第4章医学图像分割ppt课件.ppt

《第4章医学图像分割ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第4章医学图像分割4.1医学图像分割概念所谓图象分割就是根据某种均匀性（或一致性）的原则将图象分成若干个有意义的部分，使得每一部分都符合某种一致性的要求，而任意两个相邻部分的合并都会破坏这种一致性。图象的分割在很多情况下可归结为图象象素点的分类问题。图象分割的一个重要手段就是利用边缘检测技术。这种分割方法是基于物体与背景之间在灰度（或纹理）特性上存在着某种不连续性（或突变性）。边缘是指它的两侧分属于两个区域，每个区域内部的特性相对比较均匀一致，而两个区域之间在特性上则存在一定差异。图象区域的边界（或称轮廓）是指包围着连通区域的一条封闭的边缘曲线。4.2阈值分割技术4.2.1全局阈

2、值法最简单的阈值分割的例子就是从背景图象中提出目标物体图象。全局阈值法假设图像基本上是由感性趣物体与背景两部分组成。令位于点的象素灰度为，，选择灰度阈值为θ，则分得的二值图象:灰度值超过阈值的象素被分做物体，低于阈值的象素划分为背景。分割结果是对应物体（g值=1）与对应背景（g值=0）的二值图像。这里的1与0是逻辑值。显示分割结果时，可以根据黑白两色设置对应的灰度值（例如，在0-255灰度等级的图像中，白色的物体g值=200，黑色背景g值=0）。图(a)是未经处理的原始图像，白细胞在黑色背景上。图像(a)的直方图如图(b)所示，表现为一高一低两个峰值。如果用T=127作为阈值，分

3、割结果如图(c)。图(d)是用3x3Laplace微分算子得到的细胞边缘。实际上，在几乎所有的医学图象中，无论是物体还是背景的灰度值都有一个分布，而且有时在灰度直方图中呈双峰或多峰情况。采用加权变换后的直方图，可以得到比原直方图更明显的峰谷分界，从而改进图象的分割。在全局阈值法中，显然阈值的适当选取是一个十分重要的问题。不同的阈值会导致不同的分割结果。一般来说，我们希望阈值的微小变化不要引起分割出物体的面积或周长很大的改变。即分割结果对阈值的变化不敏感。如果用T表示阈值，A(T)和P(T)分别代表分割物体的面积与周长，如果所选择的阈值T*满足下式，称T*为最佳阈值。4.2.2大津

4、阈值分割大津(N.Otsu)阈值法可以自动寻找阈值，对图像进行划分，将目标物和背景区分开来。把直方图在某一阈值处分割成两组，当被分成的两组间方差为最大时，确定阈值。例如，设一幅图象的灰度值为1～m级，灰度值i的象素数为ni，此时我们得到：总象素数各灰度值的归一化值然后用将其分成两组和各组产生的概率如下：产生的概率产生的概率的平均值的平均值其中是整体图象的灰度平均值；是阈值为时的灰度平均值两组直方图的灰度平均值为两组间的方差用下式求出值便是阈值,不管图象的直方图有无明显的双峰，都能得到较满意的结果。4.3微分算子边缘检测边缘检测：基于物体与背景之间在灰度（或纹理）特性上存在着某种不

5、连续性（或突变性）。边缘是指它的两侧分属于两个区域，每个区域特性相对比较均匀一致，而两个区域之间在特性上则存在一定差异。左图是用一阶导数和二阶导数检测均匀灰度分布图像边缘的示例。圆形区域边缘灰度变化及导数4.3.1灰度梯度及Prewitt模板当图像微分运算后，每个象素处有一梯度值反映该点边缘强度。灰度梯度由幅值和方向两个因素构成。方向：幅值：对于数字化离散图像，可用差分代替微分。在左图中,也可表示为：利用梯度算子也可以得到二值分割图像：这里,tH是域值。域值的选取是面向问题的。一阶差分梯度算子：医学图像中大部分梯度算子可以用卷积运算，也就是局部象素灰度加权和来实现。权系数以与图像

6、中局部邻域对应的数值数组形式列出，也称作模板、窗口或卷积核。当对图像采用边缘检测模板卷积运算后，该值也反映此像点的边缘强度。边缘强度：这里，p是选用模板的序号。f(x,y)是象素(x,y)处的灰度值,h(k,l)与该象素位移(k,l)处的卷积模板元素值。对3×3模板k=l=-1,0,1m=n=1，Prewitt模板。垂直方向模板g1：-101-101-101if则在(x,y)处垂直方向有边缘通过。同理，111000-1-1-1if水平方向有边缘通过。如不考虑方向，则if则有边缘通过。水平方向模板g2：4.3.2Roberts交叉算子普通灰度梯度算子检测图像象素灰度在水平方向或垂直

7、方向上的变化。因此，这种算子对水平方向或垂直方向上的边缘较为敏感。与普通灰度梯度算子不同，Roberts交叉算子检测的是沿与图像坐标轴45度角或135度角方向上的灰度梯度：或对3×3模板k=l=-1,0,1m=n=1，Sobel模板。垂直方向模板g1：-101-202-101if则在(x,y)处垂直方向有边缘通过。同理，121000-1-2-1if水平方向有边缘通过。如不考虑方向，则if则有边缘通过。水平方向模板g2：4.3.3Sobel模板上图是用Sobel模板做边缘检测的实例