Optimisation d'un algorithme statistique itératif de reconstruction tomographique en utilisant les coordonnées cylindriques

La tomographie axiale calculée par ordinateur utilise ordinairement des techniques de analytiques pour la reconstruction d'image, mais les méthodes de reconstruction statistiques itératives peuvent être une bonne alternative à ces techniques conventionnelles. Le principal problème de ces méthodes est qu'elles nécessitent une amélioration de leurs performances en terme de temps de calcul et d'espace mémoire requis pour pouvoir être plus couramment utilisées. Ce projet a pour but d'optimiser les performances globales des techniques itératives de reconstruction en 3D en utilisant une représentation alternative de l'objet ainsi que d'étudier l'impacte d'une telle représentation sur la qualité des images reconstruites. L'utilisation de la représentation de l'objet en coordonnées cylindriques a un impact significatif sur les performances globales des techniques itératives car elle permet une réduction considérable de la mémoire utilisée mais ce modèle est encore très lent, son effet sur la qualité de l'image méconnu et ses possibilités d'amélioration très peu étudiées. Ce projet se concentre sur les améliorations possibles avec cette représentation cylindrique permettant l'augmentation de la vitesse de calcul en utilisant la forme particulière de la matrice de projection en 3D. Ce projet consiste dans un premier temps à développer une nouvelle technique de parallélisation efficace de la méthode qui utilise la structure bloc-circulante de la matrice de projection, mais aussi à améliorer la vitesse de convergence grâce à la création de préconditionneurs utilisant la propriété diagonale-dominante de la matrice de projection bloc-diagonalisée dans le domaine de Fourier. La deuxième partie de ce projet consiste en une étude poussée de la qualité des images résultant de l'utilisation de cette représentation cylindrique de l'objet, une application d'un modèle polychromatique afin de réduire les artéfacts métalliques et la mise en oeuvre de la méthode pour la reconstruction de données réelles.