爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 6824|回复: 5

[y已解决]如何快速的将海洋上的contour去掉

[复制链接]
发表于 2014-4-4 22:42:21 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 Jillou娄 于 2014-4-5 16:41 编辑

111.png
图片是陆面数据插值后的结果,但海洋信息是我们不需要的,所以想把上面的contour覆盖掉。尝试出图设置res,但好像没通,请教各位。
出图部分的代码就是这样了,可是眼拙,检查了半天没查来原因。
因为数据是站点信息,所以用mask不太好整。主要思路参考的http://bbs.06climate.com/forum.php?mod=viewthread&tid=11417

解决思路:先画contour,然后把地图轮廓填上去,再把海洋区域设成transparent~~~~

127 plot = new(6,graphic)
128 res=True
129 res@cnInfoLabelOn           = False
130 res@gsnDraw              = False        ; do not draw picture
131 res@gsnFrame             = False       ; do not advance frame
132
133 wks=gsn_open_wks("pdf","interpolation")
134 gsn_define_colormap(wks,"WhiteBlueGreenYellowRed")
135 res@mpMinLatF            = 20.5        ; zoom in on map
136 res@mpMaxLatF            = 40
137 res@mpMinLonF            = 104
138 res@mpMaxLonF            = 124.75
139 res@cnFillOn             = True     ; turn on color fill
140 res@cnLinesOn            = False    ; turn of contour lines
res@cnFillDrawOrder      = "Predraw"  ; draw contours firstres@mpLandFillColor       = "transparent"
res@mpOceanFillColor = 0   
142 res@cnLevelSpacingF      = 0.03        ; interval
143 res@gsnSpreadColors     = True         ; Use full range of colormap.
144 res@gsnSpreadColorStart =97      ;color to start from the last color in the color table      
145 res@gsnSpreadColorEnd   =152     ; color to end
147 res@mpOutlineOn          = True
148 res@mpLeftCornerLatF    = 20.5
149 res@mpLeftCornerLonF    = 104
150 res@mpRightCornerLatF    = 40

151 res@mpRightCornerLonF   = 124.75
152 res@mpDataBaseVersion     = "Ncarg4_1"          ; DEFAULT LOW RESOLUTION OR "HighRes"
153 res@mpLandFillColor       = 177
154 res@mpOutlineBoundarySets = "Geophysical"
155 res@gsnAddCyclic          = False               ; Not global,set False

157 res@mpDataSetName = "Earth..4"
158 res@mpDataBaseVersion = "MediumRes"
159 res@mpOutlineSpecifiers         = (/"Taiwan"/)
160 res@mpGeophysicalLineThicknessF = 0.5
161 ;;;Turn on Mask;;;;;
162 res@mpMonoFillColor = False
163 res@mpFillColor = True
164 res@mpAreaMaskingOn = True
165 res@mpMaskAreaSpecifiers = (/"China:states","Taiwan"/)   ;China:states
166 res@mpOceanFillColor = 0
QQ截图20140405164047.png





密码修改失败请联系微信:mofangbao
发表于 2014-4-4 23:22:20 | 显示全部楼层
本帖最后由 infocert 于 2014-4-4 23:26 编辑

没看清你的要求,不好意思。看能不能找个海陆地形文件对数据做mask
密码修改失败请联系微信:mofangbao
发表于 2014-4-5 08:48:01 | 显示全部楼层
加入
res@mpFillAreaSpecifiers=(/"water"/)
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-4-5 09:41:52 | 显示全部楼层
longlivehj 发表于 2014-4-5 08:48
加入
res@mpFillAreaSpecifiers=(/"water"/)

请问添加这个属性的时候,还有其他附带的属性要添加吗?对于Earth..4同样适用吗?为何我添加之后,海洋上的值还没有覆盖掉
密码修改失败请联系微信:mofangbao
发表于 2014-4-16 10:25:02 | 显示全部楼层
好东西必须顶
密码修改失败请联系微信:mofangbao
发表于 2014-4-24 08:36:09 | 显示全部楼层
正好要用,谢谢分享
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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