CN111067508B - 非临床环境下对高血压进行非干预式的监测和评估方法 - Google Patents
非临床环境下对高血压进行非干预式的监测和评估方法 Download PDFInfo
- Publication number
- CN111067508B CN111067508B CN201911412103.8A CN201911412103A CN111067508B CN 111067508 B CN111067508 B CN 111067508B CN 201911412103 A CN201911412103 A CN 201911412103A CN 111067508 B CN111067508 B CN 111067508B
- Authority
- CN
- China
- Prior art keywords
- features
- feature
- correlation coefficient
- time window
- hypertension
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/352—Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/363—Detecting tachycardia or bradycardia
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明提供一种非临床环境下对高血压进行非干预式的监测和评估方法,包括对采集的心电图数据进行分片,得到以5分钟的时间窗为单位的心电信号序列,利用交叠滑动窗口算法定位每次心跳所发生的时刻,从而计算初级心跳间期序列;采用阈值法对获取的初级心跳间期序列进行修正处理;从每个时间窗的心跳间期序列中提取线性域、非线性域和波形分布形态特征,通过皮尔逊相关性分析方法定量计算心率变异性线性域和非线性域特征之间的相关性,并基于特征之间的相关强度进行特征聚合,通过信息增益选择特征形成特征子集;融合特征子集特征,使用随机森林方法构建高血压评估模型。本发明进行了更加细粒度的分析,对高血压患者的评估准确率可以达到97.1%。
Description
技术领域
本发明属于高血压评价领域,具体涉及一种非临床环境下对高血压进行非干预式的监测和评估方法。
背景技术
高血压是指以体循环动脉血压(收缩压和/或舒张压)增高为主要特征,可伴有心、脑、肾等器官的功能或器质性损害的临床综合征,是一种典型的慢性病。高血压的患病人数在持续上升,1975年全世界约有5.94亿的高血压患者,而到了2015年,这一数字则超过了11亿,这就意味着,自1975年以来,30年的时间中全世界患高血压的成年人数量翻了近一倍。在我国,高血压同样是威胁人们健康的主要疾病之一,我国分别在1959年、1979年、1991年、2002年和2015年共进行过五次高血压患病率的调查,调查结果分别为5.10%,7.73%,11.26%,18.80%和23%。可以看出我国高血压的患病率呈快速上升趋势。
高血压因高患病率和高危害性一直是国内外重点防治的慢性病,是全社会普遍关注的重大公共卫生问题。高血压虽为慢性病,但危害性很大。短期来看,血压升高会引起头痛、头晕、失眠、记忆力衰退、注意力不集中等症状,直接影响患者的正常工作、生活;长期来看,持续的血压升高会对人的大脑、肝、肾脏等重要器官造成损伤,甚至引发脑卒中、心肌梗死、心力衰竭、肾功能衰竭等危及生命的临床并发症,是造成中风和心脏病的主要危险因素。此外,高血压致死率高,每年全世界约有750万人因高血压而死亡,其中我国因血压升高而过早死亡的达200万人,直接医疗费用约为366亿元,给家庭、社会带来极大的经济负担。
心脏的每一次收缩和舒张的过程都会产生一定的电位变化,通过心脏周围的导电组织和体液,反应到体表。心电图技术就是一种以时间为单位来记录心脏跳动过程中带来的电位变化的方法,通过放置在体表特定位置的电极来捕获和记录这种电位的变化。
人体的自主神经***能无意识地调节身体机能,如血压、心率、呼吸等,按照功能可分为交感神经和副交感神经。交感神经能使心跳加快,血压上升,当人处于紧张状态时,交感神经起主导作用,副交感神经能使心跳减慢,血压降低。高血压的形成有多种机制,交感神经兴奋是重要一环,高血压和交感神经的兴奋互为因果关系。心率变异性(HRV)能够反映自主神经***活性,是定量评估交感神经与副交感神经张力及其平衡性的生理指标。因此,可以使用HRV作为反映心血管疾病病情的一个重要指标,用来区分健康人和高血压患者。
目前,主要通过测量血压来判别是否患有高血压,一种是使用血压计或袖带测量,这种方法虽然测量结果比较准确,但无法持续测量血压值,无法监测血压的变化动态。另一种是使用可穿戴设备中的各种传感器测量脉搏波的传导时间间接计算血压值,这种方法能持续监测,但准确率较低。综上所述,早期无明显症状是高血压知晓率低的主要原因,需要一种在非临床环境下对高血压进行非干预式的监测和评估方法。
发明内容
本发明的目的是针对上述存在的问题和不足,提供一种非临床环境下基于非干预的,持续监测的心电信号实现高血压监测和评估的方法。
本发明采用的技术方案是:
非临床环境下对高血压进行非干预式的监测和评估方法,其步骤包括:
步骤S1:获取受试者的心电图数据:
S11:数据采集终端采集受试者夜间睡眠状态下的心电图数据,所述心电图数据采集时长范围为6-12h;
S12:细粒度分片:采用固定时间窗分片方法将采集的心电图数据切分成若干个时间窗,每个时间窗包括若干个心电波;
步骤S2:获取心跳间期序列:
S21:利用交叠滑动窗口算法定位每次心跳所发生的时刻,从而计算初级心跳间期序列;
S22:采用阈值法对获取的初级心跳间期序列进行修正处理,以获得最终的心跳间期序列;
步骤S3:多维心率变异特征的提取及数据分析:
S31:多维心率变异特征的提取:从每个时间窗的心跳间期序列中提取心率变异性特征,组合成原特征集,所述原特征集包括线性域特征、非线性域特征和分布形态特征;
S32:对线性域特征和非线性域特征进行相关性分析:
S321:将原特征集包括的所有线性域特征和非线性域特征排列成特征矩阵,采用皮尔逊相关性分析方法得到相关系数矩阵;
S322:基于线性域特征和非线性域特征之间的相关强度对提取的若干个线性域特征、若干个非线性域特征进行聚合,得到多组聚合后特征,采用信息增益对每组聚合后特征进行特征选择,除去冗余特征,形成新的特征子集;
S33:所述分布形态特征包括偏度、峰度及高斯混合模型的六个参数,对分布形态特征进行定量评价:
步骤S4:高血压评估模型构建:将上述步骤S3提取的特征子集进行融合,通过随机森林分类算法构建高血压评估模型;
步骤S5:通过特征子集的包含特征指标评价上述高血压评估模型的性能,验证评估模型的有效性,该特征指标包括方差、心跳间期标准差、差值均方根、间期大于50ms的比例、超低频、低频、高频、总功率、长轴标准差、短轴标准差、采样熵、近似熵。
优选的,在步骤S1中的1),所述心电图数据采集时长为8小时;
在步骤S12中,通过固定时间窗分片的方法将采集的心电图数据切分成96个时间窗,每个时间窗的时长为5min。
优选的,在步骤S31中,高血压患者心电图数据的采样频率是128HZ,则滑动窗口为128个采样点;健康人心电图数据的采样频率是200HZ,则滑动窗口为200个采样点,所述滑动窗口每次向前移动自身长度的一半。
优选的,步骤S321中,计算每个人的相关系数矩阵的具体方法为:
在特征提取部分,每个时间窗提取的线性域特征和非线性域特征共计16个特征,每个人8小时的数据划分为96个时间窗,从每个时间窗中提取同一维特征,组成特征列向量Fi,即Fi={Fi,1;Fi,2;...;Fi,j...Fi,96},其中Fi,j表示从第j个时间窗中的第i个特征,一共有16个特征,得到96行16列的特征矩阵FMk
其中,k表示的是第k个人,特征矩阵的行表示每个时间窗口中的特征,列表示每一维特征;
分别对特征矩阵FMk中的每两个特征向量计算皮尔逊相关系数,得到16行16列的相关系数矩阵CMk;
其中,PCCi,j表示表示第i,j这两个特征向量之间的相关系数;
皮尔逊相关系数,对于两个特征向量X和Y,X=(X1,X2,...,Xn),F2=(Y1,Y2,...,Yn),这两个特征的相关系数PCC的计算公式如下所示:
优选的,步骤S321中,计算每一类人的相关系数矩阵的具体方法为:
a.计算每个人的相关系数矩阵:
在特征提取部分,每个时间窗提取的线性域特征和非线性域特征共计16个特征,每个人8小时的数据划分为96个时间窗,从每个时间窗中提取同一维特征,组成特征列向量Fi,即Fi={Fi,1;Fi,2;...;Fi,j...Fi,96},其中Fi,j表示从第j个时间窗中的第i个特征,一共有16个特征,得到96行16列的特征矩阵FMk
其中,k表示的是第k个人,特征矩阵的行表示每个时间窗口中的特征,列表示每一维特征;
分别对特征矩阵FMk中的每两个特征向量计算皮尔逊相关系数,得到16行16列的相关系数矩阵CMk;
其中,PCCi,j表示表示第i,j这两个特征向量之间的相关系数;
皮尔逊相关系数,对于两个特征向量X和Y,X=(X1,X2,...,Xn),F2=(Y1,Y2,...,Yn),这两个特征的相关系数PCC的计算公式如下所示:
b.计算每类人内所有人的相关系数矩阵中的每一个元素的均值,得到该一类人的相关系数矩阵。
优选的,步骤S322中,选择特征的具体方法为:计算线性域特征和非线性域特征之间的相关强度,当线性域特征和非线性域特征之间的相关系数大于预设的相关强度阈值,其强度阈值范围值为(-1,+1),则判断该线性域特征和非线性域特征强相关,并将该线性域特征和非线性域特征分为一类,如此迭代,得到多组聚合后特征,并计算每组聚合后特征的信息增益,将每组中信息增益最大的特征代替该组的其他特征,完成特征选择。
优选的,在步骤S322中,将每个人所有时间窗单个特征的值进行聚合:由于所述时间窗为96个,则每个特征对应96个数值,每个特征对应的96个数值范围相差大,通过min-max归一化方法将每个特征的数值范围变换到(0,1)区间;当需要聚合为一个值时,则采用平方和求根的方法进行特征聚合。
优选的,在步骤S33中,由于每个时间窗内心跳间期序列的分布呈现双峰分布,运用高斯混合模型对每个时间窗心跳间期的分布进行建模的具体方法:
对于给定的心跳间期序列RR={x1,x2,…,xn},该心跳间期序列对应的高斯混合模型的概率密度表达式为
其中,fi(xiθi)表示第i个高斯分布的密度函数;未知参数θi=(μi,σi),μi和σi分别表示第i个高斯分布的均值和标准差,wi表示每个高斯分布在高斯混合模型中所占的权重系数;满足的约束条件为为了求解高斯混合模型中的未知参数wi,μi,σi采用期望最大化算法进行求解。
优选的,采用期望最大化算法求解高斯混合模型的未知参数wi,μi,σi具体包括:
高斯混合模型中的参数wi,μi,σi可以通过给定样本来确定,将高斯混合模型中的位置参数标记为ζ={ζ1,ζ2,L,ζK},对于ζi={wi,μi,σi},则表示混合高斯模型中第i个高斯分布的未知参数,分别对应每个高斯分布的权重、均值和标准差,为了得到未知参数的最优解,采用最大似然估计的方法来求解,具体为:先求样本的联合概率密度函数,因各个样本之间符合独立同分布,因此,这些样本的联合概率密度就是单个样本概率密度的乘积,表示为:
接着,求联合概率密度函数取对数,如下公式所示:
然后使用最大似然估计确定未知参数θi={μi,σi}的最优值,即:
具体步骤:
第一步,设定高斯混合模型中单个高斯分布的数量K,给每个高斯分布的未知参数wi,μi,σi赋初始值;
第二步,根据参数初值或上次迭代的模型参数计算隐变量w的后验概率,作为隐变量当前的估计值;
第三步,用上述计算得到的隐变量的值输入,将似然函数最大化计算新的参数wi,μi,σi的值;
第四步,反复迭代第二步和第三步,直至wi,μi,σi的值保持不变,这时就能得到参数wi,μi,σi的最大值。
与现有技术相比,本发明使用心电数据提取多维心率异性特征进行高血压的评估,利用心电数据采集时非干预性和可持续监测等优势,通过采集夜晚睡眠状态下的心电信号,并采用固定时间窗分片方法将原始心电信号切分成细粒度的时间片,使用滑动窗口方法检测R波,并通过R波错检和漏检算法对R波进行修正,得到心跳间期序列,进而提取心率变异性的线性域特征和非线性域特征及分布形态特征,并通过高斯混合模型对每个时间窗内心跳间期序列的分布进行建模。分析线性域特征和非线性域特征之间的相关强度,并进行特征选择,除去冗余特征,将提取的多维心率变异性特征进行融合,通过随机森林分类算法进行高血压评估模型。本发明进行了更加细粒度的分析,对高血压患者的评估准确率可以达到97.1%。
附图说明
图1为本发明的结构图;
图2为本发明中心跳间期提取过程的流程图;
图3为本发明进行特征选择的过程框图;
图4为本发明的流程框图;
图5为本发明中偏度示意图;
图6为本发明中高血压识别结果ROC曲线。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚完整的描述。显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域的普通技术人员在没有做出创造性劳动的前提下所获得的其他实施例,都属于本发明的保护范围。
参照图1-4所示,本发明具体公开了一种非临床环境下对高血压进行非干预式的监测和评估方法,其步骤包括:
步骤S1:获取受试者的心电图数据:
S11:数据采集终端采集受试者夜间睡眠状态下的心电图数据,所述心电图数据采集时长优选8小时;
S12:细粒度分片:采用固定时间窗分片的方法将采集的心电图数据切分成96个时间窗,每个时间窗的时长为5min,所述每个时间窗包括若干个心电波;
步骤S2:获取心跳间期序列:
心跳间期是指相邻两次心跳的时间间隔。血压的变化由交感神经与副交感神经共同控制,而心跳间期序列的变化被公认为是交感神经与副交感神经相互作用的最直观表现。因此分析心跳间期序列可以间接地反推某受试者是否患有高血压疾病,这也是本发明的理论基础。
本发明实施例采用以下方法进行心跳间期序列的获取:
S21:利用交叠滑动窗口算法定位每次心跳所发生的时刻,从而计算初级心跳间期序列;
具体地:
滑动窗口采用矩形窗,高血压患者心电图数据的采样频率是128HZ,则滑动窗口为128个采样点;健康人心电图数据的采样频率是200HZ,则滑动窗口为200个采样点,滑动窗口每次向前移动自身长度的一半,从而保证不遗漏任何可能的心跳信息。
每个窗口中最大值所对应的时刻即为该窗口内心跳发生的时刻,并将第i个窗口所对应的时刻记为Ti,则心跳时刻序列T可表示为:T=T1,T2,…Ti…Tn;其中,窗口的个数为n。若将第i个心跳间期记为RRi,则初级心跳间期序列RR可表示为:RR=RR1,RR2,…RRi,…RRn-1;即:RR=T2-T1,T3-T2,…Ti+1-Ti…Tn-Tn-1;
S22:采用阈值法对获取的初级心跳间期序列进行修正处理,以获得最终的心跳间期序列;
具体地:在对心跳间期进行统计分析的基础上,该阈值条件设定为:其中,表示初级心跳间期序列的平均值。不满足该条件的心跳间期使用平均值替代。
步骤S3:多维心率变异特征的提取及数据分析:
S31:多维心率变异特征的提取:从每个时间窗的心跳间期序列中提取出心率变异性特征,组合成原特征集,所述原特征集包括线性域特征、非线性域特征和分布形态特征,所述线性域特征包括时域特征和频域特征,其中,时域特征包括均值(Mean)、最大值(Max)、最小值(Min)、标准差(SD)、差值均方根(RMSSD)、差值标准差(SDSD)、相邻RR间期大于50ms的数量(NN50)和变异系数(VC);频域特征包括高频频段能量值(HF)、低频频段能量值(LF)、低频频段与高频频段能量比值(LF/HF)、归一化的高频频带能量(HFnorm)、归一化的低频频段能量(LFnorm);所述非线性域特征包括散点图的特征,分别是长轴的标准差(SD1)、短轴的标准差(SD2)及长轴短轴标准差的比值(SD1/SD2);所述分布形态特征包括偏度(SKEWNESS)、峰度(KURTOSIS)及高斯混合模型的六个参数;
S32:对线性域特征和非线性域特征进行相关性分析:
S321:
计算每个人的相关系数矩阵的具体方法为:
在特征提取部分,每个时间窗提取的线性域特征和非线性域特征共计16个特征,每个人8小时的数据划分为96个时间窗,从每个时间窗中提取同一维特征,组成特征列向量Fi,即Fi={Fi,1;Fi,2;...;Fi,j...Fi,96},其中Fi,j表示从第j个时间窗中的第i个特征,一共有16个特征,得到96行16列的特征矩阵FMk
其中,k表示的是第k个人,特征矩阵的行表示每个时间窗口中的特征,列表示每一维特征;
分别对特征矩阵FMk中的每两个特征向量计算皮尔逊相关系数,得到16行16列的相关系数矩阵CMk;
其中,PCCi,j表示表示第i,j这两个特征向量之间的相关系数;
皮尔逊相关系数,对于两个特征向量X和Y,X=(X1,X2,...,Xn),F2=(Y1,Y2,...,Yn),这两个特征的相关系数PCC的计算公式如下所示:
或
计算每一类人的相关系数矩阵的具体方法为:
a.计算每个人的相关系数矩阵:
在特征提取部分,每个时间窗提取的线性域特征和非线性域特征共计16个特征,每个人8小时的数据划分为96个时间窗,从每个时间窗中提取同一维特征,组成特征列向量Fi,即Fi={Fi,1;Fi,2;...;Fi,j...Fi,96},其中Fi,j表示从第j个时间窗中的第i个特征,一共有16个特征,得到96行16列的特征矩阵FMk
其中,k表示的是第k个人,特征矩阵的行表示每个时间窗口中的特征,列表示每一维特征;
分别对特征矩阵FMk中的每两个特征向量计算皮尔逊相关系数,得到16行16列的相关系数矩阵CMk;
其中,PCCi,j表示表示第i,j这两个特征向量之间的相关系数;
皮尔逊相关系数,对于两个特征向量X和Y,X=(X1,X2,...,Xn),F2=(Y1,Y2,...,Yn),这两个特征的相关系数PCC的计算公式如下所示:
b.计算每类人内所有人的相关系数矩阵中的每一个元素的均值,得到该一类人的相关系数矩阵。
S322:计算线性域特征和非线性域特征之间的相关强度,当线性域特征和非线性域特征之间的相关系数大于预设的相关强度阈值,其强度阈值范围值为(-1,+1),则判断该线性域特征和非线性域特征强相关,并将该线性域特征和非线性域特征分为一类,如此迭代,得到多组聚合后特征,并计算每组聚合后特征的信息增益,将每组中信息增益最大的特征代替该组的其他特征,完成特征选择,形成新的特征子集;
S33:所述分布形态特征包括偏度、峰度及高斯混合模型的六个参数,对分布形态特征进行定量评价(峰度及高斯混合模型的六个参数的定量评价过程为公开现有技术,因此在这里不再一一赘述。);
偏度系数的计算为例:给定心跳间期序列RR={RR1,RR2,...,RRn},偏度系数s的计算公式为:
考虑到每个时间窗内心跳间期序列的分布呈现双峰分布,用高斯混合模型对每个时间窗内心跳间期的分布进行建模:
对于给定的心跳间期序列RR={x1,x2,…,xn},该心跳间期序列对应的高斯混合模型的概率密度表达式为
其中,fi(xi|θi)表示第i个高斯分布的密度函数;未知参数θi=(μi,σi),μi和σi分别表示第i个高斯分布的均值和标准差,wi表示每个高斯分布在高斯混合模型中所占的权重系数;满足的约束条件为为了求解高斯混合模型中的未知参数wi,μi,σi,采用期望最大化算法进行求解。
为了求解高斯混合模型中的未知参数wi,μi,σi采用期望最大化算法进行求解,具体包括:
高斯混合模型中的参数wi,μi,σi可以通过给定样本来确定,将高斯混合模型中的位置参数标记为ζ={ζ1,ζ2,L,ζK},对于ζi={wi,μi,σi},则表示混合高斯模型中第i个高斯分布的未知参数,分别对应每个高斯分布的权重、均值和标准差,为了得到未知参数的最优解,采用最大似然估计的方法来求解,具体为:先求样本的联合概率密度函数,因各个样本之间符合独立同分布,因此,这些样本的联合概率密度就是单个样本概率密度的乘积,表示为:
接着,求联合概率密度函数取对数,如下公式所示:
然后使用最大似然估计确定未知参数θi={μi,σi}的最优值,即:
具体步骤:
第一步,设定高斯混合模型中单个高斯分布的数量K,给每个高斯分布的未知参数wi,μi,σi赋初始值;
第二步,根据参数初值或上次迭代的模型参数计算隐变量w的后验概率,作为隐变量当前的估计值;
第三步,用上述计算得到的隐变量的值输入,将似然函数最大化计算新的参数wi,μi,σi的值;
步骤S4:高血压评估模型构建:将上述步骤S3提取的特征子集进行融合,通过随机森林分类算法构建高血压评估模型;
步骤S5:通过特征子集的包含特征指标评价上述高血压评估模型的性能,验证评估模型的有效性,该特征指标包括方差var、心跳间期标准差sdnn、差值均方根rmssd、间期大于50ms的比例pnn50、超低频VLF、低频LF、高频HF、总功率TP、标准差sd1、sd2、采样熵SampEn、近似熵ApEn。
在上述步骤S322中,将每个人所有时间窗单个特征的值进行聚合:由于所述时间窗为96个,则每个特征对应96个数值,每个特征对应的96个数值范围相差大,通过min-max归一化方法将每个特征的数值范围变换到(0,1)区间;当需要聚合为一个值时,则采用平方和求根的方法进行特征聚合。
受试者的测量评估过程:
受试者由医院选择,共138人,均为高血压患者,采集受试者整晚ECG数据,受试者基本信息如下:
研究对象实验数据统计信息
将采用本发明的评估方法与已公开发表的国际同类研究进行比较,其结果如图6所示,可知,采用本发明的评估方法,其确认率远远大于其他的评估方法。
上述实施例仅用于说明本发明的技术方案,而非对其限制。本发明并不局限于上面已经描述并在附图中示意的准确结构,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,做出的各种改变和变形,都应当视为属于本发明的保护范围。
Claims (6)
1.非临床环境下对高血压进行非干预式的监测和评估方法,其特征在于,其步骤包括:
步骤S1:获取受试者的心电图数据:
S11:数据采集终端采集受试者夜间睡眠状态下的心电图数据,所述心电图数据采集时长范围为6-12h;
S12:细粒度分片:采用固定时间窗分片方法将采集的心电图数据切分成若干个时间窗,每个时间窗包括若干个心电波;
步骤S2:获取心跳间期序列:
S21:利用交叠滑动窗口算法定位每次心跳所发生的时刻,从而计算初级心跳间期序列;
S22:采用阈值法对获取的初级心跳间期序列进行修正处理,以获得最终的心跳间期序列;
步骤S3:多维心率变异特征的提取及数据分析:
S31:多维心率变异特征的提取:从每个时间窗的心跳间期序列中提取心率变异性特征,组合成原特征集,所述原特征集包括线性域特征、非线性域特征和分布形态特征;
S32:对线性域特征和非线性域特征进行相关性分析:
S321:将原特征集包括的所有线性域特征和非线性域特征排列成特征矩阵,采用皮尔逊相关性分析方法得到相关系数矩阵;
S322:基于线性域特征和非线性域特征之间的相关强度对提取的若干个线性域特征、若干个非线性域特征进行聚合,得到多组聚合后特征,采用信息增益对每组聚合后特征进行特征选择,除去冗余特征,形成新的特征子集;
S33:所述分布形态特征包括偏度、峰度及高斯混合模型的六个参数,对分布形态特征进行定量评价:
步骤S4:高血压评估模型构建:将上述步骤S3提取的特征子集进行融合,通过随机森林分类算法构建高血压评估模型;
步骤S5:通过特征子集的包含特征指标评价上述高血压评估模型的性能,验证评估模型的有效性,该特征指标包括方差、心跳间期标准差、差值均方根、间期大于50ms的比例、超低频、低频、高频、总功率、长轴标准差、短轴标准差、采样熵、近似熵;
步骤S321中,计算每个人的相关系数矩阵的具体方法为:
在特征提取部分,每个时间窗提取的线性域特征和非线性域特征共计16个特征,每个人8小时的数据划分为96个时间窗,从每个时间窗中提取同一维特征,组成特征列向量Fi,即Fi={Fi,1;Fi,2;...;Fi,j...Fi,96},其中Fi,j表示从第j个时间窗中的第i个特征,一共有16个特征,得到96行16列的特征矩阵FMk
分别对特征矩阵FMk中的每两个特征向量计算皮尔逊相关系数,得到16行16列的相关系数矩阵CMk;
其中,PCCi,j表示表示第i,j这两个特征向量之间的相关系数;
皮尔逊相关系数,对于两个特征向量X和Y,X=(X1,X2,...,Xn),F2=(Y1,Y2,...,Yn),这两个特征的相关系数PCC的计算公式如下所示:
步骤S321中,计算每一类人的相关系数矩阵的具体方法为:
a.计算每个人的相关系数矩阵:
在特征提取部分,每个时间窗提取的线性域特征和非线性域特征共计16个特征,每个人8小时的数据划分为96个时间窗,从每个时间窗中提取同一维特征,组成特征列向量Fi,即Fi={Fi,1;Fi,2;...;Fi,j...Fi,96},其中Fi,j表示从第j个时间窗中的第i个特征,一共有16个特征,得到96行16列的特征矩阵FMk
分别对特征矩阵FMk中的每两个特征向量计算皮尔逊相关系数,得到16行16列的相关系数矩阵CMk;
其中,PCCi,j表示表示第i,j这两个特征向量之间的相关系数;
皮尔逊相关系数,对于两个特征向量X和Y,X=(X1,X2,...,Xn),F2=(Y1,Y2,...,Yn),这两个特征的相关系数PCC的计算公式如下所示:
b.计算每类人内所有人的相关系数矩阵中的每一个元素的均值,得到该一类人的相关系数矩阵;
步骤S322中,选择特征的具体方法为:计算线性域特征和非线性域特征之间的相关强度,当线性域特征和非线性域特征之间的相关系数大于预设的相关强度阈值,其强度阈值范围值为(-1,+1),则判断该线性域特征和非线性域特征强相关,并将该线性域特征和非线性域特征分为一类,如此迭代,得到多组聚合后特征,并计算每组聚合后特征的信息增益,将每组中信息增益最大的特征代替该组的其他特征,完成特征选择。
2.根据权利要求1所述的非临床环境下对高血压进行非干预式的监测和评估方法,其特征在于,在步骤S11中,所述心电图数据采集时长为8小时;
在步骤S12中,通过固定时间窗分片的方法将采集的心电图数据切分成96个时间窗,每个时间窗的时长为5min。
3.根据权利要求1所述的非临床环境下对高血压进行非干预式的监测和评估方法,其特征在于,在步骤S31中,高血压患者心电图数据的采样频率是128HZ,则滑动窗口为128个采样点;健康人心电图数据的采样频率是200HZ,则滑动窗口为200个采样点,所述滑动窗口每次向前移动自身长度的一半。
4.根据权利要求3所述的非临床环境下对高血压进行非干预式的监测和评估方法,其特征在于,在步骤S322中,将每个人所有时间窗单个特征的值进行聚合:由于所述时间窗为96个,则每个特征对应96个数值,每个特征对应的96个数值范围相差大,通过min-max归一化方法将每个特征的数值范围变换到(0,1)区间;当需要聚合为一个值时,则采用平方和求根的方法进行特征聚合。
6.根据权利要求5所述的非临床环境下对高血压进行非干预式的监测和评估方法,其特征在于,采用期望最大化算法求解高斯混合模型的未知参数wi,μi,σi,具体包括:
高斯混合模型中的参数wi,μi,σi可以通过给定样本来确定,将高斯混合模型中的位置参数标记为ζ={ζ1,ζ2,L,ζK},对于ζi={wi,μi,σi},则表示混合高斯模型中第i个高斯分布的未知参数,分别对应每个高斯分布的权重、均值和标准差,为了得到未知参数的最优解,采用最大似然估计的方法来求解,具体为:先求样本的联合概率密度函数,因各个样本之间符合独立同分布,因此,这些样本的联合概率密度就是单个样本概率密度的乘积,表示为:
接着,求联合概率密度函数取对数,如下公式所示:
然后使用最大似然估计确定未知参数θi={μi,σi}的最优值,即:
具体步骤:
第一步,设定高斯混合模型中单个高斯分布的数量K,给每个高斯分布的未知参数wi,μi,σi赋初始值;
第二步,根据参数初值或上次迭代的模型参数计算隐变量w的后验概率,作为隐变量当前的估计值;
第三步,用上述计算得到的隐变量的值输入,将似然函数最大化计算新的参数wi,μi,σi的值;
第四步,反复迭代第二步和第三步,直至wi,μi,σi的值保持不变,这时就能得到参数wi,μi,σi的最大值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911412103.8A CN111067508B (zh) | 2019-12-31 | 2019-12-31 | 非临床环境下对高血压进行非干预式的监测和评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911412103.8A CN111067508B (zh) | 2019-12-31 | 2019-12-31 | 非临床环境下对高血压进行非干预式的监测和评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111067508A CN111067508A (zh) | 2020-04-28 |
CN111067508B true CN111067508B (zh) | 2022-09-27 |
Family
ID=70321158
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911412103.8A Active CN111067508B (zh) | 2019-12-31 | 2019-12-31 | 非临床环境下对高血压进行非干预式的监测和评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111067508B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111820879A (zh) * | 2020-06-23 | 2020-10-27 | 浙江清华长三角研究院 | 一种适用于慢性疾病患者的健康评价管理方法 |
CN112370036A (zh) * | 2020-10-20 | 2021-02-19 | 复旦大学 | 基于级联型rls自适应滤波的ppg心率提取装置和方法 |
CN113017650B (zh) * | 2021-03-12 | 2022-06-28 | 南昌航空大学 | 一种基于功率谱密度图像的脑电特征提取方法和*** |
CN115363598B (zh) * | 2022-09-19 | 2023-03-24 | 深圳市太美亚电子科技有限公司 | 一种心电信号处理方法及*** |
CN116821809B (zh) * | 2023-08-30 | 2023-12-12 | 济宁医学院附属医院 | 基于人工智能的生命体征数据采集*** |
CN116861369A (zh) * | 2023-09-04 | 2023-10-10 | 山东捷瑞数字科技股份有限公司 | 工业互联网异构数据融合的空压机健康评估方法及*** |
CN117591987B (zh) * | 2024-01-18 | 2024-04-26 | 北京国旺盛源智能终端科技有限公司 | 一种基于人工智能的用电设备监测方法及*** |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1682654A (zh) * | 2004-04-15 | 2005-10-19 | Ge医药***信息科技公司 | 用于确定心电图信号的交替数据的方法和装置 |
CN104382582A (zh) * | 2014-11-10 | 2015-03-04 | 哈尔滨医科大学 | 一种对动态心电数据进行分类的装置 |
CN105852807A (zh) * | 2016-03-25 | 2016-08-17 | 中山大学 | 一种无监督式自主神经功能量化评估方法 |
CN106570178A (zh) * | 2016-11-10 | 2017-04-19 | 重庆邮电大学 | 一种基于图聚类的高维文本数据特征选择方法 |
CN108491769A (zh) * | 2018-03-08 | 2018-09-04 | 四川大学 | 基于rr间期和多特征值的房颤分类方法 |
CN108577804A (zh) * | 2018-02-02 | 2018-09-28 | 西北工业大学 | 一种面向高血压病人监测的bcg信号分析方法及*** |
CN109117769A (zh) * | 2018-07-31 | 2019-01-01 | 东南大学 | 一种针对穿戴式心电采集的实时质量评估反馈方法 |
CN110151169A (zh) * | 2019-07-04 | 2019-08-23 | 中山大学 | 一种基于心电数据的睡眠状态识别分类方法 |
CN110393519A (zh) * | 2019-08-19 | 2019-11-01 | 广州视源电子科技股份有限公司 | 心电信号的分析方法、装置、存储介质和处理器 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007074600A1 (ja) * | 2005-12-26 | 2007-07-05 | Nec Corporation | 特徴抽出装置、特徴抽出方法、および特徴抽出プログラム |
JP5553303B2 (ja) * | 2010-02-18 | 2014-07-16 | 株式会社デルタツーリング | 生体状態推定装置及びコンピュータプログラム |
US12002545B2 (en) * | 2011-08-03 | 2024-06-04 | ExSano, Inc. | Technique for identifying features |
CN204515353U (zh) * | 2015-03-31 | 2015-07-29 | 深圳市长桑技术有限公司 | 一种智能手表 |
EP3153097B1 (en) * | 2015-10-07 | 2024-04-03 | Tata Consultancy Services Limited | A method and system for assessing learning experience of a person |
CN109044301A (zh) * | 2018-07-02 | 2018-12-21 | 西北工业大学 | 一种面向高血压疾病检测的bcg信号分析方法 |
CN109528197B (zh) * | 2018-11-20 | 2022-07-08 | 中国科学院脑科学与智能技术卓越创新中心 | 基于脑功能图谱进行精神疾病的个体化预测方法和*** |
CN109770859A (zh) * | 2019-03-28 | 2019-05-21 | 广州视源电子科技股份有限公司 | 心电信号的处理方法和装置、存储介质、处理器 |
CN109770861B (zh) * | 2019-03-29 | 2022-04-12 | 广州视源电子科技股份有限公司 | 心电节律模型的训练及检测方法、装置、设备和存储介质 |
CN110299205A (zh) * | 2019-07-23 | 2019-10-01 | 上海图灵医疗科技有限公司 | 基于人工智能的生物医学信号特征处理及评测方法、装置及应用 |
-
2019
- 2019-12-31 CN CN201911412103.8A patent/CN111067508B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1682654A (zh) * | 2004-04-15 | 2005-10-19 | Ge医药***信息科技公司 | 用于确定心电图信号的交替数据的方法和装置 |
CN104382582A (zh) * | 2014-11-10 | 2015-03-04 | 哈尔滨医科大学 | 一种对动态心电数据进行分类的装置 |
CN105852807A (zh) * | 2016-03-25 | 2016-08-17 | 中山大学 | 一种无监督式自主神经功能量化评估方法 |
CN106570178A (zh) * | 2016-11-10 | 2017-04-19 | 重庆邮电大学 | 一种基于图聚类的高维文本数据特征选择方法 |
CN108577804A (zh) * | 2018-02-02 | 2018-09-28 | 西北工业大学 | 一种面向高血压病人监测的bcg信号分析方法及*** |
CN108491769A (zh) * | 2018-03-08 | 2018-09-04 | 四川大学 | 基于rr间期和多特征值的房颤分类方法 |
CN109117769A (zh) * | 2018-07-31 | 2019-01-01 | 东南大学 | 一种针对穿戴式心电采集的实时质量评估反馈方法 |
CN110151169A (zh) * | 2019-07-04 | 2019-08-23 | 中山大学 | 一种基于心电数据的睡眠状态识别分类方法 |
CN110393519A (zh) * | 2019-08-19 | 2019-11-01 | 广州视源电子科技股份有限公司 | 心电信号的分析方法、装置、存储介质和处理器 |
Non-Patent Citations (1)
Title |
---|
Gaussian Mixture Model of Heart Rate Variability;Tommaso Costa等;《PLOS ONE》;20120530;第7卷(第5期);第1-9页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111067508A (zh) | 2020-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111067508B (zh) | 非临床环境下对高血压进行非干预式的监测和评估方法 | |
Acharya et al. | Comprehensive analysis of cardiac health using heart rate signals | |
Oweis et al. | QRS detection and heart rate variability analysis: A survey | |
Clifford et al. | ECG statistics, noise, artifacts, and missing data | |
Sansone et al. | Electrocardiogram pattern recognition and analysis based on artificial neural networks and support vector machines: a review | |
EP3454730B1 (en) | System for determining the frequency of a periodic physiological process of a subject | |
EP2286723A1 (en) | Method and apparatus for determining sleep apnea | |
CN105147248A (zh) | 基于生理信息的抑郁症评估***及其评估方法 | |
Acharya et al. | Heart rate variability | |
CN204931634U (zh) | 基于生理信息的抑郁症评估*** | |
US20210007621A1 (en) | Method to analyze cardiac rhythms using beat-to-beat display plots | |
Rahimi et al. | Sleep stage classification based on ECG-derived respiration and heart rate variability of single-lead ECG signal | |
CN109938719A (zh) | 一种基于生理参数的驾驶员疲劳检测方法 | |
Ihsan et al. | Study of feature extraction algorithms on photoplethysmography (ppg) signals to detect coronary heart disease | |
Aziz et al. | Classification of heart rate signals of healthy and pathological subjects using threshold based symbolic entropy | |
Chou et al. | Comparison between heart rate variability and pulse rate variability for bradycardia and tachycardia subjects | |
Fan et al. | An electrocardiogram acquisition and analysis system for detection of human stress | |
Mendez et al. | Automatic detection of sleep macrostructure based on bed sensors | |
Hikmah et al. | Sleep Quality Assessment from Robust Heart and Muscle Fatigue Estimation Using Supervised Machine Learning. | |
CN114886403A (zh) | 一种基于脉搏主波间期的恶性心律失常识别与预测*** | |
Sisodia et al. | Sleep order detection model using support vector machines and features extracted from brain ECG signals | |
Bassiouni et al. | Combination of ECG and PPG signals for smart healthcare systems: Techniques, applications, and challenges | |
Pale et al. | Heart rate variability analysis using different wavelet transformations | |
Kumar et al. | Statistical heart rate variability analysis under rest and post-exercise | |
Gnanavel et al. | GUI Base Prediction of Heart Stroke Stages by Finding the Accuracy using Machine Learning Algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |