基于B样条小波的图像边缘检测.

基于B样条小波的图像边缘检测

周何，黄山，盛贤

（四川大学电气信息学院自动化系，成都市610065；）

摘要：研究图像边缘优化检测问题。针对图像边缘信息被噪声污染影响定位精度，经典的边缘检测方法Canny算法中的高斯平滑函数边缘定位精确度较低，导致图像缓变边缘信息丢失和假边缘的现象。在Canny最优边缘检测准则下，引入了渐进最优的B样条小波函数，采用小波变换应用于图像边缘检测中的基于模极大值的方法，并结合基于Kmeans聚类的自适应双阈值方法进行图像边缘检测。实验结果表明，改进的算法改善了噪声干扰情况下图像边缘提取效果，有效提高了边缘检测的准确性，得到较高的边缘检测图像质量。

关键词：边缘检测；小波变换；定位精度；

中图法分类号: TP391.4文献标识码: A

Image edge detection based on B-spline wavelet

ZHOU He，HUANG Shan，SHENG Xian

(School of Electrical Engineering and Information, Sichuan University, Chengdu 610065, China;) Abstract:In order to solve the low positioning accuracy of image edge detection by noise, make a research on optimization of image edge detection. The Gaussian smoothing function of Canny edge detection method, the classical algorithm, causes the missing of slowly varying edge and the producing of feigned edge and the edge detection is not accurate enough. So in the Canny criteria of optimum edge detection, the introduction of the asymptotically optimal B-spline wavelet function was put forward. The method of modulus maxima of wavelet transform and Kmeans clustering method determining its duel valves automatically was used in the edge detection experiments.The experiments proved that the new algorithm was in a higher accuracy, and improved the quality of the edge detection image.

Keywords : edge detection; wavelet transform; positioning accuracy;

1 引言

传统的边缘检测Canny算法是将图像与高斯滤波器相卷积以获得平滑降噪的效果，其基本思想是在图像中找出具有局部最大梯度幅值的像素点，对边缘提取的大部分工作集中在寻找能够用于实际图像的梯度数学逼近。这种算法会造成原图像的过度光滑，缓变边缘丢失，定位精度较低，且计算量大、复杂、耗时[1]。

小波分析具有多尺度分析的特点，能较好的综合噪声抑制和边缘保持这两个特性。任意一个信号可表示成经伸缩和平移的n次B样条的加权和，即可完全由B样条系数来刻画。该系数中的分辨阶数越小对信号的平滑程度越小，边缘定位越精确，在对不同尺度下的逼近函数取一阶导数或者二阶导数时就获得了多尺度边缘提取。

本文充分利用边缘信息的多尺度特性和B 样条函数是同次样条函数空间中具有最小支撑的基底的这一特点，选取正交三次中心B样条

作为边缘提取时的平滑函数，再采用模极大值和Kmeans聚类的自适应双阈值的方法，提取出最终的边缘图像。此算法的原理与实现简单，且有较好的抗噪性能，并拥有比以Gauss函数为平滑函数的Canny算法更加出色的定位精度，提取出了更加精细的边缘，去除了虚假边缘。

2 B样条小波

在对Canny边缘检测算法的应用和研究中发现，Canny算法用Gauss函数作为滤波器，会使原图像过度光滑，缓变边缘丢失。由于Canny 算子不能直接进行Z变换，即找不到递推公式，从而只有用它进行卷积运算。但对于一个大的图像，计算时间很长。为此，在Canny最优边缘检测准则下，引入了渐进最优的B样条小波函数。

2.1 Canny边缘提取准则

John Canny于1986年在IEEE 上发表了自己的文章《A Computational Approach to Edge

1

2Detection 》，在其文章中指出了三个准则：[2]

(1) 好的检测结果（Good detection ）一个好的检测结果应该尽可能小地漏检真实存在的边缘点和误检非边缘点。

设用于边缘测定滤波器为()f x ，它的有限相应边界为[,]W W -。边缘为()G x ，边缘发生在0x =，信号中的信噪比是高斯白噪声()n x ，

其方差为2

0n ，第一个准则数学表达式定义

|()()|

W

G x f x dx SNR +-=

（1）

(2) 定位精确（Good localization ）：标记为

边缘的点应当尽可能地接近真实边缘的中心。检测精确定义为：

''

'2

0()()()W

W

W

W

G x f x dx

L n f x dx

+-+--=

⎰

⎰

（2）

(3) 对同一边缘响应次数较少（Only one response to a single edge ）理想情况下，用滤波器对噪声响应的峰值距离来近似滤波器对一个边缘的响应宽度。而两个邻近的极大值平均距离为滤波器输出导数的零交叉两倍。零交叉点平均距离为

2

'2''2

()*()zc f x dx x f x dx π+∞-∞+∞-∞⎛⎫ ⎪= ⎪ ⎪⎝⎭⎰⎰ （3）

有了这三个准则的数学表达式，寻找最优滤

波器就转化为泛函约束优化的问题。 2.2 B 样条小波定义

小波基的种类很多，如何根据实际问题选择最佳的小波基是首要关注的问题。用于边缘检测的小波基函数应为一个紧支撑的奇函数小波。根据Canny 准则的前两个准则，对阶跃边缘的最优检测函数为阶段阶跃或者差分盒函数。定义一阶B 样条函数[3]为:

01,[0,1]

()0,other x x β∈⎧=⎨

⎩

（4）

等距单重结点条件下，n 次中心样条函数

()n x β用卷积定义为:

100

001

()()()()()...()

n n n x x x x x x ββββββ-+=*=***

10111(1)!2n

n k k n n x k k n +=+⎡⎤+⎡⎤=--+⎢⎥⎢⎥⎣

⎦⎣⎦∑ (5) 其中()max{0,}x x +=。易知，()n x β是非负的，其支撑集为11,22n n ++⎡⎤

-

⎢⎥⎣⎦

。 考虑一系列相互嵌套的多项式样条函数空间(){,}

n i S i Z ∈(n 暂取为奇数，为多项式的阶数)，使得()(1){}n n i i S S +⊃。对于i Z ∈，()

n

i S 为2()L R 的子空间，属于1n C -类（即具有n-1次连续导数）。

在每一区间2,(1)2j j k k ⎡⎤+⎣⎦上()n i S 等价于n 阶

多项式。

()()

2{()()(2)}

j

n n n j i i i

k S f x C k x k β

+∞

=-∞

==

-∑

2,()i x R C k

l ∈∈ （6） 其中，21()22

j n

n j j x

ββ=

。()n x β的Fourier 变换为

1

1

0sin(/2)()()/2n n n ωβωβωω++⎡⎤

⎡⎤==⎢

⎥⎢⎥

⎣

⎦

⎣⎦

(7)

{(),}x k k Z β-∈构成空间(0)

n

S 的Riesz 基，而(){,}n i S i Z ∈构成了2

()L R 空间的多分辨率分析，即