数字图像处理 第八章

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
数字图像处理
Digital Image Processing
1
第八章 图像分割
2
8.1 概述
3
图像分割
特征
分割出来的各区域对某种性质例如灰度,纹
理而言具有相似性,区域内部是连通的的且 没有过多小孔。
区域边界是明确的 相邻区域对分割所依据的性质有明显的差异
4
图像分割
图像 预处理 图像 理解
幅值
方向角
mag ( f ) (Gx Gy )
2 2
1
2
Gy ( x, y) arctan( ) Gx
8
梯度算子
近似计算
M1 | Gx | | Gy |

M 2 Gx Gy
2
2
M Max(Gx , Gy )
数字图像处理中用差分代替微分
9
梯度算子
Z1 Z2 Z3
Z4 Z5 Z6 Z7 Z8 Z9 -1 0 0 1
Roberts算子
0
1
-1
0
Gx Z 9 Z 5 Gy Z8 Z 6
10
梯度算子
Z1 Z2 Z3
Z4 Z5 Z6 Z7 Z8 Z9 -1 -2 -1 0 1 0 2 0 1 -1 -2
Sobel算子
0 0
1 2
-1
分水岭对应于原始图像中的边缘
图8.24 分水岭形成示意图
42
自动阈值
MATLAB函数
L=watershed(f)
分水岭算法
原始图像
分水岭
分割结果
局部极小值
图8.25 不准确标记分水岭算法导致过分割
43
自动阈值
改进的分水岭算法
原始图像
原图像的距离变换
标记外部约束
标记内部约束
由标记内外部约 束重构的梯度图
有先验知识可以利用时,对含有复杂场景或自然景物等先验知识不足的图
像进行分割, 也可以取得较好的性能。但是,空间和时间开销都比较大。
46
区域生长
区域生长法主要考虑象素及其空间邻域象素之间的关系
开始时确定一个或多个象素点作为种子,然后按某种相似性准
则增长区域,逐步生成具有某种均匀性的空间区域,将相邻的具 有相似性质的象素或区域归并从而逐步增长区域,直至没有可以 归并的点或其它小区域为止。
Gx (Z 7 Z8 Z9 ) ( Z1 Z 2 Z 3 ) Gy (Z 3 Z 6 Z9 ) ( Z1 Z 4 Z 7 )
12
梯度算子
原图
Roberts算子
Sobel算子
Prewitt算子
13
拉普拉斯算子
二阶导数算子
2 2 f ( x , y ) f ( x, y) 2 f ( x, y) 2 x y 2
区域内象素的相似性度量可以包括平均灰度值、纹理、颜色等
信息。
47
区域生长
步骤
选择合适的种子点 确定相似性准则(生长准则) 确定生长停止条件
48
区域生长
实例
原始图像及 种子点位置
三个种子点区域 生长结果
原始图像及 种子点位置
四个种子点区域 生长结果
图8.29 区域生长
49
区域分裂
条件:如果区域的某些特性不满足一致性准则 开始:从图像的最大区域开始,一般情况下,是从整幅图像开始 注意:
然后再做微分运算,所不同的是平滑部分的权值有些差异,因此 对噪声具有一定的抑制能力,但不能完全排除检测结果中出现的 虚假边缘。虽然这两个算子边缘定位效果不错,但检测出的边缘 容易出现多像素宽度。
20
算子比较

Laplacian算子:是不依赖于边缘方向的二阶微分算子算子,
对图像中的阶跃型边缘点定位准确,该算子对噪声非常敏感,它使 噪声成分得到加强,这两个特性使得该算子容易丢失一部分边缘的
方向信息,造成一些不连续的检测边缘,同时抗噪声能力比较差。
21
算子比较

LOG算子:该算子首先用高斯函数对图像作平滑滤波处理,然后
才使用 Laplacian 算子检测边缘,因此克服了 Laplacian 算子抗噪声能力比 较差的缺点,但是在抑制噪声的同时也可能将原有的比较尖锐的边缘也
平滑掉了,造成这些尖锐边缘无法检被测到。应用LOG算子,高斯函数
1 ( 1 2 ) 2
(4) 重复步骤2、3,直到Ti+1和Ti的值差别小于某个给定值
36
自动阈值
迭代法
原始图像
迭代阈值二值化图象
图8.20 迭代式阈值二值化图像
37
自动阈值
Ostu法
Ostu方法是使类间方差最大而推导出的一种能自动确定阈值
的方法
百度文库
具有简单,处理速度快的特点
Matlab工具箱提供的graythresh函数求取阈值采用的正是Ostu
0
1
Gx (Z7 2Z8 Z9 ) (Z1 2Z2 Z3 ) Gy (Z3 2Z6 Z9 ) (Z1 2Z4 Z7 )
11
梯度算子
Z1 Z2 Z3
Z4 Z5 Z6 -1 -1 -1
Prewitt算子
-1
-1
0
0
1
1
0
1
0
1
0
1
Z7 Z8 Z9
-1
0
1
35
自动阈值
(1) 选择图像灰度的中值作为初始阈值Ti=T0。
迭代法
(2) 利用阈值Ti把图像分割成两部分区域, R1和R2,并计算其灰度均

u1
in
i 0 Ti
Ti
i
n
i 0
, u2
in
i Ti L 1 i Ti
L 1
i
i
n
i
(3) 计算新的阈值Ti+1
Ti 1
将这种想法应用于图像分割,就是要在灰度图像中找出不同的
吸水盆地和分水岭,由这些不同的吸引盆地和分水岭组成的区域即为
我们要分割的目标。
40
自动阈值
分水岭算法
(a)原始图像
(b)图像对应的拓扑地形图 图8.23 图像对应的拓扑表面图
41
自动阈值
分水岭算法
分水岭阈值选择算法可以看成是一种自适应的多阈值分割算法
作用
图 像
图像 识别 图像 分割
图8.1图像分割在整个图像处理过程中的作用
5
图像分割
边缘检测 边缘跟踪 阈值分割 区域分割 运动分割
本章要点
6
8.2 边缘检测
7
梯度算子
梯度算子是一阶导数算子
f Gx x f ( x, y ) Gy f y
23
8.3
边缘跟踪
24
基本步骤
从图像中一个边缘点出发,然后根据某种判别准则搜索下 一个边缘点以此跟踪出目标边界。
确定边界的起始搜索点,起始点的选择很关键,对某些图
像,选择不同的起始点会导致不同的结果。
确定合适边界判别准则和搜索准则,判别准则用于判断一
个点是不是边界点,搜索准则则指导如何搜索下一个边缘点。
微分
差分
2 f ( x, y) f ( x 1, y) f ( x 1, y) f ( x, y 1) f ( x, y 1) 4 f ( x, y)
14
拉普拉斯算子
0 1 0 1 1 1 1 -8 1 1 1 1
1 0
-4 1
1 0
图8.5 两种常用的拉普拉斯算子模板
T=210的二值化图像
图像直方图
34
自动阈值
迭代法
基本思想:开始时选择一个阈值作为初始估计值,然后按某种策略不
断地改进这一估计值,直到满足给定的准则为止。在迭代过程中,关 键之处在于选择什么样的阈值改进策略,好的阈值的改进策略应该具 备两个特征,一是能够快速收敛,二是在每一个迭代过程中,新产生 阈值优于上一次的阈值。
22
算子比较
Canny算子:Canny算子虽然是基于最优化思想推导出的边缘
检测算子,实际效果并不一定最优,原因在于理论和实际有许多不 一致的地方。该算子同样采用高斯函数对图像作平滑处理,因此具 有较强的抑制噪声能力,同样该算子也会将一些高频边缘平滑掉, 造成边缘丢失。Canny算子其后所采用用双阈值算法检测和连接边 缘,采用的多尺度检测和方向性搜索较LOG算子要好。

38
自动阈值
Ostu法
原始图像
Otsu法二值化图象
图8.21 Otsu法二值化图像
39
自动阈值
分水岭算法
分水岭算法(watershed)是一种借鉴了形态学理论的分割方法,
它将一幅图象看成为一个拓扑地形图,其中灰度值被认为是地形高度
值。高灰度值对应着山峰,低灰度值处对应着山谷。将水从任一处流 下,它会朝地势底的地方流动,直到某一局部低洼处才停下来,这个 低洼处被称为吸水盆地,最终所有的水会分聚在不同的吸水盆地, 吸 水盆地之间的山脊被称为分水岭,水从分水岭流下时,它朝不同的吸 水盆地流去的可能性是相等的。
分割结果
图8.26 准确标记的分水岭算法分割过程
44
8.5
区域分割
45
基本思想
阈值分割法由于没有或很少考虑空间关系,使多阈值选择受到限制
于区域的分割方法可以弥补这点不足,它利用的是图像的空间性质,
该方法认为分割出来的属于同一区域的像素应具有相似的性质,其概念是 相当直观的。
传统的区域分割算法有区域增长法和区域分裂合并法。该类方法在没

确定搜索的终止条件。
25
二值图像边界跟踪
26
灰度图像边界跟踪
27
Hough 变换
Hough变换可以用于将边缘像素连接起来得到边界曲线
优点在于受噪声和曲线间断的影响较小 在已知曲线形状的条件下,Hough变换实际上是利用分散
的边缘点进行曲线逼近,它也可看成是一种聚类分析技术
28
Hough 变换
15
拉普拉斯算子
16
Canny算子
出图像真实的边缘的同时要避免检测出现虚假的边缘。
基本思想
好的检测结果:对边缘的错误检测率要尽可能低,在检测
好的边缘定位精度:标记出的边缘位置要和图像上真正
边缘的位置尽量接近。
对同一边缘要有低的响应次数:有的算子会对一个边缘回
产生多个响应。也就是说图像上本来只有一个边缘点的,可是检 测出来就会出现多个边缘点。
直线检测
原始图像
二值化图像
细化图像
Hough 变换检测出的直线
29
8.4
阈值分割
30
基本原理
原始图像——f(x,y)
灰度阈值——T
阈值运算得二值图像——g(x,y)
1 如果 f [ x, y] T g[ x, y] 0 如果 f [ x, y] T
阈值选择直接影响分割效果,通常可以通过对灰度直方图
确定分裂准则(一致性准则)
确定分裂方法,即如何分裂区域,使得分裂后的子区域的特
性尽可能都满足一致性准则值。
50
区域分裂
算法
形成初始区域
对图像的每一个区域Ri,计算P(Ri),如果
P(Ri)=FALSE
则沿着某一合适的边界分裂区域
重复步骤2,当没有区域需分裂时,算法结束。
51
区域合并
对于图像中乡邻个区域,计算是否满足一致性谓词词,
若满足则合并为一个区域。
重复步骤2,直到没有区域可以合并,算法结束。
的分析来确定它的值。
31
阈值选择
利用灰度直方图求双峰或多峰 选择两峰之间的谷底作为阈值
32
人工阈值
人工选择法是通过人眼的观察,应用人对图像的知识,在分析图像 直方图的基础上,人工选出合适的阈值。也可以在人工选出阈值后 ,根据分割效果,不断的交互操作,从而选择出最佳的阈值。
33
人工阈值
原始图像
T=155的二值化图像

克服噪声的影响
17
Canny算子
用高斯滤波器平滑图像
计算滤波后图像梯度的幅值和方向
算法步骤
对梯度幅值应用非极大值抑制,其过程为找处图像梯度中的
局部极大值点,把其它非局部极大值点置零以得到得到细化的 边缘
用双阈值算法检测和连接边缘,使用两个阈值T1和
T2(T1>T2),T1用来找到每条线段,T2用来在这些线段的两个 方向上延伸寻找边缘的断裂处,并连接这些边缘。
基本思想
合并运算就是把相邻的具有相似性质的区域合成为一个区域
合并算法中最重要的运算是确定两个区域的相似性
评判区域相似性方法有许多,评判相似性的方法可以基于区
域的灰度值,也可以基于区域边界的强弱性等因素。一种简单 的方法是比较它们的灰度均值
52
区域合并
算法
使用某种方法进行图像的初始区域分割。
中方差参数的选择很关键,对图像边缘检测效果有很大的影响。高斯滤 波器为低通滤波器,越大,通频带越窄,对较高频率的噪声的抑制作用 越大,避免了虚假边缘的检出,同时信号的边缘也被平滑了,造成某些 边缘点的丢失。反之,越小,通频带越宽,可以检测到的图像更高频率 的细节,但对噪声的抑制能力相对下降,容易出现虚假边缘。因此,应 用LOG算子,为取得更佳的效果,对于不同图像应选择不同参数。
18
Canny算子
实例
19
算子比较

Roberts算子:Roberts算子利用局部差分算子寻找边缘,边
缘定位精度较高,但容易丢失一部分边缘,同时由于图像没经过 平滑处理,因此不具备能抑制噪声能力。该算子对具有陡峭边缘 且含噪声少的图像效果较好。

Sobel算子和Prewitt算子:都是对图像先做加权平滑处理,
相关文档
最新文档