將低於閾值的值設置為 .netcdf 文件中的閾值

Question

我想在 a.netcdf 文件中將所有低於常量c的值設置為c本身：file.nc

使用氣候數據運營商 (CDO) 的解決方案是

cdo mul -gec,$c file.nc file.nc t1.nc
cdo add -mulc,$c -ltc,$c file.nc t1.nc output.nc
rm -f t1.nc

但是有沒有更簡潔/更短的方法來做到這一點？

Answer 1

您可以使用NCO 的 ncap2輕松完成此操作。

例如，在file.nc x低於 100 的所有值設置為 100，並在file2.nc輸出：

>>> ncap2 -s 'where(x<100.) x=100;' file.nc -O file2.nc 

Answer 2

ncap2 的裁剪操作符最為簡潔：

ncap2 -s 'x=x>>100' in.nc out.nc

Answer 3

將 Python 與 NumPy 和 netCDF4 結合使用，您可以執行以下操作：

import numpy as np
from netCDF4 import Dataset

dataset = Dataset('/path/to/dataset','r+')

data = dataset.variables['data_variable_name'][:]

threshold = 100 # or whatever your constant c is

# np.where(condition, value if true, value if false)
new_data = np.where(data < threshold, threshold, data)

# Write your new data back to the NetCDF file
dataset.variables['data_variable_name'][:] = new_data[:]

dataset.close()

祝你好運！

Answer 4

還有一種更有效的方法可以使用expr與氣候數據運算符 (CDO) 一起執行此操作。 將這篇文章中的示例改編為這個問題中的具體示例，假設變量的名稱是x並且常量是c ，那么命令將是這樣的：

cdo -expr,'x = ((x > c)) ? x : c' infile.nc outfile.nc

在哪里? and :是三元條件運算符，而x? y: z x? y: z ， y if x not equal 0, else z 。 有關運算符和表達式的其他信息，請參見CDO 手冊的第 2.7.1 節。