JPH04215084A - 電界又は磁界の界パターンを測定する方法 - Google Patents
電界又は磁界の界パターンを測定する方法Info
- Publication number
- JPH04215084A JPH04215084A JP3046128A JP4612891A JPH04215084A JP H04215084 A JPH04215084 A JP H04215084A JP 3046128 A JP3046128 A JP 3046128A JP 4612891 A JP4612891 A JP 4612891A JP H04215084 A JPH04215084 A JP H04215084A
- Authority
- JP
- Japan
- Prior art keywords
- vector
- field
- time
- activity
- current source
- 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
- 238000000034 method Methods 0.000 title claims description 8
- 230000005684 electric field Effects 0.000 title description 2
- 239000013598 vector Substances 0.000 claims abstract description 38
- 230000000694 effects Effects 0.000 claims description 41
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 230000002123 temporal effect Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 18
- 238000013459 approach Methods 0.000 abstract description 2
- 238000005303 weighing Methods 0.000 abstract 1
- 241000238366 Cephalopoda Species 0.000 description 4
- 230000000747 cardiac effect Effects 0.000 description 4
- 230000007423 decrease Effects 0.000 description 4
- 238000011835 investigation Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000011157 data evaluation Methods 0.000 description 1
- 210000002837 heart atrium Anatomy 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 230000030074 regulation of atrial cardiomyocyte membrane repolarization Effects 0.000 description 1
- 230000002336 repolarization Effects 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 238000011410 subtraction method Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 230000002861 ventricular Effects 0.000 description 1
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
- G01R33/025—Compensating stray fields
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)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【０００１】
【産業上の利用分野】この発明は、センサ装置と界を引
き起こす電流源の空間的位置を反復して算出するデータ
処理装置とにより、電界又は磁界の界パターンを測定す
る方法に関する。
き起こす電流源の空間的位置を反復して算出するデータ
処理装置とにより、電界又は磁界の界パターンを測定す
る方法に関する。
【０００２】
【従来の技術】地磁場に比べて１０の数乗弱い磁界の測
定は特に生体磁気活動に関連して重要である。その際生
体組織内部の神経電流から出る外部界を相応のセンサに
より種々の位置で測定し、そのようにして測定された界
パターンから相応のデータ評価により等価なダイポール
の位置を算出することが問題となる。そのために必要な
測定装置は極めて高感度でなければならず、そのために
特殊な磁界勾配計が用いられ、これらの磁界勾配計が超
伝導量子干渉計（スクイド）と共に超伝導状態のもとに
ジュワー容器中に収容されている。こうして信頼のおけ
る測定を可能にするためには、特に外部妨害界をできる
だけ十分に遮蔽しなければならない。このために患者及
び測定装置は出入り可能な磁気遮蔽室中に収容される。 磁界分布をできるだけ速やかに求めることができるよう
に、例えば専門誌「エレクトロメディカ（Ｅｌｅｃｔｒ
ｏｍｅｄｉｃａ ）」第５７巻、１９８９年、第２〜７
ページに記載されているように多チャネル装置が開発さ
れている。
定は特に生体磁気活動に関連して重要である。その際生
体組織内部の神経電流から出る外部界を相応のセンサに
より種々の位置で測定し、そのようにして測定された界
パターンから相応のデータ評価により等価なダイポール
の位置を算出することが問題となる。そのために必要な
測定装置は極めて高感度でなければならず、そのために
特殊な磁界勾配計が用いられ、これらの磁界勾配計が超
伝導量子干渉計（スクイド）と共に超伝導状態のもとに
ジュワー容器中に収容されている。こうして信頼のおけ
る測定を可能にするためには、特に外部妨害界をできる
だけ十分に遮蔽しなければならない。このために患者及
び測定装置は出入り可能な磁気遮蔽室中に収容される。 磁界分布をできるだけ速やかに求めることができるよう
に、例えば専門誌「エレクトロメディカ（Ｅｌｅｃｔｒ
ｏｍｅｄｉｃａ ）」第５７巻、１９８９年、第２〜７
ページに記載されているように多チャネル装置が開発さ
れている。
【０００３】しかし外から作用する妨害界だけが測定結
果を歪曲させるわけではない。おそらく多くの場合外部
磁界を引き起こしている調査中の生体電流源は隣接する
源の活動により付加重畳される。このバックグラウンド
活動は調査中の電流源の位置決めの際に誤差を招くおそ
れがある。このような事例は特に生体磁気的心臓活動（
心磁図）の場合に生じる。この種の重畳された界が十分
に知られているときには、調査される電流源の界を得る
ために、重畳された界を測定された界から差し引くこと
ができる。しかしこの種の一定の活動は特に心臓の刺激
電流の場合には与えられない。すなわち既知の時間帯中
のこの領域における電気生理的活動は、この時間帯での
空間的な界パターンがほぼ一定にとどまりしかし振幅係
数に関しては時間に関係して変化する例えばゆっくり増
減するような電界又は磁界を発生させる。この振幅係数
の時間的変化の推移は知られていない。この種の変化は
例えば、電気生理的活動が位置及び方向について安定し
ておりしかしイオン流の電流密度がゆっくり減少すると
きに生じる。例えばこのことは心房の再分極の場合に心
臓サイクルのＰ／Ｑ間隔で起こる。この活動が所定の時
間帯（時間窓）の開始後に位置決めしようとする別の活
動の一定の界パターンと重なるときは、このバックグラ
ウンド界パターンを位置決めしようとする活動の開始前
に所定の時点で測定し、もしこの活動が一定にとどまっ
ているときは、測定しようとする界パターンから差し引
くことができる。しかし実際に起こるように、バックグ
ラウンド活動が測定しようとする活動を突き止めるため
の測定サイクルの期間中に一定にとどまっていないなら
ば、この単純な減算法は機能しない。
果を歪曲させるわけではない。おそらく多くの場合外部
磁界を引き起こしている調査中の生体電流源は隣接する
源の活動により付加重畳される。このバックグラウンド
活動は調査中の電流源の位置決めの際に誤差を招くおそ
れがある。このような事例は特に生体磁気的心臓活動（
心磁図）の場合に生じる。この種の重畳された界が十分
に知られているときには、調査される電流源の界を得る
ために、重畳された界を測定された界から差し引くこと
ができる。しかしこの種の一定の活動は特に心臓の刺激
電流の場合には与えられない。すなわち既知の時間帯中
のこの領域における電気生理的活動は、この時間帯での
空間的な界パターンがほぼ一定にとどまりしかし振幅係
数に関しては時間に関係して変化する例えばゆっくり増
減するような電界又は磁界を発生させる。この振幅係数
の時間的変化の推移は知られていない。この種の変化は
例えば、電気生理的活動が位置及び方向について安定し
ておりしかしイオン流の電流密度がゆっくり減少すると
きに生じる。例えばこのことは心房の再分極の場合に心
臓サイクルのＰ／Ｑ間隔で起こる。この活動が所定の時
間帯（時間窓）の開始後に位置決めしようとする別の活
動の一定の界パターンと重なるときは、このバックグラ
ウンド界パターンを位置決めしようとする活動の開始前
に所定の時点で測定し、もしこの活動が一定にとどまっ
ているときは、測定しようとする界パターンから差し引
くことができる。しかし実際に起こるように、バックグ
ラウンド活動が測定しようとする活動を突き止めるため
の測定サイクルの期間中に一定にとどまっていないなら
ば、この単純な減算法は機能しない。
【０００４】
【発明が解決しようとする課題】この発明の課題は、界
パターンが測定しようとする活動の開始前に所定の時点
で測定できるような生理的又は技術的な発生源のバック
グラウンド活動を、測定期間中のその振幅推移を考慮し
ながら、各測定時点で適切な量だけ全体量から差し引く
ことができるように規定することにある。
パターンが測定しようとする活動の開始前に所定の時点
で測定できるような生理的又は技術的な発生源のバック
グラウンド活動を、測定期間中のその振幅推移を考慮し
ながら、各測定時点で適切な量だけ全体量から差し引く
ことができるように規定することにある。
【０００５】
【課題を解決するための手段】この課題はこの発明に基
づき、バックグラウンド活動を補償するために、突き止
めようとする電流源の活動の開始前の時点ｔ１ とｔ２
との間の時点ｔ０ で測定された値から、バックグラ
ウンド活動の結果を示す信号ベクトルが形成され、この
信号ベクトルが時点ｔ２ とｔ３ との間の時点ｔｘ
で測定された値のそれぞれ結果を示すベクトルからベク
トル的に重みづけして差し引かれ、それにより主として
突き止めようとする電流源の電流だけにより引き起こさ
れる信号ベクトルが形成され、その際バックグラウンド
界のベクトルの長さを変更する重みづけ係数が、調査さ
れる源パラメータと共に、調査される電流源のパラメー
タが最適に表示できるように、誤差最小化法により反復
して算定されることにより解決される。
づき、バックグラウンド活動を補償するために、突き止
めようとする電流源の活動の開始前の時点ｔ１ とｔ２
との間の時点ｔ０ で測定された値から、バックグラ
ウンド活動の結果を示す信号ベクトルが形成され、この
信号ベクトルが時点ｔ２ とｔ３ との間の時点ｔｘ
で測定された値のそれぞれ結果を示すベクトルからベク
トル的に重みづけして差し引かれ、それにより主として
突き止めようとする電流源の電流だけにより引き起こさ
れる信号ベクトルが形成され、その際バックグラウンド
界のベクトルの長さを変更する重みづけ係数が、調査さ
れる源パラメータと共に、調査される電流源のパラメー
タが最適に表示できるように、誤差最小化法により反復
して算定されることにより解決される。
【０００６】
【作用効果】この発明によれば、ベクトルの振幅減衰を
時間に関係して個々の測定点で算出する重みづけ係数を
バックグラウンド界のベクトルに割り付けることにより
、測定された界と測定しようとする電流源から出る界と
の間の界の差、従って後者の界自体のパターンを最善に
測定することができるということが達成される。
時間に関係して個々の測定点で算出する重みづけ係数を
バックグラウンド界のベクトルに割り付けることにより
、測定された界と測定しようとする電流源から出る界と
の間の界の差、従って後者の界自体のパターンを最善に
測定することができるということが達成される。
【０００７】
【実施例】次に心磁図測定装置とこの発明に基づく方法
による測定に関する線図との一実施例を示す図面により
、この発明を詳細に説明する。
による測定に関する線図との一実施例を示す図面により
、この発明を詳細に説明する。
【０００８】図１に示す装置は、患者２が横たわる患者
横臥台１及び患者の上方に配置されたジュワー容器４か
ら成り、このジュワー容器は多チャネル磁界勾配計装置
及び磁界勾配計に付設された直流スクイドを内蔵する。 患者横臥台１及びジュワー容器４は外部磁界が作用しな
いように遮蔽室７中に収容されている。電子装置８はス
クイド／磁界勾配計５の制御と給電とに役立つ。等価な
ダイポール１０から放出され患者の身体を囲む磁界１１
の種種の測定位置での時間的推移を表しスクイド／磁界
勾配計により求められた測定値の信号評価は、所定の反
復計算法に基づきコンピュータ９中で行われる。
横臥台１及び患者の上方に配置されたジュワー容器４か
ら成り、このジュワー容器は多チャネル磁界勾配計装置
及び磁界勾配計に付設された直流スクイドを内蔵する。 患者横臥台１及びジュワー容器４は外部磁界が作用しな
いように遮蔽室７中に収容されている。電子装置８はス
クイド／磁界勾配計５の制御と給電とに役立つ。等価な
ダイポール１０から放出され患者の身体を囲む磁界１１
の種種の測定位置での時間的推移を表しスクイド／磁界
勾配計により求められた測定値の信号評価は、所定の反
復計算法に基づきコンピュータ９中で行われる。
【０００９】さて測定される磁界パターンはしばしば、
空間的に位置決め可能な単一の電流源により発生させら
れるとは限らない。むしろ電流源を位置決めしようとす
る活動が例えば先駆活動から由来する可能性のある別の
活動に重畳し、先駆活動が比較的大きい組織複合体にわ
たり拡散し増大又は減衰する振幅をもって伝わり、位置
決めしようとする電流源の界に付加重畳するということ
が起こる。この種のバックグラウンド界は、評価が電流
源の誤った位置決めを招くように関心のある界パターン
に作用するおそれが大きい。
空間的に位置決め可能な単一の電流源により発生させら
れるとは限らない。むしろ電流源を位置決めしようとす
る活動が例えば先駆活動から由来する可能性のある別の
活動に重畳し、先駆活動が比較的大きい組織複合体にわ
たり拡散し増大又は減衰する振幅をもって伝わり、位置
決めしようとする電流源の界に付加重畳するということ
が起こる。この種のバックグラウンド界は、評価が電流
源の誤った位置決めを招くように関心のある界パターン
に作用するおそれが大きい。
【００１０】図２には、心臓活動に由来する外部磁界（
心磁図）の磁界勾配計のうちの一つの位置で測定された
時間的推移が示されている。ここで曲線の一部（時間窓
）だけを考察しよう。時間窓は時点ｔ１ と共に始まり
時点ｔ３ で終わる。時点ｔ１ で減衰しつつある活動
Ｓ１ は心臓サイクルのＰ／Ｑ間隔での心房の再分極に
由来する。これは磁界パターンがこの時間帯でほぼ一定
にとどまりしかし振幅係数については変化するすなわち
ゆっくり増減するような活動である。振幅係数のこの時
間的変化は知られていない。この変化は例えば、電気生
理的活動が位置及び方向において安定しておりしかしイ
オン流の電流密度がゆっくりと減少するときに生じる。 この活動Ｓ１ の振幅は相応に減衰する。時点ｔ２ で
、位置決めしようとする心室活動Ｓ２ が始まる。両活
動は窓期間全体にわたり重畳する。その際時点ｔでセン
サにより測定される信号はベクトル表示で次のように記
述することができる。 ‖ＳＮ（Ｍ）（ｔ）＝ＳＮ（Ｑ）（ｔ）＋ＳＮ（Ｂ）（
ｔ）＋ＳＮ（Ｒ）（ｔ）‖これに起因する位置決め誤差
を避けるために、活動Ｓ１ を活動Ｓ２ から差し引か
なければならない。このためにまず活動Ｓ１ はまだこ
の活動だけが発生している時点ｔ０ で測定される。時
間にわたるこの活動の振幅は一定でないので、次式に基
づく減算すなわち、‖ＳＮ（Ｍ）（ｔｘ）−［ＳＮ（Ｑ
ｉ）（ｔｘ）＋ＳＮ（Ｂ）（ｔ０）＋ＳＮ（Ｒ）（ｔ）
］‖＝０は目的を達成しない。おそらくこのためには、
バックグラウンド信号のベクトルＳＮ （Ｂ）にそれぞ
れの測定時点ｔｘ で掛けなければならない重みづけ係
数Ｗの導入が必要である。そして式はベクトルＳＮ （
Ｒ）（ｔ）を無視して次のようになる。 ‖ＳＮ（Ｍ）（ｔｘ）−［ＳＮ（Ｑｉ）（ｔｘ）＋Ｗ（
ｔｘ）×ＳＮ（Ｂ）（ｔ０）］‖＝ 最小これらの式に
おいて、 ＳＮ （Ｍ）（ｔ）＝時点ｔでのＮ個のセンサの測定信
号Ｓ１ ＋Ｓ２ のベクトルＳＮ （Ｑ）（ｔ）＝時点
ｔでのＮ個のセンサの全信号における調査中の電流源の
信号Ｓ２ のベ
クトルＳＮ （Ｂ）（ｔ）＝時点ｔでのＮ個のセンサの
バックグラウンド活動Ｓ１ のベク
トルＳＮ （Ｒ）（ｔ）＝時点ｔで
のＮ個のセンサのノイズ成分のベクトルＷ（ｔ）
＝ベクトルＳＮ （Ｂ）（ｔ）の長さに対す
る重みづけ係数Ｑｉ （ｔｘ ） ＝電流源及び
従って時点ｔｘ でのその信号を一義的に決定する
パラメータである。 この式は、バックグラウンド界が固有のバックグラウン
ド活動ＳＮ （Ｂ）と全般的なノイズ成分ＳＮ （Ｒ）
とを加えて合成されているので、ゼロにはならない。
心磁図）の磁界勾配計のうちの一つの位置で測定された
時間的推移が示されている。ここで曲線の一部（時間窓
）だけを考察しよう。時間窓は時点ｔ１ と共に始まり
時点ｔ３ で終わる。時点ｔ１ で減衰しつつある活動
Ｓ１ は心臓サイクルのＰ／Ｑ間隔での心房の再分極に
由来する。これは磁界パターンがこの時間帯でほぼ一定
にとどまりしかし振幅係数については変化するすなわち
ゆっくり増減するような活動である。振幅係数のこの時
間的変化は知られていない。この変化は例えば、電気生
理的活動が位置及び方向において安定しておりしかしイ
オン流の電流密度がゆっくりと減少するときに生じる。 この活動Ｓ１ の振幅は相応に減衰する。時点ｔ２ で
、位置決めしようとする心室活動Ｓ２ が始まる。両活
動は窓期間全体にわたり重畳する。その際時点ｔでセン
サにより測定される信号はベクトル表示で次のように記
述することができる。 ‖ＳＮ（Ｍ）（ｔ）＝ＳＮ（Ｑ）（ｔ）＋ＳＮ（Ｂ）（
ｔ）＋ＳＮ（Ｒ）（ｔ）‖これに起因する位置決め誤差
を避けるために、活動Ｓ１ を活動Ｓ２ から差し引か
なければならない。このためにまず活動Ｓ１ はまだこ
の活動だけが発生している時点ｔ０ で測定される。時
間にわたるこの活動の振幅は一定でないので、次式に基
づく減算すなわち、‖ＳＮ（Ｍ）（ｔｘ）−［ＳＮ（Ｑ
ｉ）（ｔｘ）＋ＳＮ（Ｂ）（ｔ０）＋ＳＮ（Ｒ）（ｔ）
］‖＝０は目的を達成しない。おそらくこのためには、
バックグラウンド信号のベクトルＳＮ （Ｂ）にそれぞ
れの測定時点ｔｘ で掛けなければならない重みづけ係
数Ｗの導入が必要である。そして式はベクトルＳＮ （
Ｒ）（ｔ）を無視して次のようになる。 ‖ＳＮ（Ｍ）（ｔｘ）−［ＳＮ（Ｑｉ）（ｔｘ）＋Ｗ（
ｔｘ）×ＳＮ（Ｂ）（ｔ０）］‖＝ 最小これらの式に
おいて、 ＳＮ （Ｍ）（ｔ）＝時点ｔでのＮ個のセンサの測定信
号Ｓ１ ＋Ｓ２ のベクトルＳＮ （Ｑ）（ｔ）＝時点
ｔでのＮ個のセンサの全信号における調査中の電流源の
信号Ｓ２ のベ
クトルＳＮ （Ｂ）（ｔ）＝時点ｔでのＮ個のセンサの
バックグラウンド活動Ｓ１ のベク
トルＳＮ （Ｒ）（ｔ）＝時点ｔで
のＮ個のセンサのノイズ成分のベクトルＷ（ｔ）
＝ベクトルＳＮ （Ｂ）（ｔ）の長さに対す
る重みづけ係数Ｑｉ （ｔｘ ） ＝電流源及び
従って時点ｔｘ でのその信号を一義的に決定する
パラメータである。 この式は、バックグラウンド界が固有のバックグラウン
ド活動ＳＮ （Ｂ）と全般的なノイズ成分ＳＮ （Ｒ）
とを加えて合成されているので、ゼロにはならない。
【００１１】この過程は図３に示されている。測定値ベ
クトルＳＮ （Ｍ）は測定値部分空間ＭＲ１ 中の一点
を示す。すべての測定値空間は図示の三つのセンサＧ１
、Ｇ２ 、Ｇ３ のそれぞれの信号により展開されて
次元Ｎを有し、ここでＮはセンサの数を表す。実際はセ
ンサの数ははるかに多く例えば３７である。この所定の
モデル電流源（例えば電流ダイポール）の信号ＳＮ （
Ｑ）はすべての源パラメータＱｉ の変化のもとに、測
定値部分空間中の例として十字印により示され離散した
点だけを描くことができる。探究戦略は電流源のすべて
の起こり得る信号ベクトルＳＮ （Ｑ）のうちの実際の
測定値ＳＮ （Ｍ）に最も近い信号ベクトルを求めるこ
とである。もし電流源に対する解をバックグラウンド界
ＳＮ （Ｂ）に関する補正無しに求めると、確かに測定
値部分空間ＭＲ１ での解が得られる。しかしながらこ
の解は誤差を含む源測定位置を招く。バックグラウンド
界ＳＮ （Ｂ）の減算により、新しい補正された測定信
号ＳＮ （Ｍ）−ＳＮ （Ｂ）の付近の解従って測定値
部分空間ＭＲ２ 中の解を求めるように探究戦略が導か
れる。その際補正ベクトルＳＮ （Ｂ）の長さに対する
重みづけ係数Ｗは、合成ベクトルＷＳＮ （Ｂ）が正し
い解にできるだけ近づくように反復して最適化される。
クトルＳＮ （Ｍ）は測定値部分空間ＭＲ１ 中の一点
を示す。すべての測定値空間は図示の三つのセンサＧ１
、Ｇ２ 、Ｇ３ のそれぞれの信号により展開されて
次元Ｎを有し、ここでＮはセンサの数を表す。実際はセ
ンサの数ははるかに多く例えば３７である。この所定の
モデル電流源（例えば電流ダイポール）の信号ＳＮ （
Ｑ）はすべての源パラメータＱｉ の変化のもとに、測
定値部分空間中の例として十字印により示され離散した
点だけを描くことができる。探究戦略は電流源のすべて
の起こり得る信号ベクトルＳＮ （Ｑ）のうちの実際の
測定値ＳＮ （Ｍ）に最も近い信号ベクトルを求めるこ
とである。もし電流源に対する解をバックグラウンド界
ＳＮ （Ｂ）に関する補正無しに求めると、確かに測定
値部分空間ＭＲ１ での解が得られる。しかしながらこ
の解は誤差を含む源測定位置を招く。バックグラウンド
界ＳＮ （Ｂ）の減算により、新しい補正された測定信
号ＳＮ （Ｍ）−ＳＮ （Ｂ）の付近の解従って測定値
部分空間ＭＲ２ 中の解を求めるように探究戦略が導か
れる。その際補正ベクトルＳＮ （Ｂ）の長さに対する
重みづけ係数Ｗは、合成ベクトルＷＳＮ （Ｂ）が正し
い解にできるだけ近づくように反復して最適化される。
【図１】心磁図測定装置の一例を示す斜視図である。
【図２】図１に示すスクイド／磁界勾配計の一つから得
られた信号の時間的推移を示すグラフである。
られた信号の時間的推移を示すグラフである。
【図３】この発明に基づく方法による測定値信号、源信
号及びバックグラウンド信号の合成の一実施例を幾何学
的に示すベクトル線図である。
号及びバックグラウンド信号の合成の一実施例を幾何学
的に示すベクトル線図である。
１０ 電流源
Ｂ バックグラウンド活動
Ｑｉ 源パラメータ
ＳＮ （Ｂ）（ｔ） バックグラウンド界の信号ベク
トルＳＮ （Ｍ）（ｔ） 測定値の結果を示す信号ベ
クトルＳＮ （Ｑ）（ｔ） 突き止めようとする電流
源の界の信号ベクトル Ｗ 重みづけ係数
トルＳＮ （Ｍ）（ｔ） 測定値の結果を示す信号ベ
クトルＳＮ （Ｑ）（ｔ） 突き止めようとする電流
源の界の信号ベクトル Ｗ 重みづけ係数
Claims (1)
- 【請求項１】 センサ装置と、突き止めようとする電
流源により形成されかつバックグラウンド活動と合成さ
れた界の時間的振幅推移及び空間的分布に基づき、界を
発生させる電流源の空間的位置を反復して算出するため
のデータ処理装置とにより、電界又は磁界の界パターン
を測定する方法であって、二つの時点ｔ１ とｔ３ と
の間の所定の時間窓内の所定の時点ｔｘ で得られた測
定値がベクトルによる測定値空間を形成するものにおい
て、バックグラウンド活動（Ｂ）を補償するために、突
き止めようとする電流源（１０）の活動の開始前の時点
ｔ１ とｔ２ との間の時点ｔ０ で測定された値から
、バックグラウンド活動（Ｂ）の結果を示す信号ベクト
ル（ＳＮ （Ｂ）（ｔ０ ））が形成され、この信号ベ
クトルが時点ｔ２ とｔ３ との間の時点ｔｘ で測定
された値のそれぞれ結果を示すベクトル（ＳＮ （Ｍ）
（ｔＸ ））からベクトル的に重みづけして差し引かれ
、それにより主として突き止めようとする電流源（１０
）の電流だけにより引き起こされる信号ベクトル（ＳＮ
（Ｑ））が形成され、その際バックグラウンド界（Ｂ）
のベクトル（ＳＮ （Ｂ）（ｔ０ ））の長さを変更す
る重みづけ係数（Ｗ）が、調査される源パラメータ（Ｑ
ｉ ）と共に、調査される電流源（１０）のパラメータ
が最適に表示できるように、誤差最小化法により反復し
て算定されることを特徴とする電界又は磁界の界パター
ンを測定する方法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP90103479A EP0443069A1 (de) | 1990-02-22 | 1990-02-22 | Verfahren zur Messung des Feldmusters elektrischer oder magnetischer Felder mit Hilfe einer Sensoranordnung |
| EP90103479.3 | 1990-02-22 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH04215084A true JPH04215084A (ja) | 1992-08-05 |
Family
ID=8203682
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP3046128A Pending JPH04215084A (ja) | 1990-02-22 | 1991-02-18 | 電界又は磁界の界パターンを測定する方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US5136242A (ja) |
| EP (1) | EP0443069A1 (ja) |
| JP (1) | JPH04215084A (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006520892A (ja) * | 2003-03-14 | 2006-09-14 | エレクタ ニューロマグ オイ | 磁場のマルチチャネル測定値を処理するための方法とシステム |
| JP2008219455A (ja) * | 2007-03-05 | 2008-09-18 | Oki Electric Ind Co Ltd | 等化器及び等化方法 |
Families Citing this family (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0614899A (ja) * | 1991-11-06 | 1994-01-25 | Mitsui Mining & Smelting Co Ltd | 脳磁界計測装置 |
| JP3204542B2 (ja) * | 1992-07-24 | 2001-09-04 | 株式会社東芝 | 磁場源測定装置 |
| US5311867A (en) * | 1993-01-04 | 1994-05-17 | Biomagnetic Technologies, Inc. | Detection and grouping analysis of cardiac cycles |
| JP2739804B2 (ja) * | 1993-05-14 | 1998-04-15 | 日本電気株式会社 | 双極子推定装置 |
| US5526811A (en) * | 1993-06-15 | 1996-06-18 | Biomagnetic Technologies, Inc. | Apparatus and process for determining the sources of biomagnetic activity |
| US5557199A (en) * | 1994-04-29 | 1996-09-17 | The United States Of America As Represented By The Department Of Health And Human Services | Magnetic resonance monitor |
| US5638301A (en) * | 1994-06-02 | 1997-06-10 | Ford Motor Company | Method and system for inspecting die sets using free-form inspection techniques |
| US5546328A (en) * | 1994-06-02 | 1996-08-13 | Ford Motor Company | Method and system for automated alignment of free-form geometries |
| JP3519141B2 (ja) * | 1994-10-31 | 2004-04-12 | 株式会社東芝 | 生体内電流源推定方法 |
| DE4439691A1 (de) * | 1994-11-07 | 1996-05-09 | Philips Patentverwaltung | Verfahren zur Bestimmung der räumlichen Feldverteilung |
| US6697660B1 (en) | 1998-01-23 | 2004-02-24 | Ctf Systems, Inc. | Method for functional brain imaging from magnetoencephalographic data by estimation of source signal-to-noise ratio |
| AU2043699A (en) | 1998-01-23 | 1999-08-09 | Ctf Systems Inc. | System and method for measuring, estimating and displaying rms current density maps |
| DE19860037C2 (de) * | 1998-12-23 | 2001-01-25 | Siemens Ag | Verfahren zur ortsaufgelösten Messung der elektrischen Aktivität von Nervenzellen mittels magnetischer Resonanz und Vorrichtung zur Durchführung des Verfahrens |
| JP3712350B2 (ja) * | 2000-07-18 | 2005-11-02 | 独立行政法人科学技術振興機構 | 心室遅延電位の心臓磁界診断装置およびその作動方法 |
| US8015977B2 (en) | 2003-10-31 | 2011-09-13 | Medtronic, Inc. | Indicator tool for use with an implantable medical device |
| WO2011010243A1 (en) * | 2009-07-20 | 2011-01-27 | Koninklijke Philips Electronics N.V. | Apparatus and method for influencing and/or detecting magnetic particles |
| FI124427B (fi) | 2010-07-06 | 2014-08-29 | Elekta Ab | Menetelmä häiriöavaruuden tarkentamiseksi biomagneettisissa kenttämittauksissa |
| US11372128B2 (en) | 2020-05-14 | 2022-06-28 | China Petroleum & Chemical Corporation | Method and system for detecting downhole magnetic interference on measurement while drilling operations |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4739262A (en) * | 1982-08-26 | 1988-04-19 | Mobil Oil Corporation | Method for removal of bias difference between magnetometers |
| DE3247585A1 (de) * | 1982-12-22 | 1984-06-28 | Siemens AG, 1000 Berlin und 8000 München | Mehrkanalige vorrichtung zur messung von verschiedenen feldquellen hervorgerufener schwacher magnetfelder |
| DE3324208A1 (de) * | 1983-07-05 | 1985-01-17 | Siemens AG, 1000 Berlin und 8000 München | Supraleitendes gradiometerspulensystem fuer eine vorrichtung zur mehrkanaligen messung schwacher, sich aendernder magnetfelder |
| DE3515237A1 (de) * | 1985-04-26 | 1986-10-30 | Siemens AG, 1000 Berlin und 8000 München | Vorrichtung zur messung schwacher magnetfelder mit wenigstens einem dc-squid |
| JPS63246688A (ja) * | 1987-04-01 | 1988-10-13 | Hitachi Ltd | 磁束検出装置 |
| FR2619454B1 (fr) * | 1987-08-14 | 1990-06-29 | Automatisme Tech Avancees | Procedes et portiques ou sas pour detecter le passage de masses metalliques |
| FI83266C (fi) * | 1988-09-12 | 1991-06-10 | Teknillinen Korkeakoulu | Foerfarande och anordning foer lokalisering av elektroder faestade vid kroppen av en maenniska, i synnerhet huvudet. |
| DE58905135D1 (de) * | 1988-09-22 | 1993-09-09 | Siemens Ag | Supraleitendes gradiometerschleifen-system einer mehrkanaligen messeinrichtung. |
| US4977896A (en) * | 1989-05-26 | 1990-12-18 | Biomagnetic Technologies, Inc. | Analysis of biological signals using data from arrays of sensors |
-
1990
- 1990-02-22 EP EP90103479A patent/EP0443069A1/de not_active Withdrawn
-
1991
- 1991-02-11 US US07/653,509 patent/US5136242A/en not_active Expired - Fee Related
- 1991-02-18 JP JP3046128A patent/JPH04215084A/ja active Pending
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006520892A (ja) * | 2003-03-14 | 2006-09-14 | エレクタ ニューロマグ オイ | 磁場のマルチチャネル測定値を処理するための方法とシステム |
| JP2008219455A (ja) * | 2007-03-05 | 2008-09-18 | Oki Electric Ind Co Ltd | 等化器及び等化方法 |
| JP4709789B2 (ja) * | 2007-03-05 | 2011-06-22 | Ｏｋｉセミコンダクタ株式会社 | 等化器 |
Also Published As
| Publication number | Publication date |
|---|---|
| US5136242A (en) | 1992-08-04 |
| EP0443069A1 (de) | 1991-08-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH04215084A (ja) | 電界又は磁界の界パターンを測定する方法 | |
| FI98337C (fi) | Menetelmä ja laite aivojen herätevasteiden ja spontaanitoiminnan sekä sydämestä mitattujen signaalien eri komponenttien erottelemiseksi toisistaan | |
| US7463024B2 (en) | Method and device for processing a multi-channel measurement of magnetic fields | |
| US4977896A (en) | Analysis of biological signals using data from arrays of sensors | |
| Taulu et al. | Applications of the signal space separation method | |
| US8838225B2 (en) | Analysis of multi-channel measurement data using orthogonal virtual channels | |
| Taulu et al. | Suppression of interference and artifacts by the signal space separation method | |
| JP5894156B2 (ja) | 生体磁場測定における干渉信号空間を調整する方法 | |
| Erne et al. | The Berlin magnetically shielded room (BMSR), Section B: Performances | |
| Taulu et al. | Novel noise reduction methods | |
| Yang et al. | SQUID array with optimal compensating configuration for magnetocardiography measurement in different environments | |
| JPH05337094A (ja) | 生体磁界の無雑音測定のためのシステムおよび方法 | |
| JP2019010483A (ja) | 磁界計測装置および計測磁界表示方法 | |
| Kong et al. | Multi-channel magnetocardiogardiography system based on low-Tc SQUIDs in an unshielded environment | |
| Taulu et al. | The signal space separation method | |
| Pasquarelli et al. | Calibration of a vector-MEG helmet system | |
| US6734353B2 (en) | Method for attenuating interference in a magnetically shielded room | |
| US20050182313A1 (en) | Method and apparatus for noise extraction in measurements of electromagnetic activity in biological sources | |
| CN109691997A (zh) | 脑磁场检测装置及脑磁场的检测方法 | |
| US20230218218A1 (en) | Magnetoencephalography apparatus and method | |
| JP7002416B2 (ja) | 磁界計測装置 | |
| Kandori et al. | Applying high-T/sub C/superconducting quantum interference devices with a room-temperature pickup coil in the measurement of impedance magnetocardiograms | |
| Ogata et al. | Study of spatial filter for magnetocardiography measurements without a magnetically shielded room | |
| Brauer et al. | Multipole analysis in magnetopneumography-a model study | |
| Bechstein et al. | Moderately shielded high-Tc SQUID system for rat MCG |