一種改進(jìn)的色彩保持低照度圖像增強(qiáng)方法

遆曉光，曲悠楊

(哈爾濱工業(yè)大學(xué)控制與仿真中心，150080 哈爾濱)

一種改進(jìn)的色彩保持低照度圖像增強(qiáng)方法

遆曉光，曲悠楊

(哈爾濱工業(yè)大學(xué)控制與仿真中心，150080 哈爾濱)

針對現(xiàn)有低照度圖像增強(qiáng)方法恢復(fù)結(jié)果存在的色彩失真問題，采用原用于圖像去霧領(lǐng)域的暗通道先驗(yàn)方法，提出了一種基于暗通道先驗(yàn)的低照度圖像增強(qiáng)方法.通過提取圖像的暗通道信息來獲得對應(yīng)場景的光照強(qiáng)度分布圖;針對暗通道先驗(yàn)原理導(dǎo)致的白色像素區(qū)域暗通道數(shù)值計(jì)算錯(cuò)誤的問題，提出了一種通過計(jì)算像素的RGB三通道數(shù)值差值來判定錯(cuò)誤區(qū)域方法，并給出了具體的錯(cuò)誤修正方法;為了消除暗通道先驗(yàn)產(chǎn)生的塊效應(yīng)，對修改后的光照強(qiáng)度分布圖進(jìn)行平滑;基于圖像每個(gè)像素點(diǎn)光照強(qiáng)度與標(biāo)準(zhǔn)光照強(qiáng)度的比例關(guān)系對原圖像的RGB三通道進(jìn)行放大，實(shí)現(xiàn)色彩保持的圖像增強(qiáng);與現(xiàn)有方法進(jìn)行實(shí)驗(yàn)對比分析，從恢復(fù)結(jié)果的色彩失真情況以及峰值信噪比等方面驗(yàn)證了所提方法的有效性.

低照度圖像增強(qiáng);暗通道先驗(yàn);色彩保持;圖像平滑

在夜間低照度條件下，采集圖像時(shí)周圍環(huán)境光照條件不佳以及采集圖像設(shè)備本身的問題都會對采集的圖像質(zhì)量產(chǎn)生影響，使得到的圖像整體灰度值和對比度低，感興趣區(qū)域難以識別.另外，圖像中的一部分由于獲得光照不足或位于陰影區(qū)域?qū)е禄叶戎递^低，與背景融合在一起，使提取出原始圖像中這部分的信息出現(xiàn)困難.因此，增強(qiáng)低照度下圖像的亮度和對比度在民用如智能交通系統(tǒng)、智能安防系統(tǒng)、智能視頻監(jiān)控等以及軍用中如成像偵察系統(tǒng)和成像制導(dǎo)中都具有重要的應(yīng)用價(jià)值.

目前，已有很多種方法可以用來解決低照度圖像增強(qiáng)問題.如:直方圖均衡法、Retinex方法、MSRCR方法、空域彩色圖像增強(qiáng)方法［1］、梯度域圖像增強(qiáng)方法［2］、仿生彩色圖像增強(qiáng)方法［3］、基于照明-反射模型的同態(tài)濾波方法［4］、HDR方法等等.

局部直方圖均衡方法是先將圖像的每個(gè)像素的RGB通道按一個(gè)固定比例放大，然后以一個(gè)窗遍歷整個(gè)圖像，并在每個(gè)窗內(nèi)進(jìn)行直方圖均衡［5］.這種方法不僅運(yùn)算量大，而且當(dāng)圖像中同時(shí)存在較暗和較亮的區(qū)域時(shí)圖像增強(qiáng)的效果很差.Retinex方法是將圖像變換到對數(shù)域［6］，在對數(shù)域?qū)D像的照射光分量和反射光分量分離開，并根據(jù)圖像的反射光分量恢復(fù)圖像［7］，這種方法的缺點(diǎn)是很容易造成圖像的色彩失真.MSRCR［8］方法是在Retinex方法的基礎(chǔ)上加入了顏色恢復(fù)的功能.但是在實(shí)際應(yīng)用中恢復(fù)的顏色參數(shù)需要手動調(diào)節(jié).HDR方法是對同一場景連續(xù)拍攝多幅不同曝光量的圖像，利用每個(gè)曝光時(shí)間相對應(yīng)最佳細(xì)節(jié)的圖像來合成最終的HDR圖像［9］.這種方法在拍攝時(shí)需要保持照相器材的穩(wěn)定，拍攝時(shí)間較長，無法應(yīng)用于實(shí)時(shí)圖像和視頻增強(qiáng)領(lǐng)域，而且對于全局亮度很低的圖像增強(qiáng)效果較差.

本文提出的改進(jìn)色彩保持的低照度圖像增強(qiáng)方法，是基于單幅圖像，引入原用于圖像去霧領(lǐng)域的暗通道先驗(yàn)［10］的思想，可獲得低照度圖像場景的光照強(qiáng)度分布，在增強(qiáng)圖像亮度和對比度的同時(shí)可很好地保持圖像色彩信息，具有處理速度快和復(fù)原圖像質(zhì)量高的特點(diǎn).同時(shí)對目前應(yīng)用于去霧領(lǐng)域的暗通道先驗(yàn)方法中存在的對于白色區(qū)域處理效果差的缺點(diǎn)進(jìn)行了改進(jìn).

1 算法概述

低照度圖像本質(zhì)上是由于場景光照強(qiáng)度較低產(chǎn)生的.而事實(shí)上在大多數(shù)低照度圖像中，光照分布不均將導(dǎo)致圖像中出現(xiàn)較亮的區(qū)域和較暗的區(qū)域，給圖像增強(qiáng)帶來了困難.為了恢復(fù)低照度圖像，我們希望放大較暗區(qū)域的像素值，保持較亮區(qū)域的像素值.這就要先獲得圖像的光照強(qiáng)度分布圖，根據(jù)每個(gè)像素點(diǎn)不同的光照強(qiáng)度進(jìn)行不同系數(shù)的放大.準(zhǔn)確的估計(jì)圖像的光照強(qiáng)度分布圖是方法的核心問題，是實(shí)現(xiàn)色彩保持的圖像增強(qiáng)的關(guān)鍵.在本文中引入原用于去霧領(lǐng)域的暗通道先驗(yàn)方法來估計(jì)場景的光照強(qiáng)度分布，并對暗通道先驗(yàn)原理中存在的問題進(jìn)行了改進(jìn).

本文提出的色彩保持的低照度圖像增強(qiáng)方法包括以下步驟:

1)輸入圖像，并將其轉(zhuǎn)換到RGB彩色空間，獲得圖像的RGB三通道圖像;

2)計(jì)算圖像的RGB三通道圖像各像素點(diǎn)在RGB通道中的最小值，得到輸入圖像的暗通道圖像;

3)為了解決暗通道先驗(yàn)原理對白色區(qū)域不適用的問題.對圖像中的白色區(qū)域進(jìn)行判定，并修改白色區(qū)域像素的暗通道數(shù)值，得到初始光照強(qiáng)度分布圖;

4)對修改后的暗通道圖像進(jìn)行線性平滑，以去除由于暗通道先驗(yàn)原理產(chǎn)生的分塊效應(yīng)，得到平滑后的光照強(qiáng)度分布圖;

5)根據(jù)像素點(diǎn)在平滑后光照強(qiáng)度圖像上的灰度與期望光照強(qiáng)度的比例關(guān)系，對圖像上不同的像素點(diǎn)的RGB三個(gè)通道進(jìn)行增強(qiáng)，完成圖像增強(qiáng).

2 暗通道圖像的獲得

獲得準(zhǔn)確的光照強(qiáng)度分布圖是恢復(fù)低照度圖像的關(guān)鍵，本文將通過暗通道先驗(yàn)的方法來獲得圖像光照強(qiáng)度分布圖.

在自然場景中，圖像可以被描述為［10］

其中:I是指觀測到的圖像;J是期望得到的場景的真實(shí)輻射值;A是光照強(qiáng)度分布圖;t用來描述光線通過媒介透射到照相機(jī)過程中沒有被散射的部分，低照度圖像增強(qiáng)的目標(biāo)就是從I中復(fù)原J、A、t.為了獲得圖像的光照強(qiáng)度圖像，先要提取圖像的暗通道圖像.暗通道先驗(yàn)是通過對正常光照下戶外圖像的觀察得出的，即在絕大多數(shù)非天空的局部區(qū)域里，某一些像素總會有至少一個(gè)顏色通道具有很低的值.對于一幅圖像I，定義它每個(gè)像素點(diǎn)的RGB三通道的最小值構(gòu)成的圖像為Imin:

對應(yīng)的暗通道圖像B為

其中:Q(i，j)為以像素I(i，j)為中心大小為N×N的圖像塊，(x，y)為Q(i，j)內(nèi)的任意像素點(diǎn).

經(jīng)觀察得出，除了天空方位，在一片區(qū)域中總存在有一個(gè)像素的某個(gè)通道的數(shù)值很低且趨近于0.如果I是自然光照下的戶外圖像，把B稱為圖像I的暗通道圖像，并且把以上觀察得出的經(jīng)驗(yàn)性規(guī)律稱為暗通道先驗(yàn).暗通道值的產(chǎn)生主要有3個(gè)因素:

1)汽車、建筑物和城市中玻璃窗戶的陰影，或者是樹葉、樹與巖石等自然景觀的投影;

2)色彩鮮艷的物體或表面，在RGB的3個(gè)通道中有些通道的值很低;

3)顏色較暗的物體或者表面，例如灰暗色的樹干和石頭.

總之，自然景物中到處都是陰影或者彩色，這些景物的圖像的暗通道值總是很低的.可以根據(jù)這個(gè)原理估計(jì)低照度條件下的圖像光照強(qiáng)度分布.

根據(jù)式(2)、(3)可以獲得低照度圖像I的暗通道圖像B.根據(jù)式(1)可得

因?yàn)镴dark(i，j)約等于0，故公式(4)等價(jià)于

本文不討論雨霧帶來的影響，對于室內(nèi)和街道等景深基本一致的場景可以認(rèn)為t是一個(gè)不變的常數(shù).故可以獲得低照度圖像的光照強(qiáng)度部分圖 A，即

可知對于圖像I，光照強(qiáng)度分布圖A與暗通道圖像B存在大小為固定系數(shù)1-t的比例關(guān)系.如何獲得準(zhǔn)確的暗通道圖像B將是準(zhǔn)確估計(jì)圖像光照強(qiáng)度分布的關(guān)鍵.由于t會在后續(xù)計(jì)算中被約去，故不需要進(jìn)行估計(jì).

3 白色圖像區(qū)域暗通道數(shù)值的修正

上一節(jié)已經(jīng)獲得了圖像的暗通道圖像，但是在獲得的暗通道圖像中，明顯有一些區(qū)域灰度值較高，不符合暗通道圖像與光照強(qiáng)度分布圖成固定比例的情況，這將造成圖像光照分布圖的估計(jì)錯(cuò)誤.為了獲得準(zhǔn)確的暗通道圖像，在這一節(jié)將對這些區(qū)域進(jìn)行研究，并提出判定和修改暗通道數(shù)值的方法.

3.1 白色區(qū)域的定義

實(shí)驗(yàn)中發(fā)現(xiàn)暗通道先驗(yàn)原理對于絕大部分的顏色區(qū)域都是適用的，但是由于白色像素的RGB三通道數(shù)值比例大致為1∶1∶1，白色區(qū)域的像素的RGB三通道強(qiáng)度值非常接近，并不存在1個(gè)明顯的數(shù)值接近0的通道.所以暗通道先驗(yàn)原理對于圖像中的一些白色區(qū)域，如白色的墻壁、白色的物體等并不適用，這就會導(dǎo)致獲得的白色區(qū)域的暗通道圖像出現(xiàn)錯(cuò)誤.而暗通道圖像的數(shù)值的不準(zhǔn)確就會導(dǎo)致光照強(qiáng)度分布圖的估計(jì)錯(cuò)誤，從而影響低照度圖像增強(qiáng)質(zhì)量，這就需要對圖像中的白色區(qū)域進(jìn)行判定，并對白色區(qū)域?qū)?yīng)的暗通道圖像的數(shù)值進(jìn)行修正.

在圖1為原圖像與其所對應(yīng)的初始暗通道圖像，圖1(b)所示的白圈內(nèi)的區(qū)域即為暗通道圖像中暗通道數(shù)值被誤判的區(qū)域.在原圖中，顯然白圈圈出的區(qū)域和其附近區(qū)域的光照強(qiáng)度相同，而在暗通道先驗(yàn)得到的結(jié)果中則與附近區(qū)域有明顯區(qū)別.這將導(dǎo)致這部分圖像在增強(qiáng)后出現(xiàn)色彩失真的問題.通過對圖像中白色區(qū)域像素的RGB三通道數(shù)值進(jìn)行分析，可以發(fā)現(xiàn)白色區(qū)域的像素I(k，l)具有以下兩個(gè)特點(diǎn):

1)IR(k，l)、IG(k，l)、IB(k，l)3 個(gè)數(shù)之間數(shù)值較為接近;

2)B(k，l)數(shù)值較大.

圖1 原圖像對應(yīng)的暗通道圖像

3.2 白色區(qū)域的判定及暗通道數(shù)值修正

為了修正誤判的白色區(qū)域暗通道圖像的數(shù)值，必須先要對白色區(qū)域進(jìn)行判定，可采用如下方法.

對于已經(jīng)獲得的暗通道圖像B，計(jì)算B的最大灰度值Bmaxgray以及B的平均灰度值Bmeangray:

其中w、h分別為圖像I的寬和高.提取圖像的RGB三通道圖像IR、IG、IB，并對它們兩兩做差，取絕對值的最大值記為S.對于圖像的所有像素點(diǎn)，如果 S(i，j)≥10，B(i，j)≥ Bmeangray且 B(i，j)≤Bmaxgray/2，則像素點(diǎn)I(i，j)被判定為白色區(qū)域.對于所有判定為白色區(qū)域的像素，將其在暗通道圖像B上對應(yīng)的值修改為Bmeangray，修改后的暗通道圖像記為D.由于圖像I中存在的噪聲，在判定白色區(qū)域過程中出現(xiàn)的一些誤判的點(diǎn)，這些點(diǎn)會在之后的圖像平滑過程中得到修復(fù)，不會影響恢復(fù)后圖像的質(zhì)量.

綜上所述，白色區(qū)域判定及其暗通道數(shù)值修正的步驟為:

1)獲得圖像的RGB三通道圖像IR、IG、IB.

2)計(jì)算圖像中每個(gè)像素I(i，j)三通道間的最大差值 S(i，j).

3)對于圖像中每個(gè)像素I(i，j)判定是否滿足白色區(qū)域的條件.

4)修改判定為白色區(qū)域的像素所對應(yīng)的暗通道數(shù)值.

在修正了暗通道圖像中白色區(qū)域的數(shù)值以后，獲得了圖像的初始光照強(qiáng)度分布圖D.

4 初始光照強(qiáng)度分布圖的平滑

在上一節(jié)中獲得了低照度圖像初始光照強(qiáng)度分布圖，但是由于暗通道圖像的求取是對于每個(gè)像素點(diǎn)在一個(gè)以其為中心的窗內(nèi)取最小值的計(jì)算，當(dāng)存在一像素點(diǎn)其某一通道的數(shù)值極小時(shí)，包含這一像素點(diǎn)的所有窗對應(yīng)的像素點(diǎn)的暗通道數(shù)值都將為這一極小值，于是就產(chǎn)生了分塊特性.而實(shí)際上光照強(qiáng)度在某一局部區(qū)域應(yīng)該是平滑過渡的，為了使恢復(fù)的光照強(qiáng)度圖像均勻過渡還需要對初始光照強(qiáng)度分布圖D進(jìn)行線性平滑，得到平滑的光照強(qiáng)度圖像E.在文獻(xiàn)［10］中，采用了一種軟摳圖算法(soft matting)來完善透射率分布函數(shù)，這個(gè)方法的問題在于需要求解稀疏線性系統(tǒng)，導(dǎo)致計(jì)算量非常大，嚴(yán)重的影響了算法的運(yùn)行速度，本文采用線性近似的方法來代替軟摳圖算法.

4.1 圖像平滑的要求

初始光照分布圖的平滑需要滿足兩個(gè)要求:

1)對于暗通道圖像中灰度數(shù)值變化不明顯的低頻區(qū)域要起到平滑的作用;

2)對于暗通道圖像中灰度數(shù)值變化明顯的高頻區(qū)域要起到保持的作用，保證不會模糊邊緣的信息.

4.2 線性近似平滑

為滿足以上要求，對于圖像I上的每一點(diǎn)，對以其為中心，大小為 M ×M的窗 Wk(k=1，2，…，w×h)內(nèi)的所有像素以如下公式進(jìn)行線性近似［11］:

其中M的大小根據(jù)輸入圖像的大小w×h確定，對于同一個(gè)圖像，M應(yīng)略小于求取暗通道圖像時(shí)使用的窗的大小N的值，可以使用公式M=N-2來計(jì)算M的數(shù)值.在M的數(shù)值大小確定后，線性近似的結(jié)果將由近似參數(shù)ak，bk決定.當(dāng)ak，bk的數(shù)值發(fā)生變化的時(shí)候，線性近似的結(jié)果E也會發(fā)生變化.E和D之間的差由下面公式求得:

其中ε是一個(gè)非常小的常數(shù)用來防止ak過大.通過求取R(ak，bk)的最小值的方法來確定線性近似的參數(shù) ak，bk.

4.3 線性平滑參數(shù)的求取

使用線性回歸的方法解出R(ak，bk)的最小值為

當(dāng)R(ak，bk)為最小值時(shí)，對于低頻區(qū)域，D與uk近似相等.線性近似參數(shù)ak，bk的求解結(jié)果將趨近于 ak=0，bk=uk，E=uk，即對窗內(nèi)所有像素灰度值進(jìn)行平均;對于高頻區(qū)域，線性近似參數(shù)ak，bk的求解結(jié)果將趨近于ak=1，bk=0，E=D.即保持窗內(nèi)所有像素的原有灰度值，滿足了暗通道圖像線性近似的要求.

在求得所有的窗Wk的線性近似解ak、bk以后，由于一個(gè)像素可能會同時(shí)存在于很多個(gè)窗內(nèi)，在獲得最后的平滑結(jié)果時(shí)，應(yīng)該取所有包含該點(diǎn)的窗Wk的ak、bk的平均值用于計(jì)算，即

平滑后的光照分布圖如圖2所示.由圖2可以看出在保持細(xì)節(jié)的前提下，不但有效地去除了分塊效應(yīng)，而且對暗通道數(shù)值出現(xiàn)錯(cuò)誤的區(qū)域也進(jìn)行了修正，獲得了可以用于進(jìn)行低照度圖像增強(qiáng)的光照強(qiáng)度分布圖.

圖2 平滑后的光照強(qiáng)度分布

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

在獲得了平滑的光照強(qiáng)度分布圖以后，即可以對低照度圖像進(jìn)行增強(qiáng).與普通的圖像恢復(fù)不同，對于低照度圖像更關(guān)注如何使恢復(fù)的圖像具有良好的視覺效果而不僅僅是恢復(fù)所獲的圖像場景［12］.

5.1 標(biāo)準(zhǔn)光照強(qiáng)度的獲得

本文采用的方法是對圖像上每個(gè)像素點(diǎn)光照強(qiáng)度進(jìn)行調(diào)整.希望把低照度圖像恢復(fù)到原場景在均勻的強(qiáng)度較高的光照下的圖像，來減小場景光照不均勻以及光照不足的影響.恢復(fù)后的均勻照射下的場景光照強(qiáng)度稱為期望光照強(qiáng)度.

為了獲得圖像所對應(yīng)的期望光照強(qiáng)度，對大量低照度圖像的光照強(qiáng)度分布圖的直方圖進(jìn)行分析，發(fā)現(xiàn)他們具有比較類似的分布，即圖像像素主要分布于低灰度值區(qū)域.并且隨著灰度值的增大，像素分布逐漸減少.由于強(qiáng)光照區(qū)域的存在，在某些較大的灰度區(qū)域存在較多的像素分布.在大量實(shí)驗(yàn)中發(fā)現(xiàn)，圖像中細(xì)節(jié)清晰且比較容易辨識的區(qū)域所對應(yīng)的光照強(qiáng)度值有如下特點(diǎn):1)光照強(qiáng)度值大致位于光照強(qiáng)度分布圖上接近于低灰度值區(qū)域;2)低照度圖像中這些低灰度值對應(yīng)的像素?cái)?shù)量約為總像素?cái)?shù)量的1%.針對這一特點(diǎn)，使用下面的方法獲得期望光照強(qiáng)度Id，即在光照強(qiáng)度分布圖的直方圖上，所有像素個(gè)數(shù)小于w×h/100的灰度值中最小的一個(gè)灰度值被選擇作為期望光照強(qiáng)度(適用條件為w×h＞10 000).

5.2 圖像的增強(qiáng)

在獲得了期望光照強(qiáng)度以后，根據(jù)像素點(diǎn)的光照強(qiáng)度與期望光照強(qiáng)度的比例關(guān)系，對低照度圖像I中每個(gè)像素點(diǎn)的RGB 3個(gè)通道的灰度值進(jìn)行不同大小的放大，即可得到增強(qiáng)后的圖像.具體低照度圖像增強(qiáng)公式為

其中:J為增強(qiáng)后的低照度圖像，P是圖像增強(qiáng)的比例系數(shù).P的取值與原場景的光照強(qiáng)度為反比例關(guān)系.當(dāng)原場景中平均光照強(qiáng)度較低時(shí)，需要對圖像的各通道進(jìn)行較大系數(shù)的放大，這時(shí)P應(yīng)為較大的值.當(dāng)原場景中平均光照強(qiáng)度較高時(shí)，需要對圖像的各通道進(jìn)行較小系數(shù)的放大，這時(shí)P應(yīng)為較小的值，即P×Bmeangray為一常值.為了確定這一常值的大小，對8:00至16:00期間獲取的300幅室內(nèi)外自然光照圖像的P×Bmeangray值進(jìn)行了計(jì)算，發(fā)現(xiàn)實(shí)驗(yàn)中P×Bmeangray的平均值大致為9.即P的取值方法為

對于一般的夜間場景(即Bmeangray=1.1～1.2時(shí))P的大小為7～8.

6 實(shí)驗(yàn)對比

6.1 實(shí)驗(yàn)說明

為了驗(yàn)證本文方法的可行性和有效性，本節(jié)將從圖像增強(qiáng)視覺效果、峰值信噪比及運(yùn)算時(shí)間等方面與局部直方圖方法、Retinex方法、MSRCR方法和HDR方法進(jìn)行比較.實(shí)驗(yàn)用圖像大小為600×800.實(shí)驗(yàn)平臺采用處理器為Inter Core 2.20 GHz、內(nèi)存為3 GB，運(yùn)行Windows XP的PC機(jī)，編程語言為Matlab2011b.

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

為了獲得有效的實(shí)驗(yàn)結(jié)果，對于室外場景在18:00到21:00的時(shí)間范圍內(nèi)每3 min進(jìn)行一次場景圖像采集，并對這60幅圖片進(jìn)行對比實(shí)驗(yàn).對于室內(nèi)場景也采集了超過100幅圖片進(jìn)行實(shí)驗(yàn).在下面給出了幾組具有代表性的例子.圖3給出了對于低照度室外場景圖像分別采用局部直方圖均衡方法、Retinex方法、MSRCR方法以及本文提出的方法進(jìn)行增強(qiáng)后的對比結(jié)果.在圖4中給出了針對低照度室內(nèi)場景本文方法與上述三種方法進(jìn)行圖像增強(qiáng)后的實(shí)驗(yàn)對比結(jié)果.表1為針對低照度室外場景，本文方法與其他三種方法進(jìn)行圖像增強(qiáng)后的平均峰值信噪比(PSNR)對比結(jié)果.表2為針對低照度室內(nèi)場景，本文方法與其他三種方法圖像增強(qiáng)后平均峰值信噪比對比結(jié)果.表3為分別針對低照度室內(nèi)和室外場景，本文方法與其他三種方法進(jìn)行圖像增強(qiáng)平均運(yùn)行時(shí)間的對比結(jié)果.

圖3 室外場景低照度圖像增強(qiáng)的效果

表1 室外場景平均峰值信噪比對比 dB

圖4 室內(nèi)場景低照度圖像增強(qiáng)的效果

表2 室內(nèi)場景平均峰值信噪比對比 dB

表3 室內(nèi)外場景平均運(yùn)行時(shí)間對比 s

從圖3、4可以看出，局部直方圖均衡方法很大程度合并了圖像的灰度級，一定程度上造成了圖像的模糊，而且可以看出色彩失真現(xiàn)象非常明顯.采用Retinex方法進(jìn)行低照度圖像增強(qiáng)，增強(qiáng)圖像存在著明顯的色彩失真和噪聲大的問題.MSRCR方法用于低照度圖像增強(qiáng)，雖然可以突出圖像場景邊緣信息，但是也有很明顯的色彩失真情況.增強(qiáng)圖像中的大部分景物失去了原有的色彩，而且MSRCR方法中用于顏色補(bǔ)償?shù)膮?shù)需要手動設(shè)置，設(shè)置不當(dāng)也會影響圖像增強(qiáng)的效果.而采用本文的方法增強(qiáng)的圖像在色彩保持方面明顯優(yōu)于其他三種方法，丟失的色彩信息較少.

從表1和表2中可以看出本文方法的低照度圖像增強(qiáng)結(jié)果的峰值信噪比要遠(yuǎn)高于其他三種方法的峰值信噪比，說明本文方法失真程度遠(yuǎn)低于其他兩種方法.從表3可以看出，本文方法進(jìn)行低照度圖像增強(qiáng)相對其他三種方法在運(yùn)行時(shí)間方面也有了較明顯的提高.

圖5給出了原始低照度圖像分別采用HDR方法及本文提出的方法進(jìn)行增強(qiáng)后的對比結(jié)果，表4為兩種方法的峰值信噪比對比結(jié)果.

圖5 本文方法與HDR方法對比

表4 圖5中兩種方法的峰值信噪比對比 dB

從圖5以及表4中可以看出與HDR方法相比，本文方法對圖像增強(qiáng)的結(jié)果具有更高的亮度和對比度，而且圖像整體顏色更接近于真實(shí)圖像.在色彩保持的前提下，實(shí)現(xiàn)了圖像亮度和對比度的提升.同時(shí)在進(jìn)行低照度圖像增強(qiáng)時(shí)，HDR方法需要連續(xù)拍攝3到5幅圖片才能完成圖像增強(qiáng)的工作.在此期間如果攝像器材發(fā)生輕微移動都會導(dǎo)致圖像增強(qiáng)的結(jié)果出現(xiàn)模糊，從而對于攝像器材以及拍攝者的要求都比較高，而本文方法則不存在這一問題.

圖6給出了加入暗通道修正前后的對比實(shí)驗(yàn)處理結(jié)果.從圖6中可以看出，在未經(jīng)白色區(qū)域暗通道數(shù)值修正之前，白圈內(nèi)的區(qū)域在增強(qiáng)后顏色出現(xiàn)了嚴(yán)重的失真.而經(jīng)過暗通道數(shù)值修正以后，增強(qiáng)后的圖像色彩與原圖像保持一致，解決了色彩丟失的問題.

圖7給出了拍攝圖像期間攝像器材發(fā)生移動，HDR和本文方法進(jìn)行低照度圖像增強(qiáng)的對比結(jié)果.

由圖7可以看出，當(dāng)攝像器材發(fā)生移動時(shí)，使用HDR方法獲得的增強(qiáng)圖像出現(xiàn)了極大的圖像模糊和對比度低的問題，而本文方法則能在保證圖像亮度和對比度的前提下較好地實(shí)現(xiàn)圖像增強(qiáng).

圖6 加入暗通道修正前后結(jié)果對比

圖7 攝像器材移動情況下本文方法與HDR方法的對比

7 結(jié)語

針對夜間及低光照場景中存在的圖像亮度、對比度低的問題，基于暗通道先驗(yàn)原理，本文提出了一種具有色彩保持的可單幅圖像實(shí)現(xiàn)的低照度圖像增強(qiáng)的方法.實(shí)驗(yàn)結(jié)果表明，該方法不但可高質(zhì)量的增強(qiáng)低照度圖像，而且具有運(yùn)算速度快和易于實(shí)現(xiàn)的特點(diǎn)，恢復(fù)后的圖像更方便人眼以及機(jī)器進(jìn)行識別.本文方法可以用于智能視頻監(jiān)控系統(tǒng)以及可見光成像導(dǎo)彈夜間目標(biāo)捕捉，對于提高夜間環(huán)境下的目標(biāo)檢測識別準(zhǔn)確率具有重要作用.

［1］WAKIMOTO K，KANAZAWA Y，OHTA N.Color image enhancement for Dichromate by additive image noise［J］.IPSJ Transactions on Computer Vision and Applications，2013，5(1):45－49.

［2］許欣，陳強(qiáng)，孫懷江，等.基于梯度域融合的圖像視覺效果改善［J］.中國圖象圖形學(xué)報(bào)，2011，16(2):278－286.

［3］王守覺，丁興號，廖英豪，等.一種新的仿生彩色圖像增強(qiáng)方法［J］.電子學(xué)報(bào)，2008，36(10):1970－1973.

［4］張曉宇，胡士強(qiáng).基于同態(tài)濾波抑制光照變化的視頻分割算法［J］.模式識別與人工智能，2013，26(1):99－105.

［5］劉松濤，王維，殷福亮.基于動態(tài)廣義直方圖均衡的紅外圖像增強(qiáng)方法［J］.系統(tǒng)工程與電子技術(shù)，2010，32(7):1411－1414.

［6］汪榮貴，朱靜，楊萬挺，等.基于照度分割的局部多尺度 Retinex算 法 ［J］. 電 子 學(xué) 報(bào)，2010，38(5):1181－1186.

［7］JANG C Y，LIM J H，KIM Y H.A fast multi-scale Retinexalgorithm using dominantSSR in weights selection［C］//The 2012 InternationalSoC Design Conference (ISOCC). Piscataway: IEEE，2012:37－40.

［8］WANG G Z，HE G J.A modified multi scale retinex with color restoration algorithm for automatic enhancement of landsat-5 remote sensing image［J］.Advanced Materials Research，2012，341:893－897.

［9］BYOUNG-JY YuN， HEE-DONGH. A contrast enhancement method for HDR image using a modified Image formation model［J］.IEICE Transactions on Information and Systems，2012，95(4):1112－1119.

［10］HE K，SUN J，TANG X.Single image haze removal using dark channel prior［J］.IEEE Transactions on Pattern Analysis and Machine Intelligence，2011，33(12):2341－2353.

［11］HE K M，SUN J，TANG X O.Guided image filtering［J］.IEEE Transactionson Pattern Analysisand Machine Intelligence，2013，35(6):1397－1409.

An improved low illumination image enhancement algorithm with color preserving

DI Xiaoguang，QU Youyang
(Control and Simulation Center，Harbin Institute of Technology，150080 Harbin，China)

To solve the problem of color distortion existing in low illumination image enhancement method，a novel low illumination image enhancement method based on dark-channel prior with the capability of preserving scene color is presented.First，the intensity of every image pixel with dark-channel prior is estimated directly，and for dealing with wrong dark-channel values in white region，the regions with wrong grey values by calculating the minimum difference value of the RGB channels and correct the dark-channel values of the regions are founded.Then，a linear smoothing method is used to refine the block effect caused by darkchannel prior，and the enhanced image is achieved by scaling up the pixel intensity according to its value in dark-channel image.Finally，by comparing with the existing low illumination image enhancement method，the effectiveness of our method with better color preserving capability and higher PSNR is proved.

low illumination image enhancement;dark-channel prior;color preserving;image smoothing

TP394.1

A

0367－6234(2014)03－0001－07

2013－09－27.

國家自然科學(xué)基金資助項(xiàng)目(61174202);中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(HIT.NSRIF 201163);黑龍江省博士后科研啟動金資助項(xiàng)目(LBHQ10078).

遆曉光(1972—)，男，博士，副教授.

遆曉光，dixiaoguang@hit.edu.cn.

(編輯 張 宏)