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

Jillou娄

图片是陆面数据插值后的结果，但海洋信息是我们不需要的，所以想把上面的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

infocert

没看清你的要求，不好意思。看能不能找个海陆地形文件对数据做mask

longlivehj

加入
res@mpFillAreaSpecifiers=(/"water"/)

Jillou娄

longlivehj 发表于 2014-4-5 08:48
加入
res@mpFillAreaSpecifiers=(/"water"/)

请问添加这个属性的时候，还有其他附带的属性要添加吗？对于Earth..4同样适用吗？为何我添加之后，海洋上的值还没有覆盖掉

rjys1219

好东西必须顶

lbl

正好要用，谢谢分享