《matlab-空域和频域进行滤波处理甄选.docx》由会员分享，可在线阅读，更多相关《matlab-空域和频域进行滤波处理甄选.docx（15页珍藏版）》请在三一文库上搜索。
1、精品word.仅供参考matlab-空域和频域进行滤波处理(优选.)rd图像平滑处理的空域算法和频域分析1技术要求对已知图像添加高斯白噪声，弁分别用低通滤波器(频域法)和邻域平均法(空域法)对图像进行平滑处理 (去噪处理),弁分析比较两种方法处理的效果。2基本原理2.1 图像噪声噪声在理论上可以定义为“不可预测，只能用概率统计方法来认识的随机误差”O实际获得的图像一般都因受到某种干扰而含有噪声。引起噪声的原因有敏感元器件的内部噪声、相片底片上感光材料的颗粒、 传输通道的干扰及量化噪声等。噪声产生的原因决定了噪声的分布特性及它和图像信号的关系。根据噪声和信号的关系可以将其分为两种形式：(1)加性
2、噪声。有的噪声与图像信号g(x,y)无关，在这种情况下，含噪图像 f(x,y)可表示为f(x,y)=g(x,y)+n(x,y)(2)乘性噪声。有的噪声与图像信号有关。这又可以分为两种情况：一种是某像素处的噪声只与该像 素的图像信号有关，另一种是某像点处的噪声与该像点及其邻域的图像信号有关，如果噪声与信号成正比， 则含噪图像f(x,y)可表示为f(x,y)=g(x,y)+n(x,y)g(x,y)另外，还可以根据噪声服从的分布对其进行分类，这时可以分为高斯噪声、泊松噪声和颗粒噪声等。如 果一个噪声，它的幅度分布服从高斯分布，而它的功率谱密度又是均匀分布的，则称它为高斯白噪声，一 般为加性噪声。2.
3、2 图像平滑处理技术平滑技术主要用于平滑图像中的噪声。平滑噪声在空间域中进行，其基本方法是求像素灰度的平均值 或中值。为了既平滑噪声又保护图像信号，也有一些改进的技术，比如在频域中运用低通滤波技术。(1)空域法在空域中对图像进行平滑处理主要是邻域平均法。这种方法的基本思想是用几个像素灰度的平均值来代 替每个像素的灰度。假定有一幅N*N个像素的图像f(x,y),平滑处理后得到一幅图像g(x,y)o g(x,y)由下式决定f (m, n)g(x,y) (m,n) S m1 / 9doc格式可编辑精品word.仅供参考式中，x,y=0,1,2,，N-1; S是(x,y)点邻域中点的坐标的集合，但其中
4、不包括 (x,y)点；M是集合内坐标点的总 数。上式说明，平滑化的图像 g(x,y)中每个像素的灰度值均由包含在(x,y)的预定邻域中的f(x,y)的几个像素的 灰度值的平均值来决定。(2)频域法低通滤波法是一种频域处理方法。在分析图像信号的频率特性时，一幅图像的边缘、跳跃部分以及颗 粒噪声代表图像信号的高频分量，而大面积的背景区则代表图像信号的低频分量。用滤波的方法滤除其高 频部分就能去掉噪声，使图像得到平滑。由卷积定理可知G(u,v) H(u,v)?F(u,v)其中F(u,v)是含有噪声的图像的傅立叶变换，G(u,v)是平滑处理后的图像的傅立叶变换，H(u,v)是传递函数,选择传递函数H(
5、u,v),利用H(u,v)使F(u,v)的高频分量得到衰减，得到G(u,v)后再经傅立叶反变换后就可以得到所希望的平滑图像 g(x,y)了。根据前面的分析，显然H(u,v)应该具有低通滤波特性，所以这种方法叫低通滤波法平滑化处理。常用的低通滤波器有如下几种：a.理想低通滤波器一个理想的二维低通滤波器有一个钞。它是一个规定的非负的量，叫做理想低通滤波器的截止频率。所谓理想低通滤波器是指以截频为书径的圆内的所有频率都能无损地通过，而在截频之外的频率分量完全被衰减。理想低通滤波器可以用计算机模拟实现，但是却不能用电子元器件实现。b.布特沃斯(Butterworth)低通滤波器一个n阶布特沃斯低通滤波
6、器的传递函数由下式表示2nH(u,v) 1 1DDo式中Do为截止频率。布特沃斯低通滤波器又称最大平坦低通滤波器。它与理想低通滤波器不同，它的通带与阻带之间没有明显的不连续性。也就是说，在通带和阻带之间有一个平滑的过度带，通常把H(u,v)下降到某一值的那一点定为截止频率。一般情况下常常到下降到H(u,v)最大值的那一点为截止频率点，该点也常称为半功率点2这样上式可修改成2nH(u,v) 1 1.,2 1 Doc.指数低通滤波器在图像处理中常用的另一种平滑滤波器是指数低通滤波器，由于指数低通滤波器有更快的衰减率，所 以，经指数低通滤波器处理的图像比布特沃斯低通滤波器处理的图像稍模糊一些。1所示
7、3建立模型描述图像平滑处理用到两种方法，其中空域法使用邻域平均法处理，其流程图如图2 / 9doc格式可编辑精品word.仅供参考Imread()函 数读 取rgb2gray()函数将原图转1fimnoise()函数为原图加double()将加了噪声的图创建3*3模创建5*5模创建7*7模用模板对噪用模板对噪用模板对噪声图邻域平声图邻域平声图邻域平Image()显示Image()显示Image()显 示处理后的图处理后的图处理后的图图1邻域平均法模型图频域法用二阶布特沃斯低通滤波器对图像滤除噪声，其程序流程图如图2所示:设定巴特我司低通滤波器的参数，并求与噪声图的傅立叶imread()函数读S
8、ize ()函数获取图将原图转换为灰度图fft2()函数对噪声图进行傅立wgn()函数产生高斯白噪声并加设定巴特我司低通滤波债的参数户并求 9doc格式可编辑与噪声图的傅立叶精品word.仅供参考巴特沃斯低通滤波器传递函数与图ifft2()函数进行逆傅立叶 image()函数显示图像图2低通滤波法模型图4源程序代码(1)邻域平均法源程序clear;clc;close all;I=imread(tu.jpg); %读取图像f=rgb2gray(I);%转化成灰度图subplot(2,3,1);%分割2*3个窗口。取第一个窗口，下面在第一个窗口处显示图像imshow;显示灰度图colormap(g
9、ray);title(原始图)给显示的图像命名为“原始图”J=imnoise(f,gaussian,0,0.09);%合原始图力口入参数为 0.09的高斯白噪声subplot(2,3,2);%分割2*3个窗口。取第二个窗口，如果下面有绘图语句，就表示要%在第二个窗口中绘图imshow(J);%显示加了高斯白噪声的图像Jtitle(噪声图)命名为“噪声图”4 / 9doc格式可编辑精品word.仅供参考m n=size(f);幅取灰度图的大小f=double;转换f为双精度型c=1/9*1 1 1;1 1 1;1 1 1; %3*3 模板for i=1:mfor j=1:nL=f(i:i,j:j
10、).*c; % 求点积G(i,j)=sum(sum(L);% 和endendsubplot(2,3,3);image(G);% 取第三个窗 口title(3*3模板)命名为“ 3*3模板”c=1/25*1 1 1 1 1;1 1 1 1 1;1 1 1 1 1; %5*5 模板for i=1:mfor j=1:nL=f(i:i,j:j).*c; % 求点积G(i,j)=sum(sum(L);%S 和endendsubplot(2,3,4);image(G);%在第四个窗口中显示图像title(5*5模板)命名为“ 5*5模板”c=1/49*1 1 1 1 1 1 1;1 1 1 1 1 1 1
11、;1 1 1 1 1 1 1; %7*7 模板 for i=1:mfor j=1:nL=f(i:i,j:j).*c; % 求点积G(i,j)=sum(sum(L);%S 和endendsubplot(2,3,5);image(G);%在第五个窗口中显示图像title(7*7模板)命名为“ 7*7模板”5 / 9doc格式可编辑精品word.仅供参考(2)低通滤波法源程序 clear;clc;close all;I=imread(tu.jpg); %读取图像m n p=size(I);%取图像矩阵大小I=double(I);将原来的图像矩阵转换为双精度型的I=I(1:m,1:n,1).*0.3+
12、I(1:m,1:n,2).*0.51+I(1:m,1:n,1).*0.11;% 将原图转换为灰度图subplot(2,3,1);image(I);colormap(gray);%分割 2*3 个窗口。取第一个窗口，显示灰度图title(原始图)命名为“原始图”Noise=wgn(m,n,25);%产生25dBm的高斯白噪声New = Noise + I;%各高斯白噪声与原始信号叠加subplot(2,3,2);image(New);%在第二个窗口中显示加了噪声的图像title(噪声图)命名为“噪声图”g=fft2(New);% 傅立叶变换g=fftshift(g);%转换数据矩阵M,N=siz
13、e(g);% 取矩阵大小nn=2;%定义二阶巴特沃斯(Butterworth)低通滤波器d0=20;%截止频率为20m=fix(M/2); n斗ix(N/2);圆整函数for i=1:Mfor j=1:Nd=sqrt(i-m)A2+(j-n)A2);h=1/(1+0.414*(d/d0)A(2*nn);%计算低通滤波器传递函数result(i,j)=h*g(i,j);%传递函数与信号傅立叶变换函数相乘end6 / 9doc格式可编辑精品word.仅供参考end result=ifftshift(result);J2=ifft2(result);%逆傅立叶变换J3=uint8(real(J2);
14、%各输出图彳t转换为8位数据存储subplot(2,3,3);image(J3); %在第三个窗口中显示滤波后的图像title(截止频率20 )%命名为“截止频率 20”nn=2;%定义二阶巴特沃斯(Butterworth)低通滤波器d0=100;%截止频率为100m=fix(M/2); n=fix(N/2);圆整函数for i=1:Mfor j=1:Nd=sqrt(i-m)A2+(j-n)A2);h=1/(1+0.414*(d/d0)A(2*nn);%计算低通滤波器传递函数result(i,j)=h*g(i,j);%传递函数与信号傅立叶变换函数相乘endendresult=ifftshift
15、(result);J2=ifft2(result);%二维逆傅立叶变换J3=uint8(real(J2); %专换为8位存储的图像subplot(2,3,4);image(J3); %在第四个窗口中显示绿波后的图像title(截止频率100 )%命名为“截止频率100” 5调试过程及结论3所示:在MatLab中输入邻域平均法程序代码后运行程序，得到输出结果如图精品word.仅供参考图3邻域平均法输出图像由噪声图与滤波后的图像对比可看出，邻域平均法对抑制噪声有明显的效果，但随着邻域的加大，就 是随着模板的加大，图像的模糊程度也愈加严重。在MatLab中输入低通滤波法程序代码后运行程序，得到输出结
16、果如图4所示：图4低通滤波法输出图像8 / 9doc格式可编辑精品word.仅供参考巴特沃斯低通滤波器的去噪效果与所选的截止频率有关，由于一幅图像的边缘、跳跃部分以及颗粒噪 声代表图像信号的高频分量，而大面积的背景区则代表图像信号的低频分量所以，截频设的太低会使图像 变得越模糊，因为图像的许多细节信息也被滤掉了。6心得体会这次能力强化训练是图像平滑处理的空预算法和频域分析，需要有一定的对数字图像处理方面的知识，由于没有学习过数字图像处理这门课，已开始题目布置下来时感到无从下手，查阅了大量资料，对图像平滑处理有了一定的了解，对于 MatLab刚开始接触，许多语法、函数都不会使用，于是和同学互相讨
17、 论，通过参考大量资源，站在巨人的肩膀上稍加整合创新，最终得出了自己的程序，程序调试运行后，发 现最开始根本无法显示图像， 检查原因，原来是对MatLab的符号规则不清楚，MatLab中的标点符号必须是 在英文状态下输入，这在编写程序时必须小心！修改后得到了正确的显示，但是图像效果对比不明显，通 过不断摸索试探，使用不同的图像显示函数，得到了对比鲜明的图像，对分析不同方法的优缺点提供了很 大方便。这次数字图像处理课程设计虽然时间短暂，但我在短时间内学到了许多以前没有接触过的数字图像处 理方面的知识，并最终完成了程序设计，得到正确的结果，学到了许多 MatLab的知识。为了完成这次课程 设计，我
18、花费了很多精力学习 MatLab,也查阅了很多关于 MatLab的书籍，了解了许多 MatLab语法和函数 功能。在设计过程中，我也碰到了许多问题，走了很多弯路，主要是因为对数字图像处理的原理还弄得不 透彻，对MatLab还不熟练，犯了一些常规错误，但最后通我查阅相关书籍和请教同学，问题都迎刃而解。 在解决问题过程中，加深了对知识的理解，最重要的是学会了自己解决问题的能力，受益非浅。7参考文献1何东健，耿楠，张义宽.数字图像处理.西安电子科技大学出版社.20032阮秋琦.数字图像处理学.电子工业出版社.20073罗军辉，冯平，哈WA.MATLAB7.0在图像处理中的应用.机械工业出版社.2005赠人玫瑰，手留余香。感谢您使用本店文档您的满意是我们的永恒的追求！（本句可删）9 / 9doc格式可编辑