《形状分析与描述》PPT课件.ppt

第五章形状分析与描述,形状是由组成物体的轮廓线或物体表面的所有点相对位置决定。这说明我们只能通过物体的轮廓线或外表面才可知道其形状，而外轮廓和外表面是能为视觉所感知的。利用边缘检测和图象分割，能够提取关于景物的重要的形状信息。计算机视觉的另一主要难题是形状的表示。只有通过表示，才能对感兴趣的景物形状进行学习、匹配、重构与利用。把边缘连接起来就成为轮廓。轮廓可以是断开的，也可以是封闭的。封闭轮廓对应于区域的边界，而区域内的象素可以通过填充算法来填满。断开的轮廓可能是区域边界的一部分，也可能是图象线条特征。区域之间的对比度太弱或边缘检测阈值设置太高都有可能产生间断的轮廓。轮廓可以用边缘序列表或曲线来表示。曲线通常称为轮廓的数学模型。曲线表示包括线段、二次曲线、三次样条曲线等。,轮廓表示的评价标准：简单：轮廓应该是一种简洁的表示。精确：轮廓应能精确地逼近图象特征。有效：轮廓应适合于后处理阶段的计算。决定轮廓表示精确性的主要因素有以下三个方面：（1）用于轮廓建模的曲线形式；（2）曲线拟合算法的性能；（3）边缘位置估计的精度。轮廓的最简单表示形式是边缘有序表。这种表示的精度就是边缘估计的精度，但其表示的紧凑性是最差的，后处理也不方便，因此不是一种高效的图象分析方法。用适当的曲线模型来拟合边缘会提高精确度，因为曲线模型拟合边缘时往往具有均值化效应，因此可以减少边缘位置误差。曲线模型也会提高轮廓表示的经济性，为后处理提供了一种更简单、更紧凑的表示。,已知一组控制点，曲线拟合常采用内插曲线或逼近曲线来实现。内插：指使得拟合曲线通过所有的控制点。逼近：指使拟合曲线非常接近这些控制点，而无需一定通过这些控制点。平面曲线函数可表示为三种形式：（1）显式；（2）隐式：；（3）参数式：，其中u是某一参数；函数的显式表示很少用在计算机视觉中，主要原因是平面上的曲线可能卷曲，致使一个x值可能对应曲线上多个y值。,5.1数字曲线及其表示,下面讨论一组计算曲线几何元素的算法，包括轮廓长度、正切方向、曲率等。由于相邻象素间的量化增量是45，因此，精确计算斜率和曲率是很困难的。估计正切方向的基本思路是使用边缘表中非邻接的边缘点，这就允许存在一个较大可能的正切方向集合。设是边缘表中第i个边缘坐标。K斜率是在边缘表中相距K个边缘点的两个边缘点之间的方向矢量。进一步又分为左K斜率和右K斜率。K曲率是左、右K斜率之差。假定边缘表中有n个边缘。则数字曲线的长度S及轮廓端点之间的距离D可表示为：,一、链码链码是沿着轮廓记录边缘表的一种表示方法。分为4方向链码和8方向链码。如用8邻点链码表示一条曲线，即从边缘表中的第一个边缘点开始，沿着轮廓按逆时针方向行走，行走方向用8个链码中一个表示。,下图所示曲线的链码是：602222202101344444454577012其差分链码是：22000062771210000017120111,将上页图中曲线旋转90后如上图。曲线的链码是：024444424323566666676711234其差分链码不变。二、斜率表示法,5.2曲线拟合,常用的曲线模型有：直线段、圆锥曲线和三次样条曲线。一般，拟合之前应考虑如下两个问题：（1）用什么方法进行边缘点曲线模型拟合？（2）如何测量拟合的逼近程度？现假设边缘位置足够精确，不会对拟合结果产生影响。以下讨论用曲线模型拟合边缘点的方法。设di是边缘点到一条拟合曲线的距离（含正负号），在曲线同一侧时具有相同的符号。以下是一些常用的用于衡量曲线拟合效果的方法。（1）最大绝对误差（MAE）（2）均方差（MSE）（3）规范化最大误差（4）误差符号变化次数（5）曲线长度与端点距离之比,一、多直线段多直线段是指端点连接端点的直线段序列，直线段序列的连接点称为顶点。多线段算法的输入值是边缘点有序表拟合边缘表并把第一个边缘点和最后一个边缘点连接起来的直线段公式如下：上式可改写成：其中：而是边缘点和之间的距离。任给一点，设，则r的符号可用来计算符号变化次数。点与拟合直线段的距离为：规范化最大误差为：,（1）多直线段分裂自顶向下的分裂算法是将整条曲线作为初始曲线，通过反复增加顶点来用直线段拟合曲线。直到所有的直线段对应的规范化最大误差均小于某一阈值为止。该过程也称为迭代分解。,（2）线段合并线段合并是指用一直线段尽量多地拟合边缘表中的边缘点。当边缘点离直线太远而无法用该直线段拟合时，则开始新的直线段拟合。（自底而上合并的多线段拟合方法）,（3）分裂与合并将多直线分裂与线段合并方法组合起来，形成合并与分裂算法。,二、二次曲线二次曲线的一般表示如下：二次曲线也称为圆锥曲线，包括：圆、椭圆、抛物线、双曲线。（1）圆弧段（2）圆锥曲线,5.3样条曲线,样条：富有弹性的细长条。样条曲线：将样条上的若干点固定，沿样条画出的光滑曲线。在数学意义上，样条曲线是用分段多项式表示的一个函数，在其连接点处具有连续的一阶和二阶导数。样条曲线有很多应用。在数学分析中，当没有合适的函数模型时，可选用样条函数拟合数据点；在计算机图形学和计算机辅助设计中，样条函数用来表示自由曲线；在计算机视觉中，若没有表示曲线的合适模型时，样条函数可以提供曲线的通用表示形式。需指出，几何等效和参数等效是两个不同的概念。几何等效：是指它们连接相同的点集（即在空间上对应着相同的形状）。参数等效：是指两条曲线的方程相同。显然，参数等效比几何等效更稳定。两条曲线可以是几何上等效但可具有不同的参数表示式，这是计算机视觉中的一个重要概念。在计算机视觉的形状表示和物体识别中，常常基于几何等效性。,一、三次样条曲线样条函数最常见的形式是三次样条函数，它是分段三次多项式的一个序列。直线段、二次曲线序列都是样条函数的特例。三次样条函数可以用很少的几个样条段表示很复杂的曲线。已广泛用于图形学及轮廓表示。三次样条具有足够的自由度来逼近边缘段位置和方向。大多数边缘检测算子同时提供边缘方向和位置估计。在直线段、二次曲线拟合中，仅使用了边缘的位置信息。下面介绍一种在三次样条曲线拟合中如何使用有边缘检测器产生的方向信息的例子。平面三次曲线方程如下：或：参数u取值范围在0和1之间。三次曲线起始点为，终点为。三次样条是由构成的一个序列。这一序列定义在连续区间0，1，1，2，,n-1，n上，并将端点连接起来使得在端点处。样条中每一个三次曲线段称为样条段，连接样条段两端的端点称为结点。设第i个结点连接第i-1样条段和第i样条段，该结点表示为或。与前述曲线拟合算法相同，把边缘点序列分成一个个子序列，每一子序列的第一个和最后一个边缘点为样条曲线的结点，然后再用样条段拟合这些结点。由三次曲线方程可知，样条中每一个三次曲线段都需要确定四个二维矢量共计8个参数，其中，曲线段的两个端点提供4个约束，结点处的一阶连续性提供2个约束，结点处的二阶连续性提供两个约束，结点处的方向信息仅提供1个约束（由于结点由两个样条段共享），这样产生的方程数量为9个，多于三次样条段所需的8个参数。,在结点处光滑连接样条段是非常重要的。在计算机图形学中，光滑连接是通过增加二阶连续性来实现。由上述分析知，二阶连续性提供2个约束，从而对样条段产生过约束。为避免过约束，同时又要使结点处光滑，可以采用结点处二阶不连续性的极小化条件，也即将结点处的曲率差值极小化作为一个约束代替二阶连续性提供的两个约束。（推导略）同前面介绍的多线段、二次曲线拟合算法一样，结点必须从边缘表中选出。调节结点的位置和数量可以改善三次样条对整个边缘点集的拟合效果。三次样条拟合算法仅需要求解一个小的线性系统就可得到正切值的符号和量值，因此该算法十分有效。另外，也可以使用交互式图形界面，在其上可方便地调节三次样条曲线。,二、B样条曲线B样条曲线是由结点引导的逐段多项式曲线，是一种平滑和内插技术。,5.4Hough变换,Hough变换(HT)是一种用于区域边界形状描述的方法，经典HT常常被用于直线段、圆和椭圆的检测。HT是于1962年由P.V.C.Hough提出的，后经不断改进。广义HT可推广至检测任意形状。无论是HT还是广义HT其基本思想是将图象空间变换到参数空间，用大多数边界点满足的某种参数形式来描述图象中的曲线（边界）。由于HT是根据局部度量来计算全面描述参数，因此对区域边界被噪声干扰或被其它目标遮盖而引起边界发生某些间断的情况，它具有很好的容错性和鲁棒性。,一、Hough变换的原理Hough（HT）变换是一种用于区域边界形状描述的重要方法。采用Hough变换检测任意曲线的原理如下：假设为需检测曲线的参数方程。式中为形状参数，x，y为空间域的图象点坐标。对于图象空间的任一点，利用上式可将其变换为参数空间中的一条曲线。假定空间域中位于同一曲线上的n个点，对这n个点逐一进行上述变换，则在参数空间中对应地得到n条曲线，由上式知，这n条曲线必定经过同一点，找到参数空间的这个点就决定了空间域中的曲线l。传统的Hough变换将空间域中的每一个轮廓点代入上式，其计算结果对参数空间中的量化点进行投票，若票数超过某一门限值，则认为有足够多的图象点位于该参数点所决定的曲线上。即需要逐点投票、记录，故耗时长，占用存储量也大。为克服这一缺点，人们在应用中提出了许多改进算法。,二、线段检测直线的方程为：如下图，按上述思想，一种检测直线的简单方法为：首先把参数平面离散化，并建立一个参数矩阵。对于图象空间的每一个边缘点，建立方程，并对离散化后的每个a值，计算出相应的b值，然后，将参数矩阵的元素的值加1：,即重复这一过程，直到扫描完所有的边界点。在过程结束后，参数矩阵元素的值表示图象空间中满足方程的边界点的个数，如果其大于某一阈值，就表示检测到了相应的直线。这种在参数空间进行“投票”的方法体现了Hough变换抗干扰的鲁棒性。上述参数方程不适合处理垂直直线，因为此时直线的a值趋于无穷大。为此，需引进直线的极坐标形式：,空间与空间的变换,y,x,此时的参数空间为。其检测步骤如下：（1）适当量化参数空间；（2）累加器清零；（3）对每一投票，相应累加器加1；（4）对投票结果进行阈值化处理。由上述分析可知，若对参数空间量化过细，则计算量增大；反之，若量化过粗，则参数空间的集聚效果差，检测精度降低。因此，在应用Hough变换时，应根据实际选取合适的量化值。如果图象空间各点的梯度方向已知，在寻求直线边缘时，可在边缘点梯度方向的一定范围内对精细量化，其它角则粗量化，这样可提高检测直线方向角的精度，又不致增加计算量。,例：,三、圆检测对于圆，其参数方程为：其参数空间为，增加到三维。,其算法步骤如下（从灰度图象出发）：（1）应用边缘检测算子提取图象边缘，并按一定的阈值对其进行二值化处理，生成边缘图象；（2）圆的参数方程可改写为：将角度值按参数空间的大小离散化，并求出相应的值存入数据表中；（3）对边缘图象中的所有点，当r的取值在变化时，求出值，并将对应的累加器阵列中的单元加1；（4）对累加阵进行处理，当时，其参数对应为图象空间的圆形边界。显然，上述圆的Hough变换的计算量是非常大的。为降低计算复杂性，有许多改进算法，如增加边缘方向信息等。,四、椭圆检测参数空间变为5维：五、广义Hough变换与任意形状检测,5.5傅立叶描述子,由于沿着封闭轮廓的位置函数是周期性的，因此傅立叶级数可以用来逼近轮廓，轮廓逼近的分辨率取决于傅立叶级数的项数。假定物体的边界表示为一个坐标的序列，其中。用复数来表示每一个坐标对，即即，x轴为实轴，y轴是复数序列的虚轴。对于封闭边界，这一序列是周期的，其周期是N，这样，边界就可在一维空间上表示。一维序列函数u（n）的离散傅立叶变换（DFT）定义为：复数系数a（k）称为边界的傅立叶描述子。,傅立叶描述子是轮廓表示的简洁表示方法。然而，更简洁的表示形式是使用傅立叶序列的低阶项的低分辨率逼近。如果仅使用前M个系数，则对u（n）的逼近如下：尽管仅使用M项来求边界u（n）的每一项，n仍然取0到N-1的区间。换句话说，被逼近的边界具有相同的点数，但重建每一个点的项数不同。边界的简单几何变换，如平移、旋转和尺度变换，将关系到边界傅立叶描述子的简单操作。这使得傅立叶描述子在边界匹配中的应用非常具有吸引力。然而，傅立叶描述子在描述具有遮挡物体时的确存在一些问题。使用其它边界表示也可以得到类似的描述子。傅立叶描述子已成功地应用于医学图象分析领域。,上机作业：编程实现用Hough变换检测图象中直线的算法，要求检测出所有长度大于20象素的直线段。,