TIME-DOMAIN INVERSE SCATTERING USING LAGRANGE MULTIPLIERS: AN ITERATIVE FDTD-BASED OPTIMIZATION TECHNIQUE

In this paper, we present an electromagnetic inverse scattering technique in the time domain, which is based on the minimization of an augmented cost functional. This functional is composed of an error term representing the discrepancy between the measured and the estimated values of the electromagnetic field, a regularization term related to the first-order Tikhonov's regularization scheme and, finally, an equality-constraints term associated with the fulfillment of the Maxwell's curl equations. These equality constraints are introduced by applying Lagrange multipliers. The minimization of the functional is carried out by applying the Polak-Ribière nonlinear conjugate-gradient algorithm. The required Fr´echet derivatives of the augmented cost functional with respect to the functions that describe the scatterer properties are derived by an analysis based on the calculus of variations. Actually, it is proven that the Lagrange multipliers are waves satisfying the Maxwell's curl equations. Consequently, during each iteration of the minimization procedure, we apply the finite-difference time-domain method and the perfectly-matched-layer absorber to calculate both the electromagnetic fields and the Lagrange multipliers. The presented technique is successfully applied to the reconstruction of multiple two-dimensional scatterers, which can be dielectric, lossy, and magnetic.

[1]  J. March Introduction to the Calculus of Variations , 1999 .

[2]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

[3]  Aria Abubakar,et al.  Iterative reconstructions of electrical conductivity from multiexperiment low-frequency electromagnetic data , 2000 .

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

[5]  Takashi Takenaka,et al.  FDTD approach to time‐domain inverse scattering problem for stratified lossy media , 1997 .

[6]  R. Kress,et al.  Inverse Acoustic and Electromagnetic Scattering Theory , 1992 .

[7]  Weng Cho Chew,et al.  Time-domain inverse scattering using the local shape function (LSF) method , 1993 .

[8]  Theodoros D. Tsiboukis,et al.  An iterative numerical method for inverse scattering problems , 1999 .

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

[10]  W. Chew,et al.  Study of some practical issues in inversion with the Born iterative method using time-domain data , 1993 .

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

[12]  A. Tijhuis Iterative determination of permittivity and conductivity profiles of a dielectric slab in the time domain , 1981 .

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

[14]  K. Paulsen,et al.  Two-dimensional hybrid element image reconstruction for TM illumination , 1995 .

[15]  R. J. Joseph,et al.  Advances in Computational Electrodynamics: The Finite - Di erence Time - Domain Method , 1998 .

[16]  Sailing He,et al.  An optimization approach to multi-dimensional time-domain electromagnetic inverse problems , 1998 .

[17]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[18]  M. Pastorino,et al.  Two-dimensional microwave imaging by a numerical inverse scattering solution , 1990 .

[19]  W. Chew,et al.  A frequency-hopping approach for microwave imaging of large inhomogeneous bodies , 1995, IEEE Antennas and Propagation Society International Symposium. 1995 Digest.

[20]  Theodoros D. Tsiboukis,et al.  A finite element-based technique for microwave imaging of two-dimensional objects , 2000, IEEE Trans. Instrum. Meas..