Question

我有一个netCDF文件monthly_qc_data.nc，它在0.5º的边界框中代表参数调用Lai_500m的每月值。

考虑到边界框/ netCDF文件的中心为参考点。我想计算参数Lai_500m与该参数在边界框中心处的值之差。

为此，我正在使用以下内容：

##SPATIAL VARIANCE
os.chdir(inbasedir)
data = xr.open_dataset('monthly_qc_data.nc')
ref_data = data.where((data['lat'] ==  10) & (data['lon'] == 10)) #considering the poin lat:10 and lon:10 as the center of the bounding box
dif_data = data.where((data['Lai_500m'] - ref_data))

不幸的是，这将返回以下错误：

ufunc 'bitwise_and' not supported for the input types, and the inputs could not be 
safely coerced to any supported types according to the casting rule ''safe''

我也尝试使用python netCDF4：

from netCDF4 import Dataset 
os.chdir(inbasedir)
dataset = Dataset("monthly_qc_data.nc")
dif_data = dataset.variables['Lai_500m'][:,:,:] - dataset.variables['Lai_500m'][:,10,10]

谁还返回了（显而易见的）错误：

ValueError: operands could not be broadcast together with shapes (12,120,120) (12,)

有人知道如何克服吗？

Answer 1

您应该能够将ref_data作为浮点数，然后从数据集中减去

ref_data = float(data.Lai_500m.sel(lat=10.0, lon=10.0).values)
dif_data = data.Lai_500m - ref_data

Answer 2

我知道您在寻找python的答案，但是如果万一有用，这里是您可以使用cdo从命令行执行相同功能的方法：

cdo sub in.nc -remapnn,lon=10/lat=10 in.nc diff.nc

Python-netCDF文件中参数的空间差异

2 个答案: