从 GRIB 文件中提取特定坐标和日期/时间的天气数据

Question

我有一个带有 ECMWF 预测的 GRIB 文件，我很想根据坐标输入从中提取数据。如，提供坐标，并获取未来 5 天特定时间的预报（风速、阵风速度、风向、波高..）。

我认为 Python 可能是实现这一目标的最佳选择。有人可以指出我正确的方向吗？给我一些。

我猜二进制需要转换为 JSON（或其他可读格式），然后我可以根据提供的坐标解析并查找数据？

Answer 1

在原生 Python 中执行此操作的一种方法是使用 xarray 和 cfgrib。 Here is a tutorial。以下是教程中的关键代码：

    import xarray as xr
    ds = xr.tutorial.load_dataset('<your_grib>.grib', engine='cfgrib')

完成此操作后，grib 文件中的所有字段都将可用。一般形式是

    ds.<field_name>[<index>].values

请注意，与使用美国国家气象局提供的 GRIB 工具相比，此代码非常慢。查看degrib。大多数天气处理代码都是用 C 和 Fortran 编写的，因为它是 so much faster than Python。根据您可用的计算资源和数据大小，您可能无法在其包含的预测过期之前用 Python 处理整个 grib 文件。

最后，这个话题在GIS stack exchange 上得到了更广泛的讨论。 “grib”是那边的标签。