一种腹部ct序列图像肝脏的快速鲁棒自动分割方法
【技术领域】
[0001] 本发明属于图像处理技术领域，涉及腹部CT序列图像器官的分割，特别涉及腹部 CT序列图像肝脏的快速鲁棒自动分割，可用于医学影像辅助诊断与治疗。
【背景技术】
[0002] 肝脏分割是实现肝脏疾病计算机辅助诊断和肝脏移植术前规划的前提。利用分割 并重建得到的肝脏模型可以辅助肝脏病灶分析、体积测量、血管分析、肝脏分段、疾病诊断 与评估等工作。CT血管造影图像由于分辨率高，对人体损伤小，能够形象、准确地反映肝脏 及其病变位置而受到医生的普遍青睐。由于三维成像使用的图像切片数量很大（每位病人 的肝脏CT序列约为200张左右），人工分割每个切片很耗时且分割结果有很大的主观性。 因此，研究腹部CT序列图像中肝脏的快速鲁棒自动分割方法对肝脏疾病的诊断与治疗，对 提高计算机辅助诊断的精度与效率具有重要意义。
[0003] 由于肝脏器官结构复杂、形状不规则、不同个体差异较大，以及成像时受到噪音、 偏移和组织运动等的影响，得到的肝脏CT序列图像通常具有复杂性和多样性等特点。此 外，CT图像中肝脏与周围的腹部肌肉、胃、膈肌、脾脏、肾脏和心脏等器官缺乏较好的灰度对 比，这都将对CT序列图像中肝脏的准确自动分割带来很大困难。
[0004] 现有的肝脏CT序列分割方法一般可分为基于图像和统计模型两大类。单纯基于 图像的分割方法是指直接运用亮度、纹理和其他图像自身信息进行分割的方法，主要包括 阈值法、聚类、区域生长、活动轮廓模型和图割等，这些分割技术主要存在以下缺点和不足： (1)需要进行复杂的预处理，包括去除肋骨、脊椎、肾脏等周围组织或器官；（2)难以分割对 比度低、肝脏边界模糊的CT序列图像。基于统计模型的方法首先运用大量的CT序列图像 构建目标先验模型，然后将其应用于当前序列的分割，该类方法对于对比度较低的图像有 较好的分割效果，但是对于形状不规则肝脏分割效果较差，且耗时长、对数据初始化和配准 敏感。
【发明内容】
[0005] 本发明充分考虑了上述现有技术的缺点与不足，其目的在于，提供一种精确、快 速、鲁棒的腹部CT序列图像肝脏自动分割方法。并且本发明的腹部CT序列图像肝脏分割 方法能够推广到其他腹部器官的分割中。
[0006] 本发明的腹部CT序列图像肝脏鲁棒自动分割方法，包括以下步骤：
[0007] 根据肝脏区域的亮度概率分布特征，建立基于高斯曲线拟合的亮度模型，抑制复 杂背景、突出肝脏区域；
[0008] 利用图像局部信息，建立基于PCA的肝脏外观模型，进一步增强肝脏弱边界、加大 肝脏与背景的区分度；
[0009] 利用相邻切片的空间相关性，结合亮度和外观模型，构建图割能量函数，实现肝脏 区域的快速鲁棒自动分割。
[0010] 所述亮度模型建立包括：
[0011] a.从输入序列中手动选取部分肝脏区域，采用高斯函数拟合其亮度概率分布，并 根据高斯分布的置信区间获取该序列图像肝脏的大致亮度范围[1_，1_]。
[0012] b.针对肝脏亮度范围[1_，1_]，构建如下亮度模型：
[0013]
[0014] 其中，k为0. 5~5的正常数，I为图像亮度，Ln和1_分别为步骤a获取的肝脏 亮度的最小值和最大值。I越接近亮度范围的中心值，亮度模型f intOTSlty(I)的取 值越大，表明该像素点属于肝脏的概率越大。
[0015] 所述外观模型建立包括：
[0016] 以每个像素点为中心取大小为（2n+l) X (2n+l)的子图像，并以该子图像的亮度 信息表示该中心像素点的特征，其中η为1~8的自然数；
[0017] 运用PCA对用户指定的部分肝脏区域所有像素点的邻域子图像进行主成 分分析并计算其前N个主成分的均值，得到肝脏的外观特征，其中N为1~ (2 n+l) X (2n+l)的自然数；
[0018] 采用距离函数作为外观相似性度量，计算整个序列中所有像素点与肝脏的外观特 性差异：
[0019] d(p) = I |PpCA(p)-PpCA_average
[0020] 其中，pPeA(P)表示对像素点P的邻域子图像进行PCA特征提取得到的前N个主成 分，d(p)表示主成分特征P pm(P)与ppeA_avCTagf；之间的距离；
[0021] 构建基于PCA的肝脏外观模型：
[0022]
[0023] 其中，mean_dp表示对所有像素点p计算得到的距离d(p)的均值。f PeA(p)的取值 范围为0~1，fpeA(p)越大表示该像素点P属于肝脏的概率越大。
[0024] 本发明的腹部CT序列图像肝脏鲁棒自动分割方法特征还在于，所述肝脏区域的 快速鲁棒自动分割包括：
[0025] 用户指定初始切片。序列中满足肝脏区域相对较大且只包含一个肝脏连通区域的 切片图像均可被指定为初始切片，默认位于输入序列从上至下的三分之一处，例如，若输入 CT序列总数为210片，则初始切片默认为第70片，且若CT序列的切片数目不能被3整除， 则对相除结果进行四舍五入取整。
[0026] 初始切片分割。首先运用该切片的亮度和外观模型构建图割能量函数：
[0027]
[0028] 其中，P表示图像f中的所有像素集；Np表示像素点p的邻域像素集；F intanslty(fp) 和Fptt(fp)分别为亮度惩罚项和外观惩罚项，分别由亮度和外观模型计算得到；α和β分 别为控制亮度和外观惩罚项的权重，取值范围均为〇~1，且满足α +β = I ;B(fp, fq)为边 界惩罚项，由相邻像素之间的梯度计算得到，用来控制分割曲线的平滑度。然后采用最优化 算法最小化该能量函数，得到分割结果。最后取其最大连通区域作为该序列初始切片的肝 脏最终分割结果。
[0029] 序列切片分割。本发明采用迭代的方式以初始切片为起点分别向上和向下分割序 列中其他所有切片。由于相邻切片间肝脏的大小和位置不会有显著变化，因此，在迭代分割 过程中，相邻切片间的肝脏位置信息也被采用作为能量惩罚项增加分割的准确性。此时，能 量函数可表示为：
[0030]
[0031] 其中，Fwtlcin表示由上一片分割结果计算得到的位置惩罚项。本发明通过利用相 邻切片的肝脏位置信息，可以有效去除与肝脏不相关的相似组织，如脾脏、胰腺、肌肉等，得 到更准确的分割结果。
[0032] 与现有技术相比，本发明方法有以下几个优点：
[0033] 本发明提供了一种自适应的肝脏序列分割方法，针对每一位病人的数据特征建立 相应的亮度和外观模型并结合相邻切片之间的空间相关性进行自动分割，可以有效处理形 状、亮度各异的肝脏，具有较强的鲁棒性；
[0034] 不需要复杂的预处理过程，不需要对肌肉、肋骨、脊椎、肾脏等进行预先分割；
[0035] 不需要进行繁琐的训练与统计模型构建，分割速度较快；
[0036] 本发明的腹部CT序列图像肝脏分割方法能够推广到腹部CT序列图像其他器官的 分割应用中，如脾脏和肾脏的分割。
【附图说明】
[0037] 图1本发明实施方式的腹部CT序列图像肝脏鲁棒自动分割方法流程图；
[0038] 图2本发明实施方式用户指定的部分肝脏区域示意图；
[0039] 图3本发明实施方式的用户指定肝脏区域高斯曲线拟合结果图；
[0040] 图4本发明实施方式的的肝脏亮度模型示例图；
[0041] 图5本发明实施方式的的肝脏外观模型示例图；
[0042] 图6本发明实施方式的肝脏分割结果示例图。
【具体实施方式】
[0043] 实施例1
[0044] 图1所示为本发明实施方式的腹部CT序列图像肝脏鲁棒自动分割方法流程图。 首先从输入的CT序列中任意选取部分肝脏区域分别建立肝脏亮度模型和外观模型，然后 结合亮度和外观模型运用图割算法对初始切片中的肝脏进行分割，最后采用迭代的方式以 初始分割切片为起点分别向上、向下分割序列中其他所有切片的肝脏。在迭代分割过程中， 将上一片分割结果的肝脏位置信息融入当前切片的图割能量函数，以增加分割结果的准确 性。直到所有切片分割完成，程序运行结束，输出分割结果。
[0045] 下面结合图1，以一较佳的实施例详细说明本发明的腹部CT序列图像肝脏鲁棒自 动分割方法。
[0046] 1.肝脏亮度模型建立。具体实现步骤如下：
[0047] 根据CT扫描特性，某个特定器官的亮度通常大致位于一个较窄的区间范围，且在 该区间上呈高斯分布。设腹部CT序列中肝脏亮度范围为，则肝脏亮度模型可构 建如下：
[0048]
[0049] 其中，k为0. 5~5的正常数，本实施例优选为1. 5, I为图像亮度。I越接近亮度 范围的中心值，亮度模型fintenslty(I)的取值越大，表明该像素点属于肝脏的概率 越大。
[0050] 为了获取亮度范围[1_，1_]，本发明对用户指定的部分肝脏区域的亮度分布进 行高斯拟合：
[0051]
[0052] 其中，a表示高斯分布的中心位置，b控制高斯分布的宽度，c为高斯分布的峰值。 图2为原始CT序列中的任一切片图像，其中矩形区域表示用户指定的部分肝脏区域。图 3所示为利用最小二乘算法对该肝脏区域进行高斯曲线拟合的结