matlab 图像直方图处理,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('规定化函数');

效果图：