获得给定边界的网格值的替代(更快)方法？

Question

我使用pygrib提取给定边界内特定波段的所有grib值。该函数将返回一个Json对象，该对象包含每个分析数据的所有值。

现在我已经读了大约73格里布，每次运行，反应是缓慢的，现在是大约3分钟。

结果应该是：网格结果

您知道使用GDAL或ECCODES等其他工具在bbox中提取数据的另一种方法吗？

我不能使用grib2json，因为它不认识GRIB。

提前感谢

Answer 1

你看过cfgrib了吗？它是xarray的。使用xarray，您将拥有大量的grib数据集服务和应用程序。xarray遵循数据立方体原则。

例如，如果已经将数据加载到xarray.Dataset中，则可以使用where函数将数据切片到bbox：

import xarray

dataset = xarray.open_dataset('path to grib file', engine='cfgrib')

# or for multiple files 
dataset = xarray.open_dataset(['path_to_grib1', 'path_to_grib2', ...], engine='cfgrib')

dataset.where((dataset.latitude > min_lat) & 
              (dataset.latitude < max_lat) & 
              (dataset.longitude > min_lon) &
              (dataset.longitude < max_lon), drop=True)