(57)【要約】【課題】心臓の電気生理学的マッピングの方法及びシステムを提供すること。【解決手段】方法は、心腔の入力メッシュの表現、心腔の壁組織上の測定位置の組及びこれらの位置で測定された局所興奮時間（ＬＡＴ）のそれぞれの組を受信することを含む。入力メッシュは、正則化された多角形を含む規則メッシュに再メッシュ化される。測定位置の組及びそれぞれのＬＡＴは、正則化された多角形にデータ適合される。正則化された多角形について、それぞれのＬＡＴ値及び壁組織が瘢痕組織を含むそれぞれの確率は反復的に計算され、瘢痕組織を示す規則メッシュ上での電気生理学的（ＥＰ）興奮波が求められる。ＥＰ興奮波及び瘢痕組織を含む、規則メッシュ上にオーバーレイされた電気解剖学的マップが表される。【選択図】図１
2020075132.tif 000002
【特許請求の範囲】
【請求項１】
心腔の入力メッシュの表現、前記心腔の壁組織上の測定位置の組及び前記位置で測定された局所興奮時間（ＬＡＴ）のそれぞれの組を受信するように構成されたインターフェースと、
プロセッサであって、
前記入力メッシュを、正則化された多角形を含む規則メッシュに再メッシュ化し、
前記測定位置の組及びそれぞれのＬＡＴを前記正則化された多角形にデータ適合させ、
前記正則化された多角形について、（ｉ）それぞれのＬＡＴ値及び（ｉｉ）前記壁組織が瘢痕組織を含むそれぞれの確率を反復的に計算して、瘢痕組織を示す前記規則メッシュ上での電気生理学的（ＥＰ）興奮波を求め、
前記規則メッシュ上にオーバーレイされた、前記ＥＰ興奮波及び前記瘢痕組織を含む電気解剖学的マップを表す、
ように構成されている、プロセッサと、
を備える、システム。
【請求項２】
前記正則化された多角形は、正則化された三角形を含む、請求項１に記載のシステム。
【請求項３】
前記プロセッサは、前記ＬＡＴ値、スローネス値及び前記ＥＰ波の伝播についての方程式の組を反復的に解くことによって、前記ＬＡＴ値及び前記確率を反復的に計算するように構成されている、請求項１に記載のシステム。
【請求項４】
前記プロセッサは、瘢痕重み乗数を使用して再計算されたスローネスベクトルを、各正則化された多角形に再割り当てすることにより、前記確率を反復的に計算するように構成されている、請求項１に記載のシステム。
【請求項５】
前記インターフェースは、瘢痕の存在及び二重ＥＰ電位の存在の少なくとも一方を示す、前記測定位置の少なくともいくつかのバイナリタグ付けを受信するように構成されている、請求項１に記載のシステム。
【請求項６】
前記プロセッサは、前記電気解剖学的マップ上に瘢痕組織を示す伝導矢印をオーバーレイすることによって、前記電気解剖学的マップを表すように構成されている、請求項１に記載のシステム。
【請求項７】
心腔の入力メッシュの表現、前記心腔の壁組織上の測定位置の組及び前記位置で測定された局所興奮時間（ＬＡＴ）のそれぞれの組を受信することと、
前記入力メッシュを、正則化された多角形を含む規則メッシュに再メッシュ化することと、
前記測定位置の組及びそれぞれのＬＡＴを前記正則化された多角形にデータ適合させることと、
前記正則化された多角形について、（ｉ）それぞれのＬＡＴ値及び（ｉｉ）前記壁組織が瘢痕組織を含むそれぞれの確率を反復的に計算して、瘢痕組織を示す前記規則メッシュ上での電気生理学的（ＥＰ）興奮波を求めることと、
前記規則メッシュ上にオーバーレイされた、前記ＥＰ興奮波及び前記瘢痕組織を含む電気解剖学的マップを表すことと、
を含む、方法。
【請求項８】
前記正則化された多角形は、正則化された三角形を含む、請求項７に記載の方法。
【請求項９】
前記ＬＡＴ値及び前記確率を反復的に計算することは、前記ＬＡＴ値、スローネス値及び前記ＥＰ波の伝播についての方程式の組を反復的に解くことを含む、請求項７に記載の方法。
【請求項１０】
前記確率を反復的に計算することは、瘢痕重み乗数を使用して再計算されたスローネスベクトルを、各正則化された多角形に再割り当てすることを含む、請求項７に記載の方法。
【請求項１１】
前記測定位置のうちの少なくともいくつかにバイナリタグを付与して、瘢痕の存在及び二重ＥＰ電位の存在の少なくとも一方を示すことを含む、請求項７に記載の方法。
【請求項１２】
前記電気解剖学的マップを表すことは、前記電気解剖学的マップ上に瘢痕組織を示す伝導矢印をオーバーレイすることを含む、請求項７に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
（関連出願の相互参照）
本出願は、その開示が参照により本明細書に援用される、「ＩＴＥＲＡＴＩＶＥ  ＣＯＨＥＲＥＮＴ  ＭＡＰＰＩＮＧ  ＯＦ  ＣＡＲＤＩＡＣ  ＥＬＥＣＴＲＯＰＨＹＳＩＯＬＯＧＩＣＡＬ（ＥＰ）ＡＣＴＩＶＡＴＩＯＮ  ＩＮＣＬＵＤＩＮＧ  ＲＥＥＮＴＲＹ  ＥＦＦＥＣＴＳ」と題する代理人整理番号１００２−１８５４．１の同日付の米国特許出願に関連している。
【０００２】
（発明の分野）
本発明は、概して、電気生理学的マッピングに関し、詳細には心臓の電気生理学的マッピングの方法及びシステムに関する。
【背景技術】
【０００３】
電気生理学的心臓マッピングは、心臓組織の心不整脈の潜在的な原因を特定するためにしばしば使用される。例えば、米国特許出願公開第２０１７／０２８１０３１号は、マルチ電極プローブを生体被験者の心臓内に挿入することと、心臓のそれぞれの位置において同時に、電極からの電位図を記録することと、電位図中のそれぞれの各興奮時間間隔の範囲を定めることと、興奮時間間隔から、電気的伝播波のマップを生成することと、電位図の興奮時間内の局所興奮時間を調整することにより波のコヒーレンスを最大化することと、調整された局所興奮時間を報告することとにより行われる電気解剖学的マッピングを記載している。
【０００４】
別の例として、米国特許出願公開第２０１６／０１０６３７６号は、マルチ電極カテーテルを使用して、それぞれが位置データと局所興奮時間（local activation、ＬＡＴ）データとの両方を含む複数の電気生理学的（ＥＰ）データポイントを収集することによって、心臓興奮波面の局所伝導速度を計算することを記載している。任意のＥＰデータポイントについて、選択されたＥＰデータポイントと少なくとも２つの追加のＥＰデータポイントとを含むＥＰデータポイントの近傍を定義することができる。近傍内のＥＰデータポイントの位置及びＬＡＴをそれぞれ使用して、位置及びＬＡＴの複数の平面を定義することができる。位置及びＬＡＴの平面の交点から、伝導速度を計算することができる。求められた複数の伝導速度を、例えば三次元心臓モデル上に均一な格子状に配置された複数のベクトルアイコンを表示することにより、グラフィカル表現（例えば電気生理学マップ）として出力することができる。
【０００５】
米国特許第６，３０１，４９６号は、心臓などの生物学的構造の異常状態を診断する方法であって、生物学的構造の表面上の少なくとも３つのサンプリングポイントの生理学的反応を測定し、その反応に関連するベクトル関数を計算し、そのベクトル関数の表現を表示し、その表現から異常状態を推測するステップを含む方法を記載している。本発明は、心不整脈の診断に特に有用であり、この場合、生理学的反応が電圧であり、そこから局所興奮時間が推定され、ベクトル関数が局所興奮時間の勾配、具体的には伝導速度である。伝導速度の大きさは、瘢痕組織では異常に低いことが予想される。
【発明の概要】
【課題を解決するための手段】
【０００６】
本発明の一実施形態は、心腔の入力メッシュの表現、心腔の壁組織上の測定位置の組、及びこれらの位置で測定された局所興奮時間（ＬＡＴ）のそれぞれの組を受信することを含む方法を提供する。入力メッシュは、正則化された多角形を含む規則メッシュに再メッシュ化される。測定位置の組及びそれぞれのＬＡＴは、正則化された多角形にデータ適合される。正則化された多角形について、それぞれのＬＡＴ値及び壁組織が瘢痕組織を含むそれぞれの確率は反復的に計算され、瘢痕組織を示す規則メッシュ上の電気生理学的（ＥＰ）興奮波が求められる。ＥＰ興奮波及び瘢痕組織を含む、規則メッシュ上にオーバーレイされた電気解剖学的マップが表される。
【０００７】
いくつかの実施形態では、正則化された多角形は、正則化された三角形を含む。
【０００８】
いくつかの実施形態では、ＬＡＴ値及び確率を反復的に計算することは、ＬＡＴ値、スローネス値及びＥＰ波の伝播についての方程式の組を反復的に解くことを含む。
【０００９】
一実施形態では、確率を反復的に計算することは、瘢痕重み乗数を使用して再計算されたスローネスベクトルを、各正則化された多角形に再割り当てすることを含む。
【００１０】
別の一実施形態では、この方法は、測定位置の少なくともいくつかにバイナリタグを付与して、瘢痕の存在及び二重ＥＰ電位の存在の少なくとも一方を示すことを更に含む。
【００１１】
いくつかの実施形態では、電気解剖学的マップを表すことは、電気解剖学的マップ上に瘢痕組織を示す伝導矢印をオーバーレイすることを含む。
【００１２】
更に、本発明の一実施形態によれば、インターフェースとプロセッサとを含むシステムが提供される。インターフェースは、心腔の入力メッシュの表現、心腔の壁組織上の測定位置の組及びこれらの位置で測定された局所興奮時間（ＬＡＴ）のそれぞれの組を受信するように構成されている。プロセッサは、入力メッシュを、複数の正則化された多角形を含む規則メッシュに再メッシュ化し、測定位置の組及びそれぞれＬＡＴを複数の正則化された多角形にデータ適合させ、正則化された多角形について、（ｉ）それぞれＬＡＴ値及び（ｉｉ）壁組織が瘢痕組織を含むそれぞれの確率を反復的に計算して、瘢痕組織を示す規則メッシュ上での電気生理学的（ＥＰ）興奮波を求めるように構成されている。プロセッサは更に、規則メッシュ上にオーバーレイされた、ＥＰ興奮波及び瘢痕組織を含む電気解剖学的マップを表すように構成されている。
【００１３】
本発明の別の一実施形態は、心腔の入力メッシュの表現、心腔の壁組織上の測定位置の組及びこれらの位置で測定された局所興奮時間（ＬＡＴ）のそれぞれの組を受信することを含む方法を提供する。
【００１４】
入力メッシュは、正則化された多角形を含む規則メッシュに再メッシュ化される。測定位置の組及びそれぞれのＬＡＴは、正則化された多角形にデータ適合される。正則化された多角形について、それぞれのＬＡＴ値が反復的に計算されて、ＥＰ波のリエントリーを考慮した規則メッシュ上での周期的なＥＰ興奮波の解が求められる。規則メッシュ上にオーバーレイされた周期的なＥＰ興奮波を含む電気解剖学的マップが表される。
【００１５】
いくつかの実施形態では、ＬＡＴ値を反復的に計算することは、リエントリーサイクル長値を含めることによりＥＰ波のリエントリーを記述するように構成された複素数値の方程式の組を反復的に解くことを含む。
【００１６】
いくつかの実施形態では、ＬＡＴ値を反復的に計算することは、ＬＡＴ値、スローネス値及びＥＰ波の伝播についての３つの線形方程式の組を反復的に解くことを含む。
【００１７】
一実施形態では、ＬＡＴ値を反復的に計算することは、近接重み乗数を使用して再計算されたスローネスベクトルを、各正則化された多角形に再割り当てすることを含む。
【００１８】
更に、本発明の一実施形態によれば、インターフェースとプロセッサとを含むシステムが提供される。インターフェースは、心腔の入力メッシュの表現、心腔の壁組織上の測定位置の組及びこれらの位置で測定された局所興奮時間（ＬＡＴ）のそれぞれの組を受信するように構成されている。プロセッサは、入力メッシュを、正則化された多角形を含む規則メッシュに再メッシュ化し、測定位置の組及びそれぞれＬＡＴを複数の正則化された多角形にデータ適合させ、正則化された多角形ついて、それぞれのＬＡＴ値を反復的に計算して、ＥＰ波のリエントリーを考慮した規則メッシュ上での周期的なＥＰ興奮波の解を求めるように構成されている。プロセッサは更に、規則メッシュ上にオーバーレイされた周期的なＥＰ興奮波を含む電気解剖学的マップを表すように構成されている。
【００１９】
本発明は、以下の「発明を実施するための形態」を図面と併せて考慮することで、より完全に理解されよう。
【図面の簡単な説明】
【００２０】
【図１】本発明の一実施形態による、心臓の三次元（３Ｄ）ナビゲーション及び電気生理学的（ＥＰ）信号解析システムの概略的な絵画図である。
【図２】本発明の一実施形態による、コヒーレントＥＰ興奮波を計算するための方法及びアルゴリズムを概略的に示すフローチャートである。
【図３】本発明の一実施形態による、図２に記載のプロセスの一部として瘢痕確率を計算するための方法及びアルゴリズムを概略的に示すフローチャートである。
【図４】本発明の一実施形態による、図２に記載のプロセスからＥＰ興奮波マップを生成する方法及びアルゴリズムを概略的に示すフローチャートである。
【図５】本発明の一実施形態による、ＥＰ波のリエントリー及び瘢痕組織を考慮した左心房のコヒーレントＥＰ興奮マップのボリュームレンダリングの概略的な絵画図である。
【図６】本発明の一実施形態による、ＥＰ興奮波の伝播を示す伝導矢印がオーバーレイされた左心房のコヒーレントＥＰ興奮マップのボリュームレンダリングの概略的な絵画図である。
【発明を実施するための形態】
【００２１】
概論
心不整脈は、心拍が不規則な病態の一種である。この種類の病態の中でも重要なグループが、上室性頻拍症（supraventricular tachycardia、ＳＶＴ）であり、これには心房性頻拍、心房細動、心房粗動及び発作性上室性頻拍が含まれる。ＳＶＴでは、上記の名前が示すように、不整脈が心臓の心房から発生する。
【００２２】
ＳＶＴ事象の一部に、瘢痕組織に起因する電気生理学的（ＥＰ）障害がある。例えば、非典型的な心房性頻拍事象は、アブレーション関連の瘢痕を持つ患者でしばしば発生する。瘢痕組織の領域では、ＥＰの興奮が完全にブロックされるか、部分的にブロックされる場合がある（遅伝導領域）。これらの領域を、以降「ブロックライン」又は瘢痕領域と呼ぶ。
【００２３】
以下に説明する本発明の実施形態は、ＳＶＴグループに属する潜在的な不整脈を特定するための反復コヒーレント電気解剖学的（iterative coherent electroanatomical、ＩＣＥＡ）マッピング方法及びシステムを提供する。開示された説明の文脈において、「コヒーレント」とは、伝導速度が周期的かつ連続的である最も良く適合するＥＰ興奮波伝播の記述を導き出すことを意味する。
【００２４】
したがって、特定の解剖学的障壁又は瘢痕組織によって興奮がブロックされる瘢痕などの非伝導性の領域を除き、興奮波の方向の急激な変化又は興奮波速度の急激な変化は数学的に許可されない。したがって、心臓ＥＰ興奮の連続的かつ周期的な（即ち以下に記載のように興奮波がリエントリーする）フィーチャを正しく捕捉することに加えて、開示されたコヒーレントマッピング方法は、（例えば線などの）あらゆる形態の瘢痕領域を示すＥＡマップを生成する。
【００２５】
開示されたコヒーレントマッピング方法は、上記のように、ＥＰ興奮波が、心臓の一方の端から他方の端に移動した後に生理学的障壁（非伝導性の解剖学的構造）を越えて終了するのではなく、心臓内で環状に移動する（即ち周期的なＥＰ興奮波が発生する）リエントリー性頻拍（ＲＴ）などの興奮波リエントリーの発生に対処する。ＥＰ興奮波が（ＥＣＧの連続するピーク間の時間などの）一定のサイクル長で連続的に伝播する場合には、心腔内の波の伝播は、あるサイクルの「終わり近くの（late）」波面が次のサイクルの「最初の方の（early）」波面と交わるリエントリーサイクル全体にわたり得る。その結果、これらの波面が時間的に近接しているにもかかわらず、１サイクル長離れた誤った局所興奮時間（ＬＡＴ）が領域に割り当てられる可能性がある。
【００２６】
開示されたＩＣＥＡ方法の実施形態は、波のリエントリーを考慮することによって、ＥＡマップにおけるリエントリーに起因するアーチファクトを克服することから、ＩＣＥＡはほとんどの種類の不整脈で機能するように設計されているが、いくつかの実施形態では、開示された方法は心房粗動（atrial flutter、ＡＦＬ）の症例及び他のマクロリエントリーの症例向けに最適化されている。この方法で得られる、例えばＡＦＬの存在のために典型的には異常なＥＰ活動を示す心腔のＥＡマップは、そうでなければ、即ちＥＰ波リエントリーを考慮していないモデルであれば引き起こされる、リエントリーが関与する（以下の記載の）「レインボー」アーチファクトを伴わずに、そのような異常なＥＰ活動を表す。したがって、そのような「レインボー」アーチファクトは、開示されたＩＣＥＡマッピング方法を使用することによって解消される。
【００２７】
いくつかの実施形態では、開示されたモデルは、心腔の入力メッシュの表現、心腔の壁組織上の測定位置の組及びこれらの位置で測定されたＬＡＴのそれぞれの組を受信する。また、位置及びそれぞれのＬＡＴ値を、以降「データポイント」と称する。プロセッサは、入力メッシュを、以降「面」とも称する三角形などの正則化された多角形を含む規則メッシュに再メッシュ化する。いくつかの実施形態では、プロセッサは、測定されたデータポイントの組の正則化された多角形へのデータ適合を実行する。
【００２８】
最初の反復コヒーレント解が求められた後に、そうでなければ連続的であるＩＣＥＡマップの不連続性を示す情報が使用される。例えば、プロセッサは、正常組織領域、瘢痕組織領域又は二重電位の組織領域を示すデータポイントのユーザ提供のタグ付けを利用する。開示された説明においては、二重ＥＰ電位は、心房などの心腔の組織壁上で記録された細分化された多成分の電位図で特定可能な特定の二重スパイクとして定義される。
【００２９】
本発明のいくつかの実施形態では、プロセッサは、各面（即ち三角形）に、電気的に非伝導性の領域、即ち瘢痕組織である確率を更に割り当てる。確率は、マップ上の測定ポイントが正常、瘢痕又は二重電位のいずれの測定結果であるのかを考慮することにより（自動的に又はユーザによって手動で）割り当てられる。割り当ては、ＬＡＴパターン及びスローネス値を考慮することを更に含んでもよい。スローネスは、２つの隣り合う面を結ぶ測地線に沿ったＥＰ信号速度ベクトルの大きさの逆数であり、方向が以下に記載のように速度ベクトルと同じであるベクトルである。
【００３０】
ある面に非伝導性である高い確率が割り当てられると、緩やかなスローネス連続性制約条件及び緩やかな境界速度制約条件が適用され、瘢痕組織の存在の可能性を考慮したマップが生成される。（ブロックラインを含む）遅伝導の領域は、相互に対向する伝播ベクトル（即ちそれらの間にブロックラインがあることを示すＩＣＥＡマップ上の逆平行の伝導矢印）を有する興奮波によって更に特定され、これらの領域では伝播ベクトルが瘢痕の周囲を周回し、非伝導性の組織を示す。一実施形態では、ＥＰ信号の振幅が小さい位置に、伝導障壁の一部である高い確率が割り当てられる。
【００３１】
プロセッサは、ＬＡＴ値及び各面が伝導性であるのか又は瘢痕領域の一部であるのかの確率（即ち重み）を反復的に計算することで、瘢痕領域の表示を含む規則メッシュ上でのＥＰ興奮波を推定する。反復的に計算された重みは、ある面が瘢痕である可能性を増加又は減少させるデータを推定し、各面の計算された瘢痕確率を出力する。その後、瘢痕重み（即ち瘢痕確率）は、ブロックライン領域の連続性の方程式の強度を低下させる乗算重みとして線形方程式に組み込まれる。同じ乗算重みにより、スローネスに対する制約も軽減される。これにより、ブロックライン領域における解は（例えば突然になど）より自由に変化することが可能となる。プロセッサは、ＥＰ興奮波が瘢痕領域に遭遇しない限り連続的であるという仮定に基づいて反復計算を実行する。
【００３２】
次いで、プロセッサは、例えば瘢痕の存在のために典型的には少なくとも部分的に異常なＥＰ興奮波を記述するコヒーレント興奮マップを生成する。一実施形態では、プロセッサは、生成されたＩＣＥＡマップを、以下に記載の「レインボーアーチファクト」を引き起こす「Ｅａｒｌｙ  Ｍｅｅｔｓ  Ｌａｔｅ」への表示依存度を抑えるために、周期的なカラースケールを使用して表す。別の一実施形態では、コヒーレントＥＰ興奮マップは、ＥＰ興奮波の正常な伝播又は異常な伝播を示す伝導矢印を含む。
【００３３】
典型的には、プロセッサは、プロセッサが上記のプロセッサ関連のステップ及び機能のそれぞれを実行することを可能にする特定のアルゴリズムを含むソフトウェア内にプログラムされている。
【００３４】
開示された反復コヒーレントＥＰマッピング方法は、医師が複雑な不整脈を容易かつ確実にマッピング及び解釈し得るマッピングシステムを作成する。したがって、開示されたシステム及び方法は、後続のカテーテルアブレーションなどの侵襲的心臓治療の成功率を高め得るツールを提供する。
【００３５】
システムの説明
図１は、本発明の一実施形態による、三次元（３Ｄ）ナビゲーション及び電気生理学的（ＥＰ）信号解析システム２０の概略的な絵画図である。システム２０は、実質的に任意の生理学的パラメータ又はそのようなパラメータの組み合わせを解析するように構成され得る。本明細書において、例として、解析される信号は、心臓内及び／又は心臓外（体表）心電図（electrocardiogram、ＥＣＧ）の電位−時間関係であると想定される。かかる関係を完全に特性化するために、例えばＬＡＴマップ生成中に行われるように、様々な位置において、時間内に相互に信号を参照する必要がある。時間参照は、例えばＥＣＧ基準信号の各ＱＲＳ群の開始（即ち各心拍の開始）などの（例えば時点などの）基準時間と比較して測定することによって実現される。一実施形態では、基準信号は、冠状静脈洞内に配置されたカテーテルから受信される。ＳＶＴの場合、基準信号は心房活動を表し、その頻度は心室活動の頻度の２倍又は３倍になり得る。ＬＡＴマップを生成するための方法は、前述の米国特許第９，０５０，０１１号に記載されている。
【００３６】
簡潔性及び明瞭性のために、以下の説明は、特に記述のない限り、システム２０がプローブ２４を使用して心臓３４の実際の電気活動を測定する調査手順を想定する。プローブの遠位端３２は、電極２２を有すると想定される。測定した信号は、特に、患者２６の心臓３４の壁組織の少なくとも一部のＬＡＴマップを作成するために使用される。
【００３７】
典型的には、プローブ２４は、システム２０を使用する医師２８によって実施されるマッピング手技の間に患者２６の身体の中に挿入されるカテーテルを備える。手技中、患者２６は、接地電極２３に取り付けられると想定される。加えて、電極２９が、心臓３４の領域にて患者２６の皮膚に取り付けられると想定される。
【００３８】
一実施形態では、プローブ２４は、心腔の一部分の上を移動する際に局所的な心内心電図（ＥＣＧ）を取得する。異常なＥＰ興奮波がカテーテル電極の下を通過した時点で、測定された心内ＥＣＧトレースのフィーチャの一部が注記される。これらの時点で、プローブ２４の位置も記録される。
【００３９】
システム２０は、システムプロセッサ４０によって制御されてもよく、システムプロセッサ４０は、メモリ４４と通信する処理ユニット４２を含む。いくつかの実施形態では、システムプロセッサ４０に含まれているメモリ４４が、患者２６の心臓３４の壁組織の少なくとも一部のＬＡＴ及び／又は電圧マップ６２を記憶する。プロセッサ４０は、典型的にはコンソール４６内に載置されており、コンソール４６は、医師２８がプロセッサと対話するために使用する典型的にはマウス又はトラックボールなどのポインティングデバイス３９を含む操作制御部３８を備える。
【００４０】
プロセッサ４０（特に処理ユニット４２）は、不整脈をモデル化するように、プローブ追跡モジュール３０と、ＥＣＧモジュール３６と、ＥＰ興奮解析モジュール３５とを含むソフトウェアを実行して、システム２０及び／又はＥＰ興奮解析モジュール３５を動作させ、（例えばメモリ４４に記憶されているＬＡＴ又は調整されたＬＡＴマップ６２を使用して）開示された解析の少なくとも一部を実行する。
【００４１】
ＥＣＧモジュール３６は、電極２２及び電極２９からの実際の電気信号を受信するよう連結されている。モジュールは、実際の信号を解析するように構成され、ディスプレイ４８上に、典型的には経時的に変動するグラフィカル表現の標準的なＥＣＧ形式で、解析の結果を提示することができる。
【００４２】
プローブ追跡モジュール３０は、通常、患者２６の心臓内で、プローブ２４の遠位端３２の位置を追跡する。追跡モジュールは、当該技術分野において既知であるいかなるプローブ位置追跡方法を使用してもよい。例えば、モジュール３０は、磁場に基づく位置追跡サブシステムを操作し得る。（簡潔性のために、このようなサブシステムの構成要素は図１には示していない）。
【００４３】
代替的に又は追加的に、追跡モジュール３０は、電極２３と電極２９と電極２２との間のインピーダンス、並びに、プローブ上に位置し得る他の電極に対するインピーダンスを測定することによって、プローブ２４を追跡してもよい。（この場合には、電極２２及び／又は電極２９は、ＥＣＧ信号と位置追跡信号の両方を提供してもよい）。Ｂｉｏｓｅｎｓｅ  Ｗｅｂｓｔｅｒ（カリフォルニア州アーバイン）社により製造されるＣａｒｔｏ３（登録商標）システムは、磁場位置追跡と位置追跡に関するインピーダンス測定の両方を位置追跡に使用する。
【００４４】
プロセッサ４０は、追跡モジュール３０を使用して遠位端３２の位置を測定することができる。加えて、追跡モジュール３０とＥＣＧモジュール３６の両方を使用することによって、プロセッサは、遠位端の位置だけでなく、これらの特定の位置において検出される実際の電気信号のＬＡＴも測定することができる。上述のように、個々の電極２２からの電気追跡信号を、各電極の位置が記録されるように磁気追跡信号と組み合わせることができる。改良された電流位置特定（Advanced Current Location、ＡＣＬ）と名付けられたこのようなハイブリッド（即ち磁気／電気）追跡システム及び方法は、Ｂｉｏｓｅｎｓｅ−Ｗｅｂｓｔｅｒ  Ｉｎｃ．が製造し、その開示が参照により本明細書に援用される米国特許第８，４５６，１８２号に詳細に記載されているＣＡＲＴＯ（登録商標）システムなどの様々な医療用途で実装されている。
【００４５】
プロセッサ４０によって実行される操作の結果は、ディスプレイ４８上で医師２８に表示され、このディスプレイ４８は、典型的には、医師に対するグラフィックユーザインターフェース、電極２２によって感知されたＥＣＧ信号の視覚的表現、及び／又は検査中の心臓３４の画像若しくはマップを表す。
【００４６】
プロセッサ４０によって実行されるソフトウェアは、例えば、電子的な形でネットワークを介してプロセッサ４０にダウンロードされてもよいが、代替的に若しくは付加的に、ソフトウェアが磁気メモリ、光学メモリ若しくは電子メモリなどの非一時的な有形の媒体に提供及び／又は記憶されてもよい。特に、プロセッサ４０は、以下で説明するように、プロセッサ４０が開示されたステップを実行することを可能にする専用のアルゴリズムを実行する。
【００４７】
リエントリー及び瘢痕効果を含む心臓ＥＰ興奮の反復コヒーレントマッピング
図２は、本発明の一実施形態による、コヒーレントＥＰ興奮波を計算するための方法及びアルゴリズムを概略的に示すフローチャートである。本実施形態によるアルゴリズムは、プロセッサ４０によって実行される。
【００４８】
次に、図２のフローチャートを使用してＩＣＥＡ方法のステップ８０〜８４を詳細に説明する前に、これらのステップの概要を説明する。いくつかの実施形態では、開示された方法は、例えば心臓の心房をモデル化し得る入力三次元（３Ｄ）メッシュ（即ち形状）を受信することを含む。３Ｄ入力メッシュは、三角形（即ち三角形のメッシュ）などの多角形を含む。本発明のいくつかの実施形態では、プロセッサは、入力メッシュを、正三角形（以降「面」とも称する）を含む規則メッシュに再メッシュ化する。規則メッシュは、ＥＰ波の連続性を維持する線形方程式の組を定義して反復的に解くことを可能にし、更に以下で説明するように計算の労力も軽減可能とする。
【００４９】
開示された方法は、心房の壁組織上の測定位置と、これらの位置で測定された局所興奮時間（ＬＡＴ）のそれぞれの組とを含むデータポイントの組を受信することを更に含む。いくつかの実施形態では、データポイントの組は、ＰＥＮＴＡＲＡＹ（登録商標）又はＬＡＳＳＯ（登録商標）などのマッピングカテーテルを使用するＣＡＲＴＯ（登録商標）などのカテーテルベースの電気解剖学的マッピングシステムを使用して測定される（システム及びカテーテルは、両方共にカリフォルニア州アーバインのＢｉｏｓｅｎｓｅ−Ｗｅｂｓｔｅｒ製）。
【００５０】
いくつかの実施形態では、開示された方法は、入力として、高速解剖学的マッピング（Fast Anatomical Mapping、ＦＡＭ）技術を使用して得られた解剖学的マップを使用する。プロセッサ４０が解剖学的構造自体以外のマップ内の情報を破棄することにより、開示された方法がより有益な電気解剖学的マップを入力解剖学的マップとして使用してもよい（例えば、心腔の入力電気解剖学的マップのうち使用されるものは心腔形状の元のメッシュである）。
【００５１】
本発明のいくつかの実施形態では、プロセッサは、受信したデータポイントの一部（即ち測定位置及びそれぞれのＬＡＴ）を破棄するか、及び／又は、データポイントが開示された規則メッシュと一致するように、データポイントの一部に対してデータ適合を実行する。
【００５２】
一実施形態では、プロセッサ４０は、組織近接度指数（Tissue Proximity Index、ＴＰＩ）と呼ばれる技術を使用して、マッピング電極が組織と物理的に接触していない間に取得されたと判定されたデータポイントを破棄する。
【００５３】
別の一実施形態では、プロセッサは、そうでない場合に規則メッシュによって提供される解剖学的構造の現実的な記述に適合しない測定位置を調整する。このようにして、開示されたコヒーレントマッピング方法は、カテーテルが房室壁に力（ＥＰ信号を取得する電極を組織と良好に接触させるために加えられる力）を加えて形状を変形させることに起因するか、又は房室壁上のカテーテル位置の変更とは無関係の呼吸運動に起因する誤差などの、カテーテルで測定された位置の誤差を克服する。
【００５４】
次に、プロセッサは入力、即ち測定位置及びそれぞれのＬＡＴ、並びに三角形の面に再メッシュ化された入力メッシュを使用する。一実施形態では、開示されたモデルは、２つの隣り合う面間のＥＰ興奮波が、面の中心を結ぶ形状上の測地線に沿って、そのスローネスベクトルの方向に伝播すると想定する。したがって、面の中心間の時間差は、これらの間の測地線距離にスローネスベクトルを乗算したものである。
【００５５】
このモデルは更に、ＥＰ興奮波が連続的であり、隣り合う面は同様のスローネスベクトルを有する必要があると想定する。このモデルは更に、ＥＰ興奮波の伝導速度が予想される最小速度と最大速度との間（即ち正常な心筋で知られている速度範囲内）であると想定する。
【００５６】
開示されたモデルは、上記の仮定を境界条件下での線形方程式の組で定式化する。プロセッサは最初に方程式を非反復的に解いて、形状の上にオーバーレイされたコヒーレントＥＰ興奮波の初期推定値を求める。データの収集に伴ってＥＰ興奮マップが変化し続けるように、ＩＣＥＡモデルが局所速度ベクトルなどの裏付けとなるデータの蓄積に基づいて再計算されてもよい。一実施形態では、開示されたＩＣＥＡ計算は、通常はすべてのＥＰデータが取得された後にユーザによって開始される。新たなＥＰデータを追加する場合には、ユーザは最初の非反復的な解を含めて最初から計算を再開する。
【００５７】
プロセッサは次に、最初の解を使用して、各面の測定されたＬＡＴ値が、隣り合う面のＬＡＴ値に、測定ポイントから隣り合う面の中心までの距離の関数として影響を与えるという更なる仮定に基づき、「最適な」ＥＰ興奮を見出すことを目的とした最小二乗法で線形方程式の組を反復的に解く。反復計算により、導き出されたＥＰ興奮波を形状の上にオーバーレイしたものから構成される大域的最適解が得られる。典型的には数回までの反復回数は、予め定められた反復回数を超えて解に明らかな変化が生じないように、開発段階での多数のＥＡマップの目視検査に基づく予め設定されたパラメータである。ユーザが計算を開始させ、典型的には数秒以内に完全なＩＣＥＡマップを受信する。
【００５８】
プロセスは、プロセッサ４０が方法のための入力をアップロードするアップロードステップ８０で始まり、入力が（ｉ）心腔壁を表す（三角形の）３Ｄメッシュ、（ｉｉ）（例えば位置及びＬＡＴなどの）測定されたデータポイントの組、（ｉｉｉ）挿入図２８０に示す心房粗動のサイクル長ｔｃｌなどのＥＣＧトレース中に注記されたサイクル長及び、（ｉｖ）瘢痕や二重電位などのタグ付けされた位置を含む。いくつかの実施形態では、例えば、以下の計算がリエントリーを考慮しない場合又は瘢痕関連の計算を含まない場合には、ステップ（ｉｉｉ）又は（ｉｖ）は不要である。
【００５９】
次に、方程式構築ステップ８１において、プロセッサ４０は、メッシュ内の各面（即ち三角形）に対して２つの変数、ＬＡＴ値及びスローネスベクトルを割り当てる。２つの変数の方程式の組を構築するために、方法は２つの変数に影響を与える３つの仮定を行う：
１．ＥＰ興奮波は連続的であり、隣り合う面は同様のスローネスベクトルを有する必要がある。
２．所定の面で測定されたＬＡＴ値は、隣り合う面の計算されたＬＡＴに、測定位置から隣り合う面の重心までの距離の関数として影響を与える。
３．興奮波が開示された伝播の方程式を維持し、波が複数の面の間をそのスローネスベクトルの方向に伝播し、面の中心間の時間差が、これらの間の測地線距離にスローネスベクトルを乗算したものである。
【００６０】
上記の仮定から線形方程式の組が導き出される：
【００６１】
【数１】
2020075132.tif 000003
方程式１（Ｉ）中、以下「近接重み」とも称される重み
【００６２】
【数２】
2020075132.tif 000004
は、面ｊの測定されたＬＡＴ値ｉが計算されたＬＡＴ値に及ぼす相対的な効果を表し、φは、求められるＬＡＴ値を変数として有する（リエントリーを考慮した）周期的なＥＰ興奮波を表す周期的な複素関数である。
【００６３】
方程式１中、方程式（Ｉ）は、使用される測定位置と面の位置との間の距離によって指数関数的に減衰し、規則的な形状からの距離が増加するのに伴って影響が減少する
【００６４】
【数３】
2020075132.tif 000005
を使用して、面ｉでのＬＡＴ測定を面ｊのＬＡＴ値に関連付ける。
【００６５】
【数４】
2020075132.tif 000006
は、典型的には約数ミリメートルの所与の距離を超えるとゼロに設定されるため、これらの方程式のいくつかはヌル方程式である。このように、各ＥＰ測定は、規則的な形状からの距離に比例する相対的な効果を有し、各測定が隣り合う面に影響を及ぼす。規則的な形状に最も近い測定ポイントは、正確である可能性が高くなる。
【００６６】
方程式（ＩＩ）は、後述のように、回転行列Ｒ（ｊ，ｉ）によってリエントリーを処理するために複素数で与えられるスローネスベクトル
【００６７】
【数５】
2020075132.tif 000007
を接続する。
【００６８】
方程式（ＩＩＩ）は、２つの隣り合う面の間を伝播する波は、これらの面間の距離によって関連付けられているスローネス及びＬＡＴを有するべきであることを記述する伝播の方程式である。方程式は、２つの隣り合う面に関して記述されることから、形状内のすべての面に関連付けられたすべてのＬＡＴとスローネスベクトルとの間の関係を反映する。
【００６９】
方程式１（Ｉ）〜（ＩＩＩ）は、以下の関係を適用する：
【００７０】
【数６】
2020075132.tif 000008
【００７１】
方程式２の関係は、すべての時間スケールが短い
【００７２】
【数７】
2020075132.tif 000009
であり、仮定Ｓｉｎ（ｘ）≒ｘ及びＣｏｓ（ｘ）≒１が成り立つと仮定する。方程式２のＴｃｌは、リエントリーのサイクル長である。
【００７３】
時間の関数としての上記の周期的な波の複素数形式は、興奮波が一定のサイクル長で連続的に伝播する上記のリエントリーの発生を考慮している。
【００７４】
方程式１（ＩＩ）は連続伝播の方程式（相似）であり、隣り合う２つの面が同様のスローネスベクトル
【００７５】
【数８】
2020075132.tif 000010
を持つ必要がある。スローネスベクトル
【００７６】
【数９】
2020075132.tif 000011
は、
【００７７】
【数１０】
2020075132.tif 000012
である。Ｒ（ｊ，ｉ）は、隣り合う２つの面ＩとＪとの間の回転行列である。
【００７８】
方程式１（ＩＩＩ）は、面の間の周期的なＥＰ興奮波の伝播の方程式である。波面が移動した距離は、２つの面の重心間の距離である測地線距離の射影、
【００７９】
【数１１】
2020075132.tif 000013
である。
【００８０】
最初の方程式解法ステップ８２で、プロセッサ４０は、方程式１を非反復的に解いて、ＥＰ興奮波の初期推定値を求める。一実施形態では、方程式１は最小二乗法で解かれる。この解法は、ＳＩＡＭ出版社から出版されているＴｒｅｆｅｔｈｅｎ、ＬＮ、Ｂａｕ、Ｄ．著「Ｎｕｍｅｒｉｃａｌ  Ｌｉｎｅａｒ  Ａｌｇｅｂｒａ」（１９９７）の１７２〜１７８ページで説明されているように、コレスキー分解とそれに続く前方及び後方置換を含む数値線形代数からの数学的手法を使用して実行される。
【００８１】
一実施形態では、方程式１の（Ｉ）（ＩＩ）（ＩＩＩ）の項が、単一の疎行列方程式に組み込まれる。各部分方程式（即ち項）に異なる最適化重みを割り当てて、様々な影響因子のバランスをとることが可能である。これらの最適化重みが、反復間で異なっていてもよい。これらの最適化重みは、反復回数のような予め設定されたパラメータでもある。
【００８２】
初期推定値において、スローネスベクトルの物理的意味は、距離に伴う時間の複素位相の変化を表す。これは、距離に伴う時間の変化を表す通常のスローネスベクトルとは対照的である。スローネスが複素数であるため、物理的な定義が異なり、幾分不自然である。
【００８３】
複素領域内で方程式１を解くことで、自然な方法で周期的な興奮に対処し、Ｅａｒｌｙ  Ｍｅｅｔｓ  Ｌａｔｅ効果による複雑さを解消することができる。スローネスの物理的に有意な表現を得るには、反復的な方程式の組を構築して解く追加のステップが必要である。
【００８４】
複素関数を解くことで、周期的な興奮の簡単な計算が可能になるが、スローネスベクトルの物理的な意味を保持するには、複素数のスローネスを実数値に戻す必要があることに留意すべきである。したがって、最終的に、
【００８５】
【数１２】
2020075132.tif 000014
は
【００８６】
【数１３】
2020075132.tif 000015
に変換される。これは、実変数
【００８７】
【数１４】
2020075132.tif 000016
Ｒｅ（φｉ），Ｉｍ（φｉ）についての方程式１を使用して行われる。
【００８８】
本開示の方法は、反復的な方程式構築ステップ８３で、ステップ８２で求められた初期推定値を使用して方程式１を反復的に解くために、
【００８９】
【数１５】
2020075132.tif 000017
φｉについて：
【００９０】
【数１６】
2020075132.tif 000018
を使用して方程式１を線形化し、線形化された方程式１を得る：
【００９１】
【数１７】
2020075132.tif 000019
【００９２】
一実施形態では、実数の線形化された方程式１、即ち方程式４は、線形変数Δφｉ，
【００９３】
【数１８】
2020075132.tif 000020
について解かれる。定数
【００９４】
【数１９】
2020075132.tif 000021
は、前の反復で求められた解である。最初の反復では、
【００９５】
【数２０】
2020075132.tif 000022
は、複素数ＬＡＴと、ステップ８２で求められた初期推定値の実数のスローネスとの解によって初期化される。
【００９６】
各面の伝導速度は、測定された興奮時間、計算された伝播ベクトル及び各面の中心間の既知の距離に基づいて計算される。いくつかの条件下では、線形解により、ゼロに近いスローネスベクトルの領域（「ソース」又は「シンク」とも呼ばれる）が生成される場合がある。このような領域は、無限の速度を持つ位置を意味するため、現実的ではない。ＥＰ伝導速度を何らかの正の生理学的定数ａに近づけるために、制約が加えられる：
【００９７】
【数２１】
2020075132.tif 000023
【００９８】
方程式５は、ａからの大きな速度（小さなスローネス）の逸脱に「罰則を科す」。テイラー展開後に、
【００９９】
【数２２】
2020075132.tif 000024
Δφｉ上の次の制約の組が求められる：
【０１００】
【数２３】
2020075132.tif 000025
【０１０１】
方程式６による制約は、制約が
【０１０２】
【数２４】
2020075132.tif 000026
などの平方根の式を含むために、
【０１０３】
【数２５】
2020075132.tif 000027
の線形形式で記述することができないことから、線形の最初の解が見つかった後にのみ追加される。方程式をデルタ形式で記述する場合には、制約を
【０１０４】
【数２６】
2020075132.tif 000028
へとテイラー展開し、線形形式で記述することができる。
【０１０５】
次に、反復解法ステップ８４で、プロセッサ４０が、変数Ｒｅ（Δφｉ），Ｉｍ（Δφｉ），
【０１０６】
【数２７】
2020075132.tif 000029
を含む方程式４を、最初に再構成部分の表面上にスローネスを投影し、ステップ８２で説明した方法を使用して行列形式で方程式を記述する（即ち、線形ベースを方程式４をより簡単に解くことができるものへと変更する）ことによって解く。
【０１０７】
ステップ８１〜８４では、方程式４が、興奮波の連続伝播を仮定して解かれる。解に瘢痕情報を組み込むには、瘢痕を検出し、ＥＰ興奮波（即ち解）にブロックライン及び非伝導領域で不連続性を生じさせるためのフレームワークが必要である。
【０１０８】
開示されたコヒーレントマッピング方法は、解を非コヒーレント（即ち、一貫性がないか又は不連続）にする制約に、より小さい重みを割り当てる。したがって、プロセッサが無効であると判断した測定ＬＡＴ値（即ち外れ値のＬＡＴ値）は、反復中にプロセッサによって削除される。この方法は、反復の解ベクトルを受信した後に新たな重みを割り当てる。ＬＡＴ制約の場合、このプロセスは一貫性のないポイントを除外する。隣接物の制約の場合、このプロセスは一貫性のないエッジを削除する。
【０１０９】
次に、図２のフローチャートを使用してＩＣＥＡ方法のステップ８５〜９０を詳細に説明する前に、これらのステップの概要を説明する。フレームワークも反復的方法に基づいて構築されるが、ここでは瘢痕重み導入ステップ８５瘢痕（面）重み
【０１１０】
【数２８】
2020075132.tif 000030
が導入され、ステップ８７〜１００で反復的に計算される。重みの範囲は［０．０５，１］で、０．０５は瘢痕を示す。
【０１１１】
反復的に計算された重みは、ある面が瘢痕である可能性を増加又は減少させるデータを推定し、各面の計算された瘢痕確率を出力する。次に、瘢痕重み（即ち瘢痕確率）
【０１１２】
【数２９】
2020075132.tif 000031
が、ブロックライン領域の連続性の方程式の強度を低下させる乗算重みとしてステップ８７に組み込まれる。方程式６によって設定されたものなどのスローネスノルムに対する制約も、同様に、乗算重みによって軽減される。これにより、解がブロックライン領域でより自由に変化することが可能となる。
【０１１３】
まず、伝播の方程式４（ＩＩ）に、エッジ（Ｉ，ｊ）に隣り合う２つの面のそれぞれの２つの重み
【０１１４】
【数３０】
2020075132.tif 000032
の最小値を乗算する：
【０１１５】
【数３１】
2020075132.tif 000033
【０１１６】
これは、面の間の伝導が、これらの面のいずれかがすでに瘢痕として示されている場合（即ち小さな
【０１１７】
【数３２】
2020075132.tif 000034
を持っている場合）には遅いという仮定を反映している。
【０１１８】
したがって、瘢痕が検出された領域では、連続伝播の方程式４０（ＩＩ）の重みが小さくなる。これにより、隣り合う面間のスローネス類似性の要件が緩和され、瘢痕のある領域での波の伝播が不連続となり得る。
【０１１９】
次に、正規化されたスローネス及び速度の方程式（方程式６）に、計算された重みが乗算される。したがって、瘢痕面のある目標速度のサイズの要件が軽減される：
【０１２０】
【数３３】
2020075132.tif 000035
【０１２１】
瘢痕確率が方程式に統合された後に、瘢痕重み統合ステップ８７で、プロセッサ４０が修正された方程式４を解き、ＬＡＴ及びスローネス出力ステップ８９で新しい解が提示される。次に、ステップ１００で、図６で説明したコヒーレントＥＰ興奮マップなどの十分に正確なマップを生成するために検証される所定の反復回数まで、新たなより正確な瘢痕確率が計算される。
【０１２２】
図３は、本発明の一実施形態による、図２に記載のプロセスの一部として瘢痕確率を計算するための方法及びアルゴリズムを概略的に示すフローチャートである。図３は、ステップ１００に含まれる、ＥＰ興奮波の解にブロックライン及び非伝導領域のラインで不連続性を生じさせ得る複数の異なるモジュール（即ち計算サブステップ）を記載している。本実施形態によるアルゴリズムは、プロセッサ４０によって実行される。
【０１２３】
ステップ１００は、メッシュのある三角形の面が瘢痕を表す可能性を増加又は減少させる各面ごとのデータを推定する計算モジュール（即ちステップ）９１〜９８を含む。ステップ９１〜９８のそれぞれの出力は、０〜１の範囲の重みである。重みは［０，１］の範囲内にあり、１は正常な面を示し、０は瘢痕又は「遅伝導」面を示す（現在のバージョンではこれら２つに区別はない）。
【０１２４】
以下で説明するように、ステップ９２〜９６で計算された重みは、コスト関数の行列形式で配列され、ステップ９９は、コスト関数を最小化することにより、瘢痕重み
【０１２５】
【数３４】
2020075132.tif 000036
を計算する。ステップ９９において、プロセッサ４０は、ステップ９１〜９８で収集された情報を各面の複合確率を結合させたものに統合することにより、瘢痕重み
【０１２６】
【数３５】
2020075132.tif 000037
を計算する。
【０１２７】
境界領域の重み付けステップ９２で、プロセッサ４０は、近傍の測定からのＬＡＴの分散の勾配を示す境界領域の面の重み
【０１２８】
【数３６】
2020075132.tif 000038
を計算する。急勾配の分散が、瘢痕に隣接する領域を示してもよい。
【０１２９】
大きなスローネス重み付けステップ９３で、プロセッサ４０は、面が大きなスローネスベクトルを有するかどうかを示すスローネス重み
【０１３０】
【数３７】
2020075132.tif 000039
を計算する。このことは、面が瘢痕である確率を高くする。
【０１３１】
ＬＡＴ関連重み付けステップ９４で、プロセッサ４０は、測定されたＬＡＴと計算された面ＬＡＴとの間の差を重み付けする。これは測定面重みと呼ばれ、
【０１３２】
【数３８】
2020075132.tif 000040
で表される。重み
【０１３３】
【数３９】
2020075132.tif 000041
は、測定ポイントｉと面ｊとの間の「リンク強度」を反映している。
【０１３４】
瘢痕関連重み付けステップ９５で、プロセッサ４０は、瘢痕測定ポイントｉと面ｊとの間の「リンク強度」を反映する重み
【０１３５】
【数４０】
2020075132.tif 000042
を計算する。面ｉと面ｊとの間の空間距離が増加すると、重みは減少する。
【０１３６】
単一瘢痕密度重み付けステップ９６で、プロセッサ４０は、単一の瘢痕測定値が規則的な（非瘢痕）測定値に囲まれている領域を示す単一瘢痕密度重み
【０１３７】
【数４１】
2020075132.tif 000043
を計算する。この重みは、瘢痕の表示があまりにもまばらであるために現実的には有効な測定値ではない面の瘢痕確率を低くするように定義されている。このことは、近傍の瘢痕測定値と近傍の非瘢痕測定値との比を調べることによって行われる。非瘢痕の測定値及び孤立した瘢痕測定値が多数ある場合、このモジュールから生成される重みが、この面が瘢痕である確率を低くする。
【０１３８】
ステップ９２〜９６で計算された複数の異なる重みは、各面の瘢痕確率を計算する瘢痕情報推定統合モジュール９９によって処理される。次に、図２に示すように、瘢痕確率が、ブロックライン領域の連続性の式の強度を低下させるために反復計算のステップ８７に組み込まれる。更に、スローネスノルム及び速度ノルムの式の相対強度も、上述のように低減される。これにより、解がブロックライン領域でより自由に変化することが可能となる。
【０１３９】
ステップ９９は、以下の規則をコスト関数に定式化することにより、各面の瘢痕確率重みを計算する：
１．瘢痕／遅伝導領域が連続している。
２．瘢痕タグが付けられた測定値は、近くに瘢痕があることを意味する。
３．近傍に多くの測定ポイントがある領域は伝導性である。
４．大きな境界領域の重み又は二重電位のタグが付けられた測定値は、瘢痕の境界を示す。
５．周辺情報が少ない、つまり、測定ポイントが近傍にない領域は、瘢痕である確率が高い。
６．瘢痕タグが付けられた測定値に近接し、かつスローネスの大きな領域は、瘢痕である確率がより高い。
【０１４０】
求められたコストは、二次項と線形項の和として定義される：
【０１４１】
【数４２】
2020075132.tif 000044
式中、
【０１４２】
【数４３】
2020075132.tif 000045
はすべての面の重みのベクトル、Ｈ（ｉ，ｊ）Ｈは面の重みの対を含むコスト項の対称行列である：
【０１４３】
【数４４】
2020075132.tif 000046
【０１４４】
方程式９の二次コスト項は、次のように定義される：
【０１４５】
【数４５】
2020075132.tif 000047
【０１４６】
Ｈは、瘢痕確率の連続性を反映し、隣り合う面が同様の瘢痕確率を有するべきであることを意味している。これは、断片化された瘢痕領域を繋げることを目的としている。
【０１４７】
Ｌは、各面に関連する線形コスト項のベクトルである。これは、矛盾する２つの項で構成されており、一方の項には瘢痕を示すすべての情報が含まれ、他方の項には正常な面を示すすべての情報が含まれる。
【０１４８】
次に、線形コストベクトルは次のように定義される：
【０１４９】
【数４６】
2020075132.tif 000048
【０１５０】
各面の線形項
【０１５１】
【数４７】
2020075132.tif 000049
は、瘢痕面又は正常な面の存在を示すものであるか、又はステップ９２〜９６で計算された重みの関数である。ｆは、［０〜１］の間の値を持つ正規化係数である。コスト関数が面ｉに対して正の場合、これはスコアが瘢痕を示していることを意味し、プロセッサ４０はその面に低い重みを割り当てることでスコアを最小化するように動作する。コスト関数が負の場合、これは正常な面を示し、プロセッサ４０はその面に高い重みを割り当てることで重みを最大化するように動作する。これは、線形コスト項を次のように記述することにより行われる：
方程式１３  Ｃｌｉｎｅａｒ＝（ｗｓｆ）ｔＬ
【０１５２】
この項は、最終コスト関数を求めるために二次項（方程式１１）に追加される。オプティマイザが、各面ごとに定義された重みｗｓｆを求める。重みｗｓｆは、１の場合は伝導面を示し、０の場合はブロックを示す。
【０１５３】
例えば、ｆが０．５であると仮定すると、面の瘢痕スコア及び正常スコアである
【０１５４】
【数４８】
2020075132.tif 000050
が１、
【０１５５】
【数４９】
2020075132.tif 000051
が０．２の場合、Ｌｉはゼロよりも大きい０．４となる。この場合、オプティマイザは、総コスト関数を最小化するために、この面に対してゼロに近い重みｗｓｆを割り当てる。ゼロに近い重みｗｓｆは、瘢痕を示す。逆に、
【０１５６】
【数５０】
2020075132.tif 000052
が０で、
【０１５７】
【数５１】
2020075132.tif 000053
が１である場合には、Ｌｉ＝−１となり、コスト関数が負となって、オプティマイザは、正常な面を示す１に近いｗｓｆを割り当てる。負の値と正の値とを使用することで、オプティマイザは、正常な面に１に近い値を割り当て、瘢痕のある面に０に近い値を割り当てることが可能となる。
【０１５８】
プロセッサ４０は、二次計画法を使用して行われる総コスト関数
【０１５９】
【数５２】
2020075132.tif 000054
を最小化することにより、ベクトル
【０１６０】
【数５３】
2020075132.tif 000055
の重みを導き出す。次のように問題を定式化する：
【０１６１】
【数５４】
2020075132.tif 000056
【０１６２】
「Ｉｎｔｅｒｉｏｒ−ｐｏｉｎｔ  ｍｅｔｈｏｄｓ」（Ｊｏｕｒｎａｌ  ｏｆ  Ｃｏｍｐｕｔａｔｉｏｎａｌ  ａｎｄ  Ａｐｐｌｉｅｄ  Ｍａｔｈｅｍａｔｉｃｓ、Ｖｏｌｕｍｅ  １２４、ｉｓｓｕｅｓ（１−２）、２８１〜３０２ページ（２０００）でＦｌｏｒｉａｎ  Ａ．Ｐｏｔｒａ及びＳｔｅｐｈｅｎ  Ｊ．Ｗｒｉｇｈｔにより説明されているように、この問題は内点凸最適化を使用して解決される。このモジュールの出力は、面ｉが正常にステップ８７に組み込まれる確率である瘢痕面の重み
【０１６３】
【数５５】
2020075132.tif 000057
である。
【０１６４】
図３に示す例示的なフローチャートは、単に概念を明確化する目的のために選ばれたものである。例えば、表現を簡単にするために、実装の更なる詳細は省略されている。
【０１６５】
図４は、本発明の一実施形態による、図２に記載のプロセスからＥＰ興奮波マップを生成する方法及びアルゴリズムを概略的に示すフローチャートである。本実施形態によるアルゴリズムは、入力メッシュ受信ステップ２９９で、プロセッサ４０が心腔の三次元メッシュ表現２９９ａを受信することから始まるプロセスを実行する。図示するように、メッシュ２９９ａａは不規則な三角形の面を含む。典型的には前述のＦＡＭ技術によって生成される、大きい辺と小さい辺とを持つ三角形を含むメッシュ２９９ａなどのメッシュは、開示された方法によって行われる幾何学的近似でエラーを引き起こす。次に、入力メッシュ正則化ステップ３００で、プロセッサ４０は入力メッシュ２９９を再メッシュ化して、より規則的なメッシュ３００ａを生成する。再メッシュ化は、開示された方法によって用いられる計算ステップをＥＰ波推定の必要性に応じて最適化する目的でも行われる。一実施形態では、メッシュ解像度（三角形の数）は、生成される方程式の数を（例えば所定のサイズの行列で行われるように）決定する。ステップ３００は、図２及び図３で説明した計算の労力を、方程式の数を少なくすることで削減する。
【０１６６】
次に、メッシュデータ適合ステップ３０２で、プロセッサ４０は、再構成部分を損なわずに、三角形の中心の一部を仮想的に変換して、近傍の測定位置により良好に適合させる。ステップ３０２は、メッシュが作成されるデータ取得中に、マッピングカテーテルが心房壁を押して「膨張した」再構成形状を生成する傾向があることから必要である。測定の多くが再構成部分から離間しているために、測定位置を心腔の形状に関連付ける際に不確実性が生じる。
【０１６７】
次に、プロセッサ４０は、対話式計算ステップ９０で、専用アルゴリズムを使用して、図２及び図３に記載の計算を実行し、最適化された三角形メッシュの各面について、ＬＡＴ、スローネス及び瘢痕確率の値を導き出す。典型的には、瘢痕重みは、ＬＡＴ及びスローネス値が更新される都度更新されるわけではない。
【０１６８】
次に、投影ステップ３０４で、プロセッサ４０は、ステップ３０２で求められた値を元のメッシュ２９９ａ上に投影する。上述のように、ＥＰ興奮波は単純化されたメッシュについて計算される。それでも、開示された方法は、元のメッシュから求められたコヒーレントマップを提供する。マップは、低解像度のメッシュの場合には、面の各中央に定義される。プロセッサ４０は、高解像度メッシュから面の法線方向に、内向き及び外向きに光線を発射することにより、（入力マップ２９９ａの）高解像度マップの面を低解像度メッシュ上のポイントにマッピングする。
【０１６９】
低解像度メッシュに当たる光線が１つだけの場合、プロセッサ４０はそのポイントを低解像度メッシュ上の関連ポイントとして選択する。両方向の光線がメッシュに当たった場合には、プロセッサ４０は最も近いポイントを選択する。
【０１７０】
次に、プロセッサ４０は、光線が当たったポイントで低解像度メッシュの面内のＬＡＴ値を補間し、この値を高解像度メッシュの面に割り当てる。一実施形態では、プロセッサ４０は、取得されたデータポイントの数が所定の必要な値よりも少ない組織領域のＬＡＴ情報を省略する。したがって、ＥＰマッピングがもたらすマッピングがあまりにもまばらであるこのような組織領域は、得られた視覚化では「グレーアウト」される。
【０１７１】
ステップ３０４により、プロセッサ４０は、図６に記載されているコヒーレントＥＰ興奮波マップ４００を生成する。
【０１７２】
不整脈の粗動を理解するための重要なステップは、メッシュ全体の波のサイクルを表示する機能である。２Ｄ画面での３Ｄメッシュの視覚化には、メッシュの裏側を見ることや、メッシュの前方部及び遠位部の判断などの困難が伴う。等時図生成ステップ４１０で、プロセッサ４０はマップ４００を等時図５００へと画像処理する。
【０１７３】
等時図５００は等高線マップ５００ａである、各等高線は、ＬＡＴが一定の曲線である。等時図５００ａを、左心房の前−後（ＡＰ）図を一例として示す。実際には、マップは色分けされている。
【０１７４】
マップ５００ａは透明であるため、ユーザは波のサイクルのスパン全体を一望視できる一方で、前方と遠方の成分は簡単に区別できる。更に、この図により、構造全体のサイクルを簡単にたどることができるため、不整脈をより良く理解することができる。等高線間隔は時間的に等間隔であるため、波の伝播速度は曲線の密度から明らかである。等高線の密度が高い領域は伝播が遅いことを表し、密度が低い領域は伝播が速いことを表す。瘢痕領域は、等高線のない均一な領域として示される。
【０１７５】
図４に示す例示的なフローチャートは、単に概念を明確化する目的のために選ばれたものである。例えば、ステップ３０２などのいくつかのステップの例示的なグラフィック出力は、表現を簡単にするために省略されている。
【０１７６】
図５は、本発明の一実施形態による、ＥＰ波のリエントリー及び瘢痕組織を考慮した左心房のコヒーレントＥＰ興奮マップ４００のボリュームレンダリングの概略的な絵画図である。図示するように、マップ４００は、上述のレインボーアーチファクトが開示された方法によって抑制されるため、領域２１０にリアルで連続的なＥＰ興奮パターンを示す。予想されるように、ＥＰ伝導領域２１０と非伝導瘢痕領域２２０との間は、明らかに不連続的である。
【０１７７】
図５の挿入図１５０は、参照のために開示された技術を用いずに計算された、左心房上のＥＰ興奮の色表現の補間のボリュームレンダリングである。挿入図１５０に示されている例では、波のリエントリーは考慮されていない。その結果、挿入図に見られるように、「レインボー」アーチファクト１５５が現れる。挿入図に示されている例では、心腔内のＥＰ波の伝播は、サイクル内の波面の「終わり近く（late）」が次のサイクルの「最初の方の（early）」波面と交わるリエントリーサイクル全体にわたっている。この例では、組織領域に、１サイクル長離れたＬＡＴが割り当てられる。虹（レインボー）状の狭いマルチシェードパスは、色分けされると、隣り合う領域間で１サイクル長の遷移を示すため、上記のように、本発明の実施形態が排除するマップの不連続性を引き起こす。即ち、上記のレインボーアーチファクトは、開示されたＩＣＥＡ方法を使用して完全に解消される。
【０１７８】
いくつかの実施形態では、周期的なカラースケールを使用して、対象ウィンドウ（ＷＯＩ）内の「Ｅａｒｌｙ  Ｍｅｅｔｓ  Ｌａｔｅ」への依存を解消する。色は、参照電位図に関連する相対的な時間差を表し、「Ｅａｒｌｙ  Ｍｅｅｔｓ  Ｌａｔｅ」という概念を排除する。カラーサークルは、興奮パターンを表す。したがって、例えば、標準のマップの場合とは異なり、赤色が早いこと（early）や紫色が遅いこと（late）を表すのではない。カラーサークルは、虹（レインボー）状の遷移を生じない自然な周期的な方法で、興奮時間の比例関係を表す。
【０１７９】
図６は、本発明の一実施形態による、ＥＰ興奮波の伝播を示す伝導矢印３３０がオーバーレイされた左心房のコヒーレントＥＰ興奮マップ４４０のボリュームレンダリングの概略的な絵画図である。伝導矢印３３０は、すべて固定長のベクトルであり、それぞれがＥＰ活動の更なる視覚化を提供する形状上の位置においてそれぞれのスローネスベクトルの方向を有する。例えば、伝導矢印３３０は、典型的には瘢痕領域を取り囲むように表示されている。別の一実施形態（図６には示されていない）では、伝導矢印は、その方向に加えて、スローネスのサイズを表す異なる長さを有する。
【０１８０】
上に述べた実施形態は例として挙げたものであり、本発明は上記に具体的に示し説明したものに限定されない点が理解されよう。むしろ、本発明の範囲は、上述の様々な特徴の組み合わせ及びその一部の組み合わせの両方、並びに上述の説明を読むことで当業者により想到されるであろう、また従来技術において開示されていない、それらの変形及び修正を含むものである。参照により本特許出願に援用される文献は、これらの援用文献において、いずれかの用語が本明細書において明示的又は暗示的になされた定義と矛盾して定義されている場合には、本明細書における定義のみを考慮するものとする点を除き、本出願の一部とみなすものとする。
【０１８１】
〔実施の態様〕
（１）  心腔の入力メッシュの表現、前記心腔の壁組織上の測定位置の組及び前記位置で測定された局所興奮時間（ＬＡＴ）のそれぞれの組を受信することと、
前記入力メッシュを、正則化された多角形を含む規則メッシュに再メッシュ化することと、
前記測定位置の組及びそれぞれのＬＡＴを前記正則化された多角形にデータ適合させることと、
前記正則化された多角形について、（ｉ）それぞれのＬＡＴ値及び（ｉｉ）前記壁組織が瘢痕組織を含むそれぞれの確率を反復的に計算して、瘢痕組織を示す前記規則メッシュ上での電気生理学的（ＥＰ）興奮波を求めることと、
前記規則メッシュ上にオーバーレイされた、前記ＥＰ興奮波及び前記瘢痕組織を含む電気解剖学的マップを表すことと、
を含む、方法。
（２）  前記正則化された多角形は、正則化された三角形を含む、実施態様１に記載の方法。
（３）  前記ＬＡＴ値及び前記確率を反復的に計算することは、前記ＬＡＴ値、スローネス値及び前記ＥＰ波の伝播についての方程式の組を反復的に解くことを含む、実施態様１に記載の方法。
（４）  前記確率を反復的に計算することは、瘢痕重み乗数を使用して再計算されたスローネスベクトルを、各正則化された多角形に再割り当てすることを含む、実施態様１に記載の方法。
（５）  前記測定位置のうちの少なくともいくつかにバイナリタグを付与して、瘢痕の存在及び二重ＥＰ電位の存在の少なくとも一方を示すことを含む、実施態様１に記載の方法。
【０１８２】
（６）  前記電気解剖学的マップを表すことは、前記電気解剖学的マップ上に瘢痕組織を示す伝導矢印をオーバーレイすることを含む、実施態様１に記載の方法。
（７）  心腔の入力メッシュの表現、前記心腔の壁組織上の測定位置の組及び前記位置で測定された局所興奮時間（ＬＡＴ）のそれぞれの組を受信するように構成されたインターフェースと、
プロセッサであって、
前記入力メッシュを、正則化された多角形を含む規則メッシュに再メッシュ化し、
前記測定位置の組及びそれぞれのＬＡＴを前記正則化された多角形にデータ適合させ、
前記正則化された多角形について、（ｉ）それぞれのＬＡＴ値及び（ｉｉ）前記壁組織が瘢痕組織を含むそれぞれの確率を反復的に計算して、瘢痕組織を示す前記規則メッシュ上での電気生理学的（ＥＰ）興奮波を求め、
前記規則メッシュ上にオーバーレイされた、前記ＥＰ興奮波及び前記瘢痕組織を含む電気解剖学的マップを表す、
ように構成されている、プロセッサと、
を備える、システム。
（８）  前記正則化された多角形は、正則化された三角形を含む、実施態様７に記載のシステム。
（９）  前記プロセッサは、前記ＬＡＴ値、スローネス値及び前記ＥＰ波の伝播についての方程式の組を反復的に解くことによって、前記ＬＡＴ値及び前記確率を反復的に計算するように構成されている、実施態様７に記載のシステム。
（１０）  前記プロセッサは、瘢痕重み乗数を使用して再計算されたスローネスベクトルを、各正則化された多角形に再割り当てすることにより、前記確率を反復的に計算するように構成されている、実施態様７に記載のシステム。
【０１８３】
（１１）  前記インターフェースは、瘢痕の存在及び二重ＥＰ電位の存在の少なくとも一方を示す、前記測定位置の少なくともいくつかのバイナリタグ付けを受信するように構成されている、実施態様７に記載のシステム。
（１２）  前記プロセッサは、前記電気解剖学的マップ上に瘢痕組織を示す伝導矢印をオーバーレイすることによって、前記電気解剖学的マップを表すように構成されている、実施態様７に記載のシステム。
【図１】
2020075132.tif 000058
【図２】
2020075132.tif 000059
【図３】
2020075132.tif 000060
【図４】
2020075132.tif 000061
【図５】
2020075132.tif 000062
【図６】
2020075132.tif 000063
【外国語明細書】
2020075132000001.pdf