灰度图像频域处理信号课程设计.doc

信号课程设计题 目: 灰度图像频域处理1．理论目的 《信号与系统》课程设计的目的之一是为了巩固课堂理论学习，并能用所学理论知识正确分析信号处理的基本问题和解释信号处理的基本现象。2．实践目的 正确理解二维傅里叶变换及滤波的基本概念，掌握低通、高通滤波器。3、内容和要求 提供给学生一幅灰度图像（house.bmp，下载地址 HYPERLINK /s/1eQ7TXo2 /s/1eQ7TXo2 ）。使用Matlab中的imread命令读取图像，并用imshow显示图像。将图像数据变换到二维频域（采用命令fft2），判断该图像的截止频率。在频域进行低通滤波和高通滤波，恢复空域结果（采用命令ifft2），比较滤波前后的图像差异。设计一个差分滤波器，得到对该图像的二维一阶差分结果。Matlab命令：imread, imshow, fft2, fftshift一、摘要本次课程设计主要是利用MATLAB－图像处理工具箱，对二维信号（数字图像）进行快速傅立叶变换与反变换，并且设计差分滤波器、低通滤波器和高通滤波器对二维信号（数字图像）进行滤波。关键词：图像处理工具箱、二维数字图像处理、低通滤波、高通滤波、差分滤波二维信号采集1.1图像的选取从所给网址下载图像并存于D盘中命名为house，得到原图如下图 SEQ 图 \* ARABIC 1 选取的需要处理的原图1.2 图像的读取MATLAB是解释型语言，就是说MATLAB命令行中敲入的命令在当前MATLAB进程中被解释运行。利用函数imread()可完成图形图像文件的读取，语法：MATLAB程序：imshow(f); f=imread(D:\house.bmp);imshow(f);title(原图);图 2 MATLAB中读取图像的效果图2、将图像数据变换到二维频域MATLAB为计算数据的离散快速傅立叶变换，提供了一系列丰富的数学函数，主要有Fft、Ifft，Fft2、Ifft2，Fftn、Ifftn和Fftshift、Ifftshift等。当所处理的数据的长度为2的幂次时，采用基-2算法进行计算，计算速度会显著增加。所以，要尽可能使所要处理的数据长度为2的幂次或者用添零的方法来添补数据使之成为2的幂次。2.1利用命令fft2判断截止频率MATLAB程序：f=imread(D:\house.bmp);imshow(f);C= fft2(double(f));D=fftshift(fft2(double(f))); E=ifft2(D); figure,imshow(log(abs(D)+1),[]); title(变换频谱图); 图3 直接变换频谱图从图像可以看出，图中两竖线与X轴线的交点就是截止频率。3、低通滤波处理3.1 低通滤波器介绍 低通滤波器是容许低于截至频率的信号通过， 但高于截止频率的信号不能通过的电子滤波装置。对于不同滤波器而言，每个频率的信号的减弱程度不同。当使用在音频应用时，它有时被称为高频剪切滤波器, 或高音消除滤波器。　　低通滤波器概念有许多不同的形式，其中包括电子线路、图像模糊处理等等，这两个工具都通过剔除短期波动、保留长期发展趋势提供了信号的平滑形式。低通滤波器有很多种，其中，最通用的就是巴特沃斯滤波器和切比雪夫滤波器。3.1.1 巴特沃斯滤波器 巴特沃斯滤波器，又被称为“最平”的幅频响应滤波器。这是因为，该滤波器在通带内具有最大平坦的幅度特性，而且随着频率升高呈现出出单调减小的特点。巴特沃斯滤波器的特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。 在振幅的对数对角频率的波得图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大。一阶巴特沃斯滤波器的衰减率为每倍频6分贝，每十倍频20分贝。二阶巴特沃斯滤波器的衰减率为每倍频12分贝、 三阶巴特沃斯滤波器的衰减率为每倍频18分贝、如此类推。巴特沃斯滤波器的振幅对角频率单调下降，并且也是唯一的无论阶数，振幅对角频率曲线都保持同样的形状的滤波器。只不过滤波器阶数越高，在阻频带振幅衰减速度越快。其他滤波器高阶的振幅对角频率图和低级数的振幅对角频率有不同的形状。3.2 二维信号（数字图像）的低通滤波处理先对二维信号（数字图像）进行FFT，然后进行低通滤波（即去掉高频分量），然后再反变换，得到滤波图像。M=imread(D:\house.bmp); %原图figure,imshow(M)title(原图)F=fft2(M); %对原图像进行傅里叶变换F=fftshift(F);figure;imshow(log(abs(F)),[]);