JP2008307228A

JP2008307228A - 画像生成方法、位置補正方法および磁気共鳴イメージング装置

Info

Publication number: JP2008307228A
Application number: JP2007158014A
Authority: JP
Inventors: Hitoshi Ikeda; 仁池田
Original assignee: GE Medical Systems Global Technology Co LLC
Current assignee: GE Medical Systems Global Technology Co LLC
Priority date: 2007-06-14
Filing date: 2007-06-14
Publication date: 2008-12-25
Anticipated expiration: 2027-06-14
Also published as: US7791341B2; US20080309334A1; JP5355866B2

Abstract

【課題】複数の磁気共鳴信号における信号強度が最大となる部位をｋ空間における周波数軸方向の中心Ｏに変位させるための変位量に基づいて、センターリング処理を実施することにより、画像品質を向上させる。
【解決手段】複数の周波数軸方向処理データにおける信号強度の最大値を算出する。そして、全ての周波数軸方向処理データにおける信号強度の最大値のうち最大値の一番大きい周波数軸方向処理データを決定し、決定した周波数軸方向処理データにおける信号強度が最大となる点がｋ空間における周波数軸方向の中心Ｏに変位されるための変位量を算出する。そして、算出された変位量に基づいて、ｋ空間において、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、ハーフエコー処理を実施して画像を生成する。
【選択図】図６

Description

本発明は、磁気共鳴画像生成方法、磁気共鳴信号の位置補正方法および磁気共鳴イメージング装置に関する。

磁気共鳴イメージング装置は、核磁気共鳴現象を利用して磁気共鳴信号を発生させ、被検体の断層画像を撮影する装置である。

磁気共鳴イメージング装置は、撮像時間に時間を要するため、様々な方法により撮像時間の短縮が図られている。

例えば、エコー時間に対するエコー波形の対称性および周波数方向におけるｋ空間データの対称性を利用して、発生するエコー信号の半分よりやや多いデータを収集し、残りの部分を共役対称性を利用して算出するハーフエコー（half echo）法や、感度分布の異なる複数のＲＦコイルからなるフェーズドアレイコイルを用いて、位相エンコードを間引いたシーケンスを実行し、行列演算によって折り返しアーチファクトを除去する展開処理を行うことにより撮影時間を短縮化するパラレルイメージング（parallel imaging）法などが知られている（例えば、特許文献１参照。）。

ハーフエコー法において、欠落しているデータを埋め合わせるためにホモダイン処理を行っている。ホモダイン処理を行うためには、エコー信号の中心をハイパスフィルター（High Pass Filter）およびローパスフィルター（Low Pass Filter）に通して処理を行う。このハイパスフィルターおよびローパスフィルターはｋ空間の中心にあるため、ｋ空間の中心からずれているエコー信号の中心をｋ空間の中心に変位させる必要がある。通常エコー信号のピークがエコー信号の中心であるため、エコー信号のピークをｋ空間の中心に変位させる。

しかし、回転磁場（Ｂ１）不均一などの影響により、フェーズドアレイコイルにより受信した複数のエコー信号におけるピークは必ずしも一致しない。例えば、図２９に示すように、フェーズドアレイコイルを構成する複数のＲＦコイルそれぞれから受信した複数のエコー信号が、エコー信号Ｅ０１、Ｅ０２およびＥ０３であるとすると、それぞれのエコー信号のピークは、ｋ空間の中心Ｏから周波数Ｆ_Ｅ０１、Ｆ_Ｅ０２およびＦ_Ｅ０３と異なった周波数だけずれている。そのため、受信したエコー信号それぞれのピークをｋ空間における周波数軸方向の中心Ｏに変位させるためには、それぞれのエコー信号を異なった変位量で変位させる必要がある。

また、パラレルイメージング法においては、ｋ空間において、全てのエコー信号が同じ変位量で変位されていないと、キャリブレーションスキャンと本スキャンにおけるコイルの感度マップの位相関係が乱れ、アーチファクトが生じる原因となる。そのため、ハーフエコー法とパラレルイメージング法を併用する場合、基準となるチャンネルのＲＦコイルを選択し、この選択されたＲＦコイルが受信したエコー信号のピークをｋ空間における周波数軸方向の中心Ｏに変位させるための変位量を、全てのエコー信号に適用することにより全てのエコー信号に対してセンターリング処理を実施している。例えば、図３０に示すように、基準と成るチャンネルのＲＦコイルが受信したエコー信号をエコー信号Ｅ０１とすると、周波数ＦＥ０１をエコー信号Ｅ０１のピークをｋ空間における周波数軸方向の中心Ｏに変位させるための変位量として、エコー信号Ｅ０１、Ｅ０２およびＥ０３を変位させ、センターリング処理を実施する。
特開２００５−１９８７１５号公報（例えば、段落０００７から段落００１０参照。）

ここで、パラレルイメージング法による撮影が、頭などの局所的な部位から腹等の広範囲の部位まで及んできたため、フェーズドアレイコイルを構成する複数のＲＦコイルが多チャンネル化されている。そのため、フェーズドアレイコイルを構成する複数のＲＦコイルのそれぞれが受信するそれぞれのエコー信号における信号強度の差が、顕著になってきている。

図３１に示すように、信号強度が低いエコー信号は、ノイズなどの影響がある場合、エコー信号の信号強度の最大値がノイズ部分となり、エコー信号の信号強度の最大値とピークＰとが一致しない場合がある。したがって、ノイズなどの影響がある場合、信号強度の低いエコー信号を受信するＲＦコイルを基準コイルとして選択した場合、エコー信号のノイズ部分をエコー信号の最大値と判断してしまう。そのため、このエコー信号に基づいて、全てのエコー信号をセンターリング処理させるための変位量を算出し、この変位量に基づいて、全てのエコー信号を変位させると、全てのエコー信号において、ｋ空間における周波数軸方向の中心Ｏにエコー信号の実際のピークＰが変位されず、実際のピークＰからずれた部分がｋ空間における周波数軸方向の中心Ｏに変位されることとなる。そのため、生成される画像にアーチファクトが生じることがある。

したがって、本発明は、ハーフエコー法とパラレルイメージング法を併用した際に、生成した画像にアーチファクトが生じないようなエコー信号の位置補正方法、画像生成方法および磁気イメージング装置を提供する。

本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、前記複数の磁気共鳴信号のプロファイルに基づいて、当該複数の磁気共鳴信号から基準となる基準磁気共鳴信号を選択する第１ステップと、前記第１ステップにおいて選択された前記基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数との差である周波数差を算出する第２ステップと、前記第２ステップにおいて算出された前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップと、前記第３ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第４ステップとを含む。

好適には、前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて信号強度の最大値を算出し、当該最大値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する。

好適には、前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて積分値を算出し、当該積分値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する。

好適には、前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて当該複数の磁気共鳴信号同士の相関係数を算出し、前記複数の磁気共鳴信号ごとに算出された前記相関係数の平均値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する。

本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数と前記ｋ空間の中心における周波数との差である周波数差を算出し、当該周波数差の平均値を算出する第１ステップと、前記第１ステップにおいて算出された前記周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと、前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップとを含む。

本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と前記前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性を算出し、前記周波数差と前記プロファイル特性の積の平均値を算出する第１ステップと、前記第１ステップにおいて算出された前記周波数差と前記プロファイル特性の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと、前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップとを含む。

好適には、前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性は、前記複数の磁気共鳴信号のそれぞれにおける積分値である。

好適には、前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性は、前記複数の磁気共鳴信号のそれぞれにおける信号強度の最大値である。

本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出する第１ステップと、前記第１ステップにおいて算出されたそれぞれの前記周波数差の中央値を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと、前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップとを含む。

本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の標準偏差を算出する第１ステップと、前記第１ステップにおいて算出された前記標準偏差に基づいて、前記複数の磁気共鳴信号のそれぞれについて算出された前記周波数差から基準となる複数の基準周波数差を選択し、選択された当該複数の基準周波数差の平均値を算出する第２ステップと、前記第２ステップにおいて算出された前記複数の基準周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップと、前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップとを含む。

本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、前記複数の磁気共鳴信号のプロファイルに基づいて、当該複数の磁気共鳴信号から基準となる基準磁気共鳴信号を選択する第１ステップと、前記第１ステップにおいて選択された前記基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数との差である周波数差を算出する第２ステップと、前記第２ステップにおいて算出された前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップとを含む。

本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数と前記ｋ空間の中心における周波数との差である周波数差を算出し、当該周波数差の平均値を算出する第１ステップと、前記第１ステップにおいて算出された前記周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップとを含む。

本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と前記前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性を算出し、前記周波数差と前記プロファイル特性の積の平均値を算出する第１ステップと、前記第１ステップにおいて算出された前記周波数差と前記プロファイル特性の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップとを含む。

本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出する第１ステップと、前記第１ステップにおいて算出されたそれぞれの前記周波数差の中央値を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップとを含む。

本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の標準偏差を算出する第１ステップと、前記第１ステップにおいて算出された前記標準偏差に基づいて、前記複数の磁気共鳴信号のそれぞれについて算出された前記周波数差から基準となる複数の基準周波数差を選択し、選択された当該複数の基準周波数差の平均値を算出する第２ステップと、前記第２ステップにおいて算出された前記複数の基準周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップとを含む。

本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する磁気共鳴イメージング装置であって、前記被検体にＲＦパルスを照射し、前記複数の磁気共鳴信号を受信するスキャン部と、前記複数の磁気共鳴信号のプロファイルに基づいて、前記スキャン部により受信された前記複数の磁気共鳴信号から基準となる基準磁気共鳴信号を選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する変位量算出部と、前記変位量算出部により算出された前記変位量に基づいて、前記複数の磁気共鳴信号それぞれを変位させ、センターリング処理を実施するセンターリング実施部と、前記センターリング実施部により変位された前記複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像再構成部とを有する。

好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度の最大値を算出し、当該最大値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する。

好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて積分値を算出し、当該積分値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する。

好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて当該複数の磁気共鳴信号同士の相関係数を算出し、前記複数の磁気共鳴信号ごとに算出された前記相関係数の平均値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する。

好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数と前記ｋ空間の中心における周波数との差である周波数差を算出し、当該周波数差の平均値を算出し、当該周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する。

好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と、前記複数の磁気共鳴信号それぞれにおける積分値を算出し、前記周波数差と前記積分値の積の平均値を算出し、当該周波数差と前記積分値の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する。

好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と、前記複数の磁気共鳴信号それぞれにおける信号強度の最大値を算出し、前記周波数差と前記信号強度の最大値の積の平均値を算出し、当該周波数差と前記信号強度の最大値の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する。

好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の中央値を変位させる変位量として設定する。

好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の標準偏差を算出し、当該標準偏差に基づいて、前記複数の磁気共鳴信号のそれぞれについて算出された前記周波数差から基準となる複数の基準周波数差を選択し、選択された当該複数の基準周波数差の平均値を算出し、当該複数の基準周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する。

本発明によれば、パラレルイメージング法とハーフエコー法を併用した際に、生成した画像にアーチファクトが生じないような画像生成方法、エコー信号の位置補正方法、および磁気イメージング装置を提供することができる。

以下より、本発明にかかる一実施形態について図面を参照して説明する。

＜第１の実施形態＞
（装置構成）
図１は、本発明にかかる一実施形態におけるＲＦコイル部により構成される磁気共鳴イメージング装置の構成を示す構成図である。本装置は、本発明の実施形態の一例である。

図１に示すように、磁気共鳴イメージング装置１は、スキャン部２と、操作コンソール部３とを有している。ここで、スキャン部２は、静磁場マグネット部１２と、勾配コイル部１３と、ＲＦコイル部１４と、クレードル１５とを有する。そして、操作コンソール部３は、ＲＦ駆動部２２と、勾配駆動部２３と、データ収集部２４と、制御部３０と、記憶部３１と、操作部３２と、画像再構成部３３と、表示部３４とを有する。

スキャン部２について説明する。

スキャン部２は、図１に示すように、被検体４０において体動する撮影対象を含む撮影スライス領域が収容される静磁場空間１１を含んでいる。そして、スキャン部２は、操作コンソール部３からの制御信号に基づいて、その静磁場が形成される静磁場空間１１に収容した被検体４０の撮影領域にＲＦパルスを照射し、その撮影領域から生ずる磁気共鳴信号をイメージングデータとして収集するスキャンを被検体４０について実施する。

スキャン部２の各構成要素について、順次、説明する。図２は、フェーズドアレイコイルの一例を示す図である。理解促進の便宜上、表面コイルの重なり合いをずらして記載しており、またコイル線の太さを変更している。図３は、フェーズドアレイコイルにより受信したエコー信号にフーリエ変換処理を実施したデータを示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。周波数０の破線はｋ空間の中心Ｏを示す。

静磁場マグネット部１２は、被検体４０が収容される静磁場空間１１に静磁場を形成するために設けられている。静磁場マグネット部１２は、水平磁場型であって、被検体４０が収容される静磁場空間１１において載置される被検体４０の体軸方向（ｚ方向）に沿うように、超伝導磁石（図示なし）が静磁場を形成する。なお、静磁場マグネット部１２は、水平磁場型の他に、垂直磁場型であってもよく、また静磁場マグネット部１２は、永久磁石により構成されていてもよい。

勾配コイル部１３は、ＲＦコイル部１４が受信する磁気共鳴信号に３次元の位置情報を持たせるために、静磁場空間１１に勾配磁場を形成する。勾配コイル部１３は、スライス選択勾配磁場、読み取り勾配磁場、位相エンコード勾配磁場の３種類の勾配磁場を形成するために勾配コイルを３系統有する。

ＲＦコイル部１４は、たとえば、被検体４０を囲むように配置される。ＲＦコイル部１４は、静磁場マグネット部１２によって静磁場が形成される静磁場空間１１内において、制御部３０からの制御信号に基づいて、電磁波であるＲＦパルスを被検体４０に送信して高周波磁場を形成する。これにより、被検体４０の撮影スライス領域におけるプロトンのスピンを励起する。そして、ＲＦコイル部１４は、その励起された被検体４０の撮影スライス領域におけるプロトンのスピンが元の磁化ベクトルへ戻る際に生ずる電磁波を、磁気共鳴信号として受信する。

また、ＲＦコイル部１４は、例えば、図２に示す複数の独立した表面コイルからなるフェーズドアレイコイル１４１である。この複数の独立した表面コイルは、それぞれ個別に送信機および受信機と接続されており、多チャンネルで複数の磁気共鳴信号を受信するように構成されている。

図２に示すように、例えば、フェーズドアレイコイル１４１は、独立した表面コイルＣ１、Ｃ２およびＣ３からなる。そして、図３に示すように、表面コイルＣ１、Ｃ２およびＣ３により受信されたそれぞれのエコー信号についてフーリエ変換処理を実施したデータ（以下、フーリエ変換処理データとも称する）は、異なるｋ空間に表される。また、回転磁場不均一の度合いが場所によって異なるため、フーリエ変換処理データにおけるピークは、ｋ空間の中心Ｏからそれぞれ異なる周波数だけずれることとなる。例えば、図３（ａ）に示すように、ｋ空間に表された表面コイルＣ１のフーリエ変換処理データＤ１におけるピークは、ｋ空間の中心Ｏから周波数Ｆ_Ｄ1だけずれている。同様に図３（ｂ）、図３（ｃ）に示すように、フーリエ変換処理データＤ２、Ｄ３におけるピークは、ｋ空間の中心Ｏからそれぞれ周波数Ｆ_Ｄ２、Ｆ_Ｄ３だけずれており、それぞれｋ空間の中心Ｏから異なる周波数ずれている。

クレードル１５は、被検体４０を載置するテーブルを有する。クレードル１５は、制御部３０からの制御信号に基づいて、静磁場空間１１の内部と外部との間でテーブルに載置された被検体４０を移動する。

操作コンソール部３について説明する。

操作コンソール部３は、スキャン部２が被検体４０についてスキャンを実施するように制御し、そのスキャン部２が実施したスキャンによって得られた磁気共鳴信号に基づいて、被検体４０の画像を生成すると共に、その生成した画像を表示する。

操作コンソール部３を構成する各部について、順次、説明する。

ＲＦ駆動部２２は、ＲＦコイル部１４を駆動させて静磁場空間１１内に高周波磁場を形成するために、ゲート変調器（図示なし）とＲＦ電力増幅器（図示なし）とＲＦ発振器（図示なし）とを有する。ＲＦ駆動部２２は、制御部３０からの制御信号に基づいて、ＲＦ発振器からのＲＦ信号を、ゲート変調器を用いて所定のタイミングおよび所定の包絡線の信号に変調する。そして、ゲート変調器により変調されたＲＦ信号を、ＲＦ電力増幅器により増幅した後、ＲＦコイル部１４に出力する。

勾配駆動部２３は、制御部３０の制御信号に基づいて勾配コイル部１３を駆動させて、静磁場空間１１内に勾配磁場を発生させる。勾配駆動部２３は、勾配コイル部１３の３系統の勾配コイルに対応して３系統の駆動回路（図示なし）を有する。

データ収集部２４は、ＲＦコイル部１４が受信する磁気共鳴信号を収集するために、位相検波器（図示なし）とアナログ／デジタル変換器（図示なし）とを有する。データ収集部２４は、ＲＦコイル部１４からの磁気共鳴信号を、ＲＦ駆動部２２のＲＦ発振器の出力を参照信号として、位相検波器によって位相検波し、アナログ／デジタル変換器に出力する。そして、位相検波器により位相検波されたアナログ信号である磁気共鳴信号を、アナログ／デジタル変換器によってデジタル信号に変換して、画像再構成部３３に出力する。

制御部３０は、コンピュータと、コンピュータを用いて所定のスキャンに対応する動作を各部に実行させるプログラムとを有する。そして、制御部３０は、操作部３２に接続されており、操作部３２に入力された操作信号を処理し、クレードル１５とＲＦ駆動部２２と勾配駆動部２３とデータ収集部２４との各部に、制御信号を出力し制御を行う。また、制御部３０は、所望の画像を得るために、操作部３２からの操作信号に基づいて画像再構成部３３を制御する。

記憶部３１は、コンピュータにより構成されている。そして、記憶部３１は、データ収集部２４に収集された画像再構成処理前の磁気共鳴信号、画像再構成部３３で画像再構成処理された画像データ等を記憶する。

操作部３２は、キーボードやマウスなどの操作デバイスにより構成されている。操作部３２は、オペレータによって操作データ、撮像プロトコルなどが入力され、その操作データ、撮像プロトコルを制御部３０に出力する。

画像再構成部３３は、コンピュータにより構成されている。そして、画像再構成部３３は、データ収集部２４に接続されており、データ収集部２４から出力される磁気共鳴信号に対して画像再構成処理を実施して、画像を生成する。

図４は、本発明にかかる一実施形態における画像再構成部３３の構成を示すブロック図である。

図４に示すように、画像再構成部３３は、変位量算出部３０１と、センターリング実施部３０２とを有している。詳細については後述するが、制御部３０は、ＲＦコイル部１４により受信したエコー信号をｋ空間における周波数軸方向の中心Ｏにセンターリングさせるために、エコー信号を変位させる変位量を算出し、算出した変位量に基づいてセンターリング処理を施す。

画像再構成部３３の各部について順次説明する。

変位量算出部３０１は、複数の表面コイルにより受信された複数のエコー信号に基づき、エコー信号のピークをｋ空間における周波数軸方向の中心Ｏにセンターリングさせるためのエコー信号の変位量を算出する。

センターリング実施部３０２は、変位量算出部３０１が算出した変位量に基づき、受信した全てのエコー信号を変位させ、センターリング処理を実施する。

表示部３４は、ディスプレイなどの表示デバイスにより構成されており、画像再構成部３３が生成する被検体４０の画像を表示する。

（動作）
以下より、本実施形態の磁気共鳴イメージング装置１を用いて、被検体４０を撮像する際の動作について説明する。

図５は、本発明にかかる第１の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。

まず、図５に示すように、スキャンを実施する（ＳＴ１０）。

ここでは、スキャン部２が被検体４０のスキャンを実施する。

ここでは、オペレータによって操作部３２に入力された操作信号に基づいて、制御部３０は、ＲＦ駆動部２２と勾配駆動部２３とデータ収集部２４とにそれぞれ制御信号を出力し制御を行って、被検体４０をスキャンし、データ収集部２４にエコー信号を取得させる。そして、そのデータ収集部２４が取得したエコー信号を画像再構成部３３に出力させる。

つぎに、図５に示すように、画像の生成を実施する（ＳＴ２０）。

ここでは、データ収集部２４から出力されるエコー信号に対して画像再構成処理を実施して、画像を生成する。画像再構成処理の際、ｋ空間において、受信した複数のエコー信号のセンターリング処理を実施する。

図６は、本発明にかかる第１の実施形態における受信した複数の位相軸方向処理データのセンターリング処理の実施の動作を示すフロー図である。図７は、本発明にかかる第１の実施形態における周波数軸方向処理データを示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。独立した表面コイルのそれぞれにより受信されたエコー信号について高速フーリエ変換（以下、ＦＦＴとも称する）処理を実施したデータは異なるｋ空間に表されるが、理解補助の便宜上、同一のｋ空間に表している。また、理解補助の便宜上、周波数軸方向処理データＤ０１を太線、周波数軸方向処理データＤ０２を中太線、周波数軸方向処理データＤ０３を細線で記載している。図８は、本発明に係る第１の実施形態におけるノイズの影響のある周波数軸方向処理データを示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。図９および図１０は、本発明にかかる第１の実施形態において、位相軸方向処理データのセンターリング処理を示した図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。また、周波数０の破線は、ｋ空間における周波数軸方向の中心Ｏを示す。また、変位量Ｆ_Ｄ０１は、周波数軸方向処理データＤ０１の信号強度が最大となる点をｋ空間における周波数軸方向の中心Ｏに変位させるための変位量である。また、太線はホモダインハイパスフィルターＨＦ、ホモダインローパスフィルターＬＦを示し、理解促進の便宜上、ホモダインハイパスフィルターＨＦとホモダインローパスフィルターＬＦの重なり合いをずらして記載している（以下の図で同じ）。

図６に示すように、まず、iチャンネル目のエコー信号に対し、周波数方向にＦＦＴ処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０１１）。

ここでは、iチャンネル目の表面コイルＣｉにより受信されたエコー信号に対して周波数軸方向にＦＦＴ処理を実施する。

具体的には、画像再構成部３３が、iチャンネル目の表面コイルＣｉにより受信されたエコー信号に対して周波数軸方向にＦＦＴ処理を実施する。そして、エコー信号に対して周波数軸方向にＦＦＴ処理を実施することにより算出されたデータ（以下、周波数軸方向処理データとも称する）を記憶部３１が保存する。

次に、図６に示すように、信号強度の最大値を算出する（ＳＴ２０２１）。

ここでは、ステップＳＴ２０１１において算出された周波数軸方向処理データの信号強度の最大値を算出する。

具体的には、ステップＳＴ２０１１において、画像再構成部３３により算出された周波数軸方向処理データの信号強度の最大値を算出する。図７に示すように、iチャンネル目の表面コイルＣｉにより受信されたエコー信号に対して周波数軸方向にＦＦＴ処理を実施したデータが周波数軸方向処理データＤ０１の場合には、周波数軸方向処理データＤ０１における信号強度の最大値として、ピーク値Ｐ_Ｄ０１が算出される。また、データが周波数軸方向処理データＤ０２の場合、信号強度の最大値としてピーク値Ｐ_Ｄ０２が算出され、データが周波数軸方向処理データＤ０３の場合、信号強度の最大値としてピーク値Ｐ_Ｄ０３が算出される。図７の場合、それぞれの周波数軸方向処理データにおける信号強度の最大値はピーク値となっているが、例えば、図８に示すように、ノイズなどの影響のある周波数軸方向処理データの場合、信号強度の最大値がピーク値となるとは限らないが、この場合、信号強度の最大値として、ピーク値Ｐ_Ｄ０４ではなく、Ｎ_Ｄ０４が算出される。

次に、図６に示すように、位相軸方向にＦＦＴ処理を実施する（ＳＴ２０３１）。

ここでは、ステップＳＴ２０１１において周波数軸方向にＦＦＴ処理を実施すことにより算出されたデータに対して位相軸方向にＦＦＴ処理を実施する。

具体的には、ステップＳＴ２０１１において、画像再構成部３３により算出され、記憶部３１に保存されている周波数軸方向処理データに対して、位相軸方向にＦＦＴ処理を実施する。そして、位相軸方向にＦＦＴ処理を実施することにより算出されたデータ（以下、位相軸方向処理データとも称する）を記憶部３１が保存する。

次に、図６に示すように、ステップＳＴ２０３１において算出された位相軸方向処理データが、最後のチャンネルの表面コイルにより受信されたエコー信号から算出されたデータであるかの判断を行う（ＳＴ２０４１）。

ここでは、ステップＳＴ２０３１において算出された位相軸方向処理データが、最後のチャンネルの表面コイルにより受信されたエコー信号から算出されたデータであるかを、制御部３０が判断する。そして、位相軸方向処理データが最後のチャンネルの表面コイルにより受信されたエコー信号から算出されたデータでない場合（Ｎｏ）には、上記の周波数軸方向のＦＦＴ処理、信号強度の最大値の算出および位相軸方向のＦＦＴ処理を継続するように制御部３０が各部を制御する。

次に、図６に示すように、全てのチャンネルにおける信号強度の最大値を比較し、最大値の一番大きいチャンネルの変位量を算出する（ＳＴ２０５１）。

ここでは、ステップＳＴ２０２１において全てのチャンネルの表面コイルにより受信されたエコー信号から算出された周波数軸方向処理データにおける信号強度の最大値を比較する。そして、変位量算出部３０１が、ｋ空間において信号強度の最大値の一番大きいチャンネルの周波数軸方向処理データを基準周波数軸方向処理データとして、信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。

具体的には、図７に示すように、変位量算出部３０１が、周波数軸方向処理データＤ０１、Ｄ０２およびＤ０３のそれぞれの信号強度の最大値であるピーク値Ｐ_Ｄ０１、ピーク値Ｐ_Ｄ０２およびピーク値Ｐ_Ｄ０３を比較し、最大値の一番大きいデータである基準周波数軸方向処理データを周波数軸方向処理データＤ０１と決定する。ここで、図８に示す周波数軸方向処理データＤ０４のように、ピークにおける信号強度よりもノイズ部分における信号強度の方が高い周波数軸方向処理データが複数の周波数軸方向処理データの中に含まれており、ノイズ部分における信号強度Ｎ_Ｄ０４が全ての周波数軸方向処理データにおける信号強度の最大値の中で一番大きい場合には、周波数軸方向処理データＤ０４が選択されることとなる。しかし、ノイズ部分における信号強度がピークにおける信号強度よりも高くなる周波数軸方向処理データは、周波数軸方向処理データ全体として信号強度が低い。そのため、このような周波数軸方向処理データが選択されることはない。

そして、図９に示すように、変位量算出部３０１が、ｋ空間において周波数軸方向処理データＤ０１における信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、当該周波数差を変位量Ｆ_Ｄ０１とし、全ての位相軸方向処理データを変位させるための変位量として設定する。ここで、周波数軸方向処理データと位相軸方向処理データのピーク位置は変わらないため、位相軸方向処理データを変位させる変位量として、変位量Ｆ_Ｄ０１を用いることができる。

次に、図６に示すように、ｉチャンネル目のデータのセンターリング処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０６１）。

ここでは、ステップＳＴ２０５１において変位量算出部３０１が算出した変位量に基づいて、ステップＳＴ２０３１において画像再構成部３３が算出したiチャンネル目の表面コイルＣｉの位相軸方向処理データに対するセンターリング処理をセンターリング実施部３０２が実施する。

具体的には、例えば、iチャンネル目の表面コイルＣｉが表面コイルＣ１である場合には、ｋ空間において、図１０（ａ）に示す位相軸方向処理データＤ１１を、図１０（ｂ）に示すように変位量Ｆ_Ｄ０１だけ変位させるセンターリング処理をセンターリング実施部３０２が実施する。ステップＳＴ２０５１において算出された変位量Ｆ_Ｄ０１が、位相軸方向処理データＤ１１をｋ空間における周波数軸方向の中心Ｏに変位させるための変位量であるため、位相軸方向処理データＤ１１を変位量Ｆ_Ｄ０１だけ変位すると、位相軸方向処理データＤ１１の信号強度が最大となる点（ピーク）は、ｋ空間における周波数軸方向の中心Ｏに変位される。そして、ｋ空間における周波数軸方向の中心Ｏに変位された位相軸方向処理データＤ１１のピークがホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１１に適切にフィルターがかかることとなる。

また、図１０（ｃ）に示すように、ｋ空間において、位相軸方向処理データＤ１２におけるピークは、位相軸方向処理データＤ１１のピークと異なる周波数にあるため、位相軸方向処理データＤ１２を変位量Ｆ_Ｄ０１に基づいて変位させると、位相軸方向処理データＤ１２のピークはｋ空間における周波数軸方向の中心Ｏに変位されない。しかし、図１０（ｄ）に示すように、変位された位相軸方向処理データＤ１２のピークがｋ空間における周波数軸方向の中心ＯにあるホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１２に適切にフィルターがかかることとなる。したがって、アーチファクトは殆ど生じない。

また、位相軸方向処理データＤ１３についても同様に変位量Ｆ_Ｄ０１に基づいて、センターリング実施部３０２がセンターリング処理を実施する。位相軸方向処理データＤ１３も位相軸方向処理データＤ１２と同様に、変位された位相軸方向処理データＤ１３のピークがｋ空間における周波数軸方向の中心ＯにあるホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１３に適切にフィルターがかかることとなる。したがって、アーチファクトは殆ど生じない。

次に、図６に示すように、ハーフエコー処理を実施する（ＳＴ２０７１）。

ここでは、ＳＴ２０６１においてセンターリング処理を実施した位相軸方向処理データに対して、画像再構成部３３がハーフエコー処理を実施する。

次に、図６に示すように、ハーフエコー処理を実施したデータが最後のチャンネルのデータであるかの判断を行う（ＳＴ２０８１）。

ここでは、ＳＴ２０７１においてハーフエコー処理を実施することにより算出されたデータ（以下、ハーフエコー処理とも称する）が、最後のチャンネルにおける表面コイルにより受信されたエコー信号から算出されたかを、制御部３０が判断する。そして、ハーフエコー処理データが最後のチャンネルにおける表面コイルにより受信されたエコー信号から算出されたデータでない場合（Ｎｏ）には、上記のセンターリング処理およびハーフエコー処理を継続するように制御部３０が各部を制御する。

以上のように、本発明の第１の実施形態は、複数の表面コイルにより受信された全てのエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出されたそれぞれの周波数軸方向処理データにおける信号強度の最大値を算出する。そして、全ての周波数軸方向処理データにおける信号強度の最大値の中で、最大値の一番大きい基準周波数軸方向処理データを決定し、決定した基準周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、当該周波数差を全ての周波数軸方向処理データを変位させるための変位量と設定する。そして、算出された変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。

このように、最大値がピークとならない周波数軸方向処理データ、すなわちノイズ部分の方が実際のピークよりも信号強度が高くなる周波数軸方向処理データは、周波数軸方向処理データ全体として信号強度が低いため、信号強度の最大値が一番大きい基準周波数軸方向処理データとして選択されることはない。したがって、本発明にかかる第１の実施形態は、信号強度の最大値が一番大きい周波数軸方向処理データにおける信号強度が最大となる点をｋ空間における周波数軸方向の中心Ｏに変位させるための変位量に基づいて、全ての位相軸方向処理データをｋ空間において変位させることにより、生成した画像のアーチファクトを防止することができ、画像品質を向上させることができる。

＜第２の実施形態＞
以下より本発明にかかる第２の実施形態について図１１により説明する。
第２の実施形態は、動作フローにおけるステップＳＴ２０２２である積分値の算出、およびステップＳＴ２０５２である変位量の算出以外の部分は第１の実施形態と同じである。そのため、重複する箇所については、記載を省略する。

図１１は、本発明にかかる第２の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図１２は、本発明にかかる第２の実施形態における周波数軸方向処理データのプロファイルの積分値を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。また、周波数軸方向処理データのプロファイルと横軸で囲まれた斜線部が積分値であり、Ｓ_Ｄ０１、Ｓ_Ｄ０２およびＳ_Ｄ０３はそれぞれ周波数軸方向処理データＤ０１、Ｄ０２およびＤ０３の積分値を示す。

まず、図１１に示すように、積分値を算出する（ＳＴ２０２２）。

ここでは、ステップＳＴ２０１１において、変位量算出部３０１が周波数軸方向のＦＦＴ処理を実施することにより算出されたデータのｋ空間における積分値を算出する。

具体的には、ステップＳＴ２０１１において、変位量算出部３０１が画像再構成部３３により算出された周波数軸方向処理データのプロファイルの積分値を算出する。図１２（ａ）に示すように、iチャンネル目の表面コイルＣｉにより受信されたエコー信号に対して周波数軸方向にＦＦＴ処理を実施することにより算出されたデータが周波数軸方向処理データＤ０１の場合には、周波数軸方向処理データＤ０１のプロファイルの積分値Ｓ_Ｄ０１を算出する。積分値Ｓ_Ｄ０１は、周波数軸方向処理データＤ０１のデータプロファイルと横軸とで囲まれている斜線部分の面積である。また、算出されたデータが周波数軸方向処理データＤ０２、Ｄ０３の場合、積分値として、積分値Ｓ_Ｄ０２、積分値Ｓ_Ｄ０３のそれぞれが算出される。

次に、図１１に示すように、全てのチャンネルにおける積分値を比較し、積分値の一番大きいチャンネルにおける変位量を算出する（ＳＴ２０５２）。

ここでは、ＳＴ２０２２において算出した全ての周波数軸方向処理データのプロファイルの積分値を比較し、変位量算出部３０１が、ｋ空間において積分値の一番大きいチャンネルの周波数軸方向処理データを基準周波数軸方向処理データとして、信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。

具体的には、図１２に示すように、周波数軸方向処理データＤ０１、Ｄ０２、Ｄ０３のプロファイルのそれぞれの積分値Ｓ_Ｄ０１、積分値Ｓ_Ｄ０２、積分値Ｓ_Ｄ０３を比較し、積分値の一番大きいデータである基準周波数軸方向処理データを周波数軸方向処理データＤ０１と決定する。そして、図９に示すように、変位量算出部３０１が、ｋ空間において周波数軸方向処理データＤ０１の信号強度が最大となる点における周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、当該周波数差を変位量Ｆ_Ｄ０１とし、全ての位相軸方向処理データを変位させるための変位量として設定する。ここで、周波数軸方向処理データと位相軸方向処理データのピーク位置は変わらないため、位相軸方向処理データを変位させる変位量として、変位量Ｆ_Ｄ０１を用いることができる。

以上のように、本発明の第２の実施形態は、複数の表面コイルにより受信した全てのエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出されたそれぞれの周波数軸方向処理データのプロファイルの積分値を算出する。そして、全ての周波数軸方向処理データのプロファイルの積分値から積分値の一番大きい基準周波数軸方向処理データを決定し、決定した周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、当該周波数差を全ての周波数軸方向処理データを変位させるための変位量として設定する。そして、算出された変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。

このように、周波数軸方向処理データの最大値がピークとならない場合、すなわち、ノイズ部分の方が実際のピークよりも信号強度が高くなる周波数軸方向処理データは、周波数軸方向処理データのプロファイルの積分値が小さいため、本実施形態における基準周波数軸方向処理データとして選択されることはない。したがって、本発明にかかる第２の実施形態は、周波数軸方向処理データのプロファイルの積分値が一番大きい周波数軸方向処理データの信号強度が最大となる点をｋ空間における周波数軸方向の中心Ｏに変位させるための変位量に基づいて、ｋ空間において全ての位相軸方向処理データを変位させることにより、生成した画像のアーチファクトを防止することができ、画像品質を向上させることができる。

＜第３の実施形態＞
以下より本発明にかかる第３の実施形態について図１３により説明する。
第３の実施形態は、動作フローにおけるステップＳＴ２０５３である変位量の算出および、ステップＳＴ２０６３であるセンターリング処理の実施以外の部分は第１の実施形態と同じである。そのため、重複する箇所については、記載を省略する。

図１３は、本発明にかかる第３の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図１４は、本発明にかかる第３の実施形態における周波数軸方向処理データそれぞれの相関を示す相関図である。縦軸は信号強度であり、横軸は周波数である。理解補助の便宜上、周波数軸方向処理データＤ０１を太線、周波数軸方向処理データＤ０２を中太線、周波数軸方向処理データＤ０３を細線で記載している。図１５は、本発明にかかる第３の実施形態における周波数軸方向処理データの相関係数が最大である周波数軸方向処理データの周波数差に基づいて、位相軸方向処理データのセンターリング処理を示した図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。

まず、図１３に示すように、各チャンネルにおけるデータの相関をとり、相関係数の一番大きいチャンネルの変位量を算出する（ＳＴ２０５３）。

ここでは、変位量算出部３０１が、各チャンネルにおける周波数軸方向処理データごとに他のチャンネルにおける周波数軸方向処理データとの相関をとり、相関係数を算出する。そして、変位量算出部３０１が、ｋ空間において相関係数が一番大きいチャンネルの周波数軸方向処理データを基準周波数軸方向処理データとして、信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。

具体的には、例えば、図１４（ａ）に示すように、周波数軸方向処理データＤ０１と周波数軸方向処理データＤ０２との相関をとり、周波数軸方向処理データＤ０１と周波数軸方向処理データＤ０２との相関係数Ｒ_１２を変位量算出部３０１が算出する。そして、図１４（ｂ）および図１４（ｃ）に示すように、周波数軸方向処理データＤ０１と周波数軸方向処理データＤ０３との相関係数Ｒ_１３および周波数軸方向処理データＤ０２と周波数軸方向処理データＤ０３との相関係数Ｒ_２３も変位量算出部３０１が算出する。相関係数は、数式（１）により算出する。

次に、変位量算出部３０１が、それぞれの周波数軸方向処理データにおける相関係数の平均値を算出する。例えば、周波数軸方向処理データＤ０１における相関係数の平均値は（Ｒ_１２＋Ｒ_１３）／２である。同様に周波数軸方向処理データＤ０２における相関係数の平均値は（Ｒ_１２＋Ｒ_2３）／２であり、周波数軸方向処理データＤ０３における相関係数の平均値は（Ｒ_１３＋Ｒ_2３）／２である。そして、変位量算出部３０１が、それぞれの平均値を比較し、平均値の一番大きいデータである基準周波数軸方向処理データを決定する。そして、平均値の一番大きい周波数軸方向処理データにおける、信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を、変位量算出部３０１が算出する。そして、当該周波数差を全ての位相軸方向処理データを変位させるための変位量として設定する。本実施形態の場合、相関係数の一番大きいデータである基準周波数軸方向処理データは、周波数軸方向処理データＤ０３である。したがって、設定される変位量は、変位量Ｆ_Ｄ０３となる。ここで、周波数軸方向処理データと位相軸方向処理データのピーク位置は変わらないため、位相軸方向処理データを変位させる変位量として、変位量Ｆ_Ｄ０３を用いることができる。

次に、図１３に示すように、iチャンネル目のデータのセンターリング処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０６３）。

ここでは、ＳＴ２０５３において変位量算出部３０１が算出した変位量に基づいて、センターリング実施部３０２がｉチャンネル目の表面コイルＣｉの位相軸方向処理データに対してセンターリング処理を実施する。

具体的には、例えば、図１５（ａ）に示すように、ｉチャンネル目の表面コイルＣｉが表面コイルＣ１である場合には、ｋ空間において、図１５（ａ）に示す位相軸方向処理データＤ１１を、図１５（ｂ）に示すように、変位量Ｆ_Ｄ０３だけ変位させて、センターリング実施部３０２がセンターリング処理を実施する。この場合、位相軸方向処理データＤ１１におけるピークが、ｋ空間における周波数軸方向の中心Ｏから（Ｆ_Ｄ０３−Ｆ_Ｄ０１）だけずれるが、ホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１１に適切にフィルターがかかることとなる。

また、上記と同様に位相軸方向処理データＤ１２および位相軸方向処理データＤ１３に対しても変位量Ｆ_Ｄ０３に基づいてｋ空間において変位させて、センターリング実施部３０２がセンターリング処理を実施する。この場合にも、ピークが、ｋ空間における周波数軸方向の中心ＯにあるホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１２およびＤ１３に適切にかかることとなる。

以上のように、本発明の第３の実施形態は、複数の表面コイルにより受信したエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出されたそれぞれの周波数軸方向処理データについて他の周波数軸方向処理データとの相関をとり、相関係数を算出する。そして、相関係数の一番大きいチャンネルの周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、当該周波数差を全ての周波数軸方向処理データを変位させるための変位量と設定する。そして、算出された変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。

このように、周波数軸方向処理データの最大値がピークとならない場合、すなわち、ノイズ部分の方が実際のピークよりも信号強度が高くなる周波数軸方向処理データは、他の周波数軸方向処理データとの相関が低いため、本実施形態における基準周波数軸方向処理データとして選択されることはない。したがって、本発明にかかる第３の実施形態は、相関係数の一番大きい周波数軸方向処理データにおける信号強度が最大となる点をｋ空間における周波数軸方向の中心Ｏに変位させるための変位量に基づいて、全ての位相軸方向処理データをｋ空間において変位させることにより、生成した画像のアーチファクトを防止することができる。よって、画像品質を向上させることができる。

＜第４の実施形態＞
以下より本発明にかかる第４の実施形態について図１６により説明する。
第４の実施形態は、動作フローにおけるステップＳＴ２０５４である変位量の平均値の算出および、ステップＳＴ２０６４であるセンターリング処理の実施以外の部分は第１の実施形態と同じである。そのため、重複する箇所については、記載を省略する。

図１６は、本発明にかかる第４の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図１７は、本発明にかかる第４の実施形態における周波数軸方向処理データの信号強度の最大値とｋ空間における周波数軸方向の中心Ｏの周波数差を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。図１８は、本発明にかかる第４の実施形態における全ての周波数差の平均値である変位量に基づいて、位相軸方向処理データのセンターリング処理を示した図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。また、変位量Ｆ_ＡＶ１は、全ての周波数軸方向処理データの変位量の平均値である。

まず、図１６に示すように、全てのチャンネルにおけるデータの周波数差の平均値を算出する（ＳＴ２０５４）。

ここでは、ＳＴ２０２１において算出された、それぞれの周波数軸方向処理データの信号強度が最大値を示す点における周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、全ての周波数軸方向処理データにおける周波数差の平均値を算出する。

具体的には、例えば、図１７（ａ）に示すように、変位量算出部３０１が、周波数軸方向処理データＤ０１の信号強度が最大となる点における周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差Ｆ_Ｄ０１を算出する。また、同様に図１７（ｂ）、図１７（ｃ）に示すように、変位量算出部３０１が、周波数軸方向処理データＤ０２、Ｄ０３それぞれの信号強度が最大となる点における周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差Ｆ_Ｄ０２、周波数差Ｆ_Ｄ０３をそれぞれ算出する。そして、それぞれの周波数軸方向処理データの周波数差の平均値である（Ｆ_Ｄ０１＋Ｆ_Ｄ０２＋Ｆ_Ｄ０３）／３を変位量算出部３０１が算出する。この平均値を全てのチャンネルの位相軸方向処理データの変位量Ｆ_ＡＶ１として設定する。ここで、周波数軸方向処理データと位相軸方向処理データのピーク位置は変わらないため、位相軸方向処理データを変位させる変位量として、変位量Ｆ_ＡＶ１を用いることができる。

次に、図１６に示すように、iチャンネル目のデータのセンターリング処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０６４）。

ここでは、ＳＴ２０５４において変位量算出部３０１が算出した変位量Ｆ_ＡＶ１に基づいて、ステップＳＴ２０３１において画像再構成部３３が算出したiチャンネル目の表面コイルＣｉの位相軸方向処理データに対するセンターリング処理をセンターリング実施部３０２が実施する。

具体的には、例えば、iチャンネル目の表面コイルＣｉが表面コイルＣ１である場合には、ｋ空間において、図１８（ａ）に示す位相軸方向処理データＤ１１を、図１８（ｂ）に示すように変位量Ｆ_ＡＶ１だけ変位させて、センターリング実施部３０２がセンターリング処理を実施する。センターリング処理を実施した後の位相軸方向処理データＤ１１における信号強度が最大となる点は、ｋ空間における周波数軸方向の中心Ｏから（Ｆ_ＡＶ１−Ｆ_Ｄ０１）だけずれているが、位相軸方向処理データＤ１１のピークがホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１１に適切にフィルターがかかることとなる。また、上記と同様に、iチャンネル目の表面コイルＣｉが表面コイルＣ２またはＣ３である場合にも、位相軸方向処理データＤ１２またはＤ１３を変位量Ｆ_ＡＶ１だけ変位させて、センターリング処理を実施させる。そして、この場合にも、位相軸方向処理データＤ１２、Ｄ１３における信号強度が最大となる点はそれぞれ、ｋ空間における周波数軸方向の中心Ｏから（Ｆ_ＡＶ１−Ｆ_Ｄ０２）、（Ｆ_ＡＶ１−Ｆ_Ｄ０３）だけずれている。しかし、位相軸方向処理データＤ１２、Ｄ０３のピークがホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１２、Ｄ０３に適切にフィルターがかかることとなる。

上記のように、位相軸方向処理データＤ１１、Ｄ１２およびＤ１３をｋ空間において、変位量Ｆ_ＡＶ１だけ変位させると、それぞれの周波数軸方向処理データにおける信号強度が最大となる点は、ｋ空間における周波数軸方向の中心Ｏには変位しない。しかし、ｋ空間における周波数軸方向の中心ＯにあるホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦが、変位された位相軸方向処理データＤ１１、Ｄ１２およびＤ１３に適切にかかるため、アーチファクトは殆ど生じない。

以上のように、本発明の第４の実施形態は、複数の表面コイルにより受信した全てのエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出された、それぞれの周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。そして、全ての周波数軸方向処理データにおける周波数差の平均を算出する。そして、算出した周波数差の平均値を変位量として設定し、当該変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。

このように、本発明にかかる第４の実施形態は、全ての周波数軸方向処理データそれぞれについて、信号強度が最大となる点をｋ空間における周波数軸方向の中心Ｏに変位させるための変位量を算出し、その全ての変位量の平均値に基づいて、全ての位相軸方向処理データをｋ空間において変位させることにより、ノイズの影響のある周波数軸方向処理データが含まれていたとしても、殆ど影響が及ぼされることはない。したがって、生成した画像のアーチファクトを防止することができ、画像品質を向上させることができる。

＜第５の実施形態＞
以下より本発明にかかる第５の実施形態について図１９により説明する。
第５の実施形態は、動作フローにおけるステップＳＴ２０５５である積分値と変位量の積の平均値の算出および、ステップＳＴ２０６５であるセンターリングの処理の実施以外の部分は第２の実施形態と同じである。そのため、重複する箇所については、記載を省略する。

図１９は、本発明にかかる第５の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図２０は、本発明にかかる第５の実施形態における周波数軸方向処理データにおけるプロファイルの積分値と周波数差を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。図２１は、本発明にかかる第５の実施形態における周波数軸方向処理データにおけるプロファイルの積分値と周波数差との積の平均値に基づいて、位相軸方向処理データのセンターリング処理を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。変位量Ｆ_ＡＶ２は、全ての周波数軸方向処理データのプロファイルの積分値と変位量の積の平均値である。

まず、図１９に示すように、全てのチャンネルにおける積分値と周波数差の積の平均値を算出する（ＳＴ２０５５）。

ここでは、変位量算出部３０１が各チャンネルにおける周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。そして、当該周波数差とステップＳＴ２０２２で算出した各周波数軸方向処理データプロファイルの積分値との積を変位量算出部３０１が算出し、それらの平均値を算出する。

具体的には、例えば、図２０（ａ）に示すように、変位量算出部３０１が、周波数軸方向処理データＤ０１における信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差Ｆ_Ｄ０１を算出する。そして、周波数軸方向処理データＤ０１におけるプロファイルの積分値と周波数差の積である（Ｓ_Ｄ０１×Ｆ_Ｄ０１）を変位量算出部３０１が算出する。
同様に、周波数軸方向処理データＤ０２におけるプロファイルの積分値と周波数差の積（Ｓ_Ｄ０２×Ｆ_Ｄ０２）および周波数軸方向処理データＤ０３におけるプロファイルの積分値と周波数差の積（Ｓ_Ｄ０３×Ｆ_Ｄ０３）を変位量算出部３０１が算出する。

次に、全ての周波数軸方向処理データにおけるプロファイルの積分値と周波数差の積の平均値を算出する。そして、当該積分値と周波数差の積の平均値を全ての位相軸方向処理データを変位させるための変位量として設定する。
例えば、周波数軸方向処理データＤ０１、Ｄ０２およびＤ０３におけるプロファイルの積分値と周波数差の積の平均値を算出する。そして、当該平均値を全てのチャンネルにおける位相軸方向処理データの変位量Ｆ_ＡＶ２として設定する。ここで、周波数軸方向処理データと位相軸方向処理データのピーク位置は変わらないため、位相軸方向処理データを変位させる変位量として、変位量Ｆ_ＡＶ２を用いることができる。また、最大値と周波数差の積の平均値である変位量Ｆ_ＡＶ２は以下の数式（２）により算出される。

次に、図１９に示すように、iチャンネル目のデータのセンターリング処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０６５）。

ここでは、ＳＴ２０５５において変位量算出部３０１が算出した変位量Ｆ_ＡＶ２に基づいて、センターリング実施部３０２がiチャンネル目の表面コイルＣｉの位相軸方向処理データに対してセンターリング処理を実施する。

具体的には、例えば、iチャンネル目の表面コイルＣｉが表面コイルＣ１である場合には、ｋ空間において、図２１（ａ）に示す位相軸方向処理データＤ１１を、図２１（ｂ）に示すように、変位量Ｆ_ＡＶ２だけ変位させて、センターリング実施部３０２がセンターリング処理を実施する。この場合、位相軸方向処理データＤ１１におけるピークは、ｋ空間における周波数軸方向の中心Ｏから（Ｆ_ＡＶ２−Ｆ_Ｄ０１）だけずれているが、ホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１１に適切にフィルターがかかることとなる。また、iチャンネル目の表面コイルＣｉが表面コイルＣ２またはＣ３である場合にも、位相軸方向処理データＤ１１と同様に、位相軸方向処理データＤ１２またはＤ１３をそれぞれ変位量Ｆ_ＡＶ２だけ変位させて、センターリング処理を実施する。そして、この場合にも、位相軸方向処理データＤ１２またはＤ１３におけるそれぞれのピークは、ｋ空間における周波数軸方向の中心Ｏから（Ｆ_ＡＶ２−Ｆ_Ｄ０２）、（Ｆ_ＡＶ２−Ｆ_Ｄ０３）だけずれているが、ホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１２またはＤ１３に適切にフィルターがかかることとなる。

以上のように、本発明の第５の実施形態は、複数の表面コイルにより受信した全てのエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出されたそれぞれの周波数軸方向処理データにおけるプロファイルの積分値を算出する。そして、それぞれの周波数軸方向処理データにおける信号強度が最大となる周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。そして、それぞれの周波数軸方向処理データにおけるプロファイルの積分値と周波数差の積を算出し、全てのエコー信号におけるプロファイルの積分値と周波数差の積の平均値を算出する。そして、算出した平均値を変位量と設定し、当該変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。

このように、本発明にかかる第５の実施形態は、全ての周波数軸方向処理データにおけるプロファイルの積分値と周波数差との積の平均値に基づいて、全ての位相軸方向処理データをｋ空間において変位させることにより、ノイズの影響のある位相軸方向処理データが含まれていたとしても、殆ど影響が及ぼされることはない。したがって、生成した画像のアーチファクトを防止することができ、画像品質を向上させることができる。

＜第６の実施形態＞
以下より本発明にかかる第６の実施形態について図２２により説明する。
第６の実施形態は、動作フローにおけるステップＳＴ２０５６である最大値と変位量の積の平均値の算出および、ステップＳＴ２０６６であるセンターリングの処理の実施以外の部分は第１の実施形態と同じである。そのため、重複する箇所については、記載を省略する。

図２２は、本発明にかかる第６の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図２３は、本発明にかかる第６の実施形態における周波数軸方向処理データにおける信号強度の最大値と周波数差を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。
図２４は、本発明にかかる第６の実施形態における周波数軸方向処理データにおける信号強度の最大値と周波数差との積の平均値に基づいて、位相軸方向処理データのセンターリング処理を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。また、変位量Ｆ_ＡＶ３は、全ての周波数軸方向処理データの最大値と変位量の積の平均値である。

まず、図２２に示すように、全てのチャンネルにおける信号強度の最大値と周波数差の積の平均値を算出する（ＳＴ２０５６）。

ここでは、変位量算出部３０１が各チャンネルにおける周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。そして、変位量算出部３０１が当該周波数差とステップＳＴ２０２１で算出した各周波数軸方向処理データの信号強度の最大値との積を算出し、周波数差と信号強度の最大値との積の平均値を算出する。

具体的には、例えば、図２３（ａ）に示すように、変位量算出部３０１が、周波数軸方向処理データＤ０１における信号強度が最大（Ｐ_Ｄ０１）となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差Ｆ_Ｄ０１を算出する。そして、周波数軸方向処理データＤ０１における信号強度の最大値と周波数差の積である（Ｐ_Ｄ０１×Ｆ_Ｄ０１）を変位量算出部３０１が算出する。
同様に、周波数軸方向処理データＤ０２における信号強度の最大値と周波数差の積（Ｐ_Ｄ０２×Ｆ_Ｄ０２）および周波数軸方向処理データＤ０３における最大値と周波数差の積（Ｐ_Ｄ０３×Ｆ_Ｄ０３）を変位量算出部３０１が算出する。

次に、全ての周波数軸方向処理データにおける最大値と周波数差の積の平均値を算出する。そして、当該最大値と周波数差の積の平均値を全ての位相軸方向処理データを変位させるための変位量として設定する。
例えば、周波数軸方向処理データＤ０１、Ｄ０２およびＤ０３における最大値と周波数差の積の平均値を算出する。そして、当該平均値を全てのチャンネルにおける位相軸方向処理データの変位量Ｆ_ＡＶ３として設定する。ここで、周波数軸方向処理データと位相軸方向処理データのピーク位置は変わらないため、位相軸方向処理データを変位させる変位量として、変位量Ｆ_ＡＶ３を用いることができる。また、最大値と周波数差の積の平均値である変位量Ｆ_ＡＶ３は以下の数式（３）により算出される。

次に、図２２に示すように、iチャンネル目のデータのセンターリング処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０６６）。

ここでは、ＳＴ２０５６において変位量算出部３０１が算出した変位量Ｆ_ＡＶ３に基づいて、センターリング実施部３０２がiチャンネル目の表面コイルＣｉの位相軸方向処理データに対してセンターリング処理を実施する。

具体的には、例えば、iチャンネル目の表面コイルＣｉが表面コイルＣ１である場合には、ｋ空間において、図２４（ａ）に示す位相軸方向処理データＤ１１を、図２４（ｂ）に示すように、変位量Ｆ_ＡＶ３だけ変位させて、センターリング実施部３０２がセンターリング処理を実施する。この場合、位相軸方向処理データＤ１１におけるピークは、ｋ空間における周波数軸方向の中心Ｏから（Ｆ_ＡＶ３−Ｆ_Ｄ０１）だけずれているが、ホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１１に適切にフィルターがかかることとなる。また、iチャンネル目の表面コイルＣｉが表面コイルＣ２またはＣ３である場合にも、位相軸方向処理データＤ１１と同様に、位相軸方向処理データＤ１２またはＤ１３をそれぞれ変位量Ｆ_ＡＶ２だけ変位させて、センターリング処理を実施する。そして、この場合にも、位相軸方向処理データＤ１２またはＤ１３におけるそれぞれのピークは、ｋ空間における周波数軸方向の中心Ｏから（Ｆ_ＡＶ３−Ｆ_Ｄ０２）、（Ｆ_ＡＶ３−Ｆ_Ｄ０３）だけずれているが、ホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１２またはＤ１３に適切にフィルターがかかることとなる。

以上のように、本発明の第６の実施形態は、複数の表面コイルにより受信した全てのエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出されたそれぞれの周波数軸方向処理データにおける信号強度の最大値を算出する。そして、それぞれの周波数軸方向処理データにおける信号強度が最大となる周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。そして、それぞれの周波数軸方向処理データにおける信号強度の最大値と周波数差の積を算出し、全てのエコー信号における信号強度の最大値と周波数差の積の平均値を算出する。そして、算出した平均値を変位量と設定し、当該変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。

このように、本発明にかかる第６の実施形態は、全ての周波数軸方向処理データにおける信号強度の最大値と周波数差との積の平均値に基づいて、全ての位相軸方向処理データをｋ空間において変位させることにより、ノイズの影響のある位相軸方向処理データが含まれていたとしても、殆ど影響が及ぼされることはない。したがって、生成した画像のアーチファクトを防止することができ、画像品質を向上させることができる。

＜第７の実施形態＞
以下より本発明にかかる第７の実施形態について図２５により説明する。
第７の実施形態は、動作フローにおけるステップＳＴ２０５７である全てのチャンネルの変位量の中央値を算出の算出および、ステップＳＴ２０６７であるセンターリングの処理の実施以外の部分は第１の実施形態と同じである。そのため、重複する箇所については、記載を省略する。

図２５は、本発明にかかる第７の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図２６は、本発明にかかる第７の実施形態におけるそれぞれの周波数軸方向処理データを周波数軸方向の中心Ｏに変位させるための変位量の中央値に基づいて、位相軸方向処理データのセンターリング処理を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。

まず、図２５に示すように、全てのチャンネルにおけるデータのそれぞれの変位量の中央値を算出（ＳＴ２０５７）。

ここでは、変位量算出部３０１が各チャンネルにおける周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、全てのチャンネルにおける周波数軸方向処理データにおける周波数差を昇順または降順に並べたときに中央となる周波数差である中央値を算出する。

具体的には、例えば、図２０に示すように、変位量算出部３０１が、周波数軸方向処理データＤ０１における信号強度が最大（Ｐ_Ｄ０１）となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差Ｆ_Ｄ０１を算出する。
同様に、周波数軸方向処理データＤ０２における周波数差Ｆ_Ｄ０２および周波数軸方向処理データＤ０３における周波数差Ｆ_Ｄ０３を変位量算出部３０１が算出する。

次に、変位量算出部３０１が全ての周波数軸方向処理データにおけるそれぞれの周波数差の中央値を算出し、中央値を全ての位相軸方向処理データを変位させるための変位量として設定する。
例えば、図２０に示すように、周波数差の大きさがＦ_Ｄ０１＜Ｆ_Ｄ０２＜Ｆ_Ｄ０３となる場合、周波数差の中央値は周波数差Ｆ_Ｄ０２となる。そして、この中央値である周波数差Ｆ_Ｄ０２を全ての位相軸方向処理データを変位させるための変位量として設定する。

次に、図２５に示すように、iチャンネル目のデータのセンターリング処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０６７）。

ここでは、ＳＴ２０５７において変位量算出部３０１が算出した変位量に基づいて、センターリング実施部３０２がiチャンネル目の表面コイルＣｉの位相軸方向処理データに対してセンターリング処理を実施する。

具体的には、例えば、図２６（ａ）に示すように、iチャンネル目の表面コイルＣｉが表面コイルＣ１である場合には、ｋ空間において、図２６（ａ）に示す位相軸方向処理データＤ１１を、図２６（ｂ）に示すように、変位量Ｆ_Ｄ０２だけ変位させて、センターリング実施部３０２がセンターリング処理を実施する。この場合、位相軸方向処理データＤ１１におけるピークが、ｋ空間における周波数軸方向の中心Ｏから（Ｆ_Ｄ０２−Ｆ_Ｄ０１）だけずれるが、ホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１１に適切にフィルターがかかることとなる。

また、上記と同様に位相軸方向処理データＤ１２および位相軸方向処理データＤ１３に対しても変位量Ｆ_Ｄ０２に基づいてｋ空間において変位させて、センターリング処理を実施する。この場合にも、ピークが、ｋ空間における周波数軸方向の中心ＯにあるホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１２およびＤ１３に適切にかかることとなる。

以上のように、本発明の第７の実施形態は、複数の表面コイルにより受信した全てのエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出されたそれぞれの周波数軸方向処理データにおける信号強度の最大値を算出する。そして、それぞれの周波数軸方向処理データにおける信号強度が最大となる周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、全てのチャンネルにおける周波数軸方向処理データにおける周波数差の中央値を算出する。そして、当該中央値を全ての周波数軸方向処理データを変位させるための変位量と設定する。そして、算出された変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。

このように、全てのチャンネルにおける周波数軸方向処理データにおける周波数差の中央値を、周波数軸方向処理データを変位させるための変位量とすることにより、全ての周波数軸方向処理データのピーク値が、ｋ空間における周波数軸方向の中心Ｏの近傍に変位されるため、全てのチャンネルにおける周波数軸方向処理データにおける周波数差の中央値に基づいて、全ての位相軸方向処理データをｋ空間において変位させることにより、生成した画像のアーチファクトを防止することができる。よって、画像品質を向上させることができる。

＜第８の実施形態＞
以下より本発明にかかる第８の実施形態について図２７により説明する。
第８の実施形態は、動作フローにおけるステップＳＴ２０５８である最大値と変位量の積の平均値の算出および、ステップＳＴ２０６８であるセンターリングの処理の実施以外の部分は第１の実施形態と同じである。そのため、重複する箇所については、記載を省略する。

図２７は、本発明にかかる第８の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図２８は、本発明にかかる第８の実施形態における周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差の標準偏差に基づいて、選択された周波数差の平均値により位相軸方向処理データのセンターリング処理を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。また、変位量Ｆ_ＡＶ４は、標準偏差に基づいて選択された周波数差の平均値である。

まず、図２７に示すように、標準偏差に基づいて選択されたデータの変位量の平均値を算出（ＳＴ２０５８）。

ここでは、変位量算出部３０１が各チャンネルにおける周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、全てのチャンネルにおける周波数軸方向処理データにおける周波数差の標準偏差を算出する。そして、標準偏差に基づいて周波数差を選択し、選択された周波数差の平均値を算出する。

次に、変位量算出部３０１が全ての周波数軸方向処理データにおける周波数差の標準偏差を算出し、当該標準偏差に基づいて、全ての周波数差の中から周波数差を選択する。
例えば、周波数軸方向処理データＤ０１、Ｄ０２およびＤ０３における周波数差の相加平均値と標準偏差σを算出する。そして、｜周波数差−相加平均値｜＜２σとなる周波数差を選択する。そして、選択された周波数差の平均値を全てのチャンネルにおける位相軸方向処理データの変位量Ｆ_ＡＶ４として設定する。例えば、周波数軸方向処理データＤ０１、Ｄ０２およびＤ０３において、｜周波数差−相加平均値｜＜２σとなる周波数軸方向処理データがＤ０１、Ｄ０２の場合、変位量Ｆ_ＡＶ４＝（Ｆ_Ｄ０１＋Ｆ_Ｄ０２）／２となる。
ここで、周波数軸方向処理データと位相軸方向処理データのピーク位置は変わらないため、位相軸方向処理データを変位させる変位量として、変位量Ｆ_ＡＶ３を用いることができる。また、周波数差の標準偏差は以下の数式（４）により算出される。

次に、図２７に示すように、iチャンネル目のデータのセンターリング処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０６８）。

ここでは、ＳＴ２０５７において変位量算出部３０１が算出した変位量に基づいて、iチャンネル目の表面コイルＣｉの位相軸方向処理データに対してセンターリング処理を実施する。

具体的には、例えば、図２８（ａ）に示すように、iチャンネル目の表面コイルＣｉが表面コイルＣ１である場合には、ｋ空間において、図２８（ａ）に示す位相軸方向処理データＤ１１を、図２８（ｂ）に示すように、変位量Ｆ_ＡＶ３だけ変位させて、センターリング処理を実施する。この場合、位相軸方向処理データＤ１１におけるピークが、ｋ空間における周波数軸方向の中心Ｏから（Ｆ_ＡＶ４−Ｆ_Ｄ０１）だけずれるが、ホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１１に適切にフィルターがかかることとなる。

また、上記と同様に位相軸方向処理データＤ１２および位相軸方向処理データＤ１３に対しても変位量Ｆ_ＡＶ４に基づいてｋ空間において変位させて、センターリング処理を実施する。この場合にも、ピークが、ｋ空間における周波数軸方向の中心ＯにあるホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１２およびＤ１３に適切にかかることとなる。

以上のように、本発明の第８の実施形態は、複数の表面コイルにより受信した全てのエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出されたそれぞれの周波数軸方向処理データにおける信号強度の最大値を算出する。そして、それぞれの周波数軸方向処理データにおける信号強度が最大となる周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、全てのチャンネルにおける周波数軸方向処理データにおける周波数差の標準偏差を算出する。そして、当該標準偏差に基づいて、全てのチャンネルにおける周波数軸方向処理データにおける周波数差から周波数差を選択し、選択した周波数差の平均値を全ての周波数軸方向処理データを変位させるための変位量と設定する。そして、算出された変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。

このように、全てのチャンネルにおける周波数軸方向処理データにおける周波数差の標準偏差に基づいて、周波数差を選択することにより、極端に値の異なる周波数差を、平均値を算出するための周波数差から除外することができる。そのため、標準偏差に基づいて選択した周波数差における平均値に基づいて、全ての位相軸方向処理データをｋ空間において変位させることにより、生成した画像のアーチファクトを防止することができる。よって、画像品質を向上させることができる。

なお、上記の本実施形態におけるエコー信号は、本発明の磁気共鳴信号に相当する。また、本実施形態の基準周波数軸方向処理データは、本発明の基準磁気共鳴信号に相当する。また、本実施形態のスキャン部２は、本発明のスキャン部に相当する。また、本実施形態の変位量算出部３０１は、本発明の変位量算出部に相当する。また、本実施形態のセンターリング実施部３０２は、本発明のセンターリング実施部に相当する。また、本実施形態の画像再構成部３３は、本発明の画像再構成部に相当する。

なお、本発明の実施に際しては、上記した実施形態に限定されるものではなく、種々の変形形態を採用することができる。

本発明の実施形態において、エコー信号にＦＦＴ処理を実施したデータに対して、センターリング処理を実施しているが、これに限定されず、エコー信号に対してセンターリング処理を実施したデータに対して、ＦＦＴ処理を実施しても良い。また、本発明の第５の実施形態において、積分値と周波数差との積の平均値に基づいて、センターリング処理を実施しているが、これに限定されず、エコー信号のプロファイル特性と周波数差との積の平均値に基づいて、センターリング処理を実施しても良い。例えば、エコー信号のプロファイル特性をエコー信号プロファイル全体の平均値としても良い。また、その他のエコー信号のプロファイル特性であっても良い。

図１は、本発明にかかる一実施形態のＲＦコイル部により構成される磁気共鳴イメージング装置の構成を示す構成図である。図２に、フェーズドアレイコイルの一例を示す。図３は、フェーズドアレイコイルにより受信したエコー信号にフーリエ変換処理を実施したデータを示す図である。図４は、本発明にかかる一実施形態における画像再構成部３３の構成を示すブロック図である。図５は、本発明にかかる第１の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図６は、本発明にかかる第１の実施形態における受信した複数の位相軸方向処理データのセンターリング処理の実施の動作を示すフロー図である。図７は、本発明にかかる第１の実施形態における周波数軸方向処理データを示す図である。図８は、本発明に係る第１の実施形態におけるノイズの影響のある周波数軸方向処理データを示す図である。図９は、本発明にかかる第１の実施形態において、位相軸方向処理データのセンターリング処理を示した図である。図１０は、本発明にかかる第１の実施形態において、位相軸方向処理データのセンターリング処理を示した図である。図１１は、本発明にかかる第２の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図１２は、本発明にかかる第２の実施形態における周波数軸方向処理データのプロファイルにおける積分値を示す図である。図１３は、本発明にかかる第３の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図１４は、本発明にかかる第３の実施形態における周波数軸方向処理データそれぞれの相関を示す相関図である。図１５は、本発明にかかる第３の実施形態における周波数軸方向処理データの相関係数が最大である周波数軸方向処理データの周波数差に基づいて、位相軸方向処理データのセンターリング処理を示した図である。図１６は、本発明にかかる第４の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図１７は、本発明にかかる第４の実施形態における周波数軸方向処理データの信号強度の最大値と周波数差を示す図である。図１８は、本発明にかかる第４の実施形態における全ての周波数差の平均値である変位量に基づいて、位相軸方向処理データのセンターリング処理を示した図である。図１９は、本発明にかかる第５の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図２０は、本発明にかかる第５の実施形態における周波数軸方向処理データにおけるプロファイルの積分値と周波数差を示す図である。図２１は、本発明にかかる第５の実施形態における周波数軸方向処理データにおけるプロファイルの積分値と周波数差との積の平均値に基づいて、位相軸方向処理データのセンターリング処理を示す図である。図２２は、本発明にかかる第６の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図２３は、本発明にかかる第６の実施形態における周波数軸方向処理データにおける信号強度の最大値と周波数差を示す図である。図２４は、本発明にかかる第６の実施形態における周波数軸方向処理データにおける信号強度の最大値と周波数差との積の平均値に基づいて、位相軸方向処理データのセンターリング処理を示す図である。図２５は、本発明にかかる第７の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図２６は、本発明にかかる第７の実施形態におけるそれぞれの周波数軸方向処理データを周波数軸方向の中心Ｏに変位させるための変位量の中央値に基づいて、位相軸方向処理データのセンターリング処理を示す図である。図２７は、本発明にかかる第８の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図２８は、本発明にかかる第８の実施形態における周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差の標準偏差に基づいて、選択された周波数差の平均値により位相軸方向処理データのセンターリング処理を示す図である。図２９は、従来技術を説明するための図である。図３０は、従来技術を説明するための図である。図３１は、従来技術を説明するための図である。

符号の説明

１：磁気共鳴イメージング装置、２:スキャン部、３：操作コンソール部、１１：静磁場空間、１２：静磁場マグネット部、１３：勾配コイル部、１４：ＲＦコイル部、１５：クレードル、２２：ＲＦ駆動部、２３：勾配駆動部、２４：データ収集部、３０：制御部、３１：記憶部、３２：操作部、３３：画像再構成部、３４：表示部、４０：被検体、３０１：変位量算出部、３０２：センターリング実施部、Ｃ１：表面コイル、Ｃ２：表面コイル、Ｃ３：表面コイル、Ｄ０１：表面コイルＣ１の周波数軸方向処理データ、Ｄ０２：表面コイルＣ２の周波数軸方向処理データ、Ｄ０３：表面コイルＣ３の周波数軸方向処理データ、Ｆ_Ｄ０１：周波数軸方向処理データＤ０１の変位量、Ｆ_Ｄ０２：周波数軸方向処理データＤ０２の変位量、Ｆ_Ｄ０３：周波数軸方向処理データＤ０３の変位量、Ｆ_ＡＶ１：全ての周波数軸方向処理データにおける変位量の平均値、Ｆ_ＡＶ２：全ての周波数軸方向処理データにおける最大値と変位量の積の平均値

Claims

パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、
前記複数の磁気共鳴信号のプロファイルに基づいて、当該複数の磁気共鳴信号から基準となる基準磁気共鳴信号を選択する第１ステップと、
前記第１ステップにおいて選択された前記基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数との差である周波数差を算出する第２ステップと、
前記第２ステップにおいて算出された前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップと、
前記第３ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第４ステップと
を含む、
画像生成方法。
前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて信号強度の最大値を算出し、当該最大値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する、
請求項１に記載の画像生成方法。
前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて積分値を算出し、当該積分値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する、
請求項１に記載の画像生成方法。
前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて当該複数の磁気共鳴信号同士の相関係数を算出し、前記複数の磁気共鳴信号ごとに算出された前記相関係数の平均値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する、
請求項１に記載の画像生成方法。
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数と前記ｋ空間の中心における周波数との差である周波数差を算出し、当該周波数差の平均値を算出する第１ステップと、
前記第１ステップにおいて算出された前記周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと、
前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップと
を含む、
画像生成方法。
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と前記前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性を算出し、前記周波数差と前記プロファイル特性の積の平均値を算出する第１ステップと、
前記第１ステップにおいて算出された前記周波数差と前記プロファイル特性の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと、
前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップと
を含む、
画像生成方法。
前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性は、前記複数の磁気共鳴信号のそれぞれにおける積分値である、
請求項６に記載の画像生成方法。
前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性は、前記複数の磁気共鳴信号のそれぞれにおける信号強度の最大値である、
請求項６に記載の画像生成方法。
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出する第１ステップと、
前記第１ステップにおいて算出されたそれぞれの前記周波数差の中央値を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと、
前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップと
を含む、
画像生成方法。
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の標準偏差を算出する第１ステップと、
前記第１ステップにおいて算出された前記標準偏差に基づいて、前記複数の磁気共鳴信号のそれぞれについて算出された前記周波数差から基準となる複数の基準周波数差を選択し、選択された当該複数の基準周波数差の平均値を算出する第２ステップと、
前記第２ステップにおいて算出された前記複数の基準周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップと、
前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップと
を含む、
画像生成方法。
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、
前記複数の磁気共鳴信号のプロファイルに基づいて、当該複数の磁気共鳴信号から基準となる基準磁気共鳴信号を選択する第１ステップと、
前記第１ステップにおいて選択された前記基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数との差である周波数差を算出する第２ステップと、
前記第２ステップにおいて算出された前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップと
を含む、
位置補正方法。
前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて信号強度の最大値を算出し、当該最大値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する、
請求項１１に記載の位置補正方法。
前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて積分値を算出し、当該積分値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する、
請求項１１に記載の位置補正方法。
前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて当該複数の磁気共鳴信号同士の相関係数を算出し、前記複数の磁気共鳴信号ごとに算出された前記相関係数の平均値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する、
請求項１１に記載の位置補正方法。
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数と前記ｋ空間の中心における周波数との差である周波数差を算出し、当該周波数差の平均値を算出する第１ステップと、
前記第１ステップにおいて算出された前記周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと
を含む、
位置補正方法。
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と前記前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性を算出し、前記周波数差と前記プロファイル特性の積の平均値を算出する第１ステップと、
前記第１ステップにおいて算出された前記周波数差と前記プロファイル特性の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと
を含む、
位置補正方法。
前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性は、前記複数の磁気共鳴信号のそれぞれにおける積分値である、
請求項１６に記載の位置補正方法。
前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性は、前記複数の磁気共鳴信号のそれぞれにおける信号強度の最大値である、
請求項１６に記載の位置補正方法。
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出する第１ステップと、
前記第１ステップにおいて算出されたそれぞれの前記周波数差の中央値を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと
を含む、
位置補正方法。
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の標準偏差を算出する第１ステップと、
前記第１ステップにおいて算出された前記標準偏差に基づいて、前記複数の磁気共鳴信号のそれぞれについて算出された前記周波数差から基準となる複数の基準周波数差を選択し、選択された当該複数の基準周波数差の平均値を算出する第２ステップと、
前記第２ステップにおいて算出された前記複数の基準周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップと
を含む、
位置補正方法。
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する磁気共鳴イメージング装置であって、
前記被検体にＲＦパルスを照射し、前記複数の磁気共鳴信号を受信するスキャン部と、
前記複数の磁気共鳴信号のプロファイルに基づいて、前記スキャン部により受信された前記複数の磁気共鳴信号から基準となる基準磁気共鳴信号を選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する変位量算出部と、
前記変位量算出部により算出された前記変位量に基づいて、前記複数の磁気共鳴信号それぞれを変位させ、センターリング処理を実施するセンターリング実施部と、
前記センターリング実施部により変位された前記複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像再構成部と
を有する、
磁気共鳴イメージング装置。
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度の最大値を算出し、当該最大値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて積分値を算出し、当該積分値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて当該複数の磁気共鳴信号同士の相関係数を算出し、前記複数の磁気共鳴信号ごとに算出された前記相関係数の平均値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数と前記ｋ空間の中心における周波数との差である周波数差を算出し、当該周波数差の平均値を算出し、当該周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と、前記複数の磁気共鳴信号それぞれにおける積分値を算出し、前記周波数差と前記積分値の積の平均値を算出し、当該周波数差と前記積分値の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と、前記複数の磁気共鳴信号それぞれにおける信号強度の最大値を算出し、前記周波数差と前記信号強度の最大値の積の平均値を算出し、当該周波数差と前記信号強度の最大値の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の中央値を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の標準偏差を算出し、当該標準偏差に基づいて、前記複数の磁気共鳴信号のそれぞれについて算出された前記周波数差から基準となる複数の基準周波数差を選択し、選択された当該複数の基準周波数差の平均値を算出し、当該複数の基準周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。