数字图像处理作业

基于直方图的图像增强方法姓名：XXX专业：电子与通信工程学号：XXX导师：XXX基于直方图的图像增强方法一．直方图均衡化算法直方图均衡化是图像处理领域中利用图像直方图对对比度进行调整的方法，这种方法通常用来增加许多图像的局部对比度，尤其是当图像的有用数据的对比度相当接近的时候。通过这种方法，亮度可以更好地砸直方图上分布。这样就额可以拥有增强局部的对比度而不影响整体的对比度，直方图均衡化通过有效地扩展常用的亮度来实现这种功能。1.1算法介绍与仿真图像的直方图是图像的重要统计特征，是用来表征数字图像的每一灰度级与其出现的频率间的统计关系的方法。其数学公式如式（1）所示：)1,....2,1,0()r(kLkNnPk(1)式中，P（rk）为图像f（x，y）的第k级灰度出现的概率；rk为第k级灰度的灰度值级；nk为图像中灰度值为rk的像素的个数；N为图像f（x，y）的总像素数；L为总的灰度级数，需要注意的是直方图不表示图像的空间信息，且任一特定图像都有唯一直方图，但反之并不成立在图像直方图中，整体较暗的图像其直方图的组成部分集中在灰度低的一侧，而明亮的图像的直方图组成部分集中在灰度较高的一侧。当一副图像其像素占据全部可能的灰度范围且分布均匀时，该图像有较高的对比度，其图像也比较清晰。因此，可以通过改变直方图的灰度级分布达到增强图像的目的。直方图均衡化是一种以累积分布函数变换法为基础的直方图修正方法。其基本思路是将一已知灰度概率分布的图像经过变换，使之成为具有均匀灰度概率分布、输出图像的直方图近似服从均匀分布的变换算法。其计算步骤如下：1）列出原始图像的灰度级rk，k＝0，1，2，…，L－1，L为总的灰度级数；2）统计各灰度级像素数目nk，k＝0，1，2，…，L－1；3）计算原始图像直方图各灰度级的频率P（rk）＝nk／N，k＝0，1，2，…，L－1，N为图像f（x，y）的总像素数；4）计算累计分布函数1,...2,1,0,)()k(10LkrPCLk;5）计算映射后的输出图像灰度级g（i）＝INT［（gmax－gmin）C（k）＋gmin＋0．5］，i＝0，1，…，P－1，P为输出图像灰度等级个数，INT为取整符号；6）统计映射后各灰度等级的像素数目ni，i＝0，1，…，P－1；7）计算输出图像直方图P（gi）＝ni／N，i＝0，1，…，P－1；8）用rk和gi的映射关系修改原始图像的灰度级，获得均衡化后的输出图像。为了验证直方图均衡方法是否可以达到增强图像的目的，下面我们就利用在宿舍拍摄的照片来进行验证，图像直方图均衡化前后效果如图1所示。由图1中的a和b可看出均衡后的图像对比度明显增强，亮度增大，原本图中的黑暗区域也可以看到了。由c和d可看出原图的灰度级集中在数值较低的部分而直方图均衡后使得直方图变得均匀分布了，在整体灰度级上都有分布。图11.2算法存在的问题由于直方图均衡化理论来源于连续函数而数字图像的的灰度是离散值应用于数字图像处理的变换函数进行了从连续到离散的近似。如积分运算就变成了累加运算因而在直方图均衡存在以下问题。1）量化误差会造成原图像信息的丢失，原来的像素的总的灰度值为连续的而变换为离散的这样就存在量化误差如1.1中步骤5）中得取整，如变换后的灰度值为50.15和50.00由于灰度级只有256个所以只能把变换后的两个灰度值都量化到50，那这样就势必会造成原图像的信息丢失。丢失的一定是数量很少的像素。2）无法增强局部细节由于直方图均衡变换是针对整幅图像所有像素值进行统一变换，所以这种算法只能整体增强图像，而不能针对局部细节部分也进行增强。并且当灰度集中在地区域或高区域时这是会造成画面过亮，而整幅图像的灰度的范围没有提升，也达不到突出图像细节的目的。所以这样就提出了自适应直方图均衡化（AHE）。二．自适应直方图均衡化上面的介绍的直方图均衡是图像中的每一点都进行运算,也就是说,灰度变换函数运算与像素所处的位置无关,这种全局性处理的算法,它具有算法简单,计算速度较快等优点,但由于它是对所有像素点都作同样的处理,忽略了图像的局部特征,这样导致经过直方图均衡法算法的图像将丢失图像的有用信息,对图像的去噪处理及边缘检测带来损失。那么该如何提取图像的局部特征呢?采用局部法对图像进行处理,不同局部采用不同的对比度增强方法,也就是说,根据图像象素的局部统计特征来决定处理方法.每个象素的灰度值都通过一个均衡化变换函数得到的,而该变换函数是由以该象素中心的一个局部子图像的直方图得到的，称其为局部对比度增强法局部对比度增强法的计算公式为（2）其中xi,j,xi,j分布为变换前、后的中心象素,wjijijixnmm),(,,1为窗W内象素的平均灰值,从（2）可看出,当k1时,如果xi,jmi,j,那么xi,jxi,j,否则xi,jxi,j。我们并没有调节整幅图像的动态范围，而只是改变局部窗口内的对比度,同样达到了细节增强的效果,优化了图像质量.下面我们详细地分析局部对比度增强法的具体算法。为了保存图像中的细节部分,在实施直方图均衡化之前,先对原图的细节部分保存,再到直方图均衡时将这些细节添入算法过程。于是,图像均衡算法修改为:它　　　　　　　　　其　　　)(2550))((,,,,,,jijijijijijixTxmxkxTx（3）其中xi,j和xi,j表示变换前后的图像灰度值,mi,j表示为以xi,j为中心的窗邻域均值,T表示对xi,j的变换函数。从式（3）中我们可以看出,T起到了调节直方图动态范围的作用,k(xi,j-mi,j)则相当于一个高通滤波器,起到了增强细节(即局部对比度)的作用,强化细节的同时,也在增加高频噪声。自然想到,我们该如何即可增强细节又能避免增加高频噪声呢?其实,要达到这种效果比较简单,算法关键就在k这个)(,,,,jijijijimxkmx1)(2,2njikk自适应参量的选取上。考虑到参量k的自适应性,我们先对k进行分析。当窗W的中心象素xi,j位于无图像细节处时，k趋向于0，而当窗W的中心象素xi,j位于图像细节处时，k取一个较大正值。基于以上考虑我们选取窗W内邻域灰度方差作为自适应变量,那么的表达式可以写为:（4）其中k为比例系数，n2为整幅图像的噪声方差，j,i2为窗W内的灰度方差。若令njir,，则式（2）可表为k=k[r2-1],由于图像窗W内的灰度方差j,i2不小于图像的噪声方差，可知r≥1,当r=1时，k=0，也即nj,i那么图像无细节存在，因此该领域中点的像素值只能增强，当r越大，k值也越大，由此可知r1，得到nj,i，该领域内有细节存在，而k值也变大，细节被增强。综上所述,我们得出自适应直方图均衡算法具体的实现过程如下。1)给出原始图像的所有灰度级i,统计原始图像各灰度级的像素数ni2)计算原始图像的直方图与累积直方图。3)用式（4）计算k值.用式（3）进行计算局部灰度值。4)用式（2）计算局部对比度,实现均衡化。5)用p(ti)=nni计算新的直方图。三．对比度受限的自适应直方图均衡AHE有过度放大图像中相同区域的噪音的问题，另外一种自适应的直方图均衡算法即限制对比度直方图均衡（CLAHE）算法能有限的限制这种不利的放大。3.1算法简介对比度受限自适应直方图均衡法（CLAHE）通过限制局部直方图的高度来限制局部对比度的增强幅度，从而限制噪声的放大及局部对比度的过增强。对比度增强的幅度可定义为灰阶映射函数的斜率。假定自适应直方图均衡方法的滑动窗口的大小为M×M，则局部映射函数为：MMiCDFim)(255)((5)CDF(i)为滑动窗口局部直方图的累积分布函数。累积分布函数CDF(i)的导数为直方图Hist(i)，从而局部映射函数m(i)的斜率S为：MMiHistdiimdS255)()((6)因此限制直方图高度就等效于限制映射函数m(i)的斜率S，进而限制对比增强度。若限定最大斜率为Smax，则允许的最大直方图高度为：255maxmaxMMSH(7)从而，对于高度大于Hmax的直方图应截去多余的那部分，如图1所示。由图可知，实际上是从阈值T（而非Hmax）处对直方图进行截断，然后将截去的部分均匀地分布在整个灰阶范围上，以保证总的直方图面积不变，从而使整个直方图上升一个高度L。因此Hmax、T、L三者之间应满足下面关系：Hmax=T+L(8)最后，改进的直方图值为：　　　　　　　　　　　　　　　TiHistHTiHistLiHistiHist)()()()(max'通过改变最大的映射函数斜率Smax及相应的最大直方图高度Hmax，便可以获得不同增强效果的图像。CLAHE可以有效抑制局部对比度的增强及噪声放大。然而，在输出图像中仍然会产生大量的人为噪声，尤其是在灰阶突变的交界区域，这是由于灰阶突变交界区域的局部直方图剧烈变化而引起的。基本步骤：1)图像分块2)产生局部直方图3)分别对每个图像块进行限制对比度直方图均衡，每个图像块分别产生独立的变换函数4)将多个图像块的灰度查找表拼接成输出图像其基本思想是构造限制函数，限制灰度级的概率密度，并将超过限制函数的像素点在直方图内进行重整。具体算法：1）定义限制函数clipLimit，计算超过限制函数的总像素点数totalExcess2）计算totalExcess平均分到每个灰度级的平均值avgBinIncr3）定义限制函数与平均值的差upperLimit=clipLimit-avgBinIncr;4）对直方图进行限高处理●if第K灰度级像素数限制函数●降低该灰度级像素数=限制函数●else●if第K灰度级像素数upperLimit●增高该灰度级概率密度=限制函数；●并在超过限制函数的总像素点数totalExcess中减去增加的像素点数●Else第K灰度级像素数upperLimit●该灰度级像素点数+多余像素点平均值avgBinIncr；●并在超过限制函数的总像素点数totalExcess中减去该平均值。5）对直方图进行均衡3.2仿真结果图2由上图的仿真结果可以看出来，限制对比度的直方图均衡得到的图片细节更加的突出，没有过亮或过暗的情况出现，弥补了图像直方图均衡在处理低灰度层密集分布的图像时,不能有效增加图像动态范围与对比度的缺陷.该算法与直方图匹配算法相比,算法简单,计算量小,对图像处理系统要求低,容易实现.在图像实时处理领域中可以使用该算法来取代直方图匹配算法,以实现图像直方图均衡的处理。四．总结通过这次的数字图像处理学习，我也更好的理解了什么是数字图像处理，并且也学到了很多图像处理的相关知识，包括图像增强、图像切割、图像变换、图像编码、图像分析等好多的知识，并通过本次的实验也熟悉了MATLAB的编程和软件的操作，并且了解了上述三种算法的过程与优劣并且自己仿真，所以说还是收益颇丰的。五．附录下面这部分是所涉及到的英文文献的粗略翻译与文章的程序ANovelApproachforImageEnhancementbyUsingContrastLimitedAdaptiveHistogramEqualizationMethod（一种通过运用对比度受限自适应直方图均衡方法的图像增强方法）摘要一个图像增强的新方法通过运用对比度受限的自适应直方图均衡方法将铲射给一个好的对比度图像如医学图像。在本文中，我们提出了一个新的图像增强方式通过运用对比度首先的自适应直方图均衡的方法。我们提出了一个自适应直方图均衡方法总体的框架。我们将通过对比其他的对比度增强方法来证明本方法的有效性。索引词汇：、、、、、介绍图像增强的目的是产生一个重要的步骤在图像处理中，这个步骤是通过编辑原始的图像来观测到更多的图像对比度在一个特定的应用领域中。对比度增强技术将会在图像处理应用中扮演重要的角色，如移动图像，数字照片，分析医疗图像，遥感，和各种图形科学中。所有的图像会由几点因素产生不好的对比度，因为运用不好质量的成像设备或者光线和环境