В документах Pyephem есть пример, как вычислить созвездие, в котором планета находится в определенное время.
>>> m = ephem.Moon('1980/6/1')
>>> print(ephem.constellation(m))
('Sgr', 'Sagittarius')
Но как рассчитать время, когда планета сменит созвездие?
Поскольку планеты могут менять направление и скорость по небу, невозможно заранее предсказать, когда именно планета пересечет границу с другим созвездием. Простейшим подходом было бы сделать шаг вперед во времени на один день (или больше? или меньше?), пока выход constellation() не изменится для планеты, а затем вернуться назад и точно определить, в какое время дня произошел переход.
