灰度直方图均衡化及其实现

直方图均衡化

定义

当直方图中像素值集中在狭窄的灰度级范围内或分布极不均匀时，图像呈现较差的对比度。直方图均衡化的目的就是将直方图的灰度级概率分布变换为均匀分布。(如下图所。注意，p(x)是概率密度函数，P(x)是概率分布函数）

由于直方图中概率较小的灰度级合并为更少的几个或一个灰度级内，从而降低了图像的灰度级分辨率，且某些细节信息处于概率较小的灰度级中，这样的灰度级归并到其他灰度级内，从而造成图像细节信息的丢失。

公式

输入r，输出s
$$s=T(r)=(L-1)\sum _{j=0}^{k}P(r_j)=\frac{L-1}{MN}\sum _{j=0}^k{n}_j,(k=0,1...L-1)$$

对于输入图像中每个具有r值的像素值产生一个输出灰度值s.

略微证明

$$先定义两个函数s=T(r),r=T^{-1}(s)。两个互为反函数$$
接下来，F（x）是概率分布函数，P（x）概率分布。
$$F_s(s)=P_r(S<s)=P_r(T(R)<s)=P_r(R<T^{-1}(s))=P_r(R<r)=F_r(r)$$

接下来对上式两边求导s（概率分布函数的导数是概率密度函数），可得如下
$$p_s(s)=p_r(r)\frac{dr}{ds}$$
当直方图均衡化的变换函数T( r )是连续型随机变量R的概率分布函数Fr( r),也就是输入灰度级r的累计直方图时，可表示为
$$s=F_r(r)={\int }_0^r{p}_r(x)dx,$$
对上面两边求导r，得到
$$\frac{ds}{dr}=p_r(r)$$

将这个式子带入当初求的哪个，得到
$$p_s(s)=1$$
这说明什么？？

当直方图均衡化的变换函数T( r)是连续型随机变量R的概率分布函数Fr( r),也就是输入灰度级r的累计直方图时，变换后的概率密度函数是均匀分布的！！！

所以
$$s=T(r)=(L-1)\sum _{j=0}^{k}P(r_j)=\frac{L-1}{MN}\sum _{j=0}^k{n}_j,(k=0,1...L-1)$$

计算步骤

1. 计算$$P(r_k)=\frac{n_k}{MN}(k=0,1，2,3...L-1)$$
   其中nk是出现的次数。L是离散灰度级数
2. 计算$$s_k=(L-1)\sum _{j=0}^{k}P(r_j),(K=0,1,...L-1)四舍五入为整数$$
3. 确定$$r_k与s_k的对应关系$$
4. 确定概率分布
   

matlab实现

i=imread('d:\\cat.jpg');
i=rgb2gray(i);//rgb转化为灰度图
j=histeq(i);//灰度直方图均衡化
figure('name','contrast');
subplot(1,2,1);imshow(i);title('raw pic');
subplot(1,2,2);imshow(j);title('process pic');