MeteoInfoLab脚本示例：加载地图图层

应用最广泛的的地图数据应该是shape格式，网络上有很多免费下载资源。MeteoInfoLab中读取shape文件的函数是shaperead，参数即文件名，返回数据包含图形和属性信息的图层对象。矢量图层按照图元类型通常分点、线、面三种，都可以通过geoshow函数加入到地图坐标系（Axesm）中。在geoshow函数中可以指定图层的显示的Symbol（颜色等），可以设置facecolor, edgecolor, size等，这种情况下图层的所有图元都会以同一个Symbol来显示。也可以利用图层的属性信息来生成更复杂的图例（makesymbolspec函数），函数第一个参数是图元类型（point, line, polygon），还需要设置一个属性字段（field='NAME'），然后利用属性数据来设置图例。这里例子中加载了rivers.shp图层，该图层是世界主要河流，这里只显示长江、黄河，用下面的代码处理：

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函数中设置labelfield等参数来给图层添加标注（城市名）。

脚本代码如下：

#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)

运行结果：