Méthodes de Reconstruction en Tomographie de diffraction 3-D

Ce travail porte sur la tomographie micro-ondes, technique d'imagerie permettant de reconstruire une image tridimensionnelle des propriétés diélectriques d'un volume inconnu à partir d'une succession de mesures du champ électromagnétique diffracté par ce volume. Cette modalité d'imagerie est prometteuse dans un grand nombre d'applications (imagerie médicale,géophysique, contrôle non destructif de matériaux, ...) mais souffre d'un inconvénient majeur, un coût de calcul très élevé des algorithmes de reconstruction, qui freine son développement industriel. Dans cette thèse, le problème d'imagerie micro-ondes est vu comme un problème inverse où l'on cherche à minimiser l'erreur entre les données et un modèle direct simulant la diffraction de l'onde en fonction des propriétés de l'objet inconnu. Le modèle physique est discrétisé par la méthode des moments. Ce problème est mal-posé, en particulier car le nombre d'inconnues est supérieur au nombre de mesures. La solution est donc définie par la minimisation d'un critère des moindres carrés pénalisé par une fonction de régularisation, dont l'optimisation requiert la mise en oeuvre d'algorithmes itératifs d'optimisation locale. Le modèle direct est en outre non-linéaire, rendant la reconstruction plus difficile. Le calcul numérique de ce modèle direct, nécessaire au calcul du critère et de son gradient à chaque itération, concentre la majorité du coût de calcul de l'inversion. Il repose sur la résolution d'un grand nombre de systèmes linéaires. Nous proposons dans cette thèse de réduire le coût de calcul de ces méthodes de reconstruction en se focalisant sur ces résolutions de systèmes. Pour cela, les travaux sont divisés en deux contributions. Tout d'abord, nous proposons une procédure permettant de réduire le nombre de systèmes à résoudre en s'appuyant sur la configuration du montage d'acquisition. La seconde contribution est d'accélérer les résolutions de systèmes à l'aide d'algorithmes par blocs. Ce type d'approche permet de résoudre conjointement plusieurs systèmes linéaires impliquant la même matrice opérateur. Ces méthodes sont validées sur des simulations de problèmes 3D réalistes, puis appliquées à la reconstruction d'objets réels sur des données expérimentales de champs diffractés.

Abstract

This thesis is focused on microwave tomography. This imaging technique consists in estimating a three-dimensional mapping of the dielectric properties of an unknown volume from measurements of the electromagnetic field from a known incident wave and scattered by this volume. This is a promising technique that is used in various applications (medical imaging, geophysics, non- destructive testing, ...) but suffers from high computational costs. This is a reason why microwave imaging is not widely used in industry. In this thesis, microwave imaging is considered as an inverse problem, where the error between the measurements and a forward model that describes the scattered field is minimized as a function of the properties of the volume. The physical model is discretized using the method of moments. This inverse problem is ill-posed because the number of unknowns is higher than the number of measurements. It is tackled through the minimization of a regularized leastsquares cost function, which is addresed by local iterative optimization algorithms. Moreover, the forward model is non-linear. Thus, reconstruction is a difficult and costful procedure. The computation of the objective function and of its gradient requires the resolution of a high number of linear systems, which are performed at each iteration of the optimization algorithm and represent most of the computational cost. In this thesis, we propose to reduce the computational costs of the reconstruction algorithms by focusing on the resolution of these linear systems. Two contributions are presented. The first one is a procedure in order to reduce the number of linear systems depending on the configuration of the measurement setup. The second contribution offers an efficient way to speed up the resolutions of the systems. We adapt block resolution algorithms, in order to jointly solve multiple linear systems involving a common operator matrix. These methods are validated on simulated, realistic, 3D problems, and applied to the reconstruction of real objects from experimental measurements of scattered fields. satisfactory results are obtained, where the computation time can be reduced by a factor of two, in particular for the most difficult reconstruction problems.