CN110298796B - 基于改进Retinex与对数图像处理的低照度图像增强方法 - Google Patents

基于改进Retinex与对数图像处理的低照度图像增强方法 Download PDF

Info

Publication number
CN110298796B
CN110298796B CN201910431030.0A CN201910431030A CN110298796B CN 110298796 B CN110298796 B CN 110298796B CN 201910431030 A CN201910431030 A CN 201910431030A CN 110298796 B CN110298796 B CN 110298796B
Authority
CN
China
Prior art keywords
image
illumination component
pixel
image processing
logarithmic
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
CN201910431030.0A
Other languages
English (en)
Other versions
CN110298796A (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.)
Sun Yat Sen University
Original Assignee
Sun Yat Sen University
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 Sun Yat Sen University filed Critical Sun Yat Sen University
Priority to CN201910431030.0A priority Critical patent/CN110298796B/zh
Publication of CN110298796A publication Critical patent/CN110298796A/zh
Application granted granted Critical
Publication of CN110298796B publication Critical patent/CN110298796B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • G06T5/94Dynamic range modification of images or parts thereof based on local image properties, e.g. for local contrast enhancement
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection

Landscapes

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

Abstract

本发明提供了一种基于改进Retinex与对数图像处理的低照度图像增强方法,包括:计算样本图像的亮通道值,将该值作为样本图像的光照分量;利用现有的对数图像处理模型下的背景强度对光照分量进行自适应局部调整;结合Sobel边缘检测方法,对局部调整后的光照分量进行滤波细化;根据细化后的光照分量,基于Retinex理论得到增强图像。本发明提供一种基于改进Retinex与对数图像处理的低照度图像增强方法,在现有的对数图像处理模型下结合了Retinex理论,对样本图像进行了加强,有效地解决了Retinex算法可能存在的光晕效应以及过增强的问题,该方法能适应更多样的光照环境。

Description

基于改进Retinex与对数图像处理的低照度图像增强方法
技术领域
本发明涉及计算机数字图像处理技术领域,更具体地,涉及一种基于改进Retinex与对数图像处理的低照度图像增强方法。
背景技术
目前,各种光学成像装置广泛地应用在各种领域,但是,在环境光照亮度较低或光学传感器灵敏度较低的情况下,成像装置所成图像存在亮度不够、对比度降低等退化问题,进而导致图像观感下降或影响数字图像处理***的分析与识别等问题的出现。因此有必要研究低照度下图像增强算法,使光学成像***可以适用于缺乏光照的场景。
Retinex是低照度图像增强的经典算法,该算法可以处理常见的低照度图像,但难以处理光照情况比较复杂的图像。对于同时存在光照充足区域与光照不足区域的图像,使用Retinex算法处理后可能会出现过度增强、出现光晕效应等问题。
为了解决上述问题,一些改进方法为:采用多尺度的Retinex算法、基于暗通道的低照度图像增强、结合直方图均衡的增强算法,但仍不能满足复杂光照的场景,同时会出现颜色失真的问题,仍然未能满足需求。
发明内容
本发明为了克服现有改进的Retinex算法在处理图像的过程中存在无法满足复杂光照的场景及会出现颜色失真的技术缺陷,提供一种基于改进Retinex与对数图像处理的低照度图像增强方法。
为解决上述技术问题,本发明的技术方案如下:
基于改进Retinex与对数图像处理的低照度图像增强方法,其特征在于,包括以下步骤:
S1:计算样本图像的亮通道值,将该值作为样本图像的光照分量;
S2:利用现有的对数图像处理模型下的背景强度对光照分量进行自适应局部调整;
S3:结合Sobel边缘检测方法,对局部调整后的光照分量进行滤波细化;
S4:根据细化后的光照分量,基于Retinex理论得到增强图像。
其中,所述步骤S1具体包括以下步骤:
S11:基于亮通道先验理论,将样本图像进行颜色通道分离,将三个颜色通道图像逐像素点比较,将最大像素点提取并保存为最大值图像;
S12:对获得的最大值图像进行最大值滤波,从而获得样本图像的光照分量。
其中,所述步骤S2具体包括以下步骤:
S21:将样本图像的光照分量转换为对数图像处理模型下的表达,即将样本图像的灰度值表示为:
f(x,y)=M-I(x,y);
上式中,I(x,y)表示样本图像的像素值,f(x,y)为对数图像处理模型下的图像灰度值,M为图像像素值的最大值,对于8位数字图像为256;
S22:根据每个像素点的邻域像素值计算样本图像的平均背景强度;
S23:根据每个像素点的平均背景强度得到每个像素点对应的亮度调整系数;
S24:将亮度调整系数代入对数图像处理模型中的数乘运算,完成对光照分量的局部调整。
其中,步骤S22中所述的每个像素点的邻域像素值按8邻域进行检测而得到。
其中,根据人眼视觉的方向特性,人眼对水平和垂直方向上的亮度变化比斜线方向上的亮度变化更敏感,因此所述8邻域中四个邻域像素值与对角线上四个邻域像素值根据不同的比例权重计算背景强度,从而得到更符合人眼视觉特性的背景强度。
其中,所述步骤S3具体包括以下步骤:
S31:运用Sobel边缘检测方法提取局部调整后的光照分量的边缘图像;
S32:根据边缘图像求取各像素点的平滑因子;
S33:根据各像素点的平滑因子对局部调整后的光照分量进行滤波、细化,得到细化后的光照分量。
其中,步骤S31所述的Sobel边缘检测方法使用3×3的Sobel算子对局部调整后的光照分量进行边缘图像提取,分解为水平方向与垂直方向进行处理。
其中,所述步骤S4具体包括以下步骤:
S41:输入低照度的样本图像并进行颜色通道分离,得到三个通道各像素点的像素值;
S42:分别把三个通道各像素点的像素值除以经调整与细化后的光照分量中对应的像素点的像素值,得到三个新的颜色通道;
S43:将三个新的颜色通道进行合成,得到增强图像。
与现有技术相比,本发明技术方案的有益效果是:
本发明提供一种基于改进Retinex与对数图像处理的低照度图像增强方法,在现有的对数图像处理模型下结合了Retinex理论,对样本图像进行了加强,有效地解决了Retinex算法可能存在的光晕效应以及过增强的问题,该方法能适应更多样的光照环境。
附图说明
图1本发明方法的流程示意图;
图2是本发明中基于亮通道先验理论提取的光照分量;
图3是本发明中自适应调整与细化后的光照分量;
图4是本发明中的输入图像与增强后图像对比。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;
对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
如图1所示,基于改进Retinex与对数图像处理的低照度图像增强方法,其特征在于,包括以下步骤:
S1:计算样本图像的亮通道值,将该值作为样本图像的光照分量;
S2:利用现有的对数图像处理模型下的背景强度对光照分量进行自适应局部调整;
S3:结合Sobel边缘检测方法,对局部调整后的光照分量进行滤波细化;
S4:根据细化后的光照分量,基于Retinex理论得到增强图像。
在具体实施过程中,本方法在现有的对数图像处理模型下结合了Retinex理论,对样本图像进行了加强,有效地解决了Retinex算法可能存在的光晕效应以及过增强的问题,该方法能适应更多样的光照环境。
实施例2
更具体的,在实施例1的基础上,将任意图像采集设备采集到的图像作为样本图像,基于亮通道先验理论,将样本图像进行颜色通道分离,将三个颜色通道图像逐像素点比较,将最大像素点提取并保存为最大值图像;再对获得的最大值图像进行最大值滤波,从而获得样本图像的光照分量,结果如图2所示。
在具体实施过程中,基于对数图像处理模型,将样本图像的光照分量转换为对数图像处理模型下的表达,即将样本图像的灰度值表示为:
f(x,y)=M-I(x,y);
上式中,I(x,y)表示样本图像的像素值,f(x,y)为对数图像处理模型下的图像灰度值,M为图像像素值的最大值,对于8位数字图像为256;根据每个像素点的邻域像素值计算样本图像的平均背景强度,按8邻域进行检测;设四邻域比例系数为a,对角线上四邻域比例系数为b,平均背景强度根据以下公式计算:
Figure BDA0002069005890000041
上式中,a和b表示比例权重,∑QX(i,j)为像素点(i,j)的四个邻域像素的集合,∑DX(i,j)为像素点(i,j)的对角线上四个邻域像素的集合;其中,⊕与分别为对数图像处理模型中的相加与相乘操作;
根据平均背景强度B(x,y)计算变换系数γ(x,y),具体为:
γ(x,y)=a(1-cos(B(x,y)×π));
将变换系数γ(x,y)代入对数图像处理模型的数乘运算,得到以下公式:
Figure BDA0002069005890000042
式中,L(x,y)为增强后的光照分量,L0(x,y)为步骤S1所得到的样本图像的光照分量,至此得到局部调整后的光照分量。
在具体实施过程中,运用Sobel边缘检测方法提取局部调整后的光照分量的边缘图像,将得到的边缘图像s(x,y)代入下式中,求得各像素点的平滑因子ε(x,y):
Figure BDA0002069005890000043
根据各像素点的平滑因子对局部调整后的光照分量进行滤波、细化,得到细化后的光照分量,如图3所示。
在具体实施过程中,输入低照度的样本图像并进行颜色通道分离,得到三个通道各像素点的像素值,分别把三个通道各像素点的像素值除以经调整与细化后的光照分量中对应的像素点的像素值,得到三个新的颜色通道,最后将三个新的颜色通道进行合成,得到增强图像。增强后图像与原图像对比如图4所示。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。

Claims (5)

1.基于改进Retinex与对数图像处理的低照度图像增强方法,其特征在于,包括以下步骤:
S1:计算样本图像的亮通道值,将该值作为样本图像的光照分量;
S2:利用现有的对数图像处理模型下的背景强度对光照分量进行自适应局部调整;
步骤S2具体包括以下步骤:
S21:将样本图像的光照分量转换为对数图像处理模型下的表达,即将样本图像的灰度值表示为:
f(x,y)=M-I(x,y);
上式中,I(x,y)表示样本图像的像素值,f(x,y)为对数图像处理模型下的图像灰度值,M为图像像素值的最大值,对于8位数字图像为256;
S22:根据每个像素点的邻域像素值计算样本图像的平均背景强度;
S23:根据每个像素点的平均背景强度得到每个像素点对应的亮度调整系数;
S24:将亮度调整系数代入对数图像处理模型中的数乘运算,完成对光照分量的局部调整;
步骤S22中所述的每个像素点的邻域像素值按8邻域进行检测而得到;
根据人眼视觉的方向特性,人眼对水平和垂直方向上的亮度变化比斜线方向上的亮度变化更敏感,因此所述8邻域中四个邻域像素值与对角线上四个邻域像素值根据不同的比例权重计算背景强度,从而得到更符合人眼视觉特性的背景强度;
设四邻域比例系数为a,对角线上四邻域比例系数为b,平均背景强度根据以下公式计算:
Figure FDA0004165116420000011
上式中,a和b表示比例权重,∑QX(i,j)为像素点(i,j)的四个邻域像素的集合,∑DX(i,j)为像素点(i,j)的对角线上四个邻域像素的集合;其中,
Figure FDA0004165116420000013
与分别为对数图像处理模型中的相加与相乘操作;
根据平均背景强度B(x,y)计算变换系数γ(x,y),具体为:
γ(x,y)=a(1-cos(B(x,y)×π));
将变换系数γ(x,y)代入对数图像处理模型的数乘运算,得到以下公式:
Figure FDA0004165116420000012
式中,L(x,y)为增强后的光照分量,L0(x,y)为步骤S1所得到的样本图像的光照分量,至此得到局部调整后的光照分量;
S3:结合Sobel边缘检测方法,对局部调整后的光照分量进行滤波细化;
步骤S3具体包括以下步骤:
S31:运用Sobel边缘检测方法提取局部调整后的光照分量的边缘图像;
S32:根据边缘图像求取各像素点的平滑因子;
S33:根据各像素点的平滑因子对局部调整后的光照分量进行滤波、细化,得到细化后的光照分量;
S4:根据细化后的光照分量,基于Retinex理论得到增强图像;
步骤S4具体包括以下步骤:
S41:输入低照度的样本图像并进行颜色通道分离,得到三个通道各像素点的像素值;
S42:分别把三个通道各像素点的像素值除以经调整与细化后的光照分量中对应的像素点的像素值,得到三个新的颜色通道;
S43:将三个新的颜色通道进行合成,得到增强图像。
2.根据权利要求1所述的基于改进Retinex与对数图像处理的低照度图像增强方法,其特征在于:所述步骤S1具体包括以下步骤:
S11:基于亮通道先验理论,将样本图像进行颜色通道分离,将三个颜色通道图像逐像素点比较,将最大像素点提取并保存为最大值图像;
S12:对获得的最大值图像进行最大值滤波,从而获得样本图像的光照分量。
3.根据权利要求2所述的基于改进Retinex与对数图像处理的低照度图像增强方法,其特征在于:所述步骤S3具体包括以下步骤:
S31:运用Sobel边缘检测方法提取局部调整后的光照分量的边缘图像;
S32:根据边缘图像求取各像素点的平滑因子;
S33:根据各像素点的平滑因子对局部调整后的光照分量进行滤波、细化,得到细化后的光照分量。
4.根据权利要求3所述的基于改进Retinex与对数图像处理的低照度图像增强方法,其特征在于:步骤S31所述的Sobel边缘检测方法使用3×3的Sobel算子对局部调整后的光照分量进行边缘图像提取,分解为水平方向与垂直方向进行处理。
5.根据权利要求4所述的基于改进Retinex与对数图像处理的低照度图像增强方法,其特征在于,所述步骤S4具体包括以下步骤:
S41:输入低照度的样本图像并进行颜色通道分离,得到三个通道各像素点的像素值;
S42:分别把三个通道各像素点的像素值除以经调整与细化后的光照分量中对应的像素点的像素值,得到三个新的颜色通道;
S43:将三个新的颜色通道进行合成,得到增强图像。
CN201910431030.0A 2019-05-22 2019-05-22 基于改进Retinex与对数图像处理的低照度图像增强方法 Active CN110298796B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910431030.0A CN110298796B (zh) 2019-05-22 2019-05-22 基于改进Retinex与对数图像处理的低照度图像增强方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910431030.0A CN110298796B (zh) 2019-05-22 2019-05-22 基于改进Retinex与对数图像处理的低照度图像增强方法

Publications (2)

Publication Number Publication Date
CN110298796A CN110298796A (zh) 2019-10-01
CN110298796B true CN110298796B (zh) 2023-05-16

Family

ID=68027038

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910431030.0A Active CN110298796B (zh) 2019-05-22 2019-05-22 基于改进Retinex与对数图像处理的低照度图像增强方法

Country Status (1)

Country Link
CN (1) CN110298796B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112249830A (zh) * 2020-11-20 2021-01-22 中北大学 基于数据融合的电梯安全评估***
CN113096033B (zh) * 2021-03-22 2024-05-28 北京工业大学 基于Retinex模型自适应结构的低光照图像增强方法
CN113129236B (zh) * 2021-04-25 2022-07-12 中国石油大学(华东) 基于Retinex和卷积神经网络的单张低光照图像增强方法及***
CN114429426B (zh) * 2021-12-20 2022-08-16 哈尔滨理工大学 一种基于Retinex模型的低照度图像质量改善方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB201219920D0 (en) * 2012-11-05 2012-12-19 British Broadcasting Corp Method and apparatus for detection of chromatic aberration
CN104732555A (zh) * 2015-04-13 2015-06-24 南通理工学院 基于Sobel算子的图像边缘检测方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003061105A (ja) * 2001-06-07 2003-02-28 Seiko Epson Corp 画像処理方法、画像処理プログラム、画像処理装置およびそれを用いたディジタルスチルカメラ
CN103578084B (zh) * 2013-12-09 2017-04-19 西安电子科技大学 基于亮通道滤波的彩色图像增强方法
RU2622095C1 (ru) * 2016-02-05 2017-06-09 Акционерное общество "Государственный Рязанский приборный завод" Способ улучшения цифровых цветных изображений
CN106780375A (zh) * 2016-12-02 2017-05-31 南京邮电大学 一种低照度环境下的图像增强方法
CN106846282B (zh) * 2017-03-28 2019-06-04 华侨大学 一种采用自适应校正的低照度图像增强方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB201219920D0 (en) * 2012-11-05 2012-12-19 British Broadcasting Corp Method and apparatus for detection of chromatic aberration
CN104732555A (zh) * 2015-04-13 2015-06-24 南通理工学院 基于Sobel算子的图像边缘检测方法

Also Published As

Publication number Publication date
CN110298796A (zh) 2019-10-01

Similar Documents

Publication Publication Date Title
CN110298796B (zh) 基于改进Retinex与对数图像处理的低照度图像增强方法
US11127122B2 (en) Image enhancement method and system
CN109064426B (zh) 一种用于抑制低照度图像中眩光并增强图像的方法及装置
CN107680056B (zh) 一种图像处理方法及装置
WO2013168618A1 (ja) 画像処理装置及び画像処理方法
Pei et al. Effective image haze removal using dark channel prior and post-processing
CN110473152B (zh) 基于改进Retinex算法的图像增强方法
CN109214996B (zh) 一种图像处理方法及装置
CN110852956A (zh) 一种高动态范围图像的增强方法
CN104966266A (zh) 自动模糊身体部位的方法及***
CN111340717A (zh) 一种用于非制冷红外热像仪的图像预处理装置
CN114240789A (zh) 一种基于优化的保持亮度的红外图像直方图均衡增强方法
CN112991197B (zh) 一种基于暗通道的细节保留的低照度视频增强方法及装置
CN109377464A (zh) 一种红外图像的双平台直方图均衡化方法及其应用***
CN111223110A (zh) 一种显微图像增强方法、装置及计算机设备
Liu et al. Single image haze removal via depth-based contrast stretching transform
CN112435184A (zh) 一种基于Retinex和四元数的雾霾天图像识别方法
CN109859138B (zh) 一种基于人眼视觉特性的红外图像增强方法
CN117496019B (zh) 一种驱动静态图像的图像动画处理方法及***
CN113592750A (zh) 一种基于梯度直方图的红外增强方法
Han et al. Automatic illumination and color compensation using mean shift and sigma filter
Fathy et al. Benchmarking of pre-processing methods employed in facial image analysis
CN110136085B (zh) 一种图像的降噪方法及装置
CN116452447A (zh) 一种低照度高清图像处理方法
CN112532938B (zh) 一种基于大数据技术的视频监控***

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant