CN110403600A - 基于差值时间散点图的阵发性房颤智能分析方法及系统 - Google Patents

Abstract

本发明公开了一种基于差值时间散点图的阵发性房颤智能分析方法及系统，所述方法包括训练阶段和检测阶段；训练阶段包括：A1)采集心电图原始训练数据；A2)绘制训练阶段的差值时间散点图；A3)对差值时间散点图标注房颤发作标签和非房颤发作标签；A4)训练卷积神经网络模型；检测阶段包括：B1)采集心电图原始检测数据；B2)绘制检测阶段的差值时间散点图；B3)将检测阶段的差值时间散点图输入至卷积神经网络模型，输出阵发性房颤智能分析结果。本发明利用差值时间散点图分类的方法，可以实现阵发性房颤的智能诊断，从而实现对亚临床房颤患者的检出。

Description

基于差值时间散点图的阵发性房颤智能分析方法及系统

技术领域

本发明涉及医疗人工智能技术领域，具体地指一种基于差值时间散点图的阵发性房颤智能分析方法及系统。

背景技术

心房颤动(简称房颤)是最常见的快速性心律失常症状，流行病学研究显示，1990-2010年间临床房颤的患病率和发病率显著增加，2010年全球临床房颤人数估计有3350万，我国房颤患者的人数超过1000万，临床房颤相关死亡率10年间统计数据增加2倍。房颤可形成左心房附壁血栓，导致的栓塞90％是脑动脉栓塞(缺血性脑卒中)，10％是外周动脉栓塞或者肠系膜动脉栓塞等。据全球疾病负担研究(Global Burden of Disease study，GBD)2017年的统计数据，卒中所造成的伤残调整寿命年占比5.29％。亚临床房颤一般认为无或者几乎无典型临床症状的房颤，其中一部分因体检或其他原因就医行心电图检查而发现，另一部分直到房颤相关并发症如卒中出现后才被发现，甚至卒中发生后仍未能发现。亚临床房颤最早是在植入起搏器患者观测到的。研究发现在植入起搏器检出房颤的患者中，仅有17％-21％的患者表现出房颤发作的相关症状。同时阵发性房颤的患者中，使用5天动态心电图监测评估，只有10％的房颤发作会引起症状。房颤发作无症状的患者比房颤发作有症状的数量要高12倍，这使得这类患者房颤的检出面临挑战。亚临床房颤患者的人数是临床房颤患者的数倍，而亚临床房颤的危害与房颤相同。预防房颤相关血栓栓塞并发征尤其是卒中，早期检出房颤意义重大。

目前关于24小时动态心电监测在亚临床房颤的筛查中的作用的证据是有限的。通过动态心电监测亚临床房颤的功效取决于心电监护的持续时间和方法，与植入记录仪和CIEDs(敏感度91％)相比，24小时动态心电监测具有中等敏感性(44％-46％)，随着检测时间的延长可增加亚临床房颤的检出率。

植入心脏电子设备(CIEDs)的患者能够检测到房颤的发作，研究表明部分植入CIEDs的患者发作房颤时无任何不适症状，亦不会寻求治疗，血栓栓塞的风险增加的最小房颤持续时间目前还未得到精确的定义，它可能缩短至几分钟至几小时。至2012年大部分试验都支持大于6分钟的亚临床房颤负荷具有临床意义。目前植入电子设备对亚临床房颤的患者的诊断，仅限于原本已植入CIEDs的患者。有几项研究基于植入CIEDs的隐源性卒中的患者，虽然植入CIEDs不能完全自动检测房颤，但在非密集随访的条件下，这类设备房颤检出率更高。

英国国家卫生与临床技术优化研究所(NICE)评估了不同设备在房颤初级筛查中的成本效益，WachBP Home A和Alivercor heart monitor device两种可穿戴设备在房颤初级筛查和预防卒中方面比更具节省成本。最近的研究显示，使用智能手表、手环在识别房颤上技术是可行的，使用该类技术基于社区的筛查可发现高危人群中的房颤。美国心律学会HRS在2015年的远程监测共识声明中对远程心电监测在早期房颤的检测及量化的应用为ⅠA级推荐。但这些设备都不是实时连续监测心电，真正的无症状房颤的应用方面监测是否有价值是暂待商榷。动态心电监测对亚临床房颤的诊断取决于检测的时间，然而随着检测时间的增加，既增加动态心电分析的工作量，同时也增加检测的费用。而植入设备昂贵，且为侵入式检查。而目前的心电可穿戴设备多为间断监测心电，因此以上三种检测亚临床房颤的方式均有明显的缺点。

基于亚临床房颤患者心律的异常可能在长时间的检测中暴露出来的特点。长时间连续心电检测和分析需要多天、长期监测，人工分析的难度和工作强度很大，误诊率也较高，同时房颤负荷无法具体量化。因此针对亚临床房颤患者检出需制定一整套心电监测方案，能够达到智能诊断，并且能具体量化房颤负荷，迫在眉睫。

发明内容

针对现有技术的不足之处，本发明提出了一种基于差值时间散点图的阵发性房颤智能分析方法及系统，通过人工智能技术实现阵发性房颤的智能分析和判断，从而实现对亚临床房颤患者的检出预测分析。

为实现上述目的，本发明提出一种基于差值时间散点图的阵发性房颤智能分析方法，其特殊之处在于，所述方法包括训练阶段和检测阶段；

所述训练阶段包括步骤：

A1)采集若干份长时动态心电图的原始训练数据，并过滤质量过低的信号段；

A2)根据每一份所述心电图的原始训练数据绘制训练阶段的差值时间散点图；所述差值时间散点图的绘制方法为：根据原始心电图数据标记所有的R波顶点；将取得的相邻两个R波顶点之间的间隔作为一个RR间期的时间宽度；对得到的心电图数据RR间期宽度序列进行遍历，以当前RR间期宽度M减去下一个RR间期宽度N得到的数据(M-N)为纵坐标，以时间纬度t为横坐标，绘制散点图(M-N，t)；

A3)对每一份所述训练阶段的差值时间散点图标注房颤发作标签和非房颤发作标签；

A4)将所有所述训练阶段的差值时间散点图和对应的房颤发作标签和非房颤发作标签作为卷积神经网络的输入，训练得到卷积神经网络模型；

所述检测阶段包括步骤：

B1)采集待检测的长时动态心电图的原始检测数据，并过滤质量过低的信号段；

B2)根据所述心电图的原始检测数据绘制检测阶段的差值时间散点图；

B3)将所述检测阶段的差值时间散点图输入至卷积神经网络模型，计算得到检测散点图的分类结果，输出阵发性房颤智能分析结果。

优选地，所述步骤A3)中对具有房颤时间段点集的心电图标注房颤发作标签，对具有非房颤时间段点集的心电图标注非房颤发作标签。

优选地，所述房颤时间段点集为围绕时间轴上下均匀分布、无中心连续分布的点集。

优选地，所述非房颤时间段点集为具有中心线装、呈连续分布或多条线状的点集。

优选地，所述步骤A4)之后还包括验证步骤：将样本原始心电图波形数据绘制为样本的差值时间散点图，输入经过训练的卷积神经网络模型输出样本心电图的阵发性房颤智能分析结果，所述结果的正确率高于设定值则前往检测阶段，否则返回步骤A1)。

优选地，所述长时动态心电图数据为连续采集的超过6小时的动态心电图数据。

本发明还提出一种上述基于差值时间散点图的阵发性房颤智能分析方法的系统，包括存储器、处理器、数据解析模块、数据分析模块；

所述数据解析模块：用于对采集长时动态心电图数据，并对原始数据进行格式转换；

所述数据分析模块：用于根据卷积神经网络模型，将输入的原始数据绘制为差值时间散点图，进行计算、分类，输出阵发性房颤智能分析结果。

进一步地，所述数据解析模块包括数据读取模块、图像格式转换模块和采样率转换模块；

所述数据读取模块：用于与心电图数据采集设备连接，采集长时动态心电图数据；

所述图像格式转换模块：用于将图像数据转换为CSV数据格式；

所述采样率转换模块：用于将原始检测数据的采样频率转换为与原始训练数据的采样频率相一致。

更进一步地，所述数据分析模块包括数据预处理模块、差值时间散点图生成模块和标注分类模块；

所述数据预处理模块：用于对原始检测数据进行质量评估，并过滤质量过低的信号；

所述散点图生成模块：用于根据原始检测数据生成差值时间散点图；

所述标注分类模块：用于根据卷积神经网络模型利用深度学习算法对差值时间散点图进行分类，从而实现对心电图的心律失常类型的分析判断，输出阵发性房颤智能分析结果。

在机器学习中，卷积神经网络(Convolutional Neural Network,CNN)是一种前馈神经网络，其人工神经元可以响应一部分覆盖范围内的周围单元。卷积神经网络是近年发展起来，并引起广泛重视的一种高效识别方法。一般地，CNN的基本结构包括两个阶段，每个阶段都包含若干层次，其一为特征提取阶段，每个神经元的输入与前一层的局部接受域相连，并提取该局部的特征。一旦该局部特征被提取后，它与其它特征间的位置关系也随之确定下来；其二是特征映射阶段，网络的每个计算层由多个特征映射组成，每个特征映射是一个平面，平面上所有神经元的权值相等。因为CNN的特征检测层通过训练数据进行学习，所以在使用CNN时，避免了显式的特征抽取，而隐式地从训练数据中进行学习；再者因为同一特征映射面上的神经元权值相同，所以网络可以并行学习。同一特征映射面上的神经元权值相同，所以网络可以并行学习。

差值时间散点图代表随着时间变化的RR间期序列的变化，每种心律都具有特征性内在的相互关联性，房颤的心电生理特性可以做RR间期显示出来。所以本发明利用机器学习技术学习时间散点图的特性，对心律失常进行初步的预判，针对预判的结果再进行分析：一方面对海量的心电数据可以进行快速的预判，另一方面强化了动态心电的诊断，让复杂的心律失常诊断更为智能，更为准确。本发明提出的训练阶段，通过区分不同时间段的阵发性房颤的时间散点图，将发作房颤的时间段和非发颤发作的时间段进行归类，训练卷积神经网络分类模型。

本发明利用差值时间散点图分类的方法，可以对有亚临床房颤的患者进行长时间的检测分析，智能分析其有无阵发房颤的发作，从而让阵发房颤得以证实，提高阵发性房颤的检出率，同时降低假阳性率，早干预及治疗，降低全社会临床房颤的危害；另一方面也降低亚临床房颤的检测费用，同时也降低亚临床房颤并发症的经济负担和社会负担。

附图说明

图1本发明基于差值时间散点图的阵发性房颤智能分析方法的流程图。

图2为本发明基于差值时间散点图进行分类的深度学习模型训练和测试的原理图。

图3为主导节律为无房颤发作的心电图对应的差值时间散点图示例。

图4为主导节律为持续房颤发作的心电图对应的差值时间散点图示例。

图5为主导节律为阵发性房颤的心电图对应的差值时间散点图示例。

图6本发明基于差值时间散点图的阵发性房颤智能分析方法的系统的结构框图。

图7为本发明提出的用于对基于差值时间散点图进行分类的深度学习模型的一个实施案例的网络结构。

具体实施方式

下面结合附图及实施例对本发明作进一步的详细描述，但该实施例不应理解为对本发明的限制。

本发明提出的一种基于差值时间散点图的阵发性房颤智能分析方法，如图1和图2所示，包括训练阶段和检测阶段；

训练阶段包括步骤：

A1)采集若干份长时动态心电图的原始训练数据，并过滤质量过低的信号段。

获取大量长时阵发心房颤动的动态心电图的原始数据，并过滤质量过低的信号段。采集大量动态心电图样本，采样频率为500Hz，每个样本时长为24小时，基于心电信号波形的形态、频域特征以及其他非线性分析方法，分段评估心电信号的质量，并过滤质量过低的信号段。

A2)根据每一份心电图的原始训练数据绘制训练阶段的差值时间散点图：根据原始心电图数据标记所有的R波顶点。用基于阈值和小波分析的方法(或其他有效手段)提取R波顶点位置。利用R波顶点数据，得到RR间期宽度序列。相邻两个R波顶点之间的间隔就是一个RR间期的时间宽度。对得到的RR间期宽度序列进行遍历，以当前RR间期宽度(单位为毫秒，取整数)为横坐标下一个RR间期宽度为纵坐标，绘制散点图。图3、图4、图5均为心电散点图的例图。

A3)对每一份所述训练阶段的差值时间散点图标注房颤发作标签和非房颤发作标签。

涵盖的心律失常类型标签具体包括：

1)房颤发作时间段标签房颤发作标签标注于具有房颤时间段点集的心电图，即具有围绕时间轴上下均匀分布、无中心连续分布特征的点集。

2)非房颤发作时间段标签非房颤发作标签标注于具有非房颤时间段点集的心电图，即具有中心线装、呈连续分布或多条线状特征的点集。

根据差值时间散点图形态区别房颤发作时间段点集和非房颤发作时间段点集：

a)房颤发作时间段点集

如图3所示，心房颤动的时间段点集，形态为以围绕时间轴上下(y＝0)分布散在间段点集

如图4所示，非心房颤动时间段点集，通常存在中心线状点集，呈连续分布，或多条线状点集。

A4)将所有训练阶段的差值时间散点图和对应的房颤发作标签和非房颤发作标签作为卷积神经网络(如图7所示)的输入，训练得到卷积神经网络模型。

将样本原始心电图波形数据绘制为样本的差值时间散点图，输入经过训练的卷积神经网络模型输出样本心电图的阵发性房颤智能分析结果，所述结果的正确率高于设定值则前往检测阶段，否则返回步骤A1)。

检测阶段包括步骤：

B1)采集待检测的长时动态心电图的原始检测数据，并过滤质量过低的信号段；获取长时动态心电图，每个样本时长为24小时，采样频率为500Hz(同训练阶段)，若频率与训练阶段预设的采样频率不同，则采用插值法转换到相同频率。同模型训练阶段的步骤A1)，基于心电信号波形的形态、频域特征以及其他非线性分析方法，分段评估心电信号的质量，并过滤质量过低的信号段。

B2)根据心电图的原始检测数据绘制检测差值时间散点图。

使用与模型训练阶段的步骤A2)相同的方式提取R波定点、得到RR间期宽度序列，进而获得差值时间散点图。

B3)利用通过模型训练阶段的步骤A4)得到训练完成的卷积神经网络模型，将B2)中得到的差值时间散点图作为卷积神经网络模型的输入，计算得到散点图的分类结果，即心电图的心律失常类型。

在模型训练阶段，深度学习框架采用TensorFlow，卷积神经网络的学习过程如下：

1)将采集并处理的长时动态心电图绘制成差值时间散点图之后，将差值时间散点图和房颤发作标签及非房颤发作标签作为训练数据存储于计算机中，设定迭代次数上限；

2)按图7构建卷积神经网络的结构，初始化参数，生成各卷积层、池化层、全连接层的权值；

3)根据当前的权值参数按照网络结构从前往后逐层计算，得到输出结果；

4)根据现有原始标签的数据，对整个神经网络的参数进行调整优化，更新各层的权值；

5)迭代次数自增，若达到迭代次数上限，则进入步骤6)，否则进入步骤3)；

6)卷积神经网络的各层参数训练完毕，得到卷积神经网络模型。

基于上述方法，本发明还提出一种基于差值时间散点图的阵发性房颤智能分析方法的系统，如图6所示，包括存储器、处理器、数据解析模块、数据分析模块。

数据解析模块：用于对采集长时动态心电图数据，并对原始数据进行格式转换。

数据分析模块：用于根据卷积神经网络模型，将输入的原始数据绘制为差值时间散点图，进行计算、分类，输出阵发性房颤智能分析结果。

数据解析模块包括数据读取模块、图像格式转换模块和采样率转换模块：数据读取模块：用于与心电图数据采集设备连接，采集长时动态心电图数据；图像格式转换模块：用于将图像数据转换为CSV数据格式；采样率转换模块：用于将原始检测数据的采样频率转换为与原始训练数据的采样频率相一致。

数据分析模块包括数据预处理模块、差值时间散点图生成模块和标注分类模块：数据预处理模块：用于对原始检测数据进行质量评估，并过滤质量过低的信号；散点图生成模块：用于根据原始检测数据生成差值时间散点图；标注分类模块：用于根据卷积神经网络模型利用深度学习算法对差值时间散点图进行分类，从而实现对心电图的心律失常类型的分析判断，输出阵发性房颤智能分析结果。

尽管上面结合附图对本发明的优选实施例进行了描述，但是本发明并不局限于上述的具体实方式，上述的具体实施方式仅仅是示意性的，并不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可以作出很多形式的具体变换，这些均属于本发明的保护范围内。

Claims (9)

1.一种基于差值时间散点图的阵发性房颤智能分析方法，其特征在于：所述方法包括训练阶段和检测阶段；

所述训练阶段包括步骤：

A1)采集若干份长时动态心电图的原始训练数据，并过滤质量过低的信号段；

A2)根据每一份所述心电图的原始训练数据绘制训练阶段的差值时间散点图；所述差值时间散点图的绘制方法为：根据原始心电图数据标记所有的R波顶点；将取得的相邻两个R波顶点之间的间隔作为一个RR间期的时间宽度；对得到的心电图数据RR间期宽度序列进行遍历，以当前RR间期宽度M减去下一个RR间期宽度N得到的数据(M-N)为纵坐标，以时间纬度t为横坐标，绘制散点图(M-N，t)；

A3)对每一份所述训练阶段的差值时间散点图标注房颤发作标签和非房颤发作标签；

A4)将所有所述训练阶段的差值时间散点图和对应的房颤发作标签和非房颤发作标签作为卷积神经网络的输入，训练得到卷积神经网络模型；

所述检测阶段包括步骤：

B1)采集待检测的长时动态心电图的原始检测数据，并过滤质量过低的信号段；

B2)根据所述心电图的原始检测数据绘制检测阶段的差值时间散点图；

B3)将所述检测阶段的差值时间散点图输入至卷积神经网络模型，计算得到检测散点图的分类结果，输出阵发性房颤智能分析结果。

2.根据权利要求1所述的基于差值时间散点图的阵发性房颤智能分析方法，其特征在于：所述步骤A3)中对具有房颤时间段点集的心电图标注房颤发作标签，对具有非房颤时间段点集的心电图标注非房颤发作标签。

3.根据权利要求2所述的基于差值时间散点图的阵发性房颤智能分析方法，其特征在于：所述房颤时间段点集为围绕时间轴上下均匀分布、无中心连续分布的点集。

4.根据权利要求2所述的基于差值时间散点图的阵发性房颤智能分析方法，其特征在于：所述非房颤时间段点集为具有中心线装、呈连续分布或多条线状的点集。

5.根据权利要求1所述的基于差值时间散点图的阵发性房颤智能分析方法，其特征在于：所述步骤A4)之后还包括验证步骤：将样本原始心电图波形数据绘制为样本的差值时间散点图，输入经过训练的卷积神经网络模型输出样本心电图的阵发性房颤智能分析结果，所述结果的正确率高于设定值则前往检测阶段，否则返回步骤A1)。

6.根据权利要求1所述的基于差值时间散点图的阵发性房颤智能分析方法，其特征在于：所述长时动态心电图数据为连续采集的超过6小时的动态心电图数据。

7.一种根据权利要求1～6中任一权利要求所述的基于差值时间散点图的阵发性房颤智能分析方法的系统，其特征在于：包括存储器、处理器、数据解析模块、数据分析模块；

所述数据解析模块：用于对采集长时动态心电图数据，并对原始数据进行格式转换；

所述数据分析模块：用于根据卷积神经网络模型，将输入的原始数据绘制为差值时间散点图，进行计算、分类，输出阵发性房颤智能分析结果。

8.根据权利要求7所述的基于差值时间散点图的阵发性房颤智能分析方法的系统，其特征在于：所述数据解析模块包括数据读取模块、图像格式转换模块和采样率转换模块；

所述数据读取模块：用于与心电图数据采集设备连接，采集长时动态心电图数据；

所述图像格式转换模块：用于将图像数据转换为CSV数据格式；

所述采样率转换模块：用于将原始检测数据的采样频率转换为与原始训练数据的采样频率相一致。

9.根据权利要求7所述的基于差值时间散点图的阵发性房颤智能分析方法的系统，其特征在于：所述数据分析模块包括数据预处理模块、差值时间散点图生成模块和标注分类模块；

所述数据预处理模块：用于对原始检测数据进行质量评估，并过滤质量过低的信号；

所述散点图生成模块：用于根据原始检测数据生成差值时间散点图；

所述标注分类模块：用于根据卷积神经网络模型利用深度学习算法对差值时间散点图进行分类，从而实现对心电图的心律失常类型的分析判断，输出阵发性房颤智能分析结果。