基于稀疏低秩与Atlas集的3D MRI胰腺分割方法
【专利摘要】本发明提出一种基于稀疏低秩与Atlas集的3D MRI胰腺分割方法，主要解决现有技术中分割胰腺粘连现象严重、分割结果不精确的问题。其步骤如下：1)输入3D MRI数据并归一化处理；2)对归一化处理后的序列图像进行稀疏低秩分解并进行Hessian矩阵增强；3)在增强后的图像上人机交互分割出胰腺标签；4)构造粗分割集；5)输入Atlas集，依次进行刚性配准和弹性配准，并进行标签融合，输出分割结果。本发明运用矩阵增强、矩阵分解、图集法，并利用3D MRI图像的空间信息，进行胰腺分割，减少了胰腺目标与其他组织器官的粘连现象，图像轮廓清晰，提高了分割的准确度，可用于胰腺癌精确放射治疗中的胰腺目标检测。
【专利说明】
基于稀疏低秩与At I as集的3D MR I胰腺分割方法
技术领域
[0001]本发明属于医学图像处理技术领域，更进一步涉及三维核磁共振图像3D MRI胰腺 分割，可用于胰腺癌精确放射治疗中的胰腺目标检测。
【背景技术】
[0002] 胰腺癌是出现在胰脏的癌症，通常认为胰腺癌是常见肿瘤中恶性程度最高、死亡 率最高的癌症，约90%的患者无法以手术根治，五年存活率低于5%。根据世界卫生组织估 计，胰腺癌在全世界发生率排第十三，死亡率排第八。随着科技的发展，其他癌症的治疗都 有了提高，但胰腺癌却没有明显的突破，死亡率也没有明显的下降趋势，所以必须加强研 究，早日攻克胰腺癌对人类的威胁。
[0003] 传统检查胰腺癌的方式是医生根据经验和观察直接阅片，定性的判断患者的病 情，这样易受情绪影响，且随着工作量增加，阅片准确率下降。对胰腺癌的治疗方法中，放射 治疗是一种主要的方法，他试图用射线杀死肿瘤，但也会对正常组织造成一次性或永久性 伤害。精确放疗通过精确的肿瘤定位、精确的计划设计、精确的剂量计算及在治疗机上精确 执行的一种全新的肿瘤放疗技术，有助于减少放射治疗中对健康组织放射造成的伤害。精 确分割胰腺组织有助于帮助医生分析胰腺病情，也是精确放疗不可缺少的一步。
[0004] 目前，对胰腺癌进行检测和治疗的设备主要是核磁共振成像MRI系统，该设备功能 强大、技术含量高、软组织图像清晰，而且对人体检测损害小，已成为对胰腺进行分析、诊断 和治疗的重要手段，而从MRI图像中精确的分割出胰腺，就是这些手段的重要前提。
[0005] 传统医学图像分割方法包括基于区域的方法、基于边缘的方法、基于神经网络的 方法、基于模糊聚类的方法、基于图谱的方法等。由于这些方法大多数是基于2D图像序列的 分割，利用2D图像本身的特性信息完成感兴趣区域的分割，例如灰度、纹理等信息，而没有 利用体素的空间信息，因此不能进行准确的胰腺提取;医学MRI图像有噪声强、对比度低、软 组织边界不明显的特点，而这些方法大多数对图像噪声和对比度敏感，所以难以定位目标 边界，使得分割结果不精确。
【发明内容】
[0006] 本发明的目的在于针对以上现有技术的不足，提出一种基于稀疏低秩与Atlas集 的3D MRI胰腺分割方法，以提高3D MRI胰腺分割的精度。
[0007] 实现本发明目的的技术方案如下：
[0008] (1)输入三维核磁共振3D MRI序列图像，对其进行归一化处理，得到归一化后的序 列图像C;
[0009] (2)用交替迭代法对归一化处理后的序列图像C进行稀疏低秩矩阵分解，得到稀疏 图像序列A和低秩图像序列B;
[0010] (3)对稀疏图像序列A进行Hessian矩阵增强，得到二值增强图像序列A1;
[0011] (4)通过一次人机交互，在二值增强图像序列六:中确定胰腺的位置区域P，根据胰 腺位置区域P，利用形态学原理在二值增强图像序列^中分割出胰腺的标签图像序列L;
[0012] (5)根据标签图像序列L和步骤（1)中的序列图像C，得到粗分割集Q，其中包含F帧 图像；
[0013] (6)输入胰腺图像序列W及人工勾画出的胰腺标签序列E，即Atlas集；
[0014] (7)选择粗分割集Q中的第i帧图像作为参考图像，记为&，将Atlas集中的胰腺图 像序列W作为浮动图像，用基于互信息的刚性配准算法对QnW进行配准，得到刚性配准后的 胰腺图像序列I和刚性转换公式!^，根据^对胰腺标签序列E进行迀移，得到刚性配准后的 胰腺标签序列Ei，其中i的取值范围为1~F;
[0015] (8)将刚性配准后的胰腺图像序列1作为浮动图像，使用基于Demons的弹性配准 算法对参考图像&和该1进行配准，得到弹性配准后的胰腺图像序列％和弹性转换公式T 2， 根据^对刚性配准后的胰腺标签序列E1进行迀移，得到弹性配准后的胰腺标签序列Ε 2;
[0016] (9)计算参考图像&与弹性配准后的胰腺图像序列此的相关性序列，根据该相关性 序列及弹性配准后的胰腺图像序列此和弹性配准后的胰腺标签序列扮这三者是一一对应的 关系，取相关性序列最高的前五个系数对应的标签，并对这五个标签进行融合，得到参考图 像Qi精分割后的结果Ri，其中i的取值范围为1~F;
[0017] (10)重复步骤(7)~(9)，对粗分割集Q中所有图像精分割后，得到最终结果精分割 集R 〇
[0018] 本发明与现有技术相比具有以下优点：
[0019] 1.本发明采用Hessian矩阵增强的方法，将三维核磁共振图像3D MRI的空间信息 运用到胰腺目标的分割中，使分割结果更加准确。
[0020] 2.本发明采用稀疏低秩矩阵分解方法，将三维核磁共振图像3D MRI分解成低秩序 列部分和稀疏序列部分，通过对稀疏序列部分进行操作，有效减少目标与其他器官的粘连 现象。
[0021] 3.本发明运用Atlas集方法，将粗分割结果与胰腺模板库图像配准，再进行标签迀 移和融合得到最终分割结果，使得分割目标区域能保持较好的轮廓。
【附图说明】
[0022]图1是本发明的实现流程图；
[0023]图2是待分割3D MRI图像序列中的一帧及该帧的标签；
[0024]图3是用本发明与现有level set法和graph cut法的分割结果对比图。
【具体实施方式】
[0025]参照图1，本发明的实现步骤如下：
[0026]步骤1、输入数据并归一化处理。
[0027]输入三维核磁共振图像3D MRI序列Co,按下列公式对其进行归一化处理，得到归 一化处理后的序列图像C:
[0029]其中，max〇为提取矩阵最大元素。
[0030] 步骤2、对归一化处理后的序列图像进行稀疏低秩分解。
[0031] 用交替迭代法对归一化处理后的序列图像C进行稀疏低秩矩阵分解，得到稀疏图 像序列A和低秩图像序列B:
[0032 ] (2a)设置初始值，稀疏图像序列A为零矩阵，低秩图像序列B = C;
[0033] (2b)通过公式：
更新稀疏矩阵A，其中， 是欧几里得投影，γ为常数，且γ>〇，β指偏离线性约束的惩罚参数，且i3>〇，k为迭代 次数，Bk为第k次迭代后的低秩矩阵，Zk为第k次迭代后的约束矩阵；
[0034] (2c)对
进行奇异值分解产生三个不同的矩阵：左奇异向量矩 阵Uk+1、右奇异向量矩阵Vk+1及对角奇异值矩阵
[0035] (2d)通过公式：
更新低秩矩阵Β，其中，Τ 为矩阵的转置，max{}为提取最大元素，diag()为提取对角线元素；
[0036] (2e)通过公式Zk+1 = Zk-f3(Ak+1+Bk+1-C)更新约束性矩阵Z，其中，A k+1是第k+Ι次迭代 后的稀疏矩阵，zk+1是第k+Ι次迭代后的线性约束矩阵，Bk+1是第k+Ι次迭代后的低秩矩阵； [0037] (2f)设定迭代次数k=1000,重复步骤(2a)~（2e)，迭代结束后，得到稀疏图像序 列A和低秩图像序列B。
[0038]步骤3、对稀疏图像序列进行Hessian矩阵增强。
[0039]对稀疏图像序列A进行Hessian矩阵增强，得到二值增强图像序列A1:
[0040] (3a)根据下面公式，求出稀疏图像序列A中每个像素点I的Hessian矩阵Η: 4- 4, Α,
[0041 ] Η = Iyx Ιγζ 4ν '二
[0042]其中，I表示稀疏图像序列Α中每一个像素点，Ιχχ代表沿χ方向的二阶偏微分，1#代 表先沿X方向再沿y方向的二阶偏微分，Ιχζ代表先沿X方向再沿ζ方向的二阶偏微分，Iyx代表 先沿y方向再沿X方向的二阶偏微分，I yy代表沿y方向的二阶偏微分，Iyz代表先沿y方向再沿 z方向的二阶偏微分，IZXR表先沿z方向再沿X方向的二阶偏微分，I zyR表先沿z方向再沿y 方向的二阶偏微分，Izz代表沿Z方向的二阶偏微分；
[0043] (3b)根据Choiesky分解法，计算步骤(3a)中每一个Hessian矩阵Η的三个特征值， 并按绝对值从小到大顺序排序为λι，λ2，λ3;
[0044] (3c)根据下面公式和每一个像素点I的Hessian矩阵特征值，对稀疏图像序列Α的 每一个像素点进行增强；
if Xt < 0 or 1? < 〇
[0046]其中，α是起调节面状权重RA作用的常数，α >〇，τ是起调节管状权重RB作用的常数， t>0，c是起调节球状权重S作用的常数，
νι(λ)是指每个像素点I增强后的结果。
[0047]步骤4、人机交互分割出胰腺标签。
[0048] 利用matlab上的roipoly函数，在步骤3中得到的二值增强图像序列Αι上手动选择 胰腺的位置区域P，利用形态学原理，将Ai中与P组成连通域的部分提取出来，得到胰腺的标 签图像序列L:
[0049] (4a)对二值增强图像序列Μ进行形态学腐蚀操作，得到腐蚀后的二值增强图像序 列Αι' ：
[0050] (4b)根据胰腺位置区域P从腐蚀后的二值增强图像序列心中提取出连通区域I/;
[0051] (4c)对提取出来的连通区域1/进行形态学膨胀操作，得到标签图像序列L。
[0052]步骤5、构造粗分割集。
[0053] 根据标签图像序列L和步骤（1)中的序列图像C，按下列公式，得到包含F帧图像的 粗分割集Q:
[0054] Q = C.*L
[0055] 其中，.*表示矩阵对应元素相乘。
[0056] 步骤6、输入Atlas集。
[0057]输入胰腺图像序列W及人工勾画出的胰腺标签序列E，W和E构成Atlas集。
[0058]步骤7、基于互信息的刚性配准。
[0059] (7a)选择粗分割集Q中的第i帧图像作为参考图像，记为&，将Atlas集中的胰腺图 像序列W作为浮动图像，用基于互信息的刚性配准算法对QnW进行配准，得到刚性配准后的 胰腺图像序列1和刚性转换公式?\:
[0060] (7a 1)设初始状态k = 0，初始参数(α?ρ3°)为(0，0，0)，W, = W，其中，？1代表水平 变量、Ρ2代表垂直变量、Ρ3代表顺时针旋转变量，k为迭代次数，为迭代第0次后的刚性配 准胰腺图像序列；
[0061] (7a2)用PV插值法统计QdPWf的联合直方图，根据该联合直方图计算QdPWf的互 信息值，其中，Wf为迭代第k次后的刚性配准胰腺图像序列；
[0062] (7a3)用Powell算法根据步骤(7a2)中得到的互信息值判断参数(pf，g，g)是否最 优，如果不是最优，则进行步骤（7a4)，如果是最优，则停止迭代，输出最优参数 (A，巧，凡）=(A"，，乂），进行步骤( 7a5)，其中，Ρ?为迭代第k次后的水平变量，K为迭代第 k次后的垂直变量，g为迭代第k次后的顺时针旋转变量；
[0063] (7a4)对妒进行平移、旋转，得到Wt1和调整后的参数(pf+1,M +1,g+1)^k = k+l， 返回到步骤(7a2)，其中，WP1为迭代第k+1次后的刚性配准胰腺图像序列，< +1为迭代第k+1 次后的水平变量，g+1为迭代第k+1次后的垂直变量，为迭代第k+1次后的顺时针旋转变 量；
[0064] (7a5)根据最优参数(Pl，p2，p3)按下列公式得到刚性转换公式!^，根据叫挪转换成 刚性配准后的胰腺图像序列W 1: 1 Ο Ο cos ρ, sin ρ, 0
[0065] JJ = 0 1 0 -sin cos/? 3 〇 Ρι _ 〇 〇 1
[0066] (7b)根据1\对胰腺标签序列Ε进行迀移，得到刚性配准后的胰腺标签序列Eu
[0067] 步骤8、基于Demons的非刚性配准。
[0068] (8a)将刚性配准后的胰腺图像序列1作为浮动图像，使用基于Demons的弹性配准 算法对参考图像&和该1进行配准，得到弹性配准后的胰腺图像序列此和弹性转换公式T2:
[0069] (8al)设初始状态k=l，在参考图像Qi上选择所有像素点作为Demons点，，〗=啊， 给定允许误差ε，其中e>〇，k为迭代次数，其中时为第一次迭代的弹性配准胰腺图像序列；
[0070] (8a2)根据Symmetric Demons公式计算中所有像素点的形变向量：
[0072] 其中，u表示形变向量，uk表示迭代第k次后的形变向量表示迭代第k次后弹性 配准胰腺图像序列，▽表示求图像的梯度，I I表示取模操作，η表示归一化因子；
[0073] (8&3)根据形变向量1!1{按下列公式求出弹性转换公式2^，根据0将1； 1转换为 %?+1，其中，g为迭代第k次后的弹性转换公式，％i+1为迭代第k+Ι次后的弹性配准胰腺图像 序列：
[0074] Jf = ξ·ι/
[0075] 其中，ξ表示项力衰减系数；
[0076] (8a4)计算Qi与％?Η?的误差平方和，如果该误差平方和小于ε，输出弹性转换公式 Γ 2 =2t和弹性配准后的胰腺图像序列％ =甿"，结束迭代，否则k = k+l，返回步骤(8a2);
[0077] (8b)根据^对刚性配准后的胰腺标签序列E1进行迀移，得到弹性配准后的胰腺标 签序列E2。
[0078] 步骤9、进行标签融合。
[0079]计算参考图像&与弹性配准后的胰腺图像序列％的相关性序列，根据该相关性序 列及弹性配准后的胰腺图像序列此和弹性配准后的胰腺标签序列扮这三者是一一对应的关 系，取相关性序列最高的前五个系数对应的标签，并对这五个标签进行融合，得到参考图像 Qi精分割后的结果Ri。
[0080] (9a)根据相关性序列最高的前五个系数gl~g5和他们对应的标签图像6 1~65，计 算标签融合后的结果R/ :
[0082] 其中j的取值范围为1~5, Σ表示求和操作；
[0083] (9b)按下列公式对(9a)中标签融合后的结果R/进行修正，得到精分割后的结果 Ri： 0 I? ' <0.4
[0084] ' ?;、0,4
[0085] 步骤10、重复步骤(7)~(9)，对粗分割集Q中所有图像精分割，得到最终结果精分 割集R。
[0086] 以下结合仿真图对本发明的效果做进一步的说明。
[0087] 1.仿真条件：
[0088] 硬件平台为:CPU为Inter Core i7-4790K，主频为4GHz，内存为32GB;
[0089] 软件平台为:Windows 10专业版64位操作系统，MatlabR2014a。
[0090] 2.仿真数据：
[0091]本例中待分割的数据采用的是通过三维核磁共振仪获取的一名志愿者腹部三维 核磁共振图像序列，图像大小为320*284,一共有5帧图像，通过人工勾画出胰腺作为标签， 如图2所示，其中，图2a是胰腺图像，图2b是胰腺图像的标签。Atlas集采用的是其他志愿者 通过三维核磁共振仪获取的三维核磁共振图像，它们对应的标签图像通过人工勾画获得。 [0092] 3.仿真内容：
[0093]为了验证算法的有效性，一般在进行图像分割以后，会以标签为参考图像计算指 标，豪斯多夫距离、平均绝对表面距离、相似性、灵敏度、特异性五个指标作为分割精度的判 断标准，豪斯多夫距离、平均绝对表面距离这两个指标数值越大分割结果越精确，相似性、 灵敏度、特异性这三个指标数值越小分割结果越精确。
[0094]用本发明和现有的两个具有代表性的图像分割方法level set法和graph cut法 进行仿真实验对比。
[0095]仿真一：分别通过这三种方法对上述仿真数据进行胰腺分割仿真，对三种方法得 到的分割结果进行三维重建，结果如图3所示，其中图3a为本发明的分割结果，图3b为level set法的分割结果，图3c为graph cut法的分割结果，可以看出用本发明分割出的胰腺轮廓 更好。
[0096]仿真二：分别通过这三种方法对上述仿真数据进行胰腺分割仿真，计算分割结果 的五个指标:豪斯多夫距离、平均绝对表面距离、相似性、灵敏度、特异性，结果如表1所示： [0097] 表1用本发明与现有level set法和graph cut法的分割指标对比表
[0099]从表1可以看出本发明的分割结果在指标上优于其他两种算法，说明本发明分割 胰腺结果精度更高。
【主权项】
1. 一种基于稀疏低秩与Atlas集的3D MRI膜腺分割方法，包括： (1) 输入Ξ维核磁共振3D MRI序列图像，对其进行归一化处理，得到归一化后的序列图 像C; (2) 用交替迭代法对归一化处理后的序列图像C进行稀疏低秩矩阵分解，得到稀疏图像 序列A和低秩图像序列B; (3) 对稀疏图像序列A进行Hessian矩阵增强，得到二值增强图像序列Ai; (4) 通过一次人机交互，在二值增强图像序列Ai中确定膜腺的位置区域P，根据膜腺位置 区域P，利用形态学原理在二值增强图像序列Ai中分割出膜腺的标签图像序列レ (5) 根据标签图像序列L和步骤(1)中的序列图像C，得到粗分割集Q，其中包含F帖图像； (6) 输入膜腺图像序列W及人工勾画出的膜腺标签序列E，即Atlas集； (7) 选择粗分割集Q中的第i帖图像作为参考图像，记为化，将Atlas集中的膜腺图像序列 W作为浮动图像，用基于互信息的刚性配准算法对化、W进行配准，得到刚性配准后的膜腺图 像序列Wi和刚性转换公式Ti，根据Τι对膜腺标签序列E进行迁移，得到刚性配准后的膜腺标 签序列Ei，其中i的取值范围为1~F; (8) 将刚性配准后的膜腺图像序列Wi作为浮动图像，使用基于Demons的弹性配准算法对 参考图像化和该Wi进行配准，得到弹性配准后的膜腺图像序列化和弹性转换公式T2,根据T2 对刚性配准后的膜腺标签序列El进行迁移，得到弹性配准后的膜腺标签序列E2; (9) 计算参考图像化与弹性配准后的膜腺图像序列W2的相关性序列，根据该相关性序列 及弹性配准后的膜腺图像序列化和弹性配准后的膜腺标签序列E2运Ξ者是一一对应的关 系，取相关性序列最高的前五个系数对应的标签，并对运五个标签进行融合，得到参考图像 化精分割后的结果Ri，其中i的取值范围为1~F; (10) 重复步骤(7)~(9)，对粗分割集帅所有图像精分割后，得到最终结果精分割集R。2. 根据权利要求1所述的基于稀疏低秩与Atlas集的3D MRI膜腺分割方法，其中步骤 (2)用交替迭代法对归一化处理后的序列图像C进行稀疏低秩矩阵分解，得到稀疏图像序列 A和低秩图像序列B，按如下步骤进行： (2a)设置初始值，稀疏图像序列A为零矩阵，低秩图像序列B = C; (2b)通过公式：夏新稀疏矩阵A，其中， 是欧几里得投影，丫为常数，且丫〉〇，β指偏离线性约束的惩罚参数，且於〇，k为迭代次数，Bk 为第k次迭代后的低秩矩阵，Zk为第k次迭代后的约束矩阵； (2c)对.t行奇异值分解产生Ξ个不同的矩阵:左奇异向量矩阵护+1、 右奇异向量矩阵yk+i及对角奇异值矩阵; (2d)通过公式更新低秩矩阵B，其中，T为矩 阵的转置，max{}为提取最大元素，diag()为提取对角线元素； (2e)通过公式Zk" = zk-i3(Ak"+Bk"-C)更新约束性矩阵Z，其中，Ak"是第k+1次迭代后的 稀疏矩阵，Zk+i是第k+1次迭代后的线性约束矩阵，Bk+i是第k+1次迭代后的低秩矩阵； (2f)设定迭代次数k的值，重复步骤(2a)~(2e)，迭代结束后，得到稀疏图像序列A和低 秩图像序列B。3. 根据权利要求1所述的基于稀疏低秩与Atlas集的3D MRI膜腺分割方法，其中步骤 (3)对稀疏图像序列A进行化ssian矩阵增强，得到二值增强图像序列Ai，按如下步骤进行： (3a)根据下面公式，求出稀疏图像序列A中每个像素点I的Hessian矩阵H:其中，I表示稀疏图像序列A中每一个像素点，Ixx代表沿X方向的二阶偏微分，Ixy代表先 沿X方向再沿y方向的二阶偏微分，Ιχζ代表先沿X方向再沿Z方向的二阶偏微分，lyx代表先沿 y方向再沿X方向的二阶偏微分，lyy代表沿y方向的二阶偏微分，lyz代表先沿y方向再沿Z方 向的二阶偏微分，Izx代表先沿Z方向再沿X方向的二阶偏微分，Izy代表先沿Z方向再沿y方向 的二阶偏微分，Izz代表沿Z方向的二阶偏微分； (313)根据化〇16347分解法，计算步骤(3曰）中每一个化331曰]1矩阵邮勺^个特征值，并按 绝对值从小到大顺序排序为λ?，λ2，λ3 ; (3c)根据下面公式和每一个像素点I的化ssian矩阵特征值，对稀疏图像序列A的每一 个像素点进行增强；其中，α是起调节面状权重Ra作用的常数，α〉〇，τ是起调节管状权重化作用的常数，T〉〇，c 是起调节球状权重S作用的常数，c〉0，νι(λ)是 指每个像素点I增强后的结果。4. 根据权利要求1所述的基于稀疏低秩与Atlas集的3D MRI膜腺分割方法，其中步骤 (5)根据标签图像序列L和步骤（1)中的序列图像C，得到粗分割集Q，是将序列图像C与标签 图像序列L相乘，即： Q = C.礼 其中，.*表示矩阵对应元素相乘。5. 根据权利要求1所述的基于稀疏低秩与Atlas集的3D MRI膜腺分割方法，其中步骤 (7)中用基于互信息的刚性配准算法对Qi、W进行配准，得到刚性配准后的膜腺图像序列Wi和 刚性转换公式Τι，按如下步骤进行： (7曰)设初始状态1^ = 0，初始参数(片，此记：)为(0,0,0)，胖1*^ =胖，其中，91代表水平变量、 口2代表垂直变量、93代表顺时针旋转变量，4为迭代次数，￥1<^为迭代第0次后的刚性配准膜腺 图像序列； (7b)用PV插值法统计Qi和Wik的联合直方图，根据该联合直方图计算Qi和Wik的互信息 值，其中，Wik为迭代第k次后的刚性配准膜腺图像序列； (7c)用Powell算法根据步骤(7b)中得到的互信息值判断参数(如，知，是否最优，如果 不是最优，则进行步骤(7d)，如果是最优，则停止迭代，输出最优参数化.知/)，)^护如片）， 进行步骤(7e)，其中，片为迭代第k次后的水平变量，如为迭代第k次后的垂直变量，贫为迭 代第k次后的顺时针旋转变量； (7d)对Wik进行平移、旋转，得至ijWik"和调整后的参数始f"，记"，皆1)，令k = k+i，返回到 步骤(7b)，其中，Wik"为迭代第k+1次后的刚性配准膜腺图像序列，掉"为迭代第k+1次后的 水平变量，"为迭代第k+1次后的垂直变量，知"为迭代第k+i次后的顺时针旋转变量. (7e)根据最优参数(pi，p2，P3)按下列公式得到刚性转换公式Ti，根据Τι将W转换成刚性 配准后的膜腺图像序列Wl:6. 根据权利要求1所述的基于稀疏低秩与Atlas集的3D MRI膜腺分割方法，其中步骤 (8)中使用基于Demons的弹性配准算法对化和Wi进行配准，得到弹性配准后的膜腺图像序列 W2和弹性转换公式T2，按如下步骤进行： (8a)设初始状态k=l，在参考图像化上选择所有像素点作为Demons点，听=巧，给定允 许误差ε，其中e>0，k为迭代次数，其中听为第一次迭代的弹性配准膜腺图像序列； (8b)根据Symmetric Demons公式计算Wf中所有像素点的形变向量；其中，U表示形变向量，uk表示迭代第k次后的形变向量，表示迭代第k次后弹性配准 膜腺图像序列，V表示求图像的梯度，II表示取模操作，η表示归一化因子； (8c)根据形变向量uk按下列公式求出弹性转换公式巧，根据巧将巧转换为跨4,其中， 巧为迭代第k次后的弹性转换公式，为迭代第k+1次后的弹性配准膜腺图像序列： 其中，ξ表示项力衰减系数；(8d)计算化与W/"的误差平方和，如果该误差平方和小于ε，输出弹性转换公式Γ;二巧 和弹性配准后的膜腺图像序列馬=賊"，结束迭代，否则k = k+l，返回步骤(8b)。7. 根据权利要求1所述的基于稀疏低秩与Atlas集的3D MRI膜腺分割方法，步骤(9)取 相关性序列最高的前五个系数对应的标签，并对运五个标签进行融合，得到参考图像化精 分割后的结果Ri，按如下步骤进行： (9a)根据相关性序列最高的前五个系数gi~g5和他们对应的标签图像Gi~Gs，计算标 签融合后的结果R/ :其中j的取值范围为1~5, Σ表示求和操作； (9b)按下列公式对(9a)中标签融合后的结果R/进行修正，得到精分割后的结果Ri:
【文档编号】G06T7/00GK106097374SQ201610473797
【公开日】2016年11月9日
【申请日】2016年6月24日
【发明人】缑水平, 牛帅, 姚瑶, 王蓉芳, 焦李成, 盛珂
【申请人】西安电子科技大学