在做數據處理時,我們經常遇到的一個問題就是分類問題,如下圖所示,現有一份中國各省、自治區、直轄市的shp文件,現在要求根據Name字段將該shp文件分割成多個shp文件,每個shp文件代表一個省、自治區或直轄市。
熟悉ArcGIS的同志應該清楚,在ArcToolbox工具箱下有一個分割工具可以幫我們實現該功能,如下圖所示:
其實在ArcEngine中我們也可以實現這一功能,核心接口爲IFeatureDataConverter
,代碼如下:
using System;
using System.Collections;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
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();
}
// 分割
private void btnSplit_Click(object sender, EventArgs e)
{
IFeatureClass pFeatureClass = GetFeatureClass(@"E:\Users\dsf\Desktop\data\china.shp");
List<string> list = GetUniqueValue(pFeatureClass, "Name");
foreach (string item in list)
{
IQueryFilter pQueryFilter = new QueryFilter();
pQueryFilter.WhereClause = "Name='" + item + "'";
SplitByAttribute(pFeatureClass, pQueryFilter, @"E:\Users\dsf\Desktop\新建文件夾\" + item + ".shp");
}
}
// 獲取要素類
private IFeatureClass GetFeatureClass(string filePath)
{
IWorkspaceFactory pWorkspaceFactory = new ShapefileWorkspaceFactory();
IWorkspaceFactoryLockControl pWorkspaceFactoryLockControl = pWorkspaceFactory as IWorkspaceFactoryLockControl;
if (pWorkspaceFactoryLockControl.SchemaLockingEnabled)
{
pWorkspaceFactoryLockControl.DisableSchemaLocking();
}
IWorkspace pWorkspace = pWorkspaceFactory.OpenFromFile(System.IO.Path.GetDirectoryName(filePath), 0);
IFeatureWorkspace pFeatureWorkspace = pWorkspace as IFeatureWorkspace;
IFeatureClass pFeatureClass = pFeatureWorkspace.OpenFeatureClass(System.IO.Path.GetFileName(filePath));
return pFeatureClass;
}
// 獲取唯一值
private List<string> GetUniqueValue(IFeatureClass pFeatureClass, string fieldName)
{
IQueryFilter pQueryFilter = new QueryFilter();
pQueryFilter.AddField(fieldName);
// 要素遊標
IFeatureCursor pFeatureCursor = pFeatureClass.Search(pQueryFilter, true);
ICursor pCursor = pFeatureCursor as ICursor;
// 設置統計信息
IDataStatistics pDataStatistics = new DataStatistics();
pDataStatistics.Field = fieldName;
pDataStatistics.Cursor = pCursor;
// 獲取唯一值
IEnumerator uniqueValues = pDataStatistics.UniqueValues;
uniqueValues.Reset();
// 遍歷唯一值
List<string> list = new List<string>();
while (uniqueValues.MoveNext())
{
list.Add(uniqueValues.Current.ToString());
}
return list;
}
// 根據屬性分割
private void SplitByAttribute(IFeatureClass pSourceFeatureClass, IQueryFilter pQueryFilter, string filePath)
{
// 源數據工作空間
IDataset pSourceDataset = pSourceFeatureClass as IDataset;
IFeatureClassName pSourceFeatureClassName = pSourceDataset.FullName as IFeatureClassName;
IWorkspace pSourceWorkspace = pSourceDataset.Workspace;
// 目標數據工作空間
IWorkspaceFactory pTargetWorkspaceFactory = new ShapefileWorkspaceFactory();
IWorkspace pTargetWorkspace = pTargetWorkspaceFactory.OpenFromFile(System.IO.Path.GetDirectoryName(filePath), 0);
IDataset pTargetDataSet = pTargetWorkspace as IDataset;
// 設置目標數據屬性
IWorkspaceName pTargetWorkspaceName = pTargetDataSet.FullName as IWorkspaceName;
IFeatureClassName pTargetFeatureClassName = new FeatureClassName() as IFeatureClassName;
IDatasetName pTargetDatasetName = pTargetFeatureClassName as IDatasetName;
pTargetDatasetName.WorkspaceName = pTargetWorkspaceName;
pTargetDatasetName.Name = System.IO.Path.GetFileName(filePath);
// 檢查字段
IFieldChecker pFieldChecker = new FieldChecker();
pFieldChecker.InputWorkspace = pSourceWorkspace;
pFieldChecker.ValidateWorkspace = pTargetWorkspace;
// 字段轉換
IFields pSourceFields = pSourceFeatureClass.Fields;
IFields pTargetFields = null;
IEnumFieldError pEnumFieldError = null;
pFieldChecker.Validate(pSourceFields, out pEnumFieldError, out pTargetFields);
// 數據轉換
IFeatureDataConverter pFeatureDataConverter = new FeatureDataConverter();
pFeatureDataConverter.ConvertFeatureClass(pSourceFeatureClassName, pQueryFilter, null, pTargetFeatureClassName, null, pTargetFields, "", 1000, 0);
}
}
}
最終的分割結果如下圖所示: