JPH1170089A - 心臓波形特性表示方法 - Google Patents

心臓波形特性表示方法

Info

Publication number
JPH1170089A
JPH1170089A JP10198968A JP19896898A JPH1170089A JP H1170089 A JPH1170089 A JP H1170089A JP 10198968 A JP10198968 A JP 10198968A JP 19896898 A JP19896898 A JP 19896898A JP H1170089 A JPH1170089 A JP H1170089A
Authority
JP
Japan
Prior art keywords
signal
lead
waveform
sqi
sampling
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.)
Pending
Application number
JP10198968A
Other languages
English (en)
Inventor
Cheuw Lee Yeo
チュー・リー・イエオ
Jyh-Yun J Wang
ジー−ユン・ジェイ・ワン
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.)
HP Inc
Original Assignee
Hewlett Packard Co
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 Hewlett Packard Co filed Critical Hewlett Packard Co
Publication of JPH1170089A publication Critical patent/JPH1170089A/ja
Pending legal-status Critical Current

Links

Classifications

    • 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/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle
    • A61B5/35—Detecting specific parameters of the electrocardiograph cycle by template matching
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7221—Determining signal validity, reliability or quality
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S128/00—Surgery
    • Y10S128/901—Suppression of noise in electric signal

Landscapes

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

Abstract

(57)【要約】 【課題】 波形表現を心電計測定値から得る場合に、Ｑ
ＲＳ群のような心臓機能の波形表現を高精度に特性表示
する。 【解決手段】 心臓の機能を示す一つ以上の信号を測定
し、カージナル数の点で測定した一つ以上の信号のそれ
ぞれをサンプリングする。サンプリングした一つ以上の
信号のそれぞれをカージナル数のサンプリング点での最
大サンプリング値を決定する。更に、サンプリングした
一つ以上の信号のそれぞれを決定された最大値をサンプ
リングした一つ以上の信号のそれぞれの関連する除数で
除算することによって信号の良否の指数を計算する。

Description

【発明の詳細な説明】
【０００１】
【発明の属する技術分野】本発明は、ＱＲＳ群のような
心臓機能の波形表現を、心電計測定値から得る際に表示
する心臓波形特性表示技術に関する。
【０００２】本発明は、波形表現が心電計測定値の場合
に、ＱＲＳ群のような心臓の機能の波形表現を、自動的
に特性表示として行うものである。この波形表現の特性
表示のために、（１）分析される構成要素が心臓機能に
対する高精度の表現であるか又は低い精度の表現である
かを決定するために波形表現の一組の心電計測定値の要
素を分析する。（２）この後、波形表現の特性表示の高
精度の信号を使用する。
【０００３】本発明の一つの例示的な実施形態では、心
臓機能を監視する心電計から得られた一定の特定信号を
利用している。これら一定の特定信号を利用は、（１）
心臓の電気化学的及び機械的動作、（２）心臓の電気化
学的動作が、心電計で心臓の機械的動作を図式的に示す
際に用いられる電気エネルギに変換する処理、及び
（３）一定の特定信号（又は「リード」）を心電計から
得る処理である。
【０００４】
【従来の技術】従来、心臓の機械的事象は、先に心臓の
電気化学的活動（すなわち、活動電位の伝播）があって
開始される。心臓の電気化学的活動を人が目視できる形
に変換する心電計は、心臓の電気化学的活動の可視表現
が可能である。この可視表現は心電図（ＥＫＧ）であ
る。
【０００５】ＥＫＧの期間中、電極が身体の表面に取付
けられる。電極は、電極内部の電荷担体が電気化学的交
換により身体内部の電荷担体（イオン）に対応できるよ
うになっている。電極を身体表面に取付けると信号を適
切に増幅して身体内部の電圧変化を記録できるようにな
る。ＥＫＧ機械内部のガルヴァノメータが記録装置とし
て使用される。ガルヴァノメータは、二つの電極の間の
電位差を記録する。ＥＫＧは、身体表面上の二つの電極
の間の電圧差を単に時間の関数として記録するものであ
り、通常帯形チャートとして記録される。心臓が休息し
ている、すなわち弛緩期にあるとき、心臓の細胞は分極
しており、電荷移動が生じない。したがって、ＥＫＧの
ガルヴァノメータは、振れを記録しない。しかし、心臓
が活動電位を伝播し始めると、減極が生じている電極
が、心臓がまだ減極していない身体上の領域から電位差
を記録する。したがって、ガルヴァノメータが振れる。
【０００６】完全な心臓サイクルが心搏である。ＥＫＧ
上では、正常な心搏が、明確な信号を有している。最
初、ガルヴァノメータは、比較的短い持続時間の丸みの
ある正の振れのＰ波を記録するが、これは心房の減極に
よって発生する。これに続いて、小さいが鋭い負の振れ
のＱ波が存在する。また、極めて大きく、鋭い正の振れ
のＲ波が存在し、その後に鋭く、大きい負の振れのＳ波
が存在する。これらの波を一緒にしてＱＲＳ群という。
ＱＲＳ群は、心室の減極により生ずる。ＱＲＳ群に続い
て、比較的長い持続時間の丸みのある正の振れのＴ波が
あり、これは心室の再分極によって発生する。
【０００７】ＥＫＧは、多数組の電極を使用する。しか
し、これら電極は、受信した信号が前記説明と同一の波
形形状を有するように身体の表面に配置される。周知の
バイポーラ対（二つて一組）の電極が通常、患者の右腕
（ＲＡ）、左腕（ＬＡ）、右脚（ＲＬ）（これらが通常
基準として使用される）、及び、左脚（ＬＬ）に配置さ
れる。参照用の単極電極がＶリードであり、解剖学的に
確定した約束事に従って患者の胸に配置される（リード
Ｖ１−Ｖ６として以降のように名付けられる）。心臓の
監視及び診断に際しては、二つの、このような電極の間
又は一つの電極と一群の他の電極の平均との間に現われ
る電圧差が、心臓の電気的活動の特定の眺望を表し、こ
れが一般的なＥＫＧである。電極の特定の組合せがリー
ドである。例えば、標準の１２リード心電図システムに
使用することが出来るリードは、下記のとおりである。 リードＩ＝（ＬＡ−ＲＡ） リードII＝（ＬＬ−ＲＡ） リードIII＝（ＬＬ−ＬＡ） リードＶ１＝Ｖ１−（ＬＡ＋ＲＡ＋ＬＬ）／３ リードＶ２＝Ｖ２−（ＬＡ＋ＲＡ＋ＬＬ）／３ リードＶ３＝Ｖ３−（ＬＡ＋ＲＡ＋ＬＬ）／３ リードＶ４＝Ｖ４−（ＬＡ＋ＲＡ＋ＬＬ）／３ リードＶ５＝Ｖ５−（ＬＡ＋ＲＡ＋ＬＬ）／３ リードＶ６＝Ｖ６−（ＬＡ＋ＲＡ＋ＬＬ）／３ リードａＶＦ＝ＬＬ−（ＬＡ＋ＲＡ）／２ リードａＶＲ＝ＲＡ−（ＬＡ＋ＬＬ）／２ リードａＶＬ＝ＬＡ−（ＲＡ＋ＬＬ）／２ したがって、用語「リード」は物理的ワイヤを示す場合
に用いられるが、心電計では上記に示したように一定の
電極構成から取得される電気信号を意味している。
【０００８】多年にわたり、保健の専門家は、異型のＥ
ＫＧを種々の病気や心臓欠陥と関連付けている。この関
連付けプロセスが「心電図記録法」である。
【０００９】この心電図記録法は、主に心電計測定値の
波形記録の変化に基づいて、視覚的に検査する技術であ
る。このような診断では、波形の不正確が診断の不正確
となる点で心電計測定値から得た波形表現の正確な特性
表示、例えば、ＱＲＳ群が現われる場所に関する知識を
有していることが不可欠である。
【００１０】波形表現の正確な特性表示（例えば、ＱＲ
Ｓ群がどこに存在するか）の必要性は、機械化（自動
化）心電図記録法にとっては更に決定的である。すなわ
ち、通常行なわれた多数の機能を自動化して取得する装
置がある。この自動化装置は心臓専門家の診断を、一つ
以上の波形の雛型を患者の心臓機能の波形表現と比較す
る基準として使用することによって、その近似化を行う
ことである。
【００１１】このような雛型を通常下記のようにして利
用する。最初に、多数リード心電計を特定の患者に貼り
つけ、心電図の測定を開始する。第２に、ＥＫＧ波形を
一つ以上のリードから得る。第３に、関係のある一定の
心臓状態を示す一つ以上の波形雛型を、患者の保管波形
と相互に関連づける。受け取った波形を波形雛型に対し
て走査し、波形雛型と患者のＥＫＧ波形との類似を考慮
して、波形雛型を患者のＥＫＧ波形と相互に関連づけ
る。基準雛型が患者の受信データとよく相関しているこ
とが確認されれば、患者が雛型に関連する心臓状態にあ
ると判断される。
【００１２】このような相関処理で正確な診断を得るた
めには、患者のＥＫＧ波形（心臓機能の波形表現）の特
性表示（例えば、ＱＲＳ群がどこに存在するか）が可能
な限り正確であることが重要である。これは、このよう
な相関処理が患者のＥＫＧ波形内部の一定のＥＫＧ特徴
（例えば、ＥＫＧデータ内のＱＲＳ群）の形状間の形態
学的差異又は類似及び患者のＥＫＧ波形に相関している
雛型に関して診断しているために重要である。波形の特
性表示の不正確が時として不正確な診断を生ずるという
点では、機械化した装置によって心臓機能の心電計測定
値から得られた波形表現の正確な特性表示（例えば、Ｑ
ＲＳ群がどこに存在するか）を保有していることが重要
となる。
【００１３】前記の説明のように、患者のＥＫＧ波形デ
ータの正確な特性表示は、多くの理由から極めて重要で
ある。しかし、このような正確な特性表示は、様々なＥ
ＫＧリードが時として「汚れた（雑音重畳）」信号又は
歪んだ信号を発生するので、その達成が困難である。こ
れらの歪みは、ＥＫＧ電極の準備が不十分で患者との電
気的接続が不正確になる場合、及び、患者の移動又は筋
肉の痙攣から生ずる低周波雑音又は高周波雑音のよう
な、多数の要因から発生する。このような歪信号が存在
すると、一つ以上の「汚れた」又は歪信号は、このよう
な特性表示を困難にするため、ＥＫＧ波形表現を適切に
特性表示することが難しい。すでに説明したように、心
臓専門家又は機械化した装置がＥＫＧデータを利用して
心臓機能の正確な診断を行なうためには、ＥＫＧ波形の
正確な特性表示が必要である。
【００１４】多数リードに基づいて患者のＥＫＧ波形デ
ータの正確な特性表示を行う試みでは、多数リードが全
てみ合わされてＱＲＳ検出のための単独信号を形成す
る。これまでの多数のＥＣＧ信号の組み合わせに使用す
る方法には、（１）信号を加算する処理、又は、（２）
信号を掛け合わせる処理がある。しかし、多数リードを
加算することは、処理中のリードの大多数に雑音が存在
する場合、信号対雑音比（Ｓ／Ｎ比）の改善に役立つの
みである。更に、単独リード内の十分な大きさの雑音
は、加算すると最終検出信号に悪影響を及ぼす場合があ
る。他方、掛算処理で多数リードを組み合わせること
は、加算処理の問題を解決するが、単独リードに等電位
のＱＲＳ群がある際に発生する偽陰性のような他の問題
が発生する場合がある。
【００１５】
【発明が解決しようとする課題】このように上記従来例
では種々の問題がある。すなわち、これまで説明したよ
うに多種のニーズがあり、波形表現を心電計測定値から
得る場合に、ＱＲＳ群のような心臓機能の波形表現の高
精度の特性表示を得る必要がある。
【００１６】本発明は、このような従来の技術における
課題を解決するものであり、ＱＲＳ群のような、心臓機
能の波形表現を高精度に特性表示する改善された心臓波
形特性表示技術の提供を目的としている。
【００１７】更に、本発明は波形表現を心電計測定値か
ら得る場合に、ＱＲＳ群のような、心臓機能の波形表現
を高精度に特性表示する心臓波形特性表示方法及びその
システムの提供を他の目的としている。
【００１８】
【課題を解決するための手段】上記目的を達成するため
に、本発明は、波形表現を心電計測定から得る場合に、
ＱＲＳ群のような、心臓機能の波形表現を高精度に特性
表示するために、心臓の機能を示す少なくとも一つの信
号（一つ以上）を測定してサンプリングする。この少な
くとも一つの信号の良否の指数を測定してサンプリング
した信号に対して計算を行う。信号の良否の指数を、一
つ以上の測定し、かつ、サンプリングした信号に対する
最大サンプリング値を、計算した除数で除算することに
よって算出する。一つ以上の測定し、かつ、サンプリン
グした信号のそれぞれの計算を行った信号の良否の指数
に応じて、一つ以上の重み付け係数を計算する。そし
て、心臓機能を示す一つ以上の測定し、かつ、サンプリ
ングした信号のそれぞれに計算した重み付け係数を掛け
合わせた波形表現を生成する。
【００１９】本発明の、上記目的及び他の目的、特徴及
び長所は下記の詳細な説明で明らかになるものである。
【００２０】本発明の新規な特徴を特許請求の範囲に示
してある。しかし、本発明にあって使用の好適を態様、
更に他の目的及びその長所は、添付した図面と共に例示
的な実施形態の詳細な説明によって、より良く理解され
るものである。
【００２１】
【発明の実施の形態】次に、本発明の心臓波形特性表示
方法及びそのシステムの実施の形態を添付図面を参照し
て詳細に説明する。この実施形態は、波形表現が心電計
測定値から構成されている場合の、ＱＲＳ群のような心
臓機能の波形表現を高精度の特性表示として得るもので
ある。
【００２２】心電計リード（すなわち、上記説明の信
号）はアナログ信号であるが、それぞれの心電計リード
を時間サンプリングする。したがって、以下の説明の大
部分は、このような個別にサンプリング（標本化）かつ
量子化したデジタル信号を取り扱うが、特定の場合に
は、アナログ波形を参照すると概念が明瞭に理解される
ため、リードをアナログ的に説明する。
【００２３】この実施形態では、心電計測定値に基づき
心臓機能の波形表現の特性表示を生成する。波形表現は
「活動関数」を利用して生成される。「活動関数」は、
一つ以上の個別心電計リード（信号）の複合体である。
以降で説明するように「活動関数」は、ＱＲＳ群のよう
な、心臓機能の波形表現を極めて正確に特性表示する。
すなわち、「活動関数」を構成するのに使用される一つ
以上の別個の心電計リードは、（１）「良い（雑音のな
い）」信号を示す規定された基線量を超えるよう決定さ
れ、かつ、（２）信号が「良い」ほど大きい重みを与え
て、固有的の「活動関数」が得られるようにしている。
活動関数の説明前に、信号の良否を判断する処理及び重
み付け係数を、信号良否との関連について説明する。
【００２４】本発明の実施形態では、その一つの処理と
して、心電計リードの信号の良否を判断している。すな
わち、どのリードが「良い」信号であるかを決定してい
る。心電図記録法は、人心臓専門家が実施し、主として
視覚的な判断技術である。すなわち、例えば、１２個か
らなる一組の心電計リード（信号）が得られるとリード
を視覚検査に基づいて最良の質から最悪の質まで等級を
付けを行う。このような視覚検査は、主として心臓専門
家にとって主観的である。したがって、この処理技術に
おける大きな課題は、心臓専門家が検査を主観的に行な
っているものを機械化することである。
【００２５】本発明の実施形態では、「信号の良否の指
数」という量を導入し、視覚検査によって行っていた検
査処理を機械化して達成する。
【００２６】以上説明したように、本発明では、信号の
「良さ」は、その「信号の良否の指数」に対応してい
る。一般に（特殊の場合を除く。この場合については以
降で説明する）特定のｉ番目のリードに対するＳＱＩ
は、以下の数式１のように規定される。
【００２７】
【数１】
【００２８】ここで、Ｖ(i)(j)は、ｉ番目のリードに対
するｊ番目のデータサンプリングでのＥＫＧ値である。
ｉは、考察中のリードであり、１からＮまでの範囲であ
る。ただし、Ｎは存在するリードの総数である。ｊは、
考察中のリードｉに対するｊ番目のＥＫＧデータサンプ
リングである。ｎは、ＳＱＩを計算するサンプリング窓
の前縁サンプリング時間である（ここでのサンプリング
窓の幅は５０サンプリングである）。
【００２９】ＳＱＩの定義は、ＳＱＩが存在するそれぞ
れのリードに対して計算される。この数式は、５０個の
サンプリングを行う場合、それぞれのリード内で発生す
る最大絶対値が、ある時間基線にわたり決定される。こ
の決定された最大絶対値にサンプリング数の５０が乗ぜ
られる。この乗算結果は、５０サンプリング点の窓で取
得されたリードのＥＫＧサンプリング値の和で除算され
る。
【００３０】上記の数式は、あらゆる時間軸に適用でき
る。例えば、一つの例示実施形態において、ＳＱＩは、
４００ミリ秒の窓（すなわち、５０個の８ミリ秒データ
サンプリング、又は、８ミリ秒の間隔で取得された５０
個のデータサンプリング）で２００ミリ秒ごとに計算さ
れる。ＳＱＩは最初、５０個の等間隔データサンプリン
グからなる４００ミリ秒間で計算される。続いて、ＳＱ
Ｉは２００ミリ秒ごとに計算されるが、これはＳＱＩが
２５個の新しいサンプリング及び２５個の前のサンプリ
ングに対して計算される。すなわち、４００ミリ秒幅の
前縁は再計算前に２００ミリ秒前進するのみであるため
である。
【００３１】更に、ＳＱＩの上記定義は、５０サンプリ
ングを取り扱うだけであるが、上記のＳＱＩがいずれの
数のサンプリングを有する基線でも計算できる。
【００３２】以上の説明のように「活動関数」は、「良
い」信号のみであり、「活動関数」は、それらの活動関
数への寄与が「不正確な」質の信号よりも大きな「良
い」信号に重み付けして生成される。したがって、重み
付け係数は、多種のリードを判定した「良さ」又は信号
の良否によって決定される。この関係は、重み付け係数
の以下の数式２による定義で明確である。この定義は、
重み付け係数が信号の良否、又は、ＳＱＩによって決定
されることを明瞭に示している。
【００３３】
【数２】
【００３４】ここで、ｉは、考察中のリードである。Ｓ
ＱＩ(i)は、ｉ番目のリードの信号の良否の指数であ
る。Ｎは、存在するＥＫＧリードの総数である。したが
って、上記の重み付け係数の数式は、下記に規定する例
外（その効用は以降で説明する）を除いて一般的に当て
はまることになる。
【００３５】次にＳＱＩ及び重み付け係数を規定して、
「活動関数」を以下の数式３で規定するが、これはサン
プリングに基づいて算出される。
【００３６】
【数３】
【００３７】ただし、ｋは、ｋ番目のＥＫＧデータサン
プリングである。Ｗ(i)は、ｉ番目のリードに対する重
み付け係数である。Ｎは、ＥＫＧリードの総数である。
Ｖ(i)(k)は、ｉ番目のリードに対するｋ番目のデータサ
ンプリングでのＥＫＧ値である。
【００３８】ＱＲＳ群のような、心臓機能の波形表現を
修正された（すなわち、極性が重要でない）場合にも利
用できる。これは「活動関数」がサンプリングの絶対値
を利用しているためである。
【００３９】図１は、活動関数に関係する概念を説明す
るための部分概略図である。図１に示されているのは４
個の例示的な心電計信号（リード）である。心電計信号
１００は、理想的なＱＲＳ群である。心電計信号１０
２，１０４，１０６もＱＲＳ群であるが、理想的ではな
いＱＲＳ群を表している。
【００４０】活動関数は、心電計信号の絶対値を利用し
ている。したがって、図１に示したのは４個の心電計信
号１００，１０２，１０４，１０６の、それぞれ修正信
号１０８，１１０，１１２，１１４である。
【００４１】上記で説明したように、活動関数は、心電
計信号に関連して重み付け係数を利用し、「最良の」心
電計波形が最大の寄与を行なう複合関数を発生する。し
たがって、図１に示した修正ＱＲＳ群１０８（「最良
の」波形である理想ＱＲＳ群と同等に扱う）が、重み付
け係数５／８が設定される（ここに示した重み付け係数
は単に関連する概念を示すのに提示しただけであり、図
示した信号から得られた実際の重み付け係数ではな
い）。更に、修正ＱＲＳ群１１０，１１２，１１４（そ
れぞれＱＲＳ群１０２，１０４，１０６と同等に扱う）
には、それぞれ重み付け係数１／８が設定される。
【００４２】図示した修正ＱＲＳ群には全て、上記で説
明した活動関数によって、それぞれの重み付け係数が掛
けられ、更に、加算されて複合ＱＲＳ群１２０を生成す
る。図１の下部は、第１モジュール２１０と第２モジュ
ール２２０との間に接触が確立して複合ＱＲＳ群（すな
わち、活動関数）から、最良のＱＲＳ群１００により、
大きい部分に、複合ＱＲＳ群が形成されている。すなわ
ち、図１は、心電計波形の高精度の表現を得ていること
を示している。このような表現は、最良の信号を有する
波形に最も大きい重み付け係数が設定されるためであ
る。
【００４３】以上説明したように、特定のリードを判定
した「良さ」（すなわち、信号の良否）は、「信号の良
否の指数」又は「ＳＱＩ」により決定される。ＳＱＩは
一般的に、一定の時間軸上に存在するそれぞれのリード
について、その時間軸で取得される一定数のサンプリン
グを利用して算出される。しかし、「ＳＱＩ」を計算す
る前に調べる特殊事例が存在する。これら特殊事例は、
ＳＱＩが計算される通常の処理を、一つ以上のリードが
特定の不満足な信号を明示するという事実のため、修正
すべきであるか否かを決定するのに使用される。以下に
示すように、信号が特に不満足であることが決定されれ
ば、このような信号には、計算せずに、所定のＳＱＩが
設定される。
【００４４】チェックされる最初の特殊事例は、それぞ
れのリードが弱い信号を示す場合である。このような場
合には、上記で示したＳＱＩの式を無視し、代わりにそ
れぞれのリードに対するＳＱＩに０を設定する。上記か
ら、存在するリードの信号の良否の指数が全て０（すな
わち、ｉ＝１からＮ（ただしＮは存在する全てのリード
に等しい）に対してＳＱＩ(i)=０）であることが判明す
れば、重み付け係数は、ｉ＝１からＮ（ただしＮは存在
する全てのリードに等しい）に対して、Ｗ(i)=１／Ｎで
あるように規定される。換言すれば、全てのリードが弱
いか、又は、不満足であることが判明すれば、それぞれ
のリードに対する信号指数に対して、それぞれの信号に
与えられる等しい重みを発生させるための０が設定され
る。
【００４５】すなわち、このような場合にはそれぞれの
リードに等しく重みづ係数が設定される。一つの例示実
施形態では、それぞれのリードが弱い信号として決定さ
れる判定基準が以下の数式４で表される。
【００４６】
【数４】
【００４７】ここで、ｉは、考察中のリードであり、Ｎ
を存在するリードの総数として１からＮまでの範囲にあ
る。ｊは、考察中のリードｉに対するｊ番目のＥＫＧデ
ータサンプリングである。ｎは、それにわたってＳＱＩ
を計算するサンプリング窓（示した等式では、サンプリ
ング窓の幅は５０サンプリングである）の前縁サンプリ
ング時間である。Ｖ(i)(j)は、ｉ番目のリードに対する
ｊ番目のデータサンプリングでのＥＫＧ値である。
【００４８】この等式は、それぞれの信号リードにおけ
る、ある５０サンプリングのブロックの中で取得した最
大サンプリング値が２５０マイクロボルト以下であれ
ば、それぞれの信号リードは低振幅である弱い信号を示
している。したがって、このような場合、全てのリード
に対するＳＱＩが０に設定されることになる。これはそ
れぞれのリードに対して全て１／Ｎに設定される重みづ
け係数を設定する。
【００４９】図２は、全ての信号が弱い特殊事例に関係
する概念を示す部分概略図である。図２は、４個の例示
心電計信号（リード）を示している。心電計信号２０
０，２０２，２０４，２０６は全てＱＲＳ群であるが、
ＱＲＳ群の全ての構成要素の大きさが２５０マイクロボ
ルト以下であることを示している。
【００５０】活動関数は、心電計信号の絶対値を利用す
る。したがって、図２に示されているのは、４個の心電
計信号２００，２０２，２０４，２０６のそれぞれ修正
形態２０８，２１０，２１２，２１４である。
【００５１】上記の説明のように、全ての心電計信号が
弱い図２に示す特殊事例では、それぞれのリードに対す
るＳＱＩに０が設定されている。更に、上記で説明した
ように、全てのリードのＳＱＩが０である場合には、重
み付け係数は、全てのリードに１／Ｎの重み付け係数が
設定される。ただし、Ｎは存在するリードの総数であ
る。したがって、図２に示したのは、修正ＱＲＳ群２０
８，２１０，２１２，２１４（それぞれＱＲＳ群２０
０，２０２，２０４，２０６と同等に扱う）にそれぞ
れ、４個のリードが存在する（すなわち、Ｎ＝４）から
１／４の重み付け係数が設定されることである。
【００５２】図示した修正ＱＲＳ群の全てに（上に説明
した活動関数により）そのそれぞれの重み付け係数が掛
けられ、加算されて複合ＱＲＳ群２２０を発生する。複
合ＱＲＳ群（すなわち、活動関数）から、図２に示した
弱い信号により、等しい部分に、複合ＱＲＳ群が形成さ
れていることが判明する。したがって、図２は、他の信
号より「良い」信号が実質的に存在しないことが決定さ
れた場合に、全ての信号に等しく重み付けすることを示
している。
【００５３】これまではチェックされる第１事例につい
て説明した。第２特殊事例は、少なくとも一つの良質の
リード、及び一つ以上の弱い（低振幅）リードが存在す
る場合である。このような場合は、弱いリードのＳＱＩ
は０に等しく設定されるが、弱いリードの判定基準に適
合しないリードのＳＱＩが、通常、ＳＱＩの上記の数式
を利用して計算される。
【００５４】一つの例示実施形態では、少なくとも一つ
の良質のリードが存在し、一つ以上の弱い（低振幅）リ
ードが存在するかを決定する明白な判定基準が以下の数
式５で示される。
【００５５】
【数５】
【００５６】ここで、ＳＱＩ(i)は、リードについて計
算した信号の良否の指数である。ｉは、考察中のリード
であり、Ｎを存在するリードの総数として、１からＮま
での範囲にある。ｊは、考察中のリードｉに対するｊ番
目のＥＫＧデータサンプリングである。ｎは、ＳＱＩを
計算するサンプリング窓（示した等式におけるサンプリ
ング窓の幅は５０サンプリングである）の前縁サンプリ
ング時間である。Ｖ(i)(j)は、ｉ番目のリードに対する
ｊ番目のデータサンプリングでのＥＫＧ値である。
【００５７】この数式は、少なくとも一つのリードが３
より大きいＳＱＩ計算値を有し、少なくとも一つのリー
ドが３００マイクロボルト以下の最大サンプリング値を
有していれば、信号の最大サンプリング値が３００マイ
クロボルト以下であるリードに対してＳＱＩを０になる
ように設定する。その他の場合には、残りのリードに対
するＳＱＩを上記のＳＱＩの数式を用いて算出する。
【００５８】図３は、特殊事例に関係する概念を説明す
るための部分概略図である。図３は、４個の心電計信号
（リード）示し、心電計信号３００，３０２，３０４，
３０６は全てＱＲＳ群である。図示の心電計信号３０
２，３０６は、心電計信号３００，３０４より大きさ
が、はるかに小さく、形状が悪いことを示している。
【００５９】活動関数は、心電計信号の絶対値を利用し
ている。したがって、図３はそれぞれ４個の心電計信号
３００，３０２，３０４，３０６の修正形３０８，３１
０，３１２，３１４を示している。
【００６０】上記の説明のように、少なくとも一つの心
電計信号が「良い」ことが判明した場合、一つ以上の心
電計信号が「弱い」ことが判明している図３の特殊事例
では、「弱い」リードのそれぞれに対するＳＱＩに０が
設定される。したがって、全てのリードに対するＳＱＩ
は、上記のＳＱＩの数式を使用して算出される。
【００６１】図３では、心電計信号３００，３０４は
「良い」状態を示している。したがって、例示のため
に、心電計信号３００，３０４に対するＳＱＩを上記の
ＳＱＩの数式を使用して算出する際に、それぞれに対す
るＳＱＩは共に、関連している第１閾値、例えば、上記
の数式３を超えていると仮定する。この後、全てのリー
ドに関する信号強度を第２関連閾値、例えば、上記の特
殊方程式における３００マイクロボルト）と比較する。
信号強度が第２関連閾値以下であるリードについて、信
号は「弱い」と判断され、このようなリードに対するＳ
ＱＩを０に設定する。
【００６２】その後で、上記に示した標準重み付け係数
方程式を使用して全てのリードに重み付け係数を設定す
る。この標準重み付け係数方程式を検討すると、リード
のＳＱＩが０であれば、そのリードに対する重み付け係
数も０であることが判明する。したがって、図３では、
修正ＱＲＳ群３１０及び３１４に０の重み付け係数が設
定される（このような「弱い」リードのＳＱＩは０であ
ると規定されている）ことを示している。更に、図３
は、二つのリードが残っており、この残りのリードは同
等に良い信号であるＳＱＩ指数を有しているので、修正
ＱＲＳ群３０８，３１２（それぞれＱＲＳ群３００，３
０４と同等に扱う）にそれぞれ１／２の重み付け係数が
設定されることを示している。なお、ここに示した重み
付け係数は単に関連する概念を示すものであって、図示
した信号から得られた実際の重み付け係数ではない。
【００６３】図示した修正ＱＲＳ群に全て（前記で説明
した活動関数により）そのそれぞれの重み付け係数を掛
けて加算し、複合ＱＲＳ群３２０を生成する。複合ＱＲ
Ｓ群（すなわち、活動関数）から、図３に示した「強
い」信号によって、等しい部分に、複合ＱＲＳ群が形成
されていること、及び「弱い」信号の寄与が無視されて
いることが判明する。したがって、図３は、過度に「弱
い」と考えられる信号の活動関数への関与を効果的に排
除できることを示している。
【００６４】以上、第２特殊事例を説明した。更に、第
３事例は、一つ以上のリードの中の低周波雑音及び高周
波雑音の「移動平均」が、このようなリードが元来信頼
できないことを示している場合に、リードのＳＱＩを０
に等しく設定して、「活動関数」に現われないように除
去する。
【００６５】高周波雑音に対する「移動平均」は、「高
周波雑音指数」であり、ＱＲＳ群の検出に続く時間サン
プリングで更新される以下の数式６で規定される。
【００６６】
【数６】
【００６７】ここでｉは、考察中のリードであり、Ｎを
存在するリードの総数として、１からＮまでの範囲にあ
る。ｊは、考察中のリードｉに対するｊ番目のＥＫＧデ
ータサンプリングである。上バーが無いＨＦＮＩ(i)量
は現在の（後に規定する上バーの無いＨＦＮＩ(i)の数
式により）高周波雑音指数値を示し、上バーを付したＨ
ＦＮＩ(i)量は高周波雑音指数の移動平均を示してい
る。
【００６８】すでに説明したように、高周波雑音の「移
動平均」がリードが元来信頼できないことを示していれ
ば、このようなリードは活動関数を形成する構成要素か
ら除外される。「移動平均」が信頼できないことを示す
か否かは、以下の数式７で決定される。
【００６９】
【数７】
【００７０】低周波雑音に対する「移動平均」は、「低
周波雑音指数」であり、ＱＲＳ群の検出に続く時間サン
プリングで更新される以下の数式８で規定される。
【００７１】
【数８】
【００７２】ここでｉは、考察中のリードであり、Ｎを
存在するリードの総数として、１からＮまでの範囲にあ
る。ｊは、考察中のリードｉに対するｊ番目のＥＫＧデ
ータサンプリングである。上バーが無いＬＦＮＩ(i)量
は現在の計算された（後に規定する上バーの無いＬＦＮ
Ｉ(i)の数式により）低周波雑音指数値を示し、上バー
のあるＬＦＮＩ(i)量は低周波雑音指数の移動平均を示
している。
【００７３】すでに説明したように、低周波雑音の「移
動平均」がリードが信頼できないことを示していれば、
このようなリードは活動関数を形成する構成要素から除
外される。「移動平均」が信頼できないことを示すか否
かは、以下の数式９で決定される。
【００７４】
【数９】
【００７５】移動平均を計算するための上記の数式は、
高周波雑音及び低周波雑音の測定値を使用している。高
周波雑音及び低周波雑音の測定値は下記のように、心電
計リードから得られる。
【００７６】高周波雑音指数は、信号に対する高周波雑
音の比として計算される。高周波雑音を求めるために、
入力ＥＫＧ信号Ｖ(j)を、サンプリングごとに以下の数
式１０で計算される第２差分フィルタを使用して高域を
濾波（フィルタリング）する。
【００７７】
【数１０】
【００７８】高周波雑音指数をＱＲＳ群の検出に続いて
以下の数式１１で算出する。
【００７９】
【数１１】
【００８０】ここで、ｉは、考察中のリードであり、Ｎ
を存在するリードの総数として、１からＮまでの範囲に
ある。ｊは、考察中のリードｉに対するｊ番目のＥＫＧ
データサンプリングである。Ｖ(i)(j)は、ｉ番目のリー
ドに対するｊ番目のタサンプリングでのＥＫＧ値であ
る。ＲＨＴk (i)はｉ番目のリードのｋ番目のサンプリン
グで検出されたＲ波の振幅の絶対値に等しい。加算はＲ
波の中心かつ８ミリ秒間隔で１２５回を行なう。Ｘ
(i)(j)は高域の濾波後の出力である。
【００８１】低周波雑音指数は、信号に対する低周波雑
音の比として算出される。低周波雑音を求めるために、
入力ＥＫＧ信号Ｖ(j) を、サンプリングごとに計算され
る以下の数式１２を用いて低域を濾波する。
【００８２】
【数１２】
【００８３】低周波雑音指数をＱＲＳ群の検出に続いて
以下の数式１３で算出する。
【００８４】
【数１３】
【００８５】ここで、ｉは、考察中のリードであり、Ｎ
を存在するリードの総数として、１からＮまでの範囲に
ある。ｊは、考察中のリードｉに対するｊ番目のＥＫＧ
データサンプリングである。Ｖ(i)(j)は、ｉ番目のリー
ドに対するｊ番目のタサンプリングでのＥＫＧ値であ
る。ＲＨＴk (i)はｉ番目のリードのｋ番目のサンプリン
グで検出されたＲ波の振幅の絶対値に等しい。加算はＲ
波の中心で８ミリ秒間隔で１２５回を行なう。Ｙ(i)(j)
は低域の濾波後の出力である。
【００８６】図４は以上説明した特殊事例に関係する概
念説明するための部分概略図である。ここでは図示した
信号の二つが受け入れ可能であるが、図示した信号の一
方は過度の高周波雑音を有し、図示した信号の他方は過
度の低周波雑音を有している。図４は４個の例示心電計
信号（リード）を示し、心電計信号４００，４０２，４
０４，４０６は全てＱＲＳ群を示している。
【００８７】図示してあるのは、それぞれ多くの高周波
雑音成分及び低周波雑音成分を有する心電計信号４０
２，４０６である。すなわち、心電計信号４０２がきれ
いな心電計波形成分４０１及び高周波雑音成分４０３で
あることを示している。また、図示はは、心電計信号４
０６が他のきれいな心電計波形成分４０５及び低周波成
分４０７を示している。
【００８８】上記き説明のように、活動関数は心電計信
号の絶対値を利用している。しかし、図４では、心電計
信号４００，４０２，４０４，４０６の修正形態を図示
してないが、それらが存在しているものである。
【００８９】上記の説明のように、心電計信号が過度な
高周波雑音又は過度な低周波雑音を有する図４に示した
特殊事例では、それぞれのリードに対するＳＱＩに０が
設定されている。上記をもって、すでに説明したよう
に、一つ以上の心電計信号に過度に雑音が多いことが判
明している図４の特殊事例では、このような雑音の多い
それぞれのリードに対するＳＱＩに０が設定される。続
いて、雑音が多くないことが判明している全てのリード
のＳＱＩを上記のＳＱＩ数式を使用して算出する。その
後に、上記の標準重み付け係数の数式を使用して全ての
リードに重み付け係数を設定する。図示した標準重み付
け係数の式を検討すると、リードのＳＱＩが０であれ
ば、そのリードに対する重み付け係数も０であることが
判明する。
【００９０】図４は心電計信号４０２が過度の高周波成
分を有していることが示されている。これは、上記で説
明したように、高周波雑音指数を移動平均することによ
り決定される。したがって、図４は過度の高周波雑音成
分４０３のために過度に雑音の多い心電計信号４０３
に、（このように過度に雑音の多いリードのＳＱＩは０
と規定されている。かつ、このように重み付けの式を使
用すると０の重み付け係数を生ずるので）０の重み付け
係数が設定されていることを示している。
【００９１】また、図４は、過度の低周波雑音成分４０
７のため、過度に雑音の多い心電計信号４０６に、（こ
のように過度に雑音の多いリードのＳＱＩは０と規定さ
れている。かつ、このように重み付けの式を使用すると
０の重み付け係数を生ずるので）０の重み付け係数が設
定されていることを示している。
【００９２】更に、図４は、０のＳＱＩが設定されてい
ない二つのリードが残っており、その残りのそれぞれの
リードは同等に良い信号であるＳＱＩ指数を有している
ので、ＱＲＳ群４００，４０６には、それぞれ１／２の
重み付け係数が設定されている。但し、ここに示した重
み付け係数は単に関連する概念を示すのに提示したもの
であって、図示した信号から得られた実際の重み付け係
数ではない。
【００９３】図示したＱＲＳ群（すなわち、心電計信
号）は全て、（上記の説明での活動関数により）そのそ
れぞれの重み付け係数が掛けられて加算され、複合ＱＲ
Ｓ群４２０を生成する。複合ＱＲＳ群（すなわち、活動
関数）から、図４に示した「強い」信号により、等しい
部分に、複合ＱＲＳ群が形成されていること、及び、過
度に雑音の多い信号の関与が無視されていることが判明
する。したがって、図４は、過度に雑音の多い信号と考
えられる信号の活動関数への関係を排除することを示し
ている。
【００９４】図５は、本発明の心臓波形特性表示方法及
びプロセスを示す高レベルの論理流れ図である。このプ
ロセスにあって、方法ステップ５００は、プロセスの開
始を示している。方法ステップ５０２は、Ｎ個の心電計
リード（心電計信号）の取得を示している。
【００９５】方法ステップ５０４は、Ｎ個の全ての信号
が「弱い」信号（すなわち、図２に関連して上記で説明
したような）を示しているか否かの決定を示している。
全ての信号が「弱い」信号であると決定されれば、方法
ステップ５０６に進み、Ｎ個のリード全てのＳＱＩを０
に設定する。その後、方法ステップ５０８に進み、それ
ぞれのリードについて決定されたＳＱＩ値を使用して活
動関数を計算する。
【００９６】全ての信号が「弱い」信号ではないことが
決定された場合、全てのリードに対するＳＱＩを算出す
る方法ステップ５０９に進む。その後、少なくとも一つ
のリードが「良い」信号を有しているか否か、及び、一
つ以上のリードが（例えば、図３に関連して説明したよ
うな）「弱い」信号を示しているか否かの方法ステップ
５１０に進む。方法ステップ５１０の規準に合わないこ
とが決定されれた場合、方法ステップ５１４に進む。
（方法ステップ５０９で計算したＳＱＩ値を使用するこ
とにより）少なくとも一つのリードが「良い」信号を有
し、（それらのリードの振幅を使用して）一つ以上のリ
ードが「弱い」ことが決定されれば、方法ステップ５１
２に進み、全ての「弱い」リードのＳＱＩが０に設定さ
れることを示す。その後、方法ステップ５１４に進む。
【００９７】方法ステップ５１４は、一つ以上の他の場
合に「良い」又は「強い」信号が、リードに存在する高
周波雑音及び低周波雑音のいずれか又は双方の累積指数
の信号が信頼できないことを示すという点で、廃棄され
るか否かの決定を示す。一つ以上の他の場合に「良
い」、又は「強い」信号が雑音のため信頼できないこと
が決定されれば、これら過度に雑音の多いチャンネルに
対するＳＱＩが０に設定されることを示す方法ステップ
５１６に進む。その後、それぞれのリードについて決定
されたＳＱＩ値を使用する活動関数の計算を示す方法ス
テップ５０８に進む。
【００９８】他の場合に「良い」又は「強い」信号に元
来雑音が多くないことが決定された場合、それぞれのリ
ードについて決定されたＳＱＩ値を使用する活動関数の
計算を示す方法ステップ５０８に進む。
【００９９】方法ステップ５０８ に続いてプロセス
は、分析すべき更に多数のデータが存在するか否かの方
法ステップ５１７に進む。分析すべき更に多数のデータ
が存在すれば、方法ステップ５０２に戻る。しかし、分
析すべき更に多数のデータが存在しなければ、方法ステ
ップ５１８に進んで終了となる。
【０１００】図６は、本発明を実施する心臓波形特性表
示方法及びそのシステムを示す部分概略ブロック図であ
る。図６は先に示した等式及び図に関連して前に既に説
明した１からＮまでのＥＫＧリード（又はチャンネル、
又は信号）を示している。図示は、それぞれのリードが
対応する信号良否測定ブロック５５０に送出されること
を示している。それぞれの信号良否測定ブロック５５０
は、多種の公式及び特殊の場合のＳＱＩ調査資料を使用
してそれぞれのリードに対するＳＱＩを計算する。当業
者は、それぞれの信号良否測定ブロック５５０としてソ
フトウェア、ハードウェア、ファームウェア、又は、こ
れらの組合せのような、種々の態様が存在することを認
識できるものである。
【０１０１】図６は、それぞれの信号良否測定ブロック
５５０によりそれぞれのリードに対するＳＱＩを決定し
たのに続き、このようなＳＱＩを重み付け係数決定ブロ
ック５６０に送出することを示している。重み付け係数
決定ブロック５６０は、多種の公式及び特殊の場合のＳ
ＱＩ調査資料を使用してそれぞれのリード（又はチャン
ネル、又は信号）に対する重み付け係数を算出する。当
業者は、重み付け係数決定ブロック５６０をソフトウェ
ア、ハードウェア、ファームウェア、又は、これらの組
合せのような、種々の態様が存在することを認識できる
ものである。
【０１０２】それぞれのリード（又はチャンネル、又は
信号）に適切な重み付け係数の決定に続いて、このよう
な重み付け係数を、ＥＫＧリード１からＮまでに対する
活動関数の出力を計算して発生するために、上記で説明
した活動関数に関連して、その重み付け係数を利用する
ＱＲＳ検出信号発生ブロック５７０に送出する。この出
力を図６でＱＲＳ検出信号５８０と記載している。当業
者は、ＱＲＳ検出信号発生ブロック５７０が、ソフトウ
ェア、ハードウェア、ファームウェア、又は、これらの
組合せのような、種々の態様が存在することを認識でき
るものである。
【０１０３】図７は本発明の例示実施形態の心臓波形特
性表示方法及びそのシステムに従って利用することが出
来るデータ処理システムの絵画的表現を示す図である。
上記の実施形態での処理を図７に示したデータ処理シ
ステムで実行できる。この例はシステムユニット６２
２、ビデオ表示端末６２４、キーボード６２６及びマウ
ス６２８を備えたコンピュータ６２０を有している。コ
ンピュータ６２０は、市販されているメインフレームコ
ンピュータ、ミニコンピュータ、又は、マイクロコンピ
ュータなどの、どのようなコンピュータにも適用でき
る。
【０１０４】図８は、本発明の例示実施形態の方法及び
システムに従って利用できる代表的ハードウェア環境を
説明するための図である。この例は上記の実施形態を実
行するコンピュータ６２０の構成を示している。システ
ムユニット６２２は、通常のマイクロプロセッサのよう
な、中央処理装置（「ＣＰＵ」）７３１及びシステムバ
ス７３２により相互接続されている多数の他の装置を備
えている。コンピュータ６２０は、読出し書込み記憶装
置（「ＲＡＭ」）７３４、読出し専用記憶装置（「ＲＯ
Ｍ」）７３６、システムバス７３２をビデオ表示端末６
２４に接続するための表示アダプタ７３７及び周辺装置
（例えば、ディスク・テープ駆動機構７３３）をシステ
ムバス７３２に接続するためのＩ／Ｏアダプタ７３９を
備えている。ビデオ表示端末６２４は、コンピュータ６
２０の出力を画面表示するものであり、コンピュータハ
ードウェアの分野で周知のＣＲＴを用いた表示装置であ
る。しかし、携帯式又はノートブック型コンピュータで
は、ビデオ表示端末６２４をＬＣＤを用いた扁平パネル
表示装置、又は、気体プラズマを用いた扁平パネル表示
装置に置換できる。コンピュータ６２０は更に、キーボ
ード６２６、マウス６２８、スピーカ７４６、マイクロ
ホン７４８、及び／又はタッチスクリーン装置（図示せ
ず）のような他のユーザインターフェース装置をシステ
ムバス７３２に接続するためのユーザインターフェース
アダプタを備えている。通信アダプタ７４９は、コンピ
ュータ６２０をデータ処理ネットワークに接続するため
のものである。
【０１０５】ＲＡＭ７３４、ＲＯＭ７３６、磁気ディス
ケット、磁気テープ、又は、光学ディスク（最後の三つ
は、ディスク・テープ駆動装置７３３に設置される）の
ような、いずれの適切な読取りが可能な媒体も利用でき
る。ユーザインターフェースもＣＰＵ７３１を指示でき
る。また、タッチスクリーン技術又は音声制御装置のよ
うな他の技術もＣＰＵ７３１と関連して使用することが
出来る。コンピュータ６２０は、コンピュータ記憶装置
７５０に装備される制御プログラム７５１を実行する。
制御プログラム７５１には、ＣＰＵ７３１で実行される
場合に、ここで説明したように図５の論理流れ図に、及
び図１、図２、図３、図４及び図６の部分概略図に示し
た動作を行なう制御プログラムが格納されている。
【０１０６】当業者は、図８に示したハードウェアを特
定の用途について変更できることを認識できるものであ
る。例えば、光学ディスク媒体、オーディオアダプタ、
又はコンピュータハードウェアの分野で周知のＰＡＬ又
はＥＰＲＯＭプログラミング装置（チッププログラミン
グ装置）などのような他の周辺装置をハードウェアに代
えて利用できる。
【０１０７】このように例示実施形態をコンピュータに
関連して説明した。当業者は、本発明の例示実施形態が
多様な形態でプログラム製品として配布され、また、例
示実施形態が記憶保持媒体によって本発明の関係人に関
係なく配付されることは重要な問題である。このような
記憶保持媒体としては、フロッピーディスクのような記
録可能な媒体、ハードディスク駆動機構、ＣＤ−ＲＯＭ
及びディジタル通信リンク及びアナログ通信リンクのよ
うな伝送形式の媒体がある。
【０１０８】以上、例示的な実施形態を特に図示し説明
してきたが、当業者は種々の実施形態及び、その細目に
関して、例示実施形態の技術思想及び技術的な範囲から
逸脱することなく、多種の変更を行なうことが出来るこ
とを理解できるものである。
【０１０９】以下に本発明の実施の形態を要約する。 １．心電計測定値から心臓の機能を示す信号の良否特性
を表示する心臓波形特性表示方法であって、心臓の機能
を示す少なくとも一つの信号を測定するステップ（５０
２）と、カージナル数の点で測定した信号に対するサン
プリングを行うステップ（５０２）と、前記カージナル
数のサンプリング点での最大サンプリング値を決定する
ステップと、前記決定した最大サンプリング値をサンプ
リングした信号を関連する除数で除算することによっ
て、信号の良否指数を算出するステップ（５０４，５０
６，５０９，５１０，５１２，５１４，５１６）と、を
有する心臓波形特性表示方法。
【０１１０】２．心臓機能の波形表現を心電計測定値か
ら得る場合に、心臓機能の波形表現の特性を高精度に表
示する方法であって、心臓の機能を示す少なくとも一つ
の信号を測定するステップ（５０２）と、心臓の機能を
示す前記測定した少なくとも一つの信号について、信号
の良否の指数を計算するステップ（５０４，５０６，５
０９，５１０，５１２，５１４，５１６）と、心臓の機
能を示す前記測定した少なくとも一つの信号について、
前記計算した信号の良否の指数に応じて少なくとも一つ
の重み付け係数を計算するステップと、心臓の機能を示
す前記測定した少なくとも一つの信号に心臓の機能を示
す前記測定した少なくとも一つの信号に対応する前記計
算した重み付け係数を乗じて波形表現を生成するステッ
プ（５０８）と、を有する心臓波形特性表示方法。
【０１１１】３．前記信号の良否の指数を計算する前記
ステップは更に、心臓の機能を示す少なくとも一つの信
号を測定するステップ（５０２）と、カージナル数の点
で前記測定した少なくとも一つの信号をサンプリングす
るステップ（５０２）と、前記サンプリングした少なく
とも一つの信号について、前記カージナル数のサンプリ
ング点での最大サンプリング値を決定するステップと、
前記サンプリングした少なくとも一つの信号について、
前記決定された最大値を前記サンプリングした少なくと
も一つの信号に関連する除数で除算して信号の良否の指
数を算出するステップ（５０４，５０６，５０９，５１
０，５１２，５１４，５１６）と、を有する上記２記載
の心臓波形特性表示方法。
【０１１２】４．前記信号の良否の指数を計算する前記
ステップは更に、前記サンプリングした少なくとも一つ
の信号について、前記カージナル数のサンプリング点で
取得したサンプリング値を集計して、前記サンプリング
した少なくとも一つの信号に関連する前記除数を計算す
るステップを有する上記１又は３記載の心臓波形特性表
示方法。
【０１１３】５．前記少なくとも一つの重み付け係数を
算出する前記ステップは更に、前記算出した信号の良否
の指数を補正するステップと、心臓の機能を示す前記少
なくとも一つの測定した信号について、前記計算した信
号の良否の指数を集計するステップと、前記少なくとも
一つの測定した信号について、前記計算した信号の良否
の指数の総計を０と比較するステップと、前記比較する
ステップが前記計算した信号の良否の指数の前記総計が
０に等しい場合に、前記少なくとも一つの重み付け係数
を全て、前記測定した少なくとも一つの信号の総計数で
除算した値に等しくなるように計算するステップと、前
記比較するステップが前記計算した信号の良否の指数の
前記総計が前記少なくとも一つの信号について０に等し
くない場合に、前記少なくとも一つの測定した信号につ
いて、前記少なくとも一つの重み付け係数を、前記少な
くとも一つの信号について、前記信号の良否の指数を、
前記集計ステップの結果で除算するステップと、を有す
る上記２記載の心臓波形特性表示方法。
【０１１４】６．心臓の機能を示す信号を心電計測定値
から得る場合に、心臓の機能を示す信号の良否特性を表
示するシステムであって、心臓の機能を示す少なくとも
一つの信号を測定する手段と、カージナル数の点で前記
測定値を少なくとも一つの信号をサンプリングする手段
と、前記サンプリングした少なくとも一つの信号につい
て、前記カージナル数のサンプリング点で生じた最大サ
ンプリング値を決定する手段と、前記サンプリングした
少なくとも一つの信号について、前記決定された最大値
を前記サンプリングした少なくとも一つの信号に関連す
る除数で除算することにより、信号の良否の指数を計算
する手段（５５０）と、を備える心臓波形特性表示シス
テム。
【０１１５】７．心臓機能の波形表現を心電計測定値か
ら得る場合に、心臓機能の波形表現の特性を高精度に表
示するシステムであって、心臓の機能を示す少なくとも
一つの信号を測定する手段と、心臓の機能を示す前記測
定した少なくとも一つの信号について、信号の良否の指
数を計算する手段（５５０）と、心臓の機能を示す前記
測定した少なくとも一つの信号について、前記計算した
信号の良否の指数に応じて少なくとも一つの重み付け係
数を計算する手段（５６０）と、心臓の機能を示す前記
測定した少なくとも一つの信号に心臓の機能を示す前記
測定した少なくとも一つの信号に対応する前記計算した
重み付け係数を乗じた波形表現を生成する手段（５７
０）と、を備える心臓波形特性表示システム。
【０１１６】８．前記信号の良否の指数を計算する前記
手段は更に、心臓の機能を示す少なくとも一つの信号を
測定する手段と、カージナル数の点で前記測定した少な
くとも一つの信号をサンプリングする手段と、前記サン
プリングした少なくとも一つの信号について、前記カー
ジナル数のサンプリング点にわたり生じた最大サンプリ
ング値を決定する手段と、前記サンプリングした少なく
とも一つの信号について、前記決定された最大値を前記
サンプリングした少なくとも一つの信号に関連する除数
で除算することにより、信号の良否の指数を計算する手
段と、を備える上記７記載の心臓波形特性表示システ
ム。
【０１１７】９．前記信号の良否の指数を計算する前記
手段は更に、前記サンプリングした少なくとも一つの信
号について、前記カージナル数のサンプリング点で取得
したそれぞれのサンプリング値を集計する手段と、前記
少なくとも一つのサンプリングした信号に関連する前記
除数を、前記集計ステップの結果に前記サンプリングが
行なわれたカージナル点の数の逆数を乗ずることにより
計算する手段と、を備える上記６又は８記載の心臓波形
特性表示システム。
【０１１８】１０．前記少なくとも一つの重み付け係数
を計算する前記手段は更に、特別の場合の事例におい
て、前記計算した信号の良否の指数を調節する手段と、
心臓の機能を示す前記少なくとも一つの測定した信号に
ついて、前記計算した信号の良否の指数を集計する手段
と、前記少なくとも一つの測定した信号について、前記
計算した信号の良否の指数の総計を０と比較する手段
と、前記比較する手段が前記計算した信号の良否の指数
の前記総計が０に等しいことを示したことに応答して、
前記少なくとも一つの重み付け係数を全て、数１を前記
測定した少なくとも一つの信号の総計数で除算した値に
等しくなるように計算する手段と、前記比較する手段が
前記計算した信号の良否の指数の前記総計が前記少なく
とも一つの信号について０に等しくないことを示したこ
とに応答して、前記少なくとも一つの測定した信号につ
いて、前記少なくとも一つの重み付け係数を、前記少な
くとも一つの信号についての前記信号の良否の指数を、
前記集計手段の結果で除算する手段と、を備える上記７
記載の心臓波形特性表示システム。
【０１１９】
【発明の効果】以上の説明から明らかなように、本発明
の心臓波形特性表示技術によれば、心臓の機能を示す一
つ以上の信号のカージナル数の点でサンプリングし、か
つ、最大サンプリング値を決定し、サンプリングした一
つ以上の信号のそれぞれの関連する除数で除算して良否
の指数を算出している。この結果、波形表現を心電計測
定値から得る場合に、ＱＲＳ群のような心臓機能の波形
表現を高精度に特性表示できるようになる。
【図面の簡単な説明】
【図１】本発明の心臓波形特性表示方法及びそのシステ
ムにおける実施形態での活動関数に関係する概念を説明
するための図である。
【図２】本発明の全ての信号が弱い特殊事例に関係する
概念を説明するための部分概要を示す図である。
【図３】本発明の第２特殊事例に関係する概念を説明す
るための部分概要を示す図である。
【図４】本発明の第３特殊事例に関係する概念を説明す
るための部分概要を示す図である。
【図５】本発明の心臓波形特性表示方法及びそのシステ
ムを示す高レベル論理流れ図である。
【図６】本発明を実施するための心臓波形特性表示方法
及びそのシステムを示す部分概略ブロック図である。
【図７】本発明の実施形態の心臓波形特性表示方法及び
そのシステムに従って使用されるデータ処理システムの
外観構成を示す斜視図である。
【図８】本発明の実施形態の心臓波形特性表示方法及び
そのシステムに従って使用される代表的なハードウェア
構成を示すブロック図である。
【符号の説明】
１００ 理想ＱＲＳ群 １０２，１０４，１０６ 理想に届かないＱＲＳ群 １０８，１１０，１１２，１１４ 修正ＱＲＳ群 １２０ 活動関数

Claims (1)

    【特許請求の範囲】
  1. 【請求項１】 心電計測定値から心臓の機能を示す信号
    の良否特性を表示する心臓波形特性表示方法であって、 心臓の機能を示す少なくとも一つの信号を測定するステ
    ップ（５０２）と、カージナル数の点で測定した信号に
    対するサンプリングを行うステップ（５０２）と、 前記カージナル数のサンプリング点での最大サンプリン
    グ値を決定するステップと、 前記決定した最大サンプリング値をサンプリングした信
    号を関連する除数で除算することによって、信号の良否
    指数を算出するステップ（５０４，５０６，５０９，５
    １０，５１２，５１４，５１６）と、を有することを特
    徴とする心臓波形特性表示方法。
JP10198968A 1997-07-15 1998-07-14 心臓波形特性表示方法 Pending JPH1170089A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US08/893,115 US5827196A (en) 1997-07-15 1997-07-15 Method and system for providing characterizations of waveform representations of heart function
US08/893-115 1997-07-15

Publications (1)

Publication Number Publication Date
JPH1170089A true JPH1170089A (ja) 1999-03-16

Family

ID=25401057

Family Applications (1)

Application Number Title Priority Date Filing Date
JP10198968A Pending JPH1170089A (ja) 1997-07-15 1998-07-14 心臓波形特性表示方法

Country Status (3)

Country Link
US (1) US5827196A (ja)
JP (1) JPH1170089A (ja)
DE (1) DE19827697B4 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005253609A (ja) * 2004-03-10 2005-09-22 Fukushima Prefecture 生体情報処理システム
JP2014528342A (ja) * 2011-10-12 2014-10-27 コーニンクレッカ フィリップス エヌ ヴェ リアルタイムｅｃｇにおいて心電図のｓｔ部レベルを全自動で測定する方法およびシステム
JP2016506779A (ja) * 2013-12-20 2016-03-07 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． Ｅｃｇデータにおけるｑｒｓ群の発生を決定する装置及び方法
CN107580471A (zh) * 2015-04-02 2018-01-12 微软技术许可有限责任公司 可穿戴脉搏感测设备信号质量估计
JP2018528812A (ja) * 2015-08-25 2018-10-04 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． Ｅｃｇリード信号の高／低周波信号品質評価

Families Citing this family (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU4214199A (en) 1998-06-03 1999-12-20 Masimo Corporation Stereo pulse oximeter
US6452679B1 (en) 1999-12-29 2002-09-17 Kimberly-Clark Worldwide, Inc. Method and apparatus for controlling the manufacturing quality of a moving web
US6438405B1 (en) 2000-04-28 2002-08-20 Koninklijke Philips Electronics, N.V. Imaging safety device
US6766189B2 (en) 2001-03-30 2004-07-20 Cardiac Pacemakers, Inc. Method and apparatus for predicting acute response to cardiac resynchronization therapy
US6705999B2 (en) * 2001-03-30 2004-03-16 Guidant Corporation Method and apparatus for determining the coronary sinus vein branch accessed by a coronary sinus lead
US6993389B2 (en) 2001-03-30 2006-01-31 Cardiac Pacemakers, Inc. Identifying heart failure patients suitable for resynchronization therapy using QRS complex width from an intracardiac electrogram
US7587237B2 (en) * 2004-02-02 2009-09-08 Cardionet, Inc. Biological signal management
US7283864B2 (en) * 2005-02-10 2007-10-16 Cardiac Pacemakers, Inc. Method and apparatus for identifying patients with wide QRS complexes
US8942795B2 (en) * 2005-03-31 2015-01-27 Medtronic, Inc. Implantable medical device with real time T-wave oversensing detection
US8369944B2 (en) 2007-06-06 2013-02-05 Zoll Medical Corporation Wearable defibrillator with audio input/output
TWI360971B (en) * 2008-09-09 2012-03-21 Realtek Semiconductor Corp An receiving apparatus and method for receiving pa
CN103002800B (zh) 2010-05-18 2015-08-26 佐尔医药公司 具有多个感测电极的可穿戴移动医疗装置
CN102038497B (zh) * 2010-12-02 2012-07-18 广东宝莱特医用科技股份有限公司 一种心电信号噪声分析方法
WO2012135062A1 (en) * 2011-03-25 2012-10-04 Zoll Medical Corporation Selection of optimal channel for rate determination
US8600486B2 (en) 2011-03-25 2013-12-03 Zoll Medical Corporation Method of detecting signal clipping in a wearable ambulatory medical device
EP4152340A1 (en) 2011-03-25 2023-03-22 Zoll Medical Corporation System and method for adapting alarms in a wearable medical device
US9684767B2 (en) 2011-03-25 2017-06-20 Zoll Medical Corporation System and method for adapting alarms in a wearable medical device
JP2015510780A (ja) 2012-03-02 2015-04-13 ゾール メディカル コーポレイションＺＯＬＬ Ｍｅｄｉｃａｌ Ｃｏｒｐｏｒａｔｉｏｎ 着用式医療監視および／または処置装置を構成するためのシステムおよび方法
BR112014029588A2 (pt) 2012-05-31 2017-06-27 Zoll Medical Corp sistema e métodos para detectar distúrbios de saúde
WO2014020484A2 (en) * 2012-08-01 2014-02-06 Koninklijke Philips N.V. A method and system to identify motion artifacts and improve reliability of measurements and alarms in photoplethysmographic measurements
US9579516B2 (en) 2013-06-28 2017-02-28 Zoll Medical Corporation Systems and methods of delivering therapy using an ambulatory medical device
US9750423B2 (en) 2013-08-13 2017-09-05 Koninklijke Philips N.V. Method and display for long term physiological signal quality indication
CN105050493B (zh) * 2013-12-20 2016-08-17 皇家飞利浦有限公司 用于确定ecg数据中的qrs复合波的出现的装置和方法
WO2015123198A1 (en) 2014-02-12 2015-08-20 Zoll Medical Corporation System and method for adapting alarms in a wearable medical device
WO2016100906A1 (en) 2014-12-18 2016-06-23 Zoll Medical Corporation Pacing device with acoustic sensor
WO2016149583A1 (en) 2015-03-18 2016-09-22 Zoll Medical Corporation Medical device with acoustic sensor
CN105725966A (zh) * 2016-01-27 2016-07-06 浙江铭众科技有限公司 一种基于神经网络模型的心电信号质量判别方法
CN105726013A (zh) * 2016-01-27 2016-07-06 浙江铭众科技有限公司 一种具有心电信号质量判别功能的心电监测系统
US11617538B2 (en) 2016-03-14 2023-04-04 Zoll Medical Corporation Proximity based processing systems and methods
CN109770889B (zh) * 2017-11-15 2022-03-11 深圳市理邦精密仪器股份有限公司 心电数据选段方法和装置
US10960213B2 (en) 2018-03-12 2021-03-30 Zoll Medical Corporation Verification of cardiac arrhythmia prior to therapeutic stimulation
US11568984B2 (en) 2018-09-28 2023-01-31 Zoll Medical Corporation Systems and methods for device inventory management and tracking

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4751931A (en) * 1986-09-22 1988-06-21 Allegheny-Singer Research Institute Method and apparatus for determining his-purkinje activity
US4987901A (en) * 1989-06-28 1991-01-29 Kunig Horst E Method and apparatus for selecting a physiologically standardized sensor of a multi-sensor electrocardiogram sensor set
AU5205493A (en) * 1992-12-01 1994-06-16 Siemens Aktiengesellschaft Cardiac event detection in implantable medical devices
US5685315A (en) * 1992-12-01 1997-11-11 Pacesetter, Inc. Cardiac arrhythmia detection system for an implantable stimulation device

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005253609A (ja) * 2004-03-10 2005-09-22 Fukushima Prefecture 生体情報処理システム
JP4633373B2 (ja) * 2004-03-10 2011-02-16 公立大学法人会津大学 生体情報処理システム
JP2014528342A (ja) * 2011-10-12 2014-10-27 コーニンクレッカ フィリップス エヌ ヴェ リアルタイムｅｃｇにおいて心電図のｓｔ部レベルを全自動で測定する方法およびシステム
JP2016506779A (ja) * 2013-12-20 2016-03-07 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． Ｅｃｇデータにおけるｑｒｓ群の発生を決定する装置及び方法
US10058260B2 (en) 2013-12-20 2018-08-28 Koninklijke Philips N.V. Apparatus and method for determining the occurrance of a QRS complex in ECG data
US10285613B2 (en) 2013-12-20 2019-05-14 Koninklijke Philips N.V. Apparatus and method for determining the occurrence of a QRS complex in ECG data
CN107580471A (zh) * 2015-04-02 2018-01-12 微软技术许可有限责任公司 可穿戴脉搏感测设备信号质量估计
US10765331B2 (en) 2015-04-02 2020-09-08 Microsoft Technology Licensing, Llc Wearable pulse sensing device signal quality estimation
CN107580471B (zh) * 2015-04-02 2020-12-04 微软技术许可有限责任公司 可穿戴脉搏感测设备信号质量估计
JP2018528812A (ja) * 2015-08-25 2018-10-04 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． Ｅｃｇリード信号の高／低周波信号品質評価

Also Published As

Publication number Publication date
DE19827697A1 (de) 1999-02-11
US5827196A (en) 1998-10-27
DE19827697B4 (de) 2005-04-07

Similar Documents

Publication Publication Date Title
JPH1170089A (ja) 心臓波形特性表示方法
JP4768095B2 (ja) 心臓機能を表す信号の特性を決定するシステム
EP3698707B1 (en) Electrocardiogram information dynamic monitoring system, computer program and computer readable storage medium
JP6381444B2 (ja) リアルタイムｅｃｇにおいて心電図のｓｔ部レベルを全自動で測定する方法およびシステム
Xue et al. Algorithms for computerized QT analysis
EP3692904A1 (en) Method and device for self-learning dynamic electrocardiography analysis employing artificial intelligence
EP3692903A1 (en) Artificial intelligence self-learning-based automatic electrocardiography analysis method and apparatus
EP3698709A1 (en) Electrocardiogram information processing method and electrocardiogram workstation system
CN106815570B (zh) 一种基于动态模式识别的心电信号st-t段识别方法
Agostinelli et al. Noninvasive fetal electrocardiography Part I: Pan-tompkins' algorithm adaptation to fetal R-peak identification
EP3698708A1 (en) Ecg information processing method and ecg workstation
CN110693483A (zh) 一种动态心电图自动分析的方法
JPH01110344A (ja) 分娩中の胎児の監視
CN108024750A (zh) Ecg导联信号的高/低频信号质量评价
CN107811631A (zh) 心电信号质量评估方法
Gendelman et al. PhysioZoo ECG: Digital electrocardiography biomarkers to assess cardiac conduction
Zhang et al. An Improved Real‐Time R‐Wave Detection Efficient Algorithm in Exercise ECG Signal Analysis
Hermawan et al. Development of ECG signal interpretation software on Android 2.2
JP5933138B2 (ja) Ｅｃｇデータにおけるｑｒｓ群の発生を決定する装置及び方法
Palaniappan et al. R-peak detection using altered pan-Tompkins algorithm
Ostertag et al. Reconstructing ECG precordial leads from a reduced lead set using independent component analysis
Mengko et al. Design and implementation of 12 Lead ECG signals interpretation system
CN117338309B (zh) 一种身份识别方法及存储介质
Agrawal et al. Morphological Event Based Signal Quality Assessment of Electrocardiogram
CN111685754B (zh) 针对可穿戴ecg采集设备的心率计算方法及系统

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20050713

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20050713

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20070410

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20080812

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20090127