椭球数学投影变换(§4.7).ppt

1 大地测量主题解算 4 7 1大地主题解算的一般说明主题解算分为 短距离 400km 中距离 1000km 长距离 1000km以上 2 1 以大地线在大地坐标系中的微分方程为基础 直接在地球椭球面上进行积分运算 主要特点 解算精度与距离有关 距离越长 收敛越慢 因此只适用于较短的距离典型解法 高斯平均引数法 大地测量主题解算 3 2 以白塞尔大地投影为基础1 按椭球面上的已知值计算球面相应值 即实现椭球面向球面的过渡 2 在球面上解算大地问题 3 按球面上得到的数值计算椭球面上的相应数值 即实现从圆球向椭球的过渡 典型解法 白塞尔大地主题解算 特点 解算精度与距离长短无关 它既适用于短距离解算 也适用于长距离解算 可适应20000km或更长的距离 这对于国际联测 精密导航 远程导弹发射等都具有重要意义 大地测量主题解算 4 4 7 2勒让德级数式为了计算的级数展开式 关键问题是推求各阶导数 大地测量主题解算 5 一阶导数 二阶导数 大地测量主题解算 6 三阶导数 大地测量主题解算 7 大地测量主题解算 8 大地测量主题解算 9 大地测量主题解算 10 4 7 3高斯平均引数正算公式高斯平均引数正算公式推导的基本思想 首先把勒让德级数在P 点展开改在大地线长度中点M展开 以使级数公式项数减少 收敛快 精度高 其次 考虑到求定中点M的复杂性 将M点用大地线两端点平均纬度及平均方位角相对应的m点来代替 并借助迭代计算便可顺利地实现大地主题正解 大地测量主题解算 11 1 建立级数展开式 大地测量主题解算 12 同理可得 2 大地测量主题解算 13 大地测量主题解算 14 大地测量主题解算 3 由大地线微分方程依次求偏导数 15 大地测量主题解算 16 大地测量主题解算 17 同理可得 大地测量主题解算 18 注意 从公式可知 欲求 及 必先有 及 但由于 2和 21未知 故精确值尚不知 为此须用逐次趋近的迭代方法进行公式的计算 除此之外 此方法适合与200公里以下的大地问题解算 其计算经纬计算精度可达到0 0001 方位角计算精度可达到0 001 19 4 7 4高斯平均引数反算公式高斯平均引数反算公式可以依正算公式导出 上述两式的主式为 20 21 已知 求得 22 4 7 5白塞尔大地主题解算方法 白塞尔法解算大地主题的基本思想 以辅助球面为基础 将椭球面三角形转换为辅助球面的相应三角形 由三角形对应元素关系 将椭球面上的大地元素按照白塞尔投影条件投影到辅助球面上 然后在球面上进行大地主题解算 最后再将球面上的计算结果换算到椭球面上 这种方法的关键问题是找出椭球面上的大地元素与球面上相应元素之间的关系式 同时也要解决在球面上进行大地主题解算的方法 23 在球面上进行大地主题解算球面上大地主题正算 已知求解球面上大地主题反算 已知求解 24 1 球面三角元素间的相互关系 25 球面上大地主题正解 26 球面上大地主题反解方法 27 2 椭球面和球面上坐标关系式 28 在椭球面上与单位球面上的大地线微分方程为 29 白塞尔提出如下三个投影条件 1 椭球面大地线投影到球面上为大圆弧 2 大地线和大圆弧上相应点的方位角相等 3 球面上任意一点纬度等于椭球面上相应点的归化纬度 30 31 以上为白塞尔微分方程 32 3 白塞尔微分方程的积分 33 34 积分得到下式 35 反算 正算 迭代法 直接法 36 适合于反算 适合于正算 迭代法 直接法 37 38 将三角函数幂级数用倍角函数代替 合并同类项 积分 截去4倍角项 其值小于0 0001秒 39 正算 反算 40 4白塞尔法大地主题正算步骤 1 计算起点的归化纬度2 计算辅助函数值 解球面三角形可得 3 按公式计算相关系数A B C以及 41 4 计算球面长度迭代法 直接法 42 5 计算经度差改正数6 计算终点大地坐标及大地方位角 43 44 5白塞尔法大地主题反算步骤 1 辅