A Time-Domain Multigrid Solver With Higher-Order Born Approximation for Full-Wave Radar Tomography of a Complex-Shaped Target

This article introduces and evaluates numerically a multigrid solver for non-linear tomographic radar imaging. Our goal is to enable the fast and robust inversion of sparse time-domain data with a mathematical full-wave approach utilizing a higher-order Born approximation (BA). Full-wave inversion is computationally expensive, hence techniques to speed up the numerical procedures are needed. To model the wave propagation effectively, we use the finite element time-domain (FETD) method, which is equipped with a multigrid scheme to enable the rapid evaluation of the higher-order BA. As a potential application, we consider the tomography of small solar system bodies (SSSBs) and asteroid interiors particular, the latter of which can contain internal details observable by radar, e.g., layers, voids and cracks. In the numerical experiments, we investigated monostatic, bistatic and multistatic measurement configurations. The results obtained suggest that, with a relevant noise level, the tomographic reconstruction quality can be improved by applying the higher-order BA in comparison to the first-order case, which we interpret as a linearization of the inverse problem. Our open multigrid-FETD solver for Matlab (The Mathworks Inc.) is available online. It applies Matlabs features for graphics computing unit acceleration to enhance computational performance.

[1]  G. M.,et al.  Partial Differential Equations I , 2023, Applied Mathematical Sciences.

[2]  E. Ebbini,et al.  Adaptive regularization of a distorted born iterative algorithm for diffraction tomography , 1996, Proceedings of 3rd IEEE International Conference on Image Processing.

[3]  Benoit Carry,et al.  Density of asteroids , 2012, 1203.4336.

[4]  Manfred Wittig,et al.  DISCUS – The Deep Interior Scanning CubeSat mission to a rubble pile near-Earth asteroid , 2018, Advances in Space Research.

[5]  M. Oelze,et al.  A study on the reconstruction of moderate contrast targets using the distorted born iterative method , 2008, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[6]  Weizhang Huang,et al.  Conditioning of finite element equations with arbitrary anisotropic meshes , 2012, Math. Comput..

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

[8]  D. Plettemeier,et al.  The Comet Nucleus Sounding Experiment by Radiowave Transmission (CONSERT): A Short Description of the Instrument and of the Commissioning Stages , 2007 .

[9]  Alain Herique,et al.  A Direct Observation of the Asteroid's Structure from Deep Interior to Regolith: Two Radars on the AIM Mission , 2016 .

[10]  F. Harris On the use of windows for harmonic analysis with the discrete Fourier transform , 1978, Proceedings of the IEEE.

[11]  James Irving,et al.  Numerical modeling of ground-penetrating radar in 2-D using MATLAB , 2006, Comput. Geosci..

[12]  Leone,et al.  Second-order iterative approach to inverse scattering: numerical results , 2000, Journal of the Optical Society of America. A, Optics, image science, and vision.

[13]  W. Chew,et al.  Fast inverse scattering solutions using the distorted Born iterative method and the multilevel fast multipole algorithm. , 2010, The Journal of the Acoustical Society of America.

[14]  Andreas Fhager,et al.  Effects of noise on tomographic breast imaging , 2011, 2011 XXXth URSI General Assembly and Scientific Symposium.

[15]  Carretera de Valencia,et al.  The finite element method in electromagnetics , 2000 .

[16]  A. P. Annan Chapter 1 – Electromagnetic Principles of Ground Penetrating Radar , 2009 .

[17]  Mahta Moghaddam,et al.  Real-time Microwave Imaging of Differential Temperature for Thermal Therapy Monitoring , 2014, IEEE Transactions on Biomedical Engineering.

[18]  Domenico Ciuonzo,et al.  On Time-Reversal Imaging by Statistical Testing , 2017, IEEE Signal Processing Letters.

[19]  Wlodek Kofman,et al.  A two dimensional simulation of the CONSERT experiment (radio tomography of comet Wirtanen) , 1999 .

[20]  A.J. Devaney Time reversal imaging of obscured targets from multistatic data , 2005, IEEE Transactions on Antennas and Propagation.

[21]  Eero P. Simoncelli,et al.  Image quality assessment: from error visibility to structural similarity , 2004, IEEE Transactions on Image Processing.

[22]  Roberto J. Lavarello,et al.  Tomographic Reconstruction of Three-Dimensional Volumes Using the Distorted Born Iterative Method , 2009, IEEE Transactions on Medical Imaging.

[23]  H. Eom Green’s Functions: Applications , 2004 .

[24]  Mika Takala,et al.  Multigrid-Based Inversion for Volumetric Radar Imaging With Asteroid Interior Reconstruction as a Potential Application , 2017, IEEE Transactions on Computational Imaging.

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

[26]  L. Kettunen,et al.  Yee‐like schemes on a tetrahedral mesh, with diagonal lumping , 1999 .

[27]  D. Plettemeier,et al.  Properties of the 67P/Churyumov-Gerasimenko interior revealed by CONSERT radar , 2015, Science.

[28]  M. Born Quantenmechanik der Stoßvorgänge , 1926 .

[29]  Katja Bachmeier,et al.  Finite Elements Theory Fast Solvers And Applications In Solid Mechanics , 2017 .

[30]  Pierluigi Salvo Rossi,et al.  Noncolocated Time-Reversal MUSIC: High-SNR Distribution of Null Spectrum , 2017, IEEE Signal Processing Letters.

[31]  Qing Huo Liu,et al.  Microwave Imaging in Layered Media: 3-D Image Reconstruction From Experimental Data , 2010, IEEE Transactions on Antennas and Propagation.

[32]  Wlodek Kofman,et al.  Microwave imaging from experimental data within a Bayesian framework with realistic random noise , 2009 .

[33]  A. Tarantola,et al.  Deconvolution and inverse theory: Application to Geophysical Problems (Methods in Geochemistry and Geophysics, Vol. 29) by Vijay Dimiri, Elsevier, Amsterdam, 1992, xviii + 230 pp., hardback, US $131.50, ISBN 0-444-89493-4 , 1993 .

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

[35]  Tommaso Isernia,et al.  Electromagnetic inverse scattering: Retrievable information and measurement strategies , 1997 .

[36]  Giovanni Leone,et al.  Inverse scattering under the distorted Born approximation for cylindrical geometries , 1999 .

[37]  Mika Takala,et al.  Far-Field Inversion for the Deep Interior Scanning CubeSat , 2017, IEEE Transactions on Aerospace and Electronic Systems.

[38]  W. Kofman,et al.  Radar techniques to study subsurfaces and interiors of the solar system objects , 2007, 2012 19th International Conference on Microwaves, Radar & Wireless Communications.

[39]  K. Yee Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[40]  A. Litman,et al.  A large 3D target with small inner details: A difficult cocktail for imaging purposes without a priori knowledge on the scatterers geometry , 2012 .

[41]  E. Grün,et al.  Comet nucleus sounding experiment by radiowave transmission , 1998 .

[42]  Giovanni Leone,et al.  Information content of Born scattered fields: results in the circular cylindrical case , 1998 .

[43]  L. Sorsa,et al.  Bistatic Full-wave Radar Tomography Detects Deep Interior Voids, Cracks, and Boulders in a Rubble-pile Asteroid Model , 2019, The Astrophysical Journal.

[44]  C. Eyraud,et al.  Imaging the interior of a comet from bistatic microwave measurements: Case of a scale comet model , 2017, Advances in Space Research.

[45]  A. Nuttall Some windows with very good sidelobe behavior , 1981 .

[46]  Amelie Litman,et al.  Optimization of a bistatic microwave scattering measurement setup: From high to low scattering targets , 2009 .

[47]  P. Sava,et al.  3D radar wavefield tomography of comet interiors , 2018 .

[48]  C. Eyraud,et al.  Influence of the description of the scattering matrix on permittivity reconstruction with a quantitative imaging procedure: polarization effects. , 2019, Journal of the Optical Society of America. A, Optics, image science, and vision.

[49]  D. Plettemeier,et al.  Direct observations of asteroid interior and regolith structure: Science measurement requirements , 2017, Advances in Space Research.

[50]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .

[51]  A. Ghatak,et al.  Elementary Theory of Scattering , 2004 .

[52]  Francesco Soldovieri,et al.  Dielectric profiles reconstruction via the quadratic approach in 2-D geometry from multifrequency and multifrequency/multiview data , 2002, IEEE Trans. Geosci. Remote. Sens..