数字图像处理与matlab编程.doc

数字图像处理与matlab编程 1、图像的二值化 clear i=imread(lenna_2.bmp); j=double(i); [m n]=size(j); lk=j; for p=1:m; for q=1:n; if j(p,q)=125 lk(p,q)=255; else lk(p,q)=0; end end end k=uint8(lk); figure,imshow(i); figure,imshow(k) 实验结果： 处理前 处理后 2、数字图像旋转、放大及缩小 1）旋转 i=imread(Lenna.bmp); j=imrotate(i,45,bilinear); subplot(1,2,1),imshow(i); subplot(1,2,2),imshow(j); 实验结果： 2）放大 i=imread(Lenna.bmp); j=imresize(i,1.5,bilinear); figure(1),imshow(i); figure(2),imshow(j); 实验结果： 放大前 放大后 3）缩小 i=imread(Lenna.bmp); j=imresize(i,0.8,bilinear); figure(1),imshow(i); figure(2),imshow(j); 实验结果： 缩小前 缩小后 3、数字图像水平镜像、垂直镜像和平移 1）数字图像的水平镜像 i=imread(Lenna.bmp); subplot(1,2,1),imshow(i); [x,y]=size(i); p=1:x;q=1:y; j(p,y-q+1)=i(p,q); subplot(1,2,2),imshow(j); 实验结果： 2）数字图像的垂直镜像 i=imread(Lenna.bmp); subplot(1,2,1),imshow(i); [x,y]=size(i); p=1:x;q=1:y; j(x-p+1,q)=i(p,q); subplot(1,2,2),imshow(j); 实验结果： 3）数字图像的平移 i=imread(Lenna.bmp); subplot(1,2,1),imshow(i); [x,y]=size(i); p=1:x;q=1:y; j(p+20,q+20)=i(p,q); j(1:p+19,q+20)=0; j(p+20,p+20:q+19)=0; subplot(1,2,2),imshow(j); 实验结果： 4、图像的错切 clear a=imread(lenna_2.bmp); figure,imshow(a); [m n]=size(a); b(1:m+0.5*n,1:n)=0; for x=1:m for y=1:n b(fix(x+0.5*y),y)=double(a(x,y)); end end figure,imshow(uint8(b)); 实验结果： 错切前 错切后 5、中值滤波 i=imread(rice.tif); j=imnoise(i,‘salt pepper ’,0.02);加椒盐噪声 i2=medfilt2(j,[3,3]); subplot(1,3,1),imshow(i);title(原图像) subplot(1,3,2),imshow(j);title(加噪声图像) subplot(1,3,3),imshow(i2,[]);title(中值滤波后图像) 实验结果： 6、均值滤波 i=imread(flowers.tif); j=imnoise(i,‘gaussian’,0.02);加高斯噪声 h=fspecial(average,3); i2=filter2(h,j); subplot(1,3,1),imshow(i);title(原图像) subplot(1,3,2),imshow(j);title(加噪声图像) subplot(1,3,3),imshow(i2,[]);title(均值滤波后图像) 实验结果： 7、图像加波浪噪声 i=imread(rice.tif); s=size(i); [x,y]=meshgrid(1:s(1),1:s(