JPH08299295A - 生体内複数電流双極子推定用の初期パラメータ設定方法 - Google Patents
生体内複数電流双極子推定用の初期パラメータ設定方法Info
- Publication number
- JPH08299295A JPH08299295A JP7109734A JP10973495A JPH08299295A JP H08299295 A JPH08299295 A JP H08299295A JP 7109734 A JP7109734 A JP 7109734A JP 10973495 A JP10973495 A JP 10973495A JP H08299295 A JPH08299295 A JP H08299295A
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- principal component
- initial parameter
- current
- living body
- 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.)
- Pending
Links
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
(57)【要約】
【目的】 計測した生体磁界に基づいてその磁界源であ
る電流双極子の位置，大きさ，向きを時間的に連続して
推定する収束計算の際の適切な初期パラメータを設定す
る。 【構成】 生体磁界波形から磁界源である電流双極子を
推定する時間帯を設定し、設定した時間帯から主成分分
析を施す複数の時間帯を決定し（Ｓ４）、決定した各時
間帯における計測データに対して主成分分析を施し（Ｓ
５）、主成分分析を行った各時間帯において、その固有
値が最大である固有ベクトルを第一主成分として選び
（Ｓ６）、選んだ複数の第一主成分により初期パラメー
タを設定する（Ｓ７）。
る電流双極子の位置，大きさ，向きを時間的に連続して
推定する収束計算の際の適切な初期パラメータを設定す
る。 【構成】 生体磁界波形から磁界源である電流双極子を
推定する時間帯を設定し、設定した時間帯から主成分分
析を施す複数の時間帯を決定し（Ｓ４）、決定した各時
間帯における計測データに対して主成分分析を施し（Ｓ
５）、主成分分析を行った各時間帯において、その固有
値が最大である固有ベクトルを第一主成分として選び
（Ｓ６）、選んだ複数の第一主成分により初期パラメー
タを設定する（Ｓ７）。
Description
【０００１】
【産業上の利用分野】本発明は、生体内における電流双
極子を推定する演算方法に関し、特に、多チャネルの超
伝導量子干渉素子（ＳＱＵＩＤ：Superconducting Quan
tum Interference Device)を用いて計測した生体磁界に
基づいて、その磁界源である電流双極子のパラメータ
（位置，大きさ，向き）を推定する際の収束計算に必要
な初期パラメータを設定する方法に関するものである。
極子を推定する演算方法に関し、特に、多チャネルの超
伝導量子干渉素子（ＳＱＵＩＤ：Superconducting Quan
tum Interference Device)を用いて計測した生体磁界に
基づいて、その磁界源である電流双極子のパラメータ
（位置，大きさ，向き）を推定する際の収束計算に必要
な初期パラメータを設定する方法に関するものである。
【０００２】
【従来の技術】人間の脳における神経活動に伴って、脳
から微弱ではあるが磁界が発生する。超伝導体を用いた
超高感度の磁気センサであるＳＱＵＩＤを用いることに
よって、脳から発生するこの微弱磁界を頭部表面にて非
侵襲的に計測することが可能である。そして、脳の計測
磁界の分布から磁界源としての電流双極子の初期パラメ
ータを所定時刻毎に間欠的に求め、得られた初期パラメ
ータに基づいて各時刻における電流双極子の位置，大き
さ，向きを推定する。図５は、この電流双極子のパラメ
ータ推定を行う実施状態を示す模式図である。被験者Ａ
の頭部には、複数の一次微分グラジオメータＭをその球
面上に配置しているＳＱＵＩＤが被せられている。各グ
ラジオメータＭからは、検出された磁界信号がＣＰＵに
送られ、ＣＰＵにて計測磁界から電流双極子のパラメー
タを推定するようになっている。
から微弱ではあるが磁界が発生する。超伝導体を用いた
超高感度の磁気センサであるＳＱＵＩＤを用いることに
よって、脳から発生するこの微弱磁界を頭部表面にて非
侵襲的に計測することが可能である。そして、脳の計測
磁界の分布から磁界源としての電流双極子の初期パラメ
ータを所定時刻毎に間欠的に求め、得られた初期パラメ
ータに基づいて各時刻における電流双極子の位置，大き
さ，向きを推定する。図５は、この電流双極子のパラメ
ータ推定を行う実施状態を示す模式図である。被験者Ａ
の頭部には、複数の一次微分グラジオメータＭをその球
面上に配置しているＳＱＵＩＤが被せられている。各グ
ラジオメータＭからは、検出された磁界信号がＣＰＵに
送られ、ＣＰＵにて計測磁界から電流双極子のパラメー
タを推定するようになっている。
【０００３】この磁界の発生源である電流双極子につい
ての位置，大きさ，向きに関するデータは、脳における
神経活動の状態を把握するのに重要である。従って、こ
れらのデータは、脳神経手術における開頭前の手術部位
の特定または脳疾患（てんかん，脳梗塞，アルツハイマ
病など）の予知，診断，治療に対して、あるいは脳の高
次機能解明などの大脳生理学に対して、有益な情報とな
り得る。
ての位置，大きさ，向きに関するデータは、脳における
神経活動の状態を把握するのに重要である。従って、こ
れらのデータは、脳神経手術における開頭前の手術部位
の特定または脳疾患（てんかん，脳梗塞，アルツハイマ
病など）の予知，診断，治療に対して、あるいは脳の高
次機能解明などの大脳生理学に対して、有益な情報とな
り得る。
【０００４】ところで、上述したような電流双極子の推
定方法では、電流双極子の推定は可能であるが、脳の対
象領域以外の他の領域の神経活動に伴って生じる磁界
（ブレインノイズ）による影響により、時間的に不連続
な推定結果しか得られない。時間の経過とともに連続的
に変化する電流双極子を推定するためには、脳磁界を時
系列信号として捉えた初期パラメータを設定する必要が
ある。
定方法では、電流双極子の推定は可能であるが、脳の対
象領域以外の他の領域の神経活動に伴って生じる磁界
（ブレインノイズ）による影響により、時間的に不連続
な推定結果しか得られない。時間の経過とともに連続的
に変化する電流双極子を推定するためには、脳磁界を時
系列信号として捉えた初期パラメータを設定する必要が
ある。
【０００５】そこで、脳磁界を時系列信号として捉え
て、主成分分析を施すことにより得られる複数の主成分
により初期パラメータを得る方法が、本発明と同一出願
人により提案されている（特開平６−311975号公報）。
この方法では、電流双極子を推定する時間帯を設定し、
当該時間帯における各チャネルの相関行列（共分散行
列）を求め、主成分分析を施すことにより、線形な相関
を抽出する方法である。電流双極子の推定時間帯に主成
分分析を施した後に、第一主成分から順に双極子性の相
関パターンになっているかを調べ、双極子性パターンで
あれば初期パラメータを設定する動作を、寄与率の累計
値が所定の閾値を超えるまで繰り返す。この方法にて複
数の電流双極子の初期パラメータを得る場合には、第一
主成分から第一電流双極子の初期パラメータを得て、第
二以降の主成分から第二以降の電流双極子の初期パラメ
ータを得ている。
て、主成分分析を施すことにより得られる複数の主成分
により初期パラメータを得る方法が、本発明と同一出願
人により提案されている（特開平６−311975号公報）。
この方法では、電流双極子を推定する時間帯を設定し、
当該時間帯における各チャネルの相関行列（共分散行
列）を求め、主成分分析を施すことにより、線形な相関
を抽出する方法である。電流双極子の推定時間帯に主成
分分析を施した後に、第一主成分から順に双極子性の相
関パターンになっているかを調べ、双極子性パターンで
あれば初期パラメータを設定する動作を、寄与率の累計
値が所定の閾値を超えるまで繰り返す。この方法にて複
数の電流双極子の初期パラメータを得る場合には、第一
主成分から第一電流双極子の初期パラメータを得て、第
二以降の主成分から第二以降の電流双極子の初期パラメ
ータを得ている。
【０００６】
【発明が解決しようとする課題】主成分分析により得ら
れる第一主成分は主たる電流双極子の磁界を反映する
が、第二主成分以降は電流双極子の位置または方向など
の変化に伴う磁界変化を反映し、必ずしも従たる電流双
極子による磁界を反映しない（栗城他：Med. Img.Tech.
Vol. 11, No. 4, September 1993）。従って、第一主
成分は電流双極子推定の初期パラメータの設定に用いる
ことができるが、第二主成分以降は初期パラメータの設
定に用いることができない場合が多いので、磁界源とし
ての複数の電流双極子の推定を行う際の初期パラメータ
の設定にこの方法を適用することは最適ではないことが
判明した。
れる第一主成分は主たる電流双極子の磁界を反映する
が、第二主成分以降は電流双極子の位置または方向など
の変化に伴う磁界変化を反映し、必ずしも従たる電流双
極子による磁界を反映しない（栗城他：Med. Img.Tech.
Vol. 11, No. 4, September 1993）。従って、第一主
成分は電流双極子推定の初期パラメータの設定に用いる
ことができるが、第二主成分以降は初期パラメータの設
定に用いることができない場合が多いので、磁界源とし
ての複数の電流双極子の推定を行う際の初期パラメータ
の設定にこの方法を適用することは最適ではないことが
判明した。
【０００７】本発明は斯かる事情に鑑みてなされたもの
であり、生体磁界の発生源である複数の電流双極子の位
置，大きさ，向きを時間的に連続して推定できるような
適切な初期パラメータを設定することが可能である生体
内複数電流双極子推定の収束計算に用いる初期パラメー
タ設定方法を提供することを目的とする。
であり、生体磁界の発生源である複数の電流双極子の位
置，大きさ，向きを時間的に連続して推定できるような
適切な初期パラメータを設定することが可能である生体
内複数電流双極子推定の収束計算に用いる初期パラメー
タ設定方法を提供することを目的とする。
【０００８】
【課題を解決するための手段】本発明に係る生体内複数
電流双極子推定用の初期パラメータ設定方法は、超伝導
量子干渉素子を用いて計測した生体磁界波形に基づい
て、生体磁界の磁界源である生体内の複数の電流双極子
を収束計算にて推定する際の初期パラメータを設定する
方法であって、複数の超伝導量子干渉素子にて計測した
生体磁界波形から電流双極子を推定する時間帯を設定
し、設定した時間帯から、複数の計測データ間の共分散
行列を導出してその共分散行列の固有値及び固有ベクト
ルを求める主成分分析を施す複数の時間帯を決定し、決
定した時間帯における計測データに主成分分析を施し、
この主成分分析により得られる、固有値が最大となる固
有ベクトルとしての複数の第一主成分により前記初期パ
ラメータを設定することを特徴とする。
電流双極子推定用の初期パラメータ設定方法は、超伝導
量子干渉素子を用いて計測した生体磁界波形に基づい
て、生体磁界の磁界源である生体内の複数の電流双極子
を収束計算にて推定する際の初期パラメータを設定する
方法であって、複数の超伝導量子干渉素子にて計測した
生体磁界波形から電流双極子を推定する時間帯を設定
し、設定した時間帯から、複数の計測データ間の共分散
行列を導出してその共分散行列の固有値及び固有ベクト
ルを求める主成分分析を施す複数の時間帯を決定し、決
定した時間帯における計測データに主成分分析を施し、
この主成分分析により得られる、固有値が最大となる固
有ベクトルとしての複数の第一主成分により前記初期パ
ラメータを設定することを特徴とする。
【０００９】
【作用】本発明では、設定した電流双極子の推定時間帯
から複数の時間帯を選択し、選択した各時間帯の計測デ
ータに主成分分析を施して得られる複数の第一主成分か
ら複数の初期パラメータを設定する。よって、時間的相
関が最も高くて電流双極子による磁界を反映する第一主
成分のみから、複数の初期パラメータを設定するので、
生体磁界を時系列信号として捉えた初期パラメータを提
供でき、かつ、複数の信号源を推定する際の初期パラメ
ータとして適当な値を与えることができる。この結果、
各電流双極子について時間的に連続した推定結果が得ら
れる。
から複数の時間帯を選択し、選択した各時間帯の計測デ
ータに主成分分析を施して得られる複数の第一主成分か
ら複数の初期パラメータを設定する。よって、時間的相
関が最も高くて電流双極子による磁界を反映する第一主
成分のみから、複数の初期パラメータを設定するので、
生体磁界を時系列信号として捉えた初期パラメータを提
供でき、かつ、複数の信号源を推定する際の初期パラメ
ータとして適当な値を与えることができる。この結果、
各電流双極子について時間的に連続した推定結果が得ら
れる。
【００１０】
【実施例】以下、本発明をその実施例を示す図面に基づ
いて説明する。なお、以下の実施例では、脳磁界を計測
する場合の例について説明する。
いて説明する。なお、以下の実施例では、脳磁界を計測
する場合の例について説明する。
【００１１】まず、各電流双極子におけるパラメータに
ついて説明する。１個の電流双極子には、その３次元的
位置，大きさ，向きの６種類のパラメータが存在する。
電流双極子が作る磁界を計算する際には、球モデルを想
定する。球の半径方向の電流双極子は球の外部に磁界を
発生しないので、その向きについては球の径方向に対す
る垂直成分のみを考慮すれば良い。従って、球モデルの
場合、各電流双極子のパラメータは以下の５種類とす
る。これらの５種類のパラメータを図２に模式的に示
す。 ｒ，θ，φ：電流双極子を極座標で表した空間的位置 ｑ：電流双極子の大きさ α：電流双極子の向き（球の経線との角度）
ついて説明する。１個の電流双極子には、その３次元的
位置，大きさ，向きの６種類のパラメータが存在する。
電流双極子が作る磁界を計算する際には、球モデルを想
定する。球の半径方向の電流双極子は球の外部に磁界を
発生しないので、その向きについては球の径方向に対す
る垂直成分のみを考慮すれば良い。従って、球モデルの
場合、各電流双極子のパラメータは以下の５種類とす
る。これらの５種類のパラメータを図２に模式的に示
す。 ｒ，θ，φ：電流双極子を極座標で表した空間的位置 ｑ：電流双極子の大きさ α：電流双極子の向き（球の経線との角度）
【００１２】図１は、本発明の生体内複数電流双極子推
定用の初期パラメータ設定方法の手順を示すフローチャ
ートであり、脳磁界における初期パラメータ設定の例を
表している。以下、このフローチャートに従って設定手
順を説明する。
定用の初期パラメータ設定方法の手順を示すフローチャ
ートであり、脳磁界における初期パラメータ設定の例を
表している。以下、このフローチャートに従って設定手
順を説明する。
【００１３】まず、計測対象である被験者の頭部を球と
見なして、その中心座標及び半径ｒ 0 を求め、有限時間
長にわたって、図５に示すようにＳＱＵＩＤを用いて、
被験者の脳磁界をｎ点で計測する（ステップＳ１）。そ
して、得られたｎ本の脳磁界波形を、基線を合わせて一
括表示する（ステップＳ２）。表示された波形から電流
双極子を推定する時間帯を決定し、その時間帯をｍ分割
する（ステップＳ３）。時刻Ｔk におけるｊ番目のＳＱ
ＵＩＤで得られる計測磁界をＢｍj （Ｔk ）とする（ｊ
＝１，２，…，ｎ, ｋ＝１，２，…，ｍ）。
見なして、その中心座標及び半径ｒ 0 を求め、有限時間
長にわたって、図５に示すようにＳＱＵＩＤを用いて、
被験者の脳磁界をｎ点で計測する（ステップＳ１）。そ
して、得られたｎ本の脳磁界波形を、基線を合わせて一
括表示する（ステップＳ２）。表示された波形から電流
双極子を推定する時間帯を決定し、その時間帯をｍ分割
する（ステップＳ３）。時刻Ｔk におけるｊ番目のＳＱ
ＵＩＤで得られる計測磁界をＢｍj （Ｔk ）とする（ｊ
＝１，２，…，ｎ, ｋ＝１，２，…，ｍ）。
【００１４】波形データから電流双極子の数Ｄを決定
し、その数だけ主成分分析を行う時間帯（Ｔpsa 〜Ｔpe
a ）（ａ＝１，２，…，Ｄ）を決定する（ステップＳ
４）。各々の時間帯に含まれるＴk をＴａg とし、その
数をｍｐa とする（ｇ＝１，２，…，ｍｐa ）。
し、その数だけ主成分分析を行う時間帯（Ｔpsa 〜Ｔpe
a ）（ａ＝１，２，…，Ｄ）を決定する（ステップＳ
４）。各々の時間帯に含まれるＴk をＴａg とし、その
数をｍｐa とする（ｇ＝１，２，…，ｍｐa ）。
【００１５】それぞれの時間帯で得られた（ｍｐa ×
ｎ）個の脳磁界データに対して主成分分析を施す（ステ
ップＳ５）。すなわち、個々のＳＱＵＩＤで得られた脳
磁界間の共分散行列σａを導出し、共分散行列の固有値
λａi ，固有ベクトルｖａi （ｉ＝１，２，…，ｎ）を
求める。なお、共分散行列の各要素σａpqは具体的には
以下のように求める。
ｎ）個の脳磁界データに対して主成分分析を施す（ステ
ップＳ５）。すなわち、個々のＳＱＵＩＤで得られた脳
磁界間の共分散行列σａを導出し、共分散行列の固有値
λａi ，固有ベクトルｖａi （ｉ＝１，２，…，ｎ）を
求める。なお、共分散行列の各要素σａpqは具体的には
以下のように求める。
【００１６】
【数１】
【００１７】主成分分析を行ったそれぞれの時間帯にお
いて、求めた固有ベクトルｖａi の中で、その固有値λ
ａi が最大のものを選ぶ（ステップＳ６）。固有値が最
大である選び出した固有ベクトルを第一主成分ｖａ0 と
呼び、これをもとに、各々について初期パラメータ（ｒ
a ，θa ，φa ，ｑa ，αa ）を設定する（ステップＳ
７）。初期パラメータの設定は次のように行う。
いて、求めた固有ベクトルｖａi の中で、その固有値λ
ａi が最大のものを選ぶ（ステップＳ６）。固有値が最
大である選び出した固有ベクトルを第一主成分ｖａ0 と
呼び、これをもとに、各々について初期パラメータ（ｒ
a ，θa ，φa ，ｑa ，αa ）を設定する（ステップＳ
７）。初期パラメータの設定は次のように行う。
【００１８】ａ番目の時間帯のデータに主成分分析を施
した結果の第一主成分ｖａ0 （ａ＝１，２，…，Ｄ）の
すべての値における最大値及び最小値とこれらの最大
値，最小値を示す計測点（２点）の位置とに注目し、こ
の２点の位置，磁界強度により、各パラメータの初期値
を以下のようにして決定する。
した結果の第一主成分ｖａ0 （ａ＝１，２，…，Ｄ）の
すべての値における最大値及び最小値とこれらの最大
値，最小値を示す計測点（２点）の位置とに注目し、こ
の２点の位置，磁界強度により、各パラメータの初期値
を以下のようにして決定する。
【００１９】
【数２】
【００２０】θa ，φa ：２点間の中点の座標から決定 αa ：２点間を結んだ線に垂直な方向とし、２点間にお
ける磁界の符号を満たす向き ｑa ：以上の４種のパラメータにて決まる位置，向きに
電流双極子を設置した場合に、２点の磁界強度を満たす
大きさ
ける磁界の符号を満たす向き ｑa ：以上の４種のパラメータにて決まる位置，向きに
電流双極子を設置した場合に、２点の磁界強度を満たす
大きさ
【００２１】以上のように、本発明では、複数の各時間
帯において得られる第一主成分に基づいて、複数の電流
双極子の初期パラメータを設定している。第一主成分は
時間的相関が最も高く電流双極子による磁界を反映する
ので、この第一主成分を用いることにより、第２電流双
極子以降についても適切な初期パラメータを設定でき
る。
帯において得られる第一主成分に基づいて、複数の電流
双極子の初期パラメータを設定している。第一主成分は
時間的相関が最も高く電流双極子による磁界を反映する
ので、この第一主成分を用いることにより、第２電流双
極子以降についても適切な初期パラメータを設定でき
る。
【００２２】次に、本発明における初期パラメータ設定
の具体例について、図３，図４を参照して説明する。Ｓ
ＱＵＩＤの総数（計測点数）が37個である脳磁計で計測
した磁界の一括表示を図３に示している。設定した推定
時間帯（図３中Ａ）において磁界のピークが２つあるの
で、電流双極子の数を２個と決定する（Ｄ＝２）。ここ
では、時間帯Ａを前後半２分割した時間帯Ａ１，Ａ２を
主成分分析を施す時間帯として決定し、それぞれの時間
帯Ａ１，Ａ２における計測データに主成分分析を施す。
その主成分分析結果の第一主成分について等磁界線を描
いたものを図４（ａ），（ｂ）に示す。図４（ａ），
（ｂ）はそれぞれ時間帯Ａ１，Ａ２における等磁界線を
示している。また、図４（ａ），（ｂ）において、実
線，破線はそれぞれ磁界の湧き出し，磁界の吸い込みを
示し、太い実線は磁界が０であることを示し、また●は
計測点を示す。
の具体例について、図３，図４を参照して説明する。Ｓ
ＱＵＩＤの総数（計測点数）が37個である脳磁計で計測
した磁界の一括表示を図３に示している。設定した推定
時間帯（図３中Ａ）において磁界のピークが２つあるの
で、電流双極子の数を２個と決定する（Ｄ＝２）。ここ
では、時間帯Ａを前後半２分割した時間帯Ａ１，Ａ２を
主成分分析を施す時間帯として決定し、それぞれの時間
帯Ａ１，Ａ２における計測データに主成分分析を施す。
その主成分分析結果の第一主成分について等磁界線を描
いたものを図４（ａ），（ｂ）に示す。図４（ａ），
（ｂ）はそれぞれ時間帯Ａ１，Ａ２における等磁界線を
示している。また、図４（ａ），（ｂ）において、実
線，破線はそれぞれ磁界の湧き出し，磁界の吸い込みを
示し、太い実線は磁界が０であることを示し、また●は
計測点を示す。
【００２３】それぞれの第一主成分から最大値，最小値
を示す計測点（図４（ａ），（ｂ）で点Ｌ：最大値を示
す計測点，点Ｓ：最小値を示す計測点）を求め、これら
の両点のデータに基づいて１つずつ電流双極子の初期パ
ラメータを設定する（図４（ａ），（ｂ）中のＱa,Ｑ
b)。
を示す計測点（図４（ａ），（ｂ）で点Ｌ：最大値を示
す計測点，点Ｓ：最小値を示す計測点）を求め、これら
の両点のデータに基づいて１つずつ電流双極子の初期パ
ラメータを設定する（図４（ａ），（ｂ）中のＱa,Ｑ
b)。
【００２４】
【発明の効果】以上のように、本発明の生体内複数電流
双極子推定用の初期パラメータ設定方法では、設定した
電流双極子の推定時間帯から複数の時間帯を選択して主
成分分析を施し、得られた複数の第一主成分からの複数
の初期パラメータを設定するので、生体磁界を時系列信
号として捉えた初期パラメータを提供でき、かつ、複数
の電流双極子を推定する際の初期パラメータとして適当
な値を与えることができる。この結果、本発明を利用す
れば、各電流双極子について時間的に連続した推定結果
が得られ、神経活動が複雑である脳内の電流双極子につ
いてより正確なデータを得ることができ、脳疾患部位の
特定または脳の高次機能の解明に対して有効な情報を提
供することが可能となる。
双極子推定用の初期パラメータ設定方法では、設定した
電流双極子の推定時間帯から複数の時間帯を選択して主
成分分析を施し、得られた複数の第一主成分からの複数
の初期パラメータを設定するので、生体磁界を時系列信
号として捉えた初期パラメータを提供でき、かつ、複数
の電流双極子を推定する際の初期パラメータとして適当
な値を与えることができる。この結果、本発明を利用す
れば、各電流双極子について時間的に連続した推定結果
が得られ、神経活動が複雑である脳内の電流双極子につ
いてより正確なデータを得ることができ、脳疾患部位の
特定または脳の高次機能の解明に対して有効な情報を提
供することが可能となる。
【図１】本発明に係る生体内複数電流双極子推定用の初
期パラメータ設定方法の手順を示すフローチャートであ
る。
期パラメータ設定方法の手順を示すフローチャートであ
る。
【図２】電流双極子の５個のパラーメータを示す模式図
である。
である。
【図３】本発明の具体例を説明するための計測磁界波形
図である。
図である。
【図４】本発明の具体例を示すための主成分分析結果の
第一主成分の等磁界線図である。
第一主成分の等磁界線図である。
【図５】ＳＱＵＩＤを用いて脳磁界の電流双極子のパラ
メータ推定を行う実施状態を示す模式図である。
メータ推定を行う実施状態を示す模式図である。
Claims (1)
- 【請求項１】 超伝導量子干渉素子を用いて計測した生
体磁界波形に基づいて、生体磁界の磁界源である生体内
の複数の電流双極子を収束計算にて推定する際の初期パ
ラメータを設定する方法であって、複数の超伝導量子干
渉素子にて計測した生体磁界波形から電流双極子を推定
する時間帯を設定し、設定した時間帯から、複数の計測
データ間の共分散行列を導出してその共分散行列の固有
値及び固有ベクトルを求める主成分分析を施す複数の時
間帯を決定し、決定した時間帯における計測データに主
成分分析を施し、この主成分分析により得られる、固有
値が最大となる固有ベクトルとしての複数の第一主成分
により前記初期パラメータを設定することを特徴とする
生体内複数電流双極子推定用の初期パラメータ設定方
法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7109734A JPH08299295A (ja) | 1995-05-08 | 1995-05-08 | 生体内複数電流双極子推定用の初期パラメータ設定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7109734A JPH08299295A (ja) | 1995-05-08 | 1995-05-08 | 生体内複数電流双極子推定用の初期パラメータ設定方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH08299295A true JPH08299295A (ja) | 1996-11-19 |
Family
ID=14517896
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP7109734A Pending JPH08299295A (ja) | 1995-05-08 | 1995-05-08 | 生体内複数電流双極子推定用の初期パラメータ設定方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH08299295A (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6073040A (en) * | 1997-01-20 | 2000-06-06 | Nec Corporation | Electrophysiological activity estimation method |
| US6226544B1 (en) | 1997-06-13 | 2001-05-01 | Nec Corporation | Living body internal active source estimation apparatus |
-
1995
- 1995-05-08 JP JP7109734A patent/JPH08299295A/ja active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6073040A (en) * | 1997-01-20 | 2000-06-06 | Nec Corporation | Electrophysiological activity estimation method |
| US6226544B1 (en) | 1997-06-13 | 2001-05-01 | Nec Corporation | Living body internal active source estimation apparatus |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Uutela et al. | Global optimization in the localization of neuromagnetic sources | |
| US5307807A (en) | Method and system for three dimensional tomography of activity and connectivity of brain and heart electromagnetic waves generators | |
| US4977896A (en) | Analysis of biological signals using data from arrays of sensors | |
| US6741888B2 (en) | Method and apparatus for estimating degree of neuronal impairment in brain cortex | |
| Stenbacka et al. | Comparison of minimum current estimate and dipole modeling in the analysis of simulated activity in the human visual cortices | |
| Kobayashi et al. | Systematic source estimation of spikes by a combination of independent component analysis and RAP-MUSIC: I: Principles and simulation study | |
| CN114065825B (zh) | 一种基于结构相似性的脑磁meg源定位方法 | |
| Ioannides | Estimates of Brain Activity. Using Magnetic Field Tomography and | |
| JPH04210046A (ja) | 生体シグナルの解析法 | |
| JPH08299295A (ja) | 生体内複数電流双極子推定用の初期パラメータ設定方法 | |
| JP3519141B2 (ja) | 生体内電流源推定方法 | |
| Im et al. | Anatomically constrained dipole adjustment (ANACONDA) for accurate MEG/EEG focal source localizations | |
| Yetik et al. | Line-source modeling and estimation with magnetoencephalography | |
| JPH06311975A (ja) | 生体内複数電流双極子推定用の初期パラメータ設定方法 | |
| US7398120B2 (en) | Method and arrangement and computer programme with programme code means and computer programme products for the analysis of neuronal activities in neuronal areas | |
| JP3489079B2 (ja) | 三次元生体情報解析装置 | |
| Khan et al. | Encoding cortical dynamics in sparse features | |
| Shi et al. | Magnetoencephalography for epileptic focus localization based on Tucker decomposition with ripple window | |
| Yetik et al. | Surface-source modeling and estimation using biomagnetic measurements | |
| Fell et al. | Linear inverse filtering improves spatial separation of nonlinear brain dynamics: a simulation study | |
| Mino et al. | An estimation algorithm for neuromagnetic source distribution using MRI information | |
| Mino et al. | A neuromagnetic source distribution estimation using MRI information | |
| Uchikawa et al. | A 3-D vector measurement and its application on biomagnetic signals | |
| JP2003038455A (ja) | 生体活動解析装置 | |
| Saita et al. | An estimation algorithm of neuromagnetic sources in the cortical region using realistically-shaped head model |