JP2009072639A - パルスフォトメータ - Google Patents
パルスフォトメータ Download PDFInfo
- Publication number
- JP2009072639A JP2009072639A JP2009006876A JP2009006876A JP2009072639A JP 2009072639 A JP2009072639 A JP 2009072639A JP 2009006876 A JP2009006876 A JP 2009006876A JP 2009006876 A JP2009006876 A JP 2009006876A JP 2009072639 A JP2009072639 A JP 2009072639A
- Authority
- JP
- Japan
- Prior art keywords
- pulse wave
- light
- norm
- pulse
- wave data
- 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
Links
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
【課題】 同一の媒体からほぼ同時に抽出される２つの同種の信号を処
理して共通の信号成分を抽出する計算処理負担を軽減したパルスフォトメータを提供する。
【解決手段】異なる２つの波長の光を同一の生体組織に照射する発光手段と、前記発光手段から発生し前記同一の生体組織を透過または反射した各波長の光を電気信号に変換する受光手段と、前記受光手段で変換された前記電気信号より得られた２つの波長の離散的時系列脈波データをそれぞれ縦軸または横軸とする２次元直交座標に展開した、各波長の脈波データそれぞれについてノルム値を求め、さらにそのノルム値の比を求めるノルム比計算手段と、前記ノルム比計算手段により計算されたノルム比に基づいて、血中の酸素の濃度を求める血中酸素濃度演算手段とを具備することを特徴とするパルスフォトメータ。
【選択図】 図７
理して共通の信号成分を抽出する計算処理負担を軽減したパルスフォトメータを提供する。
【解決手段】異なる２つの波長の光を同一の生体組織に照射する発光手段と、前記発光手段から発生し前記同一の生体組織を透過または反射した各波長の光を電気信号に変換する受光手段と、前記受光手段で変換された前記電気信号より得られた２つの波長の離散的時系列脈波データをそれぞれ縦軸または横軸とする２次元直交座標に展開した、各波長の脈波データそれぞれについてノルム値を求め、さらにそのノルム値の比を求めるノルム比計算手段と、前記ノルム比計算手段により計算されたノルム比に基づいて、血中の酸素の濃度を求める血中酸素濃度演算手段とを具備することを特徴とするパルスフォトメータ。
【選択図】 図７
Description
本発明は、一つの媒体からほぼ同時に抽出される２つの同種の信号を処理して共通の信
号成分を抽出する信号処理に関し、特には医療の分野において、特に循環器系の診断に用
いられるパルスフォトメータにおける信号処理の改良に関する。
号成分を抽出する信号処理に関し、特には医療の分野において、特に循環器系の診断に用
いられるパルスフォトメータにおける信号処理の改良に関する。
一つの媒体からほぼ同時に抽出された２つの信号から信号成分と雑音成分に分離する方
法には様々な方法が提案されている。
それらは、一般的には周波数領域や時間領域による処理が行われている。
医療現場でも、光電脈波計と言われる脈波波形や脈拍数を測定する装置、血液に含まれ
る吸光物質の濃度測定として、酸素飽和度ＳｐＯ2の測定装置、一酸化炭素ヘモグロビン
やＭｅｔヘモグロビン等の特殊ヘモグロビンの濃度の測定装置、注入色素濃度の測定装置
などがパルスフォトメータとして知られている。
中でも酸素飽和度ＳｐＯ2の測定装置を特にパルスオキシメータと呼んでいる。
法には様々な方法が提案されている。
それらは、一般的には周波数領域や時間領域による処理が行われている。
医療現場でも、光電脈波計と言われる脈波波形や脈拍数を測定する装置、血液に含まれ
る吸光物質の濃度測定として、酸素飽和度ＳｐＯ2の測定装置、一酸化炭素ヘモグロビン
やＭｅｔヘモグロビン等の特殊ヘモグロビンの濃度の測定装置、注入色素濃度の測定装置
などがパルスフォトメータとして知られている。
中でも酸素飽和度ＳｐＯ2の測定装置を特にパルスオキシメータと呼んでいる。
パルスフォトメータの原理は、対象物質への吸光性が異なる複数の波長の光を生体組織
に透過又は反射させ、その透過光又は反射光の光量を連続的に測定することで得られる脈
波データ信号から対象物質の濃度を求めるものである。
そしてその脈波データに雑音が混入すると、正しい濃度の計算が出来ず、誤処置につな
がる危険が生じる。
パルスフォトメータにおいても従来より雑音を低減するために周波数帯域を分割して信
号成分に着目したり、２つの信号の相関を取るなどの方法が提案されてきた。
しかし、これらの方法は解析に時間がかかるなどの問題があった。
に透過又は反射させ、その透過光又は反射光の光量を連続的に測定することで得られる脈
波データ信号から対象物質の濃度を求めるものである。
そしてその脈波データに雑音が混入すると、正しい濃度の計算が出来ず、誤処置につな
がる危険が生じる。
パルスフォトメータにおいても従来より雑音を低減するために周波数帯域を分割して信
号成分に着目したり、２つの信号の相関を取るなどの方法が提案されてきた。
しかし、これらの方法は解析に時間がかかるなどの問題があった。
そこで、本出願人は、特許第３２７０９１７号（特許文献１）において、異なる２つの
波長の光を生体組織に照射して透過光から得られる２つの脈波信号のそれぞれの大きさを
縦軸、横軸としてグラフを描き、その回帰直線を求め、その回帰直線の傾きに基づいて、
動脈血中の酸素飽和度ないし吸光物質濃度を求めることを提案している。
この発明により、測定精度を高め、低消費電力化することができた。
しかし、各波長の脈波信号についての多くのサンプリングデータを用いて回帰直線ない
しその傾きを求めるためには、なお多くの計算処理を要していた。
波長の光を生体組織に照射して透過光から得られる２つの脈波信号のそれぞれの大きさを
縦軸、横軸としてグラフを描き、その回帰直線を求め、その回帰直線の傾きに基づいて、
動脈血中の酸素飽和度ないし吸光物質濃度を求めることを提案している。
この発明により、測定精度を高め、低消費電力化することができた。
しかし、各波長の脈波信号についての多くのサンプリングデータを用いて回帰直線ない
しその傾きを求めるためには、なお多くの計算処理を要していた。
更に本出願人は、特願２００１−３３２３８３号（特許文献２）においては、周波数解
析を用いてはいるが、その解析においては従来技術のように脈波信号そのものを抽出する
のではなく、脈波信号の基本周波数を求め、さらには精度を高めるためにその高調波周波
数を用いたフィルタを用いて脈波信号をフィルタリングする方法を提案している。
しかし、基本周波数を求める点に関しては更なる改善が望まれていた。
析を用いてはいるが、その解析においては従来技術のように脈波信号そのものを抽出する
のではなく、脈波信号の基本周波数を求め、さらには精度を高めるためにその高調波周波
数を用いたフィルタを用いて脈波信号をフィルタリングする方法を提案している。
しかし、基本周波数を求める点に関しては更なる改善が望まれていた。
本発明の課題（目的）は、同一の媒体からほぼ同時に抽出される２つの同種の信号を処
理して共通の信号成分を抽出する計算処理負担を軽減したパルスフォトメータを提供することにある。
また、上記信号処理方法を適用して、前記媒体の体動によるノイズが脈波信号に生じた
場合であっても、対象物質の濃度を精度よく測定することにある。
また、体動によるノイズが脈波データ信号に生じた場合であっても、脈波信号からノイ
ズを除去し、精度よく脈拍を求めることにある。
理して共通の信号成分を抽出する計算処理負担を軽減したパルスフォトメータを提供することにある。
また、上記信号処理方法を適用して、前記媒体の体動によるノイズが脈波信号に生じた
場合であっても、対象物質の濃度を精度よく測定することにある。
また、体動によるノイズが脈波データ信号に生じた場合であっても、脈波信号からノイ
ズを除去し、精度よく脈拍を求めることにある。
・前記課題を解決するために、本発明のパルスフォトメータは、異なる２つの波長の光を同一の生体組織に照射する発光手段と、前記発光手段から発生し前記同一の生体組織を透過または反射した各波長の光を電気信号に変換する受光手段と、前記受光手段で変換された前記電気信号より得られた２つの波長の離散的時系列脈波データをそれぞれ縦軸または横軸とする２次元直交座標に展開した、各波長の脈波データそれぞれについてノルム値を求め、さらにそのノルム値の比を求めるノルム比計算手段と、前記ノルム比計算手段により計算されたノルム比に基づいて、血中の酸素の濃度を求める血中酸素濃度演算手段とを具備することを特徴とする。（請求項１）
また、参考例では、
・同一の媒体から抽出される所定期間の２つの信号を２次元直交座標上に展開するステップと、前記展開された信号を、回転行列を用いて、あらかじめ決められた角度または所定条件に基づいて決められた角度に回転させるステップとにより前記２つの信号を処理する。
この信号処理方法によって、同一の媒体からほぼ同時に抽出される２つの同種の信号を
処理して共通の信号成分を抽出する計算処理負担を軽減することができる。
・前記信号処理の結果、信号主成分の基本周波数を求める。
この構成により、ノイズが低減された基本波波形を容易に取得することができる。
・同一の媒体からの２つの信号を入力する入力手段と、前記入力手段から入力された所定
期間分の２つの信号を２次元直交座標上に展開する展開手段と、該展開されたデータをあ
らかじめ決められた角度または所定条件に基づいて決められた角度に回転させる回転行列
を用いて処理する処理手段とを具えた構成とする。
この構成により、回転行列を用いて処理することでノイズが低減された波形を容易に取
得することができる。
・同一の媒体から抽出される所定期間の２つの信号を２次元直交座標上に展開するステップと、前記展開された信号を、回転行列を用いて、あらかじめ決められた角度または所定条件に基づいて決められた角度に回転させるステップとにより前記２つの信号を処理する。
この信号処理方法によって、同一の媒体からほぼ同時に抽出される２つの同種の信号を
処理して共通の信号成分を抽出する計算処理負担を軽減することができる。
・前記信号処理の結果、信号主成分の基本周波数を求める。
この構成により、ノイズが低減された基本波波形を容易に取得することができる。
・同一の媒体からの２つの信号を入力する入力手段と、前記入力手段から入力された所定
期間分の２つの信号を２次元直交座標上に展開する展開手段と、該展開されたデータをあ
らかじめ決められた角度または所定条件に基づいて決められた角度に回転させる回転行列
を用いて処理する処理手段とを具えた構成とする。
この構成により、回転行列を用いて処理することでノイズが低減された波形を容易に取
得することができる。
・異なる２つの波長の光を生体組織に照射する発光手段と、前記発光手段から発生し前記
生体組織を透過または反射した各波長の光を電気信号に変換する受光手段とを具えたパル
スフォトメータにおいて、
前記各波長の電気信号より得られた脈波データを、あらかじめ決められた角度にまたは
所定条件に基づいて決められた角度に回転させる回転行列を用いて処理する処理手段とを
具備することを特徴とする。
この構成により、脈波データを回転行列を用いて処理することでノイズが低減された波
形を容易に取得することができる。
・前記脈波データが、ほぼ同時点に受光した２つの波長の受光信号から得られる電気信号
を、前記２つの波長をそれぞれ縦軸または横軸とする２次元直交座標に展開したものであ
ることを特徴とする。
・前記脈波データは所定時間分であり、かつ経時移動して処理されることを特徴とする。
・異なる２つの波長の光を生体組織に照射する発光手段と、前記発光手段から発生し前記
生体組織を透過または反射した各波長の光を電気信号に変換する受光手段とを具えたパル
スフォトメータにおいて、前記各波長の電気信号より得られた脈波信号を、各波長をそれ
ぞれ縦軸または横軸とする２次元直交座標に展開し、縦軸または横軸に射影される領域が
最大または最小となるいずれかの条件を満足するように回転行列を用いて回転させて処理
する第２の処理手段とを具備することを特徴とする。
生体組織を透過または反射した各波長の光を電気信号に変換する受光手段とを具えたパル
スフォトメータにおいて、
前記各波長の電気信号より得られた脈波データを、あらかじめ決められた角度にまたは
所定条件に基づいて決められた角度に回転させる回転行列を用いて処理する処理手段とを
具備することを特徴とする。
この構成により、脈波データを回転行列を用いて処理することでノイズが低減された波
形を容易に取得することができる。
・前記脈波データが、ほぼ同時点に受光した２つの波長の受光信号から得られる電気信号
を、前記２つの波長をそれぞれ縦軸または横軸とする２次元直交座標に展開したものであ
ることを特徴とする。
・前記脈波データは所定時間分であり、かつ経時移動して処理されることを特徴とする。
・異なる２つの波長の光を生体組織に照射する発光手段と、前記発光手段から発生し前記
生体組織を透過または反射した各波長の光を電気信号に変換する受光手段とを具えたパル
スフォトメータにおいて、前記各波長の電気信号より得られた脈波信号を、各波長をそれ
ぞれ縦軸または横軸とする２次元直交座標に展開し、縦軸または横軸に射影される領域が
最大または最小となるいずれかの条件を満足するように回転行列を用いて回転させて処理
する第２の処理手段とを具備することを特徴とする。
・異なる２つの波長の光を生体組織に照射する発光手段と、前記発光手段で発生し前記生
体組織を透過または反射した各波長の光を電気信号に変換する受光手段と、前記受光部に
より得られる前記２つの波長の脈波信号をそれぞれ縦軸または横軸とする２次元直交座標
に展開した、各波長の各波長の脈波データを、所定角度に回転させる回転行列を用いて処
理した波形を求める波形取得手段と、前記波形取得手段により得られた波形の周波数解析
により、脈波の基本周波数または脈拍数を求める波形解析手段と、前記波形解析手段の出
力から血中の吸光物質の濃度を求める血中吸光物質濃度演算手段とを具備するパルスフォ
トメータ。
・前記血中吸光物質濃度演算手段は、動脈血中の酸素飽和度、特殊ヘモグロビン濃度、ま
たは注入色素濃度のうち少なくとも１つを演算する。
体組織を透過または反射した各波長の光を電気信号に変換する受光手段と、前記受光部に
より得られる前記２つの波長の脈波信号をそれぞれ縦軸または横軸とする２次元直交座標
に展開した、各波長の各波長の脈波データを、所定角度に回転させる回転行列を用いて処
理した波形を求める波形取得手段と、前記波形取得手段により得られた波形の周波数解析
により、脈波の基本周波数または脈拍数を求める波形解析手段と、前記波形解析手段の出
力から血中の吸光物質の濃度を求める血中吸光物質濃度演算手段とを具備するパルスフォ
トメータ。
・前記血中吸光物質濃度演算手段は、動脈血中の酸素飽和度、特殊ヘモグロビン濃度、ま
たは注入色素濃度のうち少なくとも１つを演算する。
請求項１に記載のパルスフォトメータによれば、各波長の脈波データから得られるノル
ム比から精度よく血中の酸素の濃度を求めることができる。
また、パルスフォトメータによれば、動脈血中の酸素飽和度、特殊ヘモグロビン濃度、または注入色素濃度のうち少なくとも１つを演算することができる。
ム比から精度よく血中の酸素の濃度を求めることができる。
また、パルスフォトメータによれば、動脈血中の酸素飽和度、特殊ヘモグロビン濃度、または注入色素濃度のうち少なくとも１つを演算することができる。
参考例のパルスフォトメータによれば、脈波データを回転行列を用いて処理することでノイズが低減された波形を取得することができる。
参考例のパルスフォトメータによれば、脈波データを回転行列を用いて処理することでノイズが低減された波形を取得し、その波形から精度よく脈波データの基本周波数または脈拍数を求めることができる。
参考例のパルスフォトメータによれば、脈波データを回転行列を用いてノイズを除去するのに適した回転角度を決定でき、その決定した回転角度の回転行列によって処理することでノイズが低減された波形を取得することができる。
参考例のパルスフォトメータによれば、脈波データを回転行列を用いてノイズを除去するのに適した回転角度を決定でき、その決定した回転角度の回転行列によって処理することでノイズが低減された波形を取得し、その波形から精度よく脈波データの基本周波数または脈拍数を求めることができる。
参考例のパルスフォトメータによれば、回転行列の回転角の決定を容易に決定できる。
参考例のパルスフォトメータによれば、脈波データを回転行列を用いて処理することでノイズが低減された波形による精度の良い脈波の基本周波数または脈拍数の取得と、動脈血中の酸素飽和度、特殊ヘモグロビン濃度、または注入色素濃度のうち少なくとも１つの演算とが同一の装置でほぼ同時に処理して出力することができる。
参考例のパルスフォトメータによれば、動脈血中の酸素飽和度、特殊ヘモグロビン濃度、または注入色素濃度のうち少なくとも１つを演算することができる。
参考例のパルスフォトメータによれば、動脈血中の酸素飽和度、特殊ヘモグロビン濃度、または注入色素濃度のうち少なくとも１つを演算することができる。
本発明の実施の形態を説明するにあたり、動脈血酸素飽和度を測定するパルスオキシメ
ータを例に挙げて原理を説明する。
なお、本発明の技術は、パルスオキシメータに限られず、特殊ヘモグロビン（一酸化炭
素ヘモグロビン、Ｍｅｔヘモグロビンなど）、血中に注入された色素などの血中吸光物質
をパルスフォトメトリーの原理を用いて測定する装置(パルスフォトメータ)に適用できる。
ータを例に挙げて原理を説明する。
なお、本発明の技術は、パルスオキシメータに限られず、特殊ヘモグロビン（一酸化炭
素ヘモグロビン、Ｍｅｔヘモグロビンなど）、血中に注入された色素などの血中吸光物質
をパルスフォトメトリーの原理を用いて測定する装置(パルスフォトメータ)に適用できる。
動脈血酸素飽和度を測定するパルスオキシメータの構成は、概略構成ブロック図である
図１のようになっている。
異なる波長の光を発光する発光素子１、２は、交互に発光するように駆動回路３により
駆動される。
発光素子１、２に採用する光はそれぞれ動脈血酸素飽和度による影響が少ない赤外光（
例えば９４０[nm]）、動脈血酸素飽和度の変化に対する感度が高い赤色光（例えば６６０
[nm]）がよい。
図１のようになっている。
異なる波長の光を発光する発光素子１、２は、交互に発光するように駆動回路３により
駆動される。
発光素子１、２に採用する光はそれぞれ動脈血酸素飽和度による影響が少ない赤外光（
例えば９４０[nm]）、動脈血酸素飽和度の変化に対する感度が高い赤色光（例えば６６０
[nm]）がよい。
これらの発光素子１、２からの発光は生体組織４を透過してフォトダイオード５で受光
して電気信号に変換される。
なお、反射光を受光するようにしてもよい。
そして、これらの変換された信号は増幅器６で増幅され、マルチプレクサ７によりそれ
ぞれの光波長に対応したフィルタ８−１、８−２に振り分けられる。
各フィルターに振り分けられた信号はフィルタ８−１、８−２によりフィルタリングさ
れてノイズ成分が低減され、Ａ／Ｄ変換器９によりデジタル化される。
して電気信号に変換される。
なお、反射光を受光するようにしてもよい。
そして、これらの変換された信号は増幅器６で増幅され、マルチプレクサ７によりそれ
ぞれの光波長に対応したフィルタ８−１、８−２に振り分けられる。
各フィルターに振り分けられた信号はフィルタ８−１、８−２によりフィルタリングさ
れてノイズ成分が低減され、Ａ／Ｄ変換器９によりデジタル化される。
デジタル化された赤外光、赤色光に対応する各信号列が、それぞれの脈波信号を形成し
ている。
デジタル化された各信号列は処理部１０に入力され、ＲＯＭ１２に格納されているプロ
グラムにより処理され、酸素飽和度ＳｐＯ2が測定され、その値が表示部１１に表示され
る。
ている。
デジタル化された各信号列は処理部１０に入力され、ＲＯＭ１２に格納されているプロ
グラムにより処理され、酸素飽和度ＳｐＯ2が測定され、その値が表示部１１に表示され
る。
＜回転行列によるノイズ低減と脈波の基本周波数の演算＞
先ず、血液中の吸光物質の吸光度（減光度）の変動の測定について説明する。
図９(a)及び(b)は、前記発光素子１、２からの発光された光が生体組織４を透過してフ
ォトダイオード５で受光して電気信号に変換された脈波データで、(a)は赤色光の場合を
、(b)は赤外光を示している。
図９の(a)では、横軸を時間、縦軸を受光出力とすると、フォトダイオード５での受光
出力は、赤色光の直流成分（Ｒ’）と脈動成分（ΔＲ’）が重畳された波形となっている。
また、図９の(b)では、横軸を時間、縦軸を受光出力とすると、フォトダイオード５で
の受光出力は、赤外光の直流成分（ＩＲ’）と脈動成分（ΔＩＲ’）が重畳された波形となっている。
図２は、図９に示すような脈波において、８秒間分の、直流成分（Ｒ’、ＩＲ’）に対
する脈動成分（ΔＲ’、ΔＩＲ’）の比（ＩＲ＝ΔＩＲ’／ＩＲ’）をとり、さらにその８秒間分のデータの平均値をゼロに合わせたものである。
なお、図２の如き、平均値をゼロとする処理を行わなくとも演算は可能である。
図３は、図９に示される赤外光ＩＲのデータを横軸に、赤色光Ｒのデータを縦軸にとっ
たグラフである。
先ず、血液中の吸光物質の吸光度（減光度）の変動の測定について説明する。
図９(a)及び(b)は、前記発光素子１、２からの発光された光が生体組織４を透過してフ
ォトダイオード５で受光して電気信号に変換された脈波データで、(a)は赤色光の場合を
、(b)は赤外光を示している。
図９の(a)では、横軸を時間、縦軸を受光出力とすると、フォトダイオード５での受光
出力は、赤色光の直流成分（Ｒ’）と脈動成分（ΔＲ’）が重畳された波形となっている。
また、図９の(b)では、横軸を時間、縦軸を受光出力とすると、フォトダイオード５で
の受光出力は、赤外光の直流成分（ＩＲ’）と脈動成分（ΔＩＲ’）が重畳された波形となっている。
図２は、図９に示すような脈波において、８秒間分の、直流成分（Ｒ’、ＩＲ’）に対
する脈動成分（ΔＲ’、ΔＩＲ’）の比（ＩＲ＝ΔＩＲ’／ＩＲ’）をとり、さらにその８秒間分のデータの平均値をゼロに合わせたものである。
なお、図２の如き、平均値をゼロとする処理を行わなくとも演算は可能である。
図３は、図９に示される赤外光ＩＲのデータを横軸に、赤色光Ｒのデータを縦軸にとっ
たグラフである。
次に、Ａ／Ｄ変換器９によってデジタル化した各波長の２つの脈波データ信号を回転行
列を用いてノイズを低減する演算処理について説明する。
なお、赤外光と赤色光とは交互に発光されるため厳密には同時に発光されるものではな
いが、隣り合う得られた赤外光受光値と赤色光受光値を同時刻に得られたものとして扱い
、所定時間分の赤外光の脈波信号と赤色光の脈波信号を２次元直交座標上に展開する。
すなわち図３のグラフを作成している。
また、脈波の直流成分に対する脈動成分の比をとることで脈拍による吸光度の脈動分が
近似される。
図３のグラフに見られる推移は４５度になっていないが、その理由は、赤外光脈波の脈
動成分の振幅と赤色光脈波の脈動成分の振幅とに差があるため、およびノイズが重畳して
いるためである。
列を用いてノイズを低減する演算処理について説明する。
なお、赤外光と赤色光とは交互に発光されるため厳密には同時に発光されるものではな
いが、隣り合う得られた赤外光受光値と赤色光受光値を同時刻に得られたものとして扱い
、所定時間分の赤外光の脈波信号と赤色光の脈波信号を２次元直交座標上に展開する。
すなわち図３のグラフを作成している。
また、脈波の直流成分に対する脈動成分の比をとることで脈拍による吸光度の脈動分が
近似される。
図３のグラフに見られる推移は４５度になっていないが、その理由は、赤外光脈波の脈
動成分の振幅と赤色光脈波の脈動成分の振幅とに差があるため、およびノイズが重畳して
いるためである。
次に、展開された脈波データに回転行列を用いて回転演算を施すこととする。
赤外光脈波の直流成分に対する脈動成分の比（ＩＲ）のデータ列を、
赤外光脈波の直流成分に対する脈動成分の比（ＩＲ）のデータ列を、
赤色光脈波の直流成分に対する脈動成分の比（Ｒ）のデータ列を、
とする。
同じ時刻tiに得られたＩＲとＲとのデータを次のように行列で定義する。
すなわち、
同じ時刻tiに得られたＩＲとＲとのデータを次のように行列で定義する。
すなわち、
また、θ[rad]回転させる回転行列をＡとすると、Ａは次のように表すことができる。
そうすると、ＳをＡによりθ[rad]回転させることにより次のＸが得られる。
なお、回転行列Ａは、上記のほか、
を用いてもよい。
ここで、θを０〜９π／３０[rad] までπ／３０[rad]ずつ脈波データＳを回転させて
得られるグラフを図４に示す。
図４からわかるように、横軸ゼロ、縦軸ゼロの点（赤色光脈波と赤外光脈波との両方が
平均である点）を中心として回転されており、θが９π／３０[rad]のときに、横軸(X1)
へ射影した領域が最小になり、縦軸(X2)へ射影した領域が最大となっている。
θを９π／３０[rad]よりさらにπ／２[rad]回転させ２４π／３０[rad] （＝１２π／
１５[rad]）回転させた場合には横軸(X1)へ射影した領域が最大になり、縦軸(X2)へ射影
した領域が最小となることは明らかである。
ここで、θを０〜９π／３０[rad] までπ／３０[rad]ずつ脈波データＳを回転させて
得られるグラフを図４に示す。
図４からわかるように、横軸ゼロ、縦軸ゼロの点（赤色光脈波と赤外光脈波との両方が
平均である点）を中心として回転されており、θが９π／３０[rad]のときに、横軸(X1)
へ射影した領域が最小になり、縦軸(X2)へ射影した領域が最大となっている。
θを９π／３０[rad]よりさらにπ／２[rad]回転させ２４π／３０[rad] （＝１２π／
１５[rad]）回転させた場合には横軸(X1)へ射影した領域が最大になり、縦軸(X2)へ射影
した領域が最小となることは明らかである。
次に、θを９π／３０[rad]、２４π／３０[rad]としたときの回転行列Ａにより、測定
された脈波データＳが処理されてＸとなった結果、どのような波形となるかを説明する。
図５は、図２に示した脈波データＳを、θを９π／３０[rad]として回転行列Ａにより
処理したＸの波形を示す。
横軸へ射影した領域が最小になったX1(t i)は、
された脈波データＳが処理されてＸとなった結果、どのような波形となるかを説明する。
図５は、図２に示した脈波データＳを、θを９π／３０[rad]として回転行列Ａにより
処理したＸの波形を示す。
横軸へ射影した領域が最小になったX1(t i)は、
一方、横軸へ射影した領域が最大になったX2(t i)は、
により演算される。
図５のX1の波形からはノイズが除去されたことがわかる。
なお、脈波データＳを、θを２４π／３０[rad]として回転行列Ａにより処理した場合
には、X2の波形がノイズが除去された波形となる。
横軸へ射影した領域が最大になるX1(t i)は、
図５のX1の波形からはノイズが除去されたことがわかる。
なお、脈波データＳを、θを２４π／３０[rad]として回転行列Ａにより処理した場合
には、X2の波形がノイズが除去された波形となる。
横軸へ射影した領域が最大になるX1(t i)は、
一方、横軸へ射影した領域が最小になるX2(t i)は、
により演算される。
このように横軸へ射影した領域が最小になるように回転角θを設定して、脈波データSを処理すれば、ノイズが抑制された脈波主成分波形を得ることができる。
このように横軸へ射影した領域が最小になるように回転角θを設定して、脈波データSを処理すれば、ノイズが抑制された脈波主成分波形を得ることができる。
次に、脈波の基本周波数の演算について説明する。
ノイズが低減される前の図２に示した脈波信号と、回転行列を用いてノイズが低減された
脈波主成分波形を周波数解析して得られたスペクトルをそれぞれ図６に示す。
横軸は周波数、縦軸はスペクトルである。
ノイズが低減される前の脈波(Before-rotation)信号のスペクトルは、ノイズの周波数帯
域ｆｎのスペクトルが強くでており、脈波信号の基本周波数ｆｓのスペクトルはほとんど
現れていない。
一方、回転行列を用いてノイズが低減された脈波主成分波形(After-rotation)を周波数
解析して得られたスペクトルでは、脈波信号の基本周波数ｆｓのスペクトルがノイズの周
波数帯域ｆｎのスペクトルと区別できるほど強く現れていることがわかり、脈波信号の基
本周波数ｆｓを求めることができる。
そして、脈波信号の基本周波数ｆｓ[Hz]が求まれば、脈拍数ｆｓ×６０[回／min]を容
易に求めることができる。
ノイズが低減される前の図２に示した脈波信号と、回転行列を用いてノイズが低減された
脈波主成分波形を周波数解析して得られたスペクトルをそれぞれ図６に示す。
横軸は周波数、縦軸はスペクトルである。
ノイズが低減される前の脈波(Before-rotation)信号のスペクトルは、ノイズの周波数帯
域ｆｎのスペクトルが強くでており、脈波信号の基本周波数ｆｓのスペクトルはほとんど
現れていない。
一方、回転行列を用いてノイズが低減された脈波主成分波形(After-rotation)を周波数
解析して得られたスペクトルでは、脈波信号の基本周波数ｆｓのスペクトルがノイズの周
波数帯域ｆｎのスペクトルと区別できるほど強く現れていることがわかり、脈波信号の基
本周波数ｆｓを求めることができる。
そして、脈波信号の基本周波数ｆｓ[Hz]が求まれば、脈拍数ｆｓ×６０[回／min]を容
易に求めることができる。
このように、所定角度の回転行列を用いることにより、ノイズが低減された脈波主成分
波形を得ることができ、脈波信号の基本周波数ないし脈拍数を求めることができる。
ここで、所定角度は、予め決められたものでもよく、測定期間中アダプティブに変化さ
せてもよい。
波形を得ることができ、脈波信号の基本周波数ないし脈拍数を求めることができる。
ここで、所定角度は、予め決められたものでもよく、測定期間中アダプティブに変化さ
せてもよい。
＜酸素飽和度の演算＞
図３は、上述のように図９に示される赤外光ＩＲのデータを横軸に、赤色光Ｒのデータ
を縦軸にとったグラフであるが、このグラフの傾きをノルム比を用いて求める。
まず、赤外光脈波データＩＲのL２ノルムを求める。
赤外光脈波データ列はIR = [ IR(ti) ： ti = 0, 1, 2, 3, ・・・]（１）であるから、L２ノルムは次の式で表すことができる。
図３は、上述のように図９に示される赤外光ＩＲのデータを横軸に、赤色光Ｒのデータ
を縦軸にとったグラフであるが、このグラフの傾きをノルム比を用いて求める。
まず、赤外光脈波データＩＲのL２ノルムを求める。
赤外光脈波データ列はIR = [ IR(ti) ： ti = 0, 1, 2, 3, ・・・]（１）であるから、L２ノルムは次の式で表すことができる。
次に、赤色光脈波データＲのL２ノルムを求める。
赤色光脈波データ列はR = [ R(ti) ： ti = 0, 1, 2, 3, ・・・]（２）であるから、L２ノルムは次の式で表すことができる。
赤色光脈波データ列はR = [ R(ti) ： ti = 0, 1, 2, 3, ・・・]（２）であるから、L２ノルムは次の式で表すことができる。
そこで、
とすればΦは、酸素飽和度SｐO2と相関するので、その相関を表す関数をｆとすれば、
と表され、酸素飽和度ＳｐＯ2を求めることができる。
なお、ノルム比を傾きとした直線を図３に示す。
ノルムとは、数学的概念の１つで、ユークリッドノルム（Euclidean-norm）又は２乗ノ
ルムは、ｎ個の要素を持つベクトルの大きさをスカラ量に写像するものである。
このように、所定期間の赤色光脈波データＲのL２ノルム値（２乗ノルム）と赤外光脈波データのL２ノルム値の比に基づいて、酸素飽和度SｐＯ2を求めることができる。
ここで、所定期間は逐次得られる現在の脈波から過去にさかのぼって所定期間分の赤色
光脈波データＲ、赤外光脈波データＩＲを用いるとよい。
また、ノルム値として、L２ノルムを用いたが、他の演算方法によるノルム値を用いて
もよい。
なお、ノルム比を傾きとした直線を図３に示す。
ノルムとは、数学的概念の１つで、ユークリッドノルム（Euclidean-norm）又は２乗ノ
ルムは、ｎ個の要素を持つベクトルの大きさをスカラ量に写像するものである。
このように、所定期間の赤色光脈波データＲのL２ノルム値（２乗ノルム）と赤外光脈波データのL２ノルム値の比に基づいて、酸素飽和度SｐＯ2を求めることができる。
ここで、所定期間は逐次得られる現在の脈波から過去にさかのぼって所定期間分の赤色
光脈波データＲ、赤外光脈波データＩＲを用いるとよい。
また、ノルム値として、L２ノルムを用いたが、他の演算方法によるノルム値を用いて
もよい。
また、酸素飽和度の演算に関しては、
脈波信号に対して、ノイズ信号が比較的小なる場合は上記ノルム比を用いて演算しても良
いが、比較的大なる場合は、上記のノルム比を用いて求める方法とは別に、前記特願２０
０１−３３２３８３号に記載した周波数解析を用いて求めた基本周波数に替えて、前記回
転による処理により求めた基本周波数を用いて演算することもできる。
脈波信号に対して、ノイズ信号が比較的小なる場合は上記ノルム比を用いて演算しても良
いが、比較的大なる場合は、上記のノルム比を用いて求める方法とは別に、前記特願２０
０１−３３２３８３号に記載した周波数解析を用いて求めた基本周波数に替えて、前記回
転による処理により求めた基本周波数を用いて演算することもできる。
（第１の実施例）
次に、上記原理を用いた装置を、概略構成ブロック図と処理フローにより説明する。
概略構成ブロック図は先に説明した図１と同じである。
発光素子１、２から、交互に発光するように駆動回路３により駆動されることにより、
異なる波長の光が発光される。
これらの発光素子１、２からの発光は生体組織４を透過して受光部（フォトダイオード
）５で受光され、電気信号に変換される。
そして、これらの変換された信号は増幅器６で増幅され、マルチプレクサ７によりそれ
ぞれの光波長に対応したフィルタ８−１、８−２に振り分けられる。
各フィルターに振り分けられた信号はフィルタ８−１、８−２によりフィルタリングさ
れてノイズ成分が低減され、A／Ｄ変換器９によりデジタル化される。
デジタル化された赤外光、赤色光に対応する各信号列が、それぞれの脈波を形成してい
る。
デジタル化された各信号列は処理部１０に入力され、ＲＯＭ１２に格納されているプロ
グラムにより処理され、脈拍数ＰＲ、酸素飽和度ＳｐＯ2が演算され、その値が表示部１
１に表示される。
次に、上記原理を用いた装置を、概略構成ブロック図と処理フローにより説明する。
概略構成ブロック図は先に説明した図１と同じである。
発光素子１、２から、交互に発光するように駆動回路３により駆動されることにより、
異なる波長の光が発光される。
これらの発光素子１、２からの発光は生体組織４を透過して受光部（フォトダイオード
）５で受光され、電気信号に変換される。
そして、これらの変換された信号は増幅器６で増幅され、マルチプレクサ７によりそれ
ぞれの光波長に対応したフィルタ８−１、８−２に振り分けられる。
各フィルターに振り分けられた信号はフィルタ８−１、８−２によりフィルタリングさ
れてノイズ成分が低減され、A／Ｄ変換器９によりデジタル化される。
デジタル化された赤外光、赤色光に対応する各信号列が、それぞれの脈波を形成してい
る。
デジタル化された各信号列は処理部１０に入力され、ＲＯＭ１２に格納されているプロ
グラムにより処理され、脈拍数ＰＲ、酸素飽和度ＳｐＯ2が演算され、その値が表示部１
１に表示される。
次に、脈拍数ＰＲ、酸素飽和度ＳｐＯ2を演算する処理フローを図７を用いて説明する。
測定が開始される（ステップＳ１）と、上記のように赤色光脈波、赤外光脈波が検出され
（ステップＳ２）、デジタル化された各信号列（各脈波データ）が処理部１０に取り込ま
れる。
処理部１０では、ＲＯＭ１２に格納されているプログラムにより、処理過程のデータを
ＲＡＭ１３に読み書きしながら、各脈波データを次のように処理する。
測定が開始される（ステップＳ１）と、上記のように赤色光脈波、赤外光脈波が検出され
（ステップＳ２）、デジタル化された各信号列（各脈波データ）が処理部１０に取り込ま
れる。
処理部１０では、ＲＯＭ１２に格納されているプログラムにより、処理過程のデータを
ＲＡＭ１３に読み書きしながら、各脈波データを次のように処理する。
先ず、赤外光脈波、赤色光脈波それぞれの脈波の直流成分に対する脈動成分の比を脈波ご
とに求める。（ステップＳ３）
次に、脈拍数ＰＲを求める処理（ステップＳ４〜Ｓ６）と酸素飽和度ＳｐＯ2を求める
処理（ステップＳ７〜Ｓ９）が同時に行われる。
とに求める。（ステップＳ３）
次に、脈拍数ＰＲを求める処理（ステップＳ４〜Ｓ６）と酸素飽和度ＳｐＯ2を求める
処理（ステップＳ７〜Ｓ９）が同時に行われる。
脈拍数ＰＲを求める処理（ステップＳ４〜Ｓ６）では、
予め回転角が設定された回転行列Ａにより、赤外光脈波データＩＲと赤色光脈波データ
ＲとのデータＳから、式５よりノイズが低減された波形を得る。（ステップＳ４）
設定する回転角は、図３に示すような赤外光脈波データＩＲを横軸に赤色光脈波データ
Ｒを縦軸にとったグラフを回転すると図４に示すように軸方向へ射影した領域が最小とな
るような角度である。
回転角度は例えば、９π／３０[rad]あるいは２４π／３０[rad]がよい。
ノイズが低減された波形は、射影領域が最小となる軸成分のデータより得ることができ
る。
そして、ノイズが低減された波形を図６に示すように周波数解析を行い、脈波データの基
本周波数を求める。（ステップＳ５）
そして、その基本周波数から脈拍数ｆｓをｆＳ×６０[回／min]から求め、表示部１１
に表示する。
予め回転角が設定された回転行列Ａにより、赤外光脈波データＩＲと赤色光脈波データ
ＲとのデータＳから、式５よりノイズが低減された波形を得る。（ステップＳ４）
設定する回転角は、図３に示すような赤外光脈波データＩＲを横軸に赤色光脈波データ
Ｒを縦軸にとったグラフを回転すると図４に示すように軸方向へ射影した領域が最小とな
るような角度である。
回転角度は例えば、９π／３０[rad]あるいは２４π／３０[rad]がよい。
ノイズが低減された波形は、射影領域が最小となる軸成分のデータより得ることができ
る。
そして、ノイズが低減された波形を図６に示すように周波数解析を行い、脈波データの基
本周波数を求める。（ステップＳ５）
そして、その基本周波数から脈拍数ｆｓをｆＳ×６０[回／min]から求め、表示部１１
に表示する。
酸素飽和度ＳｐＯ２を求める処理（ステップＳ７〜Ｓ９）では、所定期間分の赤外光脈
波データＩＲと赤色光脈波データＲとからそれぞれのＬ２ノルム値を式１０、式１１から
求め、さらにそれぞれのＬ２ノルム値の比を式１２から求める。
次に、ノイズの除去された赤外光と赤色光の脈波信号の比を求め、酸素飽和度を演算する。
（ステップＳ７）
そのＬ２ノルム比をΦとして、式１２より酸素飽和度ＳｐＯ２を求め（ステップＳ８）
、表示部１１に表示する。（ステップＳ９）
波データＩＲと赤色光脈波データＲとからそれぞれのＬ２ノルム値を式１０、式１１から
求め、さらにそれぞれのＬ２ノルム値の比を式１２から求める。
次に、ノイズの除去された赤外光と赤色光の脈波信号の比を求め、酸素飽和度を演算する。
（ステップＳ７）
そのＬ２ノルム比をΦとして、式１２より酸素飽和度ＳｐＯ２を求め（ステップＳ８）
、表示部１１に表示する。（ステップＳ９）
測定を継続するときはステップＳ２に戻り処理を繰り返し、測定を計測しない場合は測
定を終了する。（ステップＳ１１）
定を終了する。（ステップＳ１１）
（第２の実施例）
次に、別の第２の実施例を図８を用いて説明する。
第２の実施例が第１の実施例と相違する点はステップＳ４において、回転角は予め定め
られたものではなく、得られるデータから回転角を求める点であり、図８に示すようにス
テップＳ４−１とステップＳ４−２に分けて処理する。
他のステップは第１の実施例と同様なので、説明は省略する。
次に、別の第２の実施例を図８を用いて説明する。
第２の実施例が第１の実施例と相違する点はステップＳ４において、回転角は予め定め
られたものではなく、得られるデータから回転角を求める点であり、図８に示すようにス
テップＳ４−１とステップＳ４−２に分けて処理する。
他のステップは第１の実施例と同様なので、説明は省略する。
脈拍数ＰＲを求める処理（ステップＳ４−１〜Ｓ６）では、
先ず、所定期間分の赤外光脈波データＩＲと赤色光脈波データＲとのデータを用い、図
３に示すようなグラフを描く。
そして、いかなる回転角度であれば軸方向へ射影した領域が最小になるかを求める。（
ステップＳ４−１）
次に、求められた回転角度による回転行列により、各波長の脈波データを処理し、射影
領域が最小となる軸の成分のデータからノイズが低減された波形を得る。（ステップＳ４
−２）。
このように、第２の実施例の特徴は、回転行列の回転角度が、固定化された角度ではな
く、検出される脈波データに応じて適時変更されるアダプティブ性を有する点にある。
先ず、所定期間分の赤外光脈波データＩＲと赤色光脈波データＲとのデータを用い、図
３に示すようなグラフを描く。
そして、いかなる回転角度であれば軸方向へ射影した領域が最小になるかを求める。（
ステップＳ４−１）
次に、求められた回転角度による回転行列により、各波長の脈波データを処理し、射影
領域が最小となる軸の成分のデータからノイズが低減された波形を得る。（ステップＳ４
−２）。
このように、第２の実施例の特徴は、回転行列の回転角度が、固定化された角度ではな
く、検出される脈波データに応じて適時変更されるアダプティブ性を有する点にある。
（第３実施例）
次に、脈拍数ＰＲ、酸素飽和度ＳｐＯ2を、周波数解析を用いて求めた基本周波数に替
えて、前記回転による処理により求めた基本周波数を用いて演算する処理フローを図１０
を用いて説明する。
測定が開始される（ステップＳ１）と、上記のように赤色光脈波、赤外光脈波が検出さ
れ（ステップＳ２）、デジタル化された各信号列（各脈波データ）が処理部１０に取り込
まれる。
処理部１０では、ＲＯＭ１２に格納されているプログラムにより、処理過程のデータを
ＲＡＭ１３に読み書きしながら、各脈波データを次のように処理する。
次に、脈拍数ＰＲ、酸素飽和度ＳｐＯ2を、周波数解析を用いて求めた基本周波数に替
えて、前記回転による処理により求めた基本周波数を用いて演算する処理フローを図１０
を用いて説明する。
測定が開始される（ステップＳ１）と、上記のように赤色光脈波、赤外光脈波が検出さ
れ（ステップＳ２）、デジタル化された各信号列（各脈波データ）が処理部１０に取り込
まれる。
処理部１０では、ＲＯＭ１２に格納されているプログラムにより、処理過程のデータを
ＲＡＭ１３に読み書きしながら、各脈波データを次のように処理する。
先ず、赤外光脈波、赤色光脈波それぞれの脈波の直流成分に対する脈動成分の比を脈波
ごとに求める。（ステップＳ３）
次に、脈拍数ＰＲを求める処理（ステップＳ４〜Ｓ６）と酸素飽和度ＳｐＯ2を求める
処理（ステップＳ７〜Ｓ９）が同時に行われる。
ごとに求める。（ステップＳ３）
次に、脈拍数ＰＲを求める処理（ステップＳ４〜Ｓ６）と酸素飽和度ＳｐＯ2を求める
処理（ステップＳ７〜Ｓ９）が同時に行われる。
脈拍数ＰＲを求める処理（ステップＳ４〜Ｓ６）では、
予め回転角が設定された回転行列Ａにより、赤外光脈波データＩＲと赤色光脈波データ
ＲとのデータＳから、式５よりノイズが低減された波形を得る。（ステップＳ４）
設定する回転角は、図３に示すような赤外光脈波データＩＲを横軸に赤色光脈波データ
Ｒを縦軸にとったグラフを回転すると図４に示すように軸方向へ射影した領域が最小とな
るような角度である。
回転角度は例えば、９π／３０[rad]あるいは２４π／３０[rad]がよい。
ノイズが低減された波形は、射影領域が最小となる軸成分のデータより得ることができ
る。
そして、ノイズが低減された波形を図６に示すように周波数解析を行い、脈波データの
基本周波数を求める。（ステップＳ５）
更に、その基本周波数から脈拍数ｆｓをｆＳ×６０[回／min]から求め、表示部１１に
表示する。
予め回転角が設定された回転行列Ａにより、赤外光脈波データＩＲと赤色光脈波データ
ＲとのデータＳから、式５よりノイズが低減された波形を得る。（ステップＳ４）
設定する回転角は、図３に示すような赤外光脈波データＩＲを横軸に赤色光脈波データ
Ｒを縦軸にとったグラフを回転すると図４に示すように軸方向へ射影した領域が最小とな
るような角度である。
回転角度は例えば、９π／３０[rad]あるいは２４π／３０[rad]がよい。
ノイズが低減された波形は、射影領域が最小となる軸成分のデータより得ることができ
る。
そして、ノイズが低減された波形を図６に示すように周波数解析を行い、脈波データの
基本周波数を求める。（ステップＳ５）
更に、その基本周波数から脈拍数ｆｓをｆＳ×６０[回／min]から求め、表示部１１に
表示する。
また、酸素飽和度ＳｐＯ２を求める処理（ステップＳ７〜Ｓ９）では、
赤外光及び赤色光の脈波信号を前記基本周波数で、または基本周波数とその高調波で構
成したフィルタを通し、ノイズの除去された信号を求める。（ステップＳ７）
次に、ノイズの除去された赤外光と赤色光の脈波信号の比を求め、酸素飽和度を演算し
（ステップＳ８）表示部１１に表示する。（ステップＳ９）
赤外光及び赤色光の脈波信号を前記基本周波数で、または基本周波数とその高調波で構
成したフィルタを通し、ノイズの除去された信号を求める。（ステップＳ７）
次に、ノイズの除去された赤外光と赤色光の脈波信号の比を求め、酸素飽和度を演算し
（ステップＳ８）表示部１１に表示する。（ステップＳ９）
測定を継続するとき（ステップＳ１０）は、ステップＳ２に戻り処理を繰り返し、測定
を計測しない場合は測定を終了する。（ステップＳ１１）
を計測しない場合は測定を終了する。（ステップＳ１１）
（第４の実施例）
更に、別の第４の実施例を図１１を用いて説明する。
第４の実施例が第３の実施例と相違する点はステップＳ４において、回転角は予め定めら
れたものではなく、得られるデータから回転角を求める点であり、図１１に示すようにス
テップＳ４−１とステップＳ４−２に分けて処理する。
他のステップは第３の実施例と同様なので、説明は省略する。
更に、別の第４の実施例を図１１を用いて説明する。
第４の実施例が第３の実施例と相違する点はステップＳ４において、回転角は予め定めら
れたものではなく、得られるデータから回転角を求める点であり、図１１に示すようにス
テップＳ４−１とステップＳ４−２に分けて処理する。
他のステップは第３の実施例と同様なので、説明は省略する。
脈拍数ＰＲを求める処理（ステップＳ４−１〜Ｓ６）では、
先ず、所定期間分の赤外光脈波データＩＲと赤色光脈波データＲとのデータを用い、図
３に示すようなグラフを描く。
そして、いかなる回転角度であれば軸方向へ射影した領域が最小になるかを求める。（
ステップＳ４−１）
次に、求められた回転角度による回転行列により、各波長の脈波データを処理し、射影
領域が最小となる軸の成分のデータからノイズが低減された波形を得る。（ステップＳ４
−２）。
このように、第４の実施例の特徴は、回転行列の回転角度が、固定化された角度ではな
く、検出される脈波データに応じて適時変更されるアダプティブ性を有する点にある。
先ず、所定期間分の赤外光脈波データＩＲと赤色光脈波データＲとのデータを用い、図
３に示すようなグラフを描く。
そして、いかなる回転角度であれば軸方向へ射影した領域が最小になるかを求める。（
ステップＳ４−１）
次に、求められた回転角度による回転行列により、各波長の脈波データを処理し、射影
領域が最小となる軸の成分のデータからノイズが低減された波形を得る。（ステップＳ４
−２）。
このように、第４の実施例の特徴は、回転行列の回転角度が、固定化された角度ではな
く、検出される脈波データに応じて適時変更されるアダプティブ性を有する点にある。
以上は、動脈血酸素飽和度を測定するパルスオキシメータを例に挙げて説明したが、本
発明の技術はパルスオキシメータに限られず、特殊ヘモグロビン（一酸化炭素ヘモグロビ
ン、Ｍｅｔヘモグロビンなど）、血中に注入された色素などの血中吸光物質をパルスフォ
トメトリーの原理を用いて測定する装置(パルスフォトメトリー)にも光源の波長を選択す
ることで適用できる。
なお、本明細書中において、生体組織（媒体）を介して測定された信号を、２次元直交
座標に展開する前までは、「脈波信号」と記載し、２次元直交座標に展開中は、「脈波デ
ータ」と記載し、求められた波形は「脈波主成分波形」と区別して記載している。
発明の技術はパルスオキシメータに限られず、特殊ヘモグロビン（一酸化炭素ヘモグロビ
ン、Ｍｅｔヘモグロビンなど）、血中に注入された色素などの血中吸光物質をパルスフォ
トメトリーの原理を用いて測定する装置(パルスフォトメトリー)にも光源の波長を選択す
ることで適用できる。
なお、本明細書中において、生体組織（媒体）を介して測定された信号を、２次元直交
座標に展開する前までは、「脈波信号」と記載し、２次元直交座標に展開中は、「脈波デ
ータ」と記載し、求められた波形は「脈波主成分波形」と区別して記載している。
本発明では、同一の媒体からほぼ同時に抽出される２つの同種の信号
を処理して共通の信号成分を抽出する計算処理負担を軽減した信号処理が実現できる。
また、本発明では、前記媒体の体動によるノイズが脈波信号に生じ
た場合であっても、対象物質の濃度を精度よく測定できる。
また、体動によるノイズが脈波データ信号に生じた場合であっても、脈波信号からノイ
ズを除去し、精度よく脈拍や血中の吸光物質の濃度を求めることができるので、産業上の
利用可能性は極めて大きい。
を処理して共通の信号成分を抽出する計算処理負担を軽減した信号処理が実現できる。
また、本発明では、前記媒体の体動によるノイズが脈波信号に生じ
た場合であっても、対象物質の濃度を精度よく測定できる。
また、体動によるノイズが脈波データ信号に生じた場合であっても、脈波信号からノイ
ズを除去し、精度よく脈拍や血中の吸光物質の濃度を求めることができるので、産業上の
利用可能性は極めて大きい。
１ 発光素子
２ 発光素子
３ 駆動回路
４ 生体組織
５ フォトダイオード
６ 変換器
７ マルチプレクサ
８ フィルタ
９ Ａ／Ｄ変換器
１０ 処理部
１１ 表示部
１２ ＲＯＭ
１３ ＲＡＭ
２ 発光素子
３ 駆動回路
４ 生体組織
５ フォトダイオード
６ 変換器
７ マルチプレクサ
８ フィルタ
９ Ａ／Ｄ変換器
１０ 処理部
１１ 表示部
１２ ＲＯＭ
１３ ＲＡＭ
Claims (1)
- 異なる２つの波長の光を同一の生体組織に照射する発光手段と、
前記発光手段から発生し前記同一の生体組織を透過または反射した各波長の光を電気信号に変換する受光手段と、
前記受光手段で変換された前記電気信号より得られた２つの波長の離散的時系列脈波データをそれぞれ縦軸または横軸とする２次元直交座標に展開した、各波長の脈波データそれぞれについてノルム値を求め、さらにそのノルム値の比を求めるノルム比計算手段と、
前記ノルム比計算手段により計算されたノルム比に基づいて、血中の酸素の濃度を求める血中酸素濃度演算手段とを具備することを特徴とするパルスフォトメータ。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2009006876A JP4807598B2 (ja) | 2002-10-31 | 2009-01-15 | パルスフォトメータ |
Applications Claiming Priority (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002318278 | 2002-10-31 | ||
| JP2002318278 | 2002-10-31 | ||
| JP2002362261 | 2002-12-13 | ||
| JP2002362261 | 2002-12-13 | ||
| JP2009006876A JP4807598B2 (ja) | 2002-10-31 | 2009-01-15 | パルスフォトメータ |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2003333613A Division JP4352315B2 (ja) | 2002-10-31 | 2003-09-25 | 信号処理方法／装置及びそれを用いたパルスフォトメータ |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2009072639A true JP2009072639A (ja) | 2009-04-09 |
| JP4807598B2 JP4807598B2 (ja) | 2011-11-02 |
Family
ID=40608200
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2009006876A Expired - Fee Related JP4807598B2 (ja) | 2002-10-31 | 2009-01-15 | パルスフォトメータ |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4807598B2 (ja) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH07308299A (ja) * | 1994-05-18 | 1995-11-28 | Tetsuji Fukamizu | 脈波検出方法及び装置 |
| JPH11183377A (ja) * | 1997-12-17 | 1999-07-09 | Matsushita Electric Ind Co Ltd | 光学式成分計 |
| JP2001078990A (ja) * | 1999-09-13 | 2001-03-27 | Nippon Koden Corp | 血中吸光物質測定装置 |
| JP2001120511A (ja) * | 1999-10-26 | 2001-05-08 | Shimadzu Corp | 生体信号計測装置 |
| JP2003235819A (ja) * | 2001-12-14 | 2003-08-26 | Nippon Koden Corp | 信号処理方法および脈波信号処理方法 |
-
2009
- 2009-01-15 JP JP2009006876A patent/JP4807598B2/ja not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH07308299A (ja) * | 1994-05-18 | 1995-11-28 | Tetsuji Fukamizu | 脈波検出方法及び装置 |
| JPH11183377A (ja) * | 1997-12-17 | 1999-07-09 | Matsushita Electric Ind Co Ltd | 光学式成分計 |
| JP2001078990A (ja) * | 1999-09-13 | 2001-03-27 | Nippon Koden Corp | 血中吸光物質測定装置 |
| JP2001120511A (ja) * | 1999-10-26 | 2001-05-08 | Shimadzu Corp | 生体信号計測装置 |
| JP2003235819A (ja) * | 2001-12-14 | 2003-08-26 | Nippon Koden Corp | 信号処理方法および脈波信号処理方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP4807598B2 (ja) | 2011-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4352315B2 (ja) | 信号処理方法／装置及びそれを用いたパルスフォトメータ | |
| JP4632143B2 (ja) | 信号処理方法及びそれを適用したパルスフォトメータ | |
| JP4830693B2 (ja) | 酸素飽和度測定装置及び測定方法 | |
| JP4196209B2 (ja) | 信号処理方法及びそれを適用したパルスフォトメータ | |
| JP4454854B2 (ja) | パルス式計測システムにおける信号復調方法および装置 | |
| US6650918B2 (en) | Cepstral domain pulse oximetry | |
| JP3627214B2 (ja) | 血中吸光物質測定装置 | |
| US6101404A (en) | Optical diagnosis point detector for noninvasive diagnosis of blood constituents and noninvasive diagnostic device | |
| US20040039273A1 (en) | Cepstral domain pulse oximetry | |
| JP2004248819A (ja) | 血液分析装置 | |
| DE60016445D1 (de) | Apparat zur entfernung von bewegungs-artefakten von körperparametermessungen | |
| CN1665443A (zh) | 用于增强信噪比的信号处理方法和设备 | |
| JP5203210B2 (ja) | パルスオキシメトリーのための信号処理 | |
| JP4831111B2 (ja) | 信号処理方法及びそれを用いたパルスフォトメータ | |
| KR101656746B1 (ko) | 혈당농도 보정 장치 및 방법 | |
| JP2813240B2 (ja) | 末梢循環検出装置 | |
| JP2003235819A (ja) | 信号処理方法および脈波信号処理方法 | |
| JP4470056B2 (ja) | パルスフォトメータに用いられる信号処理方法及びそれを用いたパルスフォトメータ | |
| JP4807598B2 (ja) | パルスフォトメータ | |
| CN111631733B (zh) | 一种动脉血液光谱检测方法及装置 | |
| JP2012019967A (ja) | ノイズ除去方法及びパルスフォトメータ | |
| JP4040224B2 (ja) | 血管の画像化装置並びに脈波信号の空間分布測定方法および装置 | |
| JP2012019966A (ja) | ノイズ除去方法及びパルスフォトメータ | |
| JP2012019968A (ja) | ノイズ除去方法及びパルスフォトメータ | |
| Herzog et al. | Contactless measurement of dermal oxygen saturation with spatial resolution |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20090115 |
|
| 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: 20110712 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20110803 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140826 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 4807598 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |