【求助】添加shp地图文件后怎么除去地图外填色部分

云雾里

求助！在用ncl制作贵州省地图时，发现自带地图有边界缺失部分，叠加贵州省地图后出现部分地区并未填色，然后想着采用shp文件直接进行处理，但是在用插值后的数据画填色图中出现超出边界部分，自己尝试了很多方法，仍没有办法解决去除边界外填色部分，还请各位大神帮帮忙抽空给指导一下，万分感谢！

画图的程序代码如下：
  minlat = 24   
  maxlat = 30  
  minlon = 103  
  maxlon = 110  

  wks  = gsn_open_wks("x11","guizhou_Map")  
  gsn_merge_colormaps(wks,"MPL_Purples","MPL_Greens")

  res                            = True            
  res@gsnMaximize               = True
  res@gsnDraw                    = False
  res@gsnFrame                   = False

  res@tmXTOn                     = False
  res@tmYROn                     = False
  res@gsnAddCyclic               = False
;------------------------------------------------------
    res@mpFillOn                      = True
    res@mpDataBaseVersion      = "LowRes"  
    res@mpOutlineOn               = False   
    res@mpLandFillColor           = "white"  

;---Zoom in on area of interest   
    res@mpLimitMode           = "LatLon"  
    res@mpMinLatF              = minlat   
    res@mpMaxLatF             = maxlat   
    res@mpMinLonF             = minlon  
    res@mpMaxLonF            = maxlon  
;-----------------------------------------------------
  res@cnInfoLabelOn          = False
  res@cnLinesOn                = False
  res@cnFillOn                   = True

  res@cnLineLabelsOn                = False
  res@cnLevelSelectionMode       = "ExplicitLevels"
  res@cnLevels                          = (/145,150,155,160,165,170,175,180,185,190,195,200,205/)

  res@lbLabelBarOn                   = True

  res@gsnSpreadColorStart              = 10
  res@gsnSpreadColorEnd               = 120
  res@lbLabelFontHeightF               = 0.01   
  res@pmLabelBarOrthogonalPosF   = 0.15   ;调节色标大小
  res@pmLabelBarHeightF               = 0.05

  map_plot = gsn_csm_contour_map(wks,gridtd,res)  ；画填色图
;****************************************************************************  
; sections for shapefiles  
;****************************************************************************  
;---加上省级区域地图   
  shpfn2="Guizhou_province.shp"  

  gres   =  True   
  gres@gsEdgesOn   =True  
  gres@gsColors       = -1      ；省内填色透明

shp_plot= gsn_add_shapefile_polygons(wks,map_plot,shpfn2,gres)  

无糖

gc_inout函数

云雾里

无糖 发表于 2019-11-29 16:21
gc_inout函数

是去掉外面的吗？试过只保存里面的，但是画出的图边界是方块的，就是一个地区一个方块一个方块的填图

我头发乱了吗

我以前画中国地图的时候也遇到过，res@cnFillDrawOrder = "PreDraw"，然后发现是画图顺序的问题，这个设置解决了，不知道你是不是一样的问题

muyuzio

你的问题解决了吗？可以参考我今天的帖子来做“使用shapefile_mask_data方法”

muyuzio

map_plot = gsn_csm_contour_map(wks,gridtd,res)  ；画填色图
这里用 gsn_csm_contour 代替，然后单独在画mapid = gsn_csm_map ,最后再画gsn_add_shapefile_polygons边界线，地图叠加等值线图

Lucia

res@mpFillDrawOrder="PostDraw"