JPH05337094A - 生体磁界の無雑音測定のためのシステムおよび方法 - Google Patents
生体磁界の無雑音測定のためのシステムおよび方法Info
- Publication number
- JPH05337094A JPH05337094A JP5041586A JP4158693A JPH05337094A JP H05337094 A JPH05337094 A JP H05337094A JP 5041586 A JP5041586 A JP 5041586A JP 4158693 A JP4158693 A JP 4158693A JP H05337094 A JPH05337094 A JP H05337094A
- Authority
- JP
- Japan
- Prior art keywords
- noise
- gradiometer
- reference system
- dipole
- biomagnetic
- 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.)
- Withdrawn
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/022—Measuring gradient
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring Magnetic Variables (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
Abstract
(57)【要約】
【目的】 生体磁界を測定している領域の外側からの磁
界による雑音を抑圧する。 【構成】 磁界勾配を表す第一の信号を発生する少なく
とも１個の一次センサ勾配計１００を生体磁気源１３０
の近くに配置し、干渉磁界およびその一次勾配を測定す
る基準系１２０を該勾配計の近くに配置し、勾配計およ
び基準系からの信号を組み合わせて、雑音を除去した信
号を発生する相殺手段１５０を設ける。
界による雑音を抑圧する。 【構成】 磁界勾配を表す第一の信号を発生する少なく
とも１個の一次センサ勾配計１００を生体磁気源１３０
の近くに配置し、干渉磁界およびその一次勾配を測定す
る基準系１２０を該勾配計の近くに配置し、勾配計およ
び基準系からの信号を組み合わせて、雑音を除去した信
号を発生する相殺手段１５０を設ける。
Description
【０００１】
【産業上の利用分野】本発明は、その中で生体磁界を測
定している領域の外側から発する磁界による雑音を抑圧
するためのシステムおよび方法に関するものである。
定している領域の外側から発する磁界による雑音を抑圧
するためのシステムおよび方法に関するものである。
【０００２】
【従来の技術】脳、心臓のような生体発生源からの磁界
は大きさが背景地球磁界に比べて約７桁から８桁小さ
い。このような小さい磁界の測定は超伝導量子干渉装置
（ＳＱＵＩＤ：Ｓｕｐｅｒｃｏｎｄｕｃｔｉｎｇ Ｑｕ
ａｎｔｕｍ Ｉｎｔｅｒｆｅｒｅｎｃｅ Ｄｅｖｉｃ
ｅ）およびそれに結合された電子回路を使用して行うこ
とができる。しかし、背景地磁気雑音の抑圧に多少注意
を払わないと、より小さな生体信号はこの背景地磁気雑
音によって圧倒されるのが普通である。背景地磁気雑音
の抑圧には、二つの方法がある。第一の方法では、磁気
干渉を排除するかまたは測定領域のまわりに磁気干渉を
「向ける」高導電性や高透磁率の材料の閉じた小室を使
用して、所望の発生源およびセンサを干渉磁界から遮へ
いする。この方法により良好な結果が得られ、商用の生
体磁気システムはこの技術を用いる。しかし、このよう
な遮へいされた小室の構成は、材料としても不動産とし
ても高価である。
は大きさが背景地球磁界に比べて約７桁から８桁小さ
い。このような小さい磁界の測定は超伝導量子干渉装置
（ＳＱＵＩＤ：Ｓｕｐｅｒｃｏｎｄｕｃｔｉｎｇ Ｑｕ
ａｎｔｕｍ Ｉｎｔｅｒｆｅｒｅｎｃｅ Ｄｅｖｉｃ
ｅ）およびそれに結合された電子回路を使用して行うこ
とができる。しかし、背景地磁気雑音の抑圧に多少注意
を払わないと、より小さな生体信号はこの背景地磁気雑
音によって圧倒されるのが普通である。背景地磁気雑音
の抑圧には、二つの方法がある。第一の方法では、磁気
干渉を排除するかまたは測定領域のまわりに磁気干渉を
「向ける」高導電性や高透磁率の材料の閉じた小室を使
用して、所望の発生源およびセンサを干渉磁界から遮へ
いする。この方法により良好な結果が得られ、商用の生
体磁気システムはこの技術を用いる。しかし、このよう
な遮へいされた小室の構成は、材料としても不動産とし
ても高価である。
【０００３】第二の方法では、所望の磁界に対する感度
に比べて干渉磁界に対する感度がずっと小さいセンサを
構成することにより干渉磁界を抑圧する。このようなセ
ンサは勾配計と呼ばれる。磁界のある高次の勾配（また
は空間導関数）に対して敏感であるように設計されるか
らである。この方式が利用している事実は、生体磁気源
に比べて磁気干渉源がセンサからずっと離れているの
で、それらの磁気勾配がより小さい（少なくとも勾配計
コイルの「有限差」により測定した場合）ということで
ある。双極子状の干渉源が１５０フィート以上離れてい
れば、二次の勾配計はその干渉を臨床品質の磁気遮へい
された室のほんの数分の一のコストで、それと同程度減
衰する。
に比べて干渉磁界に対する感度がずっと小さいセンサを
構成することにより干渉磁界を抑圧する。このようなセ
ンサは勾配計と呼ばれる。磁界のある高次の勾配（また
は空間導関数）に対して敏感であるように設計されるか
らである。この方式が利用している事実は、生体磁気源
に比べて磁気干渉源がセンサからずっと離れているの
で、それらの磁気勾配がより小さい（少なくとも勾配計
コイルの「有限差」により測定した場合）ということで
ある。双極子状の干渉源が１５０フィート以上離れてい
れば、二次の勾配計はその干渉を臨床品質の磁気遮へい
された室のほんの数分の一のコストで、それと同程度減
衰する。
【０００４】都合の悪いことに、完全な勾配計を作るこ
とはできない。線を巻いたコイルの特定の機械的な構成
を使用して勾配計が作られるからである。理想的な勾配
計の性能を得るために、アライメントが非常に重要であ
る。たとえば、「作られたときの」二次の勾配計は多分
一次および零次の磁界にも応答する。これは「共通モー
ド応答」と呼ばれる。共通モード応答を最小にするため
に超伝導トリムタブが勾配計のまわりに機械的に配置さ
れるが、これらの調整は長たらしくて、時間とともに変
化しやすい。機械的トリム調整のもう一つの問題は、次
世代の生体磁気計に見られるような多数の（２０個より
多い）勾配計があると、すぐに手に負えなくなるという
ことである。
とはできない。線を巻いたコイルの特定の機械的な構成
を使用して勾配計が作られるからである。理想的な勾配
計の性能を得るために、アライメントが非常に重要であ
る。たとえば、「作られたときの」二次の勾配計は多分
一次および零次の磁界にも応答する。これは「共通モー
ド応答」と呼ばれる。共通モード応答を最小にするため
に超伝導トリムタブが勾配計のまわりに機械的に配置さ
れるが、これらの調整は長たらしくて、時間とともに変
化しやすい。機械的トリム調整のもう一つの問題は、次
世代の生体磁気計に見られるような多数の（２０個より
多い）勾配計があると、すぐに手に負えなくなるという
ことである。
【０００５】したがって、平衡とトリムの操作に伴う困
難無しに遠くからの磁気干渉を更により良く抑圧するよ
うに高次の勾配計応答を達成するため、従来の勾配計を
改善しなければならない。
難無しに遠くからの磁気干渉を更により良く抑圧するよ
うに高次の勾配計応答を達成するため、従来の勾配計を
改善しなければならない。
【０００６】
【発明の目的】したがって、本発明の一つの目的は周囲
の発生源からの干渉無しに生体信号を測定するシステム
および方法を提供することである。本発明のもう一つの
目的は時間とともに変化する干渉磁界に応動するように
実時間で有効な、このようなシステムおよび方法を提供
することである。
の発生源からの干渉無しに生体信号を測定するシステム
および方法を提供することである。本発明のもう一つの
目的は時間とともに変化する干渉磁界に応動するように
実時間で有効な、このようなシステムおよび方法を提供
することである。
【０００７】更にもう一つの目的は組み立ての容易な比
較的簡単な構成要素を使用する、このようなシステムを
提供することである。
較的簡単な構成要素を使用する、このようなシステムを
提供することである。
【０００８】
【発明の概要】本発明の上記および他の目的は、一次セ
ンサ勾配計、該勾配計の近くに配置された基準系、勾配
計および基準系からの信号を相殺することにより雑音を
除去した信号を作成するための手段をそなえた生体磁気
測定用のシステムを設けることにより達成される。セン
サ位置での雑音磁界は、基準系の測定から決定される三
つの「等価」双極子を使用して推定される。事実上の勾
配計応答は、実際の勾配計応答から勾配計での雑音磁束
差を減算することにより計算される。この事実上の勾配
計応答は生体磁界に対する、雑音を除去した所望の応答
である。
ンサ勾配計、該勾配計の近くに配置された基準系、勾配
計および基準系からの信号を相殺することにより雑音を
除去した信号を作成するための手段をそなえた生体磁気
測定用のシステムを設けることにより達成される。セン
サ位置での雑音磁界は、基準系の測定から決定される三
つの「等価」双極子を使用して推定される。事実上の勾
配計応答は、実際の勾配計応答から勾配計での雑音磁束
差を減算することにより計算される。この事実上の勾配
計応答は生体磁界に対する、雑音を除去した所望の応答
である。
【０００９】本発明の新規な特徴は特許請求の範囲に明
確に記載されているが、本発明の構成および内容は本発
明の他の目的および特徴とともに、図面を参照した以下
の詳細な説明から、より良く理解され、評価される筈で
ある。
確に記載されているが、本発明の構成および内容は本発
明の他の目的および特徴とともに、図面を参照した以下
の詳細な説明から、より良く理解され、評価される筈で
ある。
【００１０】
【詳しい説明】図１は本発明の概略図である。主検知素
子１００は一次勾配計群である。一次勾配計群の基線ｄ
（コイル間の距離）は、生体磁気源１３０に最も近い検
知コイルと生体磁気源自体との間の距離より幾分か大き
くなるように選定される。生体磁気源１３０からの磁界
はｒ-2に比例するので、その生体磁気源に対する勾配計
の応答は勾配計より磁気計の応答に似たものとなる。た
とえば、勾配計の基線は５ｃｍであり、生体磁気源に対
する距離が２ｃｍである場合には、遠方コイルにより遮
断される磁束は近傍コイルによって遮断される磁束のほ
ぼ（２／７）2 すなわち１０％未満となる。磁気干渉源
が勾配計の基線に比べて大きい距離の所に配置されてい
るとすれば、勾配計の応答はほぼ、干渉磁界の一次勾配
に勾配計基線を乗算したものとなる。
子１００は一次勾配計群である。一次勾配計群の基線ｄ
（コイル間の距離）は、生体磁気源１３０に最も近い検
知コイルと生体磁気源自体との間の距離より幾分か大き
くなるように選定される。生体磁気源１３０からの磁界
はｒ-2に比例するので、その生体磁気源に対する勾配計
の応答は勾配計より磁気計の応答に似たものとなる。た
とえば、勾配計の基線は５ｃｍであり、生体磁気源に対
する距離が２ｃｍである場合には、遠方コイルにより遮
断される磁束は近傍コイルによって遮断される磁束のほ
ぼ（２／７）2 すなわち１０％未満となる。磁気干渉源
が勾配計の基線に比べて大きい距離の所に配置されてい
るとすれば、勾配計の応答はほぼ、干渉磁界の一次勾配
に勾配計基線を乗算したものとなる。
【００１１】
【数７】
【００１２】ここでＮはコイル当たりの巻数であり、Ａ
はコイルの面積である。殆どの干渉源は磁気双極子によ
って表すことができる。これらの磁気双極子は磁界強度
がｒ-3に比例しているので、その勾配がｒ-4に比例し、
より高次の勾配の大きさはより早く減衰する。センサ勾
配計１００より上、したがって生体磁気源から更に遠く
に、基準系１２０が配置されている。基準系は空間内の
一点に於ける「雑音」磁界の測定値およびそれの一次勾
配を供給する。磁界の三つのベクトル成分があるので、
３個の直交磁気計が必要とされる。勾配テンソルの９個
の成分があるが、測定点では磁界がソレノイド状（▽・
Ｂ＝０）であり、保存性（▽×ベクトルＢ＝０）である
ので、５個の勾配成分だけを測定すればよい。したがっ
て基準系１２０には、３個の磁気計および５個の一次勾
配計が含まれている。
はコイルの面積である。殆どの干渉源は磁気双極子によ
って表すことができる。これらの磁気双極子は磁界強度
がｒ-3に比例しているので、その勾配がｒ-4に比例し、
より高次の勾配の大きさはより早く減衰する。センサ勾
配計１００より上、したがって生体磁気源から更に遠く
に、基準系１２０が配置されている。基準系は空間内の
一点に於ける「雑音」磁界の測定値およびそれの一次勾
配を供給する。磁界の三つのベクトル成分があるので、
３個の直交磁気計が必要とされる。勾配テンソルの９個
の成分があるが、測定点では磁界がソレノイド状（▽・
Ｂ＝０）であり、保存性（▽×ベクトルＢ＝０）である
ので、５個の勾配成分だけを測定すればよい。したがっ
て基準系１２０には、３個の磁気計および５個の一次勾
配計が含まれている。
【００１３】基準系１２０は雑音源のモデルを作成する
ために必要とされるデータを供給する。詳しく述べる
と、モデルは空間内の３個の双極子、すなわちｘ−ｙ平
面内のｚ方向の双極子、ｚ−ｘ平面内のｙ方向の双極
子、およびｙ−ｚ平面内のｘ方向の双極子で構成され
る。これらの双極子は測定された磁界に対する「最良の
はめ込み」に過ぎず、物理的な双極子を表すものではな
いので、ここでは「等価」双極子と呼ぶ。
ために必要とされるデータを供給する。詳しく述べる
と、モデルは空間内の３個の双極子、すなわちｘ−ｙ平
面内のｚ方向の双極子、ｚ−ｘ平面内のｙ方向の双極
子、およびｙ−ｚ平面内のｘ方向の双極子で構成され
る。これらの双極子は測定された磁界に対する「最良の
はめ込み」に過ぎず、物理的な双極子を表すものではな
いので、ここでは「等価」双極子と呼ぶ。
【００１４】双極子の大きさＭ、および位置（ｘ，ｙ，
ｚ）のパラメータは下記のように基準系測定値から決め
ることができる。ｚ双極子の場合、
ｚ）のパラメータは下記のように基準系測定値から決め
ることができる。ｚ双極子の場合、
【００１５】
【数８】
【００１６】ｙ双極子の場合、
【００１７】
【数９】
【００１８】ｘ双極子の場合、
【００１９】
【数１０】
【００２０】Ｂxy，Ｂxz，Ｂyzは独立でないので、双極
子パラメータが矛盾無く決められるようにいくつかの付
加的な拘束が必要とされた。使用される拘束の組は発見
的な（ｈｅｕｒｉｓｔｉｃ）アルゴリズムに従って得ら
れた。ｘ−ｙ平面のｚ方向の双極子は磁界勾配
子パラメータが矛盾無く決められるようにいくつかの付
加的な拘束が必要とされた。使用される拘束の組は発見
的な（ｈｅｕｒｉｓｔｉｃ）アルゴリズムに従って得ら
れた。ｘ−ｙ平面のｚ方向の双極子は磁界勾配
【００２１】
【数１１】
【００２２】によって決定され、それらに寄与する。し
かし、ベクトルＢの回転の無い性質は、ｙ−ｚ平面のｘ
方向の双極子によって
かし、ベクトルＢの回転の無い性質は、ｙ−ｚ平面のｘ
方向の双極子によって
【００２３】
【数１２】
【００２４】も決まるということを意味する。同様なこ
とは
とは
【００２５】
【数１３】
【００２６】およびｚ−ｘ平面のｙ方向の双極子につい
ても言える。一般に言えることであるが、各回転導関数
には３個の等価双極子の中のある２個の等価双極子から
の寄与があり、それらを分離する独特の方法は無い。等
価双極子の位置を決める式に対する回転導関数の寄与
は、影響を受ける双極子が作る磁界の大きさによって重
み付けされる。たとえば、ｘ−ｙ平面のｚ方向の双極子
の場合、双極子決定式（２）−（４）で量Ｂzxが
ても言える。一般に言えることであるが、各回転導関数
には３個の等価双極子の中のある２個の等価双極子から
の寄与があり、それらを分離する独特の方法は無い。等
価双極子の位置を決める式に対する回転導関数の寄与
は、影響を受ける双極子が作る磁界の大きさによって重
み付けされる。たとえば、ｘ−ｙ平面のｚ方向の双極子
の場合、双極子決定式（２）−（４）で量Ｂzxが
【００２７】
【数１４】
【００２８】に置き換えられた。残りの回転導関数に対
しても同様の置き換えが行われた。式（２）−（４）を
使用してこれらの３個の等価双極子が決定されれば、主
検知勾配計での雑音磁界を推定することができる。等価
雑音双極子の大きさと位置が与えられれば、磁気双極子
の式を使って次のように検知勾配計での磁界を計算する
ことができる。
しても同様の置き換えが行われた。式（２）−（４）を
使用してこれらの３個の等価双極子が決定されれば、主
検知勾配計での雑音磁界を推定することができる。等価
雑音双極子の大きさと位置が与えられれば、磁気双極子
の式を使って次のように検知勾配計での磁界を計算する
ことができる。
【００２９】
【数１５】
【００３０】但し、ベクトルｒp およびベクトルｒq は
それぞれ検知コイルおよび等価双極子に対する位置ベク
トルである。検知コイルの寸法および方向を使用して勾
配計コイルの磁束を計算する方法は知られている。基準
系の発散導関数出力
それぞれ検知コイルおよび等価双極子に対する位置ベク
トルである。検知コイルの寸法および方向を使用して勾
配計コイルの磁束を計算する方法は知られている。基準
系の発散導関数出力
【００３１】
【数１６】
【００３２】は等価雑音双極子が供給する外挿された磁
界に対する簡単な線形補正として使用される。これらの
補正は次式に従って作成される。
界に対する簡単な線形補正として使用される。これらの
補正は次式に従って作成される。
【００３３】
【数１７】
【００３４】補正は双極子の外挿された磁界に付加され
る。図２に示すように、基準系２１０はｘ、ｙ、および
ｚを規定する原点（０，０，０）に配置されている。基
準系２１０の８個の成分で得られる測定値に基づいて、
３個の「等価」双極子２１２、２１４、および２１６が
計算される。３個の双極子は干渉磁界を生じる干渉源に
対する近似である。ベクトル２２０は式（５）のベクト
ルｒq であり、これは基準系２１０から一つの等価双極
子２１２に向かう位置ベクトルである。図２にはベクト
ル２３０も示されている。ベクトル２３０は式（５）の
ベクトルｒp である。ベクトル２３０は基準系２１０か
ら勾配計２４０の一つのコイルに向かう位置べクトルで
ある。式（５）は一つの双極子による一つのコイルでの
磁界を与えるので、両方の勾配計コイルに対する３個の
等価双極子の影響を明らかにするため、存在する各勾配
計に対してこの式を６回使用しなければならない。
る。図２に示すように、基準系２１０はｘ、ｙ、および
ｚを規定する原点（０，０，０）に配置されている。基
準系２１０の８個の成分で得られる測定値に基づいて、
３個の「等価」双極子２１２、２１４、および２１６が
計算される。３個の双極子は干渉磁界を生じる干渉源に
対する近似である。ベクトル２２０は式（５）のベクト
ルｒq であり、これは基準系２１０から一つの等価双極
子２１２に向かう位置ベクトルである。図２にはベクト
ル２３０も示されている。ベクトル２３０は式（５）の
ベクトルｒp である。ベクトル２３０は基準系２１０か
ら勾配計２４０の一つのコイルに向かう位置べクトルで
ある。式（５）は一つの双極子による一つのコイルでの
磁界を与えるので、両方の勾配計コイルに対する３個の
等価双極子の影響を明らかにするため、存在する各勾配
計に対してこの式を６回使用しなければならない。
【００３５】勾配計２４０の一つのコイルを通る磁束は
雑音磁束と（生体磁気源からの）所望の信号磁束との合
計である。前に説明したように、式（５）ならびにコイ
ルの面積と方向の知識から雑音磁束を計算することがで
きる。測定された合計磁束からこの雑音磁束を減算すれ
ば、所望の生体磁界だけによる磁束が結果として得られ
る。勾配計はそれの２個のコイルの間の磁束の差を測定
するので、勾配計２４０の事実上の勾配計出力（生体磁
気源だけに対して勾配計応答がどうなるか）は次のよう
に計算される。
雑音磁束と（生体磁気源からの）所望の信号磁束との合
計である。前に説明したように、式（５）ならびにコイ
ルの面積と方向の知識から雑音磁束を計算することがで
きる。測定された合計磁束からこの雑音磁束を減算すれ
ば、所望の生体磁界だけによる磁束が結果として得られ
る。勾配計はそれの２個のコイルの間の磁束の差を測定
するので、勾配計２４０の事実上の勾配計出力（生体磁
気源だけに対して勾配計応答がどうなるか）は次のよう
に計算される。
【００３６】 事実上の勾配計出力＝物理的な勾配計出力−［コイル２の雑音磁束−コイル１ の雑音磁束］ （６） 再び図１を参照して説明する。物理的な勾配計信号１１
０（式（６）の物理的な勾配計出力）が相殺手段１５０
に与えられる。相殺手段１５０は簡単な電子回路とする
ことができる。信号１６０は式（６）の括弧内の部分で
あり、これを信号１１０から減算することにより雑音を
除去した信号１４０（式（６）の「事実上の」勾配計出
力）が得られる。
０（式（６）の物理的な勾配計出力）が相殺手段１５０
に与えられる。相殺手段１５０は簡単な電子回路とする
ことができる。信号１６０は式（６）の括弧内の部分で
あり、これを信号１１０から減算することにより雑音を
除去した信号１４０（式（６）の「事実上の」勾配計出
力）が得られる。
【００３７】本発明の特定の実施例を図示し、説明して
きたが、熟練した当業者が変形および変更を考えつき得
ることは明らかである。したがって、本発明の趣旨と範
囲内に入るこのような変形および変更をすべて包含する
ように、特許請求の範囲を記述してある。
きたが、熟練した当業者が変形および変更を考えつき得
ることは明らかである。したがって、本発明の趣旨と範
囲内に入るこのような変形および変更をすべて包含する
ように、特許請求の範囲を記述してある。
【図１】本発明による生体磁気測定システムの概略構成
図である。
図である。
【図２】本発明により等価双極子を「はめ込んだ」、三
次元空間内の本発明の要素の概略構成図である。
次元空間内の本発明の要素の概略構成図である。
１００ センサ勾配計 １１０ 物理的な勾配計信号 １２０ 基準系 １３０ 生体磁気源 １４０ 雑音を除去した信号 １５０ 相殺手段 １６０ 雑音磁束差信号 ２１０ 基準系 ２４０ 勾配計
Claims (8)
- 【請求項１】 生体磁界の無雑音測定のための方法に於
いて、 磁界勾配を表す第一の信号を発生する少なくとも１個の
一次センサ勾配計を生体磁気源の近くに配置するステッ
プ、 干渉磁界およびその一次の勾配を測定する基準系を上記
勾配計の近くに配置するステップ、 上記基準系に対する各平面ｘ−ｙ、ｘ−ｚ、ｙ−ｚの磁
気双極子として上記干渉磁界を表すステップ、 上記基準系から上記勾配計コイルまでの位置ベクトルを
ベクトルｒp 、上記基準系から上記磁気双極子までの位
置ベクトルをベクトルｒq として、次の磁気双極子の式 【数１】 を使って、上記勾配計のコイルでの上記磁気干渉による
雑音の推定値を計算するステップ、 上記基準系で測定された一次勾配 【数２】 を使って、上記推定値に対して線形補正を行うステッ
プ、 上記勾配計のコイルに於ける雑音磁束を計算して、上記
勾配計に対する雑音磁束勾配を表す第二の信号を発生す
るステップ、および上記第一の信号と上記第二の信号の
差を求めることにより、無雑音勾配計出力を表す信号を
得るステップを含むことを特徴とする生体磁界の無雑音
測定方法。 - 【請求項２】 上記基準系が３個の磁気計および５個の
一次勾配計を含んでいる請求項１記載の生体磁界の無雑
音測定方法。 - 【請求項３】 上記干渉磁界を表す上記ステップが、次
式に従ってｘ−ｙ平面のｚ方向の双極子を計算するステ
ップ、 【数３】 次式に従ってｘ−ｚ平面のｙ方向の双極子を計算するス
テップ、 【数４】 次式に従ってｙ−ｚ平面のｘ方向の双極子を計算するス
テップ 【数５】 が含まれる請求項２記載の生体磁界の無雑音測定方法。 - 【請求項４】 上記線形補正ステップは次式 【数６】 に従って上記雑音推定値に対する線形補正を行う請求項
３記載の生体磁界の無雑音測定方法。 - 【請求項５】 生体磁界測定システムに於いて、 少なくとも１個の一次センサ勾配計、 ３個の磁気計および５個の一次勾配計を含む、上記勾配
計の近くに配置された基準系、および上記勾配計および
上記基準系の作成する信号を相殺することにより、雑音
を除去した信号を発生する相殺手段を含むことを特徴と
する生体磁界測定システム。 - 【請求項６】 上記相殺手段がアナログ回路である請求
項５記載の生体磁界測定システム。 - 【請求項７】 上記相殺手段がコンピュータソフトウェ
アである請求項５記載の生体磁界測定システム。 - 【請求項８】 上記相殺手段がディジタル回路である請
求項５記載の生体磁界測定システム。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US07/850,725 US5187436A (en) | 1992-03-13 | 1992-03-13 | Noise cancellation method in a biomagnetic measurement system using an extrapolated reference measurement |
| US850725 | 1992-03-13 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH05337094A true JPH05337094A (ja) | 1993-12-21 |
Family
ID=25308949
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP5041586A Withdrawn JPH05337094A (ja) | 1992-03-13 | 1993-03-03 | 生体磁界の無雑音測定のためのシステムおよび方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US5187436A (ja) |
| EP (1) | EP0560621A3 (ja) |
| JP (1) | JPH05337094A (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006520892A (ja) * | 2003-03-14 | 2006-09-14 | エレクタ ニューロマグ オイ | 磁場のマルチチャネル測定値を処理するための方法とシステム |
| JP2007522465A (ja) * | 2004-02-13 | 2007-08-09 | エレクタ アクチボラゲット（パブル） | 測定装置における干渉抑制方法 |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7092352B2 (en) * | 1993-07-23 | 2006-08-15 | Aquity, Llc | Cancellation systems for multicarrier transceiver arrays |
| DE4439691A1 (de) * | 1994-11-07 | 1996-05-09 | Philips Patentverwaltung | Verfahren zur Bestimmung der räumlichen Feldverteilung |
| US5657756A (en) * | 1995-06-07 | 1997-08-19 | Ctf Systems Inc. | Method and systems for obtaining higher order gradiometer measurements with lower order gradiometers |
| US5642045A (en) * | 1995-08-18 | 1997-06-24 | International Business Machines Corporation | Magnetic field gradiometer with improved correction circuits |
| DE19733920A1 (de) * | 1997-08-06 | 1999-02-25 | Forschungszentrum Juelich Gmbh | Vorrichtung zur Störsignalunterdrückung mit frequenzabhängiger adaptiver Kompensation sowie eine solche Vorrichtung enthaltendes Gradiometer |
| FI112399B (fi) | 1998-08-28 | 2003-11-28 | Neuromag Oy | Menetelmä ja laite taustahäiriön poistamiseksi monikanavaisista ilmaisinasetelmista |
| JP3379488B2 (ja) * | 1999-09-14 | 2003-02-24 | 株式会社日立製作所 | 磁場計測装置 |
| AU2002226874A1 (en) * | 2000-08-29 | 2002-04-08 | Cardiomag Imaging, Inc. | High balance gradiometer |
| KR100400755B1 (ko) * | 2001-11-08 | 2003-10-08 | 엘지전자 주식회사 | 보조센서에 의한 ｓｑｕｉｄ 센서 |
| FI124427B (fi) * | 2010-07-06 | 2014-08-29 | Elekta Ab | Menetelmä häiriöavaruuden tarkentamiseksi biomagneettisissa kenttämittauksissa |
| US9671214B2 (en) | 2013-07-17 | 2017-06-06 | Infineon Technologies Ag | Discrete magnetic angle sensor device, a magnetic angle sensor arrangement, a method for generating an angle signal and a method for providing a sensor signal |
| US10545197B2 (en) * | 2017-09-21 | 2020-01-28 | Apple Inc. | Self-calibrating system architecture for magnetic distortion compensation |
| EP4199806B1 (en) | 2020-08-31 | 2025-06-04 | Bard Access Systems, Inc. | Magnetic field direction detection |
| GB202213839D0 (en) * | 2022-09-22 | 2022-11-09 | Univ Court Univ Of Glasgow | Biomagnetidm measurement system and method for sensing biomagnetic signals |
| CN121091072B (zh) * | 2025-11-11 | 2026-01-30 | 国网甘肃省电力公司兰州供电公司 | 一种基于磁感应的压板状态监测方法及系统 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4386318A (en) * | 1980-09-26 | 1983-05-31 | Her Majesty The Queen In Right Of Canada, As Represented By The Minister Of National Defence | Method and apparatus to compensate a gradiometer having first and second unwanted terms |
| FR2558599B1 (fr) * | 1984-01-24 | 1986-08-01 | Thomson Csf | Procede et dispositif de compensation automatique des erreurs mecaniques d'un gradientmetre magnetique |
| DE4005079A1 (de) * | 1990-02-17 | 1991-08-22 | Dornier Gmbh | Vorrichtung zur vollstaendigen bestimmung des gradiententensors eines magnetfelds |
| DE4118126A1 (de) * | 1990-07-03 | 1992-01-09 | Siemens Ag | Verfahren zur aktiven unterdrueckung von magnetischen stoerfeldern bei einer anordnung zum messen biomagnetischer signale |
| US5122744A (en) * | 1990-10-09 | 1992-06-16 | Ibm Corporation | Gradiometer having a magnetometer which cancels background magnetic field from other magnetometers |
| GB2258314B (en) * | 1991-08-02 | 1995-07-19 | Dornier Gmbh | Device for measuring local magnetic field distributions |
-
1992
- 1992-03-13 US US07/850,725 patent/US5187436A/en not_active Expired - Fee Related
-
1993
- 1993-03-03 JP JP5041586A patent/JPH05337094A/ja not_active Withdrawn
- 1993-03-12 EP EP19930301885 patent/EP0560621A3/en not_active Withdrawn
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006520892A (ja) * | 2003-03-14 | 2006-09-14 | エレクタ ニューロマグ オイ | 磁場のマルチチャネル測定値を処理するための方法とシステム |
| JP2007522465A (ja) * | 2004-02-13 | 2007-08-09 | エレクタ アクチボラゲット（パブル） | 測定装置における干渉抑制方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| EP0560621A2 (en) | 1993-09-15 |
| US5187436A (en) | 1993-02-16 |
| EP0560621A3 (en) | 1994-08-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH05337094A (ja) | 生体磁界の無雑音測定のためのシステムおよび方法 | |
| EP0966689B1 (en) | Magnetic gradiometer | |
| Taulu et al. | Applications of the signal space separation method | |
| Drung | The PTB 83-SQUID system for biomagnetic applications in a clinic | |
| US8838225B2 (en) | Analysis of multi-channel measurement data using orthogonal virtual channels | |
| US5644229A (en) | Method of determining the spatial field distribution | |
| JP7222814B2 (ja) | 測定装置および測定方法 | |
| Taulu et al. | Novel noise reduction methods | |
| ter Brake et al. | Improvement of the performance of a mu-metal magnetically shielded room by means of active compensation (biomagnetic applications) | |
| JPH04215084A (ja) | 電界又は磁界の界パターンを測定する方法 | |
| JP7525297B2 (ja) | 磁場計測装置、磁場計測方法、および、磁場計測プログラム | |
| JP2007522465A (ja) | 測定装置における干渉抑制方法 | |
| US20120032679A1 (en) | Gradient coil arrangement | |
| JP3379488B2 (ja) | 磁場計測装置 | |
| JP2019010483A (ja) | 磁界計測装置および計測磁界表示方法 | |
| WO2023079081A1 (en) | Improvements in and relating to magnetic field nulling | |
| Hilgenfeld et al. | Active magnetic shielding for biomagnetic measurement using spatial gradient fields | |
| JP3655425B2 (ja) | 生体磁場計測装置 | |
| JP2003517734A (ja) | 磁気遮蔽室内の干渉を弱める方法 | |
| US20210247471A1 (en) | Optimal methods to feedback control and estimate magnetic fields to enable a neural detection system to measure magnetic fields from the brain | |
| JP7330014B2 (ja) | 心磁計測装置、心磁計測方法、および、心磁計測プログラム | |
| Kandori et al. | Cancellation technique of external noise inside a magnetically shielded room used for biomagnetic measurements | |
| Resmer et al. | Development of an active magnetic screen to allow a biomagnetometer to be used in an unshielded environment | |
| Garachtchenko et al. | Baseline distance optimization for SQUID gradiometers | |
| JP7636933B2 (ja) | 計測装置、計測方法、およびプログラム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A300 | Application deemed to be withdrawn because no request for examination was validly filed |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20000509 |