Green's function estimation in speckle using the decomposition of the time reversal operator: application to aberration correction in medical imaging

The FDORT method (French acronym for decomposition of the time reversal operator using focused beams) is a time reversal based method that can detect point scatterers in a heterogeneous medium and extract their Green's function. It is particularly useful when focusing in a heterogeneous medium. This paper generalizes the theory of the FDORT method to random media (speckle), and shows that it is possible to extract Green's functions from the speckle signal using this method. Therefore it is possible to achieve a good focusing even if no point scatterers are present. Moreover, a link is made between FDORT and the Van Cittert-Zernike theorem. It is deduced from this interpretation that the normalized first eigenvalue of the focused time reversal operator is a well-known focusing criterion. The concept of an equivalent virtual object is introduced that allows the random problem to be replaced by an equivalent deterministic problem and leads to an intuitive understanding of FDORT in speckle. Applications to aberration correction are presented. The reduction of the variance of the Green's function estimate is discussed. Finally, it is shown that the method works well in the presence of strong interfering scatterers.     

Time reversal and the inverse filter

To focus ultrasonic waves in an unknown inhomogeneous medium using a phased array, one has to calculate the optimal set of signals to be applied on the transducers of the array. In the case of time-reversal mirrors, one assumes that a source is available at the focus, providing the Green's function of this point. In this paper, the robustness of this time-reversal method is investigated when loss of information breaks the time-reversal invariance. It arises in dissipative media or when the field radiated by the source is not entirely measured by the limited aperture of a time-reversal mirror. However, in both cases, linearity and reciprocity relations ensure time reversal to achieve a spatiotemporal matched filtering. Nevertheless, though it provides robustness to this method, no constraints are imposed on the field out of the focus and sidelobes may appear. Another approach consists of measuring the Green's functions associated to the focus but also to neighboring points. Thus, the whole information characterizing the medium is known and the inverse source problem can be solved. A matrix formalism of the propagation operator is introduced to compare the time-reversal and inverse filter techniques. Moreover, experiments investigated in various media are presented to illustrate this comparison.     

Detection of cracks in a thin air-filled hollow cylinder by application of the DORT method to elastic components of the echo

The DORT method (Decomposition de l'Opérateur de Retournement Temporel in French) is a scattering analysis technique which uses arrays of transducers. This method is efficient for detection of selective focusing on point-like scatterers. It has been also applied to analyze the scattering by an air-filled cylindrical steel shell immersed in water. It was shown that the diagonalization of the time reversal operator allows us to separate the different elastic components of the scattered field. Here, we apply the method to detect flaws in hollow cylinders. In this case, the dominant components are the three circumferential waves (A0, A1 and S0 Lamb modes). Each Lamb mode corresponds to an invariant of the time reversal operator. The dispersion curves of these waves are calculated from these invariants. Resonance frequencies of the shell are deduced from the frequency dependence of the eigenvalues of the time reversal operator. It is shown that the presence of a crack (0.2 mm in depth) affects significantly the eigenvalue distribution of the time reversal operator. Thus, the DORT method offers a new means for detecting defects in a shell.     

Optimal focusing by spatio-temporal inverse filter. I. Basic principles

A focusing technique based on the inversion of the propagation operator relating an array of transducers to a set of control points inside a medium was proposed in previous work [Tanter et al., J. Acoust. Soc. Am. 108, 223-234 (2000)] and is extended here to the time domain. As the inversion of the propagation operator is achieved both in space and time, this technique allows calculation of the set of temporal signals to be emitted by each element of the array in order to optimally focus on a chosen control point. This broadband inversion process takes advantage of the singular-value decomposition of the propagation operator in the Fourier domain. The physical meaning of this decomposition is explained in a homogeneous medium. In particular, a definition of the number of degrees of freedom necessary to define the acoustic field generated by an array of limited aperture in a focal plane of limited extent is given. This number corresponds to the number of independent signals that can be created in the focal area both in space and time. In this paper, this broadband inverse-focusing technique is compared in homogeneous media with the classical focusing achieved by simple geometrical considerations but also with time-reversal focusing. It is shown that, even in a simple medium, slight differences appear between these three focusing strategies. In the companion paper [Aubry et al., J. Acoust. Soc. Am. 110, 48-58 (2001)] the three focusing techniques are compared in heterogeneous, absorbing, or complex media where classical focusing is strongly degraded. The strong improvement achieved by the spatio-temporal inverse-filter technique emphasizes the great potential of multiple-channel systems having the ability to apply completely different signal waveforms on each transducer of the array. The application of this focusing technique could be of great interest in various ultrasonic fields such as medical imaging, nondestructive testing, and underwater acoustics.     

Time reversed reverberation focusing in a waveguide

Time reversal mirrors have been applied to focus energy at probe source locations and point scatterers in inhomogeneous media. In this paper, we investigate the application of a time reversal mirror to rough interface reverberation processing in a waveguide. The method is based on the decomposition of the time reversal operator which is computed from the transfer matrix measured on a source-receiver array [Prada et al., J. Acoust. Soc. Am. 99, 2067-2076 (1996)]. In a similar manner, reverberation data collected on a source-receiver array can be filtered through an appropriate temporal window to form a time reversal operator. The most energetic eigenvector of the time reversal operator focuses along the interface at the range corresponding to the filter delay. It is also shown that improved signal-to-noise ratio measurement of the time reversal operator can be obtained by ensonifying the water column with a set of orthogonal array beams. Since these methods do not depend upon a priori environmental information, they are applicable to complex shallow water environments. Numerical simulations with a Pekeris waveguide demonstrate this method.     

Experimental detection and focusing in shallow water by decomposition of the time reversal operator

A rigid 24-element source-receiver array in the 10-15 kHz frequency band, connected to a programmable electronic system, was deployed in the Bay of Brest during spring 2005. In this 10- to 18-m-deep environment, backscattered data from submerged targets were recorded. Successful detection and focusing experiments in very shallow water using the decomposition of the time reversal operator (DORT method) are shown. The ability of the DORT method to separate the echo of a target from reverberation as well as the echo from two different targets at 250 m is shown. An example of active focusing within the waveguide using the first invariant of the time reversal operator is presented, showing the enhanced focusing capability. Furthermore, the localization of the scatterers in the water column is obtained using a range-dependent acoustic model.     

Revisiting iterative time reversal processing: application to detection of multiple targets

The iterative time reversal processing represents a high speed and easy way to self-focus on the strongest scatterer in a multitarget medium. However, finding weaker scatterers is a more difficult task that can be solved by computing the eigenvalue and eigenvector decomposition of the time reversal operator, the so-called DORT method. Nevertheless, as it requires the measurement of the complete interelements response matrix and time-consuming computation, the separation of multiple targets may not be achieved in real time. In this study, a new real time technique is proposed for multitarget selective focusing that does not require the experimental acquisition of the time reversal operator. This technique achieves the operator decomposition using a particular sequence of filtered waves propagation instead of computational power. Due to its simplicity of implementation, this iterative process can be achieved in real time. This high speed selective focusing is experimentally demonstrated by detecting targets through a heterogeneous medium and in a speckle environment. A theoretical analysis compares this technique to the DORT formalism.     

Detection and location of cracks using loss of reciprocity in ultrasonic waves propagation

The reciprocity theorem is a general statement valid for elastic media, and it has been applied to the solution of elastic wave equations, transducers calibration, time reversal acoustics, etc. However, localized nonlinear scatterers are expected to break reciprocity even though the effect is, in several cases, negligible. Here the dependence of the reciprocity break on the presence of a localized damage and the influence of its relative position has been experimentally investigated. It will be shown that the break of reciprocity, usually considered a disadvantage, can be exploited as an imaging tool for localized cracks detection.     

Numerical time reversal of waves

A numerical time reversal of waves is proposed instead of the conventional time reversal of wave procedure used in underwater acoustics. In the numerical method, the test sound source and the receiving arrays are used, as in the conventional method, but the transmission of the received signals after their time reversal into the same medium, as well as the measurement of the field obtained in this way at the point of the test source, is replaced by computations. To use the proposed technique for obtaining the same results as those provided by the conventional time reversal of waves, the teset source should be placed at different depths. A simplified numerical algorithm with the test source operating at a single depth is proposed and justified. This version of the time reversal of waves is successfully applied to the experiment in the Barents Sea. In contrast to the conventional method, the proposed technique allows one to study the stability of the sea medium with currents.     

Optimal adaptive focusing through heterogeneous media with the minimally invasive inverse filter

The inverse filter is a technique used to adaptively focus waves through heterogeneous media. It is based on the inversion of the Green's functions matrix between the M transducers of a focusing array and N control points in the focal area. The inverse filter minimizes the pressure deposited around the focal point. However it is highly invasive, requiring the presence of N transducers or hydrophones in the focal area at the control points' locations to measure the Green's functions. This paper presents a way of reaching the inverse filter's focusing quality with a minimally invasive setup: only one transducer (at the desired focal location) is needed. This minimally invasive inverse filter takes advantage of the fact all the information about the propagation medium can be retrieved from the signals backscattered by the medium towards the focusing array, if the propagation medium is lossless. A numerical simulation is performed to test this minimally invasive inverse filter through a scattering, lossless medium. The focusing quality equals the conventional, highly invasive inverse filter's. The average spatial and temporal contrast is increased by up to 10 dB compared to the time reversal focusing.     

Time reversal and its application to tomography with diffracting sources

An exact time-domain method is proposed to time reverse a transient scalar wave using only the field measured on an arbitrary closed surface enclosing the initial source. Under certain conditions, a time-reversed field can be approximated by retransmitting the measured signals in a reversed temporal order. Exact reconstruction for three-dimensional broadband diffraction tomography (a linearized inverse scattering problem) is proposed by time-reversing the measured field back to the time when each secondary source is excited. The algorithm is verified by a numerical simulation. Extension to the case using Green's function in a heterogeneous medium is discussed.     

Including dispersion and attenuation directly in the time domain for wave propagation in isotropic media

When sound propagates in a lossy fluid, causality dictates that in most cases the presence of attenuation is accompanied by dispersion. The ability to incorporate attenuation and its causal companion, dispersion, directly in the time domain has received little attention. Szabo [J. Acoust. Soc. Am. 96, 491-500 (1994)] showed that attenuation and dispersion in a linear medium can be accounted for in the linear wave equation by the inclusion of a causal convolutional propagation operator that includes both phenomena. Szabo's work was restricted to media with a power-law attenuation. Waters et al. [J. Acoust. Soc. Am. 108, 2114-2119 (2000)] showed that Szabo's approach could be used in a broader class of media. Direct application of Szabo's formalism is still lacking. To evaluate the concept of the causal convolutional propagation operator as introduced by Szabo, the operator is applied to pulse propagation in an isotropic lossy medium directly in the time domain. The generalized linear wave equation containing the operator is solved via a finite-difference-time-domain scheme. Two functional forms for the attenuation often encountered in acoustics are examined. It is shown that the presence of the operator correctly incorporates both, attenuation and dispersion.     

Some general problems of fast electrons transport

A generalization is given of the segments method in the form of a multistep method with generalized time for computing the transport of fast particles. The integral equation for a flow with generalized time in the phase space of variables is written under the assumption that the flow cuts the generalized time surface at right angles. The Green's function for the differential flow operator is the kernel of the integral equation. It is also shown that such an integral equation which can be obtained from a nonstationary kinetic equation provides a uniform consistent algorithm for solving either nonstationary or stationary problems. Examples of Green's functions are given for an operator of differential flow of fast electrons.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 8, pp. 110–114, August, 1974.The author would like to express his thanks to A. A. Vorob'ev and B. A. Kononov for their encouragement, to A. P. Yalovets for discussing the work with him, and to A. M. Kol'chuzhkin for going through the text.     

A time reversal damage imaging method for structure health monitoring using Lamb waves

This paper investigates the Lamb wave imaging method combining time reversal for health monitoring of a metal-lic plate structure.The temporal focusing effect of the time reversal Lamb waves is investigated theoretically.It demonstrates that the focusing effect is related to the frequency dependency of the time reversal operation.Numerical simulations are conducted to study the time reversal behaviour of Lamb wave modes under broadband and narrowband excitations.The results show that the reconstructed time reversed wave exhibits close similarity to the reversed nar-rowband tone burst signal validating the theoretical model.To enhance the similarity,the cycle number of the excited signal should be increased.Experiments combining finite element model are then conducted to study the imaging method in the presence of damage like hole in the plate structure.In this work,the time reversal technique is used for the recompression of Lamb wave signals.Damage imaging results with time reversal using broadband and narrowband excitations are compared to those without time reversal.It suggests that the narrowband excitation combined time reversal can locate and determine the size of structural damage more precisely,but the cycle number of the excited signal should be chosen reasonably.     

Transmission mode time-reversal super-resolution imaging

The theory of time-reversal super-resolution imaging of point targets embedded in a reciprocal background medium [A. J. Devaney, "Super-resolution imaging using time-reversal and MUSIC," J. Acoust. Soc. Am. (to be published)] is generalized to the case where the transmitter and receiver sensor arrays need not be coincident and for cases where the background medium can be nonreciprocal. The new theory developed herein is based on the singular value decomposition of the generalized multistatic data matrix of the sensor system rather than the standard eigenvector/eigenvalue decomposition of the time-reversal matrix as was employed in the above-mentioned work and other treatments of time-reversal imaging [Prada, Thomas, and Fink, "The iterative time reversal process: Analysis of the convergence," J. Acoust. Soc. Am. 97, 62 (1995); Prada et al., "Decomposition of the time reversal operator: Detection and selective focusing on two scatterers," J. Acoust. Soc. Am. 99, 2067 (1996)]. A generalized multiple signal classification (MUSIC) algorithm is derived that allows super-resolution imaging of both well-resolved and non-well-resolved point targets from arbitrary sensor array geometries. MUSIC exploits the orthogonal nature of the scatterer and noise subspaces defined by the singular vectors of the multistatic data matrix to form scatterer images. The time-reversal/MUSIC algorithm is tested and validated in two computer simulations of offset vertical seismic profiling where the sensor sources are aligned along the earth's surface and the receiver array is aligned along a subsurface borehole. All results demonstrate the high contrast, high resolution imaging capabilities of this new algorithm combination when compared with "classical" backpropagation or field focusing. Above and beyond the application of seismo-acoustic imaging, the time-reversal super-resolution theory has applications in ocean acoustics for target location, and ultrasonic nondestructive evaluation of parts.     

Time-dependent transport: Time domain recursively solving NEGF technique

We present a numerical method in this paper to calculate time-dependent transport properties in mesoscopic systems. This method recursively computes the non-equilibrium Green's functions (NEGF) in the time domain. In the simulations, we utilize a dynamically allocated data structure to compute current densities in response to input signals of any time duration and arbitrary shape. To demonstrate the performance, we apply the proposed method to solve NEGF for time-dependent transport in resonant-tunneling devices (RTD) with a single resonant level. We also investigate how the electrodes’ energy band structure affects the transport properties.     

Focusing of low-frequency sound fields in shallow water

A numerical experiment is carried out to study the focusing of a low-frequency (100–300 Hz) sound field in a shallow-water acoustic waveguide typical of an oceanic shelf. Focusing with the use of time reversal of broadband acoustic signals, which is called time reversal mirror (TRM) of waves, is considered along with focusing by phase conjugation (PC) of a monochromatic sound field. It is demonstrated that, in the case of focusing by the TRM method in the waveguide of interest, it is sufficient to have a single source-receiving element. The use of a vertical array improves the quality of focusing. The quality achieved in the latter case proves to be approximately the same as that achieved in the case of focusing by phase conjugation of a monochromatic field at a frequency identical to the carrier frequency of the broadband signals. It is also shown that, in a range-independent waveguide, intense surface waves considerably reduce the quality of focusing. This effect is most pronounced in the case of using phase conjugation.     

Discriminating resonant targets from clutter using Lanczos iterated single-channel time reversal

Power iterated single-channel time-reversal is extended to employ Lanczos iterations. The properties of these algorithms are studied in the presence of varying levels of noise and broadband clutter. It is shown the Lanczos iterated method possesses superior convergence properties in comparison to the standard power iterated technique. Results demonstrate that such algorithms provide an efficient means through which to isolate and extract the properties of resonant scatterers in the presence of noise and coherent interference.     

On optical properties of nonspherical inhomogeneous particles

To solve the problem of light scattering by multilayer scatterers of an arbitrary axisymmetric shape, a separation of variables method that involves special scalar potentials and their expansions in spherical functions is developed. The approach is shown to yield highly exact results even for particles that have 100 layers or more. A graphic library that illustrates the optical properties of layered and homogeneous (with an effective refractive index) spheroids, spheres, and Chebyshev particles of various shapes and sizes (about 650 figures) is created and is put on the Internet. It is noted that the linear polarization of radiation transmitted forward through a polydisperse medium containing partially oriented nonspherical porous particles strongly depends on the structure of scatterers. It is shown that the difference between the degrees of polarization of layered and corresponding homogeneous scatterers can exceed 200–300%.