基于核磁共振图像的乳腺肿瘤分割方法

基于核磁共振图像的乳腺肿瘤分割方法
【专利摘要】本发明公开了一种乳腺肿瘤分割方法，包括构造乳腺组织核磁共振图像的分类与偏移场矫正的耦合框架的步骤，对乳房区域以及周边区域进行增强的步骤，以及结合形状先验对乳房肿瘤图像进行分割的步骤。本发明将偏移场信息融入到分类模型中将二者归纳为统一框架，通过快速能量极小化方法同时求解乳腺核磁共振图像分类以及矫正偏移场，使得两者在模型演化过程中相互利用对方的信息最终实现二者的准确求解；并对血管以及肿瘤的形状进行分析，抓住二者在形状上的差异，利用特征值以及特征向量等参数建立形状先验，构造基于形状的水平集驱动力并将其与基于局部信息的水平集方法相结合，使得水平集在演化时克服管状结构的干扰只捕获乳腺肿瘤区域。
【专利说明】基于核磁共振图像的乳腺肿瘤分割方法

【技术领域】
[0001] 本发明属于图像分割【技术领域】，尤其是涉及一种基于核磁共振图像的乳腺肿瘤组 织的提取和分割方法。

【背景技术】
[0002] 乳腺癌是女性排名第一的常见恶性肿瘤，全球每年有超过10000名女性死于该疾 病。2010年，美国有近207090个新的乳腺癌患者，而到2011年该项统计数据增长到230480 例，占女性新发占女性新发恶性肿瘤的30%，排名女性恶性肿瘤发病率第一位。在我国，北 京、上海、天津等大城市的统计显示乳腺癌同样是我国女性最常见的恶性肿瘤，且发病率呈 逐年上升趋势。乳腺癌的症状多种多样，常见的有：乳腺肿块，乳腺疼痛，乳头溢液、糜烂或 皮肤凹陷，腋窝淋巴结肿大等。其中，乳腺肿瘤是最常见的症状，约90%的患者是以该症状 进行就诊的。乳腺肿瘤分为良性肿瘤和恶性肿瘤。其中，乳腺良性肿瘤中较多见的有腺纤 维瘤和管内或囊内乳头状瘤，而乳腺恶性肿瘤有癌、肉瘤及癌肉瘤等。为了对症下药，需要 对乳腺肿瘤进行定性的分析和诊断。
[0003] 目前乳腺钥靶、B超以及核磁共振是检测乳腺肿瘤最常用的三种成像方式。其中， 核磁检查是软组织分辨率最高的影像检查手段，与钥靶和B超相比核磁共振图像可以更好 地显示肿瘤形态和血液动力学的特征，大大提高了乳腺癌的诊断率，发现乳腺肿瘤的敏感 性高达94%?100%，特异性为71 %，且能发现X线片上不能发现的小肿瘤，尤其对致密型 乳腺及乳腺癌术前分期都有优势。乳腺核磁共振成像在检测重度危险病人和阶段病人都起 到非常重要的作用。
[0004] 为了给临床和放射科医生提供有效的乳腺肿瘤的形状特征以及体积等信息以帮 助医生制定治疗计划，需要将乳腺肿瘤区域从核磁共振图像中分割出来。目前，对乳腺核磁 共振图像进行分析主要依赖专家手工分割与借助计算机辅助分割。手工分割可以较精确地 得到图像中乳腺肿瘤区域，但手动分割需要较强的先验知识，并且需要消耗大量时间。为了 解决该问题，很多专家致力于构造乳腺肿瘤的自动化分割方法，希望通过计算机辅助诊断 的方式将医生从海量的医学数据的分割中解放出来。但是，受成像设备的限制以及呼吸导 致的运动伪影，临床采集的乳腺核磁共振图像的灰度分布非常不均匀，此外乳腺肿瘤区域 的灰度与乳房中的血管以及肋骨的灰度非常接近，因此传统的基于图像灰度的分割方法无 法正确检测出乳腺肿瘤区域。因此，设计出一种快速且准确的乳腺肿瘤的自动化分割方法 显得非常重要。
[0005] 医学图像的自动化分割是进行计算机图像分析的关键，是制约医学图像处理中其 他相关技术发展与应用的瓶颈，也是医学图像理解的基础。医学图像分割在生物医学研究、 临床诊断、病理分析等方面具有重要的意义。近几十年来，许多研究人员在医学图像分割 领域做了大量的研究工作，提出了很多有效的分割算法。其中，统计学理论、模糊集理论、 神经网络、形态学方法、小波理论等已经在医学图像分割中受到了广泛的关注和应用。此 夕卜，近期涌现出的遗传算法、尺度空间、多分辨率方法、图论以及非线性扩散方程等新方法 也不断的背应用于医学图像分割的领域中。
[0006] 根据医学图像分析任务的要求，可将分割分为基于区域的分割和基于边界的分 割两类。基于区域的分割是对一幅图像按其特征相似性划分为若干有意义的区域。该类 方法可以进一步的分为基于阈值的分割方法，基于区域生长和分裂合并的方法，基于聚 类的方法，基于随机场的方法等。基于边界的分割是将图像中感兴趣的区域（Region of Interesting，R0I)分离出来，提取其边界，为R0I进行准确定位。该类方法主要有基于曲面 拟合的方法，基于边界曲线拟合的方法，基于反应-扩散方程的方法，串行边界查找，基于 形变模型的方法等。
[0007] 由于辅助临床医疗的需要，如何提高图像分割模型的速度、精度是个重点、难点问 题。乳腺肿瘤的分割主要存在如下三个难点问题：1)乳房核磁共振图像受偏移场的影响非 常严重，图像灰度分布非常不均匀；2)图像中除了含有乳房部分还包括肋骨、胸腔等非乳 房部分，且这两个部分并没有非常明显的边界；3)在乳房区域内乳腺肿瘤的灰度与血管的 灰度非常接近，很难从灰度上区分这两个区域。由于这几个难点问题，上述提到的医学图像 分割方法都无法直接应用于乳腺肿瘤的分割中。因而设计一种完全自动化的分割方法精确 的将乳腺肿瘤从核磁共振图像中提取出来具有较大的挑战性。
[0008] 对于核磁共振图像分割而言偏移场现象的影响最为明显。由于偏移场往往比较光 滑，许多学者在此基础上提出基于多项式的模型，基于离散正余弦变换的模型以及基于薄 板样条的模型等。这些方法的精确度由自由度确定，随着自由度的增加，导致计算复杂度的 增加，且易陷入局部最优。还有一些模型假设在图像由有限个同质区域构成，利用这个性质 将聚类模型引入到偏移场的计算中。Wells等人提出一种基于最大似然估计的方法，然而该 方法需要知道每个类别的具体信息；目前还出现了基于聚类的方法，该类方法往往依赖初 始值的选取且易陷入局部极值，当偏移场较强时该类方法往往不能得到较好的结果。
[0009] 检测乳腺肿瘤时，一些非乳房区域的部位例如胸腔，肺以及心脏等区域会影响肿 瘤检测的速度和精度，因此需要在乳腺肿瘤检测之前将乳房部位作为感兴趣区域分割出 来，去除其他非感兴趣区域。由于心脏和肺部的运动，乳房区域与其他区域的边界非常模 糊，使得感兴趣区域的提取存在困难。Hayton等人基于数学形态学方法和图搜索方法提出 一种乳房区域的分割方法，但该方法非常耗时，且无法处理病人的胸腔不是非常平坦时的 情况。Twellmann等人将中值滤波、灰度阈值法以及形态学方法相结合用以分割乳房区域。 但由于核磁共振图像的灰度不均匀，该方法的灰度阈值部分无法正确的将乳房区域与非感 兴趣区域分离开。Ertas等人利用细胞神经网络方法代替形态学操作可以有效的将乳房区 域分割出来。但该方法是一个迭代算法，需要根据不同的情况设置迭代中止条件。
[0010] 为了测量肿瘤的形状、体积以及亮度等特征，需将乳腺肿瘤从乳房部位分割出来。 全局阈值法是一个常用的图像分割方法，该方法基于全局灰度信息。该方法可以将灰度值 高于正常部位的灰度值的乳腺肿瘤区域分割出来。但某些乳腺肿瘤的灰度值经常与其他部 位叠加，因此该方法不适用于乳腺肿瘤的分割。Li等人提出了局部自适应阈值方法，与全局 阈值法相比局部阈值法的分割方法更加准确，但分割出的目标边界都不光滑影响后续的肿 瘤性质分析。Szekely等人用马尔科夫随机场将粗分割结果进一步提炼，再结合半径梯度索 引法，直方图方法以及二值化方法将乳腺肿瘤从乳房区域中分割出来，该方法得到的目标 边界比较准确但实现过程比较复杂。Zheng等人提出了自适应的拓扑区域增长方法定义肿 瘤区域的初始边界，并应用主动轮廓模型方法捕获目标边界。Petrick等人使用高斯滤波的 拉普拉斯变化以及密度权重对比增强法分割乳腺肿瘤。密度权重对比增强方法可以增强图 像中的结构信息使得基于边缘检测的方法更好的捕获目标边界。Zou等人使用梯度向量流 (GVF)方法分割乳腺肿瘤，但初始曲线的自适应选择仍有待研究。
[0011] 综上所述，由于乳房核磁共振图像灰度不均匀，非乳房区域对乳腺肿瘤的自动分 割影响较大以及乳腺肿瘤的灰度与血管灰度非常接近等问题导致自动的从核磁共振图像 中分割出乳腺肿瘤显得非常困难。


【发明内容】

[0012] 为解决上述问题，本发明公开了一种乳腺肿瘤分割方法，将偏移场信息融入到分 类模型中将二者归纳为统一框架，通过快速能量极小化方法同时求解乳腺核磁共振图像分 类以及矫正偏移场，使得两者在模型演化过程中相互利用对方的信息最终实现二者的准确 求解；并对血管以及肿瘤的形状进行分析，抓住二者在形状上的差异，利用特征值以及特征 向量等参数建立形状先验，构造基于形状的水平集驱动力并将其与基于局部信息的水平集 方法相结合，使得水平集在演化时克服管状结构的干扰只捕获乳腺肿瘤区域。
[0013] 为了达到上述目的，本发明提供如下技术方案：
[0014] 一种基于核磁共振图像的乳腺肿瘤分割方法，包括如下步骤：
[0015] 步骤1，构造乳腺组织核磁共振图像的分类与偏移场矫正的耦合框架：
[0016] 步骤1. 1，对乳房核磁共振的灰度分布进行分析，确定乳腺核磁共振图像的类别个 数；
[0017] 步骤1. 2,增加保证偏移场光滑的约束项，结合多相位的水平集方法构造乳房核磁 共振图像的分类与偏移场矫正的耦合模型：
[0018]

【权利要求】
1. 基于核磁共振图像的乳腺肿瘤分割方法，其特征在于，包括如下步骤： 步骤1，构造乳腺组织核磁共振图像的分类与偏移场矫正的耦合框架： 步骤1.1，对乳房核磁共振的灰度分布进行分析，确定乳腺核磁共振图像的类别个数； 步骤1. 2,增加保证偏移场光滑的约束项，结合多相位的水平集方法构造乳房核磁共振 图像的分类与偏移场矫正的耦合模型：
其中第一项是数据拟合项，该项包括灰度信息I(X)、偏移场信息b(X)以及类别信息Ui(X),而第二项是空间约束项； 步骤1. 3,应用快速的能量最小化方法求解图像分类和偏移场； 步骤2,利用图像边缘增强的方法对乳房区域以及周边区域进行增强，区分乳房区域和 非乳房区域，提取乳房区域； 步骤3,结合形状先验对乳房肿瘤图像进行分割： 步骤3.1，构造初始曲线； 步骤3. 2,根据形状统计的方法构造管状以及块状结构的先验信息； 步骤3. 3,将形状先验信息融入基于局部信息的水平集方法进行图像的分割。
2. 根据权利要求1所述的基于核磁共振图像的乳腺肿瘤分割方法，其特征在于，所述 步骤 I. 1包括如下步骤： Step1:统计图像的直方图并进行高斯平滑，得到平滑后的直方图H; Step2:对平滑后的直方图进行初始波峰检测，若Hi >Hh且Hi >Hi+1则为波峰，否则 为非波峰；
Step3:去除小波峰和相邻波峰，设定阈值T1，若波峰高度小于阈值时，从集合中去除 该波峰；
当两个波峰之间的宽度小于阈值T2时，将两个波峰合并，去除较小的波峰，最终得到波 峰的集合P; Step4:根据波峰集合，确定波峰的数量，以此作为图像的类别个数，波峰对应的灰度 值则为各个类别的初始中心； Step5 :用融入偏移场信息与类别的均值信息（b(X)Ci)表示每个类别。
3. 根据权利要求1或2所述的基于核磁共振图像的乳腺肿瘤分割方法，其特征在于： 所述步骤1. 2的耦合模型中第一项中的偏移场信息b建模为一组光滑的基函数的线性组 合： b=wTg 其中，W=O1,O2,…，〇n)T是参数，g= (gl,g2,…，gn)T是基函数。
4. 根据权利要求1所述的基于核磁共振图像的乳腺肿瘤分割方法，其特征在于：所述 步骤 3. 2包括如下步骤： Step1:用高斯滤波器求解每个点的Hessian矩阵； Step2:求解Hessian矩阵的特征值入丨，入2; Step3:根据局部结构特征与Hessian矩阵特征值的关系定义如下管状结构先验：
其中，Rb= 0，(：是自定义参数，
'根据多尺度的概念，最终的先验 公式表示如下：
当Vtj(X)接近1时，当前点X的结构是管状结构，反之是块状结构。
5. 根据权利要求1所述的基于核磁共振图像的乳腺肿瘤分割方法，其特征在于，所述 步骤3. 3中在模型中增加形状先验，具体能量泛函如下：
其中，M1O(y)) =H〇 (y))，M2O(y)) =l-H〇 (y))，H是Heaviside函数， 对图像中所有点x，其局部高斯概率的形式如下：
而Pr1 (I(x))和pr2 (I(x))分别表示非管状结构和管状结构的先验： Pr1(Kx)) =P(XG^non-vessel) =I-V0(x) ,pr2(I(x) ) =P(xG^vessei) =V0(x). 能量泛函公式的极小化可通过梯度下降法求解：
其中，
》通过交替迭 代上述三个公式即可得到最终的肿瘤分割结果。
【文档编号】G06T7/00GK104268873SQ201410500096
【公开日】2015年1月7日 申请日期:2014年9月25日 优先权日:2014年9月25日
【发明者】顾升华, 詹天明, 郑钰辉, 陈允杰, 罗君 申请人:南京信息工程大学