登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 MeteoInfo 于 2020-5-30 21:59 编辑
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)
|