基于ArcGIS的VBA多边形空洞去除方法研究

2013-08-29张金兰欧阳婷萍

大众科技 2013年3期

张金兰欧阳婷萍

(1.广东工贸职业技术学院，广东广州 510510；2.中国科学院广州地球化学研究所边缘海地质重点实验室，广东广州 510640)

在GIS中，矢量数据模型是地理数据的最主要表现形式。在矢量数据模型中，将二维地理实体抽象为点、线、面（多边形）三类要素对象[1]。这些几何形体对象可以精确描述具有离散特征的要素在地球上的具体位置和形状，并能够表达出地理实体之间的空间关系。正是因为空间关系，才使得GIS可以对这些要素进行不同的空间分析和运算，得出用户需要的结果，这也是GIS区别于一般地理信息系统的本质——空间查询和分析功能。空间关系主要包括方向关系、拓扑关系和度量关系。空间拓扑关系是空间实体间的一个重要而又基本的关系，是GIS空间推理、查询和分析的基础[2-7]。多边形的拓扑关系可归纳为“相离”、“相接”、“相交”、“覆盖”、“包含”、“覆盖于”、“包含于”、“相等”关系[6,7]。但对于这些关系的描述和应用都是在简单无洞多边形的基础上的。如果多边形存在空洞（内环），便会对拓扑关系造成困扰，干扰GIS的应用。因此经常需要对多边形的空洞（内环）进行删除。ESRI公司的ArcGIS系列软件是一个全面的、完善的、功能强大的GIS软件平台[8]。ArcGIS虽然能够对空洞多边形进行去除，但只能对空洞（内环）一个一个的进行删除，不能达到批量处理的效果。本文针对上述问题，利用ArcGIS中内置VBA对多边形空洞进行批量的删除处理，生成无洞多边形。经试验表明，该方法简单，效果理想，为GIS的数据处理和应用分析提高效率，更适用于广大GIS用户。

1 VBA介绍

众所周知，Arcmap中包含了内置的VBA，VBA全称Visual Basic For Application，是 Visual Basic的子集，作为ArcGIS的二次开发工具。VBA提供了一种集成的编程环境——VB编辑器（VBE）。在该环境中，可以编写并调试在Arcmap中使用的VB宏及ESRI的对象库，宏的命令可以将VB的一些或部分功能与Arcmap的功能集成起来。创建宏的过程就是编写VB的子程序，改程序的名字就是宏的名字[9]。VBA的开发使用ActiveX技术，这种技术能使用户能够从ArcGIS的内部或外部以编程方式来操作ArcGIS，同时具有一个功能强大的对象库——ArcObjects。ArcGIS与VBA的结合具有如下特点：

（1）VBA可以通过应用程序满足用户的需要，增强与用户的交互，集成ArcGIS的功能，使ArcGIS任务自动化，增强ArcGIS的功能[10,11]。

（2）VBA是VB的子集，二者结构和语法完全一致，易于掌握和快速应用，并且内嵌于ArcGIS软件中，无需安装[11,12]。

2 基于ArcGIS的VBA多边形空洞去除的实现

2.1 算法依据

在ArcObjects中， Polygon（多边形）是一个有序Ring（环）对象的集合，环有内环和外环之分，内环及多边形的空洞。内环和外环的存储方向是有区别的，内环存储方向是逆时针方向，而外环是顺时针方向，这也决定了环的面积特性有所区别。内环的面积为负值，外环的面积为正值[13-15]。这是判断内环和外环的关键，也是去除多边形空洞的算法依据。

2.2 算法实现

（1）在ArcGIS的“ArcMap—T0ols—Macros-Create”下，输入文件名，编好程序后，若想保存文件，点VBE中的“File—save Project”图标，输入文件名，则文件以后缀名为．mxd的格式保存．当下次需要运行时，打开“ArcMap—tools—Macms—VBE—Run”，运行宏的结果将在桌面显示。

（2）在本次开发中主要用到Polygon类和Ring类，首先通过ITopologicalOperator2接口的SimplifyPreserveFromTo函数确保多边形的内环和外环的方向是正确的；接着通过IArea接口的Area属性值是否小于0，来筛选内环进行删除，只保留多边形的外环，完成多边形空洞的去除。部分代码如下：