Full-Wave Three-Dimensional Microwave Imaging With a Regularized Gauss–Newton Method— Theory and Experiment

A reconstruction algorithm is detailed for three-dimensional full-vectorial microwave imaging based on Newton-type optimization. The goal is to reconstruct the three-dimensional complex permittivity of a scatterer in a homogeneous background from a number of time-harmonic scattered field measurements. The algorithm combines a modified Gauss-Newton optimization method with a computationally efficient forward solver, based on the fast Fourier transform method and the marching-on-in-source-position extrapolation procedure. A regularized cost function is proposed by applying a multiplicative-additive regularization to the least squares datafit. This approach mitigates the effect of measurement noise on the reconstruction and effectively deals with the non-linearity of the optimization problem. It is furthermore shown that the modified Gauss-Newton method converges much faster than the Broyden-Fletcher-Goldfarb-Shanno quasi-Newton method. Promising quantitative reconstructions from both simulated and experimental data are presented. The latter data are bi-static polarimetric free-space measurements provided by Institut Fresnel, Marseille, France.

[1]  N. Joachimowicz,et al.  SICS: a sensor interaction compensation scheme for microwave imaging , 2002 .

[2]  Jean-Paul Hugonin,et al.  Microwave imaging-complex permittivity reconstruction-by simulated annealing , 1991 .

[3]  Weng Cho Chew,et al.  A discrete BCG-FFT algorithm for solving 3D inhomogeneous scatterer problems , 1995 .

[4]  van Mc Martijn Beurden,et al.  Iterative Solution of Field Problems with a Varying Physical Parameter , 2002 .

[5]  Ann Franchois,et al.  Quantitative microwave imaging with a 2.45-GHz planar microwave camera , 1998, IEEE Transactions on Medical Imaging.

[6]  Tim Hopkins,et al.  The Parallel Iterative Methods (PIM) package for the solution of systems of linear equations on para , 1995 .

[7]  P. M. Berg,et al.  The three dimensional weak form of the conjugate gradient FFT method for solving scattering problems , 1992 .

[8]  W. Chew,et al.  Reconstruction of two-dimensional permittivity distribution using the distorted Born iterative method. , 1990, IEEE transactions on medical imaging.

[9]  Qing Huo Liu,et al.  Three-dimensional nonlinear image reconstruction for microwave biomedical imaging , 2004, IEEE Transactions on Biomedical Engineering.

[10]  Marc Saillard,et al.  Retrieval of inhomogeneous targets from experimental frequency diversity data , 2005 .

[11]  Paul M. Meaney,et al.  Nonactive antenna compensation for fixed-array microwave imaging. II. Imaging results , 1999, IEEE Transactions on Medical Imaging.

[12]  Anton G. Tijhuis,et al.  Transient scattering by a lossy dielectric cylinder: marching-on-in-frequency approach , 1993 .

[13]  C. Pichot,et al.  Inverse scattering: an iterative numerical method for electromagnetic imaging , 1991 .

[14]  R. Fletcher.  Practical Methods of Optimization , 1988 .

[15]  Robert H. Svenson,et al.  Three-dimensional microwave tomography. Theory and computer experiments in scalar approximation , 2000 .

[16]  A. Kirsch.  An Introduction to the Mathematical Theory of Inverse Problems , 1996, Applied Mathematical Sciences.

[17]  Ignace Bogaert,et al.  An efficient hybrid MLFMA-FFT solver for the volume integral equation in case of sparse 3D inhomogeneous dielectric scatterers , 2008, J. Comput. Phys..

[18]  T. Sarkar,et al.  Comments on "Application of FFT and the conjugate gradient method for the solution of electromagnetic radiation from electrically large and small conducting bodies" , 1986 .

[19]  Three dimensional complex permittivity reconstruction by means of Newton-type microwave imaging , 2006, 2006 First European Conference on Antennas and Propagation.

[20]  Takashi Takenaka,et al.  Inverse scattering for a three-dimensional object in the time domain. , 2003, Journal of the Optical Society of America. A, Optics, image science, and vision.

[21]  Matteo Pastorino,et al.  A computational technique based on a real-coded genetic algorithm for microwave imaging purposes , 2000, IEEE Trans. Geosci. Remote. Sens..

[22]  Ann Franchois,et al.  Testing a 3D BCGS-FFT solver against experimental data , 2005 .

[23]  Anne Sentenac,et al.  Influence of multiple scattering on three-dimensional imaging with optical diffraction tomography. , 2006, Journal of the Optical Society of America. A, Optics, image science, and vision.

[24]  P. M. Berg,et al.  The diagonalized contrast source approach: an inversion method beyond the Born approximation , 2005 .

[25]  Amélie Litman,et al.  Theoretical and computational aspects of 2-D inverse profiling , 2001, IEEE Trans. Geosci. Remote. Sens..

[26]  P. M. Berg,et al.  Imaging of biomedical data using a multiplicative regularized contrast source inversion method , 2002 .

[27]  A. Bos.  Complex gradient and Hessian , 1994 .

[28]  Roger Fletcher,et al.  Practical methods of optimization; (2nd ed.) , 1987 .

[29]  P. M. Berg,et al.  A modified gradient method for two-dimensional problems in tomography , 1992 .

[30]  Andrea Massa,et al.  Full-vectorial three-dimensional microwave imaging through the iterative multiscaling strategy-a preliminary assessment , 2005, IEEE Geoscience and Remote Sensing Letters.

[31]  C. Eyraud,et al.  Free space experimental scattering database continuation: experimental set-up and measurement precision , 2005 .

[32]  V. Morozov.  On the solution of functional equations by the method of regularization , 1966 .

[33]  Siyuan Chen,et al.  Inverse scattering of two-dimensional dielectric objects buried in a lossy earth using the distorted Born iterative method , 2001, IEEE Trans. Geosci. Remote. Sens..

[34]  AG Anton Tijhuis,et al.  A quasi‐Newton reconstruction algorithm for a complex microwave imaging scanner environment , 2003 .

[35]  Aria Abubakar,et al.  Total variation as a multiplicative constraint for solving inverse problems , 2001, IEEE Trans. Image Process..

[36]  C. Pichot,et al.  Microwave imaging-complex permittivity reconstruction with a Levenberg-Marquardt method , 1997 .

[37]  P. M. Berg,et al.  A contrast source inversion method , 1997 .

[38]  A. N. Tikhonov,et al.  Solutions of ill-posed problems , 1977 .