图像分割与特征提取

图像分割与特征提取第1页，共73页，2023年，2月20日，星期四8.1引言图像分割是将图像划分为若干互不相交的小区域的过程。小区域是某种意义下具有共同属性的象素的连通集合，如物体所占的图像区域、天空区域等。连通的概念是指集合中任意两个点之间都存在着完全属于该集合的连通路径。对于离散图像而言，连通有4连通和8连通之分

4连通8连通第2页，共73页，2023年，2月20日，星期四8.1引言图像分割是图像理解的基础，而在理论上图像分割又依赖图像理解，彼此是紧密关联的。图像分割在一般意义下是十分困难的问题，目前的图像分割一般作为图像的前期处理阶段，是针对分割对象的技术，是与问题相关的，如最常用到的利用阈值化处理进行的图像分割。

第3页，共73页，2023年，2月20日，星期四8.1引言图像分割有二种不同的途径：区域法：将各象素划归到相应物体或区域的象素聚类方法边界方法：通过直接确定区域间的边界来实现分割第4页，共73页，2023年，2月20日，星期四8.2图像的分割8.2.1、并行区域技术－基于阈值的图像分割技术8.2.2、串行区域技术8.2.3、并行边界技术－基于边界的图像分割第5页，共73页，2023年，2月20日，星期四8.2图像的分割8.2.1、并行区域技术－基于阈值的图像分割技术8.2.2、串行区域技术8.2.3、并行边界技术－基于边界的图像分割第6页，共73页，2023年，2月20日，星期四8.2.1并行区域技术－基于阈值

取阈值是最常见的并行的直接检测区域的分割方法。阈值是在分割时作为区分物体与背景象素的门限，大于或等于阈值的象素属于物体，而其它属于背景。这种方法对于在物体与背景之间存在明显差别（对比）的景物分割十分有效。实际上，在任何实际应用的图像处理系统中，都要用到阈值化技术。为了有效地分割物体与背景，人们发展了各种各样的阈值处理技术，包括全局阈值、自适应阈值、最佳阈值、局部阈值等。第7页，共73页，2023年，2月20日，星期四8.2.1并行区域技术－基于阈值单阈值分割 只用一个阈值分割。多阈值分割用多个阈值分割。在一般的多阈值情况下，多阈值分割取为：第8页，共73页，2023年，2月20日，星期四8.2.1并行区域技术－基于阈值

（1）极小值点阈值对于双峰直方图，选取两个峰之间的谷对应的灰度值作为阈值。

将直方图的包络看作一条曲线，求直方图包络线的极小值点对应的灰度值作为阈值。第9页，共73页，2023年，2月20日，星期四8.2.1并行区域技术－基于阈值

在阈值化处理之后，可以通过直接跟踪物体边界的方法将物体区域分割出来，得到其轮廓并进一步分析其几何形状特征。

第10页，共73页，2023年，2月20日，星期四8.2.1并行区域技术－基于阈值

第11页，共73页，2023年，2月20日，星期四8.2.1并行区域技术－基于阈值（2）最佳阈值

有时目标和背景的灰度值有部分交错，用1个全局阈值不能将它们绝然分开。希望减小误分割的概率，选取一个最佳阈值。阈值的选择需要根据具体问题来确定，一般通过实验来确定。第12页，共73页，2023年，2月20日，星期四8.2.1并行区域技术－基于阈值（3）直方图变换 利用象素邻域的局部性质变换原来的直方图得到一个新的直方图。新的直方图或者谷更深或者谷变成峰，更容易检测。常用的方法是用象素的梯度值。第13页，共73页，2023年，2月20日，星期四8.2.1并行区域技术－基于阈值（4）灰度和灰度平均图物体与背景各自的灰度都较均匀二者相差不大时，分割时可以这样处理：横轴取象素的灰度r，纵轴取r的邻域的平均。直方图为坐标点上的象素数目。边界上的点将远离对角线，因此选远离对角线的点的灰度作为分割的灰度门限将获得较好的分割效果。第14页，共73页，2023年，2月20日，星期四8.2.1并行区域技术－基于阈值（5）灰度值和梯度值散射图做直方图（r，g）横轴取象素的灰度r，纵轴取r的梯度g。直方图（r，g）为具有某个灰度和梯度值的象素数目。第15页，共73页，2023年，2月20日，星期四8.2图像的分割8.2.1、并行区域技术－基于阈值的图像分割技术8.2.2、串行区域技术8.2.3、并行边界技术－基于边界的图像分割第16页，共73页，2023年，2月20日，星期四8.2.2串行区域技术一、区域生长二、分裂合并第17页，共73页，2023年，2月20日，星期四一、区域生长

区域生长方法是根据同一物体区域内象素的相似性质来聚集象素点的方法，从初始区域（如小邻域或甚至于每个象素）开始，将相邻的具有同样性质的象素或其它区域归并到目前的区域中从而逐步增长区域，直至没有可以归并的点或其它小区域为止。区域内象素的相似性度量可以包括平均灰度值、纹理、颜色等信息。关键是相似性度量准则以及初始区域或象素的确定。8.2.2串行区域技术第18页，共73页，2023年，2月20日，星期四8.2.2串行区域技术图像：5586取跟踪门限T＝25586489748972283228333333333生长准则：区域任一象素与其邻点灰度差<T结果与起始点选择和门限选择有关第19页，共73页，2023年，2月20日，星期四8.2.2串行区域技术图像：558655865586489748974897228322832283333333333333生长准则：待检测象素的灰度与已检测的区域的平均灰度差<T例：取跟踪门限T＝2结果与起始点选择(如选6)和门限选择有关第20页，共73页，2023年，2月20日，星期四8.2.2串行区域技术二、区域分裂与合并

从整个图像开始不断分裂得到各个区域。实际中常常先把图像分成任意大小且不重叠的区域，然后再合并或分裂这些区域以满足要求。第21页，共73页，2023年，2月20日，星期四8.2.2串行区域技术确定均匀性准则。例如以一个区对某种特征（如灰度、彩色或纹理）的均匀性为准则，通常用门限T来约束。对满足均匀性准则的小区则合并，不满足均匀性准则的小区则采用四叉树的方式将此区等分为四个小区。第22页，共73页，2023年，2月20日，星期四8.2.2串行区域技术10008777111188880000788800017887443378103344881087111111871111011000877711118888000078880001788744337810334488108711111187111101第23页，共73页，2023年，2月20日，星期四8.2.2串行区域技术

10008777111188880000788800017887443378103344881087111111871111011000877711118888000078880001788744337810334488108711111187111101第24页，共73页，2023年，2月20日，星期四8.2图像的分割8.2.1、并行区域技术－基于阈值的图像分割技术8.2.2、串行区域技术8.2.3、并行边界技术－基于边界的图像分割第25页，共73页，2023年，2月20日，星期四8.2.3并行边界技术－基于边界的图像分割边缘检测是所有基于边界的分割方法的第一步。视觉系统对图像的边缘更敏感，而不是根据点的灰度区分出物体，人对边界的识别机理也不是设置一个灰度门限T来分割物体的。图像增强讨论了用梯度、拉普拉斯算子以及高通滤波增强图像边缘轮廓的方法，实际上，对增强后的图像边缘轮廓进行门限化处理，就可以用于边缘检测。第26页，共73页，2023年，2月20日，星期四8.2.3并行边界技术－基于边界的图像分割一、梯度边缘检测

当物体与背景有明显对比度时，物体的边界处于图像梯度最高的点上，通过跟踪图像中具有最高梯度的点的方式获得物体的边界，可以实现图像分割。这种方法容易受到噪声的影响而偏离物体边界，通常需要在跟踪前对梯度图像进行平滑等处理，再采用边界搜索跟踪算法来实现。第27页，共73页，2023年，2月20日，星期四8.2.3并行边界技术－基于边界的图像分割

梯度图像阈值化梯度图像第28页，共73页，2023年，2月20日，星期四8.2.3并行边界技术－基于边界的图像分割二、边界提取与轮廓跟踪

为了获得图像的边缘人们提出了多种边缘检测方法。在边缘图像的基础上，需要通过平滑等处理去除噪声点、毛刺、空洞等不需要的部分，再通过细化、边缘连接和跟踪等方法获得物体的轮廓边界。

例：采用光栅跟踪（顺序跟踪）的方法第29页，共73页，2023年，2月20日，星期四8.2.3并行边界技术－基于边界的图像分割步骤：1）先用高阈值Th（检测阈值）检出图像轨迹曲线的初选点。2）从第一行的初选点开始用较低阈值Tl（跟踪阈值）对初选点的下三个点进行判断。在阈值范围之内的接收，反之去除。3）对应于某个检出点，如果下一行无可接收的对象，则这条曲线跟踪即可结束。第30页，共73页，2023年，2月20日，星期四8.2.3并行边界技术－基于边界的图像分割步骤：4）对应于某个检出点，如果下一行有多个可接收的对象，则这条曲线发生分支。跟踪对各个分支同时进行。5）对于不在第一行的其他初选出来的象素，从该点开始，重新使用跟踪阈值进行跟踪，以检出不是从第一行开始的其他曲线。第31页，共73页，2023年，2月20日，星期四8.2.3并行边界技术－基于边界的图像分割原图像用阈值T＝7处理的结果900080004712533222801721053260098701716021118070326592184107303198787987787987第32页，共73页，2023年，2月20日，星期四8.2.3并行边界技术－基于边界的图像分割用阈值T＝4处理的结果用检测阈值Th＝7和跟踪阈值

Tl＝4处理的结果9875875698776876598479000800007125332228017210532600987017160211180703265921841070131第33页，共73页，2023年，2月20日，星期四8.2.3并行边界技术－基于边界的图像分割检测与跟踪时的准则不一定是灰度，也可以是梯度等反映局部性质的量。光栅顺序跟踪跟光扫描方向有关，最好用其他方向再跟踪一次，例如逆序向上，两种方法综合起来能得到更好的效果。若边缘和扫描方向平行时效果不好，最好在垂直方向跟踪一次。第34页，共73页，2023年，2月20日，星期四8.2.3并行边界技术－基于边界的图像分割三、HOUGH变换

HOUGH变换是利用图像全局特性而将边缘象素连接起来组成区域封闭边界的一种方法。在预先知道区域形状的情况下，用哈夫变换可以很方便的得到边界曲线而将不连续的边缘象素点连接起来。哈夫变换的主要优点：受噪声和曲线间断的影响较小。利用哈夫变换可以检测图像中某些符合参数模型的主导特征，如直线、圆、椭圆等，可以通过对其参数进行聚类的方法，抽取相应的特征。适用于二值图像。第35页，共73页，2023年，2月20日，星期四8.2.3并行边界技术－基于边界的图像分割哈夫变换的基本思想是：点――线对偶性。在图像空间XY中，所有过点的直线满足直线方程：

p为斜率，q为截距所以：可以认为该式是参数空间PQ中过点(p,q)的一条直线。YX0QP这条直线上所有点的p，q相同第36页，共73页，2023年，2月20日，星期四8.2.3并行边界技术－基于边界的图像分割显然，XY空间中的共线的点对应PQ空间中一组相交的线，这组线的交点就是p，q。反过来，在参数空间相交于同一个点的所有直线在图像空间中都有共线的点与之对应。这就是点―线对偶性。哈夫变换就是根据这些关系把图像空间中的检测问题转换到了参数空间中，通过在参数空间中进行简单的累加统计完成检测任务。第37页，共73页，2023年，2月20日，星期四8.2.3并行边界技术－基于边界的图像分割具体实现时，考虑p、q的可能取值范围，从大到小进行累加。

1)初始化A(p,q)=02)对XY空间中的每一点，P取遍所有可能值

3）算出对应的q，给A(p,q)单元加14）根据maxA(p,q)，求出直线第38页，共73页，2023年，2月20日，星期四8.2.3并行边界技术－基于边界的图像分割问题：p、q的取值较大时，计算量增加，可以采用极坐标系形式。直线在极坐标系下具有如下的参数方程形式：这条直线在极坐标下表示为一个点第39页，共73页，2023年，2月20日，星期四8.2.3并行边界技术－基于边界的图像分割

第40页，共73页，2023年，2月20日，星期四8.3图像的特征图像特征是图像的重要属性，它一般应具备4个特点：可区分性、可靠性、独立性、数量少图像的特征主要有：

1)颜色（灰度）特征

2)形状特征

3)纹理特征

4)空间关系特征第41页，共73页，2023年，2月20日，星期四8.4图像的描述边界描述区域或闭合边界的描述第42页，共73页，2023年，2月20日，星期四8.4.1图像的边界描述一、像元之间的连通和邻接（略）二、距离（略）第43页，共73页，2023年，2月20日，星期四8.4.1图像的边界描述三、周长和边长L

周长的定义通常有3种：

1）把区域像元看为正方形，边界只有水平和垂直两种线，相当于像元小方格的四个边的四邻码

2）把像元看做一个点，可用弧长的8邻码计算

3）用边界所占面积第44页，共73页，2023年，2月20日，星期四8.4.1图像的边界描述方法1：L＝25d；d为小正方形的边长方法2：方法3：18@,其中@为小正方形的面积第45页，共73页，2023年，2月20日，星期四8.4.1图像的边界描述四、边界链码链码是从物体边界上的任意一个点（x,y）出发，而后只记录下一个边界点的方向码，直至回到出发点为止的编码序列。在采用边界跟踪方法获取物体边界轮廓时，可以直接得到边界链码。第46页，共73页，2023年，2月20日，星期四8.4.1图像的边界描述8方向链码12345670第47页，共73页，2023年，2月20日，星期四8.4.1图像的边界描述下面的图像用链码表示为：20206644

这个结果与起点有关，为了避免这个问题，用循环码（归一化）表示：

02066442

具体方法为：首尾相连循环形成的自然数最小 22440660第48页，共73页，2023年，2月20日，星期四8.4.1图像的边界描述图像平移链码不变，但旋转时链码会发生变化，为此可用差分码表示原链码(4)20206644原链码(3)17175533差分码66266060（模8、反向）差分码662660602244066011337557第49页，共73页，2023年，2月20日，星期四8.4.1图像的边界描述五、线条的描述

若图中相继相邻的点集能拟合与某些曲线函数，那么这些曲线函数的少量参数就可做为该点集的一种描述。如：常用B样条曲线来拟合，用于从边界识别零件等场合。第50页，共73页，2023年，2月20日，星期四能用少量数据描述物体特征叫做描述子。区域的描述子应满足四点要求：

1）能精确地描述特征

2）对大小变化不敏感

3）对描述的起点不敏感

4）对平移旋转不敏感8.4.1图像的边界描述第51页，共73页，2023年，2月20日，星期四8.4.1图像的边界描述六、傅立叶描述子从任一点开始的一个N点的闭合边界序列（xk,yk）(k=0,…,N)可以认为是一个周期为N的周期函数，把它放到复平面上，就形成了一个1维的复数序列{uk+jvk}（xk,yk）

uk+jvkXUYV第52页，共73页，2023年，2月20日，星期四8.4.1图像的边界描述对这个复数序列

进行傅立叶变换

第53页，共73页，2023年，2月20日，星期四8.4.1图像的边界描述我们只利用s(w)的前M个系数，其余系数置0，反变换就得到s(k)的一个近似傅立叶变换的低频信息对应总体形状而高频信息对应一些细节，这样我们可用对应低频分量的傅立叶系数（个数M远<N）来近似描述边界的形状。（P231图）第54页，共73页，2023年，2月20日，星期四8.4.2区域的描述一、矩描述子

对数字图像f(x,y)，如果它分段连续且在XY平面上只有有限个点不为零，可证明其各阶矩存在，f(x,y)的p＋q阶矩定义为：

第55页，共73页，2023年，2月20日，星期四8.4.2区域的描述可以证明f(x,y)与mpq是一一对应的。区域的矩是用所有区域的点计算出来的，抗干扰性能较好。f(x,y)的p＋q阶中心矩定义为：其中：第56页，共73页，2023年，2月20日，星期四8.4.2区域的描述三阶以下的中心矩：u00、u01＝u10、u11

、u20、u02、u30、u12、u21、u03f(x,y)的归一化中心矩可表示为：其中第57页，共73页，2023年，2月20日，星期四8.4.2区域的描述可以导出具有平移、旋转和尺度变换不变性的7个不变矩：P241图第58页，共73页，2023年，2月20日，星期四8.4.2区域的描述二、区域面积通常面积的计算可以直接对分割出的物体内像元计数即可。从闭合边界来计算面积有类似计算周长的4连通和8连通的问题用不同方法描述的边界计算面积时有很大差异。第59页，共73页，2023年，2月20日，星期四8.4.2区域的描述A＝31（每一小块面积认为是1）第60页，共73页，2023年，2月20日，星期四8.4.2区域的描述第一层有1个白小块(标记为1的)，面积＝4m-1=43-1=16第二层有2个白小块，面积＝2×4m-2=2×43-2=8第三层有11个白小块，面积＝11×4m-3=11×43-3=11总面积＝16＋8＋11＝350000000000100000000001001111111011111110111111101111110011111100用四叉树描述区域时，只需将叶节点的面积相加即可。叶节点面积为4m－l，l为叶节点所在层数，m为图像层数，N＝2m第61页，共73页，2023年，2月20日，星期四8.4.2区域的描述三、区域重心重心是一种全局描述符，区域重心的坐标是根据所有属于区域的点计算出来的。在区域本身尺寸与各区域间的距离相对很小时，可将区域用其重心坐标的质点来代替第62页，共73页，2023年，2月20日，星期四8.4.2区域的描述四、区域的拓扑结构拓扑学是研究图形不受畸变变形影响性质的。欧拉数E＝C（连通组元个数）－H（孔）

E=1-2=-1E=1-1=0BA第63页，共73页，2023年，2月20日，星期四8.4.2区域的描述五、区域的纹理纹理是物体表面结构的模式，如木材、水泥、沙子、玻璃等等，对图像分析和理解十分重要。纹理可以看做许多基本纹理小单元的重复，这些小单元可称作纹理元，纹理元常包含多个象素。常用的纹理描述方法分为统计法和结构法。统计法用于分析木纹、沙地等不规则的纹理；结构法常用于分析印刷图