MeteoInfoLab脚本示例：添加南海脚图

添加南海脚图实际上就是增加一个地图坐标系（Axesm），在用axesm函数创建地图坐标系的时候可以指定position参数（即位置和大小），第1和第2个参数是坐标系的左下角x, y位置，后两个参数是坐标系的宽度和高度，单位都是Figure的归一化坐标（0 - 1）。由于地图坐标系在Figure缩放时需要保持地图x, y坐标比例不变（地图不发生变形），坐标系具体绘制的位置会随着Figure缩放而变化，需要提前设置好Figure大小，再调整position参数。

脚本程序：

#Set data folders
basedir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'
mapdir = os.path.join(basedir, 'map')
#Read shape files
bou2_layer = shaperead(os.path.join(mapdir, 'bou2_4p.shp'))
bou1_layer = shaperead(os.path.join(mapdir, 'bou1_4l.shp'))
river_layer = shaperead(os.path.join(mapdir, 'rivers.shp'))
city_layer = shaperead(os.path.join(mapdir, 'res1_4m.shp'))
#Plot
axesm()
geoshow(bou2_layer, edgecolor='lightgray')
geoshow(bou1_layer, facecolor=(0,0,255))
ss = makesymbolspec('line', {'value':'Yangtze', 'color':(0,255,255), 'size':1}, \
    {'value':'Huang He', 'color':(0,255,255), 'size':1}, field='NAME')
geoshow(river_layer, symbolspec=ss)
geoshow(city_layer, facecolor='r', size=4, labelfield='NAME', fontname=u'楷体', fontsize=16, yoffset=15)
xlim(72, 136)
ylim(16, 55)
#Add south China Sea
sc_layer = shaperead(os.path.join(mapdir, 'bou1_4l.shp'))
axesm(position=[0.14,0.18,0.15,0.2], axison=False)
geoshow(sc_layer, facecolor=(0,0,255))
xlim(106, 123)
ylim(2, 23)

运行结果：