CN101951523B - 一种自适应彩色图像处理方法及系统 - Google Patents

一种自适应彩色图像处理方法及系统 Download PDF

Info

Publication number
CN101951523B
CN101951523B CN201010290111A CN201010290111A CN101951523B CN 101951523 B CN101951523 B CN 101951523B CN 201010290111 A CN201010290111 A CN 201010290111A CN 201010290111 A CN201010290111 A CN 201010290111A CN 101951523 B CN101951523 B CN 101951523B
Authority
CN
China
Prior art keywords
image
function
luminance
highlights
overall situation
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
CN201010290111A
Other languages
English (en)
Other versions
CN101951523A (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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201010290111A priority Critical patent/CN101951523B/zh
Publication of CN101951523A publication Critical patent/CN101951523A/zh
Application granted granted Critical
Publication of CN101951523B publication Critical patent/CN101951523B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Processing (AREA)

Abstract

本发明公开了一种自适应彩色图像处理方法及系统,其方法为:读取一帧待处理图像;利用非线性反S型函数对待处理图像的亮度分量进行全局自适应变换,得到全局增强亮度图像;依据全局增强亮度图像中的中心像素与邻域像素的亮度的对应关系,对全局增强亮度图像的局部对比度进行调整,获得局部增强亮度图像;利用局部增强亮度图像与待处理图像的亮度分量两者之间的比值,对待处理图像进行线性颜色恢复,获取彩色增强图像。通过本发明的方法与系统,能够自适应地调节待处理图像亮度,同时增强暗区域和亮区域的对比度,且能够有效地增强待处理图像的局部对比度,并抑制图像边缘产生的过增强现象,且较好地保持了图像的颜色,取得更好的增强效果。

Description

一种自适应彩色图像处理方法及系统
技术领域
本发明涉及图像处理领域,更具体的说,涉及一种自适应彩色图像处理方法及系统。
背景技术
目前在大多数情况下,在彩色图像的获取过程中,由于成像传感器的动态范围有限导致采集的图像局部曝光不足、色彩暗淡、对比度低及细节损失,这极大限制和影响了后续图像分析与理解的研究工作。因此,有必要对采集到的图像进行增强处理,来改善图像的视觉效果、突出图像的特征。现有技术中常见的灰度图像增强技术有伽玛校正、自动色阶(autolevels)、USM锐化(UnSharp Masking,USM)、直方图均衡化等全局增强算法,以及局部直方图均衡化、同态滤波等局部增强算法。但是,由于彩色图像具有颜色信息,现阶段研究人员尚未完全了解大脑的颜色处理机制,因此,彩色图像的处理远比灰度图像复杂,很多针对灰度图像的处理方法不能直接推广到彩色图像处理。因此,在研究HVS(Human Visual System,人眼视觉系统)特性的基础上,已经提出了许多彩色图像增强算法。现有技术中主要包括两类基于视觉感知的彩色图像增强方法:第一类,从人眼对颜色的感知特性出发,提出了基于颜色恒常性的Retinex算法和ACE(Automatic Color Equalization,自动颜色均衡化)算法。其中,Retinex算法是一类应用最为广泛的方法。Retinex算法能够同时实现图像的全局和局部对比度增强、动态范围压缩以及基于灰度假设(gray world assumption)的颜色恒常性。ACE算法综合了GW(Gray World)和WP(White Patch)颜色校正机制。第二类,基于HVS的全局和局部自适应特性,提出了仿生(bio-inspired)彩色图像增强方法,用于增强图像的全局和局部对比度,而不考虑颜色恒常性的问题。
但是,采用上述现有技术中基于视觉感知的方法对彩色图像进行处理,仍然存在以下缺点:其一,采用Retinex算法,当图像整体或者局部区域违背灰度假设时,处理后的颜色将变为灰色调,从而产生灰化(graying out)效应,此外,沿图像边缘会产生Halo效应;其二,采用第二类中的一种仿生彩色图像增强方法,在全局自适应亮度调整过程中利用幂指数小于1的幂函数,虽然能够增强了图像中暗区域的对比度,但是降低了亮区域的对比度;而且利用双边滤波计算邻域像素的加权平均值,仍不可避免地在图像边缘处产生的过增强现象。
发明内容
有鉴于此,本发明提供了一种自适应彩色图像处理方法及系统,解决了现有技术中在对图像进行处理时,造成的图像亮区域对比度降低、图像边缘处产生的过增强现象、以及沿图像边缘产生的Halo效应和灰化效应的问题。
为实现上述目的,本发明提供如下技术方案:
一种自适应彩色图像处理方法,包括:
读取一帧待处理图像;
利用非线性反S型函数对所述待处理图像的亮度分量进行全局自适应变换,得到全局增强亮度图像;
获取全局增强亮度图像,利用增强函数
Figure GDA0000130445400000021
依据所述全局增强亮度图像中的中心像素与邻域像素的亮度的对应关系,对所述全局增强亮度图像的局部对比度进行调整,获得局部增强亮度图像;
其中,Ω(x)是以像素x为中心的小邻域,d(x,y)为所述待处理图像的中心像素x与邻域像素y的欧氏距离,L(x)-L(y)为中心像素x和邻域像素y的亮度差值u,r(u)为饱和函数,用于确定亮度差值u的调整范围;
利用所述局部增强亮度图像与所述待处理图像的亮度分量两者之间的比值,对所述待处理图像进行线性颜色恢复,获取彩色增强图像。
优选的,所述利用非线性反S型函数对所述待处理图像的亮度分量进行全局自适应变换,得到全局增强亮度图像的过程具体为:
确定非线性映射的区间范围,将所述亮度分量的灰度范围映射到所述区间范围内,获取灰度级压缩的图像;
利用非线性反S型函数对所述灰度级压缩的图像进行非线性灰度级变换,获得反S型变换图像;
将所述反S型变换图像的值归一化,获取全局自适应变换后的所述全局增强亮度图像。
优选的,所述依据所述对应关系对所述全局增强亮度图像的局部对比度进行调整,获得局部增强亮度图像的过程为:
获取全局增强亮度图像;
利用所述增强函数对所述全局增强亮度图像进行处理,得到所述局部增强亮度图像。
优选的,获取所述增强函数的过程包括:
确定所述中心像素与所述邻域像素的亮度差值;
获取所述亮度差值在饱和函数范围内的对应值;
确定所述中心像素与所述邻域像素的空间距离;
获取以所述中心像素为中心的小邻域内,所述亮度差值在饱和函数范围内的对应值与所述空间距离之间的比值之和,得到所述增强函数。
优选的,对所述待处理图像进行线性颜色恢复的过程为:
利用局部增强亮度图像与所述待处理图像的亮度分量两者之间的比值,获得颜色复原函数;
利用获取到的所述颜色复原函数等比例乘以所述待处理图像中的各颜色分量,得到所述彩色增强图像。
优选的,所述非线性映射的区间范围为:[smin,smax],其中,smin为s(vmin),smax为s(vmax),s(v)为S型函数,vmin值的取值范围为:-6~-3,vmax值的取值范围为:1~6。
优选的,所述饱和函数的取值范围为0~1,所述饱和函数的斜率为-3。
一种自适应彩色图像处理系统,包括:
读取单元,用于读取一帧待处理图像;
全局亮度调整单元,用于利用非线性反S型函数对所述待处理图像的亮度分量进行全局自适应变换,得到全局增强亮度图像;
局部亮度调整单元,用于依据所述全局增强亮度图像中的中心像素与邻域像素的亮度的对应关系,对所述全局增强亮度图像的局部对比度进行调整,获得局部增强亮度图像;
颜色恢复单元,用于利用所述局部增强亮度图像与所述待处理图像的亮度分量两者之间的比值,对所述待处理图像进行线性颜色恢复,获取彩色增强图像。
优选的,所述全局亮度调节单元中主要包括:
灰度级压缩单元,用于确定非线性映射的区间范围,将所述亮度分量的灰度范围映射到所述区间范围内,获取灰度级压缩的图像;
非线性映射单元,用于利用非线性反S型函数对所述灰度级压缩的图像进行非线性灰度级变换,获得反S型变换图像;
归一化单元,用于将所述反S型变换图像的值归一化,获取全局自适应变换后的所述全局增强亮度图像。
优选的,所述局部亮度调整单元中主要包括:
增强函数获取单元,用于获取所述亮度差值在饱和函数范围内的对应值,以及确定所述中心像素与所述邻域像素的亮度差值和空间距离,并依据获取以所述中心像素为中心的小邻域内,所述亮度差值在饱和函数范围内的对应值与所述空间距离之间的比值之和,得到所述增强函数;
局部调整单元,用于利用所述增强函数对所述全局增强亮度图像进行处理,得到所述局部增强亮度图像。经由上述的技术方案可知,与现有技术相比,本发明公开了一种自适应彩色图像处理方法及系统,基于HVS的亮度适应机理,利用神经元的非线性反S型函数对彩色降质视频中的待处理图像进行全局自适应变换,并获取全局增强亮度图像;然后再基于HVS的侧抑制机理对全局增强亮度图像的局部对比度进行增强,并获取局部增强亮度图像;利用所述局部增强亮度图像与待处理图像的亮度分量两者之间的比值,对所述待处理图像进行线性颜色恢复,获取彩色增强图像。通过本发明的方法与系统,能够自适应地调节待处理图像的亮度,并同时调节暗区域和亮区域的对比度和亮度,而且能够有效地增强待处理图像的局部对比度,并抑制待处理图像边缘产生的过增强现象,且较好地保持了图像的颜色。
此外也不会产生灰化效应和Halo效应,在图像处理的过程中取得更好的增强效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本发明实施例一公开的一种自适应彩色图像处理方法的流程图;
图2为非线性S型函数;
图3为本发明实施例二公开的一种自适应彩色图像处理方法的流程图;
图4为用于全局亮度自适应变换的非线性映射曲线;
图5为增强函数中所用到的饱和函数;
图6a为本发明公开的方法与现有技术中的方法进行比较时的原图的亮度分量直方图;
图6b为采用现有技术对图6a进行彩色增强处理后的亮度分量直方图;
图6c为采用本发明对图6a进行彩色增强处理后的亮度分量直方图;
图7a为本发明公开的方法与现有技术中的方法进行比较时的原图的R、G、B分量直方图;
图7b为采用现有技术对图7a进行彩色增强处理后的R、G、B分量直方图;
图7c为采用本发明对图7a进行彩色增强处理后的R、G、B分量直方图;
图8为本发明实施例二公开的一种自适应彩色图像处理的系统结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明公开了一种自适应彩色图像处理方法及系统,通过利用非线性反S型函数对降质彩色视频中的待处理图像进行全局自适应变换,获取全局增强亮度图像;然后再对所述全局增强亮度图像的局部对比度进行增强,并获取增强后的局部增强亮度图像;利用所述局部增强亮度图像与所述待处理图像的亮度分量两者之间的比值,对所述待处理图像进行线性颜色恢复,获取彩色增强图像。下面对本发明的方法下面将进行详细说明。
实施例一:
如图1所示,本发明的方法主要包括以下步骤:
步骤S101,读取一帧待处理图像。
执行步骤S101,从彩色降质视频序列中读取一帧视频图像I(x)=(IR(x),IG(x),IB(x))T,该视频图像即为待处理图像,是具有R、G、B三个分量的彩色图像。其中,该待处理图像宽为W个像素、高为H个像素,由于图像是二维空间,因此上式中的x代表向量,x也可用(m,n)表示,其中,0≤m≤W-1,0≤n≤H-1,m、n、W、H均为非负整数。步骤S102,利用非线性反S型函数对所述待处理图像的亮度分量进行全局自适应变换,得到全局增强亮度图像La(x)。
人眼具有对高动态范围和明暗亮度适应的特性,虽然能够响应的照度范围高达13个对数单位,但是并不能同时在这个范围内工作;对于人眼所具有的明暗亮度适应主要指,人眼能够根据环境光照的变化,自动调节对亮度的适应过程,而人眼的适应明暗条件变化的过程称为亮度适应。其中,从亮到暗的过程称为暗适应,从暗到亮的过程称为亮适应。
基于上述HVS的亮度适应机理,执行步骤S102,利用非线性反S型函数对待处理图像中的亮度进行全局自适应变换,即利用非线性反S型函数进行灰度级映射,从而完成全局自适应调节待处理图像亮度的过程,得到全局增强亮度图像La(x)。
在本发明中只对待处理图像的亮度分量L(x)进行处理,在HSV颜色空间中,彩色图像的亮度分量L(x)具体为:
L(x)=max{IR(x),IG(x),IB(x)};    (1)
在(1)式中,所述IR(x)、IG(x)和IB(x)分别代表颜色空间中点x处像素的R(Red,红色)、G(Green,绿色)和B(Blue,蓝色)。因此可以避免现有技术中分别处理R、G、B通道易导致颜色失真,且增加计算量带来的功耗问题。即仅对待处理图像的亮度分量进行处理,降低了执行该步骤时产生的功耗。
本发明公开的实施例一中,在HSV颜色空间的亮度分量是计算各颜色分量的最大值,但是,本发明并不仅限于此,在HSI颜色空间的亮度分量是计算3个颜色分量的平均值,也可以用PCA变换后的第一主成分。另外,也可以采用其他颜色空间转换方法。
需要说明的是,S型函数定义为:
s ( v ) = 1 1 - e - αv - - - ( 2 )
在(2)式中,α为控制函数s(v)的增长速率,S型函数的值域为(0,1)。而S型函数的示意图,如图2所示,图中横轴为v,纵轴为s。
步骤S103,依据所述全局增强亮度图像La(x)中的中心像素与邻域像素的亮度的对应关系,对所述全局增强亮度图像La(x)的局部对比度进行调整,获得局部增强亮度图像Len(x)。
经过步骤S102的全局自适应变换之后的待处理图像,即全局增强亮度图像La(x),增强了亮区域和暗区域的对比度,但是损失了中间亮度区域对比度。研究表明,人眼对全局对比度并不敏感,而对局部对比度更为敏感。HVS的中心/邻域侧抑制特性在视觉中的作用是增强反差,使人眼观察物体的边角和轮廓时会提高视敏度。
依据中心像素与邻域像素的亮度的对应关系,即当所述中心像素的亮度高于所述邻域像素的平均亮度时,增强所述中心像素的亮度;当所述中心像素的亮度低于所述邻域像素的平均亮度时,降低所述中心像素的亮度。
而在本发明公开的实施例一中,通过指数形式最终获取局部增强亮度图像Len(x)可以表示为:
Len(x)=La(x)E(x)    (3)
其中,E(x)为增强函数。若E(x)>1,则降低该中心像素的亮度;反之,若0<E(x)<1,则增加中心像素的亮度。
为避免待处理图像的边缘处产生Halo效应和过增强现象,本发明中利用ACE算法,定义增强函数E(x)为:
E ( x ) = Σ y ∈ Ω , y ≠ x r { L ( x ) - L ( y ) } d ( x , y ) - - - ( 4 )
在(4)式中,Ω(x)是以像素x为中心的小邻域,在本发明优选的实施例一中Ω(x)所取的邻域为7×7邻域,但是,本发明所公开的实施例并不仅限于此;d(x,y)为所述待处理图像的中心像素x与邻域像素y的欧氏距离,作为邻域像素的距离权重函数;L(x)-L(y)为中心像素x和邻域像素y的差值,用于模拟人眼的侧抑制机制;在中心像素与邻域像素的亮度进行比较的过程中,如图5所示(图中横轴为u,纵轴r),对于上述邻域中的差值L(x)-L(y)的范围用饱和函数r(u)进行调整,具体表示为:
r ( u ) = 1 , u &le; - 1 2 k - ku + 0.5 , | u | < 1 2 k 0 , u &GreaterEqual; 1 2 k - - - ( 5 )
在(5)式中,k为经验值,本发明的实施例中取值为-3,但是并不仅限于此。由上式中可知,k控制饱和函数r(u)的斜率,k的取值越大,图像局部对比度增强越大。该r(u)的值域为[0,1],采用约束饱和函数r(u)上下界的方式,能够达到抑制图像边缘处产生的Halo效应和过增强现象。
此外,通过(4)式与(5)式可知,获取增强函数E(x)的具体的过程为:
首先,确定所述中心像素与邻域像素的亮度差值u,以及该中心像素与所述邻域像素的空间距离;获取所述亮度差值u在饱和函数r(u)范围内的对应值;然后获取以该中心像素为中心的小邻域内,亮度差值u在饱和函数r(u)范围内的对应值与所述空间距离之间的比值之和;最后得到该增强函数E(x)。
步骤S104,利用所述局部增强亮度图像Len(x)与所述待处理图像的亮度分量L(x)两者之间的比值,对所述待处理图像进行线性颜色恢复处理,获取彩色增强图像。
通过执行上述步骤S101至S103,对待处理图像进行了全局自适应变换和局部对比度增强处理后,得到最终的局部增强亮度图像Len(x)。最后,再执行步骤S104恢复亮度增强后的图像的颜色。
其中,通过上述的全局自适应变换和局部对比度增强,获得增强后的局部增强亮度图像Len(x)。定义颜色复原函数为:
κ(x)=Len(x)/L(x)    (6)
将κ(x)等比例地乘以3个颜色分量,恢复增强图像的颜色信息,即
I′c(x)=κ(x)Ic(x)    (7)
在(7)式中,c={R,G,B}分别表示R、G、B颜色通道;IR(x)、IG(x)、IB(x)和I′R(x)、I′G(x)、I′B(x)分别表示增强前后的R、G、B颜色分量值。由于线性操作能够使R、G、B颜色分量的比例关系保持不变,因此,能够保持原图像的色度。
通过本发明实施例公开的方法,能够保证在进行图像处理的过程中,不仅可以自适应地调节图像的亮度,同时增强暗区域和亮区域的对比度;同时,还能有效地增强图像的局部对比度,并抑制边缘处的过增强现象。此外,也不会产生灰化效应和Halo效应。
需要说明的是,本发明公开的上述步骤S101至S104的执行过程,只是针对彩色降质视频序列中的一帧视频图像进行处理的过程,所以在执行完一帧待处理图像的处理之后,获取下一帧待处理图像进行循环处理,直至该彩色降质视频序列中所有待处理图像帧处理完毕,最终获取一组完整的、彩色增强的视频图像数据。
实施例二:
在上述本发明公开的实施例的基础上,本发明还公开了一种自适应彩色图像处理方法,请参阅附图3,主要包括以下步骤:
步骤S101,读取一帧待处理图像。
步骤S1021,确定非线性映射的区间范围,将所述亮度分量L(x)的灰度范围映射到所述区间范围内,获取灰度级压缩的图像L′(x)。
在S型函数中,动态确定灰度映射区间vmin,vmax]。vmin的取值由下列的线性分段函数进行表示。
v min = - 6 , M &le; 0.3 0.04 ( M - 0.3 ) - 6 , 0.3 < M &le; 0.6 - 3 , M > 0.6 - - - ( 8 )
在(8)式中,M为待处理图像的平均亮度。根据图像的平均亮度,自适应地调节暗区域。而vmax的取值为1至6之间,在本发明优选的实施例一中vmax的取值为3。
将亮度分量L(x)的灰度范围从[0,1]线性映射到非线性映射区间范围[smin,smax]中,进行灰度映射获取图像L′(x)的具体过程为:
L′(x)=(smax-smin)L(x)+Smin;    (9)
在(9)式中,smin=s(vmin),smax=s(vmax)。
在执行步骤S1021中,根据图像的平均亮度M,自适应地调节暗区域。其中,随着vmin取值的减小,待处理图像的暗区域的对比度和亮度的增强程度逐渐提升,图像的整体亮度逐渐提高。vmin是根据图像的平均亮度自适应确定的。因此,对于vmin的取值根据实际情况的不同而有所不同。
如图4所示,为vmin=-3、vmin=-4.5、vmin=-6时的全局亮度自适应变换的非线性映射曲线,从图中可知,vmin值越小,对待处理图像的暗区域的提升越明显。图中横轴为L,纵轴为La
步骤S1022,利用非线性反S型函数对所述灰度级压缩的图像L′(x)进行非线性灰度级变换,获得反S型变换图像L″(x),所述获取图像L″(x)具体为:
L &prime; &prime; ( x ) = - 1 &alpha; ln ( 1 L &prime; ( x ) - 1 ) - - - ( 10 )
式中,α的取值为1,但并不仅限于1。
执行步骤S1022,利用步骤S1021中确定的非线性反S型函数的非线性映射曲线中的一段,然后利用对图像L′(x)进行灰度级变换。
步骤S1023,将所述反S型变换图像L″(x)的值归一化至[0,1]中,获取全局自适应变换后的所述全局增强亮度图像La(x)。
在基于HVS亮度适应机理的条件下,根据待处理图像的明暗程度,利用非线性反S型函数对该待处理图像的亮度进行非线性调节,增强了该待处理图像的暗区域和亮区域的对比度,同时压缩了该待处理图像的动态范围,达到使图像符合人眼直接观察场景的目的。
步骤S103,依据所述全局增强亮度图像La(x)中的中心像素与邻域像素的亮度的对应关系,对所述全局增强亮度图像La(x)的局部对比度进行调整,获得局部增强亮度图像Len(x)。
步骤S104,利用所述局部增强亮度图像Len(x)与所述待处理图像的亮度分量L(x)两者之间的比值,对所述待处理图像进行线性颜色恢复,获取彩色增强图像。
通过上述的步骤可以在执行完一帧待处理图像的处理之后,经过判断获取下一帧待处理图像进行循环处理,直至该彩色降质的视频序列中所有待处理图像帧处理完毕,最终获取一组完整的、还原色彩的视频图像数据。
需要说明的是,上述本发明公开的对图像进行处理的过程,不仅适用于彩色图像,也适用于灰度图像。
此外,为进一步验证本发明实施例公开的方法对图像处理的效果,本发明从定性和定量两个方面,对本发明的方法和现有技术中的一种仿生彩色图像增强方法进行了分析比较。
如图6a所示,为一幅傍晚拍摄的低对比度图像的直方图,由于成像传感器的动态范围有限,背景照度较高造成前景曝光不足,暗区域的局部细节严重损失。图6b和图6c分别为现有技术和本发明的方法对图像增强的相应亮度图像的直方图。从图6a、图6b和图6c的直方图中可以看出,原图像的像素主要集中在直方图的两端;现有技术增强后的图像像素整体偏向直方图的亮端;而用本发明方法增强后的图像提高了暗区域的亮度,而亮区域的亮度没有明显变化。
图7a、图7b和图7c为局部区域放大图像的R、G、B颜色分量直方图,图7a、图7b和图7c中,R、G、B分别表示三个颜色亮度分量。
为了使显示效果较为明显,截去了70%灰度值为0的像素。直方图相似度测量两个直方图分布之间的相似程度。为了算法性能的定量评价与比较,本发明使用基于直方图相似度的4个评测准则,其中,h=(hO,·,hk,·,hN-1)T和h′=(hO,·,h′k,·,h′N-1)T分别为处理前后图像直方图中各灰度级k的出现统计值,分别为:
相关系数(Correlation):
d correl ( h &prime; , h ) = &Sigma; k ( h k &prime; - h &prime; &OverBar; ) ( h k - h &OverBar; ) &Sigma; k ( h k &prime; - h &prime; &OverBar; ) 2 &Sigma; k ( h k - h &OverBar; ) 2 - - - ( 11 )
其中, h &OverBar; = 1 N &Sigma; k h k
N表示灰度级。dcorrel越大,匹配度越高。完全匹配时,dcorrel为1,而完全不匹配时,dcorrel为-1。
卡方度量(Chi-square):
d chi - sqr ( h &prime; , h ) = &Sigma; k ( h k &prime; - h k ) 2 h k &prime; + h k - - - ( 12 )
ddchi-sqr越小,匹配度越高。完全匹配时,dchi-sqr为0。
直方图交集(Intersection):
d inter sec t ( h &prime; , h ) = &Sigma; k min ( h k &prime; , h k ) - - - ( 13 )
dintersect越大,匹配度越高。若将直方图做归一化处理,即为灰度级出现的概率,则完全匹配时,dintersect为1;而完全不匹配时,dintersect为0。
巴氏距离(Bhattacharyya distance):
d B ( h &prime; , h ) = 1 - &Sigma; k h k &prime; h k &Sigma; k h k &prime; &CenterDot; &Sigma; k h k - - - ( 14 )
dB越小,匹配度越高。完全匹配时,dB为0;而完全不匹配,dB为1。
表1列出了图6a与图6b(现有技术的方法)以及图6a与图6c(本发明的方法)的相似度,可以看出,本发明的方法更好地保持了直方图的形状。
表1:
Figure GDA0000130445400000126
此外,通过采用本发明的方法在大量室内外图像上所做的实验,取得了较好的效果。更进一步证明了采用本发明的方法增强了图像的局部对比度和细节,真实地展现了暗区域的颜色,提高了可视性,使增强图像符合人眼直接观察场景的视觉感知。
通过上述本发明实施例方法的描述,所属领域的技术人员可以清楚地了解到,本发明的方法可以依托于通用的硬件完成,因此,在上述公开的方法的基础上,本发明的实施例公开了一种自适应彩色图像处理的系统。
请参阅附图8,在本发明公开的实施例中,该系统主要包括:读取单元901、全局亮度调整单元902、局部亮度调整单元903和颜色恢复单元904。
读取单元901,用于从彩色降质的视频序列中读取一帧待处理图像。
全局亮度调整单元902,用于利用非线性反S型函数对所述待处理图像的亮度分量进行全局自适应变换,得到全局增强亮度图像。
局部亮度调整单元903,用于依据所述全局增强亮度图像中的中心像素与邻域像素的亮度的对应关系,对所述全局增强亮度图像的局部对比度进行调整,获得局部增强亮度图像。
颜色恢复单元904,用于利用所述局部增强亮度图像与所述待处理图像的亮度分量两者之间的比值,对所述待处理图像进行线性颜色恢复,获取彩色增强图像。
在本发明公开的实施例,通过读取单元901获取一帧待处理图像;然后,在全局亮度调整单元902中对该待处理图像进行全局自适应变换,获取全局增强亮度图像;然后,在局部亮度调整单元903中对该全局增强亮度图像进行局部对比度增强,获取局部增强亮度图像;然后,在颜色恢复单元904中,对经过上述全局亮度调整单元902与局部亮度调整单元903后获取的局部增强亮度图像进行颜色恢复,获取最终的增强图像,即彩色增强图像。
需要说明的是,如图8所示,在上述本发明公开的系统还包括一个判断单元905。
判断单元905,用于判断是否还存下一帧待处理图像,如果是,则返回读取单元;如果否,则停止图像处理过程。
在利用本发明公开的实施例二的系统中进行图像处理的过程,为在执行完一帧待处理图像的处理之后,经过判断单元905判断后获取下一帧待处理图像进行循环处理,直至该彩色降质的视频序列中所有待处理图像帧处理完毕,最终获取一组完整的、还原色彩的视频图像数据。
需要说明的是,在全局亮度调整单元902中主要包括:灰度级压缩单元9021、非线性映射单元9022和归一化单元9023,并在以上单元中进行全局自适应变换。
灰度级压缩单元9021,用于确定非线性映射的区间范围,将所述亮度分量的灰度范围映射到所述区间范围内,获取灰度级压缩的图像。
非线性映射单元9022,用于利用线性反S型函数对所述灰度级压缩的图像进行非线性灰度级变换,获得反S型变换图像。
归一化单元9023,用于将所述反S型变换图像的值,进行归一化处理,获取自适应调整后的所述全局增强亮度图像。
此外,需要说明的是,在局部亮度调整单元903中主要包括:增强函数获取单元9031和局部调整单元9032。
增强函数获取单元9031,用于获取所述亮度差值在饱和函数范围内的对应值,以及确定所述中心像素与所述邻域像素的亮度差值和空间距离,并依据获取以所述中心像素为中心的小邻域内,所述亮度差值在饱和函数范围内的对应值与所述空间距离之间的比值之和,得到所述增强函数。
局部调整单元9032,用于利用所述增强函数对所述全局增强亮度图像进行处理,得到所述局部增强亮度图像。
在上述单元对待处理图像进行处理的过程中,所依据的原理或机理、公式或函数对待处理图像进行处理的过程,与上述本发明公开的方法实施例中的处理过程相同,这里不再赘述。
因此,通过本发明实施例公开的系统,能够保证在进行图像处理的过程中,不仅可以自适应地调节图像的亮度,可以同时增强暗区域和亮区域的对比度,能有效地校正待处理图像的颜色;同时,还能有效地增强图像的局部对比度,并抑制边缘处的过增强现象,且不会产生灰化效应和Halo效应。
综上所述:
通过本发明各个实施例公开的方法及系统,通过本发明的方法与系统,能够自适应地调节待处理图像的亮度,并同时调节暗区域和亮区域的对比度和亮度,而且能够有效地增强待处理图像的局部对比度,并抑制待处理图像边缘产生的过增强现象,且较好地保持了图像的颜色。此外也不会产生灰化效应和Halo效应,在图像处理的过程中取得更好的增强效果。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块,或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (9)

1.一种自适应彩色图像处理方法,其特征在于,包括:
读取一帧待处理图像;
利用非线性反S型函数对所述待处理图像的亮度分量进行全局自适应变换,得到全局增强亮度图像;
获取全局增强亮度图像,利用增强函数依据所述全局增强亮度图像中的中心像素与邻域像素的亮度的对应关系,对所述全局增强亮度图像的局部对比度进行调整,获得局部增强亮度图像;
其中,Ω(x)是以像素x为中心的小邻域,d(x,y)为所述待处理图像的中心像素x与邻域像素y的欧氏距离,L(x)-L(y)为中心像素x和邻域像素y的亮度差值u,r(u)为饱和函数,用于确定亮度差值u的调整范围;
利用所述局部增强亮度图像与所述待处理图像的亮度分量两者之间的比值,对所述待处理图像进行线性颜色恢复,获取彩色增强图像。
2.根据权利要求1所述的方法,其特征在于,所述利用非线性反S型函数对所述待处理图像的亮度分量进行全局自适应变换,得到全局增强亮度图像的过程具体为:
确定非线性映射的区间范围,将所述亮度分量的灰度范围映射到所述区间范围内,获取灰度级压缩的图像;
利用非线性反S型函数对所述灰度级压缩的图像进行非线性灰度级变换,获得反S型变换图像;
将所述反S型变换图像的值归一化,获取全局自适应变换后的所述全局增强亮度图像。
3.根据权利要求2所述的方法,其特征在于,获取所述增强函数的过程包括:
确定所述中心像素与所述邻域像素的亮度差值;
获取所述亮度差值在饱和函数范围内的对应值;
确定所述中心像素与所述邻域像素的空间距离;
获取以所述中心像素为中心的小邻域内,所述亮度差值在饱和函数范围内的对应值与所述空间距离之间的比值之和,得到所述增强函数。
4.根据权利要求1所述的方法,其特征在于,对所述待处理图像进行线性颜色恢复的过程为:
利用局部增强亮度图像与所述待处理图像的亮度分量两者之间的比值,获得颜色复原函数;
利用获取到的所述颜色复原函数等比例乘以所述待处理图像中的各颜色分量,得到所述彩色增强图像。
5.根据权利要求1所述的方法,其特征在于,所述非线性映射的区间范围为:[smin,smax],其中,smin为s(vmin),smax为s(vmax),s(v)为S型函数,vmin值的取值范围为:-6~-3,vmax值的取值范围为:1~6。
6.根据权利要求5所述的方法,其特征在于,所述饱和函数的取值范围为0~1,所述饱和函数的斜率为-3。
7.一种自适应彩色图像处理系统,其特征在于,包括:
读取单元,用于读取一帧待处理图像;
全局亮度调整单元,用于利用非线性反S型函数对所述待处理图像的亮度分量进行全局自适应变换,得到全局增强亮度图像;
局部亮度调整单元,用于依据所述全局增强亮度图像中的中心像素与邻域像素的亮度的对应关系,对所述全局增强亮度图像的局部对比度进行调整,获得局部增强亮度图像;
颜色恢复单元,用于利用所述局部增强亮度图像与所述待处理图像的亮度分量两者之间的比值,对所述待处理图像进行线性颜色恢复,获取彩色增强图像。
8.根据权利要求7所述的系统,其特征在于,所述全局亮度调节单元中主要包括:
灰度级压缩单元,用于确定非线性映射的区间范围,将所述亮度分量的灰度范围映射到所述区间范围内,获取灰度级压缩的图像;
非线性映射单元,用于利用非线性反S型函数对所述灰度级压缩的图像进行非线性灰度级变换,获得反S型变换图像;
归一化单元,用于将所述反S型变换图像的值归一化,获取全局自适应变换后的所述全局增强亮度图像。
9.根据权利要求7所述的系统,其特征在于,所述局部亮度调整单元中主要包括:
增强函数获取单元,用于获取所述亮度差值在饱和函数范围内的对应值,以及确定所述中心像素与所述邻域像素的亮度差值和空间距离,并依据获取以所述中心像素为中心的小邻域内,所述亮度差值在饱和函数范围内的对应值与所述空间距离之间的比值之和,得到所述增强函数;
局部调整单元,用于利用所述增强函数对所述全局增强亮度图像进行处理,得到所述局部增强亮度图像。
CN201010290111A 2010-09-21 2010-09-21 一种自适应彩色图像处理方法及系统 Expired - Fee Related CN101951523B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201010290111A CN101951523B (zh) 2010-09-21 2010-09-21 一种自适应彩色图像处理方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201010290111A CN101951523B (zh) 2010-09-21 2010-09-21 一种自适应彩色图像处理方法及系统

Publications (2)

Publication Number Publication Date
CN101951523A CN101951523A (zh) 2011-01-19
CN101951523B true CN101951523B (zh) 2012-10-24

Family

ID=43454861

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010290111A Expired - Fee Related CN101951523B (zh) 2010-09-21 2010-09-21 一种自适应彩色图像处理方法及系统

Country Status (1)

Country Link
CN (1) CN101951523B (zh)

Families Citing this family (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102209246B (zh) * 2011-05-23 2013-01-09 北京工业大学 一种实时视频白平衡处理系统
CN103971330B (zh) * 2013-02-05 2017-10-20 腾讯科技(深圳)有限公司 图像增强方法及装置
CN104599238B (zh) * 2013-10-30 2017-09-12 腾讯科技(北京)有限公司 图像处理方法及装置
CN103745687B (zh) * 2013-12-04 2016-01-13 西安诺瓦电子科技有限公司 图像显示方法和图像显示装置
WO2015103739A1 (zh) * 2014-01-08 2015-07-16 富士通株式会社 增强图像对比度的装置、电子设备及方法
CN105338214B (zh) * 2014-08-13 2018-07-20 山东新北洋信息技术股份有限公司 图像处理方法和图像处理装置
CN104240197B (zh) * 2014-08-26 2017-12-15 浙江工商大学 一种保持对比度、颜色一致性和灰度像素特征的消色方法
CN104299208B (zh) * 2014-09-15 2017-10-13 深圳市云宙多媒体技术有限公司 一种自适应图像亮度的处理方法与装置
CN104392425B (zh) * 2014-12-04 2017-07-21 厦门美图移动科技有限公司 一种基于人脸的自动调整对比度的图像增强方法
CN104835128B (zh) * 2015-04-17 2017-05-24 南京大学 一种用于模糊视频文本信息增强的多光谱融合方法
KR20180015248A (ko) * 2015-06-09 2018-02-12 후아웨이 테크놀러지 컴퍼니 리미티드 비디오 인코딩 방법, 비디오 디코딩 방법, 비디오 인코더 및 비디오 디코더
CN105243646A (zh) * 2015-10-28 2016-01-13 上海大学 面部纹理特征增强方法
CN105744118B (zh) * 2016-02-01 2018-11-30 杭州当虹科技有限公司 一种基于视频帧自适应的视频增强方法及视频增强系统
CN106296620B (zh) * 2016-08-14 2019-06-04 遵义师范学院 一种基于直方图平移的色彩还原方法
CN106327443A (zh) * 2016-08-24 2017-01-11 电子科技大学 一种改进的遗传算法的夜间图像增强方法
CN106548465A (zh) * 2016-11-25 2017-03-29 福建师范大学 一种多光谱遥感影像的增强方法
CN107463895B (zh) * 2017-07-28 2019-09-20 中国科学院西安光学精密机械研究所 基于邻域向量pca的弱小损伤目标检测方法
CN107451974B (zh) * 2017-07-31 2020-06-02 北京电子工程总体研究所 一种高动态范围图像的自适应再现显示方法
CN109325922B (zh) * 2018-09-12 2022-03-11 深圳开阳电子股份有限公司 一种图像自适应增强方法、装置及图像处理设备
CN110995549B (zh) * 2019-03-19 2021-09-03 上海越安建筑安装工程有限公司 通信组件开关控制系统
CN112166598B (zh) * 2019-07-30 2022-03-11 深圳市大疆创新科技有限公司 图像处理方法、系统、可移动平台和存储介质
CN110473158A (zh) * 2019-08-14 2019-11-19 上海世茂物联网科技有限公司 一种车牌图像亮度的处理方法、装置及设备
CN111027399B (zh) * 2019-11-14 2023-08-22 武汉兴图新科电子股份有限公司 一种基于深度学习的遥感图像水面潜艇识别方法
CN112866802B (zh) * 2019-11-27 2022-03-15 深圳市万普拉斯科技有限公司 视频处理方法、装置、存储介质及计算机设备
CN111325680B (zh) * 2020-01-08 2023-04-25 深圳深知未来智能有限公司 一种具有抑制局部过曝的图像增亮方法
CN112541869A (zh) * 2020-12-07 2021-03-23 南京工程学院 一种基于matlab实现的Retinex图像去雾方法
CN114640836B (zh) * 2020-12-15 2023-07-25 深圳市奥拓电子股份有限公司 视频数据映射方法、led显示控制器及led显示系统
CN112383719B (zh) * 2021-01-13 2022-02-01 腾讯科技(深圳)有限公司 图像亮度的调整方法、装置、设备及可读存储介质
CN113096031B (zh) * 2021-03-17 2024-02-06 西安电子科技大学 高动态范围红外图像的压缩显示方法
CN113393391B (zh) * 2021-06-10 2022-05-20 维沃移动通信(杭州)有限公司 图像增强方法、图像增强装置、电子设备和存储介质
CN114298935B (zh) * 2021-12-27 2024-06-04 重庆港宇高科技开发有限公司 图像增强方法、装置及计算机可读存储介质
CN115937016B (zh) * 2022-10-31 2023-07-21 哈尔滨理工大学 一种保证图像细节的对比度增强方法
CN115512231B (zh) * 2022-11-15 2023-02-28 山东省国土空间数据和遥感技术研究院(山东省海域动态监视监测中心) 适用于国土空间生态修复的遥感解译方法
CN117201693B (zh) * 2023-11-01 2024-01-16 长春汽车工业高等专科学校 一种物联网图像压缩方法、装置、终端设备及介质

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100568279C (zh) * 2008-07-09 2009-12-09 北京航空航天大学 一种基于Retinex理论的快速彩色图像增强方法
JP2010193199A (ja) * 2009-02-18 2010-09-02 Victor Co Of Japan Ltd 画像処理装置および画像処理方法
CN101742339A (zh) * 2010-01-14 2010-06-16 中山大学 一种彩色图像增强的方法
CN101853490A (zh) * 2010-04-21 2010-10-06 中国科学院半导体研究所 一种基于人类视觉特性的仿生图像复原方法

Also Published As

Publication number Publication date
CN101951523A (zh) 2011-01-19

Similar Documents

Publication Publication Date Title
CN101951523B (zh) 一种自适应彩色图像处理方法及系统
CN108765336B (zh) 基于暗亮原色先验与自适应参数优化的图像去雾方法
US11127122B2 (en) Image enhancement method and system
CN113313657B (zh) 一种用于低光照图像增强的非监督学习方法和系统
CN105745914B (zh) 用于逆色调映射的方法和系统
CN103593830B (zh) 一种低照度视频图像增强方法
US9147238B1 (en) Adaptive histogram-based video contrast enhancement
CN106683080B (zh) 一种视网膜眼底图像预处理方法
CN108198155B (zh) 一种自适用色调映射方法及系统
CN104700376A (zh) 基于伽马校正和平滑滤波的图像直方图均衡化增强方法
CN104268843A (zh) 基于直方图修饰的图像自适应增强方法
CN106780417A (zh) 一种光照不均图像的增强方法及系统
CN103455979A (zh) 一种低照度视频图像增强方法
CN104200431A (zh) 图像灰度化的处理方法及处理装置
CN112116536A (zh) 一种低光照图像增强方法及系统
CN109919859A (zh) 一种户外场景图像去雾增强方法、计算设备及其存储介质
CN109978789A (zh) 一种基于Retinex算法与引导滤波的图像增强方法
CN115578284A (zh) 一种多场景图像增强方法及系统
CN110288546A (zh) 一种采用双向伽马变换的低照度图像增强方法
CN110580690B (zh) 一种识别峰值变换非线性曲线的图像增强方法
CN105426847A (zh) 低质量自然光虹膜图像非线性增强方法
CN114155173A (zh) 一种图像去雾方法、设备及非易失性存储介质
CN102789634A (zh) 一种获取光照均一化图像的方法
CN117611467A (zh) 一种能同时平衡不同区域细节和亮度的低光图像增强方法
CN117611501A (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
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20121024

Termination date: 20130921