(57)【特許請求の範囲】
【請求項１】
バリストカルジオグラム信号を分析して心拍数を決定する方法であって、
前記バリストカルジオグラム信号中の第１の心拍の最初の時間推定値を決定するステップと、
前記最初の時間推定値を用いて、前記バリストカルジオグラム信号中の後続の心拍の推定値を繰り返し計算するステップとを有し、
前記計算するステップの各繰り返しは、複数のスコアリング関数の重み付けした和を含む目標関数の評価を含み、
前記バリストカルジオグラム信号中の後続の心拍の推定値を計算するステップの各繰り返しは、その前の計算するステップの繰り返しで見つけた時間推定値後の目標区間に制限される、方法。
【請求項２】
目標関数の評価は、心拍に対応する前記目標関数のピークの特定を含む、
請求項１に記載の方法。
【請求項３】
前記スコアリング関数のうちの１つは、前記バリストカルジオグラム信号中の高周波成分の発生を反映する、
請求項１または２に記載の方法。
【請求項４】
前記バリストカルジオグラム信号中の高周波成分の発生を反映する前記スコアリング関数は、
前記バリストカルジオグラム信号をフィルタして、高周波成分を取りだし、前記フィルタしたバリストカルジオグラム信号を二乗して、
前記フィルタして二乗したバリストカルジオグラム信号をフィルタして、得られた信号を与えること、により決定する、
請求項３に記載の方法。
【請求項５】
前記バリストカルジオグラム信号中の第１の心拍の最初の時間推定値を決定するステップは、
前記バリストカルジオグラム信号中の高周波成分の発生を反映する前記スコアリング関数の前記最大点を選択するステップを有する、前記最大点は、前記バリストカルジオグラム信号の始めの区間から選択可能
請求項４に記載の方法。
【請求項６】
前記バリストカルジオグラム信号中の第１の心拍の最初の時間推定値を決定するステップは、
最初の推定値と、前記バリストカルジオグラム信号中の前記第２の心拍の可能性のある時間推定値とについて、前記目的関数を評価するステップと、
前記最初の時間推定値と可能性のある時間推定値を、前記最初の２つの心拍について前記目的関数を最大化する値として選択するステップとを有する、
請求項１ないし４いずれか一項に記載の方法。
【請求項７】
前記バリストカルジオグラム信号中の第１の心拍の最初の時間推定値を決定するステップは、
前記バリストカルジオグラム信号の始めの前記時間値を前記最初の時間推定値として選択するステップを有する、
請求項１ないし４いずれか一項に記載の方法。
【請求項８】
計算するステップを複数回繰り返した後、前記方法はさらに前記最初の時間推定値を、
改良する各目的区間を、前の計算するステップの繰り返しで見つけた時間推定値の周りの小さな近傍の目的区間に限定して、さらに計算するステップを実行することにより改良するステップを有する、
請求項１ないし７いずれか一項に記載の方法。
【請求項９】
前記スコアリング関数のうちの１つは、前記心拍数の長期的予測を反映する、
請求項１ないし８いずれか一項に記載の方法。
【請求項１０】
前記心拍数の長期的推定を反映する前記スコアリング関数は、
前の計算するステップの繰り返しで見つけた時間推定に応じた複数の時点で、前記バリストカルジオグラム信号のスペクトログラムを評価することにより決定する、
請求項９に記載の方法。
【請求項１１】
前記心拍数の長期的推定を反映する前記スコアリング関数は、前記バリストカルジオグラム信号の自己相関関数を評価することにより決定する、
請求項９に記載の方法。
【請求項１２】
前記スコアリング関数の１つは、心拍間隔の確率的振る舞いに関する情報を反映する、請求項１ないし１１いずれか一項に記載の方法。
【請求項１３】
心拍間隔の確率的振る舞いに関する情報を反映する前記スコアリング関数は、心拍間隔の確率分布を反映する確率分布を計算することにより決定する、
請求項１２に記載の方法。
【請求項１４】
心拍間隔の確率的振る舞いに関する情報を反映する前記スコアリング関数は、前に検出した心拍のタイミングを考慮することにより決定する、
請求項１２に記載の方法。
【請求項１５】
後続の心拍の時間推定値と前記バリストカルジオグラム信号とを用いて、前記心拍の自己同一性を改善するステップをさらに有する、
請求項１ないし１４いずれか一項に記載の方法。
【請求項１６】
時間推定を用いるステップは、
前記バリストカルジオグラム信号にハイパスフィルタをかけるステップと、
前記時間推定に関連する前記フィルタしたバリストカルジオグラム信号中のピークを特定するステップと、
心拍を、前記バリストカルジオグラム信号中の特定したピークに続く周波数は低いが振幅は大きい波として特定する、
請求項１５に記載の方法。
【請求項１７】
時間推定を用いるステップは、
時系列中の周期性を検出する方法を用いて、連続する２つの心拍間のビート間距離を特定するステップを有する、
請求項１５に記載の方法。
【請求項１８】
コンピュータまたはプロセッサで実行したとき、請求項１ないし１７いずれか一項に記載の方法のステップを実行するように構成された、コンピュータプログラムを含むコンピュータプログラム製品。
【請求項１９】
患者のバリストカルジオグラム信号を測定するデバイスと使用する装置であって、
前記デバイスからバリストカルジオグラム信号を受信する手段と、
前記受信したバリストカルジオグラム信号に、請求項１−１７いずれか一項に記載の方法を実行する処理手段とを有する、装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、バリストカルジオグラム信号を分析する方法と装置に関し、特に、バリストカルジオグラム信号中の単一の心拍イベントを検出する方法と装置に関する。
【背景技術】
【０００２】
バリストカルジオグラフ（ＢＣＧ）は、血液が心臓により送り出される時の血液の運動量による人体の動きを測定する。
【０００３】
ＢＣＧは、人体に電極を取り付けたり、ベルトや繊維などの特殊なセンサを身につけなくても、生体信号の測定が可能であるという点で、心電図（ＥＣＧ）より有利である。ＢＣＧは、睡眠品質、ストレス、または心臓性能を評価するために、脈拍数や脈拍数の変動データを調べるのに、特に有用である。これらは、控えめなＢＣＧモニタリングの性質が非常に重要なアプリケーションである。センサが患者と直接接触することにより、睡眠の質が悪くなることが避けられないからである。
【０００４】
現在、バリストカルジオグラム信号を分析して心拍数を決定するアルゴリズムは、例えば信号の自己相関関数を評価することにより、一定のパターンが繰り返すのを検出するスペクトル方法または時間ドメインにおける方法を用いる。これらのアプローチではすべて、複数の心拍をカバーするように、長さが数秒の信号のセグメントを考慮しなければならない。結果として、ある時間にわたる平均心拍を求めるが、正確なビート間情報（beat-to-beat  information）は得られない。
【０００５】
バリストカルジオグラム信号からビート間推定をするアルゴリズムが開示されている。しかし、これらは、正しく動作するには大きくて高価なセンサアレイを必要とするもの（"FFT  averaging  of  multichannel  BCG  signals  from  bed  mattress  sensor  to  improve  estimation  of  heart  beat  interval"  by  Kortelainen,  J.M.  and  Virkkala,  J.,  Engineering  in  Medicine  and  Biology  Society,  ２００７,  EMBS  ２００７,  ２９th  Annual  International  Conference  of  the  IEEE,  ２２-２６  August  ２００７,  pages  ６６８５-６６８８）、あるいは人間によるインターラクションをを必要とするもの（"Automatic  ballistocardiogram  (BCG)  Beat  Detection  Using  a  Template  Matching  Approach"  by  J.H.  Shin,  B.  H.  Choi,  Y.  G.  Lim,  D.U.  Joeng  and  K.S.  Park,  Engineering  in  Medicine  and  Biology  Society,  ２００８,  EMBS  ２００８,  ３０th  Annual  International  Conference  of  the  IEEE,  ２１-２４  August  ２００８）、あるいは異なるセンサモダリティを用いた精度が悪いもの（"Estimation  of  Respiratory  Waveform  and  Heart  Rate  Using  an  Accelerometer"  by  D.  H.  Phan,  S.  Bonnet,  R.  Guillemaud,  E.  Castelli,  N.Y.  Pham  Thi,  Engineering  in  Medicine  and  Biology  Society,  ２００８,  EMBS  ２００８,  ３０th  Annual  International  Conference  of  the  IEEE,  ２１-２４  August  ２００８）である。
【０００６】
これらのアルゴリズムを上市できるか、これらのアルゴリズムはバリストカルジオグラム信号の患者内の可変性や患者間の可変性を処理できるかについては疑問がある。
【０００７】
もちろん、ビート間推定は心電図（ＥＣＧ）からのデータを用いる標準的な方法で達成できるが、上記の通り、これらの方法は患者に配線をしなければならず、邪魔になる。
【０００８】
既述の通り、単一のセンサでキャプチャしたバリストカルジオグラム信号に基づいて心拍数を計算する技術水準の方法では、数秒の時間にわたる心拍推定しか求めることができない。睡眠ストレス、睡眠時無呼吸症候群、心不全などの分野への応用には、ビート間ベースの心拍数の変動に関する情報が必要である。すなわち、ある時間にわたる平均推定では不十分である。さらに、現在の方法は、正確な推定のためにはこの時間内で心拍が一定であることを要する。このため、異所性収縮や欠脈などの不整脈があると、心拍数の推定が乱れたり、単純に無視されたりする。したがって、脈拍数の変動の低周波数部分しか検出できない。
【発明の概要】
【課題を解決するための手段】
【０００９】
本発明の第１の態様による方法は、バリストカルジオグラム信号を分析して心拍数を決定する方法であって、前記バリストカルジオグラム信号中の第１の心拍の最初の時間推定値を決定する段階と、前記最初の時間推定値を用いて、前記バリストカルジオグラム信号中の後続の心拍の推定値を繰り返し計算する段階とを有し、前記計算するステップの各繰り返しは、複数のスコアリング関数の重み付けした和を含む目標関数の評価を含み、前記バリストカルジオグラム信号中の後続の心拍の推定値を計算するステップの各繰り返しは、その前の計算するステップの繰り返しで見つけた時間推定値後の目標区間に制限される。
【００１０】
本発明の第２の態様によるコンピュータプログラム製品は、コンピュータまたはプロセッサで実行したとき、上記の方法のステップを実行するように構成された、コンピュータプログラムを含むコンピュータプログラム製品である。
【００１１】
本発明の第３の態様による装置は、患者のバリストカルジオグラム信号を測定するデバイスと使用する装置であって、前記装置は、前記デバイスからバリストカルジオグラム信号を受信する手段と、前記受信したバリストカルジオグラム信号に上記の方法を実行する処理手段とを有する。
【００１２】
先行技術と対比して、本発明による方法と装置は、ビート間間隔をＥＣＧの場合と同様に計算する。それゆえ、上記のアプリケーションの必要性を満たすことができる。また、必要な信号を選るのに、単一のセンサのみを必要とする。高価なマルチセンサ機器や専門家による管理は必要ない。本方法とシステムは、一般的に非常に複雑なＢＣＧ信号に対してロバストに機能する。特徴的なイベント、長期的予測、及び心臓周期の長さに関するアプリオリな情報を用いるからである。このように、本方法と装置は、信号がマイナーな動きアーティファクトにより損なわれ、一方で不整脈のような不規則パターンを特定するくらい感度が高い状況を処理するにも十分にロバスト（robust）である。
【図面の簡単な説明】
【００１３】
添付した図面を参照して、一例としてここに本発明を説明する。
【図１】典型的なバリストカルジオグラム信号を示す図である。
【図２】本発明による方法を示すフローチャートである。
【図３】図２のステップ１０５を詳細に示す図である。
【図４】高周波成分に特徴的なスコアリング関数を示すグラフである。
【図５】バリストカルジオグラム信号のスペクトログラムである。
【図６】長期的予測をするためのスコアリング関数例を示すグラフである。
【図７】アプリオリな心拍の長さの確率密度関数のためのスコアリング関数例を示すグラフである。
【図８】合成関数とその最大値を示すグラフである。
【図９】図２のステップ１０５の結果を示すグラフである。
【図１０】図２のステップ１０５の結果とステップ１０７の改善した推定を示すグラフである。
【００１４】
〔発明の詳細な説明〕
図１は、ベッドや台に寝ている患者の胸部の下においた薄膜状センサを用いて取得した典型的なバリストカルジオグラム（ＢＣＧ）を示す。ＢＣＧは、呼吸による動きと、心拍数に関連する周期的なビートのパターンとを記録する。
【００１５】
図から分かるように、ＢＣＧは、患者の呼吸に関連する支配的な（約４秒の）低周波数成分と、心臓の機械的な動きによる、周波数が高く小さい変動とを含む。
【００１６】
言うまでもなく、このような明瞭なＢＣＧ信号を得るためには、患者は安静にしていなければならない。大きな動きがあると、ＢＣＧ信号では支配的な動きアーティファクトとなり、ＢＣＧ信号の分析を妨げることとなる。
【００１７】
ＢＣＧ信号を処理する最初のステップは、ＢＣＧ信号を、その他の動きや変動（perturbations）が心拍数や呼吸速度（breathing  rate）を邪魔する小区間（sub-intervals）と、推定が可能なエリアとに分割することである。この種の分割は信号のエネルギーレベルの評価により可能である。
【００１８】
さらに、呼吸による貢献と、心臓の機械的な動きによる貢献とは、フィルタを用いて分離できる。例えば、次数が３、低周波カットオフが０．０４ないし０．０８Ｈｚ、高周波カットオフが０．５０ないし０．７０Ｈｚのバターワースフィルタを用いたバンドパスフィルタリングにより、呼吸成分を得られる。心拍成分は、ハイパスフィルタ（例えば、次数が２、カットオフ周波数が０．８−１．２Ｈｚの範囲にあるバターワースフィルタ）を用いたフィルタリングにより得られる。以下に説明する方法は心拍成分を利用する。
【００１９】
本発明による方法を図２のフローチャートを参照して簡単に説明する。
【００２０】
ステップ１０１において、ＢＣＧにより信号を収集する。
【００２１】
バリストカルジオグラム信号の始めから第１の心拍発生時間の初期値をステップ１０３で計算する。この推定の計算は、以下にもっと詳しく説明するが、バリストカルジオグラム信号の終わりに到達するまで時間的に進みながら心拍の推定を繰り返し計算する（ステップ１０５に示した）繰り返し手順の開始点となる。ステップ１０５で実行する手順は、図３を参照してもっと詳しく説明する。
【００２２】
ステップ１０５の結果は、バリストカルジオグラム信号の心拍区間への第１のセグメント化である。
【００２３】
ステップ１０７において、ステップ１０５で求めた事前知識（foreknowledge）を用いて最終的な心拍区間の長さを計算する改善手順を実行する。この手続は、以下に詳しく説明する。
【００２４】
最初に、図３を参照して、ステップ１０３をもっと詳細に説明する。
【００２５】
心拍イベントが発生すると、次の心拍が一定の時間内に続くことが知られている。それゆえ、好ましい実施形態では、次の心拍を探す区間を生理的に合理的な値に制限することにより、この特性特性を利用する。このように、心拍推定ｔｎについて、次の心拍ｔｎ＋１は区間［ｔｎ＋ｔｍｉｎ，ｔｎ＋ｔｍａｘ］内になければならない。ｔｍｉｎとｔｍａｘの値としては、（心拍数が毎分約１２０拍を表す）約０．５秒ないし（心拍数が毎分約３０拍を表す）２秒がよいだろう。当業者には言うまでもないことだが、これとは別の生理的にありうる値を選択することもできる。
【００２６】
好ましくは、３つのスコアリング関数λ（ｔ）、μ（ｔ）、σ（ｔ）を生成し（ステップ１２３、１２５、１２７）、好ましくは特徴的な高周波成分の発生（スコアリング関数λ（ｔ））と心拍数の長期的振る舞い（μ（ｔ））と心拍の長さの確率分布（σ（ｔ））にそれぞれ関する基準に関するｔｎ＋１の推定がどのくらいよいか評価する。
【００２７】
ステップ１２９において、これらのスコアリング関数による最もよい推定が推定ｔｎ＋１となり、ｔｎ＋２の計算の基礎となる。ｔｎ＋２のサーチ区間は［ｔｎ＋１＋ｔｍｉｎ，ｔｎ＋１＋ｔｍａｘ］で与えられ、スコアリング関数λ（ｔ）、μ（ｔ）、σ（ｔ）が更新される。
【００２８】
この繰り返し手順を反復して、本方法は、バリストカルジオグラム信号の始めの初期値から開始して、バリストカルジオグラム信号の終わりになるまで、時間軸に沿ってバリストカルジオグラム信号を左から右に処理する。
【００２９】
次に、心拍推定ｔｎが与えられたときに、スコアリング関数λ（ｔ）、μ（ｔ）、σ（ｔ）をいかに計算するかと、スコアリング関数λ（ｔ）、μ（ｔ）、σ（ｔ）を用いて次の繰り返しステップｔｎ＋１をいかに計算するか、説明する。
【００３０】
特徴的高周波成分λ（ｔ）：
バリストカルジオグラム信号をバンドパスフィルタにかけて、高周波コンテンツを取り出す。得られた信号を二乗してローパスフィルタにかけると、支配的なピークを有する新しい信号が得られる。高周波成分は処理前の信号（raw  signal）にあったものである。これは、本願の出願人が本願と同時に出願した特許出願「Method  and  apparatus  for  the  analysis  of  ballistocardiogram  signals」でもっと詳しく説明した。この方法で生成した信号は、バリストカルジオグラム信号全体の特徴的高周波成分の発生のスコアリング関数であり、区間［ｔｎ＋ｔｍｉｎ，ｔｎ＋ｔｍａｘ］に制限すると、繰り返しにおけるこのステップのλ（ｔ）となる。図４は、バンドパスフィルタで２０ないし４０Ｈｚの周波数コンテンツを取り出し、ローパスフィルタのカットオフ周波数が３．５Ｈｚであるときのスコアリング関数λ（ｔ）の例を示す。
【００３１】
長期予測μ（ｔ）：
周波数が十分に集中した時間・周波数分布が関数μ（ｔ）の計算の開始点である。
【００３２】
スペクトログラムの絶対値を二乗すると、
【００３３】
【数１】
0005416218.tif 000002
ここで、ＢＣＧ（ｓ）は記録したバリストカルジオグラム信号を表し、ｈ（ｓ）は長さが１５秒のハニング窓を表す。
【００３４】
【数２】
0005416218.tif 000003
スペクトログラムは時刻ｔにおける周波数ｆの周波数コンテンツを測るものである。ｔｎと、時刻ｔにおける心拍イベントとにより、心拍数（ｔ−ｔｎ）−１が画定される。これが心拍数の正しい周波数と一致するなら、スペクトログラムの値はＳ（ｔ，（ｔ−ｔｎ）−１）で大きな値を示し、それ以外では小さい値となる。そのため、スコアリング関数μ（ｔ）を次の通りとする、
【００３５】
【数３】
0005416218.tif 000004
バリストカルジオグラム信号のスペクトログラムの一部を図５に示す。ここで、明るい部分はスペクトログラムのエネルギー値が高い部分であり、曲線は関数ｔ→（ｔ−ｔｎ）−１である。
【００３６】
図６は、スペクトログラムを用いた長期予測用のスコアリング関数μ（ｔ）を示す。図５中の曲線の関数値が、図６のスコアリング関数の計算に用いられる。
【００３７】
μ（ｔ）を決定する別の方法は、自己相関関数（ＡＣＦ）などの周期性決定方法の利用である。ここで、自己相関は、
【００３８】
【数４】
0005416218.tif 000005
により計算する。ここで、ＢＣＧはバリストカルジオグラム信号であり、２εは分析窓の長さであり、一般的には５ないし２０秒の範囲にある。
【００３９】
確率分布関数σ（ｔ）：
連続した２つの心拍間の間隔（interval  lengths）はガウス分布に従う。そのため、好ましい実施形態では、スコアリング関数σ（ｔ）を次式とする。
【００４０】
【数５】
0005416218.tif 000006
好ましい実施形態では、平均ｍと分散ｓの値はそれぞれｍ＝０．９２、ｓ＝０．４である。
【００４１】
別の実施形態では、前に検出した心拍の履歴を考慮してｍとｓを決めてもよい。あるいは、過去の心拍間隔に関する情報を用いて次の心拍間隔を予測する、当業者には知られた確率モデルを用いることもできる。
【００４２】
確率分布スコアリング関数σ（ｔ）の例を図７に示す。
【００４３】
３つの関数の重み付けをうまく制御するため、最大ノルムに対してスコアリング関数を規格化するか、適当なアフィン写像により関数の範囲を［０，１］区間に写像するとよい。次に、ステップ１２９において、αλ（ｔ）＋βμ（ｔ）＋χσ（ｔ）をｔ∈［ｔｎ＋ｔｍｉｎ，ｔｎ＋ｔｍａｘ］に関して最大化する。ここで、α、β、χはスカラーである。この多目的最適化問題を解決するため、この式の主要なピークを見つけなければならない。好ましくは、上記の和にローパスフィルタをかけてから、最大値を検索することにより、このピークを検出する。この様子を図８に示した。こうして発見したピークが心拍数ｔｎ＋１の次の推定となる。前述の通り、この推定から、信号の終わりに至るまで繰り返しを続ける。
【００４４】
言うまでもなく、ステップ１２３、１２５、及び１２７は、単独で用いると、誤差を生じ易い。例えば、個々の心拍に対するステップ１２３における高周波検索は、動きアーティファクトがあると（すなわち、患者が動き過ぎると）うまく機能しない。別の場合には、ステップ１２５における平均心拍数の計算は、例えば不整脈を適切に考慮していなかったり、まったく無視しているため、問題がある。最後に、ステップ１２７の確率的アプローチは、履歴データを入手できる場合にのみ有利である。
【００４５】
それゆえ、各ステップ１２３、１２５、及び１２７の出力をステップ１２９で１つの関数に合成して、パルスベースで心拍数のロバストな推定を行うために、この関数を最大化する。
【００４６】
再び図２を参照して、ここでステップ１０３を詳細に説明する。
【００４７】
図２の方法は、任意の開始点から始めて、少し繰り返すだけで安定し、正確な心拍数の推定を求めることができる。
【００４８】
あるいは、高周波成分は信号の始めに限定されているので、開始点をスコアリング関数の最大点とすることもできる。スコアリング関数を制限する区間は、［０ｓ，１.５ｓ］であってもよい。
【００４９】
さらに別の実施形態では、開始点ｔ１の推定と次の繰り返しステップｔ２の推定を同時に計算して、ｔ１が心拍セグメントのある開始点を表し、ｔ２が終了点を表すようにしてもよい。［０ｓ，１．５ｓ］の範囲にある各ｔ１について、Ｑ（ｔ１，ｔ２）＝αλ（ｔ２）＋βμ（ｔ２）＋χσ（ｔ２）を定義する。ｔ２∈［ｔ１＋ｔｍｉｎ，ｔ１＋ｔｍａｘ］である。言い換えると、Ｑ（ｔ１，ｔ２）は、上記の繰り返し方法に用いる同じスコアリング関数について、ｔ１とｔ２の選択がよいか測るものである。Ｑを最大化するｔ１、ｔ２は、最初の２つの心拍の時点ｔ１、ｔ２を最もよく推定するものである。
【００５０】
さらに別の実施形態では、上記の２つの方法のうち一方を用いて、予備的な開始点を計算してから、繰り返し点ｔｍ（例えば、ｍ＝１５）まで何回か通常の繰り返しを行う。しかし、そこから、アルゴリズムは「後ろ向きに」動く。すなわち、検索区間として、［ｔｍ＋ｔｍｉｎ，ｔｍ＋ｔｍａｘ］ではなく［ｔｍ−ｔｍａｘ，ｔｍ−ｔｍｉｎ］を選択する。心拍推定値の繰り返し計算を、信号の始めに至るまで、「後ろ向きに」続ける。「後ろ向き」の推定の最後の結果が開始点ｔ１となる。これは、本アルゴリズムは、開始点の選択が悪くても、数回の繰り返しの後に正しい推定値を与えることを用いている。
【００５１】
ステップ１０３と１０５を実行すると、バリストカルジオグラム信号の最初のセグメント化が得られる。各セグメントは、基本的には、１つの心拍のみの特徴的パターンを含んでいる。最終的に、このセグメント化を用いて改善方法の初期値を求め、これらの改善方法が正確な心拍イベントの時点や心拍間隔を求める。
【００５２】
一実施形態では、ステップ１０３と１０５により見つけた推定値に近い、バリストカルジオグラム信号中の特徴（characteristic  features）を見つけることにより、改善推定値を計算できる。好ましくはハイパスフィルタで処理したバリストカルジオグラム信号において、ステップ１０３、１０５で見つかった点の後には、周波数が低いが振幅は比較的大きい波が続く。ステップ１０３、１０５で見つかった点の小さい近傍にある局大値を検出することにより、心拍イベントをより正確に特定（localization）できる。心拍を特定してから、ビート間間隔を計算する。
【００５３】
さらに別の実施形態では、時系列における周期性を検出するために、自己相関（ＡＣＦ）その他の方法を用いて、連続する２つの心拍間のビート間隔を求める。この場合、ステップ１０３と１０５により、間隔に関するラフな情報が得られているが、ＡＣＦが最大値になる時間遅れΔｔを求めなければならない（すなわち、Δｔのすべての値を評価するわけではない）。Δｔの意味のある区間だけに限定するのは、ＡＣＦによる心拍数のビート間推定を可能とするために必要である。これにより、理論的に可能性のあるΔｔの検索空間全体を考慮する問題が解消される。ＡＣＦは、心拍イベントに関連しない偽の最大を頻繁に示す。
【００５４】
セグメント化手順の結果を図９に示す。図９には、（ハイパスフィルタ後の）バリストカルジオグラム信号を、特定したセグメントと、ＥＣＧ信号で求めた心拍イベント（×印で示した）とともに示した。セグメント結果後、次の０．１２秒以内に信号の最大値をマークすると、図１０に示したように、特徴的イベントがロバストに見つかる。
【００５５】
（改善なしと改善ありの）ビート間間隔推定値と、ＥＣＧとの比較を次の表１に示す。
【００５６】
【表１】
0005416218.tif 000007
当業者には、本方法の多数の意味のある構成と変形が理解できるだろう。ステップ１２９の目的関数に含まれる３つのスコアリング関数は、本方法を異なる必要性に合わせるため、異なる重み付けができる。これは、スカラーα、β、χを適当に選択することにより行える。例えば、規則的な心臓リズムが期待できる時は、ステップ１２５の長期的予測に重きをおくと、心拍推定値のロバスト性が増大する。一方、不整脈が検出された時は、ステップ１２３の高周波成分にフォーカスした方がよい。ステップ１２９の多目的関数の目的感覚ｔｍａｘ−ｔｍｉｎは、推定値と計算時間に明らかに影響し、適切に選択すべきである。
【００５７】
好ましい実施形態では、スカラーは次の値を有する：α＝１、β＝０．６、及びχ＝０．４。これにより、よい汎用心拍推定ができる。目標間隔（すなわち、ｔｍａｘ−ｔｍｉｎ）は計算時間を節約するため約１．２秒とすべきである。
【００５８】
別の実施形態では、本方法を不整脈のある患者に利用する場合、例えば、スカラーは次の値としてもよい：α＝１、β＝０（長期予測を実質的に無効にするか、もしくは無視する）、及びχ＝０．２不整脈を正しく特定するために、目標間隔は一般的な場合よりも少し長く設定できる（すなわち、ｔｍａｘ−ｔｍｉｎ＝２）。
【００５９】
別のスコアリング関数にそれぞれスカラーをかけたものを加えることにより、多目的関数αλ（ｔ）＋βμ（ｔ）＋χσ（ｔ）を拡張するのは、より多くの情報源を利用するときに便利である。
【００６０】
また、本方法は、長期的予測と計算に必要な時間だけ遅れるが、（ほぼ）リアルタイムで心拍を推定するために用いることもできる。
【００６１】
上記の方法により、バリストカルジオグラム信号からビート間間隔を取得することもできる。上記の方法を用いたビート間推定により、心不全管理、不整脈検出、心房細動の診断や管理など様々なアプリケーションで標準的なＥＣＧ装置を置き換えることもできる。さらに、パルス間分析により、正確に心拍数変動を計算できるようになる。これは睡眠及びストレスの分析に必要である。不整脈検出をできるという潜在的能力のため、ローアキュイティモニタリング（low  acuity  monitoring）のために、専門病院においてＢＣＧ法は受け入れられつつある。これは注目すべきことである。集中治療室（ＩＣＵ）のベッド数は限られており、一般病棟用の簡易で低費用のモニタリングソリューションが望まれているからである。
【００６２】
結局、上記のアルゴリズムにより、場合によっては他のセンサモダリティに付随して、心原性失神に起因する転倒を予測して、回避することができる。
【００６３】
本発明は、養護施設、病院、自宅介護監視に用いることができる。どの場合にも、ＥＣＧに対するＢＣＧの一般的な利点は、患者に電極などを取り付けなくても、目立たずに患者をモニタできることにある。
【００６４】
本発明を方法またはアルゴリズムにより説明したが、言うまでもなく、本発明をＢＣＧシステム（すなわち、ＢＣＧ信号を測定する装置と結合したコンピュータ装置）、またはスタンドアロンのコンピュータシステム、またはプログラムとして実施することができる。言うまでもなく、ＢＣＧシステムは、本発明の装置に、バリストカルジオグラム信号をアナログ形式またはデジタル形式で供給でき、本発明の装置はこの信号を受信するように適宜構成できる。例えば、ＢＣＧシステムは本装置にバリストカルジオグラム信号をアナログ形式で供給し、本装置は、その適宜プログラムされたデジタルシグナルプロセッサに、バリストカルジオグラム信号をデジタル化した信号を供給するアンチエイリアスフィルタ（anti-alias  filter）とアナログ・デジタル・コンバータを有する。あるいは、ＢＣＧシステムは、アナログ・デジタルコンバータを実装し、バリストカルジオグラム信号をデジタル形式で本装置に（そして、具体的には本装置のデジタルシグナルプロセッサに）供給してもよい。本装置は、ＢＣＧシステムへの有線または無線の接続などの適当な手段により、バリストカルジオグラム信号を受信できる。
【００６５】
本発明を、図面と上記の説明に詳しく示し説明したが、かかる例示と説明は例であり限定ではなく、本発明は開示した実施形態には限定されない。
【００６６】
請求項に記載した発明を実施する際、図面、本開示、及び添付した特許請求の範囲を研究して、開示した実施形態のバリエーションを、当業者は理解して実施することができるであろう。請求項において、「有する（comprising）」という用語は他の要素やステップを排除するものではなく、「１つの（"a"  or  "an"）」という表現は複数ある場合を排除するものではない。
【００６７】
単一のプロセッサまたはその他のアイテムが請求項に記載した複数のユニットの機能を満たすこともできる。相異なる従属クレームに手段が記載されているからといって、その手段を組み合わせて有利に使用することができないということではない。請求項に含まれる参照符号は、その請求項の範囲を限定するものと解してはならない。コンピュータプログラムは、光記憶媒体や他のハードウェアとともに、またはその一部として供給される固体媒体などの適切な媒体に記憶／配布することができ、インターネットや有線または無線の電気通信システムなどを介して他の形式で配信することもできる。
【図１】
0005416218.tif 000008
【図２】
0005416218.tif 000009
【図３】
0005416218.tif 000010
【図４】
0005416218.tif 000011
【図５】
0005416218.tif 000012
【図６】
0005416218.tif 000013
【図７】
0005416218.tif 000014
【図８】
0005416218.tif 000015
【図９】
0005416218.tif 000016
【図１０】
0005416218.tif 000017