超声造影图像的心肌提取及心肌梗塞诊断-控制科学与工程专业论文.docxVIP

Classified Index: TN911.73 U.D.C: 681.39 Dissertation for the Master Degree in Engineering FULL MYOCARDIUM SEGMENTATION OF THE ECHOCARDIOGRAPHIC IMAGES AND MYOCARDIAL INFARCTION DIAGNOSIS Candidate： Jinfu Xing Supervisor： Prof. Yan Wang Academic Degree Applied for： Doctor of Engineering Speciality： Control Science and Engineering Affiliation： School of Astronautics Date of Defence： July, 3rd, 2012 Degree-Conferring-Institution： Harbin Institute of Technology 摘 摘 要 - - I - 摘 要 本文提出了一种心肌超声造影图像心肌内膜与外膜的全自动提取方法及基 于心肌灌注时间强度曲线的希尔伯特黄变换的心肌梗塞诊断方法。心肌梗塞诊 断方法中使用的数据是心肌灌注过程中心肌区域的时间强度变化曲线，而这个 强度变化曲线的每一个数据点就是从每一幅心肌超声造影图像中分割得到的心 肌区域的平均灰度值。 心肌超声造影图像心肌内膜的分割主要使用的是基于数学形态学的方法， 心肌外膜的分割主要使用的是 Chan-Vese 活动轮廓模型的方法。内膜分割使用 数学形态学是因为心腔区域只与心肌区域相连，整个心腔有相似的灰度值且与 心肌区域有明显的差值，简单的形态学方法就取得非常好的效果；之所以外膜 的分割使用不同方法，是因为心肌外膜相邻组织的多样性、复杂性和较高的噪 声，使数学形态学难以取得理想的效果，而基于区域的 Chan-Vese 活动轮廓模 型方法由于对噪声不敏感，所以能够很好的提取初始轮廓，而对完全无法区分 心肌外膜的区域所产生的溢出，根据拟合圆的趋势使用拟合圆进行衔接，取得 了较理想的效果。 心肌梗塞的诊断主要使用的是 HHT 方法，先将心肌区域按顺时针分为 6 个区域，然后使用 EMD 方法分别对每一个区域的灰度值变化曲线进行分解， 得到不同频率的本征模态函数（IMFs）。然后对这些本征模态函数做 Hilbert 变 换，得到各种瞬时特征参数，并最终得到边际谱（marginal spectrum）。对心肌 梗塞区域和正常心肌区域分别进行不同试验，经过分析本文提出了一种根据边 际谱诊断心肌梗塞的方法，如果边际谱的能量主要集中于 1-2Hz 和 6-8Hz 的范 围，那么这个心肌区域是正常的心肌区域，且这两个频率范围分别对应实验用 的 SD 大鼠的呼吸频率和心率；如果边际谱的能量微小或者只是 1-2Hz 处的能 量较高，那么这一区域就是含有心肌梗塞部分的区域。 之所以使用 HHT 方法，是因为考虑到动物体能产生不同频率的信号，如心 跳、呼吸、脑电波等，而 HHT 的核心步骤 EMD 可以将一个信号分解为代表不 同频率的分量。在心跳的过程中，心肌区域必然会受到影响，因此猜想心肌区 域灰度值变化的时间强度曲线也会受到影响，所以使用 EMD 将时间强度曲线 分解，并最终通过反映能量在频率上分布的边际谱总结提出上述诊断方法。实 践证明，这种方法具有快速、准确且能确定病变位置的优点。 关键词：超声心肌造影，数学形态学，活动轮廓模型，希尔伯特黄变换，心肌 梗塞诊断 Abs Abstract - - II - Abstract In this paper we propose a full-automatic method of segmenting the epicardium and endocardium of echocardiography images, we also propose a diagnostic method of myocardial infarction. The data used in myocardial infarction diagnosis is acquired from the time-intensity curve of the whole myocardial perfusion process. And each data point corresponds to one echocardiographic image, because the data is the mean gray value of the myocardium region