第7章 图像分割与特征提取及MATLAB实现
图像分割：是指将图像中有意义的对象与其背景分离，并把这些对象按照不同的含义分割开来，也就是说，把图像中具有不同含义的对象提取出来。
图像分割的方法：大致可以分为基于边缘检测的方法和基于区域生成的方法两大类。 边缘检测技术：是所有基于边界分割的图像分析方法的第一步，首先检测出图像局部特性的不连续性，再将它们连成边界，这些边界把图像分成不同的区域，检测出边缘的图像就可以进行特征提取和形状分析。
图像特征：是指图像的原始特性或属性。其中有些是视觉直接感受到的自然特征，如区域的亮度、边缘的轮廓、纹理或色彩等；有些是需要通过变换或测量才能得到的人为特征，如变换频谱、直方图、矩等。图像特征提取工作的结果给出了某一具体的图像中与其他图像相区别的特征。如：描述物体表面灰度变化的纹理特征，描述物体外形的形状特征等。这些特征提取的结果需要一定的表达方式，要让计算机能懂得，这就是本章的任务。
7.1边缘检测方法
图像边缘对图像识别和计算机分析十分有用，边缘能勾画出目标物体，使观察者一目了然；边缘蕴含了丰富的内在信息（如方向、阶跃性质、形状等），是图像识别中重要的图像特征之一。从本质上说，图像边缘是图像局部特性不连续性（灰度突变、颜色突变、纹理结构突变等）的反映，它标志着一个区域的终结和另一个区域的开始。为了计算方便起见，通常选择一阶和二阶导数来检测边界，利用求导方法可以很方便地检测到灰度值的不连续效果。边缘的检测可以借助空域微分算子利用卷积来实现。常用的微分算子有梯度算子和拉普拉斯算子等，这些算子不但可以检测图像的二维边缘，还可检测图像序列的三维边缘。下面分别进行介绍。
7.1.1边缘算子法
（1）差分算子
图7.1是最常见的边缘，其灰度变化可能呈阶梯状，也可能呈脉冲状，对于图7. 1(a)和(d)所示的一阶差分和二阶差分如图7.1(b)、(c)和(e)、(f)。由图可知，差分图像中能够较精确地获得这两种类型的边缘。
边缘与差分值的关系归纳如下：边缘发生在差分最大值（如图7. 1(b)）或最小值处（如图(f)）；边缘发生在过零点处（如图(c)、(e)）。
图7.1常见边缘的一阶差分和二阶差分
(a)阶梯状边缘 （b）阶梯状边缘的一阶差分（c）阶梯状边缘的二阶差分 (d)脉冲状边缘 （e）脉冲状边缘的一阶差分（f）脉冲状边缘的二阶差分
1）梯度算子
在前面章节曾经讨论，在点f(m,n)处，梯度G[f(m,n)]的幅度为：
12???f???f??G?????? （7.1.1） ?f?m,n????????n?????m???22对于数字图像，式(7.1.1)可改写为：
1G??f?m,n???????mf其中：?mf?f?m,n??f?m?1,n?222??nf?? （7.1.2）
?nf?f?m,n??f?m,n?1? （7.1.3）
有时为了避免平方和运算，可将幅度用两个分量的绝对值之和或最大绝对值来表示，即
G??f?m,n?????mf??nf
或 G??f?m,n????ma?x?mf?,nf? （7.1.4）
取适当的门限T，如果G[f(m,n) ] > T，则（m,n）为阶跃状边缘点。
2）罗伯特（Robert）梯度
前面，在计算（m,n）点的梯度时只用到f(m,n), f(m-1,n)及f(m, n-1)的值。但实际上，任意一对相互垂直方向上的差分都可用来估计梯度。Robert梯度采用对角方向相邻两像素之差，即
?mf?f?m,n??f?m?1,n?1??nf?f?m?1,n??f?m,n?1? （7.1.5）
??12,n?1??2?有了?mf,?nf，很容易地算出Robert梯度的幅值。Robert梯度实际上是以?m?为中心的，应当把它们看成在这个中心点上连续梯度的近似。从图像处理的实际效果看，用式(7.1.5)的Robert梯度比用式(7.1.3)的梯度计算式来检测边缘要好。
3）拉普拉斯（Laplacian）算子
拉普拉斯（Laplacian）算子是一种二阶微分算子，也可用来提取图像的边缘。在数字图像处理中，其一般表示形式为：
?f2?m,n???u,v?s??f?u,v??f?m,n??? （7.1.6）
其中S是以f(m,n)为中心的邻点的集合，可以是上、下、左、右4邻点或8邻点的集合，或
者是对角线4邻点的集合（如图7.2所示），与其相对应的表达式分别为：
?f?m,n??f?m?1,n??f?m?1,n??f?m,n?1??f?m,n?1??4f?m,n? （7.1.7）
2?f?m,n??f?m?1,n?1??f?m?1,n??f?m?1,n?1??f?m,n?1??2f?m,n?1??f?m?1,n?1??f?m?1,n??f?m?1,n?1??8f?m,n?2 （7.1.8）
?f?m,n??f?m?1,n?1??f?m?1,n?1??f?m?1,n?1??f?m?1,n?1??4f?m,n? （7.1.9）
图7.2 Laplacian算子集合
下面以图7.3所示的图像边缘灰度分布，运用式（7.1.8）进行计算，以了解Laplacian算子用于边缘提取的特性。处理结果如图7.4。在图7.3和图7.4中，中心的点用作比较输人和输出图案的参考值，线框为处理范围。由图可知，
孤立点（图7.3(a)）的输出是一个略为扩大或略带模糊的点，其输出幅度是该点灰度值的4倍（图7.4(a)）；
对于线结构（图7.3(b)），输出宽度加粗，外观仍呈线型，组成初始线的各位置上的?2f值是原来幅度的2倍（图7.4(b)）；
对阶跃边缘（图7.3(c)），其检测结果是线（图7.4(c)）。
(a) (b) (c)
图7.3灰度分布
图7. 4 ?2f
以上是直接以?2f?m,n?作为边缘像素灰度，实际工作中也可以把?f?m,n??THg2或f?m,n???f?m,n?的像素作为边缘元。不同的定义方法其边缘检测出所获得的目标图
2像，在使用中可通过试验选择与实际图像相匹配的算法。
(2) Laplacian-Gauss算子
梯度算子和Laplacian算子对噪声比较敏感。对此，一方面可在运用这两种算子作边缘提取前，先用邻域平均法等作平滑处理，另一方面可先用高斯形二维低通滤波器对图像f(m,n)进行滤波，然后再对图像作Laplacian边缘提取，这种方法被称为Laplacian-Gauss算子法。具体是令g(m,n)为高斯低通滤波后的图像，?2g?m,n?表示边缘提取后图像，则有
?g?m,n??2?f?m22??f?n22?m2?n2??m2?n2???2?exp??? （7.1.10） ?2222?2???m?n2?????exp???22???g?m,n?式中，??W22，参数W为正瓣宽度（如图7.5）。
图7.5 ?g?m,n?算子的截面图
2
图7.5描述了函数?2g?m,n?的截面情况。它以坐标原点中心全方位对称，并带有与主瓣相反的拖尾。也可用Laplacian-Gauss算子对f(m,n)作卷积计算，所用窗口不宜太小，窗口尺寸与主瓣宽度有关，窗口模板内的系数之和为零。边缘与卷积计算 后的零交叉点对应。
7.1.2模板匹配法
模板是为了检测某些区域特征而设计的阵列。设有一个3×3模板窗口W，其元素Wi,j的位置如图7. 6(a)所示，一幅图像F的各元素f(m,n)的位置如图7.6(b)所示。模板匹配的过程是求乘积和的过程
lli,jg?m,n????f?m?i,n?j?Wi??lj??l （7.1.11）
其中，以g(m,n)为边缘检测模板输出，l = L/2 , L为窗口宽度，对于3×3模板窗口l=-1,0,1。在上述3×3模板例子中，可以分别把模板阵列W和局部图像堆叠成9维向量W=(W-1,-1, W0,-1,?, W1,1,)T以及F=(f(m-1,n-1), f(m,n-1),?, f(m+1,n+1))T，则式（7.1.1）的运算相当于计算两个向量W和F的内积FTW，即为
?F,W??FW?FWcos? （7.1.12）
T其中θ为两处向量之间的夹角，若|W| =1，内积就等于投影。因此，可以将式(7.1.11)的运算称作计算F在W上投影。下面分别讨论点、线条和边缘的模板。
图7.6 模板W及图像F的各元素
（1）点模板
图7.7所示为点模板，通常用于背景强度恒定、目标图像 灰度相同或基本相同的图像（如二值图像）。在检测区域时， 点模板通常是拖动模板在图像域移动，横向移动间隔取1个像
素，纵向移动间隔取1个扫描行。在每一个位置上，将模板元 图7.7匹配滤波器点模板 素分别与相应的图像灰度级相乘并求和（当和小于零时，可作 取绝对值和零对待两种处理）。
当g(m,n)＝0时，说明模板位于均衡背景或目标区域内部； 当g(m,n)≠0时，表明当前窗口内既有背景也有目标；
当g(m,n) =max时，表明模板中心正好是目标和背景的交界处； 当g(m,n)值减少时，表明模板中心离开目标和背景的交界处。 图7.8所示的结果正好与之结论相吻合。其中，
图7. 8(a)是背景灰度为1、目标灰度为4的二值图像；图7. 8(b)是点模板中间计算结果；图7.8(c)是点模板最后匹配输出。
(2)线模板
图7.9所示为线模板，它能够有效地检出线型类型图像结构。其中W1、W2、W3、W4分别为00，450，900，-450四个方向的线模板。当W1、W2、W3、W4在图像域内自上而下、自左向右移动时，在背景灰度级不变、线宽度为1个像素的情况下，W1、W2、W3、W4模板分别对水平线、450方向线、垂直线与斜线（＼）响应最佳。