MeteoInfo脚本教程（十八）：项目文件

MeteoInfo

在MeteoInfo中图层、地图、版面等的各种设置可以保存在一个项目文件中（.mip）。可以通过打开该项目文件还原项目的各种设置。项目文件也可以在脚本程序中使用，会很大程度降低脚本程序的复杂性（其实脚本程序中最复杂的部分主要是图形的各种设置）。比如上一讲中站点数据分析绘图的例子可以利用项目文件降低复杂性，并在某些方面提升图形质量。前面还提到了图例文件（.lgs），使用它也能是脚本中图例设置部分更加简洁。

1、首先在MeteoInfo桌面软件中根据需要编辑一个项目文件（这里是rain.mip），能更方便地控制图形中各要素，尤其是对省会城市标志位置的编辑能够提升图形质量（如果在脚本中控制每个标注的位置是非常繁琐的）。同时编辑出降水图层的图例（rain1.lgs）。


2、利用项目文件和图例文件简化上一讲的脚本程序。

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.layout import MapLayout
    
11. from org.meteoinfo.data import StationData
    
12. from org.meteoinfo.data.mapdata import MapDataManage
    
13. from org.meteoinfo.data.meteodata import MeteoDataInfo, DrawMeteoData
    
14. from org.meteoinfo.legend import LegendScheme, LegendType, PolygonBreak, MapFrame
    
15. from org.meteoinfo.shape import ShapeTypes
    
16. from org.meteoinfo.layout import LegendStyles, LayoutMap
    
17. from org.meteoinfo.global import Extent
    
18. from org.meteoinfo.projection import ProjectionInfo
    
19. from org.meteoinfo.geoprocess.analysis import InterpolationSetting
    
20. import os
    
21. from datetime import datetime
    
22. from datetime import timedelta
    
23. from java.awt import Color
    
24. from java.awt import Font
    
25. from javax.swing import JFrame
    
26. 
    
27. #---- Set directories
    
28. print 'Set directories...'
    
29. baseDir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'
    
30. dataDir = os.path.join(baseDir, 'sample/MICAPS')
    
31. mapDir = os.path.join(baseDir, 'map')
    
32. 
    
33. #---- Create MapLayout object
    
34. mapLayout = MapLayout()
    
35. 
    
36. #---- Load project file (*.mip)
    
37. pfn = 'D:/Temp/mip/rain.mip'
    
38. mapLayout.loadProjectFile(pfn)
    
39. mapFrame = mapLayout.getActiveMapFrame()
    
40. 
    
41. #---- Create MeteoDataInfo object
    
42. mdi = MeteoDataInfo()
    
43. 
    
44. #---- Set start/end time
    
45. stime = datetime(2010,10,14,14)
    
46. etime = datetime(2010,10,14,20)
    
47. 
    
48. #---- Loop
    
49. print 'Get sum station data...'
    
50. sdata = StationData()
    
51. atime = stime
    
52. i = 0
    
53. while atime < etime:
    
54.   #---- Open a MICAPS data file
    
55.   fn = os.path.join(dataDir, stime.strftime('%y%m%d%H') + '.000')
    
56.   mdi.openMICAPSData(fn)
    
57. 
    
58.   #---- Sum precipitation station data
    
59.   if i == 0:
    
60.     sdata = mdi.getStationData('Precipitation6h')
    
61.   else:
    
62.     sdata = sdata.add(mdi.getStationData('Precipitation6h'))
    
63. 
    
64.   atime = atime + timedelta(hours=6)
    
65.   
    
66. #---- Interpolate station data to grid data
    
67. print 'Interpolate station data to grid data...'
    
68. interpSet = InterpolationSetting(60,140,-20,60,160,160,"IDW_Radius",2,1)
    
69. #radList = [10.0, 8.0, 6.0, 4.0]
    
70. #interpSet = InterpolationSetting(60,140,-20,60,160,160,"Cressman",radList)
    
71. gdata = sdata.interpolateData(interpSet)
    
72. 
    
73. #---- Set legend scheme
    
74. ls = LegendScheme(ShapeTypes.Polygon)
    
75. lsfn = 'D:/Temp/rain1.lgs'
    
76. ls.importFromXMLFile(lsfn)
    
77. 
    
78. #---- Create shaded layer
    
79. print 'Create shaded layer...'  
    
80. layer = DrawMeteoData.createShadedLayer(gdata, ls, 'Rain_shaded', 'Rain', True)
    
81. layer.setMaskout(True)
    
82. 
    
83. #---- Add layer
    
84. mapFrame.addLayer(layer)
    
85. mapFrame.moveLayer(layer, 0)
    
86. 
    
87. #---- Set sub title
    
88. stime = stime + timedelta(hours=-6)
    
89. subTitle = mapLayout.getTexts().get(1)
    
90. subTitle.setLabelText(u'（' + stime.strftime('%Y-%m-%d %H:00') +
    
91.   u' 至 ' + etime.strftime('%Y-%m-%d %H:00') + u'）')
    
92. 
    
93. #---- Set legend
    
94. legend = mapLayout.getLegends().get(0)
    
95. legend.setLegendLayer(layer)
    
96. legend.setTitle(u'降水量（毫米）')
    
97. 
    
98. frame = JFrame('MeteoInfo Script Sample', size = (750, 530))
    
99. frame.add(mapLayout)
    
100. frame.visible = True
     
101. 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.layout import MapLayout<br /> from org.meteoinfo.data import StationData<br /> from org.meteoinfo.data.mapdata import MapDataManage<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 LegendStyles, LayoutMap<br /> from org.meteoinfo.global import Extent<br /> from org.meteoinfo.projection import ProjectionInfo<br /> from org.meteoinfo.geoprocess.analysis import InterpolationSetting<br /> import os<br /> from datetime import datetime<br /> from datetime import timedelta<br /> from java.awt import Color<br /> from java.awt import Font<br /> from javax.swing import JFrame<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 MapLayout object<br /> mapLayout = MapLayout()<br /> <br /> #---- Load project file (*.mip)<br /> pfn = 'D:/Temp/mip/rain.mip'<br /> mapLayout.loadProjectFile(pfn)<br /> mapFrame = mapLayout.getActiveMapFrame()<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 /> #---- 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 /> <br /> frame = JFrame('MeteoInfo Script Sample', size = (750, 530))<br /> frame.add(mapLayout)<br /> frame.visible = True<br /> print 'Finished!'

从运行结果图中可以看出省会城市标注的改进。

南宫海语

支持楼主，先收藏，日后再看

123349818

果断要学习啊。

冬蝉夏艾

二次开发的时候怎么设置项目保存、打开项目等一些菜单栏呢？