泰森多边形(Voronoi图)生成算法.docVIP

泰森多边形(Voronoi图)生成算法 作者 谢卫国 创建时间 2008-10-27 审核人 审核时间 2008-10-27 文档状态 草稿 文档版本 Ver 1.0 审批人 审批时间 最后修改人 谢卫国 最后修改时间 2008-10-27 文档编号 面向人员 技术开发人员 文档目的 本文描述了在geomodel模块中，生成泰森多边形所使用的算法。 概述 GIS和地理分析中经常采用泰森多边形进行快速插值，和分析地理实体的影响区域，是解决邻接度问题的又一常用工具。荷兰气候学家A·H·Thiessen提出了一种根据离散分布的气象站的降雨量来计算平均降雨量的方法，即将所有相邻气象站连成三角形，作这些三角形各边的垂直平分线，于是每个气象站周围的若干垂直平分线便围成一个多边形。用这个多边形内所包含的一个唯一气象站的降雨强度来表示这个多边形区域内的降雨强度，并称这个多边形为泰森多边形。如图1，其中虚线构成的多边形就是泰森多边形。泰森多边形每个顶点是每个三角形的外接圆圆心。泰森多边形也称为Voronoi图，或dirichlet图。 图1 泰森多边形泰森多边形的特性是：每个泰森多边形内仅含有一个离散点数据泰森多边形内的点到相应离散点的距离最近位于泰森多边形边上的点到其两边的离散点的距离相等泰森多边形可用于定性分析、统计分析、邻近分析等。例如，可以用离散点的性质来描述泰森多边形区域的性质；可用离散点的数据来计算泰森多边形区域的数据；判断一个离散点与其它哪些离散点相邻时，可根据泰森多边形直接得出，且若泰森多边形是n边形，则就与n个离散点相邻；当某一数据点落入某一泰森多边形中时，它与相应的离散点最邻近，无需计算距离。在泰森多边形的构建中，首先要将离散点构成三角网。这种三角网称为Delaunay三角网。Delaulay三角形的构建Delaunay三角网的构建也称为不规则三角网的构建，就是由离散数据点构建三角网，如图2，即确定哪三个数据点构成一个三角形，也称为自动联接三角网。即对于平面上n个离散点，其平面坐标为(xi，yi)，i＝1，2，…，n，将其中相近的三点构成最佳三角形，使每个离散点都成为三角形的顶点。?? 图2 Delaunay三角网 自动联接三角网的结果为所有三角形的三个顶点的标号，如1,2,8；2,8,3；3,8,7；…… 为了获得最佳三角形，在构三角网时，应尽可能使三角形的三内角均成锐角，即符合Delaunay三角形产生的准则：任何一个Delaunay三角形的外接圆内不能包含任何其它离散点。相邻两个Delaunay三角形构成凸四边形，在交换凸四边形的对角线之后，六个内角的最小者不再增大。该性质即为最小角最大准则。 图3 凸包下面介绍Tsai（1993）提出的在n维欧拉空间中构造Delaunay三角形的通用算法---凸包插值算法。（一）、凸包生成1、求出点集中满足min(x-y)、min(x+y)、max(x-y)、max(x+y)的四个点，并按逆时针方向组成一个点的链表。这4个点是离散点中与包含离散点的外接矩形的4个角点最近的点。这4个点构成的多边形作为初始凸包。2、对于每个凸包上的点I，设它的后续点为J，计算矢量线段IJ右侧的所有点到IJ的距离，求出距离最大的点K。3、将K插入I、J之间，并将K赋给J。4、重复2、3步，直到点集中没有在线段IJ右侧的点为止。5、将J赋给I，J取其后续点，重复2、3、4步。6、当凸包中任意相邻两点连线的右侧不存在离散点时，结束点集凸包求取过程。完成这一步后，形成了包含所有离散点的多边形（凸包），如图3所示。（二）、环切边界法凸包三角剖分在凸包链表中每次寻找一个由相邻两条凸包边组成的三角形，在该三角形的内部和边界上都不包含凸包上的任何其它点。将这个点去掉后得到新的凸包链表。重复这个过程，直到凸包链表中只剩三个离散点为止。将凸包链表中的最后三个离散点构成一个三角形，结束凸包三角剖分过程。 图4 环切边界法凸包三角剖分完成这一步后，将凸包中的点构成了若干Delaunay三角形，如图4所示。（三）、离散点内插在对凸包进行三角剖分之后，不在凸包上的其余离散点，可采用逐点内插的方法进行剖分。基本过程为：1、2、3、 4、重复1、2、3，直到所有非凸壳离散点都插入完为止。完成这一步后，就完成了Delaunay三角网的构建，如图5所示。 图5 离散点内插泰森多边形的建立步骤建立泰森多边形算法的关键是对离散数据点合理地连成三角网，即构建Delaunay三角网。建立泰森多边形的步骤为：1、离散点自动构建三角网，即构建Delaunay三角网。对离散点和形成的三角形编号，记录每个三角形是由哪三个离散点构成的。2、找出与每个离散点相邻的所有三角形的编号，并记录下来。这只要在已构建的三角网中找出具有一个相