MeteoInfo脚本示例：生成可见光云图

MeteoInfo

读取MICAPS格式云图数据并生成图片。

脚本代码：

1. # This Python file uses the following encoding: GB2312
   
2. #--------------------------------------------------------        
   
3. # Author: Yaqiang Wang                                          
   
4. # Date: 2013-8-27                                                
   
5. # Purpose: Read and plot MICAPS cloud data  
   
6. # Note: Sample                                                   
   
7. #-----------------------------------------------------------  
   
8. #---- 引入类库
   
9. import clr
   
10. clr.AddReferenceByPartialName("System")
    
11. clr.AddReferenceByPartialName("System.Windows.Forms")
    
12. clr.AddReferenceByPartialName("System.Drawing")
    
13. from System import *
    
14. from System.Windows.Forms import *
    
15. from System.Drawing import *
    
16. clr.AddReference("MeteoInfoC.dll")
    
17. from MeteoInfoC import *
    
18. from MeteoInfoC.Data.MapData import *
    
19. from MeteoInfoC.Data.MeteoData import *
    
20. from MeteoInfoC.Layout import *
    
21. from MeteoInfoC.Projections import *
    
22. 
    
23. #---- 设置路径变量
    
24. baseDir = "C:\\Program Files (x86)\\MeteoInfo\\"
    
25. mapDir = baseDir + "Map\\"
    
26. legendDir = baseDir + "pal\\"
    
27. dataDir = "D:\\Temp\\micaps\\"
    
28. 
    
29. #---- 打开图层
    
30. bou2Layer = MapDataManage.OpenLayer(mapDir + "bou2_4p.shp")
    
31. lb = bou2Layer.LegendScheme.LegendBreaks[0]
    
32. lb.Color = Color.Yellow
    
33. lb.OutlineColor = Color.Gray
    
34. lb.OutlineSize = 1
    
35. lb.DrawFill = False
    
36. bou1Layer = MapDataManage.OpenLayer(mapDir + "bou1_4l.shp")
    
37. lb = bou1Layer.LegendScheme.LegendBreaks[0]
    
38. lb.Color = Color.Red
    
39. res1Layer = MapDataManage.OpenLayer(mapDir + "res1_4m.shp")
    
40. lb = res1Layer.LegendScheme.LegendBreaks[0]
    
41. lb.Color = Color.Red
    
42. lb.OutlineColor = Color.Black
    
43. lb.Size = 5
    
44. res1Layer.LabelSet.FieldName = "NAME";
    
45. res1Layer.LabelSet.LabelFont = Font("楷体", 10);
    
46. res1Layer.LabelSet.YOffset = 15;
    
47. res1Layer.LabelSet.LabelColor = Color.Red
    
48. res1Layer.AddLabels();
    
49. 
    
50. #---- 创建MIApp类的对象
    
51. myApp = MIApp()
    
52. mapLayout = myApp.MapLayout
    
53. layoutMap = mapLayout.ActiveLayoutMap
    
54. mapFrame = layoutMap.MapFrame
    
55. mapView = mapFrame.MapView
    
56. 
    
57. #---- 设置Layout
    
58. mapLayout.SetPaperSize(540, 550)
    
59. 
    
60. #---- 添加地图图层
    
61. mapFrame.AddLayer(bou2Layer)
    
62. mapFrame.AddLayer(bou1Layer)
    
63. mapFrame.AddLayer(res1Layer)
    
64. 
    
65. #---- 设置ActiveLayoutMap（图层显示）
    
66. layoutMap.DrawGridLine = False
    
67. layoutMap.DrawNeatLine = False
    
68. layoutMap.DrawGridLabel = False
    
69. layoutMap.DrawGridTickLine = False
    
70. layoutMap.Left = 10
    
71. layoutMap.Top = 10
    
72. layoutMap.Width = 500
    
73. layoutMap.Height = 490
    
74. 
    
75. #---- 设置图形标题
    
76. aText = mapLayout.AddText("Temp",250,50,"黑体",12)
    
77. aText.Graphic.Legend.Color = Color.Red
    
78. 
    
79. #---- 设置时间
    
80. sTime = DateTime.Parse("2013-08-22 16:00")
    
81. 
    
82. #---- 设置图例
    
83. palfn = legendDir + "I-01.pal"
    
84. 
    
85. #---- 设置MeteoDataInfo
    
86. mid = MeteoDataInfo()
    
87. 
    
88. #---- 设置数据文件  
    
89. inFile = sTime.ToString("yyMMddHH") + "00.000"
    
90. print inFile
    
91. #---- 打开MICAPS数据文件
    
92. mid.OpenMICAPSData(dataDir + inFile)  
    
93. #---- 投影
    
94. mapView.ProjectLayers(mid.ProjInfo)
    
95. #---- 获取云图格点数据
    
96. gData = mid.GetGridData("var")
    
97. #---- 生成云图图层
    
98. cloudLayer = DrawMeteoData.CreateRasterLayer(gData, "Cloud", palfn)
    
99. cloudLayer.ProjInfo = mid.ProjInfo
    
100. #---- 添加图层
     
101. mapFrame.AddLayer(cloudLayer)
     
102. print "Display finished"  
     
103. #---- 调整图层顺序（以避免压盖）
     
104. mapFrame.MoveLayerToBottom(cloudLayer)  
     
105. #---- 缩放至云图范围
     
106. mapView.ZoomToExtentEx(cloudLayer.Extent)
     
107. #---- 设置标题名称
     
108. title = "可见光云图(" + sTime.ToString("yyyy-MM-dd HH:00") + ")"
     
109. aText.SetLabelText(title)
     
110. #---- 绘制图形
     
111. mapLayout.PaintGraphics()
     
112. #---- 输出图形为文件
     
113. outFile = "Cloud_" + sTime.ToString("yyyyMMddHH") + ".png"
     
114. print outFile
     
115. #mapLayout.ExportToPicture("E:\\Temp\\" + outFile)  
     
116. 
     
117. #---- 显示程序窗体（只是为了看效果，自动运行时不需要）
     
118. #Application.Run(myApp)
     
119. myApp.Show()

# This Python file uses the following encoding: GB2312<br /> #--------------------------------------------------------        <br /> # Author: Yaqiang Wang                                           <br /> # Date: 2013-8-27                                                 <br /> # Purpose: Read and plot MICAPS cloud data  <br /> # Note: Sample                                                   <br /> #-----------------------------------------------------------  <br /> #---- 引入类库<br /> import clr<br /> clr.AddReferenceByPartialName("System")<br /> clr.AddReferenceByPartialName("System.Windows.Forms")<br /> clr.AddReferenceByPartialName("System.Drawing")<br /> from System import *<br /> from System.Windows.Forms import *<br /> from System.Drawing import *<br /> clr.AddReference("MeteoInfoC.dll")<br /> from MeteoInfoC import *<br /> from MeteoInfoC.Data.MapData import *<br /> from MeteoInfoC.Data.MeteoData import *<br /> from MeteoInfoC.Layout import *<br /> from MeteoInfoC.Projections import *<br /> <br /> #---- 设置路径变量<br /> baseDir = "C:\\Program Files (x86)\\MeteoInfo\\"<br /> mapDir = baseDir + "Map\\"<br /> legendDir = baseDir + "pal\\"<br /> dataDir = "D:\\Temp\\micaps\\"<br /> <br /> #---- 打开图层<br /> bou2Layer = MapDataManage.OpenLayer(mapDir + "bou2_4p.shp")<br /> lb = bou2Layer.LegendScheme.LegendBreaks[0]<br /> lb.Color = Color.Yellow<br /> lb.OutlineColor = Color.Gray<br /> lb.OutlineSize = 1<br /> lb.DrawFill = False<br /> bou1Layer = MapDataManage.OpenLayer(mapDir + "bou1_4l.shp")<br /> lb = bou1Layer.LegendScheme.LegendBreaks[0]<br /> lb.Color = Color.Red<br /> res1Layer = MapDataManage.OpenLayer(mapDir + "res1_4m.shp")<br /> lb = res1Layer.LegendScheme.LegendBreaks[0]<br /> lb.Color = Color.Red<br /> lb.OutlineColor = Color.Black<br /> lb.Size = 5<br /> res1Layer.LabelSet.FieldName = "NAME";<br /> res1Layer.LabelSet.LabelFont = Font("楷体", 10);<br /> res1Layer.LabelSet.YOffset = 15;<br /> res1Layer.LabelSet.LabelColor = Color.Red<br /> res1Layer.AddLabels();<br /> <br /> #---- 创建MIApp类的对象<br /> myApp = MIApp()<br /> mapLayout = myApp.MapLayout<br /> layoutMap = mapLayout.ActiveLayoutMap<br /> mapFrame = layoutMap.MapFrame<br /> mapView = mapFrame.MapView<br /> <br /> #---- 设置Layout<br /> mapLayout.SetPaperSize(540, 550)<br /> <br /> #---- 添加地图图层<br /> mapFrame.AddLayer(bou2Layer)<br /> mapFrame.AddLayer(bou1Layer)<br /> mapFrame.AddLayer(res1Layer)<br /> <br /> #---- 设置ActiveLayoutMap（图层显示）<br /> layoutMap.DrawGridLine = False<br /> layoutMap.DrawNeatLine = False<br /> layoutMap.DrawGridLabel = False<br /> layoutMap.DrawGridTickLine = False<br /> layoutMap.Left = 10<br /> layoutMap.Top = 10<br /> layoutMap.Width = 500<br /> layoutMap.Height = 490<br /> <br /> #---- 设置图形标题<br /> aText = mapLayout.AddText("Temp",250,50,"黑体",12)<br /> aText.Graphic.Legend.Color = Color.Red<br /> <br /> #---- 设置时间<br /> sTime = DateTime.Parse("2013-08-22 16:00")<br /> <br /> #---- 设置图例<br /> palfn = legendDir + "I-01.pal"<br /> <br /> #---- 设置MeteoDataInfo<br /> mid = MeteoDataInfo()<br /> <br /> #---- 设置数据文件  <br /> inFile = sTime.ToString("yyMMddHH") + "00.000"<br /> print inFile<br /> #---- 打开MICAPS数据文件<br /> mid.OpenMICAPSData(dataDir + inFile)  <br /> #---- 投影<br /> mapView.ProjectLayers(mid.ProjInfo)<br /> #---- 获取云图格点数据<br /> gData = mid.GetGridData("var") <br /> #---- 生成云图图层<br /> cloudLayer = DrawMeteoData.CreateRasterLayer(gData, "Cloud", palfn)<br /> cloudLayer.ProjInfo = mid.ProjInfo<br /> #---- 添加图层<br /> mapFrame.AddLayer(cloudLayer)<br /> print "Display finished"  <br /> #---- 调整图层顺序（以避免压盖）<br /> mapFrame.MoveLayerToBottom(cloudLayer)  <br /> #---- 缩放至云图范围<br /> mapView.ZoomToExtentEx(cloudLayer.Extent)<br /> #---- 设置标题名称<br /> title = "可见光云图(" + sTime.ToString("yyyy-MM-dd HH:00") + ")"<br /> aText.SetLabelText(title) <br /> #---- 绘制图形<br /> mapLayout.PaintGraphics()<br /> #---- 输出图形为文件<br /> outFile = "Cloud_" + sTime.ToString("yyyyMMddHH") + ".png"<br /> print outFile<br /> #mapLayout.ExportToPicture("E:\\Temp\\" + outFile)  <br /> <br /> #---- 显示程序窗体（只是为了看效果，自动运行时不需要）<br /> #Application.Run(myApp)<br /> myApp.Show()

运行结果：

才仔

大神威武

La_respiration

非常感谢老师

Role

好东西 感谢共享

liang

学习一下，收藏备用，多谢分享

寻路人

大神威武

huangfucheng

老师 请问下 脚本怎么执行呢

墨墨卡托

楼主可以分享个MeteoInfo读取、生成FY—2C的可见光、红外、水汽通道的代码吗？

MeteoInfo

墨墨卡托 发表于 2015-6-18 11:23
楼主可以分享个MeteoInfo读取、生成FY—2C的可见光、红外、水汽通道的代码吗？

数据是什么？

墨墨卡托

MeteoInfo 发表于 2015-6-18 11:52
数据是什么？

老师您好，数据虽然不大，但也不小，具体数据我随后发您的邮箱，这里提前写过老师了~