怎樣判斷線段相交

2021-06-05 18:29:18 字數 3345 閱讀 1970

一.向量基本知識

因為後面的計算需要一些向量的基本知識,這裡只是簡單的列舉如下,如果需要更加詳細的資訊,可以自行搜尋wikipedia或google。

1.向量的概念

:如果一條線段的端點是有次序之分的,我們把這種線段成為有向線段(directed segment)。如果有向線段p1p2的起點p1在座標原點,我們可以把它稱為向量(vector)p2。

2.向量加減法

:設二維向量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 )。

3.向量的叉積

:計算向量叉積是與直線和線段相關演算法的核心部分。設向量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共線,但可能同向也可能反向。

4.折線段的拐向判斷

:折線段的拐向判斷方法可以直接由向量叉積的性質推出。對於有公共端點的線段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三點共線。

為了後文的敘述方便,先定義幾個結構:

struct point;

struct v;

計算兩條直線的叉積(cross production), 這裡由於定義的都是二維的情況,本質上說,在平面上兩個向量的叉積應該是垂直平面的,這裡函式返回的整數值即為z軸上的值:

int crossproduct(v* v1, v* v2)

二.判斷兩條直線相交

先來看乙個簡單的情況,即判斷兩條直線是否相交。

第乙個可能會想到的辦法,就是判斷斜率,這個在中學時代就學過了,不過斜率需要考慮垂直的特殊情況,比較麻煩。更好的辦法或許是計算兩個向量的叉積,如果為0,則是平行或者重合的,否則兩直線相交。

**就不貼了,直接呼叫上面的函式就ok了。

三.判斷兩線段相交

經典方法,就是跨立試驗了,即如果一條線段跨過另一條線段,則線段的兩個端點分別在另一條線段的兩側。但是,還需要檢測邊界情況,即兩條線段中可能某條線段的某個端點正好落在另一條線段上。這也是演算法導論中介紹的演算法。

程式模擬如下:

int direction(point* pi, point* pj, point* pk)

int onsegment(point* pi, point* pj, point* pk)

else

if (pi->y > pj->y)

else

if (minx <= pk->x && pk->x <= maxx && miny <= pk->y && pk->y <= maxy)

return 1;

else

return 0;

}int segmentintersect(point* p1, point* p2, point* p3, point* p4)

實際上,如果想改進上述演算法,還可以在跨立試驗前加一步,就是先做快速排斥試驗。那就是,先分別判斷以兩條線段為對角線的矩形是否相交,如果不相交,則兩個線段肯定不相交。

四.判斷兩條線段相交,然後計算交點

設一條線段為l0=p1p2, 另一條線段或直線為l1=q1q2, 要計算的就是l0和l1的交點。

1.首先判斷l0和l1是否相交(方法已在前文討論過), 如果不相交則沒有交點, 否則說明l0和l1一定有交點, 下面就將l0和l1都看作直線來考慮. 

2.如果p1和p2橫座標相同, 即l0平行於y軸 

a)若l1也平行於y軸 

i.若p1的縱座標和q1的縱座標相同, 說明l0和l1共線, 假如l1是直線的話他們有無窮的交點, 假如l1是線段的話可用"計算兩條共線線段的交點"的演算法求他們的交點(該方法在前文已討論過);

ii.否則說明l0和l1平行, 他們沒有交點; 

b)若l1不平行於y軸, 則交點橫座標為p1的橫座標, 代入到l1的直線方程中可以計算出交點縱座標; 

3.如果p1和p2橫座標不同, 但是q1和q2橫座標相同, 即l1平行於y軸, 則交點橫座標為q1的橫座標, 代入到l0的直線方程中可以計算出交點縱座標; 

4.如果p1和p2縱座標相同, 即l0平行於x軸

a)若l1也平行於x軸,

i.若p1的橫座標和q1的橫座標相同, 說明l0和l1共線, 假如l1是直線的話他們有無窮的交點, 假如l1是線段的話可用"計算兩條共線線段的交點"的演算法求他們的交點(該方法在前文已討論過);

ii.否則說明l0和l1平行, 他們沒有交點; 

b)若l1不平行於x軸, 則交點縱座標為p1的縱座標, 代入到l1的直線方程中可以計算出交點橫座標; 

5.如果p1和p2縱座標不同, 但是q1和q2縱座標相同, 即l1平行於x軸, 則交點縱座標為q1的縱座標, 代入到l0的直線方程中可以計算出交點橫座標; 

6.剩下的情況就是l1和l0的斜率均存在且不為0的情況 

a)計算出l0的斜率k0, l1的斜率k1;

b)如果k1 = k2 

i.如果q1在l0上, 則說明l0和l1共線, 假如l1是直線的話有無窮交點, 假如l1是線段的話可用"計算兩條共線線段的交點"的演算法求他們的交點(該方法在前文已討論過);

ii.如果q1不在l0上, 則說明l0和l1平行, 他們沒有交點.

c)聯立兩直線的方程組可以解出交點來

這個演算法並不複雜, 但是要分情況討論清楚, 尤其是當兩條線段共線的情況需要單獨考慮, 所以在前文將求兩條共線線段的演算法單獨寫出來. 另外, 一開始就先利用向量叉乘判斷線段與線段(或直線)是否相交, 如果結果是相交, 那麼在後面就可以將線段全部看作直線來考慮. 需要注意的是, 我們可以將直線或線段方程改寫為ax+by+c=0的形式, 這樣一來上述過程的部分步驟可以合併, 縮短了**長度, 但是由於先要求出引數, 這種演算法將花費更多的時間.

判斷線段相交

先判斷兩條線段是否不平行 最好同時判斷是否有交點並且不平行,因為浮點運算不精確 然後計算兩條線段的交點。以下是c語言 include include define eps 1e 8 define zero x x 0?x x eps 判兩直線平行 int parallel point u1,poin...

判斷線段相交

1.叉積的計算 為了確定兩個點的相對位置,可以使用叉積 有兩個點p1和p2,如果p1xp2 0,則p1位於p2的順時針方向,p1xp2 0,則p1位於p2的逆時針方向。p1的座標為 x1,y1 p2的座標為 x2,y2 則p1xp2 x1y2 x2y1 確定連續線段左轉還是右轉的問題,也可是使用叉積...

C 判斷線段是否相交

c 判斷線段是否相交 線段是否相交,一種是從幾何上就是判斷兩個線段有沒有交點,還有一種是通過向量叉乘 也就是向量積 來判斷。因為向量叉乘的結果是乙個垂直於原來兩個向量的新向量,可以簡單的理解為垂直於原來兩向量所在平面的向量。我們來看圖。線段是否相交 線段p1p2的p1點 線段p1p2的p2點 線段q...