ArcEngine實現多邊形內部的挖洞操作

在一般的數據編輯流程中，有些時候可能會涉及到多邊形內部的挖洞操作。在ArcEngine中，一般可以通過兩種方式實現挖洞操作。一種是疊置分析，這主要是針對圖層與圖層之間的操作，我會在後續的博客中進行介紹，另一種則是用戶自己繪製多邊形實現挖洞操作，這主要是利用ITopologicalOperator接口進行實現。先來看一下效果圖：

代碼如下：

using System;
using System.Collections;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Diagnostics;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
using System.Windows.Forms;
using ESRI.ArcGIS.Carto;
using ESRI.ArcGIS.Controls;
using ESRI.ArcGIS.DataSourcesFile;
using ESRI.ArcGIS.DataSourcesGDB;
using ESRI.ArcGIS.Display;
using ESRI.ArcGIS.esriSystem;
using ESRI.ArcGIS.Geodatabase;
using ESRI.ArcGIS.Geometry;
using ESRI.ArcGIS.Output;
using ESRI.ArcGIS.SystemUI;

namespace WindowsFormsApplication1
{
    public partial class MainForm : Form
    {
        public MainForm()
        {
            InitializeComponent();
            axMapControl1.LoadMxFile(@"E:\Users\dsf\Desktop\無標題.mxd");
        }

        // OnMouseDown事件
        private void axMapControl1_OnMouseDown(object sender, IMapControlEvents2_OnMouseDownEvent e)
        {
            IFeatureClass pFeatureClass = (axMapControl1.get_Layer(0) as IFeatureLayer).FeatureClass;
            IGeometry pGeometry = axMapControl1.TrackPolygon();
           
            // 空間過濾器
            ISpatialFilter pSpatialFilter = new SpatialFilter();
            pSpatialFilter.AddField("Shape");
            pSpatialFilter.Geometry = pGeometry;
            pSpatialFilter.SpatialRel = esriSpatialRelEnum.esriSpatialRelIntersects;

            // 要素遊標
            IFeatureCursor pFeatureCursor = pFeatureClass.Search(pSpatialFilter, true);
            IFeature pFeature = pFeatureCursor.NextFeature();
            if (pFeature == null)
            {
                return;
            }

            // 遍歷遊標
            while (pFeature != null)
            {
                ITopologicalOperator pTopologicalOperator = pFeature.ShapeCopy as ITopologicalOperator;
                pFeature.Shape = pTopologicalOperator.Difference(pGeometry);
                pFeature.Store();
                pFeature = pFeatureCursor.NextFeature();
            }

            // 刷新地圖
            System.Runtime.InteropServices.Marshal.ReleaseComObject(pFeatureCursor);
            axMapControl1.ActiveView.Refresh();
        }
    }
}