(57)【要約】【課題】心臓の機能不全と異常を表す情報を測定する。【解決手段】心臓の機能不全と異常を表す情報を測定するための器械は心臓血管の動きを表す信号から振幅の変化を検出する、ように構成されている処理装置（４０２）を含み、処理装置は心臓活動に関連する電磁現象を表す信号から心拍数を繰り返す第１波形を抽出し、且つ、前記第２信号から前記心拍数を繰り返す第２波形を抽出するように構成されている。処理装置は各タイミング値がある１つの心拍周期に属する前記第１波形の基準点から同一の心拍周期に属する前記第２波形の基準点までの時間を表すタイミングデータを作成する。処理装置は、少なくとも部分的にタイミングデータに基づき、心臓の機能不全と異常を表す情報を測定する。【選択図】図４
2018086278.tif 000002
【特許請求の範囲】
【請求項１】
心臓活動に関連する電磁現象を表す第１信号と心臓血管の動きを表す第２信号を受信するための信号インターフェース（４０１）と、
前記信号インターフェースに接続されており、前記第１信号からある心拍数を繰り返す第１波形を抽出し、且つ、前記第２信号から前記心拍数を繰り返す第２波形を抽出するように構成されている処理装置（４０２）と、
を備える器械であって、
前記処理装置が、さらに、
−各タイミング値がある１つの心拍周期に属する前記第１波形の基準点から同一の心拍周期に属する前記第２波形の基準点までの時間を表すタイミングデータを作成し、且つ、
−前記タイミングデータと前記心拍数を表すペーシングデータとの間の相関を決定する、
ように構成されており、
前記相関が心臓の機能不全と異常を表す、
ことを特徴とする前記器械。
【請求項２】
前記器械が、さらに、
前記第１信号を測定するための第１センサ素子（４０３）と、
前記第２信号を測定するための第２センサ素子（４０４）と、
を備える、ことを特徴とする請求項１に記載の器械。
【請求項３】
前記第２センサ素子が、加速度センサ、圧電センサ、傾斜センサ、圧力センサのうちの１つを備える、ことを特徴とする請求項２に記載の器械。
【請求項４】
前記処理装置が、
心室性筋肉組織の脱分極によって生じ、前記第１信号上で前記心拍数を繰り返す前記第１波形を表すＲピークを前記第１信号から抽出し、且つ、
大動脈弁の開放によって生じ、前記第２信号上で前記心拍数を繰り返す前記第２波形を表すＡＯピークを前記第２信号から抽出する、
ように構成されている、
ことを特徴とする請求項１〜３のいずれか一項に記載の器械。
【請求項５】
前記処理装置が、
心室性筋肉組織の脱分極によって生じ、前記第１信号上で前記心拍数を繰り返す前記第１波形を表すＲピークを前記第１信号から抽出し、且つ、
大動脈弁を通る左心室から大動脈への血流によって生じ、前記第２信号上で前記心拍数を繰り返す前記第２波形を表すＪピークを前記第２信号から抽出する、
ように構成されている、
ことを特徴とする請求項１〜３のいずれか一項に記載の器械。
【請求項６】
前記処理装置が、
心室性筋肉組織の脱分極によって生じ、前記第１信号上で前記心拍数を繰り返す前記第１波形を表すＲピークを前記第１信号から抽出し、且つ、
僧帽弁の閉鎖によって生じるＭＣピークと大動脈弁の開放によって生じるＡＯピークの間の谷であって、前記第２信号上で前記心拍数を繰り返す前記第２波形を表すＭＣピークとＡＯピークの間の前記の谷を前記第２信号から抽出する、
ように構成されている、
ことを特徴とする請求項１〜３のいずれか一項に記載の器械。
【請求項７】
前記処理装置が、
Ｒピークの頂点を前記第１波形の前記基準点として使用し、且つ、
ＡＯピークの頂点、Ｊピークの頂点、又はＭＣピークとＡＯピークの間の前記の谷の最低点を前記第２波形の前記基準点として使用する、
ように構成されている、
ことを特徴とする請求項４〜６のいずれか一項に記載の器械。
【請求項８】
前記処理装置が、
前記タイミングデータと前記ペーシングデータとの間の前記相関の程度を表す相関係数：
【数１】
2018086278.tif 000011
を計算するように構成されており、
式中、
Ｃ（ｊ）が前記相関係数であり、
Ｅが期待値演算子であり、
ＴＤが前記タイミングデータであり、
μTが前記タイミングデータの平均であり、
ＰＤが前記ペーシングデータであり、
μPが前記ペーシングデータの平均であり、
ｊが心拍周期内の前記タイミングデータに関して前記ペーシングデータの時間差を表す整数である、
ことを特徴とする請求項１に記載の器械。
【請求項９】
前記処理装置が、
実質的に一定の心拍数に対応する前記タイミングデータの変化を測定するように構成されており、
前記の変化の程度が心臓の機能不全と異常を表す、
ことを特徴とする請求項１〜８のいずれか一項に記載の器械。
【請求項１０】
前記処理装置が前記タイミングデータの前記の変化の程度を表す変化量：
【数２】
2018086278.tif 000012
を測定するように構成されており、
式中、
Ｖが前記変化量であり、
Ｅが期待値演算子であり、
ＴＤが前記タイミングデータであり、
μTが検討中の前記心拍数における前記タイミングデータの平均である、
ことを特徴とする請求項９に記載の器械。
【請求項１１】
前記処理装置が前記ペーシングデータを得るために前記第１信号上の連続するＲピーク間の時間を検出するように構成されており、
前記の検出された時間が前記心拍数に反比例するものとして前記の心拍数を表す、
ことを特徴とする請求項１又は８に記載の器械。
【請求項１２】
前記処理装置が、
前記ペーシングデータの各値が前記心拍数の瞬時値を表すように前記ペーシングデータを作成し、且つ、
前記相関係数が閾値よりも大きい状況に応答して心房細動を表す信号を発信する、
ように構成されている、
ことを特徴とする請求項８に記載の器械。
【請求項１３】
前記処理装置が、
前記ペーシングデータの各値が１回の心拍周期の時間の長さを表すように前記ペーシングデータを作成し、且つ、
前記相関係数が閾値未満である状況に応答して心房細動を表す信号を発信する、  
ように構成されている、
ことを特徴とする請求項８に記載の器械。
【請求項１４】
前記処理装置が０を前記の閾値として使用する、
ように構成されている、ことを特徴とする請求項１２又は１３に記載の器械。
【請求項１５】
前記処理装置が、
前記の心臓血管の動きを表す前記第２信号をローパスフィルタにかけ、且つ、
前記のローパスフィルタにかけられた信号からＡＯピークを検出する、
ように構成されている、
ことを特徴とする請求項４〜１４のいずれか一項に記載の器械。
【請求項１６】
前記ローパスフィルタリングの上限の周波数が３０Ｈｚである、ことを特徴とする請求項１５に記載の器械。
【請求項１７】
前記処理装置が、
前記の心臓血管の動きを表す前記第２信号をバンドパスフィルタにかけ、且つ、
大動脈弁の閉鎖によって生じるＡＣピークを前記のバンドパスフィルタにかけられた信号から検出する、
ように構成されている、請求項４〜１６のいずれか一項に記載の器械。
【請求項１８】
前記バンドパスフィルタリングの通過帯域が４０Ｈｚから１００Ｈｚまでである、ことを特徴とする請求項１７に記載の器械。
【請求項１９】
前記処理装置が、
ＡＣ〜ＡＯ間隔の時間の長さを検出し、且つ、
前記の検出された前記ＡＣ〜ＡＯ間隔の時間の長さの変化の強度を表す時間変化量を計算する、
ように構成されており、
前記ＡＣ〜ＡＯ間隔のそれぞれが前記ＡＣピークのうちの１つから前記ＡＯピークのうちの次の１つまでの時間間隔であり、
前記時間変化量が心臓の機能不全と異常を表す、
ことを特徴とする請求項１７又は１８に記載の器械。
【請求項２０】
前記処理装置が、
ＡＣ〜ＡＯ間隔の時間の長さとＡＯ〜ＡＯ間隔の時間の長さを検出し、且つ、
同一の心拍周期内の前記のＡＣ〜ＡＯ間隔の時間の長さと前記のＡＯ〜ＡＯ間隔の時間の長さの間の比率を表す第１比率量を計算する、
ように構成されており、
前記ＡＣ〜ＡＯ間隔のそれぞれが前記ＡＣピークのうちの１つから前記ＡＯピークのうちの次の１つまでの時間間隔であり、
前記ＡＯ〜ＡＯ間隔のそれぞれが前記ＡＯピークのうちの１つから前記ＡＯピークのうちの次の１つまでの時間間隔であり、
前記第１比率量が心臓の機能不全と異常を表す、
ことを特徴とする請求項１７又は１８に記載の器械。
【請求項２１】
前記処理装置が、
ＡＣ〜Ｒ間隔の時間の長さとＲ〜Ｒ間隔の時間の長さを検出し、且つ、
同一の心拍周期内の前記のＡＣ〜Ｒ間隔の時間の長さと前記のＲ〜Ｒ間隔の時間の長さの間の比率を表す第２比率量を計算する、
ように構成されており、
前記ＡＣ〜Ｒ間隔のそれぞれが前記ＡＣピークのうちの１つから前記Ｒピークのうちの次の１つまでの時間間隔であり、
前記Ｒ〜Ｒ間隔のそれぞれが前記Ｒピークのうちの１つから前記Ｒピークのうちの次の１つまでの時間間隔であり、
前記第２比率量が心臓の機能不全と異常を表す、
ことを特徴とする請求項１７又は１８に記載の器械。
【請求項２２】
心臓活動に関連する電磁現象を表す第１信号からある心拍数を繰り返す第１波形を抽出するステップ（２０１）と、
心臓血管の動きを表す第２信号から前記心拍数を繰り返す第２波形を抽出するステップ（２０２）と、
を含む方法であって、
前記方法が、
−各タイミング値がある１つの心拍周期に属する前記第１波形の基準点から同一の心拍周期に属する前記第２波形の基準点までの時間を表すタイミングデータを作成するステップ（２０３）と、
−前記タイミングデータと前記心拍数を表すペーシングデータとの間の相関を決定するステップ（２０４）と、
をさらに含み、
前記相関が心臓の機能不全と異常を表す、
ことを特徴とする方法。
【請求項２３】
前記方法が、
心室性筋肉組織の脱分極によって生じ、前記第１信号上で前記心拍数を繰り返す前記第１波形を表すＲピークを前記第１信号から抽出するステップと、
及び大動脈弁の開放によって生じ、前記第２信号上で前記心拍数を繰り返す前記第２波形を表すＡＯピークを前記第２信号から抽出するステップと
を含む、請求項２２に記載の方法。
【請求項２４】
前記方法が、
心室性筋肉組織の脱分極によって生じ、前記第１信号上で前記心拍数を繰り返す前記第１波形を表すＲピークを前記第１信号から抽出するステップと、
大動脈弁を通る左心室から大動脈への血流によって生じ、前記第２信号上で前記心拍数を繰り返す前記第２波形を表すＪピークを前記第２信号から抽出するステップと、
を含む、
ことを特徴とする請求項２２に記載の方法。
【請求項２５】
前記方法が、
心室性筋肉組織の脱分極によって生じ、前記第１信号上で前記心拍数を繰り返す前記第１波形を表すＲピークを前記第１信号から抽出するステップと、
僧帽弁の閉鎖によって生じるＭＣピークと大動脈弁の開放によって生じるＡＯピークの間の谷であって、前記第２信号上で前記心拍数を繰り返す前記第２波形を表すＭＣピークとＡＯピークの間の前記の谷を前記第２信号から抽出するステップと、
を含む、  
ことを特徴とする請求項２２に記載の方法。
【請求項２６】
Ｒピークの頂点が前記第１波形の前記基準点であり、且つ、
ＡＯピークの頂点、Ｊピークの頂点、又はＭＣピークとＡＯピークの間の前記の谷の最低点が前記第２波形の前記基準点である、
ことを特徴とする請求項２３〜２５のいずれか一項に記載の方法。
【請求項２７】
前記タイミングデータと前記ペーシングデータとの間の前記相関の程度を表す相関係数：
【数３】
2018086278.tif 000013
が計算され、
式中、
Ｃ（ｊ）が前記相関係数であり、
Ｅが期待値演算子であり、
ＴＤが前記タイミングデータであり、
μTが前記タイミングデータの平均であり、
ＰＤが前記ペーシングデータであり、
μPが前記ペーシングデータの平均であり、
ｊが心拍周期内の前記タイミングデータに関して前記ペーシングデータの時間差を表す整数である、
ことを特徴とする請求項２２に記載の方法。
【請求項２８】
前記心臓の機能不全と異常の指標を測定するステップ（２０４）が実質的に一定の心拍数に対応する前記タイミングデータの変化を測定するステップを含み、
前記の変化の程度が心臓の機能不全と異常を表す、
ことを特徴とする請求項２２〜２６のいずれか一項に記載の方法。
【請求項２９】
前記タイミングデータの前記の変化の程度を表す変化量：
【数４】
2018086278.tif 000014
が計算され、
式中、
Ｖが前記変化量であり、
Ｅが期待値演算子であり、
ＴＤが前記タイミングデータであり、
μTが検討中の前記心拍数における前記タイミングデータの平均である、
ことを特徴とする請求項２８に記載の方法。
【請求項３０】
前記方法が、
前記ペーシングデータを得るために前記第１信号上の連続するＲピーク間の時間を検出するステップを含み、
前記の検出された時間が前記心拍数に反比例するものとして前記の心拍数を表す、
ことを特徴とする請求項２２又は２７に記載の方法。
【請求項３１】
前記ペーシングデータの各値が前記心拍数の瞬時値を表し、且つ、
前記方法が、前記相関係数が閾値よりも大きい状況に応答して心房細動を表す信号を発信するステップを含む、
ことを特徴とする請求項２７に記載の方法。
【請求項３２】
前記ペーシングデータの各値が１回の心拍周期の時間の長さを表し、且つ、
前記方法が、前記相関係数が閾値未満である状況に応答して心房細動を表す信号を発信するステップを含む、
ことを特徴とする請求項２７に記載の方法。
【請求項３３】
前記閾値が０である、請求項３１又は３２に記載の方法。
【請求項３４】
前記方法が、
前記の心臓血管の動きを表す前記第２信号をローパスフィルタにかけるステップ（２１１）と、
前記のローパスフィルタにかけられた信号からＡＯピークを検出するステップ（２１２）と、
を含む、ことを特徴とする請求項２３〜３３のいずれか一項に記載の方法。
【請求項３５】
前記ローパスフィルタリングの上限の周波数が３０Ｈｚである、ことを特徴とする請求項３４に記載の方法。
【請求項３６】
前記方法が、
前記の心臓血管の動きを表す前記第２信号をバンドパスフィルタにかけるステップ（２１１）と、
大動脈弁の閉鎖によって生じるＡＣピークを前記のバンドパスフィルタにかけられた信号から検出するステップ（２１２）と、
を含む、ことを特徴とする請求項２３〜３５のいずれか一項に記載の方法。
【請求項３７】
前記バンドパスフィルタリングの通過帯域が４０Ｈｚから１００Ｈｚまでである、ことを特徴とする請求項３６に記載の方法。
【請求項３８】
前記方法が、
ＡＣ〜ＡＯ間隔の時間の長さを検出するステップと、
前記の検出された前記ＡＣ〜ＡＯ間隔の時間の長さの変化の強度を表す時間変化量を計算するステップと、を含み、
前記ＡＣ〜ＡＯ間隔のそれぞれが前記ＡＣピークのうちの１つから前記ＡＯピークのうちの次の１つまでの時間間隔であり、
前記時間変化量が心臓の機能不全と異常を表す、
ことを特徴とする請求項３６又は３７に記載の方法。
【請求項３９】
前記方法が、
ＡＣ〜ＡＯ間隔の時間の長さとＡＯ〜ＡＯ間隔の時間の長さを検出するステップと、
同一の心拍周期内の前記のＡＣ〜ＡＯ間隔の時間の長さと前記のＡＯ〜ＡＯ間隔の時間の長さの間の比率を表す第１比率量を計算するステップと、を含み、
前記ＡＣ〜ＡＯ間隔のそれぞれが前記ＡＣピークのうちの１つから前記ＡＯピークのうちの次の１つまでの時間間隔であり、
前記ＡＯ〜ＡＯ間隔のそれぞれが前記ＡＯピークのうちの１つから前記ＡＯピークのうちの次の１つまでの時間間隔であり、
前記第１比率量が心臓の機能不全と異常を表す、
ことを特徴とする請求項３６又は３７に記載の方法。
【請求項４０】
前記方法が、
ＡＣ〜Ｒ間隔の時間の長さとＲ〜Ｒ間隔の時間の長さを検出するステップと、
同一の心拍周期内の前記ＡＣ〜Ｒ間隔の時間の長さと前記Ｒ〜Ｒ間隔の時間の長さの間の比率を表す第２比率量を計算すステップと、を含み、
前記ＡＣ〜Ｒ間隔のそれぞれが前記ＡＣピークのうちの１つから前記Ｒピークのうちの次の１つまでの時間間隔であり、
前記Ｒ〜Ｒ間隔のそれぞれが前記Ｒピークのうちの１つから前記Ｒピークのうちの次の１つまでの時間間隔であり、
前記第２比率量が心臓の機能不全と異常を表す、
ことを特徴とする請求項３６又は３７に記載の方法。
【請求項４１】
コンピュータプログラムであって、
心臓活動に関連する電磁現象を表す第１信号からある心拍数を繰り返す第１波形を抽出し、且つ、心臓血管の動きを表す第２信号から前記心拍数を繰り返す第２波形を抽出する、プログラミング可能なプロセッサを制御するためのコンピュータが実行可能な命令を含み、
さらに、
−各タイミング値がある１つの心拍周期に属する前記第１波形の基準点から同一の心拍周期に属する前記第２波形の基準点までの時間を表すタイミングデータを作成し、且つ、
−前記タイミングデータと前記心拍数を表すペーシングデータとの間の相関を決定し、
前記相関が心臓の機能不全と異常を表す、プログラミング可能なプロセッサを制御するためのコンピュータが実行可能な命令を含む、
ことを特徴とするコンピュータプログラム。
【請求項４２】
請求項４１に記載のコンピュータプログラムが符号化されている非一時的コンピュータ可読媒体を備える、ことを特徴とするコンピュータプログラム製品。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は概して心臓の機能不全と異常、例えば、心房細動などを表す情報を測定することに関する。
より具体的には、本発明は心臓の機能不全と異常を表す情報を測定するための器械と方法に関する。
さらに、本発明は心臓の機能不全と異常を表す情報を測定するためのコンピュータプログラムに関する。
【背景技術】
【０００２】
心臓血管系において生じ得る機能不全と異常は、診断されず、適切に治療又は改善されなかった場合、特に、個体がストレスに出会った時に冠動脈の酸素需要を満足させる充分な酸素を供給する心臓血管系の能力を漸進的に低下させ得る。
現在では、心臓活動に関連する電磁現象に基づく心拍記録法、心エコー検査法、及び心臓血管の動きに基づく心拍記録法などの方法が様々な心臓の機能不全と異常の特定と評価において用いられている。
心臓活動に関連する電磁現象に基づく心拍記録法のよく知られている例は心電図記録法「ＥＣＧ」であり、心臓血管の動きに基づく心拍記録法の例はバリストカルジオグラフィー「ＢＣＧ」及びサイスモカルジオグラフィー「ＳＣＧ」である。
心エコー検査法は心臓の切片の画像を提供し、且つ、心臓の構造と機能についての包括的な情報を提供することができるが、高価な機器と特別な操作担当者を必要とする。
ＥＣＧは心臓の電気的評価をかなり速く与えるが、収縮力に関連するどのような情報も提供しない。
心臓血管の動きに基づく心拍記録法は心臓血管の動きを表す信号の測定を伴う。
以前はその信号は、動きを測定するための器械が用意されているベッドに個体が横たわっている間に、又は、脚の脛のいたる所に取り付けられた支援装置が存在する間に取得された。
現在では、その信号は、心臓の動きを表す微細な動きを測定するのに適切である小さなセンサ素子、例えば加速度センサを使用して取得され得る。
【０００３】
図１ａ、図１ｂは、律動的な電気的機能と関連の心臓血管の動きとの間の関係を示す。
図１ａはＥＣＧ波形の例を示し、図１ｂは、心臓血管の動きを表し、通常はｙ方向と呼ばれる「頭足」方向で加速度センサによって測定された例となる信号の波形を示す。
例示を目的として基本的な心臓機能の簡単な説明を下に提供する。
【０００４】
心臓は４つの小室を含む。
右心房は三尖弁によって右心室と相互接続されており、左心房は僧帽弁によって左心室によって相互接続されている。
血液は上大静脈を介して上半身から右心房に、下大静脈を介して下半身から右心房に送られる。
右心房の心筋と右心室の乳頭筋が同時に収縮して三尖弁が開き、それによって血液が右心房から右心室へ流れる。
その後、その乳頭筋が弛緩して三尖弁が閉じる。
右心室の心筋が収縮すると、血液が肺動脈弁を通って右心室から血液を肺に送る肺動脈へ押し出され、そこで血液は酸素負荷される。
次に酸素負荷された血液は肺静脈を介して左心房へ送られる。
酸素負荷された血液は、血液の左心房から左心室へ流れを可能にする左心房の心筋と左心室の乳頭筋の同時収縮により僧帽弁が開いたときに左心房から左心室へ流れる。
その後、その乳頭筋が弛緩して僧帽弁が閉じる。
次に酸素負荷された血液は大動脈弁を通って左心室から酸素負荷された血液を末梢血管系に送る大動脈へ押し出される。
【０００５】
各心拍周期は心房収縮期、心室収縮期及び心臓拡張期の主要な３つのステージを伴う。
心房収縮期は右心房と左心房を取り囲む心筋の収縮期である。
両方の心房が乳頭筋の収縮により同時に収縮してそれにより三尖弁と僧帽弁を開放させる。
心臓の小室の筋肉組織を刺激してそれらを収縮させる電気的活性、すなわち、電気的収縮は右心房に位置する洞房結節において始まる。
電気的脱分極の伝達が続いて、各心房筋細胞を次々と脱分極させる両方の心房を介して波として下方、左方、及び後方へ伝わる。
この電荷の伝播を図１ａに示されるＥＣＧ波形上のＰ波として見ることができる。
この直後に、図１ｂに示される波形のｈピーク（ｐｅａｋ）に対応するインパクトと図１ｂに示される波形のｉバレー（ｖａｌｌｅｙ）に対応する反動として検出される心房の機械的収縮が続く。
右心房と左心房が収縮し始めると、右心室と左心室への高速の血流が生じ、それは図１ｂに示される波形上のｊピークにより表される。
三尖弁が閉じ始めると、続いている心房の収縮によってさらに右心室と左心室への低速の血流が生じる。
その追加的な血流は「心房キック」と呼ばれ、図１ｂに示される波形の「ａ−ａ1」波複合に対応する。
心房が空になった後に三尖弁と僧帽弁が閉じ、それにより図１ｂに示される波形上の下行ｇ波が生じる。
心室収縮は左心室と右心室の筋肉の収縮であり、図１ａに示されるＥＣＧ波形の「Ｑ−Ｒ−Ｓ」波複合を生じる心室心筋層の電気的脱分極によって引き起こされる。
下行Ｑ波は、「ヒス束」と呼ばれる特別な細胞群に沿う中隔を通る脱分極の下降流によって引き起こされる。
Ｒピークは心室筋肉組織の脱分極によって生じ、Ｓ波は心房と心室の間の心臓組織の脱分極によって生じる。
脱分極が中隔を下り、心室心筋層中に伝わると、心房と洞房結節が分極し始める。
三尖弁と僧帽弁の閉鎖が心室収縮期の始まりであり、心臓が拍動するときに心臓によって生じるラブ・ダブ（ｌｕｂ−ｄｕｂ）音の１番目の部分を生じる。
この音は通常は「１音」として知られている。心室心筋層の電気的脱分極がピークに達すると、右心室と左心室を分ける房室「ＡＶ」中隔が収縮し、図１ｂに示される波形上のＨピークに対応するインパクトと図１ｂに示される波形上のＩバレーに対する反動を引き起こす。
心室の収縮により血液が肺動脈弁を介して右心室から肺動脈へ、及び大動脈弁を介して左心室から大動脈へ非常に高速で押し出され、それにより図１ｂに示される波形上のＪピークが生じる。
左心室から大動脈への血流の減速により、図１ｂに示される波形上の下行Ｋ波が生じる。
左心室が空になるとその圧力は大動脈の圧力よりも下に落ち、大動脈弁が閉じる。
同様に、右心室の圧力が肺動脈の圧力よりも下に落ちると肺動脈弁が閉じる。
ラブ・ダブ音の２番目の部分は通常は「２音」として知られており、心室収縮期の終わりの肺動脈弁と大動脈弁の閉鎖によって生じ、それにより図１ｂに示される波形上の上行Ｌ波が生じる。
肺動脈弁と大動脈弁の閉鎖と同時に房室「ＡＶ」中隔が弛緩及び上昇し、心室心筋層が再分極して図１ａに示されるＥＣＧ波形上のＴ波が生じる。
心臓拡張期は心房拡張期と心室拡張期を含み、心臓が収縮後に弛緩し、且つ、循環血液で再補充される準備をしている期間である。
心房拡張期は右心房と左心房が弛緩している期間であり、心室拡張期は右心室と左心室が弛緩している期間である。
心房拡張期の間に右心房は脱酸素化された血液によって再補充され、左心房は酸素負荷された血液によって再補充される。
心房の再補充によって図１ｂに示される波形上の下行Ｍ波がヒス束細胞の再分極と同時に起こる拡張期の初期に生じ、その再分極はＥＣＧ波形のＵ波として示される。
右心房と左心房がそれらの最大容量まで補充されると三尖弁と僧帽弁への血液の逆流により、図１ｂに示される波形上の上行Ｎ波が生じる。
【０００６】
特許文献１は、前胸部の加速血流シグナルを記録することにより対象の心収縮力を評価するための方法を記載する。
特許文献２は個体の生理的状態をモニターし、その異常を検出するための方法を記載する。
その方法は、ＥＣＧ信号である第１信号と心臓血管の動きを表す第２信号を同時に受信することを含む。
【０００７】
心臓血管の動きを表す波形の分析は、異常な心臓血管機能を正常な事例と区別するために通常は有資格の診断医によって視覚的に実施される。
しかしながら、多くの場合では、例えば、心房細動などのある特定の心臓の機能不全と異常を視覚的分析により見つけ出すことは困難であり得る。
したがって、心臓の機能不全と異常を表す情報を測定するための方法と器械の必要性が存在する。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】国際公開第２０１２１４９６５２号パンフレット
【特許文献２】米国特許出願公開第２００８１９４９７５号明細書
【発明の概要】
【０００９】
次に、発明の様々な実施形態のうちの幾つかの態様を基本的に理解してもらうために簡略化した概要を提示する。
本概要は本発明の広範囲にわたる概説ではない。
それは本発明の重要又は重大な要素を特定することを意図しても、本発明の範囲を詳しく説明することを意図してもいない。
次の概要は、本発明の好ましい実施形態のさらに詳しい説明の前置きとなる本発明の簡略化された幾つかの構想を提示するだけのものである。
【００１０】
本発明により、心臓の機能不全と異常、例えば、心房細動を表す情報を測定するための新しい方法が提供される。
本発明による方法は、
−心臓活動に関連する電磁現象を表す第１信号からある心拍数を繰り返す第１波形を抽出するステップと、
−心臓血管の動きを表す第２信号からその心拍数を繰り返す第２波形を抽出するステップと、
−各タイミング値がある１つの心拍周期に属する第１波形の基準点から同一の心拍周期に属する第２波形の基準点までの時間を表すタイミングデータを作成するステップと、
−前記タイミングデータと前記の心拍数を表すペーシングデータとの間の相関を決定するステップと、
を含み、
その相関が心臓の機能不全と異常を表す。
【００１１】
第２波形が第１信号上の第１波形に対する心臓の応答を表すものであるように第２波形が選択されることが好都合である。
第１信号は、例えば、心電計「ＥＣＧ」波形又は誘導的に測定された波形を表し得る。
第１波形は、必ずというわけではないが、例えば、図１ａに示されるＥＣＧ波形上の「Ｑ−Ｒ−Ｓ」波複合のＲピークであり得、第２波形は、必ずというわけではないが、例えば、図１ｂに示される波形上のＪピークであり得る。
この場合、Ｒピークの頂点は第１波形の基準点として使用され得、Ｊピークの頂点は第２波形の基準点として使用され得、各タイミング値はＲピークの頂点の瞬間からＪピークの頂点の瞬間までの時間を表し得る。
別の例では、第２波形は通常はｚ方向と呼ばれる「通胸」方向で測定された心臓血管の動きを表す波形上の大動脈弁開放「ＡＯ」ピークであり得る。
【００１２】
本方法は、心臓活動に関連する電磁現象を表す第１信号及び／又は心臓血管の動きを表す第２信号から上述のペーシングデータを検出するステップをさらに含み得る。
そのペーシングデータの検出は、例えば、ＥＣＧ波形の連続するＲピーク間の時間を検出するステップを含み得る。
心拍数を表すペーシングデータは個体の身体から測定される第３信号に基づいて作成されることもあり得る。
【００１３】
タイミングデータとペーシングデータとの間の上述の相関は心臓の機能不全と異常の指標として使用され得る。
例えば、実証的なデータに照らすと、ペーシングデータの各値が心拍数の瞬時値を、例えば、分当たりの心拍数で表すとき、正の相関は心房細動の可能性の上昇を意味する。
同様に、ペーシングデータの各値が１回の心拍周期の時間の長さを、例えば、秒で、すなわち１／心拍数で表すとき、負の相関は心房細動の可能性の上昇を意味する。
【００１４】
タイミングデータとペーシングデータとの間の相関の程度は、必ずというわけではないが、例えば、次の式：
【数１】
2018086278.tif 000003
に従って計算され得る相関係数を用いて表現され得、
式中、
Ｃ（ｊ）はその相関係数であり、
Ｅは期待値演算子であり、すなわちＥ｛変数｝は変数の期待値であり、
ＴＤはタイミングデータであり、
μTはタイミングデータの平均であり、
ＰＤはペーシングデータであり、
μPはペーシングデータの平均であり、
ｊは心拍周期内のタイミングデータに関してペーシングデータの時間差を表す整数である。
実証的な結果に照らすと、ペーシングデータＰＤがタイミングデータＴＤに関して１回の心拍周期の時間差、すなわちｊ＝１を有することが有利である。
この場合、タイミングデータＴＤが所与の心拍周期に関連するとき、対応するペーシングデータＰＤは先の心拍周期に関連する。
相関係数は、フォルムρT,Pで表され得、それは、常に−１から＋１までの範囲であり、
【数２】
2018086278.tif 000004
（σTとσPがそれぞれタイミングデータとペーシングデータの標準偏差である）である。
【００１５】
タイミングデータとペーシングデータとの間の可能性がある相関を表現する多数の方法が存在し、本発明はその相関を表現するどの特定の方法にも限定されないことに留意すべきである。
さらに、その相関は必ずしも数学量ではないが、従属を含む広範な統計学的関係のいずれかに関連しており、その相関はその広い意味で因果関係を含むこと、又は必要とすることが無いことに留意すべきである。
【００１６】
本発明の別の好ましい実施形態による方法では、心臓の機能不全と異常の指標の測定は実質的に一定の心拍数におけるタイミングデータの変化の測定を含む。
例えば、実証的なデータに照らすと、タイミングデータの標準偏差は心房細動の間では平均値の約１０％であり得、正常状態の場合では平均値の約１〜２％であり得る。
【００１７】
本発明の好ましい実施形態による方法では、心臓血管の動きを表す信号をローパスフィルタにかけ、そのローパスフィルタにかけられた信号から上述のＡＯピークを検出するステップ、及び／又は、心臓血管の動きを表す前記信号をバンドパスフィルタにかけ、大動脈弁の閉鎖により生じるＡＣピークをそのバンドパスフィルタにかけられた信号から検出するステップ、を含む。
ローパスフィルタリングの上限の周波数は、必ずというわけではないが、例えば、３０Ｈｚであり得、バンドパスフィルタリングの通過帯域は、必ずというわけではないが、例えば、４０Ｈｚから１００Ｈｚまでであり得る。ローパスフィルタリング及び／又はパスバンドフィルタリングはＡＯピーク及び／又はＡＣピークの検出を容易にする。
特に心房細動の間にバンドパスフィルタリングを用いると、バンドパスフィルタリングが無いときよりもＡＣピークの発見が容易である。
検出されたＡＯピーク及び／又はＡＣピークは、例えば上述の第２波形を心臓血管の動きを表す信号から抽出するときに利用され得る。
検出されたＡＯピーク及び／又はＡＣピークは他の多くの目的のためにも、例えば振幅の変化、時間変化、心拍数、収縮間隔、及び／又は拡張間隔を検出するために使用され得る。
【００１８】
本発明によって、心臓の機能不全と異常を表す情報を測定するための新しい器械が提供される。
本発明による器械は、
−心臓活動に関連する電磁現象を表す第１信号と心臓血管の動きを表す第２信号を受信するための信号インターフェースと、
−その信号インターフェースに接続されている処理装置であって、
ａ）前記第１信号からある心拍数を繰り返す第１波形を抽出し、且つ、前記第２信号から前記心拍数を繰り返す第２波形を抽出し、
ｂ）各タイミング値がある１つの心拍周期に属する前記第１波形の基準点から同一の心拍周期に属する前記第２波形の基準点までの時間を表すタイミングデータを作成し、且つ、
ｃ）前記タイミングデータと前記の心拍数を表すペーシングデータとの間の相関を決定する、
ように構成されている前記処理装置と、
を備えており、
その相関が心臓の機能不全と異常を表す。
【００１９】
その器械は心臓活動に関連する電磁現象を表す第１信号を測定するための第１センサ素子及び／又は心臓血管の動きを表す第２信号を測定するための第２センサ素子をさらに備え得る。
その第１センサ素子は個体の身体へ装着に適切な電極を備え得る。
その第２センサ素子は、例えば、加速度センサ、圧電センサ、傾斜センサ、圧力センサ、又は力、加速度、変位、若しくは心臓血管の動きに関連し、心臓血管の動きを表す他のあらゆる量の測定に適切な他のあらゆる素子を備え得る。
前記信号インターフェースは適切なセンサ素子を備える外部機器から前記第１信号及び／又は前記第２信号を受信することができるということもあり得、すなわち前記器械は前記第１信号及び／又は前記第２信号を測定するための手段を必ず備えるわけではないことが強調される。
【００２０】
本発明の好ましい実施形態による器械は、心臓血管の動きを表す信号をローパスフィルタにかけるためのローパスフィルタとそのローパスフィルタにかけられた信号からＡＯピークを検出するための手段、例えばプロセッサ、及び／又は心臓血管の動きを表す前記信号をバンドパスフィルタにかけるためのバンドパスフィルタとそのバンドパスフィルタにかけられた信号からＡＣピークを検出するための手段を備える。
【００２１】
本発明により、心臓の機能不全と異常を表す情報を測定するための新しいコンピュータプログラムが提供される。
そのコンピュータプログラムは、
−心臓活動に関連する電磁現象を表す第１信号からある心拍数を繰り返す第１波形を抽出し、
−心臓血管の動きを表す第２信号から前記心拍数を繰り返す第２波形を抽出し、
−各タイミング値がある１つの心拍周期に属する前記第１波形の基準点から同一の心拍周期に属する前記第２波形の基準点までの時間を表すタイミングデータを作成し、且つ、
−前記タイミングデータと前記の心拍数を表すペーシングデータとの間の相関を決定する、
プログラミング可能なプロセッサを制御するためのコンピュータが実行可能な命令を備え、
その相関が心臓の機能不全と異常を表す。
【００２２】
本発明の好ましい実施形態によるコンピュータプログラムは、
心臓血管の動きを表す信号をローパスフィルタにかけ、そのローパスフィルタにかけられた信号からＡＯピークを検出する、プログラミング可能なプロセッサを制御するためのコンピュータが実行可能な命令、及び／又は、
心臓血管の動きを表す信号をバンドパスフィルタにかけ、そのバンドパスフィルタにかけられた信号からＡＣピークを検出する、プログラミング可能なプロセッサを制御するためのコンピュータが実行可能な命令を備える。
【００２３】
本発明により、新しいコンピュータプログラム製品が提供される。そのコンピュータプログラム製品は、本発明によるコンピュータプログラムが符号化されている不揮発性のコンピュータ可読媒体、例えば、コンパクトディスク「ＣＤ」を備える。
【００２４】
本発明の多数の好ましい実施形態は添付されている従属請求項に記載される。
構成と操作方法に関する本発明の様々な好ましい実施形態はそれらの追加的な目的及び利点と共に、添付されている図面と関連させて読まれると、具体的な好ましい実施形態の次の記載から最も良く理解される。
「備える（ｔｏ  ｃｏｍｐｒｉｓｅ）」及び「含む（ｔｏ  ｉｎｃｌｕｄｅ）」という動詞は、本書類において、列挙されていない特徴の存在を排除することも必要とすることもない開放的な制限語として使用される。従属請求項において列挙される特徴は、別途明らかに示されていない限り、自由に相互に組合せ可能である。
【図面の簡単な説明】
【００２５】
本発明の好ましい実施形態とそれらの利点が、添付されている図面の参照により以下にさらに詳しく説明される。
【図１ａ】ＥＣＧ波形の例を示す図である。
【図１ｂ】通常はｙ方向と呼ばれる「頭足」方向で加速度センサにより測定された心臓血管の動きを表す例となる信号の波形を示す図である。
【図２ａ】心臓の機能不全と異常を表す情報を測定するための本発明の好ましい実施形態による方法のフローチャートを示す図である。
【図２ｂ】前記の心臓血管の動きを表す信号からＡＯデータ及び／又はＡＣデータを抽出するための本発明の好ましい実施形態による方法のフローチャートを示す図である。
【図３ａ−３ｂ】図３ａは、例となるＥＣＧ波形を示す図である。図３ｂは、通常はｚ方向と呼ばれる「通胸」方向で加速度センサにより測定された、心臓血管の動きを表す例となる信号の波形を示す図である。
【図４】心臓の機能不全と異常を表す情報を測定するための本発明の好ましい実施形態による器械の概略図を示す図である。
【図５】例となる心房細動の場合の様々な心拍数におけるＥＣＧ波形のＲピークと心臓血管の動きを表す波形の大動脈弁開放「ＡＯ」ピークの間の時間差を示す図である。
【発明を実施するための形態】
【００２６】
図１ａ、図１ｂは本発明の背景技術を述べたときに既に説明されている。
図２ａは心臓の機能不全と異常、例えば心房細動を表す情報を測定するための本発明の好ましい実施形態による方法のフローチャートを示す。
その方法は心臓活動に関連する電磁現象を表す第１信号からある心拍数を繰り返す第１波形をステップ２０１において抽出することを含む。
その方法は心臓血管の動きを表す第２信号からその心拍数を繰り返す第２波形をステップ２０２において抽出することを含む。
その第１信号と第２信号は同一の個体から測定される／測定されたものである。
その方法はタイミングデータの各タイミング値がある１つの心拍周期に属する前記第１波形の基準点から同一の心拍周期に属する前記第２波形の基準点までの時間を表すように、そのタイミングデータＴＤをステップ２０３において作成することを含む。
その方法はそのタイミングデータＴＤに少なくとも部分的に基づいて心臓の機能不全と異常の指標をステップ２０４において測定することを含む。
【００２７】
図３ａは心臓活動に関連する電磁現象を表す例となる信号の波形を示し、図３ｂは心臓血管の動きを表す例となる信号の波形を示す。
図３ａ、図３ｂに示されている波形は互いに対応し、同一の個体から同時に測定されたものである。
図３ａに示されている波形は個体の胸に装着された電極を用いて測定されたＥＣＧ波形であり、図３ｂに示されている波形は通常はｚ方向と呼ばれる「通胸」方向で加速度センサにより測定されたものである。
図４に示されている座標系４５０がｚ方向を示す。
【００２８】
図３ａ、図３ｂは共にタイミングデータＴＤを定義するための例となる方法を示す。
この例となる場合では、ＥＣＧ波形に現れ、且つ、心室性筋肉組織の脱分極により生じるＲピークは心拍数を繰り返す第１波形を表し、心臓血管の動きを表す波形の大動脈弁開放「ＡＯ」ピークは心拍数を繰り返す第２波形を表す。Ｒピークの頂点がその第１波形の基準点であり、ＡＯピークの頂点がその第２波形の基準点である。
【００２９】
本発明の別の好ましい実施形態による方法では、Ｒピークは前記第１波形を表し、頭足方向で測定された波形に現れ、且つ、大動脈弁を通る左心室から大動脈への血流によって生じるＪピークは前記第２波形を表す。
Ｒピークの頂点はその第１波形の基準点として使用され得、Ｊピークの頂点はその第２波形の基準点として使用され得る。
【００３０】
本発明の１つの好ましい実施形態による方法では、Ｒピークは前記第１波形を表し、図３ｂに示されている波形上のＭＣピークとＡＯピークの間のバレー３１０は前記第２波形を表す。
ＭＣピークは僧帽弁の閉鎖によって生じる。
Ｒピークの頂点は前記第１波形の前記基準点として使用され得、バレー３１０の最低点は前記第２波形の前記基準点として使用され得る。
【００３１】
図３ａ、図３ｂは３回の連続する心拍周期、すなわち、「ｉ」が整数であり得る心拍周期「ｉ−１」、心拍周期「ｉ」、及び心拍周期「ｉ＋１」を示す。
例えば、心拍周期「ｉ」に関連するタイミング値ＴＤ（ｉ）は、図３ａ、図３ｂに示されるように、心拍周期「ｉ」のＲピークの頂点の瞬間から心拍周期「ｉ」のＪピークの頂点の瞬間までの時間である。
【００３２】
図２ａに示されているステップ２０４は前記タイミングデータＴＤと前記の心拍数を表すペーシングデータＰＤとの間の相関を決定することを含む。
実証的な結果に照らすと、その相関は心臓の機能不全と異常の前記指標として使用され得る。
【００３３】
その方法は、前記の心拍数を表すペーシングデータＰＤを心臓活動に関連する電磁現象を表す波形及び／又は心臓血管の動きを表す波形から検出することをさらに含み得る。
そのペーシングデータの検出は、例えば、ＥＣＧ波形上の連続するＲピーク間の時間の検出を含み得る。
そのペーシングデータの検出は、例えば、心臓血管の動きを表す波形上の連続するＡＯピーク又はＪピークの間の時間の検出を含むこともあり得る。
より信頼できるペーシングデータを得るために１つより多くの波形を活用することができる。
さらに、第３信号が個体の身体から測定され、この第３信号が単独で、又は心臓活動に関連する電磁現象を表す波形及び／又は心臓血管の動きを表す波形と共に前記の心拍数を表すペーシングデータの測定のために用いられることもあり得る。
図３ａでは、連続するＲピーク間の時間の長さはＨＢ（ｉ−１）、ＨＢ（ｉ）、及びＨＢ（ｉ＋１）で示されている。
心拍周期「ｉ−１」、「ｉ」、及び「ｉ＋１」のペーシングデータはそれぞれ例えばＰＤ（ｉ−１）＝１／ＨＢ（ｉ−１）、ＰＤ（ｉ）＝１／ＨＢ（ｉ）、及びＰＤ（ｉ＋１）＝１／ＨＢ（ｉ＋１）として定義され得る。
【００３４】
タイミングデータＴＤとペーシングデータＰＤとの間の相関は、例えば、次の式：
【数３】
2018086278.tif 000005
に従って計算され得る数学的相関係数を用いて表現され得、
式中、
Ｎは検討中の心拍周期の回数であり、
ｊは心拍周期内のタイミングデータＴＤに関してペーシングデータＰＤの時間差を表す整数であり、
【数４】
2018086278.tif 000006
である。
【００３５】
上記の相関係数Ｃ（ｊ）は心臓の機能不全と異常の前記指標として使用され得る。実証的なデータに照らして、相関係数Ｃ（ｊ＝１）の正の値は心房細動の可能性の上昇を意味し、その相関の負の値は正常状態を表す。
【００３６】
図２ｂは心臓血管の動きを表す信号からＡＯデータ及び／又はＡＣデータを抽出するための本発明の好ましい実施形態による方法のフローチャートを示す。
前記の心臓血管の動きを表す信号が、通常はｚ方向と呼ばれる「通胸」方向で加速度センサにより測定されることが好都合である。
その方法は次のステップを含む：
−ステップ２１１：心臓血管の動きを表す信号をローパスフィルタにかけること、及び／又は心臓血管の動きを表す信号をバンドパスフィルタにかけること、及び
−ステップ２１２：そのローパスフィルタにかけられた信号からＡＯピークを検出すること、及び／又はそのバンドパスフィルタにかけられた信号からＡＣピークを検出すること。
【００３７】
ローパスフィルタリングの上限の周波数は、必ずというわけではないが、例えば、３０Ｈｚであり得、バンドパスフィルタリングの通過帯域は、必ずというわけではないが、例えば、４０Ｈｚから最大で１００Ｈｚまでであり得る。
ローパスフィルタリング及び／又はパスバンドフィルタリングはＡＯピーク及び／又はＡＣピークの検出を容易にする。
検出されたＡＯピーク及び／又はＡＣピークは例えば上述の第２波形を心臓血管の動きを表す信号から抽出するときに利用され得る。
検出されたＡＯピーク及び／又はＡＣピークは他の多くの目的のためにも、例えば振幅の変化、時間変化、心拍数、収縮間隔、及び／又は拡張間隔を検出するために使用され得る。
【００３８】
本発明の好ましい実施形態による方法はＡＣ〜ＡＯ間隔の時間の長さを検出すること、及び検出されたＡＣ〜ＡＯ間隔の時間の長さの変化の強度を表す時間変化量を計算することを含み、前記のＡＣ〜ＡＯ間隔のそれぞれがＡＣピークのうちの１つからＡＯピークのうちの次の１つまでの時間間隔であり、その時間変化量が心臓の機能不全と異常を表す。
この時間変化量は可能性がある心臓の機能不全と異常の検出の信頼性を上昇させるためにタイミングデータとペーシングデータとの間の上述の相関と共に使用され得る。
しかしながら、この時間変化量は単独でも使用され得る。
【００３９】
本発明の好ましい実施形態による方法はＡＣ〜ＡＯ間隔の時間の長さとＡＯ〜ＡＯ間隔の時間の長さを検出すること、及び同一の心拍周期内のＡＣ〜ＡＯ間隔の時間の長さとＡＯ〜ＡＯ間隔の時間の長さの間の比率を表す第１比率量を計算することを含む。
前記のＡＣ〜ＡＯ間隔のそれぞれはＡＣピークのうちの１つからＡＯピークのうちの次の１つまでの時間間隔であり、前記のＡＯ〜ＡＯ間隔のそれぞれはＡＯピークのうちの１つからＡＯピークのうちの次の１つまでの時間間隔であり、その第１比率量が心臓の機能不全と異常を表す。
この第１比率量は可能性がある心臓の機能不全と異常の検出の信頼性を上昇させるためにタイミングデータとペーシングデータとの間の上述の相関と共に使用され得る。
しかしながら、この第１比率量は単独でも使用され得る。
【００４０】
本発明の好ましい実施形態による方法はＡＣ〜Ｒ間隔の時間の長さとＲ〜Ｒ間隔の時間の長さを検出すること、及び同一の心拍周期内の前記ＡＣ〜Ｒ間隔の時間の長さと前記Ｒ〜Ｒ間隔の時間の長さの間の比率を表す第２比率量を計算すことを含む。
前記のＡＣ〜Ｒ間隔のそれぞれはＡＣピークのうちの１つからＲピークのうちの次の１つまでの時間間隔であり、前記のＲ〜Ｒ間隔のそれぞれはＲピークのうちの１つからＲピークのうちの次の１つまでの時間間隔であり、その第２比率量が心臓の機能不全と異常を表す。
この第２比率量は可能性がある心臓の機能不全と異常の検出の信頼性を上昇させるためにタイミングデータとペーシングデータとの間の上述の相関と共に使用され得る。しかしながら、この第２比率量は単独でも使用され得る。
【００４１】
図５は例となる心房細動の場合の様々な心拍数において得られたタイミングデータの実証的な値を示す。
図５中の各黒点は
ａ）ある特定の心拍周期内のＥＣＧ波形のＲピークと心臓血管の動きを表す波形の大動脈弁開放「ＡＯ」ピークとの間の時間差と
ｂ）先の心拍周期における心拍数と、を表す。
図５に示されている破線５０１の正の傾きはタイミングデータとペーシングデータとの間の正の相関を示す。
図５から理解することができるように、そのタイミングデータの傾きは心拍数が分当たり約８５拍を越えると心拍数とともに上昇する。
正常状態の場合であれば、その傾きは実質的に一定であるか、又は減少する。
【００４２】
本発明の好ましい実施形態による方法は、相関係数Ｃ（ｊ＝１）が閾値よりも大きい状況に応答して心房細動を表す信号を発信することを含む。
その閾値に適切な値は、一群の患者及び／又は他の人間から集められた実証的データに基づいて決定され得る。
その閾値は必ずしも一定というわけではなく、その閾値は検討中の個体に応じて、時間に応じて、及び／又は他の幾つかの因子に応じて変化し得る。
各閾値が心房細動の特定の可能性、又は他の幾つかの心臓の機能不全及び／又は異常の特定の可能性を表す一連の閾値を設定することも可能である。
幾つかの場合では、その閾値は０であり得る。
【００４３】
別の例では、心拍周期「ｉ−１」、「ｉ」、及び「ｉ＋１」のペーシングデータはそれぞれＰＤ（ｉ−１）＝ＨＢ（ｉ−１）、ＰＤ（ｉ）＝ＨＢ（ｉ）、及びＰＤ（ｉ＋１）＝ＨＢ（ｉ＋１）と定義され得る。
この場合、相関係数Ｃ（ｊ＝１）の負の値は心房細動の可能性の上昇を意味し、その相関の正の値は正常状態を表す。
本発明の好ましい実施形態による方法は、相関係数Ｃ（ｊ＝１）が閾値未満であり、それが、幾つかの場合では、０でありうる状況に応答して心房細動を表す信号を発信することを含む。
【００４４】
本発明の好ましい実施形態による方法では、図２ａに示されているステップ２０４は心拍数が実質的に一定であるときのタイミングデータの変化を測定するステップを含む。
実証的なデータに照らして、その変化の程度は心臓の機能不全と異常を表す。
図５では、心拍数が８０拍／分であるときのその変化の範囲が垂直の破線部分５０２によって示されている。
タイミングデータが図３ａ及び３ｂに示されるようにＲピークとＡＯピークの間の時間間隔を示すとき、そのタイミングデータは前駆出期「ＰＥＰ」を示す。
心房細動の期間において連続する心拍周期の間にＰＥＰの確率的な変化があることが気付かれている。
【００４５】
上述の変化の程度は、例えば、次の式：
【数５】
2018086278.tif 000007
に従って計算され得る数学的変化量を用いて表現され得、
式中、
Ｖはその変化量であり、
Ｍは検討中の心拍数における検討中のタイミングデータ値の数であり、
【数６】
2018086278.tif 000008
である。
実証的なデータに照らして、変化量Ｖは心房細動の間に約１０％であり得、正常状態の場合に約１〜２％であり得る。
【００４６】
本発明の好ましい実施形態による方法は、変化量Ｖが閾値よりも大きい状況に応答して心房細動を表す信号を発信するステップを含む。
その閾値に適切な値は、一群の患者及び／又は他の人間から集められた実証的データに基づいて決定され得る。
その閾値は必ずしも一定というわけではなく、その閾値は検討中の個体に応じて、時間に応じて、及び／又は他の幾つかの因子に応じて変化し得る。
各閾値が心房細動の特定の可能性、又は他の幾つかの心臓の機能不全及び／又は異常の特定の可能性を表す一連の閾値を設定することも可能である。
【００４７】
本発明の好ましい実施形態による方法は、１回よりも多くの心拍数におけるタイミングデータの変化を検出するステップ、及び心臓の機能不全と異常の前記指標を測定するために検出結果を使用するステップを含む。
例えば、式（３）及び（４）は数回の心拍数における変化量を得るために使用され得、最終的な変化量は、検討中の心拍数に関連する変化量から数学的論理的演算、例えば、算術平均により計算され得る。
【００４８】
各ピーク値、例えば１つのＡＯピークの高さは、極大を捜し求めることによる１つの点であると見なされ得る。
あるいは、そのピーク値は、最初に検討中のピークを含む時間ウィンドウから多数の試料を獲得し、次に、ノイズ効果を軽減するためにそれらの試料の数学的関数、例えば算術平均としてそのピーク値を計算するようにして取得され得る。
その時間ウィンドウは、例えば、１００ミリ秒であり得、その時間ウィンドウ内の試料数は、例えば、１０以上であり得る。
その時間ウィンドウに基づく方法はデジタルフィルタリングの一例である。
一般に、心臓血管の動きを表す信号及び心臓活動に関連する電磁現象を表す信号におけるノイズ効果を軽減するために用いられ得る多数のデジタル信号処理方法及びアナログ信号処理方法が存在する。
【００４９】
本発明の好ましい実施形態による方法は、個体の身体よりセンサ素子を用いて心臓活動に関連する電磁現象を表す第１信号と心臓血管の動きを表す第２信号を測定するステップを所望により含む。
本発明の別の好ましい実施形態による方法は、メモリからこれらの信号を読むステップを含み、その場合、それらの信号は前に測定されており、そのメモリに記録されている。
本発明の好ましい実施形態による方法は、外部データ転送システムからそれらの信号を受信するステップを含む。
したがって、測定は、本発明の実施形態による方法の必須であり、必要なステップというわけではない。
【００５０】
図４は心臓の機能不全と異常を表す情報を測定するための本発明の好ましい実施形態による器械の概略図を示す。
その器械は心臓活動に関連する電磁現象を表す第１信号と心臓血管の動きを表す第２信号を受信するための信号インターフェース４０１を備える。
その器械は前記信号インターフェースに接続されている処理装置４０２を備える。
その処理装置は、
−前記第１信号からある心拍数を繰り返す第１波形を抽出し、且つ、前記第２信号から前記心拍数を繰り返す第２波形を抽出し、
−各タイミング値がある１つの心拍周期に属する前記第１波形の基準点から同一の心拍周期に属する前記第２波形の基準点までの時間を表すタイミングデータを作成し、且つ、
−前記タイミングデータと前記の心拍数を表すペーシングデータとの間の相関を決定する、
ように構成されており、
その相関が心臓の機能不全と異常を表す。
【００５１】
その第１信号は、例えば、心電計「ＥＣＧ」波形又は誘導的に測定された波形を表し得る。
【００５２】
本発明の好ましい実施形態による器械では、処理装置４０２は、ＲピークをＥＣＧ波形から抽出し、且つ、ＡＯピークを心臓血管の動きを表す波形から抽出する、ように構成されている。
この場合、Ｒピークは前記第１信号上で前記心拍数を繰り返す第１波形を表し、ＡＯピークは前記第２信号上で前記心拍数を繰り返す第２波形を表す。
Ｒピークの頂点が前記第１波形の基準点であり得、ＡＯピークの頂点が前記第２波形の基準点であり得る。
【００５３】
本発明の別の好ましい実施形態による器械では、処理装置４０２は、ＲピークをＥＣＧ波形から抽出し、Ｊピークを心臓血管の動きを表す波形から抽出する、ように構成されている。
この場合、Ｒピークは前記第１信号上で前記心拍数を繰り返す第１波形を表し、Ｊピークは前記第２信号上で前記心拍数を繰り返す第２波形を表す。
Ｒピークの頂点が前記第１波形の基準点であり得、Ｊピークの頂点が前記第２波形の基準点であり得る。
【００５４】
本発明の好ましい実施形態による器械では、処理装置４０２は、ＲピークをＥＣＧ波形から抽出し、ＭＣピークとＡＯピークの間の谷を心臓血管の動きを表す波形から抽出する、ように構成されている。
この場合、Ｒピークは前記第１信号上で前記心拍数を繰り返す第１波形を表し、ＭＣピークとＡＯピークの間の前記の谷は前記第２信号上で前記心拍数を繰り返す第２波形を表す。
Ｒピークの頂点が前記第１波形の基準点であり得、前記の谷の最低点が前記第２波形の基準点であり得る。
【００５５】
本発明の好ましい実施形態による器械では、処理装置４０２は、前記タイミングデータと前記ペーシングデータとの間の前記相関の程度を表す相関係数：
【数７】
2018086278.tif 000009
を計算するように構成されており、
式中、
Ｃ（ｊ）はその相関係数であり、
Ｅは期待値演算子であり、
ＴＤはタイミングデータであり、
μTはタイミングデータの平均であり、
ＰＤはペーシングデータであり、
μPはペーシングデータの平均であり、
ｊは心拍周期内のタイミングデータに関してペーシングデータの時間差を表す整数である。
【００５６】
さらに、処理装置４０２は、計算された相関係数Ｃ（ｊ）を１つ以上の閾値と比較し、且つ、その比較の結果が心臓の機能不全と異常の存在を示す状況に応答して心臓の機能不全と異常を表す信号を発信する、ように構成され得る。
その信号は、例えば、ディスプレー画面４０６上に示されるメッセージであり得る。
それらの１つ以上の閾値は、前記器械にそのユーザーインターフェースを介して提供され得る調節可能パラメーターであることが好ましい。
【００５７】
本発明の好ましい実施形態による器械では、処理装置４０２は前記第１信号上の連続するＲピーク間の時間を検出する、ように構成されている。
検出された時間は前記心拍数に反比例するものとしてその心拍数を表す。
したがって、それらの検出された時間をペーシングデータとして使用することができ、又はそれらの検出された時間を活用してペーシングデータを作成することができる。
【００５８】
本発明の好ましい実施形態による器械では、処理装置４０２は、ペーシングデータの各値が前記心拍数の瞬時値を表すようにペーシングデータを作成し、且つ、相関係数が所定の閾値、例えば０よりも大きい状況に応答して心房細動を表す信号を発信する、ように構成されている。
【００５９】
本発明の好ましい実施形態による器械では、処理装置４０２は、ペーシングデータの各値が１回の心拍周期の時間の長さを表すようにペーシングデータを作成し、且つ、相関係数が所定の閾値、例えば０未満である状況に応答して心房細動を表す信号を発信する、ように構成されている。
【００６０】
本発明の好ましい実施形態による器械では、処理装置４０２は実質的に一定の心拍数に対応するタイミングデータの変化を測定する、ように構成されている。
その処理装置は、
その変化の程度を表す変化量：
【数８】
2018086278.tif 000010
を測定するように構成され得、
式中、
Ｖはその変化量であり、
Ｅは期待値演算子であり、
ＴＤはタイミングデータであり、
μTは検討中の前記心拍数におけるタイミングデータの平均である。
【００６１】
さらに、処理装置４０２は、計算された変化量Ｖを１つ以上の閾値と比較し、且つ、その比較の結果が心臓の機能不全と異常の存在を示す状況に応答して心臓の機能不全と異常を表す信号を発信する、ように構成され得る。
その信号は、例えば、ディスプレー画面４０６上に示されるメッセージであり得る。
それらの１つ以上の閾値は、前記器械にそのユーザーインターフェースを介して提供され得る調節可能パラメーターであることが好ましい。
【００６２】
本発明の好ましい実施形態による器械では、処理装置４０２は、１回よりも多くの心拍数におけるタイミングデータの変化を検出し、且つ、それらの１回より多くの心拍数に関連する変化に基づいて心臓の機能不全と異常の指標を測定するように構成されている。
例えば、式（３）及び（４）は数回の心拍数における変化量を得るために使用され得る。
最終的な変化量は、検討中の心拍数に関連する変化量から数学的論理的演算、例えば、算術平均により計算され得る。
【００６３】
本発明の好ましい実施形態による器械は個体の身体４０９から前記第１信号を測定するための第１センサ素子４０３とその個体の身体から前記第２信号を測定するための第２センサ素子４０４をさらに備える。
その第１センサ素子と第２センサ素子は、それぞれが、例えば、無線リンク又は有線リンクであり得る１つ以上のデータ転送リンクを介して前記信号インターフェースに接続されている。
センサ素子４０３及び４０４から信号インターフェース４０１へのデータ転送は直接的に、又は例えば電気通信ネットワークなどのデータ転送ネットワーク４０５を介して起こり得る。
図４に示されている例となる事例では、センサ素子４０３とセンサ素子４０４の両方が無線送信機４０８に接続されている。
処理装置４０６を備える器械が前記のセンサ素子と一体化していることもあり得る。
この場合、前記信号インターフェースは実際にはセンサ素子４０３及び４０４から処理装置４０２までの単なる配線である。
【００６４】
この例となる場合では、第１センサ素子４０３は個体の身体へ装着に適切な電極を備える。その第１センサ素子は、例えば、増幅器、信号フィルタ、及び／又はアナログデジタル「ＡＤ」変換器をさらに備え得る。
第２センサ素子４０４は、例えば、加速度センサ、圧電センサ、傾斜センサ、圧力センサ、又は力、加速度、変位、若しくは心臓血管の動きに関連し、心臓血管の動きを表す他のあらゆる物理量の測定に適切な他のあらゆる素子を備え得る。
その第２センサ素子は、例えば、増幅器、信号フィルタ、及び／又はアナログデジタル「ＡＤ」変換器をさらに備え得る。
加速度センサが、例えば図４に示されている座標系４５０の３つの互いに直交するｘ方向、ｙ方向、及びｚ方向で独立的に動きを測定することができる３軸加速度センサであることが有利である。
この場合、心臓血管の動きを表す第２信号は３つの成分を含み、その第２信号、例えば、そのユークリッドノルム、すなわち、心臓血管の動きを表す３成分ベクトルの絶対値を作成することにより前処理され得る。
【００６５】
本発明の好ましい実施形態による器械は、固定された開始時点と固定された終止時点を有する時間ウィンドウ内、又は固定された時間の長さを有し、且つ、時間経過と共に移動するスライディング時間ウィンドウ内で第１信号と第２信号を記録する、ように構成されている。
その器械は、前記の信号を記録するための内部メモリ４０７を備えることができ、及び／又は、その器械は外部メモリへの接続のためのデータポートを備えることができる。
【００６６】
本発明の好ましい実施形態による器械は心臓活動に関連する電磁現象を表す第１信号及び／又は心臓血管の動きを表す第２信号を前処理するための手段を備える。
その前処理は、例えば、呼吸、個体の非心臓血管の動き、外部的理由により引き起こされた振動等により生じた前記第２信号上のノイズの取り消しを含み得る。
その前処理のための手段を、例えば、処理装置４０２に装備することができ、又はその前処理のための１つ以上の別の処理装置が存在し得る。
【００６７】
本発明の好ましい実施形態による器械では、処理装置４０２は、心臓血管の動きを表す第２信号をローパスフィルタにかけ、上述のＡＯピークをそのローパスフィルタにかけられた信号から検出する、ように構成されている。
図４に示されている機能ブロック４２０はそのローパスフィルタリングを表し、機能ブロック４２２はＡＯピークの検出を表す。
【００６８】
本発明の好ましい実施形態による器械では、処理装置４０２は、心臓血管の動きを表す第２信号をバンドパスフィルタにかけ、大動脈弁の閉鎖によって生じるＡＣピークをそのバンドパスフィルタにかけられた信号から検出する、ように構成されている。
図４に示されている機能ブロック４２１はそのバンドパスフィルタリングを表し、機能ブロック４２２はＡＣピークの検出を表す。
図４に示されている例となる事例では、ローパスフィルタリングとバンドフィルタリングの両方が存在し、機能ブロック４２２がＡＯピークとＡＣピークの両方の検出を表す。
ローパスフィルタリングの上限の周波数は、必ずというわけではないが、例えば、３０Ｈｚであり得、バンドパスフィルタリングの通過帯域は、必ずというわけではないが、例えば、４０Ｈｚから最大で１００Ｈｚまでであり得る。
ローパスフィルタリング及び／又はパスバンドフィルタリングはＡＯピーク及び／又はＡＣピークの検出を容易にする。
検出されたＡＯピーク及び／又はＡＣピークは利用され得る例えば上述の第２波形を心臓血管の動きを表す信号から抽出するときに．検出されたＡＯピーク及び／又はＡＣピークは他の多くの目的のためにも、例えば振幅の変化、時間変化、心拍数、収縮間隔、及び／又は拡張間隔を検出するために使用され得る。
【００６９】
本発明の好ましい実施形態による器械では、処理装置５０２はＡＣ〜ＡＯ間隔の時間の長さを検出し、且つ、検出されたＡＣ〜ＡＯ間隔の時間の長さの変化の強度を表す時間変化量を計算する、ように構成されている。
前記のＡＣ〜ＡＯ間隔のそれぞれはＡＣピークのうちの１つからＡＯピークのうちの次の１つまでの時間間隔であり、前記の時間変化量は心臓の機能不全と異常を表す。
この時間変化量は可能性がある心臓の機能不全と異常の検出の信頼性を上昇させるためにタイミングデータとペーシングデータとの間の上述の相関と共に使用され得る。
しかしながら、この時間変化量は単独でも使用され得る。
【００７０】
本発明の好ましい実施形態による器械では、処理装置５０２はＡＣ〜ＡＯ間隔の時間の長さとＡＯ〜ＡＯ間隔の時間の長さを検出し、且つ、同一の心拍周期内のＡＣ〜ＡＯ間隔の時間の長さとＡＯ〜ＡＯ間隔の時間の長さの間の比率を表す第１比率量を計算するように構成されている。
前記のＡＣ〜ＡＯ間隔のそれぞれはＡＣピークのうちの１つからＡＯピークのうちの次の１つまでの時間間隔であり、前記のＡＯ〜ＡＯ間隔のそれぞれはＡＯピークのうちの１つからＡＯピークのうちの次の１つまでの時間間隔であり、その第１比率量が心臓の機能不全と異常を表す。
この第１比率量は可能性がある心臓の機能不全と異常の検出の信頼性を上昇させるためにタイミングデータとペーシングデータとの間の上述の相関と共に使用され得る。しかしながら、この第１比率量は単独でも使用され得る。
【００７１】
本発明の好ましい実施形態による器械では、処理装置５０２は、ＡＣ〜Ｒ間隔の時間の長さとＲ〜Ｒ間隔の時間の長さを検出し、且つ、同一の心拍周期内の前記のＡＣ〜Ｒ間隔の時間の長さと前記のＲ〜Ｒ間隔の時間の長さの間の比率を表す第２比率量を計算するように構成されている。
前記のＡＣ〜Ｒ間隔のそれぞれはＡＣピークのうちの１つからＲピークのうちの次の１つまでの時間間隔であり、前記のＲ〜Ｒ間隔のそれぞれはＲピークのうちの１つからＲピークのうちの次の１つまでの時間間隔であり、その第２比率量が心臓の機能不全と異常を表す。
この第２比率量は可能性がある心臓の機能不全と異常の検出の信頼性を上昇させるためにタイミングデータとペーシングデータとの間の上述の相関と共に使用され得る。
しかしながら、この第２比率量は単独でも使用され得る。
【００７２】
処理装置４０２に、例えば、１つ以上のプロセッサ回路を装備することができ、それらの回路のそれぞれが適切なソフトウェア、例えば特定用途向け集積回路「ＡＳＩＣ」などの専用ハードウェアプロセッサ、又は、例えばフィールドプログラマブルゲートアレイ「ＦＰＧＡ」などの構成可能ハードウェアプロセッサと共に提供されるプログラミング可能なプロセッサ回路であり得る。
【００７３】
本発明の好ましい実施形態によるコンピュータプログラムは心臓の機能不全と異常を表す情報を測定するためのソフトウェアモジュールを備える。
それらのソフトウェアモジュールは、
−心臓活動に関連する電磁現象を表す第１信号からある心拍数を繰り返す第１波形を抽出し、
−心臓血管の動きを表す第２信号から前記心拍数を繰り返す第２波形を抽出し、
−各タイミング値がある１つの心拍周期に属する前記第１波形の基準点から同一の心拍周期に属する前記第２波形の基準点までの時間を表すタイミングデータを作成し、且つ、
−前記タイミングデータと前記の心拍数を表すペーシングデータとの間の相関を決定する、
プログラミング可能なプロセッサを制御するためのコンピュータが実行可能な命令を備え、
その相関が心臓の機能不全と異常を表す。
【００７４】
本発明の好ましい実施形態によるコンピュータプログラムでは、前記のソフトウェアモジュールは次のもののうちの少なくとも１つを備える：
−心臓血管の動きを表す信号をローパスフィルタにかけ、そのローパスフィルタにかけられた信号からＡＯピークを検出する、プログラミング可能なプロセッサを制御するためのコンピュータが実行可能な命令、及び／又は
−心臓血管の動きを表す前記信号をバンドパスフィルタにかけ、そのバンドパスフィルタにかけられた信号からＡＣピークを検出する、プログラミング可能なプロセッサを制御するためのコンピュータが実行可能な命令。
【００７５】
前記のソフトウェアモジュールは、例えば、適切なプログラミング言語が実装されており、且つ、そのプログラミング言語と前記のプログラミング可能なプロセッサに適切なコンパイラーが実装されているサブルーチン又はファンクションであり得る。
【００７６】
本発明の好ましい実施形態によるコンピュータプログラム製品は発明の実施形態によるコンピュータプログラムが符号化されているコンピュータ可読媒体、例えばコンパクトディスク「ＣＤ」を備える。
【００７７】
本発明の好ましい実施形態による信号は発明の実施形態によるコンピュータプログラムを定義する情報を担持するために符号化される。
【００７８】
上記の説明において提供された特定の例は、添付されている特許請求の範囲及び／又は特許請求項の適用可能性を限定するものと解釈されるべきではない。さらに、多くの場合、本発明は心臓の機能不全と異常を検出するための他の技術と共に使用され得ることにも留意されたい。
【図１ａ】
2018086278.tif 000015
【図１ｂ】
2018086278.tif 000016
【図２ａ】
2018086278.tif 000017
【図２ｂ】
2018086278.tif 000018
【図３ａ−３ｂ】
2018086278.tif 000019
【図４】
2018086278.tif 000020
【図５】
2018086278.tif 000021
【手続補正書】
【提出日】平成30年2月15日(2018.2.15)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
心臓活動に関連する電磁現象を表す第１信号と心臓血管の動きを表す第２信号とを受信するための信号インターフェース（４０１）と、
心室性筋肉組織の脱分極によって生じ、前記第１信号上で前記心拍数を繰り返す前記第１波形を表すＲピークを前記第１信号から抽出し、且つ、大動脈弁の開放によって生じ、前記第２信号上で前記心拍数を繰り返す前記第２波形を表すＡＯピークを前記第２信号から抽出するように構成されている処理装置（４０２）と、
を備える器械であって、
前記処理装置が、さらに、
−各タイミング値がある１つの心拍周期に属する前記第１波形の基準点から同一の心拍周期に属する前記第２波形の基準点までの時間を表すタイミングデータを作成し、且つ、
−ペーシングデータの各値が前記心拍数の瞬時値を表すように当該ペーシングデータを作成し、
−前記タイミングデータと前記ペーシングデータとの間の相関係数を計算し、
−前記相関係数が閾値よりも大きい状況に応答して心房細動を表す信号を発信する、ように構成されており、
前記相関係数が、
【数１】
2018086278.tif 000022
であり、
式中、
Ｃ（ｊ）が前記相関係数であり、
Ｅが期待値演算子であり、
ＴＤが前記タイミングデータであり、
μTが前記タイミングデータの平均であり、
ＰＤが前記ペーシングデータであり、
μPが前記ペーシングデータの平均であり、
ｊが心拍周期内の前記タイミングデータに関して前記ペーシングデータの時間差を表す整数である、
ことを特徴とする器械。
【請求項２】
前記処理装置が、Ｒピークの頂点を前記第１波形の前記基準点として使用し、且つ、ＡＯピークの頂点を前記第２波形の前記基準点として使用するように構成されていることを特徴とする請求項１に記載の器械。
【請求項３】
前記処理装置が、前記心臓血管の動きを表す前記第２信号をローパスフィルタにかけ、且つ、前記ローパスフィルタにかけられた信号からＡＯピークを検出するように構成されていることを特徴とする請求項１又は２に記載の器械。
【請求項４】
前記ローパスフィルタリングの上限の周波数が３０Ｈｚであることを特徴とする請求項３に記載の器械。
【請求項５】
前記処理装置が０を前記閾値として使用するように構成されている、ことを特徴とする請求項１〜４のいずれか一項に記載の器械。
【請求項６】
心臓活動に関連する電磁現象を表す第１信号から、心室性筋肉組織の脱分極によって生じ且つ前記第１信号上で心拍数を繰り返す第１波形を表すＲピークを抽出するステップ（２０１）と、
心臓血管の動きを表す第２信号から、大動脈弁の開放によって生じ且つ前記第２信号上で前記心拍数を繰り返す第２波形を表すＡＯピークを抽出するステップ（２０２）と、
を含む方法であって、
当該方法が、
−各タイミング値がある１つの心拍周期に属する前記第１波形の基準点から同一の心拍周期に属する前記第２波形の基準点までの時間を表すタイミングデータを作成するステップ（２０３）と、
−ペーシングデータの各値が前記心拍数の瞬時値を表すように当該ペーシングデータを作成するステップと、
−前記タイミングデータと前記ペーシングデータとの間の相関係数を計算するステップ（２０４）と、
−前記相関係数が閾値よりも大きい状況に応答して心房細動を表す信号を発信するステップと、
をさらに含み、
前記相関係数が、
【数２】
2018086278.tif 000023
であり、
式中、
Ｃ（ｊ）が前記相関係数であり、
Ｅが期待値演算子であり、
ＴＤが前記タイミングデータであり、
μTが前記タイミングデータの平均であり、
ＰＤが前記ペーシングデータであり、
μPが前記ペーシングデータの平均であり、
ｊが心拍周期内の前記タイミングデータに関して前記ペーシングデータの時間差を表す整数である、
ことを特徴とする方法。
【請求項７】
Ｒピークの頂点が前記第１波形の前記基準点であり、且つ、ＡＯピークの頂点が前記第２波形の前記基準点であることを特徴とする請求項６に記載の方法。
【請求項８】
前記方法が、前記心臓血管の動きを表す前記第２信号をローパスフィルタにかけるステップ（２１１）と、前記ローパスフィルタにかけられた信号からＡＯピークを検出するステップ（２１２）とを含むことを特徴とする請求項６又は７に記載の方法。
【請求項９】
前記ローパスフィルタリングの上限の周波数が３０Ｈｚであることを特徴とする請求項８に記載の方法。
【請求項１０】
前記閾値が０である請求項６〜９のいずれか一項に記載の方法。
【請求項１１】
心臓活動に関連する電磁現象を表す第１信号から、心室性筋肉組織の脱分極によって生じ且つ前記第１信号上で心拍数を繰り返す第１波形を表すＲピークを抽出し、
心臓血管の動きを表す第２信号から、大動脈弁の開放によって生じ且つ前記第２信号上で前記心拍数を繰り返す第２波形を表すＡＯピークを抽出する、プログラミング可能なプロセッサを制御するためのコンピュータが実行可能な命令を含むコンピュータプログラムにおいて、
当該コンピュータプログラムが、
−各タイミング値がある１つの心拍周期に属する前記第１波形の基準点から同一の心拍周期に属する前記第２波形の基準点までの時間を表すタイミングデータを作成し、
−ペーシングデータの各値が前記心拍数の瞬時値を表すように当該ペーシングデータを作成し、
−前記タイミングデータと前記ペーシングデータとの間の相関係数を計算し、
−前記相関係数が閾値よりも大きい状況に応答して心房細動を表す信号を発信する、プログラミング可能なプロセッサを制御するためのコンピュータが実行可能な命令をさらに含み、
前記相関係数が、
【数３】
2018086278.tif 000024
であり、
式中、
Ｃ（ｊ）が前記相関係数であり、
Ｅが期待値演算子であり、
ＴＤが前記タイミングデータであり、
μTが前記タイミングデータの平均であり、
ＰＤが前記ペーシングデータであり、
μPが前記ペーシングデータの平均であり、
ｊが心拍周期内の前記タイミングデータに関して前記ペーシングデータの時間差を表す整数である、
ことを特徴とするコンピュータプログラム。
【外国語明細書】
2018086278000001.pdf