数字图像变换

1、第十章数字图像变换 数字图像处理的方法分为两类：空间域处理法和频域法。 空间域：计算复杂、费时，甚至难以实现。 频域：运算速度高、可用滤波技术简化运算。 图像变换可以将图像从空间域转换到频率域，然后在频率域对图像进行各种处理，再将所得到的结果进行反变换，即从频率域变换到空间域，从而达到图像处理的目的。第十章图像变换 图像处理中应用正交变换，进行如图像增强、复原、编码、描述和特征提取。 正交变换 ： 正弦型变换：正弦型变换：傅里叶变换、余弦变换和正弦变换。 方波型变换：方波型变换：哈达玛（Hadamarn）变换、沃尔什（ Walsh）变换、斜变换、小波变换。 基于特征向量的变换基于特征向量的变换

2、：主要包括Hotelling变换、KL变换和SVD变换。第十章 图像变换 10.1 10.1 频域世界与频域变换频域世界与频域变换10.2 10.2 傅立叶变换傅立叶变换10.3 10.3 离散余弦变换离散余弦变换 10.4 10.4 离散沃尔什哈达玛变换离散沃尔什哈达玛变换 10.5 10.5 小波变换简介小波变换简介 10.1 10.1 频域世界与频域变换频域世界与频域变换任意波形可分解为正弦波的加权和任意波形可分解为正弦波的加权和 y1 = Sin(x + /2) A=1, = /2, f=1/ 2 y2=0.5sin(2x+ ) A=0.5, = , f=1/ y3=0.25sin(4

3、x+ 3 /2) A=0.25, = 3 /2 , f=2/ y= Sin(x + /2) + 0.5sin(2x+ ) + 0.25sin(4x+ 3 /2) x 0,4 波形的频域表示波形的频域表示y= Sin(x + /2) + 0.5sin(2x+ ) + 0.25sin(4x+ 3 /2) x 0,4 f=w/2 幅频特性幅频特性Af0.250.510.751/2 3/2 1/ 2/ 相频特性相频特性f /2 2 3 /21/2 3/2 1/ 2/ 10.110.1频域世界与频域变换频域世界与频域变换幅频特性幅频特性Af0.250.510.751/2 3/2 1/ 2/ 相频特性相频

4、特性f /2 2 3 /21/2 3/2 1/ 2/ iiiff4log2)(2ifiiiffA212)(20)2sin()(iiiixfAxy10.1 10.1 频域世界与频域变换频域世界与频域变换710.2傅里叶变换: 112T直流分量基波分量n =1 谐波分量n11n)sincos()(11101tnbtnaatfnnn一：周期函数的傅里叶变换：一：周期函数的傅里叶变换：8100.cos).(211TttndttntfTadttntfTbTttn.sin).(210011直流系数余弦分量系数正弦分量系数10.2傅里叶变换:100).(110TttdttfTa1：周期函数的频谱分析910.

5、2傅里叶变换:2：周期函数的复指数级数l由前知l由欧拉公式l其中)sincos()(11101tnbtnaatfnnntjnnenFtf1)()(1)(21)(1nnjbanF)(21)(1nnjbanF0)0(aF引入了负频率10.2傅里叶变换:二：非周期函数傅立叶变换二：非周期函数傅立叶变换分析分析式：式：()()j w tFwfted t deFtftj. )(21)(10.2傅里叶变换: 假定以间隔x对一个连续函数f(x)均匀采样，离散化为一个序列 f(x0), f(x0+x), fx0+(N1)x（如图3.3所示），则将序列表示 f(x)=f(x0+x) 式中x假定为离散值0，1，2

6、，N1。换句话说，序列 f(0),f(1),f(2),f(N1) 表示取自该连续函数N个等间隔的抽样值。三：离散函数的傅里叶变换三：离散函数的傅里叶变换.1离散傅立叶变换离散傅立叶变换被抽样函数的离散傅里叶变换定义为反变换为在二维的情况下，离散的傅里叶变换对表示为NuxjNxexfuF/210)()(NuxjNueuFNxf/210)(1)(1010)/(2),(),(MxNyNvyMuxjeyxfvuF1010)/(2),(1),(MuNvNvyMuxjeyvuFMNyxf二维离散函数的傅立叶频谱、二维离散函数的傅立叶频谱、 相位谱和能量谱分别为相位谱和能量谱分别为 ),

7、(),(),(),(),(arctan),(),(),(| ),(|2222vuIvuRvuEvuRvuIvuvuIvuRvuF式中，式中，R R( (u u, , v v) )和和I I( (u u, , v v) )分别是分别是F F( (u u, , v v) )的实部和虚部。的实部和虚部。 .1离散傅立叶变换离散傅立叶变换10.2.2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换 离散傅立叶变换计算量非常大，运算时间长。研究离散傅离散傅立叶变换计算量非常大，运算时间长。研究离散傅立叶变换的快速算法（立叶变换的快速算法（Fast Fourier Transfor

8、mFast Fourier Transform， FFTFFT）是）是非常有必要的。非常有必要的。 介绍一种称为逐次加倍法的快速傅立叶变换算法（介绍一种称为逐次加倍法的快速傅立叶变换算法（FFTFFT），），它是它是19651965年年CooleyCooley和和TukeyTukey首先提出的。首先提出的。 二维离散傅立叶变换具有可分离性，即它可由两次一维二维离散傅立叶变换具有可分离性，即它可由两次一维离散傅立叶变换计算得到，因此，仅研究一维离散傅立叶离散傅立叶变换计算得到，因此，仅研究一维离散傅立叶变换的快速算法即可。改写公式：变换的快速算法即可。改写公式： 10)()(NxuxNWxfuF

9、式中，式中， =e=e-j2-j2N N ，称为旋转因子。，称为旋转因子。 = e e-j2-j2N N = =cos(22N N )-j )-j sin(22N N ) ( ) (以以N N为周期为周期) )式中很多式中很多 系数相同，不必进行多次重复计算。系数相同，不必进行多次重复计算。10.2.2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换uxWNWNWNWNWNW FFT FFT的推导过程：的推导过程： 设设N N为为2 2的正整数次幂，的正整数次幂， 即即 , 2 , 12nNn令令M M=N/2,=N/2,离散傅立叶变换可改写成如下形式：离散傅立叶变换可改写成如下形式： 1

10、0)12(2)2(2101202) 12()2()()(MxxuMxuMMxMxuxMWxfWxfWxfuF偶离散点偶离散点奇离散点奇离散点10.2.2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换uxMuxMjuxMjuxMWeeW)()(/222/222uMuxMMxMxuxMWWxfWxfuF21010) 12()2()( 定义定义 1, 1 ,0,)12()(1, 1 ,0,)2()(1010MxuWxfuFMxuWxfuFMxuxMoMxuxMe10)12(2)2(2101202) 12()2()()(MxxuMxuMMxMxuxMWxfWxfWxfuF10.2.2 10.2.

11、2 快速离散傅立叶变换快速离散傅立叶变换even 偶odd奇于是于是)()()(2uFWuFuFouMe 将一个将一个N N点的离散傅立叶变换分解成两个点的离散傅立叶变换分解成两个N N2 2短序列的离短序列的离散傅立叶变换，即分解为偶数和奇数序列的离散傅立叶变散傅立叶变换，即分解为偶数和奇数序列的离散傅立叶变换换F Fe e( (u u) )和和F Fo o( (u u) ) 。 )7()7()7()6()6()6()5()5()5()4()4()4()3()3()3()2()2()2() 1 () 1 () 1 ()0()0()0(7868584838281808oeoeoeoeoeoeo

12、eoeFWFFFWFFFWFFFWFFFWFFFWFFFWFFFWFF设设N=2N=23 310.2.2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换) 7 () 7 () 7 () 6 () 6 () 6 () 5 () 5 () 5 () 4() 4() 4() 3 () 3 () 3 () 2() 2() 2() 1 () 1 () 1 () 0 () 0 () 0 (7868584838281808oeoeoeoeoeoeoeoeFWFFFWFFFWFFFWFFFWFFFWFFFWFFFWFF)()()(2MuFWMuFMuFoMuMeMxxMuMMMuMMxxMuMWxfWWW

13、xf0)(220)()2()2(MxMxMuxMuMMxMxMuxMWWxfWWWxf020) 12()2(MxuxMuMMxuxMWxfWWxf020) 12()2()()(2uFWuFouMe)0()0()40(2ouMeFWFF)2()2()42(2ouMeFWFF) 1 () 1 ()41 (2ouMeFWFF)3()3()43(2ouMeFWFF10.2.2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换蝶形运算单元蝶形运算单元 Fe(1)F(1)F(5)Fo(1)18W18W) 1 () 1 ()5() 1 () 1 () 1 (1818oeoeFWFFFWFF)3()3()7

14、()2()2()6() 1 () 1 ()5()0()0()4()3()3()3()2()2()2() 1 () 1 () 1 ()0()0()0(3828180838281808oeoeoeoeoeoeoeoeFWFFFWFFFWFFFWFFFWFFFWFFFWFFFWFF10.2.2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换Fe(0)Fe(1)Fe(2)Fe(3)Fo(0)Fo(1)Fo(2)Fo(3)08W18W28W38W08W18W28W38WF(0)F(1)F(2)F(3)F(4)F(5)F(6)F(7) 3 () 3 () 7() 2() 2() 6() 1 () 1

15、() 5 () 0() 0() 4() 3 () 3 () 3 () 2() 2() 2() 1 () 1 () 1 () 0() 0() 0(3828180838281808oeoeoeoeoeoeoeoeFWFFFWFFFWFFFWFFFWFFFWFFFWFFFWFF10.2.2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换 F Fe e( (u u) )和和F Fo o( (u u) )都是都是4 4点的点的DFTDFT，对它们再按照奇偶进行分组，对它们再按照奇偶进行分组) 1 () 1 () 3()0()0()2() 1 () 1 () 1 ()0()0()0(28082808e

16、oeeeeoeeeeoeeeeoeeeFWFFFWFFFWFFFWFF) 1 () 1 ()3()0()0()2() 1 () 1 () 1 ()0()0()0(28082808oooeooooeooooeooooeoFWFFFWFFFWFFFWFFFee(0)Feo(1)08W28WFee(1)Feo(0)Fe(0)Fe(1)Fe(2)Fe(3)28W08WFoe(0)Foo(1)08W28WFoe(1)Foo(0)Fo(0)Fo(1)Fo(2)Fo(3)28W08W10.2.2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换8点点DFT的蝶形流程图的蝶形流程图 Fee(0)Feo(1)

17、08W28WFee(1)Feo(0)Fe(0)Fe(1)Fe(2)Fe(3)28W08WFoe(0)Foo(1)08W28WFoe(1)Foo(0)Fo(0)Fo(1)Fo(2)Fo(3)28W08Wf (0)f (4)08W08Wf (2)f (6)08W08Wf (1)f (5)08W08Wf (3)f (7)08W08W08W18W28W38W08W18W28W38WF(0)F(1)F(2)F(3)F(4)F(5)F(6)F(7)例：例：0 1 0 2 0 3 0 40 5 0 60 7 0 8Fe(0)Fo(1)04W14WFe(1)Fo(0)F (0)F (1)F (2)F (3)1

18、4W04W04W04W04W04Wf(0)f(2)f(1)f(3)10.2.2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换0 1 0 2 0 3 0 40 5 0 60 7 0 83 i -3 -i 0 3 0 40 5 0 60 7 0 80012003-13i-3-i i-i 1 1 110.2.2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换0034007-17i-7-i i-i 1 1 13 i -3 -i 7 i -7 -i0 5 0 60 7 0 83 i -3 -i 0 3 0 40 5 0 60 7 0 810.2.2 10.2.2 快速离散傅立叶变换快速离散傅立

19、叶变换00560011-111i-11-i i-i 1 1 13 i -3 -i 7 i -7 -i0 5 0 60 7 0 83 i -3 -i 7 i -7 -i11 i -11 -i0 7 0 810.2.2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换00780015-115i-15-i i-i 1 1 13 i -3 -i 7 i -7 -i11 i -11 -i0 7 0 83 i -3 -i 7 i -7 -i11 i -11 -i15 i -15 -i10.2.2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换31171514-822-836-8+8i-8-8-8i

20、i-i 1 1 13 i -3 -i 7 i -7 -i11 i -11 -i15 i -15 -i36 i -3 -i -8+8i i -7 -i-8 i -11 -i-8-8i i -15 -i10.2.2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换iiii2i02i04i000 i-i 1 1 136 i -3 -i -8+8i i -7 -i-8 i -11 -i-8-8i i -15 -i36 4i -3 -i -8+8i 0 -7 -i-8 0 -11 -i-8-8i 0 -15 -i10.2.2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换-3-11-7-15-14

21、8-228-368-8i88+8i i-i 1 1 136 4i -3 -i -8+8i 0 -7 -i-8 0 -11 -i-8-8i 0 -15 -i36 4i -36 -i -8+8i 0 8-8i -i-8 0 8 -i-8-8i 0 8+8i -i10.2.2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换-i-i-i-i-2i0-2i0-4i000 i-i 1 1 136 4i -36 -i -8+8i 0 8-8i -i-8 0 8 -i-8-8i 0 8+8i -i36 4i -36 -4i -8+8i 0 8-8i 0-8 0 8 0-8-8i 0 8+8i 010.2.

22、2 10.2.2 快速离散傅立叶变换快速离散傅立叶变换10.3 10.3 离散余弦变换（离散余弦变换（DCTDCT） 离散余弦变换（离散余弦变换（Discrete CosineDiscrete Cosine Transform Transform， DCTDCT）的变换核为余弦）的变换核为余弦函数。函数。DCTDCT变换被认为是一种语音信号变换被认为是一种语音信号、图像信号的变换的准最佳变换。、图像信号的变换的准最佳变换。 10.3.1 10.3.1 一维离散余弦变换一维离散余弦变换 一维一维DCTDCT定义如下：定义如下： 设设 f f( (x x)|)|x x=0,1,=0,1, ,N N

23、-1-1为离散的信号列。为离散的信号列。 102) 12(cos)(2)()(NxNuxxfNuCuF102) 12(cos)()(2)(NuNuxuFuCNxfu u, ,x x=0,1,2,=0,1,2, ,N N1 1其他1021)(uuC见课本P200FCffCFT一维离散余弦变换10.3.2 10.3.2 二维离散余弦变换二维离散余弦变换 二维二维DCTDCT定义如下：设定义如下：设f f( (x, yx, y) )为为M MN N的数字图像矩阵，则的数字图像矩阵，则 NvyMuxvCuCyxfMNvuFMxNy2) 12(cos2) 12(cos)()(),(2),(1010Nvy

24、MuxvuFvCuCMNyxfMuNv2) 12(cos2) 12(cos),()()(2),(1010 x x, ,u u=0,1,2,=0,1,2, ,M M1 1 y y, ,v v=0,1,2,=0,1,2, ,N N1 1C C( (u u) )和和C C( (v v) )的定义同前的定义同前完成p203例题10.410.4离散沃尔什离散沃尔什- -哈达玛变换（哈达玛变换（WHTWHT） 10.4.1 10.4.1 一维离散沃尔什一维离散沃尔什- -哈达玛变换哈达玛变换 1. 1. 沃尔什函数沃尔什函数 沃尔什函数是沃尔什函数是19231923年由美国数学家沃尔什提出的。年由美国数学

25、家沃尔什提出的。它是一个完备正交函数系，其值只能取它是一个完备正交函数系，其值只能取1 1和和1 1。从排列次序上可将沃尔什函数分为三种定义方法。从排列次序上可将沃尔什函数分为三种定义方法。在此只介绍哈达玛排列定义的沃尔什变换。在此只介绍哈达玛排列定义的沃尔什变换。2n 阶哈达玛矩阵有如下形式：阶哈达玛矩阵有如下形式： 1111 1 21HH111111111111111122224HHHHH10.4.1 10.4.1 一维离散沃尔什一维离散沃尔什- -哈达玛变换哈达玛变换2. 2. 离散沃尔什离散沃尔什- -哈达玛变换哈达玛变换一维离散沃尔什变换及逆变换定义为一维离散沃尔什变换及逆变换定义为

26、 若将若将Walsh(Walsh(u u, , x x) )用哈达玛矩阵表示，并将变换表达式写用哈达玛矩阵表示，并将变换表达式写成矩阵形式，则上两式分别为：成矩阵形式，则上两式分别为： 10),()(1)(NxxuWalshxfNuW10),()()(NuxuWalshuWxf10.4.1 10.4.1 一维离散沃尔什一维离散沃尔什- -哈达玛变换哈达玛变换) 1() 1 ()0(1) 1() 1 ()0(NfffHNNWWWN) 1() 1 ()0() 1() 1 ()0(NWWWHNfffNHN为为N阶阶哈达玛哈达玛矩阵矩阵 由哈达玛矩阵的由哈达玛矩阵的特点可知，沃尔什特点可知，沃尔什-

27、-哈达玛变换的本质上哈达玛变换的本质上是将离散序列是将离散序列f f( (x x) )的的各项值的符号按一定各项值的符号按一定规律改变后，进行加规律改变后，进行加减运算，它比采用复减运算，它比采用复数运算的数运算的DFTDFT和采用和采用余弦运算的余弦运算的DCTDCT要简要简单得多。单得多。 课本P208例子10.4.2 10.4.2 二维离散沃尔什变换二维离散沃尔什变换 二维二维WHTWHT的正变换核和逆变换分别为的正变换核和逆变换分别为 1010),(),(),(1),(NyMxyvWslshxuWalshyxfMNvuW1010),(),(),(),(NvMuyvWslshxuWals

28、hvuWyxfx x, ,u u=0,1,2,=0,1,2, ,M M1 1 y y, ,v v=0,1,2,=0,1,2, ,N N1 1例有两个二维数字图像信号矩阵如下，求这两个信号的二维例有两个二维数字图像信号矩阵如下，求这两个信号的二维WHTWHT。 13311331133113311f11111111111111112f根据题意，根据题意，M M= =N N=4=4，其二维，其二维WHTWHT变换核为变换核为 11111111111111114H10.4.2 10.4.2 二维离散沃尔什变换二维离散沃尔什变换00000000000010021111111111111111133113

29、311331133111111111111111114121W00000000000000011111111111111111111111111111111111111111111111114122W 从以上例子可看出，二维从以上例子可看出，二维WHTWHT具有能量集中的特性，而且具有能量集中的特性，而且原始数据中数字越是均匀分布，经变换后的数据越集中于矩阵原始数据中数字越是均匀分布，经变换后的数据越集中于矩阵的边角上。因此，二维的边角上。因此，二维WHTWHT可用于压缩图像信息。可用于压缩图像信息。 10.510.5小波变换简介小波变换简介 与傅立叶变换不同，小波变换是通过缩放母小波与傅立叶

30、变换不同，小波变换是通过缩放母小波（Mother waveletMother wavelet）的宽度来获得信号的频率特征，通过）的宽度来获得信号的频率特征，通过平移母小波来获得信号的时间信息。对母小波的缩放和平移母小波来获得信号的时间信息。对母小波的缩放和平移操作是为了计算小波系数，这些小波系数反映了小平移操作是为了计算小波系数，这些小波系数反映了小波和局部信号之间的相关程度。波和局部信号之间的相关程度。 1. 1. 连续小波变换（连续小波变换（CWTCWT） 小波分析就是把一个信号分解为将母小波经过缩放和小波分析就是把一个信号分解为将母小波经过缩放和平移之后的一系列小波。小波变换可以理解为用

31、经过缩平移之后的一系列小波。小波变换可以理解为用经过缩放和平移的一系列小波函数代替傅立叶变换的正弦波和放和平移的一系列小波函数代替傅立叶变换的正弦波和余弦波进行傅立叶变换的结果。余弦波进行傅立叶变换的结果。(a)(b) 从小波和正弦波的形状可以看出，变化剧烈的信号，从小波和正弦波的形状可以看出，变化剧烈的信号，用不规则的小波进行分析比用平滑的正弦波更好，用小波用不规则的小波进行分析比用平滑的正弦波更好，用小波更能描述信号的局部特征。更能描述信号的局部特征。 连续小波变换（连续小波变换（Continuous Wavelet TransformContinuous Wavelet Transfor

32、m， CWTCWT）用下式表示：用下式表示： dttpositionscaletfpositionscaleC),()(),( CWT CWT的变换结果是许多小波系数的变换结果是许多小波系数C C，这些系数是缩放因，这些系数是缩放因子（子（scalescale）和平移（）和平移（positonpositon）的函数。）的函数。 1. 1. 连续小波变换（连续小波变换（CWTCWT） (1) (1) 缩放就是压缩或伸展基波，缩放系数越小，则小波缩放就是压缩或伸展基波，缩放系数越小，则小波越窄。越窄。 小波的缩放操作小波的缩放操作 OOOf (t)f (t)f (t)tttf (t)(t)； sc

33、ale 1f (t)(2t)； scale 0.5f (t)(4t)； scale 0.251. 1. 连续小波变换（连续小波变换（CWTCWT） (2) (2) 平移就是小波的延迟或超前。在数学上，函数平移就是小波的延迟或超前。在数学上，函数f f( (t t) )延迟延迟k k的表达式为的表达式为f f( (t-kt-k) ) 。 小波的平移操作小波的平移操作(a) (a) 小波函数小波函数( (t t) )； ( (b b) ) 位移后的小波函数位移后的小波函数( (t-kt-k) ) Ot(t)Ot(t k)(a)(b)1. 1. 连续小波变换（连续小波变换（CWTCWT） CWT C

34、WT计算主要有如下五个步骤：计算主要有如下五个步骤： 第一步：第一步： 取一个小波，取一个小波， 将其与原始信号的开始一节进将其与原始信号的开始一节进行比较。行比较。 第二步：第二步： 计算小波与所取一节信号的相似程度计算小波与所取一节信号的相似程度C C，计算，计算结果取决于所选小波的形状。结果取决于所选小波的形状。 1. 1. 连续小波变换（连续小波变换（CWTCWT）原 始 信 号小 波 信 号第三步：向右移动小波，重复第一步和第二步，直至覆盖整第三步：向右移动小波，重复第一步和第二步，直至覆盖整个信号。个信号。1. 1. 连续小波变换（连续小波变换（CWTCWT）计算尺度后系数值计算尺

35、度后系数值C C 原 始 信 号小 波 信 号C 0.2247第四步：第四步： 伸展小波，伸展小波， 重复第一步至第三步。重复第一步至第三步。1. 1. 连续小波变换（连续小波变换（CWTCWT） 第五步：对于所有缩放，重复第一步至第四步。第五步：对于所有缩放，重复第一步至第四步。 缩放因子缩放因子scalescale越小，小波越窄，度量的是信号的细节越小，小波越窄，度量的是信号的细节变化，表示信号频率越高；缩放因子变化，表示信号频率越高；缩放因子scalescale越大，小波越宽，越大，小波越宽，度量的是信号的粗糙程度，表示信号频率越低。度量的是信号的粗糙程度，表示信号频率越低。 1. 1.

36、 连续小波变换（连续小波变换（CWTCWT） 2. 2. 离散小波变换（离散小波变换（DWTDWT） 在每个可能的缩放因子和平移参数下计算小波系数，其在每个可能的缩放因子和平移参数下计算小波系数，其计算量相当大，而且有许多数据是无用的。如果缩放因子和计算量相当大，而且有许多数据是无用的。如果缩放因子和平移参数都选择为平移参数都选择为2 2j j（j j00且为整数）的倍数，就会使分析且为整数）的倍数，就会使分析的数据量大大减少。使用这样的缩放因子和平移参数的小的数据量大大减少。使用这样的缩放因子和平移参数的小波变换称为波变换称为双尺度小波变换双尺度小波变换（Dyadic Wavelet Tra

37、nsformDyadic Wavelet Transform），），它是离散小波变换（它是离散小波变换（Discrete Wavelet TransformDiscrete Wavelet Transform， DWTDWT）的一种形式。通常离散小波变换就是指双尺度小波变换。的一种形式。通常离散小波变换就是指双尺度小波变换。 执行离散小波变换的有效方法是使用滤波器，它是一种执行离散小波变换的有效方法是使用滤波器，它是一种信号分解的方法，信号分解的方法， 又常称为又常称为双通道子带编码双通道子带编码。 SAD滤波器组低通高通小波分解示意图小波分解示意图2. 2. 离散小波变换（离散小波变换（DW

38、TDWT） 信号的低频分量是最重要的，而高频分量只起一个修饰信号的低频分量是最重要的，而高频分量只起一个修饰的作用。如同一个人的声音一样，把高频分量去掉后，听起的作用。如同一个人的声音一样，把高频分量去掉后，听起来声音会发生改变，但还能听出说的是什么内容，但如果来声音会发生改变，但还能听出说的是什么内容，但如果把低频分量删除后，就会什么内容也听不出来了。把低频分量删除后，就会什么内容也听不出来了。 2. 2. 离散小波变换（离散小波变换（DWTDWT）多级信号分解示意图多级信号分解示意图（a） 信号分解；信号分解； (b) 小波分数；小波分数； （c）小波分解树）小波分解树 cA3cD3cA2

39、cD2SLo_DHi_DA1D1Lo_DHi_DA2D2Lo_DHi_DA3D3Lo_D：低通滤波器；Hi_D：高通滤波器(a)ScA1cD1(b)(c)ScA1cD1cA2cD2cA3cD3一级分解一级分解对低频分量连续分解对低频分量连续分解,可得到信号不同分辨率下的可得到信号不同分辨率下的低频分量低频分量,也称为信号的多分也称为信号的多分辨率分析辨率分析分解的级数取决于分解的级数取决于要分析的信号数据特征要分析的信号数据特征及用户的具体需要。及用户的具体需要。 SDA1000个采样点1000个采样点1000个采样点ScDcA1000个采样点约500个DWT系数约500个DWT系数表示下采样

40、表示下采样 对于一个信号，如采用上述方法，理论上产生的数据量对于一个信号，如采用上述方法，理论上产生的数据量将是原始数据的两倍。于是，根据奈奎斯特（将是原始数据的两倍。于是，根据奈奎斯特（NyquistNyquist）采样）采样定理，定理， 可用可用下采样下采样的方法来减少数据量，即在每个通道内每的方法来减少数据量，即在每个通道内每两个样本数据取一个，便可得到离散小波变换的系数两个样本数据取一个，便可得到离散小波变换的系数（CoefficientCoefficient），）， 分别用分别用cAcA和和cDcD表示。表示。 2. 2. 离散小波变换（离散小波变换（DWTDWT） 3. 3. 小波

41、重构小波重构 利用信号的小波分解的系数还原出原始信号，这一过程称利用信号的小波分解的系数还原出原始信号，这一过程称为为小波重构小波重构（Wavelet ReconstructionWavelet Reconstruction）或叫）或叫小波合成小波合成（Wavelet SynthesisWavelet Synthesis）。这一合成过程的数学运算叫做）。这一合成过程的数学运算叫做逆离散逆离散小波变换小波变换（Inverse Discrete Wavelet TransformInverse Discrete Wavelet Transform， IDWTIDWT）。）。 SHLHL 1 1）重

42、构近似信号与细节信号）重构近似信号与细节信号 小波分解的近似系数和细节系数可以重构出原始信号。同小波分解的近似系数和细节系数可以重构出原始信号。同样，可由近似系数和细节系数分别重构出信号的近似值或细节样，可由近似系数和细节系数分别重构出信号的近似值或细节值，这时只要近似系数或细节系数置为零即可。值，这时只要近似系数或细节系数置为零即可。 重构近似和细节信号示意（a） 重构近似信号； (b) 重构细节信号 A1HL1000个 样 点0约 500个 0cA1约 500个 近 似 分 量(a)D1HL1000个 样 点(b)约 500个 0约 500个 近 似 分 量0cD1 2 2）多层重构）多层

43、重构 在上图中，重构出信号的近似值在上图中，重构出信号的近似值A A1 1与细节值与细节值D D1 1之后，则之后，则原信号可用原信号可用A A1 1D D1 1S S重构出来。对应于信号的多层小波分解，重构出来。对应于信号的多层小波分解，小波的多层重构如下图：小波的多层重构如下图：A3D3A2D2SA1D1重构过程为：重构过程为：A A3 3+ +D D3 3A A2 2 A A2 2+ +D D2 2A A1 1 A A1 1+ +D D1 1S S 2 2）多层重构）多层重构信号重构中，滤波器的选择非常重要，关系到能否重构出信号重构中，滤波器的选择非常重要，关系到能否重构出满意的原始信号

44、。低通分解滤波器（满意的原始信号。低通分解滤波器（L L）和高通分解滤波器）和高通分解滤波器（H H）及重构滤波器组（）及重构滤波器组（LL和和HH）构成一个系统，这个系）构成一个系统，这个系统称为正交镜像滤波器（统称为正交镜像滤波器（Quadrature Mirror FiltersQuadrature Mirror Filters， QMFQMF）系统。系统。 4. 4. 小波包分析小波包分析 小波分析是将信号分解为近似与细节两部分，近似部分又小波分析是将信号分解为近似与细节两部分，近似部分又可以分解成第二层近似与细节，可以这样重复下去。对于一个可以分解成第二层近似与细节，可以这样重复下去。对于一个N N层分解来说，层分解来说， 有有N N+1+1个分解信号的途径。个分解信号的途径。