MedicalImage9医学图像分割基于区域的分割课件

1、医学图像处理信自学院生医系MedicalImage9医学图像分割基于区域的分割1医学图像处理信自学院生医系MedicalImage9医学图像4.1 图像分割概述4.2 间断检测4.3 边缘连接和边界检测4.4 阈值分割法4.5 基于区域的分割4.6 数学形态学图像处理第四章 医学图像分割MedicalImage9医学图像分割基于区域的分割24.1 图像分割概述第四章 医学图像分割MedicalIma4.5 基于区域的分割4.5.1 基本概念4.5.2 区域生长4.5.3 区域分裂与合并MedicalImage9医学图像分割基于区域的分割34.5 基于区域的分割MedicalImage9医学图像

2、分割4.5.1 基本概念目标：将区域R划分为若干个子区域R1,R2,Rn，这些子区域满足5个条件：1)完备性：2)连通性：每个Ri都是一个连通区域3)独立性：对于任意ij，RiRj= MedicalImage9医学图像分割基于区域的分割44.5.1 基本概念目标：将区域R划分为若干个子区域R1,R4.5.1 基本概念4)单一性：即任何区域满足如下谓词（例如每个区域内的灰度级相等）P（Ri）= TRUE，i = 1,2,n5)互斥性：任何两个相邻区域不能合并成单一区域（例如任两个区域的灰度级不等），必满足谓词：P（RiRj）= FALSE，ij一致性谓词P 定义了在区域 R 上的所有点与区域模型

3、的相似程度。MedicalImage9医学图像分割基于区域的分割54.5.1 基本概念MedicalImage9医学图像分割基4.5.2 区域生长区域生长技术利用原理：同类象素具有相似性。对每一象素定义对应的特征向量 其中， 都是与（i, j）有关的性质，如灰度值、梯度幅值、是否为边缘点等 MedicalImage9医学图像分割基于区域的分割64.5.2 区域生长区域生长技术利用原理：同类象素具有相似性4.5.2 区域生长衡量两个象素（i, j）（k, l）是否为同类可以通过比较它们的特征向量是否具有相似性来确定.比较或若A (i, j, k, l)或 L (i, j, k, l)足够小即认为

4、象素是相似的（即为同类）。MedicalImage9医学图像分割基于区域的分割74.5.2 区域生长衡量两个象素（i, j）（k, l）是否4.5.2 区域生长1. 单一型链结的区域生长2. 混合型链结的区域生长3. 登山算法4. 分水岭算法MedicalImage9医学图像分割基于区域的分割84.5.2 区域生长1. 单一型链结的区域生长Medical4.5.2 区域生长1.单一型链结的区域生长算法实现：1）根据图像的不同应用选择一个或一组种子，它或 者是最亮或最暗的点，或者是位于点簇中心的点。2）选择一个相似性准则。(生长条件) （灰度级、 彩色、 纹理、 梯度等特性相似）3）从该种子开始

5、向外扩张，不断将与集合中各个像 素连通、且满足相似性准则的像素加入集合。4）上一过程进行到不再有满足条件的新结点加入集 合为止。（终止准则）MedicalImage9医学图像分割基于区域的分割94.5.2 区域生长1.单一型链结的区域生长MedicalI4.5.2 区域生长1. 单一型链结的区域生长 区域A 区域B 种子像素 种子像素MedicalImage9医学图像分割基于区域的分割104.5.2 区域生长1. 单一型链结的区域生长 区域A 区域4.5.2 区域生长1. 单一型链结的区域生长上图给出一个简单的例子。此例的相似性准则是邻近点（4邻域）的灰度级与物体的平均灰度级的差小于2。图中被

6、接受的点和起始点均用下划线标出， 其中（a）图是输入图像；（b）图是第一步接受的邻近点； （c）图是第二步接受的邻近点； （d）图是从6开始生成的结果。MedicalImage9医学图像分割基于区域的分割114.5.2 区域生长1. 单一型链结的区域生长上图给出一个简举例：一幅图像背景部分的均值为25，方差为625，在背景上分布着一些互不重叠的均值为150，方差为400的小目标。设所有目标合起来约占图像总面积的20%，提出1个基于区域生长的分割算法将这些目标分割出来。 MedicalImage9医学图像分割基于区域的分割12举例：一幅图像背景部分的均值为25，方差为625，在背景上分算法描述从

7、左至右，从上至下扫描图像。若扫描到灰度值大于150的象素，取为种子点，进行区域生长。生长准则为将相邻的灰度值与已有区域的平均灰度值的差小于60（3）的象素扩展进来。若不能再生长，标记已生长区域。若扫描到图像右下角，结束；否则回到继续。 说明MedicalImage9医学图像分割基于区域的分割13算法描述从左至右，从上至下扫描图像。说明MedicalIm4.5.2 区域生长2. 混合型链结的区域生长考虑象素的邻域，在此邻域上定义象素的特性矢量，特性矢量接近的两个象素相似。（1） 邻域分享技术：研究某象素的邻域，作出相似邻域表。如果两个象素在对方的相似邻域表中，且表中有足够多的象素，则可将二者连接

8、。S值小则可以连接。（2）依赖边缘确定两象素是否可以连接。定义两个象素特征矢量可用（x，a）、（y，b）表示，x，y为两者的灰度，a，b为两者的邻域灰度均值，相似性为S，则为非负的权值思考题MedicalImage9医学图像分割基于区域的分割144.5.2 区域生长2. 混合型链结的区域生长考虑象素的邻域4.5.2 区域生长3. 登山算法(1)灰度极大值点作为中心点；(2)16个等角度间隔方向上检测斜率值最大的点作为边缘点；(3)以16个边缘点作为种子点进行区域生长。MedicalImage9医学图像分割基于区域的分割154.5.2 区域生长3. 登山算法(1)灰度极大值点作为中心4.5.2

9、区域生长区域生长算法的优点： 计算简单 特别适用于分割小的结构如肿瘤和伤疤 区域生长也很少单独使用，往往是与其他分割方法一起使用。区域生长的缺点： 有时需要人工交互以获得种子点，这样使用者必须在每个需要抽取出的区域中植入一个种子点。 对噪声敏感，导致抽取出的区域有空洞或者将原本分开的区域连接起来。MedicalImage9医学图像分割基于区域的分割164.5.2 区域生长区域生长算法的优点：MedicalIma4.5.2 区域生长4. 分水岭算法分水岭（watershed，也称分水线/水线）把图像看成3-D地形的表示，即2-D的地基（对应图像空间）加上第3维的高度（对应图像灰度）计算过程是串行

10、的，得到的是目标的边界MedicalImage9医学图像分割基于区域的分割174.5.2 区域生长4. 分水岭算法分水岭（watersh建立不同目标间的分水岭谷底孔分水岭三类点：（a）属于局部性最小值的点；（b）当一滴水放在某点的位置上时，水一定会下落到一个单一的最小值点；（c）当水处于某点位置上时，水会等概率流向不止一个这样的最小值点。MedicalImage9医学图像分割基于区域的分割18建立不同目标间的分水岭谷底孔分水岭三类点：MedicalI4.5.2 区域生长4. 分水岭算法 图像最初在一个低灰度值上二值化。该灰度值把图像分割成正确数目的物体，但它们的边界偏向物体内部。随后阈值逐渐增

11、加，每一次增加一个灰度级。物体的边界将随着阈值增加而扩展。当边界相互接触时，这些物体并没有合并。因此，这些初次接触的点变成了相邻物体间的最终边界。这个过程在阈值达到背景的灰度级之前终止一 也就是说，在被恰当分割的物体的边界正确地确定时终止。MedicalImage9医学图像分割基于区域的分割194.5.2 区域生长4. 分水岭算法 图像最初在一个低灰度值待分割图像 f (x, y)，其梯度图像为g(x, y)用M1, M2, , MR表示g(x, y)中各局部极小值的象素位置，C(Mi)为与Mi对应的区域中的象素坐标集合。用n表示当前灰度阈值，Tn代表记为(u, v)的象素集合，g(u, v)

12、 n，对Mi所在的区域，其中满足条件的坐标集合Cn(Mi)可看作一幅二值图像 分水岭计算步骤MedicalImage9医学图像分割基于区域的分割20分水岭计算步骤MedicalImage9医学图像分割基于区域MedicalImage9医学图像分割基于区域的分割21MedicalImage9医学图像分割基于区域的分割21用Cn代表在灰度阈值为n时图像中所有满足条件的象素 Cmax + 1将是所有区域的并集Cn 1是Cn的子集，Cn是Tn的子集，所以Cn 1又是Tn的子集 MedicalImage9医学图像分割基于区域的分割22用Cn代表在灰度阈值为n时图像中所有满足条件的象素 M令 Q 代表Tn

13、中的连通组元集合，对每个连通组元q Qn，有3种可能性： (1)q Cn 1是1个空集 (2)q Cn 1里包含Cn 1中的一个连通组元 (3)q Cn 1里包含Cn 1中一个以上的连通组元分别处理： (1)Cn可由把连通组元 q 加到Cn 1中得到 (2)Cn可由把连通组元 q 加到Cn 1中得到 (3)需要在 q 中建水坝 MedicalImage9医学图像分割基于区域的分割23MedicalImage9医学图像分割基于区域的分割23基于形态学分水岭的分割水坝的构造（a）在第n-1个阶段淹没的汇水盆地的两个部分，（b）淹没的第n个阶段，显示出两个盆地间的水已经溢出，（c）用于膨胀的结构元素

14、（数学形态学方法），（d）扩展的结果和水坝的构造MedicalImage9医学图像分割基于区域的分割24基于形态学分水岭的分割水坝的构造（a）在第n-1个阶段淹 原始图 阈值分割 分水岭 叠加轮廓补充题MedicalImage9医学图像分割基于区域的分割25补充题MedicalImage9医学图像分割基于区域的分割2分水岭算法不是简单地将图像在最佳灰度级进行阈值处理，而是从一个偏低但仍然能正确分割各个物体的阈值开始。然后随着阈值逐渐上升到最佳值，使各个物体不会被合并。这个方法可以解决那些由于物体靠得太近而不能用全局阈值解决的问题。只要也只有所采用最初的阈值进行分割的结果是正确的，那么，最后的分

15、割也是正确的。也就是说，图像中每个实际物体都有相应的边界。 最初和最终的阈值灰度级都必须很好地选取。如果初始的阈值太低，那么低对比度的物体开始时会被丢失，然后随着阈值的增加就会和相邻的物体合并。如果初始阈值太高，物体一开始便会被合井。最终的阈值决定了最后的边界与实际物体的吻合程度。本章中讨论的阈值选择方法在确定这两个值时很有用。MedicalImage9医学图像分割基于区域的分割26分水岭算法不是简单地将图像在最佳灰度级进行阈值处理，而是从一4.5.3 区域分裂与合并算法实现：1）对图像中灰度级不同的区域，均分为四个子区域2）如果相邻的子区域所有像素的灰度级相同，则将其合并3）反复进行上两步操

16、作，直至不再有新的分裂与合并为止MedicalImage9医学图像分割基于区域的分割274.5.3 区域分裂与合并算法实现：MedicalImage4.5.3 区域分裂与合并算法实现：MedicalImage9医学图像分割基于区域的分割284.5.3 区域分裂与合并算法实现：MedicalImage4.5.3 区域分裂与合并算法实现（实际应用中还可作修改）：P(Ri)的定义为：1）区域内多于80%的像素满足不等式|zj-mi|=2i，其中：zj是区域Ri中第j个点的灰度级， mi是该区域的平均灰度级， i是区域的灰度级的标准方差。2）当P(Ri)=TRUE时，将区域内所有像素的灰度级置为mi。

17、MedicalImage9医学图像分割基于区域的分割294.5.3 区域分裂与合并算法实现（实际应用中还可作修改）：4.5.3 区域分裂与合并在区域合并方法中，输入图像往往先被分为多个相似的区 域，然后类似的相邻区域根据某种判断准则迭代地进行合并。 在区域分裂技术中，整个图像先被看成一个区域，然后区域不断被分裂为四个矩形区域，直到每个区域内部都是相似的。 在区域的分裂合并方法中，先从整幅图像进行分裂，然后将相邻的区域进行合并。 分裂合并方法不需要预先指定种子点，它的研究重点是分裂和合并规则的设计。但是，分裂可能会使分割区域的边界被破坏。MedicalImage9医学图像分割基于区域的分割304

18、.5.3 区域分裂与合并在区域合并方法中，输入图像往往先被正态分布时区间（-1,+1）的面积占总面积的68.27%；正态分布时区间（-1.96,+1.96）的面积占总面积的95%；正态分布时区间（-2.58,+2.58）的面积占总面积的99%。 MedicalImage9医学图像分割基于区域的分割31正态分布时区间（-1,+1）的面积占总面积的68.思考题MedicalImage9医学图像分割基于区域的分割32思考题MedicalImage9医学图像分割基于区域的分割3MedicalImage9医学图像分割基于区域的分割33MedicalImage9医学图像分割基于区域的分割33MedicalImage9医学图像分割基于区域的分割34MedicalImage9医学图像分割基于区域的分割34作业1、对下图用区域生长法进行分割，其中画线象素为种子点，令新生长进