删除重复的feature vba VS 删除重复的feature python

最新推荐文章于 2021-01-08 10:31:37 发布

weixin_34190136

最新推荐文章于 2021-01-08 10:31:37 发布

阅读量132

点赞数

CC 4.0 BY-SA版权

文章标签： python javascript ViewUI

原文链接：http://www.cnblogs.com/qiernonstop/p/3634358.html

本文对比了VBA与Python两种语言在ArcGIS环境中的使用情况，通过一个具体的案例——删除重复要素——展示了两种语言实现相同功能的不同方式。文章指出Python相比VBA在ArcGIS中更为简洁易用，并且具有更好的支持和发展前景。

VBA：

 1 Sub deleteDuplicatedFeature()
 2 
 3 Dim app As IApplication
 4 Set app = Application
 5 
 6 Dim pMxDocument As IMxDocument
 7 Set pMxDocument = Application.Document
 8 
 9 Dim pMap As IMap
10 Set pMap = pMxDocument.FocusMap
11 
12 Dim pFeatureLayer As IFeatureLayer
13 Set pFeatureLayer = pMap.Layer(0)
14 
15 
16 Dim pFeatureClass As IFeatureClass
17 Set pFeatureClass = pFeatureLayer.FeatureClass
18 
19 
20 Dim pFeatureCursor As IFeatureCursor
21 Set pFeatureCursor = pFeatureClass.Search(Nothing, False)
22 
23 Dim pArea As IArea
24 Dim pRelOp As IRelationalOperator
25 Dim pFeature As IFeature
26 Set pFeature = pFeatureCursor.NextFeature
27 
28 Dim findex As Integer
29 findex = pFeatureClass.FindField("ifCopyed")
30 
31 Dim fid As Integer
32 fid = pFeatureClass.FindField("FID")
33 
34 
35 Dim pFilter As ISpatialFilter
36 
37 Dim pFeatureCursorCompared As IFeatureCursor
38 Dim pFeatureCompared As IFeature
39 Dim pShapeCompared As IPolygon
40 Dim pAreaCompared As IArea
41 
42  
43 
44 While Not pFeature Is Nothing
45    
46    Set pRelOp = pFeature.Shape
47 
48    
49    Set pFilter = New SpatialFilter
50 
51    With pFilter
52     Set .Geometry = pFeature.Shape
53         .GeometryField = "SHAPE"
54         .SpatialRel = esriSpatialRelIntersects
55    End With
56   
57    Set pFeatureCursorCompared = pFeatureClass.Search(pFilter, False)
58    Set pFeatureCompared = pFeatureCursorCompared.NextFeature
59 
60 
61    
62    
63     While Not pFeatureCompared Is Nothing
64     
65         Set pShapeCompared = pFeatureCompared.Shape
66         
67         If pRelOp.Equals(pShapeCompared) And pFeature.Value(fid) <> pFeatureCompared.Value(fid) Then
68     
69             pFeatureCompared.Delete   '删除之后，查询可以可以自动调节
70             
71          End If
72    
73 
74          
75          Set pFeatureCompared = pFeatureCursorCompared.NextFeature
76     Wend
77 
78    Set pFeature = pFeatureCursor.NextFeature
79 Wend

Python：

#run this program at arcgis10.0 environment please

 1 import arcpy
 2 
 3 fc = "D:\Output.shp"  #change to your own shapefile
 4 
 5 desc = arcpy.Describe(fc)
 6 
 7 shapefieldname = desc.ShapeFieldName
 8 
 9 rows = arcpy.SearchCursor(fc)
10 
11 for row in rows:
12     rows2 = arcpy.UpdateCursor(fc) 
13     for row2 in rows2:
14         if row.getValue(shapefieldname).equals(row2.getValue(shapefieldname)) and row.getValue("FID") != row2.getValue ("FID"):
15             rows2.deleteRow(row2)
16     del row2       
17 
18 del rows