大地测量学资源课-第4.7节 大地测量主题解算.ppt

4.7.5 白塞尔大地主题解算方法 白塞尔法解算大地主题的基本思想: 以辅助球面为基础,将椭球面三角形转换为辅助球面的相应三角形,由三角形对应元素关系,将椭球面上的大地元素按照白塞尔投影条件投影到辅助球面上，然后在球面上进行大地主题解算，最后再将球面上的计算结果换算到椭球面上 关键问题是找出椭球面上的大地元素与球面上相应元素之间的关系式,同时也要解决在球面上进行大地主题解算的方法。 白塞尔大地主题解算 1、球面上进行大地主题解算 球面上大地主题正算: 已知： 求解： 球面上大地主题反算: 已知： 求解： 白塞尔大地主题解算 球面上大地主题解算公式 球面三角元素间的相互关系（球面三角公式） 正弦定理： 余弦定理： 五元素公式： 球面上大地主题正解 球面上大地主题解算公式 已知： ，求 球面上大地主题反解方法 球面上大地主题解算公式 已知： ，求 其中： 2 、椭球面和球面上坐标关系式 白塞尔大地主题解算方法 椭球面上与单位球面上的大地线微分方程为: 白塞尔大地主题解算方法 1.椭球面大地线投影到球面上为大圆弧 2.大地方位角A12投影后数值不变； 白塞尔大地主题解算方法 3.球面上任意一点纬度等于椭球面上相应点的归化纬度。 白塞尔提出如下三个投影条件： 按照上述条件，在球面极三角形PP1P2中，依正弦定理可得 白塞尔大地主题解算方法 依据大地线克莱劳方程可得 由上两式可得 至此白塞尔六个元素中，建立了4个元素之间的关系： 白塞尔大地主题解算方法 因为： 以上两式为白塞尔微分方程. 白塞尔大地主题解算方法 同理根据椭球面与球面大地微分方程的关系式可知： 3、白塞尔微分方程的积分 白塞尔大地主题解算方法 在球面直角三角形PP0Q1中： 白塞尔大地主题解算方法 令 因为上式的最后一项积分， 可忽略不计。 白塞尔大地主题解算方法 同类项合后并得到 在保持解算具有足够精度的情况下，其它函数积分后得： 积分得到下式： 白塞尔大地主题解算方法 因为 其中各系数项为： 适合于反算: 白塞尔大地主题解算方法 适合于正算: 迭代法: 直接法: 白塞尔大地主题解算方法 研究 L 的积分问题： 将被积函数按级数展开可得： 将三角函数幂级数用倍角函数代替，合并同类项，积分。截去4倍角项，其值小于0.0001秒。 对上式进行积分，则可得到经差计算公式 白塞尔大地主题解算方法 正算： 白塞尔大地主题解算方法 反算： 白塞尔法大地主题正算步骤 白塞尔法大地主题正算步骤 1.计算起点的归化纬度 2.计算辅助函数值，解球面三角形得: 3. 按公式计算相关系数 A, B, C以及α,β 白塞尔法大地主题正算步骤 4.计算球面长度 迭代法: 直接法: 循环计算满足: 白塞尔法大地主题正算步骤 5.计算经度差改正数 6.计算终点大地坐标及大地方位角 白塞尔法大地主题正算步骤 注意：上述相关元素的计算要进行象限的判断与确定，见教课书相关内容。 白塞尔法大地主题反算步骤 1.辅助计算 白塞尔法大地主题反算步骤 2. 逐次趋近法同时计算起点大地方位角、球面长度及经差，第一次趋近时，取δ＝０。 白塞尔法大地主题反算步骤 重复上述计算过程2，两次迭代满足 要求为止 白塞尔法大地主题反算步骤 3. 计算大地线长度S 4. 计算反方位角 白塞尔法大地主题反算应注意的问题 当两点在赤道时，由公式可知归化纬度 u1=0,u2 =0，方位角A12无法计算,导致后续计算无法进行。因此在此特殊情况下，还在计算步骤②之前加入如下算法： 当两点同在子午线上时, A12, A21按象限判断也得不到正确的解 另外也可采用其它计算方法 ,即计算步骤②之前，首先判断两点是否同在子午线上，如果不在子午线上，按公式计算；如果在子午线上，大地线长度依照子午线长度公式计算，大地方位角则满足： 大地主题计算算例 白塞尔法大地主题计算算例 算例 例1（80km) 例2（400km) 例3（15000km) 参考椭球 克拉索夫斯基椭球 克拉索夫斯基椭球 克拉索夫斯基椭球 已 知 数 据 参 考 值 大地主题计算程序： 白塞尔法大地主题计算程序 教学程序演示：（略） 4.7 大地测量主题解算 4.7大地测量主题解算 椭球面上的计算问题： 椭球面上点的大地经度 L、大地纬度B，两点间的大地线长度 S 以及正反大地方位角A12、A