JP2002238869A - 生体磁場計測装置 - Google Patents

Abstract

(57)【要約】 【課題】 電気生理学的興奮の伝播過程を定量化可能な
生体磁場計測装置を提供する。 【解決手段】 生体磁場計測装置は，生体から発生する
生体磁場を計測する複数のＳＱＵＩＤ磁束計，周期的に
発生する生体信号を計測して収集する生体信号計測装置
６，７，複数方向から対として同時に計測された生体磁
場の信号及び生体信号の演算処理を行なう演算処理装置
１０，演算処理の結果を表示する表示装置１１とから構
成される。 【効果】 磁場源推定をせず，多数の等磁場線図を表示
せず，少数マップを用いて，興奮の時間変化を詳細に把
握できる。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は，生体の脳の神経活
動，心臓の心筋活動等により発生する生体磁場を，高感
度な量子干渉素子（ＳＱＵＩＤ）からなる複数の磁束計
を用いて計測する生体磁場計測方法及び生体磁場計測装
置に関する。

【０００２】

【従来の技術】従来，生体磁場計測の結果は，計測され
た磁場成分の時間変化波形，任意の時刻での磁場強度の
等しい点を結ぶ等磁場線図により表示されていた。例え
ば，直交座標でのＺ成分（Ｂ_ｚ），又は，極座標での同
径成分（Ｂ_ｒ）を計測し，Ｂ_ｚ又はＢ_ｒの値を等磁場線
図として表示することことが知られている（Ｈ．Ｈｏｓ
ａｋａ ａｎｄ Ｄ．Ｃｏｈｅｎ，Ｊ．Ｅｌｅｃｔｒｏ
ｃａｒｄｉｏｌ．，９−４，４２６（１９７６））。ま
た，直交座標ので接線成分（Ｂ_ｘ，Ｂ_ｙ）を計測し，各
成分毎に等磁場線図として表示する，あるいは，√
｛（Ｂ_ｘ）^２，（Ｂ_ｙ）^２｝により２次元磁場ベクトル
を得て等磁場線図として表示することが知られている
（Ｋ．Ｔｕｋａｄａ ｅｔ ａｌ．，Ｒｅｖｉｅｗ ｏ
ｆ ｔｈｅ Ｓｃｉｅｎｔｉｆｉｃ Ｉｎｓｔｒｕｍｅ
ｎｔｓ，６６，１０（１９９５））。更に，法線成分Ｂ
_ｚを計測したのち，Ｂ_ｚから接線成分（Ｂ_ｘ，Ｂ_ｙ）と
等価的な磁場成分を解析的に求める方法が知られている
（Ｔ．Ｍｉｙａｓｈｉｔａ ｅｔａｌ．，Ｐｒｏｃｅｅ
ｄｉｎｇｓ ２０ｔｈ Ｉｎｔｅｒｎａｔｉｏｎａｌ
Ｃｏｎｆｅｒｅｎｃｅ ＩＥＥＥ／ＥＭＢＳ（Ｈｏｎｇ
Ｋｏｎｇ），５２０−５２３（１９９８））。

【０００３】従来，生体磁場成分の解析結果は，磁場の
時間波形と等磁場線図を用いて表示されていた。また，
任意時刻での生体内の電流源の位置，大きさ，方向等
は，逆問題を解き推定し，これらを用いて心臓に於ける
不整脈の初期興奮位置や，脳に於ける癲癇の焦点位置等
の推定に使用されていた。心臓の心筋の興奮伝播過程
や，脳での神経興奮伝播等のある時間帯での動的な現象
を追跡するために，各時間毎の等磁場線図を多数く並べ
て表示したり，各時間毎に推定した電流源のベクトルの
軌跡を一つの図に表示していた（Ｎ．Ｉｚｕｍｉｄａ
ｅｔ ａｌ．，Ｊａｐａｎｅｓｅ Ｈｅａｒｔ Ｊｏｕ
ｒｎａｌ，７３１−７４２（１９９８））。

【０００４】

【発明が解決しようとする課題】多くの等磁場線図を並
べて，パターン認識により心臓や脳の動的な興奮伝播を
解析するのではなく，パターン認識によらず動的な興奮
伝播を行なうための定量化されグラフ表示，図形表示が
要求されていた。各瞬間毎に電流源を推定する方法で
は，電流源が局在している場合はダイポールモデルとし
て電流源は推定できるが，一般には，電流源は広がりを
持って分布している時間帯の方が多い。また，各瞬間毎
に於ける逆問題を解く場合，解が収束するまでに多くの
演算が必要とされる。特に，推定した電流源が作る計算
上の磁場分布と実際に計測された磁場分布との一致度が
悪い場合には，電流源の推定値が悪くなってしまう。こ
のため，ある時間帯での個々の瞬間毎で電流源を推定し
た場合，推定誤差が大きくなり，電流源の位置，大き
さ，方向の時間変化の連続性が途切れた解析結果を与え
るという問題があった。

【０００５】本発明の目的は，ダイポール（磁場源）推
定，多数の等磁場線図を表示することなく電気生理学的
な興奮の伝播過程を定量化できる生体磁場計測方法及び
生体磁場計測装置を提供することにある。

【０００６】

【課題を解決するための手段】本発明では，生体磁場計
測に於ける座標系として，直交座標系（ｘ，ｙ，ｚ）
（磁場成分をＢ_ｘ，Ｂ_ｙ，Ｂ_ｚとする），極直交座標系
（ｒ，θ，φ）が用いられる。計測対象が心臓等である
場合には，胸壁をｘｙ平面とする直交座標系（ｘ，ｙ，
ｚ）が用いられる。計測対象が脳部等である場合には，
頭部が球に近い形状であるため極直交座標系（ｒ，θ，
φ）（磁場成分をＢ_ｒ，Ｂ_θ，Ｂ_φとする）が用いられ
る。頭部の体表面に垂直な磁場成分（法線成分）は
Ｂ_ｚ，Ｂ_ｒで表わされ，生体の表面に接する面に平行な
成分（接線成分）は，Ｂ_ｘ，Ｂ_ｙ，Ｂ_θ，Ｂ_φで表わさ
れる。

【０００７】以下の説明では，直交座標系（ｘ，ｙ，
ｚ）を例にとって説明するが，極直交座標系（ｒ，θ，
φ）を用いる場合には，Ｂ_ｚをＢ_ｒに，Ｂ_ｘをＢに_θ，
Ｂ_ｙをＢ_φにそれぞれ読み替えるものとする。

【０００８】本発明の生体磁場計測装置では，１組のセ
ンサアレイを使用して，異なる多方向から生体磁場を計
測する。この時，多方向から生体磁場の計測結果を解析
するために，（１）各方向での生体磁場の計測と同時
に，生体信号計測装置として，心電計，心音計，脈波
計，脳波計等の何れかを使用し，心電波形，心音波形，
脈波形，脳波形等の何れかの，生体磁場信号以外の周期
的に発生する生体信号を，生体磁場の信号と対として計
測して収集する，あるいは，（２）電気刺激装置により
生体の一部を電気刺激して神経系を刺激する，音刺激装
置により音を発生し聴覚神経系を刺激する，匂い刺激装
置により匂いを発生し嗅覚神経系を刺激する，視覚刺激
装置により光信号又は色信号を発生し視覚神経系を刺激
する，触覚刺激装置により皮膚を刺激し触覚神経系を刺
激する等の，何れかの刺激信号の印加の開始と同期する
同期信号を，各方向での生体磁場の信号と対として収集
する。

【０００９】胸部から発する生体磁場（以下，心磁場と
いう）を，例えば，胸面，背面の２方向から，又は，胸
面，背面，胸部の右側面，及び左側面の４方向から計測
する。勿論，これらの方向以外から胸部から発する生体
磁場を計測しても良い。

【００１０】上記した何れかの刺激により頭部（脳）か
ら発する生体磁場（以下，脳磁場という）を，例えば，
頭部の前方，後方の２方向から，又は，頭部の前方，後
方の頭部の右側方，及び左側方の４方向から，あるい
は，頭部の前方，後方の頭部の右側方，及び左側方，頭
頂方向からの５方向から計測する。勿論，これの方向以
外から頭部から発する生体磁場を計測しても良い。

【００１１】ｔを時間変数，直交座標系（ｘ，ｙ，ｚ）
に於いて，ｘ，ｙをセンサアレイの個々のセンサが配置
される座標位置，生体の表面に接する面に平行な面をｘ
ｙ面，生体の表面に接する面に垂直な軸をｚとする。

【００１２】異なる多方向から計測された生体磁場波形
に対して，各方向毎に以下の処理を行なう。周期的に発
生する生体信号を，生体磁場の信号と対として計測して
収集する場合には，時間変数をｔとして，ｍ＝１，２，
…，Ｍの複数方向から計測された生体信号の波形Ｗ
_ｍ（ｔ）の時間軸が，共通の原点（ｔ＝０）を持つよう
に，計測された各生体信号の波形Ｗ_ｍ（ｔ）（ｍ＝１，
２，…，Ｍ）の時間軸に対して変換Ｔ_ｍ（ｍ＝１，２，
…，Ｍ）を行なう。生体信号Ｗ_ｍ（ｔ）と対をなす生体
磁場の信号の波形Ｆ_ｍ（ｍ＝１，２，…，Ｍ）の時間軸
に対して変換Ｔ_ｍ（ｍ＝１，２，…，Ｍ）変換を行な
う。刺激信号の印加の開始に同期する同期信号を，生体
磁場の信号と対として収集する場合には，ｍ＝１，２，
…，Ｍの複数方向から計測された生体磁場の信号の波形
Ｆ_ｍ（ｍ＝１，２，…，Ｍ）の時間軸が，上記の同期信
号が収集された時点で共通の原点（ｔ＝０）を持つよう
に，時間軸の変換Ｔ_ｍ’（ｍ＝１，２，…，Ｍ）を行な
う。上記の変換Ｔ_ｍ，Ｔ_ｍ’（ｍ＝１，２，…，Ｍ）
は，時間軸の平行移動の変換である。

【００１３】共通の原点（ｔ＝０）を持つ，複数の各方
向から計測された生体磁場（心磁場又は脳磁場）の信号
の波形に対してそれぞれ，以下の演算処理を行なう。

【００１４】生体磁場として，生体の表面に接する面に
垂直な磁場成分Ｂ_ｚ（ｘ，ｙ，ｔ）を計測する場合に
は，垂直な磁場成分Ｂ_ｚ（ｘ，ｙ，ｔ）のｘ方向の変化
量，∂Ｂ_ｚ（ｘ，ｙ，ｔ）／∂ｘと，Ｂ_ｚ（ｘ，ｙ，
ｔ）のｙ方向の変化量，∂Ｂ_ｚ（ｘ，ｙ，ｔ）／∂ｙと
を求め，（数１），（数２）により，２乗和の平方根，
即ち，２次元磁場ベクトルＩ（ｘ，ｙ，ｔ）の大きさ
（以下，単に，ベクトル強度という）とその位相θ
（ｘ，ｙ，ｔ）を求める演算処理をそれぞれ行なう。

【００１５】

【数１】 Ｉ（ｘ，ｙ，ｔ）＝√｛（∂Ｂ_ｚ（ｘ，ｙ，ｔ）／∂ｘ）^２＋（∂Ｂ_ｚ（ｘ，ｙ ，ｔ）／∂ｙ）^２｝ …（数１）

【００１６】

【数２】 θ（ｘ，ｙ，ｔ）＝−ｔａｎ^−１｛（−∂Ｂ_ｚ（ｘ，ｙ，ｔ）／∂ｘ）／（∂Ｂ _ｚ （ｘ，ｙ，ｔ）／∂ｙ）｝ …（数２） 生体から発する磁場の接線成分（生体の表面に接する面
に平行な成分）Ｂ_ｘ，Ｂ_ｙを計測する場合には，（数
３），（数４）により，接線成分Ｂ_ｘとＢ_ｙの２乗和の
平方根からベクトル強度Ｉ（ｘ，ｙ，ｔ）とその位相θ
（ｘ，ｙ，ｔ）を求める演算処理をそれぞれ行なう。

【００１７】

【数３】 Ｉ（ｘ，ｙ，ｔ）＝√｛（Ｂ_ｘ（ｘ，ｙ，ｔ））^２＋（Ｂ_ｙ（ｘ，ｙ，ｔ））^２ ｝ …（数３）

【００１８】

【数４】 θ（ｘ，ｙ，ｔ）＝−ｔａｎ^−１｛−Ｂ_ｘ（ｘ，ｙ，ｔ）/Ｂ_ｙ（ｘ，ｙ，ｔ） ｝ …（数４） 次に，計測された生体磁場（心磁場又は脳磁場）の各時
間ｔに於ける最大ベクトル強度Ｉ_ｍａｘ（ｘ_ｉ，ｙ_ｊ，
ｔ）とその位相θ（ｘ_ｉ，ｙ_ｊ，ｔ）を求める。時間ｔ
に於いて，センサのｉ番目のｘ座標位置，ｊ番目のｙ座
標位置に於いて，即ち，（ｉ，ｊ）チャネルに於いて，
ベクトル強度Ｉ（ｘ，ｙ，ｔ）が最大を与える。求めら
れた各時間ｔに於ける最大ベクトル強度Ｉ
_ｍａｘ（ｘ_ｉ，ｙ_ｊ，ｔ）とその位相θ（ｘ_ｉ，ｙ_ｊ，
ｔ）を時間変数ｔに対して表示する。この表示を，時間
−強度プロット（ｔ−Ｉ_ｍａｘ），時間−位相プロット
（t−θ）という。

【００１９】以上の結果，複数の各方向から計測され，
共通の共通の時間原点（ｔ＝０）を持つ生体磁場（心磁
場又は脳磁場）の信号の波形から，時間−強度プロット
（ｔ−Ｉ_ｍａｘ），及び時間−位相プロット（t−θ）
が得られる。この結果，生体磁場の各計測面毎に時間−
強度プロット（ｔ−Ｉ_ｍａｘ），及び時間−位相プロッ
ト（t−θ）を比較表示できる。

【００２０】また，複数の各方向から計測され，共通の
共通の時間原点（ｔ＝０）を持つ生体磁場（心磁場又は
脳磁場）の信号の波形から得られる，全てのセンサの位
置（ｘ_ｉ，ｙ_ｊ），即ち，全チャネルに於ける，ベクト
ル強度Ｉ（ｘ，ｙ，ｔ）とその位相θ（ｘ，ｙ，ｔ）を
同一の表示画面に表示することもできる。この表示を時
間−位相・強度プロット（ｔ−θ・Ｉ）という。この表
示では，位相θ（ｘ，ｙ，ｔ）を時間変数ｔに対してプ
ロットし，ベクトル強度Ｉ（ｘ，ｙ，ｔ）は，プロット
する色，プロットする色の濃淡，又は，プロットする記
号の大きさを，ベクトル強度Ｉ（ｘ，ｙ，ｔ）に対応し
て変化させて表示する。

【００２１】以上のようにして，ダイポール推定をする
ことなく，また，多数の等磁場線図を表示することな
く，心磁場又は脳磁場を複数方向から計測することによ
り，電気生理学的な興奮の伝播過程を定量化して表示で
きる。

【００２２】

【発明の実施の形態】本発明の実施例の生体磁場計測装
置では，複数のＳＱＵＩＤ磁束計により生体から発生す
る生体磁場を計測する。この時，生体磁場以外の周期的
に発生する生体信号を計測して収集する生体信号計測装
置，あるいは，各種の神経系の何れかを刺激する刺激信
号とこの刺激信号の印加の開始に同期する同期信号を発
生する刺激装置を，使用する。

【００２３】複数方向から対として同時に計測され収集
された生体磁場の信号及び生体信号，あるいは，刺激信
号の印加の開始に同期する同期信号と対として，複数方
向から計測され収集された生体磁場の信号は，演算処理
装置により演算処理され，演算処理結果が表示装置に表
示される。

【００２４】演算処理装置は，複数方向から計測され収
集された生体磁場の信号が，共通原点（ｔ＝０）をもつ
ように，生体磁場の信号の時間軸の平行移動を行なう時
間軸変換を行ない，共通原点（ｔ＝０）をもつ生体磁場
の信号に対して演算処理を行なう。

【００２５】生体信号計測装置を使用する場合には，時
間軸変換は，時間変数をｔとして，複数方向から計測さ
れた生体信号の波形の時間軸が共通原点（ｔ＝０）を持
つように，各生体信号の波形の時間軸の変換を行なう，
この時，生体信号と対をなす生体磁場の信号の波形の時
間軸に対しても，同一の変換を行なう。刺激装置を使用
する場合には，複数方向から計測された生体磁場の信号
の時間軸が，同期信号が収集された時点で共通の原点
（ｔ＝０）を持つように，時間軸変換が実行される。

【００２６】また，本発明の他の実施例では，演算処理
装置は，共通原点（ｔ＝０）をもつ生体磁場（心磁場又
は脳磁場）の信号を用いて，生体の表面に接する面に垂
直な方向をｚ方向，ｚ方向に直交し生体の表面に接する
面に水平な方向をｘ方向及びｙ方向とし，生体磁場の各
計測点（ｘ，ｙ）で２次元磁場ベクトルの大きさ及び／
又はｘｙ面での方向を表わす角度を求める。表示装置
に，２次元磁場ベクトルの大きさ及び／又は角度の，共
通原点（ｔ＝０）を原点とする時間変化が表示される。

【００２７】本発明の他の実施例の生体磁場計測装置
は，生体の胸部又は頭部から発生する生体磁場（心磁場
又は脳磁場）の信号を計測する複数のＳＱＵＩＤ磁束計
と，生体磁場の信号の演算処理を行なう演算処理装置
と，演算処理の結果を表示する表示装置とを具備する。

【００２８】演算処理装置は，生体の表面に接する面に
垂直な方向をｚ方向，ｚ方向に直交し生体の表面に接す
る面に水平な方向をｘ方向及びｙ方向として，生体磁場
の信号から生体磁場の各計測点（ｘ，ｙ）に於ける２次
元磁場ベクトルを求め，複数の計測点（ｘ，ｙ）に於け
る２次元磁場ベクトルのうち最大の２次元磁場ベクトル
の大きさ及び／又はｘｙ面での方向を表わす角度を，生
体磁場が計測された複数の時点で求める。そして，表示
装置に，最大の２次元磁場ベクトルの大きさ及び／又は
角度の時間変化が表示される。

【００２９】別の実施例では，各計測点（ｘ，ｙ）に於
ける上記の２次元磁場ベクトルの大きさ及び角度を，生
体磁場が計測された複数の時点で求める。そして，表示
装置に，各計測点（ｘ，ｙ）に於ける２次元磁場ベクト
ルの大きさ及びｘｙ面での方向を表わす角度の時間変化
が，２次元磁場ベクトルの大きさをプロット点の大きさ
に比例させて，又は，カラーで区別して表示される。

【００３０】また，本発明の別の実施例の生体磁場計測
装置は，生体の頭部から発生する生体磁場の信号を計測
する複数のＳＱＵＩＤ磁束計と，生体を刺激する信号と
生体を刺激する信号の発生の開始に同期する同期信号を
発生する刺激装置と，生体の頭部の複数方向から計測さ
れた生体磁場（脳磁場）の信号の演算処理を行なう演算
処理装置と，演算処理の結果を表示する表示装置とを具
備する。

【００３１】演算処理装置は，生体の表面に接する面に
垂直な方向をｚ方向，ｚ方向に直交し生体の表面に接す
る面に水平な方向をｘ方向及びｙ方向とし，時間変数を
ｔとして，複数方向から計測された生体信号の波形の時
間軸が共通原点（ｔ＝０）を持つように，同期信号に基
づいて複数方向から計測された生体信号の波形の時間軸
の変換を行なう。

【００３２】更に，生体磁場を計測した複数方向につい
て，共通原点（ｔ＝０）をもつ生体磁場の信号から生体
磁場の各計測点（ｘ，ｙ）に於ける２次元磁場ベクトル
を求め，複数の計測点（ｘ，ｙ）に於ける２次元磁場ベ
クトルのうち最大の２次元磁場ベクトルの大きさ及び／
又はｘｙ面での方向を表わす角度を，共通原点（ｔ＝
０）を原点とする時間軸の各点で求める。そして，表示
装置に，最大の２次元磁場ベクトルの大きさ及び／又は
角度の時間変化が，生体磁場を計測した複数方向につい
て表示される。 （第１の実施例）図１は，本発明の第１の実施例の生体
磁場計測装置の概略構成を説明する斜視図である。心臓
から発する磁場（心磁場）の計測（以下，心磁場計測と
いう）を行なう生体磁場計測装置は，量子干渉素子（Ｓ
ＱＵＩＤ）からなる複数の磁場センサを用いる。環境磁
場雑音の影響を除去するために，心磁場計測は磁場シー
ルドルーム１の内部で行なわれる。被検体２はベッド４
に横たわる。

【００３３】ここで，ｘｙ面がベッドの面となるように
直交座標系（ｘ，ｙ，ｚ）を設定する。被検体２の上方
に，ＳＱＵＩＤとそのＳＱＵＩＤに接続した検出コイル
とが一体化された磁場センサを複数個収納し，液体Ｈｅ
を満たしたデュワ３が配置される。デュワ３はガントリ
５により床に固定されている。磁場センサからの出力
は，検出コイルが検出した磁場強度に比例する電圧を出
力するＦＬＬ（Flux Locked Loop）回路８に入力され
る。

【００３４】ＦＦＬ回路８は，ＳＱＵＩＤの出力を一定
に保つようＳＱＵＩＤに入力された生体磁場の変化を帰
還コイルを介してキャンセルしている。帰還コイルに流
した電流を電圧に変換することにより，生体磁場の信号
の変化に比例した電圧出力が得られる。この電圧出力
は，アンプフィルタ回路９の増幅器により増幅され，フ
イルター回路により周波数帯域が選択され，データ収録
解析装置（演算処理装置）１０でＡＤ変換されデータ収
録される。

【００３５】データ収録解析装置１０では，各種の演算
処理が実行され，演算処理結果がデイスプレイ１１に表
示され，更に，プリンタにより出力される。心磁場の計
測と同時に心電図の計測も行なう。被検体２の手首，足
首に心電計用電極６を貼りつけ，四肢誘導による電位が
心電計７に導かれる。心電計７の出力は，生体磁場計測
装置のアンプフィルタ回路９に入力され，増幅，周波数
帯域が選択され，心磁場と同様に演算処理され，心電波
形がデイスプレイ１１に表示される。心磁場の計測を正
面，及び背面から計測する時は，被検体は，それぞれ仰
向け，うつ伏せの状態となり，心磁場が計測される。

【００３６】図２は，本発明の第１の実施例に於ける磁
場センサの配置を説明する斜視図である。生体磁場の法
線成分を検出するコイルはｚ方向に垂直な面を持つ。磁
場センサ２０−ｉ（ｉ＝１１，１２，…，１８；ｉ＝２
１，２２，…，２８；ｉ＝３１，３２，…，３８；ｉ＝
４１，４２，…，４８；ｉ＝５１，５２，…，５８；ｉ
＝６１，６２，…，６８；ｉ＝７１，７２，…，７８；
ｉ＝８１，８２，…，８８）はデュワ内部の底部から垂
直に立位で設置される。各センサ間の距離は，生体磁場
の法線成分のｘ方向，ｙ方向での変化量を正確に求める
ように，ｘ方向，ｙ方向で等間隔（２５ｍｍ）に配置し
た。センサを，ｘ方向，ｙ方向にそれぞれ８個，正方格
子状に配置し，６４チャンネルとした。

【００３７】図３は，本発明の第１の実施例に於ける生
体磁場の法線成分を検出する磁場センサの構成を説明す
る斜視図である。磁場センサ２０は，体表面に対して垂
直な成分Ｂ_ｚを計測するセンサであり，超伝導線（Ｎｂ
ーＴｉ線）で作製されるコイルはｚ方向に垂直な面を持
つ。このコイルは２つの逆向きのコイルを組み合わせた
ものであり，生体に近い位置に配置される検出コイル２
２と，生体から遠い位置に配置され，外部磁場雑音を除
去する参照コイル２３とにより一次微分コイルを形成し
ている。検出コイル２２，及び参照コイル２３の各コイ
ル径を２０ｍｍφ，検出コイル２２，及び参照コイル２
３の各コイル間の距離（ベースライン）を５０ｍｍとし
た。

【００３８】外部磁場雑音は生体より遠い信号源から生
じており，これらは検出コイル２２及び参照コイル２３
により同じように検出される。一方，生体内の磁場源は
コイルに近いため，生体磁場は検出コイル２２によりよ
り強く検出される。検出コイル２２により生体磁場の信
号と外部磁場雑音が検出され，参照コイル２３では外部
磁場雑音のみが検出される。従って，両者のコイルで検
出された磁場の差からよりＳ/Ｎの高い計測ができる。
一次微分コイルはＳＱＵＩＤ２１を実装した実装基板の
超伝導配線を介してＳＱＵＩＤのインプットコイルに接
続され，コイルで検出された生体磁場の信号はＳＱＵＩ
Ｄに伝達される。

【００３９】磁場センサを内蔵したデュワは，ベッドに
横たわった被検体の上方に配置され，心臓から発生する
心磁場が計測される。ここで，体軸の方向をｙ軸とし，
ｙ軸に直交するｘ軸とする。

【００４０】図４は，本発明の第１の実施例に於ける磁
場センサの配置と被検体２の胸部の正面，背面との位置
関係を説明する図である。図４に於いて，白丸は，磁場
センサが８×８のアレイ状に配列される位置，即ち，磁
場の計測点３０示す。被検体２の胸部の正面での計測基
準点３１と背面の計測基準点３１’はｘｙ座標上で同じ
点とする。図４に示す例では，正面での計測基準点３１
は計測面の下から２行，左から３列に位置する磁場セン
サに対応し，背面では計測基準点は２行，６列に位置す
る磁場センサに対応する。但し，背面での座標系はｘ座
標を正面での座標系とは逆にとる。例えば，正面での１
行，１列に位置する磁場センサは，背面での１行，８列
に位置する磁場センサに対応し，正面での８行，８列に
位置する磁場センサは，背面での８行，１列に位置する
磁場センサに対応する。

【００４１】図５は，本発明の第１の実施例に於ける磁
場センサの配置と人体の胸部との位置を合わせを行なう
方法を説明する斜視図である。センサアレイの計測基準
点３１と被検体の基準点４０とを合わせるための各種の
機構，方法が知られている。図５に示す例では，直交座
標系のｘｚ面に平行な面内に扇状に広がるｘ軸ライン形
成用ビーム４３を発生するｘ方向レーザ光源４１，直交
座標系のｙｚ面に平行な面内に扇状に広がるｙ軸ライン
形成用ビーム４４を発生するｙ方向レーザ光源４２が，
位置合わせのために使用される。デュワ３の外周面に
は，直交座標系のｘｚ面の位置を示すｘｚ標識（マー
ク）４５，直交座標系のｙｚ面の位置を示すｙｚ標識
（マーク）４６が印されている。

【００４２】ｘ軸ライン形成用ビーム４３が，被検体２
及びデュワ３の外周面に照射され，被検体の体表に設定
された基準点４０とデュワ３のｘｚ標識（マーク）４５
とを照射するように，ｘ方向レーザ光源４１の位置を調
整する。同様にして，ｙ軸ライン形成用ビーム４４が，
被検体２及びデュワ３の外周面に照射され，被検体の体
表に設定された基準点４０とデュワ３のｙｚ標識（マー
ク）４６とを照射するように，ｙ方向レーザ光源４２の
位置を調整する。これによりセンサと生体の位置とを調
整できる。ビーム４３，４４は交差してｚ軸に平行な交
差線４９を形成する。

【００４３】図６は，本発明の第１の実施例に於ける生
体磁場の計測，計測された信号の解析の流れを説明する
図である。図６に示す解析では，計測された心磁場の各
時間に於ける最大ベクトルを選択する。まず始めに，正
面，及び背面の２方向から心磁場計測を行なう。心磁場
の法線成分Ｂ_ｚを計測する。心磁場計測と同時に心電図
も計測を行なう。次に，正面，及び背面からの心磁場計
測と同時に計測された心電図の波形が，同じ時間で同一
波形になるように，時間軸を調整する。即ち，各心電図
の時間軸を平行移動して，計測された各心電図での同じ
時相が同じ時間になるようにする。この各心電図で調整
した時間軸の平行移動を対応する心磁場の計測データに
ついても適用して，心電図データに基づく心磁場の計測
データ（以下，心磁場波形という）の時間軸の設定を行
なう。

【００４４】次に，図４で説明したように正面，背面で
の心磁場の計測データの対応をつけるため，背面でのｘ
座標を反転してセンサアレイの位置と計測データとの対
応を変える。次に，計測された生体面に垂直な磁場成分
Ｂ_ｚ（ｘ，ｙ，ｔ）のｘ方向での変化量∂Ｂ_ｚ（ｘ，
ｙ，ｔ）／∂ｘと，Ｂ_ｚ（ｘ，ｙ，ｔ）の方向での変化
量∂Ｂ_ｚ（ｘ，ｙ，ｔ）／∂ｙを計算により求め，正
面，及び背面でそれぞれ計測された観測点６４点に於け
るベクトル強度Ｉ（ｘ，ｙ，ｔ）を（数１）に基づい
て，ベクトル位相θ（ｘ，ｙ，ｔ）を（数２）に基づい
てそれぞれ求める。

【００４５】次に，心磁場の各計測面（正面，背面）毎
の各時間ｔ毎にベクトル強度Ｉ（ｘ，ｙ，ｔ）の最大ベ
クトルを抽出する。即ち，正面，及び背面のそれぞれで
の心磁場の計測データから求められたベクトル強度Ｉ
（ｘ，ｙ，ｔ）の各時間ｔに於ける，最大ベクトル強度
Ｉ_ｍａｘ（ｘ_ｉ，ｙ_ｊ，ｔ）とその位相θ（ｘ_ｉ，
ｙ_ｊ，ｔ）を求める。次に，求められた各時間ｔに於け
る最大ベクトル強度Ｉ_ｍａｘ（ｘ_ｉ，ｙ_ｊ，ｔ）とその
位相θ（ｘ_ｉ，ｙ_ｊ，ｔ）を時間変数ｔに対して表示す
る。即ち，心磁場の各計測面（正面，背面）毎に時間−
強度プロット（ｔ−Ｉ _ｍａｘ），及び時間−位相プロッ
ト（t−θ）を表示する。以下，図６の流れ図に基づい
て得られた具体例について説明する。

【００４６】図７は，健常者を被検体として，正面，及
び背面から計測された心磁場波形と，心磁場波形の計測
と同時に計測された心電図波形の例を示す図である。図
７では，心磁場波形を，６４チャンネルの全ての時間波
形を重ね合わせ表示しており，心電図波形は第II誘導の
波形を示している。正面，及び背面から計測された心磁
場波形の時間軸を調整するため，正面からの心磁場計測
時に計測された心電波形のＱＲＳ波と呼ばれる心室が脱
分極した時間帯の開始点の時間Ｔ１を０に設定する。次
に，背面からの心磁場計測時に計測された心電波形のＱ
ＲＳ波の開始点の時間Ｔ２を０に設定し，２つの心電波
形の時間軸を合わせる。正面からの心磁場波形の時間原
点はｔ＝Ｔ１＝０に，背面からの心磁場波形の時間原点
はｔ＝Ｔ２＝０に，それぞれ設定する。

【００４７】図８は，本発明の第１の実施例に於いて，
正面から計測された健常者の心磁場波形の心室が脱分極
する過程であるＱＲＳ波の開始点から３０ｍｓ経過後の
時間での電流アロー図，及び最大電流ベクトルを示す図
である。６４点の計測点での各矢印の大きさは，（数
１）に基づくベクトル強度を示し，位相は（数２）に基
づいて計算する。６４個のベクトルの中のもっとも大き
いベクトルを最大ベクトルとして選択する。

【００４８】図９は，本発明の第１の実施例に於いて，
最大ベクトルの位相θの基準を示す図である。図９に於
いて，ｘ軸の正方向（右方向）をθ＝０度として，時計
周りを位相のプラス方向，反時計周りを位相のマイナス
方向とする。

【００４９】図１０は，本発明の第１の実施例に於い
て，正面，及び背面から計測された健常者の心磁場波形
の心室が脱分極する過程であるＱＲＳ波の開始点から２
００ｍｓまでの時間帯での時間−強度プロット（ｔ−Ｉ
_ｍａｘ）と，時間−位相プロット（ｔ−θ）の表示例を
示す図である。正面，背面から計測された心磁波形から
得られた最大ベクトル強度（ｐＴ／ｃｍ），及び，位相
の時間変化は異なるパターンを示している。しかし，最
大ベクトル強度（ｐＴ／ｃｍ）は，ＱＲＳ波の開始点か
ら約２０ｍｓから約１００ｍｓの時間帯で大きな値を持
つパターンを示す点では共通する。

【００５０】次に，心室に刺激伝導障害がある脚ブロッ
クのうち右脚ブロックの患者に関する計測例を説明す
る。

【００５１】図１１は，本発明の第１の実施例に於い
て，正面，及び背面から計測された右脚ブロックの患者
の心磁波形の心室が脱分極する過程であるＱＲＳ波の開
始点から２００ｍｓまでの時間帯での時間−強度プロッ
ト（ｔ−Ｉ_ｍａｘ）と，時間−位相プロット（ｔ−θ）
の表示例を示す図である。図１０に示す健常者の場合の
結果と異なり，正面から計測された心磁場波形から得ら
れたベクトル強度が，ＱＲＳ波の開始点から約６０ｍｓ
経過後から約１６０ｍｓの長い時間帯わたり大きな値を
示しており，心臓の活動時間も長いことが分かる。ま
た，正面から計測された心磁場波形から得られた位相の
時間変化を，図１０に示す健常者の場合と比較すると，
ＱＲＳ波の開始点から初期の時間帯での変化が，健常者
の場合に比べて小さいことが分かる。このように，心臓
の興奮伝播の異常が，本発明により容易に判別できるよ
うになった。 （第２の実施例）本発明の第２の実施例では，心磁場の
各計測面（正面，背面）毎に時間−強度プロット（ｔ−
Ｉ_ｍａｘ），及び時間−位相プロット（t−θ）の表示
に加えて，６４チャンネルの全ての心磁場波形のデータ
を表示する。

【００５２】本発明の第２の実施例に於いて６４チャン
ネルの全ての各チャンネルの各時間でのベクトルの時間
−位相・強度プロット（ｔ−θ・Ｉ）を行なう場合の，
生体磁場の計測，計測された信号の解析の流れを説明す
る図である。図１２に示す流れ図では，各計測点（チャ
ネル）でのベクトル強度及びその位相を求める処理まで
は，図６と同じである。図１２に示す流れ図では，最大
ベクトルを選択せずにそのまま，正面，及び背面から計
測された心磁場のデータからそれぞれ得られた，６４チ
ャンネルの全てのチャネルに於けるベクトルの時間−位
相・強度プロット（ｔ−θ・Ｉ）を表示する。

【００５３】図１３は，本発明の第２の実施例に於い
て，健常者の正面から計測された心磁場波形の心室が脱
分極する過程であるＱＲＳの開始点から２００ｍｓまで
の時間帯での６４チャンネルの全チャンネルのベクトル
の時間−位相強度プロットの表示例を示す図である。

【００５４】図１４は，本発明の第２の実施例に於い
て，健常者の背面から計測された心磁場波形の心室が脱
分極する過程であるＱＲＳの開始点から２００ｍｓまで
の時間帯での６４チャンネルの全チャンネルのベクトル
の時間−位相強度プロットの表示例を示す図である。

【００５５】図１５は，本発明の第２の実施例に於い
て，正面から計測された右脚ブロックの患者の心磁波形
の心室が脱分極する過程であるＱＲＳ波の開始点から２
００ｍｓまでの時間帯での６４チャンネルの全チャンネ
ルのベクトルの時間−位相強度プロットの表示例を示す
図である。図１５のプロットのパターンと図１３のプロ
ットのパターンとは一見して大きな差が認められる。

【００５６】図１６は，本発明の第２の実施例に於い
て，背面から計測された右脚ブロックの患者の心磁波形
の心室が脱分極する過程であるＱＲＳ波の開始点から２
００ｍｓまでの時間帯での６４チャンネルの全チャンネ
ルのベクトルの時間−位相強度プロットの表示例を示す
図である。図１６のプロットのパターンと図１４のプロ
ットのパターンとは一見して大きな差が認められる。

【００５７】図１３から図１６では，各時間に於ける６
４点の各計測点でのベクトル強度を表示するプロット点
の濃度で区別して，位相を縦軸として表示しているの
で，容易に表示されたパターンの相違を識別できる。ベ
クトル強度の表示として，プロット点の濃淡で表示する
他に，プロット点の色の違い（カラースケール）や，プ
ロット点の大きさで表示することも可能である。このよ
うなプロットにより，健常者との相違点を容易に検出で
きる。 （第３の実施例）第３の実施例では，生体磁場のｘ方
向，及びｙ方向の磁場成分を計測する。

【００５８】図１７は，本発明の第３の実施例に於いて
使用する生体磁場の接線成分Ｂ_ｘ，及びＢ_ｙ成分を検出
する磁場センサの構成例の概略を説明する図である。図
１７に示す磁場センサ２０’は，平面型のコイルを使用
している。

【００５９】ｘ方向の磁場を計測するセンサは検出コイ
ル２２’−１と参照コイル２３’−１とが一つの平面に
並んで配置され，ｘ方向の磁場を計測する一次微分コイ
ルを構成している。検出コイル２２’−１，参照コイル
２３’−１のコイルの大きさを２０ｍｍ×２０ｍｍの正
方形とし，検出コイル２２’−１，参照コイル２３’−
１の各コイルの中心間の距離（ベースライン）を５０ｍ
ｍとした。ｘ方向の磁場を計測する一次微分コイルはＳ
ＱＵＩＤ２１−１’を実装した実装基板の超伝導配線を
介してＳＱＵＩＤのインプットコイルに接続され，コイ
ルで検出された生体磁場の信号はＳＱＵＩＤに伝達され
る。

【００６０】ｙ方向の磁場を計測するセンサは検出コイ
ル２２’−２と参照コイル２３’−２とが一つの平面に
並んで配置され，ｙ方向の磁場を計測する一次微分コイ
ルを構成している。検出コイル２２’−２，参照コイル
２３’−２のコイルの大きさを２０ｍｍ×２０ｍｍの正
方形とし，検出コイル２２’−２，参照コイル２３’−
２の各コイルの中心間の距離（ベースライン）を５０ｍ
ｍとした。ｙ方向の磁場を計測する一次微分コイルはＳ
ＱＵＩＤ２１−２’を実装した実装基板の超伝導配線を
介してＳＱＵＩＤのインプットコイルに接続され，コイ
ルで検出された生体磁場の信号はＳＱＵＩＤに伝達され
る。

【００６１】４角柱の支持体の互いに直交する２面に，
上記のｘ成分検出用の磁場センサとｙ成分検出用の磁場
センサを張り付けることにより，生体磁場のｘ成分，及
びｙ成分を計測できる磁場センサを形成している。図１
７に示すような４角柱の形状を持つ磁場センサは，図２
に示すようにセンサアレイ状に配置される。計測された
接線成分Ｂ_ｘ，Ｂ_ｙからベクトル強度Ｉ（ｘ，ｙ，ｔ）
を（数３）に基づいて，その位相θ（ｘ，ｙ，ｔ）を
（数４）に基づいて求める。

【００６２】図６，及び図１２で説明したように生体磁
場の法線成分のｘ方向及びｙ方向での微分を求める代わ
りに，計測された２つの接線成分から（数３），（数
４）に基づいて求めたベクトル強度Ｉ（ｘ，ｙ，ｔ）と
その位相θ（ｘ，ｙ，ｔ）を用いて，先に説明した第１
実施例，第２の実施例と同様にして，最大ベクトルの時
間−強度プロット（ｔ−Ｉ_ｍａｘ），及び時間−位相プ
ロット（ｔ−θ），あるいは，６４チャンネルの全ての
チャンネルについての時間−位相・強度プロット（ｔ−
θ・Ｉ）を求め，表示できる。 （第４の実施例）第４の実施例では，平面型微分コイル
を用いて法線成分のｘ方向，及びｙ方向の一次微分成分
を直接計測して解析する。

【００６３】図１８は，本発明の第４の実施例に於いて
使用する，コイルに同一の方向から磁場が入射したと
き，互いに逆向きの電流が流れる直線部分を持つ半円形
状のコイルが並んで形成され，コイル全体として円の形
状を持つ微分コイル５０の形状の例を示す図である。図
１８に示す平面微分型磁場センサ２０”では，互いに逆
向きの電流が流れる直線部分と直交する方向で磁場が微
分された値が検出できる。微分コイル５０は，ＳＱＵＩ
Ｄ２１”を実装した実装基板の超伝導配線を介してＳＱ
ＵＩＤのインプットコイルに接続され，微分コイルで検
出された生体磁場の信号はＳＱＵＩＤに伝達される。互
いに逆向きの電流が流れる直線部分を，ｘ方向に持つ平
面型微分コイルとｙ方向に持つ平面型微分コイルを準備
する。

【００６４】図２に示すようなセンサアレイの各計測点
に，ｘ方向，ｙ方向の微分成分を検出する平面型微分コ
イルを重ねて配置することにより，６４点の各計測点で
生体磁場の法線成分のｘ方向，ｙ方向での微分値が直接
計測できる。計測されたｘ方向，ｙ方向での微分値か
ら，（数１）に基づいてベクトル強度Ｉ（ｘ，ｙ，ｔ）
を求め，（数２）に基づいてその位相θ（ｘ，ｙ，ｔ）
を求める。先に説明した第１の実施例〜第３の実施例と
同様にして，最大ベクトルの時間−強度プロット（ｔ−
Ｉ_ｍａｘ），及び時間−位相プロット（ｔ−θ），ある
いは，６４チャンネルの全てのチャンネルについての時
間−位相・強度プロット（ｔ−θ・Ｉ）を求め，表示で
きる。 （第５の実施例）図１９は，本発明の第５の実施例に於
いて使用する，正面，及び背面の方向から同時に心磁場
計測を行なう生体磁場計測装置の構成例の概略を説明す
る図である。第５の実施例の生体磁場計測装置の構成例
では，被検体２が，正面，背面の方向から心磁場を計測
するために姿勢を変える必要がなく，一度に同時に２方
向から心磁場を計測できる。上部デュワ３’−１と下部
デュワ３’−２の２つデュワを備え，それぞれに，第１
の実施例〜第４の実施例で説明したような何れかのセン
サアレイが配置される。上部デュワ３’−１と下部デュ
ワ３’−２はガントリ５’に保持されている。ベッド
４’の脚には，水平方向に移動が可能な滑車１３が付け
られている。被検体２がベッド４’に乗った後，レール
１２に沿って，上部デュワ３’−１と下部デュワ３’−
２との間の定位置に，ベッド４’を移動させる。心電図
等の周期的な生体信号を計測することなく，被検体２の
上下に配置されるセンサアレイより，心磁場が２方向か
ら同時に計測できるので，２方向から計測された心磁場
のデータの時間軸を調整する必要がない。 （第６の実施例）第６の実施例では，音の刺激により聴
性誘発される磁場（脳磁場）を計測する。 図２０は，
本発明の第６の実施例に於いて聴性誘発された脳磁場を
計測する生体磁場計測装置の構成例を説明する斜視図で
ある。被検体はベッド４に横たわり，デュワ３の底面に
計測しようとする頭部の面を近づけて，脳磁場は計測さ
れる。

【００６５】図２０に示す構成では，１ｋＨｚで５０ｍ
ｓ保持時間幅を持つトーンバースト音を音刺激装置１８
３によって生成する。音刺激の間隔は０．３Ｈｚ（約
３．３秒に１回の頻度）で行なう。音刺激の印加のタイ
ミングに合わせて同期信号１８４を発生して，データ収
録解析装置（演算処理装置）１０へ同期信号１８４が入
力される。入力された同期信号１８４を利用して，脳磁
場の波形の時間軸を一致させた後，加算平均化処理を行
ない信号雑音比を向上させる。

【００６６】音刺激装置１８３によって生成されたトー
ンバースト音はエアーチューブ１８２とアダプタ１８１
を通して左耳に入力される。図２０に図示しないが，右
耳にホワイト雑音の音を常時与えて外来からの音による
影響がないように計測を行なう。脳磁場はデュワ３内の
磁場センサによって計測される。磁場センサはＦＬＬ回
路８にって駆動され，ＦＬＬ回路８の出力はアンプフィ
ルタ回路９を通ってデータ収録解析装置１０にデジタル
データとして収集され記録される。データ収録解析装置
１０，ＦＬＬ回路８，アンプフィルタ回路９等を制御す
るための画面やデータ解析結果を表示する画面等は，デ
ィスプレイ１１に表示される。図２０に示すベッド４，
デュワ３以外の構成要素は，図１に示す磁場シールドル
ーム１の外部に配置することが望ましい。

【００６７】図２１は，本発明の第６の実施例に於ける
磁場センサの配置と人体の頭部との位置関係を説明する
図である。図２１は，脳磁場の計測範囲（１７５ｍｍ×
１７５ｍｍ）を，磁場センサが配置される計測点３０に
より示している。被検体２の左側頭から脳磁場を計測す
る場合の計測範囲を上部に示し，下部に右側頭から脳磁
場計測する場合の計測範囲を示す。図２１では，聴性誘
発される脳磁場を計測するので，耳のやや上部にも計測
点３０を配置する。

【００６８】図２２は，本発明の第６の実施例に於ける
脳磁場の計測の手順を説明する流れ図である。最初に左
側頭から脳磁場を計測するため，被検体２の頭部とデュ
ワ３との位置合わせを行ない，最初に右耳刺激（対側側
の音刺激）を８分間行ない，連続して左耳刺激（同側側
の音刺激）を８分間行なう。続いて，右側頭から脳磁場
を計測するために，被検体２は左肩を下にしてベッド２
に横たわり，位置合わせを行なう。位置合わせ終了後，
左耳刺激（対側側の音刺激）を８分間行ない，連続して
右耳刺激（同側側の音刺激）を８分間行なう。

【００６９】図２３は，本発明の第６の実施例に於い
て，健常者の脳磁場の計測結果の表示例と，最大ベクト
ルの抽出する原理を説明する図である。図２３の上部
に，左耳刺激した場合の脳磁場の計測結果例を示し，図
２３の下部に右耳刺激した場合の脳磁場の計測結果を示
す。計測された脳磁場の波形２１２，２１４は７０回の
加算処理を行った後に，音刺激の前の時間帯で基線補正
を行なった６４チャネル全てのチャネルの計測された脳
磁場の波形を重ねて表示している。

【００７０】電流アローマップ２１１，２１３は，ぞれ
ぞれ，計測された脳磁場の波形２１２，２１４の最大ピ
ークが出現する時刻（Ｎ１００ｍと呼ばれるピーク）で
の電流アローマップを示す。電流アローマップ２１１内
の最大電流アローから強度Ｉcontralateral（対側側刺
激の強度）（以下，Ｉ_ｃと略記する）と最大電流アロー
の傾きθ１とを計算する。同様にして，電流アローマッ
プ２１３内の最大電流アローから強度Ｉipsilateral
（同側側刺激の強度）（以下，Ｉ_ｉと略記する）と最大
電流アローの傾きθ２とを計算する。これらの最大電流
アローの強度比Ｉ _ｉ／Ｉ_ｃ（同側側刺激の強度／対側側
刺激の強度）と位相差｜Δθ｜＝｜θ１−θ２｜とを求
める計算を行なう。右側頭から脳磁場を計測する場合に
も，Ｉ_ｃ（対側側刺激の強度）は，左耳刺激による脳磁
場の波形のピーク時の最大電流アローを用い，Ｉ_ｉ（同
側側刺激の強度）は，右耳刺激による脳磁場の波形のピ
ーク時の最大電流アローを用いる。

【００７１】図２４は，本発明の第６の実施例に於いて
得られた最大電流ベクトルの強度比Ｉ_ｉ／Ｉ_ｃ（同側側
刺激の強度／対側側刺激の強度）の結果例を示す図であ
る。

【００７２】図２５は，本発明の第６の実施例に於いて
得られた位相差｜Δθ｜の結果例を示す図である。

【００７３】図２６は，本発明の第６の実施例に於いて
得られたＮ１００ｍが出現した時間（潜時）を各側頭毎
にまとめた結果例を示す図である。

【００７４】健常者４例，右半球梗塞患者５例，慢性め
まいの患者２例，脳血流に障害が観測されたモヤモヤ病
患者２例についての解析結果を，図２４〜図２６を使用
して説明する。但し，右半球梗塞患者のうち５例中２例
については，梗塞範囲が広域であっため右側頭での脳磁
波形が検出されなかった。

【００７５】図２４に示すように，健常者では両側頭共
に電流比が１より小さく，対側側の電流アロー強度が強
いことが分かる。一方，右半球梗塞（右側脳梗塞）の１
例と慢性めまいの患者１例では，電流比が１より大きい
症例があることが分かる。モヤモヤ病の患者２例は，健
常者と変わらず電流比が１以下であった。

【００７６】図２５に示すように，健常者では両側頭共
に位相差は２０度以下であり，最大ベクトルの方向が良
く一致していることが分かる。一方，右側脳梗塞患者の
６例中３例で２０度を超える位相差が観測され，慢性め
まいの患者１例とモヤモヤ病患者１例でも，位相差が他
のケースと大幅に異なる結果が得られた。

【００７７】図２６に示すように，健常者の右側頭で
は，対側音刺激（左耳刺激）の潜時が有意に短い結果と
なっているが，左側頭では，特に何れにも有意性が認め
られない。右半球梗塞の左側側頭では健常者に比べて，
６例中３例で同側側の音刺激（左耳刺激）で潜時が延長
する結果が見られた。右半球梗塞の右側側頭では，健常
者に比べて，対側側の音刺激（左耳刺激）が延長する傾
向が見られ，同側側の音刺激（右耳刺激）に於いても，
３例中１例では延長していることが分かった。

【００７８】慢性めまいの患者では，２例中１例で左側
頭の潜時が延長する結果が得られた。一方，モヤモヤ病
の患者では，両側頭共に潜時に健常者との有意な差は認
められなかった。

【００７９】なお，詳細な説明は省略するが，生体に印
加する光，音等の刺激信号の印加の開始に同期する同期
信号を，脳磁場の信号と対として収集し，複数方向から
計測された脳磁場の波形の時間軸が，同期信号が収集さ
れた時点で共通の原点（ｔ＝０）を持つように，時間軸
の変換，即ち，時間軸の平行移動する変換を行なう。次
に，第１の実施例〜第４の実施例と同様にして，ベクト
ル強度Ｉ（ｘ，ｙ，ｔ）とその位相θ（ｘ，ｙ，ｔ）を
求め，最大ベクトルの時間−強度プロット（ｔ−Ｉ
_ｍａｘ），及び時間−位相プロット（ｔ−θ），あるい
は，全てのチャンネルについての時間−位相・強度プロ
ット（ｔ−θ・Ｉ）を求め，表示できる。

【００８０】以上，説明した本発明の各実施例では，心
磁場，又は，脳磁場の計測された磁場波形から得られる
ベクトル強度とその位相を使用するので，従来技術のよ
うに生体部位の各時刻に於ける状態を表わす多数の図
（マップ）を用いて生体現象を解析することなく，従来
技術で使用されていた数よりもはるかに少数の図（マッ
プ）を用いて，心臓や脳に於ける生体の電気生理学的な
興奮の時間変化を詳細に把握できる。

【００８１】

【発明の効果】本発明の生体磁場計測装置によれば，ベ
クトル強度とその位相を使用するので，逆問題を解いて
ダイポール（磁場源）推定することなく，また，多数の
等磁場線図を表示することなく，電気生理学的な興奮の
伝播過程を定量化でき，個人毎の疾患，異常を客観的，
定量的に把握できる。

【図面の簡単な説明】

【図１】本発明の第１の実施例の生体磁場計測装置の概
略構成を説明する斜視図。

【図２】本発明の第１の実施例に於ける磁場センサの配
置を説明する斜視図。

【図３】本発明の第１の実施例に於ける生体磁場の法線
成分を検出する磁場センサの構成を説明する斜視図。

【図４】本発明の第１の実施例に於ける磁場センサの配
置と人体の胸部の正面，背面との位置関係を説明する
図。

【図５】本発明の第１の実施例に於ける磁場センサの配
置と人体の胸部との位置を合わせを行なう方法を説明す
る斜視図。

【図６】本発明の第１の実施例に於ける生体磁場の計
測，計測された信号の解析の流れを説明する図。

【図７】本発明の第１の実施例に於いて，健常者を被検
体として，正面，及び背面から計測された心磁場波形
と，心磁場波形の計測と同時に計測された心電図波形の
例を示す図。

【図８】本発明の第１の実施例に於いて，正面から計測
された健常者の心磁場波形のＱＲＳ波の開始点から３０
ｍｓ経過後の時間での電流アロー図，及び最大電流ベク
トルを示す図。

【図９】本発明の第１の実施例に於いて，最大ベクトル
の位相の基準を示す図。

【図１０】本発明の第１の実施例に於いて，正面，及び
背面から計測された健常者の心磁場波形のＱＲＳ波の開
始点から２００ｍｓまでの時間帯での時間−強度プロッ
ト（ｔ−Ｉ_ｍａｘx）と，時間−位相プロット（ｔ−
θ）の表示例を示す図。

【図１１】本発明の第１の実施例に於いて，正面，及び
背面から計測された右脚ブロックの患者の心磁波形のＱ
ＲＳ波の開始点から２００ｍｓまでの時間帯での時間−
強度プロット（t−Imax）と，時間−位相プロット（t−
θ）の表示例を示す図。

【図１２】本発明の第２の実施例に於いて６４チャンネ
ルの全ての各チャンネルの各時間でのベクトルの時間−
位相・強度プロットを行なう場合の，生体磁場の計測，
計測された信号の解析の流れを説明する図。

【図１３】本発明の第２の実施例に於いて，健常者の正
面から計測された心磁場波形のＱＲＳの開始点から２０
０ｍｓまでの時間帯での６４チャンネルの全チャンネル
のベクトルの時間−位相強度プロットの表示例を示す
図。

【図１４】本発明の第２の実施例に於いて，健常者の背
面から計測された心磁場波形のＱＲＳの開始点から２０
０ｍｓまでの時間帯での６４チャンネルの全チャンネル
のベクトルの時間−位相強度プロットの表示例を示す
図。

【図１５】本発明の第２の実施例に於いて，正面から計
測された右脚ブロックの患者の心磁波形のＱＲＳ波の開
始点から２００ｍｓまでの時間帯での６４チャンネルの
全チャンネルのベクトルの時間−位相強度プロットの表
示例を示す図。

【図１６】本発明の第２の実施例に於いて，背面から計
測された右脚ブロックの患者の心磁波形のＱＲＳ波の開
始点から２００ｍｓまでの時間帯での６４チャンネルの
全チャンネルのベクトルの時間−位相強度プロットの表
示例を示す図。

【図１７】本発明の第３の実施例に於いて使用する生体
磁場の接線成分Ｂ_ｘ，及びＢ_ｙ成分を検出する磁場セン
サの構成例の概略を説明する図。

【図１８】本発明の第４の実施例に於いて使用する，コ
イル全体として円の形状を持つ微分コイルの例を示す
図。

【図１９】本発明の第５の実施例に於いて使用する，正
面，及び背面の方向から同時に心磁場計測を行なう生体
磁場計測装置の構成例の概略を説明する図。

【図２０】本発明の第６の実施例に於いて聴性誘発され
た脳磁場を計測する生体磁場計測装置の構成例を説明す
る斜視図。

【図２１】本発明の第６の実施例に於ける磁場センサの
配置と人体の頭部との位置関係を説明する図。

【図２２】本発明の第６の実施例に於ける脳磁場の計測
の手順を説明する流れ図。

【図２３】本発明の第６の実施例に於いて，健常者の脳
磁場の計測結果の表示例と，最大ベクトルの抽出する原
理を説明する図。

【図２４】本発明の第６の実施例に於いて得られた最大
電流ベクトルの強度比Ｉ_ｉ／Ｉ_ｃ（同側側刺激の強度／
対側側刺激の強度）の結果例を示す図。

【図２５】本発明の第６の実施例に於いて得られた位相
差｜Δθ｜の結果例を示す図。

【図２６】本発明の第６の実施例に於いて得られたＮ１
００ｍが出現した時間（潜時）を各側頭毎にまとめた結
果例を示す図。

【符号の説明】

１…磁場シールドルーム，２…被検体，３…デュワ，
３’−１…上部デュワ，３’−２…下部デュワ，４，
４’…ベッド，５，５’…ガントリ，６…心電計用電
極，７…心電計，８…ＦＬＬ回路，９…アンプフィルタ
ー回路，１０…データ収録解析装置，１１…ディスプレ
イ（表示装置），１２…レール，１３…滑車，２０，２
０−ｉ（ｉ＝１１，１２，…，１８；ｉ＝２１，２２，
…，２８；ｉ＝３１，３２，…，３８；ｉ＝４１，４
２，…，４８；ｉ＝５１，５２，…，５８；ｉ＝６１，
６２，…，６８；ｉ＝７１，７２，…，７８；ｉ＝８
１，８２，…，８８），２０’，２０”…磁場センサ，
２１，２１−１’，２１−２’２１”…ＳＱＵＩＤ，２
２，２２’−１，２２’−２…検出コイル，２３，２
３’−１，２３’−２…参照コイル，３０…計測点，３
１，３１’…計測基準点，４０…被検体の基準点，４１
…ｘ方向レーザ光源，４２…ｙ方向レーザ光源，４３…
ｘ軸ライン形成用ビーム，４４…ｙ軸ライン形成用ビー
ム，４５…ｘｚ標識（マーク），４６…ｙｚ標識（マー
ク），４９…交差線，５０…微分コイル，１８１…アダ
プタ，１８２…エアーチューブ，１８３…音刺激装置，
１８４…同期信号，２１１，２１３…電流アローマッ
プ，２１２，２１４…計測された脳磁場の波形。

─────────────────────────────────────────────────────

【手続補正書】

【提出日】平成１３年１０月３日（２００１．１０．
３）

【手続補正１】

【補正対象書類名】明細書

【補正対象項目名】特許請求の範囲

【補正方法】変更

【補正内容】

【特許請求の範囲】

───────────────────────────────────────────────────── フロントページの続き (72)発明者 神鳥 明彦 東京都国分寺市東恋ケ窪一丁目280番地 株式会社日立製作所中央研究所内 (72)発明者 鈴木 大介 東京都国分寺市東恋ケ窪一丁目280番地 株式会社日立製作所中央研究所内 (72)発明者 横澤 宏一 東京都国分寺市東恋ケ窪一丁目280番地 株式会社日立製作所中央研究所内 Ｆターム(参考） 2G017 AA08 AD32 BA15 4C027 AA00 AA02 AA03 AA10 CC00 DD02 FF03 GG16 HH11

Claims (22)

【特許請求の範囲】
1. 【請求項１】生体から発生する生体磁場を計測する複数
   のＳＱＵＩＤ磁束計と，前記生体磁場以外の周期的に発
   生する生体信号を計測して収集する生体信号計測装置
   と，複数方向から対として同時に計測された前記生体磁
   場の信号及び前記生体信号の演算処理を行なう演算処理
   装置と，演算処理の結果を表示する表示装置とを有する
   ことを特徴とする生体磁場計測装置。
2. 【請求項２】請求項１に記載の生体磁場計測装置に於い
   て，前記演算処理装置は，時間変数をｔとして，複数方
   向から計測された前記生体信号の波形の時間軸が共通原
   点（ｔ＝０）を持つように，前記各生体信号の波形の時
   間軸の変換を行ない，前記生体信号と対をなす前記生体
   磁場の信号の波形の時間軸に対して前記変換と同一の変
   換を行ない，前記表示装置に，前記共通原点（ｔ＝０）
   をもつ前記生体磁場の信号の演算処理から求めた結果が
   表示されることを特徴とする生体磁場計測装置。
3. 【請求項３】請求項１に記載の生体磁場計測装置に於い
   て，前記生体磁場及び前記生体信号が前記生体の胸面及
   び背面の２方向から計測されることを特徴とする生体磁
   場計測装置。
4. 【請求項４】請求項１に記載の生体磁場計測装置に於い
   て，前記生体信号計測装置が，心電計，心音計，脈波
   計，脳波計の何れかであり，前記生体信号が，心電波
   形，心音波形，脈波形，脳波形の何れかであることを特
   徴とする生体磁場計測装置。
5. 【請求項５】生体から発生する生体磁場を計測する複数
   のＳＱＵＩＤ磁束計と，前記生体磁場以外の周期的に発
   生する生体信号を計測して収集する生体信号計測装置
   と，複数方向から対として同時に計測された前記生体磁
   場の信号及び前記生体信号の演算処理を行なう演算処理
   装置と，演算処理の結果を表示する表示装置とを有し，
   前記演算処理装置は，時間変数をｔとして，複数方向で
   計測された前記生体信号を用いて，複数方向で計測され
   た前記生体磁場の信号が共通原点（ｔ＝０）を持つよう
   に時間軸の変換を行ない，前記生体の表面に接する面に
   垂直な方向をｚ方向，該ｚ方向に直交し前記生体の表面
   に接する面に水平な方向をｘ方向及びｙ方向とし，前記
   共通原点（ｔ＝０）をもつ前記生体磁場の信号を用い
   て，前記生体磁場の各計測点（ｘ，ｙ）で２次元磁場ベ
   クトルの大きさ及び／又はｘｙ面での方向を表わす角度
   を求め，前記表示装置に，前記２次元磁場ベクトルの大
   きさ及び／又は前記角度の，前記共通原点（ｔ＝０）を
   原点とする時間変化が表示されることを特徴とする生体
   磁場計測装置。
6. 【請求項６】請求項５に記載の生体磁場計測装置に於い
   て，前記各ＳＱＵＩＤ磁束計は前記ｚ方向の磁場成分
   （Ｂ_ｚ（ｘ，ｙ，ｔ））を計測し，前記演算処理装置は
   前記２次元磁場ベクトルを√｛（∂Ｂ_ｚ（ｘ，ｙ，ｔ）
   ／∂ｘ）^２＋（∂Ｂ_ｚ（ｘ，ｙ，ｔ）／∂ｙ）^２｝から
   求めることを特徴とする生体磁場計測装置。
7. 【請求項７】請求項５に記載の生体磁場計測装置に於い
   て，前記各ＳＱＵＩＤ磁束計は前記ｚ方向の磁場成分
   ｛（Ｂ_ｚ（ｘ，ｙ，ｔ）｝のｘ方向に関する微分｛∂Ｂ
   _ｚ（ｘ，ｙ，ｔ）／∂ｘ｝及びｙ方向に関する微分｛∂
   Ｂ_ｚ（ｘ，ｙ，ｔ）／∂ｙ｝を計測し，前記演算処理装
   置は前記２次元磁場ベクトルを√｛（∂Ｂ_ｚ（ｘ，ｙ，
   ｔ）／∂ｘ）^２＋（∂Ｂ_ｚ（ｘ，ｙ，ｔ）／∂ｙ）^２｝
   から求めることを特徴とする生体磁場計測装置。
8. 【請求項８】請求項５に記載の生体磁場計測装置に於い
   て，前記各ＳＱＵＩＤ磁束計は前記ｘ方向の磁場成分
   ｛Ｂ_ｘ（ｘ，ｙ，ｔ）｝及び前記ｙ方向の磁場成分｛Ｂ
   _ｙ（ｘ，ｙ，ｔ）｝を計測し，前記演算処理装置は前記
   ２次元磁場ベクトルを√｛（Ｂ _ｘ（ｘ，ｙ，ｔ））^２＋
   （Ｂ_ｙ（ｘ，ｙ，ｔ））^２｝から求めることを特徴とす
   る生体磁場計測装置。
9. 【請求項９】請求項５に記載の生体磁場計測装置に於い
   て，前記演算処理装置は，前記生体磁場を計測した各方
   向について，前記共通原点（ｔ＝０）をもつ前記生体磁
   場の信号から前記各計測点（ｘ，ｙ）に於ける前記２次
   元磁場ベクトルを求め，複数の前記計測点（ｘ，ｙ）に
   於ける前記２次元磁場ベクトルのうち最大の前記２次元
   磁場ベクトルの大きさ及び／又はｘｙ面での方向を表わ
   す角度を，前記共通原点（ｔ＝０）を原点とする時間軸
   の各点で求め，前記表示装置に，前記最大の前記２次元
   磁場ベクトルの大きさ及び／又は角度の時間変化が，前
   記生体磁場を計測した前記各方向について表示されるこ
   とを特徴とする生体磁場計測装置。
10. 【請求項１０】請求項９に記載の生体磁場計測装置に於
    いて，前記表示装置に，前記最大の前記２次元磁場ベク
    トルの大きさ及び／又は角度の時間変化が，前記生体磁
    場を計測した前記各方向について同一の時間軸を用いて
    表示されることを特徴とする生体磁場計測装置。
11. 【請求項１１】請求項９に記載の生体磁場計測装置に於
    いて，前記表示装置に，前記最大の前記２次元磁場ベク
    トルの大きさ及び／又は角度の時間変化が，前記生体磁
    場を計測した前記各方向について異なるカラーで同一の
    時間軸を用いて表示されることを特徴とする生体磁場計
    測装置。
12. 【請求項１２】請求項５に記載の生体磁場計測装置に於
    いて，前記演算処理装置は，前記生体磁場を計測した各
    方向について，前記共通原点（ｔ＝０）をもつ前記生体
    磁場の信号から前記各計測点（ｘ，ｙ）に於ける前記２
    次元磁場ベクトルの大きさ及びｘｙ面での方向を表わす
    角度を，前記共通原点（ｔ＝０）を原点とする時間軸の
    各点で求め，前記表示装置に，前記各計測点（ｘ，ｙ）
    に於ける前記２次元磁場ベクトルの大きさ及び角度の時
    間変化が，前記生体磁場を計測した前記各方向毎に，前
    記２次元磁場ベクトルの大きさをプロット点の大きさに
    比例させて，又は，カラーで区別して表示されることを
    特徴とする生体磁場計測装置。
13. 【請求項１３】生体の心臓から発生する生体磁場の信号
    を計測する複数のＳＱＵＩＤ磁束計と，生体信号とし
    て，心電波形，心音波形，脈波形の何れかを計測する装
    置と，前記生体の胸面及び背面の２方向から対として同
    時に計測された前記生体磁場の信号及び前記生体信号の
    演算処理を行なう演算処理装置と，演算処理の結果を表
    示する表示装置とを有し，前記演算処理装置は，前記生
    体の表面に接する面に垂直な方向をｚ方向，該ｚ方向に
    直交し前記生体の表面に接する面に水平な方向をｘ方向
    及びｙ方向とし，時間変数をｔとして，前記２方向から
    計測された前記生体信号の波形の時間軸が共通原点（ｔ
    ＝０）を持つように，前記２方向から計測された前記生
    体信号の波形の時間軸の変換を行ない，前記生体信号と
    対をなす前記生体磁場の信号の波形の時間軸に対して前
    記変換と同一の変換を行ない，前記生体磁場を計測した
    前記２方向について，前記共通原点（ｔ＝０）をもつ前
    記生体磁場の信号から前記生体磁場の各計測点（ｘ，
    ｙ）に於ける２次元磁場ベクトルを求め，複数の前記計
    測点（ｘ，ｙ）に於ける前記２次元磁場ベクトルのうち
    最大の前記２次元磁場ベクトルの大きさ及び／又はｘｙ
    面での方向を表わす角度を，前記共通原点（ｔ＝０）を
    原点とする時間軸の各点で求め，前記表示装置に，前記
    最大の前記２次元磁場ベクトルの大きさ及び／又は角度
    の時間変化が，前記生体磁場を計測した前記２方向につ
    いて表示されることを特徴とする生体磁場計測装置。
14. 【請求項１４】請求項１３に記載の生体磁場計測装置に
    於いて，前記角度が０°となる方向を前記生体の体軸に
    垂直な左水平方向とし，前記角度が０°から１８０°の
    範囲，及び，０°から−１８０°の範囲で表示されるこ
    とを特徴とする生体磁場計測装置。
15. 【請求項１５】請求項１３に記載の生体磁場計測装置に
    於いて，前記心臓の脚ブロック伝導障害に関する情報を
    検出することを特徴とする生体磁場計測装置。
16. 【請求項１６】生体の胸部又は頭部から発生する生体磁
    場の信号を計測する複数のＳＱＵＩＤ磁束計と，前記生
    体磁場の信号の演算処理を行なう演算処理装置と，演算
    処理の結果を表示する表示装置とを有し，前記演算処理
    装置は，前記生体の表面に接する面に垂直な方向をｚ方
    向，該ｚ方向に直交し前記生体の表面に接する面に水平
    な方向をｘ方向及びｙ方向として，前記生体磁場の信号
    から前記生体磁場の各計測点（ｘ，ｙ）に於ける２次元
    磁場ベクトルを求め，複数の前記計測点（ｘ，ｙ）に於
    ける前記２次元磁場ベクトルのうち最大の前記２次元磁
    場ベクトルの大きさ及び／又はｘｙ面での方向を表わす
    角度を，前記生体磁場が計測された複数の時点で求め，
    前記表示装置に，前記最大の前記２次元磁場ベクトルの
    大きさ及び／又は角度の時間変化が表示されることを特
    徴とする生体磁場計測装置。
17. 【請求項１７】生体の胸部又は頭部から発生する生体磁
    場の信号を計測する複数のＳＱＵＩＤ磁束計と，前記生
    体磁場の信号の演算処理を行なう演算処理装置と，演算
    処理の結果を表示する表示装置とを有し，前記演算処理
    装置は，前記生体の表面に接する面に垂直な方向をｚ方
    向，該ｚ方向に直交し前記生体の表面に接する面に水平
    な方向をｘ方向及びｙ方向として，前記生体磁場の信号
    から前記生体磁場の各計測点（ｘ，ｙ）に於ける２次元
    磁場ベクトルを求め，前記各計測点（ｘ，ｙ）に於ける
    前記２次元磁場ベクトルの大きさ及び角度を，前記生体
    磁場が計測された複数の時点で求め，前記表示装置に，
    前記各計測点（ｘ，ｙ）に於ける前記２次元磁場ベクト
    ルの大きさ及びｘｙ面での方向を表わす角度の時間変化
    が，前記２次元磁場ベクトルの大きさをプロット点の大
    きさに比例させて，又は，カラーで区別して表示される
    ことを特徴とする生体磁場計測装置。
18. 【請求項１８】生体の頭部から発生する生体磁場を計測
    する複数のＳＱＵＩＤ磁束計と，前記生体を刺激する信
    号と前記刺激する信号の発生の開始に同期する同期信号
    を発生する刺激装置と，前記生体の頭部の複数方向から
    対として計測された前記生体磁場の信号及び前記同期信
    号の演算処理を行なう演算処理装置と，演算処理の結果
    を表示する表示装置とを有することを特徴とする生体磁
    場計測装置。
19. 【請求項１９】請求項１８に記載の生体磁場計測装置に
    於いて，前記演算処理装置は，時間変数をｔとして，前
    記複数方向から計測された前記生体信号の波形の時間軸
    が共通原点（ｔ＝０）を持つように，前記同期信号に基
    づいて前記複数方向から計測された前記生体信号の波形
    の時間軸の変換を行ない，前記表示装置に，前記共通原
    点（ｔ＝０）をもつ前記生体磁場の信号の演算処理から
    求めた結果が表示されることを特徴とする生体磁場計測
    装置。
20. 【請求項２０】請求項１８に記載の生体磁場計測装置に
    於いて，前記刺激装置が，電気刺激装置，音刺激装置，
    匂い刺激装置，視覚刺激装置，触覚刺激装置の何れかで
    あることを特徴とする生体磁場計測装置。
21. 【請求項２１】生体の頭部から発生する生体磁場の信号
    を計測する複数のＳＱＵＩＤ磁束計と，前記生体を刺激
    する信号と前記刺激する信号の発生の開始に同期する同
    期信号を発生する刺激装置と，前記生体の頭部の複数方
    向から対として計測された前記生体磁場の信号及び前記
    同期信号の演算処理を行なう演算処理装置と，演算処理
    の結果を表示する表示装置とを有し，前記演算処理装置
    は，時間変数をｔとして，前記複数方向から計測された
    前記生体信号の波形の時間軸が共通原点（ｔ＝０）を持
    つように，前記同期信号に基づいて前記複数方向から計
    測された前記生体信号の波形の時間軸の変換を行ない，
    前記生体の表面に接する面に垂直な方向をｚ方向，該ｚ
    方向に直交し前記生体の表面に接する面に水平な方向を
    ｘ方向及びｙ方向とし，前記共通原点（ｔ＝０）をもつ
    前記生体磁場の信号を用いて，前記生体磁場の各計測点
    （ｘ，ｙ）で２次元磁場ベクトルの大きさ及び／又はｘ
    ｙ面での方向を表わす角度を求め，前記表示装置に，前
    記２次元磁場ベクトルの大きさ及び／又は前記角度の，
    前記共通原点（ｔ＝０）を原点とする時間変化が表示さ
    れることを特徴とする生体磁場計測装置。
22. 【請求項２２】生体の頭部から発生する生体磁場の信号
    を計測する複数のＳＱＵＩＤ磁束計と，前記生体を刺激
    する信号と前記刺激する信号の発生の開始に同期する同
    期信号を発生する刺激装置と，前記生体の頭部の複数方
    向から対として計測された前記生体磁場の信号及び前記
    同期信号の演算処理を行なう演算処理装置と，演算処理
    の結果を表示する表示装置とを有し，前記演算処理装置
    は，前記生体の表面に接する面に垂直な方向をｚ方向，
    該ｚ方向に直交し前記生体の表面に接する面に水平な方
    向をｘ方向及びｙ方向とし，時間変数をｔとして，前記
    複数方向から計測された前記生体信号の波形の時間軸が
    共通原点（ｔ＝０）を持つように，前記同期信号に基づ
    いて前記複数方向から計測された前記生体信号の波形の
    時間軸の変換を行ない，前記生体磁場を計測した前記複
    数方向について，前記共通原点（ｔ＝０）をもつ前記生
    体磁場の信号から前記生体磁場の各計測点（ｘ，ｙ）に
    於ける２次元磁場ベクトルを求め，複数の前記計測点
    （ｘ，ｙ）に於ける前記２次元磁場ベクトルのうち最大
    の前記２次元磁場ベクトルの大きさ及び／又はｘｙ面で
    の方向を表わす角度を，前記共通原点（ｔ＝０）を原点
    とする時間軸の各点で求め，前記表示装置に，前記最大
    の前記２次元磁場ベクトルの大きさ及び／又は角度の時
    間変化が，前記生体磁場を計測した前記複数方向につい
    て表示されることを特徴とする生体磁場計測装置。