clear
clc
image = imread('hw4.jpg');
[m,n] = size(image);
figure
subplot(221)
imshow(image)%显示原始图像
subplot(222)
imhist(image)%显示原始图像直方图
g1 = zeros(1,256);
g2 = g1;
for i = 1:m
for j = 1: n
g1(image(i,j) + 1) = g1(image(i,j) + 1) + 1;%对应灰度值像素点数量增加一
end
end
g1 = g1 / (m * n * 1.0);%得到灰度级所占比例
for i = 1:256
for j = 1:i
g2(i) = g2(i)+g1(j);
end
end
for i = 1:256 %直方图分布
g2(i) = floor(255*g2(i)+0.5);
end
for i = 1:m %对灰度值进行映射
for j = 1: n
image(i,j) = g2(image(i,j)+1);
end
end
%%%%%%%%%%%%%%结果比对%%%%%%%%
subplot(223)
imshow(image)%显示原始图像
subplot(224)
imhist(image)%显示原始图像直方图
评论0