JP2016530053A - 生物学的調律障害の分析のために信号セグメントを選択するシステム及び方法 - Google Patents

生物学的調律障害の分析のために信号セグメントを選択するシステム及び方法 Download PDF

Info

Publication number
JP2016530053A
JP2016530053A JP2016542819A JP2016542819A JP2016530053A JP 2016530053 A JP2016530053 A JP 2016530053A JP 2016542819 A JP2016542819 A JP 2016542819A JP 2016542819 A JP2016542819 A JP 2016542819A JP 2016530053 A JP2016530053 A JP 2016530053A
Authority
JP
Japan
Prior art keywords
signal
segment
correlation coefficient
msec
cardiac
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
JP2016542819A
Other languages
English (en)
Other versions
JP2016530053A5 (ja
Inventor
ケアリー ロバート ブリッグス
ケアリー ロバート ブリッグス
ルチル セラ
ルチル セラ
Original Assignee
トペラ インコーポレイテッド
トペラ インコーポレイテッド
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by トペラ インコーポレイテッド, トペラ インコーポレイテッド filed Critical トペラ インコーポレイテッド
Publication of JP2016530053A publication Critical patent/JP2016530053A/ja
Publication of JP2016530053A5 publication Critical patent/JP2016530053A5/ja
Ceased legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7246Details of waveform analysis using correlation, e.g. template matching or determination of similarity
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/363Detecting tachycardia or bradycardia
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2505/00Evaluating, monitoring or diagnosing in the context of a particular type of medical care
    • A61B2505/05Surgical care
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/25Bioelectric electrodes therefor
    • A61B5/279Bioelectric electrodes therefor specially adapted for particular uses
    • A61B5/28Bioelectric electrodes therefor specially adapted for particular uses for electrocardiography [ECG]
    • A61B5/283Invasive
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/25Bioelectric electrodes therefor
    • A61B5/279Bioelectric electrodes therefor specially adapted for particular uses
    • A61B5/28Bioelectric electrodes therefor specially adapted for particular uses for electrocardiography [ECG]
    • A61B5/283Invasive
    • A61B5/287Holders for multiple electrodes, e.g. electrode catheters for electrophysiological study [EPS]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6846Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be brought in contact with an internal body part, i.e. invasive
    • A61B5/6847Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be brought in contact with an internal body part, i.e. invasive mounted on an invasive device
    • A61B5/6852Catheters
    • A61B5/6858Catheters with a distal basket, e.g. expandable basket
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Cardiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Psychiatry (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (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)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

心臓信号のセグメントを処理するシステム及び方法。第1の信号中の複数のセグメントの各々を複数のオフセットで相関させて、第1の信号中の複数のセグメントの各々に関して最も高い相関係数を決定する。第2の信号中の複数のセグメントの各々を複数のオフセットで相関させて、第2の信号中の複数のセグメントの各々に関して最も高い相関係数を決定する。第1の信号のセグメント及び第2の信号のセグメントに関する最も高い相関係数を用いて、複数の複合相関係数を決定する。第1の信号のセグメントと第2の信号のセグメントは、ほぼ同時期である。第1の信号からのセグメント及び第2の信号からのセグメントを含む1組のセグメントが選択される。この組が、複数の複合相関係数からの最も高い複合相関係数に関連づけられる。

Description

本願は、一般に、生物学的調律(リズム)障害に関する。特に、本発明は、生物学的調律障害(例えば、心調律障害)の分析のための信号セグメントを選択するシステム及び方法に関する。
〔関連出願の説明〕
本願は、2013年9月12日に出願された米国特許出願第61/877,093号及び2014年9月11日に出願された米国特許出願第14/483,914号の優先権主張出願であり、これら米国特許出願の両方を参照により引用し、これらの記載内容を本明細書の一部とする。
生物学的調律障害、例えば心(心臓)調律(「心臓リズム」又は「心拍リズム」)障害は、よく見受けられ、世界中において罹患率及び死亡のかなりの原因となっている。心臓内の電気系の機能不全は、心調律障害の誘因となっている。心調律障害は、多くの形態で存在し、これら形態のうちの最も複雑で且つ治療困難な形態は、心房細動(AF)、心室性頻拍(VT)及び心室細動(VF)である。他の調律障害は、治療するのが比較的簡単であるが、臨床的に重要である場合があり、かかる調律障害としては、心房性頻拍(AT)、上室性頻拍(SVT)、心房粗動(AFL)、上心室性異所性棘波群/拍動(SVE)及び心室性期外(早期心室性)棘波群/拍動(PVC)が挙げられる。
従来、心調律障害−−特にAF、VF及び多形VTの複合性調律障害−−の治療は、困難であった。というのは、発生源が隠されている心臓内の場所を突き止めることができなかったからである。複合性調律障害がどのように機能するか及びこれら複合性調律障害を治療する臨床的な応用については、種々の理論が存在する。しかしながら、これら応用のうちで複合性調律障害の治療において実りをもたらしたものはなかった。
最近、複合性心調律障害と関連した発生源を初めて突き止めた(特定した)という画期的な発見があった。この画期的な技術は、世界中において心調律障害のうちの高いパーセンテージを生じさせ持続させている旋回賦活パターン(旋回源、例えばロータ(rotor ))並びに病巣源を突き止めるために患者の心臓内に導入されたカテーテルの電極から得られる信号において、心臓賦活情報(賦活開始時期)を首尾良く分析して再構成した。かくして、心調律障害の治療は、心調律障害をなくすために、患者の心臓内の旋回源及び/又は病巣源に的が絞られていると言える。かかる治療は、例えばアブレーション(焼灼)によって首尾良く実行できる。
上述したように、心臓信号は、一般に、患者の心臓内に導入されたカテーテルの電極から得られる(例えば、検出され、収集され、又は記録される)。多くのノイズ発生源は、信号が患者から得られるときに当該信号中に埋もれている又は重畳されている場合が多い。これらノイズ源としては、患者の心臓の別の部分からの電気的活動、患者の他の解剖学的構造、電極の動き及び/又は患者の動き(例えば、呼吸)に起因する運動アーチファクト、電極が互いに接することによって生じる機械的クロストーク、電子増幅器の飽和、外部システムからの無線周波(RF)エネルギー、並びに他のノイズ源が挙げられる。加うるに、電極は、信号の振幅を減少させ得ると共に最悪のケースでは結果的に電気的活動が信号中に存在しなくなり得る、患者の心臓との種々の接触レベル(又は、非接触レベル)を有する場合がある。
心臓賦活情報(賦活開始時期)の再構築は、コンピュータ計算集約的であると共に時間集約的であり得る信号の分析を必要とする。信号の全体を分析することは、これらのコンピュータ計算上の制約又は他のコンピュータ計算上の制約に起因して、有利とは言えない場合がある。さらに、これら信号の或る特定の部分、場合によって大部分は、ノイズによって悪影響を受け得る。かかる状況では、信号中に重畳されたノイズによって分析が複雑で制限されるような信号の部分を回避することが、有利な場合がある。
これら信号の全体の分析は、心調律障害の発生源を突き止める際の時間及び精度、並びに、治療及び除去のために心調律障害の発生源を標的にする際の精度、にマイナスの影響を及ぼす場合がある。
したがって、更なる分析のためにノイズの量が低減されている周期的な心臓情報を含むこれら信号の部分(例えば、信号セグメント)を識別する(特定する)ことが望ましい。これは、心調律障害の発生源を突き止める際の時間及び精度、並びに、治療及び除去のために心調律障害の発生源を標的にする際の時間及び精度、を向上させることができる。
本願の要旨は、患者の調律障害の原因又は発生源の検出及び治療のため、最小ノイズと安定した明確に規定される生物学的活動とを含み患者から得られる複数の信号から1組の信号セグメントを選択することに適用できる。信号セグメント選択は、心調律障害を含む種々の調律障害並びに他の生物学的調律障害、例えば神経性発作、食道痙攣、膀胱不安定性、過敏性腸管症候群、並びに、障害の原因又は発生源の判定、診断及び/又は治療を可能にするよう複数の信号から信号セグメントが選択され得る対象としての他の生物学的障害、に適用できる。しかしながら、本願の要旨は、結果的に複雑な賦活パターンをもたらす複合性調律障害に特に有用であり、特に心臓の複合性調律障害に有用であり、その目的は、疾患の原因又は発生源の検出を改善してかかる障害を都合良く治療できるようにすることにある。
複合性心調律障害においては、典型的には、賦活パターンは解読するのが極めて困難であり、複合性障害における心拍動の正確な賦活に関する情報を求めることは、従来可能ではなかった。本発明の利点のうちの1つは、最小ノイズと安定した明確に規定される生物学的活動とを含む複数の信号から信号セグメントを選択する能力であり、その結果、調律障害の原因及び/又は発生源の突き止めを行うことができ、そして調律障害をより迅速に且つより高い精度で治療することができる。別の利点は、検出装置、例えばセンサが設けられているカテーテル、を患者の体内又はその近くで用いている状態で迅速に実施することができ、信号セグメント選択、選択した信号セグメントに基づく調律障害の発生源の突き止め、及び、心組織の治療を行って調律障害を改善して多くの場合障害を治癒させること、に続くことができるシステム及び方法を提供することにある。
一実施形態又は一観点によれば、心臓信号のセグメントを処理する方法が開示される。
少なくとも1つの第1の心臓信号中の信号セグメントが、複数の時間オフセットによってシフトされた当該信号セグメントと相関されて、当該信号セグメントに関する最も高い関係数が決定される。
少なくとも1つの第2の心臓信号中の信号セグメントが、複数の時間オフセットによってシフトされた当該信号セグメントと相関されて、当該信号セグメントに関する最も高い相関係数が決定される。
前記少なくとも1つの第1の心臓信号及び前記少なくとも1つの第2の心臓信号中において追加の信号セグメントに関する相関工程が繰り返されて、追加の信号セグメントの各々に関する最も高い相関係数が決定される。
前記少なくとも1つの第1の心臓信号中の信号セグメント及び前記少なくとも1つの第2の心臓信号中の信号セグメントに関する最も高い相関係数を用いて、複数の複合相関係数が生じさせられる。前記少なくとも1つの第1の心臓信号中の信号セグメントは、前記少なくとも1つの第2の心臓信号中の信号セグメントと、ほぼ同時期に存在する。
前記少なくとも1つの第1の心臓信号からの少なくとも1つの信号セグメント及び前記少なくとも1つの第2の心臓信号からの少なくとも1つの信号セグメントを含む1組の信号セグメントが選択される。当該1組の信号セグメントが、複数の複合相関係数のうちで最も高い複合相関係数と関連付けられる。
この方法によれば、当該1組の信号セグメントを用いて、調律障害の発生源を突き止めることができる。
別の実施形態又は別の観点によれば、心臓信号のセグメントを処理するシステムが開示される。このシステムは、処理装置と、記憶装置と、を含む。当該記憶装置は、処理装置によって実行される時に、当該処理装置が以下のような演算を実行する、という命令を記憶している。
少なくとも1つの第1の心臓信号中の信号セグメントが、複数の時間オフセットによってシフトされた当該信号セグメントと相関されて、当該信号セグメントに関する最も高い関係数が決定される。
少なくとも1つの第2の心臓信号中の信号セグメントが、複数の時間オフセットによってシフトされた当該信号セグメントと相関されて、当該信号セグメントに関する最も高い相関係数が決定される。
前記少なくとも1つの第1の心臓信号及び前記少なくとも1つの第2の心臓信号中において追加の信号セグメントに関する相関工程が繰り返されて、追加の信号セグメントの各々に関する最も高い相関係数が決定される。
前記少なくとも1つの第1の心臓信号中の信号セグメント及び前記少なくとも1つの第2の心臓信号中の信号セグメントに関する最も高い相関係数を用いて、複数の複合相関係数が生じさせられる。前記少なくとも1つの第1の心臓信号中の信号セグメントは、前記少なくとも1つの第2の心臓信号中の信号セグメントと、ほぼ同時期に存在する。
前記少なくとも1つの第1の心臓信号からの少なくとも1つの信号セグメント及び前記少なくとも1つの第2の心臓信号からの少なくとも1つの信号セグメントを含む1組の信号セグメントが選択される。当該1組の信号セグメントが、複数の複合相関係数のうちで最も高い複合相関係数と関連付けられる。
これら演算は、調律障害の発生源を突き止めるために、当該1組の信号セグメントを用いる工程を更に含むのが良い。
別の実施形態又は別の観点によれば、複数の心臓信号のセグメントを処理する方法が開示される。
最も高い相関係数を決定するために、オフセットの第1の範囲に従って、複数の信号の中からほぼ同時期に存在する信号セグメントが自己相関させられる。最も高い相関係数を用いて、複合相関係数が生じさせられる。
しきい値よりも低い複合相関係数のうちの1つ又は2つ以上が決定される。ほぼ同時期に存在する信号セグメントに関する最も高い相関係数を決定するために、オフセットの第2の範囲に従って、決定された1つ又は2つ以上の複合相関係数と関連付けられたほぼ同時期に存在する信号セグメントが自己相関される。第2のオフセット範囲は、第1のオフセット範囲とは異なっている。
選択されたオフセットの範囲に関して最も高い相関係数と関連付けられる1組のほぼ同時期に存在する信号セグメントが、選択される。
さらに別の実施形態又は別の観点によれば、心臓信号のセグメントを処理するシステムが開示される。このシステムは、処理装置と、記憶装置と、を含む。当該記憶装置は、処理装置によって実行される時に、当該処理装置が以下のような演算を実行する、という命令を記憶している。
最も高い相関係数を決定するために、オフセットの第1の範囲に従って、複数の信号の中からほぼ同時期に存在する信号セグメントが自己相関させられる。最も高い相関係数を用いて、複合相関係数が生じさせられる。
しきい値よりも低い複合相関係数のうちの1つ又は2つ以上が決定される。ほぼ同時期に存在する信号セグメントに関する最も高い相関係数を決定するために、オフセットの第2の範囲に従って、決定された1つ又は2つ以上の複合相関係数と関連付けられたほぼ同時期に存在する信号セグメントが自己相関される。第2のオフセット範囲は、第1のオフセット範囲とは異なっている。
選択されたオフセットの範囲に関して最も高い相関係数と関連付けられる1組のほぼ同時期に存在する信号セグメントが、選択される。
本願のこれらの目的、目標及び利点並びに他の目的、目標及び利点は、添付の図面と関連して以下の詳細な説明を読めば明らかになろう。
幾つかの実施形態又は幾つかの観点が、添付の図面の各図に、例示として、本発明を限定しないで、図示されている。
例示の信号処理システムを示す図である。
図1の患者の体内の検出場所のところに位置決めされた例示のセンサからほぼ同時期に得られる複合性心調律障害の例示の電気信号図である。
例示として3つのセグメントにセグメント化された、図2の例示の第1の信号を示す図である。
例示として3つのセグメントにセグメント化された、図2の例示の第2の信号を示す図である。
例示として3つのセグメントにセグメント化された、図2の例示の第3の信号を示す図である。
図2の第1の信号のうちの図3の第1の例示の信号セグメント中の、3つの例示の時間シフト(又はオフセット)を示す図である。
図2の第1の信号のうちの図3の第1の信号セグメント中の、100msecから350msecまでの範囲の時間オフセットに関して求められた相関係数(r)の曲線を示す例示のグラフ図である。
第1の実施形態に従って、生物学的調律障害、例えば心調律障害、の分析のために信号セグメントを選択する例示の方法を示す流れ図である。
所与の時間オフセット、例えば100msec、についての例示の相関計算式を示す図である。
図2の信号の各々の信号セグメントに関する最も高い相関係数(r)の概要を示す例示の表である。
最も高い相関係数(r)の二乗平均平方根に関する例示の計算式を示す図である。
図1の心臓内の検出場所のところに位置決めされた例示のセンサと関連して得られた、複合性心調律障害の例示の信号を示す図である。
第1の選択されたオフセット範囲に従った、図12に示された信号の信号セグメント中の例示の時間シフト(又はオフセット)を示す図である。
図12に示された信号の信号セグメントにおいて、第1の選択されたオフセット範囲に従った時間オフセットに関して求められた相関係数(r)の曲線を示す例示のグラフ図である。
第2の選択されたオフセット範囲に従った、図12に示されている信号の信号セグメント中の例示の時間シフト(又はオフセット)を示す図である。
図12に示された信号の信号セグメントにおいて、第2の選択されたオフセット範囲に従った時間オフセットに関して求められた相関係数(r)の曲線を示す例示のグラフ図である。
第3の選択されたオフセット範囲に従った、図12に示されている信号の信号セグメント中の例示の時間シフト(又はオフセット)を示す図である。
図12に示された信号の信号セグメント中において、第3の選択されたオフセット範囲に従った時間オフセットに関して求められた相関係数(r)の曲線を示す例示のグラフ図である。
第2の実施形態に従って、生物学的調律障害、例えば心調律障害、の分析のために信号セグメントを選択する例示の方法を示す流れ図である。
図1の心臓内の検出場所のところに位置決めされた例示のセンサと関連して得られた、複合性心調律障害の例示の信号を示す図である。
図20に示された信号を含む複数の信号のセグメントと関連付けられた複合相関係数の例示の曲線を示す図である。
図20に示された信号を含む複数の信号の同時期に存在するセグメントを考慮した、心調律障害のサイクル長(CL)と関連付けられた複合オフセットの例示の曲線を示す図である。
汎用コンピュータシステムの例示の実施形態のブロック図である。
生物学的調律障害(例えば、心調律障害)の分析のために信号セグメントを選択するシステム及び方法が、本明細書において開示される。以下の説明において、説明の目的上、多くの特定の細部について説明されるが、その目的は、例示の実施形態又は観点の完全な理解を提供することにある。しかしながら、例示の実施形態は、開示される特定の細部の全てを備えなくとも実施できる、ということが当業者には明らかであろう。
図1は、例示の信号処理システム100を示している。例示のシステム100は、心調律障害と関連した患者の心臓から信号を検出し、低減された量のノイズを伴う更なる分析のための周期的な心臓賦活情報を含むこれら信号から、心調律障害の発生源を突き止める際の時間及び精度並びに治療及び除去のために心調律障害の発生源を標的にする際の時間及び精度を向上させることができるような、信号セグメントを選択するよう構成されている。心臓は、右心房122、左心房124、右心室126及び左心室128を有する。例示のシステム100は、カテーテル101、信号処理装置114及びコンピュータ計算装置116を含む。
カテーテル101は、心臓内の心臓賦活情報を検出し、検出した心臓賦活情報をワイヤレス接続方式かワイヤード接続方式かのいずれかにより信号処理装置114に伝送するよう構成されている。カテーテルは、患者の血管を通って心臓内に挿入できる複数の電極又はセンサ102〜112を有する。
幾つかの実施形態又は幾つかの観点では、センサ102〜112のうちの1つ又は2つ以上は、患者の心臓内に挿入されない。例えば、幾つかのセンサは、患者の表面(例えば、心電図、体表面マッピング)により又は患者と接触しないで遠隔で(例えば、心磁図)心臓賦活状態を検出することができる。別の例として、幾つかのセンサは又、非電気的検出装置(例えば、心エコー図)の心臓の動きから心臓賦活情報を導き出すことができる。種々の実施形態又は種々の観点では、これらセンサは、別々に又は種々の組み合わせ状態で使用でき、更に、これら別々の又は種々の組み合わせ状態は又、患者の心臓内に挿入されたセンサと組み合わせて使用できる。
検討対象の心臓内の或る特定の場所(センサ配置場所)に位置決めされたセンサ102〜112は、検出場所の心臓賦活情報を検出することができ、そして更に、エネルギーを送ってこれら検出場所で心臓を焼灼(ablation)することができる。検出場所は、センサが心臓賦活情報を検出する元となるセンサ配置場所の近くに位置すると共に当該センサ配置場所を含む領域である。注目されるように、センサ102〜112は又、心臓のオーバーラップした領域(例えば、右心房122と左心房124)からの心臓賦活情報を検出することもできる。
信号処理装置114は、検出場所のところでセンサ102〜112により検出された心臓賦活情報を処理して(例えば、明確にすると共に増幅して)電気図信号(以下、「信号」という)の状態にすると共に、処理した信号を本明細書において開示するセグメント選択及び更なる分析のためにコンピュータ計算装置116に提供するよう構成されている。センサ102〜112からの心臓賦活情報を処理する際、信号処理装置114は、心臓120のオーバーラップした領域からの心臓賦活情報を減算して、処理後の信号を分析のためにコンピュータ計算装置116に提供することができる。幾つかの実施形態又は幾つかの観点では、信号処理装置114は、単極信号を提供するよう構成されているが、他の実施形態又は他の観点では、信号処理装置114は、双極信号を提供することができる。
上述したように、センサ102〜112によって検出される賦活情報は、埋め込まれた又は重ね合わされた(重畳された)ノイズを含む場合がある。コンピュータ計算装置116は、信号処理装置114から信号を受け取り又はこれにアクセスするよう構成され、幾つかの実施形態では、コンピュータ計算装置116は、受け取った又はアクセスした信号から或る特定の長さのエポック(窓)を選択するよう構成されている。かかる信号も又、本明細書においては信号であるとみなされる。コンピュータ計算装置116は、更に、心調律障害の発生源を突き止めるために、開示される方法、機能又は論理に従って、低減された量のノイズを伴った状態で信号のセグメント(以下、「信号セグメント」という)を自動的に選択するよう構成されている。具体的に説明すれば、当該選択した信号セグメントを用いると、心調律障害の発生源を突き止める時間及び精度、並びに、治療及び除去のために心調律障害の発生源を標的にする際の時間及び精度を向上させることができる。例えば、選択された信号セグメントは、ブリッグス等(Briggs et al. )に付与された米国特許第8,165,666号の発明によって入力信号として処理することができる。この米国特許は、当該参照により引用され、その記載内容全体が本明細書の一部とされる。
図2は、第1の実施形態に従って、図1の心臓120内の検出場所に位置決めされた例示のセンサ102〜106と関連してほぼ同時期に(例えば、同時に又はほぼ同時に)得られる複合性心調律障害の例示の信号202〜206を示している。例えば、カテーテル101のセンサ102〜106が、図1に示されているように、右心房122内の検出場所に位置決めされ得る。一例として、心調律障害は、複合性調律障害AF、VF及び多形VT又は別の心調律障害である場合があり得る。
コンピュータ計算装置116は、例示の信号202〜206を、受け取るか、これらにアクセスするか、違ったやり方で選択する。例示の信号202〜206は、簡潔ではあるが明確な図示を提供するよう、長さが8秒であるよう図示されている。しかしながら、信号202〜206は、60秒とか、これよりも長い又は短い時間とか、であっても良い。しかしながら、コンピュータ計算装置116によって受け取られ又はアクセスされる信号は、任意の長さのものであって良く、例えば、離散した長さを有しても良く、或いは長さが連続していても良い。幾つかの実施形態では、信号202〜206は、かかる受け取った又はアクセスした信号からコンピュータ計算装置116によって選択された或る特定の時間長さのエポック(窓)を表すことができる。例えば、このようにして選択された信号202〜206は、60秒であっても良く、或いは、これよりも長い又は短い期間のものであっても良い。信号202〜206は、例えば図1に示されているように、複合性心調律障害中の心臓120内の検出場所に位置決めされた例示のセンサ102〜106と関連している。複数の信号がカテーテル101によって提供され得て、例えば、バスケットカテーテルの対応のセンサからの64個の信号が、心臓120の右心房122のパノラマ(写景)図を提供することができる。
或る特定の心調律障害(例えば、AF)における例示の観察に基づき、賦活開始時期の相互間の約100msec(ミリ秒)〜約300msecのサイクル長が、かかる心調律障害中に特定のセンサ(検出場所)のところで信号中に典型的に観察される。或る特定の他の心調律障害(例えば、AT、AFL)における例示の観察に基づくと、サイクル長は、一般に、これよりも長い(例えば、約300msec〜約500msec)。正常な洞調律(調律)では、作動開始時期相互間のサイクル長は、一般に、更に長い(例えば、約600msec〜約1000msec)。この例示の実施形態では、信号202〜206は、100msecから350msecまでの範囲のサイクル長と関連している。
信号中の賦活開始時期は、一般に、当該信号中の局所活動と間違えられ得るファーフィールドアーチファクト(far-field artifacts)をほとんど伴わない、検出場所に対してローカルである信号中に重畳されたベースラインからの信号の僅かな程度の変化を有するものとして、識別可能である。検出場所での局所活動は、信号中に、鮮明な変曲点及び急な勾配を有する賦活開始時期を呈し、典型的には例示の観察のうちの幾つか(例えば、心房細動)において約100msec〜250msec間続き、他の例示の観察(例えば、心房性頻拍)ではこれよりも長く続く、再分極を表すなだらかで偏差の少ない勾配の期間を続いて呈する、ということで特徴付けられ得る。
短い持続時間(例えば、約100msecの最も短いサイクル長よりも短い)の複数の偏向は、ファーフィールド賦活又は単なるノイズに対立するものとしての検出場所での局所賦活開始時期の識別を極めて困難にする場合がある。例えば、Sig.1(符号202)の最初の1,000msecでは、8つ程度の偏向が存在し、これら偏向は、場合によっては、賦活開始時期を示す場合がある。しかしながら、心臓120は、約100msec〜約300msecよりも短いサイクル長で検出場所で検出された賦活開始時期(偏向)後に再び生理学的に活動することができない。というのは、検出場所に対してローカルな心臓120の組織は、再分極を行わなければならないからである。偏向は、当該偏向が、当該検出場所の隣に位置する検出場所で検出された信号内にも有意に存在している時、当該検出場所に対してローカルでない。
さらに、信号202〜206の全体の処理は、極めて時間集約的であり且つコンピュータ計算集約的であり得て、心臓120の複合性心調律障害の発生源を識別するための時間を増大させ得ると共にこの能力を制限し得る。かくして、更なる分析のためにファーフィールド賦活及びノイズの量が減少した周期的心臓賦活情報を含む信号202〜206中の対応の信号セグメントを選択することが望ましい。それにより、心調律障害の発生源を突き止める際の時間及び精度、並びに、治療及び除去のために心調律障害の発生源を標的にする際の時間及び精度を向上させることができる。
本明細書において説明するシステム及び方法は、選択された時間オフセットの範囲(例えば、100msec〜350msec)にわたっての相関計算(例えば、自己相関)を用いて、信号202〜206における信号セグメント選択を自動化して、心調律障害の発生源の次の分析及び識別のために当該信号202〜206中の複数の信号セグメントから全体として最も相関のある信号202〜206内の1組の対応の信号セグメントをプログラム的に決定する。
注目されるように、選択された時間オフセットの範囲は、例示のサイクル長(例えば、この実施形態では、関連の複合性調律障害の説明のために用いられている300msec)よりも僅かに広く(例えば、350msec)、そのことが、例示のサイクル長(例えば、300msec)の終わりの近く又はその付近での相関を見落とす可能性を減少させる。同様に、例示のサイクル長(例えば、100msec)の開始時の近く又はその付近での相関を見落とす可能性を減少させるために、選択された時間オフセットの範囲は、開始時のところでは僅かに広い(例えば、50msec)よう調節されるのが良い。幾つかの場合、選択された時間オフセットの範囲は、この実施形態の心調律障害と関連した上述の例示のサイクル長(例えば、100msec〜300msec)に等しいかこれに近似している場合がある。
さらに、時間オフセットの範囲は、1つ又は2つ以上の特定の心調律障害(例えば、AF、VT、VF)並びに1つ又は2つ以上の他の生物学的調律障害(例えば、神経性発作、食道痙攣、膀胱不安定性、過敏性腸管症候群、及び他の生物学的障害)に基づいて選択的に変化させることができる。
この例示の実施形態では、相関計算のための最小及び最大時間オフセットは、それぞれ、100msec及び350msecに選択されている。最小及び最大オフセットは、この実施形態の例えば上述の複合性調律障害の100msec及び300msecの外側サイクル長に近似する。その結果、相関係数(r)に寄与する信号セグメントの周期性は、約100msec〜約300msecの間のサイクル長を有することになる。
例示の相関計算について、図9を参照して本明細書において詳細に説明する。この時点においては、以下を言えば十分である。すなわち、当該相関計算は、時間オフセットの範囲にわたって、例えば100msec〜350msecにわたって、各信号セグメントをそれ自体と相関させて(例えば、自動相関させて)、信号セグメントの複数のオフセットでの複数の相関係数(r)を求める。当該複数の相関係数から最も高い(最大の)相関係数(r)が、当該信号セグメントのために記録される。この実施形態の例示の時間オフセット範囲及び対応の相関係数(r)について、図6及び図7を参照して、本明細書において詳細に説明する。異なる相関方法又は技術を用いて、時間オフセット範囲にわたって各信号セグメントをそれ自体と相関させる(例えば、自己相関させる)こともできる。
最小ノイズ及び明確に規定される心臓活動を含む信号セグメントは、心臓活動にわたって重畳された相当な量のノイズを含む信号セグメントよりも高い相関係数(r)を有することになる。したがって、全ての信号202〜206全体にわたって最も高い複合相関係数(r)を有する信号202〜206中の対応の(ほぼ同時期に存在する)信号セグメントの選択は、その下にある心臓活動情報(賦活開始時期)を検出する更なる分析のための信号202〜206の全てのうちで最良の信号品質を有する1組の信号セグメントを提供することになる。それにより、心調律障害の発生源を突き止める際の時間及び精度、並びに、治療及び除去のために心調律障害の発生源を標的にする際の時間及び精度を向上させることができる。
図3は、例示の3つのセグメント302〜306にセグメント化された図2の例示の第1の信号202を示している。説明を簡潔にすると共に分かりやすくするために、3つのセグメント302〜306しか図示されていない。図2を参照して注目されるように、第1の信号202は、長さが8秒である。当然のことながら、異なる長さ、例えば60秒以上が可能である。
例示の信号セグメント302〜306の各々は、長さ(以下、「セグメント長」という)が2秒である。図3に示されているように、信号セグメント302〜306は、少なくとも部分的にオーバーラップしている場合がある。例えば、信号セグメント302は、2msec〜4msecであり、信号セグメント304は、3msec〜5msecであり、信号セグメント306は、4msec〜6msecである。注目されるように、第1の信号202は、上述のセグメント長、例えば0msec〜2msec、1msec〜3msec、5msec〜7msec及び6msec〜8msec、の追加のオーバーラップ状態の信号セグメントにセグメント化することができる。
さらに注目されるように、信号セグメント302〜306の異なるセグメント長が、例えば第1の信号202の全体的長さに基づいて選択され得る。一例として、第1の信号202は、少なくとも部分的にオーバーラップすることができる長さが4秒の信号セグメント、例えば、0msec〜4msec、1msec〜5msec、2msec〜6msec、3msec〜7msec、4msec〜8msec、にセグメント化することができる。
信号セグメント302〜306の各々が、当該信号セグメント302〜306の各々の複数のオフセットでの複数の相関係数(r)を決定するために、選択された時間オフセット範囲にわたって、それ自体に相関させられる(例えば、図9の相関関係式を用いて)。時間オフセット範囲は、生物学的調律障害(例えば、心調律障害)のサイクル長と関連されている。この実施形態によれば、この範囲は、100msec〜350msecであるように選択される。
複数の相関係数からの最も高い(最大)相関係数(r)が信号セグメント302〜306の各々について記録される。図3に示されているように、信号セグメント302に関する最大相関係数(r=0.73)が、オフセット201msecで得られる。信号セグメント304に関し、最大相関係数(r=0.76)が、オフセット204msecで得られる。最後に、信号セグメント306に関し、最大相関係数(r=0.73)が、オフセット205msecで得られる。
図4は、例示の3つのセグメント402〜406にセグメント化された図2の例示の第2の信号204を示している。説明を簡潔にすると共に分かりやすくするために3つのセグメント402〜406しか図示されていない。例示の信号セグメント402〜406の各々は、2秒のセグメント長を有する。
図4に示されているように、信号セグメント402〜406は、少なくとも部分的にオーバーラップしている場合がある。例えば、信号セグメント402は、2msec〜4msecであり、信号セグメント404は、3msec〜5msecであり、信号セグメント406は、4msec〜6msecである。注目されるように、第2の信号204は、上述のセグメント長、例えば、0msec〜2msec、1msec〜3msec、5msec〜7msec及び6msec〜8msec、の追加の少なくとも部分的にオーバーラップ状態の信号セグメントにセグメント化することができる。
図3の第1の信号202を参照して説明したように、信号セグメント402〜406の異なるセグメント長を選択することもできる。例えば、第2の信号204は、少なくとも部分的にオーバーラップし得る長さが4秒の信号セグメント、例えば、0msec〜4msec、1msec〜5msec、2msec〜6msec、3msec〜7msec、4msec〜8msecの信号セグメント、にセグメント化することができる。
信号セグメント402〜406の各々が、当該信号セグメント402〜406の各々の複数のオフセットでの複数の相関係数(r)を決定するために、時間オフセット範囲(例えば、100msec〜350msec)にわたって、それ自体に相関させられる(例えば、図9の相関関係式を用いて)。複数の相関係数のうちで最も高い(最大)相関係数(r)が信号セグメント402〜406の各々について記録される。図4に示されているように、信号セグメント402に関する最大相関係数(r=0.81)がオフセット200msecで得られる。信号セグメント404に関し、最大相関係数(r=0.62)がオフセット202msecで得られる。最後に、信号セグメント406に関し、最大相関係数(r=0.54)がオフセット205msecで得られる。
図5は、例示の3つのセグメント502〜506にセグメント化された図2の例示の第3の信号206を示している。説明を簡潔にすると共に分かりやすくするために3つのセグメント502〜506しか図示されていない。例示の信号セグメント502〜506の各々は、2秒のセグメント長を有する。
図5に示されているように、信号セグメント502〜506は、少なくとも部分的にオーバーラップしている場合がある。例えば、信号セグメント502は、2msec〜4msecであり、信号セグメント504は、3msec〜5msecであり、信号セグメント506は、4msec〜6msecである。注目されるように、第3の信号206は、上述のセグメント長、例えば、0msec〜2msec、1msec〜3msec、5msec〜7msec及び6msec〜8msec、の追加の少なくとも部分的にオーバーラップ状態の信号セグメントにセグメント化することができる。
図3の第1の信号202を参照して説明したように、信号セグメント502〜506の異なるセグメント長を選択することもできる。例えば、第3の信号206は、少なくとも部分的にオーバーラップし得る長さが4秒の信号セグメント、例えば、0msec〜4msec、1msec〜5msec、2msec〜6msec、3msec〜7msec、4msec〜8msecの信号セグメント、にセグメント化することができる。
信号セグメント502〜506の各々が、当該信号セグメント502〜506の各々の複数のオフセットでの複数の相関係数(r)を決定するために、時間オフセット範囲(例えば、100msec〜350msec)にわたって、それ自体に相関させられる(例えば、図9の相関関係式を用いて)。複数の相関係数のうちで最も高い(最大)相関係数(r)が、信号セグメント502〜506の各々について記録される。図5に示されているように、信号セグメント502に関する最大相関係数(r=0.43)がオフセット210msecで得られる。信号セグメント504に関し、最大相関係数(r=0.47)がオフセット199msecで得られる。最後に、信号セグメント506に関し、最大相関係数(r=0.39)がオフセット200msecで得られる。
図6は、図2の第1の信号202のうちの図3の第1の例示の信号セグメント302の3つの例示の時間シフト(又はオフセット)602〜606を示す図である。
例えば、第1の時間オフセット602は、100msecであり、第2の時間オフセット604は、200msecであり、第3の時間オフセット602は、350msecである。1msecのインクリメントを用いれば、100msecから350msecまでの範囲において、250個の時間オフセットが存在し得る。
図9の例示の相関計算式を用いて、100msecから350msecまでの範囲の時間オフセットの各々について、相関係数(r)が計算される。例えば、100msecの下時間オフセットでは、相関係数(r)=−0.47であり、350msecの上オフセットでは、相関係数(r)=0.02である。最大相関係数(r)=0.73が、図6に示されているように、200msecのオフセットで得られる。異なる相関方法又は技術を用いて例示の信号セグメント302をそれ自体に相関させる(自己相関させる)こともできる。
一例として、図9の例示の相関関係式を用いて、信号セグメントがそれ自体に相関させられる。具体的に説明すると、当該相関は、2つの別々のデータシリーズ(信号セグメントの各msecでの電圧要素を表す)の積の合計として説明することができる。相関係数(例えば、相関の数値結果)は、範囲−1〜1に標準化される。
一実施形態では、信号セグメント(例えば、信号セグメント302)が所与の場合、元の信号セグメントから2つの基準セグメント、即ち、所与の信号セグメントの開始から取り出されたt個の要素(例えば、1msecの各要素)を含む第1の基準セグメント及び所与の信号セグメントの終わりから取り出されたt個の要素(例えば、1msecの各要素)を含む他の基準セグメント、を生じさせることによって、所与の時間オフセットt(例えば、100msec)での相関を計算することができる。かくして、所与の時間オフセットtに関する相関係数(r)が、基準セグメントに基づき、例えば図9に示されている例示の相関計算式を用いて、求められる。
別の実施形態では、1つの基準信号セグメントが、所与の時間セグメント(例えば、信号セグメント302)について所与の時間オフセットt(例えば、100msec)で生成される。それは、所与の信号セグメントの開始時から、所与の信号セグメントの終わりから所与のオフセットtを差し引いた値によって表される終了時まで、延びるものである。かくして、所与の時間オフセットtに関する相関係数(r)は、時間オフセットから信号セグメントの終わりまで延びる所与の信号セグメント及び基準信号セグメントに基づき、例えば図9に示されている例示の相関関係式を用いて、求められる。
各信号202,204,206の信号セグメント302〜306,402〜406,502〜506が処理されて、信号セグメント302を参照して上述したように、100msecから350msecまでの範囲の時間オフセットの各々について相関係数(r)が計算される。
図7は、図2の第1の信号202の図3に示されている第1の信号セグメント302の100msecから350msecまでの範囲の時間オフセットについて求められた相関係数(r)の曲線702を示す例示のグラフ700である。
図7に示されているように、相関係数は、100msecから350msecまでの時間オフセット範囲において、−0.47〜0.02であり、最大相関係数(r)=0.73は、200msecの時間オフセットで起こる。
同様な曲線を生成させて、信号202の他の信号セグメント304,306並びに他の信号204,206の信号セグメント402〜406,502〜506について、100msecから350msecまでの範囲の時間オフセットについて求められた相関係数(r)を説明することができる。注目すべき重要なことは、最大相関係数(r)が信号(202,204,206)の各々の各信号セグメント(302〜306,402〜406,502〜506)について求められるということである。
図8は、第1の実施形態に従って、生物学的調律障害(例えば心調律障害)の分析のために信号セグメントを選択する例示の方法800を示す流れ図である。例示の方法800は、図1に示されているコンピュータ計算装置116によって実施され得る。
具体的に説明すると、例示の方法800は、演算ないしステップ802で始まり、この演算ないしステップ802では、複数の信号が、心臓120内に設けられたセンサから信号処理装置114を介してコンピュータ計算装置116によって受け取られる、又は、アクセスされる。かかる信号は、特定長さのものであるのが良く、例えば、60秒以上の長さのものである。例えば、信号202〜206は、図1に示されているように、心臓120の右心房122内に設けられたカテーテル101のセンサ102〜106からのものである。幾つかの実施形態又は幾つかの観点では、当該センサからの信号の少なくとも一部分が、信号処理装置114によって記録され得て、次にコンピュータ計算装置116に提供され得る。
演算ないしステップ804では、1つ又は2つ以上のフィルタ処理方法を用いて、複数の受け取られた又はアクセスされた信号がフィルタ処理される。使用され得るフィルタ処理方法としては、QRS群除去、中央値フィルタリング、及び周波数(バンドパス)フィルタリングが挙げられるが、これらには限定されない。当然のことながら、信号中のノイズを減少させると共に信号品質を高めるために、他のフィルタリング方法も使用できる。演算ないしステップ806では、当該複数の信号から、或る信号が選択される。本明細書において説明するように、選択される信号は、離散長を有しても良く、又は、これよりも長い信号のエポックであっても良い。演算ないしステップ808では、信号セグメントが当該選択された信号から選択される。例えば、信号セグメントは、特定長さ(例えば、2秒、4秒等)のものであって良い。図3乃至図5を参照して説明したように、信号セグメントは、少なくとも部分的にオーバーラップし得る。
演算ないしステップ810では、オフセット範囲(例えば、100msec〜350msec)の下オフセット(lower offset)(例えば、100msec)が選択された信号セグメントに対して選択される。演算ないしステップ812では、現在のオフセットが当該下オフセットに設定される。演算ないしステップ814では、基準信号セグメントが選択された信号セグメントから生じさせられる。当該基準信号セグメントは、選択された信号セグメントの開始時から、選択された信号セグメントの終わりから現在のオフセットを差し引いた値により表される終了時まで、延びる。
演算ないしステップ816では、現在のオフセットから選択された信号セグメントの終わりまで延びる選択された信号セグメントに基づき、且つ、演算ないしステップ814で生成された基準信号セグメントに基づき、例えば図9に記載されている相関計算式を用いて、相関値(例えば、相関係数(r))が求められる。演算ないしステップ818では、現在のオフセットが時間インクリメント(例えば、1msec、2msec、又は別の時間インクリメント)だけ増分される。演算ないしステップ820では、現在のオフセットがオフセット範囲(例えば、100msec〜350msec)の上オフセット(upper offset)(例えば、350msec)に等しいかどうかの判定が行われる。
演算ないしステップ820で現在のオフセットが上オフセットに等しくないことが判定された場合、方法800は、演算ないしステップ814〜820を再び実施することによって続行する。しかしながら、演算ないしステップ820で現在のオフセットが上オフセットに等しいことが判定されると、方法800は、演算ないしステップ822で続行して、選択された信号セグメントについて最も高い相関値(例えば、相関係数(r))を選択する。
変形実施形態では、演算808〜820は、減分されるオフセットに基づいて実行されても良い。具体的に説明すると、演算ないしステップ812での現在のオフセットを、選択された信号セグメントに関して上オフセット(例えば、350msec)に設定しても良い。演算ないしステップ818で、現在のオフセットを時間デクリメント(例えば、1msec、2msec、又は他の時間インクリメント)だけ減分することができる。同様に、演算ないしステップ820で、現在のオフセットが下オフセット(例えば、100msec)に等しいかどうかの判定が行われ得る。
演算ないしステップ824では、選択された信号の全ての信号セグメントが処理されたかどうかの判定が行われる。選択された信号の全ての信号セグメントは(未だ)処理されていないと判定された場合、方法800は、演算ないしステップ808〜824を再び実施することによって続行する。しかしながら、演算ないしステップ824で選択された信号の全ての信号セグメントが処理されたと判定されると、方法800は、演算ないしステップ826で続行する。理解を容易にするために演算ないしステップ808〜824は連続した順序で説明されているが、これら演算ないしステップは、信号の複数の又は全ての信号セグメントについて、同時に又は互い違いに(staggered manner)実施され得る。
演算ないしステップ826では、全ての信号が処理されたかどうかの判定を行う。全ての信号は処理されていないと判定されると、方法800は、次に選択された信号について演算ないしステップ806〜826を実施することによって続行する。しかしながら、演算ないしステップ826で全ての信号が処理されたと判定されると、方法800は、演算ないしステップ828で続行して、処理した信号の全ての同一の(ほぼ同時期に存在する)信号セグメントについて複合相関値(例えば、和、平均値、SMR)を生成する。
方法800は、異なる選択された時間オフセット範囲(例えば、100msec〜350msec、300msec〜500msec、450msec〜1300msec、又は別のオフセット範囲)と関連して実施できる。かくして、演算ないしステップ810及び820で、下オフセット及び上オフセットは、それぞれ、選択されたオフセット範囲のうちの開始オフセット及び終了オフセットに一致するのが良い。
演算ないしステップ830では、複数の複合相関値のうちの最も高い複合相関値と関連した1組の信号セグメント(ほぼ同時期に存在する信号セグメント)が選択される。選択された信号セグメントを用いて、治療および除去のために心調律障害の発生源を識別する時間及び精度並びに心調律障害の発生源を標的にする時間及び精度を高めるのが良い。例えば、選択された信号セグメントはブリッグス等に付与された米国特許第8,165,666号明細書によって入力信号として処理されることができる。当該米国特許は、当該参照により引用され、その記載内容全体が本明細書の一部とされる。
図9は、所与の時間オフセット(例えば、100msec)についての信号セグメントの例示の相関計算式900を示している。
相関計算式900は、2つの別々のデータシリーズの積の合計として相関係数(r)を定め、各x又はyは、n個の要素を有する別々のデータシリーズの各msecでの電圧要素を表している。図示のように、相関係数は、範囲−1〜1に標準化される。具体的に説明すると、x及びyは、図6〜図8を参照して説明したように、相関されるべき信号セグメントに基づいて長さ(サイズ)の等しいデータシリーズを表している。nを信号セグメントの長さ(msec)である(例えば、長さが2,000msecの信号セグメントについてn=2,000)としてtを時間オフセットであるとすると、xは、信号セグメントの開始時から始まってn−tに至るまでの電圧要素によってポピュレートされ、yは、信号セグメントのtから始まってnに至る電圧要素によってポピュレートされる。
Figure 2016530053
それぞれ、xデータシリーズ及びyデータシリーズの平均値を表している。x(i)及びy(i)はそれぞれ、データシリーズx及びy中のi番目の要素を表している。例示の相関計算式900は、所与の時間オフセットに関して、信号セグメントの相関係数(r)を決定する。
図10は、図2の信号202〜206の各々の信号セグメント1002(例えば、2msec〜4msec、3msec〜5msec、4msec〜6msec信号セグメント)に関する最も高い相関係数(r)をまとめた例示の表1000を示している。
表1000は、異なる信号202〜206の同一(ほぼ同時期に存在する)信号セグメント1002(例えば、2msec〜4msec信号セグメント)に関する最も高い相関係数(r)の複合相関係数1004を更に示している。複合相関係数1004は、信号202〜206の同一の信号セグメント1002(例えば、2msec〜4msec信号セグメント)に関する最も高い相関係数(r)の和1006、最も高い相関係数(r)の平均値1008、又は最も高い相関係数(r)の二乗平均平方根(SMR)1010であるのが良い。SMR1010の計算について、図11を参照して詳細に説明する。
図11は、最も高い相関係数(r)のSMR1010に関する例示の計算式1100を示している。具体的に説明すると、xは、nというサイズを有する相関係数(r)のアレイを表している。xiは、xアレイのi番目の要素を表している。例示の計算式1100は、nで除算することによって標準化された各要素iの平方根の和を決定する。標準化された和は、SMRを求めるために2乗される。SMR1010は、複合相関係数1004と関連した他の信号セグメントの通常の相関係数(r)に対して、或る特定の信号セグメントの高い相関係数(r)の影響を減少させる。
図12は、図1の心臓120内の検出場所に位置決めされた例示のセンサ102と関連して得られた複合性心調律障害の例示の信号1202を示している。この例では、複合性調律障害の信号1202は、図20に示された60秒の信号のエポック(信号)を表しており、その目的は、複合性心調律障害と関連した1つの調律から別の調律までの移行(例えば、約34秒での)を明確に示すことにある。
説明を分かりやすくすると共に簡潔にするために、たった1つの信号1202が示されているが、カテーテル101によって複数の信号を提供することができる。例えば、バスケットカテーテルの対応のセンサからの64個の信号が、例えば図2を参照して説明したように、心臓120の右心房122のパノラマ図を提供することができる。本明細書において説明するように、かかる信号は、ほぼ同時期に(例えば、同時に又はほぼ同時に)得られる。
図2〜図5を参照して説明したように、複数の信号(信号1202を含む)は、コンピュータ計算装置116を用いて、セグメント長(例えば、長さが2秒)を有する複数の信号セグメントに同様にセグメント化できる。同様に、信号セグメントは、少なくとも部分的にオーバーラップした複数の信号セグメント、例えば、24〜26秒、25〜27秒、26〜28秒、……、39〜41秒及び40〜42秒の信号セグメント、を含むのが良い。信号1202の例えばその全長に基づくセグメント化のために、別のセグメント長(例えば、4秒)を選択することもできる。
この例によれば、ほぼ同時期に存在する信号セグメント(例えば、ほぼ同一の開始時期を有する信号セグメント)の各群は、複数の信号相互間で連続して処理される。信号1202を参照すると、信号1202の開始時からの信号セグメントであって信号セグメント1204(例えば、34〜36秒)を含む信号セグメントの各々が、かくして、信号セグメントの各々の複数のオフセットでの複数の相関係数(r)を求めるために、第1の選択された時間オフセット範囲にわたって、それ自体に連続的に相関させられる(例えば、図9の相関計算式を用いて)。信号1202の各信号セグメントに関し、ほぼ同時期に存在する(例えば、ほぼ同一の開始時期を有する)複数の信号の信号セグメントが、同様に(図9の相関計算式を用いて)相関させられ、かくして、一群の信号セグメント(例えば、この群は、ほぼ同一の開始時期を有する信号セグメントを含む)が形成される。
第1の選択した時間オフセット範囲は、生物学的調律障害(例えば、心調律障害)のサイクル長と関連付けられている。上述したように、第1の選択範囲は、100msecから350msecまでであるのが良い。第1の選択したオフセット範囲にわたる上述の相関は、複数の信号のほぼ同時期に存在する信号セグメントの各連続した群について繰り返される。
各群のほぼ同時期に存在する信号セグメントの各々について、最も高い相関係数(r)が選択される。さらに、複数の信号の同一の信号セグメント(例えば、複数の信号についてほぼ同一の開始時期を有する各群の信号セグメント)について、第1の複合相関値(例えば、平均相関係数)を生じさせる。第1の選択範囲(例えば、100msec〜350msec)及び信号セグメントの平均オフセットが、第1の複合相関値と関連して記憶されるのが良い。しかる後、第1の複合相関値が選択されたしきい値(例えば0.3)と比較される。
第1の複合相関係数が同時期に存在する信号セグメントの群について選択されたしきい値よりも大きい場合、信号セグメントの次に連続して位置する群が、第1の選択範囲を用いて、上述したように相関させられ、信号セグメントの当該連続して位置する群について、第1の複合相関値(例えば、平均相関係数)が生成される。
しかしながら、第1の複合相関係数が選択されたしきい値(例えば、0.3)以下である場合、300msecから500msecまでの第2の選択した範囲を用いて、当該群の信号セグメントが再相関させられる。当該群について第2の複合相関値(例えば、平均相関係数)が生成される。第2の選択範囲(例えば、300msecから500msecまで)及び信号セグメントの平均オフセットが、第2の複合相関値と関連して記憶されるのが良い。しかる後、第2の複合相関係数が選択されたしきい値と比較される。
第2の複合相関係数が選択されたしきい値を超えている場合、信号セグメントの次に連続して位置する群が、上述したように相関させられ(例えば、100msecから350msecまでの第1の選択された範囲を用いて)、当該連続して位置する次の群について、複合相関値(例えば、平均相関係数)が生成される。
しかしながら、第2の複合相関係数が選択されたしきい値以下である場合、450msecから1300msecまでの第3の選択した範囲を用いて、当該群の信号セグメントが相関させられる。当該群について第3の複合相関値(例えば、平均相関係数)が生成される。第3の選択した範囲(例えば、450msecから1300msecまで)及び信号セグメントの平均オフセットが、第3の複合相関値と関連して記憶されるのが良い。しかる後、第3の複合相関係数が選択されたしきい値と比較される。幾つかの実施形態では、第1及び/又は第2の複合相関係数が選択されたしきい値を上回っている場合であっても、第2及び第3の複合相関係数が生成される。
第3の複合相関係数が選択されたしきい値を超えている場合、信号セグメントの次に連続して位置する群が、上述したように相関させられ(例えば、100msecから350msecまでの第1の選択された範囲を用いて)、当該連続して位置する次の群について、複合相関値(例えば、平均相関係数)が生成される。
しかしながら、第3の複合相関係数が選択したしきい値以下である場合、「用いない」表示が当該群の信号セグメントに関連して記憶され、例えば、複数の信号のほぼ同時期に存在する信号セグメント相互間に制限された一貫性が存在するので、当該群は用いられるべきでない。幾つかの実施形態では、追加の又は異なるオフセット選択範囲を用いて、4つ以上の複合相関係数を計算することができる。
複数の信号のほぼ同時期に存在する信号セグメントの残りの群が、それに応じて、例えばオフセットの範囲(例えば、100msecから350msecまで、300msecから500msecまで、450msecから1300msecまで)を用いて、処理される。例示の信号1202において、最大約34秒までの信号セグメントは、図2〜図5に示されている信号セグメントとほぼ同じである。予想できることとして、第1の選択されたオフセット範囲(100msecから350msecまで)と関連した第1の複合相関値は、約34秒までの当該群の信号セグメントでは選択されたしきい値よりも大きく、約34秒又はそれ以後の信号セグメントでは選択されたしきい値よりも小さいであろう。異なるオフセット範囲によるこれら信号セグメントの処理について、図13〜図18を参照して以下に詳細に説明する。
図13は、第1の選択されたオフセット範囲に従った、図12に示された信号1202の信号セグメント1204の例示の時間シフト(又はオフセット)1302を示す図である。
上述したように、第1の選択されたオフセット範囲(例えば、100msecから350msecまで)に従って、250個の時間オフセットが存在する。図示していないが、図9の例示の相関計算式を用いて、100msecから350msecまでの範囲の時間オフセットの各々について相関係数(r)が計算される。
信号セグメント1204に関し、122msecのオフセットで、最大相関係数(r)=0.11が得られる。異なる相関方法又は技術を用いて例示の信号セグメント1204をそれ自体と相関(例えば、自己相関)させてもよい。信号セグメント1204に関する最大相関係数は、図12を参照して説明したように、しきい値(例えば、0.3)未満である。
図14は、図12に示された信号1202の信号セグメント1204の第1の選択されたオフセット範囲1406に従った時間オフセットについて求められた相関係数(r)の曲線1402を示す例示のグラフ1400である。
図14に示されているように、当該相関係数は、第1の選択された時間オフセット範囲1406(例えば、100msec〜350msec)において−0.10〜0.11であり、最大相関係数(r)=0.11が122msecの時間オフセットで生じている。
同様な曲線を生成して、他の信号のほぼ同時期に存在する信号セグメントについて、第1の選択されたオフセット範囲内の時間オフセットについて求められた相関係数(r)を示すことができる。
図15は、第2の選択されたオフセット範囲に従った、図12に示されている信号1202信号セグメント1204の例示の時間シフト(又はオフセット)1502を示す図である。
第2の選択されたオフセット範囲(例えば、300msec〜500msec)に従って、200個の時間オフセットが存在する。図示していないが、図9の例示の相関計算式を用いて、300msecから500msecまでの範囲の時間オフセットの各々について相関係数(r)が計算される。
信号セグメント1204に関し、406msecのオフセットで、最大相関係数(r)=0.05が得られる。異なる相関方法又は技術を用いて例示の信号セグメント1204をそれ自体と相関(例えば、自己相関)させてもよい。信号セグメント1204に関する最大相関係数は、図12を参照して説明したように、しきい値(例えば、0.3)未満である。
図16は、図12に示された信号1202の信号セグメント1204の第2の選択されたオフセット範囲1606に従った時間オフセットについて求められた相関係数(r)の曲線1602を示す例示のグラフ1600である。
図16に示されているように、当該相関係数は、第2の選択された時間オフセット範囲1606(例えば、300msec〜500msec)において−0.08〜0.05であり、最大相関係数(r)=0.05が406msecの時間オフセットで生じている。
同様な曲線を生成して、他の信号のほぼ同時期に存在する信号セグメントについて、第2の選択されたオフセット範囲内の時間オフセットについて求められた相関係数(r)を示すことができる。
図17は、第3の選択されたオフセット範囲に従った、図12に示されている信号1202の信号セグメント1204の例示の時間シフト(又はオフセット)1702を示す図である。
第3の選択されたオフセット範囲(例えば、450msec〜1300msec)に従って、750個の時間オフセットが存在する。図示していないが、図9の例示の相関計算式を用いて、450msecから1300msecまでの範囲の時間オフセットの各々について相関係数(r)が計算される。
信号セグメント1204に関し、698msecのオフセットで、最大相関係数(r)=0.97が得られる。異なる相関方法又は技術を用いて例示の信号セグメント1204をそれ自体と相関(例えば、自己相関)させてもよい。信号セグメント1204に関する最大相関係数は、図12を参照して説明したように、しきい値(例えば、0.3)を超える。
図18は、図12に示された信号1202の信号セグメント1204の第3の選択されたオフセット範囲1806に従った時間オフセットについて求められた相関係数(r)の曲線1802を示す例示のグラフ1800である。
図18に示されているように、当該相関係数は、第3の選択された時間オフセット範囲1806(例えば、450msec〜1300msec)において−0.45〜0.97であり、最大相関係数(r)=0.97が698msecの時間オフセットで生じている。
同様な曲線を生成して、他の信号のほぼ同時期に存在する信号セグメントについて、第3の選択されたオフセット範囲内の時間オフセットについて求められた相関係数(r)を示すことができる。
図19は、第2の実施形態に従って、生物学的調律障害(例えば心調律障害)の分析のために信号セグメントを選択する例示の方法1900を示す流れ図である。例示の方法1900は、図1に示されているコンピュータ計算装置116によって実施され得る。
具体的に説明すると、例示の方法1900は、演算ないしステップ1902で始まり、この演算ないしステップ1902では、複数の信号が心臓120内に設けられたセンサから信号処理装置114を介してコンピュータ計算装置116によって受け取られる、又は、アクセスされる。かかる信号は、特定長さのものであるのが良く、例えば、60秒以上の長さのものである。幾つかの実施形態又は幾つかの観点では、センサからの信号の少なくとも一部分が信号処理装置114によって記録され得て、次にコンピュータ計算装置116に提供され得る。
演算ないしステップ1904では、1つ又は2つ以上のフィルタ処理方法を用いて、複数の受け取られた又はアクセスされた信号がフィルタ処理される。使用され得るフィルタ処理方法としては、QRS群除去、中央値フィルタリング、及び周波数(バンドパス)フィルタリングが挙げられるが、これらには限定されない。当然のことながら、信号中のノイズを減少させると共に信号品質を高めるために、他のフィルタリング方法も使用できる。
演算ないしステップ1906では、複数の信号の信号セグメントに関する時間の長さが設定される。本明細書において説明するように、当該時間の長さは、2秒間、4秒間、又は別の時間の長さ、であって良い。演算ないしステップ1908では、複数の信号のほぼ同時期に存在する信号セグメント(例えば、ほぼ同一の開始時期で始まる信号セグメント)を選択するために、信号セグメントに関する開始時期が設定される。当初、当該開始時期は、信号の開始時に設定されるのが良い(例えば、0の開始時期)。しかる後、複数の信号の連続した信号セグメントを選択するために、開始時期が演算ないしステップ1908〜1940の各繰り返しのために時間インクリメント(例えば、1秒)だけ増分されるのが良い。
演算ないしステップ1910では、心調律障害のサイクル長と関連付けられたオフセット範囲が選択される。図12〜図18を参照して説明したように、第1の範囲(例えば、100msecから350msecまで)が当初選択される。以下に説明するように、選択することができる他の範囲としては、第2の範囲(例えば、300msecから500msecまで)及び第3の範囲(例えば、450msecから1300msecまで)が挙げられる。他の範囲も規定することができる。
演算ないしステップ1912では、複数の信号から或る信号が選択される。本明細書において説明するように、選択される信号は、離散長を有しても良く、又は、これよりも長い信号のエポックであっても良い。演算ないしステップ1914では、信号セグメントが、当該選択された信号から選択される。選択された信号セグメントは、開始時期と時間の長さとによって定められる。
演算ないしステップ1916では、現在のオフセットが演算ないしステップ1910で選択されたオフセット範囲の下オフセットに設定される。演算ないしステップ1918では、基準信号セグメントが、選択された信号セグメントから生成される。当該基準信号セグメントは、選択された信号セグメントの開始時から、選択された信号セグメントの終わりから現在のオフセットを差し引いた値により表される終了時まで、延びる。
演算ないしステップ1920では、現在のオフセットから選択された信号セグメントの終わりまで延びる選択された信号セグメントに基づき、且つ、演算ないしステップ918で生成された基準信号セグメントに基づき、例えば図9に記載されている相関計算式を用いて、相関値(例えば、相関係数(r))が求められる。演算ないしステップ1922では、現在のオフセットが時間インクリメント(例えば、1msec、2msec、又は別の時間インクリメント)だけ増分される。演算ないしステップ1924では、現在のオフセットが選択オフセット範囲の上オフセット(例えば、第1の選択範囲について350msec)に等しいかどうかの判定が行われる。
演算ないしステップ1924で現在のオフセットが上オフセットに等しくないことが判定された場合、方法1900は、演算ないしステップ1918〜1924を再び実施することによって続行する。しかしながら、演算ないしステップ1924で現在のオフセットが上オフセットに等しいことが判定されると、方法1900は、演算ないしステップ1926で続行して、選択された信号セグメントについて最も高い相関値(例えば、相関係数(r))を選択する。
変形実施形態では、演算1914〜1924は、減分されるオフセットに基づいて実行されても良い。具体的に説明すると、演算ないしステップ1916での現在のオフセットを、選択された信号セグメントに関して上オフセットに設定しても良い。演算ないしステップ1922で、現在のオフセットを時間デクリメント(例えば、1msec、2msec、又は他の時間インクリメント)だけ減分することができる。同様に、演算ないしステップ1924で、現在のオフセットが選択範囲の下オフセット(例えば、第1の選択範囲について100msec)に等しいかどうかの判定が行われ得る。
演算ないしステップ1928では、全ての信号のほぼ同時期に存在する信号セグメントが処理された(例えば、選択されたセグメントの開始時期を有する信号セグメントについて全ての信号が処理された)かどうかの判定が行われる。全ての信号は(未だ)処理されていないと判定された場合、方法1900は、演算ないしステップ1912〜1928を再び実施することによって続行する。しかしながら、演算ないしステップ1928で全ての信号が処理されたと判定されると、方法1900は、演算ないしステップ1930で続行する。理解を容易にするために演算ないしステップ1912〜1928は連続した順序で説明されているが、これら演算ないしステップは、信号のほぼ同時期に存在する信号セグメントの異なる群について、同時に又は互い違いに実施され得る。
演算ないしステップ1930では、複合相関値(例えば、和、平均値、SMR)が、信号の全てに亘って同一の(ほぼ同時期に存在する)信号セグメントについて生成される。演算ないしステップ1932で、選択された範囲及び平均オフセット(例えば、ほぼ同時期に存在する信号セグメントのオフセットの平均値)が、複合相関値と関連して記憶される。
演算ないしステップ1934では、複合相関値がしきい値(例えば、0.3)を超えているかどうかについての判定が行われる。演算ないしステップ1934において複合相関値がしきい値以下であると判定された場合、方法1900は、演算ないしステップ1936で続行し、当該演算ないしステップ1936では、別のオフセット範囲が選択可能であるかどうかの判定が行われる。もし可能であれば、方法1900は、演算ないしステップ1910で続行し、当該演算ないしステップ1910で、別のオフセット範囲が選択され、演算ないしステップ1910〜1934が繰り返される。選択されるオフセット範囲は、第2のオフセット範囲(例えば、300msecから500msecまで)であるのが良い。上述したように、3つの規定されたオフセット範囲(例えば、100msecから350msecまで、300msecから500msecまで、450msecから1300msecまで)が、演算ないしステップ1936で選択され得る。演算ないしステップ1936で、他の範囲も規定され、かくして選択されることができる。別のオフセット範囲が選択可能ではない場合、信号セグメントは演算ないしステップ1938で使用されないものとしてマーク付けされ、方法1900は、演算ないしステップ1940で続行する。
演算ないしステップ1934において、複合相関値がしきい値を超えていると判定された場合、方法1900は、演算ないしステップ1940で続行し、当該演算ないしステップ1940において、信号中の全ての信号セグメントが処理されたかどうかについての判定が行われる。具体的に説明すると、複数の信号の信号セグメントに関する全ての開始時期が処理されたかどうかが判定される。もしそうでなければ、方法1900は、複数の信号中の追加の信号セグメントについて演算ないしステップ1908〜1940を実施し、即ち、複数の信号中でまだ処理されていないほぼ同時期に存在する信号セグメントに関する次の連続した開始時期を設定する。複数の信号中の全ての信号セグメントが処理されたことを演算ないしステップ1940で判定した後は、方法1900は、演算ないしステップ1942で続行する。演算ないしステップ1942では、選択されたオフセット範囲についての最も高い複合相関値と関連した1組のほぼ同時期に存在するセグメントが選択される。方法1900は、演算ないしステップ1944で終わる。
選択された信号セグメントを用いると、心調律障害の発生源を突き止める時間及び精度、並びに、治療及び除去のために心調律障害の発生源を標的にする際の時間及び精度を向上させることができる。例えば、選択された信号セグメントをブリッグス等(Briggs et al. )に付与された米国特許第8,165,666号の発明によって入力信号として処理することができる。当該米国特許は当該参照により引用され、その記載内容全体を本明細書の一部とする。
図20は、図1の心臓120内の検出場所に位置決めされた例示のセンサ102と関連して得られた複合性心調律障害の例示の信号2002を示す図である。信号2002は、60秒信号であり、図12に示されているエポック又は信号1202を含む。図示のように、区分2004の調律(リズム)が、区分2006の別の調律に移行している(例えば、約34秒で)。
図21は、信号2002を含む複数の信号のセグメントと関連した複合相関係数の例示の曲線2102を示している。複合相関係数は、信号2002を含む複数の信号の1秒ごとに得られる2秒長さ、例えば0〜2秒、1〜3秒、……、32〜34秒、……、57〜59秒、及び58〜60秒、のオーバーラップした信号セグメントと関連している。
区分2004の間、複数の信号のほぼ同時期に存在する信号セグメントの複合相関係数は、約0.58と約0.76との間で変化する。複数の信号のほぼ2106のところでは(例えば、信号セグメント32〜34秒)、複合相関係数は、しきい値(例えば、0.3)と関連したしきい値2104を下回る。区分2006の間、複数の信号中のほぼ同時期に存在する信号セグメントの複合相関係数は、約0.8と約0.96との間で変化する。
図示のように、区分2004の間の複数の信号の周期性の一貫性、2106のところでの周期性の一貫性の著しい低下、及び、区分2006の間での複数の信号の周期性の一貫性の再開、がみとめられる。区分2004の間における周期性の一貫性の方が、区分2006の間における周期性の一貫性よりも低い。心調律障害と関連した一貫性の低い調律は、一貫性の高い調律へと移行している。この例では、一貫性の高い調律は、正常な洞調律である。具体的に説明すると、心調律障害と関連した調律の停止、及び、正常な洞調律への変換、がみとめられる。
図22は、信号2002を含む複数の信号の同時期に存在するセグメントを考慮した、心調律障害のサイクル長(CL)と関連した複合オフセットの例示の曲線2202を示す図である。
区分2004の間、心調律障害のCLと関連した複合オフセットは、約220msecであり、これは、第1の選択可能なオフセット範囲(例えば、100msecから350msecまで)内にある。区分2006の間、正常な洞調律と関連した複合オフセットは、約700msecであり、これは、第3の選択可能なオフセット範囲(例えば、450msecから1300msecまで)内にある。
第2の実施形態によれば、選択可能なオフセット範囲の使用により、心調律障害と関連した1つ又は2つ以上の調律、並びに、調律相互間での移行、の検出が可能である。移行は、異常な調律相互間の移行(例えば、AFからATまで、その他)だけでなく、異常な調律(例えば、AF、AT、その他)から正常な洞調律への移行を含む。これにより、更なる分析のために心調律障害と関連した特定の調律と関係するほぼ同時期に存在する信号セグメント(複数の信号のほぼ同時期に存在する信号セグメント)の選択(例えば、オフセット範囲を用いて)が可能である。
図23は、汎用コンピュータシステム1200の例示の実施形態のブロック図である。コンピュータシステム2300は、図1の信号処理装置114及びコンピュータ計算装置116であり得る。コンピュータシステム1200は、本明細書において開示する方法又はコンピュータ利用機能のうちの任意の1つ又は2つ以上をコンピュータシステム2300に実施させるべく実行され得る1組の命令を含み得る。コンピュータシステム2300又はその任意の一部分は、スタンドアローン型装置として動作しても良く、或いは、例えばネットワーク又は他の接続方式を用いて他のコンピュータシステム又は周辺機器に接続されても良い。例えば、コンピュータシステム1200は、信号処理装置114に作動的に連結されても良い。
コンピュータシステム2300は、種々の装置、例えばパーソナルコンピュータ(PC)、タブレット型PC、携帯情報端末(PDA)、モバイル装置、パームトップ型コンピュータ、ラップトップ型コンピュータ、デスクトップ型コンピュータ、通信機器、制御システム、ウェブアプライアンス、又は、当該機械によってとられる(なされる)べき行為を特定する1組の命令を(逐次的に又は他の態様で)実行できる任意の他の機械、として実装されても良く、或いはこれらに組み込まれても良い。さらに、単一コンピュータシステム2300が図示されているが、「システム」という用語は、1つ又は2つ以上のコンピュータ機能を実行するための1組の又は複数組の命令を個々に又は一括に実行するシステム又はサブシステムの任意の集まりを含むものと解されるべきである。
図23に示されているように、コンピュータシステム2300は、プロセッサ2302、例えば中央処理装置(CPU)、グラフィックス処理ユニット(GPU)、又はこれら両方、を含むのが良い。さらに、コンピュータシステム2300は、バス2326を介して互いに通信可能なメインメモリ2304とスタティックメモリ2306を含むのが良い。図示のように、コンピュータシステム2300は、ビデオディスプレイユニット2310、例えば液晶ディスプレイ(LCD)、有機発光ダイオード(OLED)、フラットパネルディスプレイ、ソリッドステイトディスプレイ、又は陰極線管(CRT)を更に含むのが良い。加うるに、コンピュータシステム2300は、入力装置2312、例えばキーボード、及びカーソル制御装置2314、例えばマウス、を含むのが良い。コンピュータシステム2300は、ディスクドライブユニット2316、信号発生装置2322、例えばスピーカ又はリモートコントロール、及びネットワークインターフェース装置2308を更に含むのが良い。
一形態又は一観点では、図23に示されているように、ディスクドライブユニット2316は、コンピュータ可読媒体2318を含むのが良く、当該コンピュータ可読媒体には、1つ又は2つ以上の組の命令2320、例えばソフトウェア、が入れられているのが良い。さらに、命令2320は、本明細書で説明する方法又は論理のうちの1つ又は2つ以上を具体化するものであるのが良い。特定の実施形態又は観点では、命令2320は、コンピュータシステム1200による実行中、メインメモリ2304、スタティックメモリ2306及び/またはプロセッサ2302内に完全に又は少なくとも部分的に存在するのが良い。メインメモリ2304及びプロセッサ2302は、コンピュータ可読媒体を更に含むのが良い。
変形実施形態又は別の観点では、専用ハードウェア具体化例、例えば特定用途向け集積回路、プログラマブルロジックアレイ及び他のハードウェア装置、が本明細書において説明する方法のうちの1つ又は2つ以上を具体化するよう構成され得る。種々の実施形態又は種々の観点の装置及びシステムを含み得るアプリケーションは、広義には、種々の電子及びコンピュータシステムを含むことができる。本明細書において説明される1つ又は2つ以上の実施形態又は観点は、モジュール間で且つモジュールを介して通信できる関連制御及びデータ信号を含む2つ又は3つ以上の特定の相互接続ハードウェアモジュール又は装置を用いて、又は特定用途向け集積回路の幾つかの部分として、機能を実行することができる。したがって、本発明のシステムは、ソフトウェア、ファームウェア、及びハードウェア具体化例を含む。
種々の実施形態又は観点によれば、本明細書において説明した方法は、プロセッサ可読媒体内に有形的に具体化されたソフトウェアプログラムによって実施可能であると共に、プロセッサによって実行可能である。さらに、例示の非限定的な実施形態又は観点では、具体化例としては、分散処理、コンポーネント/オブジェクト分散処理、及びパラレル処理が挙げられる。変形例として、バーチャルコンピュータシステム処理方式が、本明細書に記載されている方法又は機能性の1つ以上を実現するために構成されても良い。
コンピュータ可読媒体が命令2320を含むか、又は伝達される信号に応答して命令2320を受信して実行する、ということも想定され、その結果、ネットワーク2324に接続されているデバイスは当該ネットワーク2324を介して音声、ビデオ又はデータを通信することができる。さらに、命令2320は、ネットワークインターフェースデバイス2308を介してネットワーク2324で送信又は受信することができる。
コンピュータ可読媒体は、単一の媒体として示されているが、「コンピュータ可読媒体」という用語は、単一の媒体又は複数の媒体、例えば、集中型ないし分散型データベース、及び/又は、1組又は2組以上の命令を記憶する関連キャッシュ及びサーバ、を含む。「コンピュータ可読媒体」という用語は又、プロセッサによって実行するための1組の命令を記憶するか、コード化するか、運ぶことができる任意の媒体、又は、コンピュータシステムに本明細書において開示される方法又は動作の任意の1つ以上を実行させる任意の媒体、を含むものとする。
特定の非限定的な例としての実施形態又は観点において、コンピュータ可読媒体は、固体メモリ(例えば、メモリカード又は、1つ以上の不揮発性読取り専用メモリを収容する他のパッケージ)を含む場合がある。さらに、コンピュータ可読媒体は、ランダムアクセスメモリ又は他の揮発性再書き込み可能メモリであっても良い。加うるに、コンピュータ可読媒体は、光磁気又は光学的媒体(例えば、伝送媒体を通じて伝達される信号のような搬送波信号を捕えるディスク又はテープ又はその他の記憶装置)を含む場合がある。電子メールに対するデジタル添付ファイル又は他の自己内蔵型情報アーカイブ又はアーカイブのセットは、有形の記憶媒体に等しい配布媒体と考えられる場合がある。したがって、コンピュータ可読媒体又は配布媒体及び他の等価物及び後継媒体(データ又は命令が格納され得る)の任意の1つ以上が、本発明の範囲に含まれる。
種々の実施形態又は観点によれば、本明細書に記載されている方法は、コンピュータプロセッサ上で動く1つ以上のソフトウェアプログラムとして実現可能である。特定用途向け集積回路、プログラマブルロジックアレイ、及び他のハードウェアデバイスを含むがこれらに限定されない専用ハードウェアの具体化も、同様に、本明細書に記載されている方法を実施するために構成可能である。さらに、分散処理又はコンポーネント/オブジェクト分散処理、パラレル処理、又はバーチャルマシン処理を含むがこれらに限定されない他のソフトウェア具体化例も、本明細書に記載されている方法を実施するために構成可能である。
開示された方法を実施するソフトウェアが、オプションとして、有形の記憶媒体、例えば、磁気媒体、例えばディスク又はテープ、光磁気又は光学的媒体、例えばディスク、固体媒体、例えばメモリカード、又は、1つ以上の読取り専用(不揮発性)メモリ、ランダムアクセスメモリ、ないし他の再書き込み可能な(揮発性)メモリを収納する他のパッケージ、上に記憶されて良い。電子メールに対するデジタル添付ファイル又は他の自己内蔵型情報アーカイブ又はアーカイブのセットは、有形の記憶媒体に等しい配布媒体と考えられる場合がある。したがって、本明細書に列記される有形の記憶媒体又は配布媒体、及び他の等価物及び後継媒体(ソフトウェア実装物が格納され得る)は、本発明の範囲に含まれる。
かくして、心臓賦活情報を再構成するシステム及び方法が説明された。特定の例示的実施形態又は観点が説明されたが、本発明の広い範囲から逸脱することなく、これらの実施形態又は観点に対して種々の改造及び変更がなされ得ることは明らかであろう。したがって、明細書及び図面は、限定的意味ではなく例示的であるとみなされるべきである。本発明の一部をなす添付の図面は、本発明の内容が実施され得る特定の実施形態又は観点を、限定ではなく説明のために示している。図示された実施形態又は観点は、当業者が本明細書において開示される教示を実施することができる程度に十分詳細に記載されている。そこから、他の実施形態又は観点を利用することができて導き出すことができ、従って、構造的及び論理的な置換及び変更が、本発明の範囲から逸脱することなく、なされ得る。したがって、この詳細な説明は、本発明を限定するものと解されてはならず、種々の実施形態又は観点の範囲は、当該請求項の記載に基づいて与えられる全均等範囲と共に、添付の特許請求の範囲の請求項の記載のみに基づいて定められる。
本発明のかかる実施形態又は観点は、本明細書において、事実上2つ以上の発明が開示されている場合でも、単なる説明の便宜上、本願の要旨を任意の単一の発明又は発明概念に自発的に限定するという意図をもたないで、「発明」という用語によって個々に及び/またはまとめて言及されている場合がある。かくして、特定の実施形態又は観点が本明細書において例示されると共に説明されたが、同一の目的を達成するために計画される任意の構成が、説明された特定の実施形態又は観点に置換され得る、ということが理解されるべきである。本発明は、種々の実施形態又は観点の任意且つ全ての改造例又は変形例を含むことが意図されている。上述の実施形態又は観点の組み合わせ、及び、本明細書に具体的に説明しなかった他の実施形態又は観点は、上記説明を読めば当業者には明らかであろう。
要約書は、37 C.F.R §1.72 (b)に準拠して提供されており、かかる要約書は、読者が技術的開示内容の本質及び要点を迅速に確認することを許容するであろう。要約書は、これが特許請求の範囲に記載された本発明の範囲又は意味を解釈し又は限定するためには使用されないという理解のもとで提出されている。
実施形態又は観点の上記説明において、種々の特徴が、本発明の説明を合理化するために、単一の実施形態の状態で集約されている。このような開示の仕方は、クレーム請求された実施形態又は観点が各請求項で明示的に記載されているよりも多くの特徴を有するということを反映するものと解されるべきではない。むしろ、以下の特許請求の範囲の記載が反映しているように、本発明の内容は、単一の開示された実施形態又は観点の全ての特徴よりも少ない状態にある。かくして、以下の特許請求の範囲の請求項の内容は、詳細な説明中に組み込まれ、各請求項の内容は、別個の例示の実施形態又は観点としてそれ自体に基づく。本明細書において説明した種々の実施形態又は観点は、詳細な説明に明示的には記載されていない異なる組み合わせ状態に組み合わせることができ又は寄せ集めることができる。さらに、かかる種々の組み合わせを包含する請求項の内容は、同様に、詳細な説明中に組み込み可能な別個の例示の実施形態または観点としてそれら自体に基づく。

Claims (32)

  1. 複数の心臓信号のセグメントを処理する方法であって、
    少なくとも1つの第1の心臓信号中の信号セグメントを、複数の時間オフセットによってシフトされた当該信号セグメントと相関させて、当該信号セグメントに関する最も高い相関係数を決定する工程と、
    少なくとも1つの第2の心臓信号中の信号セグメントを、複数の時間オフセットによってシフトされた当該信号セグメントと相関させて、当該信号セグメントに関する最も高い相関係数を決定する工程と、
    前記少なくとも1つの第1の心臓信号中及び前記少なくとも1つの第2の心臓信号中において追加の信号セグメントに対して前記相関を繰り返して、当該追加の信号セグメントの各々に関する最も高い相関係数を決定する工程と、
    前記少なくとも1つの第1の心臓信号の前記信号セグメント及び前記少なくとも1つの第2の心臓信号の前記信号セグメントに関する最も高い相関係数を用いて、複数の複合相関係数を生じさせる工程と、
    前記少なくとも1つの第1の心臓信号からの少なくとも1つの信号セグメント及び前記少なくとも1つの第2の心臓信号からの少なくとも1つの信号セグメントを含む1組の信号セグメントを選択する工程と、
    を備え、
    前記少なくとも1つの第1の心臓信号の前記信号セグメントは、前記少なくとも1つの第2の心臓信号の前記信号セグメントとほぼ同時期に存在し、
    前記1組の信号セグメントは、前記複数の複合相関係数からの最も高い複合相関係数と関連している、方法。
  2. 前記少なくとも1つの第1の心臓信号及び前記少なくとも1つの第2の心臓信号をフィルタ処理する工程を更に含む、請求項1記載の方法。
  3. 前記少なくとも1つの第1の心臓信号中の前記複数の信号セグメントは、少なくとも部分的にオーバーラップしている、請求項1記載の方法。
  4. 前記少なくとも1つの第2の心臓信号中の前記複数の信号セグメントは、少なくとも部分的にオーバーラップしている、請求項1記載の方法。
  5. 前記少なくとも1つの第1の心臓信号中の前記複数の信号セグメント及び前記少なくとも1つの第2の心臓信号中の前記複数の信号セグメントの各信号セグメントは、所定のセグメント長さのものである、請求項1記載の方法。
  6. 前記所定のセグメント長さは、長さが2秒及び4秒のうちの一方である、請求項5記載の方法。
  7. 前記複数の時間オフセットは、100msecから350msecまでの範囲にある、請求項1記載の方法。
  8. 第1の心臓信号中の第1の信号セグメントを複数の時間オフセットで相関させる前記工程は、
    前記第1の信号セグメントから第1の基準セグメントを生じさせる工程を含み、前記第1の基準セグメントは、前記第1の信号セグメントの開始時から、前記第1の信号セグメントの終わりから第1の時間オフセットを差し引いた値によって表される終了時まで、延びており、
    前記第1の信号セグメント中の前記第1の時間オフセットから前記第1の信号セグメントの前記終わりまで延びる前記第1の信号セグメントの第1の相関係数を、前記第1の基準セグメントに対して決定する工程を含み、
    前記第1の信号セグメントから第2の基準セグメントを生じさせる工程を含み、前記第2の基準セグメントは、前記第1の信号セグメントの開始時から、前記第1の信号セグメントの終わりから第2の時間オフセットを差し引いた値によって表される終了時まで、延びており、
    前記第1の信号セグメント中の前記第2の時間オフセットから前記第1の信号セグメントの前記終わりまで延びる前記第1の信号セグメントの第2の相関係数を、前記第2の基準セグメントに対して決定する工程を含み、
    前記第1の相関係数及び前記第2の相関係数の中から前記第1の心臓信号に関する前記最も高い相関係数を選択する工程を含む、請求項1記載の方法。
  9. 第2の心臓信号中の第1の信号セグメントを複数の時間オフセットで相関させる前記工程は、
    前記第1の信号セグメントから第1の基準セグメントを生じさせる工程を含み、前記第1の基準セグメントは、前記第1の信号セグメントの開始時から、前記第1の信号セグメントの終わりから第1の時間オフセットを差し引いた値によって表される終了時まで、延びており、
    前記第1の信号セグメント中の前記第1の時間オフセットから前記第1の信号セグメントの前記終わりまで延びる前記第1の信号セグメントの第1の相関係数を、前記第1の基準セグメントに対して決定する工程を含み、
    前記第1の信号セグメントから第2の基準セグメントを生じさせる工程を含み、前記第2の基準セグメントは、前記第1の信号セグメントの開始時から、前記第1の信号セグメントの終わりから第2の時間オフセットを差し引いた値によって表される終了時まで、延びており、
    前記第1の信号セグメント中の前記第2の時間オフセットから前記第1の信号セグメントの前記終わりまで延びる前記第1の信号セグメントの第2の相関係数を、前記第2の基準セグメントに対して決定する工程を含み、
    前記第1の相関係数及び前記第2の相関係数の中から前記第2の心臓信号に関する前記最も高い相関係数を選択する工程を含む、請求項8記載の方法。
  10. 複数の複合相関係数を生じさせる前記工程は、前記第1の心臓信号に関する前記最も高い相関係数と前記第2の心臓信号に関する前記最も高い相関係数との和を計算する工程を含む、請求項9記載の方法。
  11. 複数の複合相関係数を生じさせる前記工程は、前記第1の心臓信号に関する前記最も高い相関係数と前記第2の心臓信号に関する前記最も高い相関係数との平均値を計算する工程を含む、請求項9記載の方法。
  12. 複数の複合相関係数を生じさせる前記工程は、前記第1の心臓信号に関する前記最も高い相関係数と前記第2の心臓信号に関する前記最も高い相関係数との二乗平均平方根を計算する工程を含む、請求項9記載の方法。
  13. 前記複数の複合相関係数は、前記少なくとも1つの第1の心臓信号の前記信号セグメント及び前記少なくとも1つの第2の心臓信号の前記信号セグメントに関する前記最も高い相関係数の和、平均値、及び二乗平均平方根のうちの1つを用いて生成される、請求項1記載の方法。
  14. 前記方法は、前記1組の信号セグメントを処理して心調律障害の発生源を突き止める工程を更に含む、請求項1記載の方法。
  15. 心臓信号のセグメントを処理するシステムであって、
    処理装置と、
    前記処理装置によって実行されると、前記処理装置が演算を実行する、という命令を記憶する記憶装置と、
    を備え、
    前記演算は、
    少なくとも1つの第1の心臓信号中の信号セグメントを、複数の時間オフセットによってシフトされた当該信号セグメントと相関させて、当該信号セグメントに関する最も高い相関係数を決定する工程と、
    少なくとも1つの第2の心臓信号中の信号セグメントを、複数の時間オフセットによってシフトされた当該信号セグメントと相関させて、当該信号セグメントに関する最も高い相関係数を決定する工程と、
    前記少なくとも1つの第1の心臓信号中及び前記少なくとも1つの第2の心臓信号中において追加の信号セグメントに対して前記相関を繰り返して、当該追加の信号セグメントの各々に関する最も高い相関係数を決定する工程と、
    前記少なくとも1つの第1の心臓信号の前記信号セグメント及び前記少なくとも1つの第2の心臓信号の前記信号セグメントに関する最も高い相関係数を用いて、複数の複合相関係数を生じさせる工程と、
    前記少なくとも1つの第1の心臓信号からの少なくとも1つの信号セグメント及び前記少なくとも1つの第2の心臓信号からの少なくとも1つの信号セグメントを含む1組の信号セグメントを選択する工程と、
    を備え、
    前記少なくとも1つの第1の心臓信号の前記信号セグメントは、前記少なくとも1つの第2の心臓信号の前記信号セグメントとほぼ同時期に存在し、
    前記1組の信号セグメントは、前記複数の複合相関係数からの最も高い複合相関係数と関連している、システム。
  16. 前記演算は、前記少なくとも1つの第1の心臓信号及び前記少なくとも1つの第2の心臓信号をフィルタ処理する工程を更に含む、請求項15記載のシステム。
  17. 前記少なくとも1つの第1の心臓信号中の前記複数の信号セグメントは、少なくとも部分的にオーバーラップしている、請求項15記載のシステム。
  18. 前記少なくとも1つの第2の心臓信号中の前記複数の信号セグメントは、少なくとも部分的にオーバーラップしている、請求項15記載のシステム。
  19. 前記少なくとも1つの第1の心臓信号中の前記複数の信号セグメント及び前記少なくとも1つの第2の心臓信号中の前記複数の信号セグメントの各信号セグメントは、所定のセグメント長さのものである、請求項15記載のシステム。
  20. 前記所定のセグメント長さは、長さが2秒及び4秒のうちの一方である、請求項19記載のシステム。
  21. 前記複数の時間オフセットは、100msecから350msecまでの範囲にある、請求項15記載のシステム。
  22. 複数の時間オフセットでの第1の心臓信号中の第1の信号セグメントの相関を実行する前記演算は、
    前記第1の信号セグメントから第1の基準セグメントを生じさせる演算を含み、前記第1の基準セグメントは、前記第1の信号セグメントの開始時から、前記第1の信号セグメントの終わりから第1の時間オフセットを差し引いた値によって表される終了時まで、延びており、
    前記第1の信号セグメント中の前記第1の時間オフセットから前記第1の信号セグメントの前記終わりまで延びる前記第1の信号セグメントの第1の相関係数を、前記第1の基準セグメントに対して決定する演算を含み、
    前記第1の信号セグメントから第2の基準セグメントを生じさせる演算を含み、前記第2の基準セグメントは、前記第1の信号セグメントの開始時から、前記第1の信号セグメントの終わりから第2の時間オフセットを差し引いた値によって表される終了時まで、延びており、
    前記第1の信号セグメント中の前記第2の時間オフセットから前記第1の信号セグメントの前記終わりまで延びる前記第1の信号セグメントの第2の相関係数を、前記第2の基準セグメントに対して決定する演算を含み、
    前記第1の相関係数及び前記第2の相関係数の中から前記第1の心臓信号に関する前記最も高い相関係数を選択する演算を含む、請求項15記載のシステム。
  23. 複数の時間オフセットでの第2の心臓信号中の第1の信号セグメントの相関を実行する前記演算は、
    前記第1の信号セグメントから第1の基準セグメントを生じさせる演算を含み、前記第1の基準セグメントは、前記第1の信号セグメントの開始時から、前記第1の信号セグメントの終わりから第1の時間オフセットを差し引いた値によって表される終了時まで、延びており、
    前記第1の信号セグメント中の前記第1の時間オフセットから前記第1の信号セグメントの前記終わりまで延びる前記第1の信号セグメントの第1の相関係数を、前記第1の基準セグメントに対して決定する演算を含み、
    前記第1の信号セグメントから第2の基準セグメントを生じさせる演算を含み、前記第2の基準セグメントは、前記第1の信号セグメントの開始時から、前記第1の信号セグメントの終わりから第2の時間オフセットを差し引いた値によって表される終了時まで、延びており、
    前記第1の信号セグメント中の前記第2の時間オフセットから前記第1の信号セグメントの前記終わりまで延びる前記第1の信号セグメントの第2の相関係数を、前記第2の基準セグメントに対して決定する演算を含み、
    前記第1の相関係数及び前記第2の相関係数の中から前記第2の心臓信号に関する前記最も高い相関係数を選択する演算を含む、請求項22記載のシステム。
  24. 複数の複合相関係数の生成を実行する前記演算は、前記第1の心臓信号に関する前記最も高い相関係数と前記第2の心臓信号に関する前記最も高い相関係数との和を計算する演算を含む、請求項23記載のシステム。
  25. 複数の複合相関係数の生成を実行する前記演算は、前記第1の心臓信号に関する前記最も高い相関係数と前記第2の心臓信号に関する前記最も高い相関係数との平均値を計算する演算を含む、請求項23記載のシステム。
  26. 複数の複合相関係数の生成を実行する前記演算は、前記第1の心臓信号に関する前記最も高い相関係数と前記第2の心臓信号に関する前記最も高い相関係数との二乗平均平方根を計算する演算を含む、請求項23記載のシステム。
  27. 前記複数の複合相関係数は、前記少なくとも1つの第1の心臓信号の前記信号セグメント及び前記少なくとも1つの第2の心臓信号の前記信号セグメントに関する前記最も高い相関係数の和、平均値、及び二乗平均平方根のうちの1つを用いて生成される、請求項15記載のシステム。
  28. 前記演算は、前記1組の信号セグメントを処理して心調律障害の発生源を突き止める演算を更に含む、請求項15記載のシステム。
  29. 複数の心臓信号のセグメントを処理する方法であって、
    最も高い相関係数を決定するために、第1のオフセット範囲に従って、複数の信号の中からほぼ同時期に存在する信号セグメントを自己相関させる工程と、
    前記最も高い相関係数を用いて複合相関係数を生じさせる工程と、
    しきい値よりも低い前記複合相関係数のうちの1つ又は2つ以上を決定する工程と、
    前記ほぼ同時期に存在する信号セグメントに関する最も高い相関係数を決定するために、第2のオフセット範囲に従って、前記決定された1つ又は2つ以上の複合相関係数と関連した前記ほぼ同時期に存在する信号セグメントを自己相関させる工程と
    選択されたオフセット範囲に関して最も高い相関係数と関連した1組のほぼ同時期に存在する信号セグメントを選択する工程と、
    を備え、
    前記第2のオフセット範囲は、前記第1のオフセット範囲とは異なっている、方法。
  30. 前記選択されたオフセット範囲は、100msecから350msecまでの範囲、300msecから500msecまでの範囲、及び、450msecから1300msecまでの範囲から選択される、請求項29記載の方法。
  31. 心臓信号のセグメントを処理するシステムであって、
    処理装置と、
    前記処理装置によって実行されると、前記処理装置が演算を実行する、という命令を記憶する記憶装置と、
    を備え、
    前記演算は、
    最も高い相関係数を決定するために、第1のオフセット範囲に従って、複数の信号の中からほぼ同時期に存在する信号セグメントを自己相関させる工程と、
    前記最も高い相関係数を用いて複合相関係数を生じさせる工程と、
    しきい値よりも低い前記複合相関係数のうちの1つ又は2つ以上を決定する工程と、
    前記ほぼ同時期に存在する信号セグメントに関する最も高い相関係数を決定するために、第2のオフセット範囲に従って、前記決定された1つ又は2つ以上の複合相関係数と関連した前記ほぼ同時期に存在する信号セグメントを自己相関させる工程と
    選択されたオフセット範囲に関して最も高い相関係数と関連した1組のほぼ同時期に存在する信号セグメントを選択する工程と、
    を備え、
    前記第2のオフセット範囲は、前記第1のオフセット範囲とは異なっている、システム。
  32. 前記選択されたオフセット範囲は、100msecから350msecまでの範囲、300msecから500msecまでの範囲、及び、450msecから1300msecまでの範囲から選択される、請求項31記載のシステム。
JP2016542819A 2013-09-12 2014-09-12 生物学的調律障害の分析のために信号セグメントを選択するシステム及び方法 Ceased JP2016530053A (ja)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US201361877093P 2013-09-12 2013-09-12
US61/877,093 2013-09-12
US14/483,914 US10368764B2 (en) 2013-09-12 2014-09-11 System and method to select signal segments for analysis of a biological rhythm disorder
US14/483,914 2014-09-11
PCT/US2014/055419 WO2015038909A2 (en) 2013-09-12 2014-09-12 System and method to select signal segments for analysis of a biological rhythm disorder

Publications (2)

Publication Number Publication Date
JP2016530053A true JP2016530053A (ja) 2016-09-29
JP2016530053A5 JP2016530053A5 (ja) 2017-10-26

Family

ID=52626370

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016542819A Ceased JP2016530053A (ja) 2013-09-12 2014-09-12 生物学的調律障害の分析のために信号セグメントを選択するシステム及び方法

Country Status (11)

Country Link
US (1) US10368764B2 (ja)
EP (1) EP3043702A4 (ja)
JP (1) JP2016530053A (ja)
KR (1) KR20160055150A (ja)
CN (1) CN105705087A (ja)
AU (1) AU2014318636A1 (ja)
CA (1) CA2921523A1 (ja)
IL (1) IL244080A0 (ja)
MX (1) MX2016002942A (ja)
RU (1) RU2016108850A (ja)
WO (1) WO2015038909A2 (ja)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102275241B1 (ko) * 2013-12-06 2021-07-12 삼성전자주식회사 신호의 변환 및 복원 방법, 신호의 변환 및 복원 시스템, 원본 신호 변환 방법, 원본 신호 변환 장치 및 원본 신호 복원 방법
JP6389907B2 (ja) * 2015-06-23 2018-09-12 シェンチェン・インスティテューツ・オブ・アドバンスド・テクノロジー・チャイニーズ・アカデミー・オブ・サイエンシーズShenzhen Institutes Of Advanced Technology Chinese Academy Of Sciences Gpuによる心電信号の並列解析方法
CN111096740B (zh) * 2018-10-25 2022-04-01 上海微创电生理医疗科技股份有限公司 心电信号分析方法及装置、信号记录仪和三维标测系统
CN116322516A (zh) * 2020-09-30 2023-06-23 波士顿科学医学有限公司 Egm特征度量的交互式2d散点图
JP2023543842A (ja) * 2020-09-30 2023-10-18 ボストン サイエンティフィック サイムド,インコーポレイテッド 局所サイクル長及び局所デューティサイクルを決定するための電気生理学システム及び方法
JP2023544021A (ja) * 2020-09-30 2023-10-19 ボストン サイエンティフィック サイムド,インコーポレイテッド 局所サイクル長のグローバルヒストグラム

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030100923A1 (en) * 2001-11-27 2003-05-29 Anders Bjorling Method and circuit for detecting cardiac rhythm abnormalities by analyzing time differences between unipolar signals from a lead with a multi-electrode tip
JP2007531542A (ja) * 2003-07-11 2007-11-08 シー・アール・バード・インコーポレーテッド 心電信号の処理及び表示のための多色重畳システム
JP2008515486A (ja) * 2004-09-30 2008-05-15 カーディアック ペースメーカーズ,インコーポレイテッド 不整脈分類および療法選択

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5109862A (en) 1990-03-19 1992-05-05 Del Mar Avionics Method and apparatus for spectral analysis of electrocardiographic signals
US5217021A (en) 1991-07-30 1993-06-08 Telectronics Pacing Systems, Inc. Detection of cardiac arrhythmias using correlation of a cardiac electrical signals and temporal data compression
DE69514238T2 (de) * 1994-08-19 2000-05-11 Biosense Inc Medizinisches diagnose-, behandlungs- und darstellungssystem
US5954665A (en) * 1995-06-07 1999-09-21 Biosense, Inc. Cardiac ablation catheter using correlation measure
US6944495B2 (en) 2000-11-10 2005-09-13 C.R. Bard, Inc. Methods for processing electrocardiac signals having superimposed complexes
US6526313B2 (en) * 2001-06-05 2003-02-25 Cardiac Pacemakers, Inc. System and method for classifying cardiac depolarization complexes with multi-dimensional correlation
SE0200624D0 (sv) 2002-02-28 2002-02-28 St Jude Medical Medical device
US7567837B2 (en) * 2003-12-19 2009-07-28 Institute Of Critical Care Medicine Enhanced rhythm identification in compression corrupted ECG
US7907994B2 (en) * 2007-01-11 2011-03-15 Biosense Webster, Inc. Automated pace-mapping for identification of cardiac arrhythmic conductive pathways and foci
US8490626B2 (en) * 2008-12-22 2013-07-23 Abbott Cardiovascular Systems Inc. Electrode data analysis to indicate catheter stability and/or needle penetration success
KR20140070502A (ko) * 2011-05-02 2014-06-10 더 리젠츠 오브 더 유니버시티 오브 캘리포니아 형상 절제를 이용한 심장 리듬 장애를 대상으로 하는 시스템 및 방법
US8165666B1 (en) * 2011-05-02 2012-04-24 Topera, Inc. System and method for reconstructing cardiac activation information

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030100923A1 (en) * 2001-11-27 2003-05-29 Anders Bjorling Method and circuit for detecting cardiac rhythm abnormalities by analyzing time differences between unipolar signals from a lead with a multi-electrode tip
JP2007531542A (ja) * 2003-07-11 2007-11-08 シー・アール・バード・インコーポレーテッド 心電信号の処理及び表示のための多色重畳システム
JP2008515486A (ja) * 2004-09-30 2008-05-15 カーディアック ペースメーカーズ,インコーポレイテッド 不整脈分類および療法選択

Also Published As

Publication number Publication date
CN105705087A (zh) 2016-06-22
EP3043702A2 (en) 2016-07-20
WO2015038909A2 (en) 2015-03-19
IL244080A0 (en) 2016-04-21
MX2016002942A (es) 2016-06-17
RU2016108850A (ru) 2017-10-17
US20150073721A1 (en) 2015-03-12
KR20160055150A (ko) 2016-05-17
CA2921523A1 (en) 2015-03-19
US10368764B2 (en) 2019-08-06
WO2015038909A3 (en) 2015-04-30
AU2014318636A1 (en) 2016-03-03
EP3043702A4 (en) 2017-05-03

Similar Documents

Publication Publication Date Title
US10856760B2 (en) Method and system for detection of biological rhythm disorders
US10149622B2 (en) System and method for reconstructing cardiac activation information
US10271786B2 (en) System and method for reconstructing cardiac activation information
JP6408911B2 (ja) 生物学的調律障害に関連付けられた規則性を決定する方法およびシステム
US20130006131A1 (en) System and Method for Targeting Heart Rhythm Disorders Using Shaped Ablation
JP2016530053A (ja) 生物学的調律障害の分析のために信号セグメントを選択するシステム及び方法

Legal Events

Date Code Title Description
RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20170420

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170912

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20170912

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20180720

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180806

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20181101

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20190425

A045 Written measure of dismissal of application [lapsed due to lack of payment]

Free format text: JAPANESE INTERMEDIATE CODE: A045

Effective date: 20190829