图像实验目录及源程序.docVIP

图像实验目录 图像增强: 1.对比度增强(dip1.m) 2. 对数变换(dip2.m) 3.直方图均衡(dip3.m) 4.位平面分解(dip4.m) 5.均值滤波和中值滤波(dip5.m) 6.锐化滤波(dip6.m) 源程序 dip1.m X1=imread(pout.tif); figure,imshow(X1); f0=0;g0=0; f1=70;g1=30; f2=180;g2=230; f3=255;g3=255; figure,plot([f0,f1,f2,f3],[g0,g1,g2,g3]) axis tight,xlabel(f),ylabel(g); title(intensity transformation); r1=(g1-g0)/(f1-f0); b1=g0-r1*f0; r2=(g2-g1)/(f2-f1); b2=g1-r2*f1; r3=(g3-g2)/(f3-f2); b3=g2-r3*f2; [m,n]=size(X1); X2=double(X1); for i=1:m for j=1:n f=X2(i,j); g(i,j)=0; if(f=0)(f=f1) g(i,j)=r1*f+b1; elseif(f=f1)(f=f2) g(i,j)=r2*f+b2; elseif(f=f2)(f=f3) g(i,j)=r3*f+b3; end end end figure,imshow(mat2gray(g)) dip2.m I=imread(pout.tif); subplot(1,2,1); imshow(I); I=double(I); J=log(I+1); subplot(1,2,2); imshow(mat2gray(J)); dip3.m I=imread(pout.tif); [J,T]=histeq(I); imshow(I); figure,imshow(J); figure,imhist(I) figure,imhist(J,64) figure,plot((0:255)/255,T); dip4.m I=imread(saturn.tif); figure,imshow(I); I=double(I); [M,N]=size(I); for k=1:8 J=zeros(M,N); for i=1:M for j=1:N temp=I(i,j); s1=0;s2=0; range=[k:-1:1]; for d=range s1=2^(8-d)+s1; s2=2^(8-d+1); if temp=s1temps2; J(i,j)=255;break; end end end end J=uint8(J); figure;imshow(J); end figure; dip5.m I=imread(saturn.tif); J=imnoise(I,salt pepper,0.05); imshow(I); figure,imshow(J); J1=filter2(fspecial(average,3),J)/255; % J2=filter2(fspecial(average,5),J)/255; % J3=filter2(fspecial(average,7),J)/255; figure,imshow(J1); % figure,imshow(J2); % figure,imshow(J3); K1=medfilt2(J,[3,3]); figure,imshow(K1); dip6.m % li hua I=imread(saturn.tif); figure,imshow(I); I1=filter2(fspecial(gaussian),I); I2=filter2(fspecial(sobel),I); I3=filter2(fspecial(prewitt),I); I4=filter2(fspecial(log),I); I5=filter2(fspecial(laplacian),I); figure,