Cálculo y mapeo de la mejor visibilidad de ocurrencia celestial

Recientemente encontré un sitio web fantástico en-the-sky.org , desde aquí podemos simplemente recoger cualquier evento astronómico y verificar su visibilidad en cualquier parte del Globo. Para ocultaciones, por ejemplo, tenemos mapas muy buenos como este a continuación:

Hay muchas conjunciones, ocultaciones, eclipses, etc., descritas con bastante detalle. Creo que es el mejor sitio web hasta ahora para las personas que hacen calendarios Astro.

El problema es que si revisé alguna conjunción, es decir, Mercurio con Marte el 19 de agosto de 2021, dediqué mucho tiempo a verificar qué ubicación alrededor del Globo es la mejor para ver este evento. Tomé en consideración la depresión solar, no inferior a 6 grados (fin del crepúsculo civil), y una altitud de estos 2 cuerpos sobre el horizonte (no inferior a 3 grados por extinción atmosférica). Finalmente, mis mapas dibujados a mano se ven así:

trazado en el mapa de día y noche en timeanddate.com

croquis trazado en Google Maps

Circunstancias verificadas con stellarium-web.org

Entonces, básicamente, es bastante trabajo verificar qué ubicación es la más apropiada para observar el momento ideal de la conjunción (me refiero al momento de mayor acercamiento de estos 2 planetas, que dura tal vez ... 30 minutos más o menos). Para concluir, me gustaría preguntar sobre la posibilidad de calcular automáticamente cosas como esta. Si por ejemplo, pusiera en una columna la depresión solar mínima, en otra columna la altitud mínima de estos planetas, ¿alguna herramienta me traería el mapa como se presenta arriba?

¿Es posible, es decir, con Python u otro (solo soy básico con suave)?

¿O conoce los servicios web, que hacen el cálculo como este por usted?

Fondo alguna consulta aquí:

Como calcular conjunciones de 2 planetas

pero se aplica para calcular las fechas y circunstancias geométricas de las conjunciones. Me gustaría tener algo en el mapa como se presenta arriba.