JP6818464B2 - 心臓興奮の時系列マッチング及びトラッキング - Google Patents

心臓興奮の時系列マッチング及びトラッキング Download PDF

Info

Publication number
JP6818464B2
JP6818464B2 JP2016157331A JP2016157331A JP6818464B2 JP 6818464 B2 JP6818464 B2 JP 6818464B2 JP 2016157331 A JP2016157331 A JP 2016157331A JP 2016157331 A JP2016157331 A JP 2016157331A JP 6818464 B2 JP6818464 B2 JP 6818464B2
Authority
JP
Japan
Prior art keywords
group
template
transmission time
local
heart
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.)
Active
Application number
JP2016157331A
Other languages
English (en)
Other versions
JP2017035483A (ja
Inventor
ウラジミール・ルビンシュタイン
メイル・バル−タル
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Biosense Webster Israel Ltd
Original Assignee
Biosense Webster Israel Ltd
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 Biosense Webster Israel Ltd filed Critical Biosense Webster Israel Ltd
Publication of JP2017035483A publication Critical patent/JP2017035483A/ja
Application granted granted Critical
Publication of JP6818464B2 publication Critical patent/JP6818464B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/35Detecting specific parameters of the electrocardiograph cycle by template matching
    • 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/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/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/339Displays specially adapted therefor
    • 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
    • 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
    • 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/6867Arrangements 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 specially adapted to be attached or implanted in a specific body part
    • A61B5/6869Heart
    • 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
    • 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/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/361Detecting fibrillation
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Cardiology (AREA)
  • Engineering & Computer Science (AREA)
  • Surgery (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Physiology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Description

(関連出願の相互参照)
本出願は、2015年8月11日に出願された米国仮出願第62/203,422号の利益を主張する。同仮出願は、参照により本明細書に組み込まれる。
(発明の分野)
本発明は、心電図(ECG)信号の分析に関し、より詳細には当該信号の群の分析に関する。
Pascaleらが2013年にCirc Arrhythm Electrophysiolにて公開した「心房細動アブレーション後の心房頻拍の早急な診断に供されるパターンとタイミング(Pattern and Timing of the Coronary Sinus Activation to Guide Rapid Diagnosis of Atrial Tachycardia After Atrial Fibrillation Ablation)」が、本明細書に参照により組み込まれる。当該文献において、著者達は、心房細動中又は後の心房頻拍(AT)により、多くの場合診断を困難となると述べている。著者達の仮説によると、冠状静脈洞(CS)の興奮パターン及びタイミングがATマッピングに役立つ。
本明細書に参照により開示内容が組み込まれるHoubenらに付与された米国特許出願第2015/0208938において、心臓で発生する興奮波前面を、双極型又は単極型電位図により検出することが説明される。
上記説明は、当該分野における関連技術の一般的概論として記載したものであって、この説明に含まれる何らの情報が本特許出願に対する先行技術を構成することを容認するものと解釈するべきではない。
本特許出願中に参照に組み込まれている文献は、本明細書に明確に又は黙示的になされている定義と矛盾する方法で、任意の用語がこれらの組み込まれている文献において定義されている限りにおいて、本明細書中の定義のみが考慮されるべきであることを除いて、本出願の不可欠な部分であると考慮されるものとする。
本発明の実施形態は、方法を提供するものであり、その方法は、
単一の心拍から得られ、被験者の心臓内のカテーテルの複数の電極のそれぞれで生成される心電図(ECG)信号の第1群を受信することと、
前記ECG信号の前記第1群のアノテーション時間に関連するテンプレートを構築し、前記テンプレートにインデックスを割り当てることと、
次の単一の心拍から得られ、前記複数の電極で生成されるECG信号の第2群を受信することと、
前記ECG信号の前記第2群のアノテーション時間を計算することと、
前記テンプレートと、前記第2群の前記アノテーション時間を比較することと、
前記比較により、前記第2群の前記アノテーション時間が前記テンプレートに対応することが示された場合、前記インデックスを前記ECG信号の前記第2群に割り当て、前記心臓からの心拍を表すタイムラインに対する前記テンプレートの発生をグラフとしてディスプレイに提示することを含む。
開示された実施形態は、前記第1群及び前記第2群を使用して、前記心臓の共通マップを構築する。
更に開示された実施形態は、前記比較により、前記第2群の前記アノテーション時間が前記テンプレートに対応しないことが示された場合、前記ECG信号の前記第2群の前記アノテーション時間に関する更なるテンプレートを構築し、更なるインデックスを前記更なるテンプレートに割り当て、前記タイムラインに対する前記更なるテンプレートの更なる発生をグラフとして前記ディスプレイに提示することを含む。典型的には、前記方法は、前記第1群を使用して前記心臓の第1マップを構築し、前記第2群を使用して、前記心臓の前記第1マップと異なる第2マップを構築することを含む。
更に開示された実施形態は、チャネルアノテーションのストリームを受信することと、前記ストリームをそれぞれの心拍に対応する平均、加重平均、又は中央値のいずれかとして基準アノテーションを推定することと、前記基準アノテーションを前記テンプレートの構築に適用することと、を含む。典型的には、前記テンプレートは、前記第1群の前記アノテーション時間と、前記基準アノテーションとの差分を含む。
別の実施形態では、前記カテーテルは、前記心臓の冠状静脈洞内に配置されている。
更に別の実施形態では、前記複数の電極は、前記カテーテルに沿って等間隔に配置されている。
更に別の実施形態では、前記カテーテルは、直線状である。
ある実施形態では、前記テンプレートは第1種の不整脈に割り当てられ、前記方法は、ECG信号の次の群を受信することと、前記次の群の前記アノテーション時間が前記テンプレートに対応せず、したがって前記ECG信号の次の群が前記第1種と異なる第2種の不整脈に対応することを示すかを観察することと、を更に含む。
更に、本発明の実施形態に基づく装置は、
被験者の心臓内に配置されたカテーテルであって、
単一の心拍から得られ、前記カテーテルの複数の電極のそれぞれで生成される心電図(ECG)信号の第1群を受信し、
次の単一の心拍から得られ、前記複数の電極で生成されるECG信号の第2群を受信するように構成されている、カテーテルと、
プロセッサであって、
前記ECG信号の前記第1群のアノテーション時間に関連するテンプレートを構築し、前記テンプレートにインデックスを割り当て、
前記ECG信号の前記第2群のアノテーション時間を計算し、
前記テンプレートと、前記第2群の前記アノテーション時間を比較し、
前記比較により、前記第2群の前記アノテーション時間が前記テンプレートに対応することが示された場合、前記インデックスを前記ECG信号の前記第2群に割り当てるように構成されている、プロセッサと、
前記心臓からの心拍を表すタイムラインに対する前記テンプレートの発生をグラフとして提示するように構成されているディスプレイと、を含む。
本開示は、下記「発明を実施するための形態」のその実施形態から、図面と共により完全に理解されるであろう。
本発明の実施形態に係る、信号分析アルゴリズムの動作を示す概略的ブロック図である。 本発明の実施形態に係る、単極型ECG信号の例をそのアノテーションとその加重平均と共に示す。 本発明の実施形態に係る、図2の信号の時間パターンが提示されたグラフを示す。 本発明の実施形態に係る、アルゴリズムのテンプレートトラッキングモードの工程を示すフローチャートである。 本発明の実施形態に係る、アルゴリズムのテンプレートマッチングモードの工程を示すフローチャートである。 本発明の実施形態に係る、ディスプレイに提示される心拍に応じたテンプレート番号数、アノテーショングループ番号の概略的グラフを示す。 本発明の一実施形態に係る、テンプレートとアノテーショングループとの間の相関を示す。そして、 本発明の一実施形態に係る、テンプレートとアノテーショングループとの間の反相関示す。
概要
心房頻拍(AT)マッピングに対して、冠状静脈洞(CS)興奮の時間パターンが有用となり得る。CS興奮の時間パターンの分析により、マクロリエントリー型ATである可能性の高さによる階層化、順序付けが迅速に可能となり、更に分析はフォーカル型ATの推定された起源を示唆する。
本発明者らは、多チャンネル単極型ECG信号を分析することで、信号の時間パターンをマッチング、トラッキングする方法を開発した。一実施形態においては当該方法により、様々な心臓挙動に特有の所定の「興奮の時間的構成」を特定、追跡可能となる。「興奮の時間的構成」は、本明細書において、「時間パターン」とも称し、ECG信号の単極型アノテーション間の時間的相関に対応する。アノテーションは、信号の局所興奮伝達時間(LAT)に対応する。
ECG信号により導出されたLATマップは、心臓の分析機能に供されてもよく、当該マップは、カリフォルニア州ダイアモンド・バーのBiosense−Webster製のCARTO(登録商標)システムにより構築される。マップは概して、不整脈の心臓における電波生成のいくつかの機構により得られた結果を示す。本発明の実施形態では、単一種の心臓興奮を含むLATマップが構築されるよう、上述の異なる機構を自動で区別、トラッキングするようにしてもよい。本発明の実施形態により、LATマップがより早く、より正確に作製可能となる。
詳細な説明
1.緒言
図1は、本発明の一実施形態による、信号分析アルゴリズムの動作を示す概略的ブロック図である。同アルゴリズム動作は、多数の心電図(ECG)信号を冠状静脈洞(CS)内の基準カテーテル先端14から受信するプロセッサ10により実行され、動作結果はプロセッサによりシステムのディスプレイ18に提示されてもよい。動作結果は、心臓のマップ21に組み込まれてもよく、更に/或いはグラフ24の形態であってもよく、いずれの形態の結果も以下に示す。一実施形態において、カテーテル先端14は、略直線状で、先端に沿って等間隔で配置された10個の電極を有する。いくつかの実施形態では、電極は複数の組として配置される。
より詳細に後述するが、アルゴリズムにより基準カテーテルからのECG信号が分析され、当該分析により、複数の種類の同型心臓興奮が特定される。心臓興奮は、LATマッピング処理に提供され、アルゴリズムを当該マッピング処理の新しいポイントを取得する時期を示すフィルタとしても利用可能である。アルゴリズム、即ちプロセッサ、カテーテル、ディスプレイを実現する要素は、上述のCARTO(登録商標)内に存在するものであってもよい。アルゴリズムはECG信号のリアルタイム分析及び/又は記録信号の分析に使用可能である。
プロセッサは、パターントラッカーブロック22を使用して、ECG信号のアノテーション又はLAT間の特定の時間的相関を示す「活動の時間構成」又は「時間パターン」を判定及び追跡する。図1に示すように、トラッカー動作の際、プロセッサは関心ウィンドウが{ref_Anno(t)−τWOI,ref_Anno(t)+τWOI}の期間を設定する。(式中、tは期間の特定時間、τWOIは期間の開始時間、τWOIは期間の終了時間を表す。)当該期間の間、トラッカーは、プロセッサがECG信号に対して行う波面アノテーション計算の結果を受信する。当該計算は、波面アノテーションブロック26にて実行され、プロセッサは波面アルゴリズムを利用して各ECG信号のアノテーション時間を計算する。適切なアルゴリズムとしては、上述の米国特許出願2015/0208938に記載のものが挙げられる。
トラッカーは、現在の興奮と、予め設定されたテンプレート(テンプレート(n)(nは正の整数のテンプレートインデックスである))との類似性を推定する。テンプレートは、アノテーションの時間パターンの記述子である。
後述のように、プロセッサ10とトラッカーブロック22に動作されるアルゴリズムは、時間パターンが受信したECG信号から計算されて既存の(予め設定された)テンプレートと比較されるテンプレートトラッキング(TT)モードと、信号から新たなテンプレートを生成可能なテンプレートマッチング(TM)モードの2つのモードで動作してもよい。
2.時間パターン計算
複数のCS基準チャネルの時間パターンを計算するため、プロセッサは、所定の心拍(当該心拍を含む所定の期間)annUniChn(k,ti)(kはアノテーション群を形成するチャネルのインデックス)の各チャネルにアノテーションを付ける。そして、プロセッサは、以下の式(1)により、アノテーション群の加重平均wAverage(ti)を計算する。
Figure 0006818464
式中、Uni(k,ti)は、期間{ref_Anno(t)−τWOI,ref_Anno(t)+τWOI}内のk番目のECG信号であり、
annUniChn(k,ti)はk番目のECG信号のアノテーションであり、
Figure 0006818464
はアノテーションannUniChn(k,ti)近傍のECG信号の導関数であり、
そして、
kはチャネルの総数、即ちCSカテーテル先端の電極総数である。
図2は、単極ECG信号と、そのアノテーション(実線円)とその加重平均(黒縦線)の例を示す。そして、プロセッサは郡中の各アノテーションと加重平均との差分(dT(k,ti))を計算する。差分のk次元ベクトルを、当該群の「時間パターン」と称する。図3は、時間パターンのグラフを示す。
3.アルゴリズムの説明
プロセッサ10は、トラッカーブロック22(図1)を使用して、上述のテンプレートトラッキング(TT)モード及びテンプレートマッチング(TM)モードの2つのモードの1つでアルゴリズムを動作させる。通常、TTモードが最初に実施される。
図4は、本発明の実施形態に係るTTモードの行程を示すフローチャートである。第1工程30において、プロセッサは最初に基準アノテーション(ref_Anno(ti))、関心ウィンドウ({τ1WOI,τ2WOI})、時間パターンのテンプレート(Template(n))をブロック22に提供する。テンプレート(n)は、後述のTMモードで作製され、アルゴリズムのユーザにより選択されてもよい。或いは、テンプレート(n)はユーザにより作製されてもよく、その場合、ユーザは選択された単一の心拍からECG信号群を選択し、TMモードにより時間パターンを作製し、これがテンプレート(n)となる。
いくつかの実施形態では、プロセッサにより提供されるref_Anno(t)の値は、アルゴリズムのユーザにより提供される外部入力から求められる。別の実施形態では、プロセッサは自動的にref_Anno(t)の値を推定する。自動的な推定は、典型的には、プロセッサが波面アノテーションブロック26からチャネルアノテーションのストリームを受信し、ストリームをそれぞれの心拍に対応する群にスライシングし、ストリーム内の群のチャネルアノテーションの平均、加重平均、又は中心値をref_Anno(t)として推定する。
第1工程において、プロセッサは更にCS単極型アノテーションの値をブロック22に提供する。アノテーションは、ブロック26の波面アルゴリズムを使用して連続的に計算され、アノテーションのストリームからプロセッサが時間インターバル({ref_Anno(ti)−τ1WOI,ref_Anno(ti)+τ2WOI})に応じて単一の心拍に対して、群を選択する。
類似性ステップ34では、上述の第2項に記載のように、プロセッサは現在の群について時間パターンを計算する。計算により、時間パターンのk次元ベクトルが得られる。そして、プロセッサは現在のグループと、ステップ30にて受信した予め設定されたテンプレート(Template(n))との類似性を推定する。類似性は、以下の式により、判定される。
Figure 0006818464
式中、s(t)は計算された時間パターン{dT(k,ti)}と、テンプレート(n)に対応するベクトル{dT(k,n)}との類似性を表す値であり、
Mはステップ30で受信したアノテーション群のチャネル数である。いくつかの実施形態では、全てのチャネルがアノテーションを提供するわけではないのでM≦Kとなる。
判定ステップ38において、プロセッサは、s(t)が予め設定された閾値未満であるかを判定する。閾値は通常約0.5から約1.2の範囲であるが、本範囲を超えても下回ってもよい。ここで、判定結果は真となると、プロセッサは当該グループがテンプレート(n)に該当すると認識し、バッファ記憶ステップ42に進む。図6を参照に後述するように、プロセッサはその後特定されたテンプレートの結果をディスプレイ18に提示する。その後、フローチャートはステップ30に戻り、次の心拍のアノテーション群が受信される。
判定ステップ38の結果が偽であった場合、フローチャートはステップ30に戻り、プロセッサはCS単極型アノテーション値の新しいグループを処理する。
図5は、本発明の実施形態に係るTMモードのステップを示すフローチャートである。テンプレートトラッキングモードと異なり、テンプレートマッチングモードでは、テンプレートが入力パラメータとして取得されない。図5に示すように、TMモードでは、プロセッサは、上述のTTモードのステップ30と34を繰り返す。判定ステップ68において、TMモードでは、現在までに確認されている全てのテンプレートの内、現在の時間パターンに一致するものが検索される。ステップ68の結果が偽であれば、即ちプロセッサが現在のパターンに一致するテンプレートを発見できなければ、次の新テンプレートステップ76において現在の時間パターンから新たなテンプレートが開かれ、ディスプレイ18上に提示されもよい(図6を参照に後述)。
一致が確認された場合、即ちプロセッサが現在のパターンに一致するテンプレートを発見できた場合、次の再計算ステップ72において、選択されたテンプレートを再計算することで、テンプレート作製に使用された時間パターンのクラスタの新たな一部が反映される。再計算は通常、クラスタにおける複数の群に対し、dT(k,n)の値の平均をとることで、ステップ72で得られた値を時間パターンの計算に組み込むことで行われる。いくつかの実施形態では、クラスタの最初のq個の部分(qは約10であってもよい)にのみ実行される。
図6は本発明の実施形態に係る、ディスプレイ18に提示される概略的グラフである。ステップ42について上述したように、特定されたテンプレートはグラフ24としてディスプレイ18に提示されてもよい(図1)。グラフ24では、テンプレート番号対アノテーショングループ番号が描画されており、アノテーショングループ数は受信した心拍のタイムラインに対応する。グラフ内に記載されているように、約430グループが分析されるが、これは同数の心拍に対応する。更にグラフに記載されているように、タイムライン上の最も多くの心拍がテンプレート2に分類され、いくつかの心拍はテンプレート1に分類され、更にいくつかの心拍はテンプレート3に分類される。
図7は、本発明の一実施形態に係る、テンプレートとアノテーション群との間の相関を示す。同図は、実線で示す心拍に基づくテンプレート412と、破線で示す心拍からのアノテーション群463についてのチャネル番号対dTのグラフであり、テンプレートと、分析されたアノテーション群がかなり類似していることが示されている。
図7に示すようなグラフの対は、心臓の同一領域から2つの心拍が発せられた状況を典型的に示す。
図8は、本発明の一実施形態に係る、テンプレートとアノテーション群との間の反相関を示す。同図は、破線で示す心拍からのアノテーション群123と、実線で示す心拍に基づくテンプレート412との間に相関がないことを示している。むしろ、2つのグラフの傾きは符号が反対となっており、テンプレートと、分析されたアノテーショングループとの間には反相関が存在する。
図8に示すようなグラフの対は、典型的には2つの心拍が、心臓の互いに離れた2つの領域から得られた場合に生じる。例えば、当該状態は一方の領域が右心房であり、もう一方の領域が左心房である場合に生じ得る。
心腔内のLATのマップのような心臓電気生理学的マップは、典型的にはECG信号組を取得し、信号の結果、即ち測定されたLATをマップに組み込むことで作製される。ただし、ECG信号の組は典型的には1つ以上のテンプレートを有する信号を含んでいるため、作製されたマップの不精確さの原因となる。
本発明の実施形態では、異なる複数組の信号は、それぞれ異なるテンプレートに対応するサブセットに分割できる。各サブセットは、心臓の各LATマップの構築に用いられてもよい。図6の説明に戻ると、グラフに示す3つの異なるテンプレートに対して、取得した信号を3つの異なるサブセットに分割でき、各サブセットについてマップが作製される。典型的には、それぞれ異なるテンプレートに対応する異なる複数のサブセットは、異なる複数種類の不整脈に対応する。
更に、不整脈の変化はテンプレートの変化に対応するため、上述の通りにテンプレートとしてECG信号セットを分析することで。本発明の実施形態を不整脈の変化の観察に使用してもよい。
上記されている実施形態は、例として列挙されており、本発明は、上記で特に示され、記載されているものに限定されないことが認識されるであろう。むしろ、本発明の範囲は、上記されている種々の特徴の組み合わせ及び部分的組み合わせと、前述の説明を読むことに基づいて当業者が想起するであろう、先行技術に開示されていない変形例及び修飾との両方を含む。
〔実施の態様〕
(1) 方法であって、
単一の心拍から得られ、被験者の心臓内のカテーテルの複数の電極のそれぞれで生成される心電図(ECG)信号の第1群を受信することと、
前記ECG信号の前記第1群のアノテーション時間に関連するテンプレートを構築し、前記テンプレートにインデックスを割り当てることと、
次の単一の心拍から得られ、前記複数の電極で生成されるECG信号の第2群を受信することと、
前記ECG信号の前記第2群のアノテーション時間を計算することと、
前記テンプレートと、前記第2群の前記アノテーション時間を比較することと、
前記比較により、前記第2群の前記アノテーション時間が前記テンプレートに対応することが示された場合、前記インデックスを前記ECG信号の前記第2群に割り当て、前記心臓からの心拍を表すタイムラインに対する前記テンプレートの発生をグラフとしてディスプレイに提示することと、を含む、方法。
(2) 前記第1群及び前記第2群を使用して、前記心臓の共通マップを構築することを含む、実施態様1に記載の方法。
(3) 前記比較により、前記第2群の前記アノテーション時間が前記テンプレートに対応しないことが示された場合、前記ECG信号の前記第2群の前記アノテーション時間に関する更なるテンプレートを構築し、更なるインデックスを前記更なるテンプレートに割り当て、前記タイムラインに対する前記更なるテンプレートの更なる発生をグラフとして前記ディスプレイに提示することを更に含む、実施態様1に記載の方法。
(4) 前記第1群を使用して前記心臓の第1マップを構築し、前記第2群を使用して、前記心臓の前記第1マップと異なる第2マップを構築することを含む、実施態様3に記載の方法。
(5) チャネルアノテーションのストリームを受信することと、前記ストリームをそれぞれの心拍に対応する区間にスライシングすることと、前記ストリームの前記区間の前記チャネルアノテーションの平均、加重平均、又は中央値のいずれかとして基準アノテーションを推定することと、前記基準アノテーションを前記テンプレートの構築に適用することと、を含む、実施態様1に記載の方法。
(6) 前記テンプレートは、前記第1群の前記アノテーション時間と、前記基準アノテーションとの差分を含む、実施態様5に記載の方法。
(7) 前記カテーテルは、前記心臓の冠状静脈洞内に配置される、実施態様1に記載の方法。
(8) 前記複数の電極は、前記カテーテルに沿って等間隔に配置されている、実施態様1に記載の方法。
(9) 前記カテーテルは、直線状である、実施態様1に記載の方法。
(10) 前記テンプレートは第1種の不整脈に割り当てられ、前記方法は、ECG信号の次の群を受信することと、前記次の群の前記アノテーション時間が前記テンプレートに対応せず、したがって前記ECG信号の次の群が前記第1種と異なる第2種の不整脈に対応することを示すかを観察することと、を更に含む、実施態様1に記載の方法。
(11) 装置であって、
被験者の心臓内に配置されたカテーテルであって、
単一の心拍から得られ、前記カテーテルの複数の電極のそれぞれで生成される心電図(ECG)信号の第1群を受信し、
次の単一の心拍から得られ、前記複数の電極で生成されるECG信号の第2群を受信するように構成されている、カテーテルと、
プロセッサであって、
前記ECG信号の前記第1群のアノテーション時間に関連するテンプレートを構築し、前記テンプレートにインデックスを割り当て、
前記ECG信号の前記第2群のアノテーション時間を計算し、
前記テンプレートと、前記第2群の前記アノテーション時間を比較し、
前記比較により、前記第2群の前記アノテーション時間が前記テンプレートに対応することが示された場合、前記インデックスを前記ECG信号の前記第2群に割り当てるように構成されている、プロセッサと、
前記心臓からの心拍を表すタイムラインに対する前記テンプレートの発生をグラフとして提示するように構成されているディスプレイと、を含む、装置。
(12) 前記プロセッサは、前記第1群及び前記第2群を使用して、前記心臓の共通マップを構築するように構成されている、実施態様11に記載の装置。
(13) 前記比較により、前記第2群の前記アノテーション時間が前記テンプレートに対応しないことが示された場合、前記プロセッサが、前記ECG信号の前記第2群の前記アノテーション時間に関する更なるテンプレートを構築し、更なるインデックスを前記更なるテンプレートに割り当て、前記タイムラインに対する前記更なるテンプレートの更なる発生をグラフとして前記ディスプレイに提示するよう構成されていることを更に含む、実施態様11に記載の装置。
(14) 前記第1群を使用して前記心臓の第1マップを構築し、前記第2群を使用して、前記心臓の前記第1マップと異なる第2マップを構築することを含む、実施態様13に記載の装置。
(15) チャネルアノテーションのストリームを受信し、前記ストリームをそれぞれの心拍に対応する区間にスライシングし、前記ストリームの前記区間の前記チャネルアノテーションの平均、加重平均、又は中央値のいずれかとして基準アノテーションを推定し、前記基準アノテーションを前記テンプレートの構築に適用するように、前記プロセッサを構成することを含む、実施態様11に記載の装置。
(16) 前記テンプレートは、前記第1群の前記アノテーション時間と、前記基準アノテーションとの差分を含む、実施態様15に記載の装置。
(17) 前記カテーテルは、前記心臓の冠状静脈洞内に配置されている、実施態様11に記載の装置。
(18) 前記複数の電極は、前記カテーテルに沿って等間隔に配置されている、実施態様11に記載の装置。
(19) 前記カテーテルは、直線状である、実施態様11に記載の装置。
(20) 前記テンプレートは第1種の不整脈に割り当てられ、前記プロセッサは、ECG信号の次の群を受信し、前記次の群の前記アノテーション時間が前記テンプレートに対応せず、したがって前記ECG信号の次の群が前記第1種と異なる第2種の不整脈に対応することを示すかを観察するように構成されている、実施態様11に記載の装置。

Claims (8)

  1. 装置であって、
    被験者の心臓内に配置されたカテーテルであって、
    単一の心拍から得られ、前記カテーテルの複数の電極のそれぞれで生成される心電図(ECG)信号の第1群を受信し、
    前記単一の心拍以降の単一の心拍から得られ、前記複数の電極で生成されるECG信号の第2群を受信するように構成されている、カテーテルと、
    プロセッサであって、
    前記第1群のそれぞれの局所興奮伝達時間と、所定の期間における前記ECG信号の前記第1群の局所興奮伝達時間の加重平均との差分に基づき、前記局所興奮伝達時間に関連するテンプレートを構築し、前記テンプレートにインデックスを割り当て、
    前記ECG信号の前記第2群の局所興奮伝達時間を計算し、
    前記テンプレートと、前記第2群のそれぞれの局所興奮伝達時間と、前記所定の期間における前記第2群の前記局所興奮伝達時間の加重平均との差分と、を比較し、
    前記比較により、前記第2群の前記局所興奮伝達時間が前記テンプレートに対応することが示された場合、前記インデックスを前記ECG信号の前記第2群に割り当てるように構成されている、プロセッサと、
    前記心臓からの心拍の発生を表すタイムラインにおいて、前記心拍の発生と前記テンプレートのインデックスとの対応を提示するように構成されているディスプレイと、を含む、装置。
  2. 前記プロセッサは、前記第1群及び前記第2群を使用して、前記心臓の共通マップを構築するように構成されている、請求項に記載の装置。
  3. 前記比較により、前記第2群の前記局所興奮伝達時間が前記テンプレートに対応しないことが示された場合、前記プロセッサが、前記ECG信号の前記第2群の前記局所興奮伝達時間に関する更なるテンプレートを構築し、更なるインデックスを前記更なるテンプレートに割り当て、前記タイムラインにおいて、前記心拍の発生と、前記更なるインデックスとの対応を前記ディスプレイに提示するよう構成されていることを更に含む、請求項に記載の装置。
  4. 前記第1群を使用して前記心臓の第1マップを構築し、前記第2群を使用して、前記心臓の前記第1マップと異なる第2マップを構築することを含む、請求項に記載の装置。
  5. チャネルの局所興奮伝達時間のストリームを受信し、前記ストリームをそれぞれの心拍に対応する区間にスライシングし、前記ストリームの前記区間の前記チャネルの局所興奮伝達時間の平均、加重平均、又は中央値のいずれかとして基準局所興奮伝達時間を推定し、前記基準局所興奮伝達時間を前記テンプレートの構築に適用するように、前記プロセッサを構成することを含む、請求項に記載の装置。
  6. 前記複数の電極は、前記カテーテルに沿って等間隔に配置されている、請求項に記載の装置。
  7. 前記カテーテルは、直線状である、請求項に記載の装置。
  8. 記プロセッサは、前記テンプレートを第1種の不整脈に割り当て、ECG信号の次の群を受信し、前記次の群の前記局所興奮伝達時間が前記テンプレートに対応せず、したがって前記ECG信号の次の群が前記第1種と異なる第2種の不整脈に対応することを示すかを観察するように構成されている、請求項に記載の装置。
JP2016157331A 2015-08-11 2016-08-10 心臓興奮の時系列マッチング及びトラッキング Active JP6818464B2 (ja)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US201562203422P 2015-08-11 2015-08-11
US62/203,422 2015-08-11
US15/206,418 2016-07-11
US15/206,418 US10506939B2 (en) 2015-08-11 2016-07-11 Matching and tracking time sequences of heart activation

Publications (2)

Publication Number Publication Date
JP2017035483A JP2017035483A (ja) 2017-02-16
JP6818464B2 true JP6818464B2 (ja) 2021-01-20

Family

ID=56618054

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016157331A Active JP6818464B2 (ja) 2015-08-11 2016-08-10 心臓興奮の時系列マッチング及びトラッキング

Country Status (7)

Country Link
US (2) US10506939B2 (ja)
EP (1) EP3130285B1 (ja)
JP (1) JP6818464B2 (ja)
CN (1) CN106419896B (ja)
AU (1) AU2016210645A1 (ja)
CA (1) CA2938310A1 (ja)
IL (1) IL247088B (ja)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11389232B2 (en) 2006-06-28 2022-07-19 Kardium Inc. Apparatus and method for intra-cardiac mapping and ablation
US9119633B2 (en) 2006-06-28 2015-09-01 Kardium Inc. Apparatus and method for intra-cardiac mapping and ablation
US8906011B2 (en) 2007-11-16 2014-12-09 Kardium Inc. Medical device for use in bodily lumens, for example an atrium
US10827977B2 (en) 2012-05-21 2020-11-10 Kardium Inc. Systems and methods for activating transducers
US9017320B2 (en) 2012-05-21 2015-04-28 Kardium, Inc. Systems and methods for activating transducers
US9198592B2 (en) 2012-05-21 2015-12-01 Kardium Inc. Systems and methods for activating transducers
US10368936B2 (en) 2014-11-17 2019-08-06 Kardium Inc. Systems and methods for selecting, activating, or selecting and activating transducers
US10722184B2 (en) 2014-11-17 2020-07-28 Kardium Inc. Systems and methods for selecting, activating, or selecting and activating transducers
US10506939B2 (en) * 2015-08-11 2019-12-17 Biosense Webster (Israel) Ltd. Matching and tracking time sequences of heart activation
US11389116B2 (en) * 2018-01-19 2022-07-19 Biosense Webster (Israel) Ltd. Apparatus and method for heartbeat classification based on time sequence and morphology of intracardiac and body surface electrocardiogram (ECG) signals
US11109792B2 (en) * 2018-01-31 2021-09-07 Biosense Webster (Israel) Ltd. Signal color morphology
US20210169359A1 (en) * 2019-12-06 2021-06-10 Biosense Webster (Israel) Ltd. Intra-cardiac pattern matching
US11730414B2 (en) 2020-01-21 2023-08-22 Biosense Webster (Israel) Ltd. Automatic pattern acquisition

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5545186A (en) 1995-03-30 1996-08-13 Medtronic, Inc. Prioritized rule based method and apparatus for diagnosis and treatment of arrhythmias
US6449503B1 (en) 1999-07-14 2002-09-10 Cardiac Pacemakers, Inc. Classification of supraventricular and ventricular cardiac rhythms using cross channel timing algorithm
US6477406B1 (en) 1999-11-10 2002-11-05 Pacesetter, Inc. Extravascular hemodynamic acoustic sensor
US6725085B2 (en) * 2000-09-22 2004-04-20 Armin Schwartzman Method and apparatus for characterizing cardiac tissue from local electrograms
US20040019366A1 (en) 2002-07-25 2004-01-29 Rottenberg William B. System and method for classification of cardiac activity using multiple sensing locations
EP1641523A4 (en) * 2003-07-03 2007-06-27 Univ New York SYSTEM AND METHOD FOR EVALUATION OF CARDIAC ELECTROPHYSIOLOGICAL STABILITY AND MODULATION OF CARDIAC OSCILLATIONS
US7242978B2 (en) 2003-12-03 2007-07-10 Medtronic, Inc. Method and apparatus for generating a template for arrhythmia detection and electrogram morphology classification
US7430446B2 (en) 2005-01-20 2008-09-30 Cardiac Pacemakers, Inc. Methods and apparatuses for cardiac arrhythmia classification using morphology stability
US20090030307A1 (en) * 2007-06-04 2009-01-29 Assaf Govari Intracorporeal location system with movement compensation
US8396541B2 (en) 2007-10-24 2013-03-12 Siemens Medical Solutions Usa, Inc. Signal analysis of cardiac and other patient medical signals
EP2555674B1 (en) * 2010-04-08 2019-06-12 The Regents of The University of California Methods, system and apparatus for the detection, diagnosis and treatment of biological rhythm disorders
US9002442B2 (en) * 2011-01-13 2015-04-07 Rhythmia Medical, Inc. Beat alignment and selection for cardiac mapping
US8700136B2 (en) * 2011-11-11 2014-04-15 Biosense Webster (Israel), Ltd. Accurate time annotation of intracardiac ECG signals
EP2934312A1 (en) 2012-12-20 2015-10-28 Boston Scientific Scimed, Inc. Rotor identification using sequential pattern matching
GB201307211D0 (en) 2013-04-22 2013-05-29 Imp Innovations Ltd Image display interfaces
EP2996550B1 (en) 2013-05-14 2019-07-31 Boston Scientific Scimed Inc. Representation and identification of activity patterns during electro-physiology mapping using vector fields
US9554718B2 (en) * 2014-01-29 2017-01-31 Biosense Webster (Israel) Ltd. Double bipolar configuration for atrial fibrillation annotation
US9380953B2 (en) 2014-01-29 2016-07-05 Biosense Webster (Israel) Ltd. Hybrid bipolar/unipolar detection of activation wavefront
US10506939B2 (en) * 2015-08-11 2019-12-17 Biosense Webster (Israel) Ltd. Matching and tracking time sequences of heart activation

Also Published As

Publication number Publication date
IL247088B (en) 2020-10-29
CA2938310A1 (en) 2017-02-11
IL247088A0 (en) 2016-12-29
EP3130285A1 (en) 2017-02-15
CN106419896A (zh) 2017-02-22
US20200113474A1 (en) 2020-04-16
AU2016210645A1 (en) 2017-03-02
JP2017035483A (ja) 2017-02-16
CN106419896B (zh) 2021-01-12
US10506939B2 (en) 2019-12-17
EP3130285B1 (en) 2022-03-02
US20170042443A1 (en) 2017-02-16

Similar Documents

Publication Publication Date Title
JP6818464B2 (ja) 心臓興奮の時系列マッチング及びトラッキング
JP6072005B2 (ja) 信号の加算平均
JP6559439B2 (ja) 複数チャネル心電図信号からの基準注釈時間の決定
US10524679B2 (en) System and method for determining electrogram morphology recurrence patterns and rates during atrial fibrillation
CN108937915B (zh) 一种动态心电图中早搏的识别方法
US11890103B2 (en) Identifying ECG signals having the same morphology
JP6309205B2 (ja) 心臓活動時間の検出
US10335052B2 (en) Detection of pulmonary vein isolation
US20210378579A1 (en) Local noise identification using coherent algorithm
CN105708442B (zh) 心室远场衰减
KR20130091574A (ko) 심방세동 예측 모델 생성장치 및 방법과, 심방세동 예측장치 및 방법
EP3821808A1 (en) A method to find abnormal activations in intra-cardiac electrocardiograms
Thorpe et al. Comparative study of T-amplitude features for fitness monitoring using the ePatch® ECG recorder
EP4311490A1 (en) Method and system for identification of fractionated signals
Valinoti et al. Towards a repository of synthetic electrograms for atrial activation detection in atrial fibrillation
US20230347141A1 (en) Bystander atrium detection using coronary sinus (cs) signals
JP2023061917A (ja) 不整脈原性活動間の類似性を用いた電気生理学的(ep)信号のクラスタ化
CN117462142A (zh) 用于识别碎裂信号的方法和系统

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190708

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200630

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200908

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20201125

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20201203

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20201215

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20201228

R150 Certificate of patent or registration of utility model

Ref document number: 6818464

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250