如何替换NetCDF文件中的指定变量并导出为新的NetCDF文件

你的代码失效的核心原因是仅修改了本地变量precip的指向，没有对读入的Dataset对象内部的precipitation变量做任何修改，导出时自然还是原文件的内容。

import xarray as xr
import numpy as np

# 读取原NetCDF文件
ds = xr.open_dataset('file.nc')
# 直接替换Dataset内的precipitation变量
# astype转换用于保持和原变量一致的数据类型，避免精度丢失
ds['precipitation'] = xr.DataArray(
    new_array.astype(np.int16),
    dims=ds['precipitation'].dims,
    attrs=ds['precipitation'].attrs
)
# 导出新文件
ds.to_netcdf('newfile.nc')
# 关闭数据集
ds.close()

• 替换时指定dims=ds['precipitation'].dims可以确保维度和原变量完全匹配，避免因维度顺序问题报错
• 指定attrs=ds['precipitation'].attrs可以保留原变量的所有属性（比如单位、填充值、描述信息等），如果不需要保留也可以省略该参数
• 原变量是int16类型，如果你的new_array是float类型，转换时如果存在超出int16取值范围的数值要提前做异常处理，避免数据截断
• 推荐使用上下文管理器（with语句）读取文件，无需手动调用close()，避免文件被意外占用，写法参考：

import xarray as xr
import numpy as np

with xr.open_dataset('file.nc') as ds:
    ds['precipitation'] = xr.DataArray(
        new_array.astype(np.int16),
        dims=ds['precipitation'].dims,
        attrs=ds['precipitation'].attrs
    )
    ds.to_netcdf('newfile.nc')

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