matlab绘制直方图均衡化,matlab直方图均衡化

%二，绘制直方图

[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('出现概率')

%四，图像均衡化 PA=PS; for i=0:255      PA(find(PS==i))=S2(i+1);                %将各个像素归一化后的灰度值赋给这个像素 end figure,imshow(PA)                           %显示均衡化后的图像 title('均衡化后图像') imwrite(PA,'PicEqual.bmp');