Matlab数字图像处理-直方图变换

直方图均衡化

clear;
close all;
img = imread('E:\DIP\image_gray.jpg');
a = size(img);
%得到直方图
h = imhist(img);
%得到矩阵最大值
L =length(h);
figure(1);
subplot(321);
imshow(img);
title('原图');
subplot(322);
bar(h);
title('原直方图');

%得到直方图分布函数
fenbu = h/sum(h);
leiji = [];
for i = 1:1:L
    if(i==1)
        leiji(i) = fenbu(i);
    else
        leiji(i) = leiji(i-1) + fenbu(i);
    end
end
%映射
huidu = round(leiji.*(L-1));
junheng = zeros(a(1),a(2));
h2 = [];
for i = 1:1:a(1)*a(2)
    junheng(i) = huidu(img(i)+1)/L;
    h2(i) = junheng(i);
end
subplot(323);
imshow(junheng);
title('均衡化');
%得到直方图
h2 = hist(h2,L);
subplot(324);
bar(h2);
title('均衡化直方图');

效果图：

直方图规定化：（这个不对）

clear;
close all;
img = imread('E:\DIP\image_gray.jpg');
a = size(img);
%得到直方图
h = imhist(img);
%得到矩阵最大值
L =length(h);
figure(1);
subplot(221);
imshow(img);
title('原图');
subplot(222);
bar(h);
title('原直方图');

[f T] = histeq(img,256);
subplot(223);
imshow(f);
title('规定化后图像');

x = linspace(0,1,256);
subplot(224);
plot(x,T)
axis([0 1 0 1])
title('规定化函数');

效果图：