如何使用pygrib从lat，lon值中检索

Question

我从python编程开始，我想做一个小脚本，它显示“在浮雕上方2米高度诊断的局部温度”，“在10米高度诊断的东西向水平风的纬向分量”和“在10米高度诊断的水平风的子午线分量”作为经度和纬度的函数。

为此，我从Meteofrance OPEN DATA的开放数据站点下载一个文件，方法是选择：

域名："France 0.01°"，子包："SP1 -当前表面参数“，截止日期组"0h”，运行日期"2020-02-10 00UTC“

所以我有一个grib2格式的文件，我正在尝试用pygrib库来处理它

为了简单起见，我只是尝试获取给定点的温度(经度= 0.25，纬度= 49.21)

我创建了一个索引来读取文件(它似乎是最快的)

indx = pygrib.index('./AROME_0.01_SP1_00H_2020021000.grib2', 'typeOfLevel', 'level', 'name')

我选择与温度相对应的记录，然后恢复这些值：

msg = indx.select(level=2, typeOfLevel="heightAboveGround", name="2 metre temperature")
temp2m = msg[0].values

问题是，从那里开始，我摸索着找不到如何从这个变量(numpy.ma.core.MaskedArray)中检索对应于我的经度= 0.25和纬度= 49.21的值

如果有人有主意，我很感兴趣

import pygrib
indx = pygrib.index('./AROME_0.01_SP1_00H_2020021000.grib2', 'typeOfLevel', 'level', 'name') 
msg = indx.select(level=2, typeOfLevel="heightAboveGround", name="2 metre temperature")
temp2m = msg[0].values

Answer 1

您可以使用pygrib.data()方法在经纬度周围创建一个边界框，然后取该框中的值的平均值。只需调整公差来控制长方体的大小。

lat, lon, tolerence = 49.21, 0.25, 0.1
data, lats, lons = grb.data(lat1=lat-tolerence,lat2=lat+tolerence,
                            lon1=lon-tolerence,lon2=lon+tolerence)
data.mean()

Answer 2

您应该能够通过从您选择变量(而不是值)中使用.latlons()来获得您的经度和经度。从那里你可以得到与你的观点相对应的索引。

import pygrib
indx = pygrib.index('./AROME_0.01_SP1_00H_2020021000.grib2', 'typeOfLevel', 'level', 'name') 
msg = indx.select(level=2, typeOfLevel="heightAboveGround", name="2 metre temperature")

lat,lon = msg.latlons()

temp2m = msg[0].values