Ce script dessine l’aspect du disque lunaire pour une date donnée, ainsi que sa position relativement au Soleil et à la Terre, pour expliquer l’évolution de la fraction éclairée visible.
Ce script Python s’appuie sur le module Skyfield (éphéméride astronomiques).
Le script calcule l’angle de phase (angle entre la Terre et le Soleil, vus depuis la Lune), et utilise une relation simple pour en déduire la fraction éclairée du disque lunaire (le pourcentage affiché) et dessiner une ellipse pour simuler le terminateur.
Pour le graphique de droite, on récupère les coordonnées géocentriques de la Lune et du Soleil, pour tracer un graphe polaire.
Le script ne trace qu’une seule date, et le fichier image généré est numéroté dans l’ordre chronologique. Ainsi, après avoir généré une série d’images, il est facile d’en faire un GIF animé, voire une vidéo.
Le script python compressé :
merci pour le partage de vos codes, cela m’a permis d’intégrer des visuels de skyfield dans home assistant