如何使用 pygrib 从纬度、纬度值中检索

Question

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

为此，我从法国气象局开放数据的开放数据站点下载一个文件，方法是选择：

域："法国 - 0.01 °"，子包："SP1 - 当前表面参数"，截止时间组"0h"和运行日期"2020-02-10 00 UTC">

所以我有一个 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 周围创建一个边界框，然后取该框中值的平均值。只需调整容差即可控制框的大小。

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