基于各向異性雙邊濾波紅外背景抑制方法研究

王濤 陳凡勝 蘇曉鋒

摘 要：紅外復(fù)雜背景抑制作為紅外弱小目標(biāo)檢測(cè)的重要環(huán)節(jié)，直接影響著后續(xù)算法的目標(biāo)檢測(cè)概率和虛警率.針對(duì)各向同性雙邊濾波背景抑制方法因無法改變背景預(yù)測(cè)方向，在方向性的結(jié)構(gòu)化云邊緣存在較大預(yù)測(cè)殘差的問題，本文設(shè)計(jì)了一種利用原始圖像鄰域統(tǒng)計(jì)信息自適應(yīng)更新濾波尺度、掩模尺寸和濾波方向的各向異性雙邊濾波背景預(yù)測(cè)器，以最大限度地分離目標(biāo)和背景.仿真實(shí)驗(yàn)表明該方法對(duì)輸入信雜比為4的目標(biāo)的信雜比增益大于4和對(duì)原始圖像鄰域雜波的背景抑制算子大于4.

關(guān)鍵詞：弱小目標(biāo)檢測(cè)；雙邊濾波；自適應(yīng)濾波；背景抑制；各向異性濾波

中圖分類號(hào)：TP391 文獻(xiàn)標(biāo)志碼：A

Abstract：Infrared complex background suppression is an important part of infrared dim target detection， which directly affects the target detection probability and false alarm rate of subsequent algorithms. As the isotropic bilateral filtering background suppression method can not change the background prediction direction and remain a large residual in the direction of the structured cloud edge， in this paper， an anisotropic bilateral filter which adapts the filter parameters including filter scale， mask size and filter direction by using the local characteristics of the image was designed. It separates the target from background more precisely. The test results show that the proposed method has a signal to cluter ratio gain for target and a background suppression operator for clutter， and both are up to 4 with the input signal to clutter ratio of 4.

Key words：dim target detection； bilateral filter；adaptive filter； background suppression； anisotropic filter

紅外目標(biāo)檢測(cè)是紅外告警系統(tǒng)的重要研究?jī)?nèi)容，當(dāng)目標(biāo)成像距離較遠(yuǎn)時(shí)，目標(biāo)在探測(cè)器上的成像尺寸較小，極易淹沒在復(fù)雜背景之中，使得其檢測(cè)與追蹤變得非常困難.紅外背景抑制利用圖像的空域和時(shí)域相關(guān)性，以適當(dāng)?shù)念A(yù)測(cè)模型估計(jì)背景，增強(qiáng)目標(biāo)，提高紅外告警系統(tǒng)的目標(biāo)檢測(cè)性能.

傳統(tǒng)的紅外背景抑制方法有空域背景抑制方法和時(shí)域背景抑制方法.空域背景抑制方法基于圖像的空間鄰域相關(guān)性，容易在圖像邊緣等弱相關(guān)性區(qū)域殘留大量的背景，主要有各向異性擴(kuò)散[1]、小波變換[2]、數(shù)學(xué)形態(tài)學(xué)[3]、自適應(yīng)局部統(tǒng)計(jì)信息[4]、多尺度分割[5-10]等.以時(shí)域剖面[11]、時(shí)空域最小均方誤差[12]為代表的時(shí)域背景抑制方法雖然背景抑制效果要優(yōu)于空域法，但算法性能受制于衛(wèi)星平臺(tái)的穩(wěn)定性，當(dāng)衛(wèi)星平臺(tái)抖動(dòng)量較大時(shí)，幀間圖像的相關(guān)性降低，時(shí)域背景抑制方法難以應(yīng)用.

雙邊濾波是一種兼顧圖像空間鄰近度和灰度相似性的非線性濾波方法，相對(duì)于傳統(tǒng)的邊緣抑制方法，具有更好的邊緣抑制性能.李凡等人[13]將雙邊濾波運(yùn)用在紅外圖像背景抑制之中，取得了信噪比增益大于2的良好效果，曹瑛等人[14]采用模糊分類方法自適應(yīng)的選取像素雙邊濾波鄰域，降低了算法的復(fù)雜度.然而，傳統(tǒng)的雙邊濾波方法并沒有考慮到圖像邊緣的方向性，本文提出的各向異性雙邊濾波背景抑制方法，使濾波方向平行于邊緣方向，并通過圖像的局部統(tǒng)計(jì)信息自適應(yīng)地調(diào)節(jié)濾波尺度，仿真實(shí)驗(yàn)驗(yàn)證了各向異性雙邊濾波背景抑制算法的有效性，背景抑制后的信雜比增益可以達(dá)到4.

1 各向異性雙邊濾波

1.1 雙邊濾波

雙邊濾波最早可以追溯到Aurich和Weule[15]的非線性高斯濾波研究中，也曾出現(xiàn)在Smith和Brady[16]的SUSAN框架中，最終由Tomasi以及Manduchi[17]正式提出.其定義可由式（1）概括

1.2 各向異性雙邊濾波

雙邊濾波雖然能很好地保持圖像的邊緣，但傳統(tǒng)的雙邊濾波并沒有考慮到圖像邊緣的方向性信息，最差的情況是當(dāng)濾波方向與邊緣方向垂直時(shí)，會(huì)存在較大的邊緣背景殘留，Geusebroek等人[18]提出的各向異性快速高斯濾波器將邊緣方向投影到高斯濾波方向，使濾波器方向平行于邊緣方向.

1.3 濾波尺度

各向異性雙邊濾波中，不同的濾波尺度對(duì)于濾波效果的影響依然很大，Kugu[19]采用MSE（Mean Square Error）、Entropy、SSIM（Structural Similarity Based Image Quality Assessment）、Second Derivative度量不同濾波尺度對(duì)雙邊濾波遙感圖像處理的影響，Peng和Rao[20]采用信噪比最大化準(zhǔn)則推導(dǎo)出最優(yōu)的濾波尺度為：σr=1.5σi，σs=2.5σi，σi是圖像的噪聲標(biāo)準(zhǔn)差.

王懷野等人[21]充分分析了各向異性高斯濾波器σu和σv和濾波尺度設(shè)定問題，對(duì)于圖像局部區(qū)域方差較小的平緩區(qū)域，σu和σv的比例應(yīng)該盡可能地接近1，在圖像局部區(qū)域方差較大區(qū)域，濾波器退化為沿著圖像邊緣切線方向的線狀濾波器，σu和σv的比例接近為0.

1.4 掩模尺寸

為了提高雙邊濾波的背景抑制效果，Bae等人[22]介紹了一種基于邊緣濾波的掩模尺寸自適應(yīng)調(diào)節(jié)方法.該方法首先累積目標(biāo)像元鄰域內(nèi)水平、垂直和對(duì)角線方向上的灰度變化，并依據(jù)灰度變化的大小選擇合適的掩模尺寸.具體的掩模尺寸選擇方法如公式（11）所示：

2 實(shí) 驗(yàn)

2.1 性能評(píng)價(jià)標(biāo)準(zhǔn)

為了衡量背景抑制效果，評(píng)價(jià)算法性能，我們采用信雜比（SCR）、信雜比增益（GSCR）、背景抑制算子（BSF）指標(biāo)對(duì)算法進(jìn)行評(píng)價(jià)，各評(píng)價(jià)指標(biāo)定義如下：

依據(jù)弱小目標(biāo)的定義，本文將目標(biāo)區(qū)域尺寸設(shè)定為3*3，以7*7區(qū)域估算背景能量，考慮到目標(biāo)能量擴(kuò)散對(duì)背景抑制的影響，設(shè)定5*5區(qū)域?yàn)槟繕?biāo)與背景的過渡區(qū)，如圖2所示.

2.2 實(shí)驗(yàn)分析

為了驗(yàn)證本文背景抑制算法的有效性，在實(shí)際拍攝得到的紅外場(chǎng)景1、紅外場(chǎng)景2、紅外場(chǎng)景3和紅外場(chǎng)景4中添加仿真設(shè)計(jì)的弱小目標(biāo)得到實(shí)驗(yàn)輸入，同時(shí)對(duì)比各向異性雙邊濾波、各向同性雙邊濾波、R-NLM[23]以及改進(jìn)的2DLMS[24]的背景抑制效果，以驗(yàn)證本文濾波算法的先進(jìn)性.

實(shí)驗(yàn)中以圖2所示的背景區(qū)域設(shè)定方法統(tǒng)計(jì)圖像中所有像元鄰域的標(biāo)準(zhǔn)差，取k個(gè)最小標(biāo)準(zhǔn)差的均值作為最終的圖像噪聲標(biāo)準(zhǔn)差σi.σs=2.5σi，σu=1.5σi，σu=R×σv，本文所提各向異性雙邊濾波掩模的具體計(jì)算過程如圖3所示.

圖3以紅外場(chǎng)景2中的云層邊緣區(qū)域的中心像元為目標(biāo)像元，如左圖放大區(qū)域的中心點(diǎn).紅外場(chǎng)景2中統(tǒng)計(jì)的σi約為0.25，目標(biāo)像元鄰域的標(biāo)準(zhǔn)差約為30，R=0.389 2，M=3，sobel算子估算的邊緣方向角為135°.雙邊濾波中輸入圖像的動(dòng)態(tài)范圍為[0，1]，故圖中顯示的目標(biāo)鄰域像素灰度值已歸一化到[0，1].

完整的實(shí)驗(yàn)結(jié)果如圖4和圖5所示，其中（a）是原始紅外場(chǎng)景，（b），（c），（d），（e）分別是各向同性雙邊濾波、改進(jìn)的2DLMS、Nl-means和各向異性雙邊濾波的背景殘擦mesh圖.對(duì)比圖4（b）、圖5（b）與圖4（e）、圖5（e），采用各向同性的雙邊濾波背景抑制算法，濾波方向始終沿著圖像的水平和垂直方向，背景抑制后，在圖像的邊緣區(qū)域存在較強(qiáng)的背景殘留，采用各向異性雙邊濾波的背景抑制算法能夠根據(jù)圖像中的區(qū)域統(tǒng)計(jì)量和方向信息自適應(yīng)調(diào)節(jié)濾波器的濾波尺度和濾波方向，濾波后邊緣殘留較少，對(duì)依然存在的強(qiáng)邊緣殘留，可以通過后續(xù)多幀確認(rèn)的方法剔除這些虛警.

圖4（c）、圖5（c）和圖4（d）、圖5（d）是2DLMS、R-NLM濾波背景抑制輸出，相對(duì)于各向異性雙邊濾波，利用預(yù)測(cè)像元鄰域方差信息自適應(yīng)調(diào)節(jié)濾波尺寸的2DLMS和利用時(shí)域灰度值相似性的R-NLM算法同樣能夠較好地抑制邊緣，但是上述濾波算法在進(jìn)行背景抑制的同時(shí)對(duì)目標(biāo)能量也產(chǎn)生了一定的削弱作用.

表1列出了各向異性雙邊濾波、各向同性雙邊濾波、改進(jìn)的2DLMS濾波、NL-means濾波在四種場(chǎng)景下輸出的目標(biāo)區(qū)域局部統(tǒng)計(jì)量，各向異性雙邊濾波、改進(jìn)的2DLMS濾波、NL-means濾波紅外背景抑制后的圖像的信雜比增益和背景抑制算子都要高于各向同性雙邊濾波，且三種濾波算法輸出的目標(biāo)區(qū)域的背景抑制算子相差不大，但由于改進(jìn)的2DLMS濾波、NL-means濾波保留的目標(biāo)能量要弱于各向異性雙邊濾波算法，兩種濾波算法背景抑制后的信雜比增益要低于各向異性雙邊濾波.

3 總 結(jié)

本文基于遙感圖像的局部統(tǒng)計(jì)信息和云層邊緣的方向性特征設(shè)計(jì)的各向異性雙邊濾波背景抑制算法能夠自適應(yīng)地調(diào)節(jié)濾波尺度和濾波方向.仿真實(shí)驗(yàn)結(jié)果表明，該方法可以在盡可能保留目標(biāo)能量的同時(shí)更好地抑制背景邊緣，背景抑制后的目標(biāo)局部信雜比增益大于4.

參考文獻(xiàn)

[1] 孫海江，王延杰，陳小林. 基于自適應(yīng)梯度閾值各向異性濾波抑制紅外復(fù)雜背景[J]. 光學(xué)精密工程，2014，22（1）：146-151.

SUN H J， WANG Y J， CHEN X L.Suppression of infrared complex background based on adaptive gradient threshold anisotropic filtering[J].Optics and Precision Engineering，2014，22（1）：146-151.（In Chinese）

[2] 朱斌，樊祥，馬東輝，等. 天空紅外圖像非平穩(wěn)背景的小波域抑制[J]. 紅外與激光工程，2008，（S2）：634-637.

ZHU B， FAN X，MA D H，et al. Non-stable background suppress in wavelet domain of infrared image of sky[J].Infrared and Laser Engineering， 2008，（S2）：634-637.（In Chinese）

[3] 程塨，郭雷，韓軍偉，等. 基于形態(tài)學(xué)帶通濾波和尺度空間理論的紅外弱小目標(biāo)檢測(cè)[J]. 光學(xué)學(xué)報(bào)，2012，32（10）：152-159.

CHENG G， GUO L， HAN J W， et al. Infrared dim small target detection based on morphological band-pass filtering and scale space theory[J]. Acta Optica Sinica，2012，32（10）：152-159.（In Chinese）

[4] 孫立輝，王永仲，周冰.基于局部統(tǒng)計(jì)特征的自適應(yīng)紅外背景抑制算法[J].紅外與激光工程，2008，37（1）：177-180.

SUN L H， WANG Y Z， ZHOU B. Self-adaptive infrared image background suppressed algorithm based on local statistical character[J]. Infrared and Laser Engineering， 2008，37（1）：177-180.（In Chinese）

[5] 周慧鑫，秦翰林，賴睿，等. 基于多尺度截?cái)嗟娜跣∧繕?biāo)復(fù)雜背景抑制[J]. 光學(xué)學(xué)報(bào)，2010，30（10）：2812-2816.

ZHOU H X， QIN H L， LAI R， et al. Multiscale truncation for dim and small target background suppression[J]. Acta Optica Sinica，2010，30（10）：2812-2816. （In Chinese）

[6] HAN Q. Infrared dim and small target background suppression using shearlettransform[J]. Journal of Infrared and Millimeter Waves， 2011，30 （2） ：162-166.

[7] 秦翰林，李佳，周慧鑫，等. 采用剪切波變換的紅外弱小目標(biāo)背景抑制[J]. 紅外與毫米波學(xué)報(bào)，2011，30（2）：162-166.

QIN H L， LI J， ZHOU H X， et al.Infrared dim and small target background suppression using shearlettransform[J].J Infrared Millim Waves， 2011，30（2）：162-166.（In Chinese）

[8] 榮生輝，劉剛，周慧鑫，等. 基于改進(jìn)的剪切波變換和引導(dǎo)濾波的紅外弱小目標(biāo)背景抑制方法[J]. 光子學(xué)報(bào)，2015，44（2）：68-73.

RONG S H， LIU G，ZHOU H X， et al. Infrared dim and small target background suppression based on the improved shearlet transformed the guide filter[J].Acta Photonic Sincia， 2015，44（2）：68-73.（In Chinese）

[9] 秦翰林，周慧鑫，劉上乾，等.一種基于非下采樣Contourlet變換的運(yùn)動(dòng)點(diǎn)目標(biāo)背景抑制算法[J]. 光子學(xué)報(bào)，2009，38（12）：3318-3321.

QIN H L，ZHOU H X， LIU S Q，et al. An algorithm for moving point target background suppression based on NSCT[J]. Acta Photonic Sincia， 2009，38（12）：3318-3321.（In Chinese）

[10]QIN H， HAN J， YAN X， et al. Multiscale random projection based background suppression of infrared small target image[J]. Infrared Physics & Technology， 2015， 73： 255-262.

[11]GONG J， ZHANG Y， HOU Q， et al. Background suppression for cloud clutter using temporal difference projection[J]. Infrared Physics & Technology， 2014，64： 66-72.

[12]張媛，辛云宏，張春琴. 基于時(shí)空聯(lián)合濾波技術(shù)的緩慢運(yùn)動(dòng)紅外弱小目標(biāo)檢測(cè)算法[J]. 光子學(xué)報(bào)，2010，39（11）：2049-2054.

ZHANG Y， XIN Y H， ZHANG C Q. An algorithm based on temporal and spatial filters for infrared weak slow moving point target detection[J].Acta Photonic Sincia， 2010， 39（11）：2049-2054.（In Chinese）

[13]李凡， 劉上乾， 秦翰林，等. 自適應(yīng)雙邊濾波紅外弱小目標(biāo)檢測(cè)方法[J].光子學(xué)報(bào)， 2010， 39（6）： 1129-1131.

LI F， LIU S Q， QIN H L， et al. Dim infrared targets detection based on adaptive bilateral filter[J]. Acta Photonic Sincia， 2010， 39（6）： 1129-1131.（In Chinese）

[14]曹瑛， 李志永， 盧曉鵬， 等. 基于自適應(yīng)鄰域雙邊濾波的點(diǎn)目標(biāo)檢測(cè)預(yù)處理算法[J].電子與信息學(xué)報(bào)， 2008， 30（8）： 1909-1912.

CAO Y， LI Z Y， LU X P，et al.A preprocessing algorithm of point target detection based on temporal-spatial bilateral filter using adaptive neighborhoods[J]. Journal of Electronics & Information Technology， 2008， 30（8）： 1909-1912.（In Chinese）

[15]AURICH V， WEULE J. Non-linear gaussian filters performing edge preserving discussion[C]// Proceedings of the DAGM Symposium.1995： 538-545.

[16]SMITH S M， BRADY J M. SUSAN — A new approach to low level image processing[J]. International Journal of Computer Vision， 1997， 23（1）： 45-78.

[17]TOMASI C， MANDUCHI R. Bilateral filtering for gray and color images [C]//Proceedings of the IEEE International Conference on Computer Vision. 1998：839-846.

[18]GEUSEBROEK J， SMEULDERS A W， De WEIJER J V， et al. Fast anisotropic Gauss filtering[J]. IEEE Transactions on Image Processing， 2003， 12（8）： 938-943.

[19]KUGU E. Satellite image denoising using bilateral filter with SPEA2 optimized parameters[C]//Recent Advances in Space Technologies （RAST），6th International Conference on. 2013：217-223.

[20]PENG H， RAO R. Bilateral kernel parameter optimization by risk minimization[C]//International Conference on Image Processing. 2010： 3293-3296.

[21]王懷野， 張科， 李言俊， 等. 各向異性濾波在紅外圖像處理中的應(yīng)用[J].紅外與毫米波學(xué)報(bào)， 2005， 24（2）： 109-113.

WANG H Y， ZHANG K， LI Y J， et al. Anisotropic Guassian filtering for infrared image[J]. J Infrared Millim Waves， 2005， 24（2）： 109-113.（In Chinese）

[22]BAE T， SOHNG K. Small target detection using bilateral filter based on edge component[J]. Journal of Infrared， Millimeter， and Terahertz Waves， 2010， 31（6）： 735-743.

[23]HU J， YU Y， LIU F， et al. Small and dim target detection by background estimation[J]. Infrared Physics & Technology， 2015， 73（73）： 141-148.

[24]BAE T W， ZHANG F， KWEON I S， et al. Edge directional 2D LMS filter for infrared small target detection[J]. Infrared Physics & Technology， 2012， 55（1）： 137-145.