CN113855007A - 一种用于血糖预测的机器学习模型样本获取的方法及装置 - Google Patents

一种用于血糖预测的机器学习模型样本获取的方法及装置 Download PDF

Info

Publication number
CN113855007A
CN113855007A CN202110995661.2A CN202110995661A CN113855007A CN 113855007 A CN113855007 A CN 113855007A CN 202110995661 A CN202110995661 A CN 202110995661A CN 113855007 A CN113855007 A CN 113855007A
Authority
CN
China
Prior art keywords
value
characteristic value
blood pressure
pressure
electrocardiosignals
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.)
Pending
Application number
CN202110995661.2A
Other languages
English (en)
Inventor
蔡其武
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Lianwei Medical Technology Shanghai Co ltd
Original Assignee
Lianwei Medical Technology Shanghai Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Lianwei Medical Technology Shanghai Co ltd filed Critical Lianwei Medical Technology Shanghai Co ltd
Priority to CN202110995661.2A priority Critical patent/CN113855007A/zh
Publication of CN113855007A publication Critical patent/CN113855007A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/14532Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring glucose, e.g. by tissue impedance measurement
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/021Measuring pressure in heart or blood vessels
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/352Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/353Detecting P-waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/355Detecting T-waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/358Detecting ST segments
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/366Detecting abnormal QRS complex, e.g. widening
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Cardiology (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Artificial Intelligence (AREA)
  • Physiology (AREA)
  • Emergency Medicine (AREA)
  • Vascular Medicine (AREA)
  • Optics & Photonics (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本申请通过一种用于血糖预测的机器学习模型样本获取的技术方案,获取心电信号及其对应的血糖值;确定所述心电信号的多维特征值;将所述心电信号的多维特征值输入预设的血压预测模型,获得血压预测值;将所述血糖值确定为所述心电信号的多维三特征值、所述血压预测值的标签,并将所述心电信号的多维特征值、所述血压预测值及所述标签确定为所述机器学习模型样本。进一步地,基于若干所述样本,构建所述机器学习模型的样本集,并基于所述样本集训练所述机器学习模型,以确定血糖预测模型。该方案只需测量心电信号即可预测血糖,提升了用户体验,而且基于多样性的样本集确定的血糖预测模型可用于对不同人的心电信号进行血糖预测,具有普遍适用性。

Description

一种用于血糖预测的机器学习模型样本获取的方法及装置
技术领域
本申请涉及心电信号处理技术领域,尤其涉及一种用于血糖预测的机器学习模型样本获取的技术。
背景技术
心脏在跳动时可以产生生物电信号,自心脏向身体的各个部位传导。由于身体各部分组织的结构不同,与心脏之间的距离也不同,使得心电信号在体表的不同部位,表现出不同的电位电压变化。
心电信号是一种周期性信号,具有频率、振幅(或功率)和相位等要素,易于检测,具有较直观的规律性,是人类最早研究并应用于医学临床的一种生物电信号。典型的心电信号波形如图1所示,x轴为以毫秒(mS)为单位的时间(T),y轴为以毫伏特(mV)为单位的电压(V),可被认为是由一系列不同频率、振幅和相位的正弦波精确叠加而成的。
随着计算机技术迅速发展,特别是随着人工智能技术的迅速发展,出现越来越多的采用机器学习技术对心电信号进行分析的方法,提取心电信号的特征,通过机器学习模型进行自动分析,可以获取与心电信号对应的人体生理指标,比如,可通过经过学习的机器学习模型来预测与心电信号对应的血糖值。
目前常用的血糖检测技术大都是有创的,而现有通过学习后的机器学习模型来预测与心电信号对应的血糖值的方法,只基于心电信号的个别参数进行数值处理、提取特征值,构建的机器学习模型需要针对每个人的血压单独建立模型函数,不具有普遍性。
发明内容
本申请的目的是提供一种用于血糖预测的机器学习模型样本获取的技术方案,用以构建适用于不同人的一种基于心电信号的多维特征值预测血压进而来预测血糖值的机器学习模型。
根据本申请的一个方面,提供了一种用于血糖预测的机器学习模型样本获取的方法,其中,所述方法包括:
获取心电信号及其对应的血糖值,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
将所述心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
将所述血糖值标识为所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值和所述血压舒张压预测值的标签,并将所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值、所述血压舒张压预测值及所述标签确定为所述机器学习模型样本。
可选地,其中,确定所述心电信号的第一特征值包括:
确定所述心电信号中每个心跳波形;
基于每个心跳波形,获取每个心跳波形的特征值,
遍历所述心电信号的全部心跳波形,获取全部心跳波形的特征值的统计值,并将所述统计值确定为所述心电信号的第一特征值。
可选地,其中,所述每个心跳波形的特征值包括以下至少一项:
每个心跳波形的P波时长;
每个心跳波形的PR间期;
每个心跳波形的QRS时间;
每个心跳波形的ST段时长;
每个心跳波形的T波时长;
每个心跳波形的TP段时长;
每个心跳波形的QT间期;
每个心跳波形的U波时长。
可选地,所述统计值包括以下至少一项:
均方根;
中值绝对偏差;
中值变异系数。
可选地,其中,确定所述心电信号的第二特征值包括:
确定所述心电信号中每个心跳波形;
基于两两相邻的两个心跳波形,获取两两相邻的两个心跳波形的特征值,
遍历所述心电信号的全部心跳波形,获取全部两两相邻的两个心跳波形的特征值的统计值,并将所述统计值确定为所述心电信号的第二特征值。
可选地,其中,所述两两相邻的两个心跳波形的特征值包括以下至少一项:
两两相邻的心跳波形的PP间期;
两两相邻的心跳波形的QQ间期;
两两相邻的心跳波形的RR间期;
两两相邻的心跳波形的SS间期;
两两相邻的心跳波形的TT间期;
两两相邻的心跳波形的UU间期。
可选地,其中,所述统计值包括以下至少一项:
均方根;
中值绝对偏差;
中值变异系数;
大于50毫秒的连续间隔差异数。
可选地,其中,确定所述心电信号的第三特征值包括:
获取所述第二特征值对应的时域波形的功率频谱;
基于所述功率频谱,获取若干预设频域的功率值;
将所述若干预设频域的功率值之间的比值确定为所述心电信号的第三特征值。
可选地,其中,确定所述心电信号的第三特征值还包括:
基于所述功率频谱,计算所述功率频谱的熵值;
将所述功率频谱的熵值确定为所述心电信号的第三特征值。
可选地,所述一种用于血糖预测的机器学习模型样本获取的方法还包括:
基于获取的若干心电信号及其对应的血糖值,确定若干所述机器学习模型样本,以构建所述机器学习模型的样本集。
可选地,所述方法还包括:
基于所述样本集,训练所述机器学习模型,当所述机器学习模型的评估指标符合预设阈值,将训练后的机器学习模型确定为血糖预测模型。
可选地,所述方法还包括:
获取待测对象的心电信号,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
将所述心电信号的第一特征值、第二特征值及第三特征值分别输入所述血压收缩压预测模型和所述血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
将所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值及所述血压舒张压预测值输入所述血糖预测模型,预测所述待测对象的血糖值。
根据本申请的又一方面,还提供了一种用于血糖预测的机器学习模型样本获取的装置,其中,所述装置包括:
第一模块,用于获取心电信号及其对应的血糖值,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
第二模块,用于将所述心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
第三模块,用于将所述血糖值标识为所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值和所述血压舒张压预测值的标签,并将所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值和所述血压舒张压预测值及所述标签确定为所述机器学习模型样本。
可选地,所述一种用于血糖预测的机器学习模型样本获取的装置还包括:
第四模块,用于基于获取的若干心电信号及其对应的血糖值,确定若干所述机器学习模型样本,以构建所述机器学习模型的样本集。
可选地,所述装置还包括:
第五模块,用于基于所述样本集,训练所述机器学习模型,当所述机器学习模型的评估指标符合预设阈值,将训练后的机器学习模型确定为血糖预测模型。
可选地,所述装置还包括:
第六模块,用于获取待测对象的心电信号,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
第七模块,用于将所述心电信号的第一特征值、第二特征值及第三特征值分别输入所述血压收缩压预测模型和所述血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
第八模块,用于将所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值及所述血压舒张压预测值输入所述血糖预测模型,预测所述待测对象的血糖值。
根据本申请的另一方面,还提供一种用于血糖预测的机器学习模型样本获取的方法,其中,所述方法包括:
获取采集对象的状态信息、心电信号及其对应的血糖值,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
将所述心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
将所述血糖值标识为所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值和所述血压舒张压预测值的标签,并将所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值、所述血压舒张压预测值及所述标签确定为所述机器学习模型样本。
可选地,所述方法还包括:
基于获取的若干采集对象的状态信息、心电信号及其对应的血糖值,确定若干所述机器学习模型样本,以构建所述机器学习模型的样本集。
可选地,其中,所述方法还包括:
基于所述样本集,训练所述机器学习模型,当所述机器学习模型的评估指标符合预设阈值,将训练后的机器学习模型确定为血糖预测模型。
可选地,其中,所述方法还包括:
获取待测对象的状态信息及心电信号,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
将所述心电信号的第一特征值、第二特征值及第三特征值分别输入所述血压收缩压预测模型和所述血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
将所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值及所述血压舒张压预测值输入所述血糖预测模型,预测所述待测对象的血糖值。
根据本申请的再一方面,还提供了一种用于血糖预测的机器学习模型样本获取的装置,其中,所述装置包括:
第一模块,用于获取采集对象的状态信息、心电信号及其对应的血糖值,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
第二模块,用于将所述心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
第三模块,用于将所述血糖值标识为所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值和所述血压舒张压预测值的标签,并将所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值、所述血压舒张压预测值及所述标签确定为所述机器学习模型样本。
可选地,所述一种用于血糖预测的机器学习模型样本获取的装置还包括:
第四模块,用于基于获取的若干采集对象的状态信息、心电信号及其对应的血糖值,确定若干所述机器学习模型样本,以构建所述机器学习模型的样本集。
可选地,所述装置还包括:
第五模块,用于基于所述样本集,训练所述机器学习模型,当所述机器学习模型的评估指标符合预设阈值,将训练后的机器学习模型确定为血糖预测模型。
可选地,所述装置还包括:
第六模块,用于获取待测对象的状态信息及心电信号,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
第七模块,用于将所述心电信号的第一特征值、第二特征值及第三特征值分别输入所述血压收缩压预测模型和所述血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
第八模块,用于将所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值及所述血压舒张压预测值输入所述血糖预测模型,预测所述待测对象的血糖值。
与现有技术相比,本申请通过一种用于血糖预测的机器学习模型样本获取的技术方案,获取心电信号及其对应的血糖值;确定所述心电信号的第一特征值、第二特征值及第三特征值;将所述心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;将所述血糖值确定为所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值和所述血压舒张压预测值的标签,并将所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值、所述血压舒张压预测值及所述标签确定为所述机器学习模型样本。进一步地,基于获取的若干心电信号及其对应的血糖值,确定若干所述机器学习模型样本,以构建所述机器学习模型的样本集,其中,同一样本集中各样本对应的采集对象的状态相同。并基于所述样本集训练所述机器学习模型,以确定基于心电信号及其预测血压来预测待测对象相应状态血糖值的血糖预测模型。通过该技术方案,基于心电信号的多维特征值来预测血压值,并基于心电信号的多位特征值及预测血压,进一步预测待测对象相应状态下的血糖值,无需再单独测量,提升了用户体验,而且基于多样性的样本集训练确定的血糖预测模型可用于对不同人的心电信号进行血糖预测,具有普遍适用性。
进一步地,还可以将采集对象的状态信息作为样本中的参数之一,基于这样的样本训练后获得的血糖预测模型,可适用待测对象任一状态的血糖值预测,无需针对待测对象的状态采用相应的血糖预测模型来预测,提升了预测效率及用户体验。
附图说明
通过阅读参照以下附图所作的对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1示出典型心电信号波形示意图;
图2示出根据本申请一个方面的一种用于血糖预测的机器学习模型样本获取的方法流程图;
图3示出根据本申请又一个方面的一种用于血糖预测的机器学习模型样本获取的装置示意图;
图4示出根据本申请另一个方面的一种用于血糖预测的机器学习模型样本获取的方法流程图;
图5示出根据本申请再一个方面的一种用于血糖预测的机器学习模型样本获取的装置示意图;
附图中相同或相似的附图标签代表相同或相似的部件。
具体实施方式
下面结合附图对本发明作进一步详细描述。
在本申请一个典型的配置中,设备、装置各模块和/或可信方均包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括非暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
为更进一步阐述本申请所采取的技术手段及取得的效果,下面结合附图及可选实施例,对本申请的技术方案,进行清楚和完整的描述。
图2示出本申请一个方面的一种用于血糖预测的机器学习模型样本获取的方法流程图,其中,一个实施例的方法包括:
S21获取心电信号及其对应的血糖值,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
S22将所述心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
S23将所述血糖值标识为所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值和所述血压舒张压预测值的标签,并将所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值、所述血压舒张压预测值及所述标签确定为所述机器学习模型样本。
在本申请中,所述方法通过设备100执行,其中,设备100为安装有相关计算机程序的计算机设备和/或云,所述计算机设备包括但不限于个人计算机、笔记本电脑、工业计算机、网络主机、单个网络服务器、多个网络服务器集;所述云由基于云计算(CloudComputing)的大量计算机或网络服务器构成,其中,云计算是分布式计算的一种,由一群松散耦合的计算机集组成的一个虚拟超级计算机。
在此,所述计算机设备和/或云仅为举例,其他现有的或者今后可能出现的设备和/或资源共享平台如适用于本申请也应包含在本申请的保护范围内,在此,以引用的方式包含于此。
在该实施例中,在所述步骤S21中,设备100获取采集自采集对象的心电信号及其对应的血糖值,其中,获取的心电信号可以是通过专用心电信号采集设备采集后以数字格式保存的心电信号波形(电压-时间波形,即心电图),或者是通过专用心电信号采集设备采集后传输给设备100,设备100获取后转换成数字格式的心电信号波形。其中,获取的心电信号可以是单导联心电信号,也可以是多导联心电信号,其中,为了确保可获取到足够特征信息及其对应的特征值,每个导联心电信号的采样频率不少于100Hz,心电信号至少持续25秒。其中,获取的血糖值是与该心电信号同步采集确定的。
在步骤S21中,还对获取的心电信号波形进行处理,标注整段心电信号波形中的特征信息,确定该心电信号的第一特征值、第二特征值及第三特征值。其中,在特征信息标注完成后,需要对整段心电信号波形以及特征信息进行检验,排除掉无效的心电信号波形、心电信号波形中无效的心跳波形以及标注错误的特征信息,以确保后续确定的各类特征值的正确性。优选地,排除掉无效的心电信号波形、心电信号波形中无效的心跳波形以及标注错误的特征信息后,保留下来的正常心电信号波形中至少包括5组两两相邻的心跳波形,以保证可确定足够数量的各特征值。保留下来的正常心电信号波形中包括的两两相邻的心跳波形越多,相应地获得的各特征值的数量越多。
可选地,其中,确定所述心电信号的第一特征值包括:
确定所述心电信号中每个心跳波形;
基于每个心跳波形,获取每个心跳波形的特征值,
遍历所述心电信号的全部心跳波形,获取全部心跳波形的特征值的统计值,并将所述统计值确定为所述心电信号的第一特征值。
其中,设备100获取的心电信号波形通常由若干个心跳波形组成,设备100对心电信号波形进行处理,标注整段心电信号波形中的特征信息,包括确定心电信号波形中的每个心跳波形,标注出每个心跳波形的特征信息。每个心跳波形的特征信息包括:P波、PR间期、QRS、ST段、T波、U波、TP段及QT间期等。获取每个心跳波形的各类特征信息相应的特征值,在此基础上,设备100遍历整段心电信号波形的全部心跳波形,获取全部心跳波形的各类特征信息对应的特征值各自的统计值,并将各统计值确定为该段心电信号的第一特征值。
可选地,其中,所述每个心跳波形的特征值包括以下至少一项:
每个心跳波形的P波时长;
每个心跳波形的PR间期;
每个心跳波形的QRS时间;
每个心跳波形的ST段时长;
每个心跳波形的T波时长;
每个心跳波形的TP段时长;
每个心跳波形的QT间期;
每个心跳波形的U波时长。
其中,每个心跳波形的特征信息包括:P波、PR间期、QRS、ST段、T波、U波、TP段及QT间期等,每个心跳波形各特征信息对应的特征值可包括:P波对应的P波时长、PR间期对应的PR间期、QRS对应的QRS时间、ST段对应的ST段时长、T波对应的T波时长、TP段对应的TP段时长、QT间期对应的QT间期及U波对应的U波时长。整段心电信号的所有心跳波形的各特征值对应的统计值可以包括:整段心电信号的所有心跳波形的P波对应的特征值P波时长的统计值、整段心电信号的所有心跳波形的PR间期对应的特征值PR间期的统计值、整段心电信号的所有心跳波形的QRS对应的特征值QRS时间的统计值、整段心电信号的所有心跳波形的ST段对应的特征值ST段时长的统计值、整段心电信号的所有心跳波形的T波对应的特征值T波时长的统计值、整段心电信号的所有心跳波形的TP段对应的特征值TP段时长的统计值、整段心电信号的所有心跳波形的QT间期对应的特征值QT间期的统计值以及整段心电信号的所有心跳波形的U波对应的特征值U波时长的统计值。
可选地,其中,所述全部心跳波形的特征值的统计值包括以下至少一项:
均方根;
中值绝对偏差;
中值变异系数。
其中,所述统计值可以是整段心电信号的所有心跳波形的各类特征信息对应的特征值各自的均方根,和/或其中值绝对偏差,和/或中值变异系数,所述统计值也可以是整段心电信号的所有心跳波形的各类型特征信息对应的特征值各自的其它统计值,在此,不做具体限定,任何适用于本申请中的整段心电信号的所有心跳波形的各类特征信息对应的特征值的统计值形式都应处于本申请的保护范围内。
可选地,其中,确定所述心电信号的第二特征值包括:
确定所述心电信号中每个心跳波形;
基于两两相邻的两个心跳波形,获取两两相邻的两个心跳波形的特征值,
遍历所述心电信号的全部心跳波形,获取全部两两相邻的两个心跳波形的特征值的统计值,并将所述统计值确定为所述心电信号的第二特征值。
其中,设备100在对心电信号波形进行处理,标注出整段心电信号的每个心跳波形的特征信息如P波、PR间期、QRS、ST段、T波、U波、TP段及QT间期的基础上,根据两两相邻的两个心跳波形(连续的两两相邻的两个心跳波形,不允许拼接得到两两相邻的两个心跳波形),获取两两相邻的两个心跳波形同类特征信息对应的特征值,在此基础上,设备100遍历整段心电信号波形的全部心跳波形,获取全部两两相邻的两个心跳波形的各个同类特征信息对应的特征值各自的统计值,并将各统计值确定为该段心电信号的第二特征值。
可选地,其中,所述两两相邻的两个心跳波形的特征值包括以下至少一项:
两两相邻的心跳波形的PP间期;
两两相邻的心跳波形的QQ间期;
两两相邻的心跳波形的RR间期;
两两相邻的心跳波形的SS间期;
两两相邻的心跳波形的TT间期;
两两相邻的心跳波形的UU间期。
其中,心电信号的每个心跳波形中的特征信息可以包括:P波、PR间期、QRS、ST段、T波、QT间期、TP段及U波,两两相邻的两个心跳波形各特征信息对应的特征值分别为PP间期、QQ间期、RR间期、SS间期、TT间期及UU间期,整段心电信号的所有两两相邻的两个心跳波形各特征信息对应的各特征值的统计值可以包括:所有两两相邻的两个心跳波形的P波对应的PP间期的统计值、所有两两相邻的两个心跳波形的Q点对应的QQ间期的统计值、所有两两相邻的两个心跳波形的R点对应的RR间期的统计值、所有两两相邻的两个心跳波形的S点对应的SS间期的统计值、所有两两相邻的两个心跳波形的T波对应的TT间期的统计值及所有两两相邻的两个心跳波形的U波对应的UU间期的统计值。
其中,若整段心电信号对应的心电信号波形标注的特征信息中存在异常的特征信息,需删除异常的特征信息及其所在的整个心跳波形,此种情况下,整段心电信号的心电信号波形将被分割成若干段,则分别确定这若干段心电信号波形的两两相邻的两个心跳波形的各类特征信息对应的特征值,然后汇总,确定统计值。
可选地,其中,所述全部两两相邻的两个心跳波形的特征值的统计值包括以下至少一项:
均方根;
中值绝对偏差;
中值变异系数;
大于50毫秒的连续间隔差异数。
其中,所述统计值可以是整段心电信号的所有两两相邻的两个心跳波形的各类特征信息对应的特征值的均方根,和/或其中值绝对偏差,和/或中值变异系数,和/或其中大于50毫秒的连续间隔差异数,所述统计值也可以是整段心电信号的所有两两相邻的两个心跳波形的各类特征信息对应的特征值的其它统计值,在此,不做具体限定,任何适用于本申请中的整段心电信号的所有两两相邻的两个心跳波形的各类特征信息对应的特征值的统计值形式都应处于本申请的保护范围内。
可选地,其中,确定所述心电信号的第三特征值包括:
获取所述第二特征值对应的时域波形的功率频谱;
基于所述功率频谱,获取若干预设频域的功率值;
将所述若干预设频域的功率值之间的比值确定为所述心电信号的第三特征值。
在心电信号研究领域,两两相邻的心跳波形的PP间期、QQ间期、RR间期、SS间期、TT间期和/或UU间期等本申请中所述的第二特征值对应的时域波形可以反映采集对象的人体交感与副交感神经的活性,因此,第二特征值对应的时域波形,与心电信号的时域波形相比较,更能体现与血压(舒张压及收缩压)的关联性。第二特征值对应的时域波形的功率频谱反映的是第二特征值对应的时域波形速度的频率分布,可以从其功率频谱中获取更多的特征用于本申请中所述的机器学习模型。
其中,在获取所述第二特征值对应的时域波形的功率频谱之前,若整段心电信号对应的心电信号波形标注的特征信息中存在异常的特征信息,需删除异常的特征信息及其所在的整个心跳波形,此种情况下,整段心电信号的心电信号波形将被分割成若干段,分别获取这若干段心电信号波形中的第二特征值对应的时域波形的功率频谱,然后将获得的各功率频谱进行叠加。
结合第二特征值对应的时域波形的特点,将第二特征值对应的时域波形转换成功率频谱后,功率频谱的大部分能量(功率)集中在0.04~0.15赫兹(Hz)的频域范围内,为了便于后续机器学习模型的学习,可将此频域范围划分成若干预设频域,获取每个预设频域内的功率值,其中,各预设频域内的功率值之和归一化后为1。将各预设频域内的功率值的比值确定为所述心电信号的第三特征值。比如,可按如下划分预设频域:VLF(Very LowFrequency,甚低频),低于0.04Hz;LF(Low Frequency,低频),0.04~0.15Hz(含两端点);HF(High Frequency,高频),大于0.15Hz,其中,设这三个预设频域内的功率值归一化后分别为P1、P2及P3,P1+P2+P3=1,P1/P2=a,P2/P3=b,P1/P3=c,可将a、b及c确定为所述心电信号的第三特征值。
为了获得预测准确率更高的血糖预测模型,可选地,其中,确定所述心电信号的第三特征值还包括:
基于所述功率频谱,计算所述功率频谱的熵值;
将所述功率频谱的熵值确定为所述心电信号的第三特征值。
其中,在获取整段心电信号波形中所有正常心跳信号波形的功率频谱后,还计算该功率平铺的熵值,并将该功率频谱的熵值也确定为所述心电信号的第三特征值,即,第三特征值可以包括:所述心电信号的功率频谱在若干预设频域内的功率值的比值,以及该功率频谱的熵值。
继续在该实施例中,在所述步骤S22中,设备100将所述心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值。
其中,所述预设的血压收缩压预测模型和血压舒张压预测模型是事先已经完成训练的机器学习模型(比如,多元线性回归模型、随机森林回归模型等),无需单独测量采集对象的血压值(包括收缩压和舒张压),可根据获取到的采集对象的心电信号来确定该心电信号对应的上述第一特征值、第二特征值及第三特征值,将该心电信号的第一特征值、第二特征值及第三特征值输入预设的血压收缩压预测模型可以获取该心电信号对应的血压收缩压预测值;将该心电信号的第一特征值、第二特征值及第三特征值输入预设的血压舒张压预测模型可以获取该心电信号对应的血压舒张压预测值。
继续在该实施例中,在所述步骤S23中,设备100将在步骤S21中与心电信号同期获取的血糖值标识为在步骤S21中确定的该心电信号的第一特征值、第二特征值及第三特征值及在步骤S22中获得的该心电信号对应的血压收缩压预测值、述血压舒张压预测值的标签,并将该心电信号的第一特征值、第二特征值、第三特征值、该心电信号对应的血压收缩压预测值、述血压舒张压预测值及所述标签确定为所述机器学习模型样本。
可选地,所述实施例的方法还包括:
S24(未示出)获取若干心电信号及其对应的血糖值,采用上述实施例的方法来确定若干所述机器学习模型样本,以构建所述机器学习模型的样本集,用于所述机器学习模型的训练。
其中,所选的样本集中各样本的心电信号及其对应的血糖值采集自处于同样状态的采集对象(比如采集对象都是处于空腹状态、餐后1小时状态或者餐后2小时状态)。为了使得机器学习模型的样本具有更广泛的多样性,以使得训练后的机器学习模型适用不同人群,还可以基于不同性别、年龄段来选择采集对象,比如,可以分别按性别(男性和女性)、年龄段(比如:18~30岁、31~60岁及61岁以上)采集不同采集对象在空腹状态下的心电信号及其同期的血糖值;分别按性别(男性和女性)、年龄段(比如:18~30岁、31~60岁及61岁以上)采集不同采集对象在餐后1小时状态下的心电信号及其同期的血糖值;分别按性别(男性和女性)、年龄段(比如:18~30岁、31~60岁及61岁以上)采集不同采集对象在餐后2小时状态下的心电信号及其同期的血糖值。不同性别、各年龄段的样本在样本集中的数量也应该基本相同。
其中,若采集的原始心电信号和/或其对应的血糖值异常,比如,获取的心电信号持续时间少于25秒、无法根据获取的心电信号确定足够的特征点来获得其第一特征值、第二特征值,则放弃该原始心电信号及其对应的血糖值。
其中,所选的样本集中的样本的第一特征值、第二特征值中包括的特征值越多,训练机器学习模型后得到的血糖预测模型的预测准确度越高。
可选地,所述实施例的方法还包括:
S25(未示出)基于所述样本集,训练所述机器学习模型,当所述机器学习模型的评估指标符合预设阈值,将训练后的机器学习模型确定为血糖预测模型。
其中,一个可选实施例中,所述机器学习模型可以是多元线性回归模型。多元线性回归模型可用如下公式表达:
y=b0+b1x1+b2x2+…+bkxk+e
其中,b0为常数项,b1、b2…bk为偏回归系数(b1为x1、x2…xk固定时,自变量x1每增加一个单位对因变量y的效应;同理b2为x1、x2…xk固定时,自变量x2每增加一个单位对因变量y的效应,同理可知其它偏回归系数的含义),e为误差项。
在本可选实施例中,该多元线性回归模型可用如下公式表达:
s=a0+a1x1+a2x2+…+aixi+b1y1+b2y2+…+bjyj+c1z1+c2z2+…+ckzk+du+fv+e
其中,a0为常数项,a1、a2…ai为心电信号的多个第一特征值x1、x2…xi的偏回归系数;b1、b2…bj为心电信号的多个第二特征值y1、y2…yj的偏回归系数;c1、c2…ck为心电信号的多个第三特征值z1、z2…zk的偏回归系数;d为血压收缩压预测值u的偏回归系数;f为血压舒张压预测值v的偏回归系数;s为回归值,即基于该多元线性回归模型回归得到的心电信号对应的血糖值。
利用样本集中的样本,可以采用SVM(Support Vector Machines,支持向量机)、SGD(Stochastic Gradient Descent,随机梯度下降)等方法来训练上述多元线性回归模型,基于预设的多元线性回归模型的评估指标阈值,来确定该多元线性回归模型的常数项及各偏回归系数,并将符合评估指标阈值的训练后的多元线性回归模型确定为血糖预测模型。其中,基于上述样本集中的样本,采用样本中的心电信号的第一特征值、第二特征值、第三特征值、血压收缩压预测值、血压舒张压预测值及血糖值标签,对该多元线性回归模型进行训练,当训练后的多元线性回归模型的评估指标符合预设阈值,将训练后的多元线性回归模型确定为用于预测待测对象相应状态下血糖值的血糖预测模型。
其中,所述多元线性回归模型的评估指标是用来评估该多元线性回归模型拟合度的指标,该评估指标可以是回归值的MSE(Mean Squared Error,均方误差),也可以是回归值的RMSE(Root Mean Squared Error,均方根误差),当评估指标符合预设阈值,表明基于样本集训练后的多元线性回归模型的拟合度满足要求,可将训练后的多元线性回归模型确定为血糖预测模型。如果基于样本集多次迭代训练后的多元线性回归模型的评估指标无法满足预设阈值,可以将最优评估指标对应的训练后的多元线性回归模型确定为血糖预测模型,比如,如采用MSE作为评估指标,则将MSE最小的训练后的多元线性回归模型确定为血糖预测模型。
其中,另一个可选实施例中,所述机器学习模型还可以是随机森林回归模型,基于样本集做多次放回样本的抽样,获得多个子样本集,然后基于多个子样本集中各样本的第一特征值、第二特征值、第三特征值、血压收缩压预测值及血压舒张压预测值,分别对随机森林回归模型进行连续的决策与机会节点训练,生成各自对应的决策树,其中,每个决策树输出一个结果,再基于多个子样本集中的标记,对每个决策树进行投票,根据投票结果,将评估指标满足预设阈值的决策树对应的随机森林回归模型确定为血糖预测模型。其中,基于多个子样本集中血糖值标签,对每个决策树进行投票,根据投票结果,将评估指标满足预设阈值的决策树对应的随机森林回归模型确定为用于预测血糖值的血糖预测模型。
可选地,上述实施例的方法还包括:
S26(未示出)获取待测对象的心电信号,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
S27(未示出)将所述心电信号的第一特征值、第二特征值及第三特征值分别输入所述血压收缩压预测模型和所述血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
S28(未示出)将所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值及所述血压舒张压预测值输入所述血糖预测模型,预测所述心电信号对应的血糖值。
其中,设备100获取采集对象的一段持续时间不少于25秒、包含若干完整心跳波形的心电信号,参照前述实施例的方法确定其特征点,获得其第一特征值、第二特征值及第三特征值,然后将获得的该心电信号的第一特征值、第二特征值及第三特征值分别输入所述血压收缩压预测模型和所述血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值,再将该心电信号的第一特征值、第二特征值及第三特征值、获得的血压收缩压预测值及血压舒张压预测值输入上述血糖预测模型,以预测该心电信号对应的血糖值。
需要说明的是,训练后获得的血糖预测模型只适用于与其样本集中样本的采集对象处于相同状态的待测对象。比如,若样本集中的样本来自于处于空腹状态的采集对象,则基于这个样本集训练后获得的血糖预测模型只适用于预测待测对象在空腹状态下的血糖值,步骤S26应是获取处于空腹状态的待测对象的心电信号,确定其多维特征值,经过步骤S27、S28后获得该待测对象空腹状态下的血糖预测值。同样的,若样本集中的样本来自于处于餐后1小时状态的采集对象,则基于这个样本集训练后获得的血糖预测模型只适用于预测待测对象在餐后1小时状态下的血糖值,步骤S26应是获取处于餐后1小时状态的待测对象的心电信号,确定其多维特征值,经过步骤S27、S28后获得该待测对象餐后1小时状态下的血糖预测值。若样本集中的样本来自于处于餐后2小时状态的采集对象,则基于这个样本集训练后获得的血糖预测模型只适用于预测待测对象在餐后2小时状态下的血糖值,步骤S26应是获取处于餐后2小时状态的待测对象的心电信号,确定其多维特征值,经过步骤S27、S28后获得该待测对象餐后2小时状态下的血糖预测值。
图3示出根据本申请又一个方面的一种用于血糖预测的机器学习模型样本获取的装置示意图,其中,一个实施例的所述装置包括:
第一模块31,用于获取心电信号及其对应的血糖值,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
第二模块32,用于将所述心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
第三模块33,用于将所述血糖值标识为所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值和所述血压舒张压预测值的标签,并将所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值、所述血压舒张压预测值及所述标签确定为所述机器学习模型样本。
在该实施例中,所述装置用于实施前述实施例的方法,所述装置集成在与前述设备100软硬件配置环境相同的设备中。
其中,该装置的第一模块31获取采集自采集对象的心电信号及其对应的血糖值,其中,获取的心电信号可以是通过专用心电信号采集设备采集后以数字格式保存的心电信号波形(电压-时间波形,即心电图),或者是通过专用心电信号采集设备采集后传输给该装置,该装置获取后转换成数字格式的心电信号波形。其中,获取的心电信号可以是单导联心电信号,也可以是多导联心电信号,其中,为了确保可获取到足够特征信息及其对应的特征值,每个导联心电信号的采样频率不少于100Hz,至少持续25秒的心电信号。其中,获取的血糖值是与该心电信号同步采集确定的。
第一模块31还对获取的心电信号波形进行处理,标注整段心电信号波形中的特征信息,确定该心电信号的第一特征值、第二特征值及第三特征值。其中,在特征信息标注完成后,需要对整段心电信号波形以及特征信息进行检验,排除掉无效的心电信号波形、心电信号波形中无效的心跳波形以及标注错误的特征信息,以确保后续确定的各类特征值的正确性。
该装置的第二模块32将第一模块31确定的该心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与该心电信号对应的血压收缩压预测值和血压舒张压预测值。
该装置的第三模块33将在第一模块31中与心电信号同期获取的血糖值标识为在第一模块31中确定的该心电信号的第一特征值、第二特征值、第三特征值、在第二模块32中获得的血压收缩压预测值及血压舒张压预测值的标签,并将该心电信号的第一特征值、第二特征值、第三特征值、血压收缩压预测值、血压舒张压预测值及该标签确定为所述机器学习模型样本。
可选地,所述装置还包括:第四模块34(未示出),用于获取若干心电信号及其对应的血糖值,采用上述方法实施例的方法来确定若干所述机器学习模型样本,以构建所述机器学习模型的样本集,用于所述机器学习模型的学习。
可选地,所述装置还包括:第五模块35(未示出),用于基于第四模块34中构建的所述样本集,训练所述机器学习模型,当所述机器学习模型的评估指标符合预设阈值,将训练后的机器学习模型确定为血糖预测模型。其中,基于上述样本集中的样本,采用样本中的心电信号的第一特征值、第二特征值、第三特征值、血压收缩压预测值、血压舒张压预测值及血糖值标签,对该机器学习模型进行训练,当训练后的机器学习模型的评估指标符合预设阈值,将训练后的机器学习模型确定为用于预测待测对象相应状态下血糖值的血糖预测模型。
可选地,所述装置还包括:第六模块36(未示出),用于获取待测对象的心电信号,并确定该心电信号的第一特征值、第二特征值及第三特征值;第七模块37(未示出),用于将该心电信号的第一特征值、第二特征值及第三特征值分别输入所述血压收缩压预测模型和所述血压舒张压预测模型,获得与该心电信号对应的血压收缩压预测值和血压舒张压预测值;第八模块38(未示出),用于将该心电信号的第一特征值、第二特征值、第三特征值、第其模块获得的血压收缩压预测值及血压舒张压预测值输入上述血糖预测模型,预测所述心电信号对应的血糖值。
上述方法实施例需要根据采集对象的状态来将若干样本划分成针对不同状态的多个样本集,以训练机器学习模型来获得多个适用于预测待测对象不同状态下的血糖预测模型。
为了进一步提升预测效率及用户体验,图4示出根据本申请另一个方面的一种用于血糖预测的机器学习模型样本获取的方法流程图,其中,一个实施例的所述方法包括:
S41获取采集对象的状态信息、心电信号及其对应的血糖值,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
S42将所述心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
S43将所述血糖值标识为所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值和所述血压舒张压预测值的标签,并将所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值、所述血压舒张压预测值及所述标签确定为所述机器学习模型样本。
在该实施例中,在所述步骤S41中,设备100获取采集自采集对象的心电信号及其对应的血糖值,以及获取采集时采集对象所处的状态信息,其中,获取的心电信号及其对应的血糖值与前述方法实施例中相同,在此不再赘述。
在步骤S41中,还对获取的心电信号波形进行处理,标注整段心电信号波形中的特征信息,确定该心电信号的第一特征值、第二特征值及第三特征值。其中,确定该心电信号的第一特征值、第二特征值及第三特征值与前述方法实施例中相同,在此不再赘述。
继续在该实施例中,所述步骤S42与前述方法实施例中的步骤S22相同,在此不再赘述。
继续在该实施例中,在所述步骤S43中,设备100将在步骤S21中与心电信号同期获取的血糖值标识为在步骤S21中获取的采集对象的状态信息、确定的该心电信号的第一特征值、第二特征值及第三特征值及在步骤S42中获得的该心电信号对应的血压收缩压预测值、述血压舒张压预测值的标签,并将该状态信息、该心电信号的第一特征值、第二特征值、第三特征值、该心电信号对应的血压收缩压预测值、述血压舒张压预测值及所述标签确定为所述机器学习模型样本。
在该实施例中,将采集对象的状态信息作为输入参数之一,无需再将包括状态信息的样本基于采集对象的状态划分针对不同状态的样本集,采用由包括状态信息的样本的样本集训练机器学习模型,获得的血糖预测模型可适用采集对象的任一状态,可让用户的使用体验得到提升。其中,采集对象的状态可以是空腹、餐后1小时和餐后2小时,对应的状态信息可以分别用数值来标识,比如分别用0、1、2标识。
可选地,该实施例的方法还包括:
S44(未示出)获取若干采集对象的状态信息、心电信号及其对应的血糖值,采用上述步骤来确定若干所述机器学习模型样本,以构建所述机器学习模型的样本集,用于所述机器学习模型的训练。
可选地,该实施例的方法还包括:
S45(未示出),基于所述样本集,训练所述机器学习模型,当所述机器学习模型的评估指标符合预设阈值,将训练后的机器学习模型确定为血糖预测模型。其中,所述机器学习模型及其训练、评估等参考前述方法实施例的步骤S25,在此不再赘述。
可选地,该实施例的方法还包括:
S46(未示出)获取待测对象的状态信息及心电信号,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
S47(未示出)将所述心电信号的第一特征值、第二特征值及第三特征值分别输入所述血压收缩压预测模型和所述血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
S48(未示出)将所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值及所述血压舒张压预测值输入所述血糖预测模型,预测所述待测对象的血糖值。
图5示出根据本申请再一个方面的一种用于血糖预测的机器学习模型样本获取的装置,其中,一个实施例的所述装置包括:
第一模块51,用于获取采集对象的状态信息、心电信号及其对应的血糖值,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
第二模块52,用于将所述心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
第三模块53,用于将所述血糖值标识为所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值和所述血压舒张压预测值的标签,并将所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值、所述血压舒张压预测值及所述标签确定为所述机器学习模型样本。
可选地,所述装置还包括:
第四模块54(未示出),用于获取若干采集对象的状态信息、心电信号及其对应的血糖值,采用上述方法实施例的方法来确定若干所述机器学习模型样本,以构建所述机器学习模型的样本集,用于所述机器学习模型的学习。
可选地,所述装置还包括:
第五模块55(未示出),用于基于第四模块34中构建的所述样本集,训练所述机器学习模型,当所述机器学习模型的评估指标符合预设阈值,将训练后的机器学习模型确定为血糖预测模型。
其中,基于上述样本集中的样本,采用样本中的状态信息、心电信号的第一特征值、第二特征值、第三特征值、血压收缩压预测值、血压舒张压预测值及血糖值标签,对该机器学习模型进行训练,当训练后的机器学习模型的评估指标符合预设阈值,将训练后的机器学习模型确定为用于预测待测对象任一状态下血糖值的血糖预测模型。
可选地,所述装置还包括:
第六模块56(未示出),用于获取待测对象的状态信息及心电信号,并确定该心电信号的第一特征值、第二特征值及第三特征值;
第七模块57(未示出),用于将该心电信号的第一特征值、第二特征值及第三特征值分别输入所述血压收缩压预测模型和所述血压舒张压预测模型,获得与该心电信号对应的血压收缩压预测值和血压舒张压预测值;
第八模块58(未示出),用于将该状态信息、该心电信号的第一特征值、第二特征值、第三特征值、第其模块获得的血压收缩压预测值及血压舒张压预测值输入上述血糖预测模型,预测所述待测对象的血糖值。
根据本申请的又一方面,还提供了一种计算机可读介质,所述计算机可读介质存储有计算机可读指令,所述计算机可读指令可被处理器执行以实现前述方法。
需要注意的是,本申请可在软件和/或软件与硬件的组合体中被实施,例如,可采用专用集成电路(ASIC)、通用目的计算机或任何其他类似硬件设备来实现。在一个实施例中,本申请的软件程序可以通过处理器执行以实现上文所述步骤或功能。同样地,本申请的软件程序(包括相关的数据结构)可以被存储到计算机可读记录介质中,例如,RAM存储器,磁或光驱动器或软磁盘及类似设备。另外,本申请的一些步骤或功能可采用硬件来实现,例如,作为与处理器配合从而执行各个步骤或功能的电路。
另外,本申请的一部分可被应用为计算机程序产品,例如计算机程序指令,当其被计算机执行时,通过该计算机的操作,可以调用或提供根据本申请的方法和/或技术方案。而调用本申请的方法的程序指令,可能被存储在固定的或可移动的记录介质中,和/或通过广播或其他信号承载媒体中的数据流而被传输,和/或被存储在根据所述程序指令运行的计算机设备的工作存储器中。在此,根据本申请的一个实施例包括一个设备,该设备包括用于存储计算机程序指令的存储器和用于执行程序指令的处理器,其中,当该计算机程序指令被该处理器执行时,触发该设备运行基于前述根据本申请的多个实施例的方法和/或技术方案。
根据本申请的又一方面,还提供了一种用于用于血糖预测的机器学习模型样本获取的设备,其中,该设备包括:
一个或多个处理器;以及
存储有计算机可读指令的存储器,所述计算机可读指令在被执行时使所述处理器执行如前述方法的操作。
例如,计算机可读指令在被执行时使所述一个或多个处理器:获取心电信号及其对应的血糖值,并确定该心电信号的第一特征值、第二特征值及第三特征值,其中,所述心电信号的第一特征值包括所述心电信号的全部心跳波形的特征值的统计值,所述心电信号的第二特征值包括所述心电信号各段连续正常心电信号波形中,全部两两相邻的两个心跳波形的特征值的统计值,所述心电信号的第三特征值包括所述心电信号对应的功率频谱中若干预设频段的功率之间的比值,还包括所述功率频谱的熵值;将该心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与该心电信号对应的血压收缩压预测值和血压舒张压预测值;将所述血糖值标识为该心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值和所述血压舒张压预测值的标签,并将该心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值、所述血压舒张压预测值及所述标签确定为所述机器学习模型样本。进一步地,基于获取的若干心电信号及其对应的血糖值,采用上述方法确定若干所述机器学习模型样本,以构建所述机器学习模型的样本集。进一步地,基于所述样本集,训练所述机器学习模型,当所述机器学习模型的评估指标符合预设阈值,将训练后的机器学习模型确定为血糖预测模型。进一步地,确定获取到的待测对象的心电信号的第一特征值、第二特征值及第三特征值;将该述心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与该心电信号对应的血压收缩压预测值和血压舒张压预测值;将该心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值、所述血压舒张压预测值输入所述血糖预测模型,预测所述心电信号对应的血糖值。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化涵括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。此外,显然“包括”一词不排除其他单元或步骤,单数不排除复数。装置权利要求中陈述的多个单元或装置也可以由一个单元或装置通过软件和/或者硬件来实现。第一,第二等词语用来表示名称,而并不表示任何特定的顺序。

Claims (18)

1.一种用于血糖预测的机器学习模型样本获取的方法,其特征在于,所述方法包括:
获取心电信号及其对应的血糖值,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
将所述心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
将所述血糖值标识为所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值和所述血压舒张压预测值的标签,并将所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值、所述血压舒张压预测值及所述标签确定为所述机器学习模型样本。
2.根据权利要求1所述的方法,其特征在于,确定所述心电信号的第一特征值包括:
确定所述心电信号中每个心跳波形;
基于每个心跳波形,获取每个心跳波形的特征值,
遍历所述心电信号的全部心跳波形,获取全部心跳波形的特征值的统计值,并将所述统计值确定为所述心电信号的第一特征值。
3.根据权利要求2所述的方法,其特征在于,所述每个心跳波形的特征值包括以下至少一项:
每个心跳波形的P波时长;
每个心跳波形的PR间期;
每个心跳波形的QRS时间;
每个心跳波形的ST段时长;
每个心跳波形的T波时长;
每个心跳波形的TP段时长;
每个心跳波形的QT间期;
每个心跳波形的U波时长。
4.根据权利要求2所述的方法,其特征在于,所述统计值包括以下至少一项:
均方根;
中值绝对偏差;
中值变异系数。
5.根据权利要求1所述的方法,其特征在于,确定所述心电信号的第二特征值包括:
确定所述心电信号中每个心跳波形;
基于两两相邻的两个心跳波形,获取两两相邻的两个心跳波形的特征值,
遍历所述心电信号的全部心跳波形,获取全部两两相邻的两个心跳波形的特征值的统计值,并将所述统计值确定为所述心电信号的第二特征值。
6.根据权利要求5所述的方法,其特征在于,所述两两相邻的两个心跳波形的特征值包括以下至少一项:
两两相邻的心跳波形的PP间期;
两两相邻的心跳波形的QQ间期;
两两相邻的心跳波形的RR间期;
两两相邻的心跳波形的SS间期;
两两相邻的心跳波形的TT间期;
两两相邻的心跳波形的UU间期。
7.根据权利要求5所述的方法,其特征在于,所述统计值包括以下至少一项:
均方根;
中值绝对偏差;
中值变异系数;
大于50毫秒的连续间隔差异数。
8.根据权利要求5所述的方法,其特征在于,确定所述心电信号的第三特征值包括:
获取所述第二特征值对应的时域波形的功率频谱;
基于所述功率频谱,获取若干预设频域的功率值;
将所述若干预设频域的功率值之间的比值确定为所述心电信号的第三特征值。
9.根据权利要求8所述的方法,其特征在于,确定所述心电信号的第三特征值还包括:
基于所述功率频谱,计算所述功率频谱的熵值;
将所述功率频谱的熵值确定为所述心电信号的第三特征值。
10.根据权利要求1至9中任一项所述的方法,其特征在于,所述方法还包括:
基于获取的若干心电信号及其对应的血糖值,确定若干所述机器学习模型样本,以构建所述机器学习模型的样本集。
11.根据权利要求10所述的方法,其特征在于,所述方法还包括:
基于所述样本集,训练所述机器学习模型,当所述机器学习模型的评估指标符合预设阈值,将训练后的机器学习模型确定为血糖预测模型。
12.根据权利要求11所述的方法,其特征在于,所述方法还包括:
获取待测对象的心电信号,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
将所述心电信号的第一特征值、第二特征值及第三特征值分别输入所述血压收缩压预测模型和所述血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
将所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值及所述血压舒张压预测值输入所述血糖预测模型,预测所述待测对象的血糖值。
13.一种用于血糖预测的机器学习模型样本获取的方法,其特征在于,所述方法包括:
获取采集对象的状态信息、心电信号及其对应的血糖值,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
将所述心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
将所述血糖值标识为所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值和所述血压舒张压预测值的标签,并将所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值、所述血压舒张压预测值及所述标签确定为所述机器学习模型样本。
14.根据权利要求13所述的方法,其特征在于,所述方法还包括:
基于获取的若干采集对象的状态信息、心电信号及其对应的血糖值,确定若干所述机器学习模型样本,以构建所述机器学习模型的样本集。
15.根据权利要求14所述的方法,其特征在于,所述方法还包括:
基于所述样本集,训练所述机器学习模型,当所述机器学习模型的评估指标符合预设阈值,将训练后的机器学习模型确定为血糖预测模型。
16.根据权利要求15所述的方法,其特征在于,所述方法还包括:
获取待测对象的状态信息及心电信号,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
将所述心电信号的第一特征值、第二特征值及第三特征值分别输入所述血压收缩压预测模型和所述血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
将所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值及所述血压舒张压预测值输入所述血糖预测模型,预测所述待测对象的血糖值。
17.一种用于血糖预测的机器学习模型样本获取的装置,其特征在于,所述装置包括:
第一模块,用于获取心电信号及其对应的血糖值,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
第二模块,用于将所述心电信号的第一特征值、第二特征值及第三特征值分别输入预设的血压收缩压预测模型和血压舒张压预测模型,获得与所述心电信号对应的血压收缩压预测值和血压舒张压预测值;
第三模块,用于将所述血糖值标识为所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值和所述血压舒张压预测值的标签,并将所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值、所述血压舒张压预测值及所述标签确定为所述机器学习模型样本。
18.根据权利要求17所述的装置,其特征在于,所述第一模块用于获取采集对象的状态信息、心电信号及其对应的血糖值,并确定所述心电信号的第一特征值、第二特征值及第三特征值;
其中,第三模块用于将所述血糖值标识为所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值和所述血压舒张压预测值的标签,并将所述状态信息、所述心电信号的第一特征值、第二特征值、第三特征值、所述血压收缩压预测值、所述血压舒张压预测值及所述标签确定为所述机器学习模型样本。
CN202110995661.2A 2021-08-27 2021-08-27 一种用于血糖预测的机器学习模型样本获取的方法及装置 Pending CN113855007A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110995661.2A CN113855007A (zh) 2021-08-27 2021-08-27 一种用于血糖预测的机器学习模型样本获取的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110995661.2A CN113855007A (zh) 2021-08-27 2021-08-27 一种用于血糖预测的机器学习模型样本获取的方法及装置

Publications (1)

Publication Number Publication Date
CN113855007A true CN113855007A (zh) 2021-12-31

Family

ID=78988606

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110995661.2A Pending CN113855007A (zh) 2021-08-27 2021-08-27 一种用于血糖预测的机器学习模型样本获取的方法及装置

Country Status (1)

Country Link
CN (1) CN113855007A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114420301A (zh) * 2022-01-28 2022-04-29 广东工业大学 基于分段域rf建模预测血糖的方法、系统及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106413534A (zh) * 2015-08-08 2017-02-15 深圳先进技术研究院 连续血压测量装置、测量模型建立方法和系统
CN108420424A (zh) * 2017-02-14 2018-08-21 路提科技股份有限公司 血压测量方法及装置
CN108652640A (zh) * 2018-02-06 2018-10-16 北京大学深圳研究生院 一种基于心电信号的无创血糖检测方法及系统
CN109872821A (zh) * 2019-04-24 2019-06-11 北京信息科技大学 一种构建血糖预测模型的方法及装置
CN113133762A (zh) * 2021-03-03 2021-07-20 刘欣刚 一种无创血糖预测方法及装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106413534A (zh) * 2015-08-08 2017-02-15 深圳先进技术研究院 连续血压测量装置、测量模型建立方法和系统
CN108420424A (zh) * 2017-02-14 2018-08-21 路提科技股份有限公司 血压测量方法及装置
CN108652640A (zh) * 2018-02-06 2018-10-16 北京大学深圳研究生院 一种基于心电信号的无创血糖检测方法及系统
CN109872821A (zh) * 2019-04-24 2019-06-11 北京信息科技大学 一种构建血糖预测模型的方法及装置
CN113133762A (zh) * 2021-03-03 2021-07-20 刘欣刚 一种无创血糖预测方法及装置

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114420301A (zh) * 2022-01-28 2022-04-29 广东工业大学 基于分段域rf建模预测血糖的方法、系统及存储介质
CN114420301B (zh) * 2022-01-28 2022-08-05 广东工业大学 基于分段域rf建模预测血糖的方法、系统及存储介质

Similar Documents

Publication Publication Date Title
US11147500B2 (en) Electrocardiogram processing system for delineation and classification
Braun et al. Demonstration of nonlinear components in heart rate variability of healthy persons
Pathinarupothi et al. Instantaneous heart rate as a robust feature for sleep apnea severity detection using deep learning
US10980429B2 (en) Method and system for cuffless blood pressure estimation using photoplethysmogram features and pulse transit time
US20190059763A1 (en) Method of detecting abnormalities in ecg signals
US20230143594A1 (en) Systems and methods for reduced lead electrocardiogram diagnosis using deep neural networks and rule-based systems
Akhbari et al. ECG segmentation and fiducial point extraction using multi hidden Markov model
US9295397B2 (en) Method and apparatus for beat-space frequency domain prediction of cardiovascular death after acute coronary event
Pimentel et al. Heart beat detection in multimodal physiological data using a hidden semi-Markov model and signal quality indices
Richter et al. A novel approach to propagation pattern analysis in intracardiac atrial fibrillation signals
CN108652640B (zh) 一种基于心电信号的无创血糖检测方法及系统
Rodriguez et al. Detrended fluctuation analysis of heart intrabeat dynamics
Le et al. Multi-module recurrent convolutional neural network with transformer encoder for ECG arrhythmia classification
Antink et al. Detection of heart beats in multimodal data: a robust beat-to-beat interval estimation approach
Sharma et al. QRS complex detection in ECG signals using the synchrosqueezed wavelet transform
Swain et al. Automated detection of myocardial infarction in ECG using modified Stockwell transform and phase distribution pattern from time-frequency analysis
Deng et al. ST-Net: Synthetic ECG tracings for diagnosing various cardiovascular diseases
Liu et al. An efficient neural network-based method for patient-specific information involved arrhythmia detection
CN113855007A (zh) 一种用于血糖预测的机器学习模型样本获取的方法及装置
Zhang et al. Interpretable detection and location of myocardial infarction based on ventricular fusion rule features
Mukhopadhyay et al. ECG feature extraction using differentiation, Hilbert transform, variable threshold and slope reversal approach
CN113854985A (zh) 一种用于血压预测的机器学习模型样本获取的方法及装置
Sörnmo et al. Characterization of f Waves
CN116636856A (zh) 心电分类方法、系统、电子设备及存储介质
Naseri et al. A unified procedure for detecting, quantifying, and validating electrocardiogram T-wave alternans

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