《心电信号的QRS波群检测：差分运算方法.doc》由会员分享，可在线阅读，更多相关《心电信号的QRS波群检测：差分运算方法.doc（13页珍藏版）》请在金锄头文库上搜索。
1、心电信号的QRS波群检测：差分运算方法叶云迟 中国台湾中坜市320号国立中央大学电子工程系 中国台湾中坜市320号清云大学电子工程系王文君 中国台湾中坜市320号国立中央大学电子工程系 中国台湾台北路106号国立台北科技大学电子工程系 版权 2008 INIST，法国国家科学研究中心，保留所有版权。关键词：心电信号 麻省理工学院心律失常数据库 QRS波群检测 差分运算方法摘要：本文用一种简单可靠的方法即差分运算方法（DOM）来检测心电信号的QRS波群。 The proposed DOM includes two stages.DOM包括两个阶段： The first stage is to f
2、ind the point R by applying the difference equation operation to an ECG signal.第一阶段是对一心电信号通过差分运算找到R点。 The second stage looks for the points Q and S based on the point R to find the QRS complex.第二阶段根据R点确定Q点和S点以确定QRS波群。 在From the QRS complex, the T wave and P wave can be obtained by the existing metho
3、ds.在QRS波群，可以通过现有的方法找到T波和P波。 Some records (QRS complex and T and P waves) of ECG signals in MIT-BIH arrhythmia database is tested to show the DOM has a much more precise detection rate and faster speed than other methods.对麻省理工学院心律失常数据库中的心电信号（ QRS波和T和P波）的记录的测试显示的差分运算方法（DOM）比其他方法有一个更精确的检测率和更快的速度。 正文： 1
4、、引言：心电图(ECG)信号是心肌的电活动在体表的表现,它几乎是一个周期信号。众所周知，心电图(ECG)信号包含了许多关于心脏疾病的信息。一般的,心电图的一个周期波形中，每个波峰分别用字母R - P,Q,S和T标识。如图1所示，R峰总是最高。诊断基于这些基准点所对应的时间和形态特征。因此,心电图检测对医生做出正确的临床诊断1是非常重要的。对心电图(ECG)信号的QRS波群检测已经有30年的历史了 。根据医学检测，心电图(ECG)信号的最重要的信息都集中在P波，QRS波群和T波上面。这些数据包括PR间间隔，QRS间间隔, QT间间隔,ST间间隔，PR段和ST段的位置和大小(见图1)。事实上，对Q
5、RS波的检测因为某些原因是非常困难的，如心电信号的噪声，电磁干扰，T波幅度非常接近于QRS波群等。为了得到正确的信息，以上的问题必须被克服。图1 心电信号波形: (1) P波; (2)QRS波群; (3)T波; (4)PR间隔; (5)QRS间隔; (6) QT间隔; (7)ST间隔;(8)PR部分; (9)ST部分; (10)RR间隔(或拍);(11)心跳周期(包括 P wave, QRS 波群, 和 T 波)。这里有一些研究心电信号的QRS波群检测的方法。例如,帕和汤普金斯提出的一种算法(所谓PT法)表明他们利用一个特殊的数字带通滤波器(参考)来分析QRS波群的位置从而减少错误的心电信号检
6、测。一种用小波变换的方法来检测QRS波群的方法，通过小波换算，把QRS波群从高P波、T波，噪音和基线漂移中检测出来。此外，也可使用一种特殊的数字滤波器来及时的检测心电信号并将其分类。或者利用频域提出的“模板匹配方法”来找到心电信号的节拍。基于对一阶导数的估计，提出了一种通用的找出QRS点的算法。正如我们所知，要达到上述几种方法的目标，需要复杂数学运算。 本文提出了一种简单快速的方法来检测QRS波群，即差分运算方法（DOM）。The proposed DOM includes two stages.DOM包括两个阶段： The first stage is to find the point R
7、 by applying the difference equation operation to an ECG signal.第一阶段是对一心电信号通过差分运算找到R点。 The second stage looks for the points Q and S based on the point R to find the QRS complex.第二阶段根据R点确定Q点和S点以确定QRS波群。大多数P和T波检测技术要求QRS作为参考点。这种差分运算方法（DOM）能在不需要复杂的数学运算如傅里叶变换等的情况下很容易的检测出QRS波群。 本文接下来的内容如下：第二部分讲述差分方法的程序；第
8、三部分用麻省理工学院心律失常数据库的数据来评价这种差分运算方法（DOM）的有效性。第四部分为本文的结束部分。 2、差分操作方法 这部分提出一种用差分运算方法（DOM）来检测QRS波群的方案。只要QRS波群被检测出，则我们就能利用现有的文献将P波和T波检测出来。差分运算方法（DOM）包括两个过程，一个是差分运算过程（DOP），另一个是波形检测过程。21 差分运算过程差分运算过程如下：1-1下载一原始的心电信号X。该心电信号可以来源于任何来源。例如麻省理工学院心律失常数据库，如图2所示。图 2 (a) x, 麻省理工学院心律失常数据库中的第103号心电信号记录； (b) xd, 经过差分运算后的心
9、电信号即差分信号；(c) xdf, 经过低通滤波器后的差分信号; (d)阈值的设定； (e) xdf, 通过阈值后的信号。图 3 心电信号的噪音: (a)60Hz电线干扰；(b)基线漂移；(c)肌电图(d)运动伪影。1-2用一数字滤波器去除心电信号的噪音。在不同情况下有多种类型的噪音，如：1）60HZ电线干扰。尽管心电信号测量的接地线接地良好，但仍然会产生干扰，（如图3a）。60HZ的频谱伴随着谐波。60HZ滤波器可以用来抗干扰。2）基线漂移。它主要是由患者的呼吸、坏的电极和电极位置不恰当而引起的，（如图3b）。它的频率范围通常低于0.5HZ。截止频率为0.5HZ 的高通滤波器可消除由基线漂移
10、引起的干扰。 3)肌电图。它是由肌电活动产生。在心电信号中，肌电信号的干扰的产生是由于肌电信号快速波动的改变的速度快于心电图检查波(见图3c)。它的频谱范围为直流到10K HZ。一单位方波（最大宽度为0.07S）的形态学滤波器可用来消除肌电信号的干扰。4）运动伪影。它产生于电极相对于病人皮肤的运动。在心电图中它能产生较大幅度的信号 (见图3d),它的持续时间约为100 - 500ms。可用自适应滤波器来消除运动伪影造成的干扰。1-3 由（1）式获得差分信号xd 。xd(n) = x(n) x(n 1) (1)。x(n)表示在n时刻的输入信号，xd(n)表示在n时刻的差分输出信号。差分信号可能包
11、含小幅度但高频率变化波形（如图2b）。图 5 麻省理工学院心律失常数据库中的第103号心电信号记录中的0.864 s的RR间隔(311 个采样点)，其中，持续时间 = (R2 R1)(1/采样频率) = (577266)(1/360 Hz)。1-4 让xd通过任一截止频率大约为100赫兹的低通滤波器来消除这个小幅度但高频率变化的波形(见图2c)。滤波信号用xdf表示。一般的，一个典型的心电信号的频域变化范围在0.5HZ至100HZ。1-5 用从14中得到的滤波信号xdf通过（2）式得到最终的信号（如图2e）。（2） T1 =2MVp，T2 =2MVn，(见图2d)。MVp和MVn分别表示在麻省
12、理工学院心律失常数据库的记录中所有正负波幅的平均值。根据医学定义,每个R-R波(也就是说,每个节拍)间隔大约为0.4 1.2秒(即在360HZ的抽样频率有144 432个采样点)(见图1)。例如，在麻省理工学院心律失常数据库中的第103条记录的R-R间隔的持续时间大约为0.864 s(即311采样点)(见图5)。在30分钟的记录时间中采样点的总数是2084311（大约2084拍，见表3）。在这个例子中, MVp和MVn的均值分别是0.48 mV与0.47mV。确T1和T2的阈值是由多次实验确定的。设置T1和T2的阈值的目的是为了找到正确的正负级的极值点。但是，如果这两个极值点的值设置的太小，则
13、在心电信号的R-R间隔之间将会有很多非零点，所以正负极值点不容易找到。另一方面，如果这两个极值点的值设置的太大，则在心电信号中可能去除许多R波的波峰，所以不能找到正确的正负极值点。我们的算法取决于心电图中电极的位置。例如,假定R峰总是正极、Q和S峰始终是负极的(根据医学定义,见图5)。如果电极位置与图5中的不相同，就用其他值代替原T1、T2的值（如图2d）。 图 6 选择正确的极值点 (A点是间隔i1中的极值点, B点是间隔i中的极值点, B点到C 点检、间的距离叫RR 间隔。)。图 7 定义搜索间隔1和搜索间隔2。2.2波形检测我们获得滤波信号xdf后，下一步就是QRS波群检测。对一个原始心
14、电信号x，如果分别找到了R点Q点S点，那么就获得了QRS群。2.2.1 寻找R点2-1 将差分信号分成两部分：和 ，分别表示的正极和负极部分，(如图4a)。图 8 确定Q点和S点位置的例子. (a)因为 MVqq (0.22mV)Vq1 (0.1mV)+ Tv (0.18 mV)，所以Q点位于Q2点处。S1点和S2点的位置相同，所以他们的位置就是S点的位置(b)因为MVqq (0.18mV)Vq1 (0.24mV)+ Tv (0.18mV) 且 Vq2 (0.22mV)Vq1 (0.225 mV).所以Q点位于Q1点。 Vs2 (0.7mV)Vs1 (0.58 mV). 所以S点位于S2点处。
15、 注解: (1) Vq1; (2) Vq1 + Tv;(3) MVqq; (4) Vq2; (5) Vs1; (6) Vs2.图9 确定Q点2-2 在包含有50个样本点和360HZ样本频率的每个间隔中，按从左到右的顺序在每个间隔中记录和 的绝对值。因此，每个区间所需的 时间大约为0.14S。如果在和 之间有一个非零极值（绝对值），那么就将其记录下来。（如图4b中 A、B两点）。忽略没有非零极值的间隔。可以看到，在每个R-R间隔之间都有3到9个间隔（基于医学定义），而且，每个R-R间隔之间持续的时间是0.4-1.2S（如图4b和图5）。2-3 选择正确的极值点（如图4c），该步骤可分为以下两小步。测量（或）中毗邻的两正极值点间的时间差，有以下两种情况。1）、时差等于或小于间隔的持续的时间 (即0.14S或50G个采样点),则选定较大的正极值所对应的点。例如,在图6中A点和B点分别是i-1间隔和i间隔中具有最大的正极值的点，其中，B点所对应的正极值较大，所以我们选择B点。2）、时差大于间隔的持续的时间。则两个正极值都被选定。例如,在图6中，
