爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 6844|回复: 8

ncl绘图

[复制链接]
发表于 2022-10-10 14:22:04 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 LLL船猫 于 2022-10-10 14:24 编辑

数据是1°的,画相关的图时左上角出现一条横线,这是数据原因吗,看起来有点怪



图片1.png

cor_nino_ww冬.pdf

2.18 MB, 下载次数: 5, 下载积分: 金钱 -5

密码修改失败请联系微信:mofangbao
发表于 2022-10-10 14:35:41 | 显示全部楼层
如果是NCL,你检查下这个区域或者附近区域是否有nan,数据的缺省值是否设定了。PS:NCL不认识Nan 画图的时候甚至把Nan画出来
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

 楼主| 发表于 2022-10-10 15:21:42 | 显示全部楼层
AuuuTi 发表于 2022-10-10 14:35
如果是NCL,你检查下这个区域或者附近区域是否有nan,数据的缺省值是否设定了。PS:NCL不认识Nan 画图的时候 ...

我刚刚仔细检查了我的数据,因为我有一个数据是事件个数,看起来是正常的,也没找到缺测的格点
(40,45,30)       6
(40,45,31)       3
(40,45,32)       3
(40,45,33)       3
(40,45,34)       3
(40,45,35)       6
(40,45,36)       9
(40,45,37)       9
(40,45,38)       9
(40,45,39)       6
(40,45,40)       6
(40,45,41)       3
(40,45,42)       3
(40,45,43)       2
(40,45,44)       0
(40,45,45)       0
(40,45,46)       1
(40,45,47)       1
(40,45,48)       0
(40,45,49)       0
(40,45,50)       1
(40,45,51)       0
(40,45,52)       0
(40,45,53)       0
(40,45,54)       0
(40,45,55)       0
(40,45,56)       0
我又输出了一下算完的相关系数值,有一些-2.147484e+09的小值,请问这个有影响吗
(35,57) 0.09292381
(35,58) -2.147484e+09
(35,59) -0.05823411
(35,60) -0.05823411
(35,61) -2.147484e+09
(35,62) -2.147484e+09
(35,63) 0.1835381
(35,64) 0.1915401
(35,65) 0.05201443
(35,66) 0.060638
(35,67) 0.1791794
(35,68) -2.147484e+09
(35,69) -2.147484e+09
(35,70) -2.147484e+09
(36,0)  -0.05129126
(36,1)  0.09311675
(36,2)  0.03411894
(36,3)  0.1547644
(36,4)  0.02968155
(36,5)  0.1351585
(36,6)  -0.1041792
(36,7)  0.1316612
(36,8)  -0.3297988
(36,9)  -2.147484e+09
(36,10) -2.147484e+09
(36,11) -2.147484e+09
(36,12) -2.147484e+09
(36,13) -2.147484e+09
(36,14) -0.05823411
(36,15) 0.06601468
(36,16) -0.07850514
(36,17) -0.04537776
(36,18) -0.08556333
(36,19) 0.0396872
(36,20) -0.08841671
(36,21) 0.1718549
(36,22) 0.1054851
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2022-10-10 15:35:28 | 显示全部楼层
LLL船猫 发表于 2022-10-10 15:21
我刚刚仔细检查了我的数据,因为我有一个数据是事件个数,看起来是正常的,也没找到缺测的格点
(40,45,3 ...

-2.147484e+09  应该就是缺省值,PrintVarSummary() 看看缺省值是不是被定义成这个数据了。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

 楼主| 发表于 2022-10-10 20:39:19 | 显示全部楼层
AuuuTi 发表于 2022-10-10 15:35
-2.147484e+09  应该就是缺省值,PrintVarSummary() 看看缺省值是不是被定义成这个数据了。

这边是输出的结果

Variable: cor
Type: float
Total Size: 15904 bytes
            3976 values
Number of Dimensions: 2
Dimensions and sizes:   [lat | 56] x [lon | 71]
Coordinates:
            lat: [   0..  55]
            lon: [  70.. 140]
Number Of Attributes: 7
  missing_value :       nan
  standard_name :       air_temperature
  long_name :   Near-Surface Air Temperature
  units :       K
  type :        real
  time :         1
  _FillValue :  nan
warning:Attempt to reference attribute (_FillValue) which is undefined
看起来缺省值是nan?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2022-10-10 21:03:33 | 显示全部楼层
LLL船猫 发表于 2022-10-10 20:39
这边是输出的结果

Variable: cor

对啊,你需要统一下缺省值 比如你那个很小的数与nan都应该被设定为非nan的缺省值,这样画图就没问题了。
http://bbs.06climate.com/home.ph ... do=blog&id=6971
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

 楼主| 发表于 2022-10-12 15:44:52 | 显示全部楼层
AuuuTi 发表于 2022-10-10 21:03
对啊,你需要统一下缺省值 比如你那个很小的数与nan都应该被设定为非nan的缺省值,这样画图就没问题了。
...

大佬我发现我的数据是完整的没有缺测的,用了rtest之后出现了缺测值nan这是什么原因啊
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2022-10-13 09:27:49 | 显示全部楼层
LLL船猫 发表于 2022-10-12 15:44
大佬我发现我的数据是完整的没有缺测的,用了rtest之后出现了缺测值nan这是什么原因啊

https://www.ncl.ucar.edu/Document/Functions/Built-in/rtest.shtml
建议看下说明
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

 楼主| 发表于 2022-10-14 16:49:59 | 显示全部楼层
AuuuTi 发表于 2022-10-13 09:27
https://www.ncl.ucar.edu/Document/Functions/Built-in/rtest.shtml
建议看下说明

谢谢你已经解决啦! 因为有些地区研究时间内没有事件发生,因此数据的方差为0,算相关系数就会出现不正常的值
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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