JPH10295660A - 生体活動電流源推定装置 - Google Patents

Abstract

(57)【要約】 【課題】 短い診断期間で生体活動電流源を推定する
ことができる生体活動電流源推定装置を提供する。 【解決手段】 被検体の関心部位に装着した発信器から
の磁界を計測し（Ｓ１）、この磁界データから発信器の
位置を特定し（Ｓ２）、この発信器の位置に基づいた球
モデルを設定する（Ｓ３）。次に、関心部位の電流源の
磁界計測を行い（Ｓ４）、前記球モデルに対して電流源
を推定する（Ｓ５）。この電流源の位置、向き、大きさ
および計測磁界データを波形表示する（Ｓ６）。したが
って、関心部位を計測した計測磁界データが計測後、直
ちに有効か否かがわかるので、効率よく診断することが
できる。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】この発明は、被検体の関心部
位の電流源の磁界計測から電流源の位置、向き、大きさ
の推定までを行う生体活動電流源推定装置に係り、特
に、生体活動電流源の位置、向き、大きさの推定を効率
良く行なう生体活動電流源推定装置に関する。

【従来の技術】生体に刺激を与えると、細胞膜を挟んで
形成されている分極が壊れて生体活動電流が流れる。こ
の生体活動電流は、脳や心臓において現れ、脳波，心電
図として記録される。また、生体活動電流によって生じ
る磁界は、脳磁図，心磁図として記録される。

【０００２】近年、生体内の微小な磁界を計測する装置
として、ＳＱＵＩＤ（Superconduc-ting Quantum Inter
ference Device：超電導量子干渉計）を用いた、マルチ
チャンネルＳＱＵＩＤセンサが開発されている。このマ
ルチチャンネルＳＱＵＩＤセンサは、デュアーと呼ばれ
る容器内に多数の磁気センサを液体窒素などの冷媒に浸
漬して収納している。

【０００３】従来、まず、このマルチチャンネルＳＱＵ
ＩＤセンサを関心部位である例えば頭部外側に置き、脳
内に生じた生体活動電流源である電流双極子（以下、
「電流源」という）からの微小磁界をマルチチャンネル
ＳＱＵＩＤセンサで無侵襲に計測する。次に、ＭＲＩ装
置を用いて、関心部位である頭部を撮像することで、三
次元の形態情報が得られる。これらの形態情報から例え
ば、特開平４ー２２６６３５号公報に示す方法等を利用
して、関心部位の立体像が作成される。この関心部位の
立体像に基づいて最小自乗法を用いた方法により立体像
全体に合う球モデルが設定される。

【０００４】被検体の関心部位の電流源の微小磁界を計
測した計測磁界データと前記球モデルとに基づいて、最
小自乗法や最小ノルム法等を用いた推定方法を利用し、
電流源の位置、大きさ、向きを推定する。なお、電流源
の位置、向き、大きさは、球モデルを基準として推定さ
れている。また、この球モデルは、関心部位の立体像か
ら算出されているので、電流源の位置、向き、大きさを
脳の断層画像に重ねることができる。これにより脳のど
の位置に電流源があるかを特定することができる。

【０００５】

【発明が解決しようとする課題】しかしながら、このよ
うな構成を有する従来例の場合には、次のような問題が
ある。球モデルを設定するためには、ＭＲＩ装置を利用
して被検体の関心部位を撮像する必要がある。ＭＲＩ装
置は、被検体の関心部位に強い磁場をかけて、被検体の
関心部位を一時的に磁化することで撮像をしている。こ
の磁化された状態の間、被検体の関心部位に対して磁界
計測を行うことができないので、十分な時間が経過した
後に再度磁界計測を開始する必要がある。このような理
由で、磁界計測前に電流源を推定するための球モデルを
設定することができなかった。

【０００６】また、被検体の関心部位の電流源から発生
する磁界は、非常に微小な磁界である。この微小磁界
は、被検体の関心部位以外および周囲の磁場などのノイ
ズの影響を受け易い。このノイズを含む計測磁界データ
に基づき、電流源を推定すると、被検体の関心部位に対
する電流源の位置、向き、大きさを正確に推定すること
ができず、再度、磁界計測からやり直さなければならな
いので、診断期間が長期化する問題があった。

【０００７】このような事情に鑑みてなされたものであ
って、電流源の磁界計測から電流源の推定までを短い診
断期間で効率良く行う生体活動電流源推定装置の提供す
ることを目的とする。

【０００８】

【課題を解決するための手段】この発明は、このような
目的を達成するために、次のような構成をとる。すなわ
ち、請求項１に記載の発明は、被検体の関心部位の生体
活動電流源から生じる微小磁界を、前記被検体の関心部
位に近接配備された磁気センサで計測するとともに、画
像撮像装置で得られた三次元の形態情報に基づき、球モ
デルを設定し、この球モデルに対して、前記磁気センサ
で得られた計測磁界データに基づいて、前記生体活動電
流源の位置、向き、大きさを推定する生体活動電流源推
定装置において、（ａ）前記被検体の関心部位の表面に
装着された複数個の発信器と、（ｂ）前記各発信器から
発生される磁界データに基づき、前記磁気センサを基準
とした座標系における前記各発信器の位置を得る発信器
位置特定手段と、（ｃ）前記発信器位置特定手段で得ら
れた前記磁気センサを基準とした座標系における全発信
器の位置に基づき、球モデルを設定する球モデル設定手
段と、（ｄ）前記発信器の位置に基づく球モデルに対し
て、前記計測磁界データに基づき電流源の位置、向き、
大きさを推定する電流源推定手段と、（ｅ）前記推定さ
れた電流源の結果を表示する表示装置とを備えたことを
特徴とするものである。

【０００９】請求項２に記載の発明は、被検体の関心部
位の生体活動電流源から生じる微小磁界を、前記被検体
の関心部位に近接配備された磁気センサで計測するとと
もに、画像撮像装置で得られた三次元の形態情報に基づ
き、球モデルを設定し、この球モデルに対して、前記磁
気センサで得られた計測磁界データに基づいて、前記生
体活動電流源の位置、向き、大きさを推定する生体活動
電流源推定装置において、（ａ）前記被検体の関心部位
の表面に装着された複数個の発信器と、（ｂ）前記各発
信器から発生される磁界データに基づき、前記磁気セン
サを基準とした座標系における前記各発信器の位置を得
る発信器位置特定手段と、（ｃ）前記発信器位置特定手
段で得られた前記磁気センサを基準とした座標系におけ
る全発信器の位置に基づき、球モデルを設定する球モデ
ル設定手段と、（ｄ）前記発信器の位置に基づく球モデ
ルに対して、前記計測磁界データに基づき電流源の位
置、向き、大きさを推定する電流源推定手段と、（ｆ）
前記推定された電流源の位置、向きのデータから得られ
る磁界を推測し、この推測磁界データと前記計測磁界デ
ータとの一致する度合いを示す値を算出する一致度算出
手段と、（ｇ）前記一致度算出手段で算出された値と、
予め定められた値とを比較し、前記一致する度合いを示
す値が前記予め定められた値よりも大きくなるまで、上
述した被検体の関心部位の磁界計測、この磁界データに
基づく電流源の推定、および一致する度合いを示す値の
算出を繰り返し実行させる判断手段とを備えたことを特
徴とするものである。

【００１０】〔作用〕この発明の作用は次のとおりであ
る。請求項１の発明は、まず、被検体の関心部位の表面
に複数個の発信器を装着する。この各発信器を順次磁界
を発信させて、この磁界を磁気センサで計測する。この
各発信器の磁界データに基づいて、各発信器の位置をそ
れぞれ特定する。この特定された発信器の位置は、磁気
センサを基準とした座標系で表される。この特定された
全発信器の位置に基づいて、球モデルを設定する。次
に、被検体の関心部位の電流源の磁界計測を行い、この
球モデルに対して、電流源を推定する。この推定された
電流源の結果を表示装置に表示する。

【００１１】請求項２の発明は、まず、被検体の関心部
位の表面に複数個の発信器を装着する。この各発信器を
順次磁界を発信させて、この磁界を磁気センサで計測す
る。この各発信器の磁界データに基づいて、各発信器の
位置をそれぞれ特定する。この特定された発信器の位置
は、磁気センサを基準とした座標系で表される。この特
定された全発信器の位置に基づいて、球モデルを設定す
る。次に、被検体の関心部位の電流源の磁界計測を行
い、この球モデルに対して、電流源を推定する。この推
定された電流源の位置、向きのデータに基づき磁界デー
タを推測する。この推測磁界データを前記計測磁界デー
タとの一致の度合いを示す値と、予め定められた値とを
比べることで、計測磁界データが有効か否かを判断す
る。そして、ＭＲＩ装置によって撮像された被検体の関
心部位の三次元の形態情報に基づき球モデルを設定す
る。この球モデルに対し最小ノルム法などを用いた手法
でこの有効な計測磁界データに基づいて、被検体の関心
部位の電流源の位置、向き、大きさを推定する。

【００１２】

【発明の実施の形態】

＜第１実施例＞図１を参照して実施例に係る生体活動電
流源推定装置の概略構成を説明する。本実施例では、被
検体の関心部位を頭部として電流源を推定する。図中、
符号１はマルチチャンネルＳＱＵＩＤセンサを示す。こ
のマルチチャンネルＳＱＵＩＤセンサ１は、被検体Ｍの
頭部の外側に配備されている。マルチチャンネルＳＱＵ
ＩＤセンサ１は、デュアー１ａ内に図示しない多数の磁
気センサＳ₁ 〜Ｓ_m が冷媒に侵漬して収納されて構成さ
れている。これらマルチチャンネルＳＱＵＩＤセンサ１
や被検体Ｍは、磁気シールドルーム２内に配設される。

【００１３】被検体Ｍの頭部表面には、複数個の発信コ
イル２３が装着されている。この実施例では、４個の発
信コイル２３が図示されているが、これは球モデルを作
成する上で４個以上あれば良く、発信コイル２３の数は
多い程良い。また、発信コイル２３は被検体の関心部位
である頭部表面の図示する位置に限定するものではな
く、関心部位の形状に合わせて装着できればよい。した
がって、関心部位に発信コイル２３を直に貼り付けるも
のに限らず、例えば、関心部位を覆うマスクの表面にこ
の発信コイル２３を装着して、このマスクを関心部位に
被せてもよい。

【００１４】各発信コイル２３は、例えば、図２（ａ）
に示すように、セラミック板等の基板３１に金属を印刷
してコイル部３２を形成したコイルや、図２（ｂ）に示
すように、ボビン３３に金属ワイヤ３４を巻回して形成
されたコイルで構成される。なお、図２（ａ）の構成の
コイルを用いた場合には、コイル部３２をより真円に近
い形状に構成でき、また、複数個の発信コイル２３の各
コイル部３２の形状を揃えることができるので、後述す
る発信コイル２３の位置検出の計算を容易にかつ正確に
行うことができ、また、図２（ａ）の構成のコイルは図
２（ｂ）の構成のコイルよりもかさばらない。従って、
発信コイル２３としては、図２（ａ）の構成のコイルを
用いるのが好ましい。

【００１５】各発信コイル２３には、電流供給制御部７
から個別に電流が供給されるように構成されている。そ
して、この電流供給制御部７から、各発信コイル２３に
１個ずつ、電流を順次個別に供給されていき、各発信コ
イル２３から生じる磁界を、マルチチャネルＳＱＵＩＤ
センサ１の磁気センサＳ₁ 〜Ｓ_m で検出する。この各発
信コイル２３からの磁界の計測磁界データは、データ変
換部４でデジタルデータに変換された後、データ収集部
５に集められる。

【００１６】また、被検体Ｍには、刺激装置６から電気
的刺激（あるいは音、光刺激など）が与えられ、この刺
激により、被検体Ｍの脳内に発生した生体活動電流源に
より生じる微小磁界も、マルチチャネルＳＱＵＩＤセン
サ１の磁気センサＳ₁ 〜Ｓ_mで検出され、データ変換部
４でデジタルデータに変換された後、データ収集部５に
集められる。

【００１７】データ解析部９は、データ収集部５で集め
られた発信コイル２３の計測磁界データに基づいて、被
検体Ｍの頭部の球モデルを設定する。次に電流源からの
計測磁界データに基づいて、この球モデルに対する電流
源を推定するためのものである。また、この球モデルに
対して推定された電流源の位置、大きさ、向きをこの発
明の表示装置であるカラーモニタ１０やカラープリンタ
１１に出力するようになっている。

【００１８】また、この発明の指示手段であるキーボー
ド１２は、カラーモニタ１０やカラープリンタ１１に出
力された電流源の結果を確認してから、処理を続行する
か否かを指示する。処理を続行しない場合はここで処理
を終了する。一方、処理を続行する場合は、いったん処
理を中断し、ＭＲＩ装置を用いて被検体Ｍの関心部位を
撮像して得られた三次元の形態情報を光磁気ディスク１
３に記憶した後、処理を再開する。

【００１９】処理を再開したデータ解析部９は、この光
磁気ディスク１３に記憶されている被検体Ｍの頭部の三
次元の形態情報に基づいて、球モデルを設定する。続い
て、計測磁界データに基づいて、この球モデルに対する
電流源を推定する。この推定された電流源は、三次元の
形態情報の一部である断層画像上に重ね合わされてカラ
ーモニタ１０に表示されたり、カラープリンタ１１に印
字されるようになっている。なお、ＭＲＩ装置で撮像さ
れた三次元の形態情報は、通信回線１４を介してデータ
解析部８に直接伝送するように構成してもよい。

【００２０】以下、この実施例に係る処理手順を図３の
フローチャートを参照して詳細に説明する。まず、マル
チチャネルＳＱＵＩＤセンサ１（デュアー１ａに収納さ
れた図示しない磁気センサＳ₁ 〜Ｓ_m ）を基準とした座
標系における各発信コイル２３の位置を特定する。これ
は、電流供給制御部７が、各発信コイル２３に１個ず
つ、電流を順次個別に供給し、これにより、各発信コイ
ル２３から順次個別に発生される磁界を、電流供給タイ
ミングに応じて、マルチチャネルＳＱＵＩＤセンサ１内
の各磁気センサＳ₁ 〜Ｓ_m で計測する（ステップＳ
１）。

【００２１】この計測磁界データは、データ変換部４を
経てデータ収集部５で収集される。そして、データ収集
部５に収集された各発信コイル２３からの磁界の計測磁
界データは、データ解析部９に取り込まれる。データ解
析部９では、最小自乗法等によって、各発信コイル２３
の位置を求める（ステップＳ２）。

【００２２】このステップＳ２で行われる、最小自乗法
による各発信コイル２３の位置の特定手順を図４に示
す。例えば、第ｋ番目（ｋは、１〜ｎ個の発信コイル２
３のｋ番目を示す自然数）の発信コイル２３の位置は、
まず、その発信コイル２３の初期位置を仮想的に決める
（ステップＴ２）。次に、各磁気センサＳ₁ 〜Ｓ_m が、
仮想的に決めた位置にある発信コイル２３からの磁界を
計測したと仮定した場合の仮想的な仮想磁界データを計
算で求める（ステップＴ３）。この計算で求めたデータ
と、各磁気センサＳ₁ 〜Ｓ_m が第ｋ番目の発信コイル２
３から実際に計測した計測磁界データとの自乗誤差を求
める（ステップＴ４）。第ｋ番目の発信コイル２３の仮
想的な位置を、自乗誤差が小さくなる方向へ移動する
（ステップＴ５）。自乗誤差が予め決めておいた終了判
定値以下であれば、その位置を第ｋ番目の発信コイル２
３の位置とし、一方、自乗誤差が終了判定値よりも大き
ければ、上記ステップＴ３〜Ｔ５の処理を繰り返す（ス
テップＴ６）。このようにして、第ｉ番目の発信コイル
２３の位置を求め、同様に、全発信コイル２３の位置を
求める（ステップＴ１、Ｔ７、Ｔ８）。なお、各発信コ
イル２３のコイル部が真円に近く、各発信コイル２３の
コイル部の形状が揃っていると、上記ステップＴ３の計
算が容易で正確になる。このようにして求められた各発
信コイル２３の位置は、マルチチャネルＳＱＵＩＤセン
サ１（磁気センサＳ₁ 〜Ｓ_m ）を基準とした三次元の座
標系における位置である。

【００２３】次に、最小自乗法を含む関数を用いて、こ
れら各発信コイル２３の位置に基づき、被検体Ｍの関心
部位である頭部の球モデルを設定する（ステップＳ
３）。

【００２４】ステップＳ３での球モデル設定手順を説明
する。ステップＳ２で求められたｎ個の発信コイル２３
の各位置は、三次元の座標点で表されるので、その各発
信コイル２３の位置座標は、Ｐ_i ＝（ｘ_i , ｙ_i ，
ｚ_i ）（i=1,2,…,n) で表すことができる。ここで、求
めたい球モデルの中心座標をＰ＝（Ｘ，Ｙ，Ｚ）、半径
をｒとする。次に、球モデルの中心点Ｐから各発信コイ
ル２３の位置Ｐ_i までの距離をＲ_i (i=1,2, …,n) とす
ると、Ｒ_i は次式（１）で表される。

【００２５】 Ｒ_i ＝√（（Ｘ−ｘ_i ）² ＋（Ｙ−ｙ_i ）² ＋（Ｚ−ｚ_i ）² ）‥‥‥（１） 中心点Ｐから各発信コイル２３の位置Ｐ_i までの距離Ｒ
_i と、求めたい球モデルの半径ｒとの自乗差を最小とす
るような中心点の座標Ｐ＝（Ｘ，Ｙ，Ｚ）と球モデルの
半径ｒを求めることで、全発信コイル２３の位置に整合
した球モデルが設定することができる。すなわち、Ｒ_i
とｒの最小自乗差を求める。ここで、最小自乗法を用い
た関数ｅとして、次式（２）が用いられる。

【００２６】

【数１】

【００２７】（２）式において、関数ｅを最小とするた
めの条件を次式（３）〜（６）に示す。

【００２８】

【数２】

【００２９】ここで、（３）〜（６）式で与えられた連
立方程式の解を求めることで球モデルの中心点の座標Ｐ
＝（Ｘ，Ｙ，Ｚ）と、その半径ｒとが求まり、発信コイ
ル２３の位置に基づいた球モデルが設定される。

【００３０】次に、刺激装置６から被検体Ｍに刺激を与
え、脳内に生体活動電流源から生じる微小磁界をマルチ
チャネルＳＱＵＩＤセンサ１内の各磁気センサＳ₁ 〜Ｓ
_m で計測する（ステップＳ４）。

【００３１】この計測磁界データは、データ変換部４を
経て、データ収集部５で収集され、データ解析部９に送
られる。データ解析部９では、最小自乗法や最小ノルム
法などの従来から用いられてきた方法により、脳内の生
体活動電流源の位置、向き、大きさを推定する。なお、
ここで、マルチチャネルＳＱＵＩＤセンサ１（磁気セン
サＳ₁ 〜Ｓ_m ）に対する被検体Ｍの頭部の球モデルが用
いられる。この球モデルは、上述したように既に求めた
各発信コイル２３の位置座標に基づいて設定された球モ
デルを用いる。（ステップＳ５）

【００３２】この電流源の推定結果や各磁気センサで計
測された計測磁界データの波形表示などを、この発明の
表示装置であるカラーモニタ１０やカラープリンタ１１
に出力する（ステップＳ６）。

【００３３】＜第２実施例＞第２実施例において、装置
の概略構成は図１に示した第１実施例を同様であり、ま
た、第１実施例の主要部であるＳ１〜Ｓ５までの各手段
による処理は共通するので、ここでの説明は省く。以
下、図５のフローチャートに示すステップＶ６〜Ｖ９を
参照して、第２実施例に係る、データ解析部９で実行さ
れる処理について説明する。

【００３４】ステップＳ５で求められた電流源の位置、
向き、大きさのデータの内、電流源の位置と向きのデー
タに基づき、磁界データを推測する。ここでは、前記電
流源の位置と向きのデータに対してビオ・サバールの法
則を適用することでマルチチャンネルＳＱＵＩＤセンサ
１で計測されると推測される磁界データを得ることがで
きる。この推測磁界データとステップＳ４で得た電流源
の計測磁界データとが一致する度合いを示す値（以下、
「一致度」と呼ぶ）によって判断する。

【００３５】前記推定磁界データをｏ_i として前記計測
磁界データをｅ_i とする。このｏ_iとｅ_i とが互いに一
致しているか否かを一致度から判断する。この両データ
の一致を見るには種々の手法が存在するが、例えば、次
式（７）に基づき得られたγの値を指標とすることがで
きる。

【００３６】 γ＝Σ（ｏ_i ・ｅ_i ）／（√（Σｏ_i ² ）・√（Σｅ_i ² ））・・・（７） ここで、γは、一致度で、ー１≦γ≦１の値をとり、そ
の値が１に近いほど、両データは、一致する度合いが高
いことを示す。ここで、式（７）にｏ_i とｅ_iを代入し
一致度γを算出する。（ステップＶ６）

【００３７】この一致度γが、予め定められた値、例え
ば０．９以上の場合は、両データは一致する度合いが高
く、計測された計測データには、有効なデータであると
判断できる。一方、一致度γが、０．９より小さい場合
は、有効なデータでないと判断する。ここで無効なデー
タをして判断された場合は、再度、Ｓ４〜Ｖ６の各処理
を行い、有効なデータと判断された場合は、次のステッ
プＶ８に進む。（ステップＶ７）。

【００３８】従来と同様に、ＭＲＩ装置を用いて被検体
Ｍの関心部位である頭部の撮像を行う。この撮像により
得られた三次元の形態情報は、光磁気ディスクなどの記
憶媒体に保存される。データ解析部９は、処理の続行の
指示によって、この光磁気ディスク１３から三次元の形
態情報を読みだす（ステップＶ８）。

【００３９】この三次元の形態情報に基づきいた球モデ
ルを設定し、ステップＳ４で得られた電流源の磁界を計
測した計測磁界データとこの球モデルとに基づき、最小
ノルム法などの電流源推定方法を使い電流源の推定を行
う（ステップＶ９）。

【００４０】次に、光磁気ディスク１３から三次元の形
態情報の一部である断層画像を呼び出す。データ解析部
９は、この断層画像上にステップＶ９で推定された電流
源を重ねて、カラーモニタ１０に表示させたり、カラー
プリンタ１１に印字させたりする（ステップＶ１０）。

【００４１】

【発明の効果】以上の説明から明らかなように、請求項
１の発明によれば、複数個の発信器を被検体の関心部位
の表面に配置しているので、この発信器からの磁界デー
タを磁気センサで計測することで、発信器の位置を特定
することができる。この発信器の位置に基づいて、球モ
デルを設定する。その後、関心部位の電流源の磁界計測
データに基づき、電流源を推定し、その結果を表示装置
に表示しているので、ＭＲＩ装置で関心部位の撮像をし
なくてもこの磁界計測後、すぐに計測された計測磁界デ
ータが有効か否かを判断することができ、また、複数個
の推定された電流源間の相対的な距離が知りたい場合な
どのＭＲＩ装置で関心部位の撮像をする必要がない場合
の診断期間を短くすることができる。

【００４２】請求項２の発明によれば、複数個の発信器
を被検体の関心部位の表面に配置しているので、この発
信器からの磁界データを磁気センサで計測することで、
発信器の位置を特定することができる。この発信器の位
置に基づいて、球モデルを設定し、その球モデルに対し
て、関心部位の電流源の磁界計測データに基づき、電流
源を推定し、推定された電流源の位置と向きのデータか
ら推測した推測磁界データと前記計測磁界データとの一
致する度合いを示す値と、予め定められた値とを比較す
ることで、有効なデータについて、被検体の関心部位に
対する電流源の推定を行なうので、被検体の関心部位の
電流源を正確に推定することができる。さらに、電流源
の磁界計測から電流源の推定までを効率よくすることが
でき、短い診断期間で診断をすることができる。

【図面の簡単な説明】

【図１】この発明に係る生体活動電流源推定装置の概略
構成を示したブロック図である。

【図２】発信コイルの構成を示す図である。

【図３】実施例１の生体活動電流源推定装置で行なわれ
る処理手順を示すフローチャートである。

【図４】発信コイル位置を特定するステップＳ２で行な
われる処理手順を示すフローチャートである。

【図５】実施例２の生体活動電流源推定装置で行なわれ
る処理手順を示すフローチャートである。

【符号の説明】

１ … マルチチャンネルＳＱＵＩＤセンサ ７ … 電流供給制御部 ９ … データ解析部 ２３… 発信コイル Ｍ … 被検体

Claims (2)

【特許請求の範囲】
1. 【請求項１】 被検体の関心部位の生体活動電流源から
   生じる微小磁界を、前記被検体の関心部位に近接配備さ
   れた磁気センサで計測するとともに、画像撮像装置で得
   られた三次元の形態情報に基づき、球モデルを設定し、
   この球モデルに対して、前記磁気センサで得られた計測
   磁界データに基づいて、前記生体活動電流源の位置、向
   き、大きさを推定する生体活動電流源推定装置におい
   て、（ａ）前記被検体の関心部位の表面に装着された複
   数個の発信器と、（ｂ）前記各発信器から発生される磁
   界データに基づき、前記磁気センサを基準とした座標系
   における前記各発信器の位置を得る発信器位置特定手段
   と、（ｃ）前記発信器位置特定手段で得られた前記磁気
   センサを基準とした座標系における全発信器の位置に基
   づき、球モデルを設定する球モデル設定手段と、（ｄ）
   前記発信器の位置に基づく球モデルに対して、前記計測
   磁界データに基づき電流源の位置、向き、大きさを推定
   する電流源推定手段と、（ｅ）前記推定された電流源の
   結果を表示する表示装置とを備えたことを特徴とする生
   体活動電流源推定装置。
2. 【請求項２】 被検体の関心部位の生体活動電流源から
   生じる微小磁界を、前記被検体の関心部位に近接配備さ
   れた磁気センサで計測するとともに、画像撮像装置で得
   られた三次元の形態情報に基づき、球モデルを設定し、
   この球モデルに対して、前記磁気センサで得られた計測
   磁界データに基づいて、前記生体活動電流源の位置、向
   き、大きさを推定する生体活動電流源推定装置におい
   て、（ａ）前記被検体の関心部位の表面に装着された複
   数個の発信器と、（ｂ）前記各発信器から発生される磁
   界データに基づき、前記磁気センサを基準とした座標系
   における前記各発信器の位置を得る発信器位置特定手段
   と、（ｃ）前記発信器位置特定手段で得られた前記磁気
   センサを基準とした座標系における全発信器の位置に基
   づき、球モデルを設定する球モデル設定手段と、（ｄ）
   前記発信器の位置に基づく球モデルに対して、前記計測
   磁界データに基づき電流源の位置、向き、大きさを推定
   する電流源推定手段と、（ｆ）前記推定された電流源の
   位置、向きのデータから得られる磁界を推測し、この推
   測磁界データと前記計測磁界データとの一致する度合い
   を示す値を算出する一致度算出手段と、（ｇ）前記一致
   度算出手段で算出された値と、予め定められた値とを比
   較し、前記一致する度合いを示す値が前記予め定められ
   た値よりも大きくなるまで、上述した被検体の関心部位
   の磁界計測、この磁界データに基づく電流源の推定、お
   よび一致する度合いを示す値の算出を繰り返し実行させ
   る判断手段とを備えたことを特徴とする生体活動電流源
   推定装置。