Мне дан список дат в формате UTC, все часы переведены на 00:00.
Я хотел бы определить, произошло ли (лунное) затмение в данный день (то есть за последние 24 часа)
Учитывая фрагмент Python
from sykfield.api import load
eph = load('de421.bsp')
def eclipticangle(t):
moon, earth = eph['moon'], eph['earth']
e = earth.at(t)
x, y, _ = e.observe(moon).apparent().ecliptic_latlon()
return x.degrees
Я предполагаю, что можно определить, произошло ли затмение в течение 24 часов времени t на
- Проверка того, что первый угол достаточно близок к 180 (легко)
- Проверка, достаточно ли вторая степень близка к 0 (не так уж и сложно?)
Теперь, что касается ответа в комментарии, не так уж и тривиально решить вторую проблему, просто проверив, близок ли угол к 0.
Поэтому у меня вопрос
Может ли кто-нибудь предоставить функцию для определения того, произошло ли лунное затмение в данный день t?
Изменить. Этот вопрос был отредактирован с учетом отзывов Брэндона Роудса, оставленных в комментариях ниже.