MeteoInfo脚本示例：绘制同心圆

MeteoInfo

之前示例了用脚本进行一些图形的绘制：MeteoInfo脚本示例：MeteoInfo脚本示例：图形绘制http://bbs.06climate.com/forum.php?mod=viewthread&tid=10278&fromuid=106

这个帖子补充讲一下同心圆图形的绘制。圆的绘制用到了CircleShape类，该类继承自PolygonShape类，可以用圆心点和半径来构造CircleShape。通过圆心不同角度的线段可以用PolylineShape来绘制。需要MeteoInfo最新文件（见置顶帖子）。

比较简单，多看看代码应该能理解，运行效果如下（注意MapView的XYScaleFactor属性要设为1，缺省是1.2）：


脚本代码如下：

1. #--------------------------------------------------------        
   
2. # Author: Yaqiang Wang                                          
   
3. # Date: 2013-7-29                                                
   
4. # Purpose: Add cicle graphics  
   
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. from math import *
    
24. 
    
25. #---- Create circle graphics
    
26. print 'Create circle graphics...'
    
27. cpoint = PointD(120.6, 30.6)
    
28. radius = [0.5, 1, 1.5]
    
29. pLegend = PolygonBreak()
    
30. pLegend.DrawFill = False
    
31. pLegend.OutlineColor = Color.Black
    
32. 
    
33. for r in radius:
    
34.         cshape = CircleShape(cpoint, r)
    
35.         cgraphic = Graphic(cshape, pLegend)
    
36.         mipy.MapDocument.ActiveMapFrame.MapView.GraphicCollection.Add(cgraphic)
    
37. 
    
38. #---- Create polyline graphics
    
39. print 'Create polyline graphics...'
    
40. def getcirclepoint(cp, r, a):
    
41.         x = cp.X + r * cos(a * 3.1416 / 180.0)
    
42.         y = cp.Y + r * sin(a * 3.1416 / 180.0)
    
43.         p = PointD(x, y)
    
44.         return p
    
45. 
    
46. llegend = PolyLineBreak()
    
47. llegend.Color = Color.Black
    
48. r = 1.5
    
49. plist = List[PointD]()
    
50. plist.Add(PointD(cpoint.X - r, cpoint.Y))
    
51. plist.Add(PointD(cpoint.X + r, cpoint.Y))
    
52. lshape = PolylineShape()
    
53. lshape.Points = plist
    
54. lgraphic = Graphic(lshape, llegend)
    
55. mipy.MapDocument.ActiveMapFrame.MapView.GraphicCollection.Add(lgraphic)
    
56. plist = List[PointD]()
    
57. plist.Add(PointD(cpoint.X, cpoint.Y + r))
    
58. plist.Add(PointD(cpoint.X, cpoint.Y - r))
    
59. lshape = PolylineShape()
    
60. lshape.Points = plist
    
61. lgraphic = Graphic(lshape, llegend)
    
62. mipy.MapDocument.ActiveMapFrame.MapView.GraphicCollection.Add(lgraphic)
    
63. plist = List[PointD]()
    
64. angle = 45.0
    
65. p = getcirclepoint(cpoint, r, angle)
    
66. plist.Add(p)
    
67. angle = 225.0
    
68. p = getcirclepoint(cpoint, r, angle)
    
69. plist.Add(p)
    
70. lshape = PolylineShape()
    
71. lshape.Points = plist
    
72. lgraphic = Graphic(lshape, llegend)
    
73. mipy.MapDocument.ActiveMapFrame.MapView.GraphicCollection.Add(lgraphic)
    
74. plist = List[PointD]()
    
75. angle = 135.0
    
76. p = getcirclepoint(cpoint, r, angle)
    
77. plist.Add(p)
    
78. angle = 315.0
    
79. p = getcirclepoint(cpoint, r, angle)
    
80. plist.Add(p)
    
81. lshape = PolylineShape()
    
82. lshape.Points = plist
    
83. lgraphic = Graphic(lshape, llegend)
    
84. mipy.MapDocument.ActiveMapFrame.MapView.GraphicCollection.Add(lgraphic)
    
85. 
    
86. #---- Refresh MapView
    
87. mipy.MapDocument.ActiveMapFrame.MapView.PaintLayers()
    
88. print 'Finished!'

#--------------------------------------------------------        <br /> # Author: Yaqiang Wang                                           <br /> # Date: 2013-7-29                                                 <br /> # Purpose: Add cicle graphics  <br /> # Note: Sample                                                   <br /> #-----------------------------------------------------------  <br /> #---- Import clr and classes<br /> import clr<br /> clr.AddReferenceByPartialName("System")<br /> clr.AddReferenceByPartialName("System.Drawing")<br /> from System import *<br /> from System.Drawing import *<br /> from System.Collections.Generic import *<br /> clr.AddReference("MeteoInfoC.dll")<br /> from MeteoInfoC import *<br /> from MeteoInfoC.Data.MapData import *<br /> from MeteoInfoC.Geoprocess import *<br /> from MeteoInfoC.Layer import *<br /> from MeteoInfoC.Shape import *<br /> from MeteoInfoC.Legend import *<br /> from MeteoInfoC.Drawing import *<br /> import os.path<br /> from math import *<br /> <br /> #---- Create circle graphics<br /> print 'Create circle graphics...'<br /> cpoint = PointD(120.6, 30.6)<br /> radius = [0.5, 1, 1.5]<br /> pLegend = PolygonBreak()<br /> pLegend.DrawFill = False<br /> pLegend.OutlineColor = Color.Black<br /> <br /> for r in radius:<br />         cshape = CircleShape(cpoint, r)<br />         cgraphic = Graphic(cshape, pLegend)<br />         mipy.MapDocument.ActiveMapFrame.MapView.GraphicCollection.Add(cgraphic)<br /> <br /> #---- Create polyline graphics<br /> print 'Create polyline graphics...'<br /> def getcirclepoint(cp, r, a):<br />         x = cp.X + r * cos(a * 3.1416 / 180.0)<br />         y = cp.Y + r * sin(a * 3.1416 / 180.0)<br />         p = PointD(x, y)<br />         return p<br /> <br /> llegend = PolyLineBreak()<br /> llegend.Color = Color.Black<br /> r = 1.5<br /> plist = List[PointD]()<br /> plist.Add(PointD(cpoint.X - r, cpoint.Y))<br /> plist.Add(PointD(cpoint.X + r, cpoint.Y))<br /> lshape = PolylineShape()<br /> lshape.Points = plist<br /> lgraphic = Graphic(lshape, llegend)<br /> mipy.MapDocument.ActiveMapFrame.MapView.GraphicCollection.Add(lgraphic)<br /> plist = List[PointD]()<br /> plist.Add(PointD(cpoint.X, cpoint.Y + r))<br /> plist.Add(PointD(cpoint.X, cpoint.Y - r))<br /> lshape = PolylineShape()<br /> lshape.Points = plist<br /> lgraphic = Graphic(lshape, llegend)<br /> mipy.MapDocument.ActiveMapFrame.MapView.GraphicCollection.Add(lgraphic)<br /> plist = List[PointD]()<br /> angle = 45.0<br /> p = getcirclepoint(cpoint, r, angle)<br /> plist.Add(p)<br /> angle = 225.0<br /> p = getcirclepoint(cpoint, r, angle)<br /> plist.Add(p)<br /> lshape = PolylineShape()<br /> lshape.Points = plist<br /> lgraphic = Graphic(lshape, llegend)<br /> mipy.MapDocument.ActiveMapFrame.MapView.GraphicCollection.Add(lgraphic)<br /> plist = List[PointD]()<br /> angle = 135.0<br /> p = getcirclepoint(cpoint, r, angle)<br /> plist.Add(p)<br /> angle = 315.0<br /> p = getcirclepoint(cpoint, r, angle)<br /> plist.Add(p)<br /> lshape = PolylineShape()<br /> lshape.Points = plist<br /> lgraphic = Graphic(lshape, llegend)<br /> mipy.MapDocument.ActiveMapFrame.MapView.GraphicCollection.Add(lgraphic)<br /> <br /> #---- Refresh MapView<br /> mipy.MapDocument.ActiveMapFrame.MapView.PaintLayers()<br /> print 'Finished!'

又是那隻貓

谢谢王老师分享和指导

wbriver

王老师辛苦了，将脚本文件翻译成C#语言了，交作业了~

MeteoInfo

wbriver 发表于 2013-7-29 23:14

                               
登录/注册后可看大图

王老师辛苦了，将脚本文件翻译成C#语言了，交作业了~

很强悍实用的作业！

雷达数据是什么格式的？

wbriver

MeteoInfo 发表于 2013-7-29 23:19

                               
登录/注册后可看大图

很强悍实用的作业！

雷达数据是什么格式的？

王老师过奖了~转成micap4类再显示的，回波的空间位置和实际位置还有偏差，还要再调整~

♂雨已~○

wbriver 发表于 2013-7-29 23:26

                               
登录/注册后可看大图

王老师过奖了~转成micap4类再显示的，回波的空间位置和实际位置还有偏差，还要再调整~

哥们能不能共享一哈转换成micaps类数据的程序算法或者源代码呢？

wbriver

♂雨已~○ 发表于 2013-8-6 21:41

                               
登录/注册后可看大图

恩，谢谢，这个我也查到了这篇文献，我想问问你那个转成micaps格点数据是咋个插值的呢？分辨率多大呢？给 ...

双线性差值就够了，按照你的需求，个人觉得想插细了可以考虑1km*1km，不过数据量比较大，用shaded图层显示会比较慢，raster快多了，不过效果不是很好啊~

♂雨已~○

wbriver 发表于 2013-8-6 21:59

                               
登录/注册后可看大图

双线性差值就够了，按照你的需求，个人觉得想插细了可以考虑1km*1km，不过数据量比较大，用shaded图层显示 ...

恩，谢咯哈

♂雨已~○

wbriver 发表于 2013-8-6 14:45

                               
登录/注册后可看大图

好久没来了。。。我这里直接用的省里转好的micaps 4类数据~
其实你只要根据雷达站经纬计算出起止经纬度， ...

恩，谢谢，这个我也查到了这篇文献，我想问问你那个转成micaps格点数据是咋个插值的呢？分辨率多大呢？给个思路呢，兄台，然后出图的时候用的是shaded图层还是raster图层呢？

wbriver

♂雨已~○ 发表于 2013-7-30 18:11

                               
登录/注册后可看大图

哥们能不能共享一哈转换成micaps类数据的程序算法或者源代码呢？

好久没来了。。。我这里直接用的省里转好的micaps 4类数据~
其实你只要根据雷达站经纬计算出起止经纬度，以及径向和纬向的格点数就可以转了，网上有相关文献的~