大连理工数字图像处理作业Word格式.docx

1、stem（0:N-1）,y2,.）;C=fft（y2）;C1=fft（a）.*fft（b）;all（C（:）=C1（:）;圆周卷积函数：function y=yuanzhou（x1,x2,N）x1=x1,zeros（1,N-length（x1）;x2=x2,zeros（1,N-length（x2）;m=0:1:N-1;H=zeros（N,N）;for n=1:N H（n,:）=cirshift（x2,n-1,N）;endy=x1*H;结果： 4-9t1=imread（Lena.tift=rgb2gray（t1）;imshow（t）,title（原图像a,b=size（t）;x,y=meshgr

2、id（-a/2:a/2-1,-b/2:b/2-1）;z=sqrt（x.2+y.2）;c=（ztf2=tf.*c1;ty3=ifft2（tf2）;imshow（abs（ty3）,）,title（高通滤波之后图像g1=mat2gray（fspecial（gaussian,512,10）;tf3=tf.*g1;ty4=ifft2（tf3）;imshow（abs（ty4）,）,title（高斯低通滤波器滤波之后图像bl=lbutter（c,15,1）;tfbl=tf.*bl;ty5=ifft2（tfbl）;imshow（abs（ty5）,）,title（巴特沃兹低通滤波之后图像巴特沃兹滤波器函数设计：

3、function out=lbutter（im,d,n）height=size（im,1）;width=size（im,2）;x,y=meshgrid（-floor（width/2）:floor（width-1）/2）,-floor（height/2）: .floor（height-1）/2）;out=1./（1+（sqrt（2）-1）*（x.2+y.2）/d2）.n）;结果5-6f=imread（fg=rgb2gray（f）;tn=imnoise（fg,salt & pepper,0.05）;imshow（tn,）,title（a1=fspecial（averagetna1=filter2（

4、a1,tn）;imshow（tna1,）,title（均值滤波之后图像tnz=medfilt2（tn,3,3）;imshow（tnz）,title（中值滤波之后图像tnout=outlier（tn,0.2）;imshow（tnout,）,title（异常值判断滤波之后D=0.2tnord=ordfilt2（tn,3,0 1 0; 1 1 1;0 1 0）;imshow（tnord,）,title（伪中值滤波之后图像伪中值滤波函数：function res=outlier（im,d）f=0.125 0.125 0.125; 0.125 0 0.125; 0.125 0.125 0.125;imd

5、=im2double（im）;imf=filter2（f,imd）;r=abs（imd-imf）-d0;res=im2uint8（r.*imf+（1-r）.*imd6-2t=imread（1.tifx,y=meshgrid（1:a,1:b）;s1=1+sin（x/3+y/5）;s2=1+sin（x/5+y/1.5）;s3=1+sin（x/6+y/6）;ts1=（im2double（t）+s1/2）/2;ts2=（im2double（t）+s2/2）/2;ts3=（im2double（t）+s3/2）/2;tps1=fftshift（fft2（ts1）;imshow（log（abs（tps1）,

6、）,title（a组加噪声之后频域图）,impixelinfo;imshow（abs（ifft2（tps1）,）,title（）tps1（:,114:116）=0;%（中心位置115,120；143,138）,142:144）=0;tps1（119:121,:）=0;tps1（137:139,:tps1_1=ifft2（tps1）;imshow（abs（tps1_1）,）,title（a组噪声在陷波滤波器滤波之后结果tps2=fftshift（fft2（ts2）;imshow（log（abs（tps2）,）,title（b组加噪声之后频域图tps2（:,119:121）=0;%（中心位置120

7、,102；137,156）,136:138）=0;tps2（101:103,:tps2（155:157,:tps1_2=ifft2（tps2）;imshow（abs（tps1_2）,）,title（b组噪声在陷波滤波器滤波之后结果tps3=fftshift（fft2（ts3）;imshow（log（abs（tps3）,）,title（c组加噪声之后频域图tps3（:,121:123）=0;%（中心位置122,122；136,136,135:137）=0;tps3（121:123,:tps3（135:137,:tps1_3=ifft2（tps3）;imshow（abs（tps1_3）,）,tit

8、le（c组噪声在陷波滤波器滤波之后结果7-7circles.tif256,1:256）;t2=double（t）.*（x+y）/2+64）+x+y;t3=uint8（255*mat2gray（t2）;fun=inline（im2bw（x,graythresh（x）a,b=size（t3）;img1=blkproc（t3,a/2,b/2,fun）;subplot（221）;imshow（img1,）,title（分块矩阵大小a/2,b/2img2=blkproc（t3,a/4,b/4,fun）;subplot（222）;imshow（img2,）,title（分块矩阵大小a/4,b/4img3=

9、blkproc（t3,a/8,b/8,fun）;subplot（223）;imshow（img3,）,title（分块矩阵大小a/8,b/8img4=blkproc（t3,a/16,b/16,fun）;subplot（224）;imshow（img4,）,title（分块矩阵大小a/16,b/168-6t=rgb2gray（t）;t1=imnoise（t,0.1）;t2=imnoise（t,0,0.02）;tu1=edge（t1,sobeltu2=edge（t1,robertstu3=edge（t1,prewitttu4=edge（t1,logtu5=edge（t1,canny,0.01,0.

10、5）;%subplot（211）;imshow（im2bw（tu1,graythresh（double（tu1）,title（加椒盐噪声后sobel函数边缘处理之后%subplot（212）;%imshow（tu1）;imshow（im2bw（tu2,graythresh（double（tu2）,title（加椒盐噪声后roberts函数边缘处理之后imshow（im2bw（tu3,graythresh（double（tu3）,title（加椒盐噪声后prewitt函数边缘处理之后imshow（im2bw（tu4,graythresh（double（tu4）,title（加椒盐噪声后log函数

11、边缘处理之后imshow（im2bw（tu5,graythresh（double（tu5）,title（加椒盐噪声后canny函数边缘处理之后tu21=edge（t2,tu22=edge（t2,tu23=edge（t2,tu24=edge（t2,tu25=edge（t2,imshow（im2bw（tu21,graythresh（double（tu21）,title（加高斯噪声后sobel函数边缘处理之后imshow（im2bw（tu22,graythresh（double（tu22）,title（加高斯噪声后roberts函数边缘处理之后imshow（im2bw（tu23,graythresh（double（tu23）,title（加高斯噪声后prewitt函数边缘处理之后imshow（im2bw（tu24,gr