NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
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/ >
netcdf-java archives: