CN107440679B - 用于检测心电图异常的方法及对应系统 - Google Patents

用于检测心电图异常的方法及对应系统 Download PDF

Info

Publication number
CN107440679B
CN107440679B CN201611227015.7A CN201611227015A CN107440679B CN 107440679 B CN107440679 B CN 107440679B CN 201611227015 A CN201611227015 A CN 201611227015A CN 107440679 B CN107440679 B CN 107440679B
Authority
CN
China
Prior art keywords
beats
beat
detected
heart rate
test
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
CN201611227015.7A
Other languages
English (en)
Other versions
CN107440679A (zh
Inventor
B·罗西
P·弗拉尼托
D·赞邦
D·卡雷拉
G·博拉基
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.)
STMicroelectronics SRL
Original Assignee
STMicroelectronics SRL
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 STMicroelectronics SRL filed Critical STMicroelectronics SRL
Publication of CN107440679A publication Critical patent/CN107440679A/zh
Application granted granted Critical
Publication of CN107440679B publication Critical patent/CN107440679B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7221—Determining signal validity, reliability or quality
    • 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
    • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02438—Detecting, measuring or recording pulse rate or heart rate with portable devices, e.g. worn by the patient
    • 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
    • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
    • A61B5/0245—Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/1118—Determining activity level
    • 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/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle
    • A61B5/35—Detecting specific parameters of the electrocardiograph cycle by template matching
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle
    • A61B5/364—Detecting abnormal ECG interval, e.g. extrasystoles, ectopic heartbeats
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle
    • A61B5/366—Detecting abnormal QRS complex, e.g. widening
    • 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/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
    • 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/7232—Signal processing specially adapted for physiological signals or for diagnostic purposes involving compression of the physiological signal, e.g. to extend the signal recording period
    • 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/7253—Details of waveform analysis characterised by using transforms
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2505/00—Evaluating, monitoring or diagnosing in the context of a particular type of medical care
    • A61B2505/07—Home care
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/02—Details of sensors specially adapted for in-vivo measurements
    • A61B2562/0219—Inertial sensors, e.g. accelerometers, gyroscopes, tilt switches
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0002—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network
    • A61B5/0015—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network characterised by features of the telemetry system
    • A61B5/0022—Monitoring a patient using a global network, e.g. telephone networks, internet
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle
    • A61B5/352—Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A61B5/7278—Artificial waveform generation or derivation, e.g. synthesising signals from measured signals
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/70—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Cardiology (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Physiology (AREA)
  • Artificial Intelligence (AREA)
  • Signal Processing (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

心率监测器检测测试信号中的心跳。局部心率和加速度的能量与检测到的心跳相关联。基于与相应心跳相关联的局部心率和加速度的能量，将检测到的心跳包括在测试心跳集合中或将其排除在外。使用稀疏近似模型来检测心跳测试集合中的异常心跳。心率监测器可以检测训练心跳信号中的心跳。参考心率和加速度的能量与训练心跳信号的检测到的搏动相关联，并且基于与训练心跳信号中的检测到的搏动相关联的心率和加速度的能量而选择性地包括在训练数据集合中。可以使用训练数据集合来生成稀疏表示模型的字典。

Description

用于检测心电图异常的方法及对应系统
技术领域
本公开涉及用于检测心电图(ECG)信号(诸如由用于健康监测的可穿戴设备获取的信号)中的异常的技术。
背景技术
ECG信号中的异常检测是指找到不符合所监测的信号的预期行为的模式的问题。这些异常模式可以转化为显著信息，诸如心脏活动异常(心律失常)、感测过程中的故障等。
可穿戴设备可以连续地获取、存储和周期性地传输通常由主机设备实时处理的生理信号。ECG信号的形态特点可能对于不同的用户而显示出显著的变化，因此通用方法在可穿戴设备上可能表现不佳。例如，每次可穿戴设备被重新定位或受到意外移动时，ECG形态可能被修改。对于日常活动，心率和心跳的形态也可以改变。
发明内容
在正常条件下获得的ECG心跳的特征在于特定的形态结构。另一方面，不符合这种结构的ECG获取可以被认为是异常。在一个实施例中，假设正常ECG心跳由稀疏模型很好地描述，诸如基于稀疏表示理论的模型。假设正常心跳s可以被建模为s＝D·x+e，其中，D是未知字典(dictionary)，x是稀疏系数向量，e是来自感测过程的随机噪声。在实施例中，以数据驱动方式(例如，直接从可穿戴设备获取的数据)来学习模型。在学习阶段期间，确定提供感测的训练数据集合的稀疏表示的字典D。在测试阶段中，对于每个获取的心跳s，计算稀疏系数向量x，使得重建误差||s-D·x||2最小(粗略地说，评估这种心跳对学习模型的拟合优度)。对于异常心跳，假设观察到的重建误差将超过所容许的噪声能量。
在实施例中，一种方法，包括：检测心跳信号中的搏动；将心率和加速度的能量与检测到的搏动相关联；基于与检测到的搏动相关联的心率和加速度的能量来将检测到的搏动选择性地包括检测到的搏动集合中；和执行以下各项中的至少一项：使用所检测到的搏动集合作为训练数据集合来生成稀疏近似(sparse approximation)模型的字典；以及使用稀疏近似模型的字典来检测检测到的搏动集合中的异常搏动。在实施例中，该方法包括：将基于与检测到的搏动集合中的搏动相关联的心率和加速度的能量的搏动变换(beattransform)应用于所选择的搏动集合中的搏动，生成自适应的搏动(adapted beat)；基于与检测到的搏动集合中的搏动相关联的心率和加速度的能量来将字典变换应用于字典，生成自适应的字典(adapted dictionary)；生成相对于自适应的字典的自适应的搏动的稀疏近似；和基于所生成的稀疏近似来将检测到的搏动集合中的搏动分类为异常。在实施例中，基于稀疏近似来将搏动集合中的搏动分类为异常包括：确定与自适应的搏动的稀疏近似相关联的均方根误差；将均方根误差与阈值误差进行比较；以及基于比较来将搏动集合中的搏动分类为异常或正常。在实施例中，该方法包括：响应于将搏动集合中的搏动分类为异常，发起与搏动集合中的搏动有关的数据的传输。在实施例中，该方法包括：响应于搏动集合中的搏动被分类为异常，存储与搏动集合中的搏动有关的数据。在实施例中，心率是局部心率，并且将检测到的搏动选择性地包括在检测到的搏动集合中包括：确定检测到的搏动的准时心率是否是相对于与检测到的搏动相关联的局部心率的外点心率；确定检测到的搏动是否是相对于与检测到的搏动相关联的加速度的能量的外点搏动；当准时心率是相对于局部心率的外点心率并且检测到的搏动是相对于与检测到的搏动相关联的加速度的能量的外点搏动时，从检测到的搏动集合中排除检测到的搏动；当准时心率是相对于局部心率的内点心率并且检测到的搏动是相对于与检测到的搏动相关联的加速度的能量的内点搏动时，将所检测到的搏动包括在检测到的搏动集合中；和当准时心率是相对于局部心率的外点心率并且检测到的搏动是相对于与所检测到的搏动相关联的加速度的能量相关联的内点搏动时，从检测到的搏动集合中排除检测到的搏动并且将其分类为节律性异常搏动。在实施例中，该方法包括：将基于与检测到的搏动集合中的搏动相关联的心率和加速度的能量的搏动变换应用于所选择的搏动集合中的搏动，生成自适应的搏动；将基于与检测到的搏动集合中的搏动相关联的心率和加速度的能量的字典变换应用于字典，生成自适应的字典；生成相对于自适应的字典的自适应的搏动的稀疏近似；和基于自适应的搏动的所生成的稀疏近似来将检测到的搏动集合中的搏动分类为异常。在实施例中，该方法包括通过迭代地进行以下各项而使用检测到的搏动集合作为训练数据集合来生成稀疏表示字典：生成稀疏表示矩阵的系数；和优化字典和稀疏表示矩阵的系数。在实施例中，心率是与检测到的搏动集合相关联的参考心率，并且将检测到的搏动选择性地包括在检测到的搏动集合中包括：确定检测到的搏动的准时心率是否是相对于参考心率的外点心率；确定检测到的搏动是否是相对于与检测到的搏动相关联的加速度的能量的外点搏动；当准时心率是相对于参考心率的外点心率时，从检测到的搏动集合中排除检测到的搏动；当检测到的搏动是相对于与检测到的搏动相关联的加速度的能量的外点搏动时，从检测到的搏动集合中排除检测到的搏动。在实施例中，该方法包括：使用检测到的搏动集合作为训练数据集合来生成稀疏近似模型的字典，其中，心跳信号是训练心跳信号；检测测试心跳信号中的搏动；将局部心率和加速度的能量与测试心跳信号中的检测到的搏动相关联；基于与测试心跳信号中的检测到的搏动相关联的局部心率和加速度的能量来将检测到的搏动选择性地包括在测试搏动集合中；和使用稀疏近似模型的字典来检测测试搏动集合中的异常搏动。在实施例中，设备被配置成执行本文中所公开的方法中的一个或多个方法。
在实施例中，一种设备，包括：接口，被配置成接收测试心跳信号、以及一个或多个加速度信号；和信号处理电路，被配置成：检测测试心跳信号中的搏动；将心率和加速度的能量与检测到的搏动相关联；基于与检测到的搏动相关联的心率和加速度的能量来将检测到的搏动选择性地包括在测试搏动集合中；和使用稀疏近似模型的字典来检测测试搏动集合中的异常搏动。在实施例中，信号处理电路在操作中：将基于与测试搏动集合中的搏动相关联的心率和加速度的能量的搏动变换应用于所选择的搏动集合中的搏动，生成自适应的搏动；将基于与测试搏动集合中的搏动相关联的心率和加速度的能量的字典变换应用于字典，生成自适应的字典；生成相对于自适应的字典的自适应的搏动的稀疏近似；和基于所生成的稀疏近似来生成指示测试搏动集合中的搏动是否异常的信号。在实施例中，基于所生成的稀疏近似来生成指示测试搏动集合中的搏动是否异常的信号包括：确定与自适应的搏动的所生成的稀疏近似值相关联的均方根误差；将均方根误差与阈值误差进行比较；和基于比较来生成信号。在实施例中，当所生成的信号指示测试搏动集合中的搏动是异常的时，信号处理电路发起与测试搏动集合中的搏动有关的数据的传输。在实施例中，心率是局部心率，并且将检测到的搏动选择性地包括在测试搏动集合中包括：确定检测到的搏动的准时心率是否是相对于与检测到的搏动相关联的局部心率的外点心率；确定检测到的搏动是否是相对于与检测到的搏动相关联的加速度的能量的外点搏动；当准时心率是相对于局部心率的外点心率或检测到的搏动是相对于与检测到的搏动相关联的加速度的能量的外点搏动时，从检测到的搏动集合中排除检测到的搏动；和当准时心率是相对于局部心率的内点心率并且检测到的搏动是相对于与检测到的搏动相关联的加速度的能量的内点搏动时，将检测到的搏动包括在检测到的搏动集合中。在实施例中，接口被配置成接收训练心跳信号；并且信号处理电路被配置成：检测训练心跳信号的搏动；将参考心率和加速度的能量与训练心跳信号的检测到的搏动相关联；基于与训练心跳信号中的检测到的搏动相关联的心率和加速度的能量来将训练心跳信号的检测到的搏动选择性地包括在训练数据集合中；和使用训练数据集合来生成稀疏表示模型的字典。在实施例中，信号处理电路通过迭代地进行以下各项来生成稀疏表示模型的字典：生成稀疏表示矩阵的系数；以及优化字典和稀疏表示矩阵的系数。在实施例中，参考心率与训练心跳信号相关联，并且将训练心跳信号的检测到的搏动选择性地包括在训练数据集合中包括：确定训练心跳信号的检测到的搏动的准时心率是否是相对于参考心率的外点心率；确定训练心跳信号的检测到的搏动是否是相对于与训练心跳信号的检测到的搏动相关联的加速度的能量的外点搏动；当准时心率是相对于参考心率的外点心率时，从训练数据集合中排除训练心跳信号的检测到的搏动；和当训练心跳信号的检测到的搏动是相对于与训练心跳信号的检测到的搏动相关联的加速度的能量的外点搏动时，从训练数据集合中排除训练心跳信号的检测到的搏动。
在实施例中，一种系统，包括：传感器，被配置成生成心跳信号；加速度计，被配置成生成一个或多个加速度信号；以及信号处理电路，被配置成：检测测试心跳信号的搏动；将心率和加速度的能量与测试心跳信号的检测到的搏动相关联；基于与测试心跳信号的检测到的搏动相关联的心率和加速度的能量来将测试心跳信号的检测到的搏动选择性地包括在测试搏动集合中；和使用稀疏近似模型的字典来检测测试搏动集合中的异常搏动。在实施例中，信号处理电路在操作中：将基于与测试搏动集合中的搏动相关联的心率和加速度的能量的搏动变换应用于所选择的搏动集合中的搏动，生成自适应的搏动；将基于与测试搏动集合中的搏动相关联的心率和加速度的能量的字典变换应用于字典，生成自适应的字典；生成相对于自适应的字典的自适应的搏动的稀疏近似；和基于所生成的稀疏近似来生成指示测试搏动集合中的搏动是否异常的信号。在实施例中，信号处理电路被配置成：检测训练心跳信号的搏动；将参考心率和加速度的能量与训练心跳信号的检测到的搏动相关联；基于与训练心跳信号中的检测到的搏动相关联的心率和加速度的能量将训练心跳信号的检测到的搏动选择性地包括在训练数据集合中；和使用训练数据集合来生成稀疏表示模型的字典。在实施例中，信号处理电路通过迭代地进行以下各项来生成稀疏表示模型的字典：生成稀疏表示矩阵的系数；和优化字典和稀疏表示矩阵的系数。
在实施例中，一种非暂态计算机可读介质，其内容配置心率监测设备以执行一种方法，该方法包括：检测测试心跳信号的搏动；将心率和加速度的能量与测试心跳信号的检测到的搏动相关联；基于与测试心跳信号的检测到的搏动相关联的心率和加速度的能量来将测试心跳信号的检测到的搏动选择性地包括在测试搏动集合中；和使用稀疏近似模型的字典来检测测试搏动集合中的异常搏动。在实施例中，该方法包括：将基于与测试搏动集合中的搏动相关联的心率和加速度的能量的搏动变换应用于所选择的搏动集合中的搏动，生成自适应的搏动；将基于与测试搏动集合中的搏动相关联的心率和加速度的能量的字典变换应用于字典，生成自适应的字典；生成相对于自适应的字典的自适应的搏动的稀疏近似；和基于所生成的稀疏近似来生成指示测试搏动集合中的搏动是否异常的信号。在实施例中，该方法包括：检测训练心跳信号的搏动；将参考心率和加速度的能量与训练心跳信号的检测到的搏动相关联；基于与训练心跳信号中的检测到的搏动相关联的心率和加速度的能量来将训练心跳信号的检测到的搏动选择性地包括在训练数据集合中；和使用训练数据集合来生成稀疏表示模型的字典。在实施例中，该方法包括：通过迭代地进行以下各项来生成稀疏表示模型的字典：生成稀疏表示矩阵的系数；和优化字典和稀疏表示矩阵的系数。
在实施例中，一种系统，包括：用于生成心跳信号的装置；用于生成加速度信号的装置；以及用于通过以下各项来处理测试心跳信号的装置：检测测试心跳信号的搏动；将心率和加速度的能量与测试心跳信号的检测到的搏动相关联；基于与测试心跳信号的检测到的搏动相关联的心率和加速度的能量来将测试心跳信号的检测到的搏动选择性地包括在测试搏动集合中；和使用稀疏近似模型的字典来检测测试搏动集合中的异常搏动。在实施例中，处理测试心跳信号包括：将基于与测试搏动集合中的搏动相关联的心率和加速度的能量的搏动变换应用于所选择的搏动集合中的搏动，生成自适应的搏动；将基于与测试搏动集合中的搏动相关联的心率和加速度的能量的字典变换应用于字典，生成自适应的字典；生成相对于自适应的字典的自适应的搏动的稀疏近似；和基于所生成的稀疏近似来生成指示测试搏动集合中的搏动是否异常的信号。在实施例中，该系统包括：用于通过以下各项来处理训练心跳信号的装置：检测训练心跳信号的搏动；将参考心率和加速度的能量与训练心跳信号的检测到的搏动相关联；基于与训练心跳信号中的检测到的搏动相关联的心率和加速度的能量来将训练心跳信号的检测到的搏动选择性地包括在训练数据集合中；和使用训练数据集合来生成稀疏表示模型的字典。在实施例中，处理训练心跳信号包括：通过迭代地进行以下各项来生成稀疏表示模型的字典：生成稀疏表示矩阵的系数；和优化字典和稀疏表示矩阵的系数。
附图说明
图1图示了正常心跳的示例形态结构。
图2是实施例中的正交匹配追踪(OMP)算法的迭代中的信号向量、稀疏表示系数向量和字典之间的关系的图形说明。
图3图示了应用OMP算法的方法的实施例。
图4是学习的稀疏模型的实施例的图形说明。
图5图示了应用K-SVD算法以生成字典的方法的实施例。
图6图示了检测异常心跳的方法的实施例。
图7是心跳和加速度数据的图形说明。
图8是用户状态和心率与所测量的加速度能量之间的示例关系的图形表示。
图9是应用于所测量的搏动和字典以考虑实施例中用户状态的改变的示例变换配对的图形表示。
图10图示了用于检测异常心跳的系统的实施例。
图11图示了检测异常心跳的方法的实施例。
图12图示了生成用于检测异常心跳的训练字典的子例程或方法的实施例。
图13是ECG迹线的一部分的图形表示。
图14是加速度计数据的图形表示。
图15是ECG迹线的加窗搏动的图形表示。
图16图示了使用训练字典、ECG迹线和加速度计数据来检测异常心跳的子例程或方法的实施例。
图17是ECG迹线的图形表示。
图18是从ECG迹线提取的搏动的示例的图形说明。
图19是处于休息状态的用户的示例搏动和用户在努力状态后的休息中的示例搏动的图形比较。
图20是实施例中的搏动自适应的示例应用的图形说明。
图21是实施例中的字典自适应的示例应用的图形说明。
图22是实施例中的正交匹配追踪(OMP)算法的迭代中的自适应信号向量、稀疏表示系数向量和自适应字典之间的关系的图形说明。
图23A至图23D图示了来自研究中患者的示例搏动。
图24图示了以线框图的形式呈现的研究中的异常检测器的实施例的全局行为。
图25表示研究中样本患者的结果。
具体实施方式
专利或申请文件包含至少一张彩色图。具有彩色附图的本专利或专利申请公开的副本将由办公室根据请求并且支付必要的费用来提供。
在以下描述中，阐述了某些细节以便提供对设备、系统、方法和制品的各种实施例的透彻理解。然而，本领域技术人员应当理解，可以在没有这些细节的情况下实践其它实施例。在其它实例中，与例如可穿戴设备和信号处理电路(诸如晶体管、乘法器、发射器、NFC电路、集成电路等)相关联的公知结构和方法在一些附图中未详细示出或描述以避免不必要地模糊实施例的描述。
除非上下文另有要求，否则在下面的说明书和权利要求书中，词语“包括(comprise)”及其变型(诸如“包括(comprising)”和“包括(comprises)”)应以开放的、包含的意义来解释，也就是说，“包括但不限于”。
贯穿本说明书对“一个实施例”或“实施例”的引用意味着结合该实施例描述的特定特征、结构或特点包括在至少一个实施例中。因此，贯穿本说明书在各个地方出现的短语“在一个实施例中”或“在实施例中”不一定是指相同的实施例或所有实施例。更进一步地，特定特征、结构或特点可以在一个或多个实施例中以任何合适的方式组合以获得另外的实施例。
标题仅为了方便而提供，并且不解释本公开的范围或含义。
附图中的元件的尺寸和相对位置不一定按比例绘制。例如，各种元件的形状和角度不是按比例绘制的，并且这些元件中的一些元件被放大和定位以提高附图易读性。进一步地，所绘制的元件的特定形状不一定旨在传达关于特定元件的实际形状的任何信息，并且仅为了便于在附图中识别而选择。
在正常条件下获取的正常ECG信号的特征可以在于特定的形态结构或模式。另一方面，不符合这种结构的ECG获取可以被认为是异常。正常状态的形态可能是复杂的，并且分析或近似模型可能难以开发。
代替使用分析或近似模型，本公开提出采用数据驱动的学习过程，其在实施例中是用户特定的，适应于用户状态的改变，并且通过执行无监督和鲁棒的学习阶段来容易地重新配置。
在正常状态下获取的ECG心跳的特征可以在于特定的形态结构。图1图示了正常心跳的示例性形态结构，并且包括P，Q，R，S和T波。在实施例中，从提供正常数据的稀疏表示的训练数据集合开发或学习模型。在测试阶段，可以监测重建误差，即，来自学习模型的测试信号的偏差，以便检测异常。例如，当重建误差超过阈值时，搏动可以被分类为异常。用户状态可以被表征为两个特征的组合，心率HR和加速度EE中的能量的测量。
稀疏表示可以允许通过稀疏线性组合来建模复杂信号，诸如构成字典的冗余原子集合中的仅少数系数是非零的线性组合。使用正常ECG信号的形态的稀疏表示建模可以有助于简化与检测异常心跳相关联的信号处理任务。
已经提出了非稀疏表示模型。参见Li等人在Proceedings of the ACMConference on Bioinformatics，Computational Biology and Biomedicine，第210至217页，ACM，2012中的Detecting ecg abnormalities via transductive transferlearning；Lagerholm等人的Clustering ECG complexes using hermite functions andself-organizing maps，Biomedical Engineering，IEEE Transactions on，第47卷，第7期，第838至848页，2000；Sivaraks等人的Robust and accurate Anomaly detection inECG artifacts using time series motif discovery，Comp.Math.Methods in Medicine2015。Li的缺点是所有数据必须从可穿戴设备连续传送到主机，以便分析数据。Lagerholm采用基于神经网络的无监督聚类过程，并且涉及几种不同类型的心跳的分类问题，而不是异常检测。Lagerholm也可能难以在可穿戴设备上实现。Sivaraks不支持实时检测。尽管Alder等人的“Sparse coding with anomaly detection”，Proceedings of IEEE，International Workshop on Machine Learning for Signal Processing(MLSP)，2003年9月，第1至6页，提出了学习用户的特定模型，但是Alder的途径不鲁棒，并且不适应用户状态的改变。
术语“稀疏”是指涉及其非零项的数目的向量的可测量特性。
Figure BDA0001193823700000102
-伪范数(在本文中被称为
Figure BDA0001193823700000105
-范数或运算符)测量向量的稀疏度。向量
Figure BDA0001193823700000101
的
Figure BDA0001193823700000106
-范数是x中非零项的数目。符号||x||0用于表示x的
Figure BDA0001193823700000107
-范数。如果||x||0≤T，则向量x被认为是T稀疏。
Figure BDA0001193823700000108
-范数是向量中的项的绝对值的和。
Figure BDA0001193823700000104
-范数是向量的欧几里德长度(向量的元素的绝对平方和的平方根)。Frobenius范数是矩阵的元素的绝对平方和的平方根。符号||·||1、||·||2和||·||F将分别用于
Figure BDA0001193823700000103
和Frobenius范数。
给定用于i＝1，...，N的N个向量
Figure BDA0001193823700000111
矩阵形式的字典可以被定义为：
Figure BDA0001193823700000112
其中列di被称为原子。如果M<N，则D被认为是过完备字典。
令
Figure BDA0001193823700000113
是过完备字典，并且令
Figure BDA0001193823700000114
是给定的数据向量。
Figure BDA0001193823700000115
问题可以表示为：
Figure BDA0001193823700000116
在上述问题中，如果稀疏向量
Figure BDA0001193823700000117
被找到作为解，则x被称为s的稀疏表示。一般而言，最稀疏的近似不是唯一的，但是存在其中它们可以是的条件。对于关于稀疏表示的唯一性的条件的更多细节，请参考M.Elad的Sparse and Redundant Representations：From Theory to Applications in Signal and Image Processing，SpringerPub.Co.1st Ed.2010。
等式(2)的问题的更一般版本解决了数据中的噪声的存在，其是实际应用中的现实假设。更具体地，可以假设以下模型成立：
s＝Dx+e， (3)
其中，e是噪声项，其包括测量和近似模型中的误差。在这种情况下，信号s可以通过相对于字典D的稀疏表示来近似，s≈Dx。令
Figure BDA0001193823700000118
是固定的，我们将
Figure BDA0001193823700000119
问题定义为：
Figure BDA00011938237000001110
如果稀疏向量
Figure BDA00011938237000001111
是上述问题的解，则它是s的稀疏近似。在这个问题中，x的选择具有更多的灵活性，因为不再需要精确地再现s。在
Figure BDA00011938237000001112
问题中，噪声项被包括在模型中，换句话说，考虑系统s＝Dx+e，其中，e表示噪声分量。在这种情况下，代替唯一性条件，理论结果检查给定噪声幅度||e||2.上限的解x的稳定性。参见Elad。对于情况s≈Dx的明确地迫使T稀疏约束在x上的另一个有用的公式如下：
Figure BDA00011938237000001113
其中T是已知的固定值。该公式化在本文中被称为
Figure BDA0001193823700000121
问题。如果稀疏向量
Figure BDA0001193823700000122
是
Figure BDA0001193823700000123
问题的解，则x是s的稀疏近似。在这个问题中，x的选择还具有更多灵活性，因为不再需要精确地再现s。在这种情况下，代替唯一性条件，理论结果检查给定噪声幅度||e||2上限的解x的稳定性。参见Elad。
正交匹配追踪(OMP)是用于找到
Figure BDA0001193823700000124
和
Figure BDA0001193823700000125
问题的近似解的贪婪算法。OMP通过迭代过程来构建近似。在每次迭代时，计算局部最优解(例如，对于等式5)。这可以通过找到最接近地类似于残差向量r的D中的列向量来完成。残差向量开始等于被近似的向量，r(0)＝s，并且在每次迭代时被调整以考虑先前选择的向量。这种局部最优解的序列可以导致全局最优解。通常，这不是一般情况，尽管存在这样的条件，在该条件下，结果将是最优解。OMP基于被称为匹配追踪(MP)的早期算法的变型。MP在每次迭代时简单地从残差向量中去除所选择的列向量。OMP在每次迭代时使用最小二乘法步骤来更新残差向量，以便改进近似。这些迭代的通常终止准则是||r(i)||2≤2或||x(i)||0≥T或其组合。图2是在实施例中OMP算法的迭代中的向量S和X与字典D之间的关系的图形说明。
图3图示了应用OMP算法的方法300的实施例。该方法在302处开始，并且进行到304。在304处，例如通过检索、接收或生成字典D、数据向量s、噪声边界ε和/或目标稀疏度T来定义字典D、要近似的数据向量s、可选地噪声边界ε和目标稀疏度T。噪声边界ε设定限制，使得x s.t.||Dx-s||2≤ε。目标稀疏度T设定条件||x||0≤T。该方法从304进行到306。
在306处，方法300设定残差r(0)＝s，索引i＝0和索引集合
Figure BDA0001193823700000127
该方法从306进行到308。在308处，wi被定义为wi＝l，其中，dl给出max<ri，dl>的解。该方法从308进行到310。在310处，将集合Wi更新为Wi＝Wi-1∪{wi}。该方法从310进行到312。在312处，解出最小二乘问题
Figure BDA0001193823700000126
该方法从312进行到314。在314处，新残差被计算为
Figure BDA0001193823700000131
该方法进行到316。在316处，索引i递增。该方法从316进行到318。在318处，该方法确定迭代是否应该终止。如所图示的，这是通过确定是否||x(i)||0≥T来完成。当在318处确定||x(i)||0≥T时，方法300返回到308以进行下一次迭代时。当在318处未确定||x(i)||0≥T时，方法从318进行到320，其中终止迭代并且输出当前残差向量。方法300的实施例可以执行附加动作，可以执行较少动作，可以修改所执行的动作，并且可以以不同于所图示的顺序来执行动作。例如，可以修改该方法以在318确定是||r(i)||2≤2还是||x(i)||0≥T，并且可以基于该确定做出关于是否要执行另一迭代的决定。
另一OMP算法的实施例在下文的表1中阐述。
Figure BDA0001193823700000132
表1：示例OMP算法
在实施例中，异常的检测基于单次搏动的分析。从ECG记录中提取搏动。所提取的心跳可以被表示为
Figure BDA0001193823700000133
在实施例中，提供了在正常感测条件下获取的正常搏动的训练集合。为了简单起见，将首先考虑采用监督训练的实施例，其中，属于正常类别的标记训练数据是可用的(当提供正常和异常样本时，可能发生更容易的情形)。然后将讨论采用无监督训练的实施例，其中不提供训练数据的标签。
假设正常测试搏动s可以很好地近似为：
s≈Dx (6)
其中，D是未知字典，并且系数向量x是稀疏的。以下解释如何与所提出的异常检测器一起计算实施例的字典D和系数向量x的示例，即，字典学习和异常检测。
在实施例中，给出属于正常类别的训练样本的数据库
Figure BDA0001193823700000141
假设训练信号已经被生成为si＝Dxi+ei，其中，
Figure BDA0001193823700000142
是固定的但未知的字典，
Figure BDA0001193823700000143
是T-sparse(||x||0≤T)，，其中i＝1，...，n，，其中T是固定的已知值，并且
Figure BDA0001193823700000149
是具有有界功率的随机扰动噪声||ei||2≤ε，其中i＝1，…，n。确定给出给定训练样本的最佳T稀疏表示
Figure BDA0001193823700000144
和稀疏表示
Figure BDA0001193823700000145
的字典D。通过逐列级联向量si和xi获得的矩阵可以分别被表示为
Figure BDA0001193823700000146
和
Figure BDA0001193823700000147
实施例的合适的最小化问题可以被表示为：
Figure BDA0001193823700000148
||dj||2＝1 for j＝1，...，N (7)
其中F表示Frobenius范数。该条件||xi||0≤T约束训练集合中的每个向量中的表示的稀疏度，而||dj||2＝1是对字典D的每个列的范数的约束，其避免了平凡解。注意，该公式是等式(5)的矩阵版本，其中，D是未知的。这种联合优化问题可以使用K-SVD算法来求解。参见Aharon等人的Svd：An algorithm for designing overcomplete dictionaries forsparse representation，Signal Processing，IEEE Transactions on，第54卷，第11期，第4311至4322页，2006。图4中示出了学习稀疏模型的说明性示例。
K-SVD是通过在D和X上交替最小化来解决问题(诸如等式7的问题)的迭代算法。K-SVD开始于固定数目的迭代和初始字典D(0)(诸如具有完全随机项的字典、或具有训练向量的随机选择的子集的字典)。一般迭代在优化系数X和字典D之间交替。例如使用OMP算法首先计算系数，然后与系数的更新一起执行字典更新。联合字典和系数更新包括对字典中每个原子的迭代，联合更新该原子和对应的系数元素。对于每个字典原子dl，其中l＝1，...，N，残差被定义为：
Figure BDA0001193823700000151
其中，
Figure BDA0001193823700000152
是X的行l。残差评估字典性能，而不考虑所选原子的贡献。然后可以计算dl和
Figure BDA0001193823700000153
的最佳值：
Figure BDA0001193823700000154
由于dxT具有秩1，所以问题的最小元(例如，等式9的问题)可以被找到作为使用SVD的El的最优秩-1近似。给定E-L＝∪∑VT，El的SVD，d和x可以被设定为d＝u1和x＝σ1v1。注意，该更新通常不会给出稀疏的x，因此X的列通常不会是稀疏的。这可以通过将
Figure BDA0001193823700000155
定义为El对其中
Figure BDA0001193823700000156
中的对应的项非零的那些列并且应用与
Figure BDA0001193823700000157
相同的过程的限制来解决。
在图5中图示了应用K-SVD算法来定义字典D的方法500的实施例。方法500在502开始，并且进行到504。在504处，例如通过检索、接收或生成字典D、训练集合S、噪声界限ε、目标稀疏度T和/或迭代次数J来定义初始字典D(0)、训练集S＝[s1，…，sn]、目标稀疏度T、迭代次数J和噪声界限ε。该方法从504进行到506。
在506处，初始化OMP算法的索引i，并且该方法进行到508。在508处，方法500使用例如OMP算法来确定矩阵X的向量的系数。如所图示的，对于矩阵X中的第i个向量计算
Figure BDA0001193823700000158
并且方法从508进行到510，其中索引i递增。方法500进行到512，其中确定是否已经处理了矩阵X的所有向量，如所图示的，这通过确定是否i＝n+1来进行。当在512处确定并非所有的向量矩阵X已经被处理，则该方法从512进行到508以处理下一个向量。当在512处确定已经处理了矩阵X的所有向量时，该方法从512进行到514。
在514处，初始化用于优化字典D和系数X的索引k，并且方法500进行到516。在516处，方法500优化字典D和系数X。如所图示的，这包括设定l等于其表示使用
Figure BDA0001193823700000161
的S中的信号的索引；设定
Figure BDA0001193823700000162
设定
Figure BDA0001193823700000163
计算
Figure BDA0001193823700000164
并且更新
Figure BDA0001193823700000165
和
Figure BDA0001193823700000166
方法500进行到518，其中，索引k递增。该方法进行到520，其中确定是否已经处理了字典D的所有向量，如所图示的，这通过确定是否k＝n+1来进行。当在520处确定不是字典D和矩阵X的所有的向量已经被处理，则该方法从520进行到516以处理下一个向量。当在520处确定已经处理了字典D和矩阵X的所有向量时，方法从520进行到522。
在522处，方法500递增索引j并且进行到524。在524处，方法500确定是否发生了期望数量的迭代，如所图示的，这通过确定是否j＝J+1来进行。当在524处确定需要更多迭代，则方法从524进行到506以开始下一迭代。当在524处确定已经完成了期望数量的迭代时，该方法从524进行到526，其中输出字典D(J)和系数矩阵X。方法500的实施例可以执行附加动作、可以执行较少动作、可以修改动作、并且可以以不同于所图示的顺序来执行动作。
另一K-SVD算法的实施例在下文的表2中阐述。
Figure BDA0001193823700000171
表2：示例K-SVD算法
观察到优化问题(诸如上文的等式9的问题)也可以是矩阵因式分解的术语。发现潜在字典的问题可以被视为与发现矩阵S的因式分解为DX的问题相同，其中，D和X具有指示的大小，并且X具有稀疏列。这个矩阵因式分解观点将这个问题与非负矩阵因式分解的相关问题(特别是稀疏非负矩阵因式分解)相关联。
现在假设给出测试信号
Figure BDA0001193823700000172
相对于学习字典D的s的稀疏编码与计算预期为T稀疏的系数的向量x相对应。该问题可以被公式化为
Figure BDA0001193823700000173
问题的实例，并且可以例如通过使用图3的方法300的实施例或表1的算法的实施例来计算。
测试信号s的T稀疏近似x可以用作测量信号s是正常或异常的程度的工具。为了定量评估测试信号与正常信号的接近程度，可以使用测量稀疏近似精确程度的异常指示器。当通过D的几个原子的线性组合来近似异常信号时，可以预期重建误差的实质性偏差，因为D的原子被学习为表示正常数据。为此，在实施例中可行的异常指示器将是均方根误差(RMSE)：
Figure BDA0001193823700000181
该异常指示器可以被认为从
Figure BDA0001193823700000182
(信号的空间)映射到低维特征空间(在这种情况下为
Figure BDA0001193823700000183
)，其中，来自正常信号的特征遵循未知分布。产生该指示器的异常值的信号可以被认为是异常。检测异常的最直接的解决方案之一包括：构建合适的置信区域并且考虑到产生落在该区域外的RMSE(D.T.)(s)的信号s的异常。在实施例中，标量异常指示器可以是：
Iα＝[0，qα] (11)
其中，qα表示RMSE(9)的经验分布的α-分位数。特别地，如果
Figure BDA0001193823700000184
或等效RMSE(D,T)(s)>qα，则s可以被认为是异常的。
在实施例中，可以考虑四分位范围
Figure BDA0001193823700000185
并且将置信区域定义为：
Figure BDA0001193823700000186
在这种情况下，如果
Figure BDA0001193823700000187
则可以认为s是异常的。图6中描绘了示例流水线方法600的图示。在602处，获得或提供包括心跳信息的测试信号s。在604处对信号进行OMP处理(参见例如，图3的方法300)以产生均方根误差RMSE的指示。在606处，基于所定义的置信区域(例如，在上文所讨论的示例中，γ表示qα或q3/4+plQR)来确定阈值γ。在608处比较RMSE误差和阈值γ以确定RMSE误差是否大于阈值γ。当在608处确定RMSE误差大于阈值时，方法600进行到610以执行与检测到的异常相一致的处理(诸如在实施例中发起指示异常的数据的传输)，继续处理下一搏动等。当在608处没有确定RMSE误差大于阈值时，该方法进行到612以执行与正常心跳相一致的处理(诸如在实施例中不发起数据传输)，并且处理下一搏动等。在实施例中可以采用重建误差的其它指示。
如上文所提及的，用户状态的特征可以在于心率HR和所测量的加速度的能量EE。图7以图表形式表示可以从中测量心率的心跳702和可以从中确定加速度的能量的对应的加速度数据704。在日常活动期间，可以预期用户的心率和搏动的形态都响应于用户状态的改变而改变。实施例可以自动调整以补偿改变的用户状态。为了便于对这种变化的适应性，当前用户的状态的特征可以在于心率(HR)和对应的所测量的加速度(EE)的能量。心率可以以单位时间内的搏动来测量，诸如心脏每分钟的收缩数(bpm)。可以通过考虑每个例如10秒的固定大小的窗口来计算当前心率，并且对于每个窗口，计算窗口中的搏动的准时心率的鲁棒统计(在本文其它地方更详细地讨论)。加速度的能量可以被定义为与当前搏动相对应的加速度计数据的l2范数。
如图8所描绘的，组合这两个特征便于覆盖几个状态(例如，休息、行走、睡觉、努力后休息、努力后行走)，以及相对于组合来标识外点搏动。观察到，而不是被视为对用户状态进行分类的练习，该过程可以被视为勾画特征空间(HR；EE)。例如，每个配对值(HR；EE)可以与要应用于搏动s和字典D的对应的变换配对P(HR；EE)和F(HR；EE)相关联，以便适应到对应的状态。例如sA＝P(HR；EE)(s)和DA＝P(HR；EE)·F(HR；EE)(D)，其中，sA和DA分别表示自适应的搏动和字典。为了简单起见，复合P(HR；EE)·F(HR；EE)可以表示为G(HR，EE)。P(HR；EE)可以被认为是搏动变换，并且G(HR，EE)可以被认为是字典变换。
在实施例中，变换P和F可以是通用的，例如，它们不依赖于特定用户。例如，可以通过使用几个用户作为训练集合来学习变换。变换可以用于在用户状态改变时保持重建误差(分别在自适应之前和之后)基本不变。这样，可以使检测器相对于用户状态的改变而基本上不变。这样的变换可以经由例如统计分析(诸如Procrustes分析、迁移学习、动态时间扭曲等)来离线地生成。
图8以图表形式表示示例用户状态(如所图示的，休息、行走、睡觉、努力后休息、努力后行走)、以及用户状态与心率HR和所测量的加速度的能量EE之间的示例关系。如上文所讨论的，可以定义和采用用于考虑用户状态的变换。图9图示了要应用于所测量的搏动和字典以考虑用户状态的改变的示例变换配对P(HR，EE)902和F(HR，EE)904。示例变换及其应用在本文其它地方进行更详细地讨论。可以定义和采用其它用户状态和关系。
图10图示了用于检测异常心跳的系统1000的实施例。如所图示的系统1000包括一个或多个处理器P、一个或多个存储器M和分立电路DC，其单独或以各种组合被配置成提供功能性和/或实现系统1000的电路。如所图示的，系统1000通信地耦合(例如，有线或无线耦合)到一个或多个传感器1002，该传感器1002向系统1000提供诸如ECG迹线和加速度计数据之类的数据。在一些实施例中，一个或多个传感器可以集成到系统1000。系统1000包括训练电路或模块1020，其在操作中使用例如本文中所讨论的方法之一的实施例来训练字典。系统1000包括检测电路或模块1030，其在操作中使用训练的字典和本文中所讨论的方法的实施例来检测异常心跳。如所图示的，系统1000还通信地耦合(例如，有线或无线耦合)到提供患者监测功能性的一个或多个远程或本地服务器1050，诸如本地主机(例如，蜂窝电话、移动设备等)和本地或远程服务器。在一些实施例中，系统1000可以通信地耦合到本地主机(例如，蜂窝电话)，本地主机又通信地耦合到远程服务器(例如，健康监测服务器)。系统1000包括通信电路或模块1040，其在操作中控制系统1000和其它设备(诸如传感器1002和远程/本地服务器1050)之间的通信。
图11图示了例如可以由图10的系统1000的实施例采用的检测异常心跳的方法的实施例。为了方便起见，将参考图10的系统1000来讨论图11。方法1100在1102处开始并且进行到1104。在1104处，例如通过诸如图10的训练模块1020之类的训练模块使用数据(其可以是实时数据或所存储的数据，诸如ECG迹线和可选的加速度计数据)以训练稀疏模型字典D，其可以存储在例如图10的系统1000的存储器M中。数据可以是例如从诸如图10的传感器1002之类的传感器接收的实况数据，或所存储的训练数据，诸如从诸如图10的服务器1050之类的服务器接收的数据。输入可包括例如7分钟或更长的ECG迹线。七分钟的迹线通常将提供足够的搏动(例如，500或更多搏动)以生成足够的训练数据集合。在实施例中，训练数据集合可以是在训练时段期间处于相对休息状态的用户的数据集合。
方法1100从1104进行到1106。在1106处，例如通过诸如图10的系统1000的检测模块1030之类的检测模块，使用在1104处训练的稀疏模型字典，连同加速度计数据一起，以将连续的ECG迹线的搏动分类为正常或异常。方法1100从1106进行到1108，其中，如所图示的，该方法终止。方法1100可以包括图11中未示出的附加动作，可以不包括图11中所示的所有动作，并且可以将图11中所示的动作分成多个动作。例如，实施例可以包括响应于检测到异常心跳而启用或唤醒通信电路(诸如图10的通信电路1040)的附加动作。这可以通过在期望通信(诸如检测到异常心跳的指示)时启用通信电路来便于节省系统功率。
图12图示了训练用于检测异常心跳的字典的训练子例程或方法1200的实施例，其可以例如被图11的方法1100的实施例和/或被图10的系统1000的实施例来采用。训练是基于对单独搏动的分析，例如，用户的ECG迹线中的单独心跳的分析。训练阶段选择适当的样本搏动训练集合，并且基于样本搏动训练集合来生成稀疏字典D。图13是可以输入到训练子例程或方法1200的ECG迹线或流的一部分的图形表示。该阶段的输入可以是例如约七分钟的时间长度的ECG流，其可以提供约500搏动和加速度计数据，其中用户处于相对休息状态。输出是与固定稀疏度级别T相关联的字典D。
尽管心动周期在开始时以P波开始，但是实施例不需要精确地隔离心跳。相反，可以选择期望沿着时间重复并且彼此一致的追踪部分。在实施例中，R峰用作重复现象的参考点，并且在R峰周围选择固定窗口。在讨论中，心跳或简单的搏动是指窗口内部的样本的这样的集合。还有，实施例可以采用预处理和滤波步骤来分析信号。
方法1200从1202进行到1204，例如，响应于调用训练子例程。
在1204处，检测每个心跳的R峰值。这可以例如通过使用两个中值滤波器对ECG迹线进行滤波以去除基线、并且使用例如Pan-Tompkins算法检测滤波信号中的R峰来完成。参见Pan等人的A real-time qrs detection algorithm，Biomedical Engineering，IEEETransaction on，第3期，第230至236页，1985。R峰用作心跳的参考点。方法1200从1204进行到1206。
在1206处，估计心率。例如，可以估计准时心率(RR)和参考心率(RHR)。心率是心脏每分钟的收缩的数目，例如，每分钟的搏动。可以例如通过将搏动的R峰和相邻搏动的R峰之间的以秒计的时间差除以60来测量准时心率。参考心率可以被认为是休息条件下特定用户的参考心率，并且可以例如基于迹线的准时心率的统计分析来确定。例如，RHR可以被设定为准时心率的中值。实施例可以标识准时心率中的内点，并且排除外点。例如，可以通过确定单独准时心率是否在迹线的参考心率的阈值范围内(例如，[RHR-5bpm，RHR+5bpm])来标识内点。该方法从1206进行到1208。
在1208处，例如使用参考窗口大小的参考窗口对流的搏动进行窗口化，该参考窗口大小可以是固定大小，诸如例如，R-0.35秒、R+0.35秒。对于每个R峰，从ECG流中提取具有位于R峰中的中心的固定大小的窗口。例如，可以考虑大小为0：7s(围绕R峰值的0:35s)的窗口。观察到，由于RHR预期在60bpm至100bpm的范围内，0：7s的窗口便于选择当前搏动的所有PQRST波而不包括来自先前或后续搏动的波。所选择的固定大小的窗口大小可以被称为参考窗口大小(RWS)，并且被视为用于休息条件下的特定用户的参考窗口。在该水平下，心跳或简单的搏动是指在所选择的窗口内部的样本的收集。图15是ECG迹线的加窗搏动的示例的图形表示。该方法从1208进行到1210。
在1210处，执行运动估计。在实施例中，提供加速度计数据作为输入。ECG数据流和加速度计数据被同步或假定为同步。对于每个心跳，测量相关联的加速度计数据的加速度的能量(EE)。在实施例中，加速度的能量被定义为与围绕R峰的大小为RWS的窗口相对应的加速度计数据的l2范数。可以沿着迹线基于滑动窗口使用例如置信区间来标识加速度计数据分布中的外点。参见图7至图8及其描述。图14是加速度计数据的示例的图形表示。该方法从1210进行到1212。
在1212处，选择搏动以获得具有几乎恒定心率的可靠训练数据集合。这可以通过丢弃在1206(RR外点)或1210(与外点加速度计数据相关联的搏动、和/或与休息状态不一致的搏动)处被标识为外点的搏动来完成。关于相对于准时心率RR的外点搏动，可以排除这种搏动，因为外点状态可以是搏动被噪声或运动伪像高度影响的指示器，或者由于例如用户的不随意移动而导致的感测过程中的故障。关于相对于加速度计数据的外点搏动，这种搏动可以类似地被排除为可能受到噪声或运动伪像的影响。关于与休息状态不一致的搏动，可以从用于训练集合的数据中排除这种搏动、和/或大量的这种搏动可以指示该迹线不适合在训练阶段中使用。在实施例中，可以为训练阶段选择具有内点准时心率和内点加速度计数据的搏动。所选择的搏动还可以在1212处经历预处理，诸如通过减去平均值来获得零均值向量(例如，垂直地平移搏动)。经预处理的选择的搏动形成训练数据集si，其也可以表示为：
Figure BDA0001193823700000231
如果训练数据集中所选择的搏动的数目太低(例如，小于阈值数目，诸如500次搏动)，则可以分析附加搏动并且将其添加到训练集合。该方法从1212进行到1214。
在1214处，从训练数据集合si中学习字典。这可以例如使用图5的方法500的实施例和/或表2的K-SVD算法的实施例来完成。该方法从1214进行到1216，其中子例程可以返回字典D，可能发生进一步的处理等。例如，字典D可以被提供给异常检测子例程，如本文中别处所讨论的。方法1200的实施例可以执行图12中未示出的附加动作，可以以不同顺序执行图12的动作，可以省略图12中所示的动作，并且可以组合或分割图12中所示的动作。例如，在实施例中，方法1200可以被修改以检查1212处的搏动选择是否指示ECG迹线不适合用于生成字典D，并且当确定ECG迹线不适合时，收集另一ECG迹线以用于训练字典D。
图16图示了使用例如训练的字典、当前ECG迹线和加速度计数据来检测异常心跳的子例程或方法1600的实施例。图17是可以作为子例程1600的输入提供的ECG迹线的图形表示。方法1600可以应用于单独搏动，并且可以顺序地应用于一系列单独搏动，或者并行地应用于单独搏动集合。每个搏动的输出是搏动是正常还是异常的指示。子例程1600从输入中的连续ECG流中提取单个心跳。还有，在分析信号之前可以采用一些预处理和滤波步骤。方法1600在1602开始并且进行到1604。
在1604，检测每个心跳的R峰。这可以例如通过对ECG迹线进行滤波来完成，例如，使用两个中值滤波器来去除基线，并且使用例如Pan-Tompkins算法来检测经滤波的信号中的R峰。参见Pan等人的A real-time qrs detection algorithm，BiomedicalEngineering，IEEE Transaction on，第3期，第230至236页，1985。图18是从ECG迹线提取的搏动的示例的图形说明。方法1600从1604进行到1606。
在1606处，例如通过计算搏动和相邻搏动的R峰之间的时间差来测量准时心率(RR)。为参考窗口(时间间隔)确定局部心率(HR)，其可以是例如固定大小，诸如10秒。可以采用统计分析来确定局部心率HR，例如，在实施例中，可以将窗口中的搏动的准时心率的中值用作局部心率。在实施例中，也可以确定四分位范围IQR＝Q3-Q1，其中，Q3和Q1分别是窗口中的RR的第三和第一四分位数。可以例如基于阈值来标识准时心率中的外点。例如，低于Q1-1.5(IQR)或高于Q3+1.5(IQR)的RR可以被标识为外点。可以采用标识外点的其它方法，诸如中值绝对偏差(MAD)。应当指出，在实施例中可以采用不同的阈值和其它准则来标识训练和检测方法或例程中的外点。在异常检测之前可以去除被标识为外点的RR，如下文所描述的。
该方法从1606进行到1608。在1608处，例如使用参考窗口大小的参考窗口对流的搏动进行窗口化，参考窗口大小可以是固定大小，例如R-0.35秒，R+0.35秒。对于每个R峰，从ECG流中提取具有位于峰中的中心的固定大小的窗口。例如，可以考虑大小为0：7s(围绕R峰值的0:35s)的窗口。可以选择其它窗口大小，例如以便于选择当前搏动的所有PQRST波，而不包括来自先前或后续搏动的波。所选择的固定大小的窗口大小可以被称为参考窗口大小(RWS)，并且被认为是特定用户的参考窗口。在该水平下，心跳或简单的搏动是指在所选择的窗口内的样本的收集。图15是ECG迹线的加窗搏动的示例的图形表示。
该方法从1608进行到1610。在1610处，执行运动估计。在实施例中，提供加速度计数据作为输入。ECG数据流和加速度计数据被同步或假设为同步。对于每个心跳，测量相关联的加速度计数据的加速度的能量(EE)。在实施例中，加速度的能量被定义为与围绕R峰的大小为RWS的窗口相对应的加速度计数据的l2范数。可以沿着迹线基于滑动窗口使用例如置信区间来标识加速度计数据分布中的外点。参见图7至图8及其描述。应当指出，在实施例中可以采用不同的阈值和其它准则来标识训练和检测方法或子例程中的外点。图14是加速度计数据的示例的图形表示。该方法从1610进行到1612。
在1612处，选择搏动以便排除可能引起异常的假阳性指示的搏动。在实施例中，这可以通过排除具有以下两项的搏动来完成：(i)相对于局部心率被标识为外点的准时心率和(ii)被标识为外点的对应的加速度计数据。可以排除这样的搏动，因为在两个准则下的外点状态可以是搏动被噪声或运动伪像高度影响的指示器，或者由于例如用户的不随意移动而导致的感测过程中的故障。应当指出，作为相对于RR分布的外点的、但是作为相对于加速度计数据的内点的搏动可以是异常搏动(例如，心房早搏、节律性异常搏动等)。因此，在实施例中，作为相对于条件(i)的外点的、但是作为相对于条件(ii)的内点的搏动可以包括在针对异常检测而选择的搏动中。在实施例中，作为相对于条件(i)的外点的、但是作为相对于条件(ii)作为内点的搏动可以在1612处被分类为异常，其中，搏动应当被考虑为异常的指示在1622处考虑，如图12中虚线1626所示。
还应当指出，在实施例中，因为相对于条件(i)的内点状态可以是保留搏动的形态的指示器，所以可以选择作为相对于条件(i)的内点但是作为相对于条件(ii)的外点的搏动，以用于异常检测，因此可以通过在分析中包括搏动来获得有用的信息。所选择的搏动还可以在1612处经历预处理，诸如通过减去平均值(例如，垂直地平移搏动)来获得零均值向量。该方法从1612进行到1614。
在1614处，每次搏动与一个配对的变换P(HR，EE)和F(HR，EE)相关联。参见例如图7至图9及其描述。该方法从1614进行到1616。在1616处，将搏动自适应应用于搏动以考虑当前用户状态。在实施例中，使用所标识的变换P来自适应当前搏动s，如下：
sA:＝P(HR,EE)·s (13)
其中，sA是自适应的搏动。例如，如果对于当前搏动s，HR相对较高，而EE相对较低(相对于一些导出阈值)，则这可以被视为当前用户状态在努力之后休息的指示。参见图8。
在这种情况下，变换P(HR，EE)可以简单地与尺寸小于窗口RWS的窗口上的投影相重合。例如，如果心率较高，则搏动可以彼此跟随得更接近和更宽，窗口(例如，大小为RWS的窗口)可以包括来自先前或后续搏动的波。参见图19，其图示了与在努力条件(底部)之后的休息中的示例搏动形成对比的休息条件(顶部)下的示例搏动。在数学术语中，P(HR，EE)可以与单位矩阵相对应，其中对角线的第一元素和最后元素中的一个或多个元素被设定为零。这种情形在图20中图示，图20是使用变换P(HR，EE)来自适应努力之后休息状态下的用户的搏动sRWS以生成自适应的搏动sA的实施例中的搏动自适应的示例应用的图形说明。方法1600从1616进行到1618。
在1618，将字典自适应应用于搏动以考虑当前用户状态。在实施例中，使用所标识的变换F来自适应电流D，如下：
DA＝P(HR,EE)·F(HR,EE)(D)＝G(HR,EE)(D) (14)
其中，DA是自适应的字典。例如，如果对于当前搏动s，HR是相对较高，而EE是相对较低(相对于一些导出阈值)(与在1616处的示例中讨论的相同的条件)，这可以被看作是当前用户状态正在在努力之后休息。参见图8。在这种情况下，如上文所讨论的，变换P(HR，EE)可以简单地与尺寸小于窗口RWS的窗口上的投影相重合。从图19可以看出，图19图示了与在努力条件(底部)之后的休息中的示例搏动形成对比的休息条件(顶部)下的示例搏动，可以看出，在努力状态之后的休息中，搏动PQRST的形态略有伸展。因此，在实施例中，F(HR，EE)可以是包含与拉伸波相对应的权重的带矩阵。参见图9，其示出了当搏动处于努力状态之后的休息中时可以由实施例采用的矩阵P 902和矩阵F 904的示例。注意，例如，当搏动与努力之后休息用户状态相对应时，G(HR，EE)＝P(HR，EE)·F(HR，EE)可以是带矩阵，其还将字典的原子投影到大小比RWS更小的窗口上，如图21所示。方法1600从1618进行到1620。
在1620处，相对于自适应的字典DA来确定自适应的测试搏动sA的T稀疏近似x。注意，sA∈RMA，s相对于自适应的字典DA的自适应测试心跳。这可以使用OMP算法(诸如图3的方法300或表1的算法)来完成，以提供对以下问题的近似解：
Figure BDA0001193823700000271
图22是在用于在努力状态之后的休息中的用户的实施例中的OMP算法的迭代中的向量sA和x与字典DA之间的关系的图形说明。
方法1600从1620进行到1622，其中，执行异常搏动的检测。这可以通过以下各项来完成：基于自适应的测试搏动sA相对于自适应的字典DA的稀疏近似x来确定搏动s的异常指示器，将异常指示器与异常指示器的预期值或值的范围进行比较，以及基于比较来将搏动分类为正常或异常。示例异常指示器是均方根误差：
Figure BDA0001193823700000272
其中，sA是自适应的测试搏动向量，DA是自适应的字典，x是自适应的测试搏动向量的稀疏近似。可以将均方根误差与均方根误差的期望值进行比较，并且可以基于比较将搏动分类为正常或异常。可以采用一个或多个阈值、一个或多个阈值范围和一个或多个查找表来将搏动分类为正常或异常。例如，如果
Figure BDA0001193823700000281
则搏动可以被认为是异常的，其中，γ是确定对异常的响应性的阈值，并且通常可以被设定为控制假阳性率，例如，被分类为异常的正常搏动的百分比。参见异常指示器的讨论和上文的图6。如上文关于1612所讨论的，作为相对于条件(i)的外点的但是作为相对于条件(ii)的内点的搏动可以在1612处被分类为异常，其中搏动应当被分类为异常的指示在1622处被考虑，如虚线1626所示。
方法1600从1622进行到1624，其中，子例程可以返回搏动是否是异常的指示，其中，可以发生进一步的处理等。方法1600的实施例可以执行图16中未示出的附加动作，可以以不同顺序执行图16的动作，可以省略图16中所示的动作，可以修改图16中所示的动作，并且可以组合或分割图16中所示的动作。例如，在实施例中，可以修改方法1600以检查附加搏动是否要在1622处分类，并且当确定附加的搏动仍然要被分类时，返回到1604以处理下一搏动。在另一示例中，可以修改子例程以执行误差处理，诸如响应于搏动被分类为异常，从而发起与被分类为异常的搏动有关的数据的传输。在另一示例中，例程1600可以被修改以追踪在阈值时间段内从在1612处选择的搏动中排除的搏动的比例，并且当该比例指示字典可能不再有效时(例如，由于传感器的改变位置)，调用子例程来更新字典，诸如图12的子例程1200。可以采用阈值来确定字典是否可能不再有效。
通常，可穿戴设备连续地或周期性地将所有数据发送到主机设备，该主机设备可以进而将所有数据发送到远程应用服务器。这可能导致可穿戴设备和主机设备的强烈能量消耗。在实施例的一种操作模式中，训练阶段(例如，定义字典)可以在可穿戴设备上执行(参见图10的系统1000)。当检测到异常时，感测到的数据可以由可穿戴设备传输到主机(例如，蜂窝电话，参见图10的本地服务器1050)，而不是周期性地或连续地。主机可以依次将数据发送到远程服务器(参见图10的远程服务器1050)。在实施例的一种操作模式中，训练数据可以由可穿戴设备传输到主机设备。字典可以由主机设备或远程服务器生成，并且被传输到可穿戴设备。当检测到异常时，感测到的数据可以由可穿戴设备传输到主机，而不是周期性地或连续地。
实施例可以提供相对于不同的用户状态(时间和移动变化的动态)是自适应的，易于可重新配置(重新训练)的用户特定字典(模型学习表征单个用户的心跳的形状)，而不需要学习与每个不同用户状态相关联的搏动形状，并且其具有鲁棒的训练和测试(检测)阶段(例如，可以检测和移除噪声和运动伪像)及其各种组合。实施例可以在无监督训练模式下操作(可以从训练数据中滤波杂散信号)。在实施例中，训练和测试(检测)的全部或部分可以在可穿戴设备上执行，而非在主机或远程服务器上执行。在实施例中可以采用压缩数据(诸如使用压缩的感测数据)。
在实施例中，使用测试数据集合在脉冲传感器上实施测试(检测)阶段的执行，并且测试(检测)阶段的执行足够快以便于在可穿戴设备上的实时执行。下文对实验中使用的测试设置进行描述。使用MIT-BIH心律失常数据库来测试该程序。参见Goldberger等人的PhysioBank，Physio Toolkit and PhysioNet：Components of a new research resourcefor complex physiologic signals，Circulation，第101卷，第23期，第e215至e220页，2000。
由于对于MIT-BIH心律失常数据库，每个患者的可用搏动数目不是非常高(在大多数情况下从1000至2500)，通常执行k-交叉验证。在k折交叉验证中，将原始样本划分为k个相等大小的子样本。在k个子样本中，保留单个子样本作为用于测试模型的验证数据，并且剩余的k-1个子样本用作训练数据。交叉验证过程然后重复k次(折)，其中k个子样本中的每一个子样本精确地使用一次作为验证数据。然后对来自折的k个结果进行平均以产生单个估计。
异常检测性能通常使用常规接收器操作特点(ROC)曲线来评估，其中，参数是阈值γ。ROC曲线绘制了真阳性检测率(TPR)(被定义为真阳性结果(正确检测的异常)的比例)与对应的假阳性率误差(FPR)(被定义为假阳性结果(检测为异常的正常信号)的比例)。当异常检测阈值γ变化时，这两个值相对于彼此而改变。回想一下，曲线下面积(通常被称为AUC)等于检测器将排列高于随机选择的正常搏动的随机选择的异常的概率。
对于γ参数的固定值，还考虑以下因素：精确度(被定义为真实结果(真阳性和真阴性，例如，正确检测的异常和正确检测的正常搏动)在所检查的情况的总数中的比例)、精度(被定义为真阳性结果(正确检测的异常)对所有阳性结果(真阳性和假阳性，例如，正确检测的异常和被检测为异常的正常搏动)的比例)、以及假阴性率(FNR)(被定义为假阴性结果的比例(错误检测的异常))。
其它分析的参数包括稀疏度级别T、搏动M的长度和字典冗余度r：＝N/M。
除非另有说明，否则本文中所讨论的程序一次处理来自一个单个导联的信号。
包括在MIT-BIH心律失常数据库中的ECG的源是在1975年至1979年期间由BethIsrael Hospital Arrhythmia Laboratory获得的超过4000个长期Holter记录的集合。大约60％的记录是从住院病人获得的。数据库包含从该集合中随机选择的23个记录(编号从100到124，包括一些缺失的编号)，和从相同集合中选择的25个记录(编号从200到234，同样包括一些缺失的编号)，以包括各种罕见但临床上重要的现象，该现象不会很好地由Holter记录的小随机样本来表示。48个记录中的每个记录略长30多分钟。第一组旨在用作心律失常检测器在常规临床使用中可能遇到的各种波形和伪像的代表性样本。选择第二组中的记录以包括复杂的室性、交接性和室上性心律失常和传导异常。受试者为25名年龄32至89岁的男性，22名年龄23至89岁的女性。记录201和202来自同一个男性受试者。
对于每个患者，从两个导联来记录两个信号并且结合对搏动的注释。如在临床实践中预期的那样，受试者之间的导联不同。在大多数记录中，一个通道是通过将电极放置在胸部上作为动态ECG记录的标准实践获得的修改的肢体导联II(MLII)，并且另一个通道通常是V1(有时是V2，V4或V5，取决于受试者)。
采样频率为360Hz。选择这个值以适应使用简单的数字陷波滤波器以去除60Hz(电源频率)干扰。
对于每个记录，两个迹线给予独立工作的两个心脏病专家。心脏病专家添加了检测器丢失搏动的位置的附加的搏动标记，必要时删除假检测，并且改变所有异常搏动的标记。他们还添加了节奏标记、信号质量标记和评论。注释从纸图记录中转录。一旦给定记录的两组心脏病专家的注释都被转录和验证，它们就被自动地逐一搏动进行比较。每个偏差都通过协商一致进行审查和解决。一旦注释匹配，结果通常就附加在R波峰处。
在实验中，给定来自MIT-BIH心律失常数据库的ECG迹线，根据相关注释来提取和标记单独搏动。这里，简要描述实验程序。从所选择的导联中读取样本以及相关注释文件(注释与相关联的搏动的R峰的时间一起存储)。使用Pan-Tompkins's算法(参见Pan等人的Areal-time qrs detection algorithm，Biomedical Engineering，IEEE Transaction on，第3期，第230至236页，1985)检测R-峰，并且时间上最接近的标记被分配给每个峰。提取具有位于每个R峰中的中心的固定大小的窗口。选择分别与长度M＝360、M＝252或144的搏动相对应的1s、0.7s或0.4s的窗口大小。搏动被垂直平移以便具有零均值。
属于不同患者的正常心跳常常表现出显著差异。特别地，一些患者具有非常规则的迹线，而其他患者呈现几个不规则性。为了突出一些特定行为并且易于对结果的理解，在一些实验中，选择患者编号119，213和105作为样本。选择这样的患者是因为其正常搏动的形态显示不同的行为：来自患者119的正常搏动总是具有相同的规则形状，而对于患者213和105，形态变得越来越不规则。异常搏动通常从正常搏动中在视觉上可标识(如在三个所选示例中发生的)。由患者222表示特别不幸运的情况，其中不仅正常搏动显著混乱，而且异常搏动难以或不能与正常搏动区分开。
图23A至图23D示出了来自患者119，213，105和222的示例搏动。在图23A至图23D中，左侧的图像以黑色表示正常搏动，右侧的图像以黑色表示正常搏动，并且异常搏动以其它颜色表示。
如图24和图25所演示的，初步结果是令人鼓舞的。特别地，在图24中，描绘了异常检测器的全局行为，因为箱线图呈现了MIT-BIH数据库中所有患者的结果。图25替代为表示样本患者的结果。对于患者119而言获得最佳结果，对于患者105而言获得最差结果，这与预期结果一致。其他患者的结果可以看作是这两种极端情况之间的范围。
图24是MIT-BIH数据库中所有患者的ROC曲线下面积(AUC)的箱线图。箱线图的底部和顶部分别是AUC的第一四分位数和第三四分位数，并且箱内部的红色带是中值AUC值。晶须(whisker)的末端表示仍然在下四分位数的1.5IQR内的最低AUC，并且最高AUC仍然在上四分位数的1.5IQR内。不包括在晶须之间的数据被绘制为具有红色标记的外点。箱线图示出了稀疏度水平(spar)和冗余参数(红色)的不同值。在左侧列上，spar＝1并且red＝10，在中心列上，spar＝4并且red＝10，以及在右侧列上，spar＝4并且red＝25。窗口大小为0.7s。
图25图示了样本患者的检测性能。从顶部到底部，示出了来自患者119，213，105的结果。对于每个患者，在左侧列上示出了ROC曲线。在中心列上，患者数目(dataset)，分别用于训练(train)、评估(thres)和测试(test)阶段的负示例(#neg)(例如，正常波动)的数目，在测试阶段使用的正示例(#pos)(例如，异常搏动)的数目，稀疏度水平(sparsity)和字典大小(dict.size)；所分析的特征是：通过设定与95百分位相对应的参数λ获得的ROC曲线下面积(AUC)、精确度、精度、预期假阴性率(eFNR)和预期真阳性率(eTPR＝1-eFNR)，参见上文对特征的解释的讨论。在右侧列上图示了后验概率。
一些实施例可以采用或包括计算机程序产品的形式。例如，根据一个实施例，提供了一种包括适于执行上文所描述的方法或功能中的一个或多个的计算机程序的计算机可读介质。介质可以是诸如例如只读存储器(ROM)芯片的物理存储介质，或诸如数字通用盘(DVD-ROM)、压缩盘(CD-ROM)、硬盘、存储器、网络或便携式媒体物品之类的盘以由适当的驱动器或经由适当的连接来读取，包括如以一个或多个条形码编码或存储在一个或多个这样的计算机可读介质上的并且可由适当的阅读器设备读取的其它有关代码。
更进一步地，在一些实施例中，系统和/或模块和/或电路和/或框中的一些可以以其它方式实现或提供，诸如至少部分地以固件和/或硬件，包括但不限于，一个或多个专用集成电路(ASIC)、数字信号处理器、分立电路、逻辑门、标准集成电路、状态机、查找表、控制器(例如，通过执行适当的指令，并且包括微控制器和/嵌入式控制器)、现场可编程门阵列(FPGA)、复杂可编程逻辑器件(CPLD)等、以及采用RFID技术的设备、及其各种组合。
可以组合上文所描述的各种实施例以提供另外的实施例。如果必要，则可以修改实施例的各方面以采用各种专利、申请和出版物的概念，以提供另外的实施例。
根据上文所详述的描述，可以对这些实施例进行这些改变和其它改变。通常，在所附权利要求中，所使用的术语不应被解释为将权利要求限制于说明书和权利要求中公开的具体实施例，而是应被解释为包括所有可能的实施例以及这些权利要求所授权的等同权利要求的全部范围。因此，权利要求不受本公开限制。

Claims (24)

1.一种用于异常检测的方法，包括：
检测心跳信号中的搏动；
测量加速度的能量；
将心率和所测量的所述加速度的能量与检测到的搏动相关联；
基于与所述检测到的搏动相关联的所述心率和所测量的所述加速度的能量，来将所述检测到的搏动选择性地包括在检测到的搏动集合中；以及
使用稀疏近似模型的字典来检测所述检测到的搏动集合中的异常搏动，所述异常搏动的所述检测包括：
基于与所述检测到的搏动集合中的搏动相关联的所述心率和所述加速度的能量，将搏动变换应用于所选择的搏动集合中的所述搏动，生成自适应的搏动；
基于与所述检测到的搏动集合中的所述搏动相关联的所述心率和所述加速度的能量，将字典变换应用于所述字典，生成自适应的字典；
生成相对于所述自适应的字典的所述自适应的搏动的稀疏近似；和
基于所生成的所述稀疏近似来将所述检测到的搏动集合中的所述搏动分类为异常。
2.根据权利要求1所述的方法，其中，基于所述稀疏近似来将所述搏动集合中的所述搏动分类为异常包括：
确定与所述自适应的搏动的所述稀疏近似相关联的均方根误差；
将所述均方根误差与阈值误差进行比较；和
基于所述比较来将所述搏动集合中的所述搏动分类为异常或正常。
3.根据权利要求1所述的方法，包括：
响应于将所述搏动集合中的所述搏动分类为异常，发起与所述搏动集合中的所述搏动有关的数据的传输。
4.根据权利要求1所述的方法，其中所述心率是局部心率，并且将所述检测到的搏动选择性地包括在所述检测到的搏动集合中包括：
确定所述检测到的搏动的准时心率是否是相对于与所述检测到的搏动相关联的所述局部心率的外点心率；
确定所述检测到的搏动是否是相对于与所述检测到的搏动相关联的所述加速度的能量的外点搏动；
当所述准时心率是相对于所述局部心率的外点心率、并且所述检测到的搏动是相对于与所述检测到的搏动相关联的所述加速度的能量的外点搏动时，从所述检测到的搏动集合中排除所述检测到的搏动；和
当所述准时心率是相对于所述局部心率的内点心率、或者所述检测到的搏动是相对于与所述检测到的搏动相关联的所述加速度的能量的内点搏动时，将所述检测到的搏动包括在所述检测到的搏动集合中。
5.一种用于异常检测的方法，包括：
检测心跳信号中的搏动；
测量加速度的能量；
将心率和所测量的所述加速度的能量与检测到的搏动相关联，其中所述心率是局部心率；
基于与所述检测到的搏动相关联的所述心率和所测量的所述加速度的能量，来将所述检测到的搏动选择性地包括在检测到的搏动集合中，其中将所述检测到的搏动选择性地包括在所述检测到的搏动集合中包括：
确定所述检测到的搏动的准时心率是否是相对于与所述检测到的搏动相关联的所述局部心率的外点心率；
确定所述检测到的搏动是否是相对于与所述检测到的搏动相关联的所述加速度的能量的外点搏动；
当所述准时心率是相对于所述局部心率的外点心率、或者所述检测到的搏动是相对于与所述检测到的搏动相关联的所述加速度的能量的外点搏动时，从所述检测到的搏动集合中排除所述检测到的搏动；
当所述准时心率是相对于所述局部心率的内点心率、并且所述检测到的搏动是相对于与所述检测到的搏动相关联的所述加速度的能量的内点搏动时，将所述检测到的搏动包括在所述检测到的搏动集合中；和
当所述准时心率是相对于所述局部心率的外点心率、并且所述检测到的搏动是相对于与所述检测到的搏动相关联的所述加速度的能量的内点搏动时，从所述检测到的搏动集合中排除所述检测到的搏动并且将其分类为节律性异常搏动；以及
使用稀疏近似模型的字典来检测所述测试搏动集合中的异常搏动。
6.根据权利要求5所述的方法，包括：
基于与所述检测到的搏动集合中的搏动相关联的所述心率和所述加速度的能量，将搏动变换应用于所选择的搏动集合中的所述搏动，生成自适应的搏动；
基于与所述检测到的搏动集合中的所述搏动相关联的所述心率和所述加速度的能量，将字典变换应用于所述字典，生成自适应的字典；
生成相对于所述自适应的字典的所述自适应的搏动的稀疏近似；和
基于所述自适应的搏动的所生成的稀疏近似来将所述检测到的搏动集合中的所述搏动分类为异常。
7.一种电子设备，包括：
接口，所述接口在操作中接收测试心跳信号、以及一个或多个加速度信号；和
信号处理电路，所述信号处理电路被耦合到所述接口，并且所述信号处理电路在操作中：
检测所述测试心跳信号中的搏动；
测量所述一个或多个加速度信号的加速度的能量；
将心率和所测量的所述加速度的能量与检测到的搏动相关联；
基于与所述检测到的搏动相关联的所述心率和所述加速度的能量来将所述检测到的搏动选择性地包括在测试搏动集合中；和
使用稀疏近似模型的字典来检测所述测试搏动集合中的异常搏动，其中检测所述异常搏动包括：
基于与所述测试搏动集合中的搏动相关联的所述心率和所述加速度的能量，将搏动变换应用于所选择的搏动集合中的所述搏动，生成自适应的搏动；
基于与所述测试搏动集合中的所述搏动相关联的所述心率和所述加速度的能量，将字典变换应用于所述字典，生成自适应的字典；
生成相对于所述自适应的字典的所述自适应的搏动的稀疏近似；和
基于所生成的所述稀疏近似来生成指示所述测试搏动集合中的所述搏动是否异常的信号。
8.根据权利要求7所述的电子设备，其中，基于所生成的所述稀疏近似来生成指示所述测试搏动集合中的所述搏动是否异常的信号包括：
确定与所述自适应的搏动的所生成的所述稀疏近似值相关联的均方根误差；
将所述均方根误差与阈值误差进行比较；和
基于所述比较来生成所述信号。
9.根据权利要求7所述的电子设备，其中：
当所述生成的信号指示所述测试搏动集合中的所述搏动是异常的时，所述信号处理电路发起与所述测试搏动集合中的所述搏动有关的数据的传输。
10.根据权利要求7所述的电子设备，其中，所述心率是局部心率，并且将所述检测到的搏动选择性地包括在所述测试搏动集合中包括：
确定所述检测到的搏动的准时心率是否是相对于与所述检测到的搏动相关联的所述局部心率的外点心率；
确定所述检测到的搏动是否是相对于与所述检测到的搏动相关联的所述加速度的能量的外点搏动；
当所述准时心率是相对于所述局部心率的外点心率、或者所述检测到的搏动是相对于与所述检测到的搏动相关联的所述加速度的能量的外点搏动时，从所述检测到的搏动集合中排除所述检测到的搏动；和
当所述准时心率是相对于所述局部心率的内点心率、并且所述检测到的搏动是相对于与所述检测到的搏动相关联的所述加速度的能量的内点搏动时，将所述检测到的搏动包括在所述检测到的搏动集合中。
11.根据权利要求7所述的电子设备，其中所述心率是局部心率，并且将所述检测到的搏动选择性地包括在所述测试搏动集合中包括：
确定所述检测到的搏动的准时心率是否是相对于与所述检测到的搏动相关联的所述局部心率的外点心率；
确定所述检测到的搏动是否是相对于与所述检测到的搏动相关联的所述加速度的能量的外点搏动；
当所述准时心率是相对于所述局部心率的外点心率、并且所述检测到的搏动是相对于与所述检测到的搏动相关联的所述加速度的能量的外点搏动时，从所述检测到的搏动集合中排除所述检测到的搏动；以及
当所述准时心率是相对于所述局部心率的内点心率、或者所述检测到的搏动是相对于与所述检测到的搏动相关联的所述加速度的能量的内点搏动时，将所述检测到的搏动包括在所述检测到的搏动集合中。
12.根据权利要求7所述的电子设备，其中所述心跳信号是心电图ECG信号。
13.一种电子设备，包括：
接口，所述接口在操作中接收测试心跳信号、以及一个或多个加速度信号；和
信号处理电路，所述信号处理电路被耦合到所述接口，并且所述信号处理电路在操作中：
检测所述测试心跳信号中的搏动；
测量所述一个或多个加速度信号的加速度的能量；
将心率和所测量的所述加速度的能量与检测到的搏动相关联；
基于与所述检测到的搏动相关联的所述心率和所测量的所述加速度的能量来将所述检测到的搏动选择性地包括在测试搏动集合中；和
使用稀疏近似模型的字典来检测所述测试搏动集合中的异常搏动，其中：
所述接口在操作中接收训练心跳信号；并且
所述信号处理电路在操作中：
检测所述训练心跳信号的搏动；
将参考心率和加速度的能量与所述训练心跳信号的所述检测到的搏动相关联；
基于与所述训练心跳信号中的所述检测到的搏动相关联的所述心率和加速度的能量来将所述训练心跳信号的所述检测到的搏动选择性地包括在训练数据集合中；以及
使用所述训练数据集合来生成所述稀疏表示模型的所述字典。
14.根据权利要求13所述的电子设备，其中，所述信号处理电路通过迭代地进行以下各项来生成所述稀疏表示模型的所述字典：
生成稀疏表示矩阵的系数；和
优化所述字典、以及所述稀疏表示矩阵的所述系数。
15.根据权利要求14所述的电子设备，其中，所述参考心率与所述训练心跳信号相关联，并且将所述训练心跳信号的所述检测到的搏动选择性地包括在所述训练数据集合中包括：
确定所述训练心跳信号的所述检测到的搏动的准时心率是否是相对于所述参考心率的外点心率；
确定所述训练心跳信号的所述检测到的搏动是否是相对于与所述训练心跳信号的所述检测到的搏动相关联的所述加速度的能量的外点搏动；
当所述准时心率是相对于所述参考心率的外点心率时，从所述训练数据集合中排除所述训练心跳信号的所述检测到的搏动；和
当所述训练心跳信号的所述检测到的搏动是相对于与所述训练心跳信号的所述检测到的搏动相关联的所述加速度的能量的外点搏动时，从所述训练数据集合中排除所述训练心跳信号的所述检测到的搏动。
16.一种用于异常检测的系统，包括：
传感器，所述传感器在操作时生成心跳信号；
加速度计，所述加速度计在操作时生成一个或多个加速度信号；和
信号处理电路，所述信号处理电路被耦合到所述传感器和所述加速度计，其中所述信号处理电路在操作中：
检测测试心跳信号的搏动；
测量所述一个或多个加速度信号的加速度的能量；
将心率和所测量的所述加速度的能量与所述测试心跳信号的检测到的搏动相关联；
基于与所述测试心跳信号的所述检测到的搏动相关联的所述心率和所述加速度的能量来将所述测试心跳信号的所述检测到的搏动选择性地包括在测试搏动集合中；和
使用稀疏近似模型的字典来检测所述测试搏动集合中的异常搏动，其中，检测所述异常搏动包括：
基于与所述测试搏动集合中的搏动相关联的所述心率和所述加速度的能量，将搏动变换应用于所选择的搏动集合中的所述搏动，生成自适应的搏动；
基于与所述测试搏动集合中的所述搏动相关联的所述心率和所述加速度的能量，将字典变换应用于所述字典，生成自适应的字典；
生成相对于所述自适应的字典的所述自适应的搏动的稀疏近似；和
基于所生成的所述稀疏近似来生成指示所述测试搏动集合中的所述搏动是否异常的信号。
17.一种用于异常检测的系统，包括：
传感器，所述传感器在操作时生成心跳信号；
加速度计，所述加速度计在操作时生成一个或多个加速度信号；和
信号处理电路，所述信号处理电路被耦合到所述传感器和所述加速度计，其中所述信号处理电路在操作中：
检测测试心跳信号的搏动；
测量所述一个或多个加速度信号的加速度的能量；
将心率和所测量的所述加速度的能量与所述测试心跳信号的检测到的搏动相关联；
基于与所述测试心跳信号的所述检测到的搏动相关联的所述心率和所述加速度的能量来将所述测试心跳信号的所述检测到的搏动选择性地包括在测试搏动集合中；
使用稀疏近似模型的字典来检测所述测试搏动集合中的异常搏动；
检测训练心跳信号的搏动；
将参考心率和加速度的能量与所述训练心跳信号的所述检测到的搏动相关联；
基于与所述训练心跳信号中的所述检测到的搏动相关联的所述心率和加速度的能量来将所述训练心跳信号的所述检测到的搏动选择性地包括在训练数据集合中；并且
使用所述训练数据集合来生成所述稀疏表示模型的所述字典。
18.根据权利要求17所述的系统，其中，所述信号处理电路通过迭代地进行以下各项来生成所述稀疏表示模型的所述字典：
生成稀疏表示矩阵的系数；和
优化所述字典、以及所述稀疏表示矩阵的所述系数。
19.一种非暂态计算机可读介质，其具有配置心率监测设备以执行方法的内容，所述方法包括：
检测测试心跳信号的搏动；
测量加速度的能量；
将心率和所测量的所述加速度的能量与所述测试心跳信号的检测到的搏动相关联；
基于与所述测试心跳信号的所述检测到的搏动相关联的所述心率和所述加速度的能量来将所述测试心跳信号的所述检测到的搏动选择性地包括在测试搏动集合中；以及
使用稀疏近似模型的字典来检测所述测试搏动集合中的异常搏动，其中，所述异常搏动的所述检测包括：
基于与所述测试搏动集合中的搏动相关联的所述心率和所述加速度的能量，将搏动变换应用于所选择的搏动集合中的所述搏动，生成自适应的搏动；
基于与所述测试搏动集合中的所述搏动相关联的所述心率和所述加速度的能量，将字典变换应用于所述字典，生成自适应的字典；
生成相对于所述自适应的字典的所述自适应的搏动的稀疏近似；和
基于所生成的所述稀疏近似来生成指示所述测试搏动集合中的所述搏动是否异常的信号。
20.一种非暂态计算机可读介质，其具有配置心率监测设备以执行方法的内容，所述方法包括：
检测测试心跳信号的搏动；
测量加速度的能量；
将心率和所测量的所述加速度的能量与所述测试心跳信号的检测到的搏动相关联；
基于与所述测试心跳信号的所述检测到的搏动相关联的所述心率和所述加速度的能量来将所述测试心跳信号的所述检测到的搏动选择性地包括在测试搏动集合中；
使用稀疏近似模型的字典来检测所述测试搏动集合中的异常搏动；
检测训练心跳信号的搏动；
将参考心率和加速度的能量与所述训练心跳信号的所述检测到的搏动相关联；
基于与所述训练心跳信号中的所述检测到的搏动相关联的所述心率和加速度的能量来将所述训练心跳信号的所述检测到的搏动选择性地包括在训练数据集合中；和
使用所述训练数据集合来生成所述稀疏表示模型的所述字典。
21.根据权利要求20所述的非暂态计算机可读介质，其中，所述方法包括通过迭代地进行以下各项来生成所述稀疏表示模型的所述字典：
生成稀疏表示矩阵的系数；和
优化所述字典、以及所述稀疏表示矩阵的所述系数。
22.一种用于异常检测的系统，包括：
用于生成心跳信号的装置；
用于生成加速度信号的装置；和
用于通过以下各项来处理测试心跳信号的装置：
检测所述测试心跳信号的搏动；
将心率和加速度的能量与所述测试心跳信号的检测到的搏动相关联；
基于与所述测试心跳信号的所述检测到的搏动相关联的所述心率和所述加速度的能量来将所述测试心跳信号的所述检测到的搏动选择性地包括在测试搏动集合中；和
使用稀疏近似模型的字典来检测所述测试搏动集合中的异常搏动，所述异常搏动的所述检测包括：
基于与所述测试搏动集合中的搏动相关联的所述心率和所述加速度的能量，将搏动变换应用于所选择的搏动集合中的所述搏动，生成自适应的搏动；
基于与所述测试搏动集合中的所述搏动相关联的所述心率和所述加速度的能量，将字典变换应用于所述字典，生成自适应的字典；
生成相对于所述自适应的字典的所述自适应的搏动的稀疏近似；和
基于所生成的所述稀疏近似来生成指示所述测试搏动集合中的所述搏动是否异常的信号。
23.根据权利要求22所述的系统，包括：
用于通过以下各项来处理训练心跳信号的装置：
检测所述训练心跳信号的搏动；
将参考心率和加速度的能量与所述训练心跳信号的所述检测到的搏动相关联；
基于与所述训练心跳信号中的所述检测到的搏动相关联的所述心率和加速度的能量来将所述训练心跳信号的所述检测到的搏动选择性地包括在所述训练数据集合中；和
使用所述训练数据集合来生成所述稀疏表示模型的所述字典。
24.根据权利要求23所述的系统，其中，所述处理所述训练心跳信号包括通过迭代地进行以下各项来生成所述稀疏表示模型的所述字典：
生成稀疏表示矩阵的系数；和
优化所述字典、以及所述稀疏表示矩阵的所述系数。
CN201611227015.7A 2016-05-31 2016-12-27 用于检测心电图异常的方法及对应系统 Active CN107440679B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US15/169,184 2016-05-31
US15/169,184 US10595788B2 (en) 2016-05-31 2016-05-31 Method for the detecting electrocardiogram anomalies and corresponding system

Publications (2)

Publication Number Publication Date
CN107440679A CN107440679A (zh) 2017-12-08
CN107440679B true CN107440679B (zh) 2020-09-01

Family

ID=57737655

Family Applications (2)

Application Number Title Priority Date Filing Date
CN201621447161.6U Active CN207084814U (zh) 2016-05-31 2016-12-27 用于检测异常搏动的设备和系统
CN201611227015.7A Active CN107440679B (zh) 2016-05-31 2016-12-27 用于检测心电图异常的方法及对应系统

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN201621447161.6U Active CN207084814U (zh) 2016-05-31 2016-12-27 用于检测异常搏动的设备和系统

Country Status (3)

Country Link
US (1) US10595788B2 (zh)
EP (1) EP3251594B1 (zh)
CN (2) CN207084814U (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11915825B2 (en) * 2017-02-10 2024-02-27 Alivecor, Inc. Systems and methods of analyte measurement analysis
CN108852340B (zh) * 2018-05-17 2021-09-03 广东工业大学 一种心电采集系统、装置及计算机可读存储介质
CN109480825B (zh) * 2018-12-13 2021-08-06 武汉中旗生物医疗电子有限公司 心电数据的处理方法及装置
CN110502995B (zh) * 2019-07-19 2023-03-14 南昌大学 基于细微面部动作识别的驾驶员打哈欠检测方法
CN110693486B (zh) * 2019-09-27 2022-06-14 武汉中旗生物医疗电子有限公司 一种心电图的异常标注方法及装置
CN110652295B (zh) * 2019-10-16 2022-04-26 中山大学 一种量化激活时间序列个体差异性的方法及系统
CN112801434A (zh) * 2019-11-13 2021-05-14 中兴通讯股份有限公司 性能指标健康度的监测方法、装置、设备和存储介质
CN111832427B (zh) * 2020-06-22 2022-02-18 华中科技大学 基于欧氏对齐和Procrustes分析的EEG分类的迁移学习方法和系统
CN112232253B (zh) * 2020-10-26 2021-11-09 南通大学 一种粗糙合理粒度驱动的胎儿心电图异常检测方法
GB2612979A (en) * 2021-11-17 2023-05-24 Prevayl Innovations Ltd Method and system for measuring and displaying biosignal data to a wearer of a wearable article

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008065432A1 (en) * 2006-12-01 2008-06-05 Oxford Biosignals Limited Biomedical signal morphology analysis method
CN105101870A (zh) * 2013-03-29 2015-11-25 皇家飞利浦有限公司 用于ecg运动伪影去除的装置和方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8614630B2 (en) 2011-11-14 2013-12-24 Vital Connect, Inc. Fall detection using sensor fusion
US9289150B1 (en) 2012-08-17 2016-03-22 Analytics For Life Non-invasive method and system for characterizing cardiovascular systems
CA2900160C (en) 2013-02-05 2019-05-21 Vivaquant Llc Detecting fiducial points in physiological signals
US9931081B2 (en) * 2013-03-29 2018-04-03 Koninklijke Philips N.V. Apparatus and method for ECG motion artifact removal
EP3079571A4 (en) 2013-12-12 2017-08-02 Alivecor, Inc. Methods and systems for arrhythmia tracking and scoring

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008065432A1 (en) * 2006-12-01 2008-06-05 Oxford Biosignals Limited Biomedical signal morphology analysis method
CN105101870A (zh) * 2013-03-29 2015-11-25 皇家飞利浦有限公司 用于ecg运动伪影去除的装置和方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Advanced K-Means Clustering Algorithm for Large ECG Data Sets Based on K-SVD Approach;，Mohammadreza Balouchestani等;《2014 9th International Symposium on Communication Systems, Networks & Digital Sign (CSNDSP),IEEE》;20140723;第177-182页 *

Also Published As

Publication number Publication date
US20170340292A1 (en) 2017-11-30
EP3251594B1 (en) 2024-02-28
CN107440679A (zh) 2017-12-08
EP3251594A1 (en) 2017-12-06
US10595788B2 (en) 2020-03-24
CN207084814U (zh) 2018-03-13

Similar Documents

Publication Publication Date Title
CN107440679B (zh) 用于检测心电图异常的方法及对应系统
US10610162B2 (en) Method for the detecting electrocardiogram anomalies and corresponding system
US10869610B2 (en) System and method for identifying cardiac arrhythmias with deep neural networks
CN107714023B (zh) 基于人工智能自学习的静态心电图分析方法和装置
Wang et al. Arrhythmia classification algorithm based on multi-head self-attention mechanism
JP6367442B2 (ja) 診断信号の変換に基づいた疾病解析のための方法及びシステム
Deng et al. Extracting cardiac dynamics within ECG signal for human identification and cardiovascular diseases classification
US11589828B2 (en) System and methods for electrocardiogram beat similarity analysis using deep neural networks
Li et al. Signal processing and feature selection preprocessing for classification in noisy healthcare data
Bulusu et al. Transient ST-segment episode detection for ECG beat classification
Gavrovska et al. Automatic heart sound detection in pediatric patients without electrocardiogram reference via pseudo-affine Wigner–Ville distribution and Haar wavelet lifting
US8750973B2 (en) Method and system for detecting P-waves in the surface ECG signal
US20230181082A1 (en) System and methods for electrocardiogram beat similarity analysis
Zidelmal et al. Heartbeat classification using support vector machines (SVMs) with an embedded reject option
US11672464B2 (en) Electrocardiogram processing system for delineation and classification
KR20140097039A (ko) 자동 결합 뉴런 네트워크를 이용한 심장부정맥 분류 방법 및 장치
CN116636856A (zh) 心电分类方法、系统、电子设备及存储介质
Bashir et al. Highlighting the current issues with pride suggestions for improving the performance of real time cardiac health monitoring
Wannawijit et al. ECG classification with modification of higher-order Hjorth descriptors
Natarajan et al. Conditional random fields for morphological analysis of wireless ecg signals
Salah et al. Cardiac arrhythmia classification by wavelet transform
Pham Classification of ecg signals of heart beats using tf-ts lstm with augmented fuzzy recurrence eigenvalues
Liberczuk et al. Model Based Feature Extraction Method for Myocardial Infarction Detection
Hosseinpoor Intelligent Diagnosis of Heart Diseases Based on Electrocardiographic Signal
Konstantin et al. Noise-resilient Automatic Interpretation of Holter ECG Recordings

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