JP3213602B2 - 心内電位解析装置及び方法 - Google Patents
心内電位解析装置及び方法Info
- Publication number
- JP3213602B2 JP3213602B2 JP21194699A JP21194699A JP3213602B2 JP 3213602 B2 JP3213602 B2 JP 3213602B2 JP 21194699 A JP21194699 A JP 21194699A JP 21194699 A JP21194699 A JP 21194699A JP 3213602 B2 JP3213602 B2 JP 3213602B2
- Authority
- JP
- Japan
- Prior art keywords
- electrodes
- vector
- excitation wave
- selecting
- excitation
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 22
- 230000005284 excitation Effects 0.000 claims description 103
- 239000013598 vector Substances 0.000 claims description 66
- 238000004458 analytical method Methods 0.000 claims description 54
- 206010003119 arrhythmia Diseases 0.000 claims description 33
- 230000006793 arrhythmia Effects 0.000 claims description 33
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000013500 data storage Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 16
- 230000008569 process Effects 0.000 description 9
- 238000012937 correction Methods 0.000 description 8
- 230000002107 myocardial effect Effects 0.000 description 8
- 238000005259 measurement Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 6
- 210000004165 myocardium Anatomy 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000000737 periodic effect Effects 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 230000001960 triggered effect Effects 0.000 description 3
- 206010003658 Atrial Fibrillation Diseases 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000002964 excitative effect Effects 0.000 description 2
- 210000002837 heart atrium Anatomy 0.000 description 2
- 208000003663 ventricular fibrillation Diseases 0.000 description 2
- 206010001497 Agitation Diseases 0.000 description 1
- 206010003662 Atrial flutter Diseases 0.000 description 1
- 238000002679 ablation Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 206010003668 atrial tachycardia Diseases 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000004070 electrodeposition Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 230000036279 refractory period Effects 0.000 description 1
- 210000005245 right atrium Anatomy 0.000 description 1
- 238000009987 spinning Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 230000002459 sustained effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
- 206010047302 ventricular tachycardia Diseases 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6846—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be brought in contact with an internal body part, i.e. invasive
- A61B5/6847—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be brought in contact with an internal body part, i.e. invasive mounted on an invasive device
- A61B5/6852—Catheters
- A61B5/6858—Catheters with a distal basket, e.g. expandable basket
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/25—Bioelectric electrodes therefor
- A61B5/279—Bioelectric electrodes therefor specially adapted for particular uses
- A61B5/28—Bioelectric electrodes therefor specially adapted for particular uses for electrocardiography [ECG]
- A61B5/283—Invasive
- A61B5/287—Holders for multiple electrodes, e.g. electrode catheters for electrophysiological study [EPS]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/339—Displays specially adapted therefor
- A61B5/341—Vectorcardiography [VCG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/04—Arrangements of multiple sensors of the same type
- A61B2562/043—Arrangements of multiple sensors of the same type in a linear array
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Cardiology (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Physiology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Description
【０００１】
【発明の属する技術分野】本発明は，コンピュータを用
いた心臓電気生理検査における心内電位解析装置及び方
法に関し，特にリエントリータイプの不整脈（心房頻
拍，心房粗動，心房細動，心室頻拍，心室細動など）の
興奮伝播経路を解析する装置及び方法に関する。
いた心臓電気生理検査における心内電位解析装置及び方
法に関し，特にリエントリータイプの不整脈（心房頻
拍，心房粗動，心房細動，心室頻拍，心室細動など）の
興奮伝播経路を解析する装置及び方法に関する。
【０００２】
【従来の技術】不整脈は，そのメカニズムにより大きく
３種類に分類されている。それらは，異常自動能，リエ
ントリー，トリガードアクティビティーである。
３種類に分類されている。それらは，異常自動能，リエ
ントリー，トリガードアクティビティーである。
【０００３】異常自動能，トリガードアクティビティー
による不整脈は，ある場所から生じた興奮波が，心内を
広がって，境界に到達したり，違う経路で進んだ異なる
興奮波とぶつかったりして，最初の場所に戻ってくるこ
となく，消失してしまうような不整脈であると理解する
ことができる。逆に，リエントリーによる不整脈は，あ
る場所を通った興奮波が，心内を伝わっていって，また
同じ場所に戻ってくることによって持続する不整脈であ
ると理解することができる。
による不整脈は，ある場所から生じた興奮波が，心内を
広がって，境界に到達したり，違う経路で進んだ異なる
興奮波とぶつかったりして，最初の場所に戻ってくるこ
となく，消失してしまうような不整脈であると理解する
ことができる。逆に，リエントリーによる不整脈は，あ
る場所を通った興奮波が，心内を伝わっていって，また
同じ場所に戻ってくることによって持続する不整脈であ
ると理解することができる。
【０００４】現在，心内興奮伝播の解析には，基本とし
てある時刻を開始時刻とし，この時刻からある時間経過
した時刻に興奮している場所を線でつないで，興奮がど
のように伝わったかを表現する等時線図と呼ばれる方法
が用いられている。これは，地図で用いられる等高線と
まったく同じものであり，山の山頂が不整脈の開始で，
麓に向かって興奮が伝わって行くようなイメージで理解
されるような不整脈，つまり異常自動能やトリガードア
クティビティーをメカニズムとする不整脈には有用であ
るといえる。いいかえれば，ある時刻から興奮が開始し
て，ある時刻に終わるパターンの不整脈には等時線図は
有用である。
てある時刻を開始時刻とし，この時刻からある時間経過
した時刻に興奮している場所を線でつないで，興奮がど
のように伝わったかを表現する等時線図と呼ばれる方法
が用いられている。これは，地図で用いられる等高線と
まったく同じものであり，山の山頂が不整脈の開始で，
麓に向かって興奮が伝わって行くようなイメージで理解
されるような不整脈，つまり異常自動能やトリガードア
クティビティーをメカニズムとする不整脈には有用であ
るといえる。いいかえれば，ある時刻から興奮が開始し
て，ある時刻に終わるパターンの不整脈には等時線図は
有用である。
【０００５】
【発明が解決しようとする課題】しかし，リエントリー
による不整脈では，興奮が回路をぐるぐる回ることで成
立するため，興奮の開始も終了も特定できない。ある時
点より前には，必ずどこかに興奮があって，前に前に追
って行くと，元いた所に戻ってきてしまうのである。無
理に開始時刻，終了時刻を設定して等時線図解析を行う
と，設定した開始時刻，終了時刻によって解析結果が異
なってしまうという不具合が生じる。図１５Ａ，１５Ｂ
にその例を示す。図１５Ａ，１５Ｂは図２の心電図に興
奮開始時刻として異なるＴ１，Ｔ２を設定して作図した
例である。この図から明らかなように，開始時刻が異な
ると結果も異なる。
による不整脈では，興奮が回路をぐるぐる回ることで成
立するため，興奮の開始も終了も特定できない。ある時
点より前には，必ずどこかに興奮があって，前に前に追
って行くと，元いた所に戻ってきてしまうのである。無
理に開始時刻，終了時刻を設定して等時線図解析を行う
と，設定した開始時刻，終了時刻によって解析結果が異
なってしまうという不具合が生じる。図１５Ａ，１５Ｂ
にその例を示す。図１５Ａ，１５Ｂは図２の心電図に興
奮開始時刻として異なるＴ１，Ｔ２を設定して作図した
例である。この図から明らかなように，開始時刻が異な
ると結果も異なる。
【０００６】従って，リエントリーによる不整脈の心内
興奮伝播経路を容易に推定することのできる解析技術の
確立が待ち望まれていた。
興奮伝播経路を容易に推定することのできる解析技術の
確立が待ち望まれていた。
【０００７】
【課題を解決するための手段】心内興奮伝播は，ある場
所の興奮が隣の場所の興奮を引き起こすため，ちょうど
波が伝わるようなものと理解することができる。水面を
伝わる波と異なるのは，一度波が通った場所はある時間
波が通ることができないということである。（心筋が興
奮したあと，不応期と呼ばれる時間があり，この時間は
興奮することができないからである。）したがって，反
射や，重なりという現象は見られない。つまり，この波
は心筋の境界，たとえば心房と静脈の境界や心房と心室
の境界に到達すると消失する。また，興奮波同士がぶっ
かるとやはり消失してしまう。興奮伝播をこのような特
徴を持った興奮波として扱うことができる。
所の興奮が隣の場所の興奮を引き起こすため，ちょうど
波が伝わるようなものと理解することができる。水面を
伝わる波と異なるのは，一度波が通った場所はある時間
波が通ることができないということである。（心筋が興
奮したあと，不応期と呼ばれる時間があり，この時間は
興奮することができないからである。）したがって，反
射や，重なりという現象は見られない。つまり，この波
は心筋の境界，たとえば心房と静脈の境界や心房と心室
の境界に到達すると消失する。また，興奮波同士がぶっ
かるとやはり消失してしまう。興奮伝播をこのような特
徴を持った興奮波として扱うことができる。
【０００８】本発明では，十分に狭い範囲では，興奮波
は直進する，と仮定して解析することで，興奮波の進行
の性質を知るようにした。直進する波の進行速度と進行
方向は，場所のわかっている異なる３つの電極を興奮波
が通過する時間を測定し，計算することで求めることが
できる。従って，本発明は，ある不整脈にとって興奮波
の進行方向が変化しないと仮定できる十分狭い範囲内に
ある電極の３点を興奮波が通過する時刻と心電図上に任
意に設定した参照時刻との時間差（以下，興奮波通過時
刻と呼ぶ）と，その3つの双極電極の位置情報に基づい
て，興奮波前面の進行方向と速度とを計算し，ベクトル
として表示することで心内興奮伝播を表現し，リエント
リーによる不整脈の心内興奮伝播経路の推定を容易にし
た。
は直進する，と仮定して解析することで，興奮波の進行
の性質を知るようにした。直進する波の進行速度と進行
方向は，場所のわかっている異なる３つの電極を興奮波
が通過する時間を測定し，計算することで求めることが
できる。従って，本発明は，ある不整脈にとって興奮波
の進行方向が変化しないと仮定できる十分狭い範囲内に
ある電極の３点を興奮波が通過する時刻と心電図上に任
意に設定した参照時刻との時間差（以下，興奮波通過時
刻と呼ぶ）と，その3つの双極電極の位置情報に基づい
て，興奮波前面の進行方向と速度とを計算し，ベクトル
として表示することで心内興奮伝播を表現し，リエント
リーによる不整脈の心内興奮伝播経路の推定を容易にし
た。
【０００９】
【発明の実施の形態】本発明による心内電位解析は，心
筋興奮の測定データを解析するものである。心筋興奮の
測定はさまざまな方法があるが，例えば複数の電極を収
納した細いカテーテル電極を血管を通して心臓内部に挿
入し，心臓内で複数の電極を広げ，多点で心筋の興奮状
態を測定して記録することによって行われる。この場
合，心筋の興奮時刻と測定電極位置が対応付けられて記
録される。
筋興奮の測定データを解析するものである。心筋興奮の
測定はさまざまな方法があるが，例えば複数の電極を収
納した細いカテーテル電極を血管を通して心臓内部に挿
入し，心臓内で複数の電極を広げ，多点で心筋の興奮状
態を測定して記録することによって行われる。この場
合，心筋の興奮時刻と測定電極位置が対応付けられて記
録される。
【００１０】図１Ａ，１Ｂに６４極のバスケットカテー
テルを用いた例を示す。バスケットカテーテルには，Ａ
からＨの８本の弧状の電線にそれぞれ８極の電極が装着
されており，計６４極の電極が球状に位置し，相対的な
位置が固定されている。図１Ａは心臓内で球状に広げら
れた電極の様子を示す。それぞれの電極は電線記号（Ａ
−Ｈ）と電極番号（１，２，３，…，８）を組にしてＡ
１，Ａ２，．．．，Ｈ７，Ｈ８と名付けられている。
（図１Ｂ）
テルを用いた例を示す。バスケットカテーテルには，Ａ
からＨの８本の弧状の電線にそれぞれ８極の電極が装着
されており，計６４極の電極が球状に位置し，相対的な
位置が固定されている。図１Ａは心臓内で球状に広げら
れた電極の様子を示す。それぞれの電極は電線記号（Ａ
−Ｈ）と電極番号（１，２，３，…，８）を組にしてＡ
１，Ａ２，．．．，Ｈ７，Ｈ８と名付けられている。
（図１Ｂ）
【００１１】ここで，電極Ａ１とＡ２による双極電極に
Ａ１２と名前を付け，Ａ１２，Ａ３４，Ａ５
６，．．．，Ｈ５６，Ｈ７８の計３２組の電位差を双極
電極の位置に対応させて記録する。双極電極の位置は２
つの電極のどちらかの位置あるいはその中間点など，２
つの電極の近傍に，該２つの電極位置に関連付けて設定
された位置である。
Ａ１２と名前を付け，Ａ１２，Ａ３４，Ａ５
６，．．．，Ｈ５６，Ｈ７８の計３２組の電位差を双極
電極の位置に対応させて記録する。双極電極の位置は２
つの電極のどちらかの位置あるいはその中間点など，２
つの電極の近傍に，該２つの電極位置に関連付けて設定
された位置である。
【００１２】図２はポリグラフの出力であり，３２組の
双極電極によるリエントリーによる不整脈の記録例を示
す。図２の各チャンネル（Ａ１２，Ａ３４，…，Ｈ７
８）は心腔内心電図を示している。本発明では一回又は
複数回の拍動による心筋興奮を複数の測定点で測定し，
解析する。
双極電極によるリエントリーによる不整脈の記録例を示
す。図２の各チャンネル（Ａ１２，Ａ３４，…，Ｈ７
８）は心腔内心電図を示している。本発明では一回又は
複数回の拍動による心筋興奮を複数の測定点で測定し，
解析する。
【００１３】本実施例では，解析結果を２次元平面上に
ベクトルとして表すため，３次元空間に位置を占める複
数の双極電極のそれぞれの位置を二次元平面上に変換し
ている。図３に二次元平面上に変換された双極電極の相
対的位置関係を示す座標マップの一例を示す。それぞれ
の測定の位置は，ある場所を基準とすれば，一意に決定
される。例では，各双極電極の相対的位置関係に重点を
置き，各双極電極がＸ，Ｙ軸で定義される平面上にある
ものと仮定して，Ａ１２のＸ座標を０，Ｙ座標を３，と
して，その他の双極電極の座標を図３のように決定して
いる。最終的な結果もこの位置関係の上に描かれる。図
３の双極電極の座標マップはプログラムにより前もって
メモリ上に登録される。
ベクトルとして表すため，３次元空間に位置を占める複
数の双極電極のそれぞれの位置を二次元平面上に変換し
ている。図３に二次元平面上に変換された双極電極の相
対的位置関係を示す座標マップの一例を示す。それぞれ
の測定の位置は，ある場所を基準とすれば，一意に決定
される。例では，各双極電極の相対的位置関係に重点を
置き，各双極電極がＸ，Ｙ軸で定義される平面上にある
ものと仮定して，Ａ１２のＸ座標を０，Ｙ座標を３，と
して，その他の双極電極の座標を図３のように決定して
いる。最終的な結果もこの位置関係の上に描かれる。図
３の双極電極の座標マップはプログラムにより前もって
メモリ上に登録される。
【００１４】図４に本発明の解析装置により実行される
解析処理のフローチャートを示す。 ステップ１：図３に一例として示されたようにな複数の
双極電極の座標マップをメモリに登録する。
解析処理のフローチャートを示す。 ステップ１：図３に一例として示されたようにな複数の
双極電極の座標マップをメモリに登録する。
【００１５】ステップ２： 図２の心電図の時間軸上に
おいて，任意に設定した参照時刻Ｔ１（またはＴ２）か
ら，興奮の開始までの時間を興奮波通過時刻として測定
し，本発明の解析装置に入力する。
おいて，任意に設定した参照時刻Ｔ１（またはＴ２）か
ら，興奮の開始までの時間を興奮波通過時刻として測定
し，本発明の解析装置に入力する。
【００１６】参照時刻はどこに設定してもかまわない。
本発明では参照時刻に関係なく全く同じ結果が得られ
る。
本発明では参照時刻に関係なく全く同じ結果が得られ
る。
【００１７】入力された興奮波通過時刻は双極電極に対
応してメモリに記憶される。興奮波通過時刻の測定はユ
ーザーが心電図上で実測してもよいし，機械で自動的に
読み取り，入力するようにしてもよい。
応してメモリに記憶される。興奮波通過時刻の測定はユ
ーザーが心電図上で実測してもよいし，機械で自動的に
読み取り，入力するようにしてもよい。
【００１８】ステップ３：双極電極の中から，隣り合う
双極電極３点を選択する。本実施例では，初めから場所
が確定されているため，３点は，（Ａ１２，Ａ３４，Ｂ
１２），（Ａ１２，Ａ３４，Ｂ３４），．．．という９
６の組合せがあらかじめ決定されている。
双極電極３点を選択する。本実施例では，初めから場所
が確定されているため，３点は，（Ａ１２，Ａ３４，Ｂ
１２），（Ａ１２，Ａ３４，Ｂ３４），．．．という９
６の組合せがあらかじめ決定されている。
【００１９】ステップ４：選択された３つの双極電極に
関し，ステップ２で入力された興奮波通過時刻から，３
点を通過した興奮波ベクトルを下記の式により求め，記
憶する。
関し，ステップ２で入力された興奮波通過時刻から，３
点を通過した興奮波ベクトルを下記の式により求め，記
憶する。
【００２０】３つの双極電極の位置を，互いに直交する
ｘ，ｙ，ｚ軸で定義された直交３次元座標系にＡ（Ａ
ｘ，Ａｙ，Ａｚ），Ｂ（Ｂｘ，Ｂｙ，Ｂｚ），Ｃ（Ｃ
ｘ，Ｃｙ，Ｃｚ）と設定する。３点Ａ，Ｂ，Ｃに対応し
た興奮波通過時刻をそれぞれＡｔ，Ｂｔ，Ｃｔとする
と，興奮波ベクトル
ｘ，ｙ，ｚ軸で定義された直交３次元座標系にＡ（Ａ
ｘ，Ａｙ，Ａｚ），Ｂ（Ｂｘ，Ｂｙ，Ｂｚ），Ｃ（Ｃ
ｘ，Ｃｙ，Ｃｚ）と設定する。３点Ａ，Ｂ，Ｃに対応し
た興奮波通過時刻をそれぞれＡｔ，Ｂｔ，Ｃｔとする
と，興奮波ベクトル
【００２１】
【数１０】 は下記の式で求められる。
【００２２】
【数１１】 但し，
【００２３】
【数１２】 である。２次元での計算も式（１）はそのまま適用でき
る。この式から３点を通過した興奮波ベクトルを求め，
興奮波，即ちベクトルの進行方向，速度を求める。
る。この式から３点を通過した興奮波ベクトルを求め，
興奮波，即ちベクトルの進行方向，速度を求める。
【００２４】ステップ５：選択された３点で興奮波前面
の速度が変化していないと仮定できるかどうかの判定を
行う。この場合の判定基準として上限，下限の閾値とし
て所定の進行速度を設定できる。速度が変化していない
と仮定できると判断した場合，ステップ６に進む。速度
が変化していると判断された場合や電位の記録できなか
った点があれば，ステップ３に戻り，双極電極の次の組
合せを選択する。
の速度が変化していないと仮定できるかどうかの判定を
行う。この場合の判定基準として上限，下限の閾値とし
て所定の進行速度を設定できる。速度が変化していない
と仮定できると判断した場合，ステップ６に進む。速度
が変化していると判断された場合や電位の記録できなか
った点があれば，ステップ３に戻り，双極電極の次の組
合せを選択する。
【００２５】ステップ６：選択された双極電極３点に関
連して設定された所定の点をベクトルの始点として，計
算された進行方向にベクトルをステップ１で登録された
座標マップに関連付けて表示装置上に表示する。この
時，ベクトルの大きさは，進行速度の相対値となるが，
適切な表示になるように単位長を任意に設定する。ステ
ップ５で所定の範囲外と判断されたものであっても，範
囲外であることを表示するようにしてもよい。
連して設定された所定の点をベクトルの始点として，計
算された進行方向にベクトルをステップ１で登録された
座標マップに関連付けて表示装置上に表示する。この
時，ベクトルの大きさは，進行速度の相対値となるが，
適切な表示になるように単位長を任意に設定する。ステ
ップ５で所定の範囲外と判断されたものであっても，範
囲外であることを表示するようにしてもよい。
【００２６】ベクトルの始点としては例えば選択された
３点の中心点を選んでもよいし，３点の中で興奮波通過
時刻の最も短い点を選んでもよい。３点の中心点を始点
とする場合は，中心点Ｐ（Ｐｘ，Ｐｙ，Ｐｚ）を下記の
式によって計算する。
３点の中心点を選んでもよいし，３点の中で興奮波通過
時刻の最も短い点を選んでもよい。３点の中心点を始点
とする場合は，中心点Ｐ（Ｐｘ，Ｐｙ，Ｐｚ）を下記の
式によって計算する。
【００２７】
【数１３】
【００２８】図５に興奮波ベクトルを矢印として描く場
合の一例を示す。３点Ａ，Ｂ，Ｃとして，図３のＡ７８
（０，０），Ｂ５６（１，１），Ａ５６（０，１）を選
び，Ａｔ＝１０ｍｓ，Ｂｔ＝３０ｍｓ，Ｃｔ＝２０ｍ
ｓ，ベクトルの始点をＡ点とした場合の例を示す。
合の一例を示す。３点Ａ，Ｂ，Ｃとして，図３のＡ７８
（０，０），Ｂ５６（１，１），Ａ５６（０，１）を選
び，Ａｔ＝１０ｍｓ，Ｂｔ＝３０ｍｓ，Ｃｔ＝２０ｍ
ｓ，ベクトルの始点をＡ点とした場合の例を示す。
【００２９】
【数１４】 を式（１）に代入すると，
【００３０】
【数１５】 同様にνy＝０．０２５ よって
【００３１】
【数１６】 となり，図示すると図５のようになる。
【００３２】ステップ７：双極電極３点の組合せすべて
に対してステップ３−６の処理を行ったかどうか調べ，
ＮＯならステップ３に戻り，ＹＥＳなら解析処理を終了
する。
に対してステップ３−６の処理を行ったかどうか調べ，
ＮＯならステップ３に戻り，ＹＥＳなら解析処理を終了
する。
【００３３】図６に本発明による最終表示例を示す。興
奮波ベクトルは図３の双極電極の相対的位置関係を示し
た座標マップに表示されている。矢印の短い範囲は，興
奮伝播速度が小さいことを示す。矢印の方向から，興奮
波の進行方向も明らかである。この実施例の場合，Ｅ３
４のあたりから波が広がる様子が見られる。
奮波ベクトルは図３の双極電極の相対的位置関係を示し
た座標マップに表示されている。矢印の短い範囲は，興
奮伝播速度が小さいことを示す。矢印の方向から，興奮
波の進行方向も明らかである。この実施例の場合，Ｅ３
４のあたりから波が広がる様子が見られる。
【００３４】図６より想定される興奮伝播の経路を図７
に模式図として表現した。図７から，右房内に手術時の
切開による癖痕と思われる，興奮波が通ることのできな
い部分（二重線で示した）があり，その中で一部遅い速
度で通過できる部分（波線で示した）を通って，興奮波
が瘤痕の周りを旋回するというメカニズムが想定され
る。ここまで経路が明らかとなれば，経路のもっとも狭
い部分であるところの切開癖痕のすきまに，電極カテー
テルで高周波エネルギーを与えて，興奮波が通過できな
いように心筋を変性させる治療（カテーテル心筋焼灼術
と呼ばれる。）を行えば，この不整脈を治療することが
できることがわかる。
に模式図として表現した。図７から，右房内に手術時の
切開による癖痕と思われる，興奮波が通ることのできな
い部分（二重線で示した）があり，その中で一部遅い速
度で通過できる部分（波線で示した）を通って，興奮波
が瘤痕の周りを旋回するというメカニズムが想定され
る。ここまで経路が明らかとなれば，経路のもっとも狭
い部分であるところの切開癖痕のすきまに，電極カテー
テルで高周波エネルギーを与えて，興奮波が通過できな
いように心筋を変性させる治療（カテーテル心筋焼灼術
と呼ばれる。）を行えば，この不整脈を治療することが
できることがわかる。
【００３５】図８Ａ，８Ｂ，８Ｃ及び図９は一定の周期
を持つ不整脈を解析する場合の原理を説明するための図
である。電極Ａ１２，Ａ３４，Ｂ１２で測定された興奮
点をＥ１，Ｅ２，Ｅ３，興奮波の周期をＣＬ（ミリ
秒），１周期後の興奮点を「′」をつけて表示する。
を持つ不整脈を解析する場合の原理を説明するための図
である。電極Ａ１２，Ａ３４，Ｂ１２で測定された興奮
点をＥ１，Ｅ２，Ｅ３，興奮波の周期をＣＬ（ミリ
秒），１周期後の興奮点を「′」をつけて表示する。
【００３６】一定の周期を持つ不整脈を少なくとも一周
期分記録したポリグラフにおいて，任意の位置に参照時
刻を設定すると，測定された興奮波通過時刻の小さいほ
うから順に興奮が伝わるとは限らない。従って，３点の
興奮波通過時刻に補正を施す必要がある。
期分記録したポリグラフにおいて，任意の位置に参照時
刻を設定すると，測定された興奮波通過時刻の小さいほ
うから順に興奮が伝わるとは限らない。従って，３点の
興奮波通過時刻に補正を施す必要がある。
【００３７】図８Ａの場合，Ｅ１点とＥ２点の時間差，
及びＥ２点とＥ３点の時間差が共にＣＬ／２未満である
なら，参照時刻に近い点Ｅ１，Ｅ２，Ｅ３の順に興奮が
伝わっていると判断し，参照時刻からの時間をそのまま
Ｅ１，Ｅ２，Ｅ３点の興奮波通過時刻として式（１）に
代入する。
及びＥ２点とＥ３点の時間差が共にＣＬ／２未満である
なら，参照時刻に近い点Ｅ１，Ｅ２，Ｅ３の順に興奮が
伝わっていると判断し，参照時刻からの時間をそのまま
Ｅ１，Ｅ２，Ｅ３点の興奮波通過時刻として式（１）に
代入する。
【００３８】図８Ｂの場合，参照時刻に近い順にＥ１，
Ｅ２，Ｅ３の興奮が見られるが，実際の興奮は，Ｅ２，
Ｅ３，Ｅ１′の順に起きていると考えられる。この場
合，Ｅ１点とＥ２点の時間差がＣＬ／２以上かどうかを
調べ，ＣＬ／２以上の場合には，Ｅ１点の代わりに１周
期後のＥ１′点を採用し，Ｅ２，Ｅ３，Ｅ１′の３点の
興奮波通過時刻を式（１）に代入する。
Ｅ２，Ｅ３の興奮が見られるが，実際の興奮は，Ｅ２，
Ｅ３，Ｅ１′の順に起きていると考えられる。この場
合，Ｅ１点とＥ２点の時間差がＣＬ／２以上かどうかを
調べ，ＣＬ／２以上の場合には，Ｅ１点の代わりに１周
期後のＥ１′点を採用し，Ｅ２，Ｅ３，Ｅ１′の３点の
興奮波通過時刻を式（１）に代入する。
【００３９】図８Ｃの場合，参照時刻に近い順にＥ１，
Ｅ２，Ｅ３の興奮が見られるが，実際の興奮は，Ｅ３，
Ｅ１′，Ｅ２′の順に起きていると考えられる。この場
合，Ｅ１点とＥ２点の時間差及びＥ２点とＥ３点の時間
差がＣＬ／２以上かどうかを調べ，Ｅ１点とＥ２点の時
間差のみがＣＬ／２以上の場合には，Ｅ１点，Ｅ２点の
代わりに１周期後のＥ１′点，Ｅ２′点を採用し，Ｅ
３，Ｅ１′，Ｅ２′の３点の興奮波通過時刻を式（１）
に代入する。
Ｅ２，Ｅ３の興奮が見られるが，実際の興奮は，Ｅ３，
Ｅ１′，Ｅ２′の順に起きていると考えられる。この場
合，Ｅ１点とＥ２点の時間差及びＥ２点とＥ３点の時間
差がＣＬ／２以上かどうかを調べ，Ｅ１点とＥ２点の時
間差のみがＣＬ／２以上の場合には，Ｅ１点，Ｅ２点の
代わりに１周期後のＥ１′点，Ｅ２′点を採用し，Ｅ
３，Ｅ１′，Ｅ２′の３点の興奮波通過時刻を式（１）
に代入する。
【００４０】図９は，上記補正を行う場合の一実施例を
フローチャートで示したものである。周期性による補正
を除き，図４のフローと基本的に同じである。
フローチャートで示したものである。周期性による補正
を除き，図４のフローと基本的に同じである。
【００４１】ステップ１１：複数の双極電極の座標マッ
プをメモリに登録する。 ステップ１２：不整脈の周期ＣＬ及び双極電極の各チャ
ンネル毎に興奮点までの興奮波通過時刻を入力する。興
奮波通過時刻は双極電極に対応して記憶される。
プをメモリに登録する。 ステップ１２：不整脈の周期ＣＬ及び双極電極の各チャ
ンネル毎に興奮点までの興奮波通過時刻を入力する。興
奮波通過時刻は双極電極に対応して記憶される。
【００４２】ステップ１３：双極電極の中の隣り合う双
極電極３点を選択する。 ステップ１４：図１０のフローにより周期性による補正
処理を行う。図１０の処理フローを説明すると，興奮波
通過時刻の早い方から順に，Ａｔ，Ｂｔ，Ｃｔとした場
合，Ｂｔ−Ａｔ＜ＣＬ／２且つＣｔ−Ａｔ＜ＣＬ／２か
どうか調べる（ステップ３１）。ＹＥＳなら補正なしの
Ａｔ，Ｂｔ，Ｃｔをベクトル計算式（１）に代入する。
ＮＯの場合，Ｂｔ−Ａｔ≧ＣＬ／２かどうか調べる（ス
テップ３２）。ＹＥＳの場合，ＡｔにＣＬを加えて新た
なＡｔとし（ステップ３３），補正処理後のＡｔ，Ｂ
ｔ，Ｃｔをベクトル計算式（１）に代入する（ステップ
３４）。ＮＯの場合，Ｂｔ−Ａｔ＜ＣＬ／２，Ｃｔ−Ｂ
ｔ≧ＣＬ／２であるので，Ａｔ，ＢｔそれぞれにＣＬを
加えて新たなＡｔ，Ｂｔとし（ステップ３５），補正処
理後のＡｔ，Ｂｔ，Ｃｔをベクトル計算式（１）に代入
する（ステップ３４）。
極電極３点を選択する。 ステップ１４：図１０のフローにより周期性による補正
処理を行う。図１０の処理フローを説明すると，興奮波
通過時刻の早い方から順に，Ａｔ，Ｂｔ，Ｃｔとした場
合，Ｂｔ−Ａｔ＜ＣＬ／２且つＣｔ−Ａｔ＜ＣＬ／２か
どうか調べる（ステップ３１）。ＹＥＳなら補正なしの
Ａｔ，Ｂｔ，Ｃｔをベクトル計算式（１）に代入する。
ＮＯの場合，Ｂｔ−Ａｔ≧ＣＬ／２かどうか調べる（ス
テップ３２）。ＹＥＳの場合，ＡｔにＣＬを加えて新た
なＡｔとし（ステップ３３），補正処理後のＡｔ，Ｂ
ｔ，Ｃｔをベクトル計算式（１）に代入する（ステップ
３４）。ＮＯの場合，Ｂｔ−Ａｔ＜ＣＬ／２，Ｃｔ−Ｂ
ｔ≧ＣＬ／２であるので，Ａｔ，ＢｔそれぞれにＣＬを
加えて新たなＡｔ，Ｂｔとし（ステップ３５），補正処
理後のＡｔ，Ｂｔ，Ｃｔをベクトル計算式（１）に代入
する（ステップ３４）。
【００４３】補正処理後ステップ１５に進む。ステップ
１５−１８は図４のステップ４−７に対応しているので
説明は省略する。図１１は一定の周期を持たない不整脈
を解析する場合の一実施例を示す。図中，各興奮点に示
された数値は参照時刻からの時間を表わす。
１５−１８は図４のステップ４−７に対応しているので
説明は省略する。図１１は一定の周期を持たない不整脈
を解析する場合の一実施例を示す。図中，各興奮点に示
された数値は参照時刻からの時間を表わす。
【００４４】心房細動，心室細動などの様に，一定の周
期を持たない不整脈の場合，同じ場所が２度興奮する時
には，その場所を通過する興奮波の方向は異なる。従っ
て，ある時間範囲の興奮波の状況をベクトルで表現する
ためには，その時間帯のすべての興奮の時刻を測定する
必要がある。そして，すべての時刻での興奮波ベクトル
を同時に表現すれば，ベクトルが重なり，理解できなく
なる。しかし，十分短い時間，例えば１００msの時間幅
を設定して興奮を観察するなら，この間に同じ場所が２
度興奮することはない。そこで，１００msの時間幅の解
析ウインドウを設定してこのウインドウ内の興奮点に関
して興奮波ベクトルを求めて表示する。次に，このウイ
ンドウを例えばΔＷ＝１０msシフトして，ベクトルを求
め，表示する。このように，解析ウインドウをΔＷづつ
シフトしながらそのウインドウ内のベクトルを動画的に
表示するなら，興奮波の経時的変化を表現することがで
きる。解析ウインドウの時間幅は，１０−４００msの間
に設定できる。好ましくは５０−２００msの範囲に設定
するとよい。図１２は一定の周期を持たない不整脈を解
析する場合の解析フローチャートの一実施例である。解
析ウインドウに関する処理の部分を除いては，基本的に
図４のフローと同じである。
期を持たない不整脈の場合，同じ場所が２度興奮する時
には，その場所を通過する興奮波の方向は異なる。従っ
て，ある時間範囲の興奮波の状況をベクトルで表現する
ためには，その時間帯のすべての興奮の時刻を測定する
必要がある。そして，すべての時刻での興奮波ベクトル
を同時に表現すれば，ベクトルが重なり，理解できなく
なる。しかし，十分短い時間，例えば１００msの時間幅
を設定して興奮を観察するなら，この間に同じ場所が２
度興奮することはない。そこで，１００msの時間幅の解
析ウインドウを設定してこのウインドウ内の興奮点に関
して興奮波ベクトルを求めて表示する。次に，このウイ
ンドウを例えばΔＷ＝１０msシフトして，ベクトルを求
め，表示する。このように，解析ウインドウをΔＷづつ
シフトしながらそのウインドウ内のベクトルを動画的に
表示するなら，興奮波の経時的変化を表現することがで
きる。解析ウインドウの時間幅は，１０−４００msの間
に設定できる。好ましくは５０−２００msの範囲に設定
するとよい。図１２は一定の周期を持たない不整脈を解
析する場合の解析フローチャートの一実施例である。解
析ウインドウに関する処理の部分を除いては，基本的に
図４のフローと同じである。
【００４５】ステップ４０：複数の双極電極の座標マッ
プをメモリ上に登録する。 ステップ４１：トータルの解析対象領域，解析ウイドウ
幅Ｗ，解析ウインドウのシフト量ΔＷを入力し，設定す
る。
プをメモリ上に登録する。 ステップ４１：トータルの解析対象領域，解析ウイドウ
幅Ｗ，解析ウインドウのシフト量ΔＷを入力し，設定す
る。
【００４６】ステップ４２：解析対象時間内において，
双極電極の各チャンネル毎に，参照時刻から各興奮点ま
での時間を入力し，メモリにストアする。
双極電極の各チャンネル毎に，参照時刻から各興奮点ま
での時間を入力し，メモリにストアする。
【００４７】ステップ４３：解析ウインドウＷ内の興奮
点を探す。もし，一つの部位で複数の興奮が認められる
場合は，その中の一つを選択する。例えば，解析ウイン
ドウの開始点又は終了点に最も近いものを選択できる。
点を探す。もし，一つの部位で複数の興奮が認められる
場合は，その中の一つを選択する。例えば，解析ウイン
ドウの開始点又は終了点に最も近いものを選択できる。
【００４８】ステップ４４−４８は，図４のステップ３
‐７に対応するので説明は省略する。
‐７に対応するので説明は省略する。
【００４９】ステップ４９：１つの解析ウイドウ内での
処理をすべて終了したら解析ウインドウをΔＷだけシフ
トする。
処理をすべて終了したら解析ウインドウをΔＷだけシフ
トする。
【００５０】ステップ５０：解析ウインドウの開始時刻
がトータルの解析対象領域の終了時刻を過ぎたかどうか
調べる。ＮＯならステップ４３に戻り，新たな解析ウイ
ンドウ内で解析処理を実行する。ＹＥＳなら解析処理を
終了する。
がトータルの解析対象領域の終了時刻を過ぎたかどうか
調べる。ＮＯならステップ４３に戻り，新たな解析ウイ
ンドウ内で解析処理を実行する。ＹＥＳなら解析処理を
終了する。
【００５１】図１３に，本発明の心内電位解析装置の一
実施例のブロック図を示す。６０はＣＰＵであり，装置
全体の動作をを制御する。６１は入出力インターフェイ
ス部であり，心筋測定データを外部の心筋測定装置から
受信したり，解析結果をプリンター等の外部装置に出力
したり，キーボードやマウスなどのデータ入力デバイス
との接続を可能にするものである。６２は記憶部であ
り，ＲＡＭやＲＯＭを備え，図４，図９，図１２に示さ
れた本発明の解析フローの各ステップを実行するプログ
ラムや，電極の座標マップや興奮波通過時刻等の入力さ
れたデータや解析途中のデータや解析結果等を記憶す
る。６３はデータ入出力部であり，磁気ディスク，光デ
ィスク，光−磁気ディスク，あるいは半導体メモリなど
のデータ記録媒体に対するデータの入出力を行う。６４
は興奮波ベクトル演算部であり，電極の３つの組に関す
る興奮波通過時刻と電極の座標マップの位置情報から解
析フローを実行する。６５は表示装置であり，解析結果
をベクトル表示する。
実施例のブロック図を示す。６０はＣＰＵであり，装置
全体の動作をを制御する。６１は入出力インターフェイ
ス部であり，心筋測定データを外部の心筋測定装置から
受信したり，解析結果をプリンター等の外部装置に出力
したり，キーボードやマウスなどのデータ入力デバイス
との接続を可能にするものである。６２は記憶部であ
り，ＲＡＭやＲＯＭを備え，図４，図９，図１２に示さ
れた本発明の解析フローの各ステップを実行するプログ
ラムや，電極の座標マップや興奮波通過時刻等の入力さ
れたデータや解析途中のデータや解析結果等を記憶す
る。６３はデータ入出力部であり，磁気ディスク，光デ
ィスク，光−磁気ディスク，あるいは半導体メモリなど
のデータ記録媒体に対するデータの入出力を行う。６４
は興奮波ベクトル演算部であり，電極の３つの組に関す
る興奮波通過時刻と電極の座標マップの位置情報から解
析フローを実行する。６５は表示装置であり，解析結果
をベクトル表示する。
【００５２】尚，図４，図９，図１２に示された本発明
の解析フローの各ステップを実行するプログラムはコン
ピュータにより読み取り可能な記録媒体，例えば，磁気
ディスク，光ディスク，光−磁気ディスク，あるいは半
導体メモリなどに記録されて提供可能である。この場
合，本発明の解析プログラムは，コンピュータにより読
み取り可能な記録媒体を読み取る装置を有する入出イン
ターフェイス部，ＣＰＵ，メモリ，表示装置などを備え
たコンピュータにて実施できる。
の解析フローの各ステップを実行するプログラムはコン
ピュータにより読み取り可能な記録媒体，例えば，磁気
ディスク，光ディスク，光−磁気ディスク，あるいは半
導体メモリなどに記録されて提供可能である。この場
合，本発明の解析プログラムは，コンピュータにより読
み取り可能な記録媒体を読み取る装置を有する入出イン
ターフェイス部，ＣＰＵ，メモリ，表示装置などを備え
たコンピュータにて実施できる。
【００５３】実施例では，解析結果を平面上に表示する
例を示したが，図３に相当する双極電極の座標マップを
３次元空間上に設定し，解析結果を３次元的に表示する
ことも可能である。また，電極を双極電極として説明し
たが，単極電極の場合にも本発明はそのまま適応でき
る。
例を示したが，図３に相当する双極電極の座標マップを
３次元空間上に設定し，解析結果を３次元的に表示する
ことも可能である。また，電極を双極電極として説明し
たが，単極電極の場合にも本発明はそのまま適応でき
る。
【００５４】図１４は解析結果を格納した記録媒体のデ
ータ構造図である。この記録媒体を他のコンピュータに
入れてより読み出すことで，解析結果を別の場所でも見
ることが可能になる。記録媒体としては，磁気ディス
ク，光ディスク，光−磁気ディスク，あるいは半導体メ
モリなどがある。
ータ構造図である。この記録媒体を他のコンピュータに
入れてより読み出すことで，解析結果を別の場所でも見
ることが可能になる。記録媒体としては，磁気ディス
ク，光ディスク，光−磁気ディスク，あるいは半導体メ
モリなどがある。
【００５５】７１は解析結果を表示させるためのプログ
ラムを読み出すコード部，７２は電極の座標マップデー
タを格納しておく領域，７３は計算で求められた興奮波
ベクトルデータを格納しておく領域，７４はベクトルを
電極の座標マップに関連付けて表示するためのコード部
である。
ラムを読み出すコード部，７２は電極の座標マップデー
タを格納しておく領域，７３は計算で求められた興奮波
ベクトルデータを格納しておく領域，７４はベクトルを
電極の座標マップに関連付けて表示するためのコード部
である。
【００５６】
【発明の効果】本発明により，局所の興奮波の進行方
向，進行速度を示すベクトルの集合によって興奮の開始
時刻，終了時刻には全く影響されずに心内興奮伝播経路
が表現される。各興奮点の時間差に基づいて興奮波ベク
トルが計算されるため，興奮波通過時刻の測定のために
心電図上に設定した参照時刻の位置にも影響されない。
これにより，さまざまな不整脈，主にリエントリーによ
る不整脈の診断，治療が容易となる。更に本発明では，
興奮波の進行方向，進行速度がベクトルとして直接表現
されており，不整脈のメカニズムを推定することが容易
となる。
向，進行速度を示すベクトルの集合によって興奮の開始
時刻，終了時刻には全く影響されずに心内興奮伝播経路
が表現される。各興奮点の時間差に基づいて興奮波ベク
トルが計算されるため，興奮波通過時刻の測定のために
心電図上に設定した参照時刻の位置にも影響されない。
これにより，さまざまな不整脈，主にリエントリーによ
る不整脈の診断，治療が容易となる。更に本発明では，
興奮波の進行方向，進行速度がベクトルとして直接表現
されており，不整脈のメカニズムを推定することが容易
となる。
【００５７】本発明は開示した実施例に限定されるもの
ではない。特許請求の範囲の精神内に入るさまざまな変
形例も本発明に含まれる。
ではない。特許請求の範囲の精神内に入るさまざまな変
形例も本発明に含まれる。
【図１】Ａは，心筋興奮の測定に使用される電極の図。
Ｂは，図１Ａの電極の配置を平面上に表した図。
Ｂは，図１Ａの電極の配置を平面上に表した図。
【図２】双極電極による測定結果を示す図。
【図３】双極電極を平面上に表した図。
【図４】本発明による解析フローチャートの一実施例を
示す図。
示す図。
【図５】興奮波ベクトルを描く様子を示した図。
【図６】図２の測定データの解析結果を示した図。
【図７】興奮伝播の経路を推定した一例を示す図。
【図８】周期性のある不整脈を解析する場合の興奮波通
過時刻に対する補正原理を説明する図。
過時刻に対する補正原理を説明する図。
【図９】周期性のある不整脈を解析する場合の解析フロ
ーチャートの一実施例を示す図。
ーチャートの一実施例を示す図。
【図１０】補正フローチャートの一実施例を示す図。
【図１１】周期性のない不整脈を解析する場合の一実施
例を示す図。
例を示す図。
【図１２】周期性のない不整脈を解析する場合の解析フ
ローチャートの一実施例を示す図。
ローチャートの一実施例を示す図。
【図１３】本発明の心内電位解析装置の一実施例のブロ
ック図。
ック図。
【図１４】解析結果としての興奮波ベクトルの表示プロ
グラムを格納した記録媒体のデータ構造を説明するため
の図。
グラムを格納した記録媒体のデータ構造を説明するため
の図。
【図１５】図２のリエントリーによる不整脈を開始時刻
を異ならせて等時線図を描いた場合の結果を示す図。
を異ならせて等時線図を描いた場合の結果を示す図。
６０ ＣＰＵ ６１ Ｉ／Ｏインターフェイス ６２ 記憶部 ６３ データ入出力部 ６４ 興奮波ベクトル演算部 ６５ 表示部
───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.7，ＤＢ名) A61B 5/042 A61B 5/046 A61B 5/0464
Claims (25)
- 【請求項１】 コンピュータを利用した心内電位解析装
置であって，該装置は，複数の電極の相対的位置関係を
座標マップとして記憶する第１の手段と，前記複数の電
極のうち，隣り合う３つの電極の組を複数選択する手段
と，前記隣り合う３つの電極の複数の組の各々につい
て，前記隣り合う３つの電極を通過した興奮波のベクト
ルを求める手段と，求められた前記複数のベクトルを前
記座標マップに関連付けて表示する手段と，を含むこと
を特徴とした心内電位解析装置。 - 【請求項２】 請求項１の心内電位解析装置において，
前記ベクトルを求める手段は，心電図の時間軸に任意に
設定された基準となる参照時刻から前記各電極を興奮波
が通過するまでの興奮波通過時刻を記憶する第２の手段
と，前記隣り合う３つの電極に関し，前記第１の手段か
ら読み出された位置情報と前記第２の手段から読み出さ
れた前記興奮波通過時刻とから興奮波のベクトルを求め
る手段と，を含むことを特徴とした心内電位解析装置。 - 【請求項３】 請求項２の心内電位解析装置において，
前記ベクトルを求める手段は，下記の式で計算すること
を特徴とする；互いに直交するｘ，ｙ，ｚ軸で定義され
た直交３次元座標系に３つの電極の位置をＡ（Ａｘ，Ａ
ｙ，Ａｚ），Ｂ（Ｂｘ，Ｂｙ，Ｂｚ），Ｃ（Ｃｘ，Ｃ
ｙ，Ｃｚ）と設定する。前記隣り合う３つの電極におけ
る興奮波通過時刻をそれぞれＡｔ，Ｂｔ，Ｃｔとする
と，興奮波ベクトル 【数１】 は下記の式で求められる。 【数２】 但し， 【数３】 - 【請求項４】 請求項３の心内電位解析装置において，
前記ベクトルを求める手段は，前記選択された３つの電
極を興奮が実際に通過した順に対応して前記興奮波通過
時刻Ａｔ，Ｂｔ，Ｃｔを求める手段を含む。 - 【請求項５】 請求項３の心内電位解析装置において，
前記興奮波通過時刻Ａｔ，Ｂｔ，Ｃｔを求める手段は，
不整脈が周期ＣＬを有する場合，興奮波通過時刻の短い
ほうから順に，Ａｔ，Ｂｔ，Ｃｔとした場合，Ｂｔ−Ａ
ｔ≧ＣＬ／２且つＣｔ−Ｂｔ＜ＣＬ／２の場合，Ａｔに
ＣＬを加えて新たなＡｔとして前記べクトル計算式に代
入し，Ｂｔ−Ａｔ＜ＣＬ／2且つＣｔ−Ｂｔ≧ＣＬ／２
の場合，Ａｔ，ＢｔそれぞれにＣＬを加えて新たなＡ
ｔ，Ｂｔとして前記ベクトル計算式に代入する手段を含
む。 - 【請求項６】 請求項２の心内電位解析装置において，
前記電極の組を複数選択する手段は，所定の時間幅内の
興奮点に対応した複数の電極から選択する手段を含む。 - 【請求項７】 請求項６の心内電位解析装置において，
前記電極の組を複数選択する手段は，前記所定の時間幅
を該所定の時間幅より短い時間をシフト量として順次シ
フトしながら前記電極の組を選択する手段を含み，前記
表示する手段は，前記所定の時間幅毎に，求められた前
記複数のベクトルを順次動画的に表示する手段を含む。 - 【請求項８】 請求項１の心内電位解析装置において，
前記表示する手段は，計算で求められた前記興奮波の速
度が所定の範囲に入っているものを選択して表示する手
段を含む。 - 【請求項９】 コンピュータを利用した心内電位解析方
法であって，該方法は，複数の電極の相対的位置関係を
座標マップとしてメモリに登録するステップと，前記複
数の電極のうち，隣り合う３つの電極の組を複数選択す
るステップと，前記隣り合う３つの電極の複数の組の各
々について，前記隣り合う３つの電極を通過した興奮波
のベクトルを求めるステップと，求められた前記複数の
ベクトルを前記座標マップに関連付けて表示するステッ
プと，を含むことを特徴とする心内電位解析方法。 - 【請求項１０】 請求項９の心内電位解析方法におい
て，前記ベクトルを求めるステップは，前記選択された
隣り合う３つの電極の位置情報と，心電図の時間軸に任
意に設定された基準となる参照時刻から前記選択された
３つの電極を興奮波が通過するまでの興奮波通過時刻と
から，前記興奮波のベクトルを求めるステップを含むこ
とを特徴とする心内電位解析方法。 - 【請求項１１】 請求項１０の心内電位解析方法におい
て，前記ベクトルを求めるステップは，下記の式で計算
することを特徴とする；互いに直交するｘ，ｙ，ｚ軸で
定義された直交３次元座標系に３つの電極の位置をＡ
（Ａｘ，Ａｙ，Ａｚ），Ｂ（Ｂｘ，Ｂｙ，Ｂｚ），Ｃ
（Ｃｘ，Ｃｙ，Ｃｚ）と設定する。前記隣り合う３つの
電極における興奮波通過時刻をそれぞれＡｔ，Ｂｔ，Ｃ
ｔとすると，興奮波ベクトル 【数４】 は下記の式で求められる。 【数５】 但し， 【数６】 - 【請求項１２】 請求項１１の心内電位解析方法におい
て，前記ベクトルを求めるステップは，前記選択された
３つの電極を興奮が実際に通過した順に対応して前記興
奮波通過時刻Ａｔ，Ｂｔ，Ｃｔを求めるステップを含
む。 - 【請求項１３】 請求項１２の心内電位解析方法におい
て，前記興奮波通過時刻Ａｔ，Ｂｔ，Ｃｔを求めるステ
ップは，不整脈が周期ＣＬを有する場合，興奮波通過時
刻の短いほうから順に，Ａｔ，Ｂｔ，Ｃｔとした場合，
Ｂｔ−Ａｔ≧ＣＬ／２且つＣｔ−Ｂｔ＜ＣＬ／２の場
合，ＡｔにＣＬを加えて新たなＡｔとして前記ベクトル
計算式に代入し，Ｂｔ−Ａｔ＜ＣＬ／２且つＣｔ−Ｂｔ
≧ＣＬ／２の場合，Ａｔ，ＢｔそれぞれにＣＬを加えて
新たなＡｔ，Ｂｔとして前記ベクトル計算式に代入する
ステップを含む。 - 【請求項１４】 請求項９の心内電位解析方法におい
て，前記電極の組を複数選択するステップは，所定の時
間幅内の興奮点に対応した複数の電極から選択するステ
ップを含む。 - 【請求項１５】 請求項１４の心内電位解析方法におい
て，前記電極の組を複数選択するステップは，前記所定
の時間幅を該所定の時間幅より短い時間をシフト量とし
て順次シフトしながら前記電極の組を選択するステップ
を含み，前記表示するステップは，前記所定の時間幅毎
に，求められた前記複数のベクトルを順次動画的に表示
するステップを含む。 - 【請求項１６】 請求項９の心内電位解析方法におい
て，前記表示するステップは，計算で求められた前記興
奮波の速度が所定の範囲に入っているものを選択して表
示するステップを含む。 - 【請求項１７】 複数の電極の相対的位置関係を座標マ
ップとしてメモリに登録するステップと，前記複数の電
極のうち，隣り合う３つの電極の組を複数選択するステ
ップと，前記隣り合う３つの電極の複数の組の各々につ
いて，前記隣り合う３つの電極を通過した興奮波のベク
トルを求めるステップと，求められた前記複数のベクト
ルを前記座標マップに関連付けて表示するステップと，
を有する心内電位解析方法をコンピュータに実行させる
ためのプログラムを記録したコンピュータ読み取り可能
な媒体。 - 【請求項１８】 請求項１７のコンピュータ読み取り可
能な媒体において，前記ベクトルを求めるステップは，
前記選択された隣り合う３つの電極の位置情報と，心電
図の時間軸に任意に設定された基準となる参照時刻から
前記選択された３つの電極を興奮波が通過するまでの興
奮波通過時刻とから，前記興奮波のベクトルを求めるス
テップを含むことを特徴とする前記媒体。 - 【請求項１９】 請求項１８のコンピュータ読み取り可
能な媒体において，前記ベクトルを求めるステップは，
下記の式で計算することを特徴とする；互いに直交する
ｘ，ｙ，ｚ軸で定義された直交３次元座標系に３つの電
極の位置をＡ（Ａｘ，Ａｙ，Ａｚ），Ｂ（Ｂｘ，Ｂｙ，
Ｂｚ），Ｃ（Ｃｘ，Ｃｙ，Ｃｚ）と設定する。前記隣り
合う３つの電極における興奮波通過時刻をそれぞれＡ
ｔ，Ｂｔ，Ｃｔとすると，興奮波ベクトル 【数７】 は下記の式で求められる。 【数８】 但し， 【数９】 - 【請求項２０】 請求項１９のコンピュータ読み取り可
能な媒体において，前記ベクトルを求めるステップは，
前記選択された３つの電極を興奮が実際に通過した順に
対応して前記興奮波通過時刻Ａｔ，Ｂｔ，Ｃｔを求める
ステップを含む。 - 【請求項２１】 請求項２０のコンピュータ読み取り可
能な媒体において，前記興奮波通過時刻Ａｔ，Ｂｔ，Ｃ
ｔを求めるステップは，不整脈が周期ＣＬを有する場
合，興奮波通過時刻の短いほうから順に，Ａｔ，Ｂｔ，
Ｃｔとした場合，Ｂｔ−Ａｔ≧ＣＬ／２且つＣｔ−Ｂｔ
＜ＣＬ／２の場合，ＡｔにＣＬを加えて新たなＡｔとし
て前記ベクトル計算式に代入し，Ｂｔ−Ａｔ＜ＣＬ／２
且つＣｔ−Ｂｔ≧ＣＬ／２の場合，Ａｔ，Ｂｔそれぞれ
にＣＬを加えて新たなＡｔ，Ｂｔとして前記ベクトル計
算式に代入するステップを含む。 - 【請求項２２】 請求項１７のコンピュータ読み取り可
能な媒体において，前記電極の組を複数選択するステッ
プは，所定の時間幅内の興奮点に対応した複数の電極か
ら選択するステップを含む。 - 【請求項２３】 請求項２２のコンピュータ読み取り可
能な媒体において，前記電極の組を複数選択するステッ
プは，前記所定の時間幅を該所定の時間幅より短い時間
をシフト量として順次シフトしながら前記電極の組を選
択するステップを含み，前記表示するステップは，前記
所定の時間幅毎に，求められた前記複数のベクトルを順
次動画的に表示するステップを含む。 - 【請求項２４】 請求項１７のコンピュータ読み取り可
能な媒体において，前記表示するステップは，計算で求
められた前記興奮波の速度が所定の範囲に入っているも
のを選択して表示するステップを含む。 - 【請求項２５】 記録されたデータを読み出すためのプ
ログラムを記憶した第１の領域と，複数の電極の相対的
位置関係を記憶した第２の領域と，前記電極の内，隣り
合う３つの電極の複数の組に関連して求められた興奮波
のベクトルデータを記憶した第３の領域と，前記ベクト
ルデータに基づく複数のベクトルを，該ベクトルを求め
るために使用した前記３つの電極の位置に関連付けて表
示装置に表示させるためのプログラムを記憶した第４の
領域と，を含むコンピュータ読み取り可能なデータ記録
媒体。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP21194699A JP3213602B2 (ja) | 1999-07-27 | 1999-07-27 | 心内電位解析装置及び方法 |
| US09/405,806 US6360121B1 (en) | 1999-07-27 | 1999-09-24 | Intracardiac potential analyzing apparatus and method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP21194699A JP3213602B2 (ja) | 1999-07-27 | 1999-07-27 | 心内電位解析装置及び方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2001037730A JP2001037730A (ja) | 2001-02-13 |
| JP3213602B2 true JP3213602B2 (ja) | 2001-10-02 |
Family
ID=16614323
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP21194699A Expired - Fee Related JP3213602B2 (ja) | 1999-07-27 | 1999-07-27 | 心内電位解析装置及び方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US6360121B1 (ja) |
| JP (1) | JP3213602B2 (ja) |
Families Citing this family (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8676303B2 (en) | 2008-05-13 | 2014-03-18 | The Regents Of The University Of California | Methods and systems for treating heart instability |
| JP5991706B2 (ja) | 2008-10-09 | 2016-09-14 | ザ リージェンツ オブ ザ ユニバーシティ オブ カリフォルニア | バイオリズム障害の検出、診断、及び治療のための方法、システム、及び装置 |
| US9332915B2 (en) | 2013-03-15 | 2016-05-10 | The Regents Of The University Of California | System and method to identify sources associated with biological rhythm disorders |
| US9392948B2 (en) | 2011-12-09 | 2016-07-19 | The Regents Of The University Of California | System and method of identifying sources for biological rhythms |
| US10434319B2 (en) | 2009-10-09 | 2019-10-08 | The Regents Of The University Of California | System and method of identifying sources associated with biological rhythm disorders |
| US10398326B2 (en) | 2013-03-15 | 2019-09-03 | The Regents Of The University Of California | System and method of identifying sources associated with biological rhythm disorders |
| CA2795770A1 (en) | 2010-04-08 | 2011-10-13 | Topera, Inc. | Methods, system and apparatus for the detection, diagnosis and treatment of biological rhythm disorders |
| US9050006B2 (en) | 2011-05-02 | 2015-06-09 | The Regents Of The University Of California | System and method for reconstructing cardiac activation information |
| US9107600B2 (en) | 2011-05-02 | 2015-08-18 | The Regents Of The University Of California | System and method for reconstructing cardiac activation information |
| US8165666B1 (en) | 2011-05-02 | 2012-04-24 | Topera, Inc. | System and method for reconstructing cardiac activation information |
| WO2012151301A1 (en) | 2011-05-02 | 2012-11-08 | Topera, Inc. | System and method for targeting heart rhythm disorders using shaped ablation |
| US8715199B1 (en) | 2013-03-15 | 2014-05-06 | Topera, Inc. | System and method to define a rotational source associated with a biological rhythm disorder |
| CN105307565B (zh) | 2013-05-07 | 2019-03-15 | 圣犹达医疗用品电生理部门有限公司 | 利用电极空间布置以用于表征心脏传导状况 |
| EP3721796B1 (en) * | 2014-02-25 | 2023-11-15 | St. Jude Medical, Cardiology Division, Inc. | System for local electrophysiological characterization of cardiac substrate using multi-electrode catheters |
| JP2016063853A (ja) * | 2014-09-22 | 2016-04-28 | フクダ電子株式会社 | 電気生理学的検査装置 |
| EP3270776B1 (en) | 2015-05-12 | 2020-05-27 | St. Jude Medical, Cardiology Division, Inc. | Systems and methods for orientation independent sensing |
| WO2017041892A1 (en) | 2015-09-07 | 2017-03-16 | Ablacon Inc. | System for analyzing electrophysiological data and method for analyzing electrophysiological data |
| WO2018191686A1 (en) | 2017-04-14 | 2018-10-18 | St. Jude Medical, Cardiology Division, Inc. | Orientation independent sensing, mapping, interface and analysis systems and methods |
| JP7104177B2 (ja) * | 2018-04-26 | 2022-07-20 | セント・ジュード・メディカル，カーディオロジー・ディヴィジョン，インコーポレイテッド | 不整脈ドライバ部位をマッピングするためのシステム及び方法 |
| WO2020096689A1 (en) | 2018-09-10 | 2020-05-14 | St. Jude Medical, Cardiology Division, Inc. | System and method for displaying electrophysiological signals from multi-dimensional catheters |
| JP7186395B2 (ja) * | 2019-01-25 | 2022-12-09 | 国立大学法人滋賀医科大学 | 記録装置、記録方法およびプログラム |
| CN115605136A (zh) | 2020-05-19 | 2023-01-13 | 圣犹达医疗用品心脏病学部门有限公司（Us） | 用于标测电生理激活的***和方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0278215A1 (de) * | 1987-01-22 | 1988-08-17 | Studer Revox Ag | Verfahren und Vorrichtung zur Verarbeitung von Signalen aus einem physiologisch erzeugten elektrischen Feld |
| US5595183A (en) * | 1995-02-17 | 1997-01-21 | Ep Technologies, Inc. | Systems and methods for examining heart tissue employing multiple electrode structures and roving electrodes |
| WO1999005964A1 (en) * | 1997-08-01 | 1999-02-11 | Harbinger Medical, Inc. | System and method of noninvasevely determining a patient's susceptibility to arrhythmia |
| US6115630A (en) * | 1999-03-29 | 2000-09-05 | Medtronic, Inc. | Determination of orientation of electrocardiogram signal in implantable medical devices |
-
1999
- 1999-07-27 JP JP21194699A patent/JP3213602B2/ja not_active Expired - Fee Related
- 1999-09-24 US US09/405,806 patent/US6360121B1/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2001037730A (ja) | 2001-02-13 |
| US6360121B1 (en) | 2002-03-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3213602B2 (ja) | 心内電位解析装置及び方法 | |
| EP3517031B1 (en) | Intra-cardiac scar tissue identification using impedance sensing and contact measurment | |
| US9050011B2 (en) | Removal of artifacts from map data | |
| JP6866207B2 (ja) | 心房細動のマッピング | |
| CN104434299B (zh) | 用于在窦性节律期间标测室性/房性早搏的方法 | |
| JP6830805B2 (ja) | 電気解剖学的マップ上における多重興奮到達区域の表示 | |
| JP6719885B2 (ja) | 心内信号を利用した位置合わせマップ | |
| CN100393283C (zh) | 生物体组织动状态跟踪方法、使用该方法的图像诊断装置 | |
| JP6818464B2 (ja) | 心臓興奮の時系列マッチング及びトラッキング | |
| JP2016039901A (ja) | ブロックのラインの検出 | |
| US11160481B2 (en) | Atrial fibrillation mapping using atrial fibrillation cycle length (AFCL) gradients | |
| EP3530184B1 (en) | Automatic identification of multiple activation pathways | |
| JP4847684B2 (ja) | 画像診断装置 | |
| CN109068984B (zh) | 心肌兴奋判定装置 | |
| US20230355159A1 (en) | Detecting potential slow-conduction cardiac tissue areas in stable arrhythmias | |
| US20200155028A1 (en) | System and Method for Characterizing Arrhythmias | |
| JP2002282236A (ja) | 時間経過につれて変形可能な臓器の変形に従う画像処理方法 | |
| EP4311490A1 (en) | Method and system for identification of fractionated signals | |
| JP7490339B2 (ja) | 心房細動電位図における興奮の特定 | |
| JP2021094388A (ja) | 二重電位問題の解決方法 | |
| JP2022165414A (ja) | 心室頻拍（ｖｔ）に関連する遅い電気生理学的（ｅｐ）心臓経路のアノテーション | |
| JP2023038932A (ja) | シグモイド曲線を用いて投影された電気生理学的波速度の重み付け | |
| Neuzil et al. | 598 Electrophysiology–Diagnosing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 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 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080719 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090719 Year of fee payment: 8 |
|
| LAPS | Cancellation because of no payment of annual fees |