[go: up one dir, main page]

CN113925480B - Coronary heart disease patient bleeding risk assessment method based on machine learning - Google Patents

Coronary heart disease patient bleeding risk assessment method based on machine learning Download PDF

Info

Publication number
CN113925480B
CN113925480B CN202111124787.9A CN202111124787A CN113925480B CN 113925480 B CN113925480 B CN 113925480B CN 202111124787 A CN202111124787 A CN 202111124787A CN 113925480 B CN113925480 B CN 113925480B
Authority
CN
China
Prior art keywords
signal
ppg
bleeding
bleeding risk
risk assessment
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
Application number
CN202111124787.9A
Other languages
Chinese (zh)
Other versions
CN113925480A (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.)
Beijing Hezhong Sizhuang Space Time Material Union Technology Co ltd
Beijing Runmai Technology Co ltd
Chinese PLA General Hospital
Original Assignee
Beijing Hezhong Sizhuang Space Time Material Union Technology Co ltd
Beijing Runmai Technology Co ltd
Chinese PLA General Hospital
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 Beijing Hezhong Sizhuang Space Time Material Union Technology Co ltd, Beijing Runmai Technology Co ltd, Chinese PLA General Hospital filed Critical Beijing Hezhong Sizhuang Space Time Material Union Technology Co ltd
Priority to CN202111124787.9A priority Critical patent/CN113925480B/en
Publication of CN113925480A publication Critical patent/CN113925480A/en
Application granted granted Critical
Publication of CN113925480B publication Critical patent/CN113925480B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/024Measuring pulse rate or heart rate
    • A61B5/02416Measuring pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
    • 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/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • 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/7225Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • 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/7271Specific aspects of physiological measurement analysis
    • A61B5/7275Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Signal Processing (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Physics & Mathematics (AREA)
  • Physiology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Cardiology (AREA)
  • Power Engineering (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

Embodiments of the present disclosure provide a coronary heart disease patient bleeding risk assessment method, apparatus, device and computer-readable storage medium based on machine learning. The method comprises acquiring a PPG measurement signal of a coronary heart disease patient; carrying out segmentation processing on the PPG measurement signal to obtain M signal segments; m is a positive integer greater than or equal to 1; performing signal quality evaluation on the M signal segments through a preset algorithm to determine signal segments with qualified quality; extracting X characteristics from the signal segments with qualified quality, and inputting the characteristics into a bleeding risk evaluation model to obtain a bleeding risk score of a coronary heart disease patient; x is a positive integer of 1 or more. In this way, a bleeding risk assessment for patients with coronary heart disease is achieved.

Description

基于机器学习的冠心病患者出血风险评估方法Bleeding risk assessment method for patients with coronary heart disease based on machine learning

技术领域technical field

本公开的实施例一般涉及信号监测领域,并且更具体地,涉及基于机器学习的冠心病患者出血风险评估方法、装置、设备和计算机可读存储介质。Embodiments of the present disclosure generally relate to the field of signal monitoring, and more specifically, relate to a machine learning-based bleeding risk assessment method, device, device, and computer-readable storage medium for patients with coronary heart disease.

背景技术Background technique

根据世界卫生组织报告,心血管疾病在各种疾病中占比为17.9%,2016年有100万人死于心血管疾病,其中冠状动脉疾病(Coronary Artery Disease冠心病)患者占主要比例。对于冠心病患者,药物治疗(例如抗血小板治疗)是主要的治疗方式。而在冠心病患者进行抗栓治疗期间,出血事件是关注的焦点。出血事件的发生可能导致患者治疗中断,长期失能,甚至死亡。因此对冠心病患者进行出血风险评估研究以发现与预后相关的敏感因素具有重要意义。According to the report of the World Health Organization, cardiovascular diseases accounted for 17.9% of various diseases, and 1 million people died of cardiovascular diseases in 2016, of which patients with coronary artery disease (Coronary Artery Disease) accounted for the main proportion. For patients with coronary heart disease, drug therapy (such as antiplatelet therapy) is the main treatment. During antithrombotic therapy in patients with coronary heart disease, bleeding events are the focus of attention. The occurrence of a bleeding event may lead to treatment interruption, long-term disability, and even death of the patient. Therefore, it is of great significance to conduct bleeding risk assessment research in patients with coronary heart disease to find sensitive factors related to prognosis.

光电容积描记术(Photoplethysmography,PPG)是一种低成本且无创的检测技术,可应用于心血管系统评估。心血管系统状态的变化可以通过PPG波形的形态学特征评估,例如大动脉僵硬度、血容量检测等。Photoplethysmography (Photoplethysmography, PPG) is a low-cost and non-invasive detection technique that can be applied to the assessment of the cardiovascular system. Changes in the state of the cardiovascular system can be assessed by morphological features of the PPG waveform, such as aortic stiffness, blood volume detection, etc.

但目前还没有通过PPG评估冠心病患者出血风险的研究。However, there are no studies evaluating the risk of bleeding in patients with coronary heart disease by PPG.

发明内容Contents of the invention

根据本公开的实施例,提供了一种基于机器学习的冠心病患者出血风险评估方案。According to an embodiment of the present disclosure, a machine learning-based bleeding risk assessment scheme for patients with coronary heart disease is provided.

在本公开的第一方面,提供了一种基于机器学习的冠心病患者出血风险评估方法。该方法包括:In the first aspect of the present disclosure, a machine learning-based bleeding risk assessment method for patients with coronary heart disease is provided. The method includes:

获取冠心病患者的PPG测量信号;Obtaining PPG measurement signals of patients with coronary heart disease;

对所述PPG测量信号进行分段处理,得到N个信号分段;N为大于等于1的正整数;Carry out segment processing to described PPG measurement signal, obtain N signal segments; N is the positive integer greater than or equal to 1;

通过预设算法对所述N个信号分段进行信号质量评估,确定质量合格的信号分段;Carrying out signal quality evaluation on the N signal segments through a preset algorithm, and determining signal segments with qualified quality;

从所述质量合格的信号分段中,提取X个特征,输入到出血风险评估模型中,得到冠心病患者的出血风险评分;X为大于等于1的正整数。X features are extracted from the signal segments with qualified quality and input into the bleeding risk assessment model to obtain the bleeding risk score of patients with coronary heart disease; X is a positive integer greater than or equal to 1.

进一步地,所述对所述PPG测量信号进行分段处理,得到N个信号分段包括:Further, the segmenting the PPG measurement signal to obtain N signal segments includes:

通过巴特沃斯带通滤波器,去除所述PPG测量信号的基线漂移和高频噪音;Remove baseline drift and high-frequency noise of the PPG measurement signal through a Butterworth band-pass filter;

基于所述PPG测量信号的波形,对所述PPG测量信号进行分段处理,得到N个信号分段。Based on the waveform of the PPG measurement signal, segment processing is performed on the PPG measurement signal to obtain N signal segments.

进一步地,所述通过预设算法对所述N个信号分段进行信号质量评估,确定质量合格的信号分段包括:Further, performing signal quality evaluation on the N signal segments through a preset algorithm, and determining a signal segment with qualified quality includes:

对所述N个信号分段中的每一个分段进行重采样;其中,所述重采样长度为所述N个信号分段长度的中位数;Resampling each of the N signal segments; wherein the resampling length is the median of the N signal segment lengths;

将重采样后的分段信号分别代入预设公式,若计算结果大于预设阈值,则当前分段信号为质量合格信号。Substituting the resampled segmented signals into preset formulas, if the calculation result is greater than the preset threshold, the current segmented signal is a qualified signal.

进一步地,所述出血风险评估模型通过如下步骤训练得到:Further, the bleeding risk assessment model is trained through the following steps:

生成训练样本集合,其中,训练样本包括带有标注信息的PPG信号对应的特征向量,其中标注信息为出血状况,出血标注为1,未出血标注为0;Generate a training sample set, wherein the training sample includes a feature vector corresponding to a PPG signal with labeling information, wherein the labeling information is a bleeding condition, and the bleeding is marked as 1, and the non-bleeding is marked as 0;

利用所述训练样本集合中的样本对出血风险评估模型进行训练,以PPG信号对应的特征向量为输入,以出血状况为输出,当输出的出血状况与标注的出血状况的统一率满足预设阈值时,完成对出血风险评估模型的训练。Using the samples in the training sample set to train the bleeding risk assessment model, the feature vector corresponding to the PPG signal is used as input, and the bleeding condition is output, when the unity rate of the output bleeding condition and the marked bleeding condition meets a preset threshold , complete the training of the bleeding risk assessment model.

进一步地,所述PPG信号对应的特征向量包括:Further, the eigenvector corresponding to the PPG signal includes:

对采集的PPG信号进行分析,分别从时域、频域和小波包分解中提取与所述PPG信号对应的特征向量。The collected PPG signal is analyzed, and feature vectors corresponding to the PPG signal are extracted from time domain, frequency domain and wavelet packet decomposition respectively.

进一步地,所述对采集的PPG信号进行分析包括:Further, said analyzing the collected PPG signal includes:

对所述PPG信号、PPG信号的一阶导数和PPG信号的二阶导数进行分析。The PPG signal, the first derivative of the PPG signal and the second derivative of the PPG signal are analyzed.

进一步地,采用XGBoost算法对所述出血风险评估模型进行训练。Further, the XGBoost algorithm is used to train the bleeding risk assessment model.

在本公开的第二方面,提供了一种基于机器学习的冠心病患者出血风险评估装置。该装置包括:In the second aspect of the present disclosure, a device for assessing bleeding risk in patients with coronary heart disease based on machine learning is provided. The unit includes:

获取模块,用于获取冠心病患者的PPG测量信号;An acquisition module, configured to acquire PPG measurement signals of patients with coronary heart disease;

处理模块,用于对所述PPG测量信号进行分段处理,得到N个信号分段;N为大于等于1的正整数;A processing module, configured to perform segmentation processing on the PPG measurement signal to obtain N signal segments; N is a positive integer greater than or equal to 1;

评估模块,用于通过预设算法对所述N个信号分段进行信号质量评估,确定质量合格的信号分段;An evaluation module, configured to evaluate the signal quality of the N signal segments through a preset algorithm, and determine a signal segment with qualified quality;

评分模块,用于从所述质量合格的信号分段中,提取X个特征,输入到出血风险评估模型中,得到出血风险评分;X为大于等于1的正整数。The scoring module is used to extract X features from the qualified signal segment and input them into the bleeding risk assessment model to obtain a bleeding risk score; X is a positive integer greater than or equal to 1.

在本公开的第三方面,提供了一种电子设备。该电子设备包括:存储器和处理器,所述存储器上存储有计算机程序,所述处理器执行所述程序时实现如以上所述的方法。In a third aspect of the present disclosure, an electronic device is provided. The electronic device includes: a memory and a processor, where a computer program is stored in the memory, and the processor implements the method as described above when executing the program.

在本公开的第四方面,提供了一种计算机可读存储介质,其上存储有计算机程序,所述程序被处理器执行时实现如根据本公开的第一方面的方法。In a fourth aspect of the present disclosure, there is provided a computer-readable storage medium on which a computer program is stored, and when the program is executed by a processor, the method according to the first aspect of the present disclosure is implemented.

本申请实施例提供的基于机器学习的冠心病患者出血风险评估方法,通过获取PPG测量信号;对所述PPG测量信号进行分段处理,得到N个信号分段;N为大于等于1的正整数;通过预设算法对所述N个信号分段进行信号质量评估,确定质量合格的信号分段;从所述质量合格的信号分段中,提取X个时域特征,输入到出血风险评估模型中,得到冠心病患者的出血风险评分;X为大于等于1的正整数,实现了对CAD患者的出血风险评估。The bleeding risk assessment method for patients with coronary heart disease based on machine learning provided in the embodiment of the present application obtains the PPG measurement signal; performs segmentation processing on the PPG measurement signal to obtain N signal segments; N is a positive integer greater than or equal to 1 ; Carry out signal quality assessment on the N signal segments by a preset algorithm, and determine the qualified signal segments; extract X time-domain features from the qualified signal segments, and input them into the bleeding risk assessment model Among them, the bleeding risk score of patients with coronary heart disease is obtained; X is a positive integer greater than or equal to 1, which realizes the bleeding risk assessment of CAD patients.

应当理解,发明内容部分中所描述的内容并非旨在限定本公开的实施例的关键或重要特征,亦非用于限制本公开的范围。本公开的其它特征将通过以下的描述变得容易理解。It should be understood that what is described in the Summary of the Invention is not intended to limit the key or important features of the embodiments of the present disclosure, nor is it intended to limit the scope of the present disclosure. Other features of the present disclosure will be readily understood through the following description.

附图说明Description of drawings

结合附图并参考以下详细说明,本公开各实施例的上述和其他特征、优点及方面将变得更加明显。在附图中,相同或相似的附图标记表示相同或相似的元素,其中:The above and other features, advantages and aspects of the various embodiments of the present disclosure will become more apparent with reference to the following detailed description when taken in conjunction with the accompanying drawings. In the drawings, identical or similar reference numerals denote identical or similar elements, wherein:

图1示出了能够在其中实现本公开的实施例的示例性运行环境的示意图;FIG. 1 shows a schematic diagram of an exemplary operating environment in which embodiments of the present disclosure can be implemented;

图2示出了根据本公开的实施例的基于机器学习的冠心病患者出血风险评估方法的流程图;2 shows a flow chart of a machine learning-based bleeding risk assessment method for patients with coronary heart disease according to an embodiment of the present disclosure;

图3示出了根据本公开的实施例的PPG、VPG和APG信号示意图;Fig. 3 shows a schematic diagram of PPG, VPG and APG signals according to an embodiment of the present disclosure;

图4示出了根据本公开的实施例的XGBoost模型的10倍交叉验证ROC曲线、平均ROC曲线示意图;4 shows a schematic diagram of a 10-fold cross-validation ROC curve and an average ROC curve of an XGBoost model according to an embodiment of the present disclosure;

图5示出了根据本公开的实施例的SHAP框架对XGBoost模型进行特性分析的示意图;FIG. 5 shows a schematic diagram of analyzing the characteristics of the XGBoost model by the SHAP framework according to an embodiment of the present disclosure;

图6示出了根据本公开的实施例的基于机器学习的冠心病患者出血风险评估装置的方框图;FIG. 6 shows a block diagram of a device for assessing bleeding risk in patients with coronary heart disease based on machine learning according to an embodiment of the present disclosure;

图7示出了能够实施本公开的实施例的示例性电子设备的方框图。FIG. 7 shows a block diagram of an exemplary electronic device capable of implementing embodiments of the present disclosure.

具体实施方式Detailed ways

为使本公开实施例的目的、技术方案和优点更加清楚,下面将结合本公开实施例中的附图,对本公开实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本公开一部分实施例,而不是全部的实施例。基于本公开中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的全部其他实施例,都属于本公开保护的范围。In order to make the purpose, technical solutions and advantages of the embodiments of the present disclosure clearer, the technical solutions in the embodiments of the present disclosure will be clearly and completely described below in conjunction with the drawings in the embodiments of the present disclosure. Obviously, the described embodiments It is a part of the embodiments of the present disclosure, but not all of them. Based on the embodiments in the present disclosure, all other embodiments obtained by persons of ordinary skill in the art without creative efforts fall within the protection scope of the present disclosure.

另外,本文中术语“和/或”,仅仅是一种描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。另外,本文中字符“/”,一般表示前后关联对象是一种“或”的关系。In addition, the term "and/or" in this article is only an association relationship describing associated objects, which means that there may be three relationships, for example, A and/or B, which may mean: A exists alone, A and B exist at the same time, There are three cases of B alone. In addition, the character "/" in this article generally indicates that the contextual objects are an "or" relationship.

图1示出了可以应用本申请的基于机器学习的冠心病患者出血风险评估方法装置的实施例的示例性系统架构100。FIG. 1 shows an exemplary system architecture 100 of an embodiment of the method and apparatus for assessing bleeding risk in patients with coronary heart disease based on machine learning of the present application.

如图1所示,系统架构100可以包括终端设备101、102、103,网络104和服务器105。网络104用以在终端设备101、102、103和服务器105之间提供通信链路的介质。网络104可以包括各种连接类型,例如有线、无线通信链路或者光纤电缆等等。As shown in FIG. 1 , a system architecture 100 may include terminal devices 101 , 102 , 103 , a network 104 and a server 105 . The network 104 is used as a medium for providing communication links between the terminal devices 101 , 102 , 103 and the server 105 . Network 104 may include various connection types, such as wires, wireless communication links, or fiber optic cables, among others.

用户可以使用终端设备101、102、103通过网络104与服务器105交互,以接收或发送消息等。终端设备101、102、103上可以安装有各种通讯客户端应用,例如模型训练类应用、视频识别类应用、网页浏览器应用、社交平台软件等。Users can use terminal devices 101 , 102 , 103 to interact with server 105 via network 104 to receive or send messages and the like. Various communication client applications, such as model training applications, video recognition applications, web browser applications, and social platform software, can be installed on the terminal devices 101, 102, and 103.

终端设备101、102、103可以是硬件,也可以是软件。当终端设备101、102、103为硬件时,可以是具有显示屏的各种电子设备,包括但不限于智能手机、平板电脑、电子书阅读器、MP3播放器(Moving Picture Experts Group Audio Layer III,动态影像专家压缩标准音频层面3)、MP4(Moving Picture Experts Group Audio Layer IV,动态影像专家压缩标准音频层面4)播放器、PPG测量设备、膝上型便携计算机和台式计算机等等。当终端设备101、102、103为软件时,可以安装在上述所列举的电子设备中。其可以实现成多个软件或软件模块(例如用来提供分布式服务的多个软件或软件模块),也可以实现成单个软件或软件模块。在此不做具体限定。The terminal devices 101, 102, and 103 may be hardware or software. When the terminal devices 101, 102, 103 are hardware, they can be various electronic devices with display screens, including but not limited to smart phones, tablet computers, e-book readers, MP3 players (Moving Picture Experts Group Audio Layer III, Moving Picture Experts Compression Standard Audio Layer 3), MP4 (Moving Picture Experts Group Audio Layer IV, Moving Picture Experts Group Audio Layer 4) Players, PPG Measuring Devices, Laptop Portable Computers and Desktop Computers, etc. When the terminal devices 101, 102, 103 are software, they can be installed in the electronic devices listed above. It may be implemented as multiple software or software modules (for example, multiple software or software modules for providing distributed services), or as a single software or software module. No specific limitation is made here.

当终端101、102、103为硬件时,其上还可以安装有视频采集设备。视频采集设备可以是各种能实现采集视频功能的设备,如摄像头、传感器等等。用户可以利用终端101、102、103上的视频采集设备来采集视频。When the terminals 101, 102, and 103 are hardware, video acquisition equipment may also be installed on them. The video capture device may be various devices capable of capturing video, such as a camera, a sensor, and the like. Users can use the video capture devices on the terminals 101, 102, 103 to capture videos.

服务器105可以是提供各种服务的服务器,例如对终端设备101、102、103上显示的数据处理的后台服务器。后台服务器可以对接收到的数据进行分析等处理,并可以将处理结果(例如评估结果)反馈给终端设备。The server 105 may be a server that provides various services, such as a background server that processes data displayed on the terminal devices 101 , 102 , 103 . The background server can perform analysis and other processing on the received data, and can feed back the processing results (eg evaluation results) to the terminal device.

需要说明的是,服务器可以是硬件,也可以是软件。当服务器为硬件时,可以实现成多个服务器组成的分布式服务器集群,也可以实现成单个服务器。当服务器为软件时,可以实现成多个软件或软件模块(例如用来提供分布式服务的多个软件或软件模块),也可以实现成单个软件或软件模块。在此不做具体限定。It should be noted that the server may be hardware or software. When the server is hardware, it can be implemented as a distributed server cluster composed of multiple servers, or as a single server. When the server is software, it can be implemented as multiple software or software modules (for example, multiple software or software modules for providing distributed services), or can be implemented as a single software or software module. No specific limitation is made here.

应该理解,图1中的终端设备、网络和服务器的数目仅仅是示意性的。根据实现需要,可以具有任意数目的终端设备、网络和服务器。特别地,在目标数据不需要从远程获取的情况下,上述系统架构可以不包括网络,而只包括终端设备或服务器。It should be understood that the numbers of terminal devices, networks and servers in Fig. 1 are only illustrative. According to the implementation needs, there can be any number of terminal devices, networks and servers. In particular, in the case where the target data does not need to be obtained from a remote location, the above-mentioned system architecture may not include a network, but only include terminal devices or servers.

如图2所示,是本申请实施例基于机器学习的冠心病患者出血风险评估方法的流程图。从图2中可以看出,本实施例的基于机器学习的冠心病患者出血风险评估方法,包括以下步骤:As shown in FIG. 2 , it is a flow chart of a method for assessing bleeding risk in patients with coronary heart disease based on machine learning according to an embodiment of the present application. As can be seen from Fig. 2, the bleeding risk assessment method for patients with coronary heart disease based on machine learning in this embodiment includes the following steps:

S210,获取冠心病患者的PPG测量信号。S210, acquiring a PPG measurement signal of a patient with coronary heart disease.

在一些实施例中,用于基于机器学习的冠心病患者出血风险评估方法的执行主体(例如图1所示的服务器)可以通过有线方式或者无线连接的方式获取PPG测量信号。In some embodiments, the executive body (such as the server shown in FIG. 1 ) used for the machine learning-based bleeding risk assessment method for patients with coronary heart disease can acquire the PPG measurement signal through wired or wireless connection.

在一些实施例中,上述执行主体可以获取与之通信连接的电子设备(例如图1所示的终端设备)发送的PPG测量信号;其中,所述电子设备上装有PPG传感器,用于获取用户(患者)的PPG数据,例如,光学心率传感器等。In some embodiments, the above-mentioned executive body can obtain the PPG measurement signal sent by the electronic device (such as the terminal device shown in FIG. patient) PPG data, for example, optical heart rate sensor, etc.

在一些实施例中,当PPG测量设备不能单独进行上传数据时,如,部分指端PPG测量设备,PPG测量设备也可将测量结果发送至与其连接的移动设备(手机等),通过所述移动设备将相关的PPG测量信息上传至如图1所示的服务器。In some embodiments, when the PPG measurement device cannot upload data alone, such as some finger-end PPG measurement devices, the PPG measurement device can also send the measurement results to the mobile device (mobile phone, etc.) connected to it, through the mobile The device uploads the relevant PPG measurement information to the server shown in Figure 1.

在一些实施例中,为了测量的准确性,可通过多种光强(发光强度)采集PPG测量数据;采集时间通常为20s-60s;可根据实际应用场景进行设定,对此不做进一步限定。In some embodiments, for the accuracy of measurement, PPG measurement data can be collected through various light intensities (luminous intensity); the collection time is usually 20s-60s; it can be set according to the actual application scene, and no further limitation is made on this .

在一些实施例中,所述PPG的信号波形如图3所示,图3中VPG为PPG的一阶导数、APG为PPG的二阶导数;T和Y表示时间ms和相应点的幅度。In some embodiments, the signal waveform of the PPG is as shown in Figure 3, where VPG is the first derivative of PPG and APG is the second derivative of PPG in Figure 3; T and Y represent the time ms and the amplitude of the corresponding point.

需要说明的是,本公开中的PPG测量设备(电子设备)通常为,装有PPG传感器的便携设备,可在居家环境中进行使用。It should be noted that the PPG measuring device (electronic device) in the present disclosure is generally a portable device equipped with a PPG sensor, which can be used in a home environment.

S220,对所述PPG测量信号进行分段处理,得到M个信号分段;M为大于等于1的正整数。S220. Perform segmentation processing on the PPG measurement signal to obtain M signal segments; M is a positive integer greater than or equal to 1.

在一些实施例中,对所述PPG测量信号进行分段处理前,需要对所述PPG测量信号进行预处理;In some embodiments, before segmenting the PPG measurement signal, the PPG measurement signal needs to be preprocessed;

具体的,可通过巴特沃斯带通滤波器,设置截止频率为0.2和20Hz去除所述PPG测量信号的基线漂移和高频噪音。Specifically, the baseline drift and high-frequency noise of the PPG measurement signal can be removed by using a Butterworth band-pass filter with cut-off frequencies set at 0.2 and 20 Hz.

在一些实施例中,按照PPG测量信号的波形周期,对预处理后的PPG测量信号进行切分,得到M个信号分段,可记为S={S1,S2,...,SM};其中M为大于等于1的正整数;S表示信号分段的集合。In some embodiments, according to the waveform period of the PPG measurement signal, the preprocessed PPG measurement signal is segmented to obtain M signal segments, which can be recorded as S={S1, S2,...,SM}; Where M is a positive integer greater than or equal to 1; S represents a set of signal segments.

S230,通过预设算法对所述M个信号分段进行信号质量评估,确定质量合格的信号分段。S230. Perform signal quality evaluation on the M signal segments by using a preset algorithm, and determine signal segments with qualified quality.

考虑到用户在测量过程中可能产生的噪声,S中的一些分段可能会被损坏,因此,需要对每个分段的信号质量进行检测,即,对每个分段的信号质量进行评估。Considering the noise that may be generated by the user during the measurement process, some segments in S may be damaged. Therefore, it is necessary to detect the signal quality of each segment, that is, to evaluate the signal quality of each segment.

在一些实施例中,对所述M个信号分段中的每一个分段进行重采样,即,对集合S中每个分段进行重采样为相同长度,可记为RS={RS1,RS2,…,RSM};其中,采样长度为所有分段长度的中位数;In some embodiments, each of the M signal segments is resampled, that is, each segment in the set S is resampled to have the same length, which can be written as RS={RS 1 , RS 2 ,…,RS M }; Among them, the sampling length is the median of all segment lengths;

将S中的分段信号分别代入如下公式,若满足计算条件,则为质量合格的分段信号;Substitute the segmented signals in S into the following formulas respectively, if the calculation conditions are met, it is a segmented signal with qualified quality;

Svalid={Si∣Si∈S,RSi∈RS,r(RSi,RST)>0.9}S valid ={S i ∣S i ∈S,RS i ∈RS,r(RS i ,RS T )>0.9}

其中,所述RST为,所述RS的平均值;Wherein, the RS T is the average value of the RS;

所述r为,皮尔逊相关系数,取值范围为-1~1;The r is the Pearson correlation coefficient, and the value range is -1 to 1;

所述0.9为预设阈值,根据大量实验得出;所述预设阈值也可根据人工经验、大数据分析和/或实际应用场景进行设定。The 0.9 is a preset threshold, which is obtained from a large number of experiments; the preset threshold can also be set according to manual experience, big data analysis and/or actual application scenarios.

S240,从所述质量合格的信号分段中,提取X个特征,输入到出血风险评估模型中,得到冠心病患者的出血风险评分;X为大于等于1的正整数。S240, extracting X features from the qualified signal segment and inputting them into the bleeding risk assessment model to obtain the bleeding risk score of the coronary heart disease patient; X is a positive integer greater than or equal to 1.

在一些实施例中,从所述质量合格的信号分段中,提取出X个特征,如提取30个特征,参考表1;T和Y表示时间ms和相应点的幅度,参考图3。In some embodiments, X features are extracted from the qualified signal segment, for example, 30 features are extracted, refer to Table 1; T and Y represent the time ms and the amplitude of the corresponding point, refer to FIG. 3 .

Figure BDA0003278410850000081
Figure BDA0003278410850000081

Figure BDA0003278410850000091
Figure BDA0003278410850000091

表1Table 1

在一些实施例中,为了测量的准确性,每个波形(信号分段)在提取特征前,需要进行归一化处理,将其归一化为0~1的范围内。In some embodiments, for measurement accuracy, each waveform (signal segment) needs to be normalized before feature extraction, and normalized to a range of 0-1.

在一些实施例中,所述出血风险评估模型可通过如下步骤训练得到:In some embodiments, the bleeding risk assessment model can be trained through the following steps:

生成训练样本集合,其中,训练样本包括带有标注信息的PPG信号对应的特征向量,其中标注信息为出血状况,出血标注为1,未出血标注为0;Generate a training sample set, wherein the training sample includes a feature vector corresponding to a PPG signal with labeling information, wherein the labeling information is a bleeding condition, and the bleeding is marked as 1, and the non-bleeding is marked as 0;

利用所述训练样本集合中的样本对出血风险评估模型进行训练,以PPG信号对应的特征向量为输入,以出血状况为输出,当输出的出血状况与标注的出血状况的统一率满足预设阈值时,完成对出血风险评估模型的训练;所述阈值可跟实际应用场景进行设置;Using the samples in the training sample set to train the bleeding risk assessment model, the feature vector corresponding to the PPG signal is used as input, and the bleeding condition is output, when the unity rate of the output bleeding condition and the marked bleeding condition meets a preset threshold , complete the training of the bleeding risk assessment model; the threshold can be set according to the actual application scenario;

进一步地,基于所述输出,可通过人工经验、大数据分析和/或设置阈值的方法,计算出冠心病患者的出血风险评分,评分越高,则出血风险越大;Further, based on the output, the bleeding risk score of patients with coronary heart disease can be calculated through manual experience, big data analysis and/or threshold setting methods, the higher the score, the greater the bleeding risk;

其中,所述PPG信号对应的特征向量,还包括PPG的一阶导数(VPG)和二阶导数(APG)的特征向量;同步的PPG、VPG、APG信号如图3所示;Wherein, the eigenvector corresponding to the PPG signal also includes the eigenvector of the first order derivative (VPG) and the second order derivative (APG) of PPG; the synchronous PPG, VPG, APG signals are as shown in Figure 3;

所述PPG信号可以为通过大数据采集的PPG信号,如,两年内大规模采集的PPG信号;The PPG signal may be a PPG signal collected through big data, such as a PPG signal collected on a large scale within two years;

对采集PPG信号进行分析,得到PPG信号对应的特征向量,即,对所述PPG信号、PPG信号的一阶导数和PPG信号的二阶导数进行分析,从时域、频域和小波包分解(能量特征)中提取30维特征向量(可根据实际应用场景进行设定),参考表1;Analyze the collected PPG signal to obtain the corresponding eigenvector of the PPG signal, that is, analyze the PPG signal, the first-order derivative of the PPG signal and the second-order derivative of the PPG signal, and decompose ( Extract 30-dimensional feature vectors from energy features) (can be set according to actual application scenarios), refer to Table 1;

其中,所述频域可通过如下方法进行确定:Wherein, the frequency domain can be determined by the following method:

PPG信号由丰富的频率分量组成,可采用Welch算法计算PPG信号的功率谱密度,通过检测极值点来确定每个谐波的位置,将第二个到第六个谐波除以基频(一次谐波)的功率归一化处理,确定频域特征,即,表1中的频域H1~H5;The PPG signal is composed of rich frequency components. Welch algorithm can be used to calculate the power spectral density of the PPG signal, and the position of each harmonic can be determined by detecting the extreme points, and the second to sixth harmonics can be divided by the fundamental frequency ( The power normalization processing of the first harmonic) determines the frequency domain characteristics, that is, the frequency domain H1~H5 in Table 1;

通过大量的研究实验,PPG信号中99%的能量集中在1~10Hz范围内,因此,在本公开中,仅对10Hz内的频率分量进行研究、提取;Through a large number of research experiments, 99% of the energy in the PPG signal is concentrated in the range of 1-10Hz. Therefore, in this disclosure, only the frequency components within 10Hz are studied and extracted;

具体地,PPG信号采样率为500Hz,经过8次小波包分解后,每个子带的带宽为0.977Hz,保留10个分量为E1~E10,PPG信号的总能量表示为Eall。Specifically, the sampling rate of the PPG signal is 500 Hz. After 8 times of wavelet packet decomposition, the bandwidth of each subband is 0.977 Hz, and 10 components are reserved as E1-E10. The total energy of the PPG signal is expressed as Eall.

在本公开中,对采集的1683为冠心病患者的数据进行采集、分析,其中有114为患者至少有一次阳性事件(出血),将没有人口统计学的记录和信号质量较差的记录删除,得到如表2所示的患者人口学特征;In this disclosure, the collected data of 1683 patients with coronary heart disease were collected and analyzed, among which 114 patients had at least one positive event (bleeding), and the records without demographics and poor signal quality were deleted. Obtain the demographic characteristics of patients as shown in Table 2;

Figure BDA0003278410850000101
Figure BDA0003278410850000101

表2Table 2

其中,阴性组和阳性组间有10个特征具有显著的统计学差异,参见表3:Among them, there are 10 characteristics with significant statistical differences between the negative group and the positive group, see Table 3:

Figure BDA0003278410850000102
Figure BDA0003278410850000102

Figure BDA0003278410850000111
Figure BDA0003278410850000111

表3table 3

参见表3,阳性组Td小于阴性组,而RI较大,从PPG波形的几何特征角度进行分析,有可能是阳性组的舒张波延迟,导致舒张波宽度增加,Td、Rarea趋于变小;频域的H1~H5在阴性组和阳性组之间有统计学差异,归一化5个谐波逐渐减小;小波包分解计算的总能量Eall在阳性组中较小。See Table 3, the Td of the positive group is smaller than that of the negative group, but the RI is larger. From the perspective of the geometric characteristics of the PPG waveform, it is possible that the diastolic wave delay in the positive group increases the width of the diastolic wave, and Td and Rarea tend to decrease; H1-H5 in the frequency domain had statistical differences between the negative group and the positive group, and the normalized 5 harmonics gradually decreased; the total energy Eall calculated by wavelet packet decomposition was smaller in the positive group.

在一些实施例中,基于采集的冠心病患者的PPG数据,分别分别使用LR(LogisticRegression)、SVR(Support Vector Regression)、RF(Random Forest)、XGBoost算法,训练获得出血分析评估模型。90%的数据用于训练分类模型(训练集),其余10%的数据用于测试(测试集),经过10折交叉验证与网格搜索进行训练评估,网格搜索获得的最优超参数用于10折交叉验证,得到各模型平均AUC(模型评估指标),如表4所示。各模型的敏感度与特异度表现如表5所示。经过对比,多个算法模型中XGBoost表现最优,平均AUC为0.762,敏感度与特异度分别为0.679与0.714,显著高于其他模型。XGBoost模型的10折交叉验证ROC曲线、平均ROC曲线和AUC如图4所示。In some embodiments, based on the collected PPG data of patients with coronary heart disease, LR (Logistic Regression), SVR (Support Vector Regression), RF (Random Forest), and XGBoost algorithms are respectively used to train and obtain a bleeding analysis and evaluation model. 90% of the data is used to train the classification model (training set), and the remaining 10% of the data is used for testing (test set). After 10-fold cross-validation and grid search for training and evaluation, the optimal hyperparameters obtained by grid search are used After 10-fold cross-validation, the average AUC (model evaluation index) of each model is obtained, as shown in Table 4. The sensitivity and specificity performance of each model is shown in Table 5. After comparison, XGBoost performed the best among multiple algorithm models, with an average AUC of 0.762, sensitivity and specificity of 0.679 and 0.714, which were significantly higher than other models. The 10-fold cross-validation ROC curve, average ROC curve and AUC of the XGBoost model are shown in Figure 4.

Figure BDA0003278410850000112
Figure BDA0003278410850000112

表4Table 4

Figure BDA0003278410850000113
Figure BDA0003278410850000113

Figure BDA0003278410850000121
Figure BDA0003278410850000121

表5table 5

进一步地,使用SHAP(SHapley Additive interpretation)框架对XGBoost模型的特性进行分析,如图5所示,显示了其中的20个特征;SHAP值表示特征对目标的贡献(负:0,正:1)。特征值用不同的颜色表示,即特征值越大表示颜色越红,反之特征值越小表示颜色越蓝(图5中,灰度值较低的为蓝色,灰度值较高的为红色);例如,随着H4、E10、SI等特征值的增加,SHAP值趋于小于0;即,模型更倾向于将当前样本确定为目标0;相反,随着RI、VDae等值的减小,SHAP值趋于大于0。Further, the SHAP (SHapley Additive interpretation) framework is used to analyze the characteristics of the XGBoost model, as shown in Figure 5, which shows 20 features; the SHAP value indicates the contribution of the feature to the target (negative: 0, positive: 1) . The eigenvalues are represented by different colors, that is, the larger the eigenvalue, the redder the color, and the smaller the eigenvalue, the bluer the color (in Figure 5, the lower gray value is blue, and the higher gray value is red ); for example, with the increase of H4, E10, SI and other eigenvalues, the SHAP value tends to be less than 0; that is, the model is more inclined to determine the current sample as the target 0; on the contrary, as the RI, VDae and other values decrease , the SHAP value tends to be greater than 0.

此外,SHAP值的绝对值之和反映了该特征的重要性,因此,根据这一准则,可认为H4是最重要的特征。上述观察结果与表3中的结果也有很好的一致性;例如,对于阴性组,Td、Rarea、H4、H1、H5、H2、Eall一致,即当特征值增加时,SHAP值趋于负,而RI值的表现恰恰相反。In addition, the sum of the absolute values of the SHAP values reflects the importance of this feature, therefore, according to this criterion, H4 can be considered as the most important feature. The above observations are also in good agreement with the results in Table 3; for example, for the negative group, Td, Rarea, H4, H1, H5, H2, Eall are consistent, that is, when the eigenvalue increases, the SHAP value tends to be negative, The RI value behaves just the opposite.

综上所述,在本公开中可采用XGBoost模型训练所述出血风险评估模型;采用SHAP框架对所述出血风险评估模型进行描述。In summary, in this disclosure, the XGBoost model can be used to train the bleeding risk assessment model; the SHAP framework is used to describe the bleeding risk assessment model.

在一些实施例中,将提取的特征输入至训练好的出血风险评估模型,确定当前患者的出血风险评分;In some embodiments, the extracted features are input into the trained bleeding risk assessment model to determine the bleeding risk score of the current patient;

进一步地,将所述出血风险评分发送至PPG测量设备和/或如图1中所示的终端设备,形成闭环,帮助患者实现自我健康管理。Further, the bleeding risk score is sent to the PPG measurement device and/or the terminal device as shown in FIG. 1 to form a closed loop and help patients realize self-health management.

根据本公开的实施例,实现了以下技术效果:According to the embodiments of the present disclosure, the following technical effects are achieved:

本发明区别于其他冠心病患者出血风险评估研究,使用便携式PPG设备,不依赖于临床环境,用户在居家环境中即可使用设备评估可能发生的出血事件。系统可远程帮助提醒CAD患者要给予足够的重视或协助医生选择治疗方案。The present invention is different from other bleeding risk assessment studies for patients with coronary heart disease. The portable PPG device is used independently of the clinical environment, and the user can use the device to evaluate possible bleeding events in the home environment. The system can remotely help remind CAD patients to pay enough attention or assist doctors to choose a treatment plan.

需要说明的是,对于前述的各方法实施例,为了简单描述,故将其都表述为一系列的动作组合,但是本领域技术人员应该知悉,本公开并不受所描述的动作顺序的限制,因为依据本公开,某些步骤可以采用其他顺序或者同时进行。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于可选实施例,所涉及的动作和模块并不一定是本公开所必须的。It should be noted that for the foregoing method embodiments, for the sake of simple description, they are all expressed as a series of action combinations, but those skilled in the art should know that the present disclosure is not limited by the described action sequence. Because of this disclosure, certain steps may be performed in other orders or simultaneously. Secondly, those skilled in the art should also know that the embodiments described in the specification are all optional embodiments, and the actions and modules involved are not necessarily required by the present disclosure.

以上是关于方法实施例的介绍,以下通过装置实施例,对本公开所述方案进行进一步说明。The above is the introduction of the method embodiments, and the solution of the present disclosure will be further described through the device embodiments below.

图6示出了根据本公开的实施例的基于机器学习的冠心病患者出血风险评估装置600的方框图。如图6所示,装置600包括:Fig. 6 shows a block diagram of an apparatus 600 for assessing bleeding risk of patients with coronary heart disease based on machine learning according to an embodiment of the present disclosure. As shown in Figure 6, the device 600 includes:

获取模块610,用于获取冠心病患者的PPG测量信号;An acquisition module 610, configured to acquire a PPG measurement signal of a patient with coronary heart disease;

处理模块620,用于对所述PPG测量信号进行分段处理,得到M个信号分段;M为大于等于1的正整数;The processing module 620 is configured to perform segmentation processing on the PPG measurement signal to obtain M signal segments; M is a positive integer greater than or equal to 1;

评估模块630,用于通过预设算法对所述M个信号分段进行信号质量评估,确定质量合格的信号分段;An evaluation module 630, configured to evaluate the signal quality of the M signal segments through a preset algorithm, and determine a signal segment with qualified quality;

评分模块640,用于从所述质量合格的信号分段中,提取X个特征,输入到出血风险评估模型中,得到出血风险评分;X为大于等于1的正整数。The scoring module 640 is configured to extract X features from the qualified signal segment and input them into the bleeding risk assessment model to obtain a bleeding risk score; X is a positive integer greater than or equal to 1.

所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,所述描述的模块的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。Those skilled in the art can clearly understand that for the convenience and brevity of the description, the specific working process of the described modules can refer to the corresponding process in the foregoing method embodiments, which will not be repeated here.

图7示出了可以用来实施本公开的实施例的电子设备700的示意性框图。设备700可以用于实现图1的消息系统104和消息到达率确定系统106中的至少一个。如图所示,设备700包括中央处理单元(CPU)701,其可以根据存储在只读存储器(ROM)702中的计算机程序指令或者从存储单元708加载到随机访问存储器(RAM)703中的计算机程序指令,来执行各种适当的动作和处理。在RAM 703中,还可以存储设备700操作所需的各种程序和数据。CPU701、ROM 702以及RAM 703通过总线704彼此相连。输入/输出(I/O)接口705也连接至总线704。FIG. 7 shows a schematic block diagram of an electronic device 700 that can be used to implement embodiments of the present disclosure. The device 700 may be used to implement at least one of the message system 104 and the message arrival rate determination system 106 of FIG. 1 . As shown, the device 700 includes a central processing unit (CPU) 701 that can be programmed according to computer program instructions stored in a read-only memory (ROM) 702 or loaded from a storage unit 708 into a random-access memory (RAM) 703 program instructions to perform various appropriate actions and processes. In the RAM 703, various programs and data necessary for the operation of the device 700 can also be stored. The CPU 701 , ROM 702 , and RAM 703 are connected to each other via a bus 704 . An input/output (I/O) interface 705 is also connected to the bus 704 .

设备700中的多个部件连接至I/O接口705,包括:输入单元706,例如键盘、鼠标等;输出单元707,例如各种类型的显示器、扬声器等;存储单元708,例如磁盘、光盘等;以及通信单元709,例如网卡、调制解调器、无线通信收发机等。通信单元709允许设备700通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。Multiple components in the device 700 are connected to the I/O interface 705, including: an input unit 706, such as a keyboard, a mouse, etc.; an output unit 707, such as various types of displays, speakers, etc.; a storage unit 708, such as a magnetic disk, an optical disk, etc. ; and a communication unit 709, such as a network card, a modem, a wireless communication transceiver, and the like. The communication unit 709 allows the device 700 to exchange information/data with other devices over a computer network such as the Internet and/or various telecommunication networks.

处理单元701执行上文所描述的各个方法和处理。例如,在一些实施例中,方法可被实现为计算机软件程序,其被有形地包含于机器可读介质,例如存储单元708。在一些实施例中,计算机程序的部分或者全部可以经由ROM 702和/或通信单元709而被载入和/或安装到设备700上。当计算机程序加载到RAM 703并由CPU 701执行时,可以执行上文描述的方法的一个或多个步骤。备选地,在其他实施例中,CPU 701可以通过其他任何适当的方式(例如,借助于固件)而被配置为执行方法。The processing unit 701 executes the various methods and processes described above. For example, in some embodiments, the methods may be implemented as a computer software program tangibly embodied on a machine-readable medium, such as storage unit 708 . In some embodiments, part or all of the computer program may be loaded and/or installed on the device 700 via the ROM 702 and/or the communication unit 709 . When a computer program is loaded into RAM 703 and executed by CPU 701, one or more steps of the methods described above may be performed. Alternatively, in other embodiments, the CPU 701 may be configured to execute the method in any other suitable manner (eg, by means of firmware).

本文中以上描述的功能可以至少部分地由一个或多个硬件逻辑部件来执行。例如,非限制性地,可以使用的示范类型的硬件逻辑部件包括:场可编程门阵列(FPGA)、专用集成电路(ASIC)、专用标准产品(ASSP)、芯片上系统的系统(SOC)、负载可编程逻辑设备(CPLD)等等。The functions described herein above may be performed at least in part by one or more hardware logic components. For example, without limitation, exemplary types of hardware logic components that may be used include: field programmable gate array (FPGA), application specific integrated circuit (ASIC), application specific standard product (ASSP), system on a chip (SOC), load programmable logic device (CPLD), etc.

用于实施本公开的方法的程序代码可以采用一个或多个编程语言的任何组合来编写。这些程序代码可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器或控制器,使得程序代码当由处理器或控制器执行时使流程图和/或框图中所规定的功能/操作被实施。程序代码可以完全在机器上执行、部分地在机器上执行,作为独立软件包部分地在机器上执行且部分地在远程机器上执行或完全在远程机器或服务器上执行。Program codes for implementing the methods of the present disclosure may be written in any combination of one or more programming languages. These program codes may be provided to a processor or controller of a general-purpose computer, a special purpose computer, or other programmable data processing devices, so that the program codes, when executed by the processor or controller, make the functions/functions specified in the flow diagrams and/or block diagrams Action is implemented. The program code may execute entirely on the machine, partly on the machine, as a stand-alone software package partly on the machine and partly on a remote machine or entirely on the remote machine or server.

在本公开的上下文中,机器可读介质可以是有形的介质,其可以包含或存储以供指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合地使用的程序。机器可读介质可以是机器可读信号介质或机器可读储存介质。机器可读介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体系统、装置或设备,或者上述内容的任何合适组合。机器可读存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦除可编程只读存储器(EPROM或快闪存储器)、光纤、便捷式紧凑盘只读存储器(CD-ROM)、光学储存设备、磁储存设备、或上述内容的任何合适组合。In the context of the present disclosure, a machine-readable medium may be a tangible medium that may contain or store a program for use by or in conjunction with an instruction execution system, apparatus, or device. A machine-readable medium may be a machine-readable signal medium or a machine-readable storage medium. A machine-readable medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination of the foregoing. More specific examples of machine-readable storage media would include one or more wire-based electrical connections, portable computer discs, hard drives, random access memory (RAM), read only memory (ROM), erasable programmable read only memory (EPROM or flash memory), optical fiber, compact disk read only memory (CD-ROM), optical storage, magnetic storage, or any suitable combination of the foregoing.

此外,虽然采用特定次序描绘了各操作,但是这应当理解为要求这样操作以所示出的特定次序或以顺序次序执行,或者要求所有图示的操作应被执行以取得期望的结果。在一定环境下,多任务和并行处理可能是有利的。同样地,虽然在上面论述中包含了若干具体实现细节,但是这些不应当被解释为对本公开的范围的限制。在单独的实施例的上下文中描述的某些特征还可以组合地实现在单个实现中。相反地,在单个实现的上下文中描述的各种特征也可以单独地或以任何合适的子组合的方式实现在多个实现中。In addition, while operations are depicted in a particular order, this should be understood to require that such operations be performed in the particular order shown, or in sequential order, or that all illustrated operations should be performed to achieve desirable results. Under certain circumstances, multitasking and parallel processing may be advantageous. Likewise, while the above discussion contains several specific implementation details, these should not be construed as limitations on the scope of the disclosure. Certain features that are described in the context of separate embodiments can also be implemented in combination in a single implementation. Conversely, various features that are described in the context of a single implementation can also be implemented in multiple implementations separately or in any suitable subcombination.

尽管已经采用特定于结构特征和/或方法逻辑动作的语言描述了本主题,但是应当理解所附权利要求书中所限定的主题未必局限于上面描述的特定特征或动作。相反,上面所描述的特定特征和动作仅仅是实现权利要求书的示例形式。Although the subject matter has been described in language specific to structural features and/or methodological acts, it is to be understood that the subject matter defined in the appended claims is not necessarily limited to the specific features or acts described above. Rather, the specific features and acts described above are merely example forms of implementing the claims.

Claims (5)

1.一种基于机器学习的冠心病患者出血风险评估装置,其特征在于,包括:1. A bleeding risk assessment device for patients with coronary heart disease based on machine learning, characterized in that it comprises: 获取模块,用于基于预设的采集时间,实时获取冠心病患者的PPG测量信号;The acquisition module is used to acquire the PPG measurement signal of the coronary heart disease patient in real time based on the preset acquisition time; 处理模块,用于通过巴特沃斯带通滤波器,去除所述PPG测量信号的基线漂移和高频噪音;A processing module, configured to remove baseline drift and high-frequency noise of the PPG measurement signal through a Butterworth band-pass filter; 基于所述PPG测量信号的波形,对所述PPG测量信号进行分段处理,得到M个信号分段;M为大于等于1的正整数;Based on the waveform of the PPG measurement signal, the PPG measurement signal is segmented to obtain M signal segments; M is a positive integer greater than or equal to 1; 评估模块,用于对所述M个信号分段中的每一个分段进行重采样,将集合S中每个分段重采样为相同长度,记为RS={RS1,RS2,…,RSM};其中,所述重采样长度为所述M个信号分段长度的中位数;An evaluation module, configured to resample each of the M signal segments, resampling each segment in the set S to the same length, denoted as RS={RS 1 ,RS 2 ,..., RS M }; wherein, the resampling length is the median of the M signal segment lengths; 将重采样后的分段信号分别代入如下公式,若计算结果大于预设阈值,则当前分段信号为质量合格信号:Substitute the resampled segmented signal into the following formula, if the calculation result is greater than the preset threshold, the current segmented signal is a qualified signal: Svalid={Si∣Si∈S,RSi∈RS,r(RSi,SRT)>0.9}S valid ={S i ∣S i ∈S,RS i ∈RS,r(RS i ,SR T )>0.9} 其中,所述SRT为RS的平均值;Wherein, the SRT is the average value of RS; 所述r为皮尔逊相关系数,取值范围为-1~1;The r is the Pearson correlation coefficient, the value range is -1~1; 评分模块,用于从质量合格的信号分段中,提取X个特征,输入到出血风险评估模型中,得到冠心病患者的出血风险评分;X为大于等于1的正整数。The scoring module is used to extract X features from the quality-qualified signal segments and input them into the bleeding risk assessment model to obtain the bleeding risk score of patients with coronary heart disease; X is a positive integer greater than or equal to 1. 2.根据权利要求1所述的装置,其特征在于,所述出血风险评估模型通过如下步骤训练得到:2. The device according to claim 1, wherein the bleeding risk assessment model is trained through the following steps: 生成训练样本集合,其中,训练样本包括带有标注信息的PPG信号对应的特征向量,其中标注信息为出血状况,出血标注为1,未出血标注为0;Generate a training sample set, wherein the training sample includes a feature vector corresponding to the PPG signal with labeling information, wherein the labeling information is bleeding status, bleeding is marked as 1, and non-bleeding is marked as 0; 利用所述训练样本集合中的样本对出血风险评估模型进行训练,以PPG信号对应的特征向量为输入,以出血状况为输出,当输出的出血状况与标注的出血状况的统一率满足预设阈值时,完成对出血风险评估模型的训练。Using the samples in the training sample set to train the bleeding risk assessment model, the feature vector corresponding to the PPG signal is used as input, and the bleeding condition is output, when the unity rate of the output bleeding condition and the marked bleeding condition satisfies a preset threshold , complete the training of the bleeding risk assessment model. 3.根据权利要求2所述的装置,其特征在于,所述PPG信号对应的特征向量包括:3. The device according to claim 2, wherein the eigenvector corresponding to the PPG signal comprises: 对采集的PPG信号进行分析,分别从时域、频域和小波包分解中提取与所述PPG信号对应的特征向量。The collected PPG signal is analyzed, and feature vectors corresponding to the PPG signal are extracted from time domain, frequency domain and wavelet packet decomposition respectively. 4.根据权利要求3所述的装置,其特征在于,所述对采集的PPG信号进行分析包括:4. The device according to claim 3, wherein said analyzing the collected PPG signal comprises: 对所述PPG信号、PPG信号的一阶导数和PPG信号的二阶导数进行分析。The PPG signal, the first derivative of the PPG signal and the second derivative of the PPG signal are analyzed. 5.根据权利要求4所述的装置,其特征在于,采用XGBoost算法对所述出血风险评估模型进行训练。5 . The device according to claim 4 , wherein the bleeding risk assessment model is trained using an XGBoost algorithm.
CN202111124787.9A 2021-09-25 2021-09-25 Coronary heart disease patient bleeding risk assessment method based on machine learning Active CN113925480B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111124787.9A CN113925480B (en) 2021-09-25 2021-09-25 Coronary heart disease patient bleeding risk assessment method based on machine learning

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111124787.9A CN113925480B (en) 2021-09-25 2021-09-25 Coronary heart disease patient bleeding risk assessment method based on machine learning

Publications (2)

Publication Number Publication Date
CN113925480A CN113925480A (en) 2022-01-14
CN113925480B true CN113925480B (en) 2023-03-21

Family

ID=79276889

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111124787.9A Active CN113925480B (en) 2021-09-25 2021-09-25 Coronary heart disease patient bleeding risk assessment method based on machine learning

Country Status (1)

Country Link
CN (1) CN113925480B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109493933A (en) * 2018-08-08 2019-03-19 浙江大学 A kind of prediction meanss of the adverse cardiac events based on attention mechanism
CN109620198A (en) * 2019-02-21 2019-04-16 天津惊帆科技有限公司 Cardiovascular index detection and model training method and device

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110319724A1 (en) * 2006-10-30 2011-12-29 Cox Paul G Methods and systems for non-invasive, internal hemorrhage detection
US10470719B2 (en) * 2016-02-01 2019-11-12 Verily Life Sciences Llc Machine learnt model to detect REM sleep periods using a spectral analysis of heart rate and motion
US11672437B2 (en) * 2016-04-25 2023-06-13 Happy Health, Inc. Method and device for tissue monitoring and heart rate detection
EP3468457A4 (en) * 2016-06-13 2020-02-19 Flashback Technologies, Inc. FAST DETECTION OF BLEEDING AFTER AN INJURY
EP3479763B1 (en) * 2017-11-06 2023-03-01 Tata Consultancy Services Limited System and method for photoplethysmogram (ppg) signal quality assessment
KR102588906B1 (en) * 2017-12-01 2023-10-13 삼성전자주식회사 Apparatus and method for assessment biological signal quality
CN112203582B (en) * 2018-06-01 2024-04-16 深圳市长桑技术有限公司 Pulse propagation time determination method and system
CN108742594B (en) * 2018-06-23 2020-11-24 桂林医学院附属医院 Wearable coronary heart disease detection device
US20200245967A1 (en) * 2019-02-04 2020-08-06 General Electric Company Localization of bleeding
US20220183606A1 (en) * 2019-03-26 2022-06-16 University Of Maryland, College Park Electrocardiogram waveform reconstruction from photoplethysmogram
CN110364261A (en) * 2019-07-17 2019-10-22 上海派兰数据科技有限公司 A method of establishing the acute coronary syndrome bleeding risk prediction model postoperative in interventional therapy
CN110739072A (en) * 2019-10-14 2020-01-31 中国人民解放军总医院 Methods and systems for assessing the occurrence of bleeding events
US20210212582A1 (en) * 2019-12-23 2021-07-15 Analytics For Life Inc. Method and system for signal quality assessment and rejection using heart cycle variability
CN112244803A (en) * 2020-10-30 2021-01-22 中国科学院合肥物质科学研究院 Sufficient risk check out test set of diabetes
CN112971795B (en) * 2021-02-07 2023-04-18 中国人民解放军总医院 Electrocardiosignal quality evaluation method
CN112971797B (en) * 2021-02-07 2024-12-13 北京海思瑞格科技有限公司 Method for quality assessment of continuous physiological signals

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109493933A (en) * 2018-08-08 2019-03-19 浙江大学 A kind of prediction meanss of the adverse cardiac events based on attention mechanism
CN109620198A (en) * 2019-02-21 2019-04-16 天津惊帆科技有限公司 Cardiovascular index detection and model training method and device

Also Published As

Publication number Publication date
CN113925480A (en) 2022-01-14

Similar Documents

Publication Publication Date Title
US20190026430A1 (en) Discovering novel features to use in machine learning techniques, such as machine learning techniques for diagnosing medical conditions
CN107529645A (en) A kind of heart sound intelligent diagnosis system and method based on deep learning
WO2021068781A1 (en) Fatigue state identification method, apparatus and device
EP4167854A1 (en) Generating physio-realistic avatars for training non-contact models to recover physiological characteristics
Chen et al. Edge2Analysis: a novel AIoT platform for atrial fibrillation recognition and detection
Xiang et al. Research of heart sound classification using two-dimensional features
Li et al. Multi-modal cardiac function signals classification algorithm based on improved DS evidence theory
CN115281685A (en) Sleep stage identification method, device and computer-readable storage medium based on abnormal detection
WO2022179645A2 (en) Electrocardiogram analysis method and apparatus, electronic device and storage medium
Sabeenian et al. Transfer learning-based electrocardiogram classification using wavelet scattered features
WO2020258507A1 (en) X-ray film classification method and apparatus, terminal, and storage medium
Lavanya et al. [Retracted] Wearable Sensor‐Based Edge Computing Framework for Cardiac Arrhythmia Detection and Acute Stroke Prediction
CN113925480B (en) Coronary heart disease patient bleeding risk assessment method based on machine learning
Erkuş et al. A new collective anomaly detection approach using pitch frequency and dissimilarity: Pitchy anomaly detection (PAD)
Mirabet-Herranz et al. Deep learning for remote heart rate estimation: A reproducible and optimal state-of-the-art framework
CN113693611A (en) Machine learning-based electrocardiogram data classification method and device
CN108338777A (en) A kind of pulse signal determination method and device
WO2021073161A1 (en) Elderly people registration method, apparatus and device based on voice recognition, and storage medium
Masullo et al. CaloriNet: From silhouettes to calorie estimation in private environments
Chen et al. Artificial intelligence for heart sound classification: A review
CN110020597A (en) It is a kind of for the auxiliary eye method for processing video frequency examined of dizziness/dizziness and system
CN115886833A (en) Electrocardiosignal classification method and device, computer readable medium and electronic equipment
García et al. Reviewing mobile apps to control heart rate in literature and virtual stores
Misbhauddin et al. An initial framework for mobile healthcare systems using deep neural networks
Omarov et al. Digital stethoscope for early detection of heart disease on phonocardiography data

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