(57)【特許請求の範囲】
【請求項１】
被検体の関心部位に複数個の指標点を決め、この指標点にそれぞれ発振器を装着して、前記発振器から発生する磁界を、前記被検体の関心部位に近接配備させた磁束計で計測することで計測磁界情報を得て、この計測磁界情報に基づき前記磁束計を基準とする三次元の座標系における前記発振器の位置を特定する発振器位置特定手段と、前記被検体の関心部位の生体活動電流源から発生する生体磁気を計測し、この計測で得られた生体磁気情報に基づき生体活動電流源を推定する電流源推定手段とを備え、画像撮像装置で被検体の関心部位を撮像することで得た複数の断層画像上の指標点と、前記発振器位置特定手段で特定された発振器の位置とをそれぞれ重ねることで、前記断層画像に前記電流源推定手段で推定された生体活動電流源の位置を重ね、被検体の関心部位に対する生体活動電流源の位置を把握する生体磁気計測装置において、前記発振器位置特定手段が、（ａ）前記被検体の関心部位に装着する複数個の発振器と、（ｂ）前記各発振器に供給される電流を検出し、この電流の電流強度と電流周波数とを求める電流情報検出手段と、（ｃ）前記電流周波数に基づき前記計測磁界情報に対して周波数解析を行うことで解析磁界情報を得る計測情報解析手段と、（ｄ）前記磁束計を基準とする三次元の座標系の仮想の位置に仮想発振器を設定し、前記仮想発振器に前記電流強度と同じ大きさの電流強度を与えた場合に前記磁束計で得られる仮想磁界情報を理論的に求め、前記解析磁界情報と前記仮想磁界情報との最小自乗誤差を算出し、この最小自乗誤差が最小となる仮想発振器の位置を求め、この位置を発振器の位置とすることで、前記磁束計を基準とする三次元の座標系における前記複数個の発振器の位置をそれぞれ特定する位置特定手段とを備えたことを特徴とする生体磁気計測装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
この発明は、被検体の関心部位の生体活動電流源を推定する生体磁気計測装置に係り、特に、推定された生体活動電流源を被検体の関心部位の撮像画像に正確に重ねることができる技術に関する。
【０００２】
【従来の技術】
生体に刺激を与えると、細胞膜を挟んで形成されている分極が壊れて生体活動電流が流れる。この生体活動電流は、脳や心臓において現れ、脳波，心電図として記録される。また、生体活動電流によって生じる磁界は、脳磁図，心磁図として記録される。
【０００３】
近年、生体内の微小な磁界を計測する装置として、ＳＱＵＩＤ（Superconduc-ting Quantum Interference Device：超電導量子干渉計）を用いたマルチチャンネルＳＱＵＩＤセンサが開発されている。このマルチチャンネルＳＱＵＩＤセンサは、デュアーと呼ばれる容器内に多数の磁気センサを液体窒素などの冷媒に浸漬して収納している。
【０００４】
このマルチチャンネルＳＱＵＩＤセンサを被検体の関心部位である例えば頭部の外側に置き、脳内に生じた生体活動電流源である電流双極子（以下、「電流源」と呼ぶ）による微小磁界をそのマルチチャンネルＳＱＵＩＤセンサで無侵襲に計測することができる。計測された生体磁界データから被検体の関心部位の電流源の位置, 向き, 大きさが推定されている。この電流源の推定方法の一つとして、最小ノルム法を用いた手法がある（例えば、W.H.Kullmann, K.D.Jandt, K.Rehm, H.A.Schlitt, W.J.Dallas and W.E.Smith, Advances in Biomagnetism, pp.571-574, Plenum Pless, New York, 1989) 。
【０００５】
このような方法で推定された電流源は、ＭＲＩ装置等で得られた被検体の関心部位の断層画像上に表示させることで、被検体の関心部位に対する電流源の位置を把握することができる。従って、マルチチャンネルＳＱＵＩＤセンサを基準とした三次元の座標系における被検体の関心部位の位置を正確に把握することが重要である。
【０００６】
そこで、被検体の関心部位に複数個の発振コイルを取り付る。この発振コイルからの磁界をマルチチャンネルＳＱＵＩＤセンサで計測して計測磁界データを得る。この計測磁界データに基づいてマルチチャンネルＳＱＵＩＤセンサを基準とした三次元の座標系における各発振コイルの位置と、ＭＲＩ装置等で撮像して得られた発振コイルに位置である指標点とを重ねることで被検体の関心部位の電流源の位置などを把握している。ＭＲＩ装置で撮影を行う場合は、発振コイルをマーカに付けかえてから撮像する。
【０００７】
従来、被検体の関心部位の複数個の発振コイルからの計測磁界データに基づいてマルチチャンネルＳＱＵＩＤセンサを基準とした三次元の座標系における各発振コイルの位置を特定する方法として、以下のような方法が用いられている。
【０００８】
まず、被検体の関心部位上に複数の指標点を決める。被検体の関心部位である例えば頭部の診断では、両眉の中心点と、両耳の付け根の３箇所を指標点としている。この各指標点に発振コイルを取り付ける。この発振コイルに予め設定された値の電流強度と電流周波数の電流（以下、「既知の電流」と呼ぶ）を電流供給装置から送出する。この各発振コイルに個別に既知の電流を順次送出することで、発振コイルから磁界を発生させる。この磁界をマルチチャンネルＳＱＵＩＤセンサで計測することで、各発振コイルの計測磁界データを得ることができる。
【０００９】
次に、マルチチャンネルＳＱＵＩＤセンサを基準とする座標系に、既知の電流値に基づく磁界を発生させる仮想発振コイルの位置を仮想的に設定する。この設定された仮想位置の仮想発振コイルから発生する磁界をマルチチャンネルＳＱＵＩＤセンサで計測した場合の仮想磁界データを求める。
【００１０】
この仮想磁界データと実際の計測磁界データとの最小自乗誤差を求める。この最小自乗誤差の値が最小になるように、仮想発振コイルの仮想位置を動かして仮想位置を特定する。この最小自乗誤差の値が最小となる仮想発振コイルの仮想位置が、実際の発振コイルの位置である。この発振コイルの位置を特定する処理を各発振コイルについて行うことで、各発振コイルのマルチチャンネルＳＱＵＩＤセンサを基準とする三次元の座標系の位置をそれぞれ特定することができる。
【００１１】
【発明が解決しようとする課題】
しかしながら、このような構成を有する従来例の場合には、次のような問題がある。
すなわち、電流を供給する電流供給装置の性能や周囲の環境温度などの種々の原因によって、発振コイルに供給される電流は、予め設定された値の電流を供給できない場合が多い。しかし、予め設定された値の電流が発振コイルに供給されているものとして、この設定された値の電流に基づき発振コイルの位置を特定しているので、特定される発振コイルの位置が実際の発振コイルの位置と大きくずれるという問題がある。
【００１２】
この発明は、このような事情に鑑みてなされたものであって、被検体の関心部位の装着した発振コイルの位置を正確に特定することができる生体磁気計測装置を提供することを目的とする。
【００１３】
【課題を解決するための手段】
この発明は、このような目的を達成するために、次のような構成をとる。
すなわち、この発明は、被検体の関心部位に複数個の指標点を決め、この指標点にそれぞれ発振器を装着して、前記発振器から発生する磁界を、前記被検体の関心部位に近接配備させた磁束計で計測することで計測磁界情報を得て、この計測磁界情報に基づき前記磁束計を基準とする三次元の座標系における前記発振器の位置を特定する発振器位置特定手段と、前記被検体の関心部位の電流源から発生する生体磁気を計測し、この計測で得られた生体磁気情報に基づき電流源を推定する電流源推定手段とを備え、画像撮像装置で被検体の関心部位を撮像することで得た複数の断層画像上の指標点と、前記発振器位置特定手段で特定された発振器の位置とをそれぞれ重ねることで、前記断層画像に前記電流源推定手段で推定された電流源の位置を重ね、被検体の関心部位に対する電流源の位置を把握する生体磁気計測装置において、前記発振器位置特定手段が、（ａ）前記被検体の関心部位に装着する複数個の発振器と、（ｂ）前記各発振器に供給される電流を検出し、この電流の電流強度と電流周波数とを求める電流情報検出手段と、（ｃ）前記電流周波数に基づき前記計測磁界情報に対して周波数解析を行うことで解析磁界情報を得る計測情報解析手段と、（ｄ）前記磁束計を基準とする三次元の座標系の仮想の位置に仮想発振器を設定し、前記仮想発振器に前記電流強度と同じ大きさの電流強度を与えた場合に前記磁束計で得られる仮想磁界情報を理論的に求め、前記解析磁界情報と前記仮想磁界情報との最小自乗誤差を算出し、この最小自乗誤差が最小となる仮想発振器の位置を求め、この位置を発振器の位置とすることで、前記磁束計を基準とする三次元の座標系における前記複数個の発振器の位置をそれぞれ特定する位置特定手段とを備えたことを特徴とするものである。
【００１４】
〔作用〕
この発明の作用は次のとおりである。
被検体の関心部位に複数個の発振器を装着する。この発振器に電流を個別に供給して発生する磁界を磁束計で計測して計測磁界情報を得る。一方、電流情報検出手段は、発振器に供給される電流を検出することで実際に供給される電流の電流強度と電流周波数とを発振器ごとに求める。計測情報解析手段は、この電流周波数に基づいて、前記計測磁界情報に周波数解析を施すことで、電流周波数に基づく解析磁界情報を求める。位置特定手段は、磁束計を基準とする三次元の座標系に仮想発振器を設定する。前記電流強度と同じ大きさの電流強度を持つ電流をこの仮想発振器に供給した場合、仮想発振器から発生する磁界を前記磁束計で計測することで得られる仮想磁界情報を理論的に求める。この仮想磁界情報と前記解析磁界情報との最小自乗誤差を計算する。この最小自乗誤差の値が最小となるような前記仮想発振器の仮想位置を求める。この仮想位置が前記発振器の位置と特定する。この位置を特定する処理を各発振器についておこなうことで、磁束計を基準とする三次元の座標系における複数個の発振器の位置を特定する。
【００１５】
【発明の実施の形態】
以下、図面を参照してこの発明の一実施例を説明する。
＜第１実施例＞
図１は実施例に係る生体磁気計測装置の概略構成を示すブロック図である。本実施例では、被検体の関心部位を頭部として電流源を推定するものとする。
【００１６】
図中、符号１はマルチチャンネルＳＱＵＩＤセンサを示す。このマルチチャンネルＳＱＵＩＤセンサ１は、被検体Ｍの頭部の外側に配備される。マルチチャンネルＳＱＵＩＤセンサ１は、デュアー１ａ内に図示しない多数の磁気センサＳ1 〜Ｓm が冷媒に侵漬して収納されて構成されている。これらマルチチャンネルＳＱＵＩＤセンサ１や被検体Ｍは、磁気シールドルーム１２内に配設される。磁気シールドルーム１２内の被検体Ｍの頭部には、両眉の中心点と、両耳の付け根の３箇所に発振コイルＣ１、Ｃ２、Ｃ３（以下、発振コイルＣ１、Ｃ２、Ｃ３をまめて「各発振コイル」と呼ぶ）がそれぞれ装着されている。なお、この実施例では３個の発振コイルを用いているが、これに限定するものではなく３個以上あればよい。
【００１７】
まず、刺激装置７は、被検体Ｍに電気的刺激（あるいは音、光刺激など）を与えるためのものである。この刺激装置７は、被検体Ｍに電気的刺激などを与え、被検体の頭部の電流源から生体磁気を発生させる。この生体磁気をマルチチャンネルＳＱＵＩＤセンサ１のデュアー１ａ内の多数の磁気センサＳ1 〜Ｓm でそれぞれ計測する。この計測された生体磁気データは、データ収集部６に集められる。電流源推定部８は、データ収集部６に集められた生体磁気データに基づいて、最小ノルム法などを用いて被検体Ｍの頭部内に発生する電流源の位置、大きさ、向き等を推定する。
【００１８】
次に、各発振コイルは、電流供給部２から順次個別に電流が供給される。この時の電流をフィードバック部３で検出し、電流情報をデータ解析部４に送る。一方、各発振コイルは供給された電流によって磁界を発生する。この磁界をマルチチャンネルＳＱＵＩＤセンサ１内の多数の磁気センサＳ1 〜Ｓm でそれぞれ計測する。この多数の磁気センサＳ1 〜Ｓm で計測された計測磁界データ群（以下、マルチチャンネルＳＱＵＩＤセンサ１で計測された計測磁界データ群」と呼ぶ）は、データ収集部６を介して計測磁界データ解析部４に送られる。データ解析部４は、電流情報に基づいて、計測磁界データ群を解析する。コイル位置特定部は、この解析データと電流情報とに基づき被検体Ｍの頭部の装着された各発振コイルの位置をそれぞれ特定するためのものである。これら各発振コイルの位置を特定する処理の詳細は後述する。
【００１９】
電流源推定部８は、ディスクドライブ９にセットされた図示しない光磁気ディスクから被検体Ｍの頭部の断層画像を呼び出す。この断層画像上の指標点と、コイル位置特定部５で特定された各発振コイルの位置とを重ねることで、被検体Ｍの頭部に発生する電流源の位置を把握することができる。この断層画像上に重ねられた電流源の位置、向き、大きさは、モニタ１０に表示されたり、プリンタ１１に印刷出力されるよになっている。なお、光磁気ディスクには、例えばＸ線ＣＴ装置やＭＲＩ装置等を用いて、被検体Ｍの関心部位である頭部を撮影することで得られる多数の断層画像が記憶されている。
【００２０】
以下、図２のフローチャート参照して被検体Ｍの頭部に装着された各発振コイルの位置の特定について説明する。
【００２１】
ステップＳ１（供給電流値を検出）
各発振コイルには、電流供給部２からフィードバック部３を介して順次個別に電流が供給される。まず、この電流供給部２から、周波数および電流強度が予め決められた値の電流を発振コイルＣ１に送出する。この時、電流供給部２と各発振コイルとの間に設置されたフィードバック部３によって、この送出さた電流を検出する。この電流の値はデジタル変換されデータ解析部４に送られる。
【００２２】
ステップＳ２（発振コイルからの磁界を計測）
発振コイルＣ１は、電流供給部２から送られた電流によって磁界を発生する。この磁界をマルチチャンネルＳＱＵＩＤセンサ１で計測する。この計測によって得られる計測磁界データ群は、データ収集部６を介してデジタル変換されデータ解析部４に送られる。
【００２３】
ステップＳ３（電流データ解析）
データ解析部４は、フィードバック部３から送られた電流の値に対して周波数解析を行うことで、実際に発振コイルＣ１に供給された電流の周波数ｆおよび電流強度ｉを知ることができる。この周波数解析とは、検出データに対してフーリエ変換を施し、発振コイルＣ１に供給される電流の値を求めることである。前記フィードバック部３とこのステップＳ３の処理を行うデータ解析部４とが、この発明の電流情報検出手段に相当する。
【００２４】
ステップＳ４（計測磁界データ解析）
データ解析部４は、発振コイルＣ１の計測磁界データ群に対して周波数解析を行い、ステップＳ３で求められた電流の周波数ｆに基づく解析磁界データ群ｅm を求める。ここで、この計測磁界データ群に周波数解析を行うことで、幾つかの周波数ごとの計測磁界データ群に分けることができる。このうち周波数ｆに対する計測磁界データ群を求めたものが、解析磁界データ群ｅm である。データ解析部４で求められた電流強度ｉと解析磁界データ群ｅm とは、コイル位置特定部５に送られる。このステップＳ４の処理を行うデータ解析部４が、この発明の計測情報解析手段に相当する。
【００２５】
ステップＳ５（発振コイル位置を特定）
図３のフローチャートを参照してコイル特定部５で行われる処理手順を説明する。このコイル位置特定部５は、この発明の位置特定手段に相当する。
【００２６】
ステップＴ１（仮想発振コイルの位置を設定）
仮想発振コイルをマルチチャンネルＳＱＵＩＤセンサ１を基準とする三次元の座標系の仮想位置に設定する。
【００２７】
ステップＴ２（仮想磁界データ群Ｅm を算出）
この仮想発振コイルには、前記電流強度ｉの電流が供給されているものとしたとき、マルチチャンネルＳＱＵＩＤセンサ１の計測される磁界をビオ・サバールの法則を用いて理論的に求める。ここで、マルチチャンネルＳＱＵＩＤセンサ１は、デュアー１ａ内に図示しない多数の磁気センサＳ1 〜Ｓm をもっている。例えば、デュアー１ａ内のある１つの磁気センサＳk では、ビオ・サバールの法則によれば次式（１）で表すことができる。
【００２８】
Ｈ＝ｓｉｎθ・ｉ・πａ2 ／（４πｒ2 ）          ・・・・（１）
Ｈ：磁気センサＳk で検出する仮想磁界データ
ｓｉｎθ：仮想発振コイルの位置と磁気センサＳk とがなす角度
ｉ：仮想発振コイルに流れる電流値
πａ2 ：仮想発振コイルのコイル径（発振コイルＣ１の径）
ｒ：仮想発振コイルの位置から磁気センサＳk までの距離
【００２９】
式（１）において、ｉとπａ2 は既に判っているので、未知数となるのはｒである。このｒは、仮想発振コイルの位置から磁気センサＳk までの距離なので、仮想発振コイルの位置をマルチチャンネルＳＱＵＩＤセンサ１を基準とする三次元の座標系に仮想的に決めれば、このｒがきまる。以上から、磁界Ｈを求めることができる。ここで、式（１）は、デュアー１ａ内の１つの磁気センサが検出する仮想磁界データであり、マルチチャンネルＳＱＵＩＤセンサ１は、この磁気センサを多数有している。したがって、各磁気センサＳ1 〜Ｓm で得られる仮想磁界データを全て求めることで、マルチチャンネルＳＱＵＩＤセンサ１で検出される仮想磁界データ群Ｅm を求めることができる。
【００３０】
ステップＴ３（最小自乗誤差を算出）
仮想発振コイルの仮想磁界データ群Ｅm と発振コイルＣ１の解析磁界データ群ｅm との最小自乗誤差を求める。この仮想磁界データ群Ｅm と解析磁界データ群ｅm との最小自乗誤差を求めることで、両データ群の一致の度合いを把握することができ、この最小自乗誤差の値が０に近い程、両データ群は一致する度合いが高いといえる。なお、この両データ群の一致する度合いが高いということは、仮想発振コイルと発振コイルＣ１との位置が一致している度合いが高いことを示すものである。
【００３１】
ステップＴ４（判定値以下か？）
ステップＴ２で算出された最小自乗誤差の値が、予め与えられた判定値以下になるか否かを判断する。判定値以下の場合は、仮想発振コイルの仮想位置を、発振コイルＣ１の位置として決定する。判定値より大きい場合は、ステップＴ５、Ｔ２、Ｔ３で行われる処理をステップＴ４での判定値以下となるまで繰り返し行う。ここでステップＴ５は、最小自乗誤差が小さくなる方向に仮想発振コイルの仮想位置を移動させる処理を行う。
【００３２】
ステップＳ６（全発振コイル終了？）
ステップＳ１〜Ｓ５で行われる処理を、被検体Ｍの頭部に装着された各発振コイルについて行う。これにより、マルチチャンネルＳＱＵＩＤセンサ１を基準とする三次元の座標系の発振コイルＣ１、Ｃ２、Ｃ３の位置を特定することができる。
【００３３】
上述した手段で求められた発振コイルＣ１、Ｃ２、Ｃ３の位置は、各発振コイルに実際に供給される電流の電流周波数と電流強度とを検出する。この電流周波数に基づいて、発振コイルからの磁界を計測して得た計測磁界データから解析磁界データ群を算出しているので、正確な位置を特定することができる。さらに、この発振コイルの正確な位置に基づいて、被検体Ｍの頭部の断層画像と推定された電流源とを重ねるので、被検体Ｍの頭部内に発生する電流源の位置を正確に把握することができる。
【００３４】
＜第２実施例＞
第２実施例に係る生体磁気計測装置の概略構成は、第１実施例の生体磁気計測装置の概略構成を示すブロック図のフィードバック部３と電流供給部２が異なるので、この部分のみ図４に示し、他の部分に関しては図１を参照する。また、被検体Ｍの頭部にｎ個の発振コイルが装着されているものとする。
【００３５】
図４に示す電流供給部２０は、ｎ個の電流供給源部ａ1 、ａ2 、ａ3 、・・、ａn から構成されている。この電流供給源部の数は、被検体Ｍの関心部位である頭部に装着される発振コイルの数と同数で構成される。それぞれの電流供給源部ａ1 〜ａn からフィードバック部３０を介して発振コイルＣ1 〜Ｃn に送出される。この電流値をフィードバック部３０で検出し、この検出した電流値をデータ解析部に送る
【００３６】
以下、図５のフローチャートを参照して説明する。
【００３７】
ステップＵ１（全供給電流値を検出）
電流供給部２０は、発振コイルＣ1 〜Ｃn にそれぞれ同時に電流供給源部ａ1 〜ａn から予めきめられた電流値の電流を送出する。この時、電流供給部２０と発振コイルＣ1 〜Ｃn との間の設置されたフィードバック部３０で、電流供給部２０内の電流供給源部ａ1 〜ａn から送出される電流値を個別同時に検出する。この検出された電流の値は、デジタル変換されデータ解析部４に送られる。
【００３８】
ここで、各発振コイルに対応する電流供給源部ａ1 〜ａn は、それぞれ異なった周波数および電流強度の電流が送出するので、フィードバック部３０で検出される電流の値も発振コイルＣ1 〜Ｃn によって異なる。例えば、発振コイルＣ1 では、周波数ｆ1 、電流強度ｉ1 の電流が検出され、発振コイルＣ2 では、周波数ｆ2 ，電流強度ｉ2 の電流が、・・・・発振コイルＣn では、周波数ｆn 、電流強度ｉn の電流が検出される。
【００３９】
ステップＵ２（全発振コイルからの磁界を計測）
発振コイルＣ1 〜Ｃn は、電流供給部２０から送られた電流によって、一斉に磁界を発生する。この磁界をマルチチャンネルＳＱＵＩＤセンサ１で計測する。この計測された磁界の計測磁界データ群は、全発振コイルからの磁界を含んでいる。なお、計測磁界データ群は、データ収集部６を介してデジタル変換されデータ解析部４に送られる。
【００４０】
ステップＵ３（全電流データ解析）
データ解析部４は、フィードバック部３０から送られた各電流に対して周波数解析を行うことで、発振コイルＣ1 〜Ｃn に実際に供給された電流の周波数ｆ1 〜ｆn および電流強度ｉ1 〜ｉn をそれぞれ知ることができる。前記フィードバック部３０とこのステップＵ３の処理を行うデータ解析部４とが、この発明の電流情報検出手段に相当する。
【００４１】
ステップＵ４（周波数ごとの計測磁界データ解析）
データ解析部４は、マルチチャンネルＳＱＵＩＤセンサ１で計測された計測磁界データ群に対して周波数解析を行う。この周波数解析は、ステップＵ３で求められた電流の各周波数ｆ1 〜ｆn に対して、それぞれの解析磁界データ群を求める。ここで、この計測磁界データ群に各周波数ｆ1 〜ｆn に基づく周波数解析を行うことで、各周波数ごとの解析磁界データ群を求めることができる。この周波数ごとの解析磁界データ群をコイル位置特定部５に送る。このステップＵ４の処理を行うデータ解析部４が、この発明の計測情報解析手段に相当する。
【００４２】
ステップＵ５（発振コイルＣ1 〜Ｃn の位置を特定）
コイル位置特定部５は、周波数ｆ1 〜ｆn に対する各解析磁界データ群に対して、第１実施例のステップＳ５で行われる処理をそれぞれ行う。これにより周波数ｆ1 〜ｆn に対するマルチチャンネルＳＱＵＩＤセンサ１を基準とする三次元の座標系の発振コイルＣ1 〜Ｃn の位置が特定する。このコイル位置特定部５は、この発明の位置特定手段に相当する。
【００４３】
上述した手段によって、発振コイルＣ1 〜Ｃn 位置を、一度の磁界計測によって同時に求めることができるので、短時間で全発振コイルの位置を正確に把握することができる。さらに、この発振コイルの位置に基づいて、被検体Ｍの頭部の断層画像と推定された電流源とを重ねるので、被検体Ｍの頭部内に発生する電流源の位置を正確に把握することができる。
【００４４】
【発明の効果】
以上の説明から明らかなように、この発明によれば、発振器に供給する電流を検出することで、実際に供給されている電流の電流周波数と電流強度とを知ることができ、この電流周波数と電流強度とに基づいて発振器の位置を特定するので、この特定された発振器の位置は、従来のように既知の電流情報に基づいて発振器の位置を特定する場合に比べて、より正確に特定することができる。さらに、この特定された発振器の位置を基準として、被検体の関心部位の断層画像と電流源とを重ねるので、被検体の関心部位に対する電流源の位置をより正確に把握することができる。
【図面の簡単な説明】
【図１】実施例に係る生体磁気計測装置の概略構成を示すブロック図である。
【図２】第１実施例で行われる処理手順を示すフローチャートである。
【図３】第１実施例のステップＳ５で行われる処理手順を示すフローチャートである。
【図４】第２実施例の電流供給部とフィードバック部を示すブロック図である。
【図５】第２実施例で行われる処理手順を示すフローチャートである。
【符号の説明】
１  …  マルチチャンネルＳＱＵＩＤセンサ
１ａ…  デュアー
２  …  電流供給部
３  …  フィードバック部
４  …  データ解析部
５  …  コイル位置特定部
Ｍ  …  被検体
Ｃ1 …  発振コイル
【図１】
0003791119.tif 000002
【図２】
0003791119.tif 000003
【図３】
0003791119.tif 000004
【図４】
0003791119.tif 000005
【図５】
0003791119.tif 000006