[gradsusr] Issues displaying Himawari-8 NetCDF file in Grads 2.1
Mike Manning
michael at bsch.com.au
Sat Aug 15 00:31:57 EDT 2015
Thank you very much, MeteoInfoLab worked well!
Cheers, Mike
On 29/07/2015 5:32 PM, yaqiang.wang at gmail.com wrote:
> You may try to use MeteoInfoLab to display the Himawari-8 data. The
> example script can be found in HDF-EOS forum:
> http://hdfeos.org/forums/showthread.php?p=1588#post1588
>
>
>
> Regards
> Yaqiang
> ------------------------------------------------------------------------
> *************************************************
> Dr. Yaqiang Wang
> Chinese Academy of Meteorological Sciences (CAMS)
> 46, Zhong-Guan-Cun South Avenue
> Beijing, 100081
> China
>
> yaqiang.wang at gmail.com <mailto:yaqiang.wang at gmail.com>
>
> http://www.meteothinker.com
> **************************************************
>
> *From:* Mike Manning <mailto:michael at bsch.com.au>
> *Date:* 2015-07-26 12:18
> *To:* gradsusr <mailto:gradsusr at gradsusr.org>
> *Subject:* Re: [gradsusr] Issues displaying Himawari-8 NetCDF file
> in Grads 2.1
> Hi Jim,
>
> It's the sample Himawari-8 data on the BOM site:
> ftp://ftp.bom.gov.au/anon/sample/catalogue/Satellite/IDE00220.201507140300.nc
>
>
> Cheers, Mike
>
> On 25/07/2015 2:56 PM, James T. Potemra wrote:
>> Hi Mike,
>>
>> From where are you getting these files? The x,y coordinates are
>> in meters and not degrees. Without specifying the
>> lat/lon it looks like you are getting the full disk image. The
>> netcdf files at JMA appear to be in lat/lon
>> http://www.data.jma.go.jp/mscweb/en/himawari89/space_segment/sample_netcdf.html,
>> but as yet they only
>> seem to have the Japan-centric images.
>>
>> Jim
>>
>> On 7/23/15 8:00 PM, Mike Manning wrote:
>>>
>>> Hi all,
>>>
>>> I’m trying to display a Himawari-8 NetCDF file in Grads – here
>>> are the attributes of the nc file:
>>>
>>> dimensions:
>>>
>>> time = UNLIMITED ; // (1 currently)
>>>
>>> y = 22000 ;
>>>
>>> x = 22000 ;
>>>
>>> variables:
>>>
>>> double time(time) ;
>>>
>>> time:coordinates = "time" ;
>>>
>>> time:long_name = "time" ;
>>>
>>> time:standard_name = "time" ;
>>>
>>> time:units = "seconds since 1970-01-01 00:00:00" ;
>>>
>>> time:calendar = "standard" ;
>>>
>>> time:axis = "T" ;
>>>
>>> time:unlimited = "True" ;
>>>
>>> double y(y) ;
>>>
>>> y:coordinates = "y" ;
>>>
>>> y:long_name = "projection_y_coordinate" ;
>>>
>>> y:standard_name = "projection_y_coordinate" ;
>>>
>>> y:units = "m" ;
>>>
>>> y:axis = "Y" ;
>>>
>>> double x(x) ;
>>>
>>> x:coordinates = "x" ;
>>>
>>> x:long_name = "projection_x_coordinate" ;
>>>
>>> x:standard_name = "projection_x_coordinate" ;
>>>
>>> x:units = "m" ;
>>>
>>> x:axis = "X" ;
>>>
>>> float geostationary ;
>>>
>>> geostationary:proj4 = "+proj=geos +lon_0=140.7
>>> +h=35785863 +x_0=0 +y_0=0 +a=6378137 +b=6356752.3 +units=m
>>> +no_defs " ;
>>>
>>> geostationary:proj_name = "GEOS141" ;
>>>
>>> geostationary:grid_mapping_name = "geostationary" ;
>>>
>>> geostationary:longitude_of_projection_origin = 140.7 ;
>>>
>>> geostationary:perspective_point_height = 35785863. ;
>>>
>>> geostationary:satellite_height = 35785863. ;
>>>
>>> geostationary:semi_major_axis = 6378137. ;
>>>
>>> geostationary:semi_minor_axis = 6356752.3 ;
>>>
>>> geostationary:sweep_angle_axis = "y" ;
>>>
>>> geostationary:GeoTransform = -5500000., 500., 0., 5500000., 0.,
>>> -500. ;
>>>
>>> geostationary:spatial_ref = "PROJCS[\"unnamed\",GEOGCS[\"unnamed
>>> ellipse\",DATUM[\"unknown\",SPHEROID[\"unnamed\",6378137,298.2570248822722]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433]],PROJECTION[\"Geostationary_Satellite\"],PARAMETER[\"central_meridian\",140.7],PARAMETER[\"satellite_height\",35785863],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0]]"
>>> ;
>>>
>>> double scan_line_time(y) ;
>>>
>>> scan_line_time:units = "seconds since 2015-07-14
>>> 03:00:00" ;
>>>
>>> scan_line_time:long_name = "estimated scan line time" ;
>>>
>>> scan_line_time:calendar = "standard" ;
>>>
>>> scan_line_time:coordinates = "y" ;
>>>
>>> float channel_0003_brf(time, y, x) ;
>>>
>>> channel_0003_brf:_FillValue = 1.e+20f ;
>>>
>>> channel_0003_brf:least_significant_digit = 2s ;
>>>
>>> channel_0003_brf:algorithm_name = "channel_3_brf" ;
>>>
>>> channel_0003_brf:central_wavelength = 6.3915e-07 ;
>>>
>>> channel_0003_brf:comment = "bidirectional reflectance factor
>>> observed by the Advanced Himawari Imager at 0.64 um (Band 3)" ;
>>>
>>> channel_0003_brf:long_name = "bidirectional reflectance factor
>>> for channel 3 at 0.64 um" ;
>>>
>>> channel_0003_brf:units = "1" ;
>>>
>>> channel_0003_brf:valid_max = 1.2 ;
>>>
>>> channel_0003_brf:valid_min = 0. ;
>>>
>>> channel_0003_brf:coordinates = "time y x" ;
>>>
>>> channel_0003_brf:grid_mapping = "geostationary" ;
>>>
>>> I made a sample script as follows:
>>>
>>> 'sdfopen testfile.nc'
>>>
>>> 'd channel_0003_br'
>>>
>>> 'printim satimage.png'
>>>
>>> 'quit'
>>>
>>> Ideally I want to be able to set a lat/lon region, but when I
>>> specify any lat/lon region I just get a chart with a red
>>> output. When I don’t specify any lat/lon output at all, I get a
>>> bizzare chart (I’m guessing because it’s using x and y
>>> coordinates instead of lat/lon?) Is there any way to correctly
>>> get this to display in grads at all?
>>>
>>> Cheers, Mike
>>>
>>>
>>>
>>> _______________________________________________
>>> gradsusr mailing list
>>> gradsusr at gradsusr.org
>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150815/4e80519f/attachment-0001.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: image/png
Size: 128795 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20150815/4e80519f/attachment-0001.png
More information about the gradsusr
mailing list