Open AccessReview

Ultrawideband Microwave Sensing and Imaging Using Time-Reversal Techniques: A Review

Mehmet Emre Yavuz

^1,* and

Fernando L. Teixeira

^2,*

Halliburton Energy Services, 3000 N Sam Houston Pkwy E, Houston, TX, 77032, USA

ElectroScience Laboratory, The Ohio State University, 1320 Kinnear Rd, Columbus, OH, 43212, USA

Authors to whom correspondence should be addressed.

Remote Sens. 2009, 1(3), 466-495; https://doi.org/10.3390/rs1030466

Submission received: 27 May 2009 / Revised: 13 July 2009 / Accepted: 10 August 2009 / Published: 24 August 2009

(This article belongs to the Special Issue Microwave Remote Sensing)

Download

Browse Figures

Versions Notes

Abstract

This paper provides an overview of some time-reversal (TR) techniques for remote sensing and imaging using ultrawideband (UWB) electromagnetic signals in the microwave and millimeter wave range. The TR techniques explore the TR invariance of the wave equation in lossless and stationary media. They provide superresolution and statistical stability, and are therefore quite useful for a number of remote sensing applications. We first discuss the TR concept through a prototypal TR experiment with a discrete scatterer embedded in continuous random media. We then discuss a series of TR-based imaging algorithms employing UWB signals: DORT, space-frequency (SF) imaging and TR-MUSIC. Finally, we consider a dispersion/loss compensation approach for TR applications in dispersive/lossy media, where TR invariance is broken.

Keywords:

time-reversal; superresolution; ultrawideband electromagnetic waves; time-domain DORT method; space-frequency (SF) imaging; TR-MUSIC method; short-time Fourier transform; dispersion compensation

Graphical Abstract

1. Introduction

Remote sensing systems in the microwave and millimeter wave (mm-wave) frequency ranges provide unique capabilities for detection and imaging of obscured targets. For example, microwaves can penetrate through foliage to detect obscured vehicles and personnel in forest environments. Mm-waves can penetrate through non-conductive walls and many packaging. They can also penetrate through dust, fog and smoke (as opposed to optical- and infrared-based systems). Even for scenarios where there is no distinctive advantage among various sensor modalities, microwave and mm-wave systems can play an important role because of their ease of integration into multimodal sensing apparatus [1].

Challenges for the full exploitation of microwave and mm-wave frequencies include vulnerability to certain atmospheric and meteorological phenomena and attenuation inside lossy materials. In addition, intervening media are often disordered with complex constitutive properties and/or many secondary scatterers about which precise information is not available. As a result, signals from the primary target(s) are often weak and/or masked by clutter and multipath (multiple scattering) contributions, which confounds detection, causes erratic tracking and makes it difficult to extract relevant information using conventional imaging and classification algorithms.

However, recent results have shown that multipath can actually be exploited to improve detection and imaging capabilities in sensing scenarios (or, somewhat equivalently, channel capacity in wireless communication scenarios) [2,3]. One such technique is time-reversal (TR) [4,5,6,7], which exploits multipath components in the intervening media to achieve superresolution, i.e., resolution that beats the classical diffraction limit. TR originated in acoustics and relies on the TR invariance of the wave equation in stationary and lossless media. TR involves physical or synthetic retransmission of signals acquired by a set of transceivers in a time-reversed fashion, i.e., last-in first-out. The retransmitted signals propagate “backwards”, naturally reversing the path that they underwent during forward propagation, which results in (automatic) energy focusing around initial source location(s). The “source” in this case can be either active (transmit mode) or passive (scatterer or “echo” mode).

The TR techniques that rely on ultrawideband (UWB) operation are further attractive because they can exploit advantages of simultaneous operation at low (e.g., more penetration into lossy materials) and high frequencies (e.g., better resolution), and because they enable imaging techniques in random media that depend only on the statistical properties (instead of a particular realization) of the random media, i.e., they are statistically stable [8,9,10].

The first successful TR experiments were carried by Fink and his collaborators using acoustic waves [11]. The underlying properties of TR operation were later investigated in more details both theoretically [10,12,13,14,15] and numerically [8,9,16,17,18,19].

Since then, TR techniques have experienced explosive growth in applications. Among the TR applications in medicine, one can cite lithotripsy [20,21], ultrasonic focusing through the skull [22], hyperthermia [23,24], remote inspection of internal masses [25] and health monitoring [26,27]. Applications of TR in geophysics and geoscience have included probing earthquake epicenters [28] and detecting objects buried underground [29]. Underwater applications have included sonar and acoustic communications [30,31], intruder detection, and echo-to-reverberation enhancement [32,33]. A commercial phased array sound system utilizing time-reversal concept is also developed in [34].

Nondestructive testing applications have utilized it for detecting defects in materials and structures [35,36,37,38,39].

In addition, imaging using synthetic TR techniques [40,41], TR interferometry [42,43,44] and imaging for sensor networks [45] have been addressed for both acoustic and electromagnetic waves. For most TR applications involving imaging, analysis of the so-called TR operator (TRO) is fundamental. In particular, the eigenspace analysis of the TRO provides important information about the scattering scenario under study. In the case of multiple targets for example, selective focusing using successive TR waves becomes possible. Perhaps the most basic approach that exploits the eigenvalue/vector structure of the TRO is the “DORT” (a French language acronym standing for “decomposition of the time-reversal operator") method [29,46,47,48,49,50,51]. In a similar vein, the “TR-MUSIC” method (derived from the MUltiple SIgnal Classification algorithm widely used for angle-of-arrival estimation in radar applications) has also been developed [40,52,53,54] and considered for a variety of applications [53,54,55,56,57].

In the frequency-domain, TR is equivalent to a phase-conjugation which has found applications in optics [58,59,60] as well as in electromagnetics where it is used to cancel distortions in the medium and provide beam steering action by using self-phasing or retrodirective array antennas [61,62,63,64,65]. In [66], it was experimentally shown that by using a single frequency phase-conjugation approach, energy focusing around two targets is also possible.

Among the most notable TR experiments using RF and microwaves, one can cite [67,68] where it was experimentally shown that a 1 MHz wideband pulse can be focused in a cavity environment.

Various TR-based target/change detection and interference cancelation algorithms have been discussed and investigated experimentally in [69,70,71,72,73].

Breast cancer detection using microwaves has been a topic of much interest [74]. TR is quite suited for this problem because of its natural ability to focus energy on malignant tissues (stronger scatterer). TR-based approaches for detection of breast cancer using microwaves are discussed in, e.g., [75,76,77,78,79,80,81].

More generally, TR-based electromagnetic imaging and target detection in cluttered environments with discrete scatterers have been experimentally investigated in [82,83,84].

The DORT method has been discussed in connection with remote sensing of buried objects in [29,49,85]. Additionally, the TR operator for electromagnetic waves in homogeneous medium was recently analyzed in [86,87].

TR-enabled focusing in a forest environment has been investigated numerically in [88,89]. Moreover, a TR synthetic aperture radar (TR-SAR) imaging method was developed and applied to an environment with strong multipath in [90,91]. Recently filed patent applications consider TR-based techniques for SAR images with improved target focusing and ghost image removal [92] as well as beamforming [93]. Additionally, TR technique has been applied for the through wall imaging in [94].

Finally, TR techniques have also found a number of applications in wireless communications such as the development of TR-based spatio-temporal matched filters to reduce channel dispersion and inter-symbol interference, thereby increasing the channel capacity [70,95,96,97,98,99,100,101,102]. Finally, a practical implementation of time-reversal of broadband microwave signals using optics has been demonstrated in [103] as well as in [104].

The objective of this review paper is to provide a summary of some TR-based techniques for detection and imaging applications using UWB microwaves. We first discuss TR operation in continuous random media which acts as a scenario where superresolution effects are manifested [105]. Then, we examine two examples from the “signal space” family of TR methods, viz. the time-domain (TD) DORT method [106], which allows selective focusing of multiple scatterers, and the space-frequency UWB TR imaging method [107], which utilizes spatial and frequency information of the received signals in tandem. Next, we discuss the “null space” TR-MUSIC method [108], which outperforms signal space methods under weak clutter conditions. Finally, we illustrate a dispersion/loss compensation approach aimed at extending TR operation to dispersive and lossy media, where TR invariance is broken [109].

2. Time-Reversal of Electromagnetic Waves and Superresolution

The TR concept exploits the invariance of the wave equation under TR, i.e., if

\bar{E} (\bar{r}, t)

is a solution of the wave equation in lossless media,

\nabla^{2} \bar{E} (\bar{r}, t) - μ (\bar{r}) ϵ (\bar{r}) \frac{\partial^{2}}{\partial t^{2}} \bar{E} (\bar{r}, t) = 0

(1)

then

\bar{E} (\bar{r}, - t)

is also a solution. Here,

\bar{E} (\bar{r}, t)

represents the electric field,

\bar{r} = \hat{x} x + \hat{y} y + \hat{z} z

is the position vector, and μ and ϵ are the permeability and permittivity of the medium, respectively. This property of the wave equation guarantees that for every wave propagating away from the source, there exists a reversed wave that would precisely retrace the path of the original wave back to the source. This occurs regardless of the presence of scattering objects and of inhomogeneities in μ and ϵ. The reversed waves would converge and focus coherently at the original source location as if time were flowing backwards,

t \to - t

This concept has been exploited by recording the wave using a set of transceivers and sending it back to its source as if time had been reversed, thereby yielding a variety of applications as mentioned in the Introduction. In practice, it is often not possible to completely surround the source so that TR is more generally performed with a set of receive/transmit antennas, the so-called TR array (TRA), that comprises only a limited angle. This causes some “information loss”; in general, the focal spot size gets larger as the array aperture size becomes smaller. Note that even though the term “array” is predominantly used in the TR literature, the antenna elements in the TRA ideally should work as independent transceivers, more akin to a “MIMO (multiple-input multiple-output) radar”. Recovering some of the “lost information” is possible through the exploitation of multipath components, as illustrated in Figure 1.

Figure 1. Time reversal experiment using a limited aspect array. (left) Forward Propagation of the short input pulse, (right) Backward Propagation of the time-reversed signals.

Here, a point source located at

{\bar{r}}_{0}

transmits a short UWB pulse

s (t)

. The transmitted signal propagates through the medium and is received by an antenna array. The signal received at the

i^{t h}

antenna is

f_{i} (t) = s (t) *_{t} h_{{\bar{r}}_{0} {\bar{r}}_{i}} (t)

(2)

where

*_{t}

denotes time convolution and

h_{{\bar{r}}_{0} {\bar{r}}_{i}} (t)

is the impulse response between the antennas at

{\bar{r}}_{0}

and

{\bar{r}}_{i}

. Reciprocity allows us to write

h_{{\bar{r}}_{0} {\bar{r}}_{i}} (t) = h_{{\bar{r}}_{i} {\bar{r}}_{0}} (t)

. The signals received at each array element are recorded, reversed in time, and transmitted back to the same medium (Figure 1). The time-reversed signal at the original source point

{\bar{r}}_{0}

due to

i^{t h}

antenna is then given by

p_{i} (t) = \underset{f_{i} (- t)}{\underset{︸}{s (- t) *_{t} h_{{\bar{r}}_{0} {\bar{r}}_{i}} (- t)}} *_{t} h_{{\bar{r}}_{i} {\bar{r}}_{0}} (t)

(3)

where the last two terms (

h_{{\bar{r}}_{0} {\bar{r}}_{i}} (- t) *_{t} h_{{\bar{r}}_{i} {\bar{r}}_{0}} (t)

) represent a correlation filter (time-correlator). This correlation function has a maximum at

t = 0

which corresponds to the energy of

h_{{\bar{r}}_{0} {\bar{r}}_{i}} (t)

, i.e.,

\int | h_{{\bar{r}}_{0} {\bar{r}}_{i}} {(t) |}^{2} d t

. With multiple antennas, TR system performance improves since each antenna will have a maximum at the original source location and they will constructively interfere to improve the TR peak signal. This coherent interference does not occur arbitrarily, but always at the original source location. For a TRA with N elements, the received signal at the original source location becomes:

p ({\bar{r}}_{0}, t) = \sum_{i = 1}^{N} s (- t) *_{t} h_{{\bar{r}}_{0} {\bar{r}}_{i}} (- t) *_{t} h_{{\bar{r}}_{i} {\bar{r}}_{0}} (t)

(4)

In addition to being a time-correlator, TR also acts as a space-correlator. In the above analysis, the TR waveform is exactly matched to the original source point

{\bar{r}}_{0}

. However, at any other point

\bar{r}

in the domain, this signal can be written as

p (\bar{r}, t) = \sum_{i = 1}^{N} s (- t) *_{t} h_{{\bar{r}}_{0} {\bar{r}}_{i}} (- t) *_{t} h_{{\bar{r}}_{i} \bar{r}} (t)

(5)

As the probe antenna location

\bar{r}

gets further away from the original source location

{\bar{r}}_{0}

, then, similar to time-correlation analysis, uncorrelated terms tend to cancel each other. For media with rich multipath components, correlation peak gets sharper and a better (sharper) focusing both in time and space can be achieved.

It is well known that the focusing spot size is dictated by the classical diffraction limit [110] which states that in an homogeneous media, cross-range (

d_{c}

) resolution is given by

d_{c} = λ \frac{L}{a}

(6)

where λ is the wavelength, a is the antenna aperture length and L is the distance between the array and the source (Figure 1). In scenarios with multipath, a should be replaced by the effective aperture length

a_{e}

as the inhomogeneities in the intervening medium affect the receiving pattern of the TRA. As long as some of the diverging wave components are redirected toward the TRA, the effective aperture length increases (

a_{e} > a

) resulting in a better focusing resolution than the homogeneous medium case (superresolution) [4]. This is as illustrated in Figure 2. We will next focus on superresolution effects created by the volumetric scattering in a continuous random media. The following subsection introduces the random medium model used in this work.

Figure 2. Effective aperture increase in media with multipaths. (a) Homogeneous medium with no multipaths (

a_{e} = a

), (b) multipaths created by the waveguide-like structure composed of two perfect electric conductor (PEC) walls (

a_{e} > a

), (c) multipaths created by discrete scatterers (

a_{e} > a

Figure 2. Effective aperture increase in media with multipaths. (a) Homogeneous medium with no multipaths (

a_{e} = a

), (b) multipaths created by the waveguide-like structure composed of two perfect electric conductor (PEC) walls (

a_{e} > a

), (c) multipaths created by discrete scatterers (

a_{e} > a

2.1. Physics-Based Ultrawideband Clutter Models

Complex natural media such as snow, vegetation, rocks, soils, and some biological tissues often cannot be described in a deterministic manner. Therefore, statistical models (random medium models) should be employed instead [111]. A random or disordered media can be classified either as (i) discrete random media, characterized as a discrete set of scatterers (e.g., trees, obstacles, buildings) at random locations, or as (ii) continuous random media, characterized by pointwise fluctuations on its properties (e.g., some biological media, soils, smoke) described in terms of a stochastic process with given first- and second-order statistics (spatial correlation functions).

With subsurface sensing scenarios in mind, we focus our attention on the characterization of continuous random media using constitutive parameters from particular soil models. Soil presents natural variability in density, composition and moisture that affect its permittivity and conductivity. Additionally, since in a host of applications the aim is to detect man-made buried objects, the soil has usually been excavated and therefore it is not expected to have a homogeneous or even layered distribution. In the absence of experimental data to support a specific choice of random medium model, continuous random medium models with Gaussian distributions are preferred for their generality and mathematical properties requiring few statistical parameters. Such random media are also characterized by a spatial correlation function. The relative permittivity of the medium is described as

ϵ (\bar{r}) = ϵ_{m} + ϵ_{f} (\bar{r})

where

ϵ_{m}

is the mean value of the relative permittivity and

ϵ_{f} (\bar{r})

is a function of position characterizing the random fluctuation on

ϵ (\bar{r})

with

〈ϵ_{f} (\bar{r})〉 = 0

. At every point in space, the fluctuation term is a Gaussian random variable with zero mean and probability density function (pdf) given by:

P_{ϵ_{f}} (ζ) = \frac{1}{\sqrt{2 π δ}} exp (\frac{- ζ^{2}}{2 δ})

(7)

The underlying assumption with this pdf is that the fluctuating dielectric permittivity is real. However, for more general cases where the dielectric permittivity is complex and random, a similar pdf can independently be used for both the real and imaginary parts of the dielectric permittivity. But, in this paper, we only consider lossless random medium and we refer the interested readers to another paper where we have considered both random and lossy (dispersive) media [106]. The random medium is characterized by transverse and vertical correlation lengths (

l_{s}

and

l_{z}

) and variance (δ), with the correlation function between the permittivities at two points also described by a Gaussian function

C ({\bar{r}}_{1} - {\bar{r}}_{2}) = 〈 ϵ_{f} ({\bar{r}}_{1}) ϵ_{f}^{*} ({\bar{r}}_{2}) 〉 = δ exp (- \frac{| {\bar{x}}_{1} - {\bar{x}}_{2} |^{2} + {| {\bar{y}}_{1} - {\bar{y}}_{2} |}^{2}}{l_{s}^{2}} - \frac{| z_{1} - z_{2} |^{2}}{l_{z}^{2}})

(8)

Previously, similar correlation functions have been used in [111,112,113]. Throughout this paper, we assume

l_{s} = l_{z}

so that we can isolate the effect of single correlation length on the TR focusing. Following these definitions, the procedure discussed in [114] can be used to generate the random medium realizations for numerical simulations. For modeling the scattering of embedded discrete target(s) in such media, we employ the finite-difference time-domain (FDTD) algorithm to solve Maxwell equations. The FDTD computational grid is truncated by perfectly matched layers (PML) [115] via stretched coordinates [116] to provide the necessary absorption of outgoing waves reaching the computational boundaries. A FDTD computational domain with

N_{x} \times N_{y} = 200 \times 240

grid points and uniform space discretization size

Δ_{x} = Δ_{y} = Δ_{s} = 0.0137

m is utilized in the examples considered below. A linear TRA of

N = 7

uniformly spaced dipole antennas are located just above a lossless random medium with spatially fluctuating permittivity. The TRA lies parallel to the x-direction and the dipoles are separated by

λ_{c} / 2

where

λ_{c}

corresponds to the wavelength at the central frequency for the mean ground permittivity. The location of the central antenna is assumed as origin, i.e.,

{\bar{R}}_{4} = (0, 0) Δ_{s}

where

{\bar{R}}_{i}

is the location of the

i^{t h}

TRA antenna. An electric dipole located at

{\bar{r}}_{s} = (x_{s}, y_{s}) = (0, 155) Δ_{s} = (0, 2.12) m

is initially fed by the current source

\bar{J} ({\bar{r}}_{s}, t) = \hat{e} s_{} (t) δ_{d} ({\bar{r}}_{s})

where

\hat{e}

(

\hat{x}

\hat{y}

\hat{z}

) is the unit vector representing the dipole polarization,

δ_{d} (\bar{r})

is the Dirac delta function and

s_{} (t)

is the UWB time-domain excitation taken as the first derivative of the Blackmann–Harris (BH) pulse [117] that vanishes after a time period of

T = 1.55 / f_{c}

, where

f_{c} = 400

MHz is the central frequency.

2.2. Numerical Results

In this section, we illustrate the superresolution effects versus random medium statistical parameters, viz. variance (δ) and correlation length (

l_{s}

We first investigate the effect of the variance of the random medium on the refocusing properties of the time reversed signals for the

T M_{z}

polarization case. The variance (δ) is varied from

0.025 ϵ_{m}

0.125 ϵ_{m}

, while the correlation length is kept fixed at

l_{s} = 8 Δ_{s}

. The snapshots of the z-component of the electric field (

E_{z}

) distribution at the time of refocusing are shown in Figure 3.

It is observed that as the variance increases, the amplitude of the focused field increases and the spot size of the focused field is reduced (somewhat counter-intuitively), characterizing superresolution. This can be explained from the fact that as the variance increases more multipath is produced. For the time reversed signals, the multipath tends to interfere constructively only at the focusing point where coherent perfect phase conjugation exists. This can also be understood as an increase on the effective aperture of the TRA due to multipathing [4,105]. Figure 5(a) shows the spatial distribution of the field components of Figure 3 at the time of refocusing at the source plane (

y = y_{s}

) with respect to the (transverse) x-coordinate (cross-range).

Figure 3. Spatial distribution of the time-reversed

E_{z}

field component of the electric field at the time of refocusing for increasing δ and fixed

l_{s}

l_{s} = 8 Δ_{s}

. Plots are given in linear scale.

Figure 3. Spatial distribution of the time-reversed

E_{z}

field component of the electric field at the time of refocusing for increasing δ and fixed

l_{s}

l_{s} = 8 Δ_{s}

. Plots are given in linear scale.

Figure 4. Cross-range resolution for varying (a) 1st and (b) 2nd order medium statistics.

Figure 5. (a) Comparison of incoherent and coherent time-domain signals. Note that although (b) the magnitudes of the frequency domain representation of both signals are same, (c) the phases are not. Incoherency comes from the oscillations in the phase arising from the SVD or EVD algorithm.

Similar simulations are performed to assess the effect of varying the correlation length on the resolution for the

T M_{z}

polarization. The correlation length is varied from

l_{s} = 8 Δ_{s}

l_{s} = 20 Δ_{s}

, for fixed variance of

δ = 0.04 ϵ_{m}

. Cross-range resolutions are shown in Figure 5(b), where it is seen that increasing the correlation length degrades the focusing properties. This can be explained from the fact that, for larger correlation lengths, the random medium starts to behave (locally) closer to an homogeneous medium and volumetric scattering effects are reduced. Note that this observation applies only for the range of correlation values and problem size considered here. If the correlation length becomes very small versus the wavelength, then a homogenization would be applicable. In general, the effect of changes on the correlation length on the focusing resolution is more pronounced when the correlation length is on the order of the wavelength. Another point to consider is the effect of polarization. As shown in [105], increased multiple scattering tends to increase depolarization, hence fully polarimetric data should be utilized for improved refocusing resolution. We also note that, although only the results of single realizations are shown here, another interesting phenomenon that occurs in connection with time-reversed back-propagation of wavefields in a random medium is the occurrence of the so-called time-domain statistical stability [8,9]. This means that the retrofocused field is self-averaging in the time-domain for a single realization and does not depend on the particular realization of the random medium. Self-averaging occurs because the time reversed field is equivalent to a phase conjugated field. After back-propagation, the conjugate phase exactly cancels out the random phase of the initial field and different frequency components add fully coherently only nearby the original source location, hence the focal spot is not affected by the particular realization.

3. Time-Reversal based Detection and Imaging Methods

In this section, we review some TR-based signal processing techniques to achieve UWB EM detection and imaging for multiple objects (targets) embedded in inhomogeneous media. This requires the analysis of the eigenspace of the TR operator (TRO) and utilization of its signal or null spaces. The TRO is obtained from the multistatic data matrix (MDM). The MDM element

{[K (ω)]}_{i j}

corresponds to the signal received at the ith antenna when the jth antenna is transmitting.

3.1. Signal Space Methods

3.1.1. Time-Domain DORT

The DORT method considers a

N \times N

MDM

K (ω)

(where ω is the frequency) obtained by probing the medium by a TRA of N transceivers. Since in this case each MDM element is associated with a different spatial location, this type of MDM is labeled as “space-space MDM”. In the frequency domain, TR of

K (ω)

is represented by its Hermitian conjugate

K^{†} (ω)

and the TRO is defined by

T (ω) = K^{†} (ω) K (ω)

. The SVD of the MDM is given by

K (ω) = U (ω) Λ (ω) V^{†} (ω)

(9)

where

U (ω)

and

V (ω)

are unitary matrices and

Λ (ω)

is real diagonal matrix of singular values. The EVD of the TRO can be written as

T (ω) = V (ω) S (ω) V^{†} (ω)

, where

S (ω) = Λ^{†} (ω) Λ (ω)

is the diagonal matrix of eigenvalues. The columns of the unitary matrix

V (ω)

are normalized eigenvectors of the TRO (

v_{i} (ω), i = 1, . ., N

). The DORT method utilizes the eigenvectors of the TRO signal space (SS) which is formed by the eigenvectors with non-zero (significant) eigenvalues, i.e.,

S S (ω) = {v_{1} (ω), . . ., v_{M_{s}} (ω)} with λ_{1}^{2} (ω) > . . > λ_{M_{s}}^{2} (ω) > 0

(10)

where

M_{s}

is the number of significant eigenvalues. For isotropic scattering from well-resolved point-like scatterers, each significant eigenvalue of the TRO is associated with a particular scatterer. Subsequent back-propagation of the corresponding eigenvector yields a wavefront focusing on that particular scatterer only [118]. Therefore, selective focusing on the mth scatterer is achieved by exciting the TRA with

N \times 1

column vector

r_{m} (ω)

generated from

λ_{m} (ω)

and eigenvector

v_{m} (ω)

via

r_{m} (ω) = K^{†} (ω) u_{m} (ω) = λ_{m} (ω) v_{m} (ω)

(11)

This represents the single-frequency DORT method. For UWB signals, EVD can be applied at all the available frequencies and a time-domain signal can be generated by

r_{p} (t) = F^{- 1} (K^{†} (ω) u_{p} (ω)) = F^{- 1} (λ_{p} (ω) v_{p} (ω))

(12)

where

F^{- 1}

denotes the inverse Fourier transformation. Backpropagation of these time-domain signals into the probed medium characterizes the TD-DORT method [106]. In cases where the background Green’s function is not known (which is the case for most subsurface sensing scenarios), an approximate Green’s function (

G_{0}

) can be used to obtain synthetic images of the probed medium by using the following TD-DORT imaging functional:

\begin{matrix} D_{p}^{Ω} ({\bar{X}}_{s}) & = & 〈g_{0} ({\bar{X}}_{s}, t), r_{p} (t)〉 ∣_{t = 0} \\ = & \sum_{i = 1}^{N} r_{p}^{i} (- t) *_{t} G_{0} ({\bar{X}}_{s}, {\bar{R}}_{i}, t) ∣_{t = 0} \\ = & \int_{Ω} λ_{p} (ω) v_{p}^{T} (ω) g_{0} ({\bar{X}}_{s}, ω) d ω \end{matrix}

(13)

where

g_{0} ({\bar{X}}_{s}, t) = {[G_{0} ({\bar{X}}_{s}, {\bar{R}}_{1}, t), \dots, G_{0} ({\bar{X}}_{s}, {\bar{R}}_{N}, t)]}^{T}

(14)

is the approximate time-domain steering vector connecting any search point

{\bar{X}}_{s}

in the probed domain to the antenna locations

{\bar{R}}_{i}

for

i = 1, . ., N

r_{p}^{i}

and

v_{p}^{i}

are the ith element of

r_{p}

and

v_{p}

, respectively and Ω is the bandwidth of operation. If one assumes, for example, that at least the average constitutive parameters of the random medium are known with a good degree of confidence, then

G_{0}

is chosen the Green’s function of an homogeneous medium with the mean constitutive parameters of the underlying random (inhomogeneous) medium that produces

K (ω)

. Note also that when the background medium is not known, energy detector which is a non-coherent imaging method can also be utilized [72].

For well-resolved scatterers, the cross-range resolution of the DORT method is related to the classical diffraction limit (Equation 6) which is directly proportional to the wavelength and propagation distance, and inversely proportional to the effective aperture length of the TR array. When the well-resolvedness criterion is broken, eigenvectors of the signal space of space-space MDMs might become linear combinations of the Green’s function vectors connecting the targets to the TR array, which results in a degradation (blurring) of the image. Another problem with the time-domain DORT stems from the eigenvalue decomposition step which creates eigenvectors with arbitrary frequency dependent phase

ϕ_{s v d} (ω)

. Direct combination of these eigenvectors creates incoherent time-domain signals which in turn severely affect the TD-DORT excitation signals. This is particularly important in inhomogeneous background media with strong multiple scattering where multipath components need to be coherently combined over the entire bandwidth. Therefore, a pre-processing step should be applied to the space-space eigenvectors to obtain coherent time domain signals by canceling the arbitrary phase term

ϕ_{s v d} (ω)

. Such pre-processing steps can be done by projecting the incoherent eigenvectors onto the columns of the MDM [8] or by a phase smoothing algorithm which tracks the phase difference between the adjacent frequency eigenvectors [106,119,120]. The time- and frequency-domain signals obtained for a random medium case with and without phase smoothing algorithm are shown in Figure 6(a). The effectiveness of the phase smoothing algorithm is evident from these plots.

Figure 6. The first eigenvalue distribution (most and only significant one in this case) of the space-space MDM with respect to the frequency (left) and singular values of the individual (center) and full (right) space-frequency MDMs.

3.1.2. Space-Frequency TR Imaging

An alternative approach to pre-processing is based on the application of the SVD to a different type of MDM which incorporates both sensor location and UWB frequency data simultaneously. Such MDM is denoted as “space-frequency (SF)” MDM. This approach was first introduced in connection with MUSIC-type algorithms [121]. Here, we illustrate them in the context of TR-based UWB imaging algorithms.

SF-MDMs can be obtained by transmitting an UWB pulse from the

n^{t h}

antenna and recording the received signals from all the TRA receiver antennas to yield an

N \times M_{f}

matrix given by

K_{S F}^{n} = (\begin{matrix} k_{1 n} (ω_{1}) & \dots & k_{1 n} (ω_{M_{f}}) \\ ⋮ & \dots & ⋮ \\ k_{N n} (ω_{1}) & \dots & k_{N n} (ω_{M_{f}}) \end{matrix})

(15)

where each row consists of the uniform samples of the Fourier transform of the time-domain signal corresponding to the respective MDM element, and

M_{f}

is the number of frequency samples. Once N of these individual SF-MDMs are obtained, SVD is applied to each of them to yield

K_{S F}^{n} = U_{S F}^{n} Λ_{S F}^{n} {(V_{S F}^{n})}^{†}

, where

U_{S F}^{n}

is the

N \times N

matrix of left singular vectors,

V_{S F}^{n}

is the

M_{f} \times M_{f}

matrix of right singular vectors, and

Λ_{S F}^{n}

is the

N \times M_{f}

matrix of singular values. Note that

K_{S F}^{n}

maps the frequency space onto the receiver space

K_{S F}^{n} v_{S F_{i}}^{n} = λ_{S F_{i}}^{n} u_{S F_{i}}^{n}

, where

λ_{S F_{i}}^{n}

is the

i^{t h}

singular value,

v_{S F_{i}}^{n}

is the

i^{t h}

M_{f} \times 1

right singular vector that represents the frequency content of the received signals, and

u_{S F_{i}}^{n}

is the

i^{t h}

N \times 1

left singular vector containing spatial (sensor location) information. The left singular vectors

u_{S F_{i}}^{n}

for

i = 1, \dots, N

form an orthonormal set spanning the sensor location space; similarly, right singular vectors

v_{S F_{i}}^{n}

for

i = 1, \dots, M_{f}

form an orthonormal set spanning the frequency space. Inverse Fourier transformation can be applied to the right singular vectors to obtain time-domain signals which are coherent. In other words, SVD applied to the SF-MDM does not create arbitrary phase dependent term as seen in the TD-DORT implementation. Therefore, the time-domain excitation signals to be backpropagated can be approximated as

s^{n} (t) = \sum_{k = 1}^{P} λ_{S F_{k}}^{n} v_{S F_{k}}^{n} (t)

(16)

where

v_{S F_{k}}^{n} (t) = F^{- 1} (v_{S F_{k}}^{n} (ω))

and P is the total number of time-domain signals being included in the approximation. This is determined by examining the singular values and the associated singular vectors. Time-domain signals corresponding to relatively small singular values and those with erratic behavior are not included. The resulting time-domain signals provide the UWB frequency data but they do not possess any sensor location data. Therefore, we use the left singular vectors to provide the necessary amplitude and phase shifts to be applied to each TR antenna during backpropagation. To this end, we define a vector functional

f (a, z (t)) = F^{- 1} {{[A_{0} e^{j ϕ_{0}} z (ω), \dots, A_{N} e^{j ϕ_{N}} z (ω)]}^{T}}

(17)

where

z (t)

is the time-domain signal to be used and

a = {[A_{0} e^{j ϕ_{0}}, \dots, A_{N} e^{j ϕ_{N}}]}^{T}

is the

N \times 1

vector to determine relative time-delays and amplitudes. This functional yields a time-domain vector

r_{S F_{i}}^{n} (t) = f (u_{S F_{i}}^{n}, s^{n} (t))

which is to be backpropagated from the receivers. If the background medium is not known precisely, an approximate Green’s functions can be utilized, which results in the following imaging functional:

\begin{matrix} I_{{S F}_{i}}^{n} (\bar{X_{s}}) & = & 〈 g_{0} ({\bar{X}}_{s}, t), r_{S F_{i}}^{n} (t) 〉 ∣_{t = 0} \\ = & \sum_{k = 1}^{N} r_{S F_{i k}}^{n} (- t) *_{t} G_{0} ({\bar{X}}_{s}, {\bar{R}}_{k}, t) ∣_{t = 0} \\ = & \int_{Ω} {(r_{S F_{i}}^{n} (ω))}^{†} g_{0} ({\bar{X}}_{s}, ω) d ω \\ = & \int_{Ω} {(s^{n} (ω))}^{*} {(u_{S F_{i}}^{n})}^{†} g_{0} ({\bar{X}}_{s}, ω) d ω \end{matrix}

(18)

where * represents complex conjugation. This procedure is repeated for all

n = 1, \dots, N

. The resulting image is averaged via

I_{S F_{i}} (\bar{X_{s}}) = \frac{1}{N} \sum_{n = 1}^{N} I_{{S F}_{i}}^{n} (\bar{X_{s}})

(19)

and denoted as “SF-image”. Note that instead of the left singular vectors, one can also utilize the eigenvectors obtained from the DORT method at the central frequency for the spatial information [107]. A second kind of SF-MDM can be obtained by combining all the

K_{S F}^{n}

for

n = 1, . ., N

into a single

N^{2} \times M_{f}

matrix given as

K_{S F}^{f u l l} = {[K_{S F}^{1}; K_{S F}^{2}; \dots; K_{S F}^{N}]}^{T}

(20)

This is denoted as the “full SF-MDM”. Similar to before, SVD applied to this matrix yields

N^{2} \times N^{2}

left and

M_{f} \times M_{f}

right singular vectors. The right singular vectors of the full SF-MDM are very similar to the individual SF-MDM and can be utilized to obtain time-domain excitation signals for backpropagation. The sub-vectors of the left singular vectors can provide the sensor location data for the SF imaging. For further details, readers are referred to [107]. Finally, we note that the phase and magnitude of the significant left singular vectors highly depend on the spatial distribution of the TR array antennas (sensor). As a result, inter-distances between sensors affect the phase and magnitude distributions of the left singular vectors. Similarly, for the DORT method, it will affect the phase and magnitudes of the eigenvectors. Depending on the scenario, this change can have positive or negative effects on the final imaging performance and has to be studied in more details which is reserved as a future study.

Next, we apply both TD-DORT and SF-imaging to a subsurface detection scenario to illustrate the similarities and differences between the two methods.

3.1.3. Results

To demonstrate the performance of DORT and SF-imaging methods, we employ the same simulation scenario as considered before except for the following differences: (a) we assume the central antenna located at the origin (i.e.

R_{4} = (0, 0) Δ_{s}

) and (b) a single scatterer at

(30, 80) Δ_{s}

is considered in both homogeneous and random media, the latter characterized by variance

δ = 0.03

and correlation length

l s = 6 Δ_{s}

. Figure 6 shows the corresponding eigenvalues for the TD-DORT method and singular values for the SF-imaging method.

It is observed that with increasing multiple scattering (variance), the magnitude of the eigenvalues and singular values increases. Additionally, the first dominant singular values obtained in homogeneous (HM) and random medium (RM) are close to each other and exhibit similar time-domain signatures (except for fluctuations due to clutter in the random medium case) as seen in Figure 7.

Figure 7. The first two significant time-domain right singular vectors obtained in homogeneous and random media using both the individual and full SF-MDM.

Among the remaining non-dominant singular values, those obtained in RM are larger than the HM ones. The corresponding time-domain right singular vectors also behave quite differently. For the RM case, the non-dominant right singular vectors are mainly due to clutter and not included in the approximation used for the excitation signal. Therefore, the first right singular vector can be used as the excitation signal with appropriate amplitudes and phase shifts dictated by the left singular vectors or TD-DORT eigenvectors at the desired frequency (Figure 8).

Figure 8. Phase distribution of the most significant TD-DORT eigenvector obtained at the central frequency (left) and those of the left singular vectors of the individual (middle) and full (right) SF-MDMs.

It is worth reminding that the excitation signals obtained via SF-MDM processing do not require phase smoothing, as opposed to the TD-DORT method. In certain cases where multiple scattering is very pronounced, phase smoothing might fail to yield coherent time-domain excitation signals, whereas SF-MDM processing will always yield coherent time domain signals. Figure 9 depicts the images obtained by TD-DORT and SF-imaging in homogeneous and random media. Note that these images correspond to the electric field distributions at the time of focusing and plotted in linear scale. Note that in terms of cross-range resolutions both methods provide similar performance. However, it was shown in [107] that in HM, TD-DORT slightly performs better. As alluded to above, in the RM case TD-DORT might suffer from phase-smoothing limitations under strong multiple scattering whereas SF-imaging is not affected by this limitation. We also note that it is not possible to claim that the SF-imaging provides improved resolution as the randomness (or multiple scattering) in the background increases. Especially, this is not possible as long as the approximate steering vector (Equation (14) is utilized. However, even for strong clutter cases where the TD-DORT might fail to generate coherent time-domain excitation signals, SF-imaging is able to generate coherent time-domain excitation signals which is an advantage of SF-imaging over the TD-DORT method. Finally, as long as the background medium is known and the exact background steering vectors are used for the imaging functionals, both methods might provide statistical stability thanks to the frequency decorrelation achieved by UWB operation.

It should be pointed out that some other signal space-based algorithms have also appeared in the literature, as in [122,123].

Figure 9. Images (in linear scale) obtained both in homogeneous (top row) and random (bottom row) media by TD-DORT (1st column), SF-imaging using individual (2nd column), and full SF-MDMs (3rd column).

3.2. Null Space Methods

Both DORT and SF-imaging methods rely on the fact that, for well-resolved point-like scatterers, information on scatterer strength(s) and location(s) are partially encoded by the eigenvalues and associated eigenvectors in the TRO signal space [47]. Backpropagation of these eigenvectors yields the target images. However, the performance of these algorithms degrades if the well-resolvedness criterion is not met. In this case, the SS eigenvectors become linear combinations of the Green’s function vectors connecting scatterers to the TRA [40]. Backpropagation of such SS eigenvectors creates overlapped image fields which hamper target imaging and localization. On the other hand, regardless of the well-resolvedness criteria, the TRO NS is always orthogonal to the TRO SS, i.e., projection of any vector formed by the linear combination of SS eigenvectors onto the NS is (ideally) zero. This property is the basis of TR-based MUSIC methods, which provides better detection and localization properties than signal space based methods even for poorly-resolved scatterers (assuming homogeneous media). A mathematical summary of the TR-MUSIC is presented next.

3.2.1. TR-MUSIC

Recalling Equation (9), the TRO NS is formed by the eigenvectors having near zero eigenvalues as shown below:

N S (ω) = {v_{M_{s} + 1} (ω), . . ., v_{N} (ω)} with λ_{M_{s} + 1}^{2} (ω) \approx . . \approx λ_{N}^{2} (ω) \approx 0

(21)

Note that, an NS (and thus the possibility to utilize TR-MUSIC) exists as long as the number of significant eigenvalues is less than the number of TRA antennas, i.e.,

M_{s} \leq N

. This may not hold for non-isotropic scattering, where more than one eigenvalue may be associated with a single scatterer [49,85,86,124]. Similarly, with increasing clutter and/or noise, additional TRO eigenvalues may become significant, which makes it increasingly difficult to distinguish the clutter/noise contribution from that of the discrete scatterers. In this case, a threshold criteria can be set to determine the SS and NS. This threshold may depend on the frequency and clutter/noise level.

In order to form the scatterer images, TR-MUSIC method requires the knowledge of the exact background Green’s function vector

g ({\bar{X}}_{s}, ω)

(steering vector) at each search point

{\bar{X}}_{s}

in the probed domain, which is defined by

g ({\bar{X}}_{s}, ω) = {[G ({\bar{X}}_{s}, {\bar{R}}_{1}, ω), . . ., G ({\bar{X}}_{s}, {\bar{R}}_{N}, ω)]}^{T}

(22)

where

G (\bar{r}, {\bar{r}}^{'}, ω) = | G (\bar{r}, {\bar{r}}^{'}, ω) | e^{j ϕ (\bar{r}, {\bar{r}}^{'}, ω)}

is the Green’s function of the problem [16,57]. The conjugate of the steering vector provides the necessary phase and amplitude distribution for the array excitation to focus on the desired point. Since

g ({\bar{X}}_{s}, ω)

is often not known in an deterministic fashion in inhomogeneous media, approximate steering vectors

g_{0} ({\bar{X}}_{s}, ω)

can be employed instead.

As mentioned above, as the well-resolvedness criterion for the scatterers is weakened, SS eigenvectors become linear combinations of the Green’s function vectors connecting the scatterers to the TRA. This is illustrated in Figure 10.

Figure 10. Illustration of well-resolved and non-well-resolved cases for a scenario having

M_{s} = 2

scatterers and probed with

N = 3

antennas. While the SS is formed by the plane (P) formed by the first two eigenvectors (

v_{1}

and

v_{2}

), NS is formed by the eigenvector (

v_{3}

) orthogonal to P.

Figure 10. Illustration of well-resolved and non-well-resolved cases for a scenario having

M_{s} = 2

scatterers and probed with

N = 3

antennas. While the SS is formed by the plane (P) formed by the first two eigenvectors (

v_{1}

and

v_{2}

), NS is formed by the eigenvector (

v_{3}

) orthogonal to P.

Imaging using these eigenvectors creates an interference pattern that degrades scatterer(s) location(s) estimates. However, even for closely spaced scatterers, NS is still orthogonal to the SS. Therefore, inner products of the steering vectors with the NS eigenvectors would vanish only at the scatterer(s) location(s). This provides the TR-MUSIC imaging functional (pseudospectrum) as given below:

M ({\bar{x}}_{s}, ω) = {(\sum_{i = M_{s} + 1}^{N} | 〈g_{0} ({\bar{x}}_{s}, ω), v_{i}^{*} (ω)〉 |)}^{- 1}

(23)

Since this scheme uses a single (central) frequency (CF)

ω_{c}

, the above is referred to as CF-TR-MUSIC, or simply CF-MUSIC in what follows.

In inhomogeneous media,

g_{0}

does not cancel out exactly the phase of the SS eigenvectors

v_{p} (= g^{*} / ∥ g ∥

) nearby the scatterer(s) location(s) as

g

would. Similarly, the inner product of

g_{0}

with the NS eigenvectors may not necessarily produce a minimum (ideally zero) at the original scatterer locations. In addition, since

g

has phase factors that depend on the (random medium) realization, using

g_{0}

does not produce statistical stability. In order to achieve statistical stability, one needs to explore frequency decorrelation. This can be achieved by UWB operation and combining images obtained at different frequencies via the following frequency integration [40]:

M_{Ω} ({\bar{x}}_{s}) = {(\int_{Ω} \sum_{i = M_{s} (ω) + 1}^{N} | 〈g_{0} ({\bar{x}}_{s}, ω), v_{i}^{*} (ω)〉 | d ω)}^{- 1}

(24)

where Ω is the bandwidth of operation. We call this strategy as UWB-TR-MUSIC, or simply UWB-MUSIC.

We apply both these TR-MUSIC algorithms to the same scenario of the previous subsection. Figure 11 shows the images obtained.

Figure 11. Images (in dB scale) obtained both in homogeneous (HM) and random media (RM) by CF-MUSIC and UWB-MUSIC.

One of the first observations is that the dynamic range of the MUSIC images are larger than the signal space methods (Figure 9). This is due to the fact that MUSIC pseudospectrum utilizes the null subspace and the inner products around the scatterer locations should ideally be zero resulting in a very strong peak at the scatterer locations (especially for the homogeneous media). A more detailed comparison of the dynamic ranges for varying clutter conditions can be found in [108]. In a homogeneous medium case, both TR-MUSIC methods provide good co-range and cross-range resolutions and they overperform those of the TD-DORT and SF-imaging. In the random medium case, the clutter affects the image quality. CF-MUSIC fails to image the scatterer location due to multiple spurious peaks, which can be misinterpreted as discrete scatterer locations. Moreover, the images vary with different random medium realizations because the imaging functionals utilize

g_{0}

and not

g

. Therefore, CF-MUSIC does not provide stable images under strong clutter. UWB-MUSIC, on the other hand, combines images obtained at different frequencies to yield a statistically stable image. This is produced at the expense of poorer co-range and cross-range resolutions (blurring) compared to the homogeneous case. Note that for both TR-MUSIC methods, increased clutter (i.e., decreased signal-to-clutter ratio) yields wider cross-range resolutions. Also, both TR-MUSIC methods fail to work in cases where NS is null, e.g., when all clutter eigenvalues are above the threshold. On the other hand, both DORT and SF-imaging methods provide similar co-range and cross-range information even for increasing clutter. As for processing times, it was shown in [108] that narrow-band methods (i.e., CF-MUSIC or CF-DORT) are faster than their UWB counterparts, and as long as the number SS eigenvectors is less than the number of NS eigenvectors, SS based methods are faster than the NS based methods. Additionally, TR-MUSIC is more stable against increasing dispersion and medium loss as compared to the SS based methods. Further details on the stability of the methods to perturbations can be found in [108,125]. There are also several other variant TR-MUSIC algorithms as considered, e.g., in [126,127]. Finally, in this paper, we have assumed that the number of distinct scatterers are less than the total number of antennas. For scenarios where the scatterers are more than the number of antennas, the methods explained in this paper may not be sufficient for imaging. As discussed in [72], a practical approach is to suppress the clutter first via an anti-focusing step and then attempt to image the scatterers of interest. This algorithm has shown to have good performance over conventional approaches when the number of scatterers is larger than the TRA sensors. So as a future study, we plan to incorporate this approach into the TR-imaging algorithms presented here. Especially, the effects of anti-focusing on the eigenspace and singular value/vectors structures will be focused.

4. Frequency Dispersion/Loss Compensation Techniques

The methods discussed so far have assumed that the intervening medium is lossless and stationary, hence the invariance of the wave equation under TR remains valid. However, dispersive and lossy media are often encountered in nature, such as in soils, rocks, ice, and most biological tissues. In such cases, the TR invariance is broken and conventional TR operation can not be directly applied. Several compensation methods (within certain limitations) have been discussed in the literature [22,75,128,129]. In this section, we review a more robust compensation method developed for UWB signals propagating in dispersive media [106,109].

A dispersive medium acts as a filter for UWB signals that propagate in it. For such a medium, the electric field at a single frequency is represented as

\bar{E} (\bar{r}, ω) = \bar{F} (\bar{r}) e x p (- j ω \sqrt{μ ϵ (ω)} \bar{r})

(25)

where

\bar{F} (\bar{r})

is a frequency-independent amplitude coefficient and for simplicity the permittivity function

ϵ (ω)

is here assumed to be invariant over space (

ϵ (\bar{r}, ω) = ϵ (ω)

) [130, pg. 18].

The non-zero imaginary part of

ϵ (ω)

yields additional attenuation on the received signals by the TRA relative to the non-dispersive case as shown in Figure 12. Moreover, the larger real part of

ϵ (ω)

causes an additional phase shift (delay). Since TR signals are phase conjugated coherently along the entire bandwidth, any additional phase shift induced from dispersion during the forward propagation is exactly compensated by the TR process. In a dispersive medium however, the signal is attenuated during both forward- and back-propagation. As a result, refocusing of the TR signals can be significantly degraded relative to the non-dispersive (lossless) case. To overcome this degradation, a compensation technique should be applied to act as an inverse filter with respect to the attenuation. Additionally, dispersive effects are also cumulative in time, i.e., the longer a signal travels in the medium, the more it is attenuated. In other words, a dispersive medium acts as a frequency-dependent, cumulative-in-time filter to the propagating signals. This requires the use of space and frequency dependent compensation filters. The implementation of the compensation method does not change the basic TR experiment, i.e., forward and backpropagation steps stay the same. But received and recorded signals after the forward propagation are properly modified (compensated) and then retransmitted. During the forward propagation of the transmitted short pulse in the dispersive medium, different multipaths go through different (dispersive) relative attenuations depending on the total distance (or time) traveled in the medium. In other words, received signals at later times (or equivalently, at further distances) need more strong “compensation” than those at earlier times. Ideally, different filters should be applied for each particular instant (sample), but this is impractical. Instead, time-windowing can be applied to received signals and different frequency-dependent filters used on each window. The combined time-windowing/filtering process can be best achieved after applying a short-time Fourier transform (STFT) [132] to the received signals at the TRA.

Figure 12. Magnitude of the exponential term

e x p (- j ω \sqrt{μ ϵ (ω)} \bar{r})

plotted with respect to frequency and spatial distance for complex

ϵ (ω)

corresponding to the Puerto Rico soil with moisture level of %2.5 [131]. Space- and frequency-dependent attenuation is observed.

Figure 12. Magnitude of the exponential term

e x p (- j ω \sqrt{μ ϵ (ω)} \bar{r})

plotted with respect to frequency and spatial distance for complex

ϵ (ω)

corresponding to the Puerto Rico soil with moisture level of %2.5 [131]. Space- and frequency-dependent attenuation is observed.

The whole compensation process is relatively simple and starts by taking the discrete STFT,

X [k, i]

, of the (sampled) received signal at each TRA antenna

x [n]

, which is given as

X [k, i] = \sum_{n} x [n] w_{i} [n] e^{- j (2 π / L) n k} for i = 1, . ., M

(26)

where

w_{i} [n]

is the

i^{t h}

window, L is the FFT length, and M is the total number of windows. Note that M depends on the window length and overlapping factor (between windows). The product

x [n] w_{i} [n] = s_{i} [n]

is referred to as the

i^{t h}

windowed signal. Each windowed signal has traveled a different distance amount to reach the receiver and each of its frequency components has gone through different attenuation. Therefore, each windowed signal has to be filtered with a different space and frequency dependent filter

H [k, i]

obtained for the

i^{t h}

window. In the Fourier domain, the filtering is carried out by

Y [k, i] = X [k, i] H [k, i] for i = 1, . ., M

(27)

The final, dispersion-compensated signal

x_{c} [n]

is then obtained by an inverse STFT:

x_{c} [n] = \sum_{i = 1}^{M} \frac{1}{w_{i} [n]} \sum_{k} Y [k, i] e^{j (2 π / L) n k}

(28)

A block diagram of the whole procedure is given in Figure 13.

Filters play an important role here and to obtain them, the dispersive characteristics of the media should be (ideally) known or (practically) estimated. One way to obtain these filters is to compare the solution of the wave equation in a homogeneous reference medium having the exact dispersive model of interest with that in a nondispersive reference medium, as carried out in [109].

Figure 13. The block diagram of the dispersion compensation method.

Now, let us illustrate the compensation procedure by considering the received signal at one of the TRA elements is considered in more detail. The received signal is windowed using five Hamming windows [132] of length

256 Δ t

each and an overlapping factor of

0.5

, as shown in Figure 14.

Figure 14. (Left:) One of the original received signal and the employed Hamming windows with overlapping factor of 0.5; (right:) corresponding windowed signals. Note the amplitude difference for each windowed signal.

For each windowed signal, a filter is designed as follows: First, for each central point of window

w_{i} [n]

, a corresponding effective distance is found for the specific dispersive characteristics. The ratio between the frequency domain representations of the signals received in the dispersive and non-dispersive reference media at this effective distance gives the frequency-dependent attenuation due to dispersion to be compensated by

H [k, i]

. Note that, to avoid noise contamination, the amplification factors can be smoothly set to unity for frequencies where the spectral density are below acceptable levels (high-frequency ends) depending on the particular application. Then, these filters are applied to corresponding windowed signals to yield the filtered windowed signals shown in Figure 15.

Figure 15. Time and Frequency domain representations for some of the windowed signals and their compensated counterparts after space and frequency dependent filtering.

These filtered windowed signals are then inverse Short-Time Fourier transformed to obtain the final compensated signal, which is shown in Figure 16 along with the original signal received in the dispersive medium and a reference signal that would have been received in the non-dispersive case.

Figure 16. Signal received by one of the TRA antennas in dispersive medium, corresponding compensated signal and the reference signal that would be received in non-dispersive medium in (left) time and (right) frequency domains.

As seen in these plots, apart from the phase shift (which is automatically compensated by the TR process), the amplitude of the compensated signal is much closer to the reference signal amplitude than that of the original signal. Any discrepancy between the compensated and the reference signal is due to the finite window length used. A better agreement can be obtained with a larger number of windows, however, trade-off is the increased computational cost. This procedure is applied to all the signals received by the TRA. Note that, since backpropagation occurs in the same dispersive medium, additional attenuation is produced on the signals. Therefore, for each windowed signal, two different sets of compensation filters should be applied for each windowed signal. The first one compensates for propagation losses from the source (TRA) to scatterers and from scatterers back to the TRA; the second one compensates for propagation losses from TRA to scatterers. The compensation filters used for backpropagation should be different from those for the initial propagation, since the effective propagation duration for the backpropagation is approximately half of the initial propagation. The second set of filters for each windowed signal can be obtained by using half of the effective distances used for the first set of compensation filters. Once these compensated and TR signals are backpropagated into the same medium, the refocused fields around the scatterers have amplitudes much closer to those obtained in the nondispersive case, as shown in [106,109]. It should also be noted that compensation filters also amplify system noise; hence, they are of limited availability in scenarios where the medium dispersion reduces the signal to noise ratio below acceptable levels.

5. Conclusions

In this paper, we have provided a summary of some time-reversal techniques for UWB microwave remote sensing. For concreteness, we have focused most of the discussion on imaging scenarios consisting of obscured discrete targets in continuous random media, although TR techniques are applicable to many other remote sensing scenarios, as surveyed in the Introduction.

First, we considered the so-called “signal space” methods, viz., TD-DORT and SF-imaging, which respectively utilize space-space and space-frequency multistatic data matrices obtained from the TR array signals. For the TD-DORT method, the use of an eigenvector decomposition of the TR operator along the available bandwidth followed by a phase smoothing algorithm (on the eigenvectors associated with the significant eigenvalues) provides the required time-domain excitation signals to be used on the TRA for selective focusing of multiple targets. While TD-DORT performs well in relatively low clutter environments, it fails in environments with stronger clutter because of the limitations in the phase smoothing algorithm. SF-imaging, on the other hand, utilizes less conventional MDMs to provide coherent time-domain excitation signals regardless of the background clutter. Although the TD-DORT cross-range performance is slightly better than that provided by SF-imaging, the latter is favored in highly scattering environments.

Next, we considered the TR-MUSIC method as the prototypal “null space” based TR algorithm. In this case, the TR operator eigenvectors associated with near zero eigenvalues of the TR operator (as opposed to the significant eigenvalues, as in the “signal space” case) are utilized along with medium steering vectors to obtain desired target images. Under weak clutter conditions, TR-MUSIC outperforms TD-DORT and SF-imaging methods in terms of both co-range and cross-range resolutions. However, TR-MUSIC was found to be less stable under increased clutter.

Finally, we briefly discussed the use of compensation (inverse) filters as a means to extend the basic TR techniques to lossy media, where the TR invariance is broken.

Acknowledgements

This work has been supported in part by the National Science Foundation (NSF) under Grant ECS-0347502 and the Ohio Supercomputing Center (OSC) under Grants PAS-0061 and PAS-0110.

References and Notes

Varshney, P.K. Distributed Detection and Data Fusion; Springer-Verlag: Berlin, Germany, 1997. [Google Scholar]
Paulraj, A.J.; Gore, D.A.; Nabar, R.U.; Bolcskei, H. An overview of MIMO communications–a key to gigabit wireless. Proc. IEEE 2004, 92, 198–218. [Google Scholar] [CrossRef]
Foschini, G.J.; Gans, M.J. On limits of wireless communication in a fading environment when using multiple antennas. Wirel. Personal Commun. 1998, 6, 311–335. [Google Scholar] [CrossRef]
Fink, M.; Cassereau, D.; Derode, A.; Prada, C.; Roux, P.; Tanter, M.; Thomas, J.; Wu, F. Time-reversed acoustics. Rep. Prog. Phys. 2000, 63, 1933–1995. [Google Scholar] [CrossRef]
Fink, M. Method and device for localization and focusing of acoustic waves in tissues. US Patent 5,092,336, 1992. [Google Scholar]
Fink, M. Method and apparatus for acoustic examination using time reversal. US Patent 5,428,999, 1995. [Google Scholar]
Fink, M.; Lewiner, J. Method and device for detecting and locating a reflecting sound source. US Patent 6,161,434, 2000. [Google Scholar]
Berryman, J.G.; Borcea, L.; Papanicolaou, G.; Tsogka, C. Statistically stable ultrasonic imaging in random media. J. Acoust. Soc. Am. 2002, 112, 1509–1522. [Google Scholar] [CrossRef] [PubMed]
Papanicolaou, G.; Ryzhik, L.; Solna, K. Statistical stability in time reversal. SIAM J. Appl. Math. 2004, 64, 1133–55. [Google Scholar] [CrossRef]
Fouque, J.-P.; Garnier, J.; Papanicolaou, G.; Solna, K. Wave Propagation and Time Reversal in Randomly Layered Media; Springer: Berlin, Germany, 2007. [Google Scholar]
Fink, M. Time reversal mirrors. J. Phys. D.: Appl. Phys. 1993, 26, 1333–1350. [Google Scholar] [CrossRef]
Clouet, J.F.; Fouque, J.P. A time reversal method for an acoustical pulse propagating in randomly layered media. Wave Motion 1997, 25, 361–368. [Google Scholar] [CrossRef]
Dowling, D.R.; Jackson, D.R. Narrowband performance of phase conjugate arrays in dynamic random media. J. Acoust. Soc. Am. 1992, 91, 3257–3277. [Google Scholar] [CrossRef]
Asch, M.; Kohler, W.; Papanicolaou, G.; Postel, M.; White, B. Frequency content of randomly scattered signals. SIAM Review 1991, 33, 519–62. [Google Scholar] [CrossRef]
Kohler, W.; Papanicolaou, G.; White, B. Localization and mode conversion for elastic waves in randomly layered media. Wave Motion 1996, 23, 1–22 and 181–201. [Google Scholar] [CrossRef]
Borcea, L.; Papanicolaou, G.; Tsogka, C.; Berryman, J. Imaging and time reversal in random media. Inverse Prob. 2002, 18, 1247–1279. [Google Scholar] [CrossRef]
Tsogka, C.; Papanicolaou, G. Time reversal through a solid-liquid interface and super-resolution. Inverse Prob. 2002, 18, 1639–1657. [Google Scholar] [CrossRef]
Blomgren, P.; Papanicolaou, G.; Tsogka, C.; Berryman, J. Super resolution in time-reversal acoustics. J. Acoust. Soc. Am. 2002, 111, 230–48. [Google Scholar] [CrossRef] [PubMed]
Chan, T.; Jaruwatanadilok, S.; Kuga, Y.; Ishimaru, A. Numerical study of the time-reversal effects on super-resolution in random scattering media and comparison with an analytical mode. Waves in Random and Complex Media 2008, 18, 627–639. [Google Scholar] [CrossRef]
Montaldo, G.; Roux, P.; Derode, A.; Negreira, C.; Fink, M. Ultrasonic shock wave generator using 1-bit time-reversal in a dispersive medium : application to lithotripsy. Appl. Phys. Lett. 2002, 80, 897–899. [Google Scholar] [CrossRef]
Thomas, J.L.; Wu, F.; Fink, M. Time reversal focusing applied to lithotripsy. Ultrason. Imag. 1996, 18, 106–121. [Google Scholar] [CrossRef]
Tanter, M.; Thomas, J.L.; Fink, M. Focusing and steering through absorbing and aberrating layers: Application to ultrasonic propagation through the skull. J. Acoust. Soc. Am. 1998, 103, 2403–2410. [Google Scholar] [CrossRef] [PubMed]
Floch, C.L.; Tanter, M.; Fink, M. Self defocusing in ultrasonic hyperthermia: Experiment and simulation. Appl. Phys. Lett. 1999, 74, 3062–3064. [Google Scholar] [CrossRef]
Persson, M.; Trefna, H. Method and system relating to hyperthermia. World Intellectual Property Organization - International Patent Application Publication Application number: PCT/SE2007/000502, Publication No: WO/2007/136335, 2007. [Google Scholar]
Pepper, D.M.; Dunning, G.J.; Sumida, D.S. Time-reversed photoacoustic system and uses thereof. US Patent 6,973,830, 2005. [Google Scholar]
Cerwin, S.A.; Chang, D.B. Identifying and treating bodily tissues using electromagnetically induced, time-reversed, acoustic signals. United States Patent Application Publication Application number: 11/450,097, US 2007/0038060 A1, 2007. [Google Scholar]
Ihn, J.-B.; Dunne, J.P. Virtual time reversal acoustics for structural health monitoring. United States Patent Application Publication Application number: 11/861,056, US 2009/0083004 A1, 2007. [Google Scholar]
Larmat, C.; Tromp, T.; Liu, Q.; Montagner, J.-P. Time reversal location of glacial earthquakes. J. Geophys. Res. B 2008, 113, B09314. [Google Scholar] [CrossRef]
Micolau, G.; Saillard, M.; Borderies, P. DORT method as applied to ultrawideband signals for detection of buried objects. IEEE Trans. Geosci. Remote Sens. 2003, 41, 1813–1820. [Google Scholar] [CrossRef]
Edelmann, G.F.; Song, H.C.; Kim, S.; Akal, T.; Hodgkiss, W.S.; Kuperman, W.A. Underwater acoustic communications using time reversal. IEEE J.Oceanic Eng. 2005, 30, 852–864. [Google Scholar] [CrossRef]
Higley, W.J.; Roux, P.; Kuperman, W.A.; Hodgkiss, W.S.; Song, H.C.; Akal, T.; Stevenson, M. Synthetic aperture time-reversal communications in shallow water: Experimental demonstration at sea. J. Acoust. Soc. Am. 2005, 118, 2365–2372. [Google Scholar] [CrossRef]
Kim, S.; Kuperman, W.A.; Hodgkiss, W.S.; Song, H.C.; Edelman, G.F.; Akal, T. Echo-to-reverberation enhancement using a time reversal mirror. J. Acoust. Soc. Am. 2004, 115, 1525–1531. [Google Scholar] [CrossRef]
Song, H.C.; Hodgkiss, W.S.; Kuperman, W.A.; Roux, P.; Akal, T.; Stevenson, M. Experimental demonstration of adaptive reverberation nulling using time reversal. J. Acoust. Soc. Am. 2005, 118, 1382–1387. [Google Scholar] [CrossRef]
Milsap, J.P. Phased array sound system. US Patent 7,130,430, 2006. [Google Scholar]
Leutenegger, T.; Dual, J. Non-destructive testing of tubes using a time reverse numerical simulation (TRNS) method. Ultrasonics 2004, 41, 811–822. [Google Scholar] [CrossRef] [PubMed]
Kerbrat, E.; Prada, C.; Cassereau, D. Ultrasonic nondestructive testing of scattering media using the decomposition of the time-reversal operator. IEEE T. Ultrason. Ferroelectr. 2002, 49, 1103–1112. [Google Scholar] [CrossRef]
Leutenegger, T.; Dual, J. Detection of defects in cylindrical structures using a time reverse method and a finite-difference approach. Ultrasonics 2002, 40, 721–725. [Google Scholar] [CrossRef]
Chakroun, N.; Fink, M.; Wu, F. Time reversal processing in ultrasonic nondestructive testing. IEEE T. Ultrason. Ferroelectr. 1995, 42, 1087–1098. [Google Scholar] [CrossRef]
Sohn, H.; Kim, S. Methods, apparatuses, and systems for damage detection. World Intellectual Property Organization - International Patent Application Publication Application number: PCT/US2007/015090, Publication No: WO/2008/005311, 2008. [Google Scholar]
Lev-Ari, H.; Devaney, A.J. The Time Reversal Techniques Re-interpreted: Subspace-based Signal Processing for Multistatic Target Location. In Proceedings of IEEE Sensor Array and Multichannel Signal Processing Workshop, Cambridge, MA, USA; 2000; pp. 509–513. [Google Scholar]
Hou, S.; Solna, K.; Zhao, H. Imaging of location and geometry for extended targets using the response matrix. J. Comp. Phys. 2004, 199, 317–338. [Google Scholar] [CrossRef]
Borcea, L.; Papanicolaou, G.; Tsogka, C. Theory and applications of time reversal and interferometric imaging. Inverse Prob. 2003, 19, 1–26. [Google Scholar] [CrossRef]
Borcea, L.; Papanicolaou, G.; Tsogka, C. Interferometric array imaging in clutter. Inverse Prob. 2005, 21, 1419–1460. [Google Scholar] [CrossRef]
Slob, E. Interferometry by deconvolution of multicomponent multioffset gpr data. IEEE Trans. Geosci. Remote Sens. 2009, 47, 828–838. [Google Scholar] [CrossRef]
Derveaux, G.; Papanicolaou, G.; Tsogka, C. Time reversal imaging for sensor networks with optimal compensation in time. J. Acoust. Soc. Am. 2007, 121, 2071–2085. [Google Scholar] [CrossRef] [PubMed]
Prada, C.; Fink, M. Eigenmodes of the time reversal operator: A solution to selective focusing in multiple-target media. Wave Motion 1994, 20, 151–163. [Google Scholar] [CrossRef]
Prada, C.; Mannevile, S.; Spoliansky, D.; Fink, M. Decomposition of the time reversal operator: Detection and selective focusing on two scatterers. J. Acoust. Soc. Am. 1996, 99, 2067–2076. [Google Scholar] [CrossRef]
Prada, C.; Thomas, J.L. Experimental subwavelength localization of scatterers by decomposition of the time reversal operator interpreted as a covariance matrix. J. Acoust. Soc. Am. 2003, 114, 235–243. [Google Scholar] [CrossRef] [PubMed]
Micolau, G.; Saillard, M. DORT method as applied to electromagnetic sensing of buried objects. Radio Sci. 2003, 38, 4–1–12. [Google Scholar] [CrossRef]
Montaldo, G.; Tanter, M.; Fink, M. Real time inverse filter focusing through iterative time reversal. J. Acoust. Soc. Am. 2004, 115, 768–775. [Google Scholar] [CrossRef] [PubMed]
Montaldo, G.; Tanter, M.; Fink, M. Revisiting iterative time reversal processing: Application to detection of multiple targets. J. Acoust. Soc. Am. 2004, 115, 776–784. [Google Scholar] [CrossRef] [PubMed]
Devaney, A.J. Super-resolution processing of multistatic data using time reversal and MUSIC. Unpublished paper, preprint available on the author’s web site. Avalable Online: http://www.ece.neu.edu/faculty/devaney/ (accessed on August 20, 2009).
Gruber, F.K.; Marengo, E.A.; Devaney, A.J. Time-reversal imaging with multiple signal classification considering multiple scattering between the targets. J. Acoust. Soc. Am. 2004, 115, 3042–3047. [Google Scholar] [CrossRef]
Marengo, E.A.; Gruber, F.K. Subspace-based localization and inverse scattering of multiply scattering point targets. EURASIP J. Adv. Signal Proc. 2007. [Google Scholar] [CrossRef]
Lehman, S.K.; Devaney, A.J. Transmission mode time-reversal superresolution imaging. J. Acoust. Soc. Am. 2003, 113, 2742–2753. [Google Scholar] [CrossRef] [PubMed]
Devaney, A.J.; Marengo, E.A.; Gruber, F.K. Time-reversal-based imaging and inverse scattering of multiply scattering point targets. J. Acoust. Soc. Am. 2005, 118, 3129–3138. [Google Scholar] [CrossRef]
Devaney, A.J. Time reversal imaging of obscured targets from multistatic data. IEEE Trans. Antennas Propagat. 2005, 53, 1600–1610. [Google Scholar] [CrossRef]
Pepper, D.M. Nonlinear optical phase conjugation. Optical Engineer. 1982, 21, 156. [Google Scholar]
Byren, R.W.; Rockwell, D.A. Self-aligning phase conjugate laser. US Patent 4,812,639, 1989. [Google Scholar]
Byren, R.W.; Filgas, D. Phase conjugate relay mirror apparatus for high energy laser system and method. US Patent 6,961,171, 2005. [Google Scholar]
Skolnik, M.I.; King, D.D. Self-Phasing array antennas. IEEE Trans. Antennas Propagat. 1964, 12, 142–149. [Google Scholar] [CrossRef]
Sichelstiel, B.A.; Waters, W.M.; Wild, T.A. Self-Focusing array research model. IEEE Trans. Antennas Propagat. 1964, 12, 150–154. [Google Scholar] [CrossRef]
Fusco, V.F.; Karode, S.L. Self-Focusing antenna array techniques for mobile communications applications. Elect. & Comm. Eng. J. 1999, 11, 279–286. [Google Scholar]
Chang, Y.; Fetterman, H.R.; Newberg, I.L.; Panaretos, S.K. Microwave phase conjugation using antenna arrays. IEEE Trans. Microwave Theory Tech. 1998, 46, 1910–1919. [Google Scholar] [CrossRef]
Miyamoto, R.Y.; Itoh, T. Retrodirective arrays for wireless communications. IEEE Microwave Magazine 2002, 3, 71–79. [Google Scholar] [CrossRef]
Henty, B.E.; Stancil, D.D. Multipath-enabled super-resolution for rf and microwave communication using phase-conjugate arrays. Phys. Rev. Lett. 2004, 93, 243904/1–4. [Google Scholar] [CrossRef] [PubMed]
Lerosey, G.; de Rosny, J.; Tourin, A.; Derode, A.; Montaldo, G.; Fink, M. Time reversal of electromagnetic waves. Phys. Rev. Lett. 2004, 92, 193904. [Google Scholar] [CrossRef] [PubMed]
Lerosey, G.; de Rosny, J.; Tourin, A.; Derode, A.; Fink, M. Time reversal of wideband microwaves. Applied Phys. Lett. 2006, 88. [Google Scholar] [CrossRef]
Cepni, A.G. Experimental Investigation of Time-Reversal Techniques using Electromagnetic Waves. Ph D. Dissertation, Carnegie Mellon University, Pittsburgh, PA, USA, 2005. [Google Scholar]
Nguyen, H.T.; Andersen, J.B.; Pedersen, G.F.; Kyritsi, P.; Eggers, P.C.F. Time reversal in wireless communications: A measurement-based investigation. IEEE Trans. Wireless Comm. 2006, 5, 2242–2252. [Google Scholar] [CrossRef]
Moura, J.M.F.; Jin, Y. Detection by time reversal: Single antenna. IEEE Trans. Signal Process. 2007, 55, 187–201. [Google Scholar] [CrossRef]
Moura, J.M.F.; Jin, Y.W. Time reversal imaging by adaptive interference canceling. IEEE Trans. Signal Process. 2008, 56, 233–247. [Google Scholar] [CrossRef]
Jin, Y.; Moura, J.M.F. Time-reversal detection using antenna arrays. IEEE Trans. Signal Process. 2009, 57, 1396–1414. [Google Scholar]
Li, X.; Bond, E.J.; Veen, B.D.V.; Hagness, S.C. An overview of ultra-wideband microwave imaging via space-time beamforming for early-stage breast-cancer detection. IEEE Antennas Propagat. Mag. 2005, 47, 19–34. [Google Scholar]
Kosmas, P.; Rappaport, C.M. Time reversal with the FDTD method for microwave breast cancer detection. IEEE Trans. Microwave Theory Tech. 2005, 53, 2317–2323. [Google Scholar] [CrossRef]
Kosmas, P. FDTD modeling for forward and linear inverse electromagnetic problems in lossy, dispersive media. Ph. D. Dissertation, Northeastern University, Boston, MA, USA, 2005. [Google Scholar]
Kosmas, P.; Rappaport, C.M. FDTD-based time reversal for microwave breast cancer detection - localization in three dimensions. IEEE Trans. Microwave Theory Tech. 2006, 54, 1921–1927. [Google Scholar] [CrossRef]
Kosmas, P.; Rappaport, C.M. A matched-filter FDTD-based time reversal approach for microwave breast cancer detection. IEEE Trans. Antennas Propagat. 2006, 54, 1257–1264. [Google Scholar] [CrossRef]
Chen, Y.F.; Gunawan, E.; Low, K.S.; Wang, S.C.; Kim, Y.; Soh, C.B. Pulse design for time reversal method as applied to ultrawideband microwave breast cancer detection: A two-dimensional analysis. IEEE Trans. Antennas Propagat. 2007, 55, 194–204. [Google Scholar] [CrossRef]
Chen, Y.F.; Gunawan, E.; Low, K.S.; Wang, S.C.; Soh, C.B.; Putti, T.C. Time-reversal ultrawideband breast imaging: Pulse design criteria considering multiple tumors with unknown tissue properties. IEEE Trans. Antennas Propagat. 2008, 56, 3073–3077. [Google Scholar] [CrossRef]
Jin, Y.; Moura, J.M.F.; Jiang, Y.; Wahl, M.; Zhu, H.; He, Q. Breast cancer detection by time reversal imaging. In Proceedings of 5th IEEE International Symposium on Biomedical Imaging: From Nano to Macro (ISBI), Paris, France, May 14-17, 2008; pp. 816–819.
Liu, D.H.; Kang, G.; Li, L.; Chen, Y.; Vasudevan, S.; Joines, W.; Liu, Q.H.; Krolik, J.; Carin, L. Electromagnetic time-reversal imaging of a target in a cluttered environment. IEEE Trans. Antennas Propagat. 2005, 53, 3058–3066. [Google Scholar]
Liu, D.H.; Vasudevan, S.; Krolik, J.; Bal, G.; Carin, L. Electromagnetic time-reversal source localization in changing media: Experiment and analysis. IEEE Trans. Antennas Propagat. 2007, 55, 344–354. [Google Scholar] [CrossRef]
Liu, D.H.; Krolik, J.; Carin, L. Electromagnetic target detection in uncertain media: Time-reversal and minimum-variance algorithms. IEEE Trans. Geosci. Remote Sens. 2007, 45, 934–944. [Google Scholar] [CrossRef]
Tortel, H.; Micolau, G.; Saillard, M. Decomposition of the time reversal operator for electromagnetic scattering. J. Electromagn. Waves Appl. 1999, 13, 687–719. [Google Scholar] [CrossRef]
Chambers, D.H.; Berryman, J.G. Analysis of the time-reversal operator for a small spherical scatterer in an electromagnetic field. IEEE Trans. Antennas Propagat. 2004, 52, 1729–1738. [Google Scholar] [CrossRef]
Berryman, J.G. Method for distinguishing multiple targets using time-reversal acoustics. US Patent 6,755,083, 2004. [Google Scholar]
Sarabandi, K.; Koh, I.; Casciato, M.D. Demonstration of time reversal methods in a multi-path environment. Proceedings of the IEEE Antennas and Propagation Society International Symposium 2004, 4, 4436–4439. [Google Scholar]
Ziade, Y.; Roussel, H.; Lesturgie, M.; Tabbara, W. A coherent model of forest propagation–application to detection and localization of targets using the DORT method. IEEE Trans. Antennas Propagat. 2008, 56, 1048–1057. [Google Scholar] [CrossRef]
Jiang, Y.; Zhu, J.; Stancil, D.D.; Chabalko, M.J. Polarization sensitive time-reversal SAR imaging in an environment filled with trees. Proc. of the IEEE Antennas Propagat. Symp. 2007, 1, 4000–4003. [Google Scholar]
Jiang, Y.; Zhu, J.; Stancil, D.D. Synthetic aperture radar ghost image cancellation using broadband time reversal averaging techniques. IEEE/MTT-S Int. Microwave Symposium 2007, 1479–1482. [Google Scholar] [CrossRef]
Jin, Y.; Moura, J.M.F. Application of time reversal to synthetic aperture imaging. United States Patent Application Publication Application number: 12/217,839, US 2009/0033549 A1, 2008. [Google Scholar]
Jin, Y.; Moura, J.M.F. Imaging by time reversal beamforming. United States Patent Application Publication Application number: 12/217,838, US 2009/0076389 A1, 2009. [Google Scholar]
Zheng, W.; Zhao, Z.; Nie, Z.P. Application of TRM in the UWB through wall radar. Prog. Electromagnetics Res. - PIER 2008, 87, 279–296. [Google Scholar] [CrossRef]
Kyritsi, P.; Papanicolaou, G.; Eggers, P.; Oprea, A. MISO time reversal and delay-spread compression for FWA channels at 5 GHz. IEEE Antennas Wirel. Prop. Lett. 2004, 3, 96–99. [Google Scholar] [CrossRef]
Lerosey, G.; de Rosny, J.; Tourin, A.; Derode, A.; Montaldo, G.; Fink, M. Time reversal of electromagnetic waves and telecommunication. Radio Sci. 2005, 40. [Google Scholar] [CrossRef]
Montaldo, G.; Lerosey, G.; Derode, A.; Tourin, A.; de Rosny, J.; Fink, M. Telecommunication in a disordered environment with iterative time reversal. Waves in Random Media 2004, 14, 287–302. [Google Scholar] [CrossRef]
Xiao, S.Q.; Chen, J.; Wang, B.Z.; Liu, X.F. A numerical study on time-reversal electromagnetic wave for indoor ultra-wideband signal transmission. Prog. Electromagnetics Res. - PIER 2007, 77, 329–342. [Google Scholar] [CrossRef]
Blomgren, P.; Kyritsi, P.; Kim, A.D.; Papanicolaou, G. Spatial focusing and intersymbol interference in multiple-input-single-output time reversal communication systems. IEEE J. Oceanic Eng. 2008, 33, 341–355. [Google Scholar] [CrossRef]
Lindskog, E.D.; Paulraj, A.J. Time-reversal block transmit diversity system for channels with intersymbol interference and method. US Patent 7,362,815, 2008. [Google Scholar]
Candy, J.V.; Chambers, D.H.; Guidry, B.L.; Poggio, A.J.; Robbins, C.L. Multi-channel time-reversal receivers for multi and 1-bit implementations. US Patent 7,463,690, 2008. [Google Scholar]
Pajusco, P.; Pagani, P. On the use of uniform circular arrays for characterizing UWB time reversal. IEEE Trans. Antennas Propagat. 2009, 57, 102–109. [Google Scholar] [CrossRef]
Coppinger, F.; Bhushan, A.S.; Jalali, B. Time reversal of broadband microwave signals. Electron. Lett. 1999, 35, 1230–1232. [Google Scholar] [CrossRef]
Yanik, M.F.; Suh, W.; Wang, Z.; Fan, S. Stopping and time reversing light in a waveguide with an all-optical system. US Patent 7,116,864, 2006. [Google Scholar]
Yavuz, M.E.; Teixeira, F.L. A numerical study of time reversed UWB electromagnetic waves in continuous random media. IEEE Antennas Wirel.Prop. Lett. 2005, 4, 43–46. [Google Scholar] [CrossRef]
Yavuz, M.E.; Teixeira, F.L. Full time-domain DORT for ultrawideband fields in dispersive, random inhomogeneous media. IEEE Trans. Antennas Propagat. 2006, 54, 2305–2315. [Google Scholar] [CrossRef]
Yavuz, M.E.; Teixeira, F.L. Space-frequency ultrawideband time-reversal imaging. IEEE Trans. Geosci. Remote Sens. 2008, 46, 1115–1124. [Google Scholar] [CrossRef]
Yavuz, M.E.; Teixeira, F.L. On the sensitivity of time-reversal imaging techniques to model perturbations. IEEE Trans. Antennas Propagat. 2008, 56, 834–843. [Google Scholar] [CrossRef]
Yavuz, M.E.; Teixeira, F.L. Frequency dispersion compensation in time reversal techniques for UWB electromagnetic waves. IEEE Geosci. Remote Sensing Lett. 2005, 2, 233–237. [Google Scholar] [CrossRef]
Born, M.; Wolf, E. Principles of Optics; Academic Press: New York, NY, USA, 1970. [Google Scholar]
Fung, A.K.; Fung, H.S. Scattering from a vegetation layer. IEEE Trans. Geosci. Electron. 1979, 17, 1–6. [Google Scholar] [CrossRef]
Tsang, L.; Kong, J.A. Thermal microwave emission from a three-layer random medium with three dimensional variations. IEEE Trans. Geosci. Remote Sens. 1980, 18, 212–216. [Google Scholar] [CrossRef]
Wang, J.R.; Shiue, J.C.; Chuang, S.L.; Shin, R.T.; Dobrowski, M. Thermal microwave emission from vegetated fields: A comparison between theory and experiment. IEEE Trans. Geosci. Remote Sens. 1984, 22, 143–149. [Google Scholar] [CrossRef]
Moss, C.D.; Teixeira, F.L.; Yang, Y.E.; Kong, J.A. Finite-difference time-domain simulation of scattering from objects in continuous random media. IEEE Trans. Geosci. Remote Sens. 2002, 40, 178–186. [Google Scholar] [CrossRef]
Berenger, J.P. A perfectly matched layer for the absorption of electromagnetic waves. J. Comput. Phys. 1994, 114, 185–200. [Google Scholar] [CrossRef]
Teixeira, F.L.; Chew, W.C. A general approach to extended berenger’s absorbing boundary condition to anisotropic and dispersive media. IEEE Trans. Antennas Propagat. 1998, 46, 1386–1387. [Google Scholar] [CrossRef]
Harris, F.J. On the use of windows for harmonic analysis with the discrete fourier transform. Proc. IEEE 1978, 66, 51–83. [Google Scholar] [CrossRef]
Prada, C.; Wu, F.; Fink, M. The iterative time reversal mirror: A solution to self-focusing in the pulse echo mode. J. Acoust. Soc. Am. 1991, 90, 1119–1129. [Google Scholar] [CrossRef]
Yavuz, M.E. Time reversal based signal processing techniques for ultrawideband electromagnetic sensing in random media. Ph. D. Dissertation, The Ohio State University, Columbus, OH, USA, 2008. [Google Scholar]
Yavuz, M.E.; Teixeira, F.L. Time Reversal Based Signal Processing Techniques: Algorithms and Applications on Ultrawideband Electromagnetic Remote Sensing; VDM Verlag: Saarbrucken, Germany, 2009. [Google Scholar]
Scholz, B. Towards virtual electrical breast biopsy: Space-frequency MUSIC for trans-admittance data. IEEE Trans. Med. Imag. 2002, 21, 588–595. [Google Scholar] [CrossRef] [PubMed]
Wang, Z.S.; Li, J.; Wu, R.B. Time-delay- and time-reversal-based robust capon beamformers for ultrasound imaging. IEEE Tran. Medical Imag. 2005, 24, 1308–1322. [Google Scholar] [CrossRef] [PubMed]
Nguyen, D.Q.; Wan, W.S. Novel DORT method in non-well-resolved scatterer case. IEEE Signal Proc. Lett. 2008, 15, 705–708. [Google Scholar] [CrossRef]
Zhao, H. Analysis of the response matrix for an extended target. SIAM J. Appl. Math. 2004, 64, 725–745. [Google Scholar] [CrossRef]
Yavuz, M.E.; Teixeira, F.L. Effects of array restriction on time-reversal methods. Proc. of the 2007 IEEE Antennas Prop. Symp. 2007. [Google Scholar] [CrossRef]
Marengo, E.A.; Gruber, F.K.; Simonetti, F. Time-reversal MUSIC imaging of extended targets. IEEE Trans. Image Processing 2007, 16, 1967–1984. [Google Scholar] [CrossRef]
Baussard, A.; Boutin, T. Time-reversal RAP-MUSIC imaging. Waves in Random and Complex Media 2008, 18, 151–160. [Google Scholar] [CrossRef]
Folégot, T.; Prada, C.; Fink, M. Resolution enhancement and separation of reverbaration from target echo with the time reversal operator decomposition. J. Acoust. Soc. Am. 2003, 113, 3155–3160. [Google Scholar] [CrossRef] [PubMed]
Saillard, M.; Vincent, P.; Micolau, G. Reconstruction of buried objects surrounded by small inhomogeneities. Inverse Prob. 2000, 16, 1195–1208. [Google Scholar] [CrossRef]
Pozar, D.M. Microwave Engineering, 2nd ed.; John Wiley & Sons: New York, NY, USA, 1998. [Google Scholar]
Teixeira, F.L.; Chew, W.C.; Straka, M.; Oristaglio, M.L.; Wang, T. Finite-difference time-domain simulation of ground penetrating radar on dispersive, inhomogeneous and conductive soils. IEEE Trans. Geosci. Remote Sensing 1998, 36, 1928–1937. [Google Scholar] [CrossRef]
Oppenheim, A.; Schafer, R.W. Discrete-Time Signal Processing; Prentice Hall: Upper Saddle River, NJ, USA, 1999. [Google Scholar]

© 2009 by the authors; licensee Molecular Diversity Preservation International, Basel, Switzerland. This article is an open-access article distributed under the terms and conditions of the Creative Commons Attribution license ( http://creativecommons.org/licenses/by/3.0/).

Share and Cite

MDPI and ACS Style

Yavuz, M.E.; Teixeira, F.L. Ultrawideband Microwave Sensing and Imaging Using Time-Reversal Techniques: A Review. Remote Sens. 2009, 1, 466-495. https://doi.org/10.3390/rs1030466

AMA Style

Yavuz ME, Teixeira FL. Ultrawideband Microwave Sensing and Imaging Using Time-Reversal Techniques: A Review. Remote Sensing. 2009; 1(3):466-495. https://doi.org/10.3390/rs1030466

Chicago/Turabian Style

Yavuz, Mehmet Emre, and Fernando L. Teixeira. 2009. "Ultrawideband Microwave Sensing and Imaging Using Time-Reversal Techniques: A Review" Remote Sensing 1, no. 3: 466-495. https://doi.org/10.3390/rs1030466

APA Style

Yavuz, M. E., & Teixeira, F. L. (2009). Ultrawideband Microwave Sensing and Imaging Using Time-Reversal Techniques: A Review. Remote Sensing, 1(3), 466-495. https://doi.org/10.3390/rs1030466

Article Menu

Ultrawideband Microwave Sensing and Imaging Using Time-Reversal Techniques: A Review

Abstract

1. Introduction

2. Time-Reversal of Electromagnetic Waves and Superresolution

2.1. Physics-Based Ultrawideband Clutter Models

2.2. Numerical Results

3. Time-Reversal based Detection and Imaging Methods

3.1. Signal Space Methods

3.1.1. Time-Domain DORT

3.1.2. Space-Frequency TR Imaging

3.1.3. Results

3.2. Null Space Methods

3.2.1. TR-MUSIC

4. Frequency Dispersion/Loss Compensation Techniques

5. Conclusions

Acknowledgements

References and Notes

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI