パルスフォトメータ
【課題】 同一の媒体からほぼ同時に抽出される２つの同種の信号を処
理して共通の信号成分を抽出する計算処理負担を軽減したパルスフォトメータを提供する。
【解決手段】異なる２つの波長の光を同一の生体組織に照射する発光手段と、前記発光手段から発生し前記同一の生体組織を透過または反射した各波長の光を電気信号に変換する受光手段と、前記受光手段で変換された前記電気信号より得られた２つの波長の離散的時系列脈波データをそれぞれ縦軸または横軸とする２次元直交座標に展開した、各波長の脈波データそれぞれについてノルム値を求め、さらにそのノルム値の比を求めるノルム比計算手段と、前記ノルム比計算手段により計算されたノルム比に基づいて、血中の酸素の濃度を求める血中酸素濃度演算手段とを具備することを特徴とするパルスフォトメータ。
理して共通の信号成分を抽出する計算処理負担を軽減したパルスフォトメータを提供する。
【解決手段】異なる２つの波長の光を同一の生体組織に照射する発光手段と、前記発光手段から発生し前記同一の生体組織を透過または反射した各波長の光を電気信号に変換する受光手段と、前記受光手段で変換された前記電気信号より得られた２つの波長の離散的時系列脈波データをそれぞれ縦軸または横軸とする２次元直交座標に展開した、各波長の脈波データそれぞれについてノルム値を求め、さらにそのノルム値の比を求めるノルム比計算手段と、前記ノルム比計算手段により計算されたノルム比に基づいて、血中の酸素の濃度を求める血中酸素濃度演算手段とを具備することを特徴とするパルスフォトメータ。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、一つの媒体からほぼ同時に抽出される２つの同種の信号を処理して共通の信
号成分を抽出する信号処理に関し、特には医療の分野において、特に循環器系の診断に用
いられるパルスフォトメータにおける信号処理の改良に関する。
【背景技術】
【０００２】
一つの媒体からほぼ同時に抽出された２つの信号から信号成分と雑音成分に分離する方
法には様々な方法が提案されている。
それらは、一般的には周波数領域や時間領域による処理が行われている。
医療現場でも、光電脈波計と言われる脈波波形や脈拍数を測定する装置、血液に含まれ
る吸光物質の濃度測定として、酸素飽和度ＳｐＯ2の測定装置、一酸化炭素ヘモグロビン
やＭｅｔヘモグロビン等の特殊ヘモグロビンの濃度の測定装置、注入色素濃度の測定装置
などがパルスフォトメータとして知られている。
中でも酸素飽和度ＳｐＯ2の測定装置を特にパルスオキシメータと呼んでいる。
【０００３】
パルスフォトメータの原理は、対象物質への吸光性が異なる複数の波長の光を生体組織
に透過又は反射させ、その透過光又は反射光の光量を連続的に測定することで得られる脈
波データ信号から対象物質の濃度を求めるものである。
そしてその脈波データに雑音が混入すると、正しい濃度の計算が出来ず、誤処置につな
がる危険が生じる。
パルスフォトメータにおいても従来より雑音を低減するために周波数帯域を分割して信
号成分に着目したり、２つの信号の相関を取るなどの方法が提案されてきた。
しかし、これらの方法は解析に時間がかかるなどの問題があった。
【０００４】
そこで、本出願人は、特許第３２７０９１７号（特許文献１）において、異なる２つの
波長の光を生体組織に照射して透過光から得られる２つの脈波信号のそれぞれの大きさを
縦軸、横軸としてグラフを描き、その回帰直線を求め、その回帰直線の傾きに基づいて、
動脈血中の酸素飽和度ないし吸光物質濃度を求めることを提案している。
この発明により、測定精度を高め、低消費電力化することができた。
しかし、各波長の脈波信号についての多くのサンプリングデータを用いて回帰直線ない
しその傾きを求めるためには、なお多くの計算処理を要していた。
【０００５】
更に本出願人は、特願２００１−３３２３８３号（特許文献２）においては、周波数解
析を用いてはいるが、その解析においては従来技術のように脈波信号そのものを抽出する
のではなく、脈波信号の基本周波数を求め、さらには精度を高めるためにその高調波周波
数を用いたフィルタを用いて脈波信号をフィルタリングする方法を提案している。
しかし、基本周波数を求める点に関しては更なる改善が望まれていた。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第３２７０９１７号 （請求項１、２、図２、図４）
【特許文献２】特願２００３−１３５４３４号（０００６）
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明の課題（目的）は、同一の媒体からほぼ同時に抽出される２つの同種の信号を処
理して共通の信号成分を抽出する計算処理負担を軽減したパルスフォトメータを提供することにある。
また、上記信号処理方法を適用して、前記媒体の体動によるノイズが脈波信号に生じた
場合であっても、対象物質の濃度を精度よく測定することにある。
また、体動によるノイズが脈波データ信号に生じた場合であっても、脈波信号からノイ
ズを除去し、精度よく脈拍を求めることにある。
【課題を解決するための手段】
【０００８】
・前記課題を解決するために、本発明のパルスフォトメータは、異なる２つの波長の光を同一の生体組織に照射する発光手段と、前記発光手段から発生し前記同一の生体組織を透過または反射した各波長の光を電気信号に変換する受光手段と、前記受光手段で変換された前記電気信号より得られた２つの波長の離散的時系列脈波データをそれぞれ縦軸または横軸とする２次元直交座標に展開した、各波長の脈波データそれぞれについてノルム値を求め、さらにそのノルム値の比を求めるノルム比計算手段と、前記ノルム比計算手段により計算されたノルム比に基づいて、血中の酸素の濃度を求める血中酸素濃度演算手段とを具備することを特徴とする。（請求項１）
【０００９】
また、参考例では、
・同一の媒体から抽出される所定期間の２つの信号を２次元直交座標上に展開するステップと、前記展開された信号を、回転行列を用いて、あらかじめ決められた角度または所定条件に基づいて決められた角度に回転させるステップとにより前記２つの信号を処理する。
この信号処理方法によって、同一の媒体からほぼ同時に抽出される２つの同種の信号を
処理して共通の信号成分を抽出する計算処理負担を軽減することができる。
・前記信号処理の結果、信号主成分の基本周波数を求める。
この構成により、ノイズが低減された基本波波形を容易に取得することができる。
・同一の媒体からの２つの信号を入力する入力手段と、前記入力手段から入力された所定
期間分の２つの信号を２次元直交座標上に展開する展開手段と、該展開されたデータをあ
らかじめ決められた角度または所定条件に基づいて決められた角度に回転させる回転行列
を用いて処理する処理手段とを具えた構成とする。
この構成により、回転行列を用いて処理することでノイズが低減された波形を容易に取
得することができる。
【００１０】
・異なる２つの波長の光を生体組織に照射する発光手段と、前記発光手段から発生し前記
生体組織を透過または反射した各波長の光を電気信号に変換する受光手段とを具えたパル
スフォトメータにおいて、
前記各波長の電気信号より得られた脈波データを、あらかじめ決められた角度にまたは
所定条件に基づいて決められた角度に回転させる回転行列を用いて処理する処理手段とを
具備することを特徴とする。
この構成により、脈波データを回転行列を用いて処理することでノイズが低減された波
形を容易に取得することができる。
・前記脈波データが、ほぼ同時点に受光した２つの波長の受光信号から得られる電気信号
を、前記２つの波長をそれぞれ縦軸または横軸とする２次元直交座標に展開したものであ
ることを特徴とする。
・前記脈波データは所定時間分であり、かつ経時移動して処理されることを特徴とする。
・異なる２つの波長の光を生体組織に照射する発光手段と、前記発光手段から発生し前記
生体組織を透過または反射した各波長の光を電気信号に変換する受光手段とを具えたパル
スフォトメータにおいて、前記各波長の電気信号より得られた脈波信号を、各波長をそれ
ぞれ縦軸または横軸とする２次元直交座標に展開し、縦軸または横軸に射影される領域が
最大または最小となるいずれかの条件を満足するように回転行列を用いて回転させて処理
する第２の処理手段とを具備することを特徴とする。
【００１１】
・異なる２つの波長の光を生体組織に照射する発光手段と、前記発光手段で発生し前記生
体組織を透過または反射した各波長の光を電気信号に変換する受光手段と、前記受光部に
より得られる前記２つの波長の脈波信号をそれぞれ縦軸または横軸とする２次元直交座標
に展開した、各波長の各波長の脈波データを、所定角度に回転させる回転行列を用いて処
理した波形を求める波形取得手段と、前記波形取得手段により得られた波形の周波数解析
により、脈波の基本周波数または脈拍数を求める波形解析手段と、前記波形解析手段の出
力から血中の吸光物質の濃度を求める血中吸光物質濃度演算手段とを具備するパルスフォ
トメータ。
・前記血中吸光物質濃度演算手段は、動脈血中の酸素飽和度、特殊ヘモグロビン濃度、ま
たは注入色素濃度のうち少なくとも１つを演算する。
【発明の効果】
【００１２】
請求項１に記載のパルスフォトメータによれば、各波長の脈波データから得られるノル
ム比から精度よく血中の酸素の濃度を求めることができる。
また、パルスフォトメータによれば、動脈血中の酸素飽和度、特殊ヘモグロビン濃度、または注入色素濃度のうち少なくとも１つを演算することができる。
【００１３】
参考例のパルスフォトメータによれば、脈波データを回転行列を用いて処理することでノイズが低減された波形を取得することができる。
【００１４】
参考例のパルスフォトメータによれば、脈波データを回転行列を用いて処理することでノイズが低減された波形を取得し、その波形から精度よく脈波データの基本周波数または脈拍数を求めることができる。
【００１５】
参考例のパルスフォトメータによれば、脈波データを回転行列を用いてノイズを除去するのに適した回転角度を決定でき、その決定した回転角度の回転行列によって処理することでノイズが低減された波形を取得することができる。
【００１６】
参考例のパルスフォトメータによれば、脈波データを回転行列を用いてノイズを除去するのに適した回転角度を決定でき、その決定した回転角度の回転行列によって処理することでノイズが低減された波形を取得し、その波形から精度よく脈波データの基本周波数または脈拍数を求めることができる。
【００１７】
参考例のパルスフォトメータによれば、回転行列の回転角の決定を容易に決定できる。
【００１８】
参考例のパルスフォトメータによれば、脈波データを回転行列を用いて処理することでノイズが低減された波形による精度の良い脈波の基本周波数または脈拍数の取得と、動脈血中の酸素飽和度、特殊ヘモグロビン濃度、または注入色素濃度のうち少なくとも１つの演算とが同一の装置でほぼ同時に処理して出力することができる。
参考例のパルスフォトメータによれば、動脈血中の酸素飽和度、特殊ヘモグロビン濃度、または注入色素濃度のうち少なくとも１つを演算することができる。
【図面の簡単な説明】
【００１９】
【図１】本発明の概略構成を示すブロック図である。
【図２】検出された脈波を示す図である。
【図３】図９に示される赤外光ＩＲのデータを横軸に、赤色光Ｒのデータを縦軸にと ったグラフである。
【図４】図３のグラフをπ／３０[rad]ずつ回転させた図である。
【図５】回転角度９π／３０[rad]の回転行列により処理された脈波の波形を示す図 である。
【図６】図５に示すX1の波形のスペクトルを示す図である。
【図７】第１の実施例における処理フローを示すフローチャートである。
【図８】第２の実施例における処理フローを示すフローチャートである。
【図９】血液中の吸光物質の吸光度の変動の測定原理を説明する波形図である。
【図１０】第３の実施例における処理フローを示すフローチャートである。
【図１１】第４の実施例における処理フローを示すフローチャートである。
【発明を実施するための形態】
【００２０】
本発明の実施の形態を説明するにあたり、動脈血酸素飽和度を測定するパルスオキシメ
ータを例に挙げて原理を説明する。
なお、本発明の技術は、パルスオキシメータに限られず、特殊ヘモグロビン（一酸化炭
素ヘモグロビン、Ｍｅｔヘモグロビンなど）、血中に注入された色素などの血中吸光物質
をパルスフォトメトリーの原理を用いて測定する装置(パルスフォトメータ)に適用できる。
【００２１】
動脈血酸素飽和度を測定するパルスオキシメータの構成は、概略構成ブロック図である
図１のようになっている。
異なる波長の光を発光する発光素子１、２は、交互に発光するように駆動回路３により
駆動される。
発光素子１、２に採用する光はそれぞれ動脈血酸素飽和度による影響が少ない赤外光（
例えば９４０[nm]）、動脈血酸素飽和度の変化に対する感度が高い赤色光（例えば６６０
[nm]）がよい。
【００２２】
これらの発光素子１、２からの発光は生体組織４を透過してフォトダイオード５で受光
して電気信号に変換される。
なお、反射光を受光するようにしてもよい。
そして、これらの変換された信号は増幅器６で増幅され、マルチプレクサ７によりそれ
ぞれの光波長に対応したフィルタ８−１、８−２に振り分けられる。
各フィルターに振り分けられた信号はフィルタ８−１、８−２によりフィルタリングさ
れてノイズ成分が低減され、Ａ／Ｄ変換器９によりデジタル化される。
【００２３】
デジタル化された赤外光、赤色光に対応する各信号列が、それぞれの脈波信号を形成し
ている。
デジタル化された各信号列は処理部１０に入力され、ＲＯＭ１２に格納されているプロ
グラムにより処理され、酸素飽和度ＳｐＯ2が測定され、その値が表示部１１に表示され
る。
【００２４】
＜回転行列によるノイズ低減と脈波の基本周波数の演算＞
先ず、血液中の吸光物質の吸光度（減光度）の変動の測定について説明する。
図９(a)及び(b)は、前記発光素子１、２からの発光された光が生体組織４を透過してフ
ォトダイオード５で受光して電気信号に変換された脈波データで、(a)は赤色光の場合を
、(b)は赤外光を示している。
図９の(a)では、横軸を時間、縦軸を受光出力とすると、フォトダイオード５での受光
出力は、赤色光の直流成分（Ｒ’）と脈動成分（ΔＲ’）が重畳された波形となっている。
また、図９の(b)では、横軸を時間、縦軸を受光出力とすると、フォトダイオード５で
の受光出力は、赤外光の直流成分（ＩＲ’）と脈動成分（ΔＩＲ’）が重畳された波形となっている。
図２は、図９に示すような脈波において、８秒間分の、直流成分（Ｒ’、ＩＲ’）に対
する脈動成分（ΔＲ’、ΔＩＲ’）の比（ＩＲ＝ΔＩＲ’／ＩＲ’）をとり、さらにその８秒間分のデータの平均値をゼロに合わせたものである。
なお、図２の如き、平均値をゼロとする処理を行わなくとも演算は可能である。
図３は、図９に示される赤外光ＩＲのデータを横軸に、赤色光Ｒのデータを縦軸にとっ
たグラフである。
【００２５】
次に、Ａ／Ｄ変換器９によってデジタル化した各波長の２つの脈波データ信号を回転行
列を用いてノイズを低減する演算処理について説明する。
なお、赤外光と赤色光とは交互に発光されるため厳密には同時に発光されるものではな
いが、隣り合う得られた赤外光受光値と赤色光受光値を同時刻に得られたものとして扱い
、所定時間分の赤外光の脈波信号と赤色光の脈波信号を２次元直交座標上に展開する。
すなわち図３のグラフを作成している。
また、脈波の直流成分に対する脈動成分の比をとることで脈拍による吸光度の脈動分が
近似される。
図３のグラフに見られる推移は４５度になっていないが、その理由は、赤外光脈波の脈
動成分の振幅と赤色光脈波の脈動成分の振幅とに差があるため、およびノイズが重畳して
いるためである。
【００２６】
次に、展開された脈波データに回転行列を用いて回転演算を施すこととする。
赤外光脈波の直流成分に対する脈動成分の比（ＩＲ）のデータ列を、
【００２７】
【数１】
【００２８】
赤色光脈波の直流成分に対する脈動成分の比（Ｒ）のデータ列を、
【００２９】
【数２】
【００３０】
とする。
同じ時刻tiに得られたＩＲとＲとのデータを次のように行列で定義する。
すなわち、
【００３１】
【数３】
【００３２】
また、θ[rad]回転させる回転行列をＡとすると、Ａは次のように表すことができる。
【００３３】
【数４】
【００３４】
そうすると、ＳをＡによりθ[rad]回転させることにより次のＸが得られる。
【００３５】
【数５】
【００３６】
なお、回転行列Ａは、上記のほか、
【００３７】
【数６】
【００３８】
を用いてもよい。
ここで、θを０〜９π／３０[rad] までπ／３０[rad]ずつ脈波データＳを回転させて
得られるグラフを図４に示す。
図４からわかるように、横軸ゼロ、縦軸ゼロの点（赤色光脈波と赤外光脈波との両方が
平均である点）を中心として回転されており、θが９π／３０[rad]のときに、横軸(X1)
へ射影した領域が最小になり、縦軸(X2)へ射影した領域が最大となっている。
θを９π／３０[rad]よりさらにπ／２[rad]回転させ２４π／３０[rad] （＝１２π／
１５[rad]）回転させた場合には横軸(X1)へ射影した領域が最大になり、縦軸(X2)へ射影
した領域が最小となることは明らかである。
【００３９】
次に、θを９π／３０[rad]、２４π／３０[rad]としたときの回転行列Ａにより、測定
された脈波データＳが処理されてＸとなった結果、どのような波形となるかを説明する。
図５は、図２に示した脈波データＳを、θを９π／３０[rad]として回転行列Ａにより
処理したＸの波形を示す。
横軸へ射影した領域が最小になったX1(t i)は、
【００４０】
【数７】
【００４１】
一方、横軸へ射影した領域が最大になったX2(t i)は、
【００４２】
【数８】
【００４３】
により演算される。
図５のX1の波形からはノイズが除去されたことがわかる。
なお、脈波データＳを、θを２４π／３０[rad]として回転行列Ａにより処理した場合
には、X2の波形がノイズが除去された波形となる。
横軸へ射影した領域が最大になるX1(t i)は、
【００４４】
【数９】
【００４５】
一方、横軸へ射影した領域が最小になるX2(t i)は、
【００４６】
【数１０】
【００４７】
により演算される。
このように横軸へ射影した領域が最小になるように回転角θを設定して、脈波データSを処理すれば、ノイズが抑制された脈波主成分波形を得ることができる。
【００４８】
次に、脈波の基本周波数の演算について説明する。
ノイズが低減される前の図２に示した脈波信号と、回転行列を用いてノイズが低減された
脈波主成分波形を周波数解析して得られたスペクトルをそれぞれ図６に示す。
横軸は周波数、縦軸はスペクトルである。
ノイズが低減される前の脈波(Before-rotation)信号のスペクトルは、ノイズの周波数帯
域ｆｎのスペクトルが強くでており、脈波信号の基本周波数ｆｓのスペクトルはほとんど
現れていない。
一方、回転行列を用いてノイズが低減された脈波主成分波形(After-rotation)を周波数
解析して得られたスペクトルでは、脈波信号の基本周波数ｆｓのスペクトルがノイズの周
波数帯域ｆｎのスペクトルと区別できるほど強く現れていることがわかり、脈波信号の基
本周波数ｆｓを求めることができる。
そして、脈波信号の基本周波数ｆｓ[Hz]が求まれば、脈拍数ｆｓ×６０[回／min]を容
易に求めることができる。
【００４９】
このように、所定角度の回転行列を用いることにより、ノイズが低減された脈波主成分
波形を得ることができ、脈波信号の基本周波数ないし脈拍数を求めることができる。
ここで、所定角度は、予め決められたものでもよく、測定期間中アダプティブに変化さ
せてもよい。
【００５０】
＜酸素飽和度の演算＞
図３は、上述のように図９に示される赤外光ＩＲのデータを横軸に、赤色光Ｒのデータ
を縦軸にとったグラフであるが、このグラフの傾きをノルム比を用いて求める。
まず、赤外光脈波データＩＲのL２ノルムを求める。
赤外光脈波データ列はIR = [ IR(ti) ： ti = 0, 1, 2, 3, ・・・]（１）であるから、L２ノルムは次の式で表すことができる。
【００５１】
【数１１】
【００５２】
次に、赤色光脈波データＲのL２ノルムを求める。
赤色光脈波データ列はR = [ R(ti) ： ti = 0, 1, 2, 3, ・・・]（２）であるから、L２ノルムは次の式で表すことができる。
【００５３】
【数１２】
【００５４】
そこで、
【００５５】
【数１３】
【００５６】
とすればΦは、酸素飽和度SｐO2と相関するので、その相関を表す関数をｆとすれば、
【００５７】
【数１４】
【００５８】
と表され、酸素飽和度ＳｐＯ2を求めることができる。
なお、ノルム比を傾きとした直線を図３に示す。
ノルムとは、数学的概念の１つで、ユークリッドノルム（Euclidean-norm）又は２乗ノ
ルムは、ｎ個の要素を持つベクトルの大きさをスカラ量に写像するものである。
このように、所定期間の赤色光脈波データＲのL２ノルム値（２乗ノルム）と赤外光脈波データのL２ノルム値の比に基づいて、酸素飽和度SｐＯ2を求めることができる。
ここで、所定期間は逐次得られる現在の脈波から過去にさかのぼって所定期間分の赤色
光脈波データＲ、赤外光脈波データＩＲを用いるとよい。
また、ノルム値として、L２ノルムを用いたが、他の演算方法によるノルム値を用いて
もよい。
【００５９】
また、酸素飽和度の演算に関しては、
脈波信号に対して、ノイズ信号が比較的小なる場合は上記ノルム比を用いて演算しても良
いが、比較的大なる場合は、上記のノルム比を用いて求める方法とは別に、前記特願２０
０１−３３２３８３号に記載した周波数解析を用いて求めた基本周波数に替えて、前記回
転による処理により求めた基本周波数を用いて演算することもできる。
【００６０】
（第１の実施例）
次に、上記原理を用いた装置を、概略構成ブロック図と処理フローにより説明する。
概略構成ブロック図は先に説明した図１と同じである。
発光素子１、２から、交互に発光するように駆動回路３により駆動されることにより、
異なる波長の光が発光される。
これらの発光素子１、２からの発光は生体組織４を透過して受光部（フォトダイオード
）５で受光され、電気信号に変換される。
そして、これらの変換された信号は増幅器６で増幅され、マルチプレクサ７によりそれ
ぞれの光波長に対応したフィルタ８−１、８−２に振り分けられる。
各フィルターに振り分けられた信号はフィルタ８−１、８−２によりフィルタリングさ
れてノイズ成分が低減され、A／Ｄ変換器９によりデジタル化される。
デジタル化された赤外光、赤色光に対応する各信号列が、それぞれの脈波を形成してい
る。
デジタル化された各信号列は処理部１０に入力され、ＲＯＭ１２に格納されているプロ
グラムにより処理され、脈拍数ＰＲ、酸素飽和度ＳｐＯ2が演算され、その値が表示部１
１に表示される。
【００６１】
次に、脈拍数ＰＲ、酸素飽和度ＳｐＯ2を演算する処理フローを図７を用いて説明する。
測定が開始される（ステップＳ１）と、上記のように赤色光脈波、赤外光脈波が検出され
（ステップＳ２）、デジタル化された各信号列（各脈波データ）が処理部１０に取り込ま
れる。
処理部１０では、ＲＯＭ１２に格納されているプログラムにより、処理過程のデータを
ＲＡＭ１３に読み書きしながら、各脈波データを次のように処理する。
【００６２】
先ず、赤外光脈波、赤色光脈波それぞれの脈波の直流成分に対する脈動成分の比を脈波ご
とに求める。（ステップＳ３）
次に、脈拍数ＰＲを求める処理（ステップＳ４〜Ｓ６）と酸素飽和度ＳｐＯ2を求める
処理（ステップＳ７〜Ｓ９）が同時に行われる。
【００６３】
脈拍数ＰＲを求める処理（ステップＳ４〜Ｓ６）では、
予め回転角が設定された回転行列Ａにより、赤外光脈波データＩＲと赤色光脈波データ
ＲとのデータＳから、式５よりノイズが低減された波形を得る。（ステップＳ４）
設定する回転角は、図３に示すような赤外光脈波データＩＲを横軸に赤色光脈波データ
Ｒを縦軸にとったグラフを回転すると図４に示すように軸方向へ射影した領域が最小とな
るような角度である。
回転角度は例えば、９π／３０[rad]あるいは２４π／３０[rad]がよい。
ノイズが低減された波形は、射影領域が最小となる軸成分のデータより得ることができ
る。
そして、ノイズが低減された波形を図６に示すように周波数解析を行い、脈波データの基
本周波数を求める。（ステップＳ５）
そして、その基本周波数から脈拍数ｆｓをｆＳ×６０[回／min]から求め、表示部１１
に表示する。
【００６４】
酸素飽和度ＳｐＯ２を求める処理（ステップＳ７〜Ｓ９）では、所定期間分の赤外光脈
波データＩＲと赤色光脈波データＲとからそれぞれのＬ２ノルム値を式１０、式１１から
求め、さらにそれぞれのＬ２ノルム値の比を式１２から求める。
次に、ノイズの除去された赤外光と赤色光の脈波信号の比を求め、酸素飽和度を演算する。
（ステップＳ７）
そのＬ２ノルム比をΦとして、式１２より酸素飽和度ＳｐＯ２を求め（ステップＳ８）
、表示部１１に表示する。（ステップＳ９）
【００６５】
測定を継続するときはステップＳ２に戻り処理を繰り返し、測定を計測しない場合は測
定を終了する。（ステップＳ１１）
【００６６】
（第２の実施例）
次に、別の第２の実施例を図８を用いて説明する。
第２の実施例が第１の実施例と相違する点はステップＳ４において、回転角は予め定め
られたものではなく、得られるデータから回転角を求める点であり、図８に示すようにス
テップＳ４−１とステップＳ４−２に分けて処理する。
他のステップは第１の実施例と同様なので、説明は省略する。
【００６７】
脈拍数ＰＲを求める処理（ステップＳ４−１〜Ｓ６）では、
先ず、所定期間分の赤外光脈波データＩＲと赤色光脈波データＲとのデータを用い、図
３に示すようなグラフを描く。
そして、いかなる回転角度であれば軸方向へ射影した領域が最小になるかを求める。（
ステップＳ４−１）
次に、求められた回転角度による回転行列により、各波長の脈波データを処理し、射影
領域が最小となる軸の成分のデータからノイズが低減された波形を得る。（ステップＳ４
−２）。
このように、第２の実施例の特徴は、回転行列の回転角度が、固定化された角度ではな
く、検出される脈波データに応じて適時変更されるアダプティブ性を有する点にある。
【００６８】
（第３実施例）
次に、脈拍数ＰＲ、酸素飽和度ＳｐＯ2を、周波数解析を用いて求めた基本周波数に替
えて、前記回転による処理により求めた基本周波数を用いて演算する処理フローを図１０
を用いて説明する。
測定が開始される（ステップＳ１）と、上記のように赤色光脈波、赤外光脈波が検出さ
れ（ステップＳ２）、デジタル化された各信号列（各脈波データ）が処理部１０に取り込
まれる。
処理部１０では、ＲＯＭ１２に格納されているプログラムにより、処理過程のデータを
ＲＡＭ１３に読み書きしながら、各脈波データを次のように処理する。
【００６９】
先ず、赤外光脈波、赤色光脈波それぞれの脈波の直流成分に対する脈動成分の比を脈波
ごとに求める。（ステップＳ３）
次に、脈拍数ＰＲを求める処理（ステップＳ４〜Ｓ６）と酸素飽和度ＳｐＯ2を求める
処理（ステップＳ７〜Ｓ９）が同時に行われる。
【００７０】
脈拍数ＰＲを求める処理（ステップＳ４〜Ｓ６）では、
予め回転角が設定された回転行列Ａにより、赤外光脈波データＩＲと赤色光脈波データ
ＲとのデータＳから、式５よりノイズが低減された波形を得る。（ステップＳ４）
設定する回転角は、図３に示すような赤外光脈波データＩＲを横軸に赤色光脈波データ
Ｒを縦軸にとったグラフを回転すると図４に示すように軸方向へ射影した領域が最小とな
るような角度である。
回転角度は例えば、９π／３０[rad]あるいは２４π／３０[rad]がよい。
ノイズが低減された波形は、射影領域が最小となる軸成分のデータより得ることができ
る。
そして、ノイズが低減された波形を図６に示すように周波数解析を行い、脈波データの
基本周波数を求める。（ステップＳ５）
更に、その基本周波数から脈拍数ｆｓをｆＳ×６０[回／min]から求め、表示部１１に
表示する。
【００７１】
また、酸素飽和度ＳｐＯ２を求める処理（ステップＳ７〜Ｓ９）では、
赤外光及び赤色光の脈波信号を前記基本周波数で、または基本周波数とその高調波で構
成したフィルタを通し、ノイズの除去された信号を求める。（ステップＳ７）
次に、ノイズの除去された赤外光と赤色光の脈波信号の比を求め、酸素飽和度を演算し
（ステップＳ８）表示部１１に表示する。（ステップＳ９）
【００７２】
測定を継続するとき（ステップＳ１０）は、ステップＳ２に戻り処理を繰り返し、測定
を計測しない場合は測定を終了する。（ステップＳ１１）
【００７３】
（第４の実施例）
更に、別の第４の実施例を図１１を用いて説明する。
第４の実施例が第３の実施例と相違する点はステップＳ４において、回転角は予め定めら
れたものではなく、得られるデータから回転角を求める点であり、図１１に示すようにス
テップＳ４−１とステップＳ４−２に分けて処理する。
他のステップは第３の実施例と同様なので、説明は省略する。
【００７４】
脈拍数ＰＲを求める処理（ステップＳ４−１〜Ｓ６）では、
先ず、所定期間分の赤外光脈波データＩＲと赤色光脈波データＲとのデータを用い、図
３に示すようなグラフを描く。
そして、いかなる回転角度であれば軸方向へ射影した領域が最小になるかを求める。（
ステップＳ４−１）
次に、求められた回転角度による回転行列により、各波長の脈波データを処理し、射影
領域が最小となる軸の成分のデータからノイズが低減された波形を得る。（ステップＳ４
−２）。
このように、第４の実施例の特徴は、回転行列の回転角度が、固定化された角度ではな
く、検出される脈波データに応じて適時変更されるアダプティブ性を有する点にある。
【００７５】
以上は、動脈血酸素飽和度を測定するパルスオキシメータを例に挙げて説明したが、本
発明の技術はパルスオキシメータに限られず、特殊ヘモグロビン（一酸化炭素ヘモグロビ
ン、Ｍｅｔヘモグロビンなど）、血中に注入された色素などの血中吸光物質をパルスフォ
トメトリーの原理を用いて測定する装置(パルスフォトメトリー)にも光源の波長を選択す
ることで適用できる。
なお、本明細書中において、生体組織（媒体）を介して測定された信号を、２次元直交
座標に展開する前までは、「脈波信号」と記載し、２次元直交座標に展開中は、「脈波デ
ータ」と記載し、求められた波形は「脈波主成分波形」と区別して記載している。
【産業上の利用可能性】
【００７６】
本発明では、同一の媒体からほぼ同時に抽出される２つの同種の信号
を処理して共通の信号成分を抽出する計算処理負担を軽減した信号処理が実現できる。
また、本発明では、前記媒体の体動によるノイズが脈波信号に生じ
た場合であっても、対象物質の濃度を精度よく測定できる。
また、体動によるノイズが脈波データ信号に生じた場合であっても、脈波信号からノイ
ズを除去し、精度よく脈拍や血中の吸光物質の濃度を求めることができるので、産業上の
利用可能性は極めて大きい。
【符号の説明】
【００７７】
１ 発光素子
２ 発光素子
３ 駆動回路
４ 生体組織
５ フォトダイオード
６ 変換器
７ マルチプレクサ
８ フィルタ
９ Ａ／Ｄ変換器
１０ 処理部
１１ 表示部
１２ ＲＯＭ
１３ ＲＡＭ
【技術分野】
【０００１】
本発明は、一つの媒体からほぼ同時に抽出される２つの同種の信号を処理して共通の信
号成分を抽出する信号処理に関し、特には医療の分野において、特に循環器系の診断に用
いられるパルスフォトメータにおける信号処理の改良に関する。
【背景技術】
【０００２】
一つの媒体からほぼ同時に抽出された２つの信号から信号成分と雑音成分に分離する方
法には様々な方法が提案されている。
それらは、一般的には周波数領域や時間領域による処理が行われている。
医療現場でも、光電脈波計と言われる脈波波形や脈拍数を測定する装置、血液に含まれ
る吸光物質の濃度測定として、酸素飽和度ＳｐＯ2の測定装置、一酸化炭素ヘモグロビン
やＭｅｔヘモグロビン等の特殊ヘモグロビンの濃度の測定装置、注入色素濃度の測定装置
などがパルスフォトメータとして知られている。
中でも酸素飽和度ＳｐＯ2の測定装置を特にパルスオキシメータと呼んでいる。
【０００３】
パルスフォトメータの原理は、対象物質への吸光性が異なる複数の波長の光を生体組織
に透過又は反射させ、その透過光又は反射光の光量を連続的に測定することで得られる脈
波データ信号から対象物質の濃度を求めるものである。
そしてその脈波データに雑音が混入すると、正しい濃度の計算が出来ず、誤処置につな
がる危険が生じる。
パルスフォトメータにおいても従来より雑音を低減するために周波数帯域を分割して信
号成分に着目したり、２つの信号の相関を取るなどの方法が提案されてきた。
しかし、これらの方法は解析に時間がかかるなどの問題があった。
【０００４】
そこで、本出願人は、特許第３２７０９１７号（特許文献１）において、異なる２つの
波長の光を生体組織に照射して透過光から得られる２つの脈波信号のそれぞれの大きさを
縦軸、横軸としてグラフを描き、その回帰直線を求め、その回帰直線の傾きに基づいて、
動脈血中の酸素飽和度ないし吸光物質濃度を求めることを提案している。
この発明により、測定精度を高め、低消費電力化することができた。
しかし、各波長の脈波信号についての多くのサンプリングデータを用いて回帰直線ない
しその傾きを求めるためには、なお多くの計算処理を要していた。
【０００５】
更に本出願人は、特願２００１−３３２３８３号（特許文献２）においては、周波数解
析を用いてはいるが、その解析においては従来技術のように脈波信号そのものを抽出する
のではなく、脈波信号の基本周波数を求め、さらには精度を高めるためにその高調波周波
数を用いたフィルタを用いて脈波信号をフィルタリングする方法を提案している。
しかし、基本周波数を求める点に関しては更なる改善が望まれていた。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第３２７０９１７号 （請求項１、２、図２、図４）
【特許文献２】特願２００３−１３５４３４号（０００６）
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明の課題（目的）は、同一の媒体からほぼ同時に抽出される２つの同種の信号を処
理して共通の信号成分を抽出する計算処理負担を軽減したパルスフォトメータを提供することにある。
また、上記信号処理方法を適用して、前記媒体の体動によるノイズが脈波信号に生じた
場合であっても、対象物質の濃度を精度よく測定することにある。
また、体動によるノイズが脈波データ信号に生じた場合であっても、脈波信号からノイ
ズを除去し、精度よく脈拍を求めることにある。
【課題を解決するための手段】
【０００８】
・前記課題を解決するために、本発明のパルスフォトメータは、異なる２つの波長の光を同一の生体組織に照射する発光手段と、前記発光手段から発生し前記同一の生体組織を透過または反射した各波長の光を電気信号に変換する受光手段と、前記受光手段で変換された前記電気信号より得られた２つの波長の離散的時系列脈波データをそれぞれ縦軸または横軸とする２次元直交座標に展開した、各波長の脈波データそれぞれについてノルム値を求め、さらにそのノルム値の比を求めるノルム比計算手段と、前記ノルム比計算手段により計算されたノルム比に基づいて、血中の酸素の濃度を求める血中酸素濃度演算手段とを具備することを特徴とする。（請求項１）
【０００９】
また、参考例では、
・同一の媒体から抽出される所定期間の２つの信号を２次元直交座標上に展開するステップと、前記展開された信号を、回転行列を用いて、あらかじめ決められた角度または所定条件に基づいて決められた角度に回転させるステップとにより前記２つの信号を処理する。
この信号処理方法によって、同一の媒体からほぼ同時に抽出される２つの同種の信号を
処理して共通の信号成分を抽出する計算処理負担を軽減することができる。
・前記信号処理の結果、信号主成分の基本周波数を求める。
この構成により、ノイズが低減された基本波波形を容易に取得することができる。
・同一の媒体からの２つの信号を入力する入力手段と、前記入力手段から入力された所定
期間分の２つの信号を２次元直交座標上に展開する展開手段と、該展開されたデータをあ
らかじめ決められた角度または所定条件に基づいて決められた角度に回転させる回転行列
を用いて処理する処理手段とを具えた構成とする。
この構成により、回転行列を用いて処理することでノイズが低減された波形を容易に取
得することができる。
【００１０】
・異なる２つの波長の光を生体組織に照射する発光手段と、前記発光手段から発生し前記
生体組織を透過または反射した各波長の光を電気信号に変換する受光手段とを具えたパル
スフォトメータにおいて、
前記各波長の電気信号より得られた脈波データを、あらかじめ決められた角度にまたは
所定条件に基づいて決められた角度に回転させる回転行列を用いて処理する処理手段とを
具備することを特徴とする。
この構成により、脈波データを回転行列を用いて処理することでノイズが低減された波
形を容易に取得することができる。
・前記脈波データが、ほぼ同時点に受光した２つの波長の受光信号から得られる電気信号
を、前記２つの波長をそれぞれ縦軸または横軸とする２次元直交座標に展開したものであ
ることを特徴とする。
・前記脈波データは所定時間分であり、かつ経時移動して処理されることを特徴とする。
・異なる２つの波長の光を生体組織に照射する発光手段と、前記発光手段から発生し前記
生体組織を透過または反射した各波長の光を電気信号に変換する受光手段とを具えたパル
スフォトメータにおいて、前記各波長の電気信号より得られた脈波信号を、各波長をそれ
ぞれ縦軸または横軸とする２次元直交座標に展開し、縦軸または横軸に射影される領域が
最大または最小となるいずれかの条件を満足するように回転行列を用いて回転させて処理
する第２の処理手段とを具備することを特徴とする。
【００１１】
・異なる２つの波長の光を生体組織に照射する発光手段と、前記発光手段で発生し前記生
体組織を透過または反射した各波長の光を電気信号に変換する受光手段と、前記受光部に
より得られる前記２つの波長の脈波信号をそれぞれ縦軸または横軸とする２次元直交座標
に展開した、各波長の各波長の脈波データを、所定角度に回転させる回転行列を用いて処
理した波形を求める波形取得手段と、前記波形取得手段により得られた波形の周波数解析
により、脈波の基本周波数または脈拍数を求める波形解析手段と、前記波形解析手段の出
力から血中の吸光物質の濃度を求める血中吸光物質濃度演算手段とを具備するパルスフォ
トメータ。
・前記血中吸光物質濃度演算手段は、動脈血中の酸素飽和度、特殊ヘモグロビン濃度、ま
たは注入色素濃度のうち少なくとも１つを演算する。
【発明の効果】
【００１２】
請求項１に記載のパルスフォトメータによれば、各波長の脈波データから得られるノル
ム比から精度よく血中の酸素の濃度を求めることができる。
また、パルスフォトメータによれば、動脈血中の酸素飽和度、特殊ヘモグロビン濃度、または注入色素濃度のうち少なくとも１つを演算することができる。
【００１３】
参考例のパルスフォトメータによれば、脈波データを回転行列を用いて処理することでノイズが低減された波形を取得することができる。
【００１４】
参考例のパルスフォトメータによれば、脈波データを回転行列を用いて処理することでノイズが低減された波形を取得し、その波形から精度よく脈波データの基本周波数または脈拍数を求めることができる。
【００１５】
参考例のパルスフォトメータによれば、脈波データを回転行列を用いてノイズを除去するのに適した回転角度を決定でき、その決定した回転角度の回転行列によって処理することでノイズが低減された波形を取得することができる。
【００１６】
参考例のパルスフォトメータによれば、脈波データを回転行列を用いてノイズを除去するのに適した回転角度を決定でき、その決定した回転角度の回転行列によって処理することでノイズが低減された波形を取得し、その波形から精度よく脈波データの基本周波数または脈拍数を求めることができる。
【００１７】
参考例のパルスフォトメータによれば、回転行列の回転角の決定を容易に決定できる。
【００１８】
参考例のパルスフォトメータによれば、脈波データを回転行列を用いて処理することでノイズが低減された波形による精度の良い脈波の基本周波数または脈拍数の取得と、動脈血中の酸素飽和度、特殊ヘモグロビン濃度、または注入色素濃度のうち少なくとも１つの演算とが同一の装置でほぼ同時に処理して出力することができる。
参考例のパルスフォトメータによれば、動脈血中の酸素飽和度、特殊ヘモグロビン濃度、または注入色素濃度のうち少なくとも１つを演算することができる。
【図面の簡単な説明】
【００１９】
【図１】本発明の概略構成を示すブロック図である。
【図２】検出された脈波を示す図である。
【図３】図９に示される赤外光ＩＲのデータを横軸に、赤色光Ｒのデータを縦軸にと ったグラフである。
【図４】図３のグラフをπ／３０[rad]ずつ回転させた図である。
【図５】回転角度９π／３０[rad]の回転行列により処理された脈波の波形を示す図 である。
【図６】図５に示すX1の波形のスペクトルを示す図である。
【図７】第１の実施例における処理フローを示すフローチャートである。
【図８】第２の実施例における処理フローを示すフローチャートである。
【図９】血液中の吸光物質の吸光度の変動の測定原理を説明する波形図である。
【図１０】第３の実施例における処理フローを示すフローチャートである。
【図１１】第４の実施例における処理フローを示すフローチャートである。
【発明を実施するための形態】
【００２０】
本発明の実施の形態を説明するにあたり、動脈血酸素飽和度を測定するパルスオキシメ
ータを例に挙げて原理を説明する。
なお、本発明の技術は、パルスオキシメータに限られず、特殊ヘモグロビン（一酸化炭
素ヘモグロビン、Ｍｅｔヘモグロビンなど）、血中に注入された色素などの血中吸光物質
をパルスフォトメトリーの原理を用いて測定する装置(パルスフォトメータ)に適用できる。
【００２１】
動脈血酸素飽和度を測定するパルスオキシメータの構成は、概略構成ブロック図である
図１のようになっている。
異なる波長の光を発光する発光素子１、２は、交互に発光するように駆動回路３により
駆動される。
発光素子１、２に採用する光はそれぞれ動脈血酸素飽和度による影響が少ない赤外光（
例えば９４０[nm]）、動脈血酸素飽和度の変化に対する感度が高い赤色光（例えば６６０
[nm]）がよい。
【００２２】
これらの発光素子１、２からの発光は生体組織４を透過してフォトダイオード５で受光
して電気信号に変換される。
なお、反射光を受光するようにしてもよい。
そして、これらの変換された信号は増幅器６で増幅され、マルチプレクサ７によりそれ
ぞれの光波長に対応したフィルタ８−１、８−２に振り分けられる。
各フィルターに振り分けられた信号はフィルタ８−１、８−２によりフィルタリングさ
れてノイズ成分が低減され、Ａ／Ｄ変換器９によりデジタル化される。
【００２３】
デジタル化された赤外光、赤色光に対応する各信号列が、それぞれの脈波信号を形成し
ている。
デジタル化された各信号列は処理部１０に入力され、ＲＯＭ１２に格納されているプロ
グラムにより処理され、酸素飽和度ＳｐＯ2が測定され、その値が表示部１１に表示され
る。
【００２４】
＜回転行列によるノイズ低減と脈波の基本周波数の演算＞
先ず、血液中の吸光物質の吸光度（減光度）の変動の測定について説明する。
図９(a)及び(b)は、前記発光素子１、２からの発光された光が生体組織４を透過してフ
ォトダイオード５で受光して電気信号に変換された脈波データで、(a)は赤色光の場合を
、(b)は赤外光を示している。
図９の(a)では、横軸を時間、縦軸を受光出力とすると、フォトダイオード５での受光
出力は、赤色光の直流成分（Ｒ’）と脈動成分（ΔＲ’）が重畳された波形となっている。
また、図９の(b)では、横軸を時間、縦軸を受光出力とすると、フォトダイオード５で
の受光出力は、赤外光の直流成分（ＩＲ’）と脈動成分（ΔＩＲ’）が重畳された波形となっている。
図２は、図９に示すような脈波において、８秒間分の、直流成分（Ｒ’、ＩＲ’）に対
する脈動成分（ΔＲ’、ΔＩＲ’）の比（ＩＲ＝ΔＩＲ’／ＩＲ’）をとり、さらにその８秒間分のデータの平均値をゼロに合わせたものである。
なお、図２の如き、平均値をゼロとする処理を行わなくとも演算は可能である。
図３は、図９に示される赤外光ＩＲのデータを横軸に、赤色光Ｒのデータを縦軸にとっ
たグラフである。
【００２５】
次に、Ａ／Ｄ変換器９によってデジタル化した各波長の２つの脈波データ信号を回転行
列を用いてノイズを低減する演算処理について説明する。
なお、赤外光と赤色光とは交互に発光されるため厳密には同時に発光されるものではな
いが、隣り合う得られた赤外光受光値と赤色光受光値を同時刻に得られたものとして扱い
、所定時間分の赤外光の脈波信号と赤色光の脈波信号を２次元直交座標上に展開する。
すなわち図３のグラフを作成している。
また、脈波の直流成分に対する脈動成分の比をとることで脈拍による吸光度の脈動分が
近似される。
図３のグラフに見られる推移は４５度になっていないが、その理由は、赤外光脈波の脈
動成分の振幅と赤色光脈波の脈動成分の振幅とに差があるため、およびノイズが重畳して
いるためである。
【００２６】
次に、展開された脈波データに回転行列を用いて回転演算を施すこととする。
赤外光脈波の直流成分に対する脈動成分の比（ＩＲ）のデータ列を、
【００２７】
【数１】
【００２８】
赤色光脈波の直流成分に対する脈動成分の比（Ｒ）のデータ列を、
【００２９】
【数２】
【００３０】
とする。
同じ時刻tiに得られたＩＲとＲとのデータを次のように行列で定義する。
すなわち、
【００３１】
【数３】
【００３２】
また、θ[rad]回転させる回転行列をＡとすると、Ａは次のように表すことができる。
【００３３】
【数４】
【００３４】
そうすると、ＳをＡによりθ[rad]回転させることにより次のＸが得られる。
【００３５】
【数５】
【００３６】
なお、回転行列Ａは、上記のほか、
【００３７】
【数６】
【００３８】
を用いてもよい。
ここで、θを０〜９π／３０[rad] までπ／３０[rad]ずつ脈波データＳを回転させて
得られるグラフを図４に示す。
図４からわかるように、横軸ゼロ、縦軸ゼロの点（赤色光脈波と赤外光脈波との両方が
平均である点）を中心として回転されており、θが９π／３０[rad]のときに、横軸(X1)
へ射影した領域が最小になり、縦軸(X2)へ射影した領域が最大となっている。
θを９π／３０[rad]よりさらにπ／２[rad]回転させ２４π／３０[rad] （＝１２π／
１５[rad]）回転させた場合には横軸(X1)へ射影した領域が最大になり、縦軸(X2)へ射影
した領域が最小となることは明らかである。
【００３９】
次に、θを９π／３０[rad]、２４π／３０[rad]としたときの回転行列Ａにより、測定
された脈波データＳが処理されてＸとなった結果、どのような波形となるかを説明する。
図５は、図２に示した脈波データＳを、θを９π／３０[rad]として回転行列Ａにより
処理したＸの波形を示す。
横軸へ射影した領域が最小になったX1(t i)は、
【００４０】
【数７】
【００４１】
一方、横軸へ射影した領域が最大になったX2(t i)は、
【００４２】
【数８】
【００４３】
により演算される。
図５のX1の波形からはノイズが除去されたことがわかる。
なお、脈波データＳを、θを２４π／３０[rad]として回転行列Ａにより処理した場合
には、X2の波形がノイズが除去された波形となる。
横軸へ射影した領域が最大になるX1(t i)は、
【００４４】
【数９】
【００４５】
一方、横軸へ射影した領域が最小になるX2(t i)は、
【００４６】
【数１０】
【００４７】
により演算される。
このように横軸へ射影した領域が最小になるように回転角θを設定して、脈波データSを処理すれば、ノイズが抑制された脈波主成分波形を得ることができる。
【００４８】
次に、脈波の基本周波数の演算について説明する。
ノイズが低減される前の図２に示した脈波信号と、回転行列を用いてノイズが低減された
脈波主成分波形を周波数解析して得られたスペクトルをそれぞれ図６に示す。
横軸は周波数、縦軸はスペクトルである。
ノイズが低減される前の脈波(Before-rotation)信号のスペクトルは、ノイズの周波数帯
域ｆｎのスペクトルが強くでており、脈波信号の基本周波数ｆｓのスペクトルはほとんど
現れていない。
一方、回転行列を用いてノイズが低減された脈波主成分波形(After-rotation)を周波数
解析して得られたスペクトルでは、脈波信号の基本周波数ｆｓのスペクトルがノイズの周
波数帯域ｆｎのスペクトルと区別できるほど強く現れていることがわかり、脈波信号の基
本周波数ｆｓを求めることができる。
そして、脈波信号の基本周波数ｆｓ[Hz]が求まれば、脈拍数ｆｓ×６０[回／min]を容
易に求めることができる。
【００４９】
このように、所定角度の回転行列を用いることにより、ノイズが低減された脈波主成分
波形を得ることができ、脈波信号の基本周波数ないし脈拍数を求めることができる。
ここで、所定角度は、予め決められたものでもよく、測定期間中アダプティブに変化さ
せてもよい。
【００５０】
＜酸素飽和度の演算＞
図３は、上述のように図９に示される赤外光ＩＲのデータを横軸に、赤色光Ｒのデータ
を縦軸にとったグラフであるが、このグラフの傾きをノルム比を用いて求める。
まず、赤外光脈波データＩＲのL２ノルムを求める。
赤外光脈波データ列はIR = [ IR(ti) ： ti = 0, 1, 2, 3, ・・・]（１）であるから、L２ノルムは次の式で表すことができる。
【００５１】
【数１１】
【００５２】
次に、赤色光脈波データＲのL２ノルムを求める。
赤色光脈波データ列はR = [ R(ti) ： ti = 0, 1, 2, 3, ・・・]（２）であるから、L２ノルムは次の式で表すことができる。
【００５３】
【数１２】
【００５４】
そこで、
【００５５】
【数１３】
【００５６】
とすればΦは、酸素飽和度SｐO2と相関するので、その相関を表す関数をｆとすれば、
【００５７】
【数１４】
【００５８】
と表され、酸素飽和度ＳｐＯ2を求めることができる。
なお、ノルム比を傾きとした直線を図３に示す。
ノルムとは、数学的概念の１つで、ユークリッドノルム（Euclidean-norm）又は２乗ノ
ルムは、ｎ個の要素を持つベクトルの大きさをスカラ量に写像するものである。
このように、所定期間の赤色光脈波データＲのL２ノルム値（２乗ノルム）と赤外光脈波データのL２ノルム値の比に基づいて、酸素飽和度SｐＯ2を求めることができる。
ここで、所定期間は逐次得られる現在の脈波から過去にさかのぼって所定期間分の赤色
光脈波データＲ、赤外光脈波データＩＲを用いるとよい。
また、ノルム値として、L２ノルムを用いたが、他の演算方法によるノルム値を用いて
もよい。
【００５９】
また、酸素飽和度の演算に関しては、
脈波信号に対して、ノイズ信号が比較的小なる場合は上記ノルム比を用いて演算しても良
いが、比較的大なる場合は、上記のノルム比を用いて求める方法とは別に、前記特願２０
０１−３３２３８３号に記載した周波数解析を用いて求めた基本周波数に替えて、前記回
転による処理により求めた基本周波数を用いて演算することもできる。
【００６０】
（第１の実施例）
次に、上記原理を用いた装置を、概略構成ブロック図と処理フローにより説明する。
概略構成ブロック図は先に説明した図１と同じである。
発光素子１、２から、交互に発光するように駆動回路３により駆動されることにより、
異なる波長の光が発光される。
これらの発光素子１、２からの発光は生体組織４を透過して受光部（フォトダイオード
）５で受光され、電気信号に変換される。
そして、これらの変換された信号は増幅器６で増幅され、マルチプレクサ７によりそれ
ぞれの光波長に対応したフィルタ８−１、８−２に振り分けられる。
各フィルターに振り分けられた信号はフィルタ８−１、８−２によりフィルタリングさ
れてノイズ成分が低減され、A／Ｄ変換器９によりデジタル化される。
デジタル化された赤外光、赤色光に対応する各信号列が、それぞれの脈波を形成してい
る。
デジタル化された各信号列は処理部１０に入力され、ＲＯＭ１２に格納されているプロ
グラムにより処理され、脈拍数ＰＲ、酸素飽和度ＳｐＯ2が演算され、その値が表示部１
１に表示される。
【００６１】
次に、脈拍数ＰＲ、酸素飽和度ＳｐＯ2を演算する処理フローを図７を用いて説明する。
測定が開始される（ステップＳ１）と、上記のように赤色光脈波、赤外光脈波が検出され
（ステップＳ２）、デジタル化された各信号列（各脈波データ）が処理部１０に取り込ま
れる。
処理部１０では、ＲＯＭ１２に格納されているプログラムにより、処理過程のデータを
ＲＡＭ１３に読み書きしながら、各脈波データを次のように処理する。
【００６２】
先ず、赤外光脈波、赤色光脈波それぞれの脈波の直流成分に対する脈動成分の比を脈波ご
とに求める。（ステップＳ３）
次に、脈拍数ＰＲを求める処理（ステップＳ４〜Ｓ６）と酸素飽和度ＳｐＯ2を求める
処理（ステップＳ７〜Ｓ９）が同時に行われる。
【００６３】
脈拍数ＰＲを求める処理（ステップＳ４〜Ｓ６）では、
予め回転角が設定された回転行列Ａにより、赤外光脈波データＩＲと赤色光脈波データ
ＲとのデータＳから、式５よりノイズが低減された波形を得る。（ステップＳ４）
設定する回転角は、図３に示すような赤外光脈波データＩＲを横軸に赤色光脈波データ
Ｒを縦軸にとったグラフを回転すると図４に示すように軸方向へ射影した領域が最小とな
るような角度である。
回転角度は例えば、９π／３０[rad]あるいは２４π／３０[rad]がよい。
ノイズが低減された波形は、射影領域が最小となる軸成分のデータより得ることができ
る。
そして、ノイズが低減された波形を図６に示すように周波数解析を行い、脈波データの基
本周波数を求める。（ステップＳ５）
そして、その基本周波数から脈拍数ｆｓをｆＳ×６０[回／min]から求め、表示部１１
に表示する。
【００６４】
酸素飽和度ＳｐＯ２を求める処理（ステップＳ７〜Ｓ９）では、所定期間分の赤外光脈
波データＩＲと赤色光脈波データＲとからそれぞれのＬ２ノルム値を式１０、式１１から
求め、さらにそれぞれのＬ２ノルム値の比を式１２から求める。
次に、ノイズの除去された赤外光と赤色光の脈波信号の比を求め、酸素飽和度を演算する。
（ステップＳ７）
そのＬ２ノルム比をΦとして、式１２より酸素飽和度ＳｐＯ２を求め（ステップＳ８）
、表示部１１に表示する。（ステップＳ９）
【００６５】
測定を継続するときはステップＳ２に戻り処理を繰り返し、測定を計測しない場合は測
定を終了する。（ステップＳ１１）
【００６６】
（第２の実施例）
次に、別の第２の実施例を図８を用いて説明する。
第２の実施例が第１の実施例と相違する点はステップＳ４において、回転角は予め定め
られたものではなく、得られるデータから回転角を求める点であり、図８に示すようにス
テップＳ４−１とステップＳ４−２に分けて処理する。
他のステップは第１の実施例と同様なので、説明は省略する。
【００６７】
脈拍数ＰＲを求める処理（ステップＳ４−１〜Ｓ６）では、
先ず、所定期間分の赤外光脈波データＩＲと赤色光脈波データＲとのデータを用い、図
３に示すようなグラフを描く。
そして、いかなる回転角度であれば軸方向へ射影した領域が最小になるかを求める。（
ステップＳ４−１）
次に、求められた回転角度による回転行列により、各波長の脈波データを処理し、射影
領域が最小となる軸の成分のデータからノイズが低減された波形を得る。（ステップＳ４
−２）。
このように、第２の実施例の特徴は、回転行列の回転角度が、固定化された角度ではな
く、検出される脈波データに応じて適時変更されるアダプティブ性を有する点にある。
【００６８】
（第３実施例）
次に、脈拍数ＰＲ、酸素飽和度ＳｐＯ2を、周波数解析を用いて求めた基本周波数に替
えて、前記回転による処理により求めた基本周波数を用いて演算する処理フローを図１０
を用いて説明する。
測定が開始される（ステップＳ１）と、上記のように赤色光脈波、赤外光脈波が検出さ
れ（ステップＳ２）、デジタル化された各信号列（各脈波データ）が処理部１０に取り込
まれる。
処理部１０では、ＲＯＭ１２に格納されているプログラムにより、処理過程のデータを
ＲＡＭ１３に読み書きしながら、各脈波データを次のように処理する。
【００６９】
先ず、赤外光脈波、赤色光脈波それぞれの脈波の直流成分に対する脈動成分の比を脈波
ごとに求める。（ステップＳ３）
次に、脈拍数ＰＲを求める処理（ステップＳ４〜Ｓ６）と酸素飽和度ＳｐＯ2を求める
処理（ステップＳ７〜Ｓ９）が同時に行われる。
【００７０】
脈拍数ＰＲを求める処理（ステップＳ４〜Ｓ６）では、
予め回転角が設定された回転行列Ａにより、赤外光脈波データＩＲと赤色光脈波データ
ＲとのデータＳから、式５よりノイズが低減された波形を得る。（ステップＳ４）
設定する回転角は、図３に示すような赤外光脈波データＩＲを横軸に赤色光脈波データ
Ｒを縦軸にとったグラフを回転すると図４に示すように軸方向へ射影した領域が最小とな
るような角度である。
回転角度は例えば、９π／３０[rad]あるいは２４π／３０[rad]がよい。
ノイズが低減された波形は、射影領域が最小となる軸成分のデータより得ることができ
る。
そして、ノイズが低減された波形を図６に示すように周波数解析を行い、脈波データの
基本周波数を求める。（ステップＳ５）
更に、その基本周波数から脈拍数ｆｓをｆＳ×６０[回／min]から求め、表示部１１に
表示する。
【００７１】
また、酸素飽和度ＳｐＯ２を求める処理（ステップＳ７〜Ｓ９）では、
赤外光及び赤色光の脈波信号を前記基本周波数で、または基本周波数とその高調波で構
成したフィルタを通し、ノイズの除去された信号を求める。（ステップＳ７）
次に、ノイズの除去された赤外光と赤色光の脈波信号の比を求め、酸素飽和度を演算し
（ステップＳ８）表示部１１に表示する。（ステップＳ９）
【００７２】
測定を継続するとき（ステップＳ１０）は、ステップＳ２に戻り処理を繰り返し、測定
を計測しない場合は測定を終了する。（ステップＳ１１）
【００７３】
（第４の実施例）
更に、別の第４の実施例を図１１を用いて説明する。
第４の実施例が第３の実施例と相違する点はステップＳ４において、回転角は予め定めら
れたものではなく、得られるデータから回転角を求める点であり、図１１に示すようにス
テップＳ４−１とステップＳ４−２に分けて処理する。
他のステップは第３の実施例と同様なので、説明は省略する。
【００７４】
脈拍数ＰＲを求める処理（ステップＳ４−１〜Ｓ６）では、
先ず、所定期間分の赤外光脈波データＩＲと赤色光脈波データＲとのデータを用い、図
３に示すようなグラフを描く。
そして、いかなる回転角度であれば軸方向へ射影した領域が最小になるかを求める。（
ステップＳ４−１）
次に、求められた回転角度による回転行列により、各波長の脈波データを処理し、射影
領域が最小となる軸の成分のデータからノイズが低減された波形を得る。（ステップＳ４
−２）。
このように、第４の実施例の特徴は、回転行列の回転角度が、固定化された角度ではな
く、検出される脈波データに応じて適時変更されるアダプティブ性を有する点にある。
【００７５】
以上は、動脈血酸素飽和度を測定するパルスオキシメータを例に挙げて説明したが、本
発明の技術はパルスオキシメータに限られず、特殊ヘモグロビン（一酸化炭素ヘモグロビ
ン、Ｍｅｔヘモグロビンなど）、血中に注入された色素などの血中吸光物質をパルスフォ
トメトリーの原理を用いて測定する装置(パルスフォトメトリー)にも光源の波長を選択す
ることで適用できる。
なお、本明細書中において、生体組織（媒体）を介して測定された信号を、２次元直交
座標に展開する前までは、「脈波信号」と記載し、２次元直交座標に展開中は、「脈波デ
ータ」と記載し、求められた波形は「脈波主成分波形」と区別して記載している。
【産業上の利用可能性】
【００７６】
本発明では、同一の媒体からほぼ同時に抽出される２つの同種の信号
を処理して共通の信号成分を抽出する計算処理負担を軽減した信号処理が実現できる。
また、本発明では、前記媒体の体動によるノイズが脈波信号に生じ
た場合であっても、対象物質の濃度を精度よく測定できる。
また、体動によるノイズが脈波データ信号に生じた場合であっても、脈波信号からノイ
ズを除去し、精度よく脈拍や血中の吸光物質の濃度を求めることができるので、産業上の
利用可能性は極めて大きい。
【符号の説明】
【００７７】
１ 発光素子
２ 発光素子
３ 駆動回路
４ 生体組織
５ フォトダイオード
６ 変換器
７ マルチプレクサ
８ フィルタ
９ Ａ／Ｄ変換器
１０ 処理部
１１ 表示部
１２ ＲＯＭ
１３ ＲＡＭ
【特許請求の範囲】
【請求項１】
異なる２つの波長の光を同一の生体組織に照射する発光手段と、
前記発光手段から発生し前記同一の生体組織を透過または反射した各波長の光を電気信号に変換する受光手段と、
前記受光手段で変換された前記電気信号より得られた２つの波長の離散的時系列脈波データをそれぞれ縦軸または横軸とする２次元直交座標に展開した、各波長の脈波データそれぞれについてノルム値を求め、さらにそのノルム値の比を求めるノルム比計算手段と、
前記ノルム比計算手段により計算されたノルム比に基づいて、血中の酸素の濃度を求める血中酸素濃度演算手段とを具備することを特徴とするパルスフォトメータ。
【請求項１】
異なる２つの波長の光を同一の生体組織に照射する発光手段と、
前記発光手段から発生し前記同一の生体組織を透過または反射した各波長の光を電気信号に変換する受光手段と、
前記受光手段で変換された前記電気信号より得られた２つの波長の離散的時系列脈波データをそれぞれ縦軸または横軸とする２次元直交座標に展開した、各波長の脈波データそれぞれについてノルム値を求め、さらにそのノルム値の比を求めるノルム比計算手段と、
前記ノルム比計算手段により計算されたノルム比に基づいて、血中の酸素の濃度を求める血中酸素濃度演算手段とを具備することを特徴とするパルスフォトメータ。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【公開番号】特開２００９−７２６３９（Ｐ２００９−７２６３９Ａ）
【公開日】平成２１年４月９日（２００９．４．９）
【国際特許分類】
【出願番号】特願２００９−６８７６（Ｐ２００９−６８７６）
【出願日】平成２１年１月１５日（２００９．１．１５）
【分割の表示】特願２００３−３３３６１３（Ｐ２００３−３３３６１３）の分割
【原出願日】平成１５年９月２５日（２００３．９．２５）
【出願人】（０００２３０９６２）日本光電工業株式会社 (179)
【Ｆターム（参考）】
【公開日】平成２１年４月９日（２００９．４．９）
【国際特許分類】
【出願日】平成２１年１月１５日（２００９．１．１５）
【分割の表示】特願２００３−３３３６１３（Ｐ２００３−３３３６１３）の分割
【原出願日】平成１５年９月２５日（２００３．９．２５）
【出願人】（０００２３０９６２）日本光電工業株式会社 (179)
【Ｆターム（参考）】
[ Back to top ]