基于模糊C均值和随机漫步的CT肝脏图像分割算法研究

王琴琴

南京医科大学附属妇产医院（南京市妇幼保健院）放射科，江苏 南京 210004

引言

在计算机辅助手术、治疗计划制定等肝脏相关临床应用中，肝脏区域及肝肿瘤分割是必要的步骤[1]。目前临床上主要采用手工方式勾画肝脏目标区域，耗时且主观差异大。随着分割算法的产生，医学图像分割朝着交互式分割和全自动分割方向发展，其中全自动分割无需人为干预，适用于简单均匀图像分割，对于较复杂和灰度不均图像，分割结果往往不能满足临床需求。交互式结合人工干预与分割算法的优势，显著提高图像分割准确性。

交互式图像分割算法主要利用图像目标像素点与背景像素点的空间位置关系和灰度、边缘等特征信息来完成图像分割[2]，主要包含图切理论，分水岭算法，水平集算法和随机漫步算法。其中，随机漫步算法具有分割结果唯一、对模糊边界分割效果显著、对噪声鲁棒性佳等特点，逐步成为应用最广泛的交互式分割算法[3]。但是随机漫步算法分割结果受初始种子点影响较大，为探索更优的分割算法，本研究采用模糊C 均值算法先对图像进行初步分类，分割成相似灰度特征的目标与背景区域，然后采用随机漫步算法提取目标边界，完成图像分割，期望能减少对初始种子点的影响，并有效分割灰度不均和肿瘤边界模糊图像。

1 方法

1.1 图像预处理

图像预处理采用3×3 中值滤波器，不产生新像素灰度，且能保留图像边缘信息。

1.2 标准FCM算法

FCM 算法主要思想是将N 个L 维向量分为C 个模糊组，通过距离测度判断像素灰度属于各类别的隶属程度，不断更新图像像素隶属度和聚类中心，使目标函数达到最小，从而完成像素分类和图像分割[4]。其中像素隶属度是用来描述像素属于某个类别的程度，取值范围为[0,1]。FCM 算法的目标函数和约束条件见公式(1)[5]，uij=ui(xj)表示像素灰度xj属于类别i的程度，m是模糊加权系数，常取2，νi={ν1,…,νC}表示i个聚类中心，dij (xj, νi)=||xj-νi||2表示第j像素灰度到第i个聚类中心的欧式距离。

采用朗格朗日乘数法求取目标函数最小值，见公式(2)，分别对参数λ和隶属度uij求导可得隶属度和聚类中心，如公式(3)和(4)，直到|u(k+1)-uk|<∈（设定的阈值）停止迭代，完成聚类优化过程。

标准FCM 算法对简单图像分割效果显著，但不能分割对比度差、区域重叠和异质性显著的图像。主要由于FCM 算法对初始聚类中心敏感，容易陷入局部最优解。其次，并未关注图像像素空间的联结性与局部特征，受噪声影响显著。本研究先采用标准FCM 算法进行初步分类，然后采用随机漫步算法进一步提升图像分割准确度。

1.3 随机漫步算法

随机漫步算法是一种有监督算法，通过计算从每一个像素点出发的随机漫步者首次到达各类种子点的概率完成像素分类[6]。首先，将图像看作具有固定个数顶点和边组成的离散无向图G=，图像像素对应图的结点νi∈V，每条边e∈E⊆V×V连接两个端点νi和νj，表示为eij，并被赋予一定的权值wij，反映相邻像素之间的差异或相似度，与两个像素点在调整空间中的绝对距离呈负相关，权函数由公式(5)得出[7]。其中Ii表示节点νi的灰度值，Ij表示与结点νi相邻的结点νj的灰度值，β是自由度参数。

由电路理论启发，Grady 提出一个随机漫步者首次到达一个种子点的概率等于种子点边界条件下的狄利克雷问题的解[8]，狄利克雷问题求解过程如下：首先定义一个拉普拉斯矩阵Lij，当νi=νj时，Lij=di；当结点νi与νj相邻时，Lij=wij；其他情况，Lij=0，其中di=∑wij称为结点的度，为与结点相邻的所有结点的权值之和。

离散的狄利克雷积分公式见式(6)，此时需要将顶点G分成VM（已标记的种子结点）和VN（未标记的种子结点）两个集合，VM和VN满足VM∪VN=V，VM∩VN=∅，根据结点所属集合不同将拉普拉斯矩阵分解为对角阵，见公式(7)。

公式(6)根据变换后的拉普拉斯矩阵可改写如下形式，见公式(8)。其中xM和xN分别对应种子结点和非种子结点的电势值，即随机漫步算法的概率值。通过对D[xN]求关于xN的微分可得临界点满足公式(9)。

假设结点νi到标记为s的种子点的概率为xis，将各种子点的标记集合定义为函数Q(νj )=s,∀νj∈VM，其中0

1.4 图像分割评价指标

本研究采用定性与定量相结合评估分割图像质量，定性分析采用视觉评估，定量评价指标采用重叠错误率（Overlap Error，OE）[9]、相对误差（Relative Difference，RD）[10]、Dice 相似性系数（Dice Similarity Coefficient，DSC）[11]，见公式(13)～(15)。其中A代表分割所得区域，B代表真实的分割区域。OE 和RD 越小，DSC 越大，表明图像分割效果越佳。

2 方法

将本研究提出算法与文献报道的分割算法进行性能比较。选用3Dircadb 数据库和MIDAS 数据库图像进行仿真实验，纳入大小超过5 mm 的肝肿瘤，排除小于5 mm 的肝肿瘤。所有算法测试均在Matlab 2013a 编程环境下仿真实现。

2.1 本研究提出算法的图像分割

图1 表示不同权值wij对本研究图像分割结果的影响，可以看出β变小，wij变小，目标区域边界逐渐扩大，向真实边界靠近。图2 表示肿瘤图像分割，基于FCM 算法可完成粗布的肿瘤提取，随机漫步算法可进一步精确勾画肿瘤边界，可以看出技师勾画肿瘤边界较大，传统随机漫步算法对灰度不均的肿瘤分割精度显然低于本研究算法，对于肿瘤边界显著图像，分割效果相当。图3 表示不同尺寸肿瘤图像分割结果，可以看出对于相对尺寸较大的肿瘤，本研究算法能精确勾画肿瘤边界，但是对于小尺寸肿瘤，还有待改进。

图1 不同权值wij对图像分割结果的影响

2.2 不同算法肝脏图像分割结果比较

不同分割算法性能比较结果，见表1。由此可推断本研究分割算法所得评价指标优于其他文献报道算法，可获得更高的分割精度。在对MIDAS 数据库图像（包含10 个肿瘤）分割时，基于本研究算法可获得平均0.81%的DSC、15.61%的OE 和4.02%的RD。同时在对3Dircadb 数据库图像（包含117 个肿瘤）进行分割时，基于本研究算法所得的平均DSC 为0.75%，OE 为22.78%，RD 为8.59%，分割精度较MIDAS 数据库有所下降，主要由于后者包含众多小尺寸肿瘤病灶，影响分割精度。此外发现操作者手工分割与技师分割时均存在过度分割现象，即肿瘤边界大于实际轮廓。

图2 剪切后包含22 mm和7 mm肿瘤的CT肝脏图像

3 结论

本研究提联合使用FCM 算法和随机漫步算法对CT肝脏图像进行分割，其中FCM 算法可将图像初步划分为不同类的目标区域和背景区域，降低各区域内种子点灰度差异。然后随机漫步算法可快速收敛，完成各像素分类，提取肿瘤边界。仿真实验结果表明，基于本研究提出的CT 肝脏图像分割算法优于多篇文献报道的其他分割算法，能有效分割灰度不均和肿瘤边界模糊图像，是一种可行的图像分割算法。

表1 不同图像分割算法所得结果比较

图3 基于本研究算法不同尺寸肝脏肿瘤图像分割结果（红线）与真实分割结果（绿线）比较