劉艾琳

(天津電子信息職業(yè)技術(shù)學(xué)院數(shù)字藝術(shù)系，天津300350)

引 言

紅外圖像總體上對(duì)比度較低，且由于成像環(huán)境復(fù)雜以及在傳輸存儲(chǔ)等環(huán)節(jié)中，不可避免受到噪聲的干擾，導(dǎo)致所獲取的圖像一般來說是降質(zhì)圖像。KANG［1］等人基于紅外圖像的自相似統(tǒng)計(jì)特征提出一種兩階段3維濾波算法;WANG［2］等人將偽中值濾波與小波變換相結(jié)合用于紅外圖像處理;YU［3］等人提出一種基于混合雙邊濾波的紅外圖像自適應(yīng)濾波算法;HU［4］等人將同態(tài)濾波用于紅外圖像增強(qiáng)研究，取得了不錯(cuò)的效果;ZHANG［5］提出一種針對(duì)紅外圖像中的顆粒噪聲的濾波算法。近年來該領(lǐng)域的研究成果特點(diǎn)有:(1)針對(duì)紅外圖像濾波增強(qiáng)的研究，主要基于空間域進(jìn)行，較少利用多尺度分析方法，如小波變換、提升小波變換等;(2)對(duì)于紅外圖像濾波與增強(qiáng)基本是分開研究，并未作為一個(gè)整體進(jìn)行，即對(duì)圖像噪聲的濾波，較少考慮到濾波后圖像視覺效果;對(duì)圖像增強(qiáng)處理并未考慮到圖像受到噪聲干擾這一特殊情形。

提升小波變換與傳統(tǒng)建立在傅里葉變換基礎(chǔ)上的小波變換的區(qū)別主要有:(1)傳統(tǒng)小波變換通過將小波母函數(shù)伸縮或者平移來構(gòu)造小波基函數(shù)，運(yùn)算效率較為低下;而提升小波變換采用“分步驟”的思想來實(shí)現(xiàn)小波變換，即采用預(yù)測(cè)和更新兩個(gè)環(huán)節(jié)實(shí)現(xiàn)對(duì)圖像高頻和低頻信息的高效率分離，計(jì)算效率得到提高;(2)提升小波變換通過原位操作，能夠?qū)崿F(xiàn)圖像的整數(shù)小波變換，在分解和重構(gòu)過程中對(duì)圖像信號(hào)刻畫更為精確，這也是傳統(tǒng)小波變換所不具備的。近年來，提升小波變換以其優(yōu)良的特性在圖像濾波［6］、圖像復(fù)原［7］、圖像重構(gòu)［8］等方面得到成功的應(yīng)用，但就目前公開發(fā)表的研究成果來看，提升小波變換在紅外圖像處理方面應(yīng)用較少。

作者將圖像濾波和增強(qiáng)方法有機(jī)結(jié)合，借助提升小波變換理論框架，提出了一種基于提升小波變換的紅外降質(zhì)圖像的雙重濾波方法。

1 紅外圖像提升小波變換流程分析

一幅紅外圖像可表示成一個(gè)信號(hào)序列的集合Sj={Sj，k}，對(duì)該圖像采用單層提升小波變換，主要流程如下。

(1)分解。將上述信號(hào)集合分解成偶序列η2j和奇數(shù)序列φ2j-1，并且該兩類序列彼此互不相交，該步驟中圖像信號(hào)序列可表示成:

(2)預(yù)測(cè)。該步驟是基于圖像信號(hào)間存在相關(guān)性這一前提，采用相鄰信號(hào)偶數(shù)序列對(duì)奇數(shù)序列進(jìn)行預(yù)測(cè)，通過對(duì)奇數(shù)序列的真實(shí)值與預(yù)測(cè)值做差，該差值可稱為細(xì)節(jié)系數(shù)，預(yù)測(cè)步驟可定義為:

式中，D2j-1為細(xì)節(jié)系數(shù)，P［·］代表預(yù)測(cè)算子。

(3)更新。采用更新算子U［·］對(duì)第(2)步中產(chǎn)生的序列D2j-1進(jìn)行更新，其過程為:

2 一種新型小波域閾值函數(shù)模型

近年來，在經(jīng)典硬閾值、軟閾值函數(shù)模型的基礎(chǔ)上，提出了一系列的改進(jìn)函數(shù)模型。

(1)折中法的改進(jìn)閾值模型，該模型通過對(duì)增加一定的調(diào)節(jié)因子，對(duì)經(jīng)典的軟、硬閾值進(jìn)行改進(jìn)，在一定程度上提高了濾波效果，該閾值模型代表性的一類如下表示:

(2)指數(shù)型改進(jìn)閾值函數(shù)模型，該類模型通過對(duì)經(jīng)典閾值函數(shù)模型中增加指數(shù)型的調(diào)節(jié)系數(shù)，來實(shí)現(xiàn)對(duì)濾波效果的改進(jìn)，代表性的函數(shù)模型如下所示:

上述兩類函數(shù)都是通過增加調(diào)節(jié)因子來實(shí)現(xiàn)對(duì)經(jīng)典閾值函數(shù)存在缺陷進(jìn)行一定程度的修正。不足之處在于:(1)模型中的調(diào)節(jié)因子的功能是對(duì)函數(shù)模型的過度濾波現(xiàn)象進(jìn)行修正，而事實(shí)上圖像受到的噪聲強(qiáng)度是不同的，采用千篇一律的調(diào)節(jié)因子，濾波效果沒有從根本上得到提高;(2)閾值的單一性，特別是對(duì)于提升小波變換而言，圖像經(jīng)過多層小波變換后的小波系數(shù)幅值會(huì)隨著分解層數(shù)的增大而快速衰減，而相應(yīng)的閾值并未隨之進(jìn)行改變。鑒于紅外圖像噪聲強(qiáng)度大，且對(duì)比度較低，結(jié)合以上對(duì)兩類改進(jìn)閾值函數(shù)模型特性分析結(jié)果，提出一種新型提升小波變換的閾值函數(shù)模型:

式中，N為小波分解層數(shù)，T*為新型小波系數(shù)閾值。經(jīng)典軟硬閾值以及上述(4)式和(5)式描述的改進(jìn)型閾值函數(shù)均采用如下閾值選取策略:

式中，x和y分別為圖像尺寸大小，σ為小波分解系數(shù)的均方差:

可以看出，該閾值屬于全局閾值范疇，即無法根據(jù)小波分解層數(shù)的變化而自適應(yīng)進(jìn)行調(diào)整，理論上講，隨著小波分解層數(shù)的增加，小波系數(shù)幅值快速衰減，那么，對(duì)應(yīng)的閾值也應(yīng)當(dāng)呈現(xiàn)這一特征，對(duì)此，作者提出一種新型閾值計(jì)算方法:

3 改進(jìn)非均局部均值濾波算法

紅外噪聲圖像可采用如下模型表示:

式中，X(i)為紅外圖像有用信號(hào)，n(i)則為噪聲信號(hào)。對(duì)于噪聲圖像中任一像素點(diǎn)，采用非局部均值濾波進(jìn)行處理時(shí)，是通過求取該圖像中所有像素點(diǎn)灰度值的加權(quán)平均值作為該點(diǎn)的濾波結(jié)果:

式中，X′(i)為濾波后的紅外圖像信號(hào)，n(i)為噪聲信號(hào)，w(i，j)為權(quán)值，且j取任意正整數(shù)。權(quán)值w(i，j)是根據(jù)圖像上任意兩個(gè)像素點(diǎn)間的相似程度來進(jìn)行賦值，而兩點(diǎn)間的相似程度由兩者的灰度值矩陣進(jìn)行決定。一般來說，采用兩點(diǎn)間的灰度值矩陣的歐氏距離來衡量兩者的相似程度，即:式中，Ni，Nj分別表示像素 i，j對(duì)應(yīng)的灰度值矩陣，a為高斯加權(quán)標(biāo)準(zhǔn)差值。那么，權(quán)值w(i，j)則可進(jìn)行如下計(jì)算:式中是歸一化參量，f(·)為核函數(shù)，在某種程度上決定著該算法的成敗，對(duì)于經(jīng)典非局部均值濾波而言，該核函數(shù)為:

式中，h為衰減因子，其功能在于控制函數(shù)的衰減，從而影響最終的濾波效果。函數(shù)的類型多種多樣，一般來說主要有余弦型和高斯型兩類核函數(shù):

相關(guān)研究成果表明［9］，經(jīng)典非局部均值中的指數(shù)型核函數(shù)對(duì)于強(qiáng)度較低的噪聲濾波效果較為理想，當(dāng)噪聲持續(xù)增大時(shí)，效果則出現(xiàn)明顯下降;而余弦型核函數(shù)濾波性能則與之相反，兩者具有一定的互補(bǔ)作用。因此，本文中將兩者有機(jī)結(jié)合，提出一種基于新型核函數(shù)模型，該模型如下表示:

該模型對(duì)經(jīng)典非局部均值濾波算法中的指數(shù)型核函數(shù)增加了一個(gè)余弦型核函數(shù)作為其系數(shù)，以此來對(duì)原有的經(jīng)典指數(shù)核函數(shù)濾波性能進(jìn)行一定程度的提升。

4 算法基本流程及實(shí)驗(yàn)仿真分析

4．1 算法基本流程

步驟1:對(duì)紅外噪聲圖像進(jìn)行提升小波變換，獲得噪聲圖像的低頻和高頻提升小波變換子圖像。

步驟2:對(duì)低頻提升小波變換子圖像再次執(zhí)行提升小波變換，獲得次高頻分解子圖像1和次低頻分解子圖像1，鑒于紅外圖像絕大部分的背景信息主要集中于次低頻分解子圖像1中，該部分信息的存在反而會(huì)導(dǎo)致紅外圖像整體上偏暗，少量的有價(jià)值的目標(biāo)信息被湮沒，因而加以舍棄;采用(6)式所定義的自適應(yīng)小波閾值去噪算法來對(duì)抑制該部分圖像中的噪聲。

步驟3:對(duì)高頻提升小波變換高頻子圖像再次執(zhí)行提升小波變換，獲得次高頻分解子圖像2和次低頻分解子圖像2，鑒于噪聲主要集中于次高頻分解子圖像2中，因而加以舍棄;采用改進(jìn)非局部均值濾波算法去除次低頻分解子圖像2中存在的噪聲。

步驟4:將步驟2和步驟3中濾波后的提升小波變換子圖像進(jìn)行重構(gòu)，獲得去噪后的紅外圖像。

步驟5:對(duì)步驟4中獲得的去噪圖像采用直方圖均衡化算法處理，以提高去噪后圖像的視覺效果。

4．2 實(shí)驗(yàn)仿真分析

采用Lena標(biāo)準(zhǔn)測(cè)試圖像以及一幅紅外圖像來對(duì)本文中的算法進(jìn)行性能測(cè)試。此外，采用經(jīng)典非局部均值濾波算法(non-local means，NLM)、小波閾值函數(shù)濾波算法［6］以及改進(jìn)型非局部均值濾波算法(improved non-local means，INLM)［9］分別對(duì)上述圖像進(jìn)行去噪實(shí)驗(yàn)，并將其實(shí)驗(yàn)結(jié)果與本文中算法的實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比。對(duì)于濾波結(jié)果的評(píng)價(jià)，選用應(yīng)用較為廣泛的峰值信噪比(peak signal-to-noise ratio，PSNR)［10］以及結(jié)構(gòu)相似度(structural similarity，SSIM)［11］作為濾波后圖像質(zhì)量定量評(píng)價(jià)指標(biāo)。Lena標(biāo)準(zhǔn)測(cè)試圖像實(shí)驗(yàn)結(jié)果見圖1、表1和表2。

Fig.1 Comparison of the denoising results of Lena imagea—noise image(30%) b—NLM c—algorithm［6］ d—algorithm［9］e—algorithm proposed

Table 1 Quantitative comparison of PSNR value of the results of Lena image/dB

從圖1b可以看出，圖像中噪聲殘留程度較高，圖像相對(duì)于噪聲圖像(見圖1a)而言，視覺效果并未得到改善;采用參考文獻(xiàn)［6］中所定義的輪廓波變換域新型閾值函數(shù)進(jìn)行濾波，從圖1c中可以看出，噪聲得到抑制，但仍有大量殘余，相對(duì)于經(jīng)典非局部均值濾波算法而言，改進(jìn)小波閾值算法［6］性能還是較優(yōu)的;從圖1d中可以看出，改進(jìn)非局部均值濾波算法［9］性能不但優(yōu)于經(jīng)典非均值濾波算法而且優(yōu)于改進(jìn)小波閾值算法［6］。但就濾波后圖像以及表1、表2中的數(shù)據(jù)可以看出，本文中算法性能總體上優(yōu)于前幾類算法，并且圖1e相對(duì)于圖1b～圖1d來說，圖像信息得到增強(qiáng)，對(duì)比度得以改善，這主要得益于本文中算法增加了一個(gè)圖像增強(qiáng)環(huán)節(jié)。

為了更好地用本文中的算法對(duì)紅外圖像的濾波性能進(jìn)行測(cè)試，仍采用上述幾類算法分別與本文中的算法對(duì)1幅熱紅外飛機(jī)圖像進(jìn)行實(shí)驗(yàn)，實(shí)驗(yàn)結(jié)果分別如圖2、圖3和表3、表4所示。

Fig.2 Comparison of the denoising results of plane imagea—noise image(20%) b—NLM c—algorithm［6］ d—algorithm［9］e—algorithm proposed

Fig.3 Comparison of the denoising results of plane imagea—noise image(40%) b—NLM c—algorithm［6］ d—algorithm［9］e—algorithm proposed

Table 3 Quantitative comparison of PSNR value of the results of plane image/dB

T e

結(jié)合表3、表4以及圖2、圖3的實(shí)驗(yàn)結(jié)果分析可以看出，在對(duì)紅外圖像(飛機(jī))加入10% ～40%強(qiáng)度的噪聲過程中，經(jīng)典非局部均值濾波算法、改進(jìn)小波閾值算法［6］、改進(jìn)非局部均值濾波算法［9］以及本文中算法性能均有所下降，但本文中算法性能下降幅度最小(PSNR值下降約4dB)，這說明，該算法具有較好的抗噪性。從圖2e和圖3e可以看出，圖中目標(biāo)物——飛機(jī)的輪廓基本從噪聲中恢復(fù)出來，且能清晰辨認(rèn)。

5 結(jié)束語

為了有效濾除紅外圖像中的噪聲，提出了一種基于提升小變換的紅外圖像雙重濾波算法。標(biāo)準(zhǔn)測(cè)試圖像和紅外圖像濾波結(jié)果證明，本文中算法的PSNR和SSIM等指標(biāo)值均優(yōu)于經(jīng)典非局部均值濾波、改進(jìn)小波閾值法以及改進(jìn)非局部均值濾波算法。就算法計(jì)算效率而言，盡管通過舍棄相當(dāng)一部分小波分解系數(shù)以提高算法計(jì)算速度，但由于對(duì)圖像先后實(shí)現(xiàn)兩次提升小波變換，無疑在一定程度上增加了算法耗時(shí)。因此，大幅度提高算法計(jì)算效率，是本文中算法需要改進(jìn)之處。

［1］ KANG Ch Q，CAO W P，HUA L，et al.Infrared image denoising algorithm via tow-stage 3-D filtering［J］.Laser ＆ Infrared，2013，43(3):261-264(in Chinese).

［2］ WANG X W，WANG Sh L，LI K.Infrared image enhancement based on pseudo median filter and wavelet transform［J］.Laser＆ Infared，2013，43(1):90-93(in Chinese).

［3］ YU B，GAO L，ZHAO T Y，et al.Adaptive hybrid bilateral filtering algorithm for infrared iamge［J］.Infrared and Laser Engineering，2012，41(11):3102-3107(in Chinese).

［4］ HU D M，ZHAO H Sh，LI Y Ch，et al.A new approach to infrared image enhancement based on homomorphic filter［J］.Infrared Technology，2012，34(4):224-228(in Chinese).

［5］ ZHANG C T.A new adaptive filtering algorithm based on salt＆ pepper noise in infrared image［J］.Infrared Technology，2013，35(8):502-506(in Chinese).

［6］ CAO X Y，ZHANG Zh J，XING J J.Method of radar signal denoising based on lifting waveelt improved threshold［J］.Computer Engineering and Applications，2012，48(14):143-147(in Chinese).

［7］ HUANG D T，WU Zh Y.Application of lifting wavelet transform in blind restoration scheme based on NAS-RIF algorithm［J］.Journal of Computer-Aided Design ＆ Computer Graphics，2012，24(12):1614-1620(in Chinese).

［8］ HOU X S，ZHANG L，XIAO L.Areconstruction algorithm with Bayesian compressive sensing for synthetic aperture radar images［J］.Journal of Xi’an Jiaotong University，2013，47(8):74-79(in Chinese).

［9］ LIU X M，TIAN Y，HE H，et al.Improved non-local means algorithm for image denoising［J］.Computer Engineering，2012，38(4):199-201(in Chinese).

［10］ LU Zh L，LI R L，LI T，et al.Infrared image denoising based total variation theory［J］.Laser Technology，2012，36(2):194-197(in Chinese).

［11］ TIAN H N，LI S M.Objective evaluation method for image quality based on edge structure similarity［J］.Acta Photonica Sinica，2013，42(1):110-114(in Chinese).