clear;
I=imread('F:\\experiment\\emeishan.jpg');
I=I(:,:,1); % 像素值取RGB中的一层 [m,n]=size(I); %获得原图像I的行列数。 pre_mat=zeros(1,256); %为变换前后的直方图分别申请矩阵。
aft_mat=zeros(1,256); for i=1:m for j=1:n
pre_mat(I(i,j)+1)=pre_mat(I(i,j)+1)+1; %统计各灰度级的像素总数。 end end
aft_mat(1)=aft_mat(1)/(m*n); for k=2:256
s(k)=pre_mat(k)/(m*n); %计算出灰度的概率。 aft_mat(k)=aft_mat(k-1)+s(k); %累积直方图。 end
M=zeros(m,n); %申请m*n阶矩阵用于装入处理后的各灰度级概率。 aft_mat=aft_mat*255; for i=1:m for j=1:n
M(i,j)=aft_mat(I(i,j)+1); %将处理后的灰度装入新的矩阵。 end end
J=uint8(M); %转成unit8型。 figure(1); %显示原图像和直方图 imshow(I); xlabel('原始图像'); figure(2); imhist(I);
xlabel('原始图像直方图')
figure(3); %显示均衡化后的图像和直方图 imshow(J);
xlabel('均衡化图像'); figure(4); imhist(J);
xlabel('均衡化后直方图');
