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

生体磁気計測装置

Info

Publication number
JPH10211183A
JPH10211183A JP9019518A JP1951897A JPH10211183A JP H10211183 A JPH10211183 A JP H10211183A JP 9019518 A JP9019518 A JP 9019518A JP 1951897 A JP1951897 A JP 1951897A JP H10211183 A JPH10211183 A JP H10211183A
Authority
JP
Japan
Prior art keywords
oscillation
coils
coil
magnetic field
oscillation coils
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.)
Granted
Application number
JP9019518A
Other languages
English (en)
Other versions
JP3661329B2 (ja
Inventor
Shoichi Okamura
昇一 岡村
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shimadzu Corp
Original Assignee
Shimadzu Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shimadzu Corp filed Critical Shimadzu Corp
Priority to JP01951897A priority Critical patent/JP3661329B2/ja
Publication of JPH10211183A publication Critical patent/JPH10211183A/ja
Application granted granted Critical
Publication of JP3661329B2 publication Critical patent/JP3661329B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

(57)【要約】 【課題】 複数の関心領域が存在する場合であっても、
磁束計を移動させることなく、被検体の相対的な位置関
係を精度良く特定できる生体磁気計測装置を提供する。 【解決手段】 コンピュータ6から、電流供給部7に対
して、n個の発振コイルC1〜Cnに対し、順次所定の時間
間隔で所定強度の電流を出力するよう指示する。発振コ
イル位置算出部3は、磁束計S1〜Smでそれぞれ検出さ
れ、各発振コイルC1〜Cn個々に生じた磁場データをデー
タ収集ユニット2より入力し、周知の最小二乗法などに
より、磁束計S1〜Smに対する各発振コイルC1〜Cn個々相
対位置を算出する。発振コイル選択部4は、求めた各発
振コイルC1〜Cn個々の相対位置の正しさを示すgoodness
-of-fit (以下「GOF」という。)値を各発振コイル
C1〜Cn毎に求め、求めたGOF値を用いてMRI画像と
の位置合わせに用いる各発振コイルC1〜Cnを選択する。
発振コイルは、GOF値が所定値以上のもので任意の組
合せにより得られる三角形の面積が最大となる3つを選
択する。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、被検体内の生体活
動電流源に伴って発生する微小磁界を計測し、その計測
データに基づいて前記被検体内の生体活動電流源を求め
る生体磁気計測装置に関する。
【0002】
【従来の技術】近年の超伝導デバイス技術の発展に伴
い、SQUID(Superconducting QUantum Interferen
ce Device )と呼ばれる高感度な磁束計を利用した生体
磁気計測装置が、医療診断装置の一つとして実用化され
つつあり、脳機能の解明や循環器疾患の診断に役立つも
のと期待されている。
【0003】この生体磁気計測装置では、計測した磁場
データに基づき、例えば、最小自乗法や最小ノルム法等
によって、磁束計を基準とした座標系における生体活動
電流源の位置、向き、大きさなどの推定がなされる(Ju
kka Sarvas "Basic mathemtical and electromagneti
c concepts of the biomagnetic inverse problem" ,Ph
ys. Med. Biol., 1987, vol.32, No.1, 11-22, Printed
by the UK )。
【0004】一方、得られた生体活動電流源は、磁気共
鳴断層撮影装置(MRI装置)で得られるMRI画像や
X線CT装置で得られるX線断層画像などの医用画像上
に並記されることで、生体内の患部等の物理的位置を特
定することが可能となるため、磁束計を基準とした座標
系における生体活動電流源の位置情報と、医用画像との
位置関係を正確に把握することが重要である。
【0005】このため、鼻根部や両耳下などの頭部表面
の明確な位置にPROBE POSITIONINDICATOR とよばれる磁
場発生源を配置し、それによって生体活動電流源と被検
体との位置関係を求める以下のような方法 (1) S.Ahlfors, et al, "MAGNETOMETER POSITION IND
ICATOR FOR MULTICHANNEL MEG", Advances in Biomagne
tism, Edited by S.J.Williamson et al,Plenum Press,
New York 693-696, 1989 (2)Neuromag-122 Preliminary Technical Data, August
1991 (3) 「生体磁気測定を行うための装置及び方法」(特開
平1-503603 号) (4) 「生体磁場測定装置の位置検出装置」(特公平5-55
126 号) が提案されている。
【0006】これらの方法では、被検体の体表面に貼り
付けられた3つ或いはそれ以上の発振コイルの内、まず
1つ目の発振コイルに直流電流が与えられ、その発振コ
イルから発せられる磁場を互いにその位置関係が既知の
複数の磁束計によって検出され、発振コイルに与えた電
流の強さと各磁束計で検出した磁場の強さ、及び各磁束
計間の位置関係から磁束計群に対する上記1つ目の発振
コイルの位置が求められる。そして、この操作を2つ目
以降の発振コイルに順次適用し、発振コイルすべての位
置を求め、磁束計群に対する被検体の位置が決定され
る。
【0007】ここで、複数種の検査に際して、測定対象
となる脳内関心領域が複数存在する場合が多々生じる
が、かかる場合、磁束計を関心領域毎に最適な位置に移
動させ磁場計測を行うか、または、当該複数の関心領域
が感度範囲内に入るよう磁束計を配設する必要が生じ
る。
【0008】
【発明が解決しようとする課題】しかしながら、関心領
域毎に磁束計の位置合わせをしていたのでは、それぞれ
の磁束計の位置毎に発振コイルとの相対位置を算出する
必要が生じ、多大な労力と検査時間が必要となる。
【0009】また、複数の関心領域が感度範囲内に入る
よう磁束計を位置決めした場合、発振コイルと磁束計と
の相対位置を求める動作は一度で済むが、関心領域から
大きく離れた磁束計の一部では、正確な磁場データが得
られず、このため磁束計に対する発振コイル、すなわち
被検体との相対位置を精度良く求めることができない。
【0010】本発明は、上記課題を解決するために創案
されたもので、複数の関心領域が存在する場合であって
も、磁束計を移動させることなく、被検体の相対的な位
置関係を精度良く特定できる生体磁気計測装置を提供す
ることを目的とする。
【0011】
【課題を解決するための手段】上記目的を達成するため
に、本発明は、被検体内の生体活動電流源に伴って発生
する微小磁界を計測し、その計測データに基づいて前記
被検体内の生体活動電流源を求める生体磁気計測装置で
あって、前記被検体に付着される複数の発振コイルと、
この複数の発振コイルに所定電流を供給する電流供給手
段と、各発振コイルから生じた磁場を計測する磁束計
と、各発振コイル毎に得られた磁場データに基づき前記
磁束計に対する各発振コイルの相対位置を算出する磁場
解析手段と、得られた各発振コイルの相対位置からその
位置算出誤差を求め、当該誤差が発振コイルの位置算出
に際し許容範囲にあるものを選択するコイル選択手段
と、を備えることで、選択された発振コイルの位置情報
からMRI等の被検体画像との位置合わせを行うよう構
成したことを特徴とする。
【0012】前記コイル選択手段は、前記位置算出誤差
が許容範囲にある発振コイルのうち、任意の組合せによ
り得られる三角形の面積を算出し、当該面積が最大とな
る3つの発振コイルを選択することを特徴とする。
【0013】また、前記コイル選択手段は、位置算出誤
差が許容範囲にある発振コイルのうち、任意の組合せに
より得られる三角錐の体積が最大となる組合せを算出
し、当該体積が最大となる4つの発振コイルを選択する
ことを特徴とする。
【0014】前記コイル選択手段は、前記位置算出誤差
をgoodness-of-fit 値として求めることを特徴とする。
【0015】前記コイル選択手段は、前記位置算出誤差
を、「測定磁場」と「解析位置からの理論磁場」との相
関係数γとして求めることを特徴とする。
【0016】
【発明の実施の形態】本発明の一実施形態を図1〜図4
に基づいて説明する。図1は、本発明の一実施形態であ
る生体磁気計測装置の概略構成図である。同図におい
て、センサーユニット1は、ピックアップコイルとSQ
UIDからなる複数の高感度な磁束計S1〜Smが、デュア
ーの中に冷媒とともに収納されており、生体活動電流源
の検出に先立って、被検体Mの頭部に近接配備される。
【0017】発振コイルC1〜Cnは、鼻根部、両耳下等、
被検体Mを特定する上で特徴となる部分に付着されてお
り、それぞれ、例えば図2(a) のようにセラミック板な
ど絶縁体で形成された基板31に金属を印刷してコイル
部32を形成したコイルCや、図2(b) のようにボビン
33に金属ワイヤ34を巻いて形成したコイルC’が使
用される。電流供給部7は、上記各発振コイルC1〜Cn個
々に所定の既知電流を供給するが、例えば、周波数の異
なる交流電流を個別に指定した強度で同時に出力するよ
う構成しても良いし、所定の時間間隔で各発振コイルC1
〜Cnに所定電流を別々に供給するよう構成しても良い。
前者の場合、発振コイルC1〜Cnの位置計測時間を大幅に
縮小できるが、別途周波数解析を行い、個々の発振コイ
ルCiが磁束計S1〜Smに与えた磁場強度を算出する必要が
生じる。以下に示す実施形態は、所定の時間間隔で各発
振コイルC1〜Cnに所定電流を別々に供給する構成であ
る。
【0018】データ収集ユニット2は、磁束計S1〜Smで
計測された各発振コイルC1〜Cnから生じる磁場データを
A/D変換してコンピュータ6の発振コイル位置算出部
3に出力する。
【0019】コンピュータ6は、計測された磁場データ
の解析及び電流供給部の動作制御を行うもので、大きく
発振コイル位置算出部3、発振コイル選択部4、及び磁
場解析部5からなる。
【0020】発振コイル位置算出部3は、データ収集ユ
ニット2から出力された、各発振コイルC1〜Cnそれぞれ
に起因する磁束計S1〜Smで得られた磁場データから、周
知の最小二乗法等により磁束計S1〜Smに対する各発振コ
イルC1〜Cnの相対位置を算出する。
【0021】発振コイル選択部4は、発振コイル位置算
出部3で算出された磁束計S1〜Smに対する各発振コイル
C1〜Cnの相対位置についての誤差等からMRI画像など
とのマッチングに用いる発振コイルを選択する。
【0022】磁場解析部5は、上記した周知の最小二乗
法などにより生体活動電流源の推定演算を行うと共に、
発振コイル位置算出部3で求められ、発振コイル選択部
4で選択された発振コイルについての位置情報を、画像
記憶部7から読み出したMRI画像上の鼻根部、両耳下
等、被検体Mの特徴点に対応づけ、求めた生体活動電流
源に関する位置情報を、ここで対応づけた位置関係をも
とにMRI画像上に重ねてモニタ9に表示すると共に、
必要に応じてMOD(光磁気ディスク)などの外部メモ
リ8に保存し、或いは不図示のプリンタ等に出力する。
【0023】次に、本実施形態の作用をコンピュータ6
の動作を示す図3及び図4のフローチャートに基づいて
説明する。
【0024】まず、コンピュータ6から、電流供給部7
に対して、n個の発振コイルC1〜Cnに対し、順次所定の
時間間隔で所定強度の電流を出力するよう指示する(S
1)。
【0025】次に、発振コイル選択部4は、磁束計S1〜
Smでそれぞれ検出され、各発振コイルC1〜Cn個々に生じ
た磁場データをデータ収集ユニット2より入力し(S
2)、周知の最小二乗法などにより、磁束計S1〜Smに対
する各発振コイルC1〜Cn個々の相対位置を算出する(S
3)。
【0026】そして、発振コイル位置算出部3は、求め
た各発振コイルC1〜Cn個々の相対位置の正しさを示すgo
odness-of-fit (以下「GOF」という。)値を各発振
コイルC1〜Cn毎に求める(S4)。
【0027】ここで、GOFとは、ある時刻におけるす
べての磁気検出素子の磁場強度データを用いて最小二乗
法で解析した電流源位置がどれだけ信頼に値するかを示
す指標であり、次式で求められる。
【0028】goodness-of-fit 値=(1−Σ(Bmi−B
ci)2 /Σ(Bmi)2 )×100 Bci:解析した位置に所定の発振コイルがあると仮定し
た場合にi番目の磁気検出素子(磁束計の各素子)が検
出すると考えられる磁場強度 Bmi:i番目の磁気検出素子が実際に計測した磁場強度 各発振コイルC1〜Cn毎にGOF値が求められると、求め
たGOF値を用いてMRI画像との位置合わせに用いる
各発振コイルC1〜Cnが選択される(S5)。
【0029】ここで、各発振コイルC1〜Cnの選択動作を
図4に示すフローチャートに基づいて説明する。
【0030】まず、各発振コイルC1〜Cn毎にGOF値と
所定の基準値とを比較し、GOF値が基準値を超えてい
る発振コイルのみ選択する(S51)。
【0031】ここで、基準値は、MRI画像との位置合
わせを行うのに診断上支障のない状態を示す値である
が、GOF値は、信号強度とノイズ量に依存するため、
基準値は、測定の種類(聴覚誘発脳滋、体性感覚脳滋な
ど)や装置自身のシステムのノイズ量、また測定時の環
境ノイズ量によって異なる。
【0032】但し、過去の経験から特定の測定では信号
強度のおおよその値は既知となっており、例えば、体性
感覚脳滋(手首正中神経刺激)では85.0%、体性感
覚脳滋(足首後頚椎神経刺激)では80.0%、聴覚誘
発脳滋では、90.0%、さらに、視覚誘発脳滋(半視
覚、パターンリバーサル刺激)では、90.0%とな
る。
【0033】これにより、被検体Mの関心領域から離れ
た位置にあり、位置算出上誤差原因となる磁束計のデー
タが排除されるので、MRI画像との位置合わせをより
正確に行うことができる。
【0034】GOF値が基準値以上の発振コイルが選択
されると、選択された発振コイルの任意の3つの組合せ
を求め、各3つの発振コイルを結ぶ三角形の面積を求め
る(S52)。例えば、選択された発振コイル数がNの
場合、任意の3つの組合せ数は、 N3 で求められ、N
=5の場合、 53 =10となる。
【0035】任意の組合せ毎に三角形の面積が求められ
ると、それぞれの面積を比較し、三角形の面積が最大と
なる3つの発振コイルを選択する(S53)。
【0036】これにより、互いの距離が最も離れている
と考えられる3つの発振コイルが選択されるため、かか
る発振コイルの位置情報を用いることにより、最小限の
位置情報を用いてMRI画像との位置合わせをより正確
に行うことが可能となる。また、発振コイルは、通常髪
毛などによって不安定にならないよう、被検体の顔面に
付着されるため、選択する発振コイルが少ない場合、M
RI撮像時に被検体に付着されるコイル数が少なくな
り、被検体へ与える不快感を軽減することができる。
【0037】なお、三角形の面積Sは、ヘロンの公式よ
り次式で求めることができる。
【0038】 S=√(s(s−a)(s−b)(s−c)) a,b,c:それぞれ三角形の各片の長さ、 s:2
(a+b+c) 図1において、MRI画像との位置合わせに使用するた
めの発振コイルの選択が終了すると、次に、本来の目的
である被検体Mの生体磁気データの計測がなされ、磁場
解析部5は、得られたデータに基づいて生体活動電流源
の推定演算を行う(S6)。
【0039】そして、最終的に選択された3つの発振コ
イルの位置情報に基づき、直交座標変換等を行うこと
で、画像記憶部9から読み出したMRI画像との位置合
わせを行い、表示部8に求めた生体活動電流源をMRI
画像に重ねて表示する(S7)。
【0040】上述した実施形態では、各発振コイルの位
置算出誤差としてGOF値を用いたが、例えば、以下の
式で示される「測定磁場」と「解析位置からの理論磁
場」との相関係数γを用いても良い。
【0041】 γ=Σi (Bthi ×Bexi )/√(Σi (Bexi 2 )×
Σi (Bthi 2 )) Bexi :磁束計Si での測定磁場 Bthi :算出した位置から発生した時に、磁束計Si で
発生するはずの測定磁場 また、上述した実施形態では、3つの発振コイルを最終
的に選択する例を示したが、本発明はこれに限らず、例
えば、4つの発振コイルを選択する場合は、それらによ
り構成される三角錐の体積が最大となる組合せを選択
し、アフィン座標変換を用いてMRI画像などとの位置
合わせをすればよく、或いは、誤差が条件を満たすすべ
ての発振コイルに基づきMRI画像などとの位置合わせ
を行うようにしても良い。
【0042】これにより、4つの発振コイル位置に誤差
が同じだけ生じた場合でも、4つの位置が作る三角錐内
では空間対応付けの誤差が3つの位置を用いた場合に比
べて小さくなるというメリットが生じる。
【0043】さらに、上述した実施形態では、各発振コ
イルについてのGOF値や相関係数γと固定された基準
値とを比較したが、すべてのGOF値を高い順に並べ、
所定数の発振コイル、例えば上位5乃至10個を選択し
た上で、任意の三角形の組合せを求めても良い。
【0044】
【発明の効果】本発明によれば、各発振コイルの相対位
置の位置算出誤差が各発振コイルの位置算出に際し許容
範囲にあるもののみ選択されるため、複数の関心領域が
存在する場合であっても、磁束計を移動させることな
く、被検体の相対的な位置関係の高精度な特定が可能と
なる。
【図面の簡単な説明】
【図1】本発明にかかる生体磁気計測装置の一実施形態
を示す図である。
【図2】本発明にかかる発振コイルの一実施形態を示す
図である。
【図3】本発明の動作を示すフローチャートである。
【図4】本発明の動作を示すフローチャートである。
【符号の説明】
M 被検体 S1〜Sm 磁束計 C1〜Cn 発振コイル 1 センサーユニット 2 データ収集ユニット 6 コンピュータ 7 電流供給部

Claims (2)

    【特許請求の範囲】
  1. 【請求項1】 被検体内の生体活動電流源に伴って発生
    する微小磁界を計測し、その計測データに基づいて前記
    被検体内の生体活動電流源を求める生体磁気計測装置に
    おいて、 前記被検体に付着される複数の発振コイルと、 この複数の発振コイルに所定電流を供給する電流供給手
    段と、 各発振コイルから生じた磁場を計測する磁束計と、 各発振コイル毎に得られた磁場データに基づき前記磁束
    計に対する各発振コイルの相対位置を算出する磁場解析
    手段と、 得られた各発振コイルの相対位置からその位置算出誤差
    を求め、当該誤差が発振コイルの位置算出に際し許容範
    囲にあるものを選択するコイル選択手段と、 を備えることで、選択された発振コイルの位置情報から
    MRI等の被検体画像との位置合わせを行うよう構成し
    たことを特徴とする生体磁気計測装置。
  2. 【請求項2】 前記コイル選択手段は、前記位置算出誤
    差が許容範囲にある発振コイルのうち、任意の組合せに
    より得られる三角形の面積を算出し、当該面積が最大と
    なる3つの発振コイルを選択することを特徴とする請求
    項1記載の生体磁気計測装置。
JP01951897A 1997-01-31 1997-01-31 生体磁気計測装置 Expired - Fee Related JP3661329B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP01951897A JP3661329B2 (ja) 1997-01-31 1997-01-31 生体磁気計測装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP01951897A JP3661329B2 (ja) 1997-01-31 1997-01-31 生体磁気計測装置

Publications (2)

Publication Number Publication Date
JPH10211183A true JPH10211183A (ja) 1998-08-11
JP3661329B2 JP3661329B2 (ja) 2005-06-15

Family

ID=12001579

Family Applications (1)

Application Number Title Priority Date Filing Date
JP01951897A Expired - Fee Related JP3661329B2 (ja) 1997-01-31 1997-01-31 生体磁気計測装置

Country Status (1)

Country Link
JP (1) JP3661329B2 (ja)

Also Published As

Publication number Publication date
JP3661329B2 (ja) 2005-06-15

Similar Documents

Publication Publication Date Title
JP3473210B2 (ja) 生体磁気計測装置
US20020062076A1 (en) Apparatus for measuring bio-magnetic fields
JP3387236B2 (ja) 生体磁気計測装置
JP2751408B2 (ja) 脳磁計測装置
JP4791797B2 (ja) 生体磁場計測装置
Tsukada et al. Newly developed magnetocardiographic system for diagnosing heart disease
JP3661329B2 (ja) 生体磁気計測装置
JP3409551B2 (ja) 生体磁気計測装置
JPH04303416A (ja) 生体磁気計測装置
JP3591121B2 (ja) 生体磁気計測装置
JP3814923B2 (ja) 生体磁気計測装置
Adachi et al. A SQUID biomagnetometer system for measurement of spinal cord evoked magnetic fields
JP3233444B2 (ja) 生体磁気計測装置、生体磁気計測方法及び生体磁気計測装置用装着具
JP3298312B2 (ja) 生体活動電流源推定装置
JPH03251226A (ja) 生体磁気計測法
JP7002416B2 (ja) 磁界計測装置
JP3791119B2 (ja) 生体磁気計測装置
JP3651146B2 (ja) 生体磁気計測装置
JPH0435642A (ja) 生体磁気計測装置
JPH0555126B2 (ja)
JPH119567A (ja) 生体磁気計測装置
JP2000325322A (ja) 生体活動電流源推定装置
JP2795211B2 (ja) 生体磁気計測装置
JP2000325323A (ja) 生体活動電流源推定装置
JPH02180243A (ja) 医用計測装置

Legal Events

Date Code Title Description
TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20050301

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20050314

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080401

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090401

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100401

Year of fee payment: 5

LAPS Cancellation because of no payment of annual fees