CN111297351A - Motion artifact identification method and device in dynamic electrocardiogram - Google Patents
Motion artifact identification method and device in dynamic electrocardiogram Download PDFInfo
- Publication number
- CN111297351A CN111297351A CN202010091000.2A CN202010091000A CN111297351A CN 111297351 A CN111297351 A CN 111297351A CN 202010091000 A CN202010091000 A CN 202010091000A CN 111297351 A CN111297351 A CN 111297351A
- Authority
- CN
- China
- Prior art keywords
- heart beat
- motion artifact
- heartbeat
- width
- wave
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Cardiology (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The invention discloses a method and a device for identifying motion artifact in a dynamic electrocardiogram. Firstly, extracting characteristic parameters of each heart beat according to the heart beats obtained by detection; then, calculating the noise condition in a fixed time width before and after the heartbeat, and calculating to obtain indexes such as a signal-to-noise ratio, a confidence coefficient and the like of the heartbeat; then, according to a set threshold value, recognizing the heart beat with the low signal-to-noise ratio or confidence level as a motion artifact; and finally, combining the pseudo-difference conditions of the heart beats before and after each heart beat, and regarding the heart beats which are not successfully identified in the continuous pseudo-differences as motion pseudo-differences to finally obtain a motion pseudo-difference identification result. The method solves the problem that the existing method can not quickly prepare for recognizing the motion artifact of the heart beat, extracts the signal-to-noise ratio and the confidence index which can reflect the noise level according to the characteristic parameters of the heart beat on the basis of the detected heart beat, further recognizes the motion artifact condition by combining the noise conditions of the front heart beat and the back heart beat, and finally automatically eliminates the invalid heart beat.
Description
Technical Field
The invention belongs to the technical field of medical tests and discloses a method and a device for identifying motion artifacts in a dynamic electrocardiogram.
Background
The dynamic electrocardiogram is a body surface electrocardiogram recorded continuously for 24 hours or more, and the information content is much larger than that of the conventional electrocardiogram, thereby increasing the workload of electrocardiogram analysis, and making the automatic detection and analysis technology of the electrocardiogram necessary. In addition, in the process of acquiring the electrocardiogram, various interferences of daily activities, especially the influence of motion interferences, are easily received, so that various motion artifact interferences are doped in the actually acquired electrocardiogram waveform, most of the artifacts are invalid signals, and the artifacts need to be removed and analyzed, and if the artifacts are removed only by hand, the workload of doctors is undoubtedly and greatly increased. Therefore, the method for rapidly and accurately identifying the motion artifact interference is an effective means for improving the working efficiency of doctors.
Since the diagnosis of the electrocardiogram is performed on a heartbeat basis, motion artifacts are generally recognized on a heartbeat basis. At present, the existing products in the market mainly identify the motion artifact through manual deletion assisted by a histogram or a heart beat template, and the processing efficiency is generally low.
Disclosure of Invention
The invention aims to provide a motion artifact identification method in a dynamic electrocardiogram aiming at the defect that the motion artifact signals in the electrocardiogram can not be quickly and accurately identified in the prior art so as to achieve the purposes of enabling the motion artifact identification to be more accurate and reliable and effectively completing the identification and detection of the motion artifacts.
The invention relates to a motion artifact identification method in a dynamic electrocardiogram, which comprises the following steps:
(1) extracting characteristic parameters of each heartbeat on the basis of the dynamic electrocardiogram which is detected by the heartbeats; the characteristic parameters comprise area A, width W, R wave height H1Q wave height H2Height H of S wave3Etc.;
(2) according to the heart beats obtained by detection, the average area A of the 200 millisecond width area before and after each heart beat is obtainedT；
(3) According to the characteristic parameters of each heartbeat in the step (1) and the average area A of 200 millisecond width areas before and after each heartbeat in the step (2)TCalculating indexes such as signal-to-noise ratio and confidence degree of the heart beat, and considering the heart beat as motion artifact according to the condition that ① signal-to-noise ratio is smaller than a threshold value of 1.5 and ② confidence degree is smaller than a threshold value of 40;
the heart beat signal-to-noise ratio index SNR is calculated as follows:
if the SNR is less than 1.5, the detection heartbeat is considered as a motion artifact;
cardiac beat confidence index RconfThe calculation formula is as follows:
if the heartbeat width W is less than 60 milliseconds or greater than 200 milliseconds, RconfSet to 0;
height of heart beat H1Is less than (H)2+H3) Then R isconfSet to 0;
if the confidence degree R isconfIf the interference is less than 40, the heart beat artifact interference is considered to be too large, and the heart beat artifact interference is judged to be motion artifact;
(4) from the first heart beat, if the current heart beat is a non-motion artifact, calculating the ratio of the number of heart beat artifacts within N seconds before and after the current heart beat, and if the artifact ratio is greater than a certain threshold, regarding the current heart beat as a motion artifact;
(5) and (5) repeating the step (4) twice to finish the motion artifact identification processing.
A motion artifact identification device in a dynamic electrocardiogram, comprising:
the characteristic extraction unit is used for extracting the characteristic parameters of each heartbeat in the dynamic electrocardiogram which is subjected to heartbeat detection; the characteristic parameters comprise area A, width W, R wave height H1Q wave height H2Height H of S wave3Etc.;
a first calculation unit for acquiring an average area A of 200 msec width regions before and after each heartbeatT；ATIs T1point-to-Q wave and S-to-T wave2The integral of the amplitude of the dot divided by the width (400 milliseconds);
the second calculation unit is used for calculating the heart beat signal-to-noise ratio and the confidence coefficient of each heart beat;
the first motion artifact identification unit is used for preliminarily judging whether the current heart beat is motion artifact or not, and specifically, if any one of the following conditions is met, the heart beat is considered as motion artifact, if the signal-to-noise ratio of ① is smaller than the threshold value of 1.5, the heart beat is considered as motion artifact, and if the confidence coefficient of ② is smaller than the threshold value of 40, the heart beat is considered as motion artifact;
the second motion artifact identification unit is used for judging whether the non-motion artifact heartbeat identified by the first motion artifact identification unit is motion artifact again; specifically, calculating the heart beat artifact number ratio in N seconds before and after the current heart beat, and if the artifact ratio is greater than a certain threshold, regarding the current heart beat as a motion artifact;
the method solves the problem that the existing method can not quickly prepare for recognizing the motion artifact of the heart beat, extracts the signal-to-noise ratio and the confidence index which can reflect the noise level according to the characteristic parameters of the heart beat on the basis of the detected heart beat, further recognizes the motion artifact condition by combining the noise conditions of the front heart beat and the back heart beat, and finally automatically eliminates the invalid heart beat.
Drawings
FIG. 1 is a typical electrocardiographic waveform and characteristic point and parameter labels, wherein the horizontal dashed line represents the baseline, W represents the heart beat width, H1/H2/H3Respectively representing the heights of R wave/Q wave/S wave, the heart beat area A is the amplitude integral of Q wave to S wave divided by the heart beat width W, and the average area A of 200 milliseconds before and after the heart beatTIs T1point-to-Q wave and S-to-T wave2The integral of the amplitude of the dot divided by the width (400 msec).
Detailed Description
The invention will be further described with reference to the accompanying drawings in which:
the motion artifact identification method in the dynamic electrocardiogram comprises the following steps:
(1) extracting characteristic parameters of each heartbeat on the basis of the dynamic electrocardiogram which is detected by the heartbeats, wherein the characteristic parameters comprise heartbeat area A and heartbeat width W, R wave height H1Q wave height H2Height H of S wave3Etc.;
(2) according to the heart beats obtained by detection, calculating the average area A of 200 milliseconds before and after each heart beatT；
(3) Extracting indexes such as heart beat signal-to-noise ratio and confidence coefficient according to the calculated characteristic parameters, wherein the heart beat signal-to-noise ratio index SNR calculation formula is as follows:
if the SNR is less than 1.5, the detection heartbeat is considered as a motion artifact;
cardiac beat confidence index RconfThe calculation formula is as follows:
further, if the heartbeat width W is less than 60 milliseconds or greater than 200 milliseconds, then RconfSet to 0;
further, if the height of the heart beat H1Is less than (H)2+H3) Then R isconfSet to 0;
if the confidence degree R isconfIf the detected heartbeat is less than 40, judging the detected heartbeat as a motion artifact;
(4) calculating the ratio of the number of heart beat pseudo-errors in N seconds before and after the current heart beat from the first heart beat if the current heart beat is a non-motion pseudo-error, and regarding the current heart beat as a motion pseudo-error if the pseudo-error ratio is greater than 70%;
(5) and (5) repeating the step (4) twice to finish the motion artifact identification processing.
A motion artifact identification device in a dynamic electrocardiogram, comprising:
the characteristic extraction unit is used for extracting the characteristic parameters of each heartbeat in the dynamic electrocardiogram which is subjected to heartbeat detection; the characteristic parameters comprise area A, width W, R wave height H1Q wave height H2Height H of S wave3Etc.;
a first calculation unit for acquiring an average area A of 200 msec width regions before and after each heartbeatT；ATIs T1point-to-Q wave and S-to-T wave2The integral of the amplitude of the dot divided by the width (400 milliseconds);
the second calculation unit is used for calculating the heart beat signal-to-noise ratio and the confidence coefficient of each heart beat;
the first motion artifact identification unit is used for preliminarily judging whether the current heart beat is motion artifact or not, and specifically, if any one of the following conditions is met, the heart beat is considered as motion artifact, if the signal-to-noise ratio of ① is smaller than the threshold value of 1.5, the heart beat is considered as motion artifact, and if the confidence coefficient of ② is smaller than the threshold value of 40, the heart beat is considered as motion artifact;
the second motion artifact identification unit is used for judging whether the non-motion artifact heartbeat identified by the first motion artifact identification unit is motion artifact again; specifically, calculating the heart beat artifact number ratio in N seconds before and after the current heart beat, and if the artifact ratio is greater than a certain threshold, regarding the current heart beat as a motion artifact;
FIG. 1 is a typical electrocardiographic waveform and characteristic point and parameter labels, wherein the horizontal dashed line represents the baseline, W represents the heart beat width, H1/H2/H3Respectively representing the heights of R wave/Q wave/S wave, the heart beat area A is the amplitude integral of Q wave to S wave divided by the heart beat width W, and the average area A of 200 milliseconds before and after the heart beatTIs T1point-to-Q wave and S-to-T wave2The integral of the amplitude of the dot divided by the width (400 msec).
Claims (2)
1. A method for recognizing motion artifact in a dynamic electrocardiogram is characterized by comprising the following steps:
(1) extracting characteristic parameters of each heartbeat on the basis of the dynamic electrocardiogram which is detected by the heartbeats; the characteristic parameters comprise area A, width W, R wave height H1Q wave height H2Height H of S wave3Etc.;
(2) according to the heart beats obtained by detection, the average area A of the 200 millisecond width area before and after each heart beat is obtainedT；
(3) According to the characteristic parameters of each heartbeat in the step (1) and the average area A of 200 millisecond width areas before and after each heartbeat in the step (2)TCalculating indexes such as signal-to-noise ratio and confidence degree of the heart beat, and considering the heart beat as motion artifact according to the condition that ① signal-to-noise ratio is less than 1.5 of threshold value②, if the confidence coefficient is less than the threshold value 40, the motion artifact is considered;
calculating formula of heart beat signal-to-noise ratio SNR index:
cardiac beat confidence index RconfThe calculation formula is as follows:
if the heartbeat width W is less than 60 milliseconds or the heartbeat width W is more than 200 milliseconds, Rconf＝0；
Height of heart beat H1＜(H2+H3) Then R isconf＝0；
If the heartbeat width W is less than or equal to 200 milliseconds and the heartbeat height H is less than or equal to 60 milliseconds1≥(H2+H3) Then, then
(4) From the first heart beat, if the current heart beat is a non-motion artifact, calculating the ratio of the number of heart beat artifacts within N seconds before and after the current heart beat, and if the artifact ratio is greater than a certain threshold, regarding the current heart beat as a motion artifact;
(5) and (5) repeating the step (4) twice to finish the motion artifact identification processing.
2. An apparatus for recognizing motion artifact in a dynamic electrocardiogram, comprising:
the characteristic extraction unit is used for extracting the characteristic parameters of each heartbeat in the dynamic electrocardiogram which is subjected to heartbeat detection; the characteristic parameters comprise area A, width W, R wave height H1Q wave height H2Height H of S wave3Etc.;
a first calculation unit for acquiring an average area A of 200 msec width regions before and after each heartbeatT；ATIs T1point-to-Q wave and S-to-T wave2The integral of the amplitude of the dot divided by the width (400 milliseconds);
the second calculation unit is used for calculating the heart beat signal-to-noise ratio and the confidence coefficient of each heart beat;
the first motion artifact identification unit is used for preliminarily judging whether the current heart beat is motion artifact or not, and specifically, if any one of the following conditions is met, the heart beat is considered as motion artifact, if the signal-to-noise ratio of ① is smaller than the threshold value of 1.5, the heart beat is considered as motion artifact, and if the confidence coefficient of ② is smaller than the threshold value of 40, the heart beat is considered as motion artifact;
the second motion artifact identification unit is used for judging whether the non-motion artifact heartbeat identified by the first motion artifact identification unit is motion artifact again; specifically, the number ratio of the heart beat artifact within N seconds before and after the current heart beat is calculated, and if the artifact ratio is larger than a certain threshold, the current heart beat is regarded as the motion artifact.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010091000.2A CN111297351A (en) | 2020-02-13 | 2020-02-13 | Motion artifact identification method and device in dynamic electrocardiogram |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010091000.2A CN111297351A (en) | 2020-02-13 | 2020-02-13 | Motion artifact identification method and device in dynamic electrocardiogram |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111297351A true CN111297351A (en) | 2020-06-19 |
Family
ID=71161600
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010091000.2A Pending CN111297351A (en) | 2020-02-13 | 2020-02-13 | Motion artifact identification method and device in dynamic electrocardiogram |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111297351A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115399784A (en) * | 2022-09-02 | 2022-11-29 | 厦门纳龙健康科技股份有限公司 | Automatic shielding method for invalid electrocardiogram data, terminal equipment and storage medium |
| CN115778402A (en) * | 2022-12-02 | 2023-03-14 | 深圳华清心仪医疗电子有限公司 | Method and system for identifying artifact of dynamic electrocardiosignal |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002253546A (en) * | 2001-02-20 | 2002-09-10 | Ge Medical Systems Global Technology Co Llc | Artifact evaluation method and program, and x-ray ct device |
| JP2008182487A (en) * | 2007-01-24 | 2008-08-07 | Sony Corp | Image processing device, and image processing method, and program |
| CN103720468A (en) * | 2013-12-05 | 2014-04-16 | 深圳先进技术研究院 | Artifact identification method and device applied to dynamic electrocardiogram data |
| CN104545887A (en) * | 2014-12-24 | 2015-04-29 | 深圳先进技术研究院 | Method and device for identifying artifact electrocardiograph waveforms |
| CN106485213A (en) * | 2016-09-27 | 2017-03-08 | 鲁东大学 | A kind of utilization electrocardiosignal carries out the feature extracting method of automatic identification |
| CN108701219A (en) * | 2017-03-14 | 2018-10-23 | 华为技术有限公司 | The method and device of waveform signal processing |
| CN109411042A (en) * | 2018-02-24 | 2019-03-01 | 深圳市凯沃尔电子有限公司 | Ecg information processing method and electro cardio signal workstation |
| US20190298210A1 (en) * | 2013-07-01 | 2019-10-03 | Mayo Foundation For Medical Education And Research | Algorithms for managing artifact and detecting cardiac events using a patient monitoring system |
-
2020
- 2020-02-13 CN CN202010091000.2A patent/CN111297351A/en active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002253546A (en) * | 2001-02-20 | 2002-09-10 | Ge Medical Systems Global Technology Co Llc | Artifact evaluation method and program, and x-ray ct device |
| JP2008182487A (en) * | 2007-01-24 | 2008-08-07 | Sony Corp | Image processing device, and image processing method, and program |
| US20190298210A1 (en) * | 2013-07-01 | 2019-10-03 | Mayo Foundation For Medical Education And Research | Algorithms for managing artifact and detecting cardiac events using a patient monitoring system |
| CN103720468A (en) * | 2013-12-05 | 2014-04-16 | 深圳先进技术研究院 | Artifact identification method and device applied to dynamic electrocardiogram data |
| CN104545887A (en) * | 2014-12-24 | 2015-04-29 | 深圳先进技术研究院 | Method and device for identifying artifact electrocardiograph waveforms |
| CN106485213A (en) * | 2016-09-27 | 2017-03-08 | 鲁东大学 | A kind of utilization electrocardiosignal carries out the feature extracting method of automatic identification |
| CN108701219A (en) * | 2017-03-14 | 2018-10-23 | 华为技术有限公司 | The method and device of waveform signal processing |
| CN109411042A (en) * | 2018-02-24 | 2019-03-01 | 深圳市凯沃尔电子有限公司 | Ecg information processing method and electro cardio signal workstation |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115399784A (en) * | 2022-09-02 | 2022-11-29 | 厦门纳龙健康科技股份有限公司 | Automatic shielding method for invalid electrocardiogram data, terminal equipment and storage medium |
| CN115778402A (en) * | 2022-12-02 | 2023-03-14 | 深圳华清心仪医疗电子有限公司 | Method and system for identifying artifact of dynamic electrocardiosignal |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2979135C (en) | Systems, apparatus and methods for sensing fetal activity | |
| CN113827215B (en) | Automatic diagnosis method for various arrhythmias based on millimeter wave radar | |
| US6597943B2 (en) | Method of using spectral measures to distinguish among atrialfibrillation, atrial-flutter and other cardiac rhythms | |
| Rakshit et al. | An efficient wavelet-based automated R-peaks detection method using Hilbert transform | |
| CN109117730B (en) | Real-time electrocardiogram atrial fibrillation judgment method, device and system and storage medium | |
| CN109303559B (en) | Dynamic electrocardiogram and heartbeat classification method based on gradient boosting decision tree | |
| CN113057648A (en) | ECG signal classification method based on composite LSTM structure | |
| CN113907765B (en) | Noninvasive fetal electrocardiosignal quality assessment method | |
| CN108514414A (en) | A kind of signal handling equipment, human ecg signal processing and disease forecasting method | |
| CN112617849A (en) | Atrial fibrillation detection and classification method based on CNN + LSTM | |
| US20240389929A1 (en) | Fusion signal processing for maternal uterine activity detection | |
| CN115486855B (en) | Electrocardiogram and heart beat classification method based on QRS complex indefinite cycle superposition | |
| CN111297351A (en) | Motion artifact identification method and device in dynamic electrocardiogram | |
| CN105877739A (en) | Clinical examination method of electrocardio intelligent analyzing system | |
| CN116784860B (en) | Electrocardiosignal characteristic extraction system based on morphological heart beat template clustering | |
| Karvounis et al. | Detection of fetal heart rate through 3-D phase space analysis from multivariate abdominal recordings | |
| CN115486854A (en) | Single lead electrocardiogram ventricular premature beat identification method aiming at dry electrode acquisition | |
| Di Maria et al. | An algorithm for the analysis of fetal ECGs from 4-channel non-invasive abdominal recordings | |
| CN112017695A (en) | A system and method for automatically identifying physiological sounds | |
| CN110491504B (en) | Method for acquiring medical index data of heart sound signal | |
| CN110367973B (en) | Method for detecting initial point of QRS wave group of multi-lead electrocardiosignal | |
| CN108836312B (en) | Clutter rejection method and system based on artificial intelligence | |
| Manlong et al. | ST segment analysis for early detection of myocardial infarction | |
| CN111657933A (en) | Autonomous high-frequency QRS wave group analysis device and analysis method | |
| CN115105084B (en) | Heartbeat location classification method, system, medium and device based on deep learning |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200619 |