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" ; }
Attachment:
melb-small_AZE.nc
Description: Cdf file
Attachment:
melb-small_LCEA.nc
Description: Cdf file
Attachment:
melb-small_M-1SP.nc
Description: Cdf file
netcdf-java archives: