第7章-频域处理幻灯片.ppt

多媒体通信 北京科技大学 杨 扬 第7章 频域处理 频域世界与频域变换 傅立叶变换 频域变换的一般表达式 离散余弦变换 离散沃尔什哈达玛变换 小波变换简介 1、频域世界与频域变换 频域变换的理论基础是： “任意波形都可以用单纯的正弦波的和来表示” 因此图像的频域变换为： （1）将图像看成是线性叠加系统 ； （2）图像在空域上相关性很强； （3）图像变换是将图像从空域变换到其它域如频域的数学变换。 常用的变换有：傅立叶变换、离散余弦变换、小波变换、离散K-L变换 1、频域世界与频域变换 1、频域世界与频域变换 2、傅立叶变换 傅立叶变换是一种常用的正交变换，它的理论完善,应用程序多。 在图像处理应用领域，傅立叶变换起着非常重要的作用，可用它完成图像分析、图像增强以及图像压缩等工作。 当一个一维信号f(x)满足狄里赫莱条件，即f(x) （1） 具有有限个间断点；  （2） 具有有限个极值点；  （3） 绝对可积。 则其傅立叶变换对（傅立叶变换和逆变换）一定存在。 2.1 连续函数的傅立叶变换 2.1 连续函数的傅立叶变换 2.1 连续函数的傅立叶变换 若把一个一维输入信号作一维傅立叶变换，该信号就被变换到频域上的一个信号，即得到了构成该输入信号的频谱，频谱反映了该输入信号由哪些频率构成。 一维傅立叶变换对的定义为： 2.2 离散傅立叶变换 设{f(x)|f(0), f(1), f(2), …, f(N-1)}为一维信号f(x)的N个抽样， 其离散傅立叶变换（Discrete Fourier Transform，DFT)对为： 式中：x，u=0，1，2，…，N－1。 2.3 离散傅立叶变换的性质 2.3 离散傅立叶变换的性质 2.3 离散傅立叶变换的性质-平移性质 2.3 离散傅立叶变换的性质-可分离性 2.3 离散傅立叶变换的性质-旋转不变性 2.4 快速离散傅立叶变换 离散傅立叶变换计算量非常大，运算时间长。可以证明其运算次数正比于N2，特别是当N较大时，其运算时间将迅速增长，以至于无法容忍。为此，研究离散傅立叶变换的快速算法(Fast Fourier Transform，FFT)是非常有必要的。 2.4 快速离散傅立叶变换 一维离散傅立叶变换的快速算法如下式： 式中，W=e-j2π/N ，称为旋转因子。 2.4 快速离散傅立叶变换 前面所示的一维离散傅立叶变换（DFT）用矩阵的形式表示为： 式中，由Wux构成的矩阵称为W阵或系数矩阵。 2.4 快速离散傅立叶变换 从上式中的W阵，并结合W的定义，可以发现W是以N为周期的。这样，W阵中很多系数就是相同的，不必进行多次重复计算。同时从下面关系可以看出W的对称性： 对于N=4，W阵为： 2.4 快速离散傅立叶变换 由W的周期性得：W4＝W0，W6＝W2，W9＝W1；再由W的对称性可得： W3＝－W1，W2＝－W0。于是上式可变为： 2.4 快速离散傅立叶变换 可见N=4的W阵中只需计算W0和W1两个系数即可。这说明W阵的系数有许多计算工作是重复的，如果把一个离散序列分解成若干短序列， 并充分利用旋转因子W的周期性和对称性来计算离散傅立叶变换，便可以简化运算过程，这就是FFT的基本思想。 设N为2的正整数次幂， 即： 2.4 快速离散傅立叶变换 2.4 快速离散傅立叶变换 2.4 快速离散傅立叶变换 在此，以计算N=8的DFT为例，此时n=3，M=4。 2.4 快速离散傅立叶变换 2.4 快速离散傅立叶变换 2.4 快速离散傅立叶变换 2.4 快速离散傅立叶变换 2.4 快速离散傅立叶变换 3、频域变换的一般表达式--可分离变换 二维傅立叶变换可以利用变换核的可分离性， 用两次一维变换来实现，即可先对f(x, y)的每一行进行一维变换得到F(x, v)，再沿F(x, v)每一列取一维变换得到变换结果F(u, v)。对于其他的图像变换，只要其变换核是可分离的，同样也可用两次一维变换来实现。 如果先对f(x, y)的每一列进行一维变换得到F(y, u)，再沿F(y, u)每一行取一维变换得到F(u, v)，其最终结果是一样的。该结论对反变换核也适用。 3、频域变换的一般表达式--图像变换的矩阵表示 数字图像都是实数矩阵， 设f(x, y)为M×N的图像灰度矩阵，通常为了分析、推导方便，可将可分离变换写成矩阵的形式： F=PfQ F=P-1FQ-1 其中，F、f是二维M×N的矩阵；P是M×M矩阵；Q是N×N矩阵。 4、离散余弦变换（DCT） 离散余弦变