一种基于拍内特征子波复杂网络的室性早搏检测方法技术

本发明专利技术专利公开了一种基于拍内特征子波复杂网络的室性早搏检测方法，首先通过心电信号预处理对采集到的心电信号进行去噪和QRS波检测，将预处理后的心电信号进行心拍分割以获得心拍信号，对单个心拍内部子波分割并计算各部分的特征；对于得到的特征序列采用去趋势方法消除了长时程的心电信号特征在不同状态下的区域性差异，并去除微小波动，以消除了微小波动对后续序列编码的影响；利用符号动力学进行编码，并构建拍内特征子波复杂网络，计算复杂网络的网络特征参数，最终提取能够反映心电信号特点的心电特征，进一步通过机器学习分类算法实现心电信号的分类。算法实现心电信号的分类。算法实现心电信号的分类。

【技术实现步骤摘要】
一种基于拍内特征子波复杂网络的室性早搏检测方法


[0001]本专利技术涉及一种心电信号处理技术，属于长时程心电监护设备、医疗监护
，具体涉及一种基于拍内特征子波复杂网络的室性早搏检测方法。

技术介绍

[0002]心脏作为人体至关重要的器官，心脏疾病能显著地影响人体身心健康。据《中国心血管健康与疾病报告2021》报告，心血管病在中国城乡居民疾病死亡构成比中占首位，且处于持续增长阶段。目前我国心律失常患者约有2000万人，其发病率呈现出逐年增高的趋势。长期心律失常会增加死亡风险，早期发现心脏异常，结合适当治疗，有助于提高患者的康复率。室性早搏(Premature ventricular contractions,PVCs)是最常见的心律失常类型之一，常见于冠心病、心肌病和风湿性心脏病等疾病中。严重器质性心脏病患者，长期频繁的室性早搏可引起心绞痛、低血压或心力衰竭。
[0003]心电图是一种广泛用于诊断心脏相关疾病的工具，包含了大量反映心脏各部位状态的生理和病理信息，通过它可以检测心血管疾病引起的异常。长程动态心电监测可以捕获常规不易发现的偶发性、短阵性心电异常，明确心律失常的类型、发生频率以及风险程度等，但由于长时程心电数据量十分庞大，极大增加了医生的负担。随着心电信号自动分析技术的不断进步，心电信号预识别的功能大大减轻了医生诊断的工作量，帮助医务人员及时发现心律失常的异常现象。目前常用于心律失常检测算法中的心电信号特征选择包括心拍波形特征如QRS波宽度、R波高度和宽度、RR间期等，或是高阶统计量提取的非线性特征，这类手动提取特征可解释性强，但存在检测率易受噪声影响，鲁棒性差的问题，亟待解决。

技术实现思路

[0004]技术问题：
[0005]针对以上问题，本专利技术提供了一种基于拍内特征子波复杂网络的室性早搏检测方法及应用，用于实现长时程、抗噪、精准的心电信号分类，有效地提高心电信号自动分析的准确性，辅助医生诊断。本专利技术通过去趋势处理消除了长时程心电信号在不同状态下的区域性差异，通过去微小波动处理和符号动力学编码降低了心电信号分析的对噪声的敏感度；通过构建拍内特征子波复杂网络，利用复杂网络拓扑结构的特征反映心拍内部特征序列动力特性，提取更典型、反映本质的心电信号特征，提高检测准确率和鲁棒性。
[0006]技术方案：
[0007]为实现上述目的，本专利技术采用如下技术方案：一种基于拍内特征子波复杂网络的室性早搏检测方法及应用，包括以下步骤：
[0008]S1，心电信号预处理，包括去噪处理和QRS检测；
[0009]S2，心拍分割，将经过预处理后的心电信号以R波为基准点，分别向前向后以固定长度的窗口大小截取心拍信号；心拍内部子波分割，将每一段心电信号内部分为8段心电采样点序列；
[0010]S3，计算心拍内部每一段心电采样点序列的特征；
[0011]S4，特征去趋势处理，用以消除长时程的心电数据中不同状态下特征的区域性差异；并进行特征去微小波动处理，用以消除特征值的波动对后续编码的干扰；
[0012]S5，符号动力学编码，利用符号动力学对心拍内部特征序列进行编码，形成编码序列以便后续分析；
[0013]S6，构建拍内特征子波复杂网络，将心拍内部特征编码序列映射为复杂网络；
[0014]S7，计算复杂网络的特征参数，所述特征可以为复杂网络的边缘度、图能量、自我循环数量、节点平均度、大小为3的循环数量、大小为4的循环数量、图形链接密度、平均节点接近中心性和平均节点聚类系数等；
[0015]S8，采用上述提取的心电特征，基于机器学习方法训练早搏/非早搏分类模型。
[0016]作为本专利技术的一种改进，所述步骤S1中，采用FIR低通滤波器提取基线漂移，去除原始信号中的基线漂移；再采用FFT带通滤波器对去除基线漂移的信号进行滤波处理。心电信号预处理的QRS检测，采用Pan
‑
Tompkins方法进行QRS波检测获取R波位置。
[0017]作为本专利技术的一种改进，所述步骤S2中，心拍分割对经过去噪处理后的心电信号进行心拍分割，以R波位置为基准点向前截取0.1s以保证截取到完整的QRS波，向后截取0.4s长度，保留完整的T波同时避免截取到下一个心拍。
[0018]作为本专利技术的一种改进，所述步骤S2中，将上述截取的第n段心电信号内部按照波形进行8段的子波分割，包括以P波波峰为分界的P波前段和后段、PR间期、QRS波群中以R波波峰为界的QR波和RS波、ST间期、T波波峰为分界的T波前段和后段，按时间顺序记为波峰为界的QR波和RS波、ST间期、T波波峰为分界的T波前段和后段，按时间顺序记为为第n个心拍内部第i段心电采样点序列，其中1≤i≤8。
[0019]作为本专利技术的一种改进，所述步骤S3中，计算心拍内部每一段心电采样点序列的特征，可以选择形态特征如峰峰值、偏度、翘度等，也可以选择统计学特征如平均值、标准差等，记为等，记为为第n个心拍内部第i段心电采样点序列的特征，其中1≤i≤8。
[0020]作为本专利技术的一种改进，所述步骤S4中，对第n个心拍第i段心电采样点序列的特征值进行去趋势处理，计算局部特征均值(以当前特征值对应的心拍序号n分别向前、向后选取1个心拍内部对应采样点序列的特征值以连续3个特征值计算当前心拍对应采样点序列的局部特征均值)，以当前特征值减去局部特征均值作为特征值的取值，其计算公式为其中1≤i≤8。
[0021]作为本专利技术的一种改进，所述步骤S4中，采用分段函数映射后的取值来去除第n个心拍第i段心电采样点序列的特征值微小波动的干扰。计算局部特征均值局部特征标准差(以当前特征值对应的心拍序号n分别向前、向后截取1个心拍内部对应采样点序列的特征值以连续3个特征值计算当前心拍对应采样点序列的局部特征均值和局部特征标准差)。当特征值与局部均值的距离小于或等于局部标准差时，对应特征值取为局部均值大小；当特征值与局部均值的距离大于局部标准差时，对应特征值不变。上述过程可以理解为特征值通过分段函数映射后进行取
值，如图2所示，其计算公式为：其中1≤i≤8。
[0022]作为本专利技术的一种改进，所述步骤S5中，利用符号动力学对第n个心拍内部特征序列进行编码，将原始单值特征转化为仅包含有限取值的符号序列，是一个粗粒化过程，降低了信号分析的噪声敏感度。编码过程利用有限长的符号将特征序列形成编码序列以便后续分析,为第n个心拍内部第i段心电采样点序列的特征编码值，其中1≤i≤8。对于给定的数据长度为L的特征值序列可以通过下式将其转换为符号序列其计算公式为：其中l表示符号集长度，通常在2
‑
6之间选择；{X
i
}(1≤i≤l
‑
1)为设定的门限值的集合。
[0023]作为本专利技术的一种改进，所述步骤S6中，构建心拍内部特征复杂网络，将第n个心拍内部特征编码序列映射为复杂网络。首先，根据长度m对编码序列进行划分，取一系列嵌入向量其中然后，每个z...

【技术保护点】

【技术特征摘要】
1.一种基于拍内特征子波复杂网络的室性早搏检测方法，其特征在于，包括以下步骤：S1，心电信号预处理，包括去噪处理和QRS检测；S2，心拍分割，将经过预处理后的心电信号以R波为基准点，分别向前向后以固定长度的窗口大小截取心拍信号；进行心拍内部子波分割，将每一段心电信号内部按照子波分割为8段心电采样点序列；S3，计算心拍内部每一段心电采样点序列的特征；S4，特征去趋势处理，用以消除长时程的心电数据中不同状态下特征的区域性差异；并且进行特征去微小波动处理，用以消除特征值的波动对后续编码的干扰；S5，符号动力学编码，利用符号动力学对心拍内部特征序列进行编码，形成编码序列以便后续分析；S6，构建拍内特征子波复杂网络，将心拍内部特征编码序列映射为复杂网络；S7，计算复杂网络的特征参数，所述特征参数为复杂网络的边缘度、图能量、自我循环数量、节点平均度、大小为3的循环数量、大小为4的循环数量、图形链接密度、平均节点接近中心性和平均节点聚类系数；S8，采用上述提取的心电特征，基于机器学习方法训练早搏/非早搏分类模型。2.根据权利要求1所述的一种基于拍内特征子波复杂网络的室性早搏检测方法，其特征在于，所述步骤S1中，心电信号预处理的去噪处理采用FIR低通滤波器提取基线漂移，去除原始信号中的基线漂移；再采用FFT带通滤波器对去除基线漂移的信号进行滤波处理；心电信号预处理的QRS检测，采用Pan
‑
Tompkins方法进行QRS波检测获取R波位置。3.根据权利要求1所述的一种基于拍内特征子波复杂网络的室性早搏检测方法，其特征在于，所述步骤S2中，心拍分割将经过去噪处理后的心电信号进行心拍分割，以R波位置为基准点向前截取0.1s以保证截取到完整的QRS波，向后截取0.4s长度，保留完整的T波同时避免截取到下一个心拍。4.根据权利要求1所述的一种基于拍内特征子波复杂网络的室性早搏检测方法，其特征在于，所述步骤S2中，心拍内部子波分割是将上述截取的第n段心电信号内部按照波形进行8段的子波分割，包括以P波波峰为分界的P波前段和后段、PR间期、QRS波群中以R波波峰为界的QR波和RS波、ST间期、T波波峰为分界的T波前段和后段，按时间顺序记为为界的QR波和RS波、ST间期、T波波峰为分界的T波前段和后段，按时间顺序记为为第n个心拍内部第i段心电采样点序列，其中1≤i≤8。5.根据权利要求1所述的一种基于拍内特征子波复杂网络的室性早搏检测方法，其特征在于，所述步骤S3中，计算心拍内部每一段心电采样点序列的特征，选择形态特征如峰峰值、偏度、翘度，或选择统计学特征如平均值、标准差，记为f
in
为第n个心拍内部第i段心电采样点序列的特征，其中1≤i≤8。6.根据权利要求1所述的一种基于拍内特征子波复杂网络的室性早搏检测方法，其特征在于，所述步骤S4中，对第n个心拍第i段心电采样点序列的特征值f
in
(1≤i≤8)进行去趋势处理，计算局部特征均值以当前特征值f
in
对应的心拍序号n分别向前、向后选取1个心拍内部对应采样点序列的特征值f
in
‑1,f
in+1
，以连续3个特征值计算当前心拍对应采样点序列的局部特征均值以当前特征值f
in
减去局部特征均值作为特征值f
in
的取值，公式如下，其中1≤i≤8：
7.根据权利要求1所述的一种基于拍内特征子波复杂网络的室性早搏检测方法，其特征在于，所述步骤S4中，采用分段函数映射后的取值来去除第n个心拍第i段心电采样点序列的特征值f
in
(1≤i≤8)微小波动的干扰；计算局部特征均值局部特征标准差以当前特征值f
in
对应的心拍序号n分别向前、向后截取1个心拍内部对应采样点序列的特征值f
in
‑1,f
in+1
，以连续3个特征值计算当前心拍对应采样点序列的局部特征均值和局部特征标准差当特征值f
in
与局部均值的距离小于或等于局部标准差时，对应特征值f
in
取为局部均值大小；当特征值f
in
与局部均值的距离大于局部标准差时，对应特征值f
in
不变；上述过程为特征值通过分段函数映射后进行取值，公式如下，其中1≤i≤8：不变；上述过程为特征值通过分段函数映射后进行取值，公式如下，其中1≤i≤8：不变；上述过程为特征值通过分段函数映射后进行取值，公式如下，其中1≤i≤8：8.根据权利要求1所述的一种基于拍内特征子波复杂网络的室性早搏检测方法，其特征在于，...

【专利技术属性】
技术研发人员：蔡志鹏，沈毓敏，刘澄玉，李建清，
申请(专利权)人：东南大学，
类型：发明
国别省市：