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

Question

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

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

我怎么能得到那个？我宁愿只使用 python 和 wgrib，但如果有任何其他简单的方法，请告诉我！我在网上找到的似乎主要涉及 wgrib2 格式。

Answer 1

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

我发布了一个示例代码，通过方法从 grib 文件的所有消息中获取 lat 和 lon 数组latlons。提取后，我将数据转换为 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)