[论文解读] Geodesics on an ellipsoid of revolution

Charles F. F. Karney|arXiv (Cornell University)|Mar 1, 2009

Historical Geography and CartographySocial Sciences参考文献 25被引用 21

一句话总结

本文提出了一套高精度、非迭代的算法，用于求解旋转椭球面上的直接与反向大地测量问题，精度达到亚纳米级别。该文引入了一种基于辅助球面与雅可比椭圆函数的新公式，可实现对三角测量、地图投影、海洋边界划分以及多边形面积计算的稳健求解，确保收敛性与高数值稳定性。

ABSTRACT

Algorithms for the computation of the forward and inverse geodesic problems for an ellipsoid of revolution are derived. These are accurate to better than 15 nm when applied to the terrestrial ellipsoids. The solutions of other problems involving geodesics (triangulation, projections, maritime boundaries, and polygonal areas) are investigated.

研究动机与目标

• 开发一种数值稳定且高精度的旋转椭球面上直接与反向大地测量问题的解法。
• 通过消除对固定阶数扁平率的依赖并修复几乎对跖点情况下的收敛失败问题，克服现有算法（尤其是Vincenty方法）的局限性。
• 将大地测量解法推广至长扁椭球体，并拓展其在三角测量、地图投影、海洋边界划分以及球面多边形面积计算中的应用。
• 为GeographicLib提供一个全面的理论与计算框架，整合缩短长度等高级大地测量概念。

提出的方法

• 采用辅助球面法表述大地线，通过共形映射将椭球问题转化为可解的球面三角形问题。
• 利用雅可比椭圆函数（sn, dn, am）及其积分（F, E, G）精确表达大地线弧长、方位角与位置。
• 通过以扁平率 $f$、第三扁平率 $n$ 或偏心率平方 $e^2$ 为变量的级数展开，非迭代求解直接问题，确保精度优于15 nm。
• 采用牛顿法并辅以精心选择的初始猜测值求解反向问题，即使在几乎对跖点情况下也能在数次迭代内保证收敛。
• 引入缩短长度 $m_{12}$ 作为大地三角形的关键不变量，实现对三角测量与面积计算的高精度处理。
• 通过使用大地线距离与方位角，将平面投影与球面多边形面积公式推广至椭球情形，保持高精度。

实验结果

研究问题

• RQ1 如何实现对旋转椭球面上直接与反向大地测量问题的亚纳米级精度求解，并确保收敛性？
• RQ2 将大地测量算法推广至长扁椭球体需要哪些修改？其解法与扁平椭球体有何不同？
• RQ3 如何基于大地线方法高精度地计算球面多边形面积并定义海洋边界？
• RQ4 椭圆积分与雅可比函数在实现直接大地测量问题非迭代求解中起到何种作用？
• RQ5 缩短长度概念如何提升三角测量与大地网计算的精度与稳定性？

主要发现

• 通过在 $f$、$n$ 或 $e^2$ 上展开的级数方法，非迭代求解直接大地测量问题，精度优于15 nm，适用于扁平与长扁椭球体。
• 反向大地测量问题通过牛顿法实现可靠收敛，即使在几乎对跖点情况下也仅需数次迭代，克服了Vincenty算法的关键缺陷。
• 缩短长度 $m_{12}$ 是大地三角网中一个基本不变量，可实现高精度的三角测量与面积计算。
• 对于长扁椭球体，经度差 $\lambda_{12}$ 与初始方位角 $\alpha_1$ 之间不再保持单调关系，因此牛顿法需调整初始点选择策略。
• 大地多边形面积计算公式为 $2\pi n - \sum E_{i-1,i}$，其中 $n$ 表示完整向东环绕圈数，确保高精度面积积分。
• 大地面积积分通过 $G(\phi, \alpha^2, k)$ 与复椭圆函数重新表述，适用于长扁情形，其中 $c^2$ 通过 $\tan^{-1}(\sqrt{-e^2})/\sqrt{-e^2}$ 定义，以保持实值表达式。