CN106097374A - 基于稀疏低秩与Atlas集的3D MRI胰腺分割方法 - Google Patents
基于稀疏低秩与Atlas集的3D MRI胰腺分割方法 Download PDFInfo
- Publication number
- CN106097374A CN106097374A CN201610473797.6A CN201610473797A CN106097374A CN 106097374 A CN106097374 A CN 106097374A CN 201610473797 A CN201610473797 A CN 201610473797A CN 106097374 A CN106097374 A CN 106097374A
- Authority
- CN
- China
- Prior art keywords
- pancreas
- sequence
- image sequence
- image
- iteration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30096—Tumor; Lesion
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明提出一种基于稀疏低秩与Atlas集的3D MRI胰腺分割方法,主要解决现有技术中分割胰腺粘连现象严重、分割结果不精确的问题。其步骤如下:1)输入3D MRI数据并归一化处理;2)对归一化处理后的序列图像进行稀疏低秩分解并进行Hessian矩阵增强;3)在增强后的图像上人机交互分割出胰腺标签;4)构造粗分割集;5)输入Atlas集,依次进行刚性配准和弹性配准,并进行标签融合,输出分割结果。本发明运用矩阵增强、矩阵分解、图集法,并利用3D MRI图像的空间信息,进行胰腺分割,减少了胰腺目标与其他组织器官的粘连现象,图像轮廓清晰,提高了分割的准确度,可用于胰腺癌精确放射治疗中的胰腺目标检测。
Description
技术领域
本发明属于医学图像处理技术领域,更进一步涉及三维核磁共振图像3D MRI胰腺分割,可用于胰腺癌精确放射治疗中的胰腺目标检测。
背景技术
胰腺癌是出现在胰脏的癌症,通常认为胰腺癌是常见肿瘤中恶性程度最高、死亡率最高的癌症,约90%的患者无法以手术根治,五年存活率低于5%。根据世界卫生组织估计,胰腺癌在全世界发生率排第十三,死亡率排第八。随着科技的发展,其他癌症的治疗都有了提高,但胰腺癌却没有明显的突破,死亡率也没有明显的下降趋势,所以必须加强研究,早日攻克胰腺癌对人类的威胁。
传统检查胰腺癌的方式是医生根据经验和观察直接阅片,定性的判断患者的病情,这样易受情绪影响,且随着工作量增加,阅片准确率下降。对胰腺癌的治疗方法中,放射治疗是一种主要的方法,他试图用射线杀死肿瘤,但也会对正常组织造成一次性或永久性伤害。精确放疗通过精确的肿瘤定位、精确的计划设计、精确的剂量计算及在治疗机上精确执行的一种全新的肿瘤放疗技术,有助于减少放射治疗中对健康组织放射造成的伤害。精确分割胰腺组织有助于帮助医生分析胰腺病情,也是精确放疗不可缺少的一步。
目前,对胰腺癌进行检测和治疗的设备主要是核磁共振成像MRI系统,该设备功能强大、技术含量高、软组织图像清晰,而且对人体检测损害小,已成为对胰腺进行分析、诊断和治疗的重要手段,而从MRI图像中精确的分割出胰腺,就是这些手段的重要前提。
传统医学图像分割方法包括基于区域的方法、基于边缘的方法、基于神经网络的方法、基于模糊聚类的方法、基于图谱的方法等。由于这些方法大多数是基于2D图像序列的分割,利用2D图像本身的特性信息完成感兴趣区域的分割,例如灰度、纹理等信息,而没有利用体素的空间信息,因此不能进行准确的胰腺提取;医学MRI图像有噪声强、对比度低、软组织边界不明显的特点,而这些方法大多数对图像噪声和对比度敏感,所以难以定位目标边界,使得分割结果不精确。
发明内容
本发明的目的在于针对以上现有技术的不足,提出一种基于稀疏低秩与Atlas集的3D MRI胰腺分割方法,以提高3D MRI胰腺分割的精度。
实现本发明目的的技术方案如下:
(1)输入三维核磁共振3D MRI序列图像,对其进行归一化处理,得到归一化后的序列图像C;
(2)用交替迭代法对归一化处理后的序列图像C进行稀疏低秩矩阵分解,得到稀疏图像序列A和低秩图像序列B;
(3)对稀疏图像序列A进行Hessian矩阵增强,得到二值增强图像序列A1;
(4)通过一次人机交互,在二值增强图像序列A1中确定胰腺的位置区域P,根据胰腺位置区域P,利用形态学原理在二值增强图像序列A1中分割出胰腺的标签图像序列L;
(5)根据标签图像序列L和步骤(1)中的序列图像C,得到粗分割集Q,其中包含F帧图像;
(6)输入胰腺图像序列W及人工勾画出的胰腺标签序列E,即Atlas集;
(7)选择粗分割集Q中的第i帧图像作为参考图像,记为Qi,将Atlas集中的胰腺图像序列W作为浮动图像,用基于互信息的刚性配准算法对Qi、W进行配准,得到刚性配准后的胰腺图像序列W1和刚性转换公式T1,根据T1对胰腺标签序列E进行迁移,得到刚性配准后的胰腺标签序列E1,其中i的取值范围为1~F;
(8)将刚性配准后的胰腺图像序列W1作为浮动图像,使用基于Demons的弹性配准算法对参考图像Qi和该W1进行配准,得到弹性配准后的胰腺图像序列W2和弹性转换公式T2,根据T2对刚性配准后的胰腺标签序列E1进行迁移,得到弹性配准后的胰腺标签序列E2;
(9)计算参考图像Qi与弹性配准后的胰腺图像序列W2的相关性序列,根据该相关性序列及弹性配准后的胰腺图像序列W2和弹性配准后的胰腺标签序列E2这三者是一一对应的关系,取相关性序列最高的前五个系数对应的标签,并对这五个标签进行融合,得到参考图像Qi精分割后的结果Ri,其中i的取值范围为1~F;
(10)重复步骤(7)~(9),对粗分割集Q中所有图像精分割后,得到最终结果精分割集R。
本发明与现有技术相比具有以下优点:
1.本发明采用Hessian矩阵增强的方法,将三维核磁共振图像3D MRI的空间信息运用到胰腺目标的分割中,使分割结果更加准确。
2.本发明采用稀疏低秩矩阵分解方法,将三维核磁共振图像3D MRI分解成低秩序列部分和稀疏序列部分,通过对稀疏序列部分进行操作,有效减少目标与其他器官的粘连现象。
3.本发明运用Atlas集方法,将粗分割结果与胰腺模板库图像配准,再进行标签迁移和融合得到最终分割结果,使得分割目标区域能保持较好的轮廓。
附图说明
图1是本发明的实现流程图;
图2是待分割3D MRI图像序列中的一帧及该帧的标签;
图3是用本发明与现有level set法和graph cut法的分割结果对比图。
具体实施方式
参照图1,本发明的实现步骤如下:
步骤1、输入数据并归一化处理。
输入三维核磁共振图像3D MRI序列C0,按下列公式对其进行归一化处理,得到归一化处理后的序列图像C:
其中,max()为提取矩阵最大元素。
步骤2、对归一化处理后的序列图像进行稀疏低秩分解。
用交替迭代法对归一化处理后的序列图像C进行稀疏低秩矩阵分解,得到稀疏图像序列A和低秩图像序列B:
(2a)设置初始值,稀疏图像序列A为零矩阵,低秩图像序列B=C;
(2b)通过公式:更新稀疏矩阵A,其中,是欧几里得投影,γ为常数,且γ>0,β指偏离线性约束的惩罚参数,且β>0,k为迭代次数,Bk为第k次迭代后的低秩矩阵,Zk为第k次迭代后的约束矩阵;
(2c)对进行奇异值分解产生三个不同的矩阵:左奇异向量矩阵Uk+1、右奇异向量矩阵Vk+1及对角奇异值矩阵
(2d)通过公式:更新低秩矩阵B,其中,T为矩阵的转置,max{}为提取最大元素,diag()为提取对角线元素;
(2e)通过公式Zk+1=Zk-β(Ak+1+Bk+1-C)更新约束性矩阵Z,其中,Ak+1是第k+1次迭代后的稀疏矩阵,Zk+1是第k+1次迭代后的线性约束矩阵,Bk+1是第k+1次迭代后的低秩矩阵;
(2f)设定迭代次数k=1000,重复步骤(2a)~(2e),迭代结束后,得到稀疏图像序列A和低秩图像序列B。
步骤3、对稀疏图像序列进行Hessian矩阵增强。
对稀疏图像序列A进行Hessian矩阵增强,得到二值增强图像序列A1:
(3a)根据下面公式,求出稀疏图像序列A中每个像素点I的Hessian矩阵H:
其中,I表示稀疏图像序列A中每一个像素点,Ixx代表沿x方向的二阶偏微分,Ixy代表先沿x方向再沿y方向的二阶偏微分,Ixz代表先沿x方向再沿z方向的二阶偏微分,Iyx代表先沿y方向再沿x方向的二阶偏微分,Iyy代表沿y方向的二阶偏微分,Iyz代表先沿y方向再沿z方向的二阶偏微分,Izx代表先沿z方向再沿x方向的二阶偏微分,Izy代表先沿z方向再沿y方向的二阶偏微分,Izz代表沿z方向的二阶偏微分;
(3b)根据Cholesky分解法,计算步骤(3a)中每一个Hessian矩阵H的三个特征值,并按绝对值从小到大顺序排序为λ1,λ2,λ3;
(3c)根据下面公式和每一个像素点I的Hessian矩阵特征值,对稀疏图像序列A的每一个像素点进行增强;
其中,α是起调节面状权重RA作用的常数,α>0,τ是起调节管状权重RB作用的常数,τ>0,c是起调节球状权重S作用的常数,c>0, v1(λ)是指每个像素点I增强后的结果。
步骤4、人机交互分割出胰腺标签。
利用matlab上的roipoly函数,在步骤3中得到的二值增强图像序列A1上手动选择胰腺的位置区域P,利用形态学原理,将A1中与P组成连通域的部分提取出来,得到胰腺的标签图像序列L:
(4a)对二值增强图像序列A1进行形态学腐蚀操作,得到腐蚀后的二值增强图像序列A1′:
(4b)根据胰腺位置区域P从腐蚀后的二值增强图像序列A1中提取出连通区域L′;
(4c)对提取出来的连通区域L′进行形态学膨胀操作,得到标签图像序列L。
步骤5、构造粗分割集。
根据标签图像序列L和步骤(1)中的序列图像C,按下列公式,得到包含F帧图像的粗分割集Q:
Q=C.*L
其中,.*表示矩阵对应元素相乘。
步骤6、输入Atlas集。
输入胰腺图像序列W及人工勾画出的胰腺标签序列E,W和E构成Atlas集。
步骤7、基于互信息的刚性配准。
(7a)选择粗分割集Q中的第i帧图像作为参考图像,记为Qi,将Atlas集中的胰腺图像序列W作为浮动图像,用基于互信息的刚性配准算法对Qi、W进行配准,得到刚性配准后的胰腺图像序列W1和刚性转换公式T1:
(7a1)设初始状态k=0,初始参数为(0,0,0),W1 0=W,其中,p1代表水平变量、p2代表垂直变量、p3代表顺时针旋转变量,k为迭代次数,W1 0为迭代第0次后的刚性配准胰腺图像序列;
(7a2)用PV插值法统计Qi和W1 k的联合直方图,根据该联合直方图计算Qi和W1 k的互信息值,其中,W1 k为迭代第k次后的刚性配准胰腺图像序列;
(7a3)用Powell算法根据步骤(7a2)中得到的互信息值判断参数是否最优,如果不是最优,则进行步骤(7a4),如果是最优,则停止迭代,输出最优参数进行步骤(7a5),其中,为迭代第k次后的水平变量,为迭代第k次后的垂直变量,为迭代第k次后的顺时针旋转变量;
(7a4)对W1 k进行平移、旋转,得到W1 k+1和调整后的参数令k=k+1,返回到步骤(7a2),其中,W1 k+1为迭代第k+1次后的刚性配准胰腺图像序列,为迭代第k+1次后的水平变量,为迭代第k+1次后的垂直变量,为迭代第k+1次后的顺时针旋转变量;
(7a5)根据最优参数(p1,p2,p3)按下列公式得到刚性转换公式T1,根据T1将W转换成刚性配准后的胰腺图像序列W1:
(7b)根据T1对胰腺标签序列E进行迁移,得到刚性配准后的胰腺标签序列E1。
步骤8、基于Demons的非刚性配准。
(8a)将刚性配准后的胰腺图像序列W1作为浮动图像,使用基于Demons的弹性配准算法对参考图像Qi和该W1进行配准,得到弹性配准后的胰腺图像序列W2和弹性转换公式T2:
(8a1)设初始状态k=1,在参考图像Qi上选择所有像素点作为Demons点,给定允许误差ε,其中ε>0,k为迭代次数,其中为第一次迭代的弹性配准胰腺图像序列;
(8a2)根据Symmetric Demons公式计算中所有像素点的形变向量:
其中,u表示形变向量,uk表示迭代第k次后的形变向量,表示迭代第k次后弹性配准胰腺图像序列,表示求图像的梯度,||表示取模操作,η表示归一化因子;
(8a3)根据形变向量uk按下列公式求出弹性转换公式根据将转换为其中,为迭代第k次后的弹性转换公式,为迭代第k+1次后的弹性配准胰腺图像序列:
其中,ξ表示项力衰减系数;
(8a4)计算Qi与的误差平方和,如果该误差平方和小于ε,输出弹性转换公式和弹性配准后的胰腺图像序列结束迭代,否则k=k+1,返回步骤(8a2);
(8b)根据T2对刚性配准后的胰腺标签序列E1进行迁移,得到弹性配准后的胰腺标签序列E2。
步骤9、进行标签融合。
计算参考图像Qi与弹性配准后的胰腺图像序列W2的相关性序列,根据该相关性序列及弹性配准后的胰腺图像序列W2和弹性配准后的胰腺标签序列E2这三者是一一对应的关系,取相关性序列最高的前五个系数对应的标签,并对这五个标签进行融合,得到参考图像Qi精分割后的结果Ri。
(9a)根据相关性序列最高的前五个系数g1~g5和他们对应的标签图像G1~G5,计算标签融合后的结果Ri′:
其中j的取值范围为1~5,∑表示求和操作;
(9b)按下列公式对(9a)中标签融合后的结果Ri′进行修正,得到精分割后的结果Ri:
步骤10、重复步骤(7)~(9),对粗分割集Q中所有图像精分割,得到最终结果精分割集R。
以下结合仿真图对本发明的效果做进一步的说明。
1.仿真条件:
硬件平台为:CPU为Inter Core i7-4790K,主频为4GHz,内存为32GB;
软件平台为:Windows 10专业版64位操作系统,MatlabR2014a。
2.仿真数据:
本例中待分割的数据采用的是通过三维核磁共振仪获取的一名志愿者腹部三维核磁共振图像序列,图像大小为320*284,一共有5帧图像,通过人工勾画出胰腺作为标签,如图2所示,其中,图2a是胰腺图像,图2b是胰腺图像的标签。Atlas集采用的是其他志愿者通过三维核磁共振仪获取的三维核磁共振图像,它们对应的标签图像通过人工勾画获得。
3.仿真内容:
为了验证算法的有效性,一般在进行图像分割以后,会以标签为参考图像计算指标,豪斯多夫距离、平均绝对表面距离、相似性、灵敏度、特异性五个指标作为分割精度的判断标准,豪斯多夫距离、平均绝对表面距离这两个指标数值越大分割结果越精确,相似性、灵敏度、特异性这三个指标数值越小分割结果越精确。
用本发明和现有的两个具有代表性的图像分割方法level set法和graph cut法进行仿真实验对比。
仿真一:分别通过这三种方法对上述仿真数据进行胰腺分割仿真,对三种方法得到的分割结果进行三维重建,结果如图3所示,其中图3a为本发明的分割结果,图3b为levelset法的分割结果,图3c为graph cut法的分割结果,可以看出用本发明分割出的胰腺轮廓更好。
仿真二:分别通过这三种方法对上述仿真数据进行胰腺分割仿真,计算分割结果的五个指标:豪斯多夫距离、平均绝对表面距离、相似性、灵敏度、特异性,结果如表1所示:
表1用本发明与现有level set法和graph cut法的分割指标对比表
豪斯多夫距离 | 平均绝对表面距离 | 相似性 | 灵敏性 | 特异性 | |
本发明 | 18.0482mm | 1.6165mm | 0.7964 | 0.7632 | 0.9966 |
Level set | 32.3503mm | 3.9333mm | 0.5721 | 0.4726 | 0.996 |
Graph cut | 69.7327mm | 12.3378mm | 0.4476 | 0.6459 | 0.9722 |
从表1可以看出本发明的分割结果在指标上优于其他两种算法,说明本发明分割胰腺结果精度更高。
Claims (7)
1.一种基于稀疏低秩与Atlas集的3D MRI胰腺分割方法,包括:
(1)输入三维核磁共振3D MRI序列图像,对其进行归一化处理,得到归一化后的序列图像C;
(2)用交替迭代法对归一化处理后的序列图像C进行稀疏低秩矩阵分解,得到稀疏图像序列A和低秩图像序列B;
(3)对稀疏图像序列A进行Hessian矩阵增强,得到二值增强图像序列A1;
(4)通过一次人机交互,在二值增强图像序列A1中确定胰腺的位置区域P,根据胰腺位置区域P,利用形态学原理在二值增强图像序列A1中分割出胰腺的标签图像序列L;
(5)根据标签图像序列L和步骤(1)中的序列图像C,得到粗分割集Q,其中包含F帧图像;
(6)输入胰腺图像序列W及人工勾画出的胰腺标签序列E,即Atlas集;
(7)选择粗分割集Q中的第i帧图像作为参考图像,记为Qi,将Atlas集中的胰腺图像序列W作为浮动图像,用基于互信息的刚性配准算法对Qi、W进行配准,得到刚性配准后的胰腺图像序列W1和刚性转换公式T1,根据T1对胰腺标签序列E进行迁移,得到刚性配准后的胰腺标签序列E1,其中i的取值范围为1~F;
(8)将刚性配准后的胰腺图像序列W1作为浮动图像,使用基于Demons的弹性配准算法对参考图像Qi和该W1进行配准,得到弹性配准后的胰腺图像序列W2和弹性转换公式T2,根据T2对刚性配准后的胰腺标签序列E1进行迁移,得到弹性配准后的胰腺标签序列E2;
(9)计算参考图像Qi与弹性配准后的胰腺图像序列W2的相关性序列,根据该相关性序列及弹性配准后的胰腺图像序列W2和弹性配准后的胰腺标签序列E2这三者是一一 对应的关系,取相关性序列最高的前五个系数对应的标签,并对这五个标签进行融合,得到参考图像Qi精分割后的结果Ri,其中i的取值范围为1~F;
(10)重复步骤(7)~(9),对粗分割集Q中所有图像精分割后,得到最终结果精分割集R。
2.根据权利要求1所述的基于稀疏低秩与Atlas集的3D MRI胰腺分割方法,其中步骤(2)用交替迭代法对归一化处理后的序列图像C进行稀疏低秩矩阵分解,得到稀疏图像序列A和低秩图像序列B,按如下步骤进行:
(2a)设置初始值,稀疏图像序列A为零矩阵,低秩图像序列B=C;
(2b)通过公式:更新稀疏矩阵A,其中, 是欧几里得投影,γ为常数,且γ>0,β指偏离线性约束的惩罚参数,且β>0,k为迭代次数,Bk为第k次迭代后的低秩矩阵,Zk为第k次迭代后的约束矩阵;
(2c)对进行奇异值分解产生三个不同的矩阵:左奇异向量矩阵Uk+1、右奇异向量矩阵Vk+1及对角奇异值矩阵
(2d)通过公式:更新低秩矩阵B,其中,T为矩阵的转置,max{}为提取最大元素,diag()为提取对角线元素;
(2e)通过公式Zk+1=Zk-β(Ak+1+Bk+1-C)更新约束性矩阵Z,其中,Ak+1是第k+1次迭代后的稀疏矩阵,Zk+1是第k+1次迭代后的线性约束矩阵,Bk+1是第k+1次迭代后的低秩矩阵;
(2f)设定迭代次数k的值,重复步骤(2a)~(2e),迭代结束后,得到稀疏图像序列A和低秩图像序列B。
3.根据权利要求1所述的基于稀疏低秩与Atlas集的3D MRI胰腺分割方法,其中 步骤(3)对稀疏图像序列A进行Hessian矩阵增强,得到二值增强图像序列A1,按如下步骤进行:
(3a)根据下面公式,求出稀疏图像序列A中每个像素点I的Hessian矩阵H:
其中,I表示稀疏图像序列A中每一个像素点,Ixx代表沿x方向的二阶偏微分,Ixy代表先沿x方向再沿y方向的二阶偏微分,Ixz代表先沿x方向再沿z方向的二阶偏微分,Iyx代表先沿y方向再沿x方向的二阶偏微分,Iyy代表沿y方向的二阶偏微分,Iyz代表先沿y方向再沿z方向的二阶偏微分,Izx代表先沿z方向再沿x方向的二阶偏微分,Izy代表先沿z方向再沿y方向的二阶偏微分,Izz代表沿z方向的二阶偏微分;
(3b)根据Cholesky分解法,计算步骤(3a)中每一个Hessian矩阵H的三个特征值,并按绝对值从小到大顺序排序为λ1,λ2,λ3;
(3c)根据下面公式和每一个像素点I的Hessian矩阵特征值,对稀疏图像序列A的每一个像素点进行增强;
其中,α是起调节面状权重RA作用的常数,α>0,τ是起调节管状权重RB作用的常数,τ>0,c是起调节球状权重S作用的常数,c>0, v1(λ)是指每个像素点I增强后的结果。
4.根据权利要求1所述的基于稀疏低秩与Atlas集的3D MRI胰腺分割方法,其中步骤(5)根据标签图像序列L和步骤(1)中的序列图像C,得到粗分割集Q,是将序列图像C与标签图像序列L相乘,即:
Q=C.*L
其中,.*表示矩阵对应元素相乘。
5.根据权利要求1所述的基于稀疏低秩与Atlas集的3D MRI胰腺分割方法,其中步骤(7)中用基于互信息的刚性配准算法对Qi、W进行配准,得到刚性配准后的胰腺图像序列W1和刚性转换公式T1,按如下步骤进行:
(7a)设初始状态k=0,初始参数为(0,0,0),W1 0=W,其中,p1代表水平变量、p2代表垂直变量、p3代表顺时针旋转变量,k为迭代次数,W1 0为迭代第0次后的刚性配准胰腺图像序列;
(7b)用PV插值法统计Qi和W1 k的联合直方图,根据该联合直方图计算Qi和W1 k的互信息值,其中,W1 k为迭代第k次后的刚性配准胰腺图像序列;
(7c)用Powell算法根据步骤(7b)中得到的互信息值判断参数是否最优,如果不是最优,则进行步骤(7d),如果是最优,则停止迭代,输出最优参数 进行步骤(7e),其中,为迭代第k次后的水平变量,为迭代第k次后的垂直变量,为迭代第k次后的顺时针旋转变量;
(7d)对W1 k进行平移、旋转,得到W1 k+1和调整后的参数令k=k+1,返回到步骤(7b),其中,W1 k+1为迭代第k+1次后的刚性配准胰腺图像序列,为迭代第k+1次后的水平变量,为迭代第k+1次后的垂直变量,为迭代第k+1次后的顺时针旋转变量;
(7e)根据最优参数(p1,p2,p3)按下列公式得到刚性转换公式T1,根据T1将W转换成刚性配准后的胰腺图像序列W1:
6.根据权利要求1所述的基于稀疏低秩与Atlas集的3D MRI胰腺分割方法,其中步骤(8)中使用基于Demons的弹性配准算法对Qi和W1进行配准,得到弹性配准后的胰腺图像序列W2和弹性转换公式T2,按如下步骤进行:
(8a)设初始状态k=1,在参考图像Qi上选择所有像素点作为Demons点,给定允许误差ε,其中ε>0,k为迭代次数,其中为第一次迭代的弹性配准胰腺图像序列;
(8b)根据Symmetric Demons公式计算中所有像素点的形变向量:
其中,u表示形变向量,uk表示迭代第k次后的形变向量,表示迭代第k次后弹性配准胰腺图像序列,表示求图像的梯度,||表示取模操作,η表示归一化因子;
(8c)根据形变向量uk按下列公式求出弹性转换公式根据将转换为其中,为迭代第k次后的弹性转换公式,为迭代第k+1次后的弹性配准胰腺图像序列:
其中,ξ表示项力衰减系数;
(8d)计算Qi与的误差平方和,如果该误差平方和小于ε,输出弹性转换公式 和弹性配准后的胰腺图像序列结束迭代,否则k=k+1,返回步骤(8b)。
7.根据权利要求1所述的基于稀疏低秩与Atlas集的3D MRI胰腺分割方法,步骤(9)取相关性序列最高的前五个系数对应的标签,并对这五个标签进行融合,得到参考图像Qi精分割后的结果Ri,按如下步骤进行:
(9a)根据相关性序列最高的前五个系数g1~g5和他们对应的标签图像G1~G5,计 算标签融合后的结果Ri′:
其中j的取值范围为1~5,∑表示求和操作;
(9b)按下列公式对(9a)中标签融合后的结果Ri′进行修正,得到精分割后的结果Ri:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610473797.6A CN106097374B (zh) | 2016-06-24 | 2016-06-24 | 基于稀疏低秩与Atlas集的3D MRI胰腺分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610473797.6A CN106097374B (zh) | 2016-06-24 | 2016-06-24 | 基于稀疏低秩与Atlas集的3D MRI胰腺分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106097374A true CN106097374A (zh) | 2016-11-09 |
CN106097374B CN106097374B (zh) | 2018-11-30 |
Family
ID=57253668
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610473797.6A Active CN106097374B (zh) | 2016-06-24 | 2016-06-24 | 基于稀疏低秩与Atlas集的3D MRI胰腺分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106097374B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106803256A (zh) * | 2017-01-13 | 2017-06-06 | 深圳市唯特视科技有限公司 | 一种基于投影卷积网络的三维形状分割及语义标记方法 |
CN106952273A (zh) * | 2017-03-09 | 2017-07-14 | 上海联影医疗科技有限公司 | 医学图像中胰腺的分割方法及装置 |
CN107016395A (zh) * | 2017-03-18 | 2017-08-04 | 复旦大学 | 基于稀疏表示体系的原发性脑部淋巴瘤和胶质母细胞瘤的鉴别方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102622749A (zh) * | 2012-02-22 | 2012-08-01 | 中国科学院自动化研究所 | 三维磁共振图像脑子结构自动分割的方法 |
CN103150739A (zh) * | 2013-03-04 | 2013-06-12 | 上海大学 | 一种基于多特征稳健主成分分析的视频运动对象分割算法 |
JP2013143102A (ja) * | 2012-01-12 | 2013-07-22 | Nikon Corp | 移動物体検出装置、移動物体検出方法およびプログラム |
CN105335975A (zh) * | 2015-10-22 | 2016-02-17 | 西安电子科技大学 | 基于低秩分解和直方图统计的极化sar图像分割方法 |
-
2016
- 2016-06-24 CN CN201610473797.6A patent/CN106097374B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013143102A (ja) * | 2012-01-12 | 2013-07-22 | Nikon Corp | 移動物体検出装置、移動物体検出方法およびプログラム |
CN102622749A (zh) * | 2012-02-22 | 2012-08-01 | 中国科学院自动化研究所 | 三维磁共振图像脑子结构自动分割的方法 |
CN103150739A (zh) * | 2013-03-04 | 2013-06-12 | 上海大学 | 一种基于多特征稳健主成分分析的视频运动对象分割算法 |
CN105335975A (zh) * | 2015-10-22 | 2016-02-17 | 西安电子科技大学 | 基于低秩分解和直方图统计的极化sar图像分割方法 |
Non-Patent Citations (2)
Title |
---|
ZHOUBINGXU 等: "efficient multi-atlas abdominal segmentation on clinically acquired CT with SIMPLE context learning", 《MEDICAL IMAGE ANALYSIS》 * |
李涛: "基于稀疏理论的图像分割方法", 《万方在线论文》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106803256A (zh) * | 2017-01-13 | 2017-06-06 | 深圳市唯特视科技有限公司 | 一种基于投影卷积网络的三维形状分割及语义标记方法 |
CN106952273A (zh) * | 2017-03-09 | 2017-07-14 | 上海联影医疗科技有限公司 | 医学图像中胰腺的分割方法及装置 |
CN106952273B (zh) * | 2017-03-09 | 2019-11-19 | 上海联影医疗科技有限公司 | 医学图像中胰腺的分割方法及装置 |
CN107016395A (zh) * | 2017-03-18 | 2017-08-04 | 复旦大学 | 基于稀疏表示体系的原发性脑部淋巴瘤和胶质母细胞瘤的鉴别方法 |
CN107016395B (zh) * | 2017-03-18 | 2021-02-26 | 复旦大学 | 稀疏表示的原发性脑部淋巴瘤和胶质母细胞瘤的鉴别系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106097374B (zh) | 2018-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | CT synthesis from MRI using multi-cycle GAN for head-and-neck radiation therapy | |
Sjölund et al. | Generating patient specific pseudo-CT of the head from MR using atlas-based regression | |
Anders et al. | Performance of an atlas-based autosegmentation software for delineation of target volumes for radiotherapy of breast and anorectal cancer | |
CN102436584B (zh) | 基于字典学习的胃部ct图像感兴趣区域检测系统 | |
CN108629785B (zh) | 基于自步学习的三维磁共振胰腺图像分割方法 | |
CN107403201A (zh) | 肿瘤放射治疗靶区和危及器官智能化、自动化勾画方法 | |
Lei et al. | 4D-CT deformable image registration using an unsupervised deep convolutional neural network | |
AU2016339009A1 (en) | Pseudo-CT generation from MR data using tissue parameter estimation | |
Palma et al. | A voxel-based approach to explore local dose differences associated with radiation-induced lung damage | |
CN102622759B (zh) | 一种结合灰度与几何信息的医学图像配准方法 | |
Nitzken et al. | Improving full-cardiac cycle strain estimation from tagged CMR by accurate modeling of 3D image appearance characteristics | |
CN102999917B (zh) | 基于t2-mri和dw-mri的宫颈癌图像自动分割方法 | |
Paganelli et al. | Time-resolved volumetric MRI in MRI-guided radiotherapy: an in silico comparative analysis | |
CN109859817A (zh) | 食管放射治疗计划中风险器官剂量学评估方法、评估系统 | |
CN106097374A (zh) | 基于稀疏低秩与Atlas集的3D MRI胰腺分割方法 | |
Tran et al. | Evaluation of MRI-derived surrogate signals to model respiratory motion | |
CN106898000A (zh) | 一种磁共振成像大脑灰质核团的自动分割方法 | |
CN104166994B (zh) | 一种基于训练样本优化的骨骼抑制方法 | |
Xiao et al. | RefineNet‐based 2D and 3D automatic segmentations for clinical target volume and organs at risks for patients with cervical cancer in postoperative radiotherapy | |
ZongRen et al. | DenseTrans: multimodal brain tumor segmentation using swin transformer | |
CN104331864B (zh) | 基于非下采样轮廓波和视觉显著模型的乳腺影像处理 | |
Zhang et al. | Enhancing liver tumor localization accuracy by prior-knowledge-guided motion modeling and a biomechanical model | |
Nasief et al. | Predicting necessity of daily online adaptive replanning based on wavelet image features for MRI guided adaptive radiation therapy | |
Gooya et al. | An EM algorithm for brain tumor image registration: A tumor growth modeling based approach | |
Iyer et al. | Deep learning-based auto-segmentation of swallowing and chewing structures |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |