直方图均衡化实际是缩小原图像中小概率像素的灰度等级范围,而拉长大概率像素的灰度等级范围,使得主要像素(有用像素)的对比度增加并减小极少像素(往往为无用数据)的灰度等级范围。以下是维基百科对直方图均衡化的叙述,原文转载自:http://zh.wikipedia.org/wiki/直方图均衡化。
直方图均衡化用来增强图像的全局对比度,尤其是当图像的有用数据的对比度相当接近的时候。通过这种方法,亮度可以更好地在直方图上分布。这样就可以用于增强局部的对比度而不影响整体的对比度,直方图均衡化通过有效地扩展常用的亮度来实现这种功能。
这种方法对于背景和前景都太亮或者太暗的图像非常有用,这种方法尤其是可以带来X光图像中更好的骨骼结构显示以及曝光过度或者曝光不足照片中更好的细节。这种方法的一个主要优势是它是一个相当直观的技术并且是可逆操作,如果已知均衡化函数,那么就可以恢复原始的直方图,并且计算量也不大。这种方法的一个缺点是它对处理的数据不加选择,它可能会增加背景噪声的对比度并且降低有用信号的对比度。
实现
我们来看一个灰度图像,让表示灰度出现的次数,这样图像中灰度为 的像素的出现概率是
是图像中所有的灰度数, 是图像中所有的像素数, 实际上是图像的直方图,归一化到 。
把 作为对应于 的累计概率函数, 定义为:
是图像的累计归一化直方图。
我们创建一个形式为 的变化,对于原始图像中的每个值它就产生一个 ,这样 的累计概率函数就可以在所有值范围内进行线性化,转换公式定义为:
注意 T 将不同的等级映射到 域,为了将这些值映射回它们最初的域,需要在结果上应用下面的简单变换:
上面描述了灰度图像上使用直方图均衡化的方法,但是通过将这种方法分别用于图像RGB颜色值的红色、绿色和蓝色分量,从而也可以对彩色图像进行处理。
例子
下面的这个例子是一个8位的8×8灰度图像:
该灰度图像的灰度值出现次数如下表所示,为了简化表格,出现次数为0的值已经被省略
-
灰度值 出现次数 灰度值 出现次数 灰度值 出现次数 灰度值 出现次数 灰度值 出现次数 52 1 64 2 72 1 85 2 113 1 55 3 65 3 73 2 87 1 122 1 58 2 66 2 75 1 88 1 126 1 59 3 67 1 76 1 90 1 144 1 60 1 68 5 77 1 94 1 154 1 61 4 69 3 78 1 104 2 62 1 70 4 79 2 106 1 63 2 71 2 83 1 109 1
累积分布函数(cdf)如下所示,与上一表格类似,为了简化,累积分布函数值为0的灰度值已经被省略
-
灰度值 cdf 灰度值 cdf 灰度值 cdf 灰度值 cdf 灰度值 cdf 52 1 64 19 72 40 85 51 113 60 55 4 65 22 73 42 87 52 122 61 58 6 66 24 75 43 88 53 126 62 59 9 67 25 76 44 90 54 144 63 60 10 68 30 77 45 94 55 154 64 61 14 69 33 78 46 104 57 62 15 70 37 79 48 106 58 63 17 71 39 83 49 109 59
如表格所示,灰度值最小值为52,最大值为154. 通常,直方图均衡化算式如下:
累积分布函数最小值cdfmin在本例中为1, M和N分别代表了图像的长宽像素个数(本例中为8×8=64),而L则是灰度级数(如本例中,图像为8位深度,则灰度级数共有2^8=256级数,这也是最常见的灰度级数).则对于本例的直方图均衡化算式为:
例如,灰度为78的像素的累积分布函数为46,均衡化后,灰度值变化为:
直方图均衡化后,图像的灰度值变化如下表所示:
注意:原最小灰度值52变为了0,而最大灰度值154变为了255
彩色图像直方图均衡化 实际上,对彩色分量rgb分别做均衡化,会产生奇异的点,图像不和谐。一般采用的是用yuv空间进行亮度的均衡即可。