GDAL中可能有bug？

Question

给定此文件：

http://mtarchive.geol.iastate.edu/2019/02/18/mrms/ncep/GaugeCorr_QPE_01H/GaugeCorr_QPE_01H_00.00_20190218-150000.grib2.gz

…我得到不同的输出取决于我使用的软件。在

wgrib2 2019021815.grib2 -csv wgrib2.csv

以及

^{pr2}$

两种输出：

-3 (undefined) :  8,869,250 records
0              : 14,848,865 records
Other values   :    781,885 records
___________________________________
Total          : 24,500,000 records

但是：

gdal_translate.exe -of xyz 2019021815.grib2 gdal.csv

输出所有24500000条记录的-3。在

如果我使用Python+GDAL解析文件，也会得到相同的结果。在

这是GDAL中的一个bug吗？我希望能够使用Python解析这些文件。在

---

我在wgrib2上使用Windows，它可以从ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib2/Windows_7/下载

我也在使用windowsforgdal，它可以从https://gdal.org/download.html下载

我正在使用Linux进行cdo，可以从 https://code.mpimet.mpg.de/projects/cdo/

Answer 1

我对你的问题的回答是改用NOAA的天气气候工具箱。在

Download Here for Windows/Mac/Linux

它可以将网格和放射状NEXRAD/MRMS数据导出到CSV。任何需要此帮助的人都可以使用示例批处理命令：

wct-export DataSetName.grib2.gz csv Output.csv wctbatchconfig.xml

如果它抛出错误，请尝试WCT的其他/早期/测试版本： https://www1.ncdc.noaa.gov/pub/data/wct/dist/并确保检查网站上的示例命令。在

Answer 2

正如你所说，这似乎是GDAL中的一个bug。我用gdal2.4.2测试了同一个命令，它可以工作。change log for 2.4.0指向GRIB驱动程序的几个修复。在

GRIB driver:

• replace DataSource? and derived classes with VSILFILE directly for > 4GB file support on Windows
• turn printf() warning as CPLDebug() messages
• read and write missing data values correctly for complex packing when original data is integer (#1063)
• g2clib: avoid int overflow. Fixes ​https://bugs.chromium.org/p/oss-fuzz/issues/detail?id=8194
• g2clib: avoid potential out of bound access (CID 1393528)
• degrib: fix floating point division by zero. Fixes ​https://bugs.chromium.org/p/oss-fuzz/issues/detail?id=9407.
• degrib: avoid potential floating point division by zero. Fixes ​https://bugs.chromium.org/p/oss-fuzz/issues/detail?id=10291
• add UNIT[] node to SRS on reading, so that is valid