实验一傅立叶变换及图象的频域处理.doc

实验一傅立叶变换及图象的频域处理
2

———————————————————————————————— 作者：
———————————————————————————————— 日期：

个人收集 仅供参考学****勿做商业用途
实验一 傅立叶变换及图象的频域处理
一、实验目的
1、了解离散傅立叶变换的根本原理；
2、掌握应用MATLAB语言进展FFT及逆变换的方法；
3、了解图象在频域中处理方法，应用MATLAB语言作简单的低通滤波器。
二、实验原理
1、傅立叶变换的根本知识。
    在图象处理的广泛应用领域中，傅立叶变换起着非常重要的作用，具体表现在包括图象分析、图象增强及图象压缩等方面。
    假设f〔x,y〕是一个离散空间中的二维函数，那么该函数的二维傅立叶变换的定义如下：
　u=0,1…M-1　v=0,1…N-1　(1)
   离散傅立叶反变换的定义如下：
          x=0,1…M-1　y=0,1…N-1〔3〕
   F〔p,q〕称为f〔m,n〕的离散傅立叶变换系数。这个式子说明，函数f〔m,n〕可以用无数个不同频率的复指数信号和表示，而在频率〔w1，w2〕处的复指数信号的幅度和相位是F〔w1，w2〕。
例如，函数f〔m,n〕在一个矩形区域内函数值为1，而在其他区域为0.
假设f〔m,n〕为一个连续函数，那么f〔m,n〕的傅立叶变换的幅度值〔即〕显示为网格图。
将傅立叶变换的结果进展可视化的另一种方法是用图象的方式显示变换结果的对数幅值。
2、MATLAB提供的快速傅立叶变换函数
〔1〕fft2
fft2函数用于计算二维快速傅立叶变换，其语法格式为：
B = fft2(I)
B = fft2(I)返回图象I的二维fft变换矩阵，输入图象I和输出图象B大小一样。
例如，计算图象的二维傅立叶变换，并显示其幅值的结果，其命令格式如下
load imdemos saturn2
imshow(saturn2)
B = fftshift(fft2(saturn2));
imshow(log(abs(B)),[],'notruesize')
〔2〕fftshift
MATLAB提供的fftshift函数用于将变换后的图象频谱中心从矩阵的原点移到矩阵的中心，其语法格式为：
　　B = fftshift(I)
对于矩阵I，B = fftshift(I)将I的一、三象限和二、四象限进展互换。
〔2〕ifft2
ifft2函数用于计算图象的二维傅立叶反变换，其语法格式为：
B = ifft2(I)
B = ifft2(A)返回图象I的二维傅立叶反变换矩阵，输入图象I和输出图象B大小一样。其语法格式含义与fft2函数的语法格式一样，可以参考fft2函数的说明。
3、简单低通滤波器的设计
    一个图象经过傅立叶变换后，就从空域变到了频域，因此我们可以用信号处理中对于频域信号的处理方法对一幅图象进展处理。比方对图象进展低通滤波等。
    虽然在计算机中必定能够模拟一个锐截止频率的理想低通滤波器，但它们不能用电子元件来实现。实际中比拟常用的低通滤波器有：巴特沃思〔Butterworth〕滤波器、指数滤波器〔ELPF〕、梯形低通滤波器等。
4

个人收集 仅供参考学****勿做商业用途
    在实验中我们设计一个理想的低通滤波器。
    设计理想的低通滤波器由其定义可知只要设计一个与频域图象大小完全一样的矩阵。在某一个域值内该矩阵的值为1，其余为0即可。
例：假设图象的大小为128*128，那么可以这样设计一个低通滤波器：
H=zeros(128);
H(32:96,32:96)=1; %此处的范围是人为取定的，可以根据需要更改。
    假设图象矩阵I的傅立叶变换是B(已经用fftshift将频谱中心移至矩阵的中心)，那么对这幅图象做低通滤波，再做傅立叶逆变换命令为
LOWPASS=B.*H;　 %此处设变换后的矩阵为LOWPASS，另注意这儿是矩阵的点乘。
C=ifft2(LOWPASS