小波多尺度边缘检测

1、第二章小波多尺度边缘检测§1多尺度边缘检测的基本原理大多数多尺边缘检测器都是在不同的尺度平滑信号，然后由其一阶或二阶导数检测锐变点，所谓尺度实际上是计算信号变化的范围。平滑函数0(x):其积分等于1,且当xT±g时速降至零，例如高斯函数，平滑函数0(X)的一阶、二阶导数分别为(21)屮a(X)=学)，屮b(X)=d笄dxdx2显然，|7a()二j®0),|7b()二(j®)20)，由于0八(0)二1故6a(0)和9b(0)均为零,从而Va(X)和甲b(X)都是满足允许条件的小波。在本章以后的讨论中，2(X)表示将g(X)按尺度S伸缩的同时保持面积不变，即

2、sg(X)A-g(-)(22)S=SS将小波变换定义为信号f(X)与屮a(X)和屮b(X)的卷积积分，即1卡xT(23)(24)(25)(26)Waf(X)二f*屮a(X)二Jf(T)Va(ssssg1gXTWbf(X)=f*vb(X)=Jf(T)Vb(-)-TSSSSg由此可以导出如下重要结论Waf(X)二f*(S伴)二sf(f*0)sdXdXss-wbf(X)=f*(S2-s)=S2£(f*0)sdX2dX2s由上列两式可以看到，边缘检测可以通过小波变换来实现，边缘实际上是一阶导数的极值点，即二阶导数的过零点，也就是说，我们可以通过寻找Waf(X)的极值点或Wbf(X)的SS过零

3、点来确定边缘的位置，但是，下面我们将会看到，通过分析|waf(X)|的极大值和尺S度S的关系，进而确定边缘的性质，故寻找一阶导数的极值点较寻找二阶导数过零点的方法会获得更多关于边缘的信息。为了定量地描述一个函数的奇异性，我们首先引入Lipschitz指数的定义。定义：令o<a<1，当且仅当对任意的(x,x)e(a,b)，总有在一个常数k，使01f(xo)-f(X<乜-£(27)我们称f(x)在区间(a，b)是一致Lipschitzd,a的上界值称为Lipschitz指数。不难理解，如f(x)在xo点可微，则其Lipschitz指数至少为1。实际上，Lipschitz

4、指数越大，函数越光滑。如f(x)在X点不连续但在X的邻域有界，则其Lipschitz指数oo为0。我们也可以将Lipschitz指数推广到为负数的情况：如f(x)的原函数在%点的Lipschitz指数为a，则它在该点的Lipschitz指数为a-1，例如§(xo)的原函数在%为一单位阶跃，其Lipschitz指数为0,故§(xo)的Lipschitz指数为一1。下面，我们要将Lipschitz指数和小波变换联系起来。定理：令0<a<1，当且仅当对任意的xe(a,b)，总存在一个常数k>0，使|wf(x)|<ksa(28)s则函数f(x)在区间(a,b

5、)是一致Lipschitza(28)式可以写为logf(x)|<log+alog(29)s前面我们已经指出，|wf(x)的极大值点指明了边缘的位置，那么边缘的性质如何呢？s由上式可以看到，小波变换模的极大值是随尺度s而变化的,如按对数取值,logwf(x)|s与logs具有线性关系，它们之间的比例系数(直线的斜率)即为Lipschitz指数。§2二进小波变换我们已经学习过小波级数,那时是将小波变换中的伸缩参数和平移参数都离散化,而二进小波变换只是将(23)和(24)式中的尺度参数离散化为2的整次幕，即s二2j,jeZ,类似于(22)式，屮(x)表示将屮(x)作二进伸缩的同时保持

6、其面积不变，即j屮(x)屮(亍)，jeZ(210)j2j2j类似于(23)和(24)式，尺度为2j时的小波变换为1xTwf(x)f*屮(x)If(TN(-)dT(211)jj2jR2j而二进小波变换是所有尺度时小波变换的总体，即WfJf(x)，jeZ(2"上式中w为二进小波变换算子。下面我们介绍如何从信号的二进小波变换重构信号，较深入的讨论将在框架理论中，为了能从二进小波变换重构信号，或者说，信号的二进小波变换并未丢失信号的信息，则屮(X)必须满足如下充分必要条件：存在两个正数，使得下列不等式成立(213)(214)A<艺M(2j)|2<B,V®wRj=g由(2

7、11)式可以得到wf()=/()v|7()=/()屮(2j)jj也就是说，wf(x)相当于信号通过一个带通滤波器，所以二进小波变换相当于信号通j过中心频率和带宽不同的带通滤波器组，而(213)将确保带通滤波器组覆盖整个频率轴,从而不丢失f(x)的信息，利用Parseval定理，从(213)和(214)可导出(213)的等价表达式|2<Bf(215)Af112<艺”f(x)j=s上述不等式不仅确保二进小波变换是完备的(可以从二进小波变换重构信号)，而且是稳定的(重构公式将有很好的收敛性)，B/A越接近1,稳定性越好。在小波级数中，由于屮(x)经二进伸缩和整数平移后构成正交基，所以在进

8、行分解和重构时都是使用同样的小波函数族，以后我们学习框架理论时将会看到，完全可以放松正交性的要求，但这时分解和重构时将使用不同的小波函数族，设二进小波变换的重构小波为X(x)，其傅里叶变换必须满足如下条件(216)jwz将按下式重构信号£vp(2j®)X(2j®)=1j=-8f1X重构小波X()经二进伸缩形成的函数族(x)=X(亍),Ij2j2jf(x)=wf*X(x)(217)jjjs上列重构公式很容易证明，将上式取傅里叶变换，并将(214)代入，同时引用重构小波条件式(216)，便可证明(217)的正确性，需要指出的是，满足(216)式的重构小波将会有无穷多个。和在MRA的Mallat算法一样，在数字应用中，输入