樊香所，范錦龍，文良華，徐智勇

〈圖像處理與仿真〉

改進(jìn)時(shí)空濾波的紅外弱小目標(biāo)檢測

樊香所1,2，范錦龍3，文良華1，徐智勇4

（1. 宜賓學(xué)院智能制造學(xué)部，四川 宜賓 644600；2. 廣西科技大學(xué)，廣西土方機(jī)械協(xié)同創(chuàng)新中心，廣西 柳州 545006；3. 國家衛(wèi)星氣象中心，北京 100081；4. 中國科學(xué)院光電技術(shù)研究所，四川 成都 610209）

為了有效解決動(dòng)態(tài)背景變化導(dǎo)致弱小目標(biāo)檢測率低的問題，文中提出了改進(jìn)時(shí)空濾波的紅外弱小目標(biāo)檢測算法。首先在分析紅外圖像成像特性的基礎(chǔ)上，針對目標(biāo)區(qū)、背景區(qū)和邊緣輪廓區(qū)不同梯度特性的差異，提出改進(jìn)的各向異性空域?yàn)V波算法，該算法充分利用空間域的梯度信息來構(gòu)建不同方向的擴(kuò)散濾波函數(shù)，并結(jié)合圖像不同特性的梯度差異選取擴(kuò)散函數(shù)值最小的兩個(gè)方向的均值作為時(shí)域?yàn)V波結(jié)果，以最大限度地保留目標(biāo)信號；接著為有效增強(qiáng)弱小目標(biāo)的能量，針對高階累積量僅利用像元點(diǎn)時(shí)域信息來構(gòu)建能量增強(qiáng)的不足，提出了一種結(jié)合時(shí)空鄰域塊的能量增強(qiáng)算法，實(shí)驗(yàn)表明，本文提出的算法能有效提升動(dòng)態(tài)場景下的弱小目標(biāo)的檢測能力。

弱小目標(biāo)；改進(jìn)各向異性；能量增強(qiáng)；目標(biāo)檢測

0 引言

為了有效減少背景和動(dòng)態(tài)噪聲對弱小目標(biāo)檢測的影響，需要先期進(jìn)行預(yù)處理工作，以增強(qiáng)弱小信號的提取能力。當(dāng)前常用的預(yù)處理方式主要分為3種，背景建模法、圖像奇異點(diǎn)分析法和基于機(jī)器學(xué)習(xí)的估計(jì)法。背景建模法主要是采用背景估計(jì)算法來預(yù)測圖像的背景成分，緊接著將估計(jì)出的背景圖與原始圖像作差分運(yùn)算，經(jīng)過上述處理后獲取的差分圖像僅包含有目標(biāo)和少量噪點(diǎn)，達(dá)到去除大部分背景干擾的目的。這類算法包括二維最小均方（Two Dimensional Least Mean Squre，TDLMS）[1-2]濾波，改進(jìn)Top-Hat[3]和雙邊濾波[4]等。圖像奇異點(diǎn)分析法是當(dāng)目標(biāo)經(jīng)過圖像時(shí)會對圖像原有的背景特性造成“灰度擾動(dòng)”，通過尋找“灰度擾動(dòng)”所形成的“奇異性”來實(shí)現(xiàn)圖像預(yù)處理，如嚴(yán)高師等人通過建立圖像鄰域區(qū)域的奇異性表征函數(shù)來計(jì)算目標(biāo)與背景不同組分的奇異差異，實(shí)現(xiàn)目標(biāo)與背景的分離[5]；連可等人通過計(jì)算鄰域局部灰度極大值來尋找圖像中奇異點(diǎn)，并將灰度極大值與預(yù)設(shè)閾值進(jìn)行比較，達(dá)到檢測目標(biāo)的目的[6]。上述的背景建模法和奇異點(diǎn)分析法對于平穩(wěn)變化或是靜止的背景，檢測效果較好，但對于動(dòng)態(tài)變化的背景，由于目標(biāo)受外場環(huán)境的干擾，導(dǎo)致虛警率較高，檢測效果不好。為提升動(dòng)態(tài)場景下的目標(biāo)檢測能力，研究學(xué)者借助于機(jī)器學(xué)習(xí)的估計(jì)法從場景中獲取數(shù)據(jù)對模型進(jìn)行訓(xùn)練，這類方法主要利用一些數(shù)學(xué)模型將目標(biāo)檢測問題轉(zhuǎn)化為模式分類問題，主要包括主成分分析法[7]、混合高斯模型[8]、低秩表示法[9]、紅外圖像塊（Infrared patch image，IPI）模型[10]等。如Oliver等人提出的主成分分析法，先利用PCA來獲取訓(xùn)練數(shù)據(jù)的特征向量，然后將測試數(shù)據(jù)投影到特性向量中獲取背景成分，實(shí)現(xiàn)最終目標(biāo)的檢測，當(dāng)場景隨著光照的變化引入較多動(dòng)態(tài)成分時(shí)，檢測效果不佳；Bouwmans等人提出混合高斯模型，該算法需從場景數(shù)據(jù)中獲取訓(xùn)練樣本來估計(jì)模型的參數(shù)，當(dāng)受到噪聲干擾或光照變化快時(shí)，模型參數(shù)難以適應(yīng)動(dòng)態(tài)背景估計(jì)；低秩表示法和紅外圖像塊模型需要獲取場景數(shù)據(jù)進(jìn)行參數(shù)訓(xùn)練，當(dāng)場景發(fā)生變化時(shí)，預(yù)先訓(xùn)練的參數(shù)難以應(yīng)對動(dòng)態(tài)場景的變化，使得檢測結(jié)果中具有較多的噪點(diǎn)，加上遠(yuǎn)距離成像，目標(biāo)被各種噪聲雜波淹沒，圖像中目標(biāo)的形狀和紋理特征并不明顯，難以形成有效的特征來支撐對模型的訓(xùn)練，這增大了構(gòu)建訓(xùn)練模型的難度。

針對上述檢測方法存在的問題，文中提出了改進(jìn)時(shí)空濾波的弱小目標(biāo)檢測方法。首先考慮到各向異性濾波方法在面臨動(dòng)態(tài)場景具有良好地檢測結(jié)果[11]，將其引入文中，并對其進(jìn)行改進(jìn)，文獻(xiàn)[12]僅利用空域范圍內(nèi)的4個(gè)不同方向的梯度差來構(gòu)建擴(kuò)散函數(shù)，并求取4個(gè)方向的均值作為差分結(jié)果，這樣處理導(dǎo)致擴(kuò)散函數(shù)難以區(qū)分邊緣輪廓區(qū)與目標(biāo)的差異，導(dǎo)致差分圖中難以去除這部分邊緣輪廓區(qū)，文中深入分析圖像各個(gè)組分在空域4個(gè)方向的梯度差異，通過構(gòu)建改進(jìn)各向異性差分濾波來獲取差分濾波結(jié)果，這樣將極大程度地保留圖像中的目標(biāo)信號；其次，經(jīng)過上述差分濾波后，僅有效保留目標(biāo)信號，為進(jìn)一步增強(qiáng)目標(biāo)信號，提出了一種結(jié)合時(shí)空鄰域塊的能量增強(qiáng)算法，該算法采用像元塊的鄰域信息來建立一種新的結(jié)合時(shí)空信息的能量評價(jià)模型，以彌補(bǔ)高階累積量僅利用單個(gè)像元點(diǎn)進(jìn)行能量增強(qiáng)的不足，提升了弱小信號的檢測能力。

1 改進(jìn)各向異性空間濾波

圖像一般由3部分構(gòu)成，即背景+目標(biāo)+非平穩(wěn)邊緣輪廓區(qū)（噪聲）。從圖1中分析發(fā)現(xiàn)，圖像中不同組分在不同方向的梯度差異較明顯，背景成分的中心塊與4個(gè)臨近方向的梯度差異較小，目標(biāo)成分中4個(gè)方向的梯度較大，邊緣輪廓部分中處于平穩(wěn)區(qū)域到非平穩(wěn)區(qū)域的過渡地帶的方向梯度差異較大，其余方向梯度差異較小。因此可結(jié)合圖像不同組分的差異，對圖像進(jìn)行差別化處理，可有效保留目標(biāo)成分。

圖1 圖像中不同組分的梯度差異

為了獲取圖像中的目標(biāo)信號，利用文獻(xiàn)[12]提出的各向異性擴(kuò)散函數(shù)來進(jìn)行濾波處理，具體公式如下：

式中：?為圖像中不同像元位置的梯度；為常數(shù)，一般取100；1(×)和2(×)為兩個(gè)不同的擴(kuò)散函數(shù)。

文中通過結(jié)合擴(kuò)散函數(shù)來構(gòu)建目標(biāo)在4個(gè)不同方向的梯度差，實(shí)現(xiàn)差分濾波，具體表達(dá)式如下：

其中，為輸入圖像；DU，DD，DL和DR是指以(,)為中心的上、下、左、右4個(gè)不同像元的梯度差；是像元移動(dòng)的步長。根據(jù)4個(gè)不同方向的梯度差來建立各向異性濾波算法，具體如下：

式中：(,)表示像元點(diǎn)的坐標(biāo)位置；￠(,)為濾波后的結(jié)果。

分析發(fā)現(xiàn)，文獻(xiàn)[12]的各向異性擴(kuò)散函數(shù)僅采用單個(gè)像元的信息來構(gòu)建4個(gè)方向的梯度差異，將導(dǎo)致空域范圍內(nèi)的信息利用不夠充分，難以區(qū)別目標(biāo)點(diǎn)和噪聲點(diǎn)；同時(shí)僅簡單地利用4個(gè)不同方向的均值作為濾波后的結(jié)果，當(dāng)圖像位于非平穩(wěn)邊緣輪廓區(qū)域時(shí)，4個(gè)不同方向中至少有兩個(gè)方向的擴(kuò)散系數(shù)較大，如僅簡單地計(jì)算4個(gè)方向的均值，致使位于目標(biāo)區(qū)域中的像元和非平穩(wěn)輪廓區(qū)域中的像元擴(kuò)散結(jié)果值差別不明顯，差分濾波后將難以去除非平穩(wěn)邊緣輪廓區(qū)，導(dǎo)致差分圖中留存有大片的邊緣噪點(diǎn)，對后續(xù)分割提取目標(biāo)點(diǎn)造成很大的干擾。為了減少虛假噪點(diǎn)對弱小目標(biāo)檢測的影響，本文針對文獻(xiàn)[12]的各向異性濾波算法的不足，提出了相應(yīng)的改進(jìn)思路：首先采用空域范圍內(nèi)的梯度塊來構(gòu)建4個(gè)方向的梯度差，進(jìn)而依據(jù)圖像各個(gè)組分的梯度差異來建立4個(gè)方向的擴(kuò)散函數(shù)，并通過構(gòu)建4個(gè)方向的最小優(yōu)化函數(shù)來優(yōu)選擴(kuò)散函數(shù)值最小的兩個(gè)方向作均值處理后，獲取最終的濾波結(jié)果，上述處理將很好地去除圖像的背景和非平穩(wěn)邊緣輪廓區(qū)，極大限度地保留目標(biāo)成分。目標(biāo)在4個(gè)不同方向的梯度差，如下式：

式中：DU，DD，DL和DR是指以中心塊為中心的上、下、左和右4個(gè)方向圖像塊的梯度差；×為空域的范圍大??；和為空域×的坐標(biāo)位置；fix()為鄰近取整函數(shù)。

根據(jù)上述分析，改進(jìn)各向異性差分濾波表達(dá)式如下：

式中：min1和min2分別是4個(gè)方向中擴(kuò)散函數(shù)值最小的兩個(gè)參數(shù)，采用這兩個(gè)參數(shù)的均值對圖像進(jìn)行差分濾波，可有效去除圖像的背景區(qū)和邊緣輪廓區(qū)。

2 結(jié)合時(shí)空鄰域塊的能量增強(qiáng)算法

文獻(xiàn)[13]指出，序列圖像經(jīng)過差分濾波后，從時(shí)間域上來分析經(jīng)過差分處理后獲取的序列差分圖可看作是零均值的高斯過程，這樣可借助于高階累積量理論將差分圖中檢測弱小信號的問題變?yōu)閺母咚乖肼曋蟹纸獬龇歉咚剐盘柍煞郑詫?shí)現(xiàn)目標(biāo)信號與噪聲的分離[1]。因此，可借助于高階累積量來實(shí)現(xiàn)弱信號的檢測。文中依據(jù)高階累積量理論對差分圖構(gòu)建相應(yīng)的目標(biāo)與噪聲二元假設(shè)函數(shù)，具體如下:

式中：￠(,,)為第幀差分濾波后獲取的圖像；(,,)是高斯噪聲，其中(,,)與(,,)相互獨(dú)立。具體的3階累積量表達(dá)式如下：

由于(,,)為高斯噪聲，依據(jù)高階累積量理論，零均值的高階累積量（3階以上）恒為零，固有3N(1,2)等于0，因此可對上式進(jìn)行簡化，具體如下：

在實(shí)際應(yīng)用中，為有效實(shí)現(xiàn)目標(biāo)的檢測，定義了3階累積量的判斷準(zhǔn)則，具體公式為：

式中：為累加的幀數(shù)；3F()為第幀三階累積后的能量大??；為閾值。

文獻(xiàn)[13]的高階累積量僅利用像元點(diǎn)(,)在時(shí)間域上的信息來構(gòu)建判斷準(zhǔn)則，由于只依靠了像元的時(shí)域信息，這將極大地削弱弱小信號的增強(qiáng)效果，因此文中構(gòu)建了結(jié)合時(shí)空鄰域塊的增強(qiáng)能量模型，具體表達(dá)式如下：

式中：￠()為融合空域信息后獲取的第幀能量圖像；￠3F()為第幀結(jié)合時(shí)空鄰域塊信息的能量大小。

3 結(jié)果與分析

在結(jié)果分析中，為了對比不同算法的差分濾波效果，文中采用目標(biāo)平均灰度值（average grayscale value，AGV）與局部信噪比（local signal noise ratio，LSNR）[1]來對差分濾波效果進(jìn)行評價(jià)，各個(gè)指標(biāo)如式(11)：

式中：(o, o)是中心像元所在的坐標(biāo)位置；(o, o)為目標(biāo)塊組成的坐標(biāo)集合，其像元總個(gè)數(shù)為；(,)為像元(,)的灰度值；T為目標(biāo)區(qū)域的均值；B為背景區(qū)域的均值；B為背景區(qū)域的標(biāo)準(zhǔn)差。

3.1 濾波參數(shù)分析

對改進(jìn)各項(xiàng)異性差分濾波結(jié)果有關(guān)聯(lián)的指標(biāo)為梯度步長，文中取擴(kuò)散函數(shù)中的常數(shù)＝100，取公式(4)步長＝(1, 5, 10, 15)時(shí)，分別獲取對應(yīng)的差分圖像，并利用公式(11)的LSNR來計(jì)算差分圖像中目標(biāo)的局部信噪比，由此畫出了梯度步長與局部信噪比LSNR間的走勢圖，如圖2所示。依據(jù)圖2可知隨著梯度步長的增大，LSNR隨之變大，當(dāng)＝5時(shí)，LSNR達(dá)到2.35，之后LSNR開始下降，從中可知兩個(gè)圖像塊的梯度跟著梯度步長的增大而越大，使得改進(jìn)后的擴(kuò)散系數(shù)值越大，越有利于保留目標(biāo)信號，但隨著梯度步長的增大，到達(dá)一定限度后，圖像塊由目標(biāo)區(qū)跨入邊緣輪廓區(qū)，此時(shí)目標(biāo)區(qū)和邊緣輪廓區(qū)的梯度差異不大，擴(kuò)散函數(shù)值較小，難以去除邊緣輪廓噪聲，導(dǎo)致差分圖中噪聲較多，目標(biāo)局部信噪比較低。通過上述分析，實(shí)驗(yàn)中?。?。

圖2 步長k與LSNR間的關(guān)系

3.2 濾波結(jié)果分析

文中選用3個(gè)場景的序列圖像進(jìn)行實(shí)驗(yàn)，其中場景1為目標(biāo)點(diǎn)在動(dòng)態(tài)變化的背景中運(yùn)動(dòng)，目標(biāo)點(diǎn)由中部斜向上運(yùn)行，再拐彎往下運(yùn)動(dòng)；場景2中為目標(biāo)從右下角向左上角作直線運(yùn)動(dòng)；場景3為兩個(gè)目標(biāo)點(diǎn)分別作斜對角直線運(yùn)動(dòng)。具體的場景信息如表1所示。接著將本文提出的改進(jìn)各項(xiàng)異性差分濾波算法與文獻(xiàn)[3]的改進(jìn)Top-hat、文獻(xiàn)[9]的LRR模型、文獻(xiàn)[10]的IPI分解法和文獻(xiàn)[12]的各向異性濾波進(jìn)行比較，具體的結(jié)果見表2。

表1 序列圖像信息

從表2的兩個(gè)指標(biāo)可以看出，對于不同場景的圖像，本文提出的差分濾波算法比其他的算法在保留目標(biāo)信號AGV和LSNR方面較高。分析發(fā)現(xiàn)，AF（Anisotropic Filtering）算法由于僅是簡單地求取4個(gè)方向的擴(kuò)散函數(shù)作為差分濾波，難以區(qū)分目標(biāo)與邊緣輪廓的差異，導(dǎo)致差分圖含有較多的邊緣輪廓區(qū)，AGV和LSNR相較其他算法低；ITH（Improved Top Hat）算法利用不同尺度的結(jié)構(gòu)元素，能很好地消除大部分背景對提取目標(biāo)信號的影響，但其對于動(dòng)態(tài)變化的背景，由于其結(jié)構(gòu)元易受動(dòng)態(tài)背景的干擾，導(dǎo)致差分圖中含有較多的噪點(diǎn)；當(dāng)背景中含有較多的邊緣輪廓區(qū)時(shí)，這些強(qiáng)邊緣輪廓區(qū)會破壞整個(gè)背景圖像數(shù)據(jù)的低秩特性，其同樣滿足“稀疏”特征，導(dǎo)致LRR（Low Rank Representation）易將這些區(qū)域誤判為虛警目標(biāo)，降低目標(biāo)的局部信噪比；IPI算法采用分塊的思想劃分圖像，其目的是進(jìn)一步約束圖像塊的低秩特性，但當(dāng)圖像塊處于邊緣輪廓區(qū)域時(shí)，圖像塊非零的奇異值個(gè)數(shù)較多，難以滿足低秩特征，因此獲取的差分圖虛警率較高；本文提出的算法能有效抑制背景和邊緣輪廓區(qū)，差分圖中有效地保留了目標(biāo)信號。

表2 不同算法獲取的背景建模結(jié)果

3.3 增強(qiáng)結(jié)果分析

為了分析本文提出算法的增強(qiáng)結(jié)果，分別采用高階累積量（Higher order cumulant，HOC）增強(qiáng)方法與本文提出的增強(qiáng)方法來對差分圖像進(jìn)行增強(qiáng)處理，并分別計(jì)算兩種算法獲取的AGV和LSNR，如表3所示。本文提出的算法充分利用了像元的時(shí)間域和空間域的信息來增強(qiáng)目標(biāo)信號，比高階累積量獲得更好的目標(biāo)增強(qiáng)結(jié)果，AGV和LSNR分別達(dá)到225和12.67dB。

3.4 檢測結(jié)果分析

文中先采用改進(jìn)的各項(xiàng)異性來獲取差分圖像，接著利用結(jié)合時(shí)空鄰域塊的增強(qiáng)能量模型對差分圖像進(jìn)行增強(qiáng)處理，最后采用雙窗分割算法提取目標(biāo)信號。為了本文算法的檢測效果，選用3個(gè)處于動(dòng)態(tài)變化的真實(shí)場景進(jìn)行實(shí)驗(yàn)，并將本文算法與ITH、LRR、IPI和AF算法進(jìn)行比較，3個(gè)場景不同算法的檢測結(jié)果如圖3、圖4和圖5所示。

表3 兩種算法的增強(qiáng)結(jié)果

圖3 場景1不同算法的檢測結(jié)果

圖4 場景2不同算法的檢測結(jié)果

從3個(gè)場景看出，ITH算法易受結(jié)構(gòu)元的影響，當(dāng)選取結(jié)構(gòu)元小于孤立噪點(diǎn)的尺寸時(shí)，導(dǎo)致差分圖中含有較多的噪點(diǎn)，見圖3(b)、4(b)和5(b)；LRR算法易受動(dòng)態(tài)非平穩(wěn)背景的影響，這些非平穩(wěn)背景由于是動(dòng)態(tài)變化的，LRR易將這些區(qū)域誤判為虛警目標(biāo)，見圖3(c)、4(c)和5(c)；IPI算法將圖像劃分為不同的塊，這樣圖像子塊在其鄰域范圍內(nèi)具有較強(qiáng)的相關(guān)性，有效地保留圖像子塊的低秩特性，但當(dāng)圖像塊處于邊緣輪廓區(qū)域時(shí)，圖像子塊的相關(guān)度降低，因此差分圖難以去除邊緣輪廓區(qū)的噪聲，見圖3(d)、4(d)和5(d)；各向異性濾波由于僅是簡單地求取4個(gè)方向的擴(kuò)散函數(shù)的均值作為差分結(jié)果，導(dǎo)致差分結(jié)果圖中難以區(qū)分邊緣輪廓區(qū)和目標(biāo)區(qū)域的差異，檢測結(jié)果中留有較多的噪聲，見圖3(e)、4(e)和5(e)；本文算法先采用改進(jìn)的各向異性濾波去除大部分非平穩(wěn)背景的干擾，緊接著采用結(jié)合時(shí)空鄰域塊的能量增強(qiáng)算法來增強(qiáng)目標(biāo)信號，結(jié)果圖中僅含有少量的噪聲，見圖3(f)、4(f)和5(f)。

為進(jìn)一步分析這些算法的性能，繪制了3個(gè)場景的ROC曲線，其中橫軸為檢測率d，縱軸為虛警率f。3個(gè)場景的ROC曲線如圖6所示。

圖5 場景3不同算法的檢測結(jié)果

圖6 3個(gè)場景ROC曲線

從圖6三個(gè)場景的ROC曲線可看出，本文算法獲取的檢測率最好，其次是IPI和LRR，AF算法最低。分析發(fā)現(xiàn)，本文算法先利用改進(jìn)的各項(xiàng)異性差分濾波算法來構(gòu)建4個(gè)方向梯度塊的擴(kuò)散函數(shù)，依據(jù)擴(kuò)散函數(shù)在4個(gè)方向梯度塊的差異對圖像進(jìn)行差異化處理，通過構(gòu)建4個(gè)方向的最小優(yōu)化函數(shù)來優(yōu)選擴(kuò)散函數(shù)值最小的兩個(gè)方向作均值處理后，獲取最終的濾波結(jié)果，接著采用提出能量增強(qiáng)算法來充分利用像元塊的時(shí)空域信息來構(gòu)建能量的判別模型，去除圖像噪聲的同時(shí)有效檢測目標(biāo)信號，如在虛警率f＝0.032時(shí)，本文算法獲取的檢測率d大于90%，而其他算法均小于88%，見圖6(a)和6(b)；而對于動(dòng)態(tài)且具有較多非平穩(wěn)邊緣輪廓的場景，由于邊緣輪廓區(qū)會導(dǎo)致圖像各個(gè)組分間的相關(guān)度降低，破壞了整個(gè)背景圖像低秩特性，導(dǎo)致LRR算法檢測率降低；IPI雖采用分塊的思想，進(jìn)一步約束鄰域塊的低秩特性，但當(dāng)圖像塊處于邊緣輪廓區(qū)域時(shí)，圖像塊非零的奇異值個(gè)數(shù)較多，獲取的差分圖虛警率較高；傳統(tǒng)的ITH和AF也同樣易受孤立噪點(diǎn)和邊緣輪廓區(qū)的影響，這兩類算法對于動(dòng)態(tài)場景中，相比于其他算法，檢測率較低，如在P＝0.04時(shí)，它們的檢測率d都小于80%，而其他算法的檢測率d都大于80%，見圖6(c)。

4 總結(jié)

文中在深入分析圖像各個(gè)組分在4個(gè)方向的梯度差異基礎(chǔ)上，構(gòu)建了利用4個(gè)方向梯度塊的擴(kuò)散濾波函數(shù)，同時(shí)通過構(gòu)建4個(gè)方向的最小優(yōu)化函數(shù)來優(yōu)選擴(kuò)散函數(shù)值最小的兩個(gè)方向作均值處理后，作為最終的濾波結(jié)果，這樣將極大程度地保留圖像中的目標(biāo)信號；進(jìn)而為進(jìn)一步增強(qiáng)目標(biāo)信號，提出了一種結(jié)合時(shí)空鄰域塊的能量增強(qiáng)算法，該算法采用像元塊的鄰域信息來建立一種新的結(jié)合時(shí)空信息的能量判別模型，有效增強(qiáng)了弱小信號的能量。通過3個(gè)動(dòng)態(tài)場景的實(shí)驗(yàn)表明，本文算法很好地剔除了大部分噪點(diǎn)干擾的同時(shí)有效保留了目標(biāo)信號。

[1] 樊香所. 序列圖像弱小目標(biāo)檢測與跟蹤算法研究[D]. 成都: 電子科技大學(xué), 2019.

FAN Xiangsuo. Dim and Small Targets Detection and Tracking Algorithms in Sequence Image[D]. Chengdu: University of Electronic Science and Technology of China, 2019.

[2] Bae T W, Kim Y C, Ahn S H, et al. An efficient two dimensional least mean square based on block statistics for small target detection[J]., 2009, 30(10): 1092-1101.

[3] BAI X Z, ZHOU F G, JIN T. Enhancement of dim small target through modified top-hat transformation under the condition of heavy clutter[J]., 2010, 90(1): 1643-1654.

[4] 秦翰林, 周慧鑫, 劉上乾, 等. 基于雙邊濾波的弱小目標(biāo)背景抑制[J]. 強(qiáng)激光與粒子束, 2009, 21(1): 25-28.

QIN H, ZHOU H, LIU S, et al. Dim and small target background suppression using bilateral filtering[J]. High, 2009, 21(1): 25-28.

[5] 嚴(yán)高師, 畢務(wù)忠. 基于區(qū)域奇異性濾波的小目標(biāo)檢測[J]. 光學(xué)技術(shù), 2006, 33(2): 163-165, 169.

YAN G S, BI W Z. Detection algorithm of small target based on regional singularity filter[J]., 2006, 33(2): 163-165, 169.

[6] 連可, 王厚軍, 李丹. 基于紅外目標(biāo)局部灰度特性分析的管道濾波方法[J]. 彈箭與制導(dǎo)學(xué)報(bào), 2011, 31(4): 200-206.

LIAN K, WANG H J, LI D. Pipeline filtering method based on feature analysis of local grey level of small infrared target[J]., 2011, 31(4): 200-203, 206.

[7] Oliver N M, Rosario B, Pentland A P. A Bayesian computer vision system for modeling human interactions[J]., 2000, 22(8): 831-843.

[8] Bouwmans T, Baf F E, Vachon B. Background modeling using mixture of Gaussians for foreground detection - a survey [J]., 2008, 1(3): 219-237.

[9] 何玉杰, 李敏, 張金利, 等. 基于低秩三分解的紅外圖像雜波抑制[J]. 光學(xué)與精密工程, 2015, 23(7): 2069-2078.

HE Y J, LI M, ZHANG J L, et al. Clutter suppression of infrared image based on three-component low rank matrix decomposition[J]., 2015, 23(7): 2069-2078

[10] GUO J, WU Y, DAI Y. Small target detection based on reweighted infrared patch-image model[J]., 2018, 12(1): 70-79.

[11] 陸福星, 李夜金, 陳忻, 等. 基于Top-Hat變換的PM模型弱小目標(biāo)檢測[J]. 系統(tǒng)工程與電子技術(shù), 2018, 40(7): 1417-1422.

LU F X, LI Y J, CHEN X, et al. Weak target detection for PM model based on Top-hat transform[J]., 2018, 40(7): 1417-1422.

[12] 周慧鑫, 趙營, 秦翰林, 等. 多尺度各向異性擴(kuò)散方程的紅外弱小目標(biāo)檢測算法[J].光子學(xué)報(bào), 2015, 44(9): 146-150.

ZHOU H X, ZHAO Y, QIN H L, et al. Infrared dim and small target detection algorithm based on multi-scale anisotropic diffusion equation[J]., 2015, 44(9): 146-150.

[13] 唐意東, 黃樹彩, 鐘宇, 等. 基于形態(tài)學(xué)和高階統(tǒng)計(jì)量的弱小運(yùn)動(dòng)目標(biāo)檢測[J].現(xiàn)代防御技術(shù), 2016, 44(2): 151-156.

TANG Y D, HUANG S C, ZHONG Y, et al. Moving dim target detection based on morphology and high-order statistics in infrared image[J]., 2016, 44(2): 151-156.

Infrared Dim-Small Target Detection Based on Improved Spatio-Temporal Filtering

FAN Xiangsuo1,2，F(xiàn)AN Jinlong3，WEN Lianghua2，XU Zhiyong4

(1.,,644600,; 2.,,545006,; 3.,100081,; 4.,,610209,)

To effectively solve the problem of low detection rates of dim and small targets caused by dynamic background changes, a detection method based on spatio-temporal filtering is proposed in this paper. Based on an analysis of the imaging characteristics of infrared images, an improved anisotropic spatial filtering algorithm is proposed to evaluate the difference in various gradient characteristics of the target area, background area, and edge contour area. The algorithm fully utilizes the gradient information in the spatial domain to construct the diffusion filter function in different directions. According to the gradient difference in various characteristics of the image, the mean value of the two directions with the smallest value of the diffusion function is selected as the result of spatial filtering to retain the target signal to the maximum extent. To effectively enhance the energy of dim and small targets and address the shortcomings of high-order cumulants that only use the temporal domain information of pixel points for energy enhancement, an energy enhancement algorithm based on spatial-temporal neighborhood blocks is proposed. Experimental results reveal that the proposed algorithm can effectively enhance the detection of dim and small targets in dynamically changing scenes.

dim and small target, improve anisotropic filtering, energy enhancement, targets detection

TP751

A

1001-8891(2022)05-0475-08

2021-01-11；

2021-06-28.

樊香所（1987-），男，博士，研究方向?yàn)橛?jì)算機(jī)視覺與圖像處理，E-mail: 100002085@gxust.edu.cn。

范錦龍（1975-），男，副研究員，主要從事遙感圖像處理，E-mail: fanjl@ cma.gov.cn。

國家自然科學(xué)基金項(xiàng)目（62001129）；廣西自然科學(xué)基金（2021GXNSFBA075029）；國家自然科學(xué)基金項(xiàng)目（61975171）。