基于栈式自编码器和支持向量机的无创心脏梗死分类模型构建方法技术

本发明专利技术公开了一种基于栈式自编码器和支持向量机的无创心脏梗死分类模型构建方法，在训练阶段先由采集到的心脏磁共振图像求得心脏的运动信息，将图像信息与运动信息融合的图像块作为输入，对应的梗死状况作为标签训练自编码器，利用降噪自编码器对输入数据做预处理，在学习高维信息过程中利用了上一层所有的变量因素，最终学得输入数据的深层特征；之后以学得的深层特征作输入，加上对应标签，利用支持向量机对样本分类。本发明专利技术从数据驱动的角度实现心肌梗死分类预测，有效地解决了临床上通过注射造影剂预测梗死费时费力问题。

【技术实现步骤摘要】
基于栈式自编码器和支持向量机的无创心脏梗死分类模型构建方法
本专利技术属于医学图像分析
，具体涉及一种基于栈式自编码器和支持向量机的无创心脏梗死分类模型构建方法。
技术介绍
心肌梗死(myocardialinfarction，MI)是临床常见的心血管系统疾病，是冠心病的危重类型，死亡率高。心肌坏死后将进一步引起心脏重构，导致心律失常和心力衰竭，心脏重构包括心室重构、血管重构、神经重构及电重构，临床表现为心肌质量、心腔容量及形态结构的改变，心脏负荷增加，顺应性降低，内膜下心肌灌注减少，心肌耗氧量增加，代偿机制启动障碍，心肌电-机械失同步；形态学上表现为心肌细胞肥大、心肌细胞凋亡、间质纤维母细胞增生及间质纤维化。上述一系列表现造成心功能的进行性下降，增加心律失常的发生，最终导致心力衰竭，因此有效地控制心肌梗死后心脏重构对改善患者预后有重要的意义。在此之前，已经提出了许多检测心肌梗塞的方法，如JWvanOorschot等人已经证明了在没有使用外源造影剂的情况下检测慢性心肌梗塞患者的可行性；然而，该方法的灵敏度和特异性较低，并且不能提供关于心肌组织特征的足够信息；Asha等人使用三个分类器决策树，朴素贝叶斯和神经网络开发了一个智能心脏疾病预测系统来预测心脏病；AnbarasiM等提出了使用遗传算法对特征子集选择的心脏疾病的增强预测，从原来的预测心脏疾病的十三个属性降到六个属性；然而，他们都没有直接和准确地预测MI的位置和面积。CarlozOrdonez等人提出了四个约束来减少规则的数量：项目过滤，属性分组，最大项目集大小和先行/后续规则过滤；但是，将关联规则应用于医疗数据集时，将会生成大量的规则，这些规则大部分在医学上是微不足道的，找到他们所需要的时间是不切实际的。MI患者的症状因人而异，患者可能会出现症状多样性；另外，没有任何症状的人也可能是心肌梗塞的患者。因此，MI的早期预测将在诊断中起关键作用，通过早期识别症状，可以得到最有效的临床治疗。近年来，常常使用计算机辅助检测(CAD)系统来自动分析心脏左心室(LV)心肌功能，这些系统将在病人没有痛苦的情况下帮助心脏病专家检测心肌梗死。当然，CAD系统并不是要取代心脏病专家，而是为医生提供建议。
技术实现思路
鉴于上述，本专利技术提供了一种基于栈式自编码器和支持向量机的无创心脏梗死分类模型构建方法，在只得到病人的心脏磁共振电影序列的情况下，利用该方法构建得到的模型能够很好地检测出心肌梗死的位置及区域大小。一种基于栈式自编码器和支持向量机的无创心脏梗死分类模型构建方法，包括如下步骤：(1)获取受试者的完整心脏磁共振图像序列及其对应的延迟增强图像；(2)通过对完整心脏磁共振图像序列进行筛选，只选取其中心脏短轴图像，从而得到对应的心脏短轴电影序列，进而利用opticalflow计算出心脏短轴电影序列中与延迟增强图像所对应的一张心脏短轴图像M的运动信息；(3)从心脏短轴图像M中提取心肌区域对应的图像信息和运动信息，遍历心肌区域中所有n×n大小的图像块并利用延迟增强图像标记每个图像块的类别，即心肌梗死与否，n为大于1的自然数；(4)对图像块的图像信息和运动信息进行融合，将得到的融合信息以及图像块的类别信息组合成一个样本；(5)根据步骤(1)～(4)遍历一定数量的受试者，得到大量的样本，这些受试者中包含了心肌梗死患者以及正常人；(6)构建栈式自编码器，利用样本对该栈式自编码器进行训练，进而将栈式自编码器中的最深隐藏层作为支持向量机的输入，从而训练得到由栈式自编码器和支持向量机所组成的无创心脏梗死分类模型。进一步地，所述步骤(1)的具体实现方法为：利用磁共振仪器对受试者同时做冠、矢、轴三个方向的定位扫描，成像范围从心底及大血管根部到心尖部，从而得到受试者的完整心脏磁共振图像序列；同时，给受试者注射一定剂量的造影剂，通过扫描成像获得对应的一张延迟增强图像。进一步地，所述步骤(4)中对图像块的图像信息和运动信息进行融合，即将图像块的图像信息和运动信息进行归一化，进而使图像块中每一像素点归一化后的像素值和速度值相加。进一步地，所述步骤(5)中得到大量样本后，将所有样本分为训练集、验证集和测试集，其中训练集的样本用于训练栈式自编码器，验证集的样本用于对训练完成的栈式自编码器进行验证微调以确立最终的分类模型，测试集的样本用于对最终确立的分类模型进行测试验证。优选地，在训练集样本的融合信息中加入高斯噪声后作为训练栈式自编码器的输入，以使得模型具有抗噪能力，因为实际输入数据是采集到的真实磁共振图像，存在噪声，有必要在其训练过程中进行去噪。进一步地，所述步骤(6)中构建及训练栈式自编码器的具体方法如下：首先，训练多个级联的自编码器，所述自编码器由输入层、隐藏层和输出层组成，且前一个自编码器的隐藏层为后一个自编码器的输入层；对于第一个自编码器，使训练集样本的融合信息作为该自编码器的输入层，以该自编码器输出层与输入层的损失函数最小为目标，通过反向传播和梯度下降法不断计算更新该自编码器中的权重参数和偏置参数，直至第一个自编码器训练完成；对于第二个自编码器，使其前一个自编码器的隐藏层作为该自编码器的输入层，以该自编码器输出层与输入层的损失函数最小为目标，通过反向传播和梯度下降法不断计算更新该自编码器中的权重参数和偏置参数，直至第二个自编码器训练完成；以此类推，完成所有自编码器的训练后，将各自编码器的隐藏层依次堆叠在一起形成所述栈式自编码器，进而利用训练集样本的融合信息和类别信息分别作为栈式自编码器的输入以及输出真值并再次进行训练。进一步地，所述自编码器的函数模型如下：h＝g(wt+b)z＝f(w'h+b')其中：t、h和z分别为自编码器的输入层、隐藏层和输出层，w和b为输入层与隐藏层之间的权重参数和偏置参数，w'和b'为隐藏层与输出层之间的权重参数和偏置参数，g()为输入层与隐藏层之间的激活函数，f()为隐藏层与输出层之间的激活函数。本专利技术将深度学习的思想引入到心肌梗死分类中，其过程主要分为训练和预测两个阶段：在训练阶段先由采集到的心脏磁共振图像求得心脏的运动信息，将n×n的图像信息与运动信息融合的图像块作为输入，对应的梗死状况作为标签训练自编码器，利用降噪自编码器对输入数据做预处理，在学习高维信息过程中利用了上一层所有的变量因素，最终学得输入数据的深层特征；之后，以学得的深层特征作输入加上对应标签，利用支持向量机对样本分类；在测试阶段，将测试集中的数据X输入到堆栈式自编码器，提取隐层的深层特征，将深层特征输入到支持向量机，最终得到分类结果。因此，本专利技术从数据驱动的角度实现心肌梗死分类预测，有效地解决了临床上通过注射造影剂预测梗死费时费力问题。附图说明图1为本专利技术方法的整体框架示意图。图2为采集到的心脏磁共振电影序列。图3(a)为采集到的心脏磁共振延迟增强图像。图3(b)为心肌梗死区域的标注示意图。图4为第一个自编码器的结构示意图。图5为第二个自编码器的结构示意图。图6为栈式自编码器的结构示意图。具体实施方式为了更为具体地描述本专利技术，下面结合附图及具体实施方式对本专利技术的技术方案进行详细说明。本专利技术基于栈式自编码器和支持向量机的无创心脏梗死分类模型构建方法，整体框架如图1所示，具体实施步骤如下本文档来自技高网...

【技术保护点】
1.一种基于栈式自编码器和支持向量机的无创心脏梗死分类模型构建方法，包括如下步骤：(1)获取受试者的完整心脏磁共振图像序列及其对应的延迟增强图像；(2)通过对完整心脏磁共振图像序列进行筛选，只选取其中心脏短轴图像，从而得到对应的心脏短轴电影序列，进而利用optical flow计算出心脏短轴电影序列中与延迟增强图像所对应的一张心脏短轴图像M的运动信息；(3)从心脏短轴图像M中提取心肌区域对应的图像信息和运动信息，遍历心肌区域中所有n×n大小的图像块并利用延迟增强图像标记每个图像块的类别，即心肌梗死与否，n为大于1的自然数；(4)对图像块的图像信息和运动信息进行融合，将得到的融合信息以及图像块的类别信息组合成一个样本；(5)根据步骤(1)～(4)遍历一定数量的受试者，得到大量的样本，这些受试者中包含了心肌梗死患者以及正常人；(6)构建栈式自编码器，利用样本对该栈式自编码器进行训练，进而将栈式自编码器中的最深隐藏层作为支持向量机的输入，从而训练得到由栈式自编码器和支持向量机所组成的无创心脏梗死分类模型。

【技术特征摘要】
1.一种基于栈式自编码器和支持向量机的无创心脏梗死分类模型构建方法，包括如下步骤：(1)获取受试者的完整心脏磁共振图像序列及其对应的延迟增强图像；(2)通过对完整心脏磁共振图像序列进行筛选，只选取其中心脏短轴图像，从而得到对应的心脏短轴电影序列，进而利用opticalflow计算出心脏短轴电影序列中与延迟增强图像所对应的一张心脏短轴图像M的运动信息；(3)从心脏短轴图像M中提取心肌区域对应的图像信息和运动信息，遍历心肌区域中所有n×n大小的图像块并利用延迟增强图像标记每个图像块的类别，即心肌梗死与否，n为大于1的自然数；(4)对图像块的图像信息和运动信息进行融合，将得到的融合信息以及图像块的类别信息组合成一个样本；(5)根据步骤(1)～(4)遍历一定数量的受试者，得到大量的样本，这些受试者中包含了心肌梗死患者以及正常人；(6)构建栈式自编码器，利用样本对该栈式自编码器进行训练，进而将栈式自编码器中的最深隐藏层作为支持向量机的输入，从而训练得到由栈式自编码器和支持向量机所组成的无创心脏梗死分类模型。2.根据权利要求1所述的无创心脏梗死分类模型构建方法，其特征在于：所述步骤(1)的具体实现方法为：利用磁共振仪器对受试者同时做冠、矢、轴三个方向的定位扫描，成像范围从心底及大血管根部到心尖部，从而得到受试者的完整心脏磁共振图像序列；同时，给受试者注射一定剂量的造影剂，通过扫描成像获得对应的一张延迟增强图像。3.根据权利要求1所述的无创心脏梗死分类模型构建方法，其特征在于：所述步骤(4)中对图像块的图像信息和运动信息进行融合，即将图像块的图像信息和运动信息进行归一化，进而使图像块中每一像素点归一化后的像素值和速度值相加。4.根据权利要求1所述的无创心脏梗死分类模型构建方法，其特征在于：所述步骤(5)中得到大量样本后，将所有样本...

【专利技术属性】
技术研发人员：刘华锋，陈明强，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江,33