基于自适应表面形变模型的ct图像肝脏分割方法
【专利摘要】本发明提出了一种基于自适应表面形变模型的CT图像肝脏分割方法，为肝脏疾病的临床诊断提供了一种便利工具。首先利用各向异性滤波的方法对图像进行预处理，得到初始边界图像；第二步，利用球形简单形变模型(Deformable Simplex Model,DSM)描述肝脏的初始轮廓；第三步根据模型顶点及与其邻域之间的关系计算模型的内力，由原始图像的梯度和边界来计算模型的外力；第四步构建模型的内力、外力的约束模型；第五步三角网格的自适应分解模型；第六步设定迭代次数并在内力、外力和气球力的驱动下逼近目标区域；最后得到精确的分割结果。
【专利说明】基于自适应表面形变模型的CT图像肝脏分割方法
【技术领域】
[0001] 本发明涉及一种基于自适应表面形变模型的CT图像肝脏分割方法，属于肝脏分 割的图像扫描领域。
【背景技术】
[0002] 肝脏是人体内最大的代谢器官，人体所需物质和营养素的合成、分解、转化和存储 中心。而肝脏分割是计算机辅助诊疗、肿瘤切除、手术规划和活体移植的重要基础。肝脏分 割的结果直接影响到肝脏三维重建的精确度，影响医生对肝脏中病灶的诊断和治疗。所以 对肝脏分割方法的研究不仅具有重要的学术意义且具有重要的临床应用价值。
[0003] 然而，在计算机断层扫描（CT)图像中，肝脏组织的灰度值往往与其邻近的心脏、 胃等组织相近，使简单的阈值分割和区域增长等方法无法实现肝脏的精确分割，且大多 数临床图像中都承载着肝脏组织的严重疾病，如肝脏肿瘤、肝硬化、或部分肝脏切除后的 伤疤等信息。所有这些特殊情况的准确诊断和治疗都必须通过提高肝脏分割算法来实 现。近年来，许多研究学者尝试不同领域的分割算法应用到肝脏的分割应用中，其主流算 法是基于模型的分割算法，主要包括可形变模型（DeformableModels，DM)、概率图谱模型 (ProbabilisticAtlasModels,PAM)和统计形状模型（StatisticalShapeModels,SSM) 等。然而，由于肝脏组织形状的多变性和复杂性，肝脏分割的精确性问题仍然未能得到很好 地解决。当今主流的基于模型的分割算法主要存在以下方面的局限性：
[0004] 1.形变模型的初始位置严重影响最后的分割结果，假如模型初始化太小，则容易 产生模型达到局部平衡或缩小为一个点的情况。此时，不能有效分割出肝脏组织。
[0005] 2.假如模型初始化太大时，则容易使形变模型溢出分割区域，产生过分割。
[0006] 3.由于形变模型平滑力的影响，致使其不能有效的进入肝脏的细长尖角区域，严 重影响分割精度。
【发明内容】
[0007] 针对上述问题，本发明提出了一种基于自适应表面形变模型的CT图像肝脏分割 方法，能够分割出准确的肝脏区域，避免了人为干扰，显著提高了分割效率。
[0008] 该基于自适应表面形变模型的CT图像肝脏分割方法，主要包括以下步骤：
[0009] 第一步：利用各向异性扩散滤波对原始图像进行预处理，并得到初始边界图像；
[0010] 第二步：模型初始位姿设定：输入肝脏的原始图像，并初始化成基于简单形变模 型的球形描述模型，用于表示肝脏组织的初始表面；
[0011]第三步：模型内、外力计算：获得原始CT图像的梯度图并利用梯度计算模型的外 力，并利用初始化的球形描述模型内部关系获得形变模型的内力；
[0012] 第四步：构建内、外力的约束模型：利用原始图像的灰度直方图获得包含有分割 目标区域的二值图像，并利用此二值图像约束模型的外力；
[0013] 第五步：三角网格的自适应分解：在模型形变过程中，计算三角网格的面积，如果 大于某一个阈值，则在三角网格的重心插入新的模型顶点，将原始的三角网格分解为三个 面积较小的三角网格；
[0014] 第六步：模型的迭代形变：设定迭代次数，进行模型的迭代形变，使其有效逼近目 标区域，并得到最终分割结果。
[0015] 本发明的有益效果：
[0016] 1.利用DSM的球形描述模型表示肝脏的初始表面。
[0017] 2.提出了模型内力的约束方法，使其能够分割肝脏的细长尖角区域。
[0018] 3.构建了模型外力的约束方法，使其能够限制模型分割算法的过分割问题。
[0019] 4.提出了基于模型顶点法向量的气球力计算方法，并将其添加到模型的形变过程 中，使模型能够避免达到局部平衡或缩小为一个点的缺陷。
[0020] 5.提出了模型在形变过程中，三角网格可以自适应的分解，从而保证分割结果的 平滑性和准确性。
【专利附图】
【附图说明】
[0021] 图1是本发明所提出的工作流程图；
[0022] 图2是本发明所提出的内力计算方法示意图；
[0023] 图3是本发明所提出的自适应三角网格分解的示意图。
【具体实施方式】
[0024] 如附图1所示，本发明的流程图，具体包括以下几个步骤：
[0025] 步骤S101，图像预处理。
[0026] 利用各向异性扩散滤波的方法对原始图像进行预处理，得到与原始图像相对应的 边界图像。
[0027] 步骤S102,模型的初始化。
[0028] 通过输入的原始肝脏CT图像，确定肝脏组织的大概区域，并计算此区域的重心。 以此重心为球心，以重心到区域边界的最短距离为半径初始化一个基于DSM描述的球形模 型。以此作为肝脏组织的初始表面。
[0029] 步骤S103,模型内、外力的计算。
[0030] 如附图2所示，模型的内力分为切向力Ftangmt和法向力Fnmial。其中切向力的计算 公式如下：
[0031]
【权利要求】
1. 基于自适应表面形变模型的CT图像肝脏分割方法，其特征在于，主要包括以下步 骤： 第一步：利用各向异性扩散滤波对原始图像进行预处理，并得到初始边界图像； 第二步：模型初始位姿设定：输入肝脏的原始图像，并初始化成基于简单形变模型的 球形描述模型，用于表示肝脏组织的初始表面； 第三步：模型内、外力计算：获得原始CT图像的梯度图并利用梯度计算模型的外力，并 利用初始化的球形描述模型内部关系获得形变模型的内力； 第四步：构建内外力的约束模型：利用原始图像的灰度直方图获得包含有分割目标区 域的二值图像，并利用此二值图像约束模型的外力； 第五步：三角网格的自适应分解：在模型形变过程中，计算三角网格的面积，如果大于 某一个阈值，则在三角网格的重心插入新的模型顶点，将原始的三角网格分解为三个面积 较小的三角网格； 第六步：模型的迭代形变：设定迭代次数，进行模型的迭代形变，使其有效逼近目标区 域，并得到最终分割结果。
2. 如权利要求1所述的基于自适应表面形变模型的CT图像肝脏分割方法，其特征在 于，所述第三步中，根据图像的梯度和边界信息计算模型的梯度力Fgradimt和边界力Fedge,从 而构建模型的外力Fraitenal Fgradi ent+Fedge 0
3. 如权利要求2所述的基于自适应表面形变模型的CT图像肝脏分割方法，其特征在 于，所述第三步中，模型的梯度力Fgradimt的计算公式如下： =βΓ\(β-Υι)·η?)ηι 其中，;巧^是一个在[0, 1]内的权重系数，Pi是模型中的一个顶点，Gi是包含顶点Pi的mXmXm区域的中心点，€是？1点处的法线方向。
4. 如权利要求3所述的基于自适应表面形变模型的CT图像肝脏分割方法，其特征在 于，其中权重系数使模型某点的梯度力Fgradimt随着此点与二值图像Ibtoy (X，y，z)的边 界距离增大而增大，计算公式如下：
其中，K是一个大于0的常量，f为根据原始图像的二值图像Ibinmy (X，y，z)得到的边界 图，其计算公式如下：
其中，"▽"是梯度算子，是方差为σ的高斯核函数。
5. 如权利要求3所述的基于自适应表面形变模型的CT图像肝脏分割方法，其特征在 于，所述第三步中，模型的边界力^_的计算公式如下： Fedge=^nm 其中，成#是一个在[〇, 1]内的权重系数，Pi是模型中的一个顶点，Ei-Pi是与法向量?共线的向量。
6. 如权利要求5所述的基于自适应表面形变模型的CT图像肝脏分割方法，其特征在 于，权重系数及〃#使模型某点的边界力FedgJ逭着此点与二值图像户_y,z)的边界距离 增大而减小，计算公式如下：
其中，K是一个大于0的常量，f为根据原始图像的二值图像Ibinmy (X，y，z)得到的边界 图，其计算公式如下：
其中，"V"是梯度算子，是方差为σ的高斯核函数。
7. 如权利要求1所述的基于自适应表面形变模型的CT图像肝脏分割方法，其特征在 于，所述第三步中，根据模型某顶点与其邻域的位置关系计算模型的切向力Ftangmt和法向 力Fnarmai，从而构建t旲型的内力Fint ernal - Ftan gent+Fnormal °
【文档编号】G06T7/00GK104318553SQ201410543371
【公开日】2015年1月28日 申请日期:2014年10月15日 优先权日:2014年10月15日
【发明者】杨健, 王雪虎, 王涌天, 刘越, 艾丹妮 申请人:北京理工大学