基于心脏核磁共振图像的左心室分割及动力学模型研究

摘要

目前，心脏核磁共振影像技术已经成为心脏疾病临床诊断的重要辅助手段。心脏核磁共振影像分析是无创性评价心脏功能的重要方法，运用心脏核磁共振图像不仅能观察到心脏的形态结构，还可以估测心室的整体功能及局部心肌功能，使医师能对心脏的病理和生理状况做出正确的判断。 在临床心脏疾病诊断中，由于左心室是全身血液循环的泵体，在心脏功能中起着重要的作用，因此是目前研究的重点。从心脏核磁共振图像序列提取左心室在收缩期的心肌运动场、位移场、应变应力，进而评定左心室整体和局部心肌功能，在临床具有重要的应用价值。本文着重对基于心脏核磁共振图像的左心室分割及心肌运动进行深入研究与探讨。从左心室的二维分割、三维表面恢复、三维运动重建和应力应变分析等方面着手，提出了相应的二维分割及三维运动分析模型，构建了运用心脏核磁共振图像序列进行左心室三维心肌功能评定的一般框架。本文的工作主要包括以下内容： (1)提出了基于模拟退火算法的简化Snake图像分割模型。该模型对传统Snake模型进行了改进，运用简化snake的思想，使运算简单，并增加系数可变的面积项，使演化曲线不再受初始位置的限制。在求解模型时，引入了模拟退火优化算法，该算法与改进的Snake模型相结合，发挥了模拟退火算法的长处，也降低了运用Snake计算的复杂度。本文算法不仅能较好地处理弱边缘和不规则边缘区域，而且能较好地处理凹陷和拐角区域。 (2)针对标记线对心脏核磁共振图像进行分割时所产生的干扰，本文提出了一种基于细节信号能量的纹理分析算法，该算法利用标记线本身的信息特征，去除了标记线对分割的影响。由于心脏核磁共振图像边缘较弱，噪声相对较强，本文应用测地线活动轮廓模型定义演化曲线，并用水平集方法求解，能够很好地提取左心室的心内膜。在提取左心室心外膜时，为处理边界断裂、缺省的现象，添加了距离约束能量项，解决了依靠梯度和纹理信息无法准确分割的问题。 (3)本文利用带标记线的心脏核磁共振的时序图像，针对左心室建立了具有系数函数的物理可形变模型，对左心室在心脏收缩期的形状和三维运动进行了重建和分析。在对初始时刻左心室表面恢复的过程中，本文提出了反演系数函数成形法；在运动跟踪阶段，本文提出了一种有效的模型节点标记力计算方法，并利用薄板样条对其进行了优化。本文模型不仅能跟踪左心室质点级的运动轨迹，并能最终将左心室在心脏收缩期的运动和形变映射为系数函数的变化曲线，直观易懂，有助于临床应用。 (4)为直观地分析左心室在心脏收缩期的应力分布及形变情况，本文提出了一种左心室力学形态分析方法。首先，利用带标记线的心脏核磁共振图像数据，针对左心室建立系数可变的物理可形变模型；在所建模型的基础上，利用心脏收缩期各个时刻的左心室轮廓点数据进行三维重建。其次，利用心脏收缩期各相邻时刻的标记点数据计算左心室模型标记力。最后，将模型标记力转换为收缩应力分量、切向应力分量和拉伸应力分量，并将各个应力分量用彩色云图显示。本文提出的左心室力学形态分析方法，能直观有效地反映左心室内外表面在整个心脏收缩期的应力分布及形变趋势。 (5)利用时间序列的带标记线心脏核磁共振图像，提出了一种基于有限元方法的左心室运动分析模型，重建了左心室在心脏收缩期的三维位移场及应变分布。本文首先提出了一种基于Delamlay三角网的三维表面恢复算法，该算法能对表达形式多变、空间分布不一致的稀疏点进行表面恢复。然后建立左心室有限元模型，重建三维运动，获取左心室形变信息。本文提出的左心室表面恢复算法及基于有限元方法的左心室模型，能有效地获取左心室的形状信息和运动信息，并为临床诊断提供有用的参考数据。 (6)提出了基于有限元方法的左心室生物力学模型。心肌是生物物质高度优化的复合材料，心肌纤维的方向是连续变化的，本文在基于有限元方法的左心室模型的基础上，进一步考虑到心肌纤维方向，将心肌的物质材料属性添加到模型中，建立了左心室生物力学模型。在该模型基础上，本文分析了左心室在局部纤维坐标系下的应力应变分量，给出了各个分量的彩色云图，这对临床应用具有重要的意义。

目录

文摘

英文文摘

声明

1.绪论

1.1引言

1.2心脏MR成像技术

1.3心脏MR图像分析综述

1.3.1心脏MR图像分割

1.3.2基于MR的左心室运动重建

1.4课题来源

1.5本文的主要研究工作和内容安排

1.6本文工作的创新点

2.心脏MR图像数据准备

2.1图像采集

2.2左心室轮廓提取

2.3标记线跟踪

2.4坐标系间的转换

2.4.1图像坐标系到世界坐标系

2.4.2世界坐标系到局部坐标系

2.5心脏MR图像数据集

2.5.1轮廓点数据

2.5.2标记点数据

3.心脏MR图像分割

3.1引言

3.2参数活动轮廓模型

3.2.1经典Snake模型

3.2.2数值解法

3.2.3 Snake模型的改进模型

3.2.4 Snake模型的优化问题

3.3基于模拟退火的简化Snake模型

3.3.1分割模型

3.3.2模拟退火思想

3.3.3算法流程

3.3.4实验结果及分析

3.4本章小结

4.带标记线的心脏MR图像分割

4.1引言

4.2纹理分析

4.2.1清晰度

4.2.2细节信号能量

4.2.3边缘信号能量

4.2.4特征能量图的计算

4.3纹理表示的验证

4.3.1实验设计

4.3.2 MTF与清晰度

4.3.3 MTF与细节信号能量

4.3.4 MTF与边缘信号能量

4.3.5 MTF与图像参数小结

4.4基于细节信号能量图的带标记线心脏MR图像分割

4.4.1细节信号能量图的计算

4.4.2左心室心内膜分割

4.4.3左心室心外膜分割

4.4.4实验结果与分析

4.5推广应用——基于特征能量图的纹理图像分割

4.5.1 C-V模型概述

4.5.2实验结果与分析

4.6本章小结

5.基于可形变模型的左心室运动分析

5.1引言

5.2基于物理的可形变模型

5.2.1可形变模型的几何定义

5.2.2可形变模型的动力学定义

5.3基于可形变模型的左心室运动分析

5.3.1左心室表面恢复

5.3.2左心室标记力计算及优化

5.3.3实验结果及模型验证

5.4基于可形变模型的左心室力学形态分析

5.5本章小结

6.基于有限元方法的左心室运动分析

6.1引言

6.2有限元方法

6.2.1有限元方法的特性

6.2.2有限元方法的基本步骤

6.2.3三维单元

6.3左心室表面恢复

6.3.1基于Delaunay三角网的左心室表面恢复

6.3.2左心室表面恢复算法的验证

6.4基于有限元方法的左心室应变分析

6.4.1基于有限元方法的左心室模型

6.4.2实验结果及分析

6.5基于有限元方法的左心室生物力学模型

6.5.1左心室生物力学模型

6.5.2实验结果及分析

6.6本章小结

7.结束语

7.1本文工作的总结

7.2将来的工作

参考文献

致 谢

附录