小波分析基础学习资料学习教案.pptx

《小波分析基础学习资料学习教案.pptx》由会员分享，可在线阅读，更多相关《小波分析基础学习资料学习教案.pptx（76页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、会计学1小波分析小波分析(fnx)基础学习资料基础学习资料第一页，共76页。如图1所示的LENA图像(txin)f(x,y)，假设图像(txin)的大小是512x512，量化级是256，即xy第2页/共76页第二页，共76页。2、L2(R)空间的正交分解和变换 1对f(t)L2(R)，存在(cnzi)L2(R)的一组标准正交基gi(t)，tR，i=1,2,使得其中(qzhng)(1.2)(1.3)第3页/共76页第三页，共76页。对于给定信号f(t)，关键是选择合适的基gi(t)，使得f(t)在这组基下的表现呈现出我们需要的特性，但是如果某一个基不满足要求，可通过变换将函数转换到另一个基下表示

2、，才能得到(d do)我们需要的函数表示。常用的变换2有：(1)K-L变换(2)Walsh变换(3)傅立叶变换(4)小波变换 如图所示是信号f(t)的傅立叶变换示意图。信号f(t)经傅立叶变换由时域变换到频域，基底不同得到(d do)大变换也不同。在信号处理中，有两类非常重要的变换即傅立叶变换和小波变换。目前，可简单地将小波理解为满足以下两个条件的特殊信号：小波必须时振荡的；小波的振幅只能在一个很短的一段区间上非零，即是局部化的。第4页/共76页第四页，共76页。1、Daubechies小波一些(yxi)著名的小波3：第5页/共76页第五页，共76页。2、Coiflets小波3、Symlets

3、小波第6页/共76页第六页，共76页。4、Morlet小波5、MexicanHat小波6、Meyer小波SKIPSKIP第7页/共76页第七页，共76页。不是(bshi)小波的例第8页/共76页第八页，共76页。RETURNRETURN第9页/共76页第九页，共76页。3、傅立叶变换与时频分析4我们知道，任何复杂的周期信号f(t)可以用简单的调和振荡函数(hnsh)表示成如下形式：这就是(jish)著名的傅立叶级数，都是简单(jindn)的调和振荡函数，直观讲都是正弦波。是函数f(t)的傅立叶系数，可由以下公式计算：(1.4)第10页/共76页第十页，共76页。于是，周期函数f(t)就与下面(

4、ximian)的傅立叶序列产生了一一对应，即从数学上已经证明了，傅立叶级数(jsh)的前N项和是原函数f(t)在给定能量下的最佳逼近：(1.5)(1.6)(1.7)第11页/共76页第十一页，共76页。对于(duy)L2(R)上的非周期函数f(t)，有(1.8)(1.9)称为f(t)的傅立叶变换，反变换公式为(1.10)第12页/共76页第十二页，共76页。有了傅立叶变换，我们可以很容易地将时域信号 f(t)转换到频域上，于是信号的频率特性一目了然，并且与傅立叶级数一样，傅立叶变换将一段信号的主要低频能量都集中在频率信号的前面几项，这种能量集中性有利于进一步的处理。在过去200年里，傅立叶分析

5、在科学与工程领域发挥了巨大的作用，但傅立叶分析也有不足，主要表现在以下两点：q傅立叶分析不能刻画(khu)时域信号的局部特性；q傅立叶分析对非平稳信号的处理效果不好。q下面通过例子来说明这两点。第13页/共76页第十三页，共76页。例、歌声信号歌声是一种声音震荡的波函数，其傅立叶变换就是将这个波函数转化成某种乐谱(yup)。但遗憾地是，傅立叶变换无法反映信号在哪一时刻有高音，在哪一时刻有低音，因此结果是所有的音符都挤在了一起，如图所示。第14页/共76页第十四页，共76页。小波变换有效地克服(kf)了傅立叶变换的这一缺点，信号变换到小波域后，小波不仅能检测到高音与低音，而且还能将高音与低音发生

6、的位置与原始信号相对应，如图所示。第15页/共76页第十五页，共76页。因此我们需要这样一个数学工具：既能在时域很好地刻画信号的局部性，同时也能在频域反映信号的局部性，这种数学工具就是“小波”。从函数分解的角度，希望能找到另外一个基函数(t)来代替sint。(t)应满足以下三个特性：任何复杂的信号f(t)，都能由一个母函数(t)经过伸缩和平移(pny)产生的基底的线性组合表示；信号用新的基展开的系数要能反映出信号在时域上的局部化特性；新的基函数(t)及其伸缩平移(pny)要比三角基sint更好地匹配非平稳信号。历史上，Haar第一个找到了这样一个基函数，这就是非常著名但又及其简单的Haar小波

7、。(1.11)第16页/共76页第十六页，共76页。数学(shxu)上已经证明：小波级数(jsh)、信号的小波逼近构成L2(R)的一个正交基，通过规范化处理，(1.12)构成L2(R)的一个规范正交基。故任何一个能量有限信号 f(t)L2(R)可以分解为(1.13)(1.14)(1.15)第17页/共76页第十七页，共76页。二、小波变换二、小波变换(binhun)的定义及特的定义及特点点定义定义1 1函数函数(t)L2(R)称为基本小波，如果它满足称为基本小波，如果它满足(mnz)以下的以下的“允许允许”条件：条件：(2.1)如果是连续的，易得：(2.2)第18页/共76页第十八页，共76页

8、。(t)又称为母小波，因为其伸缩(shnsu)、平移可构成L2(R)的一个标准正交基：同傅立叶变换一样，连续(linx)小波变换可定义为函数与小波基的内积：将a,b离散(lsn)化，令可得离散小波变换：(2.3)(2.4)(2.5)第19页/共76页第十九页，共76页。总结：小波即小区域的波，是一种特殊(tsh)的长度有限、平均值为零的波形。它有两个特点：一是“小”，即在时域具有紧支集或近似紧支集；二是正负交替的“波动性”，也即支流分量为零。傅立叶分析是将信号分解成一系列不同频率的正弦波的叠加，同样小波分析是将信号分解为一系列小波函数的叠加，而这些小波函数都是由一个母小波函数经过平移和尺度伸缩

9、得来的。(2.6)(2.7)第20页/共76页第二十页，共76页。小波分析优于傅立叶分析的地方是，它在时域和频域同时具有良好的局部化性质。而且由于对高频成分采用(ciyng)逐渐精细的时域或频域取样步长，从而可以聚焦到对象的任何细节，所以被称为“数学显微镜”。小波分析广泛应用与信号处理、图像处理、语音识别等领域。第21页/共76页第二十一页，共76页。可以这样理解小波变换的含义：打个比喻，我们用镜头观察目标信号f(t)，(t)代表镜头所起的所用。b相当于使镜头相对于目标平行移动，a的所用相当于镜头向目标推进或远离。由此可见，小波变换有以下特点：多尺度/多分辨的特点，可以由粗及细地处理信号；可以

10、看成用基本频率特性为()的带通滤波器在不同尺度a下对信号做滤波。适当地选择小波，使(t)在时域上为有限支撑,()在频域上也比较(bjio)集中，就可以使WT在时、频域都具有表征信号局部特征的能力。第22页/共76页第二十二页，共76页。小波变换的思想来源于伸缩和平移(pny)方法。尺度伸缩对波形的尺度伸缩就是在时间轴上对信号进行压缩和伸展，如图所示。第23页/共76页第二十三页，共76页。第24页/共76页第二十四页，共76页。v时间平移v时间平移就是指小波函数在时间轴上的波形(bxn)平行移动，如图所示。第25页/共76页第二十五页，共76页。小波运算的基本步骤：(1)选择一个小波函数，并将

11、这个小波与要分析的信号起始点对齐；(2)计算(jsun)在这一时刻要分析的信号与小波函数的逼近程度，即计算(jsun)小波变换系数C，C越大，就意味着此刻信号与所选择的小波函数波形越相近，如图所示。第26页/共76页第二十六页，共76页。(3)将小波函数沿时间轴向右移动一个单位时间，然后重复步骤(1)、(2)求出此时的小波变换系数C，直到覆盖完整个信号(xnho)长度，如图所示；第27页/共76页第二十七页，共76页。(4)将所选择的小波函数尺度伸缩一个单位，然后(rnhu)重复步骤(1)、(2)、(3)，如图所示；(5)对所有的尺度伸缩重复(chngf)步骤(1)、(2)、(3)、(4)。第

12、28页/共76页第二十八页，共76页。v尺度与频率(pnl)的关系尺度与频率的关系如下：小尺度a压缩的小波快速变换的细节高频部分大尺度a拉伸(lshn)的小波缓慢变换的粗部低频部分第29页/共76页第二十九页，共76页。第30页/共76页第三十页，共76页。三、多分辨三、多分辨(fnbin)分析分析由母小波按如下方式(fngsh)的伸缩平移可构成L2(R)空间的标准正交基如何构造母小波呢？1989年，Mallat和Meyer提出了按多分辨分析的思想来构造母小波，其基本思想是：现构造一个具有特定性质的层层嵌套的闭子空间序列 VjjZ，这个(zhge)闭子空间序列充满了整个 L2(R)空间。在V0

13、子空间找一个函数g(t)，其平移g(t-k)kZ构成V0子空间的Riesz基。对函数g(t)进行正交化，得到函数称为正交尺度函数(t)。由(t)计算出小波函数(t)。1 1、多分辨分析、多分辨分析、多分辨分析、多分辨分析(MRA)(MRA)的概念的概念的概念的概念55(3.1)第31页/共76页第三十一页，共76页。Riesz基定义令H是Hilbert空间，H中的一个(y)序列gjjZ是Riesz基，如果它满足以下的条件：A和B分别称为(chnwi)Riesz基的上下界，Riesz基又称为(chnwi)稳定基。(3.2)(3.3)第32页/共76页第三十二页，共76页。定义定义1 空间L2(R

14、)中的多分辨分析是指L2(R)中的满足如下条件的一个子空间序列第33页/共76页第三十三页，共76页。多分辨空间的关系可用下图来形象(xngxing)地说明。第34页/共76页第三十四页，共76页。如果g(t-k)kZ是V0的Riesz基，可通过(tnggu)正交化得到V0空间的函数(t)V0，使得(t-k)kZ构成V0空间的规范正交基。由伸缩性和平移不变性可知，j,k(t)j,kZ构成Vj空间的一个规范正交基。于是(ysh)(3.4)(3.5)第35页/共76页第三十五页，共76页。注意：(t)并不是L2(R)空间的小波函数，而是与其紧密(jnm)相关的尺度函数，j,k(t)j,kZ称为尺度

15、基，多分辨空间序列VjjZ称为尺度空间，在MRA意义下，可由尺度基导出小波基。由MRA的单调性可以看出：Vj是Vj+1的严格子空间，设Wj是Vj关于Vj+1的正交补(子空间)，即(3.6)第36页/共76页第三十六页，共76页。对于一幅图像，量化级数决定了图像的分辨率，量化级数越高，图像就越清晰(qngx)，即图像的分辨率高。对于任意一幅图像，都可以用不同的量化空间来表示，细节比较丰富的部分用高分辨率来表示，细节比较单一的部分可用低分辨率来表示。我们可以将不同的量化级数构成的空间看成不同的多分辨空间Vj，显然这些量化空间是相互嵌套的，(3.7)第37页/共76页第三十七页，共76页。从图像处理

16、的角度，多分辨空间的分解可以理解为图像的分解，假设有一幅 256级量化的图像，不妨将它看成量化空间 Vj中的图像，则可理解为Vj空间中的图像有一部分保留在 Vj-1空间中，还有一部分放在 Wj-1空间，如图所示如图所示如图所示如图所示。与尺度函数的产生一样(yyng)，若存在(t)W0，使得(t-k)kZ构成空间W0的一个规范正交基，则构成L2(R)空间的一个规范正交基。称为小波基，(t)称为母小波。(3.8)SKIPSKIP第38页/共76页第三十八页，共76页。VjWj-1Vj-1RETURNRETURN第39页/共76页第三十九页，共76页。MRA非常抽象，但是它给出了构造小波的一般框架

17、。在实践中很难通过小波空间直接构造小波，但通过MRA可推导出一个非常重要(zhngyo)的关系：双尺度方程，通过求解该方程，使我们有可能求出尺度函数和小波函数。2、双尺度方程由前面的分析，我们知道：(3.9)(3.10)第40页/共76页第四十页，共76页。方程(fngchng)(3.9)和(3.10)称为双尺度方程(fngchng)。由(t)的正交性可得：对双尺度方程两边(lingbin)取傅立叶变换，可得频域上的的双尺度方程：(3.12)(3.11)(3.14)(3.13)第41页/共76页第四十一页，共76页。(3.16)(3.15)从信号处理的角度，h是与(t)对应的低通滤波器，g是与

18、(t)对应的高同滤波器，h，g既可以表示(biosh)为时域上的离散序列形式 hk，gkkZ，也可以表示(biosh)为频域上的2周期函数h()，g()。两者本质上是一样的。第42页/共76页第四十二页，共76页。若kN时，hk=0，这样的滤波器称为有限(yuxin)脉冲响应滤波器(FIR)，FIR滤波器具有好的局部化特性。此时，(t)只在有限(yuxin)区间0，N上取值，所以(t)是紧支的，其支集supp=0，N，(3.9)式变为：(3.17)此时(t)也是紧支的。所以只要(zhyo)滤波器的长度是有限的，我们称对应的小波(t)是紧支小波。第43页/共76页第四十三页，共76页。由(3.1

19、3)式得：(3.18)(3.19)第44页/共76页第四十四页，共76页。结论：结论：结论：结论：只要找到满足双尺度方程(3.9)的序列hkkZ，通过公式(3.15)就可以计算出2周期函数h()，再由公式(3.19)就可以计算出，经过傅立叶反变换，最终可得尺度函数(t)，有了尺度函数就可以计算出小波函数(t)。通过解双尺度方程(3.9)，我们希望得到满足MRA的尺度函数(t)，并最终构造出小波函数(t)，但有两个问题必须解决：问题1：双尺度方程(3.9)是否有解？解的唯一性如何？问题2：双尺度方程(3.9)的解是否满足MRA？关于(guny)问题1，I.Daubechies和Lagarias7

20、 在1991 年给出了证明。第45页/共76页第四十五页，共76页。解决问题2却是一件非常困难的事情。这里牵涉到尺度函数(t)与滤波器系数hkkZ之间的关系问题：如果有一个L2(R)空间的尺度函数(t)，一定能构造出双尺度方程(3.9)，从而找到一组满足(3.9)的滤波器hkkZ；反过来，如果有一组滤波器hkkZ满足某个双尺度方程，由此求解得到的函数却不一定是满足MRA的尺度函数，这样无法(wf)保证双尺度方程解的平移构成L2(R)Riesz基若(t)是正交的，则相应的滤波器h有什么性质呢？定理13若(t)是正交的，则相应的滤波器hk必须满足条件：(3.20)(3.21)但是，如果hk仅仅满足

21、(3.20)和(3.21)，并不能保证由双尺度方程构造出的函数(hnsh)(t)是正交尺度函数(hnsh)。(3.20)和(3.21)称为构造正交小波的必要条件。第46页/共76页第四十六页，共76页。仅有必要条件是不够的，即 hkkZ除了满足(mnz)条件(3.20)和(3.21)外，还应满足(mnz)其他条件。S.Mallat4，W.Lawton6等都在这方面作出了重大的贡献，并给出了一些有意义的结论。下面给出 W.Lawton的充分条件。定理x2设h()是FIR滤波器，若满足(mnz)若矩阵A的特征值1是非(shfi)退化的，则(t-k)kZ是标准正交的。第47页/共76页第四十七页，共

22、76页。算法：构造紧支小波基算法：构造紧支小波基算法：构造紧支小波基算法：构造紧支小波基步骤步骤步骤步骤1 1 寻找满足寻找满足寻找满足寻找满足(mnz)(mnz)双尺度方程双尺度方程双尺度方程双尺度方程(3.9)(3.9)和和和和(3.10)(3.10)的滤波器的滤波器的滤波器的滤波器hk,gkkhk,gkk0,1,N0,1,N步骤步骤步骤步骤2 2 利用公式利用公式利用公式利用公式(3.15)(3.15)计算计算计算计算2 2周期函数周期函数周期函数周期函数h(h()；步骤步骤步骤步骤3 3 验证验证验证验证h(h()是否满足是否满足是否满足是否满足(mnz)(mnz)条件条件条件条件通过

23、傅立叶反变换(binhun)求出(t)步骤5验证矩阵A的特征值1是否非退化；步骤6(t-k)kZ是正交的尺度函数，对应的紧支小波由公式(3.10)计算。步骤步骤步骤步骤(bzhu)4(bzhu)4(bzhu)4(bzhu)4 计算计算计算计算第48页/共76页第四十八页，共76页。3 3、小波与共轭镜像滤波器、小波与共轭镜像滤波器、小波与共轭镜像滤波器、小波与共轭镜像滤波器44 我们知道尺度函数和小波函数我们知道尺度函数和小波函数我们知道尺度函数和小波函数我们知道尺度函数和小波函数(t),(t),(t)t(t)t RR是在时域刻画是在时域刻画是在时域刻画是在时域刻画(khu)(khu)信号的性

24、质，对应的滤波器信号的性质，对应的滤波器信号的性质，对应的滤波器信号的性质，对应的滤波器h(h(),g(),g()RR从频域上刻画从频域上刻画从频域上刻画从频域上刻画(khu)(khu)信号的性质。实际上，信号的性质。实际上，信号的性质。实际上，信号的性质。实际上，(t),(t),(t)t(t)t RR大量的性质都可以由对应的大量的性质都可以由对应的大量的性质都可以由对应的大量的性质都可以由对应的h(h(),g(),g()RR从频域上反映出来，甚至离散小波变换都可以借助滤波器来实现，因此小波与滤波器具有紧密的关系。从频域上反映出来，甚至离散小波变换都可以借助滤波器来实现，因此小波与滤波器具有紧

25、密的关系。从频域上反映出来，甚至离散小波变换都可以借助滤波器来实现，因此小波与滤波器具有紧密的关系。从频域上反映出来，甚至离散小波变换都可以借助滤波器来实现，因此小波与滤波器具有紧密的关系。3.1 3.1 正交尺度函数产生共轭镜像滤波器正交尺度函数产生共轭镜像滤波器正交尺度函数产生共轭镜像滤波器正交尺度函数产生共轭镜像滤波器定义定义定义定义 若尺度函数若尺度函数若尺度函数若尺度函数(t)(t)是正交的，则它所对应的滤波器是正交的，则它所对应的滤波器是正交的，则它所对应的滤波器是正交的，则它所对应的滤波器h(h()称为共轭镜像滤波器。称为共轭镜像滤波器。称为共轭镜像滤波器。称为共轭镜像滤波器。h

26、(h()满足以下条件：满足以下条件：满足以下条件：满足以下条件：第49页/共76页第四十九页，共76页。滤波器hkkZ称为低通滤波器。所谓低通是指：当信号(xnho)f(t)被hkkZ作用后，其低频成分能被保留下来，而高频成分(=)却被滤掉了。对应的小波滤波器g()也是共轭镜像滤波器。也满足条件(3.22)另外，由于(t-k)kZ与(t-k)kZ分别(fnbi)是V0空间和W0空间的规范正交基，而 V0W0，则(3.23)公式(gngsh)(3.23)反映了低通滤波器h()和高通滤波器g()之间的关系。RETURNRETURN第50页/共76页第五十页，共76页。S.Mallat4同时给出了这

27、样的结论：若高通滤波器g()满足(mnz)公式(3.22)和(3.23)，则由公式产生的小波基(t-k)kZ构成W0空间的规范正交基。因此当尺度函数(t)已经确定时，只要能找到一个满足公式(3.22)和(3.23)的g()，就一定能找到对应(duyng)的小波(t)，但是这样的解并不是唯一的。例如可取(3.24)可以验证(ynzhng)g()满足(3.22)和(3.23)，对应的共轭镜像滤波器为：(3.25)第51页/共76页第五十一页，共76页。因此当找到低通共轭镜像滤波器hkkZ后，利用公式(gngsh)(3.25)马上可得高通共轭镜像滤波器gkkZ。总结：在一个MRA下的正交尺度函数和小

28、波函数(t),(t)tR，产生一组共轭镜像滤波器h,g，满足：(3.26)公式(3.26)还有几个(j)等价形式，下面以定理的形式给出。第52页/共76页第五十二页，共76页。定理定理定理定理 设设设设h,gh,g是由正交尺度函数和小波函数产生的共轭镜像滤波器，则以下几个条件等价是由正交尺度函数和小波函数产生的共轭镜像滤波器，则以下几个条件等价是由正交尺度函数和小波函数产生的共轭镜像滤波器，则以下几个条件等价是由正交尺度函数和小波函数产生的共轭镜像滤波器，则以下几个条件等价(dngji)(dngji)：在频域上在频域上在频域上在频域上(3.26)(3.26)式成立；式成立；式成立；式成立；在时

29、域上以下公式成立：在时域上以下公式成立：在时域上以下公式成立：在时域上以下公式成立：(3.27)q定义调制矩阵：(3.28)第53页/共76页第五十三页，共76页。则(3.29)3.2 3.2 利用共轭镜像滤波器实现快速正交小波变换利用共轭镜像滤波器实现快速正交小波变换利用共轭镜像滤波器实现快速正交小波变换利用共轭镜像滤波器实现快速正交小波变换44 L2(R)L2(R)空间空间空间空间(kngjin)(kngjin)的一个的一个的一个的一个MRAMRA产生了两个子空间产生了两个子空间产生了两个子空间产生了两个子空间(kngjin)(kngjin)：尺度空间：尺度空间：尺度空间：尺度空间(kng

30、jin)Vjj(kngjin)Vjj Z Z和小波空间和小波空间和小波空间和小波空间(kngjin)Wjj(kngjin)Wjj Z Z。j,kj,kj,kj,k Z Z和和和和 j,kj,kj,kj,k Z Z 分别是两个空间分别是两个空间分别是两个空间分别是两个空间(kngjin)(kngjin)的规范正交基，信号的规范正交基，信号的规范正交基，信号的规范正交基，信号f(t)f(t)L2(R)L2(R)在两个空间在两个空间在两个空间在两个空间(kngjin)(kngjin)上都可以做正交投影：上都可以做正交投影：上都可以做正交投影：上都可以做正交投影：(3.30)第54页/共76页第五十四

31、页，共76页。信号(xnho)在小波空间的展开为(3.31)但实践(shjin)中不可能进行无穷次逼近，不妨设f(t)VJ，则因为所以(suy)表示从尺度2-J到2-j进行了(J-j)次小波分解(jJ)第55页/共76页第五十五页，共76页。实际计算时，可以一次一次地进行小波分解，然后递推实现(J-j)次小波分解，不妨记一次小波分解的尺度(chd)系数和小波系数为(3.32)第56页/共76页第五十六页，共76页。而因为(ynwi)代入(3.32)式得故第57页/共76页第五十七页，共76页。从而(cngr)我们(wmen)得到如下的递推公式：(3.33)第58页/共76页第五十八页，共76页

32、。现在(xinzi)来求dj,k的递推公式，(3.34)而第59页/共76页第五十九页，共76页。因为(ynwi)代入(3.34)式得故第60页/共76页第六十页，共76页。从而(cngr)我们得到如下(rxi)的递推公式：(3.35)第61页/共76页第六十一页，共76页。通过公式(3.33)和(3.35)，可以很快计算出尺度(chd)系数和小波系数cj,k,dj,k，这就是著名的Mallat算法：因此，只要确定(qudng)VJ空间的初始序列cJ,kk Z，就可以算出任意空间 Vj(jJ)的所有尺度系数和小波系数。公式(3.33)和(3.35)称为离散小波变换的分解公式。第62页/共76页

33、第六十二页，共76页。又由于Vj+1=VjWj，VjWj，因此Vj上的标准正交基与Wj上的标准正交基是相互正交的。它们共同构成(guchng)Vj+1上的标准正交基，则 Vj+1上的函数j+1,nj,nZ可以由这两个基共同表示：有前面(qinmian)的计算可知：故第63页/共76页第六十三页，共76页。从而(cngr)这就是(jish)Mallat重构算法：第64页/共76页第六十四页，共76页。小波的应用小波的应用(yngyng)1,4,8,9小波的应用主要是信号的处理，其中最典型的应用是小波图象压缩。另外，小波在诸如信号去噪、特征提取等多方面均有成功的应用。下面以图象去噪为例说明小波应用

34、策略。小波的各种应用均可分为以下三步：1）对原始信号作小波变换(binhun)，将信号由空域变换(binhun)到频域；2）对小波系数做相应处理；3）对处理后的小波系数做小波逆变换(binhun)，还原原信号。第65页/共76页第六十五页，共76页。小小小小波波波波信信信信号号号号(x x n nh h o o)去去去去噪噪噪噪一一一一第66页/共76页第六十六页，共76页。小小小小波波波波信信信信号号号号(x x n nh h o o)去去去去噪噪噪噪二二二二第67页/共76页第六十七页，共76页。图图图图像像像像(t t x xi i n n)融融融融合合合合第68页/共76页第六十八页，

35、共76页。小波图像去噪8因为噪声信号(xnho)多包含在具有较高频率的细节中，所以小波去噪首先对图像信号(xnho)进行小波分解，可利用门限阈值对所分解的小波系数进行处理，然后对图像信号(xnho)进行小波重构，抑制图像信号(xnho)中的无用部分，恢复图像信号(xnho)中的有用部分。具体步骤为：图像信号(xnho)的小波分解：选择合适的小波及恰当的分解层次N，对目标图像进行N层的小波分解；对分解后的高频系数进行阈值量化：对于分解的每一层，选择恰当的阈值，对该层高频系数进行阈值量化处理；重构图像：根据小波分解后的第N层近似的低频系数和经过阈值量化处理后的细节高频系数，重构图像。第69页/共7

36、6页第六十九页，共76页。第70页/共76页第七十页，共76页。参考文献参考文献1唐远炎,王玲.小波分析与文本文字识别(shbi),科学出版社,20042李弼程,彭天强,彭波.智能图像处理技术,电子工业出版社,20043I.Daubechies,TenLecturesonWavelets.Philadelphia:SIAM,1992.4S.Mallat.Awavelettourofsignalprocessing.AcademicPress,USA,19985S.Mallat,“Atheoryformultiresolutionsignaldecomposition:Thewaveletrep

37、resentation,”IEEETrans.PatternAnal.MachineIntell.,vol.11,pp.674693,1989.6W.Lawton.Tightframesofcompactlysupportedwavelets,J.Math.Phys.,31:18981901,1990.第71页/共76页第七十一页，共76页。7I.Daubechies,J.C.Lagarias.Two-scaledifferenceoperationsI:existenceandglogalregularityofsolutions,SIAMJ.Math.Anal.,22:13881410,1

38、991.8D.Donoho,“De-noisingbysoft-thresholding,”IEEE Trans.Inform.Theory,vol.41,pp.613627,1995.9B.JAWERTH,etc.“anoverviewofwaveletbasedmultiresolutionanalyses,”SIAMREVIEW,vol.36,No.33,pp.377-412,September,1994.第72页/共76页第七十二页，共76页。精品精品(jn pn)课件课件！第73页/共76页第七十三页，共76页。精品精品(jn pn)课件课件！第74页/共76页第七十四页，共76页。第75页/共76页第七十五页，共76页。感谢您的观看感谢您的观看(gunkn)！第76页/共76页第七十六页，共76页。