GIS开发常用算法原理分析

本文详细介绍了计算几何中的基本算法,包括矢量概念、加减法、叉积,以及在GIS开发中涉及的线段、折线、多边形、矩形和圆的判断与操作。通过这些算法,可以判断点、线段是否在其他几何形状内,以及计算最近点和交点。此外,还探讨了凸包的概念及其求解方法。
摘要由CSDN通过智能技术生成

摘要:

计算机的出现使得很多原本十分繁琐的工作得以大幅度简化,但是也有一些在人们直观看来很容易的问题却需要拿出一套并不简单的通用解决方案,比如几何问题。作为计算机科学的一个分支,计算几何主要研究解决几何问题的算法。在现代工程和数学领域,计算几何在图形学、机器人技术、超大规模集成电路设计和统计等诸多领域有着十分重要的应用。在本文中,我们将对计算几何常用的基本算法做一个全面的介绍,希望对您了解并应用计算几何的知识解决问题起到帮助。

 

矢量的概念

如果一条线段的端点是有次序之分的,我们把这种线段成为有向线段(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,y1Q = (x2,y2),则矢量叉积定义为由(0,0)p1p2p1+p2所组成的平行四边形的带符号的面积,即:P× Q = x1*y2 - x2*y1,其结果是一个标量。显然有性质 P× Q = - ( Q× P ) P× ( - Q ) = - ( P× Q )。一般在不加说明的情况下,本文下述算法中所有的点都看作矢量,两点的加减法就是矢量相加减,而点的乘法则看作矢量叉积。

 

叉积的一个非常重要性质是可以通过它的符号判断两矢量相互之间的顺逆时针关系:

 

P× Q > 0 ,PQ的顺时针方向。 P × Q < 0 ,PQ的逆时针方向。 P × Q = 0 ,PQ共线,但可能同向也可能反向。

 

折线段的拐向判断

折线段的拐向判断方法可以直接由矢量叉积的性质推出。对于有公共端点的线段p0p1p1p2,通过计算(p2 - p0) × (p1 - p0)的符号便可以确定折线段的拐向:

 

(p2 - p0)× (p1 - p0) > 0,p0p1p1点拐向右侧后得到p1p2

 

(p2 - p0)× (p1 - p0) < 0,p0p1p1点拐向左侧后得到p1p2

 

(p2 - p0)× (p1 - p0) = 0,p0p1p2三点共线。

 

具体情况可参照下图:

 

判断点是否在线段上

设点为Q,线段为P1P2,判断点Q在该线段上的依据是:( Q- P1 )× ( P2 - P1 ) = 0 Q在以 P1P2为对角顶点的矩形内。前者保证Q点在直线P1P2上,后者是保证Q点不在线段P1P2的延长线或反向延长线上,对于这一步骤的判断可以用以下过程实现:

 

ON-SEGMENT(pi,pj,pk)

 

if min(xi,xj)<=xk<=max(xi,xj) andmin(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)快速排斥试验

 

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

 

(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<

评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

MingoJ

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值