一种结合心电时频域特征分析的可电击复律心律自动识别和归类方法【专利摘要】本发明提供了一种结合心电时频域分析的可电击复律心律自动识别和归类方法，具体步骤：S1.心电信号的预处理；S2.心脏停搏心律的自动识别，若不满足判别条件，则执行步骤S3；即基于整系数带通滤波器，计算输出信号的最大幅度比例值（Pa）、平均幅度比例值（Pb）和平均偏差比例值（Pc）；步骤S4、S5、S6和S7则根据Pa、Pb和Pc等频域特征值判别可电击复律心律和不可电击复律心律，若失败，则执行步骤S8；即计算心电标准化栅条投影标准差；S9.根据标准差值判别可电击复律心律和不可电击复律心律。本发明可用于根据体表心电图自动识别和归类可电击复律心律的仪器设备，提高其可电击复律心律的识别灵敏度和不可电击复律心律的特异性，增强算法的运行效率。【专利说明】-种结合心电时频域特征分析的可电击复律心律自动识别和归类方法【
技术领域：
】[0001]本发明属于医疗电子【
技术领域：
】，具体涉及一种心电信号（ECG)的自动识别和归类方法，特别涉及一种可用于心电监护仪和全自动体外除颤仪（AED)的可电击复律心律的自动识别和归类方法及算法。【
背景技术：
】[0002]美国健康统计中心的流行病学研究结果显示，心血管疾病的死亡总数50%以上是由心脏粹死（SuddenCardiacDeath,S⑶）导致的。临床流行病学研究表明，在过去的五十年中，S⑶的发生率有上升趋势，其中，80%的S⑶是由心室纤维颤动（VentricularFibrillation-VF，简称室颤）或持续性室性心动过速（VentricularTachycardia-VT，简称室速）恶化所导致的。当患者出现室颤这种最恶性的心律失常，在最短的时间内实施电击除颤（ElectricDefibrillation-ED,简称除颤）是临床上唯一可靠并被广泛使用的室颤转复方法。[0003]1986年美国心脏协会（AHA)已确信早期除颤是建立强有力的生存链中最重要的一环，于是设立了"公共场所除颤"特别工作小组，并在1994年、1997年两次召开国际会议，AED被认为是新的革命性技术。国际心脏复苏联络委员会（ILCOR)对公共场所实现心肺复苏急救中进行除颤收集了大量资料，根据世界各地的统计资料形成国际性共识:AED是构建公共场所快速除颤的最有希望的方法；急救存活率高于紧急医疗救护（EMT)(30:23%)，速度也更快（3.3分：8.1分钟）；在社区中值得培训和推广。1997年，美国心脏协会在Circulation期刊上发表了一篇与AED算法性能报道相关的建议《AutomaticExternalDefibrillatorsforPublicAccessDefibrillation:RecommendationsforSpecifyingandReportingArrhythmiaAnalysisAlgorithmPerformance，IncorporatingNewWaveforms，andEnhancingSafety》，美国心脏协会AED特别工作组及其附属的AED安全与功效小组委员会和AED制造商座谈小组一致同意将心律分为以下三大类：[0004](1)可电击复律心律（ShockableRhythms，简称ShR):包括粗颤（CoarseVF)和快速变化的室速（RapidVT);[0005](2)不可电击复律心律（Non-ShockabIeRhythms，简称NShR):包括窦性心律(NSR)、房扑和房颤（AF)、室性逸搏（PVC)、室上性心动过速（SVT)、心脏停搏（Asystole)等；[0006](3)中间心律（IntermediateRhythms):包括细颤（FineVF)和其它的室速(OtherVT)。[0007]目前，AED的相关技术已经比较成熟，国外临床应用也非常广泛。然而，其核心的心律自动识别和分类技术至今仍未公开，国内相关技术的研究亦几乎处于空白状态，由此导致现有产品完全依赖进口，价格居高不下，严重限制了AED-类惠民产品的国内临床应用和推广。在此背景下，我们必须突破国外核心技术的封锁，并研制出对核心技术具有自主知识产权的AED。[0008]在AED使用中，心律的自动识别和正确分类是其核心和关键技术，即当患者发生VF或者持续性VT时，需要即刻对其做出准确地识别和判断，并对患者立即进行电除颤心律转复等治疗措施。这其中，如果可电击复律心律被漏检或误检为不可电击复律心律（NShR)，则将错失及时救助患者生命的机会；相反，如果不可电击复律心律（NShR)被误检为VF或持续性VT，而对病人实施不必要的电击，则可能对人体心脏造成致命性的损害。因此，为了让没有专业急救经验的公众能快速、便捷地操作和使用AED对身边患者实施早期除颤，一套高效、准确的可电击复律心律自动识别和归类方法及算法是其关键。[0009]近年来，随着心脏猝死及电击除颤成功率的不断提高，用于全自动体外除颤仪的可电击复律心律自动识别方法及算法的研究受到了国内外学者的广泛重视。迄今，已有多种算法报道，包括：幅度概率密度函数算法、自相关函数算法、TCI算法、L-Z复杂度算法、快速模板匹配算法和希尔伯特（HILB)算法等。1976年，《MedicalInstrument〉〉曾报道了一篇题为"Considerationsinthedevelopmentoftheautomaticimplantabledefibrillator"的文章，即将心电幅度的分布情况定义为幅度概率密度函数（PDF)，并据此来检测VF。然而，该方法却无法区分VF和VT。随后，《ComputersinCardiologyProceedings〉〉和〈〈Medical&BiologicalEngineering&Computing〉〉进一步干丨J发了"Fibrillationrecognitionusingautocorrelationanalysis，'和"Ventricularfibrillationdetectionbyaregressiontestontheautocorrelationfunction，'等基于周期信号或具有一致性幅值分布特征检测ShR和NShR的文章。但是，却存在当心电信号的形态发生很大的变化时，基于心电图R波识别的算法无法适应于可电击复律心律的识另IJ，由此导致其识别的灵敏度和特异性较低的缺点。1994年和1997年，《IEEETransactionsonBiomedicalEngineering》先后干丨J发了"Multiwaysequentialhypothesistestingfortachyarrhythmiadiscrimination，'和"Detectingventriculartachycardiaandfibrillationbycomplexitymeasure",报道了基于各种变换和复杂度分析的算法，然而这类算法却存在计算复杂、对硬件计算负荷要求高等缺点，无法用之于诸如AED等一类便携式设备。最近，《ComputersinCardiology》（2005年）和《IEEETransactionsonBiomedicalEngineering》（2007年）先后手艮道了"Anewventricularfibrillationdetectionalgorithmforautomatedexternaldefibrillators，'和"DetectingVentricularFibrillationbyTime-DelayMethods"两篇关于HILB的新算法，此类方法及算法是基于相空间重建的算法，虽然特异性有了较大的提高，然而识别灵敏度却较低。综上所述，迄今虽已公开报道了多种可电击复律心律识别和归类方法及算法，然而却仍缺少对识别灵敏度和特异性二者的兼顾考虑，而其计算复杂度亦远远超过便携式AED-类设备的计算负荷。为此，在受制于国外各大厂商对AED核心技术封锁，而同时国内医患对AED的巨大潜在需求的背景下，公开和建立一套兼顾高灵敏度和高特异性、计算简单可靠、具有自主知识产权的可电击复律心律自动识别和归类方法及算法，用于自主研发的心电监护仪和AED是本发明的初衷。【
发明内容】[0010]本发明的目的在于提供一种识别效率高、分类准确、计算简单并满足应用要求的可电击复律心律自动识别和归类方法及算法，其核心结合了反映心电频域特征的带通滤波分析和时域特征的标准化栅条投影标准差分析，用于心电监护仪和全自动体外除颤仪等，以改进现有的需要使用可电击复律心律自动识别和归类方法的仪器设备的性能。[0011]本发明的目的主要由以下技术方案来实现：[0012](1)计算心电信号幅值绝对值的最大值（Max(AbsFS))，根据Max(AbsFS)值判别心脏停搏心律：[0013]如果满足条件：Max(AbsFS)〈停搏阈值，则判定该心律为心脏停搏心律；[0014]如果满足条件：Max(AbsFS)>=停搏阈值，则不是心脏停搏心律，则继续执行后续步骤。[0015](2)根据心电信号的频域特性，计算心电信号经过带通滤波器处理后获得的反映心电信号频域特征的输出信号最大幅度比例值（Pa)、平均幅度比例值（Pb)和平均偏差比例值（Pc)，然后根据Pa、Pb和Pc值判别ShR和NShR:[0016](3)如果满足条件：Pa〈阈值_1，Pb>阈值_2和PaXPb/Pc〈阈值_3，则判定该心律为NShR;如果满足条件：阈值_l〈Pa〈阈值_4,Pb〈阈值_5和PaXPb/Pc〈阈值_3，则判定该心律为NShR;[0017]如果满足条件：Pa>阈值_1和Pb〈阈值_2，则判定该心律为ShR;[0018]如果满足条件：Pb>阈值_6，则判定该心律为ShR。[0019](4)根据心电信号的时域特性，计算标准化栅条投影标准差，根据标准化栅条投影标准差判别NSR和NShR:[0020]如果满足条件：标准化栅条投影标准差>=阈值，则判定该心律为NShR;[0021]如果满足条件：标准化栅条投影标准差〈阈值，则判定该心律为ShR。[0022]上述根据心电信号频域特性，计算经过带通滤波器处理后输出的最大幅度比例值、平均幅度比例值和平均偏差比例值的具体过程如下：[0023]首先，将心电信号通过一个整系数带通滤波器；由于不同心律的分布频带范围不同，从而VF和VT的波形将得到抑制；[0024]其次，对通过上述整系数带通滤波器的输出信号取绝对值，用AbsFS表示；[0025]再次，计算信号绝对值小段时间内的最大值Max(AbsFS)、平均值Mean(AbsFS)和平均偏差（MD);[0026]最后,计算输出信号的Pa、Pb和Pc。最大幅度比例值（Pa):先计算心电信号小段时间内在0.5XMax(AbsFS)到Max(AbsFS)的之间的幅度比例，然后将大段时间内的相加；平均幅度比例值（Pb):先计算心电信号小段时间内在Mean(AbsFS)到Max(AbsFS)的之间的幅度比例，然后将大段时间内的相加；平均偏差比例值（Pc):先计算心电信号小段时间内在Mean(AbsFS)-MD到Mean(AbsFS)+MD的之间的幅度比例，然后将大段时间内的相应幅度比例相加。[0027]上述根据心电信号时域特性，计算标准化栅条投影标准差的具体过程如下：[0028]首先，将采集到的一段时间内的心电信号按相同时间间隔分成小段，每一小段称为一个栅条，每一个时间间隔称为栅条宽度；[0029]其次，计算心电图在每个栅条内覆盖的幅度范围，就是栅条在纵轴上的投影；[0030]再次，计算所有栅条投影的标准差；[0031]最后，对栅条投影的标准差进行标准化，得到标准化栅条投影标准差。[0032]采用以上技术方案不仅提高了可电击复律心律识别的灵敏度和不可电击复律心律识别的特异性、满足算法性能的应用要求，而且也简化了算法的计算复杂度。本发明可应用于心电监护仪和全自动体外除颤仪等仪器设备。【专利附图】【附图说明】[0033]图1是本发明的方法及算法主体流程图。[0034]图2是本发明主体流程图中"S1心电信号预处理"步骤的流程图。[0035]图3是本发明主体流程图中"S3计算带通滤波器输出的最大幅度比例值（Pa)、平均幅度比例值（Pb)和平均偏差比例值（Pc)"步骤的流程图。[0036]图4是本发明主体流程图中"S8计算标准化栅条投影标准差"步骤的流程图。[0037]图5是本发明中窦性心律（NSR)经过各个阶段处理的图形。[0038]图6是本发明中室颤心律（VF)经过各个阶段处理的图形。【具体实施方式】[0039]本实施例是在个人计算机（PC)和Matlab软件平台上的一种可能的实现，如图1?6所示，其具体实施步骤如下：[0040]1?对采集到的心电信号进行预处理：[0041](1)使用一个截止频率为IHz的高通滤波器，抑制基线漂移；[0042](2)使用一个截止频率为30Hz的巴特沃斯低通滤波器，滤除工频干扰和肌电噪声等；[0043](3)使用一个简单的滑动平均滤波器，进一步滤除无关的高频干扰，提高滤波效果。[0044]2.对心电信号进行心脏停搏心律的识别：如果满足条件：Max(AbsFS)〈150iiV，则判定该心律是心脏停搏心律；如果满足条件=Max(AbsFS)>=150yV，则判定该心律不是心脏停搏心律，继续执行后续的步骤。[0045]3.根据心电信号频域特性，计算心电信号经过带通滤波器处理后的反映心电信号频域特性的输出信号最大幅度比例值、平均幅度比例值和平均偏差比例值：[0046](1)将心电信号通过一个整系数带通滤波器，带通滤波器的数学公式如下所示。由于不同心律的分布频带范围不同，从而VF和VT的波形得到抑制；)/8[0047](2)对通过上述整系数带通滤波器的输出信号取绝对值，用AbsFS表示；[0048](3)计算信号绝对值每Is内的最大值Max(AbsFS)、平均值Mean(AbsFS)和平均方差（MD);[0049](4)计算输出信号的Pa、Pb和Pc。最大幅度比例值（Pa):先计算Is内在0.5XMax(AbsFS)到Max(AbsFS)的之间的幅度比例数，然后将10内的相加；平均幅度比例值（Pb):先计算Is内在Mean(AbsFS)到Max(AbsFS)的之间的幅度比例数，然后将10内的相加；平均偏差比例值（Pc):先计算Is内在Mean(AbsFS)-MD到Mean(AbsFS)+MD的之间的幅度比例数，然后将10内的相加。[0050]4.根据Pa、Pb和Pc的值来判别ShR和NShR:[0051]如果满足条件：Pa〈250,Pb>950和PaXPb/Pc〈210,则判定该心律为NShR;[0052]如果满足条件：250〈Pa〈400,Pb〈600和PaXPb/Pc〈210,则判定该心律为NShR;[0053]如果满足条件：Pa>250和Pb〈950,则判定该心律为ShR;[0054]如果满足条件：Pb>1100,则判定该心律为ShR。[0055]5.根据心电信号的时域特性，计算其标准化栅条投影标准差：[0056](1)将IOs内的心电信号按相同时间间隔分成小段，每一小段称为一个栅条，每一小段的时间间隔称为栅条宽度；[0057](2)计算心电图在每个栅条内覆盖的幅度范围，就是栅条在纵轴上的投影；[0058](3)计算所有的栅条投影的标准差；[0059](4)对栅条投影的标准差进行标准化，得到标准化栅条投影标准差。[0060]6.根据标准化栅条投影标准差来判别ShR和NShR:[0061]如果满足条件：标准化栅条投影标准差>=阈值，则判定该心律为NShR;[0062]如果满足条件：标准化栅条投影标准差〈阈值，则判定该心律为ShR。[0063]本实施例所使用的软硬件配置如下：[0064]硬件：戴尔台式电脑，Inter(R)Core(TM)i5处理器和4G内存；[0065]软件：MATLABR2010b。[0066]本实施例的训练数据库和测试数据库分别为My_DB02和My_DB03,包括NSR、VT、VF及其他心律。两数据库中所有的心电数据片段均来自于美国麻省理工大学心律失常数据库(MITDB)、美国克莱登大学室性心律失常数据库（CUDB)、美国麻省理工大学恶性室性心律失常数据库（VFDB)和美国心脏协会数据库（AHADB)四个标准的心电数据库。其中，每个心电数据片段按长度为l〇s，相邻两个片段间隔Is的方式连续截取获得。My_DB02和My_DB03数据库详情如下表：【权利要求】1.一种结合心电时频域特征分析的可电击复律心律自动识别和归类方法，其特征是结合心电频域特征的带通滤波分析和心电时域特征的标准化栅条投影标准差分析，适用于疾病监测和治疗的心电监护仪和全自动体外除颤仪（AED)，主要步骤包括：51.对采集到的心电信号进行预处理；52.对心电信号进行心脏停搏心律的识别，若是心脏停搏心律，则判为不可电击复律心律；否则，则继续执行后续步骤S3;53.将心电信号通过整系数带通滤波器处理，获得反映心电频域特征的输出信号最大幅度比例值（Pa)、平均幅度比例值（Pb)和平均偏差比例值（Pc)及各阈值；54.根据带通滤波器输出信号的Pa、Pb和Pc值，判别可电击复律心律和不可电击复律心律，若满足条件：Pa〈阈值_1，Pb>阈值_2和PaXPb/Pc〈阈值_3，则判定该心律为不可电击复律心律；否则，继续执行后续步骤S5;S5?若满足条件：阈值_1〈Pa〈阈值_4，Pb〈阈值_5和PaXPb/Pc〈阈值_3，则判定该心律为不可电击复律心律；否则，继续执行后续步骤S6;56.若满足条件：Pa>阈值_1和Pb〈阈值_2，则判定该心律为可电击复律心律；否则，继续执行后续步骤S7;57.若满足条件：Pb>阈值_6，则判定该心律为可电击复律心律；否则，继续执行后续步骤S8;58.计算反映心电时域特征的标准化栅条投影标准差；59.根据标准化栅条投影标准差判别可电击复律心律和不可电击复律心律，若标准化栅条投影标准差〈阈值，则判定该心律为可电击复律心律；若标准化栅条投影标准差>=阈值，则判定该心律为不可电击复律心律。2.如权利1所述的一种结合心电时频域特征分析的可电击复律心律自动识别和归类方法，其特征在于，所述的带通滤波方法中Pa、Pb和Pc的计算步骤包括：首先，将心电信号通过一个整系数带通滤波器，由于不同心律的分布频带范围不同，从而VF和VT的波形将得到抑制；其次，对通过上述整系数带通滤波器的输出信号取绝对值，用AbsFS表示；最后，计算信号绝对值小段时间内的最大值Max(AbsFS)、平均值Mean(AbsFS)和平均方差（MD)，其中Pa、Pb和Pc定义如下：最大幅度比例值（Pa):先计算心电信号小段时间内在0.5XMax(AbsFS)到Max(AbsFS)的之间的幅度比例，然后将大段时间内的相加；平均幅度比例值（Pb):先计算心电信号小段时间内在Mean(AbsFS)到Max(AbsFS)的之间的幅度比例，然后将大段时间内的相加；平均偏差比例值（Pc):先计算心电信号小段时间内在Mean(AbsFS)-MD到Mean(AbsFS)+MD的之间的幅度比例，然后将大段时间内的相加。3.如权利1所述的一种结合心电时频域特征分析的可电击复律心律自动识别和归类方法，其特征在于，所述的标准化栅条投影标准差的计算步骤包括：首先，将采集到的一段时间内的心电信号按相同的时间间隔分成小段，每一小段称为一个栅条，每小段的时间间隔称为栅条宽度；其次，计算心电图在每个栅条内覆盖的幅度范围，就是栅条在纵轴上的投影；再次，计算所有栅条投影的标准差；最后，对栅条投影的标准差进行标准化处理，得到标准化栅条投影标准差。4.如权利1所述的一种结合心电时频域特征分析的可电击复律心律自动识别和归类方法，其特征在于，所述的心电信号预处理，具体步骤包括：首先，使用一个截止频率小于或等于1Hz的高通滤波器，抑制基线漂移；其次，使用一个截止频率大于或等于30Hz的低通滤波器，滤除工频干扰和肌电噪声等；最后，使用一个滑动平均滤波器，进一步滤除与心电信号无关的高频干扰，提高滤波效果。5.如权利1所述的一种结合心电时频域特征分析的可电击复律心律自动识别和归类方法，其特征在于，所述的阈值_1?阈值_6:对所采集到的一段时间内的大量心电数据进行统计分析，由数据统计分析软件计算确定。6.如权利1所述的一种结合心电时频域特征分析的可电击复律心律自动识别和归类方法，其特征在于，步骤S2中当心电信号幅度的绝对值小于停搏阈值则判别为心脏停搏心律。【文档编号】A61B5/0452GK104382590SQ201410764262【公开日】2015年3月4日申请日期:2014年12月11日优先权日:2014年12月11日【发明者】赖大坤,张飞申请人:赖大坤