JPH0454439A - 生体計測方法および装置 - Google Patents
生体計測方法および装置Info
- Publication number
- JPH0454439A JPH0454439A JP16474990A JP16474990A JPH0454439A JP H0454439 A JPH0454439 A JP H0454439A JP 16474990 A JP16474990 A JP 16474990A JP 16474990 A JP16474990 A JP 16474990A JP H0454439 A JPH0454439 A JP H0454439A
- Authority
- JP
- Japan
- Prior art keywords
- light
- measurement
- subject
- living body
- measured
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000005259 measurement Methods 0.000 claims abstract description 73
- 239000000126 substance Substances 0.000 claims abstract description 23
- 230000008569 process Effects 0.000 claims abstract description 11
- 238000001514 detection method Methods 0.000 claims description 31
- 238000000691 measurement method Methods 0.000 claims description 17
- 230000000877 morphologic effect Effects 0.000 claims description 15
- 238000003384 imaging method Methods 0.000 claims description 10
- 210000003625 skull Anatomy 0.000 claims 2
- 230000003287 optical effect Effects 0.000 abstract description 53
- 230000002503 metabolic effect Effects 0.000 abstract description 4
- 238000003745 diagnosis Methods 0.000 abstract description 3
- 230000001678 irradiating effect Effects 0.000 abstract 1
- 238000010521 absorption reaction Methods 0.000 description 15
- 239000013307 optical fiber Substances 0.000 description 14
- 238000010586 diagram Methods 0.000 description 12
- 238000013425 morphometry Methods 0.000 description 12
- 230000003562 morphometric effect Effects 0.000 description 11
- 210000001519 tissue Anatomy 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 9
- 238000001228 spectrum Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 230000002123 temporal effect Effects 0.000 description 6
- 210000000988 bone and bone Anatomy 0.000 description 4
- 239000000470 constituent Substances 0.000 description 4
- 238000000149 argon plasma sintering Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 210000005013 brain tissue Anatomy 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 230000031700 light absorption Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 102000018832 Cytochromes Human genes 0.000 description 1
- 108010052832 Cytochromes Proteins 0.000 description 1
- 102000004190 Enzymes Human genes 0.000 description 1
- 108090000790 Enzymes Proteins 0.000 description 1
- 102000001554 Hemoglobins Human genes 0.000 description 1
- 108010054147 Hemoglobins Proteins 0.000 description 1
- 238000002835 absorbance Methods 0.000 description 1
- 239000011358 absorbing material Substances 0.000 description 1
- 230000008827 biological function Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 230000037149 energy metabolism Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 210000004884 grey matter Anatomy 0.000 description 1
- 239000002609 medium Substances 0.000 description 1
- 230000007102 metabolic function Effects 0.000 description 1
- 239000002207 metabolite Substances 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 210000004885 white matter Anatomy 0.000 description 1
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔産業上の利用分野〕
本発明は、生体の医学生理学的情報を非侵襲的に計測し
、画像として表示する生体計測方法および装置に関し、
特に光を用いて代謝循環等の生体機能の定量的な分布を
画像として提供するのに好適な生体計測方法および装置
に関する。
、画像として表示する生体計測方法および装置に関し、
特に光を用いて代謝循環等の生体機能の定量的な分布を
画像として提供するのに好適な生体計測方法および装置
に関する。
[従来の技術]
従来から、可視〜近赤外の光はヘモグロビン、チトクロ
ム等のエネルギー代謝酵棄に特異的に吸収されるが、他
の生体構成物質による吸収は小さいため、これら波長域
の光は生体のエネルギー代謝を計測するのに適している
ことが知られている(例えば、特開昭57−11523
2号公報参照)、また、XJilCT装置の計測原理を
応用し、上記波長域の光を用いて生体の断層面内の特定
物質の分布を計測し画像化する装置(光CT装置)が、
例えば、特開昭60−72542号公報に開示されてい
る。
ム等のエネルギー代謝酵棄に特異的に吸収されるが、他
の生体構成物質による吸収は小さいため、これら波長域
の光は生体のエネルギー代謝を計測するのに適している
ことが知られている(例えば、特開昭57−11523
2号公報参照)、また、XJilCT装置の計測原理を
応用し、上記波長域の光を用いて生体の断層面内の特定
物質の分布を計測し画像化する装置(光CT装置)が、
例えば、特開昭60−72542号公報に開示されてい
る。
ところで、X!ICTi置においては、被検体によるX
線減弱量の対数と吸収物質密度の関係は線形性が保たれ
ており、計測信号はいわゆるベアーランバート則に従う
、このため、XICT装置の計測信号はコンボリューシ
ョン(Convolusion)法と呼ばれる比較的単
純な手法で演算処理することで、断層画像に変換するこ
とができる。
線減弱量の対数と吸収物質密度の関係は線形性が保たれ
ており、計測信号はいわゆるベアーランバート則に従う
、このため、XICT装置の計測信号はコンボリューシ
ョン(Convolusion)法と呼ばれる比較的単
純な手法で演算処理することで、断層画像に変換するこ
とができる。
一方、可視から赤外の光は生体組織中では強い散乱を受
けるため、生体に光を照射しその透過光を計測する光C
Tにおいては、積比された計lIl!l値はベアーラン
バート則には従わず、1ICTi置と同様のコンボリュ
ーション法を用いた画像再生を行っても、被検体内の光
吸収量分布を正確に画像化することはできない。
けるため、生体に光を照射しその透過光を計測する光C
Tにおいては、積比された計lIl!l値はベアーラン
バート則には従わず、1ICTi置と同様のコンボリュ
ーション法を用いた画像再生を行っても、被検体内の光
吸収量分布を正確に画像化することはできない。
上述の問題、つまり、生体光計測における光散乱の影響
を除く手法としては、光の伝搬時間特性を用いた時間ゲ
ート法が特開昭63−206655号公報に開示されて
いる。
を除く手法としては、光の伝搬時間特性を用いた時間ゲ
ート法が特開昭63−206655号公報に開示されて
いる。
しかし、生体組織による光の散乱は強く、上述の時間ゲ
ート法によっても、その影響を十分に除くことはできな
い、特に、上記散乱光の影響は、生体内部の組織分布に
影響されるため、従来提案されている計測技法上の改善
のみでは、正確な画像を求めることはできないという問
題がある。
ート法によっても、その影響を十分に除くことはできな
い、特に、上記散乱光の影響は、生体内部の組織分布に
影響されるため、従来提案されている計測技法上の改善
のみでは、正確な画像を求めることはできないという問
題がある。
本発明は上記事情に鑑みてなされたもので、その目的と
するところは、従来の技術における上述の如き問題を解
消し、光CTM置における上述の散乱光の影響を効果的
に除き、正確な医学的診断を可能とする簡便な生体計測
方法およびそのための装置を提供することにある。
するところは、従来の技術における上述の如き問題を解
消し、光CTM置における上述の散乱光の影響を効果的
に除き、正確な医学的診断を可能とする簡便な生体計測
方法およびそのための装置を提供することにある。
本発明の上記目的は、生体に可視から近赤外の光を照射
して、その散乱透過光データから生体内における物質の
密度または化学状態の空間分布を計測し画像化する生体
計測方法において、可視から赤外の光を用いて被検体の
形状を計測し、該計測で得た形態情報を、生体内の物質
の密度または化学状態の空間分布を計測し画像化し表示
する過程において利用することを特徴とする生体計測方
法、および、上記方法を用いる生体計測装置であって、
複数の光入力部と光検出部を有し、これらの異なる組み
合せについて前記光入力部から光検出部への光伝搬時間
を計測して、前記光入力部と光検出部の相対位置関係を
計測して被検体の形状を計測する手段を有することを特
徴とする生体計測装置によって達成される。
して、その散乱透過光データから生体内における物質の
密度または化学状態の空間分布を計測し画像化する生体
計測方法において、可視から赤外の光を用いて被検体の
形状を計測し、該計測で得た形態情報を、生体内の物質
の密度または化学状態の空間分布を計測し画像化し表示
する過程において利用することを特徴とする生体計測方
法、および、上記方法を用いる生体計測装置であって、
複数の光入力部と光検出部を有し、これらの異なる組み
合せについて前記光入力部から光検出部への光伝搬時間
を計測して、前記光入力部と光検出部の相対位置関係を
計測して被検体の形状を計測する手段を有することを特
徴とする生体計測装置によって達成される。
〔作用1
本発明に係る生体計測方法においては、可視から赤外の
光を被写体に照射し、その透過光データから被写体の形
状を計測し、上記被写体の形態情報を得る。上記情報は
生体内の特定物質の濃度または化学状態の分布を画像化
する過程において。
光を被写体に照射し、その透過光データから被写体の形
状を計測し、上記被写体の形態情報を得る。上記情報は
生体内の特定物質の濃度または化学状態の分布を画像化
する過程において。
演算に必要な物理変数を規定し、上記画像処理時間の短
縮および高精度化を行うものである。
縮および高精度化を行うものである。
また、本発明に係る生体計測装置においては、複数の光
入力部と光検出部を有し、これらの異なる組み合せにつ
いて前記光入力部から光検出部への光伝搬時間を計測す
ることにより、前記光入力部と光検出部の相対位置関係
を計測して被検体の形状を計測し、この結果を利用する
ものである。
入力部と光検出部を有し、これらの異なる組み合せにつ
いて前記光入力部から光検出部への光伝搬時間を計測す
ることにより、前記光入力部と光検出部の相対位置関係
を計測して被検体の形状を計測し、この結果を利用する
ものである。
[実施例〕
以下、本発明の実施例を図面に基づいて詳細に説明する
。
。
第1図は、本発明の一実施例を示す生体計測装置(以下
、単に「装置jともいう)の全体構成を示すものである
0本装置は、光源部1.被検体計測部2、光検出部3.
検検出帯を数値化し記憶するデータ変換記憶部4.数値
化データを用いて被検体内の代謝機能分布を計算し画像
化するデータ処理部5.上記各部の動作を制御する計測
制御部6゜得られた画像を表示するための画像表示部7
から構成されている。
、単に「装置jともいう)の全体構成を示すものである
0本装置は、光源部1.被検体計測部2、光検出部3.
検検出帯を数値化し記憶するデータ変換記憶部4.数値
化データを用いて被検体内の代謝機能分布を計算し画像
化するデータ処理部5.上記各部の動作を制御する計測
制御部6゜得られた画像を表示するための画像表示部7
から構成されている。
次に、各部の機能の詳細を説明する。
光源部1は、被写体中の計測対象物質が特異的に吸収す
る波長とその近傍の、少なくとも2種以上の波長の時間
幅の狭いパルス光を、同時にまたは交互に計測制御部6
の指示に応じて放出する機能と、被検体の形態計測に最
適な上記波長とは異なる少なくとも一つの波長の光を放
出する機能とを有するものである。上述の如く、近傍の
2波長の光による減衰をそれぞれ計測し、両者による計
測値間の比を求めることにより、散乱に起因する減衰量
が除去でき、生体のような強い散乱体においても吸収量
の正確な計測が可能となる。また、光源部1は、吸収計
測に適した近傍の2個およびこれと異なる1個、計3個
以上の波長光を出力する機能を有する必要があり、例え
ば、血液中のヘモグロビンの酸素飽和度の計測において
は、酸素化ヘモグロビンと脱酸素化ヘモグロビンの等吸
収点(805nm)と、他の一点の波長(例えば、60
0n+m)と各波長の近傍の各1点の計4点の計測を行
う必要がある。このような場合には、上述の4個に、更
に形態計測用の1波長を加えた、計5個の波長光を出力
する機能を有することが必要となる。
る波長とその近傍の、少なくとも2種以上の波長の時間
幅の狭いパルス光を、同時にまたは交互に計測制御部6
の指示に応じて放出する機能と、被検体の形態計測に最
適な上記波長とは異なる少なくとも一つの波長の光を放
出する機能とを有するものである。上述の如く、近傍の
2波長の光による減衰をそれぞれ計測し、両者による計
測値間の比を求めることにより、散乱に起因する減衰量
が除去でき、生体のような強い散乱体においても吸収量
の正確な計測が可能となる。また、光源部1は、吸収計
測に適した近傍の2個およびこれと異なる1個、計3個
以上の波長光を出力する機能を有する必要があり、例え
ば、血液中のヘモグロビンの酸素飽和度の計測において
は、酸素化ヘモグロビンと脱酸素化ヘモグロビンの等吸
収点(805nm)と、他の一点の波長(例えば、60
0n+m)と各波長の近傍の各1点の計4点の計測を行
う必要がある。このような場合には、上述の4個に、更
に形態計測用の1波長を加えた、計5個の波長光を出力
する機能を有することが必要となる。
光源部1においては、第2図に示す如く、光源発光部8
から照射されたレーザ光ビームが、光走査部に入射する
。光走査部は回転ミラー9で構成されており、光源発光
部8からの出力光ビームの進行方向を偏向し、走査する
。上記光ビームは、複数の光ファイバーで構成された光
分配部10に入射し、回転ミラー9の走査角度に応じて
、複数の光ファイバーに順次、例えば、入力部番号■、
■。
から照射されたレーザ光ビームが、光走査部に入射する
。光走査部は回転ミラー9で構成されており、光源発光
部8からの出力光ビームの進行方向を偏向し、走査する
。上記光ビームは、複数の光ファイバーで構成された光
分配部10に入射し、回転ミラー9の走査角度に応じて
、複数の光ファイバーに順次、例えば、入力部番号■、
■。
■、・・・の順番で入力する。この結果、被検体計測部
2の光出力部からは、時系列的に、光ファイバーの配置
に応じて、生体への光照射が行われる。
2の光出力部からは、時系列的に、光ファイバーの配置
に応じて、生体への光照射が行われる。
この光走査角度と光源パルス発光のタイミングを制御す
ることにより、任意のファイバーに、設定されたタイミ
ングで設定数の光パルスを入力することが可能になる。
ることにより、任意のファイバーに、設定されたタイミ
ングで設定数の光パルスを入力することが可能になる。
このように、光走査と波長切り替えとを組み合わせるこ
とにより、設定されたスケジュールに従って、複数の波
長の光を順次高速に走査することができる。
とにより、設定されたスケジュールに従って、複数の波
長の光を順次高速に走査することができる。
上述の如く、光源部1から光ファイバーにより導かれた
光は、被検体計測部2において被検体に入射し、生体内
の代謝物質の光学特性に応じて吸収、散乱を受け、被検
体表面から放出される。被検体表面には複数の光検出器
が配置されており、生体の光学特性に応じた光が入力す
る。
光は、被検体計測部2において被検体に入射し、生体内
の代謝物質の光学特性に応じて吸収、散乱を受け、被検
体表面から放出される。被検体表面には複数の光検出器
が配置されており、生体の光学特性に応じた光が入力す
る。
本実施例に示す被検体計測部2は、光入力部および光検
出部を、計測断層面上に最適の位置で配置できる構成と
なっている。これは、第3図に示す如く、円形の固定リ
ング11から長さ可変の支持棒12で、被検体(斜線を
施した部分)表面に各光学系13を密着させることによ
って達成される。これにより、光学系を被写体表面に密
着できるため、人体頭部の如く、被検体ごとに形状が興
なる場合でも、最適の計測条件を実現できる。
出部を、計測断層面上に最適の位置で配置できる構成と
なっている。これは、第3図に示す如く、円形の固定リ
ング11から長さ可変の支持棒12で、被検体(斜線を
施した部分)表面に各光学系13を密着させることによ
って達成される。これにより、光学系を被写体表面に密
着できるため、人体頭部の如く、被検体ごとに形状が興
なる場合でも、最適の計測条件を実現できる。
上述の各支持部12の先端には、光出力または光入力用
の光学系が取り付けられており、光入力系に接続した光
ファイバーの一端は光源部1に、光検出系に接続した光
ファイバーの一端は光検出部3に入力している。光入力
系I3に入射した被検体からの透過散乱光は、光ファイ
バー14を通り光検出部3に導かれるが、07手法によ
り画像化を行う場合には、例えば、特開昭63−206
655号公報に記載された如き時間ゲートによる散乱線
除去を行なう必要がある。
の光学系が取り付けられており、光入力系に接続した光
ファイバーの一端は光源部1に、光検出系に接続した光
ファイバーの一端は光検出部3に入力している。光入力
系I3に入射した被検体からの透過散乱光は、光ファイ
バー14を通り光検出部3に導かれるが、07手法によ
り画像化を行う場合には、例えば、特開昭63−206
655号公報に記載された如き時間ゲートによる散乱線
除去を行なう必要がある。
このため、本実施例では、第4図に示す如く、光検出器
として、光の強度変化を高速に計測し記録できる装置1
5(以下、「光時間スペクトル計測装置」という)を用
いている0本光時間スペクトル計測装置】5は、それ自
身高速の時間応答特性を有しているため、入射光の時間
スペクトルを高速に計測、保存できるため、光シヤツタ
ーを用いることなく時間ゲートを設定することができる
。以下、その動作を説明する。
として、光の強度変化を高速に計測し記録できる装置1
5(以下、「光時間スペクトル計測装置」という)を用
いている0本光時間スペクトル計測装置】5は、それ自
身高速の時間応答特性を有しているため、入射光の時間
スペクトルを高速に計測、保存できるため、光シヤツタ
ーを用いることなく時間ゲートを設定することができる
。以下、その動作を説明する。
被検体計測部2から複数の光ファイバー14により導か
れた計測光は、光時間スペクトル計測装置15の受光面
16に適当な間隔で入射する。この入射光は受光面】6
において光電子を放出し、該光電子は印加電圧によって
螢光面18に向かい走行する。
れた計測光は、光時間スペクトル計測装置15の受光面
16に適当な間隔で入射する。この入射光は受光面】6
において光電子を放出し、該光電子は印加電圧によって
螢光面18に向かい走行する。
交流電場発生部17は、上述の電子の走行方向を高速に
変位させる。この結果、螢光面18には入力光の時間変
化に対応した分布画像が表示され、これをTVカメラ等
で読み出すことにより、光強度の時間変化を高速に計測
することができる。
変位させる。この結果、螢光面18には入力光の時間変
化に対応した分布画像が表示され、これをTVカメラ等
で読み出すことにより、光強度の時間変化を高速に計測
することができる。
ところで、高速のパルス光源8によって被検体計測部2
を介して光を照射された被検体の表面からは、光の走行
過程に応じた遅れ時間を持つ光が第5図(a)のような
時間分布で放出される。この光は、更に光ファイバー1
4を介して光時間スペクトル計測装置15の螢光面18
に、第5図(b)のような分布データを形成するにの分
布(b)のうち縦(Y)軸方向は時間を、横(X)軸方
向は各光ファイバーの光検出系13が配置されている被
検体表面近傍の位置を示すものである。このような分布
像をTVカメラ等で読み出し、計算機で処理することで
、各光検出系13に入力した光強度の時間変化が高速に
計測できる。上述の如き計測を、光パルスの照射位置、
すなわち、照射用光ファイバーを順次選択することで、
被検体の多方向からの光入射による光伝達特性を計測し
、これにより、被検体の3次元的な光学特性分布を画像
化するのに必要なすべてのデータを得ることができる。
を介して光を照射された被検体の表面からは、光の走行
過程に応じた遅れ時間を持つ光が第5図(a)のような
時間分布で放出される。この光は、更に光ファイバー1
4を介して光時間スペクトル計測装置15の螢光面18
に、第5図(b)のような分布データを形成するにの分
布(b)のうち縦(Y)軸方向は時間を、横(X)軸方
向は各光ファイバーの光検出系13が配置されている被
検体表面近傍の位置を示すものである。このような分布
像をTVカメラ等で読み出し、計算機で処理することで
、各光検出系13に入力した光強度の時間変化が高速に
計測できる。上述の如き計測を、光パルスの照射位置、
すなわち、照射用光ファイバーを順次選択することで、
被検体の多方向からの光入射による光伝達特性を計測し
、これにより、被検体の3次元的な光学特性分布を画像
化するのに必要なすべてのデータを得ることができる。
上述の光ファイバーの時系列的選択は、従来のX線CT
装置におけるエックス線源走査と同様の効果を持つもの
で、上記構成で得られた計測データをCTの画像再生理
論にもとづいて演算処理すれば、被検体の光学特性分布
を画像化することができる。この種の画像再生の原理お
よび方法に関しては1例えば、” I wage Re
construction fromP roject
ion” G 、 T 、 Herman、 1979
に詳細に述べられている。
装置におけるエックス線源走査と同様の効果を持つもの
で、上記構成で得られた計測データをCTの画像再生理
論にもとづいて演算処理すれば、被検体の光学特性分布
を画像化することができる。この種の画像再生の原理お
よび方法に関しては1例えば、” I wage Re
construction fromP roject
ion” G 、 T 、 Herman、 1979
に詳細に述べられている。
本実施例に示す被検体計測部2およびデータ処理部5は
、機能計測のための光学系のほかに被検体形状を計測す
るための光学系およびデータ処理機能が付加されている
ことを特徴とするものであることは、#J述の通りであ
る。
、機能計測のための光学系のほかに被検体形状を計測す
るための光学系およびデータ処理機能が付加されている
ことを特徴とするものであることは、#J述の通りであ
る。
例えば、本実施例の装置では、被検体上に配置した光学
素子の一部または全体を用いて形態計測を行う、形態計
測に用いる光源波長は、機能計測用の波長とは異なり、
生体構成物質による吸収の少ない波長を選択する。これ
により生体の生理変化による計測値の変動を小さくし、
更に信号の吸収による減弱を小さくしてS/Nの向上を
計ることができる。上述の形態計測用の波長としては、
生体による吸収が少なくかつ散乱減衰も比較的少ない9
00nm以上の近赤外光を用いることが望ましい、以下
、上述の形態計測部の詳細を述べる。
素子の一部または全体を用いて形態計測を行う、形態計
測に用いる光源波長は、機能計測用の波長とは異なり、
生体構成物質による吸収の少ない波長を選択する。これ
により生体の生理変化による計測値の変動を小さくし、
更に信号の吸収による減弱を小さくしてS/Nの向上を
計ることができる。上述の形態計測用の波長としては、
生体による吸収が少なくかつ散乱減衰も比較的少ない9
00nm以上の近赤外光を用いることが望ましい、以下
、上述の形態計測部の詳細を述べる。
形態計測部の機能は2つあり、その第1は被検体の形態
計測である。本計測で得られた形態情報は、機能分布の
画像化の過程で演算精度向上および演算時間短縮に利用
する。また、第2は検出器素子の位置計測であり、本実
施例のように被検体形状に密着し、各計測毎に変化する
光学系の配置を計測することができる。
計測である。本計測で得られた形態情報は、機能分布の
画像化の過程で演算精度向上および演算時間短縮に利用
する。また、第2は検出器素子の位置計測であり、本実
施例のように被検体形状に密着し、各計測毎に変化する
光学系の配置を計測することができる。
く形態計測法〉
以下、形態計測の手順を具体的に説明する。形態計測用
の光学系素子61,62,63.・・・は、第6図に示
す如く、被検体表面に密着させて配置する。上記素子の
うち6]、63.・・・は光入力用、62,64.・・
・は光出力用の光学系とする。入力部の素子61を選ん
で、これを基準点とし、以後の位置座標はすべてこれを
原点とする。例えば、3個の素子61,62.63の位
置関係を求めるには、簡単な幾何学的計算から、各素子
間の距離り、、 L、、 L、を求めればよい6そこで
、上述の基準点61から、時間幅が第7図にaで示され
るパルス光(約50psec)を入力する。
の光学系素子61,62,63.・・・は、第6図に示
す如く、被検体表面に密着させて配置する。上記素子の
うち6]、63.・・・は光入力用、62,64.・・
・は光出力用の光学系とする。入力部の素子61を選ん
で、これを基準点とし、以後の位置座標はすべてこれを
原点とする。例えば、3個の素子61,62.63の位
置関係を求めるには、簡単な幾何学的計算から、各素子
間の距離り、、 L、、 L、を求めればよい6そこで
、上述の基準点61から、時間幅が第7図にaで示され
るパルス光(約50psec)を入力する。
基準点61からの光は被検体内を伝搬し、配置した複数
の光検出部からファイバーを通し光検出部に入力する。
の光検出部からファイバーを通し光検出部に入力する。
この入射光は、生体構成物質の吸収が少ない波長を選択
しているため、被検体通過時には吸収の影響を殆んど受
けない。次に、この透過光を、先に述べた光時間スペク
トル計測装置15で計測する。上記検出器で計測される
信号は生体内で散乱され、パルス信号の時間幅が広がる
と同時に伝搬距離に応じて立ち上がり時間が遅れ、第7
図にbで示す形状となる。本計測における時間軸は、光
の伝達時間に対応するため、これに生体内の光伝搬速度
を積算すれば光の伝搬距離が求められる。生体組織はそ
の70%が水であるため、先の伝搬速度は水中とほぼ同
じ約2,3 x 10” cm/秒であり、伝搬時間を
1.(秒)とすると 伝搬距離 L=2,3xlO”xt、(cm)−°−°
(1)となる。
しているため、被検体通過時には吸収の影響を殆んど受
けない。次に、この透過光を、先に述べた光時間スペク
トル計測装置15で計測する。上記検出器で計測される
信号は生体内で散乱され、パルス信号の時間幅が広がる
と同時に伝搬距離に応じて立ち上がり時間が遅れ、第7
図にbで示す形状となる。本計測における時間軸は、光
の伝達時間に対応するため、これに生体内の光伝搬速度
を積算すれば光の伝搬距離が求められる。生体組織はそ
の70%が水であるため、先の伝搬速度は水中とほぼ同
じ約2,3 x 10” cm/秒であり、伝搬時間を
1.(秒)とすると 伝搬距離 L=2,3xlO”xt、(cm)−°−°
(1)となる。
例えば、人間の頭部断面の直径は15c+++程度であ
り、上述の関係から、人体頭部断面を透過する光の伝搬
時間は600p秒程度となり、前述の光時間スペクトル
計測装置15(時間分解能約10P秒)で計測できる。
り、上述の関係から、人体頭部断面を透過する光の伝搬
時間は600p秒程度となり、前述の光時間スペクトル
計測装置15(時間分解能約10P秒)で計測できる。
ところで、各計測点で計測された時間分解データ波形の
うち、立ち上がり部分の信号は被検体内各種経路を通過
し検出器に到達した伝搬光のうち最短の伝搬路を通過し
た光の信号であり、その遅れ時間1dは光入力部と光検
出部を結ぶ直線距離に比例している。つまり、上記各検
出器の時間分解データの立ち上がり時刻1dを求めれば
、式(1)から、入力点から検出素子への距離が求まる
ことになる。
うち、立ち上がり部分の信号は被検体内各種経路を通過
し検出器に到達した伝搬光のうち最短の伝搬路を通過し
た光の信号であり、その遅れ時間1dは光入力部と光検
出部を結ぶ直線距離に比例している。つまり、上記各検
出器の時間分解データの立ち上がり時刻1dを求めれば
、式(1)から、入力点から検出素子への距離が求まる
ことになる。
上述の方法で、まず、第6図の61を発光点として、6
1−62の距離つまりLlが求められる1次に、63を
発光点としてり、が求まる。上記手順により。
1−62の距離つまりLlが求められる1次に、63を
発光点としてり、が求まる。上記手順により。
L、およびり、は求まるが、61と63は両者とも発光
点のため、61−63の距離り、は計測できない。そこ
で、検出器63でも光検出が可能となるよう、光入力用
と光検出用の光ファイバーを一体で構成する。これによ
り、61を発光点、63を検出点としてり、が求められ
る0以上の手順で、L、、 L、、 L、が求まり、結
果的に、3点(61,62,63)のpl標が求まる。
点のため、61−63の距離り、は計測できない。そこ
で、検出器63でも光検出が可能となるよう、光入力用
と光検出用の光ファイバーを一体で構成する。これによ
り、61を発光点、63を検出点としてり、が求められ
る0以上の手順で、L、、 L、、 L、が求まり、結
果的に、3点(61,62,63)のpl標が求まる。
以下、三角測量と同様の原理で、61,62.63以外
のすべての素子位装置を求めることができる。
のすべての素子位装置を求めることができる。
上述の手順によって、配置されているすべての光検出部
と光入力部の位置座標が求められる。更に、これらの各
素子が被検体表面に密着していることから、jis図に
示す如く、各検出器の位置座標を清らかな曲線で結ぶ二
とにより、被検体の外部形状が求められる。
と光入力部の位置座標が求められる。更に、これらの各
素子が被検体表面に密着していることから、jis図に
示す如く、各検出器の位置座標を清らかな曲線で結ぶ二
とにより、被検体の外部形状が求められる。
ところで、上記計測信号(第7図のb)は、伝搬域の光
散乱による広がりに応じて、パルス時間幅が広がる特性
、つまり、検出光の時間計測信号の形状は、被検体の散
乱特性の情報を有する。そこで、上記計測信号の時間特
性から、光の散乱情報を得ることが可能で、これを利用
して被検体内部の形態情報を得ることができる。以下、
その手順を示す。
散乱による広がりに応じて、パルス時間幅が広がる特性
、つまり、検出光の時間計測信号の形状は、被検体の散
乱特性の情報を有する。そこで、上記計測信号の時間特
性から、光の散乱情報を得ることが可能で、これを利用
して被検体内部の形態情報を得ることができる。以下、
その手順を示す。
計測空間を、第10図(a)に示す如く、適当な大きさ
の画素に分割し、各画素の光学変数として、散乱係数”
1+吸収係数m、を設定する。ここで、lは素子番号を
示している。また、上記計測で得られた時間分解信号は
、被検体の吸収の影響を殆んど受けないため、m 4=
Oとみなすことができ、真の変数は散乱係数S、のみと
なる。
の画素に分割し、各画素の光学変数として、散乱係数”
1+吸収係数m、を設定する。ここで、lは素子番号を
示している。また、上記計測で得られた時間分解信号は
、被検体の吸収の影響を殆んど受けないため、m 4=
Oとみなすことができ、真の変数は散乱係数S、のみと
なる。
まず、先に求めた被検体の外部形状から、被検体内の生
体組織の構造分布を、解剖学的に仮定する。本実施例に
おいて被検体とする人体頭部ではその構成組織は、骨、
脳組R(白質、灰白質)および脳室に分けられる。具体
的には、例えば、第10図(b)に示す如く、計測断面
の表面の内側に約] Oimの厚さの骨を配置する。次
に、骨で覆われた領域の内側に、標準的な形状の脳室を
配置し、残りの部分に脳組織を配置する。これら各組織
の散乱係数の標準値は、実際の組織の計測データから求
めることが可能であり、この値を、先に仮定した組織分
布に応じて、上記各画素に設定する。これにより散乱係
数s(の初期値が設定される。上記初期条件を基に、イ
テレーション法により、正確な散乱係数分布 Siを求
めることができる。以下、具体的な手順を説明する。
体組織の構造分布を、解剖学的に仮定する。本実施例に
おいて被検体とする人体頭部ではその構成組織は、骨、
脳組R(白質、灰白質)および脳室に分けられる。具体
的には、例えば、第10図(b)に示す如く、計測断面
の表面の内側に約] Oimの厚さの骨を配置する。次
に、骨で覆われた領域の内側に、標準的な形状の脳室を
配置し、残りの部分に脳組織を配置する。これら各組織
の散乱係数の標準値は、実際の組織の計測データから求
めることが可能であり、この値を、先に仮定した組織分
布に応じて、上記各画素に設定する。これにより散乱係
数s(の初期値が設定される。上記初期条件を基に、イ
テレーション法により、正確な散乱係数分布 Siを求
めることができる。以下、具体的な手順を説明する。
まず、上記初期条件下で、先に求めた配置の光入力部を
光源としたときの、検出部での時間分解データを求める
。これは、各画素での散乱確率を上記散乱係数により与
える、モンテカルロ演算を行うことにより可能となる。
光源としたときの、検出部での時間分解データを求める
。これは、各画素での散乱確率を上記散乱係数により与
える、モンテカルロ演算を行うことにより可能となる。
このようなモンテカルロ計算の原理および具体的な方法
については、例えば、J 0M、Maarek、 ”S
imulation of lasertomosco
py in a heterogeneous bio
logical med−ium″Med & Bio
l & Comput、]]986,24.pp、40
7−41に詳細に述べられている。次に、これと先の計
測で求めた時間分解データとを比較し、その差が最小と
なるように散乱係数分布を変位させる。
については、例えば、J 0M、Maarek、 ”S
imulation of lasertomosco
py in a heterogeneous bio
logical med−ium″Med & Bio
l & Comput、]]986,24.pp、40
7−41に詳細に述べられている。次に、これと先の計
測で求めた時間分解データとを比較し、その差が最小と
なるように散乱係数分布を変位させる。
このような演算を繰り返すイテレーション法により、正
確な被検体の内部形状を求めることができる。なお、上
記イテレーション法において、変位させる数値パラメー
タとして、組織分布の形状のみを、散乱係数51として
は各組織の標準値を用いれば、上記演算時間を短縮する
ことができる。
確な被検体の内部形状を求めることができる。なお、上
記イテレーション法において、変位させる数値パラメー
タとして、組織分布の形状のみを、散乱係数51として
は各組織の標準値を用いれば、上記演算時間を短縮する
ことができる。
以上の手順により、被検体の計測面内における形状が求
まる。この形態情報と、先に求めた検出器位置座標を用
いることにより、以後の機能分布の画像化処理を、効率
良く正確に行うことが可能となる。
まる。この形態情報と、先に求めた検出器位置座標を用
いることにより、以後の機能分布の画像化処理を、効率
良く正確に行うことが可能となる。
(機能画像計測の手順〉
次に、機能分布の画像化における各計測部の機能および
データ処理系の動作を説明する、なお、形態画像および
検出器位置座標は、既に前記手順で得られているものと
する。
データ処理系の動作を説明する、なお、形態画像および
検出器位置座標は、既に前記手順で得られているものと
する。
まず、第2図の計測部の光分配部10に、順次光パルス
を照射し■、■、■、・・・に示す順番で照射する。こ
のとき、各光検出器には、第7図す、で示されるような
時間スペクトルを持った光信号が順次入力する。これら
の計測値は、一連の計測が終了するまで、記憶部に保存
される。ところで、これらのデータには生体内の散乱光
が多く含まれるため、このままでは、画像化に用いるこ
とはできない。そこで、前述の特開昭63−20665
5号公報に示される如き時間ゲートを用いて散乱光を除
く0本実施例においては、第5図(a)に示した如く、
上記時間スペクトルデータのうち、散乱光の少ない時間
スペクトル前部のデータを設定した時間幅で加算し、以
後の演算に用いることで同様の効果を得ている。
を照射し■、■、■、・・・に示す順番で照射する。こ
のとき、各光検出器には、第7図す、で示されるような
時間スペクトルを持った光信号が順次入力する。これら
の計測値は、一連の計測が終了するまで、記憶部に保存
される。ところで、これらのデータには生体内の散乱光
が多く含まれるため、このままでは、画像化に用いるこ
とはできない。そこで、前述の特開昭63−20665
5号公報に示される如き時間ゲートを用いて散乱光を除
く0本実施例においては、第5図(a)に示した如く、
上記時間スペクトルデータのうち、散乱光の少ない時間
スペクトル前部のデータを設定した時間幅で加算し、以
後の演算に用いることで同様の効果を得ている。
ここで、時間ゲート幅(本実施例ではデータの選択幅)
を狭くするほど散乱光の混入を少なくして、空間分解能
を向上することができるが、一方で、信号量が減るとい
う問題がある。このため、時間ゲート幅は被検体の大き
さ、用いた光の透過率により最適値が異なり、例えば、
人体頭部の場合は]0Opsec程度が最適な幅となる
。
を狭くするほど散乱光の混入を少なくして、空間分解能
を向上することができるが、一方で、信号量が減るとい
う問題がある。このため、時間ゲート幅は被検体の大き
さ、用いた光の透過率により最適値が異なり、例えば、
人体頭部の場合は]0Opsec程度が最適な幅となる
。
このように、スペクトルデータを設定した時間幅で時間
方向に加算し、画像再生のためのプロジェクションデー
タを作成する。この結果、光照射部の数をN8、光検出
部の数をN、とすると、上記計測でえられるデータはN
、XN、個の数値データとなる。
方向に加算し、画像再生のためのプロジェクションデー
タを作成する。この結果、光照射部の数をN8、光検出
部の数をN、とすると、上記計測でえられるデータはN
、XN、個の数値データとなる。
〈画像再生処理〉
次に、上記計測で得られた機能計測データの組(N、X
N、個の数値データ)と、事前に計測した形態画像を用
いて、機能分布の画像再生を行なう。
N、個の数値データ)と、事前に計測した形態画像を用
いて、機能分布の画像再生を行なう。
以下、画像再生処理の流れを具体的に説明する。
第10図(a)に示す如く、計測領域をNXMの画素に
分割する。この画素に、既に求めである計測部の光入力
部と光検出部の位置座標を与える。また、同図に示す如
く、上記各画素に光学的な変数として、散乱係数S□お
よび吸収係数miを与える(1は画素の番号)。なお、
形態計測部で計測した形態画像により、被検体のS工は
求められているので、上記各光学変数のうちSiは定数
として与える。この結果、求めるべき変数は脳組織内の
吸収係数mlのみになる。更に、被検体の外部形状から
、被検体の外側の吸収係数m1を Oと与えることがで
きる。
分割する。この画素に、既に求めである計測部の光入力
部と光検出部の位置座標を与える。また、同図に示す如
く、上記各画素に光学的な変数として、散乱係数S□お
よび吸収係数miを与える(1は画素の番号)。なお、
形態計測部で計測した形態画像により、被検体のS工は
求められているので、上記各光学変数のうちSiは定数
として与える。この結果、求めるべき変数は脳組織内の
吸収係数mlのみになる。更に、被検体の外部形状から
、被検体の外側の吸収係数m1を Oと与えることがで
きる。
上述の如き処理により、本来、各画素に与えられていた
変数の組み合わせのうち、多くの部分を定数として設定
できるため、以後の画像化が容易になり、演算に要する
時間も短縮できる。
変数の組み合わせのうち、多くの部分を定数として設定
できるため、以後の画像化が容易になり、演算に要する
時間も短縮できる。
次に、上述の変数mlおよび設定された定数を基に、吸
収係数の分布を求める方法を説明する。
収係数の分布を求める方法を説明する。
ある発光点Skに光パルスを入力したときの計測値は、
第10図(b)に示される如く、各検出点と発光点j(
j=1.2,3.・・・)の組み合わせA□となる。こ
のような計測をまず第]の波長で行い、このときの計測
値をA1□とする。次に、この近傍の第2の波長で得た
計測値をA、□とし、上記各データから X□=−1og(A、kJ/A、。) ・・・
・(2)を求める。ここで1時間ゲートを通過した光は
往体内をほぼ直進しているので、kとJの組み合わせか
ら光が通過する領域を、第1O図(C)のように求める
二とができる。そこで、この通過領域に含まれるすべて
の画素番号の集合をRとする。このとき、上記計測値よ
りえられた吸光度xk、は、X0=Σ(m1+ s 1
)、JCR・・・・(3)となる。
第10図(b)に示される如く、各検出点と発光点j(
j=1.2,3.・・・)の組み合わせA□となる。こ
のような計測をまず第]の波長で行い、このときの計測
値をA1□とする。次に、この近傍の第2の波長で得た
計測値をA、□とし、上記各データから X□=−1og(A、kJ/A、。) ・・・
・(2)を求める。ここで1時間ゲートを通過した光は
往体内をほぼ直進しているので、kとJの組み合わせか
ら光が通過する領域を、第1O図(C)のように求める
二とができる。そこで、この通過領域に含まれるすべて
の画素番号の集合をRとする。このとき、上記計測値よ
りえられた吸光度xk、は、X0=Σ(m1+ s 1
)、JCR・・・・(3)となる。
ところで、これらのうちSiは形態計測により既知であ
るので、各計測における真の吸収による減弱Yk−は、 Y、、=Xk、−Σ(Si)、1ER−(4)により求
めることができる。また、画素】が骨の領域に含まれる
場合には、吸収がほとんど無いためm1=Oと置くこと
ができる。このような計測を、すべてのkjの組み合わ
せに関して行う。これにより、吸収値m工のプロジェク
ションの組み合わせ(Y、、、に=]〜N、、j=1〜
N、)が求められるので、従来のX線0丁におけるC7
画像再生のアルゴリズムを適用して、吸収値m□の分布
が求められる。
るので、各計測における真の吸収による減弱Yk−は、 Y、、=Xk、−Σ(Si)、1ER−(4)により求
めることができる。また、画素】が骨の領域に含まれる
場合には、吸収がほとんど無いためm1=Oと置くこと
ができる。このような計測を、すべてのkjの組み合わ
せに関して行う。これにより、吸収値m工のプロジェク
ションの組み合わせ(Y、、、に=]〜N、、j=1〜
N、)が求められるので、従来のX線0丁におけるC7
画像再生のアルゴリズムを適用して、吸収値m□の分布
が求められる。
以上の演算によって、第1O図(d)のような特定の2
波長の光吸収係数の差の空間分布(つまり特定の代謝物
質の空間分布)を示す画像が得られる。
波長の光吸収係数の差の空間分布(つまり特定の代謝物
質の空間分布)を示す画像が得られる。
また、上記機能分布画像は、第I0図(e)に示すよう
に、既に得られている形態画像と重ねて表示することに
よりより、有効な診断情報を与える二とが可能となる。
に、既に得られている形態画像と重ねて表示することに
よりより、有効な診断情報を与える二とが可能となる。
なお、上記実施例の形態計測において、第9図に示す如
き素子を用いて発光点の近傍で光検出を行うことにより
、骨による反射光を捉え、この反射光の時間遅れから骨
の厚みを計測しても良い。
き素子を用いて発光点の近傍で光検出を行うことにより
、骨による反射光を捉え、この反射光の時間遅れから骨
の厚みを計測しても良い。
上記実施例によれば、被検体の形状および検出器位置を
求めて、これを機能分布画像の再生に利用することによ
り、簡単な処理で、正確な機能分布画像を効率良く求め
ることが可能になる。
求めて、これを機能分布画像の再生に利用することによ
り、簡単な処理で、正確な機能分布画像を効率良く求め
ることが可能になる。
なお、上記実施例は本発明の一例を示すものであり、本
発明はこれに限定されるべきものではないことは言うま
でもない。
発明はこれに限定されるべきものではないことは言うま
でもない。
以上、詳細に述べた如く、本発明によれば、生体に可視
から近赤外の光を照射して、その散乱透過光データから
生体内における物質の密度または化学状態の空間分布を
計測し画像化する生体計測方法において、可視から赤外
の光を用いて被検体の形状を計測し、該計測で得た形態
情報を、生体内の物質の密度または化学状態の空間分布
を計測し画像化し表示する過程において利用するように
したので、光CT装置における散乱光の影響を効果的に
除き、正確な医学的診断を可能とする簡便な生体計測方
法を実現できるという顕著な効果を奏するものである。
から近赤外の光を照射して、その散乱透過光データから
生体内における物質の密度または化学状態の空間分布を
計測し画像化する生体計測方法において、可視から赤外
の光を用いて被検体の形状を計測し、該計測で得た形態
情報を、生体内の物質の密度または化学状態の空間分布
を計測し画像化し表示する過程において利用するように
したので、光CT装置における散乱光の影響を効果的に
除き、正確な医学的診断を可能とする簡便な生体計測方
法を実現できるという顕著な効果を奏するものである。
また、上述の計測方法を用いる生体計測装置において、
複数の光入力部と光検出部を有し、これらの異なる組み
合せについて前記光入力部から光検出部への光伝搬時間
を計測して、前記光入力部と光検出部の相対位置関係を
計測して被検体の形状を計測する手段を設けることによ
り、正確な機能分布画像を効率良く求めることが可能な
生体計測装置を実現できるという顕著な効果を奏するも
のである。
複数の光入力部と光検出部を有し、これらの異なる組み
合せについて前記光入力部から光検出部への光伝搬時間
を計測して、前記光入力部と光検出部の相対位置関係を
計測して被検体の形状を計測する手段を設けることによ
り、正確な機能分布画像を効率良く求めることが可能な
生体計測装置を実現できるという顕著な効果を奏するも
のである。
第1図は本発明の一実施例を示す生体光計性装置の構成
図、第2図は光源部の構成および動作を示す図、第3図
は被検体計測部の構成を示す図、第4図は光検出部の動
作を示す図、第5図は計測信号の波形を示す図、第6図
は形態計測の原理を示す図、第7図は形態計測における
信号波形を示す図、第8図は被検体の外形計測法を示す
図、第9図は光学系の一実施例を示す図、第10図は画
像化の処理過程を示す図である。 1:光源部、2:被検体計測部、3:光検出部、4:デ
ータ変換記憶部、5:データ処理部、6:計測制御部、
7:表示部、8:光源発光部、9:光走査部、10:光
用1部、11:固定リング、12:支持部、】3コ光学
系、14:光ファイバー、15:光時間スペクトル計測
装置。 第 図 第 図 第 図 R’) 第 図 時間 七贈招 第 図 第 図 第 図(その2)
図、第2図は光源部の構成および動作を示す図、第3図
は被検体計測部の構成を示す図、第4図は光検出部の動
作を示す図、第5図は計測信号の波形を示す図、第6図
は形態計測の原理を示す図、第7図は形態計測における
信号波形を示す図、第8図は被検体の外形計測法を示す
図、第9図は光学系の一実施例を示す図、第10図は画
像化の処理過程を示す図である。 1:光源部、2:被検体計測部、3:光検出部、4:デ
ータ変換記憶部、5:データ処理部、6:計測制御部、
7:表示部、8:光源発光部、9:光走査部、10:光
用1部、11:固定リング、12:支持部、】3コ光学
系、14:光ファイバー、15:光時間スペクトル計測
装置。 第 図 第 図 第 図 R’) 第 図 時間 七贈招 第 図 第 図 第 図(その2)
Claims (1)
- 【特許請求の範囲】 1、生体に可視から近赤外の光を照射して、その散乱透
過光データから生体内における物質の密度または化学状
態の空間分布を計測し画像化する生体計測方法において
、可視から赤外の光を用いて被検体の形状を計測し、該
計測で得た形態情報を、生体内の物質の密度または化学
状態の空間分布を計測し画像化し表示する過程において
利用することを特徴とする生体計測方法。 2、計測した散乱透過光データから生体内の物質の密度
または化学状態の空間分布を画像化し表示する過程にお
いて、計測対象領域を複数の画素に分割し、該分割した
各画素に複数の物理量を変数として与え、前記形態情報
を用いて前記変数のうち少なくとも1個の物理変数の変
域を制限し、演算処理を行なうことを特徴とする請求項
1記載の生体計測方法。 3、前記形態情報に加えて、解剖学的形態情報をも利用
することを特徴とする請求項1または2記載の生体計測
方法。 4、前記被検体の形状計測に用いる計測光と、生体内の
物質の密度または化学状態の空間分布の計測に用いる計
測光との波長が異なるものであることを特徴とする請求
項1〜3のいずれかに記載の生体計測方法。 5、前記被検体の形状計測に用いる計測光が、波長90
0nm以上の光であることを特徴とする請求項4記載の
生体計測方法。 6、請求項1〜5のいずれかに記載の生体計測方法を用
いて生体計測を行うための生体計測装置であつて、複数
の光入力部と光検出部を有し、これらの異なる組み合せ
について前記光入力部から光検出部への光伝搬時間を計
測して、前記光入力部と光検出部の相対位置関係を計測
して被検体の形状を計測する手段を有することを特徴と
する生体計測装置。 7、前記光入力部と光検出部が、被検体表面に位置付け
可能に構成されたことを特徴とする請求6記載の生体計
測装置。 8、被検体として生体頭部を対象とし、該頭部表面から
光を入力し、前記被検体表面に配置した複数の光検出器
で前記入力光の散乱反射光を計測し演算することで頭蓋
骨の厚さを計測し、該頭蓋骨の厚さを画像化の過程にお
いて形態情報として利用することを特徴とする請求項6
または7記載の生体計測装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP16474990A JPH0454439A (ja) | 1990-06-22 | 1990-06-22 | 生体計測方法および装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP16474990A JPH0454439A (ja) | 1990-06-22 | 1990-06-22 | 生体計測方法および装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JPH0454439A true JPH0454439A (ja) | 1992-02-21 |
Family
ID=15799190
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP16474990A Pending JPH0454439A (ja) | 1990-06-22 | 1990-06-22 | 生体計測方法および装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPH0454439A (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05306994A (ja) * | 1992-04-30 | 1993-11-19 | Shimadzu Corp | 光断層像測定装置 |
JPH0643096A (ja) * | 1992-04-24 | 1994-02-18 | Becton Dickinson & Co | 光散乱および吸収の時間分解測定法による血液培養瓶のバクテリア検出方法および装置 |
JP2001521147A (ja) * | 1997-10-16 | 2001-11-06 | ザ・リサーチ・ファンデーション・オブ・ステート・ユニバーシティ・オブ・ニューヨーク | 近赤外診療光学走査装置 |
-
1990
- 1990-06-22 JP JP16474990A patent/JPH0454439A/ja active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0643096A (ja) * | 1992-04-24 | 1994-02-18 | Becton Dickinson & Co | 光散乱および吸収の時間分解測定法による血液培養瓶のバクテリア検出方法および装置 |
JPH05306994A (ja) * | 1992-04-30 | 1993-11-19 | Shimadzu Corp | 光断層像測定装置 |
JP2001521147A (ja) * | 1997-10-16 | 2001-11-06 | ザ・リサーチ・ファンデーション・オブ・ステート・ユニバーシティ・オブ・ニューヨーク | 近赤外診療光学走査装置 |
JP2010101902A (ja) * | 1997-10-16 | 2010-05-06 | Research Foundation Of State Univ Of New York | 近赤外診療光学走査装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6282438B1 (en) | Optical system for measuring metabolism in a body and imaging method | |
JP3112025B2 (ja) | 生体計測装置 | |
JP3310782B2 (ja) | 吸収物質濃度の空間分布画像化装置 | |
EP0904011B1 (en) | Apparatus for imaging microvascular blood flow | |
US8280494B2 (en) | Apparatus and method to measure a spectroscopic characteristic in an object | |
US8355131B2 (en) | Device and method for acquiring image data from a turbid medium | |
JP3142079B2 (ja) | 光ct装置 | |
JP2005114473A (ja) | 光検出方法及び生体光計測装置 | |
JP2007528500A (ja) | 蛍光タンパク質を使用する断層撮影イメージングのための方法およびシステム | |
US5774223A (en) | Optical measuring method and an optical measuring apparatus for determining the internal structure of an object | |
JP3164112B2 (ja) | 生体計測装置 | |
JPH0454439A (ja) | 生体計測方法および装置 | |
JPH10246697A (ja) | 光学的検査方法及び光学的検査装置 | |
JP2956777B2 (ja) | 生体光計測方法およびそのための装置 | |
JPH07120384A (ja) | 光計測方法および装置 | |
RU91517U1 (ru) | Устройство диффузионной оптической томографии | |
JP3276947B2 (ja) | 吸収物質濃度空間分布画像装置 | |
RU2787360C1 (ru) | Устройство дифференциальной диагностики глубины термических поражений кожного покрова человека с активным инфракрасным зондированием | |
JPH0641911B2 (ja) | 光断層像画像化装置 | |
JPH0620458B2 (ja) | 高指向性結像素子及び高指向性結像装置 | |
JPH10221248A (ja) | 光学的検査方法及び光学的検査装置 | |
Han et al. | Interferometric diffusing speckle contrast imaging (iDSCI) system for monitoring regional cerebral blood flow | |
WO2018165992A1 (zh) | 一种多通道脑功能成像装置及方法 | |
Tabatabaei et al. | Thermophotonic lock-in imaging: a novel early caries detection and imaging modality | |
Rao et al. | Images of biological soft tissue using synchrotron X-ray and laser CT systems |