第七章 图像分割

1、河北大学数学与计算机学院河北大学数学与计算机学院马颍丽马颍丽第七章第七章 图像分割图像分割n7.1 概述概述n7.2 像素的邻域和连通性像素的邻域和连通性n7.3 图像的阈值分割技术图像的阈值分割技术n7.4 图像的边缘检测图像的边缘检测n7.5 霍夫变换霍夫变换n7.6 区域生长法区域生长法n7.7 图像分割方法的比较图像分割方法的比较Ma Yingli CMC HBUn图像分割的定义及其在图像处理中的地位。图像分割的定义及其在图像处理中的地位。n像素间的关系：像素间的关系：n邻域和连通性。邻域和连通性。n阈值分割技术：阈值分割技术：n全局阈值分割和自适应阈值分割技术。全局阈值分割和自适应阈

2、值分割技术。n边缘检测法：边缘检测法：n梯度算子、拉普拉斯算子、拉普拉斯梯度算子、拉普拉斯算子、拉普拉斯- -高斯算子、高斯算子、方向算子、坎尼算子和边缘跟踪。方向算子、坎尼算子和边缘跟踪。n区域检测法。区域检测法。n霍夫变换。霍夫变换。知识要点Ma Yingli CMC HBU7.1 概概 述述u 7.1.1 图像分割的目的和任务图像分割的目的和任务u图像处理的重要任务就是对图像中的对象进行分析和理解。图像处理的重要任务就是对图像中的对象进行分析和理解。u在图像分析中，输出的结果是对图像的描述、分类或其他在图像分析中，输出的结果是对图像的描述、分类或其他的某种结论的某种结论 。u图像分析主要

3、包括以下几部分内容：图像分析主要包括以下几部分内容：（1）把图像分割成不同的区域，或把不同的目标分开（分割）。即把图像分成互不重叠的区域并提取出感兴趣目标。（2）找出各个区域的特征（特征提取）。（3）识别图像中的内容，或对图像进行分类（识别与分类）。（4）给出结论（描述、分类或其他的结论）。 Ma Yingli CMC HBU图图7.1 目标为飞机的图像目标为飞机的图像 （a）原图像 （b）分割后的图像Ma Yingli CMC HBUn令集合令集合R代表整个图像区域，对代表整个图像区域，对R的图像分割可以视为的图像分割可以视为将将R分成分成N个满足以下条件的非空子集：个满足以下条件的非空子集

4、：n（1）n（2）对于所有的对于所有的i和和j，ij ，有，有 ；n（3）对于）对于i = 1, 2, , N，有，有P(Ri) = TRUE； n（4）对于）对于ij ，有，有P(RiRj) = FALSE；n（5）对于）对于i = 1, 2, , N，Ri是连通的区域。是连通的区域。 1;NiiRRijRR 7.1.2 图像分割的集合定义图像分割的集合定义Ma Yingli CMC HBU 根据分割方法的不同，通常有两种分类方法：（1）根据图像的两种特性进行分割：根据各个像素点的灰度不连续性进行分割；根据同一区域具有相似的灰度进行分割。（2）根据分割的处理策略不同进行分割：并行算法，所有的

5、判断和决策可以独立进行；串行算法，后期的处理依赖前期的运算结果。7.1.3 图像分割的分类图像分割的分类Ma Yingli CMC HBU表表7.1 7.1 常见的图像分割算法常见的图像分割算法分类边界（不连续性）区域（相似性）并行处理并行边界类（边缘检测等）并行区域类（阈值分割、聚类等）串行处理串行边界类（边缘跟踪等）串行区域类（区域生长、分裂合并等）Ma Yingli CMC HBU7.2 像素的邻域和连通性像素的邻域和连通性1. 4邻域 n对一个坐标为 的像素p，它可以有两个水平和两个垂直的近邻像素。它们的坐标分别是n这四个像素称为p 的4邻域。n互为4邻域的像素又称为4连通的。2. 8

6、邻域 n取像素p四周的8个点作为相链接的邻域点，除掉p本身外，剩下的8个点就是p的8邻域。n互为8邻域的像素又称为8连通的 。 ),(yx) 1,(),1,(), 1(), 1(yxyxyxyxMa Yingli CMC HBUu目标和背景的连通性定义必须取不同，否则会引起矛盾。0000001110010100011000000图图7.2 目标和背景连通性目标和背景连通性Ma Yingli CMC HBU【例例7.1】根据4/8连通准则在二值图像中判断目标。 解：解：应用函数bwlabel可以根据4连通或8连通准则，在给定的二值图像矩阵BW中寻找目标。MATLAB程序：BW = 1 1 1 0

7、 0 0 0 0; 1 1 1 0 1 1 0 0; 1 1 1 0 1 1 0 0; 1 1 1 0 0 0 1 0; 1 1 1 0 0 0 1 0; 1 1 1 0 0 0 1 0; 1 1 1 0 0 1 1 0; 1 1 1 0 0 0 0 0; %给定的二值图像矩阵L4 = bwlabel(BW,4) %根据4连通准则判定目标L8 = bwlabel(BW,8) %根据8连通准则判定目标Ma Yingli CMC HBU根据4连通准则，得到的目标是3个： L4 = 1 1 1 0 0 0 0 0 1 1 1 0 2 2 0 0 1 1 1 0 2 2 0 0 1 1 1 0 0 0

8、 3 0 1 1 1 0 0 0 3 0 1 1 1 0 0 0 3 0 1 1 1 0 0 3 3 0 1 1 1 0 0 0 0 0根据8连通准则，得到目标是2个： L8 =1 1 1 0 0 0 0 0 1 1 1 0 2 2 0 0 1 1 1 0 2 2 0 0 1 1 1 0 0 0 2 0 1 1 1 0 0 0 2 0 1 1 1 0 0 0 2 0 1 1 1 0 0 2 2 0 1 1 1 0 0 0 0 0Ma Yingli CMC HBU7.3 图像的阈值分割技术图像的阈值分割技术 n灰度阈值分割方法。n若目标和背景具有不同的灰度集合，且两个灰度集合可用一个灰度级阈值T

9、进行分割。n这样就可以用阈值分割灰度级的方法在图像中分割出目标区域与背景区域。n设图像为f (x, y) ，其灰度集范围是Z1, ZK ，在Z1 和ZK 之间选择一个合适的灰度阈值T。7.3.1 基本原理基本原理Ma Yingli CMC HBUn图像分割方法可由下式描述： 这样得到的是一幅二值图像。n 图7.4给出了利用阈值分割图像的实例。n（a）是原图n（b）是对应的直方图n（c）是选择分割阈值为110的结果图。 TyxfTyxfyxg),(0),(1),(（7.1） Ma Yingli CMC HBU（a）原图像）原图像 （b）直方图）直方图 （c）已分割的图像）已分割的图像图图7.4

10、阈值分割阈值分割 Ma Yingli CMC HBUl全局阈值是最简单的图像分割方法。根据不同的目标，选用最佳的阈值。l1 1实验法实验法l需要知道图像的某些特征l2 2直方图法直方图法l适用于目标和背景的灰度差较大，直方图有明显谷底的情况。l3 3最小误差的方法最小误差的方法Ma Yingli CMC HBU7.3.2 全局阈值分割全局阈值分割 图图7.5 7.5 直方图阈值分割示意图直方图阈值分割示意图 图图7.6 7.6 灰度级分布灰度级分布Ma Yingli CMC HBUl当照明不均匀、有突发噪声或者背景灰度变化比较大的时候，可以对图像进行分块处理，对每一块分别选定一个阈值进行分割，

11、这种与坐标相关的阈值称为自适应阈值的方法。l这类算法的时间复杂度和空间复杂度比较大，但是抗噪声的能力比较强 。l任何一种分割方法都有其局限性。l实际的算法只能根据实际情况选择方法和阈值。 Ma Yingli CMC HBU 7.3.3 局部阈值分割局部阈值分割7.4 图像的边缘检测图像的边缘检测 l基于灰度不连续性进行的分割方法。l 图7.7 几种常见的边缘l用差分、梯度、拉普拉斯算子及各种高通滤波处理方法对图像边缘进行增强，只要再进行一次门限化的处理，便可以将边缘增强的方法用于边缘检测。Ma Yingli CMC HBU7.4.1 边缘检测的基本原理边缘检测的基本原理l对于一个连续函数f (

12、x,y)，其在(x,y)处的梯度:l常采用小型模板，然后利用卷积运算来近似，Gx和Gy各自使用一个模板。l1. Roberts算子 yfxfGGfyx（7.2） 100101107.4.2 梯度算子梯度算子 Ma Yingli CMC HBU2. Prewitt算子3. Sobel算子 u通过算子检测后，还需作二值处理从而找到边界点。u图7.5给出了利用这三个算子进行边缘检测的不同效果。u这三种模板中，Sobel算子的检测效果最好。 101101101111000111101202101121000121Ma Yingli CMC HBU(a) 原图像原图像 (b) Roberts算子检测算子

13、检测 (c) Prewitt算子检测算子检测 (d) Sobel算子检测算子检测 【例例7.3】利用梯度算子对图像进行边缘检测。利用梯度算子对图像进行边缘检测。Ma Yingli CMC HBUn I = imread(blood1.tif);nimshow(I);nBW1 = edge(I,roberts); %进行Roberts算子边缘检测，门限值采用默认值nBW2 = edge(I,prewitt); %进行Prewitt算子边缘检测，门限值采用默认值nBW3 = edge(I,sobel); %进行Sobel算子边缘检测，门限值采用默认值nfigure,imshow(BW1,);nfi

14、gure,imshow(BW2,);nfigure,imshow(BW3,); Ma Yingli CMC HBU Laplacian是二阶导数算子，也是借助模板来实现的。l对模板有一些基本要求：l模板中心的系数为正，其余相邻系数为负，且所有的系数之和为零。l常用的模板有：010141010010151010Ma Yingli CMC HBU7.4.3 拉普拉斯算子拉普拉斯算子图图7.9 二阶导数算子确定边缘位置二阶导数算子确定边缘位置（a a）平滑边缘的二阶算子）平滑边缘的二阶算子 （b b）斜坡边缘的二阶算子）斜坡边缘的二阶算子Ma Yingli CMC HBU【例例7.4】Robert、

15、Sobel和和Laplace算子的边缘检测。算子的边缘检测。 （a）Lena图像图像 （b）Robert算子检测结果算子检测结果 （c）Sobel算子检测结果算子检测结果 （d）Laplace算子检测结果算子检测结果图图7.10 各种算子的检测结果各种算子的检测结果Ma Yingli CMC HBU7.4.4 拉普拉斯拉普拉斯-高斯算子高斯算子n拉普拉斯-高斯算子（Laplacian of Gaussian, LoG），也称Marr算子。n思想：n先用高斯函数g(x, y)对图像f (x, y)进行滤波，n再对滤波后的图像进行拉普拉斯运算（2），结果为零的位置即为边缘点的位置。n滤波提高了抗噪

16、声的能力，但同时可能使原本比较尖锐的边缘平滑了，甚至无法检测到。Ma Yingli CMC HBU图图7.8 不同微分算子的边缘检测效果不同微分算子的边缘检测效果（a）原图像 （e）LoG算子检测 Ma Yingli CMC HBUCanny的主要工作：l推导了最优边缘检测算子。l考核边缘检测算子的指标是：l低误判率，即尽可能少地把边缘点误认为是非边缘点；l高定位精度，即准确地把边缘点定位在灰度变化最大的像素上；l抑制虚假边缘。 Ma Yingli CMC HBU7.4.5 Canny边缘检测算子边缘检测算子图图7.8 不同微分算子的边缘检测效果不同微分算子的边缘检测效果（a）原图像 （e）C

17、anny算子检测 Ma Yingli CMC HBUl判断一个像素是否为边缘点的条件为：（1）像素(i, j)的边缘强度大于沿梯度方向的两个相邻像素的边缘强度；（2）与该像素梯度方向上相邻两点的方向差小于45；（3）以该像素为中心的33邻域中的边缘强度的极大值小于某个阈值 。Ma Yingli CMC HBUl Canny算子的检测比较优越，可以减少小模板检测中边缘中断，有利于得到较完整的边缘。lMATLAB程序：nI = imread(blood1.tif);nimshow(I);nBW5 = edge(I, canny);nfigure,imshow(BW5,); Canny算子边缘检测的

18、结果图算子边缘检测的结果图 Ma Yingli CMC HBUl利用一组模板对图像中的同一像素求卷积，然后选取其中最大的值作为边缘强度，而将与之对应的方向作为边缘方向。l相对于梯度算子的优点：不仅仅只考虑水平和垂直方向，还可以检测其他方向上的边缘。l但计算量将大大增加。l常用的有8方向Kirsch（33）模板，如图7.7所示，方向间的夹角为45。 Ma Yingli CMC HBU7.4.6 方向算子方向算子图图7.11 33 Kirsch算子的八方向模板算子的八方向模板 Ma Yingli CMC HBU 上述方法仅得到处在边缘上的像素点。l噪声和不均匀的照明而产生的边缘间断的影响，使得经过

19、边缘检测后得到的边缘像素点很少能完整地描绘实际的一条边缘。l可以在使用边缘检测算法后，接着使用连接方法将边缘像素组合成有意义的边缘。Ma Yingli CMC HBU 7.4.7 边缘跟踪边缘跟踪l光栅扫描跟踪法：l一种简单的利用局部信息、通过扫描的方式将边缘点连接起来的方法。l该跟踪算法采用电视光栅行扫描顺序对遇到的像素进行分析，从而确定其是否为边缘。l由于光栅扫描跟踪和扫描方向有关，因此最好沿其他方向再跟踪一次 。 Ma Yingli CMC HBU图图7.12 光栅扫描跟踪光栅扫描跟踪（a）输入图像）输入图像 （b）阈值化处理）阈值化处理 （c）根据阈值进行跟踪）根据阈值进行跟踪Ma Y

20、ingli CMC HBU【例例7.5】利用函数对图像进行边缘跟踪。利用函数对图像进行边缘跟踪。nBW = imread(blobs.png); nimshow(BW,);ns = size(BW);nfor row = 2:55:s(1)nfor col=1:s(2) if BW(row,col), break; endnendncontour = bwtraceboundary(BW, row, col, W, 8, 50,counterclockwise);nif(isempty(contour) hold on;n plot(contour(:,2),contour(:,1),g,Li

21、neWidth,2);n hold on; plot(col, row,gx,LineWidth,2); nelsenhold on; plot(col, row,rx,LineWidth,2);nend Ma Yingli CMC HBU图图7.14 二值边缘跟踪二值边缘跟踪Ma Yingli CMC HBU 7.5 霍夫变换霍夫变换l霍夫（Hough）变换方法是利用图像全局特性而直接检测目标轮廓，将图像的边缘像素连接起来的常用方法。l1.基本原理l点线的对偶性。l当给定图像空间的一些边缘点，就可以通过霍夫变换确定连接这些点的直线方程。 Ma Yingli CMC HBU2霍夫变换的实现l实

22、际进行霍夫变换时，要在上述基本方法的基础上根据图像的具体情况采用一些措施。l使用极坐标直线方程用以提高精度和速度。 Ma Yingli CMC HBU7.5.1 直角坐标系中的霍夫变换直角坐标系中的霍夫变换图图7.15 7.15 图像空间和参数空间中点和线的对偶性图像空间和参数空间中点和线的对偶性 Ma Yingli CMC HBU图图7.16 参数空间中的累加数组参数空间中的累加数组 Ma Yingli CMC HBU7.5.2 极坐标系中的霍夫变换极坐标系中的霍夫变换 图7.17 直线的极坐标表示 图7.18 参数空间对应的曲线 Ma Yingli CMC HBU图图7.19 采用霍夫变换

23、提取图像中的直线采用霍夫变换提取图像中的直线 （a）原图像 （d）检测到的直线Ma Yingli CMC HBU 7.6 区域生长法区域生长法l将具有相似性质的像素集合起来构成区域。l先对每个需要分割的区域找一个种子像素作为生长的起点l然后将种子像素周围邻域中与种子像素具有相同或相似性质的像素合并到这一区域中。l将这些新像素当做新的种子像素继续进行上面的过程，直到再没有满足条件的像素可被包括进来。这样一个区域就长成了。Ma Yingli CMC HBU7.6.1 原理和步骤原理和步骤在实际应用区域生长法时需要解决三个问题： 选择一组能正确代表所需区域的种子像素；n种子像素的选取常可借助具体问题

24、的特点进行。 确定在生长过程中将相邻像素包括进来的准则；n生长准则的选取不仅依赖于具体问题本身，也和所用图像数据的种类有关 ， 制定让生长过程停止的条件或规则。n一般生长过程在进行到再没有满足生长准则需要的像素时停止。 Ma Yingli CMC HBU465225650255077401574015 51 155511555115551155511555115552255512555117751157511(a)原图像原图像 (b)T=3的生长结果的生长结果 (c) T=1的生长结果的生长结果图图7.20 区域生长区域生长Ma Yingli CMC HBUl区域生长的一个关键：l选择适合的生长准则，大部分区域生长准则使用图像的局部性质。l生长准则可根据不同的原则制定，而使用不同的生长准则会影响区域生长的过程。Ma Yingli CMC HBU7.6.2 生长准则和过程生长准则和过程1. 基于区域灰度差区域生长方法将图像以像素为基本单位来进行操作（1）对图像进行逐行扫描