PyEphem: RA и разделительные вопросы - PullRequest
Новая
       14

PyEphem: RA и разделительные вопросы

0 голосов
/ 29 октября 2018

Я экспериментирую с библиотекой PyEphem для астрономии и пытаюсь воссоздать / понять некоторые из базовых вычислений.

Существует функция, называемая разделением, которая, учитывая две планеты и дату / время, вычисляет расстояние (угол) между этими двумя планетами относительно их проекции "x" на плоскость.

Итак, на 2018/1/1 для планет Меркурий и Марс имеем: 

import ephem
import math 

mercury = ephem.Mercury('2018/1/1')
mars = ephem.Mars('2018/1/1')

s1 = ephem.separation(mercury, mars)
print(math.degrees(s1))

Который возвращает 

33.792384499568264

Но если бы я хотел вычислить это без функции разделения, то вычисление было бы так же просто, как «прямое восхождение» 1 минус «правильное восхождение» другого: 

math.degrees(mercury.ra) - math.degrees(mars.ra)

Что возвращает 

35.114532008671574

Почему углы разные? Так как я не включаю широту и долготу наблюдателя, все расчеты должны быть геоцентрическими, согласно PyEphem.

Кто-нибудь знаком с вычислениями, проводимыми за PyEphem, или другой библиотекой со встроенным эфемеридой, которая может давать согласованные результаты для разделения?

1 Ответ

0 голосов
/ 31 октября 2018

Хороший вопрос! Вам захочется обратиться к учебнику по сферической тригонометрии. Точно так же, как вы не можете предсказать количество миль между двумя точками просто по их долготе - потому что они могут быть как на экваторе, так и далеко друг от друга, или находиться вблизи Северного полюса и находиться всего в нескольких футах друг от друга - поэтому вы не можете предсказать разделение между двумя объектами только по прямому Вознесению (= "долгота неба"), но нужно будет добавить в (a) склонения двух объектов и (b) использовать специальную формулу, которая вычисляет дугу между двумя положениями на сфере по маршруту большого круга.

...