一种控制平均亮度的直方图均衡方法技术

本发明专利技术提供了一种控制平均亮度的直方图均衡方法，按颜色模型抽取数字图像中的亮度数值，根据统计后的量度直方图计算亮度、平均亮度、均衡系数等一系列数值后，映射图像亮度并还原数字图像。本发明专利技术提供的控制平均亮度的直方图均衡方法，可以有效直接地控制直方图均衡带来的亮度变化，节省存储空间。

【技术实现步骤摘要】

本专利技术涉及，该方法应用于静止图像处理、视频图像增强等

技术介绍
直方图是对输出数据进行分组统计图，它反映了目标数据的分布几率特征。直方图均衡是一种利用该分组统计图进行分布特征优化的技术，优化以后的数据分布趋向均匀分布。直方图均衡技术被广泛地应用于图像增强。直方图均衡技术由下列文件公开1)Two-dimensional Signal and Image Processing，Prentice hall，Englewoodcliffs，New Jersey，19902)Digital Image Processing，R.C.Gonzalez，P.Wints，Addison-Wesley，Reading，Massachusetts3)Evaluation of the Effectiveness of Adaptive Histogram Equalization forContrast Enhancement，J.Zimmerman，S.prizer，E.Staab，E.Perry，W.McCarteney and B.Brenton，IEEE Trans.on Medical Imaing，PP.304-3124)Application of Adaptive Histogram Equalization to x-ray chest Image.Y.Li，W.Wang and D.Y.Yu，Proc.of the SPIE，PP.513-514，vol.2321，1994这些公开的文件研究了直方图均衡技术在图像处理方面的应用并实践于医用图像和雷达图象的增强。直方图均衡的本质在于使得增强后数据的梯度与该数据的发生几率成正比，在一副图像中，如果某个像数值越多，那么增强后它与相邻的像数值之间的对比度就越强。一般情况下直方图均衡在扩展对比度的同时，扩大了图像的动态范围，使得图像的视觉感受变好，如图4(a)和图4(b)所示，分别是均衡前的图和均衡后的直方图，与图4(a)对比，图4(b)的像素亮度分布不再集中在某个地方，因而更均匀，动态范围更大。但直方图均衡技术也有非常明显的缺点它的均衡结果由数据的分布特征完全决定，均衡过程不受控制。直方图均衡的缺点是增强后的图像的平均亮度都趋向亮度值域的中间值，当原图的平均亮度远离该值时，增强后的图像的平均亮度变化太大，图像的总体视觉也就变化较大，在许多应用场合不符合应用需求。
技术实现思路
本专利技术提供的，通过控制平均亮度变化从而对直方图均衡过程进行控制。为了达到上述目的，本专利技术提供了，其包含3种技术方案，其包含如下步骤步骤a)、输入数字图像，若是彩色数字图像，则需要按颜色模型抽取其中的亮度数值；若是灰度图像，则直接应用灰度信息作为亮度数值；步骤b)、统计亮度直方图H，设为数组H，x∈{Xj|j＝0，1，...，N-1}，其中X0，X1，...，Xj，...，XN-1顺序为N级离散化的图像亮度值，而且X0对应图像亮度值域的极小值，XN-1对应图像亮度值域的极大值；步骤c)、根据直方图，计算亮度累计概率分布函数CDF，设为CDF，x∈{Xj|j＝0，1，...，N-1}，其中X0，X1，...，Xj，...，XN-1顺序为N级离散化的图像亮度值，而且X0对应图像亮度值域的极小值，XN-1对应图像亮度值域的极大值；步骤d)，计算输入图像的平均亮度BA；步骤e)、计算基于平均亮度移动限制的均衡系数K设DBA为预定的允许亮度移动最大范围，计算K＝DBA/|BAH-BA|，如果K＞1，则使K＝1，其中，||表示取绝对值，BAH为直方图均衡后的平均亮度，其为亮度值域的中值，即BAH＝(X0+XN-1)/2；步骤f)、计算亮度映射表，即亮度值域映射表G，x∈{Xj|j＝0，1，...，N-1}，其中X0，X1，...，Xj，...，XN-1顺序为N级离散化的图像亮度值，而且X0对应图像亮度值域的极小值，XN-1对应图像亮度值域的极大值；则G＝(1-K)×Xj+K×CDF×XN-1；步骤g)，设定输入图像的亮度值为Xin，计算增强后的原图像的亮度值为Xout，则Xout＝G；根据输入的图像性质，还原数字图像，若输入的是彩色图像，则还原彩色图像，若输入的是灰度图像，则还原灰度图像。在步骤e)中，对均衡系数K进行限制K＝min(DBA/(BAH-BA)，Kmax)，其中，Kmax为预设的最大直方图均衡系数，这样可以防止当平均值靠近或等于BAH时产生的较大的K值；在步骤e)中，当图像亮度分布较为集中时，平均亮度的移动要小一点，当图像亮度分布离散时，平均亮度的移动可以大一点；故对某帧/场的平均亮度的最大移动DBA进行动态优化设计DBA＝C×Max(Sigma-C0，0)，其中，C是预设常数，Sigma是上一帧/场的亮度标准偏差，C0是预设恒量；当Sigma小于等于该恒量时，即直方图离散程度不大时，DBA为零，K也为零，直方图均衡没有增强效果，标准偏差越大，直方图均衡效果越大；使用离散程度作为控制手段可以获得更好的图像处理鲁棒性。在步骤e)中，对平均亮度的最大移动DBA采用均方差动态设计DBA＝C×Max(MSE-C0，0)，其中，MSE为亮度均方差。在步骤e)中，对平均亮度的最大移动DBA采用简单平均差动态设计DBA＝C×Max(ME-C0，0)，其中，ME为简单平均差，其计算过程如下ME=1NumΣj=0Num-1|Bj-BA|,]]>其中，Num为像素总数，Bj为第j个像素亮度。在所述的步骤a)中，颜色模型可以是YUV颜色空间的Y值，YCbCr空间的Y值，HSV空间的V值，也可以是HSL空间的L值，或者它们的等效表达式。一种控制平均亮度的截断直方图均衡方法，其包含以下步骤步骤a)、输入数字图像，若是彩色数字图像，则需要按颜色模型抽取其中的亮度数值；若是灰度图像，则直接应用灰度信息作为亮度数值；步骤b)、统计直方图H、截断直方图CH、截断像数和CN，其中x∈{Xj|j＝0，1，...，N-1}，其中X0，X1，...，Xj，...，XN-1顺序为N级离散化的图像亮度值，而且X0对应图像亮度值域的极小值，XN-1对应图像亮度值域的极大值； 步骤b1)、初始化以上统计目标的存储空间H＝0；CH＝0；CN＝0，x∈{Xj|j＝0，1，...，N-1}；步骤b2)、按顺序遍历图像，读入图像当前像素的亮度值Xin步骤b3)、H＝H+1；步骤b4)、判断该像素的亮度值所对应的截断直方图存储数组元值是否小于预定的参数CountMax，如果是，则使CH＝CH+1，如果否，则使CN＝CN+1；步骤b5)、判断是否已经对图像所有的像素统计完毕，如果是，则执行步骤c)，如果不是，则返回，循环执行步骤b2)～步骤b5)；步骤c)、计算基于截断直方图的累计概率分布函数CDF，设为CCDF，x∈{Xj|j＝0，1，...，N-1}，其中X0，X1，...，Xj，...，XN-1顺序为N级离散化的图像亮度值，而且X0对应图像亮度值域的极小值，XN-1对应图像亮度值域的极大值；定义数组CCH，x∈{Xj|j＝0，1，...，N-1}，其中X0，X1，...，Xj，...，XN-本文档来自技高网...

【技术保护点】
一种控制平均亮度的直方图均衡方法，其特征在于，包含如下步骤：步骤ａ）、输入数字图像，若是彩色数字图像，则需要按颜色模型抽取其中的亮度数值；若是灰度图像，则直接应用灰度信息作为亮度数值；步骤ｂ）、统计亮度直方图Ｈ，设为数组Ｈ［ ｘ］，ｘ∈｛Ｘ↓［ｊ］｜ｊ＝０，１，…，Ｎ－１｝，其中Ｘ↓［０］，Ｘ↓［１］，…，Ｘ↓［ｊ］，…，Ｘ↓［Ｎ－１］顺序为Ｎ级离散化的图像亮度值，而且Ｘ↓［０］对应图像亮度值域的极小值，Ｘ↓［Ｎ－１］对应图像亮度值域的极大值；步骤ｃ）、 根据直方图，计算亮度累计概率分布函数ＣＤＦ，设为ＣＤＦ［ｘ］，ｘ∈｛Ｘ↓［ｊ］｜ｊ＝０，１，…，Ｎ－１｝，其中Ｘ↓［０］，Ｘ↓［１］，…，Ｘ↓［ｊ］，…，Ｘ↓［Ｎ－１］顺序为Ｎ级离散化的图像亮度值，而且Ｘ↓［０］对应图像亮度值域的极小值，Ｘ↓［Ｎ－１］对应图像亮度值域的极大值；步骤ｄ），计算输入图像的平均亮度ＢＡ；步骤ｅ）、计算基于平均亮度移动限制的均衡系数Ｋ：设ＤＢＡ为预定的允许亮度移动最大范围，计算：Ｋ＝ＤＢＡ／｜ＢＡＨ－ＢＡ｜，如果Ｋ＞１，则使Ｋ＝１；其 中，｜｜表示取绝对值，ＢＡＨ为直方图均衡后的平均亮度；步骤ｆ）、计算亮度映射表，即亮度值域映射表Ｇ［ｘ］，ｘ∈｛Ｘ↓［ｊ］｜ｊ＝０，１，…，Ｎ－１｝，其中Ｘ↓［０］，Ｘ↓［１］，…，Ｘ↓［ｊ］，…，Ｘ↓［Ｎ－１］顺序为Ｎ级离散化的图 像亮度值，而且Ｘ↓［０］对应图像亮度值域的极小值，Ｘ↓［Ｎ－１］对应图像亮度值域的极大值；则：Ｇ［Ｘ↓［ｊ］］＝（１－Ｋ）×Ｘ↓［ｊ］＋Ｋ×ＣＤＦ［Ｘ↓［ｊ］］×Ｘ↓［Ｎ－１］；步骤ｇ）、设定输入图像的亮度值为Ｘｉｎ，增强后的原图像 的亮度值为Ｘｏｕｔ，则Ｘｏｕｔ＝Ｇ［Ｘｉｎ］；根据输入的图像性质，还原数字图像，若输入的是彩色图像，则还原彩色图像，若输入的是灰度图像，则还原灰度图像。...

【技术特征摘要】

【专利技术属性】
技术研发人员：黄晓东，侯钢，王国中，
申请(专利权)人：上海广电集团有限公司中央研究院，
类型：发明
国别省市：31[中国|上海]