JP3952275B2 - 生体光計測装置 - Google Patents
生体光計測装置 Download PDFInfo
- Publication number
- JP3952275B2 JP3952275B2 JP2002070672A JP2002070672A JP3952275B2 JP 3952275 B2 JP3952275 B2 JP 3952275B2 JP 2002070672 A JP2002070672 A JP 2002070672A JP 2002070672 A JP2002070672 A JP 2002070672A JP 3952275 B2 JP3952275 B2 JP 3952275B2
- Authority
- JP
- Japan
- Prior art keywords
- light
- measurement
- value
- signal
- hemoglobin
- 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
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Description
【発明の属する技術分野】
この発明は光を用いて生体内部情報を計測する生体光計測装置に関し、特に計測結果の信頼性の向上を図った生体光計測装置に関する。
【0002】
【従来の技術】
生体光計測装置は、光源から特定波長の光を被検体に照射し、被検体を透過した光或いは被検体の表層で反射した光を受光素子で検出し、その光量から血液循環、血行動態、ヘモグロビン変化などの生体情報を得る装置である。近年、このような光計測を利用して、複数の測定点を含む領域を検査し、その領域についての生体情報、具体的にはヘモグロビン動向を画像として表示したり、脳の活性領域などを計測するようにした光トポグラフィ装置が提案され、実用化されている(特開平9-98972号、特開平9-149903号など)。また生体光計測の臨床的な応用として、例えば頭部を計測対象とする場合、てんかん患者の焦点位置を同定する技術や麻酔深度を定量化する技術なども提案されている(特願2001-107070号など)。
【0003】
例えば麻酔深度の定量化においては、麻酔前と麻酔中のヘモグロビン濃度から麻酔効果の指標である麻酔深度を算出し、麻酔中に表示する。また活性領域の測定では、被検体に光刺激やタッピング等の所定のタスクを与え、タスク実行前のヘモグロビン濃度に対するタスク実行中のヘモグロビン濃度の変化から脳の活性領域を特定する。このように、これら生体光計測における計算や判定は、ヘモグロビン変化信号の大きさに基づいている。
【0004】
【発明が解決しようとする課題】
しかし、一般にヘモグロビン変化信号は微弱であるのに対し、装置の光源や光ファイバに起因して無視できないノイズや誤差成分が生じる場合がある。このため必ずしも上記計算や判定を精度よく行うことができず、その結果は誤診の可能性を内包している。
【0005】
そこで本発明は、生体光計測によって得られた信号に統計的処理を加えることによって結果の信頼性を高め、精度よい計測を行うことが可能な生体光計測装置を提供することを目的とする。
【0006】
【課題を解決するための手段】
本発明の生体光計測装置は、光源部と、前記光源部からの光を被検体の複数の照射位置に照射するとともに照射された光が生体を透過した光又は生体から反射した光を受光するプローブと、前記プローブが受光した光を計測位置毎の光強度として計測する光計測部と、前記光計測部が測定した光強度に対応する信号を用いて測定領域内のヘモグロビン変化信号を算出する信号処理部とを備えた生体光計測装置において、信号処理部は、時系列的に得られたヘモグロビン変化信号から、当該ヘモグロビン変化信号の統計的有意性を判定する評価量を求める手段を備えたものである。
【0007】
この生体光計測装置では、ヘモグロビン変化信号について統計的有意性を判定することにより、その判定結果又は統計的に有意と判定されたヘモグロビン信号を用いて信頼性の高い診断を行うことができる。
【0008】
本発明の生体光計測装置で算出する評価量は、各計測位置毎に求めた、統計的有意性を示す検定関数の合計として求めることができる。生体光計測では、測定領域で変化が起こっているか否かのみを診断する場合が多いので、測定領域に含まれる各計測位置の検定関数の合計として有意性を判定することにより、簡便に且つ高精度の診断が可能となる。
【0009】
さらに具体的には、評価量は、例えば、時系列的に得られたヘモグロビン信号を所定の時間幅で切り出した複数のサンプルとし、これら複数のサンプルの平均値及び標準偏差から検定統計量(T値)及び/又は有意確率(P値)を求め、これらT値又はP値を所定の基準を用いて0又は1を値とする検定関数に変換し、各計測位置について得た検定関数を合算したものを評価量とすることができる。ここで、所定の基準として、例えば、被検体に対しタスク無負荷の状態で計測したヘモグロビン変化信号から求めたP値の最小値を用いることができる。
【0010】
また本発明の生体光計測装置は、信号処理部が、複数の領域について同時に計測したヘモグロビン変化信号について、それぞれヘモグロビン変化信号の統計的有意性を判定する評価量を求める手段と、前記評価量に基き複数の領域のうち、統計的に有意な変化を生じている領域を判定する手段とを備える。
これにより、例えば右脳と左脳の活性を精度良く且つ容易に比較することが可能となる。
【0011】
さらに本発明の生体光計測装置は、信号処理部の演算により得られた評価量或いは評価量を用いた判定結果を表示する手段を備えているものである。
尚、本明細書において、「測定領域(或いは単に領域という場合もある)」はプローブによって計測される複数の計測点或いは計測位置を含む所定の広さの領域を意味する。また計測点或いは計測位置は、一つの光照射位置とそれと隣接する受光位置との中間の位置を意味する。
【0012】
【発明の実施の形態】
以下、本発明の生体光計測装置を、図面に示す実施形態に基づきさらに説明する。
【0013】
図1は、本発明の生体光計測装置の全体概要を示す図である。
この生体光計測装置101は、所定の波長の光を発生する光源部102と、被検体の検査部位を透過した光或いは検査部位で反射、散乱した光(以下、まとめて透過光という)を検出する105と、光源部102からの光を被検体の検査部位に誘導する照射用光ファイバ103と、検査部位からの透過光を光計測部105〜107に誘導する受光用光ファイバ103と、これら光ファイバの各先端を被検体の検査部位(測定領域)に当接させるために、各先端を着脱自在に固定する装着具(プローブ)104と、光計測部105〜107で計測した信号をもとに血中ヘモグロビン量などを表すヘモグロビン変化信号を作成し、画像化する信号処理部108とを備えている。
【0014】
光源部102は、可視光から赤外の波長領域内の所定の波長、例えば780nmや830nmなどの光を放射する半導体レーザ等の発光素子と、発光素子からの光に電気的に変調を与える変調器とを備えている。
【0015】
光計測部は、受光用光ファイバ103が誘導する光を光量に対応する電気信号に変換する光電変換素子105と、光電変換素子105からの電気信号を入力し、照射位置及び波長に対応した変調信号を選択的に検出するためのロックインアンプ106と、ロックインアンプ106からの信号をA/D変換するA/D変換器107とからなる。光電変換素子としては、例えばアバランシェフォトダイオードが使用される。A/D変換器107の出力は、信号処理部108に送られる。
【0016】
信号処理部108は、汎用のパーソナルコンピュータ等に構築することができ、光計測部からの信号を記録するメモリ109と、光計測部からの信号をもとにヘモグロビン変化信号を作成する中央演算処理装置(CPU)110、CPU110における演算処理に必要な条件やパラメータを記憶するとともに演算結果を記録するハードディスク111、ヘモグロビン変化信号をもとにヘモグロビン量などの変化や分布を表す画像等を表示するモニタ112を備えている。
【0017】
CPU110は、上述したヘモグロビン変化信号の計算において、以下述べるような統計的評価量の計算を行う。この評価量は、計測されたヘモグロビンの変化が有意なものか否か、つまり何らかの変化が生じているか否かの判定に用いられる。
【0018】
次に上述のような構成における生体光計測装置を用いた生体光計測について、脳活動の賦活判定を例に説明する。図2に、信号処理部108における処理のフローチャートを示す。
【0019】
まず、図1に示すプローブ104を被検体頭部に装着した状態で、近赤外光を照射するとともに生体からの透過光を計測する。この際、脳の誘発反応を生じさせる刺激、例えば光刺激を等間隔で連続して与える(ステップ201)。これによって計測された計測位置毎の光強度は、CPU110において、図3(a)に示すようなヘモグロビン変化信号301に変換される。なお、図3は複数の計測位置のうちの一つの位置(j番目)におけるヘモグロビン信号を示すものであり、図中、四角で囲ったタスクは上述の刺激を意味する。
【0020】
図3(a)に示すように、タスク時には、信号強度が高くなっていることが認められるが、この変化は微弱であり、ノイズにまぎれやすい。逆に、信号強度に変化がない計測位置についてもノイズによって信号強度が変化する場合がある。そこで、次に時系列的に得られるヘモグロビン変化信号をタスク試行時毎に切り出し(ステップ201)、有意性判断のための評価量を求める。図3(b)は、タスク毎に切り出したヘモグロビン信号302を示す。
【0021】
評価量の計算では、まずタスク試行時毎に切り出されたヘモグロビン信号について、タスク開始時点を基準とした同一潜時tにおける信号値mj,t,n(jは、測定位置、tは時間、nは何番目のタスク試行かを表す)の平均値Mj,t及び標準偏差Sj,tを次式により求める(ステップ202)。
【0022】
【数1】
次にこれら平均値及び標準偏差から、次式により、検定統計量T値(Tj,t)を求め、さらにT値から有意確率P値を求める(ステップ203)。
【0023】
【数2】
【0024】
P値は、自由度n-1のT分布において確率変数tが|Tj,t|以上もしくは|Tj,t|以下となる確率Pj,tであり、この値が少ないほど、何かが起こっている可能性が高いことを示す。これらT値、P値を全ての計測位置について行う。各計測位置について求めたP値を図4に示す。
【0025】
一方、タスク試行前に計測されたデータを用いて上述の方法と同様に、各測定位置毎にP値を求める(ステップ204)。そしてこのP値の最小値を基準αとして、タスク試行中に計測したデータの有意性を判定する(ステップ205)。即ち、タスク試行前に計測されたデータのP値の最小値は、装置性能に起因する検出限界つまりノイズと判断し、P値がこの最小値以下である潜時のデータを有意と判定する。具体的には、図4に示すグラフにおいて点線で示す基準α(タスク試行前P値の最小値)を超える場合を0、最小値以下を1とする。図5にこのように二値化した判定結果を示す。図示するように、この判定結果は、PがP>αのときに0、P≦αのときに1の値をとる関数である。これを判定関数λと呼ぶ。
【0026】
次に、各潜時において、この判定結果が1(即ち、P値が基準値以下)である計測位置の個数(計測チャンネル数)を求める(ステップ206、207)。これは、図5に示す各計測位置についての判定関数を合算することと同義である。こうして求めた計測チャンネル数をこの計測の評価量とし、
1)評価量が0である潜時においては、脳活動の賦活が検出されなかった
2)評価量が0以外の潜時においては、何らかの脳活動の賦活が検出されたと判定する。
【0027】
評価量算定結果の一例を図6に示す。本例は、麻酔時と非麻酔時に被検体に光刺激(タスク)を呈示し、そのとき測定されたヘモグロビン変化信号について有意性を判定したものであり、図中601は非麻酔時、602は麻酔時を示している。図示するように、麻酔時に計測されたヘモグロビン信号からタスク中に統計的有意と判定された計測位置の個数に比べ、非麻酔時に計測されたヘモグロビン変化信号からタスク中に統計的有意と判定された計測位置の個数は、明らかに増加している。これにより麻酔時に測定領域の活動が行われていないことが統計的に有意であると判定でき、信頼性の高い診断を行うことができる。
【0028】
また、このような判定結果を、図6に示すようなグラフとして、モニタ(図1、112)に表示することにより(ステップ208)、容易に視覚的に判定を行うことができる。
【0029】
このように本実施形態によれば、複数の計測位置についてそれぞれ求めたT値、P値から、その測定で得られた計測データ(ヘモグロビン変化信号)の統計的有意性を判断する機能を加えたことにより、その判定結果自体を用いて信頼性の高い診断(例えば何らかの脳の賦活の有無の診断)を行うことができる。この場合、生体光計測の目的が、測定領域で何らかの変化があることを診断することであれば、この判定結果の表示のみで変化の有無を高い信頼性を持って診断することができる。また、ヘモグロビン変化信号を用いて麻酔深度の算出等の演算を行う場合には、この判定結果で有意性があると判定されたヘモグロビン変化信号のみを用いることにより正確で誤差の少ない演算が可能となる。
【0030】
尚、以上の実施形態では、統計的な処理を行うためのサンプルの切り出し方として、ヘモグロビン変化信号をタスク試行時毎に切り出した場合を説明したが、例えば図7(a)に示すように、同一反応を行っていると考えられる所定の時間幅(T2−T1)を切り出してもよい。この場合には、データのサンプリング周波数をf[Hz]、T1〜T2間の測定時間点数をK(=(T2−T1)/f)とすると、次式により、平均値Mj、標準偏差Sjを求めることができる。
【0031】
【数3】
その後、検定統計量T値、有意確率P値を求めることは前述の実施形態と同様である。
【0032】
また図7(b)に示すように、所定の時間長内に、複数の区間を設定し、各々の区間において平均値、標準偏差を求めることも可能である。この場合には、タスク試行時に生じている事象をさらに詳細に観察することが可能になる。
【0033】
また以上の実施形態では、複数の計測位置をカバーする一つの領域について、その領域を計測したデータの統計的有意性を判定する場合を説明したが、本発明は、例えば、図8に示すように、右脳と左脳の対比など、複数の領域(測定領域1、測定領域2)について脳活性を比較する場合にも適用することが可能である。
【0034】
その場合、複数の領域についてそれぞれ上述のように評価量(統計的有意と判定された計測位置の個数)を求め、さらにこの評価量について、T値(或いはP値)を用いた差の検定を行ってもよい。図8は、検定の結果を示すもので、801は測定領域1にのみ顕著な活動が生じた場合、802は測定領域1、2ともに顕著な活動が認められなかった場合を示している。
【0035】
このように右脳と左脳の活性の違いを判定する技術は、例えば「”Non-invasive assessment of language dominance with near-infrared spectroscopic mapping”、Eiju Watanabe ら, Neuroscience Letters Vol.256, pp.49-52, 1998」に記載された言語優位半球の決定などに応用することができる。
【0036】
【発明の効果】
本発明によれば、信頼性の高い診断を行うことが可能な生体光計測装置が提供される。
【図面の簡単な説明】
【図1】本発明の生体光計測装置の全体概要を示す図
【図2】本発明の生体光計測装置による計測データ処理のフローを示す図
【図3】ヘモグロビン変化信号とその統計的処理の手順を示す図
【図4】各計測位置におけるP値算出結果の一例を示す図
【図5】各計測位置におけるP値から求めた検定関数を示す図
【図6】最終的な評価量の表示例を示す図
【図7】ヘモグロビン変化信号からサンプルを切り出す方法の他の実施例を示す図
【図8】本発明を適用した複数の領域の計測例を示す図
【符号の説明】
102・・・光源部、103・・・光ファイバ、105・・・光電変換素子(光計測部)、108・・・信号処理部、110・・・中央演算処理装置(CPU)
Claims (6)
- 光源部と、前記光源部からの光を被検体の複数の照射位置に照射するとともに照射された光が生体を透過した光又は生体から反射した光を受光するプローブと、前記プローブが受光した光を計測位置毎の光強度として計測する光計測部と、前記光計測部が測定した光強度に対応する信号を用いて測定領域内のヘモグロビン変化信号を算出する信号処理部とを備えた生体光計測装置において、
前記信号処理部は、時系列的に得られたヘモグロビン変化信号から、当該ヘモグロビン変化信号の統計的有意性を判定する評価量を求める手段を備え、各計測位置毎に、統計的有意性を示す検定関数を求め、各計測位置の検定関数の合計を前記評価量として求めることを特徴とする生体光計測装置。 - 光源部と、前記光源部からの光を被検体の複数の照射位置に照射するとともに照射された光が生体を透過した光又は生体から反射した光を受光するプローブと、前記プローブが受光した光を計測位置毎の光強度として計測する光計測部と、前記光計測部が測定した光強度に対応する信号を用いて測定領域内のヘモグロビン変化信号を算出する信号処理部とを備えた生体光計測装置において、
前記信号処理部は、時系列的に得られたヘモグロビン変化信号から、当該ヘモグロビン変化信号の統計的有意性を判定する評価量を求める手段を備え、各計測位置毎に、統計的有意性を示す検定関数を求め、検定関数の値が有意と判定される計測位置の数を前記評価量として求めることを特徴とする生体光計測装置。 - 前記信号処理部は、時系列的に得られたヘモグロビン信号を所定の時間幅で切り出した複数のサンプルとし、これら複数のサンプルの平均値及び標準偏差から検定統計量(T値)及び/又は有意確率(P値)を求め、これらT値又はP値を所定の基準を用いて0又は1を値とする検定関数に変換し、各計測位置について得た検定関数を合算したものを評価量とすることを特徴とする請求項1記載の生体光計測装置。
- 前記信号処理部は、前記所定の基準として、被検体に対しタスク無負荷の状態で計測したヘモグロビン変化信号から求めたP値の最小値を用いることを特徴とする請求項3に記載の生体光計測装置。
- 前記信号処理部は、複数の領域について同時に計測したヘモグロビン変化信号について、それぞれヘモグロビン変化信号の統計的有意性を判定する評価量を求める手段と、前記評価量に基き、複数の領域のうち統計的に有意な変化を生じている領域を判定する手段とを備えたことを特徴とする請求項1記載の生体光計測装置。
- 請求項1ないし5いずれか1項記載の生体光計測装置であって、前記信号処理部は、前記評価量又は前記評価量に基き判定した結果を表示する手段を備えたことを特徴とする生体光計測装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002070672A JP3952275B2 (ja) | 2002-03-14 | 2002-03-14 | 生体光計測装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002070672A JP3952275B2 (ja) | 2002-03-14 | 2002-03-14 | 生体光計測装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2003265442A JP2003265442A (ja) | 2003-09-24 |
JP3952275B2 true JP3952275B2 (ja) | 2007-08-01 |
Family
ID=29201175
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002070672A Expired - Fee Related JP3952275B2 (ja) | 2002-03-14 | 2002-03-14 | 生体光計測装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3952275B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11033225B2 (en) | 2015-05-27 | 2021-06-15 | Shimadzu Corporation | Brain activity feedback system |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009082595A (ja) * | 2007-10-02 | 2009-04-23 | Shimadzu Corp | 光生体測定装置 |
JP5324999B2 (ja) * | 2009-04-09 | 2013-10-23 | 株式会社日立製作所 | 生体光計測装置及び生体光計測方法 |
-
2002
- 2002-03-14 JP JP2002070672A patent/JP3952275B2/ja not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11033225B2 (en) | 2015-05-27 | 2021-06-15 | Shimadzu Corporation | Brain activity feedback system |
Also Published As
Publication number | Publication date |
---|---|
JP2003265442A (ja) | 2003-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3916985B2 (ja) | 生体光計測装置 | |
US6587702B1 (en) | Classification and characterization of tissue through features related to adipose tissue | |
US6456870B1 (en) | Non-invasive method of determining skin thickness and characterizing layers of skin tissue in vivo | |
US7254425B2 (en) | Method for detecting artifacts in data | |
KR102375159B1 (ko) | 레이저 유도 붕괴 스펙트럼 분석을 이용하는 진단 방법 및 이를 수행하는 진단 장치 | |
KR100490461B1 (ko) | 청진기 | |
JP4097522B2 (ja) | 生体光計測装置 | |
JP2006521869A (ja) | 光音響分析評価方法と装置 | |
US11076784B2 (en) | System for analyzing tissue | |
JP5674093B2 (ja) | 濃度定量装置及び濃度定量方法並びにプログラム | |
JP4642279B2 (ja) | 生体光計測装置 | |
JP2008049091A (ja) | 生体成分濃度測定方法 | |
JP2008546430A (ja) | 個体の被検体レベルの非侵襲的な測定方法 | |
JP4477403B2 (ja) | 生体光計測装置 | |
US20030040664A1 (en) | Method and apparatus for path normalization of light transport in tissue | |
JP3952275B2 (ja) | 生体光計測装置 | |
JP2019201760A (ja) | 血管検知装置及びその方法 | |
JP4411208B2 (ja) | 生体光計測装置 | |
JP4470681B2 (ja) | 光生体計測装置 | |
JP2019198547A (ja) | 血液成分測定方法、血液成分測定装置、および血液成分測定プログラム | |
JPWO2007066589A1 (ja) | 近赤外分光を用いた生活習慣病に関する検査・診断法および装置 | |
JP4071506B2 (ja) | 生体光計測装置 | |
US7565249B2 (en) | Method for the determination of a light transport parameter in a biological matrix | |
JP4272024B2 (ja) | 光生体計測装置 | |
KR101051406B1 (ko) | 혈관 내피세포 기능 측정 장치 및 방법 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20041007 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20060803 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070130 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070320 |
|
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: 20070417 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070418 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110511 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120511 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120511 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130511 Year of fee payment: 6 |
|
LAPS | Cancellation because of no payment of annual fees |