JPH0749304A - 散乱吸収体の内部情報計測方法及び装置 - Google Patents
散乱吸収体の内部情報計測方法及び装置Info
- Publication number
- JPH0749304A JPH0749304A JP6021400A JP2140094A JPH0749304A JP H0749304 A JPH0749304 A JP H0749304A JP 6021400 A JP6021400 A JP 6021400A JP 2140094 A JP2140094 A JP 2140094A JP H0749304 A JPH0749304 A JP H0749304A
- Authority
- JP
- Japan
- Prior art keywords
- light
- scattering medium
- scattering
- internal information
- predetermined wavelength
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- 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/47—Scattering, i.e. diffuse reflection
- G01N21/49—Scattering, i.e. diffuse reflection within a body or fluid
Landscapes
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
報を計測する。 【構成】 散乱や吸収の影響を受けながら散乱吸収体の
内部を拡散伝搬した光を散乱吸収体の外部で測定し、そ
の測定値を演算処理して散乱吸収体の内部情報を計測す
る。このとき、散乱吸収体内部を拡散伝搬する光の振る
舞いおよびその結果である信号つまり光検出信号が、散
乱吸収体内部の散乱成分や吸収成分などの性質や濃度に
依存することを利用して、3種以上の異なる光入射位置
−光検出位置間距離で測定した3種以上の検出信号(計
測値)を演算処理する。ここで計測される内部情報は、
散乱吸収体の吸収係数や輸送散乱係数などの絶対値であ
り、さらにこれらを演算処理して、特定散乱成分や特定
吸収成分の濃度に関する情報などを計測することができ
る。
Description
るいは連続光を生体などの散乱吸収体に入射し、散乱吸
収体の内部を拡散伝搬した光を検出して、散乱吸収体の
内部情報を非侵襲に計測するもので、散乱吸収体の吸収
係数や輸送散乱係数の絶対値、さらには特定吸収成分の
濃度の定量、およびそれらの時間変化や空間分布などを
計測すること、あるいはこれらの計測の精度を高めるこ
とができる散乱吸収体の内部情報計測方法及び装置に関
するものである。
のような散乱吸収体内部の吸収係数や輸送散乱係数の精
密な測定、あるいは特定吸収成分の濃度の精密な測定に
対する要請は非常に強く、従来いくつかの報告や試みが
ある。その中の主なものを、この節の最後に参考文献
1-5)として一括して示す。なお、以降では、x)によって
参考文献の番号を示すものとする。
に散乱、吸収されるため、直進しない。吸収が零である
散乱体では、光の総量が減ることはないが、光は散乱成
分によってランダムに散乱されるため、折れ曲がりなが
ら進む。この場合、光が散乱の影響を受けずに進める距
離は、平均自由行程あるいは平均拡散長と呼ばれ、輸送
散乱係数μS ′の逆数となる。生体試料などの場合、こ
の平均自由行程は2mm程度である1)。また、散乱吸収体
では散乱成分のほかに吸収成分も含まれているので、光
が散乱して進む距離に応じて吸収が生じ、その距離に対
して光量が指数関数的に減衰する。
の精密測定の分野の従来技術としては、連続光やパルス
光入射に対する透過光あるいは反射光の光量(時間積分
値)を計測するものと、パルス光入射に対する透過光あ
るいは反射光を時間分解計測し、その波形を解析するも
のがある。これらのうち、前者はよく知られているラン
バート・ベール(Lambert-Beer)則を基本原理とする吸光
度測定を利用するもので、これと併せてさらに2波長分
光測定の原理を応用するものもある。ここで、ランバー
ト・ベール則は、試料の吸光度(absorbance,またはopti
cal density)が試料のモル吸光係数、モル濃度、および
試料の厚さの積に比例する、あるいは厚さを一定とした
とき吸光度差が試料の濃度差に比例するというものであ
る。
場合、散乱吸収体の光入射位置と光検出位置との間を拡
散伝搬する光の平均光路長(optical pathlength)が、
散乱吸収体の吸収係数μa に依存して変化する。したが
って、従来の光路長を一定と仮定する散乱吸収体に対す
る吸光度測定では、この平均光路長の吸収係数依存性が
大きな問題となり、吸収係数あるいは吸収成分の濃度な
どを正確に計測することができない、あるいはそれらし
い計測ができても測定誤差が大きくなるため実用できな
い。例えば、特願昭62−248590、特願昭62−
336197、特願平2−231378などは、光路長
が一定であると仮定して吸光度を測定することを基本原
理としているから、上記光路長変化に起因する測定誤差
を避けることができない。
方法で測定した平均光路長を利用する方法もあるが、上
記のように平均光路長が吸収係数に依存して変化するの
で、これを一定値として近似することに起因する誤差を
避けることはできない。その他、パルス光などを使用し
て吸光度差を測定する方法、これらに2波長分光の原理
を応用する方法、3種類以上の波長の光を使用して吸光
度を測定する方法2)もあるが、いずれの場合にも散乱吸
収体に対して光路長が一定であると仮定して吸光度測定
法を適用しているため、これら全てに共通して光路長変
化に起因する大きな誤差が発生し、十分な計測精度が得
られない。
ス光や変調光を使用して、透過光あるいは反射光を時間
分解計測して波形解析する方法は、時間分解計測をする
から、光量つまり光信号の時間積分値を測定する本願の
発明の計測方法や装置に比較して、計測方法や装置がき
わめて複雑になり、かつ装置が高価になるという難点が
ある。実際に、パルス光を散乱吸収体に入射したときの
反射光あるいは透過光を時間分解計測し、その波形から
内部の吸収情報を計測する試みがいくつかある3-5)。こ
のとき測定される光信号は、散乱と吸収の影響によっ
て、だらだらと減衰する長い裾をひく信号となる。Patt
erson らは、均一な散乱吸収体のモデルを考えて、パル
ス光入射に対する出力信号を解析的に求めた3)。Patter
son らが求めた式が与える光信号の強度の時間波形は、
均一な散乱吸収体を用いた実験で得られる波形とほぼ一
致する。彼らおよび本願の発明者らの実験結果によれ
ば、散乱吸収体を構成する吸収成分の吸収係数は、上記
の光信号が十分に減衰したとき、つまり十分時間が経過
したときの波形の傾斜(微分値)で与えられる。しか
し、吸収係数を求める部分の光信号が十分に減衰してい
る必要があるということは、信号が微弱であるというこ
とと等価であるから、被計測信号の信号対雑音比(SN
比)が悪くなり、必然的に測定誤差が増大するので、実
用化することが困難である。以上の他にも種々の試みが
あるが、いずれも、計測精度が不十分である。
光測定の原理を応用するものは、次の共通する問題があ
る。そもそも散乱吸収体の吸光度測定では、吸光係数(e
xtinction coefficient)が定義上、輸送散乱係数と吸収
係数の和となり、これらが同格のパラメータとして取り
扱われるから、そのままでは散乱と吸収の影響を分離し
て吸収の影響、例えば吸収係数を測定することが不可能
である。そこで通常は、2波長分光測定の原理を利用し
てこれを解決する。この2波長分光測定では、吸収成分
に対して吸収係数が異なる適宜の2種以上の光を用いて
吸光度を測定し、この際に前記2種以上の光に対する散
乱係数あるいは輸送散乱係数が同一、あるいは差があっ
ても極めて小さいものと仮定して、前記2種類以上の光
に対する吸光度差から散乱の影響を消去して、吸収係数
あるいは吸収成分の濃度を求めるものである。この方法
は、上記の測定原理から明らかなように、異なる波長に
対して散乱係数あるいは輸送散乱係数が等しいと仮定す
ることに起因する誤差が発生するという本質的な欠点が
ある。
告されている上記分野の計測法や計測装置では、いずれ
の場合にも十分な測定精度が得られていない。このよう
な状況に鑑み、本願の発明者は、上記分野の測定の基本
となる散乱吸収体内部の光の振る舞いについて、理論的
かつ実験的に詳細に解析、検討して、従来の解析方法は
散乱吸収体内部の複雑な光の振る舞いを記述するのに十
分でない、あるいは従来の計測法や計測装置の計測原理
を導く過程で導入している多くの近似が、実際の測定対
象や測定現場で成立しないため、その結果として大きな
誤差が測定値に含まれることなどを明らかにした。本願
の発明は、これらの解析、検討、実験などに基づいて、
上記の誤差を著しく小さくして測定精度を高める新しい
計測方法と計測装置を提供することを目的とするもので
ある。
ance andTransmittance of Tissues: Principle and Ap
plication, IEEE J. QuatumElectron., Vol.26, No.12,
pp.2186-2199 (1990) 2) I. Oda, Y. Ito, H. Eda, T. Tamura, M. Takada,
R. Abumi, K. Nagai,H. Nakagawa, and M. Tamura: Non
-invasive hemoglobin oxygenationmonitor and comput
ed tomography by NIR spectrophotometry, Proc.SPIE,
Vol.1431, pp.284-293 (1991) 3) M.S. Patterson, B. Chance, and B.C. Wilson: Tim
e resolvedreflectance and transmittance for the no
n-invasive measurement oftissue optical propertie
s, Applied Optics, Vol.28, No.12,pp.2331-2336 (198
9) 4) M.S. Patterson, J.D. Moulton, B.C. Wilson, and
B. Chance:Application of time-resolved light scatt
ering measurements tophotodynamic theraphy dosimet
ry, Proc. SPIE, Vol.1203, pp.62-75(1990) 5) E.M. Sevick, N.G. Wang, and B.Chance: Time-depe
ndent photon imaging,Proc. SPIE, Vol.1599, pp.273-
283 (1991)
の内部情報計測方法及び装置は、光路長変化を十分に考
慮した新しい知見、つまり所定波長の光が散乱吸収体の
内部を拡散伝搬して所定の検出位置に到達するときの光
検出信号(計測値)と複数の既知のパラメータとの間の
複数の連立関係を利用して、散乱吸収体の吸収係数や輸
送散乱係数の絶対値などの内部情報を計測する。さら
に、複数の波長の光を使用して上記計測を行って、特定
吸収成分の濃度などの吸収に関する情報を計測する。
情報計測方法は、(a)所定波長の光を散乱吸収体に入
射し、(b)3種類以上の異なる光入射位置−光検出位
置間距離に対応する検出位置で、散乱吸収体の内部を拡
散伝搬した所定波長の光を検出して3種以上の信号を取
得し、(c)所定波長の光がそれぞれ拡散伝搬して検出
位置に到達するときの散乱吸収体内の拡散伝搬経路にお
ける散乱特性および吸収特性と前記検出信号との3つ以
上の関係に基づいて、検出信号を演算処理して散乱吸収
体の内部情報を導出する、ことを特徴とする。
特定吸収成分に対して吸収係数が異なる2種類以上の光
であり、演算処理は1種類の所定波長による計測によっ
て得られる散乱吸収体の第1次内部情報の導出処理と、
2種類以上の波長の光に対する第1次内部情報を利用し
た散乱吸収体の前記特定成分に関する情報である第2次
内部情報の導出処理とである、ことを特徴としてもよ
い。
測装置は、(a)所定波長の光を発生する光発生手段
と、(b)所定波長の光を散乱吸収体に入射させる光入
射手段と、(c)3種類以上の異なる光入射位置−光検
出位置間距離に対応する検出位置で、散乱吸収体の内部
を拡散伝搬した所定波長の光を検出して、3種以上の信
号を取得する光検出手段と、(d)所定波長の光がそれ
ぞれ拡散伝搬して検出位置に到達するときの散乱吸収体
内の拡散伝搬経路における散乱特性および吸収特性と前
記検出信号との3つ以上の関係に基づいて、検出信号を
演算処理して散乱吸収体の内部情報を導出する演算処理
手段と、を備えることを特徴とする。
特定吸収成分に対して吸収係数が異なる2種類以上の光
であり、演算処理は1種類の所定波長による計測によっ
て得られる散乱吸収体の第1次内部情報の導出処理と、
2種類以上の波長の光に対する第1次内部情報を利用し
た散乱吸収体の前記特定成分に関する情報である第2次
内部情報の導出処理とを実行する、ことを特徴としても
よい。
び装置は、散乱と吸収の影響を受けながら散乱吸収体の
内部を拡散伝搬した光を散乱吸収体の外部で測定し、そ
の測定値を演算処理して散乱吸収体の内部情報を計測す
る。このとき、散乱吸収体内部を拡散伝搬する光の振る
舞いおよびその結果である信号つまり光検出信号が、散
乱吸収体内部の散乱成分や吸収成分などの性質や濃度に
依存することを利用して、3種以上の異なる光入射位置
−光検出位置間距離で測定した3種以上の検出信号(計
測値)を演算処理する。ここで計測される内部情報は、
散乱吸収体の吸収係数や輸送散乱係数などの絶対値であ
り、さらにこれらを演算処理して、特定散乱成分や特定
吸収成分の濃度に関する情報などを計測することができ
る。また、上記の計測を計測対象の多箇所で行えば内部
情報の空間分布、また、同一場所で異なる時間に計測す
れば内部情報の時間変化を計測することができる。
法および装置では、所定波長の光を散乱吸収体に入射
し、3種以上の異なる光入射位置−光検出位置間距離に
対応する散乱吸収体の内部を拡散伝搬した光を検出し
て、3種類以上の検出信号(計測値)を得る。つぎに、
前記所定波長の光が前記検出位置に拡散伝搬して到達す
るときの散乱および吸収と、前記検出信号(計測値)と
の3つ以上の連立関係に基づいて、前記3種類以上の検
出信号(計測値)と既知のパラメータとを演算処理し
て、前述したような散乱吸収体の内部情報を定量計測す
る。以上で用いられる検出信号(計測値)と各種パラメ
ータおよび散乱吸収体の内部情報との関係は、本願によ
って初めて開示されるものである。
して上記の計測を行い、その際に、散乱吸収体内部の特
定吸収成分に対する吸収係数が異なるように前記所定波
長を選ぶと、散乱吸収体内部の特定吸収成分の精密な定
量計測ができる。また、本願の発明による計測では複数
の波長に対する輸送散乱係数と吸収係数を高精度に計測
することができるから、従来の散乱吸収体に対する2波
長分光計測では不可能であった、2つの波長に対して散
乱係数が異なる試料あるいは波長の場合でも、散乱吸収
体の吸収係数を高精度で計測することができる。このよ
うな新しい2波長分光計測は、例えば特定物質の濃度、
ヘモグロビンの酸素飽和度、およびこれらの空間分布や
時間変化などの精密計測に利用される。
共通する光路長変化にともなう測定誤差の問題、および
2波長分光測定における2波長間の散乱係数の差の問題
が解決され、計測精度が著しく向上する。さらに、本願
の発明は時間分解計測しない方法であり、光量つまり光
検出信号の時間積分値を計測する方法であるから、装置
構成が極めて簡単になるとともに、光の利用率が向上し
て測定時間が短くなり、信号対雑音比(SN比)も向上
する。これらから、本願の発明によってさらに、手術中
や運動中の人の脳内の酸素濃度などのモニタリングある
いはリアルタイム計測などが可能になる。
詳細に説明する。なお、図面の説明において同一の要素
には同一の符号を付し、重複する説明を省略する。
がら拡散伝搬し、その一部は表面に出てくる。つまり、
光は散乱吸収体の内部で散乱、吸収されるが、これを透
過することができる。例えば図1は、散乱吸収体内部の
光の振る舞いを示した図であり、散乱吸収体の1点に光
を入射し、他の点で拡散光(反射光)を検出する様子を
示す。散乱吸収体に入射された光は、内部でランダムに
散乱されてほぼ全域に広がるが、この図では、光検出器
で検出される光子の飛跡のみが模式的に示してある。つ
まり図1は、実際の計測に利用される光子の飛跡を示し
ている。
は、光拡散理論(Photon DiffusionTheory)によって、
かなり厳密に記述、解析できることが知られている。こ
の光拡散理論によれば、散乱吸収体に入射された光パル
スは、散乱吸収体内部で散乱、吸収されながら拡散伝搬
する際に、パルス幅が広がっていくことなどが示され
る。図2は散乱吸収体内部における光検出量を求めるた
めの仮想光源の設定を示す図である。この図に示すよう
に、散乱吸収体の表面の1点(ρ=0,z=0)に光を
入射したとき、点P(ρ,z)で検出される光つまり光
信号は、光拡散方程式から導出することができる。この
場合、散乱吸収体の表面と外部では光拡散が存在しない
から、これが満足されるように境界条件を設定する必要
がある。Patterson らは、負極性の点光源を仮定してこ
の境界条件を満たす方法を開発し、これから得られる理
論値と実際の実験値とがよく一致することを報告してい
る。図2では、この負極性の点光源を○印点(ρ=0,
z=−z0 )に導入してある。
子1個1個の振る舞いは、コンピュータを利用したモン
テカルロシミュレーション(Monte Carlo Calculation)
によって解析、実験、検討することができる。また、実
際に散乱吸収体の物理模型や生体試料などを用いて実験
することもできる。最近得られている知見では、光拡散
理論から導かれる結果、モンテカルロシミュレーション
の結果、およびサンプルによる実験結果の間には、それ
ぞれ良い一致が得られている。したがって、散乱吸収体
の内部の光の振る舞いは光拡散方程式で十分に記述され
ると考えてよい。また、本願の発明者による種々の解
析、実験、検討、および最近のさらに詳しい散乱吸収体
内部の光の振る舞いの解析、実験などでも、散乱吸収体
の内部の光の振る舞いが光拡散方程式で正確に記述され
ることが示されている。
の振る舞いを記述する光拡散方程式は、例えば、光子流
動率(fluence rate) φを用いて次のように記述され
る。
ence rate)〔光子/mm2 ・sec 〕(但し、r はベクトル
である。) D:光拡散係数(diffusion coefficient)〔mm〕 μa :吸収係数 (absorption coefficient) 〔mm-1〕 c:散乱吸収体内部の光速度〔mm/sec〕(屈折率によっ
て定まる) S(r,t):光源〔光子/mm3 ・sec 〕 である。
ら、原点(r=0)にt=0のとき入射する光インパル
スは、次のように表される。
次のようになる。
とき入射する光インパルスである。
は、 μS :散乱係数 (scattering coefficient) 〔mm-1〕 μS ′:輸送散乱係数 (transport scattering coeffic
ient) 〔mm-1〕 μtr:輸送減衰係数 (transport attenuation coeffici
ent)〔mm-1〕 μeff :有効減衰係数 (effective attenuation coeffi
cient)〔mm-1〕 g:散乱角θに対するcosθの平均値 (mean cosine
of the scatteringangle) としたとき、 D=〔3 (μa +μS ′)〕-1=(3μtr)-1 …(1.4a) μS ′=(1−g)μS …(1.4b) μtr=μa +μS ′=μa +(1−g)μS …(1.4c) μeff =〔3μa (μa +μS ′)〕1/2 ={3μa 〔μa +(1−g)μs 〕}1/2 …(1.4d) なる関係にある。
光パルスを入射する場合の境界条件は、前述したPatter
son らが示したように、図2に示すような負極性の点光
源を位置(ρ=0,z=−z0 )に置くことで実現され
る。通常z0 は1/μS ′程度の値であるが、厳密には
光の入射方法や散乱吸収体の散乱成分の性質などによつ
て変化する。このことは、本願の発明者らのモンテカル
ロシミュレーションでも確認されている。
を解いて、散乱吸収体の表面の任意の位置(ρ,0)に
おける光信号R (ρ,0,t) [光子/mm2 ・sec ]を
求めると、 R (ρ,0,t) =(4πD)-3/2z0 t-5/2 ×exp〔−(z0 2 +ρ2 )/(4cDt)〕 ×exp(−cμa t) …(1.5) となる。
t)を時間積分したものであるから、
のように解くことができる。
場合には、計測値I[ρ]および既知値ρを(1.7)
または(1.8)式に代入して、散乱吸収体のμeff を
求めることができる。図3は、このような計測値I
[ρ]と検出距離ρとの関係を示す。ここでは、μS ′
=0.45,z0 =1/0.45とし、μa をパラメー
タとして示してある。
(1.8)式のz0 は、実際にはz0 ≠μS ′となる。
このことは、次の例を考えれば容易に理解できる。つま
り、通常の計測系では空気中から散乱吸収体に光を入射
し、散乱吸収体の内部を拡散伝搬した光を空気中で検出
する。この場合、前記の境界条件、つまり、図2に示し
たz0 は空気と散乱吸収体との屈折率の差に依存する値
となる。本願の発明者らの詳細な解析およびモンテカル
ロシミュレーションによれば、z0 =a/(μa +
μS ′)とおいたとき、通常の計測状態では、aの値が
概ね、0.3<a<3の範囲にあることが確認されてい
る。
あるいは散乱吸収体から光が出射するとき、何らかの光
減衰、例えば表面の着色層などでの吸収があると、前記
の(1.5)式で表される光信号R、および(1.7)
式や(1.8)式で表される検出光量が、上記の減衰に
比例して小さくなる。また通常は、出射光は光検出器で
検出される。したがって、検出信号(計測値)は光検出
器の量子効率に関係した値になる。しかし、この量子効
率は、上記の光減衰率に含めて考えることができる。
高い計測精度が得られる計測方法について説明する。上
記のような光減衰の減衰率をηとすると、計測値Qは
(1.8)式から、 lnQ[ρ]=ln〔ηz0 /(2π)〕−(3/2)ln(z0 2 +ρ2 ) +ln〔μeff (z0 2 +ρ2 )1/2 +1〕 −μeff (z0 2 +ρ2 )1/2 …(1.9) ただし、 z0 =a/(μa +μS ′) …(1.10) となる。ここで、aは経験的あるいは実験的に定められ
る定数であり、通常は0.3<a<3である。
9)式の未知数は、μa 、μS ′、及びηの3個であ
る。したがって、制御あるいは設定可能な既知数である
光検出距離ρを3つ以上の異なる値に設定して、3種以
上の計測値Qを取得し、3種以上のQに対する3つ以上
の連立関係に基づいて、3個の未知数つまり光減衰率
η、吸収係数μa 及び輸送散乱係数μS ′を算出するこ
とが可能になる。本願の発明では、上記のような基本的
な知見をさらに展開して、次のように応用する。
報計測方法を説明する図である。つまり、均一な散乱吸
収体に対して、所定の波長の光を用いて、図4に示した
ような異なる3種の検出距離で光検出して、前記計測値
Q[ρ]を得る。なお、以下では、光入射位置と光検出
位置との距離を単に検出距離と呼ぶ。いま、検出距離ρ
1 、ρ2 およびρ3 に対する3種類の計測値をそれぞ
れ、 Q[ρ1 ]:検出距離ρ1 での計測値 Q[ρ2 ]:検出距離ρ2 での計測値 Q[ρ3 ]:検出距離ρ3 での計測値 のように表すことにする。すると、3種の検出距離に対
する3つの計測値は、次の3連の連立方程式を満足す
る。
後、総称して(1.11)とも呼ぶ)を構成する3連の
式は、それぞれ独立であり、かつ未知数はμa、μS ′
およびηの3個である。したがって、3個の計測値と既
知あるいは別の方法で測定した値ρ1 、ρ2 、およびρ
3 を用いて、3個の未知数μa 、μS ′およびη、さら
には必要に応じてμeff やz0 を求めることができる。
なお、普通はηやz0 を求める必要はない。また、以上
のようなμa とμS ′、あるいはμeff を求めるには、
(1.6)または(1.7)式から導かれる3連のそれ
ぞれ独立した式であればどんな形の式でもよい。さら
に、ηが既知数として取り扱える場合には、(1.1
1)式の中の任意の2つの式からなる連立方程式、ある
いは(1.6)または(1.7)式から導かれる2連の
それぞれ独立な式からなる連立方程式を解くことによっ
て、μa とμS ′、さらにはμeff を求めることができ
る。以上のような連立方程式を解く計算は、コンピータ
を利用して高速に実行することができる。また、当然で
あるが、互いに独立な計測値Qの数が多いほど上記の計
算精度が向上する。
2 の光を用いると、それぞれの波長に対する吸収係数μ
a1とμa2、およびそれぞれの波長に対する輸送散乱係数
μs1′とμs2′を求めることができる。したがって、例
えばμa1/μa2の値からヘモグロビンの酸素飽和度など
を算出することができる。このような2波長分光計測で
は、散乱吸収体に含まれる吸収成分の吸収係数の波長依
存性を利用する。例えば、図5に吸収スペクトルを示す
酸化・還元ヘモグロビンや酸化・還元ミオグロビンの計
測では、酸化・還元に対して吸収係数の差が大きい波長
の光、つまり700nmから1.2μmの波長の光がよ
く使用される。
測方法や装置では、前述したように輸送散乱係数μs1′
とμs2′の絶対値を求めることができることである。し
たがって、従来の2波長分光における仮定、つまり異な
る2つの波長の光に対する散乱成分の散乱係数が同一ま
たは差があっても極めて小さいという仮定が成立しなく
ても、本願の発明による計測方法および装置によって精
度の高い計測ができる。このように、本願の発明による
計測方法や装置は、従来の2波長分光計測の適用範囲を
大幅に拡大して、計測精度を大きく向上することができ
る。
するので、入射する光パルスの波形はどんな形のもので
あっても、積分時間が規定できるものであれば、問題な
く利用することができる。例えば方形波や連続光でもよ
いことは明らかである。また、3種以上の波長の光や4
種以上の検出距離で計測すれば、さらに計測精度を向上
させること、あるいはバックグラウンド吸収のある試料
を計測することなどが可能になる。なお、図4では1つ
の光入射位置に対して異なる3つの位置で光検出した
が、異なる3つの位置に光を入射して1つの光検出位置
で検出する方法などでもよく、要は、3種以上の検出距
離を利用すればよい。
散乱吸収体である場合について説明したが、実際には、
有限の大きさの散乱吸収体を計測することが多い。この
場合、図4の散乱吸収体の表面および外で前記の境界条
件が満足され、かつ散乱吸収体内部を拡散伝搬する大部
分の光に対して光拡散の条件が成立すればよい。この条
件は、検出距離ρ1 、ρ2 およびρ3 に比較して散乱吸
収体が十分大きいと見做せるときに成立する。例えば、
図4に示した波線内の部分は、検出位置に到達する大部
分の光が通過する領域を示す。なお、このような領域が
あることは、例えば、前述したSevickらが報告した散乱
吸収体内部のバナナ状の光の束、あるいは本願の発明者
らのモンテカルロシミュレーションの結果からも明らか
である。
較して散乱吸収体の厚さ(z方向)が十分に大きいと見
做せないときは、散乱吸収体の背面側(光入射面と対向
する反対側の表面)に、境界条件を満足させるような第
2の負極性と正極性の点光源を導入すればよい。このと
き、原理的には、この第2の負極性と正極性の点光源の
影響を補正するために、さらに第3の負極性と正極性の
点光源、さらにそれを補正するために第4の負極性と正
極性の点光源などが必要となる。しかし、これらで導入
する点光源の出力光への影響は、その順序に従って急激
に減少するので、有限の点光源で近似することができ
る。この場合、前出の(1.5)式乃至(1.9)式
は、追加した負極性と正極性の点光源の影響に関する項
を追加して修正する必要がある。しかし、この修正は、
項の数が増えるだけで、未知数の数や各式の独立性など
は不変であるから、前記と同様のことが成立する。
スラブ状の試料に適用することもでき、その様子を図6
に示す。ただし、図6では3つの位置から光を入射して
いるが、1箇所から光を入射して検出距離の異なる3箇
所で検出してもよい。なお、スラブ状の試料の場合に
は、光拡散方程式(1.3)から導かれる結果、つまり
前述の(1.7)式に相当する式は、検出距離をdとし
たとき、例えば、
(1.9)式の式の形が変わっただけ、つまり未知数や
既知数の数が同じで、パラメータも同一であるから、前
記の場合と同様に演算処理することができ、前記と同様
の結果が得られる。したがって、前記と同様にして3個
の計測値と既知の値ρ1 、ρ2 およびρ3 を用いて、3
個の未知数μa 、μS ′およびη、さらには必要に応じ
てμeff を求めることができる。なお、通常はηやz0
を求める必要はない。また、ηが既知数である場合に
は、任意の2つの独立な式の連立方程式からμa 、
μS ′およびμeff を求めることができる。さらに、異
なる2つの波長の光を利用して、ヘモグロビンの酸素飽
和度など、吸収成分に関する内部情報を求めることがで
きる。また、前記と同様にして、連立方程式を解く計算
はコンピュータを利用して高速に実行される。
の内部情報は、光入射位置から光検出位置に拡散伝搬す
る光の光路に沿う情報の平均値である。したがって、内
部情報が図4示した光入射位置−光検出位置間の距離、
あるいは図6に示した光入射位置の間隔に比較して、十
分粗に分散している場合には、各種の計測値の空間分布
の画像化、つまり簡単なイメージングができる。さら
に、当然ではあるが前記の計測を異なる時間に行えば、
前記内部情報の時間変化を計測することができ、脳内酸
素量などのモニタリングや計測に応用することができ
る。以上のような場合の画像化演算などは、メモリ、デ
イスプレイなどを備えるコンピュータ装置によって高速
に実行される。
つの例として、本願の発明による方法と従来の一定の光
路長を仮定する吸光度法との計測誤差について述べてお
く。前述の図4に示したような半空間を占める散乱吸収
体を考えた場合、検出される光信号は(1.5)式で表
せる。このような信号に対して、平均光路長<L>、つ
まり時間分解波形の重心を考えると、
て、μS ′=0.45〔mm-1〕、z0 =1/μS ′、ま
た計測条件としてρ=50〔mm〕を考え、μa が0.0
1から0.02に変化するときの平均光路長をそれぞ
れ、<L(μa =0.01)>と<L(μa =0.0
2)>とすれば、平均光路長を一定と仮定するときの平
均光路長に対する誤差Δは、
と仮定して求めた吸光度は、14%程度の誤差を含むこ
とがわかる。これに対して、本願の発明の方法および装
置では、上記平均光路長の変化を考慮した関係を利用し
て演算処理しているので、上記の誤差をなくすことがで
きる。以上のことは、スラブ状の散乱吸収体を計測する
場合でも同様であり、誤差の程度の差は少しあるが、本
願の発明の有効性が明らかである。
する情報を演算処理して導出する方法について説明す
る。
法 哺乳類の脳における吸収成分の主なものは、水、チトク
ローム(cytochrom)、酸化および還元ヘモグロビンであ
る。近赤外線領域での水とチトクロームの吸収は、酸化
および還元ヘモグロビンに対して、ほぼ無視することが
できる程度に少ない。また、酸化および還元ヘモグロビ
ンは、前出の図5に示すように吸収スペクトルが異な
る。さらに、頭蓋骨は、近赤外線に対して散乱体と考え
てよい。
種の光を考える。すると、λ1 とλ2 に対する吸収係数
μa1とμa2は、ランバート・ベール則によって、次のよ
うになる。
収係数〔mm-1・M-1〕 εHbO,1 :酸化ヘモグロビンの波長λ1 に対するモル吸
収係数〔mm-1・M-1〕 εHb,2: 還元ヘモグロビンの波長λ2 に対するモル吸
収係数〔mm-1・M-1〕 εHbO,2 :酸化ヘモグロビンの波長λ2 に対するモル吸
収係数〔mm-1・M-1〕 〔Hb〕:還元ヘモグロビンのモル濃度[M] 〔HbO〕:酸化ヘモグロビンのモル濃度[M] である。
るいはμa1/μa2と既知のパラメータεHb,1、
εHbO,1 、εHb,2およびεHbO,2 とを用いて、還元ヘモ
グロビンのモル濃度〔Hb〕、酸化ヘモグロビンのモル
濃度〔HbO〕、ヘモグロビンのモル濃度〔Hb〕+
〔HbO〕および酸素飽和度Yを算出することができ
る。なお、酸化および還元ヘモグロビンに対して吸収が
同一になる波長(≒800 nm 、isosbestic wavelengt
h )を使用すれば上記の式はさらに簡単になる。
場合がある。この場合は、波長λ1 、λ2 およびλ3 に
対するバックグラウンド吸収を、それぞれa1、a2 およ
びa3 とすれば、ランバート・ベール則によって次式が
成立する。
宜に波長を選択すれば、前記と全く同様にして、3種類
の波長の光に対する計測値と既知のパラメータの値を用
いてμa1、μa2およびμa3を求めて、これらと既知のパ
ラメータの値を(2.4)式あるいは(2.5)式に代
入して、〔Hb〕、〔HbO〕、〔Hb〕+〔Hb
O〕、あるいは飽和度Yを求めることができる。なお、
上記の条件a1 ≒a2 ≒a3 は、計測対象である生体な
どに対して波長をうまく選べば実現できる。
法を実施する装置の第1実施例である。レーザーダイオ
ードなどを使用した光源2は所定波長λ1 とλ2 の光を
発生する。計測に使用する光の波長は、計測対象に応じ
て適宜に選ぶ必要がある。一般に生体では、ヘモグロビ
ンなどの吸収の関係から700〔nm〕以上の光がよい。
また、酸化ヘモグロビンと還元ヘモグロビンでは、前出
の図5に示したように、吸収係数が異なるから、適宜に
波長を選択して使用することによって、これらを分離し
て計測することができる。光源には、発光ダイオード、
レーザーダイオード、HeNeレーザーなど種々のもの
が使用できる。また、前述したようにパルス光や方形波
光を発生するものでもよい。
され、光ガイド6を通して計測対象である散乱吸収体2
0の表面に入射される。この場合、後述するような集光
レンズやピンホールを使用する方法もある。つまり、前
述したように生体試料などの散乱吸収体では平均拡散長
が2〔mm〕程度であるため、入射光は約2〔mm〕直進す
るまでに散乱し、光の方向性がなくなる。したがって、
数cm以上の散乱吸収体では平均拡散長の影響を無視する
ことができるから、スポット状に光を入射すれば良い。
また、太いビーム状の光を散乱吸収体に入射してもよ
い。この場合には、複数のスポット状光源が並んでいる
と考えればよい。
図7の実施例では微小になっている。しかし実際には、
これを大きくして、この空間に計測対象である散乱吸収
体20とほぼ等しい屈折率と散乱係数をもつ液状体やゼ
リー状物体(以下、インターフェース材と呼ぶ)を満た
しておいても良い。つまり、光はこのインターフェース
材中を拡散伝搬して計測対象に入射するから何ら問題は
生じない。また、散乱吸収体20の表面反射が問題にな
るときには、インターフェース材を適宜に選択すること
によって、表面反射などの影響を低減することができ
る。
は、前記光入射位置から距離ρ1 、ρ2 およびρ3 の位
置に置いた光ガイド8、9および10で受光される。こ
こでも前記と同様の理由によって、インターフェース材
を使用してもよい。
および第3の光検出器14は、前記受光信号を電気信号
に変換して増幅し、それぞれ検出信号a、bおよびcを
出力する。この光検出器12、13および14は、光電
子増倍管のほか、光電管、フォトダイオード、アバラン
シェフォトダイオード、PINフォトダイオードなど、
あらゆる種類の光検出器を使用することができる。光検
出器の選択に際しては、所定波長の光が検出できる分光
感度特性をもっていれば良い。また、光信号が微弱であ
るときは高感度あるいは高利得の光検出器を使用する。
光検出器の受光面以外の場所は、光を吸収あるいは遮光
する構造にしておくことが望ましい。また、散乱吸収体
の内部を拡散伝搬した光が複数の波長の光を含む場合に
は、光検出器12、13および14と散乱吸収体20と
の間の適宜の位置に波長選択フィルタ(図示せず)を入
れる。
第2の光検出器13および第3の光検出器14からの検
出信号a、bおよびcを信号処理して、それぞれ検出光
量に比例した計測値d、eおよびfに変換する。具体的
には、演算処理部16は光源2の光発生に同期した信号
gを利用して、検出信号a、bおよびcに対する時間域
での積分演算を行って、検出光量に比例する計測値d、
eおよびfを演算する。ただし、パルス光などを利用す
る場合には、この同期信号gを省略することができる。
つぎに、波長λ1 とλ2 に対して、計測値d、eおよび
fと、予め別の方法で測定あるいは設定された既知のパ
ラメータρ1 、ρ2 およびρ3 の値とを、前述した3連
の関係に基づいて演算処理し、散乱吸収体の内部情報つ
まり波長λ1 と波長λ2 とに対する吸収係数と輸送散乱
係数、つまりμa1とμs1′およびμa2とμs2′を算出す
る。この種の演算処理は、演算処理手段に組み込んだマ
イクロコンピュータなどで高速に実行することができ
る。また必要に応じて、上記で求めた内部情報つまりμ
a1とμa2の値を用いて、特定吸収成分の濃度などを算出
する。さらに、バックグラウンド吸収などがある場合に
は、既に述べたように3種以上の波長の光を使用する。
光させて使用する方法と、異なる波長の光を同時に含む
光を使用する方法とがある。具体的な波長選択手段とし
ては、ミラーを用いた光ビーム切り換え器、フィルター
を用いた波長切り換え器、光スイッチを用いた光切り換
え器などを利用する方法(図8(a))がある。また、
異なる波長の光を同軸状のビームにして、光入射点の直
前に設けた波長選択フィルタで波長選択する方法、ある
いはそのまま平行に散乱吸収体に入射して(図8b)、
前記光検出器の直前に設けた波長選択フィルタで波長選
択する方法(図9)、さらには6個の光検出器を用いて
それぞれの検出位置で2種の波長の光を並列に検出する
方法、およびこれらを適宜に組み合わせる方法などがあ
る。
としては、図7に示した光ガイドを使用する方法の他
に、集光レンズ(図10(a))、光ファイバー(図1
0(b))、ピンホール(図10(c))を利用する方
法、胃カメラのように体内から入射する方法(図10
(d))などがある。また、散乱吸収体の内部を拡散伝
搬した光を受光して検出する手段としては、図7に示し
た光ガイドを用いる方法以外に、直接検出する方法(図
11(a))、光ファイバーを用いる方法(図11
(b))、レンズを用いる方法(図11(c))などが
ある。
増幅する必要がある場合には、狭帯域アンプ(図12
(a))、ロックインアンプ(図12(b))などを利
用することができる。ロックインアンプを使用する場合
には、参照信号として前記同期信号gを利用する。この
方法は、方形波光やパルス光を用いてダイナミックレン
ジの高い計測を行うときに有効である。また、散乱吸収
体20が生体試料などで、表面部に吸収層がある場合に
は、検出距離ρ1 、ρ2 およびρ3 の場所の予め別の方
法で測定した吸収の値を用いて、前記計測値d、eおよ
びfの値を補正することもできる。
と光検出位置とを同期して走査(図示せず)させて、散
乱吸収体の各部の内部情報を求めて、フレームメモリー
(図示せず)に蓄積し、これをテレビ方式で読み出せ
ば、内部情報の分布を示す画像が得られる。また、異な
る時刻に計測すれば、内部情報の時間的変化を計測する
ことができる。前記演算処理部16は、このようにして
得た内部情報を記憶する機能をもち、図7の中の表示記
録手段18はこれらの途中経過や結果を表示あるいは記
録するものである。この際、これらの演算処理は、メモ
リ、ディスプレイなどを備えるコンピュータ装置によっ
て高速に実行することができる。
と同様の計測ができる。この場合には、拡散伝搬した透
過光を検出するように、光入射位置と光検出位置とを、
図6に示したように、散乱吸収体の対向する表面に設置
するとともに、検出距離が異なるように設定する必要が
ある。
方法を実施する装置の第2実施例の構成図である。この
装置は、散乱吸収体の光学定数を計測する装置で、容器
24の中に入れた試料20は散乱吸収体であり、その吸
収係数や輸送散乱係数は未知である。この実施例に示す
装置は、1つの所定波長の光を使用するもので、2波長
を使うこと以外の動作は第1の実施例の装置とほぼ同じ
である。また、図13では、前出の図7で説明したもの
と同一の機能をもつものに図7と同一の番号や記号が付
けてある。
通して容器24に入れた液状の試料20に入射される。
散乱吸収体である試料の内部を拡散伝搬した光は、前記
光入射位置から距離ρ1 、ρ2 およびρ3 の位置に置い
た光ガイド8、9および10で受光され、検出部22の
中の配置した光検出器(図示せず)で検出信号a、bお
よびcに変換され、信号ケーブル17を介して、次段の
演算処理部16に供給される。演算処理部16は、前記
検出信号a、bおよびcを信号処理して、それぞれ検出
光量に比例した計測値d、eおよびfに変換する。つぎ
に、計測値d、eおよびfと、予め別の方法で測定ある
いは設定された既知のパラメータρ1 、ρ2 およびρ3
の値とを、前述した3連の関係に基づいて演算処理し、
散乱吸収体の吸収係数μa1と輸送散乱係数μS1′を算出
する。この種の演算処理は、演算処理部16に組み込ま
れたマイクロコンピュータで高速に実行する。
段、光検出信号の増幅手段などは、第1の実施例で説明
したものが利用できる。このような装置は、未知の散乱
吸収体の光学定数を測定する汎用性の高い装置となり、
理工学、生物学、分光分析などの広い分野で利用するこ
とができる。
方法を実施する装置の第3の実施例の構成の一部分を示
す。この装置は、人の脳内の酸化ヘモグロビンの濃度あ
るいはヘモグロビンの酸素飽和度を計測あるいはモニタ
リングする装置である。バンド型の容器26は、鉢巻き
のように頭部30に装着され、脳内のヘモグロビンの酸
素飽和度が計測される。この実施例に示す装置は、第1
の実施例と同様に、2つの所定波長λ1 とλ2 の光を使
用し、その動作は第1の実施例の装置とほぼ同じであ
る。また、図14では、前出の図7で説明したものと同
一の機能をもつものに図7と同一の番号や記号が付けて
ある。
は、波長選択器4で波長選択され、光ガイド6を通して
計測対象である頭部30に入射される。この際、2つの
波長は適宜に選ぶ必要がある。散乱吸収体である頭部内
を拡散伝搬した光は、前記光入射位置から距離ρ1 、ρ
2 およびρ3 の位置に置いた光ガイド8、9および10
で受光され、第1の光検出器12、第2の光検出器13
および第3の光検出器14で電気信号に変換、増幅され
る。このとき、電源、信号などは容器26に付けたコネ
クタ28を介して、外部機器(図示せず)と接続され
る。信号に対する演算処理(図示せず)は、第1の実施
例と同様であり、この実施例ではヘモグロビンの酸素飽
和度を出力、表示する。また、上記で、信号を電波や光
信号に変換して、伝送することもできる。
段、光検出信号の増幅手段などは、第1の実施例で説明
したものが利用できる。また、人体頭部などでは、表面
反射や光ガイドと頭部間の隙間などが問題になることが
ある。この場合には、前述したインターフェース材を利
用するとよい。図15は、図14に示した装置の光入
射、検出部の変形例を詳細に示す。この例では、図14
に示した光ガイドを省略して、頭部と波長選択器4との
間、および頭部と光検出器12との間に、計測対象とほ
ぼ同一の散乱係数と吸収係数をもつインターフェース材
11を使用している。このようにすると、集光率が良く
なるとともに、ρ1 、ρ2 およびρ3 の位置に置いたそ
れぞれの光検出と頭部との光結合状態の差を小さくする
ことができる。
に、例えばマラソン中の人の脚の筋肉中の酸化ヘモグロ
ビン濃度の計測あるいはモニタリングなどに使用するこ
とができる。
散乱吸収体の内部情報計測方法および装置によれば、散
乱吸収体の内部情報の高精度計測、つまり輸送散乱係数
μS ′および吸収係数μa の絶対値測定が可能となる。
これらに2波長分光の原理を応用すれば、特定成分の濃
度などの内部情報の精密な計測ができる。また、前記内
部情報の空間分布や時間的変化を計測することができ
る。さらに、本発明を利用した計測装置では、光信号の
時間積分値を使用するから、光の利用率が高くなり、信
号対雑音比が大きくなって計測精度が高くなる。したが
って、人体頭部や胴体部、立ち木などの植物体などに対
して、簡単な装置で、非侵襲で内部情報を計測し、これ
らのイメージングなどを行うことが可能になる。
る。
の説明図である。
グラフである。
ある。
めの構成図である。
る。
る。
る。
ド、11…インターフェス材、12,13,14…光検
出器、16…演算処理部、17…信号ケーブル、18…
表示記録手段、20…散乱吸収体、22…検出部、24
…容器、26…バンド状容器、28…コネクター、30
…人体頭部。
Claims (4)
- 【請求項1】 所定波長の光を散乱吸収体に入射し、 3種以上の異なる光入射位置−光検出位置間距離に対応
する検出位置で、散乱吸収体の内部を拡散伝搬した前記
所定波長の光を検出して3種以上の検出信号を取得し、 前記所定波長の光がそれぞれ拡散伝搬して前記検出位置
に到達するときの前記散乱吸収体内の拡散伝搬経路にお
ける散乱特性および吸収特性と前記検出信号との3つ以
上の連立関係に基づいて、前記検出信号を演算処理して
散乱吸収体の内部情報を導出する、 ことを特徴とする散乱吸収体の内部情報計測方法。 - 【請求項2】 前記所定波長の光は、前記散乱吸収体内
部の特定成分に対して吸収係数が異なる2種類以上の波
長の光であり、 前記演算処理は、1種類の前記所定波長による計測によ
って得られる前記散乱吸収体の第1次内部情報の導出処
理と、前記2種以上の波長の光に対する前記第1次内部
情報を利用した前記散乱吸収体の前記特定成分に関する
情報である第2次内部情報の導出処理とである、 ことを特徴とする請求項1記載の散乱吸収体の内部情報
計測方法。 - 【請求項3】 所定波長の光を発生する光発生手段と、 前記所定波長の光を散乱吸収体に入射させる光入射手段
と、 3種以上の異なる光入射位置−光検出位置間距離に対応
する検出位置で、前記散乱吸収体の内部を拡散伝搬した
前記所定波長の光を検出して、3種以上の信号を取得す
る光検出手段と、 前記所定波長の光がそれぞれ拡散伝搬して前記検出位置
に到達するときの前記散乱吸収体内の拡散伝搬経路にお
ける散乱特性および吸収特性と前記検出信号との3つ以
上の連立関係に基づいて、前記検出信号を演算処理して
前記散乱吸収体の内部情報を導出する演算処理手段と、 を備えることを特徴とする散乱吸収体の内部情報計測装
置。 - 【請求項4】 前記所定波長の光は、散乱吸収体内部の
特定成分に対して吸収係数が異なる2種類以上の波長の
光であり、 前記演算処理手段は、1種類の前記所定波長による計測
によって得られる前記散乱吸収体の第1次内部情報の導
出処理と、前記2種以上の波長の光に対する前記第1次
内部情報を利用した前記散乱吸収体の前記特定成分に関
する情報である第2次内部情報の導出処理と、を実行す
る、 ことを特徴とする請求項3記載の散乱吸収体の内部情報
計測装置。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP02140094A JP3433498B2 (ja) | 1993-06-02 | 1994-02-18 | 散乱吸収体の内部情報計測方法及び装置 |
US08/252,227 US5517987A (en) | 1993-06-02 | 1994-06-01 | Method for measuring internal information in scattering medium and apparatus for the same |
DE69431326T DE69431326T2 (de) | 1993-06-02 | 1994-06-02 | Verfahren und Vorrichtung zur Messung innerer Eigenschaften eines Streumediums |
EP94303966A EP0627620B1 (en) | 1993-06-02 | 1994-06-02 | Method for measuring internal information in scattering medium and apparatus for the same |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP5-132211 | 1993-06-02 | ||
JP13221193 | 1993-06-02 | ||
JP02140094A JP3433498B2 (ja) | 1993-06-02 | 1994-02-18 | 散乱吸収体の内部情報計測方法及び装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH0749304A true JPH0749304A (ja) | 1995-02-21 |
JP3433498B2 JP3433498B2 (ja) | 2003-08-04 |
Family
ID=26358453
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP02140094A Expired - Lifetime JP3433498B2 (ja) | 1993-06-02 | 1994-02-18 | 散乱吸収体の内部情報計測方法及び装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US5517987A (ja) |
EP (1) | EP0627620B1 (ja) |
JP (1) | JP3433498B2 (ja) |
DE (1) | DE69431326T2 (ja) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09168531A (ja) * | 1995-08-24 | 1997-06-30 | Johnson & Johnson Medical Inc | 物質の1つまたはそれ以上の特性を定量的に決定する方法 |
JPH09264842A (ja) * | 1996-03-28 | 1997-10-07 | Technol Res Assoc Of Medical & Welfare Apparatus | 光散乱吸光媒質の吸光係数および光散乱係数の測定方法および測定装置 |
JP2001521147A (ja) * | 1997-10-16 | 2001-11-06 | ザ・リサーチ・ファンデーション・オブ・ステート・ユニバーシティ・オブ・ニューヨーク | 近赤外診療光学走査装置 |
US7286870B2 (en) | 1994-10-06 | 2007-10-23 | Hitachi, Ltd. | Optical system for measuring metabolism in a body and imaging method |
JP2010527635A (ja) * | 2006-05-30 | 2010-08-19 | ユニバーシティ オブ マサチューセッツ | 組織酸素化の測定 |
JP2014016235A (ja) * | 2012-07-09 | 2014-01-30 | Seiko Epson Corp | 光吸収係数分布推定装置、濃度測定装置及び光吸収係数分布推定装置の制御方法 |
JP2018084539A (ja) * | 2016-11-25 | 2018-05-31 | 株式会社分光応用技術研究所 | 光による特性測定システム及び光による特性測定方法 |
Families Citing this family (66)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5492118A (en) | 1993-12-16 | 1996-02-20 | Board Of Trustees Of The University Of Illinois | Determining material concentrations in tissues |
DE4417639A1 (de) * | 1994-05-19 | 1995-11-23 | Boehringer Mannheim Gmbh | Verfahren zur Bestimmung eines Analyten in einer biologischen Probe |
JPH10500338A (ja) * | 1994-05-19 | 1998-01-13 | ベーリンガー・マンハイム・ゲゼルシャフト・ミット・ベシュレンクテル・ハフツング | バイオロジカル・サンプル中の分析物を決定するための方法および装置 |
JP3433534B2 (ja) * | 1994-11-07 | 2003-08-04 | 浜松ホトニクス株式会社 | 散乱吸収体内の散乱特性・吸収特性の測定方法及び装置 |
DE4445683A1 (de) * | 1994-12-21 | 1996-06-27 | Boehringer Mannheim Gmbh | Verfahren zur Untersuchung eines streuenden Mediums mit intensitätsmoduliertem Licht |
US6031603A (en) * | 1995-06-09 | 2000-02-29 | Cybro Medical, Ltd. | Sensor, method and device for optical blood oximetry |
GB2304187B (en) * | 1995-08-07 | 1998-12-30 | Dia Stron Ltd | Translucency measurement |
DE19543020A1 (de) * | 1995-11-18 | 1997-05-22 | Boehringer Mannheim Gmbh | Verfahren und Vorrichtung zur Bestimmung von analytischen Daten über das Innere einer streuenden Matrix |
US5995856A (en) * | 1995-11-22 | 1999-11-30 | Nellcor, Incorporated | Non-contact optical monitoring of physiological parameters |
US6108576A (en) * | 1996-03-18 | 2000-08-22 | The Research Foundation Of City College Of New York | Time-resolved diffusion tomographic 2D and 3D imaging in highly scattering turbid media |
JP3662376B2 (ja) * | 1996-05-10 | 2005-06-22 | 浜松ホトニクス株式会社 | 内部特性分布の計測方法および装置 |
US5931779A (en) * | 1996-06-06 | 1999-08-03 | Wisconsin Alumni Research Foundation | Real-time in-vivo measurement of myoglobin oxygen saturation |
WO1998003847A2 (en) * | 1996-07-19 | 1998-01-29 | Mills Alexander K | Device for noninvasive determination of blood parameters |
JP3660761B2 (ja) * | 1996-10-03 | 2005-06-15 | 技術研究組合医療福祉機器研究所 | 散乱体の吸収情報の計測方法及び装置 |
US5935076A (en) * | 1997-02-10 | 1999-08-10 | University Of Alabama In Huntsville | Method and apparatus for accurately measuring the transmittance of blood within a retinal vessel |
US5776060A (en) * | 1997-02-20 | 1998-07-07 | University Of Alabama In Huntsville | Method and apparatus for measuring blood oxygen saturation within a retinal vessel with light having several selected wavelengths |
US6694157B1 (en) | 1998-02-10 | 2004-02-17 | Daedalus I , L.L.C. | Method and apparatus for determination of pH pCO2, hemoglobin, and hemoglobin oxygen saturation |
US6078833A (en) * | 1998-03-25 | 2000-06-20 | I.S.S. (Usa) Inc. | Self referencing photosensor |
US6199550B1 (en) * | 1998-08-14 | 2001-03-13 | Bioasyst, L.L.C. | Integrated physiologic sensor system |
US6850656B1 (en) * | 1998-10-07 | 2005-02-01 | Ecole Polytechnique Federale De Lausanne | Method and apparatus for measuring locally and superficially the scattering and absorption properties of turbid media |
US6678541B1 (en) * | 1998-10-28 | 2004-01-13 | The Governmemt Of The United States Of America | Optical fiber probe and methods for measuring optical properties |
US7047054B2 (en) * | 1999-03-12 | 2006-05-16 | Cas Medical Systems, Inc. | Laser diode optical transducer assembly for non-invasive spectrophotometric blood oxygenation monitoring |
AU7829500A (en) * | 1999-09-17 | 2001-04-17 | General Hospital Corporation, The | Calibration methods and systems for diffuse optical tomography and spectroscopy |
US20030139667A1 (en) * | 2000-04-13 | 2003-07-24 | Hewko Mark D. | Tissue viability/health monitor utilizing near infrared spectroscopy |
CA2417917C (en) * | 2000-08-04 | 2014-10-14 | Photonify Technologies, Inc. | Systems and methods for providing information concerning chromophores in physiological media |
EP1349493A2 (en) | 2001-01-12 | 2003-10-08 | The General Hospital Corporation | System and method for enabling simultaneous calibration and imaging of a medium |
JP4465271B2 (ja) | 2002-07-26 | 2010-05-19 | シーエーエス・メディカル・システムズ・インコーポレイテッド | 対象組織内の血液酸素飽和度を非侵襲的に決定する装置 |
DE10329915A1 (de) * | 2003-07-02 | 2005-03-17 | Manfred Holtkamp Elektronik Gmbh | Vorrichtungen und Verfahren für die Bestimmung einer zulässigen Bestrahlung der menschlichen Haut mit UV-Strahlung |
WO2006009910A2 (en) * | 2004-06-17 | 2006-01-26 | The Regents Of The University Of California | Time-resolved optometric fluorescence detection for skin diagnostics |
US7617628B2 (en) * | 2004-12-22 | 2009-11-17 | Smith & Wesson Corp. | Fire control mechanism for a firearm |
US7392611B2 (en) * | 2004-12-22 | 2008-07-01 | Smith & Wesson Corp. | Apparatus and method for firearm takedown |
US7506469B2 (en) * | 2004-12-22 | 2009-03-24 | Smith & Wesson Corp. | Firearm frame with configurable grip |
US7380362B2 (en) * | 2004-12-22 | 2008-06-03 | Smith & Wesson Corp. | Firearm extractor mechanism |
US7389719B2 (en) * | 2004-12-22 | 2008-06-24 | Smith & Wesson Corp. | Wire bushing for use with a firearm barrel |
US7600340B2 (en) * | 2004-12-22 | 2009-10-13 | Smith & Wesson Corp. | Locking apparatus for a firearm |
US7703230B2 (en) * | 2004-12-22 | 2010-04-27 | Smith & Wesson Corp. | Positive striker lock safety for use with a firearm |
EP1885235B1 (en) * | 2005-05-12 | 2013-12-18 | Cas Medical Systems, Inc. | Improved method for spectrophotometric blood oxygenation monitoring |
WO2006124696A1 (en) * | 2005-05-13 | 2006-11-23 | Children's Hospital Medical Center | Multi-wavelength spatial domain near infrared oximeter to detect cerebral hypoxia-ischemia |
JP4647447B2 (ja) * | 2005-09-20 | 2011-03-09 | 富士フイルム株式会社 | 試料分析装置 |
JP4662831B2 (ja) * | 2005-09-20 | 2011-03-30 | 富士フイルム株式会社 | 試料分析装置 |
WO2007060583A2 (en) * | 2005-11-28 | 2007-05-31 | Koninklijke Philips Electronics N.V. | Method and apparatus for determining concentrations of analytes in a turbid medium |
US9408538B2 (en) * | 2006-05-03 | 2016-08-09 | Triple Ring Technologies, Inc. | Method and apparatus for pressure sore detection |
WO2008118042A1 (en) * | 2007-03-23 | 2008-10-02 | St. Jude Medical Ab | Implantable medical device |
US8391942B2 (en) * | 2008-10-06 | 2013-03-05 | Cas Medical Systems, Inc. | Method and apparatus for determining cerebral desaturation in patients undergoing deep hypothermic circulatory arrest |
US20100105998A1 (en) * | 2008-10-28 | 2010-04-29 | Cas Medical Systems, Inc. | Method and apparatus for spectrophotometric based oximetry of spinal tissue |
US8296990B2 (en) | 2008-12-30 | 2012-10-30 | Smith & Wesson Corp. | Snap-on dovetail pistol sight |
US9057695B2 (en) * | 2009-09-24 | 2015-06-16 | Canon Kabushiki Kaisha | Apparatus and method for irradiating a scattering medium with a reconstructive wave |
US20110245637A1 (en) * | 2010-03-31 | 2011-10-06 | Nellcor Puritan Bennett Llc | Ambient light use in physiological sensors |
US7884933B1 (en) | 2010-05-05 | 2011-02-08 | Revolutionary Business Concepts, Inc. | Apparatus and method for determining analyte concentrations |
WO2012050847A2 (en) | 2010-09-28 | 2012-04-19 | Masimo Corporation | Depth of consciousness monitor including oximeter |
US9775545B2 (en) | 2010-09-28 | 2017-10-03 | Masimo Corporation | Magnetic electrical connector for patient monitors |
AU2011323310A1 (en) | 2010-11-03 | 2013-06-13 | University Of Washington Through Its Center For Commercialization | Determination of tissue oxygenation in vivo |
US8954130B2 (en) | 2010-12-17 | 2015-02-10 | Canon Kabushiki Kaisha | Apparatus and method for irradiating a medium |
US8817255B2 (en) | 2010-12-17 | 2014-08-26 | Canon Kabushiki Kaisha | Apparatus and method for irradiating a scattering medium |
US8976433B2 (en) | 2010-12-17 | 2015-03-10 | Canon Kabushiki Kaisha | Apparatus and method for irradiating a scattering medium |
US9304490B2 (en) | 2011-05-27 | 2016-04-05 | Canon Kabushiki Kaisha | Apparatus and method for irradiating a medium |
US9186112B2 (en) | 2012-05-03 | 2015-11-17 | Vioptix, Inc. | Tissue oximetry probe geometry for robust calibration and self-correction |
US20130310669A1 (en) * | 2012-05-20 | 2013-11-21 | Jerusalem College Of Technology | Pulmonary pulse oximetry method for the measurement of oxygen saturation in the mixed venous blood |
US8843186B2 (en) * | 2012-11-21 | 2014-09-23 | Folim G. Halaka | Non-invasive reagentless glucose determination |
US9693694B2 (en) | 2012-11-21 | 2017-07-04 | Folim G. Halaka | Cancer cell detection using dielectrophoretic dynamic light scattering (DDLS) spectroscopy |
US9848808B2 (en) | 2013-07-18 | 2017-12-26 | Cas Medical Systems, Inc. | Method for spectrophotometric blood oxygenation monitoring |
US10292654B2 (en) | 2013-12-16 | 2019-05-21 | Johnson & Johnson Consumer Inc. | Biomedical device, systems and methods having conductive elements |
USD763938S1 (en) | 2014-04-02 | 2016-08-16 | Cephalogics, LLC | Optical sensor array |
USD763939S1 (en) | 2014-04-02 | 2016-08-16 | Cephalogics, LLC | Optical sensor array liner with optical sensor array pad |
WO2016049103A1 (en) * | 2014-09-23 | 2016-03-31 | Halaka Folim G | Cancer cell detection using dielectrophoretic dynamic light scattering (ddls) spectroscopy |
WO2016057553A1 (en) | 2014-10-07 | 2016-04-14 | Masimo Corporation | Modular physiological sensors |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2612866B2 (ja) * | 1987-09-30 | 1997-05-21 | 株式会社島津製作所 | ヘモグロビン測定法 |
JPH01174940A (ja) * | 1987-12-29 | 1989-07-11 | Shimadzu Corp | ヘモグロビン測定法 |
CA1317669C (en) * | 1988-06-08 | 1993-05-11 | Randall L. Barbour | Method of imaging a random medium |
CA1331483C (en) * | 1988-11-02 | 1994-08-16 | Britton Chance | User-wearable hemoglobinometer for measuring the metabolic condition of a subject |
JPH06103257B2 (ja) * | 1988-12-19 | 1994-12-14 | 大塚電子株式会社 | 光散乱を用いた物質の吸光係数測定方法および装置 |
GB2228314B (en) * | 1989-02-16 | 1992-11-18 | Hamamatsu Photonics Kk | Examination apparatus |
EP0442011A1 (en) * | 1990-02-15 | 1991-08-21 | Hewlett-Packard GmbH | Sensor, apparatus and method for non-invasive measurement of oxygen saturation |
IE912141A1 (en) * | 1990-07-06 | 1992-01-15 | Kapsch Ag | Method and device for qualitative and quantitative¹determination of tissue-specific parameters of a biological¹tissue |
JP2932644B2 (ja) * | 1990-08-31 | 1999-08-09 | 株式会社島津製作所 | チトクロムオキシダーゼの測定方法及び測定装置 |
JP2539707B2 (ja) * | 1991-03-27 | 1996-10-02 | 大塚電子株式会社 | 吸光スペクトルの補正方法およびその方法を用いた光拡散物質の分光測定装置 |
CA2110716C (en) * | 1991-06-06 | 2003-10-28 | Gary D. Lewis | Optical cerebral oximeter |
US5277181A (en) * | 1991-12-12 | 1994-01-11 | Vivascan Corporation | Noninvasive measurement of hematocrit and hemoglobin content by differential optical analysis |
AU2245092A (en) * | 1991-12-31 | 1993-07-28 | Vivascan Corporation | Blood constituent determination based on differential spectral analysis |
US5297548A (en) * | 1992-02-07 | 1994-03-29 | Ohmeda Inc. | Arterial blood monitoring probe |
-
1994
- 1994-02-18 JP JP02140094A patent/JP3433498B2/ja not_active Expired - Lifetime
- 1994-06-01 US US08/252,227 patent/US5517987A/en not_active Expired - Lifetime
- 1994-06-02 DE DE69431326T patent/DE69431326T2/de not_active Expired - Fee Related
- 1994-06-02 EP EP94303966A patent/EP0627620B1/en not_active Expired - Lifetime
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7286870B2 (en) | 1994-10-06 | 2007-10-23 | Hitachi, Ltd. | Optical system for measuring metabolism in a body and imaging method |
US7440794B2 (en) | 1994-10-06 | 2008-10-21 | Hitachi, Ltd. | Optical system for measuring metabolism in a body and imaging method |
US7715904B2 (en) | 1994-10-06 | 2010-05-11 | Hitachi, Ltd. | Optical system for measuring metabolism in a body and imaging method |
US8050744B2 (en) | 1994-10-06 | 2011-11-01 | Hitachi, Ltd. | Optical system for measuring metabolism in a body and imaging method |
JPH09168531A (ja) * | 1995-08-24 | 1997-06-30 | Johnson & Johnson Medical Inc | 物質の1つまたはそれ以上の特性を定量的に決定する方法 |
JPH09264842A (ja) * | 1996-03-28 | 1997-10-07 | Technol Res Assoc Of Medical & Welfare Apparatus | 光散乱吸光媒質の吸光係数および光散乱係数の測定方法および測定装置 |
JP2001521147A (ja) * | 1997-10-16 | 2001-11-06 | ザ・リサーチ・ファンデーション・オブ・ステート・ユニバーシティ・オブ・ニューヨーク | 近赤外診療光学走査装置 |
JP2010101902A (ja) * | 1997-10-16 | 2010-05-06 | Research Foundation Of State Univ Of New York | 近赤外診療光学走査装置 |
JP2010527635A (ja) * | 2006-05-30 | 2010-08-19 | ユニバーシティ オブ マサチューセッツ | 組織酸素化の測定 |
JP2014016235A (ja) * | 2012-07-09 | 2014-01-30 | Seiko Epson Corp | 光吸収係数分布推定装置、濃度測定装置及び光吸収係数分布推定装置の制御方法 |
US9380968B2 (en) | 2012-07-09 | 2016-07-05 | Seiko Epson Corporation | Light absorption coefficient distribution estimation device, concentration measurement device, and method for controlling light absorption coefficient distribution estimation device |
JP2018084539A (ja) * | 2016-11-25 | 2018-05-31 | 株式会社分光応用技術研究所 | 光による特性測定システム及び光による特性測定方法 |
Also Published As
Publication number | Publication date |
---|---|
DE69431326T2 (de) | 2003-05-22 |
EP0627620B1 (en) | 2002-09-11 |
EP0627620A1 (en) | 1994-12-07 |
JP3433498B2 (ja) | 2003-08-04 |
US5517987A (en) | 1996-05-21 |
DE69431326D1 (de) | 2002-10-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3433498B2 (ja) | 散乱吸収体の内部情報計測方法及び装置 | |
JP3577335B2 (ja) | 散乱吸収体計測方法及び装置 | |
JP3433534B2 (ja) | 散乱吸収体内の散乱特性・吸収特性の測定方法及び装置 | |
EP0592200B1 (en) | Apparatus for measuring optical information in scattering medium and method therefor | |
JP3107914B2 (ja) | 散乱吸収体内部の吸収情報計測装置及び方法 | |
US10912504B2 (en) | Near-infrared spectroscopy and diffuse correlation spectroscopy device and methods | |
US5706821A (en) | Spectroscopic method for quantitatively determining the change of concentration of a light or other radiation absorbing compound in a medium which is interrogated through an intervening medium | |
US5137023A (en) | Method and apparatus for monitoring blood analytes noninvasively by pulsatile photoplethysmography | |
KR100612827B1 (ko) | 비 침습적인 헤모글로빈 농도와 산소 포화도 모니터링방법 및 장치 | |
CA1247397A (en) | Spectrophotometric method and apparatus for the non- invasive determination of glucose in body tissues | |
EP1013219B1 (en) | Living body information measuring apparatus, living body information measuring method, body fat measuring apparatus, body fat measuring method, and program recording medium | |
JPH07209177A (ja) | 散乱吸収体計測方法及び散乱吸収体計測装置 | |
US10925525B2 (en) | Combined pulse oximetry and diffusing wave spectroscopy system and control method therefor | |
JPH0961359A (ja) | 濃度測定装置 | |
EP1136811B1 (en) | Method and device for measuring internal information of scattering absorber | |
US6567165B1 (en) | Concentration measuring method and apparatus for absorption component in scattering medium | |
JP3903147B2 (ja) | 青果物の非破壊糖度測定装置 | |
JPH07120384A (ja) | 光計測方法および装置 | |
JP2690131B2 (ja) | 散乱電滋放射の減衰特性を用いて生体内の既知吸収度の組織色素の濃度を決定する方法及び装置 | |
Liu et al. | Spatial sensitivity of nirs tissue oxygenation measurement using a simplified instrument | |
JPH06181930A (ja) | 生体光計測装置 | |
Gravenstein et al. | Basic principles of optical radiation and some common applications in anesthesia | |
Chivel | New technology of cerebral oxygenation measurements by time-resolved spectroscopy |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313532 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080530 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090530 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100530 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110530 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110530 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120530 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130530 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130530 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140530 Year of fee payment: 11 |
|
EXPY | Cancellation because of completion of term |