Next Article in Journal
Influence of Moisture Content on Electromagnetic Response of Concrete Studied Using a Homemade Apparatus
Next Article in Special Issue
A MIMO-SAR Tomography Algorithm Based on Fully-Polarimetric Data
Previous Article in Journal
Slope Stability Monitoring Using Novel Remote Sensing Based Fuzzy Logic
Previous Article in Special Issue
Azimuth Phase Center Adaptive Adjustment upon Reception for High-Resolution Wide-Swath Imaging
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two-Dimensional Augmented State–Space Approach with Applications to Sparse Representation of Radar Signatures

School of Electronics and Information Engineering, Beihang University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(21), 4631; https://doi.org/10.3390/s19214631
Submission received: 16 September 2019 / Revised: 18 October 2019 / Accepted: 21 October 2019 / Published: 24 October 2019
(This article belongs to the Special Issue Recent Advancements in Radar Imaging and Sensing Technology)

Abstract

:
In this work, we focus on sparse representation of two-dimensional (2-D) radar signatures for man-made targets. Based on the damped exponential (DE) model, a 2-D augmented state–space approach (ASSA) is proposed to estimate the parameters of scattering centers on complex man-made targets, i.e., the complex amplitudes and the poles in down-range and aspect dimensions. An augmented state–space approach is developed for pole estimation of down-range dimension. Multiple-range search strategy, which applies one-dimensional (1-D) state–space approach (SSA) to the 1-D data for each down-range cell, is used to alleviate the pole-pairing problem occurring in previous algorithms. Effectiveness of the proposed approach is verified by the numerical and measured inverse synthetic aperture radar (ISAR) data.

1. Introduction

Sparse representation of two-dimensional (2-D) radar signatures has been widely used in many applications, such as super-resolution radar imaging, data compression, and target identification [1,2,3,4]. 2-D radar signatures can be reconstructed with fewer data, where a set of parameters including the locations, amplitudes, and damping factors are used to represent the returned signals from spatial distributed scattering centers. Moreover, 2-D ultra-wideband images of radar targets can be obtained through parameter interpolation or extrapolation [5].
Over the years, several model-based spectral estimation approaches have been developed for sparse representation of 2-D radar signatures, in which the common idea of these approaches is the development of a parametric model based on electromagnetic scattering mechanisms. In this way, a set of parameters may be found to represent the original signatures. Please note that most sparse representation approaches are considered to have super-resolution capabilities because they are supposed to estimate the parameters of scattering centers that cannot be distinguished by standard processing [6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21]. Typically, these approaches include the amplitude and phase estimation of a sinusoid (APES) algorithm [6], fast Fourier transform (FFT)-based technique (CLEAN) [7,8], the compressed sensing (CS)-based procedure [9,10,11] and the subspace-based approach [12,13,14,15,16,17,18,19,20,21]. The APES algorithm could provide accurate estimation of the complex amplitudes whereas it does not estimate the locations of scattering centers. The 2-D CLEAN technique is a computationally efficient procedure which uses the undamped exponential model and deconvolution algorithm for optimization [7,22]. The main process of the CS-based procedure is to represent the signatures using an over-complete dictionary and the corresponding coefficients [9,10]. The performance of this approach depends on the initial coefficients and the threshold used in terminating iterations.
Another important group of 2-D sparse representation technique is the subspace-based approaches. The representative algorithms of this kind include the 2-D multiple signal classification (MUSIC) [12], matrix enhancement and matrix pencil (MEMP) [13,14], 2-D total least squares (TLS) Prony [15], algebraically coupled matrix pencils (ACMP) [16], 2-D estimation of signal parameters via rotational invariance techniques (2-D ESPRIT) [17] and 2-D system realization technique [18,19]. The rationale and pole-pairing scheme of these algorithms are listed in Table 1.
While these subspace-based algorithms are demonstrated to be useful for the signatures of targets or simulated scattering centers, there are still some technique challenges, especially for sparse representation of the wideband radar signatures collected on complex man-made targets. First, the pairing schemes used by many existing algorithms cannot provide correct pole pairs in certain circumstances. i.e., MEMP, ACMP and 2-D system realization meet with the pole-pairing problem when there are numerous repeated poles in either the down-range or aspect dimension [17,20,23]. Second, the model order, which directly affects the result of parameter estimation such as 2-D MUSIC, MEMP, and 2-D ESPRIT et al., is difficult to be determined in the measurement environment [8]. Third, the compromises between the computational complexity and accuracy should be considered in the application of these algorithms [12]; for instance, the 2-D TLS Prony has smaller computational complexity than MEMP but with loss in accuracy [15].
In this paper, we focus on developing an approach for sparse representation of 2-D radar signatures collected on man-made targets. An augmented state–space approach is proposed for pole estimation of the down-range dimension. Multiple-range search strategy is then applied to estimate the pairing poles and corresponding amplitudes along the aspect dimension. Compared to the existing methods [14,17,18], the advantages of the 2-D augmented state–space approach (ASSA) are: (1) Computational complexity of the algorithm is much reduced since the newly defined Hankel matrix and several time-saving operations are adopted, whereas the pole-estimation accuracy is still at the same level; (2) The pole-pairing problem can be alleviated because all the poles are adaptively paired by using the multiple-range search strategy; And (3) an eigenvalue sequences transform algorithm is proposed, which could provide fast model order selection.
The remainder of the paper is organized as follows. Section 2 briefly presents the damped exponential models for 1-D and 2-D signals. In Section 3, a two-step procedure for 2-D ASSA is developed. Results for numerical and measured inverse synthetic aperture radar (ISAR) data representation are demonstrated to validate the effectiveness of the proposed procedure in Section 4. We conclude the paper in Section 5. Appendices are given to show more mathematical details of the proposed approach.

2. Data Model

2.1. 1-D Damped Exponential Model

As described in DE model [22,24], the 1-D radar signatures y f n can be expressed as a summation of K scattering centers corrupted with noise w ( n ) .
y f n = k = 1 K A k exp β k + j 2 π r k c f n + w n = k = 1 K a k p k n + w n
where n = 1,2,…,N and N denotes the number of pulses, Ak is the complex amplitude of the k-th scattering center; βk is the damping factor with respect to frequency; r k denotes the relative range; The parameter fn denotes the radar frequency f n = f c + n 1 N 2 Δ f where f c is the center frequency and N2 = ceil(N/2) denotes the smallest integer less than or equal to N/2; a k represents the amplitude of the k-th scattering center in pole form, the pole p k = exp β k + j 2 π r k / c Δ f represents the transfer function of the k-th scattering center; c = 3 × 108 m/s is the propagation velocity. It worth noting that the data models used in this paper are all considered in a stepped frequency radar [25,26].
According to the discrete-time control theory and auto-regressive moving average (ARMA) model, Piou and Naishadham proposed a one-dimensional state–space approach (1-D SSA) which use a state–space description to the 1-D radar signatures in (1) [22].
x n + 1 = A x n + B u n
y n = C x n + u n
where x n K × 1 is the state vector, y n denotes the signal sequence y f n , u n is the input vector, A K × K represents the open-loop matrix, B K × 1 and C 1 × K are the constant matrices. Thus, the 1-D noiseless radar signatures same as (1) can be expressed as
y ˜ ( 1 ) y ˜ ( 2 ) y ˜ ( N ) = C B C A B C A N 1 B
As described in [22], 1-D SSA could precisely estimate the state matrices A, B, and C. Once these three state matrices are computed, the model parameters in (1) can be estimated by using the eigen-decomposition technique.

2.2. 2-D Damped Exponential Model

As a 2-D extension of 1-D DE model, the radar signatures obtained from different aspect angles are considered to be the summation of a finite number of dispersive scattering centers [13,14,15,16,17,18]. Typically, it is applicable to modeling the 2-D radar signatures with small aspect ranges [27].
y ( θ m , f n ) = k = 1 K A k exp j 2 π r 1 k c + β 1 k f c θ m exp j 2 π r 2 k c + β 2 k f n + w ( m , n ) = k = 1 K a k s k m p k n + w ( m , n )
in matrix notation, the 2-D radar signatures in (5) can be expressed as:
Y = y ( 1 , 1 ) y ( 1 , 2 ) y ( 1 , N ) y ( 2 , 1 ) y ( 2 , 2 ) y ( 2 , N ) y ( M , 1 ) y ( M , 2 ) y ( M , N )
where m = 1 , , M and n = 1 , , N ; { r 1 k , r 2 k } give the relative locations of the k-th scattering center; { β 1 k , β 2 k } characterize the frequency and aspect dependence of scattering; θ m = θ 0 + ( m 1 ) Δ θ denotes the m-th aspect angle where θ 0 is the starting angle and Δ θ represents the angle interval; w ( m , n ) is the Gaussian noise with zero-mean; y ( m , n ) denotes the signal sequences y θ m , f n ; { s k , p k } refer to poles of the down-range and aspect dimension.
s k = exp j 2 π r 1 k c + β 1 k f c Δ θ
p k = exp j 2 π r 2 k c + β 2 k Δ f

3. Two-Dimensional Augmented State–Space Approach

From the 2-D DE model in (5), the vector Y ˜ ( n ) , which represents the n-th column of the noiseless signature matrix Y ˜ = Y W (where W is the noise matrix), can be decomposed as.
Y ˜ ( n ) = a 1 s 1 1 p 1 a 2 s 2 1 p 2 a K s K 1 p K a 1 s 1 2 p 1 a 2 s 2 2 p 2 a K s K 2 p K a 1 s 1 M p 1 a 2 s 2 M p 2 a K s K M p K p 1 0 0 0 p 2 0 0 0 p K n 1 l K
where K denotes the number of scattering centers, l K indicates the column vector of ones with length K.
Actually, we often meet the one-to-multiple matching situation, i.e., one pole p k is corresponding to more than one poles s k 1 , s k 2 , . These repeated poles can be merged and (9) could be rewritten as.
Y ˜ ( n ) = t = 1 l 1 a 1 t s 1 t p 1 t = 1 l 2 a 2 t s 2 t p 2 t = 1 l K 1 a K 1 t s K 1 t p K 1 t = 1 l 1 a 1 t s 1 t 2 p 1 t = 1 l 2 a 2 t s 2 t 2 p 2 t = 1 l K 1 a K 1 t s K 1 t 2 p K 1 t = 1 l 1 a 1 t s 1 t M p 1 t = 1 l 2 a 2 t s 2 t M p 2 t = 1 l K 1 a K 1 t s K 1 t M p K 1 p 1 0 0 0 p 2 0 0 0 p K n 1 l K 1
where K 1 represents the number of non-repeated poles in the down-range dimension. l 1 , l 2 , l K 1 denote the number of repeated poles, respectively, for the corresponding poles p 1 , p 2 , p K 1 .
We define the matrix P as.
P = Q p 1 0 0 0 p 2 0 0 0 p K 1 Q *
where * denotes the Hermitian operator; Q is a unitary matrix which satisfies Q * Q = Q Q * = E , E is the identity matrix.
Two constant matrices, i.e., S and D, are defined to simplify the expression in (10).
S = t = 1 l 1 a 1 t s 1 t p 1 t = 1 l 2 a 2 t s 2 t p 2 t = 1 l K 1 a K 1 t s K 1 t p K 1 t = 1 l 1 a 1 t s 1 t 2 p 1 t = 1 l 2 a 2 t s 2 t 2 p 2 t = 1 l K 1 a K 1 t s K 1 t 2 p K 1 t = 1 l 1 a 1 t s 1 t M p 1 t = 1 l 2 a 2 t s 2 t M p 2 t = 1 l K 1 a K 1 t s K 1 t M p K 1 Q *
D = Q l K 1
Using (11)–(13), expression in (10) is simplified as.
Y ˜ ( n ) = S P n 1 D
Thus, the noiseless signatures Y ˜ can be written as.
Y ˜ = S D S P D S P N 1 D
Comparing (15) with (4), it can be found that these two equations have extremely similar structures. Thus, the parameters a k , s k and p k can be estimated by 2-D augmentation of 1-D SSA, which is called 2-D augmentation state–space approach (2-D ASSA).
Here, the 2-D ASSA consists of two steps. An augmented state–space approach is applied to estimating the pole p k in down-range dimension. Then multiple-range search strategy is used for estimating the matching pole s k as well as the corresponding amplitude a k . Details are as follows.

3.1. Pole Estimation of the Down-Range Dimension

First, we introduce a single augmented Hankel matrix H of size M ( N L + 1 ) × L by analogy of the enhanced Hankel matrix in 1-D SSA:
H = Y 1 Y 2 Y L Y 2 Y 3 Y L + 1 Y N L + 1 Y N L + 2 Y N
where each element Y ( n ) , first mentioned in (9), represents the n-th column of the matrix Y; L is the step size of the correlation window, which is heuristically set to be L= N 2 [22].
Next, do the singular value decomposition (SVD) of the Hankel matrix.
H = [ U s n U n ] s n n V s n * V n *
where U s n , U n , V s n * and V n * are unitary matrices; s n and n are diagonal matrices; The matrices with subscript ‘sn’ refer to the components of signal space and those matrices with subscript ‘n’ denote the components of noise space; The rank of the noiseless matrix U s n s n V s n * is known as the model order which has been widely used in [14,22,27].
As proved in [14,17], the model order in the down-range dimension is equal to K 1 and should be less than the number of columns or rows of H at least. Thus, the number of non-repeated poles K 1 should satisfy the following condition.
M N M L + M K 1 L K 1
Model order selection is an inevitable problem in modeling the 2-D signatures [14,22,27]. A series of eigenvalue-based criteria, such as Akaike information criterion (AIC), minimum description length (MDL) criterion and the minimum eigenvalue (MEV) criterion, have been proposed for solving this problem [28,29]. Those criteria are useful in model order selection but with time consuming process [30]. Here an eigenvalue sequences transform algorithm is proposed for estimating the model order due to its low computational complexity. This algorithm can provide fast estimation but with loss in robustness to noise. More details of this algorithm are presented in Appendix A.
Based on the linear systems theory [31] and the matrix expression in (14) and (16), the noiseless matrix H ˜ can be further factorized as:
H ˜ = U s n s n V s n * = Ω Γ
where Ω is the ( N L + 1 ) M × K 1 observability matrix which can be further expressed by the matrices S and P.
Ω = U s n s n 1 / 2 = S S P S P N L T
and Γ is the K 1 × L controllability matrix which can be expressed by the matrices P and D.
Γ = s n 1 / 2 V s n * = D P D P L 1 D
Considering that computational complexity of (19)–(21) is enlarged observably for large data sets, operations (22)–(24) are used for alternative steps which have lower computational load but with minimal calculation error.
H ˜ * H ˜ = V s n s n 2 V s n * = V s n Z V s n *
Ω H V s n Z 1 / 4
Γ = Z 1 / 4 V s n *
where Z = s n 2 .
As an analogy of the open-loop matrix in 1-D SSA, the augmented open-loop matrix P can be derived from the observability matrix by using least square.
P = Ω r f * Ω r f 1 Ω r f * Ω r l
Or it can be computed by the controllability matrix Γ .
P = Γ c l Γ c f * Γ c f Γ c f * 1
where Ω r f is the first (NL)M rows of Ω and Ω r l denotes the last (NL)M rows of Ω ; Γ c f represents the first (L − 1) columns of Γ and Γ c l is the last (L − 1) columns of Γ . From (20) and (21), these matrices can be rewritten as
S S P S P N L 1 T = Ω r f
S P S P 2 S P N L T = Ω r l
D P D P L 2 D = Γ c f
P D P 2 D P L 1 D = Γ c l
More details of the derivation of the matrices S, P, and D are listed in Appendix B.
According to (11), the vector [ p 1 p 2 p K 1 ] is obtained by performing the SVD.
Λ = Q * P Q
d i a g ( Λ ) = p 1 p 2 p K 1
where d i a g ( Λ ) denotes the diagonal element of matrix Λ .

3.2. Pole Estimation of the Aspect Dimension

In this step, multiple-range search, which applies one-dimensional (1-D) state–space approach (SSA) to the 1-D data for each down-range cell, is used for pole estimation of the aspect dimension and pole adaptive pairing. Details are as follows.
We introduce the Vandermonde sub-matrix O.
O = p 1 1 p 1 2 p 1 N p 2 1 p 2 2 p 2 N p K 1 1 p K 1 2 p K 1 N
From (5), the noiseless signatures Y ˜ can be factorized as a product of the pairing matrix G and O, where each column of G is associated with each row of O.
Y ˜ = t = 1 l 1 a 1 t s 1 t t = 1 l 2 a 2 t s 2 t t = 1 l K 1 a K 1 t s K 1 t t = 1 l 1 a 1 t s 1 t 2 t = 1 l 2 a 2 t s 2 t 2 t = 1 l K 1 a K 1 t s K 1 t 2 t = 1 l 1 a 1 t s 1 t M t = 1 l 2 a 2 t s 2 t M t = 1 l K 1 a K 1 t s K 1 t M O = G O
Thus, the pairing matrix G in (34) can be calculated by using least square, i.e.
G = Y ˜ O * O O * 1
The k-th column of the pairing matrix G is presented as
G ( k ) = t = 1 l k a k t s k t 1 t = 1 l k a k t s k t 2 t = 1 l k a k t s k t M T
where k = 1, 2, …, K1.
Using (36), each column of the pairing matrix is constructed in a same structure to 1-D DE model defined in (1). This indicates that the pole s k t and a k t , which correspond to the k-th pole p k in (32), can be solved by using 1-D SSA [22] to each column of the pairing matrix G . For example, for the first column of the pairing matrix G corresponds to p 1 , the parameter matrices (C, A, B) can be obtained by using 1-D SSA in (4).
G ( 1 ) = C B C A B C A M 1 B T
The eigenvalue decomposition of the open-loop matrix A leads to.
y ˜ ( 1 ) y ˜ ( 2 ) y ˜ ( N ) = C B C A B C A N 1 B
The matching pole s 1 t and the corresponding amplitude a 1 t are computed as.
s 11 , s 12 , , s 1 l 1 = Λ 1 1 , 1 Λ 1 2 , 2 Λ 1 l 1 , l 1
a 11 , a 12 , , a 1 l 1 = C M 1 * · M 1 1 B
Thus, the pairing matrix G can be searched column by column until all poles s k t and the corresponding amplitudes a k t are estimated. No extra pairing scheme is required because the poles [ s k 1 , s k 2 , , s k l k ] and the amplitudes [ a k 1 , a k 2 , , a k l k ] have already been adaptively paired to the poles p k in (34). In addition, considering that the model order may be overestimated sometime, the searched pole pairs are usually checked again according to their amplitudes. Finally, locations and damping factors of all the scattering points can be obtained from (41) and (42).
( r 1 k , r 2 k ) = ( Arg { s k } 4 π f c Δ θ / c , Arg { p k } 4 π Δ f / c )
( β 1 k , β 2 k ) = ( ln ( s k ) f c Δ θ , l n ( p k ) Δ f )

4. Results and Discussion

In this Section, three examples are presented to demonstrate the usefulness of the proposed procedure, i.e., the numerical signatures with 14 point scattering centers, the numerical signatures of a sphere tipped cone-cylinder-frustum combination model, the measured ISAR data for an aircraft model. The results obtained by 2-D ESPRIT are used for comparison since it is one of the very few techniques which have been used in real radar applications [17,32].

4.1. Numerical Signatures with Point Scattering Centers

In this example, the noisy signatures composed of 14-point scattering centers are considered to be as follows:
y ( m , n ) = k = 1 14 a k exp j 2 π r 1 k / c + β 1 k f c θ n exp j 2 π r 2 k / c + β 2 k f m + w m , n
where (M, N) = (41, 41); a k is set to be 1; f c = 12 GHz and Δ f = 150 MHz ; θ 0 = 0 rad and Δ θ = 0.0125 rad; The setting of coupled ranges are shown in Figure 1; All the damping factors including β 1 k and β 2 k are set as −0.05/ ( f c Δ θ ) and −0.05/ Δ f except the far left scattering points ( β 11 , β 21 ) = (−0.02/ ( f c Δ θ ) , −0.02/ Δ f ); w ( m , n ) denotes the additive white Gaussian noise.
As shown in Figure 1, these scattering points form a missile-like shape in Cartesian coordinates. It contains 13 pole pairs which have repeated poles in either the down-range or aspect dimensions. The noisy signatures in space domain by Fourier transform-based imaging algorithm (add Taylor window) are displayed in Figure 2a,b. As can be seen, the positions and decay rates of these scattering points are consistent with the corresponding ranges and damping factors. To sparse representation of these noisy signatures, the parameter for 2-D ASSA is chosen to be L = 20; the parameters for 2-D ESPRIT are set as: (P, Q) = (20, 20), β = 0.8 , which are suggested in [17]. For each signal to noise ratio (SNR), the number of Monte Carlo simulation is 200. A series of range estimation results in different SNRs are presented in Figure 3. Please note that the model order in 2-D ESPRIT is pre-specified because the singular values-based criterion [28,29] cannot be used when there are repeated poles in either the down-range or aspect dimensions. Because of different pairing strategies, K1 and K in 2-D ASSA could be estimated by the eigenvalue sequences transform algorithm (Appendix A) when SNR > 10 dB. However, the numbers K1 and K in 2-D ASSA should be pre-specified or use the other criterion [30] when the noise level is higher (SNR 10 dB).
As can be seen in Figure 3a,b, the positions of scattering points estimated by 2-D ESPRIT are generally according to the preset poles in Figure 1, except for some missing or incorrect points (shown by the black rectangular). The possible reason of this problem is that the pairing procedure in 2-D ESPRIT may provide incorrect pole pairs when there are repeated poles in either the down-range or aspect dimension, i.e., six pole pairs ( s 1 , p 1 ) , ( s 1 , p 2 ) , ( s 1 , p 3 ) , ( s 2 , p 1 ) , ( s 2 , p 2 ) , ( s 2 , p 3 ) . In contrast, the results estimated by 2-D ASSA showed higher accuracy than 2-D ESPRIT for different pairing strategies.
The estimation accuracy of cross/down ranges and damping factors are displayed in Figure 4. Root mean square error (RMSE) is used as the evaluating indicator which is defined in (44). As we can see, estimation accuracy of these two algorithms are basically at the same level although the size of Hankel matrix used by 2-D ASSA is smaller than the block-Hankel matrix defined by 2-D ESPRIT.
δ R M S E = 10 log 10 1 M c t = 1 M c X est X real 2
where M c = 2000 denotes the number of Monte Carlo runs, X est and X real are the esimated and real parameters.
Figure 5 also presents the statistic result with 2000 Monte Carlo runs (SNR = 0 dB). The result demonstrates that the multiple-range search strategy used by 2-D ASSA is robust in pole-pairing for low SNR.

4.2. Numerical Signatures of Computer-Aided Design (CAD) Model

The numerical signatures are obtained using method of moment (MOM), where the target is from a computer-aided design (CAD) model of a sphere tipped cone-cylinder-frustum combination, (shown in Figure 6). The data was calculated from 8–12 GHz in 10 MHz frequency step size, view angle ranging from −5 to 5 deg with an increment of 0.25 deg. Figure 7a presents 2-D radar image processed using 2-D FFT. As it can be seen, strong scattering points can be observed at differential discontinuities, such as the base-edge, the body groove, and nosetip. To sparsely represent the numerical data, the parameters of this example are set as follows. For 2-D ASSA, L is set as N2. For 2-D ESPRIT, (P, Q) = (M2, N2), where M2 = ceil(M/2) denotes the smallest integer less than or equal to M/2, and β = 0.8.
Location estimation of key scattering points for these two algorithms are shown in Figure 7b–e. As can be seen, when the number of scattering centers is set to be 14, all key scattering points are accurately estimated by these two algorithms except for some minor differences. However, when K is set to be 18, 2-D ESPRIT encountered the pole-pairing problem and could not provide the right estimation. Relative reconstruction error (RRE) δ R R E (defined in (45)) of these two algorithms are shown in Figure 8. We can see that the RRE of 2-D ASSA has been falling when K increased from 10 to 30, whereas the RRE of 2-D ESPEIT stopped falling after K = 14. The result shows that 2-D ASSA tends to be more robust in pole-pairing for different numbers of scattering centers.
δ R R E = 1 M N m = 1 M n = 1 N 20 log 10 Y r e c o n ( m , n ) / Y r e a l ( m , n )
where Yrecon denotes the signatures reconstructed by the estimated parameters. Yreal is the original signatures.
Figure 9a displays the sub-band data image (8–10 GHz, from −5 to 5 deg) which was extracted from the full-band data. Compared to the full-band data image in Figure 7a, scattering centers of the base-edge and the body groove cannot be distinguished when the signal bandwidth is only 2 GHz. Figure 9b presents the scattering points estimated by 2-D ASSA and Figure 9c shows the full-band data image extrapolated by these extracted pole pairs. As can be seen, scattering centers of the base-edge and the body groove are distinguished from the extracted scattering points. The full-band data image extrapolated by 2-D ASSA is in keeping with the main scattering centers distributed in the original full-band data image (Figure 7a). It is clear from these figures that 2-D ASSA results in higher-resolution images than traditional 2-D FFT imagery.

4.3. Measured ISAR Signatures

In the third example, the measured ISAR signatures are acquired in an indoor test range where a mocked aircraft model is used. The data was collected using S-band radar with center frequency 3 GHz and 1.5 GHz bandwidth. The range of aspect angle is from −4° to 4° with an interval of 0.1°. Figure 10a shows the ISAR image of the aircraft model generated by 2-D FFT. As it can be seen, the scattering distribution of this model is more complex than the previous model, i.e., many scattering points are densely distributed in the fuselage. The estimation parameter for 2-D ASSA is set to be L = N 2 ; For 2-D ESPRIT, (P, Q) = M 2 , N 2 and β = 0.8.
Location estimation of main scattering centers for measured data are shown in Figure 10b–e. As we can see, both 2-D ESPRIT and 2-D ASSA could estimate the right positions of strong scattering centers. The scattering points estimated by 2-D ASSA are more densely distributed around the strong scattering centers than 2-D ESPRIT’s for different pole-pairing strategies. A few relatively weak scattering points, i.e., the red point in rectangle box, are wrongly estimated by 2-D ESPRIT when the number of scattering points is set to be 115, and more wrongly estimated pole pairs can be seen in K > 115 which are not displayed in Figure 10.
Figure 11a–h display a set of reconstructed 2-D images by the estimated pole pairs. Please note that the number of scattering points K in 2-D ESPRIT should be equal to the model order whereas K in 2-D ASSA could be larger than the model order K1 (the number of non-repeated poles in the down-range dimension) when there are one-to-multiple pairing poles. That is the reason K could be set to 607 which is already exceed the limitation of model order in [17]. As we can see from the figures, all those two algorithms can reconstruct those strong scattering centers of target which are in consistent with the estimated poles in Figure 10a–d. As displayed in Figure 11f–h, more and more relatively weak scattering centers have been reconstructed by 2-D ASSA with the increasing number of K and the reconstruction result with compression ratio of 91.04% is extremely similar to the original data in Figure 10a.
Evaluation of the reconstructed results by 2-D ASSA are listed in Table 2. Compression ratio (CR) ε C R and image similarity degree (ISD) are defined in (46) and (47). The RRE δ R R E is used for evaluating the reconstructed result in frequency-domain, whereas the ISD γ I S D is for the result in image-domain. From the table, we can see that the numbers of scattering points in 2-D ASSA are negatively correlated with RRE (or ISD). In contrast, 2-D ESPRIT performs higher RRE than 2-D ASSA when the number of scattering centers remains the same. Please note that part of evaluation results by 2-D ESPRIT are abnormal because 2-D ESPRIT cannot do the right reconstruction when K > 307. The possible reason is that for 2-D ESPRIT, all amplitudes factors a 1 a 2 a K (defined in (3)) are estimated together by using least square technique after all the pole pairs are confirmed. Thus, numerous incorrect pairing poles estimated by 2-D ESPRIT will lead to entirely wrong reconstruction. In comparison, 2-D ASSA performs better robustness for the measurement data.
ε C R = M × N n P × K M × N × 100 %
γ I S D = I r e c o n I r e a l I r e c o n 2 I r e a l 2
where n P = 3 represents parameter number of each paired pole a k , s k , p k , K is the number of scattering points, M × N denotes size of the original signatures. Irecon, Ireal represent the reconstructed and original 2-D images-domain data.
Running time of the processed examples are listed in Table 3. All the results are carried out by Matlab (R2016b) with same hardware platform: Intel Core i7-6900k 3.2 GHz and 128 G memory. Please note that the running time of 2-D ESPRIT does not contain the estimation of model order. According to the table, 2-D ASSA performs higher efficiency than 2-D ESPRIT in processing the same data with the following differences.
(1) The size of the augmented Hankel matrix used by 2-D ASSA is approximate 1/(M/4) of the block-Hankel matrix used by 2-D ESPRIT when L = Q = N2 and P =M2. Moreover, the data size of the Hankel matrix could be further reduced by using the operation (22)–(24);
(2) No extra pairing scheme is required by using multiple-range search strategy which can finish the pairing process once the pole estimation of the aspect dimension is confirmed. Moreover, calculation of the block-Hankel matrix for the aspect dimension in 2-D ESPRIT is simplified as calculating multiple small Hankel matrices with size M2 × (MM2 + 1).

5. Conclusions

In this paper, we have presented a two-dimensional augmented state–space approach for sparse representation of 2-D wideband radar signatures collected on man-made targets. To do this, a two-step procedure, i.e., an augmented state–space approach followed by multiple-range search strategy, is proposed to estimate the complex amplitudes and poles in down-range and aspect dimensions. In general, there are mainly two contributions provided in this paper. First, we develop a computationally efficient approach by adopting several time-saving operations, whereas the pole-estimation accuracy is still at the same level; second, the proposed approach can apparently alleviate the pole-pairing problem by using the multiple-range search strategy.
Numerical as well as measured ISAR data are processed to validate the proposed approach. Experimental results demonstrate that 2-D ASSA is robust and accurate in pole-paring for different SNRs, and is applicable for sparse representation of 2-D wideband radar signatures collected on man-made targets with low computational cost.
Future works are considered to be as follows. On the one hand, physical meanings of the extracted parameters of scattering centers by 2-D ASSA might be further studied to the possible applications in automatic target recognition (ATR). On the other hand, the extension of the proposed approach to 3-D radar signatures obtained from different azimuth and elevation angles would also be an interesting study.

Author Contributions

Methodology, K.W.; validation, K.W. and X.X.; writing—original draft preparation, K.W.; writing—review and editing, K.W. and X.X.

Funding

This research was funded by National Natural Science Foundation of China: 61371005.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In the appendix, we provide an eigenvalue sequences transform algorithm to estimate the model order K1. Details are as follows.
Suppose that the size of the input matrix H is m × n ( m n ), do the SVD of the matrix H * H .
H * H = V 2 V *
where V is unitary matrix; 2 is the diagonal matrix.
Λ = λ 1 , λ 2 , , λ n is the diagonal vector of . It can be normalized as:
w i = λ i min ( Λ ) max ( Λ ) min ( Λ )
where i = 1, 2, …, n; max ( Λ ) and min ( Λ ) denote the maximum and minimum elements of Λ .
The coordinate form of the normalized eigenvalue vector is:
Q 1 , Q 2 , , Q n = 1 , w 1 , 2 , w 2 , , n , w n
where Q n denotes the n-th points.
Based on information theory and principal factor analysis [28], the separable signals and noises in H often have the following hypothesis.
w 1 > > w K 1 > w K 1 + 1 > w n signal   space | noise space
w K 1 w K 1 + 1 > w K 1 + 1 w K 1 + 2 w n 1 w n
where K1 is the number of separable signals which also denotes the model order.
From (A5), it can be deduced that the point Q K 1 belonging to Q 1 , Q 2 , , Q n has the longest distance to line Q 1 D which has the following condition.
k Q i Q K 1 + 1 > k Q 1 D > k Q K 1 + 1 Q n
where i = 1, 2, …, K1, k Q i Q K 1 + 1 denotes the slope of line Q i Q K 1 + 1 , D is a point which satisfies the slope condition (A6).
Here, an approximate calculation is used to estimate the slop of Q 1 D .
k Q 1 D k Q p 1 Q n
g p = max g 1 , g 2 , , . g n
g i = k Q 1 Q n i 1 w i w 1 k Q 1 Q n 2 + 1
where Q p represents the point in Q 1 , Q 2 , , Q n which has the longest distance to line Q 1 Q n ; g i denotes the distance from point Q i to Q 1 Q n ; i = 1, 2, …, n.
Thus, Q K 1 can be confirmed by searching the longest distance from point Q i to line Q 1 D .
g i = k Q 1 D i 1 w i w 1 k Q 1 D 2 + 1
g K 1 + 1 = max g ( 1 ) , g ( 2 ) , , . g ( n )
where i = 1, 2, …, n; g i denotes the distance from point Q i to Q 1 D .

Appendix B

From the observability matrix Ω given by (20), controllability matrix Γ given by (21) and the matrices Ω r f , Ω r l , Γ c f , Γ c l in (27),(28),(29) and (30), it is not difficult to deduce that the augmented open-loop matrix P satisfies the following matrix equations.
Ω r f P = Ω r l
P Γ c f = Γ c l
Then the augmented open-loop matrix P can be obtained by least squares.
P = Ω r f * Ω r f 1 Ω r f * Ω r l
or
P = Γ c l Γ c f * Γ c f Γ c f * 1
Moreover, here are two ways to compute the corresponding constant matrices S and P. The first way is for matrix P computed by (A14), the corresponding matrix S is defined as the first M rows of the observability matrix Ω .
S = Ω 1 , 1 Ω 1 , K 1 Ω M , 1 Ω M , K 1
Using (A14) and (A16), we calculate the matrix Ω ˜
Ω ˜ = S S P S P N 1
From (15) and (A17), the original 2-D signatures Y can be written as
Y = Ω ˜ D
By using least square again, the matrix D is
D = Ω ˜ * Ω ˜ 1 Ω ˜ * Y
The second way is for matrix P computed by (A15), the corresponding matrix D is defined as
D = Γ 1 , 1 Γ K 1 , 1
Using (A15) and (A20), we calculate the matrix Γ ˜
Γ ˜ = D P D P N 1 D
Similar to (A18), the original 2-D signatures Y can also be rewritten as
Y = S Γ ˜
By using least square, the matrix S is
S = Y Γ ˜ * Γ ˜ Γ ˜ * 1
In this appendix, we provide two ways to estimate the parameter matrices (S, P, D). For consideration of the robustness to noise, (A23), (A14) (or (A15)) and (A19) are preferred. It is also worth mentioning that 2-D ASSA has two ways, which include the matrix form (S, P, D) and the pole form a k , s k , p k , to reconstruct the original signatures Y. The matrix form could provide more accurate reconstruction result than the pole form but with much more parameters.

References

  1. Wang, L.; Zhao, L.F.; Bi, G.; Wan, C. Sparse Representation-Based ISAR Imaging Using Markov Random Fields. IEEE J. Sel. Topics Appl. Earth Observ. Remote Sens. 2015, 8, 3941–3953. [Google Scholar] [CrossRef]
  2. Lasserre, M.; Bidon, S.; Chevalier, F.L. New Sparse-Promoting Prior for the Estimation of a Radar Scene with Weak and Strong Targets. IEEE Trans. Signal Process. 2016, 64, 4634–4643. [Google Scholar] [CrossRef] [Green Version]
  3. Tseng, N. A very Efficient RCS Data Compression and Reconstruction Technique. Master’s Thesis, The Ohio State University, Columbus, OH, USA, 1992. [Google Scholar]
  4. Dong, G.G.; Wang, N.; Kuang, G.Y. Sparse Representation of Monogenic Signal: With Application to Target Recognition in SAR Images. IEEE Signal Process. Lett. 2014, 21, 952–956. [Google Scholar]
  5. Xu, X.J.; He, F.Y. An Iterative Procedure for Ultra-Wideband Imagery of Space Objects from Distributed Multi-Band Radar Signatures. Proc. SPIE. 2014, 9227. [Google Scholar] [CrossRef]
  6. Li, J.; Stoica, P. An Adaptive Filtering Approach to Spectral Estimation and SAR Imaging. IEEE Trans. Signal Process. 1996, 44, 1469–1484. [Google Scholar]
  7. Graaf, S.R.D. Parametric Estimation of Complex 2-D Sinusoids. In Proceedings of the IEEE Fouth Annual ASSP Workshop on Spectrum Estimation and Modeling, Minneapolis, MN, USA, 3–5 August 1988; IEEE: Piscataway, NJ, USA; pp. 391–396. [Google Scholar]
  8. Choi, I.S.; Kim, H.T. Two-Dimensional Evolutionary Programming-Based CLEAN. IEEE Trans. Aerosp. Electron. Syst. 2003, 39, 373–382. [Google Scholar] [CrossRef]
  9. Varshney, K.R.; Cethin, M.; Fisher, J.W.; Willsky, A.S. Sparse Representation in Structured Dictionaries With Application to Synthetic Aperture Radar. IEEE Trans. Signal Process. 2008, 56, 3548–3561. [Google Scholar] [CrossRef]
  10. Potter, L.C.; Erthin, E.; Parker, J.T.; Cethin, M. Sparsity and Compressed Sensing in Radar Imaging. Proc. IEEE. 2010, 98, 1006–1020. [Google Scholar] [CrossRef]
  11. Liu, H.C.; Jiu, B.; Liu, H.W.; Bao, Z. Superresolution ISAR Imaging Based on Sparse Bayesian Learning. IEEE Trans. Geosci. Remote Sens. 2014, 52, 5005–5013. [Google Scholar]
  12. Odendaal, J.W.; Barnard, E.; Pistorius, C.W.I. Two-Dimensional Superresolution Radar Imaging Using the MUSIC Algorithm. IEEE Trans. Antennas Propag. 1994, 42, 1386–1391. [Google Scholar] [CrossRef]
  13. Hua, Y.; Sarkar, T.K. Matrix Pencil Method for Estimating Parameters of Exponentially Damped/undamped Sinusoids in Noise. IEEE Trans. Signal Process. 1990, 36, 814–824. [Google Scholar] [CrossRef]
  14. Hua, Y. Estimating Two-Dimensional Frequencies by Matrix Enhancement and Matrix Pencil. IEEE Trans. Signal Process. 1992, 40, 2267–2280. [Google Scholar] [CrossRef]
  15. Sacchini, J.J.; Steedly, W.M.; Moses, R.L. Two-Dimensional Prony Modeling and Parameter Estimation. IEEE Trans. Signal Process. 1993, 41, 3127–3137. [Google Scholar] [CrossRef]
  16. Vanpoucke, F.; Moonen, M.; Berthoumieu, Y. An Efficient Subspace Algorithm for 2-D Harmonic Retrieval. In Proceedings of the ICASSP ’94. IEEE International Conference on Acoustics, Speech and Signal Processing, Adelaide, SA, Australia, 19–22 April 1994; pp. 461–464. [Google Scholar]
  17. Rouquette, S.; Najim, M. Estimation of Frequencies and Damping Factors by Two-Dimensional ESPRIT Type Methods. IEEE Trans. Signal Process. 2001, 49, 237–245. [Google Scholar] [CrossRef]
  18. Piou, J.E. System Realization Using 2-D Output Measurements. In Proceedings of the 2004 American Control Conference (ACC), Boston, MA, USA, 30 June–2 July 2004. [Google Scholar]
  19. Piou, J.E.; Dumanian, A.J. Application of the Fornasini-Marchesini First Model to Data Collected on a Complex Target Model. In Proceedings of the 2014 American Control Conference (ACC), Portland, OR, USA, 4–6 June 2014. [Google Scholar]
  20. Chen, F.J.; Fung, C.C.; Kok, C.W.; Kwong, S. Estimation of Two-Dimensional Frequencies Using Modified Matrix Pencil Method. IEEE Trans. Signal Process. 2007, 55, 718–724. [Google Scholar] [CrossRef]
  21. Qian, C.; Huang, L.; So, H.C.; Sidiropoulos, N.D.; Xie, J. Unitary PUMA Algorithm for Estimating the Frequency of a Complex Sinusoid. IEEE Trans. Signal Process. 2015, 63, 5358–5368. [Google Scholar] [CrossRef]
  22. Naishadham, K.; Piou, J.E. A Robust State-Space Model for the Characterization of Extended Returns in Radar Target Signatures. IEEE Trans. Antennas Propag. 2008, 56, 1742–1750. [Google Scholar] [CrossRef]
  23. Ying, C.J.; Chiang, H.C.; Moses, R.L.; Potter, L.C. Complex SAR Phase History Modeling Using Two-dimensional Parametric Estimation Techniques. Proc. SPIE. 1996, 2757, 174–185. [Google Scholar] [CrossRef]
  24. Xu, X.J.; He, F.Y. Fast and Accurate RCS Calculation for Squat Cylinder Calibrators. IEEE Antennas Propagation Mag. 2015, 57, 33–40. [Google Scholar] [CrossRef]
  25. Augusto, A.; Vincenzo, C.; Antonio, D.M.; Luca, P. High Range Resolution Profile Estimation via a Cognitive Stepped Frequency Technique. IEEE Trans. Aerosp. Electron. Syst. 2019, 1, 444–458. [Google Scholar]
  26. Pia, A.; Augusto, A.; Antonio, D.M.; Luca, P.; Silvia, L.U. HRR profile estimation using SLIM. IET Radar Sonar Navig. 2019, 4, 512–521. [Google Scholar]
  27. Potter, L.C.; Moses, R.L. Attributed Scattering Centers for SAR ATR. IEEE Trans. Image Process. 1997, 6, 79–91. [Google Scholar] [CrossRef] [PubMed]
  28. Wax, M.; Kailath, T. Detection of Signals by Information Theoretic Criteria. IEEE Trans. Signal Process. 1995, 33, 387–392. [Google Scholar] [CrossRef]
  29. Wax, M.; Ziskind, I. Detection of the Number of Coherent Signals by the MDL Principle. IEEE Trans. Signal Process. 1989, 27, 1190–1196. [Google Scholar] [CrossRef]
  30. Aksasse, B.; Radouane, L. Two-Dimensional Autoregressive (2-D AR) Model Order Estimation. IEEE Trans. Signal Process. 1999, 47, 2072–2077. [Google Scholar] [CrossRef]
  31. Kung, S.Y.; Arun, K.S.; Rao, D.V.B. State-space and Singular-Value Decomposition based Approximation Methods for the Harmonic Retrieval Problem. J. Opt. Soc. Am. 1983, 73, 1799–1811. [Google Scholar] [CrossRef]
  32. Michael, L.B. Two-Dimensional ESPRIT With Tracking for Radar Imaging and Feature Extraction. IEEE Trans. Antennas Propag. 2004, 2, 524–532. [Google Scholar]
Figure 1. Spatial distribution of the simulated poles.
Figure 1. Spatial distribution of the simulated poles.
Sensors 19 04631 g001
Figure 2. Simulated radar image with different SNR by Fourier transform. (a) Simulated radar image with SNR = 5 dB; (b) Simulated radar image with SNR = 0 dB.
Figure 2. Simulated radar image with different SNR by Fourier transform. (a) Simulated radar image with SNR = 5 dB; (b) Simulated radar image with SNR = 0 dB.
Sensors 19 04631 g002
Figure 3. Estimation of the coupled ranges by 2-D ASSA and 2-D ESPRIT with different SNR. (a) Estimation by 2-D ESPRIT with SNR = 5 dB; (b) Estimation by 2-D ESPRIT with SNR = 0 dB; (c) Estimation by 2-D ASSA with SNR = 5 dB; (d) Estimation by 2-D ASSA with SNR = 0 dB.
Figure 3. Estimation of the coupled ranges by 2-D ASSA and 2-D ESPRIT with different SNR. (a) Estimation by 2-D ESPRIT with SNR = 5 dB; (b) Estimation by 2-D ESPRIT with SNR = 0 dB; (c) Estimation by 2-D ASSA with SNR = 5 dB; (d) Estimation by 2-D ASSA with SNR = 0 dB.
Sensors 19 04631 g003
Figure 4. Estimation accuracy of r11, r21, β11 and β21 with different SNR. (a) Estimation of r11 with different SNR; (b) Estimation of r21 with different SNR; (c) Estimation of β11 with different SNR; (d) Estimation of β21 with different SNR.
Figure 4. Estimation accuracy of r11, r21, β11 and β21 with different SNR. (a) Estimation of r11 with different SNR; (b) Estimation of r21 with different SNR; (c) Estimation of β11 with different SNR; (d) Estimation of β21 with different SNR.
Sensors 19 04631 g004
Figure 5. Estimation of the coupled ranges by 2-D ASSA with 2000 Monte Carlo runs (SNR = 0 dB).
Figure 5. Estimation of the coupled ranges by 2-D ASSA with 2000 Monte Carlo runs (SNR = 0 dB).
Sensors 19 04631 g005
Figure 6. CAD model of the sphere tipped cone-cylinder-frustum combination.
Figure 6. CAD model of the sphere tipped cone-cylinder-frustum combination.
Sensors 19 04631 g006
Figure 7. 2-D radar image of the CAD model and location estimation of main scattering points result by 2-D ESPRIT and 2-D ASSA. (a) 2-D radar image in aspect angle from −5 to 5 deg; (b) Pole estimation by 2-D ESPRIT(K = 14); (c) Pole Estimation by 2-D ASSA(K = 14); (d) Pole estimation by 2-D ESPRIT(K = 18); (e) Pole Estimation by 2-D ASSA(K = 18).
Figure 7. 2-D radar image of the CAD model and location estimation of main scattering points result by 2-D ESPRIT and 2-D ASSA. (a) 2-D radar image in aspect angle from −5 to 5 deg; (b) Pole estimation by 2-D ESPRIT(K = 14); (c) Pole Estimation by 2-D ASSA(K = 14); (d) Pole estimation by 2-D ESPRIT(K = 18); (e) Pole Estimation by 2-D ASSA(K = 18).
Sensors 19 04631 g007
Figure 8. Relative reconstruction error of the numerical signatures of CAD model with different numbers of scattering centers.
Figure 8. Relative reconstruction error of the numerical signatures of CAD model with different numbers of scattering centers.
Sensors 19 04631 g008
Figure 9. Sub-band data image and extrapolated data image by 2-D ASSA. (a) Sub-band data image by 2-D FFT (8–10 GHz); (b) Scattering points estimated by 2-D ASSA; (c) Full-band data image extrapolated by 2-D ASSA (8–12 GHz).
Figure 9. Sub-band data image and extrapolated data image by 2-D ASSA. (a) Sub-band data image by 2-D FFT (8–10 GHz); (b) Scattering points estimated by 2-D ASSA; (c) Full-band data image extrapolated by 2-D ASSA (8–12 GHz).
Sensors 19 04631 g009
Figure 10. Pole-estimation result of the measured data with different K. (a) Measured ISAR image of the aircraft model by using 2-D FFT (b) Pole estimation by 2-D ESPRIT (K = 57); (c) Pole estimation by 2-D ASSA (K = 57); (d) Pole estimation by 2-D ESPRIT(K = 115); (e) Pole estimation by 2-D ASSA (K = 115).
Figure 10. Pole-estimation result of the measured data with different K. (a) Measured ISAR image of the aircraft model by using 2-D FFT (b) Pole estimation by 2-D ESPRIT (K = 57); (c) Pole estimation by 2-D ASSA (K = 57); (d) Pole estimation by 2-D ESPRIT(K = 115); (e) Pole estimation by 2-D ASSA (K = 115).
Sensors 19 04631 g010
Figure 11. Reconstruction result of the measured data with different K. (a) Reconstruction result by 2-D ESPRIT(K = 57); (b) Reconstruction result by 2-D ESPRIT (K = 115); (c) Reconstruction result by 2-D ESPRIT (K = 304); (d) Reconstruction result by 2-D ASSA(K = 57);(e) Reconstruction result by 2-D ASSA(K = 115); (f) Reconstruction result by 2-D ASSA(K = 304); (g) Reconstruction result by 2-D ASSA(K = 507); (h) Reconstruction result by 2-D ASSA(K = 607).
Figure 11. Reconstruction result of the measured data with different K. (a) Reconstruction result by 2-D ESPRIT(K = 57); (b) Reconstruction result by 2-D ESPRIT (K = 115); (c) Reconstruction result by 2-D ESPRIT (K = 304); (d) Reconstruction result by 2-D ASSA(K = 57);(e) Reconstruction result by 2-D ASSA(K = 115); (f) Reconstruction result by 2-D ASSA(K = 304); (g) Reconstruction result by 2-D ASSA(K = 507); (h) Reconstruction result by 2-D ASSA(K = 607).
Sensors 19 04631 g011
Table 1. Synopsis of the 2-D subspace-based algorithms.
Table 1. Synopsis of the 2-D subspace-based algorithms.
TechniqueRationalePole (1) Pairing Scheme
2-D MUSICSignal-noise subspace decompositionNo need for pole-pairing
MEMPEnhanced matrix decompositionMaximizing a certain criterion
2-D TLS PronyTLS-based Prony modelMinimizing a certain distance
ACMPSignal subspace decompositionRank-restoration scheme
2-D ESPRITShift-invariance structure of the signal subspaceJoint diagonalization scheme
2-D system realizationState–space modelAlgebraic method
(1). Pole denotes a transfer function of the scattering center collected on target. Usually, the scattering centers are characterized via pairing the complex amplitudes and the poles in down-range and aspect dimensions.
Table 2. Evaluation of the reconstruction result of the measured data by 2-D ESPRIT and 2-D ASSA.
Table 2. Evaluation of the reconstruction result of the measured data by 2-D ESPRIT and 2-D ASSA.
2-D ASSA2-D ESPRIT
Number of Scattering δ R R E ( dB ) γ I S D δ R R E ( dB ) γ I S D
K = 57 ( ε C R = 99.16%)3.260.903.930.88
K = 115 ( ε C R = 98.30%)2.600.923.200.90
K = 156 ( ε C R = 97.70%)2.530.923.110.90
K = 207 ( ε C R = 96.95%)2.470.932.960.91
K = 256 ( ε C R = 96.22%)2.270.942.940.91
K = 304 ( ε C R = 95.51%)2.070.942.900.91
K = 407 ( ε C R = 93.99%)1.760.9541.280.26
K = 507 ( ε C R = 92.52%)1.600.9668.690.12
K = 607 ( ε C R = 91.04%)1.510.9643.570.24
Table 3. Running time of the numerical radar signatures and the measured ISAR data by 2-D ESPRIT and 2-D ASSA.
Table 3. Running time of the numerical radar signatures and the measured ISAR data by 2-D ESPRIT and 2-D ASSA.
MethodK Hankel   Matrix   ( p k ) Hankel   Matrix   ( s k ) Running Time
Numerical data (41 × 401)2-D ESPRIT144000 × 44444000 × 444434.81 s
2-D ASSA148282 × 200K1 × 20 × 220.94 s
Measured ISAR data (81 × 251)2-D ESPRIT575000 × 53345000 × 533453.27 s
1155000 × 53345000 × 533494.18 s
2-D ASSA5710287 × 125K1 × 40 × 421.04 s
11510287 × 125K1 × 40 × 421.19 s

Share and Cite

MDPI and ACS Style

Wu, K.; Xu, X. Two-Dimensional Augmented State–Space Approach with Applications to Sparse Representation of Radar Signatures. Sensors 2019, 19, 4631. https://doi.org/10.3390/s19214631

AMA Style

Wu K, Xu X. Two-Dimensional Augmented State–Space Approach with Applications to Sparse Representation of Radar Signatures. Sensors. 2019; 19(21):4631. https://doi.org/10.3390/s19214631

Chicago/Turabian Style

Wu, Kejiang, and Xiaojian Xu. 2019. "Two-Dimensional Augmented State–Space Approach with Applications to Sparse Representation of Radar Signatures" Sensors 19, no. 21: 4631. https://doi.org/10.3390/s19214631

APA Style

Wu, K., & Xu, X. (2019). Two-Dimensional Augmented State–Space Approach with Applications to Sparse Representation of Radar Signatures. Sensors, 19(21), 4631. https://doi.org/10.3390/s19214631

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop