On Fri, Oct 14, 2011 at 5:00 AM, <netcdf-java-request@xxxxxxxxxxxxxxxx> wrote:
> Message: 1
> Date: Thu, 13 Oct 2011 11:09:20 -0600
> From: John Caron <caron@xxxxxxxxxxxxxxxx>
> To: netcdf-java@xxxxxxxxxxxxxxxx
> Subject: Re: [netcdf-java] Can't open Gridded data in NetCDF Java
> based tools (ToolsUI, IDV) for AZE, LCEA, Mercator projections?
> Message-ID: <4E971B40.4030004@xxxxxxxxxxxxxxxx>
> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed"
>
> Hi Patrick:
>
> Are you using GDAL for this? if so, best to ask the question on their list.
>
> The CF compliance checker does not check that the file contains gridded
> data, so it is not definitive.
Hi John,
yes the files are being created by GDAL originally, but I have checked
the attributes
and they do seem to be appropriate for those projections in the CF-1.5
conventions doc.
The attributes of the required names are there, so unless some of the
values don't make sense in CF-1.5?
If there are any test NetCDF datasets available in these projections
that NetCDF-Java should be able to open
that would help as well.
thanks, Patrick.
>
> John
>
> On 9/29/2011 1:17 AM, Patrick Sunter wrote:
>> Hello all,
>>
>> I've been doing some testing with exporting to NetCDF CF-1 files from
>> GeoTiff's using GDAL ... and am having problems opening 3 files that
>> seem to be CF-compliant (IE passed the online compliance checker, and
>> their projections have the right attributes according to
>> http://cf-pcmdi.llnl.gov/documents/cf-conventions/1.5/apf.html).
>>
>> Can anyone see an issue here? The Azimuthal Equidistant and Lambert
>> Cylindrical Equal Area files refuse to be recognised as gridded data,
>> whereas the Mercator file loads but IDV/ToolsUI seem unable to
>> recognise the coordinates properly and display in the right spot.
>>
>> Files attached and ncdumps listed below.
>>
>> thanks, Patrick.
>>
>> --
>>
>> netcdf melb-small_AZE {
>> dimensions:
>> x = 193 ;
>> y = 83 ;
>> variables:
>> char azimuthal_equidistant ;
>> azimuthal_equidistant:grid_mapping_name =
>> "azimuthal_equidistant" ;
>> azimuthal_equidistant:semi_major_axis = 6378137. ;
>> azimuthal_equidistant:inverse_flattening = 298.257223563 ;
>> azimuthal_equidistant:latitude_of_projection_origin = -37. ;
>> azimuthal_equidistant:longitude_of_projection_origin = 145. ;
>> azimuthal_equidistant:false_easting = 0. ;
>> azimuthal_equidistant:false_northing = 0. ;
>> azimuthal_equidistant:proj4 = "+proj=aeqd +lat_0=-37 +lon_0=145
>> +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs " ;
>> azimuthal_equidistant:spatial_ref =
>> "PROJCS[\"unnamed\",GEOGCS[\"WGS
>> 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS
>> 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433],AUTHORITY[\"EPSG\",\"4326\"]],PROJECTION[\"Azimuthal_Equidistant\"],PARAMETER[\"latitude_of_center\",-37],PARAMETER[\"longitude_of_center\",145],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]]]"
>> ;
>> azimuthal_equidistant:GeoTransform = "75048.08411999565
>> 30.0083799812434 0 -7954.467523806079 0 -30.0083799812434 " ;
>> double x(x) ;
>> x:standard_name = "projection_x_coordinate" ;
>> x:long_name = "x coordinate of projection" ;
>> x:units = "m" ;
>> double y(y) ;
>> y:standard_name = "projection_y_coordinate" ;
>> y:long_name = "y coordinate of projection" ;
>> y:units = "m" ;
>> byte Band1(y, x) ;
>> Band1:long_name = "GDAL Band Number 1" ;
>> Band1:grid_mapping = "azimuthal_equidistant" ;
>>
>> // global attributes:
>> :GDAL_AREA_OR_POINT = "Area" ;
>> :Conventions = "CF-1.5" ;
>> :GDAL = "GDAL 1.9dev, released 2011/01/18" ;
>> }
>>
>> ===
>>
>> netcdf melb-small_LCEA {
>> dimensions:
>> x = 192 ;
>> y = 81 ;
>> variables:
>> char lambert_cylindrical_equal_area ;
>> lambert_cylindrical_equal_area:grid_mapping_name =
>> "lambert_cylindrical_equal_area" ;
>> lambert_cylindrical_equal_area:semi_major_axis = 6378137. ;
>> lambert_cylindrical_equal_area:inverse_flattening =
>> 298.257223563 ;
>> lambert_cylindrical_equal_area:longitude_of_central_meridian =
>> 145.f ;
>> lambert_cylindrical_equal_area:false_easting = 0.f ;
>> lambert_cylindrical_equal_area:false_northing = 0.f ;
>> lambert_cylindrical_equal_area:standard_parallel = -37. ;
>> lambert_cylindrical_equal_area:proj4 = "+proj=cea +lon_0=145
>> +lat_ts=-37 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs " ;
>> lambert_cylindrical_equal_area:spatial_ref =
>> "PROJCS[\"unnamed\",GEOGCS[\"WGS
>> 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS
>> 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433],AUTHORITY[\"EPSG\",\"4326\"]],PROJECTION[\"Cylindrical_Equal_Area\"],PARAMETER[\"standard_parallel_1\",-37],PARAMETER[\"central_meridian\",145],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]]]"
>> ;
>> lambert_cylindrical_equal_area:GeoTransform =
>> "75046.57819714122
>> 30.00481221042836 0 -4789456.678646537 0 -30.00481221042836 " ;
>> double x(x) ;
>> x:standard_name = "projection_x_coordinate" ;
>> x:long_name = "x coordinate of projection" ;
>> x:units = "m" ;
>> double y(y) ;
>> y:standard_name = "projection_y_coordinate" ;
>> y:long_name = "y coordinate of projection" ;
>> y:units = "m" ;
>> byte Band1(y, x) ;
>> Band1:long_name = "GDAL Band Number 1" ;
>> Band1:grid_mapping = "lambert_cylindrical_equal_area" ;
>>
>> // global attributes:
>> :GDAL_AREA_OR_POINT = "Area" ;
>> :Conventions = "CF-1.5" ;
>> :GDAL = "GDAL 1.9dev, released 2011/01/18" ;
>> }
>>
>> =====
>>
>> (This Mercator does actually load, but then doesn't project properly)
>>
>> netcdf melb-small_M-1SP {
>> dimensions:
>> x = 192 ;
>> y = 81 ;
>> variables:
>> char mercator ;
>> mercator:grid_mapping_name = "mercator" ;
>> mercator:semi_major_axis = 6378137. ;
>> mercator:inverse_flattening = 298.257223563 ;
>> mercator:longitude_of_projection_origin = 145.f ;
>> mercator:scale_factor_at_projection_origin = 1.f ;
>> mercator:false_easting = 0.f ;
>> mercator:false_northing = 0.f ;
>> mercator:proj4 = "+proj=merc +lon_0=145 +k=1 +x_0=0 +y_0=0
>> +datum=WGS84 +units=m +no_defs " ;
>> mercator:spatial_ref = "PROJCS[\"unnamed\",GEOGCS[\"WGS
>> 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS
>> 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433],AUTHORITY[\"EPSG\",\"4326\"]],PROJECTION[\"Mercator_1SP\"],PARAMETER[\"central_meridian\",145],PARAMETER[\"scale_factor\",1],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]]]"
>> ;
>> mercator:GeoTransform = "93968.49658325009 37.56813175897275 0
>> -4422925.417953242 0 -37.56813175897275 " ;
>> double x(x) ;
>> x:standard_name = "projection_x_coordinate" ;
>> x:long_name = "x coordinate of projection" ;
>> x:units = "m" ;
>> double y(y) ;
>> y:standard_name = "projection_y_coordinate" ;
>> y:long_name = "y coordinate of projection" ;
>> y:units = "m" ;
>> byte Band1(y, x) ;
>> Band1:long_name = "GDAL Band Number 1" ;
>> Band1:grid_mapping = "mercator" ;
>>
>> // global attributes:
>> :GDAL_AREA_OR_POINT = "Area" ;
>> :Conventions = "CF-1.5" ;
>> :GDAL = "GDAL 1.9dev, released 2011/01/18" ;
>> }
>>
>>
>> _______________________________________________
>> netcdf-java mailing list
>> netcdf-java@xxxxxxxxxxxxxxxx
>> For list information or to unsubscribe, visit:
>> http://www.unidata.ucar.edu/mailing_lists/
>
