， ， ， ，

( 1.陜西師范大學(xué)計(jì)算機(jī)科學(xué)學(xué)院， 陜西西安 710119； 2.西安現(xiàn)代控制技術(shù)研究所， 陜西西安 710065； 3.西安電子科技大學(xué)雷達(dá)信號(hào)處理國家重點(diǎn)實(shí)驗(yàn)室， 陜西西安 710071)

0 引言

合成孔徑雷達(dá)(Synthetic Aperture Radar，SAR)由于其全天時(shí)、全天候,能夠穿透云、雨、煙、霧等優(yōu)點(diǎn)，已被廣泛地應(yīng)用于軍事和民用等諸多領(lǐng)域[1-3]。其中，基于SAR的自動(dòng)目標(biāo)識(shí)別(Automatic Target Recognition，ATR)技術(shù)是SAR應(yīng)用的一個(gè)研究熱點(diǎn)[4-5]，它的研究對(duì)于環(huán)境監(jiān)測(cè)、敵我識(shí)別和災(zāi)害評(píng)估等方面都有著重要意義。

傳統(tǒng)的SAR目標(biāo)識(shí)別算法針對(duì)目標(biāo)類型的識(shí)別，也就是說，同一類型不同型號(hào)的樣本在識(shí)別時(shí)被認(rèn)為是一樣的[6-8]。而目標(biāo)的型號(hào)識(shí)別比類型識(shí)別更具挑戰(zhàn)性，因?yàn)橥活愋湍繕?biāo)的不同型號(hào)之間的差異非常小，如同一類型坦克上有沒有機(jī)關(guān)槍和油箱、天線是否展開等。高精度的目標(biāo)型號(hào)識(shí)別對(duì)于戰(zhàn)場(chǎng)感知、精確打擊等需要獲取目標(biāo)詳細(xì)信息的應(yīng)用領(lǐng)域是非常必要的，有著重要的應(yīng)用價(jià)值[9-10]。

稀疏描述是一種有效的模式識(shí)別、機(jī)器學(xué)習(xí)理論，已被成功應(yīng)用于雷達(dá)成像[11]、圖像處理[12]、數(shù)據(jù)壓縮[13]、模式識(shí)別[14]等領(lǐng)域。稀疏描述用于識(shí)別時(shí)，利用所有訓(xùn)練樣本構(gòu)造字典矩陣，然后求解測(cè)試樣本的稀疏向量，最終根據(jù)稀疏向量計(jì)算重構(gòu)誤差得到識(shí)別結(jié)果?；谙∈杳枋龅腟AR目標(biāo)型號(hào)識(shí)別算法可以得到較好的識(shí)別結(jié)果[10]，但是由于陰影效應(yīng)、信號(hào)和環(huán)境的相互作用、將三維場(chǎng)景投影到平面等原因[8]，SAR圖像對(duì)目標(biāo)方位角的變化較為敏感。為了實(shí)現(xiàn)較好的SAR目標(biāo)型號(hào)識(shí)別，對(duì)訓(xùn)練樣本有較大的需求，通常需要覆蓋整個(gè)方位區(qū)間的大量訓(xùn)練樣本來全面捕獲目標(biāo)特征，這就使得基于稀疏描述的SAR目標(biāo)型號(hào)識(shí)別算法的計(jì)算復(fù)雜度較高，不利于算法的實(shí)時(shí)性和高效性。

針對(duì)此問題，本文提出一種快速稀疏描述(Efficient Sparse Representation，ESR)算法，實(shí)現(xiàn)SAR目標(biāo)的型號(hào)識(shí)別。考慮到雖然SAR圖像對(duì)于目標(biāo)方位角的變化較為敏感，但是在一定的角度范圍內(nèi)目標(biāo)圖像的形態(tài)變化不大，也即在一定角度范圍內(nèi)目標(biāo)SAR圖像惰性變化的特點(diǎn)[15]，所提算法將每個(gè)型號(hào)的訓(xùn)練樣本在一定的方位區(qū)間內(nèi)分別取平均，采用不同型號(hào)在不同方位區(qū)間的平均樣本表征原來方位區(qū)間內(nèi)的若干個(gè)樣本。這樣可有效減少訓(xùn)練樣本的數(shù)目，顯著降低算法的計(jì)算復(fù)雜度，提高基于稀疏描述的SAR目標(biāo)型號(hào)識(shí)別算法的效率。

1 基于稀疏描述的SAR目標(biāo)型號(hào)識(shí)別算法

采用稀疏描述理論進(jìn)行SAR目標(biāo)型號(hào)識(shí)別時(shí)，首先利用屬于第ii=1,2,…,C個(gè)型號(hào)的Ni個(gè)訓(xùn)練樣本構(gòu)造矩陣Ai,Ni為第i個(gè)型號(hào)的訓(xùn)練樣本個(gè)數(shù)，Ai=xi1,xi2,…,xiNi,xij表示由第i個(gè)型號(hào)的第jj=1,2,…,Ni個(gè)訓(xùn)練樣本形成的列向量。那么，所有型號(hào)目標(biāo)的訓(xùn)練樣本可以形成字典矩陣D=A1,A2,…,AC,對(duì)于測(cè)試樣本y來說，它可用字典矩陣D線性描述為

y=Da

(1)

式中，a表示對(duì)應(yīng)的系數(shù)向量，若測(cè)試樣本y屬于第i個(gè)型號(hào)，則理想情形下a中的系數(shù)應(yīng)只與第i個(gè)型號(hào)訓(xùn)練樣本相對(duì)應(yīng)的為非零值，其余均為零[14]：

aN×1=0,…,0,ai1,ai2,…,aiNi,0,…,0T

(2)

式中，N=N1+N2+…+NC為所有型號(hào)目標(biāo)的訓(xùn)練樣本總數(shù)，當(dāng)Ni?N時(shí)，此向量是稀疏的。

(3)

嚴(yán)格來說，信號(hào)的稀疏性應(yīng)該用0范數(shù)來加以刻畫，但是，0范數(shù)的求解是一個(gè)NP hard問題。針對(duì)此問題，通常有兩種常用的求解算法，一種是Chen，Candes等提出的將該問題松弛到1范數(shù)加以求解[16-17]。在一定條件下，其求解結(jié)果和0范數(shù)是一致的，1范數(shù)的求解是凸優(yōu)化問題，可以采用基追蹤(Basis Pursuit，BP)算法進(jìn)行求解[16]；另一種是采用Tropp等提出的基于貪婪迭代的正交匹配追蹤(Orthogonal Matching Pursuit，OMP)算法[18]進(jìn)行求解。

(4)

2 基于快速稀疏描述的SAR目標(biāo)型號(hào)識(shí)別算法

圖1 不同目標(biāo)在不同方位角下的SAR圖像

如前所述，SAR圖像對(duì)于目標(biāo)方位角的變化較為敏感，但是對(duì)于某一確定方位的目標(biāo)圖像，在一定角度范圍內(nèi)目標(biāo)SAR圖像具有惰性變化的特點(diǎn)[15]，如圖2所示，以圖2中第二行的目標(biāo)BMP2-9566為例，其在161°,162°,163°和165°的SAR圖像差異并不大。根據(jù)SAR圖像的此特點(diǎn)，本文所提算法將每個(gè)型號(hào)的訓(xùn)練樣本在一定的方位區(qū)間內(nèi)分別取平均，采用不同方位區(qū)間的平均樣本表征原來方位區(qū)間內(nèi)的若干個(gè)樣本。具體來說，將0°～360°整個(gè)方位區(qū)間劃分成K塊，則每塊覆蓋360°/K的角度范圍，將矩陣Ai落在每個(gè)劃分塊內(nèi)的所有訓(xùn)練樣本取平均，即

(5)

(6)

(7)

圖2 不同目標(biāo)在5°方位角區(qū)間內(nèi)的SAR圖像

3 實(shí)驗(yàn)結(jié)果與分析

利用美國MSTAR數(shù)據(jù)進(jìn)行實(shí)驗(yàn)驗(yàn)證所提算法的有效性[8-10]，選取其中3類7種不同型號(hào)的目標(biāo)，目標(biāo)SAR圖像的大小為128像素×128像素，雷達(dá)工作在X波段，采用聚束模式進(jìn)行成像，分辨率為0.3 m×0.3 m，方位覆蓋范圍為0°～360°。表1給出了訓(xùn)練樣本和測(cè)試樣本的型號(hào)和樣本個(gè)數(shù)。

表1 訓(xùn)練樣本與測(cè)試樣本的型號(hào)及樣本個(gè)數(shù)

首先在BMP2數(shù)據(jù)集上進(jìn)行實(shí)驗(yàn)，以5°為一個(gè)方位區(qū)間對(duì)訓(xùn)練樣本進(jìn)行取平均操作，特征維數(shù)為1 024。選取SR-BP、SR-OMP、支持向量機(jī)(Support Vector Machine，SVM)、K近鄰分類器(K Nearest Neighbor，K-NN)四種不同的算法在此數(shù)據(jù)集上進(jìn)行實(shí)驗(yàn)，以驗(yàn)證所提算法的有效性。SR-BP為采用BP算法求解稀疏向量的SR算法，SR-OMP為采用OMP算法求解稀疏向量的SR算法，ESR-BP和ESR-OMP為對(duì)應(yīng)本文所提的快速SR算法。不同算法下的實(shí)驗(yàn)結(jié)果如表2所示。由識(shí)別結(jié)果可以看出，本文所提算法的效果要弱于SR，可見所提快速稀疏描述算法在降低算法復(fù)雜度的同時(shí)也損失了一部分有效信息，但是由于SAR圖像在一定的角度范圍內(nèi)惰性變化的特點(diǎn)，信息的損失較小，因此其識(shí)別結(jié)果仍然要優(yōu)于SVM算法和K-NN算法。由表2可見，所提快速稀疏描述算法在大幅度降低算法復(fù)雜度的前提下，能達(dá)到87.22%的平均識(shí)別率。

下面采用T72數(shù)據(jù)集驗(yàn)證所提快速稀疏描述算法的有效性，實(shí)驗(yàn)條件和上述相同，對(duì)應(yīng)的實(shí)驗(yàn)結(jié)果如表3所示?？梢钥吹?，在此條件下，所提快速稀疏描述算法的性能仍然較為優(yōu)越，可以獲得94.85%的正確識(shí)別率。由實(shí)驗(yàn)結(jié)果可見，所提算法在T72數(shù)據(jù)集上同樣可以取得滿意的識(shí)別結(jié)果。

表2 BMP2目標(biāo)型號(hào)識(shí)別結(jié)果 %

表3 T72目標(biāo)型號(hào)識(shí)別結(jié)果 %

最后，將所提算法在表1所示的全部7個(gè)型號(hào)目標(biāo)上進(jìn)行實(shí)驗(yàn)，相應(yīng)的實(shí)驗(yàn)結(jié)果如表4所示。可以看到，在多型號(hào)數(shù)據(jù)條件下，所提快速稀疏描述算法依然可以取得滿意的識(shí)別結(jié)果。由實(shí)驗(yàn)結(jié)果可見，SAR圖像雖然具有明顯的方位角敏感特性，但其在一定的范圍內(nèi)圖像具有惰性變化的特點(diǎn)，利用區(qū)間內(nèi)的平均樣本可較好地表征這一區(qū)間內(nèi)的樣本特性。雖然這樣會(huì)造成一部分識(shí)別信息的損失，但損失較小。在需要快速識(shí)別的情形下，可考慮采用本文所提的快速稀疏描述識(shí)別算法。圖3為不同算法對(duì)各個(gè)型號(hào)目標(biāo)的識(shí)別結(jié)果圖。

表4 7個(gè)型號(hào)目標(biāo)識(shí)別結(jié)果 %

圖3 不同算法對(duì)各個(gè)型號(hào)的識(shí)別結(jié)果

表5和表6分別給出了ESR-BP算法和ESR-OMP算法在不同方位區(qū)間下的算法復(fù)雜度，可以看出，樣本的方位區(qū)間取得越大，算法的計(jì)算復(fù)雜度就越低，但相應(yīng)的識(shí)別結(jié)果也會(huì)變差，對(duì)應(yīng)的識(shí)別結(jié)果如圖4所示。由圖4可見，在方位區(qū)間取為2°的時(shí)候，識(shí)別結(jié)果可達(dá)到91.28%，當(dāng)方位區(qū)間增大到12°的時(shí)候，識(shí)別精度會(huì)下降到86.08%。因此，實(shí)際應(yīng)用中需要對(duì)識(shí)別精度和算法復(fù)雜度進(jìn)行折中考慮。

表5 ESR-BP算法在不同方位區(qū)間下的算法復(fù)雜度

表6 ESR-OMP算法在不同方位區(qū)間下的算法復(fù)雜度

圖4 所提快速稀疏描述算法在不同方位區(qū)間下的實(shí)驗(yàn)結(jié)果

4 結(jié)束語

利用SAR圖像在一定角度范圍內(nèi)惰性變化的特點(diǎn)，本文提出了一種快速稀疏描述算法，有效地降低了基于稀疏描述的SAR目標(biāo)型號(hào)識(shí)別的算法復(fù)雜度。方位區(qū)間設(shè)置越小，則目標(biāo)信息保留越完整，算法復(fù)雜度也就相應(yīng)越高，但識(shí)別結(jié)果越好；方位區(qū)間設(shè)置越大，則目標(biāo)信息損失越多，識(shí)別結(jié)果有所降低，但對(duì)應(yīng)的算法復(fù)雜度也越低，在實(shí)際應(yīng)用中可根據(jù)具體要求選擇合適的識(shí)別算法。