Humboldt-Univ., Diss.--Berlin, 2003.

Ces travaux de thèse sont motivés par le développement de techniques d’imagerie alternatives pour le diagnostic précoce du cancer du sein. Parmi celles-ci, l’imagerie photoacoustique couple potentiellement les avantages de deux modalités d’imagerie non-invasives, à savoir la quantification de contrastes physiologiques du fait de l’excitation optique et la haute résolution du fait d’un sondage acoustique.Le but de ces travaux est de proposer une modélisation multiondes du phénomène photoacoustique, et d’incorporer ce modèle dans un algorithme de reconstruction efficace pour résoudre le problème inverse. Celui-ci se rapporte à la reconstruction de cartes de propriétés physiques (optique et/ou acoustiques) de l’intérieur du sein. La Méthode des Eléments Finis (MEF) a été retenue pour résoudre l’équation de propagation optique. Pour la résolution de l’équation de propagation acoustique, une méthode semi-analytique, basée sur des calculs par transformées de Fourier (méthod k-space), a été choisie. Pour la résolution du problème inverse, deux approches ont été étudiées : i) un sondage passif, permettant de remonter à la distribution de pression initiale, à l’aide de la méthode de retournement temporel ; ii) un sondage actif, où l’on interroge le milieu sélectivement sous différentes excitations, permettant de remonter quantitativement aux propriétés optiques du milieu. On appelle cette dernière approche Tomographie PhotoAcoustique Quantitative (TPAQ). Une étude spécifique sur le protocole d’illumination/détection a été conduite, prenant également en compte les contraintes expérimentales. / The present work was motivated by the development of alternative imaging techniques for breast cancer early diagnosis, that is photoacoustic imaging, which potentially couples the merits of optical imaging and ultrasound imaging, that is high optical functional contrasts brought by optical probing and high spatial resolution by ultrasound detection. Our work aims at modeling the photoacoustic multiwave phenomenon and incorporate it in an efficient reconstruction algorithm to solve the inverse problem. The inverse problem consists in the recovery of interior maps of physical properties of the breast. The forward model couples optical and acoustic propagations. The Finite Element Method (FEM) was chosen for solving the optical propagation equation, while a semi-analytical method based on Fourier transforms calculations (k-space method) was preferred for solving the acoustic propagation equation. For the inverse model, time reversal method was adopted to reconstruct the initial pressure distribution, an active approach of the inverse problem was also achieved, which decoupled the optical properties from measured photoacoustic pressure, this approach is called quantitative photoacoustic tomography (QPAT), in this approach, illumination/detection protocol was studied, and the experimental set up is also take into consideration. In the last step, photoacoustic pressure measurements obtained from experiment and simulation are studied and compared.

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal.

L'imagerie cone-beam computed tomography (CBCT) est une méthodologie de contrôle non destructif permettant l'obtention d'images volumiques d'un objet. Le système d'acquisition se compose d'un tube à rayons X et d'un détecteur plan numérique. La recherche développée dans ce manuscrit se déroule dans le contexte industriel. L'objet est placé sur une platine de rotation et une séquence d'images 2D est acquise. Un algorithme de reconstruction procure des données volumiques de l'atténuation de l'objet. Ces informations permettent de réaliser une étude métrologique et de valider ou non la conformité de la pièce imagée. La qualité de l'image 3D est dégradée par différents artefacts inhérents à la plateforme d'acquisition. L'objectif de cette thèse est de mettre au point une méthode de correction adaptée à une plateforme de micro-tomographie par rayons X d'objets manufacturés poly-matériaux. Le premier chapitre décrit les bases de la physique et de l'algorithmie propres à la technique d'imagerie CBCT par rayons X ainsi que les différents artefacts nuisant à la qualité de l'image finale. Le travail présenté ici se concentre sur deux types d'artefacts en particulier: les rayonnements secondaires issus de l'objet et du détecteur et le durcissement de faisceau. Le second chapitre présente un état de l'art des méthodes visant à corriger le rayonnement secondaire. Afin de quantifier le rayonnement secondaire, un outil de simulation basé sur des techniques de Monte Carlo hybride est développé. Il permet de caractériser le système d'acquisition installé au laboratoire de façon réaliste. Le troisième chapitre détaille la mise en place et la validation de cet outil. Les calculs Monte Carlo étant particulièrement prohibitifs en terme de temps de calcul, des techniques d'optimisation et d'accélération sont décrites. Le comportement du détecteur est étudié avec attention et il s'avère qu'une représentation 2D suffit pour modéliser le rayonnement secondaire. Le modèle de simulation permet une reproduction fidèle des projections acquises avec le système réel. Enfin, le dernier chapitre présente la méthodologie de correction que nous proposons. Une première reconstruction bruitée de l'objet imagé est segmentée afin d'obtenir un modèle voxélisé en densités et en matériaux. L'environnement de simulation fournit alors les projections associées à ce volume. Le volume est corrigé de façon itérative. Des résultats de correction d'images tomographiques expérimentales sont présentés dans le cas d'un objet mono-matériaux et d'un objet poly-matériaux. Notre routine de correction réduit les artefacts de cupping et améliore la description du volume reconstruit.

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal

L'imagerie cone-beam computed tomography (CBCT) est une méthodologie de contrôle non destructif permettant l'obtention d'images volumiques d'un objet. Le système d'acquisition se compose d'un tube à rayons X et d'un détecteur plan numérique. La recherche développée dans ce manuscrit se déroule dans le contexte industriel. L'objet est placé sur une platine de rotation et une séquence d'images 2D est acquise. Un algorithme de reconstruction procure des données volumiques de l'atténuation de l'objet. Ces informations permettent de réaliser une étude métrologique et de valider ou non la conformité de la pièce imagée. La qualité de l'image 3D est dégradée par différents artefacts inhérents à la plateforme d'acquisition. L'objectif de cette thèse est de mettre au point une méthode de correction adaptée à une plateforme de micro-tomographie par rayons X d'objets manufacturés poly-matériaux. Le premier chapitre décrit les bases de la physique et de l'algorithmie propres à la technique d'imagerie CBCT par rayons X ainsi que les différents artefacts nuisant à la qualité de l'image finale. Le travail présenté ici se concentre sur deux types d'artefacts en particulier: les rayonnements secondaires issus de l'objet et du détecteur et le durcissement de faisceau. Le second chapitre présente un état de l'art des méthodes visant à corriger le rayonnement secondaire. Afin de quantifier le rayonnement secondaire, un outil de simulation basé sur des techniques de Monte Carlo hybride est développé. Il permet de caractériser le système d'acquisition installé au laboratoire de façon réaliste. Le troisième chapitre détaille la mise en place et la validation de cet outil. Les calculs Monte Carlo étant particulièrement prohibitifs en terme de temps de calcul, des techniques d'optimisation et d'accélération sont décrites. Le comportement du détecteur est étudié avec attention et il s'avère qu'une représentation 2D suffit pour modéliser le rayonnement secondaire. Le modèle de simulation permet une reproduction fidèle des projections acquises avec le système réel. Enfin, le dernier chapitre présente la méthodologie de correction que nous proposons. Une première reconstruction bruitée de l'objet imagé est segmentée afin d'obtenir un modèle voxélisé en densités et en matériaux. L'environnement de simulation fournit alors les projections associées à ce volume. Le volume est corrigé de façon itérative. Des résultats de correction d'images tomographiques expérimentales sont présentés dans le cas d'un objet mono-matériaux et d'un objet poly-matériaux. Notre routine de correction réduit les artefacts de cupping et améliore la description du volume reconstruit. / Cone-beam computed tomography (CBCT) is a standard nondestructive imaging technique related to the acquisition of three-dimensional data. This methodology interests a wide range of applications. An industrial CBCT system comprises an X-ray source and a flat-panel detector. Radiographic images are acquired during a rotation of the object of interest. A reconstruction algorithm leads to a volumic representation of the object and a post-processing routine assesses its validity. Accurate quantitative reconstruction is needed to perform an efficient diagsnotic. However, it is challenged by the presence of different artefacts coming from the acquisition itself. This thesis aims at analyzing and correcting those artefacts in a context of industrial micro-tomography. After an introduction to the physical and algorithmic background of CBCT, the artefacts are presented. Our study adresses two major artefacts: beam hardening and scatter radiations coming from the object and the detector. The second chapter reports on the state of the art in secondary radiation correction. A simulation model of the CBCT imaging chain is developed in a Monte Carlo environment. This model is designed to be realistic in order to get an accurate insight on the processes contributing to the final image formation. The third chapter focuses on the built and validation of the simulation tool. Monte Carlo methods are exact but prohibitively slow. Consequently, acceleration and optimization techniques are used to speed-up the calculations without loss of accuracy. A layer model of the flat-panel detector gives some insight on its secondary radiation behavior. More specifically, we demonstrate that a 2D description of the detector would be sufficient to compute its contribution. Our projection tool fits well with the real system. Finally, the last chapter describes our iterative correction method. The noisy initial reconstruction is segmented into different materials and densities and fed to the simulation framework. Beam hardening and secondary radiations are corrected via the volume reconstructed from the difference between acquired and simulated projections. This correction method is shown to be effective on both mono-material and poly-material objects.