[发明专利]一种基于PQRST波位置识别的RonT信号检测方法

权利要求书

1.一种基于PQRST波位置识别的RonT信号检测方法，其特征在于，包括如下步骤：

a)建立RonT心电数据ALL_ECG，

，

式中为第i个心电信号，i∈{1,...,n}；

b)对RonT心电数据ALL_ECG进行滤波，消除高频干扰；

c)对滤波后的RonT心电数据ALL_ECG去除基线漂移；

d)将去除基线漂移后的RonT心电数据ALL_ECG进行归一化处理；

e)对的RonT心电数据ALL_ECG中的12导联心电数据振幅进行绝对值处理，根据振幅的绝对值将12导联合并成一个导联x(nT)，n∈{1,...,length}，length为RonT心电数据ALL_ECG的信号长度，T为采样率；

f)利用Pan-Tompkins算法识别合并后的导联信号x(nT)中的R波峰值；

g)根据R波峰值检测QRS复合波起止点；

h)根据R波峰值检测P波起止点；

i)根据R波峰值检测T波起止点；

j)根据PQRST波形的起止位置对RonT心电数据ALL_ECG划分为若干个心搏，每个心搏包含一个PQRST波形，通过PQRST波形的位置信息以及步骤d)中取出基线漂移后的RonT心电数据ALL_ECG进行RonT信号检测；

步骤j)包括如下步骤：

j-1)设置阈值Flag_RonT，并将阈值Flag_RonT赋值为0；

j-2)判断当前心搏中是否具有P波，如果没有P波则转至执行步骤j-4)；

j-3)如果P波存在则寻找下一个心搏并转至执行步骤j-2)；

j-4)对当前心搏所在的导联的所有心搏R波幅值求平均，得到均值R_avg，如果QRS波宽度大于阈值threshold_QRS且当前心搏R波幅值大于1.1倍的R_avg，则转至执行步骤j-5)，如果QRS波宽度小于等于阈值threshold_QRS则转至执行步骤j-2)；

j-5)计算QT间期，如果QT间期的长度小于设定的阈值threshold_QT，则阈值Flag_RonT赋值为1，如果QT间期的长度大于等于设定的阈值threshold_QT则转至步骤j-2)；

j-6)如果阈值Flag_RonT等于1则表示存在RonT信号，并转至执行步骤j-7)，阈值Flag_RonT等于0则表示不存在RonT信号，并转至执行步骤j-2)；

j-7)输出RonT识别结果。

2.根据权利要求1所述的基于PQRST波位置识别的RonT信号检测方法，其特征在于：步骤b)中采用窗宽为150、截止频率为40Hz的FIR数字低通滤波器对RonT心电数据ALL_ECG进行滤波。

3.根据权利要求1所述的基于PQRST波位置识别的RonT信号检测方法，其特征在于：步骤c)中采用中值滤波器将滤波后的RonT心电数据ALL_ECG减采样到50Hz后取窗宽为560ms的长窗，对长窗中70％的中间值求平均值得到当前基线值，将当前基线值重采样到原频率后作为基线被移除。

4.根据权利要求1所述的基于PQRST波位置识别的RonT信号检测方法，其特征在于：步骤d)中通过公式计算得到导联名称为m的归一化后的心电数据为获得导联名称为m的全部心电数据，为获得导联名称为m的全部心电数据的最小值，为获得导联名称为m的全部心电数据的最大值。

5.根据权利要求1所述的基于PQRST波位置识别的RonT信号检测方法，其特征在于，步骤g)中检测QRS复合波起止点的方法为：

g-1)以R波为基准，在R波前0.1RR间期数据段内，利用峰值位置确定该段数据内的极小值点作为Q波波峰点；

g-2)以R波为基准，在R波后0.1RR间期数据段内，利用峰值位置确定该段数据内的极小值点作为S波波峰点；

g-3)将Q波波峰点左侧15ms时间窗口内最靠近基线的点作为Q波起点，将S波波峰点右侧15ms时间窗口内最靠近基线的点作为S波终点。