ArcGis相接面补节点c#

相接（Touch）面执行切割后 新面与原相接面会缺少公共节点。

  private void AddPointToTouchesPolygon(IFeatureCursor newFeatureCursor, IFeatureClass featureClass)
         {
             IFeature newFeature = newFeatureCursor.NextFeature();
             while (newFeature != null)
             {
                 IPointCollection pointCollectionNewFeature = newFeature.Shape as IPointCollection;
                 for (int i = ; i < pointCollectionNewFeature.PointCount-; i++)
                 {
                     ITopologicalOperator pTopologicalOperator = pointCollectionNewFeature.Point[i] as IGeometry as ITopologicalOperator;
                     IGeometry pGeometryPointBuffer = pTopologicalOperator.Buffer(0.001);
                     ISpatialFilter spatialFilter0 = Utils.FilterUtils.SpatialFilter(pGeometryPointBuffer, esriSpatialRelEnum.esriSpatialRelIntersects);
                     IFeatureCursor featureCursor0 = featureClass.Update(spatialFilter0, false);
                     IFeature feature0 = featureCursor0.NextFeature();
                     while (feature0 != null)
                     {
                         IPointCollection pointCollection1 = feature0.Shape as IPointCollection;
                         IPoint pointTemp = new PointClass();
                         int count = ;
                         for (int k = ; k < pointCollection1.PointCount-; k++)
                         {
                             pointCollection1.QueryPoint(k, pointTemp);
                             if (Math.Abs(pointTemp.X - pointCollectionNewFeature.Point[i].X) < 0.001 && Math.Abs(pointTemp.Y - pointCollectionNewFeature.Point[i].Y) < 0.001)
                                 count++;
                         }
                         if (count == )
                         {
                             pointCollection1.AddPoint(pointCollectionNewFeature.Point[i]);
                             ITopologicalOperator topologicalOperator = feature0.Shape as ITopologicalOperator;
                             topologicalOperator.Simplify();
                             featureCursor0.UpdateFeature(feature0);
                             featureCursor0.Flush();
                         }
                         feature0 = featureCursor0.NextFeature();
                     }
                 }
                 newFeature = newFeatureCursor.NextFeature();
             }
         }

该方法只适用于规则的多边形，欢迎提出改进意见、bug及其解决

=============================================================================================================

更新方法：

      private void AddPointToTouchesPolygon(IFeatureCursor newFeatureCursor, IFeatureClass featureClass)
         {
             IFeature newFeature = newFeatureCursor.NextFeature();
             IFeatureCursor featureCursor0 = null;
             try
             {
                 while (newFeature != null)
                 {
                     IPointCollection pointCollectionNewFeature = newFeature.Shape as IPointCollection;
                     for (int i = ; i < pointCollectionNewFeature.PointCount - ; i++)
                     {
                         ITopologicalOperator pTopologicalOperator = pointCollectionNewFeature.Point[i] as IGeometry as ITopologicalOperator;
                         IGeometry pGeometryPointBuffer = pTopologicalOperator.Buffer(0.001);
                         ISpatialFilter spatialFilter0 = Utils.FilterUtils.SpatialFilter(pGeometryPointBuffer, esriSpatialRelEnum.esriSpatialRelIntersects);
                         featureCursor0 = featureClass.Update(spatialFilter0, false);
                         IFeature adjoinFeature = featureCursor0.NextFeature();
                         while (adjoinFeature != null)
                         {
                             IPointCollection pointCollectionAdjoinFeature = adjoinFeature.Shape as IPointCollection;
                             IPointCollection pointCollectionAdjoinFeatureCpoy = adjoinFeature.ShapeCopy as IPointCollection;
                             int count = ;
                             for (int k = ; k < pointCollectionAdjoinFeature.PointCount - ; k++)
                             {
                                 if (Math.Abs(pointCollectionAdjoinFeature.Point[k].X - pointCollectionNewFeature.Point[i].X) < 0.001 && Math.Abs(pointCollectionAdjoinFeature.Point[k].Y - pointCollectionNewFeature.Point[i].Y) < 0.001)
                                 {
                                     count++;
                                     break;
                                 }
                             }
                             if (count == )
                             {
                                 for (int j = ; j < pointCollectionAdjoinFeatureCpoy.PointCount - ; j++)
                                 {
                                     if(IsPointOnSegment(pointCollectionNewFeature.Point[i], pointCollectionAdjoinFeatureCpoy.Point[j], pointCollectionAdjoinFeatureCpoy.Point[j + ]))
                                     {
                                         //before 在这个索引之前加
                                         pointCollectionAdjoinFeature.AddPoint(pointCollectionNewFeature.Point[i], j + , Type.Missing);
                                         break;
                                     }
                                 }
                                 ITopologicalOperator topologicalOperator = adjoinFeature.Shape as ITopologicalOperator;
                                 topologicalOperator.Simplify();
                                 featureCursor0.UpdateFeature(adjoinFeature);
                                 featureCursor0.Flush();
                             }
                             adjoinFeature = featureCursor0.NextFeature();
                         }
                     }
                     newFeature = newFeatureCursor.NextFeature();
                 }
             }
             catch (Exception exp)
             {
                 throw new Exception("添加公共节点时发生错误！" + "\r\n" + exp.Message + "\r\n" + exp.StackTrace);
             }
         }