日志
画相关系数显著性区域时额外多出一条0线的解决方法
热度 3 |
RT，NCL在计算相关系数的显著性水平时有一个方便的函数rtest()，直接一条命令就可以解决。但是返回的值都是[0,1]范围内的显著性水平a，如果想区别正负相关区，就需要对结果进行重新处理下。一开始用
prob=rtest(ccr,N,0) ;ccr是相关系数
prob=where(ccr.lt.0,-prob,prob)
这样处理后画出的图多出了一条类似0线的东西，仔细看会发现在线是由正负颜色的线叠加起来的。
为了避免这种情况，可以采用下面的小方法
prob=rtest(ccr,N,0) ;ccr是相关系数
prob=1.0-prob
prob=where(ccr.lt.0,-prob,prob) ;这样的话画图的时候就可以以-0.99 , -0.95 , - 0.9 , 0.9 , 0.95 , 0.99来作为levels了，
而不是-0.1 , -0.05 , -0.01 , 0 , 0.01 , 0.05 , 0.1
图形效果如下
多出了一条线
线去掉了