JP2001517991A - データ信号適応平均化方法及び装置 - Google Patents

データ信号適応平均化方法及び装置

Info

Publication number
JP2001517991A
JP2001517991A JP54531898A JP54531898A JP2001517991A JP 2001517991 A JP2001517991 A JP 2001517991A JP 54531898 A JP54531898 A JP 54531898A JP 54531898 A JP54531898 A JP 54531898A JP 2001517991 A JP2001517991 A JP 2001517991A
Authority
JP
Japan
Prior art keywords
measurements
way
signal
measurement
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.)
Pending
Application number
JP54531898A
Other languages
English (en)
Inventor
ベイカー、クラーク・アール、ジュニア
ヨーケイ、トーマス・ジェイ
Original Assignee
ネルコー・ピューリタン・ベネット・インコーポレイテッド
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 ネルコー・ピューリタン・ベネット・インコーポレイテッド filed Critical ネルコー・ピューリタン・ベネット・インコーポレイテッド
Publication of JP2001517991A publication Critical patent/JP2001517991A/ja
Pending legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • A61B5/14551Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Molecular Biology (AREA)
  • Optics & Photonics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)

Abstract

(57)【要約】 生理的パラメータを計測するためのシステムでのノイズによる影響を低減する方法。この方法は、(a)生物組織を通じて伝達される電磁エネルギーの少なくとも一つの波長から導かれる複数の計測値を生成する工程、(b)選択計測値と、少なくとも一つの期待計測値特性とを比較する工程、(c)複数の可変重みのうちの一つを、上記の比較工程に基づいた各選択計測値に割り当てる工程であって、これにより、複数の異なった重み付き計測値を生成する、ところの工程、及び(d)生理的パラメータの予測に使用するためのフィルター計測値を得るために、複数の異なった重み付き計測値を平均化する工程、を含む。

Description

【発明の詳細な説明】 データ信号適応平均化方法及び装置 本発明は、生理的パラメータを計測するための方法及び装置に関し、特に、生 理的パラメータを計測するためのシステムにおいてノイズによる影響を低減する ための方法及び装置に関する。特に、データ信号を適応的に平均化するための方 法及び装置に関する。本発明は、パルス酸素計にフィルター化技術を用いて、動 脈血のヘモグロビンの酸素飽和度を予測する。 パルス酸素計は、典型的に、動脈血のヘモグロビンの酸素飽和度に限らず、様 々な血流特性を計測、表示する。酸素計は、光を指や耳といった血液潅流組織を 通過させ、組織での光吸収を光学的に感知する。次に、光吸収量は、計測される べき血液構成成分(例えば、オキシヘモグロビン)量を計算するために使用され る。 組織を通過する光には、一つ又はそれ以上の波長のものが選択され、この波長 の光は、血液に存在する血液構成成分量を表す量だけ血液に吸収される。組織を 通過する光の量は、組織の血液構成成分及び関係する光吸収の変化量に従って変 化する。次に、飽和度は、Beer-Lambertの法則に基づいて計算される。伝統的に 、飽和度は、 二つの波長(例えば、赤色光及び赤外光)での光吸収度を計測することによって 決定される。飽和度は、次に、米国特許第4802486号、米国特許第491 1167号、米国特許第4928692号、米国特許第4969254号、米国 特許第5078136号及び米国特許第5485847号に開示されるように、 “割合比(ratio of ratios)”の値を求めることによって計算できる。 組織を通過する光学信号は、ノイズ及びモーションアーチファクト(motion a rtifact)により劣化する。ノイズ源の一つは、光検出器に到達する周囲の光で ある。他のノイズ源は、他の電子機器からの電磁的な接続である。患者のモーシ ョンもまた、ノイズを発生し、信号に影響を与える。例えば、検出器と皮膚との 間の接続やエミッタと皮膚との間の接続が、いずれか一方が皮膚から離れるよう に動くと、一時的に中断される。また、血液が流体であることから、周囲の組織 と異なって慣性効果に応答して酸素計のプローブが取り付けられる場所で瞬間的 にその体積を変化させることになる。 モーションアーチファクトは、医者が気づくことなく、医者の頼るところのパ ルス酸素計信号を劣化させる。これは、特に、患者のモニターが遠隔にある場合 やモーションが小さすぎる場合、又は医者がセンサーの取付位置ではなく患者の 機器や他の部分をみている場合に起こる。 モーション検出のために加速度計を使用する酸素計システムが米国特許第50 25791号に開示される。モーションが検出されると、モーションによる読取 り値が、誤りとして除去されるか又は指示される。典型的な酸素計では、脈拍信 号の最上点と最下点でとられた計測値を使用して所望の特性を計算する。モーシ ョンにより最上点に誤りが生じ、この結果、計測値が、誤った値で、誤った時間 に記録されるものとなる。 EKG信号をモニターし、酸素計の読取り値と相関して、酸素計の読取り値上 のノイズ及びモーションによる影響を制限するように同期させる他のシステムが 米国特許第4802486号に開示される。これは、周期的なモーション信号に ロックする酸素計の機会を低減させる。 米国特許第5078136号に開示されるシステムは、酸素計信号を分析する 線形挿間及び変化量技術を使用して、ノイズ及びモーションアーチファクトによ る影響を制限する。 本発明は、様々な重みが異なった計測値に割り当てられ、フィルター計測値を 得るために、重み付き計測値が平均化される、ところの生理的パラメータ(例え ば、血中酸素飽和度)を計測するためのシステムにおいてノイズによる影響を低 下させるための技術を提供する。 従って、本発明の第一の態様は、生理的パラメータを計測するためのシステム においてノイズによる影響を低下させるための方法であって、 (a)生物組織を通じて伝達される電磁エネルギーの少なくとも一つの波長から 導かれる複数の計測値を生成する工程、 (b)選択計測値と、少なくとも一つの期待計測値特性とを比較する工程、 (c)複数の可変重みのうちの一つを、上記の比較工程に基づいた各選択計測値 に割り当てる工程であって、これにより、複数の異なった重み付き計測値を生成 する、ところの工程、及び (d)生理的パラメータの予測に使用するためのフィルター計測値を得るために 、複数の異なった重み付き計測値を平均化する工程、を含む。 他の態様では、本発明は、生理的パラメータを計測するためのシステムにおい てノイズによる影響を低減させるための装置を提供する。この装置は、 (a)生物組織を通じて伝達される電磁エネルギーの少なくとも一つの波長から 導かれる複数の計測値を生成するための手段、 (b)選択計測値と、少なくとも一つの期待計測値特性とを比較するための手段 、 (c)複数の可変重みのうちの一つを、上記の比較工程に基づいた各選択計測値 に割り当てるための手段であって、これにより、複数の異なった重み付き計測値 を生成する、ところの手段、及び (d)生理的パラメータの予測に使用するためのフィルター計測値を得るために 、複数の異なった重み付き計測値を平均化するための手段、 を含む。 好適に、この装置は、例えば、米国特許第4621643号、米国特許第47 00708号及び米国特許第4770179号に開示されるように、電磁エネル ギーの少なくとも一つの波長を表す信号に応答して生理的パラメータを計測する ため、システムを校正するための手段を含む。 好適に、この装置は、センサーを含み、センサーは、電磁エネルギーの上記の 少なくとも一つの波長のためのエミッタ、電磁エネルギーを感知し且つこれを表 す第一の信号を生成するための感知手段、センサーを酸素計に検出可能に連結し 且つセンサーと酸素計との間で信号を通信するための手段、及び電磁エネルギー の少なくとも一つの波長を表す第二の信号を与えるための手段、を有する。 本発明は、生理的パラメータが時間とともにどのように変化し、これらパラメ ータがノイズ環境で計測値にどのように関係するかを 説明するために、数学的モデルを使用するフィルター化技術を使用する。このよ うなフィルター化は、生理的パラメータの最適な予測値に到達するように、平均 化重み及び平均化時間の組を変更することができる。 本発明の技術は、動脈血のヘモグロビンの酸素飽和度を決定するために、パル ス酸素計に関連して使用できる。少なくとも部分的に、バンド外ノイズを除くた めに、バンドパスフィルターが使用され、0.5Hz以下、10Hz以上のデー タを減衰することができる。フィルター化データは、次に、後述する飽和度計算 アルゴリズムを使用して処理される。 本発明は、適応フィルター化技術(例えば、カルマン(Kalman)フィルター法 )を利用して、血液酸素飽和度を計算する。カルマンフィルター法は、パラメー タが時間とともに変化するとき、最小二乗的にパラメータを一致させることがで きる。伝統的な、予測量を平均化又はローパスフィルター化する古典的最小二乗 法(CLS)を利用できる。カルマンフィルター法は、実質的に同一の結果を達 成するが、カルマンフィルターは、平均化の最適量を計算する。適当なカルマン フィルターアルゴリズムは、「Introduction to Random Signals and Applield Kalman Filtering」(RG Brown及びPYCHwang著、第2版、1992年、John Willey & Sons)に開示される。カルマン フィルター定理を利用したカルマン心臓ゲート化平均化プロセッサが酸素飽和度 の計算に使用され、このプロセッサは、例えば、心拍数を計算するためのアルゴ リズムにより供給され、又はECG波形からの心拍数に基づいたトリガーにより ゲート化される。酸素計から心拍数を計算するための技術の詳細は「Method and Apparatus for measuring Pulse Rate and Saturation」(国際出願番号IB9 6/ )(代理人用整理番号P21988B)に開示される。これに開示さ れる情報は、本出願に参照文献として組み入れられる。 カルマンフィルターを使用すると、本発明は、適当な前処理データを、適当な 減光係数を使用してオキシヘモグロビン及び全ヘモグロビン濃度に対応する量に 変換することを含む。これら二つの変換量のある瞬間の比は、飽和度を与える。 ある瞬間の飽和度値は、減光係数を使用するか又は割合比から直接的に計算され 得る。この方法は、(所望の場合、最大値又は最小値は使用され、カルマンフィ ルター化されるが、)パルス探索アルゴリズムのような最大値又は最小値を探索 する必要がない。最大値/最小値の比(つまり、事象ベースのアルゴリズム)で はなく、ある瞬間の比(つまり、時間ベースのアルゴリズム)を使用することは 、それが到達すると事象を払いのけデータを修正するコードを維持する。したが って、本発明の方法は、パルス探索、事象ベースの飽和度計算アルゴリズムよ りも簡単に適用できる点に利点がある。 好適に、平均化した異なった重み付き計測値の数は、上記の割り当て工程に応 答して変化する。 好適に、異なった重み付き計測値の第一の数は、フィルター計測値を得るため に平均化され、ここで、第一の数は、割り当て工程での複数の継続する異なった 重み付き計測値に割り当てられる重みに従って変化する。次に、時間外(又は、 オーバータイム(over time))に、各波長のための複数のフィルター計測値が 生成されることが好適である。 好適に、複数の可変重みは、第一及び第二の重みに対応する計測値の第一及び 第二の組を含む。比較工程及び割り当て工程は、次に、計測値の第一の組で行わ れる。本方法が、計測値の第一の組に割り当てられるものと同一の計測値の第二 の組の重みに割り当てる工程を含み、これにより、第一及び第二の波長の各々の ための別々のフィルター計測値が得られることが好適である。 好適に、複数の計測値は、少なくとも二つの波長からのデータを組み合わせる ことによって得られる。 好適に、各選択計測値へ割り当てられた可変重みは、部分的に、選択計測値と 直前の計測値との間の変化量に基づく。 好適に、本方法の上記の生成する工程は、 (a)電磁エネルギーの少なくとも一つの波長を表す信号の対数をとる工程であ って、これにより、第一の信号を生成する、ところの工程、及び (b)第一の信号をバンドパスフィルター化する工程であって、これにより、第 二の信号を生成し、複数の計測値が第二の信号から導かれる、ところの工程、 を含む。 好適に、上記の生成する工程は、次に、第三の信号を生成するために第二の信 号を正規化する工程を含み、複数の計測値が第三の信号から導かれる。上記の生 成する工程は、次に、第四の信号を生成するために第三の信号の導関数をとる工 程を含むことができ、複数の計測値が第四の信号から導かれる。 好適に、各選択計測値と対応する直前の計測値との間の相似点は、期待値と比 較される。期待値は、生理的モデルに対応し得る。それは、複数の対応する直前 の計測値の間の変化量に対応し得る。 好適に、可変重みを割り当てるため、選択計測値と比較される対応する直前の 期待値は、フィルター計測値に対応する。 好適に、少なくとも一つの期待計測値特性は、少なくとも一つの直前のフィル ター計測値に対応する予測値を含む。 好適に、選択計測値は、デジタル信号を含む。 好適に、本方法の工程(a)で生成した複数の計測値は、一連の心拍に対応す る。 好適に、生成した重み付き計測値の可変重みは、複数の異なったゼロでない数 を含む。 好適に、本方法の工程(a)で生成した複数の計測値は、時間ベースであり、 被駆動事象(又は、イベント・ドライブン(event driven))ではない。各フィ ルター時間ベース計測値は、次に、下記のうちの少なくとも一つに対応する。 1.動脈血中のヘモグロビンの酸素飽和度を表す割合比の現在の値。 2.動脈血中のヘモグロビンの酸素飽和度の現在の値。少なくとも一つの期待計 測値特性は、次に、酸素飽和度値の変化量を表す一定値を含む。 複数の時間ベース計測値は、比較工程前に正規化され、これにより、時間ベー ス計測値上のモーションアーチファクトに対応するノイズによる影響を低減する 。 好適に、本発明の方法は、電磁エネルギーの少なくとも一つの波長を表す信号 を与える工程を含む。 好適に、本発明の方法は、電磁エネルギーの波長に対応するデータをフィルタ ー化する工程を含み、患者の心拍数の整数倍でない モーション及びノイズエネルギーが減衰される。 本発明は、生理的パラメータを計測するときに、ノイズによる影響を低減する ことを含む。それは、ノイズによる影響を低減するための装置を含み、この装置 は、 生物組織を通じて伝達される電磁エネルギーの少なくとも一つの波長から導か れる複数の計測値を生成するための手段、 電磁エネルギーの少なくとも一つの波長を表す信号を与えるための手段、 選択計測値と、少なくとも一つの期待計測値特性とを比較するための手段、 複数の可変重みのうちの一つを、上記の比較工程に基づいた各選択計測値に割 り当てるための手段であって、これにより、各波長のための複数の異なった重み 付き計測値を生成し、可変重みは、各選択計測値と対応する直前の計測値との間 の相似点に応答して部分的に割り当てられ、可変重みは、複数の異なったゼロで ない数を含む、ところの手段、 生理的パラメータの予測に使用するためのフィルター計測値を得るために、複 数の異なった重み付き計測値を平均化するための手段、及び 電磁エネルギーの少なくとも一つの波長を表す信号に応答して、 生理的パラメータを計測するために、システムを校正するための手段、 を含む。 本発明は、また、生理的パラメータをモニターするためのモニターを含み、モ ニターはセンサーと一緒に使用される。センサーは、電磁エネルギーの少なくと も一つの波長を放出するための手段、電磁エネルギーを感知し且つそれを表す第 一の信号を生成するための感知手段、センサーを酸素計に検出可能に連結し且つ センサーと酸素計との間の信号の通信を与えるための手段、及び電磁エネルギー の少なくとも一つの波長を示す第二の信号を与えるための手段、を有する。モニ ターは、 第一の信号から導かれた複数の計測値を生成するための手段、 選択計測値と、少なくとも一つの期待計測値特性とを比較するための手段、 複数の可変重みのうちの一つを、上記の比較工程に基づいた各選択計測値へ割 り当てるための手段であって、これにより、複数の異なった重み付き計測値を生 成し、可変重みは、各選択計測値と対応する直前の計測値との間の相似点に応答 して、部分的に、割り当てられ、可変重みは、複数の異なったゼロでない数を含 む、ところの手段、 生理的パラメータを予測することに使用するためのフィルター計測値を得るた めに、複数の異なった重み付き計測値を平均化するための手段、及び 第二の信号に応答して生理的パラメータを計測するために、モニターを校正す るための手段、 を含む。 次に、添付の図面を参照して、本発明の例のみを説明する。 図1は、患者のヘモグロビンの酸素飽和度のような生理的パラメータを計測す るための装置の略図である。 図2は、図1に示す装置でのデータフローを示すブロック図である。 図3は、古典的最小二乗アルゴリズムの性能とカルマンアルゴリズムの性能を 比較するグラフである。 図4は、カルマン心臓ゲート化平均化フィルターの入出力を比較するグラフで ある。 図面を参照すると、図1は、患者のヘモグロビンの酸素飽和度のような生理的 パラメータを計測するための装置をしめす。センサー/酸素計組合せ60は、セ ンサー61、及び酸素計モニター62、から構成される。センサー61は、それ ぞれ、赤外線及び赤色スペクトル範囲の波長放出特性を典型的に有するLED6 3、64を含 む。光ダイオードセンサー65が、LED63、64によって伝えられた光を受 信する。抵抗器66(又は同様の電気的インピーダンス関連)が、特定波長、又 は波長に対する対照インピーダンス値により特定されるような波長の組合せに一 致するように選択される。デコード手段67が、抵抗器66のインピーダンス値 を決定し、特定のセンサー61の送信特性に対応する適当な減光係数が生成され る。よって、酸素計は、再校正せずに様々な波長の光を放出するLEDを有する 様々なセンサーと一緒に使用できる。センサー61は、コネクターの手段によっ て酸素計モニター62に検出可能に連結される。このようなセンサー/酸素計組 合せの一例は、米国特許第4621643号に開示される。 センサーから受信したデータは、図2に示すスキームに従って処理される。そ れは、米国特許第5348004号に開示されるタイプの装置を使用して処理で きる。最初の処理工程1、2では、(通常、赤色及び赤外LEDからの)データ の自然対数がとられ、データは、バンドパスフィルター化される(工程1)。次 に、フィルターデータは、酸素飽和度の計算用のアルゴリズムによって処理され る。フィルターデータを処理するためのアルゴリズムは、カルマンフィルター化 を使用でき(工程11)、これは、心臓ゲート化平均化しても、しなくてもよい (工程9)。これらフィルター化技術 は、上記の刊行物「Introduction to Random Signals and Applied Kalman Filt ering」に開示される。 カルマンフィルター化を使用すると、予測されるべきパラメータx(例えば、 酸素飽和度又は心拍数)が、幾らか周期的に、時間的に変化する。x値が、ある 時間の幾らかのサンプルで知られている場合、xは、直前の値からの変化が少し だけ又は全くないことが期待され得る。Qは、この差の分散である。パラメータ xは、直接的に計測されない。実際に計測されるのは、パラメータzであり、こ れは、xに一定値Hを掛け、これに計測ノイズを加えたものである。Rは、この 計測ノイズの分散である。これらは、下記のように書ける。 xk=xk-1+nk Q k=Hkk+nk R zを知ってx値を予測できる能力及びxの最終的な予測は、R及びQにより量 化された二つのノイズに関連する。カルマンフィルターは、予測誤差Pと呼ばれ るパラメータで、二つのノイズを量化する。カルマンフィルターは、また、カル マン利得Kと呼ばれる中間項を使用する。Po-1は、ゼロ値で初期化される。次 に、各々の新しいデータ点Kでは、以下の工程が行われる。 Pk -1=Pk-1 -1+Hk 2k -1 k=Pkkk -1 k=xk-1+Kk(zk−Hkk-1) Pk+1=Pk+Qk カルマンフィルター(工程11)で、飽和度は変化でき、モデルは、二つの部 分に分けられる。第一の部分は、 vk=ukk+nR k である。 つまり、伝達される前処理したデータの比は、計測ノイズを除いた飽和度値で ある。データの点在は、ノイズ分散の実時間計測値を与える。第二の部分は、平 均飽和度が経時変化しないが、それが変化するならば、その変化の標準偏差は実 測可能の変化量である幾らかの一定値Q1/2であることを示す。つまり、第二の 部分は、 sk=sk-1+nQ k である。 この第二の方程式は、飽和度が2秒間に10ポイントだけ変化した場合にそれ が計測ノイズに起因するものであることを認知する能力をカルマンフィルターに 与える。次に、カルマンフィルターは、直前の値とともに、計算した飽和度をさ らに平均化し、この変化を、生理学で期待されるものにさらに一致させる。対照 的に、この変化が境界内である場合は、カルマンフィルターは、ほんの少しだ け平均化する。 Rの値は、最後のNポイント上のusとvとの間の差から予測される、ここで 、使用者は、値Nを特定する。一つの実施例では、飽和度のためのカルマンモデ ルは、また、パルスのより小さい部分に小さい重みを与え、より大きい部分に大 きい重みを与え、実際の分散に小さい増分値を加えて、計測システムの誤差固有 値(例えば、ハードウェアノイズ)を表す。 第二のカルマンフィルター(工程12)では、カルマンフィルターは、変化を 飽和度の時間微分に制限する。このフィルターの方程式は、飽和度の時間分散が 統計的に経時変化しないことを示す。 dxk/dt=dxk-1/dt+nQ k dzk/dt=dxk/dt+nR k ここで、zは、第一のカルマンフィルターからの飽和度の予測値であり、xは 、変化をその時間微分に制限した後の飽和度の予測値である。この実施例では、 パラメータnQは、実測変化量の第二の微分に対応するように選択されることが 好適である。カルマンフィルターの一般形では、これら二つの別々のフィルター は、一つのフィルターに組み合わせることができる。これらを別々に分けること によって、行列代数を使用する必要が除かれ、各カルマンフィルターを別々に試 験できる。 計測ノイズは、使用されるデータ値の周りのウィンドウを中心に置く(又はセ ンターリングする)ことによって予測される。このセンターリングは、ノイズの より正確な予測値を与えるが、半分のウィンドウ長さたけカルマンフィルターの 出力を遅らせる。1秒のウィンドウが使用されると、フィルターがモーションに 直ぐに応答し、飽和度予測値での1/2秒の遅れは、臨床的に著しいものではな い。心臓ゲート化平均化カルマンフィルター化 カルマンCGA(Kalman cardiac gated averaging)アルゴリズムは、カルマ ン飽和度アルゴリズムと直列的に適用される(工程8、9)。カルマンCGAプ ロセッサは、継続する血量計パルス又は波形を最適に平均化して、最適な血量計 波形を作り出す。下記の第一の方程式は、計測ノイズを除く、計測血量計形状と 平均化血量計波形とを相関する。 zk=xk+nk Q Q値は、データから予測される。以下の方程式に従って、新しいパルスは、許 容できる数パーセント以上、平均化血量計パルスと異なり得ない。 xk=xk-N+nk R カルマン心臓ゲート化平均化モデルは、入力波形が著しく変化し た場合、より多数のデータポイントを自動的に平均化し、波形が、期待される生 理的な変化に基づいた仮説に従う場合、直ちに更新する能力を有する。カルマン ゲート化平均化は、米国特許第4802486号に開示されるような、ネルコー ・インコーポレイテッドから入手可能の酸素計(N200及びN3000(商標 ))に使用される従来技術の飽和度計算方法を著しく改善する。図4は、本発明 の一つの実施例に従ったカルマンフィルターの入出力の例を示す。トリガー波形 100は、ECGのR波又は心拍計算方法からのものである。未修正データ波形 102は、モーションにより損なわれたものであり、可変平均化により、カルマ ン心臓ゲート化平均化技術が、規則的なフィルター化波形104を維持させるこ とができる。予測される残余106は、計測データ上のノイズと時間的によく相 関する。 カルマン心臓ゲート化平均化技術が、CLS飽和度計算技術、カルマン飽和度 計算技術又はその他の技術のうちのいずれかによって処理するための酸素計デー タパルスを形づけるために使用され得る。いずれかの実施例は、心臓ゲート化平 均化トリガーのように酸素計データを処理するアルゴリズムによって生成される 心拍数又はECG心拍数を使用できる。カルマン飽和度計算技術は、カルマン心 臓ゲート化平均化技術なしで使用できる。 図2を参照すると、二つの飽和度値が、バンドパスフィルターからのデータか ら計算される。一つの飽和度値は、ECG波形からのトリガーを使用し、ハーモ ニックフィルター(工程7)及び心臓ゲート化平均化を有するカルマンフィルタ ーを使用して得られる。 ハーモニックフィルター(工程7)は、心拍数の整数倍のエネルギーのみがフィ ルターを通過できるように、IR及び赤色波をデジタル式にフィルター化する。 ハーモニックフィルター(工程7)の応答は、それに印加される心拍数信号とと もに変化し、心拍数にないモーション及びノイズエネルギーを減衰させる。この アレンジメントでは、これに続くカルマンCGAによるフィルター化(工程9) 及び/又は後述する飽和度計算(工程11)アルゴリズムが、幾らかの重み化及 び平均化を、フィルター化データストリームに基づいて、IR及び赤色データス トリームに適用する。 本発明の装置に使用し得る、ノイズによる液強を低減するハーモニックフィル ターの使用の詳細は、本出願と同時に出願した「Method and Apparatus for Har monically Filtering Data」(国際出願番号IB96/ )(代理人用整 理番号P21977D)に開示される。本発明の装置に使用し得る、ノイズによ る影響を低減する適応組合せフィルター化の使用の詳細は、本出願と同時に出願 した「Method and Apparatus for Mcasuring Measuring Pulsc Rate and Saturation」(国際出願番号IB96/ )(代理人用整理番号P219 77B)に開示される。 本発明のこれら特徴に関する情報は、これら出願書類に開示され、参照文献と して、本出願の明細書に組み入れる。 第二の飽和度値は、カルマンフィルター(工程11)を適用することによって 得られる。事象ベースの心臓ゲート化平均化を有する第一のカルマンフィルター とは対照的に、第二のフィルターは時間ベースである。第二のフィルターは、バ ンドパスフィルターからのデータ及び第一のフィルターからのデータで動作する 。第二のフィルターの適用の前に、不可能な飽和度計算(例えば、負の飽和度) になるデータポイントが除去される(工程10)。第二のフィルターの適用後、 最も良好な飽和度値が、各値に関連した信頼レベルに従って選択される(工程1 2)。 第二のフィルター後の飽和度値は、どれをどのように表示するかを決定する適 当な後処理を行った後に表示される(工程14)。酸素飽和度の信頼レベルは、 どの飽和度が信頼できるかを決定する、酸素計データで行われるアルゴリズムか ら入手可能のメトリック(metric)から予測できる。例えば、信頼レベルは、飽 和度レベルを計算する信号のエージ(age)及び予測値からのそのレベルの分散 に従って決定できる。飽和度レベルの信頼性を評価する技術のより詳 細は、本出願と同時に出願した「Method and Apparatus for Arbitrating to ob tain Best Estimates for Blood Constituent Values and Rejecting Harmonics 」(国際出願番号IB96/ )(代理人用整理番号P21977C)に 開示される。本発明のこれら特徴に関する情報は、この出願書類に開示され、参 照文献として本出願の明細書に組み入れる。

Claims (1)

  1. 【特許請求の範囲】 1.生理的パラメータを計測するためのシステムにおいてノイズによる影響を低 減するための方法であって、 (a)生物組織を通じて伝達される電磁エネルギーの少なくとも一つの波長から 導かれる複数の計測値を生成する工程、 (b)選択計測値と、少なくとも一つの期待計測値特性とを比較する工程、 (c)複数の可変重みのうちの一つを、前記の比較する工程に基づいた各選択計 測値に割り当てる工程であって、これにより、複数の異なった重み付き計測値を 生成する、ところの工程、及び (d)生理的パラメータの予測に使用するためのフィルター計測値を得るために 、複数の異なった重み付き計測値を平均化する工程、 を含む方法。 2.請求項1の方法であって、 平均化される異なった重み付き計測値の数が、前記の割り当てる工程に応答し て変化する、ところの方法。 3.請求項1の方法であって、 異なった重み付き計測値の第一の数が、フィルター計測値を得るために平均化 され、 この第一の数は、前記の割り当てる工程で複数の引き続く異なっ た重み付き計測値に割り当てられる重みに従って変化し、時間外に、複数のフィ ルター計測値が、各波長のために生成される、 ところの方法。 4.請求項1の方法であって、 複数の計測値が、それぞれ、第一及び第二の波長に対応する計測値の第一及び 第二の組を含み、 前記の比較する工程及び割り当てる工程が、計測値の第一の組で行われ、 計測値の第一の組に割り当てられる重みと同一の重みを計測値の第二の組に割 り当てる工程をさらに含み、 これにより、第一及び第二の波長の主主のための別々のフィルター計測値を得 る、 ところの方法。 5.請求項1の方法であって、 複数の可変重みが、少なくとも二つの波長からのデータを組み合わせることに よって得られる、 ところの方法。 6.請求項1の方法であって、 各選択計測値に割り当てられる可変重みが、部分的に、選択計測値と直前の計 測値との間の変化量に基づいたものである、 ところの方法。 7.請求項1の方法であって、 前記の生成する工程が、 (a)電磁エネルギーの少なくとも一つの波長を表す信号の対数をとる工程であ って、これにより、第一の信号を生成する、ところの工程、及び (b)この第一の信号をバンドパスフィルター化する工程であって、これにより 、第二の信号を生成し、第二の信号から複数の計測値が導かれる、ところの工程 、 を含む、 ところの方法。 8.請求項7の方法であって、 前記の生成する工程が、 第二の信号を正規化する工程、 をさらに含み、 これにより、第三の信号を生成し、第三の信号から複数の計測値が導かれる、 ところの方法。 9.請求項8の方法であって、 前記の生成する工程が、 第三の信号の導関数をとる工程、 をさらに含み、 これにより、第四の信号を生成し、第四の信号から複数の計測値が生成される 、 ところの方法。 10.請求項1の方法であって、 各選択計測値と、対応する直前の計測値との間の相似点が、期待値と比較され る、 ところの方法。 11.請求項10の方法であって、 期待値は、生理的モデルに対応する、 ところの方法。 12.請求項10の方法であって、 期待値は、複数の対応する直前の計測値の間の変化量に対応する、 ところの方法。 13.請求項1の方法であって、 対応する直前の計測値は、フィルター計測値に対応する、 ところの方法。 14.請求項1の方法であって、 少なくとも一つの期待計測値特性は、少なくとも一つの直前のフィルター計測 値に対応する予測値を含む、 ところの方法。 15.請求項1の方法であって、 選択計測値は、デジタル信号を含む、 ところの方法。 16.請求項1の方法であって、 工程(a)で生成される複数の計測値は、一連の心拍に対応する、 ところの方法。 17.請求項1の方法であって、 前記の生成された重み付き計測値の可変重みは、複数のゼロでない数を含む、 ところの方法。 18.請求項1の方法であって、 工程(a)で生成される複数の計測値は、時間ベースであり、被駆動事象では ない、 ところの方法。 19.請求項18の方法であって、 各フィルター化時間ベース計測値は、動脈血のヘモグロビンの酸 素飽和度を表す割合比の現在の値に対応する、 ところの方法。 20.請求項18の方法であって、 各フィルター化時間ベース計測値は、動脈血のヘモグロビンの酸素飽和度の現 在の値に対応する、 ところの方法。 21.請求項20の方法であって、 少なくとも一つの期待計測値特性は、酸素飽和度値の変化量を表す一定値を含 む、 ところの方法。 22.請求項18の方法であって、 複数の時間ベース計測値は、前記の比較する工程の前に正規化され、これによ り、幾らかの時間ベース計測値上のモーションアーチファクトに対応するノイズ による影響を低減する、 ところの方法。 23.請求項1の方法であって、 少なくとも一つの電磁エネルギーを表す信号を与えることを含む、 ところの方法。 24.請求項1の方法であって、 電磁エネルギーの前記波長に対応するデータをフィルター化する工程を含み、 患者の心拍数の整数倍にないモーション及びノイズエネルギーが、減衰される 、 ところの方法。 25.生理的パラメータを計測するためのシステムにおいてノイズによる影響を 低減するための装置であって、 (a)生物組織を通じて伝達される電磁エネルギーの少なくとも一つの波長から 導かれる複数の計測値を生成するための手段、 (b)選択計測値と、少なくとも一つの期待計測値特性とを比較するための手段 、 (c)複数の可変重みのうちの一つを、上記の比較工程に基づいた各選択計測値 に割り当てるための手段であって、これにより、複数の異なった重み付き計測値 を生成する、ところの手段、及び (d)生理的パラメータの予測に使用するためのフィルター計測値を得るために 、複数の異なった重み付き計測値を平均化するための手段、 を含む装置。 26.請求項25の装置であって、 電磁エネルギーの少なくとも一つの波長を示す信号に応答して生 理的パラメータを計測するために、システムを校正するための手段を含む、 ところの装置。 27.請求項25の装置であって、 センサーを含み、 このセンサーは、 電磁エネルギーの前記の少なくとも一つの波長を放出するためのエミッタ、 電磁エネルギーを感知し且つそれを表す第一の信号を生成するための感知手段 、 センサーを酸素計に検出可能に連結し且つセンサーと酸素計との間の信号の通 信を与えるための手段、及び 電磁エネルギーの少なくとも一つの波長を示す第二の信号を与えるための手段 、 を有する、 ところの装置。
JP54531898A 1997-03-21 1997-03-21 データ信号適応平均化方法及び装置 Pending JP2001517991A (ja)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/IB1997/000287 WO1998042249A1 (en) 1997-03-21 1997-03-21 Method and apparatus for adaptively averaging data signals

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2009116451A Division JP4571220B2 (ja) 2009-05-13 2009-05-13 データ信号適応平均化方法及び装置

Publications (1)

Publication Number Publication Date
JP2001517991A true JP2001517991A (ja) 2001-10-09

Family

ID=11004541

Family Applications (1)

Application Number Title Priority Date Filing Date
JP54531898A Pending JP2001517991A (ja) 1997-03-21 1997-03-21 データ信号適応平均化方法及び装置

Country Status (9)

Country Link
EP (1) EP1006863B1 (ja)
JP (1) JP2001517991A (ja)
AT (1) ATE251864T1 (ja)
AU (1) AU720450B2 (ja)
CA (1) CA2283858C (ja)
DE (1) DE69725621T2 (ja)
ES (1) ES2208876T3 (ja)
HK (1) HK1031189A1 (ja)
WO (1) WO1998042249A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004503315A (ja) * 2000-06-16 2004-02-05 アイシス・イノベーション・リミテッド 異なるセンサからの測定の組み合わせ
WO2011096240A1 (ja) * 2010-02-05 2011-08-11 日本電気株式会社 生体情報計測器、携帯端末装置、生体情報計測方法およびプログラム

Families Citing this family (42)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
MX9702434A (es) 1991-03-07 1998-05-31 Masimo Corp Aparato de procesamiento de señales.
US7376453B1 (en) 1993-10-06 2008-05-20 Masimo Corporation Signal processing apparatus
US6575905B2 (en) 2000-09-22 2003-06-10 Knobbe, Lim & Buckingham Method and apparatus for real-time estimation of physiological parameters
US9282925B2 (en) 2002-02-12 2016-03-15 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US10022078B2 (en) 2004-07-13 2018-07-17 Dexcom, Inc. Analyte sensor
US7497827B2 (en) 2004-07-13 2009-03-03 Dexcom, Inc. Transcutaneous analyte sensor
US8260393B2 (en) 2003-07-25 2012-09-04 Dexcom, Inc. Systems and methods for replacing signal data artifacts in a glucose sensor data stream
US8010174B2 (en) 2003-08-22 2011-08-30 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US9247901B2 (en) 2003-08-22 2016-02-02 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US20190357827A1 (en) 2003-08-01 2019-11-28 Dexcom, Inc. Analyte sensor
US8275437B2 (en) 2003-08-01 2012-09-25 Dexcom, Inc. Transcutaneous analyte sensor
US7774145B2 (en) 2003-08-01 2010-08-10 Dexcom, Inc. Transcutaneous analyte sensor
US8845536B2 (en) 2003-08-01 2014-09-30 Dexcom, Inc. Transcutaneous analyte sensor
US7920906B2 (en) 2005-03-10 2011-04-05 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US20140121989A1 (en) 2003-08-22 2014-05-01 Dexcom, Inc. Systems and methods for processing analyte sensor data
DE10343863A1 (de) 2003-09-23 2005-04-14 Roche Diagnostics Gmbh Verfahren und Gerät zur laufenden Überwachung der Konzentration eines Analyten
US9247900B2 (en) 2004-07-13 2016-02-02 Dexcom, Inc. Analyte sensor
US8774886B2 (en) 2006-10-04 2014-07-08 Dexcom, Inc. Analyte sensor
US8364231B2 (en) 2006-10-04 2013-01-29 Dexcom, Inc. Analyte sensor
EP2316331B1 (en) 2003-12-09 2016-06-29 Dexcom, Inc. Signal processing for continuous analyte sensor
US7534212B2 (en) 2004-03-08 2009-05-19 Nellcor Puritan Bennett Llc Pulse oximeter with alternate heart-rate determination
US7194293B2 (en) 2004-03-08 2007-03-20 Nellcor Puritan Bennett Incorporated Selection of ensemble averaging weights for a pulse oximeter based on signal quality metrics
US8792955B2 (en) 2004-05-03 2014-07-29 Dexcom, Inc. Transcutaneous analyte sensor
US8886272B2 (en) 2004-07-13 2014-11-11 Dexcom, Inc. Analyte sensor
ATE548971T1 (de) * 2004-07-13 2012-03-15 Dexcom Inc Transkutaner analytsensor
US7392075B2 (en) 2005-03-03 2008-06-24 Nellcor Puritan Bennett Incorporated Method for enhancing pulse oximetry calculations in the presence of correlated artifacts
US7590439B2 (en) 2005-08-08 2009-09-15 Nellcor Puritan Bennett Llc Bi-stable medical sensor and technique for using the same
US7657294B2 (en) 2005-08-08 2010-02-02 Nellcor Puritan Bennett Llc Compliant diaphragm medical sensor and technique for using the same
US7869850B2 (en) 2005-09-29 2011-01-11 Nellcor Puritan Bennett Llc Medical sensor for reducing motion artifacts and technique for using the same
US8750953B2 (en) 2008-02-19 2014-06-10 Covidien Lp Methods and systems for alerting practitioners to physiological conditions
US9895068B2 (en) 2008-06-30 2018-02-20 Covidien Lp Pulse oximeter with wait-time indication
US8862194B2 (en) 2008-06-30 2014-10-14 Covidien Lp Method for improved oxygen saturation estimation in the presence of noise
US8855749B2 (en) 2009-09-24 2014-10-07 Covidien Lp Determination of a physiological parameter
US8923945B2 (en) 2009-09-24 2014-12-30 Covidien Lp Determination of a physiological parameter
US8571621B2 (en) 2009-09-24 2013-10-29 Covidien Lp Minimax filtering for pulse oximetry
US8840562B2 (en) 2009-09-24 2014-09-23 Covidien Lp Signal processing warping technique
US8930145B2 (en) 2010-07-28 2015-01-06 Covidien Lp Light focusing continuous wave photoacoustic spectroscopy and its applications to patient monitoring
EP3575796B1 (en) 2011-04-15 2020-11-11 DexCom, Inc. Advanced analyte sensor calibration and error detection
US9770210B2 (en) 2011-09-23 2017-09-26 Nellcor Puritan Bennett Ireland Systems and methods for analyzing a physiological sensor signal
US9833146B2 (en) 2012-04-17 2017-12-05 Covidien Lp Surgical system and method of use of the same
ITUA20161519A1 (it) * 2016-03-10 2017-09-10 Michele Gallamini Dispositivo di monitoraggio di parametri funzionali dell’organismo
CA3235408A1 (en) * 2021-10-18 2023-04-27 Vioptix, Inc. Handheld oximeter with display of real-time, average measurements and status indicator

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4407290A (en) * 1981-04-01 1983-10-04 Biox Technology, Inc. Blood constituent measuring device and method
US4911167A (en) * 1985-06-07 1990-03-27 Nellcor Incorporated Method and apparatus for detecting optical pulses
US4960126A (en) * 1988-01-15 1990-10-02 Criticare Systems, Inc. ECG synchronized pulse oximeter
WO1991011137A1 (en) * 1990-01-30 1991-08-08 Medical Physics Enhanced arterial oxygen saturation determination and arterial blood pressure monitoring
US5355880A (en) * 1992-07-06 1994-10-18 Sandia Corporation Reliable noninvasive measurement of blood gases
JP3116252B2 (ja) * 1992-07-09 2000-12-11 日本光電工業株式会社 パルスオキシメータ

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004503315A (ja) * 2000-06-16 2004-02-05 アイシス・イノベーション・リミテッド 異なるセンサからの測定の組み合わせ
JP4741168B2 (ja) * 2000-06-16 2011-08-03 オックスフォード・バイオシグナルズ・リミテッド 異なるセンサからの測定の組み合わせ
WO2011096240A1 (ja) * 2010-02-05 2011-08-11 日本電気株式会社 生体情報計測器、携帯端末装置、生体情報計測方法およびプログラム

Also Published As

Publication number Publication date
EP1006863B1 (en) 2003-10-15
ATE251864T1 (de) 2003-11-15
CA2283858C (en) 2005-07-26
ES2208876T3 (es) 2004-06-16
WO1998042249A1 (en) 1998-10-01
AU720450B2 (en) 2000-06-01
HK1031189A1 (en) 2001-06-08
EP1006863A1 (en) 2000-06-14
DE69725621T2 (de) 2004-08-19
AU2104997A (en) 1998-10-20
CA2283858A1 (en) 1998-10-01
DE69725621D1 (de) 2003-11-20

Similar Documents

Publication Publication Date Title
JP2001517991A (ja) データ信号適応平均化方法及び装置
US7392075B2 (en) Method for enhancing pulse oximetry calculations in the presence of correlated artifacts
US7991448B2 (en) Method, apparatus, and system for removing motion artifacts from measurements of bodily parameters
US6519486B1 (en) Method, apparatus and system for removing motion artifacts from measurements of bodily parameters
US6035223A (en) Method and apparatus for determining the state of an oximetry sensor
CN106889979B (zh) 一种基于心电信号和血氧容积波的连续性血压测量系统
US20130253331A1 (en) Method and apparatus for estimating pulse rate
JP2007527772A (ja) 信号品質測定基準に基づいたパルス酸素濃度計のためのアンサンブル平均重みの選択
JP2007527770A (ja) 交代性の心拍数決定を有するパルス酸素濃度計
JP2001517129A (ja) 調和的にデータをフィルター化するための方法及び装置
CA2623270C (en) Signal processing for pulse oximetry
US20100249551A1 (en) System And Method For Generating Corrective Actions Correlated To Medical Sensor Errors
JP4571220B2 (ja) データ信号適応平均化方法及び装置
EP1237470A2 (en) Method and system for improving photoplethysmographic analyte measurements by de-weighting motion-contaminated data
JP2001518822A (ja) パルスレート及び飽和度を測定するための方法及び装置
AU720238C (en) Method and apparatus for harmonically filtering data
EP4037556A1 (en) Method and system for evaluating the quality of a physiological signal

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040224

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20060117

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20060413

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20060605

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060713

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060718

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20060905

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070410

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20080425

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20080605

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20090116

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20090121

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20090409

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20090414

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090513