数字图像处理第6章

1、数数 字字 图图 像像 处处 理理第六章第六章 图像分割图像分割图像分割门限法边缘检测区域生长分裂与合并分水岭方法图像分割一般的图像处理过程图像输入光电转换数字化图像增强图像恢复图像编码预处理阈值分割边缘检测区域分割图像分割图像识别图像分析理解描述解释特征提取图像分割v在对图像的研究和应用中，往往对图像中的部分区域感兴趣，这些区域被称为前景或目标具有独特性质的区域。v定义：把图像划分为不相交的有意义的区域，同一区域内的像素反映了目标的一致性属性（同一性）。v在无任何先验知识的情况下，一致性的依据和准则取决于图像的灰度和颜色。v图像分割（Image segmentation）是图像分析中的关键技

2、术由于目标特征不均匀，场景复杂，图像中有噪声等原因分割困难。v图像分割方法一般是基于像素的两个特性之一：不连续性和相似性。不连续性：指不同区域之间像素值差别比较大，在边界上具有某种不连续性基于不连续性的分割方法相似性：指区域内部象像素值之间具备一定的相似性基于区域的分割方法图像分割门限法一、基于不连续性的分割方法1 1、门限法（阈值分割）、门限法（阈值分割）( , )( , )( , )EBLf x yTg x yLf x yT阈值T 的选择 基于直方图单目标和背景，单一门限多个目标和背景，多门限图像分割门限法图像分割门限法图像分割门限法自动计算阈值T：（1）选择一个T的初始估计值。（2） 用

3、T分割图像，G1由所有灰度值大于T的像素组成，G2由所有小于或等于T的像素组成。（3）对区域G1和G2中的所有像素计算平均灰度值1和2 。（4）计算新的门限值： （5）重复步骤2到4，直到逐次迭代得到的T值之差小于某个小量（收敛）。动态门限（自适应门限） 由于光照不均匀，导致图像整体有明暗分布（但目标与背景的相对灰度有差异）将图像划分成子区域（或用移动窗）对子区域进行局部门限分割。2121T图像分割门限法(a)原图；(b)全局门限分割结果；(c) 图像划分成子区域；(d)自适应门限分割结果。abcd门限法的缺点图像分割门限法最优门限的选取 zpPzpPzp2211 多数情况下，目标和背景的灰度

4、分布有重叠。若二者的灰度分布的概率密度函数已知，则可以选择门限使得错误概率最小（统计最优）。图像整体灰度级变化的总概率密度函数：图像中两个区域的灰度级概率密度函数（P1和P2是两类象素出现的概率）目标背景图像分割门限法 dzzpTET21 TEPTEPTE2112 TpPTpP2211 dzzpTET12将一个背景点当作目标点进行分类时，错误概率为：将一个目标点当作背景点进行分类时，错误概率为:出错率的整体概率是：对E(T)求导并令导数为0，得解出的T即为最佳门限。如果P1=P2，则最佳门限位于P1(z)和P2(z)的交点处。(6-1-1-1)图像分割门限法从T的表达式知，为了求取T，需要知道

5、两个概率密度。在现实中并不是总可以对这两个密度进行估计。通常的做法是利用参数化模型。例如常考虑使用高斯密度： 2222212122221122zzePePzp02CBTAT将该方程用于(6-1-1-1)得下列门限T的解：2221A2122212B2112222121222221ln2PPC其中图像分割边缘检测2 2、边缘检测、边缘检测 两个具有不同灰度值的相邻区域之间总存在灰度边缘，而灰度边缘是灰度值不连续（或突变）的结果。常用一阶和二阶导数来检测边缘（通过设置阈值将边缘像素与非边缘像素二值化）。边缘模型：阶跃边缘渐变边缘图像分割边缘检测一阶导数可以用于检测图像中的一个点是否是边缘点；二阶导数

6、的符号可以判断一个边缘像素的位置。图像分割边缘检测图中第一列的图像显示了4个斜坡边缘的特写，分别被均值为0且= 0,0.1,1,10 的随机高斯噪声污染。第二列是一阶导数图像和灰度级剖面线。第三列为二阶导数图像和灰度级剖面线。噪声的影响噪声的影响图像分割边缘检测边缘检测算子： Prewitt算子算子 Roberts算子算子 Sobel算子算子图像分割边缘检测用于检测对角线方向上突变的用于检测对角线方向上突变的Prewitt和和Sobel模板模板图像分割边缘检测(a)原图，(b) x方向上的梯度分量|Gx| ，(c) y方向上的梯度分量|Gy|，(d)梯度图像|Gx|+|Gy|abcd y x图

7、像分割边缘检测（采用55的均值滤波器预先进行了平滑处理的结果）图像分割边缘检测对角方向边缘检测。(a)用Sobel +45模板。(b)用Sobel -45模板。从图中可以看到两个模板对水平和垂直的边缘具有相似的响应，但要比在这两个方向上的水平和垂直模板的响应要弱。图像分割边缘检测 Laplacian算子算子一阶导数的局部最大值对应二阶导数的零交叉点（Zero crossing）。所以通过求图像的二阶导数的零交叉零交叉就能找到精确边缘点。拉氏算子对噪声敏感，因此在实际中通常不直接使用。一般结合平滑算子使用。图像分割边缘检测 Marr-Hildreth算子算子 一般情况下，需将平滑滤波与边缘检测结

8、合。先做滤波: g(x,y) = f (x,y)h(x,y)再对g(x)求一阶或二阶导数以检测边缘点。利用卷积性质，有 g(x,y) = f (x,y) h(x,y); g(x,y) = f (x,y) h(x,y);Marr边缘检测算法l平滑滤波器采用高斯滤波器；l边缘检测判据是二阶导数零交叉点；l采用线性插值的方法估计边缘的位置。l因为采用Laplacian算子，故称也为LoG（Laplacian of Gaussian）滤波器。图像分割边缘检测222221,2xyh x ye2222222222242,112xyhhLoGh x yxyxye 55Log模板图像分割边缘检测(a) 血管造

9、影图像。(b) 图像的Sobel梯度。(c)用于得到2727的空间平滑处理模板的高斯函数。(d)Laplacian模板。(e)通过对原图使用高斯型平滑模板进行平滑处理，然后使用拉普拉斯算子模板后得到的LoG图像。(f)对LoG设置门限后的结果。(g)零交点。图像分割边缘检测 Canny算子算子 图像边缘检测必须满足两个条件：能有效地抑制噪声；尽量精确确定边缘的位置。根据对信噪比与定位乘积进行测度，得到最优化逼近算子Canny边缘检测算子。类似于Marr边缘检测，也属于先平滑后求导数的方法在与边缘垂直的方向上求零交叉。性能优于Marr算子。图像分割边缘检测Marr-HildrethCanny图像

10、分割边缘检测3 3、间断检测间断检测 间断检测包括点检测，线检测和边界检测三种。寻找间断最一般的方法是模板检测计算模板所包围区域的灰度级与模板系数的乘积之和，若大于阈值，则视为间断处。可用图示模板检测孤立点。孤立点：该点的灰度级与其背景的差异相当大，并且它所在的位置是一个均匀的或者近似均匀的区域。基本思想：如果一个孤立点与它周围的点很不相同，则很容易被上述模板检测到。在灰度级为常数的区域，模板响应为0。图像分割边缘检测线模板第一个模板对水平方向的线条（单象素宽）有最强的响应。第二个模板对于45线有最强响应。第三个模板对垂直线有最强响应。第四个模板对于45 线有最强响应。每个模板系数相加总和为0

11、，保证了在灰度级恒定的区域，模板响应为0。图像分割边缘检测4 4、边缘连接和边界、边缘连接和边界检测检测 理论上讲，边缘检测算子仅得到处在边缘上的像素点。实际上，由于噪声、不均匀照明而产生的边缘间断以及其他由于引入虚假的亮度间断所带来的影响，使得到的一组像素很少能完整地描绘一条边缘典型的做法是在使用边缘检测算法之后，使用连接过程将边缘像素组合成有意义的边缘。原理分析图像中每个边缘点的一个小邻域内像素的性质，将所有相似点连接起来，就是共同满足一定准则的像素组成的一条边缘。理想边缘图像分割边缘检测Eyxfyxf00,Ayxyx00,基于梯度一致性：要求：(x,y)是(x0,y0)的邻域，(x,y)

12、是该点的梯度方向角。 (a)输入图像(b)Sobel梯度的y分量(c)梯度的x分量(d)边缘连接的结果a bc d图像分割边缘检测迭代端点拟合：1.确定一个散布在两个特定边缘点A和B之间的边缘点的集合P；2.设A,B组成一条（初始）边缘线，计算其它边缘点到直线的距离；3.若所有点距直线的距离小于阈值，则拟合此直线为确定边缘；4.否则，找出距离最大的点（记为C）为新端点，将P分为包含AC 线段的边缘点子集合P1和包含BC线段的子集合P2；5. 重复3、4步骤，直到所有子集合的点到直线的距离小于阈值， 确定最后的各线段的端点。6. 利用每个线段间的边缘点，通过拟合确定边缘。多边形分段线性拟合原理图

13、像分割边缘检测ABCABCABCDD迭代端点拟合图像分割边缘检测Hough变换：利用图像的全局特性而直接检测目标轮廓。基本思想：点-线的对偶性(duality)图像变换前在图像空间，变换后在参数空间。在图像空间里，直线方程：y = ax+b ；给定 (x,y) 后，可以写成：b =ax+y ；即图像空间中的一个点对应参数空间 a-b 中的一条直线。图像空间中的两个点，确定了参数空间中两条直线的交点。图像分割边缘检测ab将参数空间划分，将参数空间划分，每一个每一个“格子格子”所所通过的直线数量，通过的直线数量，反映了图像空间中反映了图像空间中边缘点的边缘点的“共线性共线性”可以通过可以通过“投投

14、票胜出票胜出”方式确定方式确定最终的有效格子，最终的有效格子，即图像中的有效线即图像中的有效线段。段。图像分割边缘检测 使用等式y = ax+b表示一条直线带来的一个问题是，当直线接近垂直时，斜率接近无穷大。也就是说a的最大接近于无穷。解决的方法是采用直线的标准式（在-空间中）：sincosyx 格子的大小，决定了共线性的精度，以及累加器的数目。 Hough变换也适用于任何形式为g(v,c)=0的函数，这里v是坐标向量，c是系数向量。图像分割边缘检测(a) 航拍的红外线图像。(b)设置了门限的梯度图像。(c) 梯度图像的Hough变换。(d) 依据一定准则判断为相连的像素集合。准则：像素属于3

15、个具有最高计数的累加器单元；没有大于5个像素的间隙。a bc d图像分割区域生长二、基于区域的分割方法1 1. . 区域生长区域生长对每个需要分割的区域找一个种子像素作为生长的起点，将种子像素邻域中与其有相似性质的像素合并到种子像素所在的区域中。将这些新像素当作新的种子像素重复上述过程，直到再没有满足条件的像素被包括进来。(1)选择或确定一组能正确代表所需区域的种子像素；(2)确定在生长过程中将相邻像素合并的准则（相似性判据）；(3)确定在生长过程中停止的条件或规则。三个问题：图像分割区域生长(1)人机交互、直方图峰值、聚类中心等；(2)取决于问题本身以及图像数据的类型 ；(3) 要求图像分割

16、后区域的可描述性。基于区域灰度差方法基于区域灰度差方法以种子像素为中心比较它的邻域像素，如果灰度差小于预先确定的阈值，将它们合并；以新合并的像素为中心，重复上述过程，直到该连接区域生长停止；寻找新的种子，继续此过程。直到所有像素都有归属，则结束整个生长过程。图像分割区域生长XiihXH0)()()()(max21XHXHX基于区域内灰度分布统计性质基于区域内灰度分布统计性质 把像素分成互不重叠的小区域；比较邻接区域的累积灰度直方图，根据灰度分布的相似性进行合并。灰度直方图h(x)的累积灰度直方图定义： Kolmogorov-Smirnov检测XXHXH)()(21Smoothed-Differ

17、ence 检测 TTXh先从目标（结节）上一个种子点开始生长，当所谓孤立对照值ICI（Isolation Contrast Integral）达到最大时停止生长。 种子点选取：33或55区域中灰度最大点（或人工干预） 。 边界像素第j个点的灰度值目标结节第i个点的灰度值 图像分割区域生长iejjmiPgMaxPgICI )(miPg)(ejPg自适应阈值结节检测法自适应阈值结节检测法ICI值:图像分割区域生长是否找到种子像素寻找所有新边界像素找到边界像素中灰值最大点检测点计算新ICI值新ICI旧ICI停止生长增加检测点为目标结节点，新ICI取代旧ICI实质上由灰度最大点一圈圈向周围扩大（灰度逐

18、渐降低）遇平原灰度值开始增加时停止生长。可进一步扩展，ICI稍降低些也可不停止生长，可抗噪声干扰。可用别的特征代替ICI特征，在其它条件下获得更好结果。算法流程图像分割区域生长(a)焊缝的x射线图像(b)灰度值为255的种子点(c)区域生长的结果(d)对有缺陷的焊缝区域进行分割后得到的边界。 区域生长的准则：区域生长的准则：(1)(1)任何像素和种子之间的灰任何像素和种子之间的灰度值绝对差小于度值绝对差小于6565。依据是。依据是(a)(a)的直方图分布；的直方图分布；(2)(2)像素必须与此区域中至少像素必须与此区域中至少一个像素是一个像素是8 8连通的。连通的。a bc d图像分割区域生长2. 2. 分裂与合并分裂与合并先从整幅图像开始，通过不断分裂把图像分成任意大小且不重叠的区域（实际中常用四分树），直到每个区域的像素满足一致性准则，然后再合并一致性区域。图像分割区域生长图像分割区域生长3. 3. 分水岭算法分水岭算法(a)原图；(b)拓扑视图；(c)-(d)水从汇水盆地的最低点上涨；abcd图像分割区域生长(e)水进一步上涨；(f)两个相邻汇水盆地的水开始会合；在此会合处建造水坝。(g)随着水继续上涨，多处开始会合，水坝也在延伸；(h)最终完成的水坝，即分水线，完成分割。efgh图像分割区域生长Step1:在一个低灰度