MeteoInfoLab脚本示例：站点填图

MeteoInfo

打开包含站点填图的站点数据文件（比如micaps 1）之后，用文件对象的smodeldata函数获取StationModel数据对象，然后用stationmodel函数绘制站点填图图层。

脚本程序：

1. #Set data folders
   
2. basedir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'
   
3. datadir = os.path.join(basedir, 'sample/MICAPS')
   
4. mapdir = os.path.join(basedir, 'map')
   
5. #Read shape files
   
6. lworld = shaperead(os.path.join(mapdir, 'country1.shp'))
   
7. lchina = shaperead(os.path.join(mapdir, 'bou2_4p.shp'))
   
8. #Read station data
   
9. f = addfile_micaps(os.path.join(datadir, '10101414.000'))
   
10. data = f.smodeldata()
    
11. #Plot
    
12. axesm(bgcolor=(204,255,255))
    
13. geoshow(lworld, facecolor=(255,251,195))
    
14. geoshow(lchina, edgecolor='k')
    
15. layer = stationmodel(data)
    
16. #colorbar(layer)
    
17. xlim(72, 136)
    
18. ylim(16, 55)

#Set data folders<br /> basedir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'<br /> datadir = os.path.join(basedir, 'sample/MICAPS')<br /> mapdir = os.path.join(basedir, 'map')<br /> #Read shape files<br /> lworld = shaperead(os.path.join(mapdir, 'country1.shp'))<br /> lchina = shaperead(os.path.join(mapdir, 'bou2_4p.shp'))<br /> #Read station data<br /> f = addfile_micaps(os.path.join(datadir, '10101414.000'))<br /> data = f.smodeldata()<br /> #Plot<br /> axesm(bgcolor=(204,255,255))<br /> geoshow(lworld, facecolor=(255,251,195))<br /> geoshow(lchina, edgecolor='k')<br /> layer = stationmodel(data)<br /> #colorbar(layer)<br /> xlim(72, 136)<br /> ylim(16, 55)

chongzika

前排，多谢王老师