第11章图像正交变换分析.ppt

第11章 图像正交变换 问题的提出： 视觉所感受到的是在空间域和时间域的信号。但是，往往许多问题在频域中讨论时，有其非常方便分析的一面。 图像变换的目的： 使图像处理问题简化 有利于图像特征提取 有助于从概念上增强对图像信息的理解　　　 变换问题的引入 什么是图像变换 将图像看成是线性叠加系统 图像在空域上相关性很强 图像变换是将图像从空域变换到其它域如频域的数学变换 常用的变换： 傅立叶变换、沃尔什变换、哈达玛变换、离散余弦变换、离散K-L变换、小波变换 先对行做变换： Fourier 变换示意图 Fourier变换的频率特性 Fourier变换的低通滤波 Fourier变换的高通滤波 Fourier变换的压缩原理 Fourier变换的压缩原理 快速傅里叶变换 问题的提出： 离散傅里叶变换已成为数字信号处理的重要工具。然而，它的计算量较大，运算时间长，在某种程度上却限制了它的使用范围。 11.2 离散余弦变换（DCT） Fourier变换的一个最大的问题是：它的参数都是复数，在数据的描述上相当于实数的两倍。为此，我们希望有一种能够达到相同功能但数据量又不大的变换。 在此期望下，产生了DCT变换。 一维离散余弦变换 二维离散余弦变换 离散余弦变换的计算 DCT变换的应用： 余弦变换实际上是傅立叶变换的实数部分。 余弦变换主要用于图像的压缩，如目前的国际压缩标准的JPEG格式中就用到了DCT变换。具体的做法与DFT 相似。给高频系数大间隔量化，低频部分小间隔量化。 程序生成一个矩形函数，区域内像素值为1，区域外为0。然后对矩形做二维傅立叶变换，由于图像的傅立叶变换矩阵元素一般是复数，不能直接显示，需要调用abs函数对变换后的结果求模，图中下面两幅图分别是傅立叶变换的对数谱和中心谱。 一维离散余弦反变换 二维离散反余弦变换 如果令N=4，由一维解析式定义可得如下展开式： 写成矩阵形式： [F(u)]=[A][f(x)] 同理，可得到反变换展开形式： 写成矩阵形式： [f(x)]=[A]T[F(u)] 二维离散余弦变换为： [F(u,v)]=[A][f(x,y)][A]T [f(x,y)]=[A]T[F(u,v)] [A]T 与傅立叶变换一样，离散余弦变换可以由定义出发进行计算，但这样的计算量太大，在实际应用中很不方便，寻找快速算法 首先，从定义出发，作如下推导 取实部的意思 如果把时域数据向量作系列延拓，即 (10)卷积定理： f(x,y)*h(x,y) = F(u,v)H(u,v) f(x,y)h(x,y) = F(u,v)*H(u,v) 二维傅立叶变换(幅值及相位)意义 左边一列: 上方为原始图像,下方为本图的相关说明说明; 中间一列: 上图幅值谱,下图为根据幅值谱的傅立叶逆变换(忽略相位信息,设相位为0); 右边一列: 上图相位谱,下图为根据相位谱的傅立叶逆变换(忽略幅值信息,设幅值为某一常数); 图像的说明 返回 返回 返回 另一幅图像效果 压缩率为：1.7：1 压缩率为：2.24：1 压缩率为：3.3：1 返回 压缩率为：8.1：1 压缩率为：10.77：1 压缩率为：16.1：1 二维离散傅立叶变换具有可分离性，即它可由两次一维 离散傅立叶变换计算得到，因此，仅研究一维离散傅立叶 变换的快速算法即可。改写公式： 式中，W=e-j2π／N ，称为旋转因子。 W= e-j2π／N =cos(2π／N )-j sin(2π／N ) (以N为周期) 式中很多Wux系数相同，不必进行多次重复计算。 FFT的推导过程： 设N为2的正整数次幂， 即 令M=N/2,离散傅立叶变换可改写成如下形式： 偶离散点 奇离散点 定义 于是 将一个N点的离散傅立叶变换分解成两个N／2短序列的离散傅立叶变换，即分解为偶数和奇数序列的离散傅立叶变换Fe(u)和Fo(u) 。 设N=23 7.2.2 快速离散傅立叶变换 7.2.2 快速离散傅立叶变换 蝶形运算单元 7.2.2 快速离散傅立叶变换 F e (0) F e (1) F e (2) F e (3) F o (0) F o (1) F o (2) F o (3) 0 8 W 1 8 W 2 8 W 3 8 W 0 8 W － 1 8 W － 2 8 W － 3 8 W － F (0) F (1) F (2) F (3) F (4) F (5) F (6) F (7) 7.2.2 快速离散傅立叶变换 Fe(u)和Fo(u)都是4点的DFT，对