1

特定の時点でのオブジェクトの緯度、経度、および地球までの距離 (ECEF) があるとします。J2000座標を取得する必要があります。私がやっていることは

import skyfield.api
planets = skyfield.api.load('de421.bsp')
earth = planets['earth']
ts = skyfield.api.load.timescale()
sat = earth + skyfield.api.Topos(latitude_degrees=latitude, longitude_degrees=longitude)
obs = earth.at(ts.utc(time)).observe(sat)
x, y, z = obs.km

次に、それを正規化[x, y, z]し、地球までの距離を掛けます。関連する 2 つの質問があります。

  1. そうするのは正しいですか?
  2. 「de421.bsp」をロードせずに高速な方法はありますか?
4

1 に答える 1