MeteoInfo脚本示例：创建点图层并判断点是否在某区域内

MeteoInfo

借用坛友的问题：如何判断点是否在某区域内？这里提供一个更完整的脚本程序来生成一个等间隔网格点图层，判断每个点是否在青藏高原内并在点图层中赋相应的属性（'Y' or 'N'），输出内部的点的坐标，设置点图层的LegendScheme（区域内的点红色，外面的点绿色）。可以很直观的验证判断结果。

运行结果：




脚本代码：

1. #---- Import clr and classes
   
2. import clr
   
3. clr.AddReferenceByPartialName("System")
   
4. clr.AddReferenceByPartialName("System.Drawing")
   
5. from System import *
   
6. from System.Drawing import *
   
7. clr.AddReference("MeteoInfoC.dll")
   
8. from MeteoInfoC import *
   
9. from MeteoInfoC.Data.MapData import *
   
10. from MeteoInfoC.Geoprocess import *
    
11. from MeteoInfoC.Layer import *
    
12. from MeteoInfoC.Shape import *
    
13. from MeteoInfoC.Legend import *
    
14. 
    
15. #---- Open layer
    
16. mapdir = "D:\\MyProgram\\CSharp\\MeteoInfoDev\\MeteoInfo\\bin\\Debug\\Map\"
    
17. aFile = mapdir + "tibet_ASCII.wmp"
    
18. aLayer = MapDataManage.OpenLayer(aFile)
    
19. 
    
20. #---- Add layer
    
21. mipy.MapDocument.ActiveMapFrame.AddLayer(aLayer)
    
22. mipy.MapDocument.ActiveMapFrame.MoveLayer(aLayer.Handle, 0)
    
23. 
    
24. #---- Create a grid point and determine grid points are in or out of Tibet
    
25. nLayer = VectorLayer(ShapeTypes.Point)
    
26. nLayer.LayerName = "New_Point_Layer"
    
27. nLayer.Visible = True
    
28. fieldName = "Inside"
    
29. nLayer.EditAddField(fieldName, String)
    
30. 
    
31. print 'The points inside Tibet are:'
    
32. tibetArea = aLayer.ShapeList[0]
    
33. sLon = 62.0
    
34. sLat = 25.0
    
35. delt = 1.0
    
36. for i in range(0, 45):
    
37.   lon = sLon + i * delt
    
38.   for j in range(0, 24):
    
39.     aPS = PointShape()   
    
40.     lat = sLat + j * delt
    
41.     aPoint = PointD(lon, lat)
    
42.     aPS.Point = aPoint
    
43.     shapeNum = nLayer.ShapeNum
    
44.     print "Shape Number: %d" % (shapeNum)
    
45.     nLayer.EditInsertShape(aPS, shapeNum)
    
46.     if GeoComputation.PointInPolygon(tibetArea, aPoint):
    
47.       nLayer.EditCellValue(fieldName, shapeNum, "Y")
    
48.       print "Lon: %5.1f, Lat: %5.1f" % (lon, lat)
    
49.     else:
    
50.       nLayer.EditCellValue(fieldName, shapeNum, "N")
    
51. 
    
52. aLS = nLayer.CreateLegendScheme(LegendType.UniqueValue, fieldName)
    
53. aLS.LegendBreaks[0].DrawOutline = False
    
54. aLS.LegendBreaks[0].Color = Color.Green
    
55. aLS.LegendBreaks[1].DrawOutline = False
    
56. aLS.LegendBreaks[1].Color = Color.Red
    
57. nLayer.LegendScheme = aLS
    
58. print "Layer Name: %s" % (nLayer.LayerName)
    
59. mipy.MapDocument.ActiveMapFrame.AddLayer(nLayer)   
    
60. print 'Finished!'

复制代码

jeydragon

必须顶起来

wet510

楼主的脚本编辑器不错哦，能不能共享下？

MeteoInfo

wet510 发表于 2012-9-26 08:13

                               
登录/注册后可看大图

楼主的脚本编辑器不错哦，能不能共享下？

最新版的MeteoInfo软件里就有呀，Tools -> Script 菜单。

haobang008

这个算法的核心思想是？求赐教

流云飞袖

很好，很强大

MeteoInfo

流云飞袖 发表于 2012-12-23 23:21

                               
登录/注册后可看大图

点文件的格式什么？

这个程序里的格点是程序里自动生成的，不存在格式问题。

清晨的雨

MeteoInfo 发表于 2012-9-26 08:16
最新版的MeteoInfo软件里就有呀，Tools -> Script 菜单。

王老师，自己的二次开发里面怎么运行脚本程序呢？

MeteoInfo

清晨的雨 发表于 2014-12-26 10:49
王老师，自己的二次开发里面怎么运行脚本程序呢？

好像你用C#开发，查阅一下如何将IronPython嵌入到C#中的资料。

wn8891

楼主，我想问问怎么用MeteoInfo提取某一地区的边界数据啊？