预判除颤时机的方法与流程

本发明属于心电信号处理的技术领域，具体涉及一种其于对心电信号进行时间-频率的能量分析来预判除颤时机的方法。

背景技术：

心脏骤停是造成现代人类非正常死亡的一个重要因素，而60％以上的心脏骤停是由于心室颤动引起的(name+year,)。根据大量临床数据及相关资料我们可得知患者发生室颤后抢救时间窗为10min，最佳抢救时间为最初的3—5min，而每延迟一分钟做cpr和除颤患者的生存率将以7％—10&递减。到目前为止，及时除颤是唯一有效的治疗手段，除此之外对于终止室颤还没有特别有效的治疗方法。

早期研究证实，高能量的电除颤会对心肌功能产生不利的影响，因此对于除颤时机的掌握就显得更为重要，如果盲目的对病人除颤不仅不会对病人的急救产生有效的帮助，反而会加重心肌功能的损伤，造成除颤失败。

早期医护人员通过ecg(心电图)的幅值来预测除颤时机，幅值越大除颤成功率也就越高，但是存在个体差异大的问题(1996，hans)。也有采用冠脉灌注压(cpp)的预测除颤时机的，然而cpr对于院外的心脏骤停患者以及院内没有预期的心脏骤停患者的治疗并不适用。到目前为止，对于除颤时机能够较为有效的预测，提高除颤的效率，和指导cpr的操作的只有heitor于2000年提出的amsa得到较为广泛的认同。amsa被定义为除颤前5秒，频率范围4-48hz的ecg信号的幅值谱曲线下的面积。

在室颤阶段，心肌细胞内的高能磷酸盐(atp)浓度影响着心肌细胞的新陈代谢，并且与心肌细胞的损伤程度成负相关，心肌细胞内atp浓度越高，对于除颤成功的结果越有帮助(robert，1991)。ecg信号能够反映心肌细胞内的atp浓度水平(david,2009)。

amsa能够从侧面反映心肌细胞的能量(david,2009)，然而，amsa并不能直接的代表和说明心肌细胞的能量的高低。

基于上述原因，amsa对于除颤时机的判断并不完善。

上述论述内容目的在于向读者介绍可能与下面将被描述和/或主张的本发明的各个方面相关的技术的各个方面，相信该论述内容有助于为读者提供背景信息，以有利于更好地理解本发明的各个方面，因此，应了解是以这个角度来阅读这些论述，而不是承认现有技术。

技术实现要素：

本发明的目的在于避免现有技术中的不足而提供一种预判除颤时机的方法，其通过对心电信号进行能量分析，能更精准判断心肺复苏过程中对患者除颤的最佳时机，为施救者提供一套是否除颤且何时进行除颤的判断标准。

本发明的目的通过以下技术方案实现：

提供一种预判除颤时机的方法，其特征在于，包括以下步骤：

s1.对采集获取的ecg信号作带通滤波；

s2.对步骤s1中处理后的滤波信号作快速傅里叶变换，得到滤波信号对应的心电频谱x(ω)，并列出相应的频率-幅值对(ωi,ai)，ai为心电频谱中不同频率点ωi所对应的幅值；

s3.根据步骤s2处理后的滤波信号的幅值ai，作以下计算power＝∑ai²；

s4.根据步骤s3得到的power值，计算预判除颤成功率p，设定阈值m与p值比较判断是否除颤。

作为进一步的改进，所述步骤s1中的带通滤波频段选定为12～30hz。

作为进一步的改进，所述步骤s1中，对采集获取的ecg信号放大500-1000倍。

作为进一步的改进，所述步骤s2中，ai的计算公式为：其中n为频率点ωi个数。

作为进一步的改进，所述步骤s4中，预判除颤成功率p的计算公式为:

power为步骤s3求得的滤波幅值平方和的值，a为0.018，b为-1.30。

作为进一步的改进，所述步骤s4中，根据步骤s3计算得的power值，阈值m取0.8，依据计算得的预判除颤成功率p值，若p>m，则对患者进行除颤；若p<m，则继续进行心肺复苏，若p值在区间[0.2,0.8]内，则每隔xs对患者做一次除颤。

本发明提供的预判除颤时机的方法，包括以下步骤：s1.对采集获取的ecg信号作带通滤波；s2.对步骤s1中处理后的滤波信号作快速傅里叶变换，得到滤波信号对应的心电频谱x(ω)，并列出相应的频率-幅值对(ωi,ai)，ai为心电频谱中不同频率点ωi所对应的幅值；s3.根据步骤s2处理后的滤波信号的幅值ai，作以下计算power＝∑ai²；s4.根据步骤s3得到的power值，计算预判除颤成功率p，设定阈值m与p值比较判断是否除颤。本发明通过对心电信号进行能量分析，能更精准判断心肺复苏过程中对患者除颤的最佳时机，为施救者提供一套是否除颤且何时进行除颤的判断标准。

附图说明

利用附图对本发明作进一步说明，但附图中的实施例不构成对本发明的任何限制，对于本领域的普通技术人员，在不付出创造性劳动的前提下，还可以根据以下附图获得其它的附图。

图1为本发明方法实施的流程示意图。

图2为心电监护仪工作系统流程示意图。

图3为心电监护仪led指示灯示意图。

图4为本发明方法与amsa方法初次除颤成功与失败t检验值的对比图。

图5为power方法与amsa法的roc曲线。

图6为power与amsa方法roc曲线指标数据比较图。

具体实施方式

为了使本领域的技术人员更好地理解本发明的技术方案，下面结合附图和具体实施例对本发明作进一步详细的描述，需要说明的是，在不冲突的情况下，本申请的实施例及实施例中的特征可以相互组合。

实施例一：

患者心脏骤停后应对其及时进行心肺复苏和除颤，若已判断患者出现室颤能及时对其进行电击除颤正确判能大大提高患者的救活率，而在心肺复苏过程中能判断出最佳除颤时机是患者成活的关键。本实施例以已判断患者出现室颤时除颤仪应用于人体为例进行说明。

图1为本发明中提出的除颤时机预判方法的流程图，包括以下步骤：

s1.对采集获取的ecg信号作带通滤波；

s2.对步骤s1中处理后的滤波信号作快速傅里叶变换，得到滤波信号对应的心电频谱x(ω)，并列出相应的频率-幅值对(ωi,ai)，ai为心电频谱中不同频率点ωi所对应的幅值；

s3.根据步骤s2处理后的滤波信号的幅值ai，作以下计算power＝∑ai²；

s4.根据步骤s3得到的power值，计算预判除颤成功率p，设定阈值m与p值比较判断是否除颤。

本实施例推荐使用12导联信号进行采集计算，步骤s1中使用的是4—48hz的带通滤波器，4—48hz的带通滤波器对工频干扰及基线漂移进行滤除。根据临床数据，所述步骤s1中的带通滤波段选定为12～30hz。

步骤s2中对步骤s1的滤波结果进行快速傅里叶变换，得到滤波信号的频谱x(ω)，ai为心电频谱中不同频率ωi所对应的幅值，计算方式如下：其中n为频率点ωi个数。

步骤s3中，基于帕斯瓦尔定理和临床数据的验证，本发明提出以下计算公式power＝∑ai²。

步骤s4中预判除颤成功率计算公式为，

其中a为0.018，power为上述步骤s3求得的滤波幅值平方和的值，b为-1.30。

步骤s4根据步骤s3计算得的power值，阈值m取0.8,依据计算得的预判除颤成功率p值，若p>m，则对患者进行除颤；若p<m，则继续进行心肺复苏，若p值在区间[0.2,0.8]内，则每隔xs(x秒)对患者做一次除颤,x为3至10中的整数。

结合图2所示，心电监护仪工作系统包括：心电采集模块，算法模块(图中未绘出)，电除颤模块，交互模块。心电采集模块由放大电路、4—48hz带通滤波器、数/模转换电路组成，放大电路把采集的患者心电信号放大500-1000倍，4—48hz的带通滤波器对工频干扰及基线漂移进行滤除。电除颤模块由充电电路和放电电路组成，算法模块由(心电波形分析算法及除颤时机判断算法)组成，除颤指令由算法模块发出，经过通信链路传输到电除颤模块，当心电分析确定需要进行除颤时，充电电路会自动充电为除颤做好准备。充电电路工作过程中，继续分析患者心电信号。。交互模块由显示模块和声音模块组成。本实施例基于的除颤仪中的显示模块设定为图3所示的一组led指示灯，所述的一组led指示灯根据心电分析的不同结果有不同的显示，相应地语音提示器会作出相关语音提示。声音模块为一语音提示器。

心电监护仪工作时，患者的心电信号经过放大电路，并经过4—48hz的带通滤波器的处理，得到去除了高频干扰的信号。此信号再经过算法模块分析，判断是否需要除颤以及除颤成功率p。

结合图3所示，当p>m时心电分析确定需要对患者进行除颤，充电电路开始自动充电为除颤做准备，充电进程有语音提示，“分析”led指示灯闪绿光，同时仍对心电信号继续分析，若分析仍建议进行除颤，语音提示器作出“按下电击键”的提示，当操作者按下“电击”键时，除颤器放电电路工作，进行电击除颤。

当p<m时心电分析确定无需电击，除颤仪就不充电，“分析”led指示灯闪红光，操作者对患者继续做心肺复苏，除颤仪算法模块继续分析，

当p值在区间[0.2,0.8]时，“分析”led指示灯闪黄光，操作者按指示每隔xs对患者做一次电击除颤，xs之间继续做心肺复苏，除颤仪继续对患者的心电信号分析，每一次电击除颤结束后语音提示器都会作出相关语音提示。

实施例二：

本发明还适用于动物患者，本实施例以对某种动物的临床实验为例进行说明，操作者应明白本发明提出的方法并不局限于本实施例的动物。

当已判断出患者出现心室颤可实施本方法。结合图1本方法的实施流程图，包括以下步骤：

s1.选取x实验动物分别获取心电信号，对采集模块中获取的心电信号作带通滤波，所述步骤s1中的带通滤波段选定为12～30hz，同时进行对每一个样本进行首次除颤；

s2.对步骤s1中采集的除颤前心电信号作快速傅里叶变换，得到每段滤波信号对应的频谱，并列出相应的频率-幅值对(ωi,ai)；

s3.根据步骤s3处理后的滤波信号的幅度谱，分别计算每个样本的power值和amsa值；其中power＝∑ai²，amsa＝∑ai·fi。

s4.采用t检验分析步骤s1中除颤成功和除颤失败的power值和amsa值，同时对power方法和amsa法进行线性回归，分析相关程度。

s5.绘制power法和amsa方法的roc曲线，计算它们的auc。

步骤s1中本实施例使用ⅱ导联信号进行采集计算，分别截取初次除颤前长度为1秒的波形并标记该次除颤是否成功，使用4—48hz的带通滤波器进行滤波去除胸外按压所带来的低频干扰及工频辐射、肌电等所引起的高频干扰。

根据临床数据，

步骤s2中对步骤s1的滤波结果进行快速傅里叶变换，得到滤波信号的频谱x(ω)，ai为心电频谱中不同频率ωi所对应的幅值，计算方式如下：其中n为频率点ωi个数。

步骤s3中，基于帕斯瓦尔定理和临床数据的验证，本发明提出以下计算公式power＝∑ai²。本实施例给出amsa的计算公式：amsa＝σai·fi,fi为不同的心电频率。

结合图4，步骤s4中power法和amsa方法的统计数据以“均值±标准差”表示。结合图5，步骤s5绘制了power法和amsa方法的roc曲线，分别计算两者的auc，相关数据比较如图6所示。

上面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用其他不同于在此描述的其他方式来实施，因此，不能理解为对本发明保护范围的限制。

总之，本发明虽然列举了上述优选实施方式，但是应该说明，虽然本领域的技术人员可以进行各种变化和改型，除非这样的变化和改型偏离了本发明的范围，否则都应该包括在本发明的保护范围内。