椭球数学投影变换节.pptx

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