[gdal-dev] Frank Eats Crow - or Fixing PixelIsPoint Interpretation

[Frank Warmerdam]

Joaquim Luis wrote:
> Frank,
> 
> Wouldn't this apply to the netCDF grids as well? GMT has always made the 
> distinction on the area vs point representation. The grids created by 
> GMT hold that information on the global attribute "node_offset". As for 
> other netCDF CF grids created elsewhere it is also possible to find 
> which model is used by checking the grid limits vs array size because 
> for the same limits the PixelIsArea model has one less row and column 
> than the PixelIsPoint.

Joaquim,

The gdal/frmts/netcdf/gmtdataset.cpp code already examines the node_offset
and adjusts the geotransform accordingly.  It does not currently attempt
to expose the node_offset via the AREA_OR_POINT metadata item though it
could.

I do not see the GMT format driver being pertinent to the RFC.

Best regards,
-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Programmer for Rent