替换netcdf变量并重写到新的netcdf文件

Question

有没有人可以帮我解决以下问题。我有一个具有以下维度和变量的netcdf文件:维度(大小)：y(100)，n2(2)，x(100)，变量(维度)：int64 valid_time()，int64 start_time()，float64 y(y)，float64 y_bounds(y，n2)，float64 x(x)，float64 x_bounds(x，n2)，int16沉淀(y，x)，int8 proj()

netcdf包含用于沉淀的二维数组(100,100)。现在我有了一个新的大小相同的二维数组(100,100)，它具有不同的值。我想知道如何用新的数组替换netcdf中的数组，并将其重写为新的netcdf。我尝试了以下代码，但它不能替换数组(它可以只重写file.nc并将其重命名为newfile.nc，而不替换数组)

import xarray as xr
ds=xr.open_dataset('file.nc')
precip=ds.variables['precipitation']
precip=np.array(precip)
precip=new_array
ds.to_netcdf('newfile.nc')

Answer 1

首先以2Dnumpy数组的形式提供new_array，您可以这样做：

import numpy as np
new_array = np.array(new_array)

然后你可以尝试：

import xarray as xr
ds=xr.open_dataset('file.nc')
ds.variables['precipitation'].values = new_array
ds.to_netcdf('newfile.nc)

或

import xarray as xr
ds=xr.open_dataset('file.nc')
ds['precipitation'].values = new_array
ds.to_netcdf('newfile.nc)

但要注意的是，x和y维度以及new_array映射到网格中心的方式应该与precipitation一致。