图像增强直方图均衡化PPT教案

1、图像增强直方图均衡化图像增强直方图均衡化1. . 灰度级直方图的概念灰度级直方图的概念2. . 直方图的计算直方图的计算3. . 直方图的性质直方图的性质4. . 直方图均衡化直方图均衡化第1页/共26页长江大学电子信息学院长江大学电子信息学院灰度级直方图灰度级直方图1 直方图的概念直方图的概念2p灰度级的直方图是反应一幅图像中的灰度级与出现这种灰度的概率之间的图形。p设图像的灰度范围为a,b，r为此灰度范围内的任一灰度级，p(r)为这幅图像中灰度级为r的像素出现的频率，可以看出p(r)是r的函数。该函数的图形称为这幅图像的直方图。( )rp r 灰度为 的像素数图像上的总像素数1)(1nii

2、rp第2页/共26页长江大学电子信息学院长江大学电子信息学院灰度级灰度级123456灰度级的像素数灰度级的像素数7437213图像的直方图图像的直方图第3页/共26页长江大学电子信息学院长江大学电子信息学院第4页/共26页长江大学电子信息学院长江大学电子信息学院设设r表示图像中像素的灰度级，可以用概率密度函数表示图像中像素的灰度级，可以用概率密度函数Pr(rk)表表示原始图像的灰度分布。则示原始图像的灰度分布。则( )/rkkP rnN式中，式中，N为一幅图像中像素的总数；为一幅图像中像素的总数；nk为第为第k级灰度的像素；级灰度的像素；rk为第为第k个灰度级；个灰度级；Pr(rk)表示该灰度

3、级出现的概率。因为表示该灰度级出现的概率。因为Pr(rk)给出了对给出了对rk出现概率的一个估计，所以直方图提供了原图的灰出现概率的一个估计，所以直方图提供了原图的灰度值分布情况，也可以说给出了一幅图像所有灰度值的整体描度值分布情况，也可以说给出了一幅图像所有灰度值的整体描述。述。 第5页/共26页长江大学电子信息学院长江大学电子信息学院（1 1） 它只反映该图像中不同灰度值出现的次数，而未它只反映该图像中不同灰度值出现的次数，而未反映某一灰度值像素所在的位置。也就是说，它只包含反映某一灰度值像素所在的位置。也就是说，它只包含了该图像中某一灰度值的像素出现的概率，而丢失了其了该图像中某一灰度值

4、的像素出现的概率，而丢失了其所在位置的信息。所在位置的信息。（2 2） 任一幅图像，都能唯一地确定出一幅与它对应的任一幅图像，都能唯一地确定出一幅与它对应的直方图，但不同的图像，可能有相同的直方图。直方图，但不同的图像，可能有相同的直方图。第6页/共26页长江大学电子信息学院长江大学电子信息学院第7页/共26页长江大学电子信息学院长江大学电子信息学院1. 大多数自然图像，其灰度分布集中在较窄的区间，引起图像大多数自然图像，其灰度分布集中在较窄的区间，引起图像细节不够清晰，采用直方图修正后可使图像的灰度间距拉开或使细节不够清晰，采用直方图修正后可使图像的灰度间距拉开或使灰度分布均匀，从而增大反差

5、，使图像细节清晰，达到增强的目灰度分布均匀，从而增大反差，使图像细节清晰，达到增强的目的。例如一幅过曝光的图片，其灰度级都集中在高亮度范围内，的。例如一幅过曝光的图片，其灰度级都集中在高亮度范围内，而曝光不足的图片，其灰度级集中在低亮度范围内，具有这样直而曝光不足的图片，其灰度级集中在低亮度范围内，具有这样直方图的图片其可视效果比较差。方图的图片其可视效果比较差。 第8页/共26页长江大学电子信息学院长江大学电子信息学院(1) 对给定的待处理图像统计其直方图，求出对给定的待处理图像统计其直方图，求出(2) 根据统计出的直方图采用累积分布函数作变换根据统计出的直方图采用累积分布函数作变换,0()

6、()kkkrjjST rP r (3) 用新灰度代替旧灰度，求出用新灰度代替旧灰度，求出Ps(s)，这一步是近似过程，应，这一步是近似过程，应根据处理目的尽量做到合理，同时把灰度值相等或近似地合根据处理目的尽量做到合理，同时把灰度值相等或近似地合并到一起。并到一起。( )/rkkP rnN求变换后的新灰度；求变换后的新灰度；第9页/共26页长江大学电子信息学院长江大学电子信息学院例子例子 假定有一幅像素数为6464，灰度级为8级的图像，将其进行均衡化处理，其灰度级分布表如右图所示：表表 64646464大小的图像灰度级分布大小的图像灰度级分布第10页/共26页长江大学电子信息学院长江大学电子信

7、息学院图图(a) (a) 原图灰度级直方图；图原图灰度级直方图；图(b) (b) 累积变换后的直方图；累积变换后的直方图；图图(c) (c) 均衡化后的直方图；均衡化后的直方图；第11页/共26页长江大学电子信息学院长江大学电子信息学院3. （1）imhist函数函数功能功能：计算和显示图像的色彩直方图。：计算和显示图像的色彩直方图。 格式格式：imhist(I, n) imhist(X, map) counts, x=imhist() 说明说明：imhist(I, n)计算和显示灰度图像计算和显示灰度图像I的直方图，的直方图，n为指定为指定的灰度级数目，缺省值为的灰度级数目，缺省值为256；

8、 imhist(X, map)计算和显示索引色图像计算和显示索引色图像X的直方图，的直方图，map为调为调色板；色板； counts, x=imhist(.) 返回直方图数据向量返回直方图数据向量counts和相应的和相应的色彩值向量色彩值向量x，用，用stem(x, counts) 同样可以显示直方图。同样可以显示直方图。第12页/共26页长江大学电子信息学院长江大学电子信息学院 显示灰度图像显示灰度图像baby.jpg的直方图的直方图 程序：程序：I=imread(C:baby.jpg);subplot(1,2,1),imshow(I);subplot(1,2,2),imhist(I);3

9、.第13页/共26页长江大学电子信息学院长江大学电子信息学院灰度图像的色彩直方图灰度图像的色彩直方图第14页/共26页长江大学电子信息学院长江大学电子信息学院显示索引色图像的直方图显示索引色图像的直方图程序：程序： clf;load trees; subplot(1,2,1),imshow(I); subplot(122);imhist(X,map);3.第15页/共26页长江大学电子信息学院长江大学电子信息学院索引色图像的直方图索引色图像的直方图第16页/共26页长江大学电子信息学院长江大学电子信息学院程序：程序： I=imread(trees.tif); c,x=imhist(I); su

10、bplot(121);imshow(I); subplot(1,2,2),stem(x,c);3.用用stem函数显示图像直方图函数显示图像直方图第17页/共26页长江大学电子信息学院长江大学电子信息学院图像的直方图图像的直方图第18页/共26页长江大学电子信息学院长江大学电子信息学院（2）imadjust函数函数功能：通过直方图变换调整对比度。 格式：J=imadjust(I, low high,bottom top, gamma) newmap=imadjust(map, low high, bottom top, gamma) 说明：J=imadjust(I, low high, bot

11、tom top, gamma)返回图像I经直方图调整后的图像J，gamma为校正量 ，low high为原图像中要变换的灰度范围，bottom top指定了变换后的灰度范围；newmap=imadjust(map, low highbottom top, gamma)调整索引色图像的调色板map。此时若low high和bottom top都为23的矩阵，则分别调整R、G、B 3个分量。 3.第19页/共26页长江大学电子信息学院长江大学电子信息学院 程序：程序： clear all I=imread(C:lena.bmp); J=imadjust(I,0.3 0.7, ); subplot(

12、221), imshow(I)； subplot(222), imshow(J)； subplot(223), imhist(I)； subplot(224), imhist(J)；调整图像的对比度，调整前后的图像见图调整图像的对比度，调整前后的图像见图3.第20页/共26页长江大学电子信息学院长江大学电子信息学院第21页/共26页长江大学电子信息学院长江大学电子信息学院3. （3）histeq函数函数功能功能：直方图均衡化直方图均衡化 格式格式：J=histeq(I, n)，指定均衡化后灰度级数指定均衡化后灰度级数n n，缺省，缺省 为为6464；第22页/共26页长江大学电子信息学院长江大学电子信息学院 程序：程序： I = imread(cameraman.tif); J = histeq(I); subplot(2,2,1),imshow(I); subplot(2,2,2