一种基于心电信号的特征取样方法与流程

本发明属于数据处理技术领域，涉及心电图数据的自动处理和分析，用以实现对心电图数据波形特征点的定位。

背景技术：

心脏疾病是严重危害人类健康的重大疾病，在心电信号及心率变异信号中蕴含着丰富的心脏活动信息，是心脏疾病的预防、诊断和治疗非常有价值的参考依据。

根据心电理论，对于小波变换的阀值去噪方法，易在Q、S波处出现振荡现象，消除基线漂移时可能使得心电信号低频部分的波形形态发生畸变，从而影响疾病诊断的准确性。在信号集抽样过程中，由于小波不具有平移不变性，易造成重构信号的Gibbs振荡。

RR间期即为心脏搏动周期，正常心率在60-100次/分。0.1-0.16s为正常时间。心律失常是复杂而多样的，目前对于RR间期心律失常的自动分类算法的研究还刚刚起步，心脏交感神经和迷走神经调节的变化规律不同人群分类是不同的；年轻人、老年人及房颤病人的生理数据样本分布存在差异，这就使得特征值的提取和预测恶性室性心律失常事件有一些困难。

技术实现要素：

基于上述现有技术的不足，本发明的目的在于提供一种基于心电信号的特征取样方法，该方法克服上述现有方法的缺点，充分利用非线性熵分析方法应用于HRV信号分析过程之中，来获取较好的效果，准确的计算出的心电特征点。

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

一种基于心电信号的特征取样方法，该方法可以概括为：A.读取RR间期序列；B.利用RR间期序列进行方差滤波，得到衍生截取数据；C.将衍生截取数据转换为原始IHR序列；D.将原始IHR序列进行三次重采样，以去除趋势，保留瞬时变化量，进而得到重采样IHR序列；E.利用重采样IHR序列进行HHT变换，将瞬时频率进行分解，消除低频干扰分量，获取有意义的瞬时频率；F.利用有意义的瞬时频率，Hilbert谱熵特征提取样本。

步骤E，具体是，将重采样IHR序列的信号进行小波包分解，将分解后的窄带信号继续进行EMD分解，来提取低频幅值(0.04～0.15Hz)、高频幅值(0.15～0.4Hz)、甚高频幅值(＞0.4Hz)及总幅值(0～1Hz)，作为反映生理病理变化的样本特征，再经过样本相关性判断处理后进行Hilbert变换得到Hilbert时频谱

进一步，所述读取RR间期序列，具体指，通过计算存放R峰位置序列中前后两个元素的差值就能得到RR间期序列，如当出现连续3个以上的RR间期大于1.5秒，则认为房室出现阻滞，进行重读；所述方差滤波是截取发生时长十分钟的RR间期序列和正常窦性RR间期序列进行RLS自适应算法方差，滤除干扰和异位心搏，转换为衍生截取数据。

进一步，所述将衍生截取数据转换为原始IHR序列，具体指，采样条频率为2Hz，采样的横轴为时间轴。

进一步，所述Hilbert时频谱按频率轴划分4个面积不等的频段。

进一步，所述利用瞬时频率Hilbert谱熵特征提取样本，该瞬时频段内每个能量单元为，

第k频段的Hibert时频谱熵为，

该方法的具体步骤包括如下：

(1)读取RR间期序列，通过计算存放R峰位置序列中前后两个元素的差值就能得到RR间期序列。

(2)利用RR间期序列进行方差滤波，滤除白噪音和异位干扰，得到衍生截取数据。方差滤波是截取发生时长十分钟的RR间期序列和正常窦性RR间期序列进行RLS自适应算法方差，滤除干扰和异位心搏，转换为衍生截取数据。

(3)将衍生截取数据转换为原始IHR序列。采样条频率为2HZ，采样的横轴为时间轴，具体算法为：t(n)＝RR(n)+t(n-1)n∈[1，N]，t(0)＝0，n，n∈Z。其中，t代表时间，N代表心跳数，在每个RR期间，将衍生截取数据转换为瞬时心率数列IHR(beats/min)，具体公式为：IHR＝6000/RR(ms)，即为原始IHR序列。

(4)将原始IHR序列进行三次重采样，得到均值用来去除趋势，保留瞬时变化量，进而得到重采样IHR序列，采样间隔为0.5秒。

(5)利用重采样IHR序列进行HHT变换，是根据HHT算法实现的。HHT算法具体的流程为是：将重采样IHR序列的信号进行小波包分解，将分解后的窄带信号继续进行EMD分解，来提取低频幅值(0.04～0.15Hz)、高频幅值(0.15～0.4Hz)、甚高频幅值(＞0.4Hz)及总幅值(0～1Hz)，作为反映生理病理变化的样本特征，再经过样本相关性判断处理后进行Hilbert变换，得到Hilbert时频谱。

(6)样本相关性判断处理过程为，将样本均匀等分为800点，每点表示0.001Hz，这样可将Hilbert时频谱按频率轴划分4(k＝4)个面积不等的频段。即1频段(0～0.04Hz)、2频段(0.04～0.15Hz)、3频段(0.15～0.4Hz)、4频段(＞0.4Hz)。在每个频率段中，再划分k(k＝1，2，3，4)个等面积的时频网格单元，根据Hilbert谱的时频多分辨率特性和HRV信号能量分布，表现为每个频段中时频单元的面积和划分疏密，以便准确提取生理病理特征。

(7)Hilbert变换过程是将Hilbert时频谱划分为四个不同频段的生理带，具体为，时间轴3000点，频率轴800点，各频段划分的时频网格单元数分别为1频段：40×150，2频段：6×150，3频段：5×60，4频段：2×30。从而在不同生理频带求取HRV信号的Hilbert能量棒形图，并提取各生理频带的HRV信号的能量特征为评价心率变异性的HRV时频特征，所述的HRV(Heart rate variability)是指心率节奏快慢随时间所发生的变化，即心率变异性。

与现有技术相比，本发明的有益效果是：通过HHT算法，实现了对年轻人、老年人及房颤病人的特征HRV信号的分频取样，相对传统的特征取样方法，Hilbert谱熵能更有效的区分各样本组，可更清晰地反映不同生理病理状况下的交感神经和迷走神经调节的变化规律。

附图说明

图1为本发明整体流程示意图。

图2为HHT变换的算法示意图。

图3为HHT算法中相关性判断及HRV信号的Hilbert时频谱划分示意图。

图4为样本组Hilbert边际谱谱熵分布特征示意图。

具体实施方式

为了使本发明的上述技术内容和构造特点能更容易地被本领域一般技术人员所理解，下面结合附图和技术要点对本申请进一步的说明。

实现本发明目的的技术关键是在RR间期，利用原始IHR序列，经HHT算法变换，获取不用样本组下的Hilbert谱熵来获取心电特征，其实现的具体实施步骤和解释如下：

请参阅图1～4。

读取RR间期序列，RR间期是两个连续R波之间的间隔，通过计算存放R峰位置序列中前后两个元素的差值就能得到RR间期。如当出现连续3个以上的RR间期大于1.5秒，则认为房室出现阻滞，进行重读。如，心率每4秒刷新一次，则显示的心率是当前4秒钟内的平均心率，它是以RR间期的平均值为依据的，从而消除计算过程中的随机性。

RR间期的平均值算法为：

R(i)是第i跳的RR间期的实时值。

利用RR间期序列进行方差滤波，滤除白噪音和异位干扰，得到衍生截取数据。方差滤波是截取发生时长十分钟的RR间期序列和正常窦性RR间期序列进行RLS自适应算法方差，滤除干扰和异位心搏，转换为衍生截取数据。RLS自适应算法是在噪声为50Hz浮动干扰时，利用自适应抵消处理，利用的器件为自适应工频陷波器，自适应滤波方法可根据信号的变化自动调节权矢量以获得最佳输出，其含有两个调节权值。如，原始信号为含有50Hz工频干扰的心电信号，则RR间期序列输入和心电信号同步采样，利用余弦波Bcos(ω0kTs+φ)，采样频率fs＝1/Ts，Ts为采样周期，输入r(k)的经过相移90度，变为R(k)，则r(k)＝Bcos(ω0kTs+φ)，R(k)＝Bsin(ω0kTs+φ)，得到的衍生截取数据以采样条形式记录。

将衍生截取数据转换为原始IHR序列，过程中，采样条频率为2HZ，采样的横轴为时间轴，具体算法为：t(n)＝RR(n)+t(n-1)n∈[1，N]，t(0)＝0，n，n∈Z。其中，t代表时间，N代表心跳数，在每个RR期间，将衍生截取数据转换为瞬时心率数列IHR(beats/min)，即原始IHR序列，具体公式为：IHR＝6000/RR(ms)。

将原始IHR序列进行三次重采样，得到均值用来去除趋势，保留瞬时变化量，进而得到重采样IHR序列，采样间隔设为0.5秒。

利用重采样IHR序列进行HHT变换，这是根据HHT算法实现的，所述HHT算法具体的流程为是：将重采样IHR序列的信号进行小波包分解，将分解后的窄带信号继续进行EMD分解，来提取低频幅值(0.04～0.15Hz)、高频幅值(0.15～0.4Hz)、甚高频幅值(＞0.4Hz)及总幅值(0～1Hz)，作为反映生理病理变化的样本特征，再经过样本相关性判断处理，处理后的数据进行IMF解算，再进行Hilbert变换，得到Hilbert时频谱，利用Hilert时频谱可在样本上绘出Hilbert边际谱。

样本相关性判断处理具体过程为，(参考图3)将样本均匀等分为800点，每点表示0.001Hz，这样可将Hilbert时频谱按频率轴划分4(k＝4)个面积不等频段。即1频段(0～0.04Hz)、2频段(0.04～0.15Hz)、3频段(0.15～0.4Hz)、4频段(＞0.4Hz)。在每个频率段中，再划分k(k＝1，2，3，4)个等面积的时频网格单元，这样，Hilbert时频谱的时频多分辨率特性和HRV信号能量分布，就表现为每个频段中时频单元的面积和划分疏密，以便准确提取生理病理特征。

与样本相对应的，是将Hilbert时频谱划分为四个不同频段的生理带，具体为，时间轴3000点，频率轴800点，各频段划分的时频网格单元数分别为1频段：40×150，2频段：6×150，3频段：5×60，4频段：2×30。在每个频率分段中，再划分m×k(k＝1，2，3，4)个等面积的时频网格单元，从而在不同生理频带求取HRV信号的Hilbert能量棒形图，并提取各生理频带的HRV信号的能量特征作为分析心率变异性的HRV时频特征，所述的HRV(Heart rate variability)是指心率节奏快慢随时间所发生的变化，即心率变异性。

如图4，是房颤病人，老年人，和年轻人的柱状图表示各频段谱熵的均值。由于心电理论中，第2和第3频段与交感神经和迷走神活动密切相关，年轻人和老年人的交感神经和迷走神活动水平明显不同，熵则是能更好的从复杂性的角度表明系统特征。房颤病人在整个频谱都有能量的分布，年轻人高频能量分布稀疏，老年人在频率大于0.1Hz的能量分布稀疏。

谱熵特征提取样本的具体算法是，设第k频段内的总能量为Ak，该时频段内每个能量单元p为：

再计算能量单元p的概率分布，

第k频段的Hibert时频谱熵为，

其中，β(0＜β＜1)，是不同频段Hilbert变换的权重值，∑β＝1，为了区分第2和第3频段的生理特征，在获取Hilbert时频谱熵时，可适当加大这两个频段的权重。

以上所述为本发明的实施例，并非因此限制本发明的专利范围，凡利用本发明说明书方法流程，或直接或间接运用在其他相关的技术领域，均包括在本发明的专利保护范围内。