第二章 地球坐标系和地球椭球.

Presentation on theme: "第二章 地球坐标系和地球椭球."— Presentation transcript:

1 第二章 地球坐标系和地球椭球

2 §2.1 概 述 大地测量采用的坐标系：天球坐标系、地球坐标系 地球坐标系：固定在地球上与地球一起自转和公转的 坐标系
§2.1 概 述 大地测量采用的坐标系：天球坐标系、地球坐标系 地球坐标系：固定在地球上与地球一起自转和公转的 坐标系 地球坐标系分类：参心坐标系、地心坐标系 定义坐标系的要素：原点位置、尺度与坐标轴指向；还包括一些天文、物理、地球等参数，若采用大地 坐标表述形式还需要椭球元素。

3 §2.2 地球椭球面的数学计算和有关计算 2.2.1 地球椭球的几何、物理元素 椭球方程： Z 扁率： 第一偏心率： 第二偏心率： O Y
X Y Z O

4 2.2.1 地球椭球的几何、物理元素（续1） 几个关系式： 1954年北京坐标系，克拉索夫斯基椭球元素：

5 2.2.1 地球椭球的几何、物理元素（续2） 1980年大地坐标系采用第16届 IAG—IUGG 椭球，其椭球元素为：

6 2.2.2 旋转椭球面的参数表示及数学性质 1、经线和纬线的曲线方程 在XOZ坐标面上的起始经线方程：
Y Z M1 M0 M L r A R S 1、经线和纬线的曲线方程 在XOZ坐标面上的起始经线方程： M0饶Z轴旋转，形成纬圈（平行圈），其半径： 经度为L的经线方程：

7 2.2.2 旋转椭球面的参数表示及数学性质（续1） O X Y Z M1 M0 M L r A R S 纬圈方程：

8 2.2.2 旋转椭球面的参数表示及数学性质（续2） 2、椭球面法线与子午线主法线的同一性、经纬线的Frenet标架 如图为过M点的子午面。
P O Q M P´ R T N A 如图为过M点的子午面。 子午线的主法线MP´位于 子午面内，且垂直于子午 线切线T；R为过M点的 平行圈切线，显然R垂直 于M点的子午面，因此R垂 直于MP´。所以， MP´垂直 于椭球面在M点的切平面，因此 它是椭球面的法线。 Frenet标架：曲线上任意一点处的三个相互正交的单位向量〔取切向、主法向和与该两个方向正交的第三个方向〕构成的三维直角坐标系。

9 2.2.2 旋转椭球面的参数表示及数学性质（续3） 3、旋转椭球面及经纬线的参数方程 1). 以大地经度L及归化纬度u为参数的方程
在XOZ子午面内，有 u a X Z O M´ M 在三维空间坐标系中：

10 2.2.2 旋转椭球面的参数表示及数学性质（续4） (2). 以大地经纬度L、B为参数的方程 切线M0T的斜率的导数式： 由椭圆方程求导得：
X Z K0 B 90°+ B O T M 0 切线M0T的斜率的导数式： 由椭圆方程求导得： 代入第一式得： 1

11 2.2.2 旋转椭球面的参数表示及数学性质（续5） 将 代入椭圆方程，化简后得： 1 引入辅助符号： 则有：

12 2.2.2 旋转椭球面的参数表示及数学性质（续6） 以大地经纬度两个参数表示的椭球面上一点的三维坐标为〔椭球面参数方程式〕：

13 以大地纬度为参数的经度为LC的子午线参数方程为:

14 在一点〔BC ,LC 〕处的子午线切向量

15 子午线切线单位向量

16 以大地经度为参数的大地纬度为 的BC纬线的参数方程为

17 在一点〔BC ,LC 〕处的平行圈切向量

18 平行圈切线单位向量

19 椭球面单位法向量为其矢量积：

20 2.2.2 旋转椭球面的参数表示及数学性质（续7） (3). 以大地经度L及球心纬度为参数的方程 球心纬度，向径，则对于
XOZ平面上的椭圆有： X Z  O M 0  在椭圆上，向径由球 心纬度唯一确定，将上式 代入椭圆方程，得：

21 2.2.2 旋转椭球面的参数表示及数学性质（续8） 对于XOZ平面上的椭圆有：
以经度、球心纬度两个参数表示的椭球面上一点的三维坐标为〔参数方程式〕为：

22 2.2.2 旋转椭球面的参数表示及数学性质（续9） 不难得出，u, B,  的关系为： 因此有： 由球心纬度公式，得：

23 4、旋转椭球面的几何性质 a). 对称性 b). 有界性 c). 正则性：曲面上每点都对应于唯一 确定的非零法向量。 d). 不可展性
2.2.2 旋转椭球面的参数表示及数学性质（续10） 4、旋转椭球面的几何性质 a). 对称性 b). 有界性 c). 正则性：曲面上每点都对应于唯一 确定的非零法向量。 d). 不可展性

24 2.2.3 法截线曲率及曲率半径 1、空间曲线的曲率几曲率半径 若以曲线的弧长s为参数，曲线上的点位用向量r(s)表示。则曲线的曲率为：
法截线曲率及曲率半径 1、空间曲线的曲率几曲率半径 若以曲线的弧长s为参数，曲线上的点位用向量r(s)表示。则曲线的曲率为： 若以t参数，则曲线的曲率可表示为：

25 2.2.3 法截线曲率及曲率半径（续1） 2、椭球面法截线的曲率 (1). 子午线曲率半径
法截线曲率及曲率半径（续1） 2、椭球面法截线的曲率 (1). 子午线曲率半径 不失一般性，以起始子午线为例推导。若以归化纬度u为子午线方程的参数，则有：

26 法截线曲率及曲率半径（续2） 则有： 同理，若以大地纬度为参数，得： 子午曲率半径M，就是曲率是倒数，即：

27 2.2.3 法截线曲率及曲率半径（续3） (2). 卯酉线曲率半径 定义：与子午面切线正交的法截面与椭球面的交线为卯 酉线。
法截线曲率及曲率半径（续3） (2). 卯酉线曲率半径 定义：与子午面切线正交的法截面与椭球面的交线为卯 酉线。 根据微分几何中的麦尼尔定理，卯酉圈曲率kn与平行圈曲率kr的关系为： 平行圈半径为子午面XOZ 平面内的X坐标，即： 则有，上述两式得卯酉曲率半径N为：

28 2.2.3 法截线曲率及曲率半径（续 4） (3). 任意方向法截线的曲率半径
法截线曲率及曲率半径（续 4） (3). 任意方向法截线的曲率半径 根据微分几何中的Euler公式，任意方向法截线的曲率与子午、卯酉曲率半径的关系为： 因此，任意方向的曲率半径为： 当A为0，/2，， 3/2时，取得极值。

29 法截线曲率及曲率半径（续 5） (4). 平均曲率半径 定义：所有方向法截线曲率半径的平均值。 代入上式，得：

30 法截线曲率及曲率半径（续 6） 不难得到：N  R  M 引入辅助量： 存在下列关系：

31 2.2.4 椭球面上第一基本形式及弧长面积计算 1. 椭球面的第一基本形式 椭球面上点的向量： 椭球面上的微分弧长： 其中： 对于椭球面：

32 2.2.4 椭球面上第一基本形式及弧长面积计算（续1）
2、子午线弧长 子午线微分弧长： 积分得： 用二项式展开，并逐项积分得： 常数 A、B、C、D、E、F、G的计算公式见教材

33 2.2.4 椭球面上第一基本形式及弧长面积计算（续2）
对于小于400km的弧长，可采用以下简化式。 其中： 根据： 求出导数，代入上式并化简，得： 对于小于40km的弧长，可进一步简化为：

34 2.2.4 椭球面上第一基本形式及弧长面积计算（续3）
已知B1和弧长S1~2求B2称为反算，可采用叠代法计算。 初值： 叠代格式： 要求： 其中：

35 2.2.4 椭球面上第一基本形式及弧长面积计算（续4）
3、平行圈的半径与弧长 相同经差的平行圈弧长在赤道最长，越靠近两极越小。

36 2.2.4 椭球面上第一基本形式及弧长面积计算（续5）
4、利用经纬格网计算椭球面的面积 L L+dL B B+dB MdB NcosBdL d

37 2.2.4 椭球面上第一基本形式及弧长面积计算（续6）
上式利用二项式展开并积分，得： 取 L2-L1 = 2，B2 = /2，B1 = 0 算得半球面积，乘2可以估算全球面积约为5.1亿平方公里

38 习 题 1、导出三种纬度 φ、u与B的关系。 2、导出子午曲率半径M与卯酉曲率半径N的计算公式。
习 题 1、导出三种纬度 φ、u与B的关系。 2、导出子午曲率半径M与卯酉曲率半径N的计算公式。 3、M、N、R的关系如何？在什么条件下三者相同？ 4、某点到赤道的子午弧长 ,求该点的纬度。 a= , =1/298.3 5、已知某点的纬度 ,求该点自赤道起的子午弧长。 a = , =1/298.3

39 2.2.5 大地线 1、大地线的定义与性质 法截弧：由椭球面上A点的法线与B点所确定的法截面与椭球面相割得到的曲线称为A到B的法截弧。
大地线 1、大地线的定义与性质 法截弧：由椭球面上A点的法线与B点所确定的法截面与椭球面相割得到的曲线称为A到B的法截弧。 相对法截弧： A到B的法截弧与B到A的法截弧。 由相对法截弧构成的椭球面三角形 不是闭合图形。

40 2、曲面上连接任何两点的最短直线必为 2.2.5 大地线（续1） 大地线的定义：大地线的主法线与曲面法线处处重合。
大地线（续1） 大地线的定义：大地线的主法线与曲面法线处处重合。 大地线的性质：1、大地线上任何点的密切平面就是该点 的法截面； 2、曲面上连接任何两点的最短直线必为 大地线。 3、大地线的测地曲率等于0 曲线的测地曲率：曲线的曲率在曲面切平面上的投影。 大地线的曲率： 大地线的挠率

41 2.2.5 大地线（续2） 2、大地坐标系中大地线的微分方程 (1). 大地线的二阶微分方程
大地线（续2） 2、大地坐标系中大地线的微分方程 (1). 大地线的二阶微分方程 以u,v 为参数的一般曲面的大地线微分方程可表示为： 下标为相应的偏导数。

42 大地线（续3） 对于椭球面，有： 代入前面公式，得： 则旋转椭球面上大地线的微分方程为：

43 2.2.5 大地线（续4） (2). 克莱劳定理 直角坐标系中的椭球面方程： 椭球面法向量为： 以大地线弧长为参数的大地线主法线向量为：
大地线（续4） (2). 克莱劳定理 直角坐标系中的椭球面方程： 椭球面法向量为： 以大地线弧长为参数的大地线主法线向量为： 两者指向一致，即：

44 大地线（续5） 由上式的前两个方程得： 1 将三维空间坐标与大地坐标的关系式及微分式代入： 代入 式，整理得： 1 2

45 大地线（续6） 将关系式： 即：大地线上各点的平行圈 半径与该点的大地线方位角 正弦的乘积是常数。 代入上式，即得克莱劳定理：

46 大地线（续7） (3). 大地线的一阶微分关系式 由克莱劳定理，微分得：

47 大地线（续8） 又如图所示： 代入上式，得： 三个微分关系式可整理为： 3

48 2.2.5 大地线（续9） 3、以弧长和大地方位角为参数的大地线方程
大地线（续9） 3、以弧长和大地方位角为参数的大地线方程 x´ y´ z´ 大地线始点坐标P0(B0,L0)，大地线上任何点的位置向量都可以展开成S,A的级数形式： 4 Frenet标架的坐标轴定义：x´指向大地线的切向t， y´指向大地线的主法向n，向内为正， z´指向大地线的副法向b，构成左手系。

49 大地线（续10） 显然有： 根据曲线论中的Frenet公式： 由以上两式可求出各阶导数：

50 大地线（续11） 将上式代入大地线展开式 ，得Frenet标架下的三维坐标： 4 5 顾及公式：

51 大地线（续12） 和： 求导得：

52 大地线（续13） 代入Frenet标架下的三维坐标公式 ，得： 5

53 2.2.5 大地线（续14） 将坐标系饶 y´逆时针旋转A， 得x”、y”、z”坐标系，则有： x" x´ z´ z" y´
大地线（续14） x´ y´ z´ z" x" 将坐标系饶 y´逆时针旋转A， 得x”、y”、z”坐标系，则有： 以P0点为原点的地平坐标系(站心坐标系) x、y、z，与x”、y”、z”坐标系的关系为： x y" z z" x" y

54 大地线（续15） 最后得到地平坐标系（站心系）中的大地线方程，称为Weingarten级数式。 6

55 大地线（续16） 法截弧为平面曲线，其挠率为0，同理可推得地平坐标系中的计算式为：

56 2.2.5 大地线（续17） 4、基于大地线的椭球面曲线坐标系 (1). 大地线极坐标系 大地圆：到极点具有相同大地线长
大地线（续17） 4、基于大地线的椭球面曲线坐标系 (1). 大地线极坐标系 大地圆：到极点具有相同大地线长 度的点所构成的轨迹。 由大地线长度和大地方位角可描述 曲面点的位置 。 如图所示： 对照第一基本形式，得： 由图中的微分直角三角形，得大地极坐标系中的微分关系式：

57 大地线（续18） 大地线的归化长度 m 的计算公式： 由 式求出偏导数代入得： 6

58 (1)以长度量表示的椭球面上坐标系的由来 2.2.7以长度量为坐标参数的新大地坐标系
早在1810年Soldner就提出了球面直角坐标系统 。此后Helmert，Grossmann，Heck等德国测量学者基于Soldner球面直角坐标系推广提出了椭球面直角坐标系。教材上所述的测地坐标系与其有类似之处。

59 (2)一种新型的大地坐标系 2005年提出了椭球面上一种新型的大地坐标系：它仍以经纬线作为坐标曲线，且与大地坐标系之间能进行精确的坐标转换；它所采用的坐标参数是以长度而不是以角度为单位；可简化椭球面上的繁复计算。 施一民，朱紫阳，范业明.坐标参数为长度量的一种新型的大地坐标系. 同济大学学报，2005,33（11）:

60 新型大地坐标系的定义 为构建新型的大地坐标系,可在区域中心附近选择一点作为其坐标原点,其在大地坐标系中的大地经纬度设为（B0，L0）。两族互为正交的经纬线构成坐标系的坐标格网。设过经纬度为（B，L）的任一点的经线与起始纬线的交点为PB0,点P至点PB0的经线上弧长作为纵坐标sL,点P0至点PB0的纬线上弧长取为横坐标sB,如图所示。

61 新型大地坐标系的定义（续）

62 新型大地坐标系的定义（续） 过点P的纬线上的微分弧长ds′与起始纬线上相应的微分弧长ds之比n称为按纬度变化方向的长度归化因子
n= ds′/ds 在大地坐标系中，与这两个平行圈上微分弧长相应的经差dl 相等，故可精确求得 n= Nsin B/ N0sin B0 式中N0，N分别为纬度B0，B处的卯酉线曲率半径。 基于微分几何和椭球大地测量的理论，n可用新大地坐标表示为(取至二次项)

63 4.新型大地坐标与大地坐标的转换关系

64 习 题 1. 纬度相同的两个点的相对法截弧是否重合？此线是否就是大地线？ 2. 推导大地线的三个微分式。
习 题 1. 纬度相同的两个点的相对法截弧是否重合？此线是否就是大地线？ 2. 推导大地线的三个微分式。 3. 试述测地坐标系的定义？测地平行线是否等距？测地大地线是否等距？ 4. 简述weingarten级数的推导步骤。

65 §2.3 椭球面上大地坐标的计算 2.3.1 水平方向、边长观测值归算到椭球面 1、水平方向观测值归算到参考椭球面的改正
包括三项改正，称为三差改正。 (1). 垂线偏差改正 (2). 标高差改正 用椭球半径的近似值代入得：

66 2.3.1 水平方向、边长观测值归算到椭球面 (3). 法截弧方向归算到大地线方向的改正
该项改正很小，100公里约0.03“，只有一等控制网才估计此项改正。

67 2.3.1 水平方向、边长观测值归算到椭球面 2、空间边长归算至参考椭球面的改正 测线端点的大地高为： 椭球面上弦长 d 的计算公式
省略H/R的二次项，得：

68 2.3.1 水平方向、边长观测值归算到椭球面 椭球面上的弧长为：

69 2.3.1 水平方向、边长观测值归算到椭球面 3. 工程控制网中的地面观测元素的归算
以平均高程面作投影面，范围小，可以用球代替椭球；球半径采用高斯平均曲率半径。计算公式为： 不难证明：椭球半径的误差对边长归算结果影响很小，R取6371km即可，但高差误差对边长归算比较敏感。

70 2.3.2 椭球面上三角形解算 1、球面角超 三块面积之和为： 代入球面角超定义式，得：

71 2.3.2 椭球面上三角形解算 按球面三角公式： 当边长小于40公里时，第二项影响小于0.0004“，可略去

72 2.3.2 椭球面上三角形解算 2、解算球面三角形的勒让德定理
勒让德定理：对于较小的球面三角形，可用平面三角公式来解算，只需使三个平面角等于相应的球面角减去三分之一的球面角超，而边长保持不变。 A B C a b c

73 2.3.3 大地主题解算 大地主题解算分类： 正算：已知(B1, L1)，A12，S12，计算(B2, L2)，A21
反算：已知(B1, L1)， (B2, L2)， 计算A12，S12 ，A21 短距离 中距离 长距离 中短距离解算方法：按级数展开采用Gauss平均引数公式； 长距离解算方法：贝塞尔公式

74 2.3.3 大地主题解算 1、纬度差、经度差和方位角差展开为大地线长度的级数式 由大地线的微分公式，得其一阶导数为：

75 2.3.3 大地主题解算 二阶和三阶导数采用复合函数求导法计算： 同理可求出四阶以上的导数和L、A的高阶导数，代入展开式即可。

76 2.3.3 大地主题解算 2、高斯平均引数公式 若取大地线中点展开，得： 两式相减，得： 1 类似地，有：

77 2.3.3 大地主题解算 两式相加，得： 类似地，有： 其中： 将 展开成级数，得： 2

78 2.3.3 大地主题解算 由大地线的微分公式： 求导，得： 可取： 代入 式，得 的计算公式。并取 2 代入 式，求出各阶导数后整理得： 1

79 2.3.3 大地主题解算 同理可得： 3 以上3式具有4次方精度，可用于解算200公里下的大地主题。

80 2.3.3 大地主题解算 因计算Bm , Lm要用到B2 , L2，因此需要叠代计算。其初值为： 叠代计算公式为： 直到 为止。
直到 为止。 最后计算纬度、经度和方位角：

81 2.3.3 大地主题解算 3、高斯平均引数反算公式 由正算公式，反解得： 右端第二项与第一项相比为小量，可以作近似：

82 2.3.3 大地主题解算 代入上式第二项，得： 由此可求得平均方位角和大地线长度如下：

83 2.3.3 大地主题解算 由正算公式的第三式，计算a： 最后得起终点的大地方位角为：

84 (2). 由（sL ,sB）求解（B,L）:B由下式迭代反解
2.3.3 新大地主题解算 4、新大地坐标系与大地坐标系间的坐标转换 (1). 由 （B,L）求解（sL ,sB） (2). 由（sL ,sB）求解（B,L）:B由下式迭代反解 而L由下式直接反解.

85 2.3.4 大地主题微分公式 1、大地主题正解微分公式 终点的经纬度（B2,L2）和大地线方位角A21，与起点的经纬度（B1,L1）和大地线方位角A12，以及大地线长度S的微分关系。

86 2.3.4 大地主题微分公式 2、大地主题反解微分公式 起点大地线方位角A12和大地线方位角A21，以及大地线长度S与起点和终点的经纬度（B1,L1）和（B2,L2）的微分关系。

87 习 题 1、地面观测方向归算到椭球面上需要加哪几项改正？ 2、地面观测距离归算到椭球面上二步改正的几何意义？
习 题 1、地面观测方向归算到椭球面上需要加哪几项改正？ 2、地面观测距离归算到椭球面上二步改正的几何意义？ 3、 P1与P2与为控制点，已知： 计算归算到椭球面上的长度 4、已知 利用Gauss平均引数公式正反算。

88 §2.4 空间大地直角坐标系及其转换模型 2.4.1 空间直角坐标系与相应大地坐标系的关系 1、X、Y、Z与B、L、H间的关系
§2.4 空间大地直角坐标系及其转换模型 2.4.1 空间直角坐标系与相应大地坐标系的关系 1、X、Y、Z与B、L、H间的关系 空间坐标系的定义：Z自转轴，X位于赤道面，指格林尼治天文台，Y指东，构成右手系。 大地坐标的定义：B为过一点的椭球面的法线与赤道面交角、L为过同一点的子午面与起始子午面二面角的平面角，H为点沿法线到椭球面的距离。 大地高与正高、正常高之间的关系：

89 2.4.1 空间直角坐标系与相应大地坐标系的关系 如图所示： 1

90 2.4.1 空间直角坐标系与相应大地坐标系的关系 2、由X、Y、Z计算B、L、H的迭代解法 计算L： 迭代计算B： 迭代初值为：

91 2.4.1 空间直角坐标系与相应大地坐标系的关系 3、X、Y、Z与B、L、H间的微分关系 由前面 式微分得； 1 其中：对角阵

92 2.4.1 空间直角坐标系与相应大地坐标系的关系 顾及A是正交阵，J是对角阵，得：

93 2.4.1 空间直角坐标系与相应大地坐标系的关系 4、B、L、H与椭球元素a, e2 之间的微分关系
若顾及椭球元素的变化，则前面的微分公式变为： 其中：

94 2.4.1 空间直角坐标系与相应大地坐标系的关系 由上式可得： 若空间坐标系的原点和坐标轴指向保持不变，即椭球的定位与定向不变，则：

95 2.4.1 空间直角坐标系与相应大地坐标系的关系 上式简化成大地坐标与椭球元素间的微分关系：

96 2.4.2 空间直角坐标系之间的旋转变换 两个右手旋转坐标系之间的旋转角，取逆时针旋转为正，顺时针旋转为负，旋转矩阵为正交阵，可表示为：

97 2.4.2 空间直角坐标系之间的旋转变换 方法一： 将X’、Y’、Z’转换到X、Y、Z坐标系：
先绕Z’将X’旋转到XOY平面与X’OY’平面的交线X” ，再绕X” 轴将Z’旋转到Z轴，最后再绕Z轴，将X” 旋转到X轴方向。由于三坐标轴的正交关系，经最后一次旋转的Y”’必位于Y轴上。

98 2.4.2 空间直角坐标系之间的旋转变换 坐标变换公式为： 旋转矩阵： 是正交矩阵。

99 2.4.2 空间直角坐标系之间的旋转变换 若表示Z与Z’之间的夹角，则有： 若、分别表示X与X’和Y与Y’之间的夹角，则有：

100 2.4.2 空间直角坐标系之间的旋转变换 方法二： 将X’、Y’、Z’转换到X、Y、Z坐标系：
先绕X’将Y’旋转到YOZ平面与Y’OZ’平面的交线Y” ，再绕Y” 轴将Z ”旋转到Z轴，最后再绕Z轴，将X” 旋转到X轴方向。由于三坐标轴的正交关系，经最后一次旋转的Y”必位于Y轴上。

101 2.4.2 空间直角坐标系之间的旋转变换 坐标变换公式为： 其中，旋转矩阵： 是正交矩阵。

102 2.4.2 空间直角坐标系之间的旋转变换 若、、 分别表示X与X’、Y与Y’和Z与Z’之间的夹角，则有：

103 2.4.2 空间直角坐标系之间的旋转变换 当旋转角是小角度时，可略去其二次项，取： 旋转矩阵简化为： 坐标转换模型简化为：

104 2.4.2 空间直角坐标系之间的旋转变换 当旋转角较大时,因旋转矩阵是正交阵，满足条件： 若： 则根据正交条件，得：

105 2.4.2 空间直角坐标系之间的旋转变换 旋转矩阵中只有5个独立未知数。在进行坐标转换时，可以直接以旋转矩阵中的9个元素为未知数，加上6个约束条件直接解算。求得旋转矩阵元素后，进行坐标转换，不必解算旋转角。 这样可避免大旋转角时，线性化过程的复杂形式。

106 习 题 1、若采用克拉索夫斯基椭球，已知大地坐标： 计算三维空间坐标，并反算检核。 2、在上题中，大地经纬度和大地高分别变化了
习 题 1、若采用克拉索夫斯基椭球，已知大地坐标： 计算三维空间坐标，并反算检核。 2、在上题中，大地经纬度和大地高分别变化了 用微分公式计算三维空间坐标的变化量。 3、在球近似下，给出球心经纬度和高程与三维空间坐标的微分关系式。 4、若要求相对误差小于10-7，则当旋转角超过多少时，不能采用略去二次项的线性近似。

107 §2.4 空间大地直角坐标系及其转换模型 2.4.3 站心地平坐标系及其应用 1、站心地平直角坐标系与空间大地直角坐标系的转换关系
§2.4 空间大地直角坐标系及其转换模型 2.4.3 站心地平坐标系及其应用 1、站心地平直角坐标系与空间大地直角坐标系的转换关系 定义：站心点的法线为z轴，向上为正在地平面上以子午线方向为x轴，y与x、z轴正交，指向以东为正。 x z y 将站心坐标轴 xyz 变换成与空间坐标系的指向一致，需要如下几步： (1). z 坐标轴反向； (2). 绕y轴90。+B； (3). 绕z轴旋转-L。

108 2.4.3 站心地平坐标系及其应用 x z y 将站心系坐标轴变换到与三维空间直角坐标轴指向一致时的旋转矩阵为：
顾及，站心系原点在空间坐标系中的坐标为：

109 2.4.3 站心地平坐标系及其应用 则，站心系坐标到空间直角坐标系的变换公式为：

110 2.4.3 站心地平坐标系及其应用 由上式得，空间直角坐标系到站心系的变换公式为：

111 2.4.3 站心地平坐标系及其应用 2、站心极坐标系与站心地平直角坐标系的关系
定义：以站心系原点到点的空间距离、方位角和天顶距为坐标参数来确定三维点位，称为站心极坐标系。 由上式，得：

112 2.4.3 站心地平坐标系及其应用 也可以用以下公式计算：
公式中的天顶距和方位角都归算到以法线为基准。测量时以垂线为基准的，需要作垂线偏差改正。改正公式下面将讲到。

113 2.4.3 站心地平坐标系及其应用 3、空间直角坐标系与站心地平直角坐标系的旋转矢量之间的关系
若x、y和z为空间坐标系的旋转矢量， x、 y和z为站心坐标系的旋转矢量。顾及旋转矢量是平移不变量，旋转关系与坐标矢量相同。

114 2.4.3 站心地平坐标系及其应用 4、站心地平直角坐标系的应用 (1). 计算基线向量的大地方位角
其中，B0,L0为基线始端的纬度和经度。 (2). 绕站心系坐标轴的旋转向量有特殊意义 z 为方位旋转角。因涉及大地方位角起始方向的变动

115 2.4.3 站心地平坐标系及其应用 (4). 计算卫星的高度角和方位角 卫星Q的方位角和高度角可用其站心坐标xQ、yQ计算。

116 2.4.4 两个空间大地直角坐标系间的转换模型 1、Bursa - Wolf 模型 转换参数包括三个平移参数、三个旋转参数与一个尺度参数。

117 2.4.4 两个空间大地直角坐标系间的转换模型 略去尺度参数和旋转参数的乘积项，上式可进一步简化为：
上式第二式常用于转换参数未知时，利用同名点在两个坐标系中的坐标计算转换参数。

118 2.4.4 两个空间大地直角坐标系间的转换模型 上式再应用于Pj，并与上式相减，得Pi与Pj两点坐标差的坐标变换模型如下：

119 2.4.4 两个空间大地直角坐标系间的转换模型 2、Molodensky 模型
如果旋转与尺度是相对于参考点PK，即以参考点PK作变换中心。则有Molodensky 模型。 旋转角为小角度时，上式可简化为：

120 2.4.4 两个空间大地直角坐标系间的转换模型 上式同样可以简化为求解转换参数的形式如下： 其中，
相应于Molodensky模型的坐标差的转换模型与Bursa-Wolf模型相同。

121 2.4.4 两个空间大地直角坐标系间的转换模型 3、范士转换模型
若旋转角是围绕参考点的站心地平坐标系的坐标轴，即为范士转换模型。将三维空间坐标系的旋转角与站心系旋转角的关系代入Molodensky模型，即得范士转换模型如下：

122 2.4.4 两个空间大地直角坐标系间的转换模型 4、卫星网与地面网之间的转换 卫星网精度高， 地面网平面坐标与高程点不重合。

123 2.4.5 大地坐标的微分公式 根据大地坐标与三维空间直角坐标间的微分公式：
大地直角坐标的变动是由于原点平移、坐标轴旋转和尺度变化引起。即： 代入上式，得大地坐标微分公式。

124 2.4.5 大地坐标的微分公式 大地坐标微分公式的矩阵形式可表示为：

125 习 题 1、给出站心坐标系的定义。 2、经过哪几步旋转和平移变换，可将站心系坐标变换到三维空间直角坐标系中。
习 题 1、给出站心坐标系的定义。 2、经过哪几步旋转和平移变换，可将站心系坐标变换到三维空间直角坐标系中。 3、导出两点的大地方位角、距离和天顶距与站心坐标的关系。 4、三维空间坐标变换有哪几种模型？各种模型间的差异在哪里？ 5、范士变换模型的旋转参数有什么意义？

126 §2.5 参心坐标系和参考椭球 2.5.1 垂线偏差与Laplace方程 1、天文经度、天文纬度和天文方位角
天文经度：包含测站垂线的子午面与起始子午面的夹角； 天文纬度：测站垂线的与赤道面的夹角； 天文方位角：包含测站垂线的子午面与测站垂线和照准面所张成的垂直面的夹角； 天文天顶距：测站垂线与观测方向的夹角

127 2.5.1 垂线偏差与Laplace方程 因地极移动(x,y)，观测的天文经纬度、方位角需要归算到地极原点，称为极移改正，其公式如下：
观测值在地面取得，归算到椭球面上时，天文纬度和方位角需要作如下改正：

128 2.5.1 垂线偏差与Laplace方程 2、垂线偏差和大地水准面差距 大地水准面 参考椭球面 数值积分，得：

129 2.5.1 垂线偏差与Laplace方程 3、垂线偏差公式和Laplace方位角
如图所示：xyz为大地站心坐标系，x1 y1 z1为天文站心坐标系。两者的关系为： 1

130 2.5.1 垂线偏差与Laplace方程 天文和大地坐标系分别与原点在站心，坐标轴与三维空间直角坐标系指向相同的坐标系的关系如下：

131 2.5.1 垂线偏差与Laplace方程 由上面第一式代入第二式，略去高次项，整理得： 上式与 式相比较，得： 1

132 2.5.1 垂线偏差与Laplace方程 并得出Laplace方程： 顾及天文站心系（x1，y1，z1）与大地站心系（x，y，z）的关系：
和天顶距、方位角和站心坐标的关系：

133 2.5.1 垂线偏差与Laplace方程 将第二式代入第一式，得： 将展开式：

134 2.5.1 垂线偏差与Laplace方程 代入上式，并略去二次以上的项，得： 由第三式，得： 由第一式或第二式，顾及上式，并略去高次项得：

135 2.5.1 垂线偏差与Laplace方程 如果椭球短轴不平行与地轴，大地起始子午面不平行大地起始子午面，则还要考虑三个旋转角的影响，此时，大地经纬度和方位角与天文经纬度和方位角的关系可推广为：

136 2.5.2 参考椭球的定位和定向 1、椭球定位和定向的意义和条件 椭球定位：确定椭球中心的位置，即三个平移量
参考椭球的定位和定向 1、椭球定位和定向的意义和条件 椭球定位：确定椭球中心的位置，即三个平移量 椭球定向：确定椭球坐标轴的指向，即三个旋转量 参考椭球定位、定向应满足的条件： （1）椭球短轴与指定历元的地球自转轴平行； （2）大地起始子午面与天文起始子午面平行； （3）在一定区域内椭球面与大地水准面最为密合。 相应的数学表达式为：

137 参考椭球的定位和定向 2、椭球定位和定向的方法 在大地原点Pk存在关系： 上式隐含了三个旋转角为0。

138 2.5.2 参考椭球的定位和定向 （1）单点定位 大地测量工作刚开始，没有充分资料确定垂线偏差和大地水准面差距，假设在大地原点处，
参考椭球的定位和定向 （1）单点定位 大地测量工作刚开始，没有充分资料确定垂线偏差和大地水准面差距，假设在大地原点处， k = k = Nk = 0。 则有： 表示：单点定位时大地原点的法线与垂线一致，大地高等于正常高

139 2.5.2 参考椭球的定位和定向 （2）多点定位 根据以前的天文大地测量成果，来对椭球进行定位、定向和计算椭球元素。
参考椭球的定位和定向 （2）多点定位 根据以前的天文大地测量成果，来对椭球进行定位、定向和计算椭球元素。 由前面的大地坐标微分公式，

140 参考椭球的定位和定向 顾及： 展开大地坐标微分公式，略去旋转参数项，取最后一式代入上式，得： 根据条件： 求解定位参数。

141 2.5.2 参考椭球的定位和定向 也可以根据垂线偏差关系： 将大地坐标微分公式，略去旋转参数项，展开后的前二式代入上式。根据条件：
参考椭球的定位和定向 也可以根据垂线偏差关系： 将大地坐标微分公式，略去旋转参数项，展开后的前二式代入上式。根据条件： 求解椭球的定位参数。

142 参心坐标系的建立 参考椭球定位后，采用天文大地测量方法，进行三角测量、导线测量、空间大地测量（GPS、SLR、SALT、VLBI）、精密水准测量、重力测量等手段，建立参心坐标系，求得覆盖全区域的大地控制测量数成果。

143 §2.6 协议地球参考系(CTRS)和平均地球椭球
定义：坐标轴指向BIH 系统，坐标轴定向使地 壳各方向的运动之和为0。 建立：由各板块的SLR、VLBI站以及部分IGS站确定。

144 2.6.2 当今技术条件下的平均椭球 利用低轨卫星观测数据求定平均椭球的四个参数，并对平均椭球进行定位和定向。

145 习 题 1、天文方位角归算到参考椭球上要加哪些改正？写出其改正公式。
习 题 1、天文方位角归算到参考椭球上要加哪些改正？写出其改正公式。 2、若某测站的垂线偏差为 = 7.02,  = -4.38,测得某方向的天顶距为854856.3, 已知该方向的方位角为1354325.6。求：改正后的天顶距。 3、参考椭球定位满足什么条件？有哪两种方法？