使用Python读取GRIB格式数据

前言

GRIB格式是一种应用于气象领域的高效存储格式，由世界气象组织进行标准化。当前有3个版本的GRIB格式，目前GRIB1和GRIB2在广泛使用。

Python读取grib文件可以选择安装xarray工具包和cfgrib。

安装

xarray github地址：https://github.com/pydata/xarray

cfgrib github地址：https://github.com/ecmwf/cfgrib

GRIB 指导文档： https://www.nco.ncep.noaa.gov/pmb/docs/grib2/

首先要安装一下ECMWF的cfgrib库，python模块依赖于ecmwf eccode 二进制库 必须安装在系统上并作为共享库访问。

在ubuntu 18.04上使用以下命令：

sudo apt-get install libeccodes0

使用pip工具安装：

pip install xarray
# 安装cfgrib需先安装二进制共享库
pip install cfgrib 

Python读取示例

# 读取grib1数据
ds = xr.open_dataset('example.grib', engine='cfgrib')
print(ds)

# 读取grib2数据
# 目前cfgrib库无法同时读取多个typeOfLevel，因此需要根据提示筛选我们需要的数据。
ds = xr.open_dataset('example.grb2', engine="cfgrib", backend_kwargs={'filter_by_keys': {'typeOfLevel': 'isobaricInhPa'}})
print(ds)

数据下载

下载地址：https://nomads.ncep.noaa.gov/