四、图像的频域变换——傅立叶变换

约瑟夫·傅立叶(Joseph Fourier)——法国数学家、物理学家,1807年提出傅立叶变换。

        傅立叶变换是最早研究与应用的酉变换;60年代出现快速傅立叶变换;傅立叶变换域也称为频域。

基本数学概念

        调谐信号(欧拉公式):

        傅立叶积分:

傅立叶变换的定义

一维连续

        f(x)为连续可积函数,其傅立叶变换定义为:

其反变换为:

通常f(x)的傅立叶变换为复数,可有通用表示式为:,、分别称为傅立叶变换的实部和虚部。

        可进一步写为指数形式:

其中:称之为的幅度谱、振幅谱或傅立叶谱;称之为的相位谱、相位角。

图2.1 将一个信号的波形分解为许多不同频率正弦波之和
图2.2 方波信号的分解
图2.3 方波信号的时、频域描述
图2.4 一维傅立叶变换示例
图2.5 矩形函数的傅立叶变换
图2.6 sin(x)/x类函数的傅立叶变换
图2.7 常数函数的傅立叶变换
图2.8 脉冲函数的傅立叶变换
图2.9 余弦函数的傅立叶变换(图中打错字了)
图2.10 几种特殊函数的傅立叶变换

一维离散

        一维离散傅立叶变换公式为:

逆变换为:

逆变换的另一种表达形式:

二维连续

        二维傅立叶变换由一维傅立叶变换推广而来:

逆变换:

幅度谱:

相位谱:

图2.11 傅立叶变换图例

二维离散

        对于二维傅立叶变换,其离散形式为:

逆变换为:

幅频谱、相位谱:

二维离散傅立叶变换的性质

        1)线性性质(加法定理):

图2.12 加法性质示例

        2)比例性质(相似性定理):

图2.13 比例性质示例

比例性质表明:信号在时域中压缩(k>1,变化速度加快)等效于在频域扩展(频带加宽);反之亦然。

        3)可分离性:

图2.14 可分离性示例

二维DFT可分离为两次一维DFT。

        4)空间位移(位移定理):

图2.15 空间位移示例

空间位移特性表明:信号在时域中沿时间轴平移一个常数时,等效于频谱函数的相位谱改变,而幅度谱不变。

        5)频率位移:

函数的频率位移相当于傅立叶变换的坐标原点平移,而幅度谱和相位谱不变。

        6)周期性:

图2.16 二维离散DFT的周期性扩展

离散傅立叶变换DFT和它的逆变换是以N为周期的函数。

        7)共轭对称性:若f(x,y)为实函数,F(u,v)为其傅立叶变换,则

图像的傅立叶变换结果是以原点为中心的共轭对称函数。

        8)旋转不变性:

图2.17 旋转特性示例

旋转特性描述:如果f(x,y)旋转了一个角度α,那么f(x,y)旋转后图像的傅立叶变换也旋转了相同的角度α。

结论:对图像的旋转变换和傅立叶变换的顺序是可交换的。

        9)平均值:

离散函数的均值等于该函数傅立叶变换在(0,0)点的值。

        10)卷积定理:空域中的卷积等价于频域中的相乘。

        11)相关定理:空域中f(x,y)与g(x,y)的相关等价于频域中F(u,v)的共轭与G(u,v)相乘。

互相关:

自相关:

        12)拉普拉斯函数:

其傅立叶变换为:

这个定理将在图像的边界提取中用到。

二维离散傅立叶变换的显示与计算

离散傅立叶变换的显示

        按照标准的傅立叶变换公式,其幅度谱的强度分布具有下列特性:

图3.1 幅度谱强度分布

        在光学傅立叶变换中,人们已习惯于变化领域中的低谱部分位于中央。使频域的频谱分布中间低、周围高,有利于对频谱的解释和进行各种计算与分析。

图3.2 低谱部分位于中央的频谱示例

为了达到上述要求——图像中心化,借助于傅立叶变换的周期性与频率位移性质,对频域进行换位:

        使频域的中心位移:

        相当于对原始图像f(x,y)乘以,再进行傅立叶变换:

        对应于的反变换不等于f(x,y):

图3.3 图像中心化示例

        二维傅立叶变换域分布特性:

图3.4 二维傅立叶变换域分布特性

离散傅立叶变换的幅度与相位

        图像信号的傅立叶变换包含幅度与相位两部分;幅度谱具有较明显的信号结构特征和易于解释;实验证明,幅度本身只包含有图像本身含有的周期结构,并不表示其在何处;相位谱类似随机图案,一般难以进行解释;物体在空间的移动,相当于频域的相位移动,相位谱具有同样重要的意义。

        单凭幅度或相位信息,均不足以恢复原图像。 

离散傅立叶变换的计算

图3.5 离散傅立叶变换示例-1
图3.6 离散傅立叶变换示例-2

快速傅立叶变换(FFT)原理

基本思想

        快速傅立叶变换的基本思想就是分解-征服,即将大的问题分解成诸多小问题,再一一解决这些小问题,从而最终解决大问题。

图3.7 快速傅立叶变换基本思想示意

        1)将变换公式分解为奇数项和偶数项之和。令:

DFT可表为:

令:N=2M

由于:

可得到:

进一步分析:

还可以得到:

逆向FFT算法

        算法思想:用正向变换计算逆向变换。

        设,可有:

即:对F(u)取共轭,利用正向FFT进行变换计算,其结果取共轭后再乘以N,即可得到f(x)。

二维快速傅立叶变换

        利用傅立叶变换的分离性质,对二维FFT进行2次的一维FFT变换:

你可能感兴趣的:(四、图像的频域变换——傅立叶变换)