第7章 节 频域处理 数字图像处理 .ppt

第七章 频域处理 ;7.1 频域世界与频域变换;图7-2 正弦波的振幅A和相位φ ;图7-3 图7-1（a）波形的频域表示 (a) 幅频特性； (b) 相频特性 ;时域和频域之间的变换可用数学公式表示如下： ; 以上一维傅立叶变换可以很容易地推广到二维，如果二维函数f(x, y)满足狄里赫莱条件，则它的二维傅立叶变换对为 ;7.2.2 离散傅立叶变换 要在数字图像处理中应用傅立叶变换， 还需要解决两个问题：一是在数学中进行傅立叶变换的f(x)为连续（模拟）信号， 而计算机处理的是数字信号（图像数据）；二是数学上采用无穷大概念，而计算机只能进行有限次计算。通常， 将受这种限制的傅立叶变换称为离散傅立叶变换（Discrete Fourier Transform，DFT)。 设{f(x)|f(0), f(1), f(2), …, f(N-1)}为一维信号f(x)的N个抽样， 其离散傅立叶变换对为 ;(7-7) ;由欧拉公式可知 ;通常傅立叶变换为复数形式， 即 ; 通常称|F(u) |为f(x)的频谱或傅立叶幅度谱，φ(u)为f(x)的相位谱。 频谱的平方称为能量谱或功率谱，它表示为 ; 式中：u, x=0, 1, 2, …, M-1；v, y=0, 1, 2, …, N-1；x, y为时域变量，u, v为频域变量。 像一维离散傅立叶变换一样，系数1/MN可以在正变换或逆变换中，也可以在正变换和逆变换前分别乘以系数 ，只要两式系数的乘积等于1／MN即可。 二维离散函数的傅立叶频谱、 相位谱和能量谱分别为 ;7.2.3 离散傅立叶变换的性质 ;; 1. 可分离性 由可分离性可知，一个二维傅立叶变换可分解为两步进行， 其中每一步都是一个一维傅立叶变换。先对f(x, y)按行进行傅立叶变换得到F(x, v)，再对F(x, v)按列进行傅立叶变换，便可得到f(x, y)的傅立叶变换结果，如图7-4所示。显然对f(x, y)先按列进行离散傅立叶变换， 再按行进行离散傅立叶变换也是可行的。 ; 2. 平移性质 平移性质表明，只要将f(x, y)乘以因子(－1)x+y，再进行离散傅立叶变换，即可将图像的频谱原点（0，0）移动到图像中心（M／2, N／2）处。图7-5是简单方块图像平移的结果。 ; 3. 旋转不变性 由旋转不变性可知，如果时域中离散函数旋转θ0角度，则在变换域中该离??傅立叶变换函数也将旋转同样的角度。离散傅立叶变换的旋转不变性如图7-6所示。 ;7.2.4 快速离散傅立叶变换 离散傅立叶变换计算量非常大，运算时间长。可以证明其运算次数正比于N2，特别是当N较大时，其运算时间将迅速增长， 以至于无法容忍。为此，研究离散傅立叶变换的快速算法（Fast Fourier Transform， FFT）是非常有必要的。 下面介绍一种称为逐次加倍法的快速傅立叶变换算法（FFT），它是1965年Cooley和Tukey首先提出的。采用该FFT算法，其运算次数正比于NlbN，当N很大时计算量可以大大减少。例如，FFT的运算次数和DFT的运算次数之比，当N=1024时，比值为1／102.4；当N=4096时，比值可达1／341.3。 ; 由于二维离散傅立叶变换具有可分离性， 即它可由两次一维离散傅立叶变换计算得到，因此，仅研究一维离散傅立叶变换的快速算法即可。 先将式（7-7）写成 ; 这样，可将式(7-21)所示的一维离散傅立叶变换（DFT）用矩阵的形式表示为 ; 观察DFT的W阵，并结合W的定义表达式W=e-j2π／N，可以发现系数W是以N为周期的。这样，W阵中很多系数就是相同的， 不必进行多次重复计算，且由于W的对称性，即 ; 由W的周期性得：W4＝W0，W6＝W2，W9＝W1；再由W的对称性可得： W3＝－W1，W2＝－W0。于是式(7-23)可变为 ; 可见N=4的W阵中只需计算W0和W1两个系数即可。这说明W阵的系数有许多计算工作是重复的，如果把一个离散序列分解成若干短序列， 并充分利用旋转因子W的周期性和对称性来计算离散傅立叶变换，便可以简化运算过程，这就是FFT的基本思想。 设N为2的正整数次幂， 即 ;将式（7-26）代入式（7-21），离散傅立叶变换可改写成如下形式： ;于是式(7-28)变