生体磁気計測装置及び測定方法
【課題】ＳＱＵＩＤ磁束計から得られた胎児の心臓や脳から発生する磁気信号と母体から発生する磁気信号の分離といった２つ以上の異なる部位から発生する磁気信号を分離する機能を有する生体磁気計測装置を提供する。
【解決手段】ＳＱＵＩＤ磁束計からの磁気波形や該磁束計以外の入力信号波形のいずれかのチャンネルの波形をトリガー信号として、前記トリガー信号を用いて、該SQUID磁束計で計測された磁気信号波形を加算平均し、かつ基線補正処理を施した後の波形を、前記トリガー信号を用いて、前記計測された磁気信号波形から差し引く手段と、前記差し引く手段から得られた波形を表示する手段を有する生体磁気計測装置。
【解決手段】ＳＱＵＩＤ磁束計からの磁気波形や該磁束計以外の入力信号波形のいずれかのチャンネルの波形をトリガー信号として、前記トリガー信号を用いて、該SQUID磁束計で計測された磁気信号波形を加算平均し、かつ基線補正処理を施した後の波形を、前記トリガー信号を用いて、前記計測された磁気信号波形から差し引く手段と、前記差し引く手段から得られた波形を表示する手段を有する生体磁気計測装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】本発明は、成人，小児，胎児などの心臓や脳などから発生する微弱な磁気信号を計測するＳＱＵＩＤ(Superconducting Quantum Interference Device：超伝導量子干渉素子）磁束計を含めた生体磁気計測装置に関するものであり、特に臨床検査分野において、生体磁気計測装置から得られた生体磁気信号の情報から有用な情報を得る解析手段および表示手段を有する生体磁気計測装置に関するものである。
【０００２】
【従来の技術】従来技術として、成人，小児，胎児などの心臓や脳などから発生する生体磁気信号の検出はＳＱＵＩＤ(Superconducting Quantum Interference Device：超伝導量子干渉素子）磁束計を用いて検出されている。
【０００３】特開平１０−３０５０１９号公報には、生体から発する生体磁気信号を検出する複数の磁束計と、信号の演算処理を行う演算処理手段と演算結果を表示する手段とを有し、磁気シールドルーム内で生体磁場分布を計測する装置について記載されている。
【０００４】特開平１０−３０５０１９号公報には、生体磁気計測装置において、検出された生体磁気信号の生データからノイズを除去する目的で加算平均（アベレージング）処理を有し、処理条件を視覚的に簡便化した表示手段を有する生体磁気計測装置について記載されている。
【０００５】また、特開平５−２４９２０８号公報には、生体磁気信号を加算平均したものを測定信号から引き算し、雑音のみの信号を得る。雑音成分は主に冷凍機から発生する周期的ノイズであることから、該雑音信号を数サイクル分加算平均し、テンプレートを作成する。その後、測定信号から雑音信号のテンプレートを引き算することにより、生体磁気信号のみのデータを得るという方法が開示されている。
【０００６】
【発明が解決しようとする課題】検出された生体磁気信号は、目的とする部位以外の生体磁気信号が混ざっていることがあり、目的とする部位の磁気信号との区別が困難な場合がある。例えば、胎児の心臓から発生する磁気信号と母体から発生する磁気信号が、ほぼ同じ大きさで検出された場合などは、どちらから発生した磁気信号なのかを区別することは、従来の技術では困難であった。
【０００７】本発明では、ＳＱＵＩＤ磁束計から得られた胎児の心臓や脳から発生する磁気信号と母体から発生する磁気信号の分離や、脳から発生する磁気信号と心臓から発生する磁気信号の分離というような２つ以上の異なる部位から発生する信号の分離を目的とした信号処理の手段を有する生体磁気計測装置が提供でき、心臓疾患や脳機能障害等の診断に有用な情報を提供することが可能となる。
【０００８】
【課題を解決するための手段】上述の課題を解決するために、本発明によれば、被験者を支持するための支持体と、被験者から発生する磁場を検出する超伝導量子干渉素子（ＳＱＵＩＤ）を備える磁気センサと、該磁気センサを外部磁場から遮断するシールドと、該超伝導量子干渉素子を低温に保持するためのクライオスタットと、該クライオスタットを保持する手段と、前記超伝導量子干渉素子の駆動回路を含む磁束計回路とを備えた生体磁気計測装置であって、更に、被験者の測定対象である磁気発生源を、磁気以外の物理量の変化で測定する測定手段を備え、該測定手段からの信号と前記磁気センサからの信号を演算し、前記被験者の複数の磁気発生源からの信号を分離する演算装置を備えた生体磁気計測装置が提供される。
【０００９】上記生体磁気計測装置において、前記磁気以外の物理量変化の測定手段が被験者の皮膚の電位の変動を測定する手段であっても良い。
【００１０】また、本発明では、生体から発生する磁場を計測する超伝導量子干渉素子（Superconducting Quantum Interference Device： 以後ＳＱＵＩＤと略す）を用いた１個ないし複数個の磁束計と、該ＳＱＵＩＤ磁束計によって得られた磁気信号の波形を表示する手段と、該ＳＱＵＩＤ磁束計以外の信号、例えば心電計や脳波計などの信号を取り込み表示する手段と、該磁束計によって得られた磁気信号波形、あるいは前記心電計や脳波計などから取り込んだ該磁束計以外の外部入力信号波形を収集データとして、前記収集データのいずれかのチャンネルの信号波形から基準時刻（トリガー信号）を決定し、前記トリガー信号を用いて前記磁気信号波形を加算平均する手段と、前記加算平均手段により得られた磁気信号波形の任意の２点の時刻における磁気信号強度の値を結んだ直線（基線）の傾きが０になるように補正し、かつ前記任意の２点の時刻における磁気信号強度の値が０（オフセット）になるように補正する手段（基線補正手段）を有し、前記加算平均手段と前記基線補正手段によって得られた波形（波形Ｂとする）を前記トリガー信号を用いて、前記計測された磁気信号波形（波形Ａとする）から差し引く手段と、前記差し引く手段によって得られた波形（波形Ｃとする）を表示する手段を備えることができる。
【００１１】また、前記差し引く手段によって得られた波形（波形Ｃ）を、前記トリガー信号を用いて前記計測された磁気信号波形（波形Ａ）から差し引く手段（該手段によって得られる波形を波形Ｄとする）を有し、前記計測された磁気信号波形Ａと、前記波形Ｃと、前記波形Ｄと、前記トリガー信号波形のいずれか２つ以上を時系列に上下並べて同一画面上に表示できる手段を有したものである。
【００１２】更に、前記波形Ｃ、又は前記波形Ｄを前記収集データとして置き換え、上記加算平均手段や基線補正手段、前記差し引く手段などの演算処理を繰り返し行う手段を有し、前記手段によって得られた波形を表示する手段を備えることができる。
【００１３】
【発明の実施の形態】本発明の一実施例を図を用いて説明する。
【００１４】図１は、本発明の実施例である生体磁気計測装置の構成を示す図である。図１に示すように磁気シールドルーム１内には、被験者９が横になるベッド４と複数個（複数チャンネル）のＳＱＵＩＤ磁気センサおよびＳＱＵＩＤ磁気センサを超伝導状態に保持するための冷媒（液体ヘリウムまたは液体窒素）が貯蔵されたクライオスタット２と、クライオスタット２を機械的に保持するガントリー３が配置されている。ベッド４は、Ｘ方向，Ｙ方向，Ｚ方向に移動可能である。磁気シールドルーム１の外部にはＳＱＵＩＤ磁束計の駆動回路５と、アンプ回路およびフィルタ回路ユニット６と、データ取り込みおよびデータ解析用コンピュータ７と、心電計などの外部入力信号を取るための回路８とが配置されている。
【００１５】ＳＱＵＩＤ磁気センサによって検出された生体磁気信号は、アンプ回路およびフィルタ回路ユニット６により増幅され、かつ設定周波数より低い周波数信号を通過させるローパスフィルタや設定周波数より高い周波数信号を通過させるハイパスフィルタ，商用電源周波数だけをカットするノッチフィルタなどの信号処理を経た後、パソコン７に生データとして取り込まれる。また、心電計などの外部入力信号をとるための回路８からの信号波形はパソコン７に生データとして取り込まれる。また、パソコン７は、取り込んだ生データを生データファイルに格納し、波形を画面表示したり、また波形の信号処理や等磁場線図処理などを行い、表示することもできる。
【００１６】図２は、本発明の実施例であるデータ処理の概念フロー図を示したものである。図２において、図１R>１のパソコン７で取り込んだ生体磁気計測信号２１および外部入力信号２２は、生データファイル２３に格納される。また生体磁気計測信号２１および外部入力信号２２は、生データファイル２３に格納されると同時に、ＤＳＰ(Digital Signal Processor）(図示せず）を用いて、リアルタイムで加算平均処理２５などの信号処理を行うこともできる。生データファイル２３からデータ処理するファイルを呼び出し、呼び出したファイルの信号波形データのうち、磁気信号波形あるいは、該磁束計以外の入力信号波形のいずれかのチャンネルの波形をトリガー信号として、該トリガー信号波形を用いて生体磁気計測信号２１（波形Ａ）に対して加算平均処理２５を行う。次に、加算平均処理２５により得られた波形に対して、基線補正処理２６を行う。次に、基線補正処理２６により得られた波形（波形Ｂ）を、前記トリガー信号を用いて、波形Ａとの差分処理２７を行う。波形差分処理２７により得られた波形Ｃは、信号処理波形データファイル２４に格納される。また波形Ｃは、波形Ａとの差分処理２８により、波形Ｄが得られる。波形Ｄは、信号処理波形データファイル２４に格納される。前記信号処理波形データファイル２４のデータは、再度、加算平均処理２５，基線補正処理２６，差分処理２７，差分処理２８の一連のデータ処理を繰り返し行うこともできる。また、トリガー信号として前記信号処理波形データファイル２４のデータを用いることもできる。生データファイル２３，信号処理波形データファイル２４に格納された波形データは、表示区間や表示スケール，表示する波形のチャンネルなどの入力設定に従って、波形表示処理３０および等磁場線図処理３１が行われ、表示画面部３２に表示される。
【００１７】図３は、図２の加算平均処理２５と基線補正処理２６と波形差分処理２７の信号処理方法の概念を示したものである。一実施例として、胎児の心臓磁気信号計測を目的とした場合について説明する。母親の腹部から計測された磁気信号波形３６について、母親の心電計波形３５をトリガー信号として、トリガー信号のピーク値の時刻のＴ１秒前からＴ２秒間を加算区間として、加算平均処理を行う。一般に母親の心電計波形には、胎児心電計波形はほとんど検出できないか、あるいは極めて小さい。したがって、心電計の波形に閾値を設定して、波形のピーク値を検出すれば、母親のみの心電計波形のピーク値を検出できる。前記ピーク値をトリガーとして、加算区間Ｔ２で加算平均処理を行って波形３７を得る。加算平均処理をすることで、母親の心電計波形３５に同期しない胎児の心臓の磁気信号波形は減衰し、心電計波形３５に同期した母親の心臓の磁気信号の加算平均処理波形３７を抽出することができる。前記波形３７の任意の２点の時刻における信号強度の値を直線で結ぶ。前記直線を波形３７の基線３８とする。前記母親の腹部から計測された磁気信号波形３６は、母親の呼吸振動に起因する波形のゆれも含む場合があり、加算平均処理だけでは完全に前記波形のゆれは除去できない場合もある。結果として、波形３７のように基線３８は、傾きとオフセット分をもった波形となる。前記基線３８を基線４０のように傾きを０にし、かつオフセットの値が０になるよう基線補正処理を行って、波形３９を得る。前記基線補正処理を行うことによって、以下記述する差分処理を行うときに、差分区間の波形と波形のつなぎ目を平滑にする効果がある。次に、心電計信号波形３５（トリガー信号）に同期させて、母親の腹部から計測された磁気信号波形３６から母親の心臓磁気信号の基線補正処理波形３９を差分処理すると、目的とする胎児の心臓磁気信号の差分処理波形４１を得る。更に、波形３６から波形４１を差分処理すると、母親の長時間の心臓磁気信号波形（図示せず）を得ることができる。上記処理においては心電計等の外部入力信号があった方が処理が簡単になるが、そのような外部入力信号がなくても、磁気センサからの信号波形のみを、基準時刻に基づき処理することにより可能となる。
【００１８】図４は、本発明における信号処理の一実施例として臨床例を示した図である。一例として、胎児の心臓から発生する磁気信号計測を目的として、母親の腹部から計測された磁気信号波形４５（収集データ）は、胎児の心臓から発生する磁気信号波形だけでなく母体の心臓から発生する磁気信号波形が重なり合っている。外部入力信号として母体の心電計波形４６を取り込み、前記心電計波形４６をトリガー信号として、磁気信号波形４５について加算平均処理２５（図２）し、更に基線補正処理２６（図２）を施し得られた波形（参照データ）を、磁気信号波形４５と波形差分処理２７（図２）を行うと、母体の心臓の磁気信号が除去され、波形４７が得られる。波形４７は、胎児の心臓の磁気信号波形である。更に、胎児の心臓波形でＱＲＳ波と呼ばれる波形の除去を目的として、波形４７のＲ波の繰り返しピークをトリガー信号として、波形４７について再度加算平均処理２５（図２）と基線補正処理２６（図２）を施し得られた波形（参照データ）を、磁気信号波形４７と波形差分処理２７（図２）を行うと、胎児の心臓波形のＱＲＳ波を除去した波形４８を抽出することができる。この臨床例については、波形４８の特徴として鋸歯のような周期性を持った波形（Ｆ波と呼ばれる）が見られ、症例として心房粗動と診断される。
【００１９】図５は、表示画面の一例を示した図である。一例として、胎児の心臓から発生する磁気信号計測を目的として、母親の腹部から計測された磁気信号波形５１と、外部入力信号として取り込んだ母体の心電計波形５３と、前記心電計波形をトリガー信号として、磁気信号波形５１について加算平均処理２５（図２）と基線補正処理２６（図２）と波形差分処理２７（図２）を施し得られた胎児の心臓磁気波形５２を、図４のように同一画面に同一時系列に並べて表示することで、視覚的に母体信号と胎児の信号との時間的関係などが容易に比較できる。更に、磁気信号波形５１から胎児の心臓の磁気信号波形５２を差分処理すると、母体の心臓の磁気信号波形（図示せず）を抽出することができ、波形表示ウィンドウを追加して、前記抽出した母体の心臓磁気信号波形も同一画面に同一時系列に表示することができる。表示する波形チャンネルは、チャンネル選択ボックス５７の▼印をクリックするとプルダウンメニューが開き、その中からチャンネルを選択する。スクロールバー５５で波形表示の時刻５６をスクロールすることができる。また波形５１と波形５２の波形表示の時刻も連動して同時にスクロールする。カーソル５４を移動させることで、カーソル５４の位置での波形の信号強度および時刻を表示する。また、カーソルは一箇所を移動させると連動して、その他の波形のカーソルも移動する。
【００２０】
【発明の効果】上記に示した本発明によれば、ＳＱＵＩＤ磁束計から得られた胎児の心臓や脳から発生する磁気信号と母体から発生する磁気信号との分離や、脳から発生する磁気信号と心臓から発生する磁気信号の分離といった２つ以上の異なる部位から発生する磁気信号を分離することができ、心臓疾患や脳機能障害等の診断に有用な情報を提供することができる。
【図面の簡単な説明】
【図１】本発明の実施例である生体磁気計測装置の構成を示す図。
【図２】本発明の実施例である信号処理の概念フロー図。
【図３】本発明の実施例である波形基線補正処理の概念図。
【図４】本発明の実施例である信号処理の一臨床例を示す図。
【図５】本発明の実施例である波形表示画面図。
【符号の説明】
１…磁気シールドルーム、２…ＳＱＵＩＤ磁気センサおよびクライオスタット、３…ガントリー、４…ベッド、５…ＳＱＵＩＤ磁束計の駆動回路、６…アンプ回路およびフィルタ回路ユニット、７…コンピュータ、８…心電計などの外部入力信号を取るための回路、９…被験者、２１…生体磁気計測信号、２２…外部入力信号、２３…生データファイル、２４…信号処理波形データファイル、２５…加算平均処理、２６…基線補正処理、２７…波形差分処理（波形Ａ−波形Ｂ）、２８…波形差分処理（波形Ａ−波形Ｃ）、３０…波形表示処理、３１…等磁場線図処理、３２…表示画面部、３５…トリガー信号波形、３６…計測された磁気信号波形、３７…加算平均処理波形、３８…加算平均処理波形の基線、３９…基線補正処理波形、４０…基線補正処理波形の基線、４１…差分処理波形（波形３６から波形３９を差し引いた波形）、４５…計測された生体磁気信号波形（母体の心臓と胎児の心臓の磁気信号波形）、４６…計測された心電計信号波形（母体の心臓の信号波形）、４７…母体信号除去処理後の波形（胎児の心臓の磁気信号波形）、４８…胎児の心臓のＱＲＳ波を除去した磁気信号波形、５１…計測された生体磁気信号波形（母体の心臓と胎児の心臓の磁気信号波形）、５２…信号処理後波形（胎児の磁気信号波形）、５３…トリガー信号波形（母体の心電計波形）、５４…時刻と信号強度表示用カーソルバー、５５…時間波形スクロールバー、５６…時刻と時間スケール表示、５７…表示チャンネル選択ボックス。
【０００１】
【発明の属する技術分野】本発明は、成人，小児，胎児などの心臓や脳などから発生する微弱な磁気信号を計測するＳＱＵＩＤ(Superconducting Quantum Interference Device：超伝導量子干渉素子）磁束計を含めた生体磁気計測装置に関するものであり、特に臨床検査分野において、生体磁気計測装置から得られた生体磁気信号の情報から有用な情報を得る解析手段および表示手段を有する生体磁気計測装置に関するものである。
【０００２】
【従来の技術】従来技術として、成人，小児，胎児などの心臓や脳などから発生する生体磁気信号の検出はＳＱＵＩＤ(Superconducting Quantum Interference Device：超伝導量子干渉素子）磁束計を用いて検出されている。
【０００３】特開平１０−３０５０１９号公報には、生体から発する生体磁気信号を検出する複数の磁束計と、信号の演算処理を行う演算処理手段と演算結果を表示する手段とを有し、磁気シールドルーム内で生体磁場分布を計測する装置について記載されている。
【０００４】特開平１０−３０５０１９号公報には、生体磁気計測装置において、検出された生体磁気信号の生データからノイズを除去する目的で加算平均（アベレージング）処理を有し、処理条件を視覚的に簡便化した表示手段を有する生体磁気計測装置について記載されている。
【０００５】また、特開平５−２４９２０８号公報には、生体磁気信号を加算平均したものを測定信号から引き算し、雑音のみの信号を得る。雑音成分は主に冷凍機から発生する周期的ノイズであることから、該雑音信号を数サイクル分加算平均し、テンプレートを作成する。その後、測定信号から雑音信号のテンプレートを引き算することにより、生体磁気信号のみのデータを得るという方法が開示されている。
【０００６】
【発明が解決しようとする課題】検出された生体磁気信号は、目的とする部位以外の生体磁気信号が混ざっていることがあり、目的とする部位の磁気信号との区別が困難な場合がある。例えば、胎児の心臓から発生する磁気信号と母体から発生する磁気信号が、ほぼ同じ大きさで検出された場合などは、どちらから発生した磁気信号なのかを区別することは、従来の技術では困難であった。
【０００７】本発明では、ＳＱＵＩＤ磁束計から得られた胎児の心臓や脳から発生する磁気信号と母体から発生する磁気信号の分離や、脳から発生する磁気信号と心臓から発生する磁気信号の分離というような２つ以上の異なる部位から発生する信号の分離を目的とした信号処理の手段を有する生体磁気計測装置が提供でき、心臓疾患や脳機能障害等の診断に有用な情報を提供することが可能となる。
【０００８】
【課題を解決するための手段】上述の課題を解決するために、本発明によれば、被験者を支持するための支持体と、被験者から発生する磁場を検出する超伝導量子干渉素子（ＳＱＵＩＤ）を備える磁気センサと、該磁気センサを外部磁場から遮断するシールドと、該超伝導量子干渉素子を低温に保持するためのクライオスタットと、該クライオスタットを保持する手段と、前記超伝導量子干渉素子の駆動回路を含む磁束計回路とを備えた生体磁気計測装置であって、更に、被験者の測定対象である磁気発生源を、磁気以外の物理量の変化で測定する測定手段を備え、該測定手段からの信号と前記磁気センサからの信号を演算し、前記被験者の複数の磁気発生源からの信号を分離する演算装置を備えた生体磁気計測装置が提供される。
【０００９】上記生体磁気計測装置において、前記磁気以外の物理量変化の測定手段が被験者の皮膚の電位の変動を測定する手段であっても良い。
【００１０】また、本発明では、生体から発生する磁場を計測する超伝導量子干渉素子（Superconducting Quantum Interference Device： 以後ＳＱＵＩＤと略す）を用いた１個ないし複数個の磁束計と、該ＳＱＵＩＤ磁束計によって得られた磁気信号の波形を表示する手段と、該ＳＱＵＩＤ磁束計以外の信号、例えば心電計や脳波計などの信号を取り込み表示する手段と、該磁束計によって得られた磁気信号波形、あるいは前記心電計や脳波計などから取り込んだ該磁束計以外の外部入力信号波形を収集データとして、前記収集データのいずれかのチャンネルの信号波形から基準時刻（トリガー信号）を決定し、前記トリガー信号を用いて前記磁気信号波形を加算平均する手段と、前記加算平均手段により得られた磁気信号波形の任意の２点の時刻における磁気信号強度の値を結んだ直線（基線）の傾きが０になるように補正し、かつ前記任意の２点の時刻における磁気信号強度の値が０（オフセット）になるように補正する手段（基線補正手段）を有し、前記加算平均手段と前記基線補正手段によって得られた波形（波形Ｂとする）を前記トリガー信号を用いて、前記計測された磁気信号波形（波形Ａとする）から差し引く手段と、前記差し引く手段によって得られた波形（波形Ｃとする）を表示する手段を備えることができる。
【００１１】また、前記差し引く手段によって得られた波形（波形Ｃ）を、前記トリガー信号を用いて前記計測された磁気信号波形（波形Ａ）から差し引く手段（該手段によって得られる波形を波形Ｄとする）を有し、前記計測された磁気信号波形Ａと、前記波形Ｃと、前記波形Ｄと、前記トリガー信号波形のいずれか２つ以上を時系列に上下並べて同一画面上に表示できる手段を有したものである。
【００１２】更に、前記波形Ｃ、又は前記波形Ｄを前記収集データとして置き換え、上記加算平均手段や基線補正手段、前記差し引く手段などの演算処理を繰り返し行う手段を有し、前記手段によって得られた波形を表示する手段を備えることができる。
【００１３】
【発明の実施の形態】本発明の一実施例を図を用いて説明する。
【００１４】図１は、本発明の実施例である生体磁気計測装置の構成を示す図である。図１に示すように磁気シールドルーム１内には、被験者９が横になるベッド４と複数個（複数チャンネル）のＳＱＵＩＤ磁気センサおよびＳＱＵＩＤ磁気センサを超伝導状態に保持するための冷媒（液体ヘリウムまたは液体窒素）が貯蔵されたクライオスタット２と、クライオスタット２を機械的に保持するガントリー３が配置されている。ベッド４は、Ｘ方向，Ｙ方向，Ｚ方向に移動可能である。磁気シールドルーム１の外部にはＳＱＵＩＤ磁束計の駆動回路５と、アンプ回路およびフィルタ回路ユニット６と、データ取り込みおよびデータ解析用コンピュータ７と、心電計などの外部入力信号を取るための回路８とが配置されている。
【００１５】ＳＱＵＩＤ磁気センサによって検出された生体磁気信号は、アンプ回路およびフィルタ回路ユニット６により増幅され、かつ設定周波数より低い周波数信号を通過させるローパスフィルタや設定周波数より高い周波数信号を通過させるハイパスフィルタ，商用電源周波数だけをカットするノッチフィルタなどの信号処理を経た後、パソコン７に生データとして取り込まれる。また、心電計などの外部入力信号をとるための回路８からの信号波形はパソコン７に生データとして取り込まれる。また、パソコン７は、取り込んだ生データを生データファイルに格納し、波形を画面表示したり、また波形の信号処理や等磁場線図処理などを行い、表示することもできる。
【００１６】図２は、本発明の実施例であるデータ処理の概念フロー図を示したものである。図２において、図１R>１のパソコン７で取り込んだ生体磁気計測信号２１および外部入力信号２２は、生データファイル２３に格納される。また生体磁気計測信号２１および外部入力信号２２は、生データファイル２３に格納されると同時に、ＤＳＰ(Digital Signal Processor）(図示せず）を用いて、リアルタイムで加算平均処理２５などの信号処理を行うこともできる。生データファイル２３からデータ処理するファイルを呼び出し、呼び出したファイルの信号波形データのうち、磁気信号波形あるいは、該磁束計以外の入力信号波形のいずれかのチャンネルの波形をトリガー信号として、該トリガー信号波形を用いて生体磁気計測信号２１（波形Ａ）に対して加算平均処理２５を行う。次に、加算平均処理２５により得られた波形に対して、基線補正処理２６を行う。次に、基線補正処理２６により得られた波形（波形Ｂ）を、前記トリガー信号を用いて、波形Ａとの差分処理２７を行う。波形差分処理２７により得られた波形Ｃは、信号処理波形データファイル２４に格納される。また波形Ｃは、波形Ａとの差分処理２８により、波形Ｄが得られる。波形Ｄは、信号処理波形データファイル２４に格納される。前記信号処理波形データファイル２４のデータは、再度、加算平均処理２５，基線補正処理２６，差分処理２７，差分処理２８の一連のデータ処理を繰り返し行うこともできる。また、トリガー信号として前記信号処理波形データファイル２４のデータを用いることもできる。生データファイル２３，信号処理波形データファイル２４に格納された波形データは、表示区間や表示スケール，表示する波形のチャンネルなどの入力設定に従って、波形表示処理３０および等磁場線図処理３１が行われ、表示画面部３２に表示される。
【００１７】図３は、図２の加算平均処理２５と基線補正処理２６と波形差分処理２７の信号処理方法の概念を示したものである。一実施例として、胎児の心臓磁気信号計測を目的とした場合について説明する。母親の腹部から計測された磁気信号波形３６について、母親の心電計波形３５をトリガー信号として、トリガー信号のピーク値の時刻のＴ１秒前からＴ２秒間を加算区間として、加算平均処理を行う。一般に母親の心電計波形には、胎児心電計波形はほとんど検出できないか、あるいは極めて小さい。したがって、心電計の波形に閾値を設定して、波形のピーク値を検出すれば、母親のみの心電計波形のピーク値を検出できる。前記ピーク値をトリガーとして、加算区間Ｔ２で加算平均処理を行って波形３７を得る。加算平均処理をすることで、母親の心電計波形３５に同期しない胎児の心臓の磁気信号波形は減衰し、心電計波形３５に同期した母親の心臓の磁気信号の加算平均処理波形３７を抽出することができる。前記波形３７の任意の２点の時刻における信号強度の値を直線で結ぶ。前記直線を波形３７の基線３８とする。前記母親の腹部から計測された磁気信号波形３６は、母親の呼吸振動に起因する波形のゆれも含む場合があり、加算平均処理だけでは完全に前記波形のゆれは除去できない場合もある。結果として、波形３７のように基線３８は、傾きとオフセット分をもった波形となる。前記基線３８を基線４０のように傾きを０にし、かつオフセットの値が０になるよう基線補正処理を行って、波形３９を得る。前記基線補正処理を行うことによって、以下記述する差分処理を行うときに、差分区間の波形と波形のつなぎ目を平滑にする効果がある。次に、心電計信号波形３５（トリガー信号）に同期させて、母親の腹部から計測された磁気信号波形３６から母親の心臓磁気信号の基線補正処理波形３９を差分処理すると、目的とする胎児の心臓磁気信号の差分処理波形４１を得る。更に、波形３６から波形４１を差分処理すると、母親の長時間の心臓磁気信号波形（図示せず）を得ることができる。上記処理においては心電計等の外部入力信号があった方が処理が簡単になるが、そのような外部入力信号がなくても、磁気センサからの信号波形のみを、基準時刻に基づき処理することにより可能となる。
【００１８】図４は、本発明における信号処理の一実施例として臨床例を示した図である。一例として、胎児の心臓から発生する磁気信号計測を目的として、母親の腹部から計測された磁気信号波形４５（収集データ）は、胎児の心臓から発生する磁気信号波形だけでなく母体の心臓から発生する磁気信号波形が重なり合っている。外部入力信号として母体の心電計波形４６を取り込み、前記心電計波形４６をトリガー信号として、磁気信号波形４５について加算平均処理２５（図２）し、更に基線補正処理２６（図２）を施し得られた波形（参照データ）を、磁気信号波形４５と波形差分処理２７（図２）を行うと、母体の心臓の磁気信号が除去され、波形４７が得られる。波形４７は、胎児の心臓の磁気信号波形である。更に、胎児の心臓波形でＱＲＳ波と呼ばれる波形の除去を目的として、波形４７のＲ波の繰り返しピークをトリガー信号として、波形４７について再度加算平均処理２５（図２）と基線補正処理２６（図２）を施し得られた波形（参照データ）を、磁気信号波形４７と波形差分処理２７（図２）を行うと、胎児の心臓波形のＱＲＳ波を除去した波形４８を抽出することができる。この臨床例については、波形４８の特徴として鋸歯のような周期性を持った波形（Ｆ波と呼ばれる）が見られ、症例として心房粗動と診断される。
【００１９】図５は、表示画面の一例を示した図である。一例として、胎児の心臓から発生する磁気信号計測を目的として、母親の腹部から計測された磁気信号波形５１と、外部入力信号として取り込んだ母体の心電計波形５３と、前記心電計波形をトリガー信号として、磁気信号波形５１について加算平均処理２５（図２）と基線補正処理２６（図２）と波形差分処理２７（図２）を施し得られた胎児の心臓磁気波形５２を、図４のように同一画面に同一時系列に並べて表示することで、視覚的に母体信号と胎児の信号との時間的関係などが容易に比較できる。更に、磁気信号波形５１から胎児の心臓の磁気信号波形５２を差分処理すると、母体の心臓の磁気信号波形（図示せず）を抽出することができ、波形表示ウィンドウを追加して、前記抽出した母体の心臓磁気信号波形も同一画面に同一時系列に表示することができる。表示する波形チャンネルは、チャンネル選択ボックス５７の▼印をクリックするとプルダウンメニューが開き、その中からチャンネルを選択する。スクロールバー５５で波形表示の時刻５６をスクロールすることができる。また波形５１と波形５２の波形表示の時刻も連動して同時にスクロールする。カーソル５４を移動させることで、カーソル５４の位置での波形の信号強度および時刻を表示する。また、カーソルは一箇所を移動させると連動して、その他の波形のカーソルも移動する。
【００２０】
【発明の効果】上記に示した本発明によれば、ＳＱＵＩＤ磁束計から得られた胎児の心臓や脳から発生する磁気信号と母体から発生する磁気信号との分離や、脳から発生する磁気信号と心臓から発生する磁気信号の分離といった２つ以上の異なる部位から発生する磁気信号を分離することができ、心臓疾患や脳機能障害等の診断に有用な情報を提供することができる。
【図面の簡単な説明】
【図１】本発明の実施例である生体磁気計測装置の構成を示す図。
【図２】本発明の実施例である信号処理の概念フロー図。
【図３】本発明の実施例である波形基線補正処理の概念図。
【図４】本発明の実施例である信号処理の一臨床例を示す図。
【図５】本発明の実施例である波形表示画面図。
【符号の説明】
１…磁気シールドルーム、２…ＳＱＵＩＤ磁気センサおよびクライオスタット、３…ガントリー、４…ベッド、５…ＳＱＵＩＤ磁束計の駆動回路、６…アンプ回路およびフィルタ回路ユニット、７…コンピュータ、８…心電計などの外部入力信号を取るための回路、９…被験者、２１…生体磁気計測信号、２２…外部入力信号、２３…生データファイル、２４…信号処理波形データファイル、２５…加算平均処理、２６…基線補正処理、２７…波形差分処理（波形Ａ−波形Ｂ）、２８…波形差分処理（波形Ａ−波形Ｃ）、３０…波形表示処理、３１…等磁場線図処理、３２…表示画面部、３５…トリガー信号波形、３６…計測された磁気信号波形、３７…加算平均処理波形、３８…加算平均処理波形の基線、３９…基線補正処理波形、４０…基線補正処理波形の基線、４１…差分処理波形（波形３６から波形３９を差し引いた波形）、４５…計測された生体磁気信号波形（母体の心臓と胎児の心臓の磁気信号波形）、４６…計測された心電計信号波形（母体の心臓の信号波形）、４７…母体信号除去処理後の波形（胎児の心臓の磁気信号波形）、４８…胎児の心臓のＱＲＳ波を除去した磁気信号波形、５１…計測された生体磁気信号波形（母体の心臓と胎児の心臓の磁気信号波形）、５２…信号処理後波形（胎児の磁気信号波形）、５３…トリガー信号波形（母体の心電計波形）、５４…時刻と信号強度表示用カーソルバー、５５…時間波形スクロールバー、５６…時刻と時間スケール表示、５７…表示チャンネル選択ボックス。
【特許請求の範囲】
【請求項１】被験者を支持するための支持体と、被験者から発生する磁場を検出する超伝導量子干渉素子（ＳＱＵＩＤ）を備える磁気センサと、該磁気センサを外部磁場から遮断するシールドと、該超伝導量子干渉素子を低温に保持するためのクライオスタットと、該クライオスタットを保持する手段と、前記超伝導量子干渉素子の駆動回路を含む磁束計回路とを備えた生体磁気計測装置であって、更に、被験者の測定対象である磁気発生源を、磁気以外の物理量の変化で測定する測定手段を備え、該測定手段からの信号と前記磁気センサからの信号を演算し、前記磁気センサからの信号に含まれる前記被験者の複数の磁気発生源からの信号を分離する演算装置を備えたことを特徴とする生体磁気計測装置。
【請求項２】請求項１記載の生体磁気計測装置において、前記磁気以外の物理量変化の測定手段が被験者の皮膚の電位の変動を測定する手段であることを特徴とする生体磁気計測装置。
【請求項３】請求項１または２記載の生体磁気計測装置において、予め用意された参照データと、前記磁気センサ及び前記磁気以外の物理量変化の測定手段から得られた収集データの間で演算処理を行い、前記収集データの信号波形から基準時刻を決定し、この同一の基準時刻に基づいて、前記参照データと前記収集データの間で演算処理を行う演算装置を備えたことを特徴とする生体磁気計測装置。
【請求項４】請求項３記載の参照データは収集データに基づき算出されたものであり、該参照データと前記収集データとの間の演算処理は、加算又は減算であることを特徴とする生体磁気計測装置。
【請求項５】請求項４記載の生体磁気計測装置において、繰り返し波形を有する入力信号からなる収集データについて、該繰り返し波形の一部を単位波形として、複数の単位波形を切り出して加算平均し、各々の単位波形切り出し区間又は前記加算平均区間を各チャンネル共に同一の基準時刻に基づいて決定し、且つ各々の入力信号から、各チャンネル毎の上記加算平均信号を上記基準時刻に基づいて減算する演算装置を備えたことを特徴とする生体磁気計測装置。
【請求項６】請求項５記載の生体磁気計測装置において、前記加算平均により得られた磁気信号波形について、任意の２点の時刻における磁気信号強度の値を結んだ直線（基線とする）の傾きが０、かつ前記任意の２点における磁気信号強度の値が０（オフセット）になるように補正演算する手段を有し、前記加算平均手段と前記基線補正手段によって得られた波形（波形Ｂとする）を、前記計測された磁気信号波形から差し引く手段と、前記差し引く手段によって得られた波形（波形Ｃとする）を表示する手段を備えたことを特徴とする生体磁気計測装置。
【請求項７】請求項６記載の生体磁気計測装置において、前記差し引く手段によって得られた波形データ（波形データＣ）を、前記基準時刻を用いて前記収集データから差し引く手段（該手段によって得られる波形データを波形データＤとする）を備えたことを特徴とする生体磁気計測装置。
【請求項８】請求項７記載の生体磁気計測装置において、前記波形データＣ、または前記波形データＤを収集データとして置き換え、再度演算処理を繰り返し行う手段を備えたことを特徴とする生体磁気計測装置。
【請求項９】請求項１〜８のいずれかに記載の磁気計測装置において、収集データ，前記波形データＣ，前記波形データＤ、前記基準時刻を示すデータのいずれか２つ以上を時系列に上下並べて同一画面上に表示できる手段を備えたことを特徴とする生体磁気計測装置。
【請求項１０】被験者を支持するための支持体と、被験者から発生する磁場を検出する超伝導量子干渉素子（ＳＱＵＩＤ）を備える磁気センサと、該磁気センサを外部磁場から遮断するシールドと、該超伝導量子干渉素子を低温に保持するためのクライオスタットと、該クライオスタットを保持する手段と、前記超伝導量子干渉素子の駆動回路を含む磁束計回路とを備えた生体磁気計測装置であって、前記磁気センサからの信号からなる収集データの信号波形に基づいて基準時刻を決定し、予め用意された参照データ収集データの信号波形から基準時刻を決定し、該基準時刻に基づいて、前記参照データと前記収集データの間で演算処理を行う演算装置を備えたことを特徴とする生体磁気計測装置。
【請求項１１】請求項１０記載の生体磁気計測装置において、前記収集データの繰り返し波形の一部を単位波形として、複数の単位波形を切り出して加算平均し、かつ該加算平均して得られた信号波形について、任意の２点の時刻における信号強度の値を結んだ直線の傾きが０、かつ前記任意の２点における心経強度が０になるように基線補正し、前記加算平均と前記基線補正によって得られた波形(波形Ｂとする)を、前記計測された磁気信号波形から差し引く演算を行う演算装置を備えたことを特徴とする生体磁気計測装置。
【請求項１２】被験者から発生する磁場を超伝導量子干渉素子（ＳＱＵＩＤ）を備える磁気センサで測定する生体磁場計測方法において、被験者の測定対象である磁気発生源を、磁気以外の物理量の変化で測定した測定結果と、前記磁気センサからの測定結果を演算し、前記被験者の複数の磁気発生源からの信号を分離することを特徴とする生体磁気計測方法。
【請求項１】被験者を支持するための支持体と、被験者から発生する磁場を検出する超伝導量子干渉素子（ＳＱＵＩＤ）を備える磁気センサと、該磁気センサを外部磁場から遮断するシールドと、該超伝導量子干渉素子を低温に保持するためのクライオスタットと、該クライオスタットを保持する手段と、前記超伝導量子干渉素子の駆動回路を含む磁束計回路とを備えた生体磁気計測装置であって、更に、被験者の測定対象である磁気発生源を、磁気以外の物理量の変化で測定する測定手段を備え、該測定手段からの信号と前記磁気センサからの信号を演算し、前記磁気センサからの信号に含まれる前記被験者の複数の磁気発生源からの信号を分離する演算装置を備えたことを特徴とする生体磁気計測装置。
【請求項２】請求項１記載の生体磁気計測装置において、前記磁気以外の物理量変化の測定手段が被験者の皮膚の電位の変動を測定する手段であることを特徴とする生体磁気計測装置。
【請求項３】請求項１または２記載の生体磁気計測装置において、予め用意された参照データと、前記磁気センサ及び前記磁気以外の物理量変化の測定手段から得られた収集データの間で演算処理を行い、前記収集データの信号波形から基準時刻を決定し、この同一の基準時刻に基づいて、前記参照データと前記収集データの間で演算処理を行う演算装置を備えたことを特徴とする生体磁気計測装置。
【請求項４】請求項３記載の参照データは収集データに基づき算出されたものであり、該参照データと前記収集データとの間の演算処理は、加算又は減算であることを特徴とする生体磁気計測装置。
【請求項５】請求項４記載の生体磁気計測装置において、繰り返し波形を有する入力信号からなる収集データについて、該繰り返し波形の一部を単位波形として、複数の単位波形を切り出して加算平均し、各々の単位波形切り出し区間又は前記加算平均区間を各チャンネル共に同一の基準時刻に基づいて決定し、且つ各々の入力信号から、各チャンネル毎の上記加算平均信号を上記基準時刻に基づいて減算する演算装置を備えたことを特徴とする生体磁気計測装置。
【請求項６】請求項５記載の生体磁気計測装置において、前記加算平均により得られた磁気信号波形について、任意の２点の時刻における磁気信号強度の値を結んだ直線（基線とする）の傾きが０、かつ前記任意の２点における磁気信号強度の値が０（オフセット）になるように補正演算する手段を有し、前記加算平均手段と前記基線補正手段によって得られた波形（波形Ｂとする）を、前記計測された磁気信号波形から差し引く手段と、前記差し引く手段によって得られた波形（波形Ｃとする）を表示する手段を備えたことを特徴とする生体磁気計測装置。
【請求項７】請求項６記載の生体磁気計測装置において、前記差し引く手段によって得られた波形データ（波形データＣ）を、前記基準時刻を用いて前記収集データから差し引く手段（該手段によって得られる波形データを波形データＤとする）を備えたことを特徴とする生体磁気計測装置。
【請求項８】請求項７記載の生体磁気計測装置において、前記波形データＣ、または前記波形データＤを収集データとして置き換え、再度演算処理を繰り返し行う手段を備えたことを特徴とする生体磁気計測装置。
【請求項９】請求項１〜８のいずれかに記載の磁気計測装置において、収集データ，前記波形データＣ，前記波形データＤ、前記基準時刻を示すデータのいずれか２つ以上を時系列に上下並べて同一画面上に表示できる手段を備えたことを特徴とする生体磁気計測装置。
【請求項１０】被験者を支持するための支持体と、被験者から発生する磁場を検出する超伝導量子干渉素子（ＳＱＵＩＤ）を備える磁気センサと、該磁気センサを外部磁場から遮断するシールドと、該超伝導量子干渉素子を低温に保持するためのクライオスタットと、該クライオスタットを保持する手段と、前記超伝導量子干渉素子の駆動回路を含む磁束計回路とを備えた生体磁気計測装置であって、前記磁気センサからの信号からなる収集データの信号波形に基づいて基準時刻を決定し、予め用意された参照データ収集データの信号波形から基準時刻を決定し、該基準時刻に基づいて、前記参照データと前記収集データの間で演算処理を行う演算装置を備えたことを特徴とする生体磁気計測装置。
【請求項１１】請求項１０記載の生体磁気計測装置において、前記収集データの繰り返し波形の一部を単位波形として、複数の単位波形を切り出して加算平均し、かつ該加算平均して得られた信号波形について、任意の２点の時刻における信号強度の値を結んだ直線の傾きが０、かつ前記任意の２点における心経強度が０になるように基線補正し、前記加算平均と前記基線補正によって得られた波形(波形Ｂとする)を、前記計測された磁気信号波形から差し引く演算を行う演算装置を備えたことを特徴とする生体磁気計測装置。
【請求項１２】被験者から発生する磁場を超伝導量子干渉素子（ＳＱＵＩＤ）を備える磁気センサで測定する生体磁場計測方法において、被験者の測定対象である磁気発生源を、磁気以外の物理量の変化で測定した測定結果と、前記磁気センサからの測定結果を演算し、前記被験者の複数の磁気発生源からの信号を分離することを特徴とする生体磁気計測方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図２】
【図３】
【図４】
【図５】
【公開番号】特開２００２−２９１７１３（Ｐ２００２−２９１７１３Ａ）
【公開日】平成１４年１０月８日（２００２．１０．８）
【国際特許分類】
【出願番号】特願２００１−９４５５５（Ｐ２００１−９４５５５）
【出願日】平成１３年３月２９日（２００１．３．２９）
【出願人】（０００００５１０８）株式会社日立製作所 (27,607)
【Ｆターム（参考）】
【公開日】平成１４年１０月８日（２００２．１０．８）
【国際特許分類】
【出願日】平成１３年３月２９日（２００１．３．２９）
【出願人】（０００００５１０８）株式会社日立製作所 (27,607)
【Ｆターム（参考）】
[ Back to top ]