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

Question

Let's assume that we have two different 2D np.arrays

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

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

So, the dimension of a and b is (1050,1440)

I want to calculate the cross correlation coefficient between a and b at each grid point (ie to create a new 2D array containing correlation coefficient values between a and b, with a dimension of (1050,1440).

Answer 1

You can do this with nctoolkit. Assuming the files have the same grid, the following will work:

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()

Likely the grids do not match, so you might want to regrid:

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()