《快速凸包算法ppt课件.ppt》由会员分享，可在线阅读，更多相关《快速凸包算法ppt课件.ppt（24页珍藏版）》请在三一办公上搜索。
1、Advanced Graphics,孙 晓 鹏博士 教授2011年 11月 16日,第二章 二维凸包,2.4 凸包的快速算法,主要思想点集S 的凸包是取决于凸包边界附近的点逐步丢掉凸包内部的点，只关注凸包附近的点，从而提高算法的效率最好情况O(nlogn)、 最坏情况O(n2),2.4 凸包的快速算法,算法过程取两个极端点，它们是最右最下点pdr和最左最上点pul有向直线 pdr pul将整个凸包被划分为右凸包和左凸包对右凸包和左凸包分别进行递归递归设S1是严格在直线 pdr pul右边的点集（S1可能是空集）在S1中寻找距离直线 pdr pul最远的点，作为pdr pul右边的一个极端点b连
2、接pdr和b ，及b 和pul 把pdr右侧的点集记为A， pul右侧的点集的点记为B对边pdr b和点集A、对边b pul 和点集B分别递归调用依次连接凸包上的顶点，得点集S1的凸包，即点集S的右凸包类似地，计算出点集S的左凸包，从而得到整个点集S的凸包,2.4 凸包的快速算法,算法过程取两个极端点，它们是最右最下点pdr和最左最上点pul有向直线 pdr pul将整个凸包被划分为右凸包和左凸包对右凸包和左凸包分别进行递归,2.4 凸包的快速算法,算法过程递归设S1是严格在直线 pdr pul右边的点集（S1可能是空集）在S1中寻找距离直线 pdr pul最远的点，作为 pdr pul右边的
3、一个极端点b连接 pdr和b ，及b 和 pul 把 pdrb右侧的点集记为A， bpul右侧的点集的点记为B对边 pdrb和点集A、对边bpul 和点集B分别递归调用,2.4 凸包的快速算法,最好情况出现在每次划分均是平衡的, O(nlogn)最坏情况出现在每次划分点的分布都很极端, O(n2),2.5 Graham算法,20世纪60年代末贝尔实验室需要求解10,000个点的凸包O(n2)的方法太慢1972年Graham出O(nlogn)的二维凸包算法,2.5 Graham算法,基本思想在凸包内部找到一个点o 如S 中任何三个不共线的点的重心，O(1)将o作为极坐标的中心，计算每个点的极角对
4、S中的点按升序排列（如pi ，pi+1 ， pi+2），O(nlogn)计算相邻三点转角的凹凸性，删除内凹的点O(n)当点集内不再包含内凹的点时，得到凸包,2.5 Graham算法,以极端点 pi为初始点，依次对相邻三个点pi ，pi+1和pi+2 ，计算pi pi+1pi+1pi+2如果在z 轴上的投影大于零，即(pi pi+1pi+1pi+2)z0说明在pi+1 处左转弯，多边形在该点上外凸，暂时保留这三点前进一步，同样去判断相邻三个点pi+1，pi+2和 pi+3如果(pi pi+1pi+1pi+2)z 0说明在pi+1处右转弯，多边形在该点上内凹，把pi+1点从多边形边界中删除后退一步
5、，同样去判断相邻三个点pi-1，pi和 pi+2时间复杂度为线性O(n),凸包计算方法对比,极端边算法O(n3)礼品包裹算法O(n2)快速算法最好情况O(nlogn)、 最坏情况O(n2)Graham 算法排序计算O(nlogn)、执行时间O(n)总的时间复杂度O(nlogn),第三章 凸包扩展,3.1 多面体,两个集合是同胚的(homeomorphic)指它们之间存在一个连续的一一映射并且这个映射的逆映射也是连续的两个同胚的集合允许它们各自拉伸和扭曲，但只要不撕裂，其结果仍然同胚如果任一集合被撕裂了，映射的连续性便被破坏，两个集合就不再同胚了,3.1 多面体,两个集合是同胚的(homeomo
6、rphic)三黑点的邻域都不能同胚于一个三维的半开半闭的半球a和b中黑点的邻域同胚于两个接触于一条线或一个点的三维的半开半闭的半球c中黑点的邻域是同胚于半个二维开圆盘,3.1 多面体,3.1.1 多面体定义三维空间中的多面体是由有限个平面多边形围成的区域，其边界满足下列三个条件多面体表面上的每一对面，它们或者是分离的、或者有一个公共顶点、或者有一条公共的边如果把多面体看成一个实心的实体，表面上任一点的足够小的邻域同胚于一个三维的如下定义的半开半闭的半球。半开半闭的半球定义为x2+y2+z2r2和z的交，其中r为半球的半径多面体表面是连通的，即表面上任意两点可通过完全在表面上的路径连接,3.1
7、多面体,如果把多面体看成厚度为零的多边形围成的空间，第二个条件也可写为多面体表面上任一点，它在表面上的邻域同胚于一个开圆盘，开圆盘是二维的开圆如果一个表面上的每一个点都满足这个条件，那么这个表面就被称为二维流形(2D manifold),3.1 多面体,第个条件表示顶点和边组成的图是连通的排除那些有“浮在”表面里面的顶点和边的物体,3.1 多面体,多面体可以存在“通孔(hole)”，从多面体表面的一面通到另一面多面体允许有任意多个这样的孔，孔的数目叫做这个体的亏格(genus)亏格为0的多面体在拓扑上同胚于球亏格为1的多面体在拓扑上同胚于圆环面凸多面体又叫做多胞形，为了强调它们是三维的，有时也
8、叫做三胞形多胞形是凸多面体，连接它上面的任意两个点的线段都在多胞形内部凸多边形在每一个顶点处是凸的多胞形在每一条边处的所有二面角(dihedral)都是凸的()二面角是指空间中两个相邻接的面在它们的公共边上的内夹角对于任意的多胞形，顶点处的所有多边形内角之和小于2这是每个顶点处是凸的必要条件，但不是充分条件,3.1 多面体,3.1.2 正则多面体只存在五种不同的正多面体正四面体、正六体、正八面体、正十二面体和正二十面体也叫柏拉图体(Platonic solids)，因为柏拉图在他的蒂迈欧篇 (Timaeus) 中讨论过它们,3.1 多面体,3.1.3 多面体的欧拉公式1758年，Leonard
9、 Euler亏格为0的多面体的顶点数、边数和面数规律顶点数和面数之和总是比边数多2正方体有8个顶点和6个面，8+6=14，比边数12大2用v、e、f分别表示多面体顶点、边和面的数目欧拉公式：v e + f = 2定理3.1 对于一个顶点数v=n，边数与面数分别为e和f的亏格为0的多面体，有 v e + f = 2，且e和f都为O(n),3.2 三维凸包算法,3.2.1 礼品包裹算法礼品包裹算法适用于任意维点集的凸包构建三维礼品包裹算法是二维的直接扩展二维礼品包裹算法从凸包的一条边开始三维礼品包裹算法从凸包的一个面f开始选择面f上的一条边e将f所在的平面，沿着e朝着点集折叠，直至碰到第一个点p，
10、则p，e成为凸包的一个新的三角面片重复上述操作,3.2 三维凸包算法,3.2.1 礼品包裹算法算法起始需要确定凸包上的一个面f首先把点集中的所有三维点投影到yz坐标平面上得到一个二维点集然后求这个点集在二维空间中的一条极端边，要求该边的一个端点是z坐标最小的一个点设这极端边的两个端点对应的三维点为pi和pj，则pipj为三维空间中点集凸包的一条边构造一平面通过pipj且其法线垂直x坐标轴再对绕着pipj旋转，直至碰到第一个点p则p，pipj便是平面f,3.2 三维凸包算法,3.2.1 礼品包裹算法确定一个面片需要O(n)时间设F为凸包上面片的数目三维礼品包裹算法时间复杂度为O(nF)三维礼品包裹算法也具有输出大小敏感特性由定理3.1得知，n个顶点的多面体面片数是O(n)凸包面片的数目也有可能达到O(n)最坏情况下，三维礼品包裹算法的时间复杂度O(n2),3.2 三维凸包算法,3.2.2 分而治之算法分而治之算法的时间复杂度下界也为(nlogn),
