JP4284674B2 - 血中吸光物質濃度測定装置 - Google Patents
血中吸光物質濃度測定装置 Download PDFInfo
- Publication number
- JP4284674B2 JP4284674B2 JP2003023990A JP2003023990A JP4284674B2 JP 4284674 B2 JP4284674 B2 JP 4284674B2 JP 2003023990 A JP2003023990 A JP 2003023990A JP 2003023990 A JP2003023990 A JP 2003023990A JP 4284674 B2 JP4284674 B2 JP 4284674B2
- Authority
- JP
- Japan
- Prior art keywords
- light intensity
- light
- blood
- transmitted light
- transmitted
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring 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/1455—Measuring 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/14551—Measuring 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/314—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry with comparison of measurements at specific and non-specific wavelengths
- G01N2021/3144—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry with comparison of measurements at specific and non-specific wavelengths for oxymetry
Description
【０００１】
【発明の属する技術分野】
本発明は、パルスオキシメータやパルス式色素希釈曲線測定装置のように、パルスフォトメトリーを原理とする血中吸光物質濃度測定装置の改良に関する。
【０００２】
【従来の技術】
パルスフォトメトリーは、パルスオキシメータだけに留まらず、現在はパルス式色素希釈法としてインドシアニングリーン(indocyanine green：ＩＣＧ）という色素を血中に投与して、その血中濃度を測定することにより、心拍出量、循環血液量、ＩＣＧの血漿消失率、クリアランスの測定を行なう装置も実用化されている（例えば、非特許文献１〜３参照）。また、一酸化炭素ヘモグロビンやメトヘモグロビンなどの異常ヘモグロビンの濃度測定や、ヘモグロビン濃度測定、血糖値測定などにも応用されている（例えば、特許文献１、特許文献２参照）。
【０００３】
例えば、２波長の光を用いて、ある物質の血液中の濃度を測定する場合、従来は血液の脈動により発生する２波長間の減光度変化分の比Φ12を求めて、Φ12と物質濃度の間に、ある一定の関係が有ることに基づいて、物質濃度を計算して求めていた（例えば、特許文献３参照）。すなわち、
C＝F(Φ12）
C:物質の血中濃度
F:一定の関係を示す関数
で表される。一般に、 n種の波長の光を用いた場合、各波長相互の減光度変化分の比Φは最大で n-1個用いる。たとえば、３波長であれば、第１波長と第２波長間の減光度変化分の比Φ12と、第１波長と第３波長間の減光度変化分の比Φ13を用いて、
C＝F(Φ12, Φ13）
となる。
物質の血中濃度C は、パルスオキシメータの場合、動脈血酸素飽和度SpO2（ヘモグロビン濃度に対する酸化ヘモグロビン濃度の比O2Hb／Hb）であり、パルス式色素希釈曲線測定の場合、色素濃度Cdとヘモグロビン濃度Hbの比Cd／Hbである。
【０００４】
しかし、このような測定方法によれば、減光度変化分比と物質濃度の間にはおおよそ一定の関係はあるものの、この関係には個体差があり、また同一個体であっても測定時点や測定部位によって変化し、測定誤差の原因になっている。例えば、パルスオキシメータでは、プローブを取付ける指を変えたり、手を上げ下げすることで、実際の動脈血酸素飽和度SpO2は一定であっても、計算値が１％程度は変化する。この測定誤差の主な原因は、次のようなことが考えられる。
（１）血液が光散乱性物質であるために、散乱により発生する減光度は、血液の厚みに依存して変化する。
（２）血液を通る光と血液を通らない光が存在する。
【０００５】
【特許文献１】
特公平３−７１１３５号公報
【特許文献２】
特開平２００２−２２８５７９号公報
【特許文献３】
特公昭５３−２６４３７号公報
【特許文献４】
特許２６０８８２８号公報
【非特許文献１】
Takehiko Iijima, et al. Cardiac output and circulating blood volume analysis by pulse dye-densitometry. J Clin Monit 1997;13: 81−89
【非特許文献２】
Takasuke Imai, et al. Measurement of cardic output by pulse dye-densitometry using indocyanine green.Anesthesiology 1997;87:816-822
【非特許文献３】
Takasuke Imai, et al. Measurement of blood concentration of indocyanine green by pulse dye-densitometry-Comparison with the conventional spectrophotometric method. J Clin Monit 1998;14: 477-484
【０００６】
【発明が解決しようとする課題】
以上のようにパルスフォトメトリーを用いて血中吸光物質濃度C を求める場合、従来は、減光度変化分比Φのみを変数とする関数を用いていた。このため、散乱により発生する減光度が、血液の厚み（ここでは変化分でなく全体の厚み）に依存して変化することは考慮されず、また、血液を通る光と、血液を通らない光（血液以外の生体組織のみを通る光）が存在するので、血液を通らない光は考慮されず、このため誤差が生じていた。
本発明は、このような従来の欠点に鑑みなされたものであり、その目的は、パルスフォトメトリーによる測定において、精度良く血中吸光物質の濃度を測定することである。
【０００７】
【課題を解決するための手段】
本発明では、血中吸光物質の濃度C を、関数F を用いて計算する場合、血液全体の厚みおよび血液以外の生体組織の厚みに関する、減光度の直流成分すなわち直流減光度Adc を考慮し、関数F の変数を、ΦおよびAdc とした。したがって、濃度C は次の関数F により表されるとした。
C＝F(Φ，Adc)
ここで、ｎ波長の光で測定する場合は、Φは最大n-1 個、Adc は最大n 個が用いられる。例えば、３波長で測定する場合、関数F は、各波長の直流減光度Adc1,Adc2,Adc3を用いて次のようになる。
C ＝F(Φ12，Φ13，Adc1，Adc2，Adc3）
直流減光度Adc は、入射光強度Iiと透過光強度Itとを用いて次のように表される。
Adc ＝log(Ii／It）＝logIi −logIt
ここで、Itは生体を透過した光の透過光強度であって連続的に測定できるが、入射光強度Iiは予め測定して知っておく必要がある。Iiを知る方法としては、本出願人の特許の特許文献４に記載した方法がある。これによれば、既知の光吸収特性を持つファントム（生体を模擬した試料）を挟み、これの透過光強度を測定することにより入射光強度を求めることができる。
【０００８】
６６０nm，８０５nm，９４０nmなどの波長は、ヘモグロビンに光吸収があるが、水による光吸収は殆どないので、これらの波長の光を生体組織に照射したときの直流的な減光度は、測定部位の血液量と主に関係がある。また、ヘモグロビンの光吸収が少なく水の光吸収が大きい１３００nmでは、直流減光度は主に生体組織の厚み（水分量）と主に関係がある。したがって、これらの波長における直流減光度を測定し、これらを計算式に変数として組み込むことにより、上記誤差要因（１）および（２）による誤差を補正し、血中吸光物質の濃度の測定精度を向上させることができる。
【０００９】
また、本発明では、直流減光度と直流透過光強度（透過光強度の直流分）は一定の関係があることから、血中吸光物質の濃度C を、直流透過光強度DC を考慮し、変数を、ΦおよびDC とした関数F1により次の式で表されるとした。
C＝F1( Φ，DC)
ここで、ｎ波長の光で測定する場合は、Φは最大n-1 個、DC は最大n 個が用いられる。例えば、３波長で測定する場合、関数F は、各波長の直流透過光強度DC1,DC2,DC3 を用いて次のようになる。
C ＝F1( Φ12，Φ13，DC1 ，DC2 ，DC3 ）
このように、減光度変化分比Φと直流透過光強度DCを変数とする関数を用いても、同様に、血中吸光物質の濃度の測定精度を向上させることができる。
【００１０】
そこで、請求項１に係る発明は、複数種の波長の光を発する光源を持ち、これらの光を生体組織に照射する光照射手段と、この光照射手段から発せられた各波長の光の生体組織への入射光強度を測定する入射光強度測定手段と、前記生体組織を透過した各波長の光を受光し、それぞれの透過光強度を測定する透過光強度測定手段と、この透過光強度測定手段が測定した透過光強度に基づいて、血液の厚み変化により発生する各波長の減光度変化分の相互の比を演算する第１の演算手段と、この第１の演算手段が求めた比と、前記入射光強度測定手段が測定した入射光強度と、前記透過光強度測定手段が測定した透過光強度とに基づいて血液中の吸光物質濃度を演算して求める第２の演算手段とを具備する。
【００１１】
また、請求項２に係る発明は、請求項１に記載の血中吸光物質濃度測定装置において、前記第２の演算手段は、各波長について、前記入射光強度測定手段が測定した入射光強度と、前記透過光強度測定手段が測定した透過光強度とに基づいて前記生体組織の直流減光度を演算して求め、この直流減光度と前記第１の演算手段が求めた比とに基づいて血液中の吸光物質濃度を演算して求めることを特徴とする。
【００１２】
また、請求項３に係る発明は、請求項２に記載の血中吸光物質濃度測定装置において、前記第２の演算手段は、各波長の直流減光度の相互の比を求め、この比と前記第１の演算手段が求めた比とに基づいて血液中の吸光物質濃度を演算して求めることを特徴とする。
【００１３】
また、請求項４に係る発明は、請求項１に記載の血中吸光物質濃度測定装置において、前記第２の演算手段は、各波長について、前記入射光強度測定手段が測定した入射光強度と、前記透過光強度測定手段が測定した透過光強度とに基づいて前記生体組織の直流透過光強度を演算して求め、この直流透過光強度と前記第１の演算手段が求めた比とに基づいて血液中の吸光物質濃度を演算して求めることを特徴ととする。
【００１４】
また、請求項５に係る発明は、請求項４に記載の血中吸光物質濃度測定装置において、前記第２の演算手段は、各波長の直流透過光強度の相互の比を求め、この比と前記第１の演算手段が求めた比とに基づいて血液中の吸光物質濃度を演算して求めることを特徴とする。
【００１５】
また、請求項６に係る発明は、請求項１乃至５のいずれか１に記載の血中吸光物質濃度測定装置において、前記光照射手段は、前記光源が発光素子から成り、この発光素子に供給する電流または印加する電圧を制御する制御部を備え、前記第２の演算手段は前記制御部により制御された電流または電圧に応じて、測定した入射光強度を補正し、この補正した入射光強度を用いて演算することを特徴とする。
【００１６】
また、請求項７に係る発明は、請求項１乃至６のいずれか１に記載の血中吸光物質濃度測定装置において、前記入射光強度測定手段は、前記光照射手段と前記透過光強度測定手段との間の所定位置に、減光度が既知のファントムが配置された場合に、前記透過光強度測定手段により測定される各波長の透過光強度に基づいて入射光強度を演算して求めることを特徴とする。
【００１７】
また、請求項８に係る発明は、請求項７に記載の血中吸光物質濃度測定装置において、前記入射光強度測定手段は、前記ファントムが前記所定位置に配置されたか否かを検出するセンサを有し、前記ファントムが前記所定位置に配置されたことを前記センサにより検出されると、入射光強度を演算して求めることを特徴とする。
【００１８】
また、請求項９に係る発明は、ｎ種の波長の光を発する光源を持ち、これらの光を生体組織に照射する光照射手段と、この光照射手段から発せられたｎ種の波長の光の生体組織への入射光強度を測定する入射光強度測定手段と、前記生体組織を透過したｎ種の波長の光を受光し、それぞれの透過光強度を測定する透過光強度測定手段と、この透過光強度測定手段が測定した透過光強度に基づいて、血液の厚み変化により発生する各波長の減光度変化分の相互の比のうち最大でｎ−１個の比を演算する第１の演算手段と、この第１の演算手段が求めた最大でｎ−１個の比と、前記入射光強度測定手段が求めた入射光強度と、前記透過光強度測定手段が求めた透過光強度とに基づいて、各波長に対応する最大でｎ個の直流減光度を演算して求め、これらの演算結果に基づいて血液中の吸光物質濃度を演算して求める第２の演算手段とを具備する。
【００１９】
また、請求項１０に係る発明は、ｎ種の波長の光を発する光源を持ち、これらの光を生体組織に照射する光照射手段と、この光照射手段から発せられたｎ種の波長の光の生体組織への入射光強度を測定する入射光強度測定手段と、前記生体組織を透過したｎ種の波長の光を受光し、それぞれの透過光強度を測定する透過光強度測定手段と、この透過光強度測定手段が測定した透過光強度に基づいて、血液の厚み変化により発生する各波長の減光度変化分の相互の比のうち最大でｎ−１個の比を演算する第１の演算手段と、この第１の演算手段が求めた最大でｎ−１個の比と、前記入射光強度測定手段が求めた入射光強度と、前記透過光強度測定手段が求めた透過光強度とに基づいて、各波長に対応する最大でｎ個の直流透過光強度を演算して求め、これらの演算結果に基づいて血液中の吸光物質濃度を演算して求める第２の演算手段とを具備する。
【００２０】
また、請求項１１に係る発明は、請求項１乃至１０のいずれか１に記載の血中吸光物質濃度測定装置において、前記光源は、血液に光吸収がある波長の光および／または血液以外の生体組織に光吸収がある波長の光を発することを特徴とする。
【００２１】
【発明の実施の形態】
［第１の実施の形態：ヘモグロビン濃度測定装置］
次に、本発明の第１の実施の形態を説明する。この第１の実施の形態は、ヘモグロビン濃度測定装置である。
【００２２】
図２は、本実施の形態の装置の全体構成を示すブロック図である。発光部１は相異なる２つの波長の光を発生するＬＥＤ２，３とこれらを駆動する駆動回路４とから成る。ＬＥＤ２から発生する光の波長を第１の波長とし、ＬＥＤ３から発生する光の波長を第２の波長とする。本装置において、第１の波長は１３００nm、第２の波長は８０５nmである。
【００２３】
受光部５は、ＬＥＤ２およびＬＥＤ３に対向して配置されたフォトダイオード６とこのフォトダイオード６の出力電流を電圧信号に変換する電流−電圧変換器７とアンプ８から成る。
【００２４】
マルチプレクサ９は、アンプ８から与えられる信号をフィルタ１０とフィルタ１１に振り分ける回路である。フィルタ１０とフィルタ１１は、各波長の透過光強度に対応する信号のノイズを除去する回路であり、これらの出力信号は、マルチプレクサ１２でそれぞれの出力のタイミングを制御されＡ／Ｄコンバータ１３に至るようにされている。Ａ／Ｄコンバータ１３はマルチプレクサ１２から与えられる信号をデジタル信号に変換する回路である。
【００２５】
ＣＰＵ１４は駆動回路４、マルチプレクサ９およびマルチプレクサ１２を制御すると共に、Ａ／Ｄコンバータ１３から与えられる信号に基づいて演算を行なう回路である。
【００２６】
メモリ１５は、ＣＰＵ１４が行なう処理のためのプログラムを格納されており、またＣＰＵ１４から与えられるデータを格納するものである。
【００２７】
表示部１６はＣＰＵ１４から与えられるデータを表示するものであり、入力部１７は複数のスイッチ（後述する校正スイッチおよび測定スイッチを含む）、複数のキーを備え、操作者の操作に応じた信号をＣＰＵ１４に出力するものである。
【００２８】
生体に装着される本装置のプローブは、ＬＥＤ２，３とフォトダイオード６を備え、これらのＬＥＤ２，３とフォトダイオード６の間で生体（例えば指先、耳朶など）３０を挟持するようになっている。
次に、本装置の動作を図１を参照して説明する。
【００２９】
＜ステップ１A ：プローブの入射光強度の測定（校正）＞
このステップ１A では、生体への入射光強度を測定する。すなわち、プローブのＬＥＤ２，３から生体に照射される光の強度を求める。本例では、既知の光吸収特性を持つファントム３０Ａを用いてこの入射光強度を求める。ファントム３０Ａとしては、例えば、乳白色のアクリル板が好適である。
【００３０】
操作者は、まず、ファントム３０ＡをプローブのＬＥＤ２，３とフォトダイオード６の間の所定位置に配置し、本装置の入力部１７の校正スイッチを操作して入射光強度の測定開始をＣＰＵ１４に指示する。これによりＬＥＤ２，３はそれぞれの波長の光を発生し、それらの光は上記のファントム３０Ａを透過してフォトダイオード６に至り、電気信号にされる。この信号は、以後、電流−電圧変換器７、アンプ８、マルチプレクサ９、フィルタ１０、１１、マルチプレクサ１２およびＡ／Ｄコンバータ１３により処理されて、ＣＰＵ１４に至り、ここで各波長の透過光強度Itcal1、Itcal2としてメモリ１５に格納される。そしてＣＰＵ１４は、測定されたItcal1、Itcal2を次式に代入して計算し、ファントム３０Ａに対する入射光強度Iical1、Iical2を求める。
Iical1＝Itcal1×exp(Af1) （１）
Iical2＝Itcal2×exp(Af2) （２）
【００３１】
式中、Af1 、Af2 はファントム３０Ａの各波長における既知の減光度であり、予めメモリ１５に格納されている。また、上記の計算結果もメモリ１５に格納される。計算結果のIical1、Iical2は、メモリ１５内の所定のエリアに書き込まれ、前回測定された値Iical1、Iical2があればこれらを書き替えるようにされている。そして、後述のステップ３A における計算では、最新のIical1、Iical2が用いられるようになっている。したがってこのステップ１は、入射光強度Iical を校正するステップである。
【００３２】
なお、入射光強度Iical の計算は、操作者がファントム３０Ａにプローブを取り付け、校正スイッチを押すことによって行なうようにしたが、プローブに光または機械的センサを設け、プローブをファントム３０Ａに取り付けると、そのセンサが感知し、これによってＣＰＵ１４は上記のような処理を開始して入射光強度Iical の計算を行なうようにしても良い。
【００３３】
また、ファントム３０Ａ自体が、プローブを収納しておくホルダーにされていて、このホルダーに上記のようなセンサを設け、プローブがそのホルダーに収納されると（もちろんファントム３０Ａの一部はＬＥＤ２，３とフォトダイオード６の間に挟持されるようになっている）、センサがその収納を感知し、これによってＣＰＵ１４は入射光強度Iical の計算を行なうようにしても良い。
【００３４】
ＬＥＤなどの発光素子は長期間使用すると発光強度が低下するし、プローブ表面の汚れなどによっても発光強度は変化するので、出荷時に決めたプローブの発光強度を継続して使うのは困難であり、測定の直前に入射光強度Iical を校正するのが良い。
【００３５】
＜ステップ２A ：脈波の測定と減光度変化分比Φの計算＞
次にステップ２A では、プローブを生体３０に取り付け、血液の脈動による各波長の減光度変化分を測定し、それら相互の比Φを求める。
【００３６】
このステップは、操作者が入力部１７の測定スイッチを操作することにより開始する。このときＣＰＵ１４は、Ａ／Ｄコンバータ１３から与えられる透過光強度It1,It2 を示す信号に基づいて、各波長の減光度A1,A2 の変化分ΔA1，ΔA2の比Φ12を求める。すなわち、次式を計算し、その結果をメモリ１５に格納する。
Φ12＝ΔA1／ΔA2＝ (AC1 ／DC1)／ (AC2 ／DC2) （３）
【００３７】
ここで、DC1,DC2 は、脈動する透過光強度It1,It2 のうちの直流分であって、直流透過光強度と称する。また、AC1,AC2 は脈動する透過光強度It1,It2 のうちの交流分であって、交流透過光強度と称する。
【００３８】
この式（３）は次のようにして得られるものである。図３（ａ）に示すように生体は、動脈血、静脈血、血液以外の組織、の各層からなると考えることができる。生体に入射光強度Iiの光を照射すると、透過光強度Itの光が測定される。このとき、生体による減光度A は、ランバート・ベールの法則に基づき、次式で表される。
A ＝log(Ii／It) ＝logIi −logIt （４）
【００３９】
次に、図３（ｂ）に示すように、動脈血の層の厚さが変化したとき、透過光強度は、変化分をΔItとすると、It−ΔItとなる。このとき、生体による減光度は、その変化分ΔA が加えられ、次式で表される。
A ＋ΔA ＝log{Ii／(It −ΔIt)}＝logIi −log(It−ΔIt) （５）
【００４０】
式（４）と式（５）からA を消去して、ΔA を求めると、ΔA は次のように入射光強度Iiを用いない式で表される。
ΔA ＝logIt −log(It−ΔIt) ＝log{It／(It −ΔIt)} （６）
この式（６）は、変形すれば、
ΔA ＝log{It／(It −ΔIt)}
＝log[1 ／{1− (ΔIt／It)}] （７）
となる。ここで (ΔIt／It) は１よりかなり小さい値をとるので（生体の動脈血層の脈動による透過光強度Itの変化分ΔItは、透過光強度Itに対し極めて小さいので）、式（７）は次の式で近似することができる。
ΔA ＝ΔIt／It （８）
したがって、Φ12の定義式であるΦ12＝ΔA1／ΔA2と、式（８）より次式が得られる。
Φ12＝ΔA1／ΔA2＝ (ΔIt1 ／It1)／ (ΔIt2 ／It2) （９）
これによれば対数計算が不要となる。この式（９）は、透過光強度がIt1,It2 を基準にしてΔIt1,ΔIt2 だけ変化した場合の式と考えられる。
【００４１】
一方、図４（ａ）に示すように、脈動する透過光強度Itのうちピーク値を直流透過光強度DCとし、ピーク値とボトム値の差（最大の変化分）の値を交流透過光強度ACとすると、この場合、透過光強度がDCを基準にしてACだけ変化したと考えられるから、ΔIt／It＝AC／DCとなり、式（９）は次のようになる。
Φ12＝ΔA1／ΔA2＝ (AC1 ／DC1)／ (AC2 ／DC2) （３）
すなわち、上記の式（３）が得られる。
【００４２】
なお、ここでは、透過光強度のうちピーク値を直流透過光強度DCとしたが、ACはDCに対し極めて小さいので、図４（ｂ）に示すように、直流透過光強度DCを、脈動する透過光強度Itのボトム値としても式（３）は成り立ち、また、図４（ｃ）に示すように、ピーク値とボトム値の平均値（ピーク値とボトム値の中間点の値）としても式（３）は成り立つ。
【００４３】
このように、直流透過光強度DCは、ピーク値からボトム値の間のどこを採っても良いので、交流透過光強度の測定時点またはその周辺時点の透過光強度Itをそのまま使っても良い。
【００４４】
このステップにおいて、Φ12の計算は、脈動する透過光強度It1,It2 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。
【００４５】
＜ステップ３A ：直流減光度の計算＞
次にステップ３A では、各波長の直流減光度Adc1, Adc2を求める。ここで、ＣＰＵ１４は、ステップ１A で求めた入射光強度Iical1、Iical2と、ステップ２A で求めた直流透過光強度DC1,DC2 を次式に代入して、各波長の直流減光度Adc1, Adc2を計算する。
Adc1＝ log(Iical1 ／DC1)＝logIical1 −logDC1 （１０）
Adc2＝ log(Iical2 ／DC2)＝logIical2 −logDC2 （１１）
【００４６】
本装置が、発光素子（ＬＥＤ２，３）に流す電流を変化させない回路構成を持つ場合は、ステップ１A で求めた入射光強度Iical1、Iical2をそのまま用いる。
【００４７】
しかし、本装置が、発光素子の電流を変化させて最適の透過光強度が得られるように調整する回路構成を持つ場合、ファントム３０Ａを用いて入射光強度測定を行なった時（校正時）の発光素子への電流値 Ccal1, Ccal2 と、生体の減光度測定を行なった時（測定時）の発光素子への電流値 Cmeas1, Cmeas2 が違うことがある。その場合ＣＰＵ１４は、入射光強度として、
IiA1 ＝Ical1 × Cmeas1 ／Ccal1 （１２）
IiA2 ＝Ical2 × Cmeas2 ／Ccal2 （１３）
のように、発光素子への通電電流に応じて補正した値IiA1,IiA2 を用いる。
【００４８】
このステップにおいても、直流減光度Adc1, Adc2の計算は、脈動する透過光強度It1,It2 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。
【００４９】
＜ステップ４A ：ヘモグロビン濃度の計算＞
次にステップ４A では、ＣＰＵ１４は、ステップ２A で求めたΦ12とステップ３A で求めた直流減光度Adc1、Adc2を変数とする次式を用いてヘモグロビン濃度Hbdcを計算する。
Hbdc＝a1×Φ12＋b1＋c1×Adc2／Adc1 （１４）
【００５０】
この式（１４）の係数a1，b1，c1は、ある母集団（例えばある選択された１０人のデータ）において式（１４）を計算した結果のHbdcと、採血してシアンメトヘモグロビン法などによって測定した正確なヘモグロビン濃度Hbs の差が最小になるように最小二乗法で予め決めた値である。
【００５１】
ここで、本装置の効果を示すため、ヘモグロビン濃度を、脈波の減光度変化分比Φ12だけを使って計算した場合と、本装置によってΦ12とAdc1,Adc2 を使って計算した場合とを、それぞれ採血法によって求めた場合とを比較する。
Φ12だけを使う場合、
Hbp ＝a2×Φ12＋b2 （１５）
により計算する。この場合の係数a2，b2も、上記母集団のデータにより上記と同様にして最小二乗法で予め決定されたものである。式（１５）を用いて計算したヘモグロビン濃度Hbp と、採血法により測定したヘモグロビン濃度Hbs との相関を図５（ａ）に示す。
【００５２】
これに対し、直流減光度を使い、式（１４）により計算した本装置によるヘモグロビン濃度Hbdcと、採血法により測定したヘモグロビン濃度Hbs との相関を図５（ｂ）に示す。これらの図を比較すると分かるように、計算式に直流減光度を入れることによって、採血法との一致が改善される。
【００５３】
また、上記の式（１４）の代わりに、次式を用いても良い。
Hbdc＝a3×Φ12＋b3＋c3×Adc2＋d3×Adc1 （１６）
この式の係数a3〜d3も上記と同様にして予め決定しておく。
【００５４】
ところで、直流減光度Adc1,Adc2 と直流透過光強度 DC1,DC2との間には逆相関の一定の関係があるので、直流減光度を用いて計算しなくとも、直流透過光強度をそのまま使って計算すれば同様の効果が得られる。この場合、直流透過光強度はプローブの入射光強度に比例するので、入射光強度 Iical1,Iical2を、発光素子（ＬＥＤ２，３）の特定の通電電流値で規格化した状態で補正する必要が有る。補正した各波長の直流透過光強度Itcomp1,Itcomp2 は次式により計算する。
Itcomp1 ＝DC1 ×(Iical1 ／Istd1)×(Cmeas1 ／Ccal1) （１７）
Itcomp2 ＝DC2 ×(Iical2 ／Istd2)×(Cmeas2 ／Ccal2) （１８）
DC1,DC2:測定した直流透過光強度
Iical1,Iical2:校正時に計算した入射光強度
Istd1,Istd2:標準とする入射光強度
Cmeas1,Cmeas2:生体を測定したときの発光素子の通電電流
Ccal1,Ccal2:入射光強度校正時の発光素子の通電電流
【００５５】
したがって上記の式（１４）、（１６）の代わりに、次式を用いても良い。
Hbdc＝a4×Φ12＋b4＋c4× Itcomp1／Itcomp2 （１９）
Hbdc＝a5×Φ12＋b5＋c5× Itcomp1＋d5×Itcomp2 （２０）
ここで、式（１９）における係数a4〜c4、式（２０）における係数a5〜d5は、それぞれ上記式（１４）の係数の決定方法と同様の方法で、予め決定されるものである。
【００５６】
このステップにおいても、Hbdcの計算は、脈動する透過光強度It1,It2 の波の１つづつについて、すなわち心拍の一拍づつについて行なう。そして、ＣＰＵ１４は、求めたHbdcをメモリ１５に格納すると共に、表示部１６に表示する。
【００５７】
本実施の形態の装置によれば、血液以外の生体組織に光吸収がある第１波長１３００ｎｍの直流減光度または直流透過光強度と、血液に光吸収がある第２波長８０５ｎｍ直流減光度または直流透過光強度を用いて計算したので、測定部位における血液と血液以外の組織からなる生体組織全体の厚さが考慮され、精度良くヘモグロビン濃度を測定することができる。
なお、本実施の形態において、第１の波長を赤色光（例えば６６０nm）、第２の波長を赤外光（例えば９４０nm）とし、酸化ヘモグロビン濃度を求めることで動脈血酸素飽和度測定に適用することができる。
【００５８】
［第２の実施の形態：一酸化炭素ヘモグロビン濃度測定装置］
次に、本発明の第２の実施の形態を説明する。この第２の実施の形態は、一酸化炭素ヘモグロビン濃度測定装置である。
【００５９】
図６は、本実施の形態の装置の全体構成を示すブロック図である。発光部１００は相異なる３つの波長の光を発生するＬＥＤ２０ａ，２０ｂ，２０ｃとこれらを駆動する駆動回路４０とから成る。
【００６０】
ＬＥＤ２０ａから発生する光の波長を第１の波長とし、ＬＥＤ２０ｂから発生する光の波長を第２の波長とし、ＬＥＤ２０ｃから発生する光の波長を第３の波長とする。本装置において、第１の波長は９４０nm、第２の波長は６６０nm、第３の波長は６２０nmである。
【００６１】
受光部５０は、ＬＥＤ２０ａ，２０ｂ，２０ｃに対向して配置されたフォトダイオード６０とこのフォトダイオード６０の出力電流を電圧信号に変換する電流−電圧変換器７０とアンプ８０から成る。
【００６２】
マルチプレクサ９０は、アンプ８０から与えられる信号を３個のフィルタ１０ａ，１０ｂ，１０ｃに振り分ける回路である。フィルタ１０ａ，１０ｂ，１０ｃは、各波長の透過光強度に対応する信号のノイズを除去する回路であり、これらの出力信号は、マルチプレクサ１２０でそれぞれの出力のタイミングを制御されＡ／Ｄコンバータ１３０に至るようにされている。Ａ／Ｄコンバータ１３０はマルチプレクサ１２０から与えられる信号をデジタル信号に変換する回路である。
【００６３】
ＣＰＵ１４０は駆動回路４０、マルチプレクサ９０およびマルチプレクサ１２０に制御信号を出力してそれらを制御すると共に、Ａ／Ｄコンバータ１３０から与えられる信号に基づいて演算を行なう回路である。
【００６４】
メモリ１５０は、ＣＰＵ１４０が行なう処理のためのプログラムを格納されており、またＣＰＵ１４０から与えられるデータを格納するものである。
【００６５】
表示部１６０はＣＰＵ１４０から与えられるデータを表示するものであり、入力部１７０は複数のスイッチ（後述する校正スイッチおよび測定スイッチを含む）、複数のキーを備え、操作者の操作に応じた信号をＣＰＵ１４０に出力するものである。本装置のプローブは、ＬＥＤ２０ａ，２０ｂ，２０ｃとフォトダイオード６０を備え、これらのＬＥＤ２０ａ，２０ｂ，２０ｃとフォトダイオード６０の間で生体３０を挟持するようになっている。
次に、本装置の動作を図７を参照して説明する。
【００６６】
＜ステップ１B ：プローブの入射光強度の測定＞
このステップ１B では、生体３０への入射光強度を測定する。ここではプローブの３個のＬＥＤ２０ａ，２０ｂ，２０ｃからそれぞれ生体３０に照射される光の強度、すなわち生体３０への入射光強度を求める。上記の第１の実施の形態の装置におけるステップ１A と同様にして、操作者は、減光度が既知のファントム３０ＡをプローブのＬＥＤ２０ａ，２０ｂ，２０ｃとフォトダイオード６０との間の所定位置に配置し、次に校正スイッチを操作して入射光強度測定開始をＣＰＵ１４０に指示する。これによりＣＰＵ１４０は、ファントム３０Ａを透過した各波長の光強度Itcal1,Itcal2,Itcal3から次式を計算して、３波長の入射光強度Iical1、Iical2、Iical3を求め、これらをメモリ１５０に格納する。
Iical1＝Itcal1×exp(Af1)
Iical2＝Itcal2×exp(Af2)
Iical3＝Itcal2×exp(Af3)
ここで、Af1,Af2,Af3 はファントム３０Ａの各波長の既知の減光度であり、予めメモリ１５０に格納されている。
【００６７】
＜ステップ２B ：脈波の測定と減光度比Φの計算＞
次にステップ２B では、操作者は、プローブを生体３０に取り付け、測定スイッチを操作する。これによりＣＰＵ１４０は、生体３０の血液の脈動による各波長の減光度変化分を測定し、それらの相互の比Φを求める。
【００６８】
すなわちこのステップでは、ＣＰＵ１４０は上記第１の実施の形態におけるステップ２A と同様にして各波長の直流透過光強度と交流透過光強度を求め、これらを用いて各波長の減光度の変化分の相互の比を求めるのであるが、ここでは、３波長の光を用いるので、次の式を計算してΦ12、Φ13を求め、その結果をメモリ１５０に格納する。
Φ12＝ΔA1／ΔA2＝ (AC1 ／DC1)／ (AC2 ／DC2) （２１）
Φ13＝ΔA1／ΔA3＝ (AC1 ／DC1)／ (AC3 ／DC3) （２２）
Φ12：第１波長と第２波長の減光度A1,A2 の変化分ΔA1, ΔA2の相互比
Φ13：第１波長と第３波長の減光度A1,A3 の変化分ΔA1, ΔA3の相互比
DC1,DC2,DC3:第１〜３波長の透過光強度It1,It2,It3 における直流透過光強度
AC1,AC2,AC3:第１〜３波長の透過光強度It1,It2,It3 における交流透過光強度
【００６９】
このステップにおいて、Φ12、Φ13の計算は、脈動する透過光強度It1,It2,Ii3 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。
【００７０】
＜ステップ３B ：直流減光度の計算＞
次にステップ３B では、各波長の直流減光度Adc1,Adc2,Adc3を求める。ここで、ＣＰＵ１４０は、ステップ１B で求めた入射光強度Iical1,Iical2,Iical3と、ステップ２B で求めた直流透過光強度DC1,DC2,DC3 を次式に代入して、各波長の直流減光度Adc1, Adc2を求め、その結果をメモリ１５０に格納する。
Adc1＝ log(Iical1 ／DC1)＝logIical1 −logDC1 （２３）
Adc2＝ log(Iical2 ／DC2)＝logIical2 −logDC2 （２４）
Adc3＝ log(Iical3 ／DC3)＝logIical3 −logDC3 （２５）
【００７１】
本装置が、発光素子（ＬＥＤ２０ａ，２０ｂ，２０ｃ）に流す電流を変化させない回路構成を持つ場合は、上記のようにステップ１B で求めた入射光強度Iical1,Iical2,Iical3をそのまま用いる。
【００７２】
しかし、本装置が、発光素子の電流を変化させて最適の透過光強度が得られるように調整する回路構成を持つ場合、ファントム３０Ａを用いて入射光強度測定を行なった時（校正時）の発光素子への電流値 Ccal1, Ccal2 ,Ccal3と生体の減光度測定を行なった時（測定時）の発光素子への電流値 Cmeas1, Cmeas2, Cmeas3 が違うことがある。その場合入射光強度IiA1,IiA2 ,IiA3 は、
IiA1 ＝Ical1 × Cmeas1 ／Ccal1 （２６）
IiA2 ＝Ical2 × Cmeas2 ／Ccal2 （２７）
IiA3 ＝Ical3 × Cmeas3 ／Ccal3 （２８）
のように、発光素子への通電電流に応じて補正した値を用いる。
【００７３】
このステップにおいても、直流減光度Adc1,Adc2,Adc3の計算は、脈動する透過光強度It1,It2,It3 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。
【００７４】
＜ステップ４B ：一酸化炭素ヘモグロビン濃度の計算＞
次にステップ４B では、ＣＰＵ１４０は、ステップ２B で求めたΦ12、Φ13とステップ３B で求めた直流減光度Adc1、Adc2、Adc3を変数とする次式を用いて一酸化炭素ヘモグロビン濃度COHbdcを計算する。
C0Hbdc＝a6×Φ12＋b6×Φ13＋c6＋d6×Adc2／Adc1＋e6×Adc3／Adc1
（２９）
【００７５】
この式（２９）の係数a6〜e6はある母集団（例えばある選択された１０人のデータ）において、式（２９）を計算した結果のCOHbdcと、採血して測定した正確な一酸化炭素ヘモグロビン値COHbs の差が最小になるように最小二乗法で予め決めた値である。
【００７６】
ここで、本装置の効果を示すため、一酸化炭素ヘモグロビン濃度を、脈波の減光度変化分比Φ12、Φ13だけを使って計算した場合と、本装置により、Φ12、Φ13と直流減光度Adc1,Adc2,Adc3を使って計算した場合とを、それぞれ採血法によって求めた場合とを比較する。
【００７７】
Φ12、Φ13だけを使う場合、
COHbp ＝a7×Φ12＋b7×Φ13＋c7 （３０）
により計算する。この場合の係数a7,b7,c7も、上記母集団のデータを用いて上記と同様にして予め決定されたものである。式（３０）を用いて計算した一酸化炭素ヘモグロビン濃度COHbp と、採血法により測定した一酸化炭素ヘモグロビン濃度COHbs との相関を図８（ａ）に示す。
【００７８】
これに対し、直流減光度を使い、式（２９）により計算した本装置による一酸化炭素ヘモグロビン濃度C0Hbdcと、採血法により測定したCOHbs との相関を図８（ｂ）に示す。
【００７９】
これらの図を比較すると分かるように、計算式に直流減光度を入れることによって、採血法との一致が改善される。
【００８０】
また、上記の式（２９）の代わりに、次式を用いても良い。
COHbdc＝a8×Φ12＋b8×Φ13＋c8＋d8×Adc1＋e8×Adc2＋f8×Adc3
（３１）
この式の係数a8〜f8も上記と同様にして予め決定しておく。
【００８１】
ところで、直流減光度Adc1,Adc2,Adc3と直流透過光強度 DC1,DC2,DC3との間には逆相関の一定の関係があるので、減光度で計算しなくても、直流透過光強度をそのまま使って計算すれば同様の効果が得られる。この場合、直流透過光強度はプローブの入射光強度に比例するので、入射光強度を、発光素子（ＬＥＤ２０ａ，２０ｂ，２０ｃ）の特定の通電電流値で規格化した状態で補正する必要が有る。補正した各波長の直流透過光強度Itcomp1,Itcomp2,Itcomp3 は次式により計算する。
Itcomp1 ＝DC1 ×(Iical1 ／Istd1)×(Cmeas1 ／Ccal1) （３２）
Itcomp2 ＝DC2 ×(Iical2 ／Istd2)×(Cmeas2 ／Ccal2) （３３）
Itcomp3 ＝DC3 ×(Iical3 ／Istd3)×(Cmeas3 ／Ccal3) （３４）
DC1,DC2,DC3:測定した透過光強度
Iical1,Iical2,Iical3: 校正時に計算した入射光強度
Istd1,Istd2,Istd3:標準とする入射光強度
Cmeas1,Cmeas2,Cmeas3: 生体を測定したときの発光素子の通電電流
Ccal1,Ccal2,Ccal3:入射光強度校正時の発光素子の通電電流
【００８２】
したがって、上記の式（２９）、（３１）の代わりに、次式を用いても良い。
COHbdc＝a9×Φ12＋b9×Φ13＋c9＋d9×Itcomp1 ／Itcomp2 ＋e9×Itcomp1 ／Itcomp3 （３５）
COHbdc＝a10 ×Φ12＋b10 ×Φ13＋ c10＋ d10×Itcomp1 ＋e10 ×Itcomp2 ＋f10 ×Itcomp3 （３６）
式（３５）における係数a9〜e9、式（３６）における係数a10 〜f10 は、それぞれ上記と同様の母集団のデータを用いて、上記と同様にして予め決定されるものである。
【００８３】
このステップにおいても、COHbdcの計算は、脈動する透過光強度It1,It2,It3 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。そして、ＣＰＵ１４０は、求めたCOHbdcをメモリ１５０に格納すると共に、表示部１６０に表示する。
【００８４】
本実施の形態の装置によれば、血液に光吸収がある第１波長９４０ｎｍ、第２波長６６０ｎｍおよび第３波長６２０ｎｍの直流減光度または直流透過光強度を用いて計算したので、測定部位の血液層全体の厚さが考慮され、精度良くヘモグロビン濃度を測定することができる。
【００８５】
［第３の実施の形態：色素濃度測定装置］
次に、本発明の第３の実施の形態を説明する。この第３の実施の形態は、色素濃度測定装置である。この装置は、第２の実施の形態と同様に３波長の光を用いるので、装置の構成ブロック図は図６と同じである。このため、本装置の構成ブロック図は省略し、以下の説明では図６を利用する。しかし、ＣＰＵ１４０が行なう処理のプログラムは、図９に示すように第２の実施の形態と異なる。また、第１の波長は９４０nm、第２の波長は６６０nm、第３の波長は８０５nmである。以下、図９を参照して本装置の動作を説明する。
【００８６】
＜ステップ１C ：プローブの入射光強度の測定（校正）＞
このステップ１C では、生体３０への入射光強度を測定する。上記の第２の実施の形態の装置におけるステップ１B と同様にして、ファントム３０Ａを用い、３波長の入射光強度Iical1、Iical2、Iical3を求め、これらをメモリ１５０に格納する。
【００８７】
＜ステップ２C ：脈波の測定と減光度変化分比Φの計算＞
次にステップ２C では、操作者は、プローブを生体３０に装着した後、色素を注入する。色素として例えばインドシアニングリーンを用いる。そして、上記第２の実施の形態におけるステップ２B と同様にして各波長の減光度の変化分の相互の比を求める。すなわち、次の式を計算してΦ12、Φ13を求め、その結果をメモリ１５０に格納する。
Φ12＝ΔA1／ΔA2＝ (AC1 ／DC1)／ (AC2 ／DC2) （２１）
Φ13＝ΔA1／ΔA3＝ (AC1 ／DC1)／ (AC3 ／DC3) （２２）
Φ12：第１波長と第２波長の減光度A1,A2 の変化分ΔA1, ΔA2の相互比
Φ13：第１波長と第３波長の減光度A1,A3 の変化分ΔA1, ΔA3の相互比
DC1,DC2,DC3:第１〜３波長の透過光強度It1,It2,It3 における直流透過光強度
AC1,AC2,AC3:第１〜３波長の透過光強度It1,It2,It3 における交流透過光強度
【００８８】
このステップにおいて、Φ12、Φ13の計算は、脈動する透過光強度It1,It2,Ii3 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。
【００８９】
＜ステップ３C ：直流減光度の計算＞
次にステップ３C では、各波長の直流減光度Adc1,Adc2,Adc3を求める。ここで、ＣＰＵ１４０は、第２の実施の形態のステップ３B と同様にして、各波長の入射光強度Iical1,Iical2,Iical3と、各波長の直流透過光強度DC1,DC2,DC3 と次式により、各波長の直流減光度Adc1, Adc2, Adc3を求め、その結果をメモリ１５０に格納する。
Adc1＝ log(Iical1 ／DC1)＝logIical1 −logDC1 （２３）
Adc2＝ log(Iical2 ／DC2)＝logIical2 −logDC2 （２４）
Adc3＝ log(Iical3 ／DC3)＝logIical3 −logDC3 （２５）
【００９０】
本装置が、発光素子（ＬＥＤ２０ａ，２０ｂ，２０ｃ）に流す電流を変化させない回路構成を持つ場合、上記のようにステップ１C で求めた入射光強度Iical1,Iical2,Iical3をそのまま用いる。
【００９１】
しかし、本装置が、発光素子の電流を変化させて最適の透過光強度が得られるように調整する回路構成を持つ場合、ファントム３０Ａを用いて入射光強度測定を行なった時（校正時）の発光素子への電流値 Ccal1, Ccal2 ,Ccal3と生体の減光度測定を行なった時（測定時）の発光素子への電流値 Cmeas1, Cmeas2, Cmeas3 が違うことがある。その場合入射光強度IiA1,IiA2 ,IiA3 は、
IiA1 ＝Ical1 × Cmeas1 ／Ccal1 （２６）
IiA2 ＝Ical2 × Cmeas2 ／Ccal2 （２７）
IiA3 ＝Ical3 × Cmeas3 ／Ccal3 （２８）
のように、発光素子への通電電流に応じて補正した値を用いる。
【００９２】
このステップにおいても、直流減光度Adc1,Adc2,Adc3の計算は、脈動する透過光強度It1,It2,It3 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。
【００９３】
＜ステップ４C ：色素濃度の計算＞
次にステップ４C では、ＣＰＵ１４０は、上記ステップ２C で求めたΦ12、Φ13と上記ステップ３C で求めた直流減光度Adc1、Adc2、Adc3を変数とする次式を用いて色素濃度Cdを計算する。
Cddc＝a11 ×Φ12＋b11 ×Φ13＋c11 ＋d11 ×Adc2／Adc1＋e11 ×Adc3／Adc1（３７）
【００９４】
この式（３７）の係数a11 〜e11 はある母集団（例えばある選択された１０人のデータ）において、式（３７）を計算した結果であるCddcと、採血して測定した正確なCds 値の差が最小になるように最小二乗法で予め決めた値である。
ここで、本装置の効果を示すため、色素濃度を、脈波の減光度変化分比Φ12、Φ13だけを使って計算した場合と、本装置によるΦ12、Φ13とAdc1,Adc2,Adc3を使って計算した場合とを、それぞれ採血法によって求めた場合とを比較する。
【００９５】
Φ12、Φ13だけを使う場合、
Cdp ＝a12 ×Φ12＋b12 ×Φ13＋c12 （３８）
により計算する。この場合の係数a12,b12,c12 も、上記母集団のデータを用いて上記と同様にして予め決定されたものである。式（３８）を用いて計算した色素濃度Cdp と、採血法により測定した色素濃度Cds との相関を図１０（ａ）に示す。これに対し、直流減光度を使い、式（３７）により計算した本装置による色素濃度Cddcと、採血法により測定した色素濃度Cds との相関を図１０（ｂ）に示す。
【００９６】
これらの図を比較すると分かるように、計算式に直流減光度を入れることによって、採血法との一致が改善される。
また、上記の式（３７）の代わりに、次式を用いても良い。
Cddc＝a13 ×Φ12＋b13 ×Φ13＋c13 ＋d13 ×Adc1＋ e13×Adc2
＋f12 ×Adc3 （３９）
この式の係数a13 〜f13 も上記と同様にして予め決定しておく。
【００９７】
ところで、直流減光度Adc1,Adc2,Adc3と直流透過光強度 DC1,DC2,DC3との間には逆相関の一定の関係があるので、減光度で計算しなくても、直流透過光強度をそのまま使っても同様の効果が得られる。この場合、直流透過光強度はプローブの入射光強度に比例するので、入射光強度を、発光素子（ＬＥＤ２０ａ，２０ｂ，２０ｃ）の特定の通電電流値で規格化した状態で補正する必要が有る。補正した各波長の直流透過光強度Itcompは次式により計算する。
Itcomp1 ＝DC1 ×(Iical1 ／Istd1)×(Cmeas1 ／Ccal1) （４０）
Itcomp2 ＝DC2 ×(Iical2 ／Istd2)×(Cmeas2 ／Ccal2) （４１）
Itcomp3 ＝DC3 ×(Iical3 ／Istd3)×(Cmeas3 ／Ccal3) （４２）
DC1,DC2,DC3:測定した透過光強度
Iical1,Iical2,Iical3: 校正時に計算した入射光強度
Istd1,Istd2,Istd3:標準とする入射光強度
Cmeas1,Cmeas2,Cmeas3: 生体を測定したときの発光素子の通電電流
Ccal1,Ccal2,Ccal3:入射光強度校正時の発光素子の通電電流
【００９８】
したがって、上記の式（３７）、（３９）の代わりに、次式を用いても良い。
Cddc＝a14 ×Φ12＋ b14×Φ13＋ c14＋ d14×Itcomp1 ／Itcomp2 ＋ e14×Itcomp1 ／Itcomp3 （４３）
Cddc＝a15×Φ12＋ b15×Φ13＋c15＋d15×Itcomp１＋e15 ×Itcomp2 ＋f15 ×Itcomp3 （４４）
【００９９】
ここで、式（４３）における係数a14 〜e14 、式（４４）における係数a15 〜f15 は、それぞれ上記と同様の母集団のデータを用いて、上記と同様にして予め決定されるものである。
【０１００】
このステップにおいても、Cddcの計算は、脈動する透過光強度It1,It2,It3 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。そして、ＣＰＵ１４０は、求めたCddcをメモリ１５０に格納すると共に、表示部１６０に表示する。
【０１０１】
本実施の形態の装置によれば、血液に光吸収がある第１波長９４０ｎｍ、第２波長８０５ｎｍおよび第３波長６６０ｎｍの直流減光度または直流透過光強度を用いて計算したので、測定部位の血液層全体の厚さが考慮され、精度良く色素濃度を測定することができる。
【０１０２】
なお、以上の説明においては、発光素子に流れる電流を制御して生体に照射する光強度を変える場合について説明したが、これは、発光素子に印加する電圧を制御して発光素子の光強度を変え、測定される入射光強度や直流透過光強度をその電圧値に応じて補正するようにしても良い。
【０１０３】
【発明の効果】
本発明によれば、パルスフォトメトリーを原理とする血中吸光物質濃度測定において、精度良く測定することができる。
【図面の簡単な説明】
【図１】第１の実施の形態のヘモグロビン濃度測定装置の動作を説明するための図。
【図２】第１の実施の形態のヘモグロビン濃度測定装置のブロック構成図。
【図３】生体へ光を照射したときの入射光強度と透過光強度との関係を説明する図。
【図４】直流透過光強度を説明する図。
【図５】第１の実施の形態のヘモグロビン濃度測定装置の効果を説明する図。
【図６】第２の実施の形態の一酸化炭素ヘモグロビン濃度測定装置のブロック構成図。
【図７】第２の実施の形態の一酸化炭素ヘモグロビン濃度測定装置の動作を説明するための図。
【図８】第２の実施の形態の一酸化炭素ヘモグロビン濃度測定装置の効果を説明する図。
【図９】第３の実施の形態の色素濃度測定装置の動作を説明するための図。
【図１０】第３の実施の形態の色素濃度測定装置の効果を説明する図。
【符号の説明】
１、１００ 発光部
５、５０ 受光部
１４、１４０ ＣＰＵ
【発明の属する技術分野】
本発明は、パルスオキシメータやパルス式色素希釈曲線測定装置のように、パルスフォトメトリーを原理とする血中吸光物質濃度測定装置の改良に関する。
【０００２】
【従来の技術】
パルスフォトメトリーは、パルスオキシメータだけに留まらず、現在はパルス式色素希釈法としてインドシアニングリーン(indocyanine green：ＩＣＧ）という色素を血中に投与して、その血中濃度を測定することにより、心拍出量、循環血液量、ＩＣＧの血漿消失率、クリアランスの測定を行なう装置も実用化されている（例えば、非特許文献１〜３参照）。また、一酸化炭素ヘモグロビンやメトヘモグロビンなどの異常ヘモグロビンの濃度測定や、ヘモグロビン濃度測定、血糖値測定などにも応用されている（例えば、特許文献１、特許文献２参照）。
【０００３】
例えば、２波長の光を用いて、ある物質の血液中の濃度を測定する場合、従来は血液の脈動により発生する２波長間の減光度変化分の比Φ12を求めて、Φ12と物質濃度の間に、ある一定の関係が有ることに基づいて、物質濃度を計算して求めていた（例えば、特許文献３参照）。すなわち、
C＝F(Φ12）
C:物質の血中濃度
F:一定の関係を示す関数
で表される。一般に、 n種の波長の光を用いた場合、各波長相互の減光度変化分の比Φは最大で n-1個用いる。たとえば、３波長であれば、第１波長と第２波長間の減光度変化分の比Φ12と、第１波長と第３波長間の減光度変化分の比Φ13を用いて、
C＝F(Φ12, Φ13）
となる。
物質の血中濃度C は、パルスオキシメータの場合、動脈血酸素飽和度SpO2（ヘモグロビン濃度に対する酸化ヘモグロビン濃度の比O2Hb／Hb）であり、パルス式色素希釈曲線測定の場合、色素濃度Cdとヘモグロビン濃度Hbの比Cd／Hbである。
【０００４】
しかし、このような測定方法によれば、減光度変化分比と物質濃度の間にはおおよそ一定の関係はあるものの、この関係には個体差があり、また同一個体であっても測定時点や測定部位によって変化し、測定誤差の原因になっている。例えば、パルスオキシメータでは、プローブを取付ける指を変えたり、手を上げ下げすることで、実際の動脈血酸素飽和度SpO2は一定であっても、計算値が１％程度は変化する。この測定誤差の主な原因は、次のようなことが考えられる。
（１）血液が光散乱性物質であるために、散乱により発生する減光度は、血液の厚みに依存して変化する。
（２）血液を通る光と血液を通らない光が存在する。
【０００５】
【特許文献１】
特公平３−７１１３５号公報
【特許文献２】
特開平２００２−２２８５７９号公報
【特許文献３】
特公昭５３−２６４３７号公報
【特許文献４】
特許２６０８８２８号公報
【非特許文献１】
Takehiko Iijima, et al. Cardiac output and circulating blood volume analysis by pulse dye-densitometry. J Clin Monit 1997;13: 81−89
【非特許文献２】
Takasuke Imai, et al. Measurement of cardic output by pulse dye-densitometry using indocyanine green.Anesthesiology 1997;87:816-822
【非特許文献３】
Takasuke Imai, et al. Measurement of blood concentration of indocyanine green by pulse dye-densitometry-Comparison with the conventional spectrophotometric method. J Clin Monit 1998;14: 477-484
【０００６】
【発明が解決しようとする課題】
以上のようにパルスフォトメトリーを用いて血中吸光物質濃度C を求める場合、従来は、減光度変化分比Φのみを変数とする関数を用いていた。このため、散乱により発生する減光度が、血液の厚み（ここでは変化分でなく全体の厚み）に依存して変化することは考慮されず、また、血液を通る光と、血液を通らない光（血液以外の生体組織のみを通る光）が存在するので、血液を通らない光は考慮されず、このため誤差が生じていた。
本発明は、このような従来の欠点に鑑みなされたものであり、その目的は、パルスフォトメトリーによる測定において、精度良く血中吸光物質の濃度を測定することである。
【０００７】
【課題を解決するための手段】
本発明では、血中吸光物質の濃度C を、関数F を用いて計算する場合、血液全体の厚みおよび血液以外の生体組織の厚みに関する、減光度の直流成分すなわち直流減光度Adc を考慮し、関数F の変数を、ΦおよびAdc とした。したがって、濃度C は次の関数F により表されるとした。
C＝F(Φ，Adc)
ここで、ｎ波長の光で測定する場合は、Φは最大n-1 個、Adc は最大n 個が用いられる。例えば、３波長で測定する場合、関数F は、各波長の直流減光度Adc1,Adc2,Adc3を用いて次のようになる。
C ＝F(Φ12，Φ13，Adc1，Adc2，Adc3）
直流減光度Adc は、入射光強度Iiと透過光強度Itとを用いて次のように表される。
Adc ＝log(Ii／It）＝logIi −logIt
ここで、Itは生体を透過した光の透過光強度であって連続的に測定できるが、入射光強度Iiは予め測定して知っておく必要がある。Iiを知る方法としては、本出願人の特許の特許文献４に記載した方法がある。これによれば、既知の光吸収特性を持つファントム（生体を模擬した試料）を挟み、これの透過光強度を測定することにより入射光強度を求めることができる。
【０００８】
６６０nm，８０５nm，９４０nmなどの波長は、ヘモグロビンに光吸収があるが、水による光吸収は殆どないので、これらの波長の光を生体組織に照射したときの直流的な減光度は、測定部位の血液量と主に関係がある。また、ヘモグロビンの光吸収が少なく水の光吸収が大きい１３００nmでは、直流減光度は主に生体組織の厚み（水分量）と主に関係がある。したがって、これらの波長における直流減光度を測定し、これらを計算式に変数として組み込むことにより、上記誤差要因（１）および（２）による誤差を補正し、血中吸光物質の濃度の測定精度を向上させることができる。
【０００９】
また、本発明では、直流減光度と直流透過光強度（透過光強度の直流分）は一定の関係があることから、血中吸光物質の濃度C を、直流透過光強度DC を考慮し、変数を、ΦおよびDC とした関数F1により次の式で表されるとした。
C＝F1( Φ，DC)
ここで、ｎ波長の光で測定する場合は、Φは最大n-1 個、DC は最大n 個が用いられる。例えば、３波長で測定する場合、関数F は、各波長の直流透過光強度DC1,DC2,DC3 を用いて次のようになる。
C ＝F1( Φ12，Φ13，DC1 ，DC2 ，DC3 ）
このように、減光度変化分比Φと直流透過光強度DCを変数とする関数を用いても、同様に、血中吸光物質の濃度の測定精度を向上させることができる。
【００１０】
そこで、請求項１に係る発明は、複数種の波長の光を発する光源を持ち、これらの光を生体組織に照射する光照射手段と、この光照射手段から発せられた各波長の光の生体組織への入射光強度を測定する入射光強度測定手段と、前記生体組織を透過した各波長の光を受光し、それぞれの透過光強度を測定する透過光強度測定手段と、この透過光強度測定手段が測定した透過光強度に基づいて、血液の厚み変化により発生する各波長の減光度変化分の相互の比を演算する第１の演算手段と、この第１の演算手段が求めた比と、前記入射光強度測定手段が測定した入射光強度と、前記透過光強度測定手段が測定した透過光強度とに基づいて血液中の吸光物質濃度を演算して求める第２の演算手段とを具備する。
【００１１】
また、請求項２に係る発明は、請求項１に記載の血中吸光物質濃度測定装置において、前記第２の演算手段は、各波長について、前記入射光強度測定手段が測定した入射光強度と、前記透過光強度測定手段が測定した透過光強度とに基づいて前記生体組織の直流減光度を演算して求め、この直流減光度と前記第１の演算手段が求めた比とに基づいて血液中の吸光物質濃度を演算して求めることを特徴とする。
【００１２】
また、請求項３に係る発明は、請求項２に記載の血中吸光物質濃度測定装置において、前記第２の演算手段は、各波長の直流減光度の相互の比を求め、この比と前記第１の演算手段が求めた比とに基づいて血液中の吸光物質濃度を演算して求めることを特徴とする。
【００１３】
また、請求項４に係る発明は、請求項１に記載の血中吸光物質濃度測定装置において、前記第２の演算手段は、各波長について、前記入射光強度測定手段が測定した入射光強度と、前記透過光強度測定手段が測定した透過光強度とに基づいて前記生体組織の直流透過光強度を演算して求め、この直流透過光強度と前記第１の演算手段が求めた比とに基づいて血液中の吸光物質濃度を演算して求めることを特徴ととする。
【００１４】
また、請求項５に係る発明は、請求項４に記載の血中吸光物質濃度測定装置において、前記第２の演算手段は、各波長の直流透過光強度の相互の比を求め、この比と前記第１の演算手段が求めた比とに基づいて血液中の吸光物質濃度を演算して求めることを特徴とする。
【００１５】
また、請求項６に係る発明は、請求項１乃至５のいずれか１に記載の血中吸光物質濃度測定装置において、前記光照射手段は、前記光源が発光素子から成り、この発光素子に供給する電流または印加する電圧を制御する制御部を備え、前記第２の演算手段は前記制御部により制御された電流または電圧に応じて、測定した入射光強度を補正し、この補正した入射光強度を用いて演算することを特徴とする。
【００１６】
また、請求項７に係る発明は、請求項１乃至６のいずれか１に記載の血中吸光物質濃度測定装置において、前記入射光強度測定手段は、前記光照射手段と前記透過光強度測定手段との間の所定位置に、減光度が既知のファントムが配置された場合に、前記透過光強度測定手段により測定される各波長の透過光強度に基づいて入射光強度を演算して求めることを特徴とする。
【００１７】
また、請求項８に係る発明は、請求項７に記載の血中吸光物質濃度測定装置において、前記入射光強度測定手段は、前記ファントムが前記所定位置に配置されたか否かを検出するセンサを有し、前記ファントムが前記所定位置に配置されたことを前記センサにより検出されると、入射光強度を演算して求めることを特徴とする。
【００１８】
また、請求項９に係る発明は、ｎ種の波長の光を発する光源を持ち、これらの光を生体組織に照射する光照射手段と、この光照射手段から発せられたｎ種の波長の光の生体組織への入射光強度を測定する入射光強度測定手段と、前記生体組織を透過したｎ種の波長の光を受光し、それぞれの透過光強度を測定する透過光強度測定手段と、この透過光強度測定手段が測定した透過光強度に基づいて、血液の厚み変化により発生する各波長の減光度変化分の相互の比のうち最大でｎ−１個の比を演算する第１の演算手段と、この第１の演算手段が求めた最大でｎ−１個の比と、前記入射光強度測定手段が求めた入射光強度と、前記透過光強度測定手段が求めた透過光強度とに基づいて、各波長に対応する最大でｎ個の直流減光度を演算して求め、これらの演算結果に基づいて血液中の吸光物質濃度を演算して求める第２の演算手段とを具備する。
【００１９】
また、請求項１０に係る発明は、ｎ種の波長の光を発する光源を持ち、これらの光を生体組織に照射する光照射手段と、この光照射手段から発せられたｎ種の波長の光の生体組織への入射光強度を測定する入射光強度測定手段と、前記生体組織を透過したｎ種の波長の光を受光し、それぞれの透過光強度を測定する透過光強度測定手段と、この透過光強度測定手段が測定した透過光強度に基づいて、血液の厚み変化により発生する各波長の減光度変化分の相互の比のうち最大でｎ−１個の比を演算する第１の演算手段と、この第１の演算手段が求めた最大でｎ−１個の比と、前記入射光強度測定手段が求めた入射光強度と、前記透過光強度測定手段が求めた透過光強度とに基づいて、各波長に対応する最大でｎ個の直流透過光強度を演算して求め、これらの演算結果に基づいて血液中の吸光物質濃度を演算して求める第２の演算手段とを具備する。
【００２０】
また、請求項１１に係る発明は、請求項１乃至１０のいずれか１に記載の血中吸光物質濃度測定装置において、前記光源は、血液に光吸収がある波長の光および／または血液以外の生体組織に光吸収がある波長の光を発することを特徴とする。
【００２１】
【発明の実施の形態】
［第１の実施の形態：ヘモグロビン濃度測定装置］
次に、本発明の第１の実施の形態を説明する。この第１の実施の形態は、ヘモグロビン濃度測定装置である。
【００２２】
図２は、本実施の形態の装置の全体構成を示すブロック図である。発光部１は相異なる２つの波長の光を発生するＬＥＤ２，３とこれらを駆動する駆動回路４とから成る。ＬＥＤ２から発生する光の波長を第１の波長とし、ＬＥＤ３から発生する光の波長を第２の波長とする。本装置において、第１の波長は１３００nm、第２の波長は８０５nmである。
【００２３】
受光部５は、ＬＥＤ２およびＬＥＤ３に対向して配置されたフォトダイオード６とこのフォトダイオード６の出力電流を電圧信号に変換する電流−電圧変換器７とアンプ８から成る。
【００２４】
マルチプレクサ９は、アンプ８から与えられる信号をフィルタ１０とフィルタ１１に振り分ける回路である。フィルタ１０とフィルタ１１は、各波長の透過光強度に対応する信号のノイズを除去する回路であり、これらの出力信号は、マルチプレクサ１２でそれぞれの出力のタイミングを制御されＡ／Ｄコンバータ１３に至るようにされている。Ａ／Ｄコンバータ１３はマルチプレクサ１２から与えられる信号をデジタル信号に変換する回路である。
【００２５】
ＣＰＵ１４は駆動回路４、マルチプレクサ９およびマルチプレクサ１２を制御すると共に、Ａ／Ｄコンバータ１３から与えられる信号に基づいて演算を行なう回路である。
【００２６】
メモリ１５は、ＣＰＵ１４が行なう処理のためのプログラムを格納されており、またＣＰＵ１４から与えられるデータを格納するものである。
【００２７】
表示部１６はＣＰＵ１４から与えられるデータを表示するものであり、入力部１７は複数のスイッチ（後述する校正スイッチおよび測定スイッチを含む）、複数のキーを備え、操作者の操作に応じた信号をＣＰＵ１４に出力するものである。
【００２８】
生体に装着される本装置のプローブは、ＬＥＤ２，３とフォトダイオード６を備え、これらのＬＥＤ２，３とフォトダイオード６の間で生体（例えば指先、耳朶など）３０を挟持するようになっている。
次に、本装置の動作を図１を参照して説明する。
【００２９】
＜ステップ１A ：プローブの入射光強度の測定（校正）＞
このステップ１A では、生体への入射光強度を測定する。すなわち、プローブのＬＥＤ２，３から生体に照射される光の強度を求める。本例では、既知の光吸収特性を持つファントム３０Ａを用いてこの入射光強度を求める。ファントム３０Ａとしては、例えば、乳白色のアクリル板が好適である。
【００３０】
操作者は、まず、ファントム３０ＡをプローブのＬＥＤ２，３とフォトダイオード６の間の所定位置に配置し、本装置の入力部１７の校正スイッチを操作して入射光強度の測定開始をＣＰＵ１４に指示する。これによりＬＥＤ２，３はそれぞれの波長の光を発生し、それらの光は上記のファントム３０Ａを透過してフォトダイオード６に至り、電気信号にされる。この信号は、以後、電流−電圧変換器７、アンプ８、マルチプレクサ９、フィルタ１０、１１、マルチプレクサ１２およびＡ／Ｄコンバータ１３により処理されて、ＣＰＵ１４に至り、ここで各波長の透過光強度Itcal1、Itcal2としてメモリ１５に格納される。そしてＣＰＵ１４は、測定されたItcal1、Itcal2を次式に代入して計算し、ファントム３０Ａに対する入射光強度Iical1、Iical2を求める。
Iical1＝Itcal1×exp(Af1) （１）
Iical2＝Itcal2×exp(Af2) （２）
【００３１】
式中、Af1 、Af2 はファントム３０Ａの各波長における既知の減光度であり、予めメモリ１５に格納されている。また、上記の計算結果もメモリ１５に格納される。計算結果のIical1、Iical2は、メモリ１５内の所定のエリアに書き込まれ、前回測定された値Iical1、Iical2があればこれらを書き替えるようにされている。そして、後述のステップ３A における計算では、最新のIical1、Iical2が用いられるようになっている。したがってこのステップ１は、入射光強度Iical を校正するステップである。
【００３２】
なお、入射光強度Iical の計算は、操作者がファントム３０Ａにプローブを取り付け、校正スイッチを押すことによって行なうようにしたが、プローブに光または機械的センサを設け、プローブをファントム３０Ａに取り付けると、そのセンサが感知し、これによってＣＰＵ１４は上記のような処理を開始して入射光強度Iical の計算を行なうようにしても良い。
【００３３】
また、ファントム３０Ａ自体が、プローブを収納しておくホルダーにされていて、このホルダーに上記のようなセンサを設け、プローブがそのホルダーに収納されると（もちろんファントム３０Ａの一部はＬＥＤ２，３とフォトダイオード６の間に挟持されるようになっている）、センサがその収納を感知し、これによってＣＰＵ１４は入射光強度Iical の計算を行なうようにしても良い。
【００３４】
ＬＥＤなどの発光素子は長期間使用すると発光強度が低下するし、プローブ表面の汚れなどによっても発光強度は変化するので、出荷時に決めたプローブの発光強度を継続して使うのは困難であり、測定の直前に入射光強度Iical を校正するのが良い。
【００３５】
＜ステップ２A ：脈波の測定と減光度変化分比Φの計算＞
次にステップ２A では、プローブを生体３０に取り付け、血液の脈動による各波長の減光度変化分を測定し、それら相互の比Φを求める。
【００３６】
このステップは、操作者が入力部１７の測定スイッチを操作することにより開始する。このときＣＰＵ１４は、Ａ／Ｄコンバータ１３から与えられる透過光強度It1,It2 を示す信号に基づいて、各波長の減光度A1,A2 の変化分ΔA1，ΔA2の比Φ12を求める。すなわち、次式を計算し、その結果をメモリ１５に格納する。
Φ12＝ΔA1／ΔA2＝ (AC1 ／DC1)／ (AC2 ／DC2) （３）
【００３７】
ここで、DC1,DC2 は、脈動する透過光強度It1,It2 のうちの直流分であって、直流透過光強度と称する。また、AC1,AC2 は脈動する透過光強度It1,It2 のうちの交流分であって、交流透過光強度と称する。
【００３８】
この式（３）は次のようにして得られるものである。図３（ａ）に示すように生体は、動脈血、静脈血、血液以外の組織、の各層からなると考えることができる。生体に入射光強度Iiの光を照射すると、透過光強度Itの光が測定される。このとき、生体による減光度A は、ランバート・ベールの法則に基づき、次式で表される。
A ＝log(Ii／It) ＝logIi −logIt （４）
【００３９】
次に、図３（ｂ）に示すように、動脈血の層の厚さが変化したとき、透過光強度は、変化分をΔItとすると、It−ΔItとなる。このとき、生体による減光度は、その変化分ΔA が加えられ、次式で表される。
A ＋ΔA ＝log{Ii／(It −ΔIt)}＝logIi −log(It−ΔIt) （５）
【００４０】
式（４）と式（５）からA を消去して、ΔA を求めると、ΔA は次のように入射光強度Iiを用いない式で表される。
ΔA ＝logIt −log(It−ΔIt) ＝log{It／(It −ΔIt)} （６）
この式（６）は、変形すれば、
ΔA ＝log{It／(It −ΔIt)}
＝log[1 ／{1− (ΔIt／It)}] （７）
となる。ここで (ΔIt／It) は１よりかなり小さい値をとるので（生体の動脈血層の脈動による透過光強度Itの変化分ΔItは、透過光強度Itに対し極めて小さいので）、式（７）は次の式で近似することができる。
ΔA ＝ΔIt／It （８）
したがって、Φ12の定義式であるΦ12＝ΔA1／ΔA2と、式（８）より次式が得られる。
Φ12＝ΔA1／ΔA2＝ (ΔIt1 ／It1)／ (ΔIt2 ／It2) （９）
これによれば対数計算が不要となる。この式（９）は、透過光強度がIt1,It2 を基準にしてΔIt1,ΔIt2 だけ変化した場合の式と考えられる。
【００４１】
一方、図４（ａ）に示すように、脈動する透過光強度Itのうちピーク値を直流透過光強度DCとし、ピーク値とボトム値の差（最大の変化分）の値を交流透過光強度ACとすると、この場合、透過光強度がDCを基準にしてACだけ変化したと考えられるから、ΔIt／It＝AC／DCとなり、式（９）は次のようになる。
Φ12＝ΔA1／ΔA2＝ (AC1 ／DC1)／ (AC2 ／DC2) （３）
すなわち、上記の式（３）が得られる。
【００４２】
なお、ここでは、透過光強度のうちピーク値を直流透過光強度DCとしたが、ACはDCに対し極めて小さいので、図４（ｂ）に示すように、直流透過光強度DCを、脈動する透過光強度Itのボトム値としても式（３）は成り立ち、また、図４（ｃ）に示すように、ピーク値とボトム値の平均値（ピーク値とボトム値の中間点の値）としても式（３）は成り立つ。
【００４３】
このように、直流透過光強度DCは、ピーク値からボトム値の間のどこを採っても良いので、交流透過光強度の測定時点またはその周辺時点の透過光強度Itをそのまま使っても良い。
【００４４】
このステップにおいて、Φ12の計算は、脈動する透過光強度It1,It2 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。
【００４５】
＜ステップ３A ：直流減光度の計算＞
次にステップ３A では、各波長の直流減光度Adc1, Adc2を求める。ここで、ＣＰＵ１４は、ステップ１A で求めた入射光強度Iical1、Iical2と、ステップ２A で求めた直流透過光強度DC1,DC2 を次式に代入して、各波長の直流減光度Adc1, Adc2を計算する。
Adc1＝ log(Iical1 ／DC1)＝logIical1 −logDC1 （１０）
Adc2＝ log(Iical2 ／DC2)＝logIical2 −logDC2 （１１）
【００４６】
本装置が、発光素子（ＬＥＤ２，３）に流す電流を変化させない回路構成を持つ場合は、ステップ１A で求めた入射光強度Iical1、Iical2をそのまま用いる。
【００４７】
しかし、本装置が、発光素子の電流を変化させて最適の透過光強度が得られるように調整する回路構成を持つ場合、ファントム３０Ａを用いて入射光強度測定を行なった時（校正時）の発光素子への電流値 Ccal1, Ccal2 と、生体の減光度測定を行なった時（測定時）の発光素子への電流値 Cmeas1, Cmeas2 が違うことがある。その場合ＣＰＵ１４は、入射光強度として、
IiA1 ＝Ical1 × Cmeas1 ／Ccal1 （１２）
IiA2 ＝Ical2 × Cmeas2 ／Ccal2 （１３）
のように、発光素子への通電電流に応じて補正した値IiA1,IiA2 を用いる。
【００４８】
このステップにおいても、直流減光度Adc1, Adc2の計算は、脈動する透過光強度It1,It2 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。
【００４９】
＜ステップ４A ：ヘモグロビン濃度の計算＞
次にステップ４A では、ＣＰＵ１４は、ステップ２A で求めたΦ12とステップ３A で求めた直流減光度Adc1、Adc2を変数とする次式を用いてヘモグロビン濃度Hbdcを計算する。
Hbdc＝a1×Φ12＋b1＋c1×Adc2／Adc1 （１４）
【００５０】
この式（１４）の係数a1，b1，c1は、ある母集団（例えばある選択された１０人のデータ）において式（１４）を計算した結果のHbdcと、採血してシアンメトヘモグロビン法などによって測定した正確なヘモグロビン濃度Hbs の差が最小になるように最小二乗法で予め決めた値である。
【００５１】
ここで、本装置の効果を示すため、ヘモグロビン濃度を、脈波の減光度変化分比Φ12だけを使って計算した場合と、本装置によってΦ12とAdc1,Adc2 を使って計算した場合とを、それぞれ採血法によって求めた場合とを比較する。
Φ12だけを使う場合、
Hbp ＝a2×Φ12＋b2 （１５）
により計算する。この場合の係数a2，b2も、上記母集団のデータにより上記と同様にして最小二乗法で予め決定されたものである。式（１５）を用いて計算したヘモグロビン濃度Hbp と、採血法により測定したヘモグロビン濃度Hbs との相関を図５（ａ）に示す。
【００５２】
これに対し、直流減光度を使い、式（１４）により計算した本装置によるヘモグロビン濃度Hbdcと、採血法により測定したヘモグロビン濃度Hbs との相関を図５（ｂ）に示す。これらの図を比較すると分かるように、計算式に直流減光度を入れることによって、採血法との一致が改善される。
【００５３】
また、上記の式（１４）の代わりに、次式を用いても良い。
Hbdc＝a3×Φ12＋b3＋c3×Adc2＋d3×Adc1 （１６）
この式の係数a3〜d3も上記と同様にして予め決定しておく。
【００５４】
ところで、直流減光度Adc1,Adc2 と直流透過光強度 DC1,DC2との間には逆相関の一定の関係があるので、直流減光度を用いて計算しなくとも、直流透過光強度をそのまま使って計算すれば同様の効果が得られる。この場合、直流透過光強度はプローブの入射光強度に比例するので、入射光強度 Iical1,Iical2を、発光素子（ＬＥＤ２，３）の特定の通電電流値で規格化した状態で補正する必要が有る。補正した各波長の直流透過光強度Itcomp1,Itcomp2 は次式により計算する。
Itcomp1 ＝DC1 ×(Iical1 ／Istd1)×(Cmeas1 ／Ccal1) （１７）
Itcomp2 ＝DC2 ×(Iical2 ／Istd2)×(Cmeas2 ／Ccal2) （１８）
DC1,DC2:測定した直流透過光強度
Iical1,Iical2:校正時に計算した入射光強度
Istd1,Istd2:標準とする入射光強度
Cmeas1,Cmeas2:生体を測定したときの発光素子の通電電流
Ccal1,Ccal2:入射光強度校正時の発光素子の通電電流
【００５５】
したがって上記の式（１４）、（１６）の代わりに、次式を用いても良い。
Hbdc＝a4×Φ12＋b4＋c4× Itcomp1／Itcomp2 （１９）
Hbdc＝a5×Φ12＋b5＋c5× Itcomp1＋d5×Itcomp2 （２０）
ここで、式（１９）における係数a4〜c4、式（２０）における係数a5〜d5は、それぞれ上記式（１４）の係数の決定方法と同様の方法で、予め決定されるものである。
【００５６】
このステップにおいても、Hbdcの計算は、脈動する透過光強度It1,It2 の波の１つづつについて、すなわち心拍の一拍づつについて行なう。そして、ＣＰＵ１４は、求めたHbdcをメモリ１５に格納すると共に、表示部１６に表示する。
【００５７】
本実施の形態の装置によれば、血液以外の生体組織に光吸収がある第１波長１３００ｎｍの直流減光度または直流透過光強度と、血液に光吸収がある第２波長８０５ｎｍ直流減光度または直流透過光強度を用いて計算したので、測定部位における血液と血液以外の組織からなる生体組織全体の厚さが考慮され、精度良くヘモグロビン濃度を測定することができる。
なお、本実施の形態において、第１の波長を赤色光（例えば６６０nm）、第２の波長を赤外光（例えば９４０nm）とし、酸化ヘモグロビン濃度を求めることで動脈血酸素飽和度測定に適用することができる。
【００５８】
［第２の実施の形態：一酸化炭素ヘモグロビン濃度測定装置］
次に、本発明の第２の実施の形態を説明する。この第２の実施の形態は、一酸化炭素ヘモグロビン濃度測定装置である。
【００５９】
図６は、本実施の形態の装置の全体構成を示すブロック図である。発光部１００は相異なる３つの波長の光を発生するＬＥＤ２０ａ，２０ｂ，２０ｃとこれらを駆動する駆動回路４０とから成る。
【００６０】
ＬＥＤ２０ａから発生する光の波長を第１の波長とし、ＬＥＤ２０ｂから発生する光の波長を第２の波長とし、ＬＥＤ２０ｃから発生する光の波長を第３の波長とする。本装置において、第１の波長は９４０nm、第２の波長は６６０nm、第３の波長は６２０nmである。
【００６１】
受光部５０は、ＬＥＤ２０ａ，２０ｂ，２０ｃに対向して配置されたフォトダイオード６０とこのフォトダイオード６０の出力電流を電圧信号に変換する電流−電圧変換器７０とアンプ８０から成る。
【００６２】
マルチプレクサ９０は、アンプ８０から与えられる信号を３個のフィルタ１０ａ，１０ｂ，１０ｃに振り分ける回路である。フィルタ１０ａ，１０ｂ，１０ｃは、各波長の透過光強度に対応する信号のノイズを除去する回路であり、これらの出力信号は、マルチプレクサ１２０でそれぞれの出力のタイミングを制御されＡ／Ｄコンバータ１３０に至るようにされている。Ａ／Ｄコンバータ１３０はマルチプレクサ１２０から与えられる信号をデジタル信号に変換する回路である。
【００６３】
ＣＰＵ１４０は駆動回路４０、マルチプレクサ９０およびマルチプレクサ１２０に制御信号を出力してそれらを制御すると共に、Ａ／Ｄコンバータ１３０から与えられる信号に基づいて演算を行なう回路である。
【００６４】
メモリ１５０は、ＣＰＵ１４０が行なう処理のためのプログラムを格納されており、またＣＰＵ１４０から与えられるデータを格納するものである。
【００６５】
表示部１６０はＣＰＵ１４０から与えられるデータを表示するものであり、入力部１７０は複数のスイッチ（後述する校正スイッチおよび測定スイッチを含む）、複数のキーを備え、操作者の操作に応じた信号をＣＰＵ１４０に出力するものである。本装置のプローブは、ＬＥＤ２０ａ，２０ｂ，２０ｃとフォトダイオード６０を備え、これらのＬＥＤ２０ａ，２０ｂ，２０ｃとフォトダイオード６０の間で生体３０を挟持するようになっている。
次に、本装置の動作を図７を参照して説明する。
【００６６】
＜ステップ１B ：プローブの入射光強度の測定＞
このステップ１B では、生体３０への入射光強度を測定する。ここではプローブの３個のＬＥＤ２０ａ，２０ｂ，２０ｃからそれぞれ生体３０に照射される光の強度、すなわち生体３０への入射光強度を求める。上記の第１の実施の形態の装置におけるステップ１A と同様にして、操作者は、減光度が既知のファントム３０ＡをプローブのＬＥＤ２０ａ，２０ｂ，２０ｃとフォトダイオード６０との間の所定位置に配置し、次に校正スイッチを操作して入射光強度測定開始をＣＰＵ１４０に指示する。これによりＣＰＵ１４０は、ファントム３０Ａを透過した各波長の光強度Itcal1,Itcal2,Itcal3から次式を計算して、３波長の入射光強度Iical1、Iical2、Iical3を求め、これらをメモリ１５０に格納する。
Iical1＝Itcal1×exp(Af1)
Iical2＝Itcal2×exp(Af2)
Iical3＝Itcal2×exp(Af3)
ここで、Af1,Af2,Af3 はファントム３０Ａの各波長の既知の減光度であり、予めメモリ１５０に格納されている。
【００６７】
＜ステップ２B ：脈波の測定と減光度比Φの計算＞
次にステップ２B では、操作者は、プローブを生体３０に取り付け、測定スイッチを操作する。これによりＣＰＵ１４０は、生体３０の血液の脈動による各波長の減光度変化分を測定し、それらの相互の比Φを求める。
【００６８】
すなわちこのステップでは、ＣＰＵ１４０は上記第１の実施の形態におけるステップ２A と同様にして各波長の直流透過光強度と交流透過光強度を求め、これらを用いて各波長の減光度の変化分の相互の比を求めるのであるが、ここでは、３波長の光を用いるので、次の式を計算してΦ12、Φ13を求め、その結果をメモリ１５０に格納する。
Φ12＝ΔA1／ΔA2＝ (AC1 ／DC1)／ (AC2 ／DC2) （２１）
Φ13＝ΔA1／ΔA3＝ (AC1 ／DC1)／ (AC3 ／DC3) （２２）
Φ12：第１波長と第２波長の減光度A1,A2 の変化分ΔA1, ΔA2の相互比
Φ13：第１波長と第３波長の減光度A1,A3 の変化分ΔA1, ΔA3の相互比
DC1,DC2,DC3:第１〜３波長の透過光強度It1,It2,It3 における直流透過光強度
AC1,AC2,AC3:第１〜３波長の透過光強度It1,It2,It3 における交流透過光強度
【００６９】
このステップにおいて、Φ12、Φ13の計算は、脈動する透過光強度It1,It2,Ii3 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。
【００７０】
＜ステップ３B ：直流減光度の計算＞
次にステップ３B では、各波長の直流減光度Adc1,Adc2,Adc3を求める。ここで、ＣＰＵ１４０は、ステップ１B で求めた入射光強度Iical1,Iical2,Iical3と、ステップ２B で求めた直流透過光強度DC1,DC2,DC3 を次式に代入して、各波長の直流減光度Adc1, Adc2を求め、その結果をメモリ１５０に格納する。
Adc1＝ log(Iical1 ／DC1)＝logIical1 −logDC1 （２３）
Adc2＝ log(Iical2 ／DC2)＝logIical2 −logDC2 （２４）
Adc3＝ log(Iical3 ／DC3)＝logIical3 −logDC3 （２５）
【００７１】
本装置が、発光素子（ＬＥＤ２０ａ，２０ｂ，２０ｃ）に流す電流を変化させない回路構成を持つ場合は、上記のようにステップ１B で求めた入射光強度Iical1,Iical2,Iical3をそのまま用いる。
【００７２】
しかし、本装置が、発光素子の電流を変化させて最適の透過光強度が得られるように調整する回路構成を持つ場合、ファントム３０Ａを用いて入射光強度測定を行なった時（校正時）の発光素子への電流値 Ccal1, Ccal2 ,Ccal3と生体の減光度測定を行なった時（測定時）の発光素子への電流値 Cmeas1, Cmeas2, Cmeas3 が違うことがある。その場合入射光強度IiA1,IiA2 ,IiA3 は、
IiA1 ＝Ical1 × Cmeas1 ／Ccal1 （２６）
IiA2 ＝Ical2 × Cmeas2 ／Ccal2 （２７）
IiA3 ＝Ical3 × Cmeas3 ／Ccal3 （２８）
のように、発光素子への通電電流に応じて補正した値を用いる。
【００７３】
このステップにおいても、直流減光度Adc1,Adc2,Adc3の計算は、脈動する透過光強度It1,It2,It3 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。
【００７４】
＜ステップ４B ：一酸化炭素ヘモグロビン濃度の計算＞
次にステップ４B では、ＣＰＵ１４０は、ステップ２B で求めたΦ12、Φ13とステップ３B で求めた直流減光度Adc1、Adc2、Adc3を変数とする次式を用いて一酸化炭素ヘモグロビン濃度COHbdcを計算する。
C0Hbdc＝a6×Φ12＋b6×Φ13＋c6＋d6×Adc2／Adc1＋e6×Adc3／Adc1
（２９）
【００７５】
この式（２９）の係数a6〜e6はある母集団（例えばある選択された１０人のデータ）において、式（２９）を計算した結果のCOHbdcと、採血して測定した正確な一酸化炭素ヘモグロビン値COHbs の差が最小になるように最小二乗法で予め決めた値である。
【００７６】
ここで、本装置の効果を示すため、一酸化炭素ヘモグロビン濃度を、脈波の減光度変化分比Φ12、Φ13だけを使って計算した場合と、本装置により、Φ12、Φ13と直流減光度Adc1,Adc2,Adc3を使って計算した場合とを、それぞれ採血法によって求めた場合とを比較する。
【００７７】
Φ12、Φ13だけを使う場合、
COHbp ＝a7×Φ12＋b7×Φ13＋c7 （３０）
により計算する。この場合の係数a7,b7,c7も、上記母集団のデータを用いて上記と同様にして予め決定されたものである。式（３０）を用いて計算した一酸化炭素ヘモグロビン濃度COHbp と、採血法により測定した一酸化炭素ヘモグロビン濃度COHbs との相関を図８（ａ）に示す。
【００７８】
これに対し、直流減光度を使い、式（２９）により計算した本装置による一酸化炭素ヘモグロビン濃度C0Hbdcと、採血法により測定したCOHbs との相関を図８（ｂ）に示す。
【００７９】
これらの図を比較すると分かるように、計算式に直流減光度を入れることによって、採血法との一致が改善される。
【００８０】
また、上記の式（２９）の代わりに、次式を用いても良い。
COHbdc＝a8×Φ12＋b8×Φ13＋c8＋d8×Adc1＋e8×Adc2＋f8×Adc3
（３１）
この式の係数a8〜f8も上記と同様にして予め決定しておく。
【００８１】
ところで、直流減光度Adc1,Adc2,Adc3と直流透過光強度 DC1,DC2,DC3との間には逆相関の一定の関係があるので、減光度で計算しなくても、直流透過光強度をそのまま使って計算すれば同様の効果が得られる。この場合、直流透過光強度はプローブの入射光強度に比例するので、入射光強度を、発光素子（ＬＥＤ２０ａ，２０ｂ，２０ｃ）の特定の通電電流値で規格化した状態で補正する必要が有る。補正した各波長の直流透過光強度Itcomp1,Itcomp2,Itcomp3 は次式により計算する。
Itcomp1 ＝DC1 ×(Iical1 ／Istd1)×(Cmeas1 ／Ccal1) （３２）
Itcomp2 ＝DC2 ×(Iical2 ／Istd2)×(Cmeas2 ／Ccal2) （３３）
Itcomp3 ＝DC3 ×(Iical3 ／Istd3)×(Cmeas3 ／Ccal3) （３４）
DC1,DC2,DC3:測定した透過光強度
Iical1,Iical2,Iical3: 校正時に計算した入射光強度
Istd1,Istd2,Istd3:標準とする入射光強度
Cmeas1,Cmeas2,Cmeas3: 生体を測定したときの発光素子の通電電流
Ccal1,Ccal2,Ccal3:入射光強度校正時の発光素子の通電電流
【００８２】
したがって、上記の式（２９）、（３１）の代わりに、次式を用いても良い。
COHbdc＝a9×Φ12＋b9×Φ13＋c9＋d9×Itcomp1 ／Itcomp2 ＋e9×Itcomp1 ／Itcomp3 （３５）
COHbdc＝a10 ×Φ12＋b10 ×Φ13＋ c10＋ d10×Itcomp1 ＋e10 ×Itcomp2 ＋f10 ×Itcomp3 （３６）
式（３５）における係数a9〜e9、式（３６）における係数a10 〜f10 は、それぞれ上記と同様の母集団のデータを用いて、上記と同様にして予め決定されるものである。
【００８３】
このステップにおいても、COHbdcの計算は、脈動する透過光強度It1,It2,It3 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。そして、ＣＰＵ１４０は、求めたCOHbdcをメモリ１５０に格納すると共に、表示部１６０に表示する。
【００８４】
本実施の形態の装置によれば、血液に光吸収がある第１波長９４０ｎｍ、第２波長６６０ｎｍおよび第３波長６２０ｎｍの直流減光度または直流透過光強度を用いて計算したので、測定部位の血液層全体の厚さが考慮され、精度良くヘモグロビン濃度を測定することができる。
【００８５】
［第３の実施の形態：色素濃度測定装置］
次に、本発明の第３の実施の形態を説明する。この第３の実施の形態は、色素濃度測定装置である。この装置は、第２の実施の形態と同様に３波長の光を用いるので、装置の構成ブロック図は図６と同じである。このため、本装置の構成ブロック図は省略し、以下の説明では図６を利用する。しかし、ＣＰＵ１４０が行なう処理のプログラムは、図９に示すように第２の実施の形態と異なる。また、第１の波長は９４０nm、第２の波長は６６０nm、第３の波長は８０５nmである。以下、図９を参照して本装置の動作を説明する。
【００８６】
＜ステップ１C ：プローブの入射光強度の測定（校正）＞
このステップ１C では、生体３０への入射光強度を測定する。上記の第２の実施の形態の装置におけるステップ１B と同様にして、ファントム３０Ａを用い、３波長の入射光強度Iical1、Iical2、Iical3を求め、これらをメモリ１５０に格納する。
【００８７】
＜ステップ２C ：脈波の測定と減光度変化分比Φの計算＞
次にステップ２C では、操作者は、プローブを生体３０に装着した後、色素を注入する。色素として例えばインドシアニングリーンを用いる。そして、上記第２の実施の形態におけるステップ２B と同様にして各波長の減光度の変化分の相互の比を求める。すなわち、次の式を計算してΦ12、Φ13を求め、その結果をメモリ１５０に格納する。
Φ12＝ΔA1／ΔA2＝ (AC1 ／DC1)／ (AC2 ／DC2) （２１）
Φ13＝ΔA1／ΔA3＝ (AC1 ／DC1)／ (AC3 ／DC3) （２２）
Φ12：第１波長と第２波長の減光度A1,A2 の変化分ΔA1, ΔA2の相互比
Φ13：第１波長と第３波長の減光度A1,A3 の変化分ΔA1, ΔA3の相互比
DC1,DC2,DC3:第１〜３波長の透過光強度It1,It2,It3 における直流透過光強度
AC1,AC2,AC3:第１〜３波長の透過光強度It1,It2,It3 における交流透過光強度
【００８８】
このステップにおいて、Φ12、Φ13の計算は、脈動する透過光強度It1,It2,Ii3 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。
【００８９】
＜ステップ３C ：直流減光度の計算＞
次にステップ３C では、各波長の直流減光度Adc1,Adc2,Adc3を求める。ここで、ＣＰＵ１４０は、第２の実施の形態のステップ３B と同様にして、各波長の入射光強度Iical1,Iical2,Iical3と、各波長の直流透過光強度DC1,DC2,DC3 と次式により、各波長の直流減光度Adc1, Adc2, Adc3を求め、その結果をメモリ１５０に格納する。
Adc1＝ log(Iical1 ／DC1)＝logIical1 −logDC1 （２３）
Adc2＝ log(Iical2 ／DC2)＝logIical2 −logDC2 （２４）
Adc3＝ log(Iical3 ／DC3)＝logIical3 −logDC3 （２５）
【００９０】
本装置が、発光素子（ＬＥＤ２０ａ，２０ｂ，２０ｃ）に流す電流を変化させない回路構成を持つ場合、上記のようにステップ１C で求めた入射光強度Iical1,Iical2,Iical3をそのまま用いる。
【００９１】
しかし、本装置が、発光素子の電流を変化させて最適の透過光強度が得られるように調整する回路構成を持つ場合、ファントム３０Ａを用いて入射光強度測定を行なった時（校正時）の発光素子への電流値 Ccal1, Ccal2 ,Ccal3と生体の減光度測定を行なった時（測定時）の発光素子への電流値 Cmeas1, Cmeas2, Cmeas3 が違うことがある。その場合入射光強度IiA1,IiA2 ,IiA3 は、
IiA1 ＝Ical1 × Cmeas1 ／Ccal1 （２６）
IiA2 ＝Ical2 × Cmeas2 ／Ccal2 （２７）
IiA3 ＝Ical3 × Cmeas3 ／Ccal3 （２８）
のように、発光素子への通電電流に応じて補正した値を用いる。
【００９２】
このステップにおいても、直流減光度Adc1,Adc2,Adc3の計算は、脈動する透過光強度It1,It2,It3 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。
【００９３】
＜ステップ４C ：色素濃度の計算＞
次にステップ４C では、ＣＰＵ１４０は、上記ステップ２C で求めたΦ12、Φ13と上記ステップ３C で求めた直流減光度Adc1、Adc2、Adc3を変数とする次式を用いて色素濃度Cdを計算する。
Cddc＝a11 ×Φ12＋b11 ×Φ13＋c11 ＋d11 ×Adc2／Adc1＋e11 ×Adc3／Adc1（３７）
【００９４】
この式（３７）の係数a11 〜e11 はある母集団（例えばある選択された１０人のデータ）において、式（３７）を計算した結果であるCddcと、採血して測定した正確なCds 値の差が最小になるように最小二乗法で予め決めた値である。
ここで、本装置の効果を示すため、色素濃度を、脈波の減光度変化分比Φ12、Φ13だけを使って計算した場合と、本装置によるΦ12、Φ13とAdc1,Adc2,Adc3を使って計算した場合とを、それぞれ採血法によって求めた場合とを比較する。
【００９５】
Φ12、Φ13だけを使う場合、
Cdp ＝a12 ×Φ12＋b12 ×Φ13＋c12 （３８）
により計算する。この場合の係数a12,b12,c12 も、上記母集団のデータを用いて上記と同様にして予め決定されたものである。式（３８）を用いて計算した色素濃度Cdp と、採血法により測定した色素濃度Cds との相関を図１０（ａ）に示す。これに対し、直流減光度を使い、式（３７）により計算した本装置による色素濃度Cddcと、採血法により測定した色素濃度Cds との相関を図１０（ｂ）に示す。
【００９６】
これらの図を比較すると分かるように、計算式に直流減光度を入れることによって、採血法との一致が改善される。
また、上記の式（３７）の代わりに、次式を用いても良い。
Cddc＝a13 ×Φ12＋b13 ×Φ13＋c13 ＋d13 ×Adc1＋ e13×Adc2
＋f12 ×Adc3 （３９）
この式の係数a13 〜f13 も上記と同様にして予め決定しておく。
【００９７】
ところで、直流減光度Adc1,Adc2,Adc3と直流透過光強度 DC1,DC2,DC3との間には逆相関の一定の関係があるので、減光度で計算しなくても、直流透過光強度をそのまま使っても同様の効果が得られる。この場合、直流透過光強度はプローブの入射光強度に比例するので、入射光強度を、発光素子（ＬＥＤ２０ａ，２０ｂ，２０ｃ）の特定の通電電流値で規格化した状態で補正する必要が有る。補正した各波長の直流透過光強度Itcompは次式により計算する。
Itcomp1 ＝DC1 ×(Iical1 ／Istd1)×(Cmeas1 ／Ccal1) （４０）
Itcomp2 ＝DC2 ×(Iical2 ／Istd2)×(Cmeas2 ／Ccal2) （４１）
Itcomp3 ＝DC3 ×(Iical3 ／Istd3)×(Cmeas3 ／Ccal3) （４２）
DC1,DC2,DC3:測定した透過光強度
Iical1,Iical2,Iical3: 校正時に計算した入射光強度
Istd1,Istd2,Istd3:標準とする入射光強度
Cmeas1,Cmeas2,Cmeas3: 生体を測定したときの発光素子の通電電流
Ccal1,Ccal2,Ccal3:入射光強度校正時の発光素子の通電電流
【００９８】
したがって、上記の式（３７）、（３９）の代わりに、次式を用いても良い。
Cddc＝a14 ×Φ12＋ b14×Φ13＋ c14＋ d14×Itcomp1 ／Itcomp2 ＋ e14×Itcomp1 ／Itcomp3 （４３）
Cddc＝a15×Φ12＋ b15×Φ13＋c15＋d15×Itcomp１＋e15 ×Itcomp2 ＋f15 ×Itcomp3 （４４）
【００９９】
ここで、式（４３）における係数a14 〜e14 、式（４４）における係数a15 〜f15 は、それぞれ上記と同様の母集団のデータを用いて、上記と同様にして予め決定されるものである。
【０１００】
このステップにおいても、Cddcの計算は、脈動する透過光強度It1,It2,It3 の波の１つずつについて、すなわち心拍の一拍ずつについて行なう。そして、ＣＰＵ１４０は、求めたCddcをメモリ１５０に格納すると共に、表示部１６０に表示する。
【０１０１】
本実施の形態の装置によれば、血液に光吸収がある第１波長９４０ｎｍ、第２波長８０５ｎｍおよび第３波長６６０ｎｍの直流減光度または直流透過光強度を用いて計算したので、測定部位の血液層全体の厚さが考慮され、精度良く色素濃度を測定することができる。
【０１０２】
なお、以上の説明においては、発光素子に流れる電流を制御して生体に照射する光強度を変える場合について説明したが、これは、発光素子に印加する電圧を制御して発光素子の光強度を変え、測定される入射光強度や直流透過光強度をその電圧値に応じて補正するようにしても良い。
【０１０３】
【発明の効果】
本発明によれば、パルスフォトメトリーを原理とする血中吸光物質濃度測定において、精度良く測定することができる。
【図面の簡単な説明】
【図１】第１の実施の形態のヘモグロビン濃度測定装置の動作を説明するための図。
【図２】第１の実施の形態のヘモグロビン濃度測定装置のブロック構成図。
【図３】生体へ光を照射したときの入射光強度と透過光強度との関係を説明する図。
【図４】直流透過光強度を説明する図。
【図５】第１の実施の形態のヘモグロビン濃度測定装置の効果を説明する図。
【図６】第２の実施の形態の一酸化炭素ヘモグロビン濃度測定装置のブロック構成図。
【図７】第２の実施の形態の一酸化炭素ヘモグロビン濃度測定装置の動作を説明するための図。
【図８】第２の実施の形態の一酸化炭素ヘモグロビン濃度測定装置の効果を説明する図。
【図９】第３の実施の形態の色素濃度測定装置の動作を説明するための図。
【図１０】第３の実施の形態の色素濃度測定装置の効果を説明する図。
【符号の説明】
１、１００ 発光部
５、５０ 受光部
１４、１４０ ＣＰＵ
Claims (11)
- 複数種の波長の光を発する光源を持ち、これらの光を生体組織に照射する光照射手段と、
この光照射手段から発せられた各波長の光の生体組織への入射光強度を測定する入射光強度測定手段と、
前記生体組織を透過した各波長の光を受光し、それぞれの透過光強度を測定する透過光強度測定手段と、
この透過光強度測定手段が測定した透過光強度に基づいて、血液の厚み変化により発生する各波長の減光度変化分の相互の比を演算する第１の演算手段と、
この第１の演算手段が求めた比と、前記入射光強度測定手段が測定した入射光強度と、前記透過光強度測定手段が測定した透過光強度とに基づいて血液中の吸光物質濃度を演算して求める第２の演算手段と
を具備する血中吸光物質濃度測定装置。 - 前記第２の演算手段は、
各波長について、前記入射光強度測定手段が測定した入射光強度と、前記透過光強度測定手段が測定した透過光強度とに基づいて前記生体組織の直流減光度を演算して求め、この直流減光度と前記第１の演算手段が求めた比とに基づいて血液中の吸光物質濃度を演算して求めることを特徴とする請求項１に記載の血中吸光物質濃度測定装置。 - 前記第２の演算手段は、各波長の直流減光度の相互の比を求め、この比と前記第１の演算手段が求めた比とに基づいて血液中の吸光物質濃度を演算して求めることを特徴とする請求項２に記載の血中吸光物質濃度測定装置。
- 前記第２の演算手段は、
各波長について、前記入射光強度測定手段が測定した入射光強度と、前記透過光強度測定手段が測定した透過光強度とに基づいて前記生体組織の直流透過光強度を演算して求め、この直流透過光強度と前記第１の演算手段が求めた比とに基づいて血液中の吸光物質濃度を演算して求めることを特徴ととする請求項１に記載の血中吸光物質濃度測定装置。 - 前記第２の演算手段は、各波長の直流透過光強度の相互の比を求め、この比と前記第１の演算手段が求めた比とに基づいて血液中の吸光物質濃度を演算して求めることを特徴とする請求項４に記載の血中吸光物質濃度測定装置。
- 前記光照射手段は、前記光源が発光素子から成り、この発光素子に供給する電流または印加する電圧を制御する制御部を備え、
前記第２の演算手段は前記制御部により制御された電流または電圧に応じて、測定した入射光強度を補正し、この補正した入射光強度を用いて演算することを特徴とする請求項１乃至５のいずれか１に記載の血中吸光物質濃度測定装置。 - 前記入射光強度測定手段は、前記光照射手段と前記透過光強度測定手段との間の所定位置に、減光度が既知のファントムが配置された場合に、前記透過光強度測定手段により測定される各波長の透過光強度に基づいて入射光強度を演算して求めることを特徴とする請求項１乃至６のいずれか１に記載の血中吸光物質濃度測定装置。
- 前記入射光強度測定手段は、前記ファントムが前記所定位置に配置されたか否かを検出するセンサを有し、前記ファントムが前記所定位置に配置されたことを前記センサにより検出されると、入射光強度を演算して求めることを特徴とする請求項７に記載の血中吸光物質濃度測定装置。
- ｎ種の波長の光を発する光源を持ち、これらの光を生体組織に照射する光照射手段と、
この光照射手段から発せられたｎ種の波長の光の生体組織への入射光強度を測定する入射光強度測定手段と、
前記生体組織を透過したｎ種の波長の光を受光し、それぞれの透過光強度を測定する透過光強度測定手段と、
この透過光強度測定手段が測定した透過光強度に基づいて、血液の厚み変化により発生する各波長の減光度変化分の相互の比のうち最大でｎ−１個の比を演算する第１の演算手段と、
この第１の演算手段が求めた最大でｎ−１個の比と、前記入射光強度測定手段が求めた入射光強度と、前記透過光強度測定手段が求めた透過光強度とに基づいて、各波長に対応する最大でｎ個の直流減光度を演算して求め、これらの演算結果に基づいて血液中の吸光物質濃度を演算して求める第２の演算手段とを具備する血中吸光物質濃度測定装置。 - ｎ種の波長の光を発する光源を持ち、これらの光を生体組織に照射する光照射手段と、
この光照射手段から発せられたｎ種の波長の光の生体組織への入射光強度を測定する入射光強度測定手段と、
前記生体組織を透過したｎ種の波長の光を受光し、それぞれの透過光強度を測定する透過光強度測定手段と、
この透過光強度測定手段が測定した透過光強度に基づいて、血液の厚み変化により発生する各波長の減光度変化分の相互の比のうち最大でｎ−１個の比を演算する第１の演算手段と、
この第１の演算手段が求めた最大でｎ−１個の比と、前記入射光強度測定手段が求めた入射光強度と、前記透過光強度測定手段が求めた透過光強度とに基づいて、各波長に対応する最大でｎ個の直流透過光強度を演算して求め、これらの演算結果に基づいて血液中の吸光物質濃度を演算して求める第２の演算手段とを具備する血中吸光物質濃度測定装置。 - 前記光源は、血液に光吸収がある波長の光および／または血液以外の生体組織に光吸収がある波長の光を発することを特徴とする請求項１乃至１０のいずれか１に記載の血中吸光物質濃度測定装置。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003023990A JP4284674B2 (ja) | 2003-01-31 | 2003-01-31 | 血中吸光物質濃度測定装置 |
| DE200410005086 DE102004005086B4 (de) | 2003-01-31 | 2004-02-02 | Vorrichtung zum Messen der Konzentration einer lichtabsorbierenden Substanz in Blut |
| US10/768,680 US7257433B2 (en) | 2003-01-31 | 2004-02-02 | Apparatus for measuring concentration of light-absorbing substance in blood |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003023990A JP4284674B2 (ja) | 2003-01-31 | 2003-01-31 | 血中吸光物質濃度測定装置 |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2004230000A JP2004230000A (ja) | 2004-08-19 |
| JP2004230000A5 JP2004230000A5 (ja) | 2006-02-23 |
| JP4284674B2 true JP4284674B2 (ja) | 2009-06-24 |
Family
ID=32709270
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2003023990A Expired - Fee Related JP4284674B2 (ja) | 2003-01-31 | 2003-01-31 | 血中吸光物質濃度測定装置 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US7257433B2 (ja) |
| JP (1) | JP4284674B2 (ja) |
| DE (1) | DE102004005086B4 (ja) |
Families Citing this family (108)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6018673A (en) | 1996-10-10 | 2000-01-25 | Nellcor Puritan Bennett Incorporated | Motion compatible sensor for non-invasive optical blood analysis |
| US9042952B2 (en) | 1997-01-27 | 2015-05-26 | Lawrence A. Lynn | System and method for automatic detection of a plurality of SPO2 time series pattern types |
| US8932227B2 (en) | 2000-07-28 | 2015-01-13 | Lawrence A. Lynn | System and method for CO2 and oximetry integration |
| US20060161071A1 (en) | 1997-01-27 | 2006-07-20 | Lynn Lawrence A | Time series objectification system and method |
| US9521971B2 (en) | 1997-07-14 | 2016-12-20 | Lawrence A. Lynn | System and method for automatic detection of a plurality of SPO2 time series pattern types |
| US20070191697A1 (en) | 2006-02-10 | 2007-08-16 | Lynn Lawrence A | System and method for SPO2 instability detection and quantification |
| US6675031B1 (en) | 1999-04-14 | 2004-01-06 | Mallinckrodt Inc. | Method and circuit for indicating quality and accuracy of physiological measurements |
| US20090281838A1 (en) | 2008-05-07 | 2009-11-12 | Lawrence A. Lynn | Medical failure pattern search engine |
| US20060195041A1 (en) | 2002-05-17 | 2006-08-31 | Lynn Lawrence A | Centralized hospital monitoring system for automatically detecting upper airway instability and for preventing and aborting adverse drug reactions |
| US9053222B2 (en) | 2002-05-17 | 2015-06-09 | Lawrence A. Lynn | Patient safety processor |
| US8135448B2 (en) * | 2001-03-16 | 2012-03-13 | Nellcor Puritan Bennett Llc | Systems and methods to assess one or more body fluid metrics |
| US7239902B2 (en) | 2001-03-16 | 2007-07-03 | Nellor Puritan Bennett Incorporated | Device and method for monitoring body fluid and electrolyte disorders |
| US6591122B2 (en) * | 2001-03-16 | 2003-07-08 | Nellcor Puritan Bennett Incorporated | Device and method for monitoring body fluid and electrolyte disorders |
| US7657292B2 (en) * | 2001-03-16 | 2010-02-02 | Nellcor Puritan Bennett Llc | Method for evaluating extracellular water concentration in tissue |
| US6757554B2 (en) * | 2001-05-22 | 2004-06-29 | Alfred E. Mann Institute For Biomedical Engineering At The University Of Southern California | Measurement of cardiac output and blood volume by non-invasive detection of indicator dilution |
| US6754516B2 (en) | 2001-07-19 | 2004-06-22 | Nellcor Puritan Bennett Incorporated | Nuisance alarm reductions in a physiological monitor |
| US7006856B2 (en) | 2003-01-10 | 2006-02-28 | Nellcor Puritan Bennett Incorporated | Signal quality metrics design for qualifying data for a physiological monitor |
| US7016715B2 (en) | 2003-01-13 | 2006-03-21 | Nellcorpuritan Bennett Incorporated | Selection of preset filter parameters based on signal quality |
| US7120479B2 (en) | 2004-02-25 | 2006-10-10 | Nellcor Puritan Bennett Inc. | Switch-mode oximeter LED drive with a single inductor |
| US7190985B2 (en) | 2004-02-25 | 2007-03-13 | Nellcor Puritan Bennett Inc. | Oximeter ambient light cancellation |
| US7534212B2 (en) | 2004-03-08 | 2009-05-19 | Nellcor Puritan Bennett Llc | Pulse oximeter with alternate heart-rate determination |
| US7194293B2 (en) | 2004-03-08 | 2007-03-20 | Nellcor Puritan Bennett Incorporated | Selection of ensemble averaging weights for a pulse oximeter based on signal quality metrics |
| US7277741B2 (en) | 2004-03-09 | 2007-10-02 | Nellcor Puritan Bennett Incorporated | Pulse oximetry motion artifact rejection using near infrared absorption by water |
| US7937128B2 (en) * | 2004-07-09 | 2011-05-03 | Masimo Corporation | Cyanotic infant sensor |
| WO2006094107A1 (en) | 2005-03-01 | 2006-09-08 | Masimo Laboratories, Inc. | Physiological parameter confidence measure |
| US7392075B2 (en) | 2005-03-03 | 2008-06-24 | Nellcor Puritan Bennett Incorporated | Method for enhancing pulse oximetry calculations in the presence of correlated artifacts |
| US7725146B2 (en) | 2005-09-29 | 2010-05-25 | Nellcor Puritan Bennett Llc | System and method for pre-processing waveforms |
| US7725147B2 (en) | 2005-09-29 | 2010-05-25 | Nellcor Puritan Bennett Llc | System and method for removing artifacts from waveforms |
| US20070106126A1 (en) | 2005-09-30 | 2007-05-10 | Mannheimer Paul D | Patient monitoring alarm escalation system and method |
| US8929973B1 (en) * | 2005-10-24 | 2015-01-06 | Lockheed Martin Corporation | Apparatus and method for characterizing optical sources used with human and animal tissues |
| WO2007084632A2 (en) * | 2006-01-20 | 2007-07-26 | Alfred E. Mann Institute For Biomedical Engineering At The University Of Southern California | Measurement of cardiac output and blood volume by non-invasive detection of indicator dilution |
| US7668579B2 (en) | 2006-02-10 | 2010-02-23 | Lynn Lawrence A | System and method for the detection of physiologic response to stimulation |
| JP2009527333A (ja) | 2006-02-22 | 2009-07-30 | ヴィヴァム ネクサス エルエルシー | 検体測定装置及び方法 |
| US8702606B2 (en) | 2006-03-21 | 2014-04-22 | Covidien Lp | Patient monitoring help video system and method |
| US8380271B2 (en) | 2006-06-15 | 2013-02-19 | Covidien Lp | System and method for generating customizable audible beep tones and alarms |
| US8064975B2 (en) | 2006-09-20 | 2011-11-22 | Nellcor Puritan Bennett Llc | System and method for probability based determination of estimated oxygen saturation |
| US8696593B2 (en) | 2006-09-27 | 2014-04-15 | Covidien Lp | Method and system for monitoring intracranial pressure |
| US7922665B2 (en) | 2006-09-28 | 2011-04-12 | Nellcor Puritan Bennett Llc | System and method for pulse rate calculation using a scheme for alternate weighting |
| US7925511B2 (en) | 2006-09-29 | 2011-04-12 | Nellcor Puritan Bennett Llc | System and method for secure voice identification in a medical device |
| US20080081956A1 (en) | 2006-09-29 | 2008-04-03 | Jayesh Shah | System and method for integrating voice with a medical device |
| US7706896B2 (en) | 2006-09-29 | 2010-04-27 | Nellcor Puritan Bennett Llc | User interface and identification in a medical device system and method |
| US7848891B2 (en) | 2006-09-29 | 2010-12-07 | Nellcor Puritan Bennett Llc | Modulation ratio determination with accommodation of uncertainty |
| US8116852B2 (en) * | 2006-09-29 | 2012-02-14 | Nellcor Puritan Bennett Llc | System and method for detection of skin wounds and compartment syndromes |
| US20080097175A1 (en) * | 2006-09-29 | 2008-04-24 | Boyce Robin S | System and method for display control of patient monitor |
| US8160668B2 (en) | 2006-09-29 | 2012-04-17 | Nellcor Puritan Bennett Llc | Pathological condition detector using kernel methods and oximeters |
| US7698002B2 (en) | 2006-09-29 | 2010-04-13 | Nellcor Puritan Bennett Llc | Systems and methods for user interface and identification in a medical device |
| US8068890B2 (en) | 2006-09-29 | 2011-11-29 | Nellcor Puritan Bennett Llc | Pulse oximetry sensor switchover |
| US8068891B2 (en) | 2006-09-29 | 2011-11-29 | Nellcor Puritan Bennett Llc | Symmetric LED array for pulse oximetry |
| US8175667B2 (en) | 2006-09-29 | 2012-05-08 | Nellcor Puritan Bennett Llc | Symmetric LED array for pulse oximetry |
| US8728059B2 (en) | 2006-09-29 | 2014-05-20 | Covidien Lp | System and method for assuring validity of monitoring parameter in combination with a therapeutic device |
| US8175665B2 (en) * | 2007-03-09 | 2012-05-08 | Nellcor Puritan Bennett Llc | Method and apparatus for spectroscopic tissue analyte measurement |
| US20080221418A1 (en) * | 2007-03-09 | 2008-09-11 | Masimo Corporation | Noninvasive multi-parameter patient monitor |
| US8265724B2 (en) | 2007-03-09 | 2012-09-11 | Nellcor Puritan Bennett Llc | Cancellation of light shunting |
| US8280469B2 (en) | 2007-03-09 | 2012-10-02 | Nellcor Puritan Bennett Llc | Method for detection of aberrant tissue spectra |
| US20080221416A1 (en) * | 2007-03-09 | 2008-09-11 | Nellcor Puritan Bennett Llc | System and method for detection of macular degeneration using spectrophotometry |
| US8781544B2 (en) | 2007-03-27 | 2014-07-15 | Cercacor Laboratories, Inc. | Multiple wavelength optical sensor |
| US8374665B2 (en) | 2007-04-21 | 2013-02-12 | Cercacor Laboratories, Inc. | Tissue profile wellness monitor |
| JP4569615B2 (ja) * | 2007-09-25 | 2010-10-27 | ブラザー工業株式会社 | 印刷装置 |
| JP5278854B2 (ja) | 2007-12-10 | 2013-09-04 | 富士フイルム株式会社 | 画像処理システムおよびプログラム |
| US8204567B2 (en) | 2007-12-13 | 2012-06-19 | Nellcor Puritan Bennett Llc | Signal demodulation |
| US8092993B2 (en) | 2007-12-31 | 2012-01-10 | Nellcor Puritan Bennett Llc | Hydrogel thin film for use as a biosensor |
| US20090171174A1 (en) * | 2007-12-31 | 2009-07-02 | Nellcor Puritan Bennett Llc | System and method for maintaining battery life |
| US8750953B2 (en) | 2008-02-19 | 2014-06-10 | Covidien Lp | Methods and systems for alerting practitioners to physiological conditions |
| US8275553B2 (en) | 2008-02-19 | 2012-09-25 | Nellcor Puritan Bennett Llc | System and method for evaluating physiological parameter data |
| US8140272B2 (en) | 2008-03-27 | 2012-03-20 | Nellcor Puritan Bennett Llc | System and method for unmixing spectroscopic observations with nonnegative matrix factorization |
| US8437822B2 (en) | 2008-03-28 | 2013-05-07 | Covidien Lp | System and method for estimating blood analyte concentration |
| US8364224B2 (en) | 2008-03-31 | 2013-01-29 | Covidien Lp | System and method for facilitating sensor and monitor communication |
| US8112375B2 (en) | 2008-03-31 | 2012-02-07 | Nellcor Puritan Bennett Llc | Wavelength selection and outlier detection in reduced rank linear models |
| US8292809B2 (en) | 2008-03-31 | 2012-10-23 | Nellcor Puritan Bennett Llc | Detecting chemical components from spectroscopic observations |
| USD626561S1 (en) | 2008-06-30 | 2010-11-02 | Nellcor Puritan Bennett Llc | Circular satseconds indicator and triangular saturation pattern detection indicator for a patient monitor display panel |
| US9895068B2 (en) | 2008-06-30 | 2018-02-20 | Covidien Lp | Pulse oximeter with wait-time indication |
| USD626562S1 (en) | 2008-06-30 | 2010-11-02 | Nellcor Puritan Bennett Llc | Triangular saturation pattern detection indicator for a patient monitor display panel |
| US8862194B2 (en) | 2008-06-30 | 2014-10-14 | Covidien Lp | Method for improved oxygen saturation estimation in the presence of noise |
| US20100016731A1 (en) * | 2008-07-15 | 2010-01-21 | Cardiox Corporation | Hemodynamic Detection of Circulatory Anomalies |
| US8433382B2 (en) | 2008-09-30 | 2013-04-30 | Covidien Lp | Transmission mode photon density wave system and method |
| US8386000B2 (en) | 2008-09-30 | 2013-02-26 | Covidien Lp | System and method for photon density wave pulse oximetry and pulse hemometry |
| US8417309B2 (en) | 2008-09-30 | 2013-04-09 | Covidien Lp | Medical sensor |
| US8968193B2 (en) | 2008-09-30 | 2015-03-03 | Covidien Lp | System and method for enabling a research mode on physiological monitors |
| US9078609B2 (en) * | 2008-10-02 | 2015-07-14 | Nellcor Puritan Bennett Ireland | Extraction of physiological measurements from a photoplethysmograph (PPG) signal |
| US8221319B2 (en) | 2009-03-25 | 2012-07-17 | Nellcor Puritan Bennett Llc | Medical device for assessing intravascular blood volume and technique for using the same |
| US8509869B2 (en) | 2009-05-15 | 2013-08-13 | Covidien Lp | Method and apparatus for detecting and analyzing variations in a physiologic parameter |
| US8494786B2 (en) | 2009-07-30 | 2013-07-23 | Covidien Lp | Exponential sampling of red and infrared signals |
| US8494606B2 (en) | 2009-08-19 | 2013-07-23 | Covidien Lp | Photoplethysmography with controlled application of sensor pressure |
| US8788001B2 (en) | 2009-09-21 | 2014-07-22 | Covidien Lp | Time-division multiplexing in a multi-wavelength photon density wave system |
| US8704666B2 (en) | 2009-09-21 | 2014-04-22 | Covidien Lp | Medical device interface customization systems and methods |
| US8494604B2 (en) | 2009-09-21 | 2013-07-23 | Covidien Lp | Wavelength-division multiplexing in a multi-wavelength photon density wave system |
| US8798704B2 (en) | 2009-09-24 | 2014-08-05 | Covidien Lp | Photoacoustic spectroscopy method and system to discern sepsis from shock |
| US9554739B2 (en) | 2009-09-29 | 2017-01-31 | Covidien Lp | Smart cable for coupling a medical sensor to an electronic patient monitor |
| US8376955B2 (en) | 2009-09-29 | 2013-02-19 | Covidien Lp | Spectroscopic method and system for assessing tissue temperature |
| US8515511B2 (en) | 2009-09-29 | 2013-08-20 | Covidien Lp | Sensor with an optical coupling material to improve plethysmographic measurements and method of using the same |
| 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 |
| US8930145B2 (en) | 2010-07-28 | 2015-01-06 | Covidien Lp | Light focusing continuous wave photoacoustic spectroscopy and its applications to patient monitoring |
| WO2012032536A2 (en) * | 2010-09-06 | 2012-03-15 | Abhishek Sen | System and method for non-invasive determination of hemoglobin concentration in blood |
| CA2833258A1 (en) | 2011-04-15 | 2012-10-18 | Constitution Medical, Inc. | Measuring volume and constituents of cells |
| US9833146B2 (en) | 2012-04-17 | 2017-12-05 | Covidien Lp | Surgical system and method of use of the same |
| WO2015137151A1 (ja) * | 2014-03-11 | 2015-09-17 | 公立大学法人横浜市立大学 | 支援用具、および生体情報取得システム |
| JP6385865B2 (ja) | 2014-03-28 | 2018-09-05 | 日本光電工業株式会社 | パルスフォトメータ |
| CN104068865B (zh) * | 2014-04-24 | 2017-05-31 | 辛勤 | 一种测量血氧饱和度的方法以及便携式设备 |
| US10524704B2 (en) | 2015-03-04 | 2020-01-07 | Nihon Kohden Corporation | Pulse photometer and method for calculating concentration of light absorber in blood |
| JP6713252B2 (ja) | 2015-03-30 | 2020-06-24 | 日本光電工業株式会社 | 生体情報測定システム |
| KR102109245B1 (ko) * | 2015-05-12 | 2020-05-11 | 코니카 미놀타 가부시키가이샤 | 식물 생육 지표 측정 장치 및 그의 방법 그리고 식물 생육 지표 측정 시스템 |
| CN109313115B (zh) * | 2016-05-11 | 2022-04-12 | 诺尔生物医药有限公司 | 全血so2传感器 |
| WO2019158385A1 (en) | 2018-02-13 | 2019-08-22 | Robert Bosch Gmbh | A controller and method to determine incident intensity of a light source and total hemoglobin |
| WO2020260365A1 (en) | 2019-06-24 | 2020-12-30 | Robert Bosch Gmbh | A bioanalyte monitoring device |
| WO2021008933A1 (en) | 2019-07-16 | 2021-01-21 | Robert Bosch Gmbh | A bioanalyte monitoring device |
| CN113456067B (zh) * | 2021-06-24 | 2023-06-23 | 江西科莱富健康科技有限公司 | 体外血红蛋白浓度测量方法、系统及计算机设备 |
| CN113588579B (zh) * | 2021-08-23 | 2024-03-01 | 重庆山外山血液净化技术股份有限公司 | 血液离子浓度探测装置和方法、钙离子浓度探测方法 |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5326437B2 (ja) | 1974-03-29 | 1978-08-02 | ||
| JPH02164341A (ja) | 1988-12-19 | 1990-06-25 | Nippon Koden Corp | ヘモグロビン濃度測定装置 |
| JP2608828B2 (ja) | 1992-02-06 | 1997-05-14 | 日本光電工業株式会社 | 非観血式オキシメータ |
| US5385143A (en) * | 1992-02-06 | 1995-01-31 | Nihon Kohden Corporation | Apparatus for measuring predetermined data of living tissue |
| JP3116252B2 (ja) * | 1992-07-09 | 2000-12-11 | 日本光電工業株式会社 | パルスオキシメータ |
| US5490523A (en) * | 1994-06-29 | 1996-02-13 | Nonin Medical Inc. | Finger clip pulse oximeter |
| DE19612425C2 (de) * | 1995-03-31 | 2000-08-31 | Nihon Kohden Corp | Apparat zur Messung von Hämoglobinkonzentration |
| US6070093A (en) * | 1997-12-02 | 2000-05-30 | Abbott Laboratories | Multiplex sensor and method of use |
| US6415236B2 (en) | 1999-11-30 | 2002-07-02 | Nihon Kohden Corporation | Apparatus for determining concentrations of hemoglobins |
| JP2002228579A (ja) | 1999-11-30 | 2002-08-14 | Nippon Koden Corp | ヘモグロビン濃度測定装置 |
| KR100612827B1 (ko) * | 2001-04-19 | 2006-08-14 | 삼성전자주식회사 | 비 침습적인 헤모글로빈 농도와 산소 포화도 모니터링방법 및 장치 |
| JP4603718B2 (ja) | 2001-04-26 | 2010-12-22 | 株式会社日立メディコ | 生体光計測装置 |
-
2003
- 2003-01-31 JP JP2003023990A patent/JP4284674B2/ja not_active Expired - Fee Related
-
2004
- 2004-02-02 US US10/768,680 patent/US7257433B2/en active Active
- 2004-02-02 DE DE200410005086 patent/DE102004005086B4/de not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| DE102004005086A1 (de) | 2004-08-12 |
| JP2004230000A (ja) | 2004-08-19 |
| US7257433B2 (en) | 2007-08-14 |
| US20040176670A1 (en) | 2004-09-09 |
| DE102004005086B4 (de) | 2013-04-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4284674B2 (ja) | 血中吸光物質濃度測定装置 | |
| US11647923B2 (en) | Tissue profile wellness monitor | |
| JP3925945B2 (ja) | 検体を傷つけずに、血液が供給されている組織の酸素飽和量を測定する方法 | |
| US8606342B2 (en) | Pulse and active pulse spectraphotometry | |
| US8632472B2 (en) | Method and apparatus for detection of venous pulsation | |
| CA2558643C (en) | Pulse oximetry motion artifact rejection using near infrared absorption by water | |
| US5891024A (en) | Two stage calibration and analyte measurement scheme for spectrophotomeric analysis | |
| US8515514B2 (en) | Compensation of human variability in pulse oximetry | |
| US6968221B2 (en) | Low-cost method and apparatus for non-invasively measuring blood glucose levels | |
| US7254432B2 (en) | Method and device for non-invasive measurements of blood parameters | |
| US7952692B2 (en) | Method and apparatus for determination of analyte concentration | |
| US6963767B2 (en) | Pulse oximeter | |
| JP2004135854A (ja) | 反射型光電脈波検出装置および反射型オキシメータ | |
| US20080027294A1 (en) | Physiological parameter tracking system | |
| JPWO2007139192A1 (ja) | 光学的測定装置、光学的測定方法、及び光学的測定プログラムを格納した記憶媒体 | |
| US8923944B2 (en) | Cancellation of light shunting | |
| JP2004290544A (ja) | 血液分析装置 | |
| JP2004290412A (ja) | 血液分析装置 | |
| WO2009124077A1 (en) | Detection of site oximetry degradation | |
| JP2019058210A (ja) | 酸素飽和度推定装置、酸素飽和度測定装置、酸素飽和度推定プログラム、及び酸素飽和度推定方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20051228 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20051228 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070221 |
|
| 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: 20090224 |
|
| 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: 20090313 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120403 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 4284674 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130403 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140403 Year of fee payment: 5 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |