爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 7278|回复: 6

[图形美化] 如何消除作差值后画图出现的由地形导致的颜色异常区域

[复制链接]

新浪微博达人勋

发表于 2020-5-3 13:25:36 | 显示全部楼层 |阅读模式

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

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

x

各位老师好!我在计算变量的时候设置了缺省值,可是打开两个文件做差值时右上角出现了异常的颜色块,请问如何在gs文件中消除右上角异常部分呢?请赐教,感激不尽!
部分gs文件如下:
'reinit'

'open G:\bs2\hua\hisvv_1.ctl'
'open G:\bs2\hua\vv_1.ctl'
'set parea 1 5.3 3.5 6.5'
'set gxout shaded'

'set rgb 21 0 0 133'
'set rgb 22 0 0 153'
'set rgb 23 0 0 174'
'set rgb 24 0 0 194'
'set rgb 25 0 0 214'
'set rgb 26 0 0 235'
'set rgb 27 0 0 255'
'set rgb 28 20 20 255'
'set rgb 29 41 41 255'
'set rgb 30 82 82 255'
'set rgb 31 102 102 255'
'set rgb 32 122 122 255'
'set rgb 33 163 163 255'
'set rgb 34 204 204 255'
'set rgb 35 224 224 255'
'set rgb 36 245 245 255'

'set rgb 37 255 255 255'

'set rgb 38 255 245 245'
'set rgb 39 255 224 224'
'set rgb 40 255 204 204'
'set rgb 41 255 163 163'
'set rgb 42 255 122 122'
'set rgb 43 255 102 102'
'set rgb 44 255 82 82'
'set rgb 45 255 41 41'
'set rgb 46 255 20 20'
'set rgb 47 255 0 0'
'set rgb 48 235 0 0'
'set rgb 49 214 0 0'
'set rgb 50 194 0 0'
'set rgb 51 174 0 0'
'set rgb 52 153 0 0'
'set ccols 35 31 33 35 40 41 43 45'
'set clevs -30 -20 -10 0 10 20 30'

'd vv.2-vv'
'draw title a.Win'


RCP26-HIS_MIROC-ESM-CHEM_vv_850.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2020-5-3 19:04:56 | 显示全部楼层
我认为正常情况下,GrADS会自动略过缺测值,不会参加计算的吧
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-5-3 20:46:58 | 显示全部楼层
river 发表于 2020-5-3 19:04
我认为正常情况下,GrADS会自动略过缺测值,不会参加计算的吧

可是出现了这种情况老师有没有觉得可行的解决方法呢?谢谢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-5-4 13:01:15 | 显示全部楼层
13765676798 发表于 2020-5-3 20:46
可是出现了这种情况老师有没有觉得可行的解决方法呢?谢谢

你的两个资料应该缺测值是一样的吧?两个资料是能见度资料?两个资料是时间上不同?建议把资料的描述文件也贴上
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-5-5 21:00:26 | 显示全部楼层
对的,两个资料是850hPa滤波后的北太平洋的经向风方差,缺测值应该是一样的,都是一个模式输出的历史和未来的资料
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-5-5 21:02:49 | 显示全部楼层
ctl文件是这样的:
dset G:\bs2\MRI-CGCM3-0000\rcp85\vv_850\vv_4.dat
title Daily NOAA-CIRES 20th Century Reanalysis
undef -9.96921e+36
xdef 49 linear 120.000000 2.500000
ydef 17 linear 20.000000 2.500000
tdef 1 linear 00Z01Jan2075 1dy
zdef 1 levels 850
vars 1
vv 1 99
endvars
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-5-6 08:51:32 | 显示全部楼层
有可能是你的缺测的描述不对,undef -9.96921e+36,就是这么一句,虽然你描述了缺测值,但是你图中的缺测的值不是这个。自己找出来从新定义就好了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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