Рассчитайте время, когда звезда находится выше высоты 30°

Чтобы найти лучшее время наблюдения за объектом, я хотел бы рассчитать время, когда он находится на 30 ° или более над горизонтом. Местного бокового времени будет достаточно.

Чтобы включить это в мою программу, мне нужна формула.

Пример: 4 июня Юпитер имеет координаты RA= 9h 19m 28.0s Dekl= 16° 32' 0"

Он восходит в 10:32 и заходит в 00:05.

После подъема, когда он снова на высоте 30°, и после прохождения, когда он снова на высоте 30°?


Я нашел эту формулу на http://www.stjarnhimlen.se/comp/riset.html . Хотя это для солнца, кажется, это то, что я ищу.

потому что ( ЛХА ) знак равно грех ( час ) грех ( лат ) × грех ( Децл ) потому что ( лат ) × потому что ( Децл )

Применительно к образцу с учетом широты 45° я получаю.

Пример расчета
Это правильный подход?

Проверьте catserver.ing.iac.es/staralt .
Спасибо, но мне нужна формула, а не инструмент.
Земля вращается на 15 градусов каждый час.
15 градусов в час: это одна часть формулы.
Я думаю, вам лучше объяснить, что вы знаете об объекте (объектах), который хотите наблюдать. Не совсем понятно, что именно вы не понимаете. Учитывая координаты объекта и дату, время, широту и долготу, это прямой перевод координат.
@ andy256 andy256 Спасибо за ваш комментарий. Я добавил пример, чтобы уточнить мой вопрос. Тем временем я нашел формулу, но я не уверен, что делаю это правильно.
Простой способ проверить свой ответ — использовать Stellarium или ssd.jpl.nasa.gov/?horizons Вы также можете посмотреть мой astronomy.stackexchange.com/questions/8390/… (формула закрытой формы), но уродливый.
@barrycarter Спасибо. Сравнение со Stellarium не доказывает, что это правильный подход. Однако, если здесь никто не опубликует ответ и не подтвердит формулу, я проверю все соответствующие сценарии и надеюсь на лучшее. Упомянутая вами формула, вероятно, сработает, но как нематематик я не могу поставить d1 в левую часть уравнения. Далее описанная задача требует двух результатов. Один, когда объект поднимается, и один, когда он собирается сесть. В дополнение к этому необходимо обрабатывать два особых случая: 1) Объект никогда не становится ниже 30 градусов 2) Объект никогда не поднимается выше 30 градусов.
idlastro.gsfc.nasa.gov/ftp/pro/astro/hadec2altaz.pro подтверждает, что ваш ответ правильный. Два особых случая (всегда выше 30 и никогда не выше 30) будут иметь место, когда правая часть больше 1 или меньше -1, значений, которых косинус никогда не может достичь (по крайней мере, для действительных чисел).
Возможно, вы захотите изучить PyEphem (если вы разрабатываете на Python) или его базовую библиотеку libastro (если на C). rhodesmill.org/pyephem

Ответы (3)

Да, это правильный подход. час в уравнении — это высота над горизонтом объекта, на котором вы считаете, что он восходит или заходит. Обычно это не ноль из-за атмосферной рефракции и, в случае Солнца или Луны, из-за их конечных диаметров. В вашем случае объект «поднимается», когда он поднимается выше час знак равно 30 и «устанавливается», когда он падает ниже этой высоты.

Если | потому что ( л ЧАС А ) | > 1 , решения нет, потому что ваш объект никогда не пересекает час знак равно 30 линия. л ЧАС А это местный часовой угол , и вы можете найти местное звездное время θ с использованием

л ЧАС А знак равно θ α ,

куда α прямое восхождение вашего объекта.

Ответ есть на сайте stjarnhimlen.se, а также на stargazing.net .

Теперь мы можем вычислить высоту Солнца над горизонтом:

sin(h) = sin(lat) * sin(Decl) + cos(lat) * cos(Decl) * cos(LHA)
LHA = LST - RA

ч=высота=30 , LHA = местный часовой угол, lat = ваша широта на Земле, Decl = склонение объекта, Ra = прямое восхождение объекта и LST = местное стандартное время. Просто нужно решить для LST.

Вы не говорите, какой язык программирования вы используете. Если это Python, или если вы можете связать библиотеки Python с ним, то PyEphem предоставит все, что вам нужно.

http://rhodesmill.org/pyephem/