如何覆盖和平均xarray打印图像？

Question

我是Python新手。我有两个netCDF文件，我使用xarray提取变量并相应地绘制它们

代码如下：

import xarray as xr
#Open dataset 1
dset = xr.open_dataset("dataset.nc")
print(dset)


#Extract tss variable from dataset 1
tss = dset['tss']
print(tss)


#Visualise 
dset['tss'].plot(cmap = 'rainbow',vmax = 20.0)

#Open dataset 2
dset2 = xr.open_dataset("datset2.nc")
print(dset2)

#Extract tss variable from dataset 2
tss_2 = dset2['tss']
print(tss_2)

#Visualise 
tss_2.plot(cmap ='rainbow', vmax = 20.0)

我的下一步是从这两个图像中生成一个平均图像。如何从这里向前迈进

Answer 1

有不同的方法可以做到这一点。例如

tss_avg = (tss*tss_2)/2

但我认为在您的情况下，有一种更聪明的方法来处理这个问题，因为xarray提供了很多有用的特性

您应该将所有数据加载到一个数据集中。（请不要在xarray documentation上显示DataArray和Dataset之间的差异）

这可以通过open_datasets实现

data = xr.open_datasets(['dataset.nc','dataset2.nc'])
 

现在可以计算变量的mean

mean_tss = data.tss.mean()
mean_tss.plot(cmap ='rainbow', vmax = 20.0)

我希望您与xarray一起工作愉快；）