untitled.m

图象处理中的直方图程序,有些文件名需自己改动一下

I=imread('rs.jpg'); 
[m,n]=size(I);%获得该图的行列数。 
m1=zeros(1,256);%为直方图申请矩阵。 
m2=zeros(1,256); 
for i=1:m 
 for j=1:n 
 m1(I(i,j)+1)=m1(I(i,j)+1)+1;%统计各级灰度的像素数。 
end 
end 
m2(1)=0; 
for k=2:256 
 m3(k)=m1(k)/(m*n);%计算出灰度的概率。 
 m2(k)=m2(k-1)+m3(k);%累加直方图。 
end 
M=zeros(m,n); 
for i=1:m 
 for j=1:n 
 M(i,j)=(255*m2(I(i,j))); %将处理后的灰度装入新的矩阵。 
 end 
end 
P=uint8(M);%转成unit8型。 
subplot(3,2,1);imshow(I);title('原始图像');subplot(3,2,2);imhist(I);title('原始直方图');%输出图像和直方图，加以对比。