JP3595660B2 - Measuring device for biological magnetic field - Google Patents

Measuring device for biological magnetic field Download PDF

Info

Publication number
JP3595660B2
JP3595660B2 JP27015897A JP27015897A JP3595660B2 JP 3595660 B2 JP3595660 B2 JP 3595660B2 JP 27015897 A JP27015897 A JP 27015897A JP 27015897 A JP27015897 A JP 27015897A JP 3595660 B2 JP3595660 B2 JP 3595660B2
Authority
JP
Japan
Prior art keywords
displayed
magnetic field
time
analysis
subject
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.)
Expired - Fee Related
Application number
JP27015897A
Other languages
Japanese (ja)
Other versions
JPH11104095A (en
Inventor
和久 町田
宜史 河▲崎▼
博之 鈴木
啓二 塚田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP27015897A priority Critical patent/JP3595660B2/en
Priority to US09/164,271 priority patent/US6473518B1/en
Publication of JPH11104095A publication Critical patent/JPH11104095A/en
Priority to US10/282,158 priority patent/US6711281B2/en
Priority to US10/770,525 priority patent/US7194121B2/en
Application granted granted Critical
Publication of JP3595660B2 publication Critical patent/JP3595660B2/en
Priority to US11/717,667 priority patent/US7433506B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Description

【０００１】
【発明の属する技術分野】
本発明は、生体の脳の神経活動や心臓の心筋活動等、生体内の電流などが原因で発生する生体の磁場を計測するのに適した生体の磁場図表示装置に関する。
【０００２】
【従来の技術】
従来から磁気センサである超伝導量子干渉素子（ＳＱＵＩＤ）を用いて、生体から発生する微弱な磁場の分布を測定し、その測定結果から、生体内部の活動電流の位置を推定し、その分布をイメージングする多チャンネルの生体磁気イメージング装置が知られている。そのような従来例は、たとえば特開平４ー３１９３３４号あるいは特開平５ー１４６４１６号等の公開公報に開示されている。
【０００３】
【発明が解決しようとする課題】
前記従来例は、生体磁気イメージング装置に関する動作原理に関するものであり、その開示内容には実施化する上での技術課題や解決手段が開示されていない。また、前記従来例は脳の内部に発生する生体活動電流に関するものであって、他の部位に関する具体的な開示はなされていない。
【０００４】
本発明は、操作および視認性を向上させた生体の磁場の計測方法を提供することを目的とする。
【０００５】
【課題を解決するための手段】
上記目的を達成するために、本発明は、被検者の生体磁場を複数の位置で計測してその処理結果としての磁場図を表示する生体の磁場図表示装置であって、入力手段と表示部とを有する計算機を備え、計算機は、計測された生体磁場を複数のチャンネルの行または列として表示部の解析データ部に波形表示し、被検者の他の生体信号による参照チャンネルを解析データ部に表示された生体磁場の複数のチャンネルの行または列の下部に表示し、入浴手段からのデータ解析の指示を受け付けて、解析データ部に当該指示に該当する解析磁場を表示すると共にその下部に参照チャンネルを表示し、解析磁場の下部に表示される参照チャンネル上に、表示される解析磁場の時刻に対応する分割線および分割線の位置を移動するための時刻選択用カーソルを表示し、入力手段による前記カーソルの移動を受けて、移動後の選択時刻に対応する解析磁場を表示するものである。
【０００６】
上記の生体の磁場図表示装置は、前記被検者の他の生体信号は被検者の心電波形、脳波、血流波形または血圧波形であることを特徴とする。
上記の生体の磁場図表示装置の前記解析データ部に表示される解析磁場は、等磁線図、伝播時間図または時間積分図であることを特徴とする。
上記の生体の磁場図表示装置は、前記解析データ部に表示される解析磁場は等磁線図であり、前記時刻選択用カーソルは２箇所表示され、前記分割線は前記２つのカーソル間に等間隔に表示されることを特徴とする。
上記の生体の磁場図表示装置は、前記解析データ部に表示される解析磁場は等磁線図であり、前記時刻選択用カーソルは各分割線に設定され、前記分割線は個別に移動可能であることを特徴とする。
上記の生体の磁場図表示装置は、前記解析データ部に表示される解析磁場は伝播時間図であり、前記時刻選択用カーソルの移動により伝播時間の起点位置を変更可能であることを特徴とする。
【０００７】
上記の生体の磁場図表示装置は、前記解析データに表示される解析磁場は時間積分図であり、前記計算機は、前記入力手段からの差分表示の指示を受け付けて、前記参照チャンネル上に４つの時刻選択用カーソルを表示し、前記４つのカーソル位置に対応する時間範囲で時間積分した２つの時間積分図およびその差分図を前記解析データ部に表示することを特徴とする。
上記の生体の磁場図表示装置の前記参照チャンネルは、複数のチャンネルの中から所定のチャンネルを選択可能に表示されることを特徴とする。
【０００８】
【発明の実施の形態】
図１は本発明が実施される生体磁気計測装置の一実施例の概略構成を示す。環境磁気雑音の影響を除去するために、生体磁気計測装置は磁気シ−ルドル−ム１内に設置される。生体からなる被検体である被検者２は、ベッド３に仰向け状態で計測が行われる。被検者の生体面（胸部の場合は一般に胸壁に平行な面）はベッド３の面と略平行であるとし、そしてこの面は直交座標系（ｘ，ｙ，ｚ）のｘ−ｙ平面と平行であるものとする。被検者の胸部は曲面であると共に傾いているが、説明を簡単にするために略平行とする。
【０００９】
被検者２の胸部の上方には、冷媒である液体Ｈｅで満たされたデユワ４が配置され、該デユワは超伝導量子干渉素子（ＳＱＵＩＤ＝Ｓｕｐｅｒｃｏｎｄｕｃｔｉｎｇ Ｑｕａｎｔｕｍ Ｉｎｔｅｒｆｅｒｅｎｃｅ Ｄｅｖｉｃｅ）とそのＳＱＵＩＤに接続された検出コイルとを含む複数個の磁気センサを収容している。液体Ｈｅは磁気シ−ルド１の外部にある自動補給装置５から連続的に補給される。
【００１０】
磁気センサの出力は、被検者２から発生して検出コイルにより検出される生体磁場の強度（磁束密度と考えることもできる）と特定の関係をもつ電圧を出力し、その出力がＦＬＬ（Ｆｌｕｘ Ｌｏｃｋｅｄ ｌｏｏｐ）回路６に入力される。このＦＬＬ回路６は、ＳＱＵＩＤの出力を一定に保つように、ＳＱＵＩＤに入力された生体磁場（生体磁気）の変化を帰還コイルを介してキャンセルする（これを磁場ロックと呼ぶ）。その帰還コイルに流した電流を電圧に変換することにより、生体磁場信号の変化と特定の関係にある電圧出力を得ることができる。このように帰還コイルを介して検出する方式を取っているので、微弱の磁場を高感度に検出できる。
【００１１】
上記出力電圧は増幅器・フイルタ−・増幅器（ＡＦＡ）７に入力され、その出力はサンプリングされて、Ａ／Ｄ変換され、計算機８に取り込まれる。
【００１２】
計算機８はパ−ソナルコンピュ−タからなり、８−１はそのデイスプレイ部、８−２はキ−ボ−ド、そして８−３はマウスを示す。マウス８−３は画面上でカ−ソルを移動させて処理対象を選択するのに用いられる。この操作はキ−ボ−ドを操作することによっても行うことができる。ＡＦＡ７の入力ゲイン（Ｉｇａｉｎ）及び出力ゲイン（Ｏｇａｉｎ）は調整可能であり、また、ＡＦＡ７は第１の基準周波数以下の周波数信号を通過させるロ−パスフイルタ（ＬＰＦ）、第１の基準周波数よりも低い第２の基準周波数以上の周波数信号を通過させるハイパスフイルタ（ＨＰＦ）及び商用電源周波数をカットするノッチフィルタ（ＢＥＦ）を含む。計算機８は各種の処理を行うことができ、その処理結果はデイスプレイ部８−１に表示され得る。なお、図１で示す前記計算機８は一実施例を示したものであり、これに限定されるものではない。例えば、タッチパネルを備えたデイスプレイを備えたものや、マウスに変えて他の座標指示装置、例えばトラックボールやジョスティック等を使用したものでもよい。また場合によっては公衆電話回線を介して接続される計算機でもよい。
【００１３】
ＳＱＵＩＤとしては例えば１例として直流ＳＱＵＩＤが用いられる。ＳＱＵＩＤに外部磁場が与えられたときに、それに対応する電圧（Ｖ）が発生するようにＳＱＵＩＤには直流バイアス電流（Ｉｂｉａｓ）が流される。その外部磁場を磁束Φで表すと、ＶのΦに対する特性曲線すなわちΦ−Ｖ特性曲線は周期関数で与えられる。計測に当っては、それに先立って、ＦＬＬ回路６のオフセット電圧（ＶＯＦＦ）を調整してΦ−Ｖ特性曲線の直流電圧をゼロレベルにする操作が行われる。更に、ＡＦＡ７の入力がゼロのときその出力がゼロとなるようにＡＦＡ７のオフセット電圧（ＡＯＦＦ）の調整が行われる。
【００１４】
ＳＱＵＩＤに外部から大きな磁場がかかると、その磁場がＳＱＵＩＤによってトラップされ、その正規の動作がなされなくなる。その場合は、ＳＱＵＩＤを加熱して一旦常伝導状態にし、そしてその後その加熱を止めて、トラップされた磁場を取り除くことができる。その場合の、ＳＱＵＩＤの加熱操作をヒ−トフラッシュと呼ぶ。
【００１５】
図２は磁気センサの配置構成を示す。磁気センサの検出コイルには生体磁場の接線成分（生体面すなわちｘ−ｙ平面に略平行な成分）を検出するコイルと生体磁場の法線成分（生体面すなわちｘ−ｙ平面に直交する成分）を検出するコイルがある。生体磁場の接線成分を検出するコイルとしては、コイル面がｘ方向及びｙ方向をそれぞれ向いた２つのコイルが用いられ、また、生体磁場の法線成分を検出するコイルとしてはコイル面がｚ方向を向いたコイルが用いられる。複数個の磁気センサ２０−１〜２０−８、２１−１〜２１−８、２２−１〜２２−８、２３−１〜２３−８、２４−１〜２４−８、２５−１〜２５−８、２６−１〜２６−８及び２７−１〜２７−９は、図２に示されるように、生体面すなわちｘ−ｙ平面と略平行な面上にマトリックス状に配置される。磁気センサの数は任意であってよいが、図２では、磁気センサのマトリックスは８行８列からなっているから、磁気センサの数は８×８＝６４である。各磁気センサは、図２に示されるように、その長手方向が生体面すなわちｘ−ｙ平面に対して垂直な方向（ｚ方向）と一致するように配置される。なお、この一実施例ではベッド面とセンサーのＸ−Ｙ面とを平行にしているが、測定精度を高めるには体に接近させる方が良く、傾けるようにすることができる。但し、被検者である人体は常に動いているので、人体に密着させるとこの動きが検出部を動かし、かえって高精度の検出が困難となる。
【００１６】
図３は磁気センサの各々の、生体磁場の法線成分Ｂｚを検出するセンサの構成を示す。同図において、超伝導線（Ｎｉ−Ｔｉ線）で作られたコイルはそのコイル面がｚ方向を向くように配置される。このコイルは互いに逆向きの２つのコイル１０及び１１の組み合わせからなり、被検者２に近い方のコイル１０は検出コイルとされ、遠い方のコイル１１は外部磁場雑音を検出する参照コイルとされる。外部磁場雑音は被検者よりも遠い信号源から生じており、したがって、その雑音信号は検出コイル１０及び参照コイル１１の両方によって検出される。一方、被検者からの磁場信号は微弱であり、したがって、その生体磁場信号は検出コイル１０によって検出されるが、参照コイル１１はその生体磁場信号にほとんど感応しない。このため、検出コイル１０は生体磁場信号と外部磁場雑音信号を検出し、参照コイル１１は外部磁場雑音信号を検出するから、両コイルで検出された信号の差をとることによりＳ／Ｎ比の高い生体磁場の計測が可能となる。これらのコイルはＳＱＵＩＤ１２を実装した実装基板の超伝導線を介してＳＱＵＩＤの入力コイルに接続され、これによって、検出された生体磁場信号の法線方向の成分ＢｚがＳＱＵＩＤに伝達される。
【００１７】
図４は磁気センサの各々の、生体磁場の接線成分Ｂｘ及びＢｙを検出するセンサの構成を示す。同図において、接線方向の生体磁場成分検出用のセンサでは平面コイルが用いられる。すなわち、検出コイル１０’及び１０”並びに参照コイル１１’及び１１”は平面コイルからなり、これらは互いにｚ方向において間隔づけられている第１及び第２の平面にそれぞれ配置される。これらのコイルは法線成分用と同様にＳＱＵＩＤ１２’及び１２”の実装基板の入力コイルに接続される。４角柱の互いに直交する２面に、これらのＢｘ成分検出用のセンサ１３及びＢｙ成分検出用のセンサ１４が貼付けられ、これによってＢｘ成分及びＢｙ成分を検出し得るセンサが形成される。
【００１８】
接線成分Ｂｘ、Ｂｙについては、これを図４に示される磁気センサを用いて検出する以外に、図３の磁気センサで得られた法線成分Ｂｚをｘ、ｙについて偏微分して求めてもよい。この場合は一つの磁気センサで接線成分Ｂｘ、Ｂｙと法線成分Ｂｚとの両方を検出し、測定することができる。
【００１９】
図５は磁気センサと被検者２の被計測部である胸部３０との位置関係を示す。示されている点は図２に示されるマトリックス上の行と列との交点すなわち被検者２の計測点すなわち計測位置を表す。これらの各計測位置をチャンネルとも呼ぶ。図からわかるように、この実施例では、被検者２の身長方向をｙ方向とし、被検者２の横方向をｘ方向としている。
【００２０】
図６は図５に示されるそれぞれの計測位置の生体磁場の測定結果を示す。この測定結果は各計測位置に対応する磁気センサでそれぞれ検出された信号を基に上述の処理を行うことで得られた、時間的に変化する生体磁場波形をマトリックス中の対応するチャンネル毎に示すものである。この実施例では心臓の筋肉が発する磁場を検出できる位置に各チャンネルを設けたので、図６の波形は心磁波形を示している。なお、心臓の筋肉から発せられる磁場を計測することで得られた波形を心磁波形と呼ぶ。図６に示されるように、チャンネル毎に計測された計測デ−タをチャンネル毎にその位置に対応させて表示する場合、これをグリッドマップ表示と呼ぶ。図６はある健常者についての心磁の計測結果の波形を示す。ここで（ａ）は接線成分Ｂｘの心磁波形を、（ｂ）は接線成分Ｂｙの心磁波形を、そして（ｃ）は法線成分Ｂｚの心磁波形をそれぞれ示す。
【００２１】
図７はある健常者について特定の２チャンネルに絞って計測された接線成分Ｂｘの心磁波形を示す。実線はあるチャンネルの、点線は他のチャンネルの心磁波形を示す。心臓の心室が脱分極した時間帯Ｔ１つまり収縮期のＱＲＳ波でのそれぞれの波形ピ−クの時間がｔＱ、ｔＲ及びｔＳとしてそれぞれ示されている。また、心臓の再分極過程（拡張期）であるＴ波の時間帯はＴ２として示されている。計測されたデ−タを表示する場合、グリッドマップ表示の他に、行毎に又は列毎に単一チャンネルのデ−タを表示するようにしてもよいし、全チャンネル又は複数のチャンネルのデ−タを重ねて表示するようにしてもよい。前者を単一波形表示、後者を重ね波形表示と呼ぶ。
【００２２】
得られる心磁波形デ−タについては、その平均化（アベレ−ジング）処理を行ったり、等磁線図や時間積分図等を作成したりして、その結果を表示することができる。たとえば、図７を参照するに、ＱＲＳ波の立上り部分がしきい値ＳＬと一致した時点から予め定められた時間（ｔＯＦＦ）だけ遡り、その遡った時点ｔ１から予め定められた時間Ｔ３だけ経過した時点ｔ２までの間のデ−タを予め定められた回数だけ加算することが行われる。これがアベレ−ジングで、予め定められる時間Ｔ３をアベレ−ジング時間、ｔＯＦＦをオフセット時間と呼ぶ。心磁波形デ−タは予め定められた時間範囲に亘って積分されてもよい。その時間積分値が等しい点（チャンネル）を結んで作られたマップを時間積分図と呼ぶ。また、心磁波形信号値が等しい点を結んで作られたマップを等磁線図と呼ぶ。なお、各チャンネルは粗く設定されているので、予め等磁線の間隔つまり磁場強度差を設定して各チャンネル間を直線補間して等磁線を描くことにより、より診断に適した図を作ることができる。図７に示される心磁波形について言えば、時点ｔ１からＱＲＳ波のピ−ク位置時点（ｔＲ時点）までの時間を伝播時間と呼び、その伝播時間が等しい点を結んで作られたマップを伝播時間図と呼ぶ。しきい値ＳＬの設定レベルは変更可能である。ｔ１時点については、これを、ＱＲＳ波の立上り部分がしきい値ＳＬと一致した時点を基準として決定しているが、ＱＲＳ波の立ち下がり部分がしきい値ＳＬと一致した時点を基準にして決定してもよい。ｔ１時点は更に、ＱＲＳ波のピ−ク位置時点（ｔＲ時点）を検出し、この時点を基準として決定されてもよい。計測された生体磁場信号は生体内の電気的生理現象によって生成されるもので、その発生源は電流双極子モデルによって近似される。その磁場発生源の電流双極子は等磁線図上に合成表示され、それを磁場源表示と呼ぶ。
【００２３】
被検者の登録から、その登録された被検者のデ−タ計測を行って、その計測されたデ−タの解析を行うまでの一連の操作はデイスプレイ８−１に表示される表示画面を見ながら行われる。このため、その一連の操作の説明に先立ってまずその表示画面のレイアウトを説明する。
【００２４】
図８は、図１のデイスプレイ８−１に表示される表示画面の基本的なレイアウトを示す。表示画面の上方部は上から順番に配置されたタイトルバ−部８０１、メニュ−バ−部８０２及びアイコンが配置されているツ−ルバ−部８０３によって占められる。上記各部は表示領域やエリアと考えることもできる。これらの配置は他の処理目的、たとえば、被検者の登録や読み出し、磁場の計測、計測デ−タの解析のための処理、などにおいてもその表示画面で共通して表示される。これにより使用し易さが増し、計測や処理の時間が短縮できる。
【００２５】
表示画面の中央部は左から右へ順番に配置された被検者情報部８０４、線図や波形のような解析デ−タが表示される解析デ−タ部８０５及び操作領域部８０６によって占められている。また、下部はステ−タスバ−８０７で占められ、これは、左側に配置された、次の操作に関するガイドメッセ−ジを表示するメッセ−ジバ−部８０７−１とその右側に配置された日時表示部８０７−２からなる。なお、前記メッセ−ジバ−部８０７−１と前記日時表示８０７−２を１つの表示領域としてもよい。
【００２６】
この実施例における表示画面では、常に最上部にこのシシテムの名称が表示されるタイトルバ−部８０１と、このシシテムの基本的な操作を行うメニュ−バ−部８０２、及び前記メニューバー部８０２における使用頻度の高い操作が可能なツ−ルバ−部８０３が配置されているので、使用者は表示画面が変わるたびに操作エリアを探す必要がなく、常に表示画面の上部を見れば現在動作中のシステムを知ることができる。しかも、表示画面の上部は、人が文章を読む場合を想定すると明らかなように、先ず最初に目を向ける部分であることから、このシステムの操作における基本的な事項を最上部に設けることで、自然な形で使いやすさを向上している。また、表示画面の中央部は、その中央にはこの表示画面の主体をなす解析デ−タ部８０５を大きく設けることで見易さを向上するとともに、その右側にこの表示画面に特有の操作領域部８０６を設けることで、右手操作における表示画面と操作部の配置と同様な配置としているので、違和感なく操作を行うことができる。したがって、この表示画面をタッチパネル付きの表示画面に採用しても、操作領域部８０６を操作する右手が解析デ−タ部８０５を邪魔することがない。また、同様に、解析デ−タ部８０５の左側に確認機能しかない被検者情報部８０４を設けているので、常に患者を確認しながら操作を行うことができる。しかも、この左側の位置は、右手操作における最も遠い位置となるのでタッチパネル付きの表示画面に採用しても表示の見易さに影響をきたすことがない。
【００２７】
また、解析デ−タ部８０５及び操作領域部８０６は被検者リスト及びその被検者のデ−タリストが表示されるときだけはそれらによって置き換えられる（図２４）。被検者情報部８０４には、被検者リスト画面（図２４）が表示されているときはその画面中の被検者リストにおいてカ−ソルが置かれる被検者の情報が常時表示され、また、解析デ−タ部に線図や波形のような解析デ−タが表示されているときは（図２５〜３４）、その表示されている解析デ−タが得られた被検者の情報が常時表示される。これによって、表示されている解析デ−タとその解析デ−タが得られた被検者との関係を明確に知ることができる。このように、このシステムの表示画面においては、メニュ−バ−部８０２と同様に、常に、被検者情報部８０４が表示画面の定位置（左側）に表示されているので、使用者は表示画面が変わるたびに被検者情報エリアを探す必要がなく、常に表示画面の所定位置（左側）を見れば知ることができる。
【００２８】
タイトルバ−部にはフレ−ムの名称、具体的には、「Ｍｕｌｔｉｃｈａｎｎｅｌ ＭＣＧ Ｓｙｓｔｅｍ」という名称が表示される（図２４〜３４）。操作領域部にはボタンやテキストボックスのような操作要素が配置されている。メニュ−バ−部は操作メニュ−を選択する部分で、メニュ−は「ファイル（Ｆ）」、「編集（Ｅ）」、「リスト（Ｌ）」、「デ−タ計測（Ｑ）」、「デ−タ解析（Ａ）」及び「ヘルプ（Ｈ）」からなり、操作の順序にしたがって配置されている。
【００２９】
図９は表示画面中のメニュ−部におけるそれぞれの操作メニュ−の内容を示し、これらのメニュ−の内容はそれぞれ対応するメニュ−ボタンをクリックすることによってプルダウンメニュ−として表示される。このため、操作メニューを必要としないときは、前記各メニューを呼び出すためのキーワードのみメニュ−バ−部にコンパクトに表示しているので、解析データ部や操作領域部等の各作業で必要な表示エリアを広く設定することができる。そして、操作メニューが必要なときは、操作手順にしたがって配列された前記キーワードをメニューバー部から選択することによって表示して操作指示することができる。この際、前記キーワードは、文字の配列（左から右）に準じて配列されているので、自然な形で操作指示することができる。
【００３０】
「ファイル（Ｆ）」のプルダウンメニュ−は、ペ−ジレイアウトダイアログボックス（図示せず）を開いて、ペイジレイアウトをセットする「ペ−ジレイアウト（Ｕ）」、印刷前のプレビュ−のための「プレビュ−（Ｖ）」、デ−タの印刷のための「印刷（Ｐ）」及びＭｕｌｔｉｃｈａｎｎｅｌ ＭＣＧ Ｓｙｓｔｅｍを終了させる「心磁システムの終了（Ｘ）」という項目を含む。
【００３１】
「編集（Ｅ）」のプルダウンメニュ−は「削除（Ｄ）」という項目を含む。被検者リスト画面（図２４）上の被検者リスト中の被検者の計測デ−タは同画面上のデ−タリストに表示されるが、「削除（Ｄ）」という項目はそのデ−タリスト中のカ−ソルが置かれたデ−タを削除するためのもので、予めデ−タが選択されていなければならない。このメニュ−がクリックされると、削除してよいかどうかの確認を行う確認ダイアログボックスが開かれる。削除が必要な場合は、そのダイアログボックス中の「ＯＫ」というボタンがクリックされ、削除をキャンセルしたいときは「キャンセル」というボタンがクリックされる。この確認ダイアログボックスの表示及び操作により、誤って被検者の計測デ−タを削除する誤操作を軽減することができる。
【００３２】
「リスト（Ｌ）」のプルダウンメニュ−は「登録（Ｒ）」、「リスト（Ｌ）」、「削除（Ｄ）」、「検索（Ｓ）」及び「解除（Ｘ）」という項目を含む。プルダウンメニュ−の「登録（Ｒ）」というボタンがクリックされると、図１０に示される被検者登録ダイアログボックスが開かれる。これは被検者に関するデ−タを登録するときに用いられるもので、登録できる項目は登録年月日、所定桁までの被検者のＩＤ番号、氏名、生年月日、身長、体重、性別、病気の分類情報を示す分類及び被検者コメントである。そのダイアログボックス中の「登録」ボタンがクリックされると、登録のために入力されたデ−タは登録されると共に、その入力ボックスはすべてクリアされて再入力が可能となり、「キャンセル」というボタンがクリックされれば、そのすべての入力ボックスはクリアされ、そして「終了」というボタンがクリックされれば、被検者登録ダイアログボックスが閉じられる。
【００３３】
プルダウンメニュ−の「リスト（Ｌ）」がクリックされると、被検者リスト画面（図２４）が表示される。また、プルダウンメニュ−の「削除（Ｄ）」は被検者リスト画面（図２４）上の被検者リスト中のカ−ソルが置かれた被検者を削除するためのもので、それがクリックされると、削除前に、プルダウンメニュ−の「削除（Ｄ）」がクリックされるときと同様に削除してよいかどうかの確認を行う確認ダイアログボックスが開かれ、削除が必要な場合は、そのダイアログボックス中の「ＯＫ」というボタンがクリックされ、削除をキャンセルしたいときは「キャンセル」というボタンがクリックされる。被検者が削除されると、その被検者に関するデ−タリスト中のすべてのデ−タも削除される。
【００３４】
プルダウンメニュ−の「検索（Ｓ）」がクリックされると、図１１に示される検索ダイアログボックスが表示される。被検者及びデ−タを検索して、被検者リスト画面に検索された被検者及びデ−タのみを表示する。被検者及びデ−タの検索対象は全被検者及び全デ−タとする。被検者名、登録日付、性別、デ−タの種別、計測日付、診断結果、コメント、検査技師、体位などをキ−ワ−ドとして与えてデ−タが検索される。前記のキ−ワ−ドを複数個組み合わせた複合検索を行うことができる。プルダウンメニュ−の「解除（Ｘ）」は検索された被検者及びデ−タだけの表示を解除して全被検者及び全デ−タを表示するのに用いられる。「デ−タ計測（Ｑ）」のプルダウンメニュ−は「調整値ファイル（Ｆ）」、「全自動調整（Ａ）」、「ＶＯＦＦ調整（Ｖ）」、「マニュアル調整（Ｍ）」、「調整値ファイル（Ｆ）」、「計測パネル（Ｐ）」、「自動波形診断（Ｄ）」及び「ＡＦＡオフセット調整（Ｏ）」という項目を含む。
【００３５】
「調整ファイル（Ｆ）」がクリックされると、「開く（Ｏ）」、「上書き保存（Ｓ）」及び「名前をつけて保存（Ａ）」を内容とするサブプルダウンメニュ−が表示される。サブプルダウンメニュ−中の「開く（Ｏ）」はシステム調整画面（図３４）を表示して、指定された調整値ファイルを開いてシステムの調整値すなわちＩｂｉａｓ及びＶＯＦＦをシステムに設定するのに用いられる。「上書き保存（Ｓ）」は確認ダイアログボックスを開いて、現在の調整値を現在開かれている調整値ファイルに上書き保存するのに用いられる。「名前をつけて保存（Ａ）」は現在の調整値デ−タを名称を変えて別の調整値ファイルに保存するのに用いられる。
【００３６】
プルダウンメニュ−の「全自動調整（Ａ）」が選択されると、システム調整画面（図３４）が表示され、後述の図２２のフロ−にしたがってバイアス電流Ｉｂｉａｓ及びＦＬＬ回路６のオフセット電圧ＶＯＦＦが自動的に調整される。「ＶＯＦＦ調整（Ｖ）」が選択されると、システム調整画面（図３４）が表示され、後述の図２３のフロ−にしたがってＦＬＬ回路６のオフセット電圧ＶＯＦＦが自動的に調整される。「マニュアル調整（Ｍ）」が選択されると、図１２に示されるマニュアル調整ダイアログボックスが開かれる。オペレ−タはスクロ−ルバ−とマウスを使ってチャンネルを選択し、バイアス電流Ｉｂｉａｓ及びオフセット電圧ＶＯＦＦを変更する。入力された値に問題なければ、「ＯＫ」ボタンをクリックすることによりその値の設定が行われる。「キャンセル」ボタンがクリックされると、変更は無効となり、ダイアログボックスが閉じられる。
【００３７】
「計測パネル（Ｐ）」がクリックされると、グリッドマップを含むデ−タ計測画面（図２５）が表示される。「自動波形診断（Ｄ）」がクリックされると、図１３に示される自動診断ダイアログボックスが開かれる。自動波形診断は、図３４の表示画面中のΦ−Ｖ特性曲線の振幅、中央値及び周期を自動的に調べ、Φ−Ｖ特性曲線を表示する際あるいは表示中の同曲線の表示を最新状態に更新する際に、指定された範囲外にあるチャンネルを計測に不適当な状態としてオペレ−タに通知するものである。計測に不適当な状態のチャンネルが検出された場合、エラ−ダイアログボックスを開いてエラ−メッセ−ジを表示することによってオペレ−タに通知するが、当該チャンネルのΦ−Ｖ特性曲線を別色で表示して不適当な状態であることを示してもよい。振幅最小値がチェックされている場合、Φ−Ｖ特性曲線における最大値と最小値の差が振幅最小値よりも小さいときが計測に不適当な状態であるとする。中央値が指定されている場合、前記最大値と最小値との平均値の絶対値が指定された中央値よりも大きい場合が計測に不適当な状態であるとする。周期が指定されている場合、Φ−Ｖ特性曲線におけるΦの周期が第１のテキストボックス（下限値）と第２のテキストボックス（上限値）で指定された範囲外にあるときが計測に不適当な状態であるとする。振幅最小値、中央値及び周期の自動診断を有効にするには、各項目の左側のチェックボックスをマウスでクリックして×印が表示されている状態にして、対応するテキストボックスに所望する値を入力すればよい。このようにして入力された自動診断のパラメ−タは「ＯＫ」ボタンを押すことによって有効となり、ダイアログボックスが閉じられる。「キャンセル」ボタンが押されると、入力されたパラメ−タは無効となり、ダイアログボックスは閉じられる。「ＡＦＡオフセット調整（Ｏ）」はＡＦＡ７のオフセット電圧ＡＯＦＦを調整するときに用いられるもので、「ＡＦＡオフセット調整（Ｏ）」がクリックされると、単一波形表示に属するデ−タ計測画面（図２５）が表示される。
【００３８】
「デ−タ解析（Ａ）」のプルダウンメニュ−は「アベレ−ジング（Ａ）」、「単一波形表示（Ｗ）」、「重ね波形表示（Ｍ）」、「グリッドマップ表示（Ｇ）」、「等磁線図（Ｂ）」、「時間積分図（Ｔ）」、「伝播時間図（Ｐ）」、「磁場源表示（Ｓ）」、「ラインモ−ド（Ｌ）」及び「塗潰しモ−ド（Ｆ）」を含む。
【００３９】
「アベレ−ジング（Ａ）」がクリックされるとアベレ−ジング画面（図２７）が、「単一波形表示（Ｗ）」がクリックされると単一波形表示字画面（図２８）が、「重ね波形表示（Ｍ）」がクリックされると重ね波形表示画面（図２９）が、「グリッドマップ表示（Ｇ）」がクリックされるとグリッドマップ表示画面（図３０）がそれぞれ表示され、メニュ−の左側に選択されたことを示すチェックマ−ク（×）が表示される。また、「等磁線図（Ｂ）」がクリックされると等磁線図画面（図３１）が、「時間積分図（Ｔ）」がクリックされると時間積分図画面（図３２）が、「伝播時間図（Ｐ）」がクリックされると伝播時間図画面（図３３）がそれぞれ表示され、メニュ−の左側に選択されたことを示すチェックマ−ク（×）が表示される。また、「磁場源表示（Ｓ）」がクリックされるとメニュ−の左側に逆三角マ−クが表示され、等磁線図を表示する際に電流ダイポ−ルで近似した磁場源が重ねて表示される（図なし）。等磁線図、時間積分図、伝播時間図及び磁場源表示の画面は「ラインモ−ド（Ｌ）」が選択されると、線間が塗潰しなしで表示され、「塗潰しモ−ド（Ｆ）」が選択されると、線間が塗潰された状態で表示される。現在のモ−ドを示すために、メニュ−の左側に選択を表すチェックマ−クが表示される。
【００４０】
「ヘルプ（Ｈ）」のプルダウンメユ−は「目次（Ｃ）」、「キ−ワ−ドで検索（Ｓ）」及び「バ−ジョン情報（Ａ）」を含み、それぞれヘルプウインドウを開いて目次を示すこと、キ−ワ−ドでトピックを検索すること、バ−ジョンダイアログボックスを開くことのために用いられる。
【００４１】
ツ−ルバ−部８０３には「被検者登録」（８０８）、「被検者リスト」（８０９）、「印刷」（８１０）、「プレビュ−」（８１１）、「システム調整」（８１２）、「デ−タ計測」（８１３）、「デ−タ解析」（８１４）というアイコンが配置される。これらは、その図示は省略されているが、メニュ−の機能とリンクしていて、プルダウンメニュ−の項目のうちの使用頻度の高いものを選択することができる。すなわち、「被検者登録」（８０８）は「リスト（Ｌ）」の「登録（Ｒ）」と、「被検者リスト」（８０９）は「リスト（Ｌ）」の「リスト（Ｌ）」と、「印刷」（８１０）は「ファイル（Ｆ）」の「印刷（Ｐ）」と、「プレビュ−」（８１１）は「ファイル（Ｆ）」の「プレビュ−（Ｖ）」と、「システム調整」（８１２）は「デ−タ計測（Ｑ）」の「マニュアル調整（Ｍ）」と、「デ−タ計測」（８１３）は「デ−タ計測（Ｑ）」の「計測パネル（Ｐ）」と、そして「デ−タ解析」（８１４）は「デ−タ解析（Ａ）」の「グリッドマップ表示（Ｇ）」とそれぞれ対応している。このようにアイコンとリンクしているメニュ−については、そのメニュ−をアイコンをクリックするだけで選択することもできる。したがって、使用頻度の高い操作は、解析データ部に隣接したアイコンをクリックするだけで簡単に操作することができるので、前記メニュ−バ−部の操作に比べて短時間にしかも認識し易いアイコンで操作することができる。なお、前記アイコンは利用者によって選択するようにしてもよく、また、利用頻度（回数）に準じて自動的に前記ツ−ルバ−部８０３に表示するようにしてもよい。
【００４２】
次に、システムの調整操作を含めて、被検者の登録から、その登録された被検者のデ−タ計測を行って、その計測されたデ−タの解析を行うまでの一連の操作を図１４〜３４を参照しながら説明する。
【００４３】
図１６は全体の操作のフロ−を示すもので、計算機８の電源がＯＮにされると（Ｓ−１）、オペレ−テイングシステムが立ち上げられ、プログラム起動アイコンがデイスプレイ部８−１に表示される（Ｓ−２）。そのアイコンの中からＭｕｌｔｉｃｈａｎｎｅｌ ＭＣＧ Ｓｙｓｔｅｍのプログラムのアイコンが選択されると（Ｓ−３）、図２４に示される被検者リスト画面が代わって表示される（Ｓ−４）。
【００４４】
この実施例に係るシステムにおいては、システム立ち上がりの初期画面として図２４に示す被検者リスト画面を表示する。この理由は、被検者と、該被検者の計測または解析データの関係が極めて重要であるため、このシステムでは被検者情報をキワードとしてデータ管理していることに起因する。すなわち、計測データや解析データは被検者情報がないと管理ができないためである。このため、このシステムでは、被検者リスト画面において、先ず被検者を登録または登録されているときは被検者を特定し、次に、新規計測の場合は計測に移行し、既に計測データがある場合は目的のデータを特定する。なお、本被検者リスト画面に先立ってシステム立ち上げ時の時間待ちの表示画面を備えてもよく、更に本システムの目次的な役割をする表示画面を設けてもよい。
【００４５】
図２４に示される被検者リスト画面について説明するに、その左側は被検者情報部で占められる。また、その右側全体の上部には被検者リストが、下部にはデ−タリストが表示されるようになっている。被検者情報部に表示される項目は図１０に関連して説明したのと同じである。被検者リストの項目は、ＩＤ（被検者ＩＤ番号）、氏名、登録年月日（デ−タ登録された日）、計測回数（デ−タ計測が行われた回数）、生年月日、年齢、身長、体重、コメント（被検者に関するコメント）等を含む。被検者リストについては、これを縦スクロ−ルバ−でスクロ−ルすることができ、被検者リストの項目については、これを水平（横）スクロ−ルバ−でスクロ−ルすることができる。選択された被検者の行は強調表示される。
【００４６】
選択された被検者に関するデ−タリストの項目は、ＩＤ、デ−タの種類（生（Ｒａｗ）デ−タかアベレ−ジング（Ａｖｅｒａｇｉｎｇ）か）、サンプリング間隔（デ−タ計測が行われたときの信号の、ミリ秒単位でのサンプリング間隔）、サンプリング時間（秒単位）、分類（病気の分類情報）、Ｄａｔｅ及びＴｉｍｅ（デ−タ計測が行われた日及び時刻）、コメント（デ−タに関するコメント）等を含む。デ−タリストについては、これを縦スクロ−ルバ−でスクロ−ルすることができ、デ−タリストの項目については、これを水平（横）スクロ−ルバ−でスクロ−ルすることができる。選択されたデ−タの行は強調表示される。
【００４７】
この被検者リスト画面によれば、被検者リストに各被検者の情報を１行表示する。これにより、上下に配列される各被検者の情報が明確に区分けすることができるので識別性を向上することができるから、たとえば誤って別の被検者を選択する誤操作を軽減することができる。この各被検者の情報は水平（横）スクロ−ルバ−でスクロ−ルすることができるとともに、選択された被検者の情報は縦長の被検者情報部に項目毎に上下に配列されるので視認性を損なうことがない。この場合、各被検者のデ−タリストの項目を前後（左右）に移動可能とすることにより視認性をより向上させてもよい。更に、各被検者の情報を１行表示すことにより、一度に沢山の被検者を見ることができるので、縦スクロ−ルバ−でスクロ−ルする回数を少なくすることができる。また、被検者リストの中から目的の被検者に関するデ−タを選択したい被検者リストにカーソルを合わせてクリックするだけの簡単な操作で下部のデ−タリストに表示することができる。しかも、被検者リストとデ−タリストが上下に配置されているので、目線移動が少なくできるから、その関連性を認識しやすい。また、前記デ−タリストは、そのエリアの上部にカーソルを移動してドラッグする簡単な操作でその大きさを自由に変えることができるので、デ−タリストのリスト数に合わせて自由にその大きさを設定することができる。
【００４８】
ステップＳ−５においては、被検者リスト画面上の被検者リストの中から所望の被検者の行が選択される。後述のアベレ−ジング処理の場合は、必ずデ−タリスト中の生（Ｒａｗ）デ−タの行が選択される。この後はフロ−はメニュ−により４つに分岐される（Ｓ−６）。分岐の一つによれば、「ファイル（Ｆ）」というメニュ−の「心磁システム終了（Ｘ）」というサブメニュ−が選択され、この場合はウインドウを閉じる等の終了処理が行われ（Ｓ−７）、それによって、システムの立ち下げが行われる（Ｓ−８）。その後、計算機１８の電源がＯＦＦにされ（Ｓ−９）、すべてが終了する。
【００４９】
分岐の残りによれば、アベレ−ジング処理（Ｓ−１０）、デ−タ解析（Ｓ−１１）及びデ−タ計測（Ｓ−１２）が行われる。アベレ−ジング処理は「デ−タ解析（Ａ）」というメニュ−の「アベレ−ジング（Ａ）」というサブメニュ−を選択することによって実行可能である。また、デ−タ解析は「デ−タ解析（Ａ）」というメニュ−の「単一波形表示（Ｗ）」、「重ね波形表示（Ｍ）」、「グリッドマップ表示（Ｇ）」、「等磁線図（Ｂ）」、「時間積分図（Ｔ）」、「伝播時間図（Ｐ）」及び「磁場源表示（Ｓ）」というサブメニュ−のいずれかを選択することにより実行可能である。更に、デ−タ計測は「デ−タ計測（Ｑ）」というメニュ−の「計測パネル（Ｐ）」というサブメニュ−を選択することにより実行可能である。ステップ１０、１１及び１２の終了後はフロ−はステップＳ−４に戻る。ステップＳ−５の被検者選択、ステップＳ−１０のアベレ−ジング処理、ステップＳ−１１のデ−タ解析及びステップＳ−１２のデ−タ解析については、その詳細はそれぞれ図１７〜２０に関連して以下に更に詳しく説明される。
【００５０】
図１７は図１６のステップＳ−５における被検者選択のフロ−を示す。被検者選択の場合は、フロ−はメニュ−選択又は被検者選択により４つに分岐される。分岐の一つは被検者を指定して選択する（Ｓ−５−１）ことで被検者選択が終了する場合である。分岐のもう一つによれば、「リスト（Ｌ）」というメニュ−の「検索（Ｓ）」というサブメニュ−が選択される。これによって、図１１に示される検索ダイアログボックスが開かれ（Ｓ−５−２）、このダイアログボックスを用いて被検者検索条件が入力される（Ｓ−５−３）。これによって被検者が検索され（Ｓ−５−４）、これにもとづいて図２４に示される被検者リスト画面上の被検者リストの表示内容が変更される（Ｓ−５−５）。分岐の更にもう一つによれば、「リスト（Ｌ）」というメニュ−の「解除（Ｘ）」というサブメニュ−が選択される。この場合は、選択された被検者が全被検者リストに戻され（Ｓ−５−６）、被検者リストの表示内容の変更が行われる（Ｓ−５−５）。分岐の残りの一つによれば、「リスト（Ｌ）」というメニュ−の「登録（Ｒ）」というサブメニュ−が選択される。この場合は、図１０に示される被検者登録ダイアログボックスが開かれ（Ｓ−５−７）、被検者情報が入力される（Ｓ−５−８）。これらのステップについてはすべての被検者の入力が終わるまでその入力終了の判断が行われ（Ｓ−５−９）、その入力が終了すると、被検者リストの更新が行われる（Ｓ−５−５）。
【００５１】
この実施例では、被検者の名前や住所等の文字入力を除いて、表示画面にプルダウンメニューを表示するなどして入力する複数の入力データまたは操作指示を表示し、その選択対象の中からマウスで特定の前記対象を指示することで入力操作を行うようにしている。これにより、マウス操作でほとんどの操作が可能となるのでキーボードに不慣れな作業者に快適な操作環境を提供できるとともに、入力／操作の時間短縮を図ることができる。前記プルダウンメニューの複数の選択対象は、この装置あるいはその入力／操作状態で入力／操作可能な選択対象が事前に設定され表示されるので誤入力／誤操作を軽減できる。また、この実施例では、入力エリアにカーソルを合わせてキーボードを介して入力することもできるので、操作者の入力の自由度を確保している。また、この実施例では、キーボードでの文字入力を想定しているが、入力時にキーボードのダイヤログを表示してこれをマウスで操作して入力してもよい。さらに手書入力ダイヤログを表示して、マウス操作で手書入力するようにしてもよい。更には、前記デイスプレイ部にタッチパネルを備えて入力／操作を画面に指先または入力ペンを介して操作してもよい。これらにより入力／操作の操作性を格段に向上することができる。
【００５２】
図１８は図１６のステップＳ−１２におけるデ−タ計測のフロ−を示す。まず初めに、初期画面として心磁波形のグリッドマップが図２５に示されるように表示される（Ｓ−１２−１）。同図において、操作領域部では、チャンネル選択、波形モニタのＯＮ−ＯＦＦ、ＦＬＬ回路６のロック−アンロック、ＡＦＡ７のオフセット電圧の自動調整及びヒ−トフラッシュのための操作をそれぞれ行うことができ、更に、信号のサンプリング条件設定、波形表示のスケ−ルの設定及びＡＦＡパラメ−タの設定が可能である。
【００５３】
チャンネルは８×８の６４チャンネルからなり、「全チャンネル選択」ボタンをクリックするか、又はチャンネルマトリックスを対角線に沿って端から端までドラッグすることによって、全チャンネルを選択することができる。また、行単位又は列単位でチャンネルマトリクスをドラッグすれば、チャンネルを行単位又は列単位で選択することができる。いずれにしても、選択されたチャンネルについては、その心磁波形が解析デ−タ部に表示される。行単位又は列単位での選択の場合は図２６のように表示される。この場合、選択された波形は時間軸に関してはフルスケ−ルいっぱいに拡大されて表示される。すなわち、６４チャンネルすべてが選択されている場合は、図２５に示すように解析データ部を上下左右（桝目）に分割してすべてのチャンネル表示を優先させ、行単位又は列単位での選択の場合は図２６のように、解析データ部を上下に分割して時間軸を左右にする馴染みのあるグラフの形態とすることで視認性を優先させた表示形態としている。
【００５４】
波形のモニタについては、「ＯＮ」ボタンが押されるとたとえば０．５ｓｅｃから２ｓｅｃの間で指定された時間毎に信号の取り込みと波形の更新が繰り返され、被検者の心磁信号がモニタされる。また、「ＯＦＦ」ボタンが押されると波形の更新が停止する。ＦＬＬについては、「Ｌｏｃｋ」ボタンまたは「Ｕｎｌｏｃｋ」ボタンをクリックすれば、６４個のＳＱＵＩＤセンサに対して磁場ロックを行ったり、そのロックを解除したりすることができる。その場合、一方のボタンが押されれば、他方が押されるまでそのままの状態が保たれる。これにより、選択されていない誤動作の状態を回避している。
【００５５】
また、「ＡＦＡオフセット調整」ボタンをクリックすれば、そのオフセット電圧が自動的に調整される。また、「ヒ−トフラッシュ」ボタンをクリックすれば、図１４に示されるヒ−トフラッシュ操作ダイアログボックスが開かれる。マウス又は矢印キ−でチャンネルを選択し、「ＯＫ」ボタンをクリックすれば、その選択されたチャンネルのＳＱＵＩＤについてヒ−トフラッシュ操作が実行される。「キャンセル」ボタンを押せば、ダイアログボックスが閉じて処理は終了する。
【００５６】
サンプリングの時間（計測時間）及び間隔については、逆三角印のついた対応テキストボックスをクリックすれば、選択可能な数値のプルダウンメニュ−が開かれ、その中から所望の数字を選択するこできる。その選択可能な数字は、時間についてはたとえば１ｓｅｃ、５ｓｅｃ、１０ｓｅｃ、３０ｓｅｃ、１ｍｉｎ及び２ｍｉｎであり、間隔についてはたとえば０．１ｍｓｅｃ、０．５ｍｓｅｃ、１．０ｍｓｅｃ、２．０ｍｓｅｃ、４．０ｍｓｅｃ、５．０ｍｓｅｃ及び１０．０ｍｓｅｃである。時間は必要に応じて１ｓｅｃ程度から２４ｈ程度までの間で選ばれるようにしてもよい。「スケ−ル」ボックス内の「時間」とあるのはｍｓｅｃ単位の時間スケ−ルすなわち水平方向のスケ−ルを、「信号」とあるのはＡ／Ｄ変換された信号のスケ−ルすなわち縦方向のスケ−ルを意味する。これらについても、サンプリングの時間及び間隔の選択と同様に、対応するテキストボックスをクリックすることによって開かれるプルダウンメニュ−の中から所望の数値が選択される。
【００５７】
ＡＦＡパラメ−タは入力ゲインＩｇａｉｎ、出力ゲインＯｇａｉｎ、ロ−パスフイルタ（ＬＰＦ）の周波数（基準周波数）、ノッチフイルタ（ＢＥＦ）の周波数、ハイパスフイルタの周波数（基準周波数）を含む。これらについても同様に、対応するテキスツボックスをクリックすることによって開かれるプルダウンメニュ−の中から所望の数字又は文字が選択される。その数字又は文字は、Ｉｇａｉｎについてはたとえば１、２、５、１０、２０、５０、１００、２００、５００及び１０００であり、Ｏｇａｉｎについてはたとえば１、１０及び１００であり、ＬＰＦについてはたとえば３０Ｈｚ、５０Ｈｚ、８０Ｈｚ、１００Ｈｚ、２００Ｈｚ、４００Ｈｚ及び１ｋＨｚであり、ＢＥＦについてはＯｆｆ、５０Ｈｚ及び６０Ｈｚであり、そしてＨＰＦについてはたとえば０．０５Ｈｚ、０．１Ｈｚ及びＴｈｒｕ．である。なお、これらは選択の代わりにキ−ボ−ドから入力されてもよい。
【００５８】
チャンネルとしては６４チャンネル以外に更にたとえば１６チャンネルの補助チャンネルが用意されていて、その補助チャンネルではたとえば心電波形が得られるようにされてもよい。図２５の最下段には参照チャンネルとしての第１０チャンネルの波形が表示されているが、これはその補助チャンネル中の第１０チャンネルで得られた心電波形である。心磁波形は一般に磁気雑音を含み、一方、心電波形はそのような雑音を含まない。したがって、表示されている心磁波形を参照チャンネルの心電波形と比較することにより心磁波形に磁気雑音が含まれているかどうかの情報が得られる。もちろん、その心電波形は補助チャンネルではなく、正規の６４チャンネルの中の任意のチャンネルから得られるようにされていてもよい。また、心電波形以外に脳波、血流波形、血圧波形等が用いられてもよい。更に、妊婦の心電波形とその胎児の心磁波形が比較されるようにしてもよい。また、参照波形としては、１チャンネルの参照波形だけでなく、複数チャンネルの参照波形が表示されるようにしてもよい。更に、参照チャンネルは、生体からの信号だけでなく、保守等を目的とした種々の制御信号を入力するのに用いられてもよい。
【００５９】
図１８のフロ−に戻るに、ステップＳ−１２−２では既述の要領でモニタチャンネルが選択され（Ｓ−１２−２）、ＦＬＬのロックボタンが押されると、すべてのＳＱＵＩＤの磁場ロックがなされる（ステップＳ−１２−３）。その状態で、計測パラメ−タであるサンプリングの時間と信号の設定並びにＡＦＡのパラメ−タの設定がなされる（ステップＳ−１２−４）。その設定については、この設定を、その設定された条件を使うことにして次回から省略することができる。これにより、毎回条件設定を行う必要がないので設定時間を短縮することができる。なお、前記設定条件について名称を付ける等して記録し呼び出し可能としてもよい。
【００６０】
「計測」ボックスの「開始」ボタンが押されると計測が開始され、図１５に示されるように「計測中」の表示がなされると共に計測の進行状況を示すプログレスバ−が表示される（ステップＳ−１２−５）。この実施例のプログレスバ−は、棒グラフ形式で左から右に処理進行の経過にしたがって棒が伸びるようにしているが、全体の処理内容（時間）を例えば１００として、現在の処理経過がわかればよいので、円グラフ等でもよい。また、前記プログレスバ−は、計測画面の所定の位置に周囲に前記計測画面を残して表示するようにする。これにより、表示画面の内容が大きく変わることがないので誤操作を軽減することができる。
【００６１】
計測が開始されると、表示されている信号波形はそのまま固定化され、プログレスバ−はその固定化された表示画面上に表示される。プログレスバ−の更新は設定された時間が終了するまで、たとえば毎秒繰り返される（Ｓ−１２−６）。「計測」ボックスの「中止」ボタンが押されると、計測が中止する。計測が終了すると、図２６に示される画面表示がなされ、波形の確認がなされる（Ｓ−１２−７）。その後、そのデ−タの保存の必要性が判断され（Ｓ−１２−８）、保存が必要な場合はメニュ−「ファイル（Ｆ）」−「保存（Ｓ）」が選択されて、信号は保存され、当該被検者のデ−タリストに追加される（Ｓ−１２−９）。その後、保存が必要ない場合も含めて、計測がもう一度必要かどうかの判断がなされ（Ｓ−１２−１０）、必要ならば以上のステップを繰り返し、必要なければデ−タ計測の全ステップは終了する。その場合、画面表示は図２４の表示となるようにメニュ−選択がなされる。なお、図２６はチャンネルとして第２列目の行のチャンネルが選択された例を示す。
【００６２】
図２６において、解析デ−タ部の最下部にはスクロ−ルボックス２６１が移動するスクロ−ルバ−２６２がある。スクロ−ルボックス２６１はスクロ−ルバ−２６２の左右の両端間で移動可能なるもので、そのスクロ−ルボックスの幅ｗは時間スケ−ルを表す。そのスクロ−ルバ−２６２の左右の両端間の時間幅は計測時間を表し、したがって、表示されている波形は計測時間中に生じる波形の、スクロ−ルボックス２６１の時間スケ−ルｗに相当する一部分の拡大された波形である。これにより、操作者は、現在解析デ−タ部に表示されている波形が計測時間（スクロ−ルバ−２６２の幅）の中でどのくらいの時間（スクロ−ルボックス２６１の幅）かを示し、前記波形が計測時間の中で前半を示しているのか後半を示しているのか等を人目で把握することができるので、視認性を向上させることができる。また、前記スクロールボックス２６１の位置と解析データ部の波形を連動させているので、カーソルを前記スクロールボックス２６１に合わせてドラッグしながら移動させることで、解析データ部の表示領域を移動させて所定時間の波形を見るようにしてもよい。このようにすれば、所定時間の波形を簡単に確認する等、前記スクロ−ルバ−２６２を目次的に扱いで計測内容を詳細に確認することができる。
【００６３】
図１９は図１６のステップＳ−１０におけるアベレ−ジング処理のフロ−を示す。アベレージング処理では、各チャンネルで測定された各データのノイズを取り除くために各チャンネル毎に加算されてその平均値が計算される。この各チャンネルのアベレージング処理にあたり基準となる時間を設定するために以下の処理がなされ、そして各チャンネルのアベレージング処理が実行される。
【００６４】
まず、指定された被検者の生デ−タが読み込まれる（Ｓ−１０−１）。これは、図２４のリストデ−タ中の「デ−タの種類」が生デ−タとなっている行を選択することにより行われる。これによって、初期表示として、図２７に示される第１列目のチャンネルの心磁波形が表示される（Ｓ−１０−２）。次いで、表示チャンネルが選択される（Ｓ−１０−３）。図２７は第２列目のチャンネルが選択された例を示す。その後、設定チャンネルが指定される（Ｓ−１０−４）。これは操作領域部にある「アベレ−ジング条件」というボックス内の「チャンネル」のテキストボックスをクリックするか又は所望のチャンネルの波形を表示している領域をマウスでクリックすることによって行われる。この場合、そのテキストボックス内の三角ボタンをクリックすればチャンネル数字が大きくなり、逆三角ボタンをクリックすればチャンネル数字が小さくなる。図２７は指定されたチャンネルが第２列第２行目のチャンネルである例を示す。このチャンネル指定によって、アベレージング処理の基準時間となるチャンネルが特定される。このチャンネルを特定するにあたっては、最も典型的な、あるいは分かり易い波形を備えたチャンネルを選定するとよい。仮に、選択した行または列に良い波形を備えたチャンネルがなければ、再び（Ｓ−１０−３）のステップからやり直すことができる。
【００６５】
また、前記解析データ部の１つのチャンネルが特定されると、図２７に示されるようにしきい値カ−ソル２７１が指定されたチャンネルの位置に表示される。これによって、指定されたチャンネルを視覚的に確認することができる。図２７において、解析デ−タ部８０５の上部に左右方向に移動可能な３個のスライダ−カ−ソル２７３〜２７５が表示されているが、これらは図２７の表示画面が表示されると同時に自動的に表示される。
【００６６】
ステップＳ−１０−５として、アベレ−ジング条件であるしきい値、オフセット時間及びアベレ−ジング時間が「アベレ−ジング条件」というボックス内の対応するテキストボックスをクリックすることにより設定される。三角ボタンをクリックすれば数字が大きくなり、逆三角ボタンをクリックすれば数字が小さくなる。しきい値の設定は「しきい値」テキストボックスでしきい値を表す数字を選択することにより行われる。これによって、その数字に対応する位置にしきい値カ−ソル２７１が自動的に移動する。その場合の移動位置はカ−ソル線によって明確に確認することができる。「しきい値」テキストボックスのしきい値を表す数字の変更（選択）としきい値カ−ソル２７１の移動は互いに連動しており、したがってしきい値はしきい値カ−ソル２７１を移動させることによっても設定可能である。スライダ−カ−ソル２７３は波形の立上り部分がスライダ−カ−ソル２７１によって設定されたしきい値と一致した時点（基準時点）位置を指し示すものであり、その指し示す位置はカ−ソル線によって明確に確認することができる。このスライダ−カ−ソル２７３はしきい値の変更に応じて変わる基準時点位置と一致するように該基準時点位置に追従して移動する。「オフセット」テキストボックスでオフセット時間を表す数字を選択し、「時間」テキストボックスでアベレ−ジング時間を表す数字を選択すると、それに対応する位置にスライダ−カ−ソル２７４及び２７５が移動する。その場合の移動位置はそれらのスライダ−カ−ソルのカ−ソル線によって明確に確認することができる。スライダ−カ−ソル２７４及び２７５の移動は「オフセット」テキストボックス及び「時間」テキストボックスでの数字の選択と連動している。このため、オフセット時間及びアベレ−ジング時間の設定はスライダ−カ−ソルを移動させることによっても可能である。
【００６７】
ステッップＳ−１０−６ではアベレ−ジングモ−ドとしてアベレ−ジングを自動で行うべきか、マニュアルで行うべきかが設定される（Ｓ−１０−７）。その後、「設定」ボタンをクリックすると（Ｓ−１０−８）、加算処理が行われる。「キャンセル」ボタンをクリックすれば、アベレ−ジング条件はすべてキャンセルされる。加算処理としては、各チャンネルにおいて、しきい値を越えた時刻ｔ（基準時点）が探索され（Ｓ−１０−９）、次いで時刻ｔを中心に波形が表示され（つまり、ｔ時点が画面の中心に位置付けられようにｔ−５０ｍｓｅｃからｔ＋５０ｍｓｅｃまでの波形が表示され）（Ｓ−１０−１０）、そしてマニュアルモ−ドでキャンセルが選択されたかどうかの判断がなされ（Ｓ−１０−１１）、キャンセルが選択されたのでなければ波形の加算が実行される（Ｓ−１０−１２）。ステップＳ−１０−９からステップＳ−１０−１２までのステップは各チャンネル毎に加算回数だけ繰り返され、その後加算されたデ−タは加算回数で割算され（Ｓ−１０−１３）、アベレ−ジング処理が終了する。
【００６８】
このように、この図２７に示すアベレ−ジング処理の表示画面においては、操作領域部と解析データ部８０５の双方からアベレ−ジング処理の諸条件の入力／操作を行うことができるので、たとえば、大まかな条件を解析データ部で設定して、詳細な値を操作領域部で設定する等、操作者に多様な設定方法を提供すると共に条件設定の時間短縮を図ることができる。特に、この実施例では、解析データ部８０５の複数のチャンネルの中からカーソルを介して特定のチャンネルを指定することでアベレ−ジング処理の基準となるチャンネルとすることができから、誤操作が軽減され操作性が向上される。しかも、その指定されたチャンネル近傍にしきい値カ−ソルやスライダ−カ−ソルを表示することで視覚的な認識が可能である。更に、アベレ−ジング処理の諸条件の入力／操作を前記チャンネル近傍に配置されるしきい値カ−ソルやスライダ−カ−ソルで行うことができるので、操作者の目線移動を少なくして、波形表示に合わせて視覚的に設定できるから、誤操作を軽減して操作性を向上することができる。なお、これらの設定条件、たとえば最新の設定条件を記憶させて次回の設定条件として表示させたり、あるいは各設定条件に名前を付して記憶させ呼び出すようにしてもよい。
【００６９】
図２０は図１６のステップＳ−１１におけるデ−タ解析のフロ−を示す。デ−タ解析はいろいろな種類の波形や線図を表示して診断に必要な情報を得ようとするもので、図９のメニュ−を選択することによりいろいろな種類の波形や線図の画面を選択的に表示することができる。すなわち、「デ−タ解析（Ａ）」の「単一波形表示（Ｗ）」を選択すれば図２８に示される単一波形画面が（Ｓ−１１−２）、「デ−タ解析（Ａ）」の「重ね合せ波形表示（Ｍ）」を選択すれば図２９に示される重ね波形画面が（Ｓ−１１−３）、「デ−タ解析（Ａ）」の「グリッドマップ表示（Ｇ）」を選択すれば図３０に示されるグリッドマップ波形画面が（Ｓ−１１−４）、「デ−タ解析（Ａ）」の「等磁線図（Ｂ）」を選択すれば図３１に示される等磁線図画面が（Ｓ−１１−５）、「デ−タ解析（Ａ）」の「伝播時間図（Ｐ）」を選択すれば図３２に示される伝播時間図画面が（Ｓ−１１−６）、そして「デ−タ解析（Ａ）」の「時間積分図（Ｔ）」を選択すれば図３３に示される時間積分図画面が（Ｓ−１１−７）それぞれ表示され、また、「ファイル（Ｆ）」の「心磁システムの終了（Ｘ）」を選択すればシステムが終了する。
【００７０】
それぞれの画面において、操作領域部にあるラジオボタン（図面中の円形ボタン）をクリックすればそのクリックによって指定された波形又は線図の画面が代わって表示される。図２０において、分岐の部分を「メニュ−で分岐」とせずに「メニュ−又はラジオボタンで分岐」としたのはそのためである。したがって、この実施例によれば、前記図９のメニュ−を選択することなく、前記操作領域部にあるラジオボタンをクリックするだけで多様な解析データが得られるので、操作時間の短縮が図れるとともに、誤操作を軽減して操作性を向上することができる。
【００７１】
図２８〜３０において、「スケ−ル」ボックス内の磁束密度とあるのは、ゼロレベルを基準としたプラス側及びマイナス側のフルスケ−ルの値（単位はピコテスラ（ｐＴ））であり、その値はそのテキストボックスの三角ボタンをクリックすることにより開かれるプルダウンメニュ−で選択される。図２８〜３３において、「表示成分」ボックス内のラジオボックスをクリックすることにより法線成分の波形又は接線成分の波形を選択して画面表示することができる。
【００７２】
図２８ではチャンネルで選択された各チャンネルの波形が、解析データ部の左端をオフセット時間に合せて表示される。この解析データによれば、解析データ部に上下に配列表示される各チャンネルの波形の形状や大きさを比較することができる。同様に図２９では、前記図２８で上下に配列された波形が重ねて表示され、その波形の形状や大きさを比較することができる。また、図３０では、すべてのチャンネルが前記図２８、図２９と同様にオフセット時間を基準にして表示される。したがって、操作者は、必要によりチャンネルの数を選定して解析することができる。
【００７３】
図３１において、解析デ−タ部の右端には縦に細長の磁場強度指標ボックス３１０が配置されている。その磁場強度指標ボックスは互いに色の異なる１２個の区画に区切られている。これは、図３１に示される等磁線図画面の各島模様で示される磁場の強度範囲を色の種類で区別することで視角的な（色彩的な）認識性を向上させるようにしたものである。すなわち、その磁場強度指標ボックス３１０の長手方向の中心位置３１１は磁場強度がゼロの位置で、その中心位置よりも上方の区画を中心位置に近い順番に第１〜第６区画とそれぞれ呼ぶことにすれば、たとえば第１区画は０〜２ｐＴの磁場強度範囲に、第２区画は２〜４ｐＴの磁場強度範囲に、第３区画は４〜６ｐＴの磁場強度範囲に、第４区画は６〜８ｐＴの磁場強度範囲に、第５区画は８〜１０ｐＴの磁場強度範囲に、そして第６区画は１０〜１２ｐＴの磁場強度範囲にそれぞれ対応している。中心位置よりも下方の区画についてもまったく同じである。ただし、中心位置よりも上方の区画はプラス方向の磁場強度を、下方の区画はマイナス方向の磁場強度を表している。図３１に示される等磁線図は、磁場強度指標ボックス３１０内の磁場強度範囲と色との対応関係の定めにしたがい、磁場強度に応じて色分け表示される。なお、色として、磁場強度のプラス側を暖色系、マイナス側は寒色系とし、中心部を黄色とするようにしてもよい。これにより、磁場の強弱を色彩的に認識することができるので視認性を向上することができる。しかも、この実施例によれば、解析デ−タ部の近傍に磁場強度指標ボックス３１０を設けたので、比較対象の色、すなわち、マップに付された色と磁場強度指標ボックス３１０の所定の色とを目線移動を大きく移動させることなく比べながら確認することができるので、前記磁場の強弱のレベルと色との関係を明確に判断することができる。なお、この実施例では、磁場強度指標ボックス３１０を解析デ−タ部の右端に設けているが、解析デ−タ部の近傍であればよく、例えば、上部、下部、左側でもよい。
【００７４】
図３１において、「再構成パラメ−タ」ボックス内の「マップ数」とあるのは表示される等磁線図の数を、「最大値」とあるのは磁場強度指標ボックス３１０の両端部に相当する磁場強度を、「間隔」とあるのは磁場強度指標ボックス３１０内の各区画の長さに対応する磁場範囲を意味する。その値については、これを対応するテキストボックスの三角又は逆三角ボタンをクリックして選択することができる。
【００７５】
解析デ−タ部の最下段には参照チャンネルの心電波形と２個のマップ時刻選択用カ−ソル３１１及び３１２が表示されている。その２個のマップ時刻選択用カ−ソル３１１及び３１２間には間隔が同じ分割線が表示され、この線の数はマップ数選択によって選択されたマップの数と一致する。また、２個のマップ時刻選択用カ−ソル３１１及び３１２はその位置を独立に左右方向に移動可能で、その移動によってそのカ−ソル間の間隔が変わると、分割線の間隔も変わるが、分割線の間隔は常に等間隔である。もちろん、各分割線にカーソルを備えて１本１本別々に設定するようにしてもよい。図３１では、表示されている等磁線図の数は１６個であるが、これらの線図は心電波形上の分割線が位置する時点での線図であり、各マップについては、該マップがいつの時点のものであるかがわかるようにその時刻も表示される。
【００７６】
これにより、図２６で説明したと同様に、操作者は、現在解析デ−タ部に表示されているマップが解析時間（心電波形の幅）の中でどのくらいの範囲（２個のカ−ソル３１１と３１２の幅）を示し、前記マップが示す範囲が解析時間の中でどこの範囲なのか等を一目で把握することができるので、視認性を向上させることができる。また、前記マップが示す範囲を２個のカーソルをマウス操作で簡単に移動させることで設定できるので操作が容易である。更に、各分割線の間隔を自由に設定するようにすれば、疑問のある部分を密にして他の部分を疎にする等、操作者に多様な解析環境を提供できる。
【００７７】
また、操作領域部にある「電流方向」のチェックボックスをクリックしてチェック印を表示すると、等磁線図上には矢印が表示される。この矢印が表示された等磁線図をアロ−マップと呼ぶ（図示なし）。矢印については、その位置はチャンネルの位置（磁気センサの位置）、長さは磁場の強度、そして方向は磁場の方向を電流の方向に変換した場合のその電流の方向をそれぞれ示す。
【００７８】
図３２において、伝播時間の起点位置（図７のｔ１時点）の変更は参照波形上で移動するカ−ソル３２１の位置を変えることによって可能であり、カ−ソル３２１の位置の変更はそのカ−ソルをマウスを用いてドラッグすることにより可能である。
【００７９】
図３３には２つの時間積分図と１つの差分図が表示されている。伝播時間図の「差分表示」はチェックボックスをクリックしてチェック印を表示することで簡単に表示することができる。「差分表示」がチエックされると、参照波形上に４つのカ−ソル３３１〜３３４が現れ、更に図示のように上方左右に２つの時間積分図が、下方左側に差分図がそれぞれ表示される。２つの時間積分図は心磁波形を、参照波形上でカ−ソル３３１及び３３２並びにカ−ソル３３３及び３３４を用いてそれぞれ設定された１００ｍｓｅｃ〜１４０ｍｓｅｃ及び１８０ｍｓｅｃ〜２４０ｍｓｅｃという時間範囲に亘って積分した値にもとづくもので、それぞれの時間範囲はカ−ソル３３１及び３３２並びにカ−ソル３３３及び３３４をマウスを用いてそれぞれドラッグすることで変えられ得る。差分図は２つの時間積分図の差を表すものである。チェック印がない場合は、カ−ソルについては２個のカ−ソル（たとえばカ−ソル３３１及び３３２）だけが現れ、時間積分図については１つの時間積分図だけが表示される。もちろん、積分時間の変更はカ−ソルの位置を変えることによって可能である。このように、この実施例によれば、「差分表示」にチェック印をクリックすることにより、次の操作を促す２組のカーソルが表示されるので、操作の迷いを与えず操作時間の短縮が図れ、しかも、前記２組のカーソルをマウスで移動させることで時間範囲を簡単に設定することができるから操作性を向上することができる。
【００８０】
図２１はシステム調整のフロ−を示す。メニュ−選択によってフロ−は５つに分岐される（Ｓ−１５）。この場合、現在表示されているデ−タはそのまま保存され（Ｓ−１３）、図３４に示されるΦ−Ｖ特性曲線が表示されている（Ｓ−１４）。「デ−タ計測（Ｑ）」の「全自動調整（Ａ）」が選択された場合は、指定されたチャンネルについてＩｂｉａｓ及びＶＯＦＦという調整値が自動的に計算され（Ｓ−１６）（その自動計算は後述）、その計算された調整値はＦＬＬ回路に設定され（Ｓ−１７）、フロ−はステップＳ−１４に戻る。
【００８１】
「デ−タ計測（Ｑ）」の「ＶＯＦＦ調整（Ｖ）」が選択された場合は、指定されたチャンネルについてＶｏＦＦが計算され（Ｓ−１８）（ＶＯＦＦの計算は後述）、その後、フロ−は既述のステップＳ−１７に進む。「デ−タ計測（Ｑ）」の「マニュアル調整（Ｍ）」が選択された場合は、図１２に示されるマニュアル調整ダイアログボックスが開かれる（Ｓ−１９）。オペレ−タはＩｂｉａｓ及びＶＯＦＦという調整値をチャンネル毎に入力すると、その入力された調整値は受け付けられ（Ｓ−２０）、ダイアログボックス中の「ＯＫ」ボタンが押されることによって調整値がＦＬＬ回路にセットされる（Ｓ−２１）。ダイアログボックスはこれによって閉じられ（Ｓ−２２）、フロ−はステップＳ−１７に進む。
【００８２】
「デ−タ計測（Ｑ）」の「調整値ファイル（Ｆ）」が選択されると、その後サブプルダウンメニュ−によりフロ−は更に２つに分岐される（Ｓ−２２）。すなわち、サブプルダウンメニュ−の「開く（Ｏ）」を選択することによりファイルをオペレ−タに問合わせ、オペレ−タは図１１の内容を含むファイルを入力する（Ｓ−２３）。その調整値ファイルの内容はＦＬＬ回路に設定される（Ｓ−２４）。また、「デ−タ計測（Ｑ）」−「調整値ファイル（Ｆ）」−「上書き保存」又は「名前を付けて保存（Ａ）」を選択することによりステップＳ−２３と同様のことが行われ（Ｓ−２５）、その調整値を調整値ファイルに書き込むことができる（Ｓ−２６）。
【００８３】
図２２は図２１のステップＳ−１６における自動計算のフロ−を示す。
【００８４】
Ｓ−１５−１：全自動調整の場合はすべてのＳＱＵＩＤチャンネルについて１チャンネルずつ以下の処理が行われる。処理中のチャンネルをｃｈとする。
【００８５】
Ｓ−１５−２：計算されたバイアス電流及びオフセット電圧はＩＢＩＡＳ及びＶＯＦＦという名称のメモリに保持されるものとする。ＩＢＩＡＳ及びＶＯＦＦはチャンネルの数だけ値を保持することができ、初期値はすべて０とする。Φ−Ｖ特性曲線の、一時記憶される振幅をΔＶとする。
【００８６】
Ｓ−１５−３：各チャンネルｃｈに対して、バイアス電流Ｉｂを０から図３４の「走査パラメ−タ」ボックスで指定されたＩｂｉａｓまでΔＩｂｉａｓのステップで変化（走査）させ、Φ−Ｖ特性曲線の振幅ΔＶが最も大きくなるＩｂｉａｓをチャンネルｃｈの最適なバイアス電流ＩＢＩＡＳ（ｃｈ）とする。
【００８７】
Ｓ−１５−４〜８：バイアス電流０からＩｂまでの間で、Φ−Ｖ特性曲線の振幅ΔＶは以下の処理で求められる。チャンネルｃｈのＳＱＵＩＤに与えられる外部磁場Φを０から図３４の「走査パラメ−タ」ボックスで指定されたΦｅｘｔまでΔΦのステップで変えて（走査して）行き、Ａ／Ｄ変換された信号を保存してその信号の最大値Ｍｍａｘ及び最小値Ｖｍｉｎを求める。
【００８８】
Ｓ−１５−９〜１０：ここで、最大値Ｍｍａｘと最小値Ｖｍｉｎの差がΔＶより大きければＩＢＩＡＳ（ｃｈ）の値をＩｂｉａｓで、ＶＯＦＦ（ｃｈ）の値を最大値Ｍｍａｘと最小値Ｖｍｉｎの平均で、ΔＶを最大値Ｖｍａｘと最小値Ｖｍｉｎの差でそれぞれ置き換える。もし、ΔＶの方が最大値Ｍｍａｘと最小値Ｖｍｉｎの差より大きければ前の値を保持する。
【００８９】
Ｓ−１５−１１：以上の処理をバイアス電流ＩｂがＩｂｉａｓになるまで繰り返したときのＩＢＩＡＳ（ｃｈ）、ＶＯＦＦ（ｃｈ）がＳＱＵＩＤチャンネルｃｈの最適なバイアス電流及びオフセット電圧となる。
【００９０】
Ｓ−１５−１２：また、以上の処理をすべてのＳＱＵＩＤチャンネルに対して実行することにより全自動調整は終了する。
【００９１】
図２３は図２１のステップ１８におけるＶＯＦＦ調整フロ−を示す。
【００９２】
Ｓ−１８−１：ＶＯＦＦ調整の場合はすべてのＳＱＵＩＤチャンネルについて１チャンネルずつ次の処理が行われる。
【００９３】
Ｓ−１８−２〜６：チャンネルｃｈのＳＱＵＩＤに与えられる外部磁場Φを０から図３４の「走査パラメ−タ」ボックスで指定されたΦｅｘｔまでΔΦのステップで変えて（走査して）行き、Ａ／Ｄ変換された信号の最大値Ｖｍａｘと最小値Ｖｍｉｎの差を求める。
【００９４】
Ｓ−１８−７：ＳＱＵＩＤチャンネルｃｈの最適なオフセット電圧ＶＯＦＦ（ｃｈ）は最大値Ｖｍａｘと最小値Ｖｍｉｎの差として計算される。
【００９５】
Ｓ−１８−８：以上の処理をすべてのＳＱＵＩＤチャンネルに対して実行することによってＶＯＦＦ調整は終了する。
【００９６】
【発明の効果】
本発明によれば、操作および視認性を向上させた生体の磁場図表示装置を提供することができる。
【図面の簡単な説明】
【図１】本発明が実施される生体磁場計測装置の一実施例の概略構成図。
【図２】図１の生体磁場計測装置に用いられる磁気センサの配置構成を示す斜視図。
【図３】図１の生体磁場計測装置において用いられる、磁場の法線成分を検出する磁気センサ単体の斜視図。
【図４】図１の生体磁場計測装置において用いられる、磁場の法線成分を検出する磁気センサ単体の斜視図。
【図５】図１の生体磁場計測装置における磁気センサと被検者の胸部との位置関係を示す図。
【図６】図１の生体磁場計測装置において、各磁気センサ一において計測された健常者の生体磁場（心磁）の各成分の時間波形図。
【図７】図１の生体磁場計測装置において、健常者について計測された特定の２チャンネルの心磁の接線成分の時間波形を示す図。
【図８】図１の生体磁場計測装置において、デイスプレイ部に表示される表示画面の基本的なレイアウトを示す図。
【図９】図１の生体磁場計測装置において、デイスプレイ部に表示される表示画面のメニュ−バ−部における操作メニュ−を示す図。
【図１０】図１の生体磁場計測装置において、デイスプレイ部に表示される表示画面中の操作メニュ−として「リスト（Ｌ）」−「登録（Ｒ）」が選択された場合に開かれる被検者登録ダイアログボックスの内容を示す図。
【図１１】図１の生体磁場計測装置において、デイスプレイ部に表示される表示画面中の操作メニュ−として「リスト（Ｌ）」−「検索（Ｓ）」が選択された場合に開かれる検索ダイアログボックスの内容を示す図。
【図１２】図１の生体磁場計測装置において、デイスプレイ部に表示される表示画面中の操作メニュ−として「デ−タ計測（Ｑ）」−「マニュアル調整（Ｍ）」が選択された場合に開かれるマニュアル調整ダイアログボックスの内容を示す図。
【図１３】図１の生体磁場計測装置において、デイスプレイ部に表示される表示画面中の操作メニュ−として「デ−タ計測（Ｑ）」−「計測パネル（Ｐ）」が選択された場合に開かれる自動診断ダイアログボックスの内容を示す図。
【図１４】図１の生体磁場計測装置において、デイスプレイ部に図２５又は図２６の表示画面が表示されている場合に「ヒ−トフラッシュ」ボタンが押されたときに開かれるヒ−トフラッシュ操作ダイアログボックスの内容を示す図。
【図１５】図１の生体磁場計測装置において、計測中にデイスプレイ部に表示される計測プログレスバ−を示す図。
【図１６】図１の生体磁場計測装置において行われる全体の操作のフロ−を示す図。
【図１７】図１６の操作フロ−中の被検者選択ステップにおける被検者選択のフロ−を示す図。
【図１８】図１６の操作フロ−中のデ−タ計測ステップにおけるデ−タ計測のフロ−を示す図。
【図１９】図１６の操作フロ−中のアベレ−ジング処理ステップにおけるアベレ−ジング処理のフロ−を示す図。
【図２０】図１６の操作フロ−中のデ−タ解析ステップにおけるデ−タ解析のフロ−を示す図。
【図２１】図１の生体磁場計測装置において行われるシステム調整のフロ−を示す図。
【図２２】図２１のシステム調整フロ−中の自動計算ステップにおける自動計算のフロ−を示す図。
【図２３】図２１のシステム調整フロ−中のＶＯＦＦ調整ステップにおけるＶＯＦＦ調整フロ−を示す図。
【図２４】図１の生体磁場計測装置のデイスプレイ部に表示される被検者リスト画面を示す図。
【図２５】図１の生体磁場計測装置のデイスプレイ部に表示されるデ−タ計測画面を示す図。
【図２６】図１の生体磁場計測装置のデイスプレイ部に計測が終了したときに表示される波形確認画面を示す図。
【図２７】図１の生体磁場計測装置のデイスプレイ部に表示されるアベレ−ジング画面を示す図。
【図２８】図１の生体磁場計測装置のデイスプレイ部に表示される単一波形表示画面を示す図。
【図２９】図１の生体磁場計測装置のデイスプレイ部に表示される重ね波形表示画面を示す図。
【図３０】図１の生体磁場計測装置のデイスプレイ部に表示されるグリッドマップ表示画面を示す図。
【図３１】図１の生体磁場計測装置のデイスプレイ部に表示される等磁線図表示画面を示す図。
【図３２】図１の生体磁場計測装置のデイスプレイ部に表示される伝播時間図表示画面を示す図。
【図３３】図１の生体磁場計測装置のデイスプレイ部に表示される時間積分図表示画面を示す図。
【図３４】図１の生体磁場計測装置のデイスプレイ部に表示されるシステム調整画面を示す図。
【符号の説明】
１：磁気シ−ルドル−ム、２：被検者、３：ベッド、４：デユワ、５：自動補給装置、６：ＦＬＬ回路、７：増幅器・フイルタ−・増幅器、８：計算機、８−１：デイスプレイ部、８−２：キ−ボ−ド、８−３：マウス、２０−１〜２０−８、２１−１〜２１−８、２２−１〜２２−８、２３−１〜２３−８、２４−１〜２４−８、２５−１〜２５−８、２６−１〜２６−８及び２７−１〜２７−８：磁気センサ、１０、１０’及び１０”並びに１１、１１’及び１１”：コイル、１２、１２’及び１２”：ＳＱＵＩＤ、１３及び１４：センサ、３０：胸部、２６１：スクロ−ルボックス、２６２：スクロ−ルバ−、２７１：しきい値カ−ソル、２７３〜２７５：スライダ−カ−ソル、３１１、３１２、３２１及び３３１〜３３４：カ−ソル、８０１：タイトルバ−部、８０２：メニュ−バ−部、８０３：ツ−ルバ−部、８０４：被検者情報部、８０５：解析デ−タ部、８０６：操作領域部、８０８：ステ−タスバ−部、８０７−１：メッセイジバ−部、８０７−２：日時表示部、８０８〜８１４：アイコン。[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention provides a magnetic field of a living body suitable for measuring a magnetic field of a living body generated by an electric current in the living body, such as a neural activity of a living brain or a cardiac muscle activity of the heart.Figure display deviceAbout.
[0002]
[Prior art]
Conventionally, the distribution of a weak magnetic field generated from a living body is measured using a superconducting quantum interference device (SQUID), which is a magnetic sensor, and the position of an active current inside the living body is estimated from the measurement result. 2. Description of the Related Art A multi-channel biomagnetic imaging apparatus for imaging is known. Such a conventional example is disclosed in, for example, JP-A-4-319334 or JP-A-5-146416.
[0003]
[Problems to be solved by the invention]
The above-mentioned conventional example relates to the operation principle of a biomagnetic imaging apparatus, and does not disclose a technical problem or a solution for implementing the invention in its disclosure contents. Further, the above-mentioned conventional example relates to a biological activity current generated inside the brain, and does not specifically disclose other parts.
[0004]
DepartureMing is,Improved operation and visibilityTo provide a method for measuring a magnetic field of a living bodyWith the goal.
[0005]
[Means for Solving the Problems]
In order to achieve the above object, the present invention is a biological magnetic field map display device which measures a biomagnetic field of a subject at a plurality of positions and displays a magnetic field map as a processing result thereof, comprising input means and display. And a calculator having a unit, the computer displays a waveform of the measured biomagnetic field as a row or column of a plurality of channels in the analysis data section of the display unit, and analyzes the reference channel based on the other biological signals of the subject. In the lower part of the row or column of the plurality of channels of the biomagnetic field displayed in the section, receive an instruction for data analysis from the bathing means, and display the analysis magnetic field corresponding to the instruction in the analysis data section and the lower part thereof And a dividing line corresponding to the time of the displayed analysis magnetic field and a time selection car for moving the position of the dividing line on the reference channel displayed below the analysis magnetic field. Displays Le, receives the movement of the cursor by the input means, and displays the analysis field corresponding to the selected time after the movement.
[0006]
The above-mentioned biological magnetic field diagram display device is characterized in that the other biological signal of the subject is an electrocardiographic waveform, an electroencephalogram, a blood flow waveform or a blood pressure waveform of the subject.
The analysis magnetic field displayed in the analysis data section of the above-described biological magnetic field diagram display device is a magnetic isochronogram, a propagation time diagram, or a time integration diagram.
In the above-described magnetic field diagram display apparatus for a living body, the analysis magnetic field displayed in the analysis data section is an isomagnetic map, the time selection cursor is displayed at two places, and the dividing line is displayed between the two cursors. It is displayed at intervals.
In the above-described magnetic field diagram display apparatus for a living body, the analysis magnetic field displayed in the analysis data section is an isomagnetic map, the time selection cursor is set to each division line, and the division lines can be individually moved. There is a feature.
The biological magnetic field diagram display apparatus described above is characterized in that the analysis magnetic field displayed in the analysis data section is a propagation time diagram, and the starting position of the propagation time can be changed by moving the time selection cursor. .
[0007]
In the above-described magnetic field diagram display apparatus for a living body, the analysis magnetic field displayed in the analysis data is a time integral diagram, and the computer receives an instruction for displaying a difference from the input unit, and displays four differential displays on the reference channel. A time selection cursor is displayed, and two time integration diagrams obtained by time integration in a time range corresponding to the four cursor positions and a difference diagram thereof are displayed in the analysis data section.
The reference channel of the above-mentioned magnetic field map display of a living body is characterized in that a predetermined channel can be selected from a plurality of channels and displayed.
[0008]
BEST MODE FOR CARRYING OUT THE INVENTION
FIG. 1 shows a schematic configuration of an embodiment of a biomagnetism measuring apparatus according to the present invention. In order to remove the influence of environmental magnetic noise, the biomagnetic measuring device is installed in the magnetic shield 1. The subject 2, which is a subject composed of a living body, is measured while lying on the bed 3. The body surface of the subject (in the case of the chest, generally parallel to the chest wall) is assumed to be substantially parallel to the surface of the bed 3, and this surface is defined by the xy plane of the rectangular coordinate system (x, y, z). It is assumed that they are parallel. Although the chest of the subject is curved and inclined, it is assumed to be substantially parallel for ease of explanation.
[0009]
Above the chest of the subject 2, a dewar 4 filled with liquid He as a refrigerant is arranged, and the dewar is connected to a superconducting quantum interference device (SQUID = Superconducting Quantum Interference Device) and the SQUID. And a plurality of magnetic sensors including a coil. The liquid He is continuously replenished from the automatic replenishing device 5 outside the magnetic shield 1.
[0010]
The output of the magnetic sensor outputs a voltage that has a specific relationship with the intensity of the biomagnetic field (which can also be considered as a magnetic flux density) generated from the subject 2 and detected by the detection coil, and the output is FLL (Flux). (Locked loop) circuit 6. The FLL circuit 6 cancels a change in the biomagnetic field (biomagnetism) input to the SQUID via a feedback coil so as to keep the output of the SQUID constant (this is called a magnetic field lock). By converting the current flowing through the feedback coil into a voltage, a voltage output having a specific relationship with a change in the biomagnetic field signal can be obtained. Since the detection is performed via the feedback coil as described above, a weak magnetic field can be detected with high sensitivity.
[0011]
The output voltage is input to an amplifier / filter / amplifier (AFA) 7 whose output is sampled, A / D converted, and taken into a computer 8.
[0012]
The computer 8 is composed of a personal computer, 8-1 is its display unit, 8-2 is a keyboard, and 8-3 is a mouse. The mouse 8-3 is used to move a cursor on the screen to select a processing target. This operation can also be performed by operating the keyboard. Input gain of AFA7 (Igain) And output gain (Ogain) Is adjustable, and the AFA 7 is a low-pass filter (LPF) that passes a frequency signal lower than the first reference frequency, and passes a frequency signal higher than a second reference frequency lower than the first reference frequency. Includes a high-pass filter (HPF) and a notch filter (BEF) that cuts off the commercial power frequency. The computer 8 can perform various types of processing, and the processing results can be displayed on the display unit 8-1. Note that the computer 8 shown in FIG. 1 shows one embodiment, and the present invention is not limited to this. For example, a device provided with a display having a touch panel, or a device using another coordinate pointing device, for example, a trackball or a joystick, instead of a mouse, may be used. In some cases, the computer may be connected via a public telephone line.
[0013]
For example, a DC SQUID is used as the SQUID. When an external magnetic field is applied to the SQUID, a direct current (I) is applied to the SQUID so that a voltage (V) corresponding to the external magnetic field is generated.bias) Is washed away. When the external magnetic field is represented by a magnetic flux Φ, a characteristic curve of V with respect to Φ, that is, a Φ-V characteristic curve is given by a periodic function. Prior to the measurement, the offset voltage (VOFF) Is adjusted to make the DC voltage of the Φ-V characteristic curve zero level. Further, when the input of the AFA 7 is zero, the offset voltage (AOFF) Adjustment is performed.
[0014]
When a large magnetic field is applied to the SQUID from the outside, the magnetic field is trapped by the SQUID, and the SQUID cannot perform its normal operation. In that case, the SQUID can be heated to a normal state once, and then the heating can be stopped to remove the trapped magnetic field. In this case, the heating operation of the SQUID is called a heat flash.
[0015]
FIG. 2 shows an arrangement configuration of the magnetic sensor. The detection coil of the magnetic sensor has a coil for detecting a tangential component of the biomagnetic field (a component substantially parallel to the biological surface, ie, the xy plane) and a normal component of the biomagnetic field (a component orthogonal to the biological surface, ie, the xy plane). There is a coil that detects As the coil for detecting the tangential component of the biomagnetic field, two coils whose coil surfaces are respectively oriented in the x direction and the y direction are used, and as the coil for detecting the normal component of the biomagnetic field, the coil surface is in the z direction. Are used. Plural magnetic sensors 20-1 to 20-8, 21-1 to 21-8, 22-1 to 22-8, 23-1 to 23-8, 241-1 to 24-8, 25-1 to 25 As shown in FIG. 2, -8, 26-1 to 26-8, and 27-1 to 27-9 are arranged in a matrix on a biological surface, that is, a surface substantially parallel to the xy plane. Although the number of magnetic sensors may be arbitrary, in FIG. 2, since the matrix of magnetic sensors is composed of 8 rows and 8 columns, the number of magnetic sensors is 8 × 8 = 64. As shown in FIG. 2, each magnetic sensor is arranged such that its longitudinal direction coincides with the direction perpendicular to the living body plane, that is, the xy plane (z direction). In this embodiment, the bed surface and the XY plane of the sensor are parallel to each other. However, in order to increase the measurement accuracy, it is better to approach the body and to tilt the bed. However, since the human body, which is the subject, is constantly moving, if the human body is brought into close contact with the human body, this movement will move the detection unit, and it will be difficult to perform highly accurate detection.
[0016]
FIG. 3 shows a configuration of each of the magnetic sensors for detecting a normal component Bz of the biomagnetic field. In the figure, a coil made of a superconducting wire (Ni-Ti wire) is arranged so that its coil surface faces the z direction. This coil is composed of a combination of two coils 10 and 11 having opposite directions. The coil 10 closer to the subject 2 is a detection coil, and the coil 11 farther away is a reference coil for detecting external magnetic field noise. You. External magnetic field noise originates from a signal source that is farther than the subject, so that the noise signal is detected by both the detection coil 10 and the reference coil 11. On the other hand, the magnetic field signal from the subject is weak. Therefore, the biomagnetic signal is detected by the detection coil 10, but the reference coil 11 is hardly sensitive to the biomagnetic signal. For this reason, since the detection coil 10 detects the biomagnetic field signal and the external magnetic field noise signal, and the reference coil 11 detects the external magnetic field noise signal, the difference between the signals detected by the two coils is used to obtain the S / N ratio. Measurement of a high biomagnetic field becomes possible. These coils are connected to the input coil of the SQUID via a superconducting wire of the mounting board on which the SQUID 12 is mounted, whereby the component Bz in the normal direction of the detected biomagnetic signal is transmitted to the SQUID.
[0017]
FIG. 4 shows a configuration of each of the magnetic sensors for detecting tangential components Bx and By of the biomagnetic field. In the figure, a planar coil is used in a sensor for detecting a biomagnetic component in a tangential direction. That is, the detection coils 10 ′ and 10 ″ and the reference coils 11 ′ and 11 ″ consist of planar coils, which are respectively arranged on first and second planes that are spaced apart in the z-direction. These coils are connected to the input coils of the mounting boards of SQUIDs 12 'and 12 "in the same manner as those for the normal component. The Bx component detection sensor 13 and the By component detection are provided on two mutually orthogonal surfaces of a quadrangular prism. Sensor 14 is attached, thereby forming a sensor capable of detecting the Bx component and the By component.
[0018]
The tangential components Bx and By may be obtained by partially differentiating the normal component Bz obtained by the magnetic sensor of FIG. 3 with respect to x and y in addition to the detection using the magnetic sensor shown in FIG. Good. In this case, one magnetic sensor can detect and measure both the tangential components Bx and By and the normal component Bz.
[0019]
FIG. 5 shows a positional relationship between the magnetic sensor and the chest 30 which is a measurement section of the subject 2. The points shown represent the intersections of the rows and columns on the matrix shown in FIG. 2, that is, the measurement points of the subject 2, that is, the measurement positions. Each of these measurement positions is also called a channel. As can be seen from the figure, in this embodiment, the height direction of the subject 2 is defined as the y direction, and the lateral direction of the subject 2 is defined as the x direction.
[0020]
FIG. 6 shows the measurement results of the biomagnetic field at the respective measurement positions shown in FIG. This measurement result shows a time-varying biomagnetic field waveform obtained by performing the above-described processing based on the signals detected by the magnetic sensors corresponding to the respective measurement positions, for each corresponding channel in the matrix. Things. In this embodiment, since each channel is provided at a position where the magnetic field generated by the heart muscle can be detected, the waveform in FIG. 6 shows a magnetocardiogram waveform. Note that a waveform obtained by measuring a magnetic field generated from the heart muscle is called a magnetocardiogram waveform. As shown in FIG. 6, when the measurement data measured for each channel is displayed corresponding to the position for each channel, this is called a grid map display. FIG. 6 shows the waveform of the measurement result of the magnetocardiogram for a certain healthy person. Here, (a) shows the magnetocardiogram waveform of the tangential component Bx, (b) shows the magnetocardiogram waveform of the tangential component By, and (c) shows the magnetocardiogram waveform of the normal component Bz.
[0021]
FIG. 7 shows a magnetocardiogram waveform of a tangential component Bx measured for a certain healthy person by focusing on two specific channels. The solid line shows the magnetocardiogram waveform of one channel, and the dotted line shows the magnetocardiogram waveform of another channel. Time period T during which the ventricles of the heart depolarized1In other words, the time of each waveform peak in the systolic QRS wave is tQ, TRAnd tSRespectively. The time period of the T wave, which is the heart repolarization process (diastole), is T2It is shown as When displaying the measured data, in addition to the grid map display, data of a single channel may be displayed row by row or column by row, or data of all channels or a plurality of channels may be displayed. The data may be superimposed and displayed. The former is called a single waveform display, and the latter is called an overlap waveform display.
[0022]
With respect to the obtained magnetocardiogram waveform data, an averaging (averaging) process is performed, an isomagnetic diagram, a time integration diagram, and the like are created, and the results can be displayed. For example, referring to FIG. 7, the rising portion of the QRS wave has a threshold value SLA predetermined time (tOFF), And the time t1Predetermined time T from3Time t2The data up to is added a predetermined number of times. This is averaging, which is a predetermined time T3Is the averaging time, tOFFIs called an offset time. The magnetocardiogram waveform data may be integrated over a predetermined time range. A map formed by connecting points (channels) having the same time integration value is called a time integration diagram. A map formed by connecting points having the same magnetocardiogram waveform signal value is called an isomagnetic diagram. Since each channel is roughly set, a diagram more suitable for diagnosis is created by setting an interval between magnetic lines, that is, a magnetic field intensity difference in advance, and linearly interpolating between each channel to draw magnetic lines. be able to. Referring to the magnetocardiogram waveform shown in FIG.1From the peak position of the QRS wave (t)RThe time up to the point) is called a propagation time, and a map formed by connecting points having the same propagation time is called a propagation time diagram. Threshold SLCan be changed. t1For the time point, the rising portion of the QRS waveLIs determined on the basis of the time at which the threshold value SLAlternatively, the determination may be made based on a point in time when the value matches. t1The time point is further the peak position time point of the QRS wave (tR(Time point) may be detected, and the time point may be determined as a reference. The measured biomagnetic signal is generated by an electrophysiological phenomenon in a living body, and its source is approximated by a current dipole model. The current dipole of the magnetic field source is synthesized and displayed on the isomagnetic diagram, which is called a magnetic field source display.
[0023]
A series of operations from registration of the subject to measurement of the data of the registered subject and analysis of the measured data are performed on a display screen displayed on the display 8-1. It is performed while watching. Therefore, prior to the description of the series of operations, the layout of the display screen will be described first.
[0024]
FIG. 8 shows a basic layout of a display screen displayed on the display 8-1 in FIG. The upper part of the display screen is occupied by a title bar section 801, a menu section 802, and a toolbar section 803 in which icons are arranged in order from the top. Each of the above units can be considered as a display area or area. These arrangements are commonly displayed on the display screen for other processing purposes, such as registration and readout of a subject, measurement of a magnetic field, and analysis of measurement data. This increases ease of use and reduces the time for measurement and processing.
[0025]
The central portion of the display screen is occupied by a subject information section 804 arranged in order from left to right, an analysis data section 805 for displaying analysis data such as a diagram and a waveform, and an operation area section 806. Have been. The lower part is occupied by a status bar 807, which is a message bar part 807-1 arranged on the left side for displaying a guide message for the next operation and a date and time display arranged on the right side thereof. Unit 807-2. Note that the message bar section 807-1 and the date and time display 807-2 may be one display area.
[0026]
On the display screen in this embodiment, a title bar section 801 that always displays the name of this system at the top, a menu bar section 802 for performing basic operations of the system, and a menu bar section 802 are used. Since the toolbar section 803 capable of performing frequent operations is arranged, the user does not have to search for an operation area every time the display screen changes, and if the user always looks at the top of the display screen, the currently operating system can be operated. You can know. Moreover, since the upper part of the display screen is the first place to look at, assuming that a person reads a sentence, it is possible to provide basic items in the operation of this system at the top. In a natural way, usability has been improved. In the center of the display screen, an analysis data section 805, which is a main body of the display screen, is provided at the center to improve the visibility, and on the right side, an operation area unique to the display screen is provided. By providing the unit 806, the arrangement is the same as the arrangement of the display screen and the operation unit in the right hand operation, so that the operation can be performed without discomfort. Therefore, even if this display screen is adopted as a display screen with a touch panel, the right hand operating the operation area section 806 does not disturb the analysis data section 805. Similarly, since the subject information section 804 having only a confirmation function is provided on the left side of the analysis data section 805, the operation can be performed while always confirming the patient. In addition, since the left position is the farthest position in the right hand operation, even if adopted for a display screen with a touch panel, the visibility of the display is not affected.
[0027]
The analysis data section 805 and the operation area section 806 are replaced by the subject list and the data list of the subject only when they are displayed (FIG. 24). When the subject list screen (FIG. 24) is displayed, the subject information section 804 always displays information on the subject on which the cursor is placed in the subject list in the screen. When analysis data such as a diagram or a waveform is displayed in the analysis data section (FIGS. 25 to 34), the subject of the subject who obtained the displayed analysis data is displayed. Information is always displayed. Thus, the relationship between the displayed analysis data and the subject from which the analysis data has been obtained can be clearly known. As described above, on the display screen of this system, the subject information section 804 is always displayed at the fixed position (left side) on the display screen, similar to the menu section 802, so that the user can display the screen. There is no need to search the subject information area every time the screen changes, and it is possible to know by always looking at a predetermined position (left side) of the display screen.
[0028]
The name of the frame is displayed in the title bar, specifically, the name "Multichannel MCG System" (FIGS. 24 to 34). In the operation area, operation elements such as buttons and text boxes are arranged. The menu section is used to select an operation menu. The menus are "File (F)", "Edit (E)", "List (L)", "Data measurement (Q)", " Data analysis (A) "and" Help (H) ", which are arranged in the order of operation.
[0029]
FIG. 9 shows the contents of the respective operation menus in the menu section on the display screen, and the contents of these menus are displayed as pull-down menus by clicking the corresponding menu buttons. For this reason, when the operation menu is not required, only the keywords for calling the respective menus are compactly displayed in the menu section, so that the display necessary for each operation such as the analysis data section and the operation area section is performed. The area can be set wide. When an operation menu is required, the keywords arranged in accordance with the operation procedure can be displayed by selecting the keyword from a menu bar portion to give an operation instruction. At this time, since the keywords are arranged according to an arrangement of characters (from left to right), the operation can be instructed in a natural manner.
[0030]
The “File (F)” pull-down menu opens a page layout dialog box (not shown), sets “Page layout (U)” to set the page layout, and provides a preview for the preview before printing. The items include “Preview (V)”, “Print (P)” for printing data, and “End (X) of the magnetocardiographic system” for terminating the Multichannel MCG System.
[0031]
The "edit (E)" pull-down menu includes an item "delete (D)". The measurement data of the subject in the subject list on the subject list screen (FIG. 24) is displayed in the data list on the same screen, but the item "Delete (D)" is displayed in the data list. -To delete the data on which the cursor is placed in the data list, the data must be selected in advance. When this menu is clicked, a confirmation dialog box for confirming whether or not deletion is permitted is opened. When the deletion is required, a button "OK" in the dialog box is clicked, and when the user wants to cancel the deletion, a button "Cancel" is clicked. By displaying and operating the confirmation dialog box, it is possible to reduce the erroneous operation of erroneously deleting the measurement data of the subject.
[0032]
The “List (L)” pull-down menu includes items of “Registration (R)”, “List (L)”, “Delete (D)”, “Search (S)”, and “Cancel (X)”. When the button "Register (R)" in the pull-down menu is clicked, a subject registration dialog box shown in FIG. 10 is opened. This is used when registering data about the subject, and the items that can be registered are the registration date, the subject's ID number up to a predetermined digit, name, date of birth, height, weight, and gender. , Classification indicating disease classification information, and subject comments. When the "Register" button in the dialog box is clicked, the data entered for registration is registered, and all of the input boxes are cleared so that re-entry is possible. Is clicked, all its input boxes are cleared, and if the button "Exit" is clicked, the subject registration dialog box is closed.
[0033]
When "List (L)" in the pull-down menu is clicked, a subject list screen (FIG. 24) is displayed. "Delete (D)" in the pull-down menu is for deleting the subject on which the cursor is placed in the subject list on the subject list screen (FIG. 24). When the button is clicked, a confirmation dialog box for confirming whether or not to delete the file is displayed in the same manner as when "Delete (D)" in the pull-down menu is clicked. Then, the button "OK" in the dialog box is clicked, and when the deletion is to be canceled, the button "Cancel" is clicked. When a subject is deleted, all data in the data list for that subject is also deleted.
[0034]
When "Search (S)" in the pull-down menu is clicked, a search dialog box shown in FIG. 11 is displayed. The subject and data are searched, and only the searched subject and data are displayed on the subject list screen. The subjects and data to be searched are all subjects and all data. The data is retrieved by giving the subject name, registration date, gender, data type, measurement date, diagnosis result, comment, laboratory technician, body position, etc. as keywords. A composite search combining a plurality of the above keywords can be performed. "Release (X)" in the pull-down menu is used to release the display of only the retrieved subjects and data and to display all subjects and all data. The pull-down menu of “Data measurement (Q)” is “Adjustment value file (F)”, “Fully automatic adjustment (A)”, “VOFFAdjustment (V), “Manual adjustment (M)”, “Adjustment value file (F)”, “Measurement panel (P)”, “Automatic waveform diagnosis (D)”, and “AFA offset adjustment (O)” including.
[0035]
When "Adjustment file (F)" is clicked, a sub pull-down menu containing "Open (O)", "Overwrite save (S)", and "Save as (A)" is displayed. . "Open (O)" in the sub pull-down menu displays the system adjustment screen (FIG. 34), opens the specified adjustment value file, and adjusts the system adjustment value, ie, I.biasAnd VOFFIs used to set the system. “Save (S)” is used to open a confirmation dialog box and overwrite the current adjustment value to the currently opened adjustment value file. "Save as" (A) "is used to change the name of the current adjustment value data and save it in another adjustment value file.
[0036]
When "Fully automatic adjustment (A)" is selected in the pull-down menu, a system adjustment screen (FIG. 34) is displayed, and the bias current I is adjusted according to the flow of FIG.biasAnd the offset voltage V of the FLL circuit 6OFFIs automatically adjusted. "VOFFWhen "Adjustment (V)" is selected, a system adjustment screen (FIG. 34) is displayed, and the offset voltage V of the FLL circuit 6 is set in accordance with the flow of FIG.OFFIs automatically adjusted. When "manual adjustment (M)" is selected, a manual adjustment dialog box shown in FIG. 12 is opened. The operator selects a channel using a scroll bar and a mouse, and sets a bias current IbiasAnd offset voltage VOFFTo change. If there is no problem with the input value, the value is set by clicking the "OK" button. If the "Cancel" button is clicked, the changes are invalidated and the dialog box is closed.
[0037]
When the "measurement panel (P)" is clicked, a data measurement screen including a grid map (FIG. 25) is displayed. When "automatic waveform diagnosis (D)" is clicked, an automatic diagnosis dialog box shown in FIG. 13 is opened. The automatic waveform diagnosis automatically checks the amplitude, median value, and period of the Φ-V characteristic curve in the display screen of FIG. 34, and updates the display of the Φ-V characteristic curve when displaying or displaying the same curve. At the time of updating, the operator is notified of a channel outside the specified range as an inappropriate state for measurement. If a channel in a state unsuitable for measurement is detected, an error dialog box is opened and an error message is displayed to notify the operator, but the Φ-V characteristic curve of the channel is changed to another color. May be displayed to indicate an inappropriate state. When the minimum amplitude value is checked, when the difference between the maximum value and the minimum value in the Φ-V characteristic curve is smaller than the minimum amplitude value, it is determined that the state is inappropriate for measurement. When the median value is specified, it is assumed that the case where the absolute value of the average of the maximum value and the minimum value is larger than the specified median value is inappropriate for measurement. When the period is specified, the measurement is not performed when the period of Φ in the Φ-V characteristic curve is out of the range specified by the first text box (lower limit) and the second text box (upper limit). It is assumed that it is in an appropriate state. To enable automatic diagnosis of the minimum amplitude, median and period, click the check box to the left of each item with the mouse so that the X mark is displayed, and enter the desired value in the corresponding text box. Should be input. The parameters of the automatic diagnosis input in this manner are enabled by pressing the "OK" button, and the dialog box is closed. When the "Cancel" button is pressed, the entered parameters become invalid and the dialog box is closed. "AFA offset adjustment (O)" is the offset voltage A of AFA7.OFFWhen "AFA offset adjustment (O)" is clicked, a data measurement screen (FIG. 25) belonging to the single waveform display is displayed.
[0038]
The pull-down menu of "Data analysis (A)" is "Averaging (A)", "Single waveform display (W)", "Overlap waveform display (M)", "Grid map display (G)". , "Equimagnetic field diagram (B)", "Time integration diagram (T)", "Propagation time diagram (P)", "Magnetic field source display (S)", "Line mode (L)", and "Filling" Mode (F) ".
[0039]
When "Averaging (A)" is clicked, the averaging screen (FIG. 27) is displayed, and when "Single waveform display (W)" is clicked, the single waveform display character screen (FIG. 28) is displayed. Clicking on “Overlay waveform display (M)” displays the overlay waveform display screen (FIG. 29), and clicking on “Grid map display (G)” displays the grid map display screen (FIG. 30). A check mark (x) indicating that the selection has been made is displayed on the left side of. Clicking on the “isomagnetic diagram (B)” displays the isomagnetic diagram screen (FIG. 31), and clicking on the “time integral diagram (T)” displays the time integral diagram screen (FIG. 32). When "Propagation time diagram (P)" is clicked, a propagation time diagram screen (FIG. 33) is displayed, and a check mark (x) indicating that the selection has been made is displayed on the left side of the menu. When "magnetic field source display (S)" is clicked, an inverted triangular mark is displayed on the left side of the menu, and a magnetic field source approximated by a current dipole is overlapped when displaying an isomagnetic diagram. Displayed (not shown). When the line mode (L) is selected, the lines between the lines are displayed without filling in the screen of the isomagnetic diagram, the time integration diagram, the propagation time diagram, and the magnetic field source display, and the "filling mode ( When “F)” is selected, the image is displayed in a state where the space between lines is filled. To indicate the current mode, a check mark is displayed on the left side of the menu to indicate the selection.
[0040]
The "Help (H)" pull-down menu includes "Contents (C)", "Search by keyword (S)", and "Version information (A)". Used to search for topics in a keyword, open a version dialog box.
[0041]
In the toolbar section 803, “examinee registration” (808), “examinee list” (809), “print” (810), “preview” (811), and “system adjustment” (812) , "Data measurement" (813), and "data analysis" (814). Although these are not shown, they are linked to menu functions, and the most frequently used items can be selected from the pull-down menu items. That is, “subject registration” (808) is “registration (R)” of “list (L)”, and “subject list” (809) is “list (L)” of “list (L)”. “Print” (810) is “Print (P)” of “File (F)”, “Preview” (811) is “Preview (V)” of “File (F)”, and “System”. "Adjustment" (812) is "Manual adjustment (M)" of "Data measurement (Q)", and "Data measurement" (813) is "Measurement panel (P)" of "Data measurement (Q)". )) And "data analysis" (814) correspond to "grid map display (G)" of "data analysis (A)", respectively. As for the menu linked to the icon, the menu can be selected by simply clicking the icon. Therefore, frequently-used operations can be easily operated by simply clicking the icon adjacent to the analysis data section. Therefore, compared to the operation of the menu section, the icons can be easily recognized in a shorter time. Can be operated. The icon may be selected by the user, or may be automatically displayed on the toolbar 803 according to the frequency of use (the number of times).
[0042]
Next, a series of operations from registration of the subject to data measurement of the registered subject and analysis of the measured data, including the adjustment operation of the system. Will be described with reference to FIGS.
[0043]
FIG. 16 shows the flow of the entire operation. When the power of the computer 8 is turned on (S-1), the operating system is started up, and the program start icon is displayed on the display section 8-1. Is performed (S-2). When a Multichannel MCG System program icon is selected from the icons (S-3), the subject list screen shown in FIG. 24 is displayed instead (S-4).
[0044]
In the system according to this embodiment, a subject list screen shown in FIG. 24 is displayed as an initial screen at the start of the system. The reason for this is that the relationship between the subject and the measurement or analysis data of the subject is extremely important, and therefore, this system manages the subject information as keywords. In other words, measurement data and analysis data cannot be managed without subject information. Therefore, in this system, on the subject list screen, first, when the subject is registered or registered, the subject is specified. If there is, specify the target data. It is to be noted that a display screen for waiting time at system startup may be provided prior to the subject list screen, and a display screen serving as a table of contents of the system may be provided.
[0045]
To describe the subject list screen shown in FIG. 24, the left side is occupied by the subject information section. The subject list is displayed at the top of the entire right side, and the data list is displayed at the bottom. Items displayed in the subject information section are the same as those described with reference to FIG. Items in the subject list include ID (subject ID number), name, registration date (date of data registration), number of measurements (number of times data measurement was performed), date of birth , Age, height, weight, and comments (comments on the subject). The subject list can be scrolled with a vertical scroll bar, and the subject list items can be scrolled with a horizontal (horizontal) scroll bar. . The row of the selected subject is highlighted.
[0046]
The items in the data list relating to the selected subject include ID, data type (raw data or averaging), and sampling interval (data measurement was performed). Time signal sampling interval in milliseconds), sampling time (seconds), classification (disease classification information), Date and Time (date and time when data measurement was performed), comment (data Comments on data). The data list can be scrolled with a vertical scroll bar, and the data list item can be scrolled with a horizontal (horizontal) scroll bar. The selected line of data is highlighted.
[0047]
According to this subject list screen, information of each subject is displayed on the subject list in one line. As a result, the information of the subjects arranged vertically can be clearly separated, so that the discriminability can be improved. For example, the erroneous operation of selecting another subject by mistake can be reduced. it can. The information of each subject can be scrolled by a horizontal (horizontal) scroll bar, and the information of the selected subjects is vertically arranged for each item in a vertically long subject information section. Therefore, visibility is not impaired. In this case, the visibility may be further improved by enabling the data list items of each subject to move forward and backward (left and right). Further, by displaying the information of each subject on one line, many subjects can be seen at a time, so that the number of times of scrolling with the vertical scroll bar can be reduced. Further, the data can be displayed in the lower data list by a simple operation of simply placing the cursor on the subject list to be selected from the subject list and clicking on the subject list to be selected. In addition, since the subject list and the data list are arranged vertically, movement of the line of sight can be reduced, so that the relevance can be easily recognized. The size of the data list can be freely changed by a simple operation of moving the cursor to the upper part of the area and dragging the data list. Therefore, the size of the data list can be freely adjusted according to the number of data lists. Can be set.
[0048]
In step S-5, a desired subject line is selected from the subject list on the subject list screen. In the case of averaging processing described later, a row of raw data in the data list is always selected. Thereafter, the flow is branched into four by the menu (S-6). According to one of the branches, a submenu "End of magnetocardiography (X)" of the menu "File (F)" is selected. In this case, an end process such as closing a window is performed (S-). 7), whereby the system is shut down (S-8). Thereafter, the power of the computer 18 is turned off (S-9), and all the operations are terminated.
[0049]
According to the rest of the branches, the averaging process (S-10), the data analysis (S-11), and the data measurement (S-12) are performed. The averaging process can be executed by selecting the submenu "Averaging (A)" of the menu "Data analysis (A)". The data analysis is performed in the menu "Data analysis (A)", such as "single waveform display (W)", "overlap waveform display (M)", "grid map display (G)", "etc." This can be executed by selecting any one of the submenus of “magnetic diagram (B)”, “time integral diagram (T)”, “propagation time diagram (P)”, and “magnetic field source display (S)”. Further, data measurement can be executed by selecting a submenu "measurement panel (P)" of a menu "data measurement (Q)". After the completion of steps 10, 11 and 12, the flow returns to step S-4. The details of the subject selection in step S-5, the averaging process in step S-10, the data analysis in step S-11, and the data analysis in step S-12 are shown in FIGS. This is described in more detail below in connection with
[0050]
FIG. 17 shows a flow of subject selection in step S-5 in FIG. In the case of subject selection, the flow is branched into four depending on menu selection or subject selection. One of the branches is a case where the subject selection is completed by designating and selecting the subject (S-5-1). According to another branch, the submenu "search (S)" of the menu "list (L)" is selected. Thereby, the search dialog box shown in FIG. 11 is opened (S-5-2), and the subject search condition is input using this dialog box (S-5-3). As a result, the subject is searched (S-5-4), and based on this, the display content of the subject list on the subject list screen shown in FIG. 24 is changed (S-5-5). . According to yet another branch, a submenu "release (X)" of the menu "list (L)" is selected. In this case, the selected subject is returned to the list of all subjects (S-5-6), and the display contents of the subject list are changed (S-5-5). According to the remaining one of the branches, the submenu "Registration (R)" of the menu "List (L)" is selected. In this case, the subject registration dialog box shown in FIG. 10 is opened (S-5-7), and subject information is input (S-5-8). In these steps, the end of the input is determined until all the subjects have been input (S-5-9), and when the input is completed, the subject list is updated (S-5). -5).
[0051]
In this embodiment, a plurality of input data or operation instructions to be input are displayed, such as displaying a pull-down menu on the display screen, excluding character input such as the subject's name and address, and the selection target is displayed. An input operation is performed by pointing the specific target with a mouse. Thus, most operations can be performed by mouse operation, so that a comfortable operation environment can be provided to an operator who is not accustomed to the keyboard, and input / operation time can be reduced. As for the plurality of selection targets in the pull-down menu, selection targets that can be input / operated in this device or its input / operation state are set and displayed in advance, so that erroneous input / erroneous operation can be reduced. Further, in this embodiment, since the cursor can be positioned on the input area and the input can be performed via the keyboard, the degree of freedom of the input by the operator is secured. Further, in this embodiment, it is assumed that characters are input using a keyboard. However, it is also possible to display a dialog of the keyboard at the time of input and operate the mouse with a mouse to input the dialog. Further, a handwriting input dialog may be displayed and handwriting input may be performed by operating the mouse. Further, a touch panel may be provided in the display unit, and input / operation may be performed on a screen via a fingertip or an input pen. Thus, the operability of input / operation can be remarkably improved.
[0052]
FIG. 18 shows the flow of data measurement in step S-12 in FIG. First, a grid map of the magnetocardiogram waveform is displayed as an initial screen as shown in FIG. 25 (S-12-1). In the figure, in the operation area section, operations for channel selection, ON / OFF of the waveform monitor, lock / unlock of the FLL circuit 6, automatic adjustment of the offset voltage of the AFA 7, and heat flash can be respectively performed. Further, it is possible to set signal sampling conditions, set the scale of waveform display, and set AFA parameters.
[0053]
The channels consist of 64 8 × 8 channels, and all channels can be selected by clicking the “Select All Channels” button or by dragging the channel matrix diagonally from end to end. Also, by dragging the channel matrix in units of rows or columns, channels can be selected in units of rows or columns. In any case, the magnetocardiogram waveform of the selected channel is displayed in the analysis data section. In the case of selection in units of rows or columns, the display is as shown in FIG. In this case, the selected waveform is enlarged to the full scale on the time axis and displayed. That is, when all 64 channels are selected, the analysis data portion is divided into upper, lower, left and right (cells) as shown in FIG. 26, as shown in FIG. 26, the analysis data portion is divided into upper and lower parts to form a familiar graph form in which the time axis is left and right, thereby giving a display form giving priority to visibility.
[0054]
As for the monitoring of the waveform, when the "ON" button is pressed, the acquisition of the signal and the updating of the waveform are repeated at a specified time, for example, from 0.5 sec to 2 sec, and the magnetocardiographic signal of the subject is monitored. You. When the “OFF” button is pressed, the updating of the waveform is stopped. For the FLL, by clicking the “Lock” button or the “Unlock” button, the magnetic field lock can be performed on the 64 SQUID sensors, and the lock can be released. In this case, if one button is pressed, the state is maintained until the other button is pressed. This avoids an unselected malfunction state.
[0055]
Clicking the "AFA offset adjustment" button automatically adjusts the offset voltage. When the "heat flash" button is clicked, a heat flash operation dialog box shown in FIG. 14 is opened. When a channel is selected with a mouse or an arrow key and an "OK" button is clicked, a heat flash operation is executed for the SQUID of the selected channel. If the "Cancel" button is pressed, the dialog box closes and the process ends.
[0056]
As for the sampling time (measurement time) and interval, if a corresponding text box with an inverted triangle is clicked, a pull-down menu of selectable numerical values is opened, and a desired numerical value can be selected from the pull-down menu. The selectable numbers are, for example, 1 sec, 5 sec, 10 sec, 30 sec, 1 min and 2 min for the time, and 0.1 msec, 0.5 msec, 1.0 msec, 2.0 msec, 4.0 msec, 5 m for the intervals. 0.0 msec and 10.0 msec. The time may be selected from about 1 sec to about 24 h as needed. The "time" in the "scale" box is the time scale in msec units, that is, the horizontal scale, and the "signal" is the scale of the A / D converted signal, that is, It means the vertical scale. As for these, similarly to the selection of the sampling time and interval, a desired numerical value is selected from a pull-down menu opened by clicking the corresponding text box.
[0057]
The AFA parameters include an input gain Igain, an output gain Ogain, a low-pass filter (LPF) frequency (reference frequency), a notch filter (BEF) frequency, and a high-pass filter frequency (reference frequency). Similarly, a desired number or character is selected from a pull-down menu opened by clicking a corresponding text box. The numbers or characters are for example 1, 2, 5, 10, 20, 50, 100, 200, 500 and 1000 for Igain, for example 1, 10 and 100 for Ogain, for example 30 Hz for LPF, 50 Hz, 80 Hz, 100 Hz, 200 Hz, 400 Hz and 1 kHz, Off for BEF, 50 Hz and 60 Hz, and for HPF e.g. 0.05 Hz, 0.1 Hz and Thru. It is. These may be input from a keyboard instead of selection.
[0058]
As channels, in addition to 64 channels, for example, 16 auxiliary channels may be prepared, and an electrocardiographic waveform may be obtained in the auxiliary channels. The waveform of the tenth channel as a reference channel is displayed at the bottom of FIG. 25, which is the electrocardiographic waveform obtained on the tenth channel among the auxiliary channels. The magnetocardiogram waveform generally contains magnetic noise, while the electrocardiogram waveform does not contain such noise. Therefore, by comparing the displayed magnetocardiogram waveform with the electrocardiogram waveform of the reference channel, information as to whether or not the magnetocardiogram waveform contains magnetic noise can be obtained. Of course, the electrocardiographic waveform may be obtained from any of the regular 64 channels instead of the auxiliary channel. Further, other than the electrocardiographic waveform, an electroencephalogram, a blood flow waveform, a blood pressure waveform, or the like may be used. Further, the electrocardiographic waveform of the pregnant woman and the magnetocardiographic waveform of the fetus may be compared. As the reference waveform, not only a reference waveform of one channel but also a reference waveform of a plurality of channels may be displayed. Further, the reference channel may be used to input not only a signal from a living body but also various control signals for maintenance or the like.
[0059]
Returning to the flowchart of FIG. 18, in step S-12-2, the monitor channel is selected in the manner described above (S-12-2), and when the lock button of the FLL is pressed, the magnetic field lock of all SQUIDs is performed. This is performed (step S-12-3). In this state, the setting of the sampling time and signal, which are the measurement parameters, and the setting of the AFA parameters are performed (step S-12-4). With respect to the setting, this setting can be omitted from the next time by using the set condition. As a result, it is not necessary to set the conditions every time, so that the setting time can be reduced. The setting conditions may be recorded by giving a name or the like so that the setting conditions can be called.
[0060]
When the "start" button in the "measurement" box is pressed, the measurement is started, "measurement in progress" is displayed and a progress bar indicating the progress of the measurement is displayed as shown in FIG. S-12-5). In the progress bar of this embodiment, the bar extends in the bar graph format from left to right as the process progresses. If the entire process content (time) is set to 100, for example, the progress of the process is known. For example, a pie chart or the like may be used. The progress bar displays the measurement screen at a predetermined position on the measurement screen while leaving the measurement screen around the measurement screen. As a result, the contents of the display screen do not change significantly, and erroneous operations can be reduced.
[0061]
When the measurement is started, the displayed signal waveform is fixed as it is, and the progress bar is displayed on the fixed display screen. Updating of the progress bar is repeated, for example, every second until the set time ends (S-12-6). When the "stop" button in the "measurement" box is pressed, the measurement is stopped. When the measurement is completed, the screen shown in FIG. 26 is displayed, and the waveform is confirmed (S-12-7). Thereafter, the necessity of storing the data is determined (S-12-8). If the data needs to be stored, the menu "File (F)"-"Save (S)" is selected, and the signal is output. The data is stored and added to the data list of the subject (S-12-9). Thereafter, it is determined whether measurement is necessary again, including the case where storage is not necessary (S-12-10). If necessary, the above steps are repeated. If not, all steps of data measurement are completed. I do. In this case, the menu is selected so that the screen display is as shown in FIG. FIG. 26 shows an example in which the channel in the second column is selected as the channel.
[0062]
In FIG. 26, there is a scroll bar-262 in which a scroll box 261 moves at the bottom of the analysis data section. The scroll box 261 is movable between the left and right ends of the scroll bar-262, and the width w of the scroll box 261 represents a time scale. The time width between the left and right ends of the scroll bar-262 indicates the measurement time, and therefore, the displayed waveform is a part of the waveform generated during the measurement time, which corresponds to the time scale w of the scroll box 261. 5 is an enlarged waveform of FIG. Thus, the operator indicates how long (the width of the scroll box 261) the waveform currently displayed in the analysis data section is within the measurement time (the width of the scroll bar-262). Since it is possible to visually recognize whether the waveform indicates the first half or the second half of the measurement time, the visibility can be improved. Further, since the position of the scroll box 261 and the waveform of the analysis data section are linked with each other, the display area of the analysis data section is moved for a predetermined time by moving the cursor while dragging the cursor in accordance with the scroll box 261. May be viewed. In this way, the scroll bar-262 can be treated as a table of contents, for example, the waveform of a predetermined time can be easily confirmed, and the measurement contents can be confirmed in detail.
[0063]
FIG. 19 shows the flow of the averaging process in step S-10 in FIG. In the averaging process, in order to remove noise of each data measured in each channel, addition is performed for each channel and an average value is calculated. The following processing is performed to set a reference time for the averaging processing of each channel, and the averaging processing of each channel is executed.
[0064]
First, raw data of the designated subject is read (S-10-1). This is performed by selecting a line in the list data of FIG. 24 in which "data type" is raw data. Thus, the magnetocardiogram waveform of the channel in the first column shown in FIG. 27 is displayed as an initial display (S-10-2). Next, a display channel is selected (S-10-3). FIG. 27 shows an example in which the channel in the second column is selected. Thereafter, the setting channel is designated (S-10-4). This is performed by clicking a text box of “channel” in a box of “averaging condition” in the operation area portion or by clicking a mouse on an area displaying a waveform of a desired channel. In this case, clicking the triangle button in the text box increases the channel number, and clicking the inverted triangle button decreases the channel number. FIG. 27 shows an example in which the specified channel is the channel in the second column and second row. By this channel designation, a channel serving as a reference time for the averaging process is specified. In specifying this channel, a channel having the most typical or easy-to-understand waveform may be selected. If there is no channel having a good waveform in the selected row or column, it is possible to start over from the step (S-10-3) again.
[0065]
When one channel of the analysis data section is specified, a threshold cursor 271 is displayed at the position of the specified channel as shown in FIG. Thereby, the designated channel can be visually confirmed. In FIG. 27, three slider cursors 273 to 275 movable in the left-right direction are displayed on the upper part of the analysis data unit 805. These are displayed at the same time when the display screen of FIG. 27 is displayed. Displayed automatically.
[0066]
In step S-10-5, a threshold value, an offset time, and an averaging time, which are averaging conditions, are set by clicking a corresponding text box in a box of "averaging condition". Click the triangle button to increase the number; click the inverted triangle button to decrease the number. The threshold is set by selecting a number representing the threshold in the "Threshold" text box. As a result, the threshold cursor 271 automatically moves to the position corresponding to the numeral. The moving position in that case can be clearly confirmed by the cursor line. The change (selection) of the number representing the threshold value in the “threshold value” text box and the movement of the threshold cursor 271 are linked with each other, and thus the threshold value moves the threshold cursor 271. This can also be set. The slider cursor 273 indicates the position (reference time) at which the rising portion of the waveform coincides with the threshold value set by the slider cursor 271, and the indicated position is clearly indicated by the cursor line. Can be confirmed. The slider cursor 273 moves following the reference time position so as to coincide with the reference time position that changes according to the change of the threshold value. When a number representing the offset time is selected in the "Offset" text box and a number representing the averaging time is selected in the "Time" text box, the slider cursors 274 and 275 move to the corresponding positions. The moving position in that case can be clearly confirmed by the cursor lines of those slider cursors. Movement of the slider cursors 274 and 275 is linked to the selection of numbers in the "Offset" and "Time" text boxes. For this reason, the setting of the offset time and the averaging time can also be performed by moving the slider cursor.
[0067]
In step S-10-6, whether averaging should be performed automatically or manually as the averaging mode is set (S-10-7). After that, when the "setting" button is clicked (S-10-8), an addition process is performed. Clicking the "Cancel" button cancels all the averaging conditions. As the addition processing, a time t (reference time) exceeding the threshold value is searched for in each channel (S-10-9), and a waveform is displayed around the time t (that is, the time t is displayed on the screen). A waveform from t-50 msec to t + 50 msec is displayed so as to be positioned at the center (S-10-10), and it is determined whether or not cancel is selected in the manual mode (S-10-11). If cancel is not selected, waveform addition is performed (S-10-12). Steps S-10-9 to S-10-12 are repeated by the number of additions for each channel, and the added data is divided by the number of additions (S-10-13). -The zing process ends.
[0068]
In this manner, on the display screen of the averaging process shown in FIG. 27, it is possible to input / operate various conditions of the averaging process from both the operation area portion and the analysis data portion 805. A variety of setting methods can be provided to the operator, such as setting a rough condition in the analysis data section and setting a detailed value in the operation area section, and the time for setting the condition can be reduced. In particular, in this embodiment, by specifying a specific channel from the plurality of channels of the analysis data section 805 via a cursor, the channel can be used as a reference for the averaging process, thereby reducing erroneous operations. Operability is improved. Moreover, by displaying a threshold cursor or a slider cursor near the designated channel, visual recognition is possible. Furthermore, since the input / operation of various conditions of the averaging process can be performed by a threshold cursor or a slider cursor arranged near the channel, the movement of the operator's line of sight can be reduced. Since it can be set visually according to the waveform display, erroneous operations can be reduced and operability can be improved. These setting conditions, for example, the latest setting conditions may be stored and displayed as the next setting conditions, or each setting condition may be stored with a name and stored.
[0069]
FIG. 20 shows a flow of the data analysis in step S-11 in FIG. The data analysis is to display various types of waveforms and diagrams to obtain information necessary for diagnosis. By selecting the menu in FIG. 9, various types of waveforms and diagrams are displayed. Can be selectively displayed. That is, if "single waveform display (W)" of "data analysis (A)" is selected, the single waveform screen shown in FIG. 28 is displayed (S-11-2), and "data analysis (A)" is displayed. )), The superimposed waveform screen (M) shown in FIG. 29 is displayed (S-11-3), and the grid map display (G) of the data analysis (A) is selected. Is selected, the grid map waveform screen shown in FIG. 30 is displayed (S-11-4), and if "Equimagnetic field map (B)" of "Data analysis (A)" is selected, it is shown in FIG. When the “Propagation time diagram (P)” of “Data analysis (A)” is selected in the isomagnetic diagram screen displayed (S-11-5), the propagation time diagram screen shown in FIG. 11-6) Then, if "Time integration diagram (T)" of "Data analysis (A)" is selected, the time integration diagram screen shown in FIG. 33 is displayed (S-11-7). Is, also, the system is completed by selecting the "File (F)" and "End of magnetocardiogram system (X)" of.
[0070]
In each screen, if a radio button (circular button in the drawing) in the operation area is clicked, the screen of the waveform or diagram specified by the click is displayed instead. For this reason, in FIG. 20, the branch portion is not "branched by menu" but "branched by menu or radio button". Therefore, according to this embodiment, various analysis data can be obtained only by clicking the radio button in the operation area without selecting the menu in FIG. 9, so that the operation time can be shortened. In addition, erroneous operation can be reduced and operability can be improved.
[0071]
In FIGS. 28 to 30, the magnetic flux density in the "scale" box is the value of the full scale on the plus side and the minus side with reference to the zero level (the unit is picotesla (pT)). The value is selected in a pull-down menu that is opened by clicking the triangle button in that text box. 28 to 33, by clicking a radio box in a “display component” box, a waveform of a normal component or a waveform of a tangent component can be selected and displayed on a screen.
[0072]
In FIG. 28, the waveform of each channel selected by the channel is displayed with the left end of the analysis data portion adjusted to the offset time. According to this analysis data, it is possible to compare the shapes and sizes of the waveforms of the channels arranged and displayed vertically in the analysis data section. Similarly, in FIG. 29, the waveforms arranged vertically in FIG. 28 are displayed in a superimposed manner, and the shapes and sizes of the waveforms can be compared. In FIG. 30, all the channels are displayed on the basis of the offset time as in FIGS. 28 and 29. Therefore, the operator can select and analyze the number of channels as needed.
[0073]
In FIG. 31, an elongated magnetic field strength index box 310 is vertically arranged at the right end of the analysis data section. The magnetic field strength index box is divided into twelve sections having different colors. This is to improve the visual (color) recognizability by distinguishing the intensity range of the magnetic field indicated by each island pattern on the isomagnetic diagram screen shown in FIG. 31 by the type of color. It is. That is, the center position 311 in the longitudinal direction of the magnetic field strength index box 310 is a position where the magnetic field strength is zero, and the sections above the center position are referred to as the first to sixth sections in order from the center position. Then, for example, the first section is in a magnetic field intensity range of 0 to 2 pT, the second section is in a magnetic field intensity range of 2 to 4 pT, the third section is in a magnetic field intensity range of 4 to 6 pT, and the fourth section is 6 to 8 pT. The fifth section corresponds to a magnetic field strength range of 8 to 10 pT, and the sixth section corresponds to a magnetic field strength range of 10 to 12 pT. The same applies to the section below the center position. However, a section above the center position indicates a plus direction magnetic field strength, and a section below the center position indicates a minus direction magnetic field strength. The isomagnetic diagram shown in FIG. 31 is displayed in different colors according to the magnetic field strength in accordance with the definition of the correspondence between the magnetic field strength range and the color in the magnetic field strength index box 310. In addition, as a color, the plus side of the magnetic field strength may be a warm color system, the minus side may be a cool color system, and the center may be yellow. Thereby, the strength of the magnetic field can be recognized in color, so that the visibility can be improved. In addition, according to this embodiment, the magnetic field strength index box 310 is provided near the analysis data section, so that the color to be compared, that is, the color added to the map and the predetermined color of the magnetic field strength index box 310 Can be checked while comparing the movement of the line of sight without greatly moving the line of sight, so that the relationship between the level of the magnetic field and the color can be clearly determined. In this embodiment, the magnetic field strength index box 310 is provided at the right end of the analysis data section. However, the magnetic field strength index box 310 may be provided near the analysis data section.
[0074]
In FIG. 31, “number of maps” in the “reconstruction parameter” box indicates the number of displayed isomagnetic maps, and “maximum value” indicates both ends of the magnetic field strength index box 310. The term “interval” corresponding to the magnetic field strength means a magnetic field range corresponding to the length of each section in the magnetic field strength index box 310. The value can be selected by clicking the triangle or inverted triangle button in the corresponding text box.
[0075]
The ECG waveform of the reference channel and two map time selection cursors 311 and 312 are displayed at the bottom of the analysis data section. A division line having the same interval is displayed between the two map time selection cursors 311 and 312, and the number of the division lines matches the number of maps selected by the map number selection. In addition, the two map time selection cursors 311 and 312 can move their positions independently in the left-right direction. If the movement changes the distance between the cursors, the space between the dividing lines also changes. The intervals between the dividing lines are always equal. Of course, each division line may be provided with a cursor and individually set one by one. In FIG. 31, the number of displayed isomagnetic maps is 16, but these diagrams are diagrams at the time when the dividing line on the electrocardiographic waveform is located. The time is also displayed so you can see when the map is up to date.
[0076]
Thus, as described with reference to FIG. 26, the operator can determine how much the map (currently displayed in the analysis data section) is within the analysis time (the width of the electrocardiographic waveform) (two car (The width of the sols 311 and 312), and it is possible to grasp at a glance where the range indicated by the map is within the analysis time, so that the visibility can be improved. Further, the range indicated by the map can be set by simply moving the two cursors with a mouse operation, so that the operation is easy. Furthermore, if the interval between the dividing lines is freely set, various analysis environments can be provided to the operator, for example, by making the questionable part dense and the other part sparse.
[0077]
Also, when a check mark is displayed by clicking the check box of “current direction” in the operation area, an arrow is displayed on the isomagnetic diagram. The isomagnetic map on which the arrow is displayed is called an allo map (not shown). For the arrows, the position indicates the position of the channel (position of the magnetic sensor), the length indicates the strength of the magnetic field, and the direction indicates the direction of the current when the direction of the magnetic field is converted into the direction of the current.
[0078]
In FIG. 32, the starting point of the propagation time (t in FIG. 7)1The position can be changed by changing the position of the cursor 321 moving on the reference waveform, and the position of the cursor 321 can be changed by dragging the cursor with the mouse. .
[0079]
FIG. 33 shows two time integration diagrams and one difference diagram. "Difference display" in the propagation time diagram can be easily displayed by clicking the check box and displaying a check mark. When the "differential display" is checked, four cursors 331 to 334 appear on the reference waveform, and two time integration diagrams are displayed on the upper left and right, and a difference diagram is displayed on the lower left as shown. . The two time integral diagrams integrate the magnetocardiogram over the time range of 100 msec to 140 msec and 180 msec to 240 msec set using the cursors 331 and 332 and the cursors 333 and 334 on the reference waveform, respectively. Based on the values, the respective time ranges can be changed by dragging the cursors 331 and 332 and the cursors 333 and 334, respectively, with a mouse. The difference diagram represents the difference between the two time integral diagrams. If there is no check mark, only two cursors (for example, the cursors 331 and 332) appear for the cursor, and only one time integration diagram is displayed for the time integration diagram. Of course, the integration time can be changed by changing the position of the cursor. As described above, according to this embodiment, by clicking the check mark for “differential display”, two sets of cursors for prompting the next operation are displayed, so that the operation time can be reduced without giving any doubt about the operation. The time range can be easily set by moving the two sets of cursors with a mouse, so that operability can be improved.
[0080]
FIG. 21 shows the flow of system adjustment. The flow is branched into five by menu selection (S-15). In this case, the currently displayed data is stored as it is (S-13), and the Φ-V characteristic curve shown in FIG. 34 is displayed (S-14). When "Fully automatic adjustment (A)" of "Data measurement (Q)" is selected, IbiasAnd VOFFIs automatically calculated (S-16) (the automatic calculation will be described later), the calculated adjustment value is set in the FLL circuit (S-17), and the flow returns to step S-14.
[0081]
"V" of "Data measurement (Q)"OFFIf “Adjust (V)” is selected, the VoFFIs calculated (S-18) (VOFFWill be described later), and then the flow proceeds to step S-17 described above. When "manual adjustment (M)" of "data measurement (Q)" is selected, a manual adjustment dialog box shown in FIG. 12 is opened (S-19). The operator is IbiasAnd VOFFIs input for each channel, the input adjustment value is accepted (S-20), and the adjustment value is set in the FLL circuit by pressing an "OK" button in the dialog box (S-). 21). The dialog box is thereby closed (S-22), and the flow proceeds to step S-17.
[0082]
When the "adjustment value file (F)" of the "data measurement (Q)" is selected, the flow is further branched into two by the sub pull-down menu (S-22). That is, a file is inquired to the operator by selecting "Open (O)" in the sub pull-down menu, and the operator inputs a file including the contents of FIG. 11 (S-23). The contents of the adjustment value file are set in the FLL circuit (S-24). Also, by selecting "data measurement (Q)"-"adjustment value file (F)"-"overwrite save" or "save as (A)", the same as step S-23 can be performed. The adjustment is performed (S-25), and the adjustment value can be written to the adjustment value file (S-26).
[0083]
FIG. 22 shows the flow of the automatic calculation in step S-16 in FIG.
[0084]
S-15-1: In the case of fully automatic adjustment, the following processing is performed for all SQUID channels one by one. The channel being processed is assumed to be ch.
[0085]
S-15-2: It is assumed that the calculated bias current and offset voltage are stored in memories named IBIAS and VOFF. IBIAS and VOFF can hold values as many as the number of channels, and the initial values are all set to 0. The temporarily stored amplitude of the Φ-V characteristic curve is ΔV.
[0086]
S-15-3: Bias current I for each channel chbFrom 0 to I specified in the "scanning parameter" box in FIG.biasUp to ΔIbiasIs changed (scanned) in the step of I, and the amplitude ΔV of the Φ-V characteristic curve becomes the largest IbiasIs the optimum bias current IBIAS (ch) for the channel ch.
[0087]
S-15-4 to 8: Between the bias current 0 and Ib, the amplitude ΔV of the Φ-V characteristic curve is obtained by the following processing. The external magnetic field Φ given to the SQUID of the channel ch is changed from 0 to Φ specified in the “scanning parameter” box in FIG.extThen, the A / D converted signal is stored and the maximum value Mmax and the minimum value Vmin of the signal are obtained.
[0088]
S-15-9-10: Here, maximum value MmaxAnd the minimum value VminIs larger than ΔV, the value of IBIAS (ch) is changed to IbiasAnd the value of VOFF (ch) is set to the maximum value MmaxAnd the minimum value VminΔV is the maximum value VmaxAnd the minimum value VminReplace with the difference of If ΔV is the maximum value MmaxAnd the minimum value VminIf it is larger than the difference, the previous value is retained.
[0089]
S-15-11: The above processing is performed with the bias current IbIs IbiasIBIAS (ch) and VOFF (ch) when repeated until becomes the optimum bias current and offset voltage of the SQUID channel ch.
[0090]
S-15-12: The fully automatic adjustment is completed by executing the above processing for all SQUID channels.
[0091]
FIG. 23 shows the VOFF adjustment flow in step 18 of FIG.
[0092]
S-18-1: VOFFIn the case of adjustment, the following processing is performed for all SQUID channels one by one.
[0093]
S-18-2 to 6: The external magnetic field Φ given to the SQUID of the channel ch is changed from 0 to Φ specified in the “scanning parameter” box in FIG.ext(Scanning) in steps of ΔΦ until the maximum value V of the A / D converted signalmaxAnd the minimum value VminFind the difference between
[0094]
S-18-7: The optimum offset voltage VOFF (ch) of the SQUID channel ch is the maximum value VmaxAnd the minimum value VminIs calculated as the difference between
[0095]
S-18-8: By performing the above processing on all SQUID channels, VOFFThe adjustment ends.
[0096]
【The invention's effect】
According to the present invention,It is possible to provide a biological magnetic field map display device with improved operation and visibility.
[Brief description of the drawings]
FIG. 1 is a schematic configuration diagram of an embodiment of a biomagnetic field measurement apparatus according to the present invention.
FIG. 2 is a perspective view showing an arrangement configuration of a magnetic sensor used in the biomagnetic field measuring apparatus of FIG. 1;
FIG. 3 is a perspective view of a single magnetic sensor for detecting a normal component of a magnetic field, which is used in the biomagnetic field measuring apparatus of FIG. 1;
FIG. 4 is a perspective view of a single magnetic sensor for detecting a normal component of a magnetic field used in the biomagnetic field measuring apparatus of FIG. 1;
FIG. 5 is a diagram showing a positional relationship between a magnetic sensor and a chest of a subject in the biomagnetic field measuring apparatus of FIG. 1;
6 is a time waveform diagram of each component of a healthy person's biomagnetic field (cardiac magnet) measured by each magnetic sensor in the biomagnetic field measuring apparatus of FIG. 1;
FIG. 7 is a view showing a time waveform of a tangential component of a specific two-channel magnetocardiogram measured for a healthy person in the biomagnetic field measuring apparatus of FIG. 1;
FIG. 8 is a diagram showing a basic layout of a display screen displayed on a display unit in the biomagnetic field measuring apparatus of FIG. 1;
9 is a diagram showing an operation menu in a menu section of a display screen displayed on a display unit in the biomagnetic field measuring apparatus of FIG. 1;
10 is a test object that is opened when “list (L)”-“registration (R)” is selected as an operation menu in a display screen displayed on the display unit in the biomagnetic field measurement apparatus of FIG. The figure which shows the content of the person registration dialog box.
11 is a search dialog that is opened when “list (L)”-“search (S)” is selected as an operation menu in a display screen displayed on the display unit in the biomagnetic field measurement apparatus of FIG. The figure which shows the content of a box.
FIG. 12 shows a case where “data measurement (Q)”-“manual adjustment (M)” is selected as an operation menu in a display screen displayed on the display unit in the biomagnetic field measurement apparatus of FIG. The figure which shows the content of the manual adjustment dialog box opened.
FIG. 13 shows a case where “data measurement (Q)” − “measurement panel (P)” is selected as an operation menu in a display screen displayed on the display unit in the biomagnetic field measurement apparatus of FIG. The figure which shows the content of the automatic diagnosis dialog box opened.
14 is a heat flash that is opened when a “heat flash” button is pressed when the display screen of FIG. 25 or 26 is displayed on the display unit in the biomagnetic field measurement apparatus of FIG. 1; The figure which shows the content of the operation dialog box.
FIG. 15 is a diagram showing a measurement progress bar displayed on a display unit during measurement in the biomagnetic field measurement device of FIG. 1;
FIG. 16 is a diagram showing a flow of an overall operation performed in the biomagnetic field measuring apparatus of FIG. 1;
FIG. 17 is a diagram showing a flow of subject selection in a subject selection step in the operation flow of FIG. 16;
FIG. 18 is a diagram showing a flow of data measurement in a data measurement step in the operation flow of FIG. 16;
FIG. 19 is a diagram showing a flow of an averaging process in an averaging process step in the operation flow of FIG. 16;
FIG. 20 is a diagram showing a flow of data analysis in a data analysis step in the operation flow of FIG. 16;
FIG. 21 is a diagram showing a flow of system adjustment performed in the biomagnetic field measuring apparatus of FIG. 1;
FIG. 22 is a diagram showing an automatic calculation flow in an automatic calculation step in the system adjustment flow of FIG. 21.
FIG. 23 shows V during the system adjustment flow of FIG. 21;OFFV in adjustment stepOFFThe figure which shows an adjustment flow.
FIG. 24 is a view showing a subject list screen displayed on a display unit of the biomagnetic field measuring apparatus in FIG. 1;
FIG. 25 is a diagram showing a data measurement screen displayed on a display unit of the biomagnetic field measurement apparatus of FIG. 1;
FIG. 26 is a diagram showing a waveform confirmation screen displayed on the display unit of the biomagnetic field measurement device of FIG. 1 when measurement is completed.
FIG. 27 is a view showing an averaging screen displayed on a display unit of the biomagnetic field measuring apparatus of FIG. 1;
FIG. 28 is a view showing a single waveform display screen displayed on a display unit of the biomagnetic field measuring apparatus of FIG. 1;
FIG. 29 is a diagram showing a superimposed waveform display screen displayed on the display unit of the biomagnetic field measuring apparatus in FIG. 1;
FIG. 30 is a view showing a grid map display screen displayed on a display unit of the biomagnetic field measuring apparatus of FIG. 1;
FIG. 31 is a view showing an isomagnetic map display screen displayed on the display unit of the biomagnetic field measuring apparatus in FIG. 1;
FIG. 32 is a diagram showing a propagation time diagram display screen displayed on the display unit of the biomagnetic field measurement apparatus of FIG. 1;
FIG. 33 is a diagram showing a time integration diagram display screen displayed on a display unit of the biomagnetic field measurement apparatus in FIG. 1;
FIG. 34 is a view showing a system adjustment screen displayed on a display unit of the biomagnetic field measuring apparatus of FIG. 1;
[Explanation of symbols]
1: magnetic shield room, 2: subject, 3: bed, 4: dewar, 5: automatic supply device, 6: FLL circuit, 7: amplifier / filter / amplifier, 8: computer, 8-1 : Display section, 8-2: keyboard, 8-3: mouse, 20-1 to 20-8, 21-1 to 21-8, 22-1 to 22-8, 23-1 to 23- 8, 241-1 to 24-8, 25-1 to 25-8, 26-1 to 26-8 and 27-1 to 27-8: magnetic sensors, 10, 10 'and 10 "and 11, 11' and 11 ": coil, 12, 12 'and 12": SQUID, 13 and 14: sensor, 30: chest, 261: scroll box, 262: scroll bar, 271: threshold cursor, 273-275 : Slider cursor, 311, 312, 321 and 331 to 334: cursor 801: title bar section, 802: menu section, 803: tool bar section, 804: subject information section, 805: analysis data section, 806: operation area section, 808: status bar Part, 807-1: message bar part, 807-2: date and time display part, 808 to 814: icon.

Claims (8)

被検者の生体磁場を複数の位置で計測してその処理結果としての磁場図を表示する生体の磁場図表示装置において、
入力手段と表示部とを有する計算機を備え、
前記計算機は、
計測された生体磁場を複数のチャンネルの行または列として前記表示部の解析データ部に波形表示し、
前記被検者の他の生体信号による参照チャンネルを前記解析データ部に表示された生体磁場の複数のチャンネルの行または列の下部に表示し、
前記入力手段からのデータ解析の指示を受け付けて、前記解析データ部に当該指示に該当する解析磁場を表示すると共にその下部に前記参照チャンネルを表示し、
前記解析磁場の下部に表示される参照チャンネル上に、表示される解析磁場の時刻に対応する分割線および前記分割線の位置を移動するための時刻選択用カーソルを表示し、
前記入力手段による前記カーソルの移動を受けて、移動後の選択時刻に対応する解析磁場を表示する
ことを特徴とする生体の磁場図表示装置。In a biological magnetic field map display device that measures the biomagnetic field of the subject at a plurality of positions and displays a magnetic field map as the processing result,
A computer having input means and a display unit,
The calculator is:
The waveform of the measured biomagnetic field is displayed in the analysis data section of the display section as rows or columns of a plurality of channels,
A reference channel based on another biological signal of the subject is displayed below a row or column of a plurality of channels of the biomagnetic field displayed in the analysis data unit,
Receiving a data analysis instruction from the input means, displaying the analysis magnetic field corresponding to the instruction in the analysis data section and displaying the reference channel below the analysis magnetic field,
On the reference channel to be displayed at the bottom of the analysis field, displays the time selection cursor for moving the position of the dividing line and the dividing line corresponding to the time of the analysis field to be displayed,
An apparatus for displaying a magnetic field map of a living body, wherein an analysis magnetic field corresponding to a selected time after the movement is displayed in response to movement of the cursor by the input means. 請求項１において、前記被検者の他の生体信号は被検者の心電波形、脳波、血流波形または血圧波形であることを特徴とする生体の磁場図表示装置。2. The apparatus according to claim 1, wherein the other biological signal of the subject is an electrocardiographic waveform, an electroencephalogram, a blood flow waveform, or a blood pressure waveform of the subject. 請求項１または２において、
前記解析データ部に表示される解析磁場は、等磁線図、伝播時間図または時間積分図である
ことを特徴とする生体の磁場図表示装置。In claim 1 or 2,
The analysis magnetic field displayed in the analysis data section is a magnetic field diagram, a propagation time diagram, or a time integration diagram, wherein the magnetic field map display of the living body is provided. 請求項１または２において、
前記解析データ部に表示される解析磁場は等磁線図であり、前記時刻選択用カーソルは２箇所表示され、前記分割線は前記２つのカーソル間に等間隔に表示されることを特徴とする生体の磁場図表示装置。In claim 1 or 2,
The analysis magnetic field displayed in the analysis data section is an isomagnetic diagram, the time selection cursor is displayed at two places, and the dividing line is displayed at equal intervals between the two cursors. A biological magnetic field diagram display device. 請求項１または２において、
前記解析データ部に表示される解析磁場は等磁線図であり、前記時刻選択用カーソルは各分割線に設定され、前記分割線は個別に移動可能であることを特徴とする生体の磁場図表示装置。In claim 1 or 2,
The analysis magnetic field displayed in the analysis data section is an isomagnetic map, the time selection cursor is set on each division line, and the division lines are individually movable, and the magnetic field diagram of the living body is characterized in that: Display device. 請求項１または２において、
前記解析データ部に表示される解析磁場は伝播時間図であり、前記時刻選択用カーソルの移動により伝播時間の起点位置を変更可能であることを特徴とする生体の磁場図表示装置。In claim 1 or 2,
The analysis magnetic field displayed in the analysis data section is a propagation time diagram, and the starting position of the propagation time can be changed by moving the time selection cursor. 請求項１または２において、前記解析データに表示される解析磁場は時間積分図であり、
前記計算機は、
前記入力手段からの差分表示の指示を受け付けて、前記参照チャンネル上に４つの時刻選択用カーソルを表示し、
前記４つのカーソル位置に対応する時間範囲で時間積分した２つの時間積分図およびその差分図を前記解析データ部に表示する
ことを特徴とする生体の磁場図表示装置。The analysis magnetic field displayed in the analysis data according to claim 1 or 2, which is a time integration diagram,
The calculator is:
Receiving a difference display instruction from the input unit, displaying four time selection cursors on the reference channel,
An apparatus for displaying a magnetic field map of a living body, wherein two time integration diagrams obtained by time integration in a time range corresponding to the four cursor positions and a difference diagram thereof are displayed in the analysis data section. 請求項１乃至７のいずれかにおいて、
前記参照チャンネルは、複数のチャンネルの中から所定のチャンネルを選択可能に表示されることを特徴とする生体の磁場図表示装置。In any one of claims 1 to 7,
The magnetic field map display of a living body, wherein the reference channel is displayed so that a predetermined channel can be selected from a plurality of channels.
JP27015897A 1997-10-02 1997-10-02 Measuring device for biological magnetic field Expired - Fee Related JP3595660B2 (en)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP27015897A JP3595660B2 (en) 1997-10-02 1997-10-02 Measuring device for biological magnetic field
US09/164,271 US6473518B1 (en) 1997-10-02 1998-10-01 Method of measuring a biomagnetic field, method of analyzing a measured biomagnetic field, method of displaying biomagnetic field data, and apparatus therefor
US10/282,158 US6711281B2 (en) 1997-10-02 2002-10-29 Method of measuring a biomagnetic field, method of analyzing a measured biomagnetic field, and method of displaying biomagnetic field data, and apparatus therefor
US10/770,525 US7194121B2 (en) 1997-10-02 2004-02-04 Method of measuring a biomagnetic field, method of analyzing a measured biomagnetic field, method of displaying biomagnetic field data, and an apparatus therefor
US11/717,667 US7433506B2 (en) 1997-10-02 2007-03-14 Method of measuring a biomagnetic field, method of analyzing a measured biomagnetic field, method of displaying biomagnetic field data, and an apparatus therefor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP27015897A JP3595660B2 (en) 1997-10-02 1997-10-02 Measuring device for biological magnetic field

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2002251109A Division JP3768178B2 (en) 2002-08-29 2002-08-29 Magnetic field diagram display method

Publications (2)

Publication Number Publication Date
JPH11104095A JPH11104095A (en) 1999-04-20
JP3595660B2 true JP3595660B2 (en) 2004-12-02

Family

ID=17482356

Family Applications (1)

Application Number Title Priority Date Filing Date
JP27015897A Expired - Fee Related JP3595660B2 (en) 1997-10-02 1997-10-02 Measuring device for biological magnetic field

Country Status (1)

Country Link
JP (1) JP3595660B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6473518B1 (en) 1997-10-02 2002-10-29 Hitachi, Ltd. Method of measuring a biomagnetic field, method of analyzing a measured biomagnetic field, method of displaying biomagnetic field data, and apparatus therefor
JP2001314388A (en) * 2000-05-10 2001-11-13 Toshiba Corp Biomagnetometer

Also Published As

Publication number Publication date
JPH11104095A (en) 1999-04-20

Similar Documents

Publication Publication Date Title
US6961605B2 (en) Method for processing biomagnetic field data, magnetic field contour mapping, forming their waveforms and a biomagnetic instrument using the same
JP3997224B2 (en) Display device of magnetic field measuring device
JPH11276451A (en) Display method and printing method of living body magnetic field measuring data
JP3595660B2 (en) Measuring device for biological magnetic field
JP5085155B2 (en) Biomagnetic field measurement device
JP3228702B2 (en) Measurement method of magnetic field of living body and magnetic field map display
JP3228703B2 (en) Method and apparatus for processing measurement data of magnetic field generated by living body
JP3720960B2 (en) Processing method of measurement data of magnetic field generated by living body
JP4054011B2 (en) Magnetic field measuring device
JP3514655B2 (en) Biomagnetic field measurement device
JP3992705B2 (en) Display device of magnetic field measuring device
JP3058619B2 (en) Biological magnetic field map generator
JP3722140B2 (en) Biomagnetic field measurement device
JP3140734B2 (en) Biological magnetic field map display
JP3140733B2 (en) Biological magnetic field map display
JP3768178B2 (en) Magnetic field diagram display method
JP4044529B2 (en) Display device of magnetic field measuring device
JP2004089724A (en) Organism magnetic field measuring apparatus
JP3714346B2 (en) Biomagnetic field measurement device
JP3140731B2 (en) Biological magnetic field map display
JP3518495B2 (en) Biomagnetic field measurement device
JP3140735B2 (en) How to display the magnetic field map of a living body
JP3140732B2 (en) Biological magnetic field map display
JP2004041788A (en) Biomagnetic field measuring instrument
JP3518494B2 (en) Biomagnetic field measurement device

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20040218

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20040225

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040426

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20040526

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040723

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20040818

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20040906

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20070910

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20080910

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20080910

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20090910

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20090910

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20100910

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20100910

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20110910

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20120910

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20120910

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20130910

Year of fee payment: 9

LAPS Cancellation because of no payment of annual fees