JP2006325754A - Electroencephalogram (eeg) dipole tracing apparatus, eeg dipole tracing method, program for eeg dipole tracing, and storage medium storing the program - Google Patents
Electroencephalogram (eeg) dipole tracing apparatus, eeg dipole tracing method, program for eeg dipole tracing, and storage medium storing the program Download PDFInfo
- Publication number
- JP2006325754A JP2006325754A JP2005151416A JP2005151416A JP2006325754A JP 2006325754 A JP2006325754 A JP 2006325754A JP 2005151416 A JP2005151416 A JP 2005151416A JP 2005151416 A JP2005151416 A JP 2005151416A JP 2006325754 A JP2006325754 A JP 2006325754A
- Authority
- JP
- Japan
- Prior art keywords
- dipole
- brain
- data
- scalp
- initial position
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 80
- 238000012545 processing Methods 0.000 claims abstract description 90
- 210000004556 brain Anatomy 0.000 claims abstract description 87
- 238000004458 analytical method Methods 0.000 claims abstract description 67
- 210000004761 scalp Anatomy 0.000 claims abstract description 51
- 239000011159 matrix material Substances 0.000 claims abstract description 44
- 238000005259 measurement Methods 0.000 claims abstract description 40
- 238000012546 transfer Methods 0.000 claims description 42
- 210000003625 skull Anatomy 0.000 claims description 29
- 210000003128 head Anatomy 0.000 claims description 14
- 210000001519 tissue Anatomy 0.000 claims description 8
- 230000000007 visual effect Effects 0.000 claims 2
- 230000005540 biological transmission Effects 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 description 26
- 238000013500 data storage Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 5
- 238000007781 pre-processing Methods 0.000 description 5
- 210000005013 brain tissue Anatomy 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 210000004709 eyebrow Anatomy 0.000 description 2
- 238000010223 real-time analysis Methods 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 1
- 230000007177 brain activity Effects 0.000 description 1
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000007428 craniotomy Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000004070 electrodeposition Methods 0.000 description 1
- 210000002615 epidermis Anatomy 0.000 description 1
- 230000000763 evoking effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 230000005477 standard model Effects 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
Description
本発明は、脳波双極子解析装置、脳波双極子解析方法、脳波双極子解析用プログラム及び該プログラムを記憶した記憶媒体に関するものである。 The present invention relates to an electroencephalogram dipole analyzer, an electroencephalogram dipole analysis method, an electroencephalogram dipole analysis program, and a storage medium storing the program.
脳内の起電力を少数の電流双極子で近似し、各双極子の位置とベクトル成分を頭表皮上に配置した複数の電極で実測される電位分布から最小二乗法を使用して決定する手法は等価双極子法と呼ばれ、脳内活動の非侵襲的な推定方法として以前から利用されている。これは脳内に仮定した双極子が発生する電位分布と実測の電位分布との二乗誤差を最小とするように各双極子の位置とベクトル成分を決定する手法である。 A method that approximates the electromotive force in the brain with a small number of current dipoles and determines the position and vector component of each dipole from the potential distribution measured by multiple electrodes arranged on the scalp using the least square method. Is called the equivalent dipole method and has been used for a long time as a non-invasive method of estimating brain activity. This is a method for determining the position and vector component of each dipole so as to minimize the square error between the potential distribution generated by the dipole assumed in the brain and the measured potential distribution.
一つの双極子は、6個の自由度(位置とベクトル成分に関してそれぞれ3個)を持つため、N個の双極子を使った場合には二乗誤差は6N変数の関数となるが、ベクトル成分が頭皮上の電位分布と線形関係にあることを利用すると二乗誤差を最小化する問題は双極子の位置(3N変数)のみに依存する評価関数を最小化する問題に還元される。 Since one dipole has 6 degrees of freedom (3 each for position and vector component), when N dipoles are used, the square error is a function of 6N variables, but the vector component is By utilizing the linear relationship with the potential distribution on the scalp, the problem of minimizing the square error is reduced to the problem of minimizing the evaluation function that depends only on the position of the dipole (3N variable).
評価関数の最小化にはシンプレックス法に代表される繰返し法に頼らざるを得ない。つまり、適当な初期位置から出発し、評価関数が減少する方向に双極子を少しずつ移動させて行く手法である。 In order to minimize the evaluation function, it is necessary to rely on an iterative method represented by the simplex method. That is, this is a method of starting from an appropriate initial position and moving the dipole little by little in a direction in which the evaluation function decreases.
このため、双極子を移動させる度に電位計算を繰り返す必要がある。更に、評価関数は複数の極小点を持つ可能性があり、真の最小点を高い確率で見出すため、複数の初期位置から出発して最小化を繰返し、得られる複数の結果の中から評価関数を最小とするものを選ぶ、という手法が用いられている。 For this reason, it is necessary to repeat the potential calculation every time the dipole is moved. Furthermore, the evaluation function may have a plurality of local minimum points. In order to find the true minimum point with a high probability, the evaluation function starts from a plurality of initial positions and is repeatedly minimized. The method of selecting the one that minimizes is used.
以上、二種類の要因により、電位分布計算の回数は膨大なものとなる。電位計算の高速化のためには頭蓋モデルとしては単純なものを採用することが好ましいが、電位計算の精度を向上させるには現実の頭蓋を良く近似表現する複雑なモデルが望まれる。実際には計算速度と計算精度との兼ね合いから実形状の3層頭蓋モデルが採用されることが多く、こうした3層頭蓋モデルをはじめそれ以外の多層の頭蓋モデルに特化した高速な電位計算手法が提案されている。 As described above, the number of potential distribution calculations becomes enormous due to two types of factors. In order to increase the speed of the potential calculation, it is preferable to use a simple skull model. However, in order to improve the accuracy of the potential calculation, a complex model that closely represents the actual skull is desired. Actually, a real three-layer skull model is often used because of the balance between calculation speed and calculation accuracy, and a high-speed potential calculation method specialized for such a three-layer skull model and other multi-layer skull models. Has been proposed.
以上の事柄は以下の文献に詳述されている。
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
The above matters are described in detail in the following documents.
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
脳内双極子の位置とベクトル成分を決定する従来技術として、例えば特許文献1に開示された等価双極子装置が提案されている。特許文献1の等価双極子装置は、頭蓋内の空洞部、頭蓋骨等の種々の物質による計算誤差を補正し、正確な脳内双極子の位置とベクトル成分を推定するように構成したものである。
As a conventional technique for determining the position and vector component of a brain dipole, for example, an equivalent dipole device disclosed in
しかし、この特許文献1の場合においても、電界分布計算の回数は膨大であり、更に補正演算も必要であるため、双極子解析の実時間解析を実現し得る高速処理を実現することは容易ではない。
上述のように、脳内の双極子についての双極子解析には膨大な計算が必要であることから、従来においては実時間解析を実現し得る高速処理が可能な脳波双極子解析装置、方法が存在しない点である。 As described above, since dipole analysis of a dipole in the brain requires enormous calculations, an electroencephalogram dipole analysis apparatus and method capable of high-speed processing capable of realizing real-time analysis have been conventionally used. It is a point that does not exist.
本発明は、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とするように各双極子の位置とベクトル成分を実時間で決定する脳波双極子解析装置であって、被検体の頭皮上に配置する脳波計測用の複数の電極と、予め構築した脳内の双極子から頭皮上への電位の伝達行列及び脳内に仮定した任意数の双極子の初期位置のデータを格納したファイルシステムと、前記複数の電極による電位の実測データと、前記ファイルシステムから読み込んだ伝達行列及び初期位置のデータとを使用し最小二乗法による解析処理を複数の計算機による初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する解析処理手段と、解析結果を出力する出力手段とを有することを最も特徴とする。 The present invention provides an electroencephalogram dipole analysis in which the position and vector components of each dipole are determined in real time so as to minimize the square error between the potential distribution generated by the assumed dipole in the brain and the measured potential distribution. A plurality of electrodes for electroencephalogram measurement arranged on the scalp of the subject, a transfer matrix of potential from the dipole in the brain to the scalp and an arbitrary number of dipoles assumed in the brain Using a file system storing initial position data, measured potential data of the plurality of electrodes, transfer matrix read from the file system and initial position data, and performing analysis processing by a least square method on a plurality of computers The optimal dipole position and vector components in the brain that minimize the square error between the potential distribution in which the assumed dipole is generated in the brain and the measured potential distribution are time And analysis means for determining, and most an outputting means for outputting an analysis result.
本発明によれば、以下の効果を奏する。
請求項1乃至3記載の発明によれば、頭部内への非侵襲な構成の基に、頭部における頭皮表面上に配置した複数の電極で測定される電位分布から脳内の起電力を近似する双極子の位置とベクトル成分を計算機による並列処理を実行することで実時間で高速決定(推定)し、決定結果を実時間で表示したり印刷したり、更には記憶保持することが可能な脳波双極子解析装置を提供することができる。
The present invention has the following effects.
According to the first to third aspects of the present invention, the electromotive force in the brain is obtained from the potential distribution measured by the plurality of electrodes arranged on the scalp surface in the head based on the noninvasive configuration in the head. It is possible to determine (estimate) the approximate dipole position and vector component in real time by executing parallel processing by a computer, display the result of the determination in real time, print it, and store it in memory. An electroencephalogram dipole analyzer can be provided.
請求項4乃至6記載の発明によれば、本発明の脳波双極子解析装置を使用して脳内の起電力を近似する双極子の位置とベクトル成分を、並列的な双極子解析処理を実行することで実時間で高速決定(推定)し、決定結果を実時間で表示したり印刷したり、更には記憶保持することが可能な脳波双極子解析方法を提供することができる。 According to the fourth to sixth aspects of the present invention, a parallel dipole analysis process is performed on the position and vector component of the dipole that approximates the electromotive force in the brain using the electroencephalogram dipole analyzer of the present invention. Thus, it is possible to provide an electroencephalogram dipole analysis method capable of performing high-speed determination (estimation) in real time, displaying the determination result in real time, printing it, and storing it in memory.
請求項7乃至9記載の発明によれば、伝達行列の決定、脳内に仮定する双極子の初期位置のデータの設定、脳波実測データの取得、これらに基づく高速な解析処理からなる一連の処理を実行し、脳内の起電力を近似する双極子の位置とベクトル成分を実時間で決定(推定)し、決定結果を実時間で表示したり印刷したり、更には記憶保持する処理を的確に実現できる脳波双極子解析用プログラムを提供することができる。 According to the seventh to ninth aspects of the present invention, a series of processes including determination of a transfer matrix, setting of initial position data of a dipole assumed in the brain, acquisition of brain wave measurement data, and high-speed analysis processing based on these To determine (estimate) the position and vector component of the dipole that approximates the electromotive force in the brain in real time, and to display and print the determination result in real time, and to further accurately store and hold the data. It is possible to provide an electroencephalogram dipole analysis program that can be realized in a simple manner.
請求項10乃至12記載の発明によれば、本発明の脳波双極子解析装置を使用した高速処理可能な脳波双極子解析方法を実現することができるコンピュータシステムを容易に構築することができる記憶媒体を提供することができる。
According to the invention described in
本発明は、脳内双極子の位置とベクトル成分とを脳波計測に対応して実時間で高速解析処理する脳波双極子解析装置を提供するという目的を、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とするように各双極子の位置とベクトル成分を実時間で決定する脳波双極子解析装置であって、被検体の頭皮上に配置する脳波計測用の複数の電極と、各電極を効率的に配置するための電極キャップと、各電極の計測電位を増幅する多チャンネルの脳波用増幅器と、増幅後の電位をディジタルデータに変換するA/D変換器と、予め構築した脳内の双極子から頭皮上への電位の伝達行列及び脳内に仮定した任意数の双極子の初期位置のデータを格納したファイルシステムと、前記A/D変換器から出力される電位の実測データと、前記ファイルシステムから読み込んだ伝達行列及び初期位置のデータとを使用し最小二乗法による解析処理を複数の計算機による初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する解析処理手段と、解析結果を可視的に出力する出力手段とを有する構成に実現した。 An object of the present invention is to provide an electroencephalogram dipole analyzer that analyzes a brain dipole position and a vector component in real time at high speed corresponding to electroencephalogram measurement. An electroencephalogram dipole analyzer that determines the position and vector components of each dipole in real time so as to minimize the square error between the potential distribution and the measured potential distribution. The electroencephalogram is placed on the scalp of the subject. A plurality of electrodes for measurement, an electrode cap for efficiently arranging each electrode, a multi-channel electroencephalogram amplifier that amplifies the measurement potential of each electrode, and an A / A that converts the amplified potential into digital data D converter, file system storing potential transfer matrix from dipole in brain to scalp and data of initial position of arbitrary number of dipoles assumed in brain, and A / D conversion Of the potential output from the Using the data, the transfer matrix read from the file system and the data of the initial position, analysis processing by the least square method is performed by parallel processing for each initial position by a plurality of computers, and the assumed dipole is generated in the brain. Analysis processing means for determining in real time the optimal dipole position and vector component in the brain that minimizes the square error between the potential distribution and the measured potential distribution; and output means for visually outputting the analysis result It realized in the composition which has.
以下に、本発明の実施例を詳細に説明する。
図1、図2は、本発明の実施例に係る実時間処理が可能な脳波双極子解析装置1の構成を示すものである。
本発明の実施例に係る脳波双極子解析装置1は、被験者の頭部30における頭皮上に配置された脳波計測用の複数の電極2と、各電極2を効率的に配置するための電極キャップ3と、前記各電極2の電位を増幅する多チャンネルの脳波用増幅器4と、脳波用増幅器4の出力をディジタルデータに変換するA/D変換器5と、A/D変換器5から出力される電位データを取り込み、実時間の双極子解析用のプログラム(詳細は後述する)に基く最小二乗法による双極子解析、そのための前処理、初期設定等を行う解析処理手段6を構成する主計算機7及び複数の副計算機8A、8B・・・8N(Nは正の整数)と、図2に示すように解析処理手段6の処理結果を表示するための表示装置9と、処理結果を印刷するプリンタ10と、処理結果を記憶する記憶手段であるハードディスク、CD−RAM等の補助記憶部11と、脳組織の導電率等のデータを指定するためのデータ入力を行うキーボード等の入力部12とからなる入出力機器13と、前記解析処理手段6の計算処理に必要な各種データを記憶するためのデータ格納ファイルを含むファイルシステム14とを有している。
Examples of the present invention will be described in detail below.
1 and 2 show the configuration of an
An
主計算機7は、図2に示すように、前記副計算機8A、8B・・・8Nとの間で各種データの交信を行うための副計算機用インターフェース21と、前記A/D変換器5、表示装置9、プリンタ10、補助記憶部11、入力部12からなる入出力機器13及びファイルシステム14用の入出力インターフェース22と、詳細は後述する処理を実行可能な実時間の双極子解析用のプログラムを格納したプログラム格納部23と、前記処理プログラムの実行、計算処理等を行うCPU24と具備している。
As shown in FIG. 2, the main computer 7 includes a
この主計算機7は、実時間の双極子解析のための前処理、すなわち実時間の双極子解析を実施する前にその計算に必要な各種のデータを設定したり、初期設定、すなわち計算開始の直前にそれらのデータを副計算機8A、8B・・・8Nに配信する処理を実行し、更に前記複数の電極2による実測データに基づき実時間で双極子解析処理を実行するようになっている。
This main computer 7 sets various data necessary for the pre-processing for the real-time dipole analysis, i.e., the real-time dipole analysis, or sets the initial settings, i.e., the calculation start. Immediately before, a process of distributing the data to the
次に、本実施例に係る脳波双極子解析装置1を使用し、双極子解析用のプログラムに基づいた実時間処理による脳波双極子解析方法について詳述する。
Next, an electroencephalogram dipole analysis method using real-time processing based on a program for dipole analysis using the
本実施例に係る脳波双極子解析方法は、図3に示すように、脳波双極子解析に先立ち双極子解析用のプログラムに基づいて脳内の双極子d(図1に白丸を付して示す)から頭皮上への電位の伝達状態を記述する伝達行列を求める前処理(ステップS10)、解析処理手段6による解析処理のための初期設定(ステップS20)、解析処理手段6による解析処理(ステップS30)、解析結果の出力(ステップS50)の各過程を含んで実現される。以下これら各過程について項目別に説明する。 In the electroencephalogram dipole analysis method according to the present embodiment, as shown in FIG. 3, the dipole d in the brain (shown with a white circle in FIG. 1) based on the dipole analysis program prior to the electroencephalogram dipole analysis. ) For preprocessing for obtaining a transfer matrix describing the transmission state of potential from the scalp to the scalp (step S10), initial setting for analysis processing by the analysis processing means 6 (step S20), analysis processing by the analysis processing means 6 (step S30) and an analysis result output process (step S50). Each of these processes will be described below by item.
図4を参照して前処理(ステップS10)について説明する。双極子計算に際しては脳内の双極子dから頭皮上の電位への伝達行列が必要である。伝達行列は、頭蓋モデルと頭皮上での電極配置で決まり例えばn行m列の行列で記述される。また、前記頭蓋モデルは各組織分画の形状と導電率によって決まる。 The preprocessing (step S10) will be described with reference to FIG. In the dipole calculation, a transfer matrix from the dipole d in the brain to the potential on the scalp is necessary. The transfer matrix is determined by the skull model and electrode arrangement on the scalp, and is described by, for example, a matrix of n rows and m columns. The skull model is determined by the shape and conductivity of each tissue fraction.
(a)頭蓋モデル(頭蓋形状モデル)の設定(ステップS11)
被検体(被験者)の頭部30のMRI画像、あるいはX線CT画像が利用できる場合は、これらの断層像から被験者に個人適合する頭蓋モデルを構築する。このような断層像が利用できない場合は、標準の頭蓋モデルを被検体に合わせて変形する。つまり、被検体の眉間・外後頭隆起間距離、左右両耳介前点間距離、及び頭皮沿いに眉間から頭項を経て外後頭隆に至るまでの距離を計測し、これらが一致するように標準モデルを各軸方向に伸長圧縮して頭蓋モデルを構築する。
(A) Setting of skull model (cranium shape model) (step S11)
When an MRI image or an X-ray CT image of the
(b)各組織の導電率の指定(ステップS12)
3層頭蓋モデルであれば、脳組織、骨組織及び表皮の導電率を指定し、4層頭蓋モデルではそれらに加えて脳脊髄液の導電率を指定する。必要ならば通流実験等によって得られる導電率を使用するが、通常は標準値を採用する。
(B) Designation of conductivity of each tissue (step S12)
In the case of a three-layer skull model, the conductivity of brain tissue, bone tissue, and epidermis is designated, and in the case of the four-layer skull model, the conductivity of cerebrospinal fluid is designated in addition to them. If necessary, the conductivity obtained by a flow experiment is used, but a standard value is usually adopted.
(c)電極配置の設定(ステップS13)
国際標準の10・20法に基づいて電極を配置する場合には、頭蓋モデルを基に主計算機7の内部で自動的に電極位置を算出する。独自の手法で電極を配置する場合には、個々の電極位置を計測してそのデータを主計算機7に入力する。
(C) Setting of electrode arrangement (step S13)
When the electrodes are arranged based on the international standard 10/20 method, the electrode positions are automatically calculated inside the main computer 7 based on the skull model. When the electrodes are arranged by a unique method, the positions of the individual electrodes are measured and the data is input to the main computer 7.
(d)各データのデータファイルへの格納(ステップS14)
上述した頭蓋モデル、導電率、電極配置の各データを前記ファイルシステム14のデータ格納ファイルへ格納する。
(D) Storage of each data in a data file (step S14)
The above-described skull model, conductivity, and electrode arrangement data are stored in the data storage file of the
(e)伝達行列の計算(ステップS15)
主計算機7は、ファイルシステム14のデータ格納ファイルに格納した頭蓋モデル、導電率、電極配置の各データを使用して伝達行列を計算する。通常、伝達行列の計算は、各被検体に対して1度だけ実施すれば良い。ただし、電極配置を変更した場合や、開頭術等で頭蓋の形状や導電率分布が変化した場合には再度の計算が必要となる。
(E) Transfer matrix calculation (step S15)
The main computer 7 calculates a transfer matrix using the skull model, conductivity, and electrode arrangement data stored in the data storage file of the
(f)計算結果の格納(ステップS16)
主計算機7は、計算し求めた伝達行列を前記ファイルシステム14のデータ格納ファイルへ格納する。
(F) Storing calculation results (step S16)
The main computer 7 stores the calculated transfer matrix in the data storage file of the
次に、図5を参照して初期設定(ステップS20)について説明する。
主計算機7は、ファイルシステム14のデータ格納ファイルに予め設定され格納された双極子の初期位置のリストを読み込む(ステップS21)この場合初期位置がN通りある場合には、主計算機7は高速処理を実現するためにN個の副計算機8A、8B・・・8Nにそれぞれ異なる初期位置のデータを配信する(ステップS22)。
次に、主計算機7はデータ格納ファイルから当該被検体の伝達行列のデータを読み込み(ステップS23)、ステップS22の処理で初期位置を配布した全ての副計算機8A、8B・・・8Nに伝達行列のデータを配信する(ステップS24)。その結果、主計算機7はオペレータによる計算開始指令の待機状態となる。
Next, the initial setting (step S20) will be described with reference to FIG.
The main computer 7 reads a list of initial positions of dipoles preset and stored in the data storage file of the file system 14 (step S21). In this case, if there are N initial positions, the main computer 7 performs high-speed processing. In order to realize the above, data at different initial positions are distributed to the N sub-computers 8A, 8B... 8N (step S22).
Next, the main computer 7 reads the transfer matrix data of the subject from the data storage file (step S23), and transfers the transfer matrix to all the
次に、図6を参照して双極子解析(ステップS30)について説明する。
主計算機7は、オペレータによる計算開始指令があると(ステップS31)以下の手順を繰り返し実行する。すなわち、主計算機7は、まずA/D変換器5に変換開始指令を送信する(ステップS32)。これにより、A/D変換器5は各電極2から出力される脳波の実測データをディジタルデータに変換して、主計算機7の入出力インターフェース22へ送る。
Next, dipole analysis (step S30) will be described with reference to FIG.
When there is a calculation start command from the operator (step S31), the main computer 7 repeatedly executes the following procedure. That is, the main computer 7 first transmits a conversion start command to the A / D converter 5 (step S32). As a result, the A /
主計算機7は、A/D変換器5のディジタルデータを取り込み(ステップS33)、次にオフセット、ドリフト除去の処理を行う(ステップS34)。すなわち、脳波の測定にどのような電極2を使うにせよ分極電位の変動等に起因するオフセットやドリフトの混入は避けられない。例えば、適当な刺激に対する誘発脳波をオフラインで双極子解析する場合には、刺激前の適当な時間区間における電位の平均値を零とするように脳波のベースラインを決定できるが、実時間で双極子解析では脳電位を零とすべき明確な時間帯は存在しない。このため、適当な時間長Tを指定し、時刻tで測定された電位分布を双極子解析する場合には、時刻t−Tから時刻tまでの平均値が零となるようにベースラインを決定するものである。 The main computer 7 takes in the digital data of the A / D converter 5 (step S33), and then performs offset and drift removal processing (step S34). In other words, no matter what electrode 2 is used for the measurement of the electroencephalogram, it is inevitable that offsets and drifts are caused by fluctuations in the polarization potential. For example, when the dipole analysis of the evoked electroencephalogram for an appropriate stimulus is performed offline, the baseline of the electroencephalogram can be determined so that the average value of the electric potential in the appropriate time interval before the stimulation is zero, In the child analysis, there is no clear time zone in which the brain potential should be zero. For this reason, when an appropriate time length T is designated and the potential distribution measured at time t is analyzed for dipoles, the baseline is determined so that the average value from time t-T to time t becomes zero. To do.
次に、主計算機7は、オフセット、ドリフト除去処理後の電位分布データをN個の副計算機8A、8B・・・8Nに配信し(ステップS35)、更に計算開始指令を出す(ステップS36)。 Next, the main computer 7 distributes the potential distribution data after the offset and drift removal processing to the N sub-computers 8A, 8B... 8N (step S35), and further issues a calculation start command (step S36).
各副計算機8A、8B・・・8Nは、計算開始指令に基づき予め配信された異なる初期位置のデータ及び伝達行列のデータと、オフセット、ドリフト除去処理後の電位分布データを使用し、各々双極子解析処理を並列で実行する(ステップS37)。具体的には、最小二乗法に基づき、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とするような双極子解析処理を実時間で実行する。この場合、極小値に収束するまでの繰返し回数は初期値によって異なり、最も多くの繰返しを要した初期値からの計算時間が全体の計算時間を決定するが、繰返しの上限を設定しておけば、計算完了までの最大の時間を保証することができる。
Each of the
次に、各副計算機8A、8B・・・8Nは、双極子解析の処理結果を各々主計算機7に送信する(ステップS38)。主計算機7は、各副計算機8A、8B・・・8Nからの処理結果を収集し、その中から最適な処理結果(電位分布の実測値と計算値との二乗誤差を最小とするものの中の最も小さいもの)を選択して、頭部30内の各双極子の実際の位置とベクトル成分を実時間で決定(推定)する(ステップS39)。
Next, each of the
主計算機7は、更に必要に応じて最適な処理結果を可視的に出力するための処理、すなわち表示装置9による表示のための出力処理、プリンタ10による印刷のための出力処理を行い、更には、前記補助記憶部11への書き込み処理等を行う(ステップS40)。
The main computer 7 further performs processing for visually outputting the optimum processing result as necessary, that is, output processing for display by the display device 9, output processing for printing by the
この後、オペレータによる停止命令があれば(ステップS41)、一連の処理を終了し、停止命令が無ければステップS2以降の処理を繰り返す。 Thereafter, if there is a stop command from the operator (step S41), the series of processing is terminated, and if there is no stop command, the processing from step S2 is repeated.
以上説明した本実施例の脳波双極子解析装置1を使用し、双極子解析用のプログラムに基づいた実時間処理による脳波双極子解析方法によれば、頭部30内への非侵襲な構成の基に、この頭部30における頭皮表面上に配置した複数の電極2で測定される電位分布から脳内の起電力を近似する双極子の位置とベクトル成分を副計算機8A、8B・・・8Nによる並列処理を実行することで実時間で高速決定(推定)し、決定結果を実時間で表示したり印刷したり、更には記憶保持することが可能となる。
According to the electroencephalogram dipole analysis method using real time processing based on the dipole analysis program using the electroencephalogram
また、上述したように、予め設定した被検体の頭蓋モデルと被検体の頭皮上の複数の電極2の配置により脳内の双極子から頭皮上への電位の伝達行列を求める処理と、脳内に仮定する双極子の初期位置のデータを設定する処理と、被検体の頭皮上に配置する脳波計測用の複数の電極2からの実測データを得る処理と、伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する処理とを行う脳波双極子解析用プログラムにより、脳内の起電力を近似する双極子の位置とベクトル成分を実時間で決定(推定)し、決定結果を実時間で表示したり印刷したり、更には記憶保持する処理を的確に実現できる。 Further, as described above, a process for obtaining a transfer matrix of a potential from a dipole in the brain to the scalp based on a predetermined skull model of the subject and the arrangement of the plurality of electrodes 2 on the scalp of the subject; Processing for setting the initial position data of the dipole assumed in FIG. 2, processing for obtaining measured data from a plurality of electrodes 2 for electroencephalogram measurement placed on the scalp of the subject, transfer matrix, initial position data, actual measurement Based on the data, the analysis process by the least square method is performed by parallel processing for each initial position, and the optimal brain that minimizes the square error between the potential distribution in which the assumed dipole is generated in the brain and the measured potential distribution Determines (estimates) the position and vector component of the dipole that approximates the electromotive force in the brain using a brain wave dipole analysis program that performs real-time processing to determine the position and vector component of the dipole. , Display the decision result in real time To print Ri, more can be accurately realized a process for storing and holding.
更に、前記脳波双極子解析用プログラムを、例えばCD−ROM、DVD−ROM等の記憶媒体に記憶し、コンピュータ読み取り可能とすることで、この記憶媒体を用いて、本実施例の脳波双極子解析装置1を使用した高速処理可能な脳波双極子解析方法を実現することができるコンピュータシステムを容易に構築することができる。
Furthermore, the electroencephalogram dipole analysis program is stored in a storage medium such as a CD-ROM or a DVD-ROM, and can be read by a computer, so that the electroencephalogram dipole analysis of this embodiment can be performed using this storage medium. A computer system capable of realizing a brain wave dipole analysis method capable of high-speed processing using the
本発明における上述した最適な処理結果は、BCI(Brain Computer Interface)データとして活用することも可能である。また、本発明において、主計算機は個の副計算機を兼ねることができ、その場合には、N個の初期位置に対してN個の計算機(主計算機とN−1個の副計算機)があれば良い。高速性を重要視しない場合には、一つの計算機(主計算機)が複数の副計算機の役割を演ずることで安価なシステムを構築することも可能である。 The above-described optimum processing result in the present invention can be utilized as BCI (Brain Computer Interface) data. In the present invention, the main computer can also serve as sub-computers. In this case, there are N computers (the main computer and N-1 sub-computers) for N initial positions. It ’s fine. When high speed is not important, it is possible to construct an inexpensive system by one computer (main computer) playing the role of a plurality of sub computers.
d 脳内の双極子
1 脳波双極子解析装置
2 電極
3 電極キャップ
4 脳波用増幅器
5 A/D変換器
6 解析処理手段
7 主計算機
8A 副計算機
8B 副計算機
8N 副計算機
9 表示装置
10 プリンタ
11 補助記憶部
12 入力部
13 入出力機器
14 ファイルシステム
21 副計算機用インターフェース
22 入出力インターフェース
23 プログラム格納部
24 CPU
30 頭部
d Dipole in the
30 heads
Claims (12)
被検体の頭皮上に配置する脳波計測用の複数の電極と、
予め構築した脳内の双極子から頭皮上への電位の伝達行列及び脳内に仮定した任意数の双極子の初期位置のデータを格納したファイルシステムと、
前記複数の電極による電位の実測データと、前記ファイルシステムから読み込んだ伝達行列及び初期位置のデータとを使用し最小二乗法による解析処理を複数の計算機による初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する解析処理手段と、
解析結果を出力する出力手段と、
を有することを特徴とする脳波双極子解析装置。 An electroencephalogram dipole analyzer that determines the position and vector components of each dipole in real time so as to minimize the square error between the potential distribution of the assumed dipole in the brain and the measured potential distribution. ,
A plurality of electrodes for electroencephalogram measurement placed on the scalp of the subject;
A file system that stores a potential transfer matrix from a dipole in the brain to the scalp and data of an initial position of an arbitrary number of dipoles assumed in the brain,
Using the measured data of the potentials by the plurality of electrodes, the transfer matrix and the initial position data read from the file system, analysis processing by the least square method is performed by parallel processing for each initial position by a plurality of computers, An analysis processing means for determining in real time the optimal position and vector component of the dipole in the brain that minimizes the square error between the potential distribution generated by the dipole assumed in ii and the measured potential distribution;
An output means for outputting the analysis result;
An electroencephalogram dipole analyzing apparatus characterized by comprising:
被検体の頭皮上に配置する脳波計測用の複数の電極と、
各電極を効率的に配置するための電極キャップと、
各電極の計測電位を増幅する多チャンネルの脳波用増幅器と、
増幅後の電位をディジタルデータに変換するA/D変換器と、
予め構築した脳内の双極子から頭皮上への電位の伝達行列及び脳内に仮定した任意数の双極子の初期位置のデータを格納したファイルシステムと、
前記A/D変換器から出力される電位の実測データと、前記ファイルシステムから読み込んだ伝達行列及び初期位置のデータとを使用し最小二乗法による解析処理を複数の計算機による初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する解析処理手段と、
解析結果を可視的に出力する出力手段と、
を有することを特徴とする脳波双極子解析装置。 An electroencephalogram dipole analyzer that determines the position and vector components of each dipole in real time so as to minimize the square error between the potential distribution of the assumed dipole in the brain and the measured potential distribution. ,
A plurality of electrodes for electroencephalogram measurement placed on the scalp of the subject;
An electrode cap for efficiently arranging each electrode;
A multichannel EEG amplifier that amplifies the measurement potential of each electrode;
An A / D converter that converts the amplified potential into digital data;
A file system that stores a potential transfer matrix from a dipole in the brain to the scalp and data of an initial position of an arbitrary number of dipoles assumed in the brain,
Parallel processing for each initial position by a plurality of computers is performed using the least square method using measured data of the potential output from the A / D converter and the transfer matrix and initial position data read from the file system. Analytical processing means that determines the optimal dipole position and vector component in the brain that minimizes the square error between the potential distribution of the assumed dipole generated in the brain and the measured potential distribution in real time When,
An output means for visually outputting the analysis result;
An electroencephalogram dipole analyzing apparatus characterized by comprising:
被検体の頭皮上に配置する脳波計測用の複数の電極と、
各電極を効率的に配置するための電極キャップと、
各電極の計測電位を増幅する多チャンネルの脳波用増幅器と、
増幅後の電位をディジタルデータに変換するA/D変換器と、
予め構築した脳内の双極子から頭皮上への電位の伝達行列及び脳内に仮定した任意数の双極子の初期位置のデータを格納したファイルシステムと、
前記A/D変換器から出力される電位の実測データと、前記ファイルシステムから読み込んだ伝達行列及び初期位置のデータとを使用し最小二乗法による解析処理を複数の計算機による初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する解析処理手段と、
解析結果を可視的に出力する表示装置及びプリンタ及び解析結果を記憶する記憶手段と、
を有することを特徴とする脳波双極子解析装置。 An electroencephalogram dipole analyzer that determines the position and vector components of each dipole in real time so as to minimize the square error between the potential distribution of the assumed dipole in the brain and the measured potential distribution. ,
A plurality of electrodes for electroencephalogram measurement placed on the scalp of the subject;
An electrode cap for efficiently arranging each electrode;
A multichannel EEG amplifier that amplifies the measurement potential of each electrode;
An A / D converter that converts the amplified potential into digital data;
A file system that stores a potential transfer matrix from a dipole in the brain to the scalp and data of an initial position of an arbitrary number of dipoles assumed in the brain,
Parallel processing for each initial position by a plurality of computers is performed using the least square method using measured data of the potential output from the A / D converter and the transfer matrix and initial position data read from the file system. Analytical processing means that determines the optimal dipole position and vector component in the brain that minimizes the square error between the potential distribution of the assumed dipole generated in the brain and the measured potential distribution in real time When,
A display device and printer for visually outputting the analysis result, and storage means for storing the analysis result;
An electroencephalogram dipole analyzing apparatus characterized by comprising:
脳内に仮定する双極子の初期位置のデータを設定する過程と、
被検体の頭皮上に配置する脳波計測用の複数の電極からの実測データを得る過程と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する過程と、
を有することを特徴とする脳波双極子解析方法。 A process for obtaining a transfer matrix of potential from a dipole in the brain to the scalp by arranging a plurality of electrodes on the scalp of the subject and the subject's skull,
The process of setting the initial dipole data assumed in the brain;
The process of obtaining measured data from multiple electrodes for electroencephalogram measurement placed on the subject's scalp,
Based on the transfer matrix, the initial position data, and the measured data, the analysis process by the least square method is performed by parallel processing for each initial position, and the potential distribution in which the assumed dipole is generated in the brain and the measured potential distribution The process of determining in real time the optimal dipole position and vector component in the brain that minimizes the square error;
A method for analyzing an electroencephalogram dipole, comprising:
脳内に仮定する双極子の初期位置のデータを設定する過程と、
被検体の頭皮上に配置する脳波計測用の複数の電極から出力されディジタル変換された実測データを得る過程と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、各並列処理結果の中から脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で選定する過程と、
選定結果を可視的に出力する過程と、
を有することを特徴とする脳波双極子解析方法。 The transfer matrix of potential from the dipole in the brain to the scalp is determined by the skull model determined by the shape and conductivity of each tissue fraction of the subject's head set in advance and the arrangement of multiple electrodes on the subject's scalp. The process of seeking
The process of setting the initial dipole data assumed in the brain;
A process of obtaining actual measurement data output from a plurality of electrodes for electroencephalogram measurement placed on the scalp of the subject and converted into digital data,
Based on the transfer matrix, the initial position data, and the actual measurement data, analysis processing by the least square method is performed by parallel processing for each initial position, and a potential distribution in which a dipole assumed in the brain is generated from each parallel processing result and The process of selecting the optimal dipole position and vector component in the brain that minimizes the square error with the measured potential distribution in real time,
The process of visually outputting the selection result,
A method for analyzing an electroencephalogram dipole, comprising:
脳内に仮定する双極子の初期位置のデータを設定する過程と、
被検体の頭皮上に配置する脳波計測用の複数の電極から出力されディジタル変換された実測データを得る過程と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、各並列処理結果の中から脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で選定する過程と、
選定結果の表示処理、印刷処理、記憶処理を行う過程と、
を有することを特徴とする脳波双極子解析方法。 The transfer matrix of potential from the dipole in the brain to the scalp is determined by the skull model determined by the shape and conductivity of each tissue fraction of the subject's head set in advance and the arrangement of multiple electrodes on the subject's scalp. The process of seeking
The process of setting the initial dipole data assumed in the brain;
A process of obtaining actual measurement data output from a plurality of electrodes for electroencephalogram measurement placed on the scalp of the subject and converted into digital data,
Based on the transfer matrix, the initial position data, and the actual measurement data, analysis processing by the least square method is performed by parallel processing for each initial position, and a potential distribution in which a dipole assumed in the brain is generated from each parallel processing result and The process of selecting the optimal dipole position and vector component in the brain that minimizes the square error with the measured potential distribution in real time,
The process of displaying the selection result, printing, and storing,
A method for analyzing an electroencephalogram dipole, comprising:
脳内に仮定する双極子の初期位置のデータを設定する処理と、
被検体の頭皮上に配置する脳波計測用の複数の電極からの実測データを得る処理と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する処理と、
を行うことを特徴とする脳波双極子解析用プログラム。 A process for obtaining a transfer matrix of potential from a dipole in the brain to the scalp by arranging a plurality of electrodes on the scalp of the subject and the skull model of the subject,
Processing to set the initial position data of the dipole assumed in the brain;
Processing to obtain actual measurement data from a plurality of electrodes for electroencephalogram measurement placed on the scalp of the subject;
Based on the transfer matrix, the initial position data, and the measured data, the analysis process by the least square method is performed by parallel processing for each initial position, and the potential distribution in which the assumed dipole is generated in the brain and the measured potential distribution A process for determining in real time the optimal position and vector component of the dipole in the brain that minimizes the square error;
A program for analyzing EEG dipoles.
脳内に仮定する双極子の初期位置のデータを設定する処理と、
被検体の頭皮上に配置する脳波計測用の複数の電極出力されディジタル変換された実測データを得る処理と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、各並列処理結果の中から脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で選定する過程と、
選定結果の可視的な出力処理と、
を行うことを特徴とする脳波双極子解析用プログラム。 The transfer matrix of potential from the dipole in the brain to the scalp is determined by the skull model determined by the shape and conductivity of each tissue fraction of the subject's head set in advance and the arrangement of multiple electrodes on the subject's scalp. The required processing,
Processing to set the initial position data of the dipole assumed in the brain;
A process of obtaining a plurality of electrodes for electroencephalogram measurement placed on the scalp of the subject and obtaining actual measurement data converted into digital data,
Based on the transfer matrix, the initial position data, and the actual measurement data, analysis processing by the least square method is performed by parallel processing for each initial position, and a potential distribution in which a dipole assumed in the brain is generated from each parallel processing result and The process of selecting the optimal dipole position and vector component in the brain that minimizes the square error with the measured potential distribution in real time,
Visual output processing of selection results,
A program for analyzing EEG dipoles.
脳内に仮定する双極子の初期位置のデータを設定する処理と、
被検体の頭皮上に配置する脳波計測用の複数の電極出力されディジタル変換された実測データを得る処理と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、各並列処理結果の中から脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で選定する処理と、
選定結果の表示、印刷及び記憶の各処理と、
を行うことを特徴とする脳波双極子解析用プログラム。 The transfer matrix of potential from the dipole in the brain to the scalp is determined by the skull model determined by the shape and conductivity of each tissue fraction of the subject's head set in advance and the arrangement of multiple electrodes on the subject's scalp. The required processing,
Processing to set the initial position data of the dipole assumed in the brain;
A process of obtaining a plurality of electrodes for electroencephalogram measurement placed on the scalp of the subject and obtaining actual measurement data converted into digital data,
Based on the transfer matrix, the initial position data, and the actual measurement data, analysis processing by the least square method is performed by parallel processing for each initial position, and a potential distribution in which a dipole assumed in the brain is generated from each parallel processing result and In real time, select the optimal dipole position and vector component in the brain that minimizes the square error with the measured potential distribution,
Each process of displaying, printing and storing selection results;
A program for analyzing EEG dipoles.
脳内に仮定する双極子の初期位置のデータを設定する処理と、
被検体の頭皮上に配置する脳波計測用の複数の電極からの実測データを得る処理と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で決定する処理と、
を行う脳波双極子解析用プログラムを記憶したことを特徴とするコンピュータ読み取り可能な記憶媒体。 A process for obtaining a transfer matrix of potential from a dipole in the brain to the scalp by arranging a plurality of electrodes on the scalp of the subject and the skull model of the subject,
Processing to set the initial position data of the dipole assumed in the brain;
Processing to obtain actual measurement data from a plurality of electrodes for electroencephalogram measurement placed on the scalp of the subject;
Based on the transfer matrix, the initial position data, and the measured data, the analysis process by the least square method is performed by parallel processing for each initial position, and the potential distribution in which the assumed dipole is generated in the brain and the measured potential distribution A process for determining in real time the optimal position and vector component of the dipole in the brain that minimizes the square error;
A computer-readable storage medium characterized by storing a program for analyzing an electroencephalogram dipole that performs.
脳内に仮定する双極子の初期位置のデータを設定する処理と、
被検体の頭皮上に配置する脳波計測用の複数の電極出力されディジタル変換された実測データを得る処理と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、各並列処理結果の中から脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で選定する過程と、
選定結果の可視的な出力処理と、
を行う脳波双極子解析用プログラムを記憶したことを特徴とするコンピュータ読み取り可能な記憶媒体。 The transfer matrix of potential from the dipole in the brain to the scalp is determined by the skull model determined by the shape and conductivity of each tissue fraction of the subject's head set in advance and the arrangement of multiple electrodes on the subject's scalp. The required processing,
Processing to set the initial position data of the dipole assumed in the brain;
A process of obtaining a plurality of electrodes for electroencephalogram measurement placed on the scalp of the subject and obtaining actual measurement data converted into digital data,
Based on the transfer matrix, the initial position data, and the actual measurement data, analysis processing by the least square method is performed by parallel processing for each initial position, and a potential distribution in which a dipole assumed in the brain is generated from each parallel processing result and The process of selecting the optimal dipole position and vector component in the brain that minimizes the square error with the measured potential distribution in real time,
Visual output processing of selection results,
A computer-readable storage medium characterized by storing a program for analyzing an electroencephalogram dipole that performs.
脳内に仮定する双極子の初期位置のデータを設定する処理と、
被検体の頭皮上に配置する脳波計測用の複数の電極出力されディジタル変換された実測データを得る処理と、
前記伝達行列、初期位置のデータ、実測データを基に最小二乗法による解析処理を初期位置毎の並列処理により行い、各並列処理結果の中から脳内に仮定した双極子が発生する電位分布と、実測の電位分布との二乗誤差を最小とする最適な脳内の双極子の位置とベクトル成分を実時間で選定する処理と、
選定結果の表示、印刷及び記憶の各処理と、
を行う脳波双極子解析用プログラムを記憶したことを特徴とするコンピュータ読み取り可能な記憶媒体。 The transfer matrix of potential from the dipole in the brain to the scalp is determined by the skull model determined by the shape and conductivity of each tissue fraction of the subject's head set in advance and the arrangement of multiple electrodes on the subject's scalp. The required processing,
Processing to set the initial position data of the dipole assumed in the brain;
A process of obtaining a plurality of electrodes for electroencephalogram measurement placed on the scalp of the subject and obtaining actual measurement data converted into digital data,
Based on the transfer matrix, the initial position data, and the actual measurement data, analysis processing by the least square method is performed by parallel processing for each initial position, and a potential distribution in which a dipole assumed in the brain is generated from each parallel processing result and In real time, select the optimal dipole position and vector component in the brain that minimizes the square error with the measured potential distribution,
Each process of displaying, printing and storing selection results;
A computer-readable storage medium characterized by storing a program for analyzing an electroencephalogram dipole that performs.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005151416A JP5060025B2 (en) | 2005-05-24 | 2005-05-24 | EEG dipole analysis device, EEG dipole analysis method, EEG dipole analysis program, and storage medium storing the program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005151416A JP5060025B2 (en) | 2005-05-24 | 2005-05-24 | EEG dipole analysis device, EEG dipole analysis method, EEG dipole analysis program, and storage medium storing the program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2006325754A true JP2006325754A (en) | 2006-12-07 |
JP5060025B2 JP5060025B2 (en) | 2012-10-31 |
Family
ID=37548099
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2005151416A Expired - Fee Related JP5060025B2 (en) | 2005-05-24 | 2005-05-24 | EEG dipole analysis device, EEG dipole analysis method, EEG dipole analysis program, and storage medium storing the program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5060025B2 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103142224A (en) * | 2013-01-24 | 2013-06-12 | 上海帝仪科技有限公司 | Electroencephalogram signal detector and method for detecting electroencephalogram signal |
KR20140009715A (en) * | 2012-07-12 | 2014-01-23 | 중앙대학교 산학협력단 | Apparatus and method for determining of optimal eeg channel based on pso |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62148642A (en) * | 1985-12-23 | 1987-07-02 | 日本電気三栄株式会社 | System for analysis of living body signal |
JPH06154184A (en) * | 1992-08-05 | 1994-06-03 | Siemens Ag | Method for discrimination and classification of three-dimensional field pattern |
JPH078275B2 (en) * | 1988-07-21 | 1995-02-01 | 利光 武者 | In vivo equivalent current dipole tracker |
JPH0779804B2 (en) * | 1989-09-13 | 1995-08-30 | 利光 武者 | Equivalent dipole measuring device |
JP2536717B2 (en) * | 1993-04-21 | 1996-09-18 | 日本電気株式会社 | Dipole searcher |
JP2560651B2 (en) * | 1994-11-08 | 1996-12-04 | 日本電気株式会社 | Body condition estimation device |
JP2611960B2 (en) * | 1992-07-03 | 1997-05-21 | 利光 武者 | Dipole estimation method |
JP2626712B2 (en) * | 1988-04-08 | 1997-07-02 | 利光 武者 | In vivo equivalent current dipole display |
JPH11319A (en) * | 1997-06-13 | 1999-01-06 | Nec Corp | In vivo active region estimate device |
JP2001037729A (en) * | 1999-07-29 | 2001-02-13 | Toshiba Corp | Cardiac load test system |
JP3191862B2 (en) * | 1997-01-22 | 2001-07-23 | 本間 三郎 | A dipole tracking device in the brain by fixing the dipole estimation position of the background brain potential |
-
2005
- 2005-05-24 JP JP2005151416A patent/JP5060025B2/en not_active Expired - Fee Related
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62148642A (en) * | 1985-12-23 | 1987-07-02 | 日本電気三栄株式会社 | System for analysis of living body signal |
JP2626712B2 (en) * | 1988-04-08 | 1997-07-02 | 利光 武者 | In vivo equivalent current dipole display |
JPH078275B2 (en) * | 1988-07-21 | 1995-02-01 | 利光 武者 | In vivo equivalent current dipole tracker |
JPH0779804B2 (en) * | 1989-09-13 | 1995-08-30 | 利光 武者 | Equivalent dipole measuring device |
JP2611960B2 (en) * | 1992-07-03 | 1997-05-21 | 利光 武者 | Dipole estimation method |
JPH06154184A (en) * | 1992-08-05 | 1994-06-03 | Siemens Ag | Method for discrimination and classification of three-dimensional field pattern |
JP2536717B2 (en) * | 1993-04-21 | 1996-09-18 | 日本電気株式会社 | Dipole searcher |
JP2560651B2 (en) * | 1994-11-08 | 1996-12-04 | 日本電気株式会社 | Body condition estimation device |
JP3191862B2 (en) * | 1997-01-22 | 2001-07-23 | 本間 三郎 | A dipole tracking device in the brain by fixing the dipole estimation position of the background brain potential |
JPH11319A (en) * | 1997-06-13 | 1999-01-06 | Nec Corp | In vivo active region estimate device |
JP2001037729A (en) * | 1999-07-29 | 2001-02-13 | Toshiba Corp | Cardiac load test system |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140009715A (en) * | 2012-07-12 | 2014-01-23 | 중앙대학교 산학협력단 | Apparatus and method for determining of optimal eeg channel based on pso |
KR101999271B1 (en) | 2012-07-12 | 2019-07-11 | 중앙대학교 산학협력단 | Apparatus and method for determining of optimal eeg channel based on pso |
CN103142224A (en) * | 2013-01-24 | 2013-06-12 | 上海帝仪科技有限公司 | Electroencephalogram signal detector and method for detecting electroencephalogram signal |
Also Published As
Publication number | Publication date |
---|---|
JP5060025B2 (en) | 2012-10-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hunyadi et al. | Block term decomposition for modelling epileptic seizures | |
Viola et al. | Semi-automatic identification of independent components representing EEG artifact | |
Lee et al. | Group nonnegative matrix factorization for EEG classification | |
Ravichandran et al. | Novel tool for complete digitization of paper electrocardiography data | |
Habermehl et al. | Optimizing the regularization for image reconstruction of cerebral diffuse optical tomography | |
US20060149160A1 (en) | System and method for EEG imaging of cerebral activity using small electrode sets | |
US9271679B2 (en) | Method and apparatus for processing medical image signal | |
CN117954055A (en) | Brain functional map drawing system | |
JP2019093008A (en) | System for estimating activity of intracerebral network, method for estimating activity of intracerebral network, program for estimating activity of intracerebral network, and learned brain activity estimation model | |
Dohmatob et al. | Brain topography beyond parcellations: local gradients of functional maps | |
JP2018187044A (en) | Emotion estimation device, emotion estimation method, and computer program | |
Zebende et al. | Auto-correlation in the motor/imaginary human EEG signals: A vision about the FDFA fluctuations | |
Taberna et al. | Detection of resting-state functional connectivity from high-density electroencephalography data: impact of head modeling strategies | |
JP5060025B2 (en) | EEG dipole analysis device, EEG dipole analysis method, EEG dipole analysis program, and storage medium storing the program | |
JP2017192425A (en) | Recognition state estimation system and recognition state estimation method | |
JP2020203051A (en) | Computer program, information processing device, information processing method, leaned model generation method, and learned model | |
van der Cruijsen et al. | A method to experimentally estimate the conductivity of chronic stroke lesions: a tool to individualize transcranial electric stimulation | |
CN112164027B (en) | Spinal dural sac compression detection method, equipment and storage medium | |
RU2743608C1 (en) | Method of brain segment localization | |
KR20210154695A (en) | Brain-computer interface apparatus and operating method of reducing burden of individual calibration process by clustering subjects based on brain activation | |
Filatova et al. | Human emotions valency and level change monitoring by means of the EEG and speech signals analysis | |
JP6433616B2 (en) | Mental activity state evaluation support device, mental activity state evaluation support system, and mental activity state evaluation support method | |
JP6821171B2 (en) | Brain activity state quantification method and brain activity state measurement device | |
Schlesinger et al. | Scalp surface estimation and head registration using sparse sampling and 3D statistical models | |
Im | Basics of EEG: Generation, Acquisition, and Applications of EEG |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080522 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20110121 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110802 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20111102 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120313 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120613 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20120703 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20120803 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150810 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 Ref document number: 5060025 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |