A New Value Picking Regularization Strategy—Application to the 3-D Electromagnetic Inverse Scattering Problem

The nonlinear electromagnetic inverse scattering problem of reconstructing a possibly quasi-piecewise constant inhomogeneous complex permittivity profile is solved by iterative minimization of a pixel-based data fit cost function. Because of the ill-posedness it is necessary to introduce some form of regularization. Many authors apply a smoothing constraint on the reconstructed permittivity profile, but such regularization smooths away sharp edges. In this paper, a simple yet effective regularization strategy, the value picking (VP) regularization, is proposed. This new technique is capable of reconstructing piecewise constant permittivity profiles without degrading the edges. It is based on the knowledge that only a few different permittivity values occur in such profiles, the values of which need not be known in advance. VP regularization does not impose this a priori information in a strict sense, such that it can be applied also to profiles that are only approximately piecewise constant. The VP regularization is introduced in the solution of the inverse problem by adding a choice function to the data fit cost function for every permittivity unknown in the discretized problem. When minimized, the choice function forces the corresponding permittivity unknown to be close to one member of a set of auxiliary variables, the VP values, which are continuously updated throughout the iterations. To minimize the regularized cost function, a half quadratic Gauss-Newton optimization technique is presented. Finally, a stepwise relaxed VP regularization scheme is proposed, in which the number of VP values is gradually increased. This scheme is tested with synthetic and measured scattering data, obtained from inhomogeneous 3D targets, and is shown to achieve high reconstruction quality.

[1]  T. Isernia,et al.  Inverse scattering with real data: detecting and imaging homogeneous dielectric objects , 2001 .

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

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

[4]  F. Santosa,et al.  Reconstruction of a two-dimensional binary obstacle by controlled evolution of a level-set , 1998 .

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

[6]  Aria Abubakar,et al.  The contrast source inversion method for location and shape reconstructions , 2002 .

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

[8]  Michel Barlaud,et al.  Deterministic edge-preserving regularization in computed imaging , 1997, IEEE Trans. Image Process..

[9]  Michel Barlaud,et al.  Microwave imaging: Reconstructions from experimental data using conjugate gradient and enhancement by edge‐preserving regularization , 1997 .

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

[11]  Curtis R. Vogel,et al.  Ieee Transactions on Image Processing Fast, Robust Total Variation{based Reconstruction of Noisy, Blurred Images , 2022 .

[12]  P. M. Berg,et al.  Extended contrast source inversion , 1999 .

[13]  Michel Barlaud,et al.  Conjugate-gradient algorithm with edge-preserving regularization for image reconstruction from Ipswi , 1997 .

[14]  P. M. Berg,et al.  A total variation enhanced modified gradient algorithm for profile reconstruction , 1995 .

[15]  A. Franchois,et al.  Full-Wave Three-Dimensional Microwave Imaging With a Regularized Gauss–Newton Method— Theory and Experiment , 2007, IEEE Transactions on Antennas and Propagation.

[16]  Peter Monk,et al.  A modified dual space method for solving the electromagnetic inverse scattering problem for an infinite cylinder , 1994 .

[17]  O. Bucci,et al.  Representation of electromagnetic fields over arbitrary surfaces by a finite and nonredundant number of samples , 1998 .

[18]  A. Litman,et al.  Reconstruction by level sets of n-ary scattering obstacles , 2005 .

[19]  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..

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

[21]  P. M. van den Berg,et al.  Total variation as a multiplicative constraint for solving inverse problems. , 2001, IEEE transactions on image processing : a publication of the IEEE Signal Processing Society.

[22]  A. Abubakar,et al.  A General Framework for Constraint Minimization for the Inversion of Electromagnetic Measurements , 2004 .

[23]  P. Chaumet,et al.  Validation of a 3D bistatic microwave scattering measurement setup , 2008 .

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

[25]  Pierre Sabouroux,et al.  Electromagnetic three-dimensional reconstruction of targets from free space experimental data , 2008 .

[26]  J. Hadamard Sur les problemes aux derive espartielles et leur signification physique , 1902 .

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

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

[29]  J. D. Zaeytijd On the 3D electromagnetic quantitative inverse scattering problem: algorithms and regularization , 2009 .