楊躍輪

(上海地區(qū)艦炮系統(tǒng)軍代表室，　上海 201913)

?

機(jī)載SAR對地面慢速運(yùn)動目標(biāo)的DPCA-CFAR聯(lián)合檢測方法

楊躍輪

(上海地區(qū)艦炮系統(tǒng)軍代表室，上海 201913)

摘要：在機(jī)載合成孔徑雷達(dá)(SAR)動目標(biāo)檢測中，由于地面慢速運(yùn)動目標(biāo)速度較小，因此很容易被地面強(qiáng)主瓣雜波淹沒而檢測不到。文中針對地面慢速運(yùn)動目標(biāo)的運(yùn)動特點(diǎn)，提出基于分?jǐn)?shù)階傅里葉變換(Frft)的DPCA-CFAR聯(lián)合檢測的方法，在采用天線相位中心偏置(DPCA)技術(shù)進(jìn)行雜波抑制的基礎(chǔ)上進(jìn)行恒虛警(CFAR)處理，從而檢測到運(yùn)動目標(biāo)，并且根據(jù)DPCA對消后的信號幅度和CFAR檢測門限推導(dǎo)了最小可檢測速度，說明了提出的算法對慢速運(yùn)動目標(biāo)的檢測性能，并進(jìn)一步采用Frft方法估計(jì)出目標(biāo)速度，最后通過仿真對算法的有效性進(jìn)行了驗(yàn)證。

關(guān)鍵詞：機(jī)載SAR；地面慢速運(yùn)動目標(biāo)；DPCA-CFAR；傅里葉變換

0引言

機(jī)載SAR動目標(biāo)的檢測方法主要有兩種：單通道方法和多通道方法。單通道方法是根據(jù)目標(biāo)的回波特性，估計(jì)目標(biāo)的多普勒中心頻率和調(diào)頻率，進(jìn)而得到目標(biāo)的速度，但是目標(biāo)的方位向位置不能確定，而多通道方法則能確定目標(biāo)的方位位置，但不一定能檢測到目標(biāo)速度[1]。

通常我們將速度低于12 km/h的目標(biāo)稱為慢速運(yùn)動目標(biāo)，由于地物雜波的存在，此類目標(biāo)很容易淹沒在強(qiáng)主瓣雜波中而檢測不到，所以對于地面慢速運(yùn)動目標(biāo)的檢測，關(guān)鍵在于有效抑制雜波，提高信雜比。單通道檢測方法主要有頻域?yàn)V波法，參數(shù)估計(jì)法和時(shí)頻分析法，由于慢速運(yùn)動目標(biāo)頻譜會被雜波掩蓋，因此，不能采用頻域?yàn)V波法來檢測，參數(shù)估計(jì)法和時(shí)頻分析法

的主要作用都是估計(jì)目標(biāo)的多普勒參數(shù)。多通道檢測方法有天線相位中心偏置(DPCA)、沿航跡干涉(ATI)等，DPCA是通過將多通道SAR圖像做差抑制雜波而檢測到動目標(biāo)，ATI是通過對多通道SAR圖像做干涉處理，對干涉相位進(jìn)行門限判斷而區(qū)分動靜目標(biāo)，并且確定目標(biāo)方位位置。綜合以上分析，可以得出，對于地面慢速運(yùn)動目標(biāo)的檢測，可以將多通道方法與單通道方法相結(jié)合，多通道方法主要用來進(jìn)行雜波抑制并得到運(yùn)動目標(biāo)信息，而單通道方法可以用來估計(jì)運(yùn)動目標(biāo)的參數(shù)[2]。

1基于Frft的DPCA-CFAR聯(lián)合檢測

1.1DPCA-CFAR聯(lián)合檢測原理

雙通道DPCA能夠抑制雜波并濾除靜止目標(biāo)，從而實(shí)現(xiàn)機(jī)載SAR對地面目標(biāo)的檢測，對于快速運(yùn)動目標(biāo)和靜止目標(biāo)背景(即不存在雜波)下的慢速運(yùn)動目標(biāo)具有很好的檢測效果，但是在雜波環(huán)境下，即使采用

DPCA處理，也會有殘余雜波，干擾對慢速目標(biāo)的檢測，在此采用DPCA-CFAR聯(lián)合檢測的方法，在對兩個(gè)通道回波信號進(jìn)行DPCA雜波對消，信雜比已提高的情況下，再設(shè)置合理的檢測門限進(jìn)行CFAR處理，從而檢測到運(yùn)動目標(biāo)。最后可采用時(shí)頻分析Frft方法估計(jì)目標(biāo)多普勒中心頻率和調(diào)頻率，進(jìn)一步得到目標(biāo)速度估計(jì)值。

DPCA-CFAP聯(lián)合檢測系統(tǒng)結(jié)構(gòu)，如圖1所示，飛機(jī)速度為va，在飛機(jī)運(yùn)動方向依次放置兩個(gè)相距為d的天線，坐標(biāo)分別為(vat+d,0,h)，(vat,0,h)，由天線1發(fā)射LFM信號，兩個(gè)天線同時(shí)接收回波信號。P點(diǎn)為地面動目標(biāo)，坐標(biāo)為(x0,y0,0)，方位向和地面距離向速度分別為vx和vy，R1(tm)和R2(tm)分別為tm時(shí)刻兩天線與地面動目標(biāo)的距離[3]。

圖1　系統(tǒng)結(jié)構(gòu)

(1)

(2)

兩個(gè)天線接收到的回波信號分別為

exp[jπγ(t′-2R1(tm)/c)2]

(3)

exp[jπγ(t′-2R2(tm)+R1(tm)/c)2]

(4)

對接收到的回波信號進(jìn)行距離壓縮后，得到的信號分別為

S1(t′,tm)=στPsinc[πB(t′-2R1(tm)/c)]

(5)

S2(t′,tm)=στPsinc[πB(t′-(R2(tm)+R1(tm))/c]

(6)

式中：σ為目標(biāo)反射系數(shù)；τP為發(fā)射信號脈沖寬度；B為發(fā)射信號帶寬；Tsar為合成孔徑時(shí)間。

對于慢速運(yùn)動目標(biāo)，有vx?va，vy?va，因此可認(rèn)為R1(tm)≈R2(tm)，對通道1和2的信號進(jìn)行圖像插值配準(zhǔn)，得到幅度對消信號

S12(t′,tm)=S1(t′,tm)·C1-S2(t′,tm+k·PRTr)=

(7)

(8)

(9)

對于靜止目標(biāo)，有vx=vy=0，則S12(t′,tm)=0，說明靜止目標(biāo)回波經(jīng)過DPCA處理后被抵消掉[4]。

雙通道的DPCA技術(shù)主要用于實(shí)現(xiàn)雜波的對消與動目標(biāo)的檢測，但是在實(shí)際情況中，由于各種非理想因素的存在，經(jīng)過DPCA處理后仍有雜波殘余，對于慢速運(yùn)動目標(biāo)的檢測有一定影響。在經(jīng)過DPCA處理后，信雜比已提高的情況下，再設(shè)定合適的檢測門限，采用CFAR恒虛警檢測，就可以檢測動目標(biāo)及其所在的距離單元。

(10)

1.2慢速運(yùn)動目標(biāo)檢測性能

選取合適的虛警概率，采用CFAR對經(jīng)過DPCA處理后的回波信號進(jìn)行檢測，可以檢測到運(yùn)動目標(biāo)及其所在的距離單元。我們認(rèn)為雜波服從瑞利分布，則虛警概率Pf和所需的檢測門限k具有以下關(guān)系

(11)

設(shè)定了虛警概率后，就可以得到檢測門限值。設(shè)DPCA處理后被檢測距離單元的信號幅度為A，根據(jù)式(7)得到

(12)

(13)

(14)

圖2　最小可檢測速度

1.3Frft估計(jì)目標(biāo)多普勒參數(shù)

分?jǐn)?shù)階Frft定義如下[7]

(15)

Kα(t,u)=

(16)

Frft的本質(zhì)是對信號的旋轉(zhuǎn)，選擇合適的旋轉(zhuǎn)角度對LFM信號進(jìn)行Frft變換，可使信號能量在某一特定的分?jǐn)?shù)階上聚集，幅度出現(xiàn)明顯的峰值，而在其他階數(shù)上幅度為零，因此能夠準(zhǔn)確地檢測和估計(jì)運(yùn)動目標(biāo)參數(shù)[8-9]。

對通道1和通道2的DPCA信號，忽略快時(shí)間sinc函數(shù)，只保留慢時(shí)間信號得到S12(tm)，進(jìn)行Frft得

j2π(fd-ucscα)tm]dtm

(17)

BmTsar·sinc[πTsar(fd-ucscα)]

(18)

Frft估計(jì)出目標(biāo)的多普勒中心頻率和調(diào)頻率如下

(19)

(20)

聯(lián)合式(10)、式(12)、式(19) 和式(20)，即可估計(jì)到目標(biāo)的運(yùn)動速度。

2仿真結(jié)果

2.1靜止目標(biāo)背景下的DPCA檢測

在不存在雜波，只存在靜止目標(biāo)和慢速運(yùn)動目標(biāo)的情況下，采用DPCA方法就可以實(shí)現(xiàn)對慢速運(yùn)動目標(biāo)的檢測，下面進(jìn)行仿真驗(yàn)證。仿真參數(shù)如表1所示，設(shè)置9個(gè)點(diǎn)目標(biāo)，目標(biāo)位置坐標(biāo)及速度值設(shè)置如表2所示，其中5個(gè)為慢速運(yùn)動目標(biāo)，4個(gè)為靜止目標(biāo)，DPCA仿真結(jié)果如圖3所示?？梢钥闯觯?jīng)過DPCA處理后，靜止目標(biāo)都被過濾掉，只顯示動目標(biāo)。仿真結(jié)果說明，在只有靜止目標(biāo)存在而無地物雜波的環(huán)境下，DPCA可以檢測到慢速運(yùn)動目標(biāo)。

表1　DPCA仿真參數(shù)

表2　目標(biāo)參數(shù)

圖3　DPCA檢測結(jié)果

2.2瑞利雜波環(huán)境下的DPCA-CFAR聯(lián)合檢測

在存在雜波時(shí)，對DPCA-CFAR聯(lián)合檢測的效果進(jìn)行分析，圖4a)為仿真DPCA處理后的點(diǎn)目標(biāo)圖像，可以看出，在存在瑞利雜波的環(huán)境下，經(jīng)過DPCA雜波對消后，雖然目標(biāo)能夠被檢測出來，但是仍有殘余雜波，使得檢測效果并不理想，對經(jīng)過DPCA處理后的信號再進(jìn)行CFAR處理得到如圖4b)所示，檢測效果得到了改善，能夠較為清晰地檢測到目標(biāo)。

圖4　DPCA-CFAR檢測結(jié)果

圖5　Frft檢測結(jié)果

3結(jié)束語

本文針對地面慢速運(yùn)動目標(biāo)易被地物雜波掩蓋的特性，提出在采用雙通道DPCA進(jìn)行雜波抑制后仍有雜波殘余的情況下，采用CFAR處理，通過設(shè)置較低的虛警門限可以進(jìn)一步雜抑制雜波，從而檢測到慢速目標(biāo)，最后采用Frft方法估計(jì)出目標(biāo)速度。并且推導(dǎo)了基于CFAR門限的最小可檢測速度，說明了該算法對慢速運(yùn)動目標(biāo)的檢測性能。

參 考 文 獻(xiàn)

[1]李瑞英. 合成孔徑雷達(dá)動目標(biāo)檢測與成像[D]．哈爾濱：哈爾濱工業(yè)大學(xué), 2010．

LI Ruiying. The detection and imaging of moving targets based on SAR[D]. Harbin：Harbin Institute of Technology, 2010.

[2]尉志平. 機(jī)載合成孔徑雷達(dá)地面運(yùn)動目標(biāo)檢測與成像技術(shù)[D]. 武漢：武漢理工大學(xué), 2012．

WEI Zhiping. The detection and imaging of ground moving targets with airborne synthetic aperture radar[D]. Wuhan: Wuhan University of Technology, 2012.

[3]鄭明潔, 楊汝良.一種改進(jìn)的DPCA運(yùn)動目標(biāo)檢測方法[J]. 電子學(xué)報(bào), 2004, 32(9): 1429-1432.

ZHENG Mingjie, YANG Ruliang. An improved DPCA moving targets detecting algorithm[J]. Acta Electronica Sinica, 2004, 32(9):1429-1432.

[4]陳偉志. 干涉合成孔徑雷達(dá)高度計(jì)與壓縮感知技術(shù)研究[D]. 南京：南京理工大學(xué), 2010．

CHEN Weizhi. Research on interferometry synthetic aperture radar altimeter and compressive sensing technology[D]. Nanjing: Nanjing University of Science and Technology, 2010.

[5]孫華東. 多通道SAR地面運(yùn)動目標(biāo)檢測與參數(shù)估計(jì)研究[D]. 哈爾濱：哈爾濱工業(yè)大學(xué), 2009．

SUN Huadong. Research on multi-channel SAR ground moving target detection and parameter estimations[D]. Harbin：Harbin Institute of Technology, 2009.

[6]袁昌, 楊麗春, 顏衛(wèi), 等. 一種改進(jìn)的恒虛警率SAR圖像艦船檢測算法[J]. 現(xiàn)代雷達(dá), 2012, 34(8): 29-32.

YUAN Chang, YANG Lichun, YAN Wei, et al. An improved SAR shirp detection algorithm based on constant false alarm rate[J]. Modern Radar, 2012, 34(8): 29-32.

[7]鄧兵, 陶然, 曲長文. 分?jǐn)?shù)階Fourier域中多分量chirp信號的遮蔽分析[J]. 電子學(xué)報(bào), 2007, 35(6): 1094-1098.

DENG Bing, TAO Ran, QU Changwen. Analysis of the shading between multicomponent chirp signals in the fractional Fourier domain[J]. Acta Electronica Sinica, 2007, 35(6):1094-1098.

[8]劉明敬,劉剛. 基于FRFT的機(jī)動目標(biāo)ISAR成像算法[J]. 現(xiàn)代雷達(dá), 2009, 31(12): 49-52.

LIU Mingjing, LIU Gang. ISAR imaging algorithm of maneuvering targets based on FRFT[J]. Modern Radar, 2009, 31(12): 49-52.

[9]ZHANG X P, LIAO G S, ZHU S Q, et al. Geometry-information-aided efficient motion parameter estimation for moving-target imaging and location[J]. IEEE Transactions on Geoscience and Romote Sensing Letters, 2015,12(1): 155-159．

楊躍輪男，1972年生，高級工程師。研究方向?yàn)槔走_(dá)、艦炮系統(tǒng)作戰(zhàn)使用研究。

聲明

為適應(yīng)我國信息化建設(shè)，擴(kuò)大本刊以及作者的知識信息交流渠道，本刊已被目次頁上的數(shù)據(jù)庫全文收錄，其作者文章著作權(quán)使用費(fèi)與本刊稿酬一次性付給。如作者不同意文章被收錄，請?jiān)趤砀鍟r(shí)向本刊申明，本刊將作適當(dāng)處理。

《現(xiàn)代雷達(dá)》編輯部

Detection Method for Ground Slow Moving Targets by

Airborne SAR

YANG Yuelun

(The Military Delegate Office of Shipborne Gun System Navy at Shanghai,Shanghai 201913, China)

Abstract：In airborne SAR moving target detection, ground slow moving target's speed is small, so it's easy to be submerged by ground strong main lobe clutter and undetectable. In this paper, according to the characteristics of the movement of ground slow moving targets, a DPCA-CFAR joint detection method based on Frft is proposed. After using DPCA to achieve clutter suppression, CFAR processing is made to detect the moving targets. And according to the signal amplitude after DPCA and CFAR detection threshold, the minimum detectable velocity is inferred to show the performance of this proposed algorithm for slow moving targets detection. Further Frft is used to estimate the target speed. Finally the effectiveness of the algorithm is verified by simulation.

Key words：airborne synthetic aperture radar；ground slow moving target； displaced phase center antenna-constant false alarm rate；fractional Fourier transform

DOI：·顯示技術(shù)· 10.16592/ j.cnki.1004-7859.2015.12.019

收稿日期：2015-07-19

修訂日期：2015-09-18

通信作者：楊躍輪Email：10430412@qq.com

中圖分類號：TN957

文獻(xiàn)標(biāo)志碼：A

文章編號：1004-7859(2015)12-0083-05