目录
- 1. 国际地球参考系统 (ITRS) 与参考框架 (ITRF)
- 2. 坐标框架的转换与历元归算
- 3. 地球物理因素对坐标的影响
- 4. WGS-84 坐标系统及其精化
- 5. 天球坐标系与地球坐标系的转换
- AI 总结
1. 国际地球参考系统 (ITRS) 与参考框架 (ITRF) 原片 @ 00:00*
Screenshot-[00:00]
1.1 概念引入
在介绍完天球坐标系后,课程转入地球坐标系。地球坐标系主要用于描述地面点的位置以及近地空间物体(如飞机、卫星)的位置。目前最常用的两种体系是:
- 国际地球参考系统 (ITRS, International Terrestrial Reference System):理论层面的定义。
- 国际地球参考框架 (ITRF, International Terrestrial Reference Frame):实际层面的实现。
1.2 ITRS (系统) 的定义 原片 @ 00:51*
ITRS 是一个理论上的规定,用户无法直接使用“定义”,必须通过框架来实现。其主要定义标准包括:
- 原点:包括海洋和大气的整个地球的质心。
- 尺度:采用国际单位制 (SI) 的米。
- 定向:初始指向与 BIH 1984.0 定义的方向重合。
- 旋转条件:规定整个地球地壳没有整体的旋转(即满足 No-net-rotation 条件)。
- 由于板块运动,各板块向不同方向移动,ITRS 规定所有板块运动的平均角动量为零。
- 如果存在整体旋转,则意味着地球自转速度定义的改变(这部分已被归入 UT1 时间系统处理,不应体现在坐标系旋转中)。
Screenshot-[01:12]
1.3 ITRF (框架) 的实现 原片 @ 02:24*
用户使用坐标系时,依赖的是 IERS (国际地球自转服务) 发布的坐标框架。
- 构成:由分布在全世界的观测站(VLBI, SLR, GPS 等)组成的一个网。
- 发布内容:IERS 解算并公布这些站点的坐标 (Coordinates)、坐标变化率 (Velocities) 以及地球定向参数 (EOP)。
- 用户应用:用户通过联测这些控制点,将自己的测量成果纳入该坐标框架中。
2. 坐标框架的转换与历元归算 原片 @ 02:53*
Screenshot-[02:53]
2.1 框架的精度演变
随着观测技术的进步(如 VLBI, GPS 精度提高)和数据处理方法的改进,ITRF 框架的精度越来越高。
- 早期:数据积累少(1-2年),精度较差(原点偏差可能达 20-30 厘米)。
- 现在:积累了数十年资料,框架非常稳定。
- 更新频率:早期每年发布一次,现在框架稳定后约 3-5 年发布一次(如 ITRF2000, ITRF2008, ITRF2014 等)。
2.2 坐标转换模型 (布尔沙模型) 原片 @ 04:19*
不同坐标框架之间(例如 ITRF2000 到 ITRF2005)的转换,通常采用 布尔沙模型 (Bursa Model),包含 7 个参数:
- 3 个平移参数 ($T_x, T_y, T_z$):解决原点不重合。
- 1 个尺度参数 ($D$ 或 $k$):解决尺度不一致。
- 3 个旋转参数 ($R_x, R_y, R_z$):解决坐标轴指向差异。
转换公式如下:
$$
\begin{bmatrix} X \ Y \ Z \end{bmatrix}{New} = \begin{bmatrix} T_x \ T_y \ T_z \end{bmatrix} + (1 + D) \cdot \begin{bmatrix} 1 & R_z & -R_y \ -R_z & 1 & R_x \ R_y & -R_x & 1 \end{bmatrix} \cdot \begin{bmatrix} X \ Y \ Z \end{bmatrix}{Old}
$$
Screenshot-[04:34]
2.3 历元归算 (Epoch Propagation) 原片 @ 05:12*
ITRF 框架给出的坐标是基于特定 参考历元 (Reference Epoch, $t_0$) 的。
- IERS 提供:参考历元时刻的坐标 ($X_0, Y_0, Z_0$) 和 站速度 ($\dot{X}, \dot{Y}, \dot{Z}$)。
- 计算当前历元坐标:由于板块运动,站点位置随时间变化。如果要求解当前时刻 $t$ 的坐标,必须进行历元归算:
$$
X(t) = X(t_0) + \dot{X} \cdot (t – t_0)
$$
注意:不能直接使用参考历元的坐标,必须考虑从参考历元到当前时刻的时间段内,板块运动引起的位移。
3. 地球物理因素对坐标的影响 原片 @ 06:52*
Screenshot-[06:52]
3.1 瞬时坐标与平均坐标
IERS 给出的站坐标是剔除了各种周期性地球物理影响后的“平均”坐标。但在高精度 GPS 测量中,必须考虑地球瞬时的形变。
3.2 主要影响因素
- 固体潮 (Solid Earth Tides):
- 地球不是刚体,具有弹性。在日月引潮力作用下,地壳会发生隆起和下陷。
- 量级:高程方向可达 50 cm 以上,平面方向约几厘米。
- 海潮负荷 (Ocean Loading):
- 海洋潮汐引起海水质量重新分布,压在海底地壳上,导致沿海地区地面高程随潮汐变化。
- 大气负荷 (Atmospheric Loading):大气压力的变化导致地壳形变。
- 冰雪负荷与回弹:积雪压低地面,融化后地面回弹(季节性变化)。
3.3 数据处理策略 原片 @ 09:37*
- IERS 发布的坐标是 Tide-Free (无潮) 或剔除了周期性影响的。
- 用户计算时:
- 使用 IERS 提供的坐标和速度算出当前历元的平均坐标。
- 利用地球物理模型(固体潮模型、海潮模型等)计算当前的改正数。
- 将改正数加回到平均坐标上,得到该时刻的瞬时坐标。
- 这些复杂的模型计算通常由专业的 GPS 解算软件(如 GAMIT, Bernese)自动完成。
4. WGS-84 坐标系统及其精化 原片 @ 10:30*
Screenshot-[10:42]
4.1 WGS-84 概述
- 定义:世界大地坐标系 (World Geodetic System 1984),由美国国防部 (DoD) 建立。
- 演变:WGS-60 -> WGS-66 -> WGS-72 -> WGS-84。
- 理论定义:与 ITRS 的四条定义完全一致(原点、尺度、定向、旋转)。
4.2 WGS-84 与 ITRF 的关系
- 实现方式不同:ITRF 由 IERS 组织通过全球多种技术(VLBI, SLR, GPS)综合实现;WGS-84 最初主要基于 GPS 广播星历系统,但也利用了多普勒、SLR 等资料。
- 早期差异:早期的 WGS-84 实现精度较低,与 ITRF 可能存在米级差异。
4.3 WGS-84 的精化版本 原片 @ 12:41*
为了提高精度并与国际标准接轨,美国国防部多次对 WGS-84 进行精化 (Realization),使其逼近 ITRF:
- WGS-84 (G730):
- G 代表 GPS,730 代表 GPS 第 730 周(1994年1月2日启用)。
- 方法:将 GPS 跟踪站与 ITRF 框架点联测平差,固定 ITRF 点坐标。
- 精度:与 ITRF 偏差约 10 cm。
- WGS-84 (G873):
- 启用时间:1996年9月29日。
- 精度:与 ITRF 偏差约 2 cm。
- WGS-84 (G1150):
- 启用时间:2002年1月20日。
- 精度:与 ITRF 偏差约 1 cm。
Screenshot-[15:41]
4.4 结论
对于绝大多数工程应用(非毫米级科学研究),目前的 WGS-84 与 ITRF 可以视为同一坐标系,混用不会带来显著误差。
5. 天球坐标系与地球坐标系的转换 原片 @ 16:36*
Screenshot-[17:12]
5.1 转换的必要性
- 卫星定轨:在天球坐标系 (GCRS) 中进行。因为卫星运动遵循牛顿力学,需要一个惯性系(GCRS 可视为准惯性系,坐标轴指向空间固定方向,如 J2000.0 平春分点和平北天极)。
- 用户应用:用户位于地球表面,需要地球坐标系 (ITRS) 中的坐标。
- 因此,必须在惯性系(计算轨道)和地固系(用户位置/观测站位置)之间进行转换。
5.2 转换矩阵 原片 @ 21:18*
已知地球坐标系向量 $X_{TRS}$,求天球坐标系向量 $X_{CRS}$(或反之),需要经过一系列旋转矩阵。虽然视频中提到顺序可能颠倒,但标准转换流程涉及四个主要部分:
$$
X_{TRS} = [W] \cdot [S] \cdot [N] \cdot [P] \cdot X_{CRS}
$$
其中:
- $P$ (Precession, 岁差):描述地球自转轴在空间中的长期圆锥运动。
- $N$ (Nutation, 章动):描述叠加在岁差上的短周期微小震动(多达1000多项)。
- $S$ (Spin/Rotation, 地球自转):描述地球绕自转轴的旋转(涉及格林尼治恒星时 GAST)。
- $W$ (Polar Motion, 极移):描述瞬时自转轴与地壳之间的相对运动。
Screenshot-[22:36]
5.3 计算复杂性
- 章动矩阵包含一千多项正弦/余弦项,系数取决于日月及行星位置,计算非常繁琐。
- 极移参数需要实测。
- 实际操作:对于研究生课程或工程应用,不需要手动计算这些矩阵。通常使用现成的权威软件(如 IERS 提供的 SOFA 库)或通过 GPS 软件自动处理。用户主要关注地球坐标系下的结果。
AI 总结
本节课程深入讲解了 GPS 测量中至关重要的坐标系统原理,重点涵盖了以下内容:
- ITRS 与 ITRF:区分了理论定义的系统 (ITRS) 与实际实现的框架 (ITRF)。强调了 ITRF 框架随时间演变,且包含速度场,使用时必须结合历元归算 ($X_t = X_0 + V \cdot \Delta t$)。
- 地球物理改正:指出了 IERS 坐标是“平均”坐标,实际高精度应用需通过模型恢复固体潮(高达50cm)、海潮负荷等引起的瞬时形变。
- WGS-84 的演进:WGS-84 虽然定义上与 ITRS 一致,但早期实现精度较低。经过 G730, G873, G1150 等多次精化,目前 WGS-84 与 ITRF 的差异已在厘米级以内,常规应用可不予区分。
- 坐标系转换:阐述了从天球坐标系(惯性系,用于定轨)到地球坐标系(地固系,用于定位)的转换逻辑,涉及岁差、章动、地球自转、极移四个复杂的旋转矩阵。
该笔记适合需要深入理解 GPS 坐标基准、框架转换及高精度数据处理背景的研究生复习使用。