《【精选】傅里叶变换及图像的频域处理》由会员分享，可在线阅读，更多相关《【精选】傅里叶变换及图像的频域处理（20页珍藏版）》请在金锄头文库上搜索。
1、 燕山大学课 程 设 计 说 明 书题目： 傅里叶变换及图像的频域处理 学院（系）： 里仁学院电气工程系 年级专业： 生物医学工程 10-1 班 学 号： 101203041007 学生姓名： 赵 林 静 指导教师： 孟 辉 教师职称： 讲 师 燕山大学课程设计（论文）任务书院（系）：里仁学院电气工程系 基层教学单位：生物医学工程 学 号 101203041007 学生姓名 赵林静 专业（班级） 生物医学工程 10-1 班设计题目 傅立叶变换及图像的频域处理设计技术参数在设计巴特沃思滤波器时满足的性能指标：通带截止频率 1000rad/s，通带最大衰减 3dB，阻带的截止频率 40000rad
2、/s 阻带的最小衰减 6dB设计要求理解离散傅立叶变换的基本原理；掌握应用 MATLAB 语言进行 FFT 及逆变换的方法；熟悉图像在频域中处理方法，应用 MATLAB 语言作简单的低通滤波器。工作量网上查阅相关资料图书馆查阅有关资料综合整理所学的知识完成相关程序工作计划第 1-2 天查找有关资料，开始准备；第 3-4 天编写程序，进行运行以及改进第 5 天整理程序，完成设计报告参考资料百度文库MATLAB 数据库论坛张德峰， MATLAB 数字图像处理 ，机械工业出版社。2009,1指导教师签字 孟 辉 基层教学单位主任签字 赵 勇说明：此表一式四份，学生、指导教师、基层教学单位、系部各一份
3、。2013 年 12 月 5 日 燕 山 大 学 课 程 设 计 说 明 书1目 录一 绪论 2二 设计目的4三 设计原理43.1 傅里叶变换的基本知识43.2MATLAB 提供的快速傅里叶变换63.3 简单低通滤波器的设计7四 设计要求104.1 理想低通滤波器的设计 114.2 理想高通滤波器的设计 13五 总结 17六 参考文献17燕 山 大 学 课 程 设 计 说 明 书2绪 论MATLAB 的英文全称 Matrix Laboratory(矩阵实验室).一开始它是一种专门用于矩阵数值计算的软件。 从这一点可以看出，它在矩阵运算方面有自己的特点。实际上，MATLAB 中的绝大多数运算都是
4、通过矩阵这一形式完成的。从理论上讲，图像是一种二维的连续函数，然而在计算机上对图像进行数字处理时，首先必须对其在空间和亮度上进行数字化。图像处理技术的发展大致经历了处创期，发展期，普及期和实用化期 4 个阶段。初创期开始于 20 世纪 60 年代，当时的图像采用像素型光栅进行扫描显示，大多采用中，大型机对其处理。 20 世纪 90 年代是图像处理技术的实用化时期，图像处理的信息量巨大，对处理器速度的要求极高。傅立叶变换是数字信号处理领域一种很重要的算法。要知道傅立叶变换算法的意义，首先要了解傅立叶理的意义。傅立叶原理表明：任何连续测量的时序或信号，都可以表示为不同频率的正弦波信号的无限叠加。而
5、根据该原理创立的傅立叶变换算法利用直接测量到的原始信号，以累加方式来计算该信号中不同正弦波信号的频率、振幅和相位。和傅立叶变换算燕 山 大 学 课 程 设 计 说 明 书3法对应的是反傅立叶变换算法。该反变换从本质上说也是一种累处理，这样就可以将单独改变的正弦波信号转换成一个信号。因此，可以说，傅立叶变换将原来难以处理的时域信号转换成了易于分析的频域信号（信号的频谱），可以利用一些工具对这些频域信号进行处理、加工。最后还可以利用傅立叶反变换将这些频域信号转换成时域信号。从现代数学的眼光来看，傅里叶变换是一种特殊的积分变换。它能将满足一定条件的某个函数表示成正弦基函数的线性组合或者积分。在不同的
6、研究领域，傅里叶变换具有多种不同的变体形式，如连续傅里叶变换和离散傅里叶变换。燕 山 大 学 课 程 设 计 说 明 书4傅立叶变换及图像的频域处理1、设计目的1、理解离散傅立叶变换的基本原理；2、掌握应用 MATLAB 语言进行 FFT 及逆变换的方法；3、熟悉图像在频域中处理方法，应用 MATLAB 语言作简单的低通滤波器。二、设计原理1、傅立叶变换的基本知识在图像处理的广泛应用领域中，傅立叶变换起着非常重要的作用，具体表现在包括图像分析、图像增强及图像压缩等方面。假设 f（x, y）是一个离散空间中的二维函数，则该函数的二维傅立叶变换的定义如下：p=0,1M-1 q=0,1N-1(2/)
7、(2/)0(,),MNjMpmjNqnmnFpqfe1 (6.1)或 p=0,1M-1 q=0,1N-1 121120(,)(,)MNjmjnmnfe（6.2）离散傅立叶反变换的定义如下：m=0,1M-1 1(2/)(2/)0(,),MNjMpmjNqnpqfnFen=0,1N-1 (6.3)F（p, q）称为 f（m, n）的离散傅立叶变换系数。这个式子表明，函数f（m, n）可以用无数个不同频率的复指数信号和表示，而在频率（w1，w2）处的复指数信号的幅度和相位是 F（w1，w2） 。燕 山 大 学 课 程 设 计 说 明 书5例如，函数 f（m, n）在一个矩形区域内函数值为 1，而在其
8、他区域为0，如图所示。了简便起见，假设 f（m, n）为一个连续函数，则 f（m , n）的傅立叶变换的幅度值（即 ）显示为网格图，如图所示。12(,)F将傅立叶变换的结果进行可视化的另一种方法是用图像的方式显示变换结果的对数幅值 ，如图所示。12log(,)几种简单函数的傅立叶变换的频谱可以直观的表示为图所示的样子。燕 山 大 学 课 程 设 计 说 明 书62、MATLAB 提供的快速傅立叶变换函数（1） fft2fft2 函数用于计算二维快速傅立叶变换，其语法格式为：B = fft2(I)B = fft2(I)返回图像 I 的二维 fft 变换矩阵，输入图像 I 和输出图像 B 大小相同
9、。例如，计算图像的二维傅立叶变换，并显示其幅值的结果，如图所示，其命令格式如下load imdemos saturn2imshow(saturn2)B = fftshift(fft2(saturn2);imshow(log(abs(B), , notruesize)燕 山 大 学 课 程 设 计 说 明 书7（2） fftshiftMATLAB 提供的 fftshift 函数用于将变换后的图像频谱中心从矩阵的原点移到矩阵的中心，其语法格式为：B = fftshift(I)对于矩阵 I，B = fftshift(I)将 I 的一、三象限和二、四象限进行互换。（3） ifft2ifft2 函数用于
10、计算图像的二维傅立叶反变换，其语法格式为：B = ifftn(I)B = ifftn(I)返回图像 I 的二维傅立叶反变换矩阵，输入图像 I 和输出图像 B 大小相同。其语法格式含义与 fft2 函数的语法格式相同，可以参考 fft2 函数的说明。3、简单低通滤波器的设计一个图像经过傅立叶变换后，就从空域变到了频域，因此我们可以用信号处理中对于频域信号的处理方法对一幅图像进行处理。比如对图像进行低通滤波等。一个二维的理想低通滤波器（ILPF） ，它的传递函数由下式确定：0),(01),( DvuvuH若若（6.4）式中 D0 是一个规定的非负的量，称为截止频率，虽然在计算机中必定能够模拟一个锐
11、截止频率的理想低通滤波器，但它们不能用电子元件来实现。实际中比较常用的低通滤波器有：巴特沃思（Butterworth）滤波器、指数滤波器（ELPF） 、梯形低通滤波器等。在实验中我们设计一个理想的低通滤波器。设计理想的低通滤波器由其定义可知只要设计一个与频域图像大小完全燕 山 大 学 课 程 设 计 说 明 书8相同的矩阵。在某一个域值内该矩阵的值为 1，其余为 0 即可。例：若图像的大小为 128*128，则可以这样设计一个低通滤波器：H=zeros(128);H(32:96, 32:96)=1; %此处的范围是人为取定的，可以根据需要更改。若图像矩阵 I 的傅立叶变换是 B(已经用 fft
12、shift 将频谱中心移至矩阵的中心)，则对这幅图像做低通滤波，再做傅立叶逆变换命令为LOWPASS=B.* H; %此处变换后的矩阵为 LOWPASS，另注意这儿是矩阵的点乘。C=ifft2(LOWPASS);imshow(abs(C)参考代码实现：I=imread(guzhe.jpg); %读入原图像文件imshow(I); %显示原图像fftI=fft2(I); %二维离散傅立叶变换sfftI=fftshift(fftI); %直流分量移到频谱中心RR=real(sfftI); %取傅立叶变换的实部II=imag(sfftI); %取傅立叶变换的虚部A=sqrt(RR.2+II.2);%
13、计算频谱幅值A=（A-min(min(A)）/(max(max(A)-min(min(A)*225;%归一化figure; %设定窗口imshow(A); %显示原图像的频谱2理想低通滤波器参考代码实现：I = imread( guzhe.jpg);f1,f2 = freqspace(size(I),meshgrid); Hd = ones(size(I); r = sqrt(f1.2 + f2.2); %0.1Hd(r0.1) = 0; Y=fft2(double(I); Y=fftshift(Y); Ya=Y.*Hd;Ya=ifftshift(Ya); Ia01=ifft2(Ya);%0.
14、2Hd(r0.2) = 0; 燕 山 大 学 课 程 设 计 说 明 书9Y=fft2(double(I); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia02=ifft2(Ya);%0.5Hd(r0.5) = 0; Y=fft2(double(I); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia05=ifft2(Ya);subplot(2,2,1),imshow(I),title(原图像)subplot(2,2,2),imshow(uint8(Ia01),title(r=0.1)subplot(2,2,3),imshow(uint8(Ia02),title(r=0.2)subplot(2,2,4),imshow(uint8(Ia05),title(r=0.5)3. 理想高通滤波器参考代
