图像边缘提取与分割.ppt

第3章 图像边缘提取和分割,3.1引言 3.2 图像分割处理,2019/5/17,第3章 图像边缘提取和分割,3.1引言 图像最基本的特征是边缘，边缘是指其周围像素灰度有阶跃变化或屋顶状变化的那些像素的集合，它存在于目标与背景、目标与目标、区域与区域、基元与基元之间。它是图像分割所依赖的最重要的特征，也是纹理特征中的重要信息源和形状特征的基础。而图像的纹理形状特征的提取又常常要依赖于图像分割。,2019/5/17,3.1.1 统计模式识别简介,统计模式识别认为图像可能包含一个或多个物体，并且每个物体属于若干事先定义的类型、范畴或模式之一。虽然模式识别可以用多种方法实现，但是在此只关心用数字图像处理技术对它的实现。 在给定一幅含有多个物体的数字图像的条件下，模式识别过程如图3.1所示，由三个主要阶段组成。,2019/5/17,3.1.1 统计模式识别简介,图3.1 模式识别的三个阶段,2019/5/17,图像分割：检测出各个物体，把它们的图像和其 余景物分离，这一过程也可以称为图像预处理。,2019/5/17,3.1.1 统计模式识别简介,3.2 图像分割处理,用计算机进行数字图像处理的目的有两个： 一是产生更适合人类视觉观察和识别的图像； 二是希望计算机能够自动进行识别和理解图像。 图像处理的关键问题是对图像进行分解。 分解的最终结果是各种特征的最小成分（基元）。 产生基元的过程就是图像分割的过程。,2019/5/17,3.2 图像分割处理,图像分割也可以按照如下的标准分类： 1.基于区域的分割方法 包括阈值分割法、区域生长和分裂合并法、聚类分割法等； 2.基于边界的分割方法 包括微分算子法、基于区域和边界技术相结合的分割方法。,2019/5/17,3.2.1 基于区域的分割方法,直方图分割（灰度阈值分割） 最简单的方法是建立在灰度直方图分析的基础上。如果一个图像是由明亮目标在一个暗的背景上组成的，其灰度直方图将显示两个最大值，一个是由目标点产生的峰值，另一个峰值是由背景点产生的。,2019/5/17,图 31 图像 的直方图,由直方图可以知道图像 的大部分像素灰度值较低，其余像素较均匀地分布在其他灰度级上。由此可以推断这幅图像是由有灰度级的物体叠加在一个暗背景上形成的。可以设一个阈值 T ，把直方图分成两个部分，如图所示。 T 的选择要本着如下原则:B1 应尽可能包含与背景相关连的灰度级，而B2 则应包含物体的所有灰度级。,当扫描这幅图像时，从 到 之间的灰度变化就指示出有边界存在。当然，为了找出水平方向和垂直方向上的边界，要进行两次扫描。也就是说，首先确定一个门限 ，然后执行下列步骤：,第一，对 的每一行进行检测，产生的图像的灰度将遵循如下规则,式中 是指定的边缘灰度级， 是背景灰度级。,第二，对 的每一列进行检测，产生的图像的灰度将遵循下述规则,在数字图像处理中，样板是为了检测某些不变区域特性而设计的阵列。样板可根据检测目的不同而分为点样板、线样板、梯度样板、正交样板等等。 点样板的例子如图32所示。下面用一幅具有恒定强度背景的图像来讨论。 1）、点样板,3.2.2 基于边界检测方法（样板匹配）,用点样板 的检测步 骤如下：,3-2 点样板检测,样板中心（标号为8）沿着图像从一个像素移到另一个像素，在每一个位置上，把处在样板内的图像的每一点的值乘以样板的相应方格中指示的数字，然后把结果相加。如果在样板区域内所有图像的像素有同样的值，则其和为零。否则其和不为零。,例如，设 代表33模板的权，并使 为模板内各像素的灰度值。从上述方法来看，应求两个矢量的积，即：,线检测样板如图33所示。其中，样板(a)沿一幅图像移动，它将对水平取向的线（一个像素宽度）有最强的响应。对于恒定背景，当线通过样板中间一行时出现最大响应；样板(b)对45方向的那些线具有最好响应；样板(c)对垂直线有最大响应；样板(d)则对-45方向的那些线有最好的响应。,1）、线样板,图33 线样板,设 是图33中四个样板的权值组成的九维矢量。与点样板的操作步骤一样，在图像中的任一点上，线样板的各个响应为 ，这 里 i =1、2、3、4。此处 X 是样板面积内九个像素形成的矢量。给定一个特定的 X ，希望能确定在讨论问题的区域与四个线样板中的哪一个有最相近的匹配。如果第 i 个样板响应最大，则可以断定 X 和第 i 个样板最相近。,换言之，如果对所有的 值，除 外，有：,就 可 以 说 和 第 个 样 板 最 接近 。 如果 ， =2、3、4，可以断定 代表的区域有水平线的性质。,对于边缘检测来说也同样遵循上述原理。通常采用的方法是执行某种形式的二维导数。类似于离散梯度计算，考虑大小的模板，如图34所示。,图34 3 3样板,考虑的图像区域， 及 分别用下式表示,采用绝对值的一种定义为,在 点的梯度为,梯度模板如图35所示。,图35 梯度样板,边缘检测也可以表示成矢量，其形式与线样板检测相同。如果 代表所讨论的图像区域，则：,这里 ， 是图35中的两个样板矢量。 分别代表它们的转置。,这样，梯度公式如下这种形式：,3.2.2 图像分割的一些常用基本方法,原始图像 阈值T=91 阈值T=130 阈值T=43 图3.3 不同阈值对分割结果的影响,2019/5/17,3.Sobel边缘算子 对于阶跃状边缘，Sobel提出一种检测边缘点的算子。对数值图像的每个像素考察它上、下、左、右邻点灰度的加权差，与之接近的邻点权大。据此，定义Sobel算子,2019/5/17,3.2.2 边缘检测,4.Laplacian算子 对于阶跃状边缘，二阶导数在边缘点出现零交叉，即边缘点两旁二阶导数取异号，据此，对数字图像的每个像素，取它关于轴方向和轴方向的二阶差分之和。,2019/5/17,3.2.2 边缘检测,对于屋顶状边缘，在边缘点的二阶导数取极小值。对数字图像的每个像素取它的关于方向和方向的二阶差分之和的相反数，即Laplacian算子的相反数。,2019/5/17,3.2.2 边缘检测,5 Kirsch边缘算子 图3.7所示的8个卷积核组成了Kirsch边缘算子。图像中的每个点都用8个掩模进行卷积，每个掩模对某个特定边缘方向做出最大响应。所有8个方向中的最大值作为边缘幅度图像的输出。最大响应掩模的序号构成了对边缘方向的编码。,2019/5/17,3.2.2 边缘检测,00 450 900 1350,2019/5/17,3.2.2 边缘检测,1800 2250 2700 3150 图3.7 Kirsch边缘算子,2019/5/17,3.2.2 边缘检测,6.Marr-Hildreth 边缘检测算子 Marr-Hildreth边缘检测算子是将高斯算子和拉普拉斯算子结合在一起而形成的一种新的边缘检测算子，先用高斯算子对图像进行平滑处理，然后采用拉普拉斯算子根据二阶微分过零点来检测图像边缘，因此该算子也可称为LOG（Laplacian of Gaussian）算子。,2019/5/17,3.2.2 边缘检测,在数字图像中实现图像与模块卷积运算时，运算速度与选取的模块大小有直接关系，模块越大，检测效果越明显，速度越慢，反之则效果差一点，但速度提高很多。因此在不同的条件下应选取不同大小的模块。在实际计算过程中，还可以通过分解的方法提高运算速度，即把二维滤波器分解为独立的行、列滤波器。常用的55 模块的Marr-Hildreth算子如图3.8所示。,2019/5/17,3.2.2 边缘检测,图3.8 LOG算子的55模板,2019/5/17,3.2.2 边缘检测,3.2.2 边缘检测,7Canny 边缘检测算子 Canny边缘检测算子是近年来在数字图像处理中广泛应用的边缘算子，它是应用变分原理推导出的一种用高斯模块导数逼近的最优算子。通过Canny算子的应用，可以计算出数字图像的边缘强度和边缘梯度方向，为后续边缘点的判断提供依据。,2019/5/17,Canny算子用范函求导方法推导出高斯函数的一阶导数，即为最优边缘检测算子的最佳近似。由于卷积运算可交换，可结合，故Canny算法首先采用二维高斯函数对图像进行平滑，二维高斯函数表示为 其中，为高斯滤波器参数，它控制着平滑的程度，较小的滤波器定位精度高，但信噪比低；较大的滤波器情况正好相反，因此，要根据需要选取高斯滤波器参数。,2019/5/17,3.2.2 边缘检测,传统Canny算法利用一阶微分算子来计算平滑后图像各点处的梯度幅值和梯度方向，获得相应的梯度幅值图像和梯度方向图像，其中，点处两个方向的偏导数 和分别为 则此时点处的梯度幅值和梯度方向分别表示为,2019/5/17,3.2.2 边缘检测,为了精确定位边缘，必须细化梯度幅值图像中的屋脊带，只保留幅值的局部极大值，即非极大值抑制(NMS)。Canny算法在梯度幅值图像中以点为中心的邻域内沿梯度方向进行插值，若点处的梯度幅值大于方向上与其相邻的两个插值，则将点标记为候选边缘点，反之则标记为非边缘点。这样，就得到了候选的边缘图像。,2019/5/17,3.2.2 边缘检测,传统Canny算法采用双阈值法从候选边缘点中检测和连接出最终的边缘。双阈值法首先选取高阈值和低阈值，然后开始扫描图像。对候选边缘图像中标记为候选边缘点的任一像素点进行检测，若点梯度幅值高于高阈值，则认为该点一定是边缘点，若点梯度幅值低于低阈值，则认为该点一定不是边缘点。而对于梯度幅值处于两个阈值之间的像素点，则将其看作疑似边缘点，再进一步依据边缘的连通性对其进行判断，若该像素点的邻接像素中有边缘点，则认为该点也为边缘点，否则，认为该点为非边缘点。,2019/5/17,3.2.2 边缘检测,Canny边缘检测算子的最优性与以下的三个标准有关： （1）检测标准：不丢失重要的边缘，不应有虚假的边缘； （2）定位标准：实际边缘与检测到的边缘位置之间的偏差 最小； （3）单响应标准：将多个响应降低为单个边缘响应。,2019/5/17,3.2.2 边缘检测,3.2.2 边缘检测,Canny边缘检测算子基于如下几个概念： （1）边缘检测算子是针对一维信号表达的，对检测标准和 定位标准最优； （2）如果考虑第三个标准（多个响应），需要通过数值优 化的办法得到最优解。该最优滤波器可以有效地近似 为标准差为的高斯平滑滤波器的一阶微分，为了便于 实现检测误差小于20%，与LOG边缘检测算子很相似； （3）将边缘检测算子推广到二维情况。阶跃状边缘由位置、 方向和可能的幅度来确定。,2019/5/17,3.2.2 边缘检测,8Prewitt算子 Prewitt提出了类似的计算偏微分估计值的方法，梯度计算表示为,2019/5/17,1方向 2方向 3方向 4方向 5方向 6方向 7方向 8方向 图3.9 Prewitt1-8方向边缘检测算子模板,2019/5/17,3.2.2 边缘检测,3.2.2 边缘检测,图 3.10 8个算子样板对应的边缘方向,2019/5/17,3.2.4 边缘检测算子的对比,在数字图像处理中，对边缘检测主要要求就是运算速度快，边缘定位准确，噪声抑制能力强，因此就这几方面对以上介绍的几个算子进行分析比较。首先，在运算速度方面，对于一个图像，其计算量如表3-1所示。,2019/5/17,2019/5/17,根据实际测试结果，简单介绍各个算子的特点。 1.Roberts算子 Roberts算子利用局部差分算子寻找边缘，边缘定位精度较高，但容易丢失一部分边缘信息，同时由于没经过图像平滑计算，因此不能抑制噪声。该算子对具有陡峭的低噪声图像响应最好。,2019/5/17,3.2.4 边缘检测算子的对比,2.Sobel算子和Prewitt算子 Sobel算子和Prewitt算子都是对图像进行差分和滤波运算，差别只是平滑部分的权值有些差异，因此对噪声具有一定的抑制能力，但不能完全排除检测结果中出现伪边缘。同时这2个算子边缘定位比较准确和完整。该类算子对灰度渐变和具有噪声的图像处理结果较好。 3.Krisch算子 该算子对八个方向边缘信息进行检测，因此具有较好的边缘定位能力，并且对噪声有一定的抑制作用，就边缘定位能力和抗噪声能力来说，该算子的处理效果比较理想。,2019/5/17,3.2.4 边缘检测算子的对比,4.Laplacian算子 拉普拉斯算子为二阶微分算子，对图像中的阶跃状边缘点定位准确且具有旋转不变性，即无方向性，但是该算子容易丢失一部分边缘的方向信息，造成一些不连续的检测边缘，同时抗噪声能力比较差。拉普拉斯算子比较适用于屋顶型边缘的检测。,2019/5/17,3.2.4 边缘检测算子的对比,5.Marr-Hildreth算子 Marr-Hildreth算子首先通过高斯函数对图像作平滑处理，因此对噪声的抑制作用比较明显，但同时也可能将原有的边缘也平滑了，造成某些边缘无法检测到。此外高斯函数中方差 参数的选择，对图像边缘检测效果有很大的影响。 越大，检测到的图像细节越丰富，但对噪声抑制能力相对下降，易出现伪边缘；反之，则抗噪声性能提高，但边缘定位准确性下降，易丢失许多真边缘，因此，对于不同的图像应该选择不同的参数。,2019/5/17,3.2.4 边缘检测算子的对比,6.Canny算子 Canny算子同样采用高斯函数对图像作平滑处理，因此具有较强的去噪能力，但同样存在