JPH0622915A - Ecg信号に基づく心臓不整脈検出装置 - Google Patents

Ecg信号に基づく心臓不整脈検出装置

Info

Publication number
JPH0622915A
JPH0622915A JP3270038A JP27003891A JPH0622915A JP H0622915 A JPH0622915 A JP H0622915A JP 3270038 A JP3270038 A JP 3270038A JP 27003891 A JP27003891 A JP 27003891A JP H0622915 A JPH0622915 A JP H0622915A
Authority
JP
Japan
Prior art keywords
qrs
cluster
group
normal
value
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.)
Granted
Application number
JP3270038A
Other languages
English (en)
Other versions
JP3224569B2 (ja
Inventor
Eru Ripurii Kenesu
ケネス・エル・リプリー
Teii Kapuran Danieru
ダニエル・ティー・カプラン
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.)
Nippon Colin Co Ltd
Original Assignee
Nippon Colin Co 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 Nippon Colin Co Ltd filed Critical Nippon Colin Co Ltd
Publication of JPH0622915A publication Critical patent/JPH0622915A/ja
Application granted granted Critical
Publication of JP3224569B2 publication Critical patent/JP3224569B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

(57)【要約】 【目的】 不整脈が正確に検出される不整脈検出装置を
提供する。 【構成】 不正脈検出装置が、ECG装置から少なくと
も1個の連続アナログ信号を取得する手段と、そのアナ
ログ信号に基づいて少なくとも1個のデジタル信号を発
生させる手段と、その少なくとも1個のデジタル信号か
ら複数のスカラ信号を発生させる手段と、そのスカラ信
号から特徴情報を抽出する手段と、その抽出特徴情報の
種類の数に等しい次元を持つ特徴情報空間にその抽出特
徴情報をプロットする手段とを含む。正常なQRS群
が、前記特徴情報空間においてQRS特徴情報のクラス
タ内に位置付けられたQRS群の個数に基づいて判別さ
れる。その正常QRS群の判別後に取得されるQRS群
の各々は、複数のルールと、その判別前後の両方の正常
QRS群に関してのその判別後の各QRS群の位置とに
基づいてグループ分けされる。

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は、心電計ないしECG装
置に関し、特に、生体の心臓血管系から発生する信号に
含まれるQRS群の特徴情報に基づいて不整脈を検出す
る装置に関する。
【0002】
【従来の技術】図1は、典型的なECG波形の一部であ
って、2つの正常な心拍と1つの異常な心拍を示してい
る。少数の数字や特徴情報を用いることによりECG信
号中のQRS形状を数量化し、心室性期外収縮(PV
C)等の検出や心臓不整脈の数量化を行う種々の方法お
よび装置が既に知られている。一般に、QRS波形情報
を利用して心拍を類似する群にグループ分けすることに
より、形状的に類似する心拍から得られる特徴情報を、
任意の心拍が正常なものであるかあるいはPVC等の異
常なものであるかの判別に利用することができる。既存
のQRS特徴情報抽出装置では、最初の工程としてQR
Sが検出される。すなわち、まずQRS候補が検出さ
れ、該検出期時間を含むウインドウ内に含まれるECG
データの一部が取り出される。次いで、QRSパターン
ベクトル等の当該ウインドゥ内のECGデータが、QR
S特徴情報ベクトル等のQRS形状を表す特徴情報に変
換される。これまでに、上記QRSパターンベクトルを
少ない種類の特徴情報に変換する種々の技法が提案、実
施されてきた。例えば、米国特許第4,336,810
号によって開示されるECG記録の不整脈分析方法/装
置では、磁気記録されたECGテープがスキャンされて
QRS群が検出される。スキャンされたECG信号はデ
ジタル変換され、テンプレートとしてグループ分けされ
た既知のQRSパターンと比較される。今回比較された
QRS群に類似するパターンが見つからない場合には、
新たにテンプレートが作成され、次のQRS群がスキャ
ンされグループ化される。各QRS群のグループ化は、
QRS群の幅とQRS群の基線の上下の面積とを含むQ
RS群の形状の諸特徴情報に基づいて行われる。各QR
S群はテンプレートを利用したグループ化により、正
常、上室性期外収縮、心室性期外収縮ないしは不明(un
known )と判別される。米国特許4,583,553号
が開示する装着携行型ECG分析/記録装置において
は、患者のECG信号がデジタル化された後ハードウエ
ア回路によるフィルタにかけられ、次いで、QRS群の
検出、グループ化、類似するQRS群の集団への組み入
れおよびその結果に基づく報告のためのアルゴリズムへ
の入力として使用される。このアルゴリズムは、リアル
タイムで作動し、1人の患者について43種類(event
s)の検出ができるように設計されている。また、2チ
ャネルのECGを使用することによりグループ化に伴う
エラーが最小限に抑制される。米国特許4,742,4
58号は、ECG信号をフィルタにかけハードウエア回
路によってそのECG信号に含まれるQRS群を分析す
るパターン認識分析方法/装置を開示する。それぞれが
QRS群の特徴情報である傾きの(符号)変化点とその
(符号)変化点間の間隔とを利用することにより、その
ECG信号を特定する指標(signiture )が作成され、
QRS群の検出、グループ化の際に、そのQRS群が類
似するQRS群のグループに組み込まれ、そのような適
合するグループがない場合には新しくグループが作成さ
れる。米国特許第4,589,420号は、ECG信号
に含まれるQRS群を検出しそれをグループ化して患者
の状態の特徴を決定するアルゴリズムを使用するECG
律動分析方法/装置を開示する。QRS群幅、R−R間
隔、瞬間と平均の心拍数等の幾つもの異なったQRS群
の特徴情報が、ECG信号から抽出される。QRS群の
認定とグループ化とは、そのQRS群に含まれるノイズ
の量に基づいて行われる。心拍は2つの段階を経てグル
ープ化され、まず、新しいQRS群の各々がそれ以前の
QRS群と比較され、次いでその比較結果に基づいて有
限状態機械finite state machine) 処理により検査され
る。
【0003】
【発明が解決しようとする課題】本発明の主な目的は、
QRS群を分類することにより正常なQRS群と異常な
QRS群とを区別する装置を提供することである。本発
明の他の目的は、N次元特徴情報空間にプロットされた
QRS群候補のN種類の特徴情報に基づいて正常なQR
S群を判別する装置を提供することである。本発明の他
の目的は、特徴情報空間内に位置するQRS群のクラス
タに基づいて、ECG装置から得られるQRS群を正常
なものと異常なものとに区別する装置を提供することで
ある。本発明の他の目的は、特徴情報空間に位置するク
ラスタに含まれるQRS群の個数に基づいてQRS群を
正常なものと異常なものとに区別する装置を提供するこ
とである。本発明の他の目的は、ノイズおよび基線ドリ
フトの影響を受けることなくQRS群を位置付ける装置
を提供することである。
【0004】
【課題を解決するための手段】上記のおよびその他の目
的および効果は、本発明に従い、ECG装置から得られ
る少なくとも1個の連続アナログ信号を取得し、該アナ
ロク信号に基づいて少なくとも1個のデジタル信号を発
生させ、その少なくとも1個のデジタル信号から複数の
スカラ信号を発生させ、該スカラ信号から複数の特徴情
報を抽出し、抽出された特徴情報の数に等しい次元を持
つ特徴情報空間内に該抽出特徴情報をプロットすること
により達成される。正常なQRS群は、特徴情報空間中
のQRS特徴情報のクラスタに含まれるQRS群の個数
と発生時間とに基づいて判別される。正常なQRS群の
判別の後に得られる後続のQRS群は、その正常QRS
群との相対位置関係等に基づいてグループ分けされる。
【0005】
【作用および発明の効果】本発明によれば、N次元の特
徴情報空間内に位置するQRS群のクラスタに含まれる
QRS群の個数と発生時間とに基づいて、QRS群を正
常なものと異常なものとにグループ分けされる。そし
て、正常なグループに属するQRS群を有する信号から
不整脈の有無が判断される。このため、正確に不整脈が
検出される。
【0006】
【実施例】本発明によって達成される上記のおよびその
他の目的、特徴および効果は以下の実施例の説明によっ
て開示され、明確にされる。なお、以下の実施例の説明
のために参照される図において、対応する要素は一般に
同一符号によって指示される。図2には本発明の不整脈
検出装置1が示されている。装置1は、既存のECG装
置に接続されその装置からECG信号を受け取る第1、
第2、第3のチャンネルコネクタ10、12、14を含
む。装置1はまた、取得されたECG信号をマルチプレ
クサ(MUX)16とフィルタ−バンク18を介して受
け取ってデジタル信号に変換するアナログ/デジタル変
換回路(ADC)20と、そのデシタル信号を処理し定
量化するマイクロプロセッサ−22と、注記付きのEC
Gデ−タを表示する表示器32とを含む。マイクロプロ
セッサ−22は、前記デシタル信号を受信しまた処理さ
れたECGデ−タを出力する入出力(I/O)回路24
と、中央処理装置(CPU)26と、プログラム指令お
よび予め定められた定数とを記憶するROM28と、デ
ジタル信号とプログラムの中間結果と可変の定数(vari
able constants)とを記憶するRAM30とを含む。
【0007】I/O回路24は、ECG信号を集めるタ
イミングを制御するタイミング信号TをMUX16に出
力する。信号Tにより、各チャネルコネクタ10、1
2、14から約250ヘルツ(Hz)の頻度でECG信
号がサンプリングされる。信号Tによるサンプリング頻
度はフィルタバンク18に設けられる後述のアンチエイ
リアシングフィルタ(anti-aliasing )と共働するよう
に300Hzに設定されてもよい。ADC回路20は、
単位分解能が2.4マイクロボルト未満であって公称作
動電圧が±5ボルトである少なくとも12ビットのデ−
タを発生させる従来のADC回路である。
【0008】フィルタ−バンク18は、受け取ったEC
G信号から高周波ノイズを除去するための従来のロ−パ
スアンチエイリアシングフィルタを含む。フィルタ−バ
ンク18はさらに、受け取ったECG信号からそれに無
関係なACライン周波数(extraneous AC line frequen
cies)の誘導ノイズとその高調波ノイズを取り除くため
の従来のノッチフィルタを備えてもよい。このACライ
ン周波数ノッチフィルタをソフトウエアにて構成し、マ
イクロプロセッサ22によって実行するようにしてもよ
い。表示器32は、ECGデータ、つまり、マーク付き
のECG波形をチャート紙上に記録する従来の記録装置
である。表示器32はマイクロプロセッサ22と共働し
てECGデ−タを表に表すようにすることもできる。
【0009】図2の装置1は、3個のアナログECG信
号を同時に取得する手段としてのチャネルコネクタ1
0、12、14を備える。これら3個のECG信号を供
給するECG装置は、従来と同様のECG装置(図示せ
ず)である。ECG装置から3個のECG信号を取得で
きない場合に、マイクロプロセッサ22は、1個ないし
2個の信号を出力するように使用され得る。2個のEC
G信号が利用可能な場合には、マイクロプロセッサ22
はその2個のECG信号の差に基づいて補間ECG信号
を発生させる。また、1個のECG信号のみが利用可能
な場合には、マイクロプロセッサ22はその信号の0.
5倍に等しい第1補外信号とその信号の−10倍に等し
い第2補外ECG信号と発生させる。
【0010】以下、本発明の第1実施例の不整脈検出装
置の作動を、図2乃至図5を参照して説明する。図3の
ステップS10(以下単にS10という。他のステップ
についても同じ)において装置1の作動が開始し、初期
処理が実行される。この初期処理では、特徴情報空間が
クリアされ、カウンタCOUNTの値がゼロにリセット
される。さらに、初期状態を示すためにフラグFFが1
にセットされ、また、カウンタWARM−UP−COU
NTの値が予め定められた値LEARNCOUNTと等
しい値に設定される。LEARNCOUNTの値は20
から500位の値の範囲内であればどのような値であっ
てもよいが、128位であることが望まい。
【0011】次いでS20においてECG信号が取り込
まれる。チャネルコネクタ10、12、14、MUX1
6、フィルタ−バンク18、ADC回路20は共働して
アナログECG信号を取得し、そのアナログ信号をデジ
タル信号に変換し、そのデジタル信号をRAM30に記
憶する。MUX16とI/O回路24は共働して、EC
G信号を250Hz位ないし300Hz位の所定のサン
プリング頻度でサンプリングする。こうして、コネクタ
10、12、14を介して取得されたECG信号にそれ
ぞれ対応するデジタル信号がベクトルe1(t)、e2(t)、
3(t)として取り出され、RAM30に記憶される。
【0012】S30では、記憶されたデジタル信号にQ
RS群候補が含まれているか否かがチェックされる。図
4には、このQRS群候補特定ルーチンの具体的なステ
ップが図示されている。まずS301において、デジタ
ル信号ベクトルがスカラ信号に変換される。スカラ−信
号M(t) は次式を用いて求められる。 m(t)=|e1(t) −e1(t-1)][e2(t)−e2(t-1)][e3(t)−e3(t-1)]| M(t)=|m(t-1)+m(t)+m(t+1)| 但し、M(t)がスカラ−信号であり、t は時間、数字1は
前回ないし次回のサンプリングサイクルで得られるEC
Gサンプルを示すために使用されている。例えば、サン
プリング頻度が250Hzで4ミリ秒(msec)間隔で各
サンプルが取得される場合、(t-1) の記号は、時間(t)
で取得されたサンプルよりも4msec前に取得されたサン
プルを示し、(t+1) の記号は、時間(t) で取得されたサ
ンプルよりも4msec後に取得されるサンプルを示してい
る。スカラ−信号M(t)は次式を使用して得ることもでき
る。 M(t)=|[e1(t+1)−e1(t-2)][e2(t+1)−e2(t-2)][e3(t+1)−e3(t-2)]| 信号処理関係の言葉で表現すれば、スカラ信号M(t)は、
従来の3点(12msec等)ボックスカーフィルタ(3 po
int boxcar filter )から成るフィルタで処理された3
個のECG信号の積の絶対値である。従来の整合フィル
タ(matched filter) 等の他のフィルタを使用して、以
下に詳述するように、デジタル信号中のQRS群の増幅
およびノイズの抑圧を行うことができる。
【0013】続くS302では、スカラ信号M(t)が予め
定められた閾値FIRST−THRESH(通常、最下
位桁の単位(units of the least significant digital
bit)に換算して64)と下側から交差する点を特定す
ることにより、QRS群の候補が検出される。スカラ信
号M(t)が前記所定の閾値と下側から交差した時には、そ
の交差時点を開始点とするウインドウMAX−FWD−
SEARCHの範囲(通常、264msec)内で信号M(t)
の最大値が求めらる。信号M(t)の最大値が精確に決定さ
れているか否かを証明するために、その最大値の位置の
確認を移動ウインドウを使用して行うことができる。こ
うして決定された信号M(t)の最大値がQRS群候補とさ
れ、その発生時間Tqrs は以下のステップでの利用のた
めに記憶される。次いで、図3のS40が実行され、S
30においてQRS群候補が特定されたか否かが判定さ
れる。この判定の結果が否定された場合には、再びS2
0が実行され、肯定された場合は、S50が実行され
る。
【0014】S50では、S301で得られたスカラ信
号から特徴情報が抽出される。図5のS501におい
て、各QRS群候補の周りのスカラ信号が分析されて、
そのQRS群がノイズとされるべきか否かが決定され
る。第i番目のQRS群の場合、Tqrs (i) −wnoise
≦ t ≦Tqrs (i) +wnoise の区間内の信号M(t)が取
り出される。従ってwnoise はノイズ分析に使用される
ウインドウの1/2の幅に相当する。こうして、Tqrs
(i) を中心とする信号M(t)の一部Ni (t) について、こ
れがノイズであるか否かが判断される。この分析におい
ては、信号Ni (t)がノイズ閾値候補と交差する回数が
両方向での交差について数えられる。このノイズ閾値候
補は後に説明するように変化するものであってよい。こ
の交差の回数が予め定められた値MAX−SIGN−C
HANGES(通常、24)に満たない場合にはそのQ
RS群はノイズでないと見なされ、そのノイズ閾値候補
が決定されたノイズ閾値となる。他方、交差の回数が所
定値MAX−SIGN−CHANGESと同数かまたは
それ以上である場合には、ノイズ閾値候補が予め定めら
れた値IJ−INCREASE(通常、最下位桁の単位
に換算して2)分だけ増加させられ、再びノイズ分析が
実行される。もしも、予め定められた回数ノイズ閾値候
補が増加させられたにも関わらず、交差回数が所定値M
AX−SIGN−CHANGESを下回らないうちに、
ノイズ閾値候補が予め定められた最大ノイズ閾値に達し
た場合には、その信号部分Ni (t) はノイズであると見
なされて、今回の分析が終了する。
【0015】前記ノイズ閾値候補は、小さい値から最大
限界値θnoise max まで変化する。その小さい値は、ス
カラ信号M(t)の発生時間Tqrs における評価値の所定割
合に等しい値に設定され、この値は、それ以前の心拍に
ついて決定された前のノイズ閾値との移動平均(θ)
noise ave として平均化された値である。前記最大限界
値は、スカラ信号M(t)の発生時間Tqrs における評価値
を値PEAK−FACTOR(通常、6)で除した値で
ある。これらの記載から、上記ノイズ閾値交差基準を満
たす閾値候補を探す手段として使用され得る他のノイズ
閾値変更技法が当業者には明らかであろう。
【0016】次いで、S502では、スカラ−信号の特
徴情報の1つが抽出される。まず、各QRS群kについ
ての時系列が次式によって求められる。 但し、i=1,2,...,2Wqrs-width であり、W
qrs-width はQRS群の幅分析に使用されるウインドウ
の1/2の間隔を表す。また、F( ) は、スパン約12
msecの従来のボックスカ−フィルタを示し、F5( )は、
そのボックスカ−フィルタが5回連続して適用されるこ
とを表す。前記S501においてQRS群候補がノイズ
であると判断された場合には、QRS群の幅は次式を用
いて決定される。 次いで、幅閾値θqrs-width がその幅閾値の交差が8回
未満となるような最小値に設定される。また、k番目の
QRS群の幅width[k]は、以下の式を用いて決定され
る。 width[k] = Toffset[k] - T onset [k] 但し、T offset[k] およびT onset [k] は、e
width [i] 時系列と、閾値θqrs-width から決定され
る。
【0017】幅width[k]が予め定められた定数である値
MIN−WIDTHより小さい場合、または、値θ
qrs-width が値ewidth (i) のピーク値と値PEAK−
PORTION−FRACTION/PEAK−POR
TION(通常、8分の6)との積よりも大きい場合
は、その心拍は不適当(bad )と見なされ、以下の工程
から除かれる。各リード信号は再びボックスカーフィル
タF()にかけられ、フィルタの出力信号により2次導
関数が得られる。ノイズ心拍を探すための閾値交差基準
が各リードと全てのリードの加算とにそれぞれ適用さ
れ、その心拍候補がノイズと見なされるべきか否か、つ
まり、不適当な心拍として以下の工程から除去されるべ
きか否かについて再度決定がなされる。
【0018】次に、S503において、各QRS群につ
いて複数の面積が決定される。各面積は、QRS群とそ
のQRS群の基線との差を積分することによって求めら
れる。その基線は次式によって求められる。 Bi [k] = (ei [Tonset [k]]+e i [Toffset[k]])/2 但し、i はECG信号のチャネル番号を、kはQRS群
を表す。デジタル信号チャネルについての符号を有する
前半面積は、基線とQRS群との差をQRS群の開始か
らその中央までについて積分することによって求められ
る。ここで、QRS群の中央とは、QRS群の開始(on
set )と終了(offset)の間に位置し、基線とQRS群
の差をQRS群の開始から中央まで積分して得られる符
号を持たない前半面積が、開始から終了点まで基線とQ
RS群の差を積分して得られる符号を持たない全面積の
約1/2になる時点をいう。その無符号全面積および無
符号前半面積は、全てのデジタル信号チャネルについて
の無符号全面積および無符号前半面積をそれぞれ加算す
ることにより決定される。また、QRS群の中央は、そ
の候補値を系統的に探すことによっても決定することが
できる。符号を持たない前半面積は、基線とQRS群と
の差の絶対値を積分することによって求められる。符号
を持つ後半面積と符号を持たない後半面積とが前半面積
と同様に得られる。
【0019】一定の条件下において、心拍は不適当と判
断され、以後の処理から除かれる。生のECG信号値の
QRS群の開始時と終了時の値の差の大きさが各リード
における前高さと後高さの絶対値の合計以上である場合
には、心拍は不適当と判断される。その前高さおよび後
高さの値は、処理される各リードについて2つずつの計
算された値を含む。通常、3本のリードが処理されるの
で、3個の前高さと3個の後高さが求められる。ここ
で、前高さとは、任意のリードについて、前記前半面積
の最大振幅(正である場合も負である場合もある)をい
う。後高さとは、任意のリードについて、前記後半面積
の最大振幅をいう。また、QRS群の幅決定のために選
択される最後の閾値での交差回数がゼロである場合に、
心拍は不適当と判断される。さらに、その幅が値MIN
−WIDTHよりも小さい場合に、心拍は不適当と判断
される。これまでに定義されていないパラメータについ
ては、以下において定義される。次いで、プログラムは
S60に進む。
【0020】図3のS60では、スカラ信号から抽出さ
れた特徴情報が、N次元特徴情報空間にプロットされ
る。次元Nは抽出された特徴情報の数に一致する。これ
によって、QRS群のクラスタ−が形成される。形成さ
れたクラスタ−は、後述するように、本ステップにおい
て更新される。ECG信号から3個のチャネルが受信、
補間ないし補外によって取得され、QRS群の幅に加え
て、各チャネルについて各々4個の面積が決定されるこ
とから、合計13個の次元を有する特徴情報空間が形成
される。もっとも、個々のQRS群の特徴情報自体は、
それらQRS群が正常であるか否かの指標とはならな
い。その判断のために、他のQRS群と比較される。す
なわち、各QRS群は、特徴情報空間中のクラスタに含
まれるQRS群と形状および発生時間について比較され
る。
【0021】各QRS群から抽出された特徴情報は、以
下に説明するクラスタリングまたはグルーピングと呼ば
れる手法によりそれ以前のQRS群から抽出された特徴
情報と比較される。RAM30には既存のクラスタのリ
ストが記憶されている。このリストには、特徴情報空間
中の各クラスタの位置についての情報(例えば、そのク
ラスタに含まれるQRS群の特徴情報の平均値)、その
クラスタ内に位置するQRS群の通常の前方および後方
連結間隔(例えば、その平均値)、そのクラスタ内に位
置するQRS群の数、および、そのクラスタに含まれる
QRS群のラベルないしグループ名が含まれる。ここ
で、前方および後方連結間隔とは、各QRS群とその前
後のQRS群との間の各時間間隔をいう。前記リストは
さらに、各クラスタに位置付けられた最初のQRS群と
最新のQRS群の時間Tqrs 、各クラスタに2個のQR
S群が連続して加えられた回数である連続付加回数、お
よび、UNKNOWN、NORMAL、PVC、DEL
ETE等のクラスタの種類を表すグループ名を含む。E
CG記録の分析の開始時にはクラスタは1個も存在しな
い。クラスタは、以下の態様で形成される。
【0022】まず、検出された各QRS群候補につい
て、そのQRS群候補の特徴情報が、クラスタリスト中
の各クラスタの特徴情報と比較される。このリスト中に
クラスタが全く存在しない場合は、新たにクラスタが作
られ、このクラスタにはそのQRS群候補の特徴情報に
対応する前記空間内の位置とそのQRS群候補の前方お
よび後方連結間隔が付与され、そのクラスタに位置付け
られる最初と最後のQRS群の発生時間がそのQRS群
候補の発生時間とされる。また、そのクラスタに含まれ
るQRS群の数が1とされ、グループ名UNKNOWN
が付与される。そして、この新しいクラスタがリストに
加えられる。
【0023】リストにクラスタが存在する場合には、前
記特徴情報空間におけるそのQRS群の特徴情報と各ク
ラスタの位置の間の距離が求められ、そのQRS群に最
も近接するクラスタが判別される。この最近のクラスタ
との距離が値MAX−DXより小さい場合には、そのQ
RS群は最近のクラスタに加えられる。このようにし
て、クラスタに含まれるQRS群の数、つまり、個数が
増加し、また、特徴情報空間内におけるクラスタの以前
の位置とそのQRS群の特徴情報との平均値を求めるこ
とによってその情報空間内のクラスタの位置を更新し
(例えば、クラスタの以前の平均位置には8分の7の重
みが与え、新しいQRS群の特徴情報には8分の1の重
みを与える)、同様に、前方および後方連結間隔を更新
し、連続付加回数を増加させる。そのQRS群に最も近
いクラスタが値MAX−DX(通常、400)以上であ
る場合には、前記のように、新たにクラスタが形成さ
れ、クラスタリストに加えられる。
【0024】クラスタが次々に形成されるに従い、クラ
スタリストは大きくなる。リストの長さが所定の長さに
達すると、新たなクラスタが加えられる度に、古いクラ
スタがリストから削除され、それによりリストの長さは
一定に保たれる。リスト可能なクラスタの数を最高16
位とすれば、分析能率を低下させることなく精度を向上
させることができる。削除されるクラスタ−は、QRS
群が加えられた時間がもっとも古く、かつ、リストされ
ているクラスタ−の中でQRS群の個数の少ない方の3
分の1に含まれるものである。続いて、S70が実行さ
れ、カウンタCOUNTの値が1だけに増加される。
【0025】S80において、フラグFFが0であるか
否かが判定される。YESなら後述するステップS90
が実行される。NOならS110が実行され、本装置が
作動を開始してから充分な数、すなはち、予め定められ
た数WARM−UP−COUNTのQRS群候補が集め
られて特徴情報空間中にプロットされたか否かが判断さ
れる。当初、QRS群は、特徴情報空間においてクラス
タの形成のために使用され、クラスタや個々のQRS群
に適当なグループ名が付与されることはない。。もし、
カウンタCOUNTの値が値WARM−UP−COUN
Tより小さい場合には、S20が実行され、前記のよう
に、後続のQRS群が集められる。特徴情報空間中にプ
ロットされたQRS群の総数が、値WARM−UP−C
OUNT以上になると、S110の判断が肯定され、続
くS120において各クラスタを分析しそれが正常なQ
RS群に対応するものであるか否かの判定がなされる。
【0026】S120においては、まず、最高個数のQ
RS群を含むクラスタが併合処理のために選ばれる。近
くの他のクラスタが類似の特徴情報を有することがある
ため、最大のクラスタと他のクラスタとの距離が決定さ
れる。決定された距離が予め定められた値MAX−DX
未満であれば、それら2つのクラスタは、それぞれの平
均値の平均値を取り、それぞれの有する最大値の最大値
を取り、それぞれの最小値の最小値を取り、それら2つ
のクラスタへの最後のQRS群の付加を新しいクラスタ
への最後のQRS群の付加とすることにより、併合され
る。
【0027】最大のクラスタに含まれるQRS群の個数
が他のクラスタとの比較において極めて大きい場合、例
えば、その個数が予め定められた数MIN−INI−N
SR−POP(通常、数WARM−UP−COUNTの
2分の1)より大きい場合は、そのクラスタは正常と判
別されて正常のグループ名が付与される。他方、正常心
拍律動状態が肯定される。最大のクラスタに含まれるQ
RS群の個数が数MIN−INI−NSR−POP未満
の場合には、二段脈である可能性がある。二段脈は、最
大クラスタと二番目に多くのQRS群を含むクラスタと
を判別し、各クラスタについて後方連結間隔の平均値を
決定することにより行われる。ここで、後方連結間隔と
は、今回のQRS群と前回のQRS群の間の時間間隔を
いう。2つの後方連結間隔の中の短い方の後方連結間隔
が長い方の後方連結間隔の0.75倍未満であり、か
つ、その2つの後方連結間隔の平均値(短い方と長い方
の両方の後方連結間隔の平均値)の8分の7未満であ
り、かつ、両クラスタの個数が予め定められた数MAX
−INI−BGM−TOG(通常、数WARM−UP−
COUNTの8分の1)未満であり、かつ、両クラスタ
の個数の合計が予め定められた数MIN−INI−BG
M−TOG未満である場合にも、二段脈の可能性があ
る。この場合、長い方の後方連結間隔を有するクラスタ
が正常とされ、他のクラスタはPVC心拍とされる。ま
た、二段脈状態が肯定され、正常心拍律動状態が否定さ
れる。クラスタをグループ化するための上記条件が満た
されない場合には、(i)全クラスタ−の総計個数に等
しい数POP−SUMが数WARM─UP─COUNT
より小さいか、(ii)最大のクラスタの個数が2番目
に大きなクラスタの個数の2倍より大きいか、(ii
i)最大のクラスタの個数が数POP−SUMの半数よ
り大きければ、その最大クラスタが正常であると判別さ
れる。この場合には、二段脈状態が否定され、正常心拍
律動状態が肯定される。S130の結果、正常なクラス
タが見つからなかった場合は、S150において値WA
RM−UP−COUNTに値LEARN−COUNTが
加えられた後、再びステップS20が実行され、新たに
値WARM−UP−COUNTの数のQRS群が集めら
れる。S130が肯定されればS140においてフラグ
FFがゼロにセットされ、次いでS20が実行される。
【0028】S140においてフラグFFがゼロとされ
ることにより、クラスタ“学習”が終了すると、図示し
ない複数の付加的一回初期化ステップが実行される。二
段脈がある場合を除き、正常クラスタを除いた全てのク
ラスタのQRS群の個数が1とされる。二段脈がある場
合には、正常クラスタとPVCクラスタを除いた全ての
クラスタのQRS群の個数が1とされる。正常なQRS
群のパラメータが確定される。すなはち、平均的な正常
QRS群の幅、符号付き面積、平均NN時間間隔および
平均RR時間間隔が確定される。平均的正常QRS群の
幅としては、正常クラスタの平均幅が使用される。平均
的正常QRS群の符号付き面積としては、正常クラスタ
の平均面積が使用される。平均的正常QRS群の平均N
N時間間隔は、二段脈がない場合には、正常クラスタの
平均後方連結間隔と平均前方連結間隔との平均値であ
り、他方、二段脈がある場合には、長い方のクラスタの
平均後方連結間隔と短い方のクラスタの平均後方連結間
隔との平均値である。平均的正常QRS群の平均RR時
間間隔は、平均NN時間間隔は、R波間隔であって上記
平均NN時間間隔と等しい。クラスタ“学習”の後、新
しいQRS群候補が検出されると、それらはクラスタと
されグループ名を付与される。未だグループ名を与えら
れていないクラスタに含まれる新しいQRS群の数が大
きくなってそのクラスタに信頼性を以てグループ名を付
与することができる充分な数に達する場合があり得る。
グループ名は複数あり、例えば、正常(NORMA
L)、心室性期外収縮(PVC)、削除(DELETE
D)を含む。未だグループ名を与えられていないクラス
タに属する1個のQRS群が正常とグループ化される
と、後に説明するルールによってそのクラスタ全体が正
常とグループ化される場合がある。もっともそれが行わ
れるのは、そのクラスタのQRS群の総数が値MIN−
NEW−NORMAL−POP(通常、10)より大き
く、かつ、そのクラスタに含まれる正常なQRS群の個
数がそのクラスタのQRS群の総数の半分より大きいか
もしくはそのクラスタに含まれる正常なQRS群の個数
がそのクラスタに含まれる異常なQRS群の個数の4倍
より大きいかの何れかである場合のみである。未だグル
ープ名を与えられていないクラスタに属する1個のQR
S群が削除(DELETED)とグループ化されると、
後に説明するルールによってそのクラスタ全体が“削
除”とグループ化される場合がある。もっともそれが行
われるのは、そのクラスタのQRS群の総数が値MIN
−NEW−DELETE−POP(通常、6)より大き
く、かつ、そのクラスタに含まれる“削除”QRS群の
個数がそのクラスタのQRS群の総数の半分より大きい
場合のみである。クラスタ学習の後に、患者の体動等の
外部的作用によって新しいQRS群が検出されることが
ある。この場合、患者が極めて乱雑な律動を有するとき
にデータに含まれるノイズ部のために新しく多数のクラ
スタが形成され、それらによって既存のクラスタが置換
される。この状態においては、各クラスタに含まれるQ
RS群の個数が1とされ、クラスタに付与されたクラス
タ名が消去される。その後、また新しいQRS群候補が
検出されるに従ってクラスタにグループ名を付与する作
動が行われる。クラスタ学習の後に新しいQRS群候補
が検出されてグループ名を付与されるに従い、平均NN
時間,平均NW時間,平均NSS時間,平均RR時間を
含む全体的律動/形状情報が更新される。心拍が正常と
グループ化される毎に、平均NNが更新される。平均N
N時間の値は、(7×平均NN+RR)/8に更新され
る。この式中、平均NN時間は平均NN時間の今回の値
であり、RR時間は正常心拍間の時間である。また、心
拍が正常とグループ化される毎に、平均NW時間および
平均NSS時間が更新される。平均NW時間の値は、
(7×平均NW+幅)/8に更新される。この式中、幅
は、今回の正常心拍の測定された幅である。平均NSS
の値は、(7×平均NSS+SAS)/8に更新され
る。式中、SASは、今回の正常心拍の符号を持つ面積
の総計、すなわち、すべてのリードについての前半面積
と後半面積との合計である。さらに、心拍が検出される
と、その心拍が正常であるかPVC心拍であるかに関わ
らず、平均RR時間が(7×平均RR+rr)/8に更
新される。この式中、rrは今回の心拍と前回の心拍の
間の時間である。患者の律動状態は、(1)相当に規則
的な連結間隔を持った正常なものであるか、(2)連結
間隔が変動する不規則なものであるか、(3)連結間隔
が変動してかなりのクラスタが新しく形成される稀(un
usual )なものであるか、(4)連結間隔が大きく変動
して多数のクラスタが新しく形成される乱雑(bizarre
)なものであるかの何れかに評価される。連結間隔の
変動度はNN間隔でなくRR間隔に基づいて決定され
る。個々の患者の状態を判断するための基準が満たされ
るとその状態が肯定され、それ以外の場合には否定され
る。各状態は“正常”から“乱雑”の方向に順位が付さ
れており、“乱雑”が第1順位を有する。従って“乱
雑”のための基準が満たされれば、“稀”であるかどう
かに関わらず、“乱雑”状態が宣言される。ある状態で
あることが認定される毎に、この認定が出力として示さ
れ、回数定数RHYTHM−COUNTがセットされ
る。各状態の最初の検出の際に、その状態がこの値にセ
ットされる。心拍が新たに処理される毎に、前記状態カ
ウントが減少する。各状態のための状態カウントが減少
してゼロになり、さらに他の状態のための基準が満たさ
れると、新たな出力指示が出される。グループ名を付与
された心拍のみ(グループ名の付与が先に延ばされたも
のは含まない)が処理されて、患者の律動状態の決定に
使用される。心拍の処理のためにクラスタリストが調べ
られて、クラスタの形成速度が決定される。クラスタの
形成速度は、クラスタリストに含まれるクラスタの数で
あって数RECENT−SHAPES(通常、32)よ
り小さいものを所定の期間内に数えることにより決定さ
れる。患者が普通でない律動状態(“乱雑”ないし
“稀”)にあって、普通でない律動のための基準(すな
わち、多くの新しいクラスタが形成され、間隔が不規
則)がもはや満たされなくなったならば、律動状態は正
常律動状態に変更され、再学習が開始される。“乱雑”
な状態と“稀”な状態の何れかが否定されず、カウント
RECENT−SHAPESが値shape−coun
t−LO−THRESHOLD(通常、5)より小さ
く、かつ、カウントRR−differenceが値R
Rdif−LO−THRESHOLD(通常、5)より
小さい場合には、再学習が適当とされる。再学習が適当
とされ、かつ、“乱雑”な状態が肯定されれば、再学習
が開始され、クラスタのQRS群の個数が1とされ、
“正常”,“異常”,“削除”カウントがゼロとされ
る。その後、“正常”状態が値RHYTHM−COUN
T(通常、100)とされ、“不規則”状態が否定さ
れ、“稀”状態が否定され、“乱雑”状態が否定され
る。第2順位は“乱雑”状態のためのテストである。
“乱雑”状態が否定され、カウントRECENT−SH
APEが値shape−count−HI−THRES
HOLD(通常、10)以上であり、かつ、カウントR
R−differenceが値RRdif−count
−HI−THRESHOLD(通常、9)以上である場
合には、“乱雑”形状が値RHYTHM−COUNTに
設定され、“正常”状態が否定され、“不規則”状態が
否定され、“稀”状態が否定される。その他、“稀”状
態または“乱雑”状態が否定され、カウントRECEN
T−SHAPEが値shape−count−LO−T
HRESHOLDより大きく、かつ、カウントRR−d
ifferenceが値RRdif−count−LO
−THRESHOLD以上である場合には、“稀”状態
が値RHYTHM−COUNTに設定され、“正常”状
態が否定され、“不規則”状態が否定され、“乱雑”状
態が否定される。さらに、“不規則”状態が否定され、
カウントRR−differenceが値RRdif−
count−LO−THRESHOLD以上であり、か
つ、カウントRECENT−SHAPEが値shape
−count−LO−THRESHOLDより小さいな
らば、“不規則”状態と再学習が設定され得る。“不規
則”状態と再学習が設定され得る状態にあって、かつ、
“乱雑”状態が設定されれば、再学習が開始され、クラ
スタのQRS群の個数が1とされ、“正常”,“異
常”,“削除”カウントがゼロとされる。また、“不規
則”状態が値RHYTHM−COUNTに設定され、
“正常”状態が否定され、“稀”状態が否定され、“乱
雑”状態が否定される。上記の基準テストがグループ化
後の有効な心拍の各々について実施される。そのテスト
の実施後、“不規則”状態カウントはゼロより大きけれ
ば減少し、“稀”状態カウントはゼロより大きければ減
少し、“乱雑”状態カウントはゼロより大きければ減少
する。最後に、“乱雑”状態または“稀”状態が肯定さ
れると、非正常状態(unnormal)が肯定される。クラス
タ学習の後に新しいQRS群候補(心拍)が検出される
際に、各心拍の分析の前に、律動および形状の限界が確
定される。パラメータearly,lateがそれぞれ
7×平均NN/8と3×平均NN/2と定義される。ま
た、パラメータearly−fusion,late−
fusionがそれぞれ7×平均NN/8と9×平均N
N/8と定義される。さらに、パラメータlow−c−
pause,hi−c−pauseがそれぞれ29×平
均NN/16と(平均NN+late)と定義される。
パラメータwide,narrowがそれぞれ5×平均
NW/4と5×平均NW/8と定義される。パラメータ
Twv−RRは、T波(T-wave)が発生する最大時間で
あり、通常、350msecに設定される。パラメータ
PZ−DURは休止(pause )と認定するために必要な
心拍間の時間長さであり、通常、2.5secに設定さ
れる。パラメータNN−TOL,NNZ−TOLはそれ
ぞれ非ノイズ心拍とノイズ心拍についてRR間隔を平均
NNと比較するための許容誤差であって、通常、48s
ecと96secとに設定される。新しいQRS群候補
が検出される時、その中の幾つかのQRS群候補が不適
当と判断されて以後の処理から除かれる場合がある。各
心拍についてのRR間隔が値earlyTより小さく、
アサーションLowHeight()が肯定され、ノイ
ズ検出の際に決定された符号変化の回数が値SIGNC
H−NOBEAT(通常、5)より大きく、幅決定のた
めに使用されたウインドウ内の信号の2次導関数におけ
る符号の変化回数が値SUM123−NOBEAT(通
常、13)より大きく、かつ、すべてのリードについて
の合計の無符号面積が値SUMU−NOBEAT(通
常、12500)より大きい場合には、その心拍が不適
当と認定される。面積はある時間間隔についての複数サ
ンプルの合計であり、サンプリング率は通常250Hz
である。リードから得られるアナログ信号は増幅された
後にデジタル信号に変換され、その増幅ゲインは自動ゲ
イン制御によって調整される。また、アサーションLo
wHeight()が肯定され、幅が値WIDTH2−
NOBEAT(通常、160msec)より大きい場合
に、心拍は不適当と認定される。さらに、RR間隔が2
00msecより小さく、かつ、前回の心拍が正常であ
ったならば、その心拍は不適当と認定される。また、前
回の心拍が正常であり、RR間隔が値earlyTより
小さく、かつ、幅がパラメータnarrow+NW−T
OL(通常、16msec)より小さく、アサーション
LowHeight()が肯定され、ノイズ検出の際に
決定された符号変化の回数が値SIGNCH−NOBE
AT(通常、5)より大きい場合に、心拍は不適当と認
定される。パラメータNWZ−TOLは通常、パラメー
タNW−TOLの2倍である。図3のS90において、
以下の表1乃至表4に示されるルール等の複数のル−ル
が、デジタル信号に含まれるQRS群候補に適用され
て、各QRS群をグループ化され、正常なQRS群、異
常なQRS群(すなわち、正常または削除QRS群以外
のもの)、削除QRS群等のグループ名が付与される。
【0029】
【表1】
【表2】
【表3】
【表4】
【0030】これらのル−ルは、すでにグループ名を付
与されたクラスタ−内に含まれるQRS群に優先的にそ
のクラスタ−のグループ名を与えられるように重みが付
けられている。こうして、各QRS群は受信の際にグル
ープ名が与えられる。しかし、直ぐにはQRS群のグル
ープ化を行うことができず、後続のQRS群を受け取る
までQRS群のグループ化を遅らせて待機状態とする場
合がある。待機させられ得るQRS群の最大数は5個等
の予め定められた数であり、その場合に、その後に供給
されるQRS群を加えて未クループ化QRS群は最大で
6個となる。6個以上のQRS群が待機させられる場合
には、最も古い待機QRS群が、削除QRS群と同様の
効果をもって削除される。合計6個の未クループ化QR
S群が存在する場合には、その6個のQRS群は以下の
ようにしてグループ化される。表1乃至表4のグループ
化ルールは待機中のQRS群の数に従い表中の順序で適
用される。例えば、1拍の場合には待機中の心拍はな
く、2拍の場合には今回の心拍と待機中の1個の心拍と
を含む。表1乃至表4の左側部分の各項はアサーション
(assertion )であり、()内の数字はそのアサーショ
ンが適用されるQRS群を示している。最も小さい数字
で示されるQRS群は最も古いQRS群である。例え
ば、アサーションNormL(1) は、最も古いQRS群
が正常クラスタ−に含まれるか否かを問い、アサーショ
ンNormL(5) は、第5番目の待機QRS群が正常ク
ラスタ−に含まれるか否かを問う。さらに、アサーショ
ンNormT(3-2) は、第3番目と第2番目の待機QR
S群の発生時間の差が正常範囲内であるか否かを問う。
また、NormT(2) のような表現はNormT(2-0)
を簡略形であり、ここで第ゼロ番目のQRS群とは今問
題にしている心拍列の直前に発生したQRS群である。
表1乃至表4において、記号“&&“は論理積(AN
D)を、記号“||“は論理和(OR)を表す。アサー
ションは、QRS群を含むクラスタ、QRS群のノイズ
成分、QRS群の形状、QRS群の発生時間、さらに、
デジタル信号の全体的律動を扱う。
【0031】表1乃至表4の複数のルールに従ってグル
ープ化される各心拍は、正常、T波、削除、異常、また
は、未グループ化の一つに類別される。表1乃至表4の
右側部分のグループ化指示は種々のグループ名に従って
グループ化し、待機、正常、T波、削除、未グループ化
の何れでもないグループ化指示はすべて、異常に類別さ
れる。例えば、グループSVPCWP(休止(pause)付
き上室性期外収縮)は異常に類別される。また、グルー
プSVPC(上室性期外収縮)、VentESC(心室
性逸脱心拍)、PVC(心室性期外収縮)、SupV
(super ventricular )、フュージョン(fusion)、フ
ァニー(funny )はすべて異常心拍に類別される。記号
SetQ()はグループ化指示ではなく、後に説明され
る関数である。
【0032】クラスタの位置付けに関係するアサーショ
ンは、アサーションNormL( )、NoL( ) 、PV
C( ) およびDelL( ) を含み、それらのアサーショ
ンはそれぞれ、QRS群が正常、未グループ化、PVC
または削除にグループ化されたクラスタに位置付けられ
ているか否かを問う。また、アサーションSameCl
uster(a,b)は、QRS群a,bが同じクラス
タに属するか否かを問い、アサーションSmallPo
p( ) は、QRS群を含むクラスタの個数が予め定めら
れた値MIN−NEW−NORMAL−POP未満であ
るか否かを問う。ノイズ成分に関係するアサーション
は、アサーションNoisy( ) ,NotNoisy
( ) を含み、既に詳しく説明したように、QRS群を含
む部分がノイズ性と判断されるか否かを問う。
【0033】QRS群の形状を扱うアサーションは、ア
サーションWideW( ) 、NormW( ) 、N2W
( ) およびNNoiseW( ) を含む。アサーションW
ideW( )は、QRS群の幅が(5×平均NW)/
4と定義される値Wideより大きい場合に肯定され
る。アサーションNormW( ) は、QRS群の幅が値
wideと値narrowの間にある場合に肯定され
る。但し、値narrowは最初(5×平均NW)/8
に設定され、平均NWは正常心拍が2個連続して取り込
まれた後に(7×平均NW+幅)/8に更新される。ア
サーションN2W( )とNNoiseW( ) は、QRS
群の幅がそれぞれ平均NW±NW−TOLの範囲内と平
均NW±NWZ−TOLの範囲内にある場合に肯定され
る。値NW−TOLと値NWZ−TOLはそれぞれ予め
定められた許容誤差値(通常、16msecと32ms
ec)である。アサーションLowHight( ) は、
QRS群の高さが予め定められた値MIN−HIGHT
(通常、最小A−Dステップ高さの300倍)未満であ
る時に肯定される。また、アサーションDiffAre
aSign(a,b)は、何れかのリードにおいて、Q
RS群a,bのそれぞれの符号付き前半面積ないし後半
面積の符号が異なる場合に肯定される。
【0034】アサーションAbnormalRhyth
m( ) は、まれにしか問われないアサーションであっ
て、患者が不正常な律動状態にある時(小さなクラスタ
が多数新たに形成される時)に肯定される。例えば、患
者の心臓律動状態が“乱雑”または“稀”な場合に肯定
される。発生時間に関係するアサーションは、QRS群
のRR間隔を取り扱う。図6および図7には、時間を取
り扱うアサーションを評価するために使用される平均N
Nを中心として複数の区間が表されている。図中、ライ
ンの下側の記号は平均NNに基づく端点値を示す。ライ
ンの上側の区間は、ラインの下側の値から導かれる。例
えば、アサーションNormT( ) は、RR間隔が値e
arly(7/8 ×平均NN)と値late(3/2 ×平均
NN)の間にあれば肯定される。アサーションEarl
yT( ) およびLateT( ) は、RR間隔が値ear
ly未満であり、値lateより大きい場合にそれぞれ
肯定される。アサーションGoodT( ) は、RR間隔
が値earlyとlateの間にある場合(発生時間が
良好である場合)に肯定される。アサーションFusi
onT( ) は、RR間隔が値early−fusion
と値late−fusionの間にある時に肯定され、
アサーションOKT( ) は、RR間隔が値earlyと
late−fusionの間にある場合に肯定され、ア
サーションCpause( ) は、RR間隔が値low−
c−pauseと値hi−c−pauseの間にある場
合に肯定される。また、アサーションLT−c−Pau
se( ) は、RR間隔が値low−c−pause未満
である場合に肯定される。アサーションRunT( )
は、RR間隔が平均NN+early未満である場合に
肯定される。また、幾つかの発生時間アサーションが、
RR間隔が平均NNの予め定められた許容誤差内に入っ
ているか否かを扱う。例えばアサーションNNT( ) や
ZNNT( ) がそれで、それぞれ、RR間隔が平均NN
±NNTOL(値NNTOLは約48msec)の範囲
内と平均NN±NNZTOL(値NNZTOLは約96
msec)の範囲内とにあれば肯定される。さらに、ア
サーションSlateT( ) およびSearlyT( )
は、RR間隔が平均NN±(NNTOL/4)の範囲内
であれば肯定される。最後に、アサーションLongP
auseT( ) は、RR間隔が予め定められた値PZ−
DUR(例えば、約2.5msec)より大きい場合に
肯定され、アサーションTwaveT( ) は、RR間隔
が予め定められた値Twv−RR(例えば、約350m
sec)未満である場合に肯定される。
【0035】図3のS90において、待機中のQRS群
の数が決定され、その数に適したルールが、グループ化
を必要とするQRS群に適用される。但し、6拍の場合
は例外で、それを以下に説明する。表1乃至表4の左側
部分の各ラインに示される全てのアサーションが肯定さ
れると、同じラインの右側部分に示されているグループ
名がそのQRS群に付与される。例えば、3拍の場合に
アサーションOKT(3)が肯定されると、最も新しい
QRS群は正常なQRS群、残る2つのQRS群は削除
QRS群とされる。3拍の場合に最も新しいQRS群と
は第3番目の心拍である。右側部分に示されるD( ) 、
N( ) 、PVC( ) 、F( ) 等のグループ化指示記号
は、第( ) 番目のQRS群がそれぞれ削除、正常、PV
C、フュージョンのグループ名を付与されるべきことを
示す。D( ) 等のグループ化指示記号は、当QRS群を
含むクラスタ−が以下に説明するクラスタ更新ルールに
従ってグループ化されるべきことを示す。
【0036】記号SetQ( ) は、心拍が類別されず、
その心拍の前後の心拍も類別されていない可能性がある
場合に、グループ化を指示するための一群の付加的ルー
ルである。記号SetQ(p) で示されるルールは以下の
定義を用いて理解される。すなわち、値RRは、心拍p
と心拍p−1の間の時間間隔である。また、値CRR
は、既に説明した心拍中央について補正された時間RR
である。さらに、値RR2は、心拍p−1と心拍p+1
の間の時間間隔である。アサーションNotNoisy
(p)&&NotNoL(p)が肯定されれば、心拍p
はそれを含むクラスタのグループ名を取得する。それ以
外の場合に、次のアサーション (5×平均NN/8<RR<5×平均NN/4)&&
(NNoiseW(p))&&(平均NN−NNZTO
L<CRR<平均NN+NNZTOL)&&(平均NW
−NNZTOL<心拍pの幅<平均NW+NNZTO
L) が肯定されれば、心拍pは正常とされる。それ以外の場
合に、次のアサーション (7×平均NN/4<RR2<5×平均NN/2)&&
(RR<7×平均NN/8)&&Not(平均NW−T
OL<心拍pの幅<平均NW+TOL) が肯定されれば、心拍pはPVCとされる。但し、値T
OLは前記サンプリング頻度の場合に5である。それ以
外の場合に、不正常状態が否定されれば、発生時間が極
めて早いためにT波とみなすことができるかどうかの決
定がなされる。今回の心拍候補と前回の心拍の間の時間
間隔RRが350msec未満であるならば、その心拍
候補はT波とされる。そしてT波はらば、“削除”のグ
ループ名が付与される。それ以外の場合に、 NotAbnormalRhythm が肯定されれば、心拍pは“疑いあり”のグループ名が
付与される。それ以外の場合に、 RR<7×平均NN/8 が肯定されれば、心拍がPVC心拍であるかのように心
拍のクラスタが更新される。それ以外の場合に、 RR<3×平均NN/2 が肯定されれば、心拍が正常な心拍であるかのように心
拍のクラスタを更新される。それ以外の場合には、心拍
がPVC心拍として再グループ化され、それに従って心
拍のクラスタが更新される。グループ名を付与されたク
ラスタ−が、そのグループ名と異なるグループ名を有す
るQRS群を含むことがある。
【0037】6個のQRS群のグループ化を行う場合に
は、それらQRS群のグループ化を反復して行うために
付加的ルールが適用される。最初に、チェックNorm
L(6) &&NormL(5) &&GoodT(6-5) が適用
され、これが肯定されると、それら6個の心拍にそれぞ
れSetQ(1) 、SetQ(2) 、SetQ(3) 、Set
Q(4) 、N(5) 、N(6) のグループ名を付与するように
指示される。他方、このチェックが否定された場合に
は、最新のQRS群から順をおってその心拍が既にグル
ープ名を付与さえたクラスタ内に位置付けられているか
否かが確認される。この作業は反復して行われる。今回
問題とされる心拍を心拍nとする。各QRS群(心拍
n)についてチェックNormL(n) が肯定されるか否
かが判定される。肯定される場合には、それら心拍には
グループ名N(n)が付与されて、次の心拍の処理に進
む。チェックNormL(n) が否定された場合には、そ
れら心拍が実際はECG信号のT波であるか否かの判定
がなされる。チェックNotGoodT(n-(n-1)) (但
しnは今回のQRS群、n-1 は前回のQRS群を表す)
およびチェックLowHeight(n) の双方が肯定さ
れれば、第n番目のQRS群はT波であり、“削除”に
グループ化される。
【0038】最後に、QRS群が正常の範囲内にあるか
否かの確認がなされる。(i)チェックNoisy(n)
が肯定されかつチェックPVCL(n) が否定されるか、
もしくは、第n番目のQRS群の幅が値narrowと
値(wide+NW−TOL)の間の値である場合(但
し、NW−TOLは予め定められた許容誤差値であ
る)、または、(ii)すべてのリードについての符号
を有する面積の総和が(3×平均NSS)/4と(5×
平均NSS)/4の間の値である場合(但し、平均NS
Sは正常クラスタに含まれるQRS群の符号付き前半面
積と符号付き後半面積の総和の平均値である)に、それ
らQRS群は正常の範囲内にあるとされる。心拍が正常
の範囲内にあるとされた場合には、それら心拍を“正
常”とグループ化してよいか否かが決定される。QRS
群nとQRS群mのそれぞれのRR間隔の差が値(ea
rly−NNTOL)より大きなQRS群mが探され
(但し、値NN−TOLは予め定められた値である)、
もしもQRS群nとQRS群mが隣接しておりそのため
m=n−1である場合には、グループ名N(n) が付与さ
れる。しかし、心拍m、nの間にQRS群候補が存在す
る場合には、それらの中間QRS群に対して以下のル−
ルが適用される。チェックNoisy( ) またはLow
Height( ) が肯定されれば、上記QRS群にはグ
ループ名D( ) とされる。それ以外の場合に、チェック
AbnormalRhythmが肯定されれば、グルー
プ名PVC( ) が付与される。それ以外の場合に、心拍
の属するクラスタがグループ名を有していればそれら心
拍にはそのクラスタのグループ名と同じグループ名が付
与される。それ以外の場合には、QRS群には“未決定
(suspicious)”とのグループ名が与えられる。以上の
ようにして最新の心拍から最も古い心拍まで6個の心拍
が反復的に処理されると、グループ化作業の第2段階に
入り、その6個の心拍が再び反復的に作業の対象とされ
る。今回の処理は最も古い心拍から始まって最新の心拍
で終了する。今回処理される心拍を心拍pとする。今回
の心拍が既にグループ名を付与されている場合には、今
回の心拍の処理を省略して次の心拍の処理に取りかか
る。今回の心拍pがp<6であり、チェックNotNo
rmW(p)が肯定され、チェックEarlyT(p)
が肯定され、チェックNotNoisy(p)が肯定さ
れ、かつ、心拍p,p+1の間の間隔が値early−
fusionより大きい場合には、グループ名PVC
(p)が与えられる。心拍pがp<6であり、かつ、チ
ェックNoL(p+1)が肯定された場合には、グルー
プ名N(p)が与えられる。それ以外の場合で、チェッ
クNoL(p+1)が肯定され、チェックNotNor
mW(p)が肯定され、チェックNotNormW(p
+1)が肯定され、チェックTwaveT((p+
1),p)が肯定され、かつ、チェックDiffAre
aSign(p,(p+1))が肯定されれば、グルー
プ名PVC(p)およびD(p+1)が与えられる。そ
れ以外の場合で、チェックNoisy(p)が肯定され
るか、または、チェックLowHeight(p)が肯
定されると、グループ名D(p)が与えられる。それ以
外の場合で、チェックAbnormalRhythmが
肯定されると、グループ名PVC(p)が与えられる。
それ以外の場合で、チェックNoL(p)が否定される
と、その心拍を含むクラスタのグループ名がその心拍に
与えられる。それ以外の場合で、その心拍を含むクラス
タの正常心拍の個数もしくは異常心拍の個数が値MIN
−NEW−NORMAL/2よりも大きい場合に、異常
心拍の個数が正常心拍の個数より大きければそれら心拍
を“異常”とグループ化し、異常心拍の個数が正常心拍
の個数以下であればそれら心拍を“正常”とグループ化
する。それ以外の場合には、心拍を“未決定”とグルー
プ化する。ついで作動はS100に進む。心拍グループ
化指示が決定されると、その後、心拍の属するクラスタ
がその心拍のグループ名の類別に従って更新される場合
がある。すなわち、心拍が正常とされた場合にはクラス
タも“正常”とされ、PVCを含むその他すべての真の
QRS群については“異常”とされ、T波ないしノイズ
等の偽のQRS群については“削除”とされる。正常心
拍に従ってクラスタのグループ名を更新する際には、ク
ラスタの正常心拍の個数が増加され、そのクラスタを
“正常”とすることが可能か否かのテストが適用され
る。このテストは、後に説明するようにクラスタにグル
ープ名が付与されていないものとされる場合に適用され
るテストと同じものである。異常心拍に従ってクラスタ
のグループ名を更新する際には、クラスタの異常心拍の
個数が増加され、そのクラスタを“異常”とすることが
可能か否かのテストが適用される。このテストは、後に
説明するようにクラスタにグループ名が付与されていな
いものとされる場合に適用されるテストと同じものであ
る。削除心拍に従ってクラスタのグループ名を更新する
際には、クラスタの削除心拍の個数が増加され、そのク
ラスタを“削除”とすることが可能か否かのテストが適
用される。このテストは、後に説明するようにクラスタ
にグループ名が付与されていないものとされる場合に適
用されるテストと同じものである。
【0039】S100において、1個ないし複数個のQ
RS群のグループ化の後に、前記クラスタ更新ルールに
従ってクラスタのグループ化が行われる。あるQRS群
が正常とされた場合において、クラスタに含まれる全Q
RS群の個数が予め定められた値MIN−NEW−NO
RMAL−POPより大きく、かつ、そのクラスタに含
まれる正常QRS群の個数(例えば、“正常”とグルー
プ化されたQRS群の個数)が全QRS群の個数の半数
より大きければ、そのクラスタは“正常”とされる。ま
た、クラスタ−の全QRS群の個数が値MIN−NEW
−NORMAL−POPより大きく、かつ、そのクラス
タの正常QRS群の個数が異常QRS群の個数の4倍よ
り大きければ、そのクラスタ−は“正常”とされる。あ
るQRS群がPVCとされた場合において、そのクラス
タに含まれる全QRS群の個数が予め定められた値MI
N−NEW−ABNORMAL−POPより大きく、か
つ、そのクラスタに含まれるPVCQRS群の個数が全
QRS群の個数の半数より大きければ、そのクラスタは
“PVC”とされる。また、クラスタ−の全QRS群の
個数が値MIN−NEW−ABNORMAL−POPよ
り大きく、かつ、そのPVCQRS群の個数が正常QR
S群の個数の2倍より大きければ、そのクラスタ−は
“PVC”とされる。最後に、あるQRS群が削除とさ
れた場合において、クラスタに含まれる全QRS群の個
数が予め定められた値MIN−NEW−DELETE−
POPより大きく、かつ、そのクラスタに含まれる削除
QRS群の個数が全QRS群の個数の半数より大きけれ
ば、そのクラスタは“削除”とされる。その後、プログ
ラムの作動はS20以下に戻り、新たなECG信号を処
理する。
【0040】次に、本発明の第2実施例について説明す
る。本実施例においても図2に示される構成の装置1が
使用され、その作動全般は図3のフロ−チャ−トに従っ
て行われるが、本実施例では図3のS30が図8のフロ
ーチャートに従って実行される。
【0041】図8のS311において、基準点評価関数
Q(i) が、式Q(i) = Q(fa (i), fb (i), f c (i), f
d (i)) によって求められる。但し、f a (i), f b
(i), fc (i) および f d (i)はそれぞれ、以下に詳述す
るように基礎ベクトルを発生させるフィルタにかけられ
るECG信号の出力である。本実施例の基準点評価関数
Q(i) は次式によって与えられる。 Q(i) = A|fa (i) |+B|fb (i) |+C|fc (i) |+D|fd (i) | 但し、A、B、C、Dは予め定められた定数である。
【0042】図9には、S311の作動工程における情
報の流れが図示されている。QRS検出/特徴情報抽出
装置2は、マイクロプロセッサ22のROM28内に記
憶されるソフトウエアによって構成され得るものであっ
て、第1、第2、第3、第4のフィルタ101、10
3、105、107を含む。各フィルタ−は受け取った
デジタル信号に基づいて特徴情報を出力する。装置2は
また、第1、第2、第3、第4の全波整流器109、1
11、113、115を含み、各全波整流器は対応する
フィルタの出力信号を受け取り、特徴情報信号を抽出す
る。装置2はさらに、加算器117を含み、整流された
各特徴情報信号を結合することにより単一の出力、すな
わち結合特徴情報信号ないし関数Q(i) とする。装置2
はまた、各フィルタ−101、103、105、107
に対応して出力端子119、121、123、125を
含み、前記S60と同様の処理のための特徴情報信号を
出力する。各フィルタ101、103、105、107
はそれぞれ、図10の示される基礎ベクトルのようなな
異なった基礎ベクトルを発生させる。フィルタ101、
103、105、107によって他の基礎ベクトルを発
生させることも可能であり、また、フィルタ101、1
03、105、107に幾つかのフィルタ(図示せず)
を加えることによってより大きな数の基礎ベクトルを発
生させることもできる。本実施例では、装置2はマイク
ロプロセッサ22に備えられたソフトウエアによって構
成されているが、専用の回路要素を使用して構成するこ
とも可能である。
【0043】前記基礎ベクトルは、デジタル信号に含ま
れるQRS群の形状に対応したフィルタ出力を発生させ
るように選択される。図9において、各フィルタは同一
のデジタル信号を受信するが、それぞれ他の全てのフィ
ルタの出力と異なった特性出力を発生させる。基礎ベク
トルは部分的に線型であって少数の線分を含んでいてよ
い。これらの基礎ベクトルにより、ACラインノイズや
基線ドリフト等のノイズが除去される。基礎ベクトルは
また、250Hzの頻度で発生させられ、そのウインド
は約100〜200msecの範囲である。こうして、
フィルタ101、103、105、107は、入力デシ
タル信号のあらゆる点について出力を発生させ、その特
徴情報に基づいて時間Tqrs 等の基準点の選択が可能と
される。各整流器109、111、113、115は、
対応するフィルタからの出力の整流された特徴情報信
号、すなわち、その絶対値を発生させる。全ての整流器
の出力は和算器117に集められ、関数Q(i) 等の特徴
情報信号の加算出力を発生させる。
【0044】S312において、基準点評価関数Q(i)
の極大値を特定することにより、その関数Q(i) から基
準点が決定される。本実施例の時系列Q(i) は、S30
2で与えられる関数M(t) と同様に機能し、関数Q(i)
の極大値が特定、利用される。特定された各極大値につ
いて、基準点Tqrs が特定される。こうして、基準点T
qrs が、基礎ベクトルによって得られる特徴情報信号の
極大値によって得られる。プログラムは次いで図3のS
40に進む。第2実施例における装置1の作動は図3の
残りの工程では前記と同一であるので説明を省略する。
【0045】図9の装置2は他のECG評価装置ないし
不整脈検出装置と伴に用いることが可能である。また、
装置2の出力は複数の特徴情報信号であり、前記のよう
に、N次元特徴情報空間にプロットされクラスタの形成
に使用され得る。
【0046】その他、当業者の知識、技能に基づいて種
々の変形、改良を施した態様で、本発明を実施すること
ができる。ゆえに、ここでは限られた実施例のみを説明
するが、本発明の趣旨を逸脱しない範囲内で数々の実施
例が考えられる。
【図面の簡単な説明】
【図1】典型的なECG波形の一部を示す図である。
【図2】本発明の不整脈検出装置を示す略図である。
【図3】QRS群をグループ化するために図2の装置に
よって実行される作動ステップの高レベルフローチャー
トを示す図である。
【図4】本発明の第1実施例の図3のステップS30に
おいて実行される作動ステップの低レベルフローチャー
トを示す図である。
【図5】本発明の第1実施例の図3のステップS50に
おいて実行される作動ステップの低レベルフローチャー
トを示す図である。
【図6】図3のステップS90におけるRS群評価のタ
イミングパターンを示す図である。
【図7】平均的QRS群幅と正常な範囲を示す図であ
る。
【図8】本発明の他の実施例において、図3のステップ
S30において実行される作動ステップの低レベルフロ
ーチャートを示す図である。
【図9】更に他の実施例において、図3の作動ステップ
S30における情報の流れを示す図である。
【図10】図9の基礎フィルタの出力信号を概念的に示
す略図である。

Claims (3)

    【特許請求の範囲】
  1. 【請求項1】 不整脈検出装置であって、 少なくとも1個のアナログECG信号を所定のサンプリ
    ング周期で取得する取得手段と、 前記取得手段に接続され、前記少なくとも1個のアナロ
    グECG信号を少なくとも1個のデジタル信号に変換す
    る変換手段と、 前記変換手段に接続され、前記少なくとも1つのデジタ
    ル信号に基づいて診断情報を発生させる処理手段と前記
    処理手段に接続され、前記診断情報を表示する表示手段
    とを含み、前記処理手段が、 複数の第1QRS群に基づいて正常なQRS群を判別す
    る第1手段と、 前記正常なQRS群に基づいて複数の第2QRS群の各
    々をグループ分けする第2手段とを含み、前記第1手段
    が、 前記少なくとも1個のデジタル信号に基づいて複数のス
    カラ信号を発生させる手段と、 前記複数のスカラ信号に含まれる複数のピーク値を検出
    することにより前記複数の第1QRS群の各々の発生時
    間を決定する手段と、 前記各第1QRS群についてQRS幅と複数の前半面積
    と複数の後半面積とを算出する手段とを含み、前記複数
    の前半面積が符号を有する前半面積と符号を有しない前
    半面積とを含み、前記複数の後半面積が符号を有する後
    半面積と符号を有しない後半面積とを含むことを特徴と
    する装置。
  2. 【請求項2】 ECG装置から供給される少なくとも1
    個のアナログECG信号から得られる少なくとも1個の
    デジタル信号に基づいて不整脈を検出する不整脈検出装
    置であって、該少なくとも1個のデジタル信号を取得す
    る入力回路を有するプロセッサを含む装置であって、 閾値交差基準に従いQRS群候補の発生時間に基づいて
    前記少なくとも1個のデジタル信号についてノイズ性若
    しくは非ノイズ性のノイズ属性を決定する決定手段と、 前記QRS群候補の発生時間と前記ノイズ属性とに基づ
    いて前記少なくとも1個のデジタル信号からN種類の特
    徴情報を抽出する抽出手段と、 前記QRS群候補を前記抽出されたN種類の特徴情報に
    従いN次元特徴情報空間にプロットするプロット手段
    と、 前記決定手段と前記抽出手段と前記プロット手段とを繰
    り返し作動させて複数の第1QRS群候補をプロットす
    ることによって、該複数の第1QRS群候補から少なく
    とも1個のクラスタを発生させる手段と、 前記少なくとも1個のクラスタに基づいて正常なクラス
    タを判別する手段と、 前記決定手段と前記抽出手段と前記プロット手段とを繰
    り返し作動させて複数の第2QRS群候補をプロットす
    ることによって、前記少なくとも1個のクラスタに該第
    2QRS群候補を加えるかないしは新たに少なくとも1
    個のクラスタを発生させるかの少なくとも一方を行う手
    段と、 前記デジタル信号に含まれる前記複数の第2QRS群候
    補の各々を、前記正常クラスタと該各第2QRS群候補
    から抽出された前記N種類の特徴情報とに基づいてグル
    ープ分けする手段とを含むことを特徴とする装置。
  3. 【請求項3】 ECG装置から供給される少なくとも1
    個のアナログECG信号から得られる少なくとも1個の
    デジタル信号に基づいて不整脈を検出する不整脈検出装
    置であって、該少なくとも1個のデジタル信号を取得す
    る入力回路を有するプロセッサを含む装置であって、 QRS群候補の発生時間に基づいて前記少なくとも1個
    のデジタル信号からN種類の特徴情報を抽出する抽出手
    段と、 前記QRS群候補を前記抽出されたN種類の特徴情報に
    従いN次元特徴情報空間にプロットするプロット手段
    と、 前記抽出手段と前記プロット手段とを繰り返し作動させ
    て複数の第1QRS群候補をプロットすることによっ
    て、該複数の第1QRS群候補から少なくとも1個のク
    ラスタを発生させる手段と、 前記少なくとも1個のクラスタに基づいて正常なクラス
    タを判別する手段と、 前記抽出手段と前記プロット手段とを繰り返し作動させ
    て複数の第2QRS群候補をプロットすることによっ
    て、前記少なくとも1個のクラスタに該第2QRS群候
    補を加えるかないしは新たに少なくとも1個のクラスタ
    を発生させるかの少なくとも一方を行う手段と、 真のQRS群と非QRS群信号の両方を含む前記複数の
    第2QRS群候補の最大限すべてを、前記正常クラスタ
    と該第2QRS群候補の各々から抽出された前記N種類
    の特徴情報とに基づいてグループ分けし、その際に該非
    QRS群信号を除去する手段と、 前記グループ分けされたQRS群に基づいて診断情報を
    表示する手段とを含むことを特徴とする装置。
JP27003891A 1990-09-21 1991-09-21 Ecg信号に基づく心臓不整脈検出装置 Expired - Fee Related JP3224569B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US58649290A 1990-09-21 1990-09-21
US586,492 1990-09-21

Publications (2)

Publication Number Publication Date
JPH0622915A true JPH0622915A (ja) 1994-02-01
JP3224569B2 JP3224569B2 (ja) 2001-10-29

Family

ID=24345968

Family Applications (1)

Application Number Title Priority Date Filing Date
JP27003891A Expired - Fee Related JP3224569B2 (ja) 1990-09-21 1991-09-21 Ecg信号に基づく心臓不整脈検出装置

Country Status (1)

Country Link
JP (1) JP3224569B2 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3816077A1 (de) * 1988-05-11 1989-11-23 Metallgesellschaft Ag Verfahren und vorrichtung zur abreinigung eines staubfilters
LT4756B (lt) 2000-01-11 2001-01-25 Kauno Medicinos Universiteto Psichofiziologijos Ir Reabilitacijos Institutas Sinusinio ir ektopinių širdies ritmų išskyrimo būdas
JP2013520275A (ja) * 2010-02-26 2013-06-06 ケイエヌユー−インダストリー コーポレイション ファウンデーション 心臓不整脈の診断方法及びその診断装置
US9008763B2 (en) 2009-04-06 2015-04-14 Cu Medical Systems, Inc. Arrhythmia diagnosis method and device
CN117281529A (zh) * 2023-11-23 2023-12-26 西南医科大学附属医院 一种心电图波形数据智能处理方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3816077A1 (de) * 1988-05-11 1989-11-23 Metallgesellschaft Ag Verfahren und vorrichtung zur abreinigung eines staubfilters
LT4756B (lt) 2000-01-11 2001-01-25 Kauno Medicinos Universiteto Psichofiziologijos Ir Reabilitacijos Institutas Sinusinio ir ektopinių širdies ritmų išskyrimo būdas
US9008763B2 (en) 2009-04-06 2015-04-14 Cu Medical Systems, Inc. Arrhythmia diagnosis method and device
JP2013520275A (ja) * 2010-02-26 2013-06-06 ケイエヌユー−インダストリー コーポレイション ファウンデーション 心臓不整脈の診断方法及びその診断装置
CN117281529A (zh) * 2023-11-23 2023-12-26 西南医科大学附属医院 一种心电图波形数据智能处理方法
CN117281529B (zh) * 2023-11-23 2024-01-30 西南医科大学附属医院 一种心电图波形数据智能处理方法

Also Published As

Publication number Publication date
JP3224569B2 (ja) 2001-10-29

Similar Documents

Publication Publication Date Title
US5271411A (en) Method and apparatus for ECG signal analysis and cardiac arrhythmia detection
Hossain et al. An accurate QRS complex and P wave detection in ECG signals using complete ensemble empirical mode decomposition with adaptive noise approach
Hamilton et al. Quantitative investigation of QRS detection rules using the MIT/BIH arrhythmia database
US5090418A (en) Method and apparatus for screening electrocardiographic (ECG) data
US9198634B2 (en) Medical decision support system
US7239988B2 (en) Apparatus and method for efficient representation of periodic and nearly periodic signals for analysis
Maglaveras et al. An adaptive backpropagation neural network for real-time ischemia episodes detection: development and performance analysis using the European ST-T database
US4336810A (en) Method and apparatus for arrhythmia analysis of ECG recordings
Tan et al. Detection of the QRS complex, P wave and T wave in electrocardiogram
EP2224346B1 (en) Ecg signal analysis tool
CN112528783B (zh) 心电图心搏数据聚类方法、装置、电子设备和介质
Smrdel et al. Automated detection of transient ST-segment episodes in 24h electrocardiograms
US20230309899A1 (en) Determination method and determination apparatus for beginning of t-wave, storage medium and computer program product
WO2021123271A1 (en) Device and process for ecg measurements
CN111557659A (zh) 基于多特征融合与Stacking-DWKNN的心律失常分类方法
CN111345814B (zh) 一种心电信号中心拍的分析方法、装置、设备和存储介质
Halder et al. Detection and identification of ECG waves by histogram approach
CN115062655A (zh) 一种多导联心电信号分析方法及多导联心肌梗死分析系统
JP3224569B2 (ja) Ecg信号に基づく心臓不整脈検出装置
Shahram et al. ECG beat classification based on a cross-distance analysis
Young et al. A comparative study of a hidden Markov model detector for atrial fibrillation
JPH05207985A (ja) 心電図波形認識方式
KR980008168A (ko) 심전도 신호의 큐.알.에스 컴플렉스 특성 표시 장치 및 방법
Atanasoski et al. Unsupervised classification of premature ventricular contractions based on RR interval and heartbeat morphology
Sörnmo A model-based approach to QRS delineation

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees