数字图像处理实验报告_打印..doc

实验一I = imread(1.jpg); I1=im2double(I); figure,subplot(151),imshow(I); title(Original image); I2=I1.^2; subplot(152), imshow(I2); title(a=2 b=0); imwrite(I2,1-1-1.jpg); I3=I1.^0.5; subplot(153), imshow(I3); title(a=0.5 b=0); imwrite(I3,1-1-2.jpg); I4=I1+0.3; subplot(154), imshow(I4); title( a=1 b=0.3); imwrite(I4,1-1-3.jpg); I5=I1-0.3; subplot(155), imshow(I5); title( a=1 b=-0.3); imwrite(I5,1-1-4.jpg);2.图象二值化 I = imread(1.jpg); I1=im2double(I); subplot(141),imshow(I1); title(original image); I2 = (I10.2); subplot(142),imshow(I2); title(Level=0.2); imwrite(I2,1-2-1.jpg); I3 = (I10.6); subplot(143),imshow(I3); title(Level=0.6); imwrite(I3,1-2-2.jpg); I4 = (I10.8); subplot(144),imshow(I4); title(Level=0.8); imwrite(I4,1-2-3.jpg); 3.负向变换 I = imread(1.jpg); I1=im2double(I); subplot(121),imshow(I1); title(Original image); fz=1-I1;subplot(122),imshow(fz); title(Reserved image); imwrite(fz,1-3-Reserved.jpg); 二. 图象直方图及均衡 I=imread(2.jpg); add0=[]; add1=[]; tab1=zeros(1,256); tab2=zeros(1,256); I1=I; for n=0:255 X=I==n; add0=[add0;sum(sum(X))]; end; [a,b]=size(I); final=add0/(a*b); for n=1:256 for i=1:n tab1(n)=tab1(n)+final(i); end; end; tab1=tab1*255; tab2=uint8(tab1); for n=1:a for m=1:b for t=0:255 if I(n,m)==t I1(n,m)=tab2(t+1); end; end; end; end; for n=0:255 X1=I1==n; add1=[add1;sum(sum(X1))]; end; [a1,b1]=size(I1); final1=add1/(a1*b1); imwrite(I1,2equaled.jpg); subplot(221),imshow(I); title(原始图像) subplot(223),imshow(I1);title(均衡化的图片) subplot(222),bar(0:255,final);title(原图像直方图) subplot(224),bar(0:255,final1);title(均衡化的直方图)  三：图像中值滤波 I=imread(3.jpg); [m,n]=size(I); J=zeros(m+2,n+2); for i=1:(m+2) if (j==1||j==n+2) J(i,j)=0; end end for j=1:(n+2) if (i==1||i==m+2) J(i,j)=0; end end for i=2:(m+1) for j=2:(n+1) J(i,j)=I(i-1,j-1); end end % 3×3中值滤波 for i=2:m+1 for j=2:n+1 M=[J(i-1,j-1), J(i-1,j), J(i-1,j