设地球的经纬度坐标是(x,y),x是经度,y是纬度
先把它转换成直角坐标(X,Y,Z)
设地球半径为r
则有转换式:
X=r*cos(y)*cos(x)
Y=r*cos(y)*sin(x)
Z=r*sin(y)
分别得到A(x1,y1),B(x2,y2)两点的直角坐标A(X1,Y1,Z1),B(X2,Y2,Z2)
那么AB两点间直线距离是L=√[(X2-X1)^2+(Y2-Y1)^2+(Z2-Z1)^2]
把上面转换式代入化简可以得到L=r*√[2-2cos(y1)*cos(y2)*cos(x1-x2)-2sin(y1)*sin(y2)]
设AB两点对球心所夹圆心角为θ(弧度),则有L=2r*sin(θ/2)=2r*√[(1-cosθ)/2]=r*√[2*(1-cosθ)]
因此cosθ=1-L^2/(2r^2)=cos(y1)*cos(y2)*cos(x1-x2)+sin(y1)*sin(y2)
那么弧长=r*θ
=r*arccos[cos(y1)*cos(y2)*cos(x1-x2)+sin(y1)*sin(y2)]
这里的经纬度用角度弧度均可,因为都要代入三角函数,但arccos得到的θ必须用弧度
追问这个回答我已经看过了,持怀疑态度啊。算出来不准确。方位角是与正北方的方位角。此处的
θ是圆心角?感谢能给予回答,希望能有朋友再次解答。