MeteoInfo脚本教程（十九）：miapp

MeteoInfo

从MeteoInfo软件界面中打开脚本编辑器，在其中运行脚本时有一个内置对象（miapp），该对象其实就是MeteoInfo软件的主界面，通过该对象可以对MeteoInfo软件进行控制。在很多情况下，我们写脚本程序并不是为了日常定时后台运行，这时候可以将从数据中创建的图层直接加入到MeteoInfo软件界面中，从而简化脚本程序。

还是以上两讲的例子来进行示例，脚本程序中无需再创建MapLayout对象，而是通过：
    miapp.getMapDocument()    #获取MeteoInfo的图形控制对象（LayersLegend类）
    miapp.getMapDocument().getActiveMapFrame()    #获取当前活动的MapFrame对象
    miapp.getMapDocument().getActiveMapFrame().getMapView()    #获取当前活动的MapFrame对象中的MapView对象
    miapp.getMapDocument().getMapLayout()    #获取MapLayout对象

脚本程序：

1. # encoding=utf-8
   
2. #-----------------------------------------------------
   
3. # Author: Yaqiang Wang
   
4. # Date: 2014-12-6
   
5. # Purpose: Interpolate station data to grid data
   
6. # Note: Sample
   
7. #-----------------------------------------------------
   
8. #---- Import classes
   
9. print 'Import classes...'
   
10. from org.meteoinfo.data import StationData
    
11. from org.meteoinfo.data.meteodata import MeteoDataInfo, DrawMeteoData
    
12. from org.meteoinfo.legend import LegendScheme, LegendType, PolygonBreak, MapFrame
    
13. from org.meteoinfo.shape import ShapeTypes
    
14. from org.meteoinfo.layout import LayerUpdateTypes
    
15. from org.meteoinfo.geoprocess.analysis import InterpolationSetting
    
16. import os
    
17. from datetime import datetime
    
18. from datetime import timedelta
    
19. 
    
20. #---- Set directories
    
21. print 'Set directories...'
    
22. baseDir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'
    
23. dataDir = os.path.join(baseDir, 'sample/MICAPS')
    
24. mapDir = os.path.join(baseDir, 'map')
    
25. 
    
26. #---- Create MeteoDataInfo object
    
27. mdi = MeteoDataInfo()
    
28. 
    
29. #---- Set start/end time
    
30. stime = datetime(2010,10,14,14)
    
31. etime = datetime(2010,10,14,20)
    
32. 
    
33. #---- Loop
    
34. print 'Get sum station data...'
    
35. sdata = StationData()
    
36. atime = stime
    
37. i = 0
    
38. while atime < etime:
    
39.   #---- Open a MICAPS data file
    
40.   fn = os.path.join(dataDir, stime.strftime('%y%m%d%H') + '.000')
    
41.   mdi.openMICAPSData(fn)
    
42. 
    
43.   #---- Sum precipitation station data
    
44.   if i == 0:
    
45.     sdata = mdi.getStationData('Precipitation6h')
    
46.   else:
    
47.     sdata = sdata.add(mdi.getStationData('Precipitation6h'))
    
48. 
    
49.   atime = atime + timedelta(hours=6)
    
50.   
    
51. #---- Interpolate station data to grid data
    
52. print 'Interpolate station data to grid data...'
    
53. interpSet = InterpolationSetting(60,140,-20,60,160,160,"IDW_Radius",2,1)
    
54. #radList = [10.0, 8.0, 6.0, 4.0]
    
55. #interpSet = InterpolationSetting(60,140,-20,60,160,160,"Cressman",radList)
    
56. gdata = sdata.interpolateData(interpSet)
    
57. 
    
58. #---- Set legend scheme
    
59. ls = LegendScheme(ShapeTypes.Polygon)
    
60. lsfn = 'D:/Temp/rain1.lgs'
    
61. ls.importFromXMLFile(lsfn)
    
62. 
    
63. #---- Create shaded layer
    
64. print 'Create shaded layer...'  
    
65. layer = DrawMeteoData.createShadedLayer(gdata, ls, 'Rain_shaded', 'Rain', True)
    
66. layer.setMaskout(True)
    
67. 
    
68. #---- miapp object
    
69. mapFrame = miapp.getMapDocument().getActiveMapFrame()
    
70. mapLayout = miapp.getMapDocument().getMapLayout()
    
71. 
    
72. #---- Add layer
    
73. mapFrame.addLayer(layer)
    
74. mapFrame.moveLayer(layer, 0)
    
75. 
    
76. #---- Set sub title
    
77. stime = stime + timedelta(hours=-6)
    
78. subTitle = mapLayout.getTexts().get(1)
    
79. subTitle.setLabelText(u'（' + stime.strftime('%Y-%m-%d %H:00') +
    
80.   u' 至 ' + etime.strftime('%Y-%m-%d %H:00') + u'）')
    
81. 
    
82. #---- Set legend
    
83. legend = mapLayout.getLegends().get(0)
    
84. legend.setLegendLayer(layer)
    
85. legend.setTitle(u'降水量（毫米）')
    
86. legend.setLayerUpdateType(LayerUpdateTypes.NotUpdate)
    
87. 
    
88. mapLayout.paintGraphics()
    
89. 
    
90. print 'Finished!'

# encoding=utf-8<br /> #-----------------------------------------------------<br /> # Author: Yaqiang Wang<br /> # Date: 2014-12-6<br /> # Purpose: Interpolate station data to grid data<br /> # Note: Sample<br /> #-----------------------------------------------------<br /> #---- Import classes<br /> print 'Import classes...'<br /> from org.meteoinfo.data import StationData<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo, DrawMeteoData<br /> from org.meteoinfo.legend import LegendScheme, LegendType, PolygonBreak, MapFrame<br /> from org.meteoinfo.shape import ShapeTypes<br /> from org.meteoinfo.layout import LayerUpdateTypes<br /> from org.meteoinfo.geoprocess.analysis import InterpolationSetting<br /> import os<br /> from datetime import datetime<br /> from datetime import timedelta<br /> <br /> #---- Set directories<br /> print 'Set directories...'<br /> baseDir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'<br /> dataDir = os.path.join(baseDir, 'sample/MICAPS')<br /> mapDir = os.path.join(baseDir, 'map')<br /> <br /> #---- Create MeteoDataInfo object<br /> mdi = MeteoDataInfo()<br /> <br /> #---- Set start/end time<br /> stime = datetime(2010,10,14,14)<br /> etime = datetime(2010,10,14,20)<br /> <br /> #---- Loop<br /> print 'Get sum station data...'<br /> sdata = StationData()<br /> atime = stime<br /> i = 0<br /> while atime < etime:<br />   #---- Open a MICAPS data file<br />   fn = os.path.join(dataDir, stime.strftime('%y%m%d%H') + '.000')<br />   mdi.openMICAPSData(fn)<br /> <br />   #---- Sum precipitation station data<br />   if i == 0:<br />     sdata = mdi.getStationData('Precipitation6h')<br />   else:<br />     sdata = sdata.add(mdi.getStationData('Precipitation6h'))<br /> <br />   atime = atime + timedelta(hours=6)<br />   <br /> #---- Interpolate station data to grid data<br /> print 'Interpolate station data to grid data...'<br /> interpSet = InterpolationSetting(60,140,-20,60,160,160,"IDW_Radius",2,1)<br /> #radList = [10.0, 8.0, 6.0, 4.0]<br /> #interpSet = InterpolationSetting(60,140,-20,60,160,160,"Cressman",radList)<br /> gdata = sdata.interpolateData(interpSet)<br /> <br /> #---- Set legend scheme<br /> ls = LegendScheme(ShapeTypes.Polygon)<br /> lsfn = 'D:/Temp/rain1.lgs'<br /> ls.importFromXMLFile(lsfn)<br /> <br /> #---- Create shaded layer<br /> print 'Create shaded layer...'  <br /> layer = DrawMeteoData.createShadedLayer(gdata, ls, 'Rain_shaded', 'Rain', True)<br /> layer.setMaskout(True)<br /> <br /> #---- miapp object<br /> mapFrame = miapp.getMapDocument().getActiveMapFrame()<br /> mapLayout = miapp.getMapDocument().getMapLayout()<br /> <br /> #---- Add layer<br /> mapFrame.addLayer(layer)<br /> mapFrame.moveLayer(layer, 0)<br /> <br /> #---- Set sub title<br /> stime = stime + timedelta(hours=-6)<br /> subTitle = mapLayout.getTexts().get(1)<br /> subTitle.setLabelText(u'（' + stime.strftime('%Y-%m-%d %H:00') +<br />   u' 至 ' + etime.strftime('%Y-%m-%d %H:00') + u'）')<br /> <br /> #---- Set legend<br /> legend = mapLayout.getLegends().get(0)<br /> legend.setLegendLayer(layer)<br /> legend.setTitle(u'降水量（毫米）')<br /> legend.setLayerUpdateType(LayerUpdateTypes.NotUpdate)<br /> <br /> mapLayout.paintGraphics()<br /> <br /> print 'Finished!'

运行结果：

freekiller

感谢王老师无私分享大作。 你在研发，我照抄，我都跟不上你的步伐，惭愧啊。