数字图像处理2-赵建峰

数字图像处理

(DigitalImageProcessing)内蒙古科技大学信息学院信息处理研究室赵建峰第二章、图像变换理论§2.1空域、变换域基础知识§2.2空域、变换域变换方法§2.1空域、变换域基础知识一、变化的必要性二、变换基础知识三、常用变换简介一、变换的必要性我们人类视觉所感受到的是在空间域和时间域的信号。许多问题在频域中讨论时，有其非常方便分析的一面。1、空间位置上的变化不改变信号的频域特性。2、利用频率成分和图像外表之间的对应关系，使在空间域表述困难的增强任务在频率域中执行。3、滤波在频率域更为直观，可以解释空间域滤波的某些性质。二、变换基础知识空域和频域之间的变换关系：采用复数表示法，可同时表示信号的振幅和相位。故上式可表示为：这种变换，一般为线性正交变换。以傅立叶变换为例：任意信号可分解为正弦波的加权和正弦波的振幅A和相位φ信号的频域表示

(a)幅频特性；(b)相频特性三、常用变换简介数学变换是将原定义在空间的图像以某种形式转换到另外一些空间，并利用输入图像在这些空间的特有性质有效而快速地对图像进行处理和分析。除了傅立叶变换外，常用的变换还有Gabor变换、小波变换、离散余弦变换、PCA变换等。这些变换在图像分析、滤波、增强、压缩等处理中都有着非常典型而重要的应用。§2.2空域、变换域变换方法一、频域变换的一般形式二、傅立叶变换三、离散余弦变换（DCT）四、离散沃尔什-哈达玛变换（WHT）五、小波变换一、频域变换的一般形式二维频域变换的一般形式为：上式中，1为正变换，2为逆变换。式中g(x,y,u,v)和h(x,y,u,v)分别为正、逆变换核，x,u=0,1,2,…,M－1；y,v=0,1,2,…,N－1。二维变换的特性如果：g(x,y,u,v)=g1(x,u)g2(y,v)h(x,y,u,v)=h1(x,u)h2(y,v)则称正、反变换核是可分离的。如果：g1(x,u)=g2(y,v)h1(x,u)=h2(y,v)则称该变换核是对称的。可分离性的应用利用变换核的可分离性，可用两次一维变换来实现二维变换。1、先对f(x,y)的每一行进行一维变换得到F(x,v)，再沿F(x,v)每一列取一维变换得到变换结果F(u,v)。2、先对f(x,y)的每一列进行一维变换得到F(y,u)，再沿F(y,u)每一行取一维变换得到变换结果F(u,v)。二者变换的最终结果是一样的。变换的矩阵表示数字图像都是实数矩阵，设f(x,y)为M×N的图像灰度矩阵，为了分析、推导方便，可将可分离变换写成矩阵的形式：F=PfQF=P-1FQ-1

其中，F、f是二维M×N的矩阵；P是M×M矩阵；Q是N×N矩阵。式中，u=0,1,2,…,M－1，v=0,1,2,…,N－1。二、傅立叶变换一、连续函数傅立叶变换二、离散傅立叶变换三、常用变换简介连续函数傅立叶变换对一个一维信号进行傅立叶变换，即得到构成该信号的频谱，其频谱反映了该输入信号由哪些频率构成。一个可进行傅立叶变换的一维信号f(x)需满足狄里赫莱条件，即f(x)：1、具有有限个间断点；2、具有有限个极值点；3、绝对可积。连续函数傅立叶变换定义式中： ，x称为时域变量，u称为频域变量。式中：x,y为时域变量；u,v为频域变量。一维：二维：离散傅立叶变换因为数字图像为离散信号，对其进行的傅立叶变换称为离散傅立叶变换(DiscreteFourierTransform，DFT)。设{f(x)|f(0),f(1),f(2),…,f(N-1)}为一维信号f(x)的N个抽样，其离散傅立叶变换对为：一维离散傅立叶变换设{f(x)|f(0),f(1),f(2),…,f(N-1)}为一维信号f(x)的N个抽样，其离散傅立叶变换对为：式中：x，u=0，1，2，…,N－1，其中系数1/N可以在正变换或逆变换中，也可在傅立叶正变换和逆变换前分别乘以 ，这是无关紧要的，只要正变换和逆变换前系数乘积等于1/N即可。复数表示法由欧拉公式可得：可见，离散序列的傅立叶变换仍是一个离散的序列，每一个u对应的傅立叶变换结果是所有输入序列f(x)的加权和（每一个f(x)都乘以不同频率的正弦和余弦值），u决定了每个傅立叶变换结果的频率。通常傅立叶变换为复数形式：式中，R(u)和I(u)分别是F(u)的实部和虚部。上式也可表示成指数形式：F(u)=|F(u)|

ejφ(u)

通常称|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即可。二维离散函数的傅立叶频谱、相位谱和能量谱分别为：式中，R(u,v)和I(u,v)分别是F(u,v)的实部和虚部。二维傅立叶变换的性质一、可分离性二、周期与共轭对称三、平移性四、旋转特性五、线性与相似性六、均值性七、拉普拉斯八、卷积与相关二维傅立叶变换的矩阵形式二维傅立叶变换的快速算法由于二维离散傅立叶变换具有可分离性，即它可由两次一维离散傅立叶变换计算得到，三、离散余弦变换（DCT）一、一维离散余弦变换二、二维离散余弦变换三、快速离散余弦变换离散余弦变换（DiscreteCosineTransform，DCT）的变换核为余弦函数。DCT除了具有一般的正交变换性质外，它的变换阵的基向量能很好地描述人类语音信号和图像信号的相关特征。因此，在对语音信号、图像信号的变换中，DCT变换被认为是一种准最佳变换。近年颁布的一系列视频压缩编码的国际标准建议中，都把DCT作为其中的一个基本处理模块。除此之外，DCT还是一种可分离的变换。一维离散余弦变换一维DCT的变换核定义为：式中，x,u=0,1,2,…,N－1；一维DCT定义如下：设{f(x)|x=0,1,…,N-1}为离散的信号列。式中，u,x=0,1,2,…,N－1。一维离散余弦变换的矩阵形式将变换式展开整理后，可以写成矩阵的形式，即F=Gf。一维DCT的逆变换IDCT定义为式中，x,u=0,1,2,…,N－1。可见一维DCT的逆变换核与正变换核是相同的。二维离散余弦变换二维离散余弦变换定义如下：设f(x,y)为M×N的数字图像矩阵，则：式中，C(u)和C(v)的定义同式（7-48）；x,u=0,1,2,…,M－1；y,v=0,1,2,…,N－1。式中：x,u=0,1,2,…,M－1；y,v=0,1,2,…,N－1。其中正变换核为：二维DCT逆变换定义如下：式中：x,u=0,1,2,…,M－1；y,v=0,1,2,…,N－1。二维DCT的矩阵形式如下：F=GfGT二维DCT的逆变换核与正变换核相同，且是可分离的，即：根据可分离性，二维DCT可用两次一维DCT来完成，其算法流程与DFT类似，即：快速离散余弦变换离散余弦变换(DCTforDiscreteCosineTransform)是与傅里叶变换相关的一种变换，它类似于离散傅里叶变换(DFTforDiscreteFourierTransform),但是只使用实数。离散余弦变换相当于一个长度大概是它两倍的离散傅里叶变换，这个离散傅里叶变换是对一个实偶函数进行的（因为一个实偶函数的傅里叶变换仍然是一个实偶函数），在有些变形里面需要将输入或者输出的位置移动半个单位(DCT有8种标准类型，其中4种是常见的)。四、离散沃尔什-哈达玛变换（WHT）一、沃尔什函数二、一维离散沃尔什-哈达玛变换三、二维离散沃尔什-哈达玛变换一、沃尔什函数沃尔什函数系是一个完备正交函数系，其值只能取＋1和－1，由美国数学家沃尔什在1923年提出。从排列次序上可将沃尔什函数分为三种定义方法：第一种是按照沃尔什排列来定义（按列率排序）；第一种是按照佩利排列来定义（按自然排序）；第三种是按照哈达玛排列来定义。按哈达玛排序的沃尔什函数可由2n（n=0，1，2，…）阶哈达玛矩阵（H