基于ArcEngine和VSTO的勘测定界地类自动化分析方法

朱高峰 ZHU Gao-feng

（广州市城市规划勘测设计研究院，广州 510000）

0 引言

勘测定界是一项为满足土地征收、划拨、农用地转用、违法用地处罚及土地规划等土地管理工作需要而进行的实地界定土地使用范围、调绘土地利用现状、计算用地面积的基础性工作[1]。由于土地勘测定界成果为土地管理提供了基础的地类面积数据，其数据准确性会对相应的土地管理业务正确性开展产生直接影响。

目前，很多从事土地管理工作的技术人员和研究人员都对土地勘测定界进行了研究，相关的理论及技术方法较为成熟[2，3，4]。同时，土地管理也是一个动态变化的过程，不同时期土地管理的目标和方法均会存在差异。本文在广州市启用“三调”数据作为土地管理基础数据的背景下，根据新时期土地管理工作对勘测定界的要求，讨论了土地勘测定界地类自动化分析实现方法和关键技术。

1 广州市勘测定界地类分析现状

所谓勘测定界地类分析，主要是指依托最新年份和历史年份的土地利用现状数据，对目标地块的土地利用现状地类及历史地类变化情况进行分析，最终将分析得到的地类面积数据以表格的形式呈现。目前，测绘人员主要采取ArcGIS和Excel相结合的方式进行地类面积变化分析与统计，这种地类分析方法不仅工作量大，而且容易出错[5]。

随着《土地管理法》的修订，广东省以《广东省自然资源厅关于建设用地报批地类有关问题的通知》停用了原有按照最新年份、2009年、2007年、1999年四个年份的用地审查地类分析规则，新的地类分析规则需要由最新年份逐年追溯至1999年，并新增了耕地占补平衡分析的要求。2022年4月，全国范围启用“三调”数据，用地审查所依据的土地利用现状分类需要在原有2007版、2002版和1984版的基础上新增“三调”版本。在新的地类分析规则以及启用“三调”数据的双重背景下，现有采用人工统计的方式进行土地勘测定界地类分析工作变得繁琐且易出错。

2 勘测定界地类自动化分析关键技术实现

2.1 主要技术流程

勘测定界地类分析主要流程如图1所示，主要包括空间数据分析、地类面积统计、表格格式调整三大部分。

图1 土地勘测定界作业流程图及方法对比

2.2 基于ArcEngine的空间数据自动分析

ArcEngine中封装了大量的ArcObjects接口，通过ArcEngine二次开发可以轻松实现ArcMap中的相关功能的配置。利用ArcEngine中的GP工具可以实现勘测定界中的空间数据叠加分析、地块面积计算、要素属性赋值等功能[6]。通过循环语句重复调用GP工具，即可完成各个年份土地利用现状数据的叠加分析工作，其核心代码如图2所示。

图2 空间分析功能核心代码

2.3 基于VSTO的数据统计及格式整饰

所谓VSTO，就是在Visual Studio的开发环境中，使用C#语言开发用于微软Office的插件或文档。由于ArcEngine二次开发是在Visual Studio环境下进行的，使用VSTO可以很好的将ArcEngine的空间数据分析功能与VBA的数据统计功能结合起来。在Excel中使用VBA定义实现地类面积数据统计及格式整饰的代码，再通过VSTO实现对VBA代码的调用，从而实现勘测定界地类自动化分析的目的[7]。VSTO调用VBA核心代码如图3所示。

图3 VSTO核心代码

3 应用实例与结果

软件基于Esri开发组件ArcEngine和VSTO，采用C#和VBA语言开发完成，在空间分析上主要实现了现状建设用地逐年倒推、可调整地类逐年倒推、设施农用地逐年倒推和非耕农用地逐年倒推功能，在数据导出上主要实现了耕地数据、违法用地数据和地类统计表自动导出功能。程序界面如图4所示。

图4 软件主界面

在程序界面中输入必要的数据信息，在数据分析模块和成果数据导出模块中点击相应的功能按钮，即可实现相应的地类分析和数据导出功能。图5为通过建设用地倒推功能导出的地类叠加分析结果，根据该部分结果，结合ArcMap制图功能可以方便的进行各个年份土地利用现状图绘制。表1为通过生成地类表功能生成的勘测定界土地分类面积汇总表标准表格（部分）。

图5 地类叠加分析结果（部分）

表1 地类面积表（部分） 单位：公顷

4 结语

本文根据广州市土地勘测定界现状基础数据及地类分析的实际需求，设计了一种基于ArcEngine空间数据分析功能和VSTO数据统计功能，采用C#、VBA语言实现的地类自动分析统计方法，并介绍了自动分析过程中的关键技术。实践证明，根据该方法统计的地类分析数据准确，符合《土地勘测定界规范》要求，能够满足实际生产需要。