图像分割的目的就是把目标从背景中提取出来, 分割过程主要基于图像的固有特征, 如灰度、纹理、对比度、亮度、彩色特征等将图像分成具有各自特性的同质区域。
医学成像技术包括核磁共振成像（MRI）、计算机断层扫描（CT）、超声、正电子发射断层扫描（PET）等。
医学图像分割方法和理论众多，包括基于边界、阈值、区域增长、统计学、图论、活动轮廓、信息论、模糊集理论、神经网络等。下面进行详细介绍。
1. 基于边界的医学图像分割方法
将图像划分为子图, 获得二维图像时, 不可避免会存在边界, 边界像素表现为灰度幅值的不连续。而显现边界的关键是边界两侧像素特征存在明显的差异, 因此可以通过边界进行子图像划分, 实现图像的分割。
过去, 科研医师与临床工作者采用鼠标以及光笔等工具确定边界, 这种方法效率低、可重复性差、对主观依赖性强。而利用计算机, 通过特定的数学模型选择合适边界点, 可极大的提高边界分割效率, 临床上常通过确定目标物体的边缘像素，利用一阶导数的级数或二阶导数的过零点信息提取边缘特征。
胸华晖等采集健康男性层间距为4 mm冠状面头颈部MRI图像, 采用距离正则化水平集图像分割方法分割鼻腔边界, 采用插足算法可计算任意位置的鼻腔横截面积, 可方便的进行鼻腔结构三维重建。
2. 基于阈值的医学图像分割方法
阈值是一种经典的基于区域的并行分割技术, 其将图像像素点分为若干类, 在肿瘤性疾病诊断中具有较高的应用价值。
Taxt等将基于局部或全局信息的阈值设置技术定义为基于上下文和非上下文分割方法, 按照分割区域方式又可分为局部阈值分割, 或称为自适应性阈值法。阈值分割利用图像灰度直方图信息, 计算量小, 利于实现, 但是易导致错误分割, 同时未考虑空间信息, 易受噪音等原因产生伪影, 其主要作为一种预处理方法。
A.Varga-Szemes等试构建了一种基于心肌信号强度阈值的半自动心脏MRI分割法, 并与常规基于轮廓的图像分割方法进行比较, 结果显示阈值法耗时更短, 其对EDV、ESW、SV、EF显示低于常规方法, 阈值法与主动脉血流测量的SV具有较好的一致性。
3. 基于区域增长分割方法
区域生长和分裂合并方法, 是一种典型的串行区域分割方法, 其根据用户的预定义相似性函数, 将图像像素或者子区域聚合成为更大的区域。在实际应用中, 区域增长法需要确定区域数目, 选择具有代表性的种子点, 选择有意义的特征, 选择停止的准则, 计算方法简单, 区域算法的精髓是将图像分类为多个相似的小区域, 按照邻接参考某种标准合并, 对于连续的均匀的目标, 如血管、实质性组织, 具有较好的准确性、高效性。该法对噪音较敏感, 更适合分割肿瘤、伤疤这些较小的结构。
区域生长是一种基于传统的区域生长机制、利用区域内像素的相似性进行图像分割的图像处理算法。区域生长的算法思想是, 对于给定的区域, 我们选取一个或者多个像素作为种子点, 以4邻域或8邻域的方式遍历种子点周围的像素点, 把那些与种子点相似的像素合并到该种子点所在的区域, 接着把那些新加入种子点所在区域的像素作为新的种子点继续进行上述过程, 直到所有相似的像素都合并到种子区域为止, 这样区域生长的过程就完成了并停止生长。
区域生长的特点就是初始化简单、执行速度快, 执行过程中可以人工参与调节, 终止易控制, 适用于灰度均匀的目标区域提取。但是该算法对噪声敏感, 需要人工指定种子点, 并且对于复杂图像很可能产生过分割和欠分割。
国外提出许多自适应区域生长算法, 该方法可自动优化学习分割区域形状特征同质性准则, 鲁棒性较好。代双凤等基于3D区域增长法进行全肺分割, 将连通域标记法、形态学方法相结合去掉气管和主支气管, 得到初步的肺实质掩膜, 最后采用凸包算法修补平滑, 修补后的精确度较高。
4. 基于统计学的医学图像分割方法
从统计学的角度进行图像分割是指将图像的各个像素点特征值作为一种概率分布的随机变量, 从概率的角度来看, 找出最应该得到图像的某种组合, 又可分为分类器、聚类、随机场、混合分布等。以最常用的聚类分析为例, 其余分类器算法类似, 但无需训练样本, 是一种无监督的算法, 通过迭代进行图像分类, 提取每一类的特征值, 按照算法又可分为K均值、模糊C均值、期望最大化算法等。
周显国等基于改进模糊聚类分析, 对医学脑部MRI图像进行分割, 利用图像的直方图作为模糊聚类的初始聚类中心, 减少迭代次数, 仿真算法证实其相较于标准的FCM算法以及其他算法, 分割效果显著提高。利用统计学方法图像分割方法较多, 各有优劣。
5.基于图论的医学图像分割方法
基于图论的图像分割是一种新的分割技术, 无需初始化, 其可将图像分割转换为图的最优划分问题, 利用剪切标准得到最佳的分割策略, 是一种全局准则, 图像分割的过程中是得到全局最优解的过程。
Shi等提出归一化割原则, Ng提出了NIW算法, 后者是当前常用的谱分割算法之一。崔宝霞等采用最小生成树法对人肺部CT图像进行分割, 而后采用MC算法进行三维重建, 实现了肺部的三维立体显示, 实验仿真证实其快速有效, 可明显提高肺部CT图像三维重建效率、完整度, 但该方法因为需要处理的数据较大, 运行速度较慢, 仍需要进行算法改进, 以提高运行效率。
6. 基于活动轮廓模型的图像分割方法
基于边缘的活动轮廓莫悉尼港采用边缘停止函数,可分割边缘化较明显的图像, 对噪声较敏感, 其主要优势在于无需设置图像的全局约束项, 对于某些分割目标和背景不均匀的图像分割效果较好。
最近Zhu等提出了一种基于活动轮廓模型的区域竞争分割方法, 通过构建基于目标与区域背景的拟合图像函数, 可有效同时检测目标和目标边界。赵立川等基于改进自仿射映射系统与参数活动轮廓的医学图像分割法, 采用高斯滤波器对图像进行平滑处理, 计算小波系数, 而后在每个小波尺度中定义一些自仿射映射, 将小波尺度对应的子力叠加获得自仿射能力, 采用动态力公式引导Snake模型, 仿真实验显示该算法分割性能好, 可提高Snake模型对弱边界、深凹陷区域的处理能力, 其鲁棒性较好, 对于100%的噪声级, 算法拟合误差度至少为40%。
7. 基于信息论图像分割方法
信息论是旨在采用概率论、数理统计方法并 将信息的传递视为一种概率事件, 熵是信息论最基本的度量单位, 医学图像分割便采用信息论中的熵理论进行, 其定义为
其中, P (x) 为随机变量x的概率密度函数, 处理医学影像时x表示像素特征, 如灰度、梯度等, 通过求熵的极值可获得最合适的分割阈值。目前最常用的方法为 极小极大熵值公式分割算法, 其将图像灰度直方图划分为多个相互独立的类, 使各类熵总量最大, 其具有全局客观属性, 但收敛速度慢, 计算量大, 对噪声较大、结构复杂的图像分割效果不理想。
李爱菊等采用改进布谷鸟搜索算法最大熵值进行医学图像分割, 采用最大熵值法找到图像的分割目标函数, 而后采用改进布谷鸟算法优化函数, 找到最佳的分割点, 反正模式计算其解决了传统最大熵值法缺陷, 提高了分割的精度, 具有较好的鲁棒性。
8. 基于模糊集理论图像分割方法
模糊集是指采用不准确的经验知识定义一个事件。模糊理论主要被用于模式识别领域, 用于表示信息不完备、不确定、病态表示的不确定性, 用于图像分割技术又可分为模糊阈值分割、模糊聚类分割、模糊连接度分割等。
刘会明等尝试基于模糊集理论, 进行胸片图像分割, 并与传统的区域生长法进行对比, 证实模糊C-均值聚类算法效果较好, 可准确的提取出肺部的轮廓, 无需人工干预, 避免了阈值设定问题, 在图像质量不佳的情况下, 也能够提取出感兴趣的区域, 分割性能有所提高。
9. 基于神经网络图像分割方法
神经网络通过大量的具有自适应性的简单单元组成相互连接的网络模型, 在处理医学图像时, 可作为训练分类器, 来确定权重, 其实质是一个分类以及标记图像的过程中, 按照拓扑结构又可分为前向神经网络、反馈网络、自组织神经网络。
以自组织神经网络为例, 其可通过输入的*纹理特征8共生矩阵进行学习, 客观处理目标图像的特征值。神经网络图像分割将分割问题转化为能量函数、分类等问题, 通过训练, 可修正节点之间的连接关系以及权值, 进而更好的对新的图像进行分割, 该分割方法需要大量训练样本。
李健等基于多尺度卷积神经网络尝试对MRI脑肿瘤进行分割, 通过多尺度输入、多尺度下的采样, 克服脑肿瘤个体差异, 自适应性不同图层之间肿瘤大小位置差异,对30例患者进行模拟训练结果显示, Dice系数为83.11%, 平均灵敏度为89.48%, 平均阳性预测值为78.91%, 分割精度明显提高。
参考链接：
医学图像分割方法综述-刘宇 陈胜
医学图像分割方法综述-李兰兰