CN105654437B - 一种对低照度图像的增强方法 - Google Patents
一种对低照度图像的增强方法 Download PDFInfo
- 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
- image
- enhanced
- illumination
- low
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000012545 processing Methods 0.000 claims abstract description 45
- 230000002708 enhancing effect Effects 0.000 claims abstract description 16
- 230000003044 adaptive effect Effects 0.000 claims abstract description 8
- 238000005286 illumination Methods 0.000 claims description 89
- 238000001914 filtration Methods 0.000 claims description 21
- 230000002194 synthesizing effect Effects 0.000 claims description 13
- 238000012937 correction Methods 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims 1
- 230000000007 visual effect Effects 0.000 description 4
- 210000001525 retina Anatomy 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000003710 cerebral cortex Anatomy 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10016—Video; Image sequence
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10024—Color image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image 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*γi,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。
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通道的灰度值。
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 (33)
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 | 北京华易明新科技有限公司 | 一种低照度图像增强处理方法 |
CN107656454A (zh) * | 2017-09-21 | 2018-02-02 | 深圳市晟达机械设计有限公司 | 一种高效的小区监控管理系统 |
CN107443403B (zh) * | 2017-09-21 | 2019-11-01 | 上海昕鼎网络科技有限公司 | 一种用于巡逻的机器人 |
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模型的低照度图像质量改善方法 |
CN115767059A (zh) * | 2022-11-04 | 2023-03-07 | 合肥城市学院 | 光照不均情况下的白平衡校正方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413275A (zh) * | 2013-07-26 | 2013-11-27 | 北京工业大学 | 基于梯度零范数最小化的Retinex夜间图像增强方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001197321A (ja) * | 2000-01-12 | 2001-07-19 | Nec Shizuoka Ltd | カラー画像処理方法並びに画像処理装置 |
-
2015
- 2015-12-24 CN CN201510996197.3A patent/CN105654437B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413275A (zh) * | 2013-07-26 | 2013-11-27 | 北京工业大学 | 基于梯度零范数最小化的Retinex夜间图像增强方法 |
Non-Patent Citations (2)
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) | 一种对低照度图像的增强方法 | |
WO2019091270A1 (zh) | 图像增强方法及系统 | |
EP2852152B1 (en) | Image processing method, apparatus and shooting terminal | |
KR102234092B1 (ko) | 이미지의 역 톤 매핑을 위한 방법 | |
EP2833317B1 (en) | Image display device and/or method therefor | |
US9135681B2 (en) | Image chroma noise reduction | |
CN104240194B (zh) | 一种基于抛物线函数的低照度图像增强算法 | |
CN107578383A (zh) | 一种低照度图像增强处理方法 | |
IES20080497A2 (en) | Image processing method and apparatus | |
US9830690B2 (en) | Wide dynamic range imaging method | |
CN108830800A (zh) | 一种暗光场景下图像的亮度提升增强方法 | |
CN113989147A (zh) | 基于直方图均衡和改进反锐化掩模的内窥镜图像增强方法 | |
CN111968065A (zh) | 一种亮度不均匀图像的自适应增强方法 | |
KR20160049954A (ko) | 영상 처리 장치, 영상 처리 방법, 및 영상 처리 프로그램 | |
CN104766285A (zh) | 一种水下降质图像的自适应增强方法 | |
RU2622095C1 (ru) | Способ улучшения цифровых цветных изображений | |
Albu et al. | One scan shadow compensation and visual enhancement of color images | |
JP2008511048A (ja) | 画像処理方法、及び画像処理用コンピュータソフトウェア | |
CN115809966A (zh) | 一种低照度图像的增强方法和系统 | |
JP5410378B2 (ja) | 映像信号補正装置および映像信号補正プログラム | |
WO2011021518A1 (ja) | 画像処理方法、画像処理装置及びプログラム | |
CN114037641A (zh) | 一种低照度图像增强方法、装置、设备及介质 | |
CN110992287B (zh) | 一种非均匀光照视频的清晰化方法 | |
Terai et al. | Color image contrast enhancement by retinex model | |
KR101514152B1 (ko) | 특이값 분해를 이용한 화질 개선 방법 및 장치 |
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 |