<div dir="ltr">Simon,<div><br></div><div>I skimmed the code and there is no actual datum support in the grib reader.  It does attempt to differentiate between ellipsoids and spheres so I assume the file actually did record things as a sphere.  Generally speaking it is not clear to me that the GRIB driver can actually do better based on the information in the files but I'd be happy to be proven wrong.  In the meantime, it appears you will need to overwrite the geographic coordinate system in your ingestor. </div><div><br></div><div>Best regards,</div><div>Frank</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Mar 11, 2015 at 1:25 PM, Simon (Vsevolod) Ilyushchenko <span dir="ltr"><<a href="mailto:simonf@google.com" target="_blank">simonf@google.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>If I use GDAL to read GRIB files like these:</div><div><a href="ftp://hydro1.sci.gsfc.nasa.gov/data/s4pa/NLDAS/NLDAS_NOAH0125_H.002/2015/065/NLDAS_NOAH0125_H.A20150306.1300.002.grb" target="_blank">ftp://hydro1.sci.gsfc.nasa.gov/data/s4pa/NLDAS/NLDAS_NOAH0125_H.002/2015/065/NLDAS_NOAH0125_H.A20150306.1300.002.grb</a><br></div><div><br></div><div>the CRS I get back is:</div><div><div>GEOGCS["Coordinate System imported from GRIB file",</div><div>    DATUM["unknown",</div><div>        SPHEROID["Sphere",6371200,0]],</div><div>    PRIMEM["Greenwich",0],</div><div>    UNIT["degree",0.0174532925199433]]</div></div><div><br></div><div>Plotting these files shows that they are shifted 1.5 pixels north, but setting CRS to EPSG:4326 that uses SPHEROID["WGS 84",6378137,298.257223563 instead fixes the issue.</div><div><br></div><div>I'm trying to understand whether the source files incorrectly specifies that it's a sphere, or whether the GDAL driver is wrong, but I cannot find a GRIB reader that would show me this information about the file. What's the best way to debug this?</div><div><br></div><div>Thanks,</div><div>Simon</div></div>
<br>_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">---------------------------------------+--------------------------------------<br>I set the clouds in motion - turn up   | Frank Warmerdam, <a href="mailto:warmerdam@pobox.com" target="_blank">warmerdam@pobox.com</a><br>light and sound - activate the windows | <a href="http://pobox.com/~warmerdam" target="_blank">http://pobox.com/~warmerdam</a><br>and watch the world go round - Rush    | Geospatial Software Developer<br></div>
</div>