[发明专利]一种CRP伪影抑制方法、系统、电子设备和存储介质

说明书

本发明公开一种CRP伪影抑制方法、系统、电子设备和存储介质，包括以下步骤：S1：获取待处理的包含CPR伪影的原始ECG数据；S2：对原始ECG数据进行预处理，将ECG信号和CPR伪影进行第一次分离，得到不包含ECG信号的CPR伪影信号；S3：计算S2中CPR伪影信号的功率谱密度；S4：根据功率谱密度PSD找出CPR伪影信号的分量频点，再根据分量频点对CPR伪影信号进行重构，得到CPR干扰信号；S5：采用卡尔曼自适应滤波器对CPR干扰信号进行递归求解得到CPR伪影估计量；S6：从原始ECG数据去除CPR伪影估计量，得到无CPR伪影的ECG数据。通过该方法，可以从原始ECG数据中信号中去除CPR伪影，提高除颤判断类型的准确性。

技术领域

本发明涉及心电信号处理技术领域，特别涉及一种CRP伪影抑制方法、系统、电子设备和存储介质。

背景技术

当受到外部严重撞击或自身疾病出现变化时，心肺有可能出现停止跳动，将对人体造成损害。在现代急救技术中，一般是采用按压的方式对患者进行心肺复苏，同时除颤仪对患者的心电信号进行分析显示，以供工作人员直观查看，决定是否进行除颤。

但在胸外持续按压期间，CRP(cardiopulmonary resuscitation，心肺复苏)伪影信号将对原始心电信号ECG产生极大干扰，会影响除颤仪对心电节律的分析，降低除颤的分析结果准确性，可能导致无法正确区分除颤/不可除颤类别。

发明内容

针对现有技术中CRP伪影对心电信号干扰的技术问题，本发明提出一种CRP伪影抑制方法、系统、电子设备和存储介质，通过带通滤波器对原始心电信号进行初步分离，然后利用主频分量的频点进行相位CPR干扰信号重构，再用卡尔曼自适应滤波器进行递归估计得到重建的CPR干扰信号，最后从原始心电信号中去除CPR干扰信号。

为了实现上述目的，本发明提供以下技术方案：

一种CRP伪影抑制方法，包括以下步骤：

S1：获取待处理的包含CPR伪影的原始ECG数据；

S2：对原始ECG数据进行预处理，将ECG信号和CPR伪影进行第一次分离，得到不包含ECG信号的CPR伪影信号；

S3：计算S2中CPR伪影信号的功率谱密度；

S4：根据功率谱密度PSD找出CPR伪影信号的分量频点，再根据分量频点对CPR伪影信号进行重构，得到CPR干扰信号；

S5：采用卡尔曼自适应滤波器对CPR干扰信号进行递归求解得到CPR伪影估计量；

S6：从原始ECG数据去除CPR伪影估计量，得到无CPR伪影的ECG数据。

优选地，所述S1中，原始ECG数据的类型包括至少两个波形；每个波形包括至少100个样本，每个样本至少有1000个采样点。

优选地，所述S2中，预处理方法为：

采用巴特沃斯一阶滤波器设计带通频率为[1Hz,3Hz]的带通滤波器，对原始ECG数据进行低通滤波去除其中的ECG信号，实现ECG信号和CPR伪影的第一次分离。

优选地，所述S3包括：

步骤S31：对CPR伪影信号进行自相关计算：

公式(1)中，R_x(τ)表示对x(t)进行自相关计算的自相关结果；t表示时间；τ表示时延；x(t)表示经过带通滤波器后时间t的时域信号；x(t+τ)表示经过带通滤波后时间t+τ的时域信号；

步骤S32：对自相关结果R_x(τ)进行傅里叶变换：