李連志，邢 川

(河北科技大學(xué)，河北 石家莊 050000)

1 引言

圖像在成像過(guò)程中，會(huì)因自然環(huán)境因素或人為因素導(dǎo)致出現(xiàn)缺乏光照、視覺(jué)效果不佳以及圖像清晰度不高等情況，而陰影圖像、弱光圖像的處理一直是圖像處理中難點(diǎn)之一。解決上述問(wèn)題不僅具有重要的理論意義，而且在實(shí)際應(yīng)用中也有很大應(yīng)用前景，例如航空影像、遙感影像、視頻監(jiān)控中的烏云影像等。合理運(yùn)用圖像增強(qiáng)技術(shù)，可以從根本上解決圖像存在的質(zhì)量問(wèn)題，令圖像視覺(jué)效果更清楚，更能體現(xiàn)圖像的細(xì)節(jié)。

王建衛(wèi)[1]為了解決低照度彩色圖像邊界提取不完整的問(wèn)題，依據(jù)點(diǎn)處理和顏色空間變換理論，提出了應(yīng)用非線性函數(shù)變換方法處理HSI顏色模型的強(qiáng)度分量彩色圖像增強(qiáng)方法。設(shè)計(jì)非線性函數(shù)變換進(jìn)行彩色圖像增強(qiáng)計(jì)算，并分析了算法中RGB和HSI顏色模型的轉(zhuǎn)換強(qiáng)度，進(jìn)行分量計(jì)算，獲取增強(qiáng)后圖像。趙軍輝[2]等人提出基于Lab色彩空間和色調(diào)映射的Retinex圖像對(duì)比增強(qiáng)算法。首先，在Lab色彩空間中將一個(gè)低對(duì)比度的輸入圖像分解成亮度和色度分量，并使用自適應(yīng)雙邊濾波估計(jì)照明的強(qiáng)度，以便根據(jù)亮度和顏色值來(lái)考慮合適的相鄰像素。然后利用基于拋物線的色調(diào)映射函數(shù)來(lái)提高估計(jì)光照?qǐng)D像的對(duì)比度。最后，將加強(qiáng)的亮度和原始色度結(jié)合在一起，產(chǎn)生一個(gè)增強(qiáng)彩色輸出圖像。

雖然上述方法都能使圖像色彩增強(qiáng)，但是各色塊間像素分布依舊不夠均衡，人眼視覺(jué)舒適度較差，因此，為了能夠獲取出較好的視覺(jué)效果，將圖像色彩空間從RGB轉(zhuǎn)換為HSV，并對(duì)轉(zhuǎn)換后圖像進(jìn)行亮度分量增強(qiáng)；根據(jù)需要處理圖像的構(gòu)成因素，利用壓縮低頻、提升高頻法，獲取出經(jīng)同態(tài)濾波處理后的視覺(jué)效果圖像，構(gòu)建巴特沃思型的同態(tài)濾波器；因?yàn)樵谟?jì)算同態(tài)濾波過(guò)程中需要同時(shí)在不相同的圖像塊中單獨(dú)運(yùn)行，所以導(dǎo)致圖像邊緣處會(huì)形成不同樣式的塊效應(yīng)，并會(huì)有圖像像素越平滑，塊效應(yīng)越顯著的情況，因此，提出基于重疊塊的塊效應(yīng)去除方法，實(shí)現(xiàn)平面視覺(jué)圖像色彩增強(qiáng)。

2 空間轉(zhuǎn)換及亮度增強(qiáng)

2.1 色彩空間轉(zhuǎn)換

RGB顏色模型由R、G和B原色組成，雖然RGB模型有一個(gè)舒適的色彩顯示效果，但是模型中兩個(gè)相似顏色值之間的差異太大，導(dǎo)致計(jì)算結(jié)果不完全一致，修改后模型三個(gè)分量之間的相關(guān)性很高。如果其中一個(gè)組件值更改，則其它兩個(gè)組件值也更改。與RGB模型不同，HSV(Hue， Saturation， Value)是一種基于人眼視覺(jué)系統(tǒng)的顏色模型。HSV模型主要有兩點(diǎn)明顯的特征，其一是V分量的取值變化和圖像整體的色彩信息沒(méi)有關(guān)聯(lián)；其二是H或S分量與人眼視覺(jué)所感受顏色的方法是一樣的。根據(jù)這兩點(diǎn)特征，即可進(jìn)一步得知HSV色彩空間模型比RGB模型更適用于彩色圖像處理算法，因此，在處理圖像前，首先將需要處理圖像的色彩空間由RGB模式轉(zhuǎn)換為HSV模式[3]。

假設(shè)(r，g，b)為一個(gè)顏色的紅、綠、藍(lán)坐標(biāo)，要得到HSV空間中的(h，s，v)值，那么轉(zhuǎn)換公式就可以寫(xiě)為

υ=max{r，g，b}

(1)

(2)

(3)

當(dāng)給定(h，s，v)值，對(duì)應(yīng)的(r，g，b)三原色為

(r，g，b)=(r′+λ，g′+λ，b′+λ)

(4)

其中，λ=υ-c。

2.2 亮度分量增強(qiáng)

由于在運(yùn)用同態(tài)濾波算法對(duì)圖像進(jìn)行色彩增強(qiáng)處理之前，需要將圖像預(yù)處理，獲取出與原圖相比較亮的圖像，因此，利用全局直方圖和局部直方圖均衡化方法，令圖像灰度值達(dá)到分布均勻效果，然后運(yùn)用灰度映射擴(kuò)展圖像中像素值的動(dòng)態(tài)范圍。但預(yù)處理過(guò)程中，容易發(fā)生顏色失真，產(chǎn)生新噪聲的情況，并具有難以達(dá)到滿意效果、計(jì)算量過(guò)多以及形成馬賽克效應(yīng)等缺點(diǎn)[4]。

為了解決上述問(wèn)題，將圖像分塊，然后將基于DCT變換的同態(tài)濾波器分為所使用的每個(gè)子圖像[5]。將待處理的圖像分為8×8個(gè)子圖像，圖像塊為8×8二維矩陣f(m，n)，其中0≤m≤7，0≤n≤7，則f(m，n)的DCT變換定義為

(5)

(6)

其中，0≤u≤7，0≤t≤7。

(7)

經(jīng)過(guò)色彩空間轉(zhuǎn)換和亮度分量增強(qiáng)可以提取出較亮的圖像，使圖像灰度值的分布更加均勻。

3 基于DCT變換的同態(tài)濾波

一般情況下常用的同態(tài)濾波算法有高斯型與指數(shù)型，這兩種濾波算法的不同之處在于處理圖像時(shí)傳遞函數(shù)的取值，根據(jù)計(jì)算選取出適用于解決圖像問(wèn)題的傳遞函數(shù)，就可以直接讓圖像達(dá)到最滿意的效果，如果不是最合適的函數(shù)，那么便可獲取出不同的增強(qiáng)效果，直至經(jīng)過(guò)調(diào)整令其達(dá)到最佳效果。然而在具體計(jì)算過(guò)程中，傳遞函數(shù)的控制參數(shù)太多，需要經(jīng)過(guò)大量計(jì)算來(lái)確定最終可以達(dá)到滿意效果的傳遞函數(shù)。因此，實(shí)用性不強(qiáng)。

3.1 構(gòu)建同態(tài)濾波模型

鑒于現(xiàn)有方法存在的問(wèn)題，根據(jù)構(gòu)成圖像的光學(xué)特征得知，圖像由自然光照度分量i(x，y)和圖中主要目標(biāo)的反射分量r(x，y)構(gòu)成，根據(jù)其組成條件[6]，數(shù)學(xué)模型可以表示為

f(x，y)=i(x，y)·r(x，y)

(8)

一般情況下，照度分量i(x，y)可以充分反射出光照條件是好或是壞，如果照度分量呈緩慢變化，那么頻譜就會(huì)落在低頻區(qū)域；而反射分量r(x，y)則是主要映射圖像目標(biāo)物件的內(nèi)容細(xì)節(jié)，該分量圖像中的頻譜區(qū)域?yàn)楦哳l部分[7]。根據(jù)上述分布情況，只需要在f(x，y)中將照度分量i(x，y)和反射分量r(x，y)劃分開(kāi)，并在此基礎(chǔ)上采用壓縮圖像低頻，提高圖像高頻的方法，就可以實(shí)現(xiàn)降低圖像中亮度分量，使圖像視覺(jué)效果更清晰。

針對(duì)式(8)左右兩側(cè)，計(jì)算選取出具有代表性的對(duì)數(shù)，這樣就可以將乘法運(yùn)算轉(zhuǎn)換為加減運(yùn)算組合，便有

lnf(x，y)=lni(x，y)+lnr(x，y)

(9)

根據(jù)上式計(jì)算結(jié)果，利用DCT變換代替?zhèn)鹘y(tǒng)同態(tài)濾波中的傅里葉變換[8]，經(jīng)過(guò)對(duì)式(9)變換后即可得出下式

F(u，t)=I(u，t)+R(u，t)

(10)

其中

F(u，t)=DCT[lnf(x，y)]

I(u，t)=DCT[lni(x，y)]

R(u，t)=DCT[lnr(x，y)]

(11)

用同態(tài)濾波函數(shù)H(u，t)來(lái)處理上式中的F(u，t)，得出

S(u，t)=H(u，t)F(u，t)=H(u，t)I(u，t)+H(u，t)R(u，t)

(12)

經(jīng)濾波處理后，再經(jīng)過(guò)逆DCT變換(IDCT)，得到

s(x，y)=IDCT[H(u，t)F(u，t)]

=IDCT[H(u，t)I(u，t)]+IDCT[H(u，t)R(u，t)]

(13)

經(jīng)過(guò)DCT變換計(jì)算后，針對(duì)式(13)進(jìn)行指數(shù)變換，這樣即可獲取出濾波后的圖像

g(x，y)=es(x，y)

(14)

通過(guò)DCT變換運(yùn)算的均勻?yàn)V波算法，選取圖像展開(kāi)的對(duì)數(shù)，得出對(duì)數(shù)結(jié)果后實(shí)施DCT變換處理，其次，在變換域中對(duì)適合改進(jìn)的圖像進(jìn)行濾波操作，可得出基于反向和指數(shù)變換的同態(tài)濾波器。

3.2 同態(tài)濾波器

在同態(tài)濾波計(jì)算過(guò)程中，同態(tài)濾波函數(shù)H(u，t)的取值結(jié)果會(huì)直接影響最后變換計(jì)算的結(jié)果是否精準(zhǔn)。為了能從根本上緩解圖像因光照環(huán)境導(dǎo)致的圖像色彩不均勻或是圖像視覺(jué)效果失真等情況，應(yīng)該降低I(u，t)頻率分量，同理為了更加凸顯圖像細(xì)節(jié)，提高圖像中目標(biāo)物的對(duì)比度，則需要提高R(u，t)頻率分量。通過(guò)上述分析，再融合同態(tài)濾波函數(shù)的特點(diǎn)，設(shè)計(jì)巴特沃思型的同態(tài)濾波器，其中表達(dá)式為

(15)

其中，γL<1，γH>1；L代表濾波器的平滑階數(shù)；D0代表停止過(guò)濾的頻率；D(u，t)代表點(diǎn)(u，t)到DCT變換最初點(diǎn)的距離，其計(jì)算公式如下

(16)

4 基于重疊塊的塊效應(yīng)去除方法

4.1 基于HVS的塊效應(yīng)可見(jiàn)度函數(shù)

在使用消除塊原理的過(guò)程中，需結(jié)合視覺(jué)系統(tǒng)特征。曾有相關(guān)研究表明：部分空間活動(dòng)及圖像光亮程度可遮掩消除塊效應(yīng)，但在紋理部分，區(qū)塊效應(yīng)的可見(jiàn)度略低。因此，在局部背景亮度較高的區(qū)域，塊效應(yīng)的可見(jiàn)度相對(duì)較低。如圖1所示，塊與相鄰塊之間的塊邊界被理解為塊邊界的局部區(qū)域[9]。圖1為水平相鄰塊組成的新塊。

圖1 水平相鄰塊組成新塊示意圖

塊效應(yīng)可以通過(guò)垂直和局部空間活動(dòng)來(lái)掩蓋，基于同態(tài)濾波計(jì)算結(jié)果定義水平和垂直函數(shù)

(17)

(18)

其中，將Ah和Av描述為塊c在水平或是垂直方向的特征活動(dòng)性，這樣根據(jù)其特性即可令R(u，v)代表塊c模型中反映局部活動(dòng)性的剩余塊在DCT域中的取值。

一般情況下，對(duì)于存在水平方向中的塊效應(yīng)來(lái)說(shuō)，垂直方向的活動(dòng)是導(dǎo)致塊效應(yīng)隱蔽的主要原因，因此，可根據(jù)下式判定垂直或水平方向塊效應(yīng)的所有活動(dòng)性

(19)

(20)

因此，定義了空間頻率活動(dòng)對(duì)水平方向阻滯效應(yīng)的掩蔽函數(shù)Mh，見(jiàn)式(21)

(21)

根據(jù)上式計(jì)算結(jié)果，即可寫(xiě)出亮度均值的計(jì)算式

(22)

根據(jù)上述計(jì)算結(jié)果得知，如果圖像空間活動(dòng)性的頻率不能達(dá)到掩蓋函數(shù)的正常取值范圍，那么圖像中的方塊效應(yīng)也會(huì)隨之產(chǎn)生可見(jiàn)度變化，掩蓋函數(shù)取值越小，塊效應(yīng)的可見(jiàn)度就會(huì)呈現(xiàn)出越低的情況，這樣就可以進(jìn)一步判定，垂直方向的塊效應(yīng)可見(jiàn)度函數(shù)ηv

(23)

4.2 圖像色彩增強(qiáng)的塊效應(yīng)去除

當(dāng)塊c的ηh大于等于閾值T時(shí)，對(duì)塊c進(jìn)行邊緣檢測(cè)，邊緣檢測(cè)的目的是防止邊緣被視為塊效應(yīng)。為了能夠有效防止出現(xiàn)忽略圖像邊緣像素的情況，導(dǎo)致在后期處理過(guò)程中圖像呈現(xiàn)出視覺(jué)效果不清晰的情況，提出一種新的圖像邊緣檢測(cè)方法，并且當(dāng)像素轉(zhuǎn)換塊c滿足于下列計(jì)算條件時(shí)，就有

c(i，j)-c(i，j+1)≤c(i，4)-c(i，5)

(24)

c(i，4)-c(i，5)≤2QP

(25)

其中，i=0，1…，7，j=0，1，…，6。

根據(jù)下述兩點(diǎn)條件，可以認(rèn)為轉(zhuǎn)換塊中沒(méi)有邊。

條件一：為了避免轉(zhuǎn)換過(guò)程中出現(xiàn)較大的效應(yīng)誤差，利用空間頻率結(jié)果將像素之間的差異控制在塊邊界像素可接受范圍內(nèi)[10]。

條件二：防止邊緣僅出現(xiàn)在塊邊界上。

為了消除平滑塊的塊效應(yīng)，用線性函數(shù)代替階躍函數(shù)。如果用線性函數(shù)代替梯形函數(shù)，那么這些系數(shù)的取值在本質(zhì)上就會(huì)發(fā)生變化，但只運(yùn)用部分轉(zhuǎn)換塊信息的話，就會(huì)忽略圖像塊中左右兩側(cè)部分的信息，進(jìn)而形成新的塊效應(yīng)。通過(guò)塊a和塊b左右局部信息，改變C(0，1)，C(0，3)C(0，5)和C(0，7)的值，從而改變?chǔ)碌闹?，最后再用線性函數(shù)代替階躍函數(shù)去除塊效應(yīng)。從而完成對(duì)平面視覺(jué)圖像整體色彩的增強(qiáng)。

5 仿真研究

為了驗(yàn)證所提算法的有效性，運(yùn)用MATLAB軟件進(jìn)行仿真，并將所提算法的仿真結(jié)果與文獻(xiàn)[1]、[2]方法進(jìn)行比較。

5.1 主觀質(zhì)量評(píng)價(jià)

在實(shí)驗(yàn)過(guò)程中，將塊效應(yīng)可見(jiàn)性函數(shù)的閾值設(shè)置為0.02。圖2(a)代表原始圖像，圖2(b)代表塊效應(yīng)嚴(yán)重的JPEG解碼圖像，圖2(c)代表本文算法處理的重建圖像。

圖2 效果圖像圖

分析圖2可知，從視覺(jué)角度出發(fā)，運(yùn)用所提算法增強(qiáng)后的效果圖，不論是從圖像色彩還是圖像塊效應(yīng)方面均得到了有效改善，圖2(b)中圖像具有明顯的塊效應(yīng)，圖像模糊、邊緣信息嚴(yán)重丟失，不能達(dá)到更清晰的效果。這是由于所提算法為了能夠有效防止出現(xiàn)忽略圖像邊緣像素導(dǎo)致形成塊效應(yīng)的問(wèn)題，提出了一種圖像邊緣檢測(cè)方法，該方法可以將圖像像素之間的差異控制在塊邊界像素可接受范圍內(nèi)，從而有效抑制塊效應(yīng)的產(chǎn)生。

5.2 客觀質(zhì)量評(píng)價(jià)

為了客觀地測(cè)量和評(píng)價(jià)算法性能，采用信噪比(SNR)對(duì)處理后的圖像質(zhì)量進(jìn)行測(cè)量，結(jié)果如表1所示。

表1 不同比特率下的信噪比

從表中可以看出，所提方法在不同比特率下的信噪比均高于文獻(xiàn)[1]方法和文獻(xiàn)[2]方法，說(shuō)明提出的算法可以有效地消除圖像的塊效應(yīng)，保護(hù)圖像邊緣信息，增強(qiáng)了圖像中因光照環(huán)境影響導(dǎo)致的圖像色彩。

為了進(jìn)一步驗(yàn)證所提方法的圖像色彩增強(qiáng)效果，建立圖像客觀質(zhì)量綜合評(píng)價(jià)準(zhǔn)則，運(yùn)用該準(zhǔn)則評(píng)價(jià)不同方法，將綜合評(píng)價(jià)準(zhǔn)則定義為CAC(Comprehensive Assessment Criteria Index)

CAC=EmNnCp

(26)

其中，Em表示圖像結(jié)構(gòu)相似性測(cè)度；Nn表示歸一化灰度差；Cp表示色彩加權(quán)還原度。CAC取值越大，說(shuō)明圖像色彩增強(qiáng)效果更好，不同方法的CAC值對(duì)比結(jié)果見(jiàn)表2。

表2 不同方法CAC值對(duì)比

分析表2可知，采用所提方法對(duì)不同圖像色彩進(jìn)行增強(qiáng)之后，所得的CAC值均高于1，該方法的CAC值明顯高于文獻(xiàn)[1]方法和文獻(xiàn)[2]方法，說(shuō)明所提方法的圖像增強(qiáng)效果更好，進(jìn)一步驗(yàn)證了該方法的優(yōu)勢(shì)性。

以上實(shí)驗(yàn)結(jié)果說(shuō)明同態(tài)濾波技術(shù)是在該研究領(lǐng)域中較為有效的色彩光照?qǐng)D像增強(qiáng)法，該方法可以避免圖像失真的問(wèn)題，并且可以做到在保持原有圖像色彩的基礎(chǔ)上，增強(qiáng)色彩暗區(qū)的細(xì)節(jié)，解決因光照環(huán)境導(dǎo)致的色彩分布不均或是過(guò)暗的問(wèn)題。

6 結(jié)論

圖像成像過(guò)程中，因自然或人為因素導(dǎo)致的圖像失真或色彩變暗會(huì)影響視覺(jué)效果，該問(wèn)題已經(jīng)成為該研究領(lǐng)域中需要解決的主要問(wèn)題之一，因此，提出基于同態(tài)濾波的平面視覺(jué)圖像色彩增強(qiáng)算法。根據(jù)色彩空間特征，將RGB轉(zhuǎn)換為HSV形式，以便于得出效果最佳的處理圖像效果，經(jīng)過(guò)同態(tài)濾波計(jì)算構(gòu)建出巴特沃思型的同態(tài)濾波器，利用可見(jiàn)度函數(shù)對(duì)圖像中存有的塊效應(yīng)去除，獲取出色彩增強(qiáng)后的圖像。雖然本文算法最后可以實(shí)現(xiàn)圖像色彩增強(qiáng)目的，但由于去除塊效應(yīng)計(jì)算中，可見(jiàn)度函數(shù)結(jié)果閾值需要精準(zhǔn)且有效，導(dǎo)致需要針對(duì)閾值結(jié)果反復(fù)，下一步將主要針對(duì)計(jì)算過(guò)程進(jìn)行優(yōu)化。