CN101739672B - 一种基于分区域插值的直方图均衡方法和装置 - Google Patents
一种基于分区域插值的直方图均衡方法和装置 Download PDFInfo
- Publication number
- CN101739672B CN101739672B CN200910242018.1A CN200910242018A CN101739672B CN 101739672 B CN101739672 B CN 101739672B CN 200910242018 A CN200910242018 A CN 200910242018A CN 101739672 B CN101739672 B CN 101739672B
- Authority
- CN
- China
- Prior art keywords
- histogram
- sub
- image
- num
- lim
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000012545 processing Methods 0.000 claims abstract description 23
- 238000013507 mapping Methods 0.000 claims description 40
- 239000011159 matrix material Substances 0.000 claims description 18
- 230000001932 seasonal effect Effects 0.000 claims description 8
- 238000011946 reduction process Methods 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 abstract description 3
- 230000000694 effects Effects 0.000 description 6
- 230000001276 controlling effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000007687 exposure technique Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Landscapes
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
- Picture Signal Circuits (AREA)
- Studio Devices (AREA)
Abstract
本发明涉及一种基于分区域插值的直方图均衡方法和装置,用于采集影像并数字化,以及对图像的预处理,将图像输入模块输入的数据划分成多个区域,然后对每个区域进行直方图统计,对直方图数据以及输入图像进行有限对比度调整,改变其直方图分布,重新分配溢出像素到各个直方图区间,并将有限对比度直方图调整的结果进行直方图均衡,将图像区域划分为若干子块,并对位于不同区域的相邻子块进行插值处理。本发明的思想是基于有限对比度直方图均衡技术对输入源图像的直方图进行分析,通过自适应的直方图规定化技术以及区域插值等过程,使得源图像中较暗部分变得变亮,凸现出更多的细节,同时又能够保证较亮部分不会因此过曝,丢失图像细节。
Description
技术领域
本发明属于计算机视频图像处理领域,涉及一种基于分区域插值的直方图均衡方法和装置。
背景技术
在明暗变化较大的场景中,例如倒车入库或出库时,由于车库外面有正常的日光照射,而车库里面由于顶棚的缘故,其光线明显不如外面强,对于人眼或者摄像头来说,处于充足日光照射的车库外面的细节,以及处于阴暗处的车库内的细节变得难以区分和辨认,对于通过车载图像获取装置所获取的图像而言,图像的动态范围太大而使得一些细节难以被清楚地呈现,而这些细节对于泊车者来说是有重要参考价值的,因此还是希望能看到外面的亮部细节,以及室内的暗部细节。又比如山洞内外的明暗差别较大,在山洞外面拍摄得到山洞里面的场景通常较暗且细节不明显。发生这种背光时,通常采用调整曝光时间的方式,来提高整帧图像的平均亮度,但是,在使某些区域变得可视的情况下不可避免的会使得图像中的另一些区域过曝或者细节消失。因此,有必要寻找一种局部自适应曝光技术,提高整帧图像的可视化效果。
发明内容
本发明针对现有技术中的不足提出一种局部曝光调整方案,其基于有限对比度直方图均衡技术对输入源图像的直方图进行分析,通过自适应的直方图规定化技术以及区域插值等过程,使得源图像中较暗部分变得变亮,凸现出更多的细节,同时又能够保证较亮部分不会因此过曝,丢失图像细节。
具体而言,本发明提供一种基于分区域插值的直方图均衡装置,包括以下模块:
图像输入模块,用于采集影像并数字化,以及对图像的预处理;
分区域直方图统计模块,用于将图像输入模块输入的数据划分成多个区域,然后对每个区域进行直方图统计;
有限对比度直方图调整模块,用于对直方图数据以及输入图像进行有限对比度调整,改变其直方图分布,重新分配溢出像素到各个直方图区间;
直方图均衡模块,用于将有限对比度直方图调整的结果进行直方图均衡;
区域映射插值模块,用于将图像区域划分为若干子块,并对位于不同区域的相邻子块进行插值处理。
以及一种基于分区域插值的直方图均衡方法,其特征在于包括以下步骤:
(1)将采集到的视频图像划分成多个区域,然后对每个区域进行直方图统计;
(2)对统计得到的直方图数据以及输入图像进行有限对比度调整,改变其直方图分布,重新分配溢出像素到各个直方图区间;
(3)对有限对比度直方图调整的结果进行直方图均衡;
(4)将图像区域划分为若干子块,并对位于不同区域的相邻子块进行插值处理。
本发明的效果在于:本发明基于有限对比度直方图均衡技术对输入源图像的直方图进行分析,通过自适应的直方图规定化技术以及区域插值等过程,使得源图像中较暗部分变得变亮,凸现出更多的细节,同时又能够保证较亮部分不会因此过曝,丢失图像细节。
附图说明
图1是基于分区域插值的直方图均衡装置的框图;
图2是分区域插值示意图;
图3是对角拼接处的相邻左上方子块UL、右上方子块UR、左下方子块BL、右下方子块BR进行插值处理的示意图;
图4是用于说明本发明效果的原图;
图5是本发明所述的方法对图4进行处理后的图像。
具体实施方式
下面结合说明书附图和具体实施方式对本发明作进一步详细描述。
本发明首先提供一种基于分区域插值的直方图均衡装置,如图1所示,包括图像输入模块、分区域直方图统计模块、有限对比直方图调整模块、区域映射差值模块以及输出图像模块。
其中,所述的图像输入模块用于采集影像并数字化,以及对图像的预处理,本发明的一个实施方式中,其包含一个影像获取单元和一个预处理单元,影像获取单元可以是CCD摄像机,或者包含CMOS感光元件的影像采集设备;在另一个实施方式中,其可以是一个小型摄像头,可固定在诸如车辆的移动设备之上,可以固定在车辆的背部、尾部。影像获取单元获取的视频图像经过预处理单元进行预处理,预处理包括去噪降噪处理,因为动态采集的视频图像存在模糊,噪声干扰等问题,为保证画面质量以及后续操作的效果,可以对影像获取单元获取的图像预先做降噪处理。
图像输入模块采集得到的数字图像被输入到分区域直方图统计模块,该模块所做的处理包括:
首先将输入的图像等分成M×N个区域,然后对每个区域进行直方图统计,在一个优选的实施方式中,M=N=8。
分区域直方图统计过程如下:将灰度区间分为NBin个区间,遍历整个图像,根据每个像素灰度或亮度的值(像素值)决定每个区间内的取值Hist[i],i=1,...,NBin,即如果有像素值落在该区间内,区间计数加1。更进一步,如果视频中连续两帧的差分值不大,即相邻两帧图像中像素值变化不大,则当前帧的直方图统计可以供下一帧处理使用,这种情况下,可以仅统计相邻两帧中像素值不同的位置的像素值,以快速地获得下一帧的直方图。
所获得的直方图数据以及输入图像发送到有限对比度直方图调整模块,在该模块中,首先对统计后的直方图进行有限对比度调整,改变其直方图分布。
首先,统计直方图中每个像素值(即像素值分布的区间数,共Nbin个)上超过对比度门限Clip_Lim的像素个数总和Num_TotalExcess;对比度门限Clip_Lim用于控制图像的对比度,该值越大则调整后的图像对比度越大,在一个实施方式中,该值可以通过以下方式来获得:
其中,Num_AllPixels为图像所有像素总数,NBin为直方图区间个数;α为对比度控制因子;
其次,计算得到每个直方图中每个像素值上的平均溢出像素,即待分配像素个数Num_AvgBinIncr,以及修正上限Upper_Lim;
Num_TotalExcess的计算过程如下:
待分配像素个数Num_AvgBinIncr,以及修正上限Upper_Lim的计算过程如下:
Num_AvgBinIncr=floor(Num_TotalExcess/NBin);(3)
Upper_Lim=Clip_Lim-Num_AvgBinIncr;(4)
然后修正直方图,重新分配溢出像素到各个bin内,过程如下:
如果该区间像素数Hist[i]大于Clip_Lim,那么该区间内像素数Hist[i]被钳位为Clip_Lim;否则
如果该区间像素数Hist[i]大于Upper_Lim,将Num_TotalExcess减去(Clip_Lim-Hist[i]);同时令Hist[i]为Clip_Lim;否则,将Num_TotalExcess减去(Clip_Lim-Num_AvgBinIncr);同时令Hist[i]为Hist[i]+Num_AvgBinIncr
接着分配剩余的像素,令步长step为max(floor(NBin/Num_TotalExcess),1),从1开始至NBin,以step为步长,执行操作:如果该区间像素数Hist[i]小于Clip_Lim,则令该区间像素数目加1。
有限对比度直方图调整的结果输入到直方图均衡模块,在直方图均衡模块中,
首先将上述步骤中修正后的直方图进行累加,生成递增序列HistSum;在本发明的一个实施方式中,可以采用Sn-1+an的方法得到HistSum,其中Sn-1为前n-1项和,而an为第n项,例如对于序列{0,2,5,3,5,8,12,1,4,...},累加后得到的HistSum结果为{0,2,7,10,15,23,35,36,40,...}
其次对HistSum序列进行均匀分布处理,获得每个子区域的映射曲线:
其中,Num_AllPixels为图像内像素个数;Value_High为图像像素亮度值上限,本发明中优选的Value_High取值为256;Value_Low为图像像素亮度值下限,本发明中Value_Low优选的取值范围为0-32之间的数。
得到每个子区域的映射曲线后,如果不进行插值使邻域间像素相关的话,可能在区域分界处产生明显的边界,影响图像效果。因此,必须进行邻域间插值。该功能由图1中区域映射插值模块来实现,其用于分别对各区域与其四周相邻区域边缘处的像素灰度值进行插值处理。
发明所采用的插值过程思想是分别将各区域划分为若干子块,并对位于不同区域的相邻子块进行插值处理。首先,将原子图像区域块再次二等分,可以分别将如图2实线所示各区域内等分为像素行列数相等的2×2共4个如虚线划分的子块,并对每4个对角拼接的相邻区域中,位于如图中箭头所引出的对角拼接处的4个相邻子块进行插值处理,即对图2中每4个对角拼接处的相邻左上方子块UL、右上方子块UR、左下方子块BL、右下方子块BR进行插值处理,可以按照如下方式进行插值处理:
L′i,j=rowRevW×[colRevW×Mapping(Li,j)UL+colW×Mapping(Li,j)UR]
+rowW×[colRevW×Mapping(Li,j)BL+colW×Mapping(Li,j)BR]
其中,L′i,j表示所述每个子块中第i行第j列先后经直方图均衡化和插值处理后的像素灰度值,rowW、rowRevW分别为正反方向的行插值系数矩阵,colW、colRevW分别为正反方向的列插值系数矩阵,Mapping(Li,j)UL、Mapping(Li,j)UR、Mapping(Li,j)BL、Mapping(Li,j)BR分别为左上方字块、右上方子块、左下方子块、右下方子块中第i行第j列经直方图均衡化的像素灰度值;在一个可选的实施方式中,
上述各矩阵中的rows为每个子块的行数,cols为每个子块的列数。
利用上述正反方向的行插值系数矩阵rowW和rowRevW、以及正反方向的列插值系数矩阵colW、colRevW,即可使得各子块中越靠近拼接角顶端的像素插值系数越大、越远离拼接角顶端的像素插值系数越小。
且,由于每个子块中的像素在插值运算时,以需要使用、且仅使用一次自身所在行和列的其余像素,因而正反方向的行插值系数矩阵rowW和rowRevW中,只能有一个矩阵的插值系数选取1-rows、另一个则只能选取0-rows-1,正反方向的列插值系数矩阵colW、colRevW中,只能有一个矩阵的插值系数选取1-cols、另一个则只能选取cols-1,从而能够避免在插值处理时重复使用每个像素所在行和列的其余像素灰度值。
需要说明的是,对于边缘区域中未标示出UL、UR、BL、BR的子块,可以不进行处理、或采用任意适用于边缘的处理方式进行处理。上述插值操作的更为形象但非限定性的描述可以由图3来表示。
图像输出模块用于输出处理后的图像。
与上述基于分区域插值的直方图均衡装置相应地,本发明同时还提出一种基于分区域插值的直方图均衡方法,其特征在于包括以下步骤:
将采集到的视频图像划分成多个区域,然后对每个区域进行直方图统计;
对统计得到的直方图数据以及输入图像进行有限对比度调整,改变其直方图分布,重新分配溢出像素到各个直方图区间;
对有限对比度直方图调整的结果进行直方图均衡;
将图像区域划分为若干子块,并对位于不同区域的相邻子块进行插值处理。
所述对每个区域进行直方图统计的步骤进一步包括:将输入的图像等分成M×N个区域,然后对每个区域进行直方图统计,包括将灰度区间分为NBin个区间,遍历整个图像,根据每个像素灰度或亮度的值决定每个区间内的取值Hist[i],i=1,...,NBin,如果有像素值落在该区间内,在该区间计数加1。
所述对统计得到的直方图数据以及输入图像进行有限对比度调整,改变其直方图分布,重新分配溢出像素到各个直方图区间的步骤进一步包括以下步骤:
首先,统计直方图中每个区间上像素值大于对比度门限Clip_Lim的像素个数总和Num_TotalExcess;对比度门限Clip_Lim通过以下方式来获得:
其中,Num_AllPixels为图像所有像素总数,NBin为直方图区间个数;α为对比度控制因子;
其次,计算得到每个直方图中每个像素值上的平均溢出像素,即待分配像素个数Num_AvgBinIncr,以及修正上限Upper_Lim;Num_TotalExcess的计算过程如下:
待分配像素个数Num_AvgBinIncr,以及修正上限Upper_Lim的计算过程如下:
Num_AvgBinIncr=floor(Num_TotalExcess/NBin);
Upper_Lim=Clip_Lim-Num_AvgBinIncr;
改变图像直方图分布,重新分配溢出像素到各个直方图,其步骤为:
如果该区间像素数Hist[i]大于Cklip_Lim,那么该区间内像素数Hist[i]被钳位为Clip_Lim;否则
如果该区间像素数Hist[i]大于Upper_Lim,将Num_TotalExcess减去(Clip_Lim-Hist[i]);同时令Hist[i]为Clip_Lim;否则,将Num_TotalExcess减去(Clip_Lim-Num_AvgBinIncr);同时令Hist[i]为Hist[i]+Num_AvgBinIncr;
接着分配剩余的像素,令步长step为max(floor(NBin/Num_TotalExcess),1),从1开始至NBin,以step为步长,执行操作:如果该区间像素数Hist[i]小于Clip_Lim,则令该区间像素数目加1。
所述直方图均衡模块中执行以下操作:对有限对比度直方图调整模块修正后的直方图进行累加,生成递增序列HistSum;
对HistSum进行均匀分布处理,获得每个子区域的映射曲线:
其中,Num_AllPixels为图像内像素个数;Value_High为图像像素亮度值上限。
得到每个子区域的映射曲线后,可以按照如下方式进行插值处理:
L′i,j=rowRevW×[colRevW×Mapping(Li,j)UL+colW×Mapping(Li,j)UR]
+rowW×[colRevW×Mapping(Li,j)BL+colW×Mapping(Li,j)BR]
其中,L′i,j表示所述每个子块中第i行第j列先后经直方图均衡化和插值处理后的像素灰度值,rowW、rowRevW分别为正反方向的行插值系数矩阵,colW、colRevW分别为正反方向的列插值系数矩阵,Mapping(Li,j)UL、Mapping(Li,j)UR、Mapping(Li,j)BL、Mapping(Li,j)BR分别为左上方字块、右上方子块、左下方子块、右下方子块中第i行第j列经直方图均衡化的像素灰度值。rowW、rowRevW、colW、colRevW矩阵可以采用与基于分区域插值的直方图均衡装置相同的rowW、rowRevW、colW、colRevW矩阵来实现。
图4和图5可以用来说明本发明提出的基于分区域插值的直方图均衡装置及方法的效果,图4为原图,图片中间门洞内细节部分无法看清楚;图5为分区域直方图均衡技术处理后的图片,门洞内的细节部分清晰可见。本发明提出的基于分区域直方图均衡化装置和方法,将图像分为若干子区域计算其直方图分布,相比于全局均衡化技术,引入了空间约束,算法保证了图像整体的对比度;并且所述基于分区域直方图均衡化装置的模块结构清晰,各部分分工明确,独立性强,便于提高整体的稳定性和可靠性,并且引入对比度限定参数,可以精确控制输出图像的整体对比度,算法的人机操控性强,所使用的参数少,计算方法简便,硬件代价低,可以嵌入硬件芯片中,实用性强。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (7)
1.一种基于分区域插值的直方图均衡装置,包括以下模块:
图像输入模块,用于采集影像并数字化,以及对图像的预处理;
分区域直方图统计模块,用于将图像输入模块输入的数据划分成多个区域,然后对每个区域进行直方图统计;
有限对比度直方图调整模块,用于对直方图数据以及输入图像进行有限对比度调整,改变其直方图分布,重新分配溢出像素到各个直方图区间;所述有限对比度直方图调整模块对直方图数据以及输入图像进行有限对比度调整,改变其直方图分布进一步包括以下操作:
首先,统计直方图中每个区间上像素值大于对比度门限Clip_Lim的像素个数总和Num_TotalExcess;对比度门限Clip_Lim通过以下方式来获得:
其中,Num_AllPixels为图像所有像素总数,NBin为直方图区间个数;α为对比度控制因子;
其次,计算得到每个直方图中每个像素值上的平均溢出像素,即待分配像素个数Num_AvgBinIncr,以及修正上限Upper_Lim;Num_TotalExcess的计算过程如下:
其中,Hist[k]表示第k个直方图区间的像素数;
待分配像素个数Num_AvgBinIncr,以及修正上限Upper_Lim的计算过程如下:
Num_AvgBinIncr=floor(Num_TotalExcess/NBin);Upper_Lim=Clip_Lim-Num_AvgBinIncr;
直方图均衡模块,用于将有限对比度直方图调整的结果进行直方图均衡;
区域映射插值模块,用于将图像区域划分为若干子块,并对位于不同区域的相邻子块进行插值处理;所述区域映射插值模块执行的插值操作包括:
首先,将图像区域块等分为像素行列数相等的2×2子块,并对每4个对角拼接的相邻区域中对角拼接处的4个相邻子块,即左上方子块UL、右上方子块UR、左下方子块BL、右下方子块BR进行插值处理:
L'i,j=rowRevW×[colRevW×Mapping(Li,j)UL+colW×Mapping(Li,j)UR]
+rowW×[colRevW×Mapping(Li,j)BL+colW×Mapping(Li,j)BR]
其中,L'i,j表示所述每个子块中第i行第j列先后经直方图均衡化和插值处理后的像素灰度值,rowW、rowRevW分别为正反方向的行插值系数矩阵,colW、colRevW分别为正反方向的列插值系数矩阵,Mapping(Li,j)UL、Mapping(Li,j)UR、Mapping(Li,j)BL、Mapping(Li,j)BR分别为左上方子块、右上方子块、左下方子块、右下方子块中第i行第j列经直方图均衡化的像素灰度值;
其中,
上述各矩阵中的Rows为每个子块的行数,cols为每个子块的列数。
2.如权利要求1所述的一种基于分区域插值的直方图均衡装置,其特征在于:
所述图像输入模块包含一个影像获取单元和一个预处理单元,所述影像获取单元用于获取视频图像,所获取的视频图像经过预处理单元进行预处理,预处理包括去噪降噪处理;
所述分区域直方图统计模块对每个区域进行直方图统计进一步包括:将输入的图像等分成M×N个区域,然后对每个区域进行直方图统计,包括将灰度区间分为NBin个区间,遍历整个图像,根据每个像素灰度或亮度的值决定每个区间内的取值Hist[k],k=1,...,NBin,如果有像素值落在该区间内,在该区间计数加1。
3.如权利要求1所述的一种基于分区域插值的直方图均衡装置,其特征在于:所述有限对比度直方图调整模块对直方图数据以及输入图像进行有限对比度调整,改变其直方图分布,重新分配溢出像素到各个直方图的具体操作还包括:
如果该区间像素数Hist[k]大于Clip_Lim,那么该区间内像素数Hist[k]被钳位为Clip_Lim;
否则,如果该区间像素数Hist[k]大于Upper_Lim,将Num_TotalExcess减去(Clip_Lim-Hist[k]);同时令Hist[k]为Clip_Lim;
否则,将Num_TotalExcess减去(Clip_Lim-Num_AvgBinIncr);同时令Hist[k]为Hist[k]+Num_AvgBinIncr;
接着分配剩余的像素,令步长step为max(floor(NBin/Num_TotalExcess),1),从1开始至NBin,以step为步长,执行操作:如果该区间像素数Hist[k]小于Clip_Lim,则令该区间像素数目加1。
4.如权利要求1所述的一种基于分区域插值的直方图均衡装置,其特征在于,所述直方图均衡模块中执行以下操作:
对有限对比度直方图调整模块修正后的直方图进行累加,生成递增序列HistSum;
对HistSum进行均匀分布处理,获得每个子区域的映射曲线:
其中,Num_AllPixels为图像内像素个数;Value_High为图像像素亮度值上限,Value_Low为图像像素亮度值下限。
5.一种基于分区域插值的直方图均衡方法,其包括以下步骤:
(1)将采集到的视频图像划分成多个区域,然后对每个区域进行直方图统计;
(2)对统计得到的直方图数据以及输入图像进行有限对比度调整,改变其直方图分布,重新分配溢出像素到各个直方图区间;所述对统计得到的直方图数据以及输入图像进行有限对比度调整,改变其直方图分布,重新分配溢出像素到各个直方图区间的步骤进一步包括以下步骤:
首先,统计直方图中每个区间上像素值大于对比度门限Clip_Lim的像素个数总和Num_TotalExcess;对比度门限Clip_Lim通过以下方式来获得:
其中,Num_AllPixels为图像所有像素总数,NBin为直方图区间个数;α为对比度控制因子;
其次,计算得到每个直方图中每个像素值上的平均溢出像素,即待分配像素个数Num_AvgBinIncr,以及修正上限Upper_Lim;Num_TotalExcess的计算过程如下:
其中,Hist[k]表示第k个直方图区间的像素数;
待分配像素个数Num_AvgBinIncr,以及修正上限Upper_Lim的计算过程如下:
Num_AvgBinIncr=floor(Num_TotalExcess/NBin);
Upper_Lim=Clip_Lim-Num_AvgBinIncr;
改变图像直方图分布,重新分配溢出像素到各个直方图,其步骤为:
如果该区间像素数Hist[k]大于Clip_Lim,那么该区间内像素数Hist[k]被钳位为Clip_Lim;
否则,如果该区间像素数Hist[k]大于Upper_Lim,将Num_TotalExcess减去(Clip_Lim-Hist[k]);同时令Hist[k]为Clip_Lim;
否则,将Num_TotalExcess减去(Clip_Lim-Num_AvgBinIncr);同时令Hist[k]为Hist[k]+Num_AvgBinIncr;
接着分配剩余的像素,令步长step为max(floor(NBin/Num_TotalExcess),1),从1开始至NBin,以step为步长,执行操作:如果该区间像素数Hist[k]小于Clip_Lim,则令该区间像素数目加1;
(3)对有限对比度直方图调整的结果进行直方图均衡;
(4)将图像区域划分为若干子块,并对位于不同区域的相邻子块进行插值处理;插值操作包括:
首先,将图像区域块等分为像素行列数相等的2×2子块,并对每4个对角拼接的相邻区域中对角拼接处的4个相邻子块,即左上方子块UL、右上方子块UR、左下方子块BL、右下方子块BR进行插值处理:
L'i,j=rowRevW×[colRevW×Mapping(Li,j)UL+colW×Mapping(Li,j)UR]
+rowW×[colRevW×Mapping(Li,j)BL+colW×Mapping(Li,j)BR]
其中,L'i,j表示所述每个子块中第i行第j列先后经直方图均衡化和插值处理后的像素灰度值,rowW、rowRevW分别为正反方向的行插值系数矩阵,colW、colRevW分别为正反方向的列插值系数矩阵,Mapping(Li,j)UL、Mapping(Li,j)UR、 Mapping(Li,j)BL、Mapping(Li,j)BR分别为左上方子块、右上方子块、左下方子块、右下方子块中第i行第j列经直方图均衡化的像素灰度值。
6.如权利要求5所述的一种基于分区域插值的直方图均衡方法,其特征在于:所述对每个区域进行直方图统计的步骤进一步包括:将输入的图像等分成M×N个区域,然后对每个区域进行直方图统计,包括将灰度区间分为NBin个区间,遍历整个图像,根据每个像素灰度或亮度的值决定每个区间内的取值Hist[k],k=1,...,NBin,如果有像素值落在该区间内,在该区间计数加1。
7.如权利要求5所述的一种基于分区域插值的直方图均衡方法,其特征在于:步骤(3)中,对有限对比度直方图调整的结果进行直方图均衡包括:
对进行有限对比度修正后的直方图进行累加,生成递增序列HistSum;
对HistSum进行均匀分布处理,获得每个子区域的映射曲线:
其中,Num_AllPixels为图像内像素个数;Value_High为图像像素亮度值上限,Value_Low为图像像素亮度值下限。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910242018.1A CN101739672B (zh) | 2009-12-02 | 2009-12-02 | 一种基于分区域插值的直方图均衡方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910242018.1A CN101739672B (zh) | 2009-12-02 | 2009-12-02 | 一种基于分区域插值的直方图均衡方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101739672A CN101739672A (zh) | 2010-06-16 |
CN101739672B true CN101739672B (zh) | 2015-08-05 |
Family
ID=42463124
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910242018.1A Active CN101739672B (zh) | 2009-12-02 | 2009-12-02 | 一种基于分区域插值的直方图均衡方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101739672B (zh) |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102298770B (zh) * | 2011-08-11 | 2014-03-12 | 奇瑞汽车股份有限公司 | 一种图像对比度增强方法和装置 |
CN103440635B (zh) * | 2013-09-17 | 2016-06-22 | 厦门美图网科技有限公司 | 一种基于学习的对比度受限自适应直方图均衡方法 |
CN103957362B (zh) * | 2014-05-05 | 2017-06-06 | 广东欧珀移动通信有限公司 | 一种自动调整相机曝光效果的方法及*** |
CN104361566B (zh) * | 2014-11-17 | 2018-08-10 | 厦门美图之家科技有限公司 | 一种优化暗部区域的图像处理方法 |
CN106296594A (zh) * | 2015-06-02 | 2017-01-04 | 南京普爱医疗设备股份有限公司 | 一种基于直方图的x射线机影像***的图像伽玛处理算法 |
CN106530262A (zh) * | 2015-09-15 | 2017-03-22 | 株式会社理光 | 图像合成方法及设备 |
CN105243192B (zh) * | 2015-09-16 | 2018-09-21 | 中国商用飞机有限责任公司 | 质量加载方法和*** |
CN105376498A (zh) * | 2015-10-16 | 2016-03-02 | 凌云光技术集团有限责任公司 | 一种扩展相机动态范围的图像处理方法及*** |
CN105913405B (zh) * | 2016-04-05 | 2019-03-29 | 智车优行科技(北京)有限公司 | 用于呈现图像细节的处理方法、装置及车辆 |
CN106097306B (zh) * | 2016-05-31 | 2019-01-29 | 凌云光技术集团有限责任公司 | 获取图像边缘检测算子的方法、图像边缘检测方法及装置 |
CN106570833A (zh) * | 2016-10-31 | 2017-04-19 | 湖南国科微电子股份有限公司 | 视频图像自适应对比度增强方法和装置 |
CN108805873A (zh) * | 2017-04-28 | 2018-11-13 | 南宁富桂精密工业有限公司 | 图像处理方法及装置 |
TWI622019B (zh) * | 2017-08-31 | 2018-04-21 | 元智大學 | 基於保熵的映射先驗的影像增強方法及其影像處理裝置 |
CN109658363A (zh) * | 2018-10-22 | 2019-04-19 | 长江大学 | 子块自适应融合的多层子块重叠直方图均衡方法及*** |
CN112153240B (zh) | 2019-06-27 | 2021-11-09 | 深圳Tcl数字技术有限公司 | 一种调整图像画质方法、装置及可读存储介质 |
CN110751603A (zh) * | 2019-09-26 | 2020-02-04 | 中国科学院长春光学精密机械与物理研究所 | 一种图像对比度的增强方法、***及终端设备 |
CN111325680B (zh) * | 2020-01-08 | 2023-04-25 | 深圳深知未来智能有限公司 | 一种具有抑制局部过曝的图像增亮方法 |
CN111261088B (zh) * | 2020-02-25 | 2023-12-12 | 京东方科技集团股份有限公司 | 一种图像绘制方法、装置及显示装置 |
CN113160098B (zh) * | 2021-04-16 | 2022-11-11 | 浙江大学 | 一种光照不均匀条件下密集微粒图像的处理方法 |
CN113255684B (zh) * | 2021-06-28 | 2021-10-26 | 恒银金融科技股份有限公司 | 一种基于灰度图溢出的背景分离方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101510305A (zh) * | 2008-12-15 | 2009-08-19 | 四川虹微技术有限公司 | 一种改进的自适应直方图均衡方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1840831A1 (en) * | 2006-03-31 | 2007-10-03 | Sony Deutschland Gmbh | Adaptive histogram equalization for images with strong local contrast |
-
2009
- 2009-12-02 CN CN200910242018.1A patent/CN101739672B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101510305A (zh) * | 2008-12-15 | 2009-08-19 | 四川虹微技术有限公司 | 一种改进的自适应直方图均衡方法 |
Non-Patent Citations (3)
Title |
---|
Adaptive histogram equalization and its variations;Stephen M. Pizer等;《Computer Vision, Graphics, and Image Processing》;19870930;第39卷(第3期);第355–368页 * |
一种基于改进直方图均衡的图像增强方法;高岩等;《计算机仿真》;20080430;第25卷(第4期);第198-201页 * |
基于对比度受限自适应直方图均衡的乳腺图像增强;刘轩等;《计算机工程与应用》;20080522;第44卷(第10期);第173-175页 * |
Also Published As
Publication number | Publication date |
---|---|
CN101739672A (zh) | 2010-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101739672B (zh) | 一种基于分区域插值的直方图均衡方法和装置 | |
CN101930534B (zh) | 用于人脸检测的动态图像压缩方法 | |
EP3076364B1 (en) | Image filtering based on image gradients | |
CN102077572B (zh) | 用于在成像***中防止运动模糊和重影的方法及装置 | |
CN101719989B (zh) | 一种背光补偿的方法及*** | |
JP5530563B2 (ja) | 高ダイナミックレンジ画像の捕捉及びレンダリング | |
WO2021164234A1 (zh) | 图像处理方法以及图像处理装置 | |
US9041834B2 (en) | Systems and methods for reducing noise in video streams | |
KR101426298B1 (ko) | 안개 제거율을 높인 영상보정장치 및 방법 | |
US20110019094A1 (en) | System and method for random noise estimation in a sequence of images | |
CN102542552B (zh) | 视频图像的顺逆光判断和拍摄时间检测方法 | |
CN111064904A (zh) | 一种暗光图像增强方法 | |
CN105578068A (zh) | 一种高动态范围图像的生成方法、装置及移动终端 | |
WO2019056549A1 (zh) | 图像增强方法以及图像处理装置 | |
DE102020200310A1 (de) | Verfahren und System zur Dunstreduzierung für die Bildverarbeitung | |
US20220198625A1 (en) | High-dynamic-range image generation with pre-combination denoising | |
EP3534325B1 (en) | Image processing device and method for compressing dynamic range of wide-range image | |
US11790501B2 (en) | Training method for video stabilization and image processing device using the same | |
CN103313068A (zh) | 基于灰边缘约束灰世界的白平衡校正图像处理方法及装置 | |
CN112819699A (zh) | 视频处理方法、装置及电子设备 | |
CN112991236B (zh) | 一种基于模板的图像增强方法及装置 | |
CN116389912B (zh) | 脉冲相机融合普通相机重构高帧率高动态范围视频的方法 | |
CN117611467A (zh) | 一种能同时平衡不同区域细节和亮度的低光图像增强方法 | |
CN117274107A (zh) | 低照度场景下端到端色彩及细节增强方法、装置及设备 | |
US20230394632A1 (en) | Method and image processing device for improving signal-to-noise ratio of image frame sequences |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20201210 Address after: No.602, 6 / F, shining building, 35 Xueyuan Road, Haidian District, Beijing 100191 Patentee after: BEIJING ZHONGXINGTIANSHI TECHNOLOGY Co.,Ltd. Address before: 100083, Haidian District, Xueyuan Road, Beijing No. 35, Nanjing Ning building, 15 Floor Patentee before: Vimicro Corp. |
|
TR01 | Transfer of patent right |