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, [email protected]
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Programmer for Rent

_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Reply via email to