JP2022180333A - リアルタイムでの局所的伝播速度の計算 - Google Patents

Abstract

【課題】心臓マッピングを行うこと。【解決手段】方法は、心臓の解剖学的表面上の複数の電極によって獲得されたそれぞれの信号に基づいて、電極のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算することを含む。この方法は、ＬＡＴに基づいて、これらの位置における電気伝播のそれぞれの方向を計算することを更に含む。この方法は、電極のうちの隣接するもののペアを、ペアの各々について、ペアを結合するベクトルが、所定の閾値整列度内で、ペアに属する電極のうちの１つの位置における電気伝播の方向と整列されるように選択することを更に含む。この方法は、電極のペアによって測定されたそれぞれのバイポーラ電圧を解剖学的表面のデジタルモデルと関連付けることを更に含む。他の実施形態も、記載されている。【選択図】図１

Description

（関連出願の相互参照）
本出願は、その開示が参照により本明細書に組み込まれる、２０２１年５月２４日出願の「Ｃｏｍｐｕｔｉｎｇ ｌｏｃａｌ ｐｒｏｐａｇａｔｉｏｎ ｖｅｌｏｃｉｔｉｅｓ ｆｏｒ ｃａｒｄｉａｃ ｍａｐｓ」と題する米国仮出願第６３／１９２，２２１号、及びその開示が参照により本明細書に組み込まれる、２０２１年５月２４日出願の「Ｃｏｍｐｕｔｉｎｇ ｌｏｃａｌ ｐｒｏｐａｇａｔｉｏｎ ｖｅｌｏｃｉｔｉｅｓ ｉｎ ｒｅａｌ－ｔｉｍｅ」と題する米国仮出願第６３／１９２，２３１号の利益を主張するものである。

（発明の分野）
本開示は、心臓マッピングの分野に関する。

心臓組織の任意の部分における局所興奮時間（ＬＡＴ）は、（ｉ）任意の心周期中に組織が電気興奮状態になる時間と、（ｉｉ）同じ周期中の基準時間との差である。基準時間は、例えば、身体表面心電図（ＥＣＧ）記録のＱＲＳ群の一点に設定され得る。

米国特許出願公開第２０１５／０１９６７７０号は、除細動ショックを送達するための手段と、患者の現在の心臓活動パラメータを連続的に収集するための手段と、少なくとも２つの層を有するニューラルネットワークを含むニューロン分析を用いた評価手段と、を含むアクティブ医療デバイスについて記載したものである。このニューラルネットワークは、不整脈発生因子のクラスに対応する別個のサブグループに分割されたそれぞれのパラメータを受け取る上流の３つのサブネットワークと、それら３つのサブネットワークに結合され、心室性不整脈のリスクの指標を出力することが可能な下流の出力ニューロンと、を含む。リスクの指標を所与の閾値と比較して、閾値と交差した場合、デバイスの少なくとも１つの機能を有効化又は無効化する。

米国特許出願公開第２０１０／０２６８０５９号は、患者の心臓の静脈網内の様々な位置に配置されたカテーテルを介して獲得された心臓情報にアクセスすることであって、この心臓情報は、位置情報と、電気的情報と、機械的情報とを含む、アクセスすることと、局所電気的活性化時間を解剖学的位置にマッピングして電気活性化時間マップを生成することと、局所機械的活性化時間を解剖学的位置にマッピングして機械的活性化時間マップを生成することと、局所電気的活性化時間を対応する局所機械的活性化時間から減算することによって、電気機械的遅延マップを生成することと、少なくとも電気機械的遅延マップをディスプレイにレンダリングすることと、を含む方法について記載している。

Ｃａｎｔｗｅｌｌ、Ｃｈｒｉｓ Ｄ．らの「Ｔｅｃｈｎｉｑｕｅｓ ｆｏｒ ａｕｔｏｍａｔｅｄ ｌｏｃａｌ ａｃｔｉｖａｔｉｏｎ ｔｉｍｅ ａｎｎｏｔａｔｉｏｎ ａｎｄ ｃｏｎｄｕｃｔｉｏｎ ｖｅｌｏｃｉｔｙ ｅｓｔｉｍａｔｉｏｎ ｉｎ ｃａｒｄｉａｃ ｍａｐｐｉｎｇ」（Ｃｏｍｐｕｔｅｒ ｉｎ ｂｉｏｌｏｇｙ ａｎｄ ｍｅｄｉｃｉｎｅ ６５（２０１５）：２２９～２４２）は、局所興奮時間を特定し、伝導方向及び速度を計算するように設計されたアルゴリズムについて概説している。

Ｒｏｎｅｙ、Ｃａｒｏｌｉｎｅ Ｈ．らの「Ａｎ ａｕｔｏｍａｔｅｄ ａｌｇｏｒｉｔｈｍ ｆｏｒ ｄｅｔｅｒｍｉｎｉｎｇ ｃｏｎｄｕｃｔｉｏｎ ｖｅｌｏｃｉｔｙ，ｗａｖｅｆｒｏｎｔ ｄｉｒｅｃｔｉｏｎ ａｎｄ ｏｒｉｇｉｎ ｏｆ ｆｏｃａｌ ｃａｒｄｉａｃ ａｒｒｈｙｔｈｍｉａｓ ｕｓｉｎｇ ａ ｍｕｌｔｉｐｏｌａｒ ｃａｔｈｅｔｅｒ」（２０１４ ３６ｔｈ Ａｎｎｕａｌ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｆ ｔｈｅ ＩＥＥＥ Ｅｎｇｉｎｅｅｒｉｎｇ ｉｎ Ｍｅｄｉｃｉｎｅ ａｎｄ Ｂｉｏｌｏｇｙ Ｓｏｃｉｅｔｙ，ＩＥＥＥ、２０１４）は、任意の配列の電極を用いて多極カテーテルデータから伝導速度を特定する一方で、波面の方向及び震源域の位置の推定を行う自動アルゴリズムについて記載している。

以下の本開示の実施例の詳細な説明を図面と併せ読むことで、本開示のより完全な理解が得られるであろう。
本開示のいくつかの例による、電気解剖学的マッピングのためのシステムの概略図である。 本開示のいくつかの例による、伝播速度を計算し、表示するためのアルゴリズムのフロー図である。 本開示のいくつかの例による、伝播速度計算の概略図である。 本開示のいくつかの例による、図２に示された選択ステップのフロー図である。 本開示のいくつかの例による、デジタルモデルの表面の概略図である。 本開示のいくつかの例による、図２に示された選択ステップの一例のフロー図である。 本開示のいくつかの例による、図６に示された補間ステップのフロー図である。 本開示のいくつかの例による、図７に示されたクラスタ化ステップのフロー図である。 本開示のいくつかの例による、表示されたモデルの概略図である。 本開示のいくつかの例による、表示されたモデルの概略図である。 本開示のいくつかの例による、反復凝縮アルゴリズムのフロー図である。 本開示のいくつかの例による、プローブ及び解剖学的表面の概略図である。 本開示のいくつかの例による、伝播速度のリアルタイムの視覚的表示の概略図である。 本開示のいくつかの例による、バイポーラ電圧を選択するための方法の概略図である。 本開示のいくつかの例による、バイポーラ電圧を選択するための方法の概略図である。 本開示のいくつかの例による、バイポーラ電圧測定のための電極のペアを選択するためのアルゴリズムのフロー図である。 本開示のいくつかの例による、電極の各位置におけるそれぞれのＬＡＴを計算するためのアルゴリズムのフロー図である。

概論
電気解剖学的マッピング中、心臓の表面に沿って、遠位端に複数の電極を備える体内プローブが移動される。電極によってその表面から獲得された生体電気信号に基づいて、表面上の様々な位置におけるＬＡＴなど、表面の様々な電気的特性が推定される。（１つ以上の獲得された信号に基づいてＬＡＴを推定するプロセスは、以下、ＬＡＴの「測定」とも呼ばれる。）

本開示の例は、ＬＡＴに基づいて、小さな空間分解能で伝播速度を正確に計算するためのアルゴリズムを提供する。本開示の例は、適切な診断及び治療を容易にするように、伝播速度を医師に視覚的に示すための技術を更に提供する。

特に、速度が計算される表面上の各「サンプリング位置」について、コンピュータプロセッサは、それぞれのＬＡＴが測定された近くにある「測定位置」の好適なセットを選択する。有利には、このセットは、電気的に不活性な組織によってサンプリング位置から分離された、いかなる測定位置も除外する。その後、セット内の測定位置（また、ＬＡＴがサンプリング位置において測定された場合は、サンプリング位置自体）の各々について、それぞれの４次元ベクトルが構築される。各ベクトルは、測定位置において測定されたＬＡＴから導出されたＬＡＴ値と共に、測定位置の位置座標から導出された３つの位置値を含む。

その後、プロセッサは、ベクトルに対して４×４の共分散行列の主成分分析（Principal Component Analysis、ＰＣＡ）を実行し、そのＰＣＡに基づいてサンプリング位置における伝播速度を計算する。例えば、プロセッサは、共分散行列の主成分を３つの位置次元に射影することによって、伝播の方向を計算し得る。次いで、プロセッサは、セット内の測定位置を、サンプリング位置を通過し、伝播方向に配向された線へと射影し得る。次に、プロセッサは、射影と対応するＬＡＴとの間の関係に近似する回帰関数を計算し得る。最後に、プロセッサは、この関数の勾配として速度の大きさ（すなわち、伝播の速度）を推定し得る。

いくつかの例では、伝播速度は、それぞれの測定位置において多数のＬＡＴを測定するように心臓表面を横切ってプローブが移動される、心臓表面の完全なマッピングに続いて計算される。具体的には、完全なマッピングに続いて、プロセッサは、表面のモデルを構築するが、モデルの表面上のそれぞれの測定点は測定位置に対応する。次に、プロセッサは、モデル表面上で、解剖学的表面上のそれぞれのサンプリング位置に対応する複数のサンプリング点を指定する。次いで、プロセッサは、測定位置の好適なセットに基づいて、各サンプリング位置における伝播速度を計算する。

そのような例では、伝播速度を計算することに続いて、プロセッサは、典型的には、伝播速度の１つ以上の特性を示すマーカーが重ねられた状態でモデルを表示する。例えば、各サンプリング点に、プロセッサは、サンプリング点における伝播速度の方向に配向されたマーカーを配置し得る。伝播速度の大きさが所定の閾値を下回るサンプリング点は、医師が緩徐伝導の領域を容易に特定し得るように、他のサンプリング点とは異なるようにマークされてもよい。

（本明細書におけるそのような例の説明では、モデルの表面上の点についての言及は、その点に対応する解剖学的表面上の位置についての言及に置き換えられてもよく、逆もまた同様であることに留意されたい。例えば、特定の測定位置において測定されたＬＡＴは、その測定位置に対応するモデル表面上の測定点と関連付けられているということができる。同様に、解剖学的表面上の位置の（ｘ，ｙ，ｚ）位置座標は、この位置に対応するモデル表面上の点の位置座標と呼ばれ得る。）

代替的にあるいは追加的に、伝播速度は、マッピング手順中にリアルタイムで計算され得る。特に、ＬＡＴ測定の各ラウンド（典型的には心周期ごとに１回行われる）後、プロセッサは、測定位置（すなわち、電極の位置）の全体にわたって反復し得る。各測定位置について、プロセッサは、近接する測定位置の好適なセットを特定することができる。次いで、プロセッサは、測定位置のそれぞれのベクトルを構築し、対応する共分散行列のＰＣＡを実行し、そのＰＣＡに基づいて伝播速度を計算することができる。

そのような例では、典型的には、プロセッサは、伝播速度の１つ以上の特性を示すマーカーが重ねられた状態で、電極を備えるプローブの遠位端のアイコンの表示を繰り返しリフレッシュする。例えば、各電極に、プロセッサは、電極における伝播速度の方向に配向されたマーカーを配置し得る。非リアルタイム表示の場合と同様に、マーカーの特性は、伝播速度の関数として変化し得る。

電気解剖学的マッピングを実行するときの別の課題は、バイポーラ電圧を測定するために使用される電極のペアが局所的な伝播方向に対して互いに垂直に配向されている場合に、測定されたバイポーラ電圧が心臓表面上の電気的に不活性な組織を誤って示し得ることである。

この課題に対処するために、本開示の例は、前述のリアルタイム計算を使用して、局所的な伝播方向と最も緊密に整列されている電極のペアを選定する。選定された電極のペア間のバイポーラ電圧は、心臓表面のモデルと関連付けられるが、他のバイポーラ電圧は、モデルから省略される。

本開示の例は、多数のバイポーラ電圧に基づいて、強化されたＬＡＴ計算を更に提供する。

システムの説明
最初に、本開示のいくつかの例による、電気解剖学的マッピングのためのシステム２０の概略図である図１を参照する。

図１では、心室の心内膜表面などの被験者２２の心臓２４の一部分の解剖学的表面に沿って、プローブ２６の遠位端を移動させている医師３０が示されている。プローブ２６の遠位端は、表面に沿って移動されているが、システム２０に属するプロセッサ３２は、プローブの遠位端にある電極２８を使用して、表面上の様々な測定位置におけるそれぞれの局所興奮時間（ＬＡＴ）を測定する。特に、電極２８が測定位置で電位図信号を獲得すると、プロセッサは、これらの信号を処理してＬＡＴを計算する。典型的には、電位図信号は、両方のユニポーラ信号、すなわち、電極と共通基準電極との間の信号と、バイポーラ電圧、すなわち隣接する電極のペアの間の電圧とを含む。

（典型的には、周期性不整脈の発生中にマッピングが実行される場合、ＬＡＴを計算することは、２つのステップを含む。第一に、背景技術において上述したように、標準的なＬＡＴ計算が実行される。その後、ＬＡＴの絶対値がｐ^＊ＣＬよりも大きい場合（ＣＬは周期性不整脈の各周期の長さであり、ｐ≧０．５が医師によって指定される）、ＬＡＴの絶対値がｐ^＊ＣＬ未満となるように、ＣＬをＬＡＴに加算する、又はＣＬをＬＡＴから減算する。）

いくつかの例では、プローブ２６の遠位端は、複数の平行なスプライン２９を含み、各スプライン２９は、電極の線形配列を備える。代替的に、電極のグリッドが、バルーン、拡張可能なプリント回路基板（ＰＣＢ）、又はプローブの遠位端にある任意の他の適切な構造上に配列され得る。

典型的には、プロセッサは、ポート又はソケットなどの電気的インターフェース３４を備えるコンソール４０内に収容される。プローブは、電気的インターフェース３４を介してコンソール４０に接続され、それにより、電極によって獲得された電位図信号は、電気的インターフェース３４を介してプロセッサによって受信される。典型的には、信号は、ワイヤに沿って、アナログ形式でプローブを介して搬送され、コンソールは、プロセッサ３２による処理のために信号をデジタル形式に変換するように構成されたアナログデジタル（Ａ／Ｄ）変換回路を更に備える。

マッピング手順中、プロセッサは、プローブの遠位端の位置を追跡する。追跡と電極から受信された電位図信号とに基づいて、プロセッサは、本明細書では「マップ」とも呼ばれる、心臓のこの部分のデジタルモデル３８を構築することができる。プロセッサは更に、任意の好適なタイプの揮発性メモリ若しくは不揮発性メモリを備えるメモリ３３にモデルを記憶すること、及び／又はモデル３８をディスプレイ３６上に表示することができる。

いくつかの例では、前述の追跡を容易にするために、プローブの遠位端は、１つ以上の電磁センサを備える。発生した磁場の存在下では、これらのセンサは、（例えば、電気的インターフェース３４を介して）プロセッサに、センサのそれぞれの位置を示す信号を出力する。そのような位置追跡技術は、例えば、参照によってそれぞれの開示が本明細書に組み込まれる、Ｂｅｎ－Ｈａｉｍの米国特許第５，３９１，１９９号、同第５，４４３，４８９号、及び同第６，７８８，９６７号、Ｂｅｎ－Ｈａｉｍらの米国特許第６，６９０，９６３号、Ａｃｋｅｒらの米国特許第５，５５８，０９１号、並びにＧｏｖａｒの米国特許第６，１７７，７９２号に開示されている。

他の例では、インピーダンスは、電極２８（及び／又はプローブの遠位端にある他の電極）と、被験者２２の身体に結合された電極パッチとの間で測定される。インピーダンス測定に基づいて、プロセッサは、電極の位置を確認する。典型的には、そのような例では、プロセッサは、例えば、それぞれの開示が参照により本明細書に組み込まれる、Ｇｏｖａｒｉらの米国特許第７，５３６，２１８号、及びＢａｒ－Ｔａｌらの米国特許第８，４５６，１８２号に記載されているように、事前に電磁センサを使用して較正された位置マップを利用する。

更に他の例では、電極パッチ間に電流が流れる。電極２８において測定された電圧に基づいて、プロセッサは、電極の位置を確認する。そのような技術は、例えば、それぞれの開示が参照により本明細書に組み込まれる、Ｗｉｔｔｋａｍｐｆへの米国特許第５，９８３，１２６号、Ｓｗａｎｓｏｎへの米国特許第６，４５６，８６４号、及びＮａｒｄｅｌｌａへの米国特許第５，９４４，０２２号に記載されている。

システム２０は、キーボード、マウス、又はディスプレイ３６に属するタッチスクリーンなどの１つ以上の入力デバイスを更に備え得る。医師３０は、入力デバイスを使用して、以下に記載される様々な閾値のうちのいずれかなどの任意の好適な入力を入力することができる。

概して、プロセッサ３２は、単一のプロセッサとして、又は協働的にネットワーク化若しくはクラスタ化されたプロセッサのセットとして具現化され得る。プロセッサ３２のいくつかの機能は、例えば、１つ以上の固定機能若しくは汎用集積回路、特定用途向け集積回路（ＡＳＩＣ）、及び／又はフィールドプログラマブルゲートアレイ（ＦＰＧＡ）を使用して、ハードウェアにのみ実装されてもよい。代替的に、この機能は、少なくとも部分的にソフトウェアで実装されてもよい。例えば、プロセッサ３２は、例えば、中央処理装置（ＣＰＵ）及び／又はグラフィック処理ユニット（ＧＰＵ）を備える、プログラムされたプロセッサとして具現化され得る。ソフトウェアプログラムを含むプログラムコード、及び／又はデータは、ＣＰＵ及び／又はＧＰＵによる実行並びに処理のためにＲＡＭにアップロードしてもよい。プログラムコード及び／又はデータは、例えば、ネットワークを介して、電子形態でプロセッサにダウンロード可能である。代替的に又は追加的に、プログラムコード及び／又はデータは、磁気記憶装置、光学的記憶装置、又は電子的記憶装置などの非一時的有形媒体上に提供及び／又は記憶され得る。このようなプログラムコード及び／又はデータは、プロセッサに提供されると、本明細書に記載されているタスクを行うように構成された機械又は専用コンピュータを作り出す。

伝播速度の計算
次に、本発明のいくつかの例による、伝播速度を計算し、表示するためのアルゴリズム５３のフロー図である図２を参照する。プロセッサ３２（図１）は、電位図信号が獲得されている間に、あるいはその獲得に続いて、アルゴリズム５３をリアルタイムで実行することができる。

アルゴリズム５３は、プロセッサが、心臓２４の解剖学的表面上の異なるそれぞれの測定位置において複数のＬＡＴを取得する（図１）、ＬＡＴ取得ステップ５５から始まる。例えば、プロセッサは、図１を参照して上述されたように、ＬＡＴを計算することによって、あるいは、マッピングの後にアルゴリズムを実行するときに、メモリ３３（図１）から、若しくはフラッシュドライブなどの外部記憶デバイスからＬＡＴを読み取ることによって、ＬＡＴを取得することができる。

ＬＡＴを取得することに続いて、プロセッサは、評価ステップ５７において、少なくとも１つのサンプリング位置について、伝播速度の計算を試行すべきかどうかを評価する。「はい」の場合、プロセッサは、選択ステップ５９において、ＬＡＴが計算に使用され得る測定位置のサブセットと共に、速度計算のためのサンプリング位置を選択する。（選択ステップ５９について、図４を参照して以下に更に説明する。）次に、サブセットサイズ評価ステップ６３において、プロセッサは、測定位置のサブセットが、計算するのに十分に大きいかどうか、すなわち、サブセットが閾値数の測定位置を含むかどうかを評価する。

測定位置のサブセットが十分に大きいという条件で、プロセッサは、測定位置のサブセットに対応するベクトルのセットを構築する。ベクトルの各々は、サブセット内の異なるそれぞれの測定位置について、測定位置のそれぞれの位置座標から導出された３つの位置値、及び測定位置において測定されたＬＡＴから導出されたＬＡＴ値を含む。（したがって、ベクトルの各々は、４次元である。）代替的に、サブセットが十分に大きくない場合、プロセッサは評価ステップ５７に戻る。

典型的には、ベクトルのセットを構築するために、プロセッサは最初に、スケーリング係数計算ステップ６５において、測定位置のサブセット全体にわたるＬＡＴの分散に基づいて、スケーリング係数を計算する。続いて、サブセット内の測定位置の各々について、プロセッサは、スケーリング係数によって、測定位置の位置座標又は測定位置におけるＬＡＴのいずれかをスケーリングし、スケーリングされたパラメータから測定位置に対応するベクトルを構築する。（測定位置の各々について、同じパラメータがスケーリングされる。）

例えば、プロセッサは最初に、ＬＡＴスケーリングステップ６７において、スケーリング係数によってＬＡＴをスケーリングすることができる。続いて、プロセッサは、ベクトル構築ステップ６９において、位置座標及びスケーリングされたＬＡＴからベクトルを構築することができる。例えば、任意の特定の測定位置について、プロセッサは、ベクトル［ｘ０ ｙ０ ｚ０ ｓ^＊Ｌ０］を構築することができ、（ｘ０，ｙ０，ｚ０）は、測定位置の位置座標であり、Ｌ０は測定位置におけるＬＡＴであり、ｓはスケーリング係数である。

ＬＡＴがスケーリングされる例では、スケーリング係数は、典型的には、測定位置のサブセット全体にわたる標的ＬＡＴ分散

とＬＡＴの実際の分散

との比であり、この比は、典型的には１より大きい。位置座標が代わりにスケーリングされる例では、スケーリング係数は、典型的には、

であり、この比は、典型的には１未満である。

いくつかの例では、目標ＬＡＴ分散は、図３を参照して以下に記載される所定の距離Ｄ０の所定の（増加）関数に従って計算される。例えば、目標ＬＡＴ分散は、ｃ^＊Ｄ０^２に等しくてもよく、式中、ｃは、例えば、５～１０の定数である。他の例では、測定位置のサブセットの空間的分散は、ｘ軸、ｙ軸、及びｚ軸の３つに対して計算され、目標ＬＡＴ分散は、最大の空間的分散の倍数として計算される。

ベクトルのセットを構築することに続いて、プロセッサは、ベクトルのセットの４×４の共分散行列の主成分分析（ＰＣＡ）に基づいて、選択されたサンプリング位置における電気伝播の方向を計算する。例えば、プロセッサは、共分散行列の第１の主成分（４次元ベクトル）を位置座標のそれぞれの次元へと射影することができる（したがって、３次元ベクトルをもたらす）。

実例として、以下の表１は、実験データから構築された８つのベクトルのセットを示す。（この特定の例では、ＬＡＴは、６倍にスケールアップされた。）

以下の表２は、このベクトルのセットの共分散行列を示す。

この行列の第１の主成分は、［－０．１３８ －１．１６５^＊１０^－１６ －０．０１６ －０．９９０］であり、ＸＹＺ空間上への第１の主成分の射影は［－０．１３８ －１．１６５^＊１０^－１６ －０．０１６］であり、これは、単位方向ベクトルとして表現すると、［０．９９３ ８．４０６^＊１０^－１６０．１１９］である。

伝播方向の計算に続いて、プロセッサは、速度計算ステップ７３において、選択されたサンプリング位置における伝播速度を計算する。（速度計算ステップ７３について、図３を参照して以下に説明する。）プロセッサは次いで、評価ステップ５７に戻る。

評価ステップ５７において、更に計算されるべき速度がないことを確認することに加えて、プロセッサは、表示ステップ７５において、図９Ａ～図９Ｂ及び図１２を参照して以下で更に説明するように、ディスプレイ３６（図１）上に速度を示す。

代替例では、プロセッサは、速度を計算せずに、電気伝播の方向を計算し、速度を示さずに、ディスプレイ上に方向を示す。

次に、アルゴリズム５３のいくつかの例に関する更なる詳細について、本開示のいくつかの例による伝播速度計算の概略図である図３を参照する。

図３は、解剖学的表面４２上の複数の測定位置４４をしている。ＬＡＴ取得ステップ５５を実行することによって、測定位置４４におけるＬＡＴを取得することに続いて、プロセッサは、上述したアルゴリズム５３の後続のステップを実行する。したがって、例えば、選択ステップ５９を実行する際に、プロセッサは、最初に解剖学的表面４２上のサンプリング位置４５を選択することができる。（図５を参照して以下で更に説明するように、サンプリング位置４５は、測定位置４４のうちの１つと必ずしも一致しない。）続いて、プロセッサは、図６を参照して以下で更に説明するように、測定位置４４のうちの、サンプリング位置４５から所定の距離Ｄ０内にあるものを特定し、次いで、特定された測定位置から測定位置のサブセットを選択することができる。例えば、図３の例によれば、プロセッサは、サンプリング位置４５から距離Ｄ０内にある測定位置のうち、測定位置４４ｄ及び４４ｅを選択せずに、測定位置４４ａ、４４ｂ、及び４４ｃを選択する。

図５～図１０を参照して以下で更に説明するように、アルゴリズム５３がリアルタイムで実行されない場合、Ｄ０は、典型的には６ｍｍ～１２ｍｍである。図１１～図１２を参照して以下で更に説明するように、アルゴリズム５３がリアルタイムで実行される場合、Ｄ０は、典型的には、５ｍｍ～７ｍｍなど、４ｍｍ～８ｍｍである。より一般的には、Ｄ０は、測定位置の密度の関数であり得る。

（図３では２次元に見えるが、解剖学的表面４２は３次元であることに留意されたい。したがって、典型的には、プロセッサは、表面に沿った距離を計算するために測地線距離測定を使用する。）
いくつかの例では、プロセッサが続いて、サブセットサイズ評価ステップ６３において、サブセットが小さすぎることを確認した場合、プロセッサは、毎回選択を繰り返して、距離Ｄ０を１回以上（所定の最大回数まで）増加させることができる。

サブセットが十分に大きいことを確認すると、プロセッサは、上述のように伝播方向を計算する。続いて、プロセッサは、速度計算ステップ７３において伝播速度を計算する。

いくつかの例では、伝播速度を計算するために、プロセッサは、最初に、測定位置のサブセットを仮想線４７へと射影するが、仮想線４７は、サンプリング位置４５を通過し、サンプリング位置における電気伝播の方向に配向されている。次いで、プロセッサは、射影が存在する線４７に沿ったそれぞれの距離を計算する。次に、プロセッサは、サブセット内の異なるそれぞれの測定位置について、（ｉ）測定位置の射影が存在する線に沿った距離、及び（ｉｉ）測定位置におけるＬＡＴ、を各々が含む点４９のグループを定義する。続いて、プロセッサは、点４９に（あるいは、サンプリング位置４５において測定された、又は、図５を参照して以下で説明されるようにサンプリング位置４５に対して補間されたＬＡＴ付近の点４９のサブセットに）回帰関数５１を適合し、次いで、関数５１の勾配として電気伝播速度を計算する。

例えば、関数５１は線であってもよく、電気伝播速度はその線の勾配として計算されてもよい。代替的に、関数は、２次以上の多項式、スプライン関数、又は任意の他の好適な種類の関数であり得る。そのような例では、速度は、サンプリング位置４５において測定された、又はサンプリング位置４５に対して補間されたＬＡＴにおける関数の勾配として計算され得る。代替的に、プロセッサは、任意の他の好適な方式で電気伝播速度を計算し得る。例えば、プロセッサは、点４９の２×２の共分散行列を計算し、次いで、この行列の第１の主成分に基づいて電気伝播速度を計算し得る。代替的に、プロセッサは、サンプリング位置４５の一方の側の第１のグループ、及びサンプリング位置４５の他方の側の第２のグループの２つのグループに射影を分割し得る。（したがって、図３の例によれば、第１のグループは、測定位置４４ａ及び４４ｂの射影を含むが、第２のグループは、測定位置４４ｃの射影を含む。）続いて、プロセッサは、各グループのそれぞれの重心を計算することができ、各重心は、グループ内の射影の平均位置である。次いで、２つの重心の間を通る線の勾配として速度を計算し得る。

選択ステップ５９に関する更なる詳細について、本開示のいくつかの例による、選択ステップ５９のフロー図である図４を同様に参照する。

選択ステップ５９は、プロセッサがサンプリング位置４５を選択する、サンプリング位置選択ステップ６０から始まる。続いて、プロセッサは、チェックステップ６２において、サンプリング位置から距離Ｄ０内の任意の測定位置がまだ処理されていないかどうかをチェックする。「はい」の場合、プロセッサは、測定位置特定ステップ７０において、処理のために、これらの測定位置のうちの１つを特定する。続いて、プロセッサは、差計算ステップ７４において、測定位置におけるＬＡＴとサンプリング位置におけるＬＡＴとの間の差を計算する。

計算されたＬＡＴ差に基づいて、プロセッサは、速度推定ステップ７６において、サンプリング位置と測定位置との間の伝播速度を推定する伝播速度推定値を計算する。典型的には、この推定値はＬＡＴ差と、サンプリング位置と測定位置との間の距離との商であり、これはプロセッサによってチェックステップ６２において計算される。

続いて、プロセッサは、チェックステップ７８において、伝播速度推定値が所定の速度推定閾値ｖ_Ｔ１を超えているかどうかをチェックする。「はい」の場合、サブセット増加ステップ８０において、測定位置の選択されたサブセットに測定位置を追加する。そうでない場合、サブセットに測定位置を追加しない。（伝播速度推定値がｖ_Ｔ１を超えない場合、サンプリング位置と測定位置との間にブロックがあり、それにより、サブセットに測定位置を追加すると、伝播速度の後続の計算が不正確になる可能性がある。）

サブセット増加ステップ８０に続いて、又は伝播速度推定値がｖ_Ｔ１を超えない場合、プロセッサは、チェックステップ６２に戻る。距離Ｄ０内の全ての測定位置が処理されたことを確認すると、選択ステップ５９が終了する。

モデルに基づく計算及び表示
次に、本開示のいくつかの例による、デジタルモデル３８の表面４２’の概略図である図５を参照する。

いくつかの例では、プロセッサは、解剖学的表面４２（図３）上におけるプローブ２６（図１）の遠位端の複数の位置に対応する点群から、典型的には点群のトライアングルテッセレーションを実行することによってモデル３８を構築する。測定位置４４（図３）は、モデル３８の表面４２’上の異なるそれぞれの測定点４４’に対応する。（表面４２’は、点群の全ての点を必ずしも通過しない「最良の適合」であるので、いくつかの測定点４４’は、点群から表面４２’へと点を射影することによって計算され得ることに留意されたい。）

そのような例では、伝播速度の計算は、モデルの構築後に実行され得る。具体的には、プロセッサは、例えば、表面を均一にサンプリングすることによって、表面４２’上で複数のサンプリング点４５’を指定することができる。続いて、プロセッサは、アルゴリズム５３を実行するときにサンプリング点全体にわたって反復することができる（図２）。具体的には、アルゴリズム５３の評価ステップ５７において、プロセッサは、任意のサンプリング点がまだ処理されていないかどうかを評価し得る。「はい」の場合、プロセッサは、以下でモデルに基づく選択ステップ５９’と呼ばれる選択ステップ５９の例を実行し得る。続いて、測定位置の選択したサブセットが十分に大きい場合、プロセッサは、図２～３を参照して上述したように、サンプリング点の伝播速度を計算することができる。

次に、本開示のいくつかの例による、モデルに基づく選択ステップ５９’のフロー図である図６を更に参照する。

モデルに基づく選択ステップ５９’は、プロセッサがサンプリング点を選択する、サンプリング点選択ステップ６０’から始まる。

いくつかの例では、サンプリング点選択ステップ６０’に続いて、プロセッサは、チェックステップ７１において、選択されたサンプリング点から距離Ｄ０内に任意の測定点があるかどうかをチェックする。距離Ｄ０内にない場合、モデルに基づく選択ステップ５９’が終了する。そうでない場合、プロセッサは、チェックステップ６１において、選択されたサンプリング点が測定点４４’のうちのいずれかの所定の距離Ｄ１内にあるかどうかをチェックする。典型的には、Ｄ１は、Ｄ０よりもはるかに小さく、例えば、０．５ｍｍ以下である。

サンプリング点が少なくとも１つの測定点のＤ１内にある場合、プロセッサは、このサンプリング点が最も近い測定点と一致すると見なす。（言い換えれば、プロセッサは、サンプリング点に対応するサンプリング位置が、ＬＡＴが取得された最も近い測定位置と一致すると見なす）。したがって、プロセッサは、ＬＡＴ割り当てステップ５８において、最も近い測定点のＬＡＴをサンプリング点に割り当てる。

代替的に、測定点がサンプリング点のＤ１内にない場合、プロセッサは、補間ステップ５６において、図７～図８を参照して以下で更に説明するように、測定点のうちの、チェックステップ７１において特定されたものに関連付けられたＬＡＴのうちの少なくともいくつかを補間することによって、サンプリング点の補間されたＬＡＴを計算し、その補間されたＬＡＴをサンプリング点に割り当てる。

したがって、例えば、図５に示すように、第１のサンプリング点４５’ａは、最も近い測定点４４’ａのＬＡＴを割り当てられ得るが、補間されたＬＡＴは、いかなる測定点にも十分に近くない第２のサンプリング点４５’ｂについて計算され得る。

サンプリング点にＬＡＴを割り当てることに続いて、プロセッサは、測定点特定ステップ７０’において、処理のために、未処理の測定点を特定する。続いて、プロセッサは、差計算ステップ７４’において、サンプリング点に関連付けられるＬＡＴと測定点に関連付けられるＬＡＴとの間の差を計算する。計算されたＬＡＴ差に基づいて、プロセッサは、速度推定ステップ７６’において、測定点とサンプリング点との間の伝播速度を推定する。続いて、プロセッサは、チェックステップ７８において、伝播速度推定値が所定の速度推定閾値ｖ_Ｔ１を超えているかどうかをチェックする。「はい」の場合、サブセット増加ステップ８０’において、測定位置の選択されたサブセットに測定点を追加する。そうでない場合、サブセットに測定点を追加しない。

サブセット増加ステップ８０に続いて、又は伝播速度推定値がｖ_Ｔ１を超えない場合、プロセッサは、チェックステップ６２’において、距離Ｄ０内に未処理の測定点が他にもあるかどうかをチェックする。距離Ｄ０内の全ての測定点が処理されたことを確認すると、選択ステップ５９’が終了する。

いくつかの例では、補間ステップ５６において補間されたＬＡＴを計算するために、プロセッサは最初に、クラスタの各々について、クラスタの各ペアについて、測定点のペアの間の伝播速度を推定する伝播速度推定値が、ｖ_Ｔ１とは異なり得る所定の速度推定閾値ｖ_Ｔ２を超えるように、測定点のうちの、サンプリング点から距離Ｄ０内にあるものを１つ以上のクラスタへとクラスタ化する。次いで、プロセッサは、サンプリング点とクラスタとの間のそれぞれの距離に基づいて、クラスタのうちの１つを特定する。次いで、プロセッサは、ＬＡＴのうちの、特定されたクラスタ内の測定点のうちの少なくともいくつかとそれぞれ関連付けられるＬＡＴの加重平均として、サンプリング点の補間されたＬＡＴを計算する。

これに関して、次に、本開示のいくつかの例による、補間ステップ５６のフロー図である図７を参照する。

補間ステップ５６は、プロセッサが、上述のように、測定点のうちの、サンプリング点から距離Ｄ０内にあるものをクラスタ化する、クラスタ化ステップ８４から始まる。クラスタ化ステップ８４の例について、図８を参照して以下に説明する。

クラスタ化ステップ８４に続いて、第１のクラスタ選択ステップ８６において、各クラスタが選択される。選択された各クラスタについて、プロセッサは、距離計算ステップ８８において、サンプリング点からクラスタまでの距離を計算する。一般に、この距離の任意の好適な定義が使用され得る。例えば、プロセッサは、この距離を、サンプリング点と、クラスタ内の、そのサンプリング点に最も近いＮ個の測定点との間の平均距離として計算してもよく、Ｎは、例えば、２、３、又は４である。

続いて、プロセッサは、チェックステップ９０において、選択されるべきクラスタが他にもあるかどうかをチェックする。「はい」の場合、プロセッサは、クラスタ選択ステップ８６に戻る。そうでない場合、クラスタ特定ステップ９２において、プロセッサは、サンプリング点から最小の距離を有するクラスタを特定する。クラスタを特定することに続いて、プロセッサは、加重平均化ステップ９４において、サンプリング点に最も近いＮ測定点など、特定されたクラスタ内の測定点のうちの少なくともいくつかと関連付けられたＬＡＴの加重平均として、補間されたＬＡＴを計算する。典型的には、各ｉ番目の測定点のＬＡＴに対する重みｗ_ｉは、

に等しく、式中、ｄ_ｉは、サンプリング点からのｉ番目の測定点の距離である。

上述のように補間を実行することにより、プロセッサは一般に、補間されたＬＡＴが、電気的に不活性な組織によってサンプリング点から分離された測定点に基づくことを控える。

次に、本開示のいくつかの例による、クラスタ化ステップ８４のフロー図である図８を参照する。

クラスタ化ステップ８４を実行するために、プロセッサは、第１の測定点選択ステップ９６において、サンプリング点から距離Ｄ０内の各測定点を反復的に選択する。図８で「ＭＰ１」と呼ばれる選択された各測定点について、プロセッサは、チェックステップ９８において、任意のクラスタが存在し、まだ選択されていないかどうかをチェックする。存在しない場合、プロセッサは、クラスタ初期化ステップ１１２において、ＭＰ１でクラスタを初期化する。そうでない場合、プロセッサは、第２のクラスタ選択ステップ１００で、まだ選択されていない次のクラスタを選択する。

クラスタを選択することに続いて、プロセッサは、第２の測定点選択ステップ１０２において、クラスタ内の、図８で「ＭＰ２」と呼ばれる測定点のうちの１つを選択する。次に、別の速度推定ステップ１０４において、プロセッサは、ＭＰ１とＭＰ２との間の伝播速度を推定する伝播速度推定値を計算する。プロセッサは次いで、チェックステップ１０６において、伝播速度推定値が閾値ｖ_Ｔ２を超えているかどうかをチェックする。超えない場合、プロセッサは、チェックステップ９８に戻る。そうでない場合、プロセッサは、チェックステップ１０８において、クラスタ内の任意の測定点がまだ選択されていないかどうかをチェックする。まだ選択されていない少なくとも１つの測定点が存在する場合、プロセッサは、第２の測定点選択ステップ１０２に戻り、次の測定点ＭＰ２を選択する。そうでない場合、プロセッサは、クラスタ成長ステップ１１０において、クラスタにＭＰ１を追加する。

クラスタ成長ステップ１１０又はクラスタ初期化ステップ１１２に続いて、プロセッサは、チェックステップ１１４において、サンプリング点から距離Ｄ０内の任意の測定点がまだ選択されていないかどうかをチェックする。「はい」の場合、プロセッサは、第１の測定点選択ステップ９６に戻り、次の測定点ＭＰ１を選択する。そうでない場合、クラスタ化ステップ８４が終了する。

次に、本開示のいくつかの例による、表示されたモデルの概略図である図９Ａを参照する。

伝播速度を計算することに続いて、プロセッサは、それぞれのマーカー１１６が、サンプリング点４５’においてモデルの上に重なり、電気伝播の方向に配向された状態で、モデル３８（特に、モデル表面４２’）をそれぞれ表示する。（より良好な視認性のために、各マーカー１１６は、例えば、モデル表面４２’に平行な方向に、そのサンプリング点からわずかにオフセットされ得る。）各マーカー１１６は、図９Ａに示される矢頭形状又は完全な矢印形状などの任意の好適な形状を有し得る。任意選択的に、図９Ａに示される円などの追加のマーカーにより、モデル表面４２’上にサンプリング点の位置をマークしてもよい。

電気伝播の方向を示すようにマーカー１１６を配向することに加えて、プロセッサは、電気伝播速度に従って、マーカーの色、形状、長さ、又は厚さなどの少なくとも１つの他の特性を変化させ得る。例えば、その特性は、所定の速度閾値ｖ_Ｔ３を超えない（又は以下に記載されるように、信頼度の閾値測定値ではｖ_Ｔ３を超えない）各速度の第１の値に、そそうでない場合には，第２の値に設定され得る。特定の例として、速度がｖ_Ｔ３を超えない各サンプリング点には、より太いマーカー１１６ａが配置され得るが、他のサンプリング点には、より細いマーカー１１６ｂが配置され得る。したがって、より太いマーカー１１６ａは、医師にとって関心があり得る解剖学的表面上の緩徐伝導の領域を示す。

いくつかの例では、速度閾値ｖ_Ｔ３を超えない各速度について、信頼度が計算される。信頼度が所定の信頼閾値を超える場合、マーカー特性は、第１の値に設定され、そうでない場合、特性は、第２の値に設定される。信頼度は、例えば、（ｉ）速度が速度閾値を超えない近接するサンプリング点の数と、（ｉｉ）近接するサンプリング点の総数との比として定義され得る。サンプリング点Ｓ１の近接するサンプリング点Ｓ２は、例えば、Ｓ１とＳ２との間の伝播速度推定値がｖ_Ｔ１を超えるという条件で、Ｓ１からの所定の距離内にある任意のサンプリング点として定義され得る。

いくつかの例では、モデルは、解剖学的表面の他の特性を更に示すように表示される。例えば、モデル表面４２’は、解剖学的表面上のＬＡＴを示すように着色され得る。

いくつかの例では、電気伝播の方向を視覚的に示す前に、プロセッサは、電気伝播の方向を平滑化する。例えば、プロセッサは、ラプラシアン平滑化を実行することができ、それにより、平滑化の各ｉ番目の反復中に、各サンプリング点における単位伝播方向ベクトルＶ［ｉ］は、α^＊Ｖ［ｉ－１］＋（１－α）^＊Ｖ_Ａ［ｉ］，として計算され、式中、Ｖ_Ａは、サンプリング点の近接点の平均単位伝播方向ベクトルである。（上述のように、近接点は、２つの点の間の伝播速度推定値がｖ_Ｔ１を超えるという条件で、サンプリング点からの所定の距離内にある任意の他のサンプリング点であってよい。）図９Ａは、その近接点とほぼ同じ方向に再配向されたマーカーのうちの１つを示すことによって、そのような平滑化操作の例示的な結果を示している。

次に、本開示のいくつかの例による、表示されたモデルの別の概略図である図９Ｂを参照する。

いくつかの例では、電気伝播の方向を平滑化することに代わって、あるいはそれに加えて、プロセッサは、マーカーが重ねられた状態でモデルを表示する前に、サンプリング点が１つ以上の平均伝播経路１２０にほぼ追従するように、サンプリング点を凝縮する。したがって、プロセッサは、医師による、表示されたモデルの解釈を容易にする。

典型的には、凝縮は、反復アルゴリズムを実行することによって実行される。アルゴリズムの各反復中、プロセッサは、電気伝播の方向を再計算し、次いで、電気伝播の方向に応じて、サンプリング点を互いに向かってシフトさせる。

これに関して、次に、本開示のいくつかの例による、そのような１つの反復凝縮アルゴリズムステップ１２１のフロー図であるここで図１０を参照する。

アルゴリズム１２１の各反復の開始時に、プロセッサは、再計算ステップ１２８において伝播方向を再計算する。典型的には、このステップにおいて、伝播速度の計算は、サンプリング点の近接点が測定点のサブセットに置き換わっている状態で、図２を参照して上述したように実行される。例えば、各サンプリング点について、プロセッサは、（ｉ）サンプリング点及びその近接点のそれぞれのベクトルを構築し、（ｉｉ）対応する４×４の共分散行列のＰＣＡを実行し、（ｉｉｉ）ＰＣＡに基づいて、例えば共分散行列の第１の主成分を位置座標のそれぞれの次元に射影することによって、伝播方向を計算することができる。

続いて、プロセッサは、チェックステップ１２３において、選択されるべきサンプリング点が他にもあるかどうかをチェックする。「はい」の場合、プロセッサは、サンプリング点選択ステップ１２５において、次のサンプリング点を選択する。次に、プロセッサは、平均位置計算ステップ１２２において、サンプリング点の近隣の点の平均位置を計算する。これらの点は、選択されたサンプリング点の近接点と共に、選択されたサンプリング点を含む。図９Ａを参照して上述したように、近接点は、２つのサンプリング点の間の伝播速度推定値が速度推定閾値ｖ_Ｔ１を超えるという条件で、選択されたサンプリング点からの所定の距離内にある任意の他のサンプリング点として定義されてもよい。

平均位置計算ステップ１２２に続いて、プロセッサは、射影計算ステップ１２４において、サンプリング点において伝播方向に配向され、平均位置を通過する線への、サンプリング点の射影を計算する。次いで、プロセッサは、点移動ステップ１２６において、サンプリング点を射影に向かって移動させる。例えば、アルゴリズムの各ｉ番目の反復中、プロセッサは、サンプリング点の新しい位置Ｐ［ｉ］をα^＊Ｐ［ｉ－１］＋（１－α）^＊Ｐ_Ｐ［ｉ］として計算することができ、式中、Ｐ_Ｐは、サンプリング点の射影である。サンプリング点を移動させることに続いて、プロセッサは、チェックステップ１２３に戻る。

チェックステップ１２３において、全てのサンプリング点が、この反復中に選択されたことを確認すると、プロセッサは、評価ステップ１３０において、別の反復を実行するかどうかを評価する。「はい」の場合、プロセッサは、アルゴリズムの次の反復を実行する。そうでない場合、アルゴリズムの実行が終了する。

一般に、評価ステップ１３０は、任意の好適な基準に基づき得る。例えば、プロセッサは、この反復中に、サンプリング点のいずれもが、所定の閾値距離を超えて移動していない場合、又は所定の最大反復数が実行された場合、アルゴリズムの実行を終了し得る。

リアルタイム計算及び表示
次に、本開示のいくつかの例による、プローブ２６及び解剖学的表面４２の概略図である図１１を参照する。

図１～図２を参照して上述したように、測定位置４４におけるＬＡＴは、プローブ２６に属する電極２８によって獲得された信号に基づいて、プロセッサによって計算される。いくつかの例では、プロセッサは、測定位置４４のうちの少なくともいくつかにおいてリアルタイムで、すなわち電極が測定位置にある間に、それぞれの伝播速度（又は少なくとも伝播方向）を更に計算する。（言い換えれば、リアルタイムでは、プロセッサは、測定位置４４のうちの少なくともいくつかをサンプリング位置４５として取り扱い、したがって、解剖学的表面のデジタルモデルを使用せずに、これらの測定位置における伝播速度を計算する。）例えば、伝播速度は、心周期ごとに１回計算され得る。このリアルタイム計算は、図２～図４を参照して上述したように、上述のモデルに基づく計算に代わって、あるいはそれに加えて実行され得る。

次に、本開示のいくつかの例による、伝播速度のリアルタイム視覚指示の概略図である図１２を参照する。

伝播速度の各リアルタイム計算に続いて、プロセッサは速度を示す。典型的には、電気伝播の方向を示すために、プロセッサは、プローブのアイコン２６’を表示し、電気伝播の方向に配向されたそれぞれのマーカー１１７を、電極のうちの、方向を計算されたサンプリング位置に位置する、ものに対応するアイコンの部分２８’に配置する。（簡単にするために、マーカー１１７は、少数の電極についてのみ図１２に示されている。）いくつかの例では、電気伝播の方向を示す前に、プロセッサは、図９Ａを参照して上述したように、その方向を平滑化する。

典型的には、プロセッサは、電気伝播速度に従って、マーカーの少なくとも１つの特性（例えば、色、形状、長さ、又は厚さ）を変化させる。例えば、マーカーは、（ｉ）第１の範囲（例えば、ｖよりも大きい速度の場合）に属する速度のもの（例えば、速度がｖ_Ｔ３よりも大きいもの）について、第１の形状及び第１の厚さ、（ｉｉ）第１の範囲よりも低い第２の範囲に属する速度のもの（例えば、ｖ_Ｔ３未満であるが、下限閾値ｖ_Ｔ４よりも大きいもの）について、第１の形状及び第２の厚さ、並びに（ｉｉｉ）第２の範囲よりも小さい第３の範囲に属する速度のもの（例えば、ｖ_Ｔ４未満の速度のもの）について、第２の形状を有し得る。図１２は、（ｉ）第１のマーカー１１７ａが、正常な伝播速度を示す、より長くかつより細い矢印を含み、（ｉｉ）第２のマーカー１１７ｂが、より遅い伝播速度を示すより短くかつより太い矢印を含み、（ｉｉｉ）第３のマーカー１１７ｃが、電気的に不活性な組織を示す円を含む、そのような例を示す。（典型的には、誤ったマーキングを回避するために、電気的に不活性な組織は、関連する電極が組織に接触していることが知られている場合にのみマークされる。）

任意選択的に、マーカーの別の特性は、ＬＡＴに応じて変更され得る。例えば、マーカーは、ＬＡＴに基づいて、カラースケールに従って着色され得る。

図１２に示すように、プローブのアイコンは、マッピングされている解剖学的表面の画像４２”の上に重ねられ得る。代替的に、アイコンは、モデルの表面上に重ねられてもよい。いくつかの例では、アイコン及びマーカーの表示は、心周期中のプローブの動きを考慮するために、心周期当たり複数回リフレッシュされる。

バイポーラ電圧の最適な選択
次に、本開示のいくつかの例による、バイポーラ電圧測定用の電極のペアを選択するための方法の概略図である図１３を参照する。

いくつかの例では、電極２８の位置における電気伝播のそれぞれの方向を計算することに続いて、プロセッサは、電極のうちの隣接するもののペア１３９を、各ペア１３９について、ペアを互いに結合するベクトル１４０が、所定の閾値整列度度内で、ペアに属する電極のうちの１つの位置における電気伝播の方向と整列されるように選択する。（したがって、ペア選択は、伝播速度ではなく、伝播方向に基づく。）電極のペアを選択することに続いて、プロセッサは、電極のペアによって測定されたそれぞれのバイポーラ電圧をモデル３８に関連付ける（図１）。したがって、有利には、関連性がより低いバイポーラ電圧は、モデルから省略される。

図１３に示すような電極の矩形グリッドでは、同じ行の隣接する（又は「近接する」）電極間の距離が、隣接する行間の距離と同じであり、閾値整列度は一般に４５°である。電極の各々（典型的には、コーナ電極のうちの１つ）について、プロセッサは、この電極を同じ行内の隣接する電極とペアリングするかどうか、この電極を隣接する行内の隣接する電極とペアリングするかどうか、又は電極を全くペアリングしないかどうかを決定する。特に、潜在的にペアリング可能な隣接電極の各々について、プロセッサは、電極から隣接電極（又はその逆）を指すベクトル１４０と、伝播方向に配向された別のベクトル１１８との間の角度θを計算する。θ（又は｜１８０°－θ｜）が閾値角度未満である場合、隣接電極は、その電極とペアリングされる。

（典型的には、同じ行の単一の隣接電極のみが潜在的にペアリング可能であり、この電極は、ペアが求められる電極の常に左側にあるか、又は常に右側にある。同様に、典型的には、隣接する行内の単一の隣接電極のみが潜在的にペアリング可能であり、隣接する行は、ペアが求められる電極の常に上にあるか、又は常に下にある。したがって、電極の対が複数回、選択されることはない。）

したがって、図１３に示される例では、第１の電極２８ａは、同じ列内の電極２８ａの下にある隣接電極とペアリングされており、第２の電極２８ｂは、第２の電極２８ｂの右側の列内の隣接電極とペアリングされている。

次に、本開示のいくつかの例による、電極の六角配列のための図１３の方法を示している図１４を更に参照する。

いくつかの例では、その開示が参照により本明細書に組み込まれる米国特許出願第１７／０９２，６２７号に記載されているように、電極２８は、各電極が最大６つの近接電極から等距離に離間されるように、六角格子で配列される。例えば、スプライン２９上の電極の行（図１３）は互い違いにされてもよい。そのような例では、閾値整列度は一般に３０°であり、プロセッサは、ペアリングのために最大３つの隣接電極を考慮する。

次に、更なる詳細について、本開示のいくつかの例による、バイポーラ電圧測定用の電極のペアを選択するためのアルゴリズム１４２のフロー図である図１５を参照する。

典型的には、アルゴリズム１４２は、心周期ごとに少なくとも１回、例えば厳密に１回実行される。アルゴリズム１４２ごとに、電極選択ステップ１４４において電極の各々（典型的には、コーナ電極のうちの潜在的にペアリング可能な近接電極を有さない１つを除く）が選択される。電極の選択に続いて、潜在的にペアリング可能な近接電極（例えば、選択された電極の右方近接電極又は下方近接電極）が、近接選択ステップ１４６で選択される。次に、プロセッサは、角度計算ステップ１４８で、電極のペア（すなわち、選択された電極とその選択された近接電極）の角度θを計算する（図１３～図１４）。

θを計算することに続いて、プロセッサは、角度比較ステップ１５０において、θ（又は｜１８０°－θ｜）が閾値角度（例えば、４５°又は３０°）未満であるかどうかを確認する。「はい」の場合、この電極のペアは、ペア選択ステップ１５２において、バイポーラ電圧測定用に選択される（すなわち、ペア間のバイポーラ電圧は、モデル３８（図５）との関連付け用に選択される）。そうでない場合、プロセッサは、チェックステップ１５４において、潜在的にペアリング可能な近接電極がまだ選択されていないかどうかをチェックする。「はい」の場合、プロセッサは、近接選択ステップ１４６に戻り、次の潜在的にペアリング可能な近接電極を選択する。

ペア選択ステップ１５２を実行するか、又は潜在的にペアリング可能な近接電極が選択されずに残存していることを確認することに続いて、プロセッサは、チェックステップ１５６で、選択されるべき電極が残存しているかどうかをチェックする。「はい」の場合、プロセッサは、電極選択ステップ１４４に戻り、次の電極を選択する。そうでない場合、プロセッサは、モデル増加ステップ１５８において、選択された電極ペアによって測定されたバイポーラ電圧をモデル３８と関連付ける（図５）。例えば、プロセッサは、バイポーラ電圧全体の範囲に及ぶカラースケールに従って、モデルの表面４２’を着色してもよい。代替的又は追加的に、表面４２’上の点の上でのマウスポインタの移動、又は点上でのマウスのクリックに応じて、プロセッサは、その点におけるバイポーラ電圧が獲得された電極ペアの指示、及び／又はバイポーラ電圧信号自体を表示し得る。

強化されたＬＡＴ計算
次に、本発明のいくつかの例による、電極の各位置におけるそれぞれのＬＡＴを計算するためのアルゴリズム１６０のフロー図である図１６を参照する。アルゴリズム１６０は、電気解剖学的マッピング手順中に任意の時点でプロセッサによって実行され得る。

導入として、アルゴリズム１６０は、任意の位置のＬＡＴの候補セットを、この位置から獲得されたユニポーラ電圧信号及びバイポーラ電圧信号に基づいて返すように構成された関数を利用することに留意されたい。そのような関数は、例えば、その開示が参照により本明細書に組み込まれる、Ｈｏｕｂｅｎらの米国特許第９，３８０，９５３号に記載されている。

アルゴリズム１６０の各反復は、プロセッサがプローブに属する電極Ｅ１を選択する、電極選択ステップ１４４から始まる。Ｅ１の選択に続いて、プロセッサは、近接選択ステップ１６２において、Ｅ１に近接する（すなわち、Ｅ１に隣接する）電極Ｅ２を選択する。次いで、プロセッサは、信号入力ステップ１６４において、前述の関数に２つの信号、すなわち、Ｅ１と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号と、Ｅ１とＥ２との間のバイポーラ電圧を表すバイポーラ電圧信号と、を入力する。続いて、出力受信ステップ１６６で、プロセッサは、入力に基づいて関数によって計算されたＬＡＴの候補セットを、関数からの出力として受信する。例えば、この出力は、候補ＬＡＴをマーキングする注釈を有するバイポーラ電圧信号を含み得る。

続いて、プロセッサは、チェックステップ１６８において、Ｅ１が近近接電極を他にも有するかどうかをチェックする。「はい」の場合、プロセッサは、近接選択ステップ１６２に戻り、次のＥ１の近接電極を選択する。そうでない場合、プロセッサは、ＬＡＴ選択ステップ１７０において、受信された全ての候補セットからＬＡＴを選定する。例えば、プロセッサは、他の候補ＬＡＴと比較してユニポーラ信号の導関数が最大である候補ＬＡＴを選択することができる。

したがって、例えば、各電極が最大４つの等距離近接電極を有する電極の矩形格子を与えられると、プロセッサは、最大４つの候補セットからのＬＡＴを選定することができる。各電極が最大６つの等距離近接電極を有する六角配列を与えられると、プロセッサは、最大６つの候補セットからＬＡＴを選定することができる。

以下の実施例は、本明細書の教示を組み合わせるか又は適用することができる、様々な非網羅的な方法に関する。以下の実施例は、本出願における又は本出願の後の書類提出における任意の時点で提示され得るいずれの特許請求の範囲の適用範囲も限定することを意図したものではないことを理解されたい。一切の権利放棄を意図するものではない。以下の実施例は、あくまでも例示的な目的で与えられるものに過ぎない。本明細書の様々な教示は、その他の多くの方法で構成及び適用さてもよいことが企図される。また、いくつかの変形例では、以下の実施例において言及されるある特定の特徴を省略してよいことも企図される。したがって、発明者らによって又は発明者らの権利承継人によって、後日そうである旨が明示的に示されない限り、以下に言及される態様又は特徴のいずれも重要なものとしてみなされるべきではない。いずれかの特許請求の範囲が、本出願において、又は以下に言及される特徴以外の更なる特徴を含む本出願に関連する後の書類提出において示される場合、それらの更なる特徴は、特許性に関連するいかなる理由によっても追加されたものとして仮定されるべきではない。

（実施例１）
システムは、ディスプレイと、心臓の解剖学的表面上の異なるそれぞれの測定位置において、複数の局所興奮時間（ＬＡＴ）を取得するように構成されたプロセッサと、を含む。プロセッサは、解剖学的表面上の１つ以上のサンプリング位置における電気伝播のそれぞれの方向を、サンプリング位置のうちの各サンプリング位置について、サンプリング位置に対する測定位置のそれぞれのサブセットを選択することと、ベクトルのセットを構築することであって、ベクトルのうちの少なくともいくつかの各々は、サブセット内の異なるそれぞれの測定位置について、測定位置のそれぞれの位置座標から導出された３つの位置値、及び測定位置においてＬＡＴから導出されたＬＡＴ値を含む、構築することと、ベクトルのセットの４×４の共分散行列の主成分分析（ＰＣＡ）に基づいて、サンプリング位置における電気伝播の方向を計算することと、によって計算するように更に構成されている。プロセッサは、ディスプレイ上に、電気伝播の方向を示すように更に構成されている。

（実施例２）
プロセッサは、
測定位置のサブセット全体にわたるＬＡＴの分散に基づいて、スケーリング係数を計算することと、
サブセット内の各測定位置について、
スケーリング係数によって、測定位置の位置座標及び測定位置におけるＬＡＴからなるパラメータの群から選択されるパラメータをスケーリングすることと、
スケーリングされたパラメータから、測定位置に対応するベクトルを構築することと、
によって、ベクトルのセットを構築するように構成されている、実施例１に記載のシステム。

（実施例３）
プロセッサは、共分散行列の第１の主成分を位置座標のそれぞれの次元へと射影することによって、サンプリング位置における電気伝播の方向を計算するように構成されている、実施例１又は２に記載のシステム。

（実施例４）
プロセッサは、サンプリング位置の各々における電気伝播の速度を、
サンプリング位置を通過し、サンプリング位置における電気伝播の方向に配向された仮想線について、仮想線への測定位置のサブセットのそれぞれの射影が存在する仮想線に沿ったそれぞれの距離を計算することと、
回帰点のグループに適合された回帰関数の勾配として速度を計算することであって、回帰点の各々は、サブセットに属する異なるそれぞれの測定位置について、（ｉ）測定位置の射影が存在する線に沿った距離、及び（ｉｉ）測定位置におけるＬＡＴを含む、計算することと、
によって計算するように更に構成されている、実施例１～３のいずれか１つに記載のシステム。

（実施例５）
プロセッサは、電気伝播の方向を示す前に電気伝播の方向を平滑化するように更に構成されている、実施例１～４のいずれか１つに記載のシステム。

（実施例６）
プロセッサは、サンプリング位置に対する測定位置のそれぞれのサブセットを、
測定位置のうちの、サンプリング位置から所定の距離内にあるものを特定することと、
特定された測定位置から測定位置のサブセットを選択することと、
によって選択するように構成されている、実施例１～５のいずれか１つに記載のシステム。

（実施例７）
プロセッサは、サブセットを、
サンプリング位置から所定の距離内にある各測定位置について、
サンプリング位置と測定位置との間の伝播速度を推定する伝播速度推定値を計算することと、
この伝播速度が所定の速度推定閾値を超えるという条件で、測定位置を選択することと、
によって選択するように構成されている、実施例６に記載のシステム。

（実施例８）
プロセッサは、サンプリング位置の補間されたＬＡＴを計算するように更に構成されており、プロセッサは、この補間されたＬＡＴに基づいて伝播速度推定値を計算するように構成されている、実施例７に記載のシステム。

（実施例９）
所定の速度推定閾値は第１の所定の速度推定閾値であり、伝播速度推定値は、第１の伝播速度推定値であり、
プロセッサは、補間されたＬＡＴを、
クラスタの各々について、クラスタ内の測定位置の各ペアについて、測定位置のペア間の伝播速度を推定する第２の伝播速度推定値が、第２の所定の速度推定閾値を超えるように、測定位置のうちの、サンプリング位置から所定の距離内にあるものを、１つ以上のクラスタへとクラスタ化することと、
サンプリング位置とクラスタとの間のそれぞれの距離に基づいて、クラスタのうちの１つを特定することと、
クラスタのうちの特定された１つにおける測定位置のうちの少なくともいくつかにおいて、補間されたＬＡＴをＬＡＴの加重平均として計算することと、
によって計算するように構成されている、実施例８に記載のシステム。

（実施例１０）
測定位置は、解剖学的表面を表すデジタルモデル表面上の異なるそれぞれの測定点に対応し、これらの測定点は、ＬＡＴにそれぞれ関連付けられ、
プロセッサは、デジタルモデル表面上で、複数のサンプリング点を指定するように更に構成されており、
サンプリング位置は、サンプリング点にそれぞれ対応する、
実施例１～９のいずれか１つに記載のシステム。

（実施例１１）
プロセッサは、それぞれのマーカーがサンプリング点においてモデル表面の上に重なり、電気伝播の方向に配向された状態で、モデル表面を表示することによって、電気伝播の方向をそれぞれ示すように構成されている、実施例１０に記載のシステム。

（実施例１２）
プロセッサは、
サンプリング位置における電気伝播のそれぞれの速度を計算し、
この速度に従ってマーカーの少なくとも１つの特性を変化させる、
ように更に構成されている、実施例１１に記載のシステム。

（実施例１３）
プロセッサは、マーカーがモデル表面の上に重なった状態でモデル表面を表示する前に、
反復して、
電気伝播の方向を再計算し、
電気伝播の方向に応じて、サンプリング点を互いに向かってシフトさせる、
ように更に構成されている、実施例１１又は１２に記載のシステム。

（実施例１４）
プロセッサは、体内プローブに属するそれぞれの電極によって獲得された信号に基づいて、ＬＡＴを計算することによって、ＬＡＴを取得するように構成されており、
サンプリング位置は、測定位置のうちの少なくともいくつかを含み、
プロセッサは、電極が測定位置にそれぞれある間に、電気伝播の方向を示すように構成されている、
実施例１～９のいずれか１つに記載のシステム。

（実施例１５）
プロセッサは、電気伝播の方向を、
プローブのアイコンを表示することと、
電気伝播の方向に配向されたそれぞれのマーカーを、電極のうちの、サンプリング位置に位置するものに対応するアイコンの部分にそれぞれ配置することと、
によって示すように構成されている、実施例１４に記載のシステム。

（実施例１６）
プロセッサは、
サンプリング位置における電気伝播のそれぞれの速度を計算し、
この速度に従ってマーカーの少なくとも１つの特性を変化させる、
ように更に構成されている、実施例１５に記載のシステム。

（実施例１７）
方法は、心臓の解剖学的表面上の異なるそれぞれの測定位置において、複数の局所興奮時間（ＬＡＴ）を取得することを含む。この方法は、解剖学的表面上の１つ以上のサンプリング位置における電気伝播のそれぞれの方向を、サンプリング位置のうちの各サンプリング位置について、サンプリング位置に対する測定位置のそれぞれのサブセットを選択することと、ベクトルのセットを構築することであって、ベクトルのうちの少なくともいくつかの各々は、サブセット内の異なるそれぞれの測定位置について、測定位置のそれぞれの位置座標から導出された３つの位置値、及び測定位置においてＬＡＴから導出されたＬＡＴ値を含む、構築することと、ベクトルのセットの４×４の共分散行列の主成分分析（ＰＣＡ）に基づいて、サンプリング位置における電気伝播の方向を計算することと、によって計算することを更に含む。この方法は、電気伝播の方向をディスプレイ上に示すことを更に含む。

（実施例１８）
ベクトルのセットを構築することは、
測定位置のサブセット全体にわたるＬＡＴの分散に基づいて、スケーリング係数を計算することと、
サブセット内の各測定位置について、
スケーリング係数によって、測定位置の位置座標及び測定位置におけるＬＡＴからなるパラメータの群から選択されるパラメータをスケーリングすることと、
スケーリングされたパラメータから、測定位置に対応するベクトルを構築することと、
を含む、実施例１７に記載の方法。

（実施例１９）
サンプリング位置における電気伝播の方向を計算することは、共分散行列の第１の主成分を位置座標のそれぞれの次元へと射影することによって、電気伝播の方向を計算することを含む、実施例１７又は１８に記載の方法。

（実施例２０）
サンプリング位置の各々における電気伝播の速度を、
サンプリング位置を通過し、サンプリング位置における電気伝播の方向に配向された仮想線について、仮想線への測定位置のサブセットのそれぞれの射影が存在する仮想線に沿ったそれぞれの距離を計算することと、
回帰点のグループに適合された回帰関数の勾配として速度を計算することであって、回帰点の各々は、サブセットに属する異なるそれぞれの測定位置について、（ｉ）測定位置の射影が存在する線に沿った距離、及び（ｉｉ）測定位置におけるＬＡＴを含む、計算することと、
によって計算するように更に構成されている、実施例１７～１９のいずれか１つに記載の方法。

（実施例２１）
測定位置は、解剖学的表面を表すデジタルモデル表面上の異なるそれぞれの測定点に対応し、これらの測定点は、ＬＡＴにそれぞれ関連付けられ、
この方法は、デジタルモデル表面上で、複数のサンプリング点を指定することを更に含み、
サンプリング位置は、サンプリング点にそれぞれ対応する、
実施例１７～２０のいずれか１つに記載の方法。

（実施例２２）
システム（２０）は、電気的インターフェース（３４）と、プロセッサ（３２）と、を含む。プロセッサ（３２）は、電気的インターフェース（３４）を介して、心臓（２４）の解剖学的表面（４２）上の複数の電極（２８）によって獲得されたそれぞれの信号を受信するように構成されている。プロセッサ（３２）は、これらの信号に基づいて、電極（２８）のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算するように更に構成されている。プロセッサ（３２）は、ＬＡＴに基づいて、これらの位置における電気伝播のそれぞれの方向を計算するように更に構成されている。プロセッサ（３２）は、電極（２８）のうちの隣接するもののペア（１３９）を、ペア（１３９）の各々について、ペア（１３９）を結合するベクトル（１４０）が、所定の整列閾値度内で、ペア（１３９）に属する電極（２８）のうちの１つの位置における電気伝播の方向と整列されるように選択するように更に構成されている。プロセッサ（３２）は、電極（２８）のペア（１３９）によって測定されたそれぞれのバイポーラ電圧を解剖学的表面（４２）のデジタルモデル（３８）と関連付けるように更に構成されている。

（実施例２３）
プロセッサ（３２）は、電極（２８）のうちの各第１の電極の位置におけるＬＡＴを、
その位置に対するＬＡＴの複数の候補セットを取得することであって、第１の電極に隣接する、電極（２８）のうちの各第２の電極について、
関数への入力として、（ｉ）第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）第１の電極と第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
関数からの出力として、候補セットのそれぞれ１つを受信することと、
によって取得することと、
候補セットから、ＬＡＴを選定することと、
によって、計算するように構成されている、実施例２２に記載のシステム（２０）。

（実施例２４）
方法は、心臓（２４）の解剖学的表面（４２）上の複数の電極（２８）によって獲得されたそれぞれの信号に基づいて、電極（２８）のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算することを含む。この方法は、ＬＡＴに基づいて、これらの位置における電気伝播のそれぞれの方向を計算することを更に含む。この方法は、電極（２８）のうちの隣接するもののペア（１３９）を、ペア（１３９）の各々について、ペアを結合するベクトル（１４０）が、所定の整列閾値度内で、ペア（１３９）に属する電極（２８）のうちの１つの位置における電気伝播の方向と整列されるように選択することを更に含む。この方法は、電極（２８）のペア（１３９）によって測定されたそれぞれのバイポーラ電圧を解剖学的表面（４２）のデジタルモデル（３８）と関連付けることを更に含む。

（実施例２５）
ＬＡＴを計算することは、電極（２８）のうちの各第１の電極の位置におけるＬＡＴを、
その位置に対するＬＡＴの複数の候補セットを取得することであって、第１の電極に隣接する、電極（２８）のうちの各第２の電極について、
関数への入力として、（ｉ）第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）第１の電極と第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
関数からの出力として、候補セットのそれぞれ１つを受信することと、
によって取得することと、
候補セットから、ＬＡＴを選定することと、
によって計算することを含む、実施例２４に記載の方法。

（実施例２６）
コンピュータソフトウェア製品は、プログラム命令が格納されている有形の非一時的なコンピュータ可読媒体を含む。命令は、プロセッサ（３２）によって読み取られると、プロセッサ（３２）に、心臓（２４）の解剖学的表面（４２）上の複数の電極（２８）によって獲得されたそれぞれの信号を受信させ、これらの信号に基づいて、電極（２８）のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算させ、それらのＬＡＴに基づいて、それぞれの位置における電気伝播のそれぞれの方向を計算させ、電極（２８）のうちの隣接するもののペア（１３９）を、ペア（１３９）の各々について、ペア（１３９）を結合するベクトル（１４０）が、所定の整列閾値度内で、ペア（１３９）に属する電極（２８）のうちの１つの位置における電気伝播の方向と整列されるように選択させ、電極（２８）のペア（１３９）によって測定されたそれぞれのバイポーラ電圧を解剖学的表面（４２）のデジタルモデル（３８）と関連付けさせる。

（実施例２７）
システム（２０）は、ディスプレイ（３６）及びプロセッサ（３２）を含む。プロセッサ（３２）は、心臓（２４）の解剖学的表面（４２）上の、体内プローブ（２６）に属する複数の電極（２８）によって獲得されたそれぞれの信号に基づいて、電極（２８）のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算するように構成されている。プロセッサ（３２）は、ＬＡＴに基づいて、これら位置における電気伝播のそれぞれの方向及び速度を計算するように更に構成されている。プロセッサは、電極（２８）がそれぞれ当該位置にある間に、ディスプレイ（３６）上にプローブ（２６）のアイコン（２６’）を表示するように更に構成されている。プロセッサ（３２）は、電気伝播の方向にそれぞれ配向され、速度に応じて変化する少なくとも１つの特性を有するそれぞれのマーカー（１１７）を、電極（２８）に対応するアイコン（２６’）の部分に配置するように更に構成されている。

（実施例２８）
プロセッサ（３２）は、電極（２８）のうちの各第１の電極の位置におけるＬＡＴを、
その位置に対するＬＡＴの複数の候補セットを取得することであって、第１の電極に隣接する、電極のうちの各第２の電極について、
関数への入力として、（ｉ）第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）第１の電極と第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
関数からの出力として、候補セットのそれぞれ１つを受信することと、
によって取得することと、
候補セットから、ＬＡＴを選定することと、
によって計算するように構成されている、実施例２７に記載のシステム（２０）。

（実施例２９）
特性は、マーカー（１１７）が、
第１の範囲に属する速度のものについて、第１の形状及び第１の厚さを有し、
第１の範囲よりも低い第２の範囲に属する速度のものについて、第１の形状及び第２の厚さを有し、
第２の範囲よりも低い第３の範囲に属する速度のものについて、第２の形状を有する、
ことを理由に、速度に応じて変化する、実施例２７又は２８に記載のシステム（２０）。

（実施例３０）
マーカー（１１７）は、ＬＡＴに応じて変化する少なくとも１つの他の特性を有する、実施例２９に記載のシステム（２０）。

（実施例３１）
マーカー（１１７）は、ＬＡＴに基づいて、カラースケールに従って着色される、実施例３０に記載のシステム（２０）。

（実施例３２）
方法は、心臓（２４）の解剖学的表面（４２）上の、体内プローブ（２６）に属する複数の電極（２８）によって獲得されたそれぞれの信号に基づいて、電極（２８）のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算することを含む。この方法は、ＬＡＴに基づいて、これらの位置における電気伝播のそれぞれの方向及び速度を計算することを更に含む。この方法は、電極（２８）がそれぞれ当該位置にある間に、プローブ（２６）のアイコン（２６’）を表示することを更に含む。この方法は、電気伝播の方向にそれぞれ配向され、速度に応じて変化する少なくとも１つの特性を有するそれぞれのマーカー（１１７）を、電極（２８）に対応するアイコン（２６’）の部分に配置することを更に含む。

（実施例３３）
ＬＡＴを計算することは、電極（２８）のうちの各第１の電極の位置におけるＬＡＴを、
その位置に対するＬＡＴの複数の候補セットを取得することであって、第１の電極に隣接する、電極のうちの各第２の電極について、
関数への入力として、（ｉ）第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）第１の電極と第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
関数からの出力として、候補セットのそれぞれ１つを受信することと、
によって取得することと、
候補セットから、ＬＡＴを選定することと、
によって計算することを含む、実施例３２に記載の方法。

（実施例３４）
特性は、マーカー（１１７）が、
第１の範囲に属する速度のものについて、第１の形状及び第１の厚さを有し、
第１の範囲よりも低い第２の範囲に属する速度のものについて、第１の形状及び第２の厚さを有し、
第２の範囲よりも低い第３の範囲に属する速度のものについて、第２の形状を有する、
ことを理由に、速度に応じて変化する、実施例３２又は３３に記載の方法。

（実施例３５）
マーカー（１１７）は、ＬＡＴに応じて変化する少なくとも１つの他の特性を有する、実施例３４に記載の方法。

（実施例３６）
マーカー（１１７）は、ＬＡＴに基づいて、カラースケールに従って着色される、実施例３５に記載の方法。

（実施例３７）
コンピュータソフトウェア製品は、プログラム命令が格納されている有形の非一時的なコンピュータ可読媒体を含む。命令は、プロセッサ（３２）によって読み取られると、プロセッサ（３２）に、心臓（２４）の解剖学的表面（４２）上の、体内プローブ（２６）に属する複数の電極（２８）によって獲得されたそれぞれの信号に基づいて、電極（２８）のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算させる。命令は、更にプロセッサ（３２）に、ＬＡＴに基づいて、これらの位置における電気伝播のそれぞれの方向及び速度を計算させる。命令は、更にプロセッサ（３２）に、電極（２８）がそれぞれ当該位置にある間に上にプローブ（２６）のアイコン（２６’）を表示させる。命令は、更にプロセッサ（３２）に、電気伝播の方向にそれぞれ配向され、速度に応じて変化する少なくとも１つの特性を有するそれぞれのマーカー（１１７）を、電極（２８）に対応するアイコン（２６’）の部分に配置させる。

本開示は、本明細書に具体的に示され、上述されたものに限定されない点が、当業者には理解されよう。むしろ、本開示の範囲は、本明細書において上述した様々な特徴の組み合わせ及び部分的組み合わせ、並びに上述の説明を読むことで当業者が想到するであろう、従来技術にはない特徴の変形例及び修正例を含む。参照により本特許出願に組み込まれる文献は、これらの組み込まれた文献において、いずれかの用語が本明細書において明示的又は暗示的になされた定義と矛盾する様式で定義される程度まで、本明細書における定義のみを考慮するものとする点を除き、本出願の不可欠な部分と見なすものとする。

〔実施の態様〕
（１） システムであって、
電気的インターフェースと、
プロセッサであって、
前記電気的インターフェースを介して、心臓の解剖学的表面上の複数の電極によって獲得されたそれぞれの信号を受信し、
前記信号に基づいて、前記電極のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算し、
前記ＬＡＴに基づいて、前記位置における電気伝播のそれぞれの方向を計算し、
前記電極のうちの隣接する電極のペアを、前記ペアの各々について、前記ペアを結合するベクトルが、所定の閾値整列度内で、前記ペアに属する前記電極のうちの１つの前記位置における電気伝播の方向と整列されるように選択し、
前記電極のペアによって測定されたそれぞれのバイポーラ電圧を、前記解剖学的表面のデジタルモデルと関連付ける、
ように構成されている、プロセッサと、
を含む、システム。
（２） 前記プロセッサは、前記電極のうちの各第１の電極の前記位置における前記ＬＡＴを、
前記位置に対するＬＡＴの複数の候補セットを取得することであって、前記第１の電極に隣接する、前記電極のうちの各第２の電極について、
関数への入力として、（ｉ）前記第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）前記第１の電極と前記第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
前記関数からの出力として、前記候補セットのそれぞれ１つを受信することと、
によって取得することと、
前記候補セットから、前記ＬＡＴを選定することと、
によって計算するように構成されている、実施態様１に記載のシステム。
（３） 心臓の解剖学的表面上の複数の電極によって獲得されたそれぞれの信号に基づいて、前記電極のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算することと、
前記ＬＡＴに基づいて、前記位置における電気伝播のそれぞれの方向を計算することと、
前記電極のうちの隣接する電極のペアを、前記ペアの各々について、前記ペアを結合するベクトルが、所定の閾値整列度内で、前記ペアに属する前記電極のうちの１つの前記位置における電気伝播の方向と整列されるように選択することと、
前記電極のペアによって測定されたそれぞれのバイポーラ電圧を、前記解剖学的表面のデジタルモデルと関連付けることと、
を含む方法。
（４） 前記ＬＡＴを計算することは、前記電極のうちの各第１の電極の前記位置における前記ＬＡＴを、
前記位置に対するＬＡＴの複数の候補セットを取得することであって、前記第１の電極に隣接する、前記電極のうちの各第２の電極について、
関数への入力として、（ｉ）前記第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）前記第１の電極と前記第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
前記関数からの出力として、前記候補セットのそれぞれ１つを受信することと、
によって取得することと、
前記候補セットから、前記ＬＡＴを選定することと、
によって計算することを含む、実施態様３に記載の方法。
（５） プログラム命令が格納されている有形の非一時的なコンピュータ可読媒体を備えるコンピュータソフトウェア製品であって、前記命令は、プロセッサによって読み取られると、前記プロセッサに、
心臓の解剖学的表面上の複数の電極によって獲得されたそれぞれの信号を受信させ、
前記信号に基づいて、前記電極のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算させ、
前記ＬＡＴに基づいて、前記位置における電気伝播のそれぞれの方向を計算させ、
前記電極のうちの隣接する電極のペアを、前記ペアの各々について、前記ペアを結合するベクトルが、所定の閾値整列度内で、前記ペアに属する前記電極のうちの１つの前記位置における電気伝播の方向と整列されるように選択させ、
前記電極のペアによって測定されたそれぞれのバイポーラ電圧を前記解剖学的表面のデジタルモデルと関連付けさせる、
コンピュータソフトウェア製品。

（６） 前記命令は、前記プロセッサに、前記電極のうちの各第１の電極の前記位置における前記ＬＡＴを、
前記位置に対するＬＡＴの複数の候補セットを取得することであって、前記第１の電極に隣接する、前記電極のうちの各第２の電極について、
関数への入力として、（ｉ）前記第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）前記第１の電極と前記第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
前記関数からの出力として、前記候補セットのそれぞれ１つを受信することと、
によって取得することと、
前記候補セットから前記ＬＡＴを選定することと、
によって計算させる、実施態様５に記載のコンピュータソフトウェア製品。
（７） システムであって、
ディスプレイと、
プロセッサであって、
心臓の解剖学的表面上の、体内プローブに属する複数の電極によって獲得されたそれぞれの信号に基づいて、前記電極のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算し、
前記ＬＡＴに基づいて、前記位置における電気伝播のそれぞれの方向及び速度を計算し、
前記電極がそれぞれ前記位置にある間に、前記ディスプレイ上に前記プローブのアイコンを表示し、
前記電気伝播の方向にそれぞれ配向され、前記速度に応じて変化する少なくとも１つの特性を有するそれぞれのマーカーを、前記電極に対応する前記アイコンの部分に配置する、
ように構成されている、プロセッサと、
を含むシステム。
（８） 前記プロセッサは、前記電極のうちの各第１の電極の前記位置における前記ＬＡＴを、
前記位置に対するＬＡＴの複数の候補セットを取得することであって、前記第１の電極に隣接する、前記電極のうちの各第２の電極について、
関数への入力として、（ｉ）前記第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）前記第１の電極と前記第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
前記関数からの出力として、前記候補セットのそれぞれ１つを受信することと、
によって取得することと、
前記候補セットから前記ＬＡＴを選定することと、
によって計算するように構成されている、実施態様７に記載のシステム。
（９） 前記特性は、前記マーカーが、
第１の範囲に属する速度のものについて、第１の形状及び第１の厚さを有し、
前記第１の範囲よりも低い第２の範囲に属する速度のものについて、前記第１の形状及び第２の厚さを有し、
前記第２の範囲よりも低い第３の範囲に属する速度のものについて、第２の形状を有する、
ことを理由に、前記速度に応じて変化する、実施態様７に記載のシステム。
（１０） 前記マーカーは、前記ＬＡＴに応じて変化する少なくとも１つの他の特性を有する、実施態様９に記載のシステム。

（１１） 前記マーカーは、前記ＬＡＴに基づいて、カラースケールに従って着色される、実施態様１０に記載のシステム。
（１２） 心臓の解剖学的表面上の、体内プローブに属する複数の電極によって獲得されたそれぞれの信号に基づいて、前記電極のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算することと、
前記ＬＡＴに基づいて、前記位置における電気伝播のそれぞれの方向及び速度を計算することと、
前記電極がそれぞれ前記位置にある間に、前記プローブのアイコンを表示することと、
前記電気伝播の方向にそれぞれ配向され、前記速度に応じて変化する少なくとも１つの特性を有するそれぞれのマーカーを、前記電極に対応する前記アイコンの部分に配置することと、
を含む方法。
（１３） 前記ＬＡＴを計算することは、前記電極のうちの各第１の電極の前記位置における前記ＬＡＴを、
前記位置に対するＬＡＴの複数の候補セットを取得することであって、前記第１の電極に隣接する、前記電極のうちの各第２の電極について、
関数への入力として、（ｉ）前記第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）前記第１の電極と前記第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
前記関数からの出力として、前記候補セットのそれぞれ１つを受信することと、
によって取得することと、
前記候補セットから前記ＬＡＴを選定することと、
によって計算することを含む、実施態様１２に記載の方法。
（１４） 前記特性は、前記マーカーが、
第１の範囲に属する速度のものについて、第１の形状及び第１の厚さを有し、
前記第１の範囲よりも低い第２の範囲に属する速度のものについて、前記第１の形状及び第２の厚さを有し、
前記第２の範囲よりも低い第３の範囲に属する速度のものについて、第２の形状を有する、
ことを理由に、前記速度に応じて変化する、実施態様１２に記載の方法。
（１５） 前記マーカーは、前記ＬＡＴに応じて変化する少なくとも１つの他の特性を有する、実施態様１４に記載の方法。

（１６） 前記マーカーは、前記ＬＡＴに基づいて、カラースケールに従って着色される、実施態様１５に記載の方法。
（１７） プログラム命令が格納されている有形の非一時的なコンピュータ可読媒体を備えるコンピュータソフトウェア製品であって、前記命令は、プロセッサによって読み取られると、前記プロセッサに、
心臓の解剖学的表面上の、体内プローブに属する複数の電極によって獲得されたそれぞれの信号に基づいて、前記電極のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算させ、
前記ＬＡＴに基づいて、前記位置における電気伝播のそれぞれの方向及び速度を計算させ、
前記電極がそれぞれ前記位置にある間に、前記プローブのアイコンを表示させ、
前記電気伝播の方向にそれぞれ配向され、前記速度に応じて変化する少なくとも１つの特性を有するそれぞれのマーカーを、前記電極に対応する前記アイコンの部分に配置させる、
コンピュータソフトウェア製品。
（１８） 前記命令は、前記プロセッサに、前記電極のうちの各第１の電極の前記位置における前記ＬＡＴを、
前記位置に対するＬＡＴの複数の候補セットを取得することであって、前記第１の電極に隣接する、前記電極のうちの各第２の電極について、
関数への入力として、（ｉ）前記第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）前記第１の電極と前記第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
前記関数からの出力として、前記候補セットのそれぞれ１つを受信することと、
によって取得することと、
前記候補セットから前記ＬＡＴを選定することと、
によって計算させる、実施態様１７に記載のコンピュータソフトウェア製品。
（１９） 前記特性は、前記マーカーが、
第１の範囲に属する速度のものについて、第１の形状及び第１の厚さを有し、
前記第１の範囲よりも低い第２の範囲に属する速度のものについて、前記第１の形状及び第２の厚さを有し、
前記第２の範囲よりも低い第３の範囲に属する速度のものについて、第２の形状を有する、
ことを理由に、前記速度に応じて変化する、実施態様１７に記載のコンピュータソフトウェア製品。
（２０） 前記マーカーは、前記ＬＡＴに応じて変化する少なくとも１つの他の特性を有する、実施態様１９に記載のコンピュータソフトウェア製品。

Claims (20)

1. システムであって、
   電気的インターフェースと、
   プロセッサであって、
   前記電気的インターフェースを介して、心臓の解剖学的表面上の複数の電極によって獲得されたそれぞれの信号を受信し、
   前記信号に基づいて、前記電極のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算し、
   前記ＬＡＴに基づいて、前記位置における電気伝播のそれぞれの方向を計算し、
   前記電極のうちの隣接する電極のペアを、前記ペアの各々について、前記ペアを結合するベクトルが、所定の閾値整列度内で、前記ペアに属する前記電極のうちの１つの前記位置における電気伝播の方向と整列されるように選択し、
   前記電極のペアによって測定されたそれぞれのバイポーラ電圧を、前記解剖学的表面のデジタルモデルと関連付ける、
   ように構成されている、プロセッサと、
   を含む、システム。
2. 前記プロセッサは、前記電極のうちの各第１の電極の前記位置における前記ＬＡＴを、
   前記位置に対するＬＡＴの複数の候補セットを取得することであって、前記第１の電極に隣接する、前記電極のうちの各第２の電極について、
   関数への入力として、（ｉ）前記第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）前記第１の電極と前記第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
   前記関数からの出力として、前記候補セットのそれぞれ１つを受信することと、
   によって取得することと、
   前記候補セットから、前記ＬＡＴを選定することと、
   によって計算するように構成されている、請求項１に記載のシステム。
3. 心臓の解剖学的表面上の複数の電極によって獲得されたそれぞれの信号に基づいて、前記電極のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算することと、
   前記ＬＡＴに基づいて、前記位置における電気伝播のそれぞれの方向を計算することと、
   前記電極のうちの隣接する電極のペアを、前記ペアの各々について、前記ペアを結合するベクトルが、所定の閾値整列度内で、前記ペアに属する前記電極のうちの１つの前記位置における電気伝播の方向と整列されるように選択することと、
   前記電極のペアによって測定されたそれぞれのバイポーラ電圧を、前記解剖学的表面のデジタルモデルと関連付けることと、
   を含む方法。
4. 前記ＬＡＴを計算することは、前記電極のうちの各第１の電極の前記位置における前記ＬＡＴを、
   前記位置に対するＬＡＴの複数の候補セットを取得することであって、前記第１の電極に隣接する、前記電極のうちの各第２の電極について、
   関数への入力として、（ｉ）前記第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）前記第１の電極と前記第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
   前記関数からの出力として、前記候補セットのそれぞれ１つを受信することと、
   によって取得することと、
   前記候補セットから、前記ＬＡＴを選定することと、
   によって計算することを含む、請求項３に記載の方法。
5. プログラム命令が格納されている有形の非一時的なコンピュータ可読媒体を備えるコンピュータソフトウェア製品であって、前記命令は、プロセッサによって読み取られると、前記プロセッサに、
   心臓の解剖学的表面上の複数の電極によって獲得されたそれぞれの信号を受信させ、
   前記信号に基づいて、前記電極のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算させ、
   前記ＬＡＴに基づいて、前記位置における電気伝播のそれぞれの方向を計算させ、
   前記電極のうちの隣接する電極のペアを、前記ペアの各々について、前記ペアを結合するベクトルが、所定の閾値整列度内で、前記ペアに属する前記電極のうちの１つの前記位置における電気伝播の方向と整列されるように選択させ、
   前記電極のペアによって測定されたそれぞれのバイポーラ電圧を前記解剖学的表面のデジタルモデルと関連付けさせる、
   コンピュータソフトウェア製品。
6. 前記命令は、前記プロセッサに、前記電極のうちの各第１の電極の前記位置における前記ＬＡＴを、
   前記位置に対するＬＡＴの複数の候補セットを取得することであって、前記第１の電極に隣接する、前記電極のうちの各第２の電極について、
   関数への入力として、（ｉ）前記第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）前記第１の電極と前記第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
   前記関数からの出力として、前記候補セットのそれぞれ１つを受信することと、
   によって取得することと、
   前記候補セットから前記ＬＡＴを選定することと、
   によって計算させる、請求項５に記載のコンピュータソフトウェア製品。
7. システムであって、
   ディスプレイと、
   プロセッサであって、
   心臓の解剖学的表面上の、体内プローブに属する複数の電極によって獲得されたそれぞれの信号に基づいて、前記電極のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算し、
   前記ＬＡＴに基づいて、前記位置における電気伝播のそれぞれの方向及び速度を計算し、
   前記電極がそれぞれ前記位置にある間に、前記ディスプレイ上に前記プローブのアイコンを表示し、
   前記電気伝播の方向にそれぞれ配向され、前記速度に応じて変化する少なくとも１つの特性を有するそれぞれのマーカーを、前記電極に対応する前記アイコンの部分に配置する、
   ように構成されている、プロセッサと、
   を含むシステム。
8. 前記プロセッサは、前記電極のうちの各第１の電極の前記位置における前記ＬＡＴを、
   前記位置に対するＬＡＴの複数の候補セットを取得することであって、前記第１の電極に隣接する、前記電極のうちの各第２の電極について、
   関数への入力として、（ｉ）前記第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）前記第１の電極と前記第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
   前記関数からの出力として、前記候補セットのそれぞれ１つを受信することと、
   によって取得することと、
   前記候補セットから前記ＬＡＴを選定することと、
   によって計算するように構成されている、請求項７に記載のシステム。
9. 前記特性は、前記マーカーが、
   第１の範囲に属する速度のものについて、第１の形状及び第１の厚さを有し、
   前記第１の範囲よりも低い第２の範囲に属する速度のものについて、前記第１の形状及び第２の厚さを有し、
   前記第２の範囲よりも低い第３の範囲に属する速度のものについて、第２の形状を有する、
   ことを理由に、前記速度に応じて変化する、請求項７に記載のシステム。
10. 前記マーカーは、前記ＬＡＴに応じて変化する少なくとも１つの他の特性を有する、請求項９に記載のシステム。
11. 前記マーカーは、前記ＬＡＴに基づいて、カラースケールに従って着色される、請求項１０に記載のシステム。
12. 心臓の解剖学的表面上の、体内プローブに属する複数の電極によって獲得されたそれぞれの信号に基づいて、前記電極のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算することと、
    前記ＬＡＴに基づいて、前記位置における電気伝播のそれぞれの方向及び速度を計算することと、
    前記電極がそれぞれ前記位置にある間に、前記プローブのアイコンを表示することと、
    前記電気伝播の方向にそれぞれ配向され、前記速度に応じて変化する少なくとも１つの特性を有するそれぞれのマーカーを、前記電極に対応する前記アイコンの部分に配置することと、
    を含む方法。
13. 前記ＬＡＴを計算することは、前記電極のうちの各第１の電極の前記位置における前記ＬＡＴを、
    前記位置に対するＬＡＴの複数の候補セットを取得することであって、前記第１の電極に隣接する、前記電極のうちの各第２の電極について、
    関数への入力として、（ｉ）前記第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）前記第１の電極と前記第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
    前記関数からの出力として、前記候補セットのそれぞれ１つを受信することと、
    によって取得することと、
    前記候補セットから前記ＬＡＴを選定することと、
    によって計算することを含む、請求項１２に記載の方法。
14. 前記特性は、前記マーカーが、
    第１の範囲に属する速度のものについて、第１の形状及び第１の厚さを有し、
    前記第１の範囲よりも低い第２の範囲に属する速度のものについて、前記第１の形状及び第２の厚さを有し、
    前記第２の範囲よりも低い第３の範囲に属する速度のものについて、第２の形状を有する、
    ことを理由に、前記速度に応じて変化する、請求項１２に記載の方法。
15. 前記マーカーは、前記ＬＡＴに応じて変化する少なくとも１つの他の特性を有する、請求項１４に記載の方法。
16. 前記マーカーは、前記ＬＡＴに基づいて、カラースケールに従って着色される、請求項１５に記載の方法。
17. プログラム命令が格納されている有形の非一時的なコンピュータ可読媒体を備えるコンピュータソフトウェア製品であって、前記命令は、プロセッサによって読み取られると、前記プロセッサに、
    心臓の解剖学的表面上の、体内プローブに属する複数の電極によって獲得されたそれぞれの信号に基づいて、前記電極のそれぞれの位置におけるそれぞれの局所興奮時間（ＬＡＴ）を計算させ、
    前記ＬＡＴに基づいて、前記位置における電気伝播のそれぞれの方向及び速度を計算させ、
    前記電極がそれぞれ前記位置にある間に、前記プローブのアイコンを表示させ、
    前記電気伝播の方向にそれぞれ配向され、前記速度に応じて変化する少なくとも１つの特性を有するそれぞれのマーカーを、前記電極に対応する前記アイコンの部分に配置させる、
    コンピュータソフトウェア製品。
18. 前記命令は、前記プロセッサに、前記電極のうちの各第１の電極の前記位置における前記ＬＡＴを、
    前記位置に対するＬＡＴの複数の候補セットを取得することであって、前記第１の電極に隣接する、前記電極のうちの各第２の電極について、
    関数への入力として、（ｉ）前記第１の電極と基準電極との間のユニポーラ電圧を表すユニポーラ電圧信号、及び（ｉｉ）前記第１の電極と前記第２の電極との間のバイポーラ電圧を表すバイポーラ電圧信号を提供することと、
    前記関数からの出力として、前記候補セットのそれぞれ１つを受信することと、
    によって取得することと、
    前記候補セットから前記ＬＡＴを選定することと、
    によって計算させる、請求項１７に記載のコンピュータソフトウェア製品。
19. 前記特性は、前記マーカーが、
    第１の範囲に属する速度のものについて、第１の形状及び第１の厚さを有し、
    前記第１の範囲よりも低い第２の範囲に属する速度のものについて、前記第１の形状及び第２の厚さを有し、
    前記第２の範囲よりも低い第３の範囲に属する速度のものについて、第２の形状を有する、
    ことを理由に、前記速度に応じて変化する、請求項１７に記載のコンピュータソフトウェア製品。
20. 前記マーカーは、前記ＬＡＴに応じて変化する少なくとも１つの他の特性を有する、請求項１９に記載のコンピュータソフトウェア製品。

Images