我正在与 Pygrib 合作,尝试使用 NBM grib 数据(如果有帮助,可在此处获得)获取特定纬度/经度坐标的表面温度。
我一直在试图获得一个索引值以用于特定纬度和经度的代表性数据。我能够得出一个索引,但问题是纬度和经度似乎各有 2 个坐标。我将以佛罗里达州迈阿密(25.7617° N,80.1918° W)为例来说明这一点。如果提供了 grib 文件,则格式化为最小可重复性。
def get_grib_data(self, gribfile, shortName):
grbs = pygrib.open(gribfile)
# Temp needs level specified
if shortName == '2t':
grib_param = grbs.select(shortName=shortName, level=2)
# Convention- use short name for less than 5 chars
# Else, use name
elif len(shortName) < 5:
grib_param = grbs.select(shortName=shortName)
else:
grib_param = grbs.select(name=shortName)
data_values = grib_param[0].values
# Need varying returns depending on parameter
grbs.close()
if shortName == '2t':
return data_values, grib_param
else:
return …Run Code Online (Sandbox Code Playgroud)