已知两点大地坐标，基于 Vincenty 反算公式精确求解大地方位角与椭球面距离。适用于导线测量闭合差计算、GPS 基线向量解算与大地控制网精度评定。
基于 Vincenty (1975) 反算公式，精度达 ±0.5 mm
不同椭球参数会导致计算结果产生系统性偏差，请根据实际工程需求选择
十进制度数（DD），北纬为正，南纬为负
十进制度数（DD），东经为正，西经为负
深入理解 Vincenty 反算公式的数学基础与计算流程
Vincenty 反算公式的核心思想是将椭球面上的大地线问题转换为辅助球面上的大圆弧问题。通过引入辅助变量（如修正纬度 U、球面角距 σ 等），利用迭代法逐步逼近真实解。
与正算公式不同，反算公式无法直接得到闭合解，必须通过迭代计算经差修正量 λ，直至前后两次迭代值之差小于给定阈值（通常为 1×10⁻¹² 弧度）。
步骤 1：计算辅助纬度 U₁、U₂（归化纬度）：
U = arctan[(1 − f) · tan(B)]步骤 2：计算初始经差 L = L₂ − L₁，设 λ = L。
步骤 3：迭代计算 λ 直至收敛：
sin σ = √[(cos U₂ · sin λ)² + (cos U₁ · sin U₂ − sin U₁ · cos U₂ · cos λ)²]
cos σ = sin U₁ · sin U₂ + cos U₁ · cos U₂ · cos λ
σ = arctan2(sin σ, cos σ)
sin α = (cos U₁ · cos U₂ · sin λ) / sin σ
cos²α = 1 − sin²α
cos(2σₘ) = cos σ − (2 · sin U₁ · sin U₂) / cos²α
C = (f/16) · cos²α · [4 + f · (4 − 3 · cos²α)]
λ′ = L + (1−C) · f · sin α · {σ + C · sin σ · [cos(2σₘ) + C · cos σ · (−1 + 2·cos²(2σₘ))]}
当 |λ − λ′| < 1×10⁻¹² 时停止迭代。
迭代收敛后，计算以下辅助量：
u² = cos²α · (a² − b²) / b² = cos²α · [a²/b² − 1]
A = 1 + (u²/16384) · [4096 + u² · (−768 + u² · (320 − 175·u²))]
B = (u²/1024) · [256 + u² · (−128 + u² · (74 − 47·u²))]
Δσ = B·sin σ · {cos(2σₘ) + (B/4)·[cos σ·(−1+2·cos²(2σₘ)) − (B/6)·cos(2σₘ)·(−3+4·sin²σ)·(−3+4·cos²(2σₘ))]}
最终椭球面距离：
S = b · A · (σ − Δσ)大地方位角：
A₁₂ = arctan2(cos U₂ · sin λ, cos U₁ · sin U₂ − sin U₁ · cos U₂ · cos λ)
A₂₁ = arctan2(cos U₁ · sin λ, −sin U₁ · cos U₂ + cos U₁ · sin U₂ · cos λ) + π
| 算法 | 精度 | 适用范围 | 计算复杂度 |
|---|---|---|---|
| Haversine | ~0.5% | 短距离近似 | O(1) |
| 球面余弦 | ~0.5% | 短距离近似 | O(1) |
| Andoyer-Lambert | ~10 m | 中距离 | O(1) |
| Thomas (1965) | ~1 m | 全距离 | O(1) |
| Vincenty (1975) | ~0.5 mm | 全距离（非对跖） | O(n) 迭代 |
| Karney (2013) | ~nm 级 | 全距离（含对跖） | O(n) 迭代 |
注：Karney (2013) 提出的测地线算法基于椭圆积分，可处理对跖点情况，精度可达纳米级，但计算更为复杂。本工具采用 Vincenty 公式，已满足绝大多数工程与科研需求。
不同坐标系采用的参考椭球几何参数一览
| 椭球名称 | 建立年份 | 长半轴 a (m) | 扁率 f | 短半轴 b (m) | 第一偏心率 e² | 主要使用范围 |
|---|---|---|---|---|---|---|
| WGS-84 | 1984 | 6,378,137.000 | 1 / 298.257223563 | 6,356,752.3142 | 0.00669437999013 | 全球 GPS 定位 |
| CGCS2000 | 2000 | 6,378,137.000 | 1 / 298.257222101 | 6,356,752.3141 | 0.00669438002290 | 中国国家大地坐标系 |
| GRS-80 | 1980 | 6,378,137.000 | 1 / 298.257222101 | 6,356,752.3141 | 0.00669438002290 | 国际大地测量学推荐 |
| 西安 80 | 1980 | 6,378,140.000 | 1 / 298.257 | 6,356,755.2882 | 0.00669438499959 | 中国历史大地坐标系 |
| 北京 54 | 1954 | 6,378,245.000 | 1 / 298.3 | 6,356,863.0188 | 0.00669342162297 | 中国历史大地坐标系 |
| 克拉索夫斯基 | 1940 | 6,378,245.000 | 1 / 298.3 | 6,356,863.0188 | 0.00669342162297 | 苏联及东欧历史测量 |
| 贝塞尔 | 1841 | 6,377,397.155 | 1 / 299.1528128 | 6,356,078.963 | 0.00667437223180 | 德国、奥地利、日本历史测量 |
| 国际 1924 | 1924 | 6,378,388.000 | 1 / 297.0 | 6,356,911.946 | 0.00672267002233 | 海福特椭球，北美历史测量 |
数据来源：IUGG、国家测绘地理信息局、NGS（美国国土地测量局）公开技术文档