(57)【特許請求の範囲】
【請求項１】
近赤外光を被検体に照射する光照射部と、
前記被検体に照射され、前記被検体内を経由した光を検出する光検出部と、
前記光検出部が検出した前記光から取得される情報を演算処理して出力する演算処理部とを有し、
前記演算処理部は、前記光から取得される前記被検体における血液中のオキシヘモグロビン変化量Ｏ及びデオキシヘモグロビン変化量Ｄの各情報から、下記式（６）及び（７）で表される、前記オキシヘモグロビン変化量Ｏにおける流速性変動の血流変化成分ＯＦ及び前記デオキシヘモグロビン変化量Ｄにおける流速性変動の血流変化成分ＤＦ、並びに／又は前記オキシヘモグロビン変化量Ｏにおける容積性変動の血液変化成分ＯＶ及び前記デオキシヘモグロビン変化量Ｄにおける容積性変動の血流変化成分ＤＶを算出して出力することを特徴とする生体光計測装置。
0005641355.tif 000016
ただし、前記式（６）及び（７）において、ｋＦ及びｋＶは、ｋＦ≠ｋＶを満たす数値であるとともに、ｋＦが−１≦ｋＦ＜０の負の数値であり、ｋＶが０＜ｋＶの正の数値である。
【請求項２】
ｋＦは、−０．８５＜ｋＦ＜−０．３５を満たし、ｋＶは、０＜ｋＶ≦１を満たす数値である請求項１に記載の生体光計測装置。
【請求項３】
近赤外線を被検体の頭部に照射する請求項１から２のいずれかに記載の生体光計測装置。
【請求項４】
光照射部から照射される近赤外光を被検体に導く照射プローブと、前記被検体内を経由した光を光検出部に導く受光プローブとを有する請求項１から３のいずれかに記載の生体光計測装置。
【請求項５】
光検出部が複数の受光プローブから検出される光をそれぞれ検出するとともに、
演算処理部が検出された前記各光から取得される情報をそれぞれ演算処理して出力する請求項４に記載の生体光測定装置。
【請求項６】
被検体の表面に沿って取付け可能なプローブホルダを有し、前記プローブホルダに照射プローブ及び受光プローブが保持される請求項４から５のいずれかに記載の生体光測定装置。
【請求項７】
近赤外光を被検体に照射する光照射部と、前記被検体に照射され、前記被検体内を経由した光を検出する光検出部と、前記光検出部が検出した前記光から取得される情報を演算処理して出力する演算処理部とを有する生体光計測装置に用いられ、
前記演算処理部に対して、前記光から取得される前記被検体における血液中のオキシヘモグロビン変化量Ｏ及びデオキシヘモグロビン変化量Ｄの各情報から、下記式（６）及び（７）で表される、前記オキシヘモグロビン変化量Ｏにおける流速性変動の血流変化成分ＯＦ及び前記デオキシヘモグロビン変化量Ｄにおける流速性変動の血流変化成分ＤＦ、並びに／又は前記オキシヘモグロビン変化量Ｏにおける容積性変動の血液変化成分ＯＶ及び前記デオキシヘモグロビン変化量Ｄにおける容積性変動の血流変化成分ＤＶを算出して出力させることを特徴とするプログラム。
0005641355.tif 000017
ただし、前記式（６）及び（７）において、ｋＦ及びｋＶは、ｋＦ≠ｋＶを満たす数値であるとともに、ｋＦが−１≦ｋＦ＜０の負の数値であり、ｋＶが０＜ｋＶの正の数値である。
【請求項８】
近赤外光を被検体に照射する光照射ステップと、
前記被検体に照射され、前記被検体内を経由した光を検出する光検出ステップと、
前記光検出ステップで検出した前記光から取得される情報を演算処理して出力する演算処理ステップとを有し、
前記演算処理ステップは、前記光から取得される前記被検体における血液中のオキシヘモグロビン変化量Ｏ及びデオキシヘモグロビン変化量Ｄの各情報から、下記式（６）及び（７）で表される、前記オキシヘモグロビン変化量Ｏにおける流速性変動の血流変化成分ＯＦ及び前記デオキシヘモグロビン変化量Ｄにおける流速性変動の血流変化成分ＤＦ、並びに／又は前記オキシヘモグロビン変化量Ｏにおける容積性変動の血液変化成分ＯＶ及び前記デオキシヘモグロビン変化量Ｄにおける容積性変動の血流変化成分ＤＶを算出して出力することを特徴とする生体光計測方法。
0005641355.tif 000018
ただし、前記式（６）及び（７）において、ｋＦ及びｋＶは、ｋＦ≠ｋＶを満たす数値であるとともに、ｋＦが−１≦ｋＦ＜０の負の数値であり、ｋＶが０＜ｋＶの正の数値である。
【請求項９】
流速性変動の血流変化成分ＯＦ、ＤＦの各情報に基づき、神経活動に伴う被検体の大脳灰白質層における血流変動を測定する請求項８に記載の生体光測定方法。
【請求項１０】
容積性変動の血液変化成分ＯＶ、ＤＶの各情報に基づき、被検体の頭皮層における血流変動を測定する請求項８から９のいずれかに記載の生体光測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、近赤外光を利用して生体活動由来の信号を測定する生体光計測装置、プログラム及び生体光計測方法に関する。
【背景技術】
【０００２】
近赤外光を利用して生体活動由来の信号を測定する生体光計測方法としては、ｓｉｎｇｌｅ  ｄｉｓｔａｎｃｅ  ｐｒｏｂｅ  ａｒｒａｎｇｅｍｅｎｔによるｆＮＩＲＳ計測法（以下、ｓｉｎｇｌｅ  ｄｉｓｔａｎｃｅ法）が知られている。この方法は、図１に示すように、光照射部から照射された近赤外光を照射プローブ１を介して頭部１０に導き、その導かれた光を頭部１０に近接して配した受光プローブ２を介して光検出部で検出することとしている。なお、図中に示される符号は、３が頭皮、４が頭蓋骨、５が脳髄液、６が大脳灰白質、７が大脳白質層である。
【０００３】
従来型のｓｉｎｇｌｅ  ｄｉｓｔａｎｃｅ法によるｆＮＩＲＳ計測データには、１）脳機能信号、２）脳機能信号以外の生理的信号、３）プローブ接触の不具合等に由来する信号変動などの成分が混在している（非特許文献１参照）。プローブ接触の不具合に由来する信号変動は、多くの場合、ベースラインの瞬時的なシフトやノイズの極端な多さなど、他の成分とは異なる時間的挙動を示すため比較的容易に識別できる。これに対して、脳機能信号以外の生理活動に由来する信号のうちのある部分は、脳機能信号に非常に類似した時間的挙動を示す。こうした信号成分は、具体的には、表層組織等における血流変動に由来しているが、例えば、指先の反復対向動作（タッピング）や認知課題を行っているときには、該当する機能領野で脳活動が生じると同時に全身性の血流変化が生じ、これに起因する信号変化が頭部の広汎な領域で観測されることが知られている（非特許文献２，３参照）。このような場合、チャンネル位置と脳機能部位の位置相関に関する先験的知識を持たない限り、我々は観測している信号変化が脳機能由来のものであるかどうかを判断することができないという問題がある。したがって、何らかの方法で脳活動に伴う光吸収係数の変化とそれ以外を分離することが必要である。
【０００４】
これまで脳機能信号の抽出のために幾つかの手法が試みられているが、それらは以下のような理由で多少なりとも問題を含むものであった。
先ず、時間的変化の特徴を用いて信号分離するアダプティブ・フィルタリング等の手法（非特許文献４参照）は、脳機能信号が他の成分と相関が低いことを前提としているが、実際には全身性血流信号の一部に課題遂行と同期して生じる脳機能信号以外の成分（アーティファクト）が存在することから、この手法では、脳機能信号とそれ以外を十分に分離することができない。
また、ＩＣＡ（Ｉｎｄｅｐｅｎｄｅｎｔ  Ｃｏｍｐｏｎｅｎｔ  Ａｎａｌｙｓｉｓ；独立成分分析、非特許文献５参照）も多チャンネル計測データに含まれる多成分の信号間の独立性を根拠に成分分離を行う手法であるため、やはり脳機能信号との同時性や波形の類似性の高いこうしたアーティファクトの分離には向かない。
また、同じく多チャンネル計測データにＰＣＡ（Ｐｒｉｎｃｉｐａｌ  Ｃｏｍｐｏｎｅｎｔ  Ａｎａｌｙｓｉｓ；主成分分析）を適用し、グローバルに生じる信号成分と局所成分を分離する手法が提案されているが（非特許文献２参照）、この手法は、脳活動が単一箇所の領野のみで生じる場合に限定され、かつ、その場合においても当該領野以外の広汎な場所にプローブを設置する必要があり、更に、それらのプローブのいずれかで例えばプローブコンタクトの不良などがあれば、解析精度に大きな影響を及ぼすなどの問題があった。
【０００５】
一方で、実験の課題内容や配列を工夫することによって、課題遂行により誘起されるアーティファクトと脳機能信号とを分離する試みもなされている。
例えば、上記のタッピングの例では、課題配列で「安静−タッピング−安静」のように安静区間に課題を挿入した場合には課題区間では脳活動と同時に生じる全身性血流変化の影響を観察してしまう可能性を排除できない。これに対して、「右手タッピング−左手タッピング−右手タッピング−左手タッピング」のような課題配列を用いる手法（非特許文献６参照）では、全区間を通じて全身性血流の影響が重畳するので、片側一次運動野の計測では脳機能活動による信号変化のみが検出される蓋然性が高い。
しかしながら、こうした手法にも幾つかの問題がある。即ち、このような実験デザインを通じた脳機能信号の抽出は、測定者に一般に高度な実験デザインスキルが要求される。また、この種の実験パラダイムでは、用意した対照課題で当該の脳機能活動以外のすべての生理活動がターゲット課題と同様に生じることを暗黙の前提としているが、その前提を厳密に証明することはほとんどの場合困難である。
【０００６】
我々は、こうした困難を克服するため、これまでにｍｕｌｔｉ  ｄｉｓｔａｎｃｅ  ｐｒｏｂｅ  ａｒｒａｎｇｅｍｅｎｔの方法を提案している（以下、ｍｕｌｔｉ  ｄｉｓｔａｎｃｅ法という。特許文献１及び非特許文献３参照）。
この提案方法は、図１に示すｓｉｎｇｌｅ  ｄｉｓｔａｎｃｅ法のプローブ配置に加え、更に図２に示すように、参照用プローブ８を受光プローブ２よりも照射プローブ１に近い位置に配する手法であり（例えば、照射プローブ１−受光プローブ２間の距離３０ｍｍに対し、照射プローブ１−参照用プローブ８間の距離を２０ｍｍとする）、光源−検出器間の距離が短く比較的頭皮３から浅い位置を経過する光を参照用プローブ８から検出し、その検出結果を心拍、呼吸、自律神経活動に伴う内因性の変動など脳機能信号以外の信号を含むものとして、光検出プローブ２から検出され、光源−検出器間の距離が長く比較的頭皮３から深い位置を経過する光の信号との差分をとり、脳機能由来の信号のみを検出することとする。
【０００７】
この提案方法によれば、上記問題のほとんどを克服して脳機能信号抽出を実現できる。
しかしながら、特別なプローブ配置を必要とし、市販装置でこの手法を実行しようとした場合に２つの問題を生じさせる。
一つは、市販装置のプローブ固定具の大きさの都合で、上記のプローブ配置を一般ユーザーが容易に実現できない点である。もう一つは、二つの検出プローブでの実測光量は、光源−検出器間距離の相違に基づき、相互で１０倍程度も異なるため、検出器のダイナミックレンジの調節機構等を十分に備えていない市販装置では、測定が困難となる点である。
したがって、こうした問題を解決し、市販装置のほとんどで装置改変をすることなく、一般ユーザーが取得できるデータのみから脳活動由来の信号のみを抽出する手法としては、何らの提案もされていないというのが現状である。
【発明の概要】
【発明が解決しようとする課題】
【０００８】
本発明は、従来における前記諸問題を解決し、以下の目的を達成することを課題とする。即ち、本発明は、市販のｆＮＩＲＳ装置を装置改変をすることなく、一般ユーザーが取得できるデータを用い、生体活動における特定信号のみを抽出する生体光測定装置、プログラム及び生体光測定方法を提供することを目的とする。
【０００９】
前記課題を解決するために、我々は、種々の経験的知見に基づき、生体組織中での血流変動に起因する近赤外光の信号変化を容積性の変動成分と流速性の変動成分という二つの様態に類別して眺めることを提案する。そして、この提案の有用性を以下のように検証した。
先ず、前記二つの様態に係る変動成分は、血流中のオキシ、デオキシ変化量の間に正／負の比例関係が成り立つ成分としてそれぞれ特徴づけられる。次に、ｆＮＩＲＳ装置の観測信号をこれらの変動成分の線形和として記述する。そして、これらの変動成分の比例関係の係数を決定する。そして、アーティファクト課題及びアーティファクトを伴う脳機能課題時の脳機能活動に伴うヘモグロビン変化を、前記提案手法及び大脳灰白質のヘモグロビン変化を計測できるｍｕｌｔｉ  ｄｉｓｔａｎｃｅ法（特許文献１及び非特許文献３参照）を用いて計算し、両者の一致度を比較する。この一致度が高ければ、前記提案手法を有用と考えることができる。
その結果、容積性の変動成分と流速性の変動成分という二つの様態に類別することにより、生体組織中での血流変動に起因する近赤外光の信号変化から、被検体の生体活動による特定信号のみを抽出することができ、その抽出手法に係る演算処理を市販の装置に適用することで、一般ユーザーが取得できるデータのみから生体活動における特定信号のみを抽出する生体光測定装置、プログラム及び生体光測定方法を提供することができることの知見を得た。
【課題を解決するための手段】
【００１０】
本発明は、前記知見に基づくものであり、前記課題を解決するための手段としては、以下の通りである。即ち、
＜１＞  近赤外光を被検体に照射する光照射部と、前記被検体に照射され、前記被検体内を経由した光を検出する光検出部と、前記光検出部が検出した前記光から取得される情報を演算処理して出力する演算処理部とを有し、前記演算処理部は、前記光から取得される前記被検体における血液中のオキシヘモグロビン変化量Ｏ及びデオキシヘモグロビン変化量Ｄの各情報から、下記式（６）及び（７）で表される、前記オキシヘモグロビン変化量Ｏにおける流速性変動の血流変化成分ＯＦ及び前記デオキシヘモグロビン変化量Ｄにおける流速性変動の血流変化成分ＤＦ、並びに／又は前記オキシヘモグロビン変化量Ｏにおける容積性変動の血液変化成分ＯＶ及び前記デオキシヘモグロビン変化量Ｄにおける容積性変動の血流変化成分ＤＶを算出して出力することを特徴とする生体光計測装置。
0005641355.tif 000002
ただし、前記式（６）及び（７）において、ｋＦ及びｋＶは、ｋＦ≠ｋＶを満たす数値であるとともに、ｋＦが−１≦ｋＦ＜０の負の数値であり、ｋＶが０＜ｋＶの正の数値である。
＜２＞  ｋＦは、−０．８５＜ｋＦ＜−０．３５を満たし、ｋＶは、０＜ｋＶ≦１を満たす数値である前記＜１＞に記載の生体光計測装置。
＜３＞  近赤外線を被検体の頭部に照射する前記＜１＞から＜２＞のいずれかに記載の生体光計測装置。
＜４＞  光照射部から照射される近赤外光を被検体に導く照射プローブと、前記被検体内を経由した光を光検出部に導く受光プローブとを有する前記＜１＞から＜３＞のいずれかに記載の生体光計測装置。
＜５＞  光検出部が複数の受光プローブから検出される光をそれぞれ検出するとともに、演算処理部が検出された前記各光から取得される情報をそれぞれ演算処理して出力する前記＜４＞に記載の生体光測定装置。
＜６＞  被検体の表面に沿って取付け可能なプローブホルダを有し、前記プローブホルダに照射プローブ及び受光プローブが保持される前記＜４＞から＜５＞のいずれかに記載の生体光測定装置。
＜７＞  近赤外光を被検体に照射する光照射部と、前記被検体に照射され、前記被検体内を経由した光を検出する光検出部と、前記光検出部が検出した前記光から取得される情報を演算処理して出力する演算処理部とを有する生体光計測装置に用いられ、前記演算処理部に対して、前記光から取得される前記被検体における血液中のオキシヘモグロビン変化量Ｏ及びデオキシヘモグロビン変化量Ｄの各情報から、下記式（６）及び（７）で表される、前記オキシヘモグロビン変化量Ｏにおける流速性変動の血流変化成分ＯＦ及び前記デオキシヘモグロビン変化量Ｄにおける流速性変動の血流変化成分ＤＦ、並びに／又は前記オキシヘモグロビン変化量Ｏにおける容積性変動の血液変化成分ＯＶ及び前記デオキシヘモグロビン変化量Ｄにおける容積性変動の血流変化成分ＤＶを算出して出力させることを特徴とするプログラム。
0005641355.tif 000003
ただし、前記式（６）及び（７）において、ｋＦ及びｋＶは、ｋＦ≠ｋＶを満たす数値であるとともに、ｋＦが−１≦ｋＦ＜０の負の数値であり、ｋＶが０＜ｋＶの正の数値である。
＜８＞  近赤外光を被検体に照射する光照射ステップと、前記被検体に照射され、前記被検体内を経由した光を検出する光検出ステップと、前記光検出ステップで検出した前記光から取得される情報を演算処理して出力する演算処理ステップとを有し、前記演算処理ステップは、前記光から取得される前記被検体における血液中のオキシヘモグロビン変化量Ｏ及びデオキシヘモグロビン変化量Ｄの各情報から、下記式（６）及び（７）で表される、前記オキシヘモグロビン変化量Ｏにおける流速性変動の血流変化成分ＯＦ及び前記デオキシヘモグロビン変化量Ｄにおける流速性変動の血流変化成分ＤＦ、並びに／又は前記オキシヘモグロビン変化量Ｏにおける容積性変動の血液変化成分ＯＶ及び前記デオキシヘモグロビン変化量Ｄにおける容積性変動の血流変化成分ＤＶを算出して出力することを特徴とする生体光計測方法。
0005641355.tif 000004
ただし、前記式（６）及び（７）において、ｋＦ及びｋＶは、ｋＦ≠ｋＶを満たす数値であるとともに、ｋＦが−１≦ｋＦ＜０の負の数値であり、ｋＶが０＜ｋＶの正の数値である。
＜９＞  流速性変動の血流変化成分ＯＦ、ＤＦの各情報に基づき、神経活動に伴う被検体の大脳灰白質層における血流変動を測定する前記＜８＞に記載の生体光測定方法。
＜１０＞  容積性変動の血液変化成分ＯＶ、ＤＶの各情報に基づき、被検体の頭皮層における血流変動を測定する前記＜８＞から＜９＞のいずれかに記載の生体光測定方法。
【発明の効果】
【００１１】
本発明によれば、従来技術における前記諸問題を解決することができ、市販のｆＮＩＲＳ装置を装置改変をすることなく、一般ユーザーが取得できるデータを用い、生体活動における特定信号のみを抽出する生体光測定装置、プログラム及び生体光測定方法を提供することができる。
【図面の簡単な説明】
【００１２】
【図１】ｓｉｎｇｌｅ  ｄｉｓｔａｎｃｅ法の概要を示す説明図である。
【図２】ｍｕｌｔｉ  ｄｉｓｔａｎｃｅ法の概要を示す説明図である。
【図３】生理学的諸属性のもとに生じる血流変化の諸相を示す説明図である。
【図４】プローブ配置の一例を示す説明図である。
【図５】プローブホルダの一例を示す説明図である。
【図６】プローブホルダの一例を示す分解斜視図である。
【図７】被験者１の各実験課題でのデータから求めた相互情報量Ｉ（ｋＦ＊，ｋＶ）のｋＶ依存性を示すグラフである。
【図８】全ての被験者、課題条件のデータから求められたｋＶ＊のヒストグラムを示すグラフである。
【図９】被験者１の上体傾斜課題遂行に伴うヘモグロビン変化に関する従来法と本発明の実施例に係る提案手法との結果比較を示すグラフである。なお、グラフ中の各表記は、次の通りである。Ｃｏｎｖｅｎｔｉｏｎａｌ  Ｏ＆Ｄ：従来法，Ｐｒｏｐｏｓｅｄ  ＯＦ＆ＤＦ：流速性変動成分，Ｐｒｏｐｏｓｅｄ  Ｏｖ＆Ｄｖ：容積性変動成分。灰線：オキシヘモグロビン，黒線：デオキシヘモグロビン。黒太線：課題遂行時間。また、各グラフにおける縦軸の単位は、[ｍＭ×ｃｍ]である。
【図１０】被験者１の呼吸停止課題遂行に伴うヘモグロビン変化に関する従来法と本発明の実施例に係る提案手法との結果比較を示すグラフである。なお、グラフ中の各表記は、次の通りである。Ｃｏｎｖｅｎｔｉｏｎａｌ  Ｏ＆Ｄ：従来法，Ｐｒｏｐｏｓｅｄ  ＯＦ＆ＤＦ：流速性変動成分，Ｐｒｏｐｏｓｅｄ  Ｏｖ＆Ｄｖ：容積性変動成分。灰線：オキシヘモグロビン，黒線：デオキシヘモグロビン。黒太線：課題遂行時間。また、各グラフにおける縦軸の単位は、[ｍＭ×ｃｍ]である。
【図１１】被験者１の指先運動（タッピング）課題遂行に伴うヘモグロビン変化に関する従来法と本発明の実施例に係る提案手法との結果比較を示すグラフである。なお、グラフ中の各表記は、次の通りである。Ｃｏｎｖｅｎｔｉｏｎａｌ  Ｏ＆Ｄ：従来法，Ｐｒｏｐｏｓｅｄ  ＯＦ＆ＤＦ：流速性変動成分，Ｐｒｏｐｏｓｅｄ  Ｏｖ＆Ｄｖ：容積性変動成分。灰線：オキシヘモグロビン，黒線：デオキシヘモグロビン。黒太線：課題遂行時間（０−２０秒：左手，４０−６０秒：右手）。また、各グラフにおける縦軸の単位は、[ｍＭ×ｃｍ]である。
【図１２】被験者１の各課題遂行に伴うヘモグロビン変化に関する本発明の実施例に係る提案手法とｍｕｌｔｉ  ｄｉｓｔａｎｃｅ法との結果比較を示すグラフである。なお、グラフ中の各表記は、次の通りである。Ｂｏｄｙ  ｔｉｌｔｉｎｇ：上体傾斜，Ｂｒｅａｔｈ  ｈｏｌｄｉｎｇ：呼吸停止，Ｆｉｎｇｅｒ  ｔａｐｐｉｎｇ：指先運動。黒線：提案手法の３０ｍｍ配置で得られた流速性変動成分，灰線：ｍｕｌｉｔｉ  ｄｉｓｔａｎｃｅ法で得られた大脳灰白質層での変化。黒太線：課題遂行時間（０−２０秒：左手，４０−６０秒：右手）。また、各グラフにおける縦軸の単位は、[ｍＭ×ｃｍ]である。
【図１３】本発明の実施例に係る提案手法で分離された各変動成分の光減衰への寄与の光源−検出器間距離に対する依存性を示すグラフである。図中、左部分は、容積性変動成分を示し、右部分は、流速性変動成分を示す。なお、グラフ中の各表記は、次の通りである。灰線：各被験者各課題で得られた依存性曲線，黒線：モンテカルロシミュレーションで求めた頭皮層（左）と大脳灰白質層（右）の部分光路長の依存性曲線。
【発明を実施するための形態】
【００１３】
（生体光計測装置及びプログラム）
本発明の生体光計測装置は、少なくとも光照射手段と、光検出部と、演算処理部とを有してなる。
【００１４】
＜光照射手段＞
前記光照射手段は、近赤外光を被検体に照射する手段である。
その具体的な手段としては、特に制限はなく、目的に応じて適宜選択することができ、市販のｆＮＩＲＳ装置で用いられている近赤外光の光源等、公知の光照射手段を全て挙げることができる。
なお、本明細書において、前記ｆＮＩＲＳ（Ｆｕｎｃｔｉｏｎａｌ  Ｎｅａｒ−Ｉｎｆｒａｒｅｄ  Ｓｐｅｃｔｒｏｓｃｏｐｙ）装置とは、被検体に近赤外光を照射して得られる生体信号を測定する装置を指し、機能的近赤外光装置、近赤外光イメージング装置等とも称される装置の全般を指す。
【００１５】
＜光検出部＞
前記光検出部は、前記被検体内を経由した光を検出する。
前記光検出部の具体的なものとしては、特に制限はなく、目的に応じて適宜選択することができ、市販のｆＮＩＲＳ装置で用いられているもの等、公知の装置例を全て挙げることができる。
【００１６】
＜演算処理部＞
前記演算処理部は、前記光から取得される前記被検体における血液中のオキシヘモグロビン変化量Ｏ及びデオキシヘモグロビン変化量Ｄの各情報から、下記式（６）及び（７）で表される、前記オキシヘモグロビン変化量Ｏにおける流速性変動の血流変化成分ＯＦ及び前記デオキシヘモグロビン変化量Ｄにおける流速性変動の血流変化成分ＤＦ、並びに／又は前記オキシヘモグロビン変化量Ｏにおける容積性変動の血液変化成分ＯＶ及び前記デオキシヘモグロビン変化量Ｄにおける容積性変動の血流変化成分ＤＶを算出して出力する。
0005641355.tif 000005
ただし、前記式（６）及び（７）において、ｋＦ及びｋＶは、ｋＦ≠ｋＶを満たす数値である。
【００１７】
前記演算処理部の装置例としては、特に制限はなく、目的に応じて適宜選択することができ、市販のｆＮＩＲＳ装置に配される電子計算機、該ｆＮＩＲＳ装置と外部接続可能な電子計算機などが挙げられる。
【００１８】
−プログラム−
本発明のプログラムは、本発明の前記生体光計測装置と協動して前記演算処理を実行させるプログラムであり、具体的には、前記生体光計測装置の前記演算処理部に対して、前記光から取得される前記被検体における血液中のオキシヘモグロビン変化量Ｏ及びデオキシヘモグロビン変化量Ｄの各情報から、前記式（６）及び（７）で表される、前記オキシヘモグロビン変化量Ｏにおける流速性変動の血流変化成分ＯＦ及び前記デオキシヘモグロビン変化量Ｄにおける流速性変動の血流変化成分ＤＦ、並びに／又は前記オキシヘモグロビン変化量Ｏにおける容積性変動の血液変化成分ＯＶ及び前記デオキシヘモグロビン変化量Ｄにおける容積性変動の血流変化成分ＤＶを算出して出力させることを特徴とする。
【００１９】
上述した演算処理の具体的な意義について、前記被検体の血流変動の態様を解き明かしつつ、以下に説明する。
【００２０】
［血流変動の二つの様態（モダリティ）］
−概念−
前記ｆＮＩＲＳ装置による計測では、前記被検体の生体組織の活動に伴うオキシヘモグロビン量、デオキシヘモグロビン量の変化（Ｏ，Ｄ）をそれぞれ独立の二変量として推定できる。どのようなＯ，Ｄの時間変化が脳活動に伴うｈｅｍｏｄｙｎａｍｉｃｓとして典型的なものであるか、については今までに多くの報告がなされてきた。
初期の研究では、オキシの卓越した増大を報告するものが多いが（非特許文献７，８参照）、これらは光路上に存在する脳以外の生体組織での血流変化の可能性に配慮していない。
一方で、測定誤差、体動、全身性生理活動の影響等に十分な配慮を払った比較的近年の研究ではタスクと同期してオキシの増大と同時にデオキシの減少が生じることが報告されている（非特許文献３，９参照）。ラットを用いた開頭下での脳機能信号の計測結果において、同様のオキシ、デオキシの相反的な変化、即ち負の相関が認められている（非特許文献１０及び１１参照）。
【００２１】
ＰＥＴや造影剤を用いたＭＲＩなどの研究からは、脳活動に伴う局所的な血流変動では、ｒＣＢＦが数１０％ほども増大するのに対して、ｒＣＢＶは高々１０％に満たない程度しか増大しないことが比較的古くから知られている（非特許文献１２，１３，１４参照）。
ＣＢＶ増大の実体として、ラットの細動脈が神経活動に伴って血管径拡張を生じることが顕微鏡下で確認されているものの（非特許文献１５参照）、こうした血管領域での酸素交換効率は、毛細血管と比較して無視できるほど小さい（非特許文献１６参照）。
他方、酸素供給の場となる毛細血管は、多くが静止血液中の赤血球半径よりも小さい管径を有し、赤血球は、血漿流の流速勾配で生じるずり応力で変形することによってこれらの毛細血管を通過することが知られている（非特許文献１７参照）。この機構上、毛細血管では、血流量の増減によって管径がほとんど変動しないと考えられる。
これら事項は、脳活動に伴う酸素交換に関わるｈｅｍｏｄｙｎａｍｉｃｓでは、血液量の増大が決して大きくないことを示唆している。
【００２２】
いま仮に毛細血管では、血管径がほぼ一定であるとすると、脳神経活動に随伴して血管の上流側から高い酸素飽和度をもった動脈血の流入が増大した場合、局所的な酸素消費とバランスして低い酸素飽和度を実現していたデオキシ−ヘモグロビンが押し流されて、オキシ−ヘモグロビンに置き換えられるので、血管中では、オキシ増大とデオキシ減少が相補的に生じるはずである。
即ち、そこでは、Ｄ１＝ｋ１Ｏ１なる関係が成り立ち、ｋ１は、−符号の定係数として記述できる。係数ｋ１は、流速に依存したヘマトクリット値の変化やわずかな受動的血管拡張の有無により異なるであろうが、最小で−１、大きくとも０の間の負の値をとると期待される。
本明細書では、この血流変化の様態（モダリティ）を流速性血流変化（ｆｌｏｗ  ｖｅｌｏｃｉｔｙ−ｓｅｎｓｉｔｉｖｅ  ｃｏｍｐｏｎｅｎｔ）と呼ぶ。
【００２３】
他方、一部の研究者は、体動や全身性生理活動に由来する血流変動では、オキシ、デオキシが同一符号の向きに増大・減少するケースを経験的に知っており、バルサルバ法（非特許文献１８参照）やおじぎ、息止め（非特許文献３参照）では、オキシ、デオキシが同相に変化することを報告している。こうした全身性の血流変化には、自律神経を介した動脈−細動脈の血管収縮・拡張や血圧変化・姿勢変化に伴う鬱血など動静脈の受動的容積変動などが伴うと考えられる。このような血管容積の変化に伴い、オキシ量、デオキシ量は、各々の存在比率にしたがって互いに同相で変化する。
即ち、そこでもやはりＤ２＝ｋ２Ｏ２なる関係が成り立つが、ｋ２は＋符号の定係数となる。係数ｋ２は、当該血管中の血液の酸素飽和度によって決まるため、酸素飽和度１００％の場合に０、酸素飽和度０％の場合に＋∞をとり、通常は、その間の正の値をとると期待される。
これらの太い血管領域では、ほとんど組織への酸素供給は行われないため、血液の酸素飽和度は、ほぼ一定と考えてよい。そのため、これらの血管では血流速度が変化したとしても上で述べた流速性血流変化は生じない。
そこで、本明細書では、この血流変化の様態（モダリティ）を前記流速性血流変化と独立した成分と考え、容積性血流変化（ｖｏｌｕｍｅ−ｓｅｎｓｉｔｉｖｅ  ｃｏｍｐｏｎｅｎｔ）と呼ぶ。
【００２４】
以上の血流モダリティと生理学的諸属性を図３に示す。各モダリティを特徴付ける比例係数ｋ１，ｋ２の各値は、ヘマトクリット値、血中酸素飽和度、血管壁の弾性率などによって定まる。これらは生体恒常性を伴うパラメータであるから計測時間内では、ｋ１，ｋ２を十分に一定と見なせる。
したがって、従来のｓｉｎｇｌｅ  ｄｉｓｔａｎｃｅ法によるｆＮＩＲＳ計測信号を、このように生成機序が明確に異なる二つの血流モダリティの混在として考えることが可能となる。
このような解釈に立ち、前記ｓｉｎｇｌｅ  ｄｉｓｔａｎｃｅ法によるｆＮＩＲＳ計測信号から各モダリティ成分の時間変化を各々分離する手法について以下に詳述する。
【００２５】
−モデリング−
前記ｓｉｎｇｌｅ  ｄｉｓｔａｎｃｅ法でｆＮＩＲＳ計測される情報は、オキシ−ヘモグロビン変化量及びデオキシ−ヘモグロビン変化量という二変量である。以下これをＯ，Ｄと記す。これらは、流速性変動の血流変化成分ＯＦ，ＤＦと容積性変動の血流変化成分ＯＶ，ＤＶによって構成されていると考える。即ち、
0005641355.tif 000006
【００２６】
ここで、ｋＦは、−１≦ｋＦ＜０の負の定数、ｋＶは、０＜ｋＶの正の定数である。このとき前記式（１）は、以下の形に変形できる。
0005641355.tif 000007
【００２７】
ｋＦ≠ｋＶであるので、前記式（３）の右辺の行列は逆行列をもつ。したがって、以下の二つの式が得られる。
0005641355.tif 000008
【００２８】
それぞれの血流変化の成分に関して式をまとめると以下となる。
0005641355.tif 000009
【００２９】
以上に述べたように、前記式（１），（２）で記述される二つの血流変化モダリティのモデルに基づくことにより、ｓｉｎｇｌｅ  ｄｉｓｔａｎｃｅ法のｆＮＩＲＳ計測で得られる二変量Ｏ，Ｄに含まれる容積性血流変化と流速性血流変化を前記式（６），（７）として分離することができる。
【００３０】
−各モダリティ係数ｋＶ，ｋＦの決定−
上述の成分分離に基づく演算処理を有効に実行するためには、ｋＶ，ｋＦの具体的な数値を決定しなければならない。
この数値を決定する手法としては、特に制限はなく、目的に応じて適宜選択することができるが、ここでは、より適切な数値を決定する手法として、算術的アプローチと、実証的アプローチの二つの手法を説明する。
【００３１】
−算術的アプローチ−
前記算術的アプローチは、算術的手法を用いて数値を決定する手法である。
前記算術的手法としては、前記ｋＶ，ｋＦの具体的な数値を決定できる手法であれば、特に制限はなく、目的に応じて適宜選択することができ、例えば、容積性血流変化と流速性血流変化の独立性を判別する関数を用いた手法などの算術的手法が挙げられる。
【００３２】
上述の容積性血流変化と流速性血流変化の独立性を判別する関数としては、特に制限はなく、目的に応じて適宜選択することができ、例えば、（１）相互情報量（Ｍｕｔｕｒａｌ  ｉｎｆｏｒｍａｔｉｏｎ）又は伝達情報量（Ｔｒａｎｓｉｎｆｏｒｍａｔｉｏｎ）と呼ばれる量を定義する関数、（２）相関係数（ｃｏｒｒｅｌａｔｉｏｎ  ｃｏｅｆｆｉｃｉｅｎｔ）又はピアソンの積率相関係数（Ｐｅａｒｓｏｎ  ｐｒｏｄｕｃｔ−ｍｏｍｅｎｔ  ｃｏｒｒｅｌａｔｉｏｎ  ｃｏｅｆｆｉｃｉｅｎｔ）と呼ばれる量を定義する関数が挙げられる。
【００３３】
（１）相互情報量又は伝達情報量と呼ばれる量を定義する関数
この量は、二つの変数ＯＦ（ｋＦ），ＯＶ（ｋＶ）が離散値をとる場合、以下のように与えられる。
0005641355.tif 000010
ただし、前記式（８）において、ｐ（ＯＦ，ＯＶ）は、ＯＦ，ＯＶの同時分布関数であり、ｐ（ＯＦ）とｐ（ＯＶ）は、それぞれＯＦとＯＶの周辺確率分布関数である。
【００３４】
また、前記量は、ＯＦ（ｋＦ），ＯＶ（ｋＶ）が連続値をとる場合、以下のように与えられる。
0005641355.tif 000011
ただし、前記式（９）において、ｐ（ＯＦ，ＯＶ）は、ＯＦ，ＯＶの同時分布密度関数であり、ｐ（ＯＦ）とｐ（ＯＶ）は、それぞれＯＦとＯＶの周辺確率密度関数である。
【００３５】
前記式（８）と（９）のいずれの場合も、ｐ（ＯＦ，ＯＶ），ｐ（ＯＦ）及びｐ（ＯＶ）は、ＯＦ（ｋＦ）とＯＶ（ｋＶ）の頻度分布に基づいて求めることができる。Ｉ（ｋＦ，ｋＶ）は、常に正の値をとり、二つの変数ＯＦ（ｋＦ），ＯＶ（ｋＶ）の取る値の独立性が高いほど小さな値をとる。この性質を用いてＩ（ｋＦ，ｋＶ）が最小になる（ｋＦ，ｋＶ）の組み合わせを求めることができる。
【００３６】
（２）相関係数又はピアソンの積率相関係数と呼ばれる量を定義する関数
前記式（６），（７）を用いる提案手法では、各時刻ｉごとに変数ＯＦ，ｉ（ｋＦ）とＯＶ，ｉ（ｋＶ）の組がｎ個与えられる。この場合、この量は以下のように与えられる。
0005641355.tif 000012
【００３７】
前記式（１０）において、ρ（ｋＦ，ｋＶ）は、−１から１までの値を取り、二つの変数ＯＦ（ｋＦ），ＯＶ（ｋＶ）の相関が低いほどゼロに近づく。この性質を用いてρ（ｋＦ，ｋＶ）の絶対値が最小になる（ｋＦ，ｋＶ）の組み合わせを求めることができる。
【００３８】
−実証的アプローチ−
実証的アプローチは、生体メカニズムの考察及び過去の研究結果から、妥当な数値を決定する手法である。以下、詳細に説明する。
【００３９】
先ず、ｋＦの数値について検討する。
この血流成分が生じ得る毛細血管では、その直径が時間的に変化しないと想定できる。脳神経活動に随伴して血管の上流側から高い酸素飽和度をもった動脈血の流入が増大した場合、デオキシ−ヘモグロビンは押し流されてオキシ−ヘモグロビンに置き換えられる。この効果のみを考えた場合ｋＦ＝−１となる。
しかし、実際には血流速度の増大に伴ってヘマトクリット値（血液中に占める血球の容積の割合）が上昇したり、受動的な血管拡張が生じる可能性がある。
ただし、それらは血流速度の増大に対する応答作用であるから、血流速度増大のみの効果を完全に打ち消すには至らない。即ち、ｋＦ＜０が保証される。
以上より−１≦ｋＦ＜０を論理的に妥当と考えることができる。
【００４０】
過去１０数年間に行われた内外の研究のうち、正しいｋＦの値を推測し得る条件（即ち、視覚刺激実験、触覚刺激実験、適切な対照実験を伴う運動課題実験、大脳灰白質の血流変化のみを抽出する測定法を用いた任意の実験、頭蓋除去状態で行った任意の実験）で得られたデータに基づいてｋＦの値を算出したところ、その値は、下記表１に示すとおり、−０．８５＜ｋＦ＜−０．３５の範囲にあった。
したがって、過去の研究例からは、−０．８５＜ｋＦ＜−０．３５が見出される。
【００４１】
【表１】
0005641355.tif 000013
【００４２】
脳という重要な器官とそこへの酸素・栄養の供給機構を動物が進化の途上で獲得してきたことを考えると、異なる種にわたり、脳の領野によらず共通の機構、即ち、共通のｋＦ値を広く見出せることが期待できる。上記表１で見る限り、−０．６の±１０％以内の範囲にｋＦの値を示した事例は全体の半数を超えている。
したがって、この共通のｋＦの値として、約−０．６（−０．６±１０％程度）を採用することが好ましいと考えられる。
【００４３】
次に、ｋＶについて検討する。
この血流成分が生じ得る（毛細血管以外の）動脈、細動脈、細静脈、静脈からは組織に直接酸素は分配されないので、これらの血管中の血液のオキシ−，デオキシ−ヘモグロビンの割合はそれぞれの一定値を保つと想定できる。その条件下で血液体積が増大すれば、オキシ−，デオキシ−ヘモグロビンはそれぞれが占める割合に応じて増大する。いま血中の総ヘモグロビンに占めるオキシヘモグロビンの割合をＳｘＯ２とすると、その関係は、以下で表すことができる。
0005641355.tif 000014
前記式（１１）において、ＳｘＯ２として論理上考え得る最大値は、血管内が全てオキシヘモグロビンで満たされた場合の１００％であり、最小値は、逆に全てデオキシヘモグロビンで満たされた場合の０％である。これらに対応するｋＶは、前記式（１１）から、それぞれ場合で０と＋∞である。
したがって、論理的には０＜ｋＶ＜＋∞の範囲を想定し得る。
【００４４】
実際には、生理学的に動脈血のＳｘＯ２は、１００％に非常に近く、９５％以下には容易にならないことが知られている。また静脈血のＳｘＯ２は、安静時には約７０％を示し（非特許文献３０参照）、激しい運動が持続した場合にも５０％程度以下には低下しないことが知られている。
もしも、この血流成分が動脈血のみで生じ、そのＳｘＯ２が１００％だった場合、ｋＶは最も低い値０をとる。一方、血流成分が静脈血のみで生じ、そのＳｘＯ２が生理的に考え得る最も低い５０％程度だったとするとｋＶは最も高い値１をとる。
以上の考察から、生理学的には、０＜ｋＶ≦１の範囲を示すのが妥当と考えられる。
【００４５】
実際には、この血流成分には動脈血と静脈血の両方が寄与する。激しい運動負荷をかけない通常の計測条件では、動脈血のＳｘＯ２が１００％、静脈血のＳｘＯ２が約７０％程度である。前記式（１１）から、それぞれに対応するｋＶ値は、０と０．４３である。
したがって、準安静的な条件下で行われる計測では、０＜ｋＶ＜０．５程度の範囲の値を取ると考えられる。
【００４６】
なお、前記ｋＶ，ｋＦの決定に際し、前記算術的アプローチと前記実証的アプローチは、それぞれ独立して考えることができる。ただし、これらを併用して前記ｋＶ，ｋＦを決定してもよい。
例えば、前記算術的アプローチにおいて、前記式（８）〜（１０）の演算処理を実施した結果、Ｉ（ｋＦ，ｋＶ），ρ（ｋＦ，ｋＶ）の最小点、即ち、測定波形の極小点が定かに確認されない場合には、前記ｋＶ，ｋＦのいずれかを実証的アプローチにより決定し、前記極小点が確認されるように処理してもよい。
【００４７】
＜装置例＞
前記生体光計測装置の装置例について説明する。
上述の通り、前記生体光計測装置は、前記光照射部及び前記光検出部を有する。これらは、被検体を介して光学的に接続されるが、通常、該光照射部から照射される光を被検体に導く照射プローブと、前記被検体内を経由した光を前記光検出部に導く受光プローブとを有して構成される（図１参照）。本発明は、ｍｕｌｔｉ  ｄｉｓｔａｎｃｅ法（図２参照）における参照用プローブを有する構成を排除するものではないが、前記演算処理部にて血流変動の二つの様態（モダリティ）を解析する演算処理を実施することにより、こうした参照プローブがなくても、脳機能由来の信号のみを抽出して測定を行うことができる。
【００４８】
前記照射プローブ及び前記受光プローブを有する構成に関し、これらは、前記生体光計測装置において複数有することとして構成されてもよい。
この場合、前記光検出部は、前記複数の受光プローブから導かれる光をそれぞれ検出するとともに、前記演算処理部が検出された前記各光から取得される情報をそれぞれ演算処理して出力するように構成される。
【００４９】
前記生体光計測装置としては、被検体の頭部における生体信号を測定することに適するが、例えば、前記照射プローブは、前記頭部の右脳部分と左脳部分とに配されてよい。また、前記受光プローブは、前記照射プローブからの距離を一定にして、その周囲に複数配されてよい。この様子を図４に示す。該図中、各プローブの配置箇所としては、符号１１が照射プローブを示し、符号１２が受光プローブを示している。
このようなプローブ配置により、脳の各機能領野における諸神経活動を効果的に測定することができる。
【００５０】
また、前記生体光計測装置としては、前記被検体の表面に沿って取付け可能なプローブホルダを有し、前記プローブホルダに前記照射プローブと前記受光プローブとが保持されることとして構成されてもよい。
このようなプローブホルダを有することにより、前記被検体に対して安定した状態で各プローブを配置することができる。
【００５１】
前記プローブホルダの具体的な構成としては、特に制限はなく、目的に応じて適宜選択することができ、例えば、市販のｆＮＩＲＳ装置に用いられているものを適用することができるが、個人差により様々な曲面を有する被検体の頭部に対し、前記各プローブを垂直に密着保持させ、前記各プローブの配置を一定に保持させる観点から、例えば、図５，図６に示すプローブホルダ２０のように構成してもよい。
【００５２】
プローブホルダ２０は、図５に示すように、少なくとも照射プローブ１１と受光プローブ１２とを保持する二つのプローブソケット１３と、これらのソケットを連結する長板上の上部連結子１４及び下部連結子１５とを有し、該下部連結子１５上に、被検体の表面にフィットするように固定させるとともに外光を遮断する布状部材１９を有する。
【００５３】
プローブソケット１３は、図６に示すように、上部連結子１４及び下部連結子１５の孔に嵌装させながら、ソケット本体１６に対し、下部取付部材１７及び上部取付部材１８のねじ部を螺合して配される。
ここで、上部連結子１４に配される孔は、下部連結子１５に配される丸孔と異なり、上部連結子の長さ方向に拡開された長孔として形成されており、プローブソケット１３は、下部連結子１５における下部取付部材１７の取付け位置を固定しつつ、上部取付部材１８の取付け位置を上部連結子１４の長さ方向で任意に調節することができる。
したがって、図５に示すプローブホルダ２０を被検体頭部の表面形状に沿って湾曲させて取付ける際、上部取付部材１８のを下部連結子１５の孔に対応する位置よりも上部連結子１４の長手方向外側の位置に調節して取付けることにより、被検体頭部に対して、垂直に密着保持されるように、各プローブ１１，１２を配することができる。
【００５４】
（生体光計測方法）
本発明の生体光計測方法は、少なくとも、光照射ステップと、光検出ステップと、演算処理ステップとを有する。
【００５５】
＜光照射ステップ＞
前記光照射ステップでは、近赤外光を被検体に照射する。
前記光照射ステップを実行する具体的な手段としては、特に制限はなく、目的に応じて適宜選択することができ、例えば、本発明の前記生体光計測装置における光照射手段が挙げられる。
【００５６】
＜光検出ステップ＞
前記光検出ステップでは、前記被検体に照射され、前記被検体内を経由した光を検出する。
前記光検出ステップを実行する具体的な手段としては、特に制限はなく、目的に応じて適宜選択することができ、例えば、本発明の前記生体光計測装置における光検出部が挙げられる。
【００５７】
＜演算処理ステップ＞
前記演算処理ステップでは、前記光検出ステップで検出した前記光から取得される情報を演算処理して出力する。
ここで、前記演算処理ステップは、前記光から取得される前記被検体における血液中のオキシヘモグロビン変化量Ｏ及びデオキシヘモグロビン変化量Ｄの各情報から、下記式（６）及び（７）で表される、前記オキシヘモグロビン変化量Ｏにおける流速性変動の血流変化成分ＯＦ及び前記デオキシヘモグロビン変化量Ｄにおける流速性変動の血流変化成分ＤＦ、並びに／又は前記オキシヘモグロビン変化量Ｏにおける容積性変動の血液変化成分ＯＶ及び前記デオキシヘモグロビン変化量Ｄにおける容積性変動の血流変化成分ＤＶを算出して出力することを特徴とする。
0005641355.tif 000015
ただし、前記式（６）及び（７）において、ｋＦ及びｋＶは、ｋＦ≠ｋＶを満たす数値である。
前記演算処理ステップを実行する具体的な手段としては、特に制限はなく、目的に応じて適宜選択することができ、例えば、本発明の前記生体光計測装置における演算処理部が挙げられる。
【実施例】
【００５８】
＜実験内容＞
図５及び図６に示す構成に特徴付けられる、独自に開発したプローブ及びホルダシステム（特許文献２及び非特許文献３１参照）を近赤外光イメージング装置（島津製作所製、ＯＭＭ−３０００）に接続し、検出される光信号に対して、前記式（６），（７）に基づく演算処理を行うようにした。使用する波長は、７７６ｎｍ、８０８ｎｍ、８２８ｎｍとした。サンプリングレート４０Ｈｚで得られた各波長の吸光度データを１０Ｈｚにダウンサンプリングした後、１Ｈｚのローパスフィルタに通した。
被験者の頭部へのプローブの設置は、非特許文献３に記載の設置方法に準拠し、同一のｔａｐｐｉｎｇタスクを用いたＭＲＩ計測で一次運動野を特定された被験者に対して、左右運動野活動部位の直上位置にそれぞれ１つの光源と４つの検出プローブを設置した。光源−検出器は直線上に置かれ、４つの検出器の光源からの距離は、それぞれ１０ｍｍ，２０ｍｍ，３０ｍｍ，４０ｍｍとした。
【００５９】
被験者は、以下のように選定した。即ち、産総研人間工学実験委員会に承認された実験計画に基づき、実験責任者から個別に実験について口頭・文書で説明を受け、理解した健常成人被験者複数名の協力を得た。
【００６０】
以上の実験構成をもとに、先ず、ｎｏｎ−ｆｕｎｃｔｉｏｎａｌ  ｔａｓｋの課題を被験者に実施させた。
被験者を椅子に座らせ、初期安静２０秒に続き、課題２０秒と安静２０秒を５回繰り返し実施させた。合計時間は２２０秒である。課題内容は、上体傾斜と息止めとした。
上体傾斜の課題では、視覚及び聴覚指示を合図に、上体のみを前傾させ、傾斜角度３０°になる位置に置かれたバーに鼻先を接触させて、その姿勢を保持させた。２０秒後の次の視覚及び聴覚指示を合図に初期安静時の上体正立に戻ってその姿勢を保持させた。
また、息止めの課題では、同様の指示を合図に呼吸を停止させ、そのまま２０秒間保持させた。２０秒後の指示を合図に呼吸を再開させた。
【００６１】
次に、ｆｕｎｃｔｉｏｎａｌ  ｔａｓｋの課題を被験者に実施させた。
ｎｏｎ−ｆｕｎｃｔｉｏｎａｌ  ｔａｓｋ時と同様に被験者を椅子に座らせ、上体正立姿勢を保持させた。この状態で、初期安静２０秒に続き、左手タッピング動作２０秒、安静２０秒、右手タッピング動作２０秒、安静２０秒を５回繰り返して実施させた。合計時間は４２０秒である。タッピング動作は、各々の片側の手の親指と人差し指の対向動作を４Ｈｚで行う条件の下、実施させた。
【００６２】
−ｋＦの決定−
上述の実証的アプローチに基づき、ここでは、ｋＦの値としてｋＦ＊＝０．６とした。
【００６３】
−ｋＶの収束性−
また、ｋＶに関し、上述の算術的アプローチにおける前記式（８）の相互情報量を定義する関数に基づき、その情報独立性の仮定からＩ（ｋＦ＊，ｋＶ）が最小となるようにｋＶ＊を定めた。これにより、それぞれの計測データごとに（ＯＦ，ＤＦ）と（ＯＶ，ＤＶ）とを一意に分離した。
図７に、異なる光源−検出器間距離、タスクの種類の条件に応じて、二つの血流モダリティ間の相互情報量Ｉ（ｋＦ＊，ｋＶ）のｋＶ依存性がどのように異なるかについての、被験者１のデータを示す。ほとんど全ての条件下で、Ｉ（ｋＦ＊，ｋＶ）は極小点をもち、そのときのｋＶの値ｋＶ＊は、本発明の提案モデルにとって妥当と云える範囲（０＜ｋＶ＜１）にあることが分かる。
【００６４】
次いで図８に、異なる被験者のデータを用いて得られた光源−検出器間距離が３０ｍｍにおけるｋＶ＊のヒストグラムを示す。このヒストグラムによれば、ｋＶ＊のほとんどは、０＜ｋＶ＊＜０．４の範囲に分布している。これは、今回実施したタスクが大きな運動負荷を伴わないために、静脈血の酸素飽和度が安静状態の値７０％前後であった、と考えれば、ｋＶ＊の値のほとんどが０．４以下に収まる結果となったことを理解できる。
また、ゼロに近い値を持つものがより多いことは、容積性の血流変動全体の中で、動脈血管での容積変動が寄与する割合が相対的に大きいことを示唆している。
以上の結果から、二つの血流モダリティを想定するモデルがこれらの課題、測定条件下での血流変化モデルとして概ね妥当であり、ｋＦ＊やｋＶ＊は妥当に推定されていると考えられる。
【００６５】
−課題遂行時の時間変化−
図９に、上体前傾の課題遂行に伴う従来法（前記実験内容の条件でヘモグロビン変化Ｏ，Ｄを分離しない手法）でのヘモグロビン変化Ｏ，Ｄ、並びに、本発明の実施例に係る提案手法の容積性変動成分ＯＶ，ＤＶ及び流速性変動成分ＯＦ，ＤＦの被験者１の結果を示す。すでに、我々が報告しているように（非特許文献３参照）、従来法のＯ，Ｄは、上体を傾斜することにより大きく変化する。
この変化については、光源−検出器間距離が異なる様々な条件でも同様の形を示していることが分かる。
また、その振幅は光源−検出器間距離が長いほど大きいが、これはＯ，Ｄの計算において光路長による補正を行っていないためと考えられる。提案手法の容積性変動成分の変化はどのチャンネルにおいても従来法の時間変化と非常によく似た挙動を示していることが分かる。
一方で、流速性変動成分は、他の二つに較べてその変化の振幅は非常に小さく、上体傾斜の影響をほとんど受けていないことが分かる。
【００６６】
次いで図１０に、呼吸停止の課題遂行時の各ヘモグロビン変化の被験者１での結果を示す。やはりここでも、我々が以前に報告した（非特許文献３参照）呼吸停止に伴うアーティファクトが従来法のヘモグロビン変化に確認された。
光源−検出器間距離の異なる条件でもその挙動は非常に類似している。
提案手法の容積性変動成分は、ここでも従来法の各条件のヘモグロビン変化に非常によく似た挙動を示した。
一方で、流速性変動成分は他の二つに較べて変化の振幅がやはり小さく、呼吸停止の影響をほとんど受けていないことが分かる。
【００６７】
そして図１１に、手指タッピング課題遂行に伴う各ヘモグロビン変化の被験者１での結果を示す。従来法で得られたヘモグロビン変化は、タッピングを行った指の左右の違いに対する明確な片側性は認められない。
同時に、光源−検出器間距離の異なる条件間の比較では、上体傾斜や呼吸停止のとき程には時間変化の形の類似は認められない。
提案手法の容積性変動成分は、数１０秒のタイムスパンで生じる緩やかなベースラインの変動と左右によらない課題遂行と同期した上昇が見られる。この形は光源−検出器間距離の異なる各条件間でよく類似している。
一方、流速性変動成分は、光源−検出器間距離３０ｍｍ，４０ｍｍの条件で、左タッピングに対して右側で、右タッピングに対して左側で明確なＯＦの上昇（ＤＦの減少）が見られる。なお、この変化は、光源−検出器間距離を２０ｍｍ，１０ｍｍとした条件ではほとんど認められない。
【００６８】
−流速性変動成分とｍｕｌｔｉ  ｄｉｓｔａｎｃｅ法による脳機能信号の比較−
測定データのうち、光源−検出器間距離を２０ｍｍ，３０ｍｍとした条件で測定されたものを用いることにより、我々が提案したｍｕｌｔｉ  ｄｉｓｔａｎｃｅ法を用いて大脳灰白質でのヘモグロビン変化を算出することができる。
図１２に、この方法で得た大脳灰白質におけるヘモグロビン変化と光源−検出器間距離３０ｍｍで得られた流速性変動成分ＯＦ，ＤＦの比較を示す。
このｍｕｌｔｉ  ｄｉｓｔａｎｃｅ法と、本発明の実施例に係る提案手法の二つの方法で得たヘモグロビン変化は振幅のスケーリングが異なる。このため図１２では、ここに掲げた条件での両者の時間変動が総体として最も一致するようにｍｕｌｔｉ  ｄｉｓｔａｎｃｅ法で得たヘモグロビン変化のスケーリングを調節して表示している。
上体傾斜及びタッピングの課題で、それぞれの手法で得られたヘモグロビン変化は、非常によく一致していることが分かる。
一方、呼吸停止の課題遂行での時間変化は、それほどよく一致せず、提案手法の流速性変動成分の振幅は、ｍｕｌｔｉ  ｄｉｓｔａｃｅ法のそれに較べて大きい。この意味については後述する。
以上のように、提案手法を用いて実際に従来のプローブ配置によるｆＮＩＲＳ計測データから容積性変動と流速性変動の二つの成分を分離することができ、そのうちの流速性変動成分は、ｍｕｌｔｉ  ｄｉｓｔａｎｃｅ法で得られる大脳灰白質の血流変化ときわめて類似した変化を示すことが明らかとなった。
【００６９】
＜考察＞
二つの血流モダリティは、少なくともモデル上ではオキシとデオキシの相関関係の違いとして定義されているに過ぎない。我々は今回、それらと脳機能信号またはその他の信号との対応関係を検証することができる。
図９，図１０では、同形の容積性変動の振幅が光源−検出器間距離の増大にしたがって大きくなる傾向が明瞭に確認された。
また、図１１のタッピング実験でも、容積性変動については光源−検出器間距離の異なるチャンネル間で比較的類似した挙動が認められる。一方で、流速性変動の場合、これとは異なり、光源−検出器間距離が３０ｍｍを超える場合でのみ認められる。
【００７０】
このように各血流モダリティ成分の振幅は、光源−検出器間距離依存性が異なっている。この点をより詳細に議論するために、我々はタッピング実験で得られた容積性成分と流速性成分の時間変動の振幅の比率の各光源−検出器間距離による違いを各々の被験者ごとに定量した。結果を図１３に示す。なお、容積性成分と流速性成分の振幅の大きさは、当該ヘモグロビン濃度変化と部分光路長の積で決定される。
近赤外光のヒト頭部組織での伝播はシミュレーションでも研究されており、各組織層での部分光路長もまた光源−検出器間距離によって変化し、その依存性は組織層ごとに異なることが知られている（非特許文献３２参照）。
もしも、我々が仮定した通りに流速性成分と容積性成分の生理学的起源が異なるとすると、流速性成分は、実際上では大脳灰白質のみで生じ、容積性成分の大部分は、頭皮層で生じると考えられる。
その場合、図１３に示した各成分の振幅の光源−検出器間距離依存性は、その成分が生じる組織層の部分光路長の光源−検出器間距離依存性に類似することが期待される。
【００７１】
そこで我々は、観測波長の一つである７７６ｎｍでの頭皮層と大脳灰白質層での部分光路長の光源−検出器間距離依存性を積層モデルを用いてシミュレーションした。その結果を図１３に示す。なお、このシミュレーションには、マサチューセッツ総合病院が公開している計算コードを用い（ｈｔｔｐ：／／ｗｗｗ．ｎｍｒ．ｍｇｈ．ｈａｒｖａｒｄ．ｅｄｕ／ＰＭＩ／ｉｎｄｅｘ．ｈｔｍｌ）、生体組織の光学定数は文献値（非特許文献３３参照）を用いた。
結果として、頭皮層と灰白質層の部分光路長の光源−検出器間距離依存性は明確に異なり、前者は容積性変動成分と、また後者は流速性変動成分と非常に類似していることが分かる。
更に、図７で示したｍｕｌｔｉ  ｄｉｓｔａｎｃｅ法による灰白質層のヘモグロビン変化と流速性変動成分が極めてよく類似していることを考えあわせると、容積性変動成分は、頭皮層で、流速性変動成分は灰白質層でそれぞれ主に生じていると考えることによって、以上を整合的に理解できる。
以上の考察により、提案手法によって分離された流速性変動成分は灰白質層における神経活動に伴う血流変動を相当程度反映していると考えられる。
【００７２】
本実験では、流速性変動を特徴付けるｋＦ＊を文献値から０．６と決定したが、種による違いや領野間での違いなどはそこでは考慮しなかった。それは、脳組織でこの血流変動を生じさせる機構がほ乳類で比較的よく保存されていると考えたためである。
しかし、幾つかの病理的なケース、例えば、脳虚血等では、ｋＦ＊の数値が異なる可能性も考えられる。そうした病理例をこの手法で計測する場合には、収束性を犠牲にしてもｋＦとｋＶの両方を自由変数として収束解を求めるようにアルゴリズムを変更して用いることも一策と考えられる。
【００７３】
呼吸停止の課題遂行時（図１０参照）の流速性変動成分を厳密に観察すると、徐々にオキシが減少し、デオキシが増大していくわずかな傾向が多くの被験者に認められた。
これは酸素供給の停止に伴う血中酸素飽和度のわずかな減少を反映したものと考えられるが、この流速性変動は概して小さいため、今回のような負荷の低い課題では深刻ではなかった。ただし、大きな運動負荷が伴う課題下での計測では、この点も考慮して、ｋＦとｋＶを決定することも一策と考えられる。
【００７４】
ｆＮＩＲＳ計測信号のオキシとデオキシの逆相関性がアーティファクトの混入によって低下するという仮定のもとに、主としてプローブのミスコンタクトによるアーティファクトの除去を提案した報告がある（非特許文献３４参照）。
プローブのミスコンタクトにおける見かけ上のヘモグロビン変化は、入射・受光効率の変動が原因であるので、用いる波長に応じた装置特有の性質を示す。
その変動は、上述した容積性変動と同じようにオキシ，デオキシ間で正の相関を持ち、ｋＶ＊に相当する比例係数は装置ごとに異なるが、概ね０．５以上の値を取る。
我々が行った実験は、独自のプローブホルダシステムを用いているためプローブのミスコンタクトの問題はほとんど生じないと考えられる。このことは、ｋＶ＊のヒストグラム（図８参照）において、ｋＶ＊が０．５以上の値を取る頻度が極端に小さいことからも示されている。ただし、このように十分な実験上の配慮を行わない場合のｋＶ＊の収束性やその結果である脳機能信号抽出のパフォーマンスに関しては、別途検討を加えることも一策である。
【００７５】
以上に述べたように、標準的なｆＮＩＲＳ計測条件下での計測では、提案手法で抽出した流速性変動成分は大脳灰白質層における神経活動に伴う血流変化を反映するものと考えられる。この手法を用いることにより、すでに従来型のｆＮＩＲＳ装置をもつユーザは新たにハードウェアを更新する必要なく、より信頼性の高いｆＮＩＲＳ計測を行うことが可能となる。
【符号の説明】
【００７６】
１，１１    照射プローブ
２，１２    受光プローブ
３    頭皮
４    頭蓋骨
５    脳脊髄液
６    大脳灰白質層
７    大脳白質層
８    参照用プローブ
１０    頭部
１３    プローブソケット
１４    上部連結子
１５    下部連結子
１６    ソケット本体
１７    下部取付部材
１８    上部取付部材
１９    布状部材
２０    プローブホルダ
【先行技術文献】
【特許文献】
【００７７】
【特許文献１】特開２００９−１３６４３４号公報
【特許文献２】特開２００９−１７８１９２号公報
【非特許文献】
【００７８】
【非特許文献１】Ｈｕｐｐｅｒｔ，Ｔ．Ｊ．，Ｄｉａｍｏｎｄ，Ｓ．Ｇ．，Ｆｒａｎｃｅｓｃｈｉｎｉ，Ｍ．Ａ．，Ｂｏａｓ，Ｄ．Ａ．，２００９．ＨｏｍＥＲ：ａ  ｒｅｖｉｅｗ  ｏｆ  ｔｉｍｅ−ｓｅｒｉｅｓ  ａｎａｌｙｓｉｓ  ｍｅｔｈｏｄｓ  ｆｏｒ  ｎｅａｒ−ｉｎｆｒａｒｅｄ  ｓｐｅｃｔｒｏｓｃｏｐｙ  ｏｆ  ｔｈｅ  ｂｒａｉｎ．Ａｐｐｌ．Ｏｐｔ．４８，Ｄ２８０−Ｄ２９８．
【非特許文献２】Ｍ．Ａ．Ｆｒａｎｃｅｓｃｈｉｎｉ，Ｄ．Ｋ．Ｊｏｓｅｐｈ，Ｔ．Ｊ．Ｈｕｐｐｅｒｔ，Ｓ．Ｇ．Ｄｉａｍｏｎｄ，ａｎｄ  Ｄ．Ａ．Ｂｏａｓ，“Ｄｉｆｆｕｓｅ  ｏｐｔｉｃａｌ  ｉｍａｇｉｎｇ  ｏｆ  ｔｈｅ  ｗｈｏｌｅ  ｈｅａｄ，”Ｊ．Ｂｉｏｍｅｄ．Ｏｐｔ．１１，０５４００７，２００６．
【非特許文献３】Ｙａｍａｄａ  Ｔ，Ｕｍｅｙａｍａ  Ｓ，Ｍａｔｓｕｄａ  Ｋ．Ｍｕｌｔｉｄｉｓｔａｎｃｅ  ｐｒｏｂｅ  ａｒｒａｎｇｅｍｅｎｔ  ｔｏ  ｅｌｉｍｉｎａｔｅ  ａｒｔｉｆａｃｔｓ  ｉｎ  ｆｕｎｃｔｉｏｎａｌ  ｎｅａｒ−ｉｎｆｒａｒｅｄ  ｓｐｅｃｔｒｏｓｃｏｐｙ．Ｊ  Ｂｉｏｍｅｄ  Ｏｐｔ．２００９  Ｎｏｖ−Ｄｅｃ；１４（６）：０６４０３４．
【非特許文献４】Ｚｈａｎｇ  Ｑ，Ｂｒｏｗｎ  ＥＮ，Ｓｔｒａｎｇｍａｎ  ＧＥ．Ａｄａｐｔｉｖｅ  ｆｉｌｔｅｒｉｎｇ  ｆｏｒ  ｇｌｏｂａｌ  ｉｎｔｅｒｆｅｒｅｎｃｅ  ｃａｎｃｅｌｌａｔｉｏｎ  ａｎｄ  ｒｅａｌ−ｔｉｍｅ  ｒｅｃｏｖｅｒｙ  ｏｆ  ｅｖｏｋｅｄ  ｂｒａｉｎ  ａｃｔｉｖｉｔｙ：ａ  Ｍｏｎｔｅ  Ｃａｒｌｏ  ｓｉｍｕｌａｔｉｏｎ  ｓｔｕｄｙ．Ｊ  Ｂｉｏｍｅｄ  Ｏｐｔ．２００７  Ｊｕｌ−Ａｕｇ；１２（４）：０４４０１４．
【非特許文献５】Ｓ．Ｋｏｈｎｏ，Ｉ．Ｍｉｙａｉ，Ａ．Ｓｅｉｙａｍａ，Ｉ．Ｏｄａ，Ａ．Ｉｓｈｉｋａｗａ，Ｓ．Ｔｓｕｎｅｉｓｈｉ，Ｔ．Ａｍｉｔａ，ａｎｄ  Ｋ．Ｓｈｉｍｉｚｕ，“Ｒｅ−ｍｏｖａｌ  ｏｆ  ｔｈｅ  ｓｋｉｎ  ｂｌｏｏｄ  ｆｌｏｗ  ａｒｔｉｆａｃｔ  ｉｎ  ｆｕｎｃｔｉｏｎａｌ  ｎｅａｒ−ｉｎｆｒａｒｅｄ  ｓｐｅｃｔｒｏｓｃｏｐｉｃ  ｉｍａｇｉｎｇ  ｄａｔａ  ｔｈｒｏｕｇｈ  ｉｎｄｅｐｅｎｄｅｎｔ  ｃｏｍｐｏｎｅｎｔ  ａｎａｌｙｓｉｓ，”Ｊ．Ｂｉｏｍｅｄ．Ｏｐｔ．１２，０６２１１１  ２００７．
【非特許文献６】Ｓ．Ｂｏｄｅｎ，Ｈ．Ｏｂｒｉｇ，Ｃ．Ｋｏｈｎｃｋｅ（「ｏ」はウムラウト付き），Ｈ．Ｂｅｎａｖ，Ｓ．Ｐ．Ｋｏｃｈ，ａｎｄ  Ｊ．Ｓｔｅｉｎｂｒｉｎｋ，Ｔｈｅ  ｏｘｙｇｅｎａｔｉｏｎ  ｒｅｓｐｏｎｓｅ  ｔｏ  ｆｕｎｃｔｉｏｎａｌ  ｓｔｉｍｕｌａｔｉｏｎ：Ｉｓ  ｔｈｅｒｅ  ａ  ｐｈｙｓｉｏｌｏｇｉｃａｌ  ｍｅａｎｉｎｇ  ｔｏ  ｔｈｅ  ｌａｇ  ｂｅｔｗｅｅｎ  ｐａｒａｍｅｔｅｒｓ？ＮｅｕｒｏＩｍａｇｅ  ３６（２００７）１００−１０７
【非特許文献７】Ｙａｍａｍｏｔｏ，Ｔ．，Ｋａｔｏ，Ｔ．，２００２．Ｐａｒａｄｏｘｉｃａｌ  ｃｏｒｒｅｌａｔｉｏｎ  ｂｅｔｗｅｅｎ  ｓｉｇｎａｌ  ｉｎ  ｆｕｎｃｔｉｏｎａｌ  ｍａｇｎｅｔｉｃ  ｒｅｓｏｎａｎｃｅ  ｉｍａｇｉｎｇ  ａｎｄ  ｄｅｏｘｙｇｅｎａｔｅｄ  ｈａｅｍｏｇｌｏｂｉｎ  ｃｏｎｔｅｎｔ  ｉｎ  ｃａｐｉｌｌａｒｉｅｓ：ａ  ｎｅｗ  ｔｈｅｏｒｅｔｉｃａｌ  ｅｘｐｌａｎａｔｉｏｎ．Ｐｈｙｓ．Ｍｅｄ．Ｂｉｏｌ．４７，１１２１−１１４１．
【非特許文献８】Ｓｔｒａｎｇｍａｎ，Ｇ．，Ｃｕｌｖｅｒ，Ｊ．Ｐ．，Ｔｈｏｍｐｓｏｎ，Ｊ．Ｈ．，Ｂｏａｓ，Ｄ．Ａ．，２００２．Ａ  ｑｕａｌｉｔａｔｉｖｅ  ｃｏｍｐａｒｉｓｏｎ  ｏｆ  ｓｉｍｕｌｔａｎｅｏｕｓ  ＢＯＬＤ  ｆＭＲＩ  ａｎｄ  ＮＩＲＳ  ｒｅｃｏｒｄｉｎｇｓ  ｄｕｒｉｎｇ  ｆｕｎｃｔｉｏｎａｌ  ｂｒａｉｎ  ａｃｔｉｖｉｔｙ．ＮｅｕｒｏＩｍａｇｅ  １７（２），７１９−７３１．
【非特許文献９】Ｈｕｐｐｅｒｔ  ＴＪ，Ｈｏｇｅ  ＲＤ，Ｄｉａｍｏｎｄ  ＳＧ，Ｆｒａｎｃｅｓｃｈｉｎｉ  ＭＡ，Ｂｏａｓ  ＤＡ．Ａ  ｔｅｍｐｏｒａｌ  ｃｏｍｐａｒｉｓｏｎ  ｏｆ  ＢＯＬＤ，ＡＳＬ，ａｎｄ  ＮＩＲＳ  ｈｅｍｏｄｙｎａｍｉｃ  ｒｅｓｐｏｎｓｅｓ  ｔｏ  ｍｏｔｏｒ  ｓｔｉｍｕｌｉ  ｉｎ  ａｄｕｌｔ  ｈｕｍａｎｓ．Ｎｅｕｒｏｉｍａｇｅ．２００６  Ｊａｎ１５；２９（２）：３６８−８２．Ｅｐｕｂ  ２００５  Ｎｏｖ  ２１．
【非特許文献１０】Ｓｈｅｔｈ  ＳＡ，Ｎｅｍｏｔｏ  Ｍ，Ｇｕｉｏｕ  Ｍ，Ｗａｌｋｅｒ  Ｍ，Ｐｏｕｒａｔｉａｎ  Ｎ，Ｔｏｇａ  ＡＷ．Ｌｉｎｅａｒ  ａｎｄ  Ｎｏｎｌｉｎｅａｒ  Ｒｅｌａｔｉｏｎｓｈｉｐｓ  ｂｅｔｗｅｅｎ  Ｎｅｕｒｏｎａｌ  Ａｃｔｉｖｉｔｙ，Ｏｘｙｇｅｎ  Ｍｅｔａｂｏｌｉｓｍ，ａｎｄ  Ｈｅｍｏｄｙｎａｍｉｃ  Ｒｅｓｐｏｎｓｅｓ  Ｎｅｕｒｏｎ．２００４  Ａｐｒ  ２２；４２（２）：３４７−５５．
【非特許文献１１】Ｈｕｐｐｅｒｔ  ＴＪ，Ａｌｌｅｎ  ＭＳ，Ｂｅｎａｖ  Ｈ，Ｊｏｎｅｓ  ＰＢ，Ｂｏａｓ  ＤＡ．Ａ  ｍｕｌｔｉｃｏｍｐａｒｔｍｅｎｔ  ｖａｓｃｕｌａｒ  ｍｏｄｅｌ  ｆｏｒ  ｉｎｆｅｒｒｉｎｇ  ｂａｓｅｌｉｎｅ  ａｎｄ  ｆｕｎｃｔｉｏｎａｌ  ｃｈａｎｇｅｓ  ｉｎ  ｃｅｒｅｂｒａｌ  ｏｘｙｇｅｎ  ｍｅｔａｂｏｌｉｓｍ  ａｎｄ  ａｒｔｅｒｉａｌ  ｄｉｌａｔｉｏｎ．Ｊ  Ｃｅｒｅｂ  Ｂｌｏｏｄ  Ｆｌｏｗ  Ｍｅｔａｂ．２００７  Ｊｕｎ；２７（６）：１２６２−７９．Ｅｐｕｂ  ２００７  Ｊａｎ  ３．
【非特許文献１２】Ｍａｎｄｅｖｉｌｌｅ  ＪＢ，Ｍａｒｏｔａ  ＪＪＡ，Ａｙａｔａ  Ｃ，Ｚａｈａｒｃｈｕｋ  Ｇ，Ｍｏｓｋｏｗｉｔｚ  ＭＡ，Ｒｏｓｅｎ  ＢＲ，Ｗｅｉｓｓｋｏｆｆ  ＲＭ  Ｅｖｉｄｅｎｃｅ  ｏｆ  ａ  ｃｅｒｅｂｒｏｖａｓｃｕｌａｒ  ｐｏｓｔａｒｔｅｒｉｏｌｅ  ｗｉｎｄｋｅｓｓｅｌ  ｗｉｔｈ  ｄｅｌａｙｅｄ  ｃｏｍｐｌｉａｎｃｅ．Ｊｏｕｒｎａｌ  ｏｆ  Ｃｅｒｅｂｒａｌ  Ｂｌｏｏｄ  Ｆｌｏｗ  ＆  Ｍｅｔａｂｏｌｉｓｍ（１９９９）１９，６７９−６８９
【非特許文献１３】Ｏｓｔｅｒｇａａｒｄ  Ｌ，Ｓｍｉｔｈ  ＤＦ，Ｖｅｓｔｅｒｇａａｒｄ−Ｐｏｕｌｓｅｎ  Ｐ，Ｈａｎｓｅｎ  ＳＢ，Ｇｅｅ  ＡＤ，Ｇｊｅｄｄｅ  Ａ，Ｇｙｌｄｅｎｓｔｅｄ  Ｃ，Ａｂｓｏｌｕｔｅ  ｃｅｒｅｂｒａｌ  ｂｌｏｏｄ  ｆｌｏｗ  ａｎｄ  ｂｌｏｏｄ  ｖｏｌｕｍｅ  ｍｅａｓｕｒｅｄ  ｂｙ  ｍａｇｎｅｔｉｃ  ｒｅｓｏｎａｎｃｅ  ｉｍａｇｉｎｇ  ｂｏｌｕｓ  ｔｒａｃｋｉｎｇ：Ｃｏｍｐａｒｉｓｏｎ  ｗｉｔｈ  ｐｏｓｉｔｒｏｎ  ｅｍｉｓｓｉｏｎ  ｔｏｍｏｇｒａｐｈｙ  ｖａｌｕｅｓ．Ｊｏｕｒｎａｌ  ｏｆ  Ｃｅｒｅｂｒａｌ  Ｂｌｏｏｄ  Ｆｌｏｗ  ＆  Ｍｅｔａｂｏｌｉｓｍ（１９９８）１８，４２５−４３２
【非特許文献１４】Ｓｃｈｕｍａｎｎ  Ｐ，Ｔｏｕｚａｎｉ  Ｏ，Ｙｏｕｎｇ  ＡＲ，Ｂａｒｏｎ  ＪＣ，Ｍｏｒｅｌｌｏ  Ｒ，ＭａｃＫｅｎｚｉｅ  ＥＴ  Ｅｖａｌｕａｔｉｏｎ  ｏｆ  ｔｈｅ  ｒａｔｉｏ  ｏｆ  ｃｅｒｅｂｒａｌ  ｂｌｏｏｄ  ｆｌｏｗ  ｔｏ  ｃｅｒｅｂｒａｌ  ｂｌｏｏｄ  ｖｏｌｕｍｅ  ａｓ  ａｎ  ｉｎｄｅｘ  ｏｆ  ｌｏｃａｌ  ｃｅｒｅｂｒａｌ  ｐｅｒｆｕｓｉｏｎ  ｐｒｅｓｓｕｒｅ．Ｂｒａｉｎ  （１９９８），１２１，１３６９−１３７９
【非特許文献１５】Ｚｏｎｔａ  Ｍ，Ａｎｇｕｌｏ  ＭＣ，Ｇｏｂｂｏ  Ｓ，Ｒｏｓｅｎｇａｒｔｅｎ  Ｂ，Ｈｏｓｓｍａｎｎ  ＫＡ，  Ｐｏｚｚａｎ  Ｔ，Ｃａｒｍｉｇｎｏｔｏ  Ｇ．Ｎｅｕｒｏｎ−  ｔｏ−ａｓｔｒｏｃｙｔｅ  ｓｉｇｎａｌｉｎｇ  ｉｓ  ｃｅｎｔｒａｌ  ｔｏ  ｔｈｅ  ｄｙｎａｍｉｃ  ｃｏｎｔｒｏｌ  ｏｆ  ｂｒａｉｎ  ｍｉｃｒｏｃｉｒｃｕｌａｔｉｏｎ．Ｎａｔ  Ｎｅｕｒｏｓｃｉ．２００３  Ｊａｎ；６（１）：４３−５０．
【非特許文献１６】Ｉｎｔａｇｌｉｅｔｔａ  Ｍ，Ｊｏｈｎｓｏｎ  ＰＣ，Ｗｉｎｓｌｏｗ  ＲＭ  Ｍｉｃｒｏｖａｓｃｕｌａｒ  ａｎｄ  ｔｉｓｓｕｅ  ｏｘｙｇｅｎ  ｄｉｓｔｒｉｂｕｔｉｏｎ．Ｃａｒｄｉｏｖａｓｃｕｌａｒ  Ｒｅｓｅａｒｃｈ  ３２（１９９６）６３２−６４３
【非特許文献１７】Ｐｏｐｅｌ  ＡＳ，Ｊｏｈｎｓｏｎ  ＰＣ  Ｍｉｃｒｏｃｉｒｃｕｌａｔｉｏｎ  ａｎｄ  ｈｅｍｏｒｈｅｏｌｏｇｙ．Ａｎｎｕ．Ｒｅｖ．Ｆｌｕｉｄ  Ｍｅｃｈ．２００５．３７：４３−６９
【非特許文献１８】Ｏｂｒｉｇ  Ｈ，  Ｖｉｌｌｒｉｎｇｅｒ  Ａ．Ｂｅｙｏｎｄ  ｔｈｅ  ｖｉｓｉｂｌｅｉｍａｇｉｎｇ  ｔｈｅ  ｈｕｍａｎ  ｂｒａｉｎ  ｗｉｔｈ  ｌｉｇｈｔ．Ｊ  Ｃｅｒｅｂ  Ｂｌｏｏｄ  Ｆｌｏｗ  Ｍｅｔａｂ．２００３  Ｊａｎ；２３（１）：１−１８．
【非特許文献１９】Ｊａｓｄｚｅｗｓｋｉ  Ｇ，Ｓｔｒａｎｇｍａｎ  Ｇ，Ｗａｇｎｅｒ  Ｊ，Ｋｗｏｎｇ  ＫＫ，Ｐｏｌｄｒａｃｋ  ＲＡ，Ｂｏａｓ  ＤＡ．Ｄｉｆｆｅｒｅｎｃｅｓ  ｉｎ  ｔｈｅ  ｈｅｍｏｄｙｎａｍｉｃ  ｒｅｓｐｏｎｓｅ  ｔｏ  ｅｖｅｎｔ−ｒｅｌａｔｅｄ  ｍｏｔｏｒ  ａｎｄ  ｖｉｓｕａｌ  ｐａｒａｄｉｇｍｓ  ａｓ  ｍｅａｓｕｒｅｄ  ｂｙ  ｎｅａｒ−ｉｎｆｒａｒｅｄｓｐｅｃｔｒｏｓｃｏｐｙ．Ｎｅｕｒｏｉｍａｇｅ．２００３  Ｓｅｐ；２０（１）：４７９−８８．
【非特許文献２０】Ｔａｎｇ  Ｌ，Ａｖｉｓｏｎ  ＭＪ，Ｇｏｒｅ  ＪＣ．Ｎｏｎｌｉｎｅａｒ  ｂｌｏｏｄ  ｏｘｙｇｅｎ  ｌｅｖｅｌ−ｄｅｐｅｎｄｅｎｔ  ｒｅｓｐｏｎｓｅｓ  ｆｏｒ  ｔｒａｎｓｉｅｎｔ  ａｃｔｉｖａｔｉｏｎｓ  ａｎｄ  ｄｅａｃｔｉｖａｔｉｏｎｓ  ｉｎ  Ｖ１‐ｉｎｓｉｇｈｔｓ  ｉｎｔｏ  ｔｈｅ  ｈｅｍｏｄｙｎａｍｉｃ  ｒｅｓｐｏｎｓｅ  ｆｕｎｃｔｉｏｎ  ｗｉｔｈ  ｔｈｅ  ｂａｌｌｏｏｎ  ｍｏｄｅｌ．Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ｉｍａｇｉｎｇ  ２７（２００９）４４９−４５９
【非特許文献２１】ＭｃＩｎｔｏｓｈ  ＭＡ，Ｓｈａｈａｎｉ  Ｕ，Ｂｏｕｌｔｏｎ  ＲＧ，ＭｃＣｕｌｌｏｃｈ  ＤＬ．Ａｂｓｏｌｕｔｅ  ｑｕａｎｔｉｆｉｃａｔｉｏｎ  ｏｆ  ｏｘｙｇｅｎａｔｅｄ  ｈｅｍｏｇｌｏｂｉｎ  ｗｉｔｈｉｎ  ｔｈｅ  ｖｉｓｕａｌ  ｃｏｒｔｅｘ  ｗｉｔｈ  ｆｕｎｃｔｉｏｎａｌ  ｎｅａｒ  ｉｎｆｒａｒｅｄ  ｓｐｅｃｔｒｏｓｃｏｐｙ  （ｆＮＩＲＳ）．Ｉｎｖｅｓｔ  Ｏｐｈｔｈａｌｍｏｌ  Ｖｉｓ  Ｓｃｉ．２０１０  Ｓｅｐ；５１（９）：４８５６−６０．Ｅｐｕｂ  ２０１０  Ｍａｒ  ３１．
【非特許文献２２】Ｓｉｅｇｅｌ  ＡＭ，Ｃｕｌｖｅｒ  ＪＰ，Ｍａｎｄｅｖｉｌｌｅ  ＪＢ，Ｂｏａｓ  ＤＡ．Ｔｅｍｐｏｒａｌ  ｃｏｍｐａｒｉｓｏｎ  ｏｆ  ｆｕｎｃｔｉｏｎａｌ  ｂｒａｉｎ  ｉｍａｇｉｎｇ  ｗｉｔｈ  ｄｉｆｆｕｓｅ  ｏｐｔｉｃａｌ  ｔｏｍｏｇｒａｐｈｙ  ａｎｄ  ｆＭＲＩ  ｄｕｒｉｎｇ  ｒａｔ  ｆｏｒｅｐａｗ  ｓｔｉｍｕｌａｔｉｏｎ．Ｐｈｙｓ  Ｍｅｄ  Ｂｉｏｌ．２００３  Ｍａｙ  ２１；４８（１０）：１３９１−４０３．
【非特許文献２３】Ｖｉｌｌｒｉｎｇｅｒ  Ａ，Ｐｌａｎｃｋ  Ｊ，Ｈｏｃｋ  Ｃ，Ｓｃｈｌｅｉｎｋｏｆｅｒ  Ｌ，Ｄｉｒｎａｇｌ  Ｕ．Ｎｅａｒ  ｉｎｆｒａｒｅｄ  ｓｐｅｃｔｒｏｓｃｏｐｙ（ＮＩＲＳ）：ａ  ｎｅｗ  ｔｏｏｌ  ｔｏ  ｓｔｕｄｙ  ｈｅｍｏｄｙｎａｍｉｃ  ｃｈａｎｇｅｓ  ｄｕｒｉｎｇ  ａｃｔｉｖａｔｉｏｎ  ｏｆ  ｂｒａｉｎ  ｆｕｎｃｔｉｏｎ  ｉｎ  ｈｕｍａｎ  ａｄｕｌｔｓ．Ｎｅｕｒｏｓｃｉｅｎｃｅ  Ｌｅｔｔｅｒｓ，１５４（１９９３）１０１−１０４
【非特許文献２４】Ｏｂｒｉｇ  Ｈ，Ｗｅｎｚｅｌ  Ｒ，Ｋｏｈｌ  Ｍ，Ｈｏｒｓｔ  Ｓ，Ｗｏｂｓｔ  Ｐ，Ｓｔｅｉｎｂｒｉｎｋ  Ｊ，Ｔｈｏｍａｓ  Ｆ，Ｖｉｌｌｒｉｎｇｅｒ  Ａ．Ｎｅａｒ−ｉｎｆｒａｒｅｄ  ｓｐｅｃｔｒｏｓｃｏｐｙ：ｄｏｅｓ  ｉｔ  ｆｕｎｃｔｉｏｎ  ｉｎ  ｆｕｎｃｔｉｏｎａｌ  ａｃｔｉｖａｔｉｏｎ  ｓｔｕｄｉｅｓ  ｏｆ  ｔｈｅ  ａｄｕｌｔ  ｂｒａｉｎ？Ｉｎｔ  Ｊ  Ｐｓｙｃｈｏｐｈｙｓｉｏｌ．２０００  Ｍａｒ；３５（２−３）：１２５−４２．
【非特許文献２５】Ｚｅｆｆ  ＢＷ，Ｗｈｉｔｅ  ＢＲ，Ｄｅｈｇｈａｎｉ  Ｈ，Ｓｃｈｌａｇｇａｒ  ＢＬ，Ｃｕｌｖｅｒ  ＪＰ．Ｒｅｔｉｎｏｔｏｐｉｃ  ｍａｐｐｉｎｇ  ｏｆ  ａｄｕｌｔ  ｈｕｍａｎ  ｖｉｓｕａｌ  ｃｏｒｔｅｘ  ｗｉｔｈ  ｈｉｇｈ−ｄｅｎｓｉｔｙ  ｄｉｆｆｕｓｅ  ｏｐｔｉｃａｌ  ｔｏｍｏｇｒａｐｈｙ．Ｐｒｏｃ  Ｎａｔｌ  Ａｃａｄ  Ｓｃｉ  ＵＳＡ．２００７  Ｊｕｌ  １７；１０４（２９）：１２１６９−７４．Ｅｐｕｂ  ２００７  Ｊｕｌ  ６．
【非特許文献２６】Ｄｕｎｎ  ＡＫ，Ｄｅｖｏｒ  Ａ，Ｄａｌｅ  ＡＭ，Ｂｏａｓ  ＤＡ．Ｓｐａｔｉａｌ  ｅｘｔｅｎｔ  ｏｆ  ｏｘｙｇｅｎ  ｍｅｔａｂｏｌｉｓｍ  ａｎｄ  ｈｅｍｏｄｙｎａｍｉｃ  ｃｈａｎｇｅｓ  ｄｕｒｉｎｇ  ｆｕｎｃｔｉｏｎａｌ  ａｃｔｉｖａｔｉｏｎ  ｏｆ  ｔｈｅ  ｒａｔ  ｓｏｍａｔｏｓｅｎｓｏｒｙ  ｃｏｒｔｅｘ．Ｎｅｕｒｏｉｍａｇｅ．２００５  Ａｕｇ１５；２７（２）：２７９−９０．
【非特許文献２７】Ｃｏｌｉｅｒ  ＷＮ，Ｑｕａｒｅｓｉｍａ  Ｖ，Ｗｅｎｚｅｌ  Ｒ，  ｖａｎ  ｄｅｒ  Ｓｌｕｉｊｓ  ＭＣ，Ｏｅｓｅｂｕｒｇ  Ｂ，Ｆｅｒｒａｒｉ  Ｍ，Ｖｉｌｌｒｉｎｇｅｒ  Ａ．Ｓｉｍｕｌｔａｎｅｏｕｓ  ｎｅａｒ−ｉｎｆｒａｒｅｄ  ｓｐｅｃｔｒｏｓｃｏｐｙ  ｍｏｎｉｔｏｒｉｎｇ  ｏｆ  ｌｅｆｔ  ａｎｄ  ｒｉｇｈｔ  ｏｃｃｉｐｉｔａｌ  ａｒｅａｓ  ｒｅｖｅａｌｓ  ｃｏｎｔｒａｌａｔｅｒａｌ  ｈｅｍｏｄｙｎａｍｉｃ  ｃｈａｎｇｅｓ  ｕｐｏｎ  ｈｅｍｉ−ｆｉｅｌｄ  ｐａｒａｄｉｇｍ．Ｖｉｓｉｏｎ  Ｒｅｓ．２００１  Ｊａｎ；４１（１）：９７−１０２．
【非特許文献２８】Ｂｅｒｗｉｃｋ  Ｊ，Ｍａｒｔｉｎ  Ｃ，Ｍａｒｔｉｎｄａｌｅ  Ｊ，Ｊｏｎｅｓ  Ｍ，Ｊｏｈｎｓｔｏｎ  Ｄ，Ｚｈｅｎｇ  Ｙ，Ｒｅｄｇｒａｖｅ  Ｐ，Ｍａｙｈｅｗ  Ｊ．Ｈｅｍｏｄｙｎａｍｉｃ  ｒｅｓｐｏｎｓｅ  ｉｎ  ｔｈｅ  ｕｎａｎｅｓｔｈｅｔｉｚｅｄ  ｒａｔ：ｉｎｔｒｉｎｓｉｃ  ｏｐｔｉｃａｌ  ｉｍａｇｉｎｇ  ａｎｄ  ｓｐｅｃｔｒｏｓｃｏｐｙ  ｏｆ  ｔｈｅ  ｂａｒｒｅｌ  ｃｏｒｔｅｘ．Ｊ  Ｃｅｒｅｂ  Ｂｌｏｏｄ  Ｆｌｏｗ  Ｍｅｔａｂ．２００２  Ｊｕｎ；２２（６）：６７０−９．
【非特許文献２９】Ｆｒａｎｃｅｓｃｈｉｎｉ  ＭＡ，Ｆａｎｔｉｎｉ  Ｓ，Ｔｈｏｍｐｓｏｎ  ＪＨ，Ｃｕｌｖｅｒ  ＪＰ，Ｂｏａｓ  ＤＡ．Ｈｅｍｏｄｙｎａｍｉｃ  ｅｖｏｋｅｄ  ｒｅｓｐｏｎｓｅ  ｏｆ  ｔｈｅ  ｓｅｎｓｏｒｉｍｏｔｏｒ  ｃｏｒｔｅｘ  ｍｅａｓｕｒｅｄ  ｎｏｎｉｎｖａｓｉｖｅｌｙ  ｗｉｔｈ  ｎｅａｒ−ｉｎｆｒａｒｅｄ  ｏｐｔｉｃａｌ  ｉｍａｇｉｎｇ．Ｐｓｙｃｈｏｐｈｙｓｉｏｌｏｇｙ．２００３  Ｊｕｌ；４０（４）：５４８−６０．
【非特許文献３０】Ｍａｄｓｅｎ  ＰＬ，Ｓｅｃｈｅｒ  ＮＨ．Ｎｅａｒ−ｉｎｆｒａｒｅｄ  ｏｘｉｍｅｔｒｙ  ｏｆ  ｔｈｅ  ｂｒａｉｎ．Ｐｒｏｇ  Ｎｅｕｒｏｂｉｏｌ．１９９９  Ａｕｇ；５８（６）：５４１−６０．
【非特許文献３１】Ｙａｍａｄａ  Ｔ，Ｕｍｅｙａｍａ  Ｓ，Ｍａｔｓｕｄａ  Ｋ．Ａ  ｍｕｌｔｉｄｉｓｔａｎｃｅ  ｐｒｏｂｅ  ａｒｒａｎｇｅｍｅｎｔ  ＮＩＲＳ  ｆｏｒ  ｄｅｔｅｃｔｉｎｇ  ａｂｓｏｒｐｔｉｏｎ  ｃｈａｎｇｅｓ  ｉｎ  ｃｅｒｅｂｒａｌ  ｇｒａｙ  ｍａｔｔｅｒ  ｌａｙｅｒ．Ｐｒｏｃ．  ｏｆ  ＳＰＩＥ  Ｖｏｌ．７５５７  ７５５７０Ｘ
【非特許文献３２】Ｏｋａｄａ  Ｅ，Ｆｉｒｂａｎｋ  Ｍ，Ｓｃｈｗｅｉｇｅｒ  Ｍ，Ａｒｒｉｄｇｅ  ＳＲ，Ｃｏｐｅ  Ｍ，Ｄｅｌｐｙ  ＤＴ．Ｔｈｅｏｒｅｔｉｃａｌ  ａｎｄ  ｅｘｐｅｒｉｍｅｎｔａｌ  ｉｎｖｅｓｔｉｇａｔｉｏｎ  ｏｆ  ｎｅａｒ−ｉｎｆｒａｒｅｄ  ｌｉｇｈｔ  ｐｒｏｐａｇａｔｉｏｎ  ｉｎ  ａ  ｍｏｄｅｌ  ｏｆ  ｔｈｅ  ａｄｕｌｔ  ｈｅａｄ．Ａｐｐｌ  Ｏｐｔ．１９９７  Ｊａｎ  １；３６（１）：２１−３１．
【非特許文献３３】Ｏｋｕｉ  Ｎ，Ｏｋａｄａ  Ｅ．Ｗａｖｅｌｅｎｇｔｈ  ｄｅｐｅｎｄｅｎｃｅ  ｏｆ  ｃｒｏｓｓｔａｌｋ  ｉｎ  ｄｕａｌ−ｗａｖｅｌｅｎｇｔｈ  ｍｅａｓｕｒｅｍｅｎｔ  ｏｆ  ｏｘｙ−  ａｎｄ  ｄｅｏｘｙ−ｈｅｍｏｇｌｏｂｉｎ．Ｊ  Ｂｉｏｍｅｄ  Ｏｐｔ．２００５  Ｊａｎ−Ｆｅｂ；１０（１）：１１０１５．
【非特許文献３４】Ｃｕｉ  Ｘ，Ｂｒａｙ  Ｓ，Ｒｅｉｓｓ  ＡＬ．Ｆｕｎｃｔｉｏｎａｌ  ｎｅａｒ  ｉｎｆｒａｒｅｄ  ｓｐｅｃｔｒｏｓｃｏｐｙ（ＮＩＲＳ）  ｓｉｇｎａｌ  ｉｍｐｒｏｖｅｍｅｎｔ  ｂａｓｅｄ  ｏｎ  ｎｅｇａｔｉｖｅ  ｃｏｒｒｅｌａｔｉｏｎ  ｂｅｔｗｅｅｎ  ｏｘｙｇｅｎａｔｅｄ  ａｎｄ  ｄｅｏｘｙｇｅｎａｔｅｄ  ｈｅｍｏｇｌｏｂｉｎ  ｄｙｎａｍｉｃｓ．ＮｅｕｒｏＩｍａｇｅ  ４９（２０１０）３０３９−３０４６
【図１】
0005641355.tif 000019
【図２】
0005641355.tif 000020
【図３】
0005641355.tif 000021
【図４】
0005641355.tif 000022
【図５】
0005641355.tif 000023
【図６】
0005641355.tif 000024
【図７】
0005641355.tif 000025
【図８】
0005641355.tif 000026
【図９】
0005641355.tif 000027
【図１０】
0005641355.tif 000028
【図１１】
0005641355.tif 000029
【図１２】
0005641355.tif 000030
【図１３】
0005641355.tif 000031