スカイフィールド外の惑星を見つける

Question

私は外の惑星の位置を取得しようとしています。私はコードが正しいことを確認したいのですが。

lon = 45.138967 
lat = -93.342730 

ts = load.timescale() 
t = ts.now() 
print(t.utc_jpl()) 

planets = load('de422.bsp') 

body = planets['SATURN_BARYCENTER'] 
earth = planets['Earth'] 

myLocation = earth + Topos(lat, lon) 
apparent = myLocation.at(t).observe(body).apparent() 

alt, az, distance = apparent.altaz() 
print(alt.degrees) 
print(az.degrees) 
print(distance) 

これを行う正しい方法はSATURN_BARYCENTERですか？

Answer 1

スカイフィールドはすべてうまくいきませんが、以前のpyephemを使って得た計算結果を確認したいときは、一般的にJPLのHorizo​​n Webサイトに行きます。そのインターフェイスで同じパラメータを設定する方法はかなり自明です。あなたが得る結果は、一般的にわずかに異なりますが、あなたが正しい方向にいることを確認するのに十分近いでしょう。

urlは土星の重心を使用してhttps://ssd.jpl.nasa.gov/horizons.cgi

は完全に罰金です。あなたが惑星の中心を使用した場合に得られる結果は、円弧でのみ異なる場合があります。土星系（惑星+月）の重心は、質量が月よりもずっと大きいので、おそらく惑星の内側に位置するでしょう。