申请/专利权人：上海联影医疗科技股份有限公司

申请日：2015-12-22

公开（公告）日：2020-11-27

公开（公告）号：CN106910194B

主分类号：G06T7/11(20170101)

分类号：G06T7/11(20170101)

优先权：

专利状态码：有效-授权

法律状态：2020.11.27#授权;2017.10.17#实质审查的生效;2017.06.30#公开

摘要：本发明公开了一种心脏心功能磁共振图像中的心内外膜的分割方法,其包括:在包含左心室心肌的若干片层、且在不同心动期相的心脏磁共振图像INP中定位舒张末期，并获得舒张末期的N个片层的磁共振图像的血池区域的粗分割结果；基于极坐标转换方法，利用射线扫描将舒张末期的每个片层的磁共振图像中左心室感兴趣区域内的图像数据转换为二维极坐标转换图像；基于双动态规划方法在二维极坐标转换图像中检测心内膜及心外膜；利用极坐标逆转换获取原始片层图像中的心内外膜，求凸包并平滑，完成舒张末期的N个片层的磁共振图像中的心内外膜的分割；从舒张末期的磁共振图像中的心内外膜分割结果向其它心动期相的磁共振图像中的心内外膜衍生分割结果。

主权项：1.一种心脏心功能磁共振图像中的心内外膜的分割方法,其特征在于包括以下步骤:S1.获取包含左心室心肌的若干片层、且在不同心动期相的心脏磁共振图像INP，其中N表示片层的序号，P表示心动期相序号，N、P均为大于或等于1的整数；S2.确定舒张末期的期相；S3.获得舒张末期的N个片层的磁共振图像的血池区域的粗分割结果；S4.以粗分割血池的质心为扫描起点，基于极坐标转换方法，利用射线扫描将舒张末期的每个片层的磁共振图像中左心室感兴趣区域内的图像数据转换为二维极坐标转换图像；S5.基于双动态规划方法在二维极坐标转换图像中检测心内膜及心外膜；利用极坐标逆转换获取原始片层图像中的心内、外膜，求凸包并平滑，完成舒张末期的N个片层的磁共振图像中的心内、外膜的分割；S6.从舒张末期的磁共振图像中的心内、外膜分割结果向其它心动期相的磁共振图像中的心内、外膜衍生分割结果；所述步骤S6包括：选取第K个片层的第t个心动期相的磁共振图像IKt；以第K个片层的第t-1个心动期相检测的心内膜确定的血池区域的质心为扫描起点，在第t-1个心动期相的磁共振图像中分割的心内、外膜边界附近，将第t个心动期相对应的磁共振图像IKt转换成二维极坐标转换图像PKt；统计从舒张末期至第t-1个心动期相的磁共振图像中已分割的心肌灰度的统计分布，通过原始图像IKt中像素点的灰度值与其在心肌灰度统计分布中相应灰度值的分布百分比值相乘，增强原始图像IKt中的心肌部分，并将增强的原始图像IKt转换为二维极坐标转换图像Pcard；基于双动态规划方法在二维极坐标转换图像Pkt、Pcard中检测左心室心内、外膜；第t-1个心动期相的分割的左心室心内外膜位置限定第t个心动期相的左心室心内外膜搜索的范围；基于二维极坐标转换图像中检测到的左心室心内、外膜，利用极坐标逆转换获取原始图像IKt中的左心室心内外膜，求凸包并平滑，得到检测结果。

全文数据：心脏心功能磁共振图像中的心内外膜的分割方法【技术领域】[0001]本发明涉及医学图像领域，尤其是涉及心脏心功能磁共振图像的分割。【背景技术】[0002]近年来，心血管疾病的发病率和死亡率正逐年增加。每年在美国有近600万人被检查出患有心血管疾病，而我国则有近300万人死于心血管疾病。心脏功能是诊断心血管疾病的重要指标，可以通过分析四维心功能磁共振图像中左心室相关的全局和局部特征参数实现，如：心室体积、射血分数和心肌壁厚等。心功能特征参数提取的准确性依赖于左心室心内外膜的准确定位，而计算机辅助诊断方案有助于提高在四维心功能磁共振图像中定位左心室心内外膜的可靠性和效率。[0003]由于四维心功能磁共振图像在不同病患扫描图像、相同病患的不同扫描断层图像以及相同扫描断层图像中灰度和梯度分布的巨大差异，另外，左心室心外膜周边临近组织的多样性，这些都增加了左心室心内外膜检测的难度。虽然在过去的近20年中，已有很多研究方法和结果的报道，但是在四维心功能磁共振影像中对左心室心内外膜进行准确的检测仍然没有解决。针对四维心功能磁共振图像中左心室心内外膜的检测，手动和需用户辅助的半自动检测方法比较繁复、效率低、并且存在较大的观察者之间和观察者自身的差异。在全自动检测方面，由于四维心功能磁共振图像在不同病患扫描图像、相同病患的不同扫描断层图像以及相同扫描断层图像中灰度和梯度分布的巨大差异，另外，左心室心外膜周边临近组织的多样性，导致简单的基于阈值的检测方法或者基于梯度的检测方法不足以准确地提取出左心室心内外膜。基于模型的检测方法需要大量的训练样本去训练模型以获得较为准确的检测结果，而这些精准训练样本的获取则是一项更艰巨和繁复的任务。【发明内容】[0004]本发明所要解决的技术问题是提供一种心脏心功能磁共振图像中的心内外膜的分割方法。[0005]本发明为解决上述技术问题而采用的技术方案是:一种心脏心功能磁共振图像中的心内外膜的分割方法，其包括以下步骤：[0006]SI.获取包含左心室心肌的若干片层、且在不同心动期相的心脏磁共振图像INP，其中N表示片层的序号，P表示心动期相序号，N、P均为大于或等于1的整数；[0007]S2.确定舒张末期的期相；[0008]S3.获得舒张末期的N个片层的磁共振图像的血池区域的粗分割结果；[0009]S4.以粗分割血池的质心为扫描起点，基于极坐标转换方法，利用射线扫描将舒张末期的每个片层的磁共振图像中左心室感兴趣区域内的图像数据转换为二维极坐标转换图像；[0010]S5.基于双动态规划方法在二维极坐标转换图像中检测心内膜及心外膜;利用极坐标逆转换获取原始片层图像中的心内、外膜，求凸包并平滑，完成舒张末期的N个片层的磁共振图像中的心内、外膜的分割；[0011]S6.从舒张末期的磁共振图像中的心内、外膜分割结果向其它心动期相的磁共振图像中的心内、外膜衍生分割结果。[0012]优选的，所述二维极坐标转换图像中的X轴表示了扫描射线的序号，y轴表示了射线上点距离扫描起点的距离;扫描的射线数目为180条，扫描起始半径为粗分割血池区域的边界位置减数个像素点，扫描终止半径为粗分割血池区域的边界位置加数个像素点。[0013]优选的，所述双边动态规划中，射线扫描方法顺序获得的射线二维极坐标转换图像的每一列被认为是阶段，射线上的点（二维极坐标转换图像每一列上的点）被认为是阶段上的候选点，从第一阶段到最后阶段二维极坐标转换图像的第一列到最后一列具有最小累积局部能量cost的路径被认为是最优路径，即检测到的左心室心内、外膜。[0014]优选的，所述局部能量cost包括内部能量COStint和外部能量costext，其中内部能量COStint决定了最优路径的平滑性，而外部能量COSte3xt决定了最优路径位于梯度大的位置。[0015]优选的，所述内部能量C0stin^以下公式获得：[0016][0017]为心内膜候选点在相邻图像列中在垂直方向上的归一化距离；[0018]为心外膜候选点在相邻图像列中在垂直方向上的归一化距离；[0019]为由心内外膜候选点决定的心肌厚度在相邻图像列（）中的归一化距离；[0020]为心内膜候选点在图像第i列中在垂直方向上的位置，为心内膜候选点在图像第i-Ι列中在垂直方向上的位置，为心外膜候选点在图像第i列中在垂直方向上的位置，为心外膜候选点在图像第i_l列中在垂直方向上的位置，^_为心内膜候选点在相邻图像列中在垂直方向上变化距离最大值3个像素），drange为心肌厚度在相邻图像列中的变化距离最大值5个像素）。[0021]优选的，所述外部能量为Coste3xt由以下公式获得：[0022][0023]为心内膜的归一化梯度，为心外膜的归一化梯度，心内膜梯度是基于当前层原始的二维极坐标转换图像获得，心外膜梯度是基于当前层增强的二维极坐标转换图像获得。[0024]优选的，所述步骤S6包括：[0025]选取第K个片层的第t个心动期相当前心动期相的磁共振图像IKt;[0026]以第K个片层的第t-1个心动期相前一心动期相检测的心内膜确定的血池区域的质心为扫描起点，在第t-Ι个心动期相的磁共振图像中分割的心内、外膜边界附近，将第t个心动期相对应的磁共振图像Iiit转换成二维极坐标转换图像Ρια;[0027]统计从舒张末期至第t_l个心动期相的磁共振图像中已分割的心肌灰度的统计分布，通过原始图像IKt中像素点的灰度值与其在心肌灰度统计分布中相应灰度值的分布百分比值相乘，增强原始图像IKt中的心肌部分，并将增强的原始图像IKt转换为二维极坐标转换图像Pcard;[0028]基于双动态规划方法在二维极坐标转换图像Pkt、Pcard中检测左心室心内、外膜;第t-Ι个心动期相前一心动期相的分割的左心室心内外膜位置限定第t个心动期相（当前心动期相的左心室心内外膜搜索的范围。[0029]基于二维极坐标转换图像中检测到的左心室心内、外膜，利用极坐标逆转换获取原始图像Ik冲的左心室心内外膜，求凸包并平滑，得到检测结果。[0030]优选的，所述第t个心动期相（当前心动期相的左心室心内外膜搜索的范围为：以第t-1个心动期相（前一心动期相）上分割的心内膜上20个像素点至外膜边界下10个像素点。[0031]优选的，所述步骤S2包括：[0032]在中间片层的每个心动期相所对应的磁共振图像Imp上，设定初始感兴趣区域，计算初始感兴趣区域在所有心动期相上的最大强度投影图像，聚类分割最大强度投影图像，选取圆度最大的区域A;[0033]计算中间片层的每个心动期相所对应的磁共振图像Imp在区域A中的图像平均灰度值，具有最大平均灰度值的磁共振图像Imp所在心动期相被定义为舒张末期。[0034]优选的，所述步骤S3包括:从中间片层开始向心尖和心底两个方向片层的磁共振图像分别进行分割，完成舒张末期的磁共振图像InpEd中的血池区域的粗分割。[0035]本发明对比现有技术有如下的有益效果:本发明的方法结合采用左心室心内外膜位置和梯度信息，利用双动态规划技术，同时实现心脏舒张期末期中左心室心内外膜的检测;本发明的方法利用相同片层在不同时序图像中的连续性特点，准确和有效地在四维心功能磁共振图像中实现左心室心内外膜的检测。本发明的方法可为后续的基于心脏心功能磁共振图像的计算机辅助诊断提供有力的支持和保障，并提高后续诊断的效率和准确性。【附图说明】[0036]图1为本明实施例中的一种心脏心功能磁共振图像中的心内外膜的分割方法；[0037]图2为本发明实施例的一种心脏心功能磁共振图像中舒张末期图像中血池的分割方法的流程图；[0038]图3a为中间片层磁共振图像在感兴趣区域中的最大强度投影图像；[0039]图3b为中间片层磁共振图像的在初始感兴趣区域中的图像；[0040]图3c为通过模糊C均值聚类方法对图3b进行聚类分割的结果；[0041]图3d为中间片层磁共振图像中血池区域分割的结果；[0042]图4a_4c通过射线扫描方法对过分割区域处理的过程示意图；[0043]图5为本明实施例中的一种心脏心功能磁共振图像中的心内外膜的分割结果。【具体实施方式】[0044]请参阅图1-5，本发明实施例中的一种心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于包括以下步骤：[0045]SI.获取包含左心室心肌的若干片层、且在不同心动期相的心脏磁共振图像Inp，其中N表示片层的序号，P表示心动期相序号，N、P均为大于或等于1的整数；[0046]S2.确定舒张末期的期相；[0047]S3.获得舒张末期的N个片层的磁共振图像的血池区域的粗分割结果；[0048]S4.以粗分割血池的质心为扫描起点，基于极坐标转换方法，利用射线扫描将舒张末期的每个片层的磁共振图像中左心室感兴趣区域内的图像数据转换为二维极坐标转换图像；[0049]S5.基于双动态规划方法在二维极坐标转换图像中检测心内膜及心外膜;利用极坐标逆转换获取原始片层图像中的心内、外膜，求凸包并平滑，完成舒张末期的N个片层的磁共振图像中的心内、外膜的分割；[0050]S6.从舒张末期的磁共振图像中的心内、外膜分割结果向其它心动期相P-1个心动期相的磁共振图像中的心内、外膜衍生分割结果。[0051]具体的，上述Sl中的每一片层、且在每个期相中的心脏磁共振图像是通过以下方式得到：[0052]通过采集多个心动周期的、对应相同心动期相中的磁共振数据线，将所采集到的磁共振数据线填充到一个K-SpaceK空间）中，通过傅里叶变换得到相应片层在相应心动期相的磁共振图像。通过设置较小的时间窗，例如20-50ms，并在多个心动周期采集数据，可减小心脏运动伪影对图像的影响，提高图像的质量。[0053]心动期相是指将自左心室的收缩期开始至舒张末期的完整心动周期T的划分为多个时间段phase，这些时间段可以为一系类连续的时间序列，也可以为间隔的时间序列。[0054]在心室包含心肌部分的区域沿心脏的长轴方向获取若干个片层（slice，并通过磁共振扫描获得每个片层、在每个心动期相中的心脏磁共振图像Inp，其中N表示相应片层的序号，P表示心动期相phase序号，N、P均为大于或等于1的整数;具体的，N的取值可以为3、5、6等，P的取值可以为3、5、8、10等;另外，定义第1个片层为中间片层，其中1为大于或等于2的整数。[0055]具体的，上述S2中确定舒张末期的期相，通过以下方法：[0056]在中间片层的每个心动期相所对应的磁共振图像Imp上，设定初始感兴趣区域，计算初始感兴趣区域在所有心动期相P个心动期相上的最大强度投影图像MIP，聚类分割最大强度投影图像MIP，选取圆度最大的区域A;优选的聚类分割方法为模糊C均值聚类方法；[0057]计算中间片层的每个心动期相所对应的磁共振图像Imp在区域A中的图像平均灰度值，具有最大平均灰度值的磁共振图像Imp所在心动期相被定义为舒张末期PED，并确定中间片层在舒张末期对应的磁共振图像IMPiED;[0058]从中间片层开始向心尖和心底两个方向片层的磁共振图像分别进行分割，完成Ν个片层的舒张末期的磁共振图像InpEd中的血池的粗分割。[0059]初始感兴趣区域ROIo的中心为原图像（即磁共振图像Imp的中心，半径为若干个像素点，半径的大小为60-120个像素点，优选半径大小为100个像素点。[0060]然后，计算初始感兴趣区域在所有心动期相P个心动期相上的最大强度投影图像MIP，通过模糊C均值聚类方法分割最大强度投影图像MIP，选取圆度最大的区域A;其中，聚类分类的数目可以为2、3或4等;圆度的定义为周长*周长4*PI*面积）。[0061]进一步的，获取左心室感兴趣区域，所述左心室感兴趣区域的中心为区域A的质心，所述左心室感兴趣区域的边长为区域A的长轴长度加数个像素点，具体可为10-30个像素点，优选为20个像素点。[0062]具体的，对中间片层的舒张末期的磁共振图像ImpEd中的血池的粗分割包括以下步骤：[0063]对中间片层的舒张末期的磁共振图像ImpEd在左心室感兴趣区域中的图像进行伽马校正GammaCorrection，获得校正后的图像Imped’；通过伽马校正，可增强血池与心肌的对比度，方便后续的识别和分割；[0064]通过模糊C均值聚类方法对校正后的图像ImpEd’进行分类，分类数目为2，选取最亮类作为候选区域Sm;[0065]在候选区域Sm中选取与区域A重叠率最大的区域作为中间片层的舒张末期的磁共振图像Imped中的血池区域BIoodPooI，BP的粗分割结果。[0066]进一步的，对第一个至第M-I个片层的舒张末期的磁共振图像的血池粗分害U，或第M+1个至第N个片层的舒张末期的磁共振图像的血池粗分割通过以下方式获得：以前一个片层的磁共振图像的血池区域的粗分割结果引导当前片层的磁共振图像的血池区域的粗分割。[0067]具体的，若第M-I个片层或第M+1个片层作为当前片层，第M个片层则是其前一个片层，第M-I个片层则是第M-2个片层的前一个片层，第M+1个片层则是第M+2个片层的前一个片层，依次类推。[0068]具体的，需要先在第M个片层、舒张末期的磁共振图像ImpEd中血池区域完成粗分割后，以第M个片层、舒张末期的磁共振图像ImpEd中血池区域的粗分割结果为基准，对第Μ-ΐ、第M+1个片层、舒张末期的磁共振图像IM-ipED、Im+iped中血池区域进行粗分割，然后以第M-I个片层、舒张末期的磁共振图像IM-1PED的血池区域的粗分割结果为基准，对第M-2个片层、舒张末期的磁共振图像IM-2PED的血池区域进行粗分割;或者，然后以第M+1个片层、舒张末期的磁共振图像Im+1pED的血池区域的粗分割结果为基准，对第M+2个片层、舒张末期的磁共振图像Im+2pED的血池区域进行粗分割，依次类推。[0069]设定第K-1、K+1个片层为当前片层，第K个片层为前一个片层，且第K个片层、舒张末期的磁共振图像IkpEd的血池已经完成粗分割;按照以下步骤获得第Κ_1、Κ+1个片层的舒张末期的磁共振图像IK+1PED、IK-IPED的血池区域的粗分割：[0070]以前一个K片层、舒张末期的磁共振图像IkpED中已粗分割血池区域的质心和长轴加40个像素点作为第K+1、第K-I个片层、舒张末期的磁共振图像IK+1PED、IK-IPED的左心室感兴趣区域的中心和边长，并对左心室感兴趣区域中的图像进行伽马校正，获得校正后的图像IK+1PED’、IK-lPED’；[0071]通过模糊C均值聚类方法对校正后的图像IK+1PED’、IK-1PED’进行分类，选取最亮类作为候选区域SowhSow;[0072]在第Κ_1、Κ+1个片层的候选区域中选取与第K个片层、舒张末期的磁共振图像IkpEd中血池区域的粗分割结果重叠率最大的区域作为第Κ_1、Κ+1个片层、舒张末期的磁共振图像中血池区域的粗分割结果。[0073]进一步的，若当前片层的磁共振图像中粗分割的血池区域的长轴大于前一片层的磁共振图像中粗分割的血池区域的长轴的1.2倍或者当前片层的磁共振图像中粗分割的血池区域的长宽比大于前一片层的磁共振图像中粗分割的血池区域的长宽比的1.3倍，则被判断为出现LV0Tleftventricularoutflowtract现象，其中所述当前片层是位于中间片层位置以下的片层靠近心底层）。[0074]进一步的，如果出现LVOT现象，则调整伽马校正参数增大校正参数），对当前片层的磁共振图像的左心室感兴趣区域中的图像进行伽马校正，并重新进行聚类分割。[0075]进一步的，如采用调整伽马校正参数后并进行聚类分割后的当前片层的磁共振图像中的血池区域仍然被判断为出现LVOT现象，则采用射线扫描方法去除过分割的区域，即可得到分割后的血池区域。[0076]具体的，所述射线扫描方法包括以下步骤：[0077]以当前片层的磁共振图像的血池区域质心为中心点a参图4a作径向射线扫描，提取每个角度射线上距离中心点最近的区域边界点（参图4a中的点bl、b2、b3;[0078]计算上述提取的所有边界点到中心点的径向距离均值mean和标准差std，去掉径向距离大于距离均值mean加标准差std的区域边界点b3，保留的点如图4b所示；[0079]将剩余的区域边界点进行插值后得到闭合曲线d参图4c;闭合曲线包含的区域即为当前片层的磁共振图像的粗分割血池区域。[0080]具体的，上述步骤S4中：二维极坐标转换图像中的X轴表示了扫描射线的序号，y轴表示了射线上点距离扫描起点的距离;扫描的射线数目为180条，扫描起始半径为粗分割血池区域的边界位置减数个像素点（可为5-15个像素点），扫描终止半径为粗分割血池区域的边界位置加数个像素点可为15-45个像素点）。[0081]进一步的，上述步骤S5中：所述双边动态规划中，射线扫描方法顺序获得的射线二维极坐标转换图像的每一列被认为是阶段，射线上的点（二维极坐标转换图像每一列上的点被认为是阶段上的候选点，从第一阶段到最后阶段二维极坐标转换图像的第一列到最后一列）具有最小累积局部能量COSt的路径被认为是最优路径，即检测到的左心室心内、外膜。[0082]具体的，所述局部能量COSt包括内部能量COStint和外部能量COStext，其中内部能量COStint决定了最优路径的平滑性，而外部能量COSte3xt决定了最优路径位于梯度大的位置。[0083]其中，所述内部能量Costint*以下公式获得：[0084][0085]为心内膜候选点在相邻图像列中在垂直方向上的归一化距离；[0086]为心外膜候选点在相邻图像列中在垂直方向上的归一化距离；[0087]为由心内外膜候选点决定的心肌厚度在相邻图像列和cU—D中的归一化距离；[0088]为心内膜候选点在图像第i列中在垂直方向上的位置，为心内膜候选点在图像第i-Ι列中在垂直方向上的位置，为心外膜候选点在图像第i列中在垂直方向上的位置，为心外膜候选点在图像第i-Ι列中在垂直方向上的位置，yrange为心内膜候选点在相邻图像列中在垂直方向上变化距离最大值3个像素），drange为心肌厚度在相邻图像列中的变化距离最大值5个像素）。[0089]其中，所述外部能量为Coste3xt由以下公式获得：[0090][0091]为心内膜的归一化梯度，为心外膜的归一化梯度，心内膜梯度是基于当前层原始的二维极坐标转换图像获得，心外膜梯度是基于当前层增强的二维极坐标转换图像获得。[0092]具体的，所述步骤S6包括：[0093]1选取第K个片层的第t个心动期相当前心动期相的磁共振图像IKt;[0094]2以第K个片层的第t-1个心动期相前一心动期相检测的心内膜确定的血池区域的质心为扫描起点，在第t-Ι个心动期相的磁共振图像中分割的心内、外膜边界附近，将第t个心动期相对应的磁共振图像IKt转换成二维极坐标转换图像Ρια;[0095]3统计从舒张末期至第t-Ι个心动期相的磁共振图像中已分割的心肌灰度的统计分布，通过原始图像IKt中像素点的灰度值与其在心肌灰度统计分布中相应灰度值的分布百分比值相乘，增强原始图像IKt中的心肌部分，并将增强的原始图像IKt转换为二维极坐标转换图像Pcard;[0096]4基于双动态规划方法在二维极坐标转换图像Pkt、Pcard中检测左心室心内、外膜；第t-Ι个心动期相前一心动期相）的分割的左心室心内外膜位置限定第t个心动期相（当前心动期相的左心室心内外膜搜索的范围。[0097]5基于二维极坐标转换图像中检测到的左心室心内、外膜，利用极坐标逆转换获取原始图像IKt中的左心室心内外膜，求凸包并平滑，得到检测结果。[0098]优选的，所述当前心动期相的左心室心内外膜搜索的范围为:前一心动期相上分割的心内膜上20个像素点至外膜边界下10个像素点。[00"]进一步的，在完成左心室心内外I旲的分割后，可进一步获得与左心室相关的心室体积、射血分数和心肌壁厚等特征参数的分析。[0100]本发明专利提出方法步骤中涉及到的参数根据实际医学图像特点可以进行任意设置；[0101]本发明专利中双动态规划技术相关的局部能量定义可以根据图像的灰度、梯度和形状等；[0102]本发明专利中心功能相关特征参数不局限于实施例中所述。[0103]需要说明的是，通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到本申请的部分或全部可借助软件并结合必需的通用硬件平台来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可包括其上存储有机器可执行指令的一个或多个机器可读介质，这些指令在由诸如计算机、计算机网络或其他电子设备等一个或多个机器执行时可使得该一个或多个机器根据本发明的实施例来执行操作。机器可读介质可包括，但不限于，软盘、光盘、CD-ROM紧致盘-只读存储器）、磁光盘、ROM只读存储器）、RAM随机存取存储器）、EPROM可擦除可编程只读存储器）、EEPR0M电可擦除可编程只读存储器）、磁卡或光卡、闪存、或适于存储机器可执行指令的其他类型的介质机器可读介质。[0104]本申请可用于众多通用或专用的计算系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等。[0105]本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。[0106]采用本发明的方法，可实现全自动确定心脏舒张期末期相位及左心室位置，采用聚类方法实现舒张末期图像中血池的分割。[0107]本发明的方法结合采用左心室心内外膜位置和梯度信息，利用双动态规划技术，同时实现心脏舒张期末期中左心室心内外膜的检测；[0108]本发明的方法利用相同片层在不同时序图像中的连续性特点，准确和有效地在四维心功能磁共振图像中实现左心室心内外膜的检测。[0109]本发明的方法可为后续的基于心脏心功能磁共振图像的计算机辅助诊断提供有力的支持和保障，并提高后续诊断的效率和准确性。[0110]虽然本发明已以较佳实施例揭示如上，然其并非用以限定本发明，任何本领域技术人员，在不脱离本发明的精神和范围内，当可作些许的修改和完善，因此本发明的保护范围当以权利要求书所界定的为准。

权利要求：1.一种心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于包括以下步骤：51.获取包含左心室心肌的若干片层、且在不同心动期相的心脏磁共振图像Inp，其中N表示片层的序号，P表示心动期相序号，N、P均为大于或等于1的整数；52.确定舒张末期的期相；53.获得舒张末期的N个片层的磁共振图像的血池区域的粗分割结果；54.以粗分割血池的质心为扫描起点，基于极坐标转换方法，利用射线扫描将舒张末期的每个片层的磁共振图像中左心室感兴趣区域内的图像数据转换为二维极坐标转换图像；55.基于双动态规划方法在二维极坐标转换图像中检测心内膜及心外膜;利用极坐标逆转换获取原始片层图像中的心内、外膜，求凸包并平滑，完成舒张末期的N个片层的磁共振图像中的心内、外膜的分割；56.从舒张末期的磁共振图像中的心内、外膜分割结果向其它心动期相的磁共振图像中的心内、外膜衍生分割结果。2.根据权利要求1所述的心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于，所述二维极坐标转换图像中的X轴表示了扫描射线的序号，y轴表示了射线上点距离扫描起点的距离;扫描的射线数目为180条，扫描起始半径为粗分割血池区域的边界位置减数个像素点，扫描终止半径为粗分割血池区域的边界位置加数个像素点。3.根据权利要求2所述的心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于，所述双边动态规划中，射线扫描方法顺序获得的射线被认为是阶段，射线上的点被认为是阶段上的候选点，从第一阶段到最后阶段具有最小累积局部能量cost的路径被认为是最优路径，即检测到的左心室心内、外膜。4.根据权利要求3所述的心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于，所述局部能量cost包括内部能量COStint和外部能量COSte3xt，其中内部能量COStint决定了最优路径的平滑性，而外部能量COSte3xt决定了最优路径位于梯度大的位置。5.根据权利要求4所述的心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于，所述内部能量COStint由以下公式获得：为心内膜候选点在图像第i列中在垂直方向上的位置，y:为心内膜候选点在图像第i-Ι列中在垂直方向上的位置，yf为心外膜候选点在图像第i列中在垂直方向上的位置，为心外膜候选点在图像第i-Ι列中在垂直方向上的位置，yrange为心内膜候选点在相邻图像列中在垂直方向上变化距离最大值，drange为心肌厚度在相邻图像列中的变化距离最大值。6.根据权利要求4所述的心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于，所述外部能量为cost#由以下公式获得：为心内膜的归一化梯度，为心外膜的归一化梯度，心内膜梯度Gr"是基于当前层原始的二维极坐标转换图像获得，心外膜梯度是基于当前层增强的二维极坐标转换图像获得。7.根据权利要求1所述的心脏心功能磁共振图像中的心内、外膜的分割方法，其特征在于，所述步骤S6包括：选取第K个片层的第t个心动期相的磁共振图像IKt;以第K个片层的第t-1个心动期相检测的心内膜确定的血池区域的质心为扫描起点，在第t-Ι个心动期相的磁共振图像中分割的心内、外膜边界附近，将第t个心动期相对应的磁共振图像IKt转换成二维极坐标转换图像PKt;统计从舒张末期至第t-Ι个心动期相的磁共振图像中已分割的心肌灰度的统计分布，通过原始图像IKt中像素点的灰度值与其在心肌灰度统计分布中相应灰度值的分布百分比值相乘，增强原始图像IKt中的心肌部分，并将增强的原始图像IKt转换为二维极坐标转换图像Pcard;基于双动态规划方法在二维极坐标转换图像Pkt、Pc^d中检测左心室心内、外膜;第t-1个心动期相的分割的左心室心内外膜位置限定第t个心动期相的左心室心内外膜搜索的范围。基于二维极坐标转换图像中检测到的左心室心内、外膜，利用极坐标逆转换获取原始图像Ik冲的左心室心内外膜，求凸包并平滑，得到检测结果。8.根据权利要求7所述的心脏心功能磁共振图像中的心内外膜的分割方法，所述第t个心动期相的左心室心内外膜搜索的范围为：以第t-Ι个心动期相上分割的心内膜上20个像素点至外膜边界下10个像素点。9.根据权利要求1所述的心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于，所述步骤S2包括：在中间片层的每个心动期相所对应的磁共振图像Imp上，设定初始感兴趣区域，计算初始感兴趣区域在所有心动期相上的最大强度投影图像，聚类分割最大强度投影图像，选取圆度最大的区域A;计算中间片层的每个心动期相所对应的磁共振图像Imp在区域A中的图像平均灰度值，具有最大平均灰度值的磁共振图像Imp所在心动期相被定义为舒张末期。10.根据权利要求1所述的心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于，所述步骤S3包括：从中间片层开始向心尖和心底两个方向片层的磁共振图像分别进行分割，完成舒张末期的磁共振图像Inped中的血池区域的粗分割。

百度查询： 上海联影医疗科技股份有限公司 心脏心功能磁共振图像中的心内外膜的分割方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD