计算相关系数分布

ABCD

各位高手，我用grads对一个纬度—高度剖面上的变量做相关，就是说，每个点上的值与其它所有点上的时间序列（30年）做相关，求出最小值，作为这个点上的新的变量值，但是老是提示b,c不是一个函数值或变量之类的，不知是为什么，应该做怎样的修改，谢谢大家啦！gs如下：
'reinit'
'open e:\wave\flux.ctl'
'set gxout fwrite'
'set fwrite e:\wave\corr.dat'
'set lat 0 90'
'set t 1'
iz=1
while(iz<=17)
'set z 'iz
'set t 1 30'
'define a=flux'
'set z 1 17'
'set lat 0 90'
'set t 1'
'define b=tcorr(a,flux,t=1,t=30)'
'define c=min(min(b,y=1,y=73),z=1,z=17)'
'd c'
iz=iz+1
endwhile
'disable fwrite'
;

river

什么东东啊，没看明白。不过你那个gs有够乱的，上面设置了'set lat 0 90'，循环里又设置一遍。既然你的高度是个循环的值，一次只能一层，为什么循环里还设置'set z 1 17'，这是矛盾的

mofangbao

tcorr求出来是一个数值

ABCD

river 发表于 2012-12-20 13:49

                               
登录/注册后可看大图

什么东东啊，没看明白。不过你那个gs有够乱的，上面设置了'set lat 0 90'，循环里又设置一遍。既然你的高度 ...

这个是类似求遥相关的，每个点上的flux时间序列与其它所有点上的flux时间序列（30年）做相关，求出最小值，作为这个点上的新的变量值，每个网格点都是如此，所以有了第一个循环，因为要和剖面上所有的点做相关，所以里面又设置了'set z 1 17'、'set lat 0 90'。。。

ABCD

mofangbao 发表于 2012-12-20 14:00

                               
登录/注册后可看大图

tcorr求出来是一个数值

flux是一个空间场,a是一个时间序列，用tcorr函数算出来的不应当是一个场么？然后再求出这个场中的最小值。

river

ABCD 发表于 2012-12-20 21:52

                               
登录/注册后可看大图

这个是类似求遥相关的，每个点上的flux时间序列与其它所有点上的flux时间序列（30年）做相关，求出最小值 ...

你的逻辑有问题，应该是先和所有点求相关，然后再画剖面上的相关系数，你弄反了。这个和画剖面图是一样的，你得有四维的数据才能画剖面啊。

ABCD

river 发表于 2012-12-20 22:16

                               
登录/注册后可看大图

你的逻辑有问题，应该是先和所有点求相关，然后再画剖面上的相关系数，你弄反了。这个和画剖面图是一样的 ...

是的，如果固定一个点的话，让这个点上的序列去和所有点做相关，可以得到一张剖面图，求出这个剖面图的最小值，作为这个固定点上的新的值，然后每个点都这么做，可以得到一张新的剖面图。
不用外面的循环，而是把它设置成为一个固定的点的话求出的最小值没问题，为-1，但是想对每个点都这样做，所以弄上了最外面的循环，但是这样就出错了，不知道为什么。

ABCD

各位~~求解答呀~~~

ABCD

亲们，给力呀~~我的问题问得不明白么？