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

Description

本発明は，生体の脳の神経活動や心臓の心筋活動等，生体内の電流などが原因で発生する生体磁場を計測するのに適した生体磁場計測方法及び生体磁場計測装置に関する。

従来から磁気センサである超伝導量子干渉素子（ＳＱＵＩＤ）を用いて，生体から発生する微弱な磁場の分布を測定し，その測定結果から，生体内部の活動電流の位置を推定し，その分布をイメージングする多チャンネルの生体磁気イメージング装置が知られている。そのような従来例は，例えば特開平４−３１９３３４号或いは特開平５−１４６４１６号等の公開公報に開示されている。
特開平４−３１９３３４号 特開平５−１４６４１６号

前記従来例は、生体磁気イメージング装置に関する動作原理に関するものであり、その開示内容には実施化する上での技術課題や解決手段が開示されていない。また、前記従来例は脳の内部に発生する生体活動電流に関するものであって、他の部位に関する具体的な開示はされていない。本発明の目的は、複数の計測位置の磁場強度から生体の磁場図表示が容易で操作の良好な生体磁場計測装置を提供することにあり、特にオフセット電圧の調整が容易な生体磁場計測装置を提供することにある。

本発明は、被検者の生体内から発する磁場を複数の位置で計測する生体磁場計測装置において、生体の心臓から発する磁場を計測する複数の磁気センサと，計測を含む処理機能項目を示す処理機能表示部と，操作項目を表示する操作領域表示部と，及び，前記磁気センサの超伝導量子干渉素子に外部磁場φを０からΔφの間隔で変えて与えた時に，該外部磁場Φに対応して発生する電圧Ｖと前記外部磁場φとで示されるφ−Ｖ特性曲線を表示する解析データ表示部とを具備する表示手段と，前記電圧Ｖの最大値と最小値からオフセット電圧を計算しオフセット電圧を調整する計算手段とを有することを特徴とする。

また、前記表示手段は、計測を含む処理機能項目を示す処理機能表示部と、計測された生体磁場とこれを処理して得られた処理磁場を表示する解析データ表示部と、及び生体磁場は処理磁場に対応した操作項目が表示される操作領域表示部からなる表示部から構成され、前記解析データ表示部には、各磁気センサに対応したチャンネルが表示され、前記操作領域表示部には、選択チャンネル項目が表示され、且つ該選択チャンネル項目は前記解析データ表示部に表示されるチャンネルと関連づけたチャンネル項目とされて構成され、且つ前記表示された選択チャンネル項目からチャンネル項目を選択する手段を有し、前記処理機能表示部の計測項目を選択すると、前記解析データ表示部に前記チャンネルを、そして選択チャンネル項目には前記チャンネル項目を表示する手段を有する構成としてもよい。

好ましくは，前記チャンネルとチャンネル項目とは，相関のある行及び列から構成される。好ましくは，前記選択チャンネル項目から選択された行もしくは列からなるチャンネル項目又は１つのチャンネル項目を他から識別する手段を設けてある。好ましくは，前記計測項目のプルダウンメニューは，「調整値ファイル
」，「全自動調整」「オフセット電圧Ｖ_off調整」，「マニュアル調整」，「調
整値ファイル」，「計測パネル」，「自動波形診断」，「ＡＦＡオフセット調整
」などのサブプルダウンメニューを有する。好ましくは，前記操作領域表示部には，開始と中止を指示する計測項目が表示され，開始指示によって計測が開始し
，中止指示によって計測が中止される。好ましくは，前記解析データ表示部に表示画面が指定された範囲外にあるときに当該チャンネルを計測に不適当なチャンネルとして識別信号を発する自動波形診断手段を設けてある。

また、前記表示部は、計測を含む処理機能項目を示す処理機能表示部と、計測された生体磁場とこれを処理して得られた処理磁場を表示する解析データ表示部とからなり、前記計測は、「データ計測」メニューの「計測パネル」サブメニューを選択して実行し、且つ前記解析データ表示部に「データ解析」メニューの「単一波形表示」、「グリッドマップ表示」、「等磁線図」及び「伝播時間図」を表示する構成としてもよい。

図１は本発明が実施される生体磁気計測装置の一実施例の概略構成を示す。環境磁気雑音の影響を除去するために，生体磁気計測装置は磁気シールドルーム１内に設置される。被検者（生体からなる被検体）２は，ベッド３に仰向け状態で計測が行われる。被検者の生体面（胸部の場合は一般に胸壁に平行な面）はベッド３の面と略平行であるとし，そしてこの面は直交座標系（ｘ，ｙ，ｚ）のｘ−ｙ平面と平行であるものとする。被検者の胸部は曲面であると共に傾いているが
，説明を簡単にするために略平行とする。被検者２の胸部の上方には，冷媒である液体Ｈeで満たされたデユワ４が配置され，該デユワは超伝導量子干渉素子（
ＳＱＵＩＤ＝Ｓｕｐｅｒｃｏｎｄｕｃｔｉｎｇ Ｑｕａｎｔｕｍ Ｉｎｔｅｒｆｅｒｅｎｃｅ Ｄｅｖｉｃｅ）とそのＳＱＵＩＤに接続された検出コイルとを含む複数個の磁気センサを収容している。液体Ｈｅは磁気シールド１の外部にある自動補給装置５から連続的に補給される。

磁気センサの出力は，被検者２から発生して検出コイルにより検出される生体磁場の強度（磁束密度と考えることもできる）と特定の関係をもつ電圧を出力し
，その出力がＦＬＬ（Ｆｌｕｘ Ｌｏｃｋｅｄ ｌｏｏｐ）回路６に入力される
。このＦＬＬ回路６は，ＳＱＵＩＤの出力を一定に保つように，ＳＱＵＩＤに入力された生体磁場（生体磁気）の変化を帰還コイルを介してキャンセルする（これを磁場ロックと呼ぶ）。その帰還コイルに流した電流を電圧に変換して，生体
磁場信号の変化と特定の関係にある電圧出力を得ることができる。このように帰還コイルを介して検出する方式を取っているので，微弱の磁場を高感度に検出できる。上記出力電圧は増幅器・フイルター・増幅器（ＡＦＡ）７に入力され，その出力はサンプリングされて，Ａ／Ｄ変換され，計算機８に取り込まれる。

計算機８はパーソナルコンピュータからなり，８−１はそのデイスプレイ部，８−２はキ−ボ−ド，そして８−３はマウスを示す。マウス８−３は画面上でカーソルを移動させて処理対象を選択するのに用いられる。この操作はキーボードを操作しても実行できる。ＡＦＡ７の入力ゲイン（Ｉ_gain）及び出力ゲイン（Ｏ
_gain）は調整可能であり，また，ＡＦＡ７は第１の基準周波数以下の周波数信号
を通過させるローパスフイルタ（ＬＰＦ），第１の基準周波数よりも低い第２の
基準周波数以上の周波数信号を通過させるハイパスフイルタ（ＨＰＦ）及び商用電源周波数をカットするノッチフィルタ（ＢＥＦ）を含む。計算機８は各種の処理を行うことができ，その処理結果はデイスプレイ部８−１に表示され得る。なお，図１で示す前記計算機８は一実施例を示したものであり，これに限定されるものではない。例えば，タッチパネルを備えたデイスプレイを備えたものや，マウスに変えて他の座標指示装置，例えばトラックボールやジョスティック等を使用したものでも良い。また場合によっては公衆電話回線を介して接続される計算機でも良い。

ＳＱＵＩＤとしては，例えば１例として直流ＳＱＵＩＤが用いられる。ＳＱＵＩＤに外部磁場が与えられたときに，それに対応する電圧（Ｖ）が発生するようにＳＱＵＩＤには直流バイアス電流（Ｉ_bias）が流される。その外部磁場を磁束
Φで表すと，ＶのΦに対する特性曲線，即ちΦ−Ｖ特性曲線は周期関数で与えられる。計測に当っては，それに先立って，ＦＬＬ回路６のオフセット電圧（Ｖ_OF
F）を調整してΦ−Ｖ特性曲線の直流電圧をゼロレベルにする操作が行われる。
更に，ＡＦＡ７の入力がゼロのときその出力がゼロとなるようにＡＦＡ７のオフセット電圧（Ａ_OFF）の調整が行われる。ＳＱＵＩＤに外部から大きな磁場がか
かると，その磁場がＳＱＵＩＤのデバイス上にトラップされ，その正規の動作がなされなくなる。その場合は，ＳＱＵＩＤを加熱して一旦常伝導状態にし，そしてその後その加熱を止めて，トラップされた磁場を取り除くことができる。その場合の，ＳＱＵＩＤの加熱操作をヒートフラッシュと呼ぶ。

図２は磁気センサの配置構成を示す。磁気センサの検出コイルには生体磁場の接線成分（生体面，即ちｘ−ｙ平面に略平行な成分）を検出するコイルと生体磁場の法線成分（生体面，即ちｘ−ｙ平面に直交する成分）を検出するコイルがある。生体磁場の接線成分を検出するコイルとしては，コイル面がｘ方向及びｙ方向をそれぞれ向いた２つのコイルが用いられ，また，生体磁場の法線成分を検出するコイルとしてはコイル面がｚ方向を向いたコイルが用いられる。複数個の磁気センサ２０−１〜２０−８，２１−１〜２１−８，２２−１〜２２−８，２３−１〜２３−８，２４−１〜２４−８，２５−１〜２５−８，２６−１〜２６−８及び２７−１〜２７−９は，図２に示すように，生体面，即ちｘ−ｙ平面と略平行な面上にマトリックス状に配置される。磁気センサの数は任意であって良いが，図２では，磁気センサのマトリックスは８行８列からなっているから，磁気センサの数は８×８＝６４である。各磁気センサは，図２に示すように，その長手方向が生体面，即ちｘ−ｙ平面に対して垂直な方向（ｚ方向）と一致するように配置される。なお，この一実施例では，ベッド面とセンサーのＸ−Ｙ面とを平行にしているが，測定精度を高めるには体に接近させる方が良く，傾斜させることができる。但し，被検者である人体は常に動いているので，人体に密着させるとこの動きが検出部を動かし，かえって高精度の検出が困難となる。

図３は磁気センサの各々の，生体磁場の法線成分Ｂ_zを検出するセンサの構成
を示す。同図において，超伝導線（Ｎｉ−Ｔｉ線）で作られたコイルはそのコイル面がｚ方向を向くように配置される。このコイルは互いに逆向きの２つのコイル１０及び１１の組み合わせからなり，被検者２に近い方のコイル１０は検出コイルとされ，遠い方のコイル１１は外部磁場雑音を検出する参照コイルとされる
。外部磁場雑音は被検者よりも遠い信号源から生じており，従って，その雑音信号は検出コイル１０及び参照コイル１１の両方によって検出される。一方，被検者からの磁場信号は微弱であり，従って，その生体磁場信号は検出コイル１０によって検出されるが，参照コイル１１はその生体磁場信号にほとんど感応しない
。このため，検出コイル１０は生体磁場信号と外部磁場雑音信号を検出し，参照コイル１１は外部磁場雑音信号を検出するから，両コイルで検出された信号の差をとることによりＳ／Ｎ比の高い生体磁場の計測が可能となる。これらのコイルはＳＱＵＩＤ１２を実装した実装基板の超伝導線を介してＳＱＵＩＤの入力コイルに接続され，これによって，検出された生体磁場信号の法線方向の成分Ｂ_zが
ＳＱＵＩＤに伝達される。

図４は磁気センサの各々の，生体磁場の接線成分Ｂ_x及びＢ_yを検出するセンサの構成を示す。同図において，接線方向の生体磁場成分検出用のセンサでは平面コイルが用いられる。即ち，検出コイル１０’及び１０”並びに参照コイル１１
’及び１１”は平面コイルからなり，これらは互いにｚ方向において間隔づけられている第１及び第２の平面にそれぞれ配置される。これらのコイルは法線成分用と同様にＳＱＵＩＤ１２’及び１２”の実装基板の入力コイルに接続される。４角柱の互いに直交する２面に，これらのＢ_x成分検出用のセンサ１３及びＢy成分検出用のセンサ１４が貼付けられ，これによってＢ_x成分及びＢ_y成分を検出し得るセンサが形成される。接線成分Ｂx，Ｂyについては，これを図４に示す磁気センサを用いて検出する以外に，図３の磁気センサで得られた法線成分Ｂ_zをｘ
，ｙについて偏微分して求めても良い。この場合は一つの磁気センサで接線成分Ｂ_x，Ｂ_yと法線成分Ｂ_zとの両方を検出し測定できる。

図５は磁気センサと被検者２の被計測部である胸部３０との位置関係を示す。示されている点は図２に示すマトリックス上の行と列との交点，即ち被検者２の計測点，即ち計測位置を表す。これらの各計測位置をチャンネルとも呼ぶ。図からわかるように，この実施例では，被検者２の身長方向をｙ方向とし，被検者２の横方向をｘ方向としている。

図６は図５に示すそれぞれの計測位置の生体磁場の測定結果を示す。この測定結果は各計測位置に対応する磁気センサでそれぞれ検出された信号を基に上述の処理を行うことで得られた，時間的に変化する生体磁場波形をマトリックス中の対応するチャンネル毎に示すものである。この実施例では心臓の筋肉が発する磁場を検出できる位置に各チャンネルを設けたので，図６の波形は心磁波形を示している。なお，心臓の筋肉から発する磁場を計測して得られた波形を心磁波形と呼ぶ。図６に示すように，チャンネル毎に計測された計測デ−タをチャンネル毎にその位置に対応させて表示する場合，これをグリッドマップ表示と呼ぶ。図６はある健常者についての心磁の計測結果の波形を示す。ここで（ａ）は接線成分Ｂ_xの心磁波形を，（ｂ）は接線成分Ｂ_yの心磁波形を，そして（ｃ）は法線成分
Ｂ_zの心磁波形をそれぞれ示す。

図７はある健常者について特定の２チャンネルに絞って計測された接線成分Ｂ_xの心磁波形を示す。実線はあるチャンネルの，点線は他のチャンネルの心磁波
形を示す。心臓の心室が脱分極した時間帯Ｔ１つまり収縮期のＱＲＳ波でのそれぞれの波形ピークの時間がｔ_Q，ｔ_R，及びｔ_Sとしてそれぞれ示されている。ま
た，心臓の再分極過程（拡張期）であるＴ波の時間帯はＴ２として示されている
。計測されたデータを表示する場合，グリッドマップ表示の他に，行毎に又は列毎に単一チャンネルのデータを表示するようにしても良いし，全チャンネル又は複数のチャンネルのデータを重ねて表示するようにしても良い。前者を単一波形表示，後者を重ね波形表示と呼ぶ。

得られる心磁波形データについては，その平均化（アベレージング）処理を行ったり，生体磁場に基づいて等磁線図や時間積分図等を作成したりして，その結果を表示できる。例えば，図７を参照するに，ＱＲＳ波の立上り部分がしきい値ＳＬと一致した時点から予め定められた時間（ｔ_OFF）だけ遡り，その遡った時
点ｔ１から予め定められた時間Ｔ３だけ経過した時点ｔ２までの間のデータを予め定められた回数だけ加算することが行われる。これがアベレージングで，予め定められる時間Ｔ３をアベレージング時間，ｔ_OFFをオフセット時間と呼ぶ。心
磁波形データは予め定められた時間範囲に亘って積分されても良い。その時間積分値が等しい点（チャンネル）を結んで作られたマップを等時間積分図と呼ぶ。また，心磁波形信号値が等しい点を結んで作られたマップを等磁線図と呼ぶ。なお，各チャンネルは粗く設定されているので，予め等磁線の間隔つまり磁場強度差を設定して各チャンネル間を直線補間して等磁線を描くことにより，より診断に適した図を作ることができる。図７に示す心磁波形について言えば，時点ｔ１からＱＲＳ波のピーク位置時点（ｔ_R時点）までの時間を伝播時間と呼び，その
伝播時間が等しい点を結んで作られたマップを等伝播時間図と呼ぶ。しきい値ＳＬの設定レベルは変更可能である。ｔ１時点については，これを，ＱＲＳ波の立上り部分がしきい値ＳＬと一致した時点を基準として決定しているが，ＱＲＳ波の立ち下がり部分がしきい値ＳＬと一致した時点を基準にして決定しても良い。ｔ１時点は更に，ＱＲＳ波のピーク位置時点（ｔ_R時点）を検出し，この時点を
基準として決定されても良い。計測された生体磁場信号は生体内の電気的生理現象によって生成されるもので，その発生源は電流双極子モデルによって近似される。その磁場発生源の電流双極子は等磁線図上に合成表示され，それを磁場源表示と呼ぶ。

被検者の登録から，その登録された被検者のデータ計測を行って，その計測されたデータの解析を行うまでの一連の操作はデイスプレイ８−１に表示される表示画面を見ながら行われる。このため，その一連の操作の説明に先立ってまずその表示画面のレイアウトを説明する。

図８は，図１のデイスプレイ８−１に表示される表示画面の基本的なレイアウトを示す。表示画面の上方部は上から順番に配置されたタイトルバー部８０１，メニューバー部８０２及びアイコンが配置されているツールバー部８０３によって占められる。上記各部は表示領域やエリアと考えることもできる。これらの配置は他の処理目的，例えば，被検者の登録や読み出し，磁場の計測，計測データの解析のための処理，などにおいてもその表示画面で共通して表示される。これにより使用し易さが増し，計測や処理の時間が短縮できる。表示画面の中央部は左から右へ順番に配置された被検者情報部８０４，線図や波形のような解析データが表示される解析データ部８０５及び操作領域部８０６によって占められている。また，下部はステータスバー８０７で占められ，これは，左側に配置された
，次の操作に関するガイドメッセージを表示するメッセージバー部８０７−１とその右側に配置された日時表示部８０７−２からなる。なお，前記メッセージバー部８０７−１と前記日時表示８０７−２を１つの表示領域としても良い。

この実施例における表示画面では，常に最上部にこのシシテムの名称が表示されるタイトルバー部８０１と，このシステムの基本的な操作を行うメニューバー部８０２，及び前記メニューバー部８０２における使用頻度の高い操作が可能なツールバー部８０３が配置されているので，使用者は表示画面が変わるたびに操作エリアを探す必要がなく，常に表示画面の上部を見れば現在動作中のシステムを知ることができる。しかも，表示画面の上部は，人が文章を読む場合を想定すると明らかなように，先ず最初に目を向ける部分であることから，このシステムの操作における基本的な事項を最上部に設けることで，自然な形で使いやすさを向上している。また，表示画面の中央部は，その中央にはこの表示画面の主体をなす解析データ表示部８０５を大きく設けることで見易さを向上するとともに，その右側にこの表示画面に特有の操作領域部８０６を設けることで，右手操作における表示画面と操作領域部の配置と同様な配置としているので，違和感なく操作を実行できる。従って，この表示画面をタッチパネル付きの表示画面に採用しても，操作領域部８０６を操作する右手が解析データ表示部８０５を邪魔することがない。また，同様に，解析データ表示部８０５の左側に確認機能しかない被検者情報部８０４を設けているので，常に患者を確認しながら操作を実行できる
。しかも，この左側の位置は，右手操作における最も遠い位置となるのでタッチパネル付きの表示画面に採用しても表示の見易さに影響をきたすことがない。

また，解析データ表示部８０５及び操作領域部８０６は被検者リスト及びその被検者のデータリストが表示されるときだけはそれらによって置き換えられる（図２４）。被検者情報部８０４には，被検者リスト画面（図２４）が表示され
ているときはその画面中の被検者リストにおいてカーソルが置かれる被検者の情報が常時表示され，また，解析データ部に線図や波形のような解析データが表示されているときは（図２５〜３４），その表示されている解析データが得られた
被検者の情報が常時表示される。これによって，表示されている解析データとその解析データが得られた被検者との関係を明確に知ることができる。このように
，このシステムの表示画面においては，メニューバー部８０２と同様に，常に，被検者情報部８０４が表示画面の定位置（左側）に表示されているので，使用者は表示画面が変わるたびに被検者情報エリアを探す必要がなく，常に表示画面の所定位置（左側）を見れば知ることができる。

タイトルバー部にはフレームの名称，具体的には，名称「Ｍｕｌｔｉｃｈａｎｎｅｌ ＭＣＧ Ｓｙｓｔｅｍ」が表示される（図２４〜３４）。操作領域部８
０６にはボタンやテキストボックスのような操作要素が配置されている。メニューバー部は操作メニューを選択する部分で，メニューは「ファイル（Ｆ）」，「
編集（Ｅ）」，「リスト（Ｌ）」，「データ計測（Ｑ）」，「データ解析（Ａ）
」及び「ヘルプ（Ｈ）」からなり，操作の順序に従って配置されている。

図９は表示画面中のメニュー部におけるそれぞれの操作メニューの内容を示し
，これらのメニューの内容はそれぞれ対応するメニューボタンをクリックしてプルダウンメニューとして表示される。このため，操作メニューを必要としないときは，前記各メニューを呼び出すためのキーワードのみメニューバー部にコンパクトに表示しているので，解析データ部や操作領域部等の各作業で必要な表示エリアを広く設定できる。そして，操作メニューが必要なときは，操作手順に従って配列された前記キーワードをメニューバー部から選択して表示して操作指示できる。この際，前記キーワードは，文字の配列（左から右）に準じて配列されているので，自然な形で操作指示できる。

「ファイル（Ｆ）」のプルダウンメニューは，ページレイアウトダイアログボ
ックス（図示せず）を開いて，ペイジレイアウトをセットする「ページレイアウト（Ｕ）」，印刷前のプレビューのための「プレビュー（Ｖ）」，データの印刷
のための「印刷（Ｐ）」及びＭｕｌｔｉｃｈａｎｎｅｌ ＭＣＧ Ｓｙｓｔｅｍ
を終了させる項目「心磁システムの終了（Ｘ）」を含む。「編集（Ｅ）」のプル
ダウンメニューは項目「削除（Ｄ）」を含む。被検者リスト画面（図２４）上の
被検者リスト中の被検者の計測データは同画面上のデータリストに表示されるが
，項目「削除（Ｄ）」はそのデータリスト中のカーソルが置かれたデータを削除
するためのもので，予めデータが選択されていなければならない。このメニューがクリックされると，削除して良いかどうかの確認を行う確認ダイアログボックスが開かれる。削除が必要な場合は，そのダイアログボックス中の「ＯＫ」ボタンがクリックされ，削除をキャンセルしたいときは「キャンセル」ボタンがクリックされる。この確認ダイアログボックスの表示及び操作により，誤って被検者の計測データを削除する誤操作を軽減できる。

「リスト（Ｌ）」のプルダウンメニューは「登録（Ｒ）」，「リスト（Ｌ）」
，「削除（Ｄ）」，「検索（Ｓ）」及び「解除（Ｘ）」の項目を含む。プルダウ
ンメニューの「登録（Ｒ）」ボタンがクリックされると，図１０に示す被検者登
録ダイアログボックスが開かれる。これは被検者に関するデータを登録するときに用いられるもので，登録できる項目は登録年月日，所定桁までの被検者のＩＤ番号，氏名，生年月日，身長，体重，性別，病気の分類情報を示す分類及び被検者コメントである。そのダイアログボックス中の「登録」ボタンがクリックされると，登録のために入力されたデータは登録されると共に，その入力ボックスは全てクリアされて再入力が可能となり，「キャンセル」ボタンがクリックされれ
ば，その全ての入力ボックスはクリアされ，そして「終了」ボタンがクリックされれば，被検者登録ダイアログボックスが閉じられる。プルダウンメニューの「
リスト（Ｌ）」がクリックされると，被検者リスト画面（図２４）が表示される
。また，プルダウンメニューの「削除（Ｄ）」は被検者リスト画面（図２４）上
の被検者リスト中のカーソルが置かれた被検者を削除するためのもので，それがクリックされると，削除前に，プルダウンメニューの「削除（Ｄ）」がクリック
されるときと同様に削除して良いかどうかの確認を行う確認ダイアログボックスが開かれ，削除が必要な場合は，そのダイアログボックス中の「ＯＫ」ボタンがクリックされ，削除をキャンセルしたいときは「キャンセル」ボタンがクリックされる。被検者が削除されると，その被検者に関するデータリスト中の全てのデータも削除される。

プルダウンメニューの「検索（Ｓ）」がクリックされると，図１１に示す検索
ダイアログボックスが表示される。被検者及びデータを検索して，被検者リスト画面に検索された被検者及びデータのみを表示する。被検者及びデータの検索対象は全被検者及び全データとする。被検者名，登録日付，性別，データの種別，計測日付，診断結果，コメント，検査技師，体位などをキーワードとして与えてデータが検索される。前記のキーワードを複数個組み合わせた複合検索を実行できる。プルダウンメニューの「解除（Ｘ）」は検索された被検者及びデータだけ
の表示を解除して全被検者及び全データを表示するのに用いられる。「データ計
測（Ｑ）」のプルダウンメニューは「調整値ファイル（Ｆ）」，「全自動調整（
Ａ）」，「Ｖ_OFF調整（Ｖ）」，「マニュアル調整（Ｍ）」，「調整値ファイル
（Ｆ）」，「計測パネル（Ｐ）」，「自動波形診断（Ｄ）」及び「ＡＦＡオフセ
ット調整（Ｏ）」の項目を含む。「調整ファイル（Ｆ）」がクリックされると，
「開く（Ｏ）」，「上書き保存（Ｓ）」及び「名前をつけて保存（Ａ）」を内容
とするサブプルダウンメニューが表示される。サブプルダウンメニュー中の「開く（Ｏ）」はシステム調整画面（図３４）を表示して，指定された調整値ファイ
ルを開いてシステムの調整値，即ちＩ_bias及びＶ_OFFをシステムに設定するのに
用いられる。「上書き保存（Ｓ）」は確認ダイアログボックスを開いて，現在の
調整値を現在開かれている調整値ファイルに上書き保存するのに用いられる。「
名前をつけて保存（Ａ）」は現在の調整値データを名称を変えて別の調整値ファ
イルに保存するのに用いられる。

プルダウンメニューの「全自動調整（Ａ）」が選択されると，システム調整画
面（図３４）が表示され，後述の図２２のフローに従ってバイアス電流Ｉ_bias及
びＦＬＬ回路６のオフセット電圧Ｖ_OFFが自動的に調整される。「Ｖ_OFF調整(Ｖ)」
が選択されると，システム調整画面（図３４）が表示され，後述の図２３のフローに従ってＦＬＬ回路６のオフセット電圧Ｖ_OFFが自動的に調整される。「マニ
ュアル調整（Ｍ）」が選択されると，図１２に示すマニュアル調整ダイアログボ
ックスが開かれる。オペレータはスクロールバーとマウスを使ってチャンネルを選択し，バイアス電流Ｉ_bias及びオフセット電圧Ｖ_OFFを変更する。入力された
値に問題なければ，「ＯＫ」ボタンをクリックしてその値の設定が行われる。「
キャンセル」ボタンがクリックされると，変更は無効となり，ダイアログボックスが閉じられる。

「計測パネル（Ｐ）」がクリックされると，グリッドマップを含むデータ計測
画面（図２５）が表示される。「自動波形診断（Ｄ）」がクリックされると，図
１３に示す自動診断ダイアログボックスが開かれる。自動波形診断は，図３４の表示画面中のΦ−Ｖ特性曲線の振幅，中央値及び周期を自動的に調べ，Φ−Ｖ特性曲線を表示する際或いは表示中の同曲線の表示を最新状態に更新する際に，指定された範囲外にあるチャンネルを計測に不適当な状態としてオペレータに通知するものである。計測に不適当な状態のチャンネルが検出された場合，エラーダイアログボックスを開いてエラーメッセージを表示してオペレータに通知するが
，当該チャンネルのΦ−Ｖ特性曲線を別色で表示して不適当な状態であることを示しても良い。振幅最小値がチェックされている場合，Φ−Ｖ特性曲線における最大値と最小値の差が振幅最小値よりも小さいときが計測に不適当な状態であるとする。中央値が指定されている場合，前記最大値と最小値との平均値の絶対値が指定された中央値よりも大きい場合が計測に不適当な状態であるとする。周期が指定されている場合，Φ−Ｖ特性曲線におけるΦの周期が第１のテキストボックス（下限値）と第２のテキストボックス（上限値）で指定された範囲外にあるときが計測に不適当な状態であるとする。振幅最小値，中央値及び周期の自動診断を有効にするには，各項目の左側のチェックボックスをマウスでクリックして×印が表示されている状態にして，対応するテキストボックスに所望する値を入力すれば良い。このようにして入力された自動診断のパラメータは「ＯＫ」ボタンを押すことによって有効となり，ダイアログボックスが閉じられる。「キャン
セル」ボタンが押されると，入力されたパラメータは無効となり，ダイアログボックスは閉じられる。「ＡＦＡオフセット調整（Ｏ）」はＡＦＡ７のオフセット
電圧ＡＯＦＦを調整するときに用いられるもので，「ＡＦＡオフセット調整（Ｏ
）」がクリックされると，単一波形表示に属するデータ計測画面（図２５）が表
示される。「データ解析（Ａ）」のプルダウンメニューは「アベレージング（Ａ
）」，「単一波形表示（Ｗ）」，「重ね波形表示（Ｍ）」，「グリッドマップ表
示（Ｇ）」，「等磁線図（Ｂ）」，「時間積分図（Ｔ）」，「伝播時間図（Ｐ）
」，「磁場源表示（Ｓ）」，「ラインモード（Ｌ）」，及び「塗潰しモード（Ｆ
）」を含む。

「アベレージング（Ａ）」がクリックされるとアベレージング画面（図２７）
が，「単一波形表示（Ｗ）」がクリックされると単一波形表示字画面（図２８）
が，「重ね波形表示（Ｍ）」がクリックされると重ね波形表示画面（図２９）が
，「グリッドマップ表示（Ｇ）」がクリックされるとグリッドマップ表示画面（
図３０）がそれぞれ表示され，メニューの左側に選択されたことを示すチェックマーク（×）が表示される。また，「等磁線図（Ｂ）」がクリックされると等磁
線図画面（図３１）が，「時間積分図（Ｔ）」がクリックされると等時間積分図
画面（図３２）が，「伝播時間図（Ｐ）」がクリックされると等伝播時間図画面
（図３３）がそれぞれ表示され，メニューの左側に選択されたことを示すチェックマーク（×）が表示される。また，「磁場源表示（Ｓ）」がクリックされると
メニューの左側に逆三角マークが表示され，等磁線図を表示する際に電流ダイポールで近似した磁場源が重ねて表示される（図なし）。等磁線図，等時間積分図
，等伝播時間図及び磁場源表示の画面は「ラインモード（Ｌ）」が選択されると
，線間が塗潰しなしで表示され，「塗潰しモード（Ｆ）」が選択されると，線間
が塗潰された状態で表示される。現在のモードを示すために，メニューの左側に選択を表すチェックマークが表示される。

「ヘルプ（Ｈ）」のプルダウンメユーは「目次（Ｃ）」，「キーワードで検索
（Ｓ）」及び「バージョン情報（Ａ）」を含み，それぞれヘルプウインドウを開
いて目次を示すこと，キーワードでトピックを検索すること，バージョンダイアログボックスを開くことのために用いられる。ツールバー部８０３には「被検者登録」（８０８），「被検者リスト」（８０９），「印刷」（８１０），「プレ
ビュー」（８１１），「システム調整」（８１２），「データ計測」（８１３）
，「データ解析」（８１４）のアイコンが配置される。これらは，その図示は省
略されているが，メニューの機能とリンクしていて，プルダウンメニューの項目のうちの使用頻度の高いものを選択できる。即ち，「被検者登録」（８０８）は
「リスト（Ｌ）」の「登録（Ｒ）」と，「被検者リスト」（８０９）は「リスト
（Ｌ）」の「リスト（Ｌ）」と，「印刷」（８１０）は「ファイル（Ｆ）」の「
印刷（Ｐ）」と，「プレビュー」（８１１）は「ファイル（Ｆ）」の「プレビュ
ー（Ｖ）」と，「システム調整」（８１２）は「データ計測（Ｑ）」の「マニュ
アル調整（Ｍ）」と，「データ計測」（８１３）は「データ計測（Ｑ）」の「計
測パネル（Ｐ）」と，そして「データ解析」（８１４）は「データ解析（Ａ）」
の「グリッドマップ表示（Ｇ）」とそれぞれ対応している。このようにアイコン
とリンクしているメニューについては，そのメニューをアイコンをクリックするだけで選択することもできる。従って，使用頻度の高い操作は，解析データ表示部に隣接したアイコンをクリックするだけで簡単に操作できるので，前記メニューバー部の操作に比べて短時間にしかも認識し易いアイコンで操作できる。なお
，前記アイコンは利用者によって選択するようにしてもよく，また，利用頻度（
回数）に準じて自動的に前記ツールバー部８０３に表示するようにしても良い。
次に，システムの調整操作を含めて，被検者の登録から，その登録された被検者のデータ計測を行って，その計測されたデータの解析を行うまでの一連の操作を図１４〜３４を参照しながら説明する。

図１６は全体の操作のフローを示すもので，計算機８の電源がＯＮにされると（Ｓ−１），オペレーテイングシステムが立ち上げられ，プログラム起動アイコ
ンがデイスプレイ部８−１に表示される（Ｓ−２）。そのアイコンの中からＭｕ
ｌｔｉｃｈａｎｎｅｌ ＭＣＧ Ｓｙｓｔｅｍのプログラムのアイコンが選択されると（Ｓ−３），図２４に示す被検者リスト画面が代わって表示される（Ｓ−
４）。この実施例に係るシステムにおいては，システム立ち上がりの初期画面と
して図２４に示す被検者リスト画面を表示する。この理由は，被検者と，該被検者の計測又は解析データの関係が極めて重要であるため，このシステムでは被検者情報をキワードとしてデータ管理していることに起因する。即ち，計測データや解析データは被検者情報がないと管理ができないためである。このため，このシステムでは，被検者リスト画面において，先ず被検者を登録又は登録されているときは被検者を特定し，次に，新規計測の場合は計測に移行し，既に計測データがある場合は目的のデータを特定する。なお，本被検者リスト画面に先立ってシステム立ち上げ時の時間待ちの表示画面を備えてもよく，更に本システムの目次的な役割をする表示画面を設けても良い。

図２４に示す被検者リスト画面について説明するに，その左側は被検者情報部で占められる。また，その右側全体の上部には被検者リストが，下部にはデ−タリストが表示されるようになっている。被検者情報部に表示される項目は図１０に関連して説明したのと同じである。被検者リストの項目は，ＩＤ（被検者ＩＤ番号），氏名，登録年月日（データ登録された日），計測回数（データ計測が行
われた回数），生年月日，年齢，身長，体重，コメント（被検者に関するコメン
ト）等を含む。被検者リストについては，これを縦スクロールバーでスクロールすることができ，被検者リストの項目については，これを水平（横）スクロールバーでスクロールできる。選択された被検者の行は強調表示される。選択された被検者に関するデータリストの項目は，ＩＤ，データの種類（生（Ｒａｗ）データかアベレージング（Ａｖｅｒａｇｉｎｇ）か），サンプリング間隔（データ計
測が行われたときの信号の，ミリ秒単位でのサンプリング間隔），サンプリング
時間（秒単位），分類（病気の分類情報），Ｄａｔｅ及びＴｉｍｅ（データ計測
が行われた日及び時刻），コメント（データに関するコメント）等を含む。デー
タリストについては，これを縦スクロールバーでスクロ−ルすることができ，データリストの項目については，これを水平（横）スクロールバーでスクロールできる。選択されたデータの行は強調表示される。

この被検者リスト画面によれば，被検者リストに各被検者の情報を１行表示する。これにより，上下に配列される各被検者の情報が明確に区分けできるので識別性を向上できるから，例えば誤って別の被検者を選択する誤操作を軽減できる
。この各被検者の情報は水平（横）スクロールバーでスクロールできるとともに
，選択された被検者の情報は縦長の被検者情報部に項目毎に上下に配列されるので視認性を損なうことがない。この場合，各被検者のデータリストの項目を前後（左右）に移動可能として視認性をより向上させても良い。更に，各被検者の情報を１行表示すことにより，一度に沢山の被検者を見ることができるので，縦スクロールバーでスクロールする回数を少なくできる。また，被検者リストの中から目的の被検者に関するデータを選択したい被検者リストにカーソルを合わせてクリックするだけの簡単な操作で下部のデータリストに表示できる。しかも，被検者リストとデータリストが上下に配置されているので，目線移動が少なくできるから，その関連性を認識しやすい。また，前記データリストは，そのエリアの上部にカーソルを移動してドラッグする簡単な操作でその大きさを自由に変えることができるので，データリストのリスト数に合わせて自由にその大きさを設定できる。

ステップＳ−５においては，被検者リスト画面上の被検者リストの中から所望の被検者の行が選択される。後述のアベレージング処理の場合は，必ずデータリスト中の生（Ｒａｗ）データの行が選択される。この後はフローはメニューにより４つに分岐される（Ｓ−６）。分岐の一つによれば，「ファイル（Ｆ）」のメ
ニューの「心磁システム終了（Ｘ）」のサブメニューが選択され，この場合はウ
インドウを閉じる等の終了処理が行われ（Ｓ−７），それによって，システムの
立ち下げが行われる（Ｓ−８）。その後，計算機１８の電源がＯＦＦにされ（Ｓ
−９），全てが終了する。分岐の残りによれば，アベレージング処理（Ｓ−１０
），データ解析（Ｓ−１１）及びデータ計測（Ｓ−１２）が行われる。アベレー
ジング処理は「データ解析（Ａ）」のメニューの「アベレージング（Ａ）」のサ
ブメニューを選択して実行可能である。また，データ解析は「データ解析（Ａ）
」のメニューの「単一波形表示（Ｗ）」，「重ね波形表示（Ｍ）」，「グリッド
マップ表示（Ｇ）」，「等磁線図（Ｂ）」，「時間積分図（Ｔ）」，「伝播時間
図（Ｐ）」及び「磁場源表示（Ｓ）」のサブメニューのいずれかを選択して実行
可能である。更に，データ計測は「データ計測（Ｑ）」のメニューの「計測パネ
ル（Ｐ）」のサブメニューを選択して実行可能である。ステップ１０，１１及び
１２の終了後はフローはステップＳ−４に戻る。ステップＳ−５の被検者選択，ステップＳ−１０のアベレージング処理，ステップＳ−１１のデータ解析及びステップＳ−１２のデータ解析については，その詳細はそれぞれ図１７〜２０に関連して以下に更に詳しく説明される。

図１７は図１６のステップＳ−５における被検者選択のフローを示す。被検者選択の場合は，フローはメニュー選択又は被検者選択により４つに分岐される。分岐の一つは被検者を指定して選択する（Ｓ−５−１）ことで被検者選択が終了する場合である。分岐のもう一つによれば，「リスト（Ｌ）」のメニューの「検
索（Ｓ）」のサブメニューが選択される。これによって，図１１に示す検索ダイ
アログボックスが開かれ（Ｓ−５−２），このダイアログボックスを用いて被検
者検索条件が入力される（Ｓ−５−３）。これによって被検者が検索され（Ｓ−
５−４），これにもとづいて図２４に示す被検者リスト画面上の被検者リストの
表示内容が変更される（Ｓ−５−５）。分岐の更にもう一つによれば，「リスト
（Ｌ）」のメニューの「解除（Ｘ）」のサブメニューが選択される。この場合は
，選択された被検者が全被検者リストに戻され（Ｓ−５−６），被検者リストの
表示内容の変更が行われる（Ｓ−５−５）。分岐の残りの一つによれば，「リス
ト（Ｌ）」のメニューの「登録（Ｒ）」のサブメニューが選択される。この場合
は，図１０に示す被検者登録ダイアログボックスが開かれ（Ｓ−５−７），被検
者情報が入力される（Ｓ−５−８）。これらのステップについては全ての被検者
の入力が終わるまでその入力終了の判断が行われ（Ｓ−５−９），その入力が終
了すると，被検者リストの更新が行われる（Ｓ−５−５）。

この実施例では，被検者の名前や住所等の文字入力を除いて，表示画面にプルダウンメニューを表示するなどして入力する複数の入力データ又は操作指示を表示し，その選択対象の中からマウスで特定の前記対象を指示して入力操作を行うようにしている。これにより，マウス操作でほとんどの操作が可能となるのでキーボードに不慣れな作業者に快適な操作環境を提供できるとともに，入力／操作の時間短縮を図ることができる。前記プルダウンメニューの複数の選択対象は，この装置或いはその入力／操作状態で入力／操作可能な選択対象が事前に設定され表示されるので誤入力／誤操作を軽減できる。また，この実施例では，入力エリアにカーソルを合わせてキーボードを介して入力することもできるので，操作者の入力の自由度を確保している。また，この実施例では，キーボードでの文字入力を想定しているが，入力時にキーボードのダイヤログを表示してこれをマウスで操作して入力しても良い。さらに手書入力ダイヤログを表示して，マウス操作で手書入力するようにしても良い。更には，前記デイスプレイ部にタッチパネルを備えて入力／操作を画面に指先又は入力ペンを介して操作しても良い。これらにより入力／操作の操作性を格段に向上できる。

図１８は図１６のステップＳ−１２におけるデ−タ計測のフローを示す。まず初めに，初期画面として心磁波形のグリッドマップが図２５に示すように表示される（Ｓ−１２−１）。同図において，操作領域部では，チャンネル選択，波形
モニタのＯＮ−ＯＦＦ，ＦＬＬ回路６のロック−アンロック，ＡＦＡ７のオフセット電圧の自動調整及びヒートフラッシュのための操作をそれぞれ行うことができ，更に，信号のサンプリング条件設定，波形表示のスケールの設定及びＡＦＡパラメータの設定が可能である。チャンネルは８×８の６４チャンネルからなり
，「全チャンネル選択」ボタンをクリックするか，又はチャンネルマトリックス
を対角線に沿って端から端までドラッグして，全チャンネルを選択できる。また
，行単位又は列単位でチャンネルマトリクス選択或いはその内の１つの項（以下「チャンネル項目」という。）をドラッグすれば，チャンネルを行単位又は列単
位で選択できる。いずれにしても，選択されたチャンネル項目に基づいてその項目に関連したチャンネルについて，その心磁波形が解析データ表示部に表示される。行単位又は列単位での選択の場合は図２６のように表示される。この場合，選択された波形は時間軸に関してはフルスケールいっぱいに拡大されて表示される。即ち，６４チャンネル全てが選択されている場合は，図２５に示すように解析データ表示部を上下左右（桝目）に分割して全てのチャンネル表示を優先させ
，行単位又は列単位での選択の場合は図２６のように，解析データ部を上下に分割して時間軸を左右にする馴染みのあるグラフの形態とすることで視認性を優先させた表示形態としている。

波形のモニタについては，「ＯＮ」ボタンが押されると，例えば，０．５ｓｅ
ｃから２ｓｅｃの間で指定された時間毎に信号の取り込みと波形の更新が繰り返され，被検者の心磁信号がモニタされる。また，「ＯＦＦ」ボタンが押されると
波形の更新が停止する。ＦＬＬについては，「Ｌｏｃｋ」ボタン又は「Ｕｎｌｏ
ｃｋ」ボタンをクリックすれば，６４個のＳＱＵＩＤセンサに対して磁場ロックを行ったり，そのロックを解除したりできる。その場合，一方のボタンが押されれば，他方が押されるまでそのままの状態が保たれる。これにより，選択されていない誤動作の状態を回避している。また，「ＡＦＡオフセット調整」ボタンを
クリックすれば，そのオフセット電圧が自動的に調整される。また，「ヒートフ
ラッシュ」ボタンをクリックすれば，図１４に示すヒートフラッシュ操作ダイアログボックスが開かれる。マウス又は矢印キーでチャンネルを選択し，「ＯＫ」
ボタンをクリックすれば，その選択されたチャンネルのＳＱＵＩＤについてヒートフラッシュ操作が実行される。「キャンセル」ボタンを押せば，ダイアログボ
ックスが閉じて処理は終了する。

サンプリングの時間（計測時間）及び間隔については，逆三角印のついた対応テキストボックスをクリックすれば，選択可能な数値のプルダウンメニューが開かれ，その中から所望の数字を選択するこできる。その選択可能な数字は，時間については，例えば，１ｓｅｃ，５ｓｅｃ，１０ｓｅｃ，３０ｓｅｃ，１ｍｉｎ及び２ｍｉｎであり，間隔については，例えば，０．１ｍｓｅｃ，０．５ｍｓｅｃ，１．０ｍｓｅｃ，２．０ｍｓｅｃ，４．０ｍｓｅｃ，５．０ｍｓｅｃ及び１０．０ｍｓｅｃである。時間は必要に応じて１ｓｅｃ程度から２４ｈ程度までの間で選ばれるようにしても良い。「スケール」ボックス内の「時間」はｍｓｅｃ
単位の時間スケール，即ち水平方向のスケールを，「信号」はＡ／Ｄ変換された
信号のスケール，即ち縦方向のスケールを意味する。これらについても，サンプリングの時間及び間隔の選択と同様に，対応するテキストボックスをクリックして開かれるプルダウンメニューの中から所望の数値が選択される。ＡＦＡパラメータは入力ゲインＩ_gain，出力ゲインＯ_gain，ローパスフイルタ（ＬＰＦ）の周
波数（基準周波数），ノッチフイルタ（ＢＥＦ）の周波数，ハイパスフイルタの
周波数（基準周波数）を含む。これらについても同様に，対応するテキスツボックスをクリックして開かれるプルダウンメニューの中から所望の数字又は文字が選択される。その数字又は文字は，I_gainについては，例えば，１，２，５，１
０，２０，５０，１００，２００，５００及び１０００であり，Ｏ_gainについて
は，例えば，１，１０及び１００であり，ＬＰＦについては，例えば，３０Ｈｚ
，５０Ｈｚ，８０Ｈｚ，１００Ｈｚ，２００Ｈｚ，４００Ｈｚ及び１ｋＨｚであり，ＢＥＦについてはＯｆｆ，５０Ｈｚ及び６０Ｈｚであり，そしてＨＰＦについては，例えば，０．０５Ｈｚ，０．１Ｈｚ及びＴｈｒｕ．である。なお，これらは選択の代わりにキーボードから入力されても良い。

チャンネルとしては６４チャンネル以外に更に例えば１６チャンネルの補助チャンネルが用意されていて，その補助チャンネルでは例えば心電波形が得られるようにされても良い。図２５の最下段には参照チャンネルとしての第１０チャンネルの波形が表示されているが，これはその補助チャンネル中の第１０チャンネルで得られた心電波形である。心磁波形は一般に磁気雑音を含み，一方，心電波形はそのような雑音を含まない。従って，表示されている心磁波形を参照チャンネルの心電波形と比較して心磁波形に磁気雑音が含まれているかどうかの情報が得られる。勿論，その心電波形は補助チャンネルではなく，正規の６４チャンネルの中の任意のチャンネルから得られるようにされていても良い。また，心電波形以外に脳波，血流波形，血圧波形等が用いられても良い。更に，妊婦の心電波形とその胎児の心磁波形が比較されるようにしても良い。また，参照波形としては，１チャンネルの参照波形だけでなく，複数チャンネルの参照波形が表示されるようにしても良い。更に，参照チャンネルは，生体からの信号だけでなく，保守等を目的とした種々の制御信号を入力するのに用いられても良い。

図１８のフローに戻るに，ステップＳ−１２−２では既述の要領でモニタチャンネルが選択され（Ｓ−１２−２），ＦＬＬのロックボタンが押されると，全て
のＳＱＵＩＤの磁場ロックがなされる（ステップＳ−１２−３）。その状態で，
計測パラメ−タであるサンプリングの時間と信号の設定並びにＡＦＡのパラメータの設定がなされる（ステップＳ−１２−４）。その設定については，この設定
を，その設定された条件を使うことにして次回から省略できる。これにより，毎回条件設定を行う必要がないので設定時間を短縮できる。なお，前記設定条件について名称を付ける等して記録し呼び出し可能としても良い。「計測」ボックス
の「開始」ボタンが押されると計測が開始され，図１５に示すように「計測中」の表示がなされると共に計測の進行状況を示すプログレスバーが表示される（ステップＳ−１２−５）。この実施例のプログレスバーは，棒グラフ形式で左から
右に処理進行の経過に従って棒が伸びるようにしているが，全体の処理内容（時間）を例えば１００として，現在の処理経過がわかれば良いので，円グラフ等でも良い。また，前記プログレスバーは，計測画面の所定の位置に周囲に前記計測画面を残して表示するようにする。これにより，表示画面の内容が大きく変わることがないので誤操作を軽減できる。

計測が開始されると，表示されている信号波形はそのまま固定化され，プログレスバーはその固定化された表示画面上に表示される。プログレスバーの更新は設定された時間が終了するまで，例えば毎秒繰り返される（Ｓ−１２−６）。「
計測」ボックスの「中止」ボタンが押されると，計測が中止する。計測が終了すると，図２６に示す画面表示がなされ，波形の確認がなされる（Ｓ−１２−７）
。その後，そのデータの保存の必要性が判断され（Ｓ−１２−８），保存が必要
な場合はメニュー「ファイル（Ｆ）」−「保存（Ｓ）」が選択されて，信号は保
存され，当該被検者のデータリストに追加される（Ｓ−１２−９）。その後，保
存が必要ない場合も含めて，計測がもう一度必要かどうかの判断がなされ（Ｓ−１２−１０），必要ならば以上のステップを繰り返し，必要なければデータ計測
の全ステップは終了する。その場合，画面表示は図２４の表示となるようにメニュー選択がなされる。なお，図２６はチャンネルとして第２列目の行のチャンネルが選択された例を示す。

図２６において，解析データ部の最下部にはスクロールボックス２６１が移動するスクロールバー２６２がある。スクロールボックス２６１はスクロールバー２６２の左右の両端間で移動可能なるもので，そのスクロ−ルボックスの幅ｗは時間スケールを表す。そのスクロールバー２６２の左右の両端間の時間幅は計測時間を表し，従って，表示されている波形は計測時間中に生じる波形の，スクロールボックス２６１の時間スケールｗに相当する一部分の拡大された波形である
。これにより，操作者は，現在解析データ部に表示されている波形が計測時間（
スクロールバー２６２の幅）の中でどのくらいの時間（スクロールボックス２６１の幅）かを示し，前記波形が計測時間の中で前半を示しているのか後半を示しているのか等を人目で把握できるので，視認性を向上させることができる。また
，前記スクロールボックス２６１の位置と解析データ部の波形を連動させているので，カーソルを前記スクロールボックス２６１に合わせてドラッグしながら移動させることで，解析データ部の表示領域を移動させて所定時間の波形を見るようにしても良い。このようにすれば，所定時間の波形を簡単に確認する等，前記スクロールバー２６２を目次的に扱いで計測内容を詳細に確認できる。

図１９は図１６のステップＳ−１０におけるアベレージング処理のフローを示す。アベレージング処理では，各チャンネルで測定された各データのノイズを取り除くために各チャンネル毎に加算されてその平均値が計算される。この各チャンネルのアベレージング処理にあたり基準となる時間を設定するために以下の処理がなされ，そして各チャンネルのアベレージング処理が実行される。まず，指定された被検者の生データが読み込まれる（Ｓ−１０−１）。これは，図２４の
リストデータ中の「データの種類」が生デ−タとなっている行を選択して行われる。これによって，初期表示として，図２７に示す第１列目のチャンネルの心磁波形が表示される（Ｓ−１０−２）。次いで，表示チャンネルが選択される（Ｓ
−１０−３）。図２７は第２列目のチャンネルが選択された例を示す。その後，
設定チャンネルが指定される（Ｓ−１０−４）。これは操作領域部にある「アベ
レージング条件」ボックス内の「チャンネル」のテキストボックスをクリックするか又は所望のチャンネルの波形を表示している領域をマウスでクリックして行われる。この場合，そのテキストボックス内の三角ボタンをクリックすればチャンネル数字が大きくなり，逆三角ボタンをクリックすればチャンネル数字が小さくなる。図２７は指定されたチャンネルが第２列第２行目のチャンネルである例を示す。このチャンネル指定によって，アベレージング処理の基準時間となるチャンネルが特定される。このチャンネルを特定するにあたっては，最も典型的な
，或いは分かり易い波形を備えたチャンネルを選定すると良い。仮に，選択した行又は列に良い波形を備えたチャンネルがなければ，再び（Ｓ−１０−３）のステップからやり直すことができる。

また，前記解析データ表示部の１つのチャンネルが特定されると，図２７に示すようにしきい値カ−ソル２７１が指定されたチャンネルの位置に表示される。これによって，指定されたチャンネルを視覚的に確認できる。図２７において，解析データ表示部８０５の上部に左右方向に移動可能な３個のスライダーカーソル２７３〜２７５が表示されているが，これらは図２７の表示画面が表示されると同時に自動的に表示される。ステップＳ−１０−５として，アベレージング条件であるしきい値，オフセット時間及びアベレージング時間が「アベレージング条件」ボックス内の対応するテキストボックスをクリックして設定される。三角ボタンをクリックすれば数字が大きくなり，逆三角ボタンをクリックすれば数字が小さくなる。しきい値の設定は「しきい値」テキストボックスでしきい値を表す数字を選択して行われる。これによって，その数字に対応する位置にしきい値カーソル２７１が自動的に移動する。その場合の移動位置はカーソル線によって明確に確認できる。「しきい値」テキストボックスのしきい値を表す数字の変更
（選択）としきい値カーソル２７１の移動は互いに連動しており，従ってしきい値はしきい値カーソル２７１を移動させることによっても設定可能である。スライダーカーソル２７３は波形の立上り部分がスライダーカーソル２７１によって設定されたしきい値と一致した時点（基準時点）位置を指し示すものであり，その指し示す位置はカーソル線によって明確に確認できる。このスライダーカーソル２７３はしきい値の変更に応じて変わる基準時点位置と一致するように該基準時点位置に追従して移動する。「オフセット」テキストボックスでオフセット時
間を表す数字を選択し，「時間」テキストボックスでアベレージング時間を表す
数字を選択すると，それに対応する位置にスライダーカ−ソル２７４及び２７５が移動する。その場合の移動位置はそれらのスライダーカーソルのカーソル線によって明確に確認できる。スライダーカーソル２７４及び２７５の移動は「オフセット」テキストボックス及び「時間」テキストボックスでの数字の選択と連動している。このため，オフセット時間及びアベレージング時間の設定はスライダーカーソルを移動させることによっても可能である。

ステップＳ−１０−６ではアベレージングモードとしてアベレージングを自動で行うべきか，マニュアルで行うべきかが設定される（Ｓ−１０−７）。その後
，「設定」ボタンをクリックすると（Ｓ−１０−８），加算処理が行われる。「
キャンセル」ボタンをクリックすれば，アベレージング条件は全てキャンセルされる。加算処理としては，各チャンネルにおいて，しきい値を越えた時刻ｔ（基準時点）が探索され（Ｓ−１０−９），次いで時刻ｔを中心に波形が表示され（
つまり，ｔ時点が画面の中心に位置付けられようにｔ−５０ｍｓｅｃからｔ＋５０ｍｓｅｃまでの波形が表示され）（Ｓ−１０−１０），そしてマニュアルモー
ドでキャンセルが選択されたかどうかの判断がなされ（Ｓ−１０−１１），キャ
ンセルが選択されたのでなければ波形の加算が実行される（Ｓ−１０−１２）。
ステップＳ−１０−９からステップＳ−１０−１２までのステップは各チャンネル毎に加算回数だけ繰り返され，その後加算されたデ−タは加算回数で割算され（Ｓ−１０−１３），アベレージング処理が終了する。このように，この図２７
に示すアベレージング処理の表示画面においては，操作領域部と解析データ表示部８０５の双方からアベレージング処理の諸条件の入力／操作を行うことができるので，例えば，大まかな条件を解析データ部で設定して，詳細な値を操作領域部で設定する等，操作者に多様な設定方法を提供すると共に条件設定の時間短縮を図ることができる。特に，この実施例では，解析データ表示部８０５の複数のチャンネルの中からカーソルを介して特定のチャンネルを指定してアベレージング処理の基準となるチャンネルとできから，誤操作が軽減され操作性が向上される。しかも，その指定されたチャンネル近傍にしきい値カーソルやスライダーカーソルを表示して視覚的な認識が可能である。更に，アベレージング処理の諸条件の入力／操作を前記チャンネル近傍に配置されるしきい値カーソルやスライダーカーソルで行うことができるので，操作者の目線移動を少なくして，波形表示に合わせて視覚的に設定できるから，誤操作を軽減して操作性を向上できる。なお，これらの設定条件，例えば最新の設定条件を記憶させて次回の設定条件として表示させたり，或いは各設定条件に名前を付して記憶させ呼び出すようにしても良い。

図２７に示したアベレージング処理の表示画面のアベレージング処理を行なう別の実施例を，図３５から図４２に基づいて詳細に説明する。図３５にアベレージング処理プログラムの表示画面を示す。図３５から図４２に示す表示画面は，図８に示した解析データ部８０５と操作領域部８０６とを詳細に示したものである。画面は主に５つのオブジェクトグループに分かれている。解析処理の手順に従って，オブジェクトグループを表すと，操作領域部としては，データ読み込みオブジェクト１０１，グラフ表示制御用オブジェクト１０２，アベレージング用オブジェクト１０３，全チャネル概要表示用オブジェクト１０４，ファイル保存及び終了オブジェクト１０５とが配置されている。表示画面上部には元波形（計測された磁場波形を元波形という）表示領域１０６が配置され，表示画面下部にはアベレージング表示領域１０７が配置されている。以下，操作手順に従って操作方法の実施例について説明を行う。

図３６を用いてファイル読み込みの実施例を説明する。ファイル読み込みはボタンＬＯＡＤ１０１ａを押すと，ファイル選択ダイアログ１０８が表示され，選択したいファイル名の位置にアイコンを移動させファイル名を選択して，ファイル選択の設定が行える。本実施例では，名前がｔｅｓｔ１．ｄａｔ名前のファイルが選ばれており，ハイライト１０８ａで分かりやすく表示されている。全てのファイルを選択したい時は，選択ボタン１０８ｂ中のＳｅｌｅｃｔ Ａｌｌボタン１０８ｂ−１を選択すると，全ファイルがハイライトされ全ファイルが選択可能となる。複数個のファイルを選択した場合，複数個の図１に示すウィンドウが表示される。以上のように１つ以上のファイルを選択後，ＯＫボタン１０８ｂ−２を押し次のモードに移動する。間違ってファイル選択ダイアログ１０８を開いてしまった場合は，Ｃａｎｃｅｌボタン１０８ｂ−３を選択して図３５に示した表示画面へ戻ることができる。選択したファイルの読み込むデータ数はデータ点数入力ウインドウ１０１ｂに数値を入力し，設定を行うことができる。選択したファイルのデータのサンプリング周波数をＦｒｅｑｕｅｎｃｙトグル１０１ｃで選択を行う。本実施例では，１ｋＨｚのサンプリング周波数が選択されている。Ｆｒｅｑｕｅｎｃｙトグル１０１ｃのように手動で選択する方式が面倒な場合は
，例えば選択したファイル内に測定時のサンプリング周波数の情報を測定データと共に記録しておき，自動でサンプリング周波数を決定する構成でも良い。

図３７を用いてグラフ表示制御について説明を行う。図３７の元波形表示領域１０６に表示されている波形１０６ａは，図３６中でハイライト１０８ａによって選択されたｔｅｓｔ１．ｄａｔのファイルでＴｒｉｇ．ｃｈボタン１０２ａで選択されたチャネル（本実施例では８０チャネル）の波形が表示されている。本実施例では，図２に示した６４個の磁気センサの他に１６チャネルの外部入力を有しており，８０チャネル目に心磁波形と同時に計測した第２誘導心電図の波形（図４３に後述する）を取り込んでおり，心電図の波形を用いてアベレージング処理を行った表示画面を示している。当然，Ｔｒｉｇ．ｃｈボタン１０２ａで選択されるチャネルは，８０チャネルの心電波形に限ることなく，図６に示した心磁波形の任意の１チャネルを選択することも可能である。Ｄｓｐ．ｃｈボタン１０２ｂはアベレージング表示領域１０７に表示されるチャネル（本実施例では８０チャネル）を選択できる。本実施例では，元波形表示領域１０６とアベレージング表示領域１０７とで選択されているチャネルが同じであるが，必ずしも同一にする必要はない。元波形表示領域１０６とアベレージング表示領域１０７のグラフ表示領域（時間領域幅）の変更は，テキストボックスｗｉｄｔｈ１０２ｃとテキストボックスｏｆｆｓｅｔ１０２ｄに数値を入力して行える。

テキストボックスｗｉｄｔｈ１０２ｃは，元波形表示領域１０６とアベレージング表示領域１０７の横軸１０６ｂ，１０７ｂに表示される時間幅（単位はｍｓ
）を設定し，テキストボックスｏｆｆｓｅｔ１０２ｄは，読み飛ばす時刻（単位はｍｓ）を設定できる。例えばデータの時刻０のデータ点から１０００ｍｓ読み飛ばしてデータ幅４０００ｍｓを読み出す場合には，テキストボックスｏｆｆｓｅｔ１０２ｄの数値を１０００ｍｓと入力し，テキストボックスｗｉｄｔｈ１０２ｃの数値を４０００ｍｓと入力する。本実施例ではテキストボックスｏｆｆｓｅｔ１０２ｄの数値を１ｍｓと入力し，テキストボックスｗｉｄｔｈ１０２ｃの数値を４０００ｍｓと入力した例を示す。本実施例の元波形表示領域１０６とアベレージング表示領域１０７の縦軸はＡ／Ｄコンバータによってデジタルデータに変換されたデータの値で表示されているが，磁場や電位の強度の値を用いて表示しても良い。さらに本実施例では縦軸の表示の幅は，表示区間に選ばれた時間幅内のデータの最大値と最小値が表示できる範囲に自動でスケーリングされる構成としている。

図３８を用いてピーク検出の設定方法について説明を行う。ピーク検出は，アベレージング用オブジェクト１０３によって行う。最初にテキストボックスＴｈｒｅｓｈｏｌｄ１０３ａ内に数字を入力してピーク検出の閾値（本実施例では１８０）を設定する。次に検出方法をトグルボタンＭｅｔｈｏｄ１０３ｂによって設定する。検出方法は閾値を越える値をもつ最大ピークを検出する方法（ｐｅａｋ＋，図３８記載），閾値より小さい値を持つ最小ピークを検出する方法（ｐｅ
ａｋ−），閾値より小さい値から大きい値に変化する時の閾値とのクロス点を検
出する方法（ｃｒｏｓｓ＋），閾値より大きい値から小さい値に変化する時の閾
値とのクロス点を検出する方法（ｃｒｏｓｓ−）の４種類が用意されている。本実施例では，最大ピークを検出する方法（ｐｅａｋ＋）が選択されている。さらにアベレージングする時間幅の選択はテキストボックスＦｒａｍｅＬｅｎｇｔｈ１０３ｃに数字を記入して行える。テキストボックスＤａｔａＬｅｎｇｔｈ１０３ｄに数字を入力して，検出ピーク前の時間幅の設定が行える。設定した内容（
１０３ａ，１０３ｂ，１０３ｃ，１０３ｄ）でピーク検出を行いたい場合は，ボタンＣＡＬＣ．ＰＥＡＫ１０３ｅを押し，計算を行わせる。ピーク検出が終わり次第，設定したＦｒａｍｅＬｅｎｇｔｈ１０３ｃとＤａｔａＬｅｎｇｔｈ１０３ｄの値は元波形表示領域１０６上に表示された元波形１０６ａ上にアベレージング表示範囲表示部Ｂ１，アベレージング表示範囲最大時刻Ｂ２，アベレージング表示範囲最小時刻Ｂ３によって設定範囲が元波形１０６ａ上に表示される。

また，表示画面上の波形中で抽出されたピーク点には全て白抜きの丸印Ｃ１，Ｃ２，Ｃ３，Ｃ４が表示される。丸印Ｃ１の上にはピーク検出された何番目のピークかが分かるように番号１が表示されている。本実施例では，ＦｒａｍｅＬｅｎｇｔｈ１０３ｃの値が１４００でＤａｔａＬｅｎｇｔｈ１０３ｄの値が７００で設定されている。元波形表示領域１０６上に表示された元波形１０６ａとアベレージング表示範囲表示部Ｂ１，アベレージング表示範囲最大時刻Ｂ２，アベレージング表示範囲最小時刻Ｂ３とを確認しながら，アベレージング処理したい場合にはボタンＮＥＸＴ１０３ｆを押すことによってアベレージング処理が行われ
，アベレージング結果はアベレージング表示領域１０７にアベレージング波形１０７ａとして表示される。図３８では，丸印Ｃ１の波形を観察後ボタンＮＥＸＴ１０３ｆを押してアベレージング処理をした例を示しており，アベレージング波形１０７ａはアベレージング回数１回のため，丸印Ｃ１の中心とした波形と同じ波形が表示されている。アベレージング処理は表示されている元波形表示領域１０６又はアベレージング表示領域１０７に表示されていないチャネルの全ての波形にも同時にアベレージング処理が施されている。図３８では，元波形表示領域１０６に１つの波形，アベレージング表示領域１０７に１つの波形しか表示していないが，元波形表示領域１０６，アベレージング表示領域１０７にそれぞれ複数個の波形を表示しても良い。

図３９は，ボタンＳＫＩＰ１０３ｇを押すことによって２個目のピークをアベレージングデータとして選択しなかった場合を示す。さらに図３９ではボタンＳＫＩＰ１０３ｇを押した後，図３８で４０００であった波形表示時間幅を設定するテキストボックスｗｉｄｔｈ１０２ｃの数字を８０００に変更した場合を示している。従って，アベレージング表示領域１０７に表示されるアベレージング波形１０７ａは図３８で選択された１個目の波形そのままである。ピーク検出処理を中断したくなった場合は，ボタンＣＬＥＡＲ１０３ｈを押すことにより，アベレージング表示領域１０７に表示されたアベレージング波形１０７ａは消滅し，アベレージング表示領域１０７には何も表示されない。アベレージング処理に選択されたピークは丸印Ｃ１が白抜きから黒丸に変更され，アベレージング処理に選択されなかったピークは２重丸印にマークが変更される。当然，丸印の形や色は白や黒，丸に限定されるものではない。以上のようにピーク点のアベレージング選択の状態を履歴として残すことによって，アベレージング処理終了後に再度確認が可能となる上に，本実施例では明記していないが，全処理終了後に不必要なピークのアベレージング波形を選択し除去することも可能である。

図３８又は図３９に示したアベレージング処理を続けていき，検出された６個目のピークで同様の処理を施す場合の実施例を図４０に示す。テキストボックスＴｈｒｅｓｈｏｌｄ１０３ａは１８０にトグルボタンＭｅｔｈｏｄ１０３ｂは最大ピークを検出する方法（ｐｅａｋ＋）に設定されているため，６個目より先のピークはピーク値が１８０を下回るため検出されていない。さらに次の時刻幅でピークを検出する場合にはテキストボックスｏｆｆｓｅｔ１０２ｄの値を例えば８０００と入力すれば次の８０００ｍｓから１６０００ｍｓまでを表示させて同様のアベレージング処理をすることが可能となる。また前述したが，丸印Ｃ１〜Ｃ６の中で黒丸になっているＣ１，Ｃ３，Ｃ４，Ｃ５，Ｃ６はアベレージング処理として選択されたピークのものを示しており，白抜きの２重丸印Ｃ２は，アベレージング処理に選択されていないものを示す。

図４１は，全チャネル概要表示用オブジェクト１０４中のボタンＶｉｅｗ Ａｖｅｒａｇｅ１０４ｂを押した場合の全チャネルのアベレージング処理結果の波形表示１０９を示す。本操作により全チャネルのアベレージング処理を行った波形１０９ａをセンサの位置に対応して表示できる。同様にボタンＶｉｅｗ Ｓｏｕｒｃｅ１０４ａを押すことにより，全チャネルの処理前の波形を表示できる。

図４２は，アベレージング処理後の波形のファイル保存処理を説明したものである。アベレージング処理後にファイルを保存する場合，ボタンＳＡＶＥ１０５ａを押して保存用ダイアログ１１０を表示させる。表示後，ファイル名を新規に作成したい場合には，テキストボックス１１０ａにファイル名を入力して，ボタンＯＫ１１０ｂを押すことにより任意のファイル名で保存が可能となる。既存のファイルに上書き保存したい場合は使用するファイル名をハイライトさせボタンＯＫ１１０ｂを押し，同一ファイルとしてアベレージング結果を保存できる。また，間違って保存用ダイアログ１１０を表示させた場合や，ファイルに保存しない場合には，ボタンＣａｎｃｅｌ１１０ｃを押すことによってファイル保存せずに図４１に示した表示画面に戻ることができる。全ての処理が終了した場合は，ボタンＱＵＩＴ１０５ｂを押すことによって，アベレージング処理のアプリケーションを終了させることができる。

以上がアベレージング処理を行う場合の手順及び画面表示方法である。上述のアベレージング処理をする場合以下のような前処理を施すと適切なアベレージング処理が可能となる。波形が呼吸性の雑音などの低周波成分の磁場雑音のために基線が揺らいでいる場合は，１つのチャネルの信号又は全てのチャネルの信号をアベレージング処理の前処理としてハイパスフィルタを通すと良い。１チャネルだけハイパスフィルタを通した場合，ハイパスフィルタを通した信号からピーク検出を行い，ピーク検出結果から前述したアベレージング処理を全チャネルの信号に施すことも可能である。ハイパスフィルタに用いるフィルタの種類は，例えばカットオフ周波数が０．１から３Ｈｚの範囲にあるハイパスフィルタをＩＩＲ（無限インパルスレスポンス）フィルタやＦＩＲ（有限インパルスレスポンス）フィルタといったデジタルフィルタを用いると良い。その他のノイズで商用電源ノイズ（例えば５０Ｈｚ，６０Ｈｚ）が雑音として問題であり，なお且つ低周波呼吸性の雑音が問題の場合は，くし形のノッチフィルタ（５０Ｈｚ又は６０Ｈｚの整数倍の周波数を除去するフィルタ）を用いることもできる。くし形のノッチフィルタとは，０Ｈｚ，５０Ｈｚ，１００Ｈｚ，１５０Ｈｚ…といった周波数を中心とした１Ｈｚ程のバンド幅を持つノッチフィルタのことをいう。以上のようなフィルタ処理をアベレージング処理の前処理として使用すると，より雑音の少ないアベレージング処理が可能となる。

次に上記のアベレージング処理を心電図で行う場合，心磁図と心電図とを同時に計測する構成を図４３によって説明する。心電計本体１２１は磁気雑音を避けるためシールドルームの外部に配置する。心電計本体１２１に接続された心電計の配線１２２−１はシールドルーム内まで挿入され，心電計の配線１２２−１の先には四肢誘導電線１２２−２が接続されている。四肢誘導電線１２２−２の先には磁気雑音を発生しにくいカーボン電極１２２−３が配置されている。本実施例では四肢誘導電線１２２−２やカーボン電極１２２−３で説明を行ったがこれらに限定されるものではなく，例えば胸部誘導電線（１２誘導）にカーボン電極を使っても良いし，非磁性材料を用いて電極部を構成して良い。さらに心電計本体１２１はＦＬＬ回路６や増幅器・フィルター・増幅器７などが収納されているラックの中に配置しも良い。心電計の波形は心磁の波形と同時に計算機８上にデジタルデータとして格納し，心磁と心電のデータ取り込みと同時にリアルタイムに心磁又は心電の任意の数の波形を表示可能である。リアルタイム波形の表示部は１つの計算機８の画面に限定されるものでなく，複数個の画面（コンピュータ用ディスプレイやテレビモニタなど）に選択的に同時にシールドルーム内外に表示できる。

図４４に心電計本体１２１の生体磁場計測装置との配置の関係を示した２つの実施例を示す。図４４は図１に示した生体磁場計測装置の実施例で示したシールルーム１の外部に配置したＦＬＬ回路６及び増幅器・フィルター・増幅器７と心電計本体１２１との位置関係を示している。図４４（ａ）の実施例では，ＦＬＬ回路６及び増幅器・フィルター・増幅器７を一つのユニットとしてボックス１２４に内蔵されており，電源スイッチ１２５のＯＮによってＦＬＬ回路６及び増幅器・フィルター・増幅器７部の電源が一括でＯＮされる。ボックス１２４の上部には，心電計本体１２１がおかれ，配線１２２―１に接続され，シールドルーム内に導入されている。配線１２２―１と，ＦＬＬ回路６とクライオスタット４との間に接続される配線とは，シールドルームの壁に設けられた同一の穴を通してシールドルーム内部に配線しても良いし，心電計のケーブルに誘導される電磁波などによる磁気のノイズが問題になる場合にはシールドルームの壁に設けられた異なる穴を通して配線１２２―１とＦＬＬ回路６とクライオスタット４との間に接続される配線とを別々にシールドルーム内部に配線しても良い。心電計本体１２１には心電図の波形を印刷するロール紙１２３が用意されている。ロール紙１２３に印刷されるものは，心電図の波形に限ることはなく，心磁図の波形を印刷してもよく，心磁図波形印刷専用のロール紙のユニットをボックス１２４に設置しても良い。図４４（ｂ）の実施例では，ＦＬＬ回路６と増幅器・フィルター・増幅器７と心電計本体１２１を一つのユニットとしてボックス１２４に内蔵されており，図２０は図１６のステップＳ−１１におけるデータ解析のフローを示す
。 データ解析はいろいろな種類の波形や線図を表示して診断に必要な情報を得ようとするもので，図９のメニューを選択していろいろな種類の波形や線図の画面を選択的に表示できる。即ち，「データ解析（Ａ）」の「単一波形表示（Ｗ）
」を選択すれば図２８に示さす単一波形画面が（Ｓ−１１−２），「デ−タ解析
（Ａ）」の「重ね合せ波形表示（Ｍ）」を選択すれば図２９す重ね波形画面が（
Ｓ−１１−３），「データ解析（Ａ）」の「グリッドマップ表示（Ｇ）」を選択
すれば図３０に示すグリッドマップ波形画面が（Ｓ−１１−４），「データ解析（Ａ）」の「等磁線図（Ｂ）」を選択すれば図３１に示す等磁線図画面が（Ｓ−
１１−５），「データ解析（Ａ）」の「伝播時間図（Ｐ）」を選択すれば図３２
に示す等伝播時間図画面が（Ｓ−１１−６），そして「デ−タ解析（Ａ）」の「
時間積分図（Ｔ）」を選択すれば図３３に示す等時間積分図画面が（Ｓ−１１−
７）それぞれ表示され，また，「ファイル（Ｆ）」の「心磁システムの終了（Ｘ
）」を選択すればシステムが終了する。

それぞれの画面において，操作領域部にあるラジオボタン（図面中の円形ボタン）をクリックすればそのクリックによって指定された波形又は線図の画面が代わって表示される。図２０において，分岐の部分を「メニューで分岐」とせずに「メニュー又はラジオボタンで分岐」としたのはそのためである。従って，この実施例によれば，前記図９のメニューを選択することなく，前記操作領域部にあるラジオボタンをクリックするだけで多様な解析データが得られるので，操作時間の短縮が図れるとともに，誤操作を軽減して操作性を向上できる。

図２８〜図３０において，「スケール」ボックス内の「磁束密度」は，ゼロレ
ベルを基準としたプラス側及びマイナス側のフルスケールの値（単位はピコテスラ（ｐＴ））であり，その値はそのテキストボックスの三角ボタンをクリックし
て開かれるプルダウンメニューで選択される。図２８〜３３において「表示成分
」ボックス内のラジオボックスをクリックして法線成分の波形又は接線成分の波形を選択して画面表示できる。

図２８ではチャンネルで選択された各チャンネルの波形が，解析データ部の左端をオフセット時間に合せて表示される。この解析データによれば，解析データ部に上下に配列表示される各チャンネルの波形の形状や大きさを比較できる。同様に図２９では，前記図２８で上下に配列された波形が重ねて表示され，その波形の形状や大きさを比較できる。また，図３０では，全てのチャンネルが前記図２８，図２９と同様にオフセット時間を基準にして表示される。従って，操作者は，必要によりチャンネルの数を選定して解析できる。

図３１において，解析データ部の右端には縦に細長の磁場強度指標ボックス３１０が配置されている。その磁場強度指標ボックスは互いに色の異なる１２個の区画に区切られている。これは，図３１に示す等磁線図画面の各島模様で示される磁場の強度範囲を色の種類で区別して視角的な（色彩的な）認識性を向上させるようにしたものである。即ち，その磁場強度指標ボックス３１０の長手方向の中心位置３１１は磁場強度がゼロの位置で，その中心位置よりも上方の区画を中心位置に近い順番に第１〜第６区画とそれぞれ呼ぶことにすれば例えば第１区画は０〜２ｐＴの磁場強度範囲に，第２区画は２〜４ｐＴの磁場強度範囲に，第３区画は４〜６ｐＴの磁場強度範囲に，第４区画は６〜８ｐＴの磁場強度範囲に，第５区画は８〜１０ｐＴの磁場強度範囲に，そして第６区画は１０〜１２ｐＴの磁場強度範囲にそれぞれ対応している。中心位置よりも下方の区画についても全く同じである。但し，中心位置よりも上方の区画はプラス方向の磁場強度を，下方の区画はマイナス方向の磁場強度を表している。図３１に示す等磁線図は，磁場強度指標ボックス３１０内の磁場強度範囲と色との対応関係の定めにしたがい
，磁場強度に応じて色分け表示される。なお，色として，磁場強度のプラス側を暖色系，マイナス側は寒色系とし，中心部を黄色とするようにしても良い。これにより，磁場の強弱を色彩的に認識できるので視認性を向上できる。しかも，この実施例によれば，解析データ部の近傍に磁場強度指標ボックス３１０を設けたので，比較対象の色，即ち，マップに付された色と磁場強度指標ボックス３１０の所定の色とを目線移動を大きく移動させることなく比べながら確認できるので
，前記磁場の強弱のレベルと色との関係を明確に判断できる。なお，この実施例では，磁場強度指標ボックス３１０を解析データ部の右端に設けているが，解析データ部の近傍であればよく，例えば，上部，下部，左側でも良い。

図３１において，「再構成パラメータ」ボックス内の「マップ数」は表示され
る等磁線図の数を，「最大値」は磁場強度指標ボックス３１０の両端部に相当す
る磁場強度を，「間隔」は磁場強度指標ボックス３１０内の各区画の長さに対応
する磁場範囲を意味する。その値については，これを対応するテキストボックスの三角又は逆三角ボタンをクリックして選択できる。解析データ表示部の最下段には参照チャンネルの心電波形と２個のマップ時刻選択用カ−ソル３１１及び３１２が表示されている。その２個のマップ時刻選択用カーソル３１１及び３１２間には間隔が同じ分割線が表示され，この線の数はマップ数選択によって選択されたマップの数と一致する。また，２個のマップ時刻選択用カーソル３１１及び３１２はその位置を独立に左右方向に移動可能で，その移動によってそのカーソル間の間隔が変わると，分割線の間隔も変わるが，分割線の間隔は常に等間隔である。勿論，各分割線にカーソルを備えて１本１本別々に設定するようにしても良い。図３１では，表示されている等磁線図の数は１６個であるが，これらの線図は心電波形上の分割線が位置する時点での線図であり，各マップについては，該マップがいつの時点のものであるかがわかるようにその時刻も表示される。これにより，図２６で説明したと同様に，操作者は，現在解析データ表示部に表示されているマップが解析時間（心電波形の幅）の中でどのくらいの範囲（２個のカーソル３１１と３１２の幅）を示し，前記マップが示す範囲が解析時間の中でどこの範囲なのか等を一目で把握できるので，視認性を向上させることができる
。また，前記マップが示す範囲を２個のカーソルをマウス操作で簡単に移動させることで設定できるので操作が容易である。更に，各分割線の間隔を自由に設定するようにすれば，疑問のある部分を密にして他の部分を疎にする等，操作者に多様な解析環境を提供できる。また，操作領域部にある「電流方向」のチェックボックスをクリックしてチェック印を表示すると，等磁線図上には矢印が表示される。この矢印が表示された等磁線図をアローマップと呼ぶ（図示なし）。矢印
については，その位置はチャンネルの位置（磁気センサの位置），長さは磁場の
強度，そして方向は磁場の方向を電流の方向に変換した場合のその電流の方向をそれぞれ示す。

図３２には，等伝播時間図が表示してあり，この図において，伝播時間の起点位置（図７のｔ1時点）の変更は参照波形上で移動するカーソル３２１の位置を
変えることによって可能であり，カーソル３２１の位置の変更はそのカーソルをマウスを用いてドラッグして可能である。

図３３には２つの時間積分図と１つの差分図が表示されている。等時間積分図の「差分表示」はチェックボックスをクリックしてチェック印を表示して簡単に表示できる。「差分表示」がチエックされると，参照波形上に４つのカーソル３
３１〜３３４が現れ，更に図示のように上方左右に２つの等時間積分図が，下方左側に差分図がそれぞれ表示される。２つの等時間積分図は心磁波形を，参照波形上でカーソル３３１及び３３２並びにカーソル３３３及び３３４を用いてそれぞれ設定された１００ｍｓｅｃ〜１４０ｍｓｅｃ及び１８０ｍｓｅｃ〜２４０ｍｓｅｃという時間範囲に亘って積分した値にもとづくもので，それぞれの時間範囲はカーソル３３１及び３３２並びにカーソル３３３及び３３４をマウスを用いてそれぞれドラッグして変え得る。差分図は２つの時間積分図の差を表すものである。チェック印がない場合は，カーソルについては２個のカーソル（例えばカーソル３３１及び３３２）だけが現れ，時間積分図については１つの等時間積分図だけが表示される。勿論，積分時間の変更はカーソルの位置を変えることによって可能である。このように，この実施例によれば，「差分表示」にチェック印
をクリックして，次の操作を促す２組のカーソルが表示されるので，操作の迷いを与えず操作時間の短縮が図れ，しかも，前記２組のカーソルをマウスで移動させることで時間範囲を簡単に設定できるから操作性を向上できる。

図２１はシステム調整のフローを示す。メニュー選択によってフローは５つに分岐される（Ｓ−１５）。この場合，現在表示されているデータはそのまま保存
され（Ｓ−１３），図３４に示すΦ−Ｖ特性曲線が表示されている（Ｓ−１４）
。「データ計測（Ｑ）」の「全自動調整（Ａ）」が選択された場合は，指定され
たチャンネルについてＩ_bias及びＶ_OFFという調整値が自動的に計算され（Ｓ−
１６）（その自動計算は後述），その計算された調整値はＦＬＬ回路に設定され
（Ｓ−１７），フローはステップＳ−１４に戻る。「データ計測（Ｑ）」の「Ｖ
_OFF調整（Ｖ）」が選択された場合は，指定されたチャンネルについてＶ_OFFが計
算され（Ｓ−１８）（Ｖ_OFFの計算は後述），その後，フローは既述のステップ
Ｓ−１７に進む。「データ計測（Ｑ）」の「マニュアル調整（Ｍ）」が選択され
た場合は，図１２に示すマニュアル調整ダイアログボックスが開かれる（Ｓ−１９）。オペレータはＩ_bias及びＶ_OFFという調整値をチャンネル毎に入力すると
，その入力された調整値は受け付けられ（Ｓ−２０），ダイアログボックス中の
「ＯＫ」ボタンが押されることによって調整値がＦＬＬ回路にセットされる（Ｓ−２１）。ダイアログボックスはこれによって閉じられ（Ｓ−２２），フローは
ステップＳ−１７に進む。「データ計測（Ｑ）」の「調整値ファイル（Ｆ）」が
選択されると，その後サブプルダウンメニューによりフローは更に２つに分岐される（Ｓ−２２）。即ち，サブプルダウンメニューの「開く（Ｏ）」を選択して
ファイル名をオペレータに問合わせ，オペレータは図１１の内容を含むファイル名を入力する（Ｓ−２３）。その調整値ファイルの内容はＦＬＬ回路に設定され
る（Ｓ−２４）。また，「データ計測（Ｑ）」−「調整値ファイル（Ｆ）」−「
上書き保存」又は「名前を付けて保存（Ａ）」を選択してステップＳ−２３と同
様のことが行われ（Ｓ−２５），その調整値を調整値ファイルに書き込むことが
できる（Ｓ−２６）。

図２２は図２１のステップＳ−１５における自動計算のフローを示す。

Ｓ−１５−１：全自動調整の場合は全てのＳＱＵＩＤチャンネルについて１チャンネルずつ以下の処理が行われる。処理中のチャンネルをｃｈとする。

Ｓ−１５−２：計算されたバイアス電流及びオフセット電圧はＩ_BIAS及びＶ_OF
Fという名称のメモリに保持されるものとする。Ｉ_BIAS及びＶ_OFFはチャンネルの
数だけ値を保持することができ，初期値は全て０とする。Φ−Ｖ特性曲線の，一時記憶される振幅をΔＶとする。

Ｓ−１５−３：各チャンネルｃｈに対して，バイアス電流Ｉ_bを０から図３４
の「走査パラメータ」ボックスで指定されたＩ_biasまでΔＩ_biasのステップで変
化（走査）させ，Φ−Ｖ特性曲線の振幅ΔＶが最も大きくなるＩ_biasをチャンネ
ルｃｈの最適なバイアス電流Ｉ_BIAS（ｃｈ）とする。

Ｓ−１５−４〜８：バイアス電流０からＩ_bまでの間で，Φ−Ｖ特性曲線の振
幅ΔＶは以下の処理で求められる。チャンネルｃｈのＳＱＵＩＤに与えられる外部磁場Φを０から図３４の「走査パラメータ」ボックスで指定されたΦｅｘｔまでΔΦのステップで変えて（走査して）行き，Ａ／Ｄ変換された信号を保存してその信号の最大値Ｍ_max及び最小値Ｖ_minを求める。

Ｓ−１５−９〜１０：ここで，最大値Ｍ_maxと最小値Ｖ_minの差がΔＶより大き
ければＩ_BIAS（ｃｈ）の値をＩ_biasで，ＶＯＦＦ（ｃｈ）の値を最大値Ｍ_maxと
最小値Ｖ_minの平均で，ΔＶを最大値Ｖ_maxと最小値Ｖ_minの差でそれぞれ置き換
える。もし，ΔＶの方が最大値Ｍ_maxと最小値Ｖ_minの差より大きければ前の値を
保持する。

Ｓ−１５−１１：以上の処理をバイアス電流Ｉ_bがＩ_biasになるまで繰り返し
たときのＩ_BIAS（ｃｈ），Ｖ_OFF（ｃｈ）がＳＱＵＩＤチャンネルｃｈの最適な
バイアス電流及びオフセット電圧となる。

Ｓ−１５−１２：また，以上の処理を全てのＳＱＵＩＤチャンネルに対して実行して全自動調整は終了する。

図２３は図２１のステップ１８におけるＶＯＦＦ調整フローを示す。

Ｓ−１８−１：ＶＯＦＦ調整の場合は全てのＳＱＵＩＤチャンネルについて１チャンネルずつ次の処理が行われる。

Ｓ−１８−２〜６：チャンネルｃｈのＳＱＵＩＤに与えられる外部磁場Φを０から図３４の「走査パラメータ」ボックスで指定されたΦ_extまでΔΦのステッ
プで変えて（走査して）行き，Ａ／Ｄ変換された信号の最大値Ｖ_maxと最小値Ｖ_m
inの差を求める。

Ｓ−１８−７：ＳＱＵＩＤチャンネルｃｈの最適なオフセット電圧ＶＯＦＦ（
ｃｈ）は最大値Ｖ_maxと最小値Ｖ_minの差として計算される。

Ｓ−１８−８：以上の処理を全てのＳＱＵＩＤチャンネルに対して実行してＶ_OF
F調整は終了する。以上の説明で「外部磁場Φの走査」は，計算機内部で設定さ
れた値をＤ／Ａ変換した信号をＳＱＵＩＤセンサに与えることによって実現できる。

本発明が実施される生体磁場計測装置の一実施例の概略構成図。 図１の生体磁場計測装置に用いられる磁気センサの配置構成を示す斜視図。 図１の生体磁場計測装置において用いられる，磁場の法線成分を検出する磁気センサ単体の斜視図。 図１の生体磁場計測装置において用いられる，磁場の法線成分を検出する磁気センサ単体の斜視図。 図１の生体磁場計測装置における磁気センサと被検者の胸部との位置関係を示す図。 図１の生体磁場計測装置において，各磁気センサにおいて計測された健常者の生体磁場（心磁）の各成分の時間波形図。 図１の生体磁場計測装置において，健常者について計測された特定の２チャンネルの心磁の接線成分の時間波形を示す図。 図１の生体磁場計測装置において，デイスプレイ部に表示される表示画面の基本的なレイアウトを示す図。 図１の生体磁場計測装置において，デイスプレイ部に表示される表示画面のメニューバー部における操作メニューを示す図。 図１の生体磁場計測装置において，デイスプレイ部に表示される表示画面中の操作メニューとして「リスト（Ｌ）」−「登録（Ｒ）」が選択された場合に開かれる被検者登録ダイアログボックスの内容を示す図。 図１の生体磁場計測装置において，デイスプレイ部に表示される表示画面中の操作メニューとして「リスト（Ｌ）」−「検索（Ｓ）」が選択された場合に開かれる検索ダイアログボックスの内容を示す図。 図１の生体磁場計測装置において，デイスプレイ部に表示される表示画面中の操作メニューとして「データ計測（Ｑ）」−「マニュアル調整（Ｍ）」が選択された場合に開かれるマニュアル調整ダイアログボックスの内容を示す図。 図１の生体磁場計測装置において，デイスプレイ部に表示される表示画面中の操作メニューとして「データ計測（Ｑ）」−「計測パネル（Ｐ）」が選択された場合に開かれる自動診断ダイアログボックスの内容を示す図。 図１の生体磁場計測装置において，デイスプレイ部に図２５又は図２６の表示画面が表示されている場合に「ヒートフラッシュ」ボタンが押されたときに開かれるヒートフラッシュ操作ダイアログボックスの内容を示す図。 図１の生体磁場計測装置において，計測中にデイスプレイ部に表示される計測プログレスバーを示す図。 図１の生体磁場計測装置において行われる全体の操作のフローを示す図。 図１６の操作フロー中の被検者選択ステップにおける被検者選択のフローを示す図。 図１６の操作フロー中のデータ計測ステップにおけるデータ計測のフローを示す図。 図１６の操作フロー中のアベレージング処理ステップにおけるアベレージング処理のフローを示す図。 図１６の操作フロー中のデータ解析ステップにおけるデータ解析のフローを示す図。 図１の生体磁場計測装置において行われるシステム調整のフローを示す図。 図２１のシステム調整フロー中の自動計算ステップにおける自動計算のフローを示す図。 図２１のシステム調整フロー中のＶＯＦＦ調整ステップにおけるＶＯＦＦ調整フローを示す図。 図１の生体磁場計測装置のデイスプレイ部に表示される被検者リスト画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示されるデータ計測画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に計測が終了したときに表示される波形確認画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示されるアベレージング画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示される単一波形表示画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示される重ね波形表示画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示されるグリッドマップ表示画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示される等磁線図表示画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示される伝播時間図表示画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示される時間積分図表示画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示されるシステム調整画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示されるアベレージング画面の初期画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示されるアベレージング画面のファイル選択画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示されるアベレージング画面の生波形表示画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示されるアベレージング画面のアベレージング処理設定画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示されるアベレージング画面のアベレージング処理中の画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示されるアベレージング画面のアベレージング処理中の画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示されるアベレージング画面の全体波形表示画面を示す図。 図１の生体磁場計測装置のデイスプレイ部に表示されるアベレージング画面のファイル保存処理画面を示す図。 図１の生体磁場計測装置の心電計の配置及び接続を示す図。 図１の生体磁場計測装置の心電計の配置を示す図。

符号の説明

１…磁気シールドルーム，２…被検者，３…ベッド，４…デユワ，５…自動補給装置，６…ＦＬＬ回路，７…増幅器・フイルター・増幅器，８…計算機，８−１…デイスプレイ部，８−２…キーボード，８−３…マウス，２０−１〜２０−８，２１−１〜２１−８，２２−１〜２２−８，２３−１〜２３−８，２４−１〜２４−８，２５−１〜２５−８，２６−１〜２６−８及び２７−１〜２７−８…磁気センサ，１０，１０’及び１０”並びに１１，１１’及び１１”…コイル
，１２，１２’及び１２”…ＳＱＵＩＤ，１３及び１４…センサ，３０…胸部，２６１…スクロールボックス，２６２…スクロールバー，２７１…しきい値カーソル，２７３〜２７５…スライダーカーソル，３１１，３１２及び３３１〜３３４…カーソル，８０１…タイトルバー部，８０２…メニューバー部，８０３…ツールバー部，８０４…被検者情報部，８０５…解析データ部，８０６…操作領域部，８０８…ステータスバー部，８０７−１…メッセイジバー部，８０７−２…日時表示部，８０８〜８１４…アイコン，１０１…データ読み込みオブジェクト
，１０１ａ…ボタンＬＯＡＤ，１０１ｂ…データ点数入力ウィンドウ，１０１ｃ…Ｆｒｅｑｕｅｎｃｙトグル，１０２…グラフ表示制御用オブジェクト，１０２ａ…Ｔｒｉｇ．ｃｈボタン，１０２ｂ…Ｄｓｐ．ｃｈボタン，１０２ｃ…テキストボックスｗｉｄｔｈ，１０２ｄ…テキストボックスｏｆｆｓｅｔ，１０３…アベレージング用オブジェクト，１０３ａ…テキストボックスＴｈｒｅｓｈｏｌｄ
，１０３ｂ…トグルボタンＭｅｔｈｏｄ，１０３ｃ…テキストボックスＦｒａｍｅＬｅｎｇｔｈ，１０３ｄ…テキストボックスＤａｔａＬｅｎｇｔｈ，１０３ｅ
；ボタンＣａｌｃ．Ｐｅａｋ，１０３ｆ…ボタンＮＥＸＴ，１０３ｇ…ボタンＳＫＩＰ，１０３ｈ…ボタンＣＬＥＡＲ，１０４…全チャンネル概要表示用オブジェクト，１０４ａ…ボタンＶｉｅｗ Ｓｏｕｒｃｅ，１０４ｂ…ボタンＶｉｅｗ
Ａｖｅｒａｇｅ，１０５…ファイル保存及び終了オブジェクト，１０５ａ…ボタンＳＡＶＥ，１０５ｂ…ボタンＱＵＩＴ，１０６…元波形表示領域，１０６ａ…元波形，１０６ｂ…横軸，１０７…アベレージング表示領域，１０７ａ…アベレージング波形，１０７ｂ…横軸，１０８…ファイル選択ダイアログ，１０８ａ…ハイライト，１０８ｂ…選択ボタン，１０８ｂ―１…Ｓｅｌｅｃｔ Ａｌｌボタン，１０８ｂ―２…ＯＫボタン，１０８ｂ―３…Ｃａｎｃｅｌボタン，１０９…波形表示，１０９ａ…波形，１１０…保存用ダイアログ，１１０ａ…テキストボックス，１１０ｂ…ボタンＯＫ，１１０ｃ…ボタンＣａｎｃｅｌ，１２０…ガントリー，１２１…心電計本体，１２２―１…心電計の配線，１２２―２…四肢誘導電線，１２２―３…カーボン電極，１２３…ロール紙，１２４…ボックス，１２５…電源スイッチ，Ｂ１…アベレージング範囲表示部，Ｂ２…アベレージング表示範囲最大時刻，Ｂ３…アベレージング表示範囲最小時刻，Ｃ１，Ｃ２，Ｃ３，Ｃ４，Ｃ５，Ｃ６…丸印。

Claims (1)

1. 生体の心臓から発する磁場を計測する複数の磁気センサと，計測を含む処理機能項目を示す処理機能表示部と，操作項目を表示する操作領域表示部と，及び，前記磁気センサの超伝導量子干渉素子に外部磁場φを０からΔφの間隔で変えて与えた時に，該外部磁場φに対応して発生する電圧Ｖと前記外部磁場φとで示されるφ−Ｖ特性曲線を表示する解析データ表示部とを具備する表示手段と，前記電圧Ｖの最大値と最小値からオフセット電圧を計算しオフセット電圧を調整する計算手段とを有することを特徴とする生体磁場計測装置。

Images