WO2022215191A1 - Procédé de détection de battement cardiaque et dispositif de détection de battement cardiaque - Google Patents

Procédé de détection de battement cardiaque et dispositif de détection de battement cardiaque Download PDF

Info

Publication number
WO2022215191A1
WO2022215191A1 PCT/JP2021/014756 JP2021014756W WO2022215191A1 WO 2022215191 A1 WO2022215191 A1 WO 2022215191A1 JP 2021014756 W JP2021014756 W JP 2021014756W WO 2022215191 A1 WO2022215191 A1 WO 2022215191A1
Authority
WO
WIPO (PCT)
Prior art keywords
value
threshold
time
heartbeat
sampling time
Prior art date
Application number
PCT/JP2021/014756
Other languages
English (en)
Japanese (ja)
Inventor
伸昭 松浦
啓 桑原
利彦 近藤
Original Assignee
日本電信電話株式会社
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 日本電信電話株式会社 filed Critical 日本電信電話株式会社
Priority to PCT/JP2021/014756 priority Critical patent/WO2022215191A1/fr
Priority to JP2023512570A priority patent/JP7491466B2/ja
Publication of WO2022215191A1 publication Critical patent/WO2022215191A1/fr

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
    • A61B5/352—Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval

Definitions

  • the present invention relates to a heartbeat detection method and a heartbeat detection device for detecting heartbeats (R waves) from an electrocardiogram waveform.
  • the ECG (Electrocardiogram) waveform is an observation of the electrical activity of the heart and consists of a continuous heartbeat waveform.
  • One heartbeat waveform consists of components such as P wave, Q wave, R wave, S wave, T wave, etc., which reflect the activity of the atria and ventricles.
  • the R wave accompanies the contraction of the ventricle, and since it has a large amplitude, the heartbeat is often detected based on the R wave.
  • a simple method for detecting heartbeats (R waves) from ECG waveforms is to detect peaks in a time series. That is, a certain threshold is set for the data time series, and when the ECG waveform exceeds the threshold, it is determined as an R wave.
  • FIG. 5 is a flow chart explaining a conventional heartbeat detection method.
  • the heartbeat detection device calculates an index value based on the ECG waveform value and the previous ECG waveform value (step S101 in FIG. 5).
  • the heartbeat detection device compares the index value and the threshold value (step S102 in FIG. 5), determines a heartbeat when the index value exceeds the threshold value (step S103 in FIG. 5), and determines that the index value is equal to or less than the threshold value. Occasionally, it is not determined as a heartbeat.
  • the heartbeat detection device updates the threshold using the index value and the previous index value (step S104 in FIG. 5). By adaptively updating the threshold value, it is possible to perform heartbeat detection that reflects changes in the signal level due to, for example, the wearing state of electrodes on the subject.
  • Patent Document 1 As such a conventional heartbeat detection method, there is a method disclosed in Patent Document 1.
  • the index value is not reflected in the threshold when it is unlikely that the index value is derived from the ECG waveform.
  • the threshold value is not reflected in the threshold when it is unlikely that the index value is derived from the ECG waveform.
  • the present invention has been made to solve the above problems, and an object of the present invention is to provide a heartbeat detection method and a heartbeat detection device that do not erroneously recognize peaks with too small amplitudes as heartbeats.
  • a heartbeat detection method of the present invention includes a first step of calculating an index value from a sampled data string of an electrocardiogram waveform of a living body at each sampling time, and a threshold value for heartbeat detection and a threshold derived from the electrocardiogram waveform of a living body. a second step of comparing with the lower limit of the likely range; and determining in the second step that the threshold is equal to or greater than the lower limit, and when a peak of the index value exceeding the threshold is detected, the peak and a third step of setting the sampling time to the heartbeat time, and not performing heartbeat detection by comparing the threshold and the index value when the threshold is determined to be less than the lower limit in the second step. It is characterized.
  • the heartbeat detection apparatus of the present invention includes: an index value calculation unit configured to calculate an index value from a sampled data string of an electrocardiogram waveform of a living body at each sampling time; a threshold determination unit configured to compare the lower limit of a plausible range as the threshold derived from the threshold determination unit determines that the threshold is equal to or greater than the lower limit, and the index value exceeding the threshold a heartbeat time determination unit configured to set the sampling time of the peak as the heartbeat time when a peak is detected, wherein the heartbeat time determination unit determines that the threshold value is less than the lower limit value. , heartbeat detection by comparison between the threshold value and the index value is not performed.
  • Heart rate detection can be performed on
  • FIG. 1 is a block diagram showing the configuration of a heartbeat detection device according to an embodiment of the present invention.
  • FIG. 2 is a flowchart illustrating a heart rate detection method according to an embodiment of the invention.
  • FIG. 3 is a flow chart for explaining the operation of the threshold value setting section of the heartbeat detection device according to the embodiment of the present invention.
  • FIG. 4 is a block diagram showing a configuration example of a computer that implements the heartbeat detection device according to the embodiment of the present invention.
  • FIG. 5 is a flowchart illustrating a conventional heartbeat detection method.
  • FIG. 1 is a block diagram showing the configuration of a heartbeat detection device according to an embodiment of the invention
  • FIG. 2 is a flow chart for explaining a heartbeat detection method according to an embodiment of the invention.
  • the heartbeat detection apparatus includes an electrocardiograph 1 for outputting an ECG waveform sampling data string, a storage unit 2 for storing an ECG waveform sampling data string and sampling time information, and sampling data from the ECG waveform sampling data string.
  • a positive/negative inversion value calculation unit 3 that calculates the positive/negative inversion value of the time difference for each sampling time, and a positive/negative inversion value in a certain time range before the sampling time to be processed and a certain value after the sampling time to be processed.
  • a maximum value detection unit 4 that detects the maximum value of the positive/negative reversed values in the time range at each sampling time, and an index that calculates an index value obtained by subtracting the maximum value from the positive/negative reversed value at the sampling time to be processed at each sampling time.
  • a value calculation unit 5 a threshold determination unit 6 that compares the threshold value for heartbeat detection with the lower limit value of a range that is likely to be a threshold value derived from the ECG waveform of the living body, and a threshold value that is determined to be equal to or higher than the lower limit value, and exceeds the threshold value.
  • a heartbeat time determining unit 7 for setting the sampling time of the peak when the peak of the index value is detected to be the heartbeat time, and a threshold setting unit 8 for updating the threshold based on the peak of the index value exceeding the current threshold.
  • the maximum value detection unit 4 includes a FIFO buffer (First In, First Out) 40 to which the time difference positive/negative inversion value calculated by the positive/negative inversion value calculation unit 3 is input, and a FIFO buffer to which the output value of the FIFO buffer 40 is input. 41, a FIFO buffer 42 to which the output value of the FIFO buffer 41 is input, and sampling the maximum value of the positive/negative inverted value of the time difference stored in the FIFO buffer 40 and the inverted positive/negative value of the time difference stored in the FIFO buffer 42. It is composed of a detection processing unit 43 that detects each time.
  • a detection processing unit 43 that detects each time.
  • the index value calculation unit 5 detects the maximum value detected by the maximum value detection unit 4 from the FIFO buffer 50 to which the time difference positive/negative inversion value calculated by the positive/negative value calculation unit 3 is input, and the output value of the FIFO buffer 50. and a subtraction processing unit 51 that calculates an index value obtained by subtracting , at each sampling time.
  • the heartbeat detection method of this embodiment will be described below.
  • the procedure from detecting one heartbeat to obtaining the heartbeat time will be described.
  • Time-series data of heartbeat times is obtained by repeating such calculation of heartbeat times over the period of the ECG waveform data.
  • D(i) is the data string obtained by sampling the ECG waveform.
  • the electrocardiograph 1 measures an ECG waveform of a living body (human body) (not shown) and outputs a sampling data string D(i) of the ECG waveform. At this time, the electrocardiograph 1 adds sampling time information to each sampling data and outputs the data. Since the specific method of measuring the ECG waveform is a well-known technique, detailed description thereof will be omitted.
  • the storage unit 2 stores the sampling data sequence D(i) of the ECG waveform output from the electrocardiograph 1 and the information on the sampling time.
  • the positive/negative inversion value calculation unit 3 calculates the data D(i+1) after one sampling of the sampling data D(i) and the data D(i+1) one sampling before.
  • Data D(i-1) are obtained from the storage unit 2 (step S1 in FIG. 2).
  • the positive/negative inversion value calculator 3 calculates the time difference positive/negative inversion value Y(i) of the sampling data D(i) for each sampling time as shown in the following equation (step S2 in FIG. 2).
  • Y(i) - ⁇ D(i+1)-D(i-1) ⁇ (1)
  • the positive/negative inversion value calculation unit 3 inputs the calculated time difference positive/negative inversion value Y(i) to the FIFO buffer 50 at each sampling time (step S3 in FIG. 2).
  • the input value is held in the FIFO buffer 50, and after a time corresponding to the size of the FIFO buffer 50 (the delay time from when the positive/negative time difference value is input to the FIFO buffer 50 until it is output), It will be used for subtraction processing.
  • the positive/negative inversion value calculation unit 3 inputs the calculated time difference positive/negative inversion value Y(i) to the FIFO buffer 40 at each sampling time (step S4 in FIG. 2).
  • the output of FIFO buffer 40 is input to FIFO buffer 41 (step S5 in FIG. 2), and the output of FIFO buffer 41 is input to FIFO buffer 42 (step S6 in FIG. 2).
  • the FIFO buffers 40 to 42 are for obtaining the maximum value of the time difference positive/negative inversion value within a certain time range.
  • the time interval L3 corresponding to the length of the FIFO buffer 41 (the delay time from when the time difference positive/negative inversion value is input to the FIFO buffer 41 until it is output) is the width of the peak derived from the R wave (approximately 10 ms). ), preferably about 50 ms.
  • a time interval L2 corresponding to the length of the FIFO buffer 40 (a delay time from when the time difference positive/negative inversion value is input to the FIFO buffer 40 until it is output)
  • the maximum value M can be obtained for the range from (L3/2) to (L2+L3/2), and the maximum value M can be subtracted from the output value a.
  • the detection processing unit 43 detects the maximum value M of the time difference positive/negative inversion value stored in the FIFO buffer 40 and the time difference positive/negative inversion value stored in the FIFO buffer 42 at each sampling time (step S7 in FIG. 2). .
  • the threshold determination unit 6 compares the current threshold Th for heartbeat detection with the lower limit value Thmin of the plausible range as the threshold Th derived from the ECG waveform of the living body (step S9 in FIG. 2).
  • a likely range for the threshold Th derived from the ECG waveform of the living body can be determined based on past measurement results.
  • the initial value of the threshold Th may be set in advance to an arbitrary value within a plausible range as the threshold Th.
  • the heartbeat time determination unit 7 detects a peak of the index value b(i) exceeding the threshold value Th (Fig. 2: YES at step S10), the sampling time of this peak is set as the heartbeat time (step S11 in FIG. 2).
  • the sampling time of the index value b(i) refers to the sampling time of the time difference positive/negative inversion value Y(i) on which the index value b(i) is based (the sampling time of the data D(i)).
  • the sampling time of the index value b(i) can be acquired from the storage unit 2 .
  • the heartbeat time determining unit 7 detects a peak of the index value b(i) exceeding the threshold Th when the threshold determining unit 6 determines that the current threshold Th is less than the lower limit Thmin (NO in step S9). If not (NO in step S10), the heartbeat (R wave) is not determined and the heartbeat time is not determined.
  • FIG. 3 is a flow chart for explaining the operation of the threshold setting unit 8 (step S12).
  • the threshold setting unit 8 compares the calculated threshold candidate Thc with a threshold limit value L based on the current threshold Th and the difference limit value that can be said to be plausible as the time difference positive/negative inversion value of the ECG waveform of the living body.
  • the threshold Th is not updated (YES in step S23 of FIG. 3).
  • threshold value Th (step S24 in FIG. 3).
  • the threshold limit value L will be described later.
  • the threshold setting unit 8 performs the above processing at each sampling time.
  • the process of step S22 is repeatedly executed to sequentially calculate the threshold candidate Thc, and the threshold candidate Thc becomes the threshold value. If it is equal to or less than the limit value L, the threshold Th is updated.
  • the QRS interval (from the beginning of the Q wave to the end of the S wave) is 0.06-0.1 s.
  • the time from the maximum value of the R wave to the minimum value of the S wave can be considered to be about 15 to 25 ms as 1/4 of the QRS interval. Furthermore, the sharpest change occurs near the center of the range from the highest value of the R wave to the lowest value of the S wave. Assuming that the time range is about half from the maximum value of the R wave to the minimum value of the S wave, it is estimated to be about 7.5 to 12.5 ms. If there is a change in the electrocardiographic potential corresponding to the QRS amplitude during this time, the rate of change in the electrocardiographic potential is considered to be about 350 ⁇ V/ms at maximum. If the change rate of the electrocardiographic potential is converted into the positive/negative inversion value of the time difference obtained by the equation (1), it becomes 700 ⁇ V when the sampling interval is 1 ms, and is further rounded to 1000 ⁇ V with a margin.
  • the index value b represents the height of the clearance of the time difference value in the surrounding time domain with respect to the peak derived from the R wave.
  • the noise presents a sharp unimodal peak, in the worst case the value is reflected in the index value b as it is.
  • the threshold Th is set based on the value obtained by multiplying the value of the peak to be detected by a coefficient ⁇ of 1 or less.
  • the threshold Th decreases.
  • the threshold Th is unbelievably low as a threshold Th derived from a human ECG waveform, even if there is a peak in the index value b, the probability that this peak is derived from the R wave is low.
  • the heartbeat determination (step S10) is not performed. detection can be achieved.
  • the threshold update method described with reference to FIG. 3 is an example, and the threshold may be updated by another method.
  • the storage unit 2, the positive/negative inversion value calculation unit 3, the maximum value detection unit 4, the index value calculation unit 5, the threshold determination unit 6, the heartbeat time determination unit 7, and the threshold value setting unit 8 of the heartbeat detection apparatus described in this embodiment are , a CPU (Central Processing Unit), a storage device, and an interface, and a program that controls these hardware resources.
  • a CPU Central Processing Unit
  • FIG. 1 A configuration example of this computer is shown in FIG.
  • the computer comprises a CPU 200 , a storage device 201 and an interface device (I/F) 202 .
  • the I/F 202 is connected to the electrocardiograph 1 and the like.
  • a program for implementing the heartbeat detection method of the present invention is stored in the storage device 201 .
  • the CPU 200 executes the processing described in this embodiment according to the programs stored in the storage device 201 .
  • the present invention can be applied to technology for detecting the heartbeat of a living body.
  • SYMBOLS 1 Electrocardiograph, 2... Storage part, 3... Positive/negative inversion value calculation part, 4... Maximum value detection part, 5... Index value calculation part, 6... Threshold determination part, 7... Heartbeat time determination part, 8... Threshold value setting Sections 40 to 42, 50... FIFO buffer, 43... Detection processing section, 51... Subtraction processing section.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Cardiology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

L'invention concerne un dispositif de détection de battement cardiaque comprenant : une unité de calcul de valeur d'indice (5) qui calcule, pour chaque temps d'échantillonnage, des valeurs d'indice à partir d'une chaîne de données d'échantillonnage d'une forme d'onde d'ECG d'un corps vivant ; une unité de détermination de valeur seuil (6) qui compare une valeur seuil pour une détection de battement cardiaque à la limite inférieure de la plage qui est susceptible d'être une valeur seuil dérivée de la forme d'onde d'ECG du corps vivant ; une unité de détermination de temps de battement cardiaque (7) qui détermine le temps d'échantillonnage du pic en tant que temps de battement cardiaque lorsque l'unité de détermination de valeur seuil (6) détermine que la valeur seuil n'est pas inférieure à la limite inférieure et détecte un pic de valeur d'indice dépassant la valeur seuil ; et une unité de réglage de valeur seuil (8) qui renouvelle la valeur seuil sur la base du pic de valeur d'indice dépassant la valeur seuil actuelle.
PCT/JP2021/014756 2021-04-07 2021-04-07 Procédé de détection de battement cardiaque et dispositif de détection de battement cardiaque WO2022215191A1 (fr)

Priority Applications (2)

Application Number Priority Date Filing Date Title
PCT/JP2021/014756 WO2022215191A1 (fr) 2021-04-07 2021-04-07 Procédé de détection de battement cardiaque et dispositif de détection de battement cardiaque
JP2023512570A JP7491466B2 (ja) 2021-04-07 2021-04-07 心拍検出方法および心拍検出装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2021/014756 WO2022215191A1 (fr) 2021-04-07 2021-04-07 Procédé de détection de battement cardiaque et dispositif de détection de battement cardiaque

Publications (1)

Publication Number Publication Date
WO2022215191A1 true WO2022215191A1 (fr) 2022-10-13

Family

ID=83545295

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2021/014756 WO2022215191A1 (fr) 2021-04-07 2021-04-07 Procédé de détection de battement cardiaque et dispositif de détection de battement cardiaque

Country Status (2)

Country Link
JP (1) JP7491466B2 (fr)
WO (1) WO2022215191A1 (fr)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017150156A1 (fr) * 2016-02-29 2017-09-08 日本電信電話株式会社 Procédé de détection des battements cardiaques et dispositif de détection des battements cardiaques
WO2018074145A1 (fr) * 2016-10-17 2018-04-26 日本電信電話株式会社 Procédé de détection de battement cardiaque et dispositif de détection de battement cardiaque

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017150156A1 (fr) * 2016-02-29 2017-09-08 日本電信電話株式会社 Procédé de détection des battements cardiaques et dispositif de détection des battements cardiaques
WO2018074145A1 (fr) * 2016-10-17 2018-04-26 日本電信電話株式会社 Procédé de détection de battement cardiaque et dispositif de détection de battement cardiaque

Also Published As

Publication number Publication date
JPWO2022215191A1 (fr) 2022-10-13
JP7491466B2 (ja) 2024-05-28

Similar Documents

Publication Publication Date Title
JP6652655B2 (ja) 心拍検出方法および心拍検出装置
US10945623B2 (en) Heartbeat detection method and heartbeat detection device
US10172531B2 (en) Heartbeat detection method and heartbeat detection device
JP6645926B2 (ja) 生体信号処理方法および装置
JP6404784B2 (ja) 心拍検出方法および心拍検出装置
CN112773582A (zh) 一种假肢控制方法、装置、系统及介质
JP2016123594A (ja) 自律神経活動推定装置及び自律神経活動推定方法
JP2018175152A (ja) 心拍検出方法および心拍検出装置
WO2022215191A1 (fr) Procédé de détection de battement cardiaque et dispositif de détection de battement cardiaque
US10750969B2 (en) Heartbeat detection method and heartbeat detection device
JP2015217060A (ja) 心拍検出方法および心拍検出装置
WO2023032119A1 (fr) Procédé de détection de battement cardiaque et dispositif de détection de battement cardiaque
JPWO2014141682A1 (ja) センサー端末および信号取得方法
WO2022215190A1 (fr) Procédé de détection de battement cardiaque et dispositif de détection de battement cardiaque
WO2023032120A1 (fr) Procédé de détection de battement cardiaque et dispositif de détection de battement cardiaque
EP2886044A1 (fr) Système et procédé pour l'acquisition de signaux biopotentiels avec réduction d'artéfacts de mouvement
CN106037718B (zh) 可穿戴心电图系统
JP7147866B2 (ja) 心拍検出方法、心拍検出装置およびプログラム
WO2022149382A1 (fr) Dispositif d'analyse de signal, procédé d'analyse de signal, et programme
WO2022149307A1 (fr) Dispositif d'analyse de signal, procédé d'analyse de signal et programme
JP7248150B2 (ja) 心拍検出装置、心拍検出方法、および心拍検出プログラム
JP2020096710A (ja) Ｅｃｇピーク検出装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 21935995

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2023512570

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 21935995

Country of ref document: EP

Kind code of ref document: A1