我正在尝试从 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 消息执行 .latlons() ），以便我可以绘制域的子集。
基本上，我正在尝试复制此视频中正在执行的操作，并且需要数据（使用 dataset.GetRasterBand(269).ReadAsArray() 获取），然后是纬度/经度信息。
我也尝试过使用 xarray，但 Windows 也不能很好地使用 xarray。