(57)【要約】      （修正有）【課題】脳内双極子の位置とベクトル成分とを脳波計測に対応して実時間で高速解析処理する脳波双極子解析装置を提供する。【解決手段】本発明の脳波双極子解析装置は、脳波計測用の複数の電極２と、各電極２の計測電位を増幅する多チャンネルの脳波用増幅器４と、増幅後の電位をディジタルデータに変換するＡ／Ｄ変換器５と、予め構築した脳内の双極子から頭皮上への電位の伝達行列及び脳内に仮定した任意数の双極子の初期位置のデータを格納したファイルシステム１４と、電位の実測データと、ファイルシステム１４から読み込んだ伝達行列及び初期位置のデータとを使用し最小二乗法による解析処理を複数の副計算機８Ａ、８Ｂ・・・８Ｎによる初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する解析処理手段６を有するものである。【選択図】図１
2006325754.tif 000002
【特許請求の範囲】
【請求項１】
脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とするように各双極子の位置とベクトル成分を実時間で決定する脳波双極子解析装置であって、
被検体の頭皮上に配置する脳波計測用の複数の電極と、
予め構築した脳内の双極子から頭皮上への電位の伝達行列及び脳内に仮定した任意数の双極子の初期位置のデータを格納したファイルシステムと、
前記複数の電極による電位の実測データと、前記ファイルシステムから読み込んだ伝達行列及び初期位置のデータとを使用し最小二乗法による解析処理を複数の計算機による初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する解析処理手段と、
解析結果を出力する出力手段と、
を有することを特徴とする脳波双極子解析装置。
【請求項２】
脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とするように各双極子の位置とベクトル成分を実時間で決定する脳波双極子解析装置であって、
被検体の頭皮上に配置する脳波計測用の複数の電極と、
各電極を効率的に配置するための電極キャップと、
各電極の計測電位を増幅する多チャンネルの脳波用増幅器と、
増幅後の電位をディジタルデータに変換するＡ／Ｄ変換器と、
予め構築した脳内の双極子から頭皮上への電位の伝達行列及び脳内に仮定した任意数の双極子の初期位置のデータを格納したファイルシステムと、
前記Ａ／Ｄ変換器から出力される電位の実測データと、前記ファイルシステムから読み込んだ伝達行列及び初期位置のデータとを使用し最小二乗法による解析処理を複数の計算機による初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する解析処理手段と、
解析結果を可視的に出力する出力手段と、
を有することを特徴とする脳波双極子解析装置。
【請求項３】
脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とするように各双極子の位置とベクトル成分を実時間で決定する脳波双極子解析装置であって、
被検体の頭皮上に配置する脳波計測用の複数の電極と、
各電極を効率的に配置するための電極キャップと、
各電極の計測電位を増幅する多チャンネルの脳波用増幅器と、
増幅後の電位をディジタルデータに変換するＡ／Ｄ変換器と、
予め構築した脳内の双極子から頭皮上への電位の伝達行列及び脳内に仮定した任意数の双極子の初期位置のデータを格納したファイルシステムと、
前記Ａ／Ｄ変換器から出力される電位の実測データと、前記ファイルシステムから読み込んだ伝達行列及び初期位置のデータとを使用し最小二乗法による解析処理を複数の計算機による初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する解析処理手段と、
解析結果を可視的に出力する表示装置及びプリンタ及び解析結果を記憶する記憶手段と、
を有することを特徴とする脳波双極子解析装置。
【請求項４】
予め設定した被検体の頭蓋モデルと被検体の頭皮上の複数の電極の配置により脳内の双極子から頭皮上への電位の伝達行列を求める過程と、
脳内に仮定する双極子の初期位置のデータを設定する過程と、
被検体の頭皮上に配置する脳波計測用の複数の電極からの実測データを得る過程と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する過程と、
を有することを特徴とする脳波双極子解析方法。
【請求項５】
予め設定した被検体の頭部の各組織分画の形状及び導電率によって定まる頭蓋モデルと被検体の頭皮上の複数の電極の配置により脳内の双極子から頭皮上への電位の伝達行列を求める過程と、
脳内に仮定する双極子の初期位置のデータを設定する過程と、
被検体の頭皮上に配置する脳波計測用の複数の電極から出力されディジタル変換された実測データを得る過程と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、各並列処理結果の中から脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で選定する過程と、
選定結果を可視的に出力する過程と、
を有することを特徴とする脳波双極子解析方法。
【請求項６】
予め設定した被検体の頭部の各組織分画の形状及び導電率によって定まる頭蓋モデルと被検体の頭皮上の複数の電極の配置により脳内の双極子から頭皮上への電位の伝達行列を求める過程と、
脳内に仮定する双極子の初期位置のデータを設定する過程と、
被検体の頭皮上に配置する脳波計測用の複数の電極から出力されディジタル変換された実測データを得る過程と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、各並列処理結果の中から脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で選定する過程と、
選定結果の表示処理、印刷処理、記憶処理を行う過程と、
を有することを特徴とする脳波双極子解析方法。
【請求項７】
予め設定した被検体の頭蓋モデルと被検体の頭皮上の複数の電極の配置により脳内の双極子から頭皮上への電位の伝達行列を求める処理と、
脳内に仮定する双極子の初期位置のデータを設定する処理と、
被検体の頭皮上に配置する脳波計測用の複数の電極からの実測データを得る処理と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する処理と、
を行うことを特徴とする脳波双極子解析用プログラム。
【請求項８】
予め設定した被検体の頭部の各組織分画の形状及び導電率によって定まる頭蓋モデルと被検体の頭皮上の複数の電極の配置により脳内の双極子から頭皮上への電位の伝達行列を求める処理と、
脳内に仮定する双極子の初期位置のデータを設定する処理と、
被検体の頭皮上に配置する脳波計測用の複数の電極出力されディジタル変換された実測データを得る処理と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、各並列処理結果の中から脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で選定する過程と、
選定結果の可視的な出力処理と、
を行うことを特徴とする脳波双極子解析用プログラム。
【請求項９】
予め設定した被検体の頭部の各組織分画の形状及び導電率によって定まる頭蓋モデルと被検体の頭皮上の複数の電極の配置により脳内の双極子から頭皮上への電位の伝達行列を求める処理と、
脳内に仮定する双極子の初期位置のデータを設定する処理と、
被検体の頭皮上に配置する脳波計測用の複数の電極出力されディジタル変換された実測データを得る処理と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、各並列処理結果の中から脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で選定する処理と、
選定結果の表示、印刷及び記憶の各処理と、
を行うことを特徴とする脳波双極子解析用プログラム。
【請求項１０】
予め設定した被検体の頭蓋モデルと被検体の頭皮上の複数の電極の配置により脳内の双極子から頭皮上への電位の伝達行列を求める処理と、
脳内に仮定する双極子の初期位置のデータを設定する処理と、
被検体の頭皮上に配置する脳波計測用の複数の電極からの実測データを得る処理と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する処理と、
を行う脳波双極子解析用プログラムを記憶したことを特徴とするコンピュータ読み取り可能な記憶媒体。
【請求項１１】
予め設定した被検体の頭部の各組織分画の形状及び導電率によって定まる頭蓋モデルと被検体の頭皮上の複数の電極の配置により脳内の双極子から頭皮上への電位の伝達行列を求める処理と、
脳内に仮定する双極子の初期位置のデータを設定する処理と、
被検体の頭皮上に配置する脳波計測用の複数の電極出力されディジタル変換された実測データを得る処理と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、各並列処理結果の中から脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で選定する過程と、
選定結果の可視的な出力処理と、
を行う脳波双極子解析用プログラムを記憶したことを特徴とするコンピュータ読み取り可能な記憶媒体。
【請求項１２】
予め設定した被検体の頭部の各組織分画の形状及び導電率によって定まる頭蓋モデルと被検体の頭皮上の複数の電極の配置により脳内の双極子から頭皮上への電位の伝達行列を求める処理と、
脳内に仮定する双極子の初期位置のデータを設定する処理と、
被検体の頭皮上に配置する脳波計測用の複数の電極出力されディジタル変換された実測データを得る処理と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、各並列処理結果の中から脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で選定する処理と、
選定結果の表示、印刷及び記憶の各処理と、
を行う脳波双極子解析用プログラムを記憶したことを特徴とするコンピュータ読み取り可能な記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、脳波双極子解析装置、脳波双極子解析方法、脳波双極子解析用プログラム及び該プログラムを記憶した記憶媒体に関するものである。
【背景技術】
【０００２】
脳内の起電力を少数の電流双極子で近似し、各双極子の位置とベクトル成分を頭表皮上に配置した複数の電極で実測される電位分布から最小二乗法を使用して決定する手法は等価双極子法と呼ばれ、脳内活動の非侵襲的な推定方法として以前から利用されている。これは脳内に仮定した双極子が発生する電位分布と実測の電位分布との二乗誤差を最小とするように各双極子の位置とベクトル成分を決定する手法である。
【０００３】
一つの双極子は、６個の自由度（位置とベクトル成分に関してそれぞれ３個）を持つため、Ｎ個の双極子を使った場合には二乗誤差は６Ｎ変数の関数となるが、ベクトル成分が頭皮上の電位分布と線形関係にあることを利用すると二乗誤差を最小化する問題は双極子の位置（３Ｎ変数）のみに依存する評価関数を最小化する問題に還元される。
【０００４】
評価関数の最小化にはシンプレックス法に代表される繰返し法に頼らざるを得ない。つまり、適当な初期位置から出発し、評価関数が減少する方向に双極子を少しずつ移動させて行く手法である。
【０００５】
このため、双極子を移動させる度に電位計算を繰り返す必要がある。更に、評価関数は複数の極小点を持つ可能性があり、真の最小点を高い確率で見出すため、複数の初期位置から出発して最小化を繰返し、得られる複数の結果の中から評価関数を最小とするものを選ぶ、という手法が用いられている。
【０００６】
以上、二種類の要因により、電位分布計算の回数は膨大なものとなる。電位計算の高速化のためには頭蓋モデルとしては単純なものを採用することが好ましいが、電位計算の精度を向上させるには現実の頭蓋を良く近似表現する複雑なモデルが望まれる。実際には計算速度と計算精度との兼ね合いから実形状の３層頭蓋モデルが採用されることが多く、こうした３層頭蓋モデルをはじめそれ以外の多層の頭蓋モデルに特化した高速な電位計算手法が提案されている。
【０００７】
以上の事柄は以下の文献に詳述されている。
Tbshimitsu Musha & Yoshiwo Okamoto:Forward and Inverse Problems of EEG  Dipole Localization,Critical Reviews in Biomedical Engineering,vol.27,no.3-5,189-239,1999
【０００８】
脳内双極子の位置とベクトル成分を決定する従来技術として、例えば特許文献１に開示された等価双極子装置が提案されている。特許文献１の等価双極子装置は、頭蓋内の空洞部、頭蓋骨等の種々の物質による計算誤差を補正し、正確な脳内双極子の位置とベクトル成分を推定するように構成したものである。
【０００９】
しかし、この特許文献１の場合においても、電界分布計算の回数は膨大であり、更に補正演算も必要であるため、双極子解析の実時間解析を実現し得る高速処理を実現することは容易ではない。
【特許文献１】特公平７−７９８０４号公報
【発明の開示】
【発明が解決しようとする課題】
【００１０】
上述のように、脳内の双極子についての双極子解析には膨大な計算が必要であることから、従来においては実時間解析を実現し得る高速処理が可能な脳波双極子解析装置、方法が存在しない点である。
【課題を解決するための手段】
【００１１】
本発明は、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とするように各双極子の位置とベクトル成分を実時間で決定する脳波双極子解析装置であって、被検体の頭皮上に配置する脳波計測用の複数の電極と、予め構築した脳内の双極子から頭皮上への電位の伝達行列及び脳内に仮定した任意数の双極子の初期位置のデータを格納したファイルシステムと、前記複数の電極による電位の実測データと、前記ファイルシステムから読み込んだ伝達行列及び初期位置のデータとを使用し最小二乗法による解析処理を複数の計算機による初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する解析処理手段と、解析結果を出力する出力手段とを有することを最も特徴とする。
【発明の効果】
【００１２】
本発明によれば、以下の効果を奏する。
請求項１乃至３記載の発明によれば、頭部内への非侵襲な構成の基に、頭部における頭皮表面上に配置した複数の電極で測定される電位分布から脳内の起電力を近似する双極子の位置とベクトル成分を計算機による並列処理を実行することで実時間で高速決定（推定）し、決定結果を実時間で表示したり印刷したり、更には記憶保持することが可能な脳波双極子解析装置を提供することができる。
【００１３】
請求項４乃至６記載の発明によれば、本発明の脳波双極子解析装置を使用して脳内の起電力を近似する双極子の位置とベクトル成分を、並列的な双極子解析処理を実行することで実時間で高速決定（推定）し、決定結果を実時間で表示したり印刷したり、更には記憶保持することが可能な脳波双極子解析方法を提供することができる。
【００１４】
請求項７乃至９記載の発明によれば、伝達行列の決定、脳内に仮定する双極子の初期位置のデータの設定、脳波実測データの取得、これらに基づく高速な解析処理からなる一連の処理を実行し、脳内の起電力を近似する双極子の位置とベクトル成分を実時間で決定（推定）し、決定結果を実時間で表示したり印刷したり、更には記憶保持する処理を的確に実現できる脳波双極子解析用プログラムを提供することができる。
【００１５】
請求項１０乃至１２記載の発明によれば、本発明の脳波双極子解析装置を使用した高速処理可能な脳波双極子解析方法を実現することができるコンピュータシステムを容易に構築することができる記憶媒体を提供することができる。
【発明を実施するための最良の形態】
【００１６】
本発明は、脳内双極子の位置とベクトル成分とを脳波計測に対応して実時間で高速解析処理する脳波双極子解析装置を提供するという目的を、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とするように各双極子の位置とベクトル成分を実時間で決定する脳波双極子解析装置であって、被検体の頭皮上に配置する脳波計測用の複数の電極と、各電極を効率的に配置するための電極キャップと、各電極の計測電位を増幅する多チャンネルの脳波用増幅器と、増幅後の電位をディジタルデータに変換するＡ／Ｄ変換器と、予め構築した脳内の双極子から頭皮上への電位の伝達行列及び脳内に仮定した任意数の双極子の初期位置のデータを格納したファイルシステムと、前記Ａ／Ｄ変換器から出力される電位の実測データと、前記ファイルシステムから読み込んだ伝達行列及び初期位置のデータとを使用し最小二乗法による解析処理を複数の計算機による初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する解析処理手段と、解析結果を可視的に出力する出力手段とを有する構成に実現した。
【実施例】
【００１７】
以下に、本発明の実施例を詳細に説明する。
図１、図２は、本発明の実施例に係る実時間処理が可能な脳波双極子解析装置１の構成を示すものである。
本発明の実施例に係る脳波双極子解析装置１は、被験者の頭部３０における頭皮上に配置された脳波計測用の複数の電極２と、各電極２を効率的に配置するための電極キャップ３と、前記各電極２の電位を増幅する多チャンネルの脳波用増幅器４と、脳波用増幅器４の出力をディジタルデータに変換するＡ／Ｄ変換器５と、Ａ／Ｄ変換器５から出力される電位データを取り込み、実時間の双極子解析用のプログラム（詳細は後述する）に基く最小二乗法による双極子解析、そのための前処理、初期設定等を行う解析処理手段６を構成する主計算機７及び複数の副計算機８Ａ、８Ｂ・・・８Ｎ（Ｎは正の整数）と、図２に示すように解析処理手段６の処理結果を表示するための表示装置９と、処理結果を印刷するプリンタ１０と、処理結果を記憶する記憶手段であるハードディスク、ＣＤ−ＲＡＭ等の補助記憶部１１と、脳組織の導電率等のデータを指定するためのデータ入力を行うキーボード等の入力部１２とからなる入出力機器１３と、前記解析処理手段６の計算処理に必要な各種データを記憶するためのデータ格納ファイルを含むファイルシステム１４とを有している。
【００１８】
主計算機７は、図２に示すように、前記副計算機８Ａ、８Ｂ・・・８Ｎとの間で各種データの交信を行うための副計算機用インターフェース２１と、前記Ａ／Ｄ変換器５、表示装置９、プリンタ１０、補助記憶部１１、入力部１２からなる入出力機器１３及びファイルシステム１４用の入出力インターフェース２２と、詳細は後述する処理を実行可能な実時間の双極子解析用のプログラムを格納したプログラム格納部２３と、前記処理プログラムの実行、計算処理等を行うＣＰＵ２４と具備している。
【００１９】
この主計算機７は、実時間の双極子解析のための前処理、すなわち実時間の双極子解析を実施する前にその計算に必要な各種のデータを設定したり、初期設定、すなわち計算開始の直前にそれらのデータを副計算機８Ａ、８Ｂ・・・８Ｎに配信する処理を実行し、更に前記複数の電極２による実測データに基づき実時間で双極子解析処理を実行するようになっている。
【００２０】
次に、本実施例に係る脳波双極子解析装置１を使用し、双極子解析用のプログラムに基づいた実時間処理による脳波双極子解析方法について詳述する。
【００２１】
本実施例に係る脳波双極子解析方法は、図３に示すように、脳波双極子解析に先立ち双極子解析用のプログラムに基づいて脳内の双極子ｄ（図１に白丸を付して示す）から頭皮上への電位の伝達状態を記述する伝達行列を求める前処理（ステップＳ１０）、解析処理手段６による解析処理のための初期設定（ステップＳ２０）、解析処理手段６による解析処理（ステップＳ３０）、解析結果の出力（ステップＳ５０）の各過程を含んで実現される。以下これら各過程について項目別に説明する。
【００２２】
図４を参照して前処理（ステップＳ１０）について説明する。双極子計算に際しては脳内の双極子ｄから頭皮上の電位への伝達行列が必要である。伝達行列は、頭蓋モデルと頭皮上での電極配置で決まり例えばｎ行ｍ列の行列で記述される。また、前記頭蓋モデルは各組織分画の形状と導電率によって決まる。
【００２３】
（ａ）頭蓋モデル（頭蓋形状モデル）の設定（ステップＳ１１）
被検体（被験者）の頭部３０のＭＲＩ画像、あるいはＸ線ＣＴ画像が利用できる場合は、これらの断層像から被験者に個人適合する頭蓋モデルを構築する。このような断層像が利用できない場合は、標準の頭蓋モデルを被検体に合わせて変形する。つまり、被検体の眉間・外後頭隆起間距離、左右両耳介前点間距離、及び頭皮沿いに眉間から頭項を経て外後頭隆に至るまでの距離を計測し、これらが一致するように標準モデルを各軸方向に伸長圧縮して頭蓋モデルを構築する。
【００２４】
（ｂ）各組織の導電率の指定（ステップＳ１２）
３層頭蓋モデルであれば、脳組織、骨組織及び表皮の導電率を指定し、４層頭蓋モデルではそれらに加えて脳脊髄液の導電率を指定する。必要ならば通流実験等によって得られる導電率を使用するが、通常は標準値を採用する。
【００２５】
（ｃ）電極配置の設定（ステップＳ１３）
国際標準の１０・２０法に基づいて電極を配置する場合には、頭蓋モデルを基に主計算機７の内部で自動的に電極位置を算出する。独自の手法で電極を配置する場合には、個々の電極位置を計測してそのデータを主計算機７に入力する。
【００２６】
（ｄ）各データのデータファイルへの格納（ステップＳ１４）
上述した頭蓋モデル、導電率、電極配置の各データを前記ファイルシステム１４のデータ格納ファイルへ格納する。
【００２７】
（ｅ）伝達行列の計算（ステップＳ１５）
主計算機７は、ファイルシステム１４のデータ格納ファイルに格納した頭蓋モデル、導電率、電極配置の各データを使用して伝達行列を計算する。通常、伝達行列の計算は、各被検体に対して１度だけ実施すれば良い。ただし、電極配置を変更した場合や、開頭術等で頭蓋の形状や導電率分布が変化した場合には再度の計算が必要となる。
【００２８】
（ｆ）計算結果の格納（ステップＳ１６）
主計算機７は、計算し求めた伝達行列を前記ファイルシステム１４のデータ格納ファイルへ格納する。
【００２９】
次に、図５を参照して初期設定（ステップＳ２０）について説明する。
主計算機７は、ファイルシステム１４のデータ格納ファイルに予め設定され格納された双極子の初期位置のリストを読み込む（ステップＳ２１）この場合初期位置がＮ通りある場合には、主計算機７は高速処理を実現するためにＮ個の副計算機８Ａ、８Ｂ・・・８Ｎにそれぞれ異なる初期位置のデータを配信する（ステップＳ２２）。
次に、主計算機７はデータ格納ファイルから当該被検体の伝達行列のデータを読み込み（ステップＳ２３）、ステップＳ２２の処理で初期位置を配布した全ての副計算機８Ａ、８Ｂ・・・８Ｎに伝達行列のデータを配信する（ステップＳ２４）。その結果、主計算機７はオペレータによる計算開始指令の待機状態となる。
【００３０】
次に、図６を参照して双極子解析（ステップＳ３０）について説明する。
主計算機７は、オペレータによる計算開始指令があると（ステップＳ３１）以下の手順を繰り返し実行する。すなわち、主計算機７は、まずＡ／Ｄ変換器５に変換開始指令を送信する（ステップＳ３２）。これにより、Ａ／Ｄ変換器５は各電極２から出力される脳波の実測データをディジタルデータに変換して、主計算機７の入出力インターフェース２２へ送る。
【００３１】
主計算機７は、Ａ／Ｄ変換器５のディジタルデータを取り込み（ステップＳ３３）、次にオフセット、ドリフト除去の処理を行う（ステップＳ３４）。すなわち、脳波の測定にどのような電極２を使うにせよ分極電位の変動等に起因するオフセットやドリフトの混入は避けられない。例えば、適当な刺激に対する誘発脳波をオフラインで双極子解析する場合には、刺激前の適当な時間区間における電位の平均値を零とするように脳波のベースラインを決定できるが、実時間で双極子解析では脳電位を零とすべき明確な時間帯は存在しない。このため、適当な時間長Ｔを指定し、時刻ｔで測定された電位分布を双極子解析する場合には、時刻ｔ−Ｔから時刻ｔまでの平均値が零となるようにベースラインを決定するものである。
【００３２】
次に、主計算機７は、オフセット、ドリフト除去処理後の電位分布データをＮ個の副計算機８Ａ、８Ｂ・・・８Ｎに配信し（ステップＳ３５）、更に計算開始指令を出す（ステップＳ３６）。
【００３３】
各副計算機８Ａ、８Ｂ・・・８Ｎは、計算開始指令に基づき予め配信された異なる初期位置のデータ及び伝達行列のデータと、オフセット、ドリフト除去処理後の電位分布データを使用し、各々双極子解析処理を並列で実行する（ステップＳ３７）。具体的には、最小二乗法に基づき、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とするような双極子解析処理を実時間で実行する。この場合、極小値に収束するまでの繰返し回数は初期値によって異なり、最も多くの繰返しを要した初期値からの計算時間が全体の計算時間を決定するが、繰返しの上限を設定しておけば、計算完了までの最大の時間を保証することができる。
【００３４】
次に、各副計算機８Ａ、８Ｂ・・・８Ｎは、双極子解析の処理結果を各々主計算機７に送信する（ステップＳ３８）。主計算機７は、各副計算機８Ａ、８Ｂ・・・８Ｎからの処理結果を収集し、その中から最適な処理結果（電位分布の実測値と計算値との二乗誤差を最小とするものの中の最も小さいもの）を選択して、頭部３０内の各双極子の実際の位置とベクトル成分を実時間で決定（推定）する（ステップＳ３９）。
【００３５】
主計算機７は、更に必要に応じて最適な処理結果を可視的に出力するための処理、すなわち表示装置９による表示のための出力処理、プリンタ１０による印刷のための出力処理を行い、更には、前記補助記憶部１１への書き込み処理等を行う（ステップＳ４０）。
【００３６】
この後、オペレータによる停止命令があれば（ステップＳ４１）、一連の処理を終了し、停止命令が無ければステップＳ２以降の処理を繰り返す。
【００３７】
以上説明した本実施例の脳波双極子解析装置１を使用し、双極子解析用のプログラムに基づいた実時間処理による脳波双極子解析方法によれば、頭部３０内への非侵襲な構成の基に、この頭部３０における頭皮表面上に配置した複数の電極２で測定される電位分布から脳内の起電力を近似する双極子の位置とベクトル成分を副計算機８Ａ、８Ｂ・・・８Ｎによる並列処理を実行することで実時間で高速決定（推定）し、決定結果を実時間で表示したり印刷したり、更には記憶保持することが可能となる。
【００３８】
また、上述したように、予め設定した被検体の頭蓋モデルと被検体の頭皮上の複数の電極２の配置により脳内の双極子から頭皮上への電位の伝達行列を求める処理と、脳内に仮定する双極子の初期位置のデータを設定する処理と、被検体の頭皮上に配置する脳波計測用の複数の電極２からの実測データを得る処理と、伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する処理とを行う脳波双極子解析用プログラムにより、脳内の起電力を近似する双極子の位置とベクトル成分を実時間で決定（推定）し、決定結果を実時間で表示したり印刷したり、更には記憶保持する処理を的確に実現できる。
【００３９】
更に、前記脳波双極子解析用プログラムを、例えばＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ等の記憶媒体に記憶し、コンピュータ読み取り可能とすることで、この記憶媒体を用いて、本実施例の脳波双極子解析装置１を使用した高速処理可能な脳波双極子解析方法を実現することができるコンピュータシステムを容易に構築することができる。
【産業上の利用可能性】
【００４０】
本発明における上述した最適な処理結果は、BCI(Brain Computer Interface)データとして活用することも可能である。また、本発明において、主計算機は個の副計算機を兼ねることができ、その場合には、Ｎ個の初期位置に対してＮ個の計算機（主計算機とＮ−１個の副計算機）があれば良い。高速性を重要視しない場合には、一つの計算機（主計算機）が複数の副計算機の役割を演ずることで安価なシステムを構築することも可能である。
【図面の簡単な説明】
【００４１】
【図１】本発明の実施例に係る脳波双極子解析装置の概略構成を示すブロック図である。
【図２】本発明の実施例に係る脳波双極子解析装置の主計算機、入出力機器の詳細を示すブロック図である。
【図３】本発明の実施例に係る脳波双極子解析方法の全体の処理の流れを示すフローチャートである。
【図４】本発明の実施例に係る脳波双極子解析方法の前処理の内容を示すフローチャートである。
【図５】本発明の実施例に係る脳波双極子解析方法の初期設定の内容を示すフローチャートである。
【図６】本発明の実施例に係る脳波双極子解析方法の双極子解析処理の内容を示すフローチャートである。
【符号の説明】
【００４２】
ｄ    脳内の双極子
１    脳波双極子解析装置
２    電極
３    電極キャップ
４    脳波用増幅器
５    Ａ／Ｄ変換器
６    解析処理手段
７    主計算機
８Ａ  副計算機
８Ｂ  副計算機
８Ｎ  副計算機
９    表示装置
１０    プリンタ
１１    補助記憶部
１２    入力部
１３    入出力機器
１４    ファイルシステム
２１    副計算機用インターフェース
２２    入出力インターフェース
２３    プログラム格納部
２４    ＣＰＵ
３０    頭部
【図１】
2006325754.tif 000003
【図２】
2006325754.tif 000004
【図３】
2006325754.tif 000005
【図４】
2006325754.tif 000006
【図５】
2006325754.tif 000007
【図６】
2006325754.tif 000008