JP2007527771A - 酸素飽和および心拍数のための個別の集団平均を用いるパルス酸素濃度計 - Google Patents

酸素飽和および心拍数のための個別の集団平均を用いるパルス酸素濃度計 Download PDF

Info

Publication number
JP2007527771A
JP2007527771A JP2007502899A JP2007502899A JP2007527771A JP 2007527771 A JP2007527771 A JP 2007527771A JP 2007502899 A JP2007502899 A JP 2007502899A JP 2007502899 A JP2007502899 A JP 2007502899A JP 2007527771 A JP2007527771 A JP 2007527771A
Authority
JP
Japan
Prior art keywords
pulse
waveform
metric
low pass
pass filter
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
JP2007502899A
Other languages
English (en)
Other versions
JP4700682B2 (ja
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 JP2007527771A publication Critical patent/JP2007527771A/ja
Application granted granted Critical
Publication of JP4700682B2 publication Critical patent/JP4700682B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/725Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6843Monitoring or controlling sensor contact pressure

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Biomedical Technology (AREA)
  • Medical Informatics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Artificial Intelligence (AREA)
  • Optics & Photonics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • General Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

酸素飽和および脈拍数の算出における使用のために、検出された波形を処理するための2つの別個の集団平均算出器の使用。酸素飽和を算出するために使用される集団平均算出器が、正規化された信号において動作する一方、脈拍数算出のための集団平均算出器は、正規化されていない信号において動作する。2つの集団平均算出器を介して2つの経路のために選ばれる測定基準は、酸素飽和または脈拍数算出のための集団平均を最適化するために変化され得る。

Description

本発明は、酸素濃度計に関し、より詳細にはパルス酸素濃度計からの検出された波形におけるパルスの集団平均に関する。
パルス酸素濃度計は典型的には、動脈血におけるヘモグロビンの血中酸素飽和度、組織に供給する個々の脈拍の血液量、および、患者の個々の心拍に対応する脈拍の速度を含む(それらに限定されるわけではない)血液の様々な化学的特性を計測するために使用される。これらの特性の計測は、非侵襲センサの使用によって達成される。その非侵襲センサは、血液が組織を灌流する患者の組織の一部に光を散乱し、そのような組織における光の吸収を光電測光にて感知する。様々な波長において吸収される光の量は、計測される血液成分の量を計算するために使用される。
組織に散乱された光は、血液中に存在する血液成分の量を示す量において、血液によって吸収された一つ以上の波長として選択される。組織に散乱された、透過光線の量は、組織および関連する光の吸収における血液成分の変化する量に従って変化する。血中酸素レベルを計測するために、そのようなセンサは、典型的には、血中酸素飽和度を計測するための周知の技術に従い、少なくとも二つの異なる波長の光を生成するように適合される光源、および、それらの波長の両方に感度の高い光検出器を提供される。
周知の非侵襲センサは、指、耳、または頭皮などの体の一部に固定される装置を含む。動物や人間においては、これらの体の部分の組織は血液によって灌流され、その組織の表面はセンサに対して、容易にアクセス可能である。
(N−100)。およそ1985年頃に始まったN−100技術は、パルスのサイズ、パルスの形、生じる予定時間(周波数)、およびR/IRの比率のパルス履歴に基づいて、パルスを許可または拒絶した。
特にN−100は、信号最大値、次いで最大負スロープの地点、次いで最小値を捜すことによってパルスを見い出した。処理は、「マンチ(munch)」として言及される状態機械において行われた。各最大値は、信号がノイズゲートとして言及されるノイズしきい値の下を通過するまで認定されなかった。これは、異なる、期待された信号振幅に適用するためにノイズゲートレベルが次の処理ステップからのフィードバックによって設定されたため、適応フィルタとして作用した。パルスは次いで、「レベル3」処理において許可または拒否される。ここにおいて、「レベル3」処理は、新しいパルスの振幅、周期、および比率の比率(ratio−of−ratios)(RedおよびIRがAC−DC比として表されるRed−IR比)を、履歴バッファにおける値の平均と比較し、次いで差異が信頼水準内にあるかを否かを決定することによって変化する信号に適応されるフィルタであった。新しいパルスが許可された場合、履歴バッファは新しいパルスに対する値を用いて更新された。中心周波数および帯域幅(信頼限界)がフィルタの出力からのフィードバックによって適応され、レベル3処理は、適応帯域通過フィルタとして作用した。
(N−200)。N−200は、ECGと同期でき、ECGフィルタリングを含むため、N−100を改良したものである。N−200は、パルス最大値の時間の計測と最小値の時間の計測との間のベースラインシフトを補償するために補間も追加した。N−200は、信号サンプルの変化する数の平均を計算する「ボックスカー」フィルタのような、他のフィルタリング特徴を含んだ。
様々なフィルタリングおよびスケーリングステップの後、N−200は、N個のサンプルの平均を計算する「ボックスカー(boxcar)」フィルタにデジタル化された信号を適用する。ここにおいて、Nは、フィルタされた心拍数に従って、次の処理ステップからのフィードバックによって設定される。新しいサンプルがボックスカーフィルタに平均される一方、最も古いサンプルはドロップされる。ボックスカーの長さ(N)は、パルスしきい値、絶対最小パルス、および小パルスの3つのパラメータを設定するために使用される。集団平均(別名「スライダ」)フィルタは、次いで新しいサンプルと1つ早いパルス周期(period)からの前の集団平均との重み付け平均を生成する。サンプルは、次いでN−100のような「マンチ」状態機械およびノイズゲートに渡される。補間特徴は、ベースラインレベルにおける変化を補償するために、N−100処理に追加される。最小値および最大値が異なった時間において生じるため、変化するベースラインは、最大値ではなく最小値を増やしまたは減らさせ得、その逆も同様にあり得る。
「集団平均(ensemble averaging)」は、複合パルスを形成するために複数のパルスからのサンプルを一緒に平均する処理のためのNELLCORのトレードマークであるC−Lockの整数部である。この処理は、「心臓ゲート平均」としても知られる。これは、各パルスの開始をマークするために「トリガ」イベントを要求する。
Conlonによる特許文献1は、異なる重みが異なるパルスに割り当てられ、複合、平均パルス波形が酸素飽和を算出するために使用される、集団平均を開示する。上述されたN−100は、特許文献2に記載されている。N−200の局面は、特許文献3(Corenman)および特許文献4(Stone)に記載されている。
米国特許第4,690,126号明細書 米国特許第4,802,486号明細書 米国特許第4,911,167号明細書 米国特許第5,078,136号明細書
本発明は、酸素飽和および脈拍数の算出における使用のために、検出された波形を処理するための2つの別個の集団平均算出器の使用に向けられる。酸素飽和を算出するために使用される集団平均算出器が、正規化された信号において動作する一方、脈拍数算出のための集団平均算出器は、正規化されていない信号において動作する。両方の波長に対応する波形が、酸素飽和計算のための比率の比率を保存するために、IRパルス振幅のような同等の量によって正規化される必要があることに留意されたい。
脈拍数のための正規化を用いない信号の使用は、動作アーチファクトのような生理学的パルスより実質的に大きいアーチファクトを不認定するためのソフトウェアの能力を改良する。脈拍数のための正規化を用いない信号の使用は、パルスが正規化によってミスされることを回避する。
2つの集団平均算出器を介して2つの経路のために選ばれる測定基準は、酸素飽和または脈拍数算出のための集団平均を最適化するために変化され得る。例えば、脈拍数を算出するために使用される場合、酸素飽和を算出するときと比較して、より低いしきい値が不整脈パルスを検出するために測定基準に対して使用される。更に、短期パルス振幅比のための測定基準は、動作アーチファクトがちょうど平常の状態に戻ったとき小さくなり、これは酸素飽和算出より多くの重みを脈拍数算出において与えられる(短期パルス振幅比は、現行のパルス振幅/前のパルス振幅である)。
図1は、本発明を具体化する酸素測定システムの実施形態を示す。センサ10は、赤および赤外線LED、ならびに光検出器を含む。これらは、ケーブル12によって、基板14へ接続される。LED駆動電流は、LED駆動インターフェース16によって提供される。センサからの、受信された光電流は、I−Vインターフェース18に提供される。IRおよび赤電圧は、次いで、本発明を具体化する、シグマ−デルタインターフェース20に提供される。シグマ−デルタインターフェース20の出力は、10ビットA/D変換器を含むマイクロコントローラ22に提供される。コントローラ22は、プログラムのためのフラッシュメモリ、および、データのためのEEPROMメモリを含む。プロセッサはまた、フラッシュメモリ26に接続されるコントローラチップ24を含む。最後に、クロック28が使用され、センサ10におけるデジタル較正へのインターフェース30が提供される。別個のホスト32は、処理された情報を受信し、および、アナログ表示を提供するためのライン34上にて、アナログ信号を受信する。
(設計概要)本発明の設計は、不要ノイズを解決することを意図する。信号測定基準は計測され、フィルタ重み付けを決定するために使用される。信号測定基準は、パルスがおそらく周波数(人間の心拍数の範囲内であるか)、形(人間の心臓パルスのような形であるか)、立ち上がり時間等のようなプレチスモグラフまたはノイズであるか否かを示すものである。同様な技術が、この出願の背景に記載されたNellcor N200において使用された。新しい設計は、本発明の特許請求の範囲に記載されるように、2つの集団平均算出器の使用のような多数の異なる特徴および変化を追加する。
アーキテクチャの詳細は、図2に示される。この設計は、以下に別々に説明される、酸素飽和および脈拍数の両方を算出する。
I.(酸素飽和算出)
A.(信号調整)−デジタル化された赤およびIR信号は受信され、(1)第1の導関数(derivative)がベースラインシフトを取り除くように作用し、(2)固定係数を用いて低域通過フィルタリングをし、(3)比率を保存するためにDC値で割ることによって、このブロックにおいて調整される。信号調整サブシステムの機能は、人間プレチスモグラフにおいて生じるより高周波数を強調し、動作アーチファクトが通常集中している低周波数を減衰するためである。信号調整サブシステムは、初期化の間に識別されるハードウェア特性に基づいて、そのフィルタ係数(広帯域または狭帯域)を選択する。
入力−デジタル化された赤およびIR信号。
出力−前処理された赤およびIR信号。
B.(パルス識別(Identification)および認定(Qualification))−低域通過フィルタされたおよびデジタル化された赤およびIR信号は、パルスを識別し、それらを本当らしい動脈拍として認定するためにこのブロックに提供される。これは、予め訓練されたニューラルネットワークを使用して行われ、本来はIR信号において行われる。Nellcor N−100において行われたように、パルスは、その振幅、形、および周波数を調査することによって識別される。このブロックへの入力は、ブロックDからの平均パルス周期である。この機能は、脈拍数を使用して目立つ(upfront)認定を変化させたN−100と同様である。出力は、不整脈度および個別のパルス品質を示す。
入力−(1)前処理された赤およびIR信号、(2)平均パルス周期、(3)低域通過フィルタからの低域通過波形。
出力−(1)不整脈度、(2)パルス振幅変化、(3)個別のパルス品質、(4)パルスビープ通知、(5)認定されたパルス周期および段階(エイジ)。
C.(信号品質測定基準を計算する)−このブロックは、パルスの形(導関数スキュー)周期変化性、パルス振幅および変化性、比率の比率変化性、ならびに脈拍数に関する周波数内容を決定する。
入力−(1)未処理のデジタル化された赤およびIR信号、(2)不整脈度、個別のパルス品質、パルス振幅変化、(3)前処理された赤およびIR信号、(4)平均パルス周期。
出力−(1)低域通過および集団平均フィルタ重み、(2)センサオフ検出器のための測定基準、(3)正規化され前処理された波形、(4)変調率。
D.(平均パルス周期)−このブロックは、受け取られたパルスから平均パルス周期を算出する。
入力−認定されたパルス周期および段階。
出力−平均パルス周期。
E1.(低域通過フィルタおよび集団平均)−ブロックE1は、脈拍数識別のために、ブロックAによって調整され、ブロックCによって正規化された信号を低域通過フィルタし、平均を集合させる。低域通過フィルタに対する重みは、信号測定基準ブロックCによって決定される。信号は集団平均もされ(これは、脈拍数およびその高調波の近い対象のもの以外の周波数を減衰する)、集団平均フィルタ重みも信号測定基準ブロックCによって決定される。信号が低下するとともにフラグされた場合、より少ない重みが割り当てられる。不整脈の間において集団平均が適切ではないため、信号が不整脈としてフラグされた場合、より多い重みが割り当てられる。赤およびIRは別々に処理されるが、同等のフィルタリング重みを用いて処理される。フィルタリングは、信号測定基準が最初に算出されることを可能にするため、約1秒遅らされる。
フィルタは、連続的に可変の重みを使用する。サンプルが集団平均されない場合、前のフィルタされたサンプルに対する重み付けは、重み付けされた平均においてゼロに設定され、新しいサンプルは依然としてコードを介して処理される。このブロックは、信号の段階−フィルタリングの累積された量(応答時間および処理における遅れの合計)を追跡する。古すぎる結果はフラグされる(良いパルスがしばらくの間検出されていない場合)。
入力−(1)正規化され、前処理された赤およびIR信号、(2)平均パルス周期、(3)低域通過フィルタ重みおよび集団平均フィルタ重み、(4)利用可能な場合、ECGトリガ、(5)ゼロ交差トリガに対して、IR基本波。
出力−(1)フィルタされた赤およびIR信号、(2)段階。
F.(フィルタ波形相関を推定および平均重みを計算する)−これは、上述されたN100およびN200において使用されたものと同様のノイズ測定基準を使用し、フィードバックを使用しない。フィルタに対する可変の重み付けは、比率の比率分散によって制御される。この可変の重みフィルタリングの効果は、比率の比率が、アーチファクトが増えるとともにゆっくりと変化し、アーチファクトが減るとともに速く変化することである。サブシステムは2つの応答モードを有する。ファストモードにおけるフィルタリングは、3秒の段階測定基準をターゲットする。ターゲット段階は通常モードにおいては5秒である。ファストモードにおいては、現在の値の最小重み付けは、より高いレベルにおいてクリップされる。言い換えると、ノイズが存在した場合、低い重みは最も新しい比率の比率算出に割り当てられ、ノイズが存在しない場合においては高い重みが割り当てられる。
入力−(1)フィルタされた赤およびIR信号および段階、(2)較正係数、(3)応答モード(ユーザスピード設定)。
出力−比率の比率算出に対する平均重み。
H.(飽和を算出する)−飽和は、較正係数および平均された比率の比率のアルゴリズムを使用して算出される。
入力−(1)平均された比率の比率、(2)較正係数。
出力−飽和。
II.(脈拍数算出)
E2.(低域通過フィルタおよび集団平均)−ブロックE2は、脈拍数識別のために、ブロックAによって調整された信号を低域通過フィルタし、平均を集合させる。低域通過フィルタに対する重みは、信号測定基準ブロックCによって決定される。信号は、集団平均もされ、(これは、脈拍数およびその高調波の近い対象のもの以外の周波数を減衰する)、集団平均フィルタ重みも信号測定基準ブロックCによって決定される。信号が低下するとともにフラグされた場合、より少ない重みが割り当てられる。不整脈の間においてフィルタリングが適切ではないため、信号が不整脈としてフラグされた場合、より多い重みが割り当てられる。赤およびIRは別々に処理される。このブロックの処理は、信号測定基準が最初に算出されることを可能にするため、約1秒遅らされる。
フィルタは、連続的に可変の重みを使用する。サンプルが集団平均されない場合、前のフィルタされたサンプルに対する重み付けは、重み付けされた平均においてゼロに設定され、新しいサンプルは依然としてコードを介して処理される。このブロックは、信号の段階−フィルタリングの累積された量(応答時間および処理における遅れの合計)を追跡する。古すぎる結果はフラグされる(良いパルスがしばらくの間検出されていない場合)。
入力−(1)前処理された赤およびIR信号、(2)平均パルス周期、(3)低域通過フィルタ重みおよび集団平均フィルタ重み、(4)利用可能な場合、ECGトリガ、(5)ゼロ交差トリガに対して、IR基本波。
出力−(1)フィルタされた赤およびIR信号、(2)段階。
I.(フィルタされたパルス識別および認定)−このブロックは、フィルタされた波形からパルス周期を識別および認定し、その結果はパルスがブロックBによって不認定された場合にのみ使用される。
入力−(1)フィルタされた赤およびIR信号および段階、(2)平均パルス周期、(3)ハードウェアIDまたはノイズフロア、(4)センサの種類。
出力−認定されたパルス周期および段階。
J.(平均パルス周期および脈拍数を算出する)−このブロックは、脈拍数および平均パルス周期を算出する。
入力−認定されたパルス周期および段階。
出力−(1)平均パルス周期、(2)脈拍数。
III.(静脈拍)
K.(静脈拍を検出する)−ブロックKは、ブロックAから前処理された赤およびIR信号および段階を入力として、ならびに脈拍数を受け取り、出力として静脈拍の表示を提供する。このサブシステムは、集団平均フィルタに出力される一枚歯くし型フィルタを使用して、時間ドメインにおいてIR基本波波形を生成する。
入力−(1)フィルタされた赤およびIR信号および段階、(2)脈拍数。
出力−静脈拍表示、IR基本波
IV.(センサオフ)
L.(センサオフおよびパルス振幅の損失を検出する)−パルス損失およびセンサオフ検出サブシステムは、センサが患者からオフされているか否かを決定するために予め訓練されたニューラルネットを使用する。ニューラルネットへの入力は、最後の数秒間に渡ってIR値および赤値の挙動のいくつかの局面を量化する測定基準である。サンプルが、多数のアルゴリズムのサブシステムによって無視される一方、信号状態は、パルス存在またはセンサがオフかもしれないのどちらでもない。信号状態の変数の値は、「パルス存在、未接続、パルス損失、センサがオフかもしれない、センサオフ」である。
入力−(1)測定基準、(2)フロントエンドサーボ設定およびID
出力−センサオフ表示を含む信号状態
集団平均サブシステム
(集団平均サブシステム)−集団平均サブシステムの機能は、可変の重み付けを用いてその入力ストリームをフィルタし、ノイズまたは動作アーチファクトによってあまり変形されていない波形を出力することである。フィルタされた波形におけるアーチファクトを減らすことは、動作またはノイズの間によりロバストな飽和またはレート推定を可能にする。
集団平均サブシステムは、数パルスのスパンに渡って、ゼロ平均であるあらゆるサンプルに対してIRおよび赤入力を要求する。
入力サンプルは、信号測定基準サブシステムから受け取られた重み(LPF_Weight)を用いて、最初にIIR低域通過フィルタされ、一秒遅れバッファに格納される。
このサブシステムは、現行の複合パルスのi番目のサンプルを形成するために、現行の1秒遅れ入力パルスのi番目のIRおよび赤サンプルを、前の複合パルスのi番目のサンプルと平均する。パルス周期の始まりを開始するためのトリガは、RWave_Occurred入力および平均周期入力(Optical_Period)から派生する(優先順位において)。現行のサンプルに与えられる重み対前のパルスの対応するサンプルは、信号測定基準サブシステムから受け取られるEnsemble_Averaging_Weight値の値によって決定される。
複合パルスは、フィルタの入力パルスと比較してノイズまたは動作アーチファクトによる変形が少ない。図3は、複合パルスを形成するためにパルスがどのようにして一緒に平均されるかについての概念図である。
サブシステムは、心拍と同期するべきであるトリガを受け取りおよび認定する。トリガは、利用可能な場合、R波認定サブシステムからの認定されたR波トリガである。R波トリガが利用不可能な場合、トリガは、脈拍数算出サブシステムからの平均周期入力(Optical_Period)から内面的に生成される。「パルス」は、各々の認定されたトリガにおいて開始し、次の認定されたトリガにおいて終了することが考慮される。このようにして連続トリガは、集団平均周期を規定するために使用される。
図4は、動作アーチファクトによって損なわれたパルスの一連へのサブシステムの応答の代表を示す。垂直線はR波トリガである。フィルタされた出力は、入力パルスのおよその形およびサイズを再格納する。平均の量は、信号測定基準サブシステムから受け取られるEnsemble_Averaging_Weightによって決定されるように、動作アーチファクトが増えるとともに増える。
図4は、集団平均サブシステムの可変の重み付けの実施例である。R波トリガは、各パルスの始まりをマークする。IR入力線は、動作アーチファクトによって損なわれたパルスを示す。フィルタされたIR線は、オリジナルのパルスサイズおよび形を大きく再格納する複数のパルスの複合である。フィルタされたIRは、IR入力から1秒遅らされる。
「低域通過フィルタ」−IRおよび赤入力波形は、信号測定基準サブシステムから受け取られた以下のような重み付け(LPF_Weight)を使用してIIRフィルタされる。
Lowpass_Waveforms=Lowpass_Waveformst−1+LPF_Weight(Input_Waveforms+Lowpass_Waveformst−1
サブシステム初期化の間、重みは1.0のデフォルトに設定される。
「1秒遅れバッファ」−IRおよび赤Lowpass_Waveforms(それらの関連する段階および状態とともに)ならびにRWave_Occurred入力は、1秒長いバッファにおいて格納される(IR_Inputs,Red_Inputs,Input_Valid,Age_Inputs,RWave_Inputs)。
「タイムスタンプ」−タイムスタンプ(Current_Time_ctr)は単に、0に初期化され、毎サンプル周期に増加される32ビットカウンタである。いくつかの受け取られた値は、それらの段階を再構成するために、それらのタイムスタンプと格納される。
「集団重みバッファ」−受け取られた最後の4つのEnsemble_Averaging_Weightsおよびそれらのタイムスタンプは、Ensemble_Average_BufferおよびEnsemble_Weight_Timestampに格納される。これは、1秒遅れバッファにおける毎エントリが、その正しい重みに関連することを可能にする。遅れ入力が1秒遅れバッファから引き出されるたびに、関連する集団重みは、少なくとも遅れサンプルと同じくらい最近である、このバッファにおける最も古い重みに設定される(weight timestamp+one second≧Current_Time_ctr)。
「トリガ認定」−サブシステムは、RWave_OccurredまたはOptical_Periodの2つの入力のうちの1つからそのトリガを選択する。RWave_Occurredがデフォルトトリガである。RWave_Occurredがトリガとして使用される前に、遅れIRおよび赤サンプルと同期するため、上述された1秒遅れバッファを通過する。脈拍数算出サブシステムからのOptical_Period入力に基づくトリガは、R波トリガが少なくとも5秒間受け取られていない場合にのみ認定される。この待ち時間は、R波トリガが利用不可能であることを決定するのに十分でないとされる。次いで、Optical_Periodからの第1のトリガは、サブシステムのOptical_Trigger_WF入力波形の第1のゼロ交差まで遅らされる。次に、トリガはOptical_Periodから単独で派生される。
「集団平均モデル」−サブシステムは、毎パルスの開始を示すためにIRおよび赤入力サンプルを受け取る。IRおよび赤入力が既に低域通過フィルタされ、1秒送れバッファを通過したことに留意されたい。これは、現行のパルスのi番目のサンプルを前の推定されたパルスのi番目のサンプルと平均することによって、その現行の入力サンプルの拍動成分を推定する。
フィルタ出力は、従って、複数のパルスの複合であって、各サンプルが受け取られるたびに算出される。フィルタは、脈拍数およびその高調波においてまたはその近くの周波数のみを通過するくし型フィルタの周波数応答を有する。平均の量は、くし型フィルタの「歯」の幅を決定する。
フィルタは、現行パルスのi番目のサンプル値が前のパルスのi番目のサンプルと概略で均等であることを想定する。iは、各々の認定されたトリガにおいてゼロの値を有し、各々の次のサンプルにおいて増加する、時間tのランプ関数である。
「集団平均フィルタ等式および中間変数」−集団平均フィルタは、信号測定基準サブシステムによって供給されるEnsemble_Averaging_Weight,w,を使用する。以下の式は、毎パルスのi番目のサンプルにおいて実行されるべき基本ステップを示す。
Figure 2007527771
ここで、z´は、1つ前のパルスからのzの値を意味する。zは、サブシステムの集団フィルタされた出力である。全てのzは、複合パルスバッファに格納される。xは、1秒遅れバッファからの最も最近の出力サンプルである。式1における全ての変数は、スカラである。
「複合パルスバッファ」−複合IRおよび赤パルスは、別々の複合パルスバッファに格納される。インデックス、i、は、トリガが受け取られ、各サンプルに対して増加された場合、バッファの初めにリセットされる。
バッファは、脈拍数変化性を可能にするために、1つの20BPMパルスに加えて10%のマージンを格納するのに十分に長い必要がある。したがって、バッファは、20BPMにおける少なくとも1.1連続複合パルスおよび最も多い脈拍数における2つの複合パルスを用いて更新される。iがバッファの終わりを通過される場合、フィルタの出力がその入力と同一に設定される時間の間、通常処理は次のトリガまで中断される必要がある。
トリガ間の間隔が長くなると、複合パルスバッファは、パルスの終わりにおいて平均するために、最近のサンプルを含み得ない。
複合パルスのj番目のサンプルは、算出された各サンプルであって、j=i+mである。mは、現行と前の認定されたトリガとの間、すなわちパルス周期、のサンプルの数である。jが、iのように、tのランプ関数であることに留意されたい。
jに対しては、式(1)は、
Figure 2007527771
を読み取るように修正される。
wの同等の値は、zおよびzを算出するために使用される。jがバッファの終わりに届いた場合、zの算出は、jが再び有効になるまで中断される必要がある。
パルス周期を変化させることは、各トリガが受け取られた後、サブシステム出力における小さな不連続を引き起こし得る。トリガの後の最初の4つのサンプルの間における特別な処理は、この効果を減少させる。フィルタされた出力サンプルは、標準出力(式1)と第2の複合パルス(式1b)との間に補間される。それによって、フィルタされた出力のそれぞれは、これらのサンプルにおけるzの80%、60%、40%および20%である。第2の複合パルスが利用不可能な場合、入力波形は、適したところにおいて使用される。
「初期化、再初期化、クリアリング、および無視されたサンプル」−2つのトリガに加えて1秒サンプルが受け取られるまで、jは無意味であり、zは算出されない。
サブシステムは、処理における割込みから回復するための2つ方法を含む。5秒より多い時間がトリガ間において経過した場合、サブシステムは次の認定されるトリガにおいて「クリア」される。クリアリング動作は、サブシステムの持続性変数の全てをそれらの初期値に設定するが、トリガ認定に需要であり、1秒遅れバッファを維持する以下の変数は例外とする。
1.サブシステムの現行のタイムスタンプ
2.最後のR波トリガおよび最後の認定されたトリガ以来の経過時間
3.ゼロ交差を検出するために使用される、前のOptical_Trigger_Waveformサンプル
4.トリガ認定状態機械の状態
5.1秒遅れ入力バッファ
6.Ensemble_Averaging_Weightバッファ
この「クリアリング」動作の状態がまれであり、処理における長い割込みを示し得るため、この動作は実行される。フィルタ出力は、iが複合パルスバッファの終わりをオーバーフローするたびに、フィルタ入力とも同一である。
5秒より少ない時間における有効入力の欠如によって、サンプルが無視された(処理されない)場合、複合パルスバッファは、変化しないまま残り、サブシステムの出力は無効とマークされる。処理が再開し、Optical_Periodがゼロではない場合、バッファのインデクスは、整数のパルス周期の前になるはずだったものにリセットされ、1秒遅れバッファはリセットされる。これは、複合パルスバッファが、処理における短い割込みの後に現行のパルスの正確な表示を更に含むべきであるために行われる。Optical_Period推定がゼロ(無効)またはサンプルが少なくとも5秒間無視された場合、サブシステムは「クリア」される。ほとんどの「無視されたサンプル」は、通常2秒より少なくかかる、酸素濃度計のLED明るさまたは増幅器ゲインにおける調整によるように期待される。インデクスをリセットするための式は以下の通りである。
Figure 2007527771
j=New_Idx+j−i_Idx、提供されたjはまだオーバーフローせず、複合パルスバッファをオーバーフローしない
i_Idx=New_Idx、提供されたi_Idxはまだオーバーフローせず、複合パルスバッファをオーバーフローしない
パルス損失およびセンサオフ検出サブシステムが、パルスが長期間休んだ後に再取得されたことを決定した場合、またはセンサが接続された場合、集団平均サブシステムは再初期化される。「再初期化」は、サブシステムの持続性変数の全てをそれらの初期値に設定することを意味する。これは、再初期化を呼び起こす両方のイベントが、集団平均サブシステムの前のパルス表示がもう現行のものではないだろうとするために行われる。
最後のR波トリガおよび最後の認定されたトリガ以来の経過時間間隔は、無視されたサンプルにおいて増加される。無視されたサンプルの間に生じるR波またはゼロ交差は、トリガを認定するために使用されない。
図5〜図7の状態遷移図は、複合パルスバッファおよびそのインデクスにおいて維持される2つのパルスの各々を更新するための状態機械を示す。
「段階測定基準(Age Metric)」−サブシステムは、IRおよび赤入力のサンプル、Age_Inputs、における段階を受け取り、複合パルス出力の段階、Age_Out、を出力する。各サンプルに対して、Age_Outは、最後に更新され、Age_Inputsで平均されて以来の経過時間によって、複合波形を更新するために使用される(Ensemble_Averaging_Weight)、w、フィルタ重みを使用して増加される。Age_Outのための公式は、
Figure 2007527771
であり、mはAge_Outが最後に更新されて以来のサンプルの数であり、Nは1秒におけるサンプルの数である。各Age_Outが更新されるタイムスタンプは、mを算出するために格納される必要がある。
サブシステムは、Age_Outも更新する必要がある。Age_Outのための公式は、jがiのために置き換えられる以外は上記されたものと同一である。
サブシステムは、無視されたサンプルを含み、毎サンプルごとに現行のタイムスタンプ(Current_Time_Ctr)を増加する必要がある。
サブシステムがクリアされ、または再初期化された場合、Age_Outバッファにおける全てのエントリは、Age_Inputsに再初期化される。更に、Age_Outエントリが更新されたタイムスタンプに対するバッファにおける全てのエントリは、現行のタイムスタンプから1引いて設定される。これらのステップは、Age_Out値が、サブシステムがクリアされ、再初期化されて以来の経過時間より新しいことを保証するためにとられる。
波形がパルス認定および脈拍数のために使用される集団平均インスタンスは、単独で周波数コンテンツに依存するRate_LPF_Weightを使用する。波形が比率の比率および飽和を算出するために使用される集団平均インスタンスは、RoR_Variance測定基準が低域通過フィルタリングの追加によって良く(低く)なるか否かにも依存するSat_LPF_Weightを使用する。これらの重みは、0.1から1.0の間にあり、任意の単一のステップにおいて0.05より多く増えない。
集団平均重み
サブシステムに、パルス識別および認定サブシステムがポテンシャルパルスの評価をちょうど完成したことが通知された場合、サブシステムは、集団平均サブシステムのインスタンスによって使用される集団平均重みを更新する。個別の重みは、その出力が飽和および脈拍数を計算するのに使用される、2つの集団平均インスタンスのために計算される。これらの重みは、その入力波形が集団平均されていないパルス識別および認定サブシステムのインスタンスによって提供される測定基準に部分的に基づく。
Sat_Ensemble_Averaging_Filter_Weightのための式は以下の通りである。
Figure 2007527771
ここで、バウンド(a,b,c)はmin(max(a,b),c)を意味する。
上記の式は、比率の比率分散の低い値に対する0.5のデフォルト重みという結果となる。Short_RoR_VarianceおよびPulse_Qual_RoR_Varianceは両方とも、3秒間隔に渡って計算される。Pulse_Qual_RoR_Varianceに対する間隔は、通常最も最近のサンプルを含む、最も最近のパルスの認定または拒絶を用いて終了する。重みは、高い比率の比率分散および一般的に動作アーチファクトを示すLong_Term_Pulse_Amp_Ratioの高い値によって減少される。Arr_Min_Filt_Wt_For_Satは、不整脈度を量化するPeriod_Varに本来基づいて、集団平均重み(0.05〜0.55の範囲)において最小値を組み付ける(impose)。これは、集団平均が異なる周期を有するパルスに対して効果的ではないために行われる。最も最近のパルスが良いPulse_Qual_Scoreを受け取った場合、これはSat_Ensemble_Averaging_Filter_Weightの最大値を0.5から1.0に増やし得る。
Rate_Ensemble_Averaging_Filter_Weightのための式は以下の通りである。
Figure 2007527771
これらの式は以下の通りのように、Sat_Ensemble_Averaging_Filter_Weightのためのものと異なる。
a)Arr_Probを計算するために使用されるしきい値は、不整脈パルスがパルス認定の前の集団平均によって不明瞭にならないことが望まれるため、多少低い。
b)Short_Term_Pulse_Amp_Ratioの小さい値は、動作アーチファクトがちょうど平常の状態に戻ったことを一般的に示し、それは集団平均重みが早く増やされ得ることを意味する。これは、比率の比率フィルタリングおよび飽和計算に対してではなく、パルス認定に対して都合良くするために経験的に見い出された。
c)心臓が、前の不整脈を用いてまたは用いずにビートを抜かした場合、結果として生じる平均より長いQualified_Pulse_Periodは、次のパルス認定から抜かされたビートを不明瞭にしないために、集団平均重みを増やす。
(定義)
[データ入力]
「Avg_Period」−脈拍数算出サブシステムによって報告される平均パルス周期
Long_Term_Pulse_Amp_Ratio−履歴のパルス振幅と比較された最後のパルス振幅を量化する。パルス識別および認定サブシステムによって提供される。1.0より実質的に大きい値は、動作アーチファクトを一般的に示し、より低いEnsemble_Averaging_Filter_Weightsの結果として生じる。
Period_Var−パルス識別および認定サブシステムからの周期変化性測定基準。不整脈の範囲をゲージするために使用される。
Pulse_Qual_RoR_Variance−パルス識別および認定サブシステムからのRoR_Variance測定基準。例えば、0.10の値は、連続パルス周期の間の平均差がAvg_Periodの10%であることを示す。
Pulse_Qual_Score−パルス識別および認定サブシステムにおけるパルス認定ニューラルネットによって計算されるスコア。ゼロは非常に悪く、1.0は非常に良い。
Qualified_Pulse_Period−パルス識別および認定サブシステムによって認定される最も最近のパルス周期。
Short_Term_Pulse_Amp_Ratio−前のパルス振幅と比較された最後のパルス振幅を量化する。
[出力]
Frequency_Ratio−Mean_IR_Frequency_Content−脈拍数比
LPF_RoR_Variance−比率の比率の変化性を量化する。LPF_Scaled_Waveformsから9秒ウィンドウに渡って計算される。
Rate_LPF_Weight−パルス認定および脈拍数算出のために使用される波形を前処理する集団平均サブシステムのインスタンスによって使用される低域通過フィルタ重み。
RoR_Variance−比率の比率の変化性を量化する。Scaled_Waveformsから9秒ウィンドウに渡って計算される。0.10の値は、サンプル−サンプル比率の比率値が、平均比率の比率値と、平均比率の比率値の10%平均異なる。
Sat_Ensemble_Averaging_Filter_Weight−パルス認定および脈拍数算出のために使用される波形を前処理する集団平均サブシステムのインスタンスによって使用される集団平均重み。
Sat_LPF_Weight−パルス認定および脈拍数算出のために使用される波形を前処理する集団平均サブシステムのインスタンスによって使用される低域通過フィルタ重み。
Scaled_Waveforms−IRおよび赤Pre_Processed_Waveformsのスケールされたバージョン。
Short_RoR_Variance−比率の比率の変化性を量化する。Scaled_Waveformから3秒ウィンドウに渡って計算される。
[内部変数]
Arr_Prob−集団平均の量を制限する不整脈の可能性。2つのEnsemble_Averaging_Filter_Weightsの各々に対して特有であるしきい値を用いて、Period_Varに基づく。
Arr_Min_Filt_Wt_For_Rate,Arr_Min_Filt_Wt_For_Sat−それぞれのArr_Prob値に基づく2つのEnsemble_Averaging_Filter_Weightsに対する最小値。
LPF_Scaled_Waveforms−LPF_RoR_Varianceを計算するために使用されるScaled_Waveformsの低域通過フィルタバージョン。
Mean_IR_Frequency_Content−IR入力波形の平均周波数コンテンツの推定。
RoR_Variance_Based_Filt_Wt−RoR_Variance測定基準およびLong_Term_Pulse_Amp_Ratioに基づくEnsemble_Averaging_Filter_Weightsに対する構成要素。
本発明を具体化する酸素測定システムのブロック図である。 本発明の実施形態を含む、酸素濃度計のソフトウェア処理ブロックの図である。 複合パルスの生成を示す図である。 集団平均実績のチャートである。 複合パルスバッファにおける所定の変数を更新するための状態機械の図である。 複合パルスバッファにおける所定の変数を更新するための状態機械の図である。 複合パルスバッファにおける所定の変数を更新するための状態機械の図である。

Claims (18)

  1. 酸素飽和および脈拍数を決定するためにパルス酸素濃度計における信号を処理する方法であって、該方法は、
    患者からの2つの異なる光の波長に対応する波形を受け取ることと、
    第1の集団平均算出器において、該波形を集団平均することと、
    該第1の集団平均算出器の出力に基づいて、脈拍数を算出することと、
    正規化された波形を生成するために該波形を正規化することと、
    第2の集団平均算出器において、該正規化された波形を集団平均することと、
    該第2の集団平均算出器の出力に基づいて、酸素飽和を算出することと
    を包含する、方法。
  2. 前記集団平均することが可変の重みを使用することと、
    前記脈拍数を算出することを最適化するために、前記第1の集団平均算出器のための第1の測定基準を選択することと、
    前記酸素飽和を算出することを最適化するために、前記第2の集団平均算出器のための第2の測定基準を選択することと
    をさらに包含する、請求項1に記載の方法。
  3. 前記第1および第2の測定基準の両方が、不整脈パルスを検出するための不整脈測定基準を含み、該第1の測定基準のための該不整脈測定基準が、脈拍数を算出することと関連して、不整脈を認識するために、該第2の測定基準のための該不整脈測定基準より低い関連するしきい値を有する、請求項2に記載の方法。
  4. 前記第1および第2の測定基準の両方が、パルス振幅における短期変化の計測である短期測定基準を含み、
    パルス振幅における短期減少に応答して、前記第1の集団平均算出器が、前記第2の集団平均算出器より早く集団平均重みを増やす、請求項2に記載の方法。
  5. 酸素飽和および脈拍数を決定するためのパルス酸素濃度計であって、該パルス酸素濃度計は、
    患者からの2つの異なる光の波長に対応する波形を受け取る検出器と、
    第1の集団平均算出器と、
    該第1の集団平均算出器の出力に結合された脈拍数算出器と、
    正規化された波形を生成するために該波形を正規化する、該検出器に結合された正規化器と、
    第2の集団平均算出器と、
    該第2の集団平均算出器の出力に結合された酸素飽和算出器と
    を備える、パルス酸素濃度計。
  6. 前記集団平均算出器が、可変の重みを使用して集団平均するように構成され、
    前記脈拍数を算出することを最適化するために前記第1の集団平均算出器のための第1の測定基準を提供するように、および前記酸素飽和を算出することを最適化するために前記第2の集団平均算出器のための第2の測定基準を提供するように構成される、信号品質測定基準算出器をさらに備える、請求項5に記載のパルス酸素濃度計。
  7. 酸素飽和および脈拍数を決定するためにパルス酸素濃度計における信号を処理する方法であって、該方法は、
    患者からの2つの異なる光の波長に対応する波形を受け取ることと、
    第1の低域通過フィルタにおいて、該波形を低域通過フィルタリングすることと、
    該第1の低域通過フィルタの出力に基づいて、脈拍数を算出することと、
    正規化された波形を生成するために該波形を正規化することと、
    第2の低域通過フィルタにおいて、該正規化された波形を低域通過フィルタリングすることと、
    該第2の低域通過フィルタの出力に基づいて、酸素飽和を算出することと
    を包含する、方法。
  8. 前記脈拍数を算出することを最適化するために、前記第1の低域通過フィルタのための第1の測定基準を選択することと、
    前記酸素飽和を算出することを最適化するために、前記第2の低域通過フィルタのための第2の測定基準を選択することと
    をさらに包含する、請求項7に記載の方法。
  9. 前記第1の低域通過フィルタに関連する前記低域通過フィルタリング重みが、脈拍数推定に関する前記波形の周波数コンテンツを量化する周波数比測定基準に基づく、請求項8に記載の方法。
  10. 前記第2の低域通過フィルタのための低域通過フィルタリング重みが、その測定基準の脈拍数推定に関する前記波形の周波数コンテンツを量化する周波数比測定基準および別個の比率の比率分散測定基準に基づく、請求項8に記載の方法。
  11. 酸素飽和および脈拍数を決定するためにパルス酸素濃度計における信号を処理する方法であって、該方法は、
    患者からの2つの異なる光の波長に対応する波形を受け取ることと、
    第1の低域通過フィルタおよび集団平均算出器において、該波形を低域通過フィルタリングおよび集団平均することと、
    該第1の低域通過フィルタおよび集団平均算出器の出力に基づいて、脈拍数を算出することと、
    正規化された波形を生成するために該波形を正規化することと、
    第2の低域通過フィルタおよび集団平均算出器において、該正規化された波形を低域通過フィルタリングおよび集団平均することと、
    該第2の低域通過フィルタおよび集団平均算出器の出力に基づいて、酸素飽和を算出することと
    を包含する、方法。
  12. 酸素飽和および脈拍数を決定するためのパルス酸素濃度計であって、該パルス酸素濃度計は、
    患者からの2つの異なる光の波長に対応する波形を受け取る検出器と、
    第1の低域通過フィルタと、
    該第1の低域通過フィルタの出力に結合された脈拍数算出器と、
    正規化された波形を生成するために該波形を正規化する、該検出器に結合された正規化器と、
    第2の低域通過フィルタと、
    該第2の低域通過フィルタの出力に結合された酸素飽和算出器と
    を備える、パルス酸素濃度計。
  13. 前記低域通過フィルタが、可変の重みを使用して集団平均するように構成され、
    前記脈拍数を算出することを最適化するために前記第1の低域通過フィルタのための第1の測定基準を提供するように、および前記酸素飽和を算出することを最適化するために前記第2の低域通過フィルタのための第2の測定基準を提供するように構成される、信号品質測定基準算出器をさらに備える、請求項12に記載のパルス酸素濃度計。
  14. 前記第1の低域通過フィルタに関連する前記低域通過フィルタリング重みが、脈拍数推定に関する前記波形の周波数コンテンツを量化する周波数比測定基準に基づく、請求項12に記載のパルス酸素濃度計。
  15. 前記第2の低域通過フィルタのための低域通過フィルタリング重みが、その測定基準の脈拍数推定に関する前記波形の周波数コンテンツを量化する周波数比測定基準および別個の比率の比率分散測定基準に基づく、請求項12に記載のパルス酸素濃度計。
  16. 酸素飽和および脈拍数を決定するためのパルス酸素濃度計であって、該パルス酸素濃度計は、
    患者からの2つの異なる光の波長に対応する波形を受け取る検出器と、
    第1の低域通過フィルタおよび集団平均算出器と、
    該第1の低域通過フィルタおよび集団平均算出器の出力に結合された脈拍数算出器と、
    正規化された波形を生成するために該波形を正規化する、該検出器に結合された正規化器と、
    第2の低域通過フィルタおよび集団平均算出器と、
    該第2の低域通過フィルタおよび集団平均算出器の出力に結合された酸素飽和算出器と
    を備える、パルス酸素濃度計。
  17. 酸素飽和を決定するためにパルス酸素濃度計における信号を処理する方法であって、該方法は、
    患者からの2つの異なる光の波長に対応する波形を受け取ることと、
    履歴の平均波形と集団平均するために、パルス周期トリガの後に新しい波形を処理することと、
    該新しい波形が、所定のしきい値より多く該履歴の平均波形と異なる場合、新しい複合の履歴の平均波形の最初の数サンプルに対して、該新しい波形と該履歴の平均波形との間に補間することと
    を包含する、方法。
  18. 前記最初の数サンプルが4つのサンプルであり、前記補間が、前記新しい波形と前記履歴の平均波形との間の80%、60%、40%、および20%の差における、請求項17に記載の方法。
JP2007502899A 2004-03-08 2005-03-07 酸素飽和および心拍数のための個別の集団平均を用いるパルス酸素濃度計 Active JP4700682B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US10/796,578 2004-03-08
US10/796,578 US7039538B2 (en) 2004-03-08 2004-03-08 Pulse oximeter with separate ensemble averaging for oxygen saturation and heart rate
PCT/US2005/007388 WO2005087094A1 (en) 2004-03-08 2005-03-07 Pulse oximeter with separate ensemble averaging for oxygen saturation and heart rate

Publications (2)

Publication Number Publication Date
JP2007527771A true JP2007527771A (ja) 2007-10-04
JP4700682B2 JP4700682B2 (ja) 2011-06-15

Family

ID=34912588

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007502899A Active JP4700682B2 (ja) 2004-03-08 2005-03-07 酸素飽和および心拍数のための個別の集団平均を用いるパルス酸素濃度計

Country Status (13)

Country Link
US (2) US7039538B2 (ja)
EP (1) EP1737337B1 (ja)
JP (1) JP4700682B2 (ja)
KR (1) KR20070004803A (ja)
CN (1) CN1933773A (ja)
AT (1) ATE466523T1 (ja)
AU (1) AU2005221651A1 (ja)
CA (1) CA2557950C (ja)
DE (1) DE602005021086D1 (ja)
ES (1) ES2342945T3 (ja)
MX (1) MXPA06010305A (ja)
NZ (1) NZ549512A (ja)
WO (1) WO2005087094A1 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012071018A (ja) * 2010-09-29 2012-04-12 Denso Corp 脈波解析装置および血圧推定装置
JP2016083367A (ja) * 2014-10-28 2016-05-19 サルーステック株式会社 検体情報処理装置、情報処理方法、情報処理プログラム、及び同プログラムを記録したコンピュータ読み取り可能な記録媒体
JP2016539718A (ja) * 2013-12-12 2016-12-22 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 対象のバイタルサインを決定する装置及び方法
JP2021053464A (ja) * 2014-01-27 2021-04-08 リズム ダイアグノスティック システムズ,インク. 生理パラメータをモニタリングするデバイス

Families Citing this family (92)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL155955A0 (en) * 2003-05-15 2003-12-23 Widemed Ltd Adaptive prediction of changes of physiological/pathological states using processing of biomedical signal
US7039538B2 (en) * 2004-03-08 2006-05-02 Nellcor Puritant Bennett Incorporated Pulse oximeter with separate ensemble averaging for oxygen saturation and heart rate
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
US7277741B2 (en) * 2004-03-09 2007-10-02 Nellcor Puritan Bennett Incorporated Pulse oximetry motion artifact rejection using near infrared absorption by water
WO2006082589A2 (en) * 2005-02-07 2006-08-10 Widemed Ltd. Detection and monitoring of stress events during sleep
US7578793B2 (en) * 2004-11-22 2009-08-25 Widemed Ltd. Sleep staging based on cardio-respiratory signals
US7697967B2 (en) 2005-12-28 2010-04-13 Abbott Diabetes Care Inc. Method and apparatus for providing analyte sensor insertion
WO2006094168A1 (en) 2005-03-01 2006-09-08 Masimo Laboratories, Inc. Noninvasive multi-parameter patient monitor
US7392075B2 (en) 2005-03-03 2008-06-24 Nellcor Puritan Bennett Incorporated Method for enhancing pulse oximetry calculations in the presence of correlated artifacts
US20070100220A1 (en) 2005-10-28 2007-05-03 Baker Clark R Jr Adjusting parameters used in pulse oximetry analysis
US11298058B2 (en) 2005-12-28 2022-04-12 Abbott Diabetes Care Inc. Method and apparatus for providing analyte sensor insertion
US7706852B2 (en) 2006-01-30 2010-04-27 Nellcor Puritan Bennett Llc System and method for detection of unstable oxygen saturation
US8152731B2 (en) * 2006-02-10 2012-04-10 Inovise Medical, Inc. Wavelet transform and pattern recognition method for heart sound analysis
US7885698B2 (en) 2006-02-28 2011-02-08 Abbott Diabetes Care Inc. Method and system for providing continuous calibration of implantable analyte sensors
US7653425B2 (en) 2006-08-09 2010-01-26 Abbott Diabetes Care Inc. Method and system for providing calibration of an analyte sensor in an analyte monitoring system
US8473022B2 (en) 2008-01-31 2013-06-25 Abbott Diabetes Care Inc. Analyte sensor with time lag compensation
US9339217B2 (en) 2011-11-25 2016-05-17 Abbott Diabetes Care Inc. Analyte monitoring system and methods of use
US8346335B2 (en) 2008-03-28 2013-01-01 Abbott Diabetes Care Inc. Analyte sensor calibration management
US8374668B1 (en) 2007-10-23 2013-02-12 Abbott Diabetes Care Inc. Analyte sensor with lag compensation
US8140312B2 (en) 2007-05-14 2012-03-20 Abbott Diabetes Care Inc. Method and system for determining analyte levels
US7618369B2 (en) 2006-10-02 2009-11-17 Abbott Diabetes Care Inc. Method and system for dynamically updating calibration parameters for an analyte sensor
US8123695B2 (en) * 2006-09-27 2012-02-28 Nellcor Puritan Bennett Llc Method and apparatus for detection of venous pulsation
US7680522B2 (en) 2006-09-29 2010-03-16 Nellcor Puritan Bennett Llc Method and apparatus for detecting misapplied sensors
US8280469B2 (en) 2007-03-09 2012-10-02 Nellcor Puritan Bennett Llc Method for detection of aberrant tissue spectra
US20080220512A1 (en) * 2007-03-09 2008-09-11 Nellcor Puritan Bennett Llc Tunable laser-based spectroscopy system for non-invasively measuring body water content
US8229530B2 (en) * 2007-03-09 2012-07-24 Nellcor Puritan Bennett Llc System and method for detection of venous pulsation
US20080221426A1 (en) * 2007-03-09 2008-09-11 Nellcor Puritan Bennett Llc Methods and apparatus for detecting misapplied optical sensors
US8109882B2 (en) * 2007-03-09 2012-02-07 Nellcor Puritan Bennett Llc System and method for venous pulsation detection using near infrared wavelengths
US8221326B2 (en) * 2007-03-09 2012-07-17 Nellcor Puritan Bennett Llc Detection of oximetry sensor sites based on waveform characteristics
US8265724B2 (en) 2007-03-09 2012-09-11 Nellcor Puritan Bennett Llc Cancellation of light shunting
WO2008118993A1 (en) 2007-03-27 2008-10-02 Masimo Laboratories, Inc. Multiple wavelength optical sensor
EP2146625B1 (en) 2007-04-14 2019-08-14 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in medical communication system
US20080262326A1 (en) * 2007-04-19 2008-10-23 Starr Life Sciences Corp. Signal Processing Method and Apparatus for Processing a Physiologic Signal such as a Photoplethysmography Signal
US8374665B2 (en) 2007-04-21 2013-02-12 Cercacor Laboratories, Inc. Tissue profile wellness monitor
US8600681B2 (en) 2007-05-14 2013-12-03 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in a medical communication system
US8444560B2 (en) 2007-05-14 2013-05-21 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in a medical communication system
US8260558B2 (en) * 2007-05-14 2012-09-04 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in a medical communication system
US8103471B2 (en) 2007-05-14 2012-01-24 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in a medical communication system
US8560038B2 (en) 2007-05-14 2013-10-15 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in a medical communication system
US9125548B2 (en) 2007-05-14 2015-09-08 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in a medical communication system
US8239166B2 (en) 2007-05-14 2012-08-07 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in a medical communication system
US20080300500A1 (en) * 2007-05-30 2008-12-04 Widemed Ltd. Apnea detection using a capnograph
US8160900B2 (en) 2007-06-29 2012-04-17 Abbott Diabetes Care Inc. Analyte monitoring and management device and method to analyze the frequency of user interaction with the device
US8409093B2 (en) 2007-10-23 2013-04-02 Abbott Diabetes Care Inc. Assessing measures of glycemic variability
US8105241B2 (en) * 2007-10-26 2012-01-31 Inovise Medical, Inc. Combining to a singularity plural-anatomical-site, heart-functionality measurements
US8346328B2 (en) * 2007-12-21 2013-01-01 Covidien Lp Medical sensor and technique for using the same
US8352004B2 (en) * 2007-12-21 2013-01-08 Covidien Lp Medical sensor and technique for using the same
US20090171226A1 (en) * 2007-12-31 2009-07-02 Nellcor Puritan Bennett Llc System and method for evaluating variation in the timing of physiological events
US8185190B2 (en) * 2008-01-29 2012-05-22 Inovise Medical, Inc. Assessment of ischemia, and risk of sudden cardiac death, via heart-functionality parameter and acoustic cardiographic monitoring
US8348852B2 (en) * 2008-03-06 2013-01-08 Inovise Medical, Inc. Heart-activity sound monitoring
US8924159B2 (en) 2008-05-30 2014-12-30 Abbott Diabetes Care Inc. Method and apparatus for providing glycemic control
US8591410B2 (en) 2008-05-30 2013-11-26 Abbott Diabetes Care Inc. Method and apparatus for providing glycemic control
US8862194B2 (en) 2008-06-30 2014-10-14 Covidien Lp Method for improved oxygen saturation estimation in the presence of noise
US8986208B2 (en) 2008-09-30 2015-03-24 Abbott Diabetes Care Inc. Analyte sensor sensitivity attenuation mitigation
US8577448B2 (en) * 2008-10-14 2013-11-05 Inovise Medical, Inc. Differential apneic detection in aid of diagnosis and treatment
US20090171172A1 (en) * 2008-12-19 2009-07-02 Nellcor Puritan Bennett Llc Method and system for pulse gating
WO2010121084A1 (en) 2009-04-15 2010-10-21 Abbott Diabetes Care Inc. Analyte monitoring system having an alert
US8311601B2 (en) * 2009-06-30 2012-11-13 Nellcor Puritan Bennett Llc Reflectance and/or transmissive pulse oximeter
EP3689237B1 (en) 2009-07-23 2021-05-19 Abbott Diabetes Care, Inc. Method of manufacturing and system for continuous analyte measurement
DK4070729T3 (da) 2009-08-31 2024-05-21 Abbott Diabetes Care Inc Visningsindretninger til en medicinsk indretning
WO2011041469A1 (en) 2009-09-29 2011-04-07 Abbott Diabetes Care Inc. Method and apparatus for providing notification function in analyte monitoring systems
US8409108B2 (en) * 2009-11-05 2013-04-02 Inovise Medical, Inc. Multi-axial heart sounds and murmur detection for hemodynamic-condition assessment
US9839381B1 (en) 2009-11-24 2017-12-12 Cercacor Laboratories, Inc. Physiological measurement system with automatic wavelength adjustment
DE112010004682T5 (de) 2009-12-04 2013-03-28 Masimo Corporation Kalibrierung für mehrstufige physiologische Monitore
US7884933B1 (en) 2010-05-05 2011-02-08 Revolutionary Business Concepts, Inc. Apparatus and method for determining analyte concentrations
WO2012036490A2 (ko) * 2010-09-16 2012-03-22 케이티메드 주식회사 생체신호의 디지털 신호 처리용 실험키트
US9211095B1 (en) * 2010-10-13 2015-12-15 Masimo Corporation Physiological measurement logic engine
US9072433B2 (en) 2011-02-18 2015-07-07 Covidien Lp Method and apparatus for noninvasive blood pressure measurement using pulse oximetry
US8721557B2 (en) 2011-02-18 2014-05-13 Covidien Lp Pattern of cuff inflation and deflation for non-invasive blood pressure measurement
CN102178536B (zh) * 2011-03-29 2013-04-03 苏州易寻传感网络科技有限公司 一种血氧饱和度和心率的测量方法
US8710993B2 (en) 2011-11-23 2014-04-29 Abbott Diabetes Care Inc. Mitigating single point failure of devices in an analyte monitoring system and methods thereof
US9370308B2 (en) 2012-05-01 2016-06-21 Nellcor Puritan Bennett Ireland Angle distribution technique for analyzing a physiological sensor signal
ES2397768B1 (es) 2012-06-18 2014-01-21 Tecatel, S.L. Terminal multimedia con medida de parámetros vitales
WO2014035732A1 (en) 2012-08-30 2014-03-06 Abbot Diabetes Care Inc. Dropout detection in continuous analyte monitoring data during data excursions
US9241670B2 (en) 2012-09-11 2016-01-26 Covidien Lp Methods and systems for conditioning physiological information using a normalization technique
US8548588B1 (en) 2012-09-21 2013-10-01 Inovise Medical, Inc. CRM-device ventricular-pacing blanking control
EP2901153A4 (en) 2012-09-26 2016-04-27 Abbott Diabetes Care Inc METHOD AND DEVICE FOR IMPROVING DELAY CORRECTION FUNCTION DURING IN VIVO MEASUREMENT OF ANALYZ CONCENTRATION WITH ANALYZ CONCENTRATION VARIABILITY AND RANGE DATA
US10413251B2 (en) 2012-10-07 2019-09-17 Rhythm Diagnostic Systems, Inc. Wearable cardiac monitor
EP2934668B1 (en) 2012-12-19 2018-08-22 VisCardia, Inc. Hemodynamic performance enhancement through asymptomatic diaphragm stimulation
US10335592B2 (en) 2012-12-19 2019-07-02 Viscardia, Inc. Systems, devices, and methods for improving hemodynamic performance through asymptomatic diaphragm stimulation
US20210383403A1 (en) * 2014-01-15 2021-12-09 Federal Law Enforcement Development Services, Inc. UV, SOUND POINT, iA OPERATING SYSTEM
DE102014204665B4 (de) * 2014-03-13 2019-01-24 Siemens Healthcare Gmbh Geräuschoptimierung einer Magnetresonanzanlage
KR101560269B1 (ko) 2015-03-09 2015-10-14 주식회사 휴이노 심전도 신호를 처리하기 위한 방법, 시스템 및 비일시성의 컴퓨터 판독 가능한 기록 매체
US10342466B2 (en) 2015-03-24 2019-07-09 Covidien Lp Regional saturation system with ensemble averaging
WO2017011346A1 (en) 2015-07-10 2017-01-19 Abbott Diabetes Care Inc. System, device and method of dynamic glucose profile response to physiological parameters
US10493271B2 (en) 2016-04-29 2019-12-03 Viscardia, Inc. Implantable medical devices, systems, and methods for selection of optimal diaphragmatic stimulation parameters to affect pressures within the intrathoracic cavity
WO2018175489A1 (en) 2017-03-21 2018-09-27 Abbott Diabetes Care Inc. Methods, devices and system for providing diabetic condition diagnosis and therapy
US11903700B2 (en) 2019-08-28 2024-02-20 Rds Vital signs monitoring systems and methods
US11925803B2 (en) 2019-09-26 2024-03-12 Viscardia, Inc. Implantable medical systems, devices, and methods for affecting cardiac function through diaphragm stimulation, and for monitoring diaphragmatic health
EP4126191A1 (en) 2020-03-27 2023-02-08 VisCardia, Inc. Implantable medical systems, devices and methods for delivering asymptomatic diaphragmatic stimulation
WO2021219779A1 (en) * 2020-04-29 2021-11-04 Rds Vital signs or health monitoring systems and methods
CN111529829B (zh) * 2020-07-07 2020-11-27 深圳市汇顶科技股份有限公司 一种ppg设备的信号调整方法和ppg设备

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0245041A (ja) * 1988-03-30 1990-02-15 Nellcor Inc 光パルス検出方法および装置
US4960126A (en) * 1988-01-15 1990-10-02 Criticare Systems, Inc. ECG synchronized pulse oximeter
JPH1052408A (ja) * 1996-04-15 1998-02-24 Ohmeda Inc 患者の灌流指数を生成する方法及び装置
JPH10507118A (ja) * 1991-03-07 1998-07-14 マシモ コーポレイション 信号処理装置
JPH10510440A (ja) * 1993-10-08 1998-10-13 ネルカー ピューリタン ベネット インコーポレイテッド 仮想トリガを使用する心拍同期化用パルス酸素計

Family Cites Families (16)

* 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
US4700708A (en) * 1982-09-02 1987-10-20 Nellcor Incorporated Calibrated optical oximeter probe
US4928692A (en) * 1985-04-01 1990-05-29 Goodman David E Method and apparatus for detecting optical pulses
US4802486A (en) 1985-04-01 1989-02-07 Nellcor Incorporated Method and apparatus for detecting optical pulses
US4690126A (en) * 1986-07-14 1987-09-01 Orley's Manufacturing Co., Inc. Catalytic combustion assembly for wood-burning stove
US5078136A (en) 1988-03-30 1992-01-07 Nellcor Incorporated Method and apparatus for calculating arterial oxygen saturation based plethysmographs including transients
US5853364A (en) 1995-08-07 1998-12-29 Nellcor Puritan Bennett, Inc. Method and apparatus for estimating physiological parameters using model-based adaptive filtering
US6035223A (en) * 1997-11-19 2000-03-07 Nellcor Puritan Bennett Inc. Method and apparatus for determining the state of an oximetry sensor
US6293915B1 (en) * 1997-11-20 2001-09-25 Seiko Epson Corporation Pulse wave examination apparatus, blood pressure monitor, pulse waveform monitor, and pharmacological action monitor
US6760609B2 (en) * 1999-07-14 2004-07-06 Providence Health System - Oregon Adaptive calibration pulsed oximetry method and device
US6339715B1 (en) * 1999-09-30 2002-01-15 Ob Scientific Method and apparatus for processing a physiological signal
IL138884A (en) * 2000-10-05 2006-07-05 Conmed Corp Pulse oximeter and a method of its operation
US7052466B2 (en) * 2001-04-11 2006-05-30 Cardiac Pacemakers, Inc. Apparatus and method for outputting heart sounds
US6711425B1 (en) * 2002-05-28 2004-03-23 Ob Scientific, Inc. Pulse oximeter with calibration stabilization
US6879850B2 (en) * 2002-08-16 2005-04-12 Optical Sensors Incorporated Pulse oximeter with motion detection
US7039538B2 (en) * 2004-03-08 2006-05-02 Nellcor Puritant Bennett Incorporated Pulse oximeter with separate ensemble averaging for oxygen saturation and heart rate

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4960126A (en) * 1988-01-15 1990-10-02 Criticare Systems, Inc. ECG synchronized pulse oximeter
JPH0245041A (ja) * 1988-03-30 1990-02-15 Nellcor Inc 光パルス検出方法および装置
JPH10507118A (ja) * 1991-03-07 1998-07-14 マシモ コーポレイション 信号処理装置
JPH10510440A (ja) * 1993-10-08 1998-10-13 ネルカー ピューリタン ベネット インコーポレイテッド 仮想トリガを使用する心拍同期化用パルス酸素計
JPH1052408A (ja) * 1996-04-15 1998-02-24 Ohmeda Inc 患者の灌流指数を生成する方法及び装置

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012071018A (ja) * 2010-09-29 2012-04-12 Denso Corp 脈波解析装置および血圧推定装置
JP2016539718A (ja) * 2013-12-12 2016-12-22 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 対象のバイタルサインを決定する装置及び方法
JP2021053464A (ja) * 2014-01-27 2021-04-08 リズム ダイアグノスティック システムズ,インク. 生理パラメータをモニタリングするデバイス
JP7117027B2 (ja) 2014-01-27 2022-08-12 アールディーエス エスエーエス 生理パラメータをモニタリングするデバイス
JP2016083367A (ja) * 2014-10-28 2016-05-19 サルーステック株式会社 検体情報処理装置、情報処理方法、情報処理プログラム、及び同プログラムを記録したコンピュータ読み取り可能な記録媒体

Also Published As

Publication number Publication date
KR20070004803A (ko) 2007-01-09
EP1737337B1 (en) 2010-05-05
AU2005221651A1 (en) 2005-09-22
CA2557950A1 (en) 2005-09-22
ATE466523T1 (de) 2010-05-15
US7209774B2 (en) 2007-04-24
EP1737337A1 (en) 2007-01-03
US20050197793A1 (en) 2005-09-08
JP4700682B2 (ja) 2011-06-15
DE602005021086D1 (de) 2010-06-17
CN1933773A (zh) 2007-03-21
US20060195280A1 (en) 2006-08-31
ES2342945T3 (es) 2010-07-19
CA2557950C (en) 2013-02-26
MXPA06010305A (es) 2007-04-17
WO2005087094A1 (en) 2005-09-22
NZ549512A (en) 2008-10-31
US7039538B2 (en) 2006-05-02

Similar Documents

Publication Publication Date Title
JP4700682B2 (ja) 酸素飽和および心拍数のための個別の集団平均を用いるパルス酸素濃度計
US7534212B2 (en) Pulse oximeter with alternate heart-rate determination
EP2260758B1 (en) Selection of ensemble averaging weights for a pulse oximeter based on signal quality metrics
US9351674B2 (en) Method for enhancing pulse oximetry calculations in the presence of correlated artifacts
US5853364A (en) Method and apparatus for estimating physiological parameters using model-based adaptive filtering
EP1006863A1 (en) Method and apparatus for adaptively averaging data signals

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20080229

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20080717

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20080717

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20101026

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110125

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20110215

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20110304

R150 Certificate of patent or registration of utility model

Ref document number: 4700682

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250