點到線段最短距離的運算與點到直線的最短距離的運算二者之間存在一定的差別,即求點到線段最短距離時需要考慮參考點在沿線段方向的投影點是否在線段上,若在線段上才可采用點到直線距離公式,如圖1所示。
圖1(a)最短距離為點P與其在線段AB上投影C之間的線段PC
(b)最短距離為點P與端點B(或A)所構(gòu)成的線段PB(或PA)
具體算法主要有以下三種:
該算法直接用高中時所學(xué)習(xí)到的解析幾何知識對點到線段的距離進行求解。其基本思想是先判斷點在線段端點、點在線上等等的特殊情況,逐步的由特殊到一般,當(dāng)忽略點在線段上的特殊情況時,判斷點到線段方向的垂線是否落在線段上的方法是通過比較橫縱坐標(biāo)的方式來判斷,最后把不同的判斷情況用不同的幾何方式來進行處理計算得出結(jié)果。
由上面敘述的基本思路可以知道這種算法雖然很容易理解和接受,但從算法的實用性的角度分析還是有很大的缺點的,首先是算法復(fù)雜,計算量巨大,大量的比較判斷、距離計算、角度計算等等,實際應(yīng)用中往往是需要求由大量線段組成的折線到某點的最短距離,如此用這樣的算法計算量是不能想象的。其次經(jīng)典算法中使用的一些簡化運算的函數(shù)不利于語言的重新包裝,如果想換編程語言的話,就比較麻煩了。
該方法主要是先判斷投影點是否在線段上,投影點在線段延長線上時,最短距離長度為點到端點的線段長度;當(dāng)投影點在線段上時,先使用海倫公式計算三角形面積,再計算出三角形的高,即為最短距離。
運用面積算法求解點到線段最短距離思路很清晰,也很容易理解。從效率方面考慮,比如需要多次計算平方、根號,這對于大量數(shù)據(jù)進行運算是負擔(dān)很重的。求面積就必須把三條邊長全部求出,并且用到的海倫公式也需要進行開方運算,計算過程顯得繁瑣。
矢量算法過程清晰,如果具有一定的空間幾何基礎(chǔ),則是解決此類問題時應(yīng)優(yōu)先考慮的方法。當(dāng)需要計算的數(shù)據(jù)量很大時,這種方式優(yōu)勢明顯。
由于矢量具有方向性,故一些方向的判斷直接根據(jù)其正負號就可以得知,使得其中的一些問題得以很簡單的解決。
用此方法考慮,我們只需要找到向量在
方向上的投影,具體如下:
上面的是方向上的單位向量,其意義是給所求向量確定方向。是的兩個向量的內(nèi)積,且,其中θ為向量AP與AB之間的夾角。是向量長度。
那么即為上圖中線段AC的長度值,不帶有方向性。此數(shù)值與上述表征方向的整體構(gòu)成有大小、有方向的新向量,即為在方向上的投影向量,C為投影點。
根據(jù)得到的,由向量的方向性可知:如果情況是上圖(a)所示,那么0<r<1,;如果是如圖(b)所示的情況,那么r>=1,;如果是如圖(c)所示的情況,那么得到r≤0;
特殊情況如點在線段上、點在端點、點在線段延長線上等等的情況全部適用于此公式,只是作為特殊情況出現(xiàn),無需另作討論。這也是矢量算法思想的優(yōu)勢所在。
故根據(jù)r值的不同,最短距離
C#代碼為:
publicstaticdoublePointToSegDist(doublex,doubley,doublex1,doubley1,doublex2,doubley2)
{
doublecross?=?(x2?-?x1)?*?(x?-?x1)?+?(y2?-?y1)?*?(y?-?y1);
if(cross?<=?0)returnMath.Sqrt((x?-?x1)?*?(x?-?x1)?+?(y?-?y1)?*?(y?-?y1));
doubled2?=?(x2?-?x1)?*?(x2?-?x1)?+?(y2?-?y1)?*?(y2?-?y1);
if(cross?>=?d2)returnMath.Sqrt((x?-?x2)?*?(x?-?x2)?+?(y?-?y2)?*?(y?-?y2));
doubler?=?cross?/?d2;
doublepx?=?x1?+?(x2?-?x1)?*?r;
doublepy?=?y1?+?(y2?-?y1)?*?r;
returnMath.Sqrt((x?-?px)?*?(x?-?px)?+?(py?-?y1)?*?(py?-?y1));
}
后記:編程考驗的不光是寫代碼的能力,更多的是對一個問題的優(yōu)質(zhì)解決方法,點到線段最短距離的求解,在數(shù)學(xué)上可能就是一個公式的問題,編程是深究起來卻有如此巧妙的解決方法,記錄本文方法,更多的是對以后學(xué)習(xí)的一個啟發(fā)!
return?Math.Sqrt((x-px)*(x-px)+(py-y)*(py-y))
轉(zhuǎn)載自:http://blog.sina.com.cn/s/blog_5d5c80840101bnhw.html