申请/专利权人：复旦大学

申请日：2019-07-01

公开（公告）日：2019-10-25

公开（公告）号：CN110367973A

主分类号：A61B5/0472(20060101)

分类号：A61B5/0472(20060101);A61B5/00(20060101)

优先权：

专利状态码：有效-授权

法律状态：2021.10.26#授权;2019.11.19#实质审查的生效;2019.10.25#公开

摘要：本发明涉及一种多导联心电信号QRS波群起始点的检测方法。利用心电采集技术采集多导联心电信号；将多导联心电信号预处理得到去除基线漂移及高频噪声后的信号；对该信号作归一化处理，得到归一化的多导联心电信号；计算归一化后的多导联心电信号的均方根RMS及检测函数DEF，通过学习RMS和DEF前若干秒的值计算得到阈值，将满足阈值条件的DEF峰值点视为QRS波群的起始点。本发明方法有检测精确率高、算法简单和操作实时等优点。本发明方法适用于长程或短程心电信号，可用于分析窦性心律信号或心律失常信号，同时在电生理机制研究和临床医学上均具有一定的应用价值。本发明方法可推广至心内电信号的激动分析和相关量化研究。

主权项：1.一种多导联心电信号QRS波群起始点的检测方法，其特征在于具体步骤如下：1利用心电采集技术采集多导联心电信号；2对步骤1的多导联心电信号进行预处理，得到去除基线漂移和高频噪声的多导联心电信号；3对步骤2得到的预处理后的多导联心电信号进行归一化处理，得到多导联归一化心电信号；4计算步骤3得到的多导联归一化心电信号的均方根RMSRootMeanSquare，得到RMS信号，假设有n导联心电信号n≥2，均方根计算如公式1所示： 其中：ECGi是第i导联心电信号的幅值，1≤i≤n；5利用步骤4中得到的RMS信号计算检测函数DEFDetectionFunction信号，DEF计算如公式2所示： 对于步骤4得到的心电信号的RMS信号，设o点是待检测的QRS波的起始点，width是时间参数经验值，k和m均为正常数经验值；其中：a为o点向前width时间内出现的最低点与o点的幅度差，若o点为拐点，则a值接近为0，exp-ka接近于1；若o点为上升支的某点，则a为正值，exp-ka的值接近于0，这有利于区分o点与其他处于RMS信号上升支的点；b为o点向后width时间内出现的最高点与o点的幅度差，其主要反映了RMS信号上升的速度，有利于区分o点与其前面的平坦区段，同时有利于区分o点与P波、T波的起始点；c为o点向后m*width时间内出现的RMS信号的最大值；它有利于区分QRS波、P波和T波；d为o点与向后m*width时间内出现的RMS最大值点之间的时间差；室性早搏的c值较大，但是d值较小；而窦性心律的c值较小，d值较大；因此综合考虑c和d有利于减小不同心拍DEF检测函数值的差异；同时，d值也有利于将QRS波与部分缓慢上升的T波区分开；e为o点到o点向后m*width时间内的RMS积分值；6利用步骤4得到的均方根信号的前若干秒长的数据，求取该前若干秒长内均方根信号的峰值，构成一组序列，取该序列的中值的某百分比对应的数值作为RMS阈值；利用步骤5得到的检测函数信号的前若干秒长的数据，求取该前若干秒长内检测函数信号的峰值，构成一组序列，取该序列的中值的某百分比对应的数值作为DEF阈值；将满足阈值条件，并且不在预设的不应期200～400ms可调内的DEF峰值点视为QRS起始点；若2s内未能检测到QRS起点，则放宽阈值条件，重新检测，将高于该阈值的某百分比对应的数值的DEF峰值点视为QRS的起始点。

全文数据：多导联心电信号QRS波群起始点的检测方法技术领域本发明涉及一种多导联心电信号QRS波群起始点的检测方法。背景技术在心电周期信号中，QRS波群的幅度最大、特征最为明显，故在现有的心电信号特征点检测算法中，有关QRS波群的检测方法很多，常见的方法有：模板匹配法、差分阈值法、人工神经网络法和小波变换法。模板匹配法需根据不同个体的差异性做出相应的调整，算法依赖性强，且易受噪声干扰，检测结果往往不理想。当发生心律不齐，尤其是室性早搏时，QRS波群的R波幅度变化很大，且RR间期不规律，传统的差分阈值法无法准确定位。人工神经网络的方法虽有较高的自适应性及较好的检测效果，但该方法的实现需要大量具有代表性的样本，在实际应用过程中较难实现。小波变换法计算量大、步骤复杂、耗时长，且无法定位起、止点。以上的QRS波群检测方法主要是用于提取R波，而关于检测QRS波群的起点及终点等特征点的检测却很简单，主要基于R波峰值点向前或向后搜寻得到。实际上，QRS波群的起始点与终止点是正确表征心脏电活动的关键点，对于临床上诊断心源性异常如休克、心绞痛、心室传导阻滞等有重要的意义。发明内容针对传统特征点检测算法的实时性差、算法复杂、准确率低等特点，本发明的目的在于提出一种多导联心电信号QRS波群起始点的检测方法，本发明能够通过一系列信号预处理步骤，得到归一化的多导联心电信号，并通过计算归一化心电信号的均方根和检测函数，最终检测出多导联心电信号QRS波群的平均起始点。本发明提出的一种多导联心电信号QRS波群起始点的检测方法，具体步骤如下：1利用心电采集技术采集多导联心电信号；2对步骤1的多导联心电信号进行预处理，得到去除基线漂移和高频噪声的多导联心电信号；3对步骤2得到的预处理后的多导联心电信号进行归一化处理，得到多导联归一化心电信号；4计算步骤3得到的多导联归一化心电信号的均方根RMS，得到RMS信号，假设有n导联心电信号n≥2，均方根计算如公式1所示：其中：ECGi是第i导联心电信号的幅值，1≤i≤n；5利用步骤4中得到的RMS信号计算检测函数DEFDetectionFunction，DEF计算如公式2所示：对于步骤4得到的心电信号的RMS信号，设o点是待检测的QRS波的起始点，width是时间参数经验值，k和m均为正常数经验值；其中：a为o点向前width时间内出现的最低点与o点的幅度差，若o点为拐点，则a值接近为0，exp-ka接近于1；若o点为上升支的某点，则a为正值，exp-ka的值接近于0，这有利于区分o点与其他处于RMS信号上升支的点；b为o点向后width时间内出现的最高点与o点的幅度差，其主要反映了RMS信号上升的速度，有利于区分o点与其前面的平坦区段，同时有利于区分o点与P波、T波的起始点；c为o点向后m*width时间内出现的RMS信号的最大值，它有利于区分QRS波、P波和T波；d为o点与向后m*width时间内出现的RMS最大值点之间的时间差，一般而言，室性早搏的c值较大，但是d值较小；而窦性心律的c值较小，d值较大；因此综合考虑c和d有利于减小不同心拍DEF检测函数值的差异，同时，d值也有利于将QRS波与部分缓慢上升的T波区分开；e为o点到o点向后m*width时间内的RMS积分值；6利用步骤4得到的均方根信号的前若干秒长的数据，求取该前若干秒长内均方根信号的峰值，构成一组序列，取该序列的中值的某百分比对应的数值作为RMS阈值；利用步骤5得到的检测函数信号的前若干秒长的数据，求取该前若干秒长内检测函数信号的峰值，构成一组序列，取该序列的中值的某百分比对应的数值作为DEF阈值；将满足阈值条件，并且不在预设的不应期200～400ms可调内的DEF峰值点视为QRS起始点；若2s内未能检测到QRS起点，则放宽阈值条件，重新检测，将高于该阈值的某百分比对应的数值的DEF峰值点视为QRS的起始点。本发明具有以下有益效果：1.本发明提出了一种多导联心电信号QRS波群起始点检测方法，与传统的检测方法相比，具有检测精确率较高、检测实时性强以及算法简便等优点。2.本发明对于多形态室性早搏、房颤、ST段改变、二联律或其他节律不规则的异常心律的多导联心电信号QRS波群起始点检测均具有一定的适应性。3.本方法适用于短程或长程的心电信号，同时适用于窦性心律信号和心律失常信号，为心电诊断提供了一个重要的检测方法。4.检测到起始点后，分别对多导联心电信号进行积分，得出的结果可表征心电综合矢量于一个心动周期内在各个导联的一个平均投影矢量，有助于研究室性早搏等心律失常信号。5.本方法可推广至心内电信号的激动分析和相关量化研究。附图说明为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单的介绍。需要说明的是，以下附图仅展示出了本发明的某些实施例，因此不应被看作是范围的限定。图1是体表12导联心电信号的均方根信号的示意图。其中，o点为待检测的QRS波起始点；width是设定的时间参数；a为o点向前width时间内出现的最低点与o点的幅度差；b为o点向后width时间内出现的最高点与o点的幅度差；c为o点向后4*width时间内出现的RMS信号的最大值；d为o点与向后4*width时间内出现的RMS最大值点之间的时间差；e为o点与o点向后4*width时间内的均方根积分值。图2是窦性心律的体表12导联心电信号QRS波群的起始点检测示意图。其中，a为一段10s长的12导联原始心电信号，I、II、III为标准肢体导联，aVR、aVL、aVF为加压肢体导联，V1～V6为6个胸导联；b为该段12导联心电信号滤波后的信号，图中标识与a相同；c为该段12导联心电信号的归一化信号，纵坐标为其归一化幅度，横坐标为时间；d为12导联归一化信号的均方根值，纵坐标为其幅度，横坐标为时间，横虚线表示均方根阈值；e为该段12导联心电信号的检测函数值，纵坐标为其幅度，横坐标为时间，横虚线表示检测函数阈值；f为最终的检测效果，前12行信号为图b，第13行信号为图d，第14行信号为图e，竖虚线为标记的R波，竖实线为检测到的12导联心电信号平均QRS波群起始点，每条竖实线顶端的数字表示心拍数。图3是多形态室性早搏的体表12导联心电信号QRS波群的起始点检测示意图。其中，a～f的图片布局与图2中对应的a～f相同。图4是房颤的体表12导联心电信号QRS波群的起始点检测示意图。其中，a～f的图片布局与图2中对应的a～f相同。图5是ST段压低的体表12导联心电信号QRS波群的起始点检测结果示意图。图片布局与图2f相同。图6是阵发性室速的体表12导联心电信号QRS波群的起始点检测结果示意图。图片布局与图2f相同。图7是室性早搏二联律的体表12导联心电信号QRS波群的起始点检测结果示意图。图片布局与图2f相同。具体实施方式下面基于常规的体表12导联心电数据即导联数n＝12结合附图和实施例对本发明方法及应用做进一步说明。实施例1：将本发明的QRS波群起始点的检测方法应用于窦性心律信号不同部位节律一致。本实施例利用采样率为1000Hz的体表12导联心电信号即导联数n＝12，工作流程如下：1观察同步采集到的窦性心律的体表12导联心电信号数字信号，如图2a所示。每个导联的浅灰色横线表示其零电位线。由图可知有些导联的心电信号存在严重的基线漂移情况，如I、II、III、aVR、aVL、aVF、V2、V3、V4、V5、V6。2对步骤1的窦性心律的体表12导联心电信号进行预处理。首先以sym4小波基函数，对信号作10层分解，对于该1000Hz的心电信号，将第10层的分解信号频段为0～0.98Hz视为基线，将其置为0；接着对第1～2层的心电信号125～500Hz进行硬阈值处理去除高频噪声。经上述预处理后得到的体表12导联心电信号如图2b所示。3为保证各个导联对最终检测结果的贡献相同，对步骤2得到的预处理后的12导联心电信号归一化，得到如图2c所示信号。4计算步骤3所得归一化12导联心电信号的均方根，得到如图2d所示信号。求取前10s内均方根信号的峰值，构成一个序列，取该序列的中值的二分之一作为均方根阈值，如图中横虚线所示。5计算步骤4所得均方根信号的检测函数令计算公式2中的k＝5,width＝25ms,m＝4，得到如图2e所示信号。求取前10s内检测函数的峰值，构成一个序列，取该序列的中值的二分之一作为检测函数阈值，如图中横虚线所示。6预设300ms的不应期，将同时满足以下三个条件的点视为QRS波群的起始点：a.步骤5得到的检测函数中大于检测函数阈值的峰值点；b.步骤4得到的均方根信号中，符合条件a的点后100ms内存在大于均方根阈值的均方根峰值点；c.符合条件a和b的点在预设的不应期之外。检测结果如图2f所示，图中竖实线代表着12导联心电信号QRS波群的起始点。从图中可以观察到该方法有效地避开了P波和T波的起始点，准确地检测出了有规律的12导联窦性心律信号QRS波群的平均起始点。实施例2：将本发明的QRS波群起始点的检测方法应用于室性早搏信号在窦房结冲动尚未抵达心室之前，心室中的任何一个部位或室间隔的异位节律点提早发出电冲动，从而引起心室的除极。本实施例利用采样率为1000Hz的体表12导联心电信号即导联数n＝12，工作流程如下：1观察同步采集到的多形态室性早搏的体表12导联心电信号数字信号，如图3a所示。每个导联的浅灰色横线表示其零电位线，从图可知有些导联的心电信号存在严重的基线漂移情况，如II、III、aVR、aVL、AVF、V1、V3、V4、V5、V6。2对步骤1的室性早搏的12导联心电信号进行预处理。首先以sym4小波基函数，对该心电信号作10层分解，对于该1000Hz的心电信号，将第10层的分解信号频段为0～0.98Hz视为基线，将其置为0；接着对第1～2层的心电信号125～500Hz进行硬阈值处理去除高频噪声。经上述预处理后得到的12导联心电信号如图3b所示。3为保证各个导联对最终检测结果的贡献相同，对步骤2得到的预处理后的12导联心电信号归一化，得到如图3c所示信号。4计算步骤3所得归一化12导联心电信号的均方根，得到如图3d所示信号。求取前10s内均方根信号的峰值，构成一个序列，取该序列的中值的二分之一作为均方根阈值，如图中横虚线所示。5计算步骤4所得均方根信号的检测函数令计算公式2中的k＝5,width＝25ms,m＝4，得到如图3e所示信号。求取前10s内检测函数的峰值，构成一个序列，取该序列的中值的二分之一作为检测函数阈值，如图中横虚线所示。6预设300ms的不应期，将同时满足以下三个条件的点视为QRS波群的起始点：a.步骤5得到的检测函数中大于检测函数阈值的峰值点；b.步骤4得到的均方根信号中，符合条件a的点后100ms内存在大于均方根阈值的均方根信号峰值点；c.符合条件a和b的点在预设的不应期之外。检测结果如图3f所示，图中竖实线代表着12导联心电信号QRS波群的起始点。从图中可以观察得到，尽管室性早搏的心电图中出现了高大畸形的QRS波，但该方法仍然有效准确地检测出了异常的12导联心电信号QRS波群的平均起始点。实施例3：将本发明的QRS波群起始点的检测方法应用于房颤信号存在多个节律，不同部位节律可能不一致。本实施例利用采样率为1000Hz的体表12导联心电信号即导联数n＝12，工作流程如下：1观察同步采集到的房颤的体表12导联心电信号数字信号，如图4a所示。每个导联的浅灰色横线表示其零电位线，可知有些导联的心电信号存在严重的基线漂移情况，如I、II、aVR、aVL、V3、V5、V6。2对步骤1的房颤的12导联心电信号进行预处理。首先以sym4小波基函数，对信号作10层分解，对于该1000Hz的心电信号，将第10层的分解信号频段为0～0.98Hz视为基线，将其置为0；接着对第1～2层的心电信号125～500Hz进行硬阈值处理去除高频噪声。经上述预处理后得到的12导联心电信号如图4b所示。3为保证各个导联对最终检测结果的贡献相同，对步骤2得到的预处理后的12导联心电信号归一化，得到如图4c所示信号。4计算步骤3所得的归一化12导联心电信号的均方根，得到如图4d所示信号。求取前10s内均方根信号的峰值，构成一个序列，取该序列的中值的二分之一作为均方根阈值，如图中横虚线所示。5计算步骤4所得均方根信号的检测函数令计算公式2中的k＝5,width＝25ms,m＝4，得到如图4e所示信号。求取前10s内检测函数的峰值，构成一个序列，取该序列的中值的二分之一作为检测函数阈值，如图中横虚线所示。6预设300ms的不应期，将同时满足以下三个条件的点视为QRS波群的起始点：a.步骤5得到的检测函数中大于检测函数阈值的峰值点；b.步骤4得到的均方根信号中，符合条件a的点后100ms内存在大于均方根阈值的均方根信号峰值点；c.符合条件a和b的点在预设的不应期之外。检测结果如图4f所示，图中竖实线代表着12导联心电信号QRS波群的起始点。从图中可以观察得到，尽管房颤的心电图中P波消失，代之以锯齿形状的f波，但该方法仍然有效准确地检测出了异常的12导联心电信号QRS波群的平均起始点。

权利要求：1.一种多导联心电信号QRS波群起始点的检测方法，其特征在于具体步骤如下：1利用心电采集技术采集多导联心电信号；2对步骤1的多导联心电信号进行预处理，得到去除基线漂移和高频噪声的多导联心电信号；3对步骤2得到的预处理后的多导联心电信号进行归一化处理，得到多导联归一化心电信号；4计算步骤3得到的多导联归一化心电信号的均方根RMSRootMeanSquare，得到RMS信号，假设有n导联心电信号n≥2，均方根计算如公式1所示：其中：ECGi是第i导联心电信号的幅值，1≤i≤n；5利用步骤4中得到的RMS信号计算检测函数DEFDetectionFunction信号，DEF计算如公式2所示：对于步骤4得到的心电信号的RMS信号，设o点是待检测的QRS波的起始点，width是时间参数经验值，k和m均为正常数经验值；其中：a为o点向前width时间内出现的最低点与o点的幅度差，若o点为拐点，则a值接近为0，exp-ka接近于1；若o点为上升支的某点，则a为正值，exp-ka的值接近于0，这有利于区分o点与其他处于RMS信号上升支的点；b为o点向后width时间内出现的最高点与o点的幅度差，其主要反映了RMS信号上升的速度，有利于区分o点与其前面的平坦区段，同时有利于区分o点与P波、T波的起始点；c为o点向后m*width时间内出现的RMS信号的最大值；它有利于区分QRS波、P波和T波；d为o点与向后m*width时间内出现的RMS最大值点之间的时间差；室性早搏的c值较大，但是d值较小；而窦性心律的c值较小，d值较大；因此综合考虑c和d有利于减小不同心拍DEF检测函数值的差异；同时，d值也有利于将QRS波与部分缓慢上升的T波区分开；e为o点到o点向后m*width时间内的RMS积分值；6利用步骤4得到的均方根信号的前若干秒长的数据，求取该前若干秒长内均方根信号的峰值，构成一组序列，取该序列的中值的某百分比对应的数值作为RMS阈值；利用步骤5得到的检测函数信号的前若干秒长的数据，求取该前若干秒长内检测函数信号的峰值，构成一组序列，取该序列的中值的某百分比对应的数值作为DEF阈值；将满足阈值条件，并且不在预设的不应期200～400ms可调内的DEF峰值点视为QRS起始点；若2s内未能检测到QRS起点，则放宽阈值条件，重新检测，将高于该阈值的某百分比对应的数值的DEF峰值点视为QRS的起始点。

百度查询： 复旦大学 多导联心电信号QRS波群起始点的检测方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD