双正交小波滤波器构造及其应用于图像边缘检测

  • 格式:pdf
  • 大小:224.81 KB
  • 文档页数:4

双正交小波滤波器构造及其应用于图像边缘检测 王坤鹏, 杨东勇 (浙江工业大学信息工程学院 浙江 杭州 310014) 摘要: 本文在总结边缘检测小波基选取原则的基础上,利用滤波器组技术,提出了具有对称性和正则性的双正交小波滤波器的构造方法,给出了滤波器的构造公式;进行了图像边缘检测实验,结果表明按本文方法构造的滤波器具有很好的图像边缘检测性能。 关键字:边缘检测;滤波器组;正则性;双正交小波滤波器 中图分类号:TP391.41 文献标识码:A

Construction of Biorthogonal Wavelet Filter and its Application to

Image Edge Detection WANG Kun-Peng ,YANG Dong-Yong (Information Engineering College ,Zhejiang University of Technology,Hangzhou 310014,china) Abstract:In this paper, the principle of choosing wavelet base in edge detection is summerized. Then construction of symmetric biorthogonal regularity wavelet filter is presented by biorthogonal filter banks, and filter formula is given. Simulation results for image edge detection demonstrate the effectiveness of the filter constructed by the presented method. Key words: edge detection; filter banks; regularity; biorthogonal wavelet filter

1 引言 小波变换是图像边缘检测的重要工具,小波基的构造和选取是应用小波变换进行边缘检测的重要问题。小波基可以通过小波滤波器方式加以描述,实际中常常通过滤波器组技术构造小波滤波器。本文提出了一种正则、对称的双正交小波滤波器的构造方法,并进行了实验验证。

2边缘检测小波基选取原则 用于图像边缘检测的小波常具有以下特点: (1) 好的时频局部性;(2) 对称性:若小波具有对称性,则没有相位畸变且计算方便[1],

而且最好是关于零轴对称,非零轴对称时小波分解会产生卷积相移[2],导致边缘发生重叠;

(3) 适度的正则性:正则性越高,函数越光滑,在频域中能量越集中,尺度函数的正则性与小波函数的消失矩成对应关系,在某一尺度,检测出的模极值点的数目通常与小波消失矩的数目成线性增长关系[3]。但是正则性与计算效率常是一对矛盾,正则性越高,计算越复杂。

3基于滤波器组技术的双正交小波滤波器构造 定理1:如果滤波器0H,1H满足:()1202222Hzazaazaz−−⎛⎞⎛⎞=−+++−⎜⎟⎜⎟⎜⎟⎜⎟⎝⎠⎝⎠ ;

()1212222221221221221aaaaHzzzzaaaa−−

⎛⎞⎛⎞

=−+−+−⎜⎟⎜⎟

⎜⎟⎜⎟−−−−⎝⎠⎝⎠

(24a≠);则0H,1H可以作为双正交滤波器组的分解低通和高通滤波器,且0H具有偶对称,1H具有奇对称。 证明 令()()01lGzzHz−=+− ,()()10lGzzHz−=−−,则双正交滤波器组的抗混叠条件和延迟条件:()()()()0011HzGzHzGz−+−=0,()()()()0011HzGzHzGz+=

kcz−

满足 ; 双正交滤波器组的传递函数(c=2时)和冲激响应的总和:()()()()011022lkqHzHzHzHzzz−−−−==,()02khk=∑,()10khk=∑满足;

所以0H,1H可以作为双正交滤波器组的分解低通滤波器和分解高通滤波器,且由0H,1

H

的z变换表达式的系数可以看出0H具有偶对称,1H具有奇对称。 定理2:如果低通滤波器0

H满足:

()12

0

22

22Hzazaazaz−−

⎛⎞⎛⎞

=−+++−⎜⎟⎜⎟

⎜⎟⎜⎟

⎝⎠⎝⎠

且2242a<≤,

则0H对应的尺度函数能收敛成连续曲线,即具有正则性。 证明:Daubechies提出[3]:如果0H满足()()()10122pzHzFz−+=,()012Hz==,()10Fz=−≠;()11Fz==,且|()|jFeω在0~2ωπ=范围内的上界值12p−≤,则通过

迭代得到的尺度函数能收敛成连续曲线。此结论与滤波器组是否正交无关。

对0

()Hz进行变形得:()()()()1101222122222zzzHzaa−−⎡⎤++⎢⎥=−−−⎢⎥⎣⎦

令 ()()()12212222zzFzaa−⎡⎤+

⎢⎥=−−−

⎢⎥⎣⎦

则 ()|()|221222cosjFeaaωω

⎡⎤

=−−−

⎣⎦,当2242a<≤时,|()|jFeω在

0~2ωπ=范围内的上界值112−≤,即1≤,而且()012Hz==,()10Fz=−≠,

()11Fz==,所以根据Daubechies提出的正则性条件,命题成立。

根据定理1和定理2,在2242a<≤范围内选取不同的a值(如a取9232,5216,23,328,5212,7216,22等),可以构造出相应的双正交对称正则小波滤波器,分解低通滤波器系数为22,,,22aaaa⎡⎤−−⎢⎥⎣⎦,分解高通滤波器系数为22,,,22221221221221aaaaaaaa⎡⎤−−−⎢⎥

−−−−⎣⎦

当a取23时,分解滤波器系数为222222(,,,),,2,2,633622⎡⎤⎡⎤−−⎢⎥⎢⎥

⎣⎦⎣⎦,对

应的小波记为:symm_regu,使用cascade算法[4]画出的波形如图1(a)。这里还先行给出后面实验对比中用到的非对称正则小波db2和非正则对称小波bior3.1波形,如图1(b)和(c)。

0200400600800-2-1012

0200400600800-1.5-1-0.500.511.52 0200400600800-100-50050100

(a)symm_regu (b)db2 (c)bior3.1 图1 小波波形

4小波边缘检测主要步骤 (1)对输入图像进行小波变换,选用如下算法:

令*(,)jjIRC表示:图像的每一行与一维滤波器jR作卷积,再将结果的每一列与一维

滤波器j

C作卷积。则图像1Sf分解到1(0)Jjj=+≥层的小波变换可通过以下方式实现:

122*(,)jjjjSfSfHH+=;122*(,)xjjjjWfSfHG+=;122*(,)yjjjjWfSfGH+=;

式中jH为H的相邻两项间插入j12−个零值, jG为G的相邻两项间插入j12−个零值。12jSf+为低频分量,12xjWf+为水平高频分量,12yjWf+为垂直高频分量。 (2)逐点计算出小波变换的模值和相角,并对相角做量化处理。 (3)求局部极值点和阈值处理。 (4)二值化并链接边缘像素,去除链长较短的边缘集合,输出图像。

5图像边缘检测实验 笔者进行了仿真实验,采用的图像是标准测试图像Lena,小波分解的尺度为2。 图2(a)为Lena原图,图2(b)为非对称正则小波db2检测的未链接的边缘,图2(c) 为非正则对称小波bior3.1所检测的未链接的边缘;图3(a)为文中构造的正则对称小波symm_regu检测的未链接的边缘,图3(b)为symm_regu小波检测的链接并去除短链后的边缘。 未链接的边缘图对比:图2(b)中 漏检了很多微弱边缘,图2(c) 中不仅漏检严重,而且一片模糊;图3(a)中既检测出了大量的微弱边缘,又较少产生细节,边缘清晰度高,边缘连续性好,这表明按本文方法构造的正则对称双正交小波滤波器具有很好的边缘检测性能。

(a)Lena原图 (b)db2 (c)bior3.1 图2 非对称正则小波db2和非正则对称小波bior3.1检测的未链接的边缘

(a)symm_regu未链接边缘 (b)symm_regu链接后边缘 图3 文中构造的正则对称小波symm_regu检测的未链接边缘及链接后的边缘

6结束语 本文总结了边缘检测中的小波基选取原则,提出了具有对称性和正则性的双正交小波滤波器构造公式。实验结果表明,按本文方法构造的小波滤波器具有很好的图像边缘检测性能。

参考文献 [1] Sun Mingui, Li Ching-Chung, Sclabassi Robert . Symmetrical Wavelet Transforms for Edge Localization[J]. Optical Engineering, 1994,33(7):2272~2281. [2] Sun Mingui, Sclabassi Robert . Symmetric Wavelet Edge Detector of the Minimum Length[c]. In: International Conference on Image Processing, 1995,2177~2180. [3] Stephane Mallat, Wen Liang Huang. Singularity detection and processing with wavelets[J]. IEEE Trans on Information Theory,1992, 38(2): 617~643. [4] Strang Gilbert, Truong Nguyen. Wavelets and Filter Banks[M]. Wellesley-Cambridge Press,1996. 作者简介: 王坤鹏(1976-),男,汉族,在读硕士研究生,研究方向:图象处理、计算机视觉。 杨东勇(1961-),男,博士,教授,博士生导师,研究方向:计算机智能系统,计算机视觉。 基金项目:浙江省自然科学基金(编号:601078)。 作者:王坤鹏;联系地址:浙江工业大学研究生23信箱,邮编:310014; E-mail:kunpeng8069@sina.com;