Skyfield が計算する長半径と TLE のものを比較したかったのです。
from skyfield.elementslib import osculating_elements_of as OE
from skyfield.api import EarthSatellite, load, utc
ts = load.timescale()
sat = EarthSatellite(TLE1, TLE2)
date = sat.epoch.utc
test = sat.at(ts.utc(*date))
print(OE(test).semi_major_axis.km)
print(sat_df['semimajor_axis_km'].values[0])
出力:
7087.916048058872
7080.642
それらが同じではない、または少なくとも近いのはなぜですか?