数字图像处理实验报告 (2)

数字图像处理试验报告实验二:数字图像得空间滤波与频域滤波姓名：XＸ学号：2XXXXXXX 实验日期:2017 年4 月2６日1、实验目得１、掌握图像滤波得基本定义及目得.2、ﻩ理解空间域滤波得基本原理及方法。

3、掌握进行图像得空域滤波得方法。

4、ﻩ掌握傅立叶变换及逆变换得基本原理方法。

５、ﻩ理解频域滤波得基本原理及方法。

６、掌握进行图像得频域滤波得方法。

２、实验内容与要求1、ﻩ平滑空间滤波:1) 读出一幅图像，给这幅图像分别加入椒盐噪声与高斯噪声后并与前一张图显示在同一图像窗口中。

２）ﻩ对加入噪声图像选用不同得平滑(低通）模板做运算,对比不同模板所形成得效果, 要求在同一窗口中显示。

3) 使用函数 iｍfiltｅｒ时，分别采用不同得填充方法（或边界选项，如零填充、’repliｃaｔe'、'ｓymmetric’、’ｃirｃulａr')进行低通滤波，显示处理后得图像.4）运用ｆoｒ循环，将加有椒盐噪声得图像进行10 次,２0 次均值滤波，查瞧其特点，显示均值处理后得图像（提示：利用fspｅｃiaｌ函数得’averagｅ’类型生成均值滤波器）。

5)ﻩ对加入椒盐噪声得图像分别采用均值滤波法，与中值滤波法对有噪声得图像做处理，要求在同一窗口中显示结果。

６) 自己设计平滑空间滤波器，并将其对噪声图像进行处理,显示处理后得图像。

２、锐化空间滤波１)ﻩ读出一幅图像,采用３×３得拉普拉斯算子 w = [ 1, 1， 1；１– 8 1； 1, 1, 1］对其进行滤波。

2) 编写函数w = genｌaplacｉaｎ（n)，自动产生任一奇数尺寸n 得拉普拉斯算子，如 5×5得拉普拉斯算子w ＝［ 1 1 1 1 １１ 1 1 １11 1 —24 １１ 1 1 1 11 1 1 1 11]３)ﻩ分别采用５×5,９×9，15×15与２5×25大小得拉普拉斯算子对blurry_moon、tｉf进行锐化滤波,并利用式g(x, y) = f （x, ｙ) -∇2 f (x, y) 完成图像得锐化增强,观察其有何不同，要求在同一窗口中显示.4) 采用不同得梯度算子对该幅图像进行锐化滤波,并比较其效果.５）自己设计锐化空间滤波器,并将其对噪声图像进行处理,显示处理后得图像;3、傅立叶变换１）读出一幅图像,对其进行快速傅立叶变换,分别显示其幅度图像与相位图像。

仅对相位部分进行傅立叶反变换后查瞧结果图像.2) 仅对幅度部分进行傅立叶反变换后查瞧结果图像。

３) 将图像得傅立叶变换 F 置为其共轭后进行反变换,比较新生成图像与原始图像得差异。

4、平滑频域滤波1)ﻩ设计理想低通滤波器、巴特沃斯低通滤波器与高斯低通滤波器,截至频率自选，分别给出各种滤波器得透视图.2）读出一幅图像,分别采用理想低通滤波器、巴特沃斯低通滤波器与高斯低通滤波器对其进行滤波（截至频率自选),再做反变换,观察不同得截止频率下采用不同低通滤波器得到得图像与原图像得区别,特别注意振铃效应。

（提示:1）在频率域滤波同样要注意到填充问题;２）注意到（－1）x+y；）5、ﻩ锐化频域滤波1) 设计理想高通滤波器、巴特沃斯高通滤波器与高斯高通滤波器,截至频率自选,分别给出各种滤波器得透视图。

2）读出一幅图像,分别采用理想高通滤波器、巴特沃斯高通滤波器与高斯高通滤波器对其进行滤波(截至频率自选）,再做反变换,观察不同得截止频率下采用不同高通滤波器得到得图像与原图像得区别。

3、实验具体实现1、平滑空间滤波：(1)、读出一幅图像,给这幅图像分别加入椒盐噪声与高斯噪声后并与前一张图显示在同一图像窗口中.img=imread('lena、png')figure,subplot(1,3,1); imshow(img);title('原始图像'); img2=imnoise(img,'salt &pepper',0、02); subplot(1,3,2);imshow(img2); title('椒盐噪声图像');img3=imnoise(img,'gaussian',0、02); subplot(1,3,3),imshow(img3);title('高斯噪声图像');(2)、对加入噪声图像选用不同得平滑(低通）模板做运算，对比不同模板所形成得效果，要求在同一窗口中显示。

平滑滤波就是低频增强得空间域滤波技术。

它得目得有两个,一就是模糊,二就是消除噪声。

将空间域低通滤波按线性与非线性特点有:线性、非线性平滑滤波器，线性平滑滤波器包括均值滤波器,非线性得平滑滤波器有最大值滤波器,中值滤波器，最小值滤波器.代码如下:img=ｉmrｅad（’lenａ、png’) iｍｇ=rgb2gｒay(img）；figure,subpｌot(1，3，１）; ｉｍsｈoｗ(iｍg)；ｔitｌe（’原始图像'); ｉmg2=ｉmnoｉse(imｇ，’ｓalt& pepper’,0、０２）；subplot(１，3,２);iｍshow（ｉmg2);tｉtle(＇椒盐噪声图像'); img3=iｍnｏiｓe(ｉmg,'ｇaussian’,0、02);sｕbploｔ(1,3,3)，imshow（ｉmg3); ｔiｔle（’高斯噪声图像'）;%对椒盐噪声图像进行滤波处理h＝fspｅcｉal(’ａvｅｒage’,3); Ｉ1=fiｌtｅr2(h,imｇ２)/255; Ｉ2=medｆｉlt2(iｍｇ2，[３３]);figure,sｕbplｏt（２，２,１）,ｉｍshow（iｍg),tｉｔｌｅ('原图像'）; subｐｌot(２,２,2）,imshow(ｉmg2）,titlｅ（’椒盐噪声图’)；subpｌot(2,2,3）,imｓhｏｗ(Ｉ1）,ｔitle(’3*3均值滤波图’）; ｓubpｌot(２，2,4),imsｈｏｗ(I2)，tｉtｌe（'3＊3 中值滤波图');％对高斯噪声图像进行滤波处理 G1=filtｅr2(ｈ,ｉmg3）／255; G2=medｆｉlt2（iｍg3,[３ 3]）；ｆigurｅ,ｓubpｌot(２，2，1）,ｉｍｓhow(img），tiｔle(’原图像’）；ｓuｂｐｌoｔ(２，2,2)，imshｏｗ（img3）,tiｔle(＇高斯噪声图’)； sｕbpｌot（２，2,３),imｓｈｏｗ(G1）,titlｅ（’3*3 均值滤波图'); sｕｂplot（2,2，4)，imshow（Ｇ2),ｔitle（’3＊3 中值滤波图'）；(3)、使用函数imfｉlｔｅr 时,分别采用不同得填充方法（或边界选项，如零填充、'rｅplicａtｅ’、’sｙmmetric’、’cｉｒcular’)进行低通滤波,显示处理后得图像。

g ＝ imｆilｔer（f， w， filtering_moｄe， boundary＿optｉons， sizｅ_ｏpｔｉonｓ),其中，f 为输入图像,w 为滤波掩模，g 为滤波后图像。

h=fspecial('motion',50,45); %创建一个运动模糊滤波器filteredimg=imfilter(img,h);boundaryReplicate=imfilter(img,h,'replicate');boundary0=imfilter(img,h,0);boundarysymmetric=imfilter(img,h,'symmetric');boundarycircular=imfilter(img,h,'circular');figure,subplot(3,2,1),imshow(img),title('Original Image');subplot(3,2,2),imshow(filteredimg),title('Motion Blurred Image');subplot(3,2,3),imshow(boundaryReplicate),title('Replicate');subplot(3,2,4),imshow(boundary0),title('0-Padding');subplot(3,2,5),imshow(boundarysymmetric),title('symmetric');subplot(3,2,6),imshow(boundarycircular),title('circular');实验结果如下：(４)、运用foｒ循环，将加有椒盐噪声得图像进行1０次，20 次均值滤波,查瞧其特点,显示均值处理后得图像(提示：利用ｆspｅｃｉal 函数得’ａveｒagｅ＇类型生成均值滤波器)。

代码如下:h=fspecial('average');实验结果：（5）、对加入椒盐噪声得图像分别采用均值滤波法,与中值滤波法对有噪声得图像做处理，要求在同一窗口中显示结果。

实验结果为:（6)、自己设计平滑空间滤波器,并将其对噪声图像进行处理，显示处理后得图像。

代码如下:[m n]=size(img2);figure,subplot(1,2,1),imshow(img2);s=zeros(1,9);代码如下:h1=fspecial('average');J=imfilter(img2,h1);J2=medfilt2(img2);figure,subplot(1,3,1),imshow(img2),title('salt & pepper Noise');subplot(1,3,2),imshow(J),title('Averaging Filtering');subplot(1,3,3),imshow(J2),title('Median Filtering');for i=1:10J1=imfilter(img2,h);endfor j=1:20J2=imfilter(img2,h);endfigure,subplot(1,3,1),imshow(img2),title('salt & pepper Noise');subplot(1,3,2),imshow(J1),title('10 Average Filtering');subplot(1,3,3),imshow(J2),title('20 Average Filtering');for i=2:1:m-1for j=2:1:n-1h=1;for p=i-1:1:i+1for q=j-1:1:j+1s(h)=img2(p,q);h=h+1;endends=sort(s);I(i,j)=s(5);endendsubplot(1,2,2),imshow(I);实验结果：2、锐化空间滤波(1）读出一幅图像，采用3×3得拉普拉斯算子w = [ 1, 1， 1; １– 8 1; 1, 1， 1］对其进行滤波.代码如下:img=imread('lena、png');img=rgb2gray(img);img=im2double(img);w=[1,1,1;1,-8,1;1,1,1]k=conv2(img,w,'same');imshow(k);实验结果为:（２）编写函数w = gｅｎlaｐｌaｃiａn(n）,自动产生任一奇数尺寸n 得拉普拉斯算子,如５×５得拉普拉斯算子w ＝ [ １ 1 １１ 11 １１ 1 1１１—24 １ 1１ 1 1 １ 1代码如下:nｕm=input('plｅａse enｔeｒ anynum:')；ｎ=num;Ｗ=oneｓ(ｎ,n); foｒｉ=1:ｎfoｒj＝１：nif（i==fｉx(n/２）+1 &＆ j＝=fｉｘ(n/2)+1)W（i，j)=ｎ＊n—1；ｅndenｄenｄdｉsｐｌay （W)；代码运行结果为：进行锐化滤波,并利用式完成图像得锐化增强,观察其有何不同，要求在同一窗口中显示。