(57)【要約】【課題】　ＥＣＧ波形及び複合体の表示に遅延を導入せずに、ＥＣＧ波形及び複合体の保全性を維持しつつ、基線ドリフト及び低周波雑音を補償する。【解決手段】　表示領域を判定された遅延期間にまでさかのぼる最も近い時点で受信されたデータを含む波形データの第１のセグメントを表示する第１の部分（４０４）と、波形データのその他の部分を含む波形データの第２のセグメントを表示する第２の部分（４０２）とに配分する。第１の部分には、部分基線補正のために連続変化振幅レベル調整が適用され、第２の部分には、それ以上の振幅レベル調整を伴わない基線補正振幅調整が適用される。【選択図】　図１
2004065981.tif 000002
【特許請求の範囲】
【請求項１】
波形データを表示装置に表示する方法（３００）において、
　表示領域を第１の部分（４０４）と、前記第１の部分（４０４）にすぐ隣接する第２の部分（４０２）とに配分することから成り、
　前記第１の部分（４０４）は、判定された遅延周期にまでさかのぼる最も近い時点で受信されたデータから成る波形データの第１のセグメントを表示するために使用され、
　前記第２の部分（４０２）は、波形データのその他の部分から成る波形データの第２のセグメントを表示するために使用され、
　前記第１の部分（４０４）で表示されるデータには、部分基線補正のために連続変化振幅レベル調整が適用されており、前記第２の部分（４０２）で表示されるデータには、それ以上の振幅レベル調整を伴わない基線補正振幅調整が適用されている方法（３００）。
【請求項２】
前記基線調整は対称有限インパルス応答フィルタを介して適用される請求項１記載の方法（３００）。
【請求項３】
前記対称有限インパルス応答フィルタは約０．５秒から約３．０秒の時間遅延で所定のデータサンプルを中心としている請求項２記載の方法（３００）。
【請求項４】
波形データは、波形データが初めに前記第１の部分（４０４）に現れ、その後、前記基線振幅調整を受けながら前記第１の部分（４０４）を通ってスクローリングし、その後、それ以上の基線調整を適用されずに前記第２の部分（４０２）に現れ、前記第２の部分（４０２）を通ってスクローリングするように前記表示領域に沿ってスクローリングされる請求項１記載の方法（３００）。
【請求項５】
前記表示領域の前記第１の部分（４０４）がワイパーバーの一方の縁部に沿ってワイパーバーに従い、それに隣接して動くように、前記表示領域に沿ってワイパーバーがスクローリングされる請求項１記載の方法（３００）。
【請求項６】
順次波形データサンプルをフィルタ補正し且つ表示する方法（３００）において、
　格納されている未補正データサンプルのシーケンスをシフトすること（３０２）と、
　新たな未補正データサンプルを受信し且つ格納すること（３０４）と、
　前記格納されている未補正データサンプル及び前記新たな未補正データサンプルを使用して基線推定補正を計算すること（３０６）と、
　格納されている補正済みデータサンプルのシーケンスをシフトし、その後、前記格納されている未補正データサンプルのうちの特定の１つに前記基線推定補正を適用することにより新たな補正済みデータサンプルを判定すること（３０８）と、
　前記格納されている未補正データサンプルのうちの前記特定の１つの後に続いて格納されている前記未補正データサンプルの各々と、前記新たな未補正データサンプルとに前記基線推定補正を適用することにより、一時表示データサンプルのシーケンスを作成すること（３１２）と、
　前記補正済みデータサンプルのシーケンス、前記新たな補正済みデータサンプル及び前記一時表示データサンプルのシーケンスを表示すること（３１４）から成る方法（３００）。
【請求項７】
前記基線推定補正は対称有限インパルス応答フィルタにより実現される請求項６記載の方法（３００）。
【請求項８】
前記格納されている未補正データサンプルのうちの前記特定の１つは前記格納されている未補正データサンプルの中の中央に位置している請求項７記載の方法（３００）。
【請求項９】
前記基線推定補正は、前記新たな未補正データサンプルを受信し且つ格納することに関して約０．５秒から約３．０秒の遅延で適用される請求項７記載の方法（３００）。
【請求項１０】
前記波形データサンプルは心電図データを表現する請求項６記載の方法（３００）。
【請求項１１】
被検体からＥＣＧ信号を検出する一組の電極（１１０）と、
　前記一組の電極（１１０）により検出された前記ＥＣＧ信号を調整する信号調整回路（１２５）と、
　前記信号調整回路（１２５）からの調整済み信号を処理するプロセッサ（１３０）と、
　前記プロセッサ（１３０）により生成されたＥＣＧ波形データを表示する表示装置（１３５）とを具備し、前記表示装置（１３５）は、
第１の部分（４０４）と、前記第１の部分（４０４）にすぐ隣接する第２の部分（４０２）とを有する表示領域を更に具備し、
　前記第１の部分（４０４）は、判定された遅延周期にまでさかのぼる最も近い時点で受信されたデータから成る波形データの第１のセグメントを表示するために使用され、
　前記第２の部分（４０２）は、波形データのその他の部分から成る波形データの第２のセグメントを表示するために使用され、
　前記第１の部分（４０４）で表示されるデータには、部分基線補正のために連続変化振幅レベル調整が適用されており、前記第２の部分（４０２）で表示されるデータには、それ以上の振幅レベル調整を伴わない基線補正振幅調整が適用されているＥＣＧシステム。
【請求項１２】
前記基線調整は対称有限インパルス応答フィルタを介して適用される請求項１１記載のＥＣＧシステム（１００）。
【請求項１３】
前記対称有限インパルス応答フィルタは約０．５秒から約３．０秒の時間遅延で所定のデータサンプルを中心としている請求項１２記載のＥＣＧシステム（１００）。
【請求項１４】
前記波形データは、前記波形データが初めに前記第１の部分（４０４）に現れ、その後、前記基線振幅調整を受けながら前記第１の部分（４０４）を通ってスクローリングし、その後、それ以上の基線調整を適用されずに前記第２の部分（４０２）に現れ、前記第２の部分（４０２）を通ってスクローリングするように前記表示領域に沿ってスクローリングされる請求項１１記載のＥＣＧシステム（１００）。
【請求項１５】
前記表示領域の前記第１の部分（４０４）がワイパーバーの一方の縁部に沿ってワイパーバーに従い、それに隣接して動くように、前記表示領域に沿ってワイパーバーがスクローリングされる請求項１１記載のＥＣＧシステム（１００）。
【請求項１６】
順次波形データサンプルをフィルタ補正し且つ表示するための機械読み取り可能なコンピュータプログラムコードと、
　コンピュータに、
　格納されている未補正データサンプルのシーケンスをシフトすること（３０２）と、
　新たな未補正データサンプルを受信し且つ格納すること（３０４）と、
　前記格納されている未補正データサンプル及び前記新たな未補正データサンプルを使用して基線推定補正を計算すること（３０６）と、
　格納されている補正済みデータサンプルのシーケンスをシフトし、その後、前記格納されている未補正データサンプルのうちの特定の１つに前記基線推定補正を適用することにより新たな補正済みデータサンプルを判定すること（３０８）と、
　前記格納されている未補正データサンプルのうちの前記特定の１つの後に続いて格納されている前記未補正データサンプルの各々と、前記新たな未補正データサンプルとに前記基線推定補正を適用することにより、一時表示データサンプルのシーケンスを作成すること（３１２）と、
　前記補正済みデータサンプルのシーケンス、前記新たな補正済みデータサンプル及び前記一時表示データサンプルのシーケンスを表示すること（３１４）から成る方法（３００）を実現させるための複数の命令とを具備する記憶媒体。
【請求項１７】
プロセッサ（１３０）に、
　格納されている未補正データサンプルのシーケンスをシフトすること（３０２）と、
　新たな未補正データサンプルを受信し且つ格納すること（３０４）と、
　前記格納されている未補正データサンプル及び前記新たな未補正データサンプルを使用して基線推定補正を計算すること（３０６）と、
　格納されている補正済みデータサンプルのシーケンスをシフトし、その後、前記格納されている未補正データサンプルのうちの特定の１つに前記基線推定補正を適用することにより新たな補正済みデータサンプルを判定すること（３０８）と、
　前記格納されている未補正データサンプルのうちの前記特定の１つの後に続いて格納されている前記未補正データサンプルの各々と、前記新たな未補正データサンプルとに前記基線推定補正を適用することにより、一時表示データサンプルのシーケンスを作成すること（３１２）と、
　前記補正済みデータサンプルのシーケンス、前記新たな補正済みデータサンプル及び前記一時表示データサンプルのシーケンスを表示すること（３１４）から成る方法（３００）を実現させるように構成されたコードを具備するコンピュータデータ信号。
【請求項１８】
前記基線推定補正は対称有限インパルス応答フィルタにより実現される請求項１６または１７記載の発明。
【請求項１９】
前記格納されている未補正データサンプルのうちの前記特定の１つは前記格納されている未補正データサンプルの中の中央に位置している請求項１８記載の発明。
【請求項２０】
前記基線推定補正は、前記新たな未補正データサンプルを受信し且つ格納することに関して約０．５秒から約３．０秒の遅延で適用される請求項１８記載の発明。
【請求項２１】
前記波形データサンプルは心電図データを表現する請求項１６または１７記載の発明。
【発明の詳細な説明】
【技術分野】
【０００１】
　本開示は一般にフィルタ補正波形データをリアルタイム表示する方法に関し、特に、フィルタ補正心電図データをリアルタイム表示する方法及び装置に関する。
【背景技術】
【０００２】
　心臓サイクルの心電図（ＥＣＧ）は患者の身体の表面に配置された複数対の感知電極で検出され、左右の心房及び心室の周期的なＰＱＲＳＴ電気活性化シーケンスにより特徴づけられる反復波形である。ＰＱＲＳＴシーケンスは順次起こる心房の脱分極及び収縮、並びにその後に起こる心室の脱分極及び収縮と関連しており、連続するＰＱＲＳＴ複合体が基線又は等電位領域により分離されている。
【０００３】
　図１に示すように、ＰＱＲＳＴ電気活性化シーケンスは、心房の脱分極及び収縮を示すＰ波で始まる。次に、心室の脱分極及び収縮を示すＱＲＳ複合体が続く。ＳＴセグメント時間遅延の終了時のＴ波は心室の再分極と関連している。１対の感知電極で全く損なわれていないＡ−Ｖ活性化が検出されるＰＱＲＳＴ電気活性化シーケンスの形状は正しく予測可能である。正常な心拍数の範囲内で順次起こるＰ波、Ｒ波及びＴ波は、通常、身体表面に装着され、脱分極波に関して正しく方向づけされている電極により記録される外部ＥＣＧを目視検査することにより容易に認識される。Ｐ波及びＲ波は、適切に配置された感知電極対と結合されたモニタ又は治療投与装置のセンス増幅器により容易に感知される。
【０００４】
　通常、ＥＣＧのＳＴセグメントの振幅は、感知電極対の配置に応じて、基線にごく近いか、又はＰＱＲＳＴシーケンスの間で感知される信号の等電位振幅にごく近い。心筋虚血の期間中、ＳＴセグメントの振幅は（心臓に対するＥＣＧ感知電極の位置どりに応じて）基線から高められるか、又は押し下げられる。それらのＳＴセグメント偏倚は目視検査により容易に認識できる。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
　しかし、ＥＣＧ信号は、通常、低周波雑音（例えば、心拍数より低い速度で起こる呼吸に起因する雑音など）の影響を受け、その結果、基線ドリフトを生じる。そのような影響によって、特に複数のＥＣＧ波形が同時に提示されるような表示装置においては、ＥＣＧ波形を読み取ることが困難になる。現在、基線ドリフトを強力に除去するフィルタ補正技法が存在しているが、これもＥＣＧ波形の一部（例えば、ＳＴセグメント）を歪ませる結果を招き且つ／又は表示の提示に遅延を導入する。フィルタ補正技法が歪みを最小にすること又は遅延を排除することを目標にするならば、通常はその代償として基線を強力に補正するという目標が失われてしまう。従って、ＥＣＧ波形及び複合体の表示に遅延を導入せずに、ＥＣＧ波形及び複合体の保全性を維持しつつ、基線ドリフト及び低周波雑音を補償することが可能であるのが望ましい。
【課題を解決するための手段】
【０００６】
　従来の技術の以上述べた欠点及び欠陥、並びにその他の欠点及び欠陥は、波形データを表示装置に表示する方法により克服又は軽減される。一実施例では、方法は、表示領域を第１の部分と、第１の部分にすぐ隣接する第２の部分とに配分することを含む。第１の部分は、判定された遅延周期にまでさかのぼる最も近い時点で受信されたデータから成る波形データの第１のセグメントを表示するために使用される。第２の部分は、波形データのその他の部分から成る波形データの第２のセグメントを表示するために使用される。第１の部分で表示されるデータには、部分基線補正のために連続変化振幅レベル調整が適用されており、一方、第２の部分で表示されるデータには、それ以上の振幅レベル調整を伴わない基線補正振幅調整が適用されている。
【０００７】
　別の面においては、順次波形データサンプルをフィルタ補正し且つ表示する方法は、格納されている未補正データサンプルのシーケンスをシフトすることと、次に、新たな未補正データサンプルを受信し且つ格納することを含む。格納されている未補正データサンプル及び新たな未補正データサンプルを使用して基線推定補正が計算される。次に、格納されている補正済みデータサンプルのシーケンスがシフトされ、格納されている未補正データサンプルのうちの特定の１つに基線推定補正を適用することにより新たな補正済みデータサンプルが判定される。格納されている未補正データサンプルのうちの特定の１つの後に続いて格納されている未補正データサンプルの各々と、新たな未補正データサンプルとに基線推定補正を適用することにより、一時表示データサンプルのシーケンスが作成される。次に、補正済みデータサンプルのシーケンス、新たな補正済みデータサンプル及び一時表示データサンプルのシーケンスがそれぞれ表示される。
【０００８】
　更に別の面においては、心電図（ＥＣＧ）システムは、被検体からＥＣＧ信号を検出する一組の電極と、一組の電極により検出されたＥＣＧ信号を調整する信号調整回路とを含む。プロセッサは、信号調整回路からの調整済み信号を処理するために使用される。加えて、プロセッサにより生成されたＥＣＧ波形データを表示する表示装置は、第１の部分と、第１の部分にすぐ隣接する第２の部分とを有する表示領域を更に含む。第１の部分は、判定された遅延周期にまでさかのぼる最も近い時点で受信されたデータから成る波形データの第１のセグメントを表示するために使用され、一方、第２の部分は、波形データのその他の部分から成る波形データの第２のセグメントを表示するために使用される。第１の部分で表示されるデータには、部分基線補正のために連続変化振幅レベル調整が適用されており、第２の部分で表示されるデータには、それ以上の振幅レベル調整を伴わない基線補正振幅調整が適用されている。
【発明を実施するための最良の形態】
【０００９】
　図面を参照すると、いくつかの図において同様の要素は同じ図中符号により示されている。
【００１０】
　ここでは、フィルタ補正心電図データをリアルタイム表示する方法及び装置が開示される。簡単に言えば、遅延領域の絶えず調整されているＤＣレベル部分を生成するために、遅延対称有限インパルス応答フィルタ（ＦＩＲ）が実現される。これにより、表示の第１の部分が波形データ全体のうちの、絶えずレベル調整され、部分的に補正された基線を有する、より最新の部分を特徴とする連続スクローリング表示が提供される。加えて、表示の第２の部分は、波形データ全体のうちの、一定振幅でスクロールする実質的に補正された基線を有するより以前の部分を特徴とする。この技法の結果、表示されるＥＣＧ波形のＳＴセグメントは歪みを生じないままである。
【００１１】
　一実施例では、対称ＦＩＲは、実際には、２秒をわずかに超える総幅を有する三角形インパルス応答であるように設計されている。従って、フィルタ補正の前に約１秒の遅延が実現され、それにより、表示領域の絶えず調整されているＤＣレベル部分（すなわち、第１の部分）は表示の右側にあって、約１秒の表示距離をカバーしている。しかし、この遅延は約０．５秒から約３．０秒までのより広い時間範囲にわたって実現されても良い。従来は、ＥＣＧシステムの場合、新たなデータは表示領域の右縁部に導入され、それより古いデータは左方向に向かって表示領域をスクローリングする。基線補正を実現するために、表示領域の第１の部分は画面再生時間のたびに完全にペインティングし直されるが、これに対し、表示領域の第２の部分に示されているデータは再生時間のたびに一定数の画素ずつ単純にスクローリングされる（すなわち、左へシフトされる）だけである。このような表示は、垂直再生周期ごとに画面全体をペインティングし直すことが可能である（毎秒、約７０回）最新のコンピュータビデオ表示装置では可能である。
【００１２】
　まず初めに図２を参照すると、本開示の一実施例を実施するのに適するＥＣＧシステム１００の一例が示されている。システム１００は、例えば、臨床環境において使用されるか、あるいはおそらくは患者の心臓の生理的負荷試験の間に使用されると考えられる。一実施例では、システム１００は複数の電極１１０から成る電極セット１０５を含み、それらの電極は標準形ＥＣＧ電極、すなわち、可能であれば患者の胴体の前部、側部及び後部の各領域をカバーするように装着される電極のアレイであっても良い。電極１１０は互いに別個に機能するが、可撓性バンド又はその他の構造を形成するように物理的に一体に張り付けられていても良い。
【００１３】
　加えて、システム１００は、電極をシステムコントローラ１２０に接続する一組のリード１１５を更に含む。システムコントローラ１２０は信号調整回路１２５と、プロセッサ１３０とを含む。信号調整回路１２５はリード１１５からのアナログ信号を入力として受信し、調整済みデジタル信号を出力としてプロセッサ１３０に供給する。プロセッサ１３０は調整済み信号を処理して、後にプロセッサ１３０に接続されている表示装置１３５及び／又はプリンタなどの出力装置１４０に提供される出力結果を生成する。負荷試験の用途に使用される場合、プロセッサ１３０は、傾斜及び歩行速度がプログラム可能であるランニングマシン１４５などの運動用装置をケーブル１５０を介して供給される制御信号によって更に制御しても良い。同様に、長時間（例えば、２４時間）にわたりリードからの信号を記録するために、アンビュラトリーシステムのオプションの記録装置１５５が使用されても良い。その場合、システムコントローラ１２０が記録データを処理できるように、記録装置１５５はシステムコントローラ１２０に接続される。
【００１４】
　先に述べた通り、従来の方法でフィルタ補正されたＥＣＧデータのリアルタイム表示に関連する問題が存在しており、それをかなり簡略化した例によって特定して例示する。そこで、図３（ａ）及び図３（ｂ）を参照すると、ＲＱＳ複合体の一般的形状を模倣した周期的三角形インパルスの一対の波形表示が示されている。図中に指示するように、上方の波形は方形波雑音により妨害された生信号を表し、下方の波形は０．０３７Ｈｚ高域フィルタによりフィルタ補正された後の同じ信号（同様に方形波雑音による妨害を受けている）を表す。図から見て取れるように、０．０３７Ｈｚフィルタは三角形パルスの保全性を維持しているが、波形をその基線レベルに戻すことに関して応答が遅くなっている。図３（ｂ）に示すように、方形波雑音が除去された場合、波形は基線の下まで降下しており、上方の基線まで戻るのが遅くなっている。
【００１５】
　これに対し、０．５９７Ｈｚ高域フィルタを使用した場合には、図３（ｃ）及び図３（ｄ）に示すように、基線への復帰ははるかに速くなる。しかし、このフィルタは、三角形パルスの立ち下がり端を立ち上がり端の始めのレベルより低くなるまで下げることにより、立ち下がり端を歪ませることにも気づくであろう。残念ながら、このような歪みの結果、実際のＥＣＧの読みが不正確に解釈されることも起こりうるであろう。従って、この特定の型の波形の強力なフィルタ補正は応答時間が遅い場合と同じように望ましくない。更に別の種類のフィルタ補正技法は、図３（ｅ）及び図３（ｆ）に示すような、放物線基線補正と呼ばれる技法である。図からわかる通り、放物線基線補正は０．５９７Ｈｚ高域フィルタほど多くの歪みを発生させないが、波形を基線に戻すことに関してもそれほど強力ではない。
【００１６】
　従って、本発明の一実施例によれば、心電図データなどのフィルタ補正データをリアルタイム表示する方法及び装置が開示される。そこで図４を参照すると、本発明のフィルタ補正及び表示技法の原理を示す概略図が示されている。ディスプレイバッファ２００は、図２の表示装置１３５などに表示されるべき最も近い時点のＥＣＧ波形を含む。単なる例ではあるが、ＥＣＧデータは毎秒２４０サンプルの速度の６秒ディスプレイバッファに合わせて構成されており、ある任意の瞬時に合計で１４４０のデータサンプルが表示されることになる。これより長い持続時間又は短い持続時間の波形信号が表示される場合、あるいは異なるサンプリング速度が使用される場合には、一度に表示されるデータサンプルの総数が異なってくることは自明である。
【００１７】
　図４に示すように、合わせて１４４０個のデータサンプルはＤtからＤt-1439Tと指示されている。最も新しく表示されたデータサンプルはＤtであり（ディスプレイバッファ２００の右端にある）、最も古くから表示されていたデータサンプルはＤt-1439Tである（ディスプレイバッファの左端にある）。尚、Ｔは約４．１７ミリ秒のサンプル周期を表し、ｔは現在時間を表す。この実施例では、表示は新たな入力サンプルが利用可能になるのと同じ速度で再生されると仮定する。
【００１８】
　ディスプレイバッファ２００に加えて、基線補正データを格納し且つシフトするためにフィルタ補正データバッファ２０２が使用される。フィルタ補正データバッファ２０２の基線補正データは、ディスプレイバッファ２００の対応する記憶場所へ直接に送り出される。表示装置１３５の第２の部分での表示及びスクローリングに使用されるのはこの基線補正データである。フィルタ補正データバッファ２０２のデータサンプルはＹt-256TからＹt-1439Tと指示されており、これは基線補正データのうちで最も古いデータが約１．０６秒の遅延の後に表示されることを意味している。フィルタ補正データバッファ２０２のすぐ右側に、Ｙt-255Tと指示されるサンプルデータを保持する記憶要素２０４がある。フィルタ補正データバッファ２０２に到達する前に生サンプルデータに対して重み付き平均基線調整が実行されるのはここである。最後に、ＦＩＲフィルタにより使用される最も近い時点の５１１個のフィルタ補正前のデータサンプルを格納し且つ（その後に）シフトするためにフィルタ補正前データバッファ２０６が使用される。従って、フィルタ補正前データバッファ２０６に含まれているデータサンプルはＸtからＸt-510Tと指示されており、Ｘtは最も新しい未補正入力サンプルを表し、Ｘt-510Tは最も古い未補正入力サンプルを表す。
【００１９】
　対称ＦＩＲフィルタは、画面再生のたびに平均基線推定を計算するために、サンプルＸt-255Tを中心とする５１１個のフィルタ補正前サンプルの全てを使用する。次に、計算された基線推定ＢをデータサンプルＸt-255Tから減算して、先に述べた通り、Ｄt-255Tとしてディスプレイバッファへ直接に送り出されるＹt-255Tを生成する。更に、計算された基線推定Ｂは最も近い時点の２５５個のサンプル（ＤtからＤt-254T）を表示するときにも使用される。これは、表示装置１３５の右側の表示領域の絶えず調整されているＤＣレベル部分（すなわち、第１の部分）に含まれるデータであり、約１秒の表示距離をカバーしている。
【００２０】
　次に図５を参照すると、本発明の方法が新たなＥＣＧ波形データを受信し、新たなデータを最も近い時点のデータサンプルと共に表示するプロセスを表す流れ図３００が示されている。ブロック３０２から始まり、フィルタ補正前データバッファ２０６に先に格納されていた５１１個のフィルタ補正前データサンプルをシフトする（Ｘt-510Tにある以前の最も古いサンプルは排除される）。これで、ブロック３０４に示すように、新たな未補正ＥＣＧデータサンプルがフィルタ補正前データバッファ２０６のＸtに受け入れられるための経路がクリアになる。
【００２１】
　最も新しいＥＣＧサンプルがＸtに受け入れられると、ブロック３０６に示すように、５１１個の現在フィルタ補正前サンプルから新たな基線推定が計算される。これは、コンボルーション式Ｂ＝ＦＩＲlowpass＊Ｘt-nT（ｎ＝０から５１０）により表現されても良い。次に、ブロック３０８で、フィルタ補正データバッファ２０２の先のデータサンプル（合わせて１１８４）をシフトオーバして、Ｙt-255Tと指示されるサンプルデータを保持する記憶要素２０４からシフトアウトされたデータサンプルを受け入れるための余地をあける。これにより、記憶要素２０４はクリアされて、新たに計算された基線推定Ｂに基づいて最前の補正サンプルを受け入れる。このとき、先に述べた通り、Ｙt-255T＝Ｘt-255T−Ｂである。このステップはブロック３１０に示されている。同様に先に述べた通り、方法はブロック３１２へ進み、最新の２５５個のフィルタ補正前データサンプル（すなわち、ＸtからＸt-254T）からＢを減算することにより、それらのデータサンプルから追加の「一時」表示サンプルを作成する。最後に、ブロック３１４では、ディスプレイバッファ２００にロードされた更新済みの１１８５個の補正サンプルと、２５５個の一時表示サンプルを使用して、表示領域全体を再生する。ディスプレイバッファ２００にロードされたデータサンプルＤt-nTの供給源は次のように記述される。
【００２２】
　        Ｄt-nT＝Ｘt-nT−Ｂ（ｎが０から２５４の場合）
　        Ｄt-nT＝Ｙt-nT（ｎが２５５から１４３９の場合）
　以上説明した方法の効果は、波形が表示に沿ってスクロールされている間に波形の動いている動画を観測することにより最も良く理解される。図６（ａ）及び図６（ｂ）は、図３（ａ）〜図３（ｆ）で使用した三角形インパルス／方形波妨害例にデータ表示方法を適用した場合の成果をある程度実例として説明している。図からわかるであろうが、基線補正は実際には生信号における方形波の立ち上がり端の前に始まっており、従って、最初の基線の上方へのシフトはフィルタ補正前波形の上方へのシフトのほぼ半分であるにすぎない。これは、表示の補正ポイントの前後に起こるフィルタ補正前データサンプルを使用する対称ＦＩＲフィルタを介する遅延基線補正の結果である。また、基線に戻るのに要する時間は、信号の歪みを伴わないだけで、強力な０．５９７Ｈｚフィルタ（図３（ｃ）及び図３（ｄ））と比較して好都合である。
【００２３】
　三角形インパルス波形に加えられる低周波ランダム雑音の例を更に現実に即した形で図７（ａ）〜図７（ｃ）に示す。図７（ａ）は、フィルタ補正前生信号に対する０．０３７Ｈｚフィルタの性能を示す。図示されるように、０．０３７Ｈｚフィルタは基線補正をごくわずかしか改善しておらず、フィルタ補正信号からのインパルスはフィルタ補正前生信号より基線にごくわずかに近づいているだけである。図７（ｂ）では、０．５９７Ｈｚフィルタはランダム雑音をはるかに良く基線補正しているが、この場合にも、三角形ＱＳＲインパルス自体の歪みは大きくなっている。同様に、図７（ｃ）に示す放物線補正技法は基線補正と信号の歪みを相殺させてはいるが、依然として希望するほど強力ではない。
【００２４】
　これに対し、ランダム雑音に対して本発明の方法を適用した成果を図８（ａ）から図８（ｈ）に示すが、これらの図は、波形が時間の経過に伴ってスクロールされる間にコンピュータの表示画面から取り出した順次画像である。このデモンストレーションプログラムにおいては、コンピュータ表示画面は第１の部分４０４と、破線４０６により第１の部分４０４から分離された第２の部分４０２とに分割されている。尚、コンピュータ表示画面の第１の部分と第２の部分との間に示される実際の時間分割は、その性質上、単なる例であり、先に説明した特定の１４４０サンプル表示には必ずしも対応しないことを理解すべきである。むしろ、図８（ａ）から図８（ｈ）の意義は第１の部分と第２の部分とで表示技法が異なるというところにある。
【００２５】
　図８（ａ）は、所持時間ｔ0に撮影された第１の画像であり、表示の右側に第１の生信号三角形インパルスδ（０）が現れている。表示画面の第１の部分４０４には、これに対応する、δ'（０）とラベル付けされたフィルタ補正インパルスが示されている。表示画面の第２の部分４０２は従来のように波形データをスクロールし、第１の部分４０４は連続ＤＣレベル補正を行っている。
【００２６】
　図８（ｂ）は、時間ｔ1で取り出された画面ショットである。図示されるように、基線に関するδ'（０）の位置は時間ｔ0（想像線で示されている）から下方へ移動することにより変化している。これは、表示の第１の部分４０４の中で能動的基線補正が行われていることを反映している。特に、破線４０６を中心として、対称ＦＩＲフィルタ補正により特定基線補正が判定される。しかし、表示の第２の部分４０２では、単純にスクローリングしているだけであるので、時間に関するフィルタ補正信号の振幅変化は見られない。時間ｔ2で取り出された次の順次ショットを図８（ｃ）に示す。（ここで、図の間の特定の時間間隔は必ずしも等間隔の期間として設定されておらず、図８の例は単に実例を示すという性質を持っているにすぎないことに注意すべきである。）この場合にも、δ'（０）の位置はｔ1及びｔ0における位置と比較して基線に関してわずかに下がり続けている。表示をｔ0とｔ2の間でリアルタイムで観測したならば、インパルスδ'（０）は右から左へ移動する間に降下するように見えるであろう。同時に、フィルタ補正波形が右から左へ移動する間に、第２の部分４０２に位置するフィルタ補正波形の部分には振幅変化は全く存在しないであろう。
【００２７】
　次に図８（ｄ）を参照すると、時間ｔ3で、表示の第１の部分４０４の対応するフィルタ補正インパルスδ'（１）と共に、生信号波形に新たなインパルスδ（１）が現れている。この時間までに、δ'（０）は第２の部分４０２までスクロールオーバし終わっており、表示に沿ったスクロール時間の残り部分の間にそれ以上の振幅補正を受けることはない。一方、時間ｔ4で取り出された図８（ｅ）においては、インパルスδ'（１）はｔ3とｔ4の間の基線に関して上昇している。この時間中、δ'（０）は表示の第２の部分４０２に位置しているため、その振幅位置を変化させなかったことがわかるであろう。
【００２８】
　次に時間ｔ5の図８（ｆ）に目を向けると、δ'（１）は第１の部分４０４から第２の部分４０２へほぼ移行しようとする時点で、ｔ4のときの位置から大きく降下したことが見て取れる。加えて、表示の最も右側の部分に新たなインパルスδ'（０）（及びδ'（２））が現れている。更に、δ'（０）の振幅位置は更に左へ移動する間に同じままであったことがわかるであろう。次に、ｔ6で取り出された図８（ｇ）には、第１の部分４０４にとどまっている間のδ'（２）の下方基線補正移動が示されている。この場合にも、第２の部分４０２にあるδ'（１）及びδ'（０）は共に基線に関してシフトしていない。最後に、図８（ｈ）は時間ｔ7で取り出された画面ショットであり、δ'（２）は第２の部分４０２へスクロールオーバする前に上方へ引き戻されたことがわかる。この時点で、第１の部分４０４に新たな波形δ'（３）が現れており、δ'（０）は表示の左側から完全に消え去っている。このように、図８（ａ）から図８（ｈ）を検討すると、表示の第１の部分４０４における「ホイップ状」補正動作を理解する手段が得られる。
【００２９】
　最後に、図９（ａ）から図９（ｅ）に示す実際のＥＣＧ波形に本発明の方法を適用した結果と、先に説明した従来のフィルタ補正技法を適用した結果を比較する。図９（ａ）には、３つの個別のＥＣＧ生信号の読みが示されており、ランダム雑音は実際には下の２つの波形を互いに交差させている。図９（ｂ）では、下の２つの波形は依然として交差しているので、同じランダム雑音は０．０３７Ｈｚ高域フィルタによって有効に補正されてはいない。図９（ｃ）では、０．５９７Ｈｚ高域フィルタは下の２つの波形のもつれをほぐすことにより強力な基線補正を実行しているが、この場合にも、信号の歪みによりＥＣＧの解釈ミスを引き起こす可能性がある。図９（ｄ）では、放物線基線補正は０．０３７Ｈｚフィルタと同様に相対的に有効ではない。しかし、図９（ｅ）は、本発明の方法がＥＣＧを歪ませずに基線ドリフトを除去するのに最も有効であることを示している。
【００３０】
　以上説明した方法の適用範囲は全てのリアルタイム心臓モニタ、非侵入電気生理学（ＥＰ）システム、運動負荷試験装置、除細動器、並びにリアルタイムリズムモード及びＣＲＴ又はＬＣＤ表示装置を伴うＥＣＧチャートを含む。しかし、データフィルタ補正及び表示技法はＥＣＧデータの表示に限定されず、一般に、例えば、血量計データ、血圧データ又は地質／地震データなどの波形データをリアルタイムで表示することが望まれるどのようなシステムにも適用できるであろう。
【００３１】
　更に、本発明の技法はスクローリング波形表示のみに限定されず、「ワイパーバー」表示と組み合わせても使用できる。ワイパーバー表示においては、通常、幅の狭い垂直のバーが表示画面を横切るように左から右へ移動し、最も近い時点の波形データはバーの左側にすぐ隣接してバーの動きに追随する。バーが画面の右端に達すると、バーは表示の左側に戻るように回り込み、その後、バーの右側の古いデータを消去する。この種の表示では、ワイパーバーは、（例えば）データが絶えずレベル調整されて、部分的に基線を補正する１秒トレーリングウィンドウを有することができるであろう。画面に残っているデータは基線補正され、最終的に消去される前にそれ以上はレベル調整されない。
【００３２】
　ここで開示した発明をコンピュータ又はコントローラで実現されるプロセス並びにそれらのプロセスを実施するための装置の形態で具現化できることは理解されるであろう。また、本発明は、フレキシブルディスク、ＣＤ−ＲＯＭ、ハードドライブ、又はその他の何らかのコンピュータ読み取り可能な記憶媒体などの実在する媒体において具現化される命令を含むコンピュータプログラムコードの形態でも実現可能であり、そのようなコンピュータプログラムコードがコンピュータ又はコントローラにロードされ、実行されるとき、コンピュータは本発明を実施するための装置となる。本発明は、例えば、記憶媒体に格納され、コンピュータ又はコントローラにロードされ且つ／又はコンピュータ又はコントローラにより実行されるにせよ、あるいは電気的配線又はケーブルを介する、光ファイバを介する、又は電磁放射線を介するなどの何らかの送信媒体を介して送信されるにせよ、コンピュータプログラムコード又は信号の形態で具現化されても良く、その場合、コンピュータプログラムコードがコンピュータにロードされ、実行されるとき、コンピュータは本発明を実施するための装置となる。汎用マイクロプロセッサで実現される場合、コンピュータプログラムコードセグメントが特定の論理回路を作成するためにマイクロプロセッサを構成する。
【００３３】
　本発明を好ましい一実施例を参照して説明したが、本発明の範囲から逸脱せずに様々な変更を実施できること及び実施例の要素を等価の要素と置き換えても差し支えないことは当業者には理解されるであろう。更に、本発明の本質的な範囲から逸脱せずに特定の状況又は材料を本発明の教示に適合させるために数多くの変形を実施できるであろう。従って、本発明は本発明を実施するために最良であると考えられる態様として開示された特定の実施例に限定されるのではなく、本発明は添付の特許請求の範囲の範囲内に入るあらゆる実施例を含む。
【図面の簡単な説明】
【００３４】
【図１】ＥＣＧ波形を表す概略図。
【図２】本開示の一実施例を実施するのに適するＥＣＧシステムの一例の概略ブロック線図。
【図３（ａ）】方形波雑音妨害に適用される従来のデータフィルタ補正技法を示す波形表示。
【図３（ｂ）】方形波雑音妨害に適用される従来のデータフィルタ補正技法を示す波形表示。
【図３（ｃ）】方形波雑音妨害に適用される従来のデータフィルタ補正技法を示す波形表示。
【図３（ｄ）】方形波雑音妨害に適用される従来のデータフィルタ補正技法を示す波形表示。
【図３（ｅ）】方形波雑音妨害に適用される従来のデータフィルタ補正技法を示す波形表示。
【図３（ｆ）】方形波雑音妨害に適用される従来のデータフィルタ補正技法を示す波形表示。
【図４】本発明の一実施例に従ったデータフィルタ補正及び表示技法の一例を示す概略図。
【図５】図４に示すデータフィルタ補正及び表示技法を更に示す流れ図。
【図６（ａ）】方形波雑音妨害に本発明によるデータフィルタ補正及び表示技法を適用した結果を示す波形表示。
【図６（ｂ）】方形波雑音妨害に本発明によるデータフィルタ補正及び表示技法を適用した結果を示す波形表示。
【図７（ａ）】ランダム低周波雑音妨害に適用される従来のデータフィルタ補正技法を示す波形表示。
【図７（ｂ）】ランダム低周波雑音妨害に適用される従来のデータフィルタ補正技法を示す波形表示。
【図７（ｃ）】ランダム低周波雑音妨害に適用される従来のデータフィルタ補正技法を示す波形表示。
【図８（ａ）】ランダム低周波雑音妨害に本発明によるデータフィルタ補正及び表示技法を適用した結果を示す波形表示。
【図８（ｂ）】ランダム低周波雑音妨害に本発明によるデータフィルタ補正及び表示技法を適用した結果を示す波形表示。
【図８（ｃ）】ランダム低周波雑音妨害に本発明によるデータフィルタ補正及び表示技法を適用した結果を示す波形表示。
【図８（ｄ）】ランダム低周波雑音妨害に本発明によるデータフィルタ補正及び表示技法を適用した結果を示す波形表示。
【図８（ｅ）】ランダム低周波雑音妨害に本発明によるデータフィルタ補正及び表示技法を適用した結果を示す波形表示。
【図８（ｆ）】ランダム低周波雑音妨害に本発明によるデータフィルタ補正及び表示技法を適用した結果を示す波形表示。
【図８（ｇ）】ランダム低周波雑音妨害に本発明によるデータフィルタ補正及び表示技法を適用した結果を示す波形表示。
【図８（ｈ）】ランダム低周波雑音妨害に本発明によるデータフィルタ補正及び表示技法を適用した結果を示す波形表示。
【図９（ａ）】同じ実ＥＣＧデータの反復例に本発明によるデータフィルタ補正及び表示技法を適用した結果と、従来の技法を適用した結果とを比較した波形表示。
【図９（ｂ）】同じ実ＥＣＧデータの反復例に本発明によるデータフィルタ補正及び表示技法を適用した結果と、従来の技法を適用した結果とを比較した波形表示。
【図９（ｃ）】同じ実ＥＣＧデータの反復例に本発明によるデータフィルタ補正及び表示技法を適用した結果と、従来の技法を適用した結果とを比較した波形表示。
【図９（ｄ）】同じ実ＥＣＧデータの反復例に本発明によるデータフィルタ補正及び表示技法を適用した結果と、従来の技法を適用した結果とを比較した波形表示。
【図９（ｅ）】同じ実ＥＣＧデータの反復例に本発明によるデータフィルタ補正及び表示技法を適用した結果と、従来の技法を適用した結果とを比較した波形表示。
【符号の説明】
【００３５】
　１００…ＥＣＧシステム、１０５…電極セット、１１０…電極、１１５…リード、１２０…システムコントローラ、１２５…信号調整回路、１３０…プロセッサ、１３５…表示装置、１４０…出力装置、１４５…ランニングマシン、１５０…ケーブル、１５５…記録装置、２００…ディスプレイバッファ、２０２…フィルタ補正データバッファ、２０４…記憶要素、２０６…フィルタ補正前データバッファ、４０２…（表示の）第２の部分、４０４…第１の部分
【図１】
2004065981.tif 000003
【図２】
2004065981.tif 000004
【図３（ａ）】
2004065981.tif 000005
【図３（ｂ）】
2004065981.tif 000006
【図３（ｃ）】
2004065981.tif 000007
【図３（ｄ）】
2004065981.tif 000008
【図３（ｅ）】
2004065981.tif 000009
【図３（ｆ）】
2004065981.tif 000010
【図４】
2004065981.tif 000011
【図５】
2004065981.tif 000012
【図６（ａ）】
2004065981.tif 000013
【図６（ｂ）】
2004065981.tif 000014
【図７（ａ）】
2004065981.tif 000015
【図７（ｂ）】
2004065981.tif 000016
【図７（ｃ）】
2004065981.tif 000017
【図８（ａ）】
2004065981.tif 000018
【図８（ｂ）】
2004065981.tif 000019
【図８（ｃ）】
2004065981.tif 000020
【図８（ｄ）】
2004065981.tif 000021
【図８（ｅ）】
2004065981.tif 000022
【図８（ｆ）】
2004065981.tif 000023
【図８（ｇ）】
2004065981.tif 000024
【図８（ｈ）】
2004065981.tif 000025
【図９（ａ）】
2004065981.tif 000026
【図９（ｂ）】
2004065981.tif 000027
【図９（ｃ）】
2004065981.tif 000028
【図９（ｄ）】
2004065981.tif 000029
【図９（ｅ）】
2004065981.tif 000030