2019年最新-椭球面上的测量计算-精选文档.ppt

4 、平均曲率半径 由于 R A 的数值随方位 A 的变化而变化，给测量带来不便，在测量工作中， 往往根据一定的精度要求，在一定范围内，把椭球面当作球面来处理，为此， 就要推求该球面的曲率半径 -- 平均曲率半径 [ 就是过椭球面上一点的一切法截 弧 (0 —2π），当其数目趋于无穷时，它们的曲率半径的算术平均值的极限， 就称为平均曲率半径，用 R 表示 ] 。 b c Na 2 R ? M N ? 2 ? ?? ( 1 ? e ) 2 2 W V V W . M 、 N 、 R 的关系： NR M 只有在极点上，它们才相等，且均等于极曲率半 径 c ，即： N ? R ? M ? c 9 0 9 0 9 0 31 7.4 椭球面上的弧长计算 1. 子午线弧长计算公式 子午椭圆的一半，它的端点与极点相重合；而赤道 又把子午线分成对称的两部分。 P ? dx 如下图所示，取子午线上某微分弧 P ? ，令 Ｐ 点 纬度为 Ｂ ， dB ， Ｐ 点的子午圈曲率半径为 P ? 点纬度为 B ? Ｍ ，于是有： dx ? MdB 从赤道开始到任意纬度 Ｂ 的平行圈 之间的弧长可由下列积分求出： X ? ? B 0 MdB 式中 M 可用下式表达： M ? a ? a cos 2 B ? a cos 4 B ? a cos 6 B ? a cos 8 B 0 2 4 6 8 32 a 0 ? a 2 ? 其中： a 4 ? a 6 ? a 8 ? 经积分，进行整理后得子午线弧长计算式： m 2 3 5 35 ? m 0 ? ? m 4 ? m 6 ? m 8 ? ? ? 2 8 16 128 ? m 2 m 4 15 7 ? ? ? m 6 ? m 8 ? 2 2 32 16 ? m 4 3 7 ? ? m 6 ? m 8 ? 8 16 32 ? m 6 m 8 ? ? ? 32 16 ? m 8 ? ? 128 ? a a a a 6 8 2 4 X ? a B ? si 2 B n ? si 4 B n ? si 6 B n ? si 8 B n 0 2 4 6 8 为求子午线上两个纬度 B 1 及 B 1 间的弧长，只需按上式分别算 出相应的 X 1 及 X 1 ，而后取差： ? X 2 ? X ，该 ? ? X X 即为所求的 1 弧长。 33 克拉索夫斯基椭球子午线弧长计算公式： ? X ? 111134 . 861 B ? 16036 . 480 sin 2 B ? 16 . 828 sin 4 B ? 0 . 02 si 6 B 1975 年国际椭球子午线弧长计算公式： ? X ? 111133 . 005 B ? 16038 . 528 sin 2 B ? 16 . 833 sin 4 B ? 0 . 02 si 6 B 34 2. 由子午线弧长可计算大地纬度 在高斯投影反算时，已知高斯平面直角坐标（ X ， Y ）反 求其大地坐标（ L ， B ）。首先 X 当作中央子午线上弧长，反求 其纬度，此时的纬度称为底点纬度或垂直纬度。计算底点纬度 的公式可以采用迭代解法和直接解法。 迭代法 ? X ? 111134 . 861 B ? 16036 . 480 sin 2 B ? 16 . 828 sin 4 B ? 0 . 02 si 6 B 在克拉索夫斯基椭球上计算时，迭代开始时设 1 B . 8611 f ? X / 111134 以后每次迭代按下式计算： i ? 1 i B ? ( X ? F ( B / 111134 . 861 f f )) F ( B ) ? ? 16036 . 4803 sin 2 B ? 16 . 8281 sin 4 B ? 0 . 0220 sin 6 B i ? 1 i B ? B ? ? 重复迭代直至 f f 为止。 i f i f i f i f 35 在 1975 年国际椭球上计算时，也有类似公式。 2. 平行圈弧长公式 旋转椭球体的平行圈是一个圆，其半径就是圆上任意一点的 子午面直角坐标 x ： r ? x ? N c o s B ? a c o s B 1 ? e s i n B 2 2 如果平行圈上有两点，其经 ? ? l ? L L 差 2 ? 1 ， 可写出平行圈弧长公式： ? ? l ? l ? S ? r ? ? N c o