CN105654437B - 一种对低照度图像的增强方法 - Google Patents

一种对低照度图像的增强方法 Download PDF

Info

Publication number
CN105654437B
CN105654437B CN201510996197.3A CN201510996197A CN105654437B CN 105654437 B CN105654437 B CN 105654437B CN 201510996197 A CN201510996197 A CN 201510996197A CN 105654437 B CN105654437 B CN 105654437B
Authority
CN
China
Prior art keywords
component
indicate
image
indicates
pixel
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201510996197.3A
Other languages
English (en)
Other versions
CN105654437A (zh
Inventor
李明
刘峰
干宗良
吕月圆
陈色桃
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GUANGDONG XUNTONG TECHNOLOGY Co Ltd
Original Assignee
GUANGDONG XUNTONG TECHNOLOGY Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by GUANGDONG XUNTONG TECHNOLOGY Co Ltd filed Critical GUANGDONG XUNTONG TECHNOLOGY Co Ltd
Priority to CN201510996197.3A priority Critical patent/CN105654437B/zh
Publication of CN105654437A publication Critical patent/CN105654437A/zh
Application granted granted Critical
Publication of CN105654437B publication Critical patent/CN105654437B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10016Video; Image sequence
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10024Color image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20172Image enhancement details

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种对低照度图像的增强方法,包括:获取待处理的低照度图像后,将其从RGB颜色空间转换到HSV颜色空间,获得色度分量、饱和度分量和亮度分量;基于Retinex算法,采用交替极小化方法将亮度分量分解为反射分量和照射分量;分别对照射分量和反射分量进行增强处理后,合成增强后的亮度分量;对饱和度分量做自适应调整处理后获得增强后的饱和度分量;将色度分量、增强后的亮度分量和增强后的饱和度分量合成新的HSV图像;将获得的新的HSV图像转化为RGB图像,并进行白平衡处理后,获得增强后的图像。本发明可以极大地提升低照度图像的清晰度,使细节可以再现,适用性强、鲁棒性高,可广泛应用于图像处理领域中。

Description

一种对低照度图像的增强方法
技术领域
本发明涉及图像处理领域,特别是涉及一种对低照度图像的增强方法。
背景技术
名词解释:
Retinex算法:Retinex是“Retina”(视网膜)和“Cortex”(大脑皮层)的缩写。Retine算法是一种建立在科学实验和科学分析基础上的基于人类视觉系统的图像增强算法,可以在动态范围压缩、边缘增强和颜色恒常三方面达到平衡,对各种不同类型的图像进行自适应性的增强;
RGB:包括红色分量R、绿色分量G、蓝色分量B三个分量。
HSV:包括色度分量H、饱和度分量S、亮度分量V三个分量。
计算机视觉监控设备在室内、夜间等低照度条件下,由于非自然光源的照度不充分,所以目标表面的反射光较弱,造成进入成像传感器的光线不足,导致在夜间采集的图像质量退化严重、图像可辨识性很低,且含有大量噪声,以至于难以分辨图像中的细节,大大降低了图像的应用价值,这类图像统称为低照度图像。
对低照度图像的质量提升,是目前图像质量提升领域的研究热点之一,尤其是在城市交通、监控视频等计算机视觉领域,低照度图像的质量提升具有重大意义。目前对低照度图像的增强方法主要是在RGB颜色空间下直接对图像进行增强处理例如MSRCR算法,这种处理方式容易产生颜色丢失,不符合人眼的视觉特征。
发明内容
为了解决上述的技术问题,本发明的目的是提供一种对低照度图像的增强方法。
本发明解决其技术问题所采用的技术方案是:
一种对低照度图像的增强方法,包括:
S1、获取待处理的低照度图像后,将其从RGB颜色空间转换到HSV颜色空间,获得色度分量、饱和度分量和亮度分量;
S2、基于Retinex算法,采用交替极小化方法将亮度分量分解为反射分量和照射分量;
S3、分别对照射分量和反射分量进行增强处理后,合成增强后的亮度分量;
S4、对饱和度分量做自适应调整处理后获得增强后的饱和度分量;
S5、将色度分量、增强后的亮度分量和增强后的饱和度分量合成新的HSV图像;
S6、将获得的新的HSV图像转化为RGB图像,并进行白平衡处理后,获得增强后的图像。
进一步,所述步骤S2,其具体为:
构建如下目标函数,并根据该目标函数获得最优的照射分量和反射分量:
上式中,V表示亮度分量,L表示照射分量,R表示反射分量,表示经过最大值滤波后的低照度图像,IC(y)表示低照度图像,Ω表示滤波器窗口,α、β、λ均为预设的权重因子。
进一步,所述步骤S2,包括:
S21、构建如下目标函数:
上式中,V表示亮度分量,L表示照射分量,R表示反射分量,表示经过最大值滤波后的低照度图像,IC(y)表示低照度图像,Ω表示滤波器窗口,α、β、λ均为预设的权重因子;
S22、构建的等价表达式:其中p表示 满足条件的像素点的个数;
S23、根据下式获得的辅助变量:
其中,hp表示的辅助变量,vp表示的辅助变量,(hp,vp)满足γi表示迭代变量,γ1=2β=0.05且γi+1= 2*γi
S24、将亮度分量通过高斯低通滤波器滤波后的值作为照射分量的初始值;
S25、根据目标函数获得最优的反射分量如下:
其中,L0表示照射分量的初始值,f表示FFT变换,f-1表示FFT逆变换,f()*表示FFT变换后复共轭,h表示的辅助变量;
S26、根据下式对反射变量进行矫正:R=min(max(R,0),1);
S27、保持矫正后的反射分量的值不变,根据目标函数获得最优的照射分量如下:
S28、根据下式对照射变量进行矫正:L=max(L,V);
S29、判断迭代变量γi是否大于预设阈值,若是,则结束,反之将迭代次数加1后,将矫正后的照射分量作为照射分量的初始值,返回执行步骤S25。
进一步,所述步骤S3,包括:
S31、对照射分量进行伽马校正,并进行自适应直方图均衡化处理后,获得增强后的照射分量;
S32、对反射分量进行局部滤波操作后获得增强后的反射分量;
S33、将增强后的照射分量和增强后的反射分量合成增强后的亮度分量。
进一步,所述步骤S31中所述对照射分量进行伽马校正的步骤,其具体为:
根据下式对照射分量进行伽马校正:
上式中,L表示照射分量,Lg表示伽马校正后的照射分量。
进一步,所述步骤S32,其具体为:
根据下式,对反射分量进行局部滤波操作后,获得增强后的反射分量:
上式中,RF(x,y)表示对像素点(x,y)的反射分量进行局部滤波操作后的反射分量值,Ω表示以像素点(x,y)为中心的n×n的邻域窗口,(i,j)表示该邻域窗口中的像素坐标,R(i,j)表示像素点(i,j)的反射分量值,K表示归一化常数,ws(i,j)表示像素点(i,j)的空域权重,wr(i,j)表示像素点(i,j)的值域权重,且K、ws(i,j)、wr(i,j)满足以下条件:
其中,δS表示空域标准差,δr表示值域标准差,Rm,y表示像素点(m,y)的强度值,Ri,n表示像素点(i,n)的强度值。
进一步,所述步骤S4,其具体为:
根据下式,对饱和度分量做自适应调整处理后获得增强后的饱和度分量:
S'=S+t×(V'-V)×μ
上式中,S'表示增强后的饱和度分量,S表示饱和度分量,t为常数,V表示亮度分量,V'表示增强后的亮度分量,μ表示调整系数,且对于每个像素点(x,y),调整系数为:
其中,(x,y)表示像素点位置,Ω表示以像素点(x,y)为中心的n×n的邻域窗口,(i,j)表示该邻域窗口Ω中的像素坐标,V(i,j)表示像素点(i,j)的亮度值,表示邻域窗口Ω内的亮度均值,S(i,j)表示像素点(i,j)的饱和度值,表示邻域窗口Ω内的饱和度均值,δV(x,y)表示像素点(x,y)的亮度方差,δS(x,y)表示像素点的(x,y)的饱和度方差。
进一步,所述步骤S6中所述进行白平衡处理的步骤,包括:
S61、计算获得转化后的RGB图像的R、G、B三通道的颜色均值,并根据下式计算RGB通道的总均值:
Kave=(Rave+Gave+Bave)/3
上式中,Kave表示RGB通道的总均值,Rave、GaveBave分别表示R、G、B三通道的颜色均值;
S62、判断下式是否成立,若是,则直接结束,反之,执行步骤S63:
S63、根据下式对R、G、B颜色分量进行白平衡处理:
上式中,Rmod、Gmod、Bmod分别表示白平衡处理后的R、G、B三通道的颜色值,Rchannel、Gchannel、Bchannel分别表示白平衡处理前R、G、B通道的灰度值。
本发明的有益效果是:本发明的一种对低照度图像的增强方法,包括:获取待处理的低照度图像后,将其从RGB颜色空间转换到HSV颜色空间,获得色度分量、饱和度分量和亮度分量;基于Retinex算法,采用交替极小化方法将亮度分量分解为反射分量和照射分量;分别对照射分量和反射分量进行增强处理后,合成增强后的亮度分量;对饱和度分量做自适应调整处理后获得增强后的饱和度分量;将色度分量、增强后的亮度分量和增强后的饱和度分量合成新的HSV图像;将获得的新的HSV图像转化为RGB图像,并进行白平衡处理后,获得增强后的图像。本方法可以极大地提升低照度图像的清晰度,使细节可以再现,而且基本不会产生颜色丢失,更为符合人眼的视觉特征,适用性强、鲁棒性高。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1是本发明的一种对低照度图像的增强方法的流程示意图;
图2是本发明的具体实施例中的待处理低照度图像;
图3是采用传统增强方法对图2的低照度图像进行增强处理后获得的结果;
图4是本发明的具体实施例中对图2的低照度图像进行增强处理后获得的结果。
具体实施方式
参照图1,本发明提供了一种对低照度图像的增强方法,包括:
S1、获取待处理的低照度图像后,将其从RGB颜色空间转换到HSV颜色空间,获得色度分量、饱和度分量和亮度分量;
S2、基于Retinex算法,采用交替极小化方法将亮度分量分解为反射分量和照射分量;
S3、分别对照射分量和反射分量进行增强处理后,合成增强后的亮度分量;
S4、对饱和度分量做自适应调整处理后获得增强后的饱和度分量;
S5、将色度分量、增强后的亮度分量和增强后的饱和度分量合成新的HSV图像;
S6、将获得的新的HSV图像转化为RGB图像,并进行白平衡处理后,获得增强后的图像。
进一步作为优选的实施方式,所述步骤S2,其具体为:
构建如下目标函数,并根据该目标函数获得最优的照射分量和反射分量:
上式中,V表示亮度分量,L表示照射分量,R表示反射分量,表示经过最大值滤波后的低照度图像,IC(y)表示低照度图像,Ω表示滤波器窗口,α、β、λ均为预设的权重因子。
进一步作为优选的实施方式,所述步骤S2,包括:
S21、构建如下目标函数:
上式中,V表示亮度分量,L表示照射分量,R表示反射分量,表示经过最大值滤波后的低照度图像,IC(y)表示低照度图像,Ω表示滤波器窗口,α、β、λ均为预设的权重因子;
S22、构建的等价表达式:其中p表示 满足条件的像素点的个数;
S23、根据下式获得的辅助变量:
其中,hp表示的辅助变量,vp表示的辅助变量,(hp,vp)满足γi表示迭代变量,γ1=2β=0.05且γi+1= 2*γi
S24、将亮度分量通过高斯低通滤波器滤波后的值作为照射分量的初始值;
S25、根据目标函数获得最优的反射分量如下:
其中,L0表示照射分量的初始值,f表示FFT变换,f-1表示FFT逆变换,f()*表示FFT变换后复共轭,h表示的辅助变量;
S26、根据下式对反射变量进行矫正:R=min(max(R,0),1);
S27、保持矫正后的反射分量的值不变,根据目标函数获得最优的照射分量如下:
S28、根据下式对照射变量进行矫正:L=max(L,V);
S29、判断迭代次数是是否大于预设阈值,若是,则结束,反之将迭代次数加1后,将矫正后的照射分量作为照射分量的初始值,返回执行步骤S25。
进一步作为优选的实施方式,所述步骤S3,包括:
S31、对照射分量进行伽马校正,并进行自适应直方图均衡化处理后,获得增强后的照射分量;
S32、对反射分量进行局部滤波操作后获得增强后的反射分量;
S33、将增强后的照射分量和增强后的反射分量合成增强后的亮度分量。
进一步作为优选的实施方式,所述步骤S31中所述对照射分量进行伽马校正的步骤,其具体为:
根据下式对照射分量进行伽马校正:
上式中,L表示照射分量,Lg表示伽马校正后的照射分量。
进一步作为优选的实施方式,所述步骤S32,其具体为:
根据下式,对反射分量进行局部滤波操作后,获得增强后的反射分量:
上式中,RF(x,y)表示对像素点(x,y)的反射分量进行局部滤波操作后的反射分量值,Ω表示以像素点(x,y)为中心的n×n的邻域窗口,(i,j)表示该邻域窗口中的像素坐标,R(i,j)表示像素点(i,j)的反射分量值,K表示归一化常数,ws(i,j)表示像素点(i,j)的空域权重,wr(i,j)表示像素点(i,j)的值域权重,且K、ws(i,j)、wr(i,j)满足以下条件:
其中,δS表示空域标准差,δr表示值域标准差,Rm,y表示像素点(m,y)的强度值,Ri,n表示像素点(i,n)的强度值。
进一步作为优选的实施方式,所述步骤S4,其具体为:
根据下式,对饱和度分量做自适应调整处理后获得增强后的饱和度分量:
S'=S+t×(V'-V)×μ
上式中,S'表示增强后的饱和度分量,S表示饱和度分量,t为常数,V表示亮度分量,V'表示增强后的亮度分量,μ表示调整系数,且对于每个像素点(x,y),调整系数为:
其中,(x,y)表示像素点位置,Ω表示以像素点(x,y)为中心的n×n的邻域窗口,(i,j)表示该邻域窗口Ω中的像素坐标,V(i,j)表示像素点(i,j)的亮度值,表示邻域窗口Ω内的亮度均值,S(i,j)表示像素点(i,j)的饱和度值,表示邻域窗口Ω内的饱和度均值,δV(x,y)表示像素点(x,y)的亮度方差,δS(x,y)表示像素点的(x,y)的饱和度方差。
进一步作为优选的实施方式,所述步骤S6中所述进行白平衡处理的步骤,包括:
S61、计算获得转化后的RGB图像的R、G、B三通道的颜色均值,并根据下式计算RGB通道的总均值:
Kave=(Rave+Gave+Bave)/3
上式中,Kave表示RGB通道的总均值,Rave、GaveBave分别表示R、G、B三通道的颜色均值;
S62、判断下式是否成立,若是,则直接结束,反之,执行步骤S63:
S63、根据下式对R、G、B颜色分量进行白平衡处理:
上式中,Rmod、Gmod、Bmod分别表示白平衡处理后的R、G、B三通道的颜色值,Rchannel、Gchannel、Bchannel分别表示白平衡处理前R、G、B通道的灰度值。
以下结合一具体实施例对本发明做详细说明。
参照图1,一种对低照度图像的增强方法,包括:
S1、获取待处理的低照度图像后,将其从RGB颜色空间转换到HSV颜色空间,获得色度分量、饱和度分量和亮度分量;
S2、基于Retinex算法,采用交替极小化方法将亮度分量分解为反射分量和照射分量,具体为:构建如下目标函数,并根据该目标函数获得最优的照射分量和反射分量:
上式中,V表示亮度分量,L表示照射分量,R表示反射分量,表示经过最大值滤波后的低照度图像,IC(y)表示低照度图像,Ω表示滤波器窗口,α、β、λ均为预设的权重因子。
更详细的,步骤S2包括S21~S27:
S21、构建如下目标函数:
上式中,V表示亮度分量,L表示照射分量,R表示反射分量,表示经过最大值滤波后的低照度图像,IC(y)表示低照度图像,Ω表示滤波器窗口,α、β、λ均为预设的权重因子;
S22、构建的等价表达式:其中p表示 满足条件的像素点的个数;C(R)表示求取梯度不为0的像素点P的操作,P点 的梯度为
S23、根据下式获得的辅助变量:
其中,该公式中if表示符合什么条件时,if others表示其它情况时;其中,hp表示 的辅助变量,vp表示的辅助变量,(hp,vp)满足 γi表示迭代变量,γ1=2β=0.05且每次迭代后迭代变量变化为γi+1=2*γi;步骤S29中, 预设阈值设为0.5即多次迭代直到γi>0.5;此处通过γi控制迭代次数约5次;
由于L0范数是统计向量中非0的元素个数,无法直接进行求导运算,因此将上式中(hp,vp)的表达式代入步骤S21中的目标函数,将原目标函数改为:
其中,H(|hp|+|vp|)是一个二元函数,当|hp|+|vp|≠0时,返回1;其他情况,返回0;
S24、将亮度分量通过高斯低通滤波器滤波后的值作为照射分量的初始值;
S25、根据目标函数获得最优的反射分量如下:
其中,L0表示照射分量的初始值,f表示FFT变换,f-1表示FFT逆变换,f()*表示FFT变换后复共轭,h表示的辅助变量;
S26、根据下式对反射变量进行矫正:R=min(max(R,0),1);
S27、保持矫正后的反射分量的值不变(这里相当于将矫正后的反射分量作为求取照射分量时的反射分量初始值),根据目标函数获得最优的照射分量如下:
S28、根据下式对照射变量进行矫正:L=max(L,V);
S29、判断迭代变量γi是否大于预设阈值,若是,则结束,反之将迭代次数加1后,将矫正后的照射分量作为照射分量的初始值,返回执行步骤S25。本实施例中,优选将预设阈值设为0.5,通过大约5次迭代后,可满足γi>0.5,可以获得趋于稳定的照射分量L和反射分量R。本实施例通过γi控制迭代次数,优选控制迭代次数约为5次。
S3、分别对照射分量和反射分量进行增强处理后,合成增强后的亮度分量;步骤S3包括S31~S33:
S31、根据下式对照射分量进行伽马校正,获得伽马校正后的照射分量Lg,并对其进行自适应直方图均衡化处理后,获得增强后的照射分量L':
上式中,L表示照射分量,Lg表示伽马校正后的照射分量;
由于光照不均匀因素,得到的照射分量L存在明显的亮度明暗区域,需要对其进行伽马校正。最后,采用现有技术的直方图均衡化方法进行直方图均衡化处理后获得增强后的照射分量L':L'=CLAHE(Lg),CLAHE表示自适应直方图均衡函数。
S32、对反射分量进行局部滤波操作后获得增强后的反射分量;图像的反射分量是图像中变化较快的高频信息,可以反映图像上物体的内在性质,在采用Retinex理论提取获得反射分量的过程中,由于噪声常常分布在高频,噪声被放大了,因此本步骤需要进行局部滤波实现降噪处理,具体如下:
根据下式,对反射分量进行局部滤波操作后,获得增强后的反射分量:
上式中,RF(x,y)表示对像素点(x,y)的反射分量进行局部滤波操作后的反射分量值,Ω表示以像素点(x,y)为中心的n×n的邻域窗口,(i,j)表示该邻域窗口中的像素坐标,R(i,j)表示像素点(i,j)的反射分量值,K表示归一化常数,ws(i,j)表示像素点(i,j)的空域权重,wr(i,j)表示像素点(i,j)的值域权重,且K、ws(i,j)、wr(i,j)满足以下条件:
其中,δS表示空域标准差,δr表示值域标准差,Rm,y表示像素点(m,y)的强度值,Ri,n表示像素点(i,n)的强度值。
S33、将增强后的照射分量和增强后的反射分量合成增强后的亮度分量。
S4、对饱和度分量做自适应调整处理后获得增强后的饱和度分量,具体如下:
根据下式,对饱和度分量做自适应调整处理后获得增强后的饱和度分量:
S'=S+t×(V'-V)×μ
上式中,S'表示增强后的饱和度分量,S表示饱和度分量,t为常数,V表示亮度分量,V'表示增强后的亮度分量,μ表示调整系数,且对于每个像素点(x,y),调整系数为:
其中,(x,y)表示像素点位置,Ω表示以像素点(x,y)为中心的n×n的邻域窗口,(i,j)表示该邻域窗口Ω中的像素坐标,V(i,j)表示像素点(i,j)的亮度值,表示邻域窗口Ω内的亮度均值,S(i,j)表示像素点(i,j)的饱和度值,表示邻域窗口Ω内的饱和度均值,δV(x,y)表示像素点(x,y)的亮度方差,δS(x,y)表示像素点的(x,y)的饱和度方差。
S5、处理过程中,色度分量保持不变,本步骤中,将色度分量、增强后的亮度分量和增强后的饱和度分量合成新的HSV图像;
S6、将获得的新的HSV图像转化为RGB图像,并进行白平衡处理后,获得增强后的图像。
步骤S6中所述进行白平衡处理的步骤,包括:
S61、计算获得转化后的RGB图像的R、G、B三通道的颜色均值,并根据下式计算RGB通道的总均值:
Kave=(Rave+Gave+Bave)/3
上式中,Kave表示RGB通道的总均值,Rave、GaveBave分别表示R、G、B三通道的颜色均值;
S62、判断下式是否成立,若是,则直接结束,不进行白平衡处理,反之下式不成立 表示执行步骤S63:
S63、根据下式对R、G、B颜色分量进行白平衡处理:
上式中,Rmod、Gmod、Bmod分别表示白平衡处理后的R、G、B三通道的颜色值,Rchannel、Gchannel、Bchannel分别表示白平衡处理前R、G、B通道的灰度值。
图2是采集获得的一低照度图像,采用本实施例对图2的低照度图像进行增强处理,所获得的结果如图4所示,增强后的图像细节清晰可见,图像颜色正常。而图3是采用传统的MSRCR算法对图2的图像进行增强处理后获得的结果,通过图3和图4的对比可知,本方法相对于传统的在RGB颜色空间下直接对图像进行增强处理的算法,增强效果好,而且基本不会产生颜色丢失,更为符合人眼的视觉特征。因此,本发明可以极大地提升低照度图像的清晰度,使细节可以再现,本方法适用性强、鲁棒性高。
以上是对本发明的较佳实施进行了具体说明,但本发明创造并不限于所述实施例,熟悉本领域的技术人员在不违背本发明精神的前提下还可做出种种的等同变形或替换,这些等同的变型或替换均包含在本申请权利要求所限定的范围内。

Claims (7)

1.一种对低照度图像的增强方法,其特征在于,包括:
S1、获取待处理的低照度图像后,将其从RGB颜色空间转换到HSV颜色空间,获得色度分量、饱和度分量和亮度分量;
S2、基于Retinex算法,采用交替极小化方法将亮度分量分解为反射分量和照射分量;
S3、分别对照射分量和反射分量进行增强处理后,合成增强后的亮度分量;
S4、对饱和度分量做自适应调整处理后获得增强后的饱和度分量;
S5、将色度分量、增强后的亮度分量和增强后的饱和度分量合成新的HSV图像;
S6、将获得的新的HSV图像转化为RGB图像,并进行白平衡处理后,获得增强后的图像;
所述步骤S2,其具体为:
构建如下目标函数,并根据该目标函数获得最优的照射分量和反射分量:
上式中,V表示亮度分量,L表示照射分量,R表示反射分量,表示经过最大值滤波后的低照度图像,IC(y)表示低照度图像,Ω表示滤波器窗口,α、β、λ均为预设的权重因子。
2.根据权利要求1所述的一种对低照度图像的增强方法,其特征在于,所述步骤S2,包括:
S21、构建如下目标函数:
上式中,V表示亮度分量,L表示照射分量,R表示反射分量,表示经过最大值滤波后的低照度图像,IC(y)表示低照度图像,Ω表示滤波器窗口,α、β、λ均为预设的权重因子;
S22、构建的等价表达式:其中p表示满足条件的像素点的个数,C(R)表示求取梯度不为0的像素点p的个数的操作;
S23、根据下式获得的辅助变量:
其中,hp表示的辅助变量,vp表示的辅助变量,(hp,vp)满足γi表示迭代变量,γ1=2β=0.05且γi+1=2*γiH(|hp|+|vp|)是一个二元函数,当|hp|+|vp|≠0时,返回1;其他情况,返回0;
S24、将亮度分量通过高斯低通滤波器滤波后的值作为照射分量的初始值;
S25、根据目标函数获得最优的反射分量如下:
其中,L0表示照射分量的初始值,f表示FFT变换,f-1表示FFT逆变换,f()*表示FFT变换后复共轭,h表示的辅助变量;
S26、根据下式对反射变量进行矫正:R=min(max(R,0),1);
S27、保持矫正后的反射分量的值不变,根据目标函数获得最优的照射分量如下:
S28、根据下式对照射变量进行矫正:L=max(L,V);
S29、判断迭代变量γi是否大于预设阈值,若是,则结束,反之将迭代次数加1后,将矫正后的照射分量作为照射分量的初始值,返回执行步骤S25。
3.根据权利要求1所述的一种对低照度图像的增强方法,其特征在于,所述步骤S3,包括:
S31、对照射分量进行伽马校正,并进行自适应直方图均衡化处理后,获得增强后的照射分量;
S32、对反射分量进行局部滤波操作后获得增强后的反射分量;
S33、将增强后的照射分量和增强后的反射分量合成增强后的亮度分量。
4.根据权利要求3所述的一种对低照度图像的增强方法,其特征在于,所述步骤S31中所述对照射分量进行伽马校正的步骤,其具体为:
根据下式对照射分量进行伽马校正:
γ=2.2
上式中,L表示照射分量,Lg表示伽马校正后的照射分量。
5.根据权利要求3所述的一种对低照度图像的增强方法,其特征在于,所述步骤S32,其具体为:
根据下式,对反射分量进行局部滤波操作后,获得增强后的反射分量:
上式中,RF(x,y)表示对像素点(x,y)的反射分量进行局部滤波操作后的反射分量值,Ω表示以像素点(x,y)为中心的n×n的邻域窗口,(i,j)表示该邻域窗口中的像素坐标,R(i,j)表示像素点(i,j)的反射分量值,K表示归一化常数,ws(i,j)表示像素点(i,j)的空域权重,wr(i,j)表示像素点(i,j)的值域权重,且K、ws(i,j)、wr(i,j)满足以下条件:
其中,σS表示空域标准差,σr表示值域标准差,Rm,y表示像素点(m,y)的强度值,Ri,n表示像素点(i,n)的强度值。
6.根据权利要求1所述的一种对低照度图像的增强方法,其特征在于,所述步骤S4,其具体为:
根据下式,对饱和度分量做自适应调整处理后获得增强后的饱和度分量:
S'=S+t×(V'-V)×μ
上式中,S'表示增强后的饱和度分量,S表示饱和度分量,t为常数,V表示亮度分量,V'表示增强后的亮度分量,μ表示调整系数,且对于每个像素点(x,y),调整系数为:
其中,(x,y)表示像素点位置,Ω表示以像素点(x,y)为中心的n×n的邻域窗口,(i,j)表示该邻域窗口Ω中的像素坐标,V(i,j)表示像素点(i,j)的亮度值,表示邻域窗口Ω内的亮度均值,S(i,j)表示像素点(i,j)的饱和度值,表示邻域窗口Ω内的饱和度均值,δV(x,y)表示像素点(x,y)的亮度方差,δS(x,y)表示像素点(x,y)的饱和度方差。
7.根据权利要求1所述的一种对低照度图像的增强方法,其特征在于,所述步骤S6中所述进行白平衡处理的步骤,包括:
S61、计算获得转化后的RGB图像的R、G、B三通道的颜色均值,并根据下式计算RGB通道的总均值:
Kave=(Rave+Gave+Bave)/3
上式中,Kave表示RGB通道的总均值,Rave、Gave、Bave分别表示R、G、B三通道的颜色均值;
S62、判断下式是否成立,若是,则直接结束,反之,执行步骤S63:
S63、根据下式对R、G、B颜色分量进行白平衡处理:
上式中,Rmod、Gmod、Bmod分别表示白平衡处理后的R、G、B三通道的颜色值,Rchannel、Gchannel、Bchannel分别表示白平衡处理前R、G、B通道的灰度值。
CN201510996197.3A 2015-12-24 2015-12-24 一种对低照度图像的增强方法 Expired - Fee Related CN105654437B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510996197.3A CN105654437B (zh) 2015-12-24 2015-12-24 一种对低照度图像的增强方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510996197.3A CN105654437B (zh) 2015-12-24 2015-12-24 一种对低照度图像的增强方法

Publications (2)

Publication Number Publication Date
CN105654437A CN105654437A (zh) 2016-06-08
CN105654437B true CN105654437B (zh) 2019-04-19

Family

ID=56476928

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510996197.3A Expired - Fee Related CN105654437B (zh) 2015-12-24 2015-12-24 一种对低照度图像的增强方法

Country Status (1)

Country Link
CN (1) CN105654437B (zh)

Families Citing this family (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106485674A (zh) * 2016-09-20 2017-03-08 天津大学 一种基于融合技术的弱光图像增强方法
CN106530250A (zh) * 2016-11-07 2017-03-22 湖南源信光电科技有限公司 一种基于改进的Retinex低照度彩色图像增强方法
CN106651818A (zh) * 2016-11-07 2017-05-10 湖南源信光电科技有限公司 基于改进的直方图均衡化低照度图像增强算法
CN106780379B (zh) * 2016-12-08 2019-08-23 哈尔滨工业大学 一种计量显微镜的彩色图像增强方法
CN106981052B (zh) * 2017-01-12 2020-03-31 西南电子技术研究所(中国电子科技集团公司第十研究所) 基于变分框架的自适应亮度不均变分校正方法
CN106846282B (zh) * 2017-03-28 2019-06-04 华侨大学 一种采用自适应校正的低照度图像增强方法
CN107194882A (zh) * 2017-03-29 2017-09-22 南京工程学院 一种钢丝绳芯输送带x光图像矫正与增强的方法
CN106897981A (zh) * 2017-04-12 2017-06-27 湖南源信光电科技股份有限公司 一种基于引导滤波的低照度图像增强方法
CN109255756B (zh) * 2017-07-14 2020-12-29 北京大学 低光照图像的增强方法及装置
CN107392869B (zh) * 2017-07-21 2020-12-01 长安大学 一种基于边缘保持滤波器的人脸图像滤波方法
CN107507215B (zh) * 2017-08-07 2020-07-14 广东电网有限责任公司珠海供电局 一种基于自适应量化增强的电力设备红外热图分割方法
CN107578383A (zh) * 2017-08-29 2018-01-12 北京华易明新科技有限公司 一种低照度图像增强处理方法
CN107443403B (zh) * 2017-09-21 2019-11-01 上海昕鼎网络科技有限公司 一种用于巡逻的机器人
CN107656454A (zh) * 2017-09-21 2018-02-02 深圳市晟达机械设计有限公司 一种高效的小区监控管理系统
CN109785240B (zh) * 2017-11-13 2021-05-25 中国移动通信有限公司研究院 一种低照度图像增强方法、装置及图像处理设备
CN108122210A (zh) * 2017-12-19 2018-06-05 长沙全度影像科技有限公司 一种基于Retinex和增强梯度的低照度车牌视频图像增强方法
CN108932467A (zh) * 2018-03-02 2018-12-04 于志国 跟踪区域图像分析平台
CN108492348A (zh) * 2018-03-30 2018-09-04 北京金山安全软件有限公司 图像处理方法、装置、电子设备及存储介质
CN110211049A (zh) * 2018-06-28 2019-09-06 京东方科技集团股份有限公司 基于Retinex理论的图像增强方法、装置和设备
CN109300101A (zh) * 2018-10-18 2019-02-01 重庆邮电大学 一种基于Retinex理论的多曝光图像融合方法
CN109345586A (zh) * 2018-11-02 2019-02-15 国网湖南省电力有限公司 基于紫外成像技术的电气设备放电特征提取方法
CN109472758B (zh) * 2018-11-20 2019-06-04 山东科技大学 一种地震剖面图像纹理细节增强方法
CN109859144B (zh) * 2019-02-22 2021-03-12 上海商汤智能科技有限公司 图像处理方法及装置、电子设备和存储介质
US20220051375A1 (en) * 2019-03-18 2022-02-17 Peking University Shenzhen Graduate School An inverse tone mapping method, system, device and computer readable medium
CN111028542B (zh) * 2019-04-10 2020-09-22 台州明创科技有限公司 定制对象预警控制装置
CN110796609B (zh) * 2019-07-31 2022-04-08 东南大学 基于尺度感知和细节增强模型的低光图像增强方法
CN110913195B (zh) * 2019-12-26 2022-10-04 深圳壹账通智能科技有限公司 白平衡自动调节方法、装置及计算机可读存储介质
CN113947535B (zh) * 2020-07-17 2023-10-13 四川大学 一种基于光照分量优化的低照度图像增强方法
CN112102204B (zh) * 2020-09-27 2022-07-08 苏州科达科技股份有限公司 图像增强方法、装置及电子设备
CN112308793A (zh) * 2020-10-21 2021-02-02 淮阴工学院 一种新型的非均匀光照图像对比度和细节增强方法
CN112508809B (zh) * 2020-11-27 2023-04-21 湖南傲英创视信息科技有限公司 低照度图像/视频增强方法及系统
CN114429426B (zh) * 2021-12-20 2022-08-16 哈尔滨理工大学 一种基于Retinex模型的低照度图像质量改善方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103413275A (zh) * 2013-07-26 2013-11-27 北京工业大学 基于梯度零范数最小化的Retinex夜间图像增强方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001197321A (ja) * 2000-01-12 2001-07-19 Nec Shizuoka Ltd カラー画像処理方法並びに画像処理装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103413275A (zh) * 2013-07-26 2013-11-27 北京工业大学 基于梯度零范数最小化的Retinex夜间图像增强方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
影像雾霾与不均匀光的变分校正方法研究;兰霞;《中国博士学位论文全文数据库信息科技辑》;20150615(第06期);摘要,第66-96页
进多尺度Retinex理论的低照度遥感影像增强方法;邵振峰等;《武汉大学学报(信息科学版)》;20150105;第40卷(第1期);第32-39页

Also Published As

Publication number Publication date
CN105654437A (zh) 2016-06-08

Similar Documents

Publication Publication Date Title
CN105654437B (zh) 一种对低照度图像的增强方法
CN107730475B (zh) 图像增强方法及系统
CN106780368B (zh) 一种基于前景模型的水下图像增强方法
US10134359B2 (en) Device or method for displaying image
CN106846282B (zh) 一种采用自适应校正的低照度图像增强方法
CN104240194B (zh) 一种基于抛物线函数的低照度图像增强算法
US7969480B2 (en) Method of controlling auto white balance
CN107578383A (zh) 一种低照度图像增强处理方法
CN109447910A (zh) 一种基于模糊理论的低照度彩色图像增强方法
JP2005341527A (ja) 階調補正装置、階調補正用プログラム、携帯端末機器及びパーソナルコンピュータ
CN105205794A (zh) 一种低照度图像的同步增强去噪方法
CN113989147A (zh) 基于直方图均衡和改进反锐化掩模的内窥镜图像增强方法
CN109658343B (zh) 结合颜色转换和自适应曝光的水下图像增强方法
CN111476732B (zh) 一种图像融合及去噪的方法及系统
JP6432399B2 (ja) 画像処理装置および画像処理方法
CN109639994B (zh) 嵌入式车载相机曝光时间动态调整方法
CN109587395A (zh) 一种基于图像处理和虚拟增强的防止眩光的系统及其实现方法
KR102296510B1 (ko) 영상 처리 장치, 영상 처리 방법, 및 영상 처리 프로그램
CN110969584A (zh) 一种低照度图像增强方法
Albu et al. One scan shadow compensation and visual enhancement of color images
CN107358578B (zh) 一种阴阳脸处理方法和装置
Linan et al. Video defogging based on adaptive tolerance
US20130286245A1 (en) System and method for minimizing flicker
CN110992287A (zh) 一种非均匀光照视频的清晰化方法
CN114140344A (zh) 一种低照度图像的图像增强方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190419

Termination date: 20191224