2行軌道要素形式のソースを表示

{{脚注の不足|date=2017-01-03}}
'''2行軌道要素形式'''（にぎょうきどうようそけいしき、{{lang-en-short|Two-line elements}} ： '''TLE'''）は、[[アメリカ航空宇宙局]] (NASA) と[[北アメリカ航空宇宙防衛司令部]] (NORAD) が現在でも使用している、[[人工衛星]]の[[測地系|地心座標系]]における[[人工衛星の軌道要素#ケプラー軌道要素|ケプラー軌道要素]]の[[テキストファイル|テキスト形式]]のフォーマットである<ref name="c1">[http://celestrak.com/NORAD/documentation/tle-fmt.asp CelesTrak: NORAD Two-Line Element Set Format]{{リンク切れ|date=2024年4月}}</ref><ref>[http://celestrak.com/columns/v04n03/ Celestrak Two-Line Elements FAQ]</ref>。元来は初期のコンピューターの80桁の[[パンチカード]]用としてデザインされたフォーマットであるが、さまざまな分野で非常に普及しており、また他のいかなるフォーマットと比べても遜色なく働くことから、現在でも使用されている。

応用分野や対象となる軌道にもよるが、更新から30日以上経過した2行軌道要素形式を用いて計算された値は、信頼性に欠ける可能性がある。衛星の軌道上の位置は、2行軌道要素形式から、SGP、[[SGP4]]、[[SDP4]]、SGP8、SDP8 の各アルゴリズムを用いて計算される<ref name="Spacetrack Report No. 3">{{Cite report |title=Spacetrack Report No. 3 - Models for Propagation of NORAD Element Sets |date=1988-12-31 |url=http://www.celestrak.com/NORAD/documentation/spacetrk.pdf |publisher=The U.S. Department of Defense, Aerospace Defense Center, Peterson AFB, CO. |format=PDF |authors=Felix R. Hoots, Ronald L. Roehrich, TS Kelso |website=CelesTrak |accessdate=2022-04-27}}</ref><ref>Explanatory Supplement to the Astronomical Almanac. 1992. K. P. Seidelmann, Ed., University Science Books, Mill Valley, California.</ref>。SGP4 を使用した場合の精度は、位置に関して典型的には誤差 1 km である。例えば 300 km 離れた位置からは、これは最大 0.2° の観測誤差を引き起こす。

== フォーマットの詳細 ==
2行軌道要素形式は、テキスト形式の1行69文字の2行 (Line 1 と Line 2) から成る。使用可能な文字は、英大文字 <code>A</code>–<code>Z</code>、数字 <code>0</code>–<code>9</code>、小数点<code>.</code>、空白<code> </code>および <code>+</code> と <code>-</code> の符号のみである。

実際の利用においては、分かり易いように、Line 1 の前に Line 0 として24文字以内の衛星名を付加することが広く行われている。衛星名として使用可能な文字は Line 1 と Line 2 で利用可能な文字よりもやや自由度が大きく、少なくとも英大文字 <code>A</code>–<code>Z</code>、数字 <code>0</code>–<code>9</code>、丸括弧 <code>(</code> と <code>)</code>、角括弧 <code>[</code> と <code>]</code>、<code>+</code> と <code>-</code> の符号、および空白が利用可能である（その他の[[ASCII]]文字も利用されている可能性がある）。衛星名の24文字という制限は、NORAD および NASA で使用されている衛星名と整合性を保つための慣例である（Line 0 を付加して3行の形式とした場合でも、2行軌道要素形式と呼ばれる）。以下の説明は Line 0 を付加した形式で行う。

=== 一般フォーマット ===
<pre>
AAAAAAAAAAAAAAAAAAAAAAAA
1 BBBBBC DDEEEFFF GGHHH.HHHHHHHH +.IIIIIIII +JJJJJ-J +KKKKK-K L MMMMN
2 BBBBB PPP.PPPP QQQ.QQQQ RRRRRRR SSS.SSSS TTT.TTTT UU.UUUUUUUUVVVVVW
</pre>

<code>AAAAAAAAAAAAAAAAAAAAAAAA</code> が24桁の衛星名 (Line 0) である。
また、下記以外の文字は、次で説明するためのシンボルである。

凡例（一般フォーマットの各項目共通）
* <code>1</code>、<code>2</code>、<code>.</code> ： これらはそれぞれ各桁番号および小数点を表す。これ以外の文字は来ない。
* <code>+</code> ： 続く数値の符号が来ることを表し、<code>+</code>、<code>-</code>または空白が来る。
* <code>-</code> ： 1桁の数値が続く場合は、その数値の符号であり、<code>+</code> か <code>-</code> が来る。その場合、この符号付の1桁の整数は、先行する数値の10を底とする指数部である。この2桁は空白になる場合がある。
* その他 : シンボル

{| class="wikitable"
|+ Line 1
|-
! 桁
! 桁数
! 詳細
! シンボル
|-
||1||1||要素データ行番号。 Line 1 では必ず <code>1</code>||1
|-
||3–7||5||衛星カタログ番号 (NORADカタログナンバー)||BBBBB
|-
||8||1||軍事機密種別（<code>S</code>：秘匿、<code>U</code>：公開、のどちらかの文字が来る。）||C
|-
||10–11||2||国際衛星識別符号 (打上げ年のラスト2桁)||DD
|-
||12–14||3||国際衛星識別符号 (その年の打上げの通算番号)||EEE
|-
||15–17||3||国際衛星識別符号（その打上げによる飛行体の通番）飛行体が1つの場合は3桁とも空白。飛行体が複数の場合は、各々の飛行体を区別する目的で左詰で <code>A</code>、<code>B</code>、<code>C</code>、…、<code>Z</code>、<code>AA</code>、…、<code>AZ</code>、<code>BA</code>、…、<code>ZZ</code>、<code>AAA</code>、… の要領で表す||FFF
|-
||19–20||2||この軌道要素の元期 (年のラスト2桁)||GG
|-
||21–32||12||元期 (続き) (その年の通日3桁、および該日における時刻を表す9桁の小数)。先行する <code>GG</code> が示す年の1月1日0000 UTC を1.0と数え、それからの経過時間を日を単位として加算した値。最小位の桁は0.864msに相当する。UTC を使用する場合は[[閏秒]]に注意する必要がある。閏秒の無い[[国際原子時|TAI]]、[[地球時|TT]]、[[UNIX時間]]、[[グローバル・ポジショニング・システム#GPS時刻|GPSタイム]]を用いた方が良い||HHH.HHHHHHHH
|-
||34–43||10||平均運動の1次微分値を2で割った値。単位は [回転/day{{sup|2}}]。主に地球の重力ポテンシャルの球対称からのずれにより発生する二体問題への摂動の影響を、平均運動について補正する。SGPで使用されるが、SGP4では使用されていない。||+.IIIIIIII
|-
||45–52||8||平均運動の2次微分値を6で割った値（仮数部は先頭に小数点があるものと見なされる）。単位は [回転/day{{sup|3}}]。<br>実際には使用されておらず、NASA、NORADなどの公表値は <code>00000-0</code>にセットされている。||+JJJJJ-J
|-
||54–61||8||B* (B STAR) 抗力項（仮数部は先頭に小数点があるものと見なされる）。大気による抗力が無視できない場合に用いる補正係数。SGP4で使用されるが、SGPでは使用されていない。||+KKKKK-K
|-
||63||1||この軌道要素を算出した軌道モデル (Ephemeris) の種別。<code>0</code>:情報なし、<code>1</code>:SGP、<code>2</code>:SGP4、<code>3</code>:SDP4、<code>4</code>:SGP8、<code>5</code>:SDP8。この情報は通常は使用されず、NASA、NORADなどの公表値は <code>0</code> となっている。||L
|-
||65–68||4||軌道要素通番。軌道要素情報が更新されるごとに1ずつ加算される。||MMMM
|-
||69||1||Line 1 チェックサム。<br>68桁までの各桁の数字を10進1桁の数と見なして単純に加算した総和を10で割った余り。ただし、数字以外の文字、空白、小数点、<code>+</code> は無視し、<code>-</code> は <code>1</code> と見なす。||N
|-
|}

{| class="wikitable"
|+ Line 2
|-
! 桁
! 桁数
! 詳細
! シンボル
|-
||1||1||要素データ行番号。 Line 2 では必ず <code>2</code>||2
|-
||3–7||5||衛星カタログ番号（Line 1 と同一）||BBBBB
|-
||9–16||8||軌道傾斜角。単位は [度]||PPP.PPPP
|-
||18–25||8||昇交点の赤経。単位は [度]||QQQ.QQQQ
|-
||27–33||7||離心率。先頭に小数点があるものと見なされる。||RRRRRRR
|-
||35–42||8||近地点引数。単位は [度]||SSS.SSSS
|-
||44–51||8||平均近点角（Mean Anomaly）。単位は [度]||TTT.TTTT
|-
||53–63||11||平均運動 (Mean Motion)。単位は [回転/day]||UU.UUUUUUUU
|-
||64–68||5||元期における通算周回数。軌道上に打上げられてから最初に昇交点を通過するまでを 0 周目として、昇交点を通過するごとに 1ずつ増える値。||VVVVV
|-
||69||1||Line 2 チェックサム。計算方法は Line 1 と同じ||W
|-
|}

=== 実際の形式の例 ===
以下に具体例をあげる<ref name="c1" />。

<pre>
ISS (ZARYA)
1 25544U 98067A   22095.91869325  .00012930  00000-0  23502-3 0  9991
2 25544  51.6452 334.5328 0004408 351.0413  99.6998 15.49890618333972

MIDORI (ADEOS)          
1 24277U 96046A   09116.47337938 -.00000023  00000-0  73445-5 0   432
2 24277  98.3597  83.2073 0002090  64.7512 295.3886 14.28595439661547

ORBCOMM FM08 [+]        
1 25112U 97084A   09116.51259343  .00000203  00000-0  12112-3 0  2154
2 25112  45.0199 241.1109 0010042 194.4473 165.6089 14.34380830592834
</pre>

== 軌道要素間の関係 ==
軌道要素のうち、2行軌道要素形式には次の要素が直接含まれている。
:<math style="vertical-align:-50%;">t_{0}</math> : [[元期]] ([[:en:Epoch|Epoch]])
:<math>i_0</math> : 元期における [[軌道傾斜角]] ([[:en:inclination|inclination]])
:<math>\Omega_0</math> : 元期における [[昇交点赤経]] ([[:en:longitude of the ascending node|longitude of the ascending node]])
:<math>e_0</math> : 元期における[[軌道離心率]] ([[:en:Orbital eccentricity|Orbital eccentricity]])
:<math>\omega_0</math> : 元期における[[近点引数]] ([[:en:Argument of periapsis|Argument of periapsis]])
:<math>M_{0}</math> : 元期における[[平均近点角]] ([[:en:mean anomaly|mean anomaly]])
:<math>n_0</math> : 元期における[[平均運動]] ([[:en:Mean Motion|Mean Motion]])

これらの要素は、厳密な[[二体問題]]が成立する場合は定数であるが、次のような原因による[[摂動]]で変動する。
* 地球の重力ポテンシャルの球対称からのずれ
* 地球の大気による抗力
* 太陽光圧、太陽風圧
* 地球以外の天体(月、太陽、他の惑星など)の重力

2行軌道要素形式には、主に平均運動の時間変動を補正する、次のような情報が含まれている。
:<math>\dot{n}_0 / 2</math> : 元期における平均運動の時間についての1次微分を2で割った値。
:<math>\ddot{n}_0 / 6</math> : 元期における平均運動の時間についての2次微分を6で割った値。
:{{Math|''B*''}} : 地球大気の抗力による平均運動への影響を表すパラメーター

<math>\dot{n}_0/2</math> および <math>\ddot{n}_0/6</math> はSGPで用いられているが、SGP4では用いられていない。逆に {{Math|''B*''}} はSGP4で用いられているが、SGPでは用いられていない<ref name="Spacetrack Report No. 3"/>。

SGPやSGP4では、地球の重力ポテンシャルのモデルなどを用いて、2行軌道要素形式に含まれる元期 ''t<sub>0</sub>'' における軌道要素の値から、任意の時刻 ''t'' における次の軌道要素の値を予測している<ref name="Spacetrack Report No. 3"/>。
:<math>i</math> : 時刻 <math>t</math>における軌道傾斜角
:<math>\Omega</math> : 時刻 <math>t</math>における昇交点赤経
:<math>e</math> : 時刻 <math>t</math>における軌道離心率
:<math>\omega</math> : 時刻 <math>t</math>における近点引数
:<math>M</math> : 時刻 <math>t</math>における平均近点角
:<math>n</math> : 時刻 <math>t</math>おける平均運動

以下の軌道要素は、2行軌道要素形式には元期 ''t<sub>0</sub>'' における値が直接は含まれていないが、任意の時刻 ''t'' における平均運動 ''n'' と離心率 ''e'' から、''t'' における値を計算可能である。
:<math>a</math> : 時刻 <math>t</math>における[[軌道長半径]]([[:en:Semimajor axis|Semimajor axis]])
:<math>l</math> : 時刻 <math>t</math>における[[円錐曲線#離心率による分類|半通径]](または半直弦)([[:en:Conic_section#Conic_parameters|semi-latus rectum]])
:<math>q</math> : 時刻 <math>t</math>における[[近点距離]] ([[:en:Periapsis|Periapsis]])
:<math>Q</math> : 時刻 <math>t</math>における[[遠点距離]] ([[:en:Apapsis|Ap(o)apsis]])
:<math>T</math> : 時刻 <math>t</math>における[[公転周期|軌道周期]]([[:en:Orbital period|Orbital period]])

これらには、次のような関係がある。
:<math>a = \sqrt[3]{\frac{ GM_E }{ n^2 }}</math>
:<math>l = a ( 1 - e^2 )</math>
:<math>Q = a(1+e)</math>
:<math>q = a(1-e)</math>
:<math>T = \frac{ 2\pi }{ n }</math>

ここで {{Math|''GM<sub>E</sub>''}} は[[万有引力定数#天体の質量との積|地心重力定数]]([[万有引力定数]]と[[地球質量]]の積)であり{{Math|''GM<sub>E</sub>''}} = {{Val|3.986004418|(8)|e=14|u=m{{sup|3}} s{{sup-|2}}}} である。

地球中心と衛星間の距離を ''r'' とすると、 ''r'' は[[真近点角]] {{Math|''&nu;''}} の関数として次の形に表すことができる。
:<math>r = \frac{ l }{ 1 + e \cos \nu }</math>

また、''r'' は、[[離心近点角]] ''E'' の関数として次の形に表すことができる。
:<math>r =a(1-e\cos E)</math>

真近点角 {{Math|''&nu;''}} と離心近点角 ''E'' の関係は、次のようになる。
:<math>\tan\frac{\nu}{2} = \sqrt{\frac{1+e}{1-e}} \tan\frac{E}{2}</math>

離心近点角 ''E'' と平均近点角 ''M'' の関係は次の式([[ケプラー方程式]])で表される。
:<math>M = E -e \sin E</math>

平均近点角 ''M'' は次の式で表される。
:<math>M = M_0 + {n_0}(t - t_0) + \delta M</math>

{{Math|''&delta;M''}} は摂動による ''M'' の変動を表す項であり、摂動がなければ {{Math|''&delta;M''}}  = 0 である。

時刻 ''t'' における 平均運動 ''n'' は 、時刻 ''t'' における平均近点角 ''M'' の時間微分である。
:<math>n = \dot M = n_0 + \delta {\dot M}</math>

== 脚注 ==
{{脚注ヘルプ}}
{{reflist}}

== 外部リンク ==
* [http://celestrak.com/ CelesTrak]
* [https://www.space-track.org/ Space track(United States Space Surveillance Network)]
*  [http://spaceflight.nasa.gov/realdata/sightings/SSapplications/Post/JavaSSOP/SSOP_Help/tle_def.html Human Space Flight (HSF) - Realtime Data - NASA]
* [http://www.lizard-tail.com/isana/tle/misc/what_is_tle きどうようそのひみつ]

{{軌道}}

{{DEFAULTSORT:2きようきとうようそけいしき}}
[[Category:力学]]
[[Category:軌道]]