[数字图像处理]灰度直方图均衡化

function [  ] = histChange( A )
%histChange 此处显示有关此函数的摘要
%对输入图像矩阵进行灰度直方图均衡化,若输入为RGB图像矩阵，则自动转换为灰度图像进行处理
%
[M,N,a]=size(A);
if a == 3
    B=rgb2gray(A);
else
    B=A;
end
x=0:255;
y=zeros(1,256);
s=zeros(1,256);
subplot(2,2,1);
total=0;
for k=1:M
    for p=1:N
        y(B(k,p)+1)=y(B(k,p)+1)+1;
        total=total+1;
    end
end

bar(x,y);
y=y/total;
for k=1:256
    for m=1:k
        s(k)=s(k)+y(m);
    end
end
s=s*255;
subplot(2,2,2);
imshow(B);
for k=1:M
    for p=1:N
        B((k),(p))=s(B(k,p)+1);
    end
end
subplot(2,2,4);
imshow(B);
u=zeros(1,256);
for k=1:M
    for p=1:N
        u(B(k,p)+1)=u(B(k,p)+1)+1;
    end
end 
subplot(2,2,3);
bar(x,u);
end

转载请注明来源