图像频率域低通滤波处理程序设计(常用版).doc
《图像频率域低通滤波处理程序设计(常用版).doc》由会员分享,可在线阅读,更多相关《图像频率域低通滤波处理程序设计(常用版).doc(117页珍藏版)》请在咨信网上搜索。
1、图像频率域低通滤波处理程序设计(常用版)(可以直接使用,可编辑 完整版资料,欢迎下载)专业综合课程设计任务书学生姓名: 陈德松 专业班级: 电信 0901班 指导教师: 黄朝兵 工作单位: 信息工程学院 题 目:图像频率域低通滤波处理程序设计初始条件: (1)提供实验机房及其matlab软件; (2)数字图像处理的基本理论学习。要求完成的主要任务:(包括课程设计工作量及其技术要求,以及说明书撰写等具体要求)(1)掌握频率域处理的基本原理,利用matlab设计程序完成以下功能; (2)选择一幅256级的灰度图像; (3)对该图像作Fourier变换,得到其频谱图像; (4)按照二阶Butterw
2、orth低通滤波器的表达式设计传递函数,对该频谱图像进行低通滤波,对滤波后的频谱图像作逆Fourier变换得到滤波结果,显示结果图; (5)按照高斯低通滤波器的表达式设计传递函数,对该频谱图像进行低通滤波,对滤波后的频谱图像作逆Fourier变换得到滤波结果,显示结果图; (6)对二种滤波器的滤波结果进行分析比较; (7)要求阅读相关参考文献不少于5篇; (8)根据课程设计有关规范,按时、独立完成课程设计说明书。时间安排:(1) 布置课程设计任务,查阅资料,确定方案四天; (2) 进行编程设计一周; (3) 完成课程设计报告书三天;指导教师签名: 年 月 日系主任(或责任教师)签名: 年 月
3、日目录摘 要I1 MATLAB的简介12 原理与实现22.1频率域增强基本理论22.2 傅立叶变换32.3频率域巴特沃兹(Butterworth)低通滤波器42.4频率域高斯(Gaussian)低通滤波器63 程序设计73.1算法设计(程序设计流程图)73.2 对灰度图像进行Fourier变换的程序73.3 二阶巴特沃斯(Butterworth)低通滤波程序83.4 高斯(Gaussian)低通滤波程序94结果与分析114.1选择一幅256级的灰度图像114.2 对灰度图像进行Fourier变换后的频谱图134.3 二阶巴特沃斯(Butterworth)低通滤波结果与分析144.4 高斯(Ga
4、ussian)低通滤波结果与分析174.5两种滤波器的滤波结果的比较205 心得体会21参考文献22摘 要数字图像处理的目的之一是改善图像的视觉效果, 以便人眼或机器对图像进一步理解。而图像增强技术通过采用某种技术手段可以改善图像的视觉效果。其中一种技术是频率域增强技术,它是增强技术的重要组成部分。本设计重点介绍了用Matlab对一幅256级的灰度图像作Fourier变换,然后对变换后得到的频谱图像进行分别进行二阶Butterworth低通滤波和高斯低通滤波,然后对滤波后的频谱图像作逆Fourier变换得到滤波结果,显示结果图。本程序的设计很简洁,加上注释显得通俗易懂,而且本程序中很多地方都有
5、输入提示,有一定的人机交互。关键词:Matlab ,Fourier变换,Butterworth,高斯(Gaussian),低通滤波1 MATLAB的简介MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学计算软件的先进水平。 MATLAB和Math
6、ematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等,主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。2 原理与实现2.1频率域增强基本理论不对Fourier变换(FT)和图像的频率域处理技术有所了解,就不可能完全理解图像增强这个最基本的图像处理任务。频域增强指在图像的频率域内,对图像的变换系数(频率成分)直接进行运算,然后通过Fourier逆变换以获得图像的增强效果。一般来说,图像的边缘和噪声对应Fourier变换
7、中的高频部分,所以低通滤波能够平滑图像、去除噪声。图像灰度发生聚变的部分与频谱的高频分量对应,所以采用高频滤波器衰减或抑制低频分量,能够对图像进行锐化处理。卷积理论是频域技术的基础,设函数f (x, y)与算子h(x, y)的卷积结果是g(x,y),即g(x, y) = h(x, y) * f (x, y),那么根据卷积定理在频域有: 其中G(u, v),H(u, v),F(u, v)分别是g(x, y),h(x, y),f (x, y)的傅立叶(或其它)变换,H(u, v)是转移函数。在具体增强应用中,f (x, y)是给定的(所以F(u, v)可利用变换得到),需要确定的是H(u, v),
8、这样具有所需特性的 g(x, y) 就可算出 G(u, v) 而得到: 2.2 傅立叶变换傅里叶变换是将时域信号分解为不同频率的正弦信号或余弦函数叠加之和。傅立叶变换是数字图像处理技术的基础,其通过在时空域和频率域来回切换图像,对图像的信息特征进行提取和分析,简化了计算工作量,被喻为描述图像信息的第二种语言,广泛应用于图像变换,图像编码与压缩,图像分割,图像重建等。因此,对涉及数字图像处理的工作者,深入研究和掌握傅立叶变换及其扩展形式的特性,是很有价值得。把傅立叶变换的理论通其物理解释相结合,将有助于解决大多数图像处理问题。傅里叶变换可分为连续傅里叶变换、离散傅里叶变换、快速傅里叶变换。快速傅
9、里叶变换(FFT)是计算离散傅里叶变换(DFT)的快速算法。离散傅里叶变换运算量巨大,计算时间长,即运算时间很长。而快速傅里叶变换的提出将傅里叶变换的复杂度由降到了,很大程度上减少了计算量。 令,u=0,1,2,M-1则,2.3频率域巴特沃兹(Butterworth)低通滤波器物理上可实现(理想低通滤波器在数学上定义得很清楚,在计算机模拟中也可实现,但在截断频率处直上直下的理想低通滤波器是不能用实际的电子器件实现的)减少振铃效应,高低频率间的过渡比较光滑n阶Butterworth低通滤波器的传递函数为:式中,D0为截止频率。=1时,=0.5,它的特性是传递函数比较平滑,连续衰减,而不像理想滤波
10、器那样陡峭变化,即明显的不连续。因此采用该滤波器滤波在抑制噪声的同时,图像边缘的模糊程度大大减小,没有振铃效应产生,如图2.1所示。图2.1 Butterworth低通滤波效果阶数对振铃现象的影响:阶数越高,越明显,如图2.2所示图2.2 阶数对振铃现象的影响 2.4频率域高斯(Gaussian)低通滤波器高斯(Gaussian)低通滤波器的传递函数为:其中,s为标准偏差。通过令sD0,我们可以根据截止参数D0得到表达式:当D(u,v)= D0时,滤波器H(u, v)由最大值1下降为0.607。GLPF没有振铃现象,但与阶数为2的BLPF相比,其通带要宽些,这样对应的空间滤波器的灰度级轮廓更窄
11、些,因而平滑效果要差些。以上两种滤波器,振铃现象从严重到无,但平滑效果从好到差,BLPF可以看成ILPF和GLPF的过渡,阶为1时与GLPF差不多,阶越高越接近BPLG.图2.3 高斯(Gaussian)低通滤波器对于不同的D0值的滤波效果3 程序设计程序除了有题目中要求的对该图像作Fourier变换,得到其频谱图像外,还扩展了调用了函数colorbar给一个坐标轴对象添加加一条色带,该色带可以指示出该图像中不同颜色的数据值,以显示出图像的灰度级,功能较多。程序中很多地方都有输入提示,有一定的人机交互。3.1算法设计(程序设计流程图)算法流程图如图3.1所示:图3.1程序设计流程图3.2 对灰
12、度图像进行Fourier变换的程序I=imread(Test picture.jpg);f=double(I); % 数据类型转换,MATLAB不支持图像的无符号整型的计算 g=fft2(f); figure,imshow(log(1 + abs(fftshift(g), )%将频谱显示为一幅图像title(原图Fourier变换后频谱图);3.3 二阶巴特沃斯(Butterworth)低通滤波程序此函数的功能为:使用巴特沃斯低通滤波原理实现图像柔和平滑,为便比较,通过设置不同的d0值运行程序,如果对各个d0一起运行,运行速度奇慢,所以对各个d0单独运行。I=imread(Test pictu
13、re.jpg); imshow(I);title(原图);s=fftshift(fft2(I);M,N=size(s); %分别返回s的行数到M中,列数到N中n=2; %对n赋初值%BLPF滤波,d0=15,30,80d0=input(请输入d0值:); %初始化d0,d0依次取15,30,80n1=floor(M/2); %对M/2进行取整n2=floor(N/2); %对N/2进行取整for i=1:M for j=1:N d=sqrt(i-n1)2+(j-n2)2); %点(i,j)到傅立叶变换中心的距离 h=1/(1+(d/d0)(2*n); %BLPF滤波函数 s(i,j)=h*s(
14、i,j); %BLPF滤波后的频域表示 endends=ifftshift(s); %对s进行反FFT移动%对s进行二维反离散的Fourier变换后,取复数的实部转化为无符号8位整数s=uint8(real(ifft2(s);figure; %创建图形图像对象imshow(s); %显示BLPF滤波处理后的图像title(BLPF滤波(d0=15));3.4 高斯(Gaussian)低通滤波程序此函数的功能为:使用高斯低通滤波原理实现图像越柔和平滑,为便比较,通过设置不同的d0值运行程序,如果对各个d0一起运行,运行速度奇慢,所以对各个d0单独运行。(1)高斯低通滤波的功能子程序GauseLo
15、wPass()function LowPass=GauseLowPass(ima,Do)%本函的功能为:实现高斯低通滤波原理% H(u,v)=exp(-D(u,v)2/(2*Do2); Row,Col=size(ima);ima=im2double(ima);for ix=1:Row %滤波中心化 for iy=1:Col ima(ix,iy)=(-1)(ix+iy)*ima(ix,iy); endend FourIma=fft2(ima);%中心变换后的傅立叶变换FRow,FCol=size(FourIma);for u=1:FRow%Do的高斯低通滤波器 for v=1:FCol D(u,
16、v)=(u-Row/2)2+(v-Col/2)20.5; factc=-D(u,v)2/(2*Do2); H(u,v)=exp(factc); G(u,v)=H(u,v)*FourIma(u,v); endendIFourIma=ifft2(G);%高斯低通滤波后的反傅立叶变换for ix=1:FRow %反傅立叶变换后*(-1)(x+y)并取实部 for iy=1:FCol LowPass(ix,iy)=(-1)(ix+iy)*IFourIma(ix,iy); end end LowPass=real(LowPass);(2)高斯低通滤波的主程序function main()%主程序ima=
17、imread(Test picture.jpg);imshow(ima);title(原图);ima0=im2double(ima);figure;ima1=GauseLowPass(ima, input(请输入d0值:);%选取Do=15,30,80,调用自定义高斯低通滤波函数imshow(ima1);%Do=15高斯低通滤波结果title(sprintf(%d%s,15,:高斯低通滤波结果);4结果与分析4.1选择一幅256级的灰度图像(1)在MTALAB文件夹中拷入个灰度图片,比如“Test picture.jpg”文件,在MATLAB命令输入窗中,调用函数info=imfinfo(Te
18、st picture.jpg);观查Test picture.jpg是否为灰度图像。在 ColorType: grayscale 属性行如此显示,如图4.1所示,则说明以Test picture.jpg命名的文件为灰度图像。图4.1 灰度图像的判别(2)设选择的图像的名称为”Test picture.jpg”,首先命令窗口输入f=imread(Test picture.jpg)读取图像文件,然后输入”whos f”,以显示图像的信息如图4.2所示:图4.2 图像信息从Class为uint8这一项信息可知图像的灰度级为256级。(2)另外我们还可以通过函数colorbar给一个坐标轴对象添加加一
19、条色带,该色带可以指示出该图像中不同颜色的数据值,以显示出图像的灰度级。输入以下代码:f = imread(Test picture.jpg);imshow(f);colorbar;可得到添加色带后的图像如图4.3所示,从图4.3可知该图像是数据类型为uint8的灰度图像,其灰度级范围从0255,级灰度级为256。图4.3 添加色带后的灰度图4.2 对灰度图像进行Fourier变换后的频谱图在MATLAB命令输入窗中,输入Fourier变换程序,按回车即生成Fourier变换后频谱图如图4.4所示:图4.4原图Fourier变换后频谱图4.3 二阶巴特沃斯(Butterworth)低通滤波结果
20、与分析在MATLAB命令输入窗中,输入巴特沃斯(Butterworth)低通滤波程序,按回车即生成原图如图4.5, Do=15的二阶Butterworth低通滤波结果如图4.6所示,Do=30的二阶Butterworth低通滤波结果如图4.7所示,Do=80的二阶Butterworth低通滤波结果如图4.8所示。图4.5 原图图4.6 Do=15的二阶Butterworth低通滤波结果图4.7 Do=30的二阶Butterworth低通滤波结果图4.8 Do=80的二阶Butterworth低通滤波结果结果分析: 在任何经BLPF处理过的图像中都没有明显的振铃效果,这是滤波器在低频和高频之间的
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 图像 频率 域低通 滤波 处理 程序设计 常用
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,个别因单元格分列造成显示页码不一将协商解决,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【二***】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【二***】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。