NOTE: The galeon
mailing list is no longer active. The list archives are made available for historical reasons.
Frank, It does not seem too hard, we have lambert_conformal_conic projection. I can get GeoX and GeoY with Coordinates Axis Types, read X and Y and find out if they link to the name found in the projection data (GeoX and GeoY) with their proper CoordinatesAxisType metadata. That should work! I'll give it a quick try this week... Regards, Denis 2005/11/7, Frank Warmerdam <warmerdam@xxxxxxxxx>:
Denis, The ncdump output looks like this: dimensions: y = 65 ; x = 93 ; variables: float Z_sfc(y, x) ; Z_sfc:long_name = "Geopotential height at surface of the earth" ; Z_sfc:standard_name = "geopotential_height" ; Z_sfc:units = "gp m" ; Z_sfc:GRIB_parameter_number = 7 ; Z_sfc:GRIB_level_flag = 1 ; Z_sfc:_FillValue = -9999.f ; Z_sfc:navigation = "nav" ; double x(x) ; x:units = "km" ; x:long_name = "synthesized X coord" ; x:_CoordinateAxisType = "GeoX" ; double y(y) ; y:units = "km" ; y:long_name = "synthesized Y coord" ; y:_CoordinateAxisType = "GeoY" ; char Projection ; Projection:Projection_Name = "lambert_conformal_conic" ; Projection:latitude_of_projection_origin = 25. ; Projection:longitude_of_central_meridian = -95. ; Projection:standard_parallel = 25. ; Projection:_CoordinateTransformType = "Projection" ; Projection:_CoordinateAxisTypes = "GeoX GeoY" ; // global attributes: :History = "GeoGrid extracted from dataset /data/mopt/jakarta-to mcat-5.0.28/content/thredds/wcs/testdata/eta.nc" ; :Convention = "_Coordinates" ; It seems to me there would be quite a bit of work making the GDAL netCDF driver handle this properly, though I am keen on it happening. :-) Best regards,
galeon
archives: