Writeup Hack-A-Sat : Fiddlin' John Carson

Description du challenge : Your spacecraft has provided its Cartesian ICRF position (km) and velocity (km/s). What is its orbit (expressed as Keplerian elements)?

Le challenge est disponible à cette adresse

Connectons-nous à l'instance du challenge :

On va devoir à partir d'une position cartésienne et d'une vitesse retrouver les paramètres d'orbite de notre vaisseau spatial. Il est important de comprendre de quoi on parle quand on évoque les paramètres d'orbite. Heureusement pour toi, j'ai fait un petit cours que tu peux consuler juste ici :)

Voici un petit schema de ce que représente les différentes données qu'on a (le placement du vaisseau est arbitraire, il ne correspond pas à ses coordonnées réelles, flemme de faire un truc en 3D) :

Bon, en vrai, les calculs pour trouver les paramètres d'orbite sont plutôt compliqués, donc on va juste apprendre à utiliser la bibliothèque Python poliastro pour arriver à nos fins.

Dans un premier temps, on va convertir nos listes qui contiennent les coordonnées de position et de vitesse en objets avec des unités physiques appropriées pour faire des calculs avec.

Le u.km, c'est une unité de mesure qui provient de la bibliothèque astropy. En multipliant notre liste avec, on convertit chaque coordonnée en objet Quantity. Ainsi, les coordonnées x, y et z sont désormais traitées comme ayant pour unité le km ce qui va nous permettre de faire des calculs physiques avec. Pareil pour le u.s qui représente des secondes et donc u.km / u.s "transforme" notre liste de vitesse avec comme unité le km/s. Ensuite, et c'est là que la magie se fait, on a juste à appeler la méthode Orbit.from_vectors qui va s'occuper de faire tous les calculs pour nous afin qu'on puisse récupérer tout ce dont on a besoin.

Earth, c'est juste un objet qui contient les paramètres gravitationnels et géométriques de la Terre. Ces derniers sont indispensables pour le calcul mais pas la peine de rentrer dans les détails. On oublie pas aussi de lui spécifier notre temps (time) exacte qui représente le moment dans le temps où la position et la vitesse ont été mesurés. Et oui, les paramètres orbitaux peuvent être amenés à changer avec le temps en raison de divers perturbations gravitationnelles donc faut le spécifier ce temps. Une fois, l'appel à Orbit.from_vectors fait, on a plus qu'à récupérer nos 6 paramètres d'orbite :

Avec la méthode to_value d'astropy, on peut convertir directement une valeur avec l'unité de son choix. En l'occurrence, comme les angles sortent en radians, on s'en sert pour les convertir en degrés. On print tout ça, et y a plus qu'à remplir avec les bonnes valeurs :

Et on a notre flag :)