一种基于随机游走和图割的活动轮廓模型的MR图像三维交互分割方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN106780518A(43)申请公布日2017.05.31(21)申请号201710073688.X(22)申请日2017.02.10(71)申请人苏州大学地址215123江苏省苏州市工业园区仁爱路199号(72)发明人陈新建孙敏(74)专利代理机构南京纵横知识产权代理有限公司32224代理人董建林(51)Int.Cl.G06T7/11(2017.01)G06T7/136(2017.01)G06T7/149(2017.01)权利要求书3页说明书10页附图1页(54)发明名称一种基于随机游走和图割的活动轮廓模型的MR图像三维交互分割方法(57)摘要本发明公开一种基于随机游走和图割的活动轮廓模型的MR图像三维交互分割方法，其包括步骤：S1，选取种子点截取包含垂体瘤的局部三维MR图像数据，进而获取分割的初始边界曲面；S2，基于初始边界曲面，建立混合活动轮廓模型，得到模型能量函数；S3，对能量函数进行离散化；S4，将截取的局部三维MR图像的每个体素点作为图的节点，采用每个体素点的6邻域进行图的构建；对初始边界曲面内外的体素点分别赋予初始值；并根据离散后的能量函数给节点间、节点与源点间，和节点与汇点间连接的边分别赋予对应的权值；S5，基于构建的图进行图割计算，得到分割结果，提取分割结果中的边界曲面即为分割轮廓；S6，以当前分割轮廓替代S1中的初始边界曲面，重复步骤S2至S5，直至分割结果收敛，得到最终分割轮廓。利用本发明可实现垂体瘤MR影像的CN106780518A三维分割，提高垂体瘤图像分割的准确性。CN106780518A权利要求书1/3页1.一种基于随机游走和图割的活动轮廓模型的MR图像三维交互分割方法，其特征是，包括以下步骤：S1，获取初始边界曲面，包括步骤：S11，选取随机游走算法所需的种子点，以该种子点作为三维中心点，从脑部MR三维数据中截取包含垂体瘤的局部三维MR图像数据；S12，利用三维随机游走算法对截取的局部三维MR图像数据进行处理，得到各体素点到达种子点的概率三维图；S13，利用最大期望值算法选取概率阈值，以分割出疑似目标区域，将疑似目标区域内的体素作为活动轮廓模型算法的前景点，疑似目标区域的边界曲面即初始边界曲面；S2，基于初始边界曲面，建立混合活动轮廓模型：对标准活动轮廓模型的能量函数中边界能量项采用标准测地线模型，区域能量项采用局部高斯模型描述的灰度最大后验概率，进而得到混合活动轮廓模型的能量函数；S3，模型离散化，包括：S31，给每个体素点定义一个二元变量，赋值1和0分别对应代表前景点和背景点，以将区域能量项函数离散化；S32，利用割测度将边界能量项函数即测地线模型离散化；S33，得到离散化后的能量函数；S4，构建图：对于截取的局部三维MR图像数据，将其中的每个体素点作为图的节点，采用每个体素点的6邻域进行图的构建；对初始边界曲面内和初始边界曲面外的体素点分别赋予初始值；并根据离散后的能量函数给节点间连接的边，节点与源点连接的边，和节点与汇点连接的边分别赋予对应的权值；S5，基于步骤S4构建的图进行图割计算：采用最大流和最小割算法得到分割结果，提取分割结果中的边界曲面即为分割轮廓；S6，以当前分割轮廓作为初始边界曲面，重复迭代步骤S2至步骤S5，直至分割结果收敛，输出最终分割结果。2.根据权利要求1所述的方法，其特征是，还包括步骤S7，利用三维中值滤波对步骤S5得到的演化曲面进行处理，得到滤波后分割结果；步骤S6以当前滤波后的分割轮廓作为初始边界曲面。3.根据权利要求1或2所述的方法，其特征是，步骤S11中，种子点的选取步骤为，在矢状位的二维显示图像中，选取中心位置切片，在切片中垂体瘤内近似中心位置选取一点，即作为种子点。4.根据权利要求1或2所述的方法，其特征是，步骤S11中，所述包含垂体瘤的数据立方体尺寸为71×71×41个体素。5.根据权利要求1或2所述的方法，其特征是，步骤S2中，给定图像域I(x):Ω→R为二维图像；闭合曲线C将图像分成两个独立区域，L(C)表示C的长度；设曲线C内区域为Ω1，曲线C外区域外为Ω2，对于图像域Ω的各像素p，考虑其半径为ρ的圆形邻域，定义为Ox＝{y:|x-y|≤ρ}，x为圆形邻域的中心像素点，y为圆形邻域中除x外的任意像素点；则区域Ox∩Ωi内的像素点灰度概率密度函数表示为：2CN106780518A权利要求书2/3页式中ui(x)和σi(x)分别为区域Ox∩Ωi内像素点x的局部灰度均值和方差，I(y)为像素点y的灰度值，边界能量项EEdge采用标准测地线模型，区域能量项EReigon采用局部高斯模型描述的灰度最大后验概率形式的能量函数E表示为：式中β是一个任意正常数，函数w是一个权重函数。