TLE(Two-line elements)からケプラーの軌道要素を読み取る

ALOS-2の軌道
www.n2yo.com

Two Line Element Set (TLE):
1 39766U 14029A 22113.21647048 -.00000519 00000-0 -62686-4 0 9997
2 39766 97.9191 211.1359 0001719 96.6291 263.5144 14.79467396427466

TLEの2行目からパラメータが読み取れる
フォーマットは以下。
2 カタログ番号 INC RAAN ECC AOP MA Mean motion(revolution per day)

すなわち
INC=97.9191
RAAN=211.1359
ECC=0001719
AOP=96.6291
MA=263.5144

TLEには長半径(SMA: Semi-major axis)が見当たらない。
Mean motion(revolution per day)から換算できるとのこと。
space.stackexchange.com

function a=mm_to_sma(n)
mu = 3.986004418e14;
% n = 15.5918272;
a = (mu^(1/3))/(2*n*pi/86400)^(2/3);
a = a/1000; %km
format long g;
end

MM=14.79467396427466(revolution/day)より、

>> mm_to_sma(14.79467396427466)
ans =
          7009.14280972952

SMA=7009.14280972952[km]と求まる。