如何从GNOME输出的netCDF文件中提取指定id的粒子位置数据

你之前执行命令失败的核心原因是：id是文件内的变量而非维度，你现有文件的维度只有data，所以NCO的-d参数（按维度切片）、CDO的select原生参数都不支持直接按id变量的取值筛选数据。

方案1：使用NCO实现筛选 #

直接用NCO的掩码筛选参数-m匹配id变量的取值即可，命令如下：

ncks -v longitude,latitude,id -m 'id==62001' infile.nc outputfile.nc

如果需要同时筛选多个id值，写法参考：

ncks -v longitude,latitude,id -m 'id==62001 || id==62968' infile.nc outputfile.nc

方案2：使用CDO实现筛选 #

通过expr生成条件掩码、ifthen保留匹配项的组合逻辑实现筛选，命令如下：

cdo select,name=latitude,longitude,id -ifthen -expr,'mask=(id==62968)' infile.nc infile.nc outputfile.nc

方案3：使用Python xarray库实现（适合复杂筛选场景） #

如果需要更灵活的筛选逻辑，用xarray操作更直观，代码示例：

import xarray as xr

# 读取输入文件
ds = xr.open_dataset("infile.nc")
# 筛选指定id的粒子，同时删除不匹配的维度项
filtered_ds = ds.where(ds.id == 62001, drop=True)
# 仅保留需要的经纬度、id变量
filtered_ds = filtered_ds[["longitude", "latitude", "id"]]
# 输出结果文件
filtered_ds.to_netcdf("outputfile.nc")

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