劉剛+王明靜+周珩+張喜濤

摘要： 針對紅外圖像易受噪聲污染、 成像質(zhì)量差等問題， 提出一種非下采樣Contourlet域變換系數(shù)混合統(tǒng)計建模的降噪方法。 將圖像變換到非下采樣Contourlet域， 對噪聲系數(shù)、 信號系數(shù)分別按照高斯分布和廣義拉普拉斯分布建模， 利用貝葉斯框架下的最大后驗估計理論， 推導(dǎo)原始信號系數(shù)估計公式， 并對包含噪聲的圖像系數(shù)進(jìn)行處理， 實現(xiàn)非下采樣Contourlet域紅外圖像降噪。 實驗結(jié)果表明， 該方法能夠?qū)t外成像過程中產(chǎn)生的高斯噪聲實現(xiàn)有效抑制， 較完整地保持圖像的邊緣等細(xì)節(jié)信息， 在峰值信噪比提高與視覺效果上優(yōu)于部分經(jīng)典算法。

關(guān)鍵詞： 紅外圖像； 降噪； 非下采樣Contourlet變換； 拉普拉斯分布

中圖分類號： TP391.41文獻(xiàn)標(biāo)識碼： A文章編號： 1673-5048（2016）05-0066-05

Abstract： For the denoising problem of infrared image， an algorithm is proposed based on mixed statistical model of transform coefficients in nonsubsampled Contourlet domain. The image is transformed into nonsubsampled Contourlet domain， and noise coefficients and signal coefficients are modeled respectively by Gaussian distribution and generalized Laplacian distribution. By using maximum a posteriori estimation theory under Bayesian framework， a new estimation formula for the original signal coefficients is deduced， and the noise image coefficients are processed to realize the infrared image denoising of nonsubsampled Contourlet domain. The experimental results show that the proposed method can suppress the Gaussian noise effectively and keep most of the images detail information， which has higher value of peak signal noise ratio and better visual effect than some standard algorithms.

Key words： infrared image； denoising； nonsubsampled Contourlet transform； Laplacian distribution

0引言

紅外圖像噪聲豐富， 而噪聲造成圖像退化、 圖像特征被掩蓋， 直接影響圖像分割、 特征提取等后續(xù)工作的準(zhǔn)確性， 因此， 抑制紅外圖像噪聲、 提升圖像質(zhì)量是紅外成像應(yīng)用中圖像處理和分析的前提。

小波變換具有良好的時頻局域化特性和多分辨率特性[1]， 在圖像降噪中得到了廣泛的應(yīng)用[2-5]。 小波變換在表示具有點奇異性的函數(shù)時是最優(yōu)基， 但是對于更高維的特征則顯得力不從心。 在二維圖像中， 邊緣、 輪廓和紋理等具有高維奇異性的幾何特征包含了大部分信息， 對此， 小波不是表示圖像的最優(yōu)基。

Do等人提出的Contourlet變換是一種真正的圖像二維表示方法， 其將小波的優(yōu)點延伸到高維空間， 能夠更好地刻畫高維信息的特性， 更適合處理具有超平面奇異性的信息[6-7]。 但是， Contourlet變換不具備平移不變性， 在利用其進(jìn)行圖像去噪時， 奇異點周圍會引入偽吉布斯（Gibbs）現(xiàn)象。 此外， Contourlet變換存在下采樣過程， 因此低頻子帶和高頻子帶均存在頻譜混疊現(xiàn)象。 頻譜混疊造成同一方向的信息會在幾個不同方向子帶中同時出現(xiàn)， 這在一定程度上削弱了其方向選擇性。

Cunha等人[8-9]利用非下采樣塔式分解和非下采樣濾波器組構(gòu)造出了非下采樣Contourlet變換（Nonsubsampled Contourlet Transform， NSCT）， 具有平移不變性。 NSCT繼承了Contourlet變換的良好特性， 變換后系數(shù)能量更加集中， 能夠更好地捕捉和跟蹤圖像中重要的幾何特征。 同時， 由于沒有采樣過程， 圖像的分解和重構(gòu)過程中不具有頻率混疊項， 這使得NSCT具有平移不變性。

在圖像降噪中應(yīng)用Contourlet和NSCT， 最常用的是傳統(tǒng)的閾值方法[10-14]。 該類方法的缺點在于沒有考慮到各圖像尺度內(nèi)系數(shù)的相關(guān)性、 相鄰尺度系數(shù)之間的相關(guān)性以及邊緣輪廓的關(guān)系， 圖像細(xì)節(jié)信息不能很好地保留且存在振鈴現(xiàn)象。 建立變換系數(shù)尺度內(nèi)和尺度間的統(tǒng)計模型是在變換域還原信號系數(shù)、 抑制噪聲系數(shù)的關(guān)鍵問題[15]。

紅外成像應(yīng)用中的紅外探測器是系統(tǒng)噪聲的主要來源， 是影響紅外圖像質(zhì)量的主要因素， 其強(qiáng)度一般情況下遠(yuǎn)大于其他環(huán)節(jié)產(chǎn)生的噪聲， 同時也是最難以克服的。 探測器本身的噪聲是無法避免的， 按照其產(chǎn)生的機(jī)理可分為散粒噪聲、 熱噪聲、 光子噪聲、 產(chǎn)生復(fù)合噪聲和1/f噪聲等。 其中散粒噪聲、 熱噪聲、 光子噪聲和產(chǎn)生復(fù)合噪聲所產(chǎn)生的探測器電流輸出是一個隨機(jī)過程， 由中心極限定理近似服從高斯分布。 1/f噪聲是紅外探測器低頻部分的一種電流噪聲， 顧名思義， 1/f噪聲與頻率成反比， 當(dāng)高于一定頻率時， 與其他噪聲相比可忽略不計。 因此， 可認(rèn)為影響紅外圖像質(zhì)量的噪聲服從高斯分布[16]。

航空兵器2016年第5期劉剛等： 非下采樣輪廓波域混合統(tǒng)計模型紅外圖像降噪方法基于以上分析， 在NSCT域提出了一種新降噪方法。

1NSCT域混合統(tǒng)計圖像降噪模型

1.1NSCT理論

NSCT與Contourlet變換一樣， 也是采用由拉普拉斯金字塔變換與方向濾波器組所構(gòu)成的雙迭代濾波器組結(jié)構(gòu)。 但NSCT采用的是非下采樣塔式濾波器和非下采樣方向濾波器組， 變換時首先由非下采樣塔式濾波器將圖像分解為低通（低頻）部分和帶通（高頻）部分， 然后由非下采樣方向性濾波器組將帶通部分分解為若干個方向。 NSCT去掉了塔式分解和方向性濾波器組分解中信號經(jīng)分析濾波后的下采樣（抽?。┮约熬C合濾波前的上采樣（插值）， 而改為對相應(yīng)的濾波器進(jìn)行采樣， 再對信號進(jìn)行分析濾波和綜合濾波。

1.2混合統(tǒng)計降噪模型

NSCT域信號系數(shù)尺度間和尺度內(nèi)同方向子帶內(nèi)相關(guān)性強(qiáng)， 噪聲系數(shù)層間相關(guān)性較弱、 層內(nèi)無相關(guān)性。 基于此， 將NSCT系數(shù)分為噪聲系數(shù)和信號系數(shù)兩類并分別建模， 根據(jù)貝葉斯框架下的最大后驗估計理論， 推導(dǎo)降噪方程。

信號系數(shù)包含了圖像的大部分信息及主要的邊緣和紋理， 相鄰尺度間的系數(shù)相關(guān)性較強(qiáng)， 分解系數(shù)的統(tǒng)計分布在原點處的峰值更尖銳， 并具有長拖尾， 可采用廣義的拉普拉斯分布來建模。 噪聲系數(shù)主要是噪聲和微小細(xì)節(jié)的貢獻(xiàn)， 采用具有強(qiáng)局部相關(guān)的零均值高斯分布來建模。

從表1、 圖1～3可得出如下結(jié)論：

（1） 從表1中可知， 基于NSCT域系數(shù)混合統(tǒng)計模型降噪算法的PSNR最高， 能有效去除圖像中的高斯白噪聲， 較好地保持圖像的邊緣， 在噪聲強(qiáng)度較大時， 降噪能力更加明顯；

（2） 從圖1、 圖2和圖3可知， 基于NSCT域系數(shù)混合統(tǒng)計模型的降噪算法既濾除了圖像噪聲， 又有效地保持了相關(guān)邊緣細(xì)節(jié)；

（3） 基于Contourlet域混合統(tǒng)計模型降噪算法在視覺效果上同NSCT域方法接近， 但在PSNR上比NSCT域方法稍差；

（4） 小波軟閾值法在濾除噪聲的同時造成了邊緣細(xì)節(jié)的模糊；

（5） Contourlet閾值降噪方法邊緣細(xì)節(jié)部分保持較好， 也能有效濾除噪聲， 但是出現(xiàn)了較為明顯的柵格效應(yīng)。

4結(jié)論

在NSCT域?qū)υ肼曄禂?shù)、 信號系數(shù)分別按照高斯分布和廣義拉普拉斯分布建模， 根據(jù)貝葉斯框架下的最大后驗估計理論， 推導(dǎo)原始信號系數(shù)的估計公式， 對包含噪聲的圖像系數(shù)進(jìn)行處理， 實現(xiàn)NSCT域紅外圖像降噪。 NSCT具有平移不變性和更好的方向選擇特性， 實驗結(jié)果充分表明了本文提出的基于NSCT域系數(shù)混合統(tǒng)計模型的降噪方法在PSNR提高與視覺效果上優(yōu)于部分經(jīng)典降噪方法。

參考文獻(xiàn)：

[1] 高靜， 李朝偉， 董云峰， 等.空空導(dǎo)彈導(dǎo)引頭小波降噪?yún)?shù)優(yōu)選仿真研究[J].航空兵器， 2010（5）： 48-54.

[2] 徐小軍， 王友仁. 基于離散分?jǐn)?shù)階正交小波變換圖像降噪新方法[J].電子學(xué)報， 2014， 42（2）： 280-287.

[3] 張?zhí)铮?孫延奎， 田小林.二進(jìn)小波與擴(kuò)散濾波結(jié)合的光學(xué)相干層析圖像降噪[J].吉林大學(xué)學(xué)報（工學(xué)版）， 2013， 43（S1）： 340-344.

[4] 金煒， 周亞訓(xùn)， 符冉迪， 等.采用零樹結(jié)構(gòu)分類小波系數(shù)的紅外圖像降噪[J].光電工程， 2012， 39（5）： 79-83.

[5] Donoho D L. DeNoising by SoftThresholding[J]. IEEE Transactions on Information Theory， 1995， 41（3）： 613-627.

[6] Po D DY， Do M N. Directional Multiscale Modeling of Images Using the Contourlet Transform[J].IEEE Transactions on Image Processing， 2006， 15（6）： 1610-1620.

[7] Do M N， Vetterli M. The Contourlet Transform： An Efficient Directional Multiresolution Image Representation[J]. IEEE Transastions on Image Processing， 2005， 14（12）： 2091-2106.

[8] Cunha A L， Zhou Jianping， Do M N. The Nonsubsampled Contourlet Transform： Theory， Design and Applications[J].IEEE Transactions on Image Processing， 2006， 15（10）： 3089-3101.

[9] Zhou Jianping， Cunha A L， Do M N. Nonsubsampled Contourlet Transform： Construction and Application in Enhancement[C]∥Proceedings of the IEEE International Conference on Image Processing， 2005： 469-472.

[10] 戴維， 于盛林， 孫栓. 基于Contourlet變換自適應(yīng)閾值的圖像去噪算法[J].電子學(xué)報， 2007， 35（10）： 1939-1943.

[11] 劉麗萍， 崔寧海， 王琰.一種新的Contourlet變換圖像降噪算法[J].沈陽理工大學(xué)學(xué)報， 2012， 31（1）： 31-35.

[12] 程光權(quán)， 成禮智.基于冗余Contourlet變換的圖像相關(guān)法去噪[J].中國圖象圖形學(xué)報， 2008， 13（9）： 1678-1682.

[13] 倪超， 李奇， 夏良正.基于非下采樣Contourlet擴(kuò)散濾波的紅外圖像降噪[J].東南大學(xué)學(xué)報（自然科學(xué)版）， 2010， 40（S1）： 272-276.

[14] 閆晟， 劉明剛， 原建平， 等. 一種基于NSCT變換的超聲圖像降噪新算法[J].儀器儀表學(xué)報， 2012， 33（5）： 1005-1012.

[15] 王文， 芮國勝， 王曉東， 等.圖像多尺度統(tǒng)計模型綜述[J].中國圖象圖形學(xué)報， 2007， 12（6）： 961-968.

[16] 陳維真， 周曉東， 張春華.基于探測器噪聲分析的紅外圖像增強(qiáng)算法[J].紅外與激光工程， 2008， 37（S1）： 44-47.