matlab实现灰度直方图均匀化

p=imread('test.jpg');
figure(1);
imshow(p);
p=rgb2gray(p); %转换为灰度图像
figure(11);
imshow(p);
figure(2);
imhist(p);  %显示灰度分布图
figure(1234);
histgram=zeros(1,256); %灰度统计数组
[m,n]=size(p);
for i=1:m
    for j=1:n
        k=p(i,j);
        histgram(k+1)=histgram(k+1)+1;
    end
end
x=1:255;
plot(x,histgram(x));
%以上为灰度的分布情况


%计算每一个点的概率值
all=m*n;
histgram(1)=histgram(1)/all;
for i=2:255
    histgram(i)=histgram(i)/all;
    histgram(i)=histgram(i-1)+histgram(i);
end
figure(3);
x=1:255;
plot(x,histgram(x));
%开始计算新的灰度值
for i=1:m
    for j=1:n
        k=p(i,j);
        p(i,j)=histgram(k+1)*256;
    end
end
%显示直方均匀化后的图像
figure(4);
imshow(p);

imhist(p);