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

Question

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

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

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

Answer 1

在本机Python中实现这一点的一种方法是使用xarray和cfgribHere 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”是那边的标签