MeteoInfoLab脚本示例：Maskout图形

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 = shaperead('D:/Temp/map/china.shp')
m_world = shaperead('D:/Temp/map/country1.shp')
geoshow(m_china, visible=False)
geoshow(m_world, linecolor=(0,0,255))
layer = contourfm(ps, 12)
masklayer(m_china, [layer])
title('Pressure')
ylabel('Latitude')
xlabel('Longitude')
colorbar(layer)
show()