(57)【要約】【課題】  心電図（ＥＣＧ）信号を分析する。【解決手段】  単一の心拍から得られ、被験者の心臓内のカテーテルの複数の電極のそれぞれで生成される心電図（ＥＣＧ）信号の第１群を受信することと、前記ＥＣＧ信号の前記第１群のアノテーション時間に関連するテンプレートを構築し、前記テンプレートにインデックスを割り当てることを含む、方法。前記方法は、次の単一の心拍から得られ、前記電極で生成されるＥＣＧ信号の第２群を受信することと、前記ＥＣＧ信号の前記第２群のアノテーション時間を計算することと、前記テンプレートと、前記第２群の前記アノテーション時間を比較することと、前記比較により、前記第２群の前記アノテーション時間が前記テンプレートに対応することが示された場合、前記インデックスを前記ＥＣＧ信号の前記第２群に割り当て、前記心臓からの心拍を表すタイムラインに対する前記テンプレートの発生をグラフとしてディスプレイに提示することを更に含む。【選択図】  図１
2017035483.tif 000002
【特許請求の範囲】
【請求項１】
方法であって、
単一の心拍から得られ、被験者の心臓内のカテーテルの複数の電極のそれぞれで生成される心電図（ＥＣＧ）信号の第１群を受信することと、
前記ＥＣＧ信号の前記第１群のアノテーション時間に関連するテンプレートを構築し、前記テンプレートにインデックスを割り当てることと、
次の単一の心拍から得られ、前記複数の電極で生成されるＥＣＧ信号の第２群を受信することと、
前記ＥＣＧ信号の前記第２群のアノテーション時間を計算することと、
前記テンプレートと、前記第２群の前記アノテーション時間を比較することと、
前記比較により、前記第２群の前記アノテーション時間が前記テンプレートに対応することが示された場合、前記インデックスを前記ＥＣＧ信号の前記第２群に割り当て、前記心臓からの心拍を表すタイムラインに対する前記テンプレートの発生をグラフとしてディスプレイに提示することと、を含む、方法。
【請求項２】
前記第１群及び前記第２群を使用して、前記心臓の共通マップを構築することを含む、請求項１に記載の方法。
【請求項３】
前記比較により、前記第２群の前記アノテーション時間が前記テンプレートに対応しないことが示された場合、前記ＥＣＧ信号の前記第２群の前記アノテーション時間に関する更なるテンプレートを構築し、更なるインデックスを前記更なるテンプレートに割り当て、前記タイムラインに対する前記更なるテンプレートの更なる発生をグラフとして前記ディスプレイに提示することを更に含む、請求項１に記載の方法。
【請求項４】
前記第１群を使用して前記心臓の第１マップを構築し、前記第２群を使用して、前記心臓の前記第１マップと異なる第２マップを構築することを含む、請求項３に記載の方法。
【請求項５】
チャネルアノテーションのストリームを受信することと、前記ストリームをそれぞれの心拍に対応する区間にスライシングすることと、前記ストリームの前記区間の前記チャネルアノテーションの平均、加重平均、又は中央値のいずれかとして基準アノテーションを推定することと、前記基準アノテーションを前記テンプレートの構築に適用することと、を含む、請求項１に記載の方法。
【請求項６】
前記テンプレートは、前記第１群の前記アノテーション時間と、前記基準アノテーションとの差分を含む、請求項５に記載の方法。
【請求項７】
前記カテーテルは、前記心臓の冠状静脈洞内に配置される、請求項１に記載の方法。
【請求項８】
前記複数の電極は、前記カテーテルに沿って等間隔に配置されている、請求項１に記載の方法。
【請求項９】
前記カテーテルは、直線状である、請求項１に記載の方法。
【請求項１０】
前記テンプレートは第１種の不整脈に割り当てられ、前記方法は、ＥＣＧ信号の次の群を受信することと、前記次の群の前記アノテーション時間が前記テンプレートに対応せず、したがって前記ＥＣＧ信号の次の群が前記第１種と異なる第２種の不整脈に対応することを示すかを観察することと、を更に含む、請求項１に記載の方法。
【請求項１１】
装置であって、
被験者の心臓内に配置されたカテーテルであって、
単一の心拍から得られ、前記カテーテルの複数の電極のそれぞれで生成される心電図（ＥＣＧ）信号の第１群を受信し、
次の単一の心拍から得られ、前記複数の電極で生成されるＥＣＧ信号の第２群を受信するように構成されている、カテーテルと、
プロセッサであって、
前記ＥＣＧ信号の前記第１群のアノテーション時間に関連するテンプレートを構築し、前記テンプレートにインデックスを割り当て、
前記ＥＣＧ信号の前記第２群のアノテーション時間を計算し、
前記テンプレートと、前記第２群の前記アノテーション時間を比較し、
前記比較により、前記第２群の前記アノテーション時間が前記テンプレートに対応することが示された場合、前記インデックスを前記ＥＣＧ信号の前記第２群に割り当てるように構成されている、プロセッサと、
前記心臓からの心拍を表すタイムラインに対する前記テンプレートの発生をグラフとして提示するように構成されているディスプレイと、を含む、装置。
【請求項１２】
前記プロセッサは、前記第１群及び前記第２群を使用して、前記心臓の共通マップを構築するように構成されている、請求項１１に記載の装置。
【請求項１３】
前記比較により、前記第２群の前記アノテーション時間が前記テンプレートに対応しないことが示された場合、前記プロセッサが、前記ＥＣＧ信号の前記第２群の前記アノテーション時間に関する更なるテンプレートを構築し、更なるインデックスを前記更なるテンプレートに割り当て、前記タイムラインに対する前記更なるテンプレートの更なる発生をグラフとして前記ディスプレイに提示するよう構成されていることを更に含む、請求項１１に記載の装置。
【請求項１４】
前記第１群を使用して前記心臓の第１マップを構築し、前記第２群を使用して、前記心臓の前記第１マップと異なる第２マップを構築することを含む、請求項１３に記載の装置。
【請求項１５】
チャネルアノテーションのストリームを受信し、前記ストリームをそれぞれの心拍に対応する区間にスライシングし、前記ストリームの前記区間の前記チャネルアノテーションの平均、加重平均、又は中央値のいずれかとして基準アノテーションを推定し、前記基準アノテーションを前記テンプレートの構築に適用するように、前記プロセッサを構成することを含む、請求項１１に記載の装置。
【請求項１６】
前記テンプレートは、前記第１群の前記アノテーション時間と、前記基準アノテーションとの差分を含む、請求項１５に記載の装置。
【請求項１７】
前記カテーテルは、前記心臓の冠状静脈洞内に配置されている、請求項１１に記載の装置。
【請求項１８】
前記複数の電極は、前記カテーテルに沿って等間隔に配置されている、請求項１１に記載の装置。
【請求項１９】
前記カテーテルは、直線状である、請求項１１に記載の装置。
【請求項２０】
前記テンプレートは第１種の不整脈に割り当てられ、前記プロセッサは、ＥＣＧ信号の次の群を受信し、前記次の群の前記アノテーション時間が前記テンプレートに対応せず、したがって前記ＥＣＧ信号の次の群が前記第１種と異なる第２種の不整脈に対応することを示すかを観察するように構成されている、請求項１１に記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
（関連出願の相互参照）
本出願は、２０１５年８月１１日に出願された米国仮出願第６２／２０３，４２２号の利益を主張する。同仮出願は、参照により本明細書に組み込まれる。
【０００２】
（発明の分野）
本発明は、心電図（ＥＣＧ）信号の分析に関し、より詳細には当該信号の群の分析に関する。
【背景技術】
【０００３】
Ｐａｓｃａｌｅらが２０１３年にＣｉｒｃ  Ａｒｒｈｙｔｈｍ  Ｅｌｅｃｔｒｏｐｈｙｓｉｏｌにて公開した「心房細動アブレーション後の心房頻拍の早急な診断に供されるパターンとタイミング（Pattern and Timing of the Coronary Sinus Activation to Guide Rapid Diagnosis of Atrial Tachycardia After Atrial Fibrillation Ablation）」が、本明細書に参照により組み込まれる。当該文献において、著者達は、心房細動中又は後の心房頻拍（ＡＴ）により、多くの場合診断を困難となると述べている。著者達の仮説によると、冠状静脈洞（ＣＳ）の興奮パターン及びタイミングがＡＴマッピングに役立つ。
【０００４】
本明細書に参照により開示内容が組み込まれるＨｏｕｂｅｎらに付与された米国特許出願第２０１５／０２０８９３８において、心臓で発生する興奮波前面を、双極型又は単極型電位図により検出することが説明される。
【０００５】
上記説明は、当該分野における関連技術の一般的概論として記載したものであって、この説明に含まれる何らの情報が本特許出願に対する先行技術を構成することを容認するものと解釈するべきではない。
【０００６】
本特許出願中に参照に組み込まれている文献は、本明細書に明確に又は黙示的になされている定義と矛盾する方法で、任意の用語がこれらの組み込まれている文献において定義されている限りにおいて、本明細書中の定義のみが考慮されるべきであることを除いて、本出願の不可欠な部分であると考慮されるものとする。
【発明の概要】
【課題を解決するための手段】
【０００７】
本発明の実施形態は、方法を提供するものであり、その方法は、
単一の心拍から得られ、被験者の心臓内のカテーテルの複数の電極のそれぞれで生成される心電図（ＥＣＧ）信号の第１群を受信することと、
前記ＥＣＧ信号の前記第１群のアノテーション時間に関連するテンプレートを構築し、前記テンプレートにインデックスを割り当てることと、
次の単一の心拍から得られ、前記複数の電極で生成されるＥＣＧ信号の第２群を受信することと、
前記ＥＣＧ信号の前記第２群のアノテーション時間を計算することと、
前記テンプレートと、前記第２群の前記アノテーション時間を比較することと、
前記比較により、前記第２群の前記アノテーション時間が前記テンプレートに対応することが示された場合、前記インデックスを前記ＥＣＧ信号の前記第２群に割り当て、前記心臓からの心拍を表すタイムラインに対する前記テンプレートの発生をグラフとしてディスプレイに提示することを含む。
【０００８】
開示された実施形態は、前記第１群及び前記第２群を使用して、前記心臓の共通マップを構築する。
【０００９】
更に開示された実施形態は、前記比較により、前記第２群の前記アノテーション時間が前記テンプレートに対応しないことが示された場合、前記ＥＣＧ信号の前記第２群の前記アノテーション時間に関する更なるテンプレートを構築し、更なるインデックスを前記更なるテンプレートに割り当て、前記タイムラインに対する前記更なるテンプレートの更なる発生をグラフとして前記ディスプレイに提示することを含む。典型的には、前記方法は、前記第１群を使用して前記心臓の第１マップを構築し、前記第２群を使用して、前記心臓の前記第１マップと異なる第２マップを構築することを含む。
【００１０】
更に開示された実施形態は、チャネルアノテーションのストリームを受信することと、前記ストリームをそれぞれの心拍に対応する平均、加重平均、又は中央値のいずれかとして基準アノテーションを推定することと、前記基準アノテーションを前記テンプレートの構築に適用することと、を含む。典型的には、前記テンプレートは、前記第１群の前記アノテーション時間と、前記基準アノテーションとの差分を含む。
【００１１】
別の実施形態では、前記カテーテルは、前記心臓の冠状静脈洞内に配置されている。
【００１２】
更に別の実施形態では、前記複数の電極は、前記カテーテルに沿って等間隔に配置されている。
【００１３】
更に別の実施形態では、前記カテーテルは、直線状である。
【００１４】
ある実施形態では、前記テンプレートは第１種の不整脈に割り当てられ、前記方法は、ＥＣＧ信号の次の群を受信することと、前記次の群の前記アノテーション時間が前記テンプレートに対応せず、したがって前記ＥＣＧ信号の次の群が前記第１種と異なる第２種の不整脈に対応することを示すかを観察することと、を更に含む。
【００１５】
更に、本発明の実施形態に基づく装置は、
被験者の心臓内に配置されたカテーテルであって、
単一の心拍から得られ、前記カテーテルの複数の電極のそれぞれで生成される心電図（ＥＣＧ）信号の第１群を受信し、
次の単一の心拍から得られ、前記複数の電極で生成されるＥＣＧ信号の第２群を受信するように構成されている、カテーテルと、
プロセッサであって、
前記ＥＣＧ信号の前記第１群のアノテーション時間に関連するテンプレートを構築し、前記テンプレートにインデックスを割り当て、
前記ＥＣＧ信号の前記第２群のアノテーション時間を計算し、
前記テンプレートと、前記第２群の前記アノテーション時間を比較し、
前記比較により、前記第２群の前記アノテーション時間が前記テンプレートに対応することが示された場合、前記インデックスを前記ＥＣＧ信号の前記第２群に割り当てるように構成されている、プロセッサと、
前記心臓からの心拍を表すタイムラインに対する前記テンプレートの発生をグラフとして提示するように構成されているディスプレイと、を含む。
【００１６】
本開示は、下記「発明を実施するための形態」のその実施形態から、図面と共により完全に理解されるであろう。
【図面の簡単な説明】
【００１７】
【図１】本発明の実施形態に係る、信号分析アルゴリズムの動作を示す概略的ブロック図である。
【図２】本発明の実施形態に係る、単極型ＥＣＧ信号の例をそのアノテーションとその加重平均と共に示す。
【図３】本発明の実施形態に係る、図２の信号の時間パターンが提示されたグラフを示す。
【図４】本発明の実施形態に係る、アルゴリズムのテンプレートトラッキングモードの工程を示すフローチャートである。
【図５】本発明の実施形態に係る、アルゴリズムのテンプレートマッチングモードの工程を示すフローチャートである。
【図６】本発明の実施形態に係る、ディスプレイに提示される心拍に応じたテンプレート番号数、アノテーショングループ番号の概略的グラフを示す。
【図７】本発明の一実施形態に係る、テンプレートとアノテーショングループとの間の相関を示す。そして、
【図８】本発明の一実施形態に係る、テンプレートとアノテーショングループとの間の反相関示す。
【発明を実施するための形態】
【００１８】
概要
心房頻拍（ＡＴ）マッピングに対して、冠状静脈洞（ＣＳ）興奮の時間パターンが有用となり得る。ＣＳ興奮の時間パターンの分析により、マクロリエントリー型ＡＴである可能性の高さによる階層化、順序付けが迅速に可能となり、更に分析はフォーカル型ＡＴの推定された起源を示唆する。
【００１９】
本発明者らは、多チャンネル単極型ＥＣＧ信号を分析することで、信号の時間パターンをマッチング、トラッキングする方法を開発した。一実施形態においては当該方法により、様々な心臓挙動に特有の所定の「興奮の時間的構成」を特定、追跡可能となる。「興奮の時間的構成」は、本明細書において、「時間パターン」とも称し、ＥＣＧ信号の単極型アノテーション間の時間的相関に対応する。アノテーションは、信号の局所興奮伝達時間（ＬＡＴ）に対応する。
【００２０】
ＥＣＧ信号により導出されたＬＡＴマップは、心臓の分析機能に供されてもよく、当該マップは、カリフォルニア州ダイアモンド・バーのＢｉｏｓｅｎｓｅ−Ｗｅｂｓｔｅｒ製のＣＡＲＴＯ（登録商標）システムにより構築される。マップは概して、不整脈の心臓における電波生成のいくつかの機構により得られた結果を示す。本発明の実施形態では、単一種の心臓興奮を含むＬＡＴマップが構築されるよう、上述の異なる機構を自動で区別、トラッキングするようにしてもよい。本発明の実施形態により、ＬＡＴマップがより早く、より正確に作製可能となる。
【００２１】
詳細な説明
１．緒言
図１は、本発明の一実施形態による、信号分析アルゴリズムの動作を示す概略的ブロック図である。同アルゴリズム動作は、多数の心電図（ＥＣＧ）信号を冠状静脈洞（ＣＳ）内の基準カテーテル先端１４から受信するプロセッサ１０により実行され、動作結果はプロセッサによりシステムのディスプレイ１８に提示されてもよい。動作結果は、心臓のマップ２１に組み込まれてもよく、更に／或いはグラフ２４の形態であってもよく、いずれの形態の結果も以下に示す。一実施形態において、カテーテル先端１４は、略直線状で、先端に沿って等間隔で配置された１０個の電極を有する。いくつかの実施形態では、電極は複数の組として配置される。
【００２２】
より詳細に後述するが、アルゴリズムにより基準カテーテルからのＥＣＧ信号が分析され、当該分析により、複数の種類の同型心臓興奮が特定される。心臓興奮は、ＬＡＴマッピング処理に提供され、アルゴリズムを当該マッピング処理の新しいポイントを取得する時期を示すフィルタとしても利用可能である。アルゴリズム、即ちプロセッサ、カテーテル、ディスプレイを実現する要素は、上述のＣＡＲＴＯ（登録商標）内に存在するものであってもよい。アルゴリズムはＥＣＧ信号のリアルタイム分析及び／又は記録信号の分析に使用可能である。
【００２３】
プロセッサは、パターントラッカーブロック２２を使用して、ＥＣＧ信号のアノテーション又はＬＡＴ間の特定の時間的相関を示す「活動の時間構成」又は「時間パターン」を判定及び追跡する。図１に示すように、トラッカー動作の際、プロセッサは関心ウィンドウが｛ｒｅｆ＿Ａｎｎｏ（ｔｉ）−τ１ＷＯＩ，ｒｅｆ＿Ａｎｎｏ（ｔｉ）＋τ２ＷＯＩ｝の期間を設定する。（式中、ｔｉは期間の特定時間、τ１ＷＯＩは期間の開始時間、τ２ＷＯＩは期間の終了時間を表す。）当該期間の間、トラッカーは、プロセッサがＥＣＧ信号に対して行う波面アノテーション計算の結果を受信する。当該計算は、波面アノテーションブロック２６にて実行され、プロセッサは波面アルゴリズムを利用して各ＥＣＧ信号のアノテーション時間を計算する。適切なアルゴリズムとしては、上述の米国特許出願２０１５／０２０８９３８に記載のものが挙げられる。
【００２４】
トラッカーは、現在の興奮と、予め設定されたテンプレート（テンプレート（ｎ）（ｎは正の整数のテンプレートインデックスである））との類似性を推定する。テンプレートは、アノテーションの時間パターンの記述子である。
【００２５】
後述のように、プロセッサ１０とトラッカーブロック２２に動作されるアルゴリズムは、時間パターンが受信したＥＣＧ信号から計算されて既存の（予め設定された）テンプレートと比較されるテンプレートトラッキング（ＴＴ）モードと、信号から新たなテンプレートを生成可能なテンプレートマッチング（ＴＭ）モードの２つのモードで動作してもよい。
【００２６】
２．時間パターン計算
複数のＣＳ基準チャネルの時間パターンを計算するため、プロセッサは、所定の心拍（当該心拍を含む所定の期間）ａｎｎＵｎｉＣｈｎ（ｋ，ｔｉ）（ｋはアノテーション群を形成するチャネルのインデックス）の各チャネルにアノテーションを付ける。そして、プロセッサは、以下の式（１）により、アノテーション群の加重平均ｗＡｖｅｒａｇｅ（ｔｉ）を計算する。
【００２７】
【数１】
2017035483.tif 000003
式中、Ｕｎｉ（ｋ，ｔｉ）は、期間｛ｒｅｆ＿Ａｎｎｏ（ｔｉ）−τ１ＷＯＩ，ｒｅｆ＿Ａｎｎｏ（ｔｉ）＋τ２ＷＯＩ｝内のｋ番目のＥＣＧ信号であり、
ａｎｎＵｎｉＣｈｎ（ｋ，ｔｉ）はｋ番目のＥＣＧ信号のアノテーションであり、
【数２】
2017035483.tif 000004
はアノテーションａｎｎＵｎｉＣｈｎ（ｋ，ｔｉ）近傍のＥＣＧ信号の導関数であり、
そして、
ｋはチャネルの総数、即ちＣＳカテーテル先端の電極総数である。
【００２８】
図２は、単極ＥＣＧ信号と、そのアノテーション（実線円）とその加重平均（黒縦線）の例を示す。そして、プロセッサは郡中の各アノテーションと加重平均との差分（ｄＴ（ｋ，ｔｉ））を計算する。差分のｋ次元ベクトルを、当該群の「時間パターン」と称する。図３は、時間パターンのグラフを示す。
【００２９】
３．アルゴリズムの説明
プロセッサ１０は、トラッカーブロック２２（図１）を使用して、上述のテンプレートトラッキング（ＴＴ）モード及びテンプレートマッチング（ＴＭ）モードの２つのモードの１つでアルゴリズムを動作させる。通常、ＴＴモードが最初に実施される。
【００３０】
図４は、本発明の実施形態に係るＴＴモードの行程を示すフローチャートである。第１工程３０において、プロセッサは最初に基準アノテーション（ｒｅｆ＿Ａｎｎｏ（ｔｉ））、関心ウィンドウ（｛τ１ＷＯＩ，τ２ＷＯＩ｝）、時間パターンのテンプレート（Ｔｅｍｐｌａｔｅ（ｎ））をブロック２２に提供する。テンプレート（ｎ）は、後述のＴＭモードで作製され、アルゴリズムのユーザにより選択されてもよい。或いは、テンプレート（ｎ）はユーザにより作製されてもよく、その場合、ユーザは選択された単一の心拍からＥＣＧ信号群を選択し、ＴＭモードにより時間パターンを作製し、これがテンプレート（ｎ）となる。
【００３１】
いくつかの実施形態では、プロセッサにより提供されるｒｅｆ＿Ａｎｎｏ（ｔｉ）の値は、アルゴリズムのユーザにより提供される外部入力から求められる。別の実施形態では、プロセッサは自動的にｒｅｆ＿Ａｎｎｏ（ｔｉ）の値を推定する。自動的な推定は、典型的には、プロセッサが波面アノテーションブロック２６からチャネルアノテーションのストリームを受信し、ストリームをそれぞれの心拍に対応する群にスライシングし、ストリーム内の群のチャネルアノテーションの平均、加重平均、又は中心値をｒｅｆ＿Ａｎｎｏ（ｔｉ）として推定する。
【００３２】
第１工程において、プロセッサは更にＣＳ単極型アノテーションの値をブロック２２に提供する。アノテーションは、ブロック２６の波面アルゴリズムを使用して連続的に計算され、アノテーションのストリームからプロセッサが時間インターバル（｛ｒｅｆ＿Ａｎｎｏ（ｔｉ）−τ１ＷＯＩ，ｒｅｆ＿Ａｎｎｏ（ｔｉ）＋τ２ＷＯＩ｝）に応じて単一の心拍に対して、群を選択する。
【００３３】
類似性ステップ３４では、上述の第２項に記載のように、プロセッサは現在の群について時間パターンを計算する。計算により、時間パターンのｋ次元ベクトルが得られる。そして、プロセッサは現在のグループと、ステップ３０にて受信した予め設定されたテンプレート（Ｔｅｍｐｌａｔｅ（ｎ））との類似性を推定する。類似性は、以下の式により、判定される。
【００３４】
【数３】
2017035483.tif 000005
式中、ｓｎ（ｔｉ）は計算された時間パターン｛ｄＴ（ｋ，ｔｉ）｝と、テンプレート（ｎ）に対応するベクトル｛ｄＴ（ｋ，ｎ）｝との類似性を表す値であり、
Ｍはステップ３０で受信したアノテーション群のチャネル数である。いくつかの実施形態では、全てのチャネルがアノテーションを提供するわけではないのでＭ≦Ｋとなる。
【００３５】
判定ステップ３８において、プロセッサは、ｓｎ（ｔｉ）が予め設定された閾値未満であるかを判定する。閾値は通常約０．５から約１．２の範囲であるが、本範囲を超えても下回ってもよい。ここで、判定結果は真となると、プロセッサは当該グループがテンプレート（ｎ）に該当すると認識し、バッファ記憶ステップ４２に進む。図６を参照に後述するように、プロセッサはその後特定されたテンプレートの結果をディスプレイ１８に提示する。その後、フローチャートはステップ３０に戻り、次の心拍のアノテーション群が受信される。
【００３６】
判定ステップ３８の結果が偽であった場合、フローチャートはステップ３０に戻り、プロセッサはＣＳ単極型アノテーション値の新しいグループを処理する。
【００３７】
図５は、本発明の実施形態に係るＴＭモードのステップを示すフローチャートである。テンプレートトラッキングモードと異なり、テンプレートマッチングモードでは、テンプレートが入力パラメータとして取得されない。図５に示すように、ＴＭモードでは、プロセッサは、上述のＴＴモードのステップ３０と３４を繰り返す。判定ステップ６８において、ＴＭモードでは、現在までに確認されている全てのテンプレートの内、現在の時間パターンに一致するものが検索される。ステップ６８の結果が偽であれば、即ちプロセッサが現在のパターンに一致するテンプレートを発見できなければ、次の新テンプレートステップ７６において現在の時間パターンから新たなテンプレートが開かれ、ディスプレイ１８上に提示されもよい（図６を参照に後述）。
【００３８】
一致が確認された場合、即ちプロセッサが現在のパターンに一致するテンプレートを発見できた場合、次の再計算ステップ７２において、選択されたテンプレートを再計算することで、テンプレート作製に使用された時間パターンのクラスタの新たな一部が反映される。再計算は通常、クラスタにおける複数の群に対し、ｄＴ（ｋ，ｎ）の値の平均をとることで、ステップ７２で得られた値を時間パターンの計算に組み込むことで行われる。いくつかの実施形態では、クラスタの最初のｑ個の部分（ｑは約１０であってもよい）にのみ実行される。
【００３９】
図６は本発明の実施形態に係る、ディスプレイ１８に提示される概略的グラフである。ステップ４２について上述したように、特定されたテンプレートはグラフ２４としてディスプレイ１８に提示されてもよい（図１）。グラフ２４では、テンプレート番号対アノテーショングループ番号が描画されており、アノテーショングループ数は受信した心拍のタイムラインに対応する。グラフ内に記載されているように、約４３０グループが分析されるが、これは同数の心拍に対応する。更にグラフに記載されているように、タイムライン上の最も多くの心拍がテンプレート２に分類され、いくつかの心拍はテンプレート１に分類され、更にいくつかの心拍はテンプレート３に分類される。
【００４０】
図７は、本発明の一実施形態に係る、テンプレートとアノテーション群との間の相関を示す。同図は、実線で示す心拍に基づくテンプレート４１２と、破線で示す心拍からのアノテーション群４６３についてのチャネル番号対ｄＴのグラフであり、テンプレートと、分析されたアノテーション群がかなり類似していることが示されている。
【００４１】
図７に示すようなグラフの対は、心臓の同一領域から２つの心拍が発せられた状況を典型的に示す。
【００４２】
図８は、本発明の一実施形態に係る、テンプレートとアノテーション群との間の反相関を示す。同図は、破線で示す心拍からのアノテーション群１２３と、実線で示す心拍に基づくテンプレート４１２との間に相関がないことを示している。むしろ、２つのグラフの傾きは符号が反対となっており、テンプレートと、分析されたアノテーショングループとの間には反相関が存在する。
【００４３】
図８に示すようなグラフの対は、典型的には２つの心拍が、心臓の互いに離れた２つの領域から得られた場合に生じる。例えば、当該状態は一方の領域が右心房であり、もう一方の領域が左心房である場合に生じ得る。
【００４４】
心腔内のＬＡＴのマップのような心臓電気生理学的マップは、典型的にはＥＣＧ信号組を取得し、信号の結果、即ち測定されたＬＡＴをマップに組み込むことで作製される。ただし、ＥＣＧ信号の組は典型的には１つ以上のテンプレートを有する信号を含んでいるため、作製されたマップの不精確さの原因となる。
【００４５】
本発明の実施形態では、異なる複数組の信号は、それぞれ異なるテンプレートに対応するサブセットに分割できる。各サブセットは、心臓の各ＬＡＴマップの構築に用いられてもよい。図６の説明に戻ると、グラフに示す３つの異なるテンプレートに対して、取得した信号を３つの異なるサブセットに分割でき、各サブセットについてマップが作製される。典型的には、それぞれ異なるテンプレートに対応する異なる複数のサブセットは、異なる複数種類の不整脈に対応する。
【００４６】
更に、不整脈の変化はテンプレートの変化に対応するため、上述の通りにテンプレートとしてＥＣＧ信号セットを分析することで。本発明の実施形態を不整脈の変化の観察に使用してもよい。
【００４７】
上記されている実施形態は、例として列挙されており、本発明は、上記で特に示され、記載されているものに限定されないことが認識されるであろう。むしろ、本発明の範囲は、上記されている種々の特徴の組み合わせ及び部分的組み合わせと、前述の説明を読むことに基づいて当業者が想起するであろう、先行技術に開示されていない変形例及び修飾との両方を含む。
【００４８】
〔実施の態様〕
（１）  方法であって、
単一の心拍から得られ、被験者の心臓内のカテーテルの複数の電極のそれぞれで生成される心電図（ＥＣＧ）信号の第１群を受信することと、
前記ＥＣＧ信号の前記第１群のアノテーション時間に関連するテンプレートを構築し、前記テンプレートにインデックスを割り当てることと、
次の単一の心拍から得られ、前記複数の電極で生成されるＥＣＧ信号の第２群を受信することと、
前記ＥＣＧ信号の前記第２群のアノテーション時間を計算することと、
前記テンプレートと、前記第２群の前記アノテーション時間を比較することと、
前記比較により、前記第２群の前記アノテーション時間が前記テンプレートに対応することが示された場合、前記インデックスを前記ＥＣＧ信号の前記第２群に割り当て、前記心臓からの心拍を表すタイムラインに対する前記テンプレートの発生をグラフとしてディスプレイに提示することと、を含む、方法。
（２）  前記第１群及び前記第２群を使用して、前記心臓の共通マップを構築することを含む、実施態様１に記載の方法。
（３）  前記比較により、前記第２群の前記アノテーション時間が前記テンプレートに対応しないことが示された場合、前記ＥＣＧ信号の前記第２群の前記アノテーション時間に関する更なるテンプレートを構築し、更なるインデックスを前記更なるテンプレートに割り当て、前記タイムラインに対する前記更なるテンプレートの更なる発生をグラフとして前記ディスプレイに提示することを更に含む、実施態様１に記載の方法。
（４）  前記第１群を使用して前記心臓の第１マップを構築し、前記第２群を使用して、前記心臓の前記第１マップと異なる第２マップを構築することを含む、実施態様３に記載の方法。
（５）  チャネルアノテーションのストリームを受信することと、前記ストリームをそれぞれの心拍に対応する区間にスライシングすることと、前記ストリームの前記区間の前記チャネルアノテーションの平均、加重平均、又は中央値のいずれかとして基準アノテーションを推定することと、前記基準アノテーションを前記テンプレートの構築に適用することと、を含む、実施態様１に記載の方法。
【００４９】
（６）  前記テンプレートは、前記第１群の前記アノテーション時間と、前記基準アノテーションとの差分を含む、実施態様５に記載の方法。
（７）  前記カテーテルは、前記心臓の冠状静脈洞内に配置される、実施態様１に記載の方法。
（８）  前記複数の電極は、前記カテーテルに沿って等間隔に配置されている、実施態様１に記載の方法。
（９）  前記カテーテルは、直線状である、実施態様１に記載の方法。
（１０）  前記テンプレートは第１種の不整脈に割り当てられ、前記方法は、ＥＣＧ信号の次の群を受信することと、前記次の群の前記アノテーション時間が前記テンプレートに対応せず、したがって前記ＥＣＧ信号の次の群が前記第１種と異なる第２種の不整脈に対応することを示すかを観察することと、を更に含む、実施態様１に記載の方法。
【００５０】
（１１）  装置であって、
被験者の心臓内に配置されたカテーテルであって、
単一の心拍から得られ、前記カテーテルの複数の電極のそれぞれで生成される心電図（ＥＣＧ）信号の第１群を受信し、
次の単一の心拍から得られ、前記複数の電極で生成されるＥＣＧ信号の第２群を受信するように構成されている、カテーテルと、
プロセッサであって、
前記ＥＣＧ信号の前記第１群のアノテーション時間に関連するテンプレートを構築し、前記テンプレートにインデックスを割り当て、
前記ＥＣＧ信号の前記第２群のアノテーション時間を計算し、
前記テンプレートと、前記第２群の前記アノテーション時間を比較し、
前記比較により、前記第２群の前記アノテーション時間が前記テンプレートに対応することが示された場合、前記インデックスを前記ＥＣＧ信号の前記第２群に割り当てるように構成されている、プロセッサと、
前記心臓からの心拍を表すタイムラインに対する前記テンプレートの発生をグラフとして提示するように構成されているディスプレイと、を含む、装置。
（１２）  前記プロセッサは、前記第１群及び前記第２群を使用して、前記心臓の共通マップを構築するように構成されている、実施態様１１に記載の装置。
（１３）  前記比較により、前記第２群の前記アノテーション時間が前記テンプレートに対応しないことが示された場合、前記プロセッサが、前記ＥＣＧ信号の前記第２群の前記アノテーション時間に関する更なるテンプレートを構築し、更なるインデックスを前記更なるテンプレートに割り当て、前記タイムラインに対する前記更なるテンプレートの更なる発生をグラフとして前記ディスプレイに提示するよう構成されていることを更に含む、実施態様１１に記載の装置。
（１４）  前記第１群を使用して前記心臓の第１マップを構築し、前記第２群を使用して、前記心臓の前記第１マップと異なる第２マップを構築することを含む、実施態様１３に記載の装置。
（１５）  チャネルアノテーションのストリームを受信し、前記ストリームをそれぞれの心拍に対応する区間にスライシングし、前記ストリームの前記区間の前記チャネルアノテーションの平均、加重平均、又は中央値のいずれかとして基準アノテーションを推定し、前記基準アノテーションを前記テンプレートの構築に適用するように、前記プロセッサを構成することを含む、実施態様１１に記載の装置。
【００５１】
（１６）  前記テンプレートは、前記第１群の前記アノテーション時間と、前記基準アノテーションとの差分を含む、実施態様１５に記載の装置。
（１７）  前記カテーテルは、前記心臓の冠状静脈洞内に配置されている、実施態様１１に記載の装置。
（１８）  前記複数の電極は、前記カテーテルに沿って等間隔に配置されている、実施態様１１に記載の装置。
（１９）  前記カテーテルは、直線状である、実施態様１１に記載の装置。
（２０）  前記テンプレートは第１種の不整脈に割り当てられ、前記プロセッサは、ＥＣＧ信号の次の群を受信し、前記次の群の前記アノテーション時間が前記テンプレートに対応せず、したがって前記ＥＣＧ信号の次の群が前記第１種と異なる第２種の不整脈に対応することを示すかを観察するように構成されている、実施態様１１に記載の装置。
【図１】
2017035483.tif 000006
【図２】
2017035483.tif 000007
【図３】
2017035483.tif 000008
【図４】
2017035483.tif 000009
【図５】
2017035483.tif 000010
【図６】
2017035483.tif 000011
【図７】
2017035483.tif 000012
【図８】
2017035483.tif 000013
【外国語明細書】
2017035483000001.pdf