如何用Python移除ECMWF下载的netCDF4文件中的'expver'维度

ECMWF下载的NetCDF数据里，expver维度对应两类数据版本：1是实时分析数据，5是对缺失时段的再分析补充数据。直接删除这个维度会丢失部分有效数据，正确做法是先合并两个版本的有效数据，再移除该维度。

处理步骤 #

1. 读取原始数据

import xarray as xr
ds = xr.open_dataset("your_file.nc")

2. 合并expver维度的有效数据
   使用combine_first方法自动填充两个版本的缺失值，同时自动移除expver维度：

# 合并expver=1和expver=5的有效数据
ds_processed = ds.sel(expver=1).combine_first(ds.sel(expver=5))

3. 保存处理后的文件

ds_processed.to_netcdf("your_processed_file.nc")

验证结果 #

打印处理后的数据集，会看到expver维度已完全移除，变量结构符合预期：

print(ds_processed)

输出示例：

<xarray.Dataset>
Dimensions: (longitude: 93, latitude: 67, time: 2920)
Coordinates:

• longitude (longitude) float32 67.0 67.25 67.5 67.75 ... 89.5 89.75 90.0
• latitude (latitude) float32 23.5 23.25 23.0 22.75 ... 7.75 7.5 7.25 7.0
• time (time) datetime64[ns] 2022-01-01 ... 2022-12-31T21:00:00

Data variables:
u10 (time, latitude, longitude) float32 ...
v10 (time, latitude, longitude) float32 ...
msl (time, latitude, longitude) float32 ...
Attributes:
Conventions: CF-1.6
history: 2023-02-01 06:42:17 GMT by grib_to_netcdf-2.25.1: /opt/ecmw...

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