何興宇，童寧寧，賀吉峰

(1．空軍工程大學(xué)防空反導(dǎo)學(xué)院，陜西 西安710051;2．駐航天科工集團(tuán)二院23所軍代表室，北京100854)

1 引言

彈道目標(biāo)特征提取是中段彈道目標(biāo)識別的基礎(chǔ)和關(guān)鍵。彈道目標(biāo)的特征主要包括微動特征和結(jié)構(gòu)特征，微動特征如進(jìn)動角、進(jìn)動周期等，結(jié)構(gòu)特征如尺寸、形狀等［1－3］。利用ISAR實現(xiàn)高分辨成像是進(jìn)行目標(biāo)結(jié)構(gòu)特征提取的重要途徑，而ISAR圖像的橫向定標(biāo)是提取目標(biāo)結(jié)構(gòu)特征的難點(diǎn)和關(guān)鍵。中段彈道目標(biāo)在轉(zhuǎn)臺模型下通常是均勻旋轉(zhuǎn)的，其ISAR圖像定標(biāo)的關(guān)鍵是成像積累角及轉(zhuǎn)速的估計。由于ISAR成像是小角度成像，通常成像積累角不超過10°，依靠單次ISAR像難以完成轉(zhuǎn)速的準(zhǔn)確估計［4－6］。

目前，大多數(shù)橫向定標(biāo)算法都是基于多特顯點(diǎn)法，篩選出包含特顯點(diǎn)的距離單元，通過提取這些單元內(nèi)特顯點(diǎn)的回波信號估計調(diào)頻率。文獻(xiàn)［7］通過圖像分割劃分強(qiáng)散射區(qū)域來實現(xiàn)橫向定標(biāo)，該方法利用LPFT(Local Polynomial Fourier Transform，局域多項式傅里葉變換)等算法，實現(xiàn)散射中心的提取和調(diào)頻斜率的估計，有較高的估計精度，然而圖像分割門限確定不準(zhǔn)確將直接影響定標(biāo)精度，算法穩(wěn)定性不高。文獻(xiàn)［7］指出，ISAR圖像散射中心相對位置不因目標(biāo)的三維轉(zhuǎn)動而發(fā)生變化，并據(jù)此提出了基于圖像配準(zhǔn)的彈道目標(biāo)ISAR像橫向定標(biāo)算法，該算法同樣需要尋找合適的控制點(diǎn)，通過仿射變換實現(xiàn)ISAR像橫向定標(biāo)。但是該方法敏感于橫向伸縮因子的估計，且計算較復(fù)雜。本文根據(jù)文獻(xiàn)［8］分析得到的ISAR圖像散射中心相對位置不變的特點(diǎn)，通過提取ISAR圖像序列中由散射點(diǎn)構(gòu)成的模值最大的向量，根據(jù)彈道目標(biāo)的旋轉(zhuǎn)對稱性，得到目標(biāo)對稱軸方向的單位向量表示，并估計該向量在成像序列之間的轉(zhuǎn)角及轉(zhuǎn)速，最后完成ISAR像的橫向定標(biāo)。仿真實驗表明，本文的ISAR圖像定標(biāo)算法能準(zhǔn)確地完成目標(biāo)轉(zhuǎn)速估計，定標(biāo)效果較好。

2 ISAR圖像定標(biāo)原理

逆合成孔徑雷達(dá)利用發(fā)射寬帶信號達(dá)到距離向的高分辨，依靠長的合成孔徑實現(xiàn)方位向的高分辨，而距離—多普勒算法是最常用的ISAR成像算法。

雷達(dá)工作在線性調(diào)頻體制下。設(shè)目標(biāo)到雷達(dá)的距離為Ri，參考距離為Rref。解調(diào)后的雷達(dá)回波為［9－11］:

對得到的解調(diào)后的回波以fs采樣頻率進(jìn)行采樣，設(shè)每次回波可以得到N個采樣點(diǎn)，做Nf點(diǎn)離散傅里葉變換，可得目標(biāo)的一維距離像:

設(shè)目標(biāo)平動補(bǔ)償后，轉(zhuǎn)動速度為ω，則散射點(diǎn)與參考中心的距離為:

去斜后，在慢時間域作Mf點(diǎn)離散傅里葉變換實現(xiàn)橫向壓縮，得到目標(biāo)的ISAR像為:

式中，M表示脈沖積累數(shù);fp為脈沖重復(fù)頻率;(nf，mf)表示尖峰脈沖位置。

由上式可以看出，在完成距離像和橫向壓縮后，在 ( nf，mf)域會出現(xiàn)sinc函數(shù)形式的尖峰脈沖，對應(yīng)了散射點(diǎn)的位置。尖峰脈沖位置可以用下式表示:

3 轉(zhuǎn)動角速度估計

3．1 模型建立

設(shè)ISAR像的縱向分辨率和橫向分辨率分別為ρr和ρa(bǔ)，若不發(fā)生越距離單元徙動，則目標(biāo)的最大徑向尺寸和橫向尺寸應(yīng)滿足以下條件［12］:

式中，λ為雷達(dá)發(fā)射信號波長;若取λ=4 cm，ρr和ρa(bǔ)均為0．4 m，則目標(biāo)尺寸應(yīng)不大于16 m×16 m，彈道目標(biāo)通常都能滿足這一要求，因此本文成像算法中假設(shè)散射點(diǎn)不發(fā)生越距離單元徙動。

運(yùn)動目標(biāo)回波信號經(jīng)包絡(luò)對齊和相位補(bǔ)償，將最終轉(zhuǎn)為轉(zhuǎn)臺模型進(jìn)行成像。假設(shè)目標(biāo)完成平動補(bǔ)償后，距離變化可以消除，同時，仿真假設(shè)識別的彈道目標(biāo)為錐柱體。

由于ISAR成像是小角度成像，加上干擾噪聲等因素的影響，依靠單次ISAR像進(jìn)行目標(biāo)轉(zhuǎn)速估計有較大的不確定性。本文提出利用ISAR像序列估計目標(biāo)轉(zhuǎn)動角速度，對消干擾噪聲因素的影響，提高定標(biāo)精度。其模型如圖1所示。

圖1 ISAR像序列目標(biāo)位置及轉(zhuǎn)角

其中，A、B、C、E、F為目標(biāo)上的五個散射點(diǎn)，AD為目標(biāo)的對稱軸。由目標(biāo)ISAR像的散射點(diǎn)構(gòu)成的向量有多個，考察ISAR像中模最大的向量。由于彈道目標(biāo)的旋轉(zhuǎn)對稱性，通常模值最大的向量有兩個，可以通過它們來確定對稱軸的方向，進(jìn)而估計成像序列間目標(biāo)的轉(zhuǎn)角。

3．2 算法分析與設(shè)計

假設(shè)散射點(diǎn)在ISAR像中坐標(biāo)分別為(xi，yi)，i∈［A，B，C，E，F(xiàn)］。得到各向量的坐標(biāo)表示，如AB=(xB－xA，yB－yA)。可以看出，向量只與散射點(diǎn)的相對位置有關(guān)。

從圖1可以看出，模值最大的向量在AC和AE中產(chǎn)生。如果圖像聚焦效果不好，將影響向量模值的估計和目標(biāo)轉(zhuǎn)角的估計精度，進(jìn)而影響目標(biāo)的定標(biāo)精度。因此，對AC和AE向量作以下處理:

首先，判斷式(8)是否成立。式中，δr為允許的偏差值，取δr=1，即兩向量的模值差不大于一個橫向(徑向)距離單元。

當(dāng)式(8)成立時，判定兩向量模值相同，通過計算兩向量的合成向量在成像序列中的變化來估計轉(zhuǎn)動角;但是，在目標(biāo)的一定姿態(tài)角下，式(8)并不成立，此時需提取AC和AE中模值大的向量，考察它在成像序列中的變化來估計轉(zhuǎn)動角。本文后面的分析均假設(shè)式(8)成立，其他情況可按下面分析作類似處理。

設(shè)初始時刻對稱軸AD方向的單位向量為e1，則有下式成立:

經(jīng)歷時間t后，在下個成像時刻，對稱軸方向的單位向量設(shè)為e2，則e1與e2的夾角即為在時間t內(nèi)目標(biāo)轉(zhuǎn)過的角度θ(t)。有下式成立:

通過式(10)，提取目標(biāo)的轉(zhuǎn)動角，估計目標(biāo)轉(zhuǎn)速，進(jìn)而根據(jù)式(6)完成目標(biāo)的橫向定標(biāo)。

4 仿真分析

假設(shè)目標(biāo)由五個散射點(diǎn)組成，其位置關(guān)系如圖2所示。

圖2 各散射點(diǎn)位置關(guān)系

雷達(dá)發(fā)射的線性調(diào)頻信號載頻fc=6 GHz，波長λ=0．03 m，脈沖寬度Tp=40μs，帶寬B=300 MHz，調(diào)頻斜率k=B/Tp，采樣頻率Fs=300 MHz。設(shè)轉(zhuǎn)臺模型逆時針方向旋轉(zhuǎn)為正，轉(zhuǎn)速設(shè)置為ω0=0．01 rad/s。

初始成像時刻為t1=0時，根據(jù)目標(biāo)初始時刻的位置，得到目標(biāo)ISAR像如圖3所示。

圖3 t1=0 s時目標(biāo)ISAR像

各散射點(diǎn)在ISAR像中的坐標(biāo)如下:A(207，122)，B(205，126)，C(203，128)，E(201，126)，F(xiàn)(203，124)。模值最大的向量為AC=(－4，6)及AE=(－6，4)，根據(jù)式(9)計算得

二次成像時刻為t2=40 s時，由于轉(zhuǎn)臺的轉(zhuǎn)動，各散射點(diǎn)的位置發(fā)生了變化，如圖4所示。

圖4 t2=40 s時目標(biāo)ISAR像

由t2=40 s時的ISAR像可得AC=(－6，4)，AE=(－7，1)，=0．14≤1，e2==(－0．9333，0．3590)。

目標(biāo)在40 s的時間里轉(zhuǎn)過的角度為θ(t)=arccos e1·e2=0．4182 rad，計算得目標(biāo)轉(zhuǎn)速ω==0．0105 rad/s，與仿真設(shè)置的目標(biāo)轉(zhuǎn)速ω0=0．01 rad/s相當(dāng)，相對誤差=4．5561%，轉(zhuǎn)速估計效果較好。

將得到的轉(zhuǎn)速估計值帶入式(6)，對t1=0時的ISAR像定標(biāo)，結(jié)果如圖5所示。

從圖5中的定標(biāo)結(jié)果可以看出，各散射點(diǎn)的相對位置關(guān)系與散射點(diǎn)模型基本相同，定標(biāo)效果較好，精度較高。定標(biāo)得到的散射點(diǎn)位置與原始模型關(guān)于原點(diǎn)對稱，這是由于多普勒符號引起的。

圖5 ISAR定標(biāo)后圖像

為了驗證算法性能，選取不同的成像序列間隔，得到轉(zhuǎn)角、轉(zhuǎn)速、轉(zhuǎn)速估計誤差及橫向定標(biāo)誤差如表1所示。

表1 不同序列間隔下的仿真實驗結(jié)果Tab．1 Simulation results in different time interval

從表1可看出，除了初次實驗，轉(zhuǎn)速估計誤差和橫向定標(biāo)誤差均在5%以下，估計誤差較小，橫向定標(biāo)精度較高。從表中還可以看出，在時間間隔較小時，轉(zhuǎn)角及轉(zhuǎn)速估計值誤差較高;但是時間間隔的選取也不能過大，否則會由于ISAR像序列中向量相對關(guān)系和目標(biāo)散射特性的變化而導(dǎo)致誤差急劇增大。另外，還可以通過選取多個時間間隔來估計目標(biāo)轉(zhuǎn)速，增加估計和定標(biāo)精度。

5 結(jié)論

彈道目標(biāo)ISAR圖像橫向定標(biāo)在彈道目標(biāo)識別中有重要地位。本文從經(jīng)典的距離—多普勒成像算法出發(fā)，分析了ISAR成像的特點(diǎn)及利用成像序列估計轉(zhuǎn)速進(jìn)而完成橫向定標(biāo)的優(yōu)點(diǎn)。通過對ISAR像序列中散射點(diǎn)構(gòu)成的向量的分析，得到轉(zhuǎn)角及轉(zhuǎn)速估計表達(dá)式，完成彈道目標(biāo)的橫向定標(biāo)。仿真實驗證明了本文方法簡單有效，定標(biāo)精度較高。本文研究了兩幅ISAR像實現(xiàn)橫向定標(biāo)的過程，當(dāng)然還可以通過利用多個ISAR像序列綜合來提高定標(biāo)精度。

［1］ Rihaczk A W，Hershkowitz SJ．Theory and practice of radar target identification［M］．Boston:Artech House，2000．

［2］ Dieter Mehrholz．Radar techniques for the characterization of meter－sized object in space［J］．Adv Space Res，2001，28(9):1259－1268．

［3］ MA Liang，WANG Tao，F(xiàn)ENGDejun，et al．The characteristic of range profile and micro－motion feature extraction for rotary target［J］．Acta Electronica Sinica，2008，36(12):2273－2279．(in Chinese)馬梁，王濤，馮德軍，等．旋轉(zhuǎn)目標(biāo)距離像長度特性及微運(yùn)動特征提?。跩］．電子學(xué)報，2008，36(12):2273－2279．

［4］ Thayaparan T，F(xiàn)rangos P，Stankovic L．Signal processing techniques for ISAR and feature extraction［J］．IET signal processing，2008，2(3):189－191．

［5］ JLi，H Ling．Application of adaptive chirplet representation for ISAR feature extraction from target with rotating parts［J］．IEE Proc－Radar Sonar Naving，2003，150(4):284－291．

［6］CHEN Xiuwei，ZHANG Yunhua，ZHANG Xiangkun．Ballistic missile echo simulation based on ISAR［J］．Journal of Air Force Engineering University:Natural Science Edition，2011，12(1):29－34．(in Chinese)陳秀偉，張云華，張祥坤．彈道導(dǎo)彈的ISAR回波模擬［J］．空軍工程大學(xué)學(xué)報:自然科學(xué)版，2011，12(1):29－34．

［7］ Martorella M．Novel approach for ISAR image crossrange scaling［J］．IEEETrans on Aerospace and Electronic Systems，2008，44(1):281－294．

［8］ JIN Guanghu，GAO Xunzhang，LI Xiang，et al．ISAR image cross scaling method for ballistic target based on image registration［J］．Systems Engineering and Electronics，2010，32(12):2565－2569．(in Chinese)金光虎，高勛章，黎湘，等．基于圖像配準(zhǔn)的彈道目標(biāo)ISAR圖像橫向定標(biāo)［J］．系統(tǒng)工程與電子技術(shù)，2010，32(12):2565－2569．

［9］ BAO Zheng，XING Mengdao，WANG Tong．Radar imaging technology［M］．Beijing:Electronic Industry Press，2005．(in Chinese)保錚，邢孟道，王彤．雷達(dá)成像技術(shù)［M］．北京:電子工業(yè)出版社，2005．

［10］QUN Zhang，Tat Soon Yeo，Hwee Siang Tan，et al．Imaging of a moving target with a rotating parts based on the Hough transform［J］．IEEE Transactions on Geoscience and Remote Sensing，2008，46(1):291－299．

［11］Jin Guanghu．Research on ISAR imaging and physical feature extraction of midcourse ballistic target［D］．Changsha:National University of Defense Technology，2009．(in Chinese)金光虎．中段彈道目標(biāo)ISAR成像及物理特征反演技術(shù)研究［D］．長沙:國防科學(xué)技術(shù)大學(xué)，2009．

［12］HU Jiemin，ZHANG Jun，ZHAN Ronghui，et al．Novel method for estimating rotation angle based on phase cancellation［J］．Systems Engineering and Electronics，2012，34(5):897－902．(in Chinese)胡杰民，張軍，占榮輝，等．一種基于相位對消的轉(zhuǎn)角估計新方法［J］．系統(tǒng)工程與電子技術(shù)，2012，34(5):897－902．