Version imprimable
Bonjour,
Je débute sur Python et je travaille sur la transformée de Radon dans la reconstruction d'image.
Il y a un code Python de la transformée inverse sur github, mais je n'arrive pas à le comprendre dans son intégralité, est-ce que quelqu'un serait capable de m'expliquer les différentes étapes ?
Principalement des lignes : 275 à 282 ; 287 à 298
Le code.
Merci d'avance.
-Antoine
Salut,
L'algorithme est décrit ici - comme mentionné dans le source. Si vous avez des difficultés à le comprendre, vous avez le forum algo. (où se nichent quelques experts en ce type d'algo.). Impossible de l'y retrouver dans sa traduction en Python/numpy sans l'avoir compris.
- W
Pour calculer une transformée de Fourier rapidement (et la transformée inverse aussi), on utilise l'algo 'FFT' comme Fast Fourier Transform. Pour fonctionner cet algorithme prend en entrée des "tenseurs" qui doivent impérativement avoir des dimensions multiples de 2. Si l'image n'est pas un multiple de 2, alors on ajoute des pixels noir sur les bords de façon a avoir la bonne taille. C'est les lignes 275 a 277.
La transformee de radon inverse peut etre effectuee par transformee de fourier inverse (voir wikipedia et ...). Sauf que on a rajoute des 0 sur les bords de l'image donc il faut compenser ca. J'ai bien l'impression que c'est le role du filtre ``fourier_filter``. Bref ca c'est les lignes 280 a 282.
Pour les lignes 287 a 298, c'est la reconstruction de l'objet en tant que tel. Les etapes d'avant permetent de recuperer le profil de plusieurs vues en coupe. Et bien a cette etape on combine les vues en coupe pour reconstituer une seule image.
Je ne suis pas certain de ce que j'avance: je n'ai jamais vraiment maitrise la transformee de Radon puis je ne me suis pas penche des heures dans le code source.