(57)【要約】      （修正有）【課題】ユーザにマーカーを付与することなく、ＭＲＩ画像における座標を、ＭＥＧ画像における座標に変換する座標変換装置を提供する。【解決手段】座標変換装置１００は、磁場強度を示す信号を出力する複数のセンサ１７０と、複数のセンサ１７０各々から出力された信号に基づいてＭＲＩ画像を生成する第１生成部１４１と、複数のセンサ１７０各々から出力された信号に基づいてＭＥＧ画像を生成する第２生成部１４２と、複数のＭＲＩ画像各々の同じＭＲＩ座標における磁場強度値に基づいて、ＭＲＩ座標をＭＥＧ画像におけるＭＥＧ座標に変換する変換部１４３と、を備える。【選択図】図１
2019134904.tif 000002
【特許請求の範囲】
【請求項１】
磁場強度を示す信号を出力する複数のセンサと、
前記複数のセンサ各々から出力された前記信号に基づいてＭＲＩ画像を生成する第１生成部と、
前記複数のセンサ各々から出力された前記信号に基づいてＭＥＧ画像を生成する第２生成部と、
前記複数のＭＲＩ画像各々の同じＭＲＩ座標における磁場強度値に基づいて、前記ＭＲＩ座標を前記ＭＥＧ画像におけるＭＥＧ座標に変換する変換部と、
を備える座標変換装置。
【請求項２】
前記ＭＲＩ座標の指定を受け付ける入力部を更に備え、
前記変換部は、前記入力部が受け付けたＭＲＩ座標をＭＥＧ座標に変換する
ことを特徴とする請求項１に記載の座標変換装置。
【請求項３】
前記変換部は、磁場ダイポールから前記ＭＲＩ座標までの位置ベクトルと、前記ＭＲＩ座標における磁場強度と、前記磁場ダイポールのモーメントに基づいて、前記ＭＲＩ座標を前記ＭＥＧ座標に変換する
ことを特徴とする請求項１又は２に記載の座標変換装置。
【請求項４】
前記変換部は、磁場強度をＢ、前記磁場ダイポールから測定座標までの位置ベクトルをｒ、前記磁場ダイポールのモーメントをｍとしたとき、下記式（１）に基づいて前記ＭＥＧ座標を算出することを特徴とする請求項３に記載の座標変換装置。
【数１】
2019134904.tif 000005
・・・式（１）
【請求項５】
前記変換部は、前記ＭＲＩ画像における特定の座標における前記複数のセンサそれぞれの磁場強度を前記式（１）に入力して連立方程式を解くことで、前記ＭＥＧ画像における座標を算出する
ことを特徴とする請求項４に記載の座標変換装置。
【請求項６】
前記第１生成部が生成したＭＲＩ画像の歪みを補正する補正部を更に備える
ことを特徴とする請求項１〜５のいずれか一項に記載の座標変換装置。
【請求項７】
コンピュータが、
複数のセンサ各々から出力された磁場強度を示す信号に基づいてＭＲＩ画像を生成する第１生成ステップと、
前記複数のセンサ各々から出力された前記信号に基づいてＭＥＧ画像を生成する第２生成ステップと、
前記複数のＭＲＩ画像各々の同じＭＲＩ座標における磁場強度値に基づいて、前記ＭＲＩ座標を前記ＭＥＧ画像におけるＭＥＧ座標に変換する変換ステップと、
を実行する座標変換方法。
【請求項８】
コンピュータに、
複数のセンサ各々から出力された磁場強度を示す信号に基づいてＭＲＩ画像を生成する第１生成ステップと、
前記複数のセンサ各々から出力された前記信号に基づいてＭＥＧ画像を生成する第２生成ステップと、
前記複数のＭＲＩ画像各々の同じＭＲＩ座標における磁場強度値に基づいて、前記ＭＲＩ座標を前記ＭＥＧ画像におけるＭＥＧ座標に変換する変換ステップと、
を実現させる座標変換プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、ＭＲＩと脳磁計との間の座標変換を行う座標変換装置に関する。
【背景技術】
【０００２】
従来、生体内の物理的情報や、分子的、科学的情報を取得する手法の一つとして、ＭＲＩ（Magnetic Resonance Imaging）が利用されている。また、脳内の情報を取得する手法の一つとして、脳磁計を利用することも行われている。そして、ＭＲＩと脳磁計とを組み合わせることも考えられている（例えば、特許文献１）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００５−２７０４２２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところで、特許文献１において、ＭＲＩと脳磁計（ＭＥＧ）とを組み合わせる場合には、それぞれの測定座標系が互いに異なることから、得られるデータに対して、座標合わせをする必要がある。具体的には、測定対象のユーザに多数のマーカーを付与して、測定を行うことで、座標合わせを行っている。即ち、ＭＲＩによる測定を行う際には、ユーザの所定の場所（例えば、こめかみ、額、など）にビタミン剤などの水状のものを塗布して測定を行い、その後に、水状のものを塗布した同じ個所に、マーカーコイルを貼り付けてＭＥＧによる測定を行う必要がある。しかしながら、当該手法では、ユーザに多数のマーカーを付与する手間が煩雑であるという問題がある。また、マーカーの付与は、人間の手に依るものであるところ、どうしても貼付位置の精度に誤差が出るという問題があり、ＭＲＩにおける座標と、ＭＥＧにおける座標との間に誤差が出るという問題もあった。
【０００５】
そこで、本発明は、上記問題に鑑みてなされたものであり、マーカーコイル等をユーザに付与せずともＭＲＩ座標をＭＥＧ座標に変換できる座標変換装置を提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明の一実施態様に係る座標変換装置は、磁場強度を示す信号を出力する複数のセンサと、複数のセンサ各々から出力された信号に基づいて複数のＭＲＩ画像を生成する第１生成部と、複数のセンサ各々から出力された信号に基づいてＭＥＧ画像を生成する第２生成部と、複数のＭＲＩ画像各々の同じＭＲＩ座標における磁場強度値に基づいて、ＭＲＩ座標をＭＥＧ画像におけるＭＥＧ座標に変換する変換部と、を備える。
【０００７】
また、上記座標変換装置において、ＭＲＩ座標の指定を受け付ける入力部を更に備え、変換部は、入力部が受け付けたＭＲＩ座標をＭＥＧ座標に変換することとしてもよい。
【０００８】
また、上記座標変換装置において、変換部は、磁場ダイポールからＭＲＩ座標までの位置ベクトルと、ＭＲＩ座標における磁場強度と、磁場ダイポールのモーメントに基づいて、ＭＲＩ座標をＭＥＧ座標に変換することとしてもよい。
【０００９】
また、上記座標変換装置において、変換部は、磁場強度をＢ、磁場ダイポールから測定座標までの位置ベクトルをｒ、磁場ダイポールのモーメントをｍとしたとき、下記式（１）に基づいてＭＥＧ座標を算出することとしてもよい。
【００１０】
【数１】
2019134904.tif 000003
・・・式（１）
【００１１】
また、上記座標変換装置において、変換部は、ＭＥＧ画像における特定の座標における複数のセンサそれぞれの磁場強度を式（１）に入力して連立方程式を解くことで、ＭＥＧ画像における座標を算出することとしてもよい。
【００１２】
また、上記座標変換装置において、第１生成部が生成したＭＲＩ画像の歪みを補正する補正部を更に備えることとしてもよい。
【００１３】
本発明の一実施態様に係る座標変換方法は、コンピュータが、複数のセンサ各々から出力された磁場強度を示す信号に基づいてＭＲＩ画像を生成する第１生成ステップと、複数のセンサ各々から出力された信号に基づいてＭＥＧ画像を生成する第２生成ステップと、複数のＭＲＩ画像各々の同じＭＲＩ座標における磁場強度値に基づいて、ＭＲＩ座標をＭＥＧ画像におけるＭＥＧ座標に変換する変換ステップと、を実行する。
【００１４】
本発明の一実施態様に係る座標変換プログラムは、コンピュータに、複数のセンサ各々から出力された磁場強度を示す信号に基づいてＭＲＩ画像を生成する第１生成ステップと、複数のセンサ各々から出力された信号に基づいてＭＥＧ画像を生成する第２生成ステップと、複数のＭＲＩ画像各々の同じＭＲＩ座標における磁場強度値に基づいて、ＭＲＩ座標をＭＥＧ画像におけるＭＥＧ座標に変換する変換ステップと、を実現させる。
【発明の効果】
【００１５】
本発明に係る座標変換装置、座標変換方法、座標変換プログラムによれば、マーカーコイルなどをユーザに付与せずともＭＲＩ座標をＭＥＧ座標に変換することができる。したがって、ＭＲＩ画像と、ＭＥＧ画像との位置合わせを容易に実現することができる。
【図面の簡単な説明】
【００１６】
【図１】座標変換装置の構成例を示すブロック図である。
【図２】（ａ）は、ＭＲＩ画像の一例であり、（ｂ）はＭＥＧ画像の一例である。
【図３】（ａ）〜（ｄ）は、脳磁計センサの位置によって得られるＭＲＩ画像が異なることを示す図である。
【図４】座標変換装置の動作例を示すフローチャートである。
【図５】ＭＲＩ画像の一例を示しており、勾配磁場歪み等の影響を受けずに得られた画像の例を示している。
【図６】（ａ）は、ＭＲＩ画像の一例を示しており、勾配磁場歪み等の影響を受けて歪んだ画像が得られた場合の例を示している。（ｂ）は、（ａ）のａｘ線、ａｙ線における信号強度の変化を示すグラフである。
【図７】（ａ）は、図６のＭＲＩ画像を補正したＭＲＩ画像を示している。（ｂ）は、（ａ）のｂｘ線、ｂｙ線における信号強度の変化を示すグラフである。
【発明を実施するための形態】
【００１７】
本発明の一実施形態について、図面を参照して説明する。
【００１８】
＜実施の形態＞
＜構成＞
以下、実施の形態に係る座標変換装置について、図面を参照しながら、詳細に説明する。
【００１９】
図１は、座標変換装置１００の機能構成例を示すブロック図である。図１に示すように、座標変換装置１００は、通信部１１０と、入力部１２０と、記憶部１３０と、制御部１４０と、出力部１５０と、を備える。座標変換装置１００は、ＭＲＩ画像とＭＥＧ画像との間で、ＭＲＩ画像における所定の座標値を、ＭＥＧ座標系における座標値に変換する機能を有するコンピュータシステムである。座標変換装置１００は、各脳磁計センサ１７０を含んで構成されてもよい。図１においては、座標変換装置１００と各脳磁計センサ１７０とを接続線１６０で接続した例を示している。なお、接続線１６０は、脳磁計センサ１７０の数の分だけ存在する。
【００２０】
脳磁計センサ１７０は、磁場を検出するセンサであり、例えば、ＳＱＵＩＤセンサにより実現することができる。脳磁計センサ１７０は、検出した磁場の強度を示す信号を通信部１１０に送信する。
【００２１】
通信部１１０は、各脳磁計センサ１７０と通信を実行する機能を有する通信インターフェースである。通信部１１０は、各脳磁計センサ１７０から伝達された磁場信号を、制御部１４０に伝達する。通信部１１０は、各脳磁計センサ１７０と通信ができれば、有線通信、無線通信のいずれで通信を行ってもよく、通信のための通信プロトコルも任意である。
【００２２】
入力部１２０は、ユーザからの入力を受け付けて、制御部１４０に伝達する機能を有する。入力部１２０は、例えば、タッチパネルやタッチパッドなどソフトウェア入力可能な機器により実現されてもよいし、マウスやキーボード等のハードウェアにより実現されてもよいし、音声入力を受け付けるためのマイクロフォンにより実現されてもよい。入力部１２０は、ユーザから、ＭＲＩ画像のＭＲＩ座標の指定を受け付けて、制御部１４０に伝達する。
【００２３】
記憶部１３０は、座標変換装置１００が動作上必要とする各種データやプログラムを記憶する機能を有する記録媒体である。記憶部１３０は、例えば、ＨＤＤ（Hard Disc Drive）、ＳＳＤ（Solid State Drive）、フラッシュメモリなどにより実現することができるが、これらに限定されるものではない。記憶部１３０は、各脳磁計センサ１７０の配置位置（三次元空間における座標値）を示す配置情報を記憶している。また、記憶部１３０は、磁場強度と磁場ダイポールの関係式を記憶しており、当該関係式を用いて、ＭＲＩ座標をＭＥＧ座標に変換する変換プログラムを記憶している。
【００２４】
制御部１４０は、座標変換装置１００の各部を制御する機能を有するプロセッサである。制御部１４０は、従来のＭＲＩと同様に脳磁計センサ１７０からのセンシングデータを用いて、ＭＲＩ画像を生成する第１生成部１４１として機能する。なお、制御部１４０は、一つの脳磁計センサ１７０からのセンシングデータのみでＭＲＩ画像を生成することができるが、複数のセンシングデータを利用すれば、より鮮明なＭＲＩ画像を生成することができる。また、制御部１４０は、従来の脳磁計と同様に各脳磁計センサ１７０からのセンシングデータを用いて、ＭＥＧ画像を生成する第２生成部１４２として機能する。制御部１４０は、入力部１２０からＭＲＩ画像におけるＭＲＩ座標（画像上の座標系における座標値）の指定を受け付けると、ＭＥＧ座標系における座標（ＭＥＧ画像における座標値であって、実際の３次元空間上の座標値、即ち、脳磁計センサ１７０の配置空間を基準とする座標値）に変換する変換部１４３として機能する。制御部１４０は、記憶部１３０に記憶されている変換プログラムを読み出して実行することにより、ＭＲＩ座標をＭＥＧ座標に変換する。制御部１４０は、記憶部１３０に記憶されている磁場強度と磁場ダイポールの関係式を利用して、ＭＲＩ座標をＭＥＧ座標に変換する。
【００２５】
以下、制御部１４０による座標変換方法について詳細に説明する。
【００２６】
ＭＲＩ画像については、一つのセンサから得られる情報に基づいて生成できることが知られている。
【００２７】
一方、ＭＥＧ画像については、複数のセンサから得られる情報に基づいて、脳内の磁場分布を示すＭＥＧ画像を生成する。
【００２８】
即ち、座標変換装置１００は、ＭＲＩとして、例えば、図２（ａ）に示すようなＭＲＩ画像２０１を生成する。一方で、座標変換装置１００は、脳磁計として、例えば、図２（ｂ）に示すようなＭＥＧ画像を生成する。なお、図２（ａ）においては、ユーザ１０に対して向けられる脳磁計センサ１７０（１７０ａ、１７０ｂ、１７０ｃ、１７０ｄ）の配置例も示している。
【００２９】
ＭＲＩ画像の生成や、ＭＥＧ画像の生成は、従来と同様の手法で行うものとし、ここでは、その詳細は割愛する。
【００３０】
ここで、ＭＲＩ画像は、一つの脳磁計センサ１７０のセンシングデータから、周波数エンコーディング及び位相エンコーディングを行うことで生成できることが知られているが、その脳磁計センサ１７０の配置位置によって得られるＭＲＩ画像は互いに異なってくる。
【００３１】
即ち、図３（ａ）〜（ｄ）各々の右図に示すように、ＭＲＩ画像における濃淡となって表れる。図３（ａ）には、脳磁計センサ１７０ｅからのセンシングデータを用いて生成されたＭＲＩ画像３０１を、図３（ｂ）には、脳磁計センサ１７０ｆからのセンシングデータを用いて生成されたＭＲＩ画像３０２を、図３（ｃ）には、脳磁計センサ１７０ｇからのセンシングデータを用いて生成されたＭＲＩ画像３０３を、図３（ｄ）には、脳磁計センサ１７０ｈからのセンシングデータを用いて生成されたＭＲＩ画像３０４を示している。各ＭＲＩ画像から理解できるように、各センサからの距離が近い箇所ほど鮮明な画像が得られる。座標変換装置１００は、これらの画像を重畳することで、例えば、図２（ａ）に示すような、ＭＲＩ画像２０１を生成することができる。
【００３２】
ところで、ＭＲＩ画像においては、脳の位置を示す座標は、ＭＲＩ画像の２次元座標でしか得ることができない。即ち、ユーザ１０の脳の実際の３次元座標にはなっていない。そのため、ＭＲＩ画像における座標は、ＭＥＧ画像における正確な位置を得ることができない。
【００３３】
そこで、ＭＲＩ画像における所定の座標が、ＭＥＧ画像においては、どのような座標になるかを知りたい場合には、座標変換装置１００は、磁場強度と磁場ダイポールの関係式を利用することで、ＭＲＩ座標をＭＥＧ座標に変換する。より正確には、入力部１０２がＭＲＩ画像における所定のＭＲＩ座標（ｉｘ、ｉｙ）の指定を受け付けると、変換部１４３は、受け付けたＭＲＩ座標における各脳磁計センサ１７０が検出した磁場強度を取得する。そして、変換部１４３は、取得した磁場強度を、磁場強度と磁場ダイポールの関係式に代入して複数の関係式を得る。磁場強度と磁場ダイポールの関係式を、以下の式（１）に示す。
【００３４】
【数２】
2019134904.tif 000004
・・・式（１）
【００３５】
上記式において、Ｂは、磁場強度を示し、ｒは、磁場ダイポールの中心から測定座標までの位置ベクトルを示し、ｍは、磁場ダイポールのモーメントを示している。ここで、ＭＲＩ座標（ｉｘ、ｉｙ）に相当する測定座標（算出対象のＭＥＧ座標）が、未知の値であることから、少なくとも３つの関係式が得られれば、ＭＥＧ座標（ｘ、ｙ、ｚ）を算出することができる。変換部１４３は、式（１）を利用して複数の関係式を生成し、生成した複数の関係式を連立方程式として解くことで、ＭＥＧ座標（ｘ、ｙ、ｚ）を算出する。例えば、脳磁計センサ１７０ａで測定したＭＲＩ画像における座標Ａにおける磁場強度を上記式（１）に代入して第１の関係式を生成し、脳磁計センサ１７０ｂで測定したＭＲＩ画像における座標Ａにおける磁場強度を上記式（１）に代入して第２の関係式を生成し、脳磁計センサ１７０ｃで測定したＭＲＩ画像における座標Ａにおける磁場強度を上記式（１）に代入して第３の関係式を生成することで、少なくとも３つの関係式を得る。関係式が多いほど、算出されるＭＥＧ座標の確度が向上する。即ち、変換部１４３は、ＭＲＩ座標（ｉｘ、ｉｙ）を、ＭＥＧ座標（ｘ、ｙ、ｚ）に変換できたことになる。
【００３６】
出力部１５０は、変換部１４３が算出したＭＥＧ座標を示す情報（例えば、座標を示す数値情報でもよいし、座標を示す音声情報でもよいし、指定された座標をＭＥＧ画像に示した画像情報でもよい）を出力する。出力部１５０は、例えば、モニターやスピーカなどにより実現することができる。
【００３７】
以上が、座標変換装置１００の構成である。
【００３８】
＜動作＞
図４は、座標変換装置１００の動作を示すフローチャートである。当該フローチャートは、ＭＲＩ座標をＭＥＧ座標に変換する処理を示すものである。
【００３９】
座標変換装置１００の通信部１１０は、各脳磁計センサ１７０から、それぞれのセンシングデータを取得する（ステップＳ４０１）。通信部１１０は、取得したセンシングデータを制御部１４０に伝達する。
【００４０】
制御部１４０の第１生成部１４１は、伝達されたセンシングデータそれぞれからＭＲＩ画像を生成し、これを重畳して１枚のＭＲＩ画像を生成する（ステップＳ４０２）。なお、当該処理は、一つのセンサからのセンシングデータに基づいて１枚のＭＲＩ画像を生成するものであってもよい。
【００４１】
一方、制御部１４０の第２生成部１４２は、伝達されたセンシングデータそれぞれからＭＥＧ画像を生成する（ステップＳ４０３）。なお、ステップＳ４０２とステップＳ４０３との処理は、互いに並列に実行されてもよいし、ステップＳ４０３の処理の方がステップＳ４０２の処理よりも先に処理されてもよい。
【００４２】
次に、入力部１０２は、ユーザからＭＲＩ画像におけるＭＲＩ座標の指定を受け付ける。ＭＲＩ座標は、ユーザがＭＥＧ画像上においてどのような座標になるのかを知りたい座標である。入力部１０２は、入力されたＭＲＩ座標を制御部１４０に伝達する。
【００４３】
変換部１４３は、伝達されたＭＲＩ座標における磁場強度を、各脳磁計センサ１７０が検出したセンシングデータを参照して、特定する。ここで、変換部１４３は、複数の脳磁計センサ１７０が検出したセンシングデータを参照して、複数の磁場強度を特定する。
【００４４】
変換部１４３は、特定した磁場強度を、磁場強度と磁場ダイポールの関係式（上記式（１））に代入する（ステップＳ４０６）。これにより、一つの磁場強度から、一つの関係式が得られる。したがって、複数の磁場強度を特定することで、複数の関係式を得ることができる。
【００４５】
変換部１４３は、複数の関係式を作成する。そして、作成した複数の関係式を連立方程式にみたてて、式（１）における未知数ｒを解く（ステップＳ４０７）。これにより、変換部１４３は、ＭＥＧ座標を算出する、即ち、ＭＲＩ座標をＭＥＧ座標に変換することができる。
【００４６】
変換部１４３は、算出したＭＥＧ座標を、出力部１５０に出力させて（ステップＳ４０８）、処理を終了する。
【００４７】
以上が、座標変換装置１００によるＭＲＩ座標をＭＥＧ座標に変換する処理動作である。
【００４８】
＜まとめ＞
座標変換装置１００は、ユーザが特定したいＭＥＧ座標に対応するＭＲＩ座標を、ＭＥＧ画像を生成するために利用した各脳磁計センサ１７０のセンシングデータの当該ＭＲＩ座標の磁場強度を用いて、磁場強度と磁場ダイポールの関係式に代入し、複数の関係式を得ることができる。その結果、複数の関係式を連立方程式としてとくことで、ＭＲＩ座標に対応するＭＥＧ座標を得ることができる。したがって、従来のように、マーカーコイル等をユーザに付与せずとも、ＭＲＩ座標をＭＥＧ座標に変換することができる。よって、ＭＲＩや脳磁におけるユーザの処理負荷を軽減することができる。
【００４９】
＜補足＞
本発明に係る座標変換装置は、上記実施の形態に限定されるものではなく、他の手法により実現されてもよいことは言うまでもない。以下、各種変形例について説明する。
【００５０】
（１）上記実施の形態において、座標変換装置１００は、脳磁計センサ１７０からセンシングデータを取得して、ＭＲＩ画像、ＭＥＧ画像を生成することとしたが、座標変換装置１００は、単純に、ＭＲＩ座標をＭＥＧ座標に変換するだけの機能を有するものであってもよい。即ち、通信部１１０を備えず、予め、記憶部１３０に各脳磁計センサ１７０が出力したセンシングデータおよび当該センシングデータを用いて生成されたＭＲＩ画像の情報及び、座標変換プログラムを記憶し、制御部１４０がＭＲＩ座標をＭＥＧ座標に変換することとしてもよい。
【００５１】
（２）上記実施の形態において、座標変換装置１００は、ＭＲＩ画像を補正する補正部を備えてもよい。第１生成部が生成するＭＲＩ画像は、勾配磁場歪みの影響等により、歪みが発生する可能性がある。図５〜７を用いて具体的に説明する。
【００５２】
図５は、ＭＲＩ画像の一例を示しており、勾配磁場歪み等の影響を受けずに得られた画像の例を示している。何も影響を受けずに得られたＭＲＩ画像は、歪みが発生しない。
【００５３】
これに対して、勾配磁場歪み等の影響を受けると、図６（ａ）に示すように、図５と比較して、歪んだＭＲＩ画像が得られる場合がある。図６（ｂ）は、図６（ａ）のａｘ及びａｙにおける磁場強度の推移を示すグラフである。
【００５４】
補正部は、図６（ａ）に示すような歪んだＭＲＩ画像を得た場合には、ＭＲＩ画像の各ピクセル位置の磁気ダイポール推定を行う。そして、その散布図から、４０ピクセル×４０ピクセルのメッシュにデータ補完を行う。これにより、補正部は、図６（ａ）に示すＭＲＩ画像と、当該ＭＲＩ画像を得た際に利用したセンシングデータを利用することで、図７（ａ）に示すような歪みのないＭＲＩ画像に補正することができる。具体的に説明するとＭＲＩ画像のピクセル番号には絶対的な位置情報が含まれておらず、順序情報のみが含まれている。補正部は、上記実施の形態において示した手法により、各ピクセル番号に対応する位置情報を推定する。そして、推定によって得られた位置情報を用いて各ピクセルをｘｙ平面上に再投影することで、散布図を生成する。このとき、補正部が得る位置情報はＭＥＧのセンサの配置座標で定まるものであって、ＭＲＩの勾配磁場によって定まるものではない。したがって、勾配磁場歪みの影響を受けないことになる。補正部は、生成した散布図から適当なメッシュサイズ、例えば、４０×４０の等間隔メッシュの各点の強度をデータ補完（内挿）、例えば、スプライン補完や線形補完を行って、図７（ａ）に示すような補正後の画像を得ることができる。
【００５５】
図７（ａ）は、図６（ａ）に示すＭＲＩ画像を補正した例を示している。図７（ａ）を見れば理解できるように、図６（ａ）のＭＲＩ画像における歪みが補正されている。そして、図７（ｂ）は、図７（ａ）のＭＲＩ画像におけるｂｘ及びｂｙにおける強度の推移を示すグラフである。このように、補正部は、歪んだＭＲＩ画像を補正することもできる。
【００５６】
（３）上記実施の形態において、座標変換装置１００に係る各装置の各機能部は、集積回路（ＩＣ（Integrated Circuit）チップ、ＬＳＩ（Large Scale Integration））等に形成された論理回路（ハードウェア）や専用回路によって実現してもよいし、ＣＰＵ（Central Processing Unit）およびメモリを用いてソフトウェアによって実現してもよい。また、各機能部は、１又は複数の集積回路により実現されてよく、複数の機能部の機能を１つの集積回路により実現されることとしてもよい。ＬＳＩは、集積度の違いにより、ＶＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩなどと呼称されることもある。また、当該回路は、再構築可能な回路（例えば、ＦＰＧＡ：Field Programmable Gate Away）により実現されてもよい。
【００５７】
座標変換装置１００の各機能部をソフトウェアにより実現する場合、座標変換装置１００は、各機能を実現するソフトウェアである制御プログラム（ＭＲＩ画像生成プログラム、ＭＥＧ画像生成プログラム、座標変換プログラム等）の命令を実行するＣＰＵ、制御プログラムおよび各種データがコンピュータ（又はＣＰＵ）で読み取り可能に記録されたＲＯＭ（Read Only Memory）又は記憶装置（これらを「記録媒体」と称する）、上記制御プログラムを展開するＲＡＭ（Random Access Memory）などを備えている。そして、コンピュータ（又はＣＰＵ）が上記制御プログラムを上記記録媒体から読み取って実行することにより、本発明の目的が達成される。上記記録媒体としては、「一時的でない有形の媒体」、例えば、テープ、ディスク、カード、半導体メモリ、プログラマブルな論理回路などを用いることができる。また、上記制御プログラムは、当該制御プログラムを伝送可能な任意の伝送媒体（通信ネットワークや放送波等）を介して上記コンピュータに供給されてもよい。本発明は、上記委任状送信プログラムが電子的な伝送によって具現化された、搬送波に埋め込まれたデータ信号の形態でも実現され得る。
【００５８】
なお、上記制御プログラムは、例えば、ActionScript、JavaScript（登録商標）などのスクリプト言語、Objective-C、Java（登録商標）などのオブジェクト指向プログラミング言語、HTML5などのマークアップ言語などを用いて実装できる。
【００５９】
（４）本発明を諸図面や実施形態に基づき説明してきたが、当業者であれば本開示に基づき種々の変形や修正を行うことが容易であることに注意されたい。従って、これらの変形や修正は本発明の範囲に含まれることに留意されたい。例えば、各手段、各ステップ等に含まれる機能等は論理的に矛盾しないように再配置可能であり、複数の手段やステップ等を１つに組み合わせたり、或いは分割したりすることが可能である。また、上記実施の形態に示す構成を適宜組み合わせてもよい。
【符号の説明】
【００６０】
１００  座標変換装置
１１０  通信部
１２０  入力部
１３０  記憶部
１４０  制御部
１４１  第１生成部
１４２  第２生成部
１４３  変換部
１５０  出力部
１６０  接続線
１７０  脳磁計センサ
【図１】
2019134904.tif 000006
【図２】
2019134904.tif 000007
【図３】
2019134904.tif 000008
【図４】
2019134904.tif 000009
【図５】
2019134904.tif 000010
【図６】
2019134904.tif 000011
【図７】
2019134904.tif 000012