一种脉象信号时域特征提取方法技术

本发明专利技术属于中医脉象测量技术领域，涉及一种脉象信号时域特征提取方法，包括：对采集到的脉象信号进行小波变换，去除高频噪声，并重建信号，得到去除了高频噪声的脉象信号；提取基线，并去除基线漂移，实现脉象信号的预处理；对经过预处理的脉象信号求取求一阶导数并找出每段脉象信号中的极大值，定位每个估计的脉象周期的上升支；对于每段脉象信号；确定脉象周期；寻找降中峡E和重搏波F。本发明专利技术能够提高脉象时域特征提取效率和准确性。

【技术实现步骤摘要】

本专利技术属于中医脉象测量
，涉及一种脉象信号特征提取方法。
技术介绍
目前，在脉象信号的时域特征提取方面主要应用的方法有阈值法和极值点法、平滑滤波和小波变换等。阈值法是指根据脉象波形的最大值与最小值确定脉象特征点的取值范围，再结合成对的极值点确定脉象波形的特征点。这种阈值法对于具有较大波动的脉象不能准确确定特征点的范围，同时利用极值点法的方法寻找特征点时，运算量比较大，脉象数据较多时， 不适宜采用此种方法。利用小波变换的系数法和过零点法判断脉象信号的特征点时，所选取的小波函数未必能适应所有脉象信号。小波变换的过零点法会产生多个特征点，不利于准确识别，应当寻求利用小波去噪以达到较好的特征提取效果。
技术实现思路
本专利技术的目的是克服现有技术的上述不足，提供一种适用于中医脉象信号特征提取的方法，以达到更好的特征提取效果。为此，本专利技术采用如下的技术方案。，包括下列步骤1)对采集到的脉象信号进行小波变换，去除高频噪声，并重建信号，得到去除了高频噪声的脉象信号；2)提取基线，并去除基线漂移，实现脉象信号的预处理；3)对经过预处理的脉象信号求取求一阶导数并找出每段脉象信号中的极大值，定位每个估计的脉象周期的上升支；4)对于每段脉象信号，以脉象上升支基点为起点，向前寻找占估计的脉象周期的 1/10的点数中的最小值，找到的第一个最小值即为脉象周期起始点；向后脉象周期的1/10 的点数中的最大值，找到的第一个最大值即为脉象周期主波峰。5)确定脉象周期；6)以每段脉象周期结束点为基点，向前在脉象周期的后2/3范围内依据脉象在此范围内的成对极值点来寻找降中峡E和重搏波F ；作为优选实施方式，步骤1)中，对采集到的脉象信号进行db8小波变换，分解层数为14，去除高频噪声，将低频系数与处理后的高频系数进行重建；步骤2)中，根据估计的脉象周期将脉象信号分段，取出每段最低点的横、纵坐标，对取出的数据做三阶曲线拟合，提取基线；步骤3)中，对于每段脉象信号，按照下列方法求取一阶导数极大值a.找出一阶导数的最大值，即为Imax ；b.设定阈值系数，并将Imax与阈值系数相乘，得到阈值Ith ；c.按照下列条件求取极大值点i :(l)i >= Ith ；(2) i同时分别大于左侧邻近的两个点和右侧邻近的2个点。本专利技术提出的脉象信号时域特征提取方法，结合了分段去基线和小波去噪对脉象信号进行预处理，改进了阈值法，并根据各特征点在脉象周期中的具体位置，提出在某一段时间范围内寻找该特征点，提高了特征提取效率和准确性。附图说明图1原始信号与去噪信号及其局部放大的对比图，(a)和(b)分别为原始信号及其局部信号，(c)和(d)分别为去噪信号及其局部信号。图2去噪信号与去噪去基线漂移信号的对比图，上图为去噪信号，下图为预处理后的脉象信号。图3上图为预处理后的脉象信号，下图为脉象信号的一阶导数。图4标记了起始点(*)与主波峰(ο)的脉象信号。图5.标记了降中峡(“ + ”)与重搏波(五角星)的脉象信号。图6.脉象信号的重搏前波(以正方形标记)。具体实施例方式在对本专利技术进行说明之前，首先介绍一下各种脉象参数生理含义①升支脉象形中由基线至主波峰顶的一条上升曲线，是心室的快速射血时期。②降支脉象形中由主波峰顶至基线的一条下降曲线，是心室射血后期至下一次心动周期的开始。③主波脉图的主体波幅，一般顶点为脉图的最高峰，反映动脉内压力与容积的最大值。④潮波又称重搏前波，位于下降支，主波之后，一般低于主波而高于重搏波，反映左心室停止射血，动脉扩张降压，逆向反射波。⑤降中峡或称降中波，是主波降支与重搏波升支构成的向下的切迹波谷，表示主动脉静压排空时间，为心脏收缩与舒张的分界点。⑥重搏波是降支中突出的一个上升波，为主动脉瓣关闭、主动脉弹性回缩波。下面介绍脉象测量中常用的各项参数和指标的生理含义及其测量方法(1)时间指标①t 脉动周期，即脉图起始点到终点的时值。t值对应于左心室的一个心动周期。测量方法测量U-U，点间隔的时间。②tl 升支时间，即脉图起点到主波峰点的时值。、值对应于左心室快速射血期。测定方法主波峰点到基线的垂线与基线的交点与U点之间的时值。③t4 心缩时间，即脉图起点到降中峡之间的时值。t4值对应于左心室的收缩期。测量方法降中峡切迹点到基线的垂线与基线的交点与U点之间的时值。④t5 缓降时间，即降中峡到脉图终点之间的时值。t5值对应于左心室舒张期。测量方法降中峡切迹点到基线的垂线与基线的交点与U’点之间的时值。⑤W 主波上1/3处的宽度，相当于动脉内高压力水平所维持的时间。测量方法先由主波峰顶至基线垂线上1/3与下2/3交界点作一水平线(平行于基线)再测量该水平线与主波两侧升降支相交两点间的时值。(2)波幅指标①、主波高度，即主波峰顶到脉图基线的垂线长度。Ii1值代表心收缩期动管壁承受的压力和容积，反映了左心室射血功能和大动脉的顺应性。测量方法测量主波峰顶至基线垂直距离。②Il3 潮波高度，即潮波峰顶到脉图基线的幅度。Ii3值主要反映动脉血管张力和外周阻力状态。测量方法潮波顶至基线的垂直距离。③比降中峡高度，即降中峡谷底到脉图基线的幅度，反映动脉外周阻力大小。测量方法测量降中峡最低点至基线间的垂直距离。④h5 重搏波高度，即重搏波峰顶点至过降中峡谷底水平基线的幅度，反映大动脉的弹性(顺应性)状况。测量方法测量过重搏波峰顶和过降中峡谷底所作的两条水平平行线之间的距1 O(3)角度指标用量角器测量，以“度”为单位。①α :上升角，或称U角，主波升支与基线的夹角，反映血管弹性与血液粘性。测量方法以U点为角顶点，量出升支起始段与基线夹角的度数。②θ :主波角，或称P角，是主波升支与降支的夹角，反映血管弹性和血流状况。测量方法将上升支与下降支速降段各引伸一条延长线，两条延长线相交之夹角即θ角。若θ角呈平项者，可不测量角度数值，而以文字记述为“平顶型角”。(4)面积指标脉图面积指标有5项。①Aa 收缩期面积，是收缩期脉图基线以上面积(单位mm2)。②Ab 舒张期面积，是舒张期脉图基线以上面积(单位mm2)。③AT 脉图总面积，是基线以上脉图总面积(单位mm2)。以上Aa、Ab和AT可通过积分求得，但在实际工作中，也可通过求积仪或梯形面积累加法求得。其中，AT、Ab是由脉图估价心功能某些指标的基本参数。④As 收缩期总面积，是动脉腔压力与静脉腔压力之间围成的收缩期总面积(单位 mmHg · s) ο⑤Ad 舒张期总面积，是动脉腔压力与静脉腔压力之间围成的舒张期总面积(单位 mmHg · s) οAs和Ad的计算方法分别是As = 0. 04Ad = 0. 04式中Ps、Pd分别表示左臂肱动脉收缩压值和舒张压值。(5)比值指标为了更好地反映心血管的功能状态和脉图特征，除以上所测绝对数值外，常取各参数的相对比值，在反映脉图的生理含义上更为准确和灵敏。①Vt1 升支斜率 ’②Vh1 张力系数；③Vh1 阻力系数；④Vh1 弹性系数。下面介绍本专利技术采用的预处理方法脉象信号是一种弱的生物电信号，易引入背景干扰。在采集时，脉象信号会受到人体呼吸、肌肉紧绷等原因而引起基线漂移，同时还会混入高频噪声，因此在对脉象信号进行特征点提取之前要对基线漂移进行处理，同时滤除或本文档来自技高网...

【技术保护点】
１．一种脉象信号时域特征提取方法，包括下列步骤：１）对采集到的脉象信号进行小波变换，去除高频噪声，并重建信号，得到去除了高频噪声的脉象信号；２）提取基线，并去除基线漂移，实现脉象信号的预处理；３）对经过预处理的脉象信号求取求一阶导数并找出每段脉象信号中的极大值，定位每个估计的脉象周期的上升支；４）对于每段脉象信号，以脉象上升支基点为起点，向前寻找占估计的脉象周期的１／１０的点数中的最小值，找到的第一个最小值即为脉象周期起始点；向后脉象周期的１／１０的点数中的最大值，找到的第一个最大值即为脉象周期主波峰。５）确定脉象周期；６）以每段脉象周期结束点为基点，向前在脉象周期的后２／３范围内依据脉象在此范围内的成对极值点来寻找降中峡Ｅ和重搏波Ｆ。

【技术特征摘要】

【专利技术属性】
技术研发人员：王联，周鹏，李想，
申请(专利权)人：芜湖圣美孚科技有限公司，
类型：发明
国别省市：34