[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


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


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