特許 7439599 情報処理装置及び雑音低減方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2024-02-19

(45)【発行日】2024-02-28

(54)【発明の名称】情報処理装置及び雑音低減方法

(51)【国際特許分類】

   A61B 5/242 20210101AFI20240220BHJP

   A61B 5/243 20210101ALI20240220BHJP

   A61B 5/245 20210101ALI20240220BHJP

   A61B 5/248 20210101ALI20240220BHJP

【ＦＩ】

A61B5/242

A61B5/243

A61B5/245

A61B5/248

【請求項の数】 7

(21)【出願番号】P 2020050207

(22)【出願日】2020-03-19

(65)【公開番号】P2021145987

(43)【公開日】2021-09-27

【審査請求日】2023-01-16

(73)【特許権者】

【識別番号】000006747

【氏名又は名称】株式会社リコー

(74)【代理人】

【識別番号】100089118

【弁理士】

【氏名又は名称】酒井 宏明

(72)【発明者】

【氏名】前田 亮

【審査官】藤原 伸二

(56)【参考文献】

【文献】国際公開第２０１９／０４７１６５（ＷＯ，Ａ１）

【文献】特表２０１９－５３１８１５（ＪＰ，Ａ）

【文献】特開平１０－２０９８８９（ＪＰ，Ａ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ａ６１Ｂ ５／００－５／３９８

Ｇ０１Ｒ ３３／００－３３／２６

Ｇ０１Ｎ ２４／００－２４／１４

(57)【特許請求の範囲】

【請求項1】

磁気センサによって検知された生体磁気信号と雑音とを含む磁場信号を取得する手段と、
前記磁場信号に含まれる雑音の周波数を参照周波数として、前記磁場信号に含まれる前記周波数の雑音の振幅と位相とを検出する手段と、
前記磁場信号から前記振幅と位相とを有する正弦波を減算する手段と
を備えたことを特徴とする情報処理装置。

【請求項2】

生体磁気信号と雑音とを含む磁場信号を取得する磁気センサと、
参照磁気信号を取得するリファレンスセンサと、
前記リファレンスセンサによって取得された前記参照磁気信号における雑音の周波数を参照周波数として、前記参照磁気信号に含まれる前記周波数の雑音の振幅と位相とを検出する手段と、
前記振幅と前記位相を持つ正弦波と前記磁場信号から重み係数を計算する手段と、
前記磁場信号から、前記重み係数を乗算した正弦波を減算する手段と
を備えたことを特徴とする情報処理装置。

【請求項3】

前記重み係数は、前記振幅と前記位相を持つ正弦波と、前記磁気センサで環境ノイズのみを取得した磁気信号とから計算される
ことを特徴とする請求項２に記載の情報処理装置。

【請求項4】

前記重み係数は、前記磁気センサによって検出された信号に対して複数の期間に渡って計算され、その間の区間は、線形補間して時間的に重み係数を変化させる
ことを特徴とする請求項２に記載の情報処理装置。

【請求項5】

前記生体磁気信号は、脳、神経、心臓又は筋肉からの磁気信号である
ことを特徴とする請求項１～４のいずれか１つに記載の情報処理装置。

【請求項6】

情報処理装置で実行される雑音低減方法であって、
磁気センサによって検知された生体磁気信号と雑音とを含む磁場信号を取得する工程と、
前記磁場信号に含まれる雑音の周波数を参照周波数として、前記磁場信号に含まれる前記周波数の雑音の振幅と位相とを検出する工程と、
前記磁場信号から前記振幅と位相とを有する正弦波を減算する工程と
を含むことを特徴とする雑音低減方法。

【請求項7】

情報処理装置で実行される雑音低減方法であって、
磁気センサによって検知された生体磁気信号と雑音とを含む磁場信号を取得する工程と、
リファレンスセンサによって参照磁気信号を取得する工程と、
前記リファレンスセンサによって取得された前記参照磁気信号における雑音の周波数を参照周波数として、前記参照磁気信号に含まれる前記周波数の雑音の振幅と位相とを検出する工程と、
前記振幅と前記位相を持つ正弦波と前記磁場信号から重み係数を計算する工程と、
前記磁場信号から、前記重み係数を乗算した正弦波を減算する工程と
を含むことを特徴とする雑音低減方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、情報処理装置及び雑音低減方法に関する。

【背景技術】

【0002】

従来、生体磁気計測は、生体の電気信号から発生する極めて微弱な磁場を計測する。一般に、生体磁気計測では、磁気シールドルーム（Magnetic Shield Room：ＭＳＲ）内で計測することで地磁気や商用電源からの輻射ノイズ（以下では、ノイズを雑音とも表記する場合がある）等の環境磁場の影響を低減しているが、ある程度のノイズが残る場合がある。

【0003】

例えば、非特許文献１には、残留磁気ノイズを低減するために、生体磁気計測用センサの他に磁気ノイズ測定用のセンサを使用して、生体磁気信号から重み付けされた磁気ノイズ信号を差し引くことで生体磁気信号に含まれる磁気ノイズを低減する技術が開示されている。

【発明の概要】

【発明が解決しようとする課題】

【0004】

しかしながら、上記従来技術では、磁気ノイズ信号に含まれるホワイトノイズ等の統計的ノイズも加減されるため生体磁気信号のホワイトノイズが増加してしまう問題があった。

【0005】

本発明は、上記に鑑みてなされたものであって、雑音を低減することを目的とする。

【課題を解決するための手段】

【0006】

上述した課題を解決し、目的を達成するために、本発明は、磁気センサによって検知された生体磁気信号と雑音とを含む磁場信号を取得する手段と、前記磁場信号に含まれる雑音の周波数を参照周波数として、前記磁場信号に含まれる前記周波数の雑音の振幅と位相とを検出する手段と、前記磁場信号から前記振幅と位相とを有する正弦波を減算する手段とを備えたことを特徴とする。

【発明の効果】

【0007】

本発明によれば、雑音を低減することができるという効果を奏する。

【図面の簡単な説明】

【0008】

【図1】図１は、実施形態における生体磁気計測システムの構成例を示す図である。

【図2】図２は、実施形態における磁気ノイズスペクトルを示す図である。

【図3】図３は、実施形態におけるノイズ除去処理の動作一例を示すフローチャートである。

【図4】図４は、実施形態におけるロックイン検波の動作フローの概念図である。

【図5】図５は、実施形態におけるノイズを除去した磁気ノイズスペクトルを示す図である。

【図6】図６は、実施形態における生体磁気用センサとリファレンスセンサとの位置関係を示す図である。

【図7】図７は、実施形態におけるリファレンスセンサを利用したロックイン検波による特定周波数ノイズ除去の動作フローの概念図である。

【図8】図８は、実施形態における情報処理装置の構成例を示す図である。

【図9】図９は、情報処理装置の機能を実現するコンピュータの一例を示すハードウェア構成図である。

【発明を実施するための形態】

【0009】

以下に添付図面を参照して、この発明にかかる情報処理装置及び雑音低減方法の一実施形態を詳細に説明する。

【0010】

実施形態の雑音低減装置は、所定のセンサによって検知された時間的に変動する信号と雑音とを含む物理量を取得する物理量取得手段と、物理量に含まれる雑音の周波数を取得する周波数取得手段と、周波数を参照周波数として、物理量に含まれる周波数の雑音の振幅と位相とを検出するロックイン検波手段と、物理量から振幅と位相とを有する正弦波を減算する減算手段とを備えたことを特徴とする。

【0011】

〔１．本発明の概要〕
以下で説明する第１の実施形態及び第２の実施形態では、生体磁気計測システムが商用電源から放射される磁気ノイズを低減する方法を例に挙げて説明する。これにより、本発明は、磁気ノイズ信号に含まれる特定の周波数のノイズ信号のみを効果的に抽出することで、脳、神経、心臓又は筋肉からの生体磁気信号のホワイトノイズの増加を最低限に抑えつつ、生体磁気信号に含まれる磁気ノイズを低減することができる。なお、第１の実施形態及び第２の実施形態は一例であり、本発明は、音声信号や、脳、神経、心臓又は筋肉からの生体電位信号等の時間変化する如何なる信号にも応用可能である。

【0012】

〔２．第１の実施形態〕
図１は、実施形態における生体磁気計測システム１の構成例を示す図である。図１に示すように、生体磁気計測システム１は、ＭＳＲ２内に、所定のセンサとして、生体磁気信号を検知する生体磁気用センサ１１及び環境ノイズを検知するリファレンスセンサ１２を有するクライオスタット２０と、被測定者ＰＡ１を適切な位置に保持するためのベッド３０とを含む。ここで、環境ノイズは、如何なる由来の磁場によるノイズを含む。

【0013】

例えば、生体磁気用センサ１１とリファレンスセンサ１２とは、超電導量子干渉計（Superconducting Quantum Interference Device：ＳＱＵＩＤ）センサである。この場合、生体磁気用センサ１１とリファレンスセンサ１２とは、超電導状態に保つために、極低温状態を維持するクライオスタット２０の内部に配置される。

【0014】

また、生体磁気計測システム１は、ＭＳＲ外に、駆動回路１００と、収録回路２００と、情報処理装置３００（雑音低減装置の一例）とをさらに含む。例えば、駆動回路１００は、生体磁気用センサ１１とリファレンスセンサ１２とを駆動する。収録回路２００は、駆動回路１００から出力された信号を取得する。このとき、収録回路２００は、各信号を増幅してフィルタ等をかけてもよい。

【0015】

例えば、情報処理装置３００は、収録回路２００からの出力信号を収集し、収集された出力信号を解析する。なお、情報処理装置３００は、外部に表示装置を有してもよい。例えば、情報処理装置３００は、解析した出力信号等の出力データを表示する外部の表示装置として、ディスプレイ等を有してもよい。

【0016】

図２は、実施形態における磁気ノイズスペクトルを示す図である。図２には、生体磁気用センサ１１で検知された磁気ノイズのスペクトル密度が示される。図２に示されるように、商用電源の周波数である５０Ｈｚのノイズや、複数の周波数にわたってスペクトルの線幅が狭いノイズ等が検知される。

【0017】

ここで、本手法を適用する。図３は、実施形態におけるノイズ除去処理の一例を示すフローチャートである。例えば、５０Ｈｚの商用電源由来のノイズを低減するものとする。この場合、情報処理装置３００は、ロックイン検波の参照周波数を５０Ｈｚに設定する（ステップＳ１０１）。

【0018】

また、情報処理装置３００は、生体磁気用センサ１１によって検知された磁気信号に含まれる５０Ｈｚ成分の振幅と位相とを算出する（ステップＳ１０２）。そして、情報処理装置３００は、５０Ｈｚの正弦波を復元して、磁気信号から差し引くことで５０Ｈｚのノイズのみを除去する（ステップＳ１０３）。このように、本手法は、所望のノイズのみを除去することができる。

【0019】

次に、数式を用いて、ロックイン検波の動作フローを説明する。図４は、実施形態におけるロックイン検波の動作フローの概念図である。ここで、ロックイン検波の動作フローで用いられるローパスフィルタ（Low-pass filter：ＬＰＦ）のカットオフ周波数は、０．１～３Ｈｚ程度であるものとする。例えば、磁気信号に含まれる複数の周波数のノイズを、未知の振幅Ａ_ｎと未知の位相φ_ｎとを有する余弦波の和として示す。

【0020】

【数1】

ここで、磁気ノイズに特定の角周波数ω_ｉを有する正弦波を乗算する（ステップＳ２０１）。

【0021】

【数2】

式（２）の第１項を展開する。

【0022】

【数3】

また、式（２）の第２項をｎ＝ｊとして展開する。

【0023】

【数4】

【0024】

続いて、磁気ノイズに特定の角周波数ω_ｉを有する余弦波を乗算する（ステップＳ２０２）。

【0025】

【数5】

式（５）の第１項を展開する。

【0026】

【数6】

また、式（５）の第二項をｎ＝ｊとして展開する。

【0027】

【数7】

【0028】

そして、正弦波が乗算された磁気ノイズにローパスフィルタをかける（ステップＳ２０３）と、式（２）は、以下のように示される。

【0029】

【数8】

【0030】

また、余弦波が乗算された磁気ノイズにローパスフィルタをかける（ステップＳ２０４）と、式（５）は、以下のように示される。

【0031】

【数9】

【0032】

式（８）及び（９）より、磁気ノイズ信号に含まれる特定の周波数のノイズの振幅（ステップＳ２０５）と位相（ステップＳ２０６）とは、以下のように示される。

【0033】

【数10】

【0034】

【数11】

【0035】

以上のように、磁気ノイズ信号に含まれる特定の周波数のノイズの振幅と位相とを一意に算出することができる。ここで、周波数ｆ_ｉの前後にローパスフィルタの特性分のノイズが残るが、全帯域のノイズレベルと比較すると非常に軽微なレベルである。

【0036】

そして、ノイズを除去した磁気ノイズスペクトルを示す。図５は、実施形態におけるノイズを除去した磁気ノイズスペクトルを示す図である。図５は、図３から５０Ｈｚのノイズを除去した後のスペクトル密度を示す。図５に示すように、５０Ｈｚのノイズのみ大幅に減少しているのがわかる。なお、同一の磁気信号に対して複数の参照周波数でロックイン検波することで複数の周波数のノイズを除去することも可能である。

【0037】

〔３．第２の実施形態〕
第２の実施形態では、生体磁気信号が比較的大きい場合を想定する。この場合、生体磁気信号によってノイズの波形が乱されるため、特定の周波数のノイズの検出精度が低下する。これにより、雑音低減（ノイズリダクション）の性能が落ちる場合がある。そこで、生体磁気信号を検知する磁気センサの他に、環境ノイズを検知するリファレンスセンサ１２を用いることで、上記課題を解決する。

【0038】

ます、生体磁気用センサ１１とリファレンスセンサ１２との位置関係について説明する。図６は、実施形態における生体磁気用センサ１１とリファレンスセンサ１２との位置関係を示す図である。例えば、図６の例では、生体磁気用センサ１１とリファレンスセンサ１２とは、ＳＱＵＩＤセンサである。

【0039】

例えば、生体磁気用センサ１１は、被測定者の計測部位に近接して配置される。リファレンスセンサ１２は、生体磁気用センサ１１から所定の距離以上離れた位置に配置され、可能な限り、環境ノイズのみが検知できるように配置される。

【0040】

なお、図６の例では、生体磁気用センサ１１とリファレンスセンサ１２とがＳＱＵＩＤセンサである例を挙げて説明したが、これに限定されない。例えば、生体磁気用センサ１１とリファレンスセンサ１２とは、光ポンピング磁気センサ（Optical Pumping Magnetometer：ＯＰＭ）や、磁気抵抗（Magnetic Resistance：ＭＲ）センサや、磁気インピーダンス（Magnetic Impedance：ＭＩ）センサや、フラックスゲート等でもよい。

【0041】

また、上記例では、生体磁気用センサ１１とリファレンスセンサ１２とが同一種類のセンサである例を挙げて説明したが、これに限定されない。例えば、生体磁気用センサ１１とリファレンスセンサ１２とは、それぞれ異なる種類のセンサであってもよい。

【0042】

次に、リファレンスセンサ１２を使ったロックイン検波による特定の周波数のノイズを除去する動作フローについて説明する。図７は、実施形態におけるリファレンスセンサ１２を利用したロックイン検波による特定周波数ノイズ除去の動作フローの概念図である。図７の例では、商用電源から放射される５０Ｈｚのノイズを例に挙げて説明する。例えば、情報処理装置３００は、リファレンスセンサ１２によって検知された環境ノイズに対して、第１の実施形態と同様に５０Ｈｚを参照信号としてロックイン検波を行う（ステップＳ３０１）。また、情報処理装置３００は、リファレンスセンサ１２によって検知された環境ノイズに含まれる５０Ｈｚのノイズの振幅と位相とを検出し、得られた振幅と位相とを有する正弦波を生成する（ステップＳ３０２）。

【0043】

そして、情報処理装置３００は、生成された正弦波と生体磁気信号とから重み係数を計算する（ステップＳ３０３）。ここで、重み係数が時間変動しない場合ものとする。この場合、以下の方法を採用することができる。

【0044】

まず、生体磁気用センサ１１とリファレンスセンサ１２とは、環境ノイズのみを同時に検知する。続いて、情報処理装置３００は、ロックイン検波によって５０Ｈｚのノイズの振幅と位相とを得る。

【0045】

そして、情報処理装置３００は、重み係数を計算する。例えば、生体磁気用センサ１１とリファレンスセンサ１２との環境ノイズがサンプル数ｎの離散データである場合に、ｉ番目の生体磁気用センサ１１のデータをｓ（ｉ）、リファレンスセンサ１２に含まれる５０Ｈｚのノイズデータをｒ_５０（ｉ）、重み係数をｗとすると、以下のように示すことができる。

【0046】

【数12】

【0047】

続いて、情報処理装置３００は、式（１２）が最小になるように、重み係数ｗを算出する。このように、情報処理装置３００は、事前に算出された重み係数ｗを使用することで、生体磁気信号から５０Ｈｚのノイズを低減することができる。

【0048】

なお、生体磁気信号の波形に影響を及ぼす可能性はあるが、生体磁気信号が含まれる信号に対しても同様に重みを計算することができる。この場合、重みを固定する必要はなく、ある一定期間おきに重みを計算し、その間の区間は線形補間等をして時間的に重みを変化させても良い。

【0049】

また、生体磁気用センサ１１とリファレンスセンサ１２との位置や、角度の違いによって振幅や位相が異なるため、１つの生体磁気用センサ１１に対して複数のリファレンスセンサ１２を使用してもよい。

【0050】

例えば、ｉ番目の生体磁気用センサ１１のデータをｓ（ｉ）、リファレンスセンサ１２の数ｍとし、ｊ番目のリファレンスセンサ１２に含まれる５０Ｈｚのノイズデータをｒ_{（ｊ,５０）}（ｉ）、重み係数をｗ_ｊとすると、以下のようなに示すことができる。

【0051】

【数13】

【0052】

そして、情報処理装置３００は、式（１３）が最小になるように、重み係数ｗ_ｊを算出する。これにより、情報処理装置３００は、特定の周波数のノイズを除去することができる。なお、上記実施形態では、他の重みを計算する方法を適用することも可能である。

【0053】

〔４．情報処理装置の構成〕
次に、図８を用いて、実施形態に係る情報処理装置３００の構成について説明する。図８は、実施形態に係る情報処理装置３００の構成例を示す図である。図８に示すように、情報処理装置３００は、通信手段３１０と、記憶手段３２０と、制御手段３３０とを有する。

【0054】

（通信手段３１０について）
通信手段３１０は、例えば、ＮＩＣ（Network Interface Card）等によって実現される。そして、通信手段３１０は、インターネット等の所定のネットワークにより有線または無線で接続され、各種装置との間で情報の送受信を行う。

【0055】

（記憶手段３２０について）
記憶手段３２０は、例えば、ＲＡＭ（Random Access Memory)、フラッシュメモリ（Flash Memory）等の半導体メモリ素子、または、ハードディスク、光ディスク等の記憶装置によって実現される。

【0056】

（制御手段３３０について）
制御手段３３０は、コントローラ（Controller）であり、例えば、ＣＰＵ（Central Processing Unit）やＭＰＵ（Micro Processing Unit）等によって、情報処理装置３００内部の記憶装置に記憶されている各種プログラムがＲＡＭを作業領域として実行されることにより実現される。また、制御手段３３０は、コントローラであり、例えば、ＡＳＩＣ（Application Specific Integrated Circuit）やＦＰＧＡ（Field Programmable Gate Array）等の集積回路により実現される。

【0057】

図８に示すように、制御手段３３０は、取得手段３３１（物理量取得手段及び周波数取得手段の一例）と、ロックイン検波手段３３２と、計算手段３３３と、減算手段３３４とを有し、以下に説明する情報処理の機能や作用を実現または実行する。なお、制御手段３３０の内部構成は、図８に示した構成に限られず、後述する情報処理を行う構成であれば他の構成であってもよい。また、制御手段３３０が有する各処理部の接続関係は、図８に示した接続関係に限られず、他の接続関係であってもよい。

【0058】

（取得手段３３１について）
取得手段３３１は、各種情報を取得する。具体的には、取得手段３３１は、所定のセンサによって検知された時間的に変動する信号と雑音とを含む物理量を取得する。また、取得手段３３１は、物理量に含まれる雑音の周波数を取得する。

【0059】

また、複数の取得手段を有してもよい。具体的には、第１取得手段は、複数の所定のセンサのうち、少なくとも１つによって検知された信号と雑音とを含む物理量を取得する。また、第２取得手段は、複数の所定のセンサのうち、少なくとも１つによって検知された雑音を取得する。

【0060】

例えば、取得手段３３１は、生体磁気用センサ１１によって検知された生体磁気信号を取得する。また、取得手段３３１は、生体磁気用センサ１１によって検知された環境ノイズを取得する。取得手段３３１は、リファレンスセンサ１２によって検知された環境ノイズを取得する。

【0061】

（ロックイン検波手段３３２について）
ロックイン検波手段３３２は、周波数を参照周波数として、物理量に含まれる周波数の雑音の振幅と位相とを検出するロックイン検波を行う。また、ロックイン検波手段３３２は、周波数を参照周波数として、第２取得手段によって検出された雑音に含まれる周波数の振幅と位相とを検出するロックイン検波を行う。

【0062】

例えば、ロックイン検波手段３３２は、ロックイン検波の参照周波数を５０Ｈｚに設定する。そして、ロックイン検波手段３３２は、生体磁気用センサ１１によって検知された磁気信号に含まれる５０Ｈｚ成分の振幅と位相とを検出する。

【0063】

例えば、ロックイン検波手段３３２は、リファレンスセンサ１２によって検知された環境ノイズに対して、５０Ｈｚを参照信号としてロックイン検波を行う。また、ロックイン検波手段３３２は、リファレンスセンサ１２によって検知された環境ノイズに含まれる５０Ｈｚのノイズの振幅と位相とを検出し、得られた振幅と位相とを有する正弦波を生成する。

【0064】

（計算手段３３３について）
計算手段３３３は、重み係数を計算する。例えば、計算手段３３３は、生体磁気用センサ１１とリファレンスセンサ１２との環境ノイズがサンプル数ｎの離散データである場合に、ｉ番目の生体磁気用センサ１１のデータをｓ（ｉ）、リファレンスセンサ１２に含まれる５０Ｈｚのノイズデータをｒ_５０（ｉ）、重み係数をｗとすると、式（１２）のように示すことができる。そして、計算手段３３３は、式（１２）が最小になるように、重み係数ｗを算出する。

【0065】

（減算手段３３４について）
減算手段３３４は、物理量から振幅と位相とを有する正弦波を減算する。また、減算手段３３４は、第１取得手段によって検出された物理量から、位相と重み係数とを乗算した振幅を有する正弦波を減算する。例えば、減算手段３３４は、５０Ｈｚの正弦波を復元して元の磁気信号から差し引くことで５０Ｈｚのノイズのみを除去する。

【0066】

〔５．ハードウェア構成〕
また、上述してきた実施形態に係る情報処理装置３００等は、例えば、図７に示すような構成のコンピュータ１０００によって実現される。図７は、情報処理装置３００の機能を実現するコンピュータ１０００の一例を示すハードウェア構成図である。コンピュータ１０００は、ＣＰＵ１１００、ＲＡＭ１２００、ＲＯＭ（Read Only Memory）１３００、ＨＤＤ（Hard Disk Drive）１４００、通信インターフェイス（Ｉ／Ｆ）１５００、入出力インターフェイス（Ｉ／Ｆ）１６００、及びメディアインターフェイス（Ｉ／Ｆ）１７００を有する。

【0067】

ＣＰＵ１１００は、ＲＯＭ１３００又はＨＤＤ１４００に格納されたプログラムに基づいて動作し、各部の制御を行う。ＲＯＭ１３００は、コンピュータ１０００の起動時にＣＰＵ１１００によって実行されるブートプログラムや、コンピュータ１０００のハードウェアに依存するプログラム等を格納する。

【0068】

ＨＤＤ１４００は、ＣＰＵ１１００によって実行されるプログラム、及び、かかるプログラムによって使用されるデータ等を格納する。通信インターフェイス１５００は、ネットワークＮを介して他の機器からデータを受信してＣＰＵ１１００へ送り、ＣＰＵ１１００がネットワークＮを介して生成したデータを他の機器へ送信する。

【0069】

ＣＰＵ１１００は、入出力インターフェイス１６００を介して、ディスプレイやプリンタ等の出力装置、及び、キーボードやマウス等の入力装置を制御する。ＣＰＵ１１００は、入出力インターフェイス１６００を介して、入力装置からデータを取得する。また、ＣＰＵ１１００は、入出力インターフェイス１６００を介して生成したデータを出力装置へ出力する。

【0070】

メディアインターフェイス１７００は、記録媒体１８００に格納されたプログラム又はデータを読み取り、ＲＡＭ１２００を介してＣＰＵ１１００に提供する。ＣＰＵ１１００は、かかるプログラムを、メディアインターフェイス１７００を介して記録媒体１８００からＲＡＭ１２００上にロードし、ロードしたプログラムを実行する。記録媒体１８００は、例えばＤＶＤ（Digital Versatile Disc）、ＰＤ（Phase change rewritable Disk）等の光学記録媒体、ＭＯ（Magneto-Optical disk）等の光磁気記録媒体、テープ媒体、磁気記録媒体、または半導体メモリ等である。

【0071】

例えば、コンピュータ１０００が実施形態に係る情報処理装置３００として機能する場合、コンピュータ１０００のＣＰＵ１１００は、ＲＡＭ１２００上にロードされたプログラムを実行することにより、制御手段３３０の機能を実現する。また、ＨＤＤ１４００には、記憶手段３２０内のデータが格納される。コンピュータ１０００のＣＰＵ１１００は、これらのプログラムを記録媒体１８００から読み取って実行するが、他の例として、他の装置からネットワークＮを介してこれらのプログラムを取得してもよい。

【0072】

〔６．その他〕
また、上記実施形態及び変形例において説明した各処理のうち、自動的に行われるものとして説明した処理の全部または一部を手動的に行うこともでき、あるいは、手動的に行われるものとして説明した処理の全部または一部を公知の方法で自動的に行うこともできる。この他、上記文書中や図面中で示した処理手順、具体的名称、各種のデータやパラメータを含む情報については、特記する場合を除いて任意に変更することができる。例えば、各図に示した各種情報は、図示した情報に限られない。

【0073】

また、図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部または一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的または物理的に分散・統合して構成することができる。例えば、駆動回路１００と、収録回路２００と、情報処理装置３００とを統合して、情報処理装置としてもよい。

【0074】

また、上述してきた実施形態及び変形例は、処理内容を矛盾させない範囲で適宜組み合わせることが可能である。

【0075】

また、上述してきた「手段」は、「部（section、module、unit）」や、「工程」や、「回路」などに読み替えることができる。例えば、減算手段は、減算部や、減算工程や、減算手段や、減算回路に読み替えることができる。

【0076】

また、上述の本発明の実施形態は、発明の範囲を限定するものではなく、発明の範囲に含まれる一例に過ぎない。本発明のある実施形態は、上述の実施形態に対して、例えば、具体的な用途、構造、形状、作用、及び効果の少なくとも一部について、発明の要旨を逸脱しない範囲において変更、省略、及び追加がされたものであってもよい。

【符号の説明】

【0077】

１ 生体磁気計測システム
２ ＭＳＲ
１１ 生体磁気用センサ
１２ リファレンスセンサ
２０ クライオスタット
３０ ベッド
１００ 駆動回路
２００ 収録回路
３００ 情報処理装置
３１０ 通信手段
３２０ 記憶手段
３３０ 制御手段
３３１ 取得手段
３３２ ロックイン検波手段
３３３ 計算手段
３３４ 減算手段

【先行技術文献】

【非特許文献】

【0078】

【文献】Reduction of Non-periodic Environmental Magnetic Noise in MEG Measurement by Continuously Adjusted Least Squares Method Y. Adachi, M. Shimogawara, M. Higuchi, Y. Haruta and M. Ochiai IEEETRANSACTIONS ON APPLIED SUPERCONDUCTIVITY, VOL. I I, NO. I, MARCH 2001, DOI:10.1109/77.919433

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】