基于超体素的多模态MRI脑肿瘤分割方法技术

本发明专利技术公开一种基于超体素的多模态MRI脑肿瘤分割方法，首先对传统的SLIC算法进行改进，通过计算强度距离、位置距离和梯度距离确定当前体素点到其领域内各体素块种子点的距离，由此获得当前体素块的边缘信息，从而可准确获得有明显占位效应的脑肿瘤或受侧脑室挤压导致的明显肿瘤周围水肿和瘤内出血的弱边界信息，克服噪声体素对超体素边界的影响，获得整齐紧凑的边界分割信息。其次，在进行单模态分割结果信息整合时，通过横向比较的方式进行信息整合，从而获得具有相同分割结果的四张MRI图像超体素分割边界信息。MRI图像超体素分割边界信息。MRI图像超体素分割边界信息。

【技术实现步骤摘要】
基于超体素的多模态MRI脑肿瘤分割方法


[0001]本专利技术涉及医学图像分割领域，尤其是一种基于超体素的多模态MRI脑肿瘤分割方法。

技术介绍

[0002]正电子断层显像仪(positron emission tomography,PET)通过向人体注射放射性药物从而获得人体内部代谢信息，但其无法直观的查看人体内部结构信息，且放射性药物通常对人体有害。反映人体结构的计算机断层扫描(Computer Tomography,CT)图像利用人体各部分组织对X射线的吸收系数不同的原理成像，具有较高的时效性和空间分辨率，但其对人体损伤较大。在不追求成像速度时通常采用对患者身体健康损伤较小的磁共振图像(Magnetic Resonance Images,MRI)。单模态MRI图像有FLAIR图像、T1加权图像(T1 weighted image)、T2加权图像(T2 weighted image)、T1-Contrast(T1-Gd，在T1加权图像中使用造影剂)，MRI图像具有较高的精准度，可通过不同扫描序列获得反映不同结构信息的图像，具有较高的空间分辨率和时间分辨率，已应用于脑肿瘤检测。如T1加权像中各个脑组织的对比度较为明显，对于区分脑部健康组织和病灶组织有着重要的作用；T2加权图像可以有效分离脑部水肿区域，获得异常区域边界信息；T1-Contrast图像能清楚直观地展示肿瘤边缘，从而得到核心肿瘤区域的边缘；FLAIR图像可以体现脑肿瘤区域与大脑其他组织之间有着明显的差异，较为明显区分脑部白质、灰质和脑脊液的结构。然而，由于脑肿瘤的呈现上没有固定的大小形状并且其可能出现在大脑的任何位置，再者颅内肿瘤可能会对周围脑组织产生影响，不仅损害血脑屏障，增加血管壁的通透性，还会对周围脑组织产生压迫，形成局部水肿等，故在进行单模态MRI图像脑肿瘤分割时经常会遇到很多困难。
[0003]目前，已有用多模态MRI脑肿瘤分割方法。具体方法是对每个单模态MRI图像进行分割处理，得到各自的分割结果图，然后在将所得到的单模态分割结果图进行信息整合，最终得到多模态分割结果。存在如下问题：
[0004]首先，基于梯度超体素分割算法是现有单模态MRI图像进行分割处理的方法之一，其经典算法有SLIC算法。SLIC算法的核心思想是先将一幅图像划分为多个区域块，将每个块的几何中心作为该块的初始中心，然后在该点的邻域内找到一个梯度最小的点作为体素块种子点，通过计算当前体素点与其领域内的体素块种子点之间的距离，不断迭代直至收敛阈值。然而，现有SLIC算法计算距离时只有强度距离和位置距离获得该块的边界信息，不能较好的获得有明显占位效应的脑肿瘤或受侧脑室挤压导致的明显肿瘤周围水肿和瘤内出血的弱边界信息，无法得到光滑、准确且最大程度贴合3D脑部图像边缘的体素块。
[0005]其次，在进行单模态分割结果信息整合时，传统的方法是将获得的各模态分析的结果叠加到一个MRI图像上获得最终分割结果，忽略了四个MRI图像中所存在的内在联系。Calhoun等人已经意识传统整合方法的不足之处，提出了联合独立成分分析的方法，即通过串联的方式将单模态信息融合，使用盲源信号分离图像中一些不相关的成分，从而提高检测效率。但是串联融合导致图像中各个结构联系过于紧密，易出现过拟合的情况。

技术实现思路

[0006]本专利技术是为了解决现有技术所存在的上述技术问题，提供一种基于超体素的多模态MRI脑肿瘤分割方法。
[0007]本专利技术的技术解决方案是：一种基于超体素的多模态MRI脑肿瘤分割方法，依次按照如下步骤进行：
[0008]步骤1：对图像I进行预处理
[0009]输入序号为V的单模态MRI脑部三维图像I，找到三维图像I矩阵中的最大值Max_I和最小值Min_I，得到归一化后的三维图像I'：
[0010][0011]所述V＝{1,2
…
,n}，所述n为多模态MRI脑部三维图像中的单模态MRI的总数；
[0012]步骤2：初始体素块中心
[0013]将三维图像I'划分为k个大小相等的体素块S，每个体素块内包含大约个体素点，其中，N表示三维图像I'中的体素个数，以体素块S内部体素的坐标x,y,z平均值为初始体素块中心；
[0014]步骤3：寻找梯度最小点为体素块种子点
[0015]以初始体素块中心点为中心，在其邻域内进行多次迭代，寻找到一个梯度最小的体素点，并将该点定义为体素块种子点；
[0016]步骤4：计算体素点i与其邻域内的体素块种子点j之间的距离D
[0017][0018][0019][0020][0021]所述D
c
为强度距离、D
s
为位置距离、D
T
为梯度距离；所述V
i
和V
j
分别表示第i个体素点和当前的体素块种子点j的矩阵归一化的值，所述i＝{1,2
…
N,}，所述j＝{1,2
…
c,}，所述c为体素点i邻域内的体素块种子点的总数；(x
i
,y
i
,z
i
)表示第i个体素点的坐标，R
x
、R
y
和R
z
分别表示相应方向上的体素分辨率，(x
j
,y
j
,z
j
)表示当前的体素块种子点j的坐标；T
i
和T
j
分别表示第i个体素点和当前的体素块种子点j的梯度值；其中ω表示梯度特征D
T
的紧致系数；N
c
表示图像的最大色彩，m表示紧密度系数，N
s
表示空间距离；
[0022]步骤5.判断距离D是否第一次计算，是，则进行步骤6，否，进行步骤9；
[0023]步骤6.将体素块种子点j的标签分配给体素点i；
[0024]步骤7.将计算的距离值D保存为Dmin；
[0025]步骤8.判断j＝c，是，则进行步骤10，否，则j
←
j+1后返回步骤4；
[0026]步骤9.判断距离值D是否小于Dmin，是，返回步骤6，否，则返回步骤8；
[0027]步骤10.判断i＝N，是，进行步骤11，否，则i
←
i+1后返回步骤4；
[0028]步骤11.获得序号为V的三维图像I'的分割结果图I”；
[0029]步骤12.判断V＝n，是，则进行步骤13，否，则V
←
V+1后返回步骤1；
[0030]步骤13.多模态信息整合
[0031]步骤13.1判断序号为V的分割结果图I”中是否有两个具有相同特征的体素块S1、S2，否，进行步骤13.2，是，进行步骤13.3；
[0032]步骤13.2判断V＝n，是，进行步骤14，否，则V
←
V+1返回步骤13.1；
[0033]步骤13.3判定除序号本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于超体素的多模态MRI脑肿瘤分割方法，其特征在于依次按照如下步骤进行：步骤1：对图像I进行预处理输入序号为V的单模态MRI脑部三维图像I，找到三维图像I矩阵中的最大值Max_I和最小值Min_I，得到归一化后的三维图像I'：所述V＝{1,2
…
,n}，所述n为多模态MRI脑部三维图像中的单模态MRI的总数；步骤2：初始体素块中心将三维图像I'划分为k个大小相等的体素块S，每个体素块内包含大约个体素点，其中，N表示三维图像I'中的体素个数，以体素块S内部体素的坐标x,y,z平均值为初始体素块中心；步骤3：寻找梯度最小点为体素块种子点以初始体素块中心点为中心，在其邻域内进行多次迭代，寻找到一个梯度最小的体素点，并将该点定义为体素块种子点；步骤4：计算体素点i与其邻域内的体素块种子点j之间的距离D计算体素点i与其邻域内的体素块种子点j之间的距离D计算体素点i与其邻域内的体素块种子点j之间的距离D计算体素点i与其邻域内的体素块种子点j之间的距离D所述D
c
为强度距离、D
s
为位置距离、D
T
为梯度距离；所述V
i
和V
j
分别表示第i个体素点和当前的体素块种子点j的矩阵归一化的值，所述i＝{1,2
…
N,}，所述j＝{1,2
…
c,}，所述c为体素点i邻域内的体素块种子点的总数；(x
i
,y
i
,z
i
)表示第i个体素点的坐标，R
x
、R
y
和R
z
分别表示相应方向上的体素分辨率，(x
j
,y
j
,z
j
)表示当前的体素块种子点j的坐标；T
i
和T
j
分别表示第i个体素点和当前的体素块种子点j的梯度值；其中ω表示梯度特征D

【专利技术属性】
技术研发人员：方玲玲，王欣，
申请(专利权)人：辽宁师范大学，
类型：发明
国别省市：