ITRS 与 GCRS 之间的坐标转换

1. ITRS 与 GCRS 之间的坐标转换

由于地球自转, 地球坐标系并不是一个惯性坐标系, 而轨道计算是建立在牛顿力学的基础上的, 因此定轨工作不能在地球坐标系中进行。如前所述, GCRS 是一个相当不错的准惯性坐标系, 定轨工作一般都在该坐标系中进行，但是用户利用卫星导航定位系统最终是为了求得在地球坐标系中的位置和速度, 因而还必须把 GCRS 中所求得的卫星轨道(卫星位置和速度) 转换到地球坐标系 ITRS(WGS 84) 中去。

ITRS 与 GCRS 之间有下列转换关系:
(XYZ)GCRS=[P][N][R][W](XYZ)ITRS\left(\begin{array}{l} X \\ Y \\ Z \end{array}\right)_{\mathrm{GCRS}}=[\boldsymbol{P}][\boldsymbol{N}][\boldsymbol{R}][\boldsymbol{W}]\left(\begin{array}{l} X \\ Y \\ Z \end{array}\right)_{\mathrm{ITRS}} ⎝⎛XYZ⎠⎞GCRS=[P][N][R][W]⎝⎛XYZ⎠⎞ITRS
(XYZ)ITRS =[W]−1[R]−1[N]−1[P]−1(XYZ)GCRS\left(\begin{array}{l} X \\ Y \\ Z \end{array}\right)_{\text {ITRS }}=[\boldsymbol{W}]^{-1}[\boldsymbol{R}]^{-1}[\boldsymbol{N}]^{-1}[\boldsymbol{P}]^{-1}\left(\begin{array}{l} X \\ Y \\ Z \end{array}\right)_{\mathrm{GCRS}} ⎝⎛XYZ⎠⎞ITRS =[W]−1[R]−1[N]−1[P]−1⎝⎛XYZ⎠⎞GCRS
式中, [P][\boldsymbol{P}][P] 为岁差矩阵; [N][\boldsymbol{N}][N] 为章动矩阵; [R][\boldsymbol{R}][R] 为地球自转矩阵; [W][\boldsymbol{W}][W] 为极移矩阵。考虑到 IGS 已完成了坐标转换工作, 在精密星历中直接给出了卫星质心在 ITRS 中的位置和速度, 而广播星历的精度有限, 允许采用一些近似的转换方法, 因此在下面的坐标转换中, 我们仍采用经典的转换方法与术语 (与 IS-GPS-200D 及 IS-GPS-705 中给出的方法基本一致)。高精度的严格方法可参阅 IAU 的决议文件和空间大地测量学等参考资料。

（1）把 GCRS 转换至观测时刻 tit_iti 的平天球坐标系
我们知道, GCRS 是参考时刻 t0=J2000.0t_0=\mathrm{J} 2000.0t0=J2000.0 时的平天球坐标系, 要把它转换为观测时刻 tit_iti 时的平天球坐标系, 只要考虑 [t0−ti]\left[t_0-t_i\right][t0−ti] 时间段内的岁差改正, 即乘上 [P]−1[\boldsymbol{P}]^{-1}[P]−1 矩阵即可。

（2）把 tit_iti 时的平天球坐标系转换为同一时刻的真天球坐标系
要把观测时刻 tit_iti 时的平天球坐标系转换为真天球坐标系, 只需顾及该时刻的章动, 即只需乘上 [N]−1[\boldsymbol{N}]^{-1}[N]−1 矩阵即可。

（3）把 tit_iti 时的真天球坐标系转换为同一时刻的真地球坐标系
我们知道, 真天球坐标系 XXX 轴是指向该时刻的真春分点 γ\gammaγ 的,而真地球坐标系的 XXX 轴是指向起始子午线与赤道的交点,两者之间的夹角称为格林尼治真恒星时 GAST 。其计算公式如下:
GAST =360∘24h(UT1+6h41m50.54841s+8640184.812866s⋅t+0.093104s⋅t2−6.2s×10−6⋅t3)+ΔΨcos⁡(εˉ+Δε)\begin{aligned} \text { GAST }=& \frac{360^{\circ}}{24^{\mathrm{h}}}\left(\mathrm{UT} 1+6 \mathrm{~h} 41 \mathrm{~m} 50.54841 \mathrm{~s}+8640184.812866 \mathrm{~s} \cdot \mathrm{t}+0.093104 \mathrm{~s} \cdot t^2\right.\\ &\left.-6.2 \mathrm{~s} \times 10^{-6} \cdot t^3\right)+\Delta \Psi_{\cos }(\bar{\varepsilon}+\Delta \varepsilon) \end{aligned} GAST =24h360∘(UT1+6 h41 m50.54841 s+8640184.812866 s⋅t+0.093104 s⋅t2−6.2 s×10−6⋅t3)+ΔΨcos(εˉ+Δε)
式中, ttt 为离 J2000.0\mathrm{J} 2000.0J2000.0 的儒略世纪数; εˉ\bar{\varepsilon}εˉ 为仅顾及岁差时的黄赤交角, εˉ=23∘26′21.448′′−\bar{\varepsilon}=23^{\circ} 26^{\prime} 21.448^{\prime \prime}-εˉ=23∘26′21.448′′− 46.815′′⋅t−0.00059′′⋅t2+0.001813′′⋅t3;Δψ46.815^{\prime \prime} \cdot t-0.00059^{\prime \prime} \cdot t^2+0.001813^{\prime \prime} \cdot t^3 ; \Delta \psi46.815′′⋅t−0.00059′′⋅t2+0.001813′′⋅t3;Δψ46.815′′⋅t−0.00059′′⋅t2+0.001813′′⋅t3;ΔΨ46.815^{\prime \prime} \cdot t-0.00059^{\prime \prime} \cdot t^2+0.001813^{\prime \prime} \cdot t^3 ; \Delta \Psi46.815′′⋅t−0.00059′′⋅t2+0.001813′′⋅t3;ΔΨ为黄经章动; Δε\Delta \varepsilonΔε 为交角章动; UT1\mathrm{UT1}UT1 则可据观测时的 UTC 和( UTC-UT1) 值求得。

把真天球坐标系绕 ZZZ 轴旋转 GAST 角后就能转换到真地球坐标系, 旋转矩阵 R\boldsymbol{R}R 为:
R=(cos⁡GASTsinGAST⁡0−sinGAST⁡cos⁡GAST0001)\boldsymbol{R}=\left(\begin{array}{ccc} \cos G A S T & \operatorname{sinGAST} & 0 \\ -\operatorname{sinGAST} & \cos G A S T & 0 \\ 0 & 0 & 1 \end{array}\right) R=⎝⎛cosGAST−sinGAST0sinGASTcosGAST0001⎠⎞
（4）把 tit_iti 时的真地球坐标系转换为 ITRS (WGS 84)

从上图可以看出, 只需要将 tit_iti 时的真地球坐标系绕 yyy 轴旋转 (−Xp)\left(-X_p\right)(−Xp) 角后,然后再绕 xxx 轴旋转 (−Yp)\left(-Y_p\right)(−Yp) 角后, 就可以把直地球坐标系 O−xyzO-x y zO−xyz 转换为 ITRS⁡(WGS84)\operatorname{ITRS}(\mathrm{WGS} 84)ITRS(WGS84) 坐标系 O−XYZ。 O-X Y Z_{\text {。 }}O−XYZ。即

(XYZ)=Rx(−Yp)Ry(−Xp)(xyz)=(1000cos⁡Yp−sin⁡Yp0sin⁡Ypcos⁡Yp)(cos⁡Xp0sin⁡Xp010−sin⁡Xp0cos⁡Xp)(xyz)\left(\begin{array}{l} X \\ Y \\ Z \end{array}\right)=\boldsymbol{R}_x\left(-Y_p\right) \boldsymbol{R}_y\left(-X_p\right)\left(\begin{array}{l} x \\ y \\ z \end{array}\right)=\left(\begin{array}{ccc} 1 & 0 & 0 \\ 0 & \cos Y_p & -\sin Y_p \\ 0 & \sin Y_p & \cos Y_p \end{array}\right)\left(\begin{array}{ccc} \cos X_p & 0 & \sin X_p \\ 0 & 1 & 0 \\ -\sin X_p & 0 & \cos X_p \end{array}\right)\left(\begin{array}{l} x \\ y \\ z \end{array}\right) ⎝⎛XYZ⎠⎞=Rx(−Yp)Ry(−Xp)⎝⎛xyz⎠⎞=⎝⎛1000cosYpsinYp0−sinYpcosYp⎠⎞⎝⎛cosXp0−sinXp010sinXp0cosXp⎠⎞⎝⎛xyz⎠⎞
由于极移值 XpX_pXp 和 YpY_pYp 都是小于 0.5′′0.5^{\prime \prime}0.5′′ 的微小值, 所以 cos⁡Xp=cos⁡Yp=1,sin⁡Xp=Xp\cos X_p=\cos Y_p=1, \sin X_p=X_pcosXp=cosYp=1,sinXp=Xp, sin⁡Yp=Yp\sin Y_p=Y_psinYp=Yp, 于是有 :::
(XYZ)=(10001−Yp0Yp1)(10Xp010−Xp01)(xyz)=(10Xp01−Yp−XpYp1)(xyz)=[W](xyz)\left(\begin{array}{l} X \\ Y \\ Z \end{array}\right)=\left(\begin{array}{ccc} 1 & 0 & 0 \\ 0 & 1 & -Y_p \\ 0 & Y_p & 1 \end{array}\right)\left(\begin{array}{ccc} 1 & 0 & X_p \\ 0 & 1 & 0 \\ -X_p & 0 & 1 \end{array}\right)\left(\begin{array}{l} x \\ y \\ z \end{array}\right)=\left(\begin{array}{ccc} 1 & 0 & X_p \\ 0 & 1 & -Y_p \\ -X_p & Y_p & 1 \end{array}\right)\left(\begin{array}{l} x \\ y \\ z \end{array}\right)=[\boldsymbol{W}]\left(\begin{array}{l} x \\ y \\ z \end{array}\right) ⎝⎛XYZ⎠⎞=⎝⎛10001Yp0−Yp1⎠⎞⎝⎛10−Xp010Xp01⎠⎞⎝⎛xyz⎠⎞=⎝⎛10−Xp01YpXp−Yp1⎠⎞⎝⎛xyz⎠⎞=[W]⎝⎛xyz⎠⎞

出自：《GPS测量与数据处理》第二章。

ITRS 与 GCRS 之间的坐标转换相关推荐

RTKLIB学习总结（二）时间系统、坐标系统
上篇博客: RTKLIB学习总结(一)MarkDown.VS配置RTKLIB.矩阵.最小二乘和Kalman滤波基本函数关于时间系统和坐标系统的算法概念也可以看RTKLIB的manual131面一. ...
ICRS、GCRS、CIRS、TIRS和ITRS坐标系统简介
ICRS.GCRS.CIRS.TIRS和ITRS坐标系统简介 1. 简介 2. ICRS.GCRS.CIRS.TIRS和ITRS分别介绍 2.1 ICRS详细说明 2.2 GCRS详细说明 2.3 C ...
arcgis 投影变换与坐标转换研究
arcgis 投影变换与坐标转换研究 1 ArcGIS中的投影方法投影的方法可以使带某种坐标信息数据源进行向另一坐标系统做转换,并对源数据中的X和Y值进行修改.我们生产实践中一个典型的例子是利用该方 ...
ArcGIS中的投影和坐标转换
ArcGIS中的投影和坐标转换 1 ArcGIS中坐标系统的定义一般情况下地理数据库(如Personal GeoDatabase的 Feature DataSet .Shape File等)在创建 ...
多相机坐标转换_使用KITTI数据集实现坐标转换
本文的主要内容,是讨论如何将KITTI激光雷达的点云数据,映射到图像中去. KITTI数据集是怎么采集的下面这幅图是车的俯视图,可以看到,用到了2个彩色相机,2个灰度相机,1个激光雷达.1个GPS/ ...
python坐标系转换库_转载：python库Pyproj进行坐标转换
https://www.baidu.com 利用Pyproj进行坐标转换作者:郜科科两个坐标系统的参考椭球不同,实地一个点的不同坐标系的值是不同的,不同的部门采用的坐标系统经常是不一致,所以要转换 ...
gps84转换gcj02公式_【测绘基础】关于投影坐标和坐标转换，这几点你不得不知道！...
[测绘基础]关于投影坐标和坐标转换,这几点你不得不知道! 本文要点: 1.地理坐标系是什么?常见的地理坐标系有哪些? 2.投影坐标系是什么?常见的投影坐标系有哪些? 3.关于坐标转换. 4.如何在Ar ...
坐标系之间的旋转平移变换与对应变换矩阵的关系
在摄影测量和计算机视觉中,经常会遇到空间坐标系之间的坐标转换问题,而两个坐标系之间的变换关系一般可以通过一个旋转矩阵R和一个平移向量T(或C)描述.因此,理解清楚坐标系之间旋转平移的转换过程与对应变换 ...
实现坐标转换--不同地理位置系统转换入门
最近做离线地图的程序,找了一些地图坐标的转换方法. 很多流行的应用程序可以提供基于位置的服务,但是计算机如何识别真实世界中的位置呢?很多方法都涉及到地理坐标系统,并且在实际应用中存在不同的此类系统.在 ...

ITRS 与 GCRS 之间的坐标转换

1. ITRS 与 GCRS 之间的坐标转换

ITRS 与 GCRS 之间的坐标转换相关推荐

最新文章

热门文章