《小波变换压缩算法解析ppt课件.ppt》由会员分享,可在线阅读,更多相关《小波变换压缩算法解析ppt课件.ppt(85页珍藏版)》请在三一办公上搜索。
1、第6章 小波变换压缩算法,2,主要内容,小波变换用于图像压缩的理由傅里叶变换窗口傅里叶变换小波变换的原理小波变换实例小波变换与数据压缩,3,小波变换用于图像压缩的理由,基于DCT (Discrete Cosine Transform) 的压缩标准JPEGMPEG-1,MPEG-2, H.264DCT 压缩的优点简单、 便于硬件实现,4,小波变换用于图像压缩的理由,DCT 压缩的缺点图像是分块处理,沿块的边界方向相关性被破坏,出现 “blocking artifacts”,5,傅里叶变换,信号表示多种方式信号的描述: 例如一个函数表达式,这就是信号的时域表示, 傅里叶变换1822年,傅里叶提出频
2、率的概念: 通过傅里叶正变换将信号在频域分解,获得信号的频谱,再通过反变换重建原始信号。频率仍然是傅里叶变换所定义 。,6,傅里叶变换,傅里叶变换的特点具有频域准确定位,可分析信号能量在各个频域成分中的分布情况,最常用的、最广泛的信号分析工具,并且相关的理论研究已发展为一个重要的数学分支调和分析。,7,傅里叶变换,傅里叶变换的不足缺乏时间-频率的定位功能不适于非平稳信号无法根据信号的特点自动调节时域和频域的分辨率,傅里叶变换的不足成为了推动寻找新变换的动力,8,窗口傅里叶变换,窗口傅里叶变换( short time Fourier transform )1946年Gabor提出了短时傅里叶变换
3、的概念 ,从而开始了非平稳信号的时频联合分析,9,窗口傅里叶变换,窗口傅里叶变换( short time Fourier transform ),10,窗口傅里叶变换,窗口傅里叶变换( short time Fourier transform )Gabor变换: 时窗函数Gauss函数时时窗函数的Fourier变换仍然是Gauss函数,保证了窗口傅立叶变换在频域内也有局域化的功能。,11,窗口傅里叶变换,窗口傅里叶变换( short time Fourier transform )时窗(Time Window),12,窗口傅里叶变换,窗口傅里叶变换( short time Fourier tr
4、ansform )频窗(Frequency Window)时窗函数g(t)的傅立叶变换 ,,13,窗口傅里叶变换,窗口傅里叶变换( short time Fourier transform )以上定义知,g(t)和G()分别起着时窗和频窗的作用,在时间频率坐标系中,时窗和频窗共同作用的结果就构成了时-频窗,这样就从几何上直观地描述了时频局部化。,14,窗口傅里叶变换,窗口傅里叶变换( short time Fourier transform )尽管窗式傅立叶变换能解决变换函数的局域化问题,但是,其窗口的大小和形状是固定的,即窗口面积不变,窗口没有自适应性。对于高频的信息,时间间隔要相对的小,更
5、好地确定峰值和断点,或者说需要用较窄的时域窗来反映信息的高频成分。对于低频谱的信息,时间间隔要相对的宽才能给出完整的信号信息, 或者说必须用较宽的时域窗来反映信息的低频成分。,15,小波变换原理,小波变换的( wavelet transform )发展20世纪 80年代后期发展起来的小波变换理论它是继傅里叶(Joseph Fourier)分析后信号处理与分析的强大工具无论是对古老的自然学科还是对新兴的高新技术应用学科都产生了强烈冲击。 小波理论是应用数学的一个新领域。要深入理解小波理论需要用到比较多的数学知识。从工程应用角度出发,直观的方法来介绍小波变换及其应用,为读者深入研究小波理论和应用提
6、供一些背景材料,16,小波变换原理,小波变换的( wavelet transform )发展哈尔(Alfred Haar)对在函数空间中寻找一个与傅里叶类似的基非常感兴趣。1909年他发现了小波,1910年被命名为Haar wavelets最早发现和使用了小波的名称,17,小波变换原理,小波变换的( wavelet transform )发展20世纪70年代,当时在法国石油公司工作的年轻的地球物理学家Jean Morlet提出了小波变换CWT (continuous wavelet transform)的概念。法国科学家Y.Meyer创造性地构造出具有一定衰减性的光滑函数,用缩放(dilati
7、ons)与平移(translations)均为 2的j次幂的倍数构造了平方可积的实空间L2(R)的规范正交基,使小波得到真正的发展.S.Mallat于1988年在构造正交小波基时提出了多分辨率分析(multiresolution analysis)的概念, 从空间上形象地说明了小波的多分辨率的特性,提出了正交小波的构造方法和快速算法,叫做Mallat算法。 Mallat算法地位相当于快速傅里叶变换在傅里叶分析中的地位。,18,小波变换原理,小波变换的( wavelet transform )发展1988年 Inrid Daubechies 最先揭示了小波变换和滤波器组(filter banks
8、)之间的内在关系20世纪90年代中期,Sweldens提出了小波变换提升方案- 第二代小波变换, 用于JPEG2000小波在信号(如声音信号,图像信号等)处理中得到极其广泛的应用。,19,小波变换原理,小波变换的( wavelet transform )发展小波变换具有在不同尺度下保持时频分析窗口面积不变性质自动调节对信号分析的时宽和带宽被誉为信号分析的显微镜,20,小波变换原理,连续小波变换( continuous wavelet transform )小波(Wavelet (A small wave, a ripple) 就是小的波形,所谓小,就是它具有衰减性,是存在于一个较小区域的波。,
9、21,小波变换原理,连续小波变换变换( continuous wavelet transform )小波基函数,22,小波变换原理,连续小波变换( continuous wavelet transform )小波正变换小波反变换,标注:a = scale variable 缩放因子b= time shift 时间平移在CWT中,缩放和平移是连续变化的,23,小波变换原理,连续小波变换( continuous wavelet transform )函数的伸缩,24,小波变换原理,连续小波变换( continuous wavelet transform )小波函数的伸缩,25,小波变换原理,连续小
10、波变换( continuous wavelet transform ),时窗中心 :小波 的时窗中心是其母函数 的时窗中心乘 倍再平移 个单位小波的 时窗宽度是其母函数 的时窗宽度的 倍。,26,小波变换原理,连续小波变换( continuous wavelet transform )小波的 频窗中心是其母函数 的频窗中心的 倍小波的 频窗宽度是其母函数 的频窗宽度的 倍,27,小波变换原理,连续小波变换( continuous wavelet transform )用较小 对信号做高频分析时,实际是用高频小波对信号进行细致观察用较大 对信号做低频分析时,实际是用低频小波对信号进行概貌观察,2
11、8,小波变换原理,连续小波变换( continuous wavelet transform )部分小波波形,29,小波变换原理,子带编码SBC (subband coding):把信号的频率分成几个子带,然后对每个子带 分别进行编码,并根据每个子带的重要性分配 不同的位数来表示数据20世纪70年代,子带编码开始用于语音编码20世纪80年代中期开始在图像编码中使用,30,小波变换原理,离散小波变换,图中的符号 表示频带降低1/2,HH表示频率最高的子带,LL表示频率最低的子带。这个过程可以重复,直到符合应用要求为止。这样的滤波器组称为分解滤波器树(decomposition filter tre
12、es),31,小波变换原理,离散小波变换只有离散,小波变换才能应用离散的方式有很多离散小波变换的多分辨率分析Mallat创立了多分辨率分析理论在多分辨率分析基础上,Mallat提出了基于滤波器组实现信号的小波正变换和反变换算法。执行离散小波变换的有效方法,32,小波变换原理,Mallat算法低通滤波器和高通滤波器构成双通道滤波原始的输入信号:S两个互补的滤波器A表示信号的近似值(approximations) D表示信号的细节值(detail),33,小波变换原理,Mallat算法低通滤波器和高通滤波器构成小波分解树对低频分量连续分解,34,小波变换原理,Mallat算法小波包分解树对低频分量
13、和高频分量均连续分解,35,小波变换原理,Mallat算法下采样过程 原始信号的数据样本为1000个,通过滤波之后每一个通道的数据均为1000个,总共为2000个。,36,小波变换原理,Mallat算法下采样过程 原始信号的数据样本为1000个,通过滤波之后每一个通道的数据均为1000个,总共为2000个。,37,小波变换原理,Mallat算法下采样过程 原始信号的数据样本为1000个,通过滤波之后每一个通道的数据均为1000个,总共为2000个。,38,小波变换原理,Mallat算法下采样过程 原始信号的数据样本为1000个,通过滤波之后每一个通道的数据均为1000个,总共为2000个。图中
14、的符号 表示下采样。,39,小波变换实例,一维哈尔小波变换哈尔函数定义,40,小波变换实例,一维哈尔小波变换哈尔函数定义基函数一组线性无关的函数,以用来构造任意给定的信号,41,小波变换实例,一维哈尔小波变换哈尔基函数最简单的基函数,42,小波变换实例,一维哈尔小波变换哈尔基函数,43,小波变换实例,一维哈尔小波变换尺度函数 :尺度函数尺度函数张成的空间VjVj的基的个数为2j,44,小波变换实例,一维哈尔小波变换小波函数 :与尺度函数对应哈尔小波函数与哈尔函数相对应,45,小波变换实例,一维哈尔小波变换小波函数 :与尺度函数对应哈尔小波函数与哈尔函数相对应哈尔小波基函数,46,小波变换实例,
15、一维哈尔小波变换小波基函数构成的空间:W j,47,小波变换实例,一维哈尔小波变换小波基函数构成的空间:W j,48,小波变换实例,一维哈尔小波变换小波基函数构成的空间:W j,49,小波变换实例,一维哈尔小波变换生成矢量空间W 2 的哈尔小波基函数,50,小波变换实例,一维哈尔小波变换生成矢量空间W 2 的哈尔小波基函数,51,小波变换实例,一维哈尔小波变换生成矢量空间W 2 的哈尔小波基函数,52,小波变换实例,一维哈尔小波变换实例图像9 7 3 5像素个数: 2j =22=4V2 中的哈尔基表示,53,小波变换实例,一维哈尔小波变换实例V2 中的哈尔基表示的一般形式其中的系数,54,小波
16、变换实例,一维哈尔小波变换实例用V 0, W 0和W1中的函数表示图像生成空间V 0的哈尔基函数为 生成空间W 0的哈尔小波基函数为 生成矢量空间W1的哈尔小波基函数为 和I(x)可表示成,55,小波变换实例,一维哈尔小波变换实例,56,小波变换实例,一维哈尔小波变换实例,57,小波变换实例,一维哈尔小波变换实例,58,小波变换实例,一维哈尔小波变换实例,59,小波变换实例,一维哈尔小波变换实例,60,小波变换实例,一维哈尔小波变换实例 生成其中,4个系数 , , 和 就是原始图像通过哈尔小波变换所得到的系数,用来表示整幅图像的平均值和不同分辨率下的细节系数。4个函数 , , 和 就是构成空间
17、V2的基。,61,小波变换实例,一维哈尔小波变换哈尔小波变换的快速算法 计算哈尔小波变换系数步骤1:求均值(averaging)。计算相邻像素对的平均值,得到一幅分辨率比较低的新图像,它的像素数目变成了2个,即新的图像的分辨率是原来的1/2,相应的像素值为:8 4步骤2:求差值(differencing)用2个像素表示这幅图像时,图像的信息已经部分丢失。为了能够从由2个像素组成的图像重构出由4个像素组成的原始图像,就需要存储一些图像的细节系数(detail coefficient),以便在重构时找回丢失的信息。原始图像可用下面的两个平均值和两个细节系数表示,8 4 1 -1步骤3:重复步骤1和
18、2把由第一步分解得到的图像进一步分解成分辨率更低的图像和细节系数。在这个例子中,分解到最后,就用一个像素的平均值6和三个细节系数2,1和1表示整幅图像:6 2 1 -1,62,小波变换实例,一维哈尔小波变换该算法可以推广到其他小波变换,63,小波变换实例,二维哈尔小波变换图像的二维变换,转置后继续对列实施变换相当于对行实施变换,64,小波变换实例,二维哈尔小波变换例如,65,小波变换实例,二维哈尔小波变换例如,66,小波变换实例,二维哈尔小波变换例如,67,小波变换实例,二维哈尔小波变换针对图像的小波变换的两种方法标准分解(standard decomposition)非标准分解(nonsta
19、ndard decomposition),68,小波变换实例,二维哈尔小波变换标准分解(standard decomposition)对图像的每一行进行小波变换,然后对这个经过行变换的图像的每一列进行小波变换,69,小波变换实例,二维哈尔小波变换标准分解(standard decomposition),70,小波变换实例,二维哈尔小波变换标准分解(standard decomposition)对图像每一行的像素值进行一维小波变换,再进行列变换,行变换与列变换交替进行,71,小波变换实例,二维哈尔小波变换非标准分解(nonstandard decomposition)交替地对图像的行和列进行小波
20、变换。,72,小波变换实例,二维哈尔小波变换非标准分解(nonstandard decomposition),73,小波变换实例,二维哈尔小波变换非标准分解(nonstandard decomposition),74,小波变换实例,二维哈尔小波变换非标准分解(nonstandard decomposition),75,小波变换实例,二维哈尔小波变换非标准分解(nonstandard decomposition),76,小波变换与数据压缩,原始图像,压缩图像,77,小波变换与数据压缩,压缩算法1(去除最小的系数),78,小波变换与数据压缩,压缩算法2(全局),79,小波变换与数据压缩,80,小波
21、变换与数据压缩,小波变换可以用于图像、声音(Sound)、视频压缩 小波的压缩过程通常分三部分:小波变换部分,量化部分,熵编码部分,81,小波变换与数据压缩,82,小波变换与数据压缩,为了进一步说明高频中的信息的所含内容,我们将一个经过压缩后的图像,逐步还原出来,每多一次还原过程,也就是经过一次小波变换的逆变换,图像逐渐清晰起来,见下图,83,小波变换与数据压缩,Time 2 Decomp Bowl,Time 1 Decomp basic content,Time 4 Decomp Background,Time 3Decomp Figures,84,小波变换与数据压缩,85,小波变换与数据压缩,量化,