资源描述：
《gis算法计算几何基础》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、GIS算法的计算几何基础矢量的概念：如果一条线段的端点是有次序之分的，我们把这种线段成为有向线段(directedsegment)。 如果有向线段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)。矢量叉积：计算矢量叉积是与直线和线段
2、相关算法的核心部分。设矢量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共线，但可能同
3、向也可能反向。折线段的拐向判断： 折线段的拐向判断方法可以直接由矢量叉积的性质推出。对于有公共端点的线段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，判断
4、点Q在该线段上的依据是：(Q-P1)×(P2-P1)=0且Q在以P1，P2为对角顶点的矩形内。前者保证Q点在直线P1P2上，后者是保证Q点不在线段P1P2的延长线或反向延长线上，对于这一步骤的判断可以用以下过程实现： ifmin(Xp1,Xp2)<=Xq<=max(Xp1,Xp2)andmin(Yp1,Yp2)<=Yq<=max(Yp1,Yp2) {returntrue}; else {returnfalse;} 特别要注意的是，由于需要考
5、虑水平线段和垂直线段两种特殊情况，min(xi,xj)<=xk<=max(xi,xj)和min(yi,yj)<=yk<=max(yi,yj)两个条件必须同时满足才能返回真值。 判断两线段是否相交： 我们分两步确定两条线段是否相交： (1)快速排斥试验 设以线段P1P2为对角线的矩形为R，设以线段Q1Q2为对角线的矩形为T，如果R和T不相交，显然两线段不会相交。 (2)跨立试验 如果两线段相交，则两线段必然相互跨立对方。若P1P2跨立Q1Q2，则矢量(P1-Q1)和(P2-Q1)位于矢量(Q
6、2-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
7、的依据是：(Q1-P1)×(P2-P1)*(P2-P1)×(Q2-P1)>=0。 判断线段和直线是否相交：如果线段P1P2和直线Q1Q2相交，则P1P2跨立Q1Q2，即：(P1-Q1)×(Q2-Q1)*(Q2-Q1)×(P2-Q1)>=0。 判断矩形是否包含点：只要判断该点的横坐标和纵坐标是否夹在矩形的左右边和上下边之间。 判断线段、折线、多边形是否在矩形中：因为矩形是个凸集，所以只要判断所有端点是否都在矩形中就可以了。 判断矩形是否在矩形中： 只要比较左右边界和上下边界就可以了。 判断圆是否在矩
8、形中： 很容易证明，圆在矩形中的充要条件是：圆心在矩形中且圆的半径小于等于圆心到矩形四边的距离的最小值。 判断点是否在多边形中： 判断点P是否在多边形中是计算几何中一个非常基本但是十分重要的算法。以点P为端点，向左方作射线L，由于多边形是有界的，所以射线L的左端一定在多边形外，考虑沿着L从无穷远处开始自左向右移动，遇到和多边形的第一个交点的时候，进入到了多边形的内部，遇到第二个交点的时候，离开
