Filtering Deterministic Layer Effects in Imaging

Sensor array imaging arises in applications such as nondestructive evaluation of materials with ultrasonic waves, seismic exploration, and radar. The sensors probe a medium with signals and record the resulting echoes, which are then processed to determine the location and reflectivity of remote reflectors. These could be defects in materials such as voids, fault lines or salt bodies in the earth, and cars, buildings, or aircraft in radar applica- tions. Imaging is relatively well understood when the medium through which the signals propagate is smooth, and therefore nonscattering. But in many problems the medium is heterogeneous, with numerous small inhomogeneities that scatter the waves. We refer to the collection of inhomogeneities as clutter, which introduces an uncertainty in imaging because it is unknown and impossible to estimate in detail. We model the clutter as a random process. The array data is measured in one realization of the random medium, and the challenge is to mitigate cumulative clutter scattering so as to obtain robust images that are statistically stable with respect to different realizations of the inhomogeneities. Scatterers that are not buried too deep in clutter can be imaged reliably with the coher- ent interferometric (CINT) approach. But in heavy clutter the signal-to-noise ratio (SNR) is low and CINT alone does not work. The "signal," the echoes from the scatterers to be imaged, is overwhelmed by the "noise," the strong clutter reverberations. There are two existing approaches for imaging at low SNR: The first operates under the premise that data are incoherent so that only the intensity of the scattered field can be used. The unknown coherent scatterers that we want to image are modeled as changes in the coefficients of diffusion or radiative transport equations satisfied by the intensities, and the problem be- comes one of parameter estimation. Because the estimation is severely ill-posed, the results have poor resolution, unless very good prior information is available and large arrays are used. The second approach recognizes that if there is some residual coherence in the data, that is, some reliable phase information is available, it is worth trying to extract it and use it with well-posed coherent imaging methods to obtain images with better resolution. This paper takes the latter approach and presents a first attempt at enhancing the SNR of the array data by suppressing medium reverberations. It introduces filters, or annihila-

[1]  Guillaume Bal,et al.  Transport-Based Imaging in Random Media , 2008, SIAM J. Appl. Math..

[2]  P. Sheng,et al.  Introduction to Wave Scattering, Localization and Mesoscopic Phenomena. Second edition , 1995 .

[3]  Jean-Pierre Fouque,et al.  Spreading of a Pulse Travelling in Random Media , 1994 .

[4]  D. Malacara-Hernández,et al.  PRINCIPLES OF OPTICS , 2011 .

[5]  Miguel Moscoso,et al.  Robust imaging of localized scatterers using the singular value decomposition and ℓ1 minimization , 2013 .

[6]  J. Berryman Stable iterative reconstruction algorithm for nonlinear traveltime tomography , 1990 .

[7]  Andrey Bakulin,et al.  The virtual source method: Theory and case study , 2006 .

[8]  George Papanicolaou,et al.  Frequency Content of Randomly Scattered Signals , 1991, SIAM Rev..

[9]  Benjamin S. White,et al.  Localization and backscattering spectrum of seismic waves in stratified lithology , 1990 .

[10]  William W. Symes,et al.  All stationary points of differential semblance are asymptotic global minimizers: Layered acoustics 1 , 1999 .

[11]  S. Péché,et al.  Phase transition of the largest eigenvalue for nonnull complex sample covariance matrices , 2004, math/0403022.

[12]  J. Claerbout Earth Soundings Analysis: Processing Versus Inversion , 1992 .

[13]  E. Slob,et al.  Tutorial on seismic interferometry: Part 2 — Underlying theory and new advances , 2010 .

[14]  Arnaud Derode,et al.  Random matrix theory applied to acoustic backscattering and imaging in complex media. , 2009, Physical review letters.

[15]  Patrick Joly,et al.  An Analysis of New Mixed Finite Elements for the Approximation of Wave Propagation Problems , 2000, SIAM J. Numer. Anal..

[16]  Christiaan C. Stolk,et al.  Microlocal analysis of a seismic linearized inverse problem , 2000 .

[17]  Hongkai Zhao,et al.  Super-resolution in time-reversal acoustics. , 2002, The Journal of the Acoustical Society of America.

[18]  R. Adler,et al.  Random Fields and Geometry , 2007 .

[19]  Rami Nammour,et al.  Approximate Multi-Parameter Inverse Scattering Using Pseudodifferential Scaling , 2011 .

[20]  Liliana Borcea,et al.  Source Localization in Random Acoustic Waveguides , 2010, Multiscale Model. Simul..

[21]  Josselin Garnier,et al.  Use of random matrix theory for target detection, localization, and reconstruction , 2011 .

[22]  Jinho Baik,et al.  Asymptotics of Tracy-Widom Distributions and the Total Integral of a Painlevé II Function , 2007, 0704.3636.

[23]  Gregory Beylkin,et al.  Imaging of discontinuities in the inverse scattering problem by inversion of a causal generalized Radon transform , 1985 .

[24]  Liliana Borcea,et al.  Adaptive Time-Frequency Detection and Filtering for Imaging in Heavy Clutter , 2011, SIAM J. Imaging Sci..

[25]  George Papanicolaou,et al.  Transport equations for elastic and other waves in random media , 1996 .

[26]  F. Herrmann,et al.  Sparsity- and continuity-promoting seismic image recovery with curvelet frames , 2008 .

[27]  Liliana Borcea,et al.  Optimal illumination and wave form design for imaging in random media. , 2007, The Journal of the Acoustical Society of America.

[28]  Liliana Borcea,et al.  Filtering Random Layering Effects in Imaging , 2010, Multiscale Model. Simul..

[29]  Arnaud Derode,et al.  Detection and imaging in a random medium: A matrix method to overcome multiple scattering and aberration , 2009, 0906.0532.

[30]  Liliana Borcea,et al.  Adaptive interferometric imaging in clutter and optimal illumination , 2006 .

[31]  Josselin Garnier,et al.  A statistical approach to target detection and localization in the presence of noise , 2012 .

[32]  Josselin Garnier,et al.  Wave Propagation and Time Reversal in Randomly Layered Media , 2007 .

[33]  Liliana Borcea,et al.  Edge Illumination and Imaging of Extended Reflectors , 2008, SIAM J. Imaging Sci..

[34]  Sergey Fomel,et al.  Applications of plane-wave destruction filters , 2002 .

[35]  George Papanicolaou,et al.  Ray theory for a locally layered random medium , 2000 .

[36]  A. Böttcher,et al.  Introduction to Large Truncated Toeplitz Matrices , 1998 .

[37]  Laurent Demanet,et al.  Matrix probing: A randomized preconditioner for the wave-equation Hessian , 2011, 1101.3615.

[38]  Liliana Borcea,et al.  Detection and imaging in strongly backscattering randomly layered media , 2011 .

[39]  Sergey Fomel,et al.  Poststack velocity analysis by separation and imaging of seismic diffractions , 2007 .

[40]  Jack K. Cohen,et al.  Mathematics of Multidimensional Seismic Imaging, Migration, and Inversion , 2001 .

[41]  Liliana Borcea,et al.  Coherent interferometric imaging in clutter , 2006 .

[42]  Albert Fannjiang,et al.  Exact localization and superresolution with noisy data and random illumination , 2010, ArXiv.

[43]  B. Borden Mathematical problems in radar inverse scattering , 2002 .

[44]  Eliane Bécache,et al.  Étude d'un nouvel élément fini mixte permettant la condensation de masse , 1997 .

[45]  Liliana Borcea,et al.  Coherent Interferometry in Finely Layered Random Media , 2006, Multiscale Model. Simul..

[46]  A. Derode,et al.  Singular value distribution of the propagation matrix in random scattering media , 2009, 0904.0161.

[47]  R. Baraniuk,et al.  Compressive Radar Imaging , 2007, 2007 IEEE Radar Conference.

[48]  P. Schultz,et al.  Fundamentals of geophysical data processing , 1979 .

[49]  N. A. Anstey,et al.  Reflections on AMPLITUDES , 1971 .

[50]  V. Marčenko,et al.  DISTRIBUTION OF EIGENVALUES FOR SOME SETS OF RANDOM MATRICES , 1967 .

[51]  George Papanicolaou,et al.  Statistical Stability in Time Reversal , 2004, SIAM J. Appl. Math..

[52]  Josselin Garnier,et al.  Correlation-based virtual source imaging in strongly scattering random media , 2012 .

[53]  M. V. Rossum,et al.  Multiple scattering of classical waves: microscopy, mesoscopy, and diffusion , 1998, cond-mat/9804141.

[54]  Jon F. Claerbout,et al.  Fundamentals of Geophysical Data Processing: With Applications to Petroleum Prospecting , 1985 .

[55]  Liliana Borcea,et al.  Asymptotics for the Space-Time Wigner Transform with Applications to Imaging , 2007 .

[56]  W. Brown Synthetic Aperture Radar , 1967, IEEE Transactions on Aerospace and Electronic Systems.

[57]  Josselin Garnier,et al.  A universal filter for enhanced imaging with small arrays , 2009 .

[58]  William W. Symes,et al.  Velocity inversion by differential semblance optimization , 1991 .

[59]  G. Papanicolaou,et al.  Enhanced statistical stability in coherent interferometric imaging , 2011 .

[60]  Mathias Fink,et al.  Time-Reversed Acoustics , 1999 .

[61]  Liliana Borcea,et al.  Synthetic aperture radar imaging with motion estimation and autofocus , 2012 .

[62]  Plamen Stefanov,et al.  Recent progress on the boundary rigidity problem , 2005 .