冒泡....双十一刚过~购物狂欢完还是要收心学习鸭！今天来说一说直方图。

直方图

定义：

直方图是一种统计报告图，由一系列高度不等的纵向条纹或线段表示数据分布的情况。 一般用横轴表示数据类型，纵轴表示分布情况。
灰度直方图是灰度级的函数，描述的是图像中具有该灰度级的像元的个数。以横轴表示灰度级，以纵轴表示每一灰度级具有的像元数或该像元数占总像元数的比例值，做出的条形统计图即为灰度直方图。如以下：

直方图示例.png

由图像可以得到的一些直方图的性质：

a.图形直观可以发现它首先是一个一维离散函数图像。
b.直方图表示一幅图像中不同灰度像素出现的统计信息，它只能反映该图像中不同灰度值出现的频数（概率），而不能表示出像素的位置等其他信息。
c.任何一幅特定的图像都有唯一的直方图与之对应，但不同的图像可以有相同的直方图。

根据直方图的形态可以大致推断图像质量的好坏。由于图像包含有大量像元，其像元灰度值分布应符合概率统计分布规律。假定像元的灰度值是随机分布，那么其直方图应该是正态分布。

image.png

直方图的均衡化

image.png

image.png

接下来就推导一下均衡化的过程：
首先考虑连续的灰度值，用变量 r 表示输入图像的灰度，用 s 表示输出图像的灰度。r 范围是 [ 0, L-1]，则函数： s = T(r), r 的范围：[0 , L-1] （0为全黑，L-1为全白）。
T(r)满足下列两个条件：
1.T(r)在区间0≤r≤1中为 严格单调递增函数，保证原图各灰度级在变换后仍保持从黑 到白（或从白到黑）的排列次序）
2.当0≤r≤1时,0≤T(r) ≤1，保证变换前后灰度值动态范围的一致性。
在处理的过程中需要引入重要的变换函数：

image.png

image.png

image.png

最后可以发现：得到ps(s)是一个均匀的概率密度函数，与pr(r)无关
对于离散值，我们处理其概率（直方图值）与求和来替代处理概率密度函数与积分。因此，一幅数字图像中灰度级出现的概率（近似）和离散变换形式如下：

iMN为图像的总像素

image.png

image.png

注意：sk=6的为s3和s4故nk是两个的叠加，p(sk)的值为nk/MN(4096)

image.png

Matlab代码参考

I = imread('rice.png');
[height,width] = size(I);
figure
subplot(221)
imshow(I)%显示原始图像
subplot(222)
imhist(I)%显示原始图像直方图
 
%进行像素灰度统计;
NumPixel = zeros(1,256);%统计各灰度数目，共256个灰度级
for i = 1:height
    for j = 1: width
        NumPixel(I(i,j) + 1) = NumPixel(I(i,j) + 1) + 1;%对应灰度值像素点数量增加一
    end
end
%计算灰度分布密度
ProbPixel = zeros(1,256);
for i = 1:256
    ProbPixel(i) = NumPixel(i) / (height * width * 1.0);
end
%计算累计直方图分布
CumuPixel = zeros(1,256);
for i = 1:256
    if i == 1
        CumuPixel(i) = ProbPixel(i);
    else
        CumuPixel(i) = CumuPixel(i - 1) + ProbPixel(i);
    end
end
%累计分布取整
CumuPixel = uint8(255 .* CumuPixel + 0.5);
%对灰度值进行映射（均衡化）
for i = 1:height
    for j = 1: width
        I(i,j) = CumuPixel(I(i,j));
    end
end
 
subplot(223)
imshow(I)%显示原始图像
subplot(224)
imhist(I)%显示原始图像直方图

直方图的规定化

从文字就可以看出，规定化就是把直方图转化为既定的规定化的形状。
过程如下：
先设 Pr（r）和Pz（z）分别表示原始灰度图像和目标图像的灰度分布概率密度函数.根据直方图规 定化的特点与要求，应使原始图像的直方图具有Pz（z）所表示的形状。因此建立Pr（r）和 Pz（z）之间的关系是直方图规定化必须解决的问题。

image.png

第一步：要先根据图a提供的直方图求出输出之后均衡的值

s3 = 5.67≈6 ；s4 =6.23≈6 ；s5 = 6.66≈7；
s6= 6.86≈7 ；s7 = 7.00≈7 ；
第二步：求变换函数G的所有值

image.png

image.png

image.png

image.png

image.png

参考资料：（https://blog.csdn.net/u013355826/article/details/57417771）
（https://blog.csdn.net/taoyanqi8932/article/details/52758376）
（https://blog.csdn.net/xiajun07061225/article/details/6910129）

Ending 不容易地收尾啦~