第四章 地球椭球数学投影的基本理论(郭).ppt

《第四章 地球椭球数学投影的基本理论(郭).ppt》由会员分享，可在线阅读，更多相关《第四章 地球椭球数学投影的基本理论(郭).ppt（204页珍藏版）》请在金锄头文库上搜索。

1、1 第四章地球椭球数学投影的基本理论 2 4 1地球椭球基本参数及其互相关系地球椭球是选择的旋转椭球 旋转椭球的形状和大小常用子午椭圆的五个基本几何参数 或称元素 长半轴 短半轴 椭圆的扁率椭圆的第一偏心率椭圆的第二偏心率通常用a 3 为简化书写 还常引入以下符号 椭球基本参数及其互相关系 4 4 2椭球面上常用坐标系及其关系4 2 1各种坐标系的建立1 大地坐标系大地经度B大地纬度L大地高H 大地坐标是大地测量的基本坐标系 具有如下的优点 见书P99页 5 2 空间直角坐标系坐标原点位于总地球椭球 或参考椭球 质心 Z轴与地球平均自转轴相重合 亦即指向某一时刻的平均北极点 X轴指向平均自转轴

2、与平均格林尼治天文台所决定的子午面与赤道面的交点G Y轴与此平面垂直 且指向东为正 地心空间直角系与参心空间直角坐标系之分 常用坐标系及其关关系 6 3 子午面直角坐标系设P点的大地经度为L 在过P点的子午面上 以子午圈椭圆中心为原点 建立x y平面直角坐标系 在该坐标系中 P点的位置用L x y表示 常用坐标系及其关系 7 4 地心纬度坐标系及归化纬度坐标系设椭球面上P点的大地经度L 在此子午面上以椭圆中心O为原点建立地心纬度坐标系 以椭球长半径a为半径作辅助圆 延长 与辅助圆相交 点 则OP 与x轴夹角称为P点的归化纬度u 常用坐标系及其关系 8 常用坐标系及其关系 5 大地极坐标系M是椭

3、球面上一点 MN是过M的子午线 S为连接MP的大地线长 A为大地线在M点的方位角 以M为极点 MN为极轴 P点极坐标为 S A 9 常用坐标系及其关系 4 2 2坐标系之间的相互关系子午平面坐标系同大地坐标系的关系 10 常用坐标系及其关系 令 pn N 11 常用坐标系及其关系 空间直角坐标同子午面直角坐标系的关系 12 常用坐标系及其关系 空间直角坐标系同大地坐标系 在椭球面上的点 不在椭球面上的点 13 常用坐标系及其关系 由空间直角坐标计算相应大地坐标 14 大地纬度B 归化纬度u 地心纬度 之间的关系B和u之间的关系 常用坐标系及其关系 15 常用坐标系及其关系 U 之间的关系 之间

4、的关系 大地纬度 地心纬度 归化纬度之间的差异很小 经过计算 当B 45 时 16 4 3椭球面上的几种曲率半径为了在椭球面上进行测量计算 就必须了解椭球面上有关曲线的性质 过椭球面上任意一点可作一条垂直于椭球面的法线 包含这条法线的平面叫做法截面 法截面同椭球面的截线叫法截线 可见 要研究椭球面的数学性质 就要研究法截线的性质 而法截线的曲率半径便是一个基本内容 17 包含椭球面上一点的法线 可作无数多个法截面 相应就有无数多个法截线 椭球面上法截线的曲率半径不同于球面上的法截线 大圆弧 曲率半径都等于圆球的半径 而是除两极外 椭球面上任意一点的法截线 随着它们的方向不同 每条法截线在该点的

5、曲率半径也不相同 因此 本节首先研究两个特殊方向的法截线曲率半径 卯酉圈及子午圈的曲率半径 在此基础上再研究平均曲率半径及任意方向的法截线曲率半径公式 18 子午圈曲率半径如图所示为一段子午椭圆 其中DK为一微小弧素 长度为dS 与之相应的纬度无穷小增量是dB 若弧素dx的曲率中心为n 此时线段Dn或Kn可以认为等于子午圈曲率半径 根据求任意曲线曲率半径的公式可以写出 19 椭球面上几种曲率半径 20 椭球面上几种曲率半径 21 卯酉圈曲率半径 N 卯酉圈 过椭球面上一点的法线 可作无限个法截面 其中一个与该点子午面相垂直的法截面同椭球面相截形成的闭合的圈称为卯酉圈 麦尼尔定理 假设通过曲面上

6、一点引两条截弧 一为法截弧 一为斜截弧 且在该点上这两条截弧具有公共切线 这时斜截弧在该点处的曲率半径等于法截弧的曲率半径乘以两截弧平面夹角的余弦 椭球面上几种曲率半径 22 椭球面上几种曲率半径 23 卯酉圈曲率半径的特点 卯酉圈曲率半径恰好等于法线介于椭球面和短轴之间的长度 亦即卯酉圈的曲率中心位在椭球的旋转轴上 椭球面上几种曲率半径 24 主曲率半径的计算以上讨论的子午圈曲率半径M及卯酉圈曲率半径N 是两个互相垂直的法截弧的曲率半径 这在微分几何中统称为主曲率半径 椭球面上几种曲率半径 25 椭球面上几种曲率半径 26 椭球面上几种曲率半径 27 28 任意法截弧的曲率半径 椭球面上几种

7、曲率半径 29 任意法截弧的曲率半径的变化规律 不仅与点的纬度B有关 而且还与过该点的法截弧的方位角A有关 当 时 变为计算子午圈曲率半径的 即 当 90 时 为卯酉圈曲率半径 即 主曲率半径M及N分别是 的极小值和极大值 当A由0 90 时 之值由 当A由90 180 时 值由N 可见 值的变化是以90 为周期且与子午圈和卯酉圈对称的 椭球面上几种曲率半径 30 平均曲率半径椭球面上任意一点的平均曲率半径R等于该点子午圈曲率半径M和卯酉圈曲率半径N的几何平均值 椭球面上几种曲率半径 31 M N R的关系 椭球面上几种曲率半径 32 对于克拉索夫斯基椭球 椭球面上几种曲率半径 33 4 4椭

8、球面上的弧长计算子午线弧长计算公式 34 椭球面上的弧长计算 35 椭球面上几种曲率半径 36 如果以B 90 代入 则得子午椭圆在一个象限内的弧长约为10002137m 旋转椭球的子午圈的整个弧长约为40008549 995m 即一象限子午线弧长约为10000km 地球周长约为40000km 为求子午线上两个纬度B 及 间的弧长 只需按 11 42 式分别算出相应的X 及X 而后取差 该 即为所求的弧长 当弧长甚短 例如X 40km 计算精度到0 001m 可视子午弧为圆弧 而圆的半径为该圆弧上平均纬度点的子午圈的曲率半径M 椭球面上的弧长计算 37 由子午弧长求大地纬度迭代解法 平行圈弧长

9、公式 椭球面上的弧长计算 38 椭球面上的弧长计算 子午线弧长和平行圈弧长变化的比较 39 4 5大地线两点间的最短距离 在平面上是两点间的直线 在球面上是两点间的大圆弧 那么在椭球面上又是怎样的一条线呢 它应是大地线 相对法截线 40 相对法截线 大地线 41 相对法截线的特点 当A B两点位于同一子午圈或同一平行圈上时 正反法截线则合二为一 在通常情况下 正反法截线是不重合的 因此在椭球面上A B C三个点处所测得的角度 各点上正法截线之夹角 将不能构成闭合三角形 为了克服这个矛盾 在两点间另选一条单一的大地线代替相对法截线 从而得到由大地线构成的单一的三角形 大地线 42 大地线 大地线

10、的定义和性质 椭球面上两点间的最短程曲线叫大地线 43 大地线的性质 大地线是两点间惟一最短线 而且位于相对法截线之间 并靠近正法截线 它与正法截线间的夹角在椭球面上进行测量计算时 应当以两点间的大地线为依据 在地面上测得的方向 距离等 应当归算成相应大地线的方向 距离 长度差异可忽略 方向差异需改化 大地线 44 大地线的微分方程和克莱劳方程 大地线的微分方程 45 大地线的微分方程 46 大地线的微分方程 大地线的克莱劳方程 在旋转椭球面上 大地线各点的平行圈半径与大地线在该点的大地方位角的正弦的乘积等于常数 式中常数C也叫大地线常数 47 当大地线穿越赤道时当大地线达极小平行圈时由克莱劳

11、方程可以写出 48 4 6将地面观测值归算至椭球面观测的基准线不是各点相应的椭球面的法线 而是各点的垂线 各点的垂线与法线存在着垂线偏差 归算的两条基本要求 以椭球面的法线为基准 将地面观测元素化为椭球面上大地线的相应元素 将地面观测的水平方向归算至椭球面将水平方向归算至椭球面上 包括垂线偏差改正 标高差改正及截面差改正 习惯上称此三项改正为三差改正 49 垂线偏差改正以测站A为中心作出单位半径的辅助球 u是垂线偏差 它在子午圈和卯酉圈上的分量分别以 表示 M是地面观测目标m在球面上的投影 垂线偏差对水平方向的影响是 R R1 地面观测值归算至椭球面 50 标高差改正 地面观测值归算至椭球面

12、51 截面差改正 地面观测值归算至椭球面 52 将地面观测的长度归算至椭球面基线尺量距的归算将基线尺量取的长度加上测段倾斜改正后 可以认为它是基线平均高程面上的长度 以 表示 现要把它归算至参考椭球面上的大地线长度S 1 垂线偏差对长度归算的影响 地面观测值归算至椭球面 53 2 高程对长度归算的影响 地面观测值归算至椭球面 54 电磁波测距的归算 地面观测值归算至椭球面 55 地面观测值归算至椭球面 56 大地测量主题解算 4 7 1大地主题解算的一般说明主题解算分为 短距离 400km 中距离 1000km 长距离 1000km以上 57 1 以大地线在大地坐标系中的微分方程为基础 直接在

13、地球椭球面上进行积分运算 主要特点 解算精度与距离有关 距离越长 收敛越慢 因此只适用于较短的距离典型解法 高斯平均引数法 大地测量主题解算 58 2 以白塞尔大地投影为基础1 按椭球面上的已知值计算球面相应值 即实现椭球面向球面的过渡 2 在球面上解算大地问题 3 按球面上得到的数值计算椭球面上的相应数值 即实现从圆球向椭球的过渡 典型解法 白塞尔大地主题解算 特点 解算精度与距离长短无关 它既适用于短距离解算 也适用于长距离解算 可适应20000km或更长的距离 这对于国际联测 精密导航 远程导弹发射等都具有重要意义 大地测量主题解算 59 4 7 2勒让德级数式为了计算的级数展开式 关键

14、问题是推求各阶导数 大地测量主题解算 60 一阶导数 二阶导数 大地测量主题解算 61 三阶导数 大地测量主题解算 62 大地测量主题解算 63 大地测量主题解算 64 大地测量主题解算 65 4 7 3高斯平均引数正算公式高斯平均引数正算公式推导的基本思想 首先把勒让德级数在P 点展开改在大地线长度中点M展开 以使级数公式项数减少 收敛快 精度高 其次 考虑到求定中点M的复杂性 将M点用大地线两端点平均纬度及平均方位角相对应的m点来代替 并借助迭代计算便可顺利地实现大地主题正解 大地测量主题解算 66 1 建立级数展开式 大地测量主题解算 67 同理可得 2 大地测量主题解算 68 大地测量

15、主题解算 69 大地测量主题解算 3 由大地线微分方程依次求偏导数 70 大地测量主题解算 71 大地测量主题解算 72 同理可得 大地测量主题解算 73 注意 从公式可知 欲求 及 必先有 但由于 2和 21未知 故精确值尚不知 为此须用逐次趋近的迭代方法进行公式的计算 除此之外 此方法适合与200公里以下的大地问题解算 其计算经纬计算精度可达到0 0001 方位角计算精度可达到0 001 74 4 7 4高斯平均引数反算公式高斯平均引数反算公式可以依正算公式导出 上述两式的主式为 75 76 已知 求得 77 4 7 5白塞尔大地主题解算方法 白塞尔法解算大地主题的基本思想 以辅助球面为基

16、础 将椭球面三角形转换为辅助球面的相应三角形 由三角形对应元素关系 将椭球面上的大地元素按照白塞尔投影条件投影到辅助球面上 然后在球面上进行大地主题解算 最后再将球面上的计算结果换算到椭球面上 这种方法的关键问题是找出椭球面上的大地元素与球面上相应元素之间的关系式 同时也要解决在球面上进行大地主题解算的方法 78 在球面上进行大地主题解算球面上大地主题正算 已知求解球面上大地主题反算 已知求解 79 1 球面三角元素间的相互关系 80 球面上大地主题正解 81 球面上大地主题反解方法 82 2 椭球面和球面上坐标关系式 83 在椭球面上与单位球面上的大地线微分方程为 84 白塞尔提出如下三个投影条件 1 椭球面大地线投影到球面上为大圆弧 2 大地线和大圆弧上相应点的方位角相等 3 球面上任意一点纬度等于椭球面上相应点的归化纬度 85 86 以上为白塞尔微分方程 87 3 白塞尔微分方程的积分 88 89 积分得到下式 90 反算 正算 迭代法 直接法 91 适合于反算 适合于正算 迭代法 直接法 92 93 将三角函数幂级数用倍角函数代替 合并同类项 积分 截去4倍角项 其值小于0 0