数字图象处理 第7章 图像分割与特征提取.ppt

第七章图像分割、特征提取与描述,主要内容7.1引言7.2图像分割7.3图像的特征7.4图像的描述,7.1引言,图像分割是将图像划分为若干互不相交的小区域的过程。小区域是某种意义下具有共同属性的象素的连通集合，如物体所占的图像区域、天空区域等。,7.1引言,图像分割是图像理解的基础，而在理论上图像分割又依赖图像理解，彼此是紧密关联的。图像分割在一般意义下是十分困难的问题，目前的图像分割一般作为图像的前期处理阶段，是针对分割对象的技术，是与问题相关的，如最常用到的利用阈值化处理进行的图像分割。,7.1引言,图像分割有二种不同的途径：区域法：将各象素划归到相应物体或区域的象素聚类方法边界方法：通过直接确定区域间的边界来实现分割,7.2图像的分割,7.2.1、并行区域技术基于阈值的图像分割技术7.2.2、串行区域技术7.2.3、并行边界技术基于边界的图像分割,7.2图像的分割,7.2.1、并行区域技术基于阈值的图像分割技术7.2.2、串行区域技术7.2.3、并行边界技术基于边界的图像分割,7.2.1并行区域技术基于阈值,取阈值是最常见的并行的直接检测区域的分割方法。阈值是在分割时作为区分物体与背景象素的门限，大于或等于阈值的象素属于物体，而其它属于背景。这种方法对于在物体与背景之间存在明显差别（对比）的景物分割十分有效。实际上，在任何实际应用的图像处理系统中，都要用到阈值化技术。为了有效地分割物体与背景，人们发展了各种各样的阈值处理技术，包括全局阈值、自适应阈值、最佳阈值、局部阈值等。,7.2.1并行区域技术基于阈值,单阈值分割只用一个阈值分割。多阈值分割用多个阈值分割。在一般的多阈值情况下，多阈值分割取为：,7.2.1并行区域技术基于阈值,（1）极小值点阈值对于双峰直方图，选取两个峰之间的谷对应的灰度值作为阈值。将直方图的包络看作一条曲线，求直方图包络线的极小值点对应的灰度值作为阈值。,7.2.1并行区域技术基于阈值,在阈值化处理之后，可以通过直接跟踪物体边界的方法将物体区域分割出来，得到其轮廓并进一步分析其几何形状特征。,7.2.1并行区域技术基于阈值,7.2.1并行区域技术基于阈值,（2）最佳阈值有时目标和背景的灰度值有部分交错，用1个全局阈值不能将它们绝然分开。希望减小误分割的概率，选取一个最佳阈值。阈值的选择需要根据具体问题来确定，一般通过实验来确定。,7.2.1并行区域技术基于阈值,（3）直方图变换利用象素邻域的局部性质变换原来的直方图得到一个新的直方图。新的直方图或者谷更深或者谷变成峰，更容易检测。常用的方法是用象素的梯度值。,7.2.1并行区域技术基于阈值,（4）灰度和灰度平均图物体与背景各自的灰度都较均匀二者相差不大时，分割时可以这样处理：横轴取象素的灰度r，纵轴取r的邻域的平均。直方图为坐标点上的象素数目。边界上的点将远离对角线，因此选远离对角线的点的灰度作为分割的灰度门限将获得较好的分割效果。,7.2.1并行区域技术基于阈值,（5）灰度值和梯度值散射图做直方图（r，g）横轴取象素的灰度r，纵轴取r的梯度g。直方图（r，g）为具有某个灰度和梯度值的象素数目。,7.2图像的分割,7.2.1、并行区域技术基于阈值的图像分割技术7.2.2、串行区域技术7.2.3、并行边界技术基于边界的图像分割,7.2.2串行区域技术,一、区域生长二、分裂合并,一、区域生长区域生长方法是根据同一物体区域内象素的相似性质来聚集象素点的方法，从初始区域（如小邻域或甚至于每个象素）开始，将相邻的具有同样性质的象素或其它区域归并到目前的区域中从而逐步增长区域，直至没有可以归并的点或其它小区域为止。区域内象素的相似性度量可以包括平均灰度值、纹理、颜色等信息。关键是相似性度量准则以及初始区域或象素的确定。,7.2.2串行区域技术,7.2.2串行区域技术,图像：5586取跟踪门限T25586489748972283228333333333生长准则：区域任一象素与其邻点灰度差T结果与起始点选择和门限选择有关,7.2.2串行区域技术,图像：558655865586489748974897228322832283333333333333生长准则：待检测象素的灰度与已检测的区域的平均灰度差T例：取跟踪门限T2结果与起始点选择(如选6)和门限选择有关,7.2.2串行区域技术,二、区域分裂与合并从整个图像开始不断分裂得到各个区域。实际中常常先把图像分成任意大小且不重叠的区域，然后再合并或分裂这些区域以满足要求。,7.2.2串行区域技术,确定均匀性准则。例如以一个区对某种特征（如灰度、彩色或纹理）的均匀性为准则，通常用门限T来约束。对满足均匀性准则的小区则合并，不满足均匀性准则的小区则采用四叉树的方式将此区等分为四个小区。,7.2.2串行区域技术,7.2.2串行区域技术,7.2图像的分割,7.2.1、并行区域技术基于阈值的图像分割技术7.2.2、串行区域技术7.2.3、并行边界技术基于边界的图像分割,7.2.3并行边界技术基于边界的图像分割,边缘检测是所有基于边界的分割方法的第一步。视觉系统对图像的边缘更敏感，而不是根据点的灰度区分出物体，人对边界的识别机理也不是设置一个灰度门限T来分割物体的。图像增强讨论了用梯度、拉普拉斯算子以及高通滤波增强图像边缘轮廓的方法，实际上，对增强后的图像边缘轮廓进行门限化处理，就可以用于边缘检测。,7.2.3并行边界技术基于边界的图像分割,一、梯度边缘检测当物体与背景有明显对比度时，物体的边界处于图像梯度最高的点上，通过跟踪图像中具有最高梯度的点的方式获得物体的边界，可以实现图像分割。这种方法容易受到噪声的影响而偏离物体边界，通常需要在跟踪前对梯度图像进行平滑等处理，再采用边界搜索跟踪算法来实现。,7.2.3并行边界技术基于边界的图像分割,梯度图像阈值化梯度图像,7.2.3并行边界技术基于边界的图像分割,二、边界提取与轮廓跟踪为了获得图像的边缘人们提出了多种边缘检测方法。在边缘图像的基础上，需要通过平滑等处理去除噪声点、毛刺、空洞等不需要的部分，再通过细化、边缘连接和跟踪等方法获得物体的轮廓边界。例：采用光栅跟踪（顺序跟踪）的方法,7.2.3并行边界技术基于边界的图像分割,步骤：1）先用高阈值Th（检测阈值）检出图像梯度轨迹曲线的初选点。2）从第一行的初选点开始用较低阈值Tl（跟踪阈值）对初选点的下三个点进行判断。在阈值范围之内的接收，反之去除。3）对应于某个检出点，如果下一行无可接收的对象，则这条曲线跟踪即可结束。,7.2.3并行边界技术基于边界的图像分割,步骤：4）对应于某个检出点，如果下一行有多个可接收的对象，则这条曲线发生分支。跟踪对各个分支同时进行。5）对于不在第一行的其他初选出来的象素，从该点开始，重新使用跟踪阈值进行跟踪，以检出不是从第一行开始的其他曲线。,7.2.3并行边界技术基于边界的图像分割,原图像用阈值Th7处理的结果,7.2.3并行边界技术基于边界的图像分割,用检测阈值Th7和跟踪阈值Tl4处理的结果,7.2.3并行边界技术基于边界的图像分割,检测与跟踪时的准则不一定是灰度，也可以是梯度等反映局部性质的量。光栅顺序跟踪跟光扫描方向有关，最好用其他方向再跟踪一次，例如逆序向上，两种方法综合起来能得到更好的效果。若边缘和扫描方向平行时效果不好，最好在垂直方向跟踪一次。,7.2.3并行边界技术基于边界的图像分割,三、HOUGH变换HOUGH变换是利用图像全局特性而将边缘象素连接起来组成区域封闭边界的一种方法。在预先知道区域形状的情况下，用哈夫变换可以很方便的得到边界曲线而将不连续的边缘象素点连接起来。哈夫变换的主要优点：受噪声和曲线间断的影响较小。利用哈夫变换可以检测图像中某些符合参数模型的主导特征，如直线、圆、椭圆等，可以通过对其参数进行聚类的方法，抽取相应的特征。适用于二值图像。,7.2.3并行边界技术基于边界的图像分割,哈夫变换的基本思想是：点线对偶性。在图像空间XY中，所有过点的直线满足直线方程：p为斜率，q为截距所以：可以认为该式是参数空间PQ中过点(p,q)的一条直线。,这条直线上所有点的p，q相同,7.2.3并行边界技术基于边界的图像分割,显然，XY空间中的共线的点对应PQ空间中一组相交的线，这组线的交点就是p，q。反过来，在参数空间相交于同一个点的所有直线在图像空间中都有共线的点与之对应。这就是点线对偶性。哈夫变换就是根据这些关系把图像空间中的检测问题转换到了参数空间中，通过在参数空间中进行简单的累加统计完成检测任务。,7.2.3并行边界技术基于边界的图像分割,具体实现时，考虑p、q的可能取值范围，从大到小进行累加。1)初始化A(p,q)=02)对XY空间中的每一点，P取遍所有可能值3）算出对应的q，给A(p,q)单元加14）根据maxA(p,q)，求出直线,7.2.3并行边界技术基于边界的图像分割,问题：p、q的取值较大时，计算量增加可以采用极坐标系形式。直线在极坐标系下具有如下的参数方程形式：这条直线在极坐标下表示为一个点,7.2.3并行边界技术基于边界的图像分割,7.3图像的特征,图像特征是图像的重要属性，它一般应具备4个特点：可区分性、可靠性、独立性、数量少图像的特征主要有：1)颜色（灰度）特征2)形状特征3)纹理特征4)空间关系特征,7.4图像的描述,边界描述区域或闭合边界的描述,7.4.1图像的边界描述,一、周长和边长L周长的定义通常有3种：1）把区域像元看为正方形，边界只有水平和垂直两种线，相当于像元小方格的四个边的四邻码2）把像元看做一个点，可用弧长的8邻码计算3）用边界所占面积,7.4.1图像的边界描述,方法1：L25d；d为小正方形的边长方法2：方法3：18,其中为小正方形的面积,7.4.1图像的边界描述,二、边界链码链码是从物体边界上的任意一个点（x,y）出发，而后只记录下一个边界点的方向码，直至回到出发点为止的编码序列。在采用边界跟踪方法获取物体边界轮廓时，可以直接得到边界链码。,7.4.1图像的边界描述,8方向链码,7.4.1图像的边界描述,下面的图像用链码表示为：20206644这个结果与起点有关，为了避免这个问题，用循环码（归一化）表示：02066442具体方法为：首尾相连循环形成的自然数最小,7.4.1图像的边界描述,图像平移链码不变，但旋转时链码会发生变化，为此可用差分码表示原链码(4)20206644原链码(3)17175533差分码66266060（模8、反向）差分码66266060,1,1,3,3,7,5,5,7,7.4.1图像的边界描述,三、线条的描述若图中相继相邻的点集能拟合与某些曲线函数，那么这些曲线函数的少量参数就可做为该点集的一种描述。如：迭代拟合、B样条曲线拟合等。,7.4.1图像的边界描述,迭代拟合：把曲线用分段线段近似表示出来。首先用直线连接端点A和B，然后选取到直线AB距离最远的点C，如果点C偏离AB超过了某种限度，则消去线段AB，然后分别连接AC和BC。重复上述的步骤，直到偏离值小于原先设定的限度为止。,四、傅里叶描述子用来描述物体特征的少量数据叫做描述子。描述子应满足四点要求：1）能精确地描述特征2）对大小变化不敏感（尺度变化不变性）3）对描述的起点不敏感4）对平移旋转不敏感,7.4.1图像的边界描述,7.4.1图像的边界描述,四、傅立叶描述子从任一点开始的一个N点的闭合边界序列（xk,yk）(k=0,N)可以认为是一个周期为N的周期函数，把它放到复平面上，就形成了一个1维的复数序列uk+jvk,7.4.1图像的边界描述,对这个复数序列进行傅立叶变换,7.4.1图像的边界描述,我们只利用s(w)的前M个系数，其余系数置0，反变换就得到s(k)的一个近似傅立叶变换的低频信息对应总体形状而高频信息对应一些细节，这样我们可用对应低频分量的傅立叶系数（个数M远N）来近似描述边界的形状。（P231图）,7.4.1图像的边界描述,7.4.2区域的描述,一、矩描述子对数字图像f(x,y)，如果它分段连续且在XY平面上只有有限个点不为零，可证明其各阶矩存在，f(x,y)的pq阶矩定义为：,7.4.2区域的描述,可以证明f(x,y)与mpq是一一对应的。区域的矩是用所有区域的点计算出来的，抗干扰性能较好。f(x,y)的pq阶中心矩定义为：其中：,7.4.2区域的描述,三阶以下的中心矩：u00、u01u10、u11、u20、u02、u30、u12、u21、u03f(x,y)的归一化中心矩可表示为：其中,7.4.2区域的描述,可以导出具有平移、旋转和尺度变换不变性的7个不变矩：Matlab示例,7.4.2区域的描述,二、区域面积通常面积的计算可以直接对分割出的物体内像元计数即可。从闭合边界来计算面积有类似计算周长的4连通和8连通的问题用不同方法描述的边界计算面积时有很大差异。,7.4.2区域的描述,A31（每一小块面积认为是1）,7.4.2区域的描述,第一层有1个白小块(标记为1的)，面积4m-1=43-1=16第二层有2个白小块，面积24m-2=243-2=8第三层有11个白小块，面积114m-3=1143-3=11总面积1681135,用四叉树描述区域时，只需将叶节点的面积相加即可。叶节点面积为4ml，l为叶节点所在层数，m为图像层数，N2m,7.4.2区域的描述,三、区域重心重心是一种全局描述符，区域重心的坐标是根据所有属于区域的点计算出来的。在区域本身尺寸与各区域间的距离相对很小时，可将区域用其重心坐标的质点来代替,7.4.2区域的描述,四、区域的拓扑结构拓扑学是研究图形不受畸变变形影响性质的。欧拉数EC（连通组元个数）H（孔）E=1-2=-1E=1-1=0,B,A,7.4.2区域的描述,五、区域的纹理纹理是物体表面结构的模式，如木材、水泥、沙子、玻璃等等，对图像分析和理解十分重要。纹理可以看做许多基本纹理小单元的重复，这些小单元可称作纹理元，纹理元常包含多个象素。常用的纹理