Скайфилд в поисках внешних планет

Я пытаюсь получить местоположение внешних планет, и я просто хочу убедиться, что у меня правильный код.

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 правильным способом сделать это?


person brian lee    schedule 18.05.2017    source источник


Ответы (1)


Я не очень хорошо разбираюсь в вызовах Skyfield, но когда я хочу проверить расчеты, которые я получаю, используя его предшественника, pyephem, я обычно иду на веб-сайт JPL Horizon. Довольно понятно, как настроить те же параметры в его интерфейсе. Результаты, которые вы получите, обычно будут немного отличаться, но достаточно близки, чтобы убедиться, что вы на правильном пути.

URL-адрес: https://ssd.jpl.nasa.gov/horizons.cgi.

Использование барицентра Сатурна совершенно нормально. Результат, который вы получите, если вы использовали центр планеты, будет отличаться только в угловых секундах, если это так. Барицентр системы Сатурн (планета + луны), вероятно, в любом случае расположен внутри планеты, поскольку ее масса намного, намного больше, чем ее луны.

person sizzzzlerz    schedule 19.05.2017