Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP6253879B2 - 心内心電計信号の正確な時間注釈 - Google Patents
[go: Go Back, main page]

JP6253879B2 - 心内心電計信号の正確な時間注釈 - Google Patents

心内心電計信号の正確な時間注釈 Download PDF

Info

Publication number
JP6253879B2
JP6253879B2 JP2012247240A JP2012247240A JP6253879B2 JP 6253879 B2 JP6253879 B2 JP 6253879B2 JP 2012247240 A JP2012247240 A JP 2012247240A JP 2012247240 A JP2012247240 A JP 2012247240A JP 6253879 B2 JP6253879 B2 JP 6253879B2
Authority
JP
Japan
Prior art keywords
signal
time
derivative
function
unipolar
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2012247240A
Other languages
English (en)
Other versions
JP2013103130A (ja
Inventor
ウラジミール・ルビンシュタイン
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Biosense Webster Israel Ltd
Original Assignee
Biosense Webster Israel Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Biosense Webster Israel Ltd filed Critical Biosense Webster Israel Ltd
Publication of JP2013103130A publication Critical patent/JP2013103130A/ja
Application granted granted Critical
Publication of JP6253879B2 publication Critical patent/JP6253879B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/327Generation of artificial ECG signals based on measured signals, e.g. to compensate for missing leads
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/367Electrophysiological study [EPS], e.g. electrical activation mapping or electro-anatomical mapping
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7278Artificial waveform generation or derivation, e.g. synthesizing signals from measured signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/35Detecting specific parameters of the electrocardiograph cycle by template matching

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Physiology (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Cardiology (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Description

本発明は、全般的には、信号解析に関し、具体的には、医療処置の間に生成される信号の解析に関する。
心臓などの、患者の身体器官から生成される電気信号は、典型的には雑音性である。これらの信号は典型的には、患者に対する医療処置の間に測定され、信号に対する雑音は通常は、複数の要因によって引き起こされる。要因のうちの一部は、電極と器官の一区域との接触の移動若しくは変化などの、人為的結果、測定されている領域の近位で発生している他の信号による干渉、身体器官の比較的高いインピーダンス、及び生成されている信号の固有の変化である。
身体器官からの信号に対する雑音の影響を低減するためのプロセスは、それゆえに有利である。
本発明の一実施形態は、信号を解析するための方法を提供し、この方法は、
時変心内電位信号を感知することと、
既定の振動波形に対する、時変心内電位信号の適合を見出すことと、
その適合に応答して、信号の注釈時間を推定することとを含む。
開示される実施形態では、時変心内電位信号は、単極性信号を含む。典型的には、既定の振動波形は、単一の極大値、単一の極小値、及び極小値と極大値とを分離する単一の変曲を有する、単一の完全な振動を含む。
代替的な実施形態では、既定の振動波形は、ガウス関数の一次微分を含む。典型的には、一次微分は、非対称因子によって歪められる。
開示される別の実施形態では、時変心内電位信号は、双極性信号を含む。既定の振動波形は、第1の単一の完全な振動と、第2の単一の完全な振動との差異を含み得る。典型的には、第1の単一の完全な振動は、第1の単一の極大値、第1の単一の極小値、及び第1の極大値と第1の極小値とを分離する第1の変曲を含み、第2の単一の完全な振動は、第2の単一の極大値、第2の単一の極小値、及び第2の極大値と第2の極小値とを分離する第2の変曲を含む。
第1の単一の完全な振動と、第2の単一の完全な振動とは、時間差によって分離することができる。この時間差は、双極性信号を生成する電極の、空間的分離の関数とすることができる。あるいは、又はそれに加えて、この時間差は、活性化波の伝播方向に対する、電極の配向の関数とすることができる。
開示される更なる実施形態では、既定の振動波形は、第1のガウス関数の一次微分と、第2のガウス関数の一次微分との差異を含む。典型的には、第1のガウス関数の一次微分は、第1の非対称因子によって歪められ、第2のガウス関数の一次微分は、第2の非対称因子によって歪められる。
開示される、更なる実施形態では、時変心内電位信号は、時間差を間に有する3つ又は4つ以上の単極性信号を含み、活性化波の伝播方向は、その時間差の関数である。典型的には、それぞれの位置を有するそれぞれの電極が、3つ又は4つ以上の単極性信号を生成し、それぞれの位置を、関数のパラメーターとすることができる。
本発明の一実施形態により、信号を解析するための装置が更に提供され、この装置は、
時変心内電位信号を感知するように構成されるセンサーと、
既定の振動波形に対する、時変心内電位信号の適合を見出し、その適合に応答して、信号の注釈時間を推定するように構成される、プロセッサとを含む。
本開示は、以下の本開示の実施形態の「発明を実施するための形態」を、以下の図面と併せ読むことによって、より完全に理解されるであろう。
本発明の一実施形態による、心電計(ECG)解析システムの概略図。 本発明の一実施形態による、ECG解析システムによって処理される、典型的なECG信号の概略的グラフ。 本発明の実施形態による、ECG信号に適合させるために使用される等式によって作成された概略的グラフ。 本発明の実施形態による、ECG信号に適合させるために使用される等式によって作成された概略的グラフ。 本発明の一実施形態による、心内信号を解析する際の工程を示す流れ図。 本発明の一実施形態による、図1のシステムによって得られた結果を示す概略的グラフ。
概論
本発明の一実施形態は、単極性又は双極性とすることができる、「生の」心内信号、又はフィルター処理された心内信号を処理するための方法を提供する。典型的には、この処理は、既定の波形に心内信号を適合させ、生の信号からではなく、その適合された信号から、信号の注釈時間を導き出すことを含む。
典型的には、単極性信号は、単一の完全な振動を表す等式に適合される。双極性信号は、典型的には時間差によって分離される、2つの単一の完全な振動の差異を表す等式に適合させることができる。一部の実施形態では、単一の完全な振動は、ガウス関数の微分に対応する。非対称因子を、その微分に適用することができ、一部の実施形態では、非対称因子は、ガウス関数に対応する。
本発明者らは、生の信号又はフィルター処理された信号から直接判定される注釈時間と比較して、生の信号又はフィルター処理された信号を、既定の等式に適合させ、その適合された信号から注釈時間を測定することにより、注釈時間の変動が低減されることを見出した。
システムの説明
ここで図1を参照すると、この図1は、本発明の一実施形態による、心電計(ECG)解析システム20の概略図である。システム20は、ヒト患者の器官内部に位置決めされる1つ又は2つ以上の電極から、少なくとも1つの、典型的には複数個の電気信号を受信する。典型的には、それらの信号は、器官内の1つ又は2つ以上のプローブ上に定置される、多数の電極から受信される。例えば、心臓に対する侵襲的手技の間、1つ又は2つ以上の電極を有する第1のプローブを、心臓の基準領域内に位置決めして、その領域から、基準ECG信号を感知するために使用することができる。複数の電極を有する第2のプローブを使用して、心臓の他の領域から、他のECG信号を検出して記録することができる。
簡潔性及び明瞭性のために、以下の説明は、特に既述のない限り、単一のプローブ24を使用して、心臓34からの電気信号を感知する調査手順を想定する。更には、プローブの遠位端32は、2つの実質的に同様の電極22A、22Bを有するものと想定される。電極22A、22Bは、本明細書では電極22と称される場合がある。当業者であれば、この説明を、1つ又は2つ以上の電極を有する複数のプローブ関して、並びに、心臓以外の器官によって生成される信号に関して、適応させることが可能であろう。
典型的には、プローブ24は、システム20のユーザー28によって実行されるマッピング手順の間に、被験者26の身体内に挿入される、カテーテルを含む。本明細書での説明では、ユーザー28は、一例として、医療専門家であると想定される。この手順の間、被験者26は、接地電極23に取り付けられると想定される。一部の実施形態では、電極29を、心臓34の領域内で被験者26の皮膚に取り付けることができる。
システム20は、メモリ44と通信する処理装置42を含む、システムプロセッサ40によって制御することができる。プロセッサ40は典型的には、動作制御機器38を含むコンソール46内に搭載される。制御機器38は、典型的には、専門家28がプロセッサと対話するために使用する、マウス又はトラックボールなどの、ポインティングデバイス39を含む。プロセッサは、メモリ44内に格納される、プローブ経路指定モジュール30及びECGモジュール36を含めた、ソフトウェアを使用して、システム20を動作させる。ECGモジュール36は、基準ECGサブモジュール37及びマップECGサブモジュール41を含み、それらの機能は以下で説明される。プロセッサ40によって実行される動作の結果は、ディスプレイ48上で専門家に提示され、このディスプレイ48は、典型的には、操作者に対するグラフィックユーザーインターフェース、電極22によって感知されるECG信号の視覚的表現、及び/又は調査されている間の心臓34の画像を提示する。ソフトウェアは、例えば、ネットワークを介して、電子的形態で、プロセッサ40にダウンロードすることができ、あるいは、又はそれに加えて、磁気メモリ、光メモリ、又は電子メモリなどの、非一時的な有形媒体上に提供及び/又は格納することができる。
ECGモジュール36は、電極22からの電気信号を受信するように結合される。このモジュールはまた、電極29のうちの1つ又は2つ以上からの信号を受信するように結合することもできる。ECGモジュールは、信号を解析するように構成され、ディスプレイ48上に、標準的なECG形式で、典型的には、時間と共に変動するグラフ式表現で、解析の結果を提示することができる。
プローブ経路指定モジュール30は、プローブが被験者26の内部にある間に、プローブ24の諸区域を追跡する。経路指定モジュール30は、典型的には、被験者26の心臓内部での、プローブ24の遠位端32の、場所及び配向の双方を追跡する。一部の実施形態では、モジュール30は、プローブの他の区域を追跡する。この経路指定モジュールは、当該技術分野において既知の、プローブを追跡するための任意の方法を使用することができる。例えば、モジュール30は、被験者の付近で、磁場発信器を動作させることにより、その発信器からの磁場は、追跡されているプローブの諸区域内に配置された追跡コイルと、相互作用することができる。磁場と相互作用するコイルは、モジュールに伝送される信号を生成し、モジュールは、その信号を解析して、コイルの場所及び配向を判定する。(簡潔性のために、そのようなコイル及び発信器は、図1には示されない。)Biosense Webster(Diamond Bar,CA)によって生産される、Carto(登録商標)システムは、そのような追跡方法を使用する。あるいは、又はそれに加えて、経路指定モジュール30は、電極23、電極29、及び電極22の間のインピーダンス、並びに、プローブ上に配置することができる他の電極に対するインピーダンスを測定することによって、プローブ24を追跡することができる。(この場合には、電極22及び/又は電極29は、ECG信号及び追跡信号の双方を提供することができる。)Biosense Websterによって生産される、Carto3(登録商標)システムは、磁場発信器及びインピーダンス測定の双方を、追跡に関して使用する。
図2は、本発明の一実施形態による、システム20によって処理される、典型的なECG信号の概略的グラフを示す。グラフ100、102は、「生の」(すなわち、未処理の)双極性心内ECG信号の、例示的な、電位対時間のプロットを示す。これらの信号は、電極が心臓壁に接触する間の、電極22Aと電極22Bとの電位差に由来するものと想定される。当該技術分野において既知であるように、心内ECG信号は雑音性であり、その雑音は、典型的には、線放射、他の電気機器との近接性、及び、患者の筋収縮(心筋を除く)などの患者26に由来する他の電気供給源などの、多数の要因によって生成される。雑音は、典型的には、生の信号から注釈時間の定量的測定を行なう際に、問題を引き起こす。
例えば、信号の「R」ピークの時間を含む、注釈時間Tが必要とされる場合があり、この時間は、信号の立上がりから測定される。グラフ100は、Tが、約30msであると測定されることを示すが、グラフ102は、Tが、約25msであると測定されることを示す。グラフに示されるように、Tの測定値は変動する。
上述のように、グラフ100、102は、電極22Aと電極22Bとの差信号によって作り出される、双極性のグラフを示す。各電極22A又は電極22B上の信号は、共通の基準電極に対して測定される場合、単極性信号であることにより、この双極性信号は、2つの単極性信号の差異として考慮することができる。基準電極は、接地電極23、及び/又は皮膚電極29のうちの1つ又は2つ以上、及び/又は心臓と接触する1つ又は2つ以上の他の電極などの、任意の便宜的な電極とすることができる。
図3及び図4は、本発明の実施形態による、ECG信号に適合させるために使用される等式によって作成された概略的グラフを示す。本発明の実施形態は、図2に示すECG信号などの信号に、既定の等式を適合させる。この等式は、既定の振動波形、典型的には、単一の完全な振動の形態である波形、すなわち、実質的にゼロの信号レベルを有する開始点及び終了点を有し、その2つの点の間の全ての電気的活動を包含する波形に対応する。典型的には、単一の完全な振動のグラフは、単一の極小値及び単一の極大値を有する。極大値と極小値とは、単一の変曲によって分離することができる。
一部の実施形態では、また本明細書で例示されるように、信号に適合される既定の等式は、ガウス関数の一次微分から導かれ、非対称因子によって歪められる。
それゆえ、電極22A又は電極22Bから受信される単極性ECG信号に関しては、プロセッサ40は、以下の等式(1)によって与えられる一般形態を有する等式を、信号に適合させる。
Figure 0006253879
式中、V単極性(t)は、時間tでの、電極で測定される、変動する単極性電位信号を表し、
は、時間t=0に対しての、信号の時間的変位である。tは、活性化波が電極の位置を通過する際の時間に対応し、
Aは、信号の振幅であり、
は、信号の非対称性を規定するパラメーターであり、
wは、信号の幅を規定するパラメーターである。
等式(1)の検証は、この等式によって提供される非対称因子が、ガウス関数に対応することを示す。それゆえ、等式(1)は、ガウス関数と、ガウス関数の一次微分とを合計する。
以下の説明では、ti1、A1、ts1、及びwのパラメーターはまた、等式(1)の単極性適合パラメーターとも総称される。
グラフ110、112、及びグラフ114(図3)は、等式(1)によって作り出される、波形に対するt及びwのパラメーターの値の効果を示す。簡潔性のために、各グラフの縦座標及び横座標の単位は、任意であると想定される。グラフ110に示すように、t=0に関しては、グラフは、(3,0)で対称の中心を有する、2回転対称を有する。(換言すれば、グラフの平面内の180°の回転の下で、このグラフは、それ自体へと変換する。)グラフ112は、t=3の正値に関しては、グラフが非対称になることを示す。この非対称性は、tの増大と共に増大する。
グラフ114によって示されるように、wの値は、グラフの全幅を変化させるものであり、wの値を増大させることにより、幅が低減される。
ECG信号が双極性信号である場合には、そのECG信号は、電極22A上の単極性信号V単極性(t)と、電極22B上の単極性信号V単極性(t)との差異によって生成されるものと、想定することができる。これらのような双極性信号に関しては、プロセッサは、等式(1)から導かれる等式(2)を、信号に適合させる。
Figure 0006253879
式中、V双極性(t)は、時間tでの、電極で測定される、変動する双極性電位信号を表し、
単極性(t)、V単極性(t)は、V及びVとも称され、等式(1)に関して上記で規定されるようなものであり、
i1、ti2は、V、Vの時間的変位であり、
A1、A2は、V、Vの振幅であり、
s1、ts2は、V、Vの非対称性を規定し、
、wは、V、Vの幅を規定する。
双極性信号に関しては、時間差Δt=ti1−ti2が存在し、この時間差は、2つの単極性信号、V単極性(t)とV単極性(t)との、時間的変位の差異に等しい。2つの単極性信号の間の時間差は、典型的には、双極性信号を生成する2つの電極の空間的分離、及び活性化波の伝播方向に対する電極の配向の関数である。それゆえ、2つの電極の場合には、活性化波の伝播方向の少なくとも1つの成分を、単極性信号の時間差から判定することができる。3つ以上の電極に関しては、その3つ以上の電極によって検出される、それぞれの単極性信号の間の時間差、並びに電極の位置により、典型的には、伝播方向の複数の成分を見出すことが可能になる点が理解されるであろう。それらの複数の成分から、活性化波の伝播方向(単に一成分ではない)を推定することができる。
以下の説明では、ti1、ti2、A1、A2、ts1、ts2、及びw、wのパラメーターはまた、等式(2)の双極性適合パラメーターとも総称される。
グラフ120、122、及びグラフ124(図4)は、等式(2)の適用を示す。グラフ120及びグラフ122は、電圧対時間の、2つの単極性の等式のグラフであり、それぞれ、t=3及びt=4.5の時間的変位(任意の単位での)、並びに4及び2の幅を有する。グラフ124は、それらの2つの式の差異のグラフであり、Δt=4.5−3=1.5の時間差を有する、双極性電圧対時間の関数を示す。
生成される心内単極性信号及び心内双極性信号は、特に、信号を測定するために使用される電極の位置に応じて変化する。生成される信号はまた、測定されている心臓の状態、すなわち、その心臓が、健常な方式で機能しているか、又は病的な方式で機能しているかに応じても変化する。
心臓が、特定の欠陥により病的である場合には、その心臓はまた、健常な心臓の心内信号とは異なる、標準的な心内信号を作り出す(同様の相違は、皮膚ECG信号、すなわち、身体表面信号を使用する診断で使用することができる)。特定の欠陥の場合には、病的な心臓は、標準的な欠陥性の単極性信号又は双極性信号を生成し、この信号内の欠陥は、対応する心臓欠陥によって引き起こされる。
図5は、本発明の一実施形態による、心内信号を解析する際の、プロセッサ40によって実行される工程を示す流れ図200である。以下の説明では、信号は、双極性信号を含むものと想定される。当業者であれば、この説明を、必要な変更を加えて、単極性信号に関して適応させることが可能であろう。
最初の工程202では、専門家28は、心臓34内にプローブ24を挿入することにより、電極22A及び電極22Bを、心臓壁の一区域と接触させる。プロセッサ40が、電極から心内双極性ECG信号を取得し、各ECG信号は、電位V及び時間tの順序対:{(V,t)}を含む。
心拍選択工程204では、1つの完全な心拍を選択する。それゆえ、選択された心拍の持続時間が、Tであり、工程202での取得が、サンプリング周波数SampleRateで実行される場合には、選択された心拍内には、双極性信号の約T/SampleRate個のサンプルが存在する。
解析工程206では、プロセッサは、選択された心拍に、等式(2)を適合させて、選択された心拍に対する最良適合を与える、等式(2)の適合パラメーターの値のセットを導く。
比較工程208では、プロセッサは、経路指定モジュール30を使用して、電極22A及び電極22Bが、心臓に対して同じ位置にあるか否かを確認する。この比較が、肯定的な結果を返すことにより、電極が同じ位置にある場合には、平均化工程210で、プロセッサは、その位置での全ての心拍に関する適合パラメーターを平均化して、平均化適合パラメーターのセットを生成する。流れ図は、次いで注釈工程212に進む。
比較が、否定的な結果を返すことにより、電極が移動している場合には、平均化は実行されずに、流れ図は、工程212に直接進む。
注釈時間工程212では、工程210(平均化が実施される場合)又は工程206(平均化が存在していない場合)のいずれかで導かれた適合パラメーターを使用して、注釈時間を推定する。注釈時間は、ECG信号の特性の発生の基準時間である。注釈時間は、身体表面ECGに対して、又は、例えば冠状静脈洞内に定置されるカテーテルからの、心内基準ECGに対して、規定することができる。基準注釈時間を規定するために使用される、典型的な信号特性としては、QRS群のRピーク最大値が生じる時間、QRS群の最小微分係数が生じる時間、完全な信号のエネルギーの中心が生じる時間、又は完全な信号の最初の兆候が生じる時間が挙げられるが、これらに限定されない。基準注釈時間は典型的には、信号が測定される心臓内の位置に応じて決定される。基準注釈時間に関する規定、及びそれらの注釈時間の値は、基準ECGサブモジュール37内に格納される。
マップ構築工程214では、プロセッサは、心臓34の電気解剖学的マップの点を構成する。マップの点を構成するために、プロセッサは、工程212で推定された注釈時間の差異、及び関連する基準注釈時間(サブモジュール37内に格納)を、心臓のマップ内に組み込む(経路指定モジュール30を使用して)(図1)。サブモジュール41もまた、この工程で使用される。
継続条件216が、肯定的な結果を返すことによって、工程202〜214の繰り返しが指示される。典型的には、専門家28が、工程214のマッピング手順の停止を決定することによって、条件216が否定的な結果を返す場合には、流れ図は終了する。
上述のように、工程202〜214は、典型的には、健常な心臓内、及び既知の欠陥を有する病的な心臓内での、電極の種々の位置を含む、種々の状況に関して、実行することができる。
図6は、本発明の一実施形態による、上述の方法の適用の結果を示す、概略的グラフを示す。心内ECG信号を、様々な異なる場合から記録して、データプールを作成した。約5,900の心拍を、データプールから抽出した。全ての心拍を11のグループへと編成し、各グループは、規定の閾値未満の振幅を有する心拍を含むものとした。
この閾値は、信号の雑音の尺度であるため、より低い閾値を有する信号は、より高い雑音レベルを有する。特定のグループ内の各心拍に関して、Rピーク最大値の発生時間tRk、及び活性化波の通過の発生時間tCkを推定した。kは、測定されている心拍の数を表す指標である。tCkは、流れ図200に関して説明されるものと同様の、本明細書では適合注釈方法とも称される、適合解析を使用して、推定した。tRkを推定するための方法はまた、本明細書では、最大値注釈方法とも称される。
各グループの範囲内で、以下の時間の差異を算出した。
Figure 0006253879
等式(3)から、以下の変動係数を算出した。
Figure 0006253879
式中、σ(Δt)は、全てのΔt値の標準偏差であり、
M(Δt)は、全てのΔt値の平均である。
等式(4)の式は、所定のグループの範囲内の心拍の、最大値注釈方法又は適合注釈方法による、注釈時間の変動性の尺度を与える。
グラフ300は、1つのグループの、変動性VAR対閾値をプロットし、グラフ302は、グラフ300の直線回帰である。グラフ310は、1つのグループの、変動性VAR対閾値をプロットし、グラフ312は、グラフ310の直線回帰である。2つのグラフのセットの比較によって、低い閾値の値に関しては、すなわち、高い雑音の値を有する信号に関しては、本明細書で説明される方法に従って処理された、すなわち、適合注釈方法を使用して処理された信号の変動性は、これらの方法で処理されていない信号の変動性よりも小さいことが、明らかである。
上述の実施形態は、例示として記載されるものであり、本発明は、上記で具体的に図示及び説明されているものに限定されないことが理解されるであろう。むしろ、本発明の範囲には、上記で説明された様々な特徴の組み合わせ及び部分的組み合わせ、並びに当業者であれば前述の説明を読むことで想到されるであろう、先行技術に開示されていない、それらの変型及び修正が含まれる。
〔実施の態様〕
(1) 信号を解析するための方法であって、
時変心内電位信号を感知することと、
既定の振動波形に対する、前記時変心内電位信号の適合を見出すことと、
前記適合に応答して、前記信号の注釈時間を推定することと、を含む、方法。
(2) 前記時変心内電位信号が、単極性信号を含む、実施態様1に記載の方法。
(3) 前記既定の振動波形が、単一の極大値、単一の極小値、及び前記極小値と前記極大値とを分離する単一の変曲を有する、単一の完全な振動を含む、実施態様2に記載の方法。
(4) 前記既定の振動波形が、ガウス関数の一次微分を含む、実施態様1に記載の方法。
(5) 前記一次微分が、非対称因子によって歪められる、実施態様4に記載の方法。
(6) 前記時変心内電位信号が、双極性信号を含む、実施態様1に記載の方法。
(7) 前記既定の振動波形が、第1の単一の完全な振動と、第2の単一の完全な振動との差異を含む、実施態様6に記載の方法。
(8) 前記第1の単一の完全な振動が、第1の単一の極大値、第1の単一の極小値、及び前記第1の極大値と前記第1の極小値とを分離する第1の変曲を含み、前記第2の単一の完全な振動が、第2の単一の極大値、第2の単一の極小値、及び前記第2の極大値と前記第2の極小値とを分離する第2の変曲を含む、実施態様7に記載の方法。
(9) 前記第1の単一の完全な振動と、前記第2の単一の完全な振動とが、時間差によって分離される、実施態様7に記載の方法。
(10) 前記時間差が、前記双極性信号を生成する電極の、空間的分離の関数である、実施態様9に記載の方法。
(11) 前記時間差が、活性化波の伝播方向に対する、電極の配向の関数である、実施態様9に記載の方法。
(12) 前記既定の振動波形が、第1のガウス関数の一次微分と、第2のガウス関数の一次微分との差異を含む、実施態様6に記載の方法。
(13) 前記第1のガウス関数の一次微分が、第1の非対称因子によって歪められ、前記第2のガウス関数の一次微分が、第2の非対称因子によって歪められる、実施態様12に記載の方法。
(14) 前記時変心内電位信号が、時間差を間に有する3つ又は4つ以上の単極性信号を含み、活性化波の伝播方向が、前記時間差の関数である、実施態様1に記載の方法。
(15) それぞれの位置を有するそれぞれの電極が、前記3つ又は4つ以上の単極性信号を生成し、前記それぞれの位置が、前記関数のパラメーターを含む、実施態様14に記載の方法。
(16) 信号を解析するための装置であって、
時変心内電位信号を感知するように構成されるセンサーと、
既定の振動波形に対する、前記時変心内電位信号の適合を見出し、前記適合に応答して、前記信号の注釈時間を推定するように構成される、プロセッサと、
を含む、装置。
(17) 前記時変心内電位信号が、単極性信号を含む、実施態様16に記載の装置。
(18) 前記既定の振動波形が、単一の極大値、単一の極小値、及び前記極小値と前記極大値とを分離する単一の変曲を有する、単一の完全な振動を含む、実施態様17に記載の装置。
(19) 前記既定の振動波形が、ガウス関数の一次微分を含む、実施態様16に記載の装置。
(20) 前記一次微分が、非対称因子によって歪められる、実施態様19に記載の装置。
(21) 前記時変心内電位信号が、双極性信号を含む、実施態様16に記載の装置。
(22) 前記既定の振動波形が、第1の単一の完全な振動と、第2の単一の完全な振動との差異を含む、実施態様21に記載の装置。
(23) 前記第1の単一の完全な振動が、第1の単一の極大値、第1の単一の極小値、及び前記第1の極大値と前記第1の極小値とを分離する第1の変曲を含み、前記第2の単一の完全な振動が、第2の単一の極大値、第2の単一の極小値、及び前記第2の極大値と前記第2の極小値とを分離する第2の変曲を含む、実施態様22に記載の装置。
(24) 前記第1の単一の完全な振動と、前記第2の単一の完全な振動とが、時間差によって分離される、実施態様22に記載の装置。
(25) 前記時間差が、前記双極性信号を生成する電極の、空間的分離の関数である、実施態様24に記載の装置。
(26) 前記時間差が、活性化波の伝播方向に対する、電極の配向の関数である、実施態様24に記載の装置。
(27) 前記既定の振動波形が、第1のガウス関数の一次微分と、第2のガウス関数の一次微分との差異を含む、実施態様21に記載の装置。
(28) 前記第1のガウス関数の一次微分が、第1の非対称因子によって歪められ、前記第2のガウス関数の一次微分が、第2の非対称因子によって歪められる、実施態様27に記載の装置。
(29) 前記時変心内電位信号が、時間差を間に有する3つ又は4つ以上の単極性信号を含み、活性化波の伝播方向が、前記時間差の関数である、実施態様16に記載の装置。
(30) それぞれの位置を有するそれぞれの電極が、前記3つ又は4つ以上の単極性信号を生成し、前記それぞれの位置が、前記関数のパラメーターを含む、実施態様29に記載の装置。

Claims (13)

  1. 信号を解析するための装置であって、
    時変心内電位信号を感知するように構成されるセンサーと、
    既定の振動波形に対する、前記時変心内電位信号の適合を見出し、前記適合に応答して、前記信号の特性の発生の基準時間を推定するように構成される、プロセッサと、
    を含み、
    前記時変心内電位信号が、双極性信号を含み、
    前記既定の振動波形が、第1の単一の完全な振動と、第2の単一の完全な振動との差分を含む、装置。
  2. 前記時変心内電位信号が、単極性信号を含む、請求項に記載の装置。
  3. 前記既定の振動波形が、単一の極大値、単一の極小値、及び前記極小値と前記極大値とを分離する単一の変曲を有する、単一の完全な振動を含む、請求項に記載の装置。
  4. 前記既定の振動波形が、ガウス関数の一次微分を含む、請求項に記載の装置。
  5. 前記一次微分が、非対称因子によって歪められる、請求項に記載の装置。
  6. 前記第1の単一の完全な振動が、第1の単一の極大値、第1の単一の極小値、及び前記第1の極大値と前記第1の極小値とを分離する第1の変曲を含み、前記第2の単一の完全な振動が、第2の単一の極大値、第2の単一の極小値、及び前記第2の極大値と前記第2の極小値とを分離する第2の変曲を含む、請求項に記載の装置。
  7. 前記第1の単一の完全な振動と、前記第2の単一の完全な振動とが、時間差によって分離される、請求項に記載の装置。
  8. 前記時間差が、前記双極性信号を生成する電極の、空間的分離の関数である、請求項に記載の装置。
  9. 前記時間差が、活性化波の伝播方向に対する、電極の配向の関数である、請求項に記載の装置。
  10. 前記既定の振動波形が、第1のガウス関数の一次微分と、第2のガウス関数の一次微分との差異を含む、請求項に記載の装置。
  11. 前記第1のガウス関数の一次微分が、第1の非対称因子によって歪められ、前記第2のガウス関数の一次微分が、第2の非対称因子によって歪められる、請求項10に記載の装置。
  12. 前記時変心内電位信号が、時間差を間に有する3つ又は4つ以上の単極性信号を含み、活性化波の伝播方向が、前記時間差の関数である、請求項に記載の装置。
  13. それぞれの位置を有するそれぞれの電極が、前記3つ又は4つ以上の単極性信号を生成し、前記それぞれの位置が、前記関数のパラメーターを含む、請求項12に記載の装置。
JP2012247240A 2011-11-11 2012-11-09 心内心電計信号の正確な時間注釈 Active JP6253879B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US13/294,412 2011-11-11
US13/294,412 US8700136B2 (en) 2011-11-11 2011-11-11 Accurate time annotation of intracardiac ECG signals

Publications (2)

Publication Number Publication Date
JP2013103130A JP2013103130A (ja) 2013-05-30
JP6253879B2 true JP6253879B2 (ja) 2017-12-27

Family

ID=47226014

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012247240A Active JP6253879B2 (ja) 2011-11-11 2012-11-09 心内心電計信号の正確な時間注釈

Country Status (7)

Country Link
US (1) US8700136B2 (ja)
EP (1) EP2591721B1 (ja)
JP (1) JP6253879B2 (ja)
CN (1) CN103099614B (ja)
AU (1) AU2012247057B2 (ja)
CA (1) CA2794794C (ja)
IL (1) IL222660A (ja)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10758141B2 (en) 2012-04-23 2020-09-01 Biosense Webster (Israel) Ltd. Cardiac activation time detection
US9833157B2 (en) * 2012-04-23 2017-12-05 Biosense Webster (Israel) Ltd. Cardiac activation time detection
WO2015006847A1 (en) * 2013-07-19 2015-01-22 Litezone Technologies Inc. Pressure compensated glass unit
US9554718B2 (en) 2014-01-29 2017-01-31 Biosense Webster (Israel) Ltd. Double bipolar configuration for atrial fibrillation annotation
KR102169378B1 (ko) 2014-02-18 2020-10-23 삼성전자주식회사 ECG(electrocardiogram) 센서 및 이의 동작 방법
US9259165B2 (en) 2014-02-26 2016-02-16 Biosense Webster (Israel) Ltd. Determination of reference annotation time from multi-channel electro-cardiogram signals
US10506939B2 (en) * 2015-08-11 2019-12-17 Biosense Webster (Israel) Ltd. Matching and tracking time sequences of heart activation
US11058342B2 (en) * 2016-08-11 2021-07-13 Biosense Webster (Israel) Ltd. Classifying ECG signals
US11389116B2 (en) 2018-01-19 2022-07-19 Biosense Webster (Israel) Ltd. Apparatus and method for heartbeat classification based on time sequence and morphology of intracardiac and body surface electrocardiogram (ECG) signals
US20200178831A1 (en) 2018-12-11 2020-06-11 Biosense Webster (Israel) Ltd. Display of arrhythmia type
US11766206B2 (en) * 2018-12-20 2023-09-26 Biosense Webster (Israel) Ltd. Electrode arrangement to sense cardiac wave vector
US11116435B2 (en) 2019-08-26 2021-09-14 Biosense Webster (Israel) Ltd. Automatic identification of a location of focal source in atrial fibrillation (AF)
US20210059549A1 (en) 2019-08-26 2021-03-04 Biosense Webster (Israel) Ltd. Error estimation of local activation times (lat) measured by multiple electrode catheter
US11730414B2 (en) 2020-01-21 2023-08-22 Biosense Webster (Israel) Ltd. Automatic pattern acquisition
US12239449B2 (en) * 2020-06-10 2025-03-04 Biosense Webster (Israel) Ltd. System and method to detect stable arrhythmia heartbeat and to calculate and detect cardiac mapping annotations
US20220079491A1 (en) 2020-09-14 2022-03-17 Biosense Webster (Israel) Ltd. Local activation time analysis system
US12544129B2 (en) * 2021-06-07 2026-02-10 Biosense Webster (Israel) Ltd. Bipolar electrode pair selection
US12324619B2 (en) 2021-06-07 2025-06-10 Biosense Webster (Israel) Ltd. Automatic anatomical feature identification and map segmentation

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5273049A (en) * 1992-04-09 1993-12-28 Telectronics Pacing Systems, Inc. Detection of cardiac arrhythmias using template matching by signature analysis
US5400795A (en) * 1993-10-22 1995-03-28 Telectronics Pacing Systems, Inc. Method of classifying heart rhythms by analyzing several morphology defining metrics derived for a patient's QRS complex
JP2001218747A (ja) * 2000-02-10 2001-08-14 Tokyo Denki Univ 心電図解析システム、および、そのためのプログラムを記録したコンピュータ読み取り可能な記録媒体
US6829501B2 (en) * 2001-12-20 2004-12-07 Ge Medical Systems Information Technologies, Inc. Patient monitor and method with non-invasive cardiac output monitoring
US6884218B2 (en) * 2002-12-09 2005-04-26 Charles W. Olson Three dimensional vector cardiograph and method for detecting and monitoring ischemic events
SE0301201D0 (sv) 2003-04-24 2003-04-24 St Jude Medical Apparatus for analysing cardiac events
BRPI0414345A (pt) 2003-09-12 2006-11-07 Bodymedia Inc método e aparelho para medição de parámetros relacionados com o coração
US7412282B2 (en) 2005-01-26 2008-08-12 Medtronic, Inc. Algorithms for detecting cardiac arrhythmia and methods and apparatuses utilizing the algorithms
CN1951320A (zh) * 2006-11-16 2007-04-25 上海交通大学 基于超完备特征的异常心电识别方法
CN101243972A (zh) * 2008-03-10 2008-08-20 湖南大学 用于心电信号qrs波检测的对数域模拟小波变换电路
CN101972143A (zh) * 2010-11-16 2011-02-16 电子科技大学 基于盲源提取的房颤监测方法

Also Published As

Publication number Publication date
IL222660A (en) 2016-04-21
AU2012247057A1 (en) 2013-05-30
CA2794794A1 (en) 2013-05-11
EP2591721B1 (en) 2018-03-21
CN103099614B (zh) 2016-08-03
AU2012247057B2 (en) 2014-10-30
US20130123652A1 (en) 2013-05-16
JP2013103130A (ja) 2013-05-30
EP2591721A1 (en) 2013-05-15
CN103099614A (zh) 2013-05-15
US8700136B2 (en) 2014-04-15
CA2794794C (en) 2018-07-31

Similar Documents

Publication Publication Date Title
JP6253879B2 (ja) 心内心電計信号の正確な時間注釈
JP7604180B2 (ja) 複数のグループ化されたecg信号の統計的特性を使用した一貫性のない信号の検出
JP7625400B2 (ja) 不一致を検出することによる、ecg信号のマッピングの遡及的最適化
CN104799844B (zh) 激活波阵面的混合的双极性/单极性检测
JP6559439B2 (ja) 複数チャネル心電図信号からの基準注釈時間の決定
EP2663228B1 (en) Electroanatomical mapping
JP6117433B2 (ja) 類似性ベースのパターンマッチングによる増強された活動開始時間の最適化のための解剖学的マッピングシステム
CN105708442B (zh) 心室远场衰减
JP2017514597A (ja) 心臓組織をマッピングする医療用デバイス
JP2017520286A (ja) 心臓組織をマッピングする医療用デバイス
CN108685613B (zh) 组织传导速度
CN113116359A (zh) 用于估计残余ecg噪声水平和自适应噪声阈值的方法和系统
JP2022099313A (ja) 局所異常心室興奮(lava)信号を含む遅延電位のアノテーション
US20250169765A1 (en) Respiration compensation for mapping
EP3883466A1 (en) Graph total variation for ecgi
AU2014250609B2 (en) Accurate time annotation of intracardiac ecg signals

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20150827

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20160720

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160823

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20161122

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170509

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170724

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20171121

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20171129

R150 Certificate of patent or registration of utility model

Ref document number: 6253879

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250