《第七章小波变换和多分辨率处理ppt课件.ppt》由会员分享,可在线阅读,更多相关《第七章小波变换和多分辨率处理ppt课件.ppt(114页珍藏版)》请在三一办公上搜索。
1、第七章 小波变换和多分辨率处理,张 萍,电子科技大学 光电信息学院E-mail: ,参考资料,教材: Rafael C. Gonzalez, etc,Digital Image Processing (Third Edition),电子工业出版社, 2010参考书籍:(美)多布著,李建平译,小波十讲,国防工业出版社 ,2011孙延奎著,小波变换与图像、图形处理技术,清华大学出版社,2012朱希安,曹林编著,小波分析及其在数字图像处理中的应用,电子工业出版社,2012,“小波”(wavelet)就是一种“尺度”很小的波动,并具有时间和频率特性。,时间A,时间B,什么是小波?,小波函数必须满足以下
2、两个条件:小波必须是振荡的;小波的振幅只能在一个很短的一段区间上非零,即是局部化的。如:,小波变换具有良好的局部时频聚焦特性,而被称为“数学显微镜”。小波分析是纯数学、应用数学和工程技术的完美结合。从数学来说是大半个世纪“调和分析”的结晶(包括傅里叶分析、函数空间等)。小波变换是20世纪最辉煌科学成就之一。在信号处理、图像处理、模式识别、语音识别、量子物理、地震勘探、流体力学、电磁场、CT成象、机器视觉、故障诊断、分形、数值计算等已有重大突破。,小波分析发展简史,Inrid Daubechies于1988年最先揭示了小波变换和滤波器组(filter banks)之间的内在关系,使离散小波分析变
3、成为现实。Ronald Coifman和Victor Wickerhauser等著名科学家在把小波理论引入到工程应用方面做出了极其重要贡献。在信号处理领域中,自从Inrid Daubechies完善了小波变换的数学理论和Stephane Mallat构造了小波分解和重构的快速算法后,小波变换在各个工程领域中得到了广泛的应用,典型的如语音信号处理、医学信号处理、图像信息处理等。,小波理论与工程应用,傅里叶变换与小波变换,傅里叶变换的基础函数是正弦函数。小波变换基于一些小型波,称为小波,具有变化的频率和 有限的持续时间。,傅里叶变换与小波变换,傅里叶变换反映的是图像的整体特征,其频域分析具有很好的
4、局部性,但空间(时间)域上没有局部化功能。与傅里叶变换相比,小波变换是空间(时间)和频率的局部变换,它通过伸缩平移运算对信号逐步进行多尺度细化,最终达到高频处时间细分,低频处频率细分,能自动适应时频信号分析的要求,从而可聚焦到信号的任意细节。,小波变换是基于具有变化的频率和有限持续时间的小型波进行的。它是多分辨率理论的分析基础。多分辨率理论将多种学科的技术有效地统一在一起,其优势很明显某种分辨率下所无法发现的特性在另一种分辨率下将很容易被发现。本章将从多分辨率的角度解释小波变换。,主要内容,背景多分辨率展开一维小波变换快速小波变换二维小波变换小波包,主要内容,背景图象金字塔子带编码哈尔变换多分
5、辨率展开一维小波变换快速小波变换二维小波变换小波包,1.背景,物体的尺寸很小或者对比度不高的时候,通常采用较高的分辨率观察。物体尺寸很大或者对比度很强,只需要较低的分辨率。物体尺寸有大有小,强弱对比度同时存在,则适合用不同的分辨率对其进行研究。,从数学观点看,图像是一个亮度的二维矩阵,边界和强烈变化的区域局部直方图统计特性不同。无法对整个图象定义一个简单的统计模型。,一幅自然图像及其直方图的局部变化,1.背景,(1) 图像金字塔,以多分辨率来解释图像的一种简单有效的结构。一幅图像的金字塔是一系列以金字塔形状排列的分辨率逐步降低的图像集合。金字塔的底部是带处理图像的高分辨率表示,而顶部是低分辨率
6、的近似。当向金字塔的上层移动时,尺寸和分辨率就降低。,基础级J的大小为NN (J=log2N)顶点级0的大小为11 第j级的大小为2j2j (0j J)共有J+1级,但是通常我们截短到P1级,其中1 PJ,J-1级近似输出用来建立近似值金字塔;作为金字塔基级的原始图像和它的P级减少的分辨率近似都能直接获取并调整;J级的预测残差输出用于建立预测残差金字塔;近似值和预测残差金字塔都通过迭代计算获得。,金字塔方框图,(1) 图像金字塔,(1) 图像金字塔迭代算法,初始化,原始图象大小2J2J,jJj-1级,以2为步长进行子抽样,计算输入图像减少的分辨率近似值j-1级近似值,生成子抽样金字塔。对j-1
7、级近似值进行步长为2的内插,并进行过滤,生成与输入图像等分辨率的预测图像。计算输入图像和预测图像之间的差异,产生预测残差金字塔。重复2、3、4步骤。,图象的高斯近似值金字塔,分辨率分别为:512512,256256,128128,6464。金字塔的分辨率越低,伴随的细节越少;低分辨率图像用于分析大的结构或图像的整体内容,高分辨率图像用于分析单个物体的特性。,相应拉普拉斯预测残差金字塔,分辨率分别为:512512,256256,128128,6464。从低级开始通过内插和滤波获得高级高斯金字塔的预测残差图象。,(1) 图像金字塔,两种图像金字塔和它的统计特性。(a)高斯金字塔(近似),(b)拉普
8、拉斯金字塔(预测残差),(a),(b),子带编码也是多分辨率相关的重要图像技术在子带编码中,一幅图像被分解为一系列限带分量的几何,称为子带。子带可以重组在一起无失真地重建原始图象。每个子带通过对输入进行带通滤波而得到。子带带宽小于原始图像带宽,子带可以进行无信息损失的抽样原始图象的重建可以通过内插、滤波、和叠加单个子带来完成,(2) 子带编码,系统输入是一个一维的带限时间离散信号x(n)分析滤波器h0(n)和h1(n)是半波数字滤波器,理想传输函数H0,H1如下图所示。H0低通滤波,输出x(n)的近似值H1高通滤波,输出x(n)的高频或细节部分综合滤波器g0(n)和g1(n) 为重构的结果,(
9、2) 子带编码,(a)一维子带编码和解码的两频带滤波器组,(b)频谱分离特性,(a),(b),序列x(n)的Z变换时域以2为因子的抽样对应到Z域同样,以2为因子的内插对应的变换为X(n)先抽样再内插得到,(2) 子带编码,系统输出滤波h0(n)的输出整理第二项含有z,代表了抽样-内插过程带来的混叠,(2) 子带编码,对于输入的无失真重建,假定下列条件:矩阵表达,消除混叠消除幅度失真,&,(2) 子带编码,分析调制矩阵,证明分析滤波器和综合滤波器双正交,24,由性质,以及奇次方相互抵消,(2) 子带编码,将H0和G0表示成G1和H1的函数,(2) 子带编码,满足该条件的滤波器组称为具有双正交性,
10、分析滤波器和综合滤波器满足上述条件,所以具有双正交性,(2) 子带编码,(正交镜像滤波器),(共轭正交滤波器),完美重建滤波器族,一维滤波器用于图像处理的二维可分离滤波器可分离滤波器首先应用于某一维(如水平方向),在应用于另一维(如垂直方向) a(m,n), dV(m,n), dH(m,n)和dD(m,n)分别表示近似值、垂直细节、水平细节和图象的对角线细节子带,(2) 子带编码,4个8抽头Daubechies正交滤波器的冲激响应。,(2) 子带编码,低通滤波器h0(n)的系数为:-0.010 597 40,0.032 883 01,0.030 841 38,-0.187 034 81,-0.
11、027 983 76,0.630 880 76,0.714 846 57,0.230 377 81。其余正交滤波器参数可以通过公式计算获得。,a dHdV dD,(2) 子带编码,花瓶的4频段子带编码,花瓶图像从512512到256256子带的近似子带a、水平子带dH,垂直子带dV和对角线子带dD dH和dV有混叠是由于对可分辨窗口进行抽样造成的,可以通过综合滤波器重建时消除。,它的基函数是最普遍也最简单的正交小波。哈尔变换本身对称、可分离,矩阵表示: T=HFHTF是NN图象矩阵,H是NN变换矩阵,T是NN变换的结果哈尔基函数,(3) 哈尔变换,NN哈尔变换矩阵第i行包含元素hi(z),其中
12、z = 0/N, 1/N, , (N-1)/N。例如,N=4时, k,p,q的值如右: 则,44变换矩阵H4 22变换矩阵H2,(3) 哈尔变换,离散小波变换的哈尔函数,6464,128128,256256,图示为哈尔基函数对图像的多分辨率分解,离散小波变换包含了与原始图像相同的像素数其局部统计数据相对稳定,并且容易给出模型。大多数数据接近0,可以进行大量的压缩;原始图象的粗和细分辨近似可以从中提取子图象进行重建。,主要内容,背景多分辨率展开序列展开尺度函数小波函数一维小波变换快速小波变换二维小波变换小波包,2.多分辨率展开,图像金字塔、子带编码和哈尔变换,在数学理论多分辨率分析中扮演了重要角
13、色。在多分辨率分析( MRA )中,尺度函数被用于建立某一函数或图像的一系列近似值,相邻两近似值之间的近似度相差2倍。被称为小波的附加函数用于对相邻近似值之间的差异进行编码。,信号或函数可以分解为一系列展开函数的线性组合 其中,k是有限或无限和的整数下标, k是具有实数值的展开系数, k(x)是具有实数值的展开函数 如果展开方式唯一,则任何指定的f(x)只有一个k序列与之相对应k(x)称为基函数展开序列k(x)称为可表示这一类函数的基。,(1) 序列展开,可展开的函数组成了一个函数空间,被称为展开集合的闭合跨度系数k可以通过内积得到,(1) 序列展开,由于展开集合的正交性,计算有3种形式情况1
14、:如果展开函数构成了V的一个正交基基与它的对偶相等, 则系数由基函数和原函数的内积来计算,(1) 序列展开,由于展开集合的正交性,计算有3种形式情况2: 若展开函数本身不正交,而是V的正交基基函数及其对偶函数为双正交,即系数由下式计算,(1) 序列展开,情况3:如果展开集合对V来说不是函数基,但是符合线性展开,那么它是一个跨度集合,对于任一f (x)V有一个以上k集合。展开函数及其对偶称为超完备或冗余。它们组成了一个框架,其中:对于某些A0,B,及所有f(x)V。若 A=B,(1) 序列展开,考虑整数平移和实数二值尺度、平方可积函数 (x)组成的可展开函数集合k决定了j,k(x)在x轴的位置j
15、决定了j,k(x)的宽度2j/2控制其高度或幅度j,k(x)的形状随着j发生变化, (x)被称为尺度函数通过选择适当的(x),j,k(x)可以决定跨度L2(R),所有可量度的平方可积函数的集合.,(2) 尺度函数,定义代表任何j,k上的跨度子空间增大j,用于表示子空间函数的j,k(x)范围变窄增加j将增加Vj的大小,将允许具有变化较小的变量和较细节函数包含在子空间中举例说明哈尔尺度函数,(2) 尺度函数,单位高度单位宽度的尺度函数,当j=1时与j=0相反,展开函数更窄更密集对于左下角的函数将0,0(x)分解作为V1展开函数的和,可以得到从数学角度看,V0是V1的一个子空间,记做,举例:哈尔尺度
16、函数,简单尺度函数遵循多分辨率的四个基本条件MRA要求1:尺度函数对其积分变换是正交的在哈尔函数的情况下,因为无论什么时候只要尺度函数的值是1,其积分变换就是0,所以二者的乘积是0。哈尔函数是紧支撑的,即,除被称为支撑区的有限区间外,函数值都为0。事实上,其支撑区是1,半开区间0,1)外的支撑区的值是0。必须注意,当尺度函数的支撑区大于1时。积分变换正交的要求将很难满足。,(2) 尺度函数,简单尺度函数遵循多分辨率的四个基本条件MRA要求2:由低尺度函数跨越的子空间在低尺度处嵌套在由高尺度跨越的子空间内。,(2) 尺度函数,子空间:,简单尺度函数遵循多分辨率的四个基本条件MRA要求3:唯一包含
17、在所有Vj中的函数是 f(x)=0如果考虑可能的最粗糙的展开函数(即j=-),惟一可表达的函数就是没有信息的函数,即,,(2) 尺度函数,简单尺度函数遵循多分辨率的四个基本条件MRA要求4:任何函数都可以以任意精度表示虽然在任意粗糙的分辨率下展开一个特定f(x)是几乎不可能的,但所有可度量的、平方可积函数都可以用极限j表示,即,在这些条件下,子空间Vj的展开函数可以被表述为子空间Vj+1的展开函数的加权和。,(2) 尺度函数,被称为尺度函数系数; 为尺度矢量任意子空间的展开函数都可以从它们自身的双倍分辨率拷贝中得到,即从相邻较高分辨率的空间中得到。对引用子空间V0的选择是任意的。,多分辨率分析
18、的基础,称为改善等式、MBR等式或扩张等式,(2) 尺度函数,哈尔尺度函数系数,附加的简化产生了,Haar尺度和小波函数是不连续和紧支撑的,在支撑的有限区域外是0.,举例:哈尔尺度函数系数,给定满足MRA要求的尺度函数,能够定义小波函数(x) (与它的积分变换及其二进制尺度) ,跨越了相邻两个尺度子空间Vj和Vj+1的差异。,(3) 小波函数,用尺度函数可得如果f(x)=Wj尺度函数与小波函数的关系表示空间并集Vj+1中Vj的正交补集是Wj,Vj中所有成员对于Wj中的所有成员都正交,(3) 小波函数,所有可量度的、平方可积函数空间表示为:任何小波函数可以表示为平移的双倍分辨率尺度函数的加权和
19、被称为小波函数系数; 为小波向量利用小波跨越的正交补集空间、积分小波变换是正交的条件,可得,(3) 小波函数,哈尔小波函数系数哈尔尺度向量定义相应的小波向量哈尔小波函数,举例:哈尔小波函数系数,W1比W0窄,可以标志更细微的细节;函数展开这里V0尺度函数的近似W0小波函数,举例:哈尔小波函数系数,尺度函数(Scaling function, , phi) 父小波函数; 近似空间(低频);小波函数(Wavelet function, psi) 母小波函数; 细节空间(高频),主要内容,背景多分辨率展开一维小波变换小波序列展开离散小波变换连续小波变换快速小波变换二维小波变换小波包,j0 是任意开始
20、尺度 近似值或尺度系数 细节或小波系数第一个和式是 f(x) 在尺度 j0 上近似第二个和式是在较高尺度 jj0 更细分辨率的小波函数,(1) 小波序列展开,定义小波序列展开,展开系数计算:,3.一维小波变换,举例:哈尔小波序列展开,举例:哈尔小波序列展开,如果待展开函数是一个数字序列,得到的系数就称为离散小波变换(DWT)对于 j j0通常令j00且M是2的幂(M=2J)注意:序列展开中的积分变成了求和对于双正交函数要用对偶函数 代替,近似系数细节系数,(2) 离散小波变换,考虑4点离散函数:f(0)=1, f(1)=4, f(2)=-3和f(3)=0M=4, J=2, 对于x=0, 1,
21、2, 3, j=0, 1求和,举例:一维离散小波函数,考虑4点离散函数:f(0)=1, f(1)=4, f(2)=-3和f(3)=0重构原始函数,举例:一维离散小波函数,Morlet 小波,举例:一维离散小波函数,Mexihat小波,连续平方可积函数f(x)的连续小波变换(CWT)s和分别为尺度和变换参数当满足条件 反连续小波变换,其中 是 的傅里叶变换,(2) 连续小波变换,DWT和CWT的相似性连续变换参数取代了积分变换参数连续尺度参数s与二进制尺度参数2j相反。连续尺度参数s出现在分母上,小波尺度和通常意义上的频率定义相反0 1 时,扩大或展开CWT开始展开j0=-,消除了尺度函数间的明
22、显关联,函数只包括小波项和DWT相似,CWT可以被看成是一组变换系数 ,它给出f(x)与基函数集 的相似性。 在连续情况下,两个集合都是无穷的,(2) 连续小波变换,墨西哥草帽小波,左图函数的傅里叶变换,解释了尺度化的小波和傅里叶频段之间的联系,频谱中的两个显著频段(峰值)对应函数的两个类高斯扰动,函数根据墨西哥草帽小波完成的CWT的一部分 (1 s 10且 100)它同时给出了时域和频域的信息。如:当s=1,变换在=10时达到最大值,变换绝对值 用黑白之间的灰度级显示,主要内容,背景多分辨率展开一维小波变换快速小波变换二维小波变换小波包,FWT是实现DWT的高效计算,利用相邻尺度DWT系数之
23、间的关系,采用类似的两段子带编码方案,也称为Mallat人字型算法。从多分辨率等式开始,尺度向量h可以被看成是用来将 展开成尺度为j+1的尺度函数和的“权”,4.快速小波变换FWT,4.快速小波变换FWT,注意:DWT在尺度j的细节系数是尺度在j+1是近似值系数的函数,类似的可以得到,尺度为j的近似值 和细节系数 可以通过尺度为j+1的近似系数 和时域反转的尺度与小波向量 的卷积,而后对结果进行亚取样来计算,4.快速小波变换FWT,一个FWT分析滤波器族有h0(n) = h (-n) 且 h0(n) = h(-n)其中,卷积在n=2k时进行计算,在非负偶数时刻计算卷积与以2为步长进行过滤和抽样
24、的效果相同,4.快速小波变换FWT,FWT分析滤波器可以迭代产生多阶结构。如果以高于奈奎斯特频率的采样率进行采样,该样值是该分辨率的尺度系数的良好近似,可以作为起始的高分辨率尺度系数的输入,即,无须小波系数。,二阶两尺度FWT分析滤波器以及其频率分离特性,4.快速小波变换FWT,例:计算一维小波变换,离散函数f(n)=1,4,-3,0,计算基于哈尔尺度和小波函数的变换。使用小波向量,例:计算一维小波变换,离散函数f(n)=1,4,-3,0,计算基于哈尔尺度和小波函数的变换。使用小波向量,哈尔尺度和小波向量计算序列1,4,-3,0的二尺度快速小波变换,通过近似值W (j,k)和W(j,k)细节系
25、数重建f(x)的高效反变换,称为快速小波反变换(FWT-1)使用正变换中所用的尺度和小波向量以及第j级近似值和细节系数来生成第j1级的近似值系数完备重建要求对于i=0,1, gi(n) = hi(-n)即:分析滤波器和综合滤波器在时域中是相互反转的对于双正交分析/综合滤波器,不是彼此时域反转的,4.快速小波变换FWT,FWT1的综合滤波器族,FWT1迭代:二阶或两尺度FWT1的综合滤波器族,Wup代表步长为2 内插,例:计算一维小波反变换,首先对0级近似值和细节系数进行内插,产生1,0和4,0与滤波器g0(n)和 g1(n)卷积结果相加产生W(1,n),得到一级近似值的重建迭代运算,重建f(x
26、),用哈尔尺度和小波向量计算序列 的两尺度快速小波反变换,快速小波变换与FFT的比较,运算复杂性对于FWT,长度为M=2J的序列的FWT的运算次数是 O(M) 阶,即:浮点乘法和加法(使用滤波器族)的次数与序列的长度存在这线性关系FFT需要 O(MlogM) 阶变换的基函数傅里叶的基函数(正弦函数)保证了FFT的存在FWT的存在取决于使用的小波函数的尺度函数是否存在,以及尺度函数和相应的小波函数的正交/双正交性表达函数时,时间和频率通常被作为不同的域来处理,它们之间存在这不可分割的关系例如,要得到时域有价值的信息,就要忍受频域模糊,反之亦燃-海森伯测不准原理块不重叠是正交基函数的特点,FWT和
27、FFT的比较,标准时域基给出时间发生的时刻,没有频域信息正弦基给出时间发生的频率但是没有时间分辨率FWT时间和频率分辨率是变化的低频:块短而宽,即有较好的频率分辨率,对应较差的时间分辨率高频:块窄而高,即有较高的时间分辨率,频率分辨率下降,取样数据 FFT基函数 FWT基函数,时间一频率块,主要内容,背景多分辨率展开一维小波变换快速小波变换二维小波变换小波包,二维乘积可分离的尺度函数二维可分离方向敏感小波定义尺度和平移基函数,沿列方向变化,沿行方向变化,沿对角线方向变化,5.二维小波变换,MN的函数f(x, y)的离散小波变换同一维DFT一样, 定义了尺度在j0 的f(x, y) 的近似, 系
28、数对于jj0 附加了水平,垂直对角线方向的细节离散反小波变换,5.二维小波变换,二维DWT可以用数字滤波器和抽样来实现,先取f(x,y)的行的一维FWT再用结果列的一维FWT,5.二维小波变换,图7.22 二维快速小波变换。(a)分析滤波器族,二维DWT可以用数字滤波器和抽样来实现,5.二维小波变换,图7.22 二维快速小波变换。 (b)分解结果,两尺度分解结果,综合滤波器,重建过程和一维相似每一次迭代4尺度j的近似和细节图象用两个一维滤波器内插和卷积,5.二维小波变换,图7.22 二维快速小波变换。 (c)综合滤波器族,四阶对称小波ab分解滤波器cd重建滤波器e一维小波函数f一维尺度函数,a
29、 bc de f,低通重建滤波g0(n)=h(n)的系数对于 0 n 7:0.0322,-0.0992,0.2979,0.8037,0.4976,-0.0296,-0.0758,5.二维小波变换,例:计算二维小波变换,四阶对称小波 三个二维小波之一,H(x,y),5.二维小波变换,3尺度FWT,a,b,c,d,a图:由计算机产生的128X128基于滤波器分解的一个序列b-d图:对称小波滤波器分解结果,前面分析的小波族称为“对称小波”(symlet)它们不完全对称,但具有最小不对称性和最高消失矩数DWT尺度和小波函数表现为低通和高通滤波器特性,大多数基于傅里叶滤波器的技术和小波部分是等价的,小波
30、 的第k阶矩是 。0阶矩影响小波函数和尺度函数的平滑性以及多项式表示它们的能力。一个N阶对称小波有N个消失矩,5.二维小波变换,小波在图像处理中的用途,如在傅里叶域那样,基本方法是:计算一幅图像的二维小波变换修改变换计算反变换,5.二维小波变换,例:基于小波的边缘提取,a bc d,a图消除了最低尺度近似分量,再进行反变换得到b图,效果是强调和突出了图象的边缘c图将水平细节也置为0,反变换得到重建图象d,可以孤立出垂直边缘。,图7.25 改进的边缘检测DWT。(a)(c)选择的 删去系数的两尺度分解,(b)(d)相应的重建,例:基于小波的去除噪声,基于小波的图象去噪的过程选择分解用的一个小波(
31、例如:哈尔对称小波)和级别数或尺度P。计算噪声图象的FWT门限化细节系数从尺度J-1到J-P选择应用一个门限处理细节系数硬门限实现,元素绝对值低于门限值则置0软门限实现,元素绝对值低于门限值则置0,并且标定非0的系数接近0,去除了硬门限固有的门限处的不连续性基于原始的近似系数,在J-P级执行小波重建,并对J-1到J-P级改进细节系数,对噪声去除改进DWTa图:人体头部带噪声MRI图象b图显示了门限化细节系数后的重建图像4阶对称小波2尺度(P2)全局门限94.909 3c图:最高分辨率细节置0的重建图像e图:两个分解级别的细节被置0后的DWT重建d图:在c图重建时移去的信息包含了原图象大多数噪声
32、和某些边缘信息f图:在e图重建时移去的信息,a b c de f,例:基于小波的去除噪声,主要内容,背景多分辨率展开一维小波变换快速小波变换二维小波变换小波包,快速小波变换将一个函数分解为一系列与对数相关的频段低频被组成窄频段高频被组成宽频段想要较大的控制时频平面的一部分,FWT必须有更灵活的分解小波包产生过程的代价是FWT计算复杂度增加,从O(M)到O(MlogM),6. 小波包,考虑两阶的滤波器族,分解过程用二叉树表示根节点被赋予最高的尺度近似系数,它是函数自身的取样叶子继承变换的近似细节和系数细节的输出,两尺度FWT分析族的一个系数数和分析数,6. 小波包,三尺度FWT分析族、分析数和相
33、应的频谱,三尺度FWT分析滤波框图,分解空间树,谱分离特性,6. 小波包,分析树提供了多尺度小波变换的紧凑有效的方法比对应的滤波器和基于子取样的方框图更容易画,并占有较少的空间相对容易定位有效分解三阶分析数提供了三种展开选择,6. 小波包,分析树还是表示小波包的有效机理,从3阶FWT分析树到3阶小波包树A表示近似滤波D表示细节滤波3阶小波包树几乎是3阶FWT的有效分解数目的3倍,6. 小波包,注意:平均分布的频带是完全小波包分解的特征。,分析树的滤波器组,谱分离特性,随着扩展的增加,基于包的变换改进了对被分解函数的频谱分割的控制,代价是复杂度增加,滤波器组将 分解成 输出,反复迭代生成P尺度变
34、换第一次迭代得到,为适应二维输入特性,就要有单个小波子空间对应系数,子空间分析树,6. 小波包,三尺度、全小波包分解树(部分),6. 小波包,例:二维小波包分解,指纹扫描图象,三尺度、全小波包分解为64片,上图中被分解子图象构成88阵列的子带为了达到压缩的目的,64片分解达到某种程度优化的可能性相对较低选择一种更加合理的分解方式附加代价函数测量二维函数f的熵或信息量利用代价函数将接近0值的数目最大化优化算法必须使用附加代价函数使分解树的叶子节点所付出的代价最小:熵最小的节点具有更多接近0的数值,导致更大的压缩,例:二维小波包分解,构造最小熵的有效算法对于分析树的每个节点,从根节点开始逐层构造树
35、,直到叶子节点计算节点的熵和此节点的4个子节点的熵对于二维小波包分解,父节点是一个近似值或细节系数的二维阵列,子节点是经过滤波的近似值是水平、垂直、和对角线方向上的细节如果子节点的联合熵小于父节点的熵,就将这些子节点包括到分析树中,否则去掉这些子节点,只保留父节点用上述算法修剪小波包树,根据计算最优树的大致框架设计处理程序。,例:二维小波包分解,最佳小波包分解结果及最佳小波包分析树,105,图中没有被进一步分离的子图象比较平滑,由具有中间灰度值的象素构成。近似值子图象已被标定,以灰度128表示0值系数,这些子图象几乎不包含任何信息。这些子图象的分割不会增加图象总体熵。,Cohen-Daubec
36、hies-Feauveau双正交小波族滤波器族的尺度和小波函数是对称的并且具有相近的长度,双小波及尺度函数,6. 小波包,Cohen-Daubechies-Feauveau双正交小波族,107,分解滤波器,重构滤波器,0n 17时,低通滤波器系数h0(n): 0, 0.001 9, -0.001 9, -0.017, 0.011 9, 0.049 7, -0.077 3, -0.094 1, 0.420 8, 0.825 9, 0.420 8, -0.094 1, -0.0773, 0.049 7, 0.011 9, -0.017, -0.001 9, 0.001 0高通滤波器系数h1(n):
37、 0, 0, 0, 0.0144 4, -0.014 5, -0.078 7, 0.040 4, 0.417 8, -0.758 9, 0.417 8, 0.040 4, -0.078 7, -0.014 5, 0.014 4, 0, 0, 0, 0g0(n) = (-1)n+1h1(n) 和 g1(n)=(-1)nh0(n),小 结,小波变换是强有力的时频分析工具,是在克服傅立叶变换缺点的基础上发展而来的。已成功应用于很多领域,如信号处理、图像处理、模式识别等。 小波变换的一个重要性质是它在时域和频域均具有很好的局部化特征,它能够提供目标信号各个频率子段的频率信息。这种信息对于信号分类是非常
38、有用的。 小波变换一个信号为一个小波级数,这样一个信号可由小波系数来刻画。,小波变换的应用,信号稀疏表示与重构;信号、图像去噪;图像边缘检测、目标检测;数据融合、图像融合;特征提取、模式识别;红外图像背景抑制、目标识别等。,小 结,超小波(Beyond Wavelet)分析技术,Curvelet 曲波 Ridgelet 脊波 Bandelet 带波 Beamlet 束波,小线波 Wedgelet 楔波 Surfacelet 面小波 Directionlet 方向波 Contourlet 轮廓波 Shearlet 剪切波,多尺度几何分析.,本章主要内容,背景多分辨率展开一维小波变换快速小波变换二维小波变换小波包,本章要求及作业,本章要求: 1. 熟悉小波分析理论及技术;2.掌握如何应用小波分析技术应用在数字图像上的基本方法。,本章作业,1. 书后: 7.2, 7.9, 7.212. 课后编程:(1)任意选择一幅图像,实现该图像的Haar小波的两层子带分解,显示分解图像。然后进行重构,显示重构图像。,END,