医学图像分割方法综述

医学图像分割方法综述

林瑶，田捷1

北京,中国科学院自动化研究所人工智能实验室,100080

摘要: 图像分割是一个经典难题，随着影像医学的发展，图像分割在医学应用中具有特殊的重要意义。本文从医学应用的角度出发，对医学图像分割方法，特别是近几年来图像分割领域中出现的新思路、新方法或对原有方法的新的改进给出了一个比较全面的综述，最后总结了医学图像分割方法的研究特点。

关键词：医学图像分割 综述

1．背景介绍

医学图像包括CT 、正电子放射层析成像技术（PET ）、单光子辐射断层摄像（SPECT ）、MRI （磁共振成像技术）、Ultrasound （超声）及其它医学影像设备所获得的图像。随着影像医学在临床医学的成功应用，图像分割在影像医学中发挥着越来越大的作用[1]。图像分割是提取影像图像中特殊组织的定量信息的不可缺少的手段，同时也是可视化实现的预处理步骤和前提。分割后的图像正被广泛应用于各种场合，如组织容积的定量分析，诊断，病变组织的定位，解剖结构的学习，治疗规划，功能成像数据的局部体效应校正和计算机指导手术[2]。 所谓图像分割是指将图像中具有特殊涵义的不同区域区分开来，这些区域是互相不交叉的，每一个区域都满足特定区域的一致性。

定义 将一幅图像，其中g x y (,)0≤≤x Max x _,0≤≤y Max y _，进行分割就是将图像划分为满足如下条件的子区域...：

g 1g 2g 3 (a) ，即所有子区域组成了整幅图像。

(b) 是连通的区域。

g k (c) ，即任意两个子区域不存在公共元素。 (d) 区域满足一定的均一性条件。均一性（或相似性）一般指同一区域内的像素点之间的灰度值差异较小或灰度值的变化较缓慢。

g k 如果连通性的约束被取消，那么对像素集的划分就称为分类（pixel classification），每一个像素集称为类（class）。在下面的叙述中，为了简单，我们将经典的分割和像素分类通称为分割。

医学图像分割到今天仍然没有获得解决，一个重要的原因是医学图像的复杂性和多样性。由于医学图像的成像原理和组织本身的特性差异，图像的形成受到诸如噪音、场偏移效应、局部体效应和组织运动等的影响，医学图像与普通图像比较，不可避免的具有模糊、不均匀性等特点。另外，人体的解剖组织结构和形状复杂，而且人与人之间有相当大的差别。这些都给医学图像分割的分割带来了困难。传统的分割技术或者完全失败，或者需要一些特殊的处理技术。因此，我们有必要针对医学应用这个领域，对图像分割方法进行研究。

为了解决医学图像的分割问题，近几年来，很多研究人员做了大量的工作，提出了很多实用的分割算法[2][3][4]，随着统计学理论、模糊集理论、神经网络、形态学理论、小波理论等在图像分割中的应用日渐广泛，遗传算法、尺度空间、多分辨率方法、非线性扩散方程等近期涌现的新方法和新思想也不断被用于解决分割问题，国内外学者提出了不少有针对性的好分割方法。本文将主要介绍近几年这一领域中研究人员提出的新方法或对原有方法的新改进。需要指出的是，由于从不同的角度将得到不同的分类结果，本文中所涉及方法的分类并不是绝对的，而且许多分割方法还是多种简单方法的综合体，我们只能大致将它们分为属于最能反映其特点 1x x g N k k =),(),(y g y =∪φ=(y y g j k ∩),(),x

g x 1 联系人：田捷 电话：82618465 E-mail:tian@

的某一类。

2．基于区域的分割方法

图像分割通常会用到不同对象间特征的不连续性和同一对象内部的特征相似性。基于区域的算2.1 法

是最常见的并行的直接检测区域的分割方法[5]。如果只用选取一个阈值称为单阈值分同类的物体灰度值或其他特征值相差很大时，它能很有效大的图像，

对于图像中不存在明显的灰的困难所在。至今仍有.2 区域生长和分裂合并

两种典型的串行区域分割方法。其特点是将分割过程分解为顺序的多个，该方法需要先选取一个种子用者必须在每个需要抽取出的区域中法测重于利用区域内特征的相似性。

阈值阈值分割割，它将图像分为目标和背景两大类；如果用多个阈值分割称为多阈值方法，图像将被分割为多个目标区域和背景，为区分目标，还需要对各个区域进行标记。阈值分割方法基于对灰度图像的一种假设：目标或背景内的相邻像素间的灰度值是相似的，但不同目标或背景的像素在灰度上有差异，反映在图像直方图上，不同目标和背景则对应不同的峰。选取的阈值应位于两个峰之间的谷，从而将各个峰分开。

阈值分割的优点是实现简单，对于不的对图像进行分割。阈值分割通常作为医学图像的预处理，然后应用其他一系列分割方法进行后处理。它也常被用于CT 图像中皮肤、骨骼的分割。

阈值分割的缺点是不适用于多通道图像和特征值相差不度差异或各物体的灰度值范围有较大重叠的图像分割问题难以得到准确的结果。另外，由于它仅仅考虑了图像的灰度信息而不考虑图像的空间信息，阈值分割对噪声和灰度不均匀很敏感。针对阈值分割方法的缺点，不少学者提出了许多改进方法，如基于过渡区的方法[6]，还有利用像素点空间位置信息的变化阈值法[7]，结合连通信息[8]的阈值方法。

对于多目标的图像来讲，如何选取合适的阈值实在是基于阈值分割方法不少学者针对该问题进行深入的研究，提出了许多新方法。在近年来的自动选取阈值方法中，基于最大熵原则选择阈值是最重要的方法之一，由T.Pun 首先在[9]中提出。这种方法的目的在于将图像的灰度直方图分成两个或多个独立的类，使得各类熵的总值最大，从信息论角度来说就是使这样选择阈值获得的信息量最大。J.N.Kapur 等人进一步发展了这种方法[10]，P.Sahoo等人提出了用Renyi熵代替常规熵的最大熵原则[11]。Jui-Cheng Yen等人提出用最大相关性原则选择阈值[12]，这种方法其实只是用他们定义的一个最大相关性原则取代了一般用的最大熵原则。

2区域生长和分裂合并是步骤，其中后续步骤要根据前面步骤的结果进行判断而确定。

区域生长的基本思想是将具有相似性质的像素集合起来构成区域点，然后依次将种子像素周围的相似像素合并到种子像素所在的区域中。区域生长算法的研究重点一是特征度量和区域增长规则的设计，二是算法的高效性和准确性。区域增长方式的优点是计算简单。与阈值分割类似，区域增长也很少单独使用，往往是与其他分割方法一起使用，特别适用于分割小的结构如肿瘤和伤疤[13]。

区域生长的缺点是它需要人工交互以获得种子点，这样使植入一个种子点。同时，区域增长方式也对噪声敏感，导致抽取出的区域有空洞或者在局部体效应的情况下将分开的区域连接起来。为解决这些问题，J.F. Mangin 等提出了一种同伦的（homotopic ）区域生长方式[14]，以保证初始区域和最终抽取出的区域的拓扑结构相同。Shu-Yen Wan 等提出的对称区域增长算法[15]有效地弥补了原算法对种子点敏感和占用内存多的弱点，而且对3D 连接对象标记和删除空洞的算法效率较高。另外，模糊连接度方法与区域增长相结合也是一个发展方向[16]。