<div dir="ltr"><div style="font-size:12.8000001907349px">Turns out the GDAL driver is doing the right thing. <span style="font-size:12.8000001907349px">However, many other tools that read GRIB files default to using the spheroid regardless of what the files specify. Consequently, many data providers don't bother switching this setting from the default sphere to the spheroid. We tried to use the sphere instead of the intended spheroid, which was the wrong thing to do for the dataset in question.</span></div><div class="" style="font-size:12.8000001907349px"></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Mar 11, 2015 at 5:17 PM, Frank Warmerdam <span dir="ltr"><<a href="mailto:warmerdam@pobox.com" target="_blank">warmerdam@pobox.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">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"><div><div class="h5">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></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><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></div></div>_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">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><span class="HOEnZb"><font color="#888888"><br></font></span></blockquote></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div>---------------------------------------+--------------------------------------<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>
</font></span></div>
</blockquote></div><br></div>