HowTo Perform the spatial selection 'Share a line segment with' using ArcObjects
HowTo Perform the spatial selection 'Share a line segment with' using ArcObjects
Article ID: |
26528 |
Software: |
ArcGIS - ArcEditor 8.1, 8.1.2, 8.2, 8.3 ArcGIS - ArcInfo 8.1, 8.1.2, 8.2, 8.3 ArcGIS - ArcView 8.1, 8.1.2, 8.2, 8.3 |
Platforms: |
Windows NT 4.0, 2000, XP |
Summary
Instructions provided illustrate how to duplicate the functionality of ArcMap's spatial selection 'share a line segment with' through ArcObjects code.
Procedure
- Start ArcMap.
- Open the Visual Basic Editor. -show me-
- In the Project Explorer window, expand Project.mxt and select ArcMap Objects > ThisDocument. Right-click and select .-show me-
Code in the Project's ThisDocument code module will only run in the current map document. If you want to store the code in all map documents open the Normal.mxt ThisDocument code module.
- Copy/Paste the following code into the code module.
- Option Explicit
- Public Sub SelectPolyAsSharesSegementWith()
- Dim pMxDoc As IMxDocument
- Set pMxDoc = ThisDocument
- Dim pMap As IMap
- Set pMap = pMxDoc.FocusMap
- ' Union all of the search features into one geometry and strip off
- ' it's boundary for use in the search. Since we are only interested
- ' in shared segments, not areas, the boundary is all that is needed.
- ' Get search layer
- Dim pPolygonLyr As IFeatureLayer
- Set pPolygonLyr = pMap.Layer(1)
- Dim pLFeatClass As IFeatureClass
- Set pLFeatClass = pPolygonLyr.FeatureClass
- ' Loop through the polygons, simplify them and then
- ' add them to a Geometry Bag.
- Dim pFeatureCursor As IFeatureCursor
- Set pFeatureCursor = pLFeatClass.Search(Nothing, True)
- Dim pFeature As IFeature
- Set pFeature = pFeatureCursor.NextFeature
- Dim pGeometryCollection As IGeometryCollection
- Set pGeometryCollection = New GeometryBag
- Dim pTopologicalOperator2 As ITopologicalOperator2
- Do While Not pFeature Is Nothing
- Set pTopologicalOperator2 = pFeature.ShapeCopy
- pTopologicalOperator2.IsKnownSimple = False
- pTopologicalOperator2.Simplify
- pGeometryCollection.AddGeometry pTopologicalOperator2
- Set pFeature = pFeatureCursor.NextFeature
- Loop
- ' Union all the polygons
- Dim pTopologicalOperator1 As ITopologicalOperator2
- Set pTopologicalOperator1 = New Polygon
- pTopologicalOperator1.ConstructUnion pGeometryCollection
- ' Construct the boundary of the unioned polygons
- Dim pTopologicalOperator3 As ITopologicalOperator2
- Set pTopologicalOperator3 = New Polyline
- Set pTopologicalOperator3 = pTopologicalOperator1.Boundary
- ' Use the boundary to find the candidate features
- ' Get the target FeatureClass
- Dim pFeatureLayer As IFeatureLayer
- Set pFeatureLayer = pMap.Layer(0)
- Dim pFeatureClass As IFeatureClass
- Set pFeatureClass = pFeatureLayer.FeatureClass
- ' Set up a scratch workspace to use for the selection
- Dim pScratchWorkspace As IWorkspace
- Dim pScratchWorkspaceFactory As IScratchWorkspaceFactory
- Set pScratchWorkspaceFactory = New ScratchWorkspaceFactory
- Set pScratchWorkspace = pScratchWorkspaceFactory.DefaultScratchWorkspace
- ' Setup the spatial filter
- Dim pSpatialFilter As ISpatialFilter
- Set pSpatialFilter = New SpatialFilter
- With pSpatialFilter
- Set .Geometry = pTopologicalOperator3
- .GeometryField = "SHAPE"
- .SpatialRel = esriSpatialRelIntersects
- End With
- ' Perform the search and create a selection set
- Dim pTempSelection As ISelectionSet
- Set pTempSelection = pFeatureClass.Select(pSpatialFilter, esriSelectionTypeIDSet, esriSelectionOptionNormal, pScratchWorkspace)
- ' Loop through the features and check to see if they share a segment
- ' with the search geometry. If they don't, remove them from the selection
- ' set.
- ' Create a feature cursor to allow us the examine the candidate features
- Dim pTempFeatCur As IFeatureCursor
- pTempSelection.Search Nothing, True, pTempFeatCur
- Dim iOIDList(1) As Long
- Dim iOIDListCount As Long
- iOIDListCount = 1
- Dim pGeometry As IPolyline
- Dim pTempFeat As IFeature
- ' Fetch the first feature
- Set pTempFeat = pTempFeatCur.NextFeature
- Dim pTempPolyline As IPolyline
- ' Check for the one dimentional intersection of the geometries.
- ' If there is no intersection the returned geometry will be empty. If
- ' it is, remove it from the selection set.
- Dim pTopologicalOperator4 As ITopologicalOperator2
- Set pTopologicalOperator4 = New Polygon
- Do
- iOIDList(0) = pTempFeat.OID
- Set pTopologicalOperator4 = pTempFeat.ShapeCopy
- pTopologicalOperator4.IsKnownSimple = False
- pTopologicalOperator4.Simplify
- Set pTempPolyline = pTopologicalOperator4.Boundary
- Set pGeometry = pTopologicalOperator3.Intersect(pTempPolyline, esriGeometry1Dimension)
- If pGeometry.IsEmpty Then
- pTempSelection.RemoveList iOIDListCount, iOIDList(0)
- End If
- Set pTempFeat = pTempFeatCur.NextFeature
- Loop Until pTempFeat Is Nothing
- ' Set the feature selection
- Dim pFeatureSelection As IFeatureSelection
- Set pFeatureSelection = pFeatureLayer
- Set pFeatureSelection.SelectionSet = pTempSelection
- ' Display the selected features.
- Dim pdoc As IMxDocument
- Set pdoc = ThisDocument
- pdoc.ActiveView.Refresh
End Sub
- Close the Visual Basic Editor.
Created: 4/1/2004
Last Modified: 7/1/2008
Article Rating: (1)
If you would like to post a comment, please login
Comments
By Anonymous - 07/20/2005 7:47 AM
Hello, The sample code to select polygons that share a segment with the selected polygon did not work for me. Although I tested it in 9.1, no library is explicitly defined so I don't think this is the problem. http://support.esri.com/index.cfm?fa=knowledgebase.techarticles.articleShow&d=26528 Sincerely, Valentina
HowTo Perform the spatial selection 'Share a line segment with' using ArcObjects的更多相关文章
- Division of Line Segment
Division of Line Segment /** */ void Line::EqualDivision(int nCount, QLineF fline, QList<QPointF& ...
- 目标检测之线段检测---lsd line segment detector
(1)线段检测应用背景 (2)线段检测原理简介 (3)线段检测实例 a line segment detector (4)hough 变换和 lsd 的区别 --------------------- ...
- HDU 1542/POJ 1151 Atlantis (scaning line + segment tree)
A template of discretization + scaning line + segment tree. It's easy to understand, but a little di ...
- Line segment matching
FMII2方法:FMII方法的轻微的修改.有限线段和无限线段(直线)的匹配. 求解方法: SVD分解 Unit Quaternion 协方差矩阵: 通过对C进行SVD分解得到R,根据R求得T. 算法流 ...
- Android OpenGL ES(九)绘制线段Line Segment .
创建一个DrawLine Activity,定义四个顶点: float vertexArray[] = { -0.8f, -0.4f * 1.732f, 0.0f, -0.4f, 0.4f * 1.7 ...
- 用ArcGIS?37个Arcmap常用操作技巧可能帮到您
1. 要素的剪切与延伸 实用工具 TASK 任务栏 Extend/Trim feature 剪切所得内容与你画线的方向有关. 2. 自动捕捉跟踪工具 点击Editor工具栏中Snapping来打开Sn ...
- ArcGIS中的查询
最近身体不适,静下心来看了一下以前收集的电子书.下面是<ArcGIS地理信息系统教程_第5版>(李玉龙)第5章“查询”的读书笔记. 1.查询的常见应用: 选择感兴趣的要素:查找哪些要素满足 ...
- 2016-2017 National Taiwan University World Final Team Selection Contest C - Crazy Dreamoon
题目:Statements Dreamoon likes algorithm competitions very much. But when he feels crazy because he ca ...
- 2016-2017 National Taiwan University World Final Team Selection Contest A - Hacker Cups and Balls
题目: Dreamoon likes algorithm competitions very much. But when he feels crazy because he cannot figur ...
随机推荐
- PAT (Basic Level) Practise:1023. 组个最小数
[题目链接] 给定数字0-9各若干个.你可以以任意顺序排列这些数字,但必须全部使用.目标是使得最后得到的数尽可能小(注意0不能做首位).例如:给定两个0,两个1,三个5,一个8,我们得到的最小的数就是 ...
- svn Q&A
Q1:在svn commit的时候,会出现某某文件 is missing.这是因为此次提交时:远程repository中并没有该文件,而且本地repository也没有该文件. 具体原因: 1.可能因 ...
- 16年青岛网络赛 1002 Cure
题目链接:http://acm.hdu.edu.cn/contests/contest_showproblem.php?pid=1002&cid=723 Cure Time Limit: 30 ...
- vim编辑器使用相关
alias 查看vi是否已经绑定vim 一.vim的块选择 v 字符选择 V 行选择 ctrl+v 快选择 y 复制选择的地方(p进行粘贴) d 删除选择的地方 二.vim多文件编辑 :n 编辑下一个 ...
- H20的题——[noip2003]银河英雄传(并查集)
公元五八○一年,地球居民迁移至金牛座α第二行星,在那里发表银河联邦创立宣言,同年改元为宇宙历元年,并开始向银河系深处拓展. 宇宙历七九九年,银河系的两大军事集团在巴米利恩星域爆发战争.泰山压顶集团派宇 ...
- python string
string比较连接 >>> s1="python string" >>> len(s) 13 >>> s2=" p ...
- 转:Elasticsearch TermQuery 详解
JavaClient 查询ES QueryBuilder queryBuilder = QueryBuilders.termQuery("字段","term值" ...
- make自动生成依赖文件的两种形式
最近编译源文件发现当修改头文件后,make并不会自动把包含此头文件的所有源文件重新编译,而每次都是需要把对应的中间文件清除才能重新编译,非常的麻烦.因此需要make自动对源文件所依赖的头文件进行管理, ...
- 【AT91SAM3S】英蓓特EM-SAM3S开发板例子工程中的启动文件分析
手上一块英倍特的EM-SAM3S开发板,拿到已经有一个月了.本来是做uLoong活动使用的板子,可当初由于不熟悉这个芯片,使用了STM32F4当作了替代.最近准备抽点时间折腾下这个板子. 这个板子的资 ...
- Semaphore tryAcquire release 正确的使用方法
boolean permit = false; try { permit = semaphore.tryAcquire(1, TimeUnit.SECONDS); if (permit) { Syst ...