| 
本帖最后由 MeteoInfo 于 2020-5-30 21:59 编辑
x
登录后查看更多精彩内容~您需要 登录 才可以下载或查看,没有帐号?立即注册 
  
 Maskout通常有两种类型:Maskout图形和Maskout数据。这里是Maskout图形的示例。需要用shaperead读取地图数据形成图层作为Maskout图层(这里是中国的行政区域china.shp),然后用masklayer函数对数据图层进行屏蔽,第一个参数是用来Maskout的图层,第二次参数是一个图层列表(可以同时屏蔽多个图层)。
 
   
 脚本程序:
 
  f = addfile('D:/Temp/GrADS/model.ctl')
ps = f['PS'][0,'10:60','60:140']
axesm()
m_china = geoshow('china', visible=False)
geoshow('country', linecolor=(0,0,255))
layer = contourfm(ps, 12)
masklayer(m_china, [layer])
title('Pressure')
ylabel('Latitude')
xlabel('Longitude')
colorbar(layer)
 |