ArcEngine中的叠置分析

最新推荐文章于 2021-03-01 17:55:41 发布

HerryDong

最新推荐文章于 2021-03-01 17:55:41 发布

阅读量3.5k

点赞数 5

分类专栏： ArcEngine 文章标签： ArcEngine

本文链接：https://blog.csdn.net/HerryDong/article/details/105620620

版权

ArcEngine 专栏收录该内容

59 篇文章 74 订阅

订阅专栏

对于叠置分析，想必用过ArcGIS的同志应该都很了解。ArcGIS的工具箱提供了各种叠置分析的工具，如下图所示：
在这里插入图片描述
那么在ArcEngine中该如何实现叠置分析呢？其实ArcEngine已经给我们提供了IBasicGeoprocessor接口来实现叠置分析，下面就以求两个图层的相交部分来进行说明。

测试数据

在这里插入图片描述

叠置分析

代码如下：

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");
        }

        private void btn_Click(object sender, EventArgs e)
        {
            IFeatureClass pSourceFeatureClass = (axMapControl1.get_Layer(0) as IFeatureLayer).FeatureClass;
            IFeatureClass pTargetFeatureClass = (axMapControl1.get_Layer(1) as IFeatureLayer).FeatureClass;
            IFeatureClass pFeatureClass = Intersect(pSourceFeatureClass, pTargetFeatureClass, @"E:\Users\dsf\Desktop\test\相交.shp");

            IFeatureLayer pFeatureLayer = new FeatureLayer();
            pFeatureLayer.FeatureClass = pFeatureClass;
            pFeatureLayer.Name = pFeatureClass.AliasName;

            axMapControl1.AddLayer(pFeatureLayer, 0);
            axMapControl1.Refresh();
        }

        private IFeatureClass Intersect(IFeatureClass pSourceFeatureClass, IFeatureClass pTargetFeatureClass, string filePath)
        {
            ITable pSourceTable = pSourceFeatureClass as ITable;
            ITable pTargetTable = pTargetFeatureClass as ITable;

            // IWorkspaceName
            IWorkspaceName pWorkspaceName = new WorkspaceName() as IWorkspaceName;
            pWorkspaceName.WorkspaceFactoryProgID = "esriDataSourcesGDB.ShapefileWorkspaceFactory";
            pWorkspaceName.PathName = System.IO.Path.GetDirectoryName(filePath);

            // IFeatureClassName
            IFeatureClassName pFeatureClassName = new FeatureClassName() as IFeatureClassName;
            pFeatureClassName.ShapeFieldName = "Shape";
            pFeatureClassName.ShapeType = esriGeometryType.esriGeometryPolygon;
            pFeatureClassName.FeatureType = esriFeatureType.esriFTSimple;

            // IDatasetName
            IDatasetName pDatasetName = pFeatureClassName as IDatasetName;
            pDatasetName.WorkspaceName = pWorkspaceName;
            pDatasetName.Name = System.IO.Path.GetFileNameWithoutExtension(filePath);

            // 计算相交部分
            IBasicGeoprocessor pBasicGeoprocessor = new BasicGeoprocessor();
            IFeatureClass pFeatureClass = pBasicGeoprocessor.Intersect(pSourceTable, false, pTargetTable, false, 0.001, pFeatureClassName);
            return pFeatureClass;
        }
    }
}