06_第四章_正交变换和小波变换.ppt

正交变换的一般表达 当f(x)或f(x,y)为偶函数时，变换的计算公式只有余弦项。 一个任意函数采样从0,1,2,…,N-1，若向负方向折叠形成2N采样的偶函数，就可以进行2N的偶函数傅立叶变换。 余弦变换是简化傅立叶变换的一种方法 1822年Fourier变换,在频域的定位最准确，无任何时域定位能力。函数，时域定位完全准确，频域无任何定位能力 1946年Gabor变换，STFT，窗函数的大小和形状与时间和频率无关而保持固定不变。不构成正交基。 1982年Burt提出金字塔式图像压缩编码，子带编码(subband coding),多采样率滤波器组(multirate sampling filter bank). 1910年Harr提出规范正交基。 1981年Stormberg对Harr系进行改进，证明了小波函数的存在。 1984年,Morlet提出了连续小波 1985年,Meyer,Grossmann,Daubecies提出离散的小波基 1986年,Meyer证明了不可能存在时域频域同时具有正则性的正交小波基，证明了小波的自正交性。 1987年,Mallat统一了多分辨率分析和小波变换，给出了快速算法。 1988年，Daubecies在NSF的小波专题研讨会进行了讲座。 小波的应用 J.Morlet，地震信号分析。 S.Mallat，二进小波用于图像的边缘检测、图像压缩和重构 Farge，连续小波用于涡流研究 Wickerhauser，小波包用于图像压缩。 Frisch噪声的未知瞬态信号。 Dutilleux语音信号处理 H.Kim时频分析 Beykin正交小波用于算子和微分算子的简化 本章要求 1. 了解DFT的运算，掌握其性质 2. 掌握DHT和DWT的计算及其特点 本章作业 书95-96页 ： 3, 8, 10, 12, 13 2. 离散小波变换 离散小波变换（连续小波的冗余性） 能不能完整表征信号。 任意信号是否都能表示为基本单元的加权和。系数如何确定。 2. 离散小波变换 线性变换 [Tx]=x(t),yj(t) 成为一个标架: 唯一性 正变换和反演的连续性 标架：一般标架，紧标架，正交标架 通过标架进行重建 2. 离散小波变换 小波基在离散栅格上扩展为小波标架。 基本小波满足框架条件为小波框架 满足小波框架条件就是满足连续小波容许条件。 小波标架： 正交小波 半正交小波 双正交小波 二维小波变换（二维多尺度分析） 二维小波变换是由一维小波变换扩展而来的，二维尺度函数和二维小波函数可由一维尺度函数和小波函数张量积得到，即： 图像的二维小波变换包括沿行向(水平方向)和列向(垂直方向)滤波和2-下采样，如图所示： 图像滤波采样 说明：如图所示，首先对原图像I(x,y)沿行向(水平方向)进行滤波和2-下采样，得到系数矩阵IL(x,y)和IH(x,y)，然后再对IL(x,y)和IH(x,y)分别沿列向(垂直方向)滤波和2-下采样，最后得到一层小波分解的4个子图： ILL (x,y)—I(x,y)的（粗）逼近子图 IHL(x,y) — I(x,y)的水平方向细节子图 ILH (x,y) — I(x,y)的垂直方向细节子图 IHH (x,y) — I(x,y)的对角线方向细节子图 二维金字塔分解算法 令I(x,y)表示大小为M?N的原始图像，l(i)表示相对于分析小波的低通滤波器系数，i=0,1,2,…,Nl-1， Nl表示滤波器L的支撑长度； h(i)表示相对于分析小波的高通滤波器系数，i=0,1,2,…,Nh-1， Nh表示滤波器H的支撑长度，则 对逼近子图重复此过程，直到确定的分解水平，下图是二层小波分解的示意图。 图像多尺度分解，(a)一层分解，(b)二层分解 3. 多分辨率分析和多采样滤波器 理想滤波器频带分解 带宽减半－采样减半 各带通空间恒Q性 滤波器一致性 函数空间分解 小波函数和尺度函数为基底的函数空间分解 3. 多分辨率分析和多采样滤波器 与连续小波的联系 各个分辨率进行分析的高频信号就是各个尺度的小波变换 信号上采样后通过同一滤波器等于原来信号经过一个插0值的滤波器。 滤波器组和小波基可以互求，利用各自领域的设计技术 H0 H1 2 a0(n) a1(n) 2 2 2 G0 G1 3. 多分辨率分析和多采样滤波器 多采样滤波器的重建条件 H0(-z)G0(z) + H1 (-z)G1(z) = 0 H0(z)G0(z) + H1(z)G1(z) = cz-k H0 H1 2 a0(n) a1(n) 2 2 2 G0 G1 3. 多分辨率分析和多采样滤波器 滤波器组四个滤波器，两个必须条件，所以不唯一，通过其他条件约束产生特定滤波器。 正交镜像对称滤波器组 H0,H1