覆盖netcdf中特定的单元格

Question

有没有一个netcdf的操作工具（可以是nco或者任何Python的netcdf库），可以用来覆盖netcdf文件中特定的单元格？

我想要修改一个netcdf文件中包含全球气候数据的小区域的数值。谢谢！

Answer 1

CDO可以将一个特定的经纬度“矩形”区域设置为一个值，下面是手动操作的例子：

如果你想把经度从120东到90西、纬度从20北到20南的区域内的所有值都设置为常数-999，可以使用：

cdo setclonlatbox,-999,120,-90,20,-20 infile outfile 

Answer 2

或者，NCO的ncap2可以通过一个命令来完成这个操作，具体命令是：

ncap2 -s 'var[0,0,0]=-999' in.nc out.nc

ncap2的详细介绍可以在这里找到。

Answer 3

使用 netCDF4-python 来处理这个问题非常简单。比如说，假设 nc 是你的 netCDF 文件，变量叫 var，你想要修改的单元格的索引是 (0,0,0)。那么你可以这样做：

from netCDF4 import Dataset
new_value = -999
nc.variables['var'][0,0,0] = new_value

netCDF4 使用 numpy 来表示所有的 netCDF 数组，这样就可以利用 numpy 强大的功能，比如 切片 和其他操作，来方便地处理这些数据。