图像频率域低通滤波处理程序设计

1、专业综合课程设计任务书学生姓名： 陈德松 专业班级： 电信 0901班 指导教师： 黄朝兵 工作单位： 信息工程学院 题 目：图像频率域低通滤波处理程序设计初始条件： （1）提供实验机房及其matlab软件； （2）数字图像处理的基本理论学习。要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）（1）掌握频率域处理的基本原理，利用matlab设计程序完成以下功能； （2）选择一幅256级的灰度图像； （3）对该图像作Fourier变换，得到其频谱图像； （4）按照二阶Butterworth低通滤波器的表达式设计传递函数，对该频谱图像进行低通滤波，对滤波后的频谱图像作

2、逆Fourier变换得到滤波结果，显示结果图； （5）按照高斯低通滤波器的表达式设计传递函数，对该频谱图像进行低通滤波，对滤波后的频谱图像作逆Fourier变换得到滤波结果，显示结果图； （6）对二种滤波器的滤波结果进行分析比较； （7）要求阅读相关参考文献不少于5篇； （8）根据课程设计有关规范，按时、独立完成课程设计说明书。时间安排：(1) 布置课程设计任务,查阅资料,确定方案四天； (2) 进行编程设计一周； (3) 完成课程设计报告书三天；指导教师签名： 年 月 日系主任（或责任教师）签名： 年 月 日目录摘 要I1 MATLAB的简介12 原理与实现22.1频率域增强基本理论22.2

3、 傅立叶变换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 高斯(Gaussian)低通滤波结果与分析174.5两种滤波器的滤波结果的比较205 心得体会2

4、1参考文献22摘 要数字图像处理的目的之一是改善图像的视觉效果, 以便人眼或机器对图像进一步理解。而图像增强技术通过采用某种技术手段可以改善图像的视觉效果。其中一种技术是频率域增强技术，它是增强技术的重要组成部分。本设计重点介绍了用Matlab对一幅256级的灰度图像作Fourier变换，然后对变换后得到的频谱图像进行分别进行二阶Butterworth低通滤波和高斯低通滤波，然后对滤波后的频谱图像作逆Fourier变换得到滤波结果，显示结果图。本程序的设计很简洁，加上注释显得通俗易懂，而且本程序中很多地方都有输入提示，有一定的人机交互。关键词：Matlab ，Fourier变换，Butterw

5、orth，高斯(Gaussian)，低通滤波1 MATLAB的简介MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首

6、屈一指。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)，这样具有所需特性的 g(x, y) 就可算出 G(u, v) 而得到： 2

8、.2 傅立叶变换傅里叶变换是将时域信号分解为不同频率的正弦信号或余弦函数叠加之和。傅立叶变换是数字图像处理技术的基础，其通过在时空域和频率域来回切换图像，对图像的信息特征进行提取和分析，简化了计算工作量，被喻为描述图像信息的第二种语言，广泛应用于图像变换，图像编码与压缩，图像分割，图像重建等。因此，对涉及数字图像处理的工作者，深入研究和掌握傅立叶变换及其扩展形式的特性，是很有价值得。把傅立叶变换的理论通其物理解释相结合，将有助于解决大多数图像处理问题。傅里叶变换可分为连续傅里叶变换、离散傅里叶变换、快速傅里叶变换。快速傅里叶变换(FFT)是计算离散傅里叶变换(DFT)的快速算法。离散傅里叶变换

9、运算量巨大，计算时间长，即运算时间很长。而快速傅里叶变换的提出将傅里叶变换的复杂度由降到了，很大程度上减少了计算量。 令，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 对灰度图像进行Fourier变换的程序I=imread('Test p

12、icture.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 picture.jpg')

13、; 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滤波函

14、数 s(i,j)=h*s(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单独运行

15、。（1）高斯低通滤波的功能子程序GauseLowPass()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的高斯

16、低通滤波器 for v=1:FCol D(u,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）高斯低通滤波的主程序

17、function main()%主程序ima=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)在M

18、TALAB文件夹中拷入个灰度图片，比如“Test picture.jpg”文件,在MATLAB命令输入窗中，调用函数info=imfinfo('Test picture.jpg');观查Test picture.jpg是否为灰度图像。在 ColorType: 'grayscale' 属性行如此显示，如图4.1所示，则说明以Test picture.jpg命名的文件为灰度图像。图4.1 灰度图像的判别(2)设选择的图像的名称为”Test picture.jpg”,首先命令窗口输入f=imread('Test picture.jpg')读取图像文件，

19、然后输入”whos f”，以显示图像的信息如图4.2所示：图4.2 图像信息从Class为uint8这一项信息可知图像的灰度级为256级。(2)另外我们还可以通过函数colorbar给一个坐标轴对象添加加一条色带，该色带可以指示出该图像中不同颜色的数据值，以显示出图像的灰度级。输入以下代码：f = imread('Test picture.jpg');imshow(f);colorbar;可得到添加色带后的图像如图4.3所示，从图4.3可知该图像是数据类型为uint8的灰度图像，其灰度级范围从0255，级灰度级为256。图4.3 添加色带后的灰度图4.2 对灰度图像进行Four

20、ier变换后的频谱图在MATLAB命令输入窗中，输入Fourier变换程序，按回车即生成Fourier变换后频谱图如图4.4所示：图4.4原图Fourier变换后频谱图4.3 二阶巴特沃斯（Butterworth）低通滤波结果与分析在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的二阶Butterwort

21、h低通滤波结果图4.7 Do=30的二阶Butterworth低通滤波结果图4.8 Do=80的二阶Butterworth低通滤波结果结果分析： 在任何经BLPF处理过的图像中都没有明显的振铃效果，这是滤波器在低频和高频之间的平滑过渡的结果。低通滤波是一个以牺牲图像清晰度为代价来减少噪声干扰效果的修饰过程。通过与原图相比较，从随着D0值分别取15,30,80滤波后得到的图像可以看出，D0的值越大图像越清晰。巴特沃斯低通滤波器的优点是：一、模糊大大减少。因为包含了许多高频分量；二、没有振铃现象。因为滤波器是平滑连续的。4.4 高斯(Gaussian)低通滤波结果与分析(1) 在MATLAB命令输

22、入窗中，输入高斯(Gaussian)低通滤波程序，按回车即生成原图如图4.9，Do=15的高斯低通滤波结果如图4.10所示，Do=30的高斯低通滤波结果如图4.11所示，Do=80的高斯低通滤波结果如图4.12所示图4.9 原图图4.10 Do=15的高斯低通滤波结果图4.11 Do=30的高斯低通滤波结果图4.12 Do=80的高斯低通滤波结果结果分析： 通过与原图相比较，从随着D0值分别取15,30,80滤波后得到的图像可以看出，D0的值越大图像越清晰，高斯低通滤波器中，D0越大，所加的窗就越窄，那么说D0越大，所通过的频带是较低的，所滤得的图像越柔和平滑。4.5两种滤波器的滤波结果的比较

23、（1）巴特沃斯低通滤波器在任何经巴特沃斯低通滤波处理过的图像中都没有明显的振铃效果，这是滤波器在低频和高频之间的平滑过渡的结果。低通滤波是一个以牺牲图像清晰度为代价来减少噪声干扰效果的修饰过程。 巴特沃斯低通滤波器的优点是：1、模糊大大减少。因为包含了许多高频分量；2、没有振铃现象。因为滤波器是平滑连续的。（2）高斯低通滤波器高斯低通滤波后得到的图像，高斯低通滤波器没有振铃现象，但与阶数为2的巴特沃斯低通滤波器相比，其通带要宽些，这样对应的空间滤波器的灰度级轮廓更窄些，因而平滑效果要差些。（3）巴特沃思滤波器和高斯滤波器滤去的频率和通过的频率之间没有明显的不连续性，图像的模糊程度降低，而且也没

24、有振铃效应，这是由于在低频和高频之间，滤波器平滑过渡的缘故。5 心得体会 经过两周的努力，总算把专业综合课程设计做完了。 通过该课程设计，全面系统的理解了数字图像处理的一般原理和基本实现方法。把死板的课本知识变得生动有趣，激发了学习的积极性。把学过的数字图像处理基础原理的知识强化，能够把课堂上学的知识通过自己编写的程序表示出来，加深了对理论知识的理解。 在这次课程设计中，我先是认真阅读课本上的相关知识，理解透后又翻阅关于Matlab的书籍，学习Matlab中一些函数及运算符的用法。总体来说,这次课设我学到了很多。在设计过程中,加深了对可内知识的理解就,真正懂得了学以致用，熟悉了Matlab的使用，了解了Matlab在数字图像处理中的重大应用。做课程设计我体会到了设计的艰辛的同时,更让我体会