High-Resolution Imaging Utilizing Space-Frequency DORT Combining the Extrapolated Virtual Array

Space-frequency DORT (SF-DORT) is an effective timereversal (TR) imaging method due to its immunity to the noise and adaptability of complex environment. However, some potential drawbacks, such as low range and co-range resolutions, make SFDORT inferior to the Space-frequency MUSIC (SF-MUSIC). In this paper, we propose a novel high-resolution imaging method utilizing SF-DORT combining an extrapolated virtual array (EASF-DORT). The virtual extrapolated TR array is created by autoregressive vector extrapolation (ARVE). In this way, the aperture of the TR array is extended significantly and the drawbacks associated with SF-DORT can be improved. Simulation results demonstrate that SF-DORT combining the extrapolated virtual array achieve much higher range and co-range resolutions than conventional SF-DORT.


I. INTRODUCTION
Time reversal (TR) technique has attracted enormous attention in radar imaging fields in recent years due to its spatial and temporal high-resolution focusing characteristics [1].The TR operator (TRO) is one of the most fundamental TR imaging method, and TRO form the basis of the decomposition of the time-reversal operator (DORT) and the TR multiple signal classification (TR-MUSIC) [2]- [4].Space-frequency decomposition of the time-reversal operator (SF-DORT) and SF multiple signal classification (SF-MUSIC) are faster TR imaging methods because they obtain the multistatic data matrices (MDMs) much more conveniently than DORT and TR-MUSIC [5]- [8].
Although SF-DORT works well under increased clutter and noise levels, its range and co-range resolutions are inferior to SF-MUSIC [9].
To overcome the shortcomings of SF-DORT, utilizing virtual extrapolated array is a promising solution because it is suitable for an array with limited aperture and meets the requirement of miniaturizing the imaging system [10].The autoregressive vector extrapolation (ARVE) [11], which exploits the spatial-temporal coupling of the matrix constituted by the received signals, predicts the signals on the virtual antennas more accurately than the 1-D AR extrapolation in [12]- [13].
Based on the above analysis, we propose a novel high-resolution imaging method utilizing SF-High-Resolution Imaging Utilizing Space-Frequency DORT Combining the Extrapolated Virtual Array DORT combining the virtual extrapolated array (EA-SF-DORT) in this article.The performance deterioration of imaging caused by applying the ARVE with the noise matrix can be offset by SF-DORT, and high-resolution can be achieved.Simulations show that SF-DORT utilizing the extrapolated virtual array can achieve much higher range and co-range resolutions than without utilizing.
The rest of this paper is organized as follows.Section II shows the whole processing chain of the proposed TR imaging method, including the principle of ARVE and SF-DORT.To illustrate the highresolution imaging capability, computer simulations are conducted, and the results are given in Section III.Finally, Section IV offers some conclusions drawn on the basis of simulation results.

II. THE PRINCIPLE OF PROPOSED TR IMAGING METHOD
The processing chain of this novel method is shown in Fig. 1 and could be divided into two stages.
In the first stage, the matrix constituted by the received noised-  A. The principle of ARVE Assume a TR array, which contains L real antennas, lies parallel to the x-axis.Assume the TR array receives the same linear frequency modulated (LFM) signals, so the received signals are coherent.
The received signal on the l-th real antenna is given below where ml  is the distance from the m-th target to the l-th antenna.m is the target index and With the original measurement data S in hand, the extrapolated data E S is formulated through ARVE.The principle of ARVE is shown in Fig. 2 and can be divided into three steps.The first step is to determine the number of parameters of S , i.e. select the order of the AR model.
The number of parameters is denoted by a pair ( 1 p , 2 p ), where 1 p is the order of the spatial dimension and 2 p is the order of the temporal dimension.Both 1 p and 2 p are positive integers and can be determined by S itself.This step can be done by choosing the typical value of the model order directly or using an order estimation algorithm such as MDL [14].The typical range for 1 p is between 1 and 2 L .The typical range for 2 p is between 1 and 2 N .The pair (1, 1) should be excluded [10].
With the number of parameters determined, the next step is estimating the block AR model parameters.We will adopt the approach 2-D lattice linear prediction parameter estimation in [15] because it gives a recursive solution for the 2-D AR parameter matrices and it is a fast algorithm.
In the last stage, we provide a brief summary of the vector extrapolation equations.Given a column vector of size 2 ( 1) 1 p  as follows and the forward linear vector extrapolation is given in (5) The backward linear vector extrapolation is given by where * represents the conjugate operation.
For more details about vector extrapolation equations, please see [7] and [11].
Therefore, the matrix S can be extrapolated to the matrix E S as follows: (0) where is the spatial extrapolation factor of the array.The original data set S is the central part of E S .U S and D S are extended data from the top and bottom of the matrix S respectively.

B. TR imaging based on SF-DORT
After having obtained the extrapolated data E S , the SF-MDM K can be established by applying Fourier transform to each row of E S .The matrix K has size LN   and can be expressed as below Appling SVD to K, we have  KU ΛV (10) where Λ is a N × M real diagonal matrix which contains singular values.U is an N × N unitary matrix, and V is the M × M matrix.Since the matrix K includes the spatial and frequency of the target,

III. SIMULATIONS
In this section, we compare the range and co-range resolutions utilizing conventional SF-DORT and the proposed EA-SF-DORT for two different cases.The first is a single target case while the second is the multi-target one.We utilize SF-DORT to evaluate the targets detection performance.
The received LFM signal can be expressed in (2).The parameters of the LFM signal are as follows: B=0.5GHz , N=4000, 4.5GHz . The signal to noise ratio (SNR) is 10 dB.

A. Imaging for a single target case
Consider the single point target located at (62.5m, 75m).All the antennas in the TR array receives LFM signal.
The range and co-range resolutions given by SF-DORT with and without autoregressive vector

B. Imaging for a multi-target case
In this simulation, we consider a two targets case when the targets are located at (75m,62.5m) and (25m, 37.5m) respectively.All the antennas in the TR array receive the same LFM signal.The extrapolation factor α is 4 and the extrapolated matrix E S has the size of (16×4000).Once the left singular vectors p U , p = 1, . . ., P , P = 2 have been acquired using (10), we could locate the two targets with the same procedure above.
The comparison of range and co-range resolutions for the target located at (75m, 62.5m) given by SF-DORT with and without ARVE are shown in Fig. 5, while the same comparison for the target located at (25m, 37.5m) are shown in Fig. 6.From Fig. 5 and Fig. 6 the range and co-range resolutions are much higher when extrapolation is used, especially the range resolution for the target at (25m, 37.5m).

IV. CONCLUSION
In this paper, we have presented a novel high-resolution TR imaging method utilizing SF-DORT combining virtual array created by ARVE (EA-SF-DORT).The ARVE extends the aperture of the TR array significantly, consequently overcomes the shortcomings of conventional SF-DORT.It has been shown that the SF-DORT utilizing extrapolated virtual array achieves much higher range and corange resolutions than conventional SF-DORT.Simulation results validate the satisfying performance of the proposed method.

Fig. 1 .
Fig. 1.The processing chain of the new TR imaging method.

Fig. 2 .
Fig. 2. Virtual measurement data are created by modeling the original measurements with a 2-D AR model and then using vector extrapolation

2 pJ
is a reflection matrix of dimension 22 ( 1) ( 1) pp    whose anti-diagonal elements are one and all other elements are zero, i.e., " denotes transposition.By applying inverse Fourier transform to each row of TR ()  R , we could get the time-domain TR signal TR () t r .Once the vector TR ()  R or TR () t r is obtained, we can get the following pseudon-th antenna.

extrapolation are shown in Fig 4 .Fig. 4 .
Fig. 4. Range and co-range resolutions utilizing SF-DORT and EA-SF-DORT for a single target and the target is located at (62.5m,75m) (a) Co-range resolutions (b) Range resolutions

Fig. 5 .Fig. 6 .
Fig. 5. Range and co-range resolutions utilizing SF-DORT and EA-SF-DORT for Multitargets case, the target is located at (75m,62.5m).(a) Co-range resolutions (b) Range resolutions SVD to K is called SF decomposition.Assume p U as the p-th column of the matrix U, p Vas the p-th column of the matrix V, and p  as the p-th singular value.p U contains the spatial relationships between the p-th target and the extrapolated array while p V represents the frequency content of the received signals.Assuming S(ω) as the Fourier transform of the transmitted signal s(t), a new N × 1 retransmitted vector function can be defined as below