(57)【要約】      （修正有）【課題】ヒトの心房についての正確な位相マップを生成する。【解決手段】心電図マッピングシステム１０は、プロセッサーによって実行可能な命令として非一過性の機械読取り可能な媒体に格納可能な信号処理方法１２を含む。信号処理方法１２は、非侵襲的な電気データに基づいて位相データ１４を計算するようにプログラムされる。非侵襲的な電気データ１６は、事前に非侵襲的に取得したデータ、あるいは、治療の手順（心臓アブレーションを含んでもよい）のあいだなど手技のあいだに非侵襲的な電極（たとえば、体表面電極）によって取得可能なデータに対応していてもよい。生成された位相データは、メモリーに格納可能であり、三次元の心臓包絡面における複数の特定された位置それぞれの電気活動の位相に対応する。このような心臓包絡面は、心外膜表面および心内膜表面あるいは、患者の心臓に対して三次元で既知の相対位置を有する任意の３Ｄ膜に対応してもよい。【選択図】図１
2019141619.tif 000002
【特許請求の範囲】
【請求項１】
機械読取り可能な命令およびデータを格納するためのメモリーと、
前記メモリーにアクセスし、前記命令を実行するためのプロセッサーと、を含むシステムであって、前記機械読み取り可能な命令は、
少なくとも１つの時間間隔のあいだ患者から得られる非侵襲的な電気データに対応する、処理された電気データを、あらかじめ定められた心臓包絡面（cardiac envelope）上の対応する再構成された電気信号に変換し、
前記再構成された電気信号に基づいて位相データを計算すること
を含む方法を行うためのものであり、
前記再構成された電気信号は、時空の一貫性がある、システム。
【請求項２】
前記方法は、前記計算された位相データに基づいて、位相マップを作成することをさらに含む、請求項１に記載のシステム。
【請求項３】
前記方法は、
前記位相データに基づいて、らせん波特性を計算し、
前記計算された位相データに基づいて、前記少なくとも１つの時間間隔にわたって、らせん波マップを作成すること
をさらに含む、請求項２に記載のシステム。
【請求項４】
患者の胴のあらかじめ定められた部分をカバーし、かつ、前記患者の胴の外面から電気活動を非侵襲的に検出するよう構成された電極群をさらに含み、
前記検出される電気活動は、前記非侵襲的な電気データに対応する、請求項２に記載のシステム。
【請求項５】
前記計算された位相データに基づいて特定される空間位置で心臓組織に対して治療を行うよう構成された治療デバイスをさらに含む、請求項４に記載のシステム。
【請求項６】
前記計算された位相データに基づいて特定されるらせん波の旋回中心に対応する空間位置で心臓組織を焼灼するよう構成された、アブレーションカテーテルをさらに含む、請求項４に記載のシステム。
【請求項７】
前記方法は、少なくとも１つの時間間隔にわたって、らせん波の旋回中心についての空間位置の時間周波数を計算することをさらに含む、請求項２に記載のシステム。
【請求項８】
前記方法は、
前記心臓包絡面の位相勾配の積分の指標を推定し、
前記心臓包絡面に対応する前記患者の心臓の複数の腔について時空の一貫性のある位相情報を提示するために、前記推定された指標に基づいて積分位相勾配マップを作成すること
をさらに含む、請求項２に記載のシステム。
【請求項９】
前記方法は、
前記心臓包絡面上の少なくとも１つの位相特異点の位置を決定し、前記少なくとも１つの位相特異点の前記位置を前記心臓包絡面についての前記位相マップに表示することをさらに含む、請求項２に記載のシステム。
【請求項１０】
前記方法は、前記心臓包絡面上の少なくとも１つのフォーカス源の位置を決定し、前記少なくとも１つのフォーカス源の前記位置を前記心臓包絡面についての前記位相マップに表示することをさらに含む、請求項２に記載のシステム。
【請求項１１】
前記方法は、前記変換の前に、少なくとも１つの時間間隔のあいだ前記患者から得られる非侵襲的な電気データを前処理して、不整脈に寄与すると決定された信号の特徴を抽出し、前記抽出された信号の特徴に対応する前記処理された電気データを与えることをさらに含み、
前記再構成された電気信号は、前記処理された電気データから誘導される、請求項１に記載のシステム。
【請求項１２】
前記方法は、前記変換の後で、前記再構成された電気信号を処理して処理後の再構成された電気信号を与えることをさらに含み、
前記計算された位相データは、前記処理後の再構成された電気信号に基づいて計算される、請求項１に記載のシステム。
【請求項１３】
前記方法は、前記変換の前に、少なくとも１つの時間間隔のあいだ患者から得られる非侵襲的な電気データを前処理して、あらかじめ規定された不整脈に寄与しないと決定された信号の特徴を除去するとともに、前記処理された電気データを与えることをさらに含み、
前記再構成された電気信号は、前記処理された電気データから誘導される、請求項１に記載のシステム。
【請求項１４】
前記前処理は、心室の電気活動に由来する信号の特徴を前記非侵襲的な電気データから除去し、心房タイプの不整脈に対する特異性を高めることをさらに含む、請求項１３に記載のシステム。
【請求項１５】
前記信号の特徴を除去することは、前記非侵襲的な電気データによって表される信号から、ＱＲＳ波を相殺すること、またはＴ波を相殺することのうちの少なくとも１つをさらに含む、請求項１４に記載のシステム。
【請求項１６】
前記方法は、ユーザー入力に応答して、あらかじめ規定された不整脈として心房性不整脈と心室性不整脈のいずれが選択されたかに応じて行われる前記前処理を選択的に制御することをさらに含み、
前記前処理は、前記ユーザー入力が前記心房性不整脈を選択している場合、心室の電気活動を表す信号の特徴を除去するようにプログラムされる、請求項１３に記載のシステム。
【請求項１７】
前記再構成された電気信号は、前記心臓包絡面に分散した１０００箇所を超える位置についての電気情報を与える、請求項１に記載のシステム。
【請求項１８】
前記心臓包絡面は、心外膜表面であり、前記あらかじめ規定された不整脈は、心房細動、心房頻拍、心室頻拍、および心室細動のうちの少なくとも１つである、請求項１に記載のシステム。
【請求項１９】
前記変換することが、特定の患者について取得される実際の幾何学形状データと一般的な解剖学的モデルとのうちの少なくとも１つである幾何学形状に基づいて、逆解を計算することによって行われる、請求項１に記載のシステム。
【請求項２０】
前記方法は、前記計算された位相データに基づいてサイクル長の指標を計算するようさらにプログラムされる、請求項１に記載のシステム。
【請求項２１】
プロセッサーによって実行可能な命令を有する非一過性のコンピューター読取り可能な媒体であって、前記命令は、
幾何学形状データおよび非侵襲的に測定される体表面電気信号に基づいて、あらかじめ定められた心臓包絡面上またはその付近にある多数のポイント各々について再構成されたエレクトログラム信号を生成するためのエレクトログラム再構成方法と、
前記再構成されたエレクトログラム信号に基づいて前記多数のポイントについて位相信号を計算するための位相計算部と、
前記計算された位相信号に基づいて出力を生成するための可視化エンジンと、を含む、媒体。
【請求項２２】
前記命令は、前記非侵襲的に測定される体表面電気信号を処理し、あらかじめ規定された不整脈に寄与しないと決定された信号の特徴を除去するとともに、対応する処理後の電気信号を与えるようにプログラムされる、前処理方法をさらに含み、
前記エレクトログラム再構成方法は、前記対応する処理後の電気信号から前記再構成された電気信号を生成する、請求項２１に記載の媒体。
【請求項２３】
前記前処理方法は、心室の電気活動による前記非侵襲的に測定される体表面電気信号から信号の特徴を除去し、前記対応する処理後の電気信号の心房タイプの不整脈に対する特異性を高めるようさらにプログラムされる、請求項２２に記載の媒体。
【請求項２４】
前記前処理方法は、前記非侵襲的に測定される体表面電気信号からの、ＱＲＳ波の相殺またはＴ波の相殺のうちの少なくとも１つを行い、前記非侵襲的に測定される体表面電気信号から前記信号の特徴を除去するようさらにプログラムされる、請求項２３に記載の媒体。
【請求項２５】
前記命令は、ユーザー入力に応答しての評価のために、あらかじめ定められたタイプの不整脈を指定するようにプログラムされるグラフィカルユーザーインターフェースをさらに含み、
前記前処理方法は、前記ユーザー入力に応答して前記前処理を選択的に制御するようにプログラムされる、請求項２２に記載の媒体。
【請求項２６】
前記前処理方法は、前記ユーザー入力が心房タイプの不整脈の評価を選択した場合、心室の電気活動に関連する信号の特徴を前記非侵襲的に測定される体表面電気信号から除去するようにプログラムされる、請求項２５に記載の媒体。
【請求項２７】
前記命令は、前記計算された位相信号に基づいて前記心臓包絡面上の少なくとも１つの位相特異点の位置を決定するようにプログラムされる、位相特異点計算をさらに含み、
前記少なくとも１つの位相特異点の前記位置は、前記出力の中で特定される、請求項２１に記載の媒体。
【請求項２８】
前記命令は、前記計算された位相信号に基づいて前記心臓包絡面上の少なくとも１つのフォーカス源の位置を決定するようにプログラムされる、フォーカス源計算部をさらに含み、
前記少なくとも１つのフォーカス源は、前記出力において特定される、請求項２１に記載の媒体。
【請求項２９】
前記可視化エンジンは、少なくとも１つの時空の一貫性のあるマップを作成し、複数の時間間隔にわたって患者の心臓の複数の腔について不整脈の１つ以上の機序を特徴付けするようにさらにプログラムされる、請求項２１に記載の媒体。
【請求項３０】
前記出力は、前記計算された位相信号に基づいて患者の心臓の複数の腔について時空の一貫性のある情報を提示するための、位相マップ、積分位相勾配マップ、らせん波ロードマップ、サイクル長マップのうちの少なくとも１つを含む、請求項２９に記載の媒体。
【請求項３１】
前記命令は、前記計算された位相信号に基づいてアブレーションの標的としてらせん波の旋回中心の位置を特定するようにプログラムされる、らせん波特定機能をさらに含む、請求項２１に記載の媒体。
【請求項３２】
前記エレクトログラム再構成方法によって利用される前記幾何学形状データは、特定の患者について取得される実際の幾何学形状データと一般的な解剖学的モデルとのうちの少なくとも１つである、請求項２１に記載の媒体。
【請求項３３】
前記命令は、前記計算された位相信号に基づいてサイクル長の指標を計算するようにプログラムされる、サイクル長計算機能をさらに含む、請求項２１に記載の媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、細動または頻拍などの不整脈用の生理学的マッピングに関する。
【背景技術】
【０００２】
心電図マッピング（ＥＣＭ）は、検知された電気信号から心臓の電気情報を決定し、表示するのに用いられている技術である。心臓の電気活動のマッピングは、心房細動および心室細動を含む細動などの特定タイプの不整脈が存在することから、さらに複雑化している。心臓の信号は、組織のいくつかの連続した脱分極を含む。心臓の脱分極または興奮についての解析内容を説明する、今日実用化されているほとんどの方法は、信号の周期性またはサイクル長に関する推測的な知識（たとえば、３Ｄ興奮マッピング用）に頼っている。特定の周期性またはサイクル長のない心臓の興奮では、位相マッピングなどの解析技術が開発されている。しかしながら、既存の位相マッピング手法は、心臓に挿入される１本以上のバスケットカテーテルを含んだり、心筋に染料を注射して光学イメージングを行うなど、侵襲的である。光学イメージング染料は、ヒトでは有毒であり、バスケットカテーテルで得られる分解能は、ヒトの心房についての正確な位相マップを生成するには、十分ではない。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
本開示は、細動または頻拍などの不整脈用の生理学的マッピングに関する。
【０００４】
一例として、システムは、機械読取り可能な命令およびデータを格納するメモリーと、メモリーにアクセスし、方法を行うための命令を実行するプロセッサーと、を含んでもよい。この方法は、少なくとも１つの時間間隔のあいだに患者から得られる非侵襲的な電気データに対応する、処理された電気データを、あらかじめ定められた心臓包絡面（cardiac envelope）上の対応の再構成された電気信号に変換することを含んでもよい。再構成された電気信号は、時空の一貫性のあるものであってもよい。また、この方法は、再構成された電気信号をもとに、位相データを計算することを含んでもよい。
【０００５】
別の例として、非一過性のコンピューター読取り可能な媒体は、プロセッサーによって実行可能な命令を有するものであってもよい。これらの命令は、幾何学形状データおよび非侵襲的に測定された体表面の電気信号に基づいて、あらかじめ定められた心臓包絡面上またはその付近にある多数のポイント各々について再構成されたエレクトログラム信号を生成するための、エレクトログラム再構成方法を含んでもよい。これらの命令は、再構成されたエレクトログラム信号に基づいて多数のポイントについて位相信号を計算するための位相計算部と、計算された位相信号に基づいて出力を生成するための可視化エンジンとを含んでもよい。
【０００６】
＝＝関連出願へのクロスリファレンス＝＝
本出願は、発明の名称「PHYSIOLOGICAL MAPPING FOR ARRHYTHMIA（不整脈用の生理学的マッピング）」２０１２年９月２１日にファイルされた米国仮特許出願第６１／７０４，２２７号（その内容全体を本明細書に援用する）の優先権の利益を主張する。
【図面の簡単な説明】
【０００７】
【図１】図１は、不整脈についての位相マッピングを行うためのシステムの一例を示す。
【図２】図２は、心房細動について実装可能な前逆解（pre-inverse solution）処理の一例を示す。
【図３】図３は、心室細動について実装可能な前逆解処理の一例を示す。
【図４】図４は、後逆解（post-inverse solution）処理の一例を示す。
【図５】図５は、ベースライン除去プロセスの一部として利用可能な方法論の一例を示す。
【図６】図６は、心房信号および対応のベースライン信号の一例を示す。
【図７】図７は、フィルタリング前後の心室信号のプロットを示す。
【図８】図８は、可視画像の生成に用いられる位相データを生成するために実装可能な位相計算部ツールの一例を示す。
【図９】図９は、心房信号およびその対応の変換結果の一例を示す。
【図１０】図１０は、複素平面にプロットされた心房信号の位相の一例を示す。
【図１１】図１１は、複素平面にプロットされた心室信号の一例を示す。
【図１２】図１２は、診断治療システムの一例を示す。
【図１３】図１３は、位相マップの一例を示す。
【図１４】図１４は、らせん波の旋回中心の特定を明示する位相マップの一例を示す。
【図１５】図１５は、サイクル長マップの一例を示す。
【図１６】図１６は、特異点を明示する心電図マップの一例を示す。
【図１７】図１７は、本明細書に開示のシステムおよび方法に基づいて生成可能な心臓の動画位相マップの一例を示す。
【図１８】図１８は、本明細書に開示のシステムおよび方法によって計算された位相マップデータに基づいてらせん波の旋回中心の位置と軌道を決定するのに利用可能なアルゴリズムの異なる部分を示す。
【図１９】図１９は、本明細書に開示のシステムおよび方法によって計算された位相マップデータに基づいてらせん波の旋回中心の位置と軌道を決定するのに利用可能なアルゴリズムの異なる部分を示す。
【図２０】図２０は、本明細書に開示のシステムおよび方法によって計算された位相マップデータに基づいてらせん波の旋回中心の位置と軌道を決定するのに利用可能なアルゴリズムの異なる部分を示す。
【図２１】図２１は、断片化およびそれぞれのらせん波間の相互作用を含み得る時間の関数としてのらせん波の動きを明示するグラフの一例を示す。
【図２２】図２２は、異なるらせん波特性を明示するのに位相マップに適用可能なグラフィック要素の一例を示す。
【図２３】図２３は、心房粗動の位相マップの動画表現を示す。
【図２４】図２４は、心臓包絡面の一部についての近傍図の一例を示す。
【図２５】図２５は、システムおよび方法を行うことのできる例示的なコンピューター環境を示す。
【発明を実施するための形態】
【０００８】
本開示は、細動などの不整脈についての非侵襲的な電気生理学的マッピングに関する。このマッピング技術を診断および／または治療のワークフローの一部として用いて、体表面で非侵襲的に測定される体表面での電気活動に基づく細動機序の特定と位置特定を容易にすることが可能である。非侵襲的な体表面での測定を利用して、高分解能の電気信号を、心臓の表面などの心臓包絡面上で再構成することが可能である。本明細書に開示の手法は、電気データに対する逆解を計算する前、途中、および／または後の処理を含む、信号処理の技術を採用できる。このような信号処理は、高分解能の出力を可能にし、光学的マッピングでのように染料を注射したり侵襲的な測定カテーテルを使用したりすることなく、上記の出力から、治療標的を特定できる。
【０００９】
また、本明細書に開示の手法は、リアルタイムでも（たとえば、電気生理学の研究中または治療行為のあいだにガイダンスを与えるために）利用できるし、あるいは、特定の患者で事前に取得して保存された電気データとの関連で実装することも可能である。この信号処理では、生理学的に関連する情報を非侵襲的な電気データから抽出し、１つ以上の時間間隔にわたる位相を計算することが可能である。こうして計算された位相を利用して、患者の心臓における細動の機序（たとえば、局所的ならせん波、ミクロおよびマクロリエントリー回路、誘発活動、フォーカス点などを含むが、これらに限定されるものではない）を特定し、特徴付けることができる。また、得られる位相データを利用して、グラフィックな可視画像を生成し、１つ以上のマップから時空の一貫性のある情報を提示してもよい。これらのマッピング出力を、さらに、らせん波およびこれに付随する動きや臨床標的として時空の一貫性のある不整脈の永続の他の特徴を示すような、動的な動画映像などの３Ｄ解剖学的マップとしてグラフィック表示することもできる。
【００１０】
図１は、心電図マッピングシステム１０の一例を示す。システム１０は、プロセッサーによって実行可能な命令として非一過性の機械読取り可能な媒体に格納可能な信号処理方法１２を含む。信号処理方法１２は、非侵襲的な電気データ１６に基づいて位相データ１４を計算するようにプログラムされる。非侵襲的な電気データ１６は、事前に非侵襲的に取得されて以後の処理用にメモリーに格納された保存電気データに対応していてもよい。あるいは、非侵襲的な電気データ１６は、電気生理学的研究のあいだや、治療の手順（心臓アブレーションを含んでもよい）のあいだなど、手技のあいだに非侵襲的な電極（たとえば、体表面電極）によって取得可能な、リアルタイムのデータフローに対応していてもよい。生成された位相データ１４は、メモリーに格納可能であり、三次元の心臓包絡面における複数の特定された位置それぞれの電気活動の位相に対応する。このような心臓包絡面は、心外膜表面および心内膜表面あるいは、患者の心臓に対して三次元で既知の相対位置を有する任意の３Ｄ膜に対応してもよい。
【００１１】
システム１０は、グラフィック出力１８として明示される対応の可視画像を生成するようにプログラムされる可視化エンジン１９を含んでもよい。出力１８は、たとえば、心臓包絡面上に分散しているポイントに対して計算される位相データ１４に基づいて生成される位相マップを含んでもよい。その代わりにまたはそれに加えて、位相データ１４は、非侵襲的な電気データ１６によって表される信号に基づいて計算される位相特異点（たとえば、らせん波の旋回中心）または他の位相特性に対応してもよい。たとえば、グラフィック出力１８は、らせん波の旋回中心の時空の一貫性を明示する位相勾配の積分（たとえば、積分位相勾配マップ）を示すマップを含んでもよい。既存の侵襲的な技術と異なり、積分位相マップは、（たとえば、現在の両心房または両心室の同時位相特性に関して）心臓の複数の腔について時空の一貫性のあるものとなり得る。
【００１２】
別の例として、可視化エンジン１９は、ムービーまたは映画撮影のフォーマットで位相マップを生成し、１つ以上の時間間隔の時間にわたって、らせん波の動と相互作用を動的に明示できる。さらに他の例では、可視化エンジン１９は、電位マップ、興奮マップ、優位周波数マップなどを含んでもよい他のデータに基づいてグラフィック出力およびマップを生成することも可能である。電気ポテンシャルの活動に関連したこのような追加のマップを、位相データと組み合わせ、１つのディスプレイにしてもよい。その代わりにまたはそれに加えて、可視化エンジン１９は、電位ベースの心電図マッピング（たとえば、興奮マップ、優位周波数マップなど）に隣接するウィンドウまたはディスプレイに、位相マップまたは関連する位相特性を表示する複数のディスプレイを同時に生成可能である。
【００１３】
図１の例では、信号処理１２は、位相データ１４を生成すべく非侵襲的な電気データ１６を処理するために実行可能な複数の工程（たとえば、プログラムモジュール）を含んでもよい。たとえば、以後の処理用に、非侵襲的に取得した電気データ１６に前逆解処理（本明細書では、前処理とも呼ぶ）２０を適用し、非侵襲的な電気データ１６に含まれる各信号の関連のある部分を抽出する（または関連のない部分を除去する）などしてもよい。このように、前処理２０は、逆解２２を駆動するのに用いられる前処理された電気データを与えることができる。
【００１４】
逆解法２２は、逆解に供給される前処理された電気データおよび幾何学形状データ（図示せず）に基づいて、心臓包絡面上の再構成されたエレクトログラムを計算可能である。幾何学形状データは、特定の患者について取得される実際の幾何学形状データ、一般的な解剖学的モデル、またはこれらの組み合わせに対応していてもよい。前逆解処理２０は、前処理された電気データを与えて、逆解２２に供給される入力データの全体としての感度を落としつつ位相データ１４の計算対象となる特定タイプの不整脈（たとえば、心房細動（ＡＦ）、心房頻拍（ＡＴ）、心室細動（ＶＦ）、心室の（ＶＴ））の各々の特異性を高めてもよい。前逆解処理２０による信号の特徴の除去に加えて、このような前処理では、ノッチフィルタリング、低域通過フィルタリング、不良チャネルの除外などのフィルタリングによって、各信号チャネルからノイズ（たとえば、望ましくない振動に対応するノイズ）も除去できる。
【００１５】
一例として、前逆解処理２０は、ＡＦの位相データを計算するときに、心室機能（たとえば、一方または両方の心室の脱分極、心室の再分極、またはその両方）による信号の寄与を除去するようプログラム可能である。たとえば、逆解２２を計算する前に、信号におけるＱＲＳ波の相殺および／またはＴ波の相殺（たとえば、ＱＴ間隔またはその選択された一部）によって、心室の電気活動による信号の特徴を除去するための特徴抽出を行い、心臓包絡面上に心電図を再構成することが可能である。このような相殺では、波全体または選択された他の着目信号成分の特異性を高めるのに十分な任意の一部を除去できる。信号前処理２０だけでなく、他の方法を利用して、ＡＦ位相診断のための値を計算しながら心室信号による作用を緩和してもよいことは、理解されたい。このような方法は、化学物質（たとえば、薬剤）の使用、電気刺激の使用、これらの組み合わせを含んでもよい。このタイプの前処理２０を、さらにプログラムし、どのタイプの不整脈を解析するかに応じて、非侵襲的に取得された電気信号から不要な信号成分を除去してもよい。
【００１６】
本明細書で開示するように、逆解２２は、前処理されたデータに基づいて再構成されたエレクトログラムを与えるようにプログラムされる。システム１０で利用可能な逆アルゴリズムの例は、米国特許第７，９８３，７４３号および同第６，７７２，００４号（各々本明細書に援用する）に開示されている。逆解２２は、心臓包絡面上の複数の位置（たとえば、約２０００箇所以上など、１０００箇所を上回る位置）で電気活動を再構成できる。結果として、処理によって与えられる位相データ１４を、（たとえば、バスケットカテーテルのコンタクト電極を用いて）侵襲的な方法で実現できるよりも有意に大きい、対応の高い分解能のものとすることができる。また、位相データ１４は、有毒な化学染料を使用する光学マッピング技術を用いて生成されるものよりかなり安全であり得る。非侵襲的な電気データ１６によって表される入力信号から不要な情報を削除する前処理２０を行うことにより、逆解２２を残りの信号情報に適用し、選択されたタイプの位相解析の特異性の高さと特別に適合された、心臓包絡面で再構成された対応の電気信号を与えることが可能である。
【００１７】
追加の後逆解処理２４は、再構成された電気データに適用可能である。後逆解処理２４は、１つ以上のフィルタリング方法、たとえば、ベースライン除去、帯域通過フィルタリング、および低域通過フィルタリング（たとえば、図４を参照のこと）、を含んでもよい。
【００１８】
位相計算２６は、フィルタリング後の再構成されたエレクトログラムデータから位相データ１４を計算するようプログラム可能である。たとえば、位相計算２６は、フィルタリング後のエレクトログラムデータに基づいて、信号の分解を行い、心臓包絡面の表面にある複数のポイント各々の位相信号を計算することができる。位相計算は、本明細書に開示するような様々な解に応じて実装可能である（たとえば、図８に関して説明する位相計算を参照のこと）。
【００１９】
本明細書で開示するように、可視化エンジン１９はさらに、さらに心臓包絡面の位置ごとに計算された位相信号に基づいて、位相マップおよび他の可視画像を計算可能である。位相信号計算のもとになる非侵襲的な入力電気的データ１６は、患者の心臓またはユーザーが選択した他の着目領域の表面全体にわたって時空の一貫性のある位相情報を提示できるため、位相マップを心臓包絡面（たとえば、心外膜表面全体）の概観的な真の位相マップとみなすことができる。既存の侵襲的な位相マッピング技術とは異なり、概観的な可視画像は心臓全体についての時空の一貫性のある情報を与えられるため、らせん波の空間移動と心臓での相互作用（複数の腔にわたる相互作用を含む）を特定することが可能である。たとえば、可視化エンジン１９は、位相データ１４からグラフィック出力１８を生成し、複数の心腔にわたる細動機序を特徴付けるための両心室のマップおよび／または両心房のマップを含めることができる。
【００２０】
グラフィック出力１８は、特定時刻における位相情報と対応してもよいし、時間の関数として心臓包絡面での位相変化を示すために動画化してもよい。たとえば、非侵襲的な電気データを、どのように、いつ取得したかに応じて、時間はさらに連続した時間を含んでもよいし、時間的に間隔のあいた複数の時間間隔（たとえば、数分、数時間または数日あけた別の間隔）を含んでもよい。測定間にこのような遅延があるにもかかわらず、生成される位相データ１４は依然として、各々の測定の組について空間的に一貫性のあるものであり得る。このような手法は、既存の侵襲的技術では、患者のリスクを劇的に増大させることなく実現するのは無理であろう。変化は、１つ以上の時間間隔から得られる心臓の電気活動に対応可能である。特定の時間間隔が、位相データの計算対象となった位置の各々について１つ以上の位相サイクルを含んでもよい。
【００２１】
図２は、心房細動（ＡＦ）５０の解析を容易にするよう構成可能な前逆解処理（たとえば、図１の前処理２０に対応）５０の一例を示す。方法５０は、患者の胴の表面に分散可能なセンサーのベストまたは他の電気的アレイから、非侵襲的な電気データを受信する。方法５０は、ラインフィルター周波数および高調波と、ユーザー入力に応答してユーザーが指定可能な他の周波数範囲とを除去するよう構成可能な、ノッチフィルター５２を含む。このため、ノッチフィルター５２は、プログラム可能であってもよい。いくつかの例では、ノッチフィルター５２は、高速フーリエ変換（ＦＦＴ）法を用いて実装可能である。フィルタリング後のデータは、以後の処理でどの入力チャネルを用いるのか選択するのに利用可能なチャネル選択ブロック５４に提供できる。
【００２２】
チャネル選択５４は、たとえば、不良チャネルの検出と除去を含んでもよい。不良チャネル検出では、逆解（たとえば、図１の逆解２２）の計算に有害であると決定されるチャネルを特定できる。たとえば、不良チャネルは、切断されているチャネル、欠けたチャネルのみならず、歪んだデータを持つチャンネルであってもよい。チャネル選択５４は、たとえば、信号の目視検査により、および／または自動不良チャネル検出アルゴリズムによって実装可能である。その代わりにまたはそれに加えて、チャネル選択５４は、特定タイプの解析（たとえば、ＡＦまたはＡＴ）および心臓または選択された着目領域についてなされる評価にプラスに寄与するよう、検知対象となる信号があらかじめ定められている患者の胴の１つ以上の領域に対応する一組の入力チャネルを選択するよう実装可能である。
【００２３】
選択チャネルに対するフィルタリング後のデータを心室興奮除去ブロック５６に与えて、心室の興奮に対応する信号の特徴を除去することが可能である。一例として、心室興奮除去５６は、ＱＲＳ波サブトラクション法５８を含んでもよい。たとえば、ＱＲＳサブトラクションは、仮想リードの主成分分析によって実装可能であり、この場合、ＱＲＳ波のテンプレート領域の位置を特定し、テンプレート領域を平均し、対応の着目領域テンプレートを作成し、テンプレートを動的に調整し、検出されたテンプレート領域で、フィルタリング後の入力信号から調整したテンプレートを引く。
【００２４】
その代わりにまたはそれに加えて、心室興奮除去５６は、選択されたチャネルによって与えられるフィルタリング後の信号に対して実行されるＴ波サブトラクション６０を含んでもよい。Ｔサブトラクションは、ＱＲＳ波サブトラクション５８と類似であってもよいが、Ｔ波のサブトラクションに対応するよう構成された別のテンプレートを使用できる。あるいは、ＱＲＳ波とＴ波の両方に対応する単一のテンプレートを使用して、選択されたチャネルに対して与えられるフィルタリング後の入力信号からＱＲＳ波とＴ波の両方を除去してもよい。
【００２５】
このように、心室興奮除去５６は、処理された信号を与えることができる。この処理された信号は、対応の心房解析に関連して特定された信号の一部を含む。信号からそのような特徴を差し引くことにより、実装対象となる位相解析の特異性を高めるために、多少のレベルの感度を犠牲にしてもよい。心室興奮除去５６は、低域通過フィルター６２に対応する信号を与えることができる。この低域通過フィルター６２は、約４０ヘルツなどのあらかじめ定められたカットオフ周波数未満の帯域を通過させるよう構成可能である。一例として、低域通過フィルター６２はSavitzy Golayフィルターとして実装可能である。他のタイプのフィルターを利用してもよい。
【００２６】
図３は、心室性不整脈（たとえば、ＶＦまたはＶＴ）の解析を容易にするよう実装可能な前逆解処理７０の一例を示す。これはたとえば、図１の前逆解処理２０に対応していてもよい。２つ以上のタイプの異なる前逆解処理があれば、本明細書に開示の手法では、実施する対象の診断および評価（ユーザー入力に応答して選択可能である）に従って、どのタイプの前処理を行うかを選択的に制御することが可能である。
【００２７】
図３の例では、前逆解処理７０は、ＦＦＴ法を用いて実装できるようなノッチフィルター７２を含む。ノッチフィルター７２は、以後の処理でどのチャネルが利用されるか、または除外されるかを選択するよう構成可能なチャネル選択部７４に、対応するフィルタリング後の信号を与える。たとえば、チャネル選択ブロック７４は、心房の活動より心室の活動に多く寄与していると決定されるような入力データ中のチャネルを保持するよう構成可能である。追加の信号サブトラクションまたはフィルタリングを利用して、非侵襲的な入力データ（たとえば、図１の入力データ１６）から心房の特徴を除去してもよい。チャネル選択７４は、低域通過フィルター７６に、選択されたチャネルのフィルタリング後の信号を与え、心室性不整脈評価用の位相診断を行うための着目信号成分を含む対応のフィルタリング後の出力を提供してもよい。
【００２８】
図４は、図１の後逆解処理２４に対応可能な後逆解法（post-inverse method）処理８０の一例を示す。後逆解法処理８０は、心臓包絡面（たとえば、患者の心臓の心外膜表面）での電気活動に対応する再構成されたエレクトログラムデータに適用できる。後逆解法処理８０は、指定の周波数で１つ以上の帯域を通過させるようにプログラムされる帯域通過フィルター８２を含んでもよい。一例として、帯域通過フィルター８２は、ＦＩＲフィルターとして実装可能である。
【００２９】
帯域通過フィルタリング後の信号を、フィルタリング後の信号からベースラインワンダーを除去するよう構成されたベースライン除去機能８４に与えることができる。たとえば、たとえば心房信号に対する信号平均値のサブトラクションによって、ベースラインワンダーを相殺できる。別の例では、ベースライン除去８４は、多項式関数でベースラインを内挿するようにプログラムされる方法によって実装可能である。一例として、ベースラインワンダーを除去するために、各信号について選択されるアンカーポイントから、ベースラインの内挿を計算することが可能である。
【００３０】
ベースラインを除去するのに利用可能な方法論の一例を、図５に示す。図５の方法論に従ってスプライン内挿を行い、特定の信号からのアンカーポイントをリンクさせた後、得られるベースラインを上記の特定の信号から差し引くことができる。その後、得られるベースラインの補正信号それぞれを低域通過フィルター８６に通し、（たとえば、図１の位相計算２６による）位相計算をはじめとする以後の計算用に、対応の処理信号を与えればよい。
【００３１】
さらに別の例として、ベースライン除去８４は、複数レベルの分解を用いるなどの方法で、ウェーブレット分解によって実装可能である。たとえば、フィルタリング後の信号を、高次のcoifletまたは他のウェーブレット（たとえば、Mexican hatウェーブレットまたはMorletウェーブレット）を用いるなどの方法で、直交ウェーブレットの基底に分解できる。
【００３２】
さらに別の例として、時間信号をｘとおくと、ウェーブレット分解は、以下を満たすスケール関数ｓiに至る。
2019141619.tif 000003
ここで、ｒは、剰余関数である。
フィルタリング８４は、診断に関連のない情報を模式化した合計からｓi関数を除去することで、行うことができる。このような目的に向けたひとつの方法は、各ｓi関数のフーリエ変換を計算し、優位周波数（たとえば、最高ピークの周波数）が生理学的周波数（たとえば、通常約４から１０Ｈｚ）の範囲外にある関数を除去することであってもよい。このようにして、ベースラインからの低周波数成分とノイズからの高周波数成分の両方を除去することが可能である。さらに、このようなウェーブレットフィルタリングは、高信頼度でＡＦサイクル長に対応する成分だけを用いて信号を再構成できるよう、すなわちデータが真の心房脱分極に対応する信号成分を含むように、起伏の大きな波形を分解し、不要な振動に寄与する成分を除去するよう動作する。
【００３３】
図６は、心房信号９０の一例を、本明細書に開示の手法のひとつに従うなどして計算された、対応するベースライン信号９２と一緒に示す。信号ごとに、計算されたベースライン信号９２をそれぞれの信号９０から引き、実質的にベースラインドリフトを含まない対応の信号を与えることが可能である。
【００３４】
図７は、図４の後逆解処理８０によって与えることが可能な信号プロット９４、９６、９８の一例を示す。図７では、生の心室細動信号９４が図示されている。また、図７は、ノイズによる振動を除去するための周波数フィルタリング後の信号９６も示す。図６も、本明細書で開示するようなベースラインワンダーを軽減するためのウェーブレットフィルタリング後の信号９８を示す。
【００３５】
図８は、処理後の再構成されたエレクトログラムデータ（たとえば、図１の方法２４または図４の方法８０によって与えられる信号に対応）から位相および位相特性を計算するための方法を行うようプログラム可能な、位相計算部ツール１００の一例を示す。位相計算部ツール１００は、ユーザー入力に応答して起動可能であり、そのために、ツール１００とともに示されるそれぞれの方法へのアクセスを提供するグラフィカルユーザーインターフェース（ＧＵＩ）などを用いてもよい。ツールによって実装されるそれぞれの方法は、逆解後のマッピング計算時に行うことができ、その出力は、対応のグラフィックな可視画像（たとえば、図１のグラフィック出力１８）を生成するのに利用可能である。
【００３６】
たとえば、ツール１００は、対応する位相データ１０４（たとえば、図１の位相データ１４に対応）を計算するようにプログラムされる位相計算方法１０２を含んでもよい。位相計算１０２は、処理後の再構成されたエレクトログラム（たとえば、図１の後逆解処理２４または図４の８０によって与えられるものなど）で計算可能であり、そこでは、再構成されたエレクトログラム信号から、着目される生理学的周波数（たとえば、約４〜１０ヘルツ）の範囲外の信号成分が除去されている。
【００３７】
一例として、位相計算１０２は、心臓包絡面全体で信号ごとに対応の位相データ１０４を得るようヒルベルト変換を用いて実装可能である。他のタイプの信号分解法を利用してもよい。このように、位相データ１０４は、再構成されたエレクトログラムが構築されている１つ以上の時間間隔を表すことが可能である。
【００３８】
さらに別の例として、−πをサイクルの任意の開始点とすると、πは、次のサイクルの開始点となる。そのあいだにある各々の時点には、［−π，π］の間の位相値が増分で割り当てられている。たとえば、得られた位相が、大きさ１の複素数の位相であると仮定する。そのように、それぞれのサイクルを複素空間における０を中心とする１つのサイクルに変換できる。
【００３９】
たとえば、心臓包絡面の特定位置における再構成された心電図をｘとする。特定のｘについて上記の要件を満たす位相信号を見つけるために、次のような位相空間「理論」を用いることができる。
実時間信号はいずれも、複素信号の実軸に正射影したものとして見られる。
【００４０】
よって、ｘを与えると、
2019141619.tif 000004
を満たすような複素信号Ｘを見いだすことが可能である。本発明者らは、心臓包絡面で再構成およびフィルタリングされた特定の電気信号についての計算された位相を表すために、
2019141619.tif 000005
とする。Ｘを与えると、
2019141619.tif 000006
は一意であるが、ｘを与えると、無限数の信号Ｘが存在する。以下は、
2019141619.tif 000007
に整合するＸに対するいくつかの考えられる解である（ここで、
2019141619.tif 000008
は、本明細書に開示の処理技術を用いて信号から抽出されるｘの、診断の関心対象の（たとえば、関連のある）部分（本明細書では
2019141619.tif 000009
で示す）に対応する。
2019141619.tif 000010
ここで、τは、推定される必要がある。
2019141619.tif 000011
このような場合、Ｘは、
2019141619.tif 000012
の解析信号と称することが可能であり、よって、以下のように表すことが可能である。
2019141619.tif 000013
2019141619.tif 000014
あるいは、ｘの複素ウェーブレットに対応しているなど、
2019141619.tif 000015
を満たすような任意のＸであってもよい。
【００４１】
位相データ１０４を計算し、一意の解を与える（１サイクルを［−π，π］の間の値に変換）位相マッピングを行うために、Ｘに求められる要件は、サイクルごとに（０，０）のまわりに１つの円を得ることである。これは、ノイズすなわち、Ｘの目立った特徴を抽出するための信号の無関係な振動を除去するための、本明細書に開示の（たとえば、図３〜図５）前処理によって達成され、それにより、位相マッピングの上述した技術の精度と再現性が高められる。さまざまな時点に、いくつかのセグメントのデータについて位相情報を計算し、時空の一貫性という意味で解析をしっかりしたものとすることが可能である。同一患者における複数のデータセグメントからの位相情報を、本明細書で開示するようなさまざまな計算技術および可視化技術を用いて組み合わせてもよい。
【００４２】
図９は、処理後の心房信号１２０の例とともに、２２で明示した当該信号のヒルベルト変換（これは、計算された位相データ１０４に対応する）を示す。このように、位相データは、電気データが得られている１つ以上の間隔で、心臓包絡面上における各々の位置（たとえば、約２０００以上のポイント）について対応する位相信号を与えることができる。
【００４３】
図１０は、複素平面にプロットされた心房信号の一例を示す。図中、実軸は信号自体に対応し、虚軸は当該信号のヒルベルト変換に対応する。図１１は、図１０と同じようにしてプロットされた心室細動信号の一例である。
【００４４】
たとえばＡＦまたはＶＦなどの不整脈の診断および／または治療を容易にするための評価を行うために、再構成されたエレクトログラム上から心臓包絡面におけるそれぞれの位置で計算された各々の位相信号について、対応する位相マップまたは他の情報を生成して表示することが可能である。対応する位相信号は、患者の心臓の（たとえば、心外膜表面などの）心臓包絡面に対応する３Ｄ幾何学形状のエレクトログラム位置ごとにマップ上に表示可能である。位相の性質は、−πがπと等しいものであるため、位相を可視化するのに利用されるカラーコードの範囲または他の尺度は、位相信号のこうした循環性を反映するように実装する必要がある。位相マップは、表面における位置各々について同時に、特定の時点における積分位相として表示可能である。また、このマップを動画位相マップとして表示し、表面全体に空間的にある位相の時間パターンを明示してもよい。３Ｄ表面は、本明細書で開示するようにして計算された位相信号に従って、表面の他の部分を表示するためのユーザー入力に応答して回転可能である。作成可能な異なる位相マップの例が、本明細書では図１４〜図２３に開示されている。
【００４５】
再び図８に戻ると、位相計算部ツール１００は、位相特異点の位置を計算するようプログラムされる位相特異点計算１０６を含んでもよい。位相特異点は、心房内または心室内のすべての位相が合流する位置として診断に有用となり得る。計算された位相がカラースケールによって空間的に表現される位相マップでは、位相特異点は、位相マップの対応するカラースケールのすべての色が同じ位置に生じる領域に対応する。たとえば、図１３の位相マップ２００において、位相特異点は「ＰＳ」を付した領域として示されている。これは、特定の患者の心臓の右心房前面上の位置である。
【００４６】
上述したように、位相特異点は、すべての位相が合流する空間位置に生じる。これは、位相勾配の積分として決定できる。たとえば、特定の時刻ｔにおける数学的定義から、以下の式６の場合に、位相特異点は位置ｘに生じる。
2019141619.tif 000016
ここで、Φ（ｔ）は、ｘ周囲の閉じた輪郭線Ｌについて計算される位相値である。
よって、式６を満たす位置および時刻は各々、特異点（singularity point）（すなわち、らせん波の旋回中心）とみなすことができる。ここで、心臓の表面（または心臓の他の膜）は、メッシュ構造で表すことが可能であり、式６の積分値は、位置Ｘ周囲の閉じた輪郭線に沿った位相変化の総和によって、たとえば以下のように推定可能である。
2019141619.tif 000017
ここで、図２４の近傍図に示すように、Ｎは、ｘに対して時計回りで順に並べられた、メッシュ表面におけるｘの近くにある位置であり、ｘiおよびｘjは、隣接する位置である。
図２４の例では、位置ｘは、それぞれ位相値Φ１、Φ２、Φ３、Φ４、Φ５を有する一組の近傍要素に囲まれている。らせん波の旋回中心の情報は、空間内でのそれらの位置と軌道を経時的に追跡するためのダイナミックマップ上に表現可能である。あるいは、らせん波の旋回中心の情報を時間平均し、エピソード全体について単一のマップを得るようにしてもよい。
【００４７】
位相計算部ツール１００は、位相データ１０４に基づいて１つ以上のフォーカス点の位置を決定するようにプログラムされる、フォーカス源計算部１０７を含んでもよい。このため、フォーカス源計算部１０７は、データ１０４によって表される位相信号に基づいて、局所的興奮を計算し、１つ以上の波面のフォーカス源を特定することが可能である。（図２４の近傍図を使用して）Ｘからその近傍要素に波面が広がっている場合に、特定の位置ｘをフォーカス源と定義できる。特定の時刻ｔにおいて、興奮波面は位相がπ／２と等しくなる場所に位置する。したがって、波面がＸ位置を通るときのＸの位相値が近傍要素の同じ瞬間における位相より早い場合、Ｘがフォーカスの位置である。図２４の近傍図では、位置Ｘは、以下が満たされる場合に、（フォーカス源計算部１０７によって）フォーカス源として決定できる。
2019141619.tif 000018
図２４の例では５つの近傍要素があるが、他の例では、近傍要素がいくつあってもよい。上記のフォーカス源の特定を利用して、たとえば１つ以上のフォーカスの発火数の指標をマップ状にグラフィック表示するなど、１つ以上の対応する可視画像を作成することが可能である。このようなマップは各々、心臓全体についての時空の一貫性のあるデータを与えるために、１つ以上の時間間隔にわたって計算可能である。いくつかの例では、（計算部１０７によって与えられる）フォーカスの活動と、（計算１０６および／または１０８によって与えられる）らせん波の活動の指標とを、別々のマップで可視化することが可能である。他の例では、フォーカスの活動とらせん波の活動をまとめ、１つのマップで同時に可視化して、細動機序の理解を容易にしてもよい。
【００４８】
三次元マップで位相特異点の位置を特定することに加えて、図８に示されるツール１００は、らせん波特定１０８およびらせん波特徴付け１１２も含んでもよい。らせん波特定１０８は、時間空間的に特定のらせん波の位置および軌道に関する統計量を計算するようにプログラム可能である。このような統計量は、たとえば１つ以上のエピソードをカバーできる特定のらせん波について単一のマップを作成するために、１つの時間間隔または複数の時間間隔にわたって、らせん波の旋回中心の位置を平均することを含んでもよい。たとえば、らせん波特徴付け機能１１２は、さまざまな統計量を計算するようにプログラム可能であり、たとえば、三次元の心臓包絡面における再構成されたエレクトログラム位置の各々についての興奮波面の統計量（たとえば、平均値、標準偏差、中央値）を含むことが可能である。
【００４９】
その代わりにまたはそれに加えて、らせん波特徴付け機能１１２は、心臓包絡面における各々の位置について２つの連続した興奮間の統計量（たとえば、平均値、標準偏差、中央値など）を計算可能である。さらに別の例として、らせん波特徴付け機能１１２は、心臓包絡面における各々の位置について、らせん波特徴付けの統計量（たとえば、平均値、標準偏差、中央値など）を計算するようにプログラム可能である。
【００５０】
らせん波特徴付け機能１１２は、ユーザー入力に応答して１つ以上の時間間隔を選択できる間隔選択部を含んでもよい。特定の間隔を特徴付けることに加えて、（たとえば、ユーザー入力に応答して）選択されている複数の間隔に対して他の行為をしてもよい。たとえば、らせん波特徴付けは、複数の間隔を組み合わせ、選択された複数の間隔についての位相マップを表示するようにプログラム可能である。らせん波特徴付け機能１１２は、さらに、ユーザー入力に応答して、ユーザーが選択している異なる時間間隔の間で、計算された統計量およびデータを比較できる。たとえば、間隔が異なれば長さも異なる場合がある。しかしながら、単位時間あたりの興奮波面の数を各々について計算し、比較して、らせん波の旋回中心の波面と数に関する、追加の情報を確認するといったことも可能である。このような他の情報は、らせん波の旋回中心の数、フォーカス源の興奮回数のみならず、心臓包絡面におけるフォーカス源の経時的な位置を含んでもよい。
【００５１】
また、ツール１００は、計算された位相信号からサイクル長を計算するようにプログラムされるサイクル長計算機能を含んでもよい。さらに、サイクル長の計算は、細動または他の不整脈が生じているとき、ユーザーが選択した時間間隔の間に変化し得るサイクル長の統計量を計算することを含んでもよい。前逆解処理および後逆解処理を含む、実行されたフィルタリングに基づいて、これらの信号は真の心房脱分極または心室脱分極に対応するものを反映している。たとえば、サイクル長の計算では、選択された時間間隔にわたる心房細動のサイクル長を計算できる。対応するマップを作成し、たとえば心臓の三次元表現上に表示される平均サイクル長に対応する、計算されたサイクル長の統計量を明示することが可能である。計算されたサイクル長から作成可能な、平均サイクル長をカラースケールで示すサイクル長マップの一例を、図１５に示す。
【００５２】
らせん波特徴付け１１２をはじめとして、ツール１００によって計算可能な他の情報には、心臓の３Ｄ表現上でらせん波の旋回中心の空間一貫性を明示するなどの目的で、積分位相勾配（たとえば式６または式７を参照のこと）を含んでもよい。たとえば、積分位相は、上限と下限の閾値を設定することにより、位相マップから計算可能であり、この閾値上で、表面における各々の位置（たとえば、複数の腔を同時に含む）について対応する積分位相を決定できる。対応する積分位相勾配マップを作成（たとえば、図１４を参照のこと）して、位相特異点を明示することが可能である。
【００５３】
らせん波特徴付け１１２のさらに別の例は、らせん波の旋回中心位置の時間周波数を計算することを含んでもよい。このような時間周波数は、電気活動データが取得されている単一の時間間隔について計算可能である。また、らせん波特徴付け機能１１２は、複数の間隔を組み合わせて、複数の間隔を合わせた全体の組に対する時間周波数を計算することも可能である。間隔を組み合わせることに加えて、らせん波特徴付け機能は、別々のマップの複数の間隔をユーザーが視覚的に比較できるようにするために、これらの間隔を比較するようにプログラム可能である。
【００５４】
図１２は、患者の診断および／または治療を行うために利用可能なシステム１５０の一例を示す。いくつかの例では、システム１５０は、診断または治療手順の一環として、心臓１５２をリアルタイムに評価できる。この診断または治療手順は、たとえば、医師が患者に対する治療行為をするためのパラメーター（たとえば、施術位置ならびに、治療の量とタイプなど）を決定しやすくするためのものである。たとえば、１つ以上の治療行為実施デバイス１５６が固定されたペーシングカテーテルまたはアブレーションカテーテルなどのカテーテルを、患者の心臓１５２に心内膜または心外膜で接触するように、身体１５４に挿入できる。当業者であれば、手順と治療のタイプに応じて変わり得る利用可能なさまざまなタイプと構成の治療行為実施デバイス１５６を理解し、評価できるであろう。
【００５５】
いくつかの例では、治療行為実施デバイス１５６は、治療システム１５８によって供給される電気信号（たとえば、無線周波数エネルギー）に応答して、組織を焼灼するための熱を発生させるよう構成されたアブレーションカテーテルの先端に配置された１つ以上の電極を含んでもよい。他の例では、治療行為実施デバイス１５６は、アブレーションを行うための冷却剤の送達（たとえば、冷凍アブレーション）、化学物質（たとえば、薬剤）の送達、超音波アブレーション、高周波アブレーション、あるいは、これらのアブレーションメカニズムまたは他のアブレーションメカニズムの組み合わせを行うよう構成されてもよい。さらに他の例では、治療行為実施デバイス１５６は、たとえば心臓のペーシング用などで、治療システム１５８によって供給される電気信号（たとえば、ペーシングパルス）に応答して電気刺激を行うためのペーシングカテーテルの先端に配置された１つ以上の電極を含んでもよい。治療システム１５８と、体内に配置される侵襲的な治療行為実施デバイス１５６によって、他のタイプの治療を行うこともできる。
【００５６】
治療システム１５８は、患者の身体１５４の外に位置してもよく、デバイス１５６によって行われる治療を制御するよう構成可能である。たとえば、治療システム１５８は、制御回路１６０を含む。この制御回路１６０は、デバイス（たとえば、電極）１５６と治療システム１５８との間に電気的に接続された導電性リンクを経由して電気信号を通信（たとえば、供給）できる。制御システム１６０は、心臓１５２の１箇所以上の位置に、１つまたは複数の電極１５４によって治療（たとえば、アブレーションまたは刺激）を施すための、デバイス１５６に供給される信号パラメーター（たとえば、電流、電圧、繰り返し数、トリガー遅延、検知用トリガー振幅）を制御できる。制御回路１６０は、治療パラメーターを設定し、自動、手動（たとえば、ユーザー入力）、または自動と手動との組み合わせ（たとえば、半自動制御など）で刺激を印加できる。１つ以上のセンサー（図示せず）も、治療システム１５８と通信してセンサー情報を戻すことができる。心臓１５２に対するデバイス１５６の位置は、画像診断技術（たとえば、蛍光透視、Ｘ線）、マッピングシステム１６２、直接的な目視などによって術中に決定し、追跡できる。このように、デバイス１５６の位置と治療パラメーターとを組み合わせて、対応する治療パラメーターデータを与えることができる。
【００５７】
セラピーシステム１５８での治療の前、治療の途中、および／または治療の後に、別のシステムまたはサブシステムを利用して、患者の電気生理学的情報を取得できる。図１２の例では、センサーアレイ１６４は、患者の活動を記録するのに利用できる１つ以上の電極を含む。一例として、センサーアレイ１６４は、（たとえば、心電図マッピング手順の一環として）患者の心臓に付随する電気活動を測定するために患者の胴の一部に分散される高密度配置の体表面センサー（たとえば、電極数が２００を上回るなど）に対応してもよい。使用可能な非侵襲センサーアレイの一例が、２００９年１１月１０日にファイルされた国際特許出願ＰＣＴ／ＵＳ２００９／０６３８０３（本明細書に援用する）に図示され、説明されている。他の配置の他の検知用電極をセンサーアレイ１６４として用いることもできる。このアレイは、患者の胴全体を覆わずに特定の目的で電気活動を測定するよう設計された、数を減らした電極の一組（たとえば、ＡＦおよび／またはＶＦの解析用に特別に設計された電極アレイ）であってもよい。
【００５８】
患者の体内に挿入されるデバイス１５６に１つ以上のセンサーが配置されてもよい。このような電極は、心腔の壁などの心内膜表面の電気活動をマッピングするためのセンサーアレイ１６４とともに利用できる。また、そのような電極は、心臓内でのデバイス１５６の位置の特定を助けるために利用でき、システム１５０によって作成される画像またはマップ上に表示されてもよい。
【００５９】
侵襲センサー、非侵襲センサーあるいは、侵襲センサーと非侵襲センサーの組み合わせを含む、患者の電気情報を取得するための、このような例としてあげた方法ではそれぞれ、１つまたは複数のセンサーアレイ１６４によって、検知された電気情報が、対応の測定システム１６６に供給される。測定システム１６６は、センサーアレイ１６４のセンサーによって検出される電気活動を示す対応の測定データ１７０を与えるための適当な制御・信号処理回路１６８を含んでもよい。測定データ１７０は、アナログ情報を含んでもよいし、デジタル情報を含んでもよい。
【００６０】
また、制御１６８は、電気活動を測定して測定データ１７０を供給するためのデータ取得プロセスを制御するようにも構成できる。測定データ１７０は、治療システムによって治療を行うのと同時に取得できる。この測定データ１７０は、たとえば、特定の治療を（たとえば治療パラメーターに従うなどして）適用することなどに反応して生じる心臓１５２の電気活動を検出するためのものである。たとえば、それぞれのデータ１７０と治療パラメーターとの時間的な関係をインデックス化するのに、適当なタイムスタンプを利用して、その評価と解析を容易にすることができる。
【００６１】
マッピングシステム１６２は、適当な処理および計算（たとえば、図１〜図４に関して開示したようなもの）を適用することによって、心臓１５２の電気活動に対応する測定データ１７０と幾何学形状データ１７２とを組み合わせ、対応する出力データ１７４を与えるようにプログラムされる。出力データ１７４は、心臓包絡面（たとえば、心臓１５２の表面上）の位相を表すまたは特徴付けるものであってもよい。
【００６２】
一例として、出力データ１７４は、患者の心臓１５２の心外膜表面について計算される位相データに基づいて（たとえば患者の身体１５４の表面に分散されたセンサー１６４を用いて非侵襲的に取得される電気データに基づいて）、位相マップまたは別の特徴付けに対応していてもよい。その代わりにまたはそれに加えて、出力データ１７４は、たとえば位相データを生成するために体表面センサー１６４を用いて記録されるものと同じ電気活動に基づいて、患者の心臓の表面について再構成される電位マップを含んでもよい。
【００６３】
いくつかの例では、測定システム１６６は心臓全体の電気活動を同時に測定できるため、得られる出力データ（たとえば、位相の特徴付けおよび／または他の心電図マップ）も、時空の一貫性がある形で心臓の同時データを表すことが可能である。出力データ／マップを計算する時間間隔は、ユーザー入力に基づいて選択可能である。その代わりにまたはそれに加えて、選択される間隔を、治療システム１５８による治療の適用に同期させることもできる。
【００６４】
マッピングシステム１６２は、１７８で示される信号処理方法を含む。信号処理方法１７８は、本明細書に開示の前逆解処理および後逆解処理（たとえば、図２の前処理５０、図３の前処理７０、図４の後処理８０）を含んでもよい。このため、信号処理は、ユーザー入力に応答して関連する信号成分を抽出することにより、電気的な測定データ１７０を前処理できる。たとえば、ユーザー入力は、マッピングを行う対象となる細動のタイプ（ＡＦまたはＶＦなど）を特定してもよい。選択されたタイプの位相評価に応じて、本明細書で開示するようにな外来信号成分と望ましくない振動とを除去するための適切なフィルタリングおよびベースライン除去を行うよう、信号処理を適宜構成することが可能である。
【００６５】
そうすると、エレクトログラム再構成１８０は、処理信号および幾何学形状データ１７２に基づいて、処理後の信号の逆解を計算し、対応する再構成されたエレクトログラムを与えることができる。よって、再構成されたエレクトログラムは、心臓包絡面での心電図活動に対応できるし、静的なもの（特定の時刻において三次元）であっても動的なもの（たとえば、時間の経過とともに変化する四次元マップ）であってもよい。
【００６６】
本明細書で開示する心臓包絡面（cardiac envelope）は、患者の心臓の実際の三次元表面に対応してもよい。この表面は、心外膜であっても心内膜であってもよい。その代わりにまたはそれに加えて、心臓包絡面は、患者の心臓の心外膜表面とセンサーアレイ１６４が配置されている患者の体表面との間に存在する、ある表面に対応してもよい。また、エレクトログラム再構成１８０によって利用される幾何学形状データ１７２は、患者の実際の幾何学形状、一般的なモデルまたはそれらの組み合わせに対応してもよい。
【００６７】
一例として、幾何学形状データ１７２は、患者について得られる画像データなど、患者の胴のグラフィック表現の形であってもよい。このような画像処理は、デジタル画像群から解剖学的特徴（１つ以上の臓器や他の構造を含む）を抽出およびセグメント化することを含んでもよい。また、たとえば、電極が患者に装着されている間に画像を取得し、適当な抽出およびセグメント化によって座標系で電極の位置を特定するなどして、センサーアレイ１６４における電極各々の位置が幾何学形状データ１７２に含まれていても構わない。こうしてセグメント化された画像データを、患者の着目領域を含む二次元または三次元のグラフィック表現に変換してもよい。
【００６８】
その代わりに、幾何学形状データ１７２は、数学モデルに対応してもよく、このモデルは、一般的なモデルであってもよいし、患者の画像データに基づいて構成されたものであってもよい。センサーアレイ１６４における電極の位置を含む、適切な解剖学的ランドマークまたは他のランドマークを幾何学形状データ１７２の中で特定し、電気的な測定データ１７０を合わせやすくして、そこでのインバース法を行いやすくしてもよい。このようなランドマークの特定は、（たとえば、画像編集ソフトウェアを使って人間が）手動で実施してもよいし、（たとえば、画像処理技術によって）自動的に実施してもよい。
【００６９】
さらに別の例として、幾何学形状データ１７２は、本明細書で説明する対応の表現を構成できるような、ほぼどのような画像診断技術（たとえば、ｘ線、超音波、コンピューター断層撮影、磁気共鳴イメージングなど）を用いても取得できる。このようなイメージングは、患者の電気的な測定データ１７０の生成に用いられる電気活動の記録と同時に行ってもよいし、イメージングを別に（たとえば、測定データを取得する前に）行ってもよい。
【００７０】
ＥＧＭ再構成法（たとえば、逆解）１８０によって計算される再構成されたエレクトログラムデータに対して、たとえば本明細書に開示の後処理方法（たとえば、図４の後逆解処理８０）に従って、さらに信号処理１７８をほどこしてもよい。信号処理機能１７８によって行われる後逆解処理に続いて、位相計算部１８２は、本明細書で開示するような各々の信号について（たとえば、１０００箇所を超えるポイントについて）位相データを計算することができる。得られる位相は、位相データとしてメモリーに格納し、マップ生成部１８４によって対応する高分解能の出力データ１７４を生成するのに利用可能である。
【００７１】
可視化エンジン１８６は、位相情報のグラフィック表現に対応する出力データ１７４を与えることができる。時間間隔や可視化して提示する情報のタイプの選択を含むことなど、可視化に付随するパラメーターは、対応する可視化ＧＵＩ１８８を介したユーザー入力に応答して選択可能である。このように、マッピングシステム１７４は、たとえば心電図位相マップ１９４を含む、ディスプレイ１９２で可視化エンジンによって対応するグラフィック出力として提供される対応する出力データ１７４を生成できる。
【００７２】
細動間隔ごとに位相マップおよび位相特徴付けマップを作成するマッピングシステム１６２に加えて、たとえば興奮マップ、優位周波数マップなどを含む、他のタイプの心電図マッピングを利用してもよい。たとえば、ディスプレイ１９２は、位相マップデータを対応の興奮または優位周波数マップと同時に表示して、ＡＦまたはＶＦの診断および評価を容易にするための、１つ以上の領域を含んでもよい。
【００７３】
図１３から図２３は、本明細書に開示のシステムおよび方法システムによって作成可能な位相マップおよび関連のグラフィック表現の例を示す。
【００７４】
図１３は、心臓のグラフィック表現上にレンダリングされる位相マップ２００の一例を示す。マップ２００は、マップに表示される三次元表面における各々の位置について計算される位相（たとえば、±πの間）に対応するカラースケール２０２を含む。図１３の例では、らせん波の旋回中心がＰＳで示されている。これは、すべての位相が合流する単一の位置（たとえば、二次元空間領域などの領域）に対応する。他の例では、特定の時点において、心臓の表面に２つ以上のらせん波が存在する場合もあり、本明細書に開示の手法では、特定の時間間隔について心臓の表面全体から同時に得られる電気データに基づいて位相データが計算されるため、上記のような位相特異点の各々を単一の３Ｄグラフィック表示で明示できることは、自明であろう。
【００７５】
図１４は、５００ミリ秒の時間間隔にわたって計算された、ＡＦ症例についてのらせん波の旋回中心の位置を示す心臓のグラフィックマップ（たとえば、積分位相勾配マップ）の一例を示す。図１４の例では、らせん波の旋回中心は、各３Ｄマップにおいて２１０と２１２で示されている。このマップの暗色の領域は、らせん波の旋回中心がアンカーされると決定されなかった位置に対応する（たとえば、式６または式７による）。図１４のマップを作成するのに用いられるデータは、位相特異点またはらせん波特定機能（たとえば、図８の１０６または１０８に対応する実行可能なコード）によって決定できる。
【００７６】
図１５は、マッピングシステムによって作成可能なサイクル長マップの一例を示す。図１５の例では、サイクル長は、５００ミリ秒の時間間隔にわたって生じる位相についての平均サイクル長として計算されており、この時間間隔は、ユーザーが選択可能である。平均サイクル長は、心臓の三次元幾何学形状にプロットされている。カラーコードスケール２２０を利用して、サイクル長に対して計算された統計量（たとえば、中央値、標準偏差）を特徴付けることが可能であるし、カラースペクトル（たとえば１２０ミリ秒から約２３０ミリ秒の範囲）における優位周波数を特徴付けることも可能である。図１５の例では、マップ２１６および２１８は各々、三次元空間内に同じ心臓について異なる向きで示されている。これはたとえば、心臓の異なる表面領域についての時間的に一貫性のあるサイクル長を示すためのものである。
【００７７】
図１６は、ある時間間隔のあいだ、心臓包絡面（たとえば、心外膜表面全体などであってもよい、複数の腔など）に対して計算された積分位相に対応する位相マップを示す心電図マップ２３０の一例を示す。図１６の例では、マップで提示される情報は、心臓包絡面におけるらせん波の旋回中心（領域２３２で示される）の空間的な一貫性を表すことができる。マップ２３０における位相特性は、計算された位相が、すでに確立されているあらかじめ定められた位相閾値よりも大きい場合にデータをプロットするなどして、１つ以上の間隔におよぶ積分位相の総和を表現できる。この閾値は、提示されている情報の感度を変更するために（たとえば、図１２のＧＵＩ１８８を用いて）ユーザーが選択するものであってもよい。
【００７８】
図１７は、ある時間間隔にわたる複数の異なる場合における動画またはムービーのフォーマットを明示する位相マップの一例を示す。図１７において、１７Ａで示す上側の部分は、１０５０ミリ秒から１０８０ミリ秒の位相マップを明示している。１７Ｂで示す下側の部分では、２３００ミリ秒から２６６０ミリ秒の時間の範囲。経時的な位相マップの動画レンダリングは、ユーザーがらせん波の動きや、らせん波間の相互作用を特定する一助となり得る。
【００７９】
計算された位相を、図１７に示すものと同じように静的または動画的にプロットすることに加えて、本明細書で開示するようならせん波の軌道や他のらせん波特性を特定するなどの目的で、他の可視化マッピング方法を利用してもよい。
【００８０】
一例として、心臓包絡面（たとえば、心臓の心外膜表面に対応するもの）上の各ポイントを用いて、対応する表面を、図１８の例に示すような円柱の上に投影することができる。対応する心臓包絡面または心臓の表面は、たとえば特定の解析と関連する解剖学的位置に基づいて、複数の着目領域にセグメント化することが可能である。図１９は、異なるカラーコードで複数の空間領域にセグメント化されているポイントを示す。（たとえば、図１９からの）セグメント化された各々の空間領域は、図１８の投影された円柱のポイントから二次元平面上にマッピング可能である。
【００８１】
別の言い方をすると、図１８〜図１９は、図２０に示されるプロットにどのように到達するかの中間工程を示している。すなわち、図１８および図１９は、らせん波の軌道と特徴付けの解析のために、幾何学形状の各領域を、どのように分割し、２Ｄ表面に投影できるかを示している。図１８は、幾何学形状の投影プロセスを示し、図１９は、２Ｄ表面に対して心臓の各領域をどのように特定できるかを示している。２Ｄ投影後、これらの領域は、単一の領域として表示することもできるし、複数の領域の組み合わせとして表示することもできる。
【００８２】
図２０は、１つまたは複数のセグメント化された領域が、２Ｄ平面２５１に投影されたプロットを示し、たとえば図１８および図１９に明示された表面の空間領域の解析に基づいて決定可能なものである。他の形態の解析、たとえばらせん波の軌道追跡などもまた、このような解析から導出することができる。さらに図２０に示されているのは、時間に対応する垂直軸である。プロット２５０から、特定されるらせん波の旋回中心２５２、２５４および２５６の各々について、いくつかのらせん波の特性（たとえばキラリティーや輪運動速など）を計算できる。この特徴付けの方法では、らせん波の旋回中心を特定でき、ある時間間隔にわたって旋回中心の軌道を測定できる。たとえば、各旋回中心は二次元の表面投影で特定されるため、特定の表面は、その近くにある表面に対して空間的に表現される。このため、プロットならびに、時間の経過にともなう対応するらせん波の動きとらせん波の軌道を、それぞれの表面投影で追跡することが可能である。らせん波間の相互作用や、表面領域におけるそれぞれのらせん波間でのらせん波の断片化も特定可能である。計算を重ねることによって、図２０に示される表面投影２５１は空間的に関連しており、これを利用して、対応の表面投影を移動するらせん波の旋回中心を特定できる。
【００８３】
図２０の例では、時間の経過にともなう３つのらせん波の旋回中心が、２５２、２５４、２５６で明示されている。らせん波２５６は、開始時（たとえば、０から約４００ミリ秒の時点まで）には空間的にかなり静止してとどまっているが、サイクル時間の後のほうになると、空間投影２５１のまわりを回転する。３００から３５０ミリ秒付近では、らせん波２５２および２５４は、２５８で示すように、それぞれの領域をまたがって相互作用する。
【００８４】
らせん波の旋回中心が空間的に特定され、相互作用が決定されたら、特定されたらせん波の旋回中心軌道を、３Ｄ表面の対応する領域に戻して投射し、経時的に変化するマップで対応するらせん波の動きを可視化することが可能である（たとえば、図１７を参照のこと）。らせん波の旋回中心が心臓の解剖のグラフィックマップで空間的に特定されたら、らせん波の旋回中心が生じる組織を焼灼する際に使用できるよう、カテーテルのデザインを作成することが可能である。このカテーテルのデザインは、電極の数と空間的配置の特定を含んでもよい。いくつかの例では、カテーテルのデザインは、あらかじめ構成された複数のカテーテルのうちの１つを、型番などによって特定できる。他の例では、カテーテルのデザインは、１つ以上のらせん波の旋回中心が特定されている組織領域を焼灼するための専用の電極群を有する構成と寸法で作られた、カスタム構成に対応していてもよい。
【００８５】
また、らせん波相互作用マップを作成し、グラフィカルな方法でらせん波の相互作用を可視化してもよい。たとえば、心臓の３Ｄ可視画像上に異なるカラーコードまたは他のマーキングをおき、そのような相互作用を特定してもよい。らせん波の旋回中心のスピン方向を示すそれぞれの円または四角の中にシンボルを示すことによって、らせん波相互作用マップでキラリティーを見ることもできる。輪運動速も、マップ上のらせん波の相互作用を示す線の色の強度としてマップに表示可能である。また、輪運動速を１秒あたりの回転として表示し、同じく３Ｄグラフィックマップに提示可能な対応するカラースケールにフィットさせてもよい。
【００８６】
たとえば図２１および図２２に示すような、追加のマップを作成してもよい。図２１には、たとえば図２０に示す計算およびらせん波の旋回中心間の相互作用に基づいて作成可能であるような、らせん波の旋回中心相互作用マップが示されている。図２２は、グラフィックマップに提示して、キラリティー（たとえば、特定のらせん波の回転方向）を示すことが可能な情報、ならびに、特定されている各らせん波の輪運動速を示すためにマップに適用可能な色の変化の、別の例を示す。
【００８７】
ＡＦおよびＶＦに加えて、本明細書に開示の位相計算は、たとえば図２３の例に示されるような、心房粗動についてのグラフィック表現を提供するのにも利用できる。図２３では、小さな積分位相内の位相だけがカラーコード化されている。これはたとえば、心臓の弁に対する位相の反時計回りの回転に対応するように、心房粗動の伝播経路の追跡を容易にするために利用できる。一例として、心房粗動の場合について図２３に示すような、計算された位相Φ（ｘ）＝一定であるという状況を明示するために、１つ以上のマップまたは動画マップを示すことが可能である。この図のカラー版では、たとえば、条件Φ（ｘ）＝π／２を紫でグラフィカルに表現でき、他のすべての条件を灰色（または他の異なる色）で表現できる。他のタイプの不整脈について、同様の診断情報を可視化することが可能であり、これはたとえば、心房細動または心室細動の位相が一定である場合を可視化するためのものであってもよい。
【００８８】
上記の構造的または機能的な説明に鑑みて、本発明の一部を、方法、データ処理システムまたはコンピュータープログラム製品として実施してもよいことは、当業者であれば自明であろう。したがって、本発明のこのような一部は、完全にハードウェアである実施形態、完全にソフトウェアである実施形態、あるいは、図２５のコンピューターシステムで図示し、説明するような、ソフトウェアとハードウェアとを組み合わせた実施形態のうちの、いずれの形態をとってもよい。さらに、本発明の一部が、媒体上にコンピューター読取り可能なプログラムコードを有する、コンピューター使用可能な記憶媒体上のコンピュータープログラム製品であってもよい。好適なコンピューター読取り可能な媒体であれば、どのようなものを利用してもよく、一例として、静的記憶装置および動的記憶装置、ハードディスク、光学記憶装置および磁気記憶装置があげられるが、これらに限定されるものではない。
【００８９】
本明細書では、方法、システム、およびコンピュータープログラム製品のブロック図を参照して、本発明の特定の実施形態についても説明してきた。図示のブロックならびに、図示のブロックの組み合わせを、コンピューターで実行可能な命令によって実装してもよいことは、理解できよう。これらのコンピューターで実行可能な命令を、汎用コンピューター、専用コンピューター、あるいは、マシンを製造するための他のプログラム可能なデータ処理装置（または装置と回路の組み合わせ）の１つ以上のプロセッサーに対して、プロセッサーによって実行する命令が１つまたは複数のブロックに示される機能を実装するように供給してもよい。
【００９０】
これらのコンピューターで実行可能な命令は、コンピューターまたは他のプログラム可能なデータ処理装置を、コンピューター読取り可能なメモリーに格納された命令が、フローチャートの１つまたは複数のブロックに示される機能を実装する命令を含む製品につながるように、特定の方法で機能させることのできるコンピューター読取り可能なメモリーに格納されてもよい。また、コンピュータープログラムの命令を、コンピューターまたは他のプログラム可能なデータ処理装置にロードして、このコンピューターまたは他のプログラム可能な装置上で一連の動作ステップを実施させ、コンピューターまたは他のプログラム可能な装置上で実行する命令によって、フローチャートの１つまたは複数のブロックに示される機能を実装するためのステップが得られるように、コンピューターで実装されるプロセスを生成してもよい。
【００９１】
この点に鑑みて、図２５は、本発明の１つ以上の実施形態を実行するのに用いることのできるコンピューターシステム３００の一例を示す。これらの実施形態は、たとえば、センサーデータの取得および処理、画像データの処理のみならず、心臓の電気活動の解析に関連した、変換されたセンサーデータおよび画像データの解析を含む。コンピューターシステム３００は、１つ以上のネットワーク接続された汎用コンピューターシステム、組込型のコンピューターシステム、ルーター、スイッチ、サーバー装置、クライアント装置、さまざまな中間装置／ノード、またはスタンドアロンのコンピューターシステムに実装可能である。また、コンピューターシステム３００は、たとえば、携帯情報端末（ＰＤＡ）、ラップトップコンピューター、ポケットベルなどのさまざまなモバイルクライアントに、それが十分な処理能力を持つという前提で、実装可能である。
【００９２】
コンピューターシステム３００は、処理装置３０１と、システムメモリー３０２と、さまざまなシステム構成要素（システムメモリーを含む）を処理装置３０１に接続するシステムバス３０３とを含む。デュアルマイクロプロセッサーなどのマルチプロセッサーアーキテクチャを、処理装置３０１として使用してもよい。システムバス３０３は、多岐にわたるバスアーキテクチャーのうちのいずれかを用いる、メモリーバスまたはメモリーコントローラー、周辺機器用バス、ローカルバスをはじめとする、さまざまなタイプのバス構造のうち、どのようなものであってもよい。システムメモリー３０２は、リードオンリーメモリー（ＲＯＭ）３０４と、ランダムアクセスメモリー（ＲＡＭ）３０５とを含む。コンピューターシステム３００内の要素間で情報を転送しやすくする基本ルーチンを含むＲＯＭ３０４に、基本入出力システム（ＢＩＯＳ）３０６を常駐させることが可能である。
【００９３】
コンピューターシステム３００は、ハードディスクドライブ３０７と、たとえばリムーバブルディスク３０９との間で読み書きするための磁気ディスクドライブ３０８と、たとえばＣＤ−ＲＯＭディスク３１１を読み取ったり、他の光学媒体との間で読み書きしたりするための光ディスクドライブ３１０と、を含んでもよい。ハードディスクドライブ３０７、磁気ディスクドライブ３０８、光ディスクドライブ３１０はそれぞれ、ハードディスクドライブインターフェース３１２、磁気ディスクドライブインターフェース３１３、光ドライブインターフェース３１４を介して、システムバス３０３に接続されている。ドライブと、これに関連するコンピューター読取り可能な媒体が、コンピューターシステム３００用のデータ、データ構造、コンピューターで実行可能な命令を格納する、不揮発性の記憶装置となる。上述したコンピューター読取り可能な媒体の説明では、ハードディスク、リムーバブル磁気ディスク、ＣＤについて言及したが、この動作環境で、磁気カセット、フラッシュメモリーカード、デジタルビデオディスクなど、さまざまな形態でのコンピューター読取り可能な他のタイプの媒体を使用してもよい。さらに、このような媒体のいずれも、本明細書の１つ以上の部分を実現するためのコンピューター読取り可能な命令を含んでもよい。
【００９４】
オペレーティングシステム３１５、１つ以上のアプリケーションプログラム３１６、他のプログラムモジュール３１７、プログラムデータ３１８をはじめとして、多数のプログラムモジュールを、ドライブおよびＲＡＭ３０５に格納してもよい。アプリケーションプログラムおよびプログラムデータは、本明細書にて図示し、説明するような１つ以上のセンサーからの電気データを取得し、処理し、表示するようにプログラムされる機能および方法を含んでもよい。アプリケーションプログラムおよびプログラムデータは、本明細書で開示するように信号を処理し、位相データを計算するようにプログラムされる機能と方法を含んでもよい。また、アプリケーションプログラムおよびプログラムデータは、本明細書で開示するような位相マップまたは他の心電図マップを作成するようにプログラムされる機能と方法を含んでもよい。
【００９５】
ユーザーは、１つ以上の入力装置３２０を介して、コマンドおよび情報をコンピューターシステム３００に入力してもよい。入力装置３２０は、たとえば、ポインティングデバイス（たとえば、マウス、タッチスクリーン）、キーボード、マイク、ジョイスティック、ゲームパッド、スキャナーなどである。たとえば、ユーザーは、ドメインモデルを編集または修正するのに、入力装置３２０を用いることができる。これらの入力装置および他の入力装置３２０は、システムバスに接続された対応のポートインターフェース３２２を介して処理装置３０１に接続されることが多いが、パラレルポート、シリアルポート、またはユニバーサルシリアルバス（ＵＳＢ）などの他のインターフェースを介して接続されることもある。１つまたは２つ以上の出力装置３２４（たとえば、ディスプレイ、モニター、プリンター、プロジェクターまたは他のタイプの表示装置）も、ビデオアダプターなどのインターフェース３２６を介してシステムバス３０３に接続される。
【００９６】
コンピューターシステム３００は、リモートコンピューター３２８などの１台以上のリモートコンピューターとの論理接続を使用して、ネットワーク環境で動作してもよい。リモートコンピューター３２８は、ワークステーション、コンピューターシステム、ルーター、ピアデバイス、または他の一般的なネットワークノードであってもよく、一般に、コンピューターシステム３００に関して説明した多くの要素またはすべての要素を含む。３３０で概略的に示す論理接続は、ローカルエリアネットワーク（ＬＡＮ）およびワイドエリアネットワーク（ＷＡＮ）を含んでもよい。
【００９７】
ＬＡＮネットワーク環境で用いる場合、コンピューターシステム３００は、ネットワークインターフェースまたはアダプター３３２を介してローカルネットワークに接続されてもよい。ＷＡＮネットワーク環境で用いる場合、コンピューターシステム３００は、モデムを含んでもよいし、ＬＡＮ上の通信サーバーに接続されてもよい。モデム（内蔵式であっても外付けであってもよい）は、適当なポートインターフェース経由でシステムバス３０３に接続できる。ネットワーク環境では、コンピューターシステム３００との関連で示したアプリケーションプログラム３１６またはプログラムデータ３１８、あるいはこれらの一部を、リモートメモリー記憶装置３４０に格納してもよい。
【００９８】
以上説明してきたものは、本発明の例および実施形態である。もちろん、本発明を説明する目的で、構成要素または方法論の考えられるすべての組み合わせを説明するのは不可能であるが、本発明の多くの他の組み合わせおよび入れ替えが可能であることは、当業者であれば認識するであろう。したがって、本発明は、添付の請求項の範囲に入るこのような変更、修正および改変をすべて包含することを意図している。請求項において、特に明記しないかぎり、冠詞「ａ」は、「１またはそれより多く」を示す。
【図１】
2019141619.tif 000019
【図２】
2019141619.tif 000020
【図３】
2019141619.tif 000021
【図４】
2019141619.tif 000022
【図５】
2019141619.tif 000023
【図６】
2019141619.tif 000024
【図７】
2019141619.tif 000025
【図８】
2019141619.tif 000026
【図９】
2019141619.tif 000027
【図１０】
2019141619.tif 000028
【図１１】
2019141619.tif 000029
【図１２】
2019141619.tif 000030
【図１３】
2019141619.tif 000031
【図１４】
2019141619.tif 000032
【図１５】
2019141619.tif 000033
【図１６】
2019141619.tif 000034
【図１７】
2019141619.tif 000035
【図１８】
2019141619.tif 000036
【図１９】
2019141619.tif 000037
【図２０】
2019141619.tif 000038
【図２１】
2019141619.tif 000039
【図２２】
2019141619.tif 000040
【図２３】
2019141619.tif 000041
【図２４】
2019141619.tif 000042
【図２５】
2019141619.tif 000043
【手続補正書】
【提出日】令和1年5月7日(2019.5.7)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
媒体。
【外国語明細書】
2019141619000001.pdf