High Range Resolution Profile Automatic Target Recognition Using Sparse Representation
High Range Resolution Profile Automatic Target Recognition Using Sparse Representation
Journal of
Aeronautics
Chinese Journal of Aeronautics 23(2010) 556-562 www.elsevier.com/locate/cja
Abstract
Sparse representation is a new signal analysis method which is receiving increasing attention in recent years. In this article, a
novel scheme solving high range resolution profile automatic target recognition for ground moving targets is proposed. The
sparse representation theory is applied to analyzing the components of high range resolution profiles and sparse coefficients are
used to describe their features. Numerous experiments with the target type number ranging from 2 to 6 have been implemented.
Results show that the proposed scheme not only provides higher recognition preciseness in real time, but also achieves more
robust performance as the target type number increases.
Keywords: automatic target recognition; high range resolution profile; sparse representation; feature extraction; dictionary gen-
eration
In this article, we explore the use of sparse repre- quency properties. Moreover, because the freedom of
sentation into HRRP feature analysis and expect its the selection of s, u and [ is unlimited, the number of
performance in HRRP recognition. gc is infinite. That is, the created dictionary can be as
redundant as it is needed. However, in practical
2. Sparse Representation application, it is sometimes hard to decide the suitable
number of gc for a dictionary. Furthermore, as all at-
2.1. Theory oms in the dictionary derive from the same base func-
tion, the adaptability in signal component representa-
Sparse representation is originally introduced by S. tion is relatively limited.
G. Mallat and Z. F. Zhang in 1993[17], and developed Another category for dictionary generation is a
by D. L. Donoho[18-20] later. It is a method to analyze training sample-based method[14]. The dictionary is
signals under a redundant dictionary. combined with training samples from several target
types. These training samples are collected under a
ª x1 º wide scope of real scenes or their simulations where
«x » N the measured signals, or also called test samples re-
y AX [ A1c A2c " ANc ] « 2 » ¦ Aicxi (1) lated to training samples, probably take place. Due to
« # » i 1
« » the similar collecting conditions, training samples have
¬ xN ¼ similar characteristics with their associated test sam-
ples. The dictionary generated by these training sam-
where the measured signal y=[y1 y2 … yd]TęRd ples is able to represent the test samples more effi-
is a d-dimensional vector, the redundant dictionary ciently and precisely. Additionally, since the real
A [ A1c A2c " ANc ] ęRdhN is a matrix with size of scenes where signals take place are known or predict-
dhN, d˘N, the atom Aic ęA is a d-dimensional vec- able, the atom number in the dictionary can be decided
tor, i=1,2,…,N, the coefficient vector X=[x1 x2 … easily.
xN]TęRN is a N-dimensional vector obtained by y ex-
2.3. Approaches
panding over A. As d˘N, equations y=AX are under
determined. There are infinite X satisfying y=AX. The Finding the sparse representation by the minimum l0
freedom of selection of X facilitates the adaption to norm in Eq.(2) is complex and extremely hard[18].
various practical applications. Several algorithms have been designed to find subop-
In all the possible solutions, X with fewest nonzero timal or approximate solutions instead. There are three
coefficients, X0, is the sparse representation[21]: approaches widely used in sparse representation over
X0 min X 0
subject to y AX (2) redundant dictionaries. They are matching pursuit
X
(MP), basis pursuit (BP) and least absolute shrinkage
In fact, X0 is greatly attractive in actual demands and selection operation (LASSO) [17-23].
since it always means fewer sensors in a radar system MP[17-21], proposed by S. G. Mallat, is a greedy pur-
or less prerequisites for a machine. suit algorithm. It finds the sparse solution by an itera-
tive matching process. Compared with BP and LASSO,
2.2. Dictionary MP has the advantage of easy implementing and short
running time. However, since the solutions produced
Dictionary design is crucial for signal component by MP are just local optimal, they are sometimes not
analysis in sparse representation. Generally, atoms in a precise enough.
dictionary should meet two basic requirements: con- D. L. Donoho[19] simplified problems minimum l0
taining a wide range of types of atoms to better adapt norm into problems minimum l1 norm:
to different local structures of signals, and having suf- y AX
X1 min X subject to (4)
ficient number of atoms to make the dictionary redun- X 1
3.4. Recognition
Fig.1 Sparse coefficients obtained by sparse analysis. Fig.2 Performance comparison of SRBR and TBR in 2 to 6
recognition experiments.
and 270q aspect angles. It is seen that there are few
Experiments are implemented using MATLAB on a
sparse coefficients appearing in Fig.1. Significant
PC with 2.8 GHz CPU. The average running time for
sparse coefficients concentrate on the atoms related to
SRBR and TBR is 4.7 s and 2.7 s, respectively.
the type of test samples while sparse coefficients of
atoms corresponding to other target types are relatively 4.3. Discussion
much smaller, or even equal to zero.
From Fig.2, it is known that both of TBR and SRBR
4.2. Recognition comparison provide satisfied recognition results. However, SRBR
is superior to TBR throughout 2 to 6 recognition ex-
Two approaches are considered in recognition com- periments. The improvement in SRBR primarily bene-
parison. One is the currently widely used template fits from the insight of sparse representation into signal
based recognition (TBR) using HRRP magnitudes[3-5]. decomposition. For example, Fig.3 contains a template
The other is the sparse representation based recogni- matching result that a test sample of 2nd target type at
tion (SRBR) employing sparse coefficients proposed 14q aspect angle. The maximum correlation coefficient
in this article. Both approaches are implemented in 5 appears at the range between 659th template and 795th
sorts of recognition experiments. The target type template which corresponds to 5th target type. Hence,
number in these 5 sorts of recognition experiments the test sample is judged as the 5th target type accord-
No.5 Zhou Nuo et al. / Chinese Journal of Aeronautics 23(2010) 556-562 · 561 ·
signature classification using syntactic pattern recogni- [20] Donoho D L. For most large underdetermined systems
tion. IEEE Aerospace Conference. 2008: 1-9. of linear equations the minimal 1-norm solution is also
[6] Liu J, Zhang J Y, Zhao F. A new time-shift invariant the sparsest solution. Communications on Pure and
feature of radar HRRPs. Journal of electronics & In- Applied Mathematics 2006; 59(6): 797-829.
formation Technology 2008; 30(8): 1949-1953. [in [21] Yang J Y, Peng Y G, Xu W L, et al. Ways to sparse
Chinese] representation: an overview. Science in China Series F:
[7] Wang D W, Li Y J, Zhang K. Time-frequency signature Information Sciences 2009; 52(4): 695-703.
extraction and recognition of target based on [22] Amaldi E, Kann V. On the approximability of mini-
S-transform. Acta Aeronautica et Astronautica Sinica mizing nonzero variables or unsatisfied relations in
2009; 30(2): 305-309. [in Chinese] linear systems. Theoretical Computer Science 1998;
[8] Morris H C, de Pass M M. Wavelet feature extraction 209(1-2): 237-260.
of HRR radar profiles using generalized Gaussian dis- [23] Osborne M R, Presnell B, Turlach B A. A new ap-
tributions for automatic target recognition. Proceedings proach to variable selection in least squares problems.
of SPIE Signal Processing, Sensor Fusion, and Target IMA Journal of Numerical Analysis 2000; 20(3):
Recognition XIV. 2005; 5809: 165-175. 389-403.
[9] Liu H W, Yang Z, He K, et al. Radar high range resolu- [24] Ross T D, Worrell S W, Velten V J, et al. Standard SAR
tion profiles recognition based on wavelet packet and ATR evaluation experiments using the MSTAR public
subband fusion. IEEE International Conference on release data set. Proceedings of SPIE Algorithms for
Acoustics, Speech, and Signal Processing. 2005; 5: Synthetic Aperture Radar Imagery V. 1998; 3370:
445-448. 566-573.
[10] Chen B, Liu H W, Bao Z. An efficient kernel optimiza- [25] Gross D, Oppenheimer M, Kaler B, et al. Preliminary
tion method for high range resolution profile recogni- comparison of high range resolution signatures of
tion. Radar International Conference. 2006: 1-4. moving and stationary ground vehicles. Proceedings of
[11] Liu H W, Su H T, Bao Z. Radar high range resolution SPIE Algorithms for Synthetic Aperture Radar Imagery
profiles feature extraction based on kernel PCA and IX. 2002; 4727: 205-212.
kernel ICA. Advances in Neural Networks 2005; 3496: [26] Skolnik M l. Introduction to radar systems. 3rd ed.
913-918. Beijing ˖ Publishing House of Electronics Industry,
[12] Wang J Y, Yin Z K, Zhang C M. Sparse analysis and 2007.
preliminary application in signal and image. Sichuan: [27] Huang P K, Yin H C, Xu X J. Characteristics of radar
Southwest Jiaotong University Press, 2006. [in Chi- targets. Beijing˖Publishing House of Electronics In-
nese] dustry, 2005. [in Chinese]
[13] Wright J, Ganesh A, Zhou Z H, et al. Demo: robust [28] Shirman Y D. Computer simulation of aerial target
face recognition via sparse representation. IEEE Inter- radar scattering, recognition, detection, and tracking.
national Conference on Automatic Face and Gesture Boston: Artech House, 2002.
Recognition. 2008: 1-2. [29] Du L. Study on radar HRRP target recognition. PhD
[14] Wright J, Yang A Y, Ganesh A, et al. Robust face rec- thesis, Xi’an:Xidian University, 2007. [in Chinese]
ognition via sparse representation. IEEE Transactions [30] Li H J, Yang S H. Using range profiles as feature vec-
on Pattern Analysis and Machine Intelligence 2009; tors to identify aerospace objects. IEEE Transactions
31(2): 210-227. on Antennas and Propagation 1993; 41(3): 261-268.
[15] Gkalelis N, Tefas A, Pitas I. Sparse human movement
representation and recognition. IEEE 10th Workshop
Biographies:
Multimedia Signal Processing. 2008: 165-169.
[16] Hang X Y. Cancer classification by sparse representa-
tion using microarray gene expression data. IEEE In- Zhou Nuo Born in 1981, she is currently a Ph.D. candidate
ternational Conference on Bioinformatics and in Beijing University of Aeronautics and Astronautics. Her
Biomeidcine Workshops. 2008: 174-177. research interests include signal processing and automatic
[17] Mallat S G, Zhang Z F. Matching pursuits with target recognition.
time-frequency dictionaries. IEEE Transactions on Sig- E-mail: [email protected]
nal Processing 1993; 41(12): 3397-3415.
[18] Chen S S, Donoho D L, Saunders M A. Atomic de- Chen Wei Born in 1963ˈhe is an associate professor in
composition by basis pursuit. SIAM review 2001; Beijing University of Aeronautics and Astronautics. His
43(1): 129-159. research interests include signal and information processing,
[19] Donoho D L, Tsaig Y. Fast solution of l1-norm information fusion, digital image processing and target rec-
minimization problems when the solution may ognition.
be sparse. <http://www-dsp.rice.edu/files/cs/FastL1. E-mail: [email protected]
pdf >. 2006.