脚本如何实现MICAPS第一类多要素的叠加显示？【已解决】

MeteoInfo

孤蓝et 发表于 2011-12-1 18:10

                               
登录/注册后可看大图

辛苦斑竹了。。。。
继续支持中。

下载最新的MeteoInfo 1.0.1版，脚本中加入：
myApp.SetDrawType("barb")
myApp.MeteoDataInfo.MeteoUVStr.IsUV = Flase
myApp.Display("WindDirection","WindSpeed")

孤蓝et

MeteoInfo 发表于 2011-12-1 22:08

                               
登录/注册后可看大图

下载最新的MeteoInfo 1.0.1版，脚本中加入：
myApp.SetDrawType("barb")
myApp.MeteoDataInfo.MeteoUVS ...

感谢版主啊，辛苦您了。问题解决。
成功实现Micaps第一类数据的多要素叠加。

孤蓝et

效果图如下，脚本程序同样附上。

1. # This Python file uses the following encoding: GB2312
   
2. #---- 引入类库
   
3. import clr
   
4. clr.AddReferenceByPartialName("System")
   
5. clr.AddReferenceByPartialName("System.Windows.Forms")
   
6. clr.AddReferenceByPartialName("System.Drawing")
   
7. from System import *
   
8. from System.Windows.Forms import *
   
9. from System.Drawing import *
   
10. clr.AddReference("MeteoInfoC.dll")
    
11. from MeteoInfoC import *
    
12. from MeteoInfoC.Layout import *
    
13. from MeteoInfoC.Data import *
    
14. 
    
15. 
    
16. 
    
17. 
    
18. 
    
19. 
    
20. #---- 创建MIApp类的对象
    
21. myApp = MIApp()
    
22. #---- 添加图层
    
23. myApp.OpenLayer("..\\Map\\bou2_4p.shp")
    
24. #---- 设置图层渲染
    
25. myApp.SetLegendBreak("bou2_4p.shp",0,Color.Yellow,Color.Gray,1,True,False,True)
    
26. #myApp.OpenLayer("..\\Map\\bou1_4l.shp")
    
27. #myApp.SetLegendBreak("bou1_4l.shp",0,Color.Blue)
    
28. myApp.OpenLayer("..\\Map\\china.shp")
    
29. myApp.SetLayerVisible("china.shp", False)
    
30. myApp.OpenLayer("..\\Map\\res1_4m.shp")
    
31. myApp.SetLegendBreak("res1_4m.shp",0,5,Color.Red,Color.Black,True,True,True)
    
32. #---- 获取矢量图层对象，并设置标注
    
33. aLayer = myApp.GetVectorLayer("res1_4m.shp")
    
34. aLayer.LabelSet.FieldName = "NAME";
    
35. aLayer.LabelSet.LabelFont = Font("楷体", 10);
    
36. aLayer.LabelSet.Offset = 15;
    
37. aLayer.AddLabels();
    
38. 
    
39. #---- Lambert投影
    
40. myApp.ProjectLayers("+proj=lcc+lat_1=25+lat_2=47+lon_0=84")
    
41. 
    
42. #---- 按照经纬度范围缩放地图
    
43. myApp.ZoomLonLat(76,94,35,50)
    
44. 
    
45. #---- 设置屏蔽图层（只绘制中国境内图形）
    
46. myApp.SetMaskout("china.shp")
    
47. 
    
48. #---- 站点数据插值为格点数据的设置
    
49. myApp.SetInterpolation(60,140,-20,60,160,160,"IDW_Radius",1,1)
    
50. #---- 设置图例文件
    
51. #myApp.SetLegendScheme("..\\Legend\\tem.lgs")
    
52. 
    
53. #---- 设置DefalutLayoutMap（图层显示）
    
54. myApp.MapLayout.DefaultLayoutMap.DrawGridLine = True
    
55. myApp.MapLayout.DefaultLayoutMap.DrawNeatLine = True
    
56. myApp.MapLayout.DefaultLayoutMap.Left = 10
    
57. myApp.MapLayout.DefaultLayoutMap.Top = 10
    
58. myApp.MapLayout.DefaultLayoutMap.Width = 620
    
59. myApp.MapLayout.DefaultLayoutMap.Height = 450
    
60. 
    
61. 
    
62. 
    
63. #---- 设置DefaultTitle（图形标题）
    
64. myApp.MapLayout.DefaultTitle.SetFont("黑体", 12)
    
65. myApp.MapLayout.DefaultTitle.Top = 30
    
66. myApp.MapLayout.DefaultTitle.Left = 80
    
67. 
    
68. #---- 设置DefaultLegend（图例）
    
69. myApp.MapLayout.DefaultLegend.LegendStyle = LegendStyleEnum.Normal               
    
70. myApp.MapLayout.DefaultLegend.DrawNeatLine = False
    
71. myApp.MapLayout.DefaultLegend.Title = "温度"
    
72. myApp.MapLayout.DefaultLegend.BackColor = Color.White
    
73. myApp.MapLayout.DefaultLegend.Font = Font("Arial", 8)
    
74. myApp.MapLayout.DefaultLegend.Left = 640
    
75. myApp.MapLayout.DefaultLegend.Top = 150
    
76. 
    
77. 
    
78. #---- 设置起始结束时间
    
79. sTime = DateTime.Parse("2011-8-3 18:00")
    
80. eTime = DateTime.Parse("2011-8-3 18:00")
    
81. 
    
82. 
    
83. 
    
84. #---- 循环
    
85. while sTime <= eTime:
    
86.   myApp.RemoveDataLayers()
    
87.   inFile = sTime.ToString("yyMMddHH") + ".000"
    
88.   print inFile
    
89.   #---- 打开MICAPS数据文件
    
90.   myApp.OpenMICAPSData("..\\Sample\\aws\" + inFile)  
    
91. 
    
92.   #---- 设置绘图类型为shaded
    
93.   myApp.SetDrawType("Shaded")
    
94.   myApp.Display("Temperature")
    
95. 
    
96.   myApp.SetDrawType("Barb")
    
97.   myApp.MeteoDataInfo.MeteoUVStr.IsUV = False
    
98.   myApp.Display("WindDirection","WindSpeed")
    
99. 
    
100.   print "Display finished"  
     
101.   #---- 调整图层顺序（以避免压盖）
     
102.   myApp.MoveLayerToTop("bou2_4p.shp")  
     
103.   #myApp.MoveLayerToTop("bou1_4l.shp")
     
104.   myApp.MoveLayerToTop("res1_4m.shp")
     
105.   #---- 设置标题名称
     
106.   title = "                     新疆区域温度分布图 (" + sTime.ToString("yyyy-MM-dd HH:00") + ")"
     
107.   myApp.SetTitle(title)     
     
108.   #---- 设置图例名称
     
109.   myApp.MapLayout.DefaultLegend.Title = "温度(℃)"
     
110.   #---- 绘制图形
     
111.   myApp.MapLayout.PaintGraphics()
     
112.   #---- 输出图形为文件
     
113.   outFile = "tem_" + sTime.ToString("yyyyMMddHH") + ".png"
     
114.   print outFile
     
115.   myApp.SaveFigure("E:\" + outFile)  
     
116.   #---- 删除6小时降水量图层
     
117.   if sTime < eTime:
     
118.           myApp.RemoveLastLayer()
     
119.   #---- 时间加6小时
     
120.   sTime = sTime.AddHours(1)
     
121. 
     
122. #---- 显示程序窗体（只是为了看效果，自动运行时不需要）
     
123. Application.Run(myApp)

复制代码

孤蓝et

    刚好需要用一个软件自动批量生成常用的天气图或者综合图，MICAPS3的做天气图的还是需要人工干预，用楼主的这个就无敌啦，做好脚本，自动执行出图，然后我再用程序批量得到天气讲解的PPT。完美了。

赞楼主一个。

喜上眉梢

楼主和版主真厉害，此方法我也想学用。

孤蓝et

喜上眉梢 发表于 2011-12-2 11:44

                               
登录/注册后可看大图

楼主和版主真厉害，此方法我也想学用。

我那VB写交互界面，自动生成、执行meteoinfo Script，批量得到图形文件，而后使用VB自动将图形、标题等信息导入PPT，而后提供给预报员。
现在只有个思路。等我花点时间，写好好给大家分享。

MeteoInfo

孤蓝et 发表于 2011-12-2 12:19

                               
登录/注册后可看大图

我那VB写交互界面，自动生成、执行meteoinfo Script，批量得到图形文件，而后使用VB自动将图形、标题等信 ...

既然是想自动运行，似乎没必要再写交互界面，用IronPython应该能完成所有工作