劉黎明

(中國(guó)船舶重工集團(tuán)公司第七二三研究所，江蘇 揚(yáng)州 225001)

多尺度Retinex圖像增強(qiáng)應(yīng)用中的彩色空間選擇

劉黎明

(中國(guó)船舶重工集團(tuán)公司第七二三研究所，江蘇 揚(yáng)州 225001)

研究了基于Retinex理論的圖像增強(qiáng)算法，介紹了圖像和視頻處理中常用的RGB、HSV、YIQ、YCbCr和Lab5種顏色空間，給出了均方誤差(MSE)、峰值信噪比(PSNR)、信息熵、彩色增強(qiáng)因子(CEF)和結(jié)構(gòu)相似度(SSIM)5個(gè)圖像質(zhì)量評(píng)價(jià)指標(biāo)的計(jì)算公式，并用這5個(gè)指標(biāo)評(píng)價(jià)了CLAHE算法在5種顏色空間上的增強(qiáng)效果。實(shí)驗(yàn)結(jié)果表明，相比于其它4種顏色空間，多尺度Retinex算法在HSV顏色空間上取得了最好的增強(qiáng)效果。

圖像增強(qiáng)；Retinex理論；多尺度；顏色空間

0 引 言

Retinex理論是Land等人在1971年提出的一種關(guān)于人類視覺(jué)的光亮度和色彩感知模型[1]，Retinex是由Retina (視網(wǎng)膜)和Cortex (腦皮層)2個(gè)單詞組成，全稱是視網(wǎng)膜腦皮層理論。該理論認(rèn)為人類在感知物體表面顏色時(shí)，與物體表面的反射特性密切相關(guān)，而與投射到人眼的光譜特性關(guān)系不大，也就是說(shuō)人眼感知物體的表面顏色不受光照變化的影響，具有色彩常恒性。目前，Retinex理論已廣泛應(yīng)用于遙感、醫(yī)學(xué)、紅外、可見(jiàn)光等圖像的對(duì)比度增強(qiáng)處理和去霧處理中[2-5]。

經(jīng)學(xué)者們的不斷研究，現(xiàn)已出現(xiàn)多種不同形式的Retinex算法。第1個(gè)Retinex算法是Land等人提出的一種隨機(jī)路徑算法[1]，其基本思路是比較像素的選擇,是從當(dāng)前圖像像素的鄰域中隨機(jī)產(chǎn)生的，該算法缺陷在于參數(shù)不易確定。此后，學(xué)者們分別提出了基于泊松方程的Retinex算法[6]、McCann's Retinex算法[7]、單尺度Retinex算法[8]、多尺度Retinex算法[9]、帶色彩恢復(fù)的多尺度Retinex算法[10]、變分Retinex算法[11]等。其中，多尺度Retinex算法具有較大的動(dòng)態(tài)壓縮范圍和較高的色彩保真度，故最引人關(guān)注。不過(guò)，大部分學(xué)者關(guān)注的是如何選擇多尺度Retinex算法中的參數(shù)，包括尺度個(gè)數(shù)、尺度大小和權(quán)重，而少有人研究多尺度Retinex算法在什么顏色空間中能獲得更好的增強(qiáng)效果。為此，本文的研究目標(biāo)就是通過(guò)某些準(zhǔn)則來(lái)客觀評(píng)估多尺度Retinex算法在各種彩色空間中的圖像增強(qiáng)效果。

1 Retinex算法

本節(jié)將介紹Retinex理論基礎(chǔ)、單尺度Retinex和多尺度Retinex。

1.1 Retinex理論基礎(chǔ)

根據(jù)Land 提出的Retinex 模型，原始圖像S(x，y)可分解為光照?qǐng)D像L(x，y)與反射圖像R(x，y)的乘積，數(shù)學(xué)表達(dá)式為：

(1)

原始圖像S(x，y)是人眼觀察到的圖像，光照?qǐng)D像L(x，y)決定了圖像的動(dòng)態(tài)范圍，而反射圖像R(x，y)包含了圖像中物體的細(xì)節(jié)特征。因此，在原始圖像中去除或降低照射圖像的影響，并保留能反映物體本質(zhì)特征的反射屬性是Retinex理論的基本思想。針對(duì)公式(1)而言，為了將復(fù)雜的乘積形式變成簡(jiǎn)單的加減運(yùn)算，可將公式(1)放入對(duì)數(shù)域處理，即有：

(2)

式中：s(x,y)=ln(S(x,y))；l(x,y)=ln(L(x,y))；r(x,y)=ln(R(x,y))。

Retinex算法流程如圖1所示。

1.2 單尺度Retinex(SSR)

SSR在二維路徑的前提下，選取高斯函數(shù)作為卷積核函數(shù)，估計(jì)出光照分量，再用對(duì)數(shù)運(yùn)算去除光照分量。SSR的數(shù)學(xué)表達(dá)式為：

r(x,y)=ln(R(x,y))=ln(S(x,y))- ln(S(x,y)*G(x,y))

(3)

式中：r(x,y)為反射分量；S(x,y)*G(x,y)為光照分量的估計(jì),符號(hào)“*”表示卷積運(yùn)算,G(x,y)為高斯環(huán)繞函數(shù)，相當(dāng)于一個(gè)高斯模板。

G(x,y)的表達(dá)式為：

(4)

SSR無(wú)法同時(shí)在動(dòng)態(tài)范圍壓縮、細(xì)節(jié)保持、色彩保真度、邊緣銳化上取得最佳效果，進(jìn)而有人提出了多尺度Retinex。

1.3 多尺度Retinex(MSR)

MSR與SSR的本質(zhì)相同：都是用高斯環(huán)繞函數(shù)與原圖像做卷積運(yùn)算，得到對(duì)光照?qǐng)D像的估計(jì)，進(jìn)而得到反應(yīng)物體本質(zhì)特征的反射圖像。MSR就是將多個(gè)SSR的結(jié)果進(jìn)行加權(quán)求和，其數(shù)學(xué)表達(dá)式為：

(5)

式中：Gk(x,y)為第k個(gè)高斯環(huán)繞函數(shù)；wk為第k個(gè)高斯的權(quán)重；K為尺度個(gè)數(shù)。

當(dāng)K=1時(shí)，MSR退化為SSR。一般情況下，取K=3，這樣3個(gè)尺度參數(shù)分別取小尺度、中尺度和大尺度，從而使得MSR兼?zhèn)淞薙SR低、中、高3個(gè)尺度的優(yōu)點(diǎn)，而權(quán)重一般取平均值，即有w1=w2=w3=1/3。至于尺度參數(shù)的具體取值，對(duì)于百萬(wàn)像素內(nèi)的圖像，文獻(xiàn)[11]建議3個(gè)尺度參數(shù)分別取15、80和250，故本文實(shí)驗(yàn)中也使用這3個(gè)尺度參數(shù)。

2 測(cè)試實(shí)驗(yàn)結(jié)果與分析

本文選用了5種常用的顏色空間來(lái)測(cè)試多尺度Retinex的性能，它們分別是RGB顏色空間[12]、HSV顏色空間[13]、YCbCr顏色空間[14]、YIQ顏色空間[15]、Lab顏色空間[16]。此外，選用了均方差(MSE)、峰值信噪比(PSNR)、信息熵(ER)、彩色增強(qiáng)因子(CEF)[15]、結(jié)構(gòu)相似度(SSIM)[17]和計(jì)算耗時(shí)共6個(gè)指標(biāo)來(lái)評(píng)價(jià)各顏色空間輸出圖像的增強(qiáng)效果。

圖2～圖5給出了5種顏色空間在不同圖像上的增強(qiáng)效果對(duì)比圖，即主觀評(píng)價(jià)；圖6～圖10為5種顏色空間對(duì)應(yīng)MSE、PSNR、ER、CEF、SSIM這5個(gè)指標(biāo)在20幅圖像上的分布，表1給出5種顏色空間的6種指標(biāo)在20幅圖像上的平均值，即客觀評(píng)價(jià)。

從主觀評(píng)價(jià)來(lái)看，多尺度Retinex算法在5種顏色空間上處理后都得到了比原圖像好的結(jié)果，即都有所增強(qiáng)。基于RGB、YCbCr、Lab這3種顏色空間的增強(qiáng)圖像在對(duì)比度增強(qiáng)方面優(yōu)于其他2種顏色空間，但有時(shí)過(guò)度增強(qiáng)使得圖像有點(diǎn)泛白，導(dǎo)致色彩失真比較嚴(yán)重，還會(huì)引入噪聲。在HSV和YIQ顏色空間上得到的結(jié)果具有較高的色彩保真度。

對(duì)于客觀評(píng)價(jià)，MSE值越小，表示失真越小，5種顏色空間中對(duì)應(yīng)HSV的MSE值最小，其次是YIQ；PSNR值越大，表示信噪比提升得越多，5種顏色空間中對(duì)應(yīng)HSV的PSNR值最大，其次是YIQ；ER值越大，表示從圖形中可以獲得信息內(nèi)容越多，5種顏色空間中對(duì)應(yīng)RGB的ER值最大，其次是HSV；CEF值越大，表示色彩越豐富，5種顏色空間中對(duì)應(yīng)RGB的CEF值最大，其次是HSV；SSIM越大，表示與原圖像越相似，5種顏色空間中對(duì)應(yīng)HSV的SSIM值最大，其次是YIQ；對(duì)于計(jì)算耗時(shí)來(lái)說(shuō)，越小越好，5種顏色空間中對(duì)應(yīng)YCbCr的計(jì)算耗時(shí)最小，其次是YIQ，雖然RGB顏色空間不存在轉(zhuǎn)換問(wèn)題，但由于需要在每個(gè)顏色通道做一次多尺度Retinex算法，而其他顏色空間僅僅在亮度通道上實(shí)施多尺度Retinex算法，故基于RGB顏色空間的計(jì)算耗時(shí)反而大于其他顏色空間。故從客觀評(píng)價(jià)來(lái)看，總體上來(lái)說(shuō)，多尺度Retinex算法在HSV彩色空間中取得的圖像增強(qiáng)效果好于其他4種顏色空間。

表1 5種顏色空間的6種指標(biāo)的平均值

3 結(jié)束語(yǔ)

本文首先研究了Retinex理論、單尺度Retinex算法和多尺度Retinex算法，接著介紹了RGB顏色空間、HSV顏色空間、YIQ顏色空間、YCbCr顏色空間、Lab顏色空間等5種顏色空間，給出了MSE、PSNR、信息熵、CEF、SSIM這5種圖像質(zhì)量評(píng)價(jià)指標(biāo)的計(jì)算公式，并用這5個(gè)指標(biāo)和計(jì)算耗時(shí)共6個(gè)指標(biāo)對(duì)多尺度Retinex算法在5種顏色空間上的增強(qiáng)效果做評(píng)價(jià)。實(shí)驗(yàn)結(jié)果表明，總體來(lái)說(shuō)，多尺度Retinex算法在HSV顏色空間上取得的增強(qiáng)效果好于其他4種彩色圖像。

雖然多尺度Retinex算法對(duì)霧天圖像和夜色圖像都有一定的增強(qiáng)，但由于在實(shí)驗(yàn)中的參數(shù)是固定的，故不能取得很好的增強(qiáng)效果，下一步將研究如何利用梯度和紋理等信息來(lái)自動(dòng)設(shè)置算法中的某些參數(shù)，以獲取更好的增強(qiáng)效果。

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

[2] 趙忠明,朱重光.遙感圖像中薄云的去除方法[J].遙感學(xué)報(bào),1996,11(3):195-199．

[3] 王彥臣,李樹(shù)杰,黃康卿.基于多尺度Retinex的數(shù)字X光圖像增強(qiáng)方法研究[J].光學(xué)精密工程,2006,14(1):70-77.

[4] 劉瑞劍, 陳樹(shù)越, 張甲杰.多尺度Retinex算法在紅外圖像增強(qiáng)中的應(yīng)用[J].彈箭與制導(dǎo)學(xué)報(bào),2008,28(4):193-195.

[5] RAHMAN Z U,JOBSON D J,WOODELL G A.Retinex processing for automatic image enhancement [J].Journal of Electronic Imaging,2004,13(1):100-110．

[6] HORN B K P.Determining lightness from an image [J].Computer graphics and image processing,1974,3(4):277-299.

[7] FUNT B,CIUREA F,MECANN J.Retinex in matlab [J].Journal of Electronic Imaging,2004,13(1):48- 57.

[8] JOBSON D J,RAHMAN Z U,WOODELL G A.Properties and performance of a center/surround retinex [J].IEEE Transactions on Image Processing,1995(6):451-462.

[9] RAHMAN Z U,JOBSON D J,WOODELL G A.Multi-scale retinex for color image enhancement [C]//Proceedings of The International Conference on Image Processing,Lausanne:IEEE,1996:1003-1006.

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

[11]ELAD M,KILNMEL R,SHAKED D,et al.Reduced complexity Retinex algorithm via the Variational approach [J].Journal of Visual Communication and Image Representation,2003,14(l):369-388.

[12]MURAHIRA K.Hue-preserving color image enhancement in RGB color space with rich saturation [C]// International Symposium on Intelligent Signal Processing and Communications Systems (ISPACS).New Taipei:IEEE,2012:266-269.

[13]REN C C,YANG J G.A novel color microscope image enhancement method based on HSV color space and curvelet transform [J].International Journal of Computer Science Issues,2012,9(2):272-277.

[14]LEE S,KWAK Y,KIM Y J,et al.Contrast-preserved chroma enhancement technique using YCbCr color space[J].IEEE Transactions on Consumer Electronics,2012,58(2):641-645.

[15]LAL S,NARASIMHADHAN A V,KUMAR R.Automatic method for contrast enhancement of natural color images [J].Journal of Electrical Engineer & Technology,2015,10(3):1233-1243.

[16]CHO Y H,KIM Y T,LEE H K,et al.Color image enhancement technique using gamut mapping based on color space division [C]//Electronic Imaging 2003.International Society for Optics and Photonics,2003:81-91.

[17]WANG Z,BOVIK A C,SHEIKH H R,et al.Image quality assessment:from error visibility to structural similarity[J].IEEE Transactions on Image Processing,2004,13(4):600-612.

Color Space Selection for Multi-scale Retinex Image Enhancement Application

LIU Li-ming

(The 723 Institute of CSIC,Yangzhou 225001,China)

This paper studies the image enhancement algorithm based on Retinex theory,introduces five color spaces:RGB,HSV,YIQ,YCbCrandLab,which are often used in image and video processing domain,gives the calculation formulas of five image quality evaluation indexes:mean square error (MSE),peak signal to noise ratio (PSNR),entropy,color enhancement factor (CEF),structural similarity (SSIM),which are used to evaluate the enhancement effect of CLAHE algorithm on five color spaces.The experimental results show:comparied with other four color spaces,multi-scale Retinex algorithm inHSVcolor space has achieved the best enhancement effect.

image enhancement;Retinex theory;multi-scale;color space

2016-05-12

TN911.73

A

CN32-1413(2017)01-0068-06

10.16426/j.cnki.jcdzdk.2017.01.015