一看就懂的小波变换

设有信号f(t):其傅里叶变换为F(jΩ):即：第1页/共38页第一页，共39页。=++Ψ（t）1/2Ψ（2t-t0）2/3Ψ（4t-t1）第2页/共38页第二页，共39页。像Ψ(t)这样，有限长且均值为0的函数称为小波函数。常用的小波函数如下图：第3页/共38页第三页，共39页。小波函数必须满足以下两个条件的函数：小波必须是振荡的；小波的振幅只能在一个很短的一段区间上非零，即是局部化的。如：图1小波例1图2小波例2第4页/共38页第四页，共39页。不是小波的例子图4图3第5页/共38页第五页，共39页。平均与细节设一维信号{x1，x2}

平均细节则一维信号可以表示成{a，d},且原信号可以恢复如下：

当x1与x2非常接近时，一维信号{x1，x2}可近似的用{a}表示，可实现信号压缩。

a可以看成信号的整体信息

d可看成原信号用a表示时丢失的细节信息第6页/共38页第六页，共39页。平均与细节对多元素信号{x1，x2，x3，x4}

信号可以表示为：{a1,0,a1,1,d1,0,d1,1}

丢失细节信号压缩为：{a1,0,a1,1}信号可进一步表示为：{a0,0,d0,0}丢失细节信号压缩为：{a0,0}第7页/共38页第七页，共39页。平均与细节{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页第八页，共39页。金字塔算法一维信号{3,1,-2,4}的小波变换为{1.5,0.5,1,-3}{1.5}：最低分辨率低频信息{0.5}：最低分辨率细节信息{2,1}：次高分辨率低频信息{1,-3}：次高分辨率细节信息{3,1,-2,4}：最高分辨率信息第9页/共38页第九页，共39页。尺度函数与小波函数信号序列{x1，x2，x3，x4}看成单位区间上的一个函数平移伸缩第10页/共38页第十页，共39页。引入记号：定义：可得：第11页/共38页第十一页，共39页。其它其它函数可以由一个尺度函数的伸缩与平移的线性组合表示第12页/共38页第十二页，共39页。同理，对小波变换其它伸缩和平移第13页/共38页第十三页，共39页。序列的多分辨率表示：第14页/共38页第十四页，共39页。4×4图像的二维Harr小波变换行小波变换列小波变换左上角二维小波变换第15页/共38页第十五页，共39页。1.1一维小波变换（一维多尺度分析）设有L2(R)空间的子空间序列：Vj的正交基函数是由一个称为尺度函数的函数(x)经伸缩平移得到的设Wj是Vj

相对于Vj+1的正交补空间，

Wj的正交基函数是由一个称为小波函数的函数(x)经伸缩平移得到的第16页/共38页第十六页，共39页。构成Vj+1的正交基。满足下列关系式(二尺度方程)：第17页/共38页第十七页，共39页。信号的多尺度分解：第18页/共38页第十八页，共39页。第19页/共38页第十九页，共39页。

求得小波系数的算式就是小波正变换。

该式也可以理解为f(x)和Ψa,b(x)内积，小波系数表示二者的相似程度，或f(x)中含有Ψa,b(x)成分的多少。第20页/共38页第二十页，共39页。

小波系数有a和b两个自变量，分别代表不同的尺度（时间）和频率，所以小波分析属于时频分析。第21页/共38页第二十一页，共39页。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页第二十二页，共39页。离散小波变换离散小波变换就是做向量的内积。例：对（64，2，3，61，60，6，7，57）做Haar小波变换第23页/共38页第二十三页，共39页。Haar小波变换第二种做法：第24页/共38页第二十四页，共39页。Haar小波反变换：第25页/共38页第二十五页，共39页。Haar小波反变换第二种做法：第26页/共38页第二十六页，共39页。1.2二维小波变换（二维多尺度分析）二维小波变换是由一维小波变换扩展而来的，二维尺度函数和二维小波函数可由一维尺度函数和小波函数张量积得到，即：图像的二维小波变换包括沿行向(水平方向)和列向(垂直方向)滤波和2-下采样，如图所示：第27页/共38页第二十七页，共39页。图5图像滤波采样第28页/共38页第二十八页，共39页。说明：如图所示，首先对原图像I(x,y)沿行向(水平方向)进行滤波和2->1下采样，得到系数矩阵IL(x,y)和IH(x,y)，然后再对IL(x,y)和IH(x,y)分别沿列向(垂直方向)滤波和2->1下采样，最后得到一层小波分解的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)的对角线方向细节子图第29页/共38页第二十九页，共39页。二维金字塔分解算法令I(x,y)表示大小为MN的原始图像，l(i)表示相对于分析小波的低通滤波器系数，i=0,1,2,…,Nl-1，Nl表示滤波器L的支撑长度；h(i)表示相对于分析小波的高通滤波器系数，i=0,1,2,…,Nh-1，Nh表示滤波器H的支撑长度，则第30页/共38页第三十页，共39页。第31页/共38页第三十一页，共39页。对逼近子图重复此过程，直到确定的分解水平，下图是二层小波分解的示意图。图6图像多尺度分解，(a)一层分解，(b)二层分解第32页/共38页第三十二页，共39页。

图像的小波特征提取首先对输入图像做J层二维小波分解；因为小波变换具有很好的时频局部化特性，所以可以将图像的不同底层特征变换为不同的小波系数；输入图像经过经一层小波分解后，被分成4个子图：

LL1—逼近子图，它代表输入图像水平和垂直两个方向的低频成分；

HL1—细节子图，它代表输入图像水平方向的高频成分和垂直方向的低频成分；第33页/共38页第三十三页，共39页。LH1—细节子图，它代表输入图像水平方向的低频成分和垂直方向的高频成分；

HH1—细节子图，它代表输入图像水平和垂直方向高频成分。在逼近子图LL1上重复二维小波分解过程，进行二层小波分解，如此继续分解，得到子图序列{LLJ，[HLk，LHk，HHk](k=1,2,…,J)}。小波基与分解层次的选取是非常重要的，目前还没有一个统一的标准。第34页/共38页第三十四页，共39页。I(x,y)[128128]I1(x,y)[6464]I1H(x,y)[6464]I1V(x,y)[6464]I1D(x,y)[6464]I2(x,y)[3232]I2H(x,y)[3232]I2V(x,y)[3232]I2D(x,y)[3232]I3(x,y)[1616]I3H(x,y)[1616]I3V(x,y)[1616]I3D(x,y)[1616]I4(x,y)[88]I4H(x,y)[88]I4V(x,y)[88]I4D(x,y)[88]图7图像I(x,y)的多尺度分解第35页/共38页第三十五页，共39页。小波基的选取一般考虑下列因素：线性相位：如果小波具有线性相位或至少具有广义线性相位，则可以避免小波分解和重构时的图像失真，尤其是图像在边缘处的失真；紧支性和衰减性：紧支性和衰减性是小波的重要性质，紧支宽度越窄或衰减越快，小波的局部化特性越好。计算复杂度越低，便于快速实现；

正交性：用正交小波基对图像做多尺度分解，可得一正交的镜像滤波器。低通子带数据和高通子带数据分别落在相互正交的L2(R2)的子空间中，使个子带数据相关性减少；其他

第36页/共38页第三十六页，共39页