地球运转、时间系统与坐标系统：GNSS大地测量学基础总结

在学习RTKLIB和GNSS定位的过程中，坐标系统和时间系统是最容易被忽略却又最基础的知识。不理解地球的运动规律，就无法理解为什么要定义这么多时间系统和坐标框架；不理解这些框架之间的转换关系，就无法真正读懂卫星位置计算和定位解算的代码。本文系统整理了地球运转、时间系统和坐标系统的核心知识，希望能为GNSS学习者提供一份清晰的参考。

一、地球的运转

1.1 地球自转

地球绕自转轴自西向东旋转，周期约为24小时。自转产生以下效应：

昼夜交替
恒星日与太阳日的差异
地球自转线速度 ：赤道处约 465 m/s

地球自转是GNSS中最重要的运动之一，因为用户在地球上随地球一起转动，而卫星在惯性空间中运动，两者之间的相对运动必须通过坐标转换来处理。

1.2 地球公转

地球绕太阳公转，轨道周期约 365.25天。公转产生以下效应：

四季变化
太阳在星空背景中每天东移约1°
黄道：地球公转轨道面在天球上的投影

公转导致太阳相对于恒星背景的位置每天变化约1°，这直接造成了恒星日与太阳日相差约4分钟。

1.3 岁差

地球自转轴在空间中并非固定不动，而是像陀螺一样缓慢圆锥运动，周期约 26000年 。这种现象称为岁差。

岁差导致：

春分点沿黄道缓慢西退（每年约50.3角秒）
天极位置变化
天球坐标系随时间缓慢旋转

1.4 章动

在岁差的基础上，地球自转轴还有短周期的微小摆动，称为章动。主周期约 18.6年，由日月引力引起。

岁差和章动共同决定了瞬时天极 相对于协议天极的位置。

1.5 极移

地球自转轴相对于地壳本身也在微小移动，称为极移。极移的周期约1年，幅度约几米。

极移导致瞬时地球坐标系 与协议地球坐标系之间存在微小差异。

二、时间系统

2.1 为什么需要多种时间系统？

不同的应用需要不同的时间基准：

需求	对应时间系统
天文观测	恒星时
日常生活	世界时
行星运动计算	历书时
高精度计量	原子时
GNSS	GPST、BDT

2.2 恒星时

定义：以春分点为参考，地球自转一周的时间。

恒星日 ：约 23h 56m 4.09s
用途：天文观测、惯性系与地固系的转换
特点：不均匀（受地球自转不均匀影响）

2.3 世界时

定义：以太阳为参考，地球自转一周的时间。

平太阳日：24小时
UT1：实际世界时，反映地球真实自转
UTC：协调世界时，用原子时秒长，通过闰秒与UT1保持接近
关系：1 平太阳日 = 1.0027379 恒星日

2.4 历书时

定义：以太阳系天体的牛顿运动为基准的时间。

提出原因：地球自转不均匀，世界时不能作为精确时间尺度
被取代：1984年后被地球时（TT）和质心力学时（TDB）取代
用途：历史星历表、行星轨道计算

2.5 原子时

定义：以铯-133原子基态超精细跃迁频率定义秒长。

秒长：9,192,631,770个周期
TAI：国际原子时，由全球原子钟加权平均
TT：地球时，TT = TAI + 32.184s
均匀性：极高，是当前最均匀的时间尺度

2.6 GNSS时间系统

时间系统	基准	与UTC关系
GPST	原子时	GPST = UTC + 闰秒（当前18s）
BDT	原子时	BDT = GPST - 14s
GLONASS时	UTC	随闰秒调整

2.7 时间系统转换关系

复制代码

原子时 (TAI)
    │ +32.184s
    ▼
地球时 (TT) ──用于星历、定轨
    │
    │ +闰秒
    ▼
协调世界时 (UTC) ──民用、GNSS时间标签
    │
    │ UT1-UTC (IERS)
    ▼
世界时 UT1 ──地球真实自转
    │
    │ 恒星时 = UT1 × 1.0027379 + ...
    ▼
恒星时 ──惯性系/地固系转换

三、坐标系统

3.1 天球坐标系 vs. 地球坐标系

对比项	天球坐标系	地球坐标系
参考对象	遥远恒星/河外射电源	地球本身
是否随地球自转	不随（惯性系）	随地球自转
Z轴指向	北天极	北极
X轴指向	春分点	本初子午线
用途	卫星轨道计算	地面定位

3.2 四个层次的概念

坐标系	基准	是否随时间变	用途
协议天球系（GCRS）	河外射电源	固定	惯性基准
瞬时天球系（CIRS）	真赤道/真春分点	随岁差章动	中间站
瞬时地球系（TIRS）	瞬时自转轴	随极移	中间站
协议地球系（ITRS）	地球质心/本初子午线	固定	定位输出

3.3 转换链条

从协议天球系到协议地球系，需要经过四个旋转：

复制代码

协议天球坐标系（GCRS）
    │ 岁差 + 章动
    ▼
瞬时天球坐标系（CIRS）
    │ 地球自转（GAST / ERA）
    ▼
瞬时地球坐标系（TIRS）
    │ 极移
    ▼
协议地球坐标系（ITRS）

矩阵形式：

\\text{ITRS}\] = \[W\] \\cdot \[R\] \\cdot \[N\] \\cdot \[P\] \\cdot \[\\text{GCRS}\]\[ITRS\]=\[W\]⋅\[R\]⋅\[N\]⋅\[P\]⋅\[GCRS

其中：

P：岁差矩阵
N：章动矩阵
R：地球自转矩阵
W：极移矩阵

3.4 地球定向参数（EOP）

连接天球系和地球系的桥梁是地球定向参数（EOP），包括：

岁差：自转轴长期漂移
章动：自转轴短周期摆动
极移：自转轴相对于地壳的移动
UT1-UTC：地球自转角度与UTC的差异

EOP由IERS发布，GNSS高精度定位需要使用。

四、参考框架

4.1 ICRF（国际天球参考框架）

定义：天球参考系的具体实现
基准：河外射电源（类星体）
维护者：IAU / IVS
实现技术：VLBI
精度：微角秒级

4.2 ITRF（国际地球参考框架）

定义：地固参考系的具体实现
基准：地球质心、地极、本初子午线
维护者：IERS
实现技术：VLBI + SLR + GNSS + DORIS
精度：毫米级

4.3 两者的关系

ICRF和ITRF通过EOP相互转换，是GNSS定位和地球科学研究的基础。

五、在GNSS/RTKLIB中的体现

概念	RTKLIB中的对应
GCRS	广播星历计算卫星位置时，隐含使用惯性系
ITRS	`ecef2pos`、`sol->rr` 输出的ECEF坐标
岁差章动	广播星历参数已吸收，SPP中不需显式计算
地球自转	`OMGE`、`omge*toes` 等项
极移	`navs.erp` 存储ERP数据，精密定位使用
时间转换	`utc2gpst`、`gpst2time`、`bdt2gpst`
闰秒	`leaps` 数组

关键认知 ：RTKLIB的广播星历（SPP/RTK）已经隐含完成了GCRS→ITRS的转换，用户不需要显式计算岁差、章动、极移。但在精密定轨和PPP中，如果需要高精度坐标转换，则需要显式使用ERP数据。

六、精度与时间的关系

距离误差 = 速度 × 时间误差。因此：

场景	速度	1cm精度对应时间误差
地球自转（赤道）	465 m/s	2.15 \times 10^{-5}2.15×10−5 s
导航卫星轨道	3870 m/s	2.58 \times 10^{-6}2.58×10−6 s
电磁波测距	3 \times 10^83×108 m/s	3.33 \times 10^{-11}3.33×10−11 s

这解释了为什么GNSS需要原子钟、皮秒级时间测量和相对论修正。

七、总结

地球的运转------自转、公转、岁差、章动、极移------决定了时间和坐标系统的定义。天球坐标系描述卫星在惯性空间中的运动，地球坐标系描述用户在地球上的位置，两者通过地球定向参数相互转换。

时间系统：恒星时、世界时、历书时、原子时，各有用途。
坐标系统：协议天球系、瞬时天球系、瞬时地球系、协议地球系，层层递进。
参考框架：ICRF和ITRF，分别是天球和地球参考系的具体实现。
GNSS：GPST、BDT基于原子时，定位结果输出在ITRF下。

理解这些基础，才能读懂RTKLIB中每一个坐标转换、时间函数和精度指标背后的物理意义。