JPH11221198A - 生体磁気計測装置 - Google Patents
生体磁気計測装置Info
- Publication number
- JPH11221198A JPH11221198A JP10025412A JP2541298A JPH11221198A JP H11221198 A JPH11221198 A JP H11221198A JP 10025412 A JP10025412 A JP 10025412A JP 2541298 A JP2541298 A JP 2541298A JP H11221198 A JPH11221198 A JP H11221198A
- Authority
- JP
- Japan
- Prior art keywords
- biomagnetism
- information
- specific position
- magnetic sensor
- interest
- 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
【０００１】
【発明の属する技術分野】この発明は、被検体の関心部
位内の生体活動電流の状態を把握する生体磁気計測装置
に係り、特に、前記被検体の関心部位内の特定位置にお
ける生体活動電流の状態を正確に把握する技術に関す
る。
位内の生体活動電流の状態を把握する生体磁気計測装置
に係り、特に、前記被検体の関心部位内の特定位置にお
ける生体活動電流の状態を正確に把握する技術に関す
る。
【０００２】
【従来の技術】生体内に生じる生体活動電流に伴って微
小な生体磁気が生体から発生する。例えば、脳内の生体
活動電流に伴って発生する生体磁気は脳磁と呼ばれ、生
体に刺激を与えることにより発生する誘発脳磁や、α波
やてんかんのスパイク波のように脳から自然に発生する
自発脳磁などがある。
小な生体磁気が生体から発生する。例えば、脳内の生体
活動電流に伴って発生する生体磁気は脳磁と呼ばれ、生
体に刺激を与えることにより発生する誘発脳磁や、α波
やてんかんのスパイク波のように脳から自然に発生する
自発脳磁などがある。
【０００３】近年、生体内の微小な生体磁気を計測する
磁気センサとして、ＳＱＵＩＤ（Superconduc-ting Qua
ntum Interference Device：超電導量子干渉計）を用い
たマルチチャンネルＳＱＵＩＤセンサが開発されてい
る。このマルチチャンネルＳＱＵＩＤセンサは、デュア
ーと呼ばれる容器内に多数個の磁気センサを液体窒素な
どの冷媒に浸漬して収納している。この磁気センサは、
磁界を計測するコイル部と検出した磁界を電圧に変換す
るＳＱＵＩＤセンサ部とから構成されている。
磁気センサとして、ＳＱＵＩＤ（Superconduc-ting Qua
ntum Interference Device：超電導量子干渉計）を用い
たマルチチャンネルＳＱＵＩＤセンサが開発されてい
る。このマルチチャンネルＳＱＵＩＤセンサは、デュア
ーと呼ばれる容器内に多数個の磁気センサを液体窒素な
どの冷媒に浸漬して収納している。この磁気センサは、
磁界を計測するコイル部と検出した磁界を電圧に変換す
るＳＱＵＩＤセンサ部とから構成されている。
【０００４】このマルチチャンネルＳＱＵＩＤセンサを
被検体の関心部位である例えば頭部の外側に置き、脳内
に生じた生体活動電流から発生する微小な生体磁気をそ
のマルチチャンネルＳＱＵＩＤセンサ内の各磁気センサ
で無侵襲に計測することができる。この計測された計測
生体磁気情報から被検体の関心部位の生体活動電流の発
生した位置や方向などの状態を把握することができる。
被検体の関心部位である例えば頭部の外側に置き、脳内
に生じた生体活動電流から発生する微小な生体磁気をそ
のマルチチャンネルＳＱＵＩＤセンサ内の各磁気センサ
で無侵襲に計測することができる。この計測された計測
生体磁気情報から被検体の関心部位の生体活動電流の発
生した位置や方向などの状態を把握することができる。
【０００５】従来、被検体の関心部位内に発生した生体
活動電流の状態を把握するため次のように行われる。ま
ず、マルチチャンネルＳＱＵＩＤセンサと被検体の関心
部位との位置関係を捉える。この位置関係を捉える手法
として、特開平５ー２３７０６５号、特開平６ー７８８
９２５号などの種々の手法がある。
活動電流の状態を把握するため次のように行われる。ま
ず、マルチチャンネルＳＱＵＩＤセンサと被検体の関心
部位との位置関係を捉える。この位置関係を捉える手法
として、特開平５ー２３７０６５号、特開平６ー７８８
９２５号などの種々の手法がある。
【０００６】次に、被検体の関心部位内の特定位置に流
れる生体活動電流から発生する微小磁界に基づく生体磁
気情報が、強調されるような重み係数を算出する。この
重み係数は、マルチチャンネルＳＱＵＩＤセンサ内に収
納された磁気センサごとに算出される。この重み係数を
算出する手法として、S.E.Robinson,D.F.Rose,”Curren
t Source Image Estimation by Spatially Filtered MG
E" Biomaguetisum:Clinical aspects(1992) などの手法
がある。
れる生体活動電流から発生する微小磁界に基づく生体磁
気情報が、強調されるような重み係数を算出する。この
重み係数は、マルチチャンネルＳＱＵＩＤセンサ内に収
納された磁気センサごとに算出される。この重み係数を
算出する手法として、S.E.Robinson,D.F.Rose,”Curren
t Source Image Estimation by Spatially Filtered MG
E" Biomaguetisum:Clinical aspects(1992) などの手法
がある。
【０００７】被検体の関心部位からの生体磁気を計測し
て得た計測生体磁気情報と、特定位置からの生体磁気だ
けを強調する重み係数とのたたみ込み演算を行う。これ
によって、関心部位内の特定位置の生体活動電流から発
生した生体磁気だけが強調されるので、被検体の関心部
位内の任意の位置での生体活動電流の状態を把握するこ
とができる。
て得た計測生体磁気情報と、特定位置からの生体磁気だ
けを強調する重み係数とのたたみ込み演算を行う。これ
によって、関心部位内の特定位置の生体活動電流から発
生した生体磁気だけが強調されるので、被検体の関心部
位内の任意の位置での生体活動電流の状態を把握するこ
とができる。
【０００８】被検体の関心部位内の多数の位置につい
て、生体活動電流の状態を捉えれば、被検体の関心部位
全体に発生している生体活動電流の状態を把握すること
ができる。この生体活動電流を、ＭＲＩ装置等で得られ
た被検体の関心部位の断層画像上に重ねて表示させるこ
とで、関心部位の断層画像上の生体活動電流の状態を把
握することができる。
て、生体活動電流の状態を捉えれば、被検体の関心部位
全体に発生している生体活動電流の状態を把握すること
ができる。この生体活動電流を、ＭＲＩ装置等で得られ
た被検体の関心部位の断層画像上に重ねて表示させるこ
とで、関心部位の断層画像上の生体活動電流の状態を把
握することができる。
【０００９】
【発明が解決しようとする課題】しかしながら、このよ
うな構成を有する従来例の場合には、次のような問題が
ある。すなわち、従来の被検体の関心部位内の特定位置
に存在する生体活動電流からの生体磁気を強調するため
には、関心部位内の特定位置に設定されたδ関数で表さ
れる生体活動電流と、特定位置から各磁気センサまでを
磁界が伝達する伝達係数の総和とが等しくなるような、
伝達係数ごとに作用させる重み係数を求める必要があ
る。この重み係数を算出する際には、磁気センサごとの
伝達係数に基づく行列の逆行列と、特定位置にδ関数で
表される生体活動電流を設定した場合に各磁気センサで
仮想的に計測される仮想生体磁気情報とから算出する。
しかし、伝達係数に基づく行列の解が零に極めて近い場
合には、その逆行列の解が極端に大きくなる。この逆行
列に基づいて算出された重み係数が作用された計測生体
磁気情報は、その計測生体磁気情報に含まれるノイズま
でもが強調されたものとなるので、特定位置における生
体活動電流を正確に把握することができないという問題
がある。
うな構成を有する従来例の場合には、次のような問題が
ある。すなわち、従来の被検体の関心部位内の特定位置
に存在する生体活動電流からの生体磁気を強調するため
には、関心部位内の特定位置に設定されたδ関数で表さ
れる生体活動電流と、特定位置から各磁気センサまでを
磁界が伝達する伝達係数の総和とが等しくなるような、
伝達係数ごとに作用させる重み係数を求める必要があ
る。この重み係数を算出する際には、磁気センサごとの
伝達係数に基づく行列の逆行列と、特定位置にδ関数で
表される生体活動電流を設定した場合に各磁気センサで
仮想的に計測される仮想生体磁気情報とから算出する。
しかし、伝達係数に基づく行列の解が零に極めて近い場
合には、その逆行列の解が極端に大きくなる。この逆行
列に基づいて算出された重み係数が作用された計測生体
磁気情報は、その計測生体磁気情報に含まれるノイズま
でもが強調されたものとなるので、特定位置における生
体活動電流を正確に把握することができないという問題
がある。
【００１０】この発明は、このような事情に鑑みてなさ
れたものであって、計測された生体磁気情報に含まれる
ノイズに影響されることなく、被検体の関心部位内の特
定位置における生体活動電流の状態を、常に正確に求め
ることができる生体磁気計測装置を提供することを目的
とする。
れたものであって、計測された生体磁気情報に含まれる
ノイズに影響されることなく、被検体の関心部位内の特
定位置における生体活動電流の状態を、常に正確に求め
ることができる生体磁気計測装置を提供することを目的
とする。
【００１１】
【課題を解決するための手段】この発明は、このような
目的を達成するために、次のような構成をとる。すなわ
ち、請求項１に記載の発明は、被検体の関心部位に近接
配備した複数個の磁気センサによって、前記被検体の関
心部位内に流れる生体活動電流に伴って発生する生体磁
気を計測し、この計測で得られた計測生体磁気情報に基
づいて、被検体の関心部位内に流れる生体活動電流の状
態を把握する生体磁気計測装置において、（ａ）前記各
磁気センサで各々得られた計測生体磁気情報を収集する
計測情報収集手段と、（ｂ）前記被検体の関心部位内の
特定位置から磁気センサまでの間を磁界が伝達する度合
いを示す伝達関数に基づいて、磁気センサごとに伝達係
数を算出する伝達係数算出手段と、（ｃ）前記伝達係数
算出手段によって算出された磁気センサごとの伝達係数
に基づいた、行列の各要素が２つの伝達係数から算出さ
れる値で表される行列に所定値を加算する所定値加算手
段と、（ｄ）前記特定位置に既知の生体活動電流を設定
した際に前記各磁気センサで仮想的に計測される仮想生
体磁気情報と前記所定値を加算後の行列の逆行列とに基
づいて、前記特定位置から発生した生体磁気だけを強調
するための磁気センサごとの重み係数を算出する重み係
数算出手段と、（ｅ）前記重み係数算出手段で算出され
た磁気センサごとの重み係数を、前記計測情報収集手段
によって収集された磁気センサごとの計測生体磁気情報
に作用させて、特定位置からの生体磁気が強調された、
磁気センサごとの強調生体磁気情報を算出する強調生体
磁気情報算出手段とを備えたことを特徴とするものであ
る。
目的を達成するために、次のような構成をとる。すなわ
ち、請求項１に記載の発明は、被検体の関心部位に近接
配備した複数個の磁気センサによって、前記被検体の関
心部位内に流れる生体活動電流に伴って発生する生体磁
気を計測し、この計測で得られた計測生体磁気情報に基
づいて、被検体の関心部位内に流れる生体活動電流の状
態を把握する生体磁気計測装置において、（ａ）前記各
磁気センサで各々得られた計測生体磁気情報を収集する
計測情報収集手段と、（ｂ）前記被検体の関心部位内の
特定位置から磁気センサまでの間を磁界が伝達する度合
いを示す伝達関数に基づいて、磁気センサごとに伝達係
数を算出する伝達係数算出手段と、（ｃ）前記伝達係数
算出手段によって算出された磁気センサごとの伝達係数
に基づいた、行列の各要素が２つの伝達係数から算出さ
れる値で表される行列に所定値を加算する所定値加算手
段と、（ｄ）前記特定位置に既知の生体活動電流を設定
した際に前記各磁気センサで仮想的に計測される仮想生
体磁気情報と前記所定値を加算後の行列の逆行列とに基
づいて、前記特定位置から発生した生体磁気だけを強調
するための磁気センサごとの重み係数を算出する重み係
数算出手段と、（ｅ）前記重み係数算出手段で算出され
た磁気センサごとの重み係数を、前記計測情報収集手段
によって収集された磁気センサごとの計測生体磁気情報
に作用させて、特定位置からの生体磁気が強調された、
磁気センサごとの強調生体磁気情報を算出する強調生体
磁気情報算出手段とを備えたことを特徴とするものであ
る。
【００１２】請求項２に記載の発明は、請求項１に記載
の生体磁気計測装置において、（ｃ）前記所定値加算手
段は、前記計測生体磁気情報に含まれるノイズの大きさ
に応じた値を所定値として加算するものである。
の生体磁気計測装置において、（ｃ）前記所定値加算手
段は、前記計測生体磁気情報に含まれるノイズの大きさ
に応じた値を所定値として加算するものである。
【００１３】〔作用〕請求項１に記載の発明の作用は次
のとおりである。被検体の関心部位に近接配備された複
数個の磁気センサは、被検体の関心部位内に生じた生体
活動電流から発生する生体磁気を検出する。計測情報収
集手段は、磁気センサごとに検出された生体磁気を計測
生体磁気情報として収集する。伝達係数算出手段は、ビ
オ・サバールの法則に基づいて関心部位内の特定位置と
磁気センサとの間で磁界が伝達する度合いを示す伝達関
数を導きだすとともに、この伝達関数に特定位置と各磁
気センサとの位置関係を代入して、伝達係数を算出す
る。所定値加算手段は、それぞれ組み合わされた２つの
伝達係数によって求まる値を各要素とする行列に、所定
の値を加算して、この行列の解を零よりもさらに離れた
値にする。重み係数算出手段は、関心部位内の特定位置
に既知の生体活動電流を設定した場合に各磁気センサで
仮想的に計測される仮想生体磁気情報と、前記行列の逆
行列とを算出するとともに、この仮想生体磁気情報と逆
行列とに基づいて、磁気センサごとの重み係数を算出す
る。強調生体磁気情報算出手段は、この重み係数を各計
測生体磁気情報に作用させることで、特定位置から発生
した生体磁気だけが強調された強調生体磁気情報を算出
する。この強調生体磁気情報に基づく生体活動電流も強
調された状態となるので、関心部位内の特定位置におけ
る生体活動電流を正確に把握することができる。
のとおりである。被検体の関心部位に近接配備された複
数個の磁気センサは、被検体の関心部位内に生じた生体
活動電流から発生する生体磁気を検出する。計測情報収
集手段は、磁気センサごとに検出された生体磁気を計測
生体磁気情報として収集する。伝達係数算出手段は、ビ
オ・サバールの法則に基づいて関心部位内の特定位置と
磁気センサとの間で磁界が伝達する度合いを示す伝達関
数を導きだすとともに、この伝達関数に特定位置と各磁
気センサとの位置関係を代入して、伝達係数を算出す
る。所定値加算手段は、それぞれ組み合わされた２つの
伝達係数によって求まる値を各要素とする行列に、所定
の値を加算して、この行列の解を零よりもさらに離れた
値にする。重み係数算出手段は、関心部位内の特定位置
に既知の生体活動電流を設定した場合に各磁気センサで
仮想的に計測される仮想生体磁気情報と、前記行列の逆
行列とを算出するとともに、この仮想生体磁気情報と逆
行列とに基づいて、磁気センサごとの重み係数を算出す
る。強調生体磁気情報算出手段は、この重み係数を各計
測生体磁気情報に作用させることで、特定位置から発生
した生体磁気だけが強調された強調生体磁気情報を算出
する。この強調生体磁気情報に基づく生体活動電流も強
調された状態となるので、関心部位内の特定位置におけ
る生体活動電流を正確に把握することができる。
【００１４】請求項２に記載の発明の作用は次のとおり
である。前記所定値加算手段は、生体磁気を計測したと
きに含まれるノイズの大きさに応じた値を所定値として
加算する。この所定値が加算された行列に基づいた重み
係数が作用された計測生体磁気情報は、前記ノイズの影
響を受けた分が緩和されるように強調される。
である。前記所定値加算手段は、生体磁気を計測したと
きに含まれるノイズの大きさに応じた値を所定値として
加算する。この所定値が加算された行列に基づいた重み
係数が作用された計測生体磁気情報は、前記ノイズの影
響を受けた分が緩和されるように強調される。
【００１５】
【発明の実施の形態】以下、図面を参照してこの発明の
一実施例を説明する。図１は実施例に係る生体磁気計測
装置の概略構成を示すブロック図である。図中、符号１
１は磁気シールドルームであり、この磁気シールドルー
ム１１内に被検体Ｍが仰臥されるベッド２と、被検体Ｍ
の関心部位である例えば頭部に近接配備され、頭部内に
生じた生体活動電流から発生する微小な生体磁気を無侵
襲に計測するためのｍ個の磁気センサＳ1 〜Ｓm をデュ
アー１ａ内に収納するマルチチャンネルＳＱＵＩＤセン
サ１（以下、「磁束計１」と呼ぶ）とが設けられてい
る。
一実施例を説明する。図１は実施例に係る生体磁気計測
装置の概略構成を示すブロック図である。図中、符号１
１は磁気シールドルームであり、この磁気シールドルー
ム１１内に被検体Ｍが仰臥されるベッド２と、被検体Ｍ
の関心部位である例えば頭部に近接配備され、頭部内に
生じた生体活動電流から発生する微小な生体磁気を無侵
襲に計測するためのｍ個の磁気センサＳ1 〜Ｓm をデュ
アー１ａ内に収納するマルチチャンネルＳＱＵＩＤセン
サ１（以下、「磁束計１」と呼ぶ）とが設けられてい
る。
【００１６】図３に示すように、磁束計１は、デュアー
１ａ内にｍ個の磁気センサＳ1 〜Ｓ m が液体窒素等の冷
媒に侵漬して収納されて構成されている。この磁気セン
サＳ 1 〜Ｓm は、磁界を検出するコイル部と検出した磁
界を電圧や電流などに変換する図示しないＳＱＵＩＤセ
ンサ部とから構成されている。この磁気センサＳ1 〜Ｓ
m は、例えば３軸方向の磁界を正確に検出するため、３
つの磁気センサのコイル面が所定の角度で各々交わるよ
うに構成される、いわゆるベクトル型磁気センサであ
る。
１ａ内にｍ個の磁気センサＳ1 〜Ｓ m が液体窒素等の冷
媒に侵漬して収納されて構成されている。この磁気セン
サＳ 1 〜Ｓm は、磁界を検出するコイル部と検出した磁
界を電圧や電流などに変換する図示しないＳＱＵＩＤセ
ンサ部とから構成されている。この磁気センサＳ1 〜Ｓ
m は、例えば３軸方向の磁界を正確に検出するため、３
つの磁気センサのコイル面が所定の角度で各々交わるよ
うに構成される、いわゆるベクトル型磁気センサであ
る。
【００１７】磁気センサＳ1 〜Ｓ3 は、例えば、図４
（ａ）、（ｂ）に示す構成となる。図４（ａ）は、フイ
ルムＦ上にそれぞれ電気的に絶縁分離して形成された３
つのコイル対３１、３２、３３を備える。例えば、コイ
ル対３１は、三角角関数で表される２つの導電パターン
３１ａ、３１ｂの両端を直線の導電パターンＡ、Ｂで接
続したもので、導電パターン３１ｂに接続端子３１ｃが
形成されている。他のコイル対３２、３３も同様に導電
パターン３２ａ、３２ｂ、３３ａ、３３ｂで構成されて
いる。各コイル対３１、３２、３３の導電パターンはそ
れぞれ横方向に（２／３）πａだけシフトして配置され
ている。ここで、符号ａは後述するコイルボビンの半径
である。
（ａ）、（ｂ）に示す構成となる。図４（ａ）は、フイ
ルムＦ上にそれぞれ電気的に絶縁分離して形成された３
つのコイル対３１、３２、３３を備える。例えば、コイ
ル対３１は、三角角関数で表される２つの導電パターン
３１ａ、３１ｂの両端を直線の導電パターンＡ、Ｂで接
続したもので、導電パターン３１ｂに接続端子３１ｃが
形成されている。他のコイル対３２、３３も同様に導電
パターン３２ａ、３２ｂ、３３ａ、３３ｂで構成されて
いる。各コイル対３１、３２、３３の導電パターンはそ
れぞれ横方向に（２／３）πａだけシフトして配置され
ている。ここで、符号ａは後述するコイルボビンの半径
である。
【００１８】３つのコイル対３１、３２、３３が形成さ
れたフイルムＦを、直線の導電パターンＡ、Ｂが近接す
るように、コイルボビンに巻き付け固定する。その状態
を図４（ｂ）に示す。これにより、コイル対３１の導電
パターン３１ａは検出コイルを、導電パターン３１ｂは
補償コイルをそれぞれ形成し、各コイルは差動結合され
た状態になる。各コイル対３１、３２、３３の端子３１
ｃ、３２ｃ、３３ｃはＳＱＵＩＤセンサ部に接続され
る。なお、この発明は、このコイル対に限定するもので
はなく、例えば、各コイル面が、それぞれ直交する構成
をするものを使用してもよい。
れたフイルムＦを、直線の導電パターンＡ、Ｂが近接す
るように、コイルボビンに巻き付け固定する。その状態
を図４（ｂ）に示す。これにより、コイル対３１の導電
パターン３１ａは検出コイルを、導電パターン３１ｂは
補償コイルをそれぞれ形成し、各コイルは差動結合され
た状態になる。各コイル対３１、３２、３３の端子３１
ｃ、３２ｃ、３３ｃはＳＱＵＩＤセンサ部に接続され
る。なお、この発明は、このコイル対に限定するもので
はなく、例えば、各コイル面が、それぞれ直交する構成
をするものを使用してもよい。
【００１９】刺激装置４は、被検体Ｍに任意の電気的刺
激（あるいは音、光刺激など）を与え、被検体Ｍの関心
部位内に生体活動電流を発生させるためのものである。
この生体活動電流に伴って発生する微小磁界を計測生体
磁気情報として磁気センサＳ 1 〜Ｓm で個別に検出す
る。この計測生体磁気情報は、検出情報変換部５でデジ
タルデータに変換された後、計測情報収集部６に集めら
れる。計測情報収集部６は、本発明における計測情報収
集手段に相当する。
激（あるいは音、光刺激など）を与え、被検体Ｍの関心
部位内に生体活動電流を発生させるためのものである。
この生体活動電流に伴って発生する微小磁界を計測生体
磁気情報として磁気センサＳ 1 〜Ｓm で個別に検出す
る。この計測生体磁気情報は、検出情報変換部５でデジ
タルデータに変換された後、計測情報収集部６に集めら
れる。計測情報収集部６は、本発明における計測情報収
集手段に相当する。
【００２０】ポジショニング部３は、磁気センサＳ1 〜
Ｓm と被検体Ｍの関心部位との位置関係を把握するため
の装置である。例えば、被検体Ｍの関心部位の複数箇所
に小コイルを取付け、これらの小コイルにポジショニン
グ部３から給電する。そして、各コイルから発生した磁
界を磁気センサＳ1 〜Ｓm で検出することでコイル位置
情報を得る。このコイル位置情報を解析することで、磁
気センサＳ1 〜Ｓm と被検体Ｍの関心部位との位置関係
を把握する。なお、磁気センサＳ1 〜Ｓm と被検体Ｍの
関心部位との位置関係を把握するための手法は、これ以
外に、光ビームを被検体Ｍに照射して両者の位置関係を
把握するものや、あるいは、特開平５ー２３７０６５
号、特開平６ー７８８９２５号などに開示された種々の
手法が用いられる。
Ｓm と被検体Ｍの関心部位との位置関係を把握するため
の装置である。例えば、被検体Ｍの関心部位の複数箇所
に小コイルを取付け、これらの小コイルにポジショニン
グ部３から給電する。そして、各コイルから発生した磁
界を磁気センサＳ1 〜Ｓm で検出することでコイル位置
情報を得る。このコイル位置情報を解析することで、磁
気センサＳ1 〜Ｓm と被検体Ｍの関心部位との位置関係
を把握する。なお、磁気センサＳ1 〜Ｓm と被検体Ｍの
関心部位との位置関係を把握するための手法は、これ以
外に、光ビームを被検体Ｍに照射して両者の位置関係を
把握するものや、あるいは、特開平５ー２３７０６５
号、特開平６ー７８８９２５号などに開示された種々の
手法が用いられる。
【００２１】演算処理部７は、被検体Ｍの関心部位と磁
気センサＳ1 〜Ｓm の位置関係に基づいて、重み係数を
算出するとともに、この重み係数をメモリ８に記憶す
る。この重み係数を磁気センサＳ1 〜Ｓm で検出された
計測生体磁気情報に与えることで、被検体Ｍの関心部位
内の特定位置での生体活動電流の方向と大きさ等の情報
を得る。また、これらの情報を被検体Ｍの関心部位の断
層画像等に重ねた画像や生体活動電流の方向と大きさを
示す図をカラーモニタ９に表示したり、あるいはカラー
プリンタ１０に印刷出力する。なお、後述する動作説明
で明らかになるように、演算処理部７は、この発明にお
ける伝達係数算出手段、所定値加算手段、重み係数算出
手段、強調生体磁気情報算出手段に相当する機能を備え
る。
気センサＳ1 〜Ｓm の位置関係に基づいて、重み係数を
算出するとともに、この重み係数をメモリ８に記憶す
る。この重み係数を磁気センサＳ1 〜Ｓm で検出された
計測生体磁気情報に与えることで、被検体Ｍの関心部位
内の特定位置での生体活動電流の方向と大きさ等の情報
を得る。また、これらの情報を被検体Ｍの関心部位の断
層画像等に重ねた画像や生体活動電流の方向と大きさを
示す図をカラーモニタ９に表示したり、あるいはカラー
プリンタ１０に印刷出力する。なお、後述する動作説明
で明らかになるように、演算処理部７は、この発明にお
ける伝達係数算出手段、所定値加算手段、重み係数算出
手段、強調生体磁気情報算出手段に相当する機能を備え
る。
【００２２】以下、生体磁気計測装置で行われる各部の
処理を図２のフローチャートを参照して説明する。 ステップＳ１（関心部位の位置を検出） まず、頭部に装着された図示しない複数個の小コイルに
ポジショニング部３から電流を供給する。そして、各小
コイルから発生した磁界を磁気センサＳ1 〜Ｓ m で検出
する。検出情報変換部５は、この検出された磁界をデジ
タル変換して、コイル位置情報として計測情報収集部６
に送る。
処理を図２のフローチャートを参照して説明する。 ステップＳ１（関心部位の位置を検出） まず、頭部に装着された図示しない複数個の小コイルに
ポジショニング部３から電流を供給する。そして、各小
コイルから発生した磁界を磁気センサＳ1 〜Ｓ m で検出
する。検出情報変換部５は、この検出された磁界をデジ
タル変換して、コイル位置情報として計測情報収集部６
に送る。
【００２３】ステップＳ２（伝達関数を導出） 演算処理部７は、計測情報収集部６に送られてきたコイ
ル位置情報に基づいて、被検体Ｍと磁気センサＳ1 〜Ｓ
m との位置関係を求める。この位置関係を求める手法
は、従来と同様に最小自乗法を用いて算出することがで
きる。磁気センサＳ1 〜Ｓm と被検体Ｍとの位置関係が
求まると、演算処理部７は、例えば、図５に示すように
被検体Ｍの関心部位内における任意の位置を原点Ｏとす
る３次元の座標系を設定する。この座標系において、被
検体Ｍの関心部位内の例えば特定位置Ｐx から発生した
生体磁気が磁気センサＳ1 〜Ｓm まで伝達する度合いを
示す伝達関数を導出する。
ル位置情報に基づいて、被検体Ｍと磁気センサＳ1 〜Ｓ
m との位置関係を求める。この位置関係を求める手法
は、従来と同様に最小自乗法を用いて算出することがで
きる。磁気センサＳ1 〜Ｓm と被検体Ｍとの位置関係が
求まると、演算処理部７は、例えば、図５に示すように
被検体Ｍの関心部位内における任意の位置を原点Ｏとす
る３次元の座標系を設定する。この座標系において、被
検体Ｍの関心部位内の例えば特定位置Ｐx から発生した
生体磁気が磁気センサＳ1 〜Ｓm まで伝達する度合いを
示す伝達関数を導出する。
【００２４】演算処理部７は、図５に示すように、原点
Ｏからの距離ｒの特定位置Ｐx に生体活動電流が存在す
る場合に、特定位置Ｐx から磁気センサＳi までの間の
生体磁気の伝達の度合いを示す伝達関数Ｇi(r)を次のよ
うにして導出する。例えば、特定位置Ｐx に生体活動電
流Ｊ(r) （Ｊ(r) は電流密度を示す）が存在する場合に
磁気センサＳi で検出される生体磁気は、ビオ・サバー
ルの法則から次式（１）で表すことができる。
Ｏからの距離ｒの特定位置Ｐx に生体活動電流が存在す
る場合に、特定位置Ｐx から磁気センサＳi までの間の
生体磁気の伝達の度合いを示す伝達関数Ｇi(r)を次のよ
うにして導出する。例えば、特定位置Ｐx に生体活動電
流Ｊ(r) （Ｊ(r) は電流密度を示す）が存在する場合に
磁気センサＳi で検出される生体磁気は、ビオ・サバー
ルの法則から次式（１）で表すことができる。
【００２５】
【数１】
【００２６】式（１）において、μ0 は真空中の透磁率
であり、ｒi は原点Ｏから磁気センサＳi までの距離で
あり、ｅi は磁気センサＳi のコイル面を通過する単位
法線ベクトルである。したがって、Ｂi は、ｉ番目の磁
気センサで検出される磁束密度を示す。上記式（１）を
変形して次式（２）、（３）を導く。
であり、ｒi は原点Ｏから磁気センサＳi までの距離で
あり、ｅi は磁気センサＳi のコイル面を通過する単位
法線ベクトルである。したがって、Ｂi は、ｉ番目の磁
気センサで検出される磁束密度を示す。上記式（１）を
変形して次式（２）、（３）を導く。
【００２７】
【数２】
【００２８】式（２）に示すように、磁束密度Ｂi は、
生体活動電流Ｊ(r) とＧi(r)との内積を積分した式で表
すことができる。上記式（３）で示すＧi(r)が、特定位
置Ｐx から磁気センサＳi までの伝達関数である。ステ
ップＳ２は、本発明における伝達係数算出手段の機能に
相当する。なお、被検体Ｍの関心部位内の特定位置Ｐx
の原点Ｏからの距離ｒや磁気センサＳ1 〜Ｓm までの距
離の具体的な数値を式（３）に示す伝達係数Ｇi(r)に代
入して算出した値が、本発明における、関心部位内の特
定位置から各磁気センサまでの伝達係数である。以下、
被検体Ｍの関心部位内の特定位置Ｐx について、伝達関
数Ｇi(r)を用いて説明する。
生体活動電流Ｊ(r) とＧi(r)との内積を積分した式で表
すことができる。上記式（３）で示すＧi(r)が、特定位
置Ｐx から磁気センサＳi までの伝達関数である。ステ
ップＳ２は、本発明における伝達係数算出手段の機能に
相当する。なお、被検体Ｍの関心部位内の特定位置Ｐx
の原点Ｏからの距離ｒや磁気センサＳ1 〜Ｓm までの距
離の具体的な数値を式（３）に示す伝達係数Ｇi(r)に代
入して算出した値が、本発明における、関心部位内の特
定位置から各磁気センサまでの伝達係数である。以下、
被検体Ｍの関心部位内の特定位置Ｐx について、伝達関
数Ｇi(r)を用いて説明する。
【００２９】ステップＳ３（重み係数算出式を導出） 被検体Ｍの関心部位全体からの生体磁気を計測して得た
計測生体磁気情報の中から、特定位置Ｐx から発生した
生体磁気だけを強調できれば、関心部位内のあらゆる位
置に存在する生体活動電流の状態を正確に把握すること
ができる。つまり、特定位置Ｐx から磁気センサＳ1 〜
Ｓm にまで伝達してくる伝達関数Ｇi(r)に重み付けをす
れば、この伝達関数Ｇi(r)に基づいて伝達する生体磁気
が強調されることとなる。そこで、演算処理部７は、δ
関数（単位大きさの電流）で表される電流δ(r) を特定
位置Ｐx に仮想的に設定する。さらに、この電流δ(r)
と、特定位置Ｐx から磁気センサＳ1 〜Ｓm までの伝達
関数Ｇ1(r)〜Ｇm(r)の総和とを等しくするような重み係
数Ｗi を含む、次式（４）で表される評価関数ｆを設定
する。なお、この電流δ(r) を３軸方向の任意の方向に
設定することで、被検体Ｍの関心部位内に発生した生体
活動電流の任意の方向成分を強調することができる。
計測生体磁気情報の中から、特定位置Ｐx から発生した
生体磁気だけを強調できれば、関心部位内のあらゆる位
置に存在する生体活動電流の状態を正確に把握すること
ができる。つまり、特定位置Ｐx から磁気センサＳ1 〜
Ｓm にまで伝達してくる伝達関数Ｇi(r)に重み付けをす
れば、この伝達関数Ｇi(r)に基づいて伝達する生体磁気
が強調されることとなる。そこで、演算処理部７は、δ
関数（単位大きさの電流）で表される電流δ(r) を特定
位置Ｐx に仮想的に設定する。さらに、この電流δ(r)
と、特定位置Ｐx から磁気センサＳ1 〜Ｓm までの伝達
関数Ｇ1(r)〜Ｇm(r)の総和とを等しくするような重み係
数Ｗi を含む、次式（４）で表される評価関数ｆを設定
する。なお、この電流δ(r) を３軸方向の任意の方向に
設定することで、被検体Ｍの関心部位内に発生した生体
活動電流の任意の方向成分を強調することができる。
【００３０】
【数３】
【００３１】次に、この評価関数ｆを最小にするような
重み係数Ｗi を算出するために、式（４）を微分変形す
ることで、次式（５）に示す線形方程式を導き出す。
重み係数Ｗi を算出するために、式（４）を微分変形す
ることで、次式（５）に示す線形方程式を導き出す。
【００３２】Ｗ＝Γ-1Ｂ ・・・（５）
【００３３】式（５）において、Ｗは重み係数Ｗi を要
素とする行列、Γ-1はΓij＝∫Ｇi(r)Ｇj(r)dr3 を要素
とする行列、ＢはＢ＝∫Ｇi(r)δ(r)dr3を要素とする行
列を示す。
素とする行列、Γ-1はΓij＝∫Ｇi(r)Ｇj(r)dr3 を要素
とする行列、ＢはＢ＝∫Ｇi(r)δ(r)dr3を要素とする行
列を示す。
【００３４】演算処理部７は、予めオペレータによって
図示しない入力手段によって入力された、生体磁気を測
定するときに重畳するノイズに相当する値、例えば、λ
Ｉを上記式（４）のΓに加算する。λは入力されたノイ
ズの大きさを示し、Ｉは単位行列を示す。λＩを加算後
の線形方程式を重み係数算出式Ｗとして次式（６）に示
す。
図示しない入力手段によって入力された、生体磁気を測
定するときに重畳するノイズに相当する値、例えば、λ
Ｉを上記式（４）のΓに加算する。λは入力されたノイ
ズの大きさを示し、Ｉは単位行列を示す。λＩを加算後
の線形方程式を重み係数算出式Ｗとして次式（６）に示
す。
【００３５】Ｗ＝（Γ＋λＩ）-1Ｂ ・・・（６）
【００３６】上記値λＩをΓに加算することで、行列Γ
の解が極めて零に近い値になる場合であっても、行列
（Γ＋λＩ）の解は零よりも離れた値となる。つまり、
正則化された行列（Γ＋λＩ）は、いわゆる行列の性質
が安定した状態となる。演算処理部７は、この重み係数
算出式Ｗをメモリ８に記憶する。ステップＳ３は、本発
明における所定値加算手段の機能に相当する。
の解が極めて零に近い値になる場合であっても、行列
（Γ＋λＩ）の解は零よりも離れた値となる。つまり、
正則化された行列（Γ＋λＩ）は、いわゆる行列の性質
が安定した状態となる。演算処理部７は、この重み係数
算出式Ｗをメモリ８に記憶する。ステップＳ３は、本発
明における所定値加算手段の機能に相当する。
【００３７】ステップＳ４（生体磁気を計測） 演算処理部７は重み係数算出式Ｗの導出が終了すると、
刺激装置４に刺激開始の指示をする。刺激装置４は、被
検体Ｍに刺激を与えて、関心部位内に生体活動電流を生
じさせる。被検体Ｍの関心部位からは、この生体活動電
流に伴って微小磁界が発生する。検出情報変換部５は、
磁気センサＳ1 〜Ｓm で各々検出された生体磁気をデジ
タル変換後、計測生体磁気情報Ｍ1 〜Ｍm として計測情
報収集部６に送る。
刺激装置４に刺激開始の指示をする。刺激装置４は、被
検体Ｍに刺激を与えて、関心部位内に生体活動電流を生
じさせる。被検体Ｍの関心部位からは、この生体活動電
流に伴って微小磁界が発生する。検出情報変換部５は、
磁気センサＳ1 〜Ｓm で各々検出された生体磁気をデジ
タル変換後、計測生体磁気情報Ｍ1 〜Ｍm として計測情
報収集部６に送る。
【００３８】ステップＳ５（重み係数を算出） 演算処理部７は、オペレータから予め指定されている被
検体Ｍの指定位置Ｐ1〜Ｐn の３次元座標上の位置を重
み係数算出式Ｗに代入して、各指定位置Ｐ1 〜Ｐn から
の生体磁気を強調するための重み係数を算出する。例え
ば、指定位置Ｐ1 については、各磁気センサＳ1 〜Ｓm
の計測生体磁気情報に重み付けをするための重み係数Ｐ
1W1 〜Ｐ1Wm が求まる。同様にして、指定位置Ｐ2 〜Ｐ
n についても重み係数Ｐ2W1 〜Ｐ2Wm ，…，ＰnW1 〜Ｐ
nWm が求まる。この重み係数をメモリ８に記憶する。ス
テップＳ５は、本発明における重み係数算出手段の機能
に相当する。
検体Ｍの指定位置Ｐ1〜Ｐn の３次元座標上の位置を重
み係数算出式Ｗに代入して、各指定位置Ｐ1 〜Ｐn から
の生体磁気を強調するための重み係数を算出する。例え
ば、指定位置Ｐ1 については、各磁気センサＳ1 〜Ｓm
の計測生体磁気情報に重み付けをするための重み係数Ｐ
1W1 〜Ｐ1Wm が求まる。同様にして、指定位置Ｐ2 〜Ｐ
n についても重み係数Ｐ2W1 〜Ｐ2Wm ，…，ＰnW1 〜Ｐ
nWm が求まる。この重み係数をメモリ８に記憶する。ス
テップＳ５は、本発明における重み係数算出手段の機能
に相当する。
【００３９】ステップＳ６（重み係数を重畳） 演算処理部７は、指定位置Ｐk(k=1 〜n)の重み係数ＰkW
1 〜ＰkWm をメモリ８から呼び出すとともに、この重み
係数ＰkW1 〜ＰkWm を、計測情報収集部６に収集された
計測生体磁気情報Ｍ1 〜Ｍm にそれぞれ重畳する。これ
によって、指定位置Ｐk からの生体磁気だけが強調され
た強調生体磁気情報ＤPk1 〜ＤPkm が求まる。さらに、
強調生体磁気情報ＤPk1 〜ＤPkm の総和を算出すること
で、指定位置Ｐk における強調された生体活動電流Ｄk
を算出することができる。ステップＳ６で行われる計算
は、次式（７）で表される。
1 〜ＰkWm をメモリ８から呼び出すとともに、この重み
係数ＰkW1 〜ＰkWm を、計測情報収集部６に収集された
計測生体磁気情報Ｍ1 〜Ｍm にそれぞれ重畳する。これ
によって、指定位置Ｐk からの生体磁気だけが強調され
た強調生体磁気情報ＤPk1 〜ＤPkm が求まる。さらに、
強調生体磁気情報ＤPk1 〜ＤPkm の総和を算出すること
で、指定位置Ｐk における強調された生体活動電流Ｄk
を算出することができる。ステップＳ６で行われる計算
は、次式（７）で表される。
【００４０】Ｄk ＝Σ（ＰkWm ・Ｍm ）・・・（７）
【００４１】式（７）による強調生体磁気情報の算出を
指定位置Ｐ1 〜Ｐn について行うことで、各指定位置Ｐ
1 〜Ｐn での生体活動電流Ｄ1 〜Ｄn が求まる。ステッ
プＳ６は、本発明における強調生体磁気情報算出手段の
機能に相当する。
指定位置Ｐ1 〜Ｐn について行うことで、各指定位置Ｐ
1 〜Ｐn での生体活動電流Ｄ1 〜Ｄn が求まる。ステッ
プＳ６は、本発明における強調生体磁気情報算出手段の
機能に相当する。
【００４２】ステップＳ７（アローマップを表示） ステップＳ６で求められた生体活動電流Ｄ1 〜Ｄn は、
電流の向き及び大きさを持つベクトルとして扱うことが
できるので、被検体Ｍの関心部位内の指定位置Ｐ1 〜Ｐ
n に存在する生体活動電流として表示することができ
る。例えば、演算処理部７は、図６に示すように指定位
置Ｐ1 〜Ｐn を平面上に表示するとともに、各指定位置
Ｐ1 〜Ｐn にベクトル８１を表示する、いわゆるアロー
マップ８０をカラーモニタ９に表示したり、カラープリ
ンタ１０に印刷出力する。
電流の向き及び大きさを持つベクトルとして扱うことが
できるので、被検体Ｍの関心部位内の指定位置Ｐ1 〜Ｐ
n に存在する生体活動電流として表示することができ
る。例えば、演算処理部７は、図６に示すように指定位
置Ｐ1 〜Ｐn を平面上に表示するとともに、各指定位置
Ｐ1 〜Ｐn にベクトル８１を表示する、いわゆるアロー
マップ８０をカラーモニタ９に表示したり、カラープリ
ンタ１０に印刷出力する。
【００４３】上述した手段で求められたアローマップを
観察することで、被検体Ｍの関心部位内の任意の位置で
の生体活動電流の向き、大きさの状態を正確に把握する
ことができる。また、被検体Ｍの関心部位の断層画像上
に重ねることで、断層画像上の生体活動電流の位置、向
き、大きさを把握することもできる。
観察することで、被検体Ｍの関心部位内の任意の位置で
の生体活動電流の向き、大きさの状態を正確に把握する
ことができる。また、被検体Ｍの関心部位の断層画像上
に重ねることで、断層画像上の生体活動電流の位置、向
き、大きさを把握することもできる。
【００４４】この実施例は以下のように変形実施するこ
とが可能である。 （１）上記実施例では、３つの磁気センサのコイル面が
それぞれ所定の角度で交わるように構成されたベクトル
型磁気センサについて説明したが、本発明はこれに限定
されるものではなく、例えば、各々の磁気センサのコイ
ル面が交わることなく配備されたものや、２つの磁気セ
ンサのコイル面が交わるように構成されたものでもよ
い。
とが可能である。 （１）上記実施例では、３つの磁気センサのコイル面が
それぞれ所定の角度で交わるように構成されたベクトル
型磁気センサについて説明したが、本発明はこれに限定
されるものではなく、例えば、各々の磁気センサのコイ
ル面が交わることなく配備されたものや、２つの磁気セ
ンサのコイル面が交わるように構成されたものでもよ
い。
【００４５】（２）上述した式（６）では、所定値とし
て大きさλの単位行列をΓに加算した後、逆行列を計算
したが、例えば、所定値として大きさλの対象行列、い
わゆるペナルティ行列にすることもできる。
て大きさλの単位行列をΓに加算した後、逆行列を計算
したが、例えば、所定値として大きさλの対象行列、い
わゆるペナルティ行列にすることもできる。
【００４６】
【発明の効果】以上の説明から明らかなように、請求項
１の発明によれば、被検体の関心部位内の特定位置から
発生する生体磁気が、磁気センサまで伝達する度合いを
示す伝達関数によって、特定位置から磁気センサごとの
伝達係数を求める。これらの伝達係数に基づいた行列に
所定値を加算しているので、所定値を加算する前の行列
の解が極めて零に近い場合に、この行列の逆行列に基づ
く重み係数を作用させた計測生体磁気情報が、その計測
生体磁気情報に含まれるノイズまでもが強調されたよう
な不適切な値になるのを防止することができる。つま
り、特定位置から発生する生体磁気を強調するのに適切
な重み係数を常に求めることができる。その結果、この
重み係数によって重み付けされた強調計測生体磁気情報
によって、被検体の関心部位内の特定位置での生体活動
電流を常に正確に把握することができる。
１の発明によれば、被検体の関心部位内の特定位置から
発生する生体磁気が、磁気センサまで伝達する度合いを
示す伝達関数によって、特定位置から磁気センサごとの
伝達係数を求める。これらの伝達係数に基づいた行列に
所定値を加算しているので、所定値を加算する前の行列
の解が極めて零に近い場合に、この行列の逆行列に基づ
く重み係数を作用させた計測生体磁気情報が、その計測
生体磁気情報に含まれるノイズまでもが強調されたよう
な不適切な値になるのを防止することができる。つま
り、特定位置から発生する生体磁気を強調するのに適切
な重み係数を常に求めることができる。その結果、この
重み係数によって重み付けされた強調計測生体磁気情報
によって、被検体の関心部位内の特定位置での生体活動
電流を常に正確に把握することができる。
【００４７】また、請求項２の発明によれば、生体磁気
の計測時に含まれるノイズの大きさを所定値としている
ので、重み係数によって重み付けされた計測生体磁気情
報は、そのノイズの大きさ分だけが低減されたように強
調される。この強調計測生体磁気情報によって、被検体
の関心部位内の特定位置での生体活動電流を常により正
確に把握することができる。
の計測時に含まれるノイズの大きさを所定値としている
ので、重み係数によって重み付けされた計測生体磁気情
報は、そのノイズの大きさ分だけが低減されたように強
調される。この強調計測生体磁気情報によって、被検体
の関心部位内の特定位置での生体活動電流を常により正
確に把握することができる。
【図１】実施例に係る生体磁気計測装置の概略構成を示
すブロック図である。
すブロック図である。
【図２】実施例で行われる処理手順を示すフローチャー
トである。
トである。
【図３】実施例の磁束計と被検体の関心部位との拡大を
示す概略構成拡大図である。
示す概略構成拡大図である。
【図４】実施例の磁気センサのコイル部の構成を示す図
である。
である。
【図５】各磁気センサと特定位置との関係を示す図であ
る。
る。
【図６】実施例で表示されるアローマップを示す図であ
る。
る。
１ … 磁束計 １ａ… デュアー ３ … ポジショニング部 ４ … 刺激装置 ５ … 検出情報変換部 ６ … 計測情報収集部 ７ … 演算処理部 ８ … メモリ Ｍ … 被検体 Ｓi … 磁気センサ Ｐx … 特定位置 ｒi-ｒ … 特定位置から磁気センサまでの距離
Claims (2)
- 【請求項１】 被検体の関心部位に近接配備した複数個
の磁気センサによって、前記被検体の関心部位内に流れ
る生体活動電流に伴って発生する生体磁気を計測し、こ
の計測で得られた計測生体磁気情報に基づいて、被検体
の関心部位内に流れる生体活動電流の状態を把握する生
体磁気計測装置において、（ａ）前記各磁気センサで各
々得られた計測生体磁気情報を収集する計測情報収集手
段と、（ｂ）前記被検体の関心部位内の特定位置から磁
気センサまでの間を磁界が伝達する度合いを示す伝達関
数に基づいて、磁気センサごとに伝達係数を算出する伝
達係数算出手段と、（ｃ）前記伝達係数算出手段によっ
て算出された磁気センサごとの伝達係数に基づいた、行
列の各要素が２つの伝達係数から算出される値で表され
る行列に所定値を加算する所定値加算手段と、（ｄ）前
記特定位置に既知の生体活動電流を設定した際に前記各
磁気センサで仮想的に計測される仮想生体磁気情報と前
記所定値を加算後の行列の逆行列とに基づいて、前記特
定位置から発生した生体磁気だけを強調するための磁気
センサごとの重み係数を算出する重み係数算出手段と、
（ｅ）前記重み係数算出手段で算出された磁気センサご
との重み係数を、前記計測情報収集手段によって収集さ
れた磁気センサごとの計測生体磁気情報に作用させて、
特定位置からの生体磁気が強調された、磁気センサごと
の強調生体磁気情報を算出する強調生体磁気情報算出手
段とを備えたことを特徴とする生体磁気計測装置。 - 【請求項２】 請求項１に記載の生体磁気計測装置にお
いて、（ｃ）前記所定値加算手段は、前記計測生体磁気
情報に含まれるノイズの大きさに応じた値を所定値とし
て加算する生体磁気計測装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP10025412A JPH11221198A (ja) | 1998-02-06 | 1998-02-06 | 生体磁気計測装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP10025412A JPH11221198A (ja) | 1998-02-06 | 1998-02-06 | 生体磁気計測装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH11221198A true JPH11221198A (ja) | 1999-08-17 |
Family
ID=12165229
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP10025412A Pending JPH11221198A (ja) | 1998-02-06 | 1998-02-06 | 生体磁気計測装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH11221198A (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017176464A (ja) * | 2016-03-30 | 2017-10-05 | 株式会社日立ハイテクノロジーズ | 磁界計測装置、磁界計測方法、磁界計測プログラム |
-
1998
- 1998-02-06 JP JP10025412A patent/JPH11221198A/ja active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017176464A (ja) * | 2016-03-30 | 2017-10-05 | 株式会社日立ハイテクノロジーズ | 磁界計測装置、磁界計測方法、磁界計測プログラム |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6735460B2 (en) | Biomagnetic field measuring method and apparatus | |
| US6681131B2 (en) | Apparatus for measuring bio-magnetic fields | |
| US5269325A (en) | Analysis of biological signals using data from arrays of sensors | |
| US6374131B1 (en) | Biomagnetism measuring method and apparatus | |
| JP3473210B2 (ja) | 生体磁気計測装置 | |
| JP2002320596A (ja) | 生体磁場計測装置 | |
| WO2013088413A1 (en) | Removal of background in mpi | |
| Weinberg et al. | Magnetic localisation of intracranial dipoles: simulation with a physical model | |
| Allescher et al. | Biomagnetic 3-dimensional spatial and temporal characterization of electrical activity of human stomach | |
| JPH11221198A (ja) | 生体磁気計測装置 | |
| Rosen et al. | A study of the vector magnetocardiographic waveform | |
| JPH04319334A (ja) | 生体磁気イメージング装置 | |
| JP3407520B2 (ja) | 生体磁気計測装置 | |
| WO2010077390A1 (en) | Whole body electromagnetic detection system | |
| JPH08266499A (ja) | 生体磁気計測装置 | |
| Adachi et al. | A SQUID biomagnetometer system for measurement of spinal cord evoked magnetic fields | |
| JP4006543B2 (ja) | 生体活動電流源推定装置 | |
| JPH04303416A (ja) | 生体磁気計測装置 | |
| JP3409551B2 (ja) | 生体磁気計測装置 | |
| JP3298312B2 (ja) | 生体活動電流源推定装置 | |
| JPH03251226A (ja) | 生体磁気計測法 | |
| JP3814923B2 (ja) | 生体磁気計測装置 | |
| JP2795211B2 (ja) | 生体磁気計測装置 | |
| JP3591121B2 (ja) | 生体磁気計測装置 | |
| JPH10211181A (ja) | 生体活動電流源推定装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040603 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20060331 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060418 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20060808 |