[go: up one dir, main page]

WO2021260783A1 - 心拍数検知方法、装置およびプログラム - Google Patents

心拍数検知方法、装置およびプログラム Download PDF

Info

Publication number
WO2021260783A1
WO2021260783A1 PCT/JP2020/024528 JP2020024528W WO2021260783A1 WO 2021260783 A1 WO2021260783 A1 WO 2021260783A1 JP 2020024528 W JP2020024528 W JP 2020024528W WO 2021260783 A1 WO2021260783 A1 WO 2021260783A1
Authority
WO
WIPO (PCT)
Prior art keywords
heart rate
index
interval
instantaneous
reference value
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.)
Ceased
Application number
PCT/JP2020/024528
Other languages
English (en)
French (fr)
Inventor
伸昭 松浦
啓 桑原
雄一 樋口
隆子 石原
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone Corp
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 Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to US18/007,939 priority Critical patent/US20230233130A1/en
Priority to JP2022531269A priority patent/JP7439924B2/ja
Priority to PCT/JP2020/024528 priority patent/WO2021260783A1/ja
Publication of WO2021260783A1 publication Critical patent/WO2021260783A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient; User input means
    • A61B5/742—Details of notification to user or communication with user or patient; User input means using visual displays

Definitions

  • the present invention relates to a heart rate detection method, device and program that presents an accurate heart rate.
  • Heart rate or its fluctuation is biometric information obtained from ECG (Electrocardiography), and is used for evaluation of autonomic nervous function and index of exercise intensity in daily life and resting state.
  • measuring the heart rate over a relatively long time such as daily life using a wearable ECG waveform acquisition / heart rate measurement device is useful for managing the physical condition of the user.
  • Patent Document 1 discloses a method capable of appropriately detecting a heartbeat when noise is superimposed on an ECG waveform, such as measurement by a wearable device.
  • the heartbeat detection may be erroneous. For example, when an R wave to be detected is missed, a wave that is not an R wave may be detected as an R wave.
  • the instantaneous heart rate caused by the error in heart rate detection is performed by moving averaging the heart rate over a certain period of time or using an IIR filter.
  • a method of calculating and presenting an average heart rate that suppresses fluctuations in the heart rate may be adopted (Patent Document 2).
  • the heart rate detection method has a step of measuring an ECG waveform, a step of calculating an RR interval from the ECG waveform, and an instantaneous step from the RR interval.
  • the heart rate detection device calculates the RR interval from the measurement unit for measuring the ECG waveform and the ECG waveform, calculates the instantaneous heart rate from the RR interval, and calculates the instantaneous heart rate. It is provided with a calculation unit for calculating an index based on the above, and a determination unit for determining whether or not to display the heart rate obtained from the RR interval by comparing the index with a reference value.
  • the heart rate detection program includes a step of measuring an ECG waveform, a step of calculating an RR interval from the ECG waveform, a step of calculating an instantaneous heart rate from the RR interval, and the above-mentioned step.
  • a process including a step of calculating an index based on an instantaneous heart rate and a step of comparing the index with a reference value and determining whether or not to display a heart rate obtained from the RR interval is executed. It features a heart rate detector to function.
  • the present invention it is possible to provide a heart rate detection method, a device and a program for presenting an accurate heart rate.
  • FIG. 1 is a diagram showing an ECG waveform and an RR interval for explaining the heart rate detection method according to the first embodiment of the present invention.
  • FIG. 2 is a diagram showing an instantaneous heart rate and an average heart rate for explaining a heart rate detection method according to the first embodiment of the present invention.
  • FIG. 3 is a diagram showing changes over time of an index in the heart rate detection method according to the first embodiment of the present invention.
  • FIG. 4 is a block diagram of the heart rate detecting device according to the first embodiment of the present invention.
  • FIG. 5 is a flowchart showing a heart rate detection method according to the first embodiment of the present invention.
  • FIG. 6 is a diagram showing a configuration example of a computer according to an embodiment of the present invention.
  • FIG. 1 is a part of an ECG waveform measured using a wearable device.
  • FIG. 1 also shows the RR interval [ms] corresponding to the heart rate at the timing when the R wave is detected (black diamond in the figure).
  • the method disclosed in Patent Document 1 can be used for detecting the RR interval.
  • HR Heart rate
  • FIG. 2 shows an instantaneous heart rate (hereinafter referred to as “IHR” (Instantaneous Heart rate)) and an average heart rate HR calculated based on the instantaneous heart rate (hereinafter referred to as “IHR”).
  • IHR Instantaneous Heart rate
  • IHR Average heart rate
  • the average heart rate HR [n] is calculated by the following formula with respect to the time series data IHR [n] of the instantaneous heart rate.
  • the average heart rate shows a stable value, it may be unreliable.
  • the average heart rate is used as an index for accurate heart rate detection, the heart rate cannot always be detected accurately.
  • the index based on the instantaneous heart rate is the difference between the time-series data IHR [n] of the instantaneous heart rate, and the square of (IHR [n] -IHR [n-1]) is the latest 10 beats in the past. Calculated by adding (heart rate) minutes.
  • IHR [n] is an instantaneous heart rate at an arbitrary time t [n].
  • IHR [n-1] is the instantaneous heart rate at the latest time t [n-1] before t [n].
  • FIG. 3 shows the time course of the index based on the above-mentioned instantaneous heart rate for the waveform shown in FIG.
  • the index shows a high value at 15:53:20 to 15:53:25, decreases at 15:53:25 to 15:53:30, and a low value at 15:53:30 to 15:53:35. show.
  • the index exceeds the standard value at 15:53:20 to 15:53:30 and does not meet the standard. Therefore, it is determined that the heart rate measured from 15:53:30 to 15:53:35 is not accurate.
  • the index is below the standard value and meets the standard. Therefore, it is determined that the heart rate measured from 15:53:30 to 15:53:35 is accurate.
  • the index shows a value exceeding the reference value when the ECG waveform noise shown in FIG. 1 is large and accurate heart rate detection cannot be performed.
  • the index indicates a value equal to or less than the reference value.
  • an index based on the instantaneous heart rate it can be determined whether or not the measured heart rate is accurate, in other words, the heart rate can be accurately detected (presented).
  • the present invention is based on the finding that the heart rate can be accurately detected (presented) by an index based on the instantaneous heart rate.
  • FIG. 4 shows a block diagram of an example of the configuration of the heart rate detection device 1 according to the present embodiment.
  • the heart rate detection device 1 includes a measurement unit 11, a storage unit 12, a calculation unit 13, a determination unit 14, and an output unit 15.
  • the measuring unit 11 measures the ECG waveform, and as an example, an electrocardiograph in a wearable device is used.
  • the storage unit 12 stores the instantaneous heart rate calculated based on the ECG potential measured by the electrocardiograph. It is also possible to store the time measured by the electrocardiograph and the ECG potential.
  • the calculation unit 13 acquires the time and ECG potential measured by the electrocardiograph from the storage unit 12, and performs calculation based on the ECG waveform in order to detect the RR interval. In addition, the instantaneous heart rate is calculated to calculate a numerical value as an index.
  • the determination unit 14 compares the index with the reference value and determines whether or not the detected heart rate is accurate (validity).
  • the output unit 15 outputs the calculated heart rate. In addition, it is possible to output that the heart rate was not accurately measured, the measured ECG waveform, the time difference value, and the integrated value.
  • FIG. 5 shows a flowchart for explaining the heart rate detection method according to the present embodiment.
  • step 21 measure the electrocardiogram (ECG waveform) (step 21).
  • the RR interval is detected using the ECG waveform data (step 22).
  • the RR interval can be calculated (detected), for example, by comparing the ECG waveform with a predetermined threshold value based on the change in the time-difference value. Further, for example, the detection accuracy can be improved by comparing a plurality of ECG waveforms with a predetermined threshold value by using a change in a value obtained by time-difference (see Patent Document 1).
  • the instantaneous heart rate is acquired as time-series data IHR [n] from the RR interval detected at an arbitrary time t [n].
  • step 24 the time when the RR interval is detected and the instantaneous heart rate IHR calculated from the RR interval are stored (step 24).
  • IHR [n] is an instantaneous heart rate at an arbitrary time t [n].
  • IHR [n-1] is an instantaneous heart rate at the latest time t [n-1] before t [n], and the instantaneous heart rate stored in the storage unit is read out and used.
  • IHR [n] -IHR [n-1] is a difference value of the instantaneous heart rate IHR between the sampling times.
  • step 26 the calculated index is compared with the reference value, and the determination is made as follows (step 26).
  • the index is less than or equal to the reference value, it is determined that the heart rate calculated from the ECG waveform is accurate.
  • the heart rate obtained from the RR interval detected at time t [n] is output (presented) (step 27).
  • the index exceeds the reference value, it is determined that the heart rate calculated from the ECG waveform is not accurate.
  • the heart rate is not calculated, and it is output (displayed) that the accurate heart rate is not measured (step 28).
  • the most recent accurate heart rate before the time t [n] may be displayed, or a blank (a state in which nothing is displayed) may be displayed.
  • the reference value is set to 20,000 [bpm 2 ].
  • the index shown in FIG. 3 is calculated.
  • the reference value exceeds 20,000 [bpm 2 ], so it is determined that the heart rate calculated from the ECG waveform is not accurate. .. Therefore, it is output (presented) that the accurate heart rate is not measured.
  • the most recent accurate heart rate before the measurement time is presented.
  • a blank nothing is displayed
  • the present embodiment it is possible to avoid presenting an abnormal value to the user even though the heart rate is not measured as an accurate value, and it is accurate. Only heart rate numbers can be presented to the user.
  • the reference value of 20,000 [bpm 2 ] is used, but the present invention is not limited to this, and other values may be used.
  • the index is (IHR [n] -IHR [ The square of n-1]) was used, but the present invention is not limited to this.
  • the difference between the instantaneous heart rate and the average heart rate may be used.
  • a step of calculating the average heart rate after calculating the instantaneous heart rate is required.
  • or (IHR [n] -IHR [n-1]) may be used as the index.
  • the absolute value of the instantaneous heart rate may be used.
  • the present embodiment it is possible to avoid presenting an abnormal value to the user even though the heart rate is not measured as an accurate value, and it is accurate. Only heart rate numbers can be presented to the user.
  • the heart rate detection device may be worn on the user's body as an integrated wearable device.
  • the heart rate detection device may be equipped with a storage unit, a calculation unit, and a determination unit on a smartphone, a server, or the like outside the wearable device by attaching the measurement unit to the user's body as a wearable device. ..
  • the heart rate detection device is provided with a transmission / reception unit for each of the wearable device and an external server, etc., and the ECG waveform measured by the wearable device is transmitted to the server or the like, and is stored, calculated, determined, or the like by the server or the like.
  • the heart rate and the like may be output to a server or the like, or may be transmitted to a wearable device or the like and output.
  • ECG potential ECG potential
  • the ECG waveforms may be acquired from the storage unit to determine the validity of the heart rate.
  • FIG. 6 shows a configuration example of the computer 30 in the heart rate detection device according to the embodiment of the present invention.
  • the heart rate detection device can be realized by a computer 30 including a CPU (Central Processing Unit) 33, a storage device (storage unit) 32, and an interface device 31, and a program for controlling these hardware resources.
  • the interface device 31 is connected to the measurement unit and the output unit of the heart rate detection device according to the embodiment of the present invention.
  • the CPU 33 executes the process according to the embodiment of the present invention according to the heart rate detection program stored in the storage device 32. In this way, the heart rate detection program activates the heart rate detection device.
  • a computer may be provided inside the device, or at least one part of the functions of the computer may be realized by using an external computer.
  • the storage unit may also use a storage medium outside the device, or may read out and execute a heart rate detection program stored in the storage medium.
  • the storage medium includes various magnetic recording media, optical magnetic recording media, CD-ROMs, CD-Rs, and various memories.
  • the heart rate detection program may be supplied to the computer via a communication line such as the Internet.
  • the present invention can be applied to a technique for analyzing a biological signal obtained from an ECG waveform.

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)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本発明の心拍数検知方法は、ＥＣＧ波形を測定するステップと、ＥＣＧ波形よりＲ－Ｒ間隔を算出するステップと、Ｒ－Ｒ間隔より瞬時心拍数を算出するステップと、瞬時心拍数に基づく指標を算出するステップと、指標を基準値と比較して、Ｒ－Ｒ間隔より求まる心拍数を表示するか否かを判定するステップとを備える。　これにより、本発明は、正確な心拍数を提示する心拍数検知方法を提供することができる。

Description

心拍数検知方法、装置およびプログラム
　本発明は、正確な心拍数を提示する心拍数検知方法、装置およびプログラムに関する。
　心拍数あるいはその変動はＥＣＧ（Ｅｌｅｃｔｒｏｃａｒｄｉｏｇｒａｍ、心電図）から得られる生体情報であり、日常生活や安静状態においては自律神経機能の評価、運動強度の指標などに活用される。
　とくに、ウェアラブルなＥＣＧ波形取得・心拍計測デバイスを用いて、日常生活等の比較的長い時間にわたって心拍数を計測することは、ユーザの体調管理等に有用である。
　特許文献１には、ウェアラブルデバイスによる計測のように、ＥＣＧ波形にノイズが重畳している場合に適切に心拍を検出できる方法が開示されている。
　しかしながら、この方法において、体動によるＥＣＧ波形のノイズが激しい場合には、心拍検出を誤ることがある。例えば、検出すべきＲ波を取りこぼす場合、Ｒ波でないものをＲ波として検出してしまう場合などがある。
　ウェアラブル心拍計測デバイスでは、このような心拍検出の誤りを回避するために、心拍数を一定の時間にわたって移動平均することや、ＩＩＲフィルタを用いることなどにより、心拍検出の誤りに起因する瞬時心拍数の変動を抑えた、平均心拍数を算出・提示する手法が採られることがある（特許文献２）。
特許第６５２７２８６号公報 特許第６６４５９２６号公報
　しかしながら、特許文献２に開示された心拍の検出方法においては、どのようなＥＣＧ波形データに対しても、平均心拍数等の値が算出される。したがって、心拍検出に誤りが含まれ、正確な値として計測されていない状況にもかかわらず、異常な数値をユーザに正確な値として提示してしまうことがある。
　上述したような課題を解決するために、本発明に係る心拍数検知方法は、ＥＣＧ波形を測定するステップと、前記ＥＣＧ波形よりＲ－Ｒ間隔を算出するステップと、前記Ｒ－Ｒ間隔より瞬時心拍数を算出するステップと、前記瞬時心拍数に基づく指標を算出するステップと、前記指標を基準値と比較して、前記Ｒ－Ｒ間隔より求まる心拍数を表示するか否かを判定するステップとを備える。
　また、本発明に係る心拍数検知装置は、ＥＣＧ波形を測定する測定部と、前記ＥＣＧ波形よりＲ－Ｒ間隔を算出し、前記Ｒ－Ｒ間隔より瞬時心拍数を算出し、前記瞬時心拍数に基づく指標を算出する算出部と、前記指標を基準値と比較して前記Ｒ－Ｒ間隔より求まる心拍数を表示するか否かを判定する判定部とを備える。
　また、本発明に係る心拍数検知プログラムは、ＥＣＧ波形を測定するステップと、前記ＥＣＧ波形よりＲ－Ｒ間隔を算出するステップと、前記Ｒ－Ｒ間隔より瞬時心拍数を算出するステップと、前記瞬時心拍数に基づく指標を算出するステップと、前記指標を基準値と比較して、前記Ｒ－Ｒ間隔より求まる心拍数を表示するか否かを判定するステップとを備える処理を実行させることを特徴とし、心拍数検知装置を機能させる。
　本発明によれば、正確な心拍数を提示する心拍数検知方法、装置およびプログラムを提供できる。
図１は、本発明の第１の実施の形態に係る心拍数検知方法を説明するためのＥＣＧ波形とＲ－Ｒ間隔を示す図である。 図２は、本発明の第１の実施の形態に係る心拍数検知方法を説明するための瞬時心拍数と平均心拍数を示す図である。 図３は、本発明の第１の実施の形態に係る心拍数検知方法における指標の経時変化を示す図である。 図４は、本発明の第１の実施の形態に係る心拍数検知装置のブロック図である。 図５は、本発明の第１の実施の形態に係る心拍数検知方法を示すフローチャート図である。 図６は、本発明の実施の形態におけるコンピュータの構成例を示す図である。
＜第１の実施の形態＞
　本発明の第１の実施の形態について図１～図５を参照して説明する。
　図１は、ウェアラブルなデバイスを用いて計測されたＥＣＧ波形の一部である。図１には、Ｒ波が検出されたタイミングに、心拍数に相当するＲ－Ｒ間隔[ｍｓ]も併せて示す（図中、黒ひし形）。ここで、例えば、Ｒ－Ｒ間隔の検出には、特許文献１に開示された方法を用いることができる。
　この波形では、１５：５３：３０～１５：５３：３３では、ＥＣＧ波形のノイズが小さくＲ波はほぼ適切に検出されるが、１５：５３：２０～１５：５３：２５では、ＥＣＧ波形のノイズが大きく、Ｒ波が的確に検出されていないことがわかる。
　次に、従来法の１つである、正確な心拍数検知のための指標に平均心拍数（以下、「ＨＲ」（Ｈｅａｒｔ　ｒａｔｅ）という。）を用いる手法を説明する。
　図２に、瞬時心拍数（以下、「ＩＨＲ」（Ｉｎｓｔａｎｔａｎｅｏｕｓ　Ｈｅａｒｔ　ｒａｔｅ）という。）と、それを基に算出された平均心拍数ＨＲを示す。図２において、黒ひし形は瞬時心拍数を示し、白丸は平均心拍数を示す。
　瞬時心拍数ＩＨＲは、瞬時心拍数ＩＨＲ[ｂｐｍ]＝６０，０００÷（Ｒ－Ｒ間隔[ｍｓ]）より算出される。
　また、ＩＩＲフィルタを用いて、瞬時心拍数の時系列データＩＨＲ[ｎ]に対し、平均心拍数ＨＲ[ｎ]は、以下の式で算出される。
　ＨＲ[ｎ]＝（１－ａ）×ＨＲ[ｎ－１]＋ａ×ＩＨＲ[ｎ]　
　ここで、ａ＝０．１である。
　図２に示すように、１５：５３：３０～１５：５３：３３では、瞬時心拍数のばらつきが少なく、平均心拍数も安定して観測され、信頼性の高いデータが得られる。
　しかしながら、１５：５３：２０～１５：５３：２５では、瞬時心拍数のばらつきは大きい。一方、平均心拍数はほぼ一定の値で観測される。この観測された平均心拍数は、ばらつきの大きい瞬時心拍数に基づき得られているので信頼性に欠ける。
　このことは、平均心拍数は安定した値を示しても、信頼性に欠ける可能性があることを示唆している。このように、平均心拍数を正確な心拍数検知のための指標としても、必ずしも正確に心拍数を検知できない。
　そこで、瞬時心拍数のばらつきに着目して、瞬時心拍数に基づく指標について検討する。ここでは、瞬時心拍数に基づく指標を、瞬時心拍数の時系列データＩＨＲ[ｎ]の差分に関するものとして、（ＩＨＲ[ｎ]―ＩＨＲ[ｎ－１]）の２乗を過去直近の１０拍（心拍数）分加算して算出する。ここで、ＩＨＲ[ｎ]は任意の時刻ｔ［ｎ］における瞬時心拍数である。ＩＨＲ[ｎ－１]はｔ［ｎ］以前の直近の時刻ｔ［ｎ－１］における瞬時心拍数である。
　図３に、図１に示す波形について、上述の瞬時心拍数に基づく指標の経時変化を示す。指標は、１５：５３：２０～１５：５３：２５で高い値を示し、１５：５３：２５～１５：５３：３０で減少し、１５：５３：３０～１５：５３：３５で低い値を示す。
　この指標の経時変化において、基準値を２０，０００[ｂｐｍ２]とすれば、１５：５３：２０～１５：５３：３０で指標は基準値を上回り、基準を満たさない。したがって、１５：５３：３０～１５：５３：３５で測定される心拍数が正確でないと判定される。
　一方、１５：５３：３０～１５：５３：３５で指標は基準値以下となり、基準を満たす。したがって、１５：５３：３０～１５：５３：３５で測定される心拍数が正確であると判定される。
　このように、瞬時心拍数に基づく指標を用いれば、図１に示すＥＣＧ波形ノイズが大きく正確な心拍数検知ができないときには、指標は基準値を上回る値を示す。一方、ＥＣＧ波形ノイズが小さく正確な心拍数検知ができるときには、指標は基準値以下の値を示す。
　したがって、瞬時心拍数に基づく指標を用いれば、測定された心拍数が正確か否かを判定でき、換言すれば、正確に心拍数を検知（提示）できる。
　本発明は、瞬時心拍数に基づく指標により正確に心拍数を検知（提示）できるという知見に基づくものである。
＜心拍数検知装置の構成＞
　図４に、本実施の形態に係る心拍数検知装置１の構成の一例のブロック図を示す。心拍数検知装置１は測定部１１と、記憶部１２と、算出部１３と、判定部１４と、出力部１５を備える。
　測定部１１はＥＣＧ波形を測定し、一例として、ウェアラブルデバイスにおける心電計を用いる。
　記憶部１２は、心電計で測定したＥＣＧ電位に基づき算出する瞬時心拍数を記憶する。
また、心電計で測定した時刻とＥＣＧ電位を記憶することもできる。
　算出部１３では、記憶部１２から心電計で測定した時刻とＥＣＧ電位を取得して、Ｒ－Ｒ間隔を検知するためにＥＣＧ波形に基づき計算を行う。また、瞬時心拍数を算出して、指標となる数値を算出する。
　判定部１４は、指標と基準値を比較して、検知される心拍数が正確か否か（正当性）を判定する。
　出力部１５では、計算された心拍数を出力する。また、心拍数が正確に測定されなかったことや、計測されたＥＣＧ波形や時間差分値、積算値も出力できる。
＜心拍数検知方法＞
　図５に、本実施の形態に係る心拍数検知方法を説明するためのフローチャート図を示す。
　初めに、心電図（ＥＣＧ波形）を測定する（ステップ２１）。
　次に、ＥＣＧ波形のデータを用いて、Ｒ－Ｒ間隔を検出する（ステップ２２）。ここで、Ｒ－Ｒ間隔は、例えば、ＥＣＧ波形を時間差分した値の変化を基に所定の閾値と比較して算出（検出）することができる。また、例えば、複数のＥＣＧ波形を時間差分した値の変化を用いて所定の閾値と比較することにより検出精度を向上させることができる（特許文献１参照）
　次に、検出したＲ－Ｒ間隔から、瞬時心拍数ＩＨＲを、ＩＨＲ＝６０，０００÷（Ｒ－Ｒ間隔[ｍｓ]）より、算出する（ステップ２３）。ここで、瞬時心拍数は、任意の時刻ｔ［ｎ］に検出されるＲ－Ｒ間隔から、時系列データＩＨＲ[ｎ]として取得される。
　次に、Ｒ－Ｒ間隔が検出された時刻と、Ｒ－Ｒ間隔から算出された瞬時心拍数ＩＨＲを記憶する（ステップ２４）。
　次に、瞬時心拍数の時系列データＩＨＲ[ｎ]の差分に関する指標を、（ＩＨＲ[ｎ]―ＩＨＲ[ｎ－１]）の２乗を過去直近の１０拍（心拍数）分加算して算出する（ステップ２５）。ここで、ＩＨＲ[ｎ]は任意の時刻ｔ［ｎ］における瞬時心拍数である。ＩＨＲ[ｎ－１]はｔ［ｎ］以前の直近の時刻ｔ［ｎ－１］における瞬時心拍数であり、記憶部に記憶されている瞬時心拍数を読み出して用いる。このように、ＩＨＲ[ｎ]―ＩＨＲ[ｎ－１]は、サンプリング時刻間での瞬時心拍数ＩＨＲの差分値である。
　ここで、指標を算出する際に、過去直近の１０拍（心拍数）分加算する例を示したが、これに限らず、所定の拍数分を加算すればよい。
　次に、この算出された指標を基準値と比較して、以下の通り判定する（ステップ２６）。
　指標が基準値以下の場合には、ＥＣＧ波形から算出される心拍数は正確ではあると判定される。
　この場合には、時刻ｔ［ｎ］に検出されるＲ－Ｒ間隔から求まる心拍数を出力（提示）する（ステップ２７）。
　一方、指標が基準値を上回る場合には、ＥＣＧ波形から算出される心拍数は正確ではないと判定される。
　この場合は、心拍数を計算せず、正確な心拍数が測定されないことが出力（表示）される（ステップ２８）。または、時刻ｔ［ｎ］以前の直近の正確な心拍数を表示してもよいし、ブランク（何も表示されない状態）を表示してもよい。
　次に、本実施の形態に係る心拍数検知方法および装置における動作（出力）の一例を説明する。ここで、基準値を基準値を２０，０００[ｂｐｍ２]とする。
　図１に示すＥＣＧ波形が計測される場合には、図３に示す指標が算出される。次に、図３において時刻１５：５３：２０～１５：５３：３０では、基準値は２０，０００[ｂｐｍ２]を上回るので、ＥＣＧ波形から算出される心拍数は正確ではないと判定される。そこで、正確な心拍数が測定されないことが出力（提示）される。または、測定時刻以前の直近の正確な心拍数を提示される。または、ブランク（何も表示されない状態）が提示される。
　一方、図３において時刻１５：５３：３０～１５：５３：３５では、基準値は２０，０００[ｂｐｍ２]以下なので、ＥＣＧ波形から算出される心拍数は正確ではあると判定される。そこで、図１に示すＥＣＧ波形より算出された心拍数が出力（提示）される。
　以上のように、本実施の形態によれば、心拍数が正確な値として計測されていない状況にもかかわらず、異常な数値をユーザに提示してしまうことを回避することができ、正確な心拍数の数値のみをユーザに提示することができる。
　本実施の形態では、基準値の２０，０００[ｂｐｍ２]を用いたが、これに限らず、他の値でもよい。
　本実施の形態では、検知された心拍数データの信頼性が低い場合には、前後の瞬時心拍数の時系列データの差が大きいことに着目して、指標に（ＩＨＲ[ｎ]―ＩＨＲ[ｎ－１]）の２乗を用いたが、これに限らない。指標には、瞬時心拍数と平均心拍数の差を用いてもよい。この場合は、図５に示す第１の実施の形態に係る心拍数検知方法を示すフローチャート図において、瞬時心拍数を算出した後に、平均心拍数を算出するステップを必要とする。
　または、指標には、｜ＩＨＲ[ｎ]―ＩＨＲ[ｎ－１]｜や（ＩＨＲ[ｎ]―ＩＨＲ[ｎ－１]）のｎ乗値を用いてもよい。また、瞬時心拍数の絶対値を用いてもよい。
　このように、指標には、瞬時心拍数に基づく数値を用いることができる。
　以上のように、本実施の形態によれば、心拍数が正確な値として計測されていない状況にもかかわらず、異常な数値をユーザに提示してしまうことを回避することができ、正確な心拍数の数値のみをユーザに提示することができる。
　本実施の形態に係る心拍数検知装置は、ウェアラブルデバイス一体としてユーザの身体に装着してもよい。
　または、本実施の形態に係る心拍数検知装置は、測定部をウェアラブルデバイスとしてユーザの身体に装着して、ウェアラブルデバイス外部のスマートフォンやサーバ等に記憶部、算出部、判定部を備えてもよい。この場合、心拍数検知装置はウェアラブルデバイスと外部のサーバ等それぞれに送受信部を備え、ウェアラブルデバイスで測定されるＥＣＧ波形をサーバ等に送信し、サーバ等で記憶、算出、判定などを行う。最後に、心拍数等（心拍数は測定されないことの表示等を含む）はサーバ等に出力されてもよいし、ウェアラブルデバイス等に送信されて出力されてもよい。
　本発明に係る実施の形態では、測定されるタイミングで、順次記憶部より直近のＥＣＧ波形（ＥＣＧ電位）を取得して心拍数が正確か否かの正当性を判定する例を示したが、ＥＣＧ波形をまとめて記憶してから、記憶部よりＥＣＧ波形を取得して心拍数の正当性を判定してもよい。
　図６に、本発明の実施の形態に係る心拍数検知装置におけるコンピュータ３０の構成例を示す。心拍数検知装置は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）３３、記憶装置（記憶部）３２およびインタフェース装置３１を備えたコンピュータ３０と、これらのハードウェア資源を制御するプログラムによって実現することができる。ここで、インタフェース装置３１に、本発明の実施の形態に係る心拍数検知装置における測定部と、出力部が接続される。ＣＰＵ３３は、記憶装置３２に格納された心拍数検知プログラムに従って本発明の実施の形態における処理を実行する。このように、心拍数検知プログラムは心拍数検知装置を機能させる。
　本発明の実施の形態に係る心拍数検知装置では、コンピュータを装置内部に備えてもよいし、コンピュータの機能の少なくとも１部を外部コンピュータを用いて実現してもよい。また、記憶部も装置外部の記憶媒体を用いてもよく、記憶媒体に格納された心拍数検知プログラムを読み出して実行してもよい。記憶媒体には、各種磁気記録媒体、光磁気記録媒体、ＣＤ－ＲＯＭ、ＣＤ－Ｒ、各種メモリを含む。また、心拍数検知プログラムはインターネットなどの通信回線を介してコンピュータに供給されてもよい。
　本発明の実施の形態では、心拍数検知装置の構成、心拍数検知方法などにおいて、各構成部の構造、寸法、材料等の一例を示したが、これに限らない。心拍数検知装置および方法の機能を発揮し効果を奏するものであればよい。
  本発明は、ＥＣＧ波形波形から得られる生体信号を分析する技術に適用することができる。
１　心拍数検知装置
１１　測定部
１２　記憶部
１３　算出部
１４　判定部

Claims (8)

  1. 　ＥＣＧ波形を測定するステップと、
    　前記ＥＣＧ波形よりＲ－Ｒ間隔を算出するステップと、
    　前記Ｒ－Ｒ間隔より瞬時心拍数を算出するステップと、
    　前記瞬時心拍数に基づく指標を算出するステップと、
    　前記指標を基準値と比較して、前記Ｒ－Ｒ間隔より求まる心拍数を表示するか否かを判定するステップと
    を備える心拍数検知方法。
  2. 　前記指標がサンプリング時刻間での前記瞬時心拍数の差分値に基づくことを特徴とする請求項１に記載の心拍数検知方法。
  3. 　前記瞬時心拍数より平均心拍数を算出するステップを備え、
    　前記指標が前記瞬時心拍数と前記平均心拍数の差分値に基づくことを特徴とする請求項１に記載の心拍数検知方法。
  4. 　前記指標が基準値以下の場合に、前記Ｒ－Ｒ間隔より求まる心拍数を表示して、前記指標が基準値より大きい場合に、ブランクを表示することを特徴とする請求項１から請求項３のいずれか一項に記載の心拍数検知方法。
  5. 　前記指標が基準値以下の場合に、前記Ｒ－Ｒ間隔より求まる心拍数を表示して、前記指標が基準値より大きい場合に、正確な心拍数が測定されなかったことを表示することを特徴とする請求項１から請求項４のいずれか一項に記載の心拍数検知方法。
  6. 　前記指標が基準値以下の場合に、前記Ｒ－Ｒ間隔より求まる心拍数を表示して、前記指標が基準値より大きい場合に、前記Ｒ－Ｒ間隔が検出された時刻以前の直近の正確な心拍数を表示することを特徴とする請求項１から請求項５のいずれか一項に記載の心拍数検知方法。
  7. 　ＥＣＧ波形を測定する測定部と、
    　前記ＥＣＧ波形よりＲ－Ｒ間隔を算出し、前記Ｒ－Ｒ間隔より瞬時心拍数を算出し、前記瞬時心拍数に基づく指標を算出する算出部と、
    　前記指標を基準値と比較して前記Ｒ－Ｒ間隔より求まる心拍数を表示するか否かを判定する判定部と
    を備える心拍数検知装置。
  8. 　ＥＣＧ波形を測定するステップと、
    　前記ＥＣＧ波形よりＲ－Ｒ間隔を算出するステップと、
    　前記Ｒ－Ｒ間隔より瞬時心拍数を算出するステップと、
    　前記瞬時心拍数に基づく指標を算出するステップと、
    　前記指標を基準値と比較して、前記Ｒ－Ｒ間隔より求まる心拍数を表示するか否かを判定するステップと
    を備える処理を実行させることを特徴とする、心拍数検知装置を機能させるための心拍数検知プログラム。
PCT/JP2020/024528 2020-06-23 2020-06-23 心拍数検知方法、装置およびプログラム Ceased WO2021260783A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US18/007,939 US20230233130A1 (en) 2020-06-23 2020-06-23 Heart Rate Detection Method, Device, and Program
JP2022531269A JP7439924B2 (ja) 2020-06-23 2020-06-23 心拍数検知方法、装置およびプログラム
PCT/JP2020/024528 WO2021260783A1 (ja) 2020-06-23 2020-06-23 心拍数検知方法、装置およびプログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2020/024528 WO2021260783A1 (ja) 2020-06-23 2020-06-23 心拍数検知方法、装置およびプログラム

Publications (1)

Publication Number Publication Date
WO2021260783A1 true WO2021260783A1 (ja) 2021-12-30

Family

ID=79282240

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2020/024528 Ceased WO2021260783A1 (ja) 2020-06-23 2020-06-23 心拍数検知方法、装置およびプログラム

Country Status (3)

Country Link
US (1) US20230233130A1 (ja)
JP (1) JP7439924B2 (ja)
WO (1) WO2021260783A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115444387A (zh) * 2022-09-02 2022-12-09 武汉中旗生物医疗电子有限公司 心率检测装置、设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006271731A (ja) * 2005-03-29 2006-10-12 Toshiba Corp 心拍計測装置及び心拍計測方法
WO2016024495A1 (ja) * 2014-08-11 2016-02-18 日本電信電話株式会社 生体信号測定システム、生体情報測定装置および生体情報抽出アルゴリズム変更方法
WO2017199597A1 (ja) * 2016-05-20 2017-11-23 ソニー株式会社 生体情報処理装置、生体情報処理方法、及び情報処理装置
JP2019097831A (ja) * 2017-11-30 2019-06-24 パラマウントベッド株式会社 異常判定装置及びそれに用いるプログラム
WO2019172111A1 (ja) * 2018-03-06 2019-09-12 日本電信電話株式会社 拍動数算出装置および方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006271731A (ja) * 2005-03-29 2006-10-12 Toshiba Corp 心拍計測装置及び心拍計測方法
WO2016024495A1 (ja) * 2014-08-11 2016-02-18 日本電信電話株式会社 生体信号測定システム、生体情報測定装置および生体情報抽出アルゴリズム変更方法
WO2017199597A1 (ja) * 2016-05-20 2017-11-23 ソニー株式会社 生体情報処理装置、生体情報処理方法、及び情報処理装置
JP2019097831A (ja) * 2017-11-30 2019-06-24 パラマウントベッド株式会社 異常判定装置及びそれに用いるプログラム
WO2019172111A1 (ja) * 2018-03-06 2019-09-12 日本電信電話株式会社 拍動数算出装置および方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115444387A (zh) * 2022-09-02 2022-12-09 武汉中旗生物医疗电子有限公司 心率检测装置、设备及存储介质

Also Published As

Publication number Publication date
JP7439924B2 (ja) 2024-02-28
US20230233130A1 (en) 2023-07-27
JPWO2021260783A1 (ja) 2021-12-30

Similar Documents

Publication Publication Date Title
US8602986B2 (en) System and method for detecting signal artifacts
US9943237B2 (en) Analysis of direct and indirect heartbeat data variations
EP2962633B1 (en) Pulse wave propagation time calculation device
JP3770204B2 (ja) 脈波解析装置及び生体状態監視装置
US9044147B2 (en) Detection of noise during heart beat variation evaluation
JP2017042386A (ja) 生体情報処理システム及びプログラム
JP3729143B2 (ja) 脈波計測装置
KR20100130108A (ko) 심전도 신호 분석 방법 및 그를 이용한 심전도 장치
JP4458436B2 (ja) 心拍計および心拍波形のノイズ除去方法
US11253162B2 (en) Method and system for heart rate estimation
TWM578864U (zh) 用於穿戴式裝置之健康管理系統
WO2021260783A1 (ja) 心拍数検知方法、装置およびプログラム
US8086303B2 (en) Cardiac arrhythmias analysis of electrophysiological signals based on symbolic dynamics
JP7343055B2 (ja) 心拍数検知方法、装置およびプログラム
JPH08112270A (ja) 睡眠状態判定方法
JP7472741B2 (ja) 脈拍数測定装置及びプログラム
JP4451297B2 (ja) 生体情報検出装置
JPH05212006A (ja) 心拍間隔計測装置
US20220044816A1 (en) Information processing apparatus and non-transitory computer readable medium
JP7137550B2 (ja) 不整脈判定装置及び不整脈判定用プログラム
US12186112B2 (en) Blood pressure-related display device, blood pressure-related information display method, and computer-readable recording medium
Lee et al. Stability of watch-based blood pressure measurements analyzed by pre-post calibration differences
JP6537744B2 (ja) メンタルストレス検出装置及びメンタルストレス検出プログラム
TWI682768B (zh) 心率偵測方法
US20180132742A1 (en) Electrocardiogram (ekg) indication saturation correction

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: 20941959

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2022531269

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: 20941959

Country of ref document: EP

Kind code of ref document: A1