使用WGS84-椭圆体的距离

Question

考虑地球表面上的点P1(60°N,20°E,0)和P2(60°N,22°E,0)

当使用WGS-84椭圆体模拟地球的形状时,点P1和P2之间的最短距离是多少？

Answer 1

不幸的是，Vincenty 的算法无法对某些输入收敛。 GeographicLib提供了一种始终收敛的替代方案（并且也更准确）。提供了在 C++、C、Fortran、Javascript、Python、Java 和 Matlab 中的实现。例如，使用 Matlab 包：

format long;
geoddistance(60,20,60,22)
->
111595.753650629

Answer 2

常用半正矢公式（误差 < 0.5%）