大地问题正解与反解

从解析意义上看，大地问题的正解与反解，就是研究大地坐标与大地极坐标 相互变换的问题。

按大地线的长度，大地主题解算可分为短距离（小于 400km）、中距离（400～ 1000km）和长距离（大于 1000km）。短距离的大地主题解算，大多应用于传统 的国家一等三角测量中。中、长距离大地问题解算主要应用于洲际联测、导弹与 火箭发射、导航和宇宙航行等。

椭球面上的大地问题解算比平面坐标计算复杂的多，多年来大量数学家和测 量学家致力于大地问题解算的研究，并提出了 70 多种解算方法和公式。根据这 些不同解法的理论基础，可以归纳为以下四种：

（1）以大地线在大地坐标系中的微分方程为基础 直接在椭球面上进行积分运算。大地线微分方程

（1-1）式把大地线的长度S作为独立变量，将B，L，A，S联系在了一起。

沿P1点和P2点的大地线S积分，得：

 A2  − A1  ± 180︒ = P1 N dS

但是上述积分在初等函数中无法计算，难以求得精确值。所以需要进行趋近

解算，将上式积分进行变换。方法之一是运用勒让德级数，将其展开成大地线长

度S的升幂级数，再逐项计算以完成主题解算。高斯平均引数公式是这种解法的 代表。但是其解算精度与距离相关，距离越长，收敛越慢，所以只适用于较短的

距离。

(2) 以白塞尔大地投影为基础 地球椭球的形状与圆球相似，所以在球面上进行大地主题解算可利用球面三

角学函数。将椭球面上的大地线投影到球面上为大圆弧，大地线上的点与大圆弧

上相应的点一致，就实现了大地投影。如果解决了大地线上某点的数值B、L、A、S 和球面大圆弧对应的数值ϕ、λ、α、σ的关系，下式的微分方程就能实现：

由此可以得到大地主题解算的步骤：

①由椭球面上已知值计算球面上相应值，即实现从椭球过渡到圆球；

②在球面上解算大地问题；

③由球面上得到的值计算椭球面上的相应值，即实现从圆球过渡到椭球。

计算公式展开e2或e︐2的幂级数，解算精度与距离的长短无关。因此，解算 短距离和长距离问题都能适用。依据白塞尔的解法可以派生出各种公式，有的是

直接解法，有的是逐渐趋近的解法。这些公式大多可以适应 2 万公里或更长的距 离，这对于国际联测，精密导航远程导弹发射，导航和宇宙航行等具有重要意义。

（3）对大地线微分方程进行数值积分的解法 这种解法不采用辅助面或运用勒让德级数，而是直接对（1-2）式进行数值论文网

积分计算。常用的方法有高斯法，牛顿法，龙格—库塔法以及契巴雪夫法。这种 解法便于程序的编写，而且适用于任意距离。但是，这种方法随着距离的增长， 计算量增大且精度降低，在近极地区几乎无法使用。

（4）依据大地线外的其他线为基础 除大地线外，还有其他一些线连接椭球面两点，比如弦线、法截线等。这种

大地主题解算实质是三维大地测量问题。运用该方法进行三边测量的大地主题解 算有其优点。但是，解算结果还应加上归化至大地线的改正，比较麻烦，所以应 用得较少。

2 白塞尔大地主题解算公式推导

白塞尔大地主题解算的基本思想是：先将椭球面上的大地主题元素按照白塞 尔投影条件投影到辅助球面上，再在辅助球面上进行解算，最后再将球面上的计 算结果换算到椭球面上。找出椭球面上的大地元素与球面上相应元素的关系是这 种方法的关键；同时还要在球面上进行大地主题解算。 基于CGCS2000椭球的白塞尔大地主题解算公式推导及程序设计(2):http://www.youerw.com/shuxue/lunwen_80262.html