一看就懂的小波变换.pptx

《一看就懂的小波变换.pptx》由会员分享，可在线阅读，更多相关《一看就懂的小波变换.pptx（38页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、设有信号f(t):其傅里叶变换为F(j):即：第1页/共38页=+（t）1/2（2t-t0）2/3（4t-t1）第2页/共38页像(t)这样，有限长且均值为0的函数称为小波函数。常用的小波函数如下图：第3页/共38页小波函数必须满足以下两个条件的函数：(1)小波必须是振荡的；(2)小波的振幅只能在一个很短的一段区间上非零，即是局部化的。如：图1 小波例1图2 小波例2第4页/共38页不是小波的例子图4图3第5页/共38页平均与细节设一维信号x1，x2 平均 细节 则一维信号可以表示成a，d,且原信号可以恢复如下：当x1与x2非常接近时，一维信号x1，x2可近似的用a表示，可实现信号压缩。a可以

2、看成信号的整体信息 d可看成原信号用a表示时丢失的细节信息 第6页/共38页平均与细节对多元素信号x1，x2，x3，x4 信号可以表示为：a1,0,a1,1,d1,0,d1,1 丢失细节信号压缩为：a1,0,a1,1信号可进一步表示为：a0,0,d0,0丢失细节信号压缩为：a0,0第7页/共38页平均与细节x1，x2，x3，x4最高分辨率信息a1,0,a1,1次高分辨率低频信息d1,0,d1,1次高分辨率细节信息a0,0最低分辨率低频信息d0,0最低分辨率细节信息x1，x2，x3，x4的小波变换a0,0,d0,0,d1,0,d1,1由整体平均和两个不同分辨率的细节信息构成第8页/共38页金字塔

3、算法一维信号3,1,-2,4的小波变换为1.5,0.5,1,-31.5：最低分辨率低频信息0.5：最低分辨率细节信息2,1：次高分辨率低频信息1,-3：次高分辨率细节信息3,1,-2,4：最高分辨率信息第9页/共38页尺度函数与小波函数信号序列x1，x2，x3，x4看成单位区间上的一个函数平移伸缩第10页/共38页引入记号：定义：可得：第11页/共38页其它其它函数可以由一个尺度函数的伸缩与平移的线性组合表示第12页/共38页同理，对小波变换其它伸缩和平移第13页/共38页序列的多分辨率表示：第14页/共38页44图像的二维Harr小波变换行小波变换列小波变换左上角二维小波变换第15页/共38

4、页1.1 一维小波变换（一维多尺度分析）设有L2(R)空间的子空间序列：Vj 的正交基函数是由一个称为尺度函数的函数(x)经伸缩平移得到的设Wj 是Vj 相对于Vj+1的正交补空间，Wj 的正交基函数是由一个称为小波函数的函数(x)经伸缩平移得到的第16页/共38页构成Vj+1的正交基。满足下列关系式(二尺度方程)：第17页/共38页信号的多尺度分解：第18页/共38页第19页/共38页 求得小波系数的算式就是小波正变换。该式也可以理解为f(x)和a,b(x)内积，小波系数表示二者的相似程度，或f(x)中含有a,b(x)成分的多少。第20页/共38页 小波系数有a和b两个自变量，分别代表不同的

5、尺度（时间）和频率，所以小波分析属于时频分析。第21页/共38页Haar小波（1/8,1/8,1/8,1/8,1/8,1/8,1/8,1/8）（1/8,1/8,1/8,1/8,-1/8,-1/8,-1/8,-1/8）（1/4,1/4,-1/4,-1/4,0,0,0,0）（0,0,0,0,1/4,1/4,-1/4,-1/4）（1/2,-1/2,0,0,0,0,0,0）（0,0,1/2,-1/2,0,0,0,0）（0,0,0,0,1/2,-1/2,0,0）（0,0,0,0,0,0,1/2,-1/2）连续Haar小波对应的离散Haar小波第22页/共38页离散小波变换离散小波变换就是做向量的内积。例

6、：对（64，2，3，61，60，6，7，57）做Haar小波变换第23页/共38页Haar小波变换第二种做法：第24页/共38页Haar小波反变换：第25页/共38页Haar小波反变换第二种做法：第26页/共38页1.2 二维小波变换（二维多尺度分析）二维小波变换是由一维小波变换扩展而来的，二维尺度函数和二维小波函数可由一维尺度函数和小波函数张量积得到，即：图像的二维小波变换包括沿行向(水平方向)和列向(垂直方向)滤波和2-下采样，如图所示：第27页/共38页图5 图像滤波采样第28页/共38页说明：如图所示，首先对原图像I(x,y)沿行向(水平方向)进行滤波和2-1下采样，得到系数矩阵IL(

7、x,y)和IH(x,y)，然后再对IL(x,y)和IH(x,y)分别沿列向(垂直方向)滤波和2-1下采样，最后得到一层小波分解的4个子图：q ILL(x,y)I(x,y)的（粗）逼近子图q IHL(x,y)I(x,y)的水平方向细节子图q ILH(x,y)I(x,y)的垂直方向细节子图q IHH(x,y)I(x,y)的对角线方向细节子图第29页/共38页二维金字塔分解算法令I(x,y)表示大小为MN的原始图像，l(i)表示相对于分析小波的低通滤波器系数，i=0,1,2,Nl-1，Nl表示滤波器L的支撑长度；h(i)表示相对于分析小波的高通滤波器系数，i=0,1,2,Nh-1，Nh表示滤波器H的

8、支撑长度，则第30页/共38页第31页/共38页对逼近子图重复此过程，直到确定的分解水平，下图是二层小波分解的示意图。图6 图像多尺度分解，(a)一层分解，(b)二层分解第32页/共38页q 图像的小波特征提取首先对输入图像做J层二维小波分解；q 因为小波变换具有很好的时频局部化特性，所以可以将图像的不同底层特征变换为不同的小波系数；q 输入图像经过经一层小波分解后，被分成4个子图：LL1逼近子图，它代表输入图像水平和垂直两个方向的低频成分；HL1细节子图，它代表输入图像水平方向的高频成分和垂直方向的低频成分；第33页/共38页 LH1细节子图，它代表输入图像水平方向的低频成分和垂直方向的高频

9、成分；HH1细节子图，它代表输入图像水平和垂直方向高频成分。q 在逼近子图LL1上重复二维小波分解过程，进行二层小波分解，如此继续分解，得到子图序列LLJ，HLk，LHk，HHk(k=1,2,J)。q 小波基与分解层次的选取是非常重要的，目前还没有一个统一的标准。第34页/共38页I(x,y)128128I1(x,y)6464I1H(x,y)6464I1V(x,y)6464I1D(x,y)6464I2(x,y)3232I2H(x,y)3232I2V(x,y)3232I2D(x,y)3232I3(x,y)1616I3H(x,y)1616I3V(x,y)1616I3D(x,y)1616I4(x,y

10、)88I4H(x,y)88I4V(x,y)88I4D(x,y)88图7 图像I(x,y)的多尺度分解第35页/共38页小波基的选取一般考虑下列因素：q 线性相位：如果小波具有线性相位或至少具有广义线性相位，则可以避免小波分解和重构时的图像失真，尤其是图像在边缘处的失真；q 紧支性和衰减性：紧支性和衰减性是小波的重要性质，紧支宽度越窄或衰减越快，小波的局部化特性越好。计算复杂度越低，便于快速实现；q 正交性：用正交小波基对图像做多尺度分解，可得一正交的镜像滤波器。低通子带数据和高通子带数据分别落在相互正交的L2(R2)的子空间中，使个子带数据相关性减少；q 其他 第36页/共38页第37页/共38页Thanks！第38页/共38页