CN102722872B - 数字图像局部高光溢出消除方法 - Google Patents

数字图像局部高光溢出消除方法 Download PDF

Info

Publication number
CN102722872B
CN102722872B CN201210172911.3A CN201210172911A CN102722872B CN 102722872 B CN102722872 B CN 102722872B CN 201210172911 A CN201210172911 A CN 201210172911A CN 102722872 B CN102722872 B CN 102722872B
Authority
CN
China
Prior art keywords
light
roi
value
brightness
idea
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.)
Active
Application number
CN201210172911.3A
Other languages
English (en)
Other versions
CN102722872A (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.)
Synthesis Electronic Technology Co Ltd
Original Assignee
SHANDONG SYNTHESIS ELECTRONIC 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 SHANDONG SYNTHESIS ELECTRONIC TECHNOLOGY Co Ltd filed Critical SHANDONG SYNTHESIS ELECTRONIC TECHNOLOGY Co Ltd
Priority to CN201210172911.3A priority Critical patent/CN102722872B/zh
Publication of CN102722872A publication Critical patent/CN102722872A/zh
Application granted granted Critical
Publication of CN102722872B publication Critical patent/CN102722872B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Processing (AREA)

Abstract

本发明提供一种数字图像局部高光溢出消除方法,可以有效处理常见的数字摄影图像在被摄体表面光滑无纹理区域的高光溢出现象,包括处理被摄体表面曲率较大或点光源这两类特殊情况。本发明采用以下技术方案:一种数字图像局部高光溢出消除方法,包括以下步骤:1)根据数字图像亮度曲线中未溢出部分的数据估算溢出部分的理论亮度数值;2)对于待调整区间内亮度做非线性变换,压缩最大亮度值规定范围,得到较好的修复效果。

Description

数字图像局部高光溢出消除方法
技术领域
本发明涉及一种数字图像处理方法,具体涉及一种数字图像局部高光溢出消除方法。 
背景技术
由于被摄物体的反光特性,通过摄影手段获取的数字图像有时出现局部高光溢出现象。局部高光溢出是指数字图像因被摄体某些部位反射光线过强,照成摄影手段获得的数字图像局部区域像素的采样值超出允许范围。 
例如普通的可见光彩色摄影所获得的JPEG格式的图片,通常每个像素由R、G、B三个分量构成,每个分量为0~255范围的整数。当被摄物体局部反光过强,就会出现R、G、B三个分量中的一个或多个取值超出255的情况,在这种情况下,数字图像只能将其亮度值设置为255。一旦出现这类亮度溢出现象,就会造成照片局部亮度变化失真或者颜色失真。参考图1是数码照片中一条水平线上的像素亮度值的分布曲线,展示了高光溢出的后果。 
当高光溢出区域位于被摄体无纹理的表面区域内时,通常用下述方法消除高光溢出:1、单色填充法,利用局部高光区域附近像素的数值填充高光区域,使得高光溢出区域看上与附近未出现高光溢出的区域颜色接近。该方法的缺点是被消除高光溢出的区域视觉效果不自然,参见图2。2、扫描线线性平滑过渡法 
对高光溢出区域进行水平或竖直扫描,扫描线段两端的像素颜色分量亮度值取其相邻的非高光溢出像素的对应颜色分量的亮度值,扫描线其它部分的像素颜色分量值利用线性插值计算得到。该方法视觉效果有所改善。由于线性插值无法表现表面呈曲面的被摄体或者点光源的反光情况,因此该方法适用于高光溢出区域范围很小,或者被摄体表面曲率较小而且光源为漫射光的情形,参见图3。
发明内容
本发明的目的是提供一种数字图像局部高光溢出消除方法,可以有效处理常见的数字摄影图像在被摄体表面光滑无纹理区域的高光溢出现象,包括处理被摄体表面曲率较大或点光源这两类特殊情况。 
为了实现上述目的,本发明采用以下技术方案: 
一种数字图像局部高光溢出消除方法,包括以下步骤:
1)根据数字图像亮度曲线中未溢出部分的数据估算溢出部分的理论亮度数值;
2)对于待调整区间内亮度做非线性变换,压缩最大亮度值规定范围,得到较好的修复效果。
进一步,所述数字图像局部高光溢出消除方法,包括以下步骤: 
1)在待处理数字图像Image亮度曲线中确定一个内部包含待处理的高光溢出区域的且长宽分别与两坐标轴平行的矩形区域作为感兴趣区域;
2)确定待处理数字图像Image亮度曲线中在水平扫描线段上各通道亮度的预测值,并保存到感兴趣区域水平扫描用的临时浮点数组中;
3)确定待处理数字图像Image在竖直扫描线段上各通道亮度的预测值,并保存到感兴趣区域竖直扫描用的临时浮点数组中;
4)压缩所保存在临时数组中的数字图像Image亮度曲线中在水平扫描线段上各通道亮度的预测值和竖直扫描线段上各通道亮度的预测值,使其范围至不超过输出待处理数字图像Image通道的高光溢出阈值;
5)根据压缩后的数字图像Image亮度曲线中在水平扫描线段上各通道亮度的预测值和竖直扫描线段上各通道亮度的预测值,计算出亮度值混合函数;
6)通过亮度值混合函数确定感兴趣区域内每个像素在各通道的亮度值,亮度值混合函数使得感兴趣区域边界上的像素颜色值与原始图像相同;
7)根据上述6)中确定好的亮度值调整数字图像Image亮度曲线中溢出部分的亮度值,以消除高光溢出处理好图像。
进一步,所述数字图像局部高光溢出消除方法中,水平扫描线段上各通道亮度的预测值的步骤进一步包括: 
1)定义待处理数字图像通道的亮度函数Light(x) ,使亮度函数Light(x)等于数字图像通道在亮度曲线中坐标(x,y)处的像素通道的像素值Pixel(x, y, c);
2)确定待处理数字图像Image的理想亮度函数Idea(x),
Idea(ROI_X1) = Light(ROI_X1)
Idea(ROI_X2) = Light(ROI_X2)
Idea’(ROI_X1) = Light’(ROI_X1)
Idea’(ROI_X2) = Light’(ROI_X2)
用三次多项式构造待处理数字图像通道Image在水平扫描线段(x1,y)-(x2,y)通道c上的理想亮度函数Idea(x),具体步骤如下:
设 
Figure 141645DEST_PATH_IMAGE001
,根据Idea(x)构造要求可列出下面方程组:
Figure 621692DEST_PATH_IMAGE002
                           (1)
                          (2)
Figure 238061DEST_PATH_IMAGE004
                (3)
Figure 986574DEST_PATH_IMAGE005
               (4)
用高斯消元法解上面方程组,求得a,b,c,d的值,完成Idea(x)的构造;
3)对于感兴趣区域中每一个值,TmpH(x, y, c) =Idea(x),TmpH(x, y, c)是对待处理数字图像Image中像素(x,y)在通道c的水平扫描方向的亮度预测值。
进一步,所述数字图像局部高光溢出消除方法中,竖直扫描线段上各通道亮度的预测值的步骤进一步包括: 
1)定义通道竖直的亮度函数Light(y) = Pixel(x, y, c);
2)构造通道理想亮度函数Idea(y),满足条件:
Idea(ROI_Y1) = Light(ROI_Y1)
Idea(ROI_Y2) = Light(ROI_Y2)
Idea’(ROI_Y1) = Light’(ROI_Y1)
Idea’(ROI_Y2) = Light’(ROI_Y2)
按照如下方法构造理想亮度函数:
用三次多项式构造Image在水平扫描线段(x,y1)-(x,y2)通道c上的理想亮度函数Idea(y),具体步骤如下:
Figure 318198DEST_PATH_IMAGE006
,根据Idea(y)构造要求可列出下面方程组:
Figure 618598DEST_PATH_IMAGE007
                           (1)
Figure 481512DEST_PATH_IMAGE008
                          (2)
Figure 895700DEST_PATH_IMAGE009
                (3)
Figure 160460DEST_PATH_IMAGE010
               (4)
用高斯消元法解上面方程组,求得a,b,c,d的值,完成Idea(y)的构造;
3)对于感兴趣区域中每一个值,TmpV(x, y, c) = Idea(y),TmpV(x, y, c) 
是对待处理数字图像Image中像素(x,y)在通道c的竖直扫描方向的亮度预测值
进一步,所述数字图像局部高光溢出消除方法中,分别压缩所保存的图像各通道按水平扫描和竖直扫描所得的亮度预测值,使其取值范围至不超过输出待处理数字图像Image通道的高光溢出阈值Level(c),进一步包括以下步骤:
1)把通道c亮度预测值中的最大值存入变量MaxLight;
2)把待处理数字图像Image在感兴趣区域即ROI区域边界上的各像素通道c的亮度的最大值存入变量MinLight;
3)构造亮度转换函数ConvertLight满足下述条件:
ConvertLight(MinLight) = MinLight
ComvertLight(MaxLight) ≤ Level(c);
4)针对ROI内的每一个像素(x, y),执行如下操作:
v = Tmp(x, y, c),V为待处理数字图像Image在坐标(x,y)处像素通道C的像素值;
w = ConvertLight(v),
Tmp(x, y, c) = w,
其中w,v是临时变量,上述操作实现各像素(x,y)在通道c内按照水平或者竖直方向扫描的预测值压缩到规定范围。
进一步地,所述数字图像局部高光溢出消除方法中,按照下述方法构造亮度转换函数ConvertLight: 
设置K = 2×Level(c) – MaxLight;
判断K > MinLight成立,执行:
light < K时,ConvertLight(light) = light;
light≥K时,ConvertLight(light) = (light + K) / 2;
判断K > MinLight不成立,执行:
ConvertLight(light) = MinLight +(light - MinLight) × (Level(c) – MinLight) / (MaxLight – MinLight)。
进一步地,所述数字图像局部高光溢出消除方法中,构造亮度值混合函数MixValue(d1, d2, v1, v2)满足如下条件: 
d1 = min(x - ROI_X1, ROI_X2 – x) ,
d2 = min(y - ROI_Y1, ROI_Y2 – y),
当d1 = 0 时,MixValue (d1, d2, v1, v2) = v1;
当d1 > 0 且d2 = 0时,MixValue(d1, d2, v1, v2) = v2, 
当d1>0且D2>0时,
Min(v1, v2) ≤ MixValue (d1, d2, v1, v2) ≤ Max(v1, v2)
MixLight = MixValue (d1, d2, v1, v2)。
   进一步地,所述数字图像局部高光溢出消除方法中,按照下属方法构造亮度值混合函数MixValue(d1, d2, v1, v2): 
MixValue(d1, d2, v1, v2) = v1, (d1 = 0且 d2 = 0)
MixValue(d1, d2, v1, v2) = (d1*v2 + d2*v1)/(d1 + d2), (d1 > 0或d2 > 0)。
本发明的有益效果: 
可以有效处理常见的数字摄影图像在被摄体表面光滑无纹理区域的高光溢出现象,包括处理被摄体表面曲率较大或点光源这两类特殊情况,且处理图形的结果更加自然。
附图说明
图1为数码照片中一条水平线上的像素亮度值的分布曲线,展示了高光溢出的后果; 
图2单色填充法消除高光溢出后结果示意图。
图3为扫面线线性平滑过渡法消除高光溢出后结果示意图。 
图4为理论上预测溢出部分的理论数值。 
图5为根据本发明原理消除高光溢出后结果示意图。 
图6为K小于level大于minlight时,亮度转换函数图像构造原理示意图。 
图7为K小于minlight时,亮度转换函数图像构造原理示意图。 
具体实施方式
下面结合附图和具体实施例对本发明做进一步详细说明,下述说明仅是示例性的,不限定本发明的保护范围。 
一种数字图像局部高光溢出消除方法,包括以下步骤: 
1)参考附图4,根据数字图像亮度曲线中未溢出部分的数据估算溢出部分的理论亮度数值;
2)参考附图5,对于待调整区间内亮度做非线性变换,压缩最大亮度值规定范围,得到较好的修复效果。
本发明约定: 
Figure DEST_PATH_RE-2651DEST_PATH_IMAGE002
一种数字图像局部高光溢出消除方法,具体包括以下步骤:
1、在待处理数字图像Image亮度曲线中确定一个内部包含待处理的高光溢出区域的且长宽分别与两坐标轴平行的矩形区域作为感兴趣区域;
2、确定待处理数字图像Image亮度曲线中在水平扫描线段上各通道亮度的预测值,并保存到感兴趣区域水平扫描用的临时浮点数组中;
水平扫描线段上各通道亮度的预测值的步骤进一步包括:
1)定义待处理数字图像通道的亮度函数Light(x) ,使亮度函数Light(x)等于数字图像通道在亮度曲线中坐标(x,y)处的像素通道的像素值Pixel(x, y, c);
2)确定待处理数字图像通道Image的理想亮度函数Idea(x),
Idea(ROI_X1) = Light(ROI_X1)
Idea(ROI_X2) = Light(ROI_X2)
Idea’(ROI_X1) = Light’(ROI_X1)
Idea’(ROI_X2) = Light’(ROI_X2)
用三次多项式构造待处理数字图像通道Image在水平扫描线段(x1,y)-(x2,y)通道c上的理想亮度函数Idea(x),具体步骤如下:
Figure DEST_PATH_RE-112602DEST_PATH_IMAGE003
,根据Idea(x)构造要求可列出下面方程组:
Figure DEST_PATH_RE-830022DEST_PATH_IMAGE004
                           (1)
Figure DEST_PATH_RE-412182DEST_PATH_IMAGE005
                          (2)
Figure DEST_PATH_RE-164237DEST_PATH_IMAGE006
                (3)
Figure DEST_PATH_RE-173651DEST_PATH_IMAGE007
               (4)
用高斯消元法解上面方程组,求得a,b,c,d的值,完成Idea(x)的构造;
3)对于感兴趣区域中每一个值,TmpH(x, y, c) =Idea(x),TmpH(x, y, c)即为水平扫描段和竖直扫描段中临时浮点数;
3、确定待处理数字图像Image在竖直扫描线段上各通道亮度的预测值,并保存到感兴趣区域竖直扫描用的临时浮点数中;
竖直扫描线段上各通道亮度的预测值的步骤进一步包括:
1)定义通道竖直的亮度函数Light(y) = Pixel(x, y, c);
2)构造通道理想亮度函数Idea(y),满足条件:
Idea(ROI_Y1) = Light(ROI_Y1)
Idea(ROI_Y2) = Light(ROI_Y2)
Idea’(ROI_Y1) = Light’(ROI_Y1)
Idea’(ROI_Y2) = Light’(ROI_Y2)
按照如下方法构造理想亮度函数:
用三次多项式构造Image在水平扫描线段(x,y1)-(x,y2)通道c上的理想亮度函数Idea(y),具体步骤如下:
Figure DEST_PATH_RE-440684DEST_PATH_IMAGE008
,根据Idea(y)构造要求可列出下面方程组:
                           (1)
Figure DEST_PATH_RE-433097DEST_PATH_IMAGE010
                          (2)
Figure DEST_PATH_RE-364143DEST_PATH_IMAGE011
                (3)
Figure DEST_PATH_RE-305424DEST_PATH_IMAGE012
               (4)
用高斯消元法解上面方程组,求得a,b,c,d的值,完成Idea(y)的构造;
3)对于感兴趣区域中每一个值,TmpV(x, y, c) = Idea(y),TmpV(x, y, c) 为临时浮点数组;
4、压缩所保存各通道的亮度值取值范围至不超过输出待处理数字图像Image通道的高光溢出阈值;压缩所保存各通道的亮度值取值范围至不超过输出待处理数字图像Image通道的高光溢出阈值进一步包括以下步骤:
1)把临时浮点数组中通道c中的最大值存入变量MaxLight;
2)把待处理数字图像Image在感兴趣区域即ROI区域边界上的各像素通道c的亮度的最大值存入变量MinLight;
3)构造亮度转换函数ConvertLight满足下述条件:
ConvertLight(MinLight) = MinLight
ComvertLight(MaxLight) ≤ Level(c);
4)针对ROI内的每一个像素(x, y),执行如下操作:
v = Tmp(x, y, c),V为待处理数字图像Image在坐标(x,y)处像素通道C的像素值;
w = ConvertLight(v),
Tmp(x, y, c) = w;
按照下述方法构造亮度转换函数ConvertLight:
设置K = 2×Level(c) – MaxLight,即求K值使得Level(c)是K与MaxLight的平均值;
参见附图6,判断K > MinLight成立,执行:
light < K时, ConvertLight (light) = light,也就是说此时压缩后的亮度与原始图片一样,压缩函数图像对应附图OA段;
light≥K时, ConvertLight (light) = (light + K) / 2,此时函数图象对应AB’段,实现把最大亮度MaxLight压缩至Level(c)。虚线AB表示亮度不压缩时函数的变化趋势。
判断K > MinLight不成立,执行:
ConvertLight(light) = MinLight +(light - MinLight) × (Level(c) – MinLight) / (MaxLight – MinLight),对应的函数图象如下图7,虚线AB表示亮度不压缩时函数的变化趋势;
5、根据压缩后的数字图像Image亮度曲线中在水平扫描线段上各通道亮度的预测值和竖直扫描线段上各通道亮度的预测值,计算出亮度值混合函数;构造数据混合函数MixValue(d1, d2, v1, v2)满足如下条件:
d1 = min(x - ROI_X1, ROI_X2 – x) ,
d2 = min(y - ROI_Y1, ROI_Y2 – y),
当d1 = 0 时,MixValue (d1, d2, v1, v2) = v1;
当d1 > 0 且d2 = 0时,MixValue(d1, d2, v1, v2) = v2, 
Min(v1, v2) ≤ MixValue (d1, d2, v1, v2) ≤ Max(v1, v2)
MixLight = MixValue (d1, d2, v1, v2);
按照下属方法构造数据混合函数MixValue(d1, d2, v1, v2):
当d1 > 0 且d2 = 0时,MixValue(d1, d2, v1, v2) = v1, 
当d1>0且d2>0时,MixValue(d1, d2, v1, v2) = (d1*v2 + d2*v1)/(d1 + d2)。
6、如果需要感兴趣区域中某点像素离左右区域边缘近,就采用水平扫描段上亮度的预测值,如果某点像素离上下区域边缘近,就采取竖直扫描端上亮度的预测值;如果某点像素在区域中心部分,通过亮度值混合函数确定某点亮度值;
7、根据上述6中确定好的亮度值调整数字图像Image亮度曲线中溢出部分的亮度值,以消除高光溢出处理好图像。
实施例
身份证人像照片在拍照后被剪裁为大小为358×441,RGB彩色模式的JIEG照片。对身份证照片的处理环节之一是消除皮肤高光区域。本实施实例在标出感兴趣区域ROI后,采用水平扫描的方法消除ROI内高光溢出。具体步骤如下: 
一、身份证照片Image的颜色通道数Channels = 3,三个通道分别用R、G、B表示。设置R、G、B各通道阈值:Level(R) = 250;Level(G) = 240;Level(B) = 225。
二、构造通道理想亮度函数Idea(x), 
本实例用三次多项式构造Image在水平扫描线短(x1,y)-(x2,y)通道c上的理想亮度函数Idea(x),具体步骤如下:
Figure 833614DEST_PATH_IMAGE023
,根据Idea(x)构造要求可列出下面方程组:
                           (1)
 
Figure 364139DEST_PATH_IMAGE003
                         (2)
                (3)
Figure 220416DEST_PATH_IMAGE017
               (4)
用高斯消元法解上面方程组,求得a,b,c,d的值,完成Idea(x)的构造。因为本实例图像像素坐标为离散值,上述方程式(3)、(4)右端用差分代替导数计算,完成Idea(y)的构造;
三、按照下述方法构造亮度转换函数ConvertLight:
设置K = 2×Level(c) – MaxLight;
判断K > MinLight成立,执行:
light < K时,ConvertLight(light) = light;
light≥K时,ConvertLight(light) = (light + K) / 2;
判断K > MinLight不成立,执行:
ConvertLight(light) = MinLight +(light - MinLight) × (Level(c) – MinLight) / (MaxLight – MinLight)。
四、数据混合函数MixValue的构造方法 
d1 = min(x - ROI_X1, ROI_X2 – x) ,
d2 = min(y - ROI_Y1, ROI_Y2 – y),
MixValue(d1, d2, v1, v2) = v1, (d1 = 0且 d2 = 0)
MixValue(d1, d2, v1, v2) = (d1×v2 + d2×v1)/(d1 + d2), (d1 > 0或d2 > 0)。
五、根据上述中确定好的亮度值调整数字图像Image亮度曲线中溢出部分的亮度值,以消除高光溢出处理好图像。 
以上所述仅是本发明的优选实施例而已,并非对本发明做任何形式的限制,虽然本发明已经通过上述实施例披露了,但是不能用于限定本发明,本领域人员可利用上述揭示的技术内容作出许多变化或改进的等效实施例,这些等效实施例应属于本发明的技术方案的范围。 

Claims (7)

1.一种数字图像局部高光溢出消除方法,首先根据数字图像亮度曲线中未溢出部分的数据估算溢出部分的理论亮度数值;其次对于待调整区间内亮度做非线性变换,压缩最大亮度值规定范围,得到较好的修复效果;其特征在于包括以下步骤:
1)在待处理数字图像亮度曲线中确定一个内部包含待处理的高光溢出区域的且长宽分别与两坐标轴平行的矩形区域作为感兴趣区域;
2)确定待处理数字图像亮度曲线中在水平扫描线段上各通道亮度的预测值,并保存到感兴趣区域水平扫描用的临时浮点数组中;
3)确定待处理数字图像在竖直扫描线段上各通道亮度的预测值,并保存到感兴趣区域竖直扫描用的临时浮点数组中;
4)压缩所保存在临时数组中的数字图像亮度曲线中在水平扫描线段上各通道亮度的预测值和竖直扫描线段上各通道亮度的预测值,使其范围至不超过输出待处理数字图像通道的高光溢出阈值;
5)根据压缩后的数字图像亮度曲线中在水平扫描线段上各通道亮度的预测值和竖直扫描线段上各通道亮度的预测值,计算出亮度值混合函数;
6)通过亮度值混合函数确定感兴趣区域内每个像素在各通道的亮度值,亮度值混合函数使得感兴趣区域边界上的像素颜色值与原始图像相同;
7)根据上述6)中确定好的亮度值调整数字图像亮度曲线中溢出部分的亮度值,以消除高光溢出处理好图像。
2.根据权利要求1所述数字图像局部高光溢出消除方法,其特征在于:水平扫描线段上各通道亮度的预测值的步骤进一步包括:
1)定义待处理数字图像通道的亮度函数Light(x) ,使亮度函数Light(x)等于数字图像通道在亮度曲线中坐标(x,y)处的像素通道的像素值Pixel(x, y, c);
2)确定待处理数字图像的理想亮度函数Idea(x),
Idea(ROI_X1) = Light(ROI_X1)
Idea(ROI_X2) = Light(ROI_X2)
Idea’(ROI_X1) = Light’(ROI_X1)
Idea’(ROI_X2) = Light’(ROI_X2)
用三次多项式构造待处理数字图像通道Image在水平扫描线段(x1,y)-(x2,y)通道c上的理想亮度函数Idea(x),具体步骤如下:
设                                                
Figure 502300DEST_PATH_IMAGE001
,根据Idea(x)构造要求可列出下面方程组:
Figure 668708DEST_PATH_IMAGE002
                        (1)
Figure 323812DEST_PATH_IMAGE003
                         (2)
             (3)
                (4)
用高斯消元法解上面方程组,求得a,b,c,d的值,完成Idea(x)的构造;
3)对于感兴趣区域中每一个值,TmpH(x, y, c) =Idea(x),TmpH(x, y, c)是对待处理数字图像中像素(x,y)在通道c的水平扫描方向的亮度预测值。
3.根据权利要求1所述数字图像局部高光溢出消除方法,其特征在于:竖直扫描线段上各通道亮度的预测值的步骤进一步包括:
1)定义通道竖直的亮度函数Light(y) = Pixel(x, y, c);
2)构造通道理想亮度函数Idea(y),满足条件:
Idea(ROI_Y1) = Light(ROI_Y1)
Idea(ROI_Y2) = Light(ROI_Y2)
Idea’(ROI_Y1) = Light’(ROI_Y1)
Idea’(ROI_Y2) = Light’(ROI_Y2)
按照如下方法构造理想亮度函数:
用三次多项式构造Image在水平扫描线段(x,y1)-(x,y2)通道c上的理想亮度函数Idea(y),具体步骤如下:
Figure 792074DEST_PATH_IMAGE006
,根据Idea(y)构造要求可列出下面方程组:
                        (1)
Figure 450818DEST_PATH_IMAGE008
                        (2)
             (3)
Figure 50743DEST_PATH_IMAGE010
             (4)
用高斯消元法解上面方程组,求得a,b,c,d的值,完成Idea(y)的构造;
3)对于感兴趣区域中每一个值,TmpV(x, y, c) = Idea(y),TmpV(x, y, c) 
是对待处理数字图像中像素(x,y)在通道c的竖直扫描方向的亮度预测值。
4.根据权利要求1所述数字图像局部高光溢出消除方法,其特征在于,分别压缩所保存的图像各通道按水平扫描和竖直扫描所得的亮度预测值,使其取值范围至不超过输出待处理数字图像通道的高光溢出阈值Level(c),进一步包括以下步骤:
1)把通道c亮度预测值中的最大值存入变量MaxLight;
2)把待处理数字图像在感兴趣区域即ROI区域边界上的各像素通道c的亮度的最大值存入变量MinLight;
3)构造亮度转换函数ConvertLight满足下述条件:
ConvertLight(MinLight) = MinLight
ComvertLight(MaxLight) ≤ Level(c);
4)针对ROI内的每一个像素(x, y),执行如下操作:
v = Tmp(x, y, c),V为待处理数字图像在坐标(x,y)处像素通道C的像素值;
w = ConvertLight(v),
Tmp(x, y, c) = w,
其中w,v是临时变量,上述操作实现各像素(x,y)在通道c内按照水平或者竖直方向扫描的预测值压缩到规定范围。
5.根据权利要求4所述数字图像局部高光溢出消除方法,其特征在于:按照下述方法构造亮度转换函数ConvertLight:
设置K = 2×Level(c) – MaxLight;
判断K > MinLight成立,执行:
light < K时,ConvertLight(light) = light;
light≥K时,ConvertLight(light) = (light + K) / 2;
判断K > MinLight不成立,执行:
ConvertLight(light) = MinLight +(light - MinLight) × (Level(c) – MinLight) / (MaxLight – MinLight)。
6.根据权利要求1所述数字图像局部高光溢出消除方法,其特征在于:构造亮度值混合函数MixValue(d1, d2, v1, v2)满足如下条件:
d1 = min(x - ROI_X1, ROI_X2 – x) ,
d2 = min(y - ROI_Y1, ROI_Y2 – y),
当d1 = 0 时,MixValue (d1, d2, v1, v2) = v1;
当d1 > 0 且d2 = 0时,MixValue(d1, d2, v1, v2) = v2, 
当d1>0且d2>0时,
Min(v1, v2) ≤ MixValue (d1, d2, v1, v2) ≤ Max(v1, v2)
MixLight = MixValue (d1, d2, v1, v2)。
7.根据权利要求6所述数字图像局部高光溢出消除方法,其特征在于:按照下属方法构造亮度值混合函数MixValue(d1, d2, v1, v2):
当d1 > 0 且d2 = 0时,MixValue(d1, d2, v1, v2) = v1, 
当d1>0且d2>0时,MixValue(d1, d2, v1, v2) = (d1*v2 + d2*v1)/(d1 + d2)。
CN201210172911.3A 2012-05-30 2012-05-30 数字图像局部高光溢出消除方法 Active CN102722872B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210172911.3A CN102722872B (zh) 2012-05-30 2012-05-30 数字图像局部高光溢出消除方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210172911.3A CN102722872B (zh) 2012-05-30 2012-05-30 数字图像局部高光溢出消除方法

Publications (2)

Publication Number Publication Date
CN102722872A CN102722872A (zh) 2012-10-10
CN102722872B true CN102722872B (zh) 2014-07-09

Family

ID=46948619

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210172911.3A Active CN102722872B (zh) 2012-05-30 2012-05-30 数字图像局部高光溢出消除方法

Country Status (1)

Country Link
CN (1) CN102722872B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105023249B (zh) * 2015-06-26 2017-11-17 清华大学深圳研究生院 基于光场的高光图像修复方法及装置
CN106296617B (zh) * 2016-08-22 2019-03-05 腾讯科技(深圳)有限公司 人脸图像的处理方法及装置
CN107256376A (zh) * 2017-03-22 2017-10-17 武汉布偶猫科技有限公司 一种基于序贯主分量分析的人脸高光区处理方法
CN107392841B (zh) * 2017-06-16 2020-04-24 Oppo广东移动通信有限公司 人脸区域中黑眼圈消除方法、装置及终端
CN107392858B (zh) * 2017-06-16 2020-09-29 Oppo广东移动通信有限公司 图像高光区域处理方法、装置和终端设备
CN112312001B (zh) * 2019-07-30 2022-08-23 阿波罗智能技术(北京)有限公司 一种图像检测的方法、装置、设备和计算机存储介质
CN110443237B (zh) * 2019-08-06 2023-06-30 北京旷视科技有限公司 证件识别方法、装置、电子设备及计算机可读存储介质
CN110536131B (zh) * 2019-08-29 2021-01-05 中国科学院长春光学精密机械与物理研究所 Tdi探测器的高光溢出性能测试系统及测试方法
CN111882495B (zh) * 2020-07-05 2022-08-02 东北林业大学 一种基于自定义模糊逻辑与gan的图像高光处理方法
CN112837243B (zh) * 2021-03-05 2023-05-30 华侨大学 联合整体与局部信息的阴道镜图像高光消除的方法及装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4325703B2 (ja) * 2007-05-24 2009-09-02 ソニー株式会社 固体撮像装置、固体撮像装置の信号処理装置および信号処理方法、ならびに撮像装置
US8144985B2 (en) * 2007-12-21 2012-03-27 Sony Corporation Method of high dynamic range compression with detail preservation and noise constraints
JP5399739B2 (ja) * 2009-02-25 2014-01-29 ルネサスエレクトロニクス株式会社 画像処理装置
CN101783963B (zh) * 2010-02-10 2012-07-18 西安理工大学 一种具有高光抑制的夜间图像增强方法

Also Published As

Publication number Publication date
CN102722872A (zh) 2012-10-10

Similar Documents

Publication Publication Date Title
CN102722872B (zh) 数字图像局部高光溢出消除方法
EP3341913B1 (en) Inverse tone mapping based on luminance zones
CN111311482B (zh) 背景虚化方法、装置、终端设备及存储介质
US7751641B2 (en) Method and system for digital image enhancement
US9401027B2 (en) Method and apparatus for scene segmentation from focal stack images
KR101812341B1 (ko) 이미지의 에지 향상 방법
CN111127318A (zh) 一种机场环境下的全景图像拼接方法
CN104240211A (zh) 用于视频拼接的图像亮度与色彩均衡方法及系统
CN109767408B (zh) 图像处理方法、装置、存储介质及计算机设备
EP3203439B1 (en) Method and device for reducing noise in a component of a picture
CN103886565A (zh) 基于目的性优化及直方图均衡化的夜间彩色图像增强方法
CN107424137B (zh) 一种文本增强方法及装置、计算机装置、可读存储介质
CN105763768A (zh) 一种图像处理方法、装置和系统
CN111476735B (zh) 人脸图像处理方法、装置、计算机设备及可读存储介质
CN111311481A (zh) 背景虚化方法、装置、终端设备及存储介质
CN110175967B (zh) 图像去雾处理方法、系统、计算机设备和存储介质
SM et al. Underwater image enhancement using single scale retinex on a reconfigurable hardware
CN112767278B (zh) 基于非均匀大气光先验的图像去雾方法及相关设备
JP2004523839A (ja) 動的なチェーンベースの閾値化
EP2469834B1 (en) Image processing method, image processing device, and program
CN110648297B (zh) 图像去雾方法、系统、电子设备和存储介质
CN202183124U (zh) 一种基于fpga的图像景物边缘提取装置
CN105991937A (zh) 一种基于Bayer格式图像的虚拟曝光方法及装置
Vazquez-Corral et al. Angular-based preprocessing for image denoising
CN116263942A (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
C56 Change in the name or address of the patentee
CP03 Change of name, title or address

Address after: Shun high tech Zone of Ji'nan City, Shandong province 250101 China West Road No. 699

Patentee after: SYNTHESIS ELECTRONIC TECHNOLOGY CO., LTD.

Address before: 250101, No. 699, Shun Hua West Road, Lixia hi tech Zone, Lixia District, Shandong, Ji'nan

Patentee before: Shandong Synthesis Electronic Technology Co., Ltd.