基于数据分治与双层索引的并行点面叠加分析方法研究.pdfVIP

第31卷 第 2期 地理与地理信 息科 学 Vo1．31 No．2 2015年 3月 GeographyandGeo-lnformationScience March2015 doi：10．3969／j_issn．1672—0504．2015．02．001 基于数据分治与双层索引的并行点面叠加分析方法研究 周玉科 ，周成虎 ，马 廷 ，高锡章 ，范俊甫 ，许 涛 ，季 民。 (1．中国科学院地理科学与资源研究所，资源与环境信息系统国家重点实验室，北京 100101； 2．山东理工大学建筑工程学院，山东 淄博 255049；3．山东科技大学测绘工程学院，山东 青岛 266510) 摘要：地图叠加分析是一种计算密集型算法，并行化计算是加快算法执行速度的一种有效方法。该文研究分布式 环境下的点面图层并行化叠加分析方法与实现 。首先根据点面叠加的特点设置并行数据分解的方式 ，基于分治法 分解空间数据 ，在并行系统下将地理要素分而治之。然后引入双层索引的并行叠加机制，一是对面图层根据 Hil— bert空间索引的排序方式分发数据 ，二是对点图层建立四叉树索引，对每一个进行相交运算的多边形进行快速过 滤和求交。最后在 Linux集群系统下实现该并行算法，其一利用MP1分布式计算环境实现在整体计算框架下的消 息通讯模式的并行，其二在每个子节点中实现基于多核OpenMP工具的本地并行化。结果表明，利用双层空间索 引分治的方法可实现并行数据分块 ，各子节点实现独立计算 ，减少并行系统中的I／O冲突，并行加速比明显。该方 法对矢量地图运算的并行化进行了有益的尝试，为大数据时代的空间数据分析提供一种有效的途径。 关键词 ：地图叠加分析；并行计算；空间索引；MPI；OpenMP 中图分类号：P208 文献标识码：A 文章编号：1672--0504(2015)02--0001--06 果判断点被多边形包含关系_l5]。已有学者利用各 0 引言 种并行计算工具对地图并行叠加分析进行探索，如 地图叠加分析(MapOverlay)是 GIS空间分析 利用GPU加速_7]和集群并行方法F89]。本文利用并 中最基础、使用最频繁的一种操作_l1]，指同一地域范 行计算的基础理论方法，重点研究在并行化地图叠 围的两个或多个地图图层在同样空问参考系下进行 加分析中的静态负载均衡方法，以数据并行为主线， 叠加，获取具有新属性的空间区域并最终生成一个 利用 Hilbert空间索引分解地理数据，达到各任务节 新的图层，图层要素属性由叠加运算符决定，其本质 点计算量基本平衡，各个节点再次进行细粒度并行 过程是一系列计算几何布尔操作和属性传递过程的 判断包含关系。 集合。图层叠加中几何对象之间的操作已被证明为 1 双层索引方法与分治并行策略 时间复杂度最少为 O(nlogn)的操作 ]，因此涉及全 幅地图图层的叠加分析属于计算密集型操作。点与 1．1 并行叠加特征预分析 多边形图层的叠加分析是其 中基本的操作类型之 本文针对点图层与面图层的叠加分析，运算对 一 ，其实质是点包含问题(单点或复杂对象组成点)。 象之间属于多对多的关系，因此研究对象的粒度应 随着数据规模的日益膨胀、实时性应用的增加， 该从点与边的关系上升到点集合与多边形集合之间 对地图叠加分析功能在计算效率、性能、处理能力等 的关系、点集 内部 的空间分布特征关系和多边形集 方面的要求也越来越高。并行计算技术从 20世纪 合之间的空间近邻关系。传统的并行算法设计考虑 7O年代随着计算机体系架构的进步和成熟发展起 每个原子级别的操作并保证每个节点处理时间大致 来，由于并行技术发展时间不长，因此 GIS的核心空 相同。如果按照传统设计方法将点和多边形相交降 间分析功能的并行化仍然处在迅速发展时期F3,4]。 低维度将增加数据组织和计算的复杂度。假设点和 近年来，并行计算机 的强大硬件加速特性和新的并 多边形数量均为 ，采用直接求交方法的时间