一些计算几何

1）向量的外积

    计算矢量叉积是与直线和线段相关算法的核心部分。设矢量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共线，但可能同向也可能反向。

 

 2）判断两条线段是否相交，分两步判断：

    (1)快速排斥试验

    设以线段 P1P2 为对角线的矩形为R， 设以线段 Q1Q2 为对角线的矩形为T，如果R和T不相交，显然两线段不会相交。

     (2)跨立试验
    如果两线段相交，则两线段必然相互跨立对方。若P1P2跨立Q1Q2 ，则矢量 ( P1 - Q1 ) 和( P2 - Q1 )位于矢量( Q2 - Q1 ) 的两侧，即( P1 - Q1 ) × ( Q2 - Q1 ) * ( P2 - Q1 ) × ( Q2 - Q1 ) < 0。上式可改写成( P1 - Q1 ) × ( Q2 - Q1 ) * ( Q2 - Q1 ) × ( P2 - Q1 ) > 0。当 ( P1 - Q1 ) × ( Q2 - Q1 ) = 0 时，说明 ( P1 - Q1 ) 和 ( Q2 - Q1 )共线，但是因为已经通过快速排斥试验，所以 P1 一定在线段 Q1Q2上；同理，( Q2 - Q1 ) ×(P2 - Q1 ) = 0 说明 P2 一定在线段 Q1Q2上。所以判断P1P2跨立Q1Q2的依据是：( P1 - Q1 ) × ( Q2 - Q1 ) * ( Q2 - Q1 ) × ( P2 - Q1 ) >= 0。同理判断Q1Q2跨立P1P2的依据是：( Q1 - P1 ) × ( P2 - P1 ) * ( P2 - P1 ) × ( Q2 - P1 ) >= 0。

    3）已知多边形各个点的坐标，计算任意多边形面积。

    计算p[0],p[i],p[i+1]三点构成的三角形面积(i从1到最后一个点之前)。

     A（x1,y1）B（x2,y2）C（x3,y3）三点的面积公式为   
                           |x1     x2       x3|   
    S(A,B,C)   =     |y1     y2     y3|   *   0.5     (当三点为逆时针时为正，顺时针则为负的)   
                           |1       1       1   |