请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 21086|回复: 35

[作图] 画相关系数显著性区域时额外多出一条0线的解决方法

[复制链接]

新浪微博达人勋

发表于 2013-11-26 14:29:31 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
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
图形效果如下



                               
登录/注册后可看大图
多出了一条线


                               
登录/注册后可看大图
线去掉了


评分

参与人数 1金钱 +10 收起 理由
kongfeng0824 + 10

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-1-28 22:06:54 | 显示全部楼层
金玉姻缘 发表于 2015-1-27 15:21
楼主你好,我有模式转出来的温度数据,自己又下载了实测数据,都是二维经纬度数据,t(lat,lon),我想计算两 ...

如果只是做同期相关的话,可以用 escorc
如果要做超前滞后相关的话,可以用esccr
另外也可以自己根据相关系数公式编写一个小程序处理。
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2013-11-26 15:35:20 | 显示全部楼层
图画得不错,顶。再把地形线颜色调一下就更好了。,感觉有点乱
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-11-26 16:38:54 | 显示全部楼层
freekiller 发表于 2013-11-26 15:35
图画得不错,顶。再把地形线颜色调一下就更好了。,感觉有点乱

地形线啊。。一般我都是用的默认。。
你觉得怎么调比较好呢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-26 17:07:29 | 显示全部楼层
我喜欢风 发表于 2013-11-26 16:38
地形线啊。。一般我都是用的默认。。
你觉得怎么调比较好呢

可以用粗分辨率的地形,这样会光滑一点,另外可以用暗一点的颜色,目的就是不要让他那么明显,明显的东西应该是你关注的,比如填色部分
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-27 10:09:10 | 显示全部楼层
恩,学习了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-27 10:11:07 | 显示全部楼层
经验很好  
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-28 14:42:59 | 显示全部楼层
楼主可以把脚本贴出来吗 谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-29 13:06:16 | 显示全部楼层
多谢分享,学习一下
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-4-11 16:21:24 | 显示全部楼层
脚本呢??
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2014-4-24 08:34:55 | 显示全部楼层
正好要用,谢谢分享
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表