如何使用MetPy在Python中正确打开NARR GRIB1文件

针对NARR的GRIB1文件，以下几种方法可以解决读取问题，适配248个变量的场景，且无需转存为NetCDF：

1. 优化xarray+cfgrib的读取参数 #

cfgrib默认会按GRIB的消息类型过滤变量，NARR的GRIB1文件可能存在多类消息，导致默认读取丢失变量。可以通过指定backend_kwargs调整读取规则：

import xarray as xr

# 读取所有变量，关闭默认过滤
ds = xr.open_dataset('your_narr_file.grib', engine='cfgrib',
                     backend_kwargs={'filter_by_keys': {}, 'read_all_fields': True})

# 若遇到索引缓存冲突错误，添加indexpath参数避免
ds = xr.open_dataset('your_narr_file.grib', engine='cfgrib',
                     backend_kwargs={'filter_by_keys': {}, 'read_all_fields': True, 'indexpath': ''})

# 也可按层级类型过滤（比如仅读取等压面变量）
ds = xr.open_dataset('your_narr_file.grib', engine='cfgrib',
                     backend_kwargs={'filter_by_keys': {'levelType': 'isobaricInhPa'}})

2. 使用MetPy的专用GRIB读取接口 #

MetPy对NARR数据有原生支持，能自动处理坐标系统和单位，结合GRIB读取接口可直接解析变量：

import metpy.calc as mpcalc
from metpy.io import GribFile
from metpy.units import units

# 读取GRIB文件并转为xarray数据集
with GribFile('your_narr_file.grib') as grb:
    ds = grb.to_xarray()

# 直接用MetPy计算平流、涡度等参数
u = ds['u-component_of_wind_isobaric'] * units('m/s')
v = ds['v-component_of_wind_isobaric'] * units('m/s')
vor = mpcalc.vorticity(u, v)

3. 底层用pygrib读取后构建xarray数据集 #

如果前两种方法无效，可通过pygrib直接遍历所有GRIB消息，手动构建数据集，完全控制变量读取逻辑：

import pygrib
import xarray as xr
import numpy as np

grbs = pygrib.open('your_narr_file.grib')
data_vars = {}
coords = {'time': [], 'level': []}

# 遍历每个GRIB消息，提取变量、时间和层级信息
for grb in grbs:
    var_name = grb.name.replace(' ', '_').lower()
    data = grb.values
    coords['time'].append(grb.validDate)
    coords['level'].append(grb.level)
    
    if var_name not in data_vars:
        data_vars[var_name] = []
    data_vars[var_name].append(data)

# 转换为xarray数据集
for var in data_vars:
    data_vars[var] = xr.DataArray(np.stack(data_vars[var]), dims=['time', 'level', 'y', 'x'])

ds = xr.Dataset(data_vars, coords=coords)

注意：NARR的GRIB1文件包含多类层级（地面、等压面、对流层顶等），读取时若遇到维度不一致的问题，可按层级分组处理变量。

内容的提问来源于stack exchange，提问作者user8229029