(57)【要約】【課題】  類似の波形を有するＥＣＧ信号を検出する。【解決手段】  ヒト被験者の１つの心拍に対して取られた、心電図（ＥＣＧ）信号の初期セットを選択することであって、このセットが、被験者の不整脈のテンプレートとして使用されるそれぞれの波形を有する、ことと、このヒト被験者の後続心拍に対して取られたＥＣＧ信号の後続セットを受け取ることと、を含む、方法。この方法は更に、初期セットと後続セットとの間の相互相関解析を行うことにより、初期セットの幾何学的形状と後続セットの幾何学的形状との間の適合度の尺度である相関係数を生成することと、相関係数がある閾値係数を超えている場合に、後続心拍が不整脈により生じたものであるとして受け入れることと、を含む。【選択図】  図１
2018023783.tif 000002
【特許請求の範囲】
【請求項１】
方法であって、
ヒト被験者の１つの心拍に対して取られた心電図（ＥＣＧ）信号の初期セットを選択することであって、前記セットが、前記被験者の不整脈のテンプレートとして使用されるそれぞれの波形を有する、ことと、
前記ヒト被験者の後続心拍に対して取られたＥＣＧ信号の後続セットを受け取ることと、
前記初期セットと前記後続セットとの間の相互相関解析を行うことにより、前記初期セットの幾何学的形状と前記後続セットの幾何学的形状との間の適合度の尺度である相関係数を生成することと、
前記相関係数がある閾値係数を超えている場合、前記後続心拍が前記不整脈により生じたものであるとして受け入れることと、
を含む、方法。
【請求項２】
前記初期セットの発生の初期セット推定時間周辺の初期セット時間関心領域（初期セットＷＯＩ）を線で描くことと、前記後続セットの発生の後続セット推定時間周辺の後続セット時間関心領域（後続セットＷＯＩ）を線で描くことと、前記初期セットＷＯＩ内の前記初期セットの信号と、前記後続セットＷＯＩ内の前記後続セットの信号と、を使用して、前記相互相関解析を実施することと、を含む、請求項１に記載の方法。
【請求項３】
前記初期セットＷＯＩ及び前記後続セットＷＯＩが、共通の時間幅を有する、請求項２に記載の方法。
【請求項４】
前記相互相関解析を実施する前に、前記初期セットと前記後続セットとの間に位相シフトを適用することを含む、請求項１に記載の方法。
【請求項５】
前記位相シフトを繰り返し変化させて最大相関係数を決定することと、該最大相関係数が前記閾値係数を超えている場合、前記後続心拍が前記不整脈により生じたものであるとして受け入れることと、を含む、請求項４に記載の方法。
【請求項６】
前記後続心拍が前記不整脈により生じたものであるとして受け入れることが、前記不整脈の位置の指標を、前記ヒト被験者の心臓マップ内に取り込むことを含む、請求項１に記載の方法。
【請求項７】
前記ＥＣＧ信号が体表（ＢＳ）ＥＣＧ信号である、請求項１に記載の方法。
【請求項８】
前記ＥＣＧ信号が心内（ＩＣ）ＥＣＧ信号である、請求項１に記載の方法。
【請求項９】
装置であって、
ヒト被験者の１つの心拍に対して取られた、前記被験者の不整脈のテンプレートとして使用されるそれぞれの波形を有する心電図（ＥＣＧ）信号の初期セットを受け取り、かつ、前記ヒト被験者の後続心拍に対して取られたＥＣＧ信号の後続セットを受け取るよう構成されている、一組の電極と、
前記初期セットと前記後続セットとの間の相互相関解析を行うことにより、前記初期セットの幾何学的形状と前記後続セットの幾何学的形状との間の適合度の尺度である相関係数を生成するよう構成され、前記相関係数がある閾値係数を超えている場合、前記後続心拍が前記不整脈により生じたものであるとして受け入れるよう構成されている、プロセッサと、
を含む、装置。
【請求項１０】
前記初期セットの発生の初期セット推定時間周辺の初期セット時間関心領域（初期セットＷＯＩ）を線で描くことと、前記後続セットの発生の後続セット推定時間周辺の後続セット時間関心領域（後続セットＷＯＩ）を線で描くことと、前記初期セットＷＯＩ内の前記該初期セットの信号と、前記後続セットＷＯＩ内の前記後続セットの信号と、を使用して、前記相互相関解析を実施することと、を含む、請求項９に記載の装置。
【請求項１１】
前記初期セットＷＯＩ及び前記後続セットＷＯＩが、共通の時間幅を有する、請求項１０に記載の装置。
【請求項１２】
前記プロセッサが、前記相互相関解析を実施する前に、前記初期セットと前記後続セットとの間に位相シフトを適用するよう構成されている、請求項９に記載の装置。
【請求項１３】
前記プロセッサが、前記位相シフトを繰り返し変化させて最大相関係数を決定するよう構成され、かつ、前記最大相関係数が前記閾値係数を超えている場合、前記後続心拍が前記不整脈により生じたものであるとして受け入れるよう構成されている、請求項１２に記載の装置。
【請求項１４】
前記後続心拍が前記不整脈により生じたものであるとして受け入れることが、前記不整脈の位置の指標を、前記ヒト被験者の心臓マップ内に取り込むことを含む、請求項９に記載の装置。
【請求項１５】
前記ＥＣＧ信号が体表（ＢＳ）ＥＣＧ信号である、請求項９に記載の装置。
【請求項１６】
前記ＥＣＧ信号が心内（ＩＣ）ＥＣＧ信号である、請求項９に記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
（関連出願の相互参照）
本出願は、２０１６年８月１０日に出願された米国仮出願第６２／３７２，９６９号の利益を主張し、これは、参照により本明細書に組み込まれる。
【０００２】
（発明の分野）
本発明は、概ね心電図（ＥＣＧ）信号に関するものであり、より詳細には類似の波形を有するＥＣＧ信号の検出に関するものである。
【背景技術】
【０００３】
不整脈を生じている心内腔の領域を正確にマッピングするためには、その特定の不整脈を表している信号（すなわち心拍）のみを捕えることが重要となる。異所性心拍、組織の機械的刺激、及び同じ周期長さの代替活性化パターンによる波形の不整脈変化などの影響による信号は、無視されるべきである。そのような信号からの結果をマップに組み込むと、局所活性化マップに不正確さが生じ、不整脈の視覚化の変形により、不整脈メカニズムを明確に識別することが困難になる。
【０００４】
本特許出願に参照により組み込まれる文献は、いずれかの用語がこれらの組み込まれた文献において、本明細書で明確に又は暗示的になされる定義と矛盾する形で定義されている場合には、本明細書中の定義のみが考慮されるべきである点を除いて、本出願の一部とみなされるものとする。
【発明の概要】
【課題を解決するための手段】
【０００５】
本発明の一実施形態は、１つの方法を提供するものであり、その方法は、
ヒト被験者の１つの心拍に対して取られた心電図（ＥＣＧ）信号の初期セットを選択することであって、このセットは、被験者の不整脈のテンプレートとして使用されるそれぞれの波形を有する、ことと、
ヒト被験者の後続心拍に対して取られたＥＣＧ信号の後続セットを受け取ることと、
初期セットと後続セットとの間の相互相関解析を行うことにより、初期セットの幾何学的形状と後続セットの幾何学的形状との間の適合度の尺度である相関係数を生成することと、
相関係数がある閾値係数を超えている場合、後続心拍が不整脈により生じたものであるとして受け入れることと、を含む。
【０００６】
開示されている一実施形態は、初期セットの発生の初期セット推定時間周辺の初期セット時間関心領域（初期セットＷＯＩ）を線で描くことと、後続セットの発生の後続セット推定時間周辺の後続セット時間関心領域（後続セットＷＯＩ）を線で描くことと、この初期セットＷＯＩ内の初期セットの信号と、この後続セットＷＯＩ内の後続セットの信号と、を使用して、相互相関解析を実施することと、を含む。典型的に、初期セットＷＯＩ及び後続セットＷＯＩは、共通の時間幅を有する。
【０００７】
更に開示されている一実施形態は、相互相関解析を実施する前に、初期セットと後続セットとの間に位相シフトを適用することを含む。更に開示されている実施形態は、更に、位相シフトを繰り返し変化させて最大相関係数を決定することと、この最大相関係数がある閾値係数を超えている場合、後続心拍が不整脈により生じたものであるとして受け入れることと、を含む。
【０００８】
また更に開示されている一実施形態において、その後続心拍が不整脈により生じたものであるとして受け入れることは、不整脈の位置の指標を、ヒト被験者の心臓マップ内に取り込むことを含む。
【０００９】
ＥＣＧ信号は、体表（ＢＳ）ＥＣＧ信号であり得る。
【００１０】
代替的に又は付加的に、ＥＣＧ信号は心内（ＩＣ）ＥＣＧ信号であり得る。
【００１１】
本発明の一実施形態により、更に装置が提供され、この装置は、
ヒト被験者の１つの心拍に対して取られた、被験者の不整脈のテンプレートとして使用されるそれぞれの波形を有する心電図（ＥＣＧ）信号の初期セットを受け取り、かつ、このヒト被験者の後続心拍に対して取られたＥＣＧ信号の後続セットを受け取るよう構成されている、一組の電極と、
その初期セットと後続セットとの間の相互相関解析を行うことにより、初期セットの幾何学的形状と後続セットの幾何学的形状との間の適合度の尺度である相関係数を生成するよう構成され、かつ、この相関係数がある閾値係数を超えている場合に、この後続心拍が不整脈により生じたものであるとして受け入れるよう構成されている、プロセッサと、を含む。
【００１２】
以下の本開示の実施形態の詳細な説明を図面と併せ読むことで、本開示のより完全な理解が得られるであろう。
【図面の簡単な説明】
【００１３】
【図１】本発明の一実施形態による侵襲性医療処置の概略図である。
【図２】本発明の一実施形態によるＥＣＧ波形合致アルゴリズムの入力及び出力を示す概略ブロック図である。
【図３】本発明の一実施形態によるＥＣＧ波形合致アルゴリズムを示す概略ブロック図である。
【図４】本発明の一実施形態によるアルゴリズムのブロック動作を示す概略図である。
【図５】本発明の一実施形態によるアルゴリズムのブロック動作を示す概略図である。
【図６】本発明の一実施形態によるアルゴリズムのブロック動作を示す概略図である。
【発明を実施するための形態】
【００１４】
概説
本発明の実施形態は、ＥＣＧ信号において同じ波形を表すすべての心拍を識別することを目的とした、ＥＣＧ波形合致アルゴリズムを提供する。このアルゴリズムは、体表（ＢＳ）（典型的に１２誘導）信号及び／又は心内（ＩＣ）信号に適合する。このアルゴリズムは、ＥＣＧ信号の波形パターンを入力として受け取り、連続ＥＣＧ信号中で同じ波形を探す。
【００１５】
本発明のユーザーは、特定の注釈周囲で関心領域（ＷＯＩ）を規定することにより、入力波形パターンを選択する。このアルゴリズムは、選択されたパターンの波形を、入ってくるＥＣＧ信号の波形と比較する。所定の重み付き相関閾値範囲内にある心拍は、同じ波形を表すとみなされる。このアルゴリズムは、心拍信号が心臓内のプローブにより取得される際に、リアルタイムで作動する。
【００１６】
このアルゴリズムの結果を使用することにより、合致した心拍の発生源である心臓領域を、心臓マップ上に自動的に示すことができる。
【００１７】
本発明の一実施形態は、ヒト被験者の１つの心拍に対して取られた心電図（ＥＣＧ）信号の初期セットを選択することであって、このセットが、被験者の不整脈のテンプレートとして使用されるそれぞれの波形を有する、ことと、このヒト被験者の後続心拍に対して取られたＥＣＧ信号の後続セットを受け取ることと、を含む方法を提供する。この方法は更に、初期セットと後続セットとの間の相互相関解析を行うことにより、初期セットの幾何学的形状と後続セットの幾何学的形状との間の適合度の尺度である相関係数を生成することを含む。この相関係数がある閾値係数を超えている場合、後続心拍が不整脈により生じたものであるとして受け入れられる。
【００１８】
図１は、本発明の一実施形態による、装置２０を使用する侵襲性医療処置の概略図である。本処置は、医療専門家２２により行われ、一例として、本明細書の以下の説明における処置は、ヒトの患者２６の心臓２４からの心内心電図（ＩＣ  ＥＣＧ）信号の取得を含むと想定される。本発明の実施形態は、ＩＣ  ＥＣＧ又はＢＳ（体表）ＥＣＧ信号のいずれかを解析するが、以下の説明を簡潔かつ明快にするために、別途記載のない限り、ＩＣ  ＥＣＧ信号が解析されるものとして想定される。
【００１９】
ＩＣ  ＥＣＧ信号を取得するために、医療専門家２２は、患者の内腔内に事前に位置付けられているシース３０の中へプローブ２８を挿入する。シース３０は、シースの遠位端３４から出たプローブの遠位端３２が患者の心臓内に入り心臓の組織に接触できるように位置付けられている。
【００２０】
プローブ２８は、患者の心臓内に挿入でき、かつ典型的には磁気追跡システム及び／又はインピーダンス測定システムを用いて追跡可能であるような、任意のタイプのカテーテルを含み得る。例えば、プローブ２８は、ラッソーカテーテル、シャフト状カテーテル、又はｐｅｎｔａＲａｙカテーテル（Ｂｉｏｓｅｎｓｅ  Ｗｅｂｓｔｅｒ（Ｄｉａｍｏｎｄ  Ｂａｒ，ＣＡ）製造）、あるいはこれらのカテーテルに概ね類似のカテーテルを含み得る。Ｂｉｏｓｅｎｓｅ  Ｗｅｂｓｔｅｒは、本発明の実施形態に使用可能な磁気追跡システム及びインピーダンス測定システムも製造している。
【００２１】
プローブ２８は１つ以上の電極３６を含み、これらは、本明細書に記述されるアルゴリズムの実施において、装置２０に含まれるプロセッサ４０が使用するＥＣＧ信号を取得するのに使用される。プロセッサ４０は、中央処理装置としての働きに加えて、リアルタイムノイズ低減回路４４（典型的には、フィールドプログラマブルゲートアレイ（ＦＰＧＡ）、次いでアナログデジタル（Ａ／Ｄ）信号変換集積回路４６として構成される）を含み得る。このプロセッサはＡ／Ｄ回路４６からの信号を別のプロセッサへと伝達することができ、かつ、本明細書に記述されるアルゴリズムを実施するようプログラムすることができる。
【００２２】
プロセッサ４０は、装置の操作コンソール６０内に位置する。コンソール６０は、医療専門家２２がプロセッサと通信するために使用するコントロール６２を備える。処置中に、プロセッサ４０は、ＥＣＧ信号を取得し、また本明細書に記述されるアルゴリズムを実施するために、モジュールバンク７０内のＥＣＧモジュール６６と通信を行う。
【００２３】
ＥＣＧモジュール６６は電極３６からＥＣＧ信号を受け取る。一実施形態において、この信号は、モジュール６６内で、低ノイズ前置増幅器６８を通り、ローパス及びハイパスフィルタ７１Ａ、７１Ｂを通って、主増幅器７２に伝達される。モジュール４３６は更にアナログデジタル変換器（ＡＤＣ）７４を含み、これは、ＥＣＧ信号のデジタル化された値をプロセッサ４０に伝達し、これによって、本明細書に記述されるアルゴリズムをプロセッサが実施する。典型的に、プロセッサ４０は前置増幅器６８、フィルタ７１Ａ、７１Ｂ、増幅器７２、及びＡＤＣ  ７４の動作を制御する。
【００２４】
簡潔化のために、図１は、電極３６からの信号を受信するためのチャネルを１つを有するものとして、ＥＣＧモジュール６６を示す。しかしながら、このモジュールは典型的に、図示されるものと実質的に同様のチャネルを複数含むことが理解されよう。例えば、モジュール６６はそのようなチャネルを１２個含んでよく、これは１２個の体表電極からの信号を受け取るのに使用され得る。
【００２５】
ＥＣＧモジュール６６により、プロセッサ４０は、電極３６が受け取ったＥＰ（電気生理学的）信号（本明細書で参照されるＥＣＧ信号を含む）を取得及び解析することができるようになる。この信号は典型的に電圧−時間グラフとして専門家２２に示され、これはディスプレイ画面８０上にリアルタイムで更新される。
【００２６】
プロセッサ４０及びモジュールバンク７０のソフトウェアは、例えば、ネットワークで、電子的形態でプロセッサにダウンロードすることができる。あるいは又は更には、ソフトウェアは、例えば、光学的、磁気的、又は電子的記憶媒体のような非一時的有形媒体上で提供され得る。
【００２７】
装置２０を操作するために、モジュールバンク７０は典型的に、上述のＥＣＧモジュール以外のモジュール（例えば、プローブ２８の遠位端をプロセッサが追跡できるようにする１つ以上の追跡モジュール）を含む。簡潔化のために、かかる他のモジュールは、図１に図示されていない。すべてのモジュールは、ハードウェア要素並びにソフトウェア要素を含み得る。
【００２８】
電極４１１が取得するＥＣＧ信号をディスプレイ画面８０に表示することに加えて、本明細書に記載されるアルゴリズムの結果も、アルゴリズムのユーザーに対してディスプレイ画面上に示され得る。例えば、結果は心臓２４のマップ８２に組み込むことができる。
【００２９】
図２は、本発明の一実施形態によるＥＣＧ波形合致アルゴリズムの入力及び出力を示す概略ブロック図である。アルゴリズムはプロセッサ４０により実施され、入力として受け入れる波形合致フィルタ１００として機能する：
・ＥＣＧモジュール６６の複数のチャネルにより上述のように処理された、ＥＣＧ信号のセット１０２。
・プロセッサ４０により計算された、信号の参照注釈１０４。信号の注釈は、信号の発生の想定時間である。一実施形態において、この注釈は、選択された１つのＥＣＧ信号に対する最大の正の値の発生時間に対応する。参照注釈（正の値、負の値、最大の負の傾き、最大の正の傾き）及びＩＣ  ＥＣＧ信号（発生時間は典型的に、この信号を生成する心筋の部分の活性化時間に対応する）について、いくつかの基準オプションが存在する。上述又は他の基準に対応する注釈のためのＥＣＧ信号を選択するための基準は、医療専門家２２により規定され得る。医療専門家は更に、解析する信号を取得するのに使用されるＥＣＧチャネルも選択する。ＢＳ  ＥＣＧ信号については、典型的に１２個のチャネルがある。ＩＣ  ＥＣＧ信号については、チャネルの数は、使用される電極３６の数に対応する。
・波形パターン１０６。これは、医療専門家２２が選択するＥＣＧ信号のセットであり、医療専門家が規定した信号の注釈周辺の関心領域（ＷＯＩ）で、特定の時点で捕捉される。ＷＯＩは波形合致アルゴリズムの期間を画定する。波形パターン１０６は、他のＥＣＧ信号を比較するテンプレートとして機能し、このパターンは、本明細書ではテンプレートとも呼ばれることもある。
・医療専門家２２が設定する相関閾値１０８；これは、ある心拍がその波形パターンに合致するかどうかをアルゴリズムが判定するのに使用される。
【００３０】
アルゴリズムの出力は以下の通りである。
・心拍状況１１０。すなわち、波形パターン１０６と同じ波形を有するとして受け入れられた心拍、又は入力パターンとは異なる波形を有するとして却下された心拍。一実施形態において、受け入れられた心拍の位置が、ＥＣＧ信号を生成する心臓のマップに組み込まれる。
・信号セット１０２の各チャネルに対する相関スコア１１２（いくつかの実施形態においては、スコアは医療専門家２２に提示されなくともよい）。
・各心拍について、すべてのチャネルに対して計算された、重み付き相関スコア１１４。
【００３１】
図３は、ＥＣＧ波形合致アルゴリズムを示す概略ブロック図であり、図４、５、及び６は、本発明の一実施形態によるアルゴリズムのブロック動作を示す概略図である。例として、プロセッサ４０はアルゴリズムを実施すると想定される。他の実施形態において、プロセッサはスタンドアロンプロセッサであってよく、及び／又は、典型的にコンピュータを作動させている汎用プロセッサであってもよい。
【００３２】
アルゴリズムの第１の工程（これは「単一チャネル相関」ブロック１２０に対応する）において、プロセッサは、入ってくる心拍の各チャネルについて、医療専門家２２が規定するＷＯＩ期間内で、格納されている波形パターン１０６との相関解析を実施する。図４は、この相関解析がどのように実施されるかを示している。図４に示すように、プロセッサは入力として以下を使用する：
・図２を参照して上述された波形パターン１０６（パターンｉ）。
・試験されるＥＣＧ信号（ＥＣＧ  ｉ，ｊ）。試験される信号は、医療専門家２２により規定された波形パターンのＷＯＩ幅に対応するＷＯＩ時間幅を有する。ＷＯＩの時間的位置は、信号のリアルタイム注釈（プロセッサにより計算される）を含むよう選択される。
【００３３】
ｉは、パターンのチャネルを規定する数値的指数であり（典型的に、ＢＳ  ＥＣＧの場合はｉ＝１、２、．．．１２）、ｊは、ＥＣＧ信号の注釈の位置を規定する数値的指数である。
【００３４】
プロセッサは各チャネルについて、以下の式に従って相関係数を計算する：
【００３５】
【数１】
2018023783.tif 000003
式中、
ｘはテンプレート参照ＥＣＧデータのサンプル値であり、
【数２】
2018023783.tif 000004
はテンプレート参照ＥＣＧデータの平均値であり、
ｙは試験される現在の心拍ＥＣＧデータのサンプル値であり、
【数３】
2018023783.tif 000005
は試験される現在の心拍ＥＣＧデータの平均値であり、
ｋは、ＥＣＧ信号のどのデータサンプルが解析されているかを規定する数値的指数である。例えば、ＷＯＩが−５０ｍｓｅｃ（参照注釈の前）から＋７０ｍｓｅｃ（参照注釈の後）までの１２０ｍｓｅｃであり、毎ｍｓｅｃごとにサンプリングを行うとき、ｋは１２０個のサンプルに対して１２０個の値のセットとなる。
【００３６】
式（１）により実施される相関解析は、テンプレート参照ＥＣＧデータの幾何学的形状又は形を現在の心拍ＥＣＧデータと比較するものであることが理解されよう。相関（ｘ，ｙ）が高い値、すなわち、合致に近い場合は、テンプレートと現在の心拍の２つの幾何学的形状が類似していることを意味する。
【００３７】
アルゴリズムの第２の工程（これは図３の「全体重み付き相関」ブロック１２４に対応する）において、プロセッサは、第１の工程で計算された（すなわち式（１）に従って計算された）相関係数の値を用いて、特定の心拍の全体相関を計算する。
【００３８】
図５は、この相関解析がどのように実施されるかを示している。図５に示すように、プロセッサは入力として以下を使用する：
・相関スコア、すなわち、試験される心拍の各チャネルについての、式（１）の出力。試験される心拍は、現在の注釈周辺のＷＯＩ内にあるＥＣＧ信号（その特定のチャネルのＥＣＧ信号）である。
・試験されるＥＣＧ信号（心拍）。
・図１を参照して上述された波形パターン１０６（パターンｉ）。
【００３９】
更に図５に示すように、プロセッサは、試験されるＥＣＧ信号の絶対最大振幅Ａｉ，ｊ、及び波形パターンの絶対最大振幅Ｂｉを計算する。
【００４０】
プロセッサはＡｉ，ｊとＢｉとの合計を重み付けとして使用して、式（２）に従い全体相関を計算する：
【００４１】
【数４】
2018023783.tif 000006
式中、Ｃｏｒｒｉ，ｊは式（１）で計算された相関係数であり、Ｎは解析されるＥＣＧチャネルの数である。ＢＳ信号の場合、Ｎは典型的には１２である。
【００４２】
式（２）により計算される全体相関係数は、波形パターンの位相に対する、試験されるＥＣＧ信号の位相に依存する。
【００４３】
アルゴリズムの第３の工程（これは図３の「位相シフト」ブロック１２８に対応する）において、プロセッサは、波形パターンの位相に対して、試験されるＥＣＧ信号の位相を繰り返し変更する。プロセッサは入力として以下を使用する：
・式（２）から得た全体相関の値。
・医療専門家２２により設定される相関閾値１０８（図２）。この閾値は０〜１であり得、典型的な値は０．９である。
【００４４】
図６は、ブロックの繰り返しセットにおいて、単一チャネル相関ブロック１２０、全体重み付き相関ブロック１２４、及び位相シフトブロック１２８により実施される繰り返しプロセスを説明している。図６に示すように、各繰り返しにおいて、「シフトされた単一チャネル相関」ブロック１２０’及び「シフトされた全体相関」ブロック１２４’において、プロセッサは上述の最初の２つの工程を繰り返す。この繰り返し中に、プロセッサは、式（２）の結果としての全体相関の最大値を決定する。
【００４５】
図６は、解析される心拍の注釈から測定して±４０ｍｓｅｃの時間枠で、位相シフトの繰り返しが１ｍｓｅｃごとに実施されているのを示す。指数ｋは、評価されている位相シフトを規定し、ｋ＝｛−４０，．．．０，．．．＋４０｝である。繰り返しのたびに、比較ブロック１３０で、シフトされた全体相関の値が、それまでの最大相関と比較され、もし比較１３０が「はい」を返す場合、全体相関は、更新ブロック１３４で更新される。
【００４６】
「いいえ」が返された場合は、コントロールは引き続き比較ブロック１３８に進み、繰り返すべき指数ｋの値がまだあるかどうかをチェックする。まだある場合、ｋは増分ブロック１４２で増加され、ｋの新たな値が信号ブロック１４６でＥＣＧ信号に適用され、フローチャートはブロック１２０’に戻る。
【００４７】
繰り返しが完了すると、コントロールは引き続き最終比較ブロック１５２に進み、ここで、ブロックの繰り返しセットの出力である、ブロック１３４内の全体相関の最大値を、入力閾値と比較する。この比較結果が「はい」を返す場合、その心拍は、波形パターンと同じ不整脈を表すと想定される。この場合、プロセッサ４０はこの心拍情報を、心臓のマップ８２の集合的情報に加えることができる（図１）。この比較結果が「いいえ」を返す場合、この心拍は波形パターンとは異なる不整脈を表すと想定され、よって、この情報はマップの集合的情報には追加されない。
【００４８】
上に述べた実施形態は例として挙げたものであり、本発明は上記に具体的に示し説明したものに限定されない点が理解されよう。むしろ、本発明の範囲は、上記されている種々の特徴の組み合わせ及び部分的組み合わせ、並びに前述の説明を読むことで当業者が想起するであろう、先行技術に開示されていない変形及び変更の両方を含む。
【００４９】
〔実施の態様〕
（１）  方法であって、
ヒト被験者の１つの心拍に対して取られた心電図（ＥＣＧ）信号の初期セットを選択することであって、前記セットが、前記被験者の不整脈のテンプレートとして使用されるそれぞれの波形（morphologies）を有する、ことと、
前記ヒト被験者の後続心拍に対して取られたＥＣＧ信号の後続セットを受け取ることと、
前記初期セットと前記後続セットとの間の相互相関解析を行う（performing a cross-correlation）ことにより、前記初期セットの幾何学的形状と前記後続セットの幾何学的形状との間の適合度の尺度である相関係数を生成することと、
前記相関係数がある閾値係数を超えている場合、前記後続心拍が前記不整脈により生じたものであるとして受け入れることと、
を含む、方法。
（２）  前記初期セットの発生の初期セット推定時間周辺の初期セット時間関心領域（初期セットＷＯＩ）を線で描くことと、前記後続セットの発生の後続セット推定時間周辺の後続セット時間関心領域（後続セットＷＯＩ）を線で描くことと、前記初期セットＷＯＩ内の前記初期セットの信号と、前記後続セットＷＯＩ内の前記後続セットの信号と、を使用して、前記相互相関解析を実施することと、を含む、実施態様１に記載の方法。
（３）  前記初期セットＷＯＩ及び前記後続セットＷＯＩが、共通の時間幅を有する、実施態様２に記載の方法。
（４）  前記相互相関解析を実施する前に、前記初期セットと前記後続セットとの間に位相シフトを適用することを含む、実施態様１に記載の方法。
（５）  前記位相シフトを繰り返し変化させて最大相関係数を決定することと、該最大相関係数が前記閾値係数を超えている場合、前記後続心拍が前記不整脈により生じたものであるとして受け入れることと、を含む、実施態様４に記載の方法。
【００５０】
（６）  前記後続心拍が前記不整脈により生じたものであるとして受け入れることが、前記不整脈の位置の指標を、前記ヒト被験者の心臓マップ内に取り込むことを含む、実施態様１に記載の方法。
（７）  前記ＥＣＧ信号が体表（ＢＳ）ＥＣＧ信号である、実施態様１に記載の方法。
（８）  前記ＥＣＧ信号が心内（ＩＣ）ＥＣＧ信号である、実施態様１に記載の方法。
（９）  装置であって、
ヒト被験者の１つの心拍に対して取られた、前記被験者の不整脈のテンプレートとして使用されるそれぞれの波形を有する心電図（ＥＣＧ）信号の初期セットを受け取り、かつ、前記ヒト被験者の後続心拍に対して取られたＥＣＧ信号の後続セットを受け取るよう構成されている、一組の電極と、
前記初期セットと前記後続セットとの間の相互相関解析を行うことにより、前記初期セットの幾何学的形状と前記後続セットの幾何学的形状との間の適合度の尺度である相関係数を生成するよう構成され、前記相関係数がある閾値係数を超えている場合、前記後続心拍が前記不整脈により生じたものであるとして受け入れるよう構成されている、プロセッサと、
を含む、装置。
（１０）  前記初期セットの発生の初期セット推定時間周辺の初期セット時間関心領域（初期セットＷＯＩ）を線で描くことと、前記後続セットの発生の後続セット推定時間周辺の後続セット時間関心領域（後続セットＷＯＩ）を線で描くことと、前記初期セットＷＯＩ内の前記該初期セットの信号と、前記後続セットＷＯＩ内の前記後続セットの信号と、を使用して、前記相互相関解析を実施することと、を含む、実施態様９に記載の装置。
【００５１】
（１１）  前記初期セットＷＯＩ及び前記後続セットＷＯＩが、共通の時間幅を有する、実施態様１０に記載の装置。
（１２）  前記プロセッサが、前記相互相関解析を実施する前に、前記初期セットと前記後続セットとの間に位相シフトを適用するよう構成されている、実施態様９に記載の装置。
（１３）  前記プロセッサが、前記位相シフトを繰り返し変化させて最大相関係数を決定するよう構成され、かつ、前記最大相関係数が前記閾値係数を超えている場合、前記後続心拍が前記不整脈により生じたものであるとして受け入れるよう構成されている、実施態様１２に記載の装置。
（１４）  前記後続心拍が前記不整脈により生じたものであるとして受け入れることが、前記不整脈の位置の指標を、前記ヒト被験者の心臓マップ内に取り込むことを含む、実施態様９に記載の装置。
（１５）  前記ＥＣＧ信号が体表（ＢＳ）ＥＣＧ信号である、実施態様９に記載の装置。
【００５２】
（１６）  前記ＥＣＧ信号が心内（ＩＣ）ＥＣＧ信号である、実施態様９に記載の装置。
【図１】
2018023783.tif 000007
【図２】
2018023783.tif 000008
【図３】
2018023783.tif 000009
【図４】
2018023783.tif 000010
【図５】
2018023783.tif 000011
【図６】
2018023783.tif 000012
【外国語明細書】
2018023783000001.pdf