求问空间自相关系数

979162804

oyy. 发表于 2020-3-10 10:24
是对的

大佬程序在下面，我时间相关系数不知道为啥都是-9999缺省
begin





File1 = systemfunc("ls yb*.nc")
;print(dimsizes(File1))
files1=addfiles(File1,"r")
ListSetType(files1,"join")

yubao=files1[:]->P

ave_yubao = dim_avg_n_Wrap(yubao,0)
;print(ave_yubao)
;printVarSummary(yubao)


File2 = systemfunc("ls gc*.nc")
;print(dimsizes(File2))
files2=addfiles(File2,"r")
ListSetType(files2,"join")

guance=files2[:]->P


ave_guance = dim_avg_n_Wrap(guance,0)
;print(ave_guance)

yubao_1=yubao(lat |:,lon |:,ncl_join |:)
guance_1=guance(lat |:,lon |:,ncl_join |:)

ra = pattern_cor(yubao(0,:,:),guance(0,:,:),1.0,0)!求空间相关系数
print(ra)

es=escorc_n(yubao,guance,0,0)！求时间相关系数
print(es)


end

979162804

oyy. 发表于 2020-3-10 10:24
是对的

奇怪的是，我把join改为cat方式后，都可以算出来值，但是时间相关系数算出来是2个，这是正确的吗（我只选了2天的数据），

oyy.

http://www.ncl.ucar.edu/Document ... t-in/escorc_n.shtml
你可以看一下官网escorc_n的介绍。
The following will calculate the cross-correlations for x3(time,lat,lon) and y3(time,lat,lon)

     ccr3 = escorc_n(x3,y3,0,0)      ; ccr3(nlat,mlon)

天海人

oyy. 发表于 2020-3-10 10:24
是对的

我想问问我是对模式资料和再分析资料分别做了mveof分解后对二者得到的eof_ts时间系数做ACC分析可否也用这个函数

oyy.

天海人 发表于 2020-4-7 20:42
我想问问我是对模式资料和再分析资料分别做了mveof分解后对二者得到的eof_ts时间系数做ACC分析可否也用这 ...

一维数组直接用escorc就可以啦

Xiaoooooooo-

oyy. 发表于 2020-3-9 12:25
pattern_cor 或者把数据变成1维

您好！想请问您：我想算模式资料与观测资料之间的距平相关系数（Anomaly Correlation Coefficient），已经利用双线性插值方法，将格点资料插值到了站点上，但是根据公式进行计算时，不知该如何对站点数据进行区域平均，请问这个情况是要一个站一个站的算吗？
谢谢您！