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()
原文：http://www.cnblogs.com/yaqiang/p/4585969.html