当前位置: 首页 > news >正文

「多模态融合(ECG+PPG联合变异性分析)——ECG信号处理-第十六课」2025年6月10日

目录

1、引言

2、ECG 信号处理与 HRV 分析

2.1 ECG 信号预处理

2.2 心率变异性(HRV)分析

2.3 结果展示

3. PPG 信号处理与 PRV 分析

3.1 PPG 信号预处理

3.2 PRV 间期提取

3.3 脉搏率变异性(PRV)分析

3.4 结果展示

4. ECG 与 PPG 联合变异性分析

4.1 相关性分析

4.2 一致性分析

5.结论


1、引言

        心电信号(ECG)和光电容积脉搏波(PPG)是两种重要的生理信号,它们分别反映了心脏的电活动和外周血管的容积变化。近年来,多模态融合技术在生理信号分析领域取得了显著进展,通过联合分析ECG和PPG信号,可以获得更全面的心血管系统信息。本文将深入探讨基于ECG和PPG的联合变异性分析方法,包括信号预处理、特征提取以及多模态融合分析技术。

2、ECG 信号处理与 HRV 分析

2.1 ECG 信号预处理

  • 基线漂移滤除

基线漂移会干扰ECG信号的分析,采用低通滤波器提取基线漂移后将其从原始信号中减去。

  • 肌电干扰滤除

肌电干扰通常存在于高频段,使用带通滤波器保留1-45Hz的信号,从而去除肌电干扰,保留ECG信号中的有效成分。

  • 工频干扰滤除

工频干扰通常为50Hz的正弦波信号,使用陷波滤波器去除50Hz工频干扰,有效抑制该频率附近的干扰信号。

图1:原始心电信号与去噪后信号对比

相关内容略,详情见我博客👇

「心电信号中基线漂移、肌电干扰与工频干扰联合滤除——ECG信号处理-第十一课」2025年3月31日_肌电信号滤除心电信号-CSDN博客

2.2 心率变异性(HRV)分析

(1)HRV 时域分析

心率变异性(HRV)是指逐次心跳周期(RR间期)的差异变化。时域分析是 HRV 分析的常用方法之一,主要指标包括:

  • SDNN :所有正常RR间期的标准差,反映自主神经系统的整体调节能力。

  • NN50 :相邻RR间期差值大于 50ms 的次数。

  • pNN50 :NN50 占总 RR 间期数的百分比。

  • RMSSD :相邻 RR 间期差值的均方根,主要反映副交感神经活动。

  • 等等。。。

(2)HRV 非线性分析 - Poincaré 图

Poincaré 图是一种直观的非线性分析方法,将每个 RR 间期与其后一个 RR 间期进行二维散点绘图。主要分析指标为 SD1 和 SD2,其中 SD1 反映短期波动,SD2 反映长期波动。

  • SD1 和 SD2 的计算公式

    • SD1 = std(diff(rr_intervals)) / sqrt(2)

    • SD2 = sqrt(2 * std(rr_intervals).^2 - SD1^2)

  • 椭圆拟合 :在 Poincaré 图上绘制椭圆,椭圆的长轴和短轴分别对应 SD2 和 SD1 的方向。

(3)HRV 频域分析

频域分析通过傅里叶变换等方法,将 RR 间期序列转换到频域,分析不同频率成分的能量分布。主要频段包括:

  • VLF(0.0005 - 0.04Hz) :超低频成分,可能与长期调节机制相关。

  • LF(0.04 - 0.15Hz) :低频成分,与交感神经活动和压力反应有关。

  • HF(0.15 - 0.4Hz) :高频成分,主要由副交感神经调节。

图2:ECG 信号RR间期序列

图3:ECG 信号RR间期直方图

图4:ECG 信号非线性分析-庞加莱图

相关内容略,详情见我博客👇


「Pan-Tompkins算法(R波定位与检测)——ECG信号处理-第十课」2025年3月17日-CSDN博客

「心率变异性分析——ECG信号处理-第十二课」2025年4月1日_hrv(心率变异性)所有健康人群的散点图,如何绘制,值是rr间期-CSDN博客

2.3 结果展示

===== HRV时间域分析 =====
SDNN: 36.80 ms
NN50: 6 次
pNN50: 4.80 %
RMSSD: 27.42 ms

=====  HRV频域分析 =====
LF: 222.22 ms^2
HF: 553.61 ms^2
TP: 775.83 ms^2
LF/HF: 0.40

=====  Poincare图指标 =====
SD1: 19.46 ms
SD2: 48.26 ms
SD1/SD2: 0.40
EA: 2951.26 ms^2
VAI: 0.0560 rad
VLI: 37.19 ms

3. PPG 信号处理与 PRV 分析

3.1 PPG 信号预处理

        PPG(光电容积描记)信号是一种通过光电传感器测量血液容积变化的生理信号,广泛用于监测心率、血氧饱和度等生理参数。然而,PPG信号容易受到多种噪声的干扰,如基线漂移、肌电干扰(EMG)和工频干扰(通常为50Hz或60Hz)。为了提高信号质量和分析准确性,需要对PPG信号进行去噪处理。

图5:PPG 信号与去噪后信号对比

3.2 PRV 间期提取

        脉搏波间期(PPI)是 PPG 信号相邻波峰之间的时间间隔,类似于 ECG 的 RR 间期。通过检测 PPG 信号的波峰,提取 PPI 序列:

  • 波峰检测 :对于PPG信号,采用查找局部极大值的方法检测脉搏波的波峰,设置合适的最小峰值距离和最小峰值高度参数,确保检测到的波峰对应于真实的脉搏波搏动。

  • PPI 计算 :对检测到的波峰位置进行差分运算,得到 PPI 间期,单位转换为毫秒。

图6:PPG 信号与峰值检测

3.3 脉搏率变异性(PRV)分析

  • 时域分析:类似于HRV的时域分析,计算PPI(脉搏-脉搏间期)的标准差、NN50、pNN50和RMSSD等指标。

  • 频域分析:对PPI间期时间序列进行频谱分析,得到LF、HF和TP等频域指标。

  • 非线性分析:通过Poincaré图分析PPI间期的非线性特征,计算SD1、SD2及其比值。

图7:PPG 信号与去噪后信号对比

3.4 结果展示

===== PRV 时域分析 =====
SDNN: 36.35 ms
NN50: 6 次
pNN50: 4.84 %
RMSSD: 27.54 ms

===== PRV 频域分析 =====
LF: 215.89 ms^2
HF: 486.74 ms^2
TP: 702.63 ms^2
LF/HF: 0.44

===== PRV Poincare图指标 =====
SD1: 19.5489 ms
SD2: 47.5430 ms
SD1/SD2 Ratio: 0.4112
EA: 2919.84 ms^2
VAI: 0.0562 rad
VLI: 36.49 ms

4. ECG 与 PPG 联合变异性分析

4.1 相关性分析

        相关性系数计算:计算HRV和PRV各特征指标之间的相关性系数,如SD1、SD2、TP、SDNN、NN50、pNN50、RMSSD、LF和HF等指标之间的相关性。根据相关性系数的大小和正负,分析HRV和PRV指标之间的相关性。以散点图的形式展示HRV和PRV各特征指标之间的相关性关系。

图8:ECG和PPG 信号 Spearman分析

(1)结果展示:

===== (a) SD1 相关性 =====
Spearman 相关系数: 0.38

===== (d) SDNN 相关性 =====
Spearman 相关系数: 0.35

===== (g) RMSSD 相关性 =====
Spearman 相关系数: 0.03

===== (b) SD2 相关性 =====
Spearman 相关系数: 0.32

===== (e) NN50 相关性 =====
Spearman 相关系数: 0.79

===== (h) LF 相关性 =====
Spearman 相关系数: 0.72

===== (c) TP 相关性 =====
Spearman 相关系数: 0.38

===== (f) pNN50 相关性 =====
Spearman 相关系数: 0.01

===== (i) HF 相关性 =====
Spearman 相关系数: 0.31

(2)结果分析:

从输出的相关系数来看,不同HRV和PRV特征之间的相关性存在较大差异:

强相关(相关系数>0.7)

  • NN50:Spearman相关系数为0.79,表明HRV和PRV的NN50特征之间存在较强的正相关关系。这意味着在心率变异性增加时,脉搏波变异性也在增加,且这种变化在时间域的NN50指标上表现得较为明显。

  • LF(低频功率):相关系数为0.72,同样表现出较强的相关性。低频功率通常与交感神经活动有关,这表明HRV和PRV在反映交感神经调节方面具有一致性。

中等相关(0.4<相关系数<0.7)

  • SD1:相关系数为0.38,表明HRV和PRV的SD1特征之间存在中等程度的正相关。SD1主要反映短期心率变异性,在短期波动方面,ECG和PPG信号具有一定的相关性。

  • TP(总功率):相关系数也为0.38,显示HRV和PRV的总功率之间存在中等相关性。总功率反映了信号的整体变异性,这表明两种信号在整体变异性上具有一致的趋势。

  • SD2:相关系数为0.32,表明HRV和PRV的SD2特征之间存在一定的正相关关系。SD2主要反映长期心率变异性,这一相关性反映了ECG和PPG信号在长期波动模式上的相似性。

  • HF(高频功率):相关系数为0.31,处于中等相关水平。高频功率通常与副交感神经活动相关,这提示在副交感神经调节方面,HRV和PRV具有一定的一致性。

弱相关或无相关(相关系数<0.4)

  • RMSSD:相关系数仅为0.03,表明HRV和PRV的RMSSD特征之间几乎不存在相关性。RMSSD主要用于评估副交感神经活动的短期变化,在ECG和PPG信号中可能受到不同因素的影响,导致其相关性较低。

  • pNN50:相关系数为0.01,同样表现出无相关性。pNN50是基于NN50的相对指标,用于评估心率变异性的变化比例。在ECG和PPG信号中,这一指标的变化可能受到多种因素干扰,导致相关性不明显。

4.2 一致性分析

        Bland-Altman分析:绘制Bland-Altman图,评估HRV和PRV指标之间的一致性。通过计算两者的差异均值和标准差,绘制差异与平均值的关系图,并确定一致性界限。观察Bland-Altman图中数据点的分布情况,评估HRV和PRV指标之间的一致性程度。

图9:ECG和PPG 信号 Bland-Altman分析

(1)结果展示

===== (a) SD1 Bland-Altman =====
均值差异: 13.385;标准差: 40.251;一致性上限: 92.277;一致性下限: -65.507;95%一致性范围内的数据点百分比: 93.86%

===== (d) SDNN Bland-Altman =====
均值差异: 11.011;标准差: 38.002;一致性上限: 85.495;一致性下限: -63.474;95%一致性范围内的数据点百分比: 94.74%

===== (g) RMSSD Bland-Altman =====
均值差异: 621.583;标准差: 4751.322;一致性上限: 9934.175;一致性下限: -8691.009;95%一致性范围内的数据点百分比: 99.12%

===== (b) SD2 Bland-Altman =====
均值差异: 11.725;标准差: 39.541;一致性上限: 89.225;一致性下限: -65.776;95%一致性范围内的数据点百分比: 94.74%

===== (e) NN50 Bland-Altman =====
均值差异: 5.328;标准差: 24.052;一致性上限: 52.470;一致性下限: -41.815;95%一致性范围内的数据点百分比: 93.10%

===== (h) LF Bland-Altman =====
均值差异: 7.103;标准差: 18.116;一致性上限: 42.610;一致性下限: -28.404;95%一致性范围内的数据点百分比: 94.83%

===== (c) TP Bland-Altman =====
均值差异: 18.866;标准差: 56.739;一致性上限: 130.074;一致性下限: -92.342;95%一致性范围内的数据点百分比: 93.86%

===== (f) pNN50 Bland-Altman =====
均值差异: 228.590;标准差: 1891.950;一致性上限: 3936.812;一致性下限: -3479.632;95%一致性范围内的数据点百分比: 99.12%

===== (i) HF Bland-Altman =====
均值差异: 179.162;标准差: 904.552;一致性上限: 1952.085;一致性下限: -1593.761;95%一致性范围内的数据点百分比: 95.61%

(2)结果分析

总体而言,大部分特征的Bland-Altman分析显示HRV和PRV之间具有较好的一致性,尤其在SD1、SDNN、NN50、LF和TP等特征上,数据点的绝大部分落在一致性区间内,表明这些特征在两种信号间有较好的可比性。不过,RMSSD和pNN50的离散度较高,一致性稍差。

5.结论

        本文详细介绍了 ECG 和 PPG 信号的处理流程,以及基于这两种信号的变异性分析方法,包括时间域、非线性分析(Poincaré 图)和频域分析。通过联合分析 HRV 和 PRV,探讨了多模态数据融合在心血管评估中的应用。相关性分析和一致性分析结果表明,HRV 和 PRV 指标之间存在一定的关联性,这为进一步开展多模态融合研究提供了基础。

Tips:下一讲,我们将进一步探讨,心电信号处理与应用的其他内容。

以上就是多模态融合(ECG+PPG联合变异性分析)的全部内容啦~

我们下期再见,拜拜(⭐v⭐) ~

(Ps:有代码实现需求,请见主页信息,谢谢支持!~)

http://www.xdnf.cn/news/963199.html

相关文章:

  • 屋顶变身“发电站” ,中天合创屋面分布式光伏发电项目顺利并网!
  • 校园二手交易平台(微信小程序版)
  • AtCoder 第409​场初级竞赛 A~E题解
  • pycharm最近遇到的一些问题
  • Linux 阻塞非阻塞
  • Ubuntu24.04常见问题
  • AbMole推荐:Cisplatin,从机制到应用的全面了解
  • Hive实现值列表横向展示(非列转行)
  • 汽车扭杆综合试验机
  • 机器人模型文件urdf介绍
  • C++基础学习:深入理解类中的构造函数、析构函数、this指针与new关键字
  • JWT令牌如何在FastAPI中实现安全又高效的生成与验证?
  • 解锁数据库简洁之道:FastAPI与SQLModel实战指南
  • 1.2 git使用
  • 小程序中的状态管理库-mobx-miniprogram
  • 【深尚想】SX1276IMLTRT LoRa射频收发器芯片 物联网 电子元器件解析
  • 在鸿蒙HarmonyOS 5中实现抖音风格的草稿箱功能
  • Kafka 的容错与持久性:副本复制机制的工作原理与实践
  • 实现安全、经济、节约、环保的智慧交通开源了
  • JBank:Jucoin 推出的 Web3 去中心化自托管银行金融协议
  • error: Sandbox: rsync(17136) deny(1) file-write-create
  • 汽车生产虚拟实训中的技能提升与生产优化​
  • 机器学习的可解释性
  • 项目课题——智能花盆系统设计
  • “机器学习中的‘Hello World‘:为什么我们总用MNIST数据集,以及何时该放弃它“
  • 机器学习中的优化问题描述
  • 在多云环境透析连接ngx_stream_proxy_protocol_vendor_module
  • ffmpeg 新版本转码设置帧率上限
  • 搭建gitlab ci/cd runner实现对c++项目的自动编译和打包
  • 51c嵌入式※~电路~合集32~PWM