2

我正在尝试获取uvIndexgrib2 文件中存在的所有 lat,lng 。这是我从中获取文件的链接。问题是我无法理解文件的结构，因此我无法获取数据。我pygrib用来读取文件。

这是我尝试过的代码：

grbs = pygrib.open('uv.t12z.grbf01.grib2')
grb = grbs.select(name='UV index')[0]
print grb.data(23.5,55.5)

我想要实现的是遍历所有 lat long 并打印相应的 uvIndex 值，或者输入 lat long 并获取相应的值。阅读文档，pygrib但找不到任何适合我目的的命令。请帮忙。

4

1 回答 1

3

您必须遍历 GRIB 文件并找到所需的记录，然后获取数据，如下所示：

for g in grbs:
    print g.shortName, g.typeOfLevel, g.level # print info about all GRIB records and check names 
    if (g.shortName == shortName and g.typeOfLevel == typeOfLevel and g.level == level):
        tmp = np.array(g.values)
        # now work with tmp as numpy array

要获取 lat 和 lon 数组，请使用：lt, ln = g.latlons(), g- 的元素grbs。

在https://software.ecmwf.int/wiki/display/GRIB/GRIB+API+examples阅读 Python 部分中的示例（pygrib 使用此库来阅读 GRIB）。

从大型 GRIB 文件中获取数据的最快方法是创建索引：

# use attributes what you want to build index
indx = pygrib.index(gribfile,'typeOfLevel','level','parameterName') 

# important: msg is an array and may have more then one record
# get U wind component on 10 m above ground
msg = indx.select(level = 10, typeOfLevel = "heightAboveGround",
 parameterName = "U U-component of wind m s**-1")
u10 = np.array(msg[0].values)
# get V wind component on 10 m above ground
msg = indx.select(level = 10, typeOfLevel = "heightAboveGround",
 parameterName = "V V-component of wind m s**-1")
v10 = np.array(msg[0].values)
于 2016-06-03T10:53:01.017 回答