CN106910194B - 心脏心功能磁共振图像中的心内外膜的分割方法 - Google Patents

Abstract

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

Description

心脏心功能磁共振图像中的心内外膜的分割方法

【技术领域】

本发明涉及医学图像领域，尤其是涉及心脏心功能磁共振图像的分割。

【背景技术】

近年来，心血管疾病的发病率和死亡率正逐年增加。每年在美国有近600万人被检查出患有心血管疾病，而我国则有近300万人死于心血管疾病。心脏功能是诊断心血管疾病的重要指标，可以通过分析四维心功能磁共振图像中左心室相关的全局和局部特征参数实现，如：心室体积、射血分数和心肌壁厚等。心功能特征参数提取的准确性依赖于左心室心内外膜的准确定位，而计算机辅助诊断方案有助于提高在四维心功能磁共振图像中定位左心室心内外膜的可靠性和效率。

由于四维心功能磁共振图像在不同病患扫描图像、相同病患的不同扫描断层图像以及相同扫描断层图像中灰度和梯度分布的巨大差异，另外，左心室心外膜周边临近组织的多样性，这些都增加了左心室心内外膜检测的难度。虽然在过去的近20年中，已有很多研究方法和结果的报道，但是在四维心功能磁共振影像中对左心室心内外膜进行准确的检测仍然没有解决。针对四维心功能磁共振图像中左心室心内外膜的检测，手动和需用户辅助的半自动检测方法比较繁复、效率低、并且存在较大的观察者之间和观察者自身的差异。在全自动检测方面，由于四维心功能磁共振图像在不同病患扫描图像、相同病患的不同 扫描断层图像以及相同扫描断层图像中灰度和梯度分布的巨大差异，另外，左心室心外膜周边临近组织的多样性，导致简单的基于阈值的检测方法或者基于梯度的检测方法不足以准确地提取出左心室心内外膜。基于模型的检测方法需要大量的训练样本去训练模型以获得较为准确的检测结果，而这些精准训练样本的获取则是一项更艰巨和繁复的任务。

【发明内容】

本发明所要解决的技术问题是提供一种心脏心功能磁共振图像中的心内外膜的分割方法。

本发明为解决上述技术问题而采用的技术方案是：一种心脏心功能磁共振图像中的心内外膜的分割方法,其包括以下步骤:

S1.获取包含左心室心肌的若干片层、且在不同心动期相的心脏磁共振图像I_NP，其中N表示片层的序号，P表示心动期相序号，N、P均为大于或等于1的整数；

S2.确定舒张末期的期相；

S3.获得舒张末期的N个片层的磁共振图像的血池区域的粗分割结果；

S4.以粗分割血池的质心为扫描起点，基于极坐标转换方法，利用射线扫描将舒张末期的每个片层的磁共振图像中左心室感兴趣区域内的图像数据转换为二维极坐标转换图像；

S5.基于双动态规划方法在二维极坐标转换图像中检测心内膜及心外膜；利用极坐标逆转换获取原始片层图像中的心内、外膜，求凸包并平滑，完成舒张末期的N个片层的磁共振图像中的心内、外膜的分割；

S6.从舒张末期的磁共振图像中的心内、外膜分割结果向其它心动期相的磁共振图像中的心内、外膜衍生分割结果。

优选的，所述二维极坐标转换图像中的x轴表示了扫描射线的序号，y轴表示了射线上点距离扫描起点的距离；扫描的射线数目为180条，扫描起始半径为粗分割血池区域的边界位置减数个像素点，扫描终止半径为粗分割血池区域的边界位置加数个像素点。

优选的，所述双边动态规划中，射线扫描方法顺序获得的射线(二维极坐标转换图像的每一列)被认为是阶段，射线上的点(二维极坐标转换图像每一列上的点)被认为是阶段上的候选点，从第一阶段到最后阶段(二维极坐标转换图像的第一列到最后一列)具有最小累积局部能量cost的路径被认为是最优路径，即检测到的左心室心内、外膜。

优选的，所述局部能量cost包括内部能量cost_int和外部能量cost_ext，其中内部能量cost_int决定了最优路径的平滑性，而外部能量cost_ext决定了最优路径位于梯度大的位置。

优选的，所述内部能量cost_int由以下公式获得:

为心内膜候选点在相邻图像列中在垂直方向上的归一化距离；

为心外膜候选点在相邻图像列中在垂直方向上的归一化距离；

为由心内外膜候选点决定的心肌厚度在相邻图像列(

和d_i-1)中的归一化距离；

为心内膜候选点在图像第i列中在垂直方向上的位置，

为心内膜候选点在图像第i-1列中在垂直方向上的位置，

为心外膜候选点在图像第i列中在垂直方向上的位置，

为心外膜候选点在图像第i-1列中在垂直方向上的位置，y_range为心内膜候选点在相邻图像列中在垂直方向上变化距离最大值(3个像素)，d_range为心肌厚度在相邻图像列中的变化距离最大值(5个像素)。

优选的，所述外部能量为cost_ext由以下公式获得:

为心内膜的归一化梯度，

为心外膜的归一化梯度，心内膜梯度

是基于当前层原始的二维极坐标转换图像获得，心外膜梯度

是基于当前层增强的二维极坐标转换图像获得。

优选的，所述步骤S6包括：

选取第K个片层的第t个心动期相(当前心动期相)的磁共振图像I_Kt；

以第K个片层的第t-1个心动期相(前一心动期相)检测的心内膜确定的血池区域的质心为扫描起点，在第t-1个心动期相的磁共振图像中分割的心内、外膜边界附近，将第t个心动期相对应的磁共振图像I_Kt转换成二维极坐标转换图像P_Kt；

统计从舒张末期至第t-1个心动期相的磁共振图像中已分割的心肌灰度的统计分布，通过原始图像I_Kt中像素点的灰度值与其在心肌灰度统计分布中相应灰度值的分布百分比值相乘，增强原始图像I_Kt中的心肌部分，并将增强的原始图像I_Kt转换为二维极坐标转换图像P_card；

基于双动态规划方法在二维极坐标转换图像P_kt、P_card中检测左心室心内、外膜；第t-1个心动期相(前一心动期相)的分割的左心室心内外膜位置限定第t个心动期相(当前心动期相)的左心室心内外膜搜索的范围。

基于二维极坐标转换图像中检测到的左心室心内、外膜，利用极坐标逆转换获取原始图像I_Kt中的左心室心内外膜，求凸包并平滑，得到检测结果。

优选的，所述第t个心动期相(当前心动期相)的左心室心内外膜搜索的范围为：以第t-1个心动期相(前一心动期相)上分割的心内膜上20个像素点至外膜边界下10个像素点。

优选的，所述步骤S2包括：

在中间片层的每个心动期相所对应的磁共振图像I_MP上，设定初始感兴趣区域，计算初始感兴趣区域在所有心动期相上的最大强度投影图像，聚类分割最大强度投影图像，选取圆度最大的区域A；

计算中间片层的每个心动期相所对应的磁共振图像I_MP在区域A中的图像平均灰度值，具有最大平均灰度值的磁共振图像I_MP所在心动期相被定义为舒张末期。

优选的，所述步骤S3包括：从中间片层开始向心尖和心底两个方向片层的磁共振图像分别进行分割，完成舒张末期的磁共振图像I_NP(ED)中的血池区域的粗分割。

本发明对比现有技术有如下的有益效果：本发明的方法结合采用左心室心内外膜位置和梯度信息，利用双动态规划技术，同时实现心脏舒张期末期中左心室心内外膜的检测；本发明的方法利用相同片层在不同时序图像中的连续性特点，准确和有效地在四维心功能磁共振图像中实现左心室心内外膜的检测。本发明的方法可为后续的基于心脏心功能磁共振图像的计算机辅助诊断提供有力的支持和保障，并提高后续诊断的效率和准确性。

【附图说明】

图1为本明实施例中的一种心脏心功能磁共振图像中的心内外膜的分割方法；

图2为本发明实施例的一种心脏心功能磁共振图像中舒张末期图像中血池的分割方法的流程图；

图3a为中间片层磁共振图像在感兴趣区域中的最大强度投影图像；

图3b为中间片层磁共振图像的在初始感兴趣区域中的图像；

图3c为通过模糊C均值聚类方法对图3b进行聚类分割的结果；

图3d为中间片层磁共振图像中血池区域分割的结果；

图4a-4c通过射线扫描方法对过分割区域处理的过程示意图；

图5为本明实施例中的一种心脏心功能磁共振图像中的心内外膜的分割结果。

【具体实施方式】

请参阅图1-5，本发明实施例中的一种心脏心功能磁共振图像中的心内外膜的分割方法,其特征在于包括以下步骤:

S1.获取包含左心室心肌的若干片层、且在不同心动期相的心脏磁共振图像I_NP，其中N表示片层的序号，P表示心动期相序号，N、P均为大于或等于1的整数；

S2.确定舒张末期的期相；

S3.获得舒张末期的N个片层的磁共振图像的血池区域的粗分割结果；

S4.以粗分割血池的质心为扫描起点，基于极坐标转换方法，利用射线扫描将舒张末期的每个片层的磁共振图像中左心室感兴趣区域内的图像数据转换为二维极坐标转换图像；

S5.基于双动态规划方法在二维极坐标转换图像中检测心内膜及心外膜；利用极坐标逆转换获取原始片层图像中的心内、外膜，求凸包并平滑，完成舒张末期的N个片层的磁共振图像中的心内、外膜的分割；

S6.从舒张末期的磁共振图像中的心内、外膜分割结果向其它心动期相(P-1个心动期相)的磁共振图像中的心内、外膜衍生分割结果。

具体的，上述S1中的每一片层、且在每个期相中的心脏磁共振图像是通过以下方式得到：

通过采集多个心动周期的、对应相同心动期相中的磁共振数据线，将所采集到的磁共振数据线填充到一个K-Space(K空间)中，通过傅里叶变换得到相应片层在相应心动期相的磁共振图像。通过设置较小的时间窗，例如20-50ms，并在多个心动周期采集数据，可减小心脏运动伪影对图像的影响，提高图像的质量。

心动期相是指将自左心室的收缩期开始至舒张末期的完整心动周期T的划分为多个时间段(phase),这些时间段可以为一系类连续的时间序列，也可以为间隔的时间序列。

在心室包含心肌部分的区域沿心脏的长轴方向获取若干个片层(slice)，并通过磁共振扫描获得每个片层、在每个心动期相中的心脏磁共振图像I_NP，其中N表示相应片层的序号，P表示心动期相(phase)序号，N、P均为大于或等于1的整数；具体的，N的取值可以为3、5、6等，P的取值可以为3、5、8、10等；另外，定义第M个片层为中间片层，其中，M为大于或等于2的整数。

具体的，上述S2中确定舒张末期的期相，通过以下方法：

在中间片层的每个心动期相所对应的磁共振图像I_MP上，设定初始感兴趣区域，计算初始感兴趣区域在所有心动期相(P个心动期相)上的最大强度投影图像(MIP)，聚类分割最大强度投影图像(MIP)，选取圆度最大的区域A；优选的聚类分割方法为模糊C均值聚类方法；

计算中间片层的每个心动期相所对应的磁共振图像I_MP在区域A中的图像平均灰度值，具有最大平均灰度值的磁共振图像I_MP所在心动期相被定义为舒张 末期P(ED)，并确定中间片层在舒张末期对应的磁共振图像I_MP(ED)；

从中间片层开始向心尖和心底两个方向片层的磁共振图像分别进行分割，完成(N个片层的)舒张末期的磁共振图像I_NP(ED)中的血池的粗分割。

初始感兴趣区域ROI₀的中心为原图像(即磁共振图像I_MP)的中心，半径为若干个像素点，半径的大小为60-120个像素点，优选半径大小为100个像素点。

然后，计算初始感兴趣区域在所有心动期相(P个心动期相)上的最大强度投影图像(MIP)，通过模糊C均值聚类方法分割最大强度投影图像(MIP)，选取圆度最大的区域A；其中，聚类(分类)的数目可以为2、3或4等；圆度的定义为(周长*周长)/(4*PI*面积)。

进一步的，获取左心室感兴趣区域，所述左心室感兴趣区域的中心为区域A的质心，所述左心室感兴趣区域的边长为区域A的长轴长度加数个像素点，具体可为10-30个像素点，优选为20个像素点。

具体的，对中间片层的舒张末期的磁共振图像I_MP(ED)中的血池的粗分割包括以下步骤：

对中间片层的舒张末期的磁共振图像I_MP(ED)在左心室感兴趣区域中的图像进行伽马校正(Gamma Correction)，获得校正后的图像I_MP(ED)’；通过伽马校正，可增强血池与心肌的对比度，方便后续的识别和分割；

通过模糊C均值聚类方法对校正后的图像I_MP(ED)’进行分类，分类数目为2，选取最亮类作为候选区域S_M；

在候选区域S_M中选取与区域A重叠率最大的区域作为中间片层的舒张末期的磁共振图像I_MP(ED)中的血池区域(Blood Pool，BP)的粗分割结果。

进一步的，对第一个至第M-1个片层的舒张末期的磁共振图像的血池粗分 割，或第M+1个至第N个片层的舒张末期的磁共振图像的血池粗分割通过以下方式获得：以前一个片层的磁共振图像的血池区域的粗分割结果引导当前片层的磁共振图像的血池区域的粗分割。

具体的，若第M-1个片层或第M+1个片层作为当前片层，第M个片层则是其前一个片层，第M-1个片层则是第M-2个片层的前一个片层，第M+1个片层则是第M+2个片层的前一个片层，依次类推。

具体的，需要先在第M个片层、舒张末期的磁共振图像I_MP(ED)中血池区域完成粗分割后，以第M个片层、舒张末期的磁共振图像I_MP(ED)中血池区域的粗分割结果为基准，对第M-1、第M+1个片层、舒张末期的磁共振图像I_(M-1)P(ED)、I _(M+1)P(ED)中血池区域进行粗分割，然后以第M-1个片层、舒张末期的磁共振图像I_(M-1)P(ED)的血池区域的粗分割结果为基准，对第M-2个片层、舒张末期的磁共振图像I_(M-2)P(ED)的血池区域进行粗分割；或者，然后以第M+1个片层、舒张末期的磁共振图像I_(M+1)P(ED)的血池区域的粗分割结果为基准，对第M+2个片层、舒张末期的磁共振图像I_(M+2)P(ED)的血池区域进行粗分割，依次类推。

设定第K-1、K+1个片层为当前片层，第K个片层为前一个片层，且第K个片层、舒张末期的磁共振图像I_KP(ED)的血池已经完成粗分割；按照以下步骤获得第K-1、K+1个片层的舒张末期的磁共振图像I_{(K+1)P(ED)、}I_(K-1)P(ED)的血池区域的粗分割：

以前一个(K)片层、舒张末期的磁共振图像I_KP(ED)中已粗分割血池区域的质心和长轴加40个像素点作为第K+1、第K-1个片层、舒张末期的磁共振图像I_{(K+1)P(ED)、}I_(K-1)P(ED)的左心室感兴趣区域的中心和边长，并对左心室感兴趣区域中的图像进行伽马校正，获得校正后的图像I_(K+1)P(ED)’、I_(K-1)P(ED)’；

通过模糊C均值聚类方法对校正后的图像I_(K+1)P(ED)’、I_(K-1)P(ED)’进行分 类，选取最亮类作为候选区域S_(K-1)、S_(K+1)；

在第K-1、K+1个片层的候选区域S_(K-1)、S_(K+1)中选取与第K个片层、舒张末期的磁共振图像I_KP(ED)中血池区域的粗分割结果重叠率最大的区域作为第K-1、K+1个片层、舒张末期的磁共振图像中血池区域的粗分割结果。

进一步的，若当前片层的磁共振图像中粗分割的血池区域的长轴大于前一片层的磁共振图像中粗分割的血池区域的长轴的1.2倍或者当前片层的磁共振图像中粗分割的血池区域的长宽比大于前一片层的磁共振图像中粗分割的血池区域的长宽比的1.3倍，则被判断为出现LVOT(left ventr icular outflow tract)现象，其中所述当前片层是位于中间片层位置以下的片层(靠近心底层)。

进一步的，如果出现LVOT现象，则调整伽马校正参数(增大校正参数)，对当前片层的磁共振图像的左心室感兴趣区域中的图像进行伽马校正，并重新进行聚类分割。

进一步的，如采用调整伽马校正参数后并进行聚类分割后的当前片层的磁共振图像中的血池区域仍然被判断为出现LVOT现象，则采用射线扫描方法去除过分割的区域，即可得到分割后的血池区域。

具体的，所述射线扫描方法包括以下步骤：

以当前片层的磁共振图像的血池区域质心为中心点a(参图4a)作径向射线扫描,提取每个角度射线上距离中心点最近的区域边界点(参图4a中的点b1、b2、b3)；

计算上述提取的所有边界点到中心点的径向距离均值mean和标准差std，去掉径向距离大于距离均值(mean)加标准差(std)的区域边界点(b3)，保留的点如图4b所示；

将剩余的区域边界点进行插值后得到闭合曲线d(参图4c)；闭合曲线包含的区域即为当前片层的磁共振图像的粗分割血池区域。

具体的，上述步骤S4中：二维极坐标转换图像中的x轴表示了扫描射线的序号，y轴表示了射线上点距离扫描起点的距离；扫描的射线数目为180条，扫描起始半径为粗分割血池区域的边界位置减数个像素点(可为5-15个像素点)，扫描终止半径为粗分割血池区域的边界位置加数个像素点(可为15-45个像素点)。

进一步的，上述步骤S5中：所述双边动态规划中，射线扫描方法顺序获得的射线(二维极坐标转换图像的每一列)被认为是阶段，射线上的点(二维极坐标转换图像每一列上的点)被认为是阶段上的候选点，从第一阶段到最后阶段(二维极坐标转换图像的第一列到最后一列)具有最小累积局部能量cost的路径被认为是最优路径，即检测到的左心室心内、外膜。

具体的，所述局部能量cost包括内部能量cost_int和外部能量cost_ext，其中内部能量cost_int决定了最优路径的平滑性，而外部能量cost_ext决定了最优路径位于梯度大的位置。

其中，所述内部能量cost_int由以下公式获得:

为心内膜候选点在相邻图像列中在垂直方向上的归一化距离；

为心外膜候选点在相邻图像列中在垂直方向上的归一化距离；

为由心内外膜候选点决定的心肌厚度在相邻图像列(

和d_i-1)中的归一化距离；

为心内膜候选点在图像第i列中在垂直方向上的位置，

为心内膜候选点在图像第i-1列中在垂直方向上的位置，

为心外膜候选点在图像第i 列中在垂直方向上的位置，

为心外膜候选点在图像第i-1列中在垂直方向上的位置，y_range为心内膜候选点在相邻图像列中在垂直方向上变化距离最大值(3个像素)，d_range为心肌厚度在相邻图像列中的变化距离最大值(5个像素)。

其中，所述外部能量为cost_ext由以下公式获得:

为心内膜的归一化梯度，

为心外膜的归一化梯度，心内膜梯度

是基于当前层原始的二维极坐标转换图像获得，心外膜梯度

是基于当前层增强的二维极坐标转换图像获得。

具体的，所述步骤S6包括：

1)选取第K个片层的第t个心动期相(当前心动期相)的磁共振图像I_Kt；

2)以第K个片层的第t-1个心动期相(前一心动期相)检测的心内膜确定的血池区域的质心为扫描起点，在第t-1个心动期相的磁共振图像中分割的心内、外膜边界附近，将第t个心动期相对应的磁共振图像I_Kt转换成二维极坐标转换图像P_Kt；

3)统计从舒张末期至第t-1个心动期相的磁共振图像中已分割的心肌灰度的统计分布，通过原始图像I_Kt中像素点的灰度值与其在心肌灰度统计分布中相应灰度值的分布百分比值相乘，增强原始图像I_Kt中的心肌部分，并将增强的原始图像I_Kt转换为二维极坐标转换图像P_card；

4)基于双动态规划方法在二维极坐标转换图像P_kt、P_card中检测左心室心内、外膜；第t-1个心动期相(前一心动期相)的分割的左心室心内外膜位置限定第t个心动期相(当前心动期相)的左心室心内外膜搜索的范围。

5)基于二维极坐标转换图像中检测到的左心室心内、外膜，利用极坐标 逆转换获取原始图像I_Kt中的左心室心内外膜，求凸包并平滑，得到检测结果。

优选的，所述当前心动期相的左心室心内外膜搜索的范围为：前一心动期相上分割的心内膜上20个像素点至外膜边界下10个像素点。

进一步的，在完成左心室心内外膜的分割后，可进一步获得与左心室相关的心室体积、射血分数和心肌壁厚等特征参数的分析。

本发明专利提出方法步骤中涉及到的参数根据实际医学图像特点可以进行任意设置；

本发明专利中双动态规划技术相关的局部能量定义可以根据图像的灰度、梯度和形状等；

本发明专利中心功能相关特征参数不局限于实施例中所述。

需要说明的是，通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到本申请的部分或全部可借助软件并结合必需的通用硬件平台来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可包括其上存储有机器可执行指令的一个或多个机器可读介质，这些指令在由诸如计算机、计算机网络或其他电子设备等一个或多个机器执行时可使得该一个或多个机器根据本发明的实施例来执行操作。机器可读介质可包括，但不限于，软盘、光盘、CD-ROM(紧致盘-只读存储器)、磁光盘、ROM(只读存储器)、RAM(随机存取存储器)、EPROM(可擦除可编程只读存储器)、EEPROM(电可擦除可编程只读存储器)、磁卡或光卡、闪存、或适于存储机器可执行指令的其他类型的介质/机器可读介质。

本申请可用于众多通用或专用的计算系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基 于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等。

本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。

采用本发明的方法，可实现全自动确定心脏舒张期末期相位及左心室位置，采用聚类方法实现舒张末期图像中血池的分割。

本发明的方法结合采用左心室心内外膜位置和梯度信息，利用双动态规划技术，同时实现心脏舒张期末期中左心室心内外膜的检测；

本发明的方法利用相同片层在不同时序图像中的连续性特点，准确和有效地在四维心功能磁共振图像中实现左心室心内外膜的检测。

本发明的方法可为后续的基于心脏心功能磁共振图像的计算机辅助诊断提供有力的支持和保障，并提高后续诊断的效率和准确性。

虽然本发明已以较佳实施例揭示如上，然其并非用以限定本发明，任何本领域技术人员，在不脱离本发明的精神和范围内，当可作些许的修改和完善，因此本发明的保护范围当以权利要求书所界定的为准。

Claims (9)

1.一种心脏心功能磁共振图像中的心内外膜的分割方法,其特征在于包括以下步骤:

S1.获取包含左心室心肌的若干片层、且在不同心动期相的心脏磁共振图像I_NP，其中N表示片层的序号，P表示心动期相序号，N、P均为大于或等于1的整数；

S2.确定舒张末期的期相；

S3.获得舒张末期的N个片层的磁共振图像的血池区域的粗分割结果；

S4.以粗分割血池的质心为扫描起点，基于极坐标转换方法，利用射线扫描将舒张末期的每个片层的磁共振图像中左心室感兴趣区域内的图像数据转换为二维极坐标转换图像；

S5.基于双动态规划方法在二维极坐标转换图像中检测心内膜及心外膜；利用极坐标逆转换获取原始片层图像中的心内、外膜，求凸包并平滑，完成舒张末期的N个片层的磁共振图像中的心内、外膜的分割；

S6.从舒张末期的磁共振图像中的心内、外膜分割结果向其它心动期相的磁共振图像中的心内、外膜衍生分割结果；所述步骤S6包括：

选取第K个片层的第t个心动期相的磁共振图像I_Kt；

以第K个片层的第t-1个心动期相检测的心内膜确定的血池区域的质心为扫描起点，在第t-1个心动期相的磁共振图像中分割的心内、外膜边界附近，将第t个心动期相对应的磁共振图像I_Kt转换成二维极坐标转换图像P_Kt；

统计从舒张末期至第t-1个心动期相的磁共振图像中已分割的心肌灰度的统计分布，通过原始图像I_Kt中像素点的灰度值与其在心肌灰度统计分布中相应灰度值的分布百分比值相乘，增强原始图像I_Kt中的心肌部分，并将增强的原始图像I_Kt转换为二维极坐标转换图像P_card；

基于双动态规划方法在二维极坐标转换图像P_kt、P_card中检测左心室心内、外膜；第t-1个心动期相的分割的左心室心内外膜位置限定第t个心动期相的左心室心内外膜搜索的范围；

基于二维极坐标转换图像中检测到的左心室心内、外膜，利用极坐标逆转换获取原始图像I_Kt中的左心室心内外膜，求凸包并平滑，得到检测结果。

2.根据权利要求1所述的心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于，所述二维极坐标转换图像中的x轴表示了扫描射线的序号，y轴表示了射线上点距离扫描起点的距离；扫描的射线数目为180条，扫描起始半径为粗分割血池区域的边界位置减数个像素点，扫描终止半径为粗分割血池区域的边界位置加数个像素点。

3.根据权利要求2所述的心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于，所述双动态规划中，射线扫描方法顺序获得的射线被认为是阶段，射线上的点被认为是阶段上的候选点，从第一阶段到最后阶段具有最小累积局部能量cost的路径被认为是最优路径，即检测到的左心室心内、外膜。

4.根据权利要求3所述的心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于，所述局部能量cost包括内部能量cost_int和外部能量cost_ext，其中内部能量cost_int决定了最优路径的平滑性，而外部能量cost_ext决定了最优路径位于梯度大的位置。

5.根据权利要求4所述的心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于，所述内部能量cost_int由以下公式获得:

为心内膜候选点在相邻图像列中在垂直方向上的归一化距离；

为心外膜候选点在相邻图像列中在垂直方向上的归一化距离；

为由心内外膜候选点决定的心肌厚度在相邻图像列

和

中的归一化距离；

为心内膜候选点在图像第i列中在垂直方向上的位置，

为心内膜候选点在图像第i-1列中在垂直方向上的位置，

为心外膜候选点在图像第i列中在垂直方向上的位置，

为心外膜候选点在图像第i-1列中在垂直方向上的位置，y_range为心内膜候选点在相邻图像列中在垂直方向上变化距离最大值，d_range为心肌厚度在相邻图像列中的变化距离最大值。

6.根据权利要求4所述的心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于，所述外部能量为cost_ext由以下公式获得:

为心内膜的归一化梯度，

为心外膜的归一化梯度，心内膜梯度

是基于当前层原始的二维极坐标转换图像获得，心外膜梯度

是基于当前层增强的二维极坐标转换图像获得。

7.根据权利要求1所述的心脏心功能磁共振图像中的心内外膜的分割方法，所述第t个心动期相的左心室心内外膜搜索的范围为：以第t-1个心动期相上分割的心内膜上20个像素点至外膜边界下10个像素点。

8.根据权利要求1所述的心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于，所述步骤S2包括：

在中间片层的每个心动期相所对应的磁共振图像I_MP上，设定初始感兴趣区域，计算初始感兴趣区域在所有心动期相上的最大强度投影图像，聚类分割最大强度投影图像，选取圆度最大的区域A；

计算中间片层的每个心动期相所对应的磁共振图像I_MP在区域A中的图像平均灰度值，具有最大平均灰度值的磁共振图像I_MP所在心动期相被定义为舒张末期。

9.根据权利要求1所述的心脏心功能磁共振图像中的心内外膜的分割方法，其特征在于，所述步骤S3包括：

从中间片层开始向心尖和心底两个方向片层的磁共振图像分别进行分割，完成舒张末期的磁共振图像I_NP(ED)中的血池区域的粗分割。

Images