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

问题描述

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

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

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

解决方法

在原生 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”是那边的标签。