[gdal-dev] GRIB format GeoTransform error

Frank Warmerdam warmerdam at pobox.com
Fri Nov 13 09:59:16 EST 2009


Graziano Giuliani wrote:
> Hi,
> 
> I am playing in this days with the GRIB IO layer, using the API for reading GRIB
> into my application.
> 
> I have find a possible bug (I am referring to V 1.6.2, but is still there in
> SVN) on the Metadata read from file in SetGribMetaData, at lines 604-605 of the
> gribdataset.cpp source file.
> 
>            if (meta->gds.scan == GRIB2BIT_2) // Y is minY, GDAL wants maxY
>              rMaxY += (meta->gds.Ny - 1) * meta->gds.Dy; // -1 because we GDAL
> needs the coordinates of the centre of the pixel
...

Graziano,

This looks to be the same problem reported by Erick Jones within:

   http://trac.osgeo.org/gdal/ticket/2550

It would be helpful if you could add your notes there as well.  I have been
hesitant to act without knowing if I'm just breaking things for most files.

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



More information about the gdal-dev mailing list