专利名称：脉搏波形特征点提取方法
技术领域：
本发明涉及一种基于同步采集的脉搏心电信号，提取脉搏波形特征点的方法。
背景技术：
目前，脉搏波检测与波形特征的提取主要采用小波分析的模极大值法和微分法。 小波分析的模极大值法利用小波变换突出局部特征的能力，使它成为检测瞬态突变及信号边沿的有力手段，但是此方法对不规则脉搏波信号检测率效果不理想；微分法运算速度快， 对脉搏波主波的检测率较高，但对重搏前波，重搏波等特征参数的提取不是很理想。脉搏信号的有用信号频段较窄，信号突变不明显，各个波峰没有明显特征。正常状态下的规则脉搏波，通常脉搏主波峰值最大，但在受试者处于疲劳等状态时，波峰、重搏前波和重搏波的幅值非常接近，甚至出现重搏波峰值高于主波峰值的情况，同时重搏波差分值和主波差分值比较接近，这些都给脉搏波波形特征点精确检测带来极大困难。
发明内容
本发明的目的是准确检测脉搏波的重搏波和脉搏主波。脉搏波形特征点提取方法，其方法的步骤为
(1)同步采集心电信号和脉搏信号；
(2)对采集到的心电信号和脉搏信号进行去基线漂移，去工频干扰，去肌电干扰，消除
噪声；
(3)检测分析心电信号，提取心电特征信息，准确判断出心电信号的R波和T波位置；
(4)基于心电信号R波和T波位置确定脉搏信号主波、降中峡和重搏波位置。由于心电信号突变明显，比较容易检测识别，本发明在同步采集脉搏信号和心电信号的基础上，对心电信号进行分析，检测出心电信号R波和T波，利用R波和T波位置来定位检测脉搏波重搏波和主波。该方法可以准确检测规则和不规则脉搏波，提取波形特征。这种方法不但能够减弱强干扰和不规则脉搏对脉搏波检测带来的不利影响，还能定量描述心电脉搏特征点在时域上的关系。
下面结合附图和实施方式对本发明作进一步详细的说明。图1是本发明同步心电脉搏采集系统框图。图2是本发明基于心电信号的脉搏波形特征点提取波形图。
具体实施例方式本发明是脉搏波形特征点提取方法，其方法的步骤为
(1)同步采集心电信号和脉搏信号；
(2)对采集到的心电信号和脉搏信号进行去基线漂移，去工频干扰，去肌电干扰，消除噪声；
(3)检测分析心电信号，提取心电特征信息，准确判断出心电信号的R波和T波位置；
(4)基于心电信号R波和T波位置确定脉搏信号主波、降中峡和重搏波位置。下面结合附图进一步展开本发明。脉搏波形特征点提取的方法，首先同步采集到心电信号和脉搏信号，图1所示为心电脉搏信号同步采集系统框图；数据采集系统由心电电极、光电脉搏传感器、信号调理电路、MP425数据采集卡和LabVIEW虚拟仪器平台构成，实现对信号的实时采集、显示与存储。 同步采集被试者ECG (Electrocardio -gram)信号和光电脉搏波信号(Pulse wave)，脉搏和心电信号的采样频率为1000Hz。然后对实验采集到的心电信号和脉搏信号进行预处理。先确定出低通滤波器的初始条件，将心电信号通过滤波器，所得结果反转后再次通过滤波器，最后将所得结果再反转，得到心电信号的精确零相位低频噪声序列，用原始心电信号减去这个精确零相位低频噪声序列，所得结果就是心电信号去基线漂移后的结果，截止频率f。为0. 5Hz,既消除频率很低的基线漂移，又保持ECG信号其它成份的形态不受影响。去基线漂移后的心电信号通过Notch滤波器进行去工频干扰。然后对得到的心电信号进行小波分解，剔除与噪声相关的系数，重构心电信号，完成去除肌电干扰。将得到的心电信号进行二次小波包降噪，使用启发式阈值对心电信号进行消噪处理。对采集到的脉搏信号同样进行预处理，先确定出低通滤波器的初始条件，将脉搏信号通过滤波器，所得结果反转后再次通过滤波器，最后将所得结果再反转，得到脉搏信号的精确零相位低频噪声序列，用原始脉搏信号减去这个精确零相位低频噪声序列，所得结果就是脉搏信号去基线漂移后的结果。去基线漂移后的脉搏信号通过Notch滤波器进行去工频干扰，然后对得到的脉搏信号进行小波分解，剔除与噪声相关的系数，重构脉搏信号，完成去除肌电干扰。对含有噪声的脉搏信号进行小波分解，通过估计的最大风险最小化的极大极小阈值法进行小波去噪。采用TEO (Teager Energy Operator)能量算子法检测心电信号R波；在精确找出 R波的位置后，利用差分法找出Q、S波位置和，以及RR间期；由向前滑动0. 55取最大值即可得出T波位置，如图2 (b)中所示。对降噪后的脉搏信号y2(η)进行后向差分运算
权利要求
1.脉搏波形特征点提取方法，其步骤为(1)同步采集心电信号和脉搏信号；(2)对采集到的心电信号和脉搏信号进行去基线漂移，去工频干扰，去肌电干扰，消除噪声；(3)检测分析心电信号，提取心电特征信息，准确判断出心电信号的R波和T波位置；(4)基于心电信号R波和T波位置确定脉搏信号主波、降中峡和重搏波位置。
2.根据权利要求1所述的脉搏波形特征点提取方法，其特征在于对采集到的心电信号进行预处理，先确定出低通滤波器的初始条件，将心电信号通过滤波器，所得结果反转后再次通过滤波器，最后将所得结果再反转，得到心电信号的精确零相位低频噪声序列，用原始心电信号减去这个精确零相位低频噪声序列，所得结果就是心电信号去基线漂移后的结果；去基线漂移后的心电信号通过Notch滤波器进行去工频干扰，然后对得到的心电信号进行小波分解，剔除与噪声相关的系数，重构心电信号，完成去除肌电干扰；将得到的心电信号进行二次小波包降噪，使用启发式阈值对心电信号进行消噪处理。
3.根据权利要求1所述的脉搏波形特征点提取方法，其特征在于对采集到的脉搏信号进行预处理，先确定出低通滤波器的初始条件，将脉搏信号通过滤波器，所得结果反转后再次通过滤波器，最后将所得结果再反转，得到脉搏信号的精确零相位低频噪声序列，用原始脉搏信号减去这个精确零相位低频噪声序列，所得结果就是脉搏信号去基线漂移后的结果；去基线漂移后的脉搏信号通过Notch滤波器进行去工频干扰，然后对得到的脉搏信号进行小波分解，剔除与噪声相关的系数，重构脉搏信号，完成去除肌电干扰。
4.根据权利要求1所述的脉搏波形特征点提取方法，其特征在于对降噪后的脉搏信号进行后向差分运算，然后对差分运算结果进行平滑滤波，用于进一步消除高频脉冲噪声、 窄带噪声和工频干扰，得到消噪后的脉搏信号差分序列。
5.根据权利要求1所述的脉搏波形特征点提取方法，其特征在于基于心电信号和脉搏信号在时域上的同步性，在消噪后的脉搏信号差分序列上找出T点对应位置，分别前向后向滑动若干个点，找到脉搏波信号差分信号的局部极大值位置，检测出这个局部极大值前后最近的两个过零点，这两个过零点分别是脉搏波的波谷和波峰。
6.根据权利要求1所述的脉搏波形特征点提取方法，其特征在于基于心电信号和脉搏信号在时域上的同步性，在消噪后的脉搏信号差分序列上找出R点对应位置，分别前向后向滑动若干个点，找到脉搏波信号差分信号的局部极大值位置，检测出这个局部极大值前后最近的两个过零点，这两个过零点分别对应脉搏波的降中峡位置和脉搏波的重搏波位置。
全文摘要
脉搏波形特征点提取方法，其目的是准确检测脉搏波的重搏波和脉搏主波；其方法的步骤为(1)同步采集心电信号和脉搏信号；(2)对采集到的心电信号和脉搏信号进行去基线漂移，去工频干扰，去肌电干扰，消除噪声；(3)检测分析心电信号，提取心电特征信息，准确判断出心电信号的R波和T波位置；(4)基于心电信号R波和T波位置确定脉搏信号主波、降中峡和重搏波位置；其系统包括一种生理信号数据采集卡，用于同步采集心电信号和脉搏信号；一种用于检测人体心电信号的心电检测模块；一种用于检测人体脉搏信号的脉搏检测模块。
文档编号A61B5/02GK102440768SQ201010505380
公开日2012年5月9日 申请日期2010年10月13日 优先权日2010年10月13日
发明者张爱华, 柳倩, 赵治月 申请人:兰州理工大学