种基于栅格方法的空间谓词判断方法及其系统.docVIP

一种基于栅格方法的空间谓词判断方法及其系统 董 慧* ** ，程振林*，方金云*， 赵红超* （* 中国科学院计算技术研究所 北京市海淀区科学院南路6号 100190） （** 中国科学院研究生院 北京市石景山区玉泉路19号（甲） 100049） 摘要：地理信息获得了越来越广泛与深入的应用。空间分析是最的计算之一。矢量方法是面向物体的描述，物体间的几何关系4CRS（four color raster signature）空间是GIS的核心，是GIS区别于一般的信息系统、CAD或者电子地图系统的主要标志之一。Open Geospatial Consortium，开放地理信息联盟Disjoint/ Intersect、Equals、Within/Contains、Overlaps、BBOX等多种空间谓词过滤器。 空间，配合空间数据的属性信息，能提供强大、丰富的空间数据查询功能。随着网络地图服务的流行，如何在网络地图服务器上提供空间功能成为。作为计算机科学的一个分支，计算几何主要研究解决几何问题的算法。针对两个多边形进行，在的多边形计算面前无能为力。采用“暴力”算法，通过反复调用两个多边形的算法来完成算法实现复杂度。如何确定交点的进点、出点属性在实际的图形中会遇到众多的特殊情况。特别是在发生了线段与线段交在端点、线段与线段重叠的情况下，如何区分交点的出点、入点情况非常复杂。这类中采用的求交算法优点是结果比较精确，缺点是要进行频繁的坐标排序、角度计算等操作。four color raster signature --4CRS)[3][4]可以实现空间关系的判断，但是4CRS结果不够精确，只能以四种栅格属性区分栅格面积空、满、一半以上、一半以下），只能判断出强*强是确定相交的，而强*弱、弱*强、弱*弱则都是不确定的结果，还需要精确计算。在计算近似面积、置信度区间等都是用数学期望和概率等公式来估计，而真实数据可能并不符合这种规律。本文使用的基于亚像素精度的方法能准确记录边界格子的面积，因此可以根据两个对象占用的单元格的面积来确定是否相交。如同一个位置的栅格，第一个图层占有此栅格为49%，第二个图层为52%，相加100%,所以确定相交，而4CRS则不能确定。而且，4CRS只能实现两个多边形的判断，而不能处理多个多边形的图层。而本文采用的方法不但能精确记录栅格面积百分比，提高结果的准确性，并能把百分比和图层属性信息等都保存下来，给用户返回更多有效的结果信息（结果多边形的ID号等）。整个栅格索引过程算法框架如图1所示。首先，对地理要素的矢量点进行坐标转换，便于利用亚像素精度进行的计算。绘制控制器计算填充单元跨段 图1. 步骤一的算法框架 经过底图索引生成这个步骤之后，即完成了栅格化过程，以栅格矩阵形式存储，其中坐标值为索引，所属的多边形ID号和所占的栅格面积百分比等为图层信息，以文件形式保存下来。 这样，两个图层的空间谓词判断的时候，只需要找到同一栅格，对其面积等进行判断，从而获得符合条件的多边形ID号返回结果。示意图如图2。 图2 步骤二的流程图 三．基于栅格方法的空间谓词实现 栅格索引生成 栅格索引生成部分，用于对于输入的图层生成对应的栅格底图，栅格底图的栅格单元以压盖所述栅格单元的图层的图形的要素ID为要素索引，栅格单元以坐标值为索引 图3 (a) 地理要素的几何表示 (b) 栅格近似 实现图片像素与矢量点关联，具体指在绘制的同时，生成一张与匹配的栅格底图，该底图就是对应的索引图每个都对应了内容，因此完全不存在搜索存在重叠区域这种情况。地理要素的矢量点进行坐标转换坐标值为索引将double类型的矢量点坐标都乘以256相当于将该坐标的二进制表示左移8位。这种坐标转换的优点是考虑了小数部分对像素的栅格单元权值(cover)的影响便于利用亚像素精度进行的计算。 轮廓线的每一权值和覆盖面积 （1） （2） 其中，（fx1，fy1）为经过该栅格单元的轮廓的线段起始点的小数坐标部分，（fx2，fy2）为经过该栅格单元的轮廓的线段终止点的小数坐标部分，cover为权值，area为面积。 运用了亚像素精度（subpixel accuracy）的布兰森汉姆Bresenham）[6]生成直线算法，Bresenham生成直线算法是一种基于误差判别式来生成直线的方法。与传统Bresenham所不同的是，该算法利用误差判别选择像素的过程是基于亚像素的，把一个像素分成NN个小像素，将单位栅格平均分成了256256个子像素。绘制控