1、5）计算映射后的输出图像灰度级g（i）INT（gmaxgmin）C（k）gmin05，i0，1， ，P1，P 为输出图像灰度等级个数，INT 为取整符号；6）统计映射后各灰度等级的像素数目n i，i0，1，P1；7）计算输出图像直方图P（g i）n iN，i0，1，P1；8）用r k和g i的映射关系修改原始图像的灰度级，获得均衡化后的输出图像。为了验证直方图均衡方法是否可以达到增强图像的目的，下面我们就利用在宿舍拍摄的照片来进行验证，图像直方图均衡化前后效果如图1所示。由图1中的a和b可看出均衡后的图像对比度明显增强，亮度增大，原本图中的黑暗区域也可以看到了。由c和d可看出原图的灰度级集中在
2、数值较低的部分而直方图均衡后使得直方图变得均匀分布了，在整体灰度级上都有分布。 图11.2算法存在的问题由于直方图均衡化理论来源于连续函数而数字图像的的灰度是离散值应用于数字图像处理的变换函数进行了从连续到离散的近似。如积分运算就变成了累加运算因而在直方图均衡存在以下问题。1）量化误差会造成原图像信息的丢失，原来的像素的总的灰度值为连续的而变换为离散的这样就存在量化误差如1.1中步骤5）中得取整，如变换后的灰度值为50.15和50.00由于灰度级只有256个所以只能把变换后的两个灰度值都量化到50，那这样就势必会造成原图像的信息丢失。丢失的一定是数量很少的像素。2）无法增强局部细节由于直方图均
3、衡变换是针对整幅图像所有像素值进行统一变换，所以这种算法只能整体增强图像，而不能针对局部细节部分也进行增强。并且当灰度集中在地区域或高区域时这是会造成画面过亮，而整幅图像的灰度的范围没有提升，也达不到突出图像细节的目的。所以这样就提出了自适应直方图均衡化（AHE）。二自适应直方图均衡化 上面的介绍的直方图均衡是图像中的每一点都进行运算,也就是说, 灰度变换函数运算与像素所处的位置无关,这种全局性处理的算法, 它具有算法简单, 计算速度较快等优点, 但由于它是对所有像素点都作同样的处理, 忽略了图像的局部特征, 这样导致经过直方图均衡法算法的图像将丢失图像的有用信息, 对图像的去噪处理及边缘检测
4、带来损失。 那么该如何提取图像的局部特征呢? 采用局部法对图像进行处理, 不同局部采用不同的对比度增强方法, 也就是说, 根据图像象素的局部统计特征来决定处理方法. 每个象素的灰度值都通过一个均衡化变换函数得到的, 而该变换函数是由以该象素中心的一个局部子图像的直方图得到的，称其为局部对比度增强法局部对比度增强法的计算公式为 （2）其中x i, j , xi, j 分布为变换前、后的中心象素,为窗W 内象素的平均灰值,从（2）可看出, 当k 1 时, 如果x i, j mi, j , 那么x i, jx i, j,否则x i, j 1，得到，该领域内有细节存在，而k值也变大，细节被增强。综上所
5、述, 我们得出自适应直方图均衡算法具体的实现过程如下。1） 给出原始图像的所有灰度级i,统计原始图像各灰度级的像素 数ni2） 计算原始图像的直方图与累积直方图。3） 用式（4）计算k值. 用式（3）进行计算局部灰度值。4） 用式（2） 计算局部对比度, 实现均衡化。5） 用p （ti） =计算新的直方图。三对比度受限的自适应直方图均衡 AHE有过度放大图像中相同区域的噪音的问题，另外一种自适应的直方图均衡算法即限制对比度直方图均衡（CLAHE）算法能有限的限制这种不利的放大。3.1算法简介对比度受限自适应直方图均衡法（CLAHE）通过限制局部直方图的高度来限制局部对比度的增强幅度，从而限制噪
6、声的放大及局部对比度的过增强。对比度增强的幅度可定义为灰阶映射函数的斜率。假定自适应直方图均衡方法的滑动窗口的大小为MM，则局部映射函数为： （5）CDF（i）为滑动窗口局部直方图的累积分布函数。累积分布函数CDF（i）的导数为直方图Hist（i），从而局部映射函数m（i）的斜率S为： （6）因此限制直方图高度就等效于限制映射函数m（i）的斜率S，进而限制对比增强度。若限定最大斜率为Smax，则允许的最大直方图高度为： （7）从而，对于高度大于Hmax的直方图应截去多余的那部分，如图1所示。由图可知，实际上是从阈值T（而非Hmax）处对直方图进行截断，然后将截去的部分均匀地分布在整个灰阶范围上
7、，以保证总的直方图面积不变，从而使整个直方图上升一个高度L。因此Hmax、T、L 三者之间应满足下面关系： Hmax=T+L （8）最后，改进的直方图值为：通过改变最大的映射函数斜率Smax及相应的最大直方图高度Hmax，便可以获得不同增强效果的图像。CLAHE可以有效抑制局部对比度的增强及噪声放大。然而，在输出图像中仍然会产生大量的人为噪声，尤其是在灰阶突变的交界区域，这是由于灰阶突变交界区域的局部直方图剧烈变化而引起的。基本步骤：1）图像分块2）产生局部直方图3）分别对每个图像块进行限制对比度直方图均衡，每个图像块分别产生独立的变换函数4）将多个图像块的灰度查找表拼接成输出图像其基本思想是
8、构造限制函数，限制灰度级的概率密度，并将超过限制函数的像素点在直方图内进行重整。具体算法：1）定义限制函数clipLimit，计算超过限制函数的总像素点数totalExcess2）计算totalExcess平均分到每个灰度级的平均值avgBinIncr3）定义限制函数与平均值的差upperLimit = clipLimit - avgBinIncr; 4）对直方图进行限高处理 if 第K灰度级像素数限制函数 降低该灰度级像素数=限制函数 else if第K灰度级像素数 upperLimit 增高该灰度级概率密度=限制函数； 并在超过限制函数的总像素点数totalExcess中减去增加的像素点数
9、 Else 第K灰度级像素数upperLimit 该灰度级像素点数+多余像素点平均值avgBinIncr； 并在超过限制函数的总像素点数totalExcess中减去该平均值。5）对直方图进行均衡3.2仿真结果 图2由上图的仿真结果可以看出来，限制对比度的直方图均衡得到的图片细节更加的突出，没有过亮或过暗的情况出现，弥补了图像直方图均衡在处理低灰度层密集分布的图像时,不能有效增加图像动态范围与对比度的缺陷.该算法与直方图匹配算法相比,算法简单,计算量小,对图像处理系统要求低,容易实现.在图像实时处理领域中可以使用该算法来取代直方图匹配算法,以实现图像直方图均衡的处理。4总结通过这次的数字图像处理
10、学习，我也更好的理解了什么是数字图像处理，并且也学到了很多图像处理的相关知识，包括图像增强、图像切割、图像变换、图像编码、图像分析等好多的知识，并通过本次的实验也熟悉了MATLAB的编程和软件的操作，并且了解了上述三种算法的过程与优劣并且自己仿真，所以说还是收益颇丰的。5附录下面这部分是所涉及到的英文文献的粗略翻译与文章的程序A Novel Approach for Image Enhancement by Using Contrast Limited Adaptive Histogram Equalization Method（一种通过运用对比度受限自适应直方图均衡方法的图像增强方法）摘要一
11、个图像增强的新方法通过运用对比度受限的自适应直方图均衡方法将铲射给一个好的对比度图像如医学图像。在本文中，我们提出了一个新的图像增强方式通过运用对比度首先的自适应直方图均衡的方法。我们提出了一个自适应直方图均衡方法总体的框架。我们将通过对比其他的对比度增强方法来证明本方法的有效性。索引词汇：、介绍图像增强的目的是产生一个重要的步骤在图像处理中，这个步骤是通过编辑原始的图像来观测到更多的图像对比度在一个特定的应用领域中。对比度增强技术将会在图像处理应用中扮演重要的角色，如移动图像，数字照片，分析医疗图像，遥感，和各种图形科学中。所有的图像会由几点因素产生不好的对比度，因为运用不好质量的成像设备或
12、者光线和环境所以我们的提出了一个新的对比度增强技术意在消除这些问题。不同的对比度增强技术被用在提高图像的对比度如直方图均衡，直方图修正，贪婪算法，自适应直方图均衡等等。本文提出了一个新的一个基于对比度受限自适应直方图均衡的直方图均衡方法。退出对比度增强方法现在有很多的对比度增强技术可使用在不清晰图像上，所以我们先讨论下面几中对比度增强技术直方图均衡直方图均衡的目的是分配一个图像的灰度级，从而使没给灰度级的发生概率近乎相等。直方图均衡将会增加亮度和一幅暗的和低对比度图像的对比度。使原始图像中观察到的特征不可见。它还用来规范图像的亮度和对比度，直方图均衡进程是找到一个映射函数来映射，图像输入直方图
13、函数和使分布式输出直方图函数相一致。直方图均衡化还被用在生物神经网络上来是神经元输出功率达到最大做为输入统计函数。这已经在fly retina 中被证明了。直方图均衡化是在更一般的直方图映射方法中较为特别的。这些方法是通过调整图像来时期更容易被分析和增强其显示质量。B）直方图平滑为了避免尖峰带来的强烈的排斥固定点问题一个平滑度约束可以完成这个目标。直方图的向后方差被用来确保其平滑度。平滑可以修改直方图是其尖峰更少因为他们在直方图中基本上都是突然变化的。C）自适应直方图均衡普通的直方图均衡运用的是从图像直方图之得到的一样的变换源来去变换所有得像素。这个在像素值分布相似的图像中可以工作的很好。然而
14、，当图像包含区域比大部分图像部分更暗或更亮时，这些区域的对比度将得不到足够的增强。自适应直方图均衡可以通过运用来自于临近区域的变换函数来变换每一个像素值来改善这一点。它的首次开发是运用在飞机驾驶舱显示器中。当图像区域包含的一个像素的临域是相当均匀的，它的直方图将会特别的尖锐，转移函数将只会映射整个结果图像中的一个狭窄域的像素值。这导致了自适应直方图均衡过于放大了很大程度上同质的的图像上的少量信号。这个方法被用来增强图像的对比度。它通过自适应的方式改变了直方图均衡化而这种方法是通过计算几个直方图来完成的，每个直方图对应着一个图像的不同区域，从而重新分配图像的亮度值。因此它便于增减图像的局部对比度
15、和显示更多的细节。3提出的算法（a对比度受限的自适应直方图均衡一个提出算法为了医学图像特别的开发出来它提供了很好的图像增强相比于原始图像。CLAHE算法把图像分割成些相关区域然后对没个区域进行直方图均衡化。这些事件产生了使用过的灰度值并重新分配因此可以使图像隐藏的特征更加的可见。CLAHE是AHE的改进算法。我们运用上述提到的算法来去增强图像，直方图均衡化，直方图平滑，自适应直方图均衡化和加强了的对比度首先的自适应直方图均衡化。这些提到的加强技术馋了的如图一的结果，第一次测试结果的可视化表示（乳房肿瘤）。从观察到的可视化结果分析所有的算法都能够加强其对比度但是我们提到的算法中对比度受限的自适应
16、直方图可以更准切的显示出图像。（3对于图像质量来说人类的视觉不被认为是基准，所以估计上述提到的算法的质量指标已经被计算出来了通过对比原始图像和输出图像。下图2表示了直方图增强图像等级的映射。为了评价上述提到的图像技术增强方法产生了更好的增强图像。表1揭示了CLAHE对比与其他算法的效果。得出的结果再次更好的揭示了这些增强算法的优劣。其它的算法也产生了不错的图像增强，但是还是没有CLAHE好。算法步骤获得输入：制定行和列的数量然后设定限幅界限和用于直方图的统计特征数量（bin）用于构建对比度增强的转换。更高的数值引起更大的动态范围代价是处理速度变慢。对比度增强技术的限幅界限是从0到1来限制对比度
17、增强。更高的数值可以得到更高的对比度。处理输入：制定限幅界限为归一化值如果必要的话，在图像分裂成区域前补充图像（ pad the imagebefore splitting it into regions）对比每一个过程行和列区域（瓦）从而产生灰度映射值做一个此区域的直方图并用制定的数量的统计特征数量（bin）用限幅界限来限制此直方图，对比度限制程序（由变换公式推出）必须使用与每一个临近区。CLAHE为了防止噪声过度的被放大自适应直方图均衡可以解决它。差值允许效率的显著提高而不影响结果的质量上图4代表的是运用对比度受限自适应直方图均衡时不同限幅界限所产生的图。CLAHE的限幅界限归一化从0到1
18、.更高的数值产生更大的对比度。表1.1代表了不同限幅界限时对比度水平。四结论在此文中，对比度受限的自适应直方图均衡已经运用于乳房癌症图像的对比度增强。对比与现在已经提出的流行的的方法直方图均衡，直方图平滑，自适应直方图均衡他它已经被证可以得到更好的结果了。下面是仿真程序和参考程序直方图均衡：1.i1=imread（f:11.jpg）;i2=rgb2gray（i1）;figure;subplot（131）;imshow（i1）;title（原图subplot（132）;imshow（i2）;灰度图subplot（133）;imhist（i2）;原图直方图i3=histeq（i2）;subplot
19、（221）;灰度图subplot（222）;imshow（i3）;b直方图均衡后的图subplot（223）;c原图直方图subplot（224）;imhist（i3）;d直方图均衡后直方图2.%zhifangtuqunhenghuaclear;x1=imread（1.tifm,n=size（x1）;x2=double（x1）;lenna=zeros（m,n）;lenna_equ=zeros（n,n）;histgram=zeros（256）;cdf=zeros（256）;lenna,map=imread（%get histogramfor i=1:m for j=1:n k=lenna（i,j
20、）; histgram（k）=histgram（k）+1; end end %get cdf cdf（1）=histgram（1）; for i=2:256 cdf（i）=cdf（i-1）+histgram（i）; %run point operation for i=1: for j =1: lenna_equ（i,j）=cdf（k）*256/（m*n）; %生成直方图均衡化后的lenna图 g=mat2gray（lenna_equ）; figure; subplot（2,2,1）;imshow（lenna）; subplot（2,2,2）;imhist（lenna）; ylim（auto
21、subplot（2,2,3）;imshow（g）; subplot（2,2,4）;imhist（g）;） f=imread（ figture;imshow（f）; title（the orignal imageimhist（f）; g=histeq（f,256）;image after equalization %zhifangtupipeichengxu I = imread（pout.tifIeq=histeq（I）;c,x=imhist（Ieq）;I1=imread（tire.tifJ=histeq（I1,c）;subplot（2,2,1）;imshow（I）;Image Asubplo
22、t（2,2,2）;imshow（Ieq）;Image A after Equalizationsubplot（2,2,3）;imshow（I1）;Image Bsubplot（2,2,4）;imshow（J）;Image B after Matching受限的自适应直方图均衡oimg=imread（oimg1=rgb2gray（oimg）;i3=adapthisteq（oimg1）;i4=histeq（oimg1）;imshow（oimg1）;imshow（i4）;b直方图均衡c限制对比度的直方图均衡function CEImage = runCLAHE（Image,Min,Max,NrX,N
23、rY,NrBins,Cliplimit）% Contrast Limited Adaptive Histogram Equalization% by Karel Zuiderveld% The program was reproduced by Alireza Saberi in April 2010% These functions implement Contrast Limited Adaptive Histogram Equalization.% The main routine （CLAHE） expects an input image that is stored contigu
24、ously in% memory; the CLAHE output image overwrites the original input image and has the% same minimum and maximum values （which must be provided by the user）.% This implementation assumes that the X- and Y image resolutions are an integer% multiple of the X- and Y sizes of the contextual regions. A check on various other% error conditions is performed.% % Image - The input/output image% XRes - Image resolution in the X direction% YRes - Image resolution in the Y direction% Min - Minimum greyvalue of input image （also becomes minimum of output image）% Max - Ma
