数字图像处理直方图均衡化MATLAB源程序

%zhifangtuqunhenghua clear x1=imread('1.tif'); [m,n]=size(x1); x2=double(x1); lenna=zeros(m,n); lenna_equ=zeros(n,n); histgram=zeros(256); cdf=zeros(256); [lenna,map]=imread('1.tif'); %get histogram for i=1:m for j=1:n k=lenna(i,j) histgram(k)=histgram(k)+1 end end %get cdf cdf(1)=histgram(1) for i=2:256 cdf(i)=cdf(i-1)+histgram(i); end %run point operation for i=1:m for j =1:n k=lenna(i,j); lenna_equ(i,j)=cdf(k)*256/(m*n); end end %生成直方图均衡化后的lenna图 g=mat2gray(lenna_equ); figure subplot(2,2,1);imshow(lenna) subplot(2,2,2);imhist(lenna) ylim('auto'); subplot(2,2,3);imshow(g) subplot(2,2,4);imhist(g) ylim('auto') f=imread('1.tif'); figture subplot(2,2,1);imshow(f); title('the orignal image'); subplot(2,2,2);imhist(f); ylim('auto') g=histeq(f,256); subplot(2,2,3);imshow(g); title('image after equalization') subplot(2,2,4);imhist(g); ylim('auto'); %zhifangtupipeichengxu I = imread('pout.tif'); Ieq=histeq(I); [c,x]=imhist(Ieq); I1=imread('tire.tif'); J=histeq(I1,c); subplot(2,2,1); imshow(I); title('Image A'); subplot(2,2,2); imshow(Ieq); title('Image A after Equalization'); subplot(2,2,3); imshow(I1); title('Image B'); subplot(2,2,4); imshow(J); title('Image B after Matching');