ITRS 与 GCRS 之间的坐标转换
Posted 大作家佚名
tags:
篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了ITRS 与 GCRS 之间的坐标转换相关的知识，希望对你有一定的参考价值。
1. ITRS 与 GCRS 之间的坐标转换
由于地球自转, 地球坐标系并不是一个惯性坐标系, 而轨道计算是建立在牛顿力学的基础上的, 因此定轨工作不能在地球坐标系中进行。如前所述, GCRS 是一个相当不错的准惯性坐标系, 定轨工作一般都在该坐标系中进行，但是用户利用卫星导航定位系统最终是为了求得在地球坐标系中的位置和速度, 因而还必须把 GCRS 中所求得的卫星轨道(卫星位置和速度) 转换到地球坐标系 ITRS(WGS 84) 中去。
ITRS 与 GCRS 之间有下列转换关系:
(
X
Y
Z
)
G
C
R
S
=
[
P
]
[
N
]
[
R
]
[
W
]
(
X
Y
Z
)
I
T
R
S
\\left(\\beginarrayl X \\\\ Y \\\\ Z \\endarray\\right)_\\mathrmGCRS=[\\boldsymbolP][\\boldsymbolN][\\boldsymbolR][\\boldsymbolW]\\left(\\beginarrayl X \\\\ Y \\\\ Z \\endarray\\right)_\\mathrmITRS
⎝
⎛XYZ⎠
⎞GCRS=[P][N][R][W]⎝
⎛XYZ⎠
⎞ITRS
(
X
Y
Z
)
ITRS
=
[
W
]
−
1
[
R
]
−
1
[
N
]
−
1
[
P
]
−
1
(
X
Y
Z
)
G
C
R
S
\\left(\\beginarrayl X \\\\ Y \\\\ Z \\endarray\\right)_\\text ITRS =[\\boldsymbolW]^-1[\\boldsymbolR]^-1[\\boldsymbolN]^-1[\\boldsymbolP]^-1\\left(\\beginarrayl X \\\\ Y \\\\ Z \\endarray\\right)_\\mathrmGCRS
⎝
⎛XYZ⎠
⎞ITRS =[W]−1[R]−1[N]−1[P]−1⎝
⎛XYZ⎠
⎞GCRS
式中,
[
P
]
[\\boldsymbolP]
[P] 为岁差矩阵;
[
N
]
[\\boldsymbolN]
[N] 为章动矩阵;
[
R
]
[\\boldsymbolR]
[R] 为地球自转矩阵;
[
W
]
[\\boldsymbolW]
[W] 为极移矩阵。 考虑到 IGS 已完成了坐标转换工作, 在精密星历中直接给出了卫星质心在 ITRS 中的位置和速度, 而广播星历的精度有限, 允许采用一些近似的转换方法, 因此在下面的坐标转换中, 我们仍采用经典的转换方法与术语 (与 IS-GPS-200D 及 IS-GPS-705 中给出的方法基本 一致)。高精度的严格方法可参阅 IAU 的决议文件和空间大地测量学等参考资料。
（1）把 GCRS 转换至观测时刻
t
i
t_i
ti 的平天球坐标系
我们知道, GCRS 是参考时刻
t
0
=
J
2000.0
t_0=\\mathrmJ 2000.0
t0=J2000.0 时的平天球坐标系, 要把它转换为观测时刻
t
i
t_i
ti 时的平天球坐标系, 只要考虑
[
t
0
−
t
i
]
\\left[t_0-t_i\\right]
[t0−ti] 时间段内的岁差改正, 即乘上
[
P
]
−
1
[\\boldsymbolP]^-1
[P]−1 矩阵即可。
（2）把
t
i
t_i
ti 时的平天球坐标系转换为同一时刻的真天球坐标系
要把观测时刻
t
i
t_i
ti 时的平天球坐标系转换为真天球坐标系, 只需顾及该时刻的章动, 即只需乘上
[
N
]
−
1
[\\boldsymbolN]^-1
[N]−1 矩阵即可。
（3）把
t
i
t_i
ti 时的真天球坐标系转换为同一时刻的真地球坐标系
我们知道, 真天球坐标系
X
X
X 轴是指向该时刻的真春分点
γ
\\gamma
γ 的,而真地球坐标系的
X
X
X 轴 是指向起始子午线与赤道的交点,两者之间的夹角称为格林尼治真恒星时 GAST 。其计算 公式如下:
GAST
=
36
0
∘
2
4
h
(
U
T
1
+
6
h
41
m
50.54841
s
+
8640184.812866
s
⋅
t
+
0.093104
s
⋅
t
2
−
6.2
s
×
1
0
−
6
⋅
t
3
)
+
Δ
Ψ
cos
(
ε
ˉ
+
Δ
ε
)
\\beginaligned \\text GAST =& \\frac360^\\circ24^\\mathrmh\\left(\\mathrmUT 1+6 \\mathrm~h 41 \\mathrm~m 50.54841 \\mathrm~s+8640184.812866 \\mathrm~s \\cdot \\mathrmt+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) \\endaligned
GAST =24h360∘(UT1+6 h41 m50.54841 s+8640184.812866 s⋅t+0.093104 s⋅t2−6.2 s×10−6⋅t3)+ΔΨcos(εˉ+Δε)
式中,
t
t
t 为离
J
2000.0
\\mathrmJ 2000.0
J2000.0 的儒略世纪数;
ε
ˉ
\\bar\\varepsilon
εˉ 为仅顾及岁差时的黄赤交角,
ε
ˉ
=
2
3
∘
2
6
′
21.44
8
′
′
−
\\bar\\varepsilon=23^\\circ 26^\\prime 21.448^\\prime \\prime-
εˉ=23∘26′21.448′′−
46.81
5
′
′
⋅
t
−
0.0005
9
′
′
⋅
t
2
+
0.00181
3
′
′
⋅
t
3
;
Δ
ψ
46.815^\\prime \\prime \\cdot t-0.00059^\\prime \\prime \\cdot t^2+0.001813^\\prime \\prime \\cdot t^3 ; \\Delta \\psi
46.815′′⋅t−0.00059′′⋅t2+0.001813′′⋅t3;Δψ
46.81
5
′
′
⋅
t
−
0.0005
9
′
′
⋅
t
2
+
0.00181
3
′
′
⋅
t
3
;
Δ
Ψ
46.815^\\prime \\prime \\cdot t-0.00059^\\prime \\prime \\cdot t^2+0.001813^\\prime \\prime \\cdot t^3 ; \\Delta \\Psi
46.815′′⋅t−0.00059′′⋅t2+0.001813′′以上是关于ITRS 与 GCRS 之间的坐标转换的主要内容，如果未能解决你的问题，请参考以下文章