tcorr函数显著性区域打点

小蜡笔没芯

刚才发了个帖子格式乱掉了= =
最近在用tcorr做相关性分析 但是显著的地方不多 就想把整个区域都shaded之后把显著相关（比如gs里的0.4）的地方打上点 但是显示的却是整个屏幕都是点点 = =然后grads显示的value result=2e+20，这个值好像是Cru_60-11-3-8这个文件里缺省值的值  不知道怎么回事TAT这是gs文件

'reinit '
'clear'
'open C:\Users\pc.pc-PC\Desktop\xg\qbjsz.ctl'
'open C:\Users\pc.pc-PC\Desktop\xg\Cru_60-11-3-8.ctl'
'set x 1'
'set y 1'
'set z 1'
'set t 1 19'
'define js=p1'
'set dfile 2'
'set lon -177.5 177.5'
'set lat -87.5 87.5'
'set z 1'
'set t 1'
'set gxout shaded'
'd tcorr(js,t3,t=1,t=19)'
'cbarn'
'cor=tcorr(js,t3,t=1,t=19)'

ilat=-87.5
while(ilat<=87.5)
'set lat 'ilat

ilon=-177.5
while(ilon<=177.5)
'set lon 'ilon

'd COR'
say result
b=subwrd(result,4)


if( b >= 0.4 | b <= -0.4 )
'q w2xy 'ilon' 'ilat''
  x=subwrd(result,3)
  y=subwrd(result,6)
'set line 1 1 4'
'draw mark 3 'x' 'y' 0.03'
endif

ilon=ilon+5
endwhile

ilat=ilat+5
endwhile

river

试试这个吧http://bbs.06climate.com/forum.php?mod=viewthread&tid=18475

小蜡笔没芯

river 发表于 2015-5-24 08:05
试试这个吧http://bbs.06climate.com/forum.php?mod=viewthread&tid=18475

谢谢 我已经弄好了 不知道为什么缺省的地方都默认相关系数成了缺省值

任我行

thank you!!!!!!!!!!!!!!!!!!!!!!!!!!

小蜡笔没芯

任我行 发表于 2015-5-24 15:57
thank you!!!!!!!!!!!!!!!!!!!!!!!!!!

linan404x

新手上路，围观学习

lj20695463

新手上路~~~~~~~~~~~~~~~~~~

梁小梁

新手上路，围观学习

Atlantic

因为我在做青藏高原冬季积雪深度与我国西南地区夏季降水两者之间的相关性，但是用grads里的tcorr函数对两者求相关时，显示的是格点坐标无效（我的方法是把高原积雪深度处理成区域平均，再和夏季6、7、8月分别作相关，格点坐标的范围是西南地区那一块），看到楼主能画出的相关场，我想问的是，楼主选取的两个变量的资料之间所在的区域是否有公共的部分。