我对GRIB格式还很陌生,但是我现在需要处理它。使用wgrib,我设法提取了想要的字段,对应的时间等。
但是我想念网格点的坐标,或者网格的纬度和经度范围。
我怎么能得到那个? 我更喜欢只使用python和wgrib,但是如果还有其他简单的方法,请告诉我! 我在网上发现的内容似乎主要与wgrib2格式有关。
答案 0 :(得分:2)
如注释中所述,您必须使用模块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)