python 如何读取GFS数据（.grib文件）

python 可以通过xarray 包对.grib文件进行读取,以及数据处理。

环境准备：

• 1、安装anaconda
• 2、安装xarray包：conda install -c conda-forge xarray
• 3、安装解码库"eccodes"：conda install -c conda-forge eccodes
• 4、安装cfgrib 库：conda install -c conda-forge cfgrib

上述的包都可以在windows环境下安装

安装完成后直接，进行读取：

import xarray as xr
path=r"D://gfs_data//225.f000"
wind = xr.open_dataset(path, engine='cfgrib')

注意的是，这里提示报错了：

这是因为，你的.grib文件里包含了太多了关键字，也就是typeOfLevel，因此需要根据提示筛选我们需要的数据。增加一行命令，读取你需要的数据即可。

wind = xr.open_dataset(path, engine='cfgrib',\
backend_kwargs={'filter_by_keys': {'typeOfLevel': 'isobaricInhPa'}})

这样就显示读取成功了：

下面的一些操作就和nc文件一样啦，这里不做赘述。

一些测试数据可以在下面的网址进行下载：

测试数据