一、叉积
叉积的计算是线段方法的核心。对于向来p1和p2,叉积是由点(0,0)、p1、p2和p1+p2构成的平行四边形的有向面积。另一种与之等价但更有效的的叉积定义方式是将其看做矩阵行列式:
p1×p2 = x1y2 - x2y1 = - p2×p1
若p1×p2为正,则相对于原点(0,0)来说,p1位于p2顺时针方向;若p1×p2为负,p1位于p2逆时针方向;若为0则方向相同,或相反。
若是相对于点p0(x0,y0)而非原点,则p0p1和p0p2的叉积为(p1-p0)×(p2-p0) = (x1-x0)(y2-y0)-(x2-x0)(y1-y0)。
确定连续线段是向左转还是向右转
对于线段p0p1和p1p2,采用叉积可以避免计算角度,只需简单的计算一下p0p2是位于p0p1的顺时针还是逆时针方向。计算叉积****
(p2-p0)×(p1-p0) = (x2-x0)(y1-y0) - (x1-x0)(y2-y0)
若结果为负,p0p2在p0p1的逆时针方向,在p1处左转;结果为正则右转;为0表示三点共线。
判断两条线段是否相交
要判断两条线段是否相交,则需要检查每条线段是否跨越了另一条线段的直线。如果点p1位于某直线的一边,而点p2位于该直线的另一边,则称p1p2跨越了这条直线。两条线段相交,当且仅当下面两个条件至少成立一个:****
每条线段都跨越了包含另一条线段的直线
一条线段的一个端点落在另一条线段上
二、确定任意一对线段是否相交
给定一个线段的集合,仅仅判断是否有两个线段相交,不必输出所有相交的线段对。此处我们假设,线段均不垂直。
我们使用“扫除”算法来解决这个问题。在扫除过程中,一条假想的扫除线穿过一个给定的几何物体集合,会与集合中的部分线段相交。下图中扫除线r与线段a、c相交,且与a的交点的y坐标值大于c的,则认为此处a>c。
在图(a)中,在r处,a>c;在t处,a还是>c,那么我们认为a和c没有相交。而在图(b)中,在v处,e>f,而在图w处,f>e,故e和f相交。
上图描述了一个算法,我们按照线段端点的x坐标,从小到达,进行“扫除”。当在一个线段的左端点扫除时,将所有相交的线段序列按序放入一个“完全前序关系T”中。当扫到线段的右端点时,从T中去除该线段。但是,如果位于在该线段上方的线段集合,与位于该线段下方的线段集合有交集,则表明有线段相交。
关键在于,如果存储T,如果求交集。可以用红黑树来存储T。
算法步骤
- 初始化T为空集
- 对线段的端点排序
- 在端点p处,开始扫除,从最左边的处开始
如果p是线段s的左端点
Insert(T, s);
如果有线段在扫除线处相交
return true;
如果p是线段的右端点
如果above(T, s) 和 below(T, s)有交集,则
return true;//即线段集中存在线段相交
无交集则,Delete(T, s); - return false
三、寻找凸包
点集Q的凸包,是一个最小的凸多边形P,满足Q中的每个点都在P的边界上,或者在P的内部。****
Graham扫描法: 复杂度O(nlogn)
选取y最小的点,多个y最小的话,选取其中x最小的点,作为p0****
剩余的点,按照p0和pi的极角的逆时针排序,编号为p1,p2,...,pm
如果m小于2,表示点数小于3,形不成多边形
设定以空栈S,将p0、p1、p2压入栈中。
for i=3 to m
得到栈顶的2个点pi-1和pi-2,如果t1t0向t0pi转的时候,不是左转,就把顶点t0出栈;
如果出栈了t0,就继续a,直到栈的顶点不再出栈位置
将pi入栈
return S
图中,从a到f是一步一步选择的过程。
Jarvis步进法:复杂度O(nh),h是凸包顶点数
先找到最下边结点里最左边的点p0,然后寻找使得p0p1极角最小的点,则p1也是凸包顶点;继续寻找使得p1p2极角最小的点,直到达到最高点pk,上图是p3,此时已经构造好了CH(Q)的右链。为了构造左链,寻找pk+1使得pkpk+1极角最小,但此时x轴啊原x轴的负方向。