ArcGis 拓扑检查——缺顶点、悬挂检查代码 C#

看了些源码，效率挺垃圾的，折腾了一个垃圾得不太彻底的代码，还是慢。

不会折腾底层直接怼COM的悲伤……

实现思路是这样的：

1、把面层的点都塞进List，去重，取坐标4位，后边的检查使用容差0.001

2、遍历点，通过点在面层寻相交的面

3、如果结果是1，那么这个面在这个点处没有毗邻面，把点缓冲区一下给定距离，如果能找到面了，那么悬挂悬挂。

如果结果>1,那么遍历所有相交面，如果面的PointCollection里有这个点，那么计数+1；如果存在PointCollection里不包含这个点的面 ，那么缺顶点缺顶点

下面贴代码

取点集，去个重：

 1   class UserPoints
 2     {
 3         public static  List<IPoint> FeatureLayer2PointList(IFeatureLayer pFeatureLayer)
 4         {
 5             List<IPoint> pointList = new List<IPoint>();
 6             IFeatureCursor pFeatureCursor = pFeatureLayer.FeatureClass.Search(null, true);
 7             try
 8             {
 9                 IFeature pFeatuare = pFeatureCursor.NextFeature();
10                 while (pFeatuare != null)
11                 {
12                     IPointCollection pcol = pFeatuare.Shape as IPointCollection;
13                     for (int i = 0; i < pcol.PointCount - 1; i++)
14                     {
15                         pointList.Add(pcol.Point[i]);
16                     }
17                     pFeatuare = pFeatureCursor.NextFeature();
18                 }
19                 pointList = pointList.Distinct(new Compare()).ToList<IPoint>();
20             }
21             catch (Exception exp)
22             {
23                 ErrorF err = new ErrorF(exp.Message + "\r\n" + exp.StackTrace);
24                 err.Show();
25             }
26             finally
27             {
28                 System.Runtime.InteropServices.Marshal.ReleaseComObject(pFeatureCursor);             
29             }
30             return pointList;
31         }
32     }
33 
34     class Compare : IEqualityComparer<IPoint>
35     {
36 
37         bool IEqualityComparer<IPoint>.Equals(IPoint a, IPoint b)
38         {
39             if (a == null && b == null)
40                 return false;
41             else
42                 return Math.Round(a.X, 4) == Math.Round(b.X, 4) && Math.Round(a.Y, 4) == Math.Round(b.Y, 4);
43         }
44 
45         int IEqualityComparer<IPoint>.GetHashCode(IPoint obj)
46         {
47             return obj.ToString().GetHashCode();
48         }
49     }

View Code

拓扑一下：

 1      public static List<IPoint> CheckLackJunctionPointOrSuspendedPoint(IFeatureLayer pFeatureLayer, double distance)
 2         {
 3             IGraphicsContainer pGraphicsContainer = (IGraphicsContainer)m_hookHelper.FocusMap.ActiveGraphicsLayer;
 4             pGraphicsContainer.DeleteAllElements();
 5             IColor innerColor = new RgbColorClass();
 6             innerColor.NullColor = true;
 7             IColor outLineColor = DisplayUtils.RGBColor(255, 0, 0);
 8             IElement pElement;
 9             List<IPoint> listError = new List<IPoint>();
10             IFeatureCursor pFeatureCursor=null;
11             IFeature pFeature;
12             try
13             {
14                 IGeometry pGeometry;
15                 foreach (IPoint point in UserPoints.FeatureLayer2PointList(pFeatureLayer))
16                 {
17                     ITopologicalOperator pTopologicalOperator = point as ITopologicalOperator;
18                     ISpatialFilter pSpatialFilter = FilterUtil.SpatialFilter(point as IGeometry, esriSpatialRelEnum.esriSpatialRelIntersects);  
19                     int count = pFeatureLayer.FeatureClass.FeatureCount(pSpatialFilter);
20                     if (count == 1)
21                     {
22                         IGeometry pGeometryPointBuffer =pTopologicalOperator.Buffer(distance);
23                         pGeometry = pTopologicalOperator.Buffer(distance) ;
24                         pSpatialFilter = FilterUtil.SpatialFilter(pGeometry, esriSpatialRelEnum.esriSpatialRelIntersects);
25                         if (pFeatureLayer.FeatureClass.FeatureCount(pSpatialFilter) > 1)
26                         {
27                             pElement = DisplayUtils.CircleMarkElement(point, innerColor, outLineColor, 8.0);
28                             pGraphicsContainer.AddElement(pElement, 0);
29                             listError.Add(point);
30                         }
31                     }
32                     else if (count > 1)
33                     {
34                         pFeatureCursor = pFeatureLayer.FeatureClass.Search(pSpatialFilter, true);
35                         pFeature = pFeatureCursor.NextFeature();
36                         int count2 = 0;
37                         while (pFeature != null)
38                         {
39                             IPointCollection pPointCollection = pFeature.Shape as IPointCollection;
40                             IPoint pPointtemp = new PointClass();
41                             for (int k = 0; k < pPointCollection.PointCount - 1; k++)
42                             {
43                                 pPointCollection.QueryPoint(k, pPointtemp);
44                                 if (Math.Abs(pPointtemp.X - point.X) < 0.001 && Math.Abs(pPointtemp.Y - point.Y) < 0.001)
45                                 {
46                                     count2++;
47                                     break;
48                                 }
49                             }
50                             pFeature = pFeatureCursor.NextFeature();
51                         }
52                         if (count2 < count)
53                         {
54                             pElement = DisplayUtils.CircleMarkElement(point, innerColor, outLineColor, 8.0);
55                             pGraphicsContainer.AddElement(pElement, 0);
56                             listError.Add(point);
57                         }
58                     }
59                 }
60             }
61             catch (Exception exp)
62             {
63                 ErrorF err = new ErrorF(exp.Message + "\r\n" + exp.StackTrace);
64                 err.Show();
65             }
66             finally
67             {
68                 Marshal.FinalReleaseComObject(pFeatureCursor);
69             }
70             return listError;
71         }

View Code

 

哪位有高效率的代码，求侮辱！

转载于:https://www.cnblogs.com/yzhyingcool/p/10026859.html