MeteoInfo脚本示例：台风

MeteoInfo

从台风数据文件中读取每个点的信息，创建台风路径点图层和线图层（包含了相应的属性信息，比如时间、风速等）。数据文件： 198214track.dat (2.28 KB, 下载次数: 206)

2012-9-26 12:59 上传

点击文件名下载附件





代码如下：

1. #--------------------------------------------------------        
   
2. # Author: Yaqiang Wang                                          
   
3. # Date: 2012-9-26                                                
   
4. # Purpose: Read and plot typhoon data  
   
5. # Note: Sample                                                   
   
6. #-----------------------------------------------------------  
   
7. #---- Import clr and classes
   
8. import clr
   
9. clr.AddReferenceByPartialName("System")
   
10. clr.AddReferenceByPartialName("System.Drawing")
    
11. from System import *
    
12. from System.Drawing import *
    
13. from System.Collections.Generic import *
    
14. clr.AddReference("MeteoInfoC.dll")
    
15. from MeteoInfoC import *
    
16. from MeteoInfoC.Data.MapData import *
    
17. from MeteoInfoC.Geoprocess import *
    
18. from MeteoInfoC.Layer import *
    
19. from MeteoInfoC.Shape import *
    
20. from MeteoInfoC.Legend import *
    
21. import os.path
    
22. 
    
23. #---- Create a point and line layers to plot typhoon pathway
    
24. print 'Create typhoon point layer...'
    
25. pLayer = VectorLayer(ShapeTypes.Point)
    
26. pLayer.LayerName = "Typhoon_Point"
    
27. pLayer.Visible = True
    
28. pLayer.EditAddField("Time", String)
    
29. pLayer.EditAddField("Wind", Int32)
    
30. pLayer.EditAddField("Status", String)
    
31. 
    
32. print 'Create typhoon polyline layer...'
    
33. lLayer = VectorLayer(ShapeTypes.Polyline)
    
34. lLayer.LayerName = "Typhoon_Line"
    
35. lLayer.Visible = True
    
36. lLayer.EditAddField("ID", String)
    
37. 
    
38. #---- Read typhoon data file
    
39. print 'Read typhoon data file...'
    
40. fn = "D:\\Temp\\ascii\\198214track.dat"
    
41. tf = open(fn)
    
42. tf.readline()
    
43. aline = tf.readline()
    
44. id = aline.split()[1]
    
45. tf.readline()
    
46. aPolyline = PolylineShape()
    
47. plist = List[PointD]()
    
48. for aline in tf:
    
49.         print aline
    
50.         datalist = aline.split()
    
51.         lat = float(datalist[1])
    
52.         lon = float(datalist[2])
    
53.         t = datalist[3]
    
54.         wind = int(datalist[4])
    
55.         stat = datalist[6]
    
56.         if len(datalist) == 8:
    
57.                 stat = stat + ' ' + datalist[7]
    
58.         aPS = PointShape()
    
59.         aPoint = PointD(lon, lat)
    
60.         aPS.Point = aPoint
    
61.         plist.Add(aPoint)
    
62.         shapeNum = pLayer.ShapeNum
    
63.         if pLayer.EditInsertShape(aPS, shapeNum):
    
64.                 pLayer.EditCellValue("Time", shapeNum, t)
    
65.                 pLayer.EditCellValue("Wind", shapeNum, wind)
    
66.                 pLayer.EditCellValue("Status", shapeNum, stat)
    
67. 
    
68. pLayer.UpdateLegendScheme(LegendType.GraduatedColor, "Wind")
    
69. for legend in pLayer.LegendScheme.LegendBreaks:
    
70.         legend.Size = 10
    
71.         
    
72. aPolyline.Points = plist
    
73. shapeNum = lLayer.ShapeNum
    
74. if lLayer.EditInsertShape(aPolyline, shapeNum):
    
75.         lLayer.EditCellValue("ID", shapeNum, id)
    
76. lLayer.UpdateLegendScheme(LegendType.SingleSymbol, "ID")
    
77. lLayer.LegendScheme.LegendBreaks[0].Size = 2
    
78. 
    
79. mipy.MapDocument.ActiveMapFrame.AddLayer(lLayer)
    
80. mipy.MapDocument.ActiveMapFrame.AddLayer(pLayer)   
    
81. print 'Finished!'

复制代码

表现力还可以再丰富一些，将设置点图层LegendSchem的代码修改如下：

1. for legend in pLayer.LegendScheme.LegendBreaks:
   
2.         legend.MarkerType = MarkerType.Character
   
3.         legend.FontName = "Weather"
   
4.         legend.CharIndex = 170
   
5.         legend.Size = 12

复制代码

如果想台风路径线条不同段用不同颜色表示，可以每一段生成一根线条，然后分段渲染。


代码如下：

1. #--------------------------------------------------------        
   
2. # Author: Yaqiang Wang                                          
   
3. # Date: 2012-11-18                                                
   
4. # Purpose: Read and plot typhoon data - polyline with different color
   
5. # Note: Sample                                                   
   
6. #-----------------------------------------------------------  
   
7. #---- Import clr and classes
   
8. import clr
   
9. clr.AddReferenceByPartialName("System")
   
10. clr.AddReferenceByPartialName("System.Drawing")
    
11. from System import *
    
12. from System.Drawing import *
    
13. from System.Collections.Generic import *
    
14. clr.AddReference("MeteoInfoC.dll")
    
15. from MeteoInfoC import *
    
16. from MeteoInfoC.Data.MapData import *
    
17. from MeteoInfoC.Geoprocess import *
    
18. from MeteoInfoC.Layer import *
    
19. from MeteoInfoC.Shape import *
    
20. from MeteoInfoC.Legend import *
    
21. from MeteoInfoC.Drawing import *
    
22. import os.path
    
23. 
    
24. #---- Create a point and line layers to plot typhoon pathway
    
25. print 'Create typhoon point layer...'
    
26. pLayer = VectorLayer(ShapeTypes.Point)
    
27. pLayer.LayerName = "Typhoon_Point"
    
28. pLayer.Visible = True
    
29. pLayer.EditAddField("Time", String)
    
30. pLayer.EditAddField("Wind", Int32)
    
31. pLayer.EditAddField("Status", String)
    
32. 
    
33. print 'Create typhoon polyline layer...'
    
34. lLayer = VectorLayer(ShapeTypes.Polyline)
    
35. lLayer.LayerName = "Typhoon_Line"
    
36. lLayer.Visible = True
    
37. lLayer.EditAddField("ID", String)
    
38. lLayer.EditAddField("Time", String)
    
39. lLayer.EditAddField("Wind", Int32)
    
40. lLayer.EditAddField("Status", String)
    
41. 
    
42. #---- Read typhoon data file
    
43. print 'Read typhoon data file...'
    
44. fn = "D:\\Temp\\ascii\\198214track.dat"
    
45. tf = open(fn)
    
46. tf.readline()
    
47. aline = tf.readline()
    
48. id = aline.split()[1]
    
49. tf.readline()
    
50. plist = List[PointD]()
    
51. i = 0
    
52. for aline in tf:
    
53.         print aline
    
54.         datalist = aline.split()
    
55.         lat = float(datalist[1])
    
56.         lon = float(datalist[2])
    
57.         t = datalist[3]
    
58.         wind = int(datalist[4])
    
59.         stat = datalist[6]
    
60.         if len(datalist) == 8:
    
61.                 stat = stat + ' ' + datalist[7]
    
62.         aPS = PointShape()
    
63.         aPoint = PointD(lon, lat)
    
64.         aPS.Point = aPoint
    
65.         plist.Add(aPoint)
    
66.         shapeNum = pLayer.ShapeNum
    
67.         if pLayer.EditInsertShape(aPS, shapeNum):
    
68.                 pLayer.EditCellValue("Time", shapeNum, t)
    
69.                 pLayer.EditCellValue("Wind", shapeNum, wind)
    
70.                 pLayer.EditCellValue("Status", shapeNum, stat)
    
71. 
    
72.         if plist.Count > 1:
    
73.                 aPolyline = PolylineShape()
    
74.                 aPolyline.Points = plist
    
75.                 shapeNum = lLayer.ShapeNum
    
76.                 if lLayer.EditInsertShape(aPolyline, shapeNum):
    
77.                         lLayer.EditCellValue("ID", shapeNum, id)
    
78.                         lLayer.EditCellValue("Time", shapeNum, t)
    
79.                         lLayer.EditCellValue("Wind", shapeNum, wind)
    
80.                         lLayer.EditCellValue("Status", shapeNum, stat)
    
81.                 aPoint = PointD(plist[1].X, plist[1].Y)
    
82.                 plist = List[PointD]()
    
83.                 plist.Add(aPoint)
    
84.         i += 1
    
85.        
    
86. pLayer.UpdateLegendScheme(LegendType.GraduatedColor, "Wind")
    
87. for legend in pLayer.LegendScheme.LegendBreaks:
    
88.         legend.MarkerType = MarkerType.Character
    
89.         legend.FontName = "Weather"
    
90.         legend.CharIndex = 170
    
91.         legend.Size = 12
    
92. 
    
93. lLayer.UpdateLegendScheme(LegendType.GraduatedColor, "Wind")
    
94. for legend in lLayer.LegendScheme.LegendBreaks:
    
95.         legend.Size = 2
    
96. 
    
97. mipy.MapDocument.ActiveMapFrame.AddLayer(lLayer)   
    
98. mipy.MapDocument.ActiveMapFrame.AddLayer(pLayer)
    
99. print 'Finished!'

复制代码

良辰

非常感谢王老师热情的帮助，我还是第一次用脚本绘图，挺有意思的，但是有点不懂这个语言，好像是Python，没研究过，想问一下王老师，那个设置点图层LegendSchem的代码加在什么位置，万分感激！

MeteoInfo

良辰 发表于 2012-9-26 16:28

                               
登录/注册后可看大图

非常感谢王老师热情的帮助，我还是第一次用脚本绘图，挺有意思的，但是有点不懂这个语言，好像是Python，没 ...

是Python的一个分支：IronPython，可以和.Net很好地衔接。第二段代码可以替换第一段代码的69和70行。

良辰

MeteoInfo 发表于 2012-9-26 16:31

                               
登录/注册后可看大图

是Python的一个分支：IronPython，可以和.Net很好地衔接。第二段代码可以替换第一段代码的69和70行。

谢谢王老师，我如果批处理能做吗？就是画几十条台风轨迹，fn= "F:\\taifeng\\his_typhoon_track\\194806track.dat"
tf = open(fn)，怎么修改，可以把fn定义成数组吗，进行循环打开。希望你先帮帮忙，以后要好好学习Python这个语言了，很不错，主要是metoinfo的脚本处理！

良辰

F:\未命名.jpg出错了

良辰

MeteoInfo

良辰 发表于 2012-9-26 16:38

                               
登录/注册后可看大图

谢谢王老师，我如果批处理能做吗？就是画几十条台风轨迹，fn= "F:\\taifeng\\his_typhoon_track\\194806t ...

脚本程序很适合来做批处理。你可以把所有的台风文件放在一个目录中，在脚本中可以列出目录中所有的文件，获得一个文件列表，通过循环列表中每个文件就可以了。

dir = 'F:\\taifeng\\his_typhoon_track'
list = os.listdir(dir)  #列出目录下的所有文件和目录
for line in list:
    fn = os.path.join(dir,line)
    ....    #例子里的代码， 注意代码缩进。

MeteoInfo

良辰 发表于 2012-9-26 16:42

                               
登录/注册后可看大图

是代码缩进的问题，你具体在网上查查Python语言代码缩进。

良辰

MeteoInfo 发表于 2012-9-26 16:53

                               
登录/注册后可看大图

是代码缩进的问题，你具体在网上查查Python语言代码缩进。

王老师，万分感谢，谢谢谢谢

良辰

良辰 发表于 2012-9-26 16:56

                               
登录/注册后可看大图

王老师，万分感谢，谢谢谢谢

它说legend.MarkerType = MarkerType.Character，没有定义，这又是为嘛？