【數學公式及推導】根據經緯度計算地面兩點間的距離

1.假設：地球是正球體。地面兩點A和B的經緯度座標分別為(Aj,Aw)和(Bj,Bw),地球半徑R取平均值6371km。

2.建立三維直角座標系：

地球球心為原點O，地軸為Z軸，北極方向為Z軸正方向，赤道平面為X軸和Y軸所在平面，在該平面上地心到零度經線的方向為X軸正方向，根據右手定則確定Y軸正方向。

設點A的三維座標為(Ax,Ay,Az)，點B的三維座標為(Bx,By,Bz)

3.思路：

A、B、O三點所在平面與地球相交形成一個半徑為R的圓，求AB間的地面距離就是求該圓上圓弧AB的長度。可由弧長等於半徑乘以圓心角公式求得。

由於R是確定的，只要獲得OA與OB的夾角θ就可以獲得弧AB的長度。弧AB=R*θ。
角θ可通過向量公式求得：

向量OA*向量OB=|OA||OB|cosθ。
則

cosθ=向量OA*向量OB/|OA||OB|
     =(Ax*Bx+Ay*By+Az*Bz)/R*R

4.用經緯度座標表示三維直角座標：

Ax=R*cosAw*cosAj
Ay=R*cosAw*sinAj
Az=R*sinAw
Bx=R*cosBw*cosBj
By=R*cosBw*sinBj
Bz=R*sinBw

代入可得

cosθ=cosAw*cosAj*cosBw*cosBj+cosAw*sinAj*cosBw*sinBj+sinAw*sinBw
     =cosAw*cosBw(cosAj*cosBj+sinAj*sinBj)+sinAw*sinBw
     =cosAw*cosBw*cos(Aj-Bj)+sinAw*sinBw
θ=arccos[cosAw*cosBw*cos(Aj-Bj)+sinAw*sinBw]

5.綜上可得根據經緯度計算地面兩點間距離的公式：

弧AB=R*arccos[cosAw*cosBw*cos(Aj-Bj)+sinAw*sinBw]

說明：

類似的公式推導大家以前都做過，時間久了可能會忘記一些東西，於是我把它記了下來，以備查閱。

由於三角函數變換多端，該公式的表現形式不止一種，只要前提假設是一樣的，那麼在數學上應該是等價並可相互推導的。
在實際程式程式碼中用的時候需要注意單位問題，比如從定位器獲得的經緯度單位是度，而三角計算的方法很可能用的單位是弧度，再比如長度單位用的是km還是m等。