A Preconditioned Inexact Newton Method for Nonlinear Sparse Electromagnetic Imaging

A nonlinear inversion scheme for the electromagnetic microwave imaging of domains with sparse content is proposed. Scattering equations are constructed using a contrast-source (CS) formulation. The proposed method uses an inexact Newton (IN) scheme to tackle the nonlinearity of these equations. At every IN iteration, a system of equations, which involves the Frechet derivative (FD) matrix of the CS operator, is solved for the IN step. A sparsity constraint is enforced on the solution via thresholded Landweber iterations, and the convergence is significantly increased using a preconditioner that levels the FD matrix's singular values associated with contrast and equivalent currents. To increase the accuracy, the weight of the regularization's penalty term is reduced during the IN iterations consistently with the scheme's quadratic convergence. At the end of each IN iteration, an additional thresholding, which removes small “ripples” that are produced by the IN step, is applied to maintain the solution's sparsity. Numerical results demonstrate the applicability of the proposed method in recovering sparse and discontinuous dielectric profiles with high contrast values.

[1]  Clifford H. Thurber,et al.  Parameter estimation and inverse problems , 2005 .

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

[3]  George V. Frisk,et al.  A comparison between the Born and Rytov approximations for the inverse backscattering problem , 1989 .

[4]  I. Daubechies,et al.  An iterative thresholding algorithm for linear inverse problems with a sparsity constraint , 2003, math/0307152.

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

[6]  José M. Bioucas-Dias,et al.  A New TwIST: Two-Step Iterative Shrinkage/Thresholding Algorithms for Image Restoration , 2007, IEEE Transactions on Image Processing.

[7]  R. Kress,et al.  Inverse Acoustic and Electromagnetic Scattering Theory , 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]  M. Pastorino,et al.  An Inexact Newton-Based Approach to Microwave Imaging Within the Contrast Source Formulation , 2009, IEEE Transactions on Antennas and Propagation.

[10]  Hakan Bagci,et al.  Shrinkage-Thresholding Enhanced Born Iterative Method for Solving 2D Inverse Electromagnetic Scattering Problem , 2014, IEEE Transactions on Antennas and Propagation.

[11]  Weng Cho Chew,et al.  An iterative solution of the two‐dimensional electromagnetic inverse scattering problem , 1989, Int. J. Imaging Syst. Technol..

[12]  M. Saillard,et al.  Special section: Testing inversion algorithms against experimental data , 2001 .