心电信号的处理.ppt

biomedical signal processing,1,nankai university, cy li, 2019/11/1,ch4: 心电信号的处理,心电图的产生, 心电图处理的基本思路, 时间上: 动态和静态, 空间上: 心肌电特性的空间离散度, 心电图处理和分析的发展,biomedical signal processing,2,nankai university, cy li, 2019/11/1,心电信号的形成（1）,心肌细胞的动作电位 心肌细胞除极和复极的电生理现象, 极化状态(polarised): 静息电位(resting potential), 动作电位(action potential, ap): 除极(depolarisation)和复极(repolarisation),biomedical signal processing,3,nankai university, cy li, 2019/11/1,心电信号的形成（2）,电兴奋的传导(conduction or spread of electrical excitation) 窦房结: 心脏的起搏兴奋点, 其细胞自发产生50-100次/分的可传导ap, 心房的传导: 心房传导束右心房左心房 传导系统的传播: 房室束希氏束左、右分支普金野（purkinje）纤维网-心室肌, 电兴奋通过purkinje网使心室肌细胞兴奋，,biomedical signal processing,4,nankai university, cy li, 2019/11/1,心脏传导系统,biomedical signal processing,5,nankai university, cy li, 2019/11/1,english terms,窦房结: sinus node, 心房: atrium, 心室: ventricle, 房室结: atrioventricular node (junction), 希氏束: bundle of his, 普金野纤维: pukinje fibres, 心脏表面(心外膜): epicardium, 心内膜：endocardium 体表：body surface,biomedical signal processing,6,nankai university, cy li, 2019/11/1,心电信号的形成（3）,心电图(electrocardiogram, ecg)的产生 电流源：每个心肌细胞的除极和复极过程等效于一个偶极子层（dipole layer） 容积导体(volume conductor):人体组织是导电的，看作是一个容积导体， 心电向量: 所有心肌细胞的偶极子场的向量和, 心电图：所有心肌细胞的偶极子场在容积导体内产生电场，从而有电位差产生，即心电图。 体表心电图, 心外膜、心内膜电图, 希氏束电图，,biomedical signal processing,7,nankai university, cy li, 2019/11/1,典型心电信号波形,p, qrs, st, t,s-t,段,qrs,t,p,qt,biomedical signal processing,8,nankai university, cy li, 2019/11/1,心电图的记录 (recording of ecg),标准导联, 肢体导联, 胸导联, 体表多部位标测(body surface mapping), 32-512通道的数据采集系统(multi-channel data acquisition system), 电极背心, (electrode vest) 心脏表面的多部位标测(epicardial mapping), 多通道的数据采集系统 电极套(electrode sock) 心脏内膜的多部位标测(endocardial mapping), 多通道的数据采集系统 导管(catheter)和伞电极,biomedical signal processing,9,nankai university, cy li, 2019/11/1,body surface mapping,biomedical signal processing,10,nankai university, cy li, 2019/11/1,epicardial electrode sock,biomedical signal processing,11,nankai university, cy li, 2019/11/1,标准导联(1),12标准导联(standard leads) 标准i， ii，iii肢体导联(bipolar limb leads),r: right arm l: left arm f: left foot i=el-er, ii=ef-er, iii=ef-el,biomedical signal processing,12,nankai university, cy li, 2019/11/1,标准导联(2),加压肢体导联标准avr, avl, avf (augmented unipolar limb leads),l,wilson terminal: c点的电位, avl和avf的连接相似,biomedical signal processing,13,nankai university, cy li, 2019/11/1,标准导联(3),单极胸导联v1-v6 (unipolar precordial leads), precordial: chest wall, 中心电端(wilson terminal) v1-v6, 胸前电极分别与中心电端的电位差,biomedical signal processing,14,nankai university, cy li, 2019/11/1,心电信号的畸变,心脏的病变： 传导阻滞，早搏，室颤、房颤, 心肌缺血、梗塞等, qrs变宽, st段位移出现, 心率变化等, 来自心脏外的干扰信号： 50hz工频干扰 肌电干扰， 10-300hz, 呼吸的干扰，使基线漂移加剧，,biomedical signal processing,15,nankai university, cy li, 2019/11/1,心电信号的预处理,抑制工频干扰， 基线纠漂，,biomedical signal processing,16,nankai university, cy li, 2019/11/1,噪声抑制和基线漂(detrending),低通滤波器 可以滤掉心电信号中的肌电信号的高频干扰， 抑制工频干扰 自适应滤波抑制工频干扰， 基线漂移的纠正 抵消法纠漂， 基线纠漂滤波器， 0.7hz的高通截止频率，,biomedical signal processing,17,nankai university, cy li, 2019/11/1,移动平滑滤波器(低通滤波),biomedical signal processing,18,nankai university, cy li, 2019/11/1,h(z)和h(z)*h(z)的比较,h(z): 旁瓣太大, 13db; 相频特性虽在通带内保持线性, 但在进入阻带后有突变, 有可能造成心电信号的高频相位失真; h(z)*h(z): 旁瓣有较大的衰减, 26.7db, 有真正的线性相位;,biomedical signal processing,19,nankai university, cy li, 2019/11/1,移动平滑滤波器(n=8),fs=1000hz n=8,h(z), -3db, 55hz, the first side peak: -13db, nonlinear phase angle h(z)*h(z), -3db, 55hz, the first side peak: -26.7db, linear phase angle,biomedical signal processing,20,nankai university, cy li, 2019/11/1,移动平滑滤波器(n=33),h(z): -3db: 0.08hz, the first side peak: -13db h(z)*h(z): -3db: 0.08hz, the first side peak: -26.7db,biomedical signal processing,21,nankai university, cy li, 2019/11/1,移动平滑滤波器的效果,biomedical signal processing,22,nankai university, cy li, 2019/11/1,ma的h(z)和h(z)*h(z)随n的变化,biomedical signal processing,23,nankai university, cy li, 2019/11/1,自适应滤波,根据输入信号自动调节滤波器的参数, 使其性能指标最优化， 适用于对信号和噪声无先验知识(频谱)的或非平稳信号, 基本结构: 滤波器、优化指标算法、滤波器参数修改算法， 滤波器： fir, iir， 优化准则：信噪比最高、输出误差均方差最小， 参数修改：递归、非递归，,biomedical signal processing,24,nankai university, cy li, 2019/11/1,y(t),自适应消噪声,+,滤波器,修改参数,信号源,噪声源,x(t)=s(t)+n0(t),x(t),+,-,n1(t),n0(t): 噪声源（如50hz工频信号），n1(t)：通过某一未知网络的 同一噪声源，n0(t)：滤波器输出的估值， x(t)=s(t)+ n0(t)- y(t)=s(t)+ n0(t)- n0(t) 要使得x(t)在最小均方差意义上与信号s(t)最佳匹配，,biomedical signal processing,25,nankai university, cy li, 2019/11/1,基线纠漂,抵消法：从信号中减去基线的估计值, 心电基线漂移可看成某种超低频干扰信号, x(n)=s(n)+w(n), 用w(n)来逼近w(n), 并从s(n)中减去, s(n)=x(n)-w(n), pq段为基线，常以pq段的中点作为基准点，用多项式拟合基线的估计函数， 如最简单的线性拟合、三次样条函数拟合等，,biomedical signal processing,26,nankai university, cy li, 2019/11/1,基线纠漂滤波器,基线纠漂滤波器(hp) 高通滤波器，高通截止频率：0.7hz， 即以40次/分为心脏搏动过缓的下限， 基线纠漂及50hz陷波滤波器, 理想的心电信号预处理滤波器 基线纠漂+50hz限波, 可用频率抽样法设计这种滤波器,f,h(f),50,100,150,fk,biomedical signal processing,27,nankai university, cy li, 2019/11/1,纠漂和50hz陷波的效果,biomedical signal processing,28,nankai university, cy li, 2019/11/1,整系数基线纠漂滤波器,梳状滤波器的数学表示;,biomedical signal processing,29,nankai university, cy li, 2019/11/1,整系数基线纠漂滤波器(2),梳状滤波器的数学表示;,biomedical signal processing,30,nankai university, cy li, 2019/11/1,整系数filter的效果,biomedical signal processing,31,nankai university, cy li, 2019/11/1,心电信号的分析,特征波形的检测 qrs波的检测，t波的检测，st段位移， 心电图的频谱分析 fft, ar model, 时频分析, wavelet transformation 心电信号的非线性分析 相关维数， hrv (heart rate variability)， t波交替（twa）， 心电信号的空间信息， qt离散度，,biomedical signal processing,32,nankai university, cy li, 2019/11/1,qrs波的检测,qrs的特点: 其能量在心电信号中占很大的比例, 其频谱分布在中高频区, 峰值落在10-20hz之间, 二阶导数算法 心电信号的一阶和二阶导数的平方和作为qrs波标记的脉冲信号, 移动平均算法 其求导平方运算和上相同, 并对求导平方数据进行移动平均, 从而突出qrs波的特征信息, 正交滤波算法,biomedical signal processing,33,nankai university, cy li, 2019/11/1,qrs波的频谱,qrs的特点: 其能量在心电信号中占很大的比例, 其频谱分布在中高频区, 峰值落在10-20hz之间,biomedical signal processing,34,nankai university, cy li, 2019/11/1,r波峰点的检测,斜率变号 双边阈值检测法 取一个固定的阈值ra, t1,t2分别为r波上升和下降通过这个阈值的时刻，则r基准点的位置t=(t1+t2)/2 固定宽度检测法 选一个固定宽度, 寻找为一个r波的t1, t2, t=(t1+t2)/2, 此法不受波形幅度变动和基线漂移的影响。,biomedical signal processing,35,nankai university, cy li, 2019/11/1,r波峰点的检测,dx/dt0,dx/dt0,ra,t2,t1,t=(t1+t2)/2,t1,t2,r,r,r,biomedical signal processing,36,nankai university, cy li, 2019/11/1,t波的检测,t波的检出 幅度和形状， 小波变换, t波的终点 零点， 最负的斜率切线和基线的交点， 最小二乘法拟合波, 用最小二乘法在t波后半支写率最大点附近的区域拟合一直线, 它于等电位的交点,biomedical signal processing,37,nankai university, cy li, 2019/11/1,心电图的小波变换,biomedical signal processing,38,nankai university, cy li, 2019/11/1,t波的检出四种方法,技术阈值法(th): t波与阈值水平的交点, 微分阈值(dth): t波的微分与阈值水平的交点, 技术斜率交点(si): t波最大斜率与等位线交点, 技术峰斜率交点(psi): t波高峰和t波最大斜率的脸限于等位线的交点, threshold: 5-15% of the peak,biomedical signal processing,39,nankai university, cy li, 2019/11/1,qt间隔的确定,qt duration: t(q)-t(t) ecg waveforms,biomedical signal processing,40,nankai university, cy li, 2019/11/1,心电频谱分析,fft ar model,biomedical signal processing,41,nankai university, cy li, 2019/11/1,频率心电图 fcg (frequency-domain correlative cardiograph),设为一个系统, 输入: v5, 输出: ii, 几个参数: 自功率谱, px, py, 自相关函数, vxx, vyy, 同一信号在时间上前后的相关性, 互相关函数, vxy, 反映两个导联的信号在不同时刻的相互关系, 传递函数相移, v5与ii信号间的相频对应关系,biomedical signal processing,42,nankai university, cy li, 2019/11/1,心电信号的时间特性,hrv是指心率的波动，或者心动周期之间的差异性。 正常人的心率有相当于平均值10%的波动，这种波动程度的降低是心脏异常的表现。 hrv的分析方法： 时域分析，rr间期的标准差，一般为100-140ms, 50ms异常, 频域分析，lf(0.04-0.15hz)和hf(0.15-0.4hz)的功率谱的比较， 非线性分析：散点图，相关维数，复杂度等。,biomedical signal processing,43,nankai university, cy li, 2019/11/1,心室晚电位(vlp),定义: 心室肌损伤后使局部心肌延迟除极引起的破裂电位, 位置：出现在qrs波末端和st段上的高频、低幅的微小电活动。 疾病：vlp常见于心肌缺血所致的室性心动过速，尤其是心肌更死后的室速，被认为是预测室性心律失常的信号。,biomedical signal processing,44,nankai university, cy li, 2019/11/1,心室晚电位(vlp),三个指标： qrs终末40ms 的电压: rms40=40ms, qrs 时限: qrs总时限110-119ms, vlp的检测, 微伏级(2-20 v), 信号平均技术, 叠加100-200次, 采用正交心电图x, y, z双极导联, rms:,biomedical signal processing,45,nankai university, cy li, 2019/11/1,twa (t wave alternans),t波交替 t 波交替现象是一种心电变异性现象, 指心电信号t 波段出现幅度、形状甚至极性的逐拍交替变化。 临床研究表明, v量级的twa 现象与室性心动过速和室颤易发性有关。检测体表心电图中是否存在微伏量级的twa 现象, 也成为临床上预测心脏性猝死的一个可靠而无创的指标。 目前tw