12-大地坐标系与大地极坐标系的关系2概论.ppt

大地坐标系与大地极坐标系的关系（大地问题主题解算）（2） 3、贝塞耳大地问题解算公式 将被积函数展开成σ的级数并积分可得S和σ关系式 思考讨论 思考讨论 思考讨论 U是σ和λ的函数 思路：设P1P2弧上有一动点P,它到P1 的距离σ及P的纬度u也都相应的产生变动，由此可建立u和σ的关系。 按贝塞尔公式精密解算大地问题，400～17000km的距离内，其终点大地坐标和方位角误差都小于±0.003“，距离误差小于±0.1m。 3）将球面元素换算到椭球面上 4、贝塞尔大地问题正解 Bessel’s Direct Solution 最大0.0002″ 米级近似是可仅取红色部分，百米级近似时忽略β，γ项，仅取 ② 象限判断 m、M、λ1均在第Ⅰ象限λ2与（M+σ）同象限 A1在第Ⅰ象限 A2第IV 或第III象限，取决于 tanA2的符号 ② 象限判断 A1在第Ⅲ象限 m、 M、 λ1分别在IV 、 Ⅱ和III象限，λ2与 同象限； A2第I 或第Ⅱ象限，取决于 tanA2的符号。 ② 象限判断 IV I IV IV tanA2为正在Ⅰ象限 tanA2为负在Ⅱ象限 与360°－ （M+σ）同象限 III II IV III II II I II tanA2为正在III象限 tanA2为负在IV象限 与（M+σ）同象限 I I I I A2 λ2 λ1 M m A1 ③ 计算示例 5、贝塞尔大地问题反解 Bessel’s Inverse Solution P1(L1, B1) S A1 P2(L2, B2) A2 N L1，B1，L2, B2 S, A1, A2 ① 解算步骤 1） 将椭球面元素投影到球面上 由B求u 5、贝塞尔大地问题反解 Bessel’s Inverse Solution b. 由l求λ b. 由l求λ 2）解算球面三角形 求σ 求A1、A2 ① 解算步骤 5、贝塞尔大地问题反解 Bessel’s Inverse Solution ② 象限判断 3）将球面元素换算至椭球面上 ① 解算步骤 5、贝塞尔大地问题反解 Bessel’s Inverse Solution ③ 改进算法－存在问题 a）不允许大地方位角为π/2的整数倍； b）解算过程中要多次用到繁琐的象限判定，很容易出错； 5、贝塞尔大地问题反解 Bessel’s Inverse Solution ③ 改进算法－原理 A1 90o-u1 x 90o-m M 5、贝塞尔大地问题反解 Bessel’s Inverse Solution ③ 改进算法－原理 ③ 改进算法－原理 ③ 改进算法－正解步骤 b. 将S化为σ a. 将B1规范在 [-π/2,π/2]并求相关量的正余弦 ③ 改进算法－正解步骤 c. 求B2、λ ③ 改进算法－正解步骤 d. 求A2、L2 ③ 改进算法－反解步骤 5、贝塞尔大地问题反解 Bessel’s Inverse Solution a. 将椭球面上的经差l规范至[-π,π] b. 求归化纬度的正余弦 ③ 改进算法－反解步骤 c. 迭代 ③ 改进算法－反解步骤 d. 求正反方位角 内容小结 1、贝塞尔微分方程的解 2、贝塞尔大地问题正解 3、贝塞尔大地问题反解 重点： 迭代解法； 象限判断。 4、改进算法 计算实习 1、编程实现贝塞尔大地问题正反解。 大地测量学基础程序框架.cpp subroutine.h 再见！ 对于大地问题正解可直接解算，但是对于反解则需要迭代求解。 思考讨论 思考讨论 思考讨论 * * * * * * * 测绘学院一系大地测量教研室 《大地测量学基础》（FOUNDATION OF GEODESY） The Solution of the Geodetic problem Polar Geodetic Coordinate System Concept of the Solution of the Geodetic problem Bessel’s Solution of the Geodetic Problem Bessel’s Direct Solution of the Geodetic Problem Bessel’s Inverse Solution of the Geodetic Problem 1、什么是大地极坐标，它是如何表示一点的位置的？ Preview （Some Question） 椭球面上点P(S,A) N P1 P S A 2、什么是大地问题计算？ 在椭球面上推算点的大地坐标，或者根据两点的大地坐标计算大地线长和大地方位角，这样的计算问题就叫做大地问题解算（The solution of the Geode