星空背景的弱小目標(biāo)差分增強(qiáng)方法

魏 敏，文 武，張江楠

(成都信息工程大學(xué)計(jì)算機(jī)學(xué)院，四川成都610225)

0 引言

弱小目標(biāo)檢測(cè)預(yù)處理的主要目的是抑制背景和增強(qiáng)目標(biāo)信號(hào)，其方法繁多，常見的有:灰度拉伸［1］、直方圖變換［2］、高通濾波［3-4］目標(biāo)能量累積［5-7］、波原子變換［8］、四元數(shù)傅里葉變換［9］和形態(tài)濾波［10-11］等。圖像直接差分［12-13］或變形形式(采用基于時(shí)空域的雙邊濾波器估計(jì)背景，再將原圖與背景做差分運(yùn)算［14］)也廣泛應(yīng)用到弱小目標(biāo)檢測(cè)當(dāng)中，差分的主要目的是抑制背景，減小背景干擾，但是，對(duì)于信噪比較低的小目標(biāo)，經(jīng)過差分處理后有可能降低目標(biāo)的信噪比，因此，普通差分并不適合復(fù)雜背景的弱小目標(biāo)檢測(cè)。如果能夠找到一種差分方法，既能抑制背景，又能提高弱小目標(biāo)的信噪比，這對(duì)擴(kuò)大差分的應(yīng)用范圍和提高弱小目標(biāo)的檢測(cè)都是極其有意義的。

針對(duì)星空背景的弱小目標(biāo)檢測(cè)這一應(yīng)用場(chǎng)景，提出一種既能抑制背景星，又能增強(qiáng)弱小目標(biāo)的差分增強(qiáng)方法。星空背景的圖像可以分為兩大類，一類是背景星位置不變或緩動(dòng)(定義為做差分運(yùn)算的兩幀圖像中的背景星位移不大于1像素);另一類為背景星位置有明顯變化。這兩類圖像均可應(yīng)用提出的差分方法進(jìn)行背景星抑制(或完全去除)，增強(qiáng)弱小目標(biāo)，方法應(yīng)用到這兩類圖像時(shí)，對(duì)于背景星運(yùn)動(dòng)的情況需要先將背景星進(jìn)行配準(zhǔn)，再運(yùn)用該方法。為了算法描述簡(jiǎn)單，僅以第一類圖像的情形進(jìn)行描述。

1 算法描述

圖像差分運(yùn)算指的是將需要做差分運(yùn)算圖像的兩幀圖像(根據(jù)需要，這兩幀圖像可以不相鄰)進(jìn)行逐像素做差運(yùn)算，這種運(yùn)算可以表示為

其中，fn(i，j)表示第 n幀圖像(i，j)處的灰度，fn-1(i，j)與第n幀做差分運(yùn)算的那幀圖像(i，j)處的灰度。如前面所述，圖像差分的優(yōu)勢(shì)在于抑制背景而非增強(qiáng)弱小目標(biāo)，若直接將圖像差分運(yùn)用到弱小目標(biāo)檢測(cè)中，很難取得良好的效果，特別是當(dāng)目標(biāo)信噪比很低時(shí)。分析圖像差分運(yùn)算可以發(fā)現(xiàn)，差分時(shí)直接將對(duì)應(yīng)像素相減，并沒有考慮對(duì)應(yīng)兩個(gè)像素的灰度關(guān)系，提出的差分充分考慮對(duì)應(yīng)兩個(gè)像素的灰度值，并利用這兩個(gè)像素的關(guān)系，在差分時(shí)將他們的差值(差分結(jié)果)進(jìn)行放大，從而達(dá)到增強(qiáng)弱小目標(biāo)的效果。改進(jìn)后的差分表示為

其中，Sat-Value是圖像的飽和值，對(duì)于16 Bit的圖像，Sat-Value的值為65535。式(2)充分考慮了幀間灰度特征，同時(shí)將差分后的結(jié)果進(jìn)行了灰度拉伸，這種拉伸只受幀間對(duì)應(yīng)像素值的影響，換句話說，只要幀間對(duì)應(yīng)像素灰度值存在差異，這種拉伸方法就將這種差異放大，差異越大放大效果越明顯，這一特性對(duì)于弱小目標(biāo)差分增強(qiáng)是極其關(guān)鍵的，原因在于:差分時(shí)兩幀對(duì)應(yīng)的背景星灰度基本一致，經(jīng)過差分后背景星被抑制;差分時(shí)兩幀背景(沒有星的地方)灰度值也基本一致，經(jīng)過差分后也被抑制;而弱小目標(biāo)比背景灰度值略高，經(jīng)過差分后會(huì)拉大與周圍背景的差異，從而實(shí)現(xiàn)背景星抑制和弱小目標(biāo)增強(qiáng)。此外，文中差分增強(qiáng)方法與其他差分后再增強(qiáng)的方法相比較有還有兩個(gè)優(yōu)勢(shì):其一是增強(qiáng)效果不受周圍像素或全圖像素的統(tǒng)計(jì)值影響;其二是增強(qiáng)方法簡(jiǎn)單且運(yùn)算量小。

雖然經(jīng)過式(2)差分增強(qiáng)能夠很大程度抑制背景星的影響，但并不能完全將背景星剔除，通過實(shí)驗(yàn)分析發(fā)現(xiàn)，造成這一現(xiàn)象的原因在于:

(1)由于大氣擾動(dòng)或傳動(dòng)晃動(dòng)的影響，背景星的形狀在幀間會(huì)有細(xì)微的變化(位置或形狀)，此外，亮度飽和的背景星，其飽和溢出的拖尾［15］在幀間變化有時(shí)還比較大，這種變化通過改進(jìn)的差分處理后更加明顯，目標(biāo)形狀的這種變化如所示。

圖1 飽和星幀間形狀變化圖

從圖1可以看出，飽和背景星的形狀在幀間變化劇烈，這時(shí)，差分后必定會(huì)留下邊緣殘留。圖2是倒數(shù)第二與倒數(shù)第三幀圖像差分后殘留的效果圖。

圖2 差分邊緣殘留示意圖

(2)由于大氣擾動(dòng)，雖然背景星在幀間位置和形狀上幾乎沒有變化，但是背景星的灰度分布在幀間卻存在較大的變化，灰度的差異經(jīng)過差分及拉伸會(huì)造成背景星的殘留。這種幀間目標(biāo)位置不變，灰度分布存在變化的情況如圖3所示。

圖3 背景星幀間灰度變化示意圖

從圖3可以看出，背景星左上角的灰度有輕微變化，經(jīng)過差分增強(qiáng)后，這些微小變化被放大了10倍以上，這勢(shì)必導(dǎo)致背景星殘留。

為了克服以上兩種情況導(dǎo)致的背景星殘留，將式(2)的差分修改為

其中OSE(fn-1(i，j))是對(duì)上一幀圖像選擇合適的結(jié)構(gòu)元SE進(jìn)行膨脹，膨脹的作用是將上一幀背景星的形狀擴(kuò)大或?qū)⑸弦粠尘靶堑幕叶确植歼M(jìn)行調(diào)整，并取局部最大值，這個(gè)操作如所示。

使用式(3)進(jìn)行差分運(yùn)算，根據(jù)實(shí)際情況選取合適的結(jié)構(gòu)元能夠極大地消除差分后的背景星殘留。

對(duì)圖2中的兩幀圖像，采用式(3)的方法，其中結(jié)構(gòu)元SE大小為3×3，對(duì)前一幀圖像進(jìn)行膨脹后差分，其結(jié)果如圖4所示，從圖4可以看出，幀間變化的飽和星邊緣被有效去除。

圖4 差分前膨脹示意圖

2 實(shí)驗(yàn)結(jié)果

實(shí)驗(yàn)數(shù)據(jù)為某光電望遠(yuǎn)鏡采集的圖像序列，在圖像中背景星在幀間靜止，圖像中有多顆運(yùn)動(dòng)弱小目標(biāo)，從圖像中無法用肉眼區(qū)分出感興趣的目標(biāo)和背景星。為了更好地展示文中算法在弱小目標(biāo)檢測(cè)預(yù)處理的效果，將其中最弱的兩個(gè)目標(biāo)所在區(qū)域進(jìn)行裁剪，這兩個(gè)目標(biāo)從左上角運(yùn)動(dòng)到右下角，其中一個(gè)目標(biāo)的信噪比為4.4左右，對(duì)比度1.03，另一目標(biāo)信噪比為1.4左右，對(duì)比度1.008。信噪比和對(duì)比度的定義為

其中，mt是目標(biāo)的均值，mt是背景的均值，σb是背景的方差。圖5是序列圖像中具有代表性的兩幀圖像的一部分(全圖分辨率為2 K×2 K，為方便顯示和突出小目標(biāo)，因此只選取其中一部分顯示)，其中小方框中心是感興趣的目標(biāo)，這兩顆小目標(biāo)在形狀和亮度上與背景星相似，甚至更弱，采用普通的差分很難檢測(cè)出這些目標(biāo)。

圖5 序列圖像中的相鄰兩幀原始局部圖像

圖6是圖5采用文中算法的效果圖，圖6(a)為圖5的右邊圖像，矩形框內(nèi)為感興趣的小目標(biāo)，圖6(b)為普通差分的效果圖，從圖中可以看出，雖然大部分背景星被抑制，但是殘留了部分較亮背景星的邊緣，圖6(c)為文中差分公式(2)處理的結(jié)果，效果與普通差分類似，但目標(biāo)大小較普通差分大，背景噪聲也比普通差分大，這正好驗(yàn)證差分能增強(qiáng)目標(biāo)的特性(副作用是噪聲也增強(qiáng)，較弱的目標(biāo)幾乎被周圍的噪聲點(diǎn)淹沒)，圖6(d)是差分公式(3)處理的效果，背景星被完全去除，且噪聲點(diǎn)較圖6(c)明顯減少，噪聲點(diǎn)減少的原因在于上一幀圖像的隨機(jī)噪聲在膨脹時(shí)面積擴(kuò)大，當(dāng)幀隨機(jī)噪聲與上一幀圖像噪聲相減的機(jī)會(huì)增大，從而抑制本幀的噪聲點(diǎn)。

圖6 局部圖像及處理效果

圖7是序列圖像中兩個(gè)較弱目標(biāo)文中算法處理前后信噪比和對(duì)比度變化的對(duì)比圖，總體來說，通過文中算法處理后目標(biāo)的信噪比和對(duì)比度都得到極大提高，在序列圖像中，目標(biāo)一在第4、5幀穿過一顆背景星，目標(biāo)二在第21幀被噪聲淹沒，因此其信噪比和對(duì)比度為“0”，目標(biāo)二在第 5、16、20、21、23、29 幀差分后被背景淹沒，噪聲和對(duì)比度為“0”。

圖7 差分處理后目標(biāo)信噪比和對(duì)比度變化對(duì)比圖

經(jīng)過提出的差分增強(qiáng)處理后，目標(biāo)一的信噪比(平均值)由4.38提高到17.42，對(duì)比度(平均值)由1.026提高到88.888，目標(biāo)二的信噪比(平均值)由1.39提高到4.81，對(duì)比度(平均值)由1.008提高到23.882。

3 結(jié)束語

提出的改進(jìn)差分方法具有以下優(yōu)點(diǎn):

(1)運(yùn)算簡(jiǎn)單，滿足實(shí)時(shí)圖像處理需求。

(2)差分的時(shí)候同時(shí)增強(qiáng)弱小目標(biāo)，克服普通差分削弱目標(biāo)信號(hào)的缺點(diǎn)。

(3)差分前膨脹可以有效去除背景星邊緣，減少背景星邊緣對(duì)后續(xù)處理的干擾。

結(jié)合實(shí)驗(yàn)結(jié)果表明:文中方法對(duì)于提高目標(biāo)信噪比和改善對(duì)比度具有優(yōu)異的性能，是一種簡(jiǎn)單有效的弱小目標(biāo)檢測(cè)預(yù)處理方法。

［1］ 李曉冰.一種紅外測(cè)量圖像自適應(yīng)分段線性灰度變換方法［J］.光電子技術(shù)，2011，31(4):236-239.

［2］ 王炳健，劉上乾，拜麗萍.紅外圖像實(shí)時(shí)增強(qiáng)新算法［J］.光電工程，2006，33(1):46-49.

［3］ 候潔，辛云宏.基于高通濾波和圖像增強(qiáng)的紅外小目標(biāo)檢測(cè)方法［J］.紅外技術(shù)，2013，35(5):279-284.

［4］ L Yang，J Yang，K Yang.Adaptive detection for infrared small target under sea-sky complex background［J］.Electronics Letters，2004，40(17):1083-1085.

［5］ 葉斌，彭嘉雄.基于能力累計(jì)與順序形態(tài)濾波的紅外小目標(biāo)檢測(cè)［J］.中國(guó)圖象圖形學(xué)報(bào)，2002，7A(3):45-49.

［6］ 廖斌，楊衛(wèi)平，沈振康.基于多幀移位疊加的紅外小目標(biāo)檢測(cè)方法［J］.紅外與激光工程，2002，31(2):150-163.

［7］ 陳尚鋒，陳華明，盧煥章.基于加權(quán)動(dòng)態(tài)規(guī)劃和航跡關(guān)聯(lián)的小目標(biāo)檢測(cè)技術(shù)［J］.國(guó)防科技大學(xué)學(xué)報(bào)，2003，25(2):46-50.

［8］ 秦翰林，黃洋，姚柯柯，等.多尺度核局部歸一化的紅外圖像背景抑制［J］.強(qiáng)激光與離子束，2012，24(5):1063-1066.

［9］ Shenxiang Qi，Jie Ma，Hang Li，etc.Infrared small target enhancement via phase spectrum of Quaternion Fourier Transform［J］.Infrared Physics and Technology，2014，62:50-58.

［10］ XiangZhi Bai.Morphological operator for infrared dim small target enhancement using dilation and erosion through structuring element construction［J］.Optik.2013，Vol.124:6163-6166.

［11］ Jian-Nan Chi，Ping Fu，Dong-Shu Wang，.et al.A detection method of infrared image small target based on order morphology transformation and image entropy difference［J］.2005，8:5111-5116.

［12］ 陳朝陽(yáng)，張桂林.基于圖像對(duì)稱差分運(yùn)算的運(yùn)動(dòng)小目標(biāo)檢測(cè)方法［J］.華中理工大學(xué)學(xué)報(bào)，1998，26(9):34-38.

［13］ Jeff Houchard，Paul Kervin，John Africano.et al.Orbital debris detection program highlights from the Air Force Maui Optical Station［J］.SPIE Space Instrumentation and Dual-Use Technologies Vol.2214:7-20.

［14］ Tae Wuk Bae.Spatial and temporal bilateral filter for infrared small target enhancement［J］.Infrared Physics and Technology，2014，63:42-53.

［15］ 傅平.CCD圖像傳感器拖尾的研究［J］.壓電與聲光，2004，26(1):72-75.