MeteoInfoLab脚本示例：地图投影

MeteoInfo

在用axesm函数创建地图坐标系的时候可以指定地图投影（设置projinfo参数），地图投影可以通过projinfo函数来创建，里面的参数依据proj4投影字符串，可以参考此网页：http://remotesensing.org/geotiff/proj_list。投影后的地图坐标系显示范围的设置可以用axism函数，参数是一个列表（list），列表里的4个元素分别是起始、结束经度和纬度。

脚本程序：

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. proj = projinfo(proj='lcc', lon_0=105, lat_1=25, lat_2=47)
    
11. axesm(projinfo=proj, axison=False)
    
12. geoshow(bou2_layer, edgecolor='lightgray')
    
13. geoshow(bou1_layer, facecolor=(0,0,255))
    
14. ss = makesymbolspec('line', {'value':'Yangtze', 'color':(0,255,255), 'size':1}, \
    
15.     {'value':'Huang He', 'color':(0,255,255), 'size':1}, field='NAME')
    
16. geoshow(river_layer, symbolspec=ss)
    
17. geoshow(city_layer, facecolor='r', size=4, labelfield='NAME', fontname=u'楷体', fontsize=16, yoffset=15)
    
18. axism([78, 130, 14, 53])

#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 /> proj = projinfo(proj='lcc', lon_0=105, lat_1=25, lat_2=47)<br /> axesm(projinfo=proj, axison=False)<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 /> axism([78, 130, 14, 53])

运行结果：

oubiwang

mark!!!!!!

lichenlu

老师，请问打开网格线以后，也就是gridline=True以后，如何更改网格线的风格和粗细以及颜色呢