GIS算法的计算几何基础1综述.ppt

* 5 判断两直线相交 算法1： （1）快速排除：以两条直线为端点的矩形不相交。（方法？）若矩形不相交，则直线不会相交。 * 5 判断两直线相交 （2）跨立试验：如果两线段相交，则必然跨立对方。即一直线的两端点必然位于另一直线两侧。 算法2： 定义A,B,C,D为二维空间点，则有向线段AB和CD的参数方程为： * 5 判断两直线相交 如果AB与CD相交，则： 解方程得： 设P为直线AB和CD的交点，则： * 5 判断两直线相交 如果 且 ，则有向线段AB与CD相交。 如果(Bx-AX)(Dy-Cy)-(By-Ay)(Dx-CX)=0，则AB与CD平行。 如果(By-Ay)(Dx-Cx)-(Bx-Ax)(Dy-Cy)=0，则AB与CD共线。 如果直线AB和CD相交，而交点不位于线段AB和CD之间，则交点位置可通过如下条件判断： r1，则P位于有向线段AB的延长线上； r0，则P位于有向线段BA的延长线上； s1，则P位于有向线段CD的延长线上； s0，则P位于有向线段DC的延长线上； * 6 矩形是否包含点 只要判断点的横坐标与纵坐标是否夹在矩形的左右边和上下边之间。 * 7判断线段、折线、多边形是否在矩形中 矩形是凸集，所以只需判断所有的端点是否在矩形中。 * 8 判断矩形是否在矩形中 比较左右边界和上下边界 * 9 判断圆是否在矩形中 圆心在矩形中，且圆的半径小于等于圆心到矩形四边的距离的最小值。 * 常用算法： ①射线法（奇偶法） ②转角法：环绕数（多边形环绕点的次数）为0，则点在多边形外，否则，点在多边形内。 10 判断点是否在多边形内 * 10.1 射线法 射线法计算从点P开始的射线穿过多边形边界的次数，多边形的边界将多边形内部与外部分开。如果结果为偶数，则点在多边形外部，否则，若为奇数，则点在多边形内部。 * 10.1 射线法 必须决定在多边形边界上的点是在多边形内部还是外部。一个标准的约定是在左边界或者下边界上的点认为是在多边形内部，在右边界或者上边界的点认为是在多边形外部。在这种约定下，如果两个不同的多边形共享一个边，那么在这条边上的点会在一个多边形的内部而在另一个多边形的外部。 * 10.1 射线法 一个简单的射线法是选择一条水平的、向点P的右侧延伸的、平行于X轴的射线。 为了计算总的交点的数目，算法简单的遍历多边形的所有边，测试是否穿越边，穿越时增加交点的数目。 穿越测试必须处理好一些特殊的情形。完全规则如下： （1）方向向上的边包括其开始点，不包括其终止点； （2）方向向下的边不包括其开始点，包括其终止点； （3）水平边不参加穿越测试； （4）射线和多边形的边的交点必须严格在点P的右边。 * 10.1 射线法 * 10.1 射线法 射线法的有效性是基于一个简单的闭合曲线将二维平面分成两个相连接的部分：有边界的内部和无边界的外部。 曲线必须是简单的（没有自相交），否则平面会被分成两个以上的部分，并且不能保证穿越边界时不会改变出入特性。 对于一个闭合的曲线，可能将二维平面分成多个部分，但是其中只有一个没有边界且在曲线外部的部分。 有边界的部分可能在曲线内部也可能在曲线外部。 两个有共同边界的有边界部分可能都在曲线内部，那么穿越过共享的边界并不会改变出入特性。 * 10 判断点是否在多边形内 判断出现错误！ * 10.2 转角法 转角法能够很精确的判定一个点是否在多边形内部。需要计算多边形绕点多少次，环绕数为0时，点在多边形外部。 * 10.2 转角法 方法：定义二维平面中某个封闭曲线关于某点的环绕数。令C(u)=(x(u),y(u))，0≤u≤1，且C(0)=C(1)，是二维连续曲线。P是不在C(u)上的点。 令Cp(u)= C(u)-P为从点P到C(u)的矢量，并且它的单位方向矢量为W(u)=CP(u)/| CP(u)|。 W(u)坐标形式为： W(u)=(cos(u),sin(u))， 其中(u)是正的逆时针方向的角。环绕的次数(wn)就等于W(u)环绕S1的次数的整数部分。计算公式为： * 10.2 转角法 若曲线C是由顶点V0,V1,…,Vn=V0构成的多边形，那么积分可以简化为计算带符号的角度的总和。这些角PVi与PVi+1的夹角。因此，如果i=angle(PVi,PVi+1),则 这个公式效率不高，原因在于arccos函数很耗时。 改进：在S1中任取一点Q，因为曲线W(u)环绕，则可能多次经过Q点。当W(u)按逆时针方向经过Q点时，记为+1次，顺时针经过Q点时，记为-1次，那么总次数和就是W环绕S1的次数，刚好等于环绕数（wn）。 * 10.2 转角法 用一个射线R，R的起点