ArcEngine中的叠置分析

对于叠置分析,想必用过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;
        }
    }
}

运行结果如下图所示,黄色区域就是两个图层的相交区域。
在这里插入图片描述

  • 5
    点赞
  • 30
    收藏
    觉得还不错? 一键收藏
  • 9
    评论
评论 9
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值