姚敏图像处理课件第九章.ppt

1,Digital Image Processing,数字图像处理, E-MAIL:myao@,姚 敏,2,第九章 图像分割,3,9.1 概述,4,图像分割,图像,图像识别,图像 预处理,图像理解,,图9.1图像分割在整个图像处理过程中的作用,,图像分割,,,作用,5, 分割出来的各区域对某种性质例如灰度，纹理而言具有相似性，区域内部是连通的的且没有过多小孔  区域边界是明确的  相邻区域对分割所依据的性质有明显的差异,图像分割,特征,6, 边缘检测  边缘跟踪  阈值分割  区域分割  运动分割,图像分割,本章要点,7,9.2 边缘检测,8,梯度算子,梯度算子是一阶导数算子,幅值,方向角,9,梯度算子, 数字图像处理中用差分代替微分,,,,近似计算,10,梯度算子,Roberts算子,11,梯度算子,Sobel算子,12,梯度算子,Prewitt算子,13,梯度算子,原图,Prewitt算子,Sobel算子,Roberts算子,14,拉普拉斯算子,差分,微分,二阶导数算子,15,拉普拉斯算子,,图9.5 两种常用的拉普拉斯算子模板,16,拉普拉斯算子,,,,17,Canny算子,,,好的检测结果：对边缘的错误检测率要尽可能低，在检测出图像真实的边缘的同时要避免检测出现虚假的边缘。

 好的边缘定位精度：标记出的边缘位置要和图像上真正边缘的位置尽量接近 对同一边缘要有低的响应次数：有的算子会对一个边缘回产生多个响应也就是说图像上本来只有一个边缘点的，可是检测出来就会出现多个边缘点 克服噪声的影响,基本思想,18,Canny算子,,,算法步骤,用高斯滤波器平滑图像 计算滤波后图像梯度的幅值和方向 对梯度幅值应用非极大值抑制，其过程为找处图像梯度中的局部极大值点，把其它非局部极大值点置零以得到得到细化的边缘 用双阈值算法检测和连接边缘，使用两个阈值T1和T2(T1>T2)，T1用来找到每条线段，T2用来在这些线段的两个方向上延伸寻找边缘的断裂处，并连接这些边缘19,Canny算子,,,实例,,,20,算子比较,,,,, Roberts算子：Roberts算子利用局部差分算子寻找边缘，边缘定位精度较高，但容易丢失一部分边缘，同时由于图像没经过平滑处理，因此不具备能抑制噪声能力该算子对具有陡峭边缘且含噪声少的图像效果较好  Sobel算子和Prewitt算子：都是对图像先做加权平滑处理，然后再做微分运算，所不同的是平滑部分的权值有些差异，因此对噪声具有一定的抑制能力，但不能完全排除检测结果中出现的虚假边缘。

虽然这两个算子边缘定位效果不错，但检测出的边缘容易出现多像素宽度21,算子比较,,,,, Laplacian算子：是不依赖于边缘方向的二阶微分算子算子，对图像中的阶跃型边缘点定位准确，该算子对噪声非常敏感，它使噪声成分得到加强，这两个特性使得该算子容易丢失一部分边缘的方向信息，造成一些不连续的检测边缘，同时抗噪声能力比较差22,算子比较,,,,, LOG算子：该算子首先用高斯函数对图像作平滑滤波处理，然后才使用Laplacian算子检测边缘，因此克服了Laplacian算子抗噪声能力比较差的缺点，但是在抑制噪声的同时也可能将原有的比较尖锐的边缘也平滑掉了，造成这些尖锐边缘无法检被测到应用LOG算子，高斯函数中方差参数的选择很关键，对图像边缘检测效果有很大的影响高斯滤波器为低通滤波器，越大，通频带越窄，对较高频率的噪声的抑制作用越大，避免了虚假边缘的检出，同时信号的边缘也被平滑了，造成某些边缘点的丢失反之，越小，通频带越宽，可以检测到的图像更高频率的细节，但对噪声的抑制能力相对下降，容易出现虚假边缘因此，应用LOG算子，为取得更佳的效果，对于不同图像应选择不同参数23,算子比较,,,,, Canny算子：Canny算子虽然是基于最优化思想推导出的边缘检测算子，实际效果并不一定最优，原因在于理论和实际有许多不一致的地方。

该算子同样采用高斯函数对图像作平滑处理，因此具有较强的抑制噪声能力，同样该算子也会将一些高频边缘平滑掉，造成边缘丢失Canny算子其后所采用用双阈值算法检测和连接边缘，采用的多尺度检测和方向性搜索较LOG算子要好24,9.3 边缘跟踪,25,基本步骤,从图像中一个边缘点出发，然后根据某种判别准则搜索下一个边缘点以此跟踪出目标边界确定边界的起始搜索点，起始点的选择很关键，对某些图像，选择不同的起始点会导致不同的结果确定合适边界判别准则和搜索准则，判别准则用于判断一个点是不是边界点，搜索准则则指导如何搜索下一个边缘点 确定搜索的终止条件26,二值图像边界跟踪,,,27,灰度图像边界跟踪,,,28,Hough 变换,,,Hough变换可以用于将边缘像素连接起来得到边界曲线 优点在于受噪声和曲线间断的影响较小 在已知曲线形状的条件下，Hough变换实际上是利用分散的边缘点进行曲线逼近，它也可看成是一种聚类分析技术,29,Hough 变换,,,,原始图像,二值化图像,细化图像,Hough 变换检测出的直线,直线检测,30,9.4 阈值分割,31,基本原理,原始图像——f(x,y) 灰度阈值——T 阈值运算得二值图像——g(x,y), 阈值选择直接影响分割效果，通常可以通过对灰度直方图 的分析来确定它的值。

32,阈值选择,利用灰度直方图求双峰或多峰 选择两峰之间的谷底作为阈值,33,人工阈值,人工选择法是通过人眼的观察，应用人对图像的知识，在分析图像直方图的基础上，人工选出合适的阈值也可以在人工选出阈值后，根据分割效果，不断的交互操作，从而选择出最佳的阈值34,,,T=155的二值化图像 T=210的二值化图像,原始图像,图像直方图,人工阈值,35,自动阈值,迭代法,,基本思想：开始时选择一个阈值作为初始估计值，然后按某种策略不断地改进这一估计值，直到满足给定的准则为止在迭代过程中，关键之处在于选择什么样的阈值改进策略，好的阈值的改进策略应该具备两个特征，一是能够快速收敛，二是在每一个迭代过程中，新产生阈值优于上一次的阈值36,自动阈值,迭代法,,(1) 选择图像灰度的中值作为初始阈值Ti=T0 (2) 利用阈值Ti把图像分割成两部分区域， R1和R2，并计算其灰度均值,(3) 计算新的阈值Ti+1,(4) 重复步骤2、3，直到Ti+1和Ti的值差别小于某个给定值,37,自动阈值,迭代法,,原始图像,迭代阈值二值化图象,图9.20 迭代式阈值二值化图像,38,自动阈值,Ostu法,, Ostu方法是使类间方差最大而推导出的一种能自动确定阈值的方法具有简单，处理速度快的特点 Matlab工具箱提供的graythresh函数求取阈值采用的正是Ostu法,39,自动阈值,Ostu法,,原始图像,Otsu法二值化图象,图9.21 Otsu法二值化图像,40,自动阈值,分水岭算法,,分水岭算法(watershed)是一种借鉴了形态学理论的分割方法，它将一幅图象看成为一个拓扑地形图，其中灰度值被认为是地形高度值。

高灰度值对应着山峰，低灰度值处对应着山谷将水从任一处流下，它会朝地势底的地方流动，直到某一局部低洼处才停下来，这个低洼处被称为吸水盆地，最终所有的水会分聚在不同的吸水盆地, 吸水盆地之间的山脊被称为分水岭，水从分水岭流下时，它朝不同的吸水盆地流去的可能性是相等的 将这种想法应用于图像分割，就是要在灰度图像中找出不同的吸水盆地和分水岭，由这些不同的吸引盆地和分水岭组成的区域即为我们要分割的目标41,自动阈值,分水岭算法,,(a)原始图像 (b)图像对应的拓扑地形图图9.23 图像对应的拓扑表面图,,42,自动阈值,分水岭算法,,分水岭阈值选择算法可以看成是一种自适应的多阈值分割算法,图9.24 分水岭形成示意图,分水岭对应于原始图像中的边缘,,,,43,自动阈值,分水岭算法,,L=watershed(f),MATLAB函数,图9.25 不准确标记分水岭算法导致过分割,,原始图像,分水岭,分割结果,局部极小值,44,自动阈值,改进的分水岭算法,,图9.26 准确标记的分水岭算法分割过程,,原始图像,原图像的距离变换,标记外部约束,标记内部约束,由标记内外部约束重构的梯度图,分割结果,45,9.5 区域分割,46,基本思想,阈值分割法由于没有或很少考虑空间关系，使多阈值选择受到限制 于区域的分割方法可以弥补这点不足，它利用的是图像的空间性质，该方法认为分割出来的属于同一区域的像素应具有相似的性质，其概念是相当直观的。

传统的区域分割算法有区域增长法和区域分裂合并法该类方法在没有先验知识可以利用时，对含有复杂场景或自然景物等先验知识不足的图像进行分割, 也可以取得较好的性能但是，空间和时间开销都比较大47,区域生长,区域生长法主要考虑象素及其空间邻域象素之间的关系 开始时确定一个或多个象素点作为种子，然后按某种相似性准则增长区域，逐步生成具有某种均匀性的空间区域，将相邻的具有相似性质的象素或区域归并从而逐步增长区域，直至没有可以归并的点或其它小区域为止 区域内象素的相似性度量可以包括平均灰度值、纹理、颜色等信息48,区域生长,选择合适的种子点 确定相似性准则(生长准则) 确定生长停止条件,步骤,49,区域生长,实例,图9.29 区域生长,,原始图像及种子点位置,三个种子点区域生长结果,原始图像及种子点位置,四个种子点区域生长结果,50,区域分裂,,条件：如果区域的某些特性不满足一致性准则 开始：从图像的最大区域开始，一般情况下，是从整幅图像开始 注意： 确定分裂准则(一致性准则)确定分裂方法，即如何分裂区域，使得分裂后的子区域的特性尽可能都满足一致性准则值51,区域分裂,算法, 形成初始区域  对图像的每一个区域Ri，计算P(Ri)，如果 P(Ri)=FALSE则沿着某一合适的边界分裂区域  重复步骤2，当没有区域需分裂时，算法结束。

52,区域合并,,合并运算就是把相邻的具有相似性质的区域合成为一个区域 合并算法中最重要的运算是确定两个区域的相似性 评判区域相似性方法有许多，评判相似性的方法可以基于区域的灰度值，也可以基于区域边界的强弱性等因素一种简单的方法是比较它们的灰度均值,基本思想,53,区域合并,,使用某种方法进行图像的初始区域分割 对于图像中乡邻个区域，计算是否满足一致性谓词词,若满足则合并为一个区域 重复步骤2，直到没有区域可以合并，算法结束算 法,54,区域分裂合并,,数据结构,55,区域分裂合并,,数据结构,图 9.31图像分裂示例,56,区域分裂合并,,算法,设整幅图像为初始区域 对每一区域R，如果P(R)=FLASE，则把该区域分裂成四个子区域 重复上一步，直到没有区域可以分裂 对图像中任意两个相邻的R1和R2，如果P(R1UR2)=TRUE，则把这两个区域合并成一个区域 重复上一步，直到没有相邻区域可以合并，算法结束,57,区域分裂合并,,,,实例,原始图像,分裂最小子区域大小1616,分裂最小子区域大小88,分裂最小子区域大小88,图9.32图像分裂合并,。