第四章椭球数学投影变换(7节)教学幻灯片

大地测量主题解算4.7.1大地主题解算的一般说明

主题解算分为:短距离(<400km)中距离(<1000km)长距离(1000km以上)

11.以大地线在大地坐标系中的微分方程为基础，直接在地球椭球面上进行积分运算。主要特点：解算精度与距离有关，距离越长，收敛越慢，因此只适用于较短的距离典型解法：高斯平均引数法

大地测量主题解算22.以白塞尔大地投影为基础1)按椭球面上的已知值计算球面相应值，即实现椭球面向球面的过渡；2)在球面上解算大地问题；3)按球面上得到的数值计算椭球面上的相应数值，即实现从圆球向椭球的过渡。典型解法：白塞尔大地主题解算

特点：解算精度与距离长短无关，它既适用于短距离解算，也适用于长距离解算。可适应20000km或更长的距离，这对于国际联测，精密导航，远程导弹发射等都具有重要意义。

大地测量主题解算3一阶导数：二阶导数：

大地测量主题解算5三阶导数

大地测量主题解算6

大地测量主题解算7

大地测量主题解算8

大地测量主题解算94.7.3高斯平均引数正算公式

高斯平均引数正算公式推导的基本思想：

首先把勒让德级数在P１点展开改在大地线长度中点M展开，以使级数公式项数减少，收敛快，精度高；其次，考虑到求定中点M的复杂性，将M点用大地线两端点平均纬度及平均方位角相对应的m点来代替，并借助迭代计算便可顺利地实现大地主题正解。

大地测量主题解算10(1)建立级数展开式:

大地测量主题解算11同理可得:

(2)

大地测量主题解算12

大地测量主题解算13

大地测量主题解算(3)由大地线微分方程依次求偏导数:14

大地测量主题解算15大地测量主题解算16同理可得：大地测量主题解算17

注意：从公式可知，欲求ΔＬ，ΔＢ及ΔＡ，必先有Ｂｍ及Ａｍ。但由于Ｂ2和Ａ21未知，故精确值尚不知，为此须用逐次趋近的迭代方法进行公式的计算。除此之外，此方法适合与200公里以下的大地问题解算，其计算经纬计算精度可达到0.0001”,方位角计算精度可达到0.001”。184.7.4高斯平均引数反算公式

高斯平均引数反算公式可以依正算公式导出:上述两式的主式为:1920已知：求得：214.7.5白塞尔大地主题解算方法

白塞尔法解算大地主题的基本思想:

以辅助球面为基础,将椭球面三角形转换为辅助球面的相应三角形,由三角形对应元素关系,将椭球面上的大地元素按照白塞尔投影条件投影到辅助球面上，然后在球面上进行大地主题解算，最后再将球面上的计算结果换算到椭球面上。

这种方法的关键问题是找出椭球面上的大地元素与球面上相应元素之间的关系式,同时也要解决在球面上进行大地主题解算的方法。

22在球面上进行大地主题解算

球面上大地主题正算:已知求解

球面上大地主题反算:已知

求解231、球面三角元素间的相互关系24

球面上大地主题正解25

球面上大地主题反解方法

262

、椭球面和球面上坐标关系式27在椭球面上与单位球面上的大地线微分方程为:28白塞尔提出如下三个投影条件：

1.椭球面大地线投影到球面上为大圆弧2.大地线和大圆弧上相应点的方位角相等；3.球面上任意一点纬度等于椭球面上相应点的归化纬度。

2930以上为白塞尔微分方程.313、白塞尔微分方程的积分3233积分得到下式：34适合于反算:适合于正算:

迭代法:

直接法:3536将三角函数幂级数用倍角函数代替，合并同类项，积分。截去4倍角项，其值小于0.0001秒。其纬差计算公式37正算：反算：384白塞尔法大地主题正算步骤

1.计算起点的归化纬度2.计算辅助函数值，解球面三角形可得:3.按公式计算相关系数A,B,C以及α,β

394.计算球面长度迭代法:

直接法:405.计算经度差改正数6.计算终点大地坐标及大地方位角

41注意：上述相关元素的计算要进行象限的判断与确定。425白塞尔法大地主题反算步骤

1.辅助计算432.用逐次趋近法同时计算起点大地方