(57)【特許請求の範囲】
【請求項１】
計測された脳の信号に関する情報である計測脳信号情報を格納し得る計測脳信号情報格納部と、
脳活動に関する事前の情報である脳事前情報を格納し得る脳事前情報格納部と、
アーチファクト源に関する事前の情報である１以上のアーチファクト事前情報を格納し得るアーチファクト事前情報格納部と、
前記計測脳信号情報と前記脳事前情報とを用いて、脳活動の情報である脳活動情報を推定し、取得する脳活動情報取得部と、
前記計測脳信号情報と、前記１以上のアーチファクト事前情報とを用いて、アーチファクト源の情報である１以上のアーチファクト情報を推定し、取得するアーチファクト情報取得部と、
前記計測脳信号情報と前記脳活動情報と前記１以上のアーチファクト情報とを用いて、前記計測脳信号情報からアーチファクトの影響を取り除く処理を行い、出力される脳活動の情報である出力脳活動情報を取得する出力脳活動情報取得部と、
前記出力脳活動情報を出力する出力部とを具備し、
前記脳活動情報取得部は、
複数のダイポール位置での電流についての情報である脳活動情報を、前記計測脳信号情報と前記脳事前情報とを用いた階層ベイズ推定を用いて推定し、取得し、
前記アーチファクト情報取得部は、
アーチファクト源のダイポール位置での電流についての情報である１以上のアーチファクト情報を、前記計測脳信号情報と、前記１以上のアーチファクト事前情報とを用いた階層ベイズ推定を用いて推定し、取得し、
前記出力脳活動情報取得部は、
前記計測脳信号情報から、前記脳活動情報のうちの脳の活動が無いと考えられる位置に配置されたダイポールに対応する脳活動情報と前記アーチファクト情報とを除く処理を行い、出力脳活動情報を取得する脳活動情報出力装置。
【請求項２】
前記脳事前情報は、
ｆＭＲＩ強度情報を電流分散とし、当該電流分布に対する事前知識を表す事前分布であり、
前記１以上の各アーチファクト事前情報は、
ｆＭＲＩ強度情報と比較して大きな電流分散を用い、かつ当該電流分布に対する事前知識を表す事前分布である請求項１記載の脳活動情報出力装置。
【請求項３】
前記計測脳信号情報は、
頭部周辺で観測された磁場に関する情報であり、
前記出力脳活動情報取得部は、
前記脳活動情報のうちの脳の活動が無いと考えられる位置に配置されたダイポールに対応する脳活動情報である電流の情報から磁場に関する情報である背景脳磁場情報を取得し、かつ、前記アーチファクト情報である電流の情報から磁場に関する情報であるアーチファクト磁場情報を取得し、
前記計測脳信号情報から、前記背景脳磁場情報と前記アーチファクト磁場情報とを除く処理を行い、出力脳活動情報を取得する請求項１から請求項２いずれか記載の脳活動情報出力装置。
【請求項４】
前記出力脳活動情報取得部が取得した出力脳活動情報と、前記計測脳信号情報とを用いて、当該出力脳活動情報におけるノイズ除去の精度に関する情報であるノイズ除去精度情報を取得するノイズ除去精度情報取得部をさらに具備し、
前記出力部は、
前記ノイズ除去精度情報が、予め格納されている閾値と比較して高い精度である場合のみ、前記出力脳活動情報を出力する請求項１から請求項３いずれか記載の脳活動情報出力装置。
【請求項５】
脳活動情報取得部、アーチファクト情報取得部、出力脳活動情報取得部、および出力部により実現される脳活動情報出力方法であって、
前記脳活動情報取得部により、記憶媒体に格納されている計測脳信号情報と脳事前情報とを用いて、脳活動の情報である脳活動情報を推定し、取得する脳活動情報取得ステップと、
前記アーチファクト情報取得部により、前記計測脳信号情報と、記憶媒体に格納されている１以上のアーチファクト事前情報とを用いて、アーチファクト源の情報である１以上のアーチファクト情報を推定し、取得するアーチファクト情報取得ステップと、
前記出力脳活動情報取得部により、前記計測脳信号情報と前記脳活動情報と前記１以上のアーチファクト情報とを用いて、前記計測脳信号情報からアーチファクトの影響を取り除く処理を行い、出力される脳活動の情報である出力脳活動情報を取得する出力脳活動情報取得ステップと、
前記出力部により、前記出力脳活動情報を出力する出力ステップとを具備し、
前記脳活動情報取得ステップにおいて、
複数のダイポール位置での電流についての情報である脳活動情報を、前記計測脳信号情報と前記脳事前情報とを用いた階層ベイズ推定を用いて推定し、取得し、
前記アーチファクト情報取得ステップにおいて、
アーチファクト源のダイポール位置での電流についての情報である１以上のアーチファクト情報を、前記計測脳信号情報と、前記１以上のアーチファクト事前情報とを用いた階層ベイズ推定を用いて推定し、取得し、
前記出力脳活動情報取得ステップにおいて、
前記計測脳信号情報から、前記脳活動情報のうちの脳の活動が無いと考えられる位置に配置されたダイポールに対応する脳活動情報と前記アーチファクト情報とを除く処理を行い、出力脳活動情報を取得する脳活動情報出力方法。
【請求項６】
コンピュータを、
記憶媒体に格納されている計測脳信号情報と脳事前情報とを用いて、脳活動の情報である脳活動情報を推定し、取得する脳活動情報取得部と、
前記計測脳信号情報と、記憶媒体に格納されている１以上のアーチファクト事前情報とを用いて、アーチファクト源の情報である１以上のアーチファクト情報を推定し、取得するアーチファクト情報取得部と、
前記計測脳信号情報と前記脳活動情報と前記１以上のアーチファクト情報とを用いて、前記計測脳信号情報からアーチファクトの影響を取り除く処理を行い、出力される脳活動の情報である出力脳活動情報を取得する出力脳活動情報取得部と、
前記出力脳活動情報を出力する出力部、として機能させるためのプログラムであって、
前記脳活動情報取得部は、
複数のダイポール位置での電流についての情報である脳活動情報を、前記計測脳信号情報と前記脳事前情報とを用いた階層ベイズ推定を用いて推定し、取得し、
前記アーチファクト情報取得部は、
アーチファクト源のダイポール位置での電流についての情報である１以上のアーチファクト情報を、前記計測脳信号情報と、前記１以上のアーチファクト事前情報とを用いた階層ベイズ推定を用いて推定し、取得し、
前記出力脳活動情報取得部は、
前記計測脳信号情報から、前記脳活動情報のうちの脳の活動が無いと考えられる位置に配置されたダイポールに対応する脳活動情報と前記アーチファクト情報とを除く処理を行い、出力脳活動情報を取得するものとして、コンピュータを機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、脳活動の信号である脳活動情報を出力する脳活動情報出力装置等に関するものである。
【背景技術】
【０００２】
従来、主成分分析や独立成分分析を用いた統計的な手法により、ＥＥＧまたはＭＥＧのアーチファクトを除去するアーチファクト除去の方法が提案されていた。しかし、かかる方法は、アーチファクトと同じような波形の脳活動が存在する場合、その脳活動の成分まで取り除いてしまう恐れがあり、有効な方法とは言えなかった。なお、アーチファクトとは、いわゆるノイズのことである。
【０００３】
また、アーチファクト源（眼球）にダイポールをおいて電流を推定し、それをＥＥＧから差し引くことでアーチファクト除去を行なう方法が提案されている（非特許文献１）。なお、アーチファクト源とは、いわゆるノイズ源のことである。
【非特許文献１】P. Berg and M. Scherg, Electroencephalogr Clin. Neurophysiol.,79, 1, pp. 36-44 (１99１).
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、非特許文献１における方法では、観測値にはアーチファクトだけでなく脳活動も含まれているので、必要な成分まで観測値から除いてしまう可能性があり、精度の高い脳活動情報が取得できなかった。
【課題を解決するための手段】
【０００５】
本第一の発明の脳活動情報出力装置は、計測された脳の信号に関する情報である計測脳信号情報を格納し得る計測脳信号情報格納部と、脳活動に関する事前の情報である脳事前情報を格納し得る脳事前情報格納部と、アーチファクト源に関する事前の情報である１以上のアーチファクト事前情報を格納し得るアーチファクト事前情報格納部と、前記計測脳信号情報と前記脳事前情報とを用いて、脳活動の情報である脳活動情報を推定し、取得する脳活動情報取得部と、前記計測脳信号情報と、前記１以上のアーチファクト事前情報とを用いて、アーチファクト源の情報である１以上のアーチファクト情報を推定し、取得するアーチファクト情報取得部と、前記計測脳信号情報と前記脳活動情報と前記１以上のアーチファクト情報とを用いて、前記計測脳信号情報からアーチファクトの影響を取り除く処理を行い、出力される脳活動の情報である出力脳活動情報を取得する出力脳活動情報取得部と、前記出力脳活動情報を出力する出力部とを具備する脳活動情報出力装置である。
【０００６】
かかる構成により、脳活動とアーチファクト源の事前情報を用いて、アーチファクトの影響を取り除き、精度の高い脳活動の情報を出力できる。
【０００７】
また、本第二の発明の脳活動情報出力装置は、第一の発明に対して、前記アーチファクト事前情報格納部は、複数のアーチファクト事前情報を格納しており、前記アーチファクト情報取得部は、複数のアーチファクト源ごとに、アーチファクト情報を推定し、複数のアーチファクト情報を取得し、前記出力脳活動情報取得部は、前記計測脳信号情報と前記脳活動情報と前記複数のアーチファクト情報とを用いて、前記計測脳信号情報から複数のアーチファクトの影響を取り除く処理を行い、出力される脳活動の情報である出力脳活動情報を取得する脳活動情報出力装置である。
【０００８】
かかる構成により、脳活動と複数のアーチファクト源の事前情報を用いて、アーチファクトの影響を取り除き、さらに精度の高い脳活動の情報を出力できる。
【０００９】
また、本第三の発明の脳活動情報出力装置は、第一、第二いずれかの発明に対して、前記脳事前情報、および前記１以上の各アーチファクト事前情報は、信号の強度に関する情報である強度情報を含み、前記脳事前情報の強度情報と、前記１以上の各アーチファクト事前情報の各強度情報は、異なることを特徴とする脳活動情報出力装置である。
【００１０】
かかる構成により、脳活動とアーチファクト源の事前情報であり、信号の強度に関する情報である強度情報を含む事前情報を用いて、アーチファクトの影響を取り除くことにより、さらに精度の高い脳活動の情報を出力できる。
【００１１】
また、本第四の発明の脳活動情報出力装置は、第一から第三いずれかの発明に対して、前記１以上のアーチファクト事前情報は、被験者の生体に含まれるアーチファクト源に関する事前の情報である生体アーチファクト事前情報を含む脳活動情報出力装置である。
【００１２】
かかる構成により、脳活動と生体に含まれるアーチファクト源の事前情報を用いて、アーチファクトの影響を取り除き、精度の高い脳活動の情報を出力できる。
【００１３】
また、本第五の発明の脳活動情報出力装置は、第四の発明に対して、前記１以上のアーチファクト事前情報は、被験者の生体に含まない外部のアーチファクト源に関する事前の情報である外部アーチファクト事前情報を含む脳活動情報出力装置である。
【００１４】
かかる構成により、脳活動と生体外のアーチファクト源の事前情報を用いて、アーチファクトの影響を取り除き、精度の高い脳活動の情報を出力できる。
【００１５】
また、本第六の発明の脳活動情報出力装置は、第一から第五いずれかの発明に対して、前記脳事前情報は、ｆＭＲＩ強度情報を電流分散とし、当該電流分布に対する事前知識を表す事前分布であり、前記１以上の各アーチファクト事前情報は、ｆＭＲＩ強度情報と比較して大きな電流分散を用い、かつ当該電流分布に対する事前知識を表す事前分布であり、前記脳活動情報取得部は、複数のダイポール位置での電流についての情報である脳活動情報を、前記計測脳信号情報と前記脳事前情報とを用いた階層ベイズ推定を用いて推定し、取得し、前記アーチファクト情報取得部は、アーチファクト源のダイポール位置での電流についての情報である１以上のアーチファクト情報を、前記計測脳信号情報と、前記１以上のアーチファクト事前情報とを用いた階層ベイズ推定を用いて推定し、取得し、前記出力脳活動情報取得部は、前記計測脳信号情報から、前記脳活動情報のうちの脳の活動が無いと考えられる位置に配置されたダイポールに対応する脳活動情報と前記アーチファクト情報とを除く処理を行い、出力脳活動情報を取得する脳活動情報出力装置である。
【００１６】
かかる構成により、脳活動の情報を電流レベルで推定できる。
【００１７】
また、本第七の発明の脳活動情報出力装置は、第一から第五いずれかの発明に対して、前記計測脳信号情報は、頭部周辺で観測された磁場に関する情報であり、前記出力脳活動情報取得部は、前記脳活動情報のうちの脳の活動が無いと考えられる位置に配置されたダイポールに対応する脳活動情報である電流の情報から磁場に関する情報である背景脳磁場情報を取得し、かつ、前記アーチファクト情報である電流の情報から磁場に関する情報であるアーチファクト磁場情報を取得し、前記計測脳信号情報から、前記背景脳磁場情報と前記アーチファクト磁場情報とを除く処理を行い、出力脳活動情報を取得する脳活動情報出力装置である。
【００１８】
かかる構成により、センサレベルで推定できる。
【００１９】
また、本第八の発明の脳活動情報出力装置は、第一から第七いずれかの発明に対して、前記出力脳活動情報取得部が取得した出力脳活動情報と、前記計測脳信号情報とを用いて、当該出力脳活動情報におけるノイズ除去の精度に関する情報であるノイズ除去精度情報を取得するノイズ除去精度情報取得部をさらに具備し、前記出力部は、前記ノイズ除去精度情報が、予め格納されている閾値と比較して高い精度である場合のみ、前記出力脳活動情報を出力する脳活動情報出力装置である。
【００２０】
かかる構成により、一定以上の精度であることを確認した上で、脳活動の情報を出力できる。
【発明の効果】
【００２１】
本発明による脳活動情報出力装置によれば、脳活動とアーチファクト源の事前情報を用いて、アーチファクトの影響を取り除き、精度の高い脳活動の情報を出力できる。
【発明を実施するための最良の形態】
【００２２】
以下、本発明に係る脳活動情報出力装置等の実施形態について図面を参照して説明する。なお、実施の形態において同じ符号を付した構成要素は同様の動作を行うので、再度の説明を省略する場合がある。
（実施の形態１）
【００２３】
本実施の形態において、脳活動とアーチファクト源の事前情報を用いて、アーチファクトの影響を取り除き、精度の高い脳活動の情報を出力できる脳活動情報出力システムについて説明する。また、本実施の形態において、電流レベルで脳活動の情報を推定する脳活動情報出力システムについて説明する。また、本実施の形態において、アーチファクト源が複数存在する場合について説明する。また、本実施の形態において、脳事前情報とアーチファクト事前情報との信号の大きさが異なることを想定している脳活動情報出力システムについて説明する。また、アーチファクト源は、生体である眼球、心臓、筋電などであったり、ＴＶの電波などの外部のものであったりする。さらに、本実施の形態において、脳活動信号をチェックし、精度が一定以上であることを確認後、脳信号を出力について説明する。なお、アーチファクト事前情報とはノイズ事前情報と同意義である。
【００２４】
図１は、本実施の形態における脳活動情報出力システムのブロック図である。脳活動情報出力システムは、脳信号計測装置１１、脳活動情報出力装置１２を具備する。
【００２５】
脳信号計測装置１１は、脳信号を計測する装置である。脳信号計測装置１１は、例えば、脳磁計測装置である。脳磁計測装置は、脳活動に関する事象を観測する観測手段としてのブレインキャップ１１０１を含んでも良い。ブレインキャップ１１０１は、例えば、２０８チャンネルセンサを有する。ただし、ブレインキャップ１１０１のチャンネル数は問わない。なお、脳信号とは、磁場でも良いし、電流等でも良い。
【００２６】
脳活動情報出力装置１２は、計測脳信号情報格納部１２１、脳事前情報格納部１２２、アーチファクト事前情報格納部１２３、脳活動情報取得部１２４、アーチファクト情報取得部１２５、出力脳活動情報取得部１２６、ノイズ除去精度情報取得部１２７、出力部１２８を具備する。
【００２７】
脳活動情報出力装置１２は、脳信号計測装置１１が取得した脳の信号に関する情報である計測脳信号情報から、アーチファクト除去を行い、精度の高い脳活動情報を出力する装置である。
【００２８】
計測脳信号情報格納部１２１は、計測された脳の信号に関する情報である計測脳信号情報を格納し得る。計測脳信号情報は、脳磁計測装置（脳信号計測装置１１の一例）の測定結果でも良いし、当該計測結果から、低周波のドリフト除去を行った情報でも良い。また、計測脳信号情報は、ＭＥＧデータでも良いし、電流の情報でも良い。また、通常、計測脳信号情報は、複数地点（例えば、２０８チャンネル）の情報であり、計測脳信号情報格納部１２１には、複数の脳信号情報が格納されている。
【００２９】
計測脳信号情報は、例えば、被験者の脳神経活動の電気的変化に伴って、被験者の頭部周辺で観測された微弱な磁場である。例えば、神経電流を位置と方向が固定されたＮ個の電流ダイポールでモデル化し、その電流モーメントの強度をＮ次元ベクトルＪで表すＭ個のＭＥＧセンサ上での観測磁場をＭ次元ベクトルＢで表すと、次式が成り立つ。
【数１】
0005266597.tif 000002
【００３０】
このＭ次元ベクトルＢが、計測脳信号情報の例である。また、上記のＮ次元ベクトルＪを計測脳信号情報の例としても良い。ここで、Ｇはリードフィールド行列と呼ばれ、単位電流モーメントによって発生する磁場を表す行列である。脳を囲む頭蓋内部が均一な媒質の球であると仮定すると、リードフィールドの成分はＳａｒｖａｓの式によって解析的に計算できる。ＭＥＧ信号を生成する主要な脳活動成分は錐体細胞に流れる電流であるという知見に基づき、皮質表面上に垂直な向きを持つ電流ダイポールを仮定する、とする。なお、Ｓａｒｖａｓの式は、公知の式であるので、詳細な説明は省略する。
【００３１】
計測脳信号情報格納部１２１は、不揮発性の記録媒体が好適であるが、揮発性の記録媒体でも実現可能である。計測脳信号情報格納部１２１に計測脳信号情報が記憶される過程は問わない。例えば、記録媒体を介して計測脳信号情報が計測脳信号情報格納部１２１で記憶されるようになってもよく、脳信号計測装置１１から通信回線等を介して送信された計測脳信号情報が計測脳信号情報格納部１２１で記憶されるようになってもよく、あるいは、入力デバイスを介して入力された計測脳信号情報が計測脳信号情報格納部１２１で記憶されるようになってもよい。
【００３２】
脳事前情報格納部１２２は、脳活動に関する事前の情報である脳事前情報を格納し得る。脳事前情報は、脳の電流に対する事前分布であり、事前知識に基づいて定められる情報である。脳事前情報は、例えば、後述するＰ０（Ｊ）である。脳事前情報は、例えば、階層ベイズ推定を用いて、図示しない導出手段により導出されたものである。階層ベイズ推定では、次のようなダイポールに対する階層事前分布を取得する。また、脳事前情報は、信号の強度に関する情報である強度情報を含むことは好適である。さらに、脳事前情報は、ｆＭＲＩ強度情報を電流分散とし、当該電流分布に対する事前知識を表す事前分布であっても良い。
【００３３】
以下、階層事前分布の取得方法を説明する前に、階層ベイズ推定を用いて、ダイポール位置での電流を推定する方法について説明する。
【００３４】
つまり、脳表面近くのＭ点での測定した磁場をＢ＝（Ｂ１，...，ＢＭ）Ｔ、脳表面上に仮定したダイポールに流れる電流分布をＪ＝｛Ｊｉ｜ｉ＝１，...，Ｎ｝と表す、とする。また、階層ベイズ推定において、電流分布Ｊが与えられたときに磁場Ｂが観測される確率Ｐ（Ｂ｜Ｊ）が、数式２の正規分布に従うと仮定する。
【数２】
0005266597.tif 000003
【００３５】
数式２において、βはノイズ分散の逆数である。また、磁場Ｂを観測したときの電流Ｊの事後分布（Ｐ（Ｊ｜Ｂ））は次の数式３で表される。
【数３】
0005266597.tif 000004
【００３６】
数式３において、Ｐ（Ｂ）は周辺尤度である。Ｐ０（Ｊ）は電流に対する事前分布であり、所定の事前知識に基づいて定められる。階層ベイズ推定では、以下の数式４から数式６のようなダイポールに対する階層事前分布が導入される。
【数４】
0005266597.tif 000005
【数５】
0005266597.tif 000006
【数６】
0005266597.tif 000007
【００３７】
数式４において、Ａはその対角成分がα＝｛αｉ｜ｉ＝１，...，Ｎ｝である対角行列であり、α?１ｉがＪｉの事前電流分散に対応する。ここで、事前電流分散を表すパラメータを、以下の数式７に示すものとする。
【数７】
0005266597.tif 000008
【００３８】
階層事前分布のパラメータであるν０ｉバー（バーは、νの上の−を有する文字を意味する）とγ０ｉはハイパーパラメータと呼ばれている。これらの値はｆＭＲＩの活動強度情報（ｔ値）、倍率のパラメータｍ０、信頼度のパラメータγ０・γ０ｉから定義されている。ν０ｉバーはｆＭＲＩのｔ値に応じて設定され、背景電流分散をν０，ｂａｓｅとすると、ν０ｉバーは、数式８のように表される。なお、ハイパーパラメータは、信号の強度に関する情報である強度情報の例である。
【数８】
0005266597.tif 000009
【００３９】
ここで、ｔｉハット（ハットは、ｔの上に＾が存在することを示す）は最大値が１となるようにｆＭＲＩのｔ値を正規化した値である。
【００４０】
νｉバーが大きいと、電流の分散が大きくなる確率が高くなる（ｉ番目のダイポールに電流が流れやすい、という事前情報を与えることを意味する）。ハイパーパラメータγ０ｉは階層事前分布の信頼度を制御する。γ０ｉが小さいほど信頼度が低く、大きいほど信頼度が高いことを表す。
【００４１】
脳事前情報は、例えば、図示しない脳事前情報取得手段により、ｆＭＲＩ活動強度情報を数式８に代入し，ｍ０＝１００、γ０ｉ＝１００として事前分散の推定を行い、取得された情報である。なお、脳事前情報は、皮質上に配置したダイポールに対して、事前分散の推定を行い、取得された情報である。
【００４２】
脳事前情報格納部１２２は、不揮発性の記録媒体が好適であるが、揮発性の記録媒体でも実現可能である。脳事前情報格納部１２２に脳事前情報が記憶される過程は問わない。例えば、記録媒体を介して脳事前情報が脳事前情報格納部１２２で記憶されるようになってもよく、通信回線等を介して送信された脳事前情報が脳事前情報格納部１２２で記憶されるようになってもよく、あるいは、入力デバイスを介して入力された脳事前情報が脳事前情報格納部１２２で記憶されるようになってもよい。
【００４３】
アーチファクト事前情報格納部１２３は、アーチファクト源に関する事前の情報である１以上のアーチファクト事前情報を格納し得る。アーチファクト事前情報格納部１２３は、複数のアーチファクト事前情報を格納していることは好適である。アーチファクト源とは、例えば、生体の一部である眼球、咀嚼筋、心拍などである。また、アーチファクト源とは、例えば、生体以外の外部のアーチファクト源であるテレビ、輸送物（自動車、トラック、電車等）の通過などでも良い。つまり、アーチファクト事前情報は、被験者の生体に含まれるアーチファクト源に関する事前の情報である生体アーチファクト事前情報を含む。また、アーチファクト事前情報は、被験者の生体に含まない外部のアーチファクト源に関する事前の情報である外部アーチファクト事前情報を含んでも良い。アーチファクト事前情報は、脳事前情報と同様に、例えば、階層ベイズ推定を用いて得られる情報であり、数式４から数式６のようなダイポールに対する階層事前分布である。ただし、アーチファクト事前情報は、脳事前情報とは、電流分散の大きさ（信号の強度である強度情報）が異なる。例えば、各種のアーチファクト事前情報は、数式９、１０のように表される。また、アーチファクト事前情報は、強度情報を含むことは好適である。また、脳事前情報の強度情報と、１以上の各アーチファクト事前情報の各強度情報は、異なることは好適である。さらに、各アーチファクト事前情報は、ｆＭＲＩ強度情報と比較して大きな電流分散を用い、かつ当該電流分布に対する事前知識を表す事前分布であることは好適である。
【数９】
0005266597.tif 000010
【数１０】
0005266597.tif 000011
【００４４】
ここで，νｅｙｅバー（バーとは、νの上の−を有することを示す）、νｃａｒｄｉａｃバー、νｍａｓｔｉｃａｔｏｒｙバーはそれぞれ眼球中心、心臓位置、咀嚼筋に配置したダイポールに関するハイパーパラメータν０ｉバーである。また、γｅｙｅ、γｃａｒｄｉａｃ、γｍａｓｔｉｃａｔｏｒｙは、ハイパーパラメータγ０ｉに対して同様の値を表す。
【００４５】
なお、νｂｒａｉｎバー：νｅｙｅバー：νｃａｒｄｉａｃバー：νｍａｓｔｉｃａｔｏｒｙバーは、例えば、「１：５０：５０００：１００」などの簡易な比でも良い。
【００４６】
また、これらの事前電流分散は、例えば、図示しないアーチファクト事前情報取得手段が、観測磁場と復元磁場のあてはまりの度合（Ｇｏｏｄｎｅｓｓ−ｏｆ−ｆｉｔ，ＧｏＦ）である第一の指標、ｆＭＲＩの活動強度情報と平均電流分散の相関係数である第二の指標、の２種類の指標をもとに決定したものであることは好適である。
【００４７】
具体的には、ＧｏＦは数式１１により定義される量とした。
【数１１】
0005266597.tif 000012
【００４８】
数式１１において、チャンネルｊでの計測磁場をＢｏｂｓ（ｊ）、チャンネルｊにおける復元磁場をＢｅｓｔ（ｊ）とする。Ｍはセンサ数とする。推定電流から計算される磁場（復元磁場）と観測磁場が一致しているならばＧｏＦは大きくなる（１００に近づく）。推定結果と計測結果が一致するようなハイパーパラメータを選択するための判断材料としてＧｏＦを用いた。第二の指標については、ｆＭＲＩの活動が期待される位置では、脳内での神経活動も大きくなることが期待される。そのため、ｆＭＲＩの活動強度情報と平均電流分散には対応関係があると考えられる（ふたつの相関係数が最も大きくなることが期待される）。第一の指標と第二の指標の条件をできるだけ満たすようなハイパーパラメータの組を選択することでアーチファクト源の電流強度に関する粗い情報をあらかじめ設定し、アーチファクト事前情報取得手段は、ダイポールに関する事前分散（アーチファクト事前情報）を推定し、取得するものとする。アーチファクト事前情報格納部１２３は、不揮発性の記録媒体が好適であるが、揮発性の記録媒体でも実現可能である。アーチファクト事前情報格納部１２３にアーチファクト事前情報が記憶される過程は問わない。例えば、記録媒体を介してアーチファクト事前情報がアーチファクト事前情報格納部１２３で記憶されるようになってもよく、通信回線等を介して送信されたアーチファクト事前情報がアーチファクト事前情報格納部１２３で記憶されるようになってもよく、あるいは、入力デバイスを介して入力されたアーチファクト事前情報がアーチファクト事前情報格納部１２３で記憶されるようになってもよい。
【００４９】
脳活動情報取得部１２４は、計測脳信号情報格納部１２１に格納されている計測脳信号情報と、脳事前情報格納部１２２に格納されている脳事前情報とを用いて、脳活動の情報である脳活動情報を推定し、取得する。脳活動情報取得部１２４は、複数のダイポール位置での電流についての情報である脳活動情報を、計測脳信号情報（例えば、Ｂ）と脳事前情報（例えば、Ｐ０（Ｊ））とを用いた階層ベイズ推定を用いて推定し、取得することは好適である。例えば、脳活動情報取得部１２４は、上述した数式３の情報を格納しており、当該数式３を用いて、電流Ｊの事後分布（Ｐ（Ｊ｜Ｂ））を取得する。電流Ｊの事後分布（Ｐ（Ｊ｜Ｂ））は、脳活動情報である。脳活動情報取得部１２４は、通常、ＭＰＵやメモリ等から実現され得る。脳活動情報取得部１２４の処理手順は、通常、ソフトウェアで実現され、当該ソフトウェアはＲＯＭ等の記録媒体に記録されている。但し、ハードウェア（専用回路）で実現しても良い。
【００５０】
アーチファクト情報取得部１２５は、計測脳信号情報と、１以上のアーチファクト事前情報とを用いて、アーチファクト源の情報である１以上のアーチファクト情報を推定し、取得する。アーチファクト情報取得部１２５は、複数のアーチファクト源ごとに、アーチファクト情報を推定し、複数のアーチファクト情報を取得することは好適である。アーチファクト情報取得部１２５は、アーチファクト源のダイポール位置での電流についての情報である１以上のアーチファクト情報を、計測脳信号情報と、１以上のアーチファクト事前情報とを用いた階層ベイズ推定を用いて推定し、取得することは好適である。アーチファクト情報取得部１２５は、事前情報が異なるが、脳活動情報取得部１２４と同様の動作を行う。つまり、アーチファクト情報取得部１２５は、数式３を用いて、各アーチファクト源に対応する電流Ｊの事後分布（Ｐ（Ｊ｜Ｂ））を取得する。アーチファクト情報取得部１２５は、通常、ＭＰＵやメモリ等から実現され得る。アーチファクト情報取得部１２５の処理手順は、通常、ソフトウェアで実現され、当該ソフトウェアはＲＯＭ等の記録媒体に記録されている。但し、ハードウェア（専用回路）で実現しても良い。
【００５１】
出力脳活動情報取得部１２６は、計測脳信号情報と脳活動情報と１以上のアーチファクト情報とを用いて、計測脳信号情報からアーチファクトの影響を取り除く処理を行い、出力される脳活動の情報である出力脳活動情報を取得する。つまり、アーチファクトの影響を取り除くために、脳皮質のダイポールとアーチファクト源のダイポールの推定を同時に行うこととなる。なお、脳皮質のダイポールの推定とは、脳活動情報取得部１２４が脳活動情報を取得する処理である。アーチファクト源のダイポールの推定とは、アーチファクト情報取得部１２５が、アーチファクト源の情報である１以上のアーチファクト情報を取得する処理である。
【００５２】
また、出力脳活動情報取得部１２６は、計測脳信号情報と脳活動情報と複数のアーチファクト情報とを用いて、計測脳信号情報から複数のアーチファクトの影響を取り除く処理を行い、出力される脳活動の情報である出力脳活動情報を取得することは好適である。
【００５３】
さらに、出力脳活動情報取得部１２６は、計測脳信号情報から、脳活動情報のうちの脳の活動が無いと考えられる位置に配置されたダイポールに対応する脳活動情報とアーチファクト情報とを除く処理を行い、出力脳活動情報を取得することは好適である。なお、脳の活動が無いと考えられる位置に配置されたダイポールの情報は、予め設定されていても良いし、ユーザが手入力により与えるなどしても良い。
【００５４】
以下、出力脳活動情報取得部１２６等の具体的な動作例について説明する。今、脳皮質とアーチファクト源のダイポールをそれぞれＪｂｒａｉｎ，Ｊｎｏｉｓｅ，リードフィールド行列をＧｂｒａｉｎ，Ｇｎｏｉｓｅ，ダイポールの事前分散をＡ?１ｂｒａｉｎ，Ａ?１ｎｏｉｓｅとすると、以下の式１２から１４を用いて、脳活動情報取得部１２４、およびアーチファクト情報取得部１２５は、それぞれ脳皮質とアーチファクト源のダイポールＪを推定する。
【数１２】
0005266597.tif 000013
【数１３】
0005266597.tif 000014
【数１４】
0005266597.tif 000015
【００５５】
つまり、脳活動情報取得部１２４とアーチファクト情報取得部１２５により、階層ベイズ推定法でアーチファクト源と脳皮質のダイポールを同時に推定した後に、出力脳活動情報取得部１２６は、ｆＭＲＩの活動がない位置に配置されたダイポールのみが作る磁場を、数式１５にしたがって観測磁場Ｂから取り除き、Ｂｐｒｕｎｅｄを算出し、メモリ上に配置する。なお、観測磁場Ｂは、計測脳信号情報の一例であり、Ｂｐｒｕｎｅｄは、出力脳活動情報の一例である。
【数１５】
0005266597.tif 000016
【００５６】
出力脳活動情報取得部１２６は、通常、ＭＰＵやメモリ等から実現され得る。出力脳活動情報取得部１２６の処理手順は、通常、ソフトウェアで実現され、当該ソフトウェアはＲＯＭ等の記録媒体に記録されている。但し、ハードウェア（専用回路）で実現しても良い。
【００５７】
ノイズ除去精度情報取得部１２７は、出力脳活動情報取得部１２６が取得した出力脳活動情報と、計測脳信号情報とを用いて、当該出力脳活動情報におけるノイズ除去の精度に関する情報であるノイズ除去精度情報を取得する。ノイズ除去精度情報取得部１２７の具体的な処理の例は、後述する。ノイズ除去精度情報取得部１２７は、通常、ＭＰＵやメモリ等から実現され得る。ノイズ除去精度情報取得部１２７の処理手順は、通常、ソフトウェアで実現され、当該ソフトウェアはＲＯＭ等の記録媒体に記録されている。但し、ハードウェア（専用回路）で実現しても良い。
【００５８】
出力部１２８は、出力脳活動情報を出力する。出力部１２８は、ノイズ除去精度情報が、予め格納されている閾値と比較して高い精度である場合のみ、出力脳活動情報を出力しても良い。ここで、出力とは、ディスプレイへの表示、プロジェクターを用いた投影、プリンタへの印字、音出力、外部の装置への送信、記録媒体への蓄積、他の処理装置や他のプログラム等への処理結果の引渡し等を含む概念である。出力部１２８は、ディスプレイやスピーカー等の出力デバイスを含むと考えても含まないと考えても良い。出力部１２８は、出力デバイスのドライバーソフトまたは、出力デバイスのドライバーソフトと出力デバイス等で実現され得る。
【００５９】
次に、脳活動情報出力システムを構成する脳活動情報出力装置１２の動作について、図２のフローチャートを用いて説明する。
【００６０】
（ステップＳ２０１）脳活動情報取得部１２４は、計測脳信号情報格納部１２１から計測脳信号情報を読み出す。
【００６１】
（ステップＳ２０２）脳活動情報取得部１２４は、脳事前情報格納部１２２から脳事前情報を読み出す。
【００６２】
（ステップＳ２０３）脳活動情報取得部１２４は、予め格納している数式であり、脳活動情報等の推定のための数式を、メモリ上に読み出す。かかる数式は、例えば、階層ベイズ推定を実現するための式である。
【００６３】
（ステップＳ２０４）脳活動情報取得部１２４は、ステップＳ２０３で読み出した式に対して、ステップＳ２０１で読み出した計測脳信号情報と、ステップＳ２０２で読み出した脳事前情報を代入し、当該式を実行し、脳活動情報を取得し、メモリ上に配置する。
【００６４】
（ステップＳ２０５）アーチファクト情報取得部１２５は、カウンタｉに１を代入する。
【００６５】
（ステップＳ２０６）アーチファクト情報取得部１２５は、アーチファクト事前情報格納部１２３に、ｉ番目のアーチファクト事前情報が存在するか否かを判断する。ｉ番目のアーチファクト事前情報が存在すればステップＳ２０７に行き、存在しなければステップＳ２１０に行く。
【００６６】
（ステップＳ２０７）アーチファクト情報取得部１２５は、アーチファクト事前情報格納部１２３から、ｉ番目のアーチファクト事前情報を読み出す。
【００６７】
（ステップＳ２０８）アーチファクト情報取得部１２５は、ステップＳ２０３でメモリ上に読み出された数式に、ステップＳ２０１で読み出された計測脳信号情報と、ステップＳ２０７で読み出したアーチファクト事前情報を代入し、当該式を実行し、アーチファクト情報を取得し、メモリ上に配置する。
【００６８】
（ステップＳ２０９）アーチファクト情報取得部１２５は、カウンタｉを１、インクリメントし、ステップＳ２０６に戻る。
【００６９】
（ステップＳ２１０）出力脳活動情報取得部１２６は、ステップＳ２０４で取得された脳活動情報から、ｆＭＲＩ活動のない位置のダイポールに対する情報を取り除く。
【００７０】
（ステップＳ２１１）出力脳活動情報取得部１２６は、ステップＳ２１０で残った情報から、ステップＳ２０８で取得された１以上のアーチファクト情報を取り除く。そして、出力脳活動情報取得部１２６は、出力脳活動情報を得る。
【００７１】
（ステップＳ２１２）出力部１２８は、ステップＳ２１１で取得された出力脳活動情報を出力し、処理を終了する。
【００７２】
なお、図２のフローチャートにおいて、ステップＳ２１０、ステップＳ２１１における処理の順序、タイミングは問わない。ステップＳ２１０、ステップＳ２１１における処理は、例えば、上述した数式１５を用いて行われる。
【００７３】
また、図２のフローチャートにおいて、ノイズ除去精度情報取得部１２７の処理は行わなかった。しかし、図２のフローチャートにおいて、出力部１２８は、ノイズ除去精度情報が、予め格納されている閾値と比較して高い精度である場合のみ、出力脳活動情報を出力しても良い。
【００７４】
以下、本実施の形態における脳活動情報出力システムを用いた実験結果について説明する。被験者は、正常な視覚を有する健常男性３名である。被験者３名（Ｘ，Ｙ，Ｚ）が、ＭＥＧ計測実験とｆＭＲＩ計測実験の両方に参加した。実験に先だって、被験者には実験内容についての説明を行い、同意を得た上で実験を行った。
【００７５】
ＭＥＧ計測実験とは、以下のような実験である。脳信号計測装置１１は、全頭型の２０８チャンネルセンサを有する市販の脳磁計測装置である。そして、脳信号計測装置１１により、１０００ＨｚでＭＥＧデータが記録した。なお、ＭＥＧデータは、計測脳信号情報の一例であるとも考えられる。
【００７６】
そして、実験課題として、被験者は運動する視標を、固視点を注視したまま、心的に追従する課題を行った（ｃｏｖｅｒｔ  ｐｕｒｓｕｉｔ  ｔａｓｋ）。具体的には、以下の（ａ）から（ｄ）の流れで課題を遂行した。
【００７７】
（ａ）被験者は、準備が出来たら、タスク開始のボタンを押す。（ｂ）固視点と視標が画面中央に現れる（１０００±５００ｍｓ間待つ）。ターゲットは画面中央で静止している。被験者は固視点を見つめ続けるように指示された。見つめる時間は４０００ｍｓとした（ｐｒｅｔｒｉｇｇｅｒ期間）。（ｃ）その後、視標が水平方向に動き始める。被験者は、固視点を見つめながら視標を心の中で追従するように指示された。視標は４０００ｍｓ間動き続けた。（ｄ）視標が消え、被験者は休憩する。そして、（ａ）から（ｄ）の流れを１試行とした。指標の運動パターンは３種類あり、それぞれ３０試行ずつランダムに提示した（３種類×３０試行＝合計９０試行）。これを１セッションとして、被験者は合計で３セッション（１セッション９０試行×３セッション＝合計２７０試行）課題を遂行した。
【００７８】
次に、ｆＭＲＩ計測実験について説明する。ｆＭＲＩ計測には市販の１．５ＴのＭＲＩスキャナを用いた。視覚刺激は眼前から３７０ｍｍ離れた位置に設置したスクリーンに提示した。被験者は、正弦波の動きをする指標を、目を動かさずに心的に追従するテストブロックと、静止する固視点を注視しつづけるコントロールブロックを、それぞれ１８秒間ずつ繰り返した。テストブロックの指標は、振幅４ｄｅｇ、周波数０．５Ｈｚの正弦波運動する赤色の点とした。
【００７９】
まず、最初に、脳信号計測装置１１が計測し、取得したＭＥＧデータは、リファレンスセンサの値を用いて、図示しない低周波ドリフト除去手段（脳信号計測装置１１が有する）により、低周波のドリフトの影響が除去された。次に、図示しないデータ切出手段（脳信号計測装置１１が有する）により、指標の提示時刻を基準にして−４００〜＋４０００ｍｓのデータを切り出し、これを１試行のデータとした。データ切出手段により、切り出されたＭＥＧデータの最初の４００ｍｓ間は指標が動き出す前の状態であり（ｐｒｅｔｒｉｇｇｅｒ期間）、図示しない補正手段（脳信号計測装置１１が有する）により、この期間の平均値を用いて基線補正を行った。上記の前処理をしたＭＥＧデータが、計測脳信号情報格納部１２１に格納された計測脳信号情報の例であり、後述する電流推定に用いられる。
【００８０】
また、解析の際、タスク遂行中に瞬目の影響のあるデータは除外した。また、実験中に被験者が熟睡してしまったセッションのデータも除外した。そのため、被験者Ｘ、Ｙ、Ｚから計測されたデータのうち、それぞれ５３．３％、９３．０％、７６．３％のデータが解析に用いられた。
【００８１】
ｆＭＲＩデータの解析はＳＰＭ５（Ｔｈｅ  Ｗｅｌｌｃｏｍｅ  Ｄｅｐａｒｔｍｅｎｔ  ｏｆ  Ｃｏｇｎｉｔｉｖｅ  Ｎｅｕｒｏｌｏｇｙ）を用いて行われた。静止する固視点のみが提示され、それを注視している条件と、固視点を注視したままその上側で運動する指標を心の中で追従している条件での脳活動の有意差検定をｔ検定で行った。
【００８２】
あらかじめ計測しておいた被験者のＭＲＩ構造画像をもとに、Ｂｒａｉｎ  Ｖｏｙａｇｅｒ（Ｂｒａｉｎ  Ｉｎｎｏｖａｔｉｏｎ）を用いて大脳表面モデルを作成した。その全脳モデル上に、被験者Ｘは１８８６個、被験者Ｙは３２２９個、被験者Ｚは２７８９個のダイポールが配置された。このとき、タスクに関連した脳活動が期待される領域（ｆＭＲＩでのｔ検定の有意差がｐ＜０．０１［ｕｎｃｏｒｒｅｃｔｅｄ］）は密になるように、ダイポールが配置された。その結果，ｆＭＲＩの活動が期待される領域では２．１５±１．３１（ｍｅａｎ±ＳＤ）［ｍｍ］，それ以外の領域では、３．９１±３．１９［ｍｍ］の間隔でダイポールが配置された。
【００８３】
また、本実験では、主なアーチファクト源として、眼球付近、咀嚼筋、心拍を仮定し、それぞれの場所にダイポールを仮定した。眼球の中心位置、咀嚼筋の位置は、ＭＲＩの構造画像から割り出し、左右に一つずつダイポールを配置した。心臓の位置は、前交連（ａｎｔｅｒｉｏｒ  ｃｏｍｍｉｓｓｕｒｅ，ＡＣ）を原点として、前交連と後交連（ｐｏｓｔｅｒｉｏｒ  ｃｏｍｍｉｓｓｕｒｅ，ＰＣ）を結ぶ線をｙ軸とする座標系を考え（ｘ軸は正中矢状面に垂直な左右方向，ｚ軸はこれらに垂直な上下方向），（ｘ，ｙ，ｚ）＝（−５，０，−２３）［ｃｍ］の位置に配置した。
【００８４】
計測脳信号情報格納部１２１に格納されているＭＥＧデータから、上述した階層ベイズ推定によって電流推定を行う際、ｆＭＲＩ活動強度情報（ｔ値）を階層ベイズ推定の事前知識として用いた。つまり、脳事前情報格納部１２２には、ｆＭＲＩ活動強度情報（ｔ値）を含む脳事前情報が格納されている。また、アーチファクト事前情報格納部１２３には、ｆＭＲＩ活動強度情報（ｔ値）を含む、３つ（眼球付近、咀嚼筋、心拍に対応）のアーチファクト事前情報が格納されている。このとき、ｃｏｖｅｒｔ  ｐｕｒｓｕｉｔ  ｔａｓｋ  に関連のある領域に絞って事前分布を与えることはせず、ｆＭＲＩの解析によって得られた活動強度情報に基づいて事前分布が設定された。
【００８５】
以上により、計測脳信号情報格納部１２１、脳事前情報格納部１２２、およびアーチファクト事前情報格納部１２３に、それぞれの情報が格納された。
【００８６】
次に、出力脳活動情報取得部１２６等は、数式１２から数式１５を用いて、上述したようなアーチファクト除去を行った。そして、例えば、図示いない平滑化手段により、アーチファクトを除去されたデータに対して、１００ｍｓの時間幅で移動平均を行い、データを平滑化した。
【００８７】
なお、本実験において、計測したＭＥＧデータにＰＣＡをかけて、その第一主成分がＭＣＧ（心磁図）波形に類似していることを目視で確認した上で、その成分を除去した。なお、ＰＣＡとは、主成分分析のことである。
（実験結果１）
【００８８】
脳活動情報出力装置１２は、脳活動とアーチファクトの振舞いを同時に予測し、ノイズ除去を行うが、その精度を調べるために、被験者を伴うＭＥＧ計測実験を行った。被験者は、目の前に提示された動く視標を心的に追従する課題を実行するように指示され，そのときのＭＥＧが、脳信号計測装置１１により計測された。
【００８９】
このようなタスクを遂行する際に、ＭＥＧデータは、心拍や眼球アーチファクト、筋活動の影響を受けてしまう。そこで、それらのアーチファクト源と大脳皮質表面に対してダイポールを配置し、脳活動情報出力装置１２により、それらの電流を推定することを試みた。その際、アーチファクト源や脳活動の信号強度は、上述したように、それぞれ異なることを考慮し、信号強度に関する情報を事前知識として、脳事前情報格納部１２２またはアーチファクト事前情報格納部１２３に、個別に蓄積し、出力脳活動情報取得部１２６が電流の推定を行った。そして、アーチファクト情報取得部１２５は、推定した電流からアーチファクトの磁場を復元し、出力脳活動情報取得部１２６が、計測データからアーチファクトの磁場を差し引くことによってアーチファクト除去を行い、出力脳活動情報を取得した。
【００９０】
図３のＡは、頭部のどの位置で、ＭＥＧデータがどのような時間波形になるかを表している。アーチファクトを除去する前は頭部外側付近や眼球付近で波形が乱れている。しかし、アーチファクト除去後は、その乱れが小さくなっている（図３のＢ）。また、頭部中心付近は、外側部分に比べればアーチファクトの影響が小さいと予想されるが、そのような箇所では，ノイズ除去の前後で波形の特徴は変わらなかった（図３のＣ）。
【００９１】
次に、ノイズ除去精度情報取得部１２７は、ノイズ除去したデータと、（脳活動のみを反映した）真の値とのずれを、例えば、以下のように検知する。つまり、計測したデータを試行全体にわたって加算平均することによって、タスクに関係のない成分の影響を抑えることができる。ノイズ除去精度情報取得部１２７は、この値を、真の値と仮定し、その値からどの程度ずれが生じているかを検出する。具体的には、以下の数式１６を用いて、ノイズ除去精度情報取得部１２７は、計測データ、ＰＣＡ、ＶＢＭＥＧでクリーニングしたデータの各センサでのＲＭＳＥ（ｒｏｏｔ  ｍｅａｎｓｑｕａｒｅ  ｅｒｒｏｒ）を計算する。なお、ＶＢＭＥＧとは、階層ベイズ推定を行なうためのソフトウェアのことである。
【数１６】
0005266597.tif 000017
【００９２】
数式１６において、Ｂは単一試行のＭＥＧデータ、Ｂハット（ハットはＢの上に＾）はＢの試行平均を表す。数式１６において、Ｔはサンプル数、Ｋは試行回数を表す。ノイズ除去精度情報取得部１２７が行った、数式１６を用いた計算結果を図４に示す。図４のＡ〜Ｃは典型的な被験者の視標パターン１のＲＭＳＥの試行平均値を頭部マップ上に表示した図である。
【００９３】
計測データ（図４のＡ）、ＰＣＡ（図４のＢ）では頭部外側からくるアーチファクトの影響が見られるが、脳活動情報出力装置１２における拡張ダイポールを用いた方法では、その影響が小さいことが分かる（図４のＣ）。
【００９４】
図４のＤは、頭部外周部分（７１個）のセンサのＲＭＳＥが観測値とＰＣＡ、ＶＢＭＥＧを用いた方法でＲＭＳＥの大きさを比較した図である。３つの解析方法について、分散分析を行うと、有意差がみられた（Ｆ（２，６３８）＝７６２．６１）。多重比較（ＴｕｋｅｙのＨＳＤ検定）を行うと、すべてのデータ群の間に有意差があり（［Ｏｂｓｅｒｖｅｄ］−［ＰＣＡ］：ｐ＜０．０１，［Ｏｂｓｅｒｖｅｄ］−［ＶＢＭＥＧ］：ｐ＜０．０１，［ＰＣＡ］−［ＶＢＭＥＧ］：ｐ＜０．０１），ＶＢＭＥＧを使ってアーチファクト除去をした場合にＲＭＳＥが最も小さかった（図４のＤ）。
【００９５】
以上より、脳活動情報出力装置１２は、頭部外側からのアーチファクトの影響をより良く除去できたことが分かる。
（実験結果２）
【００９６】
実験結果２では、脳内の格子点で推定された電流レベルでのアーチファクト除去が行われ、出力される脳活動情報が、電流レベルで推定精度が向上していることを示す。
【００９７】
図５は、視標を心的に追従しているときの推定電流波形の時間的・空間的パターンを表す（典型的な被験者の単一試行の結果の例）。図５のＦは、視標位置の時間パターンである。図５のＡは、脳表面にのみダイポールを配置して電流推定した場合の時間波形である。図５のＢは、想定されるアーチファクト源にもダイポールを配置して、脳活動とアーチファクト源の振舞いを同時に推定した場合の推定電流波形である。そして、灰色の領域の時刻での平均電流分散の空間的なパターンを図５のＣ，Ｄに示す。
【００９８】
図５のＣ〜Ｅは、アーチファクト源に配置したダイポールの推定電流波形を表す。図５のＤは、心臓のダイポールに関してであるが、心拍に似た特徴をもつ波形が得られた。脳内の推定電流が、アーチファクトによって大きく歪められているならば、脳内とアーチファクトの推定電流の間に何らかの相関があると考えられる。しかし、脳電流は、アーチファクト源に配置したダイポールでの推定電流と脳電流（Ｌａｔｅｒａｌ  Ｏｃｉｐｐｉｔｏ  Ｔｅｍｐｏｒａｌ  Ｃｏｒｔｅｘ，  ＬＯＴＣ）との相関はいずれも小さかった（被験者毎に相関係数を求め、すべての被験者の平均値と標準偏差は，［ｅｙｅ（ｌｅｆｔ）］：（ｒｘ，ｒｙ，ｒｚ）＝（０．００±０．０７，０．００±０．０７，０．０２±０．０８），［ｅｙｅ（ｒｉｇｈｔ）］：（ｒｘ，ｒｙ，ｒｚ）＝（０．００±０．０８，０．００±０．００，０．０±０．００），［ｃａｒｄｉａｃ］：（ｒｘ，ｒｙ，ｒｚ）＝（０．００±０．０８，０．００±０．０８，０．００±０．０７），［ｍａｓｔｉｃａｔｏｒｙ  ｍｕｓｃｌｅ（ｌｅｆｔ）］：（ｒｘ，ｒｙ，ｒｚ）＝（０．００±０．０７，０．００±０．０８，０．００±０．０７）），［ｍａｓｔｉｃａｔｏｒｙ  ｍｕｓｃｌｅ（ｒｉｇｈｔ）］：（ｒｘ，ｒｙ，ｒｚ）＝（０．００±０．０７，０．０１±０．０７，０．００±０．０８））．推定された脳電流がアーチファクトによって歪められている可能性は低いと考えられる。
【００９９】
脳表面だけにダイポールを配置した場合と、アーチファクト源と脳表面にダイポールを配置した場合で、推定結果を比較した。まず、磁場に関して、計測磁場と復元磁場について数式１７で表されるあてはまりの度合（ＧｏＦ）を調べた。
【数１７】
0005266597.tif 000018
【０１００】
２つの条件で求めたＧｏＦを比較した結果、アーチファクト源にダイポールを置くことによって、ＧｏＦは有意に大きくなり（ｐａｉｒｅｄ  ｔ−ｔｅｓｔ（８）＝１５．３２，ｐ＜０．０００００１，図５Ｆ）、値は１２．７２±２．３５％改善した。脳表面に配置したダイポールの数は１８８６〜３２２９個と多数であるのに対して、アーチファクト源に追加したダイポールの数はわずか１５個（［アーチファクト源５箇所］×［３軸方向］）である。少ない数であるが、アーチファクト源に適切に配置することによって、ＧｏＦを効果的に改善することができた。
【０１０１】
次に、推定電流に関して、ｆＭＲＩで活動が期待される位置では、推定電流も大きくなることが予想される。そこで、ｆＭＲＩの活動強度情報と推定電流の平均電流分散の相関関係を調べた。その結果、拡張ダイポールを置くことによって、相関が有意に大きくなり（ｐａｉｒｅｄ  ｔ−ｔｅｓｔ（８）＝２．２８，ｐ＜０．０５，図６のＢ），ｆＭＲＩ活動強度情報と平均電流分散の相関においても改善がみられた。
【０１０２】
以上、本実施の形態における脳活動情報出力装置１２によれば、脳活動とアーチファクトの振舞いを同時に予測し、ノイズ除去を行うことができる。そのため、精度の高い脳活動の情報を出力できる。
【０１０３】
また、本実施の形態によれば、アーチファクト源や脳活動の信号強度は、それぞれ異なることを考慮し、信号強度に関する情報を事前知識として個別に設定し、電流の推定を行うことにより、より精度の高い脳活動の情報を出力できる。
【０１０４】
また、本実施の形態によれば、推定した電流からアーチファクトの磁場を復元し、計測脳信号情報から差し引くことによってアーチファクト除去を行うことができ、精度の高い脳活動の情報を出力できる。
また、本実施の形態によれば、計測脳信号情報は、主として、脳磁図（ＭＥＧ）であった。ただし、計測脳信号情報は、脳波（ＥＥＧ）でも良いことは言うまでもない。かかることは、他の実施の形態においても同様である。
【０１０５】
さらに、本実施の形態によれば、拡張ダイポールを用いてアーチファクト除去を行い、従来の方法と比べて合理的かつ効果的にセンサレベルでノイズ除去ができる。また、アーチファクトも含めて脳電流を推定することによって、電流レベルでのクリーニングもできる。そのため、精度の高い脳活動の情報を出力できる。
【０１０６】
なお、本実施の形態において、上述した数式の情報は、通常、当該数式を用いて演算する部、手段等が、予め保持しており、当該数式を読み出して、実行する。
【０１０７】
さらに、本実施の形態における処理は、ソフトウェアで実現しても良い。そして、このソフトウェアをソフトウェアダウンロード等により配布しても良い。また、このソフトウェアをＣＤ−ＲＯＭなどの記録媒体に記録して流布しても良い。なお、このことは、本明細書における他の実施の形態においても該当する。なお、本実施の形態における脳活動情報出力装置を実現するソフトウェアは、以下のようなプログラムである。つまり、このプログラムは、コンピュータを、記憶媒体に格納されている計測脳信号情報と脳事前情報とを用いて、脳活動の情報である脳活動情報を推定し、取得する脳活動情報取得部と、前記計測脳信号情報と、記憶媒体に格納されている１以上のアーチファクト事前情報とを用いて、アーチファクト源の情報である１以上のアーチファクト情報を推定し、取得するアーチファクト情報取得部と、前記計測脳信号情報と前記脳活動情報と前記１以上のアーチファクト情報とを用いて、前記計測脳信号情報からアーチファクトの影響を取り除く処理を行い、出力される脳活動の情報である出力脳活動情報を取得する出力脳活動情報取得部と、前記出力脳活動情報を出力する出力部として機能させるためのプログラム、である。
【０１０８】
また、上記プログラムにおいて、前記アーチファクト情報取得部は、複数のアーチファクト源ごとに、アーチファクト情報を推定し、複数のアーチファクト情報を取得し、前記出力脳活動情報取得部は、前記計測脳信号情報と前記脳活動情報と前記複数のアーチファクト情報とを用いて、前記計測脳信号情報から複数のアーチファクトの影響を取り除く処理を行い、出力される脳活動の情報である出力脳活動情報を取得するように、コンピュータを、機能させるためのプログラム、であることは好適である。
【０１０９】
また、上記プログラムにおいて、前記脳事前情報、および前記１以上の各アーチファクト事前情報は、信号の強度に関する情報である強度情報を含み、前記脳事前情報の強度情報と、前記１以上の各アーチファクト事前情報の各強度情報は、異なることは好適である。
【０１１０】
また、上記プログラムにおいて、前記１以上のアーチファクト事前情報は、被験者の生体に含まれるアーチファクト源に関する事前の情報である生体アーチファクト事前情報を含むことは好適である。
【０１１１】
また、上記プログラムにおいて、前記１以上のアーチファクト事前情報は、被験者の生体に含まない外部のアーチファクト源に関する事前の情報である外部アーチファクト事前情報を含むことは好適である。
【０１１２】
また、上記プログラムにおいて、前記脳事前情報は、ｆＭＲＩ強度情報を電流分散とし、当該電流分布に対する事前知識を表す事前分布であり、前記１以上の各アーチファクト事前情報は、ｆＭＲＩ強度情報と比較して大きな電流分散を用い、かつ当該電流分布に対する事前知識を表す事前分布であり、前記脳活動情報取得部は、複数のダイポール位置での電流についての情報である脳活動情報を、前記計測脳信号情報と前記脳事前情報とを用いた階層ベイズ推定を用いて推定し、取得し、前記アーチファクト情報取得部は、アーチファクト源のダイポール位置での電流についての情報である１以上のアーチファクト情報を、前記計測脳信号情報と、前記１以上のアーチファクト事前情報とを用いた階層ベイズ推定を用いて推定し、取得し、前記出力脳活動情報取得部は、前記計測脳信号情報から、前記脳活動情報のうちの脳の活動が無いと考えられる位置に配置されたダイポールに対応する脳活動情報と前記アーチファクト情報とを除く処理を行い、出力脳活動情報を取得するように、コンピュータを、機能させるためのプログラム、であることは好適である。
【０１１３】
また、上記プログラムにおいて、前記計測脳信号情報は、頭部周辺で観測された磁場に関する情報であり、前記出力脳活動情報取得部は、前記脳活動情報のうちの脳の活動が無いと考えられる位置に配置されたダイポールに対応する脳活動情報である電流の情報から磁場に関する情報である背景脳磁場情報を取得し、かつ、前記アーチファクト情報である電流の情報から磁場に関する情報であるアーチファクト磁場情報を取得し、前記計測脳信号情報から、前記背景脳磁場情報と前記アーチファクト磁場情報とを除く処理を行い、出力脳活動情報を取得するように、コンピュータを、機能させるためのプログラム、であることは好適である。
【０１１４】
さらに、上記プログラムにおいて、前記出力脳活動情報取得部が取得した出力脳活動情報と、前記計測脳信号情報とを用いて、当該出力脳活動情報におけるノイズ除去の精度に関する情報であるノイズ除去精度情報を取得するノイズ除去精度情報取得部をさらに具備し、前記出力部は、前記ノイズ除去精度情報が、予め格納されている閾値と比較して高い精度である場合のみ、前記出力脳活動情報を出力するように、コンピュータを、機能させるためのプログラム、であることは好適である。
【０１１５】
また、図７は、本明細書で述べたプログラムを実行して、上述した実施の形態の脳活動情報出力装置等を実現するコンピュータの外観を示す。上述の実施の形態は、コンピュータハードウェア及びその上で実行されるコンピュータプログラムで実現され得る。図７は、このコンピュータシステム３４０の概観図であり、図８は、コンピュータシステム３４０のブロック図である。
【０１１６】
図７において、コンピュータシステム３４０は、ＦＤドライブ３４１１、ＣＤ−ＲＯＭドライブ３４１２を含むコンピュータ３４１と、キーボード３４２と、マウス３４３と、モニタ３４４とを含む。
【０１１７】
図８において、コンピュータ３４１は、ＦＤドライブ３４１１、ＣＤ−ＲＯＭドライブ３４１２に加えて、ＭＰＵ３４１３と、ＣＤ−ＲＯＭドライブ３４１２及びＦＤドライブ３４１１に接続されたバス３４１４と、ブートアッププログラム等のプログラムを記憶するためのＲＯＭ３４１５と、ＭＰＵ３４１３に接続され、アプリケーションプログラムの命令を一時的に記憶するとともに一時記憶空間を提供するためのＲＡＭ３４１６と、アプリケーションプログラム、システムプログラム、及びデータを記憶するためのハードディスク３４１７とを含む。ここでは、図示しないが、コンピュータ３４１は、さらに、ＬＡＮへの接続を提供するネットワークカードを含んでも良い。
【０１１８】
コンピュータシステム３４０に、上述した実施の形態の脳活動情報出力装置等の機能を実行させるプログラムは、ＣＤ−ＲＯＭ３５０１、またはＦＤ３５０２に記憶されて、ＣＤ−ＲＯＭドライブ３４１２またはＦＤドライブ３４１１に挿入され、さらにハードディスク３４１７に転送されても良い。これに代えて、プログラムは、図示しないネットワークを介してコンピュータ３４１に送信され、ハードディスク３４１７に記憶されても良い。プログラムは実行の際にＲＡＭ３４１６にロードされる。プログラムは、ＣＤ−ＲＯＭ３５０１、ＦＤ３５０２またはネットワークから直接、ロードされても良い。
【０１１９】
プログラムは、コンピュータ３４１に、上述した実施の形態の脳活動情報出力装置等の機能を実行させるオペレーティングシステム（ＯＳ）、またはサードパーティープログラム等は、必ずしも含まなくても良い。プログラムは、制御された態様で適切な機能（モジュール）を呼び出し、所望の結果が得られるようにする命令の部分のみを含んでいれば良い。コンピュータシステム３４０がどのように動作するかは周知であり、詳細な説明は省略する。
【０１２０】
なお、上記プログラムにおいて、情報を送信する送信ステップや、情報を受信する受信ステップなどでは、ハードウェアによって行われる処理、例えば、送信ステップにおけるモデムやインターフェースカードなどで行われる処理（ハードウェアでしか行われない処理）は含まれない。
【０１２１】
また、上記プログラムを実行するコンピュータは、単数であってもよく、複数であってもよい。すなわち、集中処理を行ってもよく、あるいは分散処理を行ってもよい。
【０１２２】
また、上記各実施の形態において、各処理（各機能）は、単一の装置（システム）によって集中処理されることによって実現されてもよく、あるいは、複数の装置によって分散処理されることによって実現されてもよい。
【０１２３】
本発明は、以上の実施の形態に限定されることなく、種々の変更が可能であり、それらも本発明の範囲内に包含されるものであることは言うまでもない。
【産業上の利用可能性】
【０１２４】
以上のように、本発明にかかる脳活動情報出力装置は、精度の高い脳活動の情報を出力できる、という効果を有し、脳活動情報出力装置等として有用である。
【図面の簡単な説明】
【０１２５】
【図１】実施の形態１における脳活動情報出力システムのブロック図
【図２】同脳活動情報出力システムの動作について説明するフローチャート
【図３】同アーチファクト除去前後でのＭＥＧ波形の比較を示す図
【図４】同観測値、ＰＣＡ、ＶＢＭＥＧを用いた方法での平均からのずれを示す図
【図５】同推定電流の時間波形と空間マップを示す図
【図６】同拡張ダイポールの有無での推定結果の違いを示す図
【図７】同コンピュータシステムの概観図
【図８】同コンピュータシステムのブロック図
【符号の説明】
【０１２６】
１１  脳信号計測装置
１２  脳活動情報出力装置
１２１  計測脳信号情報格納部
１２２  脳事前情報格納部
１２３  アーチファクト事前情報格納部
１２４  脳活動情報取得部
１２５  アーチファクト情報取得部
１２６  出力脳活動情報取得部
１２７  ノイズ除去精度情報取得部
１２８  出力部
【図２】
0005266597.tif 000019
【図７】
0005266597.tif 000020
【図８】
0005266597.tif 000021
【図１】
0005266597.tif 000022
【図３】
0005266597.tif 000023
【図４】
0005266597.tif 000024
【図５】
0005266597.tif 000025
【図６】
0005266597.tif 000026