基于小波变换的边缘检测技术(完整

1、第一章 图像边缘的定义引 言 在实际的图像处理问题中，图像的边缘作为图像的一种基本特征，被经常用于到较高层次的特征描述，图像识别。图像分割，图像增强以及图像压缩等的图像处理和分析中，从而可以对图像进行进一步的分析和理解。由于信号的奇异点或突变点往往表现为相邻像素点处的灰度值发生了剧烈的变化，我们可以通过相邻像素灰度分布的梯度来反映这种变化。根据这一特点，人们提出了多种边缘检测算子：Roberts算子 Prewitt算子 Laplace算子等。经典的边缘检测方法是构造出像素灰度级阶跃变化敏感的微分算子。这些算子毫无例外地对噪声较为敏感。由于原始图像往往含有噪声、而边缘和噪声在空间域表现为灰度有大

2、的起落，在频域则反映为同是主频分量，这就给真正的边缘检测到来困难。于是发展了多尺度分析的边缘检测方法。小波分析与多尺度分析有着密切的联系，而且在小波变换这一统一理论框架下，可以更深刻地研究多尺度分析的边缘检测方法，Mallat S提出了一小波变换多尺度分析为基础的局部极大模方法进行边缘检测。小波变换有良好的时频局部转化及多尺度分析能力，因此比其他的边缘检测方法更实用和准确。小波边缘检测算子的基本思想是取小波函数作为平滑函数的一阶导数或二阶导数。利用信号的小波变换的模值在信号突变点处取局部极大值或过零点的性质来提取信号的边缘点。常用的小波算子有Marr 算子Canny算子 和Mallat算子等。

3、1.1信号边缘特征人类的视觉研究表明，信号知觉不是信号各部分简单的相加，而是各部分有机组成的。人类的信号识别（这里讨论二维信号即图像）具有以下几个特点：边缘与纹理背景的对比鲜明时,图像知觉比较稳定；图像在空间上比较接近的部分容易形成一个整体；在一个按一定顺序组成的图像中，如果有新的成份加入，则这些新的成份容易被看作是原来图像的继续；在视觉的初级阶段，视觉系统首先会把图像边缘与纹理背景分离出来，然后才能知觉到图像的细节，辨认出图像的轮廓，也就是说，首先识别的是图像的大轮廓；知觉的过程中并不只是被动地接受外界刺激，同时也主动地认识外界事物，复杂图像的识别需要人的先验知识作指导；图像的空间位置、方向

4、角度影响知觉的效果。从以上这几点，可以总结出待识别的图像边缘点应具有下列特征即要素：具有较强的灰度突变，也就是与背景的对比度鲜明；边缘点之间可以形成有意义的线形关系，即相邻边缘点之间存在一种有序性；具有方向特征；在图像中的空间相对位置；边缘的类型，即边缘是脉冲型、阶跃型、斜坡型、屋脊型中哪一种。1.2图像边缘的定义 边缘检测是图像处理中的重要内容。而边缘是图像中最基本的特征，也是指周围像素灰度有变化的那些像素的集合。主要表现为图像局部特征的不连续性，也就是通常说的信号发生奇异变化的地方。奇异信号沿边缘走向的灰度变化剧烈，通常分为阶跃边缘和屋顶边缘两种类型。阶跃边缘在阶跃的两边的灰度值有明显的变

5、化；屋顶边缘则位于灰度增加与减少的交界处。我们可以利用灰度的导数来刻画边缘点的变化，分别求阶跃边缘和屋顶边缘的一阶，二阶导数。如图可见，对于边缘点A，阶跃边缘的一阶导数在A点到最大值，二阶导数在A点过零点；屋顶边缘的一阶导数在A点过零点，二阶导数在A点有最大值。第二章 传统的边缘检测算子 2.1传统的边缘检测算子 边缘检测的实质是采用某种算法来提取出图像中对像与背景间的交界线。我们将边缘定义为图像中灰度发生急剧变化的区域边界。而灰度变化的情况可以用图像灰度分布的梯度来反映，所以我们可以用局部图像微分技术来获得边缘检测算子。以下对比较经典的边缘检测算子进行了理论分析，并做出了比较和评价。我们记为

6、图像的梯度，中包含局部灰度的变化信息。记： 为梯度的幅度，可以用做边缘检测算子。常用的边缘检测方法有：差分边缘检测，Roberts边缘检测算子，Sobel边缘检测算子，Prewitt边缘检测算子,Robinson边缘检测算子,Lapalce边缘检测算子等等。2.2 差分边缘检测方法利用像素灰度的一阶导数算子在灰度迅速变化处得到高值来进行奇异点的检测。它在某一点的值就代表该点的“边缘强度”，可以通过对这些值设置阈值来进一步得到边缘图像。但用差分边缘检测必须使差分的方向与边缘方向垂直，这就需要对图像的不同方向都进行差分运算，增加了实际运算的繁琐性。一般为垂直边缘、水平边缘、对角线边缘检测：图2-1

7、 差分算法检测边缘的方向模板2.3 Roberts边缘检测算子Roberts边缘检测算子根据任意一对互相垂直方向上的差分可以用来计算梯度的原理，采用对角线方向相邻两像素之差，即：他们的卷积算子为： 有了，之后，很容易计算出Roberts的梯度幅值，适当的取门限TH，作如下判断：，为阶跃边缘点。为边缘图像。Roberts算子采用对角线方向相邻两像素之差近似梯度幅值边缘检测。检测水平和垂直边缘的效果好于斜向边缘，定位精度高，对噪声敏感。 图2-2：用Roberts算子进行边缘检测的Lena图与原图像2.4 Sobel边缘检测算子对数字图像的每一个像素，考察它上，下，左，右邻点灰度的加权差，与之接近

8、的邻点的权大。据此，定义Sobel算子如下：卷积算子为：图2-3： Sobel边缘检测算子方向模板适当的取门限TH，作如下判断：为阶跃边缘点，为边缘图像。Sobel算子很容易在空间上实现，Sobel边缘检测器不但产生较好的边缘检测效果，而且受噪声的影响也比较小。当使用大的领域时，抗噪声特性会更好，但这样做会增加计算量，并得出的边缘也比较粗。Sobel算子利用像素点上下，左右邻点的灰度加权算法，根据在边缘点出达到极值这一现象进行边缘的检测。Sobel算子对噪声具有平滑作用，提供较为精确的边缘方向信息，但它同时会检测出许多的伪边缘，边缘定位精度不高。当对精度要求不是很高时，是一种较为常用的方法。

9、图2-4：用Sobel算子进行边缘检测的Lena图与原图像2.5 Prewitt边缘检测算子Prewitt算子是一种边缘样板算子。这些算子样板由理想的边缘子图像构成。依次用边缘样板去检测图像，与被检测区域为相似的样板给出最大值。用这个最大值作为算子的输出值，这样就可以将边缘像素检测出来。定义Prewitt边缘算子模板如下：图2-5： Prewitt边缘检测算子模板8个算子样板对应的边缘方向如下图所示： 图2-6： 样板方向 适当取门限TH，作如下判断：为阶跃边缘点。为边缘图像。 图2-7： 用Prewitt算子进行边缘检测的Lena图与原图2.6 Robinson 边缘检测算子Robinson

10、 边缘检测算子也是一种边缘样板算子，其算法和Prewitt边缘检测算子相似，只是8个边缘样板不同。如下所示： 图2-8：Robinson 边缘检测算子模板2.7 Laplace边缘检测算子Laplace算子是二阶微分算子，是一个标量，属于各向同性的运算，对灰度突变敏感。在数字图像中，可以用差分来近似微分运算，的Laplace算子为Laplace算子的二种估算模板： 图2-9：Laplace的两种估算模板对阶跃边缘，二阶导数在边缘点出现零交叉，即边缘点两边二阶导函数取异号。Laplace算子就是据此对的每个像素取它关于方向和方向的二阶差分之和，这是一个与边缘方向无关的边缘检测算子。而对屋顶状边缘

11、，在边缘点的二阶导数取极小值，这时对的每个像素取它关于方向和方向的二阶差分之和的相反数。 Laplace算子有两个缺点：其一是边缘的方向信息丢失，其二是Laplace算子为二阶差分，双倍加强了图像中的噪声影响：优点是各向同性，即具有旋转不变性。因为在微分学中有：一个只包含偶次阶导数和取偶次幂的奇次阶导数的线形组合算子，一定是各向同性的。Laplace算子实际二阶微分算子，利用边缘点处二阶导函数出现零交叉原理检测边缘。不具有方向性，对灰度突变敏感，定位精度高，不但检测出了绝大部分的边缘，同时基本上没有出现伪边缘。但他的检测也存在一些缺点，如丢失一些边缘，有一些边缘不够连续，对噪声敏感且不能获得边

12、缘方向的功能信息。 图2-10：用Laplace算子进行边缘检测的Lena图与原图2.8 检测结果与结论通过对以上介绍的几种边缘检测算子的算法公式和检测的结果可以看出，Roberts算子简单直观，但边缘检测图里存在有伪边缘；Sobel算子、Prewitt算子和Robinson的检测结果图能检测出更多的边缘，但也存在有伪边缘且检测出来的边缘线比较粗，并放大了噪声；Lapalce算子和改进的Laplace算子利用二阶差分运算来进行检测，但不可以检测出较多的边缘，而且还在很大程度上消除了伪边缘的存在，定位精度高。但受噪声的影响比较大。 第三章 小波变换在图像边缘检测中的应用3.1 小波思想的引入虽然

13、边缘提取已有梯度算子、Laplace算子、Sobel算子等方法，但这些算法都没有自动变焦的思想。而事实上，由于物理和光照等原因，每幅图像中的边缘通常产生在不同的尺度范围内，形成不同类型的边缘，这些信息是未知的。另外图像中还存在有噪声，因此，根据图像特性自适应地正确检测出图像的边缘是非常困难的。可以肯定，用单一尺度的边缘算子不可能检测出所有的边缘，同时，为避免在滤除噪声是影响边缘检测的正确性，用多尺度的方法检测边缘越来越引起人们的重视。由于小波变换具有良好的时频局部化特性及多尺度分析能力，在不同尺度上具有“变焦”的功能，适合于检测突变信号。是检测突变信号强有力的工具，得到了广泛的应用。3.1.1

14、小波变换在图像边缘检测中的优势 用小波变换对信号做多分辨率分析非常适合提取信号的局部特征。这是因为小波变换的尺度因子和平移因子构成了一个滑动的时间-频率窗，小尺度下的变换系数对应信号的高频分量，大尺度下的变换系数对应信号的低频分量。于是信号被分解成各个频率下的分量，这样就可以检测对应不同频率的信号局部特征。而图像中的突变信息和噪声都属于高频信号，可以利用小波变换后的高频分量进行去噪和得到边缘图像。由于各种原因，图像常常受到随机噪声的干扰。经典的边缘检测方法由于引入了各种形式的微分运算，从而必然引起对噪声的极度敏感，在其上执行边缘检测的结果常常是把噪声当作边缘点检测出来，而真正的边缘也由于受到噪

15、声干扰而没有被检测出来。因而对于有噪声图像来说，一种好的边缘检测方法应具有良好的各种噪声抑制能力，同时又有完备的边缘检测保持特性。小波变换可以提供一种很好的去噪方法。当取小波函数为平滑函数的一阶导数时，信号的小波变换的模的信号突变点出取局部极大值，边缘与噪声的区别在于，随着尺度的增加，噪声引起的小波变换的模的极大值迅速减小；而边缘的滤波模值不随尺度的变化，故小波变换可以在低信噪比的信号中检测噪声和边缘。通过计算在尺度和尺度上每一个在位置上最接近，且具有相同符号的系数最大值，我们可以找出不同尺度下小波幅度的变化，消除那些系数极值孤独随尺度减小而小波系数Mf在平均值上增加的序列。这些极值对应图像噪

16、声奇异点，从而就得到了图像边缘的真正奇异点。同时，边缘点具有较强的方向性。即在边缘点上，其滤波值表现为较强的模值，而在垂直边缘方向上，滤波模值较小；而噪声点在各个方向的滤波模值相似。因而，在对信号进行不同方向的滤波时，边缘点滤波模值会有较强与较弱的差异，而噪声点滤波模值不会产生过大的差异，因此可以通过构造任意方向的小波滤波器来检测出奇异信号和噪声点。对于很大一类图像来说，急剧变化点所对应的边缘在图像平面内是一些规则的曲线，沿着这些曲线，图像在一个方向是奇异的，而和该方向相垂直的方向却是平滑的。由于实际得到的值是一些离散点，为了检测出边缘，可以利用沿曲线方向梯度矢量的模变化平缓特点以及不同尺度下

17、梯度矢量幅角的信息，将位置及幅角十分接近的模极大值点连接起来形成模极大值链，这些链构成了图像的边缘。3.1.2小波变换定义及特点小波(wavelet)，即小区域的波，是一种特殊的长度有限、平均值为0的波形。它有两个特点：一是“小”，即在时域都具有紧支集或近似紧支集；尺度二是正负交替的“波动性”，也即直流分量为零。我们可以用小波和构成傅立叶分析基础的正弦波做一个比较，傅立叶分析所用的正弦波杂时间上没有限制，从负无穷到正无穷，但小波倾向于不规则也不对称。傅立叶分析是将信号分解成一系列不同频率的正弦波的叠加，同样小波分析是将信号分解成一系列小波函数的叠加，而这些小波函数都是有一个母小波函数经过平移与

18、尺度伸缩得来的。我们知道，用不同规则的小波函数来逼近尖锐变化的信号显然要比光滑的正弦曲线要好，同样，信号局部的特性用小波函数来逼近显然要比光滑的正弦函数来逼近要好。小波变换的定义是把某一被称为基本小波（也叫母小波）的函数做位移后，再在不同尺度下与待分析的信号做内积：， 等效的频域表示是：可以这样理解上面表达式的意义：打个比喻我们用镜头观察目标（即待分析信号），代表镜头所起的作用（例如滤波或卷积）。相当于使镜头对于目标平行移动，的作用相当于镜头向目标推进或远离。由此可见小波变换有以下特点： 有多分辨率，也叫多尺度的特点，可以有粗及细地逐步观察信号。 可以看成用基本频率特性为的带通滤波器在不同尺度

19、下对信号滤波。由于傅立叶变换的尺度特性可知这组滤波器具有品质因数恒定，即相对带宽（带宽与中心频率之比）恒定的特点。注意，越大相当频率越低。 适当的选择基小波，使在频域上也比较集中，就可以使在在时、频域都具有表征信号局部特征的能力，因此有利于检测信号的瞬态或奇异点。正是由于上述特性，有人把小波变换称为分析信号的数学显微镜。如上所述，小波分析的一个重要优点就是能够分析信号的局部特征。利用小波分析可以检测出许多其他分析方法忽略的信号特性，例如，信号的趋势，信号的高阶不连续点、自相似特性。总之，小波变换作为一种数学理论和方法在科学技术和工程界引起了越来越多的关注和重视。3.2连续小波变换3.2.1连续

20、小波的定义及其变换定义3 . 1 设,其傅立叶变换为满足允许条件（Admissible Condition ，完全重构条件或恒等分辨条件） 我们称为一个基本小波或母小波（Mother Wavelet）。将母函数按此方式经平移和伸缩后得到函数族： 称其为一个小波序列或连续小波。其中a 为伸缩因子，b 为平移因子。我们首先来看看连续小波在相空间中的局部化格式。假定是双窗函数，记由于正频和负频讨论类似，只对正频进行讨论。注意到 因此的正频窗口中心为其时窗宽，频窗宽。由此可以看出，b 仅仅影响窗口在相平面时间轴上的位置，而a 不仅仅影响窗口在频率轴上的位置，也影响窗口的形状。当越小时，时宽越小，而频宽

21、越大，且窗口中心向增大方向移动，这表明高频时连续小波的时间分辨率较高，频率分辨率较低；反之，当越大时，时宽越大，而频宽越小，窗口中心向减小方向移动，这表明低频时连续小波的时间分辨率较差，频率分辨率较高。这一特性正好符合低频信号变化缓慢而高频信号变化迅速的特点，决定了小波变换在突变信号处理上的特殊地位及功能。在这个意义上，小波变换被誉为数学显微镜。对于任意函数的连续小波变换为 其重构公式（逆变换）为有时为了方便起见，小波变换也常用如下定义 其中 式和式定义形式有所不同，主要在于1 ）伸缩系数不同。2 ）由卷积代替了相关。但事实上，它们之间是可以相互转换的。 由基本小波生成的小波在小波变换中对被分

22、析的信号起着信号观测窗的作用，所以还应满足一般函数的约束条件 故是一个连续函数。这意味着，为了满足允许条件，在原点必须等于零，即 小波只有在实轴上取值有正有负才能保证上式的积分为零，所以应具有振荡性，也就是说是一个“波”。同时，的定义域是具有紧支撑的，即在一个很小的区间之外，函数值为零，应是一个迅速衰减的“小”波。小波因此得名。如果使下式 成立，则说为阶消失距小波（Vanishing moments）这时在0 处是次可微的，即。而时就是容许条件，随着的增加，小波的振荡性就越来越强烈。 连续小波变换具有以下重要性质：1、线形性：一个信号的小波变换等于各个分量的小波变换之和。2、平移不变性：若的小

23、波变换为，则的小波变换为3、伸缩共变性：若的小波变换为，则的小波变换为 4、自相似性：对应不同尺度参数和不同平移参数的连续小波变换之间是自相似的。5、冗余性：连续小波变换中存在信息表述的冗余度。3.2.2 离散小波变换在实际运用中，特别是在计算机实现上，往往需要把上面提到的连续小波及其变换离散化。需要注意的是，这一离散化是针对连续的尺度参数和连续的平移参数的，而不是针对时间变量的。这一点和我们在离散傅里叶变换中熟悉的时间离散化是不同的。在连续小波中，考虑函数：这里，。为方便起见，在离散化中，总限制只取正值，即且，是容许的，这样相容性条件就变成通常，把连续小波变换中尺度参数和平移参数的离散公式分

24、别取作，这里，扩展步长是固定值。为了方便起见，总假定。所以对应的离散小波函数即可写作而离散化的小波变换系数则可表示为其重构公式为 是一个与信号无关的常数。显然，和决定了信号重构的精度，即网格点应尽可能地密（和尽可能小），才能保证一定的重构精度。如果网格点越稀疏，使用的小波函数和离散小波系数就越少，信号重构的精确度也就会降低。3.3常用小波函数介绍3.3.1 Haar小波Haar小波是在小波分析中最早用到的一个具有紧支撑的正交小波函数，同时也是最简单的一个函数。Haar函数的定义为 3.3.2 Daubechies小波系 Daubechies函数是由世界著名的小波分析学者Inrid Daubec

25、hies构造的小波函数，提供了比Haar小波更有效的分析和综合。除了dbl(即Haar小波)外，其它小波没有明确的表达式。但转换函数h的平方模是很明确的。 假设 其中，为二项式的系数，则有 其中： 小波函数和尺度函数的有效支撑长度为，小波函数的消失矩阶数为 N,且具有正交性。常用的Daubechies小波函数有db4和db8两种小波。3.3.3 Morlet小波 Morlet函数定义为 它的尺度函数不存在，且不具有正交性。3.3.4 Mexican Hat小波Mexican Hat函数为 它是Gauss函数的二阶导数，像墨西哥帽的截面，也称为墨西哥函数。它在时间域和频率域都有很好的局部化，并且

26、满足。由于它的尺度函数不存在，所以不具有正交性。3.3.5 Meyer小波Meyer小波的小波函数和尺度函数都是在频率域中进行定义的，是具有紧支集的正交小波。 其中， 为构造Meyer小波的辅助函数，具有 图3-1：几种常见的小波函数图形第四章 小波分析中的多尺度思想和边缘检测算子4.1 多尺度边缘检测思想图像在不同尺度上的小波变换都提供了一定的边缘信息。当小尺度时，图像的边缘细节信息较为丰富，边缘定位精度较高，但易受到噪声的干扰；大尺度时，图像的边缘稳定，抗噪性好，但定位精度差。在实际应用中经常存在着去除噪声和准确定位之间的矛盾。多尺度边缘检测的基本思想就是沿梯度方向，分别用几个不同尺度的边

27、缘检测算子在相应点上检测模极大值的变换情况，并通过对阈值的选取，再在不同尺度上进行综合得到最终边缘图像，可以较好的解决噪音和定位精度之间的矛盾。在信号的时频局部化分析中，虽然Fourier变换能较好地刻画信号的频率特性，但不能提供信号在时域上的信息，为此，Gabor首先提出加窗Fourier分析，但加窗不能敏感地反映信号的突变。小波理论与变换方法能弥补上述不足，主要表现在高频处的时间分辨率高，低频处的频率分辨率高，即具有变焦特性，小波变换是一种非平稳信号分析方法，它是通过一个基本小波函数的不同尺度的平移和伸缩构成一族小波函数系表示来逼近一信号。在小波多尺度分析中，引入了尺度函数，其伸缩与平移系

28、数构成矢量空间的正交基，构成空间在分辨率上组成逼近空间。尺度函数具有低通滤波作用，并满足双尺度方程 而与尺度函数相对应的小波函数的平移与伸缩构成矢量空间的正交补空间，小波函数具有高通滤波的作用，并满足方程 对任意一平方可积函数或信号，在尺度下的离散小波变换为 平滑信号为 小波变换与平滑分量满足递推公式其中； 信号的多分辨率分解与重构可由上式完成，推广至二维图像空间，在图像被分解成1个轮廓图像和3个细节图像。其中 代表了图像的低频成分，低频图像还可以进一步分解成4个子带，设分解层数为，则总的子带数为代表了图像的垂直高频成分（水平边界）；代表了图像的水平高频成分（垂直边界）；代表了图像的对角线上的

29、高频成分。Mallat S提出了离散小波变换的一种快速算法，即用低通滤波器H和高通滤波器G对原始数据进行逐步分层分解。其结果是产生一个低频图像，1个水平方向的细节图像，1个垂直方向的细节图像和1个对角方向的细节图像。另外小波变换的多尺度分析(或多分辨率分析) 是建立在函数概念上的，其创建者Mallat S 正是在研究数字图像处理的问题时建立此理论的.随着尺度由大到小变化,在各尺度上可以由粗及精地观察图像的目标.大尺度时,观察到的是图像的基本特征;在小尺度的空间里,则可以看到目标的细节.把二维图像信号所占据的总频带定义为空间,用理想的低通滤波器和高通滤波器在行、列方向将它分别分解成低频部分和高频

30、部分,每一方向的两部分分别反映出该图像信号在剖分方向上的概貌和细节;对于经第二级( a = 21) 分解后又被剖分成低频、水平方向的高频、垂直方向的高频、以及对角线方向的高频在这种空间剖分过程中, 反映的是图像信号在空间中沿方向概貌的低频子空间，反映的是图像信号在空间中沿方向细节的高频子空间。从多分辨率分析可以看出,空间的每次剖分包含两部分:一部分是图像信号通过低通滤波后得到的低频概貌;另一部分是通过带通滤波(小波变换) 得到的图像高频细节;对于低频概貌,重复以上过程,最终把图像信号分解成多个等级的高频细节与最后一次低通滤波后的低频概貌之和。在剖分过程中,这些子空间具有以下的特性：2 单调性:

31、 2) 逼近性: 3) 伸缩规则性: 4) 平行不变性: 满足的上述性质称为多尺度分析即:任意函数应用多尺度分析将其分解为细节部分或是某一方向上的细节部分和的基本特征部分然后将进一步分解,可得到任意尺度下的的基本特征部分以及细节部分之和。按照多尺度分析理论,描述函数分解过程如图所示图4-14.2 常用的几种边缘检测算子近年来受到广泛关注的边缘检测方法是Canny算子以及Mallat边缘检测算子。下面就简单地作以介绍：1 Canny边缘检测算子Canny算子使用一阶导数的极大值表示边缘。其基本思想是先将图像使用Gauss函数进行平滑，再由一阶微分的极大值确定边缘点。二阶导数的零交叉点不仅对应着一

32、阶导数的极大值也对应着一阶导数的极小值，也就是说，灰度变化剧烈的点与灰度变化缓慢的点都对应着二阶导数零交叉点不仅对应着一阶导数的极大值也对应着一阶导数的极小值。 在用高斯函数对图像进行滤波得到，其中为相应的尺度因子。计算其梯度矢量的模和方向：图像边缘点即为在方向上使取得局部极大值的点。2 Mallat小波边缘检测算子 基于Canny边缘检测算子的主要思想，Mallat在1992年提出了小波变换极大值方法，并用于分析信号的奇异性和图像的边缘检测，使其成为图像检测的重要工具。在边缘提取中，一般取小波函数为其中为平滑函数，满足，且。则相应的二进小波变换为4.2 小波变换中尺度参数的确定及有效边缘判定

33、进行边缘检测时，滤波尺度参数的选取取决于当前象素点所处的区域是纹理区（平滑区）还是边缘区，如果是纹理区（平滑区）,则取大的滤波尺度,以抑制噪声；如果是边缘区，则取小的滤波尺度，以准确定位边缘。（1）判断纹理区（平滑区）和边缘区计算两个尺度下的小波变换结果和，对于坐标为的象素点，在两个尺度下的小波变换模值分别为和,若满足,则判断为脉冲噪声,其中是阈值，的取值可以去掉一部分噪声和一些弱的细节，其取值受图像中噪声强度的影响，若图像噪声强度较弱，则较小；否则较大。假设图像中噪声分布为，可取,这样可以去掉较多噪声点。将整幅图像中所有这样的脉冲去掉，然后将剩下的,中的局部极值点（候选边缘点）进行线形连接，

34、得到各候选边缘点的链长，将链长 大于某一阈值的候选边缘点定义为大轮廓的边缘点，然后计算其它候选边缘点距离最近的大轮廓边缘点的距离，得到边缘点要素中的参量，候选边缘点的方向参数可使用文献中的方法来计算，这样，可以得到整幅图像中所有候选边缘点的边缘特征。（2）滤波尺度参数的选取范围边缘准确定位时，滤波尺度应尽可能小，取最小尺度。最大滤波尺度的选取应考虑抑噪能力和计算复杂度的折衷,选取能够充分抑制噪声的所有尺度中最小者为。由于使用的小波函数是平滑函数的一阶导数，检测到的边缘点是小波变换的局部极值点，在满足下式时则判定为局部极值点，其中是局部极值点，是的邻域，是阈值。理论上是不需要阈值的，但具体实现时

35、，阈值可以抑制噪声的影响，同时阈值在一定范围内时（不引起有效边缘的丢失），决定着相应的最大滤波尺度，阈值越大， 越小；阈值越小，越大。（3）利用有效边缘测度调整滤波尺度若判定当前象素点处在纹理区（平滑区），则应取最大的滤波尺度；若判定为边缘区，则应取最小滤波尺度；若不能确定时，选取尺度为，这样在一定程度上可以抑制噪声，同时又不会使边缘定位误差过大，如果图像噪声的强度较大，这个尺度可以取，以进一步抑制噪声。第五章 利用小波变换进行边缘检测 5.1利用小波变换进行边缘检测的基本原理 设上一适当光滑的二元函数，并满足条件： 引进记号： 用函数对进行光滑得。现在取二维小波基本函数如下： ， 作相应的小

36、波变换如下： ， 由，的定义，易得： 显然，向量的模取极大值的点对应函数（即按尺度光滑后的图像）的突变点。当很小时，用对光滑化的结果对的突变部位的位置形态影响不大，而当稍大时，则此光滑过程会将的一些细小突变（有噪声引起的不是图像的真正突变）消去，剩下尺寸较大的突变（图像真正的突变）。对参数取二进数列，则二维函数序列为的二维二进小波变换。在尺度下，的梯度的模正比于：它与水平方向的夹角等于：显然，取极大值的点对应的突变点，则给出了梯度方向。根据前面的分析可以知道，用小波系数的局部极大值点可刻画出图像信号突变点的位置，即图像边缘的位置。因此，利用小波变换的极大值检测可以探测出图像的边缘点。5.2 利

37、用小波变换进行边缘提取的两种方法根据上面所述的利用小波变换进行边缘提取的基本原理，对于数字图像可设计出下述计算机实现方法：1、 选择MRA滤波器、边界拓展方式、分解尺度J和阈值T2、 小波变换，设原始图像为，变换得到、3、 计算： 4、 确定与O点相邻的两点、。如图： 5、 同时满足以下三式，则O点为模极大值点：（1）；（2）；（3）或。若则令，否则。6、 将显示出来即为所求边缘图像。7、 重复步骤2-5，直至分解尺度J。 对于上述方法还有以下简化算法，该算法把二维边缘检测问题转化为小波变换来处理，简化了算法的复杂度，其实现方法如下：1、 选择MRA滤波器、边界拓展方式、分解尺度J和阈值T。2

38、、 对图像的每一行像素做一维小波变换，并取绝对值。3、 对图像的每一列像素做一维小波变换重复步骤2。4、 分别求取行列小波变换的模值极大值点，其余点置0。5、 若，则令，否则。6、 将显示出来即为所求边缘图像。7、 重复步骤2-5，直至分解尺度J。按上述方法分别对Lena图进行边缘检测，使用B样条小波，得到下图的边缘图像。 图5-1:用小波变换得到的图像第六章 结 论这种多尺度不仅改变了以前的高斯算法计算量大的缺陷, 而且以较少的运算便可得到理想的图像处理结果. 与传统的高斯尺度空间算法相比较, 得到如下结论:(1) 与传统的经典的边缘检测算法相比, 此小波变换具有定位精确和计算量小的优点;

39、而传统的和经典的边缘检测算法不仅不够精确, 而计算量大, 对噪声的敏感程度也远远高于小波变换的算法;(2) 利用小波的多尺度空间信息, 识别图像边缘特性的能力比高斯尺度空间强. 由于小波变换有多尺度的特点, 可以利用多尺度特性, 通过细节和粗节的逼近, 将图像的空间信息较好的描述, 与传统的和经典的边缘检测算法相比, 在识别边缘特征方面, 强于它们.(3) 在边缘和噪声的取舍中, 由于二者均为高频信号, 很难用频带划分. 小波变换的方法, 使得在大尺度下抑制噪声, 小尺度下, 得到边缘的真实位置; 而传统的和经典的边缘检测算法则在此问题上不能提供有效的解决办法.致 谢：为期三个月左右的毕业设计

40、结束了，在这次设计中使我对数字图像处理有了进一步的认识与研究，充实了自己的专业知识，并且在本次设计中得到了桂志国老师的大力帮助，给我提出了宝贵的意见和修改建议。在此衷心感谢桂志国老师对我的帮助，使我的毕业论文能够顺利完成。参 考 文 献：【1】 李勇，徐震等，MATLAB辅助现代工程数字信号处理 ，西安电子科技大学出版社 2002【2】 飞思科技产品研发中心，MATLAB 6.5辅助小波分析与应用，电子工业出版社 2003【3】 楼顺天， 李博菡 ，基于MATLAB的系统分析与设计信号处理，西安电子科技大学出版社 1998【4】 胡昌华，张军波，夏军，张伟，基于MATLAB的系统分析与设计小波

41、分析 1999【5】 胡昌华，张军波，夏军，张伟，基于MATLAB的系统分析与设计小波分析，西安电子科技大学出版社 1999【6】 杨福生，小波变换的工程分析与应用，科学出版社 2001译文：一种改进的边缘检测与边缘等级排列技术摘要本文说明了一种精简的Canny边缘检测的方法,随后从边界的角度排列了他们相对距离上的边缘.这促进了计算机视觉的合适细节选择,特别是对有连接的物体.第一个议题是对有错误的缝隙与伪边缘中用已经存在的Canny算法只选择一条边缘相互交叉的路径.说明了一种有效,实用的Canny技术处理过程,这种算法能够系统、合适的连接边缘，最后一个命题是边缘的等级排列。1、动机与介绍 这项

42、工作的原理与基础理论需要一个可以连接的具有多种边缘细节的周界提取器。这需要产生一个图像预处理的可以进行手动跟踪系统的平台，为了实现这个结论，使用精确模板的几何学技术来浏览图像。他能够来为图像重要部分的边界提取服务，并能够很好的分辨出周界上的边缘细节的区别于不同。第二部分简明的说明了Canny边缘检测的算法。第三部分说明了周界提取出的边缘排列与缝隙连接。第四部分提出了这种技术市怎样为后期精确的3D物体图像的处理与重现，第五部分是总结与结论。2、Canny边缘检测一种寻找物体边界的标准后简单方法是得到它的轮廓，并且对单一的，已知的可控制的背景轮廓能够很容易的被计算出来。然而轮廓边界的精确位置随着计

43、算阈值的选择，源方向与边界附近的光线水平有关。一种更加精简的方法-几乎可以克服所有的缺点-是用边缘检测，正如我们处理低分辨率图像时，精确的算法用来抵消粗糙的不足，另外一种优点是能够计算出内部细节而不是一个单点，这对多关节的物体追踪是十分重要的。在计算机视觉基础上的边缘检测，被定义用来计算图像的每一个像素值，并按比例计算两种不同边界值的像素Canny边缘检测算法，被用于openCV，选择它是因为他的快速运算结果（与重要的边缘检测技术相比较），与它的高效率实行技术。2.1 如何工作的第一长图是被高斯回旋所平滑的，然后用2D的衍生操作器对平滑图象的亮光部分检测。边缘表示在图像的梯度极大值上，算子跟踪

44、检测那些梯度变化的点，并将不变化的值赋0，所以输出即为一条线。两种阈值限制了这种算法，T1与T2T2，计算时只从值大于T2的点开始。然后想两个方向计算直到值小于T1。这种滞后现象使得噪声不会干扰多边缘的细节。2.2 Y-边缘连接点效果基于Canny操作器的一个问题是：“Y-连接” ：这种情况发生在梯度图象的3个边缘或者更多的边缘相交时。当一个边缘被另外一个物体所融合或交合时，这种算法会把两个边缘认作为同一条线，并且第三条极大值边缘作为一条直线来接近他，但是并没有连接上另外两个边缘。图1：Canny检测算子和它的“Y-连接”效果左：影像图像，右：Canny边缘检测。注意到两条垂直的显示怎么样接近

45、的但是并没有与中心的水平线相接，产生了水平的两个伪边缘。2.3 边缘等级边缘等级与下面的法则有关。0 周界上的边缘1 并不是0等级的并且截止在周界上的边缘。N 不是小于N等级的，并且截止在N-1的边缘后来平台的处理含有图像的各种边缘等级构造起来的框架，有益于选择合适的边缘细节。3.1 周界提取从上面的边缘检测发现,周长被提取出来。虽然比较琐碎，可能会（通常都会）在物体重要部分上有伪造的噪声，而且周长大体上都不再连续。首先图像从内部框架的外边缘进行扫描，由于任何有界图形的像素将会从外部接近（尽管图形没有连接）（实际上，有一种间接的可能性，即这种扫描可能先连接内部物体的重要部分，但事实上并没有干扰

46、到运算的成功。最坏打算它只是对低等级的边缘减慢了检测进程。然后，检测边界的基础图形，当考虑到其他的边界像素时对图形的外界一定要注意，这个可以用下面的实现：1、 当我们对每一个像素沿着周长追踪时，要保持“最后外部像素”的观点。2、 从与痕迹的方向一致的图形最后的最外面的像素开始旋转。换句话说，每一次我们找出另外的一个周长像素, 我们记录最后的非周长像素。(在外部) 为了要找下一个周长像素，我们替换了立即相邻像素（从最后外部像素开始，与旋转的方向相同）。我们假设物体重要的外部边缘在处理的过程中已经被移动或者忽视。尽管它可能与其他边缘相接，也可以用测量长度等技术来实现。这不在是一个没有原形的问题。3

47、.2 缝隙连接现在如果我们试着向前追踪我们已经发现的形状外部，我们会发现通常痕迹只出现在物体周界周围的一部分，因为Canny边缘检测留下了缝隙。图 3: 维持“最后的外部图素”。 周长的每个图素被一个正方形表现，借由从它的中心到它的最后外部的图素中心的一个直线部分,逆时针旋转。这种缝隙以各种不同的形式出现，所以需要解决他们的一般的技术。从图4中拿出一些现实的可观察的数据来，这种技术按照以下各项来解决各种问题：1、 从第一个边界像素开始，按逆时针方向持续追踪。（a）如果痕迹描述一个完全的线路,那么我们完成追踪。 （b）否则纪录结束像素和最后最外部的像素。2、 从起始像素，按照顺时针方向追踪。3、

48、 对于所有的结束点，试着连接他们。（a） 从外部到内部替换，按照一个合适的尺寸的半径来扫描。（b） 如果另外一个片断被缩定，连接到与它相连的最近的像素。（c） 如果这个新的片段已经是在结束的目录上，用他的新的边缘合并它，否则追踪新的片断范围，更新他的结束数据。4、 如果有任何的片断剩余（暗示所有的片断结束在伪边界），试着连接他们。（这一个情形在实际中很少地发生）（a） 背面- 痕迹沿着片段远离结束部分，从外部的缝隙到内部开始测试，否则按照3的步骤进行。（b） 如果没有找到连接，便删除这些结束点（例如噪音）。 图4：连接周界的缝隙。内部的缝隙是相似的。图5：一个实例解释了为什么必须向后旋转，而不

49、能通过边缘向前旋转（会在洞中附着）。3.3 等级排列如果所需要的等级比零等级要好，那么所有的周长等级首先被设置为他们的零等级。然后使用上面的缝隙连接技术（但是不用注意去保持一个外部的意识，并且把所有的求得的半径内的边缘连接起来），所有的邻近边缘被给等级一这一个程序被反复使用直到达到所需要的等级。所有的边缘与其他的现有N等级相连，那些没有被分配一个正确的等级，现在被给予等级N+1。4 有关节物体的3D重建试着用照相机的图像重构一个有关节的物体，对有折痕的关节连接处是非常有用的，这确定了连接点的位置，并且促进了连接点的弯曲点角度变化的察觉。然而，较好的水平细节，在这个阶段通常是不受欢迎的，如同它只

50、是引导噪声或者将不需要的复杂部分成为进一步的处理的平台。然而在关节连接的三维物体，二维图像来看，一些（或更多）的三维物体的外部边缘在投射之后成为内部的边缘。图6：注意手指 (对三度空间的手外部) 的边缘如何变成内在的边缘的在那二维空间图像。在这一个图像中，手指的所有看得见的边缘是等级 3,并且四指的末尾的白色光辉连接处与拇指末尾的关节黑暗连接处的折痕是2或者3。因此等级3是这个图的最佳等级。边缘的两边可以用阈值（T1，T2）精确的设定。这样的话，用这项技术可以正确的给定周长的边缘和等级，而且知道了粗糙物体的位置/方向（比如以前的框架，速度和模型限制），他相对的变为选择那些三维物体的边缘使用技术

51、，和其他的有关节折痕的物体的重要特征。5 摘要与结论这种技术可以简洁高效的论证附加在已存在的边缘检测算子上的错误的实时表现。借由提供这些边缘等级，适当程度的细节部分立刻被带到焦点里去（包括联系的方式，周长提取），因此，直接促进了有关节物体模型的协调，重要的减少了将来计算机视觉处理过程中的冗余。图7：使用这项技术进行图像预处理片断的截图。注意在左边（Canny边缘）上的图像如何被清除，只表现在右边次图像的周长上。 C+源代码可以自由的从网上下载，并且手式图形的追踪报告也可以在那里发现。译文原件：An Improved Edge Detectionand Ranking TechniqueDecl

52、an MurphyComputer Science Dept, Copenhagen Universitydeclandiku.dkAbstractA method is presented for refining the Canny edge detection algorithm, and subsequently ranking the resulting edges relative to their distance (in a topological sense) from the perimeter. This facilitates the selection of the

53、appropriate level of detail for computer vision, particularly of articulated objects. The first issue is the way in which the existing Canny algorithm chooses only one path where edges cross each other, resulting in erroneous gaps and spurious doubling of edges. An efficient pragmatic post-processin

54、g of the Canny technique is presented, which systematically bridges these gaps appropriately. The next issue is the automatic ranking of the edges.1 Motivation and IntroductionThe rationale and motivation behind this work was the requirement for a connected perimeter extractor with a variable amount of edge detail, tracing in from the perimeter. This need aro