原创 图像的直方图均衡

1，MATLAB代码：

clear all
%一，图像的预处理，读入彩色图像将其灰度化
PS=imread('1.jpg');                 %读入JPG彩色图像文件
imshow(PS)                                  %显示出来
title('输入的彩色JPG图像')
imwrite(rgb2gray(PS),'PicSampleGray.bmp'); %将彩色图片灰度化并保存
PS=rgb2gray(PS);                            %灰度化后的数据存入数组

%二，绘制直方图
[m,n]=size(PS);                             %测量图像尺寸参数
GP=zeros(1,256);                            %预创建存放灰度出现概率的向量
for k=0:255
     GP(k+1)=length(find(PS==k))/(m*n);      %计算每级灰度出现的概率，将其存入GP中相应位置
end
figure,bar(0:255,GP,'g')                    %绘制直方图
title('原图像直方图')
xlabel('灰度值')
ylabel('出现概率')

%三，直方图均衡化
S1=zeros(1,256);
for i=1:256
     for j=1:i
          S1(i)=GP(j)+S1(i);                 %计算Sk
     end
end
S2=round((S1*256)+0.5);                          %将Sk归到相近级的灰度
for i=1:256
     GPeq(i)=sum(GP(find(S2==i)));           %计算现有每个灰度级出现的概率
end
figure,bar(0:255,GPeq,'b')                  %显示均衡化后的直方图
title('均衡化后的直方图')
xlabel('灰度值')
ylabel('出现概率')

  自己设计程序实现图像的直方图均衡；

自己设计的均值代码：

1. clear,close all;
2. I=imread('image.jpg');
3. subplot(2,2,1);
4. imshow(I);                                 
5. title('原始图像');
6. PS=rgb2gray(I);                           
7. [m,n]=size(I);                             
8. GP=zeros(1,256);                          
9. for k=0:255
10.      GP(k+1)=length(find(I==k))/(m*n);     
11. end
12. subplot(2,2,2)
13. bar(0:255,GP,'g')                   
14. title('原图像直方图')
15. xlabel('灰度值')
16. ylabel('出现概率')
17. N=zeros(1,256);
18. for i=1:256
19.      for j=1:i
20.           N(i)=GP(j)+N(i);               
21.      end
22. end
23. S2=round((N*256)+0.5);                         
24. for i=1:256
25.      GPeq(i)=sum(GP(find(S2==i)));          
26. end
27. subplot(2,2,3)
28. bar(0:255,GPeq,'b')                
29. title('均衡化后直方图')
30. xlabel('灰度值')
31. ylabel('出现概率')
32. X=PS;
33. for i=0:255
34.      X(find(PS==i))=S2(i+1);              
35. end
36. subplot(2,2,4)
37. imshow(X);                         
38. title('均衡化后图像');        画出均衡后的图像直方图；