CN101510306B - 一种视频图像照度分布估计方法 - Google Patents
一种视频图像照度分布估计方法 Download PDFInfo
- Publication number
- CN101510306B CN101510306B CN2008101478621A CN200810147862A CN101510306B CN 101510306 B CN101510306 B CN 101510306B CN 2008101478621 A CN2008101478621 A CN 2008101478621A CN 200810147862 A CN200810147862 A CN 200810147862A CN 101510306 B CN101510306 B CN 101510306B
- Authority
- CN
- China
- Prior art keywords
- image
- gradient
- illumination
- imgl
- gradient vector
- 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
- 238000005286 illumination Methods 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims description 29
- 238000012216 screening Methods 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000010354 integration Effects 0.000 claims description 8
- 238000003708 edge detection Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 abstract description 6
- 238000013507 mapping Methods 0.000 abstract description 6
- 238000002310 reflectometry Methods 0.000 abstract description 5
- 238000012545 processing Methods 0.000 abstract description 4
- 230000002708 enhancing effect Effects 0.000 abstract description 3
- 238000011156 evaluation Methods 0.000 abstract 2
- 230000000694 effects Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000002950 deficient Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 108091008695 photoreceptors Proteins 0.000 description 2
- 230000001154 acute effect Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种视频图像照度分布估计方法,经过梯度筛选,去掉梯度幅度大于梯度阀值或边缘存在部分的梯度向量影响,经对筛选后的梯度向量图像进行照度图像重建,得到照度图像,可以有效去除场景表面反射率剧烈变化对图像照度分布统计带来的影响,从而准确地估计出视频图像采集场景中光源照度的特征和分布变化。本发明得到的视频图像照度分布估计可以为后续图像增强处理,如直方图均衡、色域映射、gamma校正等结合,进而达到视频图像质量的增强效果。
Description
技术领域
本发明涉及视频增强技术领域,具体来讲,涉及一种视频图像照度分布估计方法。
背景技术
在视频图像增强处理系统中,对图像的亮度分布进行均衡、映射、校正等处理是用以增强图像质量很重要的一类技术手段。如直方图均衡、色域映射、gamma校正等常用技术。然而对图像进行增强处理,往往需要对当前待处理的视频图像亮度分布进行估计。
现有视频图像亮度分布估计技术主要有两大类:一类是局部亮度分布估计技术,另一类是全局亮度分布估计技术。
局部亮度分布估计的典型代表是自适应直方图均衡技术,为自适应图像各个局部的亮度分布特征,局部地统计图像亮度灰度分布,从而计算在该图像局部所使用的均衡参数,从而达到自适应的直方图均衡目标。
全局亮度分布估计的一种典型方法是如图1所示,将图像分为5个区域一次编号为0、1、2、3、4。对于彩色图像需要先应用以下常用亮度分量公式计算出亮度:
Luminance=0.27red+0.67green+0.06blue (1)
其中,red、green、blue分别代表彩色图像红绿蓝三色分量。
依据公式(1)计算出每个像素的亮度值后,按图1所示分区依次计算出每个分区的亮度算数平均值。一种可能图像亮度分布的计算结果如图2所示。依据图2所示的图像亮度均值分布结果,可以大致判断出图像左上角明亮,右下角较暗。进而可以依据这些数据均定随后使用图像增强技术,如直方图均衡、色域映射等处理时的具体参数。
依据现有技术的视频图像亮度分布估计进行图像亮度分布校正的技术方案有一个根本的缺陷。具体来说,现有技术方案亮度分布估计的基础是依据公式(1)进行的逐像素点亮度值计算,其成像过程如图3、4所示。在图3中,光源101发出的光束102入射到场景物体104的表面105后,经过场景物体104吸纳一部分入射光,漫反射一部分入射光,最后一束反射光103进入摄像机。在图4中,物体反射光103经过摄像机镜头后的光束201投射到摄像机感光器表面202上,感光信号经过模数转换203和必要的摄像机图像处理流程204后最终输出场景104的数字图像形式205。场景物体104最终的数字图像一般以每个像素R/G/B三个通道逐像素表示的方式输出和存储。
由图3、4可以看出,公式(1)计算出的亮度值只反应了物体表面反射光束103的亮度水平,却并不能直接反应光源101投射到物体104的入射光束102,即光照的亮度状态。事实上,就一般自然图像成像场景来说,场景物体表面反射率在空间上会剧烈跃式变化,而现实中的光源照度一般只会在空间发生缓慢连续变化。所以现有传统视频图像亮度分布估计技术的缺陷是只反应料场景物体表面反射率在空间上剧烈的不规则变化,而场景物体表面反射率在空间上的变化是需要保留记录的场景物体有用信息,而由于光源照度变化导致的图像亮度不均衡才是需要处理的部分
发明内容
本发明的目的在于克服现有视频图像亮度分布估计的不足,提供一种能更好地估计出视频图像采集场景的照度分布和变化特征,从而使随后的图像增强处理的参数调整更为准确的视频图像照度分布估计方法。
为达到上述发明目的,本发明的视频图像照度分布估计方法,包括以下步骤:
(1)、计算出原始图像相应的亮度灰度图像;
(2)、使用基于图像梯度值的边缘检测方法检测并绘制出亮度灰度图像的边缘图像;
(3)、在亮度灰度图像的基础上,计算出梯度图像和梯度向量图像;
(4)、根据梯度图像,计算出梯度图像的梯度阀值;根据边缘图像以及梯度阀值对梯度向量图像进行梯度筛选,如果梯度向量图像中的某像素点相应梯度图像位置上的梯度幅度大于梯度阀值或边缘图像指示该位置存在于边缘,则把梯度向量图像中的某像素点的梯度向量设置为零,否则,保持;
(5)、对筛选后的梯度向量图像进行照度图像重建,得到照度图像;
(6)、根据照度图像进行照度分布估计,得到视频图像照度分布估计。
在本发明中,经过梯度筛选,去掉梯度幅度大于梯度阀值或边缘存在部分的梯度向量影响,经对筛选后的梯度向量图像进行照度图像重建,得到照度图像,可以有效去除场景表面反射率剧烈变化对图像照度分布统计带来的影响,从而准确地估计出视频图像采集场景中光源照度的特征和分布变化。本发明得到的视频图像照度分布估计可以为后续图像增强处理,如直方图均衡、色域映射、gamma校正等结合,进而达到视频图像质量的增强效果。
附图说明
图1是现有技术全局亮度分布估计示意图;
图2是图1所示分区计算出的亮度均值分布图;
图3是场景物体成像原理示意图;
图4是视频图像成像原理示意图;
图5是本发明一种具体实施方式流程图;
图6是图5所示的积分重构顺序图。
具体实施方式
为更好地理解本发明,下面结合具体实施方式对本发明进行更为详细描述。在以下的描述中,当已有的现有技术的详细描述也许会淡化本发明的主题内容时,这些描述在这儿将被忽略。
图1~4所示为现有技术,在背景技术中已经进行了描述和说明,在此不再赘述。
图5是本发明一种具体实施方式流程图。在本实施例中,包括以下步骤:
步骤ST1:原始图像Img每个像素根据以下常用亮度分量公式计算出亮度值:
Luminance=0.27red+0.67green+0.06blue (1)
其中,red、green、blue分别代表彩色图像红绿蓝三色分量。
得到原始图像相应的亮度灰度图像ImgL。
步骤ST2:使用基于图像梯度值的边缘检测方法检测并绘制出边缘图像ImgE。需要注意的是,本技术领域人员知道检测与绘制边缘图像并不限于这里提到的基于图像梯度值的方法,可以使用众多其它传统或新兴边缘检测方法。
步骤ST3:在亮度灰度图像ImgL的基础上逐像素计算其对数,得到对数图像ImgL’,计算公式如下:
ImgL′=log ImgL (2)
其中,ImgL′对数图像值,ImgL为亮度灰度图像值
步骤ST4:得到对数图像ImgL’后,在其基础上计算其梯度向量图像ImgGV和梯度图像ImgG,具体来说,使用空间差分计算,计算公式如下:
h(x,y)=ImgL′(x+1,y)-ImgL′(x,y);
V(x,y)=ImgL′(x,y+1)-ImgL′(x,y) (3)
其中(x,y)表示图像的空间坐标,ImgGV(x,y)=(h(x,y),V(x,y))构成该像素位置的梯度向量,ImgG(x,y)为该像素位置的梯度幅度值,从得到梯度向量图像ImgGV和梯度图像ImgG。
步骤ST5:根据梯度图像ImgG计算梯度阀值GradT,具体方法为:
GradMax=max(ImgG);
GradMin=min(imgG); (4)
GradT=GradMin+(GradMax-GradMin)*GradS
其中,GradS为一经验数据,在本实施中,取值为0.3。GradMax表示梯度图像ImgG中的最大梯度幅度值,GradMin表示梯度图像ImgG中的最小梯度幅度值
在边缘图像ImgE以及梯度阀值GradT的基础上,进行梯度筛选处理。具体来讲,就是逐像素判别处理,每个像素位置上如果梯度图像ImgG的梯度幅度ImgG(x,y)大于梯度阀值GradT或边缘图像ImgE指示该像素位置该位置存在于边缘,则把该位置梯度向量设置为零。具体计算公式如下:
步骤ST6:在筛选后的梯度向量图像ImgGV’的基础上,进行积分重构原始图像的对数照度图像ImgIL。图6是图5所示的积分重构顺序图,图中,P11表示图像左上角像素,逐行一次是P12,P13等,逐列依次是P21,P31等。在进行重构是P11位置的初始值使用从对数图像ImgL’对应位置P11的值,具体照度对数图像ImgIL的积分重构计算过程为:
ImgIL(1,1)=ImgL′(1,1);
ImgIL(x+1,1)=ImgIL(x,1)-ImgGV′(x,1)v; (6)
ImgIL(x,y+1)=ImgIL(x,y)-ImgGV′(x,y)h;
其中,ImgIL(x,y)表示照度对数图像在位置(x,y)上的值,ImgL′(1,1)表示对数图像在位置(x,y)上的值,ImgGV′(x,y)v表示梯度向量图像在位置(x,y)的垂直分量,ImgGV′(x,y)h表示梯度向量图像在位置(x,y)的水平分量,具体的积分重构计算线为图6所示的之字形逐行扫描。
步骤ST7:在得到照度对数图像ImgIL后,需要进行公式(2)的逆运算计算原始图像最终的照度图像ImgI,具体计算方法为:
ImgI=log′ImgIL (7)
其中,ImgI照度图像值,ImgIL为照度对数图像值。
步骤ST8:得到照度图像ImgI后,进行照度分布估计,得到视频图像照度分布估计。图像照度分布估计可以采用现有多种统计方法进行。
本实施例中,通过计算原始图像的边缘图像、对数梯度图像,结合梯度阀值的计算有效地去除对数梯度图像中梯度幅度值高于阀值或存在边缘的像素位置的梯度变化,随后经过积分重构和反对数计算得到原始图像的照度图像,然后进行照度分布估计,可以为后续图像增强处理,如直方图均衡、色域映射、gamma校正等结合,进而达到视频图像质量的增强效果。
尽管上面对本发明说明性的具体实施方式进行了描述,但应当清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (3)
1.一种视频图像照度分布估计方法,其特征在于,包括以下步骤:
(1)、计算出原始图像相应的亮度灰度图像;
(2)、使用基于图像梯度值的边缘检测方法检测并绘制出亮度灰度图像的边缘图像;
(3)、在亮度灰度图像的基础上,计算出梯度图像和梯度向量图像;
(4)、根据梯度图像,计算出梯度图像的梯度阀值;根据边缘图像以及梯度阀值对梯度向量图像进行梯度筛选,如果梯度向量图像中的某像素点相应梯度图像位置上的梯度幅度大于梯度阀值或边缘图像指示该位置存在于边缘,则把梯度向量图像中的某像素点的梯度向量设置为零,否则,保持;
(5)、对筛选后的梯度向量图像进行照度图像重建,得到照度图像;
(6)、根据照度图像进行照度分布估计,得到视频图像照度分布估计。
2.根据权利要求1所述的视频图像照度分布估计方法,其特征在于,步骤(3)所述的计算出梯度图像和梯度向量图像步骤为:
(31)、在亮度灰度图像的基础上逐像素计算其对数,得到对数图像,计算公式如下:
ImgL′=log ImgL
其中,ImgL′为对数图像值,ImgL为亮度灰度图像值
(32)、得到对数图像后,使用空间差分计算出其梯度向量图像和梯度图像,计算公式如下:
h(x,y)=ImgL′(x+1,y)-ImgL′(x,y);
V(x,y)=ImgL′(x,y+1)-ImgL′(x,y)
其中(x,y)表示图像的空间坐标,ImgGV(x,y)=(h(x,y),V(x,y))构成该像素位置的梯度向量,ImgG(x,y)为该像素位置的梯度幅度值,从得到梯度向量图像和梯度图像。
3.根据权利要求2所述的视频图像照度分布估计方法,其特征在于,步骤 (5)所述的对筛选后的梯度向量图像进行照度图像重建,步骤为:
(51)、在筛选后的梯度向量图像的基础上,进行积分重构原始图像的照度对数图像,照度对数图像的积分重构计算过程为:
ImgIL(1,1)=ImgL′(1,1);
ImgIL(x+1,1)=ImgIL(x,1)-ImgGV′(x,1)v;
ImgIL(x,y+1)=ImgIL(x,y)-ImgGV′(x,y)h;
其中,ImgIL(x,y)表示照度对数图像在位置(x,y)上的值,ImgL′(1,1)表示对数图像在位置(x,y)上的值,ImgGV′(x,y)v表示梯度向量图像在位置(x,y)的垂直分量,ImgGV′(x,y)h表示梯度向量图像在位置(x,y)的水平分量;
(52)、在得到照度对数图像后,需要进行逆运算计算原始图像最终的照度图像,具体计算方法为:
ImgI=log′ImgIL
其中,ImgI照度图像值,ImgIL为照度对数图像值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101478621A CN101510306B (zh) | 2008-12-15 | 2008-12-15 | 一种视频图像照度分布估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101478621A CN101510306B (zh) | 2008-12-15 | 2008-12-15 | 一种视频图像照度分布估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101510306A CN101510306A (zh) | 2009-08-19 |
CN101510306B true CN101510306B (zh) | 2010-12-22 |
Family
ID=41002698
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008101478621A Expired - Fee Related CN101510306B (zh) | 2008-12-15 | 2008-12-15 | 一种视频图像照度分布估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101510306B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102456222B (zh) * | 2010-10-29 | 2013-12-04 | 深圳迈瑞生物医疗电子股份有限公司 | 图像中组织均衡的方法和装置 |
CN103593847A (zh) * | 2013-11-25 | 2014-02-19 | 中国航天科工集团第三研究院第八三五七研究所 | 基于机器视觉的雨滴检测分析方法 |
CN108805220A (zh) * | 2018-07-25 | 2018-11-13 | 广东奥普特科技股份有限公司 | 一种基于梯度积分的快速模板匹配算法 |
CN109118541B (zh) * | 2018-08-24 | 2022-04-19 | 浙江大丰实业股份有限公司 | 舞台定向车台导轨位置分析系统 |
CN114565535B (zh) * | 2022-03-02 | 2022-09-13 | 北京锐影医疗技术有限公司 | 一种基于自适应梯度gamma校正的图像增强方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1139669A1 (en) * | 2000-03-28 | 2001-10-04 | STMicroelectronics S.r.l. | Coprocessor for motion estimation in digitised video sequence encoders |
CN101001377A (zh) * | 2007-01-09 | 2007-07-18 | 浙江大学 | 一种基于多参考帧的h.264快速运动估计的方法 |
CN101193303A (zh) * | 2007-09-14 | 2008-06-04 | 四川虹微技术有限公司 | 快速运动估计方法 |
-
2008
- 2008-12-15 CN CN2008101478621A patent/CN101510306B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1139669A1 (en) * | 2000-03-28 | 2001-10-04 | STMicroelectronics S.r.l. | Coprocessor for motion estimation in digitised video sequence encoders |
CN101001377A (zh) * | 2007-01-09 | 2007-07-18 | 浙江大学 | 一种基于多参考帧的h.264快速运动估计的方法 |
CN101193303A (zh) * | 2007-09-14 | 2008-06-04 | 四川虹微技术有限公司 | 快速运动估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101510306A (zh) | 2009-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9852499B2 (en) | Automatic selection of optimum algorithms for high dynamic range image processing based on scene classification | |
CN109754377B (zh) | 一种多曝光图像融合方法 | |
Park et al. | Single image dehazing with image entropy and information fidelity | |
CN106030653B (zh) | 用于生成高动态范围图像的图像处理系统和图像处理方法 | |
CN104067611B (zh) | 图像处理设备、图像处理方法及图像处理装置 | |
CN106991663B (zh) | 一种基于暗通道理论的水下彩色图像增强方法 | |
CN112785534B (zh) | 一种动态场景下去鬼影多曝光图像融合方法 | |
CN111161170A (zh) | 一种用于目标识别的水下图像综合增强方法 | |
CN108765342A (zh) | 一种基于改进暗通道的水下图像复原方法 | |
TWI518635B (zh) | 影像修復方法及使用該方法的影像處理裝置 | |
CN110473152B (zh) | 基于改进Retinex算法的图像增强方法 | |
CN107292830A (zh) | 低照度图像增强及评价方法 | |
CN112640413B (zh) | 用于显示周围环境的模型的方法、控制设备和车辆 | |
CN101908210A (zh) | 彩色图像去雾处理方法和系统 | |
US9165345B2 (en) | Method and system for noise reduction in video systems | |
CN105959510B (zh) | 一种视频快速去雾方法 | |
US20190068891A1 (en) | Method and apparatus for rapid improvement of smog/low-light-level image using mapping table | |
US9881358B2 (en) | Method and system for adaptive pixel replacement | |
CN101510306B (zh) | 一种视频图像照度分布估计方法 | |
CN107194900A (zh) | 图像处理方法、装置、计算机可读存储介质和移动终端 | |
CN104766307A (zh) | 一种图像处理的方法及设备 | |
TW201308251A (zh) | 水下影像強化系統 | |
US20250030952A1 (en) | Intensity separated local white balance correction | |
CN110023957B (zh) | 用于估计图像中的投射阴影区域和/或加亮区域的方法和设备 | |
CN110675351A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20101222 Termination date: 20151215 |
|
EXPY | Termination of patent right or utility model |