(57)【要約】【課題】ｆＭＲＩを用いた脳研究に対する新たな可能性を切り開くとともに、究極的には脳機能のより正確な把握を可能にする。【解決手段】脳神経活動に一次的に関連づけられる神経活動信号の時系列データと、ＭＲＩ装置からのＭＲ信号の変動データとの相関を分析し、その結果に基づいて脳機能を測定する脳機能測定装置１において、前記神経活動信号の時系列データから、予め定めたｏｘｙＨｂレスポンス関数及びｄｅｏｘｙＨｂレスポンス関数を介してｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度の時系列データの時系列データをそれぞれ算出するｏｘｙＨｂ濃度算出部２及びｄｅｏｘｙＨｂ算出部３と、前記ｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度の時系列データを少なくともパラメータとして、前記変動データとの相関を比較し得る比較データを算出する比較データ算出部４とを設けた。【選択図】図２
2005006774.tif 000002
【特許請求の範囲】
【請求項１】
脳神経活動に一次的に関連づけられる信号である神経活動信号の時系列データと、ＭＲＩ装置から得られるＭＲ信号の変動データとの相関を分析し、その結果に基づいて脳機能を測定するものであって、
前記神経活動信号の時系列データから、予め定めたｏｘｙＨｂレスポンス関数を介してｏｘｙＨｂ濃度の時系列データを算出するｏｘｙＨｂ濃度算出部と、
前記脳神経活動の時系列データから、予め定めたｄｅｏｘｙＨｂレスポンス関数を介してｄｅｏｘｙＨｂ濃度の時系列データを算出するｄｅｏｘｙＨｂ算出部と、
前記ｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度の時系列データを少なくともパラメータとして、前記ＭＲ信号の変動データとの相関を比較し得る比較データを算出する比較データ算出部とを備えていることを特徴とする脳機能測定装置。
【請求項２】
前記ｏｘｙＨｂレスポンス関数が、神経活動信号のｏｘｙＨｂ濃度に対する空間変動影響を表す空間変動関数をさらに含むものであり、
前記ｄｅｏｘｙＨｂレスポンス関数が、神経活動信号のｄｅｏｘｙＨｂ濃度に対する空間変動影響を表す空間変動関数をさらに含むものである請求項１記載の脳機能測定装置。
【請求項３】
脳神経活動に一次的に関連づけられる信号である神経活動信号の時系列データと、ＭＲＩ装置から得られるＭＲ信号の変動データとの相関を分析し、その結果に基づいて脳機能を測定する脳機能測定装置に適用されるものであって、
前記神経活動信号の時系列データから、予め定めたｏｘｙＨｂレスポンス関数を介してｏｘｙＨｂ濃度の時系列データを算出するｏｘｙＨｂ濃度算出部と、
前記脳神経活動の時系列データから、予め定めたｄｅｏｘｙＨｂレスポンス関数を介してｄｅｏｘｙＨｂ濃度の時系列データを算出するｄｅｏｘｙＨｂ算出部と、
前記ｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度の時系列データを少なくともパラメータとして、前記ＭＲ信号の変動データとの相関を比較し得る比較データを算出する比較データ算出部としての機能を前記脳機能測定装置に備えさせることを特徴とする脳機能測定プログラム。
【請求項４】
脳神経活動に一次的に関連づけられる信号である神経活動信号の時系列データと、ＭＲＩ装置から得られるＭＲ信号の変動データとの相関に基づいて脳の機能部位を特定するものであって、
前記神経活動信号の時系列データから、予め定めた第１ｏｘｙＨｂレスポンス関数を介して、光吸収変化から測定し得るｏｘｙＨｂ濃度の時系列データを算出するとともに、予め定めた第２ｏｘｙＨｂレスポンス関数を介して、光散乱変化から測定し得るｏｘｙＨｂ濃度の時系列データを算出するｏｘｙＨｂ濃度算出部と、
前記神経活動信号の時系列データから、予め定めた第１ｄｅｏｘｙＨｂレスポンス関数を介して、光吸収変化から測定し得るｄｅｏｘｙＨｂ濃度の時系列データを算出するとともに、予め定めた第２ｄｅｏｘｙＨｂレスポンス関数を介して、光散乱変化から測定し得るｄｅｏｘｙＨｂ濃度の時系列データを算出するｄｅｏｘｙＨｂ濃度算出部と、
前記ｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度の各時系列データをパラメータとして、前記ＭＲ信号の変動データとの相関を比較し得る比較データを算出する比較データ算出部とを備えていることを特徴とする脳機能測定装置。
【請求項５】
前記各ｏｘｙＨｂレスポンス関数が、神経活動信号のｏｘｙＨｂ濃度に対する空間変動影響を表す空間変動関数をさらに含むものであり、
前記各ｄｅｏｘｙＨｂレスポンス関数が、神経活動信号のｄｅｏｘｙＨｂ濃度に対する空間変動影響を表す空間変動関数をさらに含むものである請求項４記載の脳機能測定装置。
【請求項６】
脳神経活動に一次的に関連づけられる信号である神経活動信号の時系列データと、ＭＲＩ装置から得られるＭＲ信号の変動データとの相関を分析し、その結果に基づいて脳機能を測定する脳機能測定装置に適用されるものであって、
前記神経活動信号の時系列データから、予め定めた第１ｏｘｙＨｂレスポンス関数を介して、光吸収変化から測定し得るｏｘｙＨｂ濃度の時系列データを算出するとともに、予め定めた第２ｏｘｙＨｂレスポンス関数を介して、光散乱変化から測定し得るｏｘｙＨｂ濃度の時系列データを算出するｏｘｙＨｂ濃度算出部と、
前記神経活動信号の時系列データから、予め定めた第１ｄｅｏｘｙＨｂレスポンス関数を介して、光吸収変化から測定し得るｄｅｏｘｙＨｂ濃度の時系列データを算出するとともに、予め定めた第２ｄｅｏｘｙＨｂレスポンス関数を介して、光散乱変化から測定し得るｄｅｏｘｙＨｂ濃度の時系列データを算出するｄｅｏｘｙＨｂ濃度算出部と、
前記ｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度の各時系列データをパラメータとして、前記ＭＲ信号の変動データとの相関を比較し得る比較データを算出する比較データ算出部とを前記脳機能測定装置に備えさせることを特徴とする脳機能測定プログラム。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、血行動態変動と脳神経活動との相関から脳の活動部位を測定する脳機能測定装置等に関するものである。
【０００２】
【従来の技術】
脳機能測定（機能的磁気共鳴イメージング、以下ｆＭＲＩという）で得られる情報は、脳の神経活動に付随して起こるｏｘｙＨｂの増加を反映した二次的なものであり、ＨＲＦ（ｈｅｍｏｄｙｎａｍｉｃ  ｒｅｓｐｏｎｓｅ  ｆｕｎｃｔｉｏｎ）と称される時間項を含む一定の波形を有した関数によって神経活動と関係づけられている（非特許文献１、２参照）。実際の解析例を説明すると、まずｆＭＲＩの元の信号であるＴ２＊信号の変動とＨＲＦとの畳み込み積分を計算することによってＨＲＦに相当する信号を強調して取り出す時間的な処理を行い、その後、ランダムフィールド理論による空間的な処理、及び、統計検定を経て、ｆＭＲＩの脳機能画像を得ている。
【０００３】
【非特許文献１】
Ｂｕｃｈｅｌ，  Ｃ．，  Ｈｏｌｍｅｓ，  Ａ．  Ｐ．，  Ｒｅｅｓ，  Ｇ．  ａｎｄ  Ｆｒｉｓｔｏｎ，  Ｋ．  Ｊ．  （１９９８）．  Ｃｈａｒａｃｔｅｒｉｚｉｎｇ  ｓｔｉｍｕｌｕｓ−ｒｅｓｐｏｎｓｅ  ｆｕｎｃｔｉｏｎｓ  ｕｓｉｎｇ  ｎｏｎｌｉｎｅａｒ  ｒｅｇｒｅｓｓｏｒｓ  ｉｎ  ｐａｒａｍｅｔｒｉｃ  ｆＭＲＩ  ｅｘｐｅｒｉｍｅｎｔｓ．  Ｎｅｕｒｏｉｍａｇｅ  ８，  １４０−１４８．
【非特許文献２】
Ｃｌａｒｋ，  Ｖ．  Ｐ．  （２００２）．  Ｏｒｔｈｏｇｏｎａｌ  Ｐｏｌｙｎｏｍｉａｌ  Ｒｅｇｒｅｓｓｉｏｎ  ｆｏｒ  ｔｈｅ  Ｄｅｔｅｃｔｉｏｎ  ｏｆ  Ｒｅｓｐｏｎｓｅ  Ｖａｒｉａｂｉｌｉｔｙ  ｉｎ  Ｅｖｅｎｔ−Ｒｅｌａｔｅｄ  ｆＭＲＩ．  ＮｅｕｒｏＩｍａｇｅ  １７，  ３４４−３６３．
【０００４】
【発明が解決しようとする課題】
ところが、前記ＨＲＦは（物理量であるＴ２＊の変動）と（生体の神経活動電位）との比であり、その単位が生理学的に明確な意味を持たない。加えて脳神経活動であるニューロン発火から血行動態が変動するまでの過程は完全に解明されているわけではなく、前記ＨＲＦという時間関数のみでニューロン発火と血行動態変動とを関連づけることが適切であるかどうかは疑問である。事実、「ＨＲＦは脳の部位等により異なる」という結果が近年報告されており、前記ＨＲＦという一定の波形を有した時間関数による脳機能測定では、刺激や部位によっては正確な脳機能測定を為し得ないおそれがあることは十分に考えられる。
【０００５】
また、今後、ニューロン発火から血行動態変動までの過程を解明していくにあたって、前述したように生理学的に明確な意味を持たないＨＲＦという単一関数では、前記Ｔ２＊信号の変動の原因となるｄｅｏｘｙＨｂの変動が、果たして血流の変動によるものなのか、あるいはｄｅｏｘｙＨｂ濃度の変動によるものなのか等を深く追求していくことはできず、このままではｆＭＲＩを用いた脳研究のさらなる発展を阻害することにもなりかねない。
【０００６】
そこで、本発明は、神経活動とｆＭＲＩで得られる情報との関係を、種々のパラメータからなりしかもそれらがより生理学的に意味のある関数、すなわち少なくともｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度という血液状態の関数を介して示すことによって前記問題点を解決し、脳研究に対する新たな可能性を切り開くとともに、究極的には脳機能のより正確な把握を可能にすべく図ったものである。
【０００７】
【課題を解決するための手段】
すなわち、図１に示すように、本発明にかかる脳機能測定装置１は、脳神経活動に一次的に関連づけられる信号である神経活動信号の時系列データと、ＭＲＩ装置から得られるＭＲ信号の変動データとの相関を分析し、その結果に基づいて脳機能を測定するものであって、
前記神経活動信号の時系列データから、予め定めたｏｘｙＨｂレスポンス関数を介してｏｘｙＨｂ濃度の時系列データを算出するｏｘｙＨｂ濃度算出部２と、
前記脳神経活動の時系列データから、予め定めたｄｅｏｘｙＨｂレスポンス関数を介してｄｅｏｘｙＨｂ濃度の時系列データを算出するｄｅｏｘｙＨｂ算出部３と、
前記ｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度の時系列データを少なくともパラメータとして、前記ＭＲ信号の変動データとの相関を比較し得る比較データを算出する比較データ算出部４とを備えていることを特徴とする。
【０００８】
このようなものであれば、ニューロン発火という神経活動から、ｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度が時間的かつ空間的に変動し、その結果Ｔ２＊信号が変動するという、より生理学的な過程に近いと考えられる仮定の下、神経活動信号とＴ２＊信号の変動データとの相関を判断するようにしていることから、従来のＨＲＦという一定の時間関数による脳機能測定では全く測定し得なかった刺激や部位等に対して、新たな測定結果をもたらし得る。
【０００９】
またこのように多くのパラメータによって神経活動信号とＴ２＊信号の変動データとを結びつけているため、今後行われるであろう種々の実験や研究結果に対応できる発展性の高い脳機能測定装置１を提供できる。
【００１０】
さらにその演算過程において、ｏｘｙＨｂ濃度とｄｅｏｘｙＨｂ濃度との変動が算出されるため、そのことから逆に、ニューロン発火から血行動態変動までの過程を解明していく研究の一助にもなりうる。
【００１１】
血液状態の時間変動に加えて、空間変動の関数も組み込むことで、より正確な把握を可能にするには、前記ｏｘｙＨｂレスポンス関数が、神経活動信号のｏｘｙＨｂ濃度に対する空間変動影響を表す空間変動関数をさらに含むものであり、前記ｄｅｏｘｙＨｂレスポンス関数が、神経活動信号のｄｅｏｘｙＨｂ濃度に対する空間変動影響を表す空間変動関数をさらに含むものであることが望ましい。
【００１２】
本装置１による脳機能測定において必要な予め定めておくべき各レスポンス関数を簡単な方法で特定するには、図２に示すように、前記神経活動信号の時系列データから、予め定めた第１ｏｘｙＨｂレスポンス関数を介して、光吸収変化から測定し得るｏｘｙＨｂ濃度の時系列データを算出するとともに、予め定めた第２ｏｘｙＨｂレスポンス関数を介して、光散乱変化から測定し得るｏｘｙＨｂ濃度の時系列データを算出するｏｘｙＨｂ濃度算出部２と、前記神経活動信号の時系列データから、予め定めた第１ｄｅｏｘｙＨｂレスポンス関数を介して、光吸収変化から測定し得るｄｅｏｘｙＨｂ濃度の時系列データを算出するとともに、予め定めた第２ｄｅｏｘｙＨｂレスポンス関数を介して、光散乱変化から測定し得るｄｅｏｘｙＨｂ濃度の時系列データを算出するｄｅｏｘｙＨｂ濃度算出部３と、前記ｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度の各時系列データをパラメータとして、前記ＭＲ信号の変動データとの相関を比較し得る比較データを算出する比較データ算出部４とを備えているものが好ましい。
【００１３】
このようなものであれば、各レスポンス関数を特定するにあたって、近赤外分光法（ＮＩＲＳ）で実際に測定したｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度と、演算過程で算出されるｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度とが一致するように設定すればよいため、その信頼性を向上させることができる。また近赤外分光法（ＮＩＲＳ）は測定が容易であるという利点を有するため、本装置１と組み合わせたシステムとし、被験者毎に前記各レスポンス関数を設定し、より正確な脳機能測定を行うようにすることも無理なくできる。
【００１４】
【発明の実施の形態】
以下に本発明の一実施形態について図面を参照して説明する。
【００１５】
本実施形態に係る脳機能測定装置１は、図３に示すように、脳神経活動に一次的に関連づけられる信号である神経活動信号の時系列データと、ＭＲＩ装置から得られるＭＲ信号であるＴ２＊信号の変動データとの相関をｏｘｙＨｂ濃度とｄｅｏｘｙＨｂ濃度の関数を介して分析し、その結果に基づいて脳機能を測定するものであり、より具体的には、前記神経活動信号の時系列データからｏｘｙＨｂ濃度に関する時系列データを算出するｏｘｙＨｂ濃度算出部２と、前記脳神経活動の時系列データからｄｅｏｘｙＨｂ濃度に関する時系列データを算出するｄｅｏｘｙＨｂ濃度算出部３と、前記ｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度に関する時系列データをパラメータとして、前記Ｔ２＊信号の変動データとの相関を比較し得る比較データを算出する比較データ算出部４とを備えてなる。
【００１６】
なお、この脳機能測定装置１は、ハードウェアとしてＣＰＵや記憶装置、入出力装置等を備えており、前記記憶装置に記憶させた所定のプログラムにしたがって、前記ＣＰＵ等が動作することにより、前記各部としての機能を発揮するように構成している。
【００１７】
しかして前記ｏｘｙＨｂ濃度算出部２では、前記神経活動信号の時系列データと、予め定めたｏｘｙＨｂレスポンス関数との畳み込み積分とを行うことにより、ｏｘｙＨｂ濃度が算出されるようにしてある。
【００１８】
また、前記ｄｅｏｘｙＨｂ濃度算出部３では、ｏｘｙＨｂ濃度算出部２同様、前記神経活動信号の時系列データと、予め定めたｄｅｏｘｙＨｂレスポンス関数との畳み込み積分とを行うことにより、ｄｅｏｘｙＨｂ濃度が算出されるようにしてある。
【００１９】
より詳細に説明すると、神経活動信号の時系列データとは、例えば時間ｔと位置座標及び向きｒ（三次元座標と三次元ベクトル）とをパラメータとして表されるダイポール（ニューロン発火が生じる電流源）強度Ｄｉｐｏｌｅ（ｒ，ｔ）のことであり、ＭＥＧ等から所定の演算により求めることができる。
【００２０】
一方、この実施形態では、前記ｏｘｙＨｂレスポンス関数として、第１ｏｘｙＨｂレスポンス関数と第２ｏｘｙＨｂレスポンス関数とを用意しており、また前記ｄｅｏｘｙＨｂレスポンス関数として、第１ｄｅｏｘｙＨｂレスポンス関数と、第２ｄｅｏｘｙＨｂレスポンス関数とを用意している。
【００２１】
第１ｏｘｙＨｂレスポンス関数及び第１ｄｅｏｘｙＨｂレスポンス関数は、いずれも、近赤外分光法（ＮＩＲＳ）等を利用して測定した場合の光吸収度から求められるｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度に対応するものであり、第２ｏｘｙＨｂレスポンス関数及び第２ｄｅｏｘｙＨｂレスポンス関数は、近赤外分光法（ＮＩＲＳ）、パルス光を用いた時間分解計測法、強度変調光を用いた周波数領域の光計測法等を利用して測定した場合の光散乱変化から求められるｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度に対応するものである。
【００２２】
ここではまず第１ｏｘｙＨｂレスポンス関数及び第１ｄｅｏｘｙＨｂレスポンス関数について説明する。
【００２３】
第１ｏｘｙＨｂレスポンス関数は、
ＨＲＦ＿ａｂｓ＿ｎ２ｏｘｙＨｂ（ｔ）  ｘ  Ｖｏｌ＿ａｂｓ＿ｎ２ｏｘｙＨｂ（ｔ）ｔ・・・（１）
で表されるもので、予め所定の記憶領域に例えば数値列として格納してある。もちろん演算式として格納しておいてもよい。
ここで、ＨＲＦ＿ａｂｓ＿ｎ２ｏｘｙＨｂ（ｔ）  は、神経活動と、前記光吸収変化から求め得るｏｘｙＨｂ濃度（ｍｏｌ／ｌ）の変動との関係を表す時間変動関数である。
またＶｏｌ＿ａｂｓ＿ｎ２ｏｘｙＨｂ（ｔ）  は、神経活動と、前記光吸収変化に関するｏｘｙＨｂの体積（ｍｍ３）の時間変動との関係、すなわち神経活動によってｏｘｙＨｂ濃度が空間的にどのように影響を受けるかを表す空間変動関数である。
なおＶｏｌ＿ａｂｓ＿ｎ２ｏｘｙＨｂ（ｔ）ｔ  は、Ｖｏｌ＿ａｂｓ＿ｎ２ｏｘｙＨｂ（ｔ）の転置行列を示す。
【００２４】
同様に第１ｄｅｏｘｙＨｂレスポンス関数は、
ＨＲＦ＿ａｂｓ＿ｎ２ｄｅｏｘｙＨｂ（ｔ）  ｘ  Ｖｏｌ＿ａｂｓ＿ｎ２ｄｅｏｘｙＨｂ（ｔ）ｔ・・・（２）
で表されるものである。各関数は前記第１ｏｘｙＨｂレスポンス関数のものに準ずるため説明は省略する。
【００２５】
第２ｏｘｙＨｂレスポンス関数及び第２ｄｅｏｘｙＨｂレスポンス関数については、前記式（１）及び（２）におけるａｂｓをｓｃａｔに置き換えて同様に表されるものであるため、説明を省略する。
【００２６】
このような各レスポンス関数が規定されている条件の下、この実施形態での前記ｏｘｙＨｂ濃度算出部２は、光吸収に関するｏｘｙＨｂ濃度の時系列データｏｘｙＨｂ＿ａｂｓ（ｔ）を、
ｏｘｙＨｂ＿ａｂｓ（ｔ）  ＝
Ｄｉｐｏｌｅ（ｒ，ｔ）  ｃｏｎｖ  ｛  ＨＲＦ＿ａｂｓ＿ｎ２ｏｘｙＨｂ（ｔ）  ｘ  Ｖｏｌ＿ａｂｓ＿ｎ２ｏｘｙＨｂ（ｔ）ｔ  ｝・・・（３）
なる式から算出するとともに、
光散乱に関するｏｘｙＨｂ濃度の時系列データｏｘｙＨｂ＿ｓｃａｔ（ｔ）を
ｏｘｙＨｂ＿ｓｃａｔ（ｔ）  ＝
Ｄｉｐｏｌｅ（ｒ，ｔ）  ｃｏｎｖ  ｛  ＨＲＦ＿ｓｃａｔ＿ｎ２ｏｘｙＨｂ（ｔ）  ｘ  Ｖｏｌ＿ｓｃａｔ＿ｎ２ｏｘｙＨｂ（ｔ）ｔ  ｝・・・（４）
なる式から算出する。
【００２７】
また、前記ｄｅｏｘｙＨｂ濃度算出部３は、光吸収に関するｄｅｏｘｙＨｂ濃度の時系列データｄｅｏｘｙＨｂ＿ａｂｓ（ｔ）を、
ｄｅｏｘｙＨｂ＿ａｂｓ（ｔ）  ＝
Ｄｉｐｏｌｅ（ｒ，ｔ）  ｃｏｎｖ  ｛  ＨＲＦ＿ａｂｓ＿ｎ２ｄｅｏｘｙＨｂ（ｔ）  ｘ  Ｖｏｌ＿ａｂｓ＿ｎ２ｄｅｏｘｙＨｂ（ｔ）ｔ  ｝・・・（５）
なる式から算出するとともに、
光散乱に関するｄｅｏｘｙＨｂ濃度の時系列データｄｅｏｘｙＨｂ＿ｓｃａｔ（ｔ）を
ｄｅｏｘｙＨｂ＿ｓｃａｔ（ｔ）  ＝
Ｄｉｐｏｌｅ（ｒ，ｔ）  ｃｏｎｖ  ｛  ＨＲＦ＿ｓｃａｔ＿ｎ２ｄｅｏｘｙＨｂ（ｔ）  ｘ  Ｖｏｌ＿ｓｃａｔ＿ｎ２ｄｅｏｘｙＨｂ（ｔ）ｔ  ｝・・・（６）
なる式から算出する。
ここでｃｏｎｖは畳み込み積分を表す。
【００２８】
一方、前記比較データ算出部４は、前記ｏｘｙＨｂ濃度算出部２及びｄｅｏｘｙＨｂ濃度算出部３において、式（３）〜式（６）から求めた各時系列データｏｘｙＨｂ＿ａｂｓ（ｔ）、ｄｅｏｘｙＨｂ＿ａｂｓ（ｔ）、ｏｘｙＨｂ＿ｓｃａｔ（ｔ）、ｄｅｏｘｙＨｂ＿ｓｃａｔ（ｔ）の線形和により、前記Ｔ２＊信号の変動データとの相関を比較し得る比較データＴ２＊＿ＣＡＬを算出する。
【００２９】
具体的には、比較データＴ２＊＿ＣＡＬを
2005006774.tif 000003
なる式から算出する。
ここでＡ〜Ｄは係数であり、その値は例えばＴ２＊を目標値とする重回帰演算を施すことにより設定する。
【００３０】
このように本実施形態によれば、ニューロン発火という神経活動から、ｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度が時間的かつ空間的に変動し、その結果Ｔ２＊信号が変動するという、より生理学的な過程に近いと考えられる仮定の下、神経活動信号とＴ２＊信号の変動データとの相関を、ｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度を介在させて判断するようにしていることから、従来のＨＲＦという一定の波形を有した時間関数による脳機能測定では全く測定し得なかった刺激や部位に対して、新たな測定結果をもたらし得る。
【００３１】
またこのように多くのパラメータによって神経活動信号とＴ２＊信号の変動データとを結びつけているため、今後行われるであろう種々の実験や研究結果に対応できる発展性の高い脳機能測定装置１を提供できる。
【００３２】
さらにその演算過程において、ｏｘｙＨｂ濃度とｄｅｏｘｙＨｂ濃度との変動が算出されるため、そのことから逆に、ニューロン発火から血行動態変動までの過程を解明していく研究の一助にもなりうる。
【００３３】
加えて、本装置による脳機能測定において必要な予め定めておくべき各レスポンス関数は、近赤外分光法（ＮＩＲＳ）で実際に測定したｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度と、演算過程で算出されるものとが一致するように設定可能であるため、信頼度の非常に高いものとなる。具体的には、既にかなりの部分で解明されている脳の低次機能を利用して予め設定しておけばよい。
【００３４】
また、この近赤外分光法（ＮＩＲＳ）は測定が容易であるという利点を有するため、本装置１と組み合わせたシステムとし、被験者毎に前記各レスポンス関数を設定してより正確な脳機能測定を行うようにすることも無理なくできる。
【００３５】
なお、本発明は前記実施形態に限られるものではない。例えばｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度を、前記実施形態のように近赤外分光法（ＮＩＲＳ）に対応させるべく光の吸収と散乱に関する形式で算出する必要はなく、その他の形式で算出するようにしても構わない。
【００３６】
その他本発明は、上記図示例に限られず、その趣旨を逸脱しない範囲で種々の変更が可能である。
【００３７】
【発明の効果】
以上に詳述したように、本発明は、神経活動信号とＴ２＊信号の変動データとの相関を、ｏｘｙＨｂ濃度及びｄｅｏｘｙＨｂ濃度を介在させて求めるという、より生理学的な過程に近いものであると考えられるため、従来のＨＲＦという一定の波形を有した時間関数による脳機能測定では全く測定し得なかった刺激や部位に対して、新たな測定結果をもたらし得る。
【００３８】
また多くのパラメータによって神経活動信号とＴ２＊信号の変動データとを結びつけているため、今後行われるであろう種々の実験や研究結果に対応できる発展性の高い脳機能測定装置を提供できるうえ、その演算過程において、ｏｘｙＨｂ濃度とｄｅｏｘｙＨｂ濃度との変動が算出されるため、そのことから逆に、神経活動であるニューロン発火から血行動態変動までの過程を解明していく研究の一助にもなりうる。
【図面の簡単な説明】
【図１】請求項１に係る発明の機能構成を示す模式的機能構成図。
【図２】請求項４に係る発明の機能構成を示す模式的機能構成図。
【図３】本発明の一実施形態における脳機能測定装置の機能構成を示す模式的機能構成図。
【符号の説明】
１・・・脳機能測定装置
２・・・ｏｘｙＨｂ濃度算出部
３・・・ｄｅｏｘｙＨｂ算出部
４・・・比較データ算出部
【図１】
2005006774.tif 000004
【図２】
2005006774.tif 000005
【図３】
2005006774.tif 000006