Python：两个二维数组之间的相关系数

Question

假设我们有两个不同的 2D np.arrays

a = (np.random.randint(25,30,1512000)).reshape(1050,1440)

b = (np.random .randint(20,25,1512000)).reshape(1050,1440)

所以，维度a 和 b 是 (1050,1440)

我想计算每个网格点 a 和 b 之间的互相关系数（即创建一个新的二维数组，其中包含 a 和 b 之间的相关系数值，维度为 (1050,1440)。

Answer 1

您可以使用 nctoolkit 来完成此操作。假设文件具有相同的网格，则以下内容将起作用：

import nctoolkit as nc
ds1 = nc.open_data("file1.nc")
ds2 = nc.open_data("file2.nc")
ds_cor = nc.cor_space(ds1, ds2)
ds_cor.plot()

可能网格不匹配，因此您可能需要重新网格：

import nctoolkit as nc
ds1 = nc.open_data("file1.nc")
ds2 = nc.open_data("file2.nc")
# regrid ds1 to ds2's grid
ds1.regrid(ds2)
ds_cor = nc.cor_space(ds1, ds2)
ds_cor.plot()