我正试图绘制GFS模型数据的GRIB2文件中的字段(示例文件:https://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gfs.20220202/12/atmos/gfs.t12z.pgrb2.0p25.f006(。通常情况下,我只会使用PyGRIB,昨天我就解决了这个问题,但我使用的是Windows(因为这是我雇主使用的,所以我一直在使用它,必须在Windows环境中工作(,Windows和PyGRIB玩得不好。我可以打开GRIB2文件,甚至可以使用GDAL在整个域中绘制变量。唯一的问题是,我需要一种方法来获得每个网格点的纬度和经度值的数组(类似于PyGRIB在GRIB消息上执行.llatlons(((,这样我就可以绘制域的子集。
基本上,我试图复制这个视频中正在做的事情,并且需要数据(使用dataset.GetRasterBand(269(.ReadAsArray((获得(,然后是lat/lon信息。
我也尝试过使用xarray,但Windows对xarray也不太满意。
如果您对PyGRIB感到满意,我认为解决方案是使用Conda并将其安装在Windows上。您可以使用conda-forge的miniorge来安装conda。然后,无论你得到Conda,安装pygrib与:
conda install -c conda-forge pygrib