CN116584952B - 心电图异常监测方法、装置、设备和存储介质 - Google Patents
心电图异常监测方法、装置、设备和存储介质 Download PDFInfo
- Publication number
- CN116584952B CN116584952B CN202310881402.6A CN202310881402A CN116584952B CN 116584952 B CN116584952 B CN 116584952B CN 202310881402 A CN202310881402 A CN 202310881402A CN 116584952 B CN116584952 B CN 116584952B
- Authority
- CN
- China
- Prior art keywords
- data
- convolution
- user
- electrocardiographic
- ecg
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 142
- 238000000034 method Methods 0.000 title claims abstract description 62
- 230000005856 abnormality Effects 0.000 title claims abstract description 34
- 238000012549 training Methods 0.000 claims abstract description 62
- 230000002159 abnormal effect Effects 0.000 claims abstract description 39
- 238000012806 monitoring device Methods 0.000 claims abstract description 20
- 238000013186 photoplethysmography Methods 0.000 claims description 141
- 238000011176 pooling Methods 0.000 claims description 50
- 238000005259 measurement Methods 0.000 claims description 39
- 230000006870 function Effects 0.000 claims description 33
- 238000012360 testing method Methods 0.000 claims description 28
- 238000007781 pre-processing Methods 0.000 claims description 27
- 230000004913 activation Effects 0.000 claims description 26
- 238000004590 computer program Methods 0.000 claims description 14
- 230000004217 heart function Effects 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 9
- 238000010606 normalization Methods 0.000 claims description 8
- 208000037170 Delayed Emergence from Anesthesia Diseases 0.000 claims description 5
- 238000011156 evaluation Methods 0.000 claims description 5
- 238000012795 verification Methods 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 2
- 230000000007 visual effect Effects 0.000 claims description 2
- 230000005986 heart dysfunction Effects 0.000 claims 4
- 230000000875 corresponding effect Effects 0.000 description 42
- 230000002547 anomalous effect Effects 0.000 description 18
- 239000011159 matrix material Substances 0.000 description 16
- 239000012634 fragment Substances 0.000 description 15
- 239000000523 sample Substances 0.000 description 11
- 238000005070 sampling Methods 0.000 description 11
- 230000000717 retained effect Effects 0.000 description 10
- 230000000747 cardiac effect Effects 0.000 description 8
- 238000001514 detection method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 238000012216 screening Methods 0.000 description 7
- 238000001914 filtration Methods 0.000 description 4
- 230000003993 interaction Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000010200 validation analysis Methods 0.000 description 3
- 206010003658 Atrial Fibrillation Diseases 0.000 description 2
- 206010003671 Atrioventricular Block Diseases 0.000 description 2
- 206010006578 Bundle-Branch Block Diseases 0.000 description 2
- 208000002330 Congenital Heart Defects Diseases 0.000 description 2
- 208000000418 Premature Cardiac Complexes Diseases 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000004064 dysfunction Effects 0.000 description 2
- 230000005831 heart abnormality Effects 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000003062 neural network model Methods 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 230000033764 rhythmic process Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 210000000707 wrist Anatomy 0.000 description 2
- HPTJABJPZMULFH-UHFFFAOYSA-N 12-[(Cyclohexylcarbamoyl)amino]dodecanoic acid Chemical compound OC(=O)CCCCCCCCCCCNC(=O)NC1CCCCC1 HPTJABJPZMULFH-UHFFFAOYSA-N 0.000 description 1
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 1
- 208000020446 Cardiac disease Diseases 0.000 description 1
- 206010019280 Heart failures Diseases 0.000 description 1
- 208000003443 Unconsciousness Diseases 0.000 description 1
- 206010047295 Ventricular hypertrophy Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 210000003423 ankle Anatomy 0.000 description 1
- 206010003119 arrhythmia Diseases 0.000 description 1
- 230000006793 arrhythmia Effects 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013145 classification model Methods 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000003205 diastolic effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 210000003414 extremity Anatomy 0.000 description 1
- 208000019622 heart disease Diseases 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 210000003141 lower extremity Anatomy 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 208000010125 myocardial infarction Diseases 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001314 paroxysmal effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000000718 qrs complex Methods 0.000 description 1
- 238000001303 quality assessment method Methods 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 230000004213 regulation of atrial cardiomyocyte membrane depolarization Effects 0.000 description 1
- 230000034225 regulation of ventricular cardiomyocyte membrane depolarization Effects 0.000 description 1
- 230000013577 regulation of ventricular cardiomyocyte membrane repolarization Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000036387 respiratory rate Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/746—Alarms related to a physiological condition, e.g. details of setting alarm thresholds or avoiding false alarms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Cardiology (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本申请涉及一种心电图异常监测方法、装置、计算机设备和存储介质。方法包括:获取用户的光电容积脉搏波数据;将光电容积脉搏波数据输入预先训练好的心电图异常监测模型中,得到心电图异常监测模型的输出值;其中,心电图异常监测模型是利用包括光电容积脉搏波数据以及标记光电容积脉搏波数据对应的心电图数据是否异常的标签的数据对的训练数据集训练得到的,其用于基于输入的光电容积脉搏波数据,输出指示所输入的光电容积脉搏波数据对应的心电图数据是否估计存在异常的输出值;当心电图异常监测模型的输出值指示用户的心电图数据估计存在异常时,使与用户相关联的用户设备发出警报。本方法能够方便地对心电图可能存在异常的情况进行持续监测。
Description
技术领域
本申请涉及心电图监测技术领域,特别是涉及一种心电图异常监测方法、装置、设备和存储介质。
背景技术
心电图(Electrocardiogram,ECG)是一种广泛用于监测/分类心脏异常(如心肌梗死、心室肥厚、心力衰竭等)的医学检查手段。
12导联ECG是在基层医疗中获得心脏功能信息的标准方法,其是通过在胸部和四肢的皮肤表面放置10个电极来测量,以记录心脏的电活动。ECG波形由代表心室去极化的QRS波群、代表心房去极化的P波以及代表心室复极化的T波构成。这样的波形是反映心脏功能的信息性和可靠的测量方法,因此广泛用于心脏疾病诊断的临床实践中。12导联ECG设备能够较为完善地测量ECG数据,但是其通常体积庞大,配有电极、中央单元以及显示器和键盘等附件,普通用户难以方便地使用12导联ECG设备对自身ECG进行日常持续监测。
相比之下,智能手表和健身追踪器等便携式ECG设备体积较小,但只能测量一条ECG导联,使得测量的ECG数据不够完善。此外,它们还需要一些用户发起的动作,例如,握持探测器以闭合导电电路,以便进行测量,因此用户需要选择合适的时机来主动发起对ECG的测量。而心律失常等心脏异常具有阵发性性质,对ECG测量时机的选择即使对于临床专家而言也是一项挑战。因此,用便携式ECG设备也难以持续地采集和监测ECG。
因此,上述现有的ECG测量方法仍然存在着改进的空间。
发明内容
基于此,有必要针对上述技术问题,提供一种能够方便地对心电图可能存在异常的情况进行持续监测的心电图(ECG)异常监测方法、装置、设备和存储介质。
在一方面,提供一种ECG异常监测方法,所述方法包括:
获取用户的光电容积脉搏波(Photoplethysmography,PPG)数据;
将所述PPG数据输入预先训练好的ECG异常监测模型中,以得到所述ECG异常监测模型的输出值;其中,所述ECG异常监测模型是利用包括PPG数据以及标记所述PPG数据对应的ECG数据是否异常的标签的数据对的训练数据集训练得到的,所述ECG异常监测模型用于基于输入的PPG数据,输出指示所输入的所述PPG数据对应的ECG数据是否估计存在异常的输出值;
当所述ECG异常监测模型的输出值指示所述ECG数据估计存在异常时,使与所述用户相关联的用户设备发出警报。
在另一方面,提供一种ECG异常监测装置,包括:
数据获取模块,用于获取用户的PPG数据;
模型监测模块,用于将所述PPG数据输入预先训练好的ECG异常监测模型中,以得到所述ECG异常监测模型的输出值;其中,所述ECG异常监测模型是利用包括PPG数据以及标记所述PPG数据对应的ECG数据是否异常的标签的数据对的训练数据集训练得到的,所述ECG异常监测模型用于基于输入的PPG数据,输出指示所输入的所述PPG数据对应的ECG数据是否估计存在异常的输出值;
异常警报模块,用于当所述ECG异常监测模型的输出值指示所述ECG数据估计存在异常时,使与所述用户相关联的用户设备发出警报。
在另一方面,提供一种ECG异常监测设备,包括PPG数据探测器、输出装置、存储器和处理器;
所述PPG数据探测器用于监测用户的PPG数据,并将测得的所述PPG数据传输至所述处理器;
所述存储器存储有计算机程序;并且
所述处理器执行所述计算机程序时,用于从所述PPG数据探测器接收所述PPG数据,并实现如上所述的方法,以在当所述ECG异常监测模型的输出值指示所述用户的ECG数据估计存在异常时,利用所述输出装置使与所述用户相关联的用户设备发出警报。
在另一方面,提供一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如上所述的方法。
上述ECG异常监测方法、装置、设备和存储介质,由于PPG数据能够以被动的方式通过例如智能手表等便携式可穿戴设备对用户进行测量得到,因此通过获取用户的PPG数据,能够方便而持续地对用户进行监测,通过巧妙地使用PPG数据对应的ECG数据是否异常的标签来对PPG数据进行标记形成的数据对进行训练,能够得到可以基于输入的PPG数据,输出指示所输入的所述PPG数据对应的ECG数据是否估计存在异常的输出值的ECG异常监测模型,从而利用该模型可以基于测得PPG数据监测到对应的ECG数据是否很可能发生异常时,并在监测到发生异常时,使与所述用户相关联的用户设备发出警报,以建议用户主动进行ECG测量获知ECG是否异常的准确结果。
附图说明
图1为一个实施例中ECG异常监测方法的应用环境图;
图2为一个实施例中ECG异常监测方法的流程示意图;
图3为一个实施例中通过预处理而被保留的PPG数据片段的示例;
图4为一个实施例中未通过预处理而被筛除的数据片段的示例;
图5为一个实施例中ECG异常监测模型的架构示意图;
图6为一个实施例中残差-双卷积注意力块的架构示意图;
图7为一个实施例中双卷积注意力块的架构示意图;
图8为一个实施例中卷积注意力块的架构示意图;
图9为一个实施例中导联II、V和AVR的ECG原始数据片段以及对应的PPG原始数据片段的示例;
图10为一个实施例中被标记为正常的ECG预处理数据片段及其对应的PPG预处理数据片段的示例;
图11为一个实施例中被标记为异常的ECG预处理数据片段及其对应的PPG预处理数据片段的示例;
图12为一个实施例中ECG异常监测装置的结构示意图;
图13为一个实施例中ECG异常监测设备的结构示意图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
由于心脏功能与ECG之间的密切关联,且ECG数据包含了较为丰富的信息内容,尤其是ECG的波形结构和蕴涵的时间信息,因此12导联ECG是临床中用来监测心脏功能异常的标准手段。然而,现有的ECG测量设备均存在着使用不便的缺陷,例如12导联ECG设备庞大,不便于用户日常监测,而便携式ECG设备通常仅能测量一条ECG导联,且通常需要由用户通过设备“主动地”采样ECG,难以把握合适的测量时机。
PPG是一种通过检测血容量变化来测量心脏周期的光学技术。PPG可以通过例如脉搏血氧计等来“被动地”测量,并且通常作为一种功能嵌入到可穿戴设备中,如健身腕带和智能手表。PPG的波形由一个收缩波和一个舒张波组成,比ECG的波形更简单、更平滑。PPG可方便地用于长期持续监测例如脉搏和呼吸率等生理参数。然而,根据PPG对心脏异常进行分类或监测非常困难,因为它是一种间接监测心脏的运作的外周测量信号,且与ECG相比,其相对平滑的形态反映的心脏信息较少。因此,现有的PPG多用于测量脉搏频率、呼吸频率、血压等生理参数。
本申请提出了一种ECG异常监测方法,其使用连续监测的、普遍存在的PPG信号,来监测ECG可能存在异常的情况,所述异常包括但不限于房室传导阻滞、窦性心律异常、束支传导阻滞、房颤、早搏等常见的心脏功能异常。以提示用户进行ECG检查以得知心脏风险。
本申请提供的ECG异常监测方法,可以应用于如图1所示的应用环境中。其中,终端102通过网络与服务器104进行通信。终端102可以是具有测量PPG数据的终端设备,例如终端102可以设置有PPG数据探测器,例如脉搏血氧计。示例地,终端102可以是例如内置有PPG数据探测器的智能手表、智能手环、智能耳机等便携式可穿戴终端设备。终端102也可以具有其他形态,例如配置有脉搏血氧计的医用监仪,又例如内置有脉搏血氧计的智能手环和与该智能手环无线通信的智能手机组成的系统等等,本领域技术人员将能想到许多其他变形,本申请对此不作限定。服务器104可以用独立的服务器或者是多个服务器组成的服务器集群来实现。
本申请的ECG异常监测方法可以包括训练阶段和测试阶段。在训练阶段,可以利用训练数据集对ECG异常监测模型进行训练,以得到训练好的ECG异常监测模型,该异常包括但不限于房室传导阻滞、窦性心律异常、束支传导阻滞、房颤、早搏等常见的心电导联异常。该训练过程可以在服务器104上执行,并在完成训练得到训练好的ECG异常监测模型后,将训练好的ECG异常监测模型从服务器104部署到终端102上以备使用。在完成训练得到训练好的ECG异常监测模型后,在测试阶段,终端102可以利用其具有的PPG探测器监测用户的PPG数据,并执行本申请实施例的ECG异常监测方法,以当ECG异常监测模型的输出值指示用户的ECG数据估计存在异常时,使与该用户相关联的用户设备发出警报,例如使终端102自身发出警报和/或使终端102之外的其他关联设备发出警报。可以理解,上述ECG异常监测方法的部分和全部步骤在何处执行,可以在终端102和服务器104之间进行适当调整,例如,训练好的ECG异常监测模型也可以部署在服务器104上,终端102将采集的PPG数据发送至服务器104,由服务器104利用模型进行监测,并将监测结果发回给终端102,等等,本申请对此不作限定。
进一步地,上述警报可以用于指示对用户进一步进行ECG测量。该ECG测量可以由终端102进行,例如,终端102上可以部署有ECG测量装置,例如ECG电极或ECG传感器,从而利用终端102对用户进行ECG测量。该ECG测量也可以由终端102之外的其他任何ECG测量设备测量。
在一个实施例中,如图2所示,提供了一种ECG异常监测方法,在测试阶段,该方法包括以下步骤S210-S230:
S210,获取用户的PPG数据。
PPG描记法是借用光电手段在活体组织中检测血液容积变化的一种无创检测方法。PPG可以通过PPG数据探测器进行检测,PPG数据探测器是能够以无创方式检测用户的PPG的任何探测器设备,其可以包括光源和光电检测器,光源向被测组织发射预定波长的光,光电检测器接收被测组织反射或透射的光,基于被测组织反射或透射的光的强度得到检测的PPG信号。PPG数据探测器例如可以是脉搏血氧计。
在一个示例中,脉搏血氧计例如可以内置于终端102中,该终端102例如可以是智能手表,当用户佩戴该智能手表时,位于智能手表表盘底部(靠近用户皮肤侧)的脉搏血氧计检测用户手腕内侧的PPG数据。在其他示例中,终端102也可以具有其他形态,PPG数据的测试部位也不限于用户手腕内侧,而可以是例如下肢踝内侧、胸部等其他部位。可以理解,在测试阶段采集的测试数据的测试部位与在训练阶段所使用的训练数据的测试部位相同。
PPG数据是表征用户的脉搏波幅值随着时间变化的数据,其可以表示为一维数据,即,将脉搏波的采样值按采样时间顺序排列得到的PPG一维数据序列。
在一个实施例中,获取用户的PPG数据是由便携终端设备实时监测的用户的PPG数据。
在一个示例中,终端102以预定采样频率对被测用户的测试部位进行持续采样监测。终端102从其持续监测的PPG信号中,截取预定时长的PPG数据片段,作为待输入模型的PPG数据。该预定采样频率和预定时长可以根据实际需要设定,例如预定采样频率可以是125Hz,预定时长可以是10s,则一个10s长度的PPG数据片段大约可以包括1250个采样值。图3中的深灰色曲线表示的数据示出了一个10s长度的PPG原始数据片段的示例,其中,横坐标表示时间,其用采样点序数表示,纵坐标表示幅值,本申请其余附图中显示的ECG和PPG的数据片段中的横纵坐标也具有同理的含义。
采集得到的原始PPG数据可能存在噪音或质量不高等问题,而影响监测结果的准确性。因此,可以对采集得到的原始数据进行预处理,再使用预处理后的数据输入模型进行监测。
在一个实施例中,上述步骤S210可以包括:获取用户的PPG原始数据,对PPG原始数据进行预处理,以得到PPG预处理数据。从而利用PPG预处理数据作为要输入到模型中的PPG数据。相应地,在后续步骤S220中,可以将PPG预处理数据输入预先训练好的ECG异常监测模型中。
在一个实施例中,对PPG原始数据进行预处理,以得到PPG预处理数据,可以包括如下步骤中的一者或多者:
1)平坦度检测和信号筛选:
在此步骤中,检测信号的平坦度是否达到预定标准,若是,则保留该信号,否则,丢弃该信号。
示例地,对于截取的10s时长的PPG原始数据片段,如果其中存在任意60个连续采样点(约0.5秒)的采样值未超过阈值(1e-5),则将其视为平坦信号,并丢弃这整个10s的数据片段。否则,保留该数据片段。
2)归一化:
在此步骤中,将信号归一化为零平均值和单位方差。
示例地,可以将前述平坦度检测和信号筛选步骤保留的数据片段中的各个采样值归一化为零平均值和单位方差,得到归一化后的数据片段。
3)滤波:
在此步骤中,对信号进行滤波处理。
示例地,可以采用三阶带通巴特沃斯滤波器,对前述步骤归一化后的数据片段进行滤波,以得到滤波后的数据片段。其中,该滤波器的低频带截止为0.5Hz,高频带截止为8Hz。
4)峰值检测和信号筛选:
在此步骤中,检测信号的单位长度内有效峰值数量是否达到预定标准,若是,则保留该信号,否则,丢弃该信号。
示例地,可以使用Python工具箱“HeartpPy”来检测前述步骤滤波后的数据片段内的有效峰值的数量。如果该10s数据片段内的有效峰值的数量低于5个(对应于低于30bpm),则丢弃该数据片段。否则,保留该数据片段。
5)偏度SQI检测和信号筛选:
在此步骤中,检测信号的偏度是否达到预定标准,若是,则保留该信号,否则,丢弃该信号。
示例地,偏度以滑动窗口的方式来计算,窗口宽度为250个采样值(两秒),步幅为125个采样值(一秒)。对于前述峰值检测和信号筛选步骤中保留的10s的数据片段,即可计算0~2s、1~3s、2~4s、3~5s、4~6s、5~7s、6~8s、7~9s、8~10s共9个时间窗口中每个时间窗口内采样值的偏度。如果计算得到的9个偏度中,大多数偏度(超过50%)为负,则认为该数据片段是低质量信号,丢弃该数据片段。否则,保留该数据片段。
6)离群采样值替换:
在此步骤中,检测信号中的离群采样值,将检测到的离群采样值替换为采样值中位数,以得到消除离群采样值后的信号。
示例地,可以应用Hampel滤波器来检测离群点。在前述偏度SQI检测和信号筛选步骤中保留的数据片段中,对于每连续10个采样值,计算中位绝对偏差(MAD)。然后根据假设正态分布的MAD值来估计该数据片段的标准差(std)。如果采样值与该数据片段的中位数相差3 std,则该采样值被检测为离群点。将检测到的离群点替换为该数据片段的中位数。如此,能够将数据片段开头和结尾的较多离群点替换掉。
步骤1)可视为质量筛查步骤。步骤2)和3)是清洁和去噪,步骤4)至6)是质量控制步骤。参见图3所示,图3中示出了一个通过上述预处理步骤1)至6)而被保留的PPG数据片段,其中深灰色曲线代表原始数据片段,浅灰色曲线代表预处理后得到的预处理数据片段。参见图4所示,图4中示出了由于未能通过上述预处理步骤中的不同质量控制标准而被筛除的数据片段的示例。其中,图4中的(a)表示因平坦度未达标而被筛除的数据片段,(b)表示因有效峰值数量未达标而被筛除的数据片段,(c)表示因偏度未达标而被筛除的数据片段。
在本示例中,对于任意原始数据片段,经过上述预处理步骤后得到的预处理数据片段,可用于后续步骤输入模型中进行监测;而如果任意原始数据片段在上述任一步骤中被丢弃,则不再使用该原始数据片段进行监测。即,仅使用通过预处理步骤删选后的满足要求的预处理数据片段进行监测,从而保障监测的准确性。可以理解,在测试阶段对原始数据执行的预处理操作可以与在训练阶段对原始数据执行的预处理操作相同,以保障模型的适用性。
S220,将PPG数据输入预先训练好的ECG异常监测模型中,以得到ECG异常监测模型的输出值。
其中,本申请的ECG异常监测模型是利用包括PPG数据以及标记PPG数据对应的ECG数据是否异常的标签的数据对的训练数据集训练得到的,该ECG异常监测模型的架构和训练过程将在后续描述。
训练好的ECG异常监测模型可以用于基于输入的PPG数据,输出指示所输入的PPG数据对应的ECG数据是否估计存在异常的输出值。
示例地,当将前述步骤中取得的10s的数据片段输入该ECG异常监测模型,该ECG异
常监测模型可以输出监测类别概率作为输出值,将该监测类别概率与设定的阈值T进行
比较,即可确定对应的ECG数据是否估计存在异常。例如,当≥T时指示ECG数据估计存在
异常,而当<T时指示ECG数据估计不存在异常。
S230,当ECG异常监测模型的输出值指示ECG数据估计存在异常时,使与用户相关联的用户设备发出警报。
在一个实施例中,S230中的与用户相关联的用户设备,可以包括上述的终端102,例如是用户佩戴的便携终端设备,以告知用户自身存在ECG异常风险,应及时进行ECG测量。
在另一个实施例中,S230中的与用户相关联的用户设备还可以是与用户关联的监护者的监护终端设备,例如用户的家属、看护者所用的终端设备,或者是用户的家庭医生、主治医生所用的终端设备,或者是用户所属的监护机构的终端设备或服务器等等,以告知监护者该用户存在ECG异常风险,应及时进行ECG测量。
在一个实施例中,使与用户相关联的用户设备发出警报包括以下各项中的一项或多项组合:在用户设备的显示屏上显示指示用户的ECG数据估计存在异常的视觉信息、使用户设备的指示灯发出指示用户的ECG数据估计存在异常的亮灯提示、使用户设备的振动器做出指示用户的ECG数据估计存在异常的振动提示、使用户设备的声音输出装置发出指示用户的ECG数据估计存在异常的声音提示。
示例地,使与用户相关联的用户设备发出警报的内容除了包括告知ECG数据估计存在异常的监测结果外,还包括建议用户主动进行ECG测量的建议信息。
进一步地,在一个实施例中,上述方法还可以包括:当所述ECG异常监测模型的输出值指示所述ECG数据估计存在异常时,向预定警报接收设备拨打警报电话。
其中,该警报电话可以由警报发起设备发起,该警报发起设备可以是任何合适的终端或服务器。例如警报电话可以由上述用户所佩戴的便携终端设备发起,在此情况下,该便携终端设备需要是具有拨打电话功能的终端设备,例如是智能手表或者是智能手环和与智能手环配对的手机;或者也可以由例如上述服务器104或其他服务器或终端设备来发起。
该警报电话可以是由终端或服务器自动发起的或是基于用户的指令来决定是否发起,例如终端或服务器可以经由用户佩戴的便携终端设备向用户提供是否发起警报电话的选项,然后根据用户对选项的选择输入而决定是否向预定警报接收设备拨打警报电话,又或者终端或服务器可以在向用户提供上述选项后,如果预定时段内未收到用户的选择输入,则认为用户可能处于失去意识等危险状况中,则自动向预定警报接收设备拨打警报电话。
用于接收警报电话的预定警报接收设备可以根据需要预先设定,例如该预定警报接收设备可以包括上述与用户关联的监护者的监护终端设备,和/或指定医疗机构的警报终端(例如可通过拨号112而拨打的医院急救电话终端)等等。
取决于预先设定或根据用户的选择输入,警报电话可以是接通实时通话以使用户与预定警报接收设备处的接收人通话,警报电话也可以是包括预先设定的语音播报内容,例如可以包括告知用户的ECG数据估计存在异常的监测结果有关的语音播报内容,此外还可以包括例如用户的姓名、性别、年龄、居所等个人信息以及用户的实时位置等信息的语音播报内容,以方便警报电话接收者快速知晓用户的当前状况和实时位置。
在一个实施例中,上述方法还可以包括:当ECG异常监测模型的输出值指示ECG数据估计存在异常时,对用户进行ECG测量,以得到指示用户的ECG是否异常的测量结果。
其中,该ECG测量可由通过设置有ECG测量装置的终端102进行,也可以通过与终端102不同的其他ECG测量设备,例如设置于医疗机构中的12导联ECG,或者家用ECG测量设备等进行。例如,在步骤S230中终端102发出警报提醒用户进行ECG测量后,终端102可以基于用户输入来确定是否触发由终端102对用户的ECG测量,又例如,终端102可以在输出值指示ECG数据估计存在异常时,自动触发对用户的ECG测量,又例如,在步骤S230中终端102发出警报提醒用户/监护者进行ECG测量后,用户/监护者可以自行选择利用终端102之外的其他ECG测量设备对用户进行ECG测量,例如用户可以自行前往医疗机构利用12导联ECG进行更专业全面地ECG测量,并得到ECG是否异常的测量结果。其中,当由终端102执行上述ECG测量时,终端102可以由自身或经由服务器等外部设备对测得的ECG数据进行分析,并将测得的ECG数据和分析得到的ECG是否异常的测量结果通过终端102推送给用户。
上述ECG异常监测方法,由于PPG数据能够以被动的方式通过例如智能手表等便携式可穿戴设备对用户进行测量得到,因此通过获取用户的PPG数据,能够方便而持续地对用户进行监测,通过巧妙地使用PPG数据对应的ECG数据是否异常的标签来对PPG数据进行标记形成的数据对进行训练,能够得到可以基于输入的PPG数据,输出指示所输入的所述PPG数据对应的ECG数据是否估计存在异常的输出值的ECG异常监测模型,从而利用该模型可以基于测得PPG数据监测到对应的ECG数据是否很可能发生异常时,并在监测到发生异常时,使与所述用户相关联的用户设备发出警报,以建议用户主动进行ECG测量获知ECG是否异常的准确结果。
以上描述了在测试阶段,终端102或服务器104可能执行的方法步骤,其中在测试阶段的步骤S220中,使用了预先训练好的ECG异常监测模型,这要求在测试阶段之前的训练阶段中,预先构造ECG异常监测模型的模型架构,并对该构造好的ECG异常监测模型进行训练,优化确定ECG异常监测模型的模型参数,以得到训练好的信号质量评估模型。
本申请提出了一种双卷积注意力网络(Dual-Convolutional AttentionNetwork,DCA-Net)作为上述ECG异常监测模型。该双卷积注意力网络的特征是,包括有一个或多个残差-双卷积注意力(Res-DAC)块。每个残差-双卷积注意力块是在残差块中加入双卷积注意力(dual-convolutional-attention,DCA)块而形成。其中,该双卷积注意力块包括用于在通道维度上对数据进行卷积的第一卷积注意力块以及用于在时间维度上对数据进行卷积的第二卷积注意力块。
在一个实施例中,参见图5所示,该ECG异常监测模型至少可以包括串行设置的多个残差-双卷积注意力块1~N。
在一个实施例中,参见图5所示,该ECG异常监测模型还可以包括位于多个残差-双卷积注意力块之前,且与多个残差-双卷积注意力块串行设置的输入卷积层和池化层,以及位于多个残差-双卷积注意力块之后,且与多个残差-双卷积注意力块串行设置的平均池化层和全连接层。
其中,输入卷积层接收对ECG异常监测模型输入的PPG数据,输入卷积层的输出馈入池化层,池化层的输出馈入多个残差-双卷积注意力块中的首个残差-双卷积注意力块,多个残差-双卷积注意力块中的前一者的输出馈入后一者,多个残差-双卷积注意力块中的最后一个残差-双卷积注意力块的输出馈入平均池化层,平均池化层的输出馈入全连接层,全连接层输出ECG异常监测模型的输出值。
其中,全连接层的最后一层可以具有单个神经元,该单个神经元输出一个1*1维的
输出值,即概率。将该概率与设定的阈值T进行比较,即可得到ECG是否估计存在异常的
监测结果。
示例地,模型中的各个卷积层,例如输入卷积层,可以为一维卷积层,以适应对PPG一维数据序列的处理需求。
可以理解,图5中的模型架构仅为示例,在模型中包括有一个或多个如本申请提出的残差-双卷积注意力块的基本前提下,本领域技术人员还可以对模型的架构做出许多变形,例如残差-双卷积注意力块的数量和串行/并行方式可以改变,模型中的卷积层、池化层数量和位置等也可以改变。
在一个实施例中,参见图6和图7所示,该多个残差-双卷积注意力块中的每个残差-双卷积注意力块可以包括第一主路径以及第一短路(shortcut)分支路径。该第一主路径上串行设置有一个或多个卷积层以及双卷积注意力块。其中,双卷积注意力块包括串行或并行设置的第一卷积注意力块和第二卷积注意力块,第一卷积注意力块用于在通道维度上进行卷积,第二卷积注意力块用于在时间维度上进行卷积。该第一短路分支路径并联于第一主路径上串行设置的一个或多个卷积层以及双卷积注意力块的两端。
其中,一个或多个卷积层以及双卷积注意力块中的最前者以及第一短路分支路径接收先前路径的输出,一个或多个卷积层以及双卷积注意力块中的后一者接收前一者的输出,一个或多个卷积层以及双卷积注意力块中的最后者的输出与第一短路分支路径的输出矩阵相加后馈入后续路径。
在一个实施例中,参见图6所示,每个残差-双卷积注意力块中的各个卷积层的输出侧还分别设置有批量归一化层以及激活函数层。其中,在每个残差-双卷积注意力块中,除了最后一个卷积层之外的每个卷积层的输出侧的批量归一化层以及激活函数层设置于下一个卷积层之前,而最后一个卷积层的输出侧的批量归一化层设置于该卷积层与双卷积注意力块之间,双卷积注意力块设置于与第一短路分支路径的输出的相加运算之前,最后一个卷积层的输出侧设置的批量归一化层设置于与第一短路分支路径的输出的相加运算之后。
在一个具体示例中,参见图6所示,每个残差-双卷积注意力块中设置有两个卷积层,且它们为一维卷积层。激活函数层为ReLU激活函数层。在其他示例中,残差-双卷积注意力块中的卷积层的数量和维度可以根据情况而更改,激活函数层也可以为其他类型的激活函数层。
在一个实施例中,参见图7所示,第一卷积注意力块和第二卷积注意力块串行设置,且第一卷积注意力块位于第二卷积注意力块之前。将通道方面的第一卷积注意力块置于时间方面的第二卷积注意力块之前,使得双卷积注意力块先学习通道间交互,然后基于通道参与输出来学习时间方面交互,这比两者并行或通道方面的第一卷积注意力块在后的顺序表现更好。在其他示例中,也可以是第二卷积注意力块位于第一卷积注意力块之前或者第一卷积注意力块和第二卷积注意力块并行设置。
在一个实施例中,参见图8所示,第一卷积注意力块和第二卷积注意力块中的每一者包括:第二主路径以及第二短路分支路径。第二主路径上设置有最大池化层和平均池化层、卷积层以及激活函数层,最大池化层和平均池化层并行设置,并行设置的最大池化层和平均池化层依次与卷积层和激活函数层串行设置。第二短路分支路径并联于第二主路径上设置的最大池化层和平均池化层、卷积层以及激活函数层的两端。
其中,第二短路分支路径、最大池化层和平均池化层分别接收先前路径的输出,最大池化层的输出与平均池化层的输出联合(concatenation)所得的输出馈入卷积层,卷积层的输出馈入激活函数层,激活函数层的输出与第二短路分支路径的输出相乘后馈入后续路径。
在一个具体示例中,参见图8所示,第一/第二卷积注意力块中具有一个卷积层,且为一维卷积层。激活函数层为Sigmoid激活函数层。在其他示例中,第一/第二卷积注意力块中的卷积层的数量和维度可以根据情况而更改,激活函数层也可以为其他类型的激活函数层。
其中,卷积注意力块依赖于将卷积运算应用到数据的某个维度,以学习沿着该维度的数据点之间的交互。与MLP注意力模块相比,卷积注意力在计算上更高效,因为卷积层中不需要参数。从图8中可见,第一卷积注意力块和第二卷积注意力块具有类似的架构,但是他们应用于不同的方向,第一卷积注意力块应用于卷积通道方向,第二卷积注意力块应用于时域方向。
在上述实施例中,各个层/块的参数,例如激活函数(activation function)、卷积核尺寸(kernel size)、步长、补丁(padding)等的具体配置可以根据实际需要设定。
具体来说,对于图7和图8所示的具体示例,以图8中的卷积层为卷积核大小为7的
一维卷积层为例,设通道方面的第一卷积注意力块的输入(也是DCA模块的输入)为,其中N是采样点/批量大小,C是从原始残差块输出的卷积通道的数量,D是信
号长度(记为时间维度)。第一卷积注意力块的输出为:
(1)
其中, (2)
其中, (3)
并且, (4)
等式(1)中,表示Hadamard乘积(元素方面矩阵乘法);表示Sigmoid激活函
数,表示卷积核大小为7且在通道方向上卷积的1D卷积层。等式(2)中的表示
沿时间方向的矩阵联合,等式(3)和(4)中的下标D表示池化层也在时间维度上应用。函数输出通道方面注意力权重,这些权重进一步扩展到相同形状的X以缩放X。然
后将缩放后的X(表示为)用作向时间方面注意力的输入。
同样,设时间方面的第二卷积注意力块的输出为。
(5)
其中, (6)
其中, (7)
并且, (8)
等式(7)-(8)表示首先在通道维度上操作最大池化和平均池化,因此,和 。然后,将池化后的矩阵沿着通道维度联合,(等式(6))。
一维卷积层在时间维度上卷积,并输出单个通道,使得。最后,这
个卷积后的矩阵由Sigmoid激活函数激活并用于缩放输入的时域。
在一个具体示例中,可以通过将已知的ResNet-34模型中的二维卷积层均替换为一维卷积层,将ResNet-34模型中最后的全连接层替换为输入为512维输出为1维的全连接层,并在ResNet-34模型的每个残差块中增加残差之前增加一个双卷积注意力块,而得到一个具体示例的DCA-Net模型。可以理解,该具体示例的DCA-Net模型将包括依次串行的输入卷积层、池化层、16个残差-双卷积注意力块、平均池化层以及全连接层。该具体示例的DCA-Net模型的架构可以参考图5所示,其中N取值N=16,并且其中每个残差-双卷积注意力块的架构可以参考图6至图8所示。其中,示例地,在每个残差-双卷积注意力块中,第一卷积注意力块和第二卷积注意力块的卷积核大小均可以为7,步长均为1,平均池化层可以为一维自适应平均池化层(AdaptiveAvgPool1d),模型中其他层和块的参数均示例地可以参照已知ResNet-34模型的参数设置。
示例地,以通过批次训练(batch training)来执行训练为例,来说明该具体示例
的ResNet-34模型的前向传播过程。当批次大小为64,向该具体示例的DCA-Net模型输入一
个64x1x1250的一维PPG数据时,经过输入卷积层(步长为2)后输出64×64×625的数据,经
过(最大)池化层(步长为2)后输出64×64×313的数据,然后该数据输入到第一个残差-双
卷积注意力块中,在依次经过图6中所示的卷积层、批量归一化层、ReLU激活函数层、卷积
层、批量归一化层后,输入到双卷积注意力层中,输入数据首先进入通道方面的第一卷积注
意力块,在第一卷积注意力块内,参见图7,平均池化层对输入数据在时间维度进行平均池
化而输出一个64×64×1矩阵,同时最大池化层对输入数据在时间维度进行最大池化而输
出另一个64×64×1矩阵,这两个矩阵沿时间维度联合(类似于拼接操作)在一起得到64×
64×2的矩阵。该64×64×2的矩阵经过卷积层而被沿通道维度进行卷积,以学习64×64×1
的注意力权重,再沿时间维度被复制展开成64×64×313的注意力分数。将此注意力分数与
经过第二短路分支路径输出的原输入数据进行Hadamard乘积而输出64×64×313矩阵。接
着,第一卷积注意力块输出的64×64×313矩阵输入时间方面的第二卷积注意力块,在第二
卷积注意力块内,参见图7,平均池化层对64×64×313矩阵在通道维度进行平均池化而输
出一个64×1×313的矩阵,同时最大池化层对输入数据在通道维度进行最大池化而输出另
一个64×1×313的矩阵。这两个矩阵联合在一起得到64×2×313的矩阵。该64×2×313的
矩阵经过卷积层而被沿时间维度进行卷积,得到64×1×313的注意力权重,再沿通道维度
被复制展开成64×64×313的注意力分数。将此注意力分数与经过第二短路分支路径输出
的原输入64×64×313矩阵进行Hadamard乘积,而输出64×64×313的矩阵。该64×64×313
的矩阵作为第二卷积注意力块的输出接着被输入到下一个残差-双卷积注意力块中,在经
过16个残差-双卷积注意力块后,输出数据然后经过平均池化层和全连接层后,输出模型的
输出值。
在训练阶段中,通过对上述任一实施例中构造的ECG异常监测模型进行训练,能够确定训练好的ECG异常监测模型。具有上述构造的ECG异常监测模型可以利用各种常规的或新型的训练方法来训练。
其中,为了执行模型训练,需要构建训练数据集,在本申请中,训练数据集可以包括多组PPG数据-标记对应的ECG数据是否异常的标签的数据对。
训练数据集的构建过程可以包括如下步骤:
1)获取多个导联的ECG信号及其对应的PPG信号。
其中,ECG信号与PPG信号对应是指,该ECG信号与该PPG信号是在同一受试者身上同时测量的,以确保该ECG信号与该PPG信号之间时间对准并彼此关联。
示例地,可以从MIMIC-III波形数据库(Waveform Database,WDB)匹配子集中获取多个导联的ECG信号及其对应的PPG信号。MIMIC-III WDB与临床记录相联系——匹配子集中的所有患者均被识别并被匹配到医疗记录系统。MIMIC-III WDB包括多个导联的ECG信号及其配对的PPG信号。对于每例患者,在其住院期间的不同时间点采集不同长度的波形数据。这些信号的采样频率为125 Hz。由于不同的ECG记录具有不同的导联数量和类型,在本示例中,选择了3种最常见和最有代表性的ECG导联的数据,即导联II、V和AVR。因此,将使用这三根导联的ECG信号以及对应的PPG信号记录纳入到本分析中。然而在其他示例中,也可能选择更多或更少的ECG导联,所选导联的类型也可能变更。此外,为了降低噪声水平,只选取时长超过预定时长,例如两分钟的信号的记录。时长不到两分钟的信号将被筛除。
2)从获取的多个导联的ECG信号及其对应的PPG信号中截取预定时长的多个导联的ECG原始数据片段及其对应的PPG原始数据片段。
训练过程中训练数据要截取的数据片段的预定时长与测试过程中截取的数据片段的预定时长相同,例如也为10秒。如此,可以从每个ECG/PPG信号中截取出第二个10秒的数据片段,以避免初始时段采样值的不稳定。如此,截取得到的10秒时长的导联II、V和AVR的ECG原始数据片段以及对应的10秒时长的PPG原始数据片段的示例如图9中所示。
3)对所截取的预定时长的多个导联的ECG原始数据片段及其对应的PPG原始数据片段分别进行预处理,以预处理并筛选出满足要求的预定时长的多个导联的ECG预处理数据片段及其对应的PPG预处理数据片段。
对于10s的PPG原始数据片段,可以执行如上所述测试过程中相同的包括步骤1)至6)的预处理过程,以筛选出满足要求的PPG预处理数据片段。
对于10s的导联II、V和AVR的ECG原始数据片段,可以使用0.5Hz的5阶高通巴特沃兹滤波器,随后使用电力线频率为50Hz的电力线滤波,以得到10s的导联II、V和AVR的ECG预处理数据片段。然后,使用预定信号质量评估方法来评估每个ECG预处理数据片段的信号质量。该评估输出三类质量中的一类——不可接受、几乎不可接受或优秀。当某一导联的ECG预处理数据片段被评估为优秀时,该导联的ECG预处理数据片段被保留,反之,当某一导联的ECG预处理数据片段被评估为不可接受、几乎不可接受时,该导联的ECG预处理数据片段被丢弃。
如此,对于任意一组10s的多个导联的ECG原始数据片段及其对应的PPG原始数据片段,如果任意一个导联的ECG原始数据片段不能通过上述预处理过程被保留,或者PPG原始数据片段不能通过上述预处理过程被保留,则该组数据片段被整体丢弃。只有一组中的多个导联的ECG原始数据片段及其对应的PPG原始数据片段中的每一者均通过上述预处理过程被保留,才能得到该组的满足要求的预定时长的多个导联的ECG原始数据片段及其对应的PPG原始数据片段。
4)针对每组预定时长的多个导联的ECG预处理数据片段及其对应的PPG预处理数据片段,确定标记该多个导联的ECG预处理数据片段是否异常的标签,将该标签与该PPG预处理数据片段组成数据对,从而得到包括多组PPG数据以及标记所述PPG数据对应的ECG数据是否异常的标签的数据对的训练数据集。
由于MIMIC-III WDB中没有给出ECG标签,示例地,可以使用已知的具有较高准确性的ECG分类模型,例如AutoNet来高效地标记ECG预处理数据片段是否异常的标签。参见图10和图11所示,图10中示出了被标记为正常的ECG预处理数据片段及其对应的PPG预处理数据片段的示例,图11中示出了被标记为异常的ECG预处理数据片段及其对应的PPG预处理数据片段的示例。在其他替代示例中,也可以由专业技术人员,例如临床医生来手动标记ECG预处理数据片段是否异常的标签,以确保标签的准确性。如此,可以得到多组PPG预处理数据片段-标记对应的预处理数据片段ECG是否异常的标签这样的数据对,以构成所需的训练数据集。
在构建好上述训练数据集后,即可利用该训练数据集对上述ECG异常监测模型进行训练以确定训练好的ECG异常监测模型。
在一个实施例中,ECG异常监测模型的训练过程包括:
获取训练数据集;训练数据集包括多组PPG数据-标记对应的ECG数据是否异常的标签的数据对;
将训练数据集划分为训练集、验证集和测试集;
利用训练集对ECG异常监测模型进行训练以得到拟合的模型参数,利用验证集对拟合的模型参数确定的ECG异常监测模型进行验证并调节模型超参数,利用测试集对调节后的超参数和训练得到模型参数确定的训练好的ECG异常监测模型进行性能评估。
示例地,可以将上述一个具体示例的DCA-Net模型应用基于MIMIC-III WDB的处理的PPG数据片段,使用标记对应的ECG数据片段是否异常的标签进行二进制分类,来进行训练。将38320组PPG数据片段-标记对应的ECG数据片段是否异常的标签的数据对分成训练集、验证集和测试集,比例为64%:16%:20%。使用训练集学习模型参数,并在验证集中进行验证,以避免过度拟合。测试集仅用于报告结果。使用特异性、灵敏度/召回率、精确率、准确率、AUROC(接收器操作特征曲线下面积)和AUPRC(精确率-召回率曲线下面积)评价了模型性能。进一步对模型进行了10次随机重复,并采用了上述评估指标的平均值和标准差。
DCA-Net模型可以使用PyTorch 1.11.0版和CUDA 11.7版实施。损失函数为二元交叉熵。进一步实验通过增加损失权重来调整类失衡。模型由Adam优化器使用默认超参数设置进行优化,并使用64的批次大小进行训练。应用了步幅方面的学习率调度器,以将学习率降低到每10历元0.1倍。还应用了提前停止——如果验证损失与5个连续历元的当前最低损失相比没有减少,则终止训练。始终保存损失最低的模型,并将其应用于测试集以进行评估。
应该理解的是,虽然图2的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图2中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
性能测试
为了验证本申请提出的DCA-Net模型在架构上的优越性能,将本申请上述具体示例的DCA-Net模型与现有的神经网络模型ResNet-18、ResNet-34、ResNet-50、ECA-Net同样用于基于PPG数据监测ECG异常,以进行性能测试与比较,比较结果如下表1所示。表1中,所显示的数字是随机重复10次的平均分数百分比,其标准差在括号中显示。最佳平均准确率、AUROC和AUPRC用粗体突出显示。
表1 本申请一个具体示例的DCA-Net模型与现有神经网络模型的性能比较结果。
从上表可见,本申请提出的DCA-Net具有最佳的平均准确率、AUROC和AUPRC评分。在单侧T检验中,DCA-Net的AUPRC显著优于ECA-Net,其中p值= 0.01。此外,由于所有评估测量的标准偏差较小,DCA-Net的性能比ECA-Net更稳定。
在一个实施例中,如图12所示,提供一种ECG异常监测装置1200,包括:
数据获取模块1210,用于获取用户的PPG数据;
模型监测模块1220,用于将PPG数据输入预先训练好的ECG异常监测模型中,以得到ECG异常监测模型的输出值;其中,ECG异常监测模型是利用包括PPG数据以及标记PPG数据对应的ECG数据是否异常的标签的数据对的训练数据集训练得到的,ECG异常监测模型用于基于输入的PPG数据,输出指示所输入的PPG数据对应的ECG数据是否估计存在异常的输出值;
异常警报模块1230,用于当ECG异常监测模型的输出值指示ECG数据估计存在异常时,使与用户相关联的用户设备发出警报。
在一个实施例中,ECG异常监测装置1200还包括模型训练模块1240,模型监测模块1220还用于从模型训练模块1240获取训练好的ECG异常监测模型,该模型训练模块1240用于对ECG异常监测模型进行训练,以得到训练好的ECG异常监测模型。
示例地,上述数据获取模块1210、模型监测模块1220、异常警报模块1230可以设置于终端102中,而上述模型训练模块1300则可以设置于服务器104中。
关于ECG异常监测装置1200的具体限定可以参见上文中对于ECG异常监测方法的限定,在此不再赘述。上述ECG异常监测装置1200中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中,如图13所示,提供一种ECG异常监测设备1300,包括PPG数据探测器1310、存储器1320和处理器1330;
PPG数据探测器1310用于监测用户的PPG数据,并将测得的PPG数据传输至处理器;
存储器1320存储有计算机程序;并且
处理器1330执行计算机程序时,用于从PPG数据探测器1310接收PPG数据,并实现如上任一实施例的方法。
在一个实施例中,PPG数据探测器可以是脉搏血氧计。
在一个实施例中,ECG异常监测设备1300还包括输出装置1340,输出装置1340可以包括显示屏、指示灯、振动器和/或声音输出装置。显示屏可以是液晶显示屏或者电子墨水显示屏等。声音输出装置可以是例如扬声器、蜂鸣器等具有声音输出功能的装置。如此,当利用ECG异常监测设备1300自身来发出警报时,可以通过该输出装置1340来发出警报。
在一个实施例中,ECG异常监测设备1300还具有拨打电话功能,例如ECG异常监测设备1300还可以包括具有拨打电话功能的拨号装置,该拨号装置例如可以通过与用户关联的智能手表、手机或服务器等来实现,以在需要向预定警报接收设备拨打警报电话时,通过该拨号装置来拨打该警报电话。
在一个实施例中,ECG异常监测设备1300还具有ECG测量功能,例如ECG异常监测设备1300还可以包括ECG测量装置,以在利用PPG数据估计存在ECG异常时,进一步利用ECG测量装置对用户进行ECG测量。
可以理解,ECG异常监测设备1300可以是能够执行本申请上述实施例方法的任意设备,其可能具有各种不同的形态,例如其可以包括前述的终端102,或者其也可以包括前述的终端102和服务器104。
在一个实施例中,提供一种计算机可读存储介质,其上存储有计算机程序,其特征在于,计算机程序被处理器执行时实现以下步骤:
获取用户的PPG数据;
将PPG数据输入预先训练好的ECG异常监测模型中,以得到ECG异常监测模型的输出值;其中,ECG异常监测模型是利用包括PPG数据以及标记PPG数据对应的ECG数据是否异常的标签的数据对的训练数据集训练得到的,ECG异常监测模型用于基于输入的PPG数据,输出指示所输入的PPG数据对应的ECG数据是否估计存在异常的输出值;
当ECG异常监测模型的输出值指示ECG数据估计存在异常时,使与用户相关联的用户设备发出警报。
在其他实施例中,计算机程序被处理器执行时还实现如上任一实施例的方法的步骤。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink) DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (14)
1.一种心电图异常监测方法,其特征在于,所述方法包括:
获取用户的光电容积脉搏波数据;
将所述光电容积脉搏波数据输入预先训练好的心电图异常监测模型中,以得到所述心电图异常监测模型的输出值;其中,所述心电图异常监测模型是利用训练数据集训练得到的,所述训练数据集包括光电容积脉搏波数据,以及包括所述光电容积脉搏波数据对应的心电图数据指示是否存在心脏功能异常的标签,所述心电图异常监测模型用于基于输入的光电容积脉搏波数据,输出指示所输入的所述光电容积脉搏波数据对应的心电图数据是否估计存在心脏功能异常的输出值;
当所述心电图异常监测模型的输出值指示所述心电图数据估计存在心脏功能异常时,使与所述用户相关联的用户设备发出警报;
其中,所述心电图异常监测模型包括串行设置的一个或多个残差-双卷积注意力块,所述一个或多个残差-双卷积注意力块中的每个残差-双卷积注意力块是在残差块中加入双卷积注意力块而形成,所述双卷积注意力块包括用于在通道维度上进行卷积的第一卷积注意力块以及用于在时间维度上进行卷积的第二卷积注意力块,所述一个或多个残差-双卷积注意力块中的每个残差-双卷积注意力块包括:
第一主路径,所述第一主路径上串行设置有一个或多个卷积层以及所述双卷积注意力块;其中,所述双卷积注意力块包括串行或并行设置的所述第一卷积注意力块和所述第二卷积注意力块;以及
第一短路分支路径,所述第一短路分支路径并联于所述第一主路径上串行设置的所述一个或多个卷积层以及双卷积注意力块的两端;
其中,所述一个或多个卷积层以及双卷积注意力块中的最前者以及所述第一短路分支路径接收先前路径的输出,所述一个或多个卷积层以及双卷积注意力块中的后一者接收前一者的输出,所述一个或多个卷积层以及双卷积注意力块中的最后者的输出与所述第一短路分支路径的输出相加后馈入后续路径。
2.根据权利要求1所述的心电图异常监测方法,其特征在于,所述获取用户的光电容积脉搏波数据,将所述光电容积脉搏波数据输入预先训练好的心电图异常监测模型中,包括:
获取用户的光电容积脉搏波原始数据;
对所述光电容积脉搏波原始数据进行预处理,以得到光电容积脉搏波预处理数据;
将所述光电容积脉搏波预处理数据输入预先训练好的心电图异常监测模型中。
3.根据权利要求1所述的心电图异常监测方法,其特征在于,所述第一卷积注意力块和所述第二卷积注意力块串行设置,且所述第一卷积注意力块位于所述第二卷积注意力块之前。
4.根据权利要求1所述的心电图异常监测方法,其特征在于,所述第一卷积注意力块和所述第二卷积注意力块中的每一者包括:
第二主路径,所述第二主路径上设置有最大池化层和平均池化层、卷积层以及激活函数层,所述最大池化层和平均池化层并行设置,所述并行设置的最大池化层和平均池化层依次与卷积层和激活函数层串行设置;以及
第二短路分支路径,所述第二短路分支路径并联于所述第二主路径上设置的所述最大池化层和平均池化层、卷积层以及激活函数层的两端;
其中,所述第二短路分支路径、所述最大池化层和所述平均池化层分别接收先前路径的输出,所述最大池化层的输出与所述平均池化层的输出联合所得的输出馈入卷积层,所述卷积层的输出馈入激活函数层,所述激活函数层的输出与所述第二短路分支路径的输出相乘后馈入后续路径。
5.根据权利要求1所述的心电图异常监测方法,其特征在于,每个残差-双卷积注意力块中的各个卷积层的输出侧还分别设置有批量归一化层以及激活函数层。
6.根据权利要求1所述的心电图异常监测方法,其特征在于,所述心电图异常监测模型还包括:
输入卷积层和池化层,位于所述多个残差-双卷积注意力块之前,且与所述多个残差-双卷积注意力块串行设置,以及
平均池化层和全连接层,位于所述多个残差-双卷积注意力块之后,且与所述多个残差-双卷积注意力块串行设置;
其中,所述输入卷积层接收对所述心电图异常监测模型输入的光电容积脉搏波数据,所述输入卷积层的输出馈入池化层,所述池化层的输出馈入所述多个残差-双卷积注意力块中的首个残差-双卷积注意力块,所述多个残差-双卷积注意力块中的最后一个残差-双卷积注意力块的输出馈入所述平均池化层,所述平均池化层的输出馈入所述全连接层,所述全连接层输出所述心电图异常监测模型的输出值。
7.根据权利要求1至6中任一项所述的心电图异常监测方法,其特征在于,所述心电图异常监测模型的训练过程包括:
获取训练数据集;所述训练数据集包括多组光电容积脉搏波数据-标记对应的心电图数据指示是否存在心脏功能异常的标签的数据对;
将所述训练数据集划分为训练集、验证集和测试集;
利用所述训练集对所述心电图异常监测模型进行训练以得到拟合的模型参数,利用所述验证集对拟合的所述模型参数确定的心电图异常监测模型进行验证并调节模型超参数,利用所述测试集对调节后的超参数和训练得到的模型参数确定的训练好的心电图异常监测模型进行性能评估。
8.根据权利要求1至6中任一项所述的心电图异常监测方法,其特征在于,与所述用户相关联的用户设备包括所述用户佩戴的便携终端设备和/或与所述用户关联的监护者的监护终端设备,所述获取用户的光电容积脉搏波数据包括由所述便携终端设备实时监测所述用户的光电容积脉搏波数据;
所述使与所述用户相关联的用户设备发出警报包括以下各项中的一项或多项:在所述用户设备的显示屏上显示指示所述用户的心电图数据估计存在心脏功能异常的视觉信息、使所述用户设备的指示灯发出指示所述用户的心电图数据估计存在心脏功能异常的亮灯提示、使所述用户设备的振动器做出指示所述用户的心电图数据估计存在心脏功能异常的振动提示、使所述用户设备的声音输出装置发出指示所述用户的心电图数据估计存在心脏功能异常的声音提示。
9.根据权利要求1至6中任一项所述的心电图异常监测方法,其特征在于,所述方法还包括:
当所述心电图异常监测模型的输出值指示所述心电图数据估计存在心脏功能异常时,向预定警报接收设备拨打警报电话。
10.根据权利要求1至6中任一项所述的心电图异常监测方法,其特征在于,所述方法还包括:
当所述心电图异常监测模型的输出值指示所述心电图数据估计存在心脏功能异常时,对所述用户进行心电图测量,以得到指示所述用户的心电图是否指示存在心脏功能异常的测量结果。
11.一种心电图异常监测装置,其特征在于,所述装置包括:
数据获取模块,用于获取用户的光电容积脉搏波数据;
模型监测模块,用于将所述光电容积脉搏波数据输入预先训练好的心电图异常监测模型中,以得到所述心电图异常监测模型的输出值;其中,所述心电图异常监测模型是利用训练数据集训练得到的,所述训练数据集包括光电容积脉搏波数据,以及包括所述光电容积脉搏波数据对应的心电图数据指示是否存在心脏功能异常的标签,所述心电图异常监测模型用于基于输入的光电容积脉搏波数据,输出指示所输入的所述光电容积脉搏波数据对应的心电图数据是否估计存在心脏功能异常的输出值;其中,所述心电图异常监测模型包括串行设置的一个或多个残差-双卷积注意力块,所述一个或多个残差-双卷积注意力块中的每个残差-双卷积注意力块是在残差块中加入双卷积注意力块而形成,所述双卷积注意力块包括用于在通道维度上进行卷积的第一卷积注意力块以及用于在时间维度上进行卷积的第二卷积注意力块,所述一个或多个残差-双卷积注意力块中的每个残差-双卷积注意力块包括第一主路径以及第一短路分支路径,所述第一主路径上串行设置有一个或多个卷积层以及所述双卷积注意力块;其中,所述双卷积注意力块包括串行或并行设置的所述第一卷积注意力块和所述第二卷积注意力块;所述第一短路分支路径并联于所述第一主路径上串行设置的所述一个或多个卷积层以及双卷积注意力块的两端;其中,所述一个或多个卷积层以及双卷积注意力块中的最前者以及所述第一短路分支路径接收先前路径的输出,所述一个或多个卷积层以及双卷积注意力块中的后一者接收前一者的输出,所述一个或多个卷积层以及双卷积注意力块中的最后者的输出与所述第一短路分支路径的输出相加后馈入后续路径;
异常警报模块,用于当所述心电图异常监测模型的输出值指示所述心电图数据估计存在心脏功能异常时,使与所述用户相关联的用户设备发出警报。
12.一种心电图异常监测设备,包括光电容积脉搏波数据探测器、存储器和处理器;
所述光电容积脉搏波数据探测器用于监测用户的光电容积脉搏波数据,并将监测到的所述光电容积脉搏波数据传输至所述处理器;
所述存储器存储有计算机程序;并且
所述处理器执行所述计算机程序时,用于从所述光电容积脉搏波数据探测器接收所述光电容积脉搏波数据,并实现权利要求1至10中任一项所述的心电图异常监测方法。
13.根据权利要求12所述的心电图异常监测设备,其特征在于,还包括输出装置,所述输出装置包括显示屏、指示灯、振动器和/或声音输出装置。
14.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至10中任一项所述的心电图异常监测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310881402.6A CN116584952B (zh) | 2023-07-18 | 2023-07-18 | 心电图异常监测方法、装置、设备和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310881402.6A CN116584952B (zh) | 2023-07-18 | 2023-07-18 | 心电图异常监测方法、装置、设备和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116584952A CN116584952A (zh) | 2023-08-15 |
CN116584952B true CN116584952B (zh) | 2023-12-08 |
Family
ID=87612065
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310881402.6A Active CN116584952B (zh) | 2023-07-18 | 2023-07-18 | 心电图异常监测方法、装置、设备和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116584952B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111839495A (zh) * | 2020-07-30 | 2020-10-30 | 深圳前海微众银行股份有限公司 | 检测方法、设备和存储介质 |
KR102373778B1 (ko) * | 2021-07-29 | 2022-03-15 | 서울대학교산학협력단 | 심전도 데이터로부터 광용적맥파 데이터를 생성하는 장치 및 방법 |
CN114391824A (zh) * | 2022-01-18 | 2022-04-26 | 安徽大学 | 一种基于ppg信号的心率估计方法、电子设备及存储介质 |
CN115270848A (zh) * | 2022-06-17 | 2022-11-01 | 合肥心之声健康科技有限公司 | 一种ppg与ecg自动转换智能算法、存储介质和计算机系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3669768B1 (en) * | 2018-12-20 | 2024-02-21 | IMEC vzw | A method for generating a model for generating a synthetic ecg and a method and system for analysis of heart activity |
-
2023
- 2023-07-18 CN CN202310881402.6A patent/CN116584952B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111839495A (zh) * | 2020-07-30 | 2020-10-30 | 深圳前海微众银行股份有限公司 | 检测方法、设备和存储介质 |
KR102373778B1 (ko) * | 2021-07-29 | 2022-03-15 | 서울대학교산학협력단 | 심전도 데이터로부터 광용적맥파 데이터를 생성하는 장치 및 방법 |
CN114391824A (zh) * | 2022-01-18 | 2022-04-26 | 安徽大学 | 一种基于ppg信号的心率估计方法、电子设备及存储介质 |
CN115270848A (zh) * | 2022-06-17 | 2022-11-01 | 合肥心之声健康科技有限公司 | 一种ppg与ecg自动转换智能算法、存储介质和计算机系统 |
Also Published As
Publication number | Publication date |
---|---|
CN116584952A (zh) | 2023-08-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108577830B (zh) | 一种面向用户的体征信息动态监护方法和动态监护系统 | |
WO2019161608A1 (zh) | 多参数监护数据分析方法和多参数监护系统 | |
McManus et al. | A novel application for the detection of an irregular pulse using an iPhone 4S in patients with atrial fibrillation | |
US20170215756A1 (en) | Devices and methods for real-time denoising of electrocardiograms | |
US20240164685A1 (en) | Health monitoring device | |
CN103038772B (zh) | 预测患者的存活性的系统和装置 | |
WO2019161609A1 (zh) | 多参数监护数据分析方法和多参数监护仪 | |
EP1075214A1 (en) | Reduced lead set for detecting acute cardiac ischemic conditions | |
Pipitprapat et al. | The validation of smartphone applications for heart rate measurement | |
US11103145B1 (en) | Physiological signal monitoring and apparatus therefor | |
US20210338174A1 (en) | Method and system for assessing emergency risk for patients | |
US11179046B2 (en) | Method and system for detection of atrial fibrillation | |
Pinge et al. | A comparative study between ECG-based and PPG-based heart rate monitors for stress detection | |
US20210307624A1 (en) | Non-invasive device and methods for monitoring muscle tissue condition | |
EP3107447B1 (en) | Method and system for evaluating a noise level of a biosignal | |
CN114246569B (zh) | 检测早搏类型的方法、装置及介质 | |
Tanasković et al. | A new algorithm for fetal heart rate detection: Fractional order calculus approach | |
CN116584952B (zh) | 心电图异常监测方法、装置、设备和存储介质 | |
Shamim et al. | Diagnostic accuracy of smartphone-connected electrophysiological biosensors for prediction of blood glucose level in a type-2 diabetic patient using machine learning: A pilot study | |
US20230200746A1 (en) | System for acquisition and analysis of maternal and/or fetal physiological signals | |
US20230009430A1 (en) | Systems and methods to detect cardiac events | |
CN116548935A (zh) | 基于柔性有机光探测器和深度学习算法的血压测量系统 | |
Wisana et al. | Smartband for Heartbeat and Oxygen Saturation Monitoring with Critical Warning to Paramedic via IoT | |
Hassanuzzaman et al. | End to end solution for continuous monitoring and real-time analysis of vital signs from ecg signal | |
Skrivanos et al. | Home Healthcare Technologies and Services: Heart-Rate Fetus Monitoring System Using an MCU ESP8266 Node |
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 |