低信噪比運動(dòng)紅外點(diǎn)目標的檢測

低信噪比運動(dòng)紅外點(diǎn)目標的檢測

為解決高空背景條件下紅外低信噪比運動(dòng)點(diǎn)目標的檢測問(wèn)題，本文提出了一種基于膨脹累加的檢測方法.運用形態(tài)膨脹運算能夠使同一目標處于不同幀上的能量仍然能夠實(shí)現有效的累加，從而達到目標增強的目的.同時(shí)本文還采用小波變換預處理的方法，對圖像中相關(guān)的1/f噪聲進(jìn)行白化.本文對這一方法進(jìn)行了大量的仿真實(shí)驗和詳細的理論分析.實(shí)驗證明該算法能夠快速穩定地檢測出信噪比為2的運動(dòng)點(diǎn)目標.
　　關(guān)鍵詞:紅外圖像處理;紅外目標檢測;小波變換

Moving Infrared Low SNR Target Detection algorithm

XIONG Hui,SHEN Zhen-kang,WEI Ji-bo,LI Ji-cheng
(Institute of Electronic Science and Engineering,National University of Defense Technology,Changsha 410073,China)

　　Abstract:This paper is presented to resolve the moving infrared target detection problem under the condition of low SNR.The scheme of multi-frame accumulation is often adopted in low SNR infrared target detection,but for moving low SNR target,this simple accumulation scheme maybe degrades greatly.In this paper,a detection algorithm based on morphological dilation before accumulation and tracking before detection is proposed.The dilation operation enables the energy of the same target in different frames can be accumulated effectively in spite of the movement of target,thus we can reach the aim of enhancing the target.A preprocessing algorithm using wavelet transform is also adopted to remove the correlation of 1/f noise in parallel scanned infrared image.The whole algorithm is well analyzed and simulated with infrared image sequence.A high performance is reached and the experimental result indicates that the algorithm can effectively detect the moving target with SNR=2.
　　Key words:infrared image processing；infrared target detection；wavelet transform；mathematical morphology

一、引　　言
　　在紅外尋的制導過(guò)程中，需要能夠盡快地截獲并鎖定跟蹤目標，由于紅外傳感器與目標之間的距離很遠（通常為幾十公里甚至上百公里），目標在像平面上的成像僅為點(diǎn)目標，并且淹沒(méi)在各種雜波背景中.通常的點(diǎn)目標的檢測方法采用多幀累加來(lái)提高信噪比，這種簡(jiǎn)單的累加方法在目標幀間移動(dòng)很小時(shí)可以有效地抑制噪聲并起到增強目標的作用.但對于運動(dòng)速度較大的目標(平均幾幀圖像移動(dòng)一個(gè)象像素)，這種簡(jiǎn)單的累加方法并不能夠實(shí)現有效的能量積累，很難檢測出低信噪比并處于快速運動(dòng)中的點(diǎn)目標.
　　為解決高空背景中低信噪比運動(dòng)點(diǎn)目標的檢測問(wèn)題，本文提出了一種基于膨脹累加、檢測前跟蹤的算法.將經(jīng)過(guò)白化預處理的圖像作形態(tài)膨脹運算之后再進(jìn)行累加，這樣仍然可以使同一點(diǎn)目標在不同幀上的能量累加到同一像素點(diǎn)上，起到突出目標的作用.邊檢測邊跟蹤利用目標候選點(diǎn)在連續幀上的信息剔除虛警假目標，采用航跡關(guān)聯(lián)判決的方法來(lái)提高算法的抗干擾性能.
　　本文采用紅外線(xiàn)列掃描圖像序列對算法進(jìn)行模擬. 紅外線(xiàn)列掃描圖像主要有兩種噪聲：圖像平面上的白噪聲和線(xiàn)列掃描方向上的非平穩1/f噪聲.本文采用小波變換的方法對線(xiàn)列掃描方向上的1/f噪聲進(jìn)行白化預處理.實(shí)驗證明，該算法獲得了很好的預處理效果.膨脹累加算法能夠快速檢測出信噪比為2的點(diǎn)目標.

二、低信噪比運動(dòng)點(diǎn)目標的檢測算法
　　1.圖像預處理
　　高空背景可近似看作均勻輻射背景，傳感器的內部噪聲將起主要作用.在線(xiàn)列紅外圖像中主要存在不相關(guān)的白噪聲和線(xiàn)列掃描方向上相關(guān)的1/f噪聲.可以合理地把噪聲看作加性噪聲，圖像序列可用下面的模型來(lái)表示：

x(i,j,t)=hδ(i-it,j-jk,t)+fn(i,j,t)+(i,j,t)　(1)

這里目標被看作位于點(diǎn)(it,jt)處幅度為h的沖激函數，圖像預處理就是對1/f噪聲fn(i,j)去相關(guān).實(shí)際上，1/f噪聲是一種統計自相似的分形隨機過(guò)程［1］.小波變換方法能夠刻劃分形過(guò)程的精細結構，對于非平穩過(guò)程的分析具有一定的優(yōu)越性.因此本文采用這一方法對噪聲進(jìn)行白化處理.
　　小波變換將信號在一個(gè)由小波函數ψ(t)的整數平移和尺度伸縮所構成的規范正交基上展開(kāi)，對于信號x(t),它能夠用小波正交基來(lái)表示，同時(shí)還可由其小波變換系數來(lái)恢復：

xmn=∫+∞-∞x(t)ψmn(t)dt　(2)
　(3)

{xmn}即為x(t)在平移伸縮系{ψmn(t)}上的小波變換，對于二進(jìn)小波變換，小波基函數可表示為：

ψmn(t)=2m/2ψ(2mt-n)　(4)

　　文獻［1］中提出了一種基于小波變換的1/f噪聲表示方法，并證明了1/f噪聲的小波變換系數可近似地看作彼此互不相關(guān)，通過(guò)各個(gè)尺度上方差的規一化，就能夠對1/f噪聲進(jìn)行白化.對于1/f噪聲信號x(t)，其小波變換系數的方差具有如下形式：

Var［xmn］=σ22-γm　(5)

其中γ為1/f噪聲的功率譜參數.由于小波變換系數近似互不相關(guān)，小波變換系數可以用下式表示：

xmn=［σ2-γm/2］υmn　(6)

其中小波變換系數υmn為零均值，單位方差并且互不相關(guān)的隨機變量，它對應著(zhù)一個(gè)平穩的白噪聲υ(t),如果用(σ2-γm/2)-1對xmn歸一化,就能夠通過(guò)重構方程(3)得到一個(gè)白化的隨機過(guò)程.
　　由于紅外線(xiàn)列掃描圖像只在掃描方向存在1/f噪聲，為減少計算量，只考慮掃描方向的白化，這樣可以只進(jìn)行一維的處理，考慮到加性的白噪聲，線(xiàn)列掃描序列的小波變換系數的方差為：

Var(xmi,j)=σ22-γm+σ2　(7)

其中σ2為白噪聲的方差，經(jīng)(σ22-γm+σ2)-1/2對小波變換系數xmi,j尺度化后，再用式(3)進(jìn)行信號重構，就能夠實(shí)現線(xiàn)列掃描方向的白化.通過(guò)上述的推導，可以得知預處理后的輸出為零均值，單位方差的白噪聲過(guò)程.
　　2.膨脹累加平均
　　當目標在幀間移動(dòng)很小時(shí)，直接采用多幀圖像累加的方法可以有效地抑制噪聲，提高信噪比.但是幀間運動(dòng)較大的點(diǎn)目標，直接累加并不能夠實(shí)現目標能量的有效積累.在本文中，首先對圖像進(jìn)行膨脹運算，將目標由一個(gè)點(diǎn)膨脹成一個(gè)斑點(diǎn)，從而使得不同幀上的目標能量仍然能夠實(shí)現累加，起到目標增強的作用.結構函數g(x)對函數f(x)的膨脹運算記為，它可以由下式計算：

　(8)

其中D(g)表示結構函數g(x)的定義域，進(jìn)一步可以得到二維膨脹運算：

　(9)

本文選擇平頂形的結構函數，即對于所有(x,y)∈D(g)，都有g(shù)(x,y)=1成立，此時(shí)某一個(gè)像素點(diǎn)的膨脹運算結果相當于以該點(diǎn)為原點(diǎn)的結構函數窗口內的所有像素灰度的最大值，則有：

　(10)

其中結構函數尺寸大小根據目標在幀間運動(dòng)的大小來(lái)選取.累加平均后的圖像可表示為：

　(11)

其中N為累加幀數，經(jīng)N幀累加平均后，目標點(diǎn)和噪聲點(diǎn)的均值保持不變，而均方差變?yōu)閱螏瑫r(shí)的1/N.最后將累加的結果取門(mén)限即得到一組目標候選點(diǎn).
　　3.數據關(guān)聯(lián)與檢測判決
　　數據關(guān)聯(lián)采用最近鄰關(guān)聯(lián)方法將目標候選觀(guān)測點(diǎn)連成多條航跡.把每一次累加輸出后的目標候選點(diǎn)視為一組，在得到T組連續候選點(diǎn)后，如果某一個(gè)目標航跡中的觀(guān)測點(diǎn)數不少于S個(gè)，判決該航跡為目標航跡.假設第i條航跡觀(guān)測數為nc(i)個(gè)，檢測判決即為如下的判決過(guò)程：

　(12)

三、模擬實(shí)驗結果
　　圖(1)給出了原始的紅外線(xiàn)列掃描圖像，從圖像中可以看到圖像在列方向(掃描方向)具有很強的相關(guān)性.圖(2)為經(jīng)小波變換預處理后的圖像，圖像列方向上的相關(guān)性得以濾除，圖(3)為圖像經(jīng)膨脹累加取平均的結果，從中可以看到目標得到增強，圖(6)為最終檢測判決結果.實(shí)驗證明該算法能夠有效地檢測出信噪比為2的運動(dòng)點(diǎn)目標.