Fit d’une image¶
Le fichier data/double_star.txt contient une image de 64 par 64 pixels d’une étoile double. L’objectif de cette partie est d’ajuster cette image par la somme de deux Gaussiennes afin de déterminer la distance entre les étoiles.
Charger et afficher le fichier à l’aide de la fonction imshow
Définir une fonction de fit et tracer une image qui ressemble à celle ci. On pourra utiliser la fonction
np.meshgrid.
Il n’est pas possible de travailler directement sur des tableaux 2D pour le fit : l’ensemble des pixels de l’image doit être sous forme d’un tableau (taille \(N\)) ainsi que l’ensemble de coordonnées \((x,y)\) des ces points (tailles \(N\times 2\)).
Ecrire une fonction qui s’adapte au contraite du fit et effectuer le fit.
Quelle est la distance et son incertitude (en pixel) entre les deux étoiles ?