知財求人 - 知財ポータルサイト「IP Force」
▶ バイオセンス・ウエブスター・（イスラエル）・リミテッドの特許一覧
特許7625400不一致を検出することによる、ＥＣＧ信号のマッピングの遡及的最適化
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2025-01-24
(45)【発行日】2025-02-03
(54)【発明の名称】不一致を検出することによる、ＥＣＧ信号のマッピングの遡及的最適化
(51)【国際特許分類】
A61B 5/367 20210101AFI20250127BHJP
A61B 5/352 20210101ALI20250127BHJP
【ＦＩ】
A61B5/367
A61B5/352
【請求項の数】
16
【外国語出願】
(21)【出願番号】P 2020184241
(22)【出願日】2020-11-04
(65)【公開番号】P2021074543
(43)【公開日】2021-05-20
【審査請求日】2023-09-20
(31)【優先権主張番号】16/674,921
(32)【優先日】2019-11-05
(33)【優先権主張国・地域又は機関】US
(73)【特許権者】
【識別番号】511099630
【氏名又は名称】バイオセンス・ウエブスター・（イスラエル）・リミテッド
【氏名又は名称原語表記】Ｂｉｏｓｅｎｓｅ Ｗｅｂｓｔｅｒ （Ｉｓｒａｅｌ）， Ｌｔｄ．
(74)【代理人】
【識別番号】100088605
【弁理士】
【氏名又は名称】加藤 公延
(74)【代理人】
【識別番号】100130384
【弁理士】
【氏名又は名称】大島 孝文
(72)【発明者】
【氏名】アサフ・ゴバリ
(72)【発明者】
【氏名】バディム・グリナー
【審査官】▲高▼原 悠佑
(56)【参考文献】
【文献】特開２０１５－１３９７０７（ＪＰ，Ａ）
【文献】米国特許出願公開第２０１９／０２６１９２７（ＵＳ，Ａ１）
【文献】特表２０１６－５３０００８（ＪＰ，Ａ）
【文献】米国特許出願公開第２０１８／０３６８７１３（ＵＳ，Ａ１）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ６１Ｂ ５／２４－５／３９８
(57)【特許請求の範囲】
【請求項1】
患者の心臓内の心臓内プローブの複数の電極によって取得された、複数の心内信号を受信するように構成されている信号取得回路と、プロセッサであって、
後続の時点で、一連のアノテーション可視化操作を、各操作において、
（ｉ）前記心内信号から複数のアノテーション値を抽出することと、
（ｉｉ）前記心内信号の群を選択することと、
（ｉｉｉ）前記群内で、逸脱の既定の尺度よりも統計的に逸脱している１つ以上のアノテーション値を特定することと、
（ｉｖ）前記アノテーション値をユーザに対して可視化する一方で、前記統計的に逸脱しているアノテーション値を省略して、可視化するのを防止することと、
を実行することによって、実施し、
前記アノテーション可視化操作のうちの１つ以上にわたって、アノテーション値の省略率を評価し、
前記省略率が既定の閾値を超えることを検出したことに応じて、
省略されていない前記アノテーション値と、省略された前記アノテーション値のうちの１つ以上を用いて、前記統計的に逸脱しているアノテーション値を再び特定する、
ことによって修正措置を取る、
ように構成されているプロセッサと、
を備えるシステム。
【請求項2】
前記プロセッサが、前記逸脱の前記尺度を、前記アノテーション値の標準スコアを用いて定義するように構成されている、請求項１に記載のシステム。
【請求項3】
前記プロセッサが、前記逸脱の前記尺度を、前記アノテーション値の１つ以上の百分位数を用いて定義するように構成されている、請求項１に記載のシステム。
【請求項4】
所与のアノテーション可視化操作において、前記プロセッサが、前記心臓内で互いから既定の距離以内に位置する、空間的に関連する電極の選択されたサブセットによって取得された心内信号に対する前記アノテーション値の逸脱を、計算するように構成されている、請求項１に記載のシステム。
【請求項5】
所与のアノテーション可視化操作に対する前記アノテーション値の逸脱を計算する際に、前記プロセッサが、既定の時間期間内に発生し、互いに時間的に近接する複数の心周期にわたって、前記心内信号の平均をとるように構成されている、請求項１に記載のシステム。
【請求項6】
所与のアノテーション可視化操作において、前記プロセッサが、前記群内の所与の電極によって取得された所与の心内信号内の前記アノテーション値のうちの１つ以上を修正して、前記群内のその他の電極に対する前記所与の電極の変位を補償するように構成されている、請求項１に記載のシステム。
【請求項7】
前記アノテーション値が、局所的興奮時間（ＬＡＴ）を含む、請求項１に記載のシステム。
【請求項8】
前記プロセッサが、前記統計的に逸脱しているアノテーション値以外の前記アノテーション値を、前記心臓のモデル上に重ね合わせることで、前記アノテーション値を可視化するように構成されている、請求項１に記載のシステム。
【請求項9】
プロセッサの作動方法であって、
前記プロセッサが、患者の心臓内の心臓内プローブの複数の電極によって取得された、複数の心内信号を受信することと、
後続の時点で、一連のアノテーション可視化操作を、各操作において、
（ｉ）前記プロセッサが、前記心内信号から複数のアノテーション値を抽出することと、
（ｉｉ）前記プロセッサが、前記心内信号の群を選択することと、
（ｉｉｉ）前記プロセッサが、前記群内で、逸脱の既定の尺度よりも統計的に逸脱している１つ以上のアノテーション値を特定することと、
（ｉｖ）前記プロセッサが、前記アノテーション値をユーザに対して可視化する一方で、前記プロセッサが、前記統計的に逸脱しているアノテーション値を省略して、可視化するのを防止することと、
を実行することによって、実施することと、
前記アノテーション可視化操作のうちの１つ以上にわたって、前記プロセッサが、アノテーション値の省略率を評価することと、
前記省略率が既定の閾値を超えることを検出したことに応じて、
前記プロセッサが、省略されていない前記アノテーション値と、省略された前記アノテーション値のうちの１つ以上を用いて、前記統計的に逸脱しているアノテーション値を再び特定することと、
によって修正措置を取ることと、
を含む作動方法。
【請求項10】
前記逸脱の前記尺度が、前記アノテーション値の標準スコアを用いて定義される、請求項９に記載の作動方法。
【請求項11】
前記逸脱の前記尺度が、前記アノテーション値の１つ以上の百分位数を用いて定義される、請求項９に記載の作動方法。
【請求項12】
所与のアノテーション可視化操作において、前記心臓内で互いから既定の距離以内に位置する、前記プロセッサが、空間的に関連する電極の選択されたサブセットによって取得された心内信号に対する前記アノテーション値の逸脱を計算することを含む、請求項９に記載の作動方法。
【請求項13】
所与のアノテーション可視化操作に対する前記アノテーション値の逸脱を計算する際に、既定の時間期間内に発生し、互いに時間的に近接する複数の心周期にわたって、前記プロセッサが、前記心内信号の平均をとることを含む、請求項９に記載の作動方法。
【請求項14】
所与のアノテーション可視化操作において、前記プロセッサが、前記群内の所与の電極によって取得された所与の心内信号内の前記アノテーション値のうちの１つ以上を修正して、前記群内のその他の電極に対する前記所与の電極の変位を補償することを含む、請求項９に記載の作動方法。
【請求項15】
前記アノテーション値が、局所的興奮時間（ＬＡＴ）を含む、請求項９に記載の作動方法。
【請求項16】
前記プロセッサが、前記アノテーション値を可視化することが、前記プロセッサが、前記統計的に逸脱しているアノテーション値以外の前記アノテーション値を、前記心臓のモデル上に重ね合わせることを含む、請求項９に記載の作動方法。
【発明の詳細な説明】
【技術分野】
【0001】
（関連出願の相互参照）
本出願は、同日付で出願された「Ｕｓｉｎｇ Ｓｔａｔｉｓｔｉｃａｌ Ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ Ｍｕｌｔｉｐｌｅ Ｇｒｏｕｐｅｄ ＥＣＧ Ｓｉｇｎａｌｓ ｔｏ Ｄｅｔｅｃｔ Ｉｎｃｏｎｓｉｓｔｅｎｔ Ｓｉｇｎａｌｓ」と題された米国特許出願（代理人整理番号第ＢＩＯ６１９６ＵＳＮＰ１号）に関連し、その開示は参照により本明細書に組み込まれる。
本出願は、同日付で出願された「Ｕｓｉｎｇ Ｓｔａｔｉｓｔｉｃａｌ Ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ Ｍｕｌｔｉｐｌｅ Ｇｒｏｕｐｅｄ ＥＣＧ Ｓｉｇｎａｌｓ ｔｏ Ｄｅｔｅｃｔ Ｉｎｃｏｎｓｉｓｔｅｎｔ Ｓｉｇｎａｌｓ」と題された米国特許出願（代理人整理番号第ＢＩＯ６１９６ＵＳＮＰ１号）に関連し、その開示は参照により本明細書に組み込まれる。
【0002】
（発明の分野）
本発明は、概して、体内医療処置及び器具に関し、特に心臓の心電図（ＥＣＧ）の体内検知に関する。
本発明は、概して、体内医療処置及び器具に関し、特に心臓の心電図（ＥＣＧ）の体内検知に関する。
【背景技術】
【0003】
多数の電極によって生成される心内心電図（ｉＥＣＧ）信号を測定及びアノテーションするとき、信号に埋め込まれたノイズを低減するために、信号を処理（例えば、コンピュータ処理）することが望ましい場合がある。
【0004】
かかるｉＥＣＧ信号処理のための様々な方法が存在する。例えば、米国特許出願公開第２００９／００８９０４８号には、心室チャネル、マッピングチャネル、及び複数の基準チャネルを含む、４本以上のチャネルによる多チャネル心電図信号から、局所的興奮時間（local activation time＝ＬＡＴ）を判定するための自動化された方法が記載されている。
【発明の概要】
【課題を解決するための手段】
【0005】
本明細書に記載の本発明の一実施形態は、信号取得回路及びプロセッサを含むシステムを提供する。信号取得回路は、患者の心臓内の心臓内プローブの複数の電極によって取得された、複数の心内信号を受信するように構成される。プロセッサは、一連のアノテーション可視化操作を、後続の時点で実行するように構成されるが、そのためには、各動作において、心内信号から複数のアノテーション値を抽出し、心内信号の群を選択し、その群から、既定の逸脱の尺度を超えて統計的に逸脱している１つ以上のアノテーション値を特定し、統計的に逸脱しているアノテーション値を省略して、可視化しないようにする一方で、アノテーション値をユーザに可視化するようになっている。プロセッサは、アノテーション可視化操作のうちの１つ以上にわたって、アノテーション値が省略される割合を評価して、省略される割合が既定の閾値を超えることを検出したことに応じて修正措置を取るように、更に構成されている。
【0006】
一部の実施形態では、プロセッサが、省略されたアノテーション値のうちの１つ以上を再適用し、統計的に逸脱しているアノテーション値を再び特定することによって、修正措置を取るように構成されている。一部の実施形態では、プロセッサが、心内信号からアノテーション値のうちの１つ以上を再抽出することによって、修正措置を行うように構成されている。
【0007】
一実施形態では、プロセッサが、逸脱の尺度を、アノテーション値の標準スコアを用いて定義するように構成されている。別の実施形態では、プロセッサが、逸脱の尺度を、アノテーション値の１つ以上の百分位数を用いて定義するように構成されている。
【0008】
例示的実施形態では、所与のアノテーション可視化操作において、プロセッサが、心臓内で互いから既定の距離以内に位置する、空間的に関連する電極の選択されたサブセットによって取得された心内信号に対するアノテーション値の逸脱を、計算するように構成されている。別の実施形態では、所与のアノテーション可視化操作に対するアノテーション値の逸脱を計算する際に、プロセッサが、既定の時間期間内に発生する複数の時間的に関連する心周期にわたって、心内信号の平均をとるように構成されている。
【0009】
更に別の実施形態では、所与のアノテーション可視化操作において、プロセッサが、群内の所与の電極によって取得された所与の心内信号内のアノテーション値のうちの１つ以上を修正して、群内のその他の電極に対する所与の電極の変位を補償するように構成されている。
【0010】
一部の実施形態では、アノテーション値は、局所的興奮時間（ＬＡＴ）を含む。一部の実施形態では、プロセッサが、統計的に逸脱しているアノテーション値以外のアノテーション値を、心臓のモデル上に重ね合わせることで、アノテーション値を可視化するように構成されている。
【0011】
本発明の一実施形態によれば、患者の心臓内の心臓内プローブの複数の電極によって取得された、複数の心内信号を受信することを含む方法が、更に提供される。後続の時点で、一連のアノテーション可視化操作を、各操作において、（ｉ）心内信号から複数のアノテーション値を抽出することと、（ｉｉ）心内信号の群を選択することと、（ｉｉｉ）群内で、逸脱の既定の尺度よりも統計的に逸脱している１つ以上のアノテーション値を特定することと、（ｉｖ）アノテーション値をユーザに対して可視化する一方で、統計的に逸脱しているアノテーション値を省略して、可視化するのを防止することとを実行することによって、実施する。アノテーション値の省略率は、アノテーション可視化操作のうちの１つ以上にわたって評価される。省略率が既定の閾値を超えていることを検出したことに応じて、修正措置が取られる。
【0012】
本発明は、以下の「発明を実施するための形態」を図面と併せて考慮することで、より完全に理解されよう。
【図面の簡単な説明】
【0013】
【図1】本発明の一実施形態による、心内ＥＣＧ信号の多チャネル測定のための電気解剖学的システムの概略描写図である。
【図2】本発明の一実施形態による、複数の心周期における複数の電極による信号の取得を、概略的に示す図である。
【図3A】本発明の一実施形態による、アノテーション値の信頼性を高めるための第１の方法を模式的に示すフローチャートである。
【図3B】本発明の一実施形態による、アノテーション値の信頼性を高めるための第２の方法を模式的に示すフローチャートである。
【図4】本発明の一実施形態による、アノテーション値の信頼性を高めるための改良された方法を概略的に示すフローチャートである。
【図5A】本発明の実施形態による、電気解剖学的システム内の空間的に関連する電極群からの、一連のＬＡＴ値のアノテーションを、概略的に示す図である。
【図5C】本発明の実施形態による、省略値を再適用することによる誤り伝播の回避を、概略的に示す図である。
【図6A】本発明の実施形態による、電気解剖学的システム内の空間的に関連する電極群によって取得される一連の信号のＬＡＴ値の確率を、概略的に示す図である。
【図6B】本発明の一実施形態による、群－１の信号に対するＬＡＴ値のアノテーションを、概略的に示す図である。
【図6C】本発明の一実施形態による、群－２の信号に対するＬＡＴ値のアノテーションを、概略的に示す図である。
【図6D】本発明の実施形態による、群－１の信号の再アノテーションを、概略的に示す図である。
【図6E】本発明の一実施形態による、再アノテーション後に行われる、群－２のアノテーション値の再計算を、概略的に示す図である。
【図7】本発明の実施形態による、再アノテーションによる、誤り伝播の回避方法を、概略的に示すフローチャートである。
【発明を実施するための形態】
【0014】
概論
心臓内プローブベース（例えば、カテーテルベース）の心臓診断及び治療システムは、侵襲的処置中に心電図（ＥＣＧ）などの複数の心内信号を測定し得る。このようなシステムは、プローブの遠位端に取り付けられた電極（以下、「遠位電極」とも呼ばれる）を使用して、複数の心内信号を取得し得る。測定された信号を用いて、患者の心臓内の病理学的電気パターンの発生源の３Ｄマッピングなどの、心臓の視覚的情報を医師に提供し、アブレーションなどの矯正的医療処置を支援することができる。
心臓内プローブベース（例えば、カテーテルベース）の心臓診断及び治療システムは、侵襲的処置中に心電図（ＥＣＧ）などの複数の心内信号を測定し得る。このようなシステムは、プローブの遠位端に取り付けられた電極（以下、「遠位電極」とも呼ばれる）を使用して、複数の心内信号を取得し得る。測定された信号を用いて、患者の心臓内の病理学的電気パターンの発生源の３Ｄマッピングなどの、心臓の視覚的情報を医師に提供し、アブレーションなどの矯正的医療処置を支援することができる。
【0015】
測定された信号は、典型的には弱い信号で、信号対雑音比（ＳＮＲ）が低い信号である。更に、一部の電極と組織とのガルバニック接続は、不十分であるか、又は存在しなくてもよい。一方、多くの電極が使用されるため、システムが電極から受信するデータにはある程度の冗長性が存在し得る。
【0016】
本明細書に開示される本発明の実施形態は、遠位電極が収集する信号の統計的特性を使用して、収集されたデータの品質及び信頼性を改善する、心臓内プローブベースの電気解剖学的測定及び分析システム及び方法を提供する。
【0017】
以下の説明において、本明細書では、局所的興奮時間（ＬＡＴ）のアノテーション値に言及する。しかしながら、開示される技術は、ＬＡＴに限定されず、本発明の様々な実施形態では、様々な他の好適な信号パラメータのアノテーション値を使用し得る。
【0018】
本発明による一部の実施形態では、プロセッサは、信号のアノテーション値（例えば、ＬＡＴ）を抽出し、次いで、対応する電極群（電極の全て又は一部を含み得る）によって取得される信号群のＬＡＴ値の統計的特性を計算する。一実施形態では、統計的特性は、信号群のＬＡＴ値の平均（例えば、
【0019】
【数1】
【0020】
【数2】
【0021】
別の実施形態では、統計的特性は、ＬＡＴ値群の四分位値を含む。プロセッサは、第１の四分位値Ｑ１及び第３の四分位値Ｑ３を計算し、次いで、Ｑ１未満又はＱ３超の全ての値を無視する（なお、第１の四分位値（Ｑ１）は、データセットの最小の数と中央値との間の中間の数として定義され、第３の四分位値（Ｑ３）は、データセットの中央値と最高値との間の中間の値である）。あるいは、プロセッサは、ＬＡＴ値の逸脱の尺度を、ＬＡＴ値の任意の他の好適な百分位数（又は複数の百分位数）によって定義してもよい。更に別の方法としては、外れＬＡＴ値を廃棄する任意の他の好適な処理法を使用することができる。
【0022】
本明細書に開示される技術は、ノイズがないこと及びガルバニック接続が不良でなく、そのため各群の電極どうしは、互いに類似のアノテーション値を示すものと想定している。通常、電極どうしが互いに遠く離れている場合には、それらの電極によって取得されるアノテーション値どうしは、実質的に異なる値をとり得る。しかも、各電極からの信号は、心拍毎に（「心周期毎に」）、周期的にアノテーションされてもよく、心周期どうしが互いに時間的に遠く離れている場合には、それらの心周期から導出されるアノテーション値どうしは互いに異なる値をとり得る。一実施形態では、群内の信号は相互に関係している。一部の実施形態では、追跡システムが、電極の幾何学的位置を測定し、群は、互いに隣接する電極（「空間的に関連する」、すなわち、互いに所定の距離以内に位置する電極）のみから導出されるアノテーション値を含む。他の実施形態では、群は、互いに隣接する心周期（「時間的に関連する」、すなわち、全てが既定の継続した時間以内に発生する心周期）のみから導出されるアノテーション値を含み、また、一実施形態では、群は、空間及び時間の両面で互いに関連する値を含む（「関連値」と略記される）。
【0023】
一部の実施形態では、プロセッサは、互いに関連するＬＡＴ値の群の統計的特性を計算した後、群内で、統計的に逸脱している、例えば、群内の値の平均値と実質的に異なるＬＡＴ値を省く（残ったＬＡＴ値の群は、有効なＬＡＴ値の群と呼ぶ）。したがって、接続不良の電極に対応するＬＡＴ値、又は極度のノイズを受ける電極に対応するＬＡＴ値は、有効なＬＡＴ値の群から排除され得る。
【0024】
実施形態では、ＬＡＴ値が、一群の信号の平均ＬＡＴから統計的に逸脱しているか否かを判定するために、プロセッサは、アノテーションされたＬＡＴ値の、ＬＡＴ値群の平均値からの偏差を測定する。一実施形態では、逸脱の尺度は、値の標準スコア（当該値と平均との差を、標準偏差で割ったものとして定義される）であり、そのスコアは、事前に設定された限界値と比較される。例えば、平均値よりも、標準偏差３．５超大きい値（標準スコア＝３．５）、又は平均値よりも、標準偏差１．５超小さい値（標準スコア＝－１．５）は、統計的に逸脱していると考えられてもよく、したがって省かれてもよい。別の実施形態では、プロセッサは、第１の四分位値未満の値又は第３四分位値超の値を省く。
【0025】
本発明の一部の実施形態では、プロセッサは、心臓内での心臓信号伝播の異なる時間遅延に起因して、空間的に関連する電極のＬＡＴ値どうしの不一致を緩和し得る。実施形態によれば、プロセッサは、伝播遅延の差を相殺するために、所与の電極によって取得されたＬＡＴアノテーションを、他の電極に対する所与の電極の変位を補償することによって修正してもよい。
【0026】
心電図信号は、時には曖昧に解釈されてもよく、２つのＬＡＴ値は、異なる確率で、アノテーションされてもよい（正しい値及び誤った値）。例えば、信号は、互いに近接する２つのピークを有してもよい。その結果、一部の場合においては、サブセットの計算された平均値は、正しいＬＡＴ値よりも誤ったＬＡＴ値に近いことがあり得る。
【0027】
一部の実施形態では、プロセッサは、ある時間にわたって群の選択を繰り返す（例えば、３２回の心周期毎に１回）が、その際に後続の群どうしの間には部分的な重複がある。プロセッサは、古い値（例えば、最も古い１６心周期で抽出された値）の一部（ただし全てではない）を除去し、新たな値（例えば、最も新しい１６心周期で抽出された値）を追加することによって、新しい群を形成する（この文脈において、「古い」及び「新しい」は、信号が獲得された心周期の連続する番号を指す）。プロセッサが計算する第１の群の平均が誤っている場合（例えば、信号の曖昧さによる）には、プロセッサは、次の群から、重複している正しいＬＡＴ値を省略し、平均計算を中断し得るが、そのため、第１の平均計算における誤りが、後の群に伝播し得るものの、新たにアノテーションされたＬＡＴ値が正しいものとなり得る。
【0028】
本明細書に開示される本発明による実施形態は、そのような誤り伝播を回避する。一実施形態では、プロセッサは、省略されたＬＡＴ値の数を監視し、値が省略される割合に応じて、プロセッサは、新たな平均の計算において省略される値を考慮することを決定し、誤りが誤りのある群を超えて伝播しないようにしてもよい。他の実施形態では、プロセッサが多数の省略された値を検出する場合、プロセッサは、対応する信号を代替的に解釈することによって値の一部を再びアノテーションして、続いて省略値を再度カウントしてもよい。プロセッサは、省略された値の数が減少する場合には、代替的アノテーションを選択する。
【0029】
要約すると、本発明の実施形態によるプロセッサは、空間的及び／又は時間的に関連する心臓間信号のアノテーション値の群の品質及び信頼性を改善することができるが、そのためには、アノテーション値の統計的特性を計算し、アノテーション値を群の平均値と比較し、有効な値の群から、平均値から離れた値を除外する。一部の実施形態では、統計的特性を計算する前に、プロセッサは、アノテーション値の群を変更して、信号の伝播遅延を修正してもよい。誤り伝播を回避するために、プロセッサは、省略された値の数を監視して、省略が発生する割合に応じて、当該群を、省略された値を元に戻した状態で再評価してもよい。他の実施形態では、省略された値の割合に応じて、プロセッサは、値を再びアノテーションして、代替的なＬＡＴ解釈を探してもよい。
【0030】
システムの説明
図１は、本発明の一実施形態による、心内ＥＣＧ信号の多チャネル測定のための電気解剖学的システム２１の概略描写図である。一部の実施形態では、システム２１は、心臓の電気解剖学的マッピングに使用される。
図１は、本発明の一実施形態による、心内ＥＣＧ信号の多チャネル測定のための電気解剖学的システム２１の概略描写図である。一部の実施形態では、システム２１は、心臓の電気解剖学的マッピングに使用される。
【0031】
図１は、医師２２が電気解剖学的カテーテル２３を使用して、患者２５の心臓２４の電気解剖学的マッピングを実施しているのを描写している。カテーテル２３は、その遠位端に、機械的に柔軟性であり得る１本以上のアーム２６を備え、それらのアームのそれぞれに、１つ以上の遠位電極２７が連結されている。理解されるように、図１は、５本のアームを有するカテーテルを示しているが、他の種類のカテーテルを、本発明による代替的な実施形態で使用してもよい。電極は、インターフェース３２を介して、プロセッサ３４に連結される。
【0032】
電気解剖学的マッピング処置中に、トラッキングシステムを使用して、遠位電極２７の心臓内位置を追跡し、それによって、取得された電気生理学的信号のそれぞれを、既知の心臓内位置と関連付け得る。トラッキングシステムの一例が、米国特許第８，４５６，１８２号に記載されている、能動的現在位置特定（Active Current Location＝ＡＣＬ）である。ＡＣＬシステムでは、プロセッサは、遠位電極２７のそれぞれと患者２５の皮膚に連結されている複数の表面電極２８との間で測定されたインピーダンスに基づいて、遠位電極のそれぞれの位置を推定する（図示を容易にするため、図１には１つの表面電極のみが示されている）。次いでプロセッサは、遠位電極２７から受信した任意の電気生理学的信号と、信号が取得された位置とを関係付け得る。
【0033】
一部の実施形態では、複数の遠位電極２７は、心臓２４の心房心室の組織から心内ＥＣＧ信号を取得する。プロセッサは、インターフェース３２から心内信号を受信するように連結された信号取得回路３６と、データ及び／又は命令を記憶するためのメモリ３８と、処理ユニット４２（例えば、ＣＰＵ又は他のプロセッサ）とを備える。
【0034】
信号取得回路３６は、心内信号をデジタル化して、複数のデジタル信号を生成する。次いで、取得回路は、デジタル化された信号をプロセッサ３４に含まれる処理ユニット４２に伝達する。
【0035】
処理ユニット４２は、他のタスクの中でもとりわけ、信号からアノテーションパラメータを抽出し、互いに類似している可能性がある隣接する信号群のアノテーションされたパラメータの平均値などの、統計的特性を計算するように構成されている。（なお、この文脈では、隣接する信号とは、互いに近接して位置する電極からの信号（「空間的に関連する」）、及び／又は時間的に互いに近接する心周期から抽出されたアノテーション値（「時間的に関連する」）を指す）。
【0036】
処理ユニットは、統計的特性を計算した後、群から無効である可能性が高いアノテーション値（例えば、ガルバニック接続が不十分であるか又は強力な時間的ノイズを受ける電極からのアノテーション値など）を除外する（すなわち、省略する）ように更に構成される。残りのアノテーション値は、これ以降、「有効アノテーション値」と呼ぶものとする。
【0037】
プロセッサ３４は、有効アノテーション値、すなわち、省略された統計的に逸脱しないアノテーション値を除外したアノテーション値を、ユーザに対して可視化する。一部の実施形態では、プロセッサ３４は、例えば、有効アノテーション値を心臓の電気解剖学的マップ５０上に重ね合わせて、そのマップを医師２２に対して、画面５２上で表示することで、有効アノテーション値を可視化する。あるいは、プロセッサ３４は、任意の他の好適な方法で、有効アノテーション値（無効なアノテーション値を省略した後の値）を可視化してもよい。
【0038】
図１に示す例示的な図は、純粋に、概念を分かりやすくするために選択されたものである。例えば、本発明の代替的な実施形態では、位置測定は、表面電極２８の対どうしの間に電圧勾配を印加し、得られた電位を、遠位電極２７を用いて測定することによって（すなわち、Ｂｉｏｓｅｎｓｅ－Ｗｅｂｓｔｅｒ社（カリフォルニア州、Ｉｒｖｉｎｅ）により製造されたＣＡＲＴＯ（登録商標）４の技術を使用して）も行われ得る。したがって、本発明の実施形態は、任意の位置検知方法に適用される。
【0039】
Ｌａｓｓｏ（登録商標）Ｃａｔｈｅｔｅｒ（Ｂｉｏｓｅｎｓｅ－Ｗｅｂｓｔｅｒ製）又はバスケットカテーテルなどの他の種類のカテーテルが、同等に使用されてもよい。接触センサは、電気解剖学的カテーテル２３の遠位端に取り付けられてもよい。例えばアブレーションに使用されるもののような他の種類の電極を、電極２７と同様に使用して、心内の電気生理学的信号を取得してもよい。
【0040】
図１は、本発明の実施形態に関連する部分を主に示す。外部ＥＣＧ記録電極及びそれらの接続部などの、他のシステム要素は省略されている。様々なＥＣＧ記録システム要素も、フィルタリング、デジタル化、回路保護、及びその他の要素も省略されている。
【0041】
任意の実施形態では、心内ＥＣＧ信号を測定するために、読み出し用の、特定用途向け集積回路（ＡＳＩＣ）が使用される。信号取得回路３６をルーティングする様々な要素は、例えば、１つ以上の個別の構成要素、フィールドプログラマブルゲートアレイ（ＦＰＧＡ）又はＡＳＩＣなどを使用してハードウェアに実装されてもよい。一部の実施形態では、信号取得回路３６及び／又は処理ユニット４２のいくつかの要素は、ソフトウェアに実装されてもよく、又はソフトウェア要素及びハードウェア要素の組み合わせを使用して実装されてもよい。
【0042】
処理ユニット４２は、通常、本明細書に記載されている機能を実行するようにプログラムされたソフトウェアを有する汎用プロセッサを備える。ソフトウェアは、例えば、ネットワーク上で、電子形態でダウンロードすることができるか、あるいはそれは、代替的に又は追加的に、例えば、磁気メモリ、光学メモリ、又は電子メモリなどの、非一時的な有形媒体で提供でき、かつ／又はそのような媒体に記憶できる。
【0043】
関連するアノテーション値
関連するアノテーション値は、空間的に関連する電極（例えば、互いに幾何学的に近接している電極、すなわち、互いから、既定の距離以下だけ離れて位置している電極）から、かつ／又は時間的に関連する信号（例えば、互いに近接している心周期から抽出された値、すなわち、既定の継続時間以内で発生する値）から導出される。より正確には、関連するアノテーション値とは、電極間の幾何学的距離及び心周期間の時間的距離を含む、組み合わされた距離が、なんらかの既定の閾値を下回るアノテーション値である。
関連するアノテーション値は、空間的に関連する電極（例えば、互いに幾何学的に近接している電極、すなわち、互いから、既定の距離以下だけ離れて位置している電極）から、かつ／又は時間的に関連する信号（例えば、互いに近接している心周期から抽出された値、すなわち、既定の継続時間以内で発生する値）から導出される。より正確には、関連するアノテーション値とは、電極間の幾何学的距離及び心周期間の時間的距離を含む、組み合わされた距離が、なんらかの既定の閾値を下回るアノテーション値である。
【0044】
図２は、複数の心周期における複数の電極による信号の取得を、概略的に示す図２００である。水平軸２０２は、心周期（各垂直線は１つの心周期である）を示し、垂直軸２０４は、基準点からの電極の距離（１つの空間的次元のみが示されている）を示す。（理解されるように、２つ又は３つの次元が実際に使用されてもよいが、明確にするために図示されていない）。図２に示す例示的な実施形態によれば、全ての水平線に電極が存在し、ＬＡＴアノテーション値が、水平線と垂直線との全ての交点に対して登録されている（各交点は、これ以降、ＬＡＴ点と呼ぶものとする）。
【0045】
曲線２０６は、示されたＬＡＴ値の位置を示す等ＬＡＴ線であり、電極は、対応する心周期において、隣接する等ＬＡＴ曲線から補間された値を測定する可能性が高い。例えば、ＬＡＴ点２０８の期待される登録値（垂直方向で、等ＬＡＴ線７１４と等ＬＡＴ線７１６との間の途中にある）が７１５であるのに対し、ＬＡＴ点２１０の期待される登録値は７０８．５である。
【0046】
理解されるように、隣接する垂直線どうしのＬＡＴ値どうし、及び隣接する水平線どうしのＬＡＴ値どうしは、互いに似た値である。円２１２は、幾何学的（垂直方向）及び時間的（水平方向）距離に関して、互いに近くにある、関連するＬＡＴ値２１４の群を表す。
【0047】
図２に示されている例示的な図は、単に概念を分かりやすくする目的のために簡略化されて図示されている。例えば、代替的な実施形態では、電極間の距離は均一ではなく、関連する信号の群は円でなくてもよい。
【0048】
図３Ａは、本発明の実施形態による、アノテーション値の信頼性を高めるための第１の方法を、概略的に示すフローチャート３００である。フローは、プロセッサ３４（図１）によって実行される。フローは、信号記録工程３０２で始まるが、その工程では、プロセッサが、電極２７によって監視され、取得回路３６によって取得されたＥＣＧ信号を記録する（図１）。次に、アノテーション値抽出工程３０４において、プロセッサは、各電極及び各心周期のアノテーション値を計算する。
【0049】
次いで、プロセッサは、電極位置取得工程３０６に移行し、その工程では、電極の位置が、（例えば、ＡＣＬ技術を使用して）取得され、かつ各電極の空間位置が登録される。次いで、プロセッサは、群選択工程３０８に移行する。
【0050】
工程３０８では、プロセッサは、関連するアノテーション値の群を選択する。既に説明したように、群は、空間的及び／又は時間的に関連する信号からの、互いに類似している可能性が高いアノテーション値を含む。
【0051】
次に、平均値及び標準偏差計算工程３１０では、プロセッサが、その群の全てのアノテーション値の平均値と標準偏差とを計算する。この文脈では、例えば、算術平均、幾何平均、中央値、２乗平均平方根値（ＲＭＳ）値、質量中心、又は任意の他の平均値など、任意の好適な種類の平均値を使用することができる。
【0052】
次いで、プロセッサは、その群の各アノテーション値に対して繰り返し、工程３１２と、工程３１４と、工程３１６又は工程３１８のいずれかとに、順次移行する。標準スコア計算工程３１２では、プロセッサが、アノテーション値の標準スコアを（例えば、アノテーション値と平均との差を、標準偏差で除算することによって）計算する。標準スコア比較工程３１４では、プロセッサが、工程３１２で計算された標準スコアを、事前に設定された限界値と比較する。標準スコアが、事前に設定された限界値を超える場合には、値除外工程３１６に移行し、その工程では、プロセッサが、統計的に逸脱しているアノテーション値を除外する。標準スコアが、事前に設定された限界値の範囲内にある場合には、値追加工程３１８に移行し、その工程では、プロセッサが、アノテーション値を有効アノテーション値の群に追加する。
【0053】
プロセッサは、その群の全てのアノテーション値について、工程３１２と、工程３１４と、工程３１６又は工程３１８のいずれかとのシーケンスを繰り返す。次いで、フローチャートは、関連する電極の他の群について、（工程３０８から）工程を繰り返してよい。
【0054】
フローが終了すると、有効アノテーション値の群が、元の群と置き換わっており、その新たな群では、極端な値（例えば、ガルバニック接続が不十分な電極からの値）が省略されているため、その信頼性がより高くなっている。
【0055】
図３Ｂは、本発明の実施形態による、アノテーション値の信頼性を高めるための第２の方法を、概略的に示すフローチャート３５０である。図３Ｂに示す方法は、統計的特性及び省略される値の選択においてのみ、図３Ａに示す方法と異なっている。したがって、図３Ａに示す工程３０２～３１８は、それぞれ図３Ｂの工程３５２～３６８とそれぞれ同一である。ただし、例外は、図３Ｂの工程３６０及び３６４であり、これらは、図３Ａの工程３１０及び３１４とそれぞれ異なっており、それゆえ、以下に説明する。
【0056】
四分位値計算工程３６２では、処理ユニット４２（図１）が、ＬＡＴ値群の第１の四分位値（Ｑ１）及び第３の四分位値（Ｑ３）を計算する（なお、Ｑ１は、最小数とＬＡＴ値群の中央値との間の中間の数として定義され、Ｑ３は、中央値とＬＡＴ値群の最大値との間の中間の数である）。
【0057】
アノテーション値比較工程３６４では、処理ユニットが、アノテーションされたＬＡＴ値を、Ｑ１及びＱ３とそれぞれ比較する。その値が、Ｑ１よりも小さいか又はＱ３より大きいという場合には、処理ユニットは、アノテーション値除外工程３６６に移行し、その値がＱ１とＱ３との間にある場合には、処理ユニットは、値追加工程３６８に移行する。
【0058】
図３Ａ及び図３Ｂに示される例示的なフローチャートは、純粋に概念を明確化する目的で選択されている。例えば、代替的な実施形態では、（信号が記録された後ではなく）信号が取得されると、アノテーション値が抽出されてもよい。一実施形態では、その群の信号の選択は、医師によって行われてもよく、他の実施形態では、プロセッサが、医師が指示する空間範囲及び／又は時間範囲に従って、その群を選択する。
【0059】
一部の実施形態では、工程３１８（図３Ｂでは工程３６８）は必要ではなく、プロセッサは、工程３１６（同３６６）において、その群から極端な値を除外して、フローが完了したときには、良好な値のみが残っている。他の実施形態では、全てのアノテーション値は、最初は無効としてマークされ、工程３１６（３６６）は必要ではない。
【0060】
一部の実施形態では、使用される他の統計的特性は、既に説明したものとは異なる。例えば、一実施形態では、四分位値ではなく八分位値を使用することができ、処理ユニットは、第１の八分位値よりも小さい値又は最後の八分位値よりも大きい値を省略してよい。更に代替的に、任意の他の好適な百分位数が用いられ得る。
【0061】
極端な値を検出しかつ省略するための任意の他の好適な統計的方法を、代替的実施形態では使用してよい。
【0062】
伝播遅延補償
一部の実施形態では、既に説明した技術が、電極の異なる空間的位置による値の予想される変化のために、抽出されたＬＡＴ値を、統計的特性計算の前に修正することによって改善され得る。例えば、心臓を通る波は、所与の速度（例えば、１ｍ／秒）で移動すると想定することができる。信号を取得する電極の既知の位置を使用して、ＬＡＴ値の理論的な差を、平均値を計算する際に適用することができる。
一部の実施形態では、既に説明した技術が、電極の異なる空間的位置による値の予想される変化のために、抽出されたＬＡＴ値を、統計的特性計算の前に修正することによって改善され得る。例えば、心臓を通る波は、所与の速度（例えば、１ｍ／秒）で移動すると想定することができる。信号を取得する電極の既知の位置を使用して、ＬＡＴ値の理論的な差を、平均値を計算する際に適用することができる。
【0063】
図４は、本発明の実施形態による、アノテーション値の信頼性を高めるための改善された方法を、概略的に示すフローチャート４００である。フローは、プロセッサ３４（図１）によって実行される。フローは、信号記録工程４０２で開始され、アノテーション値計算工程４０４がそれに続き、更に電極位置取得工程４０６がそれに続き、また更に群選択工程４０８がそれに続くが、工程４０２、４０４、４０６、及び４０８はそれぞれ、工程３０２、３０４、３０６、及び３０８（図３）と同一であってもよい。
【0064】
次に、プロセッサは、ＬＡＴ値修正工程４１０に移行し、そこでは、その群の各ＬＡＴ値に対して、プロセッサが、電極の空間位置及び想定される波移動速度に従って、推定される修正値を計算し、適用する。工程４１０の後、フローは、図３に、平均値及び標準偏差計算工程３１０で戻る。
【0065】
このように、伝播遅延によって生じる偏差の推定値を、その群から除去することができるので、アノテーション信号の信頼性を更に向上させることができる。
【0066】
図４に示す例示的なフローチャートは、純粋に概念を明確にする目的で選択されている。例えば、代替的な実施形態では、予想される信号遅延の修正は、平均値及び標準偏差計算工程に統合することができる。他の実施形態では、修正は、群が選択される前に行われる（したがって、群は、より多くの関連するＬＡＴ値を含み得る）。
【0067】
誤り伝播及びその防止
本発明による一部の実施形態では、プロセッサ３４が、後続の心拍からの信号の受信に応じて、関連するＬＡＴ値の群を連続的に選択する。場合によっては、プロセッサは、時間的に重複する信号に対応する値の群を選択してもよい。例えば、第１の群が、心周期ｎ～ｍに対応する信号から抽出されたＬＡＴ値を含み、第２の群が、周期ｘ～ｙに対応する値を含むとすると、ｘがｍより大きく（ｘ＞ｍ）かつｙがｎより大きい（ｙ＞ｎ）場合には、群どうしが部分的に重複することになる。
本発明による一部の実施形態では、プロセッサ３４が、後続の心拍からの信号の受信に応じて、関連するＬＡＴ値の群を連続的に選択する。場合によっては、プロセッサは、時間的に重複する信号に対応する値の群を選択してもよい。例えば、第１の群が、心周期ｎ～ｍに対応する信号から抽出されたＬＡＴ値を含み、第２の群が、周期ｘ～ｙに対応する値を含むとすると、ｘがｍより大きく（ｘ＞ｍ）かつｙがｎより大きい（ｙ＞ｎ）場合には、群どうしが部分的に重複することになる。
【0068】
前述の説明においては、図５Ａ～図５Ｃ、図６Ａ～図６Ｅ、及び図７を参照して、異なる心周期からアノテーションされたＬＡＴ値を参照する。ＬＡＴ値のそれぞれは、実際には、複数の空間的に関連するＬＡＴ値（信号伝播遅延について修正されたＬＡＴ値を含む）を含んでもよい。しかしながら、簡潔にするために、心周期ごとに１つのＬＡＴ値を示す。
【0069】
図５Ａ～図５Ｃは、水平時間軸（図示せず）に沿ったイベントを説明するタイミング図である。図５Ａは、本発明の実施形態による、電気解剖学的システムにおける一連のＬＡＴ値のアノテーションを、概略的に示す図５００である。図５Ａに示す例示的な実施形態では、各群は、１０個の連続する心周期に対応する１０個のＬＡＴ値を含み、連続する群のそれぞれの対は、５つの共通のＬＡＴ値を共有する。
【0070】
明確にするために、図５Ａに示される例では、アノテーションされたＬＡＴ値は、「Ｘ」に近いか又は「Ｙ」に近いかのいずれかであり得る（短く、Ｘ値及びＹ値と呼ぶものとする）。理解されるように、この単純化は、断じて本開示の範囲を限定することを意味するものではない。本発明の実施形態は、ＬＡＴ値の任意の好適なセットをアノテーションし得る。
【0071】
図５Ａに示す例では、第１の５つのＬＡＴ値は「Ｘ」であるが、後続の心内信号が受信されるとき、ＬＡＴ値の過半数が「Ｙ」であり、これは、最初のアノテーション及び第１の群について取られた判断は恐らく誤っていたことを意味する。以下に（図５Ｂを参照して）説明されるように、本発明による一部の実施形態では、修正措置は採用されず、第１の群について取られた誤った判断は、更なる群に伝播し得る。更に（図５Ｃを参照して）説明されるように、他の実施形態では、修正措置が取られて、誤り伝播が防止される。
【0072】
図５Ｂは、本発明の実施形態による、修正措置が採用されていないときにプロセッサが（図５Ａで）受信するＬＡＴ値の統計的特性計算を、概略的に示す図５１０である。図５１０は、群－１の平均化スキーム５１２と、群－２の平均化スキーム５１４と、群－３の平均化スキーム５１６とを含む（示されるように、群－４の平均化は群－２及び群－３と同一であり、図示されていない）。
【0073】
群－１の平均化スキーム５１２は、最初の１０個のＬＡＴ値の平均化を表す。７つのＸ値及び３つのＹ値が存在するので、平均値はＸに近いことになる。そのため、最後の５つのサンプルのうちのＹ値が、後続の群から省略されることになる。
【0074】
ＬＡＴ値６～１５の平均化を表す群－２の平均化スキーム５１４は、４つのＸ値及び３つのＹ値（２つの元々のＹ値は省略済み）を含む。この平均値は、やはりＸに近いものとなる。そのため、最後の５つのサンプルのうちのＹ値が、後続の群から省略される。理解されるように、第１の群からの誤りが、第２の群に伝播している。
【0075】
ＬＡＴ値１１～２０の平均化を表す群－３の平均化スキーム５１６もまた、２つの元々のＹ値が省略されたので、４つのＸ値及び３つのＹ値を含み、ここでも、平均値はＸに近いものとなり、Ｙ値が省略されることになる。したがって、第１の群からの誤りが、この群に伝播し、更なる群にも続くことになる。（図５Ａ、図５Ｂに示す例によれば、後続のＬＡＴ値の大多数（例えば８０％）がＹであるという場合にのみ伝播が停止する。）
【0076】
図５Ｃは、本発明の実施形態による、省略済みの値を再適用することによって誤り伝播を回避することを、概略的に示す図５２０である。図５２０は、群－１の平均化スキーム５２２と、群－２の予備的平均化スキーム５２４と、群－２の再評価済平均化スキーム５２６と、群－３の平均化スキーム５２８とを含む。
【0077】
群－１の平均化スキーム５２２は、最初の１０個のＬＡＴ値の平均化を表し、群－１の平均化スキーム５１２（図５Ｂ）と同一である。ここでもまた、平均値はＸに近いものとなり、Ｙ値が次の群から省略されることになる。
【0078】
群－２の予備的平均化スキーム５２４は、群－２の平均化スキーム５１４（図５Ｂ）と同一であり、ここでもまた、平均値はＸに近いものとなり、Ｙ値は次の群から省略されることになる。しかしながら、図５Ｃに示す例示的な実施形態によれば、ＬＡＴ値が省略される頻度が監視される。ここで、１０個のサンプル群から６つの省略済の値が存在するため、プロセッサは、統計的特性計算において、省略済みの値を再び適用することを決定する。省略済の値を再適用する決定（例えば、参照符号５２４によって指定された図５Ｃの群２）は、実行されることになる「修正措置」である。
【0079】
群－２の再計算された平均化スキーム５２６は、省略済みのＹ値を、統計的計算において含む。ここで、群－２には、６個のＹ値と４つのＸ値とが存在し、平均値はＹに近いものとなり、Ｘ値が次の群から省略されることになる。ここで、群－３の平均化スキーム５２８は、２つのＸ値と６つのＹ値とを含み、平均値はＹに近いものとなり、Ｘ値が省略されることになる。このパターンは、群－４に続くことになる。
【0080】
したがって、図５Ｃに示す例示的な実施形態によれば、プロセッサは、値が省略される頻度を監視することによって、値の省略に起因して生じ得る誤りの伝播を、停止させることができ、値が省略される頻度に応じて、省略済の値なしで平均値を再計算させることを含み得る修正措置を取ることができる。
【0081】
本発明による他の実施形態は、値が省略される頻度に応じて、ＬＡＴ値の少なくとも一部を再びアノテーションすることによって、誤り伝播を防止することができる。これは、心内信号が曖昧であって、２通り以上に解釈され得るという場合に、有効であり得る。
【0082】
【0083】
図６Ａは、本発明の実施形態による、電気解剖学的システム内の空間的に関連する電極群によって取得される一連の信号のＬＡＴ値の確率を、概略的に示す図である。図６Ａ中の各矩形は、対応する信号のＬＡＴ値が、Ｘである確率（上段の数字）又はＹである確率（下段の数字）を表す。最初の５つの信号が、ＬＡＴ値としてＸ値を表す確率は６０％であり、Ｙ値を表す確率は４０％である。後の信号（後の心周期に関連する信号）は、ほとんどの場合、Ｘ値である確率よりもＹ値である確率が高くなっている。信号は、部分的に重複する群（１～３）に分割され、各群が１０個のサンプルを有しており、しかも、後続の群との間で、５つのサンプルが重複している。
【0084】
図６Ｂは、本発明の一実施形態による、群－１の信号に対するＬＡＴ値のアノテーションを、概略的に示す図である。プロセッサは、サンプルのうちの６つに対してＸ値を、４つに対してＹ値をアノテーションする。平均値はＸに近いものとなり、プロセッサは４つのＹ値を省略する（省略された信号のうちの、点線の矩形によって指定された３つは、群－２と重複する）。
【0085】
図６Ｃは、本発明の一実施形態による、群－２の信号に対するＬＡＴ値のアノテーションを、概略的に示す図である。群－１と共通の５つの信号から、プロセッサは、より高いｙ確率を有する３つの信号を省略済みである。残りの７つの値のうちの４つがＸ値であり、３つがＹ値である。ここでも、平均値はＸに近いものとなり、プロセッサは、次の群からＹ値を省略してもよい。しかしながら、プロセッサは、ここで、省略済の値が多い（例えば、事前に設定された限界値を超える）と判断し、ＬＡＴ値を再アノテーションして、省略される値の数を減らすように試みる。
【0086】
図６Ｄは、本発明の実施形態による、群－１の信号の再アノテーションを、概略的に示す図である。ここで、プロセッサは、最初の５つのサンプルについて、確率の２番目に高い値（信号１～５に対して、それぞれＹ、Ｙ、Ｙ、Ｘ、Ｙ）をアノテーションする。ここで第１の群に対してアノテーションされたＬＡＴ値は、３つのＸ値と７つのＹ値とを含むことになり、平均値はＹに近いものとなる。ここで、第７の値及び第９の値は、次の群から省略されることとなる。
【0087】
図６Ｅは、本発明の一実施形態による、群－１の再アノテーションに続く、群－２の値の再アノテーションを、概略的に示す図である。群－１からの５つの重複信号からの２つのＸ値は省略されていないので、群－２は、６つのＹ値及び２つのＸ値を含む。ここでは平均値がＹに近くなり、誤り伝播が防止される。
【0088】
図５Ａ～図５Ｃ及び図６Ａ～図６Ｅを参照して説明した、誤り伝播防止方法は、単に概念的明確性のために引用された例示的な方法である。本発明による実施形態は、誤り伝播の防止のための代替的な方法を採用することができる。例えば、群の長さ及び重複の長さは、任意の好適な数であってよく、例えば、時間とともに変化する群サイズ及び／又は重複を含み得る。一部の実施形態では、閾値は、動的なものであってもよい。一実施形態では、プロセッサは確率が２番目に高いＬＡＴ値を、その確率が事前に設定された限界値未満である場合には、あるいはＬＡＴ値が省略される頻度の関数である限界値未満である場合には、アノテーションしなくてもよい。
【0089】
最後に、一部の実施形態では、省略済の値を元に戻すことと、アノテーションを再評価することとを組み合わせて用いてもよい。
【0090】
【0091】
フローは、重複アノテーション受信工程７０２で開始するが、そこではプロセッサがメモリから、前の群と現在の群とに共通である心内ＥＣＧ信号のアノテーションを取得する（既に述べたように、プロセッサは、前の群の平均値から実質的に逸脱する値を省略してもよい）。
【0092】
次に、新しいＥＣＧシグナルの受信工程７０４で、プロセッサは、電極から、後続の心周期に関係するＥＣＧ信号を（取得回路３６（図１）を介して）受信し、対応するＬＡＴ値をアノテーションする。次いで、プロセッサは、統計計算工程７０６に移行し、値の群の統計的特性（例えば、平均値及び標準偏差）を計算する。なお、その群には、前の群からの重複値と、新たにアノテーションされた値とを含む。
【0093】
次いで、プロセッサは、群スキャン工程７０８に移行し、そこでは、プロセッサは、その群のＬＡＴ値のそれぞれを群の平均値（工程７０６で計算されたもの）と比較して、事前に設定された閾値を超えて平均値から逸脱した値（又は、あるいは、標準スコアが事前に設定された限界値内でない場合の値）を省略する。プロセッサはまた、省略済のＬＡＴ値の数をカウントする。
【0094】
次に、プロセッサは、カウンタ比較工程７１０に移行し、そこでは、プロセッサは、省略済の値の数を、事前に設定された限界値と比較する。省略済の値の数が限界値を超えない場合には、プロセッサは、次の群に進む（例えば、工程７０２に戻る）。
【0095】
工程７１０で、省略済の値の数が、事前に設定された限界値を超えた場合には、プロセッサは、以前のＬＡＴ値の誤ったアノテーションの結果として誤り伝播が発生したと想定する。次いで、プロセッサは、重複信号受信工程７１２に移行し、メモリから重複信号を取得する。次に、プロセッサは、代替値抽出工程７１４に移行し、信号からＬＡＴ値を再びアノテーションする。プロセッサは、ＬＡＴ値を既に（前の群の工程７０４で）アノテーションしているので、プロセッサは、ここでは、代替的アノテーション（例えば、各信号の２番目に高い確率のＬＡＴ値）を探す。
【0096】
工程７１４の後、プロセッサは工程７０６に再び移行し、統計値を、代替アノテーション値を用いて再計算する。ここで省略済の値の数が、工程７０８での以前のカウント数よりも少ない場合に、新たな値が効力を発することになる。
【0097】
要約すると、図７に示す例示的な実施形態によれば、関連する心内信号の連続する群のＬＡＴ値が、アノテーションされる。第１の群内で誤りが発生する場合、プロセッサは、後続の群内の値が頻繁に省略されているということを検出し得る。次いで、プロセッサは、信号の代替アノテーションを試し、省略されるアノテーションがより少なくなる場合に、代替アノテーションに対応するＬＡＴ値を選択する。
【0098】
図７に示す例示的なフローチャートは、純粋に概念を明確にする目的で選択されている。例えば、代替的な実施形態では、プロセッサは、最も確率の高いアノテーション値の第１のセットと、代替アノテーション値の第２のセットとを抽出し、その結果、工程７１４は、代替値選択工程によって置き換えられる。一部の実施形態では、３つ以上のＬＡＴ値候補が抽出され、フローチャートはそれに応じて変更される。一実施形態では、省略済の値の数が閾値と比較される工程７２０は、カウント数ではなく頻度を比較する他の手段によって置き換えられてもよい。実施形態では、代替アノテーションを抽出する代わりに、プロセッサは、省略済の値を用いて、群の統計的特性を再計算する。更に他の実施形態では、プロセッサは、統計的再計算と代替アノテーション値との両方を試みる。
【0099】
上記の実施形態は例として挙げたものであり、本発明は上記に具体的に示し記載したものに限定されない点が理解されよう。むしろ本発明の範囲は、上で説明される様々な特徴の組み合わせ及びその部分的組み合わせの両方、並びに上述の説明を読むことで当業者に想到されるであろう、従来技術において開示されていないそれらの変形例及び修正例を含むものである。参照により本特許出願に組み込まれる文献は、これらの組み込まれる文献において、いずれかの用語が本明細書において明示的又は暗示的になされた定義と矛盾する様式で定義されている場合には、本明細書における定義のみを考慮するものとする点を除き、本出願の一部とみなすものとする。
【0100】
〔実施の態様〕
（１） 患者の心臓内の心臓内プローブの複数の電極によって取得された、複数の心内信号を受信するように構成されている信号取得回路と、
プロセッサであって、
後続の時点で、一連のアノテーション可視化操作を、各操作において、
（ｉ）前記心内信号から複数のアノテーション値を抽出することと、
（ｉｉ）前記心内信号の群を選択することと、
（ｉｉｉ）前記群内で、逸脱の既定の尺度よりも統計的に逸脱している１つ以上のアノテーション値を特定することと、
（ｉｖ）前記アノテーション値をユーザに対して可視化する一方で、前記統計的に逸脱しているアノテーション値を省略して、可視化するのを防止することと、
を実行することによって、実施し、
前記アノテーション可視化操作のうちの１つ以上にわたって、アノテーション値の省略率を評価し、
前記省略率が既定の閾値を超えることを検出したことに応じて、修正措置を取る
ように構成されているプロセッサと、
を備えるシステム。
（２） 前記プロセッサが、
省略された前記アノテーション値のうちの１つ以上を再適用し、
前記統計的に逸脱しているアノテーション値を再び特定する、
ことによって、前記修正措置を取るように構成されている、実施態様１に記載のシステム。
（３） 前記プロセッサが、前記心内信号から前記アノテーション値のうちの１つ以上を再抽出することによって、前記修正措置を行うように構成されている、実施態様１に記載のシステム。
（４） 前記プロセッサが、前記逸脱の前記尺度を、前記アノテーション値の標準スコアを用いて定義するように構成されている、実施態様１に記載のシステム。
（５） 前記プロセッサが、前記逸脱の前記尺度を、前記アノテーション値の１つ以上の百分位数を用いて定義するように構成されている、実施態様１に記載のシステム。
（１） 患者の心臓内の心臓内プローブの複数の電極によって取得された、複数の心内信号を受信するように構成されている信号取得回路と、
プロセッサであって、
後続の時点で、一連のアノテーション可視化操作を、各操作において、
（ｉ）前記心内信号から複数のアノテーション値を抽出することと、
（ｉｉ）前記心内信号の群を選択することと、
（ｉｉｉ）前記群内で、逸脱の既定の尺度よりも統計的に逸脱している１つ以上のアノテーション値を特定することと、
（ｉｖ）前記アノテーション値をユーザに対して可視化する一方で、前記統計的に逸脱しているアノテーション値を省略して、可視化するのを防止することと、
を実行することによって、実施し、
前記アノテーション可視化操作のうちの１つ以上にわたって、アノテーション値の省略率を評価し、
前記省略率が既定の閾値を超えることを検出したことに応じて、修正措置を取る
ように構成されているプロセッサと、
を備えるシステム。
（２） 前記プロセッサが、
省略された前記アノテーション値のうちの１つ以上を再適用し、
前記統計的に逸脱しているアノテーション値を再び特定する、
ことによって、前記修正措置を取るように構成されている、実施態様１に記載のシステム。
（３） 前記プロセッサが、前記心内信号から前記アノテーション値のうちの１つ以上を再抽出することによって、前記修正措置を行うように構成されている、実施態様１に記載のシステム。
（４） 前記プロセッサが、前記逸脱の前記尺度を、前記アノテーション値の標準スコアを用いて定義するように構成されている、実施態様１に記載のシステム。
（５） 前記プロセッサが、前記逸脱の前記尺度を、前記アノテーション値の１つ以上の百分位数を用いて定義するように構成されている、実施態様１に記載のシステム。
【0101】
（６） 所与のアノテーション可視化操作において、前記プロセッサが、前記心臓内で互いから既定の距離以内に位置する、空間的に関連する電極の選択されたサブセットによって取得された心内信号に対する前記アノテーション値の逸脱を、計算するように構成されている、実施態様１に記載のシステム。
（７） 所与のアノテーション可視化操作に対する前記アノテーション値の逸脱を計算する際に、前記プロセッサが、既定の時間期間内に発生する複数の時間的に関連する心周期にわたって、前記心内信号の平均をとるように構成されている、実施態様１に記載のシステム。
（８） 所与のアノテーション可視化操作において、前記プロセッサが、前記群内の所与の電極によって取得された所与の心内信号内の前記アノテーション値のうちの１つ以上を修正して、前記群内のその他の電極に対する前記所与の電極の変位を補償するように構成されている、実施態様１に記載のシステム。
（９） 前記アノテーション値が、局所的興奮時間（ＬＡＴ）を含む、実施態様１に記載のシステム。
（１０） 前記プロセッサが、前記統計的に逸脱しているアノテーション値以外の前記アノテーション値を、前記心臓のモデル上に重ね合わせることで、前記アノテーション値を可視化するように構成されている、実施態様１に記載のシステム。
（７） 所与のアノテーション可視化操作に対する前記アノテーション値の逸脱を計算する際に、前記プロセッサが、既定の時間期間内に発生する複数の時間的に関連する心周期にわたって、前記心内信号の平均をとるように構成されている、実施態様１に記載のシステム。
（８） 所与のアノテーション可視化操作において、前記プロセッサが、前記群内の所与の電極によって取得された所与の心内信号内の前記アノテーション値のうちの１つ以上を修正して、前記群内のその他の電極に対する前記所与の電極の変位を補償するように構成されている、実施態様１に記載のシステム。
（９） 前記アノテーション値が、局所的興奮時間（ＬＡＴ）を含む、実施態様１に記載のシステム。
（１０） 前記プロセッサが、前記統計的に逸脱しているアノテーション値以外の前記アノテーション値を、前記心臓のモデル上に重ね合わせることで、前記アノテーション値を可視化するように構成されている、実施態様１に記載のシステム。
【0102】
（１１） 患者の心臓内の心臓内プローブの複数の電極によって取得された、複数の心内信号を受信することと、
後続の時点で、一連のアノテーション可視化操作を、各操作において、
（ｉ）前記心内信号から複数のアノテーション値を抽出することと、
（ｉｉ）前記心内信号の群を選択することと、
（ｉｉｉ）前記群内で、逸脱の既定の尺度よりも統計的に逸脱している１つ以上のアノテーション値を特定することと、
（ｉｖ）前記アノテーション値をユーザに対して可視化する一方で、前記統計的に逸脱しているアノテーション値を省略して、可視化するのを防止することと、
を実行することによって、実施することと、
前記アノテーション可視化操作のうちの１つ以上にわたって、アノテーション値の省略率を評価することと、
前記省略率が既定の閾値を超えることを検出したことに応じて、修正措置を取ることと、
を含む方法。
（１２） 前記修正措置を取ることが、
省略された前記アノテーション値のうちの１つ以上を再適用することと、
前記統計的に逸脱しているアノテーション値を再び特定することと、
を含む、実施態様１１に記載の方法。
（１３） 前記修正措置を取ることが、前記心内信号から前記アノテーション値のうちの１つ以上を再抽出することを含む、実施態様１１に記載の方法。
（１４） 前記逸脱の前記尺度が、前記アノテーション値の標準スコアを用いて定義される、実施態様１１に記載の方法。
（１５） 前記逸脱の前記尺度が、前記アノテーション値の１つ以上の百分位数を用いて定義される、実施態様１１に記載の方法。
後続の時点で、一連のアノテーション可視化操作を、各操作において、
（ｉ）前記心内信号から複数のアノテーション値を抽出することと、
（ｉｉ）前記心内信号の群を選択することと、
（ｉｉｉ）前記群内で、逸脱の既定の尺度よりも統計的に逸脱している１つ以上のアノテーション値を特定することと、
（ｉｖ）前記アノテーション値をユーザに対して可視化する一方で、前記統計的に逸脱しているアノテーション値を省略して、可視化するのを防止することと、
を実行することによって、実施することと、
前記アノテーション可視化操作のうちの１つ以上にわたって、アノテーション値の省略率を評価することと、
前記省略率が既定の閾値を超えることを検出したことに応じて、修正措置を取ることと、
を含む方法。
（１２） 前記修正措置を取ることが、
省略された前記アノテーション値のうちの１つ以上を再適用することと、
前記統計的に逸脱しているアノテーション値を再び特定することと、
を含む、実施態様１１に記載の方法。
（１３） 前記修正措置を取ることが、前記心内信号から前記アノテーション値のうちの１つ以上を再抽出することを含む、実施態様１１に記載の方法。
（１４） 前記逸脱の前記尺度が、前記アノテーション値の標準スコアを用いて定義される、実施態様１１に記載の方法。
（１５） 前記逸脱の前記尺度が、前記アノテーション値の１つ以上の百分位数を用いて定義される、実施態様１１に記載の方法。
【0103】
（１６） 所与のアノテーション可視化操作において、前記心臓内で互いから既定の距離以内に位置する、空間的に関連する電極の選択されたサブセットによって取得された心内信号に対する前記アノテーション値の逸脱を計算することを含む、実施態様１１に記載の方法。
（１７） 所与のアノテーション可視化操作に対する前記アノテーション値の逸脱を計算する際に、既定の時間期間内に発生する複数の時間的に関連する心周期にわたって、前記心内信号の平均をとることを含む、実施態様１１に記載の方法。
（１８） 所与のアノテーション可視化操作において、前記群内の所与の電極によって取得された所与の心内信号内の前記アノテーション値のうちの１つ以上を修正して、前記群内のその他の電極に対する前記所与の電極の変位を補償することを含む、実施態様１１に記載の方法。
（１９） 前記アノテーション値が、局所的興奮時間（ＬＡＴ）を含む、実施態様１１に記載の方法。
（２０） 前記アノテーション値を可視化することが、前記統計的に逸脱しているアノテーション値以外の前記アノテーション値を、前記心臓のモデル上に重ね合わせることを含む、実施態様１１に記載の方法。
（１７） 所与のアノテーション可視化操作に対する前記アノテーション値の逸脱を計算する際に、既定の時間期間内に発生する複数の時間的に関連する心周期にわたって、前記心内信号の平均をとることを含む、実施態様１１に記載の方法。
（１８） 所与のアノテーション可視化操作において、前記群内の所与の電極によって取得された所与の心内信号内の前記アノテーション値のうちの１つ以上を修正して、前記群内のその他の電極に対する前記所与の電極の変位を補償することを含む、実施態様１１に記載の方法。
（１９） 前記アノテーション値が、局所的興奮時間（ＬＡＴ）を含む、実施態様１１に記載の方法。
（２０） 前記アノテーション値を可視化することが、前記統計的に逸脱しているアノテーション値以外の前記アノテーション値を、前記心臓のモデル上に重ね合わせることを含む、実施態様１１に記載の方法。
【図1】
【図2】
【図3A】
【図3B】
【図4】
【図5A】
【図5B】
【図5C】
【図6A】
【図6B】
【図6C】
【図6D】
【図6E】
【図7】