MeteoInfoLab脚本示例：Maskout图形

时间：2015-06-18 17:06:27 收藏：0 阅读：331

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