JP3387236B2 - 生体磁気計測装置 - Google Patents

Description

【発明の詳細な説明】 【０００１】 【産業上の利用分野】この発明は、生体活動電流源の位
置、向き、大きさを推定するための生体磁気計測装置に
係り、特に、推定した生体活動電流源の位置、向き、大
きさの良否の評価を正確に行うための技術に関する。 【０００２】 【従来の技術】生体に刺激を与えると、細胞膜を挟んで
形成されている分極が壊れて、生体活動電流が流れる。
この生体活動電流は、脳や心臓などにおいて現れ、脳
波、心電図として記録される。また、生体活動電流によ
って生じる磁界は、脳磁図、心磁図として記録される。 【０００３】近年、生体内の微小な磁界を計測する装置
として、ＳＱＵＩＤ（Super-conducting Quantum Inter
face Device ：超電導量子干渉計）を用いたセンサが開
発されている。このＳＱＵＩＤセンサは、デュワーと呼
ばれる容器内に、磁気センサが液体窒素などの冷媒に浸
漬して収納されている。また、このＳＱＵＩＤセンサに
は、シングルチャネルＳＱＵＩＤセンサとマルチチャネ
ルＳＱＵＩＤセンサとがあり、シングルチャネルＳＱＵ
ＩＤセンサは、上記磁気センサを１個備え、マルチチャ
ンルＳＱＵＩＤセンサは、磁気センサを複数個備えてい
る。 【０００４】例えば、脳内の生体活動電流による微小磁
界の計測を、シングルチャネルＳＱＵＩＤセンサで行う
際には、被検体の頭部の近辺の複数の計測位置に、係る
ＳＱＵＩＤセンサを移動して、各計測位置における磁界
データを無浸襲に検出して収集し、複数個の計測データ
を得る。また、同様の計測を、マルチチャネルＳＱＵＩ
Ｄセンサで行う際には、係るＳＱＵＩＤセンサを被検体
の頭部に近接配置し（このとき、デュワー内の各磁気セ
ンサは各計測位置に位置される）、各磁気センサで磁界
データを無浸襲に検出して収集し、複数個の計測データ
を得る。 【０００５】そして、得られた複数個の計測データに基
づき、病巣に関連した生体活動電流源の位置、向き、大
きさを推定し、推定した電流源をＸ線ＣＴ装置や磁気共
鳴断層撮影装置（ＭＲＩ装置）等で得られた断層画像上
に重ねて表示させるなどして、患部などの物理的位置の
特定などに用いている。 【０００６】ところで、従来の電流源の推定は、最小自
乗法や最小ノルム法、格子点移動最小ノルム法等で行わ
れている。 【０００７】最小自乗法は、電流源を仮想的に設定し、
各計測位置における磁気センサが、上記仮想した電流源
により生じる微小磁界を計測したと仮定したときの仮想
的な磁場データを計算によって求め、この計算で求めた
データと、実際に計測したデータとの自乗誤差を求め、
この自乗誤差が最小になるように、仮想した電流源を移
動させ、真の電流源に近づける方法である。 【０００８】また、最小ノルム法は、検査対象領域内に
格子点を設定し、格子点上の電流源と各計測データとの
間の関係を示す行列を求め、この行列の一般逆行列と、
計測データとの積をとって電流源（格子点上の分布）を
求める方法である（例えば、W.H.Kullmann，K.D.Jandt,
K.Rehm，H.A.Schlitt, W.J.Dallas and W.E.Smith,Adva
nces in Biomagnetism pp.571-574, Plenum Pless, New
York, 1989 や、次の特願平5-160450号、特願平5-1604
51号等の従来技術等を参照）。 【０００９】さらに、格子点移動最小ノルム法は、本願
出願人によって、特願平5-160450号や特願平5-160451号
等で提案している方法であり、上記最小ノルム法の格子
点の数を変えずに、上記最小ノルム法で得られた分布の
うち、大きい電流源の近くに各格子点を移動し、その格
子点の配置における電流源の分布を最小ノルム法で求
め、以下、同様の処理を、最小の格子点間隔が予め設定
された収束値以下になるまで繰り返し、電流源を推定す
る方法である。 【００１０】 【発明が解決しようとする課題】しかしながら、このよ
うな構成を有する従来例の場合には、次のような問題が
ある。すなわち、最小自乗法で電流源を推定した場合、
推定結果は自乗誤差によりその良否が評価できる。すな
わち、自乗誤差が「０」に近いほど推定された電流源は
真の電流源に近いと評価できる。しかしながら、計測デ
ータにノイズが含まれている場合には、自乗誤差の大小
と評価の良否とは必ずしも一致せず、正確な評価が行え
ないという問題がある。 【００１１】また、最小ノルム法では、推定結果の良否
の判断を適正に行うことができないという問題がある。 【００１２】さらに、格子点移動最小ノルム法でも、各
格子点の配置における電流源の推定は、最小ノルム法で
行うので、各格子点の配置における推定結果の良否の判
断を適正に行えず、また、最終の電流源を推定するため
に、最小の格子点間隔が予め設定された収束値以下にな
ったか否かを収束判定条件としているため、この設定さ
れた収束判定値の値によって推定結果が異なり、推定法
に一般性をもたせることができないという問題がある。 【００１３】この発明は、このような事情に鑑みてなさ
れたものであって、各推定方法で推定した結果の良否を
適正、正確に評価することができる生体磁気計測装置を
提供することを目的とする。 【００１４】 【課題を解決するための手段】この発明は、このような
目的を達成するために、次のような構成をとる。すなわ
ち、この発明は、被検体の診断対象領域内の生体活動電
流源により生じる微小磁界を、磁気センサで前記診断対
象領域の近傍の複数の計測位置から計測して複数個の計
測データを収集し、それら複数個の計測データに基づ
き、前記生体活動電流源の位置、向き、大きさを推定す
るための生体磁気計測装置において、（ａ）前記複数個
の計測データのうちの一部の計測データ群（以下、これ
らの計測データを「第１のグループの計測データ」とい
う）を用いて生体活動電流源の位置、向き、大きさを推
定する推定手段と、（ｂ）前記第１のグループの計測デ
ータに含まれない計測データを含む計測データ群（以
下、これらの計測データを「第２のグループの計測デー
タという）の計測データを計測したときの各計測位置に
おける磁気センサが、前記推定手段で推定された生体活
動電流源により生じる微小磁界を計測したと仮定したと
きの仮想的な磁場データを計算によって求め、この計算
で求めたデータと、前記第２のグループの計測データと
の自乗誤差を求めて、前記推定手段で推定した生体活動
電流源の位置、向き、大きさの良否を判断するための評
価値を得る評価値算出手段とを備えたものである。 【００１５】 【作用】この発明の作用は次のとおりである。推定手段
は、複数個の計測データのうちの一部の計測データ群
（第１のグループの計測データ）を用いて生体活動電流
源の位置、向き、大きさを、最小自乗法や最小ノルム法
等で推定する。そして、評価値算出手段は、第１のグル
ープの計測データに含まれない計測データを含む計測デ
ータ群（第２のグループの計測データ）の計測データを
計測したときの各計測位置における磁気センサが、推定
手段で推定された生体活動電流源により生じる微小磁界
を計測したと仮定したときの仮想的な磁場データを計算
によって求め、この計算で求めたデータと、第２のグル
ープの実際の計測データとの自乗誤差を求めて、推定手
段で推定した生体活動電流源の位置、向き、大きさの良
否を判断するための評価値を求める。 【００１６】なお、第１のグループの計測データと第２
のグループの計測データの採り方としては、例えば、全
計測データを２分し、一方の計測データ群を第１のグル
ープの計測データ、他方の計測データ群を第２のグルー
プの計測データに設定したり、各グループの計測データ
に含まれる計測データの一部が互いに重なるように設定
するなど、種々の採り方が選択できる。しかし、第２の
グループの計測データが、完全に第１のグループの計測
データに含まれる採り方はできない。 【００１７】 【実施例】以下、図面を参照してこの発明の一実施例を
説明する。図１は、この発明の実施例に係る生体磁気計
測装置の概略構成を示した図であり、図２は実施例で用
いられるマルチチャネルＳＱＵＩＤセンサの概略構成を
示す図である。 【００１８】図中、符号２は磁気シールドルームであ
り、この磁気シールドルーム２内に被検体Ｍが仰臥され
るベッド３と、被検体Ｍの例えば頭部に近接配備され、
被検体の脳内に生じた生体活動電流源による微小磁界を
無浸襲に計測するためのマルチチャネルＳＱＵＩＤセン
サ１とが設けられている。マルチチャネルＳＱＵＩＤセ
ンサ１は、図２に示すように、デュアー１ａ内に複数個
の磁気センサＳ₁ 〜Ｓ_mが冷媒に浸漬して収納されてい
る。デュアー１ａが被検体の頭部の近傍に配置されたと
きの各磁気センサＳ₁ 〜Ｓ_m の配置位置を計測位置とし
て、各磁気センサＳ₁ 〜Ｓ_m は被検体Ｍの脳内の生体活
動電流源により生じる磁界データをそれぞれ計測する。 【００１９】マルチチャネルＳＱＵＩＤセンサ１内の各
磁気センサＳ₁ 〜Ｓ_m で計測された複数個の計測データ
はデータ変換ユニット４に与えられてディジタルデータ
に変換された後、データ収集ユニット５に集められる。
刺激装置６は、被検体Ｍに電気的刺激（あるいは音、光
刺激など）を与えるためのものである。ポジショニング
ユニット７は、マルチチャネルＳＱＵＩＤセンサ１を基
準とした３次元座標系に対する被検体Ｍの頭部の位置関
係を把握するための装置である。例えば、被検体Ｍの頭
部の複数個所（例えば、両眉の中間と両耳の付け根の３
個所）に小コイルを取付け、これら小コイルにポジショ
ニングユニット７から給電する。そして、各小コイルか
ら発生した磁界をマルチチャネルＳＱＵＩＤセンサ１で
検出することにより、マルチチャネルＳＱＵＩＤセンサ
１に対する被検体Ｍの頭部の位置関係を把握する。な
お、ＳＱＵＩＤセンサ１に対する被検体Ｍの頭部の位置
関係を把握するための手法は、これ以外に、デュアー１
ａに投光器を取付けて光ビームで被検体Ｍの頭部を照射
して両者の位置関係を把握するものや、あるいは、特開
平5-237065号公報や特開平6-788925号公報などに開示さ
れた種々の手法が用いられる。 【００２０】データ解析ユニット８は、推定手段および
評価値算出手段に相当するものであり、データ収集ユニ
ット５に集められた計測データに基づいて、被検体Ｍの
検査対象領域内（脳内）の電流源を推定するとともに、
その推定結果の評価値を求めるためのものである。デー
タ解析ユニット８に関連して設けられた光磁気ディスク
９には、例えばＸ線ＣＴ装置やＭＲＩ装置で得られた断
層画像が記録されており、データ解析ユニット８で推定
された電流源が、これら断層画像上に重ね合わされてカ
ラーモニタ１０に表示されたり、あるいはカラープリン
タ１１に印字出力され、また、評価値は、例えば、操作
モニタ１２などに表示されて操作者に通知される。な
お、Ｘ線ＣＴ装置やＭＲＩ装置で得られた断層画像は、
図１に示した通信回線１３を介してデータ解析ユニット
８に直接伝送するように構成してもよい。 【００２１】また、図１の装置では、マルチチャネルＳ
ＱＵＩＤセンサ１を用いて電流源から生じる磁界データ
を計測するように構成しているが、例えば、シングルチ
ャネルＳＱＵＩＤセンサを用いて電流源から生じる磁界
データを計測する場合には、図３に示すように、デュア
ー１ａと被検体Ｍとの位置関係を変位して、各位置にお
ける磁気センサＳ_S の配置位置を計測位置として、磁気
センサＳ_S は複数の計測位置において、電流源から生じ
る磁界データを計測し、複数個の計測データをデータ収
集ユニット５に集めるように構成される。 【００２２】次に、上記データ解析ユニット８で実行さ
れる処理について、図４に示すフローチャートを参照し
て説明する。図４のフローチャートは、最小自乗法また
は最小ノルム法で電流源を推定する場合の手順を示す。 【００２３】まず、データ収集ユニット５に集められた
複数個の計測データをグループ分けする（ステップＳ
１）。ここでは、図５に示すように、全計測データＤ₁
〜Ｄ_m（各磁気センサＳ₁ 〜Ｓ_m で計測されたデータ）
を２分し、一方の計測データ群（Ｄ₁ 〜Ｄ_i ：但し、ｉ
は２以上、（ｍ−２）以下の自然数）を第１のグループ
ＧＲ１の計測データとし、残りの計測データ群（Ｄ_i+1
〜Ｄ_m ）を第２のグループＧＲ２の計測データとしてグ
ループ分けしている。 【００２４】次に、第１のグループＧＲ１の計測データ
（Ｄ₁ 〜Ｄ_i ）を用いて、電流源を推定する（ステップ
Ｓ２）。この推定方法は、最小自乗法や最小ノルム法で
あり、これら方法は、全計測データを用いて電流源を推
定（従来例）するか、一部の計測データを用いて電流源
を推定（この実施例）するかが異なるだけで、その基本
的な処理内容は従来例と同様であるので、ここでの詳述
は省略する。 【００２５】そして、第２のグループＧＲ２の計測デー
タ（Ｄ_i+1 〜Ｄ_m ）を計測したときの各計測位置におけ
る磁気センサ（この実施例では各磁気センサＳ_i+1 〜Ｓ
_m ）が、ステップＳ２で推定された電流源により生じる
微小磁界を計測したと仮定したときの仮想的な磁場デー
タを計算によって求める（ステップＳ３）。この計算方
法は、例えば、従来例の最小自乗法において、各磁気セ
ンサが、仮想的に設定した電流源により生じる微小磁界
を計測したと仮定したときの仮想的な磁場データを求め
る際の計算と同様であるので、ここでの詳述は省略す
る。 【００２６】次に、ステップＳ３で計算したデータ（Ｔ
_i+1 〜Ｔ_m とする：但し、添字は、磁気センサＳ_i+1 〜
Ｓ_m に対応する）と、第２のグループＧＲ２の実際の計
測データ（Ｄ_i+1 〜Ｄ_m ）との自乗誤差を求め、評価値
ｆを求める（ステップＳ４）。以下に、その計算式を示
す。 【００２７】 【数１】 【００２８】そして、ステップＳ４で求めた評価値ｆ
（ステップＳ２で求めた電流源の推定結果に対する評価
値）を、操作モニタ１２に表示出力する（ステップＳ
５）。この評価値ｆは、値が小さい方がステップＳ２で
求めた電流源の推定結果が良いことを示す。 【００２９】なお、所定の判断値を予め決めておき、上
記評価値ｆとその判断値とを比較して、評価値ｆが判断
値よりも大きいときには、ステップＳ２の推定条件（例
えば、最小自乗法であれば、仮想的に設定する電流源の
初期値や電流源の移動のやり方等、また、最小ノルム法
であれば、格子点の数や配置等）を変えて、ステップＳ
２からもう一度処理をやり直し、判断値よりも小さい評
価値ｆが得られるまで、その処理を繰り返して、真の電
流源により近い電流源を自動的に求められるように構成
してもよい。 【００３０】このように、この実施例では、全計測デー
タ（Ｄ₁ 〜Ｄ_m ）の一部で構成される第１のグループＧ
Ｒ１の計測データ（Ｄ₁ 〜Ｄ_i ）を用いて電流源を推定
し、その推定に用いた計測データ（Ｄ₁ 〜Ｄ_i ）以外の
計測データ（Ｄ_i+1 〜Ｄ_m ）で構成される第２のグルー
プＧＲ２の計測データを用いて、第１のグループＧＲ１
の計測データを用いて推定した結果（電流源）を評価す
るための評価値ｆを求めるので、第１のグループＧＲ１
の計測データを用いた推定方法が最小ノルム法の場合で
あっても、その推定結果を適正に評価することが可能と
なった。 【００３１】また、評価値ｆは、推定に用いた計測デー
タ以外の計測データを用いて求められるので、従来の最
小自乗法において、全計測データを用いて得られる自乗
誤差よりも客観性が高まり、従って、従来の最小自乗法
に比べて推定の良否を正確に評価することができる。 【００３２】ところで、この上述の実施例のステップＳ
１におけるグループ分けでは、全計測データを２分して
行ったが、この発明における第１のグループの計測デー
タと第２のグループの計測データの採り方は、上記図５
の採り方に限定されるものではない。すなわち、第２の
グループＧＲ２の計測データは、第１のグループＧＲ１
の計測データに含まれない計測データを含めた計測デー
タ群で構成されていればよく、例えば、図６（ａ）に示
すように、双方のグループＧＲ１、ＧＲ２の計測データ
の一部が重なるように各グループＧＲ１、ＧＲ２の計測
データを決めてもよいし、図６（ｂ）に示すように、第
１のグループＧＲ１の計測データの全てが第２のグルー
プＧＲ２の計測データに完全に含まれるように各グルー
プＧＲ１、ＧＲ２の計測データを決めてもよい。しか
し、図６（ｃ）に示すように、第２のグループＧＲ２の
計測データの全てが第１のグループＧＲ１の計測データ
に完全に含まれるように各グループＧＲ１、ＧＲ２の計
測データを決めることはできない。これは、図５、図６
（ａ）、（ｂ）では、推定に用いた計測データ以外の計
測データを含めて表価値ｆを求めており、推定に用いた
計測データ以外の計測データが反映された評価値ｆが得
えれ、推定結果を客観的に評価することができるが、図
６（ｃ）では、推定に用いた計測データ以外の計測デー
タが評価値ｆに反映されず、推定結果を客観的に評価す
ることができないからである。なお、推定結果を、より
客観的に評価（すなわち、より正確に評価）するために
は、図５に示すように各グループＧＲ１、ＧＲ２の計測
データが完全に分かれていることが好ましい。 【００３３】次に、データ解析ユニット８により、格子
点移動最小ノルム法で電流源の推定等を行う場合の処理
について、図７に示すフローチャートを参照して説明す
る。 【００３４】まず、全計測データをグループ分けする
（ステップＳ１１）。これは、上記ステップＳ１と同じ
であるので、重複する説明は省略する。 【００３５】次に、診断対象領域に格子点を設定し（ス
テップＳ１２）、最小ノルム解（最小ノルム法による電
流源の推定）を求める（ステップＳ１３）。これら処理
は、従来例と同じであるので、その詳述は省略する。 【００３６】次に、ステップＳ１４、Ｓ１５で評価値ｆ
を求める。これは、上記ステップＳ３、Ｓ４と同様であ
る。 【００３７】そして、予め設定しておいた終了判定値
と、求めた評価値ｆとを比較し、評価値ｆが終了判定値
よりも小さければ処理を終了し、逆に大きければステッ
プＳ１７、Ｓ１８を実行して、再び、ステップＳ１３に
戻る（ステップＳ１６）。 【００３８】ステップＳ１７では、評価値ｆが最小とな
る最適な移動パラメータを求める。この移動パラメータ
は、特願平5-160451号等においてαやβ、γで説明され
ている、格子点の移動を制御するためのパラメータであ
り、上記出願の明細書にいおても説明しているように、
従来はこの移動パラメータは経験的に決めている。しか
し、この実施例では、評価値ｆを用いて、最適な移動パ
ラメータを自動的に決めることができる。 【００３９】具体的には、移動パラメータに具体的な数
値を代入して、格子点を移動し、移動後の格子点におけ
る評価値ｆを計算し、評価値ｆが最小となったときの移
動パラメータの値を最適な移動パラメータとして選択す
る。 【００４０】そして、ステップＳ１８では、上記移動パ
ラメータを用いて格子点を移動し、ステップＳ１３に戻
る。 【００４１】上述したように、この実施例によれば、格
子点移動最小ノルム法で電流源を推定する場合において
も、推定結果の評価を適正に行うことができるし、終了
判定にも一般性を持たせることができる。また、上記し
たように、最適な移動パラメータを自動的に決めること
ができる。 【００４２】なお、最適な移動パラメータを決める場合
の評価値ｆは、以下のようにして求めてもよい。 【００４３】すなわち、第２のグループＧＲ２の計測デ
ータを用いて求めた評価値ｆをｆ₂とし、これと逆に、
第２のグループＧＲ２の計測データを用いて電流源を推
定し、第１のグループＧＲ１の計測データを用いて：第
２のグループＧＲ２の計測データを用いて推定した電流
源の良否を評価するための評価値ｆを求めて、これをｆ
₁ とする。そして、（ｆ₁ ＋ｆ₂ ）を評価値ｆとして、
最適な移動パラメータを選択する。このように構成する
ことにより、各グループＧＲ１、ＧＲ２の計測データを
用いた評価値を加味させることができ、客観性が高くな
る。 【００４４】また、図８に示すように、全計測データを
３組（ｇｒ１、ｇｒ２、ｇｒ３）に分割し、まず、第１
組ｇｒ１の計測データを上記第１のグループＧＲ１の計
測データに、第２組ｇｒ２、第３組ｇｒ３を合わせた計
測データを上記第２のグループＧＲ２の計測データに割
り当てて、第１のグループＧＲ１の計測データを用いて
推定した推定結果に対する評価値ｆ₂₃（第２組ｇｒ２、
第３組ｇｒ３の計測データを用いて算出）を求める。次
に、第２組ｇｒ２の計測データを上記第１のグループＧ
Ｒ１の計測データに、第１組ｇｒ１、第３組ｇｒ３を合
わせた計測データを上記第２のグループＧＲ２の計測デ
ータに割り当てて、上記と同様の手順で求めた評価値ｆ
₁₃、および、第３組ｇｒ３の計測データを上記第１のグ
ループＧＲ１の計測データに、第１組ｇｒ１、第２組ｇ
ｒ２を合わせた計測データを上記第２のグループＧＲ２
の計測データに割り当てて求めた評価値ｆ₁₂をそれぞれ
求める。そして、（ｆ₂₃＋ｆ₁₃＋ｆ₁₂）を評価値ｆとし
て、最適な移動パラメータを選択するようにしてもよ
い。 【００４５】さらに、全計測データをｎ組（ｎは４以上
の自然数）に分割し、上記３組に分割した場合と同様に
方法で、最適な移動パラメータを選択するための評価値
ｆを求めてもよい。 【００４６】 【発明の効果】以上の説明から明らかなように、この発
明によれば、全計測データの一部で構成される第１のグ
ループの計測データを用いて電流源を推定し、その推定
に用いた計測データ以外の計測データを含む第２のグル
ープの計測データを用いて、第１のグループの計測デー
タを用いて推定した結果（電流源）を評価するための評
価値を求めるように構成したので、第１のグループの計
測データを用いた推定方法が最小ノルム法や格子点移動
最小ノルム法であっても、その推定結果を適正に評価す
ることが可能となった。また、この評価値を格子点移動
最小ノルム法の収束判定条件にも利用することができ、
格子点移動最小ノルム法に一般性を持たせることが可能
となった。 【００４７】さらに、評価値は、推定に用いた計測デー
タ以外のデータを含めて求めるので、従来の最小自乗法
において、全計測データを用いて得られる自乗誤差より
も客観性が高まり、従って、従来の最小自乗法に比べて
推定結果の良否を正確に評価することができる。

【図面の簡単な説明】 【図１】この発明の実施例に係る生体磁気計測装置の概
略構成を示した図である。 【図２】実施例で用いられるマルチチャネルＳＱＵＩＤ
センサの概略構成を示す図である。 【図３】シングルチャネルＳＱＵＩＤセンサで電流源か
らの磁界を計測する場合の構成を説明するための図であ
る。 【図４】実施例のデータ解析ユニットが、最小自乗法、
最小ノルム法で電流源を推定する場合の処理手順を示す
フローチャートである。 【図５】計測データのグループ分けを示す図である。 【図６】計測データのグループ分けの変形例等を示す図
である。 【図７】実施例のデータ解析ユニットが、格子点移動最
小ノルム法で電流源を推定する場合の処理手順を示すフ
ローチャートである。 【図８】評価値の算出の変形例を説明するための図であ
る。 【符号の説明】 １ … マルチチャネルＳＱＵＩＤセンサ ８ … データ解析ユニット Ｍ … 被検体 Ｓ₁ 〜Ｓ_m … 磁気センサ Ｄ₁ 〜Ｄ_m … 計測データ ＧＲ１ … 第１のグループ ＧＲ２ … 第２のグループ

Claims (1)

1. (57)【特許請求の範囲】 【請求項１】 被検体の診断対象領域内の生体活動電流
   源により生じる微小磁界を、磁気センサで前記診断対象
   領域の近傍の複数の計測位置から計測して複数個の計測
   データを収集し、それら複数個の計測データに基づき、
   前記生体活動電流源の位置、向き、大きさを推定するた
   めの生体磁気計測装置において、（ａ）前記複数個の計
   測データのうちの一部の計測データ群（以下、これらの
   計測データを「第１のグループの計測データ」という）
   を用いて生体活動電流源の位置、向き、大きさを推定す
   る推定手段と、（ｂ）前記第１のグループの計測データ
   に含まれない計測データを含む計測データ群（以下、こ
   れらの計測データを「第２のグループの計測データとい
   う）の計測データを計測したときの各計測位置における
   磁気センサが、前記推定手段で推定された生体活動電流
   源により生じる微小磁界を計測したと仮定したときの仮
   想的な磁場データを計算によって求め、この計算で求め
   たデータと、前記第２のグループの計測データとの自乗
   誤差を求めて、前記推定手段で推定した生体活動電流源
   の位置、向き、大きさの良否を判断するための評価値を
   得る評価値算出手段とを備えたことを特徴とする生体磁
   気計測装置。