实验三数字图像的空间域滤波模板

1、实验三、四 数字图像的空间域滤波和频域滤波1实验目的1.掌握图像滤波的基本定义及目的。2.理解空间域滤波的基本原理及方法。3.掌握进行图像的空域滤波的方法。4.掌握傅立叶变换及逆变换的基本原理方法。5.理解频域滤波的基本原理及方法。6.掌握进行图像的频域滤波的方法。2实验基本原理处理图像每一个像素的取值 空域滤波基本上是让图像在 从而改变输出图像的频率分1. 空间域增强空间域滤波是在图像空间中借助模板对图像进行领域操作， 都是根据模板对输入像素相应领域内的像素值进行计算得到的。 频域空间内某个范围的分量受到抑制， 同时保证其他分量不变， 布，达到增强图像的目的。空域滤波一般分为线性滤波和非线性

2、滤波两类。 线性滤波器的设计常基于对傅立叶变换 的分析， 非线性空域滤波器则一般直接对领域进行操作。 各种空域滤波器根据功能主要分为 平滑滤波器和锐化滤波器。平滑可用低通来实现， 平滑的目的可分为两类： 一类是模糊， 目 的是在提取较大的目标前去除太小的细节或将目标内的小肩端连接起来；另一类是消除噪 声。锐化可用高通滤波来实现， 锐化的目的是为了增强被模糊的细节。 结合这两种分类方法， 可将空间滤波增强分为四类：线性平滑滤波器（低通） 非线性平滑滤波器（低通） 线性锐化滤波器（高通） 非线性锐化滤波器（高通）将模板在图中移动，并将模板中心与图中某个像素位置重合； 将模板上的系数与模板下对应的像

3、素相乘； 将所有乘积相加；将和（模板的输出响应）赋给图中对应模板中心位置的像素。 平滑滤波器线性平滑滤波器空间滤波器都是基于模板卷积，其主要工作步骤是：1）2）3）4）2.1）线性低通平滑滤波器也称为均值滤波器，这种滤波器的所有系数都是正数，对3X39。的模板来说，最简单的是取所有系数为1，为了保持输出图像任然在原来图像的灰度值范围内，模板与象素邻域的乘积都要除以MATLAB 提供了 fspecial 函数生成滤波时所用的模板，并提供 filter2 函数用指定的滤波器模板对图像进行运算。函数fspecial的语法格式为：h=fs pecial(t yp e);h=fs pecial(t yp

4、e,p arameters);其中参数type指定滤波器的种类，parameters是与滤波器种类有关的具体参数。Typepa珂砂:eiw说明旳値溥诫”躺舉黑以片芍*剧前工亍为标甌否则由两元如网捷XlM#指;M 邻域的?1和列薮，d:kradru;有(珂11：乜+讣令!边的0睢均值谑粧體标准嵋筮为A小为址伽的格斷限逓滤谍器bplaciiiDalpha毎数由汨pheeLO)抉址們二细拉普拉斷滤盜lOF标Rr -In寸対齐黄小対tiJz#的高斷誌很啸转对称过氏即卡udticaL3. tbita按腐曖凶E饲桟訪1 r竜的运动谑浚器piSWItT杞近似汁訝壷捧度的水平边労强调薛子sobl近似汁舅垂宜梯

5、嗖光潸散应的水平Ji缘强调S子lusluipalplu根据ilpLa决定的拉氏算子Ih注的掩接谑世辭表2.1 MATLAB中预定义的滤波器种类MATLAB 提供了一个函数imnoise来给图像增添噪声，其语法格式为：J=i mno ise(l,typ e);J=i mno ise(l,type,p arameters);参数type指定噪声的种类，parameters是与噪声种类有关的具体参数。参数的种类见 表 2.2。参薮说明gaussiaiim. V肉旳为帕万状为V的崗斯呀向localvarV肉他为 方差为Y的高斯白嗓声passiou无泊松噪芦M pepperspeckleV肉社勺乩 方茎

6、为的沟匀分布随机噪晴表22噪声种类及参数说明2) 非线性平滑滤波器其滤波原理与均值滤波器方法类似，但计中值滤波器是一种常用的非线性平滑滤波器,然后选择改组的中间值作算的非加权求和，而是把领域中的图像的象素按灰度级进行排序， 为输出象素值。MATLAB提供了 medfilt2函数来实现中值滤波，其语法格式为：B=medfilt2(A,m n);B=medfilt2(A);其中，A是原图象，B是中值滤波后输出的图像。m n指定滤波模板的大小，默认模 板为3X 3。3.锐化滤波器图像平滑往往使图像中的边界、轮廓变得模糊，为了减少这类不利效果的影响,用图像锐化技术，使图像的边缘变得清晰。1)线性锐化滤

7、波器线性高通滤波器是最常用的线性锐化滤波器。 的系数都是负的，所有的系数之和为0。对3X 3的模板来说，典型的系数取值为：这种滤波器的中心系数都是正的,需要利而周围-1-1-1;-1 8 -1;-1-1-1laplacian,0.5）得到的拉普拉斯算子所以利用微分可以锐化图像。图像处理中sobel 算子、Prewitt事实上这是拉普拉斯算子。语句h=-fs pecial(为：h =-0.3333 -0.3333 -0.3333-0.3333 2.6667 -0.3333-0.3333 -0.3333 -0.33332) 非线性锐化滤波邻域平均可以模糊图像，因为平均对应积分，最常用的微分方法是利

8、用梯度。常用的空域非线性锐化滤波微分算子有算子、log算子等。4.频域增强频域增强是利用图像变换方法将原来的图像空间中的图像以某种形式转换到其他空间 中，然后利用该空间的特有性质方便地进行图像处理，最后再转换回原来的图像空间中，从而得到处理后的图像。频域增强的主要步骤是：选择变换方法，将输入图像变换到频域空间。在频域空间中，根据处理目的设计一个转移函数，并进行处理。将所得结果用反变换得到增强的图像。常用的频域增强方法有低通滤波和高通滤波。5. 低通滤波图像的能量大部分集中在幅度谱的低频和中频部分，而图像的边缘和噪声对应于高频部分。因此能降低高频成分幅度的滤波器就能减弱噪声的影响。由卷积定理，在

9、频域实现低通滤波的数学表达式：G(u,v) =H(u,v)F(u,v)1) 理想低通滤波器(ILPF )|o2)巴特沃斯低通滤波器（BLPF ）3)指数型低通滤波器（ELPF ）高通滤波6.由于图像中的细节部分与其高频分量相对应，所以高通滤波可以对图像进行锐化处理。 高通滤波与低通滤波相反，它是高频分量顺利通过，通滤波器相似，其转移函数分别为：理想高通滤波器（IHPF）使低频分量受到削弱。高通滤波器和低1)D(“) D2)巴特沃斯高通滤波器（BLPF ）3)指数型高通滤波器（ELPF ）图像经过高通滤波处理后，会丢失许多低频信息， 所以图像的平滑区基本上会消失。所以，可以采用高频加强滤波来弥补

10、。高频加强滤波就是在设计滤波传递函数时，加上一 个大于0小于1的常数C,即：H (u,v) =H(u,v)+c3.实验内容与要求1. 平滑空间滤波：1）读出eight.tif这幅图像，给这幅图像分别加入椒盐噪声和高斯噪声后并与前一张图 显示在同一图像窗口中。2）对加入噪声图像选用不同的平滑（低通）模板做运算，对比不同模板所形成的效果，要求在同一窗口中显示。3）使用函数imfilter时，分别采用不同的填充方法（或边界选项，如零填 充、rep licate ” symmetric、 circular）进行低通滤波，显示处理后的图像。4）运用for循环，将加有椒盐噪声的图像进行10次，20次均值滤

11、波，查看其特点，显示均值处理后的图像(提示：利用fspecial函数的averag类型生成均值滤波器)。5) 对加入椒盐噪声的图像分别采用均值滤波法，和中值滤波法对有噪声的图像做处 理，要求在同一窗口中显示结果。6)2.1)自己设计平滑空间滤波器，并将其对噪声图像进行处理，显示处理后的图像。锐化空间滤波读出blurry_moon.tif这幅图像，采用3X3的拉普拉斯算子 w = 1, 1, 1; 1 -8 1; 1, 1,1对其进行滤波。2) 编写函数 的拉普拉斯算子w = genlap lacia n(n)，自动产生任一奇数尺寸n的拉普拉斯算子，如5X5111111111111-241111

12、111111113)分别采用9X 9,15X15和25 X 25大小的拉普拉斯算子对blurry_moon.tif进行锐化滤波，并利用式g(x, y) = f(X, y) -2 f(x, y)完成图像的锐化增强，观察其有何不同，要求在同一窗口中显示。4)5)采用不同的梯度算子对 blurry_moon.tif进行锐化滤波，并比较其效果。 自己设计锐化空间滤波器，并将其对噪声图像进行处理，显示处理后的图像;3.傅立叶变换读出woman.tif这幅图像，对其进行快速傅立叶变换，分别显示其幅度图像和相位1)图像。仅对相位部分进行傅立叶反变换后查看结果图像。仅对幅度部分进行傅立叶反变换后查看结果图像。

13、将图像的傅立叶变换 F置为其共轭后进行反变换，比较新生成图像与原始图像的差2)3) 异。4.平滑频域滤波设计理想低通滤波器、巴特沃斯低通滤波器和高斯低通滤波器，截至频率自选，分1)别给出各种滤波器的透视图。2)读出test_pattern.tif这幅图像，分别采用理想低通滤波器、巴特沃斯低通滤波器和高 斯低通滤波器对其进行滤波(截至频率自选)，再做反变换，观察不同的截止频率下采用不同低通滤波器得到的图像与原图像的区别，特别注意振铃效应。(提示:1)在频率域滤波同样要注意到填充问题；2)注意到(-1)x+y;)5.锐化频域滤波1) 设计理想高通滤波器、巴特沃斯高通滤波器和高斯高通滤波器，截至频率

14、自选，分 别给出各种滤波器的透视图。2) 读出test_pattern.tif这幅图像，分别采用理想高通滤波器、巴特沃斯高通滤波器和高斯高通滤波器对其进行滤波(截至频率自选)，再做反变换，观察不同的截止频率下采用不同高通滤波器得到的图像与原图像的区别。1.实验具体实现1.平滑空间滤波：1）读出eight.tif这幅图像，给这幅图像分别加入椒盐噪声和高斯噪声后并与前一张图 显示在同一图像窗口中。l=imreadCeight.tif);imshow(l);%no ise den sity=0.05J = imno ise(l,salt & pepp er,0.05);K= imn oise(l,g

15、aussia n,0.01,0.01);sub plot(131),imshow(l);sub plot(132),imshow(J);sub plot(133),imshow(K);Original ImageSalt&Pepper Noise Gaussian Noise图2.1初始图像及椒盐噪声图像、高斯噪声污染图2）对加入噪声图像选用不同的平滑（低通）模板做运算，对比不同模板所形成的效果,要求在同一窗口中显示。l=imread(eight.tif);H = fsp ecial(sobel);Sobel = imfilter(I,H,re pl icate);H = fsp ecial(

16、la pl acia n,0.4); lap = imfilter(l,H,re plicate);H = fsp ecial(gaussia n,3 3,0.5); gaussia n = imfilter(l,H,re plicate); sub plot(221),imshow(l);sub plot(222),imshow(Sobel);sub plot(223),imshow(la p); sub pl ot(224),imshow(gaussia n);图2.2原图像及各类低通滤波处理图像3）使用函数imfilter时，分别采用不同的填充方法（或边界选项，如零填 充、 replic

17、ate ” symmetric circular）进行低通滤波，显示处理后的图像。orig in alRGB = imread( peppers.pn g);h = fsp ecial(moti on, 50, 45);%moti on blurredfilteredRGB = imfilter(origi nalRGB, h);bou ndaryRe plicateRGB = imfilter(origi nalRGB, h, re pl icate); bou ndaryORGB = imfilter(origi nalRGB, h, x);bou ndary0RGB = imfilter

18、(origi nalRGB, h, 0);bou ndarysymmetricRGB = imfilter(origi nalRGB, h, symmetric); bou ndarycircularRGB = imfilter(origi nalRGB, h, circular);Original imageMoticn Blurred Image图2.3原图像及运动模糊图像图2.4函数imfilter各填充方式处理图像4) 运用for循环，将加有椒盐噪声的图像进行10次，20次均值滤波，查看其特点，显示均值处理后的图像。J = imno ise(l,salt & pepp er,0.05)

19、;h=fs pecial(average);%Averag ing Filteri ngJ1=imfilter(J,h);for i=1:10J1=imfilter(J,h); end for i=1:20J2=imfilter(J,h);Salt&LPeppur Noise10 Averaging Filtering 20 Averagirg Filtering图2.5椒盐噪声污染图像经10次、20次均值滤波图像由图2.5可得，20次滤波后的效果明显好于 10次滤波，但模糊程度也更强。5) 对加入椒盐噪声的图像分别采用均值滤波法，和中值滤波法对有噪声的图像做处理，要求在同一窗口中显示结果。%

20、Media n Filteri ngh=fs pecial(average);J1=imfilter(J,h);J2=medfilt2(J);%Averag ing Filteri ngAveraging FilteringMedian FiltenngOrigiral Image图2.6椒盐噪声污染图像及均值、中值滤波图像从图2.6中可以看出，对于椒盐噪声污染的图像处理，中值滤波效果要明显好于均值滤波。经均值滤波器处理后的图像比均值滤波器中结果图像更加模糊。6）设计平滑空间滤波器，并将其对噪声图像进行处理，显示处理后的图像。 domai n=0 0 8 0 0;0 0 8 0 0;8 8 8

21、 8 8;0 0 8 0 0;0 0 8 0 0;K1= ordfilt2(J,5,domai n);Salt&Pepper Noise Image 5*5 Smoothing Spatial Fitered Image 讦猜:以弗总处紀-孑 y 痴沁晋逼f墨翳;、 -I -IJ二:沁辽冷碎y.-F仇.:mhv；.爛5礙丁:- 品图2.7椒盐噪声污染图像及5*5平滑滤波器掩模掩模值为 w=1/25*1 1 1 1 1;1 1 1 1 1;1 1 1 1 1;1 1 1 1 1;1 1 1 1 1SaltgiPepper Noised Image:沆:芯竺堀质吃菩能:洽*： ； J f:右Wi

22、Il .* ij. rj|i T I -i 丹；严-*飞品-,-3羅S醪:吃*出曲 M；空务 -. 丹礙卫应.曲霸SSSL.:- 叱沁.d%總轟U申:丁:丁沈沦T上- V词禅娜:.-朋VLow pass Filtered Image图2.8椒盐噪声污染图像及5*5平滑滤波器掩模掩模值为 w= 0 0 8 0 0;0 0 8 0 0;8 8 8 8 8; 0 0 8 0 0;0 0 8 0 02. 锐化空间滤波1) 采用3X 3的拉普拉斯算子 w = 1, 1, 1; 1 -8 1; 1, 1, 1滤波 I=imread(blun7_ moon .tif);T=double(l);sub plo

23、t(1,2,1),imshow(T,);title(Origi nal Image);w =1,1,1;1,1,1;K=c on v2(T,w,same);1,-8,1;图2.9初始图像与拉普拉斯算子锐化图像，自动产生任奇数尺寸 n的拉普拉斯算子，如 5X 5的拉2)编写函数 w = gen la placia n(n)普拉斯算子：-241 11111fun cti on w = genlap lacia n(n) %Comp utes the Lap lacia n op erator w = ones(n);x = ceil( n/2);w(x, x) = -1 * (n * n - 1)

24、;3） 分别采用5X 5, 9X 9, 15X 15和25X 25大小的拉普拉斯算子对blurry_moon.tif进行锐2化滤波，并利用式g（x，y） =f（X，y） f（X，y）完成图像的锐化增强，观察其有何不同，要求在同一窗口中显示。不同尺寸拉普拉斯算子滤波以及图像增强w1 = genlap lacia n( 5); I=imreadCblurry_ moon .tif); T=double(l);K=co nv2( T,w1,same);J=T-K;图2.10初始图像与不同拉普拉斯算子锐化图像图像锐化的实质是将原图像与梯度信息叠加，相当于对目标物的边缘进行了增强。图2.11拉普拉斯算子

25、锐化与锐化增强图像4)采用不同的梯度算子对blurry_moon.tif进行锐化滤波，并比较其效果l,ma p=imread(blurry_ moon .tif);匸double(l);Gx,Gy=gradie nt(l);G=sqrt(Gx.*Gx+Gy.*Gy);% gradie nt calculati on% matrixJ1=G;% gradie nt1imshow(J1,ma p);J2=I;% gradie nt2K=fin d(G=7);J2(K)=G(K); imshow(J2,ma p);J3=I;K=fin d(G=7);J3(K)=255; imshow(J3,ma p

26、);J4=I;% gradie nt3% gradie nt4K=fin d(G=7);J4(K)=255; imshow(J4,ma p);J5=I;K=fin d(G=7);J5(Q)=255; imshow(J5,ma p);% gradie nt5图2.12原始图像与不同梯度子锐化图像图 2.13一 i丿曦號_刃” -I B原始图像与不同边缘锐化图像jSc m作为二阶微分算子，拉普拉斯变换在图像细节的增强处理上有明显的优点，但会产生更多的噪声。梯度变换在灰度变化区域的响应更强，但对噪声和细节的响应比拉普拉斯变换弱。5)自己设计锐化空间滤波器，并将其对噪声图像进行处理，显示处理后的图像;

27、3.1)domain=8 8 0 8 8;8 8 0 8 8;傅立叶变换读出woman.tif这幅图像，对其进行快速傅立叶变换，分别显示其幅度图像和相位图像。0 0 0 0 0;8 8 0 8 8;8 8 0 8 8;K1= ordfilt2(J,5,domai n);F=imread(woma n.tif);F1=fft2(F);%amp litude sp ectrumF2=log(1+abs(F1);F3=fftshift(F1); imshow(log(1+abs(F3),);F4=a ngle(F1);%p hase sp ectrumOriginal ImageAmplitude

28、SpectrumP base Ep ectrumL-A. H MA亠ZeroTreqiuency Centralized Spectrum图2.14原始图像与快速傅立叶变换图像2)仅对相位部分进行傅立叶反变换后查看结果图像。F1=fft2(F);i=sqrt(-1);f2=ifft2(ex p( i*a ngle(F1);imshow(real(f2),);Orginal ImageRestoration Image图2.15原始图像与对全部信息进行傅立叶逆变换结果Phase-based Restoration Image图2.16仅对相位信息进行傅立叶逆变换结果相位谱决定了图像信号中各频率分

29、量的位置。3) 仅对幅度部分进行傅立叶反变换后查看结果图像。f1=ifft2(abs(F1);imshow(log(1+abs(f1),);图2.17仅对幅度信息进行傅立叶逆变换结果4)将图像的傅立叶变换F1=fft2(F);F置为其共轭后进行反变换，比较新生成图像与原始图像的差异。F2=log(1+abs(F1);F3=fftshift(F1);F4=a ngle(F1);%amp litude sp ectrum%p hase sp ectrumF5=-F4F6= double(F3*ex p(F4);F7=ifft2(F6);imshow(real(F7),);%thecomp lexc

30、on jugateof the fouriertra nsform%in verse fourier tra nsformOriginal ImageComjugte image图2.18共轭傅立叶逆变换结果傅立叶变换的相位谱为对称的，原变换与其共轭变换间仅频率谱互为相反。4.平滑频域滤波1）设计理想低通滤波器、巴特沃斯低通滤波器和高斯低通滤波器 理想低通滤波器l=imread(test_ patter n.tif);f=double(l);% chage into double as MATLAB doesn% of image in un sig ned int typeg=fft2(f)

31、;g=fftshift(g);M,N=size(g);d0=100;m=fix(M/2); n=fix(N/2);for i=1:Mt suppor calculation% fourier tra nsform% zero-freque ncy area cen tralized%cutoff freque ncyforj=1:Nd=sqrt(i-m)A2+(j-门)人2);if(d=d0)h=1;else h=0;endresult(i,j)=h*g(i,j);end end result=ifftshift(result);J1=ifft2(result);J2=ui nt8(real(

32、J1); imshow(J2)巴特沃斯低通滤波器（二阶）l=imread(test_ patter n.tif); f=double(l);g=fft2(f);g=fftshift(g);M,N=size(g);% 2-grade Butterworth low pss filternn=2;d0=100;m=fix(M/2); n=fix(N/2);for i=1:Mforj=1:Nd=sqrt(i-m)2+(j-n)人2);h=1/(1+0.414*(d/d0)A(2* nn); % filter tran sform fun ction%h=1.心+(d./d0)A(2* n)%h=ex

33、 p(-人2)./(2*90人2); result(i,j)=h*g(i,j);endend result=ifftshift(result);J1=ifft2(result);J2=ui nt8(real(J1); imshow(J2);高斯低通滤波l=imread(test_ patter n.tif); f=double(l);g=fft2(f);g=fftshift(g);M,N=size(g);d0=100;m=fix(M/2); n=fix(N/2);for i=1:Mforj=1:N% gaussia n filter tran sformd=sqrt(i-m)A2+(j-门)人

34、2); h=ex p(-(dA2)./(2*(d0A2);result(i,j)=h*g(i,j);endend result=ifftshift(result);J1=ifft2(result);J2=ui nt8(real(J1);图2.19理想低通滤波器透视图图2.20巴特沃斯低通滤波器透视图图2.21高斯低通滤波器透视图由各类低通滤波器透视图可见，高斯滤波器剖面线最平滑， 二阶巴特沃斯低通滤波器函数剖面线较为紧凑，而理想滤波器完全为圆筒状结构，未考虑选择范围内不同信息频率的差别化处理。Original imageIdeal Lowpass Filter da=15 VL r.yI I

35、I ；遜岭I I I滴泯M d. H 3 d H H.alldeal Lowpass Filter d0=30Ideal Lowpags Filter dQ=1OO - a IIIIIIIII i 14 ii a 3 1ai I* i*i-.a iiiiiiii图 2.22理想低通滤波器滤波效果(d0=15,30,100)当截止频率do = 15时，当do = 30时，滤波后的图像模糊减弱，能分辨出字母与图形轮廓，但由于理想低通滤 波器在频率域的锐截止特性，滤波后的图像仍有较明显的振铃现象。当do = 100时，滤波后的图像比较清晰，但高频分量损失后，图像边沿与文字变的有 些模糊，在图像的边框(如条带和矩形轮廓)附近仍有轻微振铃现象。滤波后的图像模糊，难以分辨，振铃现象明显。Original Image