<div dir="ltr">Even,<br><br>Thank you for your work and diligence. Perhaps what I meant by the warp "working" was that in 3.1.0 it is no longer throwing an error about being able to compute the bounds. <br><br>I can confirm that NWS's NBM v3.2 grib2 files are now translating to TIF successfully when I try using the latest docker images in the repo. Translating to TIF and then warping to 4326 looks correct when opening in QGIS. Thanks for the tip about band interleaving too.<br><br>Thomas</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jun 5, 2020 at 2:26 PM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><u></u>
<div style="font-family:monospace;font-size:9pt;font-weight:400;font-style:normal">
<p style="margin:0px;text-indent:0px">Thomas,</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">> A follow up with some more information. According to Zachary Flamig, there</p>
<p style="margin:0px;text-indent:0px">> are no issues in GDAL 2.2.0. I encountered these issues in 3.0.4.</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">I've just pushed a fix for that regression. was an added security check that was incorrectly triggered when reading band 23.</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">> In 3.1.0,</p>
<p style="margin:0px;text-indent:0px">> the warp works correctly, but gdal_translate still fails with an out of</p>
<p style="margin:0px;text-indent:0px">> memory issue.</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">Hum, are you really sure about that.. ? if gdal_translate cannot work, I can't see how gdalwarp could (if you use -overwrite to make sure creating a new file, and not updating something potentially unrelated)</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">I get now gdal_translate and gdalwarp to work, but due to how GRIB decompression works, for performant operation, I strongly suggest to first gdal_translate to a band interleaved GeoTIFF file, and then gdalwarp it:</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">gdal_translate blend.t16z.master.f001.co.grib2 tmp.tif \</p>
<p style="margin:0px;text-indent:0px">   -co tiled=yes -co compress=zstd -co interleave=band</p>
<p style="margin:0px;text-indent:0px">gdalwarp tmp.tif out.tif -t_srs EPSG:4326 -overwrite -co tiled=yes</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">Even</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">-- </p>
<p style="margin:0px;text-indent:0px">Spatialys - Geospatial professional services</p>
<p style="margin:0px;text-indent:0px"><a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a></p></div></blockquote></div>