(57)【要約】【課題】複数の磁気共鳴信号における信号強度が最大となる部位をｋ空間における周波数軸方向の中心Ｏに変位させるための変位量に基づいて、センターリング処理を実施することにより、画像品質を向上させる。【解決手段】複数の周波数軸方向処理データにおける信号強度の最大値を算出する。そして、全ての周波数軸方向処理データにおける信号強度の最大値のうち最大値の一番大きい周波数軸方向処理データを決定し、決定した周波数軸方向処理データにおける信号強度が最大となる点がｋ空間における周波数軸方向の中心Ｏに変位されるための変位量を算出する。そして、算出された変位量に基づいて、ｋ空間において、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、ハーフエコー処理を実施して画像を生成する。【選択図】図６
2008307228.tif 000002
【特許請求の範囲】
【請求項１】
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、
前記複数の磁気共鳴信号のプロファイルに基づいて、当該複数の磁気共鳴信号から基準となる基準磁気共鳴信号を選択する第１ステップと、
前記第１ステップにおいて選択された前記基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数との差である周波数差を算出する第２ステップと、
前記第２ステップにおいて算出された前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップと、
前記第３ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第４ステップと
を含む、
画像生成方法。
【請求項２】
前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて信号強度の最大値を算出し、当該最大値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する、
請求項１に記載の画像生成方法。
【請求項３】
前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて積分値を算出し、当該積分値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する、
請求項１に記載の画像生成方法。
【請求項４】
前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて当該複数の磁気共鳴信号同士の相関係数を算出し、前記複数の磁気共鳴信号ごとに算出された前記相関係数の平均値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する、
請求項１に記載の画像生成方法。
【請求項５】
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数と前記ｋ空間の中心における周波数との差である周波数差を算出し、当該周波数差の平均値を算出する第１ステップと、
前記第１ステップにおいて算出された前記周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと、
前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップと
を含む、
画像生成方法。
【請求項６】
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と前記前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性を算出し、前記周波数差と前記プロファイル特性の積の平均値を算出する第１ステップと、
前記第１ステップにおいて算出された前記周波数差と前記プロファイル特性の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと、
前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップと
を含む、
画像生成方法。
【請求項７】
前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性は、前記複数の磁気共鳴信号のそれぞれにおける積分値である、
請求項６に記載の画像生成方法。
【請求項８】
前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性は、前記複数の磁気共鳴信号のそれぞれにおける信号強度の最大値である、
請求項６に記載の画像生成方法。
【請求項９】
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出する第１ステップと、
前記第１ステップにおいて算出されたそれぞれの前記周波数差の中央値を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと、
前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップと
を含む、
画像生成方法。
【請求項１０】
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の標準偏差を算出する第１ステップと、
前記第１ステップにおいて算出された前記標準偏差に基づいて、前記複数の磁気共鳴信号のそれぞれについて算出された前記周波数差から基準となる複数の基準周波数差を選択し、選択された当該複数の基準周波数差の平均値を算出する第２ステップと、
前記第２ステップにおいて算出された前記複数の基準周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップと、
前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップと
を含む、        
画像生成方法。
【請求項１１】
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、
前記複数の磁気共鳴信号のプロファイルに基づいて、当該複数の磁気共鳴信号から基準となる基準磁気共鳴信号を選択する第１ステップと、
前記第１ステップにおいて選択された前記基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数との差である周波数差を算出する第２ステップと、
前記第２ステップにおいて算出された前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップと
を含む、
位置補正方法。
【請求項１２】
前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて信号強度の最大値を算出し、当該最大値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する、
請求項１１に記載の位置補正方法。
【請求項１３】
前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて積分値を算出し、当該積分値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する、
請求項１１に記載の位置補正方法。
【請求項１４】
前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて当該複数の磁気共鳴信号同士の相関係数を算出し、前記複数の磁気共鳴信号ごとに算出された前記相関係数の平均値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する、
請求項１１に記載の位置補正方法。
【請求項１５】
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数と前記ｋ空間の中心における周波数との差である周波数差を算出し、当該周波数差の平均値を算出する第１ステップと、
前記第１ステップにおいて算出された前記周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと
を含む、
位置補正方法。
【請求項１６】
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と前記前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性を算出し、前記周波数差と前記プロファイル特性の積の平均値を算出する第１ステップと、
前記第１ステップにおいて算出された前記周波数差と前記プロファイル特性の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと
を含む、
位置補正方法。
【請求項１７】
前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性は、前記複数の磁気共鳴信号のそれぞれにおける積分値である、
請求項１６に記載の位置補正方法。
【請求項１８】
前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性は、前記複数の磁気共鳴信号のそれぞれにおける信号強度の最大値である、
請求項１６に記載の位置補正方法。
【請求項１９】
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出する第１ステップと、
前記第１ステップにおいて算出されたそれぞれの前記周波数差の中央値を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと
を含む、
位置補正方法。
【請求項２０】
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、
前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の標準偏差を算出する第１ステップと、
前記第１ステップにおいて算出された前記標準偏差に基づいて、前記複数の磁気共鳴信号のそれぞれについて算出された前記周波数差から基準となる複数の基準周波数差を選択し、選択された当該複数の基準周波数差の平均値を算出する第２ステップと、
前記第２ステップにおいて算出された前記複数の基準周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップと
を含む、
位置補正方法。
【請求項２１】
パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する磁気共鳴イメージング装置であって、
前記被検体にＲＦパルスを照射し、前記複数の磁気共鳴信号を受信するスキャン部と、
前記複数の磁気共鳴信号のプロファイルに基づいて、前記スキャン部により受信された前記複数の磁気共鳴信号から基準となる基準磁気共鳴信号を選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する変位量算出部と、
前記変位量算出部により算出された前記変位量に基づいて、前記複数の磁気共鳴信号それぞれを変位させ、センターリング処理を実施するセンターリング実施部と、
前記センターリング実施部により変位された前記複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像再構成部と
を有する、
磁気共鳴イメージング装置。
【請求項２２】
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度の最大値を算出し、当該最大値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。
【請求項２３】
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて積分値を算出し、当該積分値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。
【請求項２４】
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて当該複数の磁気共鳴信号同士の相関係数を算出し、前記複数の磁気共鳴信号ごとに算出された前記相関係数の平均値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。
【請求項２５】
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数と前記ｋ空間の中心における周波数との差である周波数差を算出し、当該周波数差の平均値を算出し、当該周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。
【請求項２６】
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と、前記複数の磁気共鳴信号それぞれにおける積分値を算出し、前記周波数差と前記積分値の積の平均値を算出し、当該周波数差と前記積分値の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。
【請求項２７】
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と、前記複数の磁気共鳴信号それぞれにおける信号強度の最大値を算出し、前記周波数差と前記信号強度の最大値の積の平均値を算出し、当該周波数差と前記信号強度の最大値の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。
【請求項２８】
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の中央値を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。
【請求項２９】
前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の標準偏差を算出し、当該標準偏差に基づいて、前記複数の磁気共鳴信号のそれぞれについて算出された前記周波数差から基準となる複数の基準周波数差を選択し、選択された当該複数の基準周波数差の平均値を算出し、当該複数の基準周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する、
請求項２１に記載の磁気共鳴イメージング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、磁気共鳴画像生成方法、磁気共鳴信号の位置補正方法および磁気共鳴イメージング装置に関する。
【背景技術】
【０００２】
磁気共鳴イメージング装置は、核磁気共鳴現象を利用して磁気共鳴信号を発生させ、被検体の断層画像を撮影する装置である。
【０００３】
磁気共鳴イメージング装置は、撮像時間に時間を要するため、様々な方法により撮像時間の短縮が図られている。
【０００４】
例えば、エコー時間に対するエコー波形の対称性および周波数方向におけるｋ空間データの対称性を利用して、発生するエコー信号の半分よりやや多いデータを収集し、残りの部分を共役対称性を利用して算出するハーフエコー（half echo）法や、感度分布の異なる複数のＲＦコイルからなるフェーズドアレイコイルを用いて、位相エンコードを間引いたシーケンスを実行し、行列演算によって折り返しアーチファクトを除去する展開処理を行うことにより撮影時間を短縮化するパラレルイメージング（parallel imaging）法などが知られている（例えば、特許文献１参照。）。
【０００５】
ハーフエコー法において、欠落しているデータを埋め合わせるためにホモダイン処理を行っている。ホモダイン処理を行うためには、エコー信号の中心をハイパスフィルター（High Pass Filter）およびローパスフィルター（Low Pass Filter）に通して処理を行う。このハイパスフィルターおよびローパスフィルターはｋ空間の中心にあるため、ｋ空間の中心からずれているエコー信号の中心をｋ空間の中心に変位させる必要がある。通常エコー信号のピークがエコー信号の中心であるため、エコー信号のピークをｋ空間の中心に変位させる。
【０００６】
しかし、回転磁場（Ｂ１）不均一などの影響により、フェーズドアレイコイルにより受信した複数のエコー信号におけるピークは必ずしも一致しない。例えば、図２９に示すように、フェーズドアレイコイルを構成する複数のＲＦコイルそれぞれから受信した複数のエコー信号が、エコー信号Ｅ０１、Ｅ０２およびＥ０３であるとすると、それぞれのエコー信号のピークは、ｋ空間の中心Ｏから周波数ＦＥ０１、ＦＥ０２およびＦＥ０３と異なった周波数だけずれている。そのため、受信したエコー信号それぞれのピークをｋ空間における周波数軸方向の中心Ｏに変位させるためには、それぞれのエコー信号を異なった変位量で変位させる必要がある。
【０００７】
また、パラレルイメージング法においては、ｋ空間において、全てのエコー信号が同じ変位量で変位されていないと、キャリブレーションスキャンと本スキャンにおけるコイルの感度マップの位相関係が乱れ、アーチファクトが生じる原因となる。そのため、ハーフエコー法とパラレルイメージング法を併用する場合、基準となるチャンネルのＲＦコイルを選択し、この選択されたＲＦコイルが受信したエコー信号のピークをｋ空間における周波数軸方向の中心Ｏに変位させるための変位量を、全てのエコー信号に適用することにより全てのエコー信号に対してセンターリング処理を実施している。例えば、図３０に示すように、基準と成るチャンネルのＲＦコイルが受信したエコー信号をエコー信号Ｅ０１とすると、周波数ＦＥ０１をエコー信号Ｅ０１のピークをｋ空間における周波数軸方向の中心Ｏに変位させるための変位量として、エコー信号Ｅ０１、Ｅ０２およびＥ０３を変位させ、センターリング処理を実施する。
【特許文献１】特開２００５−１９８７１５号公報（例えば、段落０００７から段落００１０参照。）
【発明の開示】
【発明が解決しようとする課題】
【０００８】
ここで、パラレルイメージング法による撮影が、頭などの局所的な部位から腹等の広範囲の部位まで及んできたため、フェーズドアレイコイルを構成する複数のＲＦコイルが多チャンネル化されている。そのため、フェーズドアレイコイルを構成する複数のＲＦコイルのそれぞれが受信するそれぞれのエコー信号における信号強度の差が、顕著になってきている。
【０００９】
図３１に示すように、信号強度が低いエコー信号は、ノイズなどの影響がある場合、エコー信号の信号強度の最大値がノイズ部分となり、エコー信号の信号強度の最大値とピークＰとが一致しない場合がある。したがって、ノイズなどの影響がある場合、信号強度の低いエコー信号を受信するＲＦコイルを基準コイルとして選択した場合、エコー信号のノイズ部分をエコー信号の最大値と判断してしまう。そのため、このエコー信号に基づいて、全てのエコー信号をセンターリング処理させるための変位量を算出し、この変位量に基づいて、全てのエコー信号を変位させると、全てのエコー信号において、ｋ空間における周波数軸方向の中心Ｏにエコー信号の実際のピークＰが変位されず、実際のピークＰからずれた部分がｋ空間における周波数軸方向の中心Ｏに変位されることとなる。そのため、生成される画像にアーチファクトが生じることがある。
【００１０】
したがって、本発明は、ハーフエコー法とパラレルイメージング法を併用した際に、生成した画像にアーチファクトが生じないようなエコー信号の位置補正方法、画像生成方法および磁気イメージング装置を提供する。
【課題を解決するための手段】
【００１１】
本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、前記複数の磁気共鳴信号のプロファイルに基づいて、当該複数の磁気共鳴信号から基準となる基準磁気共鳴信号を選択する第１ステップと、前記第１ステップにおいて選択された前記基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数との差である周波数差を算出する第２ステップと、前記第２ステップにおいて算出された前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップと、前記第３ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第４ステップとを含む。
【００１２】
好適には、前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて信号強度の最大値を算出し、当該最大値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する。
【００１３】
好適には、前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて積分値を算出し、当該積分値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する。
【００１４】
好適には、前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて当該複数の磁気共鳴信号同士の相関係数を算出し、前記複数の磁気共鳴信号ごとに算出された前記相関係数の平均値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する。
【００１５】
本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数と前記ｋ空間の中心における周波数との差である周波数差を算出し、当該周波数差の平均値を算出する第１ステップと、前記第１ステップにおいて算出された前記周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと、前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップとを含む。
【００１６】
本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と前記前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性を算出し、前記周波数差と前記プロファイル特性の積の平均値を算出する第１ステップと、前記第１ステップにおいて算出された前記周波数差と前記プロファイル特性の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと、前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップとを含む。
【００１７】
好適には、前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性は、前記複数の磁気共鳴信号のそれぞれにおける積分値である。
【００１８】
好適には、前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性は、前記複数の磁気共鳴信号のそれぞれにおける信号強度の最大値である。
【００１９】
本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出する第１ステップと、前記第１ステップにおいて算出されたそれぞれの前記周波数差の中央値を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップと、前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップとを含む。
【００２０】
本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する画像生成方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の標準偏差を算出する第１ステップと、前記第１ステップにおいて算出された前記標準偏差に基づいて、前記複数の磁気共鳴信号のそれぞれについて算出された前記周波数差から基準となる複数の基準周波数差を選択し、選択された当該複数の基準周波数差の平均値を算出する第２ステップと、前記第２ステップにおいて算出された前記複数の基準周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップと、前記第２ステップにおいて変位された前記複数の磁気共鳴信号についてハーフエコー処理を施す第３ステップとを含む。
【００２１】
本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、前記複数の磁気共鳴信号のプロファイルに基づいて、当該複数の磁気共鳴信号から基準となる基準磁気共鳴信号を選択する第１ステップと、前記第１ステップにおいて選択された前記基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数との差である周波数差を算出する第２ステップと、前記第２ステップにおいて算出された前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップとを含む。
【００２２】
好適には、前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて信号強度の最大値を算出し、当該最大値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する。
【００２３】
好適には、前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて積分値を算出し、当該積分値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する。
【００２４】
好適には、前記第１ステップは、前記複数の磁気共鳴信号のそれぞれについて当該複数の磁気共鳴信号同士の相関係数を算出し、前記複数の磁気共鳴信号ごとに算出された前記相関係数の平均値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択する。
【００２５】
本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数と前記ｋ空間の中心における周波数との差である周波数差を算出し、当該周波数差の平均値を算出する第１ステップと、前記第１ステップにおいて算出された前記周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップとを含む。
【００２６】
本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と前記前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性を算出し、前記周波数差と前記プロファイル特性の積の平均値を算出する第１ステップと、前記第１ステップにおいて算出された前記周波数差と前記プロファイル特性の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップとを含む。
【００２７】
好適には、前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性は、前記複数の磁気共鳴信号のそれぞれにおける積分値である。
【００２８】
好適には、前記複数の磁気共鳴信号のそれぞれにおけるプロファイル特性は、前記複数の磁気共鳴信号のそれぞれにおける信号強度の最大値である。
【００２９】
本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出する第１ステップと、前記第１ステップにおいて算出されたそれぞれの前記周波数差の中央値を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第２ステップとを含む。
【００３０】
本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像生成方法における、前記複数の磁気共鳴信号のそれぞれを変位させる位置補正方法であって、前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の標準偏差を算出する第１ステップと、前記第１ステップにおいて算出された前記標準偏差に基づいて、前記複数の磁気共鳴信号のそれぞれについて算出された前記周波数差から基準となる複数の基準周波数差を選択し、選択された当該複数の基準周波数差の平均値を算出する第２ステップと、前記第２ステップにおいて算出された前記複数の基準周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定し、当該変位量に基づいて前記複数の磁気共鳴信号それぞれを変位させる第３ステップとを含む。
【００３１】
本発明は、パラレルイメージング法によって、複数のＲＦコイルによって静磁場空間内の被検体をスキャンすることにより、前記複数のＲＦコイルのそれぞれが受信して得られるそれぞれの複数の磁気共鳴信号に基づいて、画像を生成する磁気共鳴イメージング装置であって、前記被検体にＲＦパルスを照射し、前記複数の磁気共鳴信号を受信するスキャン部と、前記複数の磁気共鳴信号のプロファイルに基づいて、前記スキャン部により受信された前記複数の磁気共鳴信号から基準となる基準磁気共鳴信号を選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する変位量算出部と、前記変位量算出部により算出された前記変位量に基づいて、前記複数の磁気共鳴信号それぞれを変位させ、センターリング処理を実施するセンターリング実施部と、前記センターリング実施部により変位された前記複数の磁気共鳴信号にハーフエコー処理を実施して、画像を生成する画像再構成部とを有する。
【００３２】
好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度の最大値を算出し、当該最大値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する。
【００３３】
好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて積分値を算出し、当該積分値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する。
【００３４】
好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて当該複数の磁気共鳴信号同士の相関係数を算出し、前記複数の磁気共鳴信号ごとに算出された前記相関係数の平均値が全ての前記複数の磁気共鳴信号の中で最大となる磁気共鳴信号を前記基準磁気共鳴信号として選択し、選択された当該基準磁気共鳴信号における信号強度が最大となる周波数とｋ空間の中心における周波数の差である周波数差を算出し、前記周波数差を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する。
【００３５】
好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数と前記ｋ空間の中心における周波数との差である周波数差を算出し、当該周波数差の平均値を算出し、当該周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する。
【００３６】
好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と、前記複数の磁気共鳴信号それぞれにおける積分値を算出し、前記周波数差と前記積分値の積の平均値を算出し、当該周波数差と前記積分値の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する。
【００３７】
好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差と、前記複数の磁気共鳴信号それぞれにおける信号強度の最大値を算出し、前記周波数差と前記信号強度の最大値の積の平均値を算出し、当該周波数差と前記信号強度の最大値の積の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する。
【００３８】
好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の中央値を変位させる変位量として設定する。
【００３９】
好適には、前記変位量算出部は、前記スキャン部により受信された前記複数の磁気共鳴信号のそれぞれについて信号強度が最大である周波数とｋ空間の中心における周波数との周波数差を算出し、当該周波数差の標準偏差を算出し、当該標準偏差に基づいて、前記複数の磁気共鳴信号のそれぞれについて算出された前記周波数差から基準となる複数の基準周波数差を選択し、選択された当該複数の基準周波数差の平均値を算出し、当該複数の基準周波数差の平均値を全ての前記複数の磁気共鳴信号を変位させる変位量として設定する。
【発明の効果】
【００４０】
本発明によれば、パラレルイメージング法とハーフエコー法を併用した際に、生成した画像にアーチファクトが生じないような画像生成方法、エコー信号の位置補正方法、および磁気イメージング装置を提供することができる。
【発明を実施するための最良の形態】
【００４１】
以下より、本発明にかかる一実施形態について図面を参照して説明する。
【００４２】
＜第１の実施形態＞
（装置構成）
図１は、本発明にかかる一実施形態におけるＲＦコイル部により構成される磁気共鳴イメージング装置の構成を示す構成図である。本装置は、本発明の実施形態の一例である。
【００４３】
図１に示すように、磁気共鳴イメージング装置１は、スキャン部２と、操作コンソール部３とを有している。ここで、スキャン部２は、静磁場マグネット部１２と、勾配コイル部１３と、ＲＦコイル部１４と、クレードル１５とを有する。そして、操作コンソール部３は、ＲＦ駆動部２２と、勾配駆動部２３と、データ収集部２４と、制御部３０と、記憶部３１と、操作部３２と、画像再構成部３３と、表示部３４とを有する。
【００４４】
スキャン部２について説明する。
【００４５】
スキャン部２は、図１に示すように、被検体４０において体動する撮影対象を含む撮影スライス領域が収容される静磁場空間１１を含んでいる。そして、スキャン部２は、操作コンソール部３からの制御信号に基づいて、その静磁場が形成される静磁場空間１１に収容した被検体４０の撮影領域にＲＦパルスを照射し、その撮影領域から生ずる磁気共鳴信号をイメージングデータとして収集するスキャンを被検体４０について実施する。
【００４６】
スキャン部２の各構成要素について、順次、説明する。図２は、フェーズドアレイコイルの一例を示す図である。理解促進の便宜上、表面コイルの重なり合いをずらして記載しており、またコイル線の太さを変更している。図３は、フェーズドアレイコイルにより受信したエコー信号にフーリエ変換処理を実施したデータを示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。周波数０の破線はｋ空間の中心Ｏを示す。
【００４７】
静磁場マグネット部１２は、被検体４０が収容される静磁場空間１１に静磁場を形成するために設けられている。静磁場マグネット部１２は、水平磁場型であって、被検体４０が収容される静磁場空間１１において載置される被検体４０の体軸方向（ｚ方向）に沿うように、超伝導磁石（図示なし）が静磁場を形成する。なお、静磁場マグネット部１２は、水平磁場型の他に、垂直磁場型であってもよく、また静磁場マグネット部１２は、永久磁石により構成されていてもよい。
【００４８】
勾配コイル部１３は、ＲＦコイル部１４が受信する磁気共鳴信号に３次元の位置情報を持たせるために、静磁場空間１１に勾配磁場を形成する。勾配コイル部１３は、スライス選択勾配磁場、読み取り勾配磁場、位相エンコード勾配磁場の３種類の勾配磁場を形成するために勾配コイルを３系統有する。
【００４９】
ＲＦコイル部１４は、たとえば、被検体４０を囲むように配置される。ＲＦコイル部１４は、静磁場マグネット部１２によって静磁場が形成される静磁場空間１１内において、制御部３０からの制御信号に基づいて、電磁波であるＲＦパルスを被検体４０に送信して高周波磁場を形成する。これにより、被検体４０の撮影スライス領域におけるプロトンのスピンを励起する。そして、ＲＦコイル部１４は、その励起された被検体４０の撮影スライス領域におけるプロトンのスピンが元の磁化ベクトルへ戻る際に生ずる電磁波を、磁気共鳴信号として受信する。
【００５０】
また、ＲＦコイル部１４は、例えば、図２に示す複数の独立した表面コイルからなるフェーズドアレイコイル１４１である。この複数の独立した表面コイルは、それぞれ個別に送信機および受信機と接続されており、多チャンネルで複数の磁気共鳴信号を受信するように構成されている。
【００５１】
図２に示すように、例えば、フェーズドアレイコイル１４１は、独立した表面コイルＣ１、Ｃ２およびＣ３からなる。そして、図３に示すように、表面コイルＣ１、Ｃ２およびＣ３により受信されたそれぞれのエコー信号についてフーリエ変換処理を実施したデータ（以下、フーリエ変換処理データとも称する）は、異なるｋ空間に表される。また、回転磁場不均一の度合いが場所によって異なるため、フーリエ変換処理データにおけるピークは、ｋ空間の中心Ｏからそれぞれ異なる周波数だけずれることとなる。例えば、図３（ａ）に示すように、ｋ空間に表された表面コイルＣ１のフーリエ変換処理データＤ１におけるピークは、ｋ空間の中心Ｏから周波数ＦＤ1だけずれている。同様に図３（ｂ）、図３（ｃ）に示すように、フーリエ変換処理データＤ２、Ｄ３におけるピークは、ｋ空間の中心Ｏからそれぞれ周波数ＦＤ２、ＦＤ３だけずれており、それぞれｋ空間の中心Ｏから異なる周波数ずれている。
【００５２】
クレードル１５は、被検体４０を載置するテーブルを有する。クレードル１５は、制御部３０からの制御信号に基づいて、静磁場空間１１の内部と外部との間でテーブルに載置された被検体４０を移動する。
【００５３】
操作コンソール部３について説明する。
【００５４】
操作コンソール部３は、スキャン部２が被検体４０についてスキャンを実施するように制御し、そのスキャン部２が実施したスキャンによって得られた磁気共鳴信号に基づいて、被検体４０の画像を生成すると共に、その生成した画像を表示する。
【００５５】
操作コンソール部３を構成する各部について、順次、説明する。
【００５６】
ＲＦ駆動部２２は、ＲＦコイル部１４を駆動させて静磁場空間１１内に高周波磁場を形成するために、ゲート変調器（図示なし）とＲＦ電力増幅器（図示なし）とＲＦ発振器（図示なし）とを有する。ＲＦ駆動部２２は、制御部３０からの制御信号に基づいて、ＲＦ発振器からのＲＦ信号を、ゲート変調器を用いて所定のタイミングおよび所定の包絡線の信号に変調する。そして、ゲート変調器により変調されたＲＦ信号を、ＲＦ電力増幅器により増幅した後、ＲＦコイル部１４に出力する。
【００５７】
勾配駆動部２３は、制御部３０の制御信号に基づいて勾配コイル部１３を駆動させて、静磁場空間１１内に勾配磁場を発生させる。勾配駆動部２３は、勾配コイル部１３の３系統の勾配コイルに対応して３系統の駆動回路（図示なし）を有する。
【００５８】
データ収集部２４は、ＲＦコイル部１４が受信する磁気共鳴信号を収集するために、位相検波器（図示なし）とアナログ／デジタル変換器（図示なし）とを有する。データ収集部２４は、ＲＦコイル部１４からの磁気共鳴信号を、ＲＦ駆動部２２のＲＦ発振器の出力を参照信号として、位相検波器によって位相検波し、アナログ／デジタル変換器に出力する。そして、位相検波器により位相検波されたアナログ信号である磁気共鳴信号を、アナログ／デジタル変換器によってデジタル信号に変換して、画像再構成部３３に出力する。
【００５９】
制御部３０は、コンピュータと、コンピュータを用いて所定のスキャンに対応する動作を各部に実行させるプログラムとを有する。そして、制御部３０は、操作部３２に接続されており、操作部３２に入力された操作信号を処理し、クレードル１５とＲＦ駆動部２２と勾配駆動部２３とデータ収集部２４との各部に、制御信号を出力し制御を行う。また、制御部３０は、所望の画像を得るために、操作部３２からの操作信号に基づいて画像再構成部３３を制御する。
【００６０】
記憶部３１は、コンピュータにより構成されている。そして、記憶部３１は、データ収集部２４に収集された画像再構成処理前の磁気共鳴信号、画像再構成部３３で画像再構成処理された画像データ等を記憶する。
【００６１】
操作部３２は、キーボードやマウスなどの操作デバイスにより構成されている。操作部３２は、オペレータによって操作データ、撮像プロトコルなどが入力され、その操作データ、撮像プロトコルを制御部３０に出力する。
【００６２】
画像再構成部３３は、コンピュータにより構成されている。そして、画像再構成部３３は、データ収集部２４に接続されており、データ収集部２４から出力される磁気共鳴信号に対して画像再構成処理を実施して、画像を生成する。
【００６３】
図４は、本発明にかかる一実施形態における画像再構成部３３の構成を示すブロック図である。
【００６４】
図４に示すように、画像再構成部３３は、変位量算出部３０１と、センターリング実施部３０２とを有している。詳細については後述するが、制御部３０は、ＲＦコイル部１４により受信したエコー信号をｋ空間における周波数軸方向の中心Ｏにセンターリングさせるために、エコー信号を変位させる変位量を算出し、算出した変位量に基づいてセンターリング処理を施す。
【００６５】
画像再構成部３３の各部について順次説明する。
【００６６】
変位量算出部３０１は、複数の表面コイルにより受信された複数のエコー信号に基づき、エコー信号のピークをｋ空間における周波数軸方向の中心Ｏにセンターリングさせるためのエコー信号の変位量を算出する。
【００６７】
センターリング実施部３０２は、変位量算出部３０１が算出した変位量に基づき、受信した全てのエコー信号を変位させ、センターリング処理を実施する。
【００６８】
表示部３４は、ディスプレイなどの表示デバイスにより構成されており、画像再構成部３３が生成する被検体４０の画像を表示する。
【００６９】
（動作）
以下より、本実施形態の磁気共鳴イメージング装置１を用いて、被検体４０を撮像する際の動作について説明する。
【００７０】
図５は、本発明にかかる第１の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。
【００７１】
まず、図５に示すように、スキャンを実施する（ＳＴ１０）。
【００７２】
ここでは、スキャン部２が被検体４０のスキャンを実施する。
【００７３】
ここでは、オペレータによって操作部３２に入力された操作信号に基づいて、制御部３０は、ＲＦ駆動部２２と勾配駆動部２３とデータ収集部２４とにそれぞれ制御信号を出力し制御を行って、被検体４０をスキャンし、データ収集部２４にエコー信号を取得させる。そして、そのデータ収集部２４が取得したエコー信号を画像再構成部３３に出力させる。
【００７４】
つぎに、図５に示すように、画像の生成を実施する（ＳＴ２０）。
【００７５】
ここでは、データ収集部２４から出力されるエコー信号に対して画像再構成処理を実施して、画像を生成する。画像再構成処理の際、ｋ空間において、受信した複数のエコー信号のセンターリング処理を実施する。
【００７６】
図６は、本発明にかかる第１の実施形態における受信した複数の位相軸方向処理データのセンターリング処理の実施の動作を示すフロー図である。図７は、本発明にかかる第１の実施形態における周波数軸方向処理データを示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。独立した表面コイルのそれぞれにより受信されたエコー信号について高速フーリエ変換（以下、ＦＦＴとも称する）処理を実施したデータは異なるｋ空間に表されるが、理解補助の便宜上、同一のｋ空間に表している。また、理解補助の便宜上、周波数軸方向処理データＤ０１を太線、周波数軸方向処理データＤ０２を中太線、周波数軸方向処理データＤ０３を細線で記載している。図８は、本発明に係る第１の実施形態におけるノイズの影響のある周波数軸方向処理データを示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。図９および図１０は、本発明にかかる第１の実施形態において、位相軸方向処理データのセンターリング処理を示した図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。また、周波数０の破線は、ｋ空間における周波数軸方向の中心Ｏを示す。また、変位量ＦＤ０１は、周波数軸方向処理データＤ０１の信号強度が最大となる点をｋ空間における周波数軸方向の中心Ｏに変位させるための変位量である。また、太線はホモダインハイパスフィルターＨＦ、ホモダインローパスフィルターＬＦを示し、理解促進の便宜上、ホモダインハイパスフィルターＨＦとホモダインローパスフィルターＬＦの重なり合いをずらして記載している（以下の図で同じ）。
【００７７】
図６に示すように、まず、iチャンネル目のエコー信号に対し、周波数方向にＦＦＴ処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０１１）。
【００７８】
ここでは、iチャンネル目の表面コイルＣｉにより受信されたエコー信号に対して周波数軸方向にＦＦＴ処理を実施する。
【００７９】
具体的には、画像再構成部３３が、iチャンネル目の表面コイルＣｉにより受信されたエコー信号に対して周波数軸方向にＦＦＴ処理を実施する。そして、エコー信号に対して周波数軸方向にＦＦＴ処理を実施することにより算出されたデータ（以下、周波数軸方向処理データとも称する）を記憶部３１が保存する。
【００８０】
次に、図６に示すように、信号強度の最大値を算出する（ＳＴ２０２１）。
【００８１】
ここでは、ステップＳＴ２０１１において算出された周波数軸方向処理データの信号強度の最大値を算出する。
【００８２】
具体的には、ステップＳＴ２０１１において、画像再構成部３３により算出された周波数軸方向処理データの信号強度の最大値を算出する。図７に示すように、iチャンネル目の表面コイルＣｉにより受信されたエコー信号に対して周波数軸方向にＦＦＴ処理を実施したデータが周波数軸方向処理データＤ０１の場合には、周波数軸方向処理データＤ０１における信号強度の最大値として、ピーク値ＰＤ０１が算出される。また、データが周波数軸方向処理データＤ０２の場合、信号強度の最大値としてピーク値ＰＤ０２が算出され、データが周波数軸方向処理データＤ０３の場合、信号強度の最大値としてピーク値ＰＤ０３が算出される。図７の場合、それぞれの周波数軸方向処理データにおける信号強度の最大値はピーク値となっているが、例えば、図８に示すように、ノイズなどの影響のある周波数軸方向処理データの場合、信号強度の最大値がピーク値となるとは限らないが、この場合、信号強度の最大値として、ピーク値ＰＤ０４ではなく、ＮＤ０４が算出される。
【００８３】
次に、図６に示すように、位相軸方向にＦＦＴ処理を実施する（ＳＴ２０３１）。
【００８４】
ここでは、ステップＳＴ２０１１において周波数軸方向にＦＦＴ処理を実施すことにより算出されたデータに対して位相軸方向にＦＦＴ処理を実施する。
【００８５】
具体的には、ステップＳＴ２０１１において、画像再構成部３３により算出され、記憶部３１に保存されている周波数軸方向処理データに対して、位相軸方向にＦＦＴ処理を実施する。そして、位相軸方向にＦＦＴ処理を実施することにより算出されたデータ（以下、位相軸方向処理データとも称する）を記憶部３１が保存する。
【００８６】
次に、図６に示すように、ステップＳＴ２０３１において算出された位相軸方向処理データが、最後のチャンネルの表面コイルにより受信されたエコー信号から算出されたデータであるかの判断を行う（ＳＴ２０４１）。
【００８７】
ここでは、ステップＳＴ２０３１において算出された位相軸方向処理データが、最後のチャンネルの表面コイルにより受信されたエコー信号から算出されたデータであるかを、制御部３０が判断する。そして、位相軸方向処理データが最後のチャンネルの表面コイルにより受信されたエコー信号から算出されたデータでない場合（Ｎｏ）には、上記の周波数軸方向のＦＦＴ処理、信号強度の最大値の算出および位相軸方向のＦＦＴ処理を継続するように制御部３０が各部を制御する。
【００８８】
次に、図６に示すように、全てのチャンネルにおける信号強度の最大値を比較し、最大値の一番大きいチャンネルの変位量を算出する（ＳＴ２０５１）。
【００８９】
ここでは、ステップＳＴ２０２１において全てのチャンネルの表面コイルにより受信されたエコー信号から算出された周波数軸方向処理データにおける信号強度の最大値を比較する。そして、変位量算出部３０１が、ｋ空間において信号強度の最大値の一番大きいチャンネルの周波数軸方向処理データを基準周波数軸方向処理データとして、信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。
【００９０】
具体的には、図７に示すように、変位量算出部３０１が、周波数軸方向処理データＤ０１、Ｄ０２およびＤ０３のそれぞれの信号強度の最大値であるピーク値ＰＤ０１、ピーク値ＰＤ０２およびピーク値ＰＤ０３を比較し、最大値の一番大きいデータである基準周波数軸方向処理データを周波数軸方向処理データＤ０１と決定する。ここで、図８に示す周波数軸方向処理データＤ０４のように、ピークにおける信号強度よりもノイズ部分における信号強度の方が高い周波数軸方向処理データが複数の周波数軸方向処理データの中に含まれており、ノイズ部分における信号強度ＮＤ０４が全ての周波数軸方向処理データにおける信号強度の最大値の中で一番大きい場合には、周波数軸方向処理データＤ０４が選択されることとなる。しかし、ノイズ部分における信号強度がピークにおける信号強度よりも高くなる周波数軸方向処理データは、周波数軸方向処理データ全体として信号強度が低い。そのため、このような周波数軸方向処理データが選択されることはない。
【００９１】
そして、図９に示すように、変位量算出部３０１が、ｋ空間において周波数軸方向処理データＤ０１における信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、当該周波数差を変位量ＦＤ０１とし、全ての位相軸方向処理データを変位させるための変位量として設定する。ここで、周波数軸方向処理データと位相軸方向処理データのピーク位置は変わらないため、位相軸方向処理データを変位させる変位量として、変位量ＦＤ０１を用いることができる。
【００９２】
次に、図６に示すように、ｉチャンネル目のデータのセンターリング処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０６１）。
【００９３】
ここでは、ステップＳＴ２０５１において変位量算出部３０１が算出した変位量に基づいて、ステップＳＴ２０３１において画像再構成部３３が算出したiチャンネル目の表面コイルＣｉの位相軸方向処理データに対するセンターリング処理をセンターリング実施部３０２が実施する。
【００９４】
具体的には、例えば、iチャンネル目の表面コイルＣｉが表面コイルＣ１である場合には、ｋ空間において、図１０（ａ）に示す位相軸方向処理データＤ１１を、図１０（ｂ）に示すように変位量ＦＤ０１だけ変位させるセンターリング処理をセンターリング実施部３０２が実施する。ステップＳＴ２０５１において算出された変位量ＦＤ０１が、位相軸方向処理データＤ１１をｋ空間における周波数軸方向の中心Ｏに変位させるための変位量であるため、位相軸方向処理データＤ１１を変位量ＦＤ０１だけ変位すると、位相軸方向処理データＤ１１の信号強度が最大となる点（ピーク）は、ｋ空間における周波数軸方向の中心Ｏに変位される。そして、ｋ空間における周波数軸方向の中心Ｏに変位された位相軸方向処理データＤ１１のピークがホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１１に適切にフィルターがかかることとなる。
【００９５】
また、図１０（ｃ）に示すように、ｋ空間において、位相軸方向処理データＤ１２におけるピークは、位相軸方向処理データＤ１１のピークと異なる周波数にあるため、位相軸方向処理データＤ１２を変位量ＦＤ０１に基づいて変位させると、位相軸方向処理データＤ１２のピークはｋ空間における周波数軸方向の中心Ｏに変位されない。しかし、図１０（ｄ）に示すように、変位された位相軸方向処理データＤ１２のピークがｋ空間における周波数軸方向の中心ＯにあるホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１２に適切にフィルターがかかることとなる。したがって、アーチファクトは殆ど生じない。
【００９６】
また、位相軸方向処理データＤ１３についても同様に変位量ＦＤ０１に基づいて、センターリング実施部３０２がセンターリング処理を実施する。位相軸方向処理データＤ１３も位相軸方向処理データＤ１２と同様に、変位された位相軸方向処理データＤ１３のピークがｋ空間における周波数軸方向の中心ＯにあるホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１３に適切にフィルターがかかることとなる。したがって、アーチファクトは殆ど生じない。
【００９７】
次に、図６に示すように、ハーフエコー処理を実施する（ＳＴ２０７１）。
【００９８】
ここでは、ＳＴ２０６１においてセンターリング処理を実施した位相軸方向処理データに対して、画像再構成部３３がハーフエコー処理を実施する。
【００９９】
次に、図６に示すように、ハーフエコー処理を実施したデータが最後のチャンネルのデータであるかの判断を行う（ＳＴ２０８１）。
【０１００】
ここでは、ＳＴ２０７１においてハーフエコー処理を実施することにより算出されたデータ（以下、ハーフエコー処理とも称する）が、最後のチャンネルにおける表面コイルにより受信されたエコー信号から算出されたかを、制御部３０が判断する。そして、ハーフエコー処理データが最後のチャンネルにおける表面コイルにより受信されたエコー信号から算出されたデータでない場合（Ｎｏ）には、上記のセンターリング処理およびハーフエコー処理を継続するように制御部３０が各部を制御する。
【０１０１】
以上のように、本発明の第１の実施形態は、複数の表面コイルにより受信された全てのエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出されたそれぞれの周波数軸方向処理データにおける信号強度の最大値を算出する。そして、全ての周波数軸方向処理データにおける信号強度の最大値の中で、最大値の一番大きい基準周波数軸方向処理データを決定し、決定した基準周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、当該周波数差を全ての周波数軸方向処理データを変位させるための変位量と設定する。そして、算出された変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。
【０１０２】
このように、最大値がピークとならない周波数軸方向処理データ、すなわちノイズ部分の方が実際のピークよりも信号強度が高くなる周波数軸方向処理データは、周波数軸方向処理データ全体として信号強度が低いため、信号強度の最大値が一番大きい基準周波数軸方向処理データとして選択されることはない。したがって、本発明にかかる第１の実施形態は、信号強度の最大値が一番大きい周波数軸方向処理データにおける信号強度が最大となる点をｋ空間における周波数軸方向の中心Ｏに変位させるための変位量に基づいて、全ての位相軸方向処理データをｋ空間において変位させることにより、生成した画像のアーチファクトを防止することができ、画像品質を向上させることができる。
【０１０３】
＜第２の実施形態＞
以下より本発明にかかる第２の実施形態について図１１により説明する。
第２の実施形態は、動作フローにおけるステップＳＴ２０２２である積分値の算出、およびステップＳＴ２０５２である変位量の算出以外の部分は第１の実施形態と同じである。そのため、重複する箇所については、記載を省略する。
【０１０４】
（動作）
以下より、本実施形態の磁気共鳴イメージング装置１を用いて、被検体４０を撮像する際の動作について説明する。
【０１０５】
図１１は、本発明にかかる第２の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図１２は、本発明にかかる第２の実施形態における周波数軸方向処理データのプロファイルの積分値を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。また、周波数軸方向処理データのプロファイルと横軸で囲まれた斜線部が積分値であり、ＳＤ０１、ＳＤ０２およびＳＤ０３はそれぞれ周波数軸方向処理データＤ０１、Ｄ０２およびＤ０３の積分値を示す。
【０１０６】
まず、図１１に示すように、積分値を算出する（ＳＴ２０２２）。
【０１０７】
ここでは、ステップＳＴ２０１１において、変位量算出部３０１が周波数軸方向のＦＦＴ処理を実施することにより算出されたデータのｋ空間における積分値を算出する。
【０１０８】
具体的には、ステップＳＴ２０１１において、変位量算出部３０１が画像再構成部３３により算出された周波数軸方向処理データのプロファイルの積分値を算出する。図１２（ａ）に示すように、iチャンネル目の表面コイルＣｉにより受信されたエコー信号に対して周波数軸方向にＦＦＴ処理を実施することにより算出されたデータが周波数軸方向処理データＤ０１の場合には、周波数軸方向処理データＤ０１のプロファイルの積分値ＳＤ０１を算出する。積分値ＳＤ０１は、周波数軸方向処理データＤ０１のデータプロファイルと横軸とで囲まれている斜線部分の面積である。また、算出されたデータが周波数軸方向処理データＤ０２、Ｄ０３の場合、積分値として、積分値ＳＤ０２、積分値ＳＤ０３のそれぞれが算出される。
【０１０９】
次に、図１１に示すように、全てのチャンネルにおける積分値を比較し、積分値の一番大きいチャンネルにおける変位量を算出する（ＳＴ２０５２）。
【０１１０】
ここでは、ＳＴ２０２２において算出した全ての周波数軸方向処理データのプロファイルの積分値を比較し、変位量算出部３０１が、ｋ空間において積分値の一番大きいチャンネルの周波数軸方向処理データを基準周波数軸方向処理データとして、信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。
【０１１１】
具体的には、図１２に示すように、周波数軸方向処理データＤ０１、Ｄ０２、Ｄ０３のプロファイルのそれぞれの積分値ＳＤ０１、積分値ＳＤ０２、積分値ＳＤ０３を比較し、積分値の一番大きいデータである基準周波数軸方向処理データを周波数軸方向処理データＤ０１と決定する。そして、図９に示すように、変位量算出部３０１が、ｋ空間において周波数軸方向処理データＤ０１の信号強度が最大となる点における周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、当該周波数差を変位量ＦＤ０１とし、全ての位相軸方向処理データを変位させるための変位量として設定する。ここで、周波数軸方向処理データと位相軸方向処理データのピーク位置は変わらないため、位相軸方向処理データを変位させる変位量として、変位量ＦＤ０１を用いることができる。
【０１１２】
以上のように、本発明の第２の実施形態は、複数の表面コイルにより受信した全てのエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出されたそれぞれの周波数軸方向処理データのプロファイルの積分値を算出する。そして、全ての周波数軸方向処理データのプロファイルの積分値から積分値の一番大きい基準周波数軸方向処理データを決定し、決定した周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、当該周波数差を全ての周波数軸方向処理データを変位させるための変位量として設定する。そして、算出された変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。
【０１１３】
このように、周波数軸方向処理データの最大値がピークとならない場合、すなわち、ノイズ部分の方が実際のピークよりも信号強度が高くなる周波数軸方向処理データは、周波数軸方向処理データのプロファイルの積分値が小さいため、本実施形態における基準周波数軸方向処理データとして選択されることはない。したがって、本発明にかかる第２の実施形態は、周波数軸方向処理データのプロファイルの積分値が一番大きい周波数軸方向処理データの信号強度が最大となる点をｋ空間における周波数軸方向の中心Ｏに変位させるための変位量に基づいて、ｋ空間において全ての位相軸方向処理データを変位させることにより、生成した画像のアーチファクトを防止することができ、画像品質を向上させることができる。
【０１１４】
＜第３の実施形態＞
以下より本発明にかかる第３の実施形態について図１３により説明する。
第３の実施形態は、動作フローにおけるステップＳＴ２０５３である変位量の算出および、ステップＳＴ２０６３であるセンターリング処理の実施以外の部分は第１の実施形態と同じである。そのため、重複する箇所については、記載を省略する。
【０１１５】
（動作）
以下より、本実施形態の磁気共鳴イメージング装置１を用いて、被検体４０を撮像する際の動作について説明する。
【０１１６】
図１３は、本発明にかかる第３の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図１４は、本発明にかかる第３の実施形態における周波数軸方向処理データそれぞれの相関を示す相関図である。縦軸は信号強度であり、横軸は周波数である。理解補助の便宜上、周波数軸方向処理データＤ０１を太線、周波数軸方向処理データＤ０２を中太線、周波数軸方向処理データＤ０３を細線で記載している。図１５は、本発明にかかる第３の実施形態における周波数軸方向処理データの相関係数が最大である周波数軸方向処理データの周波数差に基づいて、位相軸方向処理データのセンターリング処理を示した図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。
【０１１７】
まず、図１３に示すように、各チャンネルにおけるデータの相関をとり、相関係数の一番大きいチャンネルの変位量を算出する（ＳＴ２０５３）。
【０１１８】
ここでは、変位量算出部３０１が、各チャンネルにおける周波数軸方向処理データごとに他のチャンネルにおける周波数軸方向処理データとの相関をとり、相関係数を算出する。そして、変位量算出部３０１が、ｋ空間において相関係数が一番大きいチャンネルの周波数軸方向処理データを基準周波数軸方向処理データとして、信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。
【０１１９】
具体的には、例えば、図１４（ａ）に示すように、周波数軸方向処理データＤ０１と周波数軸方向処理データＤ０２との相関をとり、周波数軸方向処理データＤ０１と周波数軸方向処理データＤ０２との相関係数Ｒ１２を変位量算出部３０１が算出する。そして、図１４（ｂ）および図１４（ｃ）に示すように、周波数軸方向処理データＤ０１と周波数軸方向処理データＤ０３との相関係数Ｒ１３および周波数軸方向処理データＤ０２と周波数軸方向処理データＤ０３との相関係数Ｒ２３も変位量算出部３０１が算出する。相関係数は、数式（１）により算出する。
【０１２０】
【数１】
2008307228.tif 000003
【０１２１】
次に、変位量算出部３０１が、それぞれの周波数軸方向処理データにおける相関係数の平均値を算出する。例えば、周波数軸方向処理データＤ０１における相関係数の平均値は（Ｒ１２＋Ｒ１３）／２である。同様に周波数軸方向処理データＤ０２における相関係数の平均値は（Ｒ１２＋Ｒ2３）／２であり、周波数軸方向処理データＤ０３における相関係数の平均値は（Ｒ１３＋Ｒ2３）／２である。そして、変位量算出部３０１が、それぞれの平均値を比較し、平均値の一番大きいデータである基準周波数軸方向処理データを決定する。そして、平均値の一番大きい周波数軸方向処理データにおける、信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を、変位量算出部３０１が算出する。そして、当該周波数差を全ての位相軸方向処理データを変位させるための変位量として設定する。本実施形態の場合、相関係数の一番大きいデータである基準周波数軸方向処理データは、周波数軸方向処理データＤ０３である。したがって、設定される変位量は、変位量ＦＤ０３となる。ここで、周波数軸方向処理データと位相軸方向処理データのピーク位置は変わらないため、位相軸方向処理データを変位させる変位量として、変位量ＦＤ０３を用いることができる。
【０１２２】
次に、図１３に示すように、iチャンネル目のデータのセンターリング処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０６３）。
【０１２３】
ここでは、ＳＴ２０５３において変位量算出部３０１が算出した変位量に基づいて、センターリング実施部３０２がｉチャンネル目の表面コイルＣｉの位相軸方向処理データに対してセンターリング処理を実施する。
【０１２４】
具体的には、例えば、図１５（ａ）に示すように、ｉチャンネル目の表面コイルＣｉが表面コイルＣ１である場合には、ｋ空間において、図１５（ａ）に示す位相軸方向処理データＤ１１を、図１５（ｂ）に示すように、変位量ＦＤ０３だけ変位させて、センターリング実施部３０２がセンターリング処理を実施する。この場合、位相軸方向処理データＤ１１におけるピークが、ｋ空間における周波数軸方向の中心Ｏから（ＦＤ０３−ＦＤ０１）だけずれるが、ホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１１に適切にフィルターがかかることとなる。
【０１２５】
また、上記と同様に位相軸方向処理データＤ１２および位相軸方向処理データＤ１３に対しても変位量ＦＤ０３に基づいてｋ空間において変位させて、センターリング実施部３０２がセンターリング処理を実施する。この場合にも、ピークが、ｋ空間における周波数軸方向の中心ＯにあるホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１２およびＤ１３に適切にかかることとなる。
【０１２６】
以上のように、本発明の第３の実施形態は、複数の表面コイルにより受信したエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出されたそれぞれの周波数軸方向処理データについて他の周波数軸方向処理データとの相関をとり、相関係数を算出する。そして、相関係数の一番大きいチャンネルの周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、当該周波数差を全ての周波数軸方向処理データを変位させるための変位量と設定する。そして、算出された変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。
【０１２７】
このように、周波数軸方向処理データの最大値がピークとならない場合、すなわち、ノイズ部分の方が実際のピークよりも信号強度が高くなる周波数軸方向処理データは、他の周波数軸方向処理データとの相関が低いため、本実施形態における基準周波数軸方向処理データとして選択されることはない。したがって、本発明にかかる第３の実施形態は、相関係数の一番大きい周波数軸方向処理データにおける信号強度が最大となる点をｋ空間における周波数軸方向の中心Ｏに変位させるための変位量に基づいて、全ての位相軸方向処理データをｋ空間において変位させることにより、生成した画像のアーチファクトを防止することができる。よって、画像品質を向上させることができる。
【０１２８】
＜第４の実施形態＞
以下より本発明にかかる第４の実施形態について図１６により説明する。
第４の実施形態は、動作フローにおけるステップＳＴ２０５４である変位量の平均値の算出および、ステップＳＴ２０６４であるセンターリング処理の実施以外の部分は第１の実施形態と同じである。そのため、重複する箇所については、記載を省略する。
【０１２９】
（動作）
以下より、本実施形態の磁気共鳴イメージング装置１を用いて、被検体４０を撮像する際の動作について説明する。
【０１３０】
図１６は、本発明にかかる第４の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図１７は、本発明にかかる第４の実施形態における周波数軸方向処理データの信号強度の最大値とｋ空間における周波数軸方向の中心Ｏの周波数差を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。図１８は、本発明にかかる第４の実施形態における全ての周波数差の平均値である変位量に基づいて、位相軸方向処理データのセンターリング処理を示した図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。また、変位量ＦＡＶ１は、全ての周波数軸方向処理データの変位量の平均値である。
【０１３１】
まず、図１６に示すように、全てのチャンネルにおけるデータの周波数差の平均値を算出する（ＳＴ２０５４）。
【０１３２】
ここでは、ＳＴ２０２１において算出された、それぞれの周波数軸方向処理データの信号強度が最大値を示す点における周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、全ての周波数軸方向処理データにおける周波数差の平均値を算出する。
【０１３３】
具体的には、例えば、図１７（ａ）に示すように、変位量算出部３０１が、周波数軸方向処理データＤ０１の信号強度が最大となる点における周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差ＦＤ０１を算出する。また、同様に図１７（ｂ）、図１７（ｃ）に示すように、変位量算出部３０１が、周波数軸方向処理データＤ０２、Ｄ０３それぞれの信号強度が最大となる点における周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差ＦＤ０２、周波数差ＦＤ０３をそれぞれ算出する。そして、それぞれの周波数軸方向処理データの周波数差の平均値である（ＦＤ０１＋ＦＤ０２＋ＦＤ０３）／３を変位量算出部３０１が算出する。この平均値を全てのチャンネルの位相軸方向処理データの変位量ＦＡＶ１として設定する。ここで、周波数軸方向処理データと位相軸方向処理データのピーク位置は変わらないため、位相軸方向処理データを変位させる変位量として、変位量ＦＡＶ１を用いることができる。
【０１３４】
次に、図１６に示すように、iチャンネル目のデータのセンターリング処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０６４）。
【０１３５】
ここでは、ＳＴ２０５４において変位量算出部３０１が算出した変位量ＦＡＶ１に基づいて、ステップＳＴ２０３１において画像再構成部３３が算出したiチャンネル目の表面コイルＣｉの位相軸方向処理データに対するセンターリング処理をセンターリング実施部３０２が実施する。
【０１３６】
具体的には、例えば、iチャンネル目の表面コイルＣｉが表面コイルＣ１である場合には、ｋ空間において、図１８（ａ）に示す位相軸方向処理データＤ１１を、図１８（ｂ）に示すように変位量ＦＡＶ１だけ変位させて、センターリング実施部３０２がセンターリング処理を実施する。センターリング処理を実施した後の位相軸方向処理データＤ１１における信号強度が最大となる点は、ｋ空間における周波数軸方向の中心Ｏから（ＦＡＶ１−ＦＤ０１）だけずれているが、位相軸方向処理データＤ１１のピークがホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１１に適切にフィルターがかかることとなる。また、上記と同様に、iチャンネル目の表面コイルＣｉが表面コイルＣ２またはＣ３である場合にも、位相軸方向処理データＤ１２またはＤ１３を変位量ＦＡＶ１だけ変位させて、センターリング処理を実施させる。そして、この場合にも、位相軸方向処理データＤ１２、Ｄ１３における信号強度が最大となる点はそれぞれ、ｋ空間における周波数軸方向の中心Ｏから（ＦＡＶ１−ＦＤ０２）、（ＦＡＶ１−ＦＤ０３）だけずれている。しかし、位相軸方向処理データＤ１２、Ｄ０３のピークがホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１２、Ｄ０３に適切にフィルターがかかることとなる。
【０１３７】
上記のように、位相軸方向処理データＤ１１、Ｄ１２およびＤ１３をｋ空間において、変位量ＦＡＶ１だけ変位させると、それぞれの周波数軸方向処理データにおける信号強度が最大となる点は、ｋ空間における周波数軸方向の中心Ｏには変位しない。しかし、ｋ空間における周波数軸方向の中心ＯにあるホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦが、変位された位相軸方向処理データＤ１１、Ｄ１２およびＤ１３に適切にかかるため、アーチファクトは殆ど生じない。
【０１３８】
以上のように、本発明の第４の実施形態は、複数の表面コイルにより受信した全てのエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出された、それぞれの周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。そして、全ての周波数軸方向処理データにおける周波数差の平均を算出する。そして、算出した周波数差の平均値を変位量として設定し、当該変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。
【０１３９】
このように、本発明にかかる第４の実施形態は、全ての周波数軸方向処理データそれぞれについて、信号強度が最大となる点をｋ空間における周波数軸方向の中心Ｏに変位させるための変位量を算出し、その全ての変位量の平均値に基づいて、全ての位相軸方向処理データをｋ空間において変位させることにより、ノイズの影響のある周波数軸方向処理データが含まれていたとしても、殆ど影響が及ぼされることはない。したがって、生成した画像のアーチファクトを防止することができ、画像品質を向上させることができる。
【０１４０】
＜第５の実施形態＞
以下より本発明にかかる第５の実施形態について図１９により説明する。
第５の実施形態は、動作フローにおけるステップＳＴ２０５５である積分値と変位量の積の平均値の算出および、ステップＳＴ２０６５であるセンターリングの処理の実施以外の部分は第２の実施形態と同じである。そのため、重複する箇所については、記載を省略する。
【０１４１】
（動作）
以下より、本実施形態の磁気共鳴イメージング装置１を用いて、被検体４０を撮像する際の動作について説明する。
【０１４２】
図１９は、本発明にかかる第５の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図２０は、本発明にかかる第５の実施形態における周波数軸方向処理データにおけるプロファイルの積分値と周波数差を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。図２１は、本発明にかかる第５の実施形態における周波数軸方向処理データにおけるプロファイルの積分値と周波数差との積の平均値に基づいて、位相軸方向処理データのセンターリング処理を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。変位量ＦＡＶ２は、全ての周波数軸方向処理データのプロファイルの積分値と変位量の積の平均値である。
【０１４３】
まず、図１９に示すように、全てのチャンネルにおける積分値と周波数差の積の平均値を算出する（ＳＴ２０５５）。
【０１４４】
ここでは、変位量算出部３０１が各チャンネルにおける周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。そして、当該周波数差とステップＳＴ２０２２で算出した各周波数軸方向処理データプロファイルの積分値との積を変位量算出部３０１が算出し、それらの平均値を算出する。
【０１４５】
具体的には、例えば、図２０（ａ）に示すように、変位量算出部３０１が、周波数軸方向処理データＤ０１における信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差ＦＤ０１を算出する。そして、周波数軸方向処理データＤ０１におけるプロファイルの積分値と周波数差の積である（ＳＤ０１×ＦＤ０１）を変位量算出部３０１が算出する。
同様に、周波数軸方向処理データＤ０２におけるプロファイルの積分値と周波数差の積（ＳＤ０２×ＦＤ０２）および周波数軸方向処理データＤ０３におけるプロファイルの積分値と周波数差の積（ＳＤ０３×ＦＤ０３）を変位量算出部３０１が算出する。
【０１４６】
次に、全ての周波数軸方向処理データにおけるプロファイルの積分値と周波数差の積の平均値を算出する。そして、当該積分値と周波数差の積の平均値を全ての位相軸方向処理データを変位させるための変位量として設定する。
例えば、周波数軸方向処理データＤ０１、Ｄ０２およびＤ０３におけるプロファイルの積分値と周波数差の積の平均値を算出する。そして、当該平均値を全てのチャンネルにおける位相軸方向処理データの変位量ＦＡＶ２として設定する。ここで、周波数軸方向処理データと位相軸方向処理データのピーク位置は変わらないため、位相軸方向処理データを変位させる変位量として、変位量ＦＡＶ２を用いることができる。また、最大値と周波数差の積の平均値である変位量ＦＡＶ２は以下の数式（２）により算出される。
【０１４７】
【数２】
2008307228.tif 000004
【０１４８】
次に、図１９に示すように、iチャンネル目のデータのセンターリング処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０６５）。
【０１４９】
ここでは、ＳＴ２０５５において変位量算出部３０１が算出した変位量ＦＡＶ２に基づいて、センターリング実施部３０２がiチャンネル目の表面コイルＣｉの位相軸方向処理データに対してセンターリング処理を実施する。
【０１５０】
具体的には、例えば、iチャンネル目の表面コイルＣｉが表面コイルＣ１である場合には、ｋ空間において、図２１（ａ）に示す位相軸方向処理データＤ１１を、図２１（ｂ）に示すように、変位量ＦＡＶ２だけ変位させて、センターリング実施部３０２がセンターリング処理を実施する。この場合、位相軸方向処理データＤ１１におけるピークは、ｋ空間における周波数軸方向の中心Ｏから（ＦＡＶ２−ＦＤ０１）だけずれているが、ホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１１に適切にフィルターがかかることとなる。また、iチャンネル目の表面コイルＣｉが表面コイルＣ２またはＣ３である場合にも、位相軸方向処理データＤ１１と同様に、位相軸方向処理データＤ１２またはＤ１３をそれぞれ変位量ＦＡＶ２だけ変位させて、センターリング処理を実施する。そして、この場合にも、位相軸方向処理データＤ１２またはＤ１３におけるそれぞれのピークは、ｋ空間における周波数軸方向の中心Ｏから（ＦＡＶ２−ＦＤ０２）、（ＦＡＶ２−ＦＤ０３）だけずれているが、ホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１２またはＤ１３に適切にフィルターがかかることとなる。
【０１５１】
以上のように、本発明の第５の実施形態は、複数の表面コイルにより受信した全てのエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出されたそれぞれの周波数軸方向処理データにおけるプロファイルの積分値を算出する。そして、それぞれの周波数軸方向処理データにおける信号強度が最大となる周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。そして、それぞれの周波数軸方向処理データにおけるプロファイルの積分値と周波数差の積を算出し、全てのエコー信号におけるプロファイルの積分値と周波数差の積の平均値を算出する。そして、算出した平均値を変位量と設定し、当該変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。
【０１５２】
このように、本発明にかかる第５の実施形態は、全ての周波数軸方向処理データにおけるプロファイルの積分値と周波数差との積の平均値に基づいて、全ての位相軸方向処理データをｋ空間において変位させることにより、ノイズの影響のある位相軸方向処理データが含まれていたとしても、殆ど影響が及ぼされることはない。したがって、生成した画像のアーチファクトを防止することができ、画像品質を向上させることができる。
【０１５３】
＜第６の実施形態＞
以下より本発明にかかる第６の実施形態について図２２により説明する。
第６の実施形態は、動作フローにおけるステップＳＴ２０５６である最大値と変位量の積の平均値の算出および、ステップＳＴ２０６６であるセンターリングの処理の実施以外の部分は第１の実施形態と同じである。そのため、重複する箇所については、記載を省略する。
【０１５４】
（動作）
以下より、本実施形態の磁気共鳴イメージング装置１を用いて、被検体４０を撮像する際の動作について説明する。
【０１５５】
図２２は、本発明にかかる第６の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図２３は、本発明にかかる第６の実施形態における周波数軸方向処理データにおける信号強度の最大値と周波数差を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。
図２４は、本発明にかかる第６の実施形態における周波数軸方向処理データにおける信号強度の最大値と周波数差との積の平均値に基づいて、位相軸方向処理データのセンターリング処理を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。また、変位量ＦＡＶ３は、全ての周波数軸方向処理データの最大値と変位量の積の平均値である。
【０１５６】
まず、図２２に示すように、全てのチャンネルにおける信号強度の最大値と周波数差の積の平均値を算出する（ＳＴ２０５６）。
【０１５７】
ここでは、変位量算出部３０１が各チャンネルにおける周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。そして、変位量算出部３０１が当該周波数差とステップＳＴ２０２１で算出した各周波数軸方向処理データの信号強度の最大値との積を算出し、周波数差と信号強度の最大値との積の平均値を算出する。
【０１５８】
具体的には、例えば、図２３（ａ）に示すように、変位量算出部３０１が、周波数軸方向処理データＤ０１における信号強度が最大（ＰＤ０１）となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差ＦＤ０１を算出する。そして、周波数軸方向処理データＤ０１における信号強度の最大値と周波数差の積である（ＰＤ０１×ＦＤ０１）を変位量算出部３０１が算出する。
同様に、周波数軸方向処理データＤ０２における信号強度の最大値と周波数差の積（ＰＤ０２×ＦＤ０２）および周波数軸方向処理データＤ０３における最大値と周波数差の積（ＰＤ０３×ＦＤ０３）を変位量算出部３０１が算出する。
【０１５９】
次に、全ての周波数軸方向処理データにおける最大値と周波数差の積の平均値を算出する。そして、当該最大値と周波数差の積の平均値を全ての位相軸方向処理データを変位させるための変位量として設定する。
例えば、周波数軸方向処理データＤ０１、Ｄ０２およびＤ０３における最大値と周波数差の積の平均値を算出する。そして、当該平均値を全てのチャンネルにおける位相軸方向処理データの変位量ＦＡＶ３として設定する。ここで、周波数軸方向処理データと位相軸方向処理データのピーク位置は変わらないため、位相軸方向処理データを変位させる変位量として、変位量ＦＡＶ３を用いることができる。また、最大値と周波数差の積の平均値である変位量ＦＡＶ３は以下の数式（３）により算出される。
【０１６０】
【数３】
2008307228.tif 000005
【０１６１】
次に、図２２に示すように、iチャンネル目のデータのセンターリング処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０６６）。
【０１６２】
ここでは、ＳＴ２０５６において変位量算出部３０１が算出した変位量ＦＡＶ３に基づいて、センターリング実施部３０２がiチャンネル目の表面コイルＣｉの位相軸方向処理データに対してセンターリング処理を実施する。
【０１６３】
具体的には、例えば、iチャンネル目の表面コイルＣｉが表面コイルＣ１である場合には、ｋ空間において、図２４（ａ）に示す位相軸方向処理データＤ１１を、図２４（ｂ）に示すように、変位量ＦＡＶ３だけ変位させて、センターリング実施部３０２がセンターリング処理を実施する。この場合、位相軸方向処理データＤ１１におけるピークは、ｋ空間における周波数軸方向の中心Ｏから（ＦＡＶ３−ＦＤ０１）だけずれているが、ホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１１に適切にフィルターがかかることとなる。また、iチャンネル目の表面コイルＣｉが表面コイルＣ２またはＣ３である場合にも、位相軸方向処理データＤ１１と同様に、位相軸方向処理データＤ１２またはＤ１３をそれぞれ変位量ＦＡＶ２だけ変位させて、センターリング処理を実施する。そして、この場合にも、位相軸方向処理データＤ１２またはＤ１３におけるそれぞれのピークは、ｋ空間における周波数軸方向の中心Ｏから（ＦＡＶ３−ＦＤ０２）、（ＦＡＶ３−ＦＤ０３）だけずれているが、ホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１２またはＤ１３に適切にフィルターがかかることとなる。
【０１６４】
以上のように、本発明の第６の実施形態は、複数の表面コイルにより受信した全てのエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出されたそれぞれの周波数軸方向処理データにおける信号強度の最大値を算出する。そして、それぞれの周波数軸方向処理データにおける信号強度が最大となる周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出する。そして、それぞれの周波数軸方向処理データにおける信号強度の最大値と周波数差の積を算出し、全てのエコー信号における信号強度の最大値と周波数差の積の平均値を算出する。そして、算出した平均値を変位量と設定し、当該変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。
【０１６５】
このように、本発明にかかる第６の実施形態は、全ての周波数軸方向処理データにおける信号強度の最大値と周波数差との積の平均値に基づいて、全ての位相軸方向処理データをｋ空間において変位させることにより、ノイズの影響のある位相軸方向処理データが含まれていたとしても、殆ど影響が及ぼされることはない。したがって、生成した画像のアーチファクトを防止することができ、画像品質を向上させることができる。
【０１６６】
＜第７の実施形態＞
以下より本発明にかかる第７の実施形態について図２５により説明する。
第７の実施形態は、動作フローにおけるステップＳＴ２０５７である全てのチャンネルの変位量の中央値を算出の算出および、ステップＳＴ２０６７であるセンターリングの処理の実施以外の部分は第１の実施形態と同じである。そのため、重複する箇所については、記載を省略する。
【０１６７】
（動作）
以下より、本実施形態の磁気共鳴イメージング装置１を用いて、被検体４０を撮像する際の動作について説明する。
【０１６８】
図２５は、本発明にかかる第７の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図２６は、本発明にかかる第７の実施形態におけるそれぞれの周波数軸方向処理データを周波数軸方向の中心Ｏに変位させるための変位量の中央値に基づいて、位相軸方向処理データのセンターリング処理を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。
【０１６９】
まず、図２５に示すように、全てのチャンネルにおけるデータのそれぞれの変位量の中央値を算出（ＳＴ２０５７）。
【０１７０】
ここでは、変位量算出部３０１が各チャンネルにおける周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、全てのチャンネルにおける周波数軸方向処理データにおける周波数差を昇順または降順に並べたときに中央となる周波数差である中央値を算出する。
【０１７１】
具体的には、例えば、図２０に示すように、変位量算出部３０１が、周波数軸方向処理データＤ０１における信号強度が最大（ＰＤ０１）となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差ＦＤ０１を算出する。
同様に、周波数軸方向処理データＤ０２における周波数差ＦＤ０２および周波数軸方向処理データＤ０３における周波数差ＦＤ０３を変位量算出部３０１が算出する。
【０１７２】
次に、変位量算出部３０１が全ての周波数軸方向処理データにおけるそれぞれの周波数差の中央値を算出し、中央値を全ての位相軸方向処理データを変位させるための変位量として設定する。
例えば、図２０に示すように、周波数差の大きさがＦＤ０１＜ＦＤ０２＜ＦＤ０３となる場合、周波数差の中央値は周波数差ＦＤ０２となる。そして、この中央値である周波数差ＦＤ０２を全ての位相軸方向処理データを変位させるための変位量として設定する。
【０１７３】
次に、図２５に示すように、iチャンネル目のデータのセンターリング処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０６７）。
【０１７４】
ここでは、ＳＴ２０５７において変位量算出部３０１が算出した変位量に基づいて、センターリング実施部３０２がiチャンネル目の表面コイルＣｉの位相軸方向処理データに対してセンターリング処理を実施する。
【０１７５】
具体的には、例えば、図２６（ａ）に示すように、iチャンネル目の表面コイルＣｉが表面コイルＣ１である場合には、ｋ空間において、図２６（ａ）に示す位相軸方向処理データＤ１１を、図２６（ｂ）に示すように、変位量ＦＤ０２だけ変位させて、センターリング実施部３０２がセンターリング処理を実施する。この場合、位相軸方向処理データＤ１１におけるピークが、ｋ空間における周波数軸方向の中心Ｏから（ＦＤ０２−ＦＤ０１）だけずれるが、ホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１１に適切にフィルターがかかることとなる。
【０１７６】
また、上記と同様に位相軸方向処理データＤ１２および位相軸方向処理データＤ１３に対しても変位量ＦＤ０２に基づいてｋ空間において変位させて、センターリング処理を実施する。この場合にも、ピークが、ｋ空間における周波数軸方向の中心ＯにあるホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１２およびＤ１３に適切にかかることとなる。
【０１７７】
以上のように、本発明の第７の実施形態は、複数の表面コイルにより受信した全てのエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出されたそれぞれの周波数軸方向処理データにおける信号強度の最大値を算出する。そして、それぞれの周波数軸方向処理データにおける信号強度が最大となる周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、全てのチャンネルにおける周波数軸方向処理データにおける周波数差の中央値を算出する。そして、当該中央値を全ての周波数軸方向処理データを変位させるための変位量と設定する。そして、算出された変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。
【０１７８】
このように、全てのチャンネルにおける周波数軸方向処理データにおける周波数差の中央値を、周波数軸方向処理データを変位させるための変位量とすることにより、全ての周波数軸方向処理データのピーク値が、ｋ空間における周波数軸方向の中心Ｏの近傍に変位されるため、全てのチャンネルにおける周波数軸方向処理データにおける周波数差の中央値に基づいて、全ての位相軸方向処理データをｋ空間において変位させることにより、生成した画像のアーチファクトを防止することができる。よって、画像品質を向上させることができる。
【０１７９】
＜第８の実施形態＞
以下より本発明にかかる第８の実施形態について図２７により説明する。
第８の実施形態は、動作フローにおけるステップＳＴ２０５８である最大値と変位量の積の平均値の算出および、ステップＳＴ２０６８であるセンターリングの処理の実施以外の部分は第１の実施形態と同じである。そのため、重複する箇所については、記載を省略する。
【０１８０】
（動作）
以下より、本実施形態の磁気共鳴イメージング装置１を用いて、被検体４０を撮像する際の動作について説明する。
【０１８１】
図２７は、本発明にかかる第８の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。図２８は、本発明にかかる第８の実施形態における周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差の標準偏差に基づいて、選択された周波数差の平均値により位相軸方向処理データのセンターリング処理を示す図である。縦軸は信号強度Ｉであり、横軸は周波数ｆである。また、変位量ＦＡＶ４は、標準偏差に基づいて選択された周波数差の平均値である。
【０１８２】
まず、図２７に示すように、標準偏差に基づいて選択されたデータの変位量の平均値を算出（ＳＴ２０５８）。
【０１８３】
ここでは、変位量算出部３０１が各チャンネルにおける周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、全てのチャンネルにおける周波数軸方向処理データにおける周波数差の標準偏差を算出する。そして、標準偏差に基づいて周波数差を選択し、選択された周波数差の平均値を算出する。
【０１８４】
具体的には、例えば、図２０に示すように、変位量算出部３０１が、周波数軸方向処理データＤ０１における信号強度が最大（ＰＤ０１）となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差ＦＤ０１を算出する。
同様に、周波数軸方向処理データＤ０２における周波数差ＦＤ０２および周波数軸方向処理データＤ０３における周波数差ＦＤ０３を変位量算出部３０１が算出する。
【０１８５】
次に、変位量算出部３０１が全ての周波数軸方向処理データにおける周波数差の標準偏差を算出し、当該標準偏差に基づいて、全ての周波数差の中から周波数差を選択する。
例えば、周波数軸方向処理データＤ０１、Ｄ０２およびＤ０３における周波数差の相加平均値と標準偏差σを算出する。そして、｜周波数差−相加平均値｜＜２σとなる周波数差を選択する。そして、選択された周波数差の平均値を全てのチャンネルにおける位相軸方向処理データの変位量ＦＡＶ４として設定する。例えば、周波数軸方向処理データＤ０１、Ｄ０２およびＤ０３において、｜周波数差−相加平均値｜＜２σとなる周波数軸方向処理データがＤ０１、Ｄ０２の場合、変位量ＦＡＶ４＝（ＦＤ０１＋ＦＤ０２）／２となる。
ここで、周波数軸方向処理データと位相軸方向処理データのピーク位置は変わらないため、位相軸方向処理データを変位させる変位量として、変位量ＦＡＶ３を用いることができる。また、周波数差の標準偏差は以下の数式（４）により算出される。
【０１８６】
【数４】
2008307228.tif 000006
【０１８７】
次に、図２７に示すように、iチャンネル目のデータのセンターリング処理を実施する（i＝０，・・・，Ｎ）（ＳＴ２０６８）。
【０１８８】
ここでは、ＳＴ２０５７において変位量算出部３０１が算出した変位量に基づいて、iチャンネル目の表面コイルＣｉの位相軸方向処理データに対してセンターリング処理を実施する。
【０１８９】
具体的には、例えば、図２８（ａ）に示すように、iチャンネル目の表面コイルＣｉが表面コイルＣ１である場合には、ｋ空間において、図２８（ａ）に示す位相軸方向処理データＤ１１を、図２８（ｂ）に示すように、変位量ＦＡＶ３だけ変位させて、センターリング処理を実施する。この場合、位相軸方向処理データＤ１１におけるピークが、ｋ空間における周波数軸方向の中心Ｏから（ＦＡＶ４−ＦＤ０１）だけずれるが、ホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１１に適切にフィルターがかかることとなる。
【０１９０】
また、上記と同様に位相軸方向処理データＤ１２および位相軸方向処理データＤ１３に対しても変位量ＦＡＶ４に基づいてｋ空間において変位させて、センターリング処理を実施する。この場合にも、ピークが、ｋ空間における周波数軸方向の中心ＯにあるホモダインハイパスフィルターＨＦおよびホモダインローパスフィルターＬＦを通るため、位相軸方向処理データＤ１２およびＤ１３に適切にかかることとなる。
【０１９１】
以上のように、本発明の第８の実施形態は、複数の表面コイルにより受信した全てのエコー信号に対して、周波数軸方向にＦＦＴ処理を実施することにより算出されたそれぞれの周波数軸方向処理データにおける信号強度の最大値を算出する。そして、それぞれの周波数軸方向処理データにおける信号強度が最大となる周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差を算出し、全てのチャンネルにおける周波数軸方向処理データにおける周波数差の標準偏差を算出する。そして、当該標準偏差に基づいて、全てのチャンネルにおける周波数軸方向処理データにおける周波数差から周波数差を選択し、選択した周波数差の平均値を全ての周波数軸方向処理データを変位させるための変位量と設定する。そして、算出された変位量に基づいて、位相軸方向処理データを変位させ、センターリング処理を実施する。そして、センターリング処理を実施したデータに対してハーフエコー処理を実施して画像を生成する。
【０１９２】
このように、全てのチャンネルにおける周波数軸方向処理データにおける周波数差の標準偏差に基づいて、周波数差を選択することにより、極端に値の異なる周波数差を、平均値を算出するための周波数差から除外することができる。そのため、標準偏差に基づいて選択した周波数差における平均値に基づいて、全ての位相軸方向処理データをｋ空間において変位させることにより、生成した画像のアーチファクトを防止することができる。よって、画像品質を向上させることができる。
【０１９３】
なお、上記の本実施形態におけるエコー信号は、本発明の磁気共鳴信号に相当する。また、本実施形態の基準周波数軸方向処理データは、本発明の基準磁気共鳴信号に相当する。また、本実施形態のスキャン部２は、本発明のスキャン部に相当する。また、本実施形態の変位量算出部３０１は、本発明の変位量算出部に相当する。また、本実施形態のセンターリング実施部３０２は、本発明のセンターリング実施部に相当する。また、本実施形態の画像再構成部３３は、本発明の画像再構成部に相当する。
【０１９４】
なお、本発明の実施に際しては、上記した実施形態に限定されるものではなく、種々の変形形態を採用することができる。
【０１９５】
本発明の実施形態において、エコー信号にＦＦＴ処理を実施したデータに対して、センターリング処理を実施しているが、これに限定されず、エコー信号に対してセンターリング処理を実施したデータに対して、ＦＦＴ処理を実施しても良い。また、本発明の第５の実施形態において、積分値と周波数差との積の平均値に基づいて、センターリング処理を実施しているが、これに限定されず、エコー信号のプロファイル特性と周波数差との積の平均値に基づいて、センターリング処理を実施しても良い。例えば、エコー信号のプロファイル特性をエコー信号プロファイル全体の平均値としても良い。また、その他のエコー信号のプロファイル特性であっても良い。
【図面の簡単な説明】
【０１９６】
【図１】図１は、本発明にかかる一実施形態のＲＦコイル部により構成される磁気共鳴イメージング装置の構成を示す構成図である。
【図２】図２に、フェーズドアレイコイルの一例を示す。
【図３】図３は、フェーズドアレイコイルにより受信したエコー信号にフーリエ変換処理を実施したデータを示す図である。
【図４】図４は、本発明にかかる一実施形態における画像再構成部３３の構成を示すブロック図である。
【図５】図５は、本発明にかかる第１の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。
【図６】図６は、本発明にかかる第１の実施形態における受信した複数の位相軸方向処理データのセンターリング処理の実施の動作を示すフロー図である。
【図７】図７は、本発明にかかる第１の実施形態における周波数軸方向処理データを示す図である。
【図８】図８は、本発明に係る第１の実施形態におけるノイズの影響のある周波数軸方向処理データを示す図である。
【図９】図９は、本発明にかかる第１の実施形態において、位相軸方向処理データのセンターリング処理を示した図である。
【図１０】図１０は、本発明にかかる第１の実施形態において、位相軸方向処理データのセンターリング処理を示した図である。
【図１１】図１１は、本発明にかかる第２の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。
【図１２】図１２は、本発明にかかる第２の実施形態における周波数軸方向処理データのプロファイルにおける積分値を示す図である。
【図１３】図１３は、本発明にかかる第３の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。
【図１４】図１４は、本発明にかかる第３の実施形態における周波数軸方向処理データそれぞれの相関を示す相関図である。
【図１５】図１５は、本発明にかかる第３の実施形態における周波数軸方向処理データの相関係数が最大である周波数軸方向処理データの周波数差に基づいて、位相軸方向処理データのセンターリング処理を示した図である。
【図１６】図１６は、本発明にかかる第４の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。
【図１７】図１７は、本発明にかかる第４の実施形態における周波数軸方向処理データの信号強度の最大値と周波数差を示す図である。
【図１８】図１８は、本発明にかかる第４の実施形態における全ての周波数差の平均値である変位量に基づいて、位相軸方向処理データのセンターリング処理を示した図である。
【図１９】図１９は、本発明にかかる第５の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。
【図２０】図２０は、本発明にかかる第５の実施形態における周波数軸方向処理データにおけるプロファイルの積分値と周波数差を示す図である。
【図２１】図２１は、本発明にかかる第５の実施形態における周波数軸方向処理データにおけるプロファイルの積分値と周波数差との積の平均値に基づいて、位相軸方向処理データのセンターリング処理を示す図である。
【図２２】図２２は、本発明にかかる第６の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。
【図２３】図２３は、本発明にかかる第６の実施形態における周波数軸方向処理データにおける信号強度の最大値と周波数差を示す図である。
【図２４】図２４は、本発明にかかる第６の実施形態における周波数軸方向処理データにおける信号強度の最大値と周波数差との積の平均値に基づいて、位相軸方向処理データのセンターリング処理を示す図である。
【図２５】図２５は、本発明にかかる第７の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。
【図２６】図２６は、本発明にかかる第７の実施形態におけるそれぞれの周波数軸方向処理データを周波数軸方向の中心Ｏに変位させるための変位量の中央値に基づいて、位相軸方向処理データのセンターリング処理を示す図である。
【図２７】図２７は、本発明にかかる第８の実施形態において、被検体４０を撮像する際の動作を示すフロー図である。
【図２８】図２８は、本発明にかかる第８の実施形態における周波数軸方向処理データにおける信号強度が最大となる点の周波数とｋ空間における周波数軸方向の中心Ｏにおける周波数との周波数差の標準偏差に基づいて、選択された周波数差の平均値により位相軸方向処理データのセンターリング処理を示す図である。
【図２９】図２９は、従来技術を説明するための図である。
【図３０】図３０は、従来技術を説明するための図である。
【図３１】図３１は、従来技術を説明するための図である。
【符号の説明】
【０１９７】
１：磁気共鳴イメージング装置、２:スキャン部、３：操作コンソール部、１１：静磁場空間、１２：静磁場マグネット部、１３：勾配コイル部、１４：ＲＦコイル部、１５：クレードル、２２：ＲＦ駆動部、２３：勾配駆動部、２４：データ収集部、３０：制御部、３１：記憶部、３２：操作部、３３：画像再構成部、３４：表示部、４０：被検体、３０１：変位量算出部、３０２：センターリング実施部、Ｃ１：表面コイル、Ｃ２：表面コイル、Ｃ３：表面コイル、Ｄ０１：表面コイルＣ１の周波数軸方向処理データ、Ｄ０２：表面コイルＣ２の周波数軸方向処理データ、Ｄ０３：表面コイルＣ３の周波数軸方向処理データ、ＦＤ０１：周波数軸方向処理データＤ０１の変位量、ＦＤ０２：周波数軸方向処理データＤ０２の変位量、ＦＤ０３：周波数軸方向処理データＤ０３の変位量、ＦＡＶ１：全ての周波数軸方向処理データにおける変位量の平均値、ＦＡＶ２：全ての周波数軸方向処理データにおける最大値と変位量の積の平均値
【図１】
2008307228.tif 000007
【図２】
2008307228.tif 000008
【図３】
2008307228.tif 000009
【図４】
2008307228.tif 000010
【図５】
2008307228.tif 000011
【図６】
2008307228.tif 000012
【図７】
2008307228.tif 000013
【図８】
2008307228.tif 000014
【図９】
2008307228.tif 000015
【図１０】
2008307228.tif 000016
【図１１】
2008307228.tif 000017
【図１２】
2008307228.tif 000018
【図１３】
2008307228.tif 000019
【図１４】
2008307228.tif 000020
【図１５】
2008307228.tif 000021
【図１６】
2008307228.tif 000022
【図１７】
2008307228.tif 000023
【図１８】
2008307228.tif 000024
【図１９】
2008307228.tif 000025
【図２０】
2008307228.tif 000026
【図２１】
2008307228.tif 000027
【図２２】
2008307228.tif 000028
【図２３】
2008307228.tif 000029
【図２４】
2008307228.tif 000030
【図２５】
2008307228.tif 000031
【図２６】
2008307228.tif 000032
【図２７】
2008307228.tif 000033
【図２８】
2008307228.tif 000034
【図２９】
2008307228.tif 000035
【図３０】
2008307228.tif 000036
【図３１】
2008307228.tif 000037