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

1、医学图像处理信自学院生医系24.1 图像分割概述4.2 间断检测4.3 边缘连接和边界检测4.4 阈值分割法4.5 基于区域的分割4.6 数学形态学图像处理第四章第四章 医学图像分割医学图像分割34.5 基于区域的分割基于区域的分割44.5.1 基本概念基本概念n目标：将区域r划分为若干个子区域r1,r2,rn，这些子区域满足5个条件：niirr154.5.1 基本概念基本概念即任何区域满足如下谓词（例如任何两个相邻区域不能合并成单一区域（例如必满足谓词：n一致性谓词p 定义了在区域 r 上的所有点与区域模型的相似程度。64.5.2 区域生长区域生长n区域生长技术利用原理：同类象素具有相似性。

2、n对每一象素定义对应的特征向量 其中， 都是与（i, j）有关的性质，如灰度值、梯度幅值、是否为边缘点等 tnjipjipjipjip),(,),(),(),(21nppp,2174.5.2 区域生长区域生长n衡量两个象素（i, j）（k, l）是否为同类可以通过比较它们的特征向量是否具有相似性来确定.n比较或n若a (i, j, k, l)或 l (i, j, k, l)足够小即认为象素是相似的（即为同类）。nmmmlkpjiplkjia1),(),(),(nmmmlkpjiplkjil12),(),(),(84.5.2 区域生长区域生长n1. 单一型链结的区域生长n2. 混合型链结的区域生

3、长n3. 登山算法n4. 分水岭算法94.5.2 区域生长区域生长1.单一型链结的区域生长算法实现：104.5.2 区域生长区域生长n1. 单一型链结的区域生长 区域a 区域b 种子像素 种子像素114.5.2 区域生长区域生长n1. 单一型链结的区域生长上图给出一个简单的例子。此例的相似性准则是邻近点（4邻域）的灰度级与物体的平均灰度级的差小于2。图中被接受的点和起始点均用下划线标出， 其中（a）图是输入图像；（b）图是第一步接受的邻近点； （c）图是第二步接受的邻近点； （d）图是从6开始生成的结果。12n举例：一幅图像背景部分的均值为25，方差为625，在背景上分布着一些互不重叠的均值为

4、150，方差为400的小目标。设所有目标合起来约占图像总面积的20%，提出1个基于区域生长的分割算法将这些目标分割出来。 13算法描述算法描述n从左至右，从上至下扫描图像。n若扫描到灰度值大于150的象素，取为种子点，进行区域生长。n生长准则为将相邻的灰度值与已有区域的平均灰度值的差小于60（3）的象素扩展进来。n若不能再生长，标记已生长区域。n若扫描到图像右下角，结束；否则回到继续。 说明144.5.2 区域生长区域生长n2. 混合型链结的区域生长232221)()()(aywbxwyxws考虑象素的邻域，在此邻域上定义象素的特性矢量，特性矢量接近的两个象素相似。（1） 邻域分享技术：研究某

5、象素的邻域，作出相似邻域表。如果两个象素在对方的相似邻域表中，且表中有足够多的象素，则可将二者连接。s值小则可以连接。（2）依赖边缘确定两象素是否可以连接。定义两个象素特征矢量可用（x，a）、（y，b）表示，x，y为两者的灰度，a，b为两者的邻域灰度均值，相似性为s，则321,www为非负的权值思考题154.5.2 区域生长区域生长n3. 登山算法(1)灰度极大值点作为中心点；(2)16个等角度间隔方向上检测斜率值最大的点作为边缘点；(3)以16个边缘点作为种子点进行区域生长。164.5.2 区域生长区域生长区域生长算法的优点：n 计算简单n 特别适用于分割小的结构如肿瘤和伤疤n 区域生长也很

6、少单独使用，往往是与其他分割方法一起使用。区域生长的缺点：n 有时需要人工交互以获得种子点，这样使用者必须在每个需要抽取出的区域中植入一个种子点。n 对噪声敏感，导致抽取出的区域有空洞或者将原本分开的区域连接起来。174.5.2 区域生长区域生长n4. 分水岭算法n分水岭（watershed，也称分水线/水线）把图像看成3-d地形的表示，即2-d的地基（对应图像空间）加上第3维的高度（对应图像灰度）n计算过程是串行的，得到的是目标的边界18建立不同目标间的分水岭谷底孔分水岭三类点：（a）属于局部性最小值的点；（b）当一滴水放在某点的位置上时，水一定会下落到一个单一的最小值点；（c）当水处于某点

7、位置上时，水会等概率流向不止一个这样的最小值点。194.5.2 区域生长区域生长4. 分水岭算法20待分割图像 f (x, y)，其梯度图像为g(x, y)用m1, m2, , mr表示g(x, y)中各局部极小值的象素位置，c(mi)为与mi对应的区域中的象素坐标集合。用n表示当前灰度阈值，tn代表记为(u, v)的象素集合，g(u, v) n，对mi所在的区域，其中满足条件的坐标集合cn(mi)可看作一幅二值图像 ),(| ),(nvugvunt)()(ntmcmciin分水岭计算步骤分水岭计算步骤21)()(ntmcmciin22用cn代表在灰度阈值为n时图像中所有满足条件的象素 cma

8、x + 1将是所有区域的并集cn 1是cn的子集，cn是tn的子集，所以cn 1又是tn的子集 riinmcnc1)(riimcc11max)( 1max23令 q 代表tn中的连通组元集合，对每个连通组元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 中建水坝 24基于形态学分水岭的分割基于形态学分水岭的分割水坝的构造水坝的构造（a）在第n-1个阶段淹没的

9、汇水盆地的两个部分，（b）淹没的第n个阶段，显示出两个盆地间的水已经溢出，（c）用于膨胀的结构元素（数学形态学方法），（d）扩展的结果和水坝的构造25 原始图 阈值分割 分水岭 叠加轮廓补充题26n分水岭算法不是简单地将图像在最佳灰度级进行阈值处理，而是从一个偏低但仍然能正确分割各个物体的阈值开始。然后随着阈值逐渐上升到最佳值，使各个物体不会被合并。这个方法可以解决那些由于物体靠得太近而不能用全局阈值解决的问题。只要也只有所采用最初的阈值进行分割的结果是正确的，那么，最后的分割也是正确的。也就是说，图像中每个实际物体都有相应的边界。n 最初和最终的阈值灰度级都必须很好地选取。如果初始的阈值太低

10、，那么低对比度的物体开始时会被丢失，然后随着阈值的增加就会和相邻的物体合并。如果初始阈值太高，物体一开始便会被合井。最终的阈值决定了最后的边界与实际物体的吻合程度。本章中讨论的阈值选择方法在确定这两个值时很有用。274.5.3 区域分裂与合并区域分裂与合并n算法实现：284.5.3 区域分裂与合并区域分裂与合并n算法实现：294.5.3 区域分裂与合并区域分裂与合并n算法实现（实际应用中还可作修改）：304.5.3 区域分裂与合并区域分裂与合并n在区域合并方法中，输入图像往往先被分为多个相似的区 域，然后类似的相邻区域根据某种判断准则迭代地进行合并。n 在区域分裂技术中，整个图像先被看成一个区域，然后区域不断被分裂为四个矩形区域，直到每个区域内部都是相似的。n 在区域的分裂合并方法中，先从整幅图像进行分裂，然后将相邻的区域进行合并。n 分裂合并方法不需要预先指定种子点，它的研究重点是分裂和合并规则的设计。但是，分裂可能会使分割区域的边界被破坏。31n正态分布时区间（-1,+1）的面积占总面积的68.27%；n正态分布时区间（-1.96,+1.96）的面积占总面积的95%；n正态分布时区间（-2.58,+2.58）的面积占总面积的99%。 32思考