(57)【要約】【課題】従来の脳活動部位の推定法では、簡便であるが、電流双極子の位置および方向について妥当性が高くないという問題があった。【解決手段】脳のＭＲＩデータを用いて電流双極子の存在するべき領域および方向を予め限定して脳の活動部位を推定することにより、電流双極子の位置を正確に、かつ効率よく推定することが可能となる。また、その電流双極子の位置および方向の推定方法を用いて、正確な脳活動部位を推定することが可能となる。脳のＭＲＩデータを用いて脳の３次元輪郭を抽出するため、個々の脳の形に適した電流双極子を推定することが可能である。【選択図】    図１
2006061321.tif 000002
【特許請求の範囲】
【請求項１】
脳の神経活動の有無を示す電流双極子の領域を大脳皮質に限定し、前記電流双極子の方向を前記大脳皮質の表面に対して垂直に限定することを特徴とする電流双極子の位置および方向の推定方法。
【請求項２】
脳の画像を取得して脳表画素を特定するステップと、
脳表から脳外に向かう単位法線ベクトルを算出するステップと、
を具備することを特徴とする電流双極子の位置および方向の推定方法。
【請求項３】
脳の画像を取得するステップと、
前記画像を、脳内と脳外の境界を含む画素群に分割するステップと、
前記画素群から脳表画素を特定するステップと、
前記脳表画素を中心とした行×列×段の画素からなる立方画素群を抽出し、前記立方画素群のうち脳外の画素を選択するステップと、
前記脳表画素から前記立方画素群内の前記脳外の画素に対する少なくとも１つの単位ベクトルを算出するステップと、
前記各単位ベクトルの和を正規化して前記脳表画素の単位法線ベクトルを算出し、電流双極子の位置および方向として推定するステップと、
を具備することを特徴とする電流双極子の位置および方向の推定方法。
【請求項４】
前記脳表画素は、前記脳外の画素に隣接する脳内の画素で構成されることを特徴とする請求項３に記載の電流双極子の位置および方向の推定方法。
【請求項５】
前記立方画素群は３行×３列×３段の２７画素からなることを特徴とする請求項３に記載の電流双極子の位置および方向の推定方法。
【請求項６】
脳の神経活動により発生する電流源の位置を大脳皮質に限定し、前記電流源の方向を前記大脳皮質の表面に対して垂直方向に限定することにより脳活動部位を推定することを特徴とする脳活動部位推定方法。
【請求項７】
被験者の刺激部位に電気刺激を与えるステップと、
前記電気刺激により発生する脳磁界もしくは脳電位の計測を行うステップと、
予め取得した脳の画像を基に、脳表画素ごとの単位法線ベクトルを算出するステップと、
前記単位法線ベクトルを用いて電流源の値を算出するステップと、
前記電流源の値の分布に基づき脳活動部位を推定するステップと、
を具備することを特徴とする脳活動部位推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、脳の神経活動の有無を示す電流双極子の存在するべき領域および方向を限定した電流双極子の位置および方向の推定方法、およびそれを用いた脳活動部位推定方法に関するものである。
【背景技術】
【０００２】
脳における活動部位は、脳の活動に伴い頭皮上に現れる電位を電極で計測した脳電位から推定することができる。前記脳電位測定を行うために、電流双極子の位置を推定する必要がある。推定方法は、現在脳を球として近似したものや、有限要素に分解したものを用いて電流双極子の位置を推定する方法が多用されている。
【０００３】
図１４に従来の電流双極子の推定方法を示す。図１４は実形状４層頭蓋モデルの斜視図である。図示のごとく、４層頭蓋モデルは脳質２１、脳髄液２２、頭蓋骨２３、頭皮２４の４層領域よりなる。脳質２１、脳髄液２２、頭蓋骨２３、頭皮２４の領域をそれぞれΩ１、Ω２、Ω３、Ω４としたとき、電流双極子は領域Ω１内にのみ発生する。
【０００４】
また、脳質２１、脳髄液２２、頭蓋骨２３、頭皮２４の導電率をそれぞれσ１、σ２、σ３、σ４として、境界要素法を用いて行列連立方程式を解く。次に、伝達行列の計算を行い、その結果を用いて電極位置の電位を求め、電流双極子の真の位置とベクトル成分を推定する（例えば特許文献１参照。）。
【特許文献１】特開平８−２５７００４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
従来の脳活動部位の推定方法では、脳を球として近似したものや有限要素に分解したものを用いて電流双極子の位置を推定する場合が多いため、簡便であるという利点を有する。
【０００６】
しかしながら、電流双極子は神経細胞の位置に存在し、その方向は神経細胞の走行方向に出現しなければならず、前記推定方法ではこの要件を満たしていない。そのため、電流双極子の正確な位置の推定が困難であるだけでなく、本来電流双極子が存在するはずのない部位に電流双極子を推定してしまう場合があるという問題があった。
【０００７】
また、一般的に脳電位から推定される脳の活動部位は、電流双極子として推定されるが、脳の神経細胞の並び方を考慮していない推定法では電流双極子の方向を予め制限していないため計算量が膨大となる。さらに、脳の形は人によってそれぞれ異なるため、各々の脳の形に適した電流双極子の位置を推定することが困難であるという問題があった。
【課題を解決するための手段】
【０００８】
本発明は係る課題に鑑みてなされ、脳の神経活動の有無を示す電流双極子の領域を大脳皮質に限定し、前記電流双極子の方向を前記大脳皮質の表面に対して垂直に限定することを特徴とするものである。
【０００９】
また、脳の画像を取得するステップと、
前記画像を、脳内と脳外の境界を含む画素群に分割するステップと、
前記画素群から脳表画素を特定するステップと、
前記脳表画素を中心とした行×列×段の画素からなる立方画素群を抽出し、前記立方画素群のうち脳外の画素を選択するステップと、
前記脳表画素から前記立方画素群内の前記脳外の画素に対する少なくとも１つの単位ベクトルを算出するステップと、
前記各単位ベクトルの和を正規化して前記脳表画素の単位法線ベクトルを算出し、電流双極子の位置および方向として推定するステップと、
を具備することを特徴とするものである。
【００１０】
また、前記脳表画素は、前記脳外の画素に隣接する脳内の画素で構成されることを特徴とするものである。
【００１１】
また、前記立方画素群は３行×３列×３段の２７画素からなることを特徴とするものである。
【００１２】
また、脳の神経活動により発生する電流源の位置を大脳皮質に限定し、前記電流源の方向を前記大脳皮質の表面に対して垂直方向に限定することにより脳活動部位を推定することを特徴とするものである。
【００１３】
また、被験者の刺激部位に電気刺激を与えるステップと、
前記電気刺激により発生する脳磁界もしくは脳電位の計測を行うステップと、
予め取得した脳の画像を基に、脳表画素ごとの単位法線ベクトルを算出するステップと、
前記単位法線ベクトルを用いて電流源の値を算出するステップと、
前記電流源の値の分布に基づき脳活動部位を推定するステップと、
を具備することを特徴とするものである。
【発明の効果】
【００１４】
本発明によれば、第１に、電流双極子の領域を大脳皮質に限定し、かつ電流双極子の方向を大脳皮質の表面（脳表）に対して垂直方向に限定することにより、電流双極子の位置および方向を効率よくかつ正確に推定することが可能となる。
【００１５】
第２に、脳のＭＲＩデータを用いて脳の３次元輪郭を抽出するため、脳を球として近似する必要がなく、個々の脳の形に適した電流双極子を推定することが可能となる。
【００１６】
第３に、対象となる画素を中心とした立方画素群から脳表画素を特定することによって正確な脳表の位置を特定することができる。
【００１７】
第４に、脳表から脳外に向かう単位法線ベクトルを求めるため、脳表に対する垂直方向を正確に得ることができる。
【００１８】
第５に、対象となる脳表画素の単位法線ベクトルは、その脳表画素に隣接した画素群、すなわち脳表画素を中心とした３行×３列×３段の立方画素群を用いて求めることができる。つまり、必要最小限の画素群からなる立方画素群により脳表画素の単位法線ベクトルを求めることができるため、それにかかる計算量が最小ですむ。
【００１９】
第６に、電流源の位置を大脳皮質に限定し、かつ電流源の方向を大脳皮質の表面（脳表）に対して垂直方向に限定することにより、脳活動部位の位置および方向を効率よく正確に推定することが可能となる。
【発明を実施するための最良の形態】
【００２０】
図１から図６を参照して本発明の第１の実施の形態について詳細に説明する。
【００２１】
本発明の実施の形態では、脳の神経活動の有無を示す電流双極子の領域を大脳皮質に限定し、その方向を大脳皮質の表面に対して垂直に限定して電流双極子の位置および方向の推定を行う。
【００２２】
図１は、電流双極子の位置および方向の推定方法のフローチャート図である。
【００２３】
本発明の電流双極子の位置および方向の推定方法は、脳の画像を取得するステップと、前記画像を、脳内と前記脳外の境界を含む画素群に分割するステップと、前記画素群から脳表画素を特定するステップと、前記脳表画素を中心とした行×列×段の画素からなる立方画素群を抽出し、前記立方画素群のうち脳外画素を選択するステップと、前記脳表画素から前記立方画素群内の脳外画素に対する少なくとも１つの単位ベクトルを算出して、それを電流双極子の位置および方向として推定するステップと、前記各単位ベクトルの和を正規化して前記脳表画素の単位法線ベクトルを算出するステップと、から構成される。
【００２４】
第１のステップは、脳の画像を取得するステップである（ステップＳ１）。
【００２５】
まず、ＭＲＩ（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ｉｍａｇｉｎｇ：磁気共鳴画像）により被験者のＭＲ（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ）画像を取得する。図２は取得した１枚のＭＲ画像データ図である。図示のごとく、ＭＲ画像１はグレースケールの情報をもつ２次元画像で構成されている。脳の３次元画像は、ＭＲ画像１のような頭部の水平断層像を約１ｍｍ間隔で撮影し、連続的に取得することによって構築される。なお、図２、図４から図６、図１１および図１２（後述）では、脳外を黒色、脳内を白色および灰色でそれぞれ表示している。
【００２６】
第２のステップは、前記画像を、脳内と脳外の境界を含む画素群に分割するステップである（ステップＳ２）。
【００２７】
脳磁界の発生源となるのは主に大脳皮質の表面に対して垂直に配列している大脳皮質内の錐体細胞の活動であると考えられている。
【００２８】
図３は、大脳皮質の概要図である。図３（Ａ）は大脳皮質４内に流れる電流双極子の概要図である。大脳皮質４においては多数の錐体細胞が存在することが知られている。この錐体細胞が集団的に活動すれば、その活動電流１３の流れは図３（Ａ）に示すように大脳皮質４の表面、すなわち脳表１２に向かって垂直な方向となる。この時間的に同期して活動する神経細胞群の電流は集合すると、まとまった値となる。
【００２９】
図３（Ｂ）は電流から発生する磁界の概要図である。図示のごとく、大脳皮質４内の神経細胞群に集中した一方向の電流が流れれば、その周りに右ねじの法則に従って磁界１５が発生する。つまり、脳磁界は主に脳表１２に対して垂直に配列した大脳皮質４内の錐体細胞の活動により生じるものだと考え、本実施形態では電流源の位置を大脳皮質４に、方向を脳表１２に対して垂直な方向に限定する。
【００３０】
そのためにまず脳表１２の特定を行う。脳表１２はＭＲ画像を基に脳内および脳外を判別することにより特定する。従って、まずＭＲ画像を脳内および脳外を含む画素群に分割する。
【００３１】
図４は図２の一部の画素群２を模式的に表した拡大図である。本実施の形態では、ＭＲＩにより取得したＭＲ画像１（図２参照）のような水平断層像を一辺の長さが約１ｍｍの正方形の画素ごとに分割する。そしてその操作を全てのＭＲ画像に対して行う。そのＭＲ画像の各画素のグレースケールの濃淡により、脳外および脳内が判別できる。
【００３２】
第３のステップは前記画素群から脳表画素を特定するステップである（ステップＳ３）。
【００３３】
脳表画素の特定は、任意の脳内画素を選択し、その脳内画素の周囲にある画素により、その脳内画素が脳表１２であるか否かを判断する。脳表１２を特定するためには任意の脳内画素に対して上下の画素にも着目する必要がある。そのため本実施の形態では、上下に間隔約１ｍｍで連続する３枚のＭＲ画像を用いて脳表画素の特定を行う。
【００３４】
図５は上下に連続する３枚のＭＲ画像を拡大した画像データ図である。画素群Ａ、Ｂ、Ｃは上下に連続する３枚のＭＲ画像の一部分であり、例えば任意の６×６の画素からなる画素群から構成されている。その画素群Ａ、Ｂ、Ｃの中から任意の脳内画素を選択し、前記脳内画素の周囲にある計２６個の画素のうち１画素でも黒色の画素が存在すると、その画素を脳表１２とする。
【００３５】
例えば、図５に示す画素群Ｂの５行４列目の脳内画素（Ｂ［４］［３］と表記する。以下同様）を選択すると、その脳内画素と隣接する画素は、画素群ＡでＢ［４］［３］の上方に位置するＡ［４］［３］を中心とした３×３（＝９）画素と、画素群ＢでＢ［４］［３］を中心とした３×３（＝９）画素のうち、Ｂ［４］［３］を除いた８画素と、画素群ＣでＢ［４］［３］の下方に位置するＣ［４］［３］を中心とした３×３（＝９）画素の合計２６画素である。そのうちＡ［３］［２］、Ａ［３］［３］、Ａ［４］［４］、Ｂ［３］［３］、Ｂ［３］［４］、Ｂ［４］［４］、Ｃ［３］［３］、Ｃ［４］［４］が黒色であるため、Ｂ［４］［３］は脳表１２と見なすことができる。
【００３６】
第４のステップは、前記脳表画素を中心とした行×列×段の画素からなる立方画素群を抽出し、前記立方画素群のうち脳外画素を選択するステップである（ステップS４）。
【００３７】
本ステップでは脳表画素に隣接する画素のうち脳外画素を選択するため、ステップＳ３で特定された脳表画素を中心として、立方画素群１４を抽出する。すなわち、立方画素群は、図５に示すごとく、Ｂ［４］［３］を中心とした３×３×３（＝２７）画素の立方画素群１４が抽出される。立方画素群１４の中で、Ａ［３］［２］、Ａ［３］［３］、Ａ［４］［４］、Ｂ［３］［３］、Ｂ［３］［４］、Ｂ［４］［４］、Ｃ［３］［３］、Ｃ［４］［４］が黒色であるため、これらの画素が１つの脳表画素に対する脳外画素として特定される。
【００３８】
第５のステップは、前記脳表画素から前記立方画素群内の脳外画素に対する少なくとも１つの単位ベクトルを算出して、それを電流双極子の位置および方向として推定するステップである（ステップＳ５）。
【００３９】
前述したように、脳磁界の発生源となるのは主に大脳皮質４の表面に対して垂直に配列している大脳皮質４内の錐体細胞の活動であると考えられている。そこで、本実施形態では、電流双極子の方向を錐体細胞が神経を伸ばしている方向、すなわち大脳皮質４の表面に対して垂直な方向に限定する。そのために、脳表画素に隣接する全ての脳外画素を求め（ステップＳ４）、本ステップでは脳表画素から全ての脳外画素へのベクトルを算出する。
【００４０】
図６は脳表画素から脳外画素へのベクトルを示す図である。
【００４１】
前記脳表画素からステップＳ４で選択された脳外画素へのベクトル群５をそれぞれ計算し、各ベクトル方向の単位ベクトル群を算出する。
【００４２】
第６のステップは、前記各単位ベクトルの和を正規化して前記脳表画素の単位法線ベクトルを算出するステップである（ステップＳ６）。
【００４３】
図７は、脳表画素Ｂ［４］［３］の単位法線ベクトルを示す図である。単位法線ベクトル６は、算出した各単位ベクトル群の和ベクトルを求め、この和ベクトル方向の単位ベクトルとして求める。
【００４４】
この手順により、全ての脳表画素に対して単位法線ベクトルを算出する。この単位法線ベクトルの位置および方向が電流双極子の位置および方向となる。
【００４５】
上述した手順を行うことにより、全ての脳表１２に対する電流双極子の位置および方向を推定する。
【００４６】
一般に、測定した電界分布、磁界分布からその発生源となる電流双極子の位置、方向、大きさを推定することは逆問題を解くことになる。しかしながら、逆問題の解は一意に決定されず、計算量が膨大となるだけでなく、本来電流双極子が存在するはずのない部位に電流双極子を推定してしまうことがあるため、電流双極子の位置および方向を正確に推定することが困難であった。
【００４７】
しかしながら、上述したように電流双極子の位置を大脳皮質４に、方向を脳表１２に対して垂直方向に限定することにより、従来の電流双極子推定方法と比較して計算負荷が軽減され、電流双極子の位置および方向を効率よくかつ正確に推定することができる。
【００４８】
次に、本発明の第２の実施の形態について詳細に説明する。
【００４９】
本実施の形態では、第１の実施の形態により求めた電流双極子の推定方法を用いて、電流源の位置および方向を推定することによって脳活動部位の推定を行う。
【００５０】
脳の活動により発生する脳磁界は、電流双極子の位置が発生源となり、電流双極子の方向に発生する。その脳の活動により発生した脳磁界の発生源が電流源となる。つまり、電流源の位置は電流双極子の位置に発生し、その方向は電流双極子の方向となる。そのため、電流源の位置および方向を電流双極子の位置および方向に限定して脳活動部位を推定する。
【００５１】
電流源の値の推定は、上記で算出した電流双極子の列ベクトルを用いて行う。
【００５２】
電流源の値を推定する式を求める方法として、脳内を空間離散化し、全ての格子点に電流要素をあらかじめ配置して、それらの強度分布を推定することで脳内の神経活動を解析することができるＬ２ノルム法を用いた。
【００５３】
Ｍ個のセンサーによって得られた計測磁界ベクトルｂは（１）式によって表される。
【００５４】
【数１】
2006061321.tif 000003
また、各格子点における実際の電流双極子の列ベクトルｐは（２）式によってあらわされる。
【００５５】
【数２】
2006061321.tif 000004
また、Ｌを脳磁計のセンサーと脳内の各格子点における電流双極子との位置関係によって一意に決まるＭ×Ｎの係数行列とすると、脳磁計によって頭表面で計測される脳磁界は（３）式によって表すことができる。
【００５６】
【数３】
2006061321.tif 000005
従って、電流源ベクトルｐ´は（４）式で表すことができる。
【００５７】
【数４】
2006061321.tif 000006
しかし、一般にＬは正則行列でないため、逆行列Ｌ−１は存在しない。そこで、（４）式に対して解の二乗ノルムの最小化をはかり、逆行列Ｌ−１を求める目的関数である（５）式を用いると、逆行列Ｌ−１は（６）式によって表すことができる。
【００５８】
【数５】
2006061321.tif 000007
電流源は（４）式と（６）式より求めることができる。しかし、実測ではノイズの影響を考慮する必要がある。そこで、γを電流の推定値に及ぼすノイズの影響を制御するパラメータとし、ＩをＭ×Ｍの単位行列としてチホノフの適切化を用いると、逆行列Ｌ−１は（７）式のように表すことができる。
【００５９】
【数６】
2006061321.tif 000008
プログラムでは、γの値を計測した磁界データからの信号成分の分散値とノイズ成分の分散値により決定している。
【００６０】
（４）式と（７）式より、導かれる電流源の値は（８）式によって求めることができる。
【００６１】
【数７】
2006061321.tif 000009
上述した手法により電流源の位置および方向を特定し、電流源の値を求める。そして、電流源の大きな値を脳活動部位として推定できる。以下、この手法の妥当性について検証を行う。
【００６２】
図８から図１３を参照して本推定方法の妥当性についての検証のための具体的な実験方法を説明する。本実験では、脳磁界計測および脳電位計測により脳活動部位の推定を行い、本発明の妥当性について検証を行う。
【００６３】
まず、脳磁界計測による脳活動部位の推定について説明する。脳磁界は、神経活動によって発生する電流源の周囲に発生する磁界であり、その位置および値は刺激部位によって変化する。
【００６４】
図８は、脳磁界計測に用いる脳磁界計測システムの構成図である。
【００６５】
図８（Ａ）はＮｅｕｒｏｍａｇ社製の脳磁界計測システムＶｅｃｔｏｒｖｉｅｗ（登録商標）である。図示のごとく、脳磁界計測システム３は被験者の頭部を覆い、脳磁界の値の検出を行うデュワー部８を備える。
【００６６】
図８（Ｂ）はデュワー部８の拡大図である。デュワー部８は、液体ヘリウムが充填されたデュワー７とセンサー部９とから構成されている。また、デュワー７には図示しない磁界検出コイルが配置されており、脳から湧き出す磁界をほぼ頭部に垂直な向きで検出する。センサー部９は伝導量子干渉素子（ＳＱＵＩＤ：ｓｕｐｅｒｃｏｎｄｕｃｔｉｎｇ  ｑｕａｎｔｕｍ  ｉｎｔｅｒｆｅｒｅｎｃｅ  ｄｅｖｉｃｅ）を用いているため、微弱な脳磁界の値を検出することができる。
【００６７】
図８（Ｃ）は、センサー部９のセンサー配置図である。図示のごとく、２０４チャネルのコイル１０が被験者の頭部全体を覆うように配置されており、全頭部における磁界を同時に計測することができる。
【００６８】
被験者は２２〜２６歳の健常な男性３名とした。刺激部位は左右の手首とした。刺激電極は手首に装着し、動かないように固定した。
【００６９】
図９は、被験者に与える電気刺激に使用する刺激パルス図である。図示のごとく、刺激装置から送られる刺激パルスは刺激時間を５ｍｓとし、刺激間隔を５００ｍｓとした。
【００７０】
なお、被験者によって刺激の感じ方が異なるため、刺激強度は６〜８ｍＡに設定した。また、電気刺激によって脳から発生する磁界は非常に小さいため、頭皮上から記録すると様々なノイズの影響を受ける。そのため、刺激によって誘発された信号成分を検出することが困難である。そこで、刺激開始時にトリガー１１をかけて反復計測を行ない、得られた磁界データを加算平均する加算平均法を用いて磁界データを得る。本実験では、加算平均は約２００回以上行い、その結果を解析データとして用いた。
【００７１】
図１０は、脳磁界計測および脳電位計測による脳活動部位推定を行うフローチャート図である。
【００７２】
第１のステップは、被験者の刺激部位に電気刺激を与えるステップである（ステップＳ１１）。上述したように、本実施の形態では被験者の刺激部位を左右の手首とした。
【００７３】
第２のステップは、前記電気刺激により発生する脳磁界もしくは脳電位の計測を行うステップである（ステップＳ１２）。本実験では脳磁界の計測を行う。脳磁界の値は脳磁界計測システム３（図８（Ａ）参照）を用いて計測する。
【００７４】
第３のステップは、予め取得したＭＲ画像（図１のステップＳ１）を基に、脳表画素ごとの単位法線ベクトルを算出するステップである（ステップＳ１３）。上述したように、脳内画素を中心とした３行×３列×３段の立方画素群のうち、１画素でも脳外画素が存在した場合、その脳内画素を脳表画素と特定する（図１のステップＳ２からステップＳ４）。
【００７５】
また、単位法線ベクトルは脳表画素から立方画素群中の全ての脳外画素への単位ベクトルの和ベクトルを正規化して算出する（図１のステップＳ５からステップＳ６）。
【００７６】
第４のステップは、前記単位法線ベクトルを用いて電流源の値を算出するステップである（ステップＳ１４）。本実施形態では推定アルゴリズムとしてＬ２ノルム法を用いる。電流源の値は、式（３）において、ベクトルｂにステップＳ１２で計測した脳磁界データを代入し、ＬにステップＳ１３で算出した各脳表画素の単位法線ベクトルから計算される各センサーでの磁界を代入することによって求めることができる。
【００７７】
第５のステップは、前記電流源の値の大きな値となる部位を脳活動部位とするステップである（ステップＳ１５）。ステップＳ１４により算出した電流源の値のうち、大きな値を示した電流源の部位を脳活動部位として推定する。
【００７８】
以上の手順により、脳磁界計測による脳活動部位の推定を行う。
【００７９】
図１１は、各被験者の右手首刺激における電流源推定結果である。図１１の右図は左図の四角で囲まれている脳を拡大した図である。本実施の形態では一例として０から電流源の最大値までを１０段階に等分割し、各段階を色分けして表示することにした。そして、他の部位より電流源が大きい値の部位を脳活動部位として推定する。
【００８０】
図１１の右図の丸で囲まれている部分は例えば８段階以上の電流源の値が推定された場所であり、これを本実施形態による脳活動部位として推定する。図示のごとく、全ての被験者において、この領域は左半球の中心溝と中心後溝に挟まれているため体性感覚野に確当する部分であるということがわかる。
【００８１】
この電流源推定結果の脳の左半球部分と図１２に示すペンフィールドの脳の左半球部分の機能局在図を比較すると、本発明による脳活動部位推定方法では、左半球の手首の刺激による電流源が推定されており、ペンフィールドの脳の機能局在図とほぼ一致する結果となった。ゆえに、本実施形態の脳活動部位推定方法による右手首刺激における脳活動部位は正確に推定することができる。
【００８２】
また、左手首刺激における脳活動部位推定を行った結果、全ての被験者において電流源が体性感覚野に推定されていたことから、左手首刺激における脳活動部位推定結果も正確に推定することができる。
【００８３】
図１３は、一人の被験者に対して行った右足首、右手首、右手中指、右手人差し指、右手親指の刺激における脳活動部位推定結果である。図１３は、３次元化した脳の左半球部分を後方から示している。図１３中の星形で囲まれた部分が右足首刺激における脳活動部位推定結果であり、三角で囲まれた部分が右手首刺激における脳活動部位推定結果であり、四角で囲まれた部分が右手中指刺激における脳活動部位推定結果であり、五角で囲まれた部分が右手人差し指刺激における脳活動部位推定結果であり、丸で囲まれた部分が右手親指刺激における脳活動部位推定結果である。いずれも電流源の値が８段階以上である等、他の部位より大きい値を示した部位である。
【００８４】
この電流源推定結果の脳の左半球部分と図１２に示すペンフィールドの脳の左半球部分の機能局在図を比較すると、本発明による脳活動部位推定方法では、大脳の中央を縦に走る大脳縦列付近の体性感覚野に足首刺激による電流源をはじめ、体性感覚野を下に向かって手首、中指、人差し指、親指刺激による電流源が推定されており、ペンフィールドの脳の機能局在図とほぼ一致する結果となった。ゆえに、本実施形態の脳活動部位推定方法による右足首、右手首、右手中指、右手人差し指、右手親指刺激における脳活動部位推定結果は妥当であると考えられる。
【００８５】
このことより、脳磁界計測における本実施形態の脳活動部位推定方法は、脳活動部位を正確に推定できるといえる。
【００８６】
次に、脳電位計測による脳活動部位の推定について説明を行う。脳電位は神経活動によって電流源が発生する際、媒質表面である頭皮の皮膚上に分布する電位であり、その位置および値は刺激部位によって変化する。
【００８７】
脳電位計測による脳活動部位の推定は、脳磁界計測による脳活動部位の推定と同様の手順で行う。（図１０参照）。
【００８８】
体性感覚誘発電位の計測は２１チャネルの電極を有する脳波計を用いて行った。被験者は眼瞼を閉じた状態で計測を行う。電気刺激は手首正中神経および手指部に与え、刺激時間を５ｍｓとし、刺激間隔を５００ｍｓとした。
【００８９】
刺激部位に電気刺激を与え（ステップＳ１１）、それによって発生する体性感覚誘発電位の計測を行う（ステップＳ１２）。
【００９０】
次に、本発明の推定方法により脳表１２の画素ごとに単位法線ベクトルを求め（ステップＳ１３）、ステップＳ１２により計測したデータを基に、電流源の方向に垂直な方向に限定した電流源の値を算出する（ステップＳ１４）。電流源の値は、式（３）において、ベクトルｂにステップＳ１２で計測した体性感覚誘発電位の計測データを代入し、ＬにステップＳ１３で算出した各脳表画素の単位法線ベクトルから計算される各電極での電位を代入することによって求めることができる。
【００９１】
そして、０から求めた電流源の値の最大値までを１０段階に等分割し、他の部位より大きい値が測定された場所を脳活動部位として推定する。（ステップＳ１５）。
以上の手順により推定を行った脳電位計測による脳活動部位と、図１２に示すペンフィールドの脳の機能局在図を比較すると、ペンフィールドの脳の機能局在図とほぼ一致する結果となった。
【００９２】
このことより、脳電位計測における本実施形態の脳活動部位推定方法は、脳活動部位を正確に推定することができる。
【００９３】
以上より、本発明の脳活動部位推定方法は、脳磁界計測および脳電位計測のどちらの計測方法を用いても脳活動部位を高い精度で推定することが可能である。
【００９４】
本実施形態では本推定方法の妥当性を検証するため、一例として電流源の最大値までの値を１０段階に等分割し、そのうち他の部位より大きい部位（例えば８段階以上）を脳活動部位として推定し説明をした。しかしこれに限らず、電流源の大きさが視覚的に判別できればよく、他の部位より大きい値を示した部位を活動部位とする。このように、本推定法によれば、電流源の大きさの分布を測定することにより、脳の活動の状態を調べることができる。
【図面の簡単な説明】
【００９５】
【図１】本発明の電流双極子の位置および方向の推定方法を示すフローチャート図である。
【図２】本発明の１枚の脳のＭＲ画像である。
【図３】本発明の電流双極子の位置および方向について説明する概要図である。
【図４】本発明の脳表画素の分割について説明する画像データ図である。
【図５】本発明の脳表画素を中心とする立方画素群を抽出した画像データ図である。
【図６】本発明の脳表画素から脳外画素へのベクトルを算出する方法について説明する画像データ図である。
【図７】本発明の電流源の方向の算出方法について説明する画像データ図である。
【図８】本発明の脳磁界の計測を行うための脳磁界計測システムの構成図である。
【図９】本発明の電気刺激に使用する刺激パルス図である。
【図１０】本発明の脳活動部位推定を行うフローチャート図である。
【図１１】本発明の各被験者の右手首刺激における電流源推定結果を示す図である。
【図１２】ペンフィールドの脳の左半球部分の機能局在図である。
【図１３】本発明の一人の被験者に対して行った右足首、右手首、右手中指、右手人差し指、右手親指の刺激における電流源推定結果を示す脳の左半球部分の図である。
【図１４】従来の実形状４層頭蓋モデルの斜視図である。
【符号の説明】
【００９６】
１      ＭＲ画像
２      画素群
３      脳磁界計測システム
４      大脳皮質
５      ベクトル群
６      単位法線ベクトル
７      デュワー
８      デュワー部
９      センサー部
１０    コイル
１１    トリガー
１２    脳表
１３    活動電流
１４    立方画素群
１５    磁界
２１    脳質
２２    脳髄液
２３    頭蓋骨
２４    頭皮
【図１】
2006061321.tif 000010
【図２】
2006061321.tif 000011
【図３】
2006061321.tif 000012
【図４】
2006061321.tif 000013
【図５】
2006061321.tif 000014
【図６】
2006061321.tif 000015
【図７】
2006061321.tif 000016
【図８】
2006061321.tif 000017
【図９】
2006061321.tif 000018
【図１０】
2006061321.tif 000019
【図１１】
2006061321.tif 000020
【図１２】
2006061321.tif 000021
【図１４】
2006061321.tif 000022
【図１３】
2006061321.tif 000023