从轨道平面到ECEF：卫星坐标旋转的RTKLIB实现

从轨道平面到ECEF：卫星坐标旋转的RTKLIB实现

在阅读RTKLIB源码时，卫星位置计算中有一个关键步骤：把卫星在轨道平面内的坐标 (x, y)转换到地固系（ECEF）的(X,Y,Z)。这个转换涉及两次旋转——从轨道平面到惯性系，再从惯性系到地固系。RTKLIB 巧妙地把这两次旋转合并成了一步，用修改后的升交点赤经 O 一次性完成。本文尝试把这一过程彻底讲清楚。

一、为什么需要坐标旋转？

广播星历计算卫星位置时，先得到的是卫星在轨道平面内的坐标：

其中 r是向径，u 是纬度幅角。

但用户最终需要的是卫星在地固系（ECEF）中的坐标 (X,Y,Z)，用于计算接收机的位置。因此必须通过旋转，把卫星从轨道平面“搬”到地球坐标系中。

二、涉及的三个坐标系

三、第一步：轨道平面 → 惯性系

从轨道平面 (x,y) 到惯性系 (Xi​,Yi​,Zi​)，需要两次旋转：

旋转1：绕 x 轴转倾角 i

轨道平面与赤道面有夹角 i，所以先绕 x 轴旋转：

旋转2：绕 Z 轴转升交点赤经 

把轨道平面从“升交点方向”转到惯性系中的实际方位：

这一步得到的是惯性系中的卫星坐标。 注意 OMG0 就是在惯性系中定义的升交点赤经。

四、第二步：惯性系 → 地固系

地固系相对于惯性系，绕 Z 轴旋转了地球自转角度 θ：

其中：

• toes：toe 时刻的周内秒（从本周起点累计的自转）

• tk：从 toe 到当前时刻的时间差

• omge：地球自转角速度

从惯性系到地固系的旋转矩阵为：

五、RTKLIB 的合并技巧：一步到位

如果把第一步的 (Xi​,Yi​) 代入第二步的旋转公式，并用三角恒等式合并，会得到：

令：

这正是 RTKLIB 中的：

O = OMG0 + (OMGd – omge) * tk – omge * toes;

于是最终公式变成：

rs[0] = x * cos(O) – y * cosi * sin(O);
rs[1] = x * sin(O) + y * cosi * cos(O);
rs[2] = y * sin(i);

RTKLIB 没有显式地先转惯性系再转地固系，而是直接把“地球自转角度”从升交点赤经里扣掉，得到地固系中的升交点赤经 O，然后用同一套旋转公式一次性算出 ECEF 坐标。数学上完全等价，但代码更简洁、计算量更小。

六、总结：两步 vs. 一步

核心公式：

它把“惯性系中的升交点赤经”转换成了“地固系中的升交点赤经”，让后续的坐标旋转一步到位。