JP2019201997A - 瞬時心拍信頼性評価装置、方法およびプログラム - Google Patents
瞬時心拍信頼性評価装置、方法およびプログラム Download PDFInfo
- Publication number
- JP2019201997A JP2019201997A JP2018099962A JP2018099962A JP2019201997A JP 2019201997 A JP2019201997 A JP 2019201997A JP 2018099962 A JP2018099962 A JP 2018099962A JP 2018099962 A JP2018099962 A JP 2018099962A JP 2019201997 A JP2019201997 A JP 2019201997A
- Authority
- JP
- Japan
- Prior art keywords
- divided
- evaluation
- signals
- evaluated
- signal
- 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.)
- Granted
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 266
- 238000000034 method Methods 0.000 title description 61
- 238000005259 measurement Methods 0.000 claims abstract description 215
- 238000004364 calculation method Methods 0.000 claims abstract description 63
- 230000002159 abnormal effect Effects 0.000 claims abstract description 60
- 238000000605 extraction Methods 0.000 claims abstract description 21
- 230000028161 membrane depolarization Effects 0.000 claims abstract description 12
- 230000007704 transition Effects 0.000 claims description 12
- 239000000284 extract Substances 0.000 claims description 5
- 230000005856 abnormality Effects 0.000 abstract description 51
- 230000008569 process Effects 0.000 description 35
- 238000007796 conventional method Methods 0.000 description 23
- 238000012545 processing Methods 0.000 description 21
- 230000008035 nerve activity Effects 0.000 description 14
- 210000001186 vagus nerve Anatomy 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 230000002889 sympathetic effect Effects 0.000 description 8
- 210000000467 autonomic pathway Anatomy 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 210000000056 organ Anatomy 0.000 description 7
- 238000001514 detection method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 239000000470 constituent Substances 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 230000000747 cardiac effect Effects 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 3
- 230000007774 longterm Effects 0.000 description 3
- 210000005036 nerve Anatomy 0.000 description 3
- 238000004904 shortening Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000002565 electrocardiography Methods 0.000 description 2
- 238000009532 heart rate measurement Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 210000005037 parasympathetic nerve Anatomy 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000001515 vagal effect Effects 0.000 description 2
- 206010040047 Sepsis Diseases 0.000 description 1
- 230000006793 arrhythmia Effects 0.000 description 1
- 206010003119 arrhythmia Diseases 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000002301 combined effect Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 210000003792 cranial nerve Anatomy 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003340 mental effect Effects 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003183 myoelectrical effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 239000002243 precursor Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000004938 stress stimulation Effects 0.000 description 1
- 230000035900 sweating Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/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/7217—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise originating from a therapeutic or surgical apparatus, e.g. from a pacemaker
-
- 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/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/333—Recording apparatus specially adapted therefor
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Engineering & Computer Science (AREA)
- Surgery (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Signal Processing (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
【課題】適切に生体信号の異常判別を実現する。【解決手段】実施形態における瞬時心拍信頼性評価装置は、被験者の生体信号における心臓の脱分極に対応する極大値をとる波形を抽出する抽出手段と、抽出された、時系列で隣接する2つの波形の間隔を算出する第1の算出手段と、計測手段から出力される信号を予め設定された期間ごとの信号に分割する分割手段と、分割された各々の信号の電位の特徴量を算出する第2の算出手段と、分割された各々の信号の計測状態が正常であるか異常であるかを特徴量に基づいて評価する第1の評価手段と、抽出された、隣接する2つの波形のそれぞれの計測状態を、第1の評価手段による評価結果に基づいて評価し、評価された波形の計測状態の種別に応じて、第1の算出手段により算出された波形の間隔の計測状態の信頼性を評価する第2の評価手段と、を有する。【選択図】図1
Description
本発明の実施形態は、瞬時心拍信頼性評価装置、方法およびプログラムに関する。
自律神経には、交感神経と迷走神経の2つがある。交感神経および迷走神経は各臓器などに広く分布し、循環、代謝をはじめとする不随意な身体機能を制御する。多くの場合、交感神経および迷走神経が1つの臓器を拮抗的に支配すると言われている。
自律神経活動の1つである交感神経活動は、暗算負荷などのストレス刺激によって亢進することが知られている。
一方、もう1つの自律神経活動である迷走神経活動は、迷走神経は各臓器において主に副交換性の神経活動を担うことから、副交感神経活動と同等に理解されることも多い。なお、迷走神経とは、厳密には脳神経の1つである第X神経の名称であり、脳から各臓器などに至る当該神経すべてを指す。このため、例えば心臓迷走神経のように、支配対象となっている臓器の名称を付記することで、対象臓器における副交感神経活動を示す場合がある。
自律神経が支配する臓器の1つに心臓がある。心臓は交感神経および迷走神経によって拮抗的に支配されており、交感神経活動と迷走神経活動の静的なバランスを反映すると言われている。
図9は、瞬時心拍(RRI)の一例を示す図である。図9に示すような、時系列で隣接する(以下、単に隣接すると称することがある)2つのR波の間隔である瞬時心拍(RRI:R-R interval)のゆらぎは交感神経活動と迷走神経活動の両方によって変化することが知られている。なお、R波とは、心電図(ECG(Electrocardiogram))計測によって得られる心電波形のひとつであり、心臓の脱分極活動を反映している。
実環境で自律神経活動を推定する手法として、瞬時心拍変動の周波数スペクトル解析が挙げられる。この手法によれば、不等間隔である瞬時心拍を周波数スペクトルで解析した際の低周波成分(以降HRVLF)は、交感神経活動と心臓迷走神経活動とを反映する指標として解釈され、高周波成分(以降HRVHF)は、心臓迷走神経活動を反映する指標として解釈される。
心電図を計測する手段の1つとして、ホルター(Holter)心電計などのウェアラブルなデバイスがある。これらのデバイスを用いて取得される心電図は、電極の変形、ズレをはじめとする電極異常、あるいは、体動、発汗、静電気など様々な要因によって計測異常が生じる。図10は、心電図における計測異常の一例を示す図である。この計測異常は、心電図では図10に示すようなノイズ(W1)、アーチファクト(W2)(異常計測状態)という形で確認できる。なお、ノイズ、アーチファクトともに、その持続時間は計測異常の継続時間によって変化する。
図10に示すように、アーチファクトとして観測される波形は、R波と酷似した周波数特性を持つため、一般的なフィルタリングで完全に除去するのは非常に困難である。このため、心電図を解析してR波を抽出するアルゴリズムが、アーチファクトをR波と誤判断して不必要に抽出してしまう場合がある。
HRVLFおよびHRVHFは、解析対象となるすべてのデータが正常な瞬時心拍を示す場合でのみ、自律神経活動を反映する。ここでいう正常な瞬時心拍とは、計測対象と計測器両方において異常がないときの瞬時心拍を意味する。計測対象の異常とは、被験者の不整脈などを指し、計測器の異常とは、心電図で計測異常が生じている状態を指す。
計測異常の1つであるアーチファクトをR波と誤判断した波形(以降では、このような波形を計測異常R波と表記する)は、その発生機序から心臓の脱分極活動を一切反映しない。このため、解析対象となる瞬時心拍を構成するR波のうち、少なくとも1つがアーチファクトをR波と誤判断したものである場合、HRVLFおよびHRVHFのいずれも自律神経活動を反映するとは言えない。
心電図の計測異常への対処法は、フィルタリングなどによる計測異常の影響低減と、計測異常発生部分の検出との二つに大別される。後者の手法である、計測異常発生部分の検出は、前者のフィルタリングなどによる計測異常の影響低減が、特にアーチファクト発生時に適切に行えないことに着目し、心電図の電位情報の統計情報を使って計測異常発生部分の検出を図る手法である(以降、従来手法a)。
瞬時心拍の異常値を除外する手法として、瞬時心拍の時間特徴量を使用する手法(非特許文献1乃至3参照)(以降、従来手法b)と、瞬時心拍の計測状態と時間特徴量に基づく手法(非特許文献4参照)(以降、従来手法c)とがある。
従来手法bの具体的な例としては、瞬時心拍の下限値・上限値と、前後の瞬時心拍との差分値に閾値を設定し、当該閾値を逸脱したものを除外する手法、および瞬時心拍の正規分布から外れた波形を除外する手法がある。この、瞬時心拍の正規分布から外れた波形を除外する手法の中では、「瞬時心拍の平均±標準偏差」による異常値検出が最も簡便なものであり、一般的には2σあるいは3σルールが用いられることが多い。
従来手法cの具体的な例としては、瞬時心拍の計測状態をR波の電位振幅情報に基づいて判別し、異常と思われる瞬時心拍を除外した上で、従来手法bの下限値・上限値除外、3σルールを順に適用するものがある。この時、瞬時心拍の計測状態は、瞬時心拍を構成するR波の電位振幅に基づいて評価される。
佐久間大輝,神田尚子,吉見真聡,吉永努,入江英嗣,「座位状態での心拍測定を用いたリアルタイムなストレス緩和システム」,マルチメディア、分散協調とモバイルシンポジウム2013論文集,pp.1188-1195,2013
横田康成,河村洋子,松丸直樹,白井邦博,「心拍変動時系列変化を用いた敗血症の前駆症状モニタリング」,第54回自動制御連合講演会,pp.1258-1261,2011
江口佳那,角田啓介,籔内勉,吉田和広,渡邉智樹,水野理,「ウェアラブル心電計を対象としたノイズ混入領域の検出方法」,信学技報,Vol.115,No.345,pp.27-32,2015
江口佳那,青木良輔,吉田和広,山田智広,「QRS群電位特性を用いたウェアラブル心電計用のRRI計測信頼度評価」,信学技報Vol.116,No.412,pp.171-176,2017
従来手法aは、心電図の計測異常がおもに電位の違いとして見られることに着目した手法であるが、長期・短期の二種類の時間長で分割した心電図の統計情報を用いる。このため、任意の分割心電図の計測状態は、その分割心電図が属する長期分割心電図の次の長期分割心電図の統計情報を用いて判別されるため、計測状態判別には長期時間長で設定した時間分がかかる。
従来手法aにおいて、以前計測異常が生じていないと判定された複数の分割心電図との比較以外については、ある分割心電図と別の分割心電図の統計情報による比較評価が行われる。この比較は、参考情報として設定した分割心電図の定数倍であるかどうかを基準とするため、当該心電図において計測異常が生じていた場合、その影響を低減することが困難なおそれがある。
瞬時心拍は時間の経過とともに変化する値であるため、動的に閾値を設定することが難しい。すなわち、一般に従来手法bでは、瞬時心拍の異常値が病理的に正常な範囲内である場合には、当該瞬時心拍を異常として検出できない。
上記の、瞬時心拍の正規分布から外れた波形を除外する手法においては、正規分布から逸脱した値を検出するため、平均値、標準偏差、中央値、四分位数をはじめとする各種統計値に基づいて閾値を設定する。異常値の判別を行うためには、正常な瞬時心拍に対するこれらの統計値を算出する必要があるが、この段階において正常な瞬時心拍を特定することはできない。よって、異常な瞬時心拍が多くなると、その存在が正常とみなされてしまうため、当該瞬時心拍を異常として検出できなくなる。
従来手法bにおいては、計測器の異常は考慮されていない。すなわち、アーチファクトをR波と誤判断したものを少なくともひとつ以上含む瞬時心拍について、時間的な長さが病理的に正常な場合、あるいは正常な瞬時心拍によって仮定される正規分布の範囲内である場合、従来手法による異常な瞬時心拍の検出は不可能である。
従来手法cは,電位情報を用いることで、従来手法bの問題点の一部および一般的なフィルタリングの問題点を克服している。しかし、心電図に生じる計測異常は、必ずしもR波と異なる電位振幅特性を有するとは限らない。特に筋電アーチファクトなどR波と同等の電位振幅特性を誤判定したR波については、異常として検出することができない。
本発明は、上記の事情に鑑みてなされたものであり、その目的は、適切に生体信号の異常判別を実現できる瞬時心拍信頼性評価装置、方法およびプログラムを提供することにある。
上記目的を達成するために、この発明の一実施形態における瞬時心拍信頼性評価装置の第1の態様は、瞬時心拍信頼性評価装置が、被験者の生体信号を計測する計測手段から出力される信号に基づいて前記生体信号における心臓の脱分極に対応する極大値をとる波形を抽出する抽出手段と、前記抽出手段により抽出された、時系列で隣接する2つの前記波形の間隔を算出する第1の算出手段と、前記計測手段から出力される信号を予め設定された期間ごとの信号に分割する分割手段と、前記分割手段により分割された各々の信号の電位の特徴量を算出する第2の算出手段と、前記分割手段により分割された各々の信号の計測状態が正常であるか異常であるかを、前記第2の算出手段により算出された特徴量に基づいて評価する第1の評価手段と、前記抽出手段により抽出された、前記隣接する2つの前記波形のそれぞれの計測状態を、前記第1の評価手段による評価結果に基づいて評価し、前記評価された前記波形の計測状態の種別に応じて、前記第1の算出手段により算出された前記波形の間隔の計測状態の信頼性を評価する第2の評価手段と、を備えるようにしたものである。
この発明の瞬時心拍信頼性評価装置の第2の態様は、第1の態様において、前記第2の算出手段は、前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の大きさを表す第1の特徴量、前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位のばらつきを表す第2の特徴量、および前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の時系列での推移量を表す第3の特徴量のうち少なくとも1つを算出し、前記第1の評価手段は、前記分割手段により分割された信号の1つである評価対象の信号について前記第2の算出手段により算出された特徴量、および前記評価対象の信号を含んで時系列で連続する信号の集合に属する他の信号について前記第2の算出手段により算出された特徴量に基づいて、前記評価対象の信号の計測状態が正常である可能性があるか又は異常である可能性があるかを、前記第2の算出手段により算出された特徴量の種別ごとに評価し、前記特徴量の種別ごとの評価結果に基づいて、前記分割手段により分割された各々の信号の計測状態が正常であるか異常であるかを評価するようにしたものである。
この発明の瞬時心拍信頼性評価装置の第3の態様は、第1の態様において、前記第2の算出手段は、前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の大きさを表す第1の特徴量、前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位のばらつきを表す第2の特徴量、および前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の時系列での推移量を表す第3の特徴量のうち少なくとも1つを算出し、前記第1の評価手段は、前記分割手段により分割された信号の1つである評価対象の信号について前記第2の算出手段により算出された特徴量、および前記評価対象の信号に対し時系列で隣接する信号について前記第2の算出手段により算出された特徴量に基づいて、前記評価対象の信号の計測状態が正常である可能性があるか又は異常である可能性があるかを前記第2の算出手段により算出された特徴量の種別ごとに評価し、前記特徴量の種別ごとの評価結果に基づいて、前記分割手段により分割された各々の信号の計測状態が正常であるか異常であるかを評価するようにしたものである。
この発明の瞬時心拍信頼性評価装置の第4の態様は、第1の態様において、前記第2の算出手段は、前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の大きさを表す第1の特徴量、前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位のばらつきを表す第2の特徴量、および前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の時系列での推移量を表す第3の特徴量のうち少なくとも1つを算出し、前記第1の評価手段は、前記分割手段により分割された信号の1つである評価対象の信号について前記第2の算出手段により算出された特徴量、および前記分割手段により分割された信号のうち計測状態が正常であると予め評価された信号について前記第2の算出手段により算出された特徴量に基づいて、前記評価対象の信号の計測状態が正常である可能性があるか又は異常である可能性があるかを前記第2の算出手段により算出された特徴量の種別ごとに評価し、前記特徴量の種別ごとの評価結果に基づいて、前記分割手段により分割された各々の信号の計測状態が正常であるか異常であるかを評価するようにしたものである。
この発明の瞬時心拍信頼性評価装置の第5の態様は、第1の態様において、前記第2の算出手段は、前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の大きさを表す第1の特徴量、前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位のばらつきを表す第2の特徴量、および前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の時系列での推移量を表す第3の特徴量のうち少なくとも1つを算出し、前記第1の評価手段は、前記分割手段により分割された信号の1つである評価対象の信号について前記第2の算出手段により算出された特徴量、および前記評価対象の信号を含んで時系列で連続する信号の集合に属する他の信号について前記第2の算出手段により算出された特徴量に基づいて、前記評価対象の信号の計測状態が正常である可能性があるか又は異常である可能性があるかを示す第1の評価結果を前記第2の算出手段により算出された特徴量の種別ごとに得て、前記分割手段により分割された信号の1つである評価対象の信号について前記第2の算出手段により算出された特徴量、および前記評価対象の信号に対し時系列で隣接する信号について前記第2の算出手段により算出された特徴量に基づいて、前記評価対象の信号の計測状態が正常である可能性があるか又は異常である可能性があるかを示す第2の評価結果を前記第2の算出手段により算出された特徴量の種別ごとに得て、前記評価対象の信号について前記第2の算出手段により算出された特徴量、および計測状態が正常であると予め評価された信号について前記第2の算出手段により算出された特徴量に基づいて、前記評価対象の信号の計測状態が正常である可能性があるか又は異常である可能性があるかを示す第3の評価結果を前記第2の算出手段により算出された特徴量の種別ごとに得て、前記第1、第2、および第3の評価結果の少なくとも一つが、前記評価対象の信号の計測状態が異常である可能性があることを示すときに、前記評価対象の信号の計測状態が異常であると評価するようにしたものである。
この発明の瞬時心拍信頼性評価装置の第6の態様は、第1の態様において、前記第2の評価手段は、前記第1の評価手段により計測状態が異常であると評価された信号から前記抽出手段により抽出された前記波形の計測状態が異常であると評価し、前記第1の評価手段により計測状態が正常であると評価された信号から前記抽出手段により抽出された前記波形の計測状態が正常であると評価するようにしたものである。
本発明の一実施形態における、瞬時心拍信頼性評価装置が行う瞬時心拍信頼性評価方法の一つの態様は、被験者の生体信号における心臓の脱分極に対応する極大値をとる波形を抽出し、前記抽出された、時系列で隣接する2つの前記波形の間隔を算出し、前記生体信号を予め設定された期間ごとの信号に分割し、前記分割された各々の信号の電位の特徴量を算出し、前記分割された各々の信号の計測状態が正常であるか異常であるかを、前記算出された特徴量に基づいて評価し、前記抽出された、前記隣接する2つの前記波形のそれぞれの計測状態を、前記分割された各々の信号の計測状態の評価結果に基づいて評価し、前記評価された前記波形の計測状態の種別に応じて、前記算出された前記波形の間隔の計測状態の信頼性を評価するようにしたものである。
本発明の一実施形態における瞬時心拍信頼性評価処理プログラムの一つの態様は、第1乃至第6の態様のいずれか1つにおける瞬時心拍信頼性評価装置の前記各手段としてプロセッサを機能させるものである。
この発明の一実施形態における瞬時心拍信頼性評価装置の第1の態様によれば、心臓の脱分極に対応する極大値をとる波形の間隔の計測状態を、分割した生体信号の計測状態の評価結果に基づいて評価するので、生体信号の計測状態が、心臓の脱分極に対応する極大値をとる波形の間隔の計測状態の評価に及ぼす影響を低減することができる。
この発明の一実施形態における瞬時心拍信頼性評価装置の第2の態様によれば、分割された各々の信号の計測状態が正常又は異常である可能性があるかを特徴量の種別ごとに評価し、特徴量の種別ごとの評価結果に基づいて、分割された各々の信号の計測状態が正常であるか異常であるかを評価するので、生体信号の計測状態の評価の精度を向上させることができる。
この発明の一実施形態における瞬時心拍信頼性評価装置の第3乃至第5の態様によれば、評価対象の信号に時系列で前後に位置する信号、または計測状態が正常であると予め評価された信号の特徴量に基づいて、評価対象の信号の計測状態を評価するので、分割した信号の1つである評価対象の信号について算出した計測状態の評価にかかる時間を短縮できる。
この発明の一実施形態における瞬時心拍信頼性評価装置の第6の態様によれば、評価対象の生体信号の計測状態の評価結果により、心臓の脱分極に対応する極大値をとる波形の計測状態を評価するので、生体信号を分割した信号の電位の特徴量に基づいて心臓の脱分極に対応する極大値をとる波形の計測状態を適切に評価することができる。
すなわち、本発明によれば、適切に生体信号の異常判別を実現することが可能になる。
すなわち、本発明によれば、適切に生体信号の異常判別を実現することが可能になる。
以下、図面を参照しながら、この発明に係わる一実施形態を説明する。
本発明の一実施形態では、瞬時心拍を構成する2つのR波の電位情報に基づいて当該瞬時心拍の信頼性を評価する。このとき、心電図の計測状態をあらかじめ別途解析しておき、瞬時心拍を構成するR波に当該計測状態の情報を付与することで、対象とする瞬時心拍の異常判別を実現する。さらに、このとき、上記の従来手法aで課題となっていた、計測異常判別にかかる待機時間を削減しつつ、計測異常の影響を受けにくい手法とする。なお、瞬時心拍の評価においては、上記の従来手法cで述べたR波の電位振幅情報、または上記の従来手法bで述べた時間情報を考慮してもよい。
本発明の一実施形態では、瞬時心拍を構成する2つのR波の電位情報に基づいて当該瞬時心拍の信頼性を評価する。このとき、心電図の計測状態をあらかじめ別途解析しておき、瞬時心拍を構成するR波に当該計測状態の情報を付与することで、対象とする瞬時心拍の異常判別を実現する。さらに、このとき、上記の従来手法aで課題となっていた、計測異常判別にかかる待機時間を削減しつつ、計測異常の影響を受けにくい手法とする。なお、瞬時心拍の評価においては、上記の従来手法cで述べたR波の電位振幅情報、または上記の従来手法bで述べた時間情報を考慮してもよい。
(構成)
図1は、本発明の一実施形態における瞬時心拍信頼性評価装置の構成例を示す図である。図1に示す瞬時心拍信頼性評価装置10は、心電図計測部11および瞬時心拍信頼性評価部12を備える。一例として、瞬時心拍信頼性評価装置10は、心電図計測部11を被験者(ユーザ)に装着可能なウェアラブルデバイスとし、瞬時心拍信頼性評価部12をスマートフォン、タブレット型端末、パーソナルコンピュータ(PC)などのコンピュータデバイスを用いたシステムにより実現される。例えば、コンピュータデバイスは、CPU(Central Processing Unit)などのプロセッサと、プロセッサに接続されるメモリと、心電図計測部11と(例えば無線で)通信するための通信インタフェースと、を備える。メモリは、不揮発性メモリなどの記憶媒体を有する記憶装置により構成される。
図1は、本発明の一実施形態における瞬時心拍信頼性評価装置の構成例を示す図である。図1に示す瞬時心拍信頼性評価装置10は、心電図計測部11および瞬時心拍信頼性評価部12を備える。一例として、瞬時心拍信頼性評価装置10は、心電図計測部11を被験者(ユーザ)に装着可能なウェアラブルデバイスとし、瞬時心拍信頼性評価部12をスマートフォン、タブレット型端末、パーソナルコンピュータ(PC)などのコンピュータデバイスを用いたシステムにより実現される。例えば、コンピュータデバイスは、CPU(Central Processing Unit)などのプロセッサと、プロセッサに接続されるメモリと、心電図計測部11と(例えば無線で)通信するための通信インタフェースと、を備える。メモリは、不揮発性メモリなどの記憶媒体を有する記憶装置により構成される。
なお、瞬時心拍信頼性評価装置10の実現形態は、この例に限定されるものではない。例えば、瞬時心拍信頼性評価装置10は1つのデバイスとして実現されてもよい。また、心電図計測部11は瞬時心拍信頼性評価装置10の外部に設けられてもよい。言い換えると、瞬時心拍信頼性評価装置10は、心電図計測部11に相当する外部の心電計測装置から被験者の心電を計測した結果を取得してもよい。
本発明の一実施形態と従来技術との相違点は、瞬時心拍を構成する2つのR波の電位情報に基づいて当該瞬時心拍を評価する際に、あらかじめ解析しておいた心電図の計測状態を用いることが挙げられる。
心電図計測部11は、被験者の心電図を計測し、計測結果を瞬時心拍信頼性評価部12に送る。心電は、循環器系の生体信号であり、例えば、心室の収縮と同期した周期的な信号を含む。心電図計測部11は、少なくとも2極の電極によって心電の計測を行う。計測結果は、心電図におけるR波相当の心電を抽出可能なデータを含む。例えば、計測結果は心電図のデータを含む。心電図計測部11は、R波相当の心電波形を計測することができればよく、その実現形態は問わない。例えば、心電図計測部11はホルター心電計からなる。
瞬時心拍信頼性評価部12は、心電図計測部11から計測結果を受け取り、受け取った計測結果に基づいて被験者の瞬時心拍の信頼性を評価する。例えば、瞬時心拍信頼性評価部12は、R波抽出部121、R波関連情報記録部122、瞬時心拍算出部123、瞬時心拍記録部124、瞬時心拍評価部125、心電図計測状態評価部131、および瞬時心拍異常値処理部132を備える。
R波抽出部121、R波関連情報記録部122、瞬時心拍算出部123、瞬時心拍記録部124、瞬時心拍評価部125、心電図計測状態評価部131、瞬時心拍異常値処理部132の機能は、例えば、プロセッサがメモリに格納されているプログラムを読み出して実行することにより実現される。なお、これらの機能の一部または全部は、特定用途向け集積回路(ASIC)などの回路によって実現されてもよい。
心電図計測状態評価部131は、心電図の計測状態を評価する。図2は、本発明の一実施形態における瞬時心拍信頼性評価装置の心電図計測状態評価部の構成例を示す図である。図2に示すように、心電図計測状態評価部131は、信号取得部131a、信号前処理部131b、信号分割部131c、特徴量算出部131d、評価部131e、心電図記録部131fを有する。各部の具体的な実現例は後述する。
瞬時心拍信頼性評価部12の各部について説明する。R波抽出部121は、心電図計測部11で取得した心電図を解析し、R波を抽出する。また、R波抽出部121は、抽出したR波の情報を、R波関連情報記録部122に記録する。
瞬時心拍信頼性評価部12の各部について説明する。R波抽出部121は、心電図計測部11で取得した心電図を解析し、R波を抽出する。また、R波抽出部121は、抽出したR波の情報を、R波関連情報記録部122に記録する。
R波関連情報記録部122は、R波の電位振幅情報に基づいたR波の計測状態判別を併用する場合には必須でなる。R波関連情報記録部122は、不揮発性メモリなどの記憶媒体を含み、この記憶媒体に対し、R波抽出部121で抽出したR波の情報を記録する。本発明の一実施形態は、R波の計測状態として、正常計測状態とアーチファクトとの少なくとも2種類の計測状態が区別可能な情報をR波関連情報記録部122への記録対象とするが、それ以外の情報については特に指定しない。例えば、抽出されたR波が出現した時間の情報をR波関連情報記録部122への記録対象としても良い。また、具体的な記録形式については特に指定しない。
瞬時心拍算出部123は、R波抽出部121で抽出したR波に基づき、隣接するR波間の間隔である瞬時心拍を算出する。瞬時心拍算出部123は、算出した瞬時心拍の情報を瞬時心拍記録部124に記録する。
瞬時心拍記録部124は、不揮発性メモリなどの記憶媒体を含み、この記憶媒体に対し、瞬時心拍算出部123で算出した瞬時心拍の情報を記録する。瞬時心拍記録部124への具体的な記録形式については特に指定しないが、例えば、(1)瞬時心拍の行列、および(2)瞬時心拍を構成する1つ目のR波の時刻情報と瞬時心拍との2つから構成されるデータ行列、が挙げられる。
なお、瞬時心拍記録部124の機能は一実施形態において必須の機能ではない。この機能は、心電図の計測状態、R波の電位情報に加え、瞬時心拍の時間情報も考慮して瞬時心拍を評価する場合のみに必要となる。
瞬時心拍評価部125は、R波関連情報記録部122に記録される情報を基に、瞬時心拍算出部123で算出した瞬時心拍の信頼性を評価する。具体的な評価手法については、後述する動作例で述べる。
(動作例)
本発明の一実施形態の動作例について説明する。この動作例では、R波の電位振幅特性では判別が難しいアーチファクトを対象として、このアーチファクトを誤判定したものを異常値として除外する方法について説明する。
本発明の一実施形態の動作例について説明する。この動作例では、R波の電位振幅特性では判別が難しいアーチファクトを対象として、このアーチファクトを誤判定したものを異常値として除外する方法について説明する。
具体的な処理手順を以下に示す。図3は、本発明の一実施形態における瞬時心拍信頼性評価装置による処理手順の一例を示すフローチャートである。
心電図計測部11」、最低でも2極の電極によって心電図の計測を行なう(S1)。
R波抽出部121は、心電図計測部11で計測した心電図からR波を抽出する(S2)。上記の従来手法cのようにR波の電位振幅を用いる場合には、R波抽出部121は、抽出した各R波の電位情報をR波関連情報記録部122に記録する。
心電図計測部11」、最低でも2極の電極によって心電図の計測を行なう(S1)。
R波抽出部121は、心電図計測部11で計測した心電図からR波を抽出する(S2)。上記の従来手法cのようにR波の電位振幅を用いる場合には、R波抽出部121は、抽出した各R波の電位情報をR波関連情報記録部122に記録する。
瞬時心拍算出部123は、R波抽出部121で抽出したR波に基づき、隣接する2つのR波から瞬時心拍を算出する(S3)。上記の従来手法bのように、瞬時心拍の時間情報も考慮に入れて瞬時心拍の計測状態の評価を行なう場合は、瞬時心拍算出部123は、算出した瞬時心拍の情報を、瞬時心拍記録部124に記録する。
心電図計測状態評価部131は、心電図計測部11で計測した心電図を元に、この心電図の計測状態を評価する(S4)。以下にS4の詳細について説明する。図4は、本発明の一実施形態における瞬時心拍信頼性評価装置による心電図の計測状態の評価の処理手順の一例を示すフローチャートである。
心電図計測状態評価部131の信号取得部131aは、心電図計測部11で計測した心電図の信号を取得する。信号前処理部131bは、信号取得部131aで取得した心電図の低周波ノイズを除去する(S4a)。この処理は、図9にノイズとして示す低周波の基線変動を除去できる手法であれば実現手段は問わない。一例としては、低周波ノイズを除去できる手法は、固定的なカットオフ周波数を持つハイパスフィルタ、ケプストラム解析によるスペクトル包絡の除去などによる可変バンドパスフィルタなどがある。
心電図計測状態評価部131の信号分割部131cは、心電図の信号を少なくとも1周期を含む複数の期間毎に分割する。図5は、本発明の一実施形態における瞬時心拍信頼性評価装置による心電図の分割例、および心電図の計測状態の評価の第1の例を示す図である。具体的には、信号分割部131cは、評価対象となる心電図を任意の時間間隔である任意間隔Aで分割する(S4b)。具体的には、信号分割部131cは、心電図の信号を、任意間隔Aのデータを有する複数の分割ECG、例えば図5に示す分割ECG1、分割ECG2、分割ECG3、分割ECG4、分割ECG5、分割ECG6、…、分割ECGn、分割ECGn+1、分割ECGn+2、分割ECGn+3、分割ECGn+4、分割ECGn+5に分割する。ここで、複数の分割ECGは、それぞれ心電図の信号の1単位となるパケットとみなすこともできる。
心電図計測状態評価部131においては、任意間隔Aごとに区切られた心電図を、以降における最小評価単位とする。すなわち、ある任意間隔Aの心電図が計測異常と判定されたときには、心電図計測状態評価部131は、心電図の分割領域に含まれるすべての心電図を計測異常とみなす。ここでは、任意間隔Aの具体的な長さについて特に規定しないが、任意間隔Aの中で観測されるR波の数が各分割心電図でおよそ同等となることが望ましい。たとえば、設定値を、健常者の最大RRIとされる1.5[sec]などとして、分割した心電図の中に最低でも一つのR波が観測されるようにすることが考えられる。
心電図計測状態評価部131の特徴量算出部131dは、任意間隔Aで分割した各心電図について、各種の特徴量を算出する。なお、以下に示す各種の特徴量は一例であり、異なる種類の特徴量を加えても良いし、同等の特性を持つ他の種類の特徴量を代わりに用いてもよいが、計測異常時と正常時とで異なる値となる特徴量である必要がある。
このように分割した各心電図について、各種の特徴量を算出することは、心電は、類似した波形が繰り返し出現する生体信号であることに着目したものである。このため、任意間隔ごとに心電図を区切って、任意間隔ごとに各心電図の計測状態を評価する場合、対象となる心電の周波数特性と類似または同一の周波数特性、例えば類似または同一の周波数帯域、を有するノイズ(以下、「類似周波数特性ノイズ」と称す。)が混入していない理想的な状態であれば、以下に示す各種の特徴量(心電位の大きさ、ばらつき、および波形の推移量)は、どの任意間隔の領域でも近い値を得ると考えられる。
特徴量の算出の詳細について説明する。
第1に、特徴量算出部131dは、特徴量aとして、分割心電図に含まれる心電の絶対値について、平均値を算出する(S4c)。なお、この特徴量aは、「心電波形の大きさ」の大小関係を大局的に表現可能な正の値であれば、他の特徴量であっても構わない。この特徴量aによって、心電の電位の大きさに基づいた計測異常の判別を可能とする。
第1に、特徴量算出部131dは、特徴量aとして、分割心電図に含まれる心電の絶対値について、平均値を算出する(S4c)。なお、この特徴量aは、「心電波形の大きさ」の大小関係を大局的に表現可能な正の値であれば、他の特徴量であっても構わない。この特徴量aによって、心電の電位の大きさに基づいた計測異常の判別を可能とする。
第2に、特徴量算出部131dは、特徴量bとして、分割心電図に含まれる心電について、標準偏差を算出する(S4d)。なお、この本特徴量bは、「心電波形のばらつき」を大局的に表現可能な値であれば、例えば一定間隔ごとに算出可能な分散値などの他の統計量であっても構わない。この特徴量bによって、類似周波数特性ノイズによって平均から外れた値を示す心電データの数、または、この平均からの外れ量を評価することができる。これにより、心電の電位のばらつきに基づいた計測異常の判別を可能とする。
第3に、特徴量算出部131dは、特徴量cとして、分割心電図に含まれる心電(または心電の絶対値)について、差分の累積値(以降,累積差分)を算出する(S4e)。なお、この特徴量cは、時系列における「心電波形の推移量」を大局的に表現可能な値であれば、例えば一定間隔ごとに算出可能な移動平均または、微分の積分値などの他の統計量であっても構わない。特徴量算出部131dは、この特徴量cによって、例えば、値が小さく長期的に発生する異常など、局所的には検出が困難な異常を評価することができる。これにより、心電の電位の推移量に基づいた計測異常の判別を可能とする。
[0059]〜[0061]では、「心電波形の大きさ」,「心電波形のばらつき」,「心電波形の推移量」を表す各特徴量を順に算出する場合について述べたが、本実施形態では各特徴量の算出手順については特に規定しない。すなわち、各特徴量を並列に算出しても良いし、[0059]〜[0061]の記載と異なる順序で特徴量を算出しても良い。
心電図計測状態評価部131の評価部131eは、S4c〜S4eで算出した特徴量a〜cを用いて、当該分割心電図の計測状態を評価する。ここでは、S4c〜S4eに示す平均、標準偏差、および累積差分を特徴量として用いる場合について説明する。評価部131eは、ある分割心電図(被評価分割心電図)(例えば図5の分割ECGn+1)を評価するとき、この心電図に対し時系列に近傍に位置する分割心電図などを参照して、次の処理i〜iiiで評価を行う。
なお、以下に示す処理i〜iiiは特徴量の種別ごとに独立して行う。また、本実施形態では処理i〜iiiの実施順序については問わないため、並列処理などで各処理を実現しても良い。
i)処理iについて説明する。処理iでは、評価部131eは、ある種別の特徴量について、評価対象となる被評価分割心電図(例えば図5に示す分割ECGn+1)の近傍のi個の分割心電図の特徴量を参照する。従来手法よりも待機時間を短縮しつつ、分割心電図の計測異常を適切に検出するために、評価部131eは、図5に示すように、任意の被評価分割心電図(例えば図5の分割ECGn+1)の計測状態を、当該被評価分割心電図が属するクラスタの分割心電図(参照分割心電図)(例えば図5に示す分割ECGn、分割ECGn+2、分割ECGn+3、分割ECGn+4、分割ECGn+5)を参照して評価する(S4f)。つまり、処理iでは、被評価分割心電図の計測状態を、この被評価分割心電図を含んで時系列で連続した分割心電図の集合に属する他の分割心電図を参照して評価する。
このとき、S4c、S4dで算出した値について、参照分割心電図の電位の平均値μ、標準偏差σで表すとき、評価部131eは、以下の式(1)で被評価分割心電図の計測状態を評価する。式(1)の左辺は、被評価分割心電図について求めた特徴量である。式(1)の右辺は、参照分割心電図の電位の平均値μ、標準偏差σである。これによって、いずれかの参照分割心電図の計測状態が異常であった場合についても、この影響を低減することが可能となる。
Statisticstarget≦μref+k×σref …式(1)
式(1)の条件を満たす場合は、評価部131eは、被評価分割心電図について算出した特徴量の変化が大局的に見て所定の範囲内の値に収まっているとみなし、被評価分割心電図が正常なデータであると判定する。そして、評価部131eは、当該特徴量を算出するために用いた被評価分割心電図にTRUEiを付与する。
式(1)の条件を満たす場合は、評価部131eは、被評価分割心電図について算出した特徴量の変化が大局的に見て所定の範囲内の値に収まっているとみなし、被評価分割心電図が正常なデータであると判定する。そして、評価部131eは、当該特徴量を算出するために用いた被評価分割心電図にTRUEiを付与する。
一方、式(1)の条件を満たさない場合は、評価部131eは、被評価分割心電図で計測異常が発生した可能性があると評価し、特徴量を算出するために用いた被評価分割心電図にFALSEiを付与する。TRUEi、FALSEiの付与については、上記のように、特徴量の種別ごとに独立して行われる。
ii)処理iiについて説明する。図6は、本発明の一実施形態における瞬時心拍信頼性評価装置による心電図の計測状態の評価の第2の例を示す図である。従来手法よりも待機時間を短縮しつつ、分割心電図の計測異常を適切に検出するために、処理iiでは、評価部131eは、ある種別の特徴量について、図6に示すように、評価対象となる被評価分割心電図(例えば図6の分割ECGn+1)に時系列で隣接するj個の分割心電図(例えば図5の分割ECGn、分割ECGn+2)の特徴量を参照して被評価分割心電図の計測状態を評価することで(S4g)、計測異常の開始・終了近傍の適切な判別を行う。このとき、処理iと同じく、評価部131eは、式(1)で被評価分割心電図の計測状態を評価することによって、いずれかの参照分割心電図の計測状態が異常であった場合についても、この影響の低減化を図る。
式(1)の条件を満たす場合は、評価部131eは、被評価分割心電図について算出した特徴量の変化が大局的に見て所定の範囲内の値に収まっているとみなし、被評価分割心電図が正常なデータであると判定する。そして、評価部131eは、特徴量を算出するために用いた被評価分割心電図にTRUEiiを付与する。
一方、式(1)の条件を満たさない場合は、評価部131eは、被評価分割心電図で計測異常が発生した可能性があると評価し、特徴量を算出するために用いた被評価分割心電図にFALSEiiを付与する。TRUEii、FALSEiiの付与については、上記のように、特徴量の種別ごとに独立して行われる。
iii)処理iiiについて説明する。図7は、本発明の一実施形態における瞬時心拍信頼性評価装置による心電図の計測状態の評価の第3の例を示す図である。従来手法よりも待機時間を短縮しつつ、分割心電図の計測異常を適切に検出するために、処理iiiでは、評価部131eは、ある種別の特徴量について、後述のS4jの処理であらかじめ「正常」とみなされたk個の過去の分割心電図の特徴量を心電図記録部131fから参照し、処理i、処理iiと同じく、評価部131eは、式(1)で被評価分割心電図の計測状態を評価する(S4h)。
式(1)の条件を満たす場合は、評価部131eは、被評価分割心電図について算出した特徴量の変化が大局的に見て所定の範囲内の値に収まっているとみなし、被評価分割心電図が正常なデータであると判定する。そして、評価部131eは、特徴量を算出するために用いた被評価分割心電図にTRUEiiiを付与する。
一方、式(1)の条件を満たさない場合は、評価部131eは、被評価分割心電図で計測異常が発生した可能性があると評価し、特徴量を算出するために用いた被評価分割心電図にFALSEiiiを付与する。TRUEiii、FALSEiiiの付与については、上記のように、特徴量の種別ごとに独立して行われる。
次に、S4f〜S4h(処理i,ii,iii)にて、各種の特徴量について全ての分割心電図に評価結果が付与された後の総合評価について説明する。図8は、各評価結果と総合評価結果との関係を表形式で示す図である。
総合評価は、図8に示すように、同じデータから算出された特徴量についての処理i,ii,iiiでの各評価結果がAnd条件を満たすか否かにより決定する場合について説明する。この総合評価は、特徴量の種別ごとに独立してなされる。
Data#(j)は、j番目の被評価分割心電図を示している。例えば、「1:1−300」は1番目の被評価分割心電図に関するデータを示している。そして、1−300は、例えば、1番目の被評価分割心電図における心電図のサンプリング番号を示している。
Data#(j)は、j番目の被評価分割心電図を示している。例えば、「1:1−300」は1番目の被評価分割心電図に関するデータを示している。そして、1−300は、例えば、1番目の被評価分割心電図における心電図のサンプリング番号を示している。
そして、総合評価結果rを示している総合評価結果r項目、処理iの評価結果を示している処理i項目、処理iiの評価結果を示している処理ii項目、および処理iiiの評価結果を示している処理iii項目、の各々において、各Data#(j)に対応するように、TRUE評価結果またはFALSE評価結果を示している。
例えば、「2:301−600」は2番目の被評価分割心電図に関するデータを示しており、処理i項目、処理ii項目、および処理iii項目の何れにおいてもTRUE評価結果であることを示している。したがって、And条件を満たすため、総合評価結果r項目には、TRUE評価結果であることが示されている。
一方、例えば、「5:1201−1500」は5番目の被評価分割心電図に関するデータを示しており、処理i項目は、TRUE評価結果であることを示しており、処理ii項目および処理iii項目は、FALSE評価結果であることを示している。したがって、And条件を満たさないため、総合評価結果r項目には、FALSE評価結果であることが示されている。
つまり、ある被評価分割心電図における、ある種別の特徴量について、処理i,ii,iii項目のうち1つ以上の処理についての評価結果がFALSEである場合、心電図計測状態評価部131の評価部131eは、当該種別の特徴量にかかる被評価分割心電図に計測異常が発生したと総合的に評価し、この心電図にFALSE特徴量名を付与する。
特徴量a,b,cのすべてについて、被評価分割心電図に対してS4f〜S4hで総合的な評価結果が付与された後で、心電図計測状態評価部131の評価部131eは、当該被評価分割心電図の最終評価を行う(S4j)。例えば、ある被評価分割心電図の各種の特徴量a,b,cの総合的な評価結果のうち1つ以上での総合的な評価結果がFALSEである場合、心電図計測状態評価部131の評価部131eは、当該被評価分割心電図を計測異常と最終評価する。一方、ある被評価分割心電図の各特徴量a,b,cの総合的な評価結果が全てTRUEである場合、心電図計測状態評価部131の評価部131eは、当該被評価分割心電図を正常と最終評価する。
また、ある被評価分割心電図における、各種の特徴量について、処理i,ii,iiiの全てについての評価結果が計測異常ではない場合、つまり、S4f〜S4h(処理i,ii,iii)のすべての評価結果がTRUEである場合、心電図計測状態評価部131の評価部131eは、当該被評価分割心電図を正常と最終評価し、この心電図を時間的に最も新しいデータである正解値として、心電図記録部131fに保存されている特徴量のうち、時間的に最も古い心電図に関するデータと置き換える。このデータは上記の処理iiiにおける参照に用いられる。
ある被評価分割心電図について、S4f〜S4hあるいはS4iで評価結果がFALSEとなった場合、S4f〜S4hに代わって、現在の被評価分割心電図よりも後に評価の対象となる分割心電図に対する評価結果がTRUEとなるまで、心電図計測状態評価部131の評価部131eは、以下の式(2)に基づいて被評価分割心電図の評価を行う。
Statisticstarget≦μref+w×k×σref …式(2)
電極ずれが生じ、その電極が別の場所で再度固定されたような場合、すなわち、始めの電位と異なる電位で安定的に生体信号が取得出来るようになった場合、電極がずれる前と同等の電位で生体信号を正常に取得できない場合がある。そこで、上記の式(1)の右辺の第2項の標準偏差の重みを意味する係数kをw倍した式(2)を用いることで、正常な状態へ復帰することができる。
電極ずれが生じ、その電極が別の場所で再度固定されたような場合、すなわち、始めの電位と異なる電位で安定的に生体信号が取得出来るようになった場合、電極がずれる前と同等の電位で生体信号を正常に取得できない場合がある。そこで、上記の式(1)の右辺の第2項の標準偏差の重みを意味する係数kをw倍した式(2)を用いることで、正常な状態へ復帰することができる。
式(2)を用いることによって、正常値であるとみなす範囲が広くなり、正常な状態であるか否かの評価条件を緩和することができる。そして、異常な状態から復帰したことの検出漏れを抑制することもできる。
瞬時心拍評価部125は、S4での心電図計測状態評価部131による評価結果に基づいて、瞬時心拍を評価する(S5)。なお、この処理に加えて、上記の従来手法cを組み合わせて実施しても良いが、その際はS5の処理よりも先行して上記の従来手法cで述べた、R波の電位振幅特性に基づく瞬時心拍の計測状態評価を実施しておく。
本実施形態では、上記の心電図計測状態評価部131による評価結果のみが示す値を用いる場合について説明する。R波の計測状態の例として、正常計測状態とアーチファクトの二種を考慮するとき、瞬時心拍評価部125は、心電図計測状態評価部131で計測異常と判別された分割心電図において検出したすべてのR波を計測異常とみなし、当該R波をアーチファクトとみなす。
また、瞬時心拍評価部125は、心電図計測状態評価部131で正常と判別された分割心電図において検出したすべてのR波については、正常とみなす。
また、瞬時心拍評価部125は、心電図計測状態評価部131で正常と判別された分割心電図において検出したすべてのR波については、正常とみなす。
瞬時心拍評価部125は、各R波に付与ついて判別された計測状態に基づき、隣接する2つのR波から構成される瞬時心拍の計測状態を評価する。正常計測状態、アーチファクトの2種を考慮する場合、瞬時心拍を構成するR波の計測状態の判別結果の組み合わせは、以下の表1に示す通し番号#1〜#4で示されるパタンのいずれかとなる。なお、通し番号の形式は上記に限られない。
表1における判別結果の「R」は正常計測状態を示し、「A」はアーチファクトを示す。つまり、表1における通し番号#1に対応する判別結果の「R、R」は、隣接する1つ目および2つ目のR波の計測状態の判別結果がともに正常計測状態であることを示す。
表1における通し番号#2に対応する判別結果の「R、A」は、隣接する1つ目のR波の計測状態の判別結果が正常状態で、2つ目のR波の計測状態の判別結果がアーチファクトであることを示す。表1における通し番号#3に対応する判別結果の「A、R」は、隣接する1つ目のR波の計測状態の判別結果がアーチファクトで、2つ目のR波の計測状態の判別結果が正常状態であることを示す。表1における通し番号#4に対応する判別結果の「A、A」は、隣接する1つ目および2つ目のR波の計測状態の判別結果がともにアーチファクトであることを示す。
表1における「状態の詳細」は、表1における同じ行の判別結果に基づく、隣接するR波の計測状態の詳細を示す。
「状態の詳細」の表現では、瞬時心拍を構成する2つのR波の計測状態の判別結果の組み合わせのみを区別し、時系列の前後を区別しない。つまり、表1に示した例では、計測状態の判別結果の組み合わせは通し番号#1〜#4に対応する4通りであるが、#2、#3に対応する「状態の詳細」は同じ「片方は正常計測状態、他方はアーチファクト」となる。このため、表1における「状態の詳細」は3通りである。
「状態の詳細」の表現では、瞬時心拍を構成する2つのR波の計測状態の判別結果の組み合わせのみを区別し、時系列の前後を区別しない。つまり、表1に示した例では、計測状態の判別結果の組み合わせは通し番号#1〜#4に対応する4通りであるが、#2、#3に対応する「状態の詳細」は同じ「片方は正常計測状態、他方はアーチファクト」となる。このため、表1における「状態の詳細」は3通りである。
本実施形態では、表1の「状態の詳細」ごとの評価を行う場合について説明するが、それ以外の評価基準を設けても良い。例えば、同じ「状態の詳細」であっても、2つのR波の前後の情報が判別可能な評価基準を設けても良い。
瞬時心拍評価部125は、状態の詳細がユーザに容易に区別できるよう、各状態に別個の評価値を割り当てる。評価値の一例を表1の「評価値」として示す。なお、この評価値はあくまでも一例であり、本実施形態では評価値の決め方を特に限定しない。次に、表1における「評価値」について説明する。この「評価値」は、表1における同じ行の「状態の詳細」で示される状態である、瞬時心拍を構成する2つのR波のそれぞれの計測状態の信頼性を、0から1の間の数値で表現し、「状態の詳細」で示される各状態に対して、任意の評価値を割り当てたものである。なお、評価値の範囲、および各状態に対する評価値の刻み方は特に限られず、例えば1から10の間で各状態に対し1刻みで異なる評価値を割り当ててもよいし、また、各状態の間で評価値の刻み幅が異なっていてもよい。また、評価値に代えて、例えば、横棒グラフの長さ(例えば横棒グラフの長さが長いほど信頼性が高い)などを用いてもよい。
評価値の具体例を説明する。表1に示すように、通し番号#1における「状態の詳細」が「2つとも正常計測状態」のとき、この通し番号#1における「評価値」は、最高値の「1」となる。
また、「状態の詳細」の表現と同じく、「評価値」の表現では、瞬時心拍を構成する2つのR波の計測状態の判別結果の組み合わせのみを区別し、時系列の前後を区別しない。つまり、通し番号#2、#3における「状態の詳細」は、共通した「片方は正常計測状態、他方はアーチファクト」であり、これら通し番号#2、#3における「評価値」は、#1における「評価値」に対して0.6減じた、共通した「0.4」である。
また、「状態の詳細」の表現と同じく、「評価値」の表現では、瞬時心拍を構成する2つのR波の計測状態の判別結果の組み合わせのみを区別し、時系列の前後を区別しない。つまり、通し番号#2、#3における「状態の詳細」は、共通した「片方は正常計測状態、他方はアーチファクト」であり、これら通し番号#2、#3における「評価値」は、#1における「評価値」に対して0.6減じた、共通した「0.4」である。
通し番号#4における「状態の詳細」は、「2つともアーチファクト」であり、この通し番号#4における「評価値」は、#2、#3における「評価値」に対して0.4減じた、最低値の「0」である。
瞬時心拍異常値処理部132は、瞬時心拍評価部125による評価結果に基づいて、異常値判別に設定された評価値よりも低い値の瞬時心拍を異常値とみなして、後段処理に引き渡す瞬時心拍から除外する(S6)。例えば、心拍特徴量の算出など、アーチファクトを後段処理に引き渡す瞬時心拍に含めてはならない場合には評価値「1」が必要となるため、瞬時心拍異常値処理部132は、この値以下の評価値を有する瞬時心拍を異常値とみなして、後段処理に引き渡す瞬時心拍から除外する。
以上が本実施形態の概要である。本実施形態では、上記処理i〜iiiでの評価に用いるi,j,kの設定値については詳細に規定しないが、例えば、処理iiにj=2(すなわち、被評価分割心電図の前後の分割心電図)を用いることが考えられる。このとき、処理i,iiiでの評価に用いる分割心電図の個数をj<i,kとすることで、処理iiのみでは完全な検出が困難な長期的な計測異常を検出することが可能となる。
なお、上記のS4iでの総合評価、S4jでの最終評価については、評価部131eは、ある被評価分割心電図について、上記のS4f〜S4hのような、1つ以上のFALSE評価の有無をTRUEまたはFALSEを付与する基準とする、つまりTRUEのAND条件で評価するのではなく、ある被評価分割心電図について、上記のS4f〜S4hのうち、FALSEが付与された評価の個数を基準としてTRUEまたはFALSEを付与してもよい。
また、本実施形態では、周波数特性および振幅特性からはR波との判別が困難であるアーチファクトのみを異常値除外の対象とする場合について説明したが、低周波の基線変動を心電図の計測異常に含め、これを異常値除外の対象に含めてもよい。このような場合には、心電図計測状態評価部131は、S4aで行なっていたノイズ除去を割愛し、S3の直後にS4bを行なってもよい。
また、上記に示した計測異常判別については、心電図に限らず、脈波などの心電図と同様の特徴を有する循環器系の信号に対して、計測異常判別を行なってもよい。脈波について計測異常判別を行なうときは、心電で用いたR波に代えて、脈波のP波に対応して極大値をとる波形(心臓の脱分極に対応する極大値をとる波形)を用いればよい。
次に、本発明の一実施形態によって生じる効果について説明する。
本発明の一実施形態では、心電図の計測状態をあらかじめ評価しておき、この評価結果を用いて瞬時心拍の異常値を除外することで、周波数特性および電位振幅特性のみでは判別が困難なアーチファクトが心拍変動解析に及ぼす影響を低減する。
また、本発明の一実施形態では、心電図の計測状態評価において、評価方法に応じて複数種類の分割心電図を参照した上で、この評価結果の統計を用いた最終評価を行なうことによって、心電図の計測状態判別にかかる待機時間と、参照分割心電図の計測状態が評価に及ぼす影響との双方を低減する。
本発明の一実施形態では、心電図の計測状態をあらかじめ評価しておき、この評価結果を用いて瞬時心拍の異常値を除外することで、周波数特性および電位振幅特性のみでは判別が困難なアーチファクトが心拍変動解析に及ぼす影響を低減する。
また、本発明の一実施形態では、心電図の計測状態評価において、評価方法に応じて複数種類の分割心電図を参照した上で、この評価結果の統計を用いた最終評価を行なうことによって、心電図の計測状態判別にかかる待機時間と、参照分割心電図の計測状態が評価に及ぼす影響との双方を低減する。
なお、本願発明は、上記実施形態に限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で種々に変形することが可能である。また、各実施形態は適宜組み合わせて実施してもよく、その場合組み合わせた効果が得られる。更に、上記実施形態には種々の発明が含まれており、開示される複数の構成要件から選択された組み合わせにより種々の発明が抽出され得る。例えば、実施形態に示される全構成要件からいくつかの構成要件が削除されても、課題が解決でき、効果が得られる場合には、この構成要件が削除された構成が発明として抽出され得る。
また、各実施形態に記載した手法は、計算機(コンピュータ)に実行させることができるプログラム(ソフトウエア手段)として、例えば磁気ディスク(フロッピー(登録商標)ディスク、ハードディスク等)、光ディスク(CD−ROM、DVD、MO等)、半導体メモリ(ROM、RAM、フラッシュメモリ等)等の記録媒体に格納し、また通信媒体により伝送して頒布することもできる。なお、媒体側に格納されるプログラムには、計算機に実行させるソフトウエア手段(実行プログラムのみならずテーブル、データ構造も含む)を計算機内に構成させる設定プログラムをも含む。本装置を実現する計算機は、記録媒体に記録されたプログラムを読み込み、また場合により設定プログラムによりソフトウエア手段を構築し、このソフトウエア手段によって動作が制御されることにより上述した処理を実行する。なお、本明細書でいう記録媒体は、頒布用に限らず、計算機内部あるいはネットワークを介して接続される機器に設けられた磁気ディスク、半導体メモリ等の記憶媒体を含むものである。
10…瞬時心拍信頼性評価装置、11…心電図計測部、12…瞬時心拍信頼性評価部、121…R波抽出部、122…R波関連情報記録部、123…瞬時心拍算出部、124…瞬時心拍記録部、125…瞬時心拍評価部、131…心電図計測状態評価部、132…瞬時心拍異常値処理部。
Claims (8)
- 被験者の生体信号を計測する計測手段から出力される信号に基づいて前記生体信号における心臓の脱分極に対応する極大値をとる波形を抽出する抽出手段と、
前記抽出手段により抽出された、時系列で隣接する2つの前記波形の間隔を算出する第1の算出手段と、
前記計測手段から出力される信号を予め設定された期間ごとの信号に分割する分割手段と、
前記分割手段により分割された各々の信号の電位の特徴量を算出する第2の算出手段と、
前記分割手段により分割された各々の信号の計測状態が正常であるか異常であるかを、前記第2の算出手段により算出された特徴量に基づいて評価する第1の評価手段と、
前記抽出手段により抽出された、前記隣接する2つの前記波形のそれぞれの計測状態を、前記第1の評価手段による評価結果に基づいて評価し、前記評価された前記波形の計測状態の種別に応じて、前記第1の算出手段により算出された前記波形の間隔の計測状態の信頼性を評価する第2の評価手段と、
を備えた瞬時心拍信頼性評価装置。 - 前記第2の算出手段は、
前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の大きさを表す第1の特徴量、前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位のばらつきを表す第2の特徴量、および前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の時系列での推移量を表す第3の特徴量のうち少なくとも1つを算出し、
前記第1の評価手段は、
前記分割手段により分割された信号の1つである評価対象の信号について前記第2の算出手段により算出された特徴量、および前記評価対象の信号を含んで時系列で連続する信号の集合に属する他の信号について前記第2の算出手段により算出された特徴量に基づいて、前記評価対象の信号の計測状態が正常である可能性があるか又は異常である可能性があるかを、前記第2の算出手段により算出された特徴量の種別ごとに評価し、
前記特徴量の種別ごとの評価結果に基づいて、前記分割手段により分割された各々の信号の計測状態が正常であるか異常であるかを評価する、
請求項1に記載の瞬時心拍信頼性評価装置。 - 前記第2の算出手段は、
前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の大きさを表す第1の特徴量、前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位のばらつきを表す第2の特徴量、および前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の時系列での推移量を表す第3の特徴量のうち少なくとも1つを算出し、
前記第1の評価手段は、
前記分割手段により分割された信号の1つである評価対象の信号について前記第2の算出手段により算出された特徴量、および前記評価対象の信号に対し時系列で隣接する信号について前記第2の算出手段により算出された特徴量に基づいて、前記評価対象の信号の計測状態が正常である可能性があるか又は異常である可能性があるかを前記第2の算出手段により算出された特徴量の種別ごとに評価し、
前記特徴量の種別ごとの評価結果に基づいて、前記分割手段により分割された各々の信号の計測状態が正常であるか異常であるかを評価する、
請求項1に記載の瞬時心拍信頼性評価装置。 - 前記第2の算出手段は、
前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の大きさを表す第1の特徴量、前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位のばらつきを表す第2の特徴量、および前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の時系列での推移量を表す第3の特徴量のうち少なくとも1つを算出し、
前記第1の評価手段は、
前記分割手段により分割された信号の1つである評価対象の信号について前記第2の算出手段により算出された特徴量、および前記分割手段により分割された信号のうち計測状態が正常であると予め評価された信号について前記第2の算出手段により算出された特徴量に基づいて、前記評価対象の信号の計測状態が正常である可能性があるか又は異常である可能性があるかを前記第2の算出手段により算出された特徴量の種別ごとに評価し、
前記特徴量の種別ごとの評価結果に基づいて、前記分割手段により分割された各々の信号の計測状態が正常であるか異常であるかを評価する、
請求項1に記載の瞬時心拍信頼性評価装置。 - 前記第2の算出手段は、
前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の大きさを表す第1の特徴量、前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位のばらつきを表す第2の特徴量、および前記分割手段により分割された信号の各々に対して、前記分割された各々の信号の電位の時系列での推移量を表す第3の特徴量のうち少なくとも1つを算出し、
前記第1の評価手段は、
前記分割手段により分割された信号の1つである評価対象の信号について前記第2の算出手段により算出された特徴量、および前記評価対象の信号を含んで時系列で連続する信号の集合に属する他の信号について前記第2の算出手段により算出された特徴量に基づいて、前記評価対象の信号の計測状態が正常である可能性があるか又は異常である可能性があるかを示す第1の評価結果を前記第2の算出手段により算出された特徴量の種別ごとに得て、
前記分割手段により分割された信号の1つである評価対象の信号について前記第2の算出手段により算出された特徴量、および前記評価対象の信号に対し時系列で隣接する信号について前記第2の算出手段により算出された特徴量に基づいて、前記評価対象の信号の計測状態が正常である可能性があるか又は異常である可能性があるかを示す第2の評価結果を前記第2の算出手段により算出された特徴量の種別ごとに得て、
前記評価対象の信号について前記第2の算出手段により算出された特徴量、および計測状態が正常であると予め評価された信号について前記第2の算出手段により算出された特徴量に基づいて、前記評価対象の信号の計測状態が正常である可能性があるか又は異常である可能性があるかを示す第3の評価結果を前記第2の算出手段により算出された特徴量の種別ごとに得て、
前記第1、第2、および第3の評価結果の少なくとも一つが、前記評価対象の信号の計測状態が異常である可能性があることを示すときに、前記評価対象の信号の計測状態が異常であると評価する、
請求項1に記載の瞬時心拍信頼性評価装置。 - 前記第2の評価手段は、
前記第1の評価手段により計測状態が異常であると評価された信号から前記抽出手段により抽出された前記波形の計測状態が異常であると評価し、
前記第1の評価手段により計測状態が正常であると評価された信号から前記抽出手段により抽出された前記波形の計測状態が正常であると評価する、
請求項1に記載の瞬時心拍信頼性評価装置。 - 瞬時心拍信頼性評価装置が行う瞬時心拍信頼性評価方法であって、
被験者の生体信号における心臓の脱分極に対応する極大値をとる波形を抽出し、
前記抽出された、時系列で隣接する2つの前記波形の間隔を算出し、
前記生体信号を予め設定された期間ごとの信号に分割し、
前記分割された各々の信号の電位の特徴量を算出し、
前記分割された各々の信号の計測状態が正常であるか異常であるかを、前記算出された特徴量に基づいて評価し、
前記抽出された、前記隣接する2つの前記波形のそれぞれの計測状態を、前記分割された各々の信号の計測状態の評価結果に基づいて評価し、
前記評価された前記波形の計測状態の種別に応じて、前記算出された前記波形の間隔の計測状態の信頼性を評価する、
瞬時心拍信頼性評価方法。 - 請求項1乃至6のいずれか1項に記載の瞬時心拍信頼性評価装置の前記各手段としてプロセッサを機能させる瞬時心拍信頼性評価処理プログラム。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018099962A JP6839680B2 (ja) | 2018-05-24 | 2018-05-24 | 瞬時心拍信頼性評価装置、方法およびプログラム |
US17/056,269 US11510605B2 (en) | 2018-05-24 | 2019-05-22 | Instantaneous heartbeat reliability evaluation device, method, and program |
PCT/JP2019/020321 WO2019225662A1 (ja) | 2018-05-24 | 2019-05-22 | 瞬時心拍信頼性評価装置、方法およびプログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018099962A JP6839680B2 (ja) | 2018-05-24 | 2018-05-24 | 瞬時心拍信頼性評価装置、方法およびプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2019201997A true JP2019201997A (ja) | 2019-11-28 |
JP6839680B2 JP6839680B2 (ja) | 2021-03-10 |
Family
ID=68616018
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2018099962A Active JP6839680B2 (ja) | 2018-05-24 | 2018-05-24 | 瞬時心拍信頼性評価装置、方法およびプログラム |
Country Status (3)
Country | Link |
---|---|
US (1) | US11510605B2 (ja) |
JP (1) | JP6839680B2 (ja) |
WO (1) | WO2019225662A1 (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021245902A1 (ja) * | 2020-06-05 | 2021-12-09 | 三菱電機株式会社 | 生体信号処理装置および見守りシステム、見守り方法 |
JP7162232B1 (ja) | 2022-08-05 | 2022-10-28 | 公立大学法人会津大学 | 心拍分類装置及び入浴者見守りシステム |
WO2023089832A1 (ja) * | 2021-11-22 | 2023-05-25 | 日本電信電話株式会社 | 評価装置、評価方法、及び評価プログラム |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017080026A (ja) * | 2015-10-27 | 2017-05-18 | 日本電信電話株式会社 | ノイズ判定装置、方法、およびプログラム |
JP2018011819A (ja) * | 2016-07-22 | 2018-01-25 | 日本電信電話株式会社 | 生体信号処理方法および装置 |
JP2018020018A (ja) * | 2016-08-05 | 2018-02-08 | 日本電信電話株式会社 | 生体信号処理装置、異常判別方法およびプログラム |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4589420A (en) * | 1984-07-13 | 1986-05-20 | Spacelabs Inc. | Method and apparatus for ECG rhythm analysis |
-
2018
- 2018-05-24 JP JP2018099962A patent/JP6839680B2/ja active Active
-
2019
- 2019-05-22 WO PCT/JP2019/020321 patent/WO2019225662A1/ja active Application Filing
- 2019-05-22 US US17/056,269 patent/US11510605B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017080026A (ja) * | 2015-10-27 | 2017-05-18 | 日本電信電話株式会社 | ノイズ判定装置、方法、およびプログラム |
JP2018011819A (ja) * | 2016-07-22 | 2018-01-25 | 日本電信電話株式会社 | 生体信号処理方法および装置 |
JP2018020018A (ja) * | 2016-08-05 | 2018-02-08 | 日本電信電話株式会社 | 生体信号処理装置、異常判別方法およびプログラム |
Non-Patent Citations (1)
Title |
---|
江口 佳那 他3名: "QRS群電位特性を用いたウェアラブル心電計用のRRI計測信頼度評価", 電子情報通信学会技術研究報告, vol. 116, no. 412, JPN6019030670, 12 January 2017 (2017-01-12), JP, pages 171 - 176, ISSN: 0004440135 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021245902A1 (ja) * | 2020-06-05 | 2021-12-09 | 三菱電機株式会社 | 生体信号処理装置および見守りシステム、見守り方法 |
JP7433428B2 (ja) | 2020-06-05 | 2024-02-19 | 三菱電機株式会社 | 生体信号処理装置および見守りシステム、見守り方法 |
WO2023089832A1 (ja) * | 2021-11-22 | 2023-05-25 | 日本電信電話株式会社 | 評価装置、評価方法、及び評価プログラム |
JP7162232B1 (ja) | 2022-08-05 | 2022-10-28 | 公立大学法人会津大学 | 心拍分類装置及び入浴者見守りシステム |
JP2024022122A (ja) * | 2022-08-05 | 2024-02-16 | 公立大学法人会津大学 | 心拍分類装置及び入浴者見守りシステム |
Also Published As
Publication number | Publication date |
---|---|
WO2019225662A1 (ja) | 2019-11-28 |
JP6839680B2 (ja) | 2021-03-10 |
US20210244338A1 (en) | 2021-08-12 |
US11510605B2 (en) | 2022-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2010510850A (ja) | 改良された自動セグメント分割を含む生物医学信号解析方法 | |
WO2012051320A2 (en) | Motion and noise artifact detection for ecg data | |
JP6689183B2 (ja) | 状態推定装置、方法およびプログラム | |
WO2019225662A1 (ja) | 瞬時心拍信頼性評価装置、方法およびプログラム | |
Ng et al. | Iterative method to detect atrial activations and measure cycle length from electrograms during atrial fibrillation | |
CA2990877A1 (en) | Data processing apparatus for assessing a condition of a myocardium | |
US20220125328A1 (en) | Time features calculation apparatus, calculation method and its program | |
Lee et al. | A real-time abnormal beat detection method using a template cluster for the ECG diagnosis of IoT devices | |
US8073218B2 (en) | Method for detecting bio signal features in the presence of noise | |
JP6527830B2 (ja) | 生体信号処理装置、方法、およびプログラム | |
JP6692283B2 (ja) | 瞬時心拍信頼性評価装置、方法およびプログラム | |
WO2019185768A1 (en) | Method to analyze cardiac rhythms using beat-to-beat display plots | |
WO2019100564A1 (zh) | 检测报告数据生成方法 | |
JP2020517337A (ja) | アーチファクト耐性脈拍変動測定 | |
Magrupov et al. | ECG signal processing algorithms to determine heart rate | |
JP6764830B2 (ja) | 不整脈判定装置、不整脈判定方法および不整脈判定処理プログラム | |
JP6857582B2 (ja) | 瞬時心拍の時系列データの補完装置、補完方法及びそのプログラム | |
JP6468986B2 (ja) | ノイズ判定装置、方法、およびプログラム | |
Tarvainen et al. | Analysis and preprocessing of hrv—kubios HRV software | |
JP6535575B2 (ja) | 生理指標検出装置、生理指標検出方法、生理指標検出処理プログラム | |
JP6807279B2 (ja) | 期外収縮判別装置、期外収縮判別方法及びプログラム | |
Pulavskyi et al. | Evaluation of the effectiveness of post-filtration smoothing using lossless compression for heart rate variability obtained from a very noisy ECG | |
Anisimov et al. | Extrasystole Episodes Detection Algorithm in Real Time Mode | |
JP2019069044A (ja) | 瞬時心拍の時系列データの補完装置、補完方法及びそのプログラム | |
JP6846288B2 (ja) | 瞬時心拍の時系列データの補完装置、補完方法及びそのプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200820 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20210209 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210215 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6839680 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |