【图像处理算法】直方图均衡化 |
您所在的位置:网站首页 › 直方图均衡化实验中遇到的问题与解决方法 › 【图像处理算法】直方图均衡化 |
参考文献: Rafael C. Gonzalez, Richard E. Woods,Digital Image Processing (Third Edition)胡学龙. 数字图像处理(第三版)左飞. 图像处理中的数学修炼目录 直方图均衡化的介绍 直方图的概念 直方图均衡化的理论基础 手工实现直方图均衡化 MATLAB上实现直方图均衡化 直方图均衡化的缺点 直方图均衡化的介绍直方图均衡化是一种简单有效的图像增强技术,通过改变图像的直方图来改变图像中各像素的灰度,主要用于增强动态范围偏小的图像的对比度。原始图像由于其灰度分布可能集中在较窄的区间,造成图像不够清晰。例如,过曝光图像的灰度级集中在高亮度范围内,而曝光不足将使图像灰度级集中在低亮度范围内。采用直方图均衡化,可以把原始图像的直方图变换为均匀分布(均衡)的形式,这样就增加了像素之间灰度值差别的动态范围,从而达到增强图像整体对比度的效果。换言之,直方图均衡化的基本原理是:对在图像中像素个数多的灰度值(即对画面起主要作用的灰度值)进行展宽,而对像素个数少的灰度值(即对画面不起主要作用的灰度值)进行归并,从而增大对比度,使图像清晰,达到增强的目的。举个例子,如图1所示,左图为原始图像,右图为直方图均衡化后的图像。 图1 直方图均衡化前后效果对比 直方图的概念对一幅灰度图像,其直方图反映了该图像中不同灰度级出现的统计情况。图2给出了一个直方图的示例,其中图(a)是一幅图像,其灰度直方图可表示为图(b),其中横轴表示图像的各灰度级,纵轴表示图像中各灰度级像素的个数。(需要注意,灰度直方图表示了在图像中各个单独灰度级的分布,而图像对比度则取决于相邻近像素之间灰度级的关系。) 图 2 图像及其直方图在MATLAB中可以使用函数imhist()来计算灰度图像的直方图。我们计算一下图1中左图的灰度直方图,结果如图3所示。 I = imread('pout.tif'); figure; subplot(121), imshow(I); subplot(122), imhist(I); axis([0, 256, 0, 4000]); 图 3 灰度图像及其直方图严格地说,图像的灰度直方图是一个一维的离散函数,可写成: (公式1) 式中,是图像中灰度级为的像素的个数。直方图的每一列(称为bin)的高度对应。直方图提供了原图中各种灰度值分布的情况,也可以说直方图给出了一幅图像所有灰度值的整体描述。直方图的均值和方差也是图像灰度的均值和方差。图像的视觉效果与其直方图有对应关系,或者说,直方图的形状和改变对图像有很大的影响。 在直方图的基础上,进一步定义归一化的直方图为灰度级出现的相对频率。即: (公式2) 式中,表示图像的像素的总数,是图像中灰度级为的像素的个数。 我们以图1中左图为例,计算其归一化的直方图,结果如图4所示。 I = imread('pout.tif'); N = numel(I); % 求图像像素的总数 Pr = imhist(I) / N; k = 0 : 255; figure; subplot(121), imshow(I); subplot(122), bar(k, Pr); 图 4 灰度图像及其归一化直方图 直方图均衡化的理论基础为讨论方便起见,以 r 和 s 分别表示归一化了的原图像灰度和经直方图均衡化后的图像灰度(因为归一化了,所以 r 和 s 的取值在0到1之间)。当 r = s = 0时,表示黑色;当 r = s = 1时,表示白色;当 r, s ∈(0, 1)时,表示像素灰度在黑白之间变化。(所谓直方图均衡化,其实是根据直方图对像素点的灰度值进行变换,属于点操作范围。换言之,即:已知r,求其对应的s。) 在 [0,1] 区间内的任何一个 r ,经变换函数 T(r) 都可以产生一个对应的 s ,且 (公式3) 式中,T(r) 应当满足以下两个条件: 在 0 ≤ r ≤ 1 内,T(r) 为单调递增函数;(此条件保证了均衡化后图像的灰度级从黑到白的次序不变)在 0 ≤ r ≤ 1 内有 0 ≤ T(r) ≤ 1。(此条件保证了均衡化后图像的像素灰度值在允许的范围内)公式3的逆变换关系为: (公式4) 式中, 对 s 同样满足上述的两个条件。 由概率论可知,如果已知随机变量 r 的概率密度是 ,而随机变量 s 是 r 的函数,则 s 的概率密度 可以由 求出。假定随机变量 s 的分布函数用 表示,根据分布函数的定义有 (公式5) 又因为概率密度函数是分布函数的导数,因此公式5两边对 s 求导可得: (公式6) 从公式6可以看出,通过变换函数 T(r) 可以控制图像灰度级的概率密度函数 ,从而改善图像的灰度层次,这就是直方图均衡化的理论基础。 又有:从人眼视觉特性来考虑,一幅图像的灰度直方图如果是均匀分布的,那么该图像看上去效果比较好(参考冈萨雷斯数字图像处理3.3节)。因此要做直方图均衡化,这里的 应当是均匀分布的概率密度函数。 由概率论知识可知,对于区间 [a,b]上的均匀分布,其概率密度函数等于 。 如果原图像没有进行归一化,即 , 那么,归一化之后 ,所以这里的 。
由公式6可以知道 ,又因为 ,所以有 。对这个式子两边积分得: (公式7) 公式7就是我们所求的变换函数 。它表明当变换函数 是原图像直方图的累积分布概率时,能达到直方图均衡化的目的。 对于灰度级为离散的数字图像,用频率来代替概率,则变换函数 的离散形式可以表示为: (公式8) 式中,, (注:这里的 ,表示归一化后的灰度级;k表示归一化前的灰度级)。由公式8可以知道,均衡化后各像素的灰度级 可直接由原图像的直方图算出来。需要说明的是,这里的 也是归一化后的灰度级,其值在 0 到 1 之间;有时需要将其乘以再取整,使其灰度级范围在 0 到 L-1之间,与原图像一致。 手工实现直方图均衡化了解直方图均衡化的原理之后,我们以一个简单的例子来手工计算均衡化后的图像。这里我们假设存在以下这张图像(假定图像的灰度级范围是 [0, 9]): 图 5 原始图像计算过程如下: 第一步,计算原始图像的灰度直方图 。 ,即原始图像中灰度级为0的像素的个数是3,直接从图5中可以数出来。按这样的方式分别得出 , , ..., 。这里我们可以得到 。 图 6 原始图像的灰度直方图第二步,计算原始图像的像素总个数。 这里 。 第三步,计算原始图像的灰度分布频率。 , 第四步,计算原始图像的灰度累积分布频率。 , 第五步,将归一化的 乘以 再四舍五入,以使得均衡化后图像的灰度级与归一化前的原始图像一致。 ,四舍五入之后其值为1,也就是说原始图像中灰度级0对应均衡化后的灰度级1,即0→1。 同理,,四舍五入之后为2,即1→2;,四舍五入之后为3,即2→3;,四舍五入之后为5,即3→5;,四舍五入之后为5,即4→5;,四舍五入之后为5,即5→5;,四舍五入之后为7,即6→7;,四舍五入之后为7,即7→7;,四舍五入之后为8,即8→8;,四舍五入之后为9,即9→9。 以上的映射关系,就是变换函数 的作用。 第六步,根据以上映射关系,参照原始图像中的像素,可以写出直方图均衡化之后的图像,如图7所示。 图 7 直方图均衡化之后的图像以上过程即为手动计算直方图均衡化。接着,我们看一下均衡化后图像的灰度直方图,如图8所示。 图 8 均衡化后图像的灰度直方图根据图6和图8,说明一下直方图均衡化是如何增强图像对比度的。在图6中,原始图像灰度值为4,5,7的像素的个数为1,因此在图8中,这三个像素值点分别归并到相邻的灰度值中。因为有三个灰度值归并,因此在均衡化处理后,出现了三个空位,由这些空位将原来相邻的灰度值展开(举个例子:5和6相邻,均衡化后,变成5和7相邻),故而展宽了对比度,但是归并也带来了某些相邻像素对比度的降低(举个例子:4和5相邻,均衡化后,变成5和5相邻)。这也说明了直方图均衡化方法对灰度分布比较集中的图像的处理效果比较明显。 MATLAB上实现直方图均衡化在MATLAB中提供了现成的函数histeq()来实现灰度图像的直方图均衡化,如下例所示: close all; clear; clc; I = imread('pout.tif'); J = histeq(I,256); imshowpair(I,J,'montage'); 图 9 使用函数histeq()进行直方图均衡化但为了演示说明算法的原理,下面将在MATLAB中自行编码实现灰度图像的直方图均衡化。通过代码来演示这个算法显然更加直观,更加易懂。 close all; clear; clc; % 首先读入灰度图像,并提取图像的高度和宽度 image = imread('pout.tif'); [height, width] = size(image); % 然后统计每个灰度的像素值的累计数目 NumPixel = zeros(1,256); % 建立一个256列的行向量,以统计各灰度级的像素个数 for i = 1 : height for j = 1 : width k = image(i,j); % k是像素点(i,j)的灰度值 % 因为NumPixel数组的下标是从1开始的,但是图像像素的取值范围是0~255 % 所以用NumPixel(k+1) NumPixel(k+1) = NumPixel(k+1) + 1; % 对应灰度值像素点数量加1 end end % % 这里我们将数组NumPixel显示出来,以观测效果 % figure; % subplot(121), imshow(image); % subplot(122), bar(NumPixel); % 灰度图像的直方图可以正确显示出来 % 接下来,将频数值算为频率 ProbPixel = zeros(1,256); % 统计各灰度级出现的频率 for i = 1 : 256 ProbPixel(i) = NumPixel(i) / (height * width); end % % 这里我们将数组ProbPixel显示出来,以观测效果 % figure; % subplot(121), imshow(image); % subplot(122), bar(ProbPixel); % 灰度图像的归一化直方图可以正确显示出来 % 再用函数cumsum()来计算累积分布函数(CDF),并将频率(取值范围是0~1)映射到0~255的无符号整数 CumPixel = cumsum(ProbPixel); % 这里的数组CumPixel大小也是1×256 CumPixel = uint8((256-1) .* CumPixel + 0.5); % % 这里我们将数组CumPixel显示出来,以观测效果 % figure; % subplot(121), imshow(image); % subplot(122), bar(CumPixel); % 数组CumPixel可以正确显示出来 % 在下列用作直方图均衡化实现的赋值语句右端,image(i,j)被用来作为CumPixel的索引 % 例如,image(i,j)=120,则从CumPixel中取出第120个值作为image(i,j)的新像素值 outImage = uint8(zeros(height, width)); % 预分配数组 for i = 1 : height for j = 1 : width outImage(i,j) = CumPixel(image(i,j)); end end % 显示直方图均衡化前后的图像,可以发现,与调用函数histeq()的效果一致 imshowpair(image, outImage, 'montage');实验效果如图10所示。 图 10 自行编码实现的直方图均衡化效果上述讨论的是灰度图像的直方图均衡化。对于彩色图像而言,可以分别对R、G、B三个分量来做直方图均衡化,这也确实是一种方法。但有些时候,这样做很有可能会导致结果图像色彩失真。因此有人建议将RGB空间转换为HSV之后,对V分量进行直方图均衡化,以保证图像色彩不失真。(HSV分别指色调、饱和度、亮度)。 下面我们采用图像处理工具箱中的测试用图(图11)分别做RGB空间和HSV空间的直方图均衡化。 图 11 MATLAB自带的图像首先在RGB空间进行直方图均衡化处理。这里为了简便,直接调用MATLAB函数histeq()。 close all; clear; clc; I = imread('baby.jpg'); % figure,imshow(I); % 分别提取R、G、B三个分量 R = I(:, :, 1); G = I(:, :, 2); B = I(:, :, 3); % 分别对三个分量进行直方图均衡化 R = histeq(R, 256); G = histeq(G, 256); B = histeq(B, 256); J = I; J(:, :, 1) = R; J(:, :, 2) = G; J(:, :, 3) = B; imshowpair(I, J, 'montage');在RGB空间进行直方图均衡化的前后图像如图12所示。 图 12 RGB空间直方图均衡化前后效果对比接下来,我们在HSV空间对V分量进行直方图均衡化处理。这里的代码可以采用调用MATLAB函数histeq()的方式,但是我们使用自行编码的方式进行处理。 close all; clear; clc; I = imread('baby.jpg'); % 将RGB空间转换为HSV空间 hsvImage = rgb2hsv(I); % 提取V分量 v = hsvImage(:, :, 3); % 这里的v是double类型的矩阵 % 以下代码与前面介绍基本一致,这里不再做过多注释 [height, width] = size(v); v = uint8(v .* 255 + 0.5); % 这里的0.5有必要加上,以免矩阵v中出现0 N = zeros(1, 256); for i = 1 : height for j = 1 : width k = v(i,j); N(k+1) = N(k+1) + 1; end end ProbPixel = zeros(1, 256); for i = 1 : 256 ProbPixel(i) = N(i) / (height * width); end CumPixel = cumsum(ProbPixel); CumPixel = uint8(255 .* CumPixel + 0.5); % 四舍五入 for i = 1 : height for j = 1 : width v(i,j) = CumPixel(v(i,j)); % 这里的v(i,j)不能为0,否则数组索引出错 end end v = im2double(v); hsvImage(:, :, 3) = v; outImage = hsv2rgb(hsvImage); imshowpair(I, outImage, 'montage');在HSV空间进行直方图均衡化的前后图像如图13所示。 图 13 HSV空间直方图均衡化前后效果对比 直方图均衡化的缺点如果一幅图像整体偏暗或者偏亮,那么直方图均衡化的方法很适用。但直方图均衡化是一种全局处理方式,它对处理的数据不加选择,可能会增加背景干扰信息的对比度并且降低有用信号的对比度(如果图像某些区域对比度很好,而另一些区域对比度不好,那采用直方图均衡化就不一定适用)。此外,均衡化后图像的灰度级减少,某些细节将会消失;某些图像(如直方图有高峰),经过均衡化后对比度不自然的过分增强。针对直方图均衡化的缺点,已经有局部的直方图均衡化方法出现。 |
今日新闻 |
推荐新闻 |
CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |