数字图像处理 第八讲 图像分割课件

1、图像分割第八讲 图像分割8.1 概述8.2 图像边缘检测8.3 阈值法分割图像分割8.1 概述图像分割的目的: 将一幅图像中的各成分分离成若干与景物中的实际物体相对应的子集。 例：地貌遥感图像，平原、山区、水域、森林等；医学图像，细胞图像中分离出单个细胞。图像分割的基础: 像素间的相似性和不连续性。图像分割图像分割的方法： 从分割依据角度出发，图像分割方法可分为相似性分割和非连续性分割。一、相似性分割： 将具有同一灰度级或相同组织结构的像素聚集在一起，形成图像中的不同区域。这种基于相似性原理的方法也称为基于区域相关的分割技术。图像分割二、非连续性分割： 首先检测局部不连续性，然后将它们连接起来

2、形成边界，这些边界把图像分以不同的区域。这种基于不连续性原理检出物体边缘的方法称为基于点相关的分割技术。 这两种方法是互补的。有时将它们有机地结合起来，以求得到更好的分割效果。 本讲主要介绍边缘检测、阈值分割。图像分割8.2 图像边缘检测 图像边缘是图像最基本的特征之一，携带着一幅图像的大部分信息，因此在图像分析时，通过图像检测提取出它的边缘。 图像边缘是以图像局部特性的不连续性的形式出现的，从本质上说，边缘线意味着一个区域的终结和另一个区域的开始，是不同目标区域或组成部分的分界线。图像分割图81 边缘灰度变化的几种类型图像分割 边缘检测的实质：采用某种算法来提取出图像中对象与背景间的交界线。

3、 边缘图像灰度的变化情况可以用图像灰度分布的梯度来反映，根据这一特点，提出多种边缘检测算子，如Robert算子、Sobel算子、Laplacian算子等。 最早的边缘检测方法都是基于像素的数值导数的，在数字图像中应用差分代替导数运算。图像分割对于数字图像，可用一阶差分替代一阶微分：此时梯度的幅度可表示为：图像分割为了避免平方和、开方运算，可将上式表示为： 利用图像灰度的一阶导数算子在灰度迅速变化处得到的极值来进行奇异点的检测。它在某一点的值就代表该点的边缘强度，可以通过对这些值设置阈值来进一步得到边缘图像。但该算子具有方向性，计算不便，很少采用。图像分割 在图像边缘提取中，常使用罗伯特(Rob

4、ert)梯度算子和Sobel梯度算子检测边缘。一、罗伯特(Robert)梯度算子（22个像素） Robert梯度采用对角方向相邻两像素之差，即： 将上面二者代入到梯度的幅度表达式中，即可求得图像的梯度幅度。图像分割 选取适当的门限TH，作如下判断：Gf(x,y)TH，(x,y)为阶跃状边缘点。二、Sobel梯度算子（33个像素） 先做加权平均，再作微分，即图像分割 代入到梯度的幅度表达式中，即可求得图像的梯度幅度。 选取适当的门限TH，作如下判断：Gf(x,y)TH，(x,y)为阶跃状边缘点。 Sobel算子是边缘检测器中最常见的算子之一。该算子不但产生较好的边缘检测效果，同时因为引入了局部平

5、均，使其受到噪声的影响也比较小，也正因为平均，有时会检测出许多伪边缘。图像分割 拉普拉斯(Laplacian)算子是不依赖于边缘方向的二阶微分算子。它是一个标量而不是向量，具有旋转不变即各向同性的性质，在图像处理中经常被用来提取图像的边缘。其表示式为 在数字图像中使用差分有：图像分割 将二者合并为一个算子有：图像增强 f(x-1,y) f(x,y-1) f(x,y)f(x,y+1) f(x+1,y) 图像分割 当拉普拉斯算子输出出现过零点时就表明有边缘存在。该算子有两个缺点：其一就是边缘方向信息的丢失，其二它是二阶差分，双倍加强了图像中噪声的影响。 改进的LOG算法： 在进行拉普拉斯运算前先进

6、行平滑去噪，然后再提取边缘。平滑去噪采用高斯滤波器，然后与拉普拉斯边缘检测合并在一起，形成LOG（Laplacian Of Gaussian）。图像分割8.3 阈值法分割 利用图像中要提取的目标物与其背景在灰度特性上的差异，把图像视为具有不同灰度级的两类区域(目标和背景)的组合，选取一个合适的阈值，以确定图像中每一个像素点应该属于目标还是背景区域，从而产生相应的二值图像。 优点：大量压缩数据、减少存储容量、大大简化其后的分析和处理步骤。图像分割 要从复杂的景物中分辨出目标并将其形状完整地提取出来，阈值的选取是阈值分割技术的关键。如果阈值选取过高，则过多的目标点被误认为背景；阈值选得过低，则会出

7、现相反的情况。至今还未能找到一种对所有图像都能有效分割的阈值选取方法。 阈值分割的方法：二值化分割、半阈值分割、区域生长法分割。图像分割 针对不同的处理对象，可分别选用以下定义式：图像分割 二值分割时的处理流程有两种方法：一、一次比较法 顺序读取各像素点灰度级与阈值比较，凡符合约定关系的用1代替，否则以0代替。（最常用）二、二次扫描法 先按行的顺序扫描比较，得到过渡图像g1(x,y)，再按列的顺序扫描比较，得到过渡图像g2(x,y)，然后将这两个结果相或，便为二次扫描图像分割输出g(x,y)。图像分割 B1灰度区为背景区，B2灰度区为目标区，阈值确定在背景与物体图像灰度交界处的物体图像灰度一侧

8、。图85 直方图双峰图像分割 2、将图像中的每一行扫描，产生的图像的灰度将遵循如下规则： 3、将图像中的每一列扫描，产生的图像的灰度将遵循如下规则：图像分割 4、将g1(x,y)和g2(x,y)相或，得到最后输出g(x,y)。 二次扫描法分割图像，可以检测出水平和垂直方向上的边界。例：下图是一幅1216，背景灰度为12、图像灰 度为37的离散图像。现令灰度12为B1区、灰度37为B2区，vc=8,vb=0。图像分割1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 11 1 1 1 2 2 3 3 1 1 1 12 2 1 1 1 3 3 3 3 2 2

9、 12 2 1 2 3 4 5 3 4 3 1 21 2 2 4 4 5 6 5 6 5 3 11 1 5 5 6 7 6 5 3 4 4 4 1 2 4 4 5 5 5 4 5 5 4 31 2 2 3 4 5 6 7 7 6 3 21 1 2 4 4 4 3 6 7 6 1 22 1 1 1 5 5 6 6 6 2 1 11 1 2 2 1 3 4 4 3 2 2 11 1 1 1 1 1 3 5 1 1 1 11 1 1 2 2 1 1 4 1 1 1 11 1 1 1 2 2 1 1 1 1 1 11 1 1 1 1 2 2 1 1 1 1 11 1 1 1 1 2 2 1 1 1 1

10、 1 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 8 8 0 0 0 00 0 0 0 0 8 0 0 8 0 0 00 0 0 0 8 0 0 0 0 8 0 00 0 0 8 0 0 0 0 0 0 8 00 0 8 0 0 0 0 0 0 0 0 8 0 0 0 0 0 0 0 0 0 0 0 00 0 8 0 0 0 0 0 0 0 0 80 0 0 0 0 0 0 0 0 0 8 00 0 0 8 0 0 0 0 0 8 0 00 0 0 0 8 0 0 0 0 0 0 00 0 0 0 0 8 0 0 8 0 0 00 0 0 0 0 0 8 0 0 0

11、0 00 0 0 0 0 0 0 8 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0原图像第一次扫描图像分割0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 8 8 8 0 0 00 0 0 0 0 8 0 0 8 8 0 00 0 0 0 8 0 0 0 0 8 8 00 0 0 8 0 0 0 0 0 0 8 80 0 8 0 0 0 0 0 0 0 0 8 0 0 8 0 0 0 0 0 0 0 0 80 0 8 8 0 0 0 0 0 0 0 80 0 0 8 0 0 0 0 0 0 8 00 0 0 8

12、 8 0 0 0 0 8 0 00 0 0 0 8 8 0 0 0 8 0 00 0 0 0 0 8 8 8 8 0 0 00 0 0 0 0 0 8 8 8 0 0 00 0 0 0 0 0 0 8 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0最后输出结果1 0 0 0 0 0 0 0 0 0 0 01 0 0 0 0 0 0 0 0 0 0 01 0 0 0 0 0 8 0 8 0 0 02 0 0 0 0 8 0 0 0 8 0 02 0 0 0 8 0 0 0 0 0 8 01 0 0 8 0 0 0 0 0 0 0 81 0 8 0 0 0 0 0 0 0 0 0

13、1 0 8 0 0 0 0 0 0 0 0 01 0 0 0 8 0 0 0 0 0 0 81 0 0 8 0 0 0 0 0 0 8 02 0 0 0 8 0 0 0 0 8 0 01 0 0 0 0 8 0 0 0 8 0 01 0 0 0 0 0 8 0 8 0 0 01 0 0 0 0 0 0 8 8 0 0 01 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0第二次扫描图像分割8.3.2 半阈值分割 该分割方法是将灰度级低于某一阈值的像素灰度级变换为0，其他灰度级像素保持不变的分割技术，表达式为：图像分割8.3.3 最佳阈值 前面介绍了二值

14、化分割和半阈值分割，知道选择最佳阈值是阈值法分割图像的首要工作，如何选择最佳阈值呢？ 令图像由背景和物体两部分组成，背景所对应的灰度级概率分布密度函数为P1(f) ，其均值为1 ， 方差为12；物体所对应的灰度级概率分布密度函数为P2(f)，其均值为2，方差为22。背景和图像的面积比分别为P1、P2(P1+P2=1)。图像分割 假定暗区相当于背景，亮区相当于物体。凡灰度级低于阈值T的像素均作为背景来对待，凡灰度级高于阈值T的像素均作为物体。 把物体上的点错分类成为背景上的点的概率为： 把背景上的点错分类成为物体上的点的概率为：图像分割 总的误差概率为： 为了获得最佳阈值，将E(T)对T求微分，

15、并使结果为零。从而有： 设P1(f) ，P2(f)均为正态分布，且二者的方差12、22相等，同时面积比P1P21/2，则最佳阈值T为：图像分割 可见，在背景和目标图像的先验概率相等这一特定条件下，最佳阈值就是背景灰度均值与目标图像灰度均值之均值。（求解过程忽略）图像分割8.3.4 区域生长法分割 区域生长法又叫区域扩张法。该方法从被称为“生长点”（可以是单个像素，也可以是某个小区域）地方开始，搜索其邻域，把符合接收准则的点或子域归并进来，形成新的生长点，反复操作，直到当前区域不能再合并为止。 简单连接、子区合并法图像分割一、简单连接 基于单个像素的区域生长法。 从满足检验点（以下简称连接核）开

16、始，考虑其周围（4像素邻域或8像素邻域）的不属于任何一个域的点。如果其特性符合接收准则，就把它作为同一个区域加以合并，形成连接核A1，既而检测A1周围的点，把符合接收准则的点并入A1，产生新的连接核A2。重复上述过程，直到没有可并入的点时，生长过程结束。图像分割 这里的接收准则可设置为邻近像素与Ai灰度均值之差作为门限T，凡灰度小于该门限的邻点均认为是区域附属点而并入。 例：设接收准则为可并入点的灰度与连接核灰度均值差小于2。图像分割二、子区合并法 以小区域内的灰度分布的相似性进行区域合并。 步骤如下：1、将图像分割成MN个，大小为LxLy的子区。2、从左上角第一个子区开始，分别计算子区和相邻子区的灰度统计特量（直方图），然后作相似性判定。若两者的灰度分布相似且符合接收准则，相邻子区并入当前子区，形成下一轮判定合并时的当前子区。如果某个相邻子区不符合接收则，图像分割 就打上“未分割标记”。继续新的一轮判定，使当前子区不断生长，直到没有可