小波分析及其应用(精).ppt

1,小波分析及其应用WaveletAnalysisandItsApplications,2,小波分析及其应用,1、小波变换简介2、小波分析在一维信号处理中的应用3、小波分析在图象分析中的应用图象特征抽取图象压缩数据隐藏和图象水印,3,小波变换简介,1.1小波变换的理论基础,信号分析是为了获得时间和频率之间的相互关系。傅立叶变换提供了有关频率域的信息，但有关时间的局部化信息却基本丢失。与傅立叶变换不同，小波变换是通过缩放母小波（Motherwavelet）的宽度来获得信号的频率特征，通过平移母小波来获得信号的时间信息。对母小波的缩放和平移操作是为了计算小波系数，这些小波系数反映了小波和局部信号之间的相关程度。,4,（a）正弦波曲线；(b)小波曲线,5,6,从小波和正弦波的形状可以看出，变化剧烈的信号，用不规则的小波进行分析比用平滑的正弦波更好，即用小波更能描述信号的局部特征。连续小波变换（ContinuousWaveletTransform，CWT）用下式表示：,(1.1),式（1.1）表示小波变换是信号f(x)与被缩放和平移的小波函数()之积在信号存在的整个期间里求和的结果。CWT的变换结果是许多小波系数C，这些系数是缩放因子（scale）和平移（positon）的函数。,7,基本小波函数()的缩放和平移操作含义如下：(1)缩放。简单地讲，缩放就是压缩或伸展基本小波，缩放系数越小，则小波越窄，如图1.2所示。,图1.2小波的缩放操作,8,(2)平移。简单地讲，平移就是小波的延迟或超前。在数学上，函数f(t)延迟k的表达式为f(t-k)，如图1.3所示。,图1.3小波的平移操作(a)小波函数(t)；(b)位移后的小波函数(t-k),9,图1.4计算系数值C,10,图1.5计算平移后系数值C,11,图1.6计算尺度后系数值C,12,图1.7小波分解示意图,13,图1.12多层小波重构示意图,14,小波的时间和频率特性,运用小波基，可以提取信号中的“指定时间”和“指定频率”的变化。时间：提取信号中“指定时间”（时间A或时间B）的变化。顾名思义，小波在某时间发生的小的波动。频率：提取信号中时间A的比较慢速变化，称较低频率成分；而提取信号中时间B的比较快速变化，称较高频率成分。,时间A,时间B,15,多分辨度分析（MRA）,1988年Mallat提出的多分辨度分析理论，统一了几个不相关的领域：包括语音识别中的镜向滤波，图象处理中的金字塔方法，地震分析中短时波形处理等。当在某一个分辨度检测不到的现象，在另一个分辨度却很容易观察处理。例如：,16,17,参考：M.Vetterli,”WaveletsandSubbandCoding“,PrenticeHallPTR,1995p.11,18,小波的3个特点,小波变换，既具有频率分析的性质，又能表示发生的时间。有利于分析确定时间发生的现象。（傅里叶变换只具有频率分析的性质）小波变换的多分辨度的变换，有利于各分辨度不同特征的提取（图象压缩，边缘抽取，噪声过滤等）小波变换比快速Fourier变换还要快一个数量级。信号长度为M时，Fourier变换（左）和小波变换（右）计算复杂性分别如下公式：,19,小波基表示发生的时间和频率,“时频局域性”图解：Fourier变换的基（上）小波变换基（中）和时间采样基（下）的比较,傅里叶变换(Fourier)基小波基时间采样基,20,Haar小波基母函数,（a）Haar“近似”基函数（b）Haar“细节”基函数低频滤波系数高频滤波系数H0=11qH1=1-1q=qq=q-q其中：,21,Haar小波的基函数,第1行基函数是取平均（近似），第2-8行基函数是取变化（细节）。细节包括变化速率和发生的时间。,H0=11qH1=1-1q,尺度函数近似基函数,小波函数细节基函数,22,小波分析发展历史,1807年Fourier提出傅里叶分析，1822年发表“热传导解析理论”论文1910年Haar提出最简单的小波1980年Morlet首先提出平移伸缩的小波公式，用于地质勘探。1985年Meyer和稍后的Daubeichies提出“正交小波基”，此后形成小波研究的高潮。1988年Mallat提出的多分辨度分析理论（MRA），统一了语音识别中的镜向滤波，子带编码，图象处理中的金字塔法等几个不相关的领域。,23,小波基可以通过给定滤波系数生成,小波基（尺度函数和小波函数）可以通过给定滤波系数生成。有的小波基是正交的，有的是非正交的。有的小波基是对称的，有的是非对称的。小波的近似系数和细节系数可以通过滤波系数直接导出，而不需要确切知道小波基函数，这是I.Daubechies等的重要发现，使计算简化，是快速小波分解和重建的基础。,24,小波基函数和滤波系数(Haar-正交，对称),“近似”基函数,“反变换”低频和高频“滤波系数”,“细节”基函数,Haar小波,“正变换”低频和高频“滤波系数”,25,小波基函数和滤波系数(db2-正交，不对称),“近似”基函数,“细节”基函数,db小波,“反变换”低频和高频“滤波系数”,“正变换”低频和高频“滤波系数”,26,小波基函数和滤波系数(db4-正交，不对称),27,小波基函数和滤波系数(sym4-正交，近似对称),28,小波基函数和滤波系数(bior2.4双正交，对称),29,小波基函数和滤波系数(bior6.8双正交，对称),30,2、小波分析在一维信号处理中的应用,小波变换就是将“原始信号s”变换成“小波系数w”，w=wa,wd包括近似(approximation)系数wa与细节(detail)系数wd近似系数wa-平均成分（低频）细节系数wd-变化成分（高频）,31,小波原始信号分解过程：,原始信号s可分解成小波近似a与小波细节d之和。s=a+d小波系数w=wa,wd的分量，乘以基函数，形成小波分解：小波近似系数wa基函数A=近似分解a-平均小波细节系数wd基函数D=细节分解d-变化,32,小波分解和小波基,正变换：原始信号在小波基上，获得“小波系数”分量反变换：所有“小波分解”合成原始信号例如：小波分解a=小波系数wa小波基A,33,离散小波变换公式,正变换反变换其中：是小波基函数参考“数字图象处理”英文版,电子工业出版社,2002年（R.C.Gonzalaz,”DigitalImageProcessing”，p.375),信号s有M个样本，J级小波变换：,小波分解,小波系数,34,一维信号小波变换例子,Haar小波，例子：16点信号:6598378565981339通过MATLAB实现(wavemenu)波形图小波正变换：小波系数：小波近似系数（加）；小波细节系数（减）小波反变换：可以由分解信号恢复原始信号。有2种：近似分解；细节分解,35,一维信号的二级小波变换系数,原始信号2级小波系数w2=wa2,wd2,wd1*Haar是正交变换。除以常数，目的使变换后平方和不变。例如：,16位,2级近似系数,2级细节系数,1级细节系数,16位,36,一维信号的二级小波变换分解,2级近似分解（原始信号每4个平均值）2级细节分解（原始信号每2个平均的差值）1级细节分解（原始信号单数和双数的差值）恢复信号,37,一维信号的二级小波变换系数和分解,原始信号2级小波系数w2=wa2,wd2,wd12级近似分解（原始信号每4个平均值）2级细节分解（原始信号每2个平均的差值）1级细节分解（原始信号单数和双数的差值）恢复信号,38,原始信号16点,39,两级小波系数16点,原始信号小波系数,40,16点信号的Haar小波近似值和细节分解,两级分解,41,小波去噪声,一般噪声特点：（1）高频成分（细节），（2）幅度小：用阈值；去噪声过程：去除原始信号高频成分（细节）中幅度小于阈值部分。对2级小波，设定2个阈值，称“阈值2”和“阈值1”。去除1级噪声：去除1级小波细节分解中小于“阈值1”部分。去除2级噪声：去除2级小波细节分解中小于“阈值2”部分。恢复：将小波近似分解，加上去噪声后小波细节分解，即获得去除噪声的信号,42,噪声去除,两级分解噪声去除，括号内保留部分数据,43,小波去噪声16点6598378565981339,原始信号(红),去噪后(黄),44,Haar小波去噪声（16点信号）,16点原始信号6598378565981339,小波去噪声,两级分解,45,一维信号的小波变换例子2（电压曲线）,通过MATLAB实现(wavemenu)波形图(MATLABtoolboxwaveletwavedemoleleccum.mat)是“电网监视的电压曲线”，有4570个点Haar小波变换,46,haar小波(s=a2+d2+d1)(wavemenu)leleccumLevel2(s-原始信号，a2-近似，d1-d2细节),47,48,图-5haar(s=a5+d5+.+d1)(wavemenu)leleccumLevel5a5-近似，d5-d1细节,附录-5(wavemenu)leleccumhaarLevel5,leleccum.mat是有36560个点的一维电压信号(s-原始信号，a1-近似，d1-细节),信号前2和后2的差-细节2,信号奇偶数值的差-细节1,原始信号,信号-近似值,5级小波分解,49,小波去噪声leleccumhaar小波,两级小波系数,1级细节小波系数,2级细节小波系数,黄虚线表示阈值,wd1,wd2,原始信号(红),去噪后(黄),50,小波压缩leleccumhaar,51,小波压缩效果leleccumhaar,52,3、小波分析在图象处理中的应用,图象是二维信号，其小波变换相当于二次一维信号的小波变换：。（1）第一次一维信号的小波变换相当于图象的行变换。（2）第二次一维信号的小波变换相当于图象的列变换。小波变换用于图象压缩有良好的效果，已形成图象压缩的标准如JPEG2000。,53,小波变换用于图象特征抽取,54,第1级L1斜线细节,第1级L1水平细节,第1级L1垂直细节,第2级L2细节,近似图象,第3级L3,小波系数分级方块表示法,55,第3级L3分辨率,第2级L2分辨率,第1级L1分辨率,小波系数分级树形表示法,56,小波变换用于图象压缩,采用小波进行压缩。作“小波变换”后，统计特性有改善，消除行和列之间的相关关系。有损压缩：根据视觉原理，不同分辨率小波系数进行比特分配。然后转换到一维作熵编码，如算术编码或霍夫曼编码。无损压缩：选择“整数小波变换”，无舍入误差。但不能进行比特分配。,57,小波变换用于图象压缩,第3级L3水平、斜线、垂直细节,第2级L2水平、斜线、垂直细节,第1级L1水平、斜线、垂直细节,两阈值线之间的直方图被去除（有损压缩）,58,小波变换用于无损数据隐藏,无损数据隐藏：是基于无损压缩：选择“整数小波变换”，无舍入误差。例如可以采用第二代小波。无损数据隐藏：避免在嵌入数据后小波反变换时图象灰度的溢出。小波变换前要作预处理，作直方图调整，将图象中灰度出现少的数据，合并入隐藏数据。第一个无损数据隐藏是1999年科达公司发表的一个专利。由于法律上原因，医学图象数据隐藏必须是无损的。此外、无损数据隐藏在电子银行、电子政务、电子商务、图象建档等有广泛的用途。,59,数据嵌入核磁共振医学图象(可无损恢复)(水印图象见下页)(a)原始(5125128)(b)小波域嵌入水印图象,60,水印图象(1921202二值图象),61,小波变换用于无损数据隐藏（交通图象）,原始图象（1024768）信息隐藏后的伪装图象（1024768）,同时隐藏5张（320280）图象（见下页）,62,同时隐藏的5张（320280）交通图象，可完全恢复,63,小波变换用于图象水印,指纹原始图象嵌入水印(取款密码等)后图象,指纹传感器：标准的Veridicom指纹鼠标指纹开发工具：VeridicomAuthenticationSDK以Windows的DLL库方式提供指纹库：(FingerprintVerificationCompetition,FVC)。FVC2000db1是由光学设备采集；FVC2000db2是由电容设备采集。,银行取款密码嵌入指纹，网上进行身份认证,64,小波变换用于图象水印,小波正变换,小波反变换,小波正变换,小波反变换,数据嵌入,数据提取,原始图象,加水印后图象输入,原始图象,加水印后图象输出,隐藏数据,隐藏数据,65,小波分析最新进展,（1）第二代小波，称提升算