python - 提取 grib 文件的 Lon 和 Lat 范围

Question

我对 GRIB 格式还很陌生，但我目前需要处理它。使用 wgrib，我设法提取了我想要的字段、相应的时间等。

但是我错过了网格点的坐标，或者网格的经纬度范围。

我怎样才能得到它？我宁愿只使用 python 和 wgrib，但如果有任何其他简单的方法，请告诉我!我在网上找到的内容似乎主要与 wgrib2 格式有关。

Answer 1

如评论中所述，您必须使用模块 pygrib。

我发布了一个示例代码，通过 latlons 方法从 grib 文件的所有消息中获取经纬度数组。提取后我将数据转换为 numpy 数组:

import pygrib
import numpy as np

grib = pygrib.open('test.grib')
for g in grib:
    lt, ln = g.latlons() # lt - latitude, ln - longitude
    lt, ln = np.array(lt), np.array(ln)
    print (lt, ln)