怎么给无维度信息的数组赋予经纬度信息

Divine

#==========相关性计算==============================================
r = np.zeros((preS.shape[1],preS.shape[2]))
p = np.zeros((preS.shape[1],preS.shape[2]))
for i in range(preS.shape[1]):
    for j in range(preS.shape[2]):
       if np.isnan(preS[0,i,j]) != True:
            r[i,j], p[i,j] = pearsonr(SIOI , preS[:,i,j])

在计算相关系数的时候，计算出的r和q没有继承preS里面原有维度信息，请问各位大佬有谁知道怎么将维度信息后面再加进去呢？

biubiubiu123

import xarray as xr
使用xr.DataArray()函数