python - GDAL 中可能的错误？

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

和

cdo outputtab,date,time,lat,lon,value 2019021815.grib2 > cdo.txt

两个输出:

-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

输出 -3所有 24,500,000 条记录。

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

这是GDAL中的错误吗？我希望能够使用 Python 解析这些文件。

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

我也在使用 Windows 进行 GDAL，可以从 https://gdal.org/download.html 下载

我使用 Linux 作为 cdo，可以从
https://code.mpimet.mpg.de/projects/cdo/

Answer 1

我对您的问题的回答是改用 NOAA 的 Weather-Climate Toolkit。

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/并确保检查网站上的示例命令。