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

MeteoInfo

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

脚本程序：

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

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

运行结果：

karenlk

这怎么没人抢沙发，好棒！

aas314

好好好  学习 感谢

兔子屎

实用。。。。。。。。。。

chenfei0430

学习了，好东西要分享

lorraine

好东西，学习了。

慧仔333

非常有用，谢谢！

黎子

老师您好，请问要怎么添加横纵坐标轴标题~~

MeteoInfo

黎子 发表于 2017-9-6 10:41
老师您好，请问要怎么添加横纵坐标轴标题~~

参考这里：http://www.meteothinker.com/docs ... nctions/xlabel.html

黎子

MeteoInfo 发表于 2017-9-6 10:57
参考这里：http://www.meteothinker.com/docs/meteoinfolab/plot/functions/xlabel.html

这个明白，我的意思是直接导数据画站点分布图时，可以添加坐标轴标题吗？