(57)【要約】較正作業を不要にし、若しくは最小限に留めることのできる脳機能計測装置と脳機能計測方法を提供することを課題とする。  本課題を解決するため、被検体の脳に光を照射する第一の光照射プローブＳ１と、第一の光照射プローブS1に隣接して配設され、第一の光照射プローブS1から照射された光のうち脳で反射された光を検出する第一の光検出プローブD1と、第一の光検出プローブD1に隣接して配設され、被検体の脳に光を照射する第二の光照射プローブS2と、第二の光照射プローブS2に隣接して配設され、第二の光照射プローブS2から照射された光のうち脳で反射された光を検出する第二の光検出プローブD2と、第一の光検出プローブD1と第二の光照射プローブS2との間のチャンネルにつき計測される光量が、第一の光照射プローブS1と第一の光検出プローブD1との間のチャンネルにつき観測された観測値になるよう、第二の光照射プローブS2により照射される光量を調節すると共に、第二の光照射プローブS2と第二の光検出プローブD2との間のチャンネルにつき計測される光量が上記観測値になるよう、第二の光検出プローブD2により検出される光量を調節する制御部３を備えた脳機能計測装置を提供する。
2018190130.tif 000002
【特許請求の範囲】
【請求項１】
被検体の脳に光を照射する第一の照射手段と、
前記第一の照射手段に隣接して配設され、前記第一の照射手段から照射された光のうち前記脳で反射された光を検出する第一の検出手段と、
前記第一の検出手段に隣接して配設され、前記被検体の脳に光を照射する第二の照射手段と、
前記第二の照射手段に隣接して配設され、前記第二の照射手段から照射された光のうち前記脳で反射された光を検出する第二の検出手段と、
前記第一の検出手段と前記第二の照射手段との間の計測点につき計測される光量が、前記第一の照射手段と前記第一の検出手段との間の計測点につき観測された観測値になるよう、前記第二の照射手段により照射される光量につき第一の調節を行うと共に、前記第二の照射手段と前記第二の検出手段との間の計測点につき計測される光量が前記観測値になるよう、前記第二の検出手段により検出される光量につき第二の調節を行う制御手段とを備えた脳機能計測装置。
【請求項２】
前記脳機能計測装置は、複数の前記照射手段及び前記検出手段を備え、  
前記制御手段は、複数の前記照射手段の中から最大の実効入射光量を持つ前記照射手段を選択すると共に、複数の前記検出手段の中から最大の実効検出効率を持つ前記検出手段を選択し、選択された前記照射手段と前記検出手段を結ぶ最短経路に位置する各々の前記計測点について、前記選択された前記照射手段から前記選択された前記検出手段に向かう順序で前記第一及び第二の調節を実行した後、前記最短経路に位置する各々の前記計測点について、前記選択された前記検出手段から前記選択された前記照射手段に向かう順序で前記第二及び第一の調節を実行する、請求項１に記載の脳機能計測装置。
【請求項３】
前記制御手段は、全ての使用波長毎に前記調節を行う、請求項２に記載の脳機能計測装置。
【請求項４】
前記最大の実効入射光量を持つ照射手段が照射する前記計測点、若しくは前記最大の実効検出効率を持つ検出手段が検出する前記計測点の少なくとも一方は、前記被検体の無毛部に配設された、請求項２に記載の脳機能計測装置。
【請求項５】
記憶手段をさらに備え、
前記制御手段は、さらに、
前記第一の検出手段と前記第二の照射手段との間の計測点につき計測される光量が、前記第一の照射手段と前記第一の検出手段との間の計測点につき観測された観測値になるよう、前記第二の照射手段により照射される光量を調節すると共に、前記第二の照射手段と前記第二の検出手段との間の計測点につき計測される光量が前記観測値になるよう、前記第二の検出手段により検出される光量を調節した後、前記観測値を基準値として前記記憶手段に保存し、かつ、少なくとも一つの前記検出手段へ入射する光の透過率を百パーセントとしたときに、前記少なくとも一つの前記検出手段が計測する前記計測点につき観測された最大光量値を前記記憶手段に保存した上で、
平準化対象とする前記計測点につき計測するための少なくとも一つの前記検出手段へ入射する光の透過率を百パーセントとし、前記平準化対象とする前記計測点につき観測される光量が前記最大光量値となるように、前記平準化対象とする前記計測点につき計測するための前記照射手段から照射される光量を調節した後に、前記平準化対象とする前記計測点につき観測される光量の値が前記基準値となるように、前記平準化対象とする前記計測点につき計測するための前記検出手段で検出される光量を調節し、前記平準化対象とする前記計測点を、前記第一の検出手段と第一の照射手段との間の計測点とみなして前記照射手段の調節を行う、請求項４に記載の脳機能計測装置。
【請求項６】
記憶手段をさらに備え、
前記制御部は、さらに、
前記第一の検出手段と前記第二の照射手段との間の計測点につき計測される光量が、前記第一の照射手段と前記第一の検出手段との間の計測点につき観測された観測値になるよう、前記第二の照射手段により照射される光量を調節すると共に、前記第二の照射手段と前記第二の検出手段との間の計測点につき計測される光量が前記観測値になるよう、前記第二の検出手段により検出される光量を調節した時点における前記第一の検出手段の実効検出効率の値と前記第一の照射手段の実効入射光量の値とを前記記憶手段に保存した上で、
前記第一の検出手段における検出効率を保存された前記実効検出効率とし、前記第一の照射手段における照射光量を保存された前記実効入射光量とし、再度前記照射手段に対する前記調節を行う、請求項２に記載の脳機能計測装置。
【請求項７】
被検体の脳に光を照射し、前記脳において反射された光を検出することにより前記脳の機能を計測する脳機能計測方法であって、
隣接する第一の検出手段と第二の照射手段との間の計測点につき計測される光量が、前記第一の検出手段に隣接する第一の照射手段と前記第一の検出手段との間の計測点につき観測された観測値になるよう、前記第二の照射手段により照射される光量を調節する第一のステップと、
前記第二の照射手段と前記第二の照射手段に隣接する第二の検出手段との間の計測点につき計測される光量が、前記観測値になるよう、前記第二の検出手段により検出される光量を調節する第二のステップとを有する脳機能計測方法。
【請求項８】
前記脳機能計測方法は、複数の前記計測点において前記脳の機能を計測する方法であって、
複数の前記照射手段の中から最大の実効入射光量を持つ照射手段を選択すると共に、複数の前記検出手段の中から最大の実効検出効率を持つ検出手段を選択する第三のステップと、
前記第三のステップで選択された前記照射手段と前記検出手段を結ぶ最短経路に位置する各々の前記計測点について、前記選択された前記照射手段から前記選択された前記検出手段に向かう順序で前記第一及び第二のステップを実行した後、前記最短経路に位置する各々の前記計測点について、前記選択された前記検出手段から前記選択された前記照射手段に向かう順序で前記第二及び第一のステップを実行する第四のステップと、
前記最短経路に位置する各々の前記計測点と、前記各々の前記計測点に隣接する前記計測点との間において、前記第一のステップ及び前記第二のステップを実行する第五のステップとをさらに有する、請求項７に記載の脳機能計測方法。
【請求項９】
前記第三から第五のステップを、全ての使用波長毎に実行する、請求項８に記載の脳機能計測方法。
【請求項１０】
前記最大の実効入射光量を持つ照射手段が照射する前記計測点、若しくは前記最大の実効検出効率を持つ検出手段が検出する前記計測点の少なくとも一方は、前記被検体の無毛部に配設される、請求項８に記載の脳機能計測方法。
【請求項１１】
前記第二のステップを実行した後に前記観測値を基準値として保存する第六のステップと、
少なくとも一つの前記検出手段へ入射する光の透過率を百パーセントとしたときに、前記少なくとも一つの前記検出手段が計測する前記計測点につき観測された最大光量値を保存する第七のステップと、
平準化対象とする前記計測点につき計測するための少なくとも一つの前記検出手段へ入射する光の透過率を百パーセントとし、前記平準化対象とする前記計測点につき観測される光量が前記最大光量値となるように、前記平準化対象とする前記計測点につき計測するための前記照射手段から照射される光量を調節する第八のステップと、
前記第八のステップを実行した後に、前記平準化対象とする前記計測点につき観測される光量の値が前記基準値となるように、前記平準化対象とする前記計測点につき計測するための前記検出手段で検出される光量を調節する第九のステップと、
前記平準化対象とする前記計測点を、前記第一の検出手段と第一の照射手段との間の計測点として前記第一のステップから前記第五のステップを実行する第十のステップとをさらに有する、請求項１０に記載の脳機能計測方法。
【請求項１２】
前記第五のステップを終えた時点における前記実効検出効率の値と前記実効入射光量の値とを保存する第六のステップと、
基準となる隣接した二つの前記計測点を計測するための前記検出手段における検出効率を前記第六のステップで保存された前記実効検出効率とし、前記二つの前記計測点を計測するための前記照射手段における照射光量を前記第六のステップで保存された前記実効入射光量とする第七のステップと、
前記第七のステップの後に、前記基準とした前記計測点を基点として再度前記第一から第五のステップを実行する、請求項８に記載の脳機能計測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、脳の機能を計測するための装置と方法に関するものである。
【背景技術】
【０００２】
近年においては、簡便な脳機能計測手法として機能的近赤外分光法（functional near infrared spectroscopy: fNIRS）が知られている。このfNIRSを用いた計測は、プローブを頭皮に当てて行われるが、プローブ直下の毛髪の多寡や、プローブと頭皮の密着度によって、頭皮とプローブの間で生じる光減衰の大きさが異なることになる。従って、計測点（チャンネル）毎に信号対雑音比も著しく異なるため、計測された信号に厳格な統計解析を適用することはできない。【０００３】
このような背景から、特許文献１に示されるように、各チャンネルで計測された全ての信号の雑音分散を等しくする技術が考案されている。すなわち、本技術は、各光源や各検出器と生体との間に光減衰器を導入し、これらの光透過率を制御することで全チャンネルにおける観測光量を平準化し、雑音分散の平準化を図るものである。なお、各光減衰器の光透過率は、計測値に基づいて設定される。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】ＷＯ２０１６１３２９８９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
特許文献１に示された技術においては、全ての光減衰器について、できる限り正確な較正を行う必要があるという問題がある。【０００６】
本発明は、このような問題を解決するためになされたもので、較正作業を不要にし、若しくは最小限に留めることのできる脳機能計測装置と脳機能計測方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記課題を解決するため、本発明は、被検体の脳に光を照射する第一の照射手段と、第一の照射手段に隣接して配設され、第一の照射手段から第一の方向へ照射された光のうち脳で反射された光を検出する第一の検出手段と、第一の検出手段に隣接して配設され、被検体の脳に光を照射する第二の照射手段と、第二の照射手段に隣接して配設され、第二の照射手段から照射された光のうち脳で反射された光を検出する第二の検出手段と、第一の検出手段と第二の照射手段との間の計測点につき計測される光量が、第一の照射手段と第一の検出手段との間の計測点につき観測された観測値になるよう、第二の照射手段により照射される光量を調節すると共に、第二の照射手段と第二の検出手段との間の計測点につき計測される光量が上記観測値になるよう、第二の検出手段により検出される光量を調節する制御手段を備えた脳機能計測装置を提供する。【０００８】
また、上記課題を解決するため、本発明は、被検体の脳に光を照射し、脳において反射された光を検出することにより脳の機能を計測する脳機能計測方法であって、隣接する第一の検出手段と第二の照射手段との間の計測点につき計測される光量が、第一の検出手段に隣接する第一の照射手段と第一の検出手段との間の計測点につき観測された観測値になるよう、第二の照射手段により照射される光量を調節する第一のステップと、第二の照射手段と第二の照射手段に隣接する第二の検出手段との間の計測点につき計測される光量が上記観測値になるよう、第二の検出手段により検出される光量を調節する第二のステップを有する脳機能計測方法を提供する。
【発明の効果】
【０００９】
本発明によれば、較正作業を不要にし、若しくは最小限に留めることのできる脳機能計測装置と脳機能計測方法を提供することができる。
【図面の簡単な説明】
【００１０】
【図１】本発明の実施の形態に係る脳機能計測システムの構成を示す図である。
【図２】図１に示されたホルダ２０上における光照射プローブS1,S2及び光検出プローブD1,D2の第一の配置例を示す図である。
【図３】図２に示された光照射プローブS及び光検出プローブDの配置を二次元平面に記した模式図である。
【図４】本発明の実施の形態に係る第一の脳機能計測方法を図１に示された制御部３による制御により実現する方法を示すフローチャートである。
【図５】図１に示されたホルダ２０上における光照射プローブS1,S2及び光検出プローブD1,D2の第二の配置例を示す図である。
【図６】本発明の実施の形態に係る第二の脳機能計測方法を図１に示された制御部３による制御により実現する方法を示すフローチャートである。
【図７】本発明の実施の形態に係る第三の脳機能計測方法を図１に示された制御部３による制御により実現する方法を示すフローチャートである。
【図８】図１に示されたホルダ２０上における光照射プローブS1,S2及び光検出プローブD1,D2の第三の配置例を示す図である。S1
【図９】図８に示された光照射プローブ及び光検出プローブの配置を二次元平面に記した模式図である。
【発明を実施するための形態】
【００１１】
以下において、本発明の実施の形態を図面を参照しつつ詳しく説明する。なお、図中同一符号は同一又は相当部分を示す。
【００１２】
図１は、本発明の実施の形態に係る脳機能計測システムの構成を示す図である。図１に示されるように、脳機能計測システムは、制御部３と光源部５、データ計測部７、及び記憶部９を含む脳機能計測装置１と、光減衰器１１〜１４と、被検体に装着され光照射プローブS1,S2及び光検出プローブD1,D2を含むホルダ２０とを備える。
【００１３】
ここで、制御部３は光源部５、データ計測部７、記憶部９、及び光減衰器１１〜１４に接続され、光照射プローブS1,S2はそれぞれ光減衰器１１，１２を介して光源部５に接続される。また、光検出プローブD1,D2はそれぞれ光減衰器１３，１４を介してデータ計測部７に接続される。また、データ計測部７は記憶部９にも接続される。
【００１４】
一方で、例えば、図１に示されたホルダ２０上における光照射プローブS1,S2は、図２の白丸で示された光照射プローブSの位置に配置され、図１に示されたホルダ２０上における光検出プローブD1,D2は、図２の黒丸で示された光検出プローブDの位置に配置される。なお、図２においてグレーの部分は毛髪を示し、図５及び図８でも同様である。
【００１５】
ここで、図２に示された光照射プローブS及び光検出プローブDは、例えば図３に示された模式図で示されるように交互に配置される。
【００１６】
［連結チャンネル群内の雑音平準化］
図４は、本発明の実施の形態に係る第一の脳機能計測方法を図１に示された制御部３による制御により実現する方法を示すフローチャートである。以下において、図３を参照しつつ、図４に示された第一の脳機能計測方法を詳しく説明する。なお、図４に示された第一の脳機能計測方法は、図１に示された制御部３による制御により実現する場合に限られないことは言うまでもない。
【００１７】
ステップＳ１において、制御部３は、隣接する第一の光検出プローブD1と第二の光照射プローブS2との間のチャンネルCH2につき第一の光検出プローブD1で計測される光量が、第一の光検出プローブD1に隣接する第一の光照射プローブS1と第一の光検出プローブD1との間のチャンネルCH1につき第一の光検出プローブD1で観測された光量の値（観測値）になるよう、第二の光照射プローブS2に接続された光減衰器１２を制御することにより、第二の光検出プローブD2で検出される光量を調節する。
【００１８】
次に、ステップＳ２において、制御部３は、第二の光照射プローブS2と第二の光照射プローブS2に隣接する第二の光検出プローブD2との間のチャンネルCH3につき第二の光検出プローブD2で計測される光量が上記観測値になるよう、光減衰器１４を制御することにより、第二の光検出プローブD2で検出される光量を調節する。
【００１９】
次に、ステップＳ３において、制御部３は、複数の光照射プローブの中から最大の実効入射光量を持つ光照射プローブを選択すると共に、複数の光検出プローブの中から最大の実効検出効率を持つ光検出プローブを選択する。
【００２０】
次に、ステップＳ４において、制御部３は、ステップＳ３で選択された光照射プローブと光検出プローブを結ぶ最短経路に位置する各々のチャンネルについて、上記選択された光照射プローブから上記選択された光検出プローブにむかう順序でステップＳ１及びステップＳ２を実行したのち、前記最短経路の各々のチャンネルについて、上記選択された光検出プローブから上記選択された光照射プローブにむかう順序でステップＳ２及びステップＳ１を実行する。
【００２１】
次に、ステップＳ５において、制御部３は、上記最短経路に位置する各々のチャンネルと、当該チャンネルに隣接するチャンネルとの間において、ステップＳ１及びステップＳ２を実行する。以下において、図４に示された第一の脳機能計測方法の具体例を詳しく説明する。
【００２２】
チャンネルkに関し各波長λにおける観測光量Jは、以下の式（１）で示される。
【００２３】
【数１】
2018190130.tif 000003
【００２４】
ここで、ri,0,λは光照射プローブと頭皮の間における光透過率の時間平均値、rj,0,λは光検出プローブと頭皮の間における光透過率の時間平均値、Rは組織透過率の時間平均値を意味する。
【００２５】
また、多チャンネルを構成するｎ本の光照射プローブとｍ本の光検出器プローブには、それぞれ、いずれか一つの図１に示された光減衰器１１〜１４が接続されているところ、ai（但し、iは１以上n以下で、aiは０以上１以下）は光照射プローブに接続された光減衰器の透過率、aｊ（jは(n+1)以上(n+m)以下で、aｊは０以上１以下）は光検出プローブに接続された光減衰器の透過率を示す。このとき、チャンネルkでのfNIRS信号に含まれた雑音の大きさhは、以下の式（２）で表される。
【００２６】
【数２】
2018190130.tif 000004
【００２７】
ここで、nj,λ(t)は光検出器プローブjの計測で生じる検出器雑音を表すが、同じ規格の検出器を用いる場合にはプローブ間での較差はほぼ無視できる。このため、式（２）を参照すると、チャンネル間でfNIRS信号に含まれた雑音の大きさhを平準化することは、各チャンネルにおける観測光量Jk,λを平準化することを意味することが分かる。
【００２８】
ところで、式（１）の両辺の対数をとると、次式（３）となる。
【００２９】
【数３】
2018190130.tif 000005
【００３０】
ただし、式（３）の係数C1と係数C2は次式（４）の関係を満たす。
【００３１】
【数４】
2018190130.tif 000006
【００３２】
なお、以下においては、式（３）の第一項におけるIi,λlairi,0,λlを実効入射光量、第二項におけるrj,0,λlajを実効検出効率と呼ぶ。
【００３３】
ここで、N個の計測チャンネルとこれを構成するプローブの関係をプローブ配置行列と呼ぶ行列Gで表現する。行列GはN×(n+m)の行列であり、その要素gは光照射プローブiと光検出器プローブjがチャンネルkを構成するとき、次式（５）で定義される。
【００３４】
【数５】
2018190130.tif 000007
【００３５】
どのような規模やパターンのチャンネル若しくはプローブの配置であっても、これに対応して行列Gは必ず一つだけ定まり、また行列Gの擬似逆行列G+もこれに対応して必ず一つ定まる。このようなプローブ配置行列を用いると、任意のチャンネルに関する式（３）の関係は次式（６）のような行列演算にまとめることができる。
【００３６】
【数６】
2018190130.tif 000008
【００３７】
ここで、sλは観測光量の対数logJk,λ（但し、kは１以上N以下）を要素とする列ベクトルであり、実測によって決められる。一方、ρλは推定したい実効入射光量に関する項log(Ii,λairi,0,λC1)（但し、iは１以上n以下）及び、実効検出効率に関する項log(rj,0,λlajC2)（但し、jは(n+1)以上(n+m)以下）を要素とする以下の式（７）で示される列ベクトルである。
【００３８】
【数７】
2018190130.tif 000009
【００３９】
ここで、式（７）の左から擬似逆行列G+を掛けることによって、ρλの特殊解が以下の式（８）のように求まる。
【００４０】
【数８】
2018190130.tif 000010
【００４１】
ここで、式（７）と式（８）を比較すると、以下の式（９）及び式（１０）の関係が得られる。
【００４２】
【数９】
2018190130.tif 000011
【００４３】
【数１０】
2018190130.tif 000012
【００４４】
このようにして各チャンネルの観測光量Jk,λが得られ、チャンネルの配置を指定する行列Gが分かっているときには、各光照射プローブの実効入射光量と各光検出プローブの実効検出効率は、それぞれebi,λ/C1、ebj,λ/C2と求めることができる。
【００４５】
ところで、観測光量を平準化する上で、観測光量のみならず実効入射光量と実効検出効率を知る必要があるが、その理由は以下の通りである。観測光量は実効入射光量と実効検出効率の積である。ここで、実効入射光量が著しく異なるチャンネルがあったとしても、それに応じた実効検出効率を持てば、見かけ上等しい観測光量を示すことが起こり得る。一方、このときの実効入射光量が一つでも安全な基準光量を超えてしまっていれば、この多チャンネル計測は安全に行われているとは言えなくなる。従って、安全な基準光量以下の光量しか生体組織に照射されていないことを確認するために実効入射光量をモニタリングすることが必要である。
【００４６】
いま、図３のような多チャンネルプローブ配置において、式（８）に基づいて最大の実効入射光量を持つ光照射プローブimaxと最大の実効検出効率を持つ光検出プローブjmaxを選び出すことができる。ここでは、光照射プローブimaxと光検出プローブjmaxを結ぶ最短の経路がν個の光照射プローブとν個の光検出プローブによって構成されているものとして、その経路上にある(2ν−1)個のチャンネルにおける観測光量の平準化について説明する。
【００４７】
なお、本経路内の光照射プローブを示す符号としてS1,S2,…,Sνを、光検出プローブを表す符号としてD1,D2,…,Dνを、各チャンネルを表す符号としてCH1,CH2,…,CHκ（但し、κ＝2ν−1）を用いる。従って、図３においては、S1が上記の光照射プローブimaxを表し、Dνが上記の光検出プローブjmaxを表している。
【００４８】
以下では、説明を平易にするため、チャンネルCH1〜CH3のみを考える。このとき、式（６）の内容を次式（１１）のように書き下すことができる。
【００４９】
【数１１】
2018190130.tif 000013
【００５０】
なお、式（１１）におけるSCH1からSCH3は逐次、計測値を算出することができる。ここで、SCH2とSCH3をSCH1の観測値/SCH1に揃えることを考える。まず、チャンネルCH1とチャンネルCH2の観測光量をモニタしながら、SCH2が上記観測値/SCH1となるように光照射プローブS2に接続された光減衰器１２を制御する。これにより、次式（１２）で示される状態を実現することができる。
【００５１】
【数１２】
2018190130.tif 000014
【００５２】
次に、チャンネルCH2とチャンネルCH3における観測光量をモニタしながら、SCH3が上記観測値/SCH1となるように光検出プローブD2に接続された光減衰器１４を制御する。このような制御により、次式（１３）で示される状態を実現することができる。
【００５３】
【数１３】
2018190130.tif 000015
【００５４】
ここで、式（１３）が成り立つための必要十分条件は、次式（１４）で示される。
【００５５】
【数１４】
2018190130.tif 000016
【００５６】
従って、上記一連の操作によってチャンネルCH1〜CH3で観測光量がe(ρS1+ρD1)に平準化されると同時に、全ての光照射プローブの実効入射光量と全ての検出プローブの実効検出効率がそれぞれeρS1、eρD1に平準化されていることが分かる。本操作は、透過率ai,ajを指定した値に制御することなく雑音平準化を達成するための中核的な手段であり、以下ではこれを隣接チャンネル平準化操作と呼ぶ。図３においては、この隣接チャンネル平準化操作をチャンネルCH4（図示していない）からチャンネルCHνに敷衍することにより、全てのチャンネルCH1〜CHνで平準化することができる。
【００５７】
このように、隣接チャンネル間における平準化の操作では、基点となるチャンネルの観測光量、またそのチャンネルを構成する光照射プローブ若しくは光検出プローブでの実効入射光量若しくは実効検出効率に合わせて隣接したチャンネルやプローブの平準化が行える。また、隣接チャンネル平準化操作は、光照射プローブ又は光検出プローブを介して連結されるチャンネルに対して敷衍してゆくことが可能である。このように連結されて構成される多チャンネルを、以下では連結チャンネル群と呼ぶ。
【００５８】
ところで、上記の操作では、一連の光検出プローブの実効検出効率はeρD1に平準化される。しかし、チャンネルCHνを構成する光検出プローブDνの実効検出効率は全ての光検出プローブ中で最大であるため、eρDνがeρD1より大きいことは自明である。そこで、観測光量を最大化し、実現しうる最大のS/N比に雑音を平準化するためには、各光検出プローブの実効検出効率をeρDνに平準化する必要がある。以下、この操作について説明する。
【００５９】
本目的を達成するためには、上記隣接チャンネル平準化操作をチャンネルCH(κ-1)まで行った時点で一旦終わらせる。この時点での状態について、煩雑を避けるためにチャンネルCH(κ-2)からチャンネルCHκ間についてのみ記述すると次式（１５）となる。
【００６０】
【数１５】
2018190130.tif 000017
【００６１】
ここで、チャンネルCH(κ-1)とチャンネルCHκの観測光量をモニタしながら、SCH(κ-1)が観測値/SCHκとなるように光検出プローブD(ν-1)に接続された光減衰器を制御する。これにより、次式（１６）で示される状態が実現される。
【００６２】
【数１６】
2018190130.tif 000018
【００６３】
式（１６）からは、このとき次式（１７）も成立していることがわかる。
【００６４】
【数１７】
2018190130.tif 000019
【００６５】
従って、チャンネルCH(κ-2)からチャンネルCHκで観測光量がe(ρS1+ρDν)に平準化されると同時に、全ての光照射プローブの実効入射光量と全ての光検出プローブの実効検出効率がそれぞれeρS1、eρDνに平準化されていることが分かる。本操作をチャンネルCH(ν-3)からチャンネルCH1に遡行してゆけば、チャンネルCHκからチャンネルCH1までの全チャンネルを観測光量e(ρS1+ρDν)に平準化できる。
【００６６】
以上のような手順により、任意の連結チャンネル群内の最大の実効入射光量を持つ光照射プローブimaxと最大の実効検出効率を持つ光検出プローブjmax間の観測光量を平準化することができ、さらに群内の全チャンネルについて隣接チャンネル平準化操作によって、観測光量の平準化を図ることができる。以上を踏まえて連結チャンネル群内を平準化する手順を整理すると、以下のようになる。
【００６７】
第一に、プローブを被検体の頭部に装着した後に、照射側と検出側の全ての光減衰器の透過率を例えば約５０％に設定し、全ての光源出力を生体（被検体）に対して安全な光照射の最大光量Isafeになるように設定する。
【００６８】
第二に、波長λを選び、式（９）及び式（１０）により各プローブのebi,λとebj,λを推定し、ebi,λが最大値となる光照射プローブimaxとebj,λが最大値となる光検出プローブjmaxを選び出す。
【００６９】
第三に、光照射プローブimaxと光検出プローブjmaxを結ぶ最短の経路上にあるチャンネルで隣接チャンネル平準化操作を行い、観測光量をe(bimax,λ+bjmax,λ)に揃える。
【００７０】
第四に、第三の手順で平準化されたチャンネルに隣接するチャンネルとの間でも同様に隣接チャンネル平準化操作を行う。
【００７１】
第五に、第二から第四の手順を全ての使用波長について行う。そして、第六に、第五までの手順により決められた設定を保持した状態でfNIRS計測を行う。さらに、第七として、本計測終了と同時に第一の手順に戻るか、若しくは全ての光源出力を遮断する。
【００７２】
［雑音分散の最小化］
雑音分散をなるべく低い水準に合わせて平準化するためには、計測される連結群の少なくとも一つのチャンネルを、図５に示されるように、前額部などの無毛頭部に配置するための光照射プローブSa及び光検出プローブDaと、光照射プローブS及び光検出プローブDを配設することが望ましい。その理由は以下の通りである。
【００７３】
無毛部のチャンネルにおける実効入射光量と実効検出効率は有毛部と比べて著しく大きく時間的にも恒常的であるため、当該チャンネルの雑音分散は著しく小さい。従って、無毛部を含んでプローブ装着を行った上でそのチャンネルを基点として上記平準化を行えば、安全な光照射強度の範囲内で最も良好な水準の信号対雑音比に全チャンネルを平準化することができる。
【００７４】
次に、異なる日時にfNIRS計測を行って得られたデータや異なる対象者について得られたデータを統計的に比較したい場合に有用な雑音分散標準化手法について説明する。
【００７５】
本手法は、独立した二つの連結チャンネル群（非連結二群）の雑音分散を平準化することに帰結するが、上記の連結チャンネル群内の平準化手順では実現することはできない。そこで以下では、時間的にのみ非連結な二群間の雑音平準化と、時間的かつ空間的に非連結な二群間の雑音平準化の方法について説明する。
【００７６】
［時間的にのみ非連結なチャンネル群間の雑音平準化］
時間的にのみ非連結な二群とは、例えば同一の対象者や、同一の部位に対して異なる日時にfNIRS計測を行う場合をいう。以下においては、このときの二群の雑音平準化の方法を、図１に示された制御部３による制御により実現する場合を例として、図６を参照しつつ詳しく説明する。なお、図６に示された脳機能計測方法は、図１に示された制御部３による制御により実現する場合に限られないことは言うまでもない。
【００７７】
図６に示された脳機能計測方法では、前提として、図５に示されるように、被検体の頭部へホルダ２０を装着することにより、前額部などの無毛頭部の形態上位置を同定しやすい部位にチャンネルを配置するための光照射プローブSa及び光検出プローブDaと、光照射プローブS及び光検出プローブDを配設する。
【００７８】
そして、最初にステップＳ６では、制御部３は、図４に示されたステップＳ１及びステップＳ２を実行した後に、ステップＳ１で得られた観測値を基準値（基準観測光量）として記憶部９に保存する。
【００７９】
次に、ステップＳ７では、制御部３は、前額部などの無毛頭部における基準となる計測において、少なくとも一つの光検出プローブに接続された光減衰器の機能をオフすることにより、当該光検出プローブへ入射する光の透過率を百パーセントとし、少なくとも一つの光検出プローブが計測するチャンネルにつき観測された最大光量値（最大観測光量）を記憶部９に保存する。
【００８０】
次に、ステップＳ８では、制御部３は、平準化を行いたい独立したチャンネル群のうち前額部などの無毛頭部における少なくとも一つの光検出プローブに接続された光減衰器の機能をオフにした状態で、当該光検出プローブへ入射する光の透過率を百パーセントとした上で、平準化対象とする当該チャンネルにつき観測される光量が上記最大光量値に一致するように、本チャンネルにつき計測するための光照射プローブに接続された光減衰器を調整することにより、本光照射プローブから照射される光量を調節する。
【００８１】
次に、ステップＳ９では、制御部３は、ステップＳ８を実行した後に、平準化対象とするチャンネルにつき観測される光量の値が上記基準値に一致するように、本チャンネルにつき計測するための光検出プローブに接続された光減衰器を調整することにより、本光検出プローブで検出される光量を調節する。なお、この段階で、本チャンネルの実効入射光量と実効検出効率はそれぞれ、平準化基準となるチャンネル群における実効入射光量及び実効検出効率の値と一致することになる。
【００８２】
次に、ステップＳ１０では、制御部３は、平準化対象とした上記チャンネルを、図４に示された第一の光検出プローブと第一の光照射プローブとの間のチャンネルに相当するものとして、図４に示されたステップＳ１からステップＳ５を実行する。なお、本ステップにより、当該チャンネル群の雑音分散は平準化基準をなすチャンネル群の雑音分散と一致することになる。
【００８３】
以下においては、図６に示された脳機能計測方法が、非連結二群の平準化を実現することを説明する。いま、図６に示されたステップＳ６及びステップＳ７で得られる基準観測光量と最大観測光量をそれぞれJ(1)、J(2)とし、log J(1)をS(1)、log J(2)をS(2)とする。このとき、次の式（１８）及び式（１９）が成立する。
【００８４】
【数１８】
2018190130.tif 000020
【００８５】
【数１９】
2018190130.tif 000021
【００８６】
ここで、式（１９）のρj0は光検出プローブに接続された光減衰器の機能をオフにして光透過率を百パーセントとしたときの、最大観測光量の計測時における実効検出効率の対数である。いま、基準となる連結群とは独立した別の連結群の無毛部に装着されたチャンネルをCHとし、これを構成する光照射プローブをs、光検出プローブをdとすると、次式（２０）が成立する。
【００８７】
【数２０】
2018190130.tif 000022
【００８８】
ここで、ステップＳ７に従って、チャンネルCHを計測するための光検出プローブdに接続された光減衰器の機能をオフにすると、段落［００７３］で説明した無毛頭部での観測の恒常性から式（２０）のρdはρj0に一致するため、次式（２１）が成立する。
【００８９】
【数２１】
2018190130.tif 000023
【００９０】
さらに、ステップＳ８で説明したように、SCHがS(2)となるように光照射プローブに接続された光減衰器を調整することにより、式（１９）と式（２１）の比較からρsをρimaxとすることができることが分かる。このとき、次式（２２）が成立する。
【００９１】
【数２２】
2018190130.tif 000024
【００９２】
さらに、ステップＳ９で説明したように、SCHがS(1)となるように光検出プローブに接続された光減衰器を調整することにより、式（１８）と式（２２）との比較からρj0をρjmaxとすることができることが分かる。このとき、次式（２３）が成立する。
【００９３】
【数２３】
2018190130.tif 000025
【００９４】
そして、本チャンネルを基点として図４に示された方法により平準化を行うことにより、一つの連結群内のチャンネルを、基準となる連結群の雑音分散に平準化することができる。
【００９５】
［時間的かつ空間的に非連結なチャンネル群間の雑音平準化］
運動関連領野などが位置する頭頂部の計測を行う場合には、上記のように前額部等に基点となるチャンネルを配置すると、計測したい位置から離れているため、しばしば非効率となる。一方、空間的に共通基点を持たない場合には、時間的のみならず空間的にも当該二群は非連結になる。
【００９６】
このようなことから、以下では、時間的かつ空間的に非連結の二群の雑音平準化を行う方法を、図１に示された制御部３による制御により実現する場合を例として、図７を参照しつつ詳しく説明する。
【００９７】
なお、図７に示された脳機能計測方法は、図１に示された制御部３による制御により実現する場合に限られないことは言うまでもない。
【００９８】
図７に示された脳機能計測方法は、前提として、図８に示されるように、連結群内の特定の隣接した二つの基準チャンネルを構成する３つのプローブ、すなわち基準照射プローブSS1,SS2及び基準検出プローブSDに対してのみ指定値制御可能な光減衰器を接続し、
これらの光減衰器の透過率が指定値になるよう、制御部３が各々の光減衰器を調整する。
【００９９】
なお、図９は、図８に示された光照射プローブ及び光検出プローブの配置を二次元平面に記した模式図であるが、それぞれ図９に示された光照射プローブS1が基準照射プローブSS1に対応し、光照射プローブS2が基準照射プローブSS2に対応する。また、図９に示された光検出プローブD1が基準検出プローブSDに対応する。また、図９に示されたチャンネルCH1,CH2が、上記二つの基準チャンネルに対応する。
【０１００】
最初に、ステップＳ１６において、制御部３は、基準チャンネルSCH1,SCH2を含む連結多チャンネル計測につき、図４に示された方法を適用して雑音平準化を実行し、これにより得られた実効検出効率と実効入射光量の値を記憶部９に保存する。
【０１０１】
次に、別の日時や別の位置で連結多チャンネル計測を行った際、ステップＳ１７において、制御部３は、基準検出プローブSDに接続された光減衰器を制御することにより、基準検出プローブSDにおける検出効率をステップＳ１６で保存された実効検出効率とし、基準照射プローブSS1,SS2に接続された光減衰器を制御することにより、基準照射プローブSS1,SS2における照射光量をステップＳ１６で保存された実効入射光量とする。
【０１０２】
次に、ステップＳ１８において、制御部３は、ステップＳ１７の後に、基準チャンネルSCH1,SCH2を基点として再度図４に示されたステップＳ１からステップＳ５、すなわち隣接チャンネル平準化操作を実行し、連結チャンネル群内の全てのチャンネルについて雑音平準化を行う。
【０１０３】
以下において、図９を参照しつつ、上記ステップＳ１７における制御方法を具体的に説明する。図９に示されたチャンネルCH1,CH2、すなわち光照射プローブS1,S2及び光検出プローブD1に関して次式（２４）が成立する。
【０１０４】
【数２４】
2018190130.tif 000026
【０１０５】
ここで、式（２４）にプローブ配置行列の一般化逆行列を演算することにより、次式（２５）が得られる。
【０１０６】
【数２５】
2018190130.tif 000027
【０１０７】
ところで、平準化を終えている基準となる連結群の各チャンネルでは、式（２３）が成立している。また、上記のような指定値制御が可能な素子を用いている場合には、ρxを実現している当該プローブxに接続された光減衰器の透過率axを具体的に式（７）により得ることができる。これらのことから、ρS1をρimaxとし、かつρD1をρjmaxとするためには、透過率aS1をe(ρimax-ρS1)倍にし、透過率aD1をe(ρjmax-ρD1)倍に変更すればよいことが分かる。
【０１０８】
また、式（２４）及び図９では、上記の基準検出プローブ及び基準照射プローブとして、一つの光検出プローブD1とこれを挟む二つの光照射プローブS1,S2が用いられているが、これらの替わりに、一つの光照射プローブとこれを挟む二つの光検出プローブを用いた場合にも、同じ方法により当該連結群の平準化が可能である。
【０１０９】
以上のような本発明の実施の形態に係る脳機能計測装置及び本装置を含む脳機能計測システムや脳機能計測方法によれば、当該装置の較正作業を不要にし、若しくは最小限に留めることができる。すなわち、光減衰器の透過率を指定値に調整せずに、連結チャンネル群内の全てのチャンネルの雑音を平準化することができる。
【０１１０】
より具体的には、図４及び図６に示された方法によれば、当該光減衰器の較正は全く不要となる。一方、図７に示された方法においては、３つのプローブで指定値制御が可能な光減衰器を用いる必要があるが、この場合にもいかにチャンネル数が大規模化しても、指定値制御が可能な光減衰器は上記３つのプローブに接続される３個だけで足りることになる。従って、脳機能計測装置や本装置を含むシステムの調整メンテナンスを大幅に簡略化することができ、小型化も可能になる。
【０１１１】
また、他面においては、本発明の実施の形態に係る脳機能計測装置及び本装置を含む脳機能計測システムや脳機能計測方法によれば、複数の独立した連結群の間の雑音分散を平準化することができる。その方法には、図６に示されたように、時間的に雑音分散が不変と想定できる基点を群間で共通に設定する方法と、図７に示されたように、一つの群内の隣接した二つのチャンネルを構成するプローブにのみ指定値制御が可能な光減衰器を用いる方法がある。
【０１１２】
これら双方の方法は、一人の被験者におけるデータの経時的変化の統計的比較を可能にするものである。加えて、後者の方法は、複数の被験者や複数の頭部部位での計測に対しての統計的比較も可能にするものである。
【０１１３】
さらに、他面においては、本発明の実施の形態に係る脳機能計測装置及び本装置を含む脳機能計測システムや脳機能計測方法によれば、無毛部に置かれたチャンネルを平準化の基点とすることにより、連結化が可能なあらゆる部位や日時の計測における雑音分散を無毛部の水準に合わせて平準化することができる。これにより、安全な光照射強度の範囲を遵守した上で達成しうる最大限の信号対雑音比を定常的に実現することができる。
【符号の説明】
【０１１４】
１  脳機能計測装置、３  制御部、５  光源部、７  データ計測部、９  記憶部、１１〜１４  光減衰器、Ｓ，Ｓ１，Ｓ２  光照射プローブ、Ｄ，Ｄ１〜Ｄ３  光検出プローブ、ＳＳ１，ＳＳ２  基準照射プローブ、ＳＤ  基準検出プローブ。
【図１】
2018190130.tif 000028
【図２】
2018190130.tif 000029
【図３】
2018190130.tif 000030
【図４】
2018190130.tif 000031
【図５】
2018190130.tif 000032
【図６】
2018190130.tif 000033
【図７】
2018190130.tif 000034
【図８】
2018190130.tif 000035
【図９】
2018190130.tif 000036
【国際調査報告】
2018190130.tif 000037
2018190130.tif 000038