《小波分析及其应用》(孙延奎,2005)第5章 可分离二维小波变换算法实现的问题与讨论
孙延奎
摘要:总结分析可分离二维小波变换算法实现的细节,澄清一些概念与问题,纠正例5-2中“转置”的错误,回答读者的问题。
教材中三个方向小波的定义:
??1(x,y)??(x)?(y)?2 ??(x,y)??(x)?(y)
??3(x,y)??(x)?(y)?
1j,22j,33gj(x,y)??dkj,,1m?j,k,m??dk,m?j,k,m??dk,m?j,k,m (5.5)
k,mk,mk,m经简单计算可得,
1?ckj,?m??j,1?d?k,m??j,2?dk,m??j,3??dk,m?f,?j?1,k,mf,?1j,k,mf,?2j,k,m
f,?3j,k,m我们称序列
?c,djj,1,dj,2,dj,3?为cj?1的(一级)二维小波变换。下面讨论二维小波变换
的快速算法。
设一维多分辨分析
?V?的两尺度方程和小波方程为:
j??t??2?hk??2t?k?k??t??2?gk??2t?k?k
其中,
?h?k为实滤波器,gk???1?h1?k。则类似一维正交多分辨分析的推导,由
k1j,11j,22j,33fj?1(x,y)??ckj,?m?j?1,k,m??dk,m?j,k,m??dk,m?j,k,m??dk,m?j,k,m
k,mk,mk,mk,mckj,m?f,?j,k,m??????RR?Rf?x,y??*x,y?dxdyj,k,m???Rf?x,y?2j?*(2jx?k)?*(2jy?m)dxdy????f?x,y?2j?2?hl*?*(2j?1x?2k?l)??2?hn*?*(2j?1y?2m?n)?dxdyln????RRR??hl*hn*?l,n?R2j?1f?x,y??*(2j?1x?2k?l)?*(2j?1y?2m?n)dxdy1**j?1**j?1??hl*hn*c2jk???hhc?hhc??l,2m?n2k?l2m?nl,nl?2kn?2ml,nl,nl,nl,ndj,1k,m?f,?1j,k,m??R?f?x,y???R1j,k,m?x,y??*dxdy????R??Rf?x,y?2j?*(2jx?k)?*(2jy?m)dxdy????**f?x,y?2j?2?hl*?*(2j?1x?2k?l)??2?gn?(2j?1y?2m?n)?dxdyln????RRR*??hl*gn?l,n?R2j?1f?x,y??*(2j?1x?2k?l)?*(2j?1y?2m?n)dxdy*j?1*j?1**j?1??hl*gnc2k?l,2m?n??hl*?2kgnc?hgc?2k?l2m?nl,n?2ml,nl,nl,nl,n
得二维Mallat算法如下(假设h是实滤波器):
?jj?1j?1c?hhc?hhc?k,m?l?2kn?2ml,n?2k?l2m?nl,nl,nl,n??j,1j?1j?1d?hgc?hgc?k,m?l?2kn?2ml,n?2k?l2m?nl,nl,nl,n???dj,2?ghcj?1?ghcj?1?l?2kn?2ml,n2k?l2m?nl,n (5.6)
?k,m?l,nl,n??dj,3?ggj?1j?1c?ggc?l?2kn?2ml,n2k?l2m?nl,n?k,m?l,nl,n?重构算法:
1jj,1j,2j,3ckj,?m??hk?2lhm?2ncl,n??hk?2lgm?2ndl,n??gk?2lhm?2ndl,n??gk?2lgm?2ndl,n(5.7)
l,nl,nl,nl,n
评注:
1) 在关履泰 编著“小波方法与应用”中,空间分解表示与我教材中是一致的,但二维可
分离小波?1与?2的意思正好与我的相反。其?与?j对应我教材中的d理论分析及试验验证,Matlab中的用法与关履泰书中的一致。
jj,2和dj,1。经
2) 在Mallat 编著“信号处理的小波导引”中,空间分解表示及二维可分离小波?1,?2,
?2,dj,1,dj,2等都与教材中的一致。只是最后对变换的结果的表示略有不同:
教材中采用的方法是:
?cj ?j,2?d?cj ?j,1?d如对下面一幅图像,
dj,1? j,3?d?dj,2?? dj,3?而Mallat著作中,最后的表示为
图1
教材中的表示如下右图所示,其中左上角表示低频系数,左下表示垂直边缘; 右上角表示水平边缘;右下角表示对角边缘。
图2
而在Mallat的著作“信号处理的小波导引”中(图7-26)中,各个分辨率下的图像中水平边缘与竖直边缘的位置正好相反。其中,左上角表示低频系数,左下表示水平细节系数(水平边缘);右上角表示垂直细节系数;右下角表示对角细节系数。
可分离双正交基:
可将一维双正交小波基推广到L(R)的可分离正交小波基。令?,?和?,?是生成L(R)的双正交小波基的两对对偶的尺度函数与小波。则
222??1(x,y)??(x)?(y)?2??(x,y)??(x)?(y) ??3(x,y)??(x)?(y)?所定义的?1,?2,?3的对偶小波是:
??1(x,y)??(x)?(y)?2??(x,y)??(x)?(y) ??3(x,y)??(x)?(y)?可以证明:
??1j,n3,?2,?j,nj,n?j,n?Z和?j,n,?j,n,?j,n3?123?j,n?Z3是L2(R2)的双正交Riesz基。
二维可分离小波的频率特性分析:
二维可分离小波?1,?2,?3在不同的尺度和方向上提取图像细节。即用?1和?2计算出的小波系数分别在水平和垂直边缘上取得大的值。分别提取图像的水平与垂直特征。小波?3在角点上产生大的系数,提取对角边缘。具体的,
?1: 垂直高频,提取水平边缘
?2: 水平高频,提取垂直边缘
?和??的能量分别集中在?0,??和??,2??上。由可分离小波?,?,?的在正频率上,?123表达式可推出:
?????????,??1??x,?y??????y?,??2??x,?y??????x???3??x,?y??????x?????y? ?xy?因此,在低的水平频率?x处和高的垂直频率?y处??2?x,?y和低的垂直频率?y处?1??,??大;在高的水平频率?xyx处
??大;而在高的水平频率?x和垂直频率?y处
?3??x,?y?大。可显示出正交小波对应的可分离小波的傅立叶变换图。 ?在王大凯,彭进业编著的“小波分析及其在信号处理中的应用” 中,塔式分解写法中低频子带与三个高频子带的位置正好反映它们对频域的内在划分,注意,该书中?,?,?的定义与教材相同。
即
123? ??2 c d ?jj,2
?1 ??3 d ?j,1 dj,3
若用L表示低通滤波器,用H表示高通滤波器,则滤波器LL,LH,HL和HH构成4个具有
不同频率特性和方向特性的滤波器。 LL用于检索图像中的低频分量,LH用于检测水平方向的边缘、细节分量,HL用于检测垂直方向的边缘、细节分量,HH用于检测主对角与副对角方向的分量。【这里的LH,意思是先用L做行变换,再H做列变换。与教材中的意思相反。】
可见,采用以下表示
?cj ?j,1?ddj,2? j,3?d?能够与它们对频域的划分一致。
其它著作中有关问题的描述:
1. 在杨福生著的“小波变换的工程分析与应用”中(P117),三个小波?1,?2,?3与本教材的是一致的。它指出:三个小波?1,?2,?3中都至少包含一个带通的??x?或??y?,因此,它们都是带通的,也就是说,这三部分反映的都是高通细节。
指出,对f?x,y?,先沿x方向分别用??x?和??x?作分析,把f?x,y?分成平滑逼近和细节这两部分,然后对这两部分再分别用??y?和??y?作类似分析。得到一个平滑逼近与三个细节函数。
注意:该著作中(5.13a), (5.13b)的写法不够严谨。 其中,Ajf,Dj,Dj,Dj的含义同本教材中的c,d?1??2??3??1??2??3?jj,1,dj,2,dj,3。
图5.4中的b,c,d是否与Dj,Dj,Dj相对应,没有指出。是对应的。 图5.6中Dj,Dj的位置搞错了。
2) 李弼程等编著的“小波分析及其应用”,在P40页指出:一幅图像可分解为一个低频子图
?1??2?cj和水平、垂直与对角线3个方向的高频子图dj,1,dj,2,dj,3。但没有对“水平、垂直与对
角线”做进一步的解释。
以下通过例子分析可分离二维小波变换的计算过程,主要用于理解概念。
第一种方法,很直观,将可分离二维小波变换看成是:先对各行进行小波变换,然后对各列进行小波变换。反之也可以。
例5.1 一个2x2图象的二维Haar小波变换。这里采用非标准的Haar小波滤波器。
小波分析及其应用孙延奎2005
![](/skin/haowen/images/icon_star.png)
![](/skin/haowen/images/icon_star.png)
![](/skin/haowen/images/icon_star.png)
![](/skin/haowen/images/icon_star.png)