计算机几何题算法.docVIP

计算几何相关算法介绍一 2009-08-03 14:55 【声明】本文内容均源于互联网。 算法介绍 　　矢量的概念： 　　如果一条线段的端点是有次序之分的，我们把这种线段称为有向线段(directed segment)。如果有向线段p1p2的起点p1在坐标原点，我们可以把它称为矢量(vector)p2。 　　矢量加减法： 　　设 二维矢量P = ( x1, y1 )，Q = ( x2 , y2 )，则矢量加法定义为： P + Q = ( x1 + x2 , y1 + y2 )，同样的，矢量减法定义为： P - Q = ( x1 - x2 , y1 - y2 )。显然有性质 P + Q = Q + P，P - Q = - ( Q - P )。 　　矢量叉积： 　　计算矢量叉积是与直线和线段相关算法的核心部分。设矢量P = ( x1, y1 )，Q = ( x2, y2 )，则矢量叉积定义为由(0,0)、p1、p2和p1+p2所组成的平行四边形的带符号的面积，即：P × Q = x1*y2 - x2*y1，其结果是一个标量。显然有性质 P × Q = - ( Q × P ) 和 P × ( - Q ) = - ( P × Q )。一般在不加说明的情况下，本文下述算法中所有的点都看作矢量，两点的加减法就是矢量相加减，而点的乘法则看作矢量叉积。 　　叉积的一个非常重要性质是可以通过它的符号判断两矢量相互之间的顺逆时针关系： 　　若 P × Q 0 , 则P在Q的顺时针方向。　　若 P × Q 0 , 则P在Q的逆时针方向。　　若 P × Q = 0 , 则P与Q共线，但可能同向也可能反向。 　　折线段的拐向判断： 　　折线段的拐向判断方法可以直接由矢量叉积的性质推出。对于有公共端点的线段p0p1和p1p2，通过计算(p2 - p0) × (p1 - p0)的符号便可以确定折线段的拐向： 　　若(p2 - p0) × (p1 - p0) 0,则p0p1在p1点拐向右侧后得到p1p2。 　　若(p2 - p0) × (p1 - p0) 0,则p0p1在p1点拐向左侧后得到p1p2。 　　若(p2 - p0) × (p1 - p0) = 0,则p0、p1、p2三点共线。 　　具体情况可参照下图： 判断点是否在线段上： 　　设点为Q，线段为P1P2 ，判断点Q在该线段上的依据是：( Q - P1 ) × ( P2 - P1 ) = 0 且 Q 在以 P1，P2为对角顶点的矩形内。前者保证Q点在直线P1P2上，后者是保证Q点不在线段P1P2的延长线或反向延长线上，对于这一步骤的判断可以用以下过程实现： 　　ON-SEGMENT(pi,pj,pk) 　　if min(xi,xj) = xk = max(xi,xj) and min(yi,yj) = yk = max(yi,yj) 　　then return true; 　　else return false; 　　特别要注意的是，由于需要考虑水平线段和垂直线段两种特殊情况，min(xi,xj)=xk=max(xi,xj)和min(yi,yj)=yk=max(yi,yj)两个条件必须同时满足才能返回真值。 1. 矢量减法设二维矢量 P = （x1,y1） ，Q = (x2,y2)则矢量减法定义为： P - Q = ( x1 - x2 , y1 - y2 )显然有性质 P - Q = - ( Q - P )如不加说明，下面所有的点都看作矢量，两点的减法就是矢量相减；2.矢量叉积设矢量P = （x1,y1） ，Q = (x2,y2)则矢量叉积定义为： P × Q = x1*y2 - x2*y1?? 得到的是一个标量显然有性质 P × Q = - ( Q × P )?? P × ( - Q ) = - ( P × Q )如不加说明，下面所有的点都看作矢量，点的乘法看作矢量叉积；叉乘的重要性质：??? 若 P × Q 0 , 则P 在Q的顺时针方向??? 若 P × Q 0 , 则P 在Q的逆时针方向??? 若 P × Q = 0 , 则P 与Q共线，但可能同向也可能反向3.判断点在线段上设点为Q，线段为P1P2 ，判断点Q在该线段上的依据是：( Q - P1 ) × ( P2 - P1 ) = 0 且 Q 在以 P1，P2为对角顶点的矩形内4.判断两线段是否相交我们分两步确定两条线段是否相交：(1)．?? 快速排斥试验设以线段 P1P2 为对角线的矩形为R， 设以线段 Q1Q2 为对角线的矩形为T，如果R和T不相交，显然两线段不会相交；(2)．?? 跨立试验如果两线段相交，则两线段必然相互跨立对方，如图1所示。在图1中，