張 杰，周浦城，張 謙

?

基于迭代多尺度引導(dǎo)濾波Retinex的低照度圖像增強(qiáng)

張 杰1,2，周浦城1,2，張 謙1,2

(1. 陸軍炮兵防空兵學(xué)院，安徽 合肥 230031； 2. 偏振光成像探測技術(shù)安徽省重點(diǎn)實(shí)驗(yàn)室，安徽 合肥 230031)

為提高低照度圖像的視覺感知，提出了基于迭代多尺度引導(dǎo)濾波Retinex的低照度圖像增強(qiáng)算法。首先，將圖像轉(zhuǎn)換到Y(jié)CbCr空間；然后對Y分量用迭代多尺度引導(dǎo)濾波算法估計(jì)照度圖像，以指數(shù)形式迭代增加尺度參數(shù)與平滑參數(shù)，在保持邊緣的同時(shí)消除各種對比度細(xì)節(jié)；接著用照度圖像求得反射圖像，并進(jìn)行對比度校正；最后將圖像轉(zhuǎn)換到RGB空間，得到最終增強(qiáng)后的圖像。實(shí)驗(yàn)表明，該方法能更好地估計(jì)照度圖像，從而使增強(qiáng)后的圖像視覺效果獲得明顯改善。

低照度圖像；Retinex；引導(dǎo)濾波；保邊平滑；多尺度

在夜晚、黃昏等低照度條件下拍攝的圖像，其亮度、對比度較低，多含有噪聲，有些還具有很高的動態(tài)范圍，給后續(xù)圖像處理與分析帶來困難，對低照度圖像增強(qiáng)是目前解決這個(gè)問題的重要途徑。

低照度圖像增強(qiáng)的傳統(tǒng)方法主要有基于空域的增強(qiáng)算法和基于頻域的增強(qiáng)算法，如文獻(xiàn)[1-2]提出的基于灰度水平來實(shí)現(xiàn)圖像的增強(qiáng)算法。其可提高亮度與對比度，但往往會丟失信息，而且對于高動態(tài)范圍的低照度圖像，也不能壓縮動態(tài)范圍。LAND和MCCANN[3]提出的基于Retinex理論的圖像增強(qiáng)方法具有動態(tài)范圍壓縮大和顏色恒常性等特點(diǎn)。在Retinex理論中，圖像表示為照度圖像與反射圖像的乘積，照度圖像反映的是光照的情況，而反射圖像才是圖像的本征屬性。由于照度圖像對應(yīng)于圖像的低頻部分，往往通過將圖像卷積一個(gè)低通濾波器來求取照度圖像。

最經(jīng)典的Retinex算法是美國國家航空航天局的JOBSON等[4]提出的單尺度Retinex(single scale retinex，SSR)算法和多尺度Retinex(multi-scale retinex，MSR)算法[5]，這兩種算法以高斯濾波器作為低通濾波器，但增強(qiáng)后的圖像邊界處有較嚴(yán)重的光暈現(xiàn)象。文獻(xiàn)[6-7]先后提出了基于雙邊濾波的Retinex算法，雖然在保持邊緣、克服光暈上有了長足的進(jìn)步，但未徹底解決光暈的問題。文獻(xiàn)[8]提出一種局部線性濾波器——引導(dǎo)濾波器，其具有良好的邊緣保持性。為此，文獻(xiàn)[9-10]也提出基于引導(dǎo)濾波的Retinex算法，能夠較好地保留邊緣細(xì)節(jié)，并且效率很高。但引導(dǎo)濾波器在維持保邊能力的同時(shí)不能對圖像充分平滑，所以得不到對照度圖像的準(zhǔn)確估計(jì)，因此，增強(qiáng)后的圖像往往細(xì)節(jié)恢復(fù)不好。針對此，LI等[11]提出了一種加權(quán)引導(dǎo)濾波，以自適應(yīng)的因子作為權(quán)重來提高引導(dǎo)濾波的保邊平滑能力，可以在較充分的平滑同時(shí)不產(chǎn)生邊緣模糊，但是會有增加噪聲的問題出現(xiàn)。WANG等[12]提出將高斯平滑后的圖像作為引導(dǎo)濾波器的引導(dǎo)圖像，雖然能更好地平滑，但得到的增強(qiáng)圖像仍存在光暈。

影響引導(dǎo)濾波器保邊平滑能力的因素是圖像中存在各種對比度的細(xì)節(jié)。為了解決以上問題，本文提出一種基于迭代多尺度引導(dǎo)濾波的Retinex算法，通過迭代地增加引導(dǎo)濾波的尺度，在保邊的同時(shí)平滑各種對比度細(xì)節(jié)，能夠得到更好的照度圖像，從而有效地實(shí)現(xiàn)低照度圖像的增強(qiáng)。

1 Retinex理論

Retinex理論的基本思想是入射光決定了一幅圖像中的所有像素點(diǎn)的動態(tài)范圍的大小，而物體自身所固有不變的反射系數(shù)決定了圖像的內(nèi)在固有屬性。所以將一幅圖像分為照度圖像與反射圖像兩部分，通過去除或降低照度圖像對反射圖像的影響，就能得到圖像的本質(zhì)視覺感知，從而達(dá)到增強(qiáng)圖像的目的。根據(jù)Retinex理論，可以得到

其中，為輸入圖像；為反射圖像；為照度圖像。

為了符合人眼在獲取信息的非線性特性與簡化運(yùn)算，將式(1)改為對數(shù)形式

其中，(,)=log2(,),(,)=log2(,),(,)= log2(,)。

反射圖像對應(yīng)于圖像的高頻部分，照度圖像對應(yīng)于圖像的低頻部分。因此，求取照度圖像就是圖像與一個(gè)低通濾波器進(jìn)行卷積的結(jié)果。在Retinex理論中，光照是分片平滑的，所以光照區(qū)域分的越準(zhǔn)，區(qū)域內(nèi)細(xì)節(jié)平滑的越充分，照度圖像就會估計(jì)的越準(zhǔn)確，反射圖像也就越能反映圖像的本質(zhì)。因此，低通濾波器的選擇對于Retinex算法至關(guān)重要。SSR算法是最經(jīng)典的基于中心環(huán)繞的Retinex算法，其可表示為

其中，?表示卷積符號；表示高斯函數(shù)，即

為了兼顧顏色保真與細(xì)節(jié)保持，又發(fā)展了MSR算法。MSR算法是通過將不同尺度的SSR增強(qiáng)結(jié)果以線性加權(quán)的方法結(jié)合起來，同時(shí)實(shí)現(xiàn)細(xì)節(jié)增強(qiáng)和色彩保真?？杀硎緸?/p>

其中，G(,)為高斯函數(shù)；為尺度的總個(gè)數(shù)；W為第個(gè)尺度SSR所占的權(quán)重，且∑W=1。一般，MSR可選擇小、中和大3個(gè)尺度，令W=1/3，通過結(jié)合3個(gè)尺度的特性來達(dá)到壓縮動態(tài)范圍、增強(qiáng)細(xì)節(jié)和恢復(fù)色彩的目的。

SSR與MSR以高斯濾波器作為估計(jì)照度圖像的低通濾波器，缺乏保邊平滑能力，會導(dǎo)致光照分片不準(zhǔn)。如圖1所示，增強(qiáng)后的圖像雖然亮度與對比度有所提高，但存在顏色失真與光暈偽影的現(xiàn)象，使得視覺效果非常差。其中光暈偽影是因?yàn)樵诠烙?jì)照度圖像時(shí)，在高對比度的邊緣處發(fā)生了平滑，導(dǎo)致產(chǎn)生了光暈偽影現(xiàn)象。

圖1 SSR與MSR增強(qiáng)結(jié)果

2 照度圖像估計(jì)

針對傳統(tǒng)Retinex算法的不足，本文采用迭代的多尺度引導(dǎo)濾波代替高斯濾波，以便更好地估計(jì)照度圖像。

2.1 引導(dǎo)濾波及其性能分析

引導(dǎo)濾波采用最小二乘法的思想，在引導(dǎo)圖像和輸出圖像之間是一個(gè)局部線性濾波，如下式

其中，ω為尺寸為(2+1)×(2+1)的方形濾波窗口，為ω的半徑，為ω的中心像素；U為圖像中的像素；q為U線性映射到圖像中的像素；a與b為線性關(guān)系的系數(shù)，可由下式求出

其中，(a,b)為窗口ω的代價(jià)函數(shù)；為原圖像；為一個(gè)正則化參數(shù)，由式(7)解出

其中，μ和σ分別為引導(dǎo)圖像在窗口ω中的均值和標(biāo)準(zhǔn)差，||為ω中像素的個(gè)數(shù)。

與高斯濾波類似，引導(dǎo)濾波也具有核權(quán)重，考慮到一個(gè)像素點(diǎn)與所有覆蓋像素的重疊窗口ω相關(guān)，于是對所有情況取平均，核權(quán)重如下式表示[8]

那么，式(6)用核權(quán)重可表示為

下面對引導(dǎo)濾波器的保邊平滑性能進(jìn)行分析。當(dāng)引導(dǎo)圖像為原圖像時(shí)，濾波后圖像會保留原圖像的邊緣梯度，因此引導(dǎo)濾波表現(xiàn)出保邊平滑的特性。當(dāng)σ2>>時(shí)，該區(qū)域?yàn)楦叻讲顓^(qū)域，區(qū)域中像素U與另一個(gè)像素U有較高的對比度，U–μ與U–μ異號，則核權(quán)重W()接近于0，U被認(rèn)為是邊緣上的像素，而得到a≈1，b≈0，因此q≈U，邊緣得到了保留；當(dāng)σ2<，窗口ω中像素U與U相差較小，則W()為1，U被視為細(xì)節(jié)區(qū)域內(nèi)的像素，而得到a=0，b=μ，因此q=μ，細(xì)節(jié)得到了平滑。所以，對于引導(dǎo)濾波，可被視為尺度因子，決定了濾波器的大小。可被視為平滑因子，決定了平滑的程度。

為進(jìn)一步分析與對濾波效果的影響，可通過改變其值，以平滑后圖像的平均梯度作為平滑效果的評價(jià)。平均梯度表示圖像的邊界附近灰度的差異，反映了圖像微小細(xì)節(jié)反差變化的速率，表征了圖像的相對清晰程度。如下式

其中，與為圖像尺寸；ΔI與ΔI為圖像中像素的水平方向與垂直方向的一階導(dǎo)數(shù)。平均梯度G越小，說明平滑效果越充分。

對圖像Tree計(jì)算不同尺度參數(shù)與平滑參數(shù)下平滑后的平均梯度G，圖2表示在一定平滑參數(shù)、不同尺度參數(shù)下的G，隨著尺度變大，G呈現(xiàn)出負(fù)指數(shù)減少的趨勢，相對的，平滑效果則為對數(shù)增加的趨勢。圖3表示在一定尺度參數(shù)、不同平滑參數(shù)下的G，隨著平滑參數(shù)變大，G也呈現(xiàn)出負(fù)指數(shù)減少的趨勢，那么平滑效果則也為對數(shù)增加的趨勢。

圖2 不同尺度參數(shù)下的Gr

圖3 不同平滑參數(shù)下的Gr

雖然引導(dǎo)濾波具有較好的保邊平滑能力，但仍有不足之處?；赟UBR等[13]的分析，圖像中存在各種對比度的紋理細(xì)節(jié)與邊緣。對滿足邊緣為高對比度、細(xì)節(jié)為低對比度的圖像，引導(dǎo)濾波能夠較好地實(shí)現(xiàn)保邊平滑。但對于低照度圖像，往往存在具有較高對比度的細(xì)節(jié)，而這些細(xì)節(jié)會被誤認(rèn)為邊緣，因此在求照度圖像會被保留下來，于是無法估計(jì)出正確的照度圖像，也無法實(shí)現(xiàn)低照度圖像的有效增強(qiáng)。

2.2 迭代多尺度引導(dǎo)濾波

尺度是衡量不同對比度的重要性質(zhì)，尺度變大，一些低對比度的細(xì)節(jié)就會消失，因此可以在不同尺度下平滑不同對比度的細(xì)節(jié)。由引導(dǎo)濾波求取照度圖像如下

其中，為低照度圖像；為照度圖像?？芍?，濾波結(jié)果與不同尺度窗口有關(guān)。

當(dāng)較小時(shí)，引導(dǎo)濾波對低對比度細(xì)節(jié)區(qū)域的像素可以較好地平滑。而高對比度細(xì)節(jié)區(qū)域的像素之間相差較大，當(dāng)σ2>時(shí)，這些像素被誤認(rèn)為邊緣，很難被平滑。只有很大時(shí)，才能讓核權(quán)重W()接近1，但很大的不滿足式(7)代價(jià)函數(shù)的最小化，所以小尺度無法實(shí)現(xiàn)高對比度細(xì)節(jié)的平滑。而當(dāng)變大時(shí)，高對比度細(xì)節(jié)區(qū)域的像素方差σ2減小，在一定大小的下，就可以讓核權(quán)重W()接近于1，從而像素得到平滑。所以大尺度能夠在一定程度上實(shí)現(xiàn)高對比度細(xì)節(jié)的平滑，但相對較大時(shí)，需付出性能降低的代價(jià)。如圖4所示，圖4(b)~(d)為在小、中、大3種尺度下引導(dǎo)濾波估計(jì)出的照度圖像的掃描線，不難看出，隨著尺度增加，更高對比度的細(xì)節(jié)得到了平滑。圖4(d)在大尺度下，邊緣發(fā)生了偏離，表明引導(dǎo)濾波的保邊性能有所下降。

為了解決單尺度下引導(dǎo)濾波不能很好的保邊平滑的問題，本文提出采用迭代多尺度引導(dǎo)濾波，即

其中，(n)為第步照度圖像；(n–1)為第–1步照度圖像；(n)為第步尺度因子；(n)為第步平滑因子；G表示引導(dǎo)濾波，由與決定。該式表述了迭代濾波的方法，即由第–1步照度圖像，進(jìn)行尺度因子為(n)、平滑因子為(n)的引導(dǎo)濾波，得到第步照度圖像。當(dāng)滿足一定條件時(shí)，迭代終止，從而獲得最終的照度圖像。其計(jì)算過程步驟如下：

(1) 初始化。設(shè)(0)為原始圖像，(0)的平均梯度記為G(0)，初始(1)取2，(1)取0.042，迭代終止閾值為，最大迭代次數(shù)為max，令=1；

(2) 根據(jù)式(13)計(jì)算(n)，轉(zhuǎn)步驟(3)；

(3) 根據(jù)式(11)計(jì)算(n)的平均梯度G(n)，計(jì)算平均梯度G(n)與G(n–1)之差ΔG，轉(zhuǎn)步驟(4)

(4) 迭代終止條件判斷。若滿足ΔG<或>max，則終止迭代，轉(zhuǎn)步驟(5)；否則令←+1，(n)=f((n–1))，(n)=f((n–1))，返回步驟(2)；

(5) 輸出。(n)為最終的照度圖像。

由圖2~3分析，隨與的增大，濾波呈現(xiàn)出對數(shù)增加的效果。為此，與在迭代中選擇指數(shù)增加可以加快迭代，這里取(n)=f((n–1))=(1)×2，(n)=f((n–1))=(1)×2。如圖5所示，當(dāng)閾值取10–3、最大迭代次數(shù)max取10時(shí)，對圖像Tree兩次迭代之間照度圖像的平均梯度差值ΔG隨迭代次數(shù)的變化關(guān)系?？梢钥闯?，在第6次迭代后，ΔG=0.0009<，就可以完成迭代濾波過程。

圖4 不同尺度濾波結(jié)果

圖5 迭代之間的ΔGr隨迭代次數(shù)的關(guān)系

迭代多尺度濾波的保邊平滑性能體現(xiàn)在隨著尺度參數(shù)與平滑參數(shù)迭代增加，在保持邊緣的同時(shí)，逐步平滑掉各種對比度的細(xì)節(jié)。對像素I與I，滿足

其中，W()(n)表示在第步迭代時(shí)的核權(quán)重；I(n–1)表示第–1步迭代時(shí)像素I的值，I(n)表示第步迭代時(shí)像素I的值。若I(0)與I(0)為原圖像高對比度細(xì)節(jié)區(qū)域的兩個(gè)像素，由于其有較高的對比度，在某個(gè)小尺度下核權(quán)重W()較小，不易被平滑，但小程度的平滑也會讓兩個(gè)像素趨向于均值μ。隨著尺度增大，像素I與I會逐漸被平滑，平滑的結(jié)果使二者中高灰度像素越來越小，低灰度像素越來越大，也就是二者的差值逐漸減小，即|I(n)–I(n)|< |I(n–1)–I(n–1)|。由式(9)可知，|I(n)–I(n)|減小，則核權(quán)重變大，即W()(n)>W()(n–1)，這時(shí)就可以更大程度地平滑。因此，尺度增大到一定程度，存在的最高對比度的細(xì)節(jié)也會被平滑掉。而在整個(gè)迭代過程中，邊緣像素核權(quán)重一直接近0，不會被平滑。所以，迭代結(jié)束時(shí)，各種對比度的細(xì)節(jié)均得到了平滑，而邊緣一直較好地保持著。

為驗(yàn)證本文估計(jì)照度圖像算法的有效性，選取引導(dǎo)濾波和近些年比較有效的保邊平滑濾波器如域變換濾波器(domain transfer filter，DT)[14]和局部拉普拉斯金字塔濾波器(local Laplacian filter，LLF)[15]做比較。分別對圖像Tree求照度圖像，計(jì)算其掃描線。如圖6所示，引導(dǎo)濾波、DT濾波器和LLF估計(jì)的照度圖像中還存在樹葉等高對比度的細(xì)節(jié)，而墻面等邊緣卻發(fā)生了模糊，本文算法則很好地估計(jì)出光照情況，即平滑掉各種細(xì)節(jié)，又較準(zhǔn)確地劃分了光照區(qū)域。從掃描線也可以看出，本文算法對邊緣有很強(qiáng)地跟蹤能力，而對各種對比度的細(xì)節(jié)又有很強(qiáng)地平滑能力。

圖6 幾種保邊平滑濾波器的比較結(jié)果

3 低照度圖像增強(qiáng)

3.1 顏色空間轉(zhuǎn)換

如果在R、G、B 3個(gè)通道分別進(jìn)行Retinex處理，再合成彩色圖像，往往會導(dǎo)致色彩的失真。為此，本文將待增強(qiáng)的圖像從RGB空間轉(zhuǎn)為YCbCr空間，Y是亮度，Cb是RGB輸入信號藍(lán)色部分與RGB信號亮度值之間的差異，Cr是RGB輸入信號紅色部分與RGB信號亮度值之間的差異。在進(jìn)行Retinex處理時(shí)，只處理Y分量，不處理Cb與Cr分量，就可以避免色彩的失真，實(shí)現(xiàn)對色彩的保護(hù)。

RGB空間轉(zhuǎn)換到Y(jié)CbCr空間如下

YCbCr空間轉(zhuǎn)換到RGB空間如下

3.2 反射分量校正

由迭代多尺度引導(dǎo)濾波算法估計(jì)出照度圖像(,)，然后根據(jù)式(3)即可得到反射圖像。反射圖像中包含大量圖像細(xì)節(jié)，需要增強(qiáng)其對比度，以獲得更好的視覺效果。利用下式實(shí)現(xiàn)全局對比度提高

其中，r(,)為歸一化后的反射圖像；final為增強(qiáng)后的反射圖像；為調(diào)整參數(shù)。

如圖7(a)表示隨著的變化的函數(shù)曲線，越小，拉伸越明顯，增強(qiáng)效果也就會越好；圖7(b)顯示了本文曲線與伽馬曲線的區(qū)別，伽馬曲線對低亮度區(qū)域拉伸過大，會導(dǎo)致該區(qū)域的細(xì)節(jié)過增強(qiáng)，而本文的曲線則拉伸比較柔和，既實(shí)現(xiàn)了細(xì)節(jié)的增強(qiáng)，又可以避免失真。

3.3 本文增強(qiáng)算法流程

綜上，本文提出的低照度圖像增強(qiáng)算法可分為：

(1)根據(jù)式(16)，將待增強(qiáng)的低照度圖像從RGB空間轉(zhuǎn)換到Y(jié)CbCr空間；

(2)提取Y分量作為初始圖像，根據(jù)式(13)進(jìn)行迭代多尺度引導(dǎo)濾波，得到照度圖像；

(3)根據(jù)式(2)，得到反射圖像；

(4)由式(18)對反射圖像進(jìn)行全局對比度校正，得到增強(qiáng)后的反射圖像final，即為Y分量的增強(qiáng)結(jié)果；

(5)由式(17)，可將處理后的Y分量與Cb、Cr分量轉(zhuǎn)化為R、G、B分量，合成為低照度圖像的最終增強(qiáng)結(jié)果final。

算法流程如圖8所示。

圖7 校正曲線效果

圖8 本文算法流程示意圖

4 實(shí)驗(yàn)結(jié)果及分析

選取多組低照度圖像進(jìn)行試驗(yàn)，所有的試驗(yàn)均在操作系統(tǒng)為Windows 7、CPU為2.6 GHz酷睿i5，內(nèi)存為4 G的PC機(jī)上運(yùn)行，編程環(huán)境為MatlabR2013a。

為驗(yàn)證本文算法的可行性，首先與經(jīng)典Retinex算法SSR與MSR算法做比較。SSR的參數(shù)為=60；MSR的參數(shù)為1=15，2=80，3=120；本文算法迭代終止閾值取10–3，最大迭代次數(shù)max取10。對圖像Boat、Flower兩幅圖像進(jìn)行處理，結(jié)果如圖9~10所示。可以看出，SSR與MSR增強(qiáng)后的結(jié)果，邊緣出現(xiàn)了嚴(yán)重的光暈偽影，還存在一定的顏色失真，破壞了視覺效果；相比之下，本文算法增強(qiáng)后的結(jié)果邊緣清晰，沒有光暈，且未出現(xiàn)顏色失真的現(xiàn)象，有效克服了SSR與MSR算法存在的問題。

圖9 不同算法對圖像Boat增強(qiáng)結(jié)果

圖10 不同算法對圖像Flower增強(qiáng)結(jié)果

為進(jìn)一步驗(yàn)證本文算法的有效性，將與SSR算法、MSR算法、文獻(xiàn)[1]算法、文獻(xiàn)[16]算法進(jìn)行比較，文獻(xiàn)[1]、[16]算法參數(shù)按原文中的參數(shù)。對圖像Chair、Car、Street 3幅圖像進(jìn)行處理，結(jié)果如圖11~ 13所示?？梢钥闯?，SSR算法與MSR算法增強(qiáng)后的結(jié)果有較嚴(yán)重的光暈問題，尤其對于圖Chair，在椅子邊緣處有強(qiáng)烈的光暈偽影，導(dǎo)致視覺效果很差。另外，SSR算法與MSR算法的顏色失真問題也很大程度上破壞了增強(qiáng)效果。而文獻(xiàn)[1]算法、文獻(xiàn)[16]算法以及本文算法與其相比，在抑制光暈與克服色彩失真上表現(xiàn)較好。但文獻(xiàn)[1]作為一種基于灰度直方圖的算法，不能有效壓縮動態(tài)范圍，如圖11(d)中，椅子表面的高光沒有得到抑制；圖12(d)中，車燈處出現(xiàn)過亮。文獻(xiàn)[16]的算法對信息細(xì)節(jié)恢復(fù)不夠，從增強(qiáng)結(jié)果的放大圖中可以看出，如車燈、廣告牌上的數(shù)字經(jīng)過處理后較為模糊。本文算法增強(qiáng)后的圖像細(xì)節(jié)明顯改善，對各種對比度的細(xì)節(jié)與邊緣都得到了較好的增強(qiáng)，清晰度較高。

圖11 不同算法對圖像Chair增強(qiáng)結(jié)果

圖12 不同算法對圖像Car增強(qiáng)結(jié)果

圖13 不同算法對圖像Street增強(qiáng)結(jié)果

通常情況下，對于低照度圖像增強(qiáng)算法的評價(jià)沒有真實(shí)圖像可做參考，屬于無參考客觀質(zhì)量評測范疇。本文主要在細(xì)節(jié)信息恢復(fù)、對比度、是否有光暈偽影這3個(gè)方面進(jìn)行評價(jià)。其中，對比度采用文獻(xiàn)[17]的CPBD(cumulative probability of blur detection)評價(jià)指標(biāo)，其值越大，代表圖像對比度越高；是否有光暈偽影采用文獻(xiàn)[18]提出的PQM (perceptual quality metric)指標(biāo)，其值越接近10，說明圖像存在的光暈越少，質(zhì)量也就越好；對于細(xì)節(jié)信息恢復(fù)程度，這里采用清晰度表示[19]

其中，max為圖像中像素最大的值。清晰度的值越大，代表圖像的細(xì)節(jié)與邊緣越清晰。

不同算法增強(qiáng)圖像的定量評價(jià)結(jié)果見表1?？梢钥闯?，與SSR與MSR算法相比，本文算法的客觀質(zhì)量評價(jià)結(jié)果有較大地提高，表明本文算法在性能上確實(shí)有一定的改善。與文獻(xiàn)[1]和[16]的算法相比，本文算法的清晰度與PQM值要高于文獻(xiàn)[1]和[16]的算法，說明本文算法在恢復(fù)細(xì)節(jié)、增強(qiáng)邊緣、克服光暈等方面取得了更好的效果；而CPBD值對于圖像Chair與文獻(xiàn)[1]算法相同，對于圖像Street略小于文獻(xiàn)[16]的算法，這是因?yàn)楸疚乃惴ú捎萌謱Ρ榷刃Ｕ栽趯Ρ榷确矫姹憩F(xiàn)還有待于提升。

表1 不同算法增強(qiáng)圖像的客觀評價(jià)結(jié)果

5 結(jié) 論

為解決低照度圖像的增強(qiáng)問題，本文在Retinex理論的基礎(chǔ)上提出一種改進(jìn)的Retinex算法。在求取照度圖像時(shí)利用迭代的多尺度引導(dǎo)濾波，很好地保留邊緣與平滑各種細(xì)節(jié)，可以得到較為準(zhǔn)確的照度圖像。以此照度圖像求得的反射圖像可以很好地增強(qiáng)圖像細(xì)節(jié)信息，提高清晰度。實(shí)驗(yàn)結(jié)果表明，本文算法對低照度圖像，能夠有效提高亮度，增強(qiáng)細(xì)節(jié)，提高視覺感知。但本算法增強(qiáng)后的圖像還存在“泛灰”的現(xiàn)象，在色彩與對比度校正方面仍有不足，需要在下一步工作中研究解決。

[1] SINGH K, KAPOOR R. Image enhancement using exposure based sub image histogram equalization[J]. Pattern Recognition Letters, 2014, 36(1): 10-14.

[2] LEE C, KIM C. Contrast enhancement based on layered difference representation[J]. IEEE Transactions on Image Processing, 2013, 22(12): 5372-5384.

[3] LAND E H, MCCANN J. Lightness and retinex theory[J]. Journal of the Optical Society of America, 1971, 61(1): 1-11.

[4] JOBSON D, RAHMAN Z, WOODEL G A. Properties and performance of a center/surround Retinex[J]. IEEE Transactions on Image Processing, 1997, 6(3): 451-462.

[5] JOBSON D, RAHMAN Z, WOODEL G A. Amulti-scale Retinex for bridging the gap between color and the human observation of scenes[J]. IEEE Transactions on Image Processing, 1997, 6(7): 965-976.

[6] ELADE M. Retinex by two bilateral filters[C]// Proceedings of the Scale-Space Conference 2005. Berlin: Springer Press, 2005: 217-229.

[7] 胡韋偉, 汪榮貴, 方帥, 等. 基于雙邊濾波的Retinex圖像增強(qiáng)算法[J]. 工程圖學(xué)學(xué)報(bào), 2010, 38(2): 217-223.

[8] HE K M, SUN J, TANG X. Guided image filtering[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2013, 35(6): 1397-1409.

[9] AHN H, KEUM B, KIM D, et al. Adaptive local tone mapping based on retinex for high dynamic range images[C]//IEEE International Conference on Consumer Electronics (ICCE). New York: IEEE Press, 2013: 153-156.

[10] 肖創(chuàng)柏, 趙宏宇, 禹晶. 基于引導(dǎo)濾波的Retinex快速夜間彩色圖像增強(qiáng)技術(shù)[J]. 北京工業(yè)大學(xué)學(xué)報(bào), 2013, 39(12): 1868-1873.

[11] LI Z, ZHENG J, ZHU Z, et al. Weighted guided image filtering[J]. IEEE Transactions on Image Processing, 2015, 24(1): 120-129.

[12] WANG Y, WANG H, YIN C, et al. Biologically inspired image enhancement based on Retinex[J]. Neurocomputing, 2015, 177: 373-384.

[13] SUBR K, SOLER C, DURAND F. Edge-preserving multi-scale image decomposition based on local extrema[J]. ACM Transactions on Graphics, 2009, 28(5): 147-154.

[14] GASTAL E S, OLIVEIRA M M. Domin transform for edge-aware image and video processing [J]. ACM Transactions on Graphics, 2011, 30(3): 69-77.

[15] AUBRY M, SYLVAIN R. Fast local Laplacian filters: theory and applications[J]. ACM Transactions on Graphics, 2014, 33(5): 167-175.

[16] GUO X J, LI Y, LING H B. LIME: low-light image enhancement via illumination map estimation[J]. IEEE Transactions on Image Processing, 2016, 26(2): 982-993.

[17] NARVEKAR N D, KARAM L J. A no-reference perceptual image sharpness metric based on a cumulative probability of blur detection[C]//International Workshop on Quality of Multimedia Experience. New York: IEEE Press, 2009: 87-91.

[18] GUPTA N, JHA R K. Enhancement of dark images using dynamics to chastic resonance with an is otropic diffusion[J]. Journal of Electronic Imaging, 2016, 25(2): 023017.

[19] FENG B B, TANG Y B, Zhou L, et al. Image enhancement under low luminancewith strong light weakening [C]//2016 8th International Conferenceon Wireless Communications & Signal Processing (WCSP). New York: IEEE Press, 2016: 1-5.

Low-Light Image Enhancement Based on Iterative Multi-Scale Guided Filter Retinex

ZHANG Jie1,2, ZHOU Pucheng1,2, ZHANG Qian1,2

(1. Army Academy of Artillery and Air Defence Forces, Hefei Anhui 230031, China; 2. Anhui Province Key Laboratory of Polarized Imaging Detection Technology, Hefei Anhui 230031, China)

In order to improve the visual perception of low-light image, a kind of low-light image enhancement algorithm based on iterative multi-scale guided filter Retinex was proposed. Firstly, the low-light image was converted into YCbCr color space. Secondly, the illumination image of the Y channel was estimated with iterative multi-scale guided filter. By increasing the scale parameter and the smoothing parameter iteratively in exponential form, the edge was preserved and details of different contrasts was eliminate. And then, the reflected image was obtained by the estimated illumination image and the contrast correction procedure. Finally, the image was converted into RGB color space, and this is the final enhanced image. Experimental results show that our method can better obtain the illumination image so that the visual effect of the enhanced image can be improved obviously.

low-light image; Retinex; guided filter; edge-preserving and smoothing; multi-scale

TP391

10.11996/JG.j.2095-302X.2018010001

A

2095-302X(2018)01-0001-11

2017-06-05；

2017-07-09

國家自然科學(xué)基金項(xiàng)目(61379105)

張 杰(1994–)，男，河北張家口人，碩士研究生。主要研究方向?yàn)閳D像分析與處理。E-mail：zj199402@163.com

周浦城(1977–)，男，江西宜春人，副教授，博士。主要研究方向?yàn)閳D像分析與處理、智能計(jì)算。E-mail：zhoupc@hit.edu.cn