地球运转、时间系统与坐标系统：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中每一个坐标转换、时间函数和精度指标背后的物理意义。
