凸包算法设计.pptVIP

凸包问题算法设计 制作者：地信*** *** 2014.05.10 凸包概念 1.点集Q的凸包（convex hull）是指一个最小凸多边形，满足Q中的点或者在多边形边上或者在其内。右图中由红色线段表示的多边形就是点集Q={p0,p1,...p12}的凸包。 2.一组平面上的点，求一个包含所有点的最小的凸多边形，这就是凸包问题了。这可以形象地想成这样：在地上放置一些不可移动的木桩，用一根绳子把他们尽量紧地圈起来，并且为凸边形，这就是凸包了。 算法：增量式算法,包裹法（Jarvis步进法）,葛立恒（Graham）扫描法,单调链,分治法,快包法（Akl-Toussaint启发式）等。 1.1 分治法的概念 所谓分治法就是把问题划分成多个子问题来进行处理。这些子问题，在结构上与原来的问题一样，但在规模上比原来的小。如果得到的子问题相对来说还大，可以反复地使用分治策略，把这些子问题再划分成更小的、结构相同的子问题。这样就可以使用递归的方法，分别求解这些子问题，把这些子问题的解结合起来，从而获得原来问题的解。 化整为零，各个击破。 1.2 问题分析图（PFD图） 1.3 过程描述 一般来说，分治法的求解过程由以下三个阶段组成： （1）划分：把规模为n 的原问题划分为k 个规模较小且 规模大致相同的子问题。 （2）求解子问题：使用递归方法分别求解子问题。 （3）合并：把各子问题的解合并起来，合并的代价因情况而异，分治算法的有效性很大程度上依赖于合并的实现。 1.4 凸包问题的分治思想： 第一步：把给定点集中的点在横坐标方向上按照大小排序。 如下图所示，p1 和pn 必定是凸多边形的两个顶点。 第二步：在上凸包点集合s1 中找到一个距离直线最远点pmax，如下 图所示。显然直线段p1pmax与直线段pmaxpn把点集s1分成了三个集合。由凸包的性质可知p1，pmax，pn三点围成的三角形中的点不可能作为凸包的顶点，所以只需考虑直线p1pmax左边的点s11以及直线pmaxpn右边的点s12。 第三步：递归求解得到凸多边形的边。 第四步：合并这些边即得所求凸包。 1.5 凸包问题的分治算法 算法： convex_divide() findtwo(MIN_Y,MAX_Y); //查找Y值最小和最大的两个点，作为初始凸包顶点 initialmain(leftList,rightList,side,MIN_Y,MAX_Y);//由这两点确定的线段side ，将平面的点分成相互独立的两部分，左边区域内的散点放在leftList中，右边的散点放在rightList中 dealwith(side,leftList,resultList); dealwith(side,rightList,resultList); write(resultList); //输出凸包的顶点 dealwith(VAR side,list， resultList) //处理某一条边side及其所属的离散点集list，并把找到的凸包顶点放到凸包顶点集resultList中 if not is Empty(list) then { node:=LongNode(side,list);//找出list中离这条边最远的点，并把它从List中删除 insert(resultList,side, node);//该点肯定是凸包的顶点，把它放到凸包顶点集的相应位置，该位置就在该边起点的后面 Triangle:=create Triangle(side,node); //该点和原先的边生成一个三角形 while(not is Empty(list)) do { Node:=GetNextNode(list); //从list中取出一个点，进行归边处理 If in Triangle(node,Triangle) //如果点在三角形内，则抛弃 then delete(node) else //点在三角形外，则判断其所属的边 if (attachSide(node,Lside,Rside)=left) then insert(LeftList,node); //该点属于Lside边放入leftList散点集中 else insert(RightList,node); //处理属于Rside 边，放入rightList 集中 }//list 中所有点处理完毕 dealwith(Lside,Leftlist,resultList);// 处理属于Lside 边及其所属的点 dealwith(rsid