(57)【特許請求の範囲】
【請求項１】
被検体の体動によって変位する第１の部位のｋ空間のデータを、複数のスキャンに分けて収集する磁気共鳴イメージング装置であって、
前記複数のスキャンのうちの第１のスキャンでは、第１の補正用ｋ空間の所定のデータと、画像用ｋ空間の第１の領域のデータとが収集され、
前記複数のスキャンのうちの第２のスキャンでは、第２の補正用ｋ空間の所定のデータと、前記画像用ｋ空間の第２の領域のデータとが収集され、
前記第１の補正用ｋ空間の所定のデータと、前記第２の補正用ｋ空間の所定のデータとに基づいて、前記画像用ｋ空間の第１の領域および第２の領域のうちの少なくともいずれか一方の領域のデータを補正する補正手段と、
前記補正手段により補正された後の前記画像用ｋ空間のデータを用いて、前記第１の部位の画像を再構成する再構成手段と、
を有し、
前記補正手段は、前記少なくともいずれか一方の領域のデータの強度を補正する、磁気共鳴イメージング装置。
【請求項２】
前記補正手段は、
前記第１の補正用ｋ空間の所定のデータの強度と、前記第２の補正用ｋ空間の所定のデータの強度とに基づいて、前記少なくともいずれか一方の領域のデータの強度を補正する、請求項１に記載の磁気共鳴イメージング装置。
【請求項３】
前記補正手段は、
前記少なくともいずれか一方の領域のデータの強度を補正した後、強度が補正されたデータの位相を補正する、請求項２に記載の磁気共鳴イメージング装置。
【請求項４】
前記補正手段は、
前記第１の補正用ｋ空間の所定のデータの位相と前記第２の補正用ｋ空間の所定のデータの位相との差に基づいて、強度が補正されたデータの位相を補正する、請求項３に記載の磁気共鳴イメージング装置。
【請求項５】
被検体の体動によって変位する第１の部位のｋ空間のデータを、複数のスキャンに分けて収集する磁気共鳴イメージング装置であって、
前記複数のスキャンのうちの第１のスキャンでは、第１の補正用ｋ空間の所定のデータと、画像用ｋ空間の第１の領域のデータとが収集され、
前記複数のスキャンのうちの第２のスキャンでは、第２の補正用ｋ空間の所定のデータと、前記画像用ｋ空間の第２の領域のデータとが収集され、
前記第１の補正用ｋ空間の所定のデータと、前記第２の補正用ｋ空間の所定のデータとに基づいて、前記画像用ｋ空間の第１の領域および第２の領域のうちの少なくともいずれか一方の領域のデータを補正する補正手段と、
前記補正手段により補正された後の前記画像用ｋ空間のデータを用いて、前記第１の部位の画像を再構成する再構成手段と、
を有し、
前記補正手段は、前記少なくともいずれか一方の領域のデータの位相を補正し、位相が補正されたデータの強度を補正する、磁気共鳴イメージング装置。
【請求項６】
前記補正手段は、
前記第１の補正用ｋ空間の所定のデータの位相と前記第２の補正用ｋ空間の所定のデータの位相との差に基づいて、前記少なくともいずれか一方の領域のデータの位相を補正する、請求項５に記載の磁気共鳴イメージング装置。
【請求項７】
前記補正手段は、
前記第１の補正用ｋ空間の所定のデータの強度と、前記第２の補正用ｋ空間の所定のデータの強度とに基づいて、位相が補正されたデータの強度を補正する、請求項５又は６に記載の磁気共鳴イメージング装置。
【請求項８】
前記第１の補正用ｋ空間の所定のデータは、前記第１の補正用ｋ空間の中心のデータであり、前記第２の補正用ｋ空間の所定のデータは、前記第２の補正用ｋ空間の中心のデータである、請求項１〜７のうちのいずれか一項に記載の磁気共鳴イメージング装置。
【請求項９】
前記第１のスキャンでは、前記第１の補正用ｋ空間の中心のデータを収集するための第１のシーケンスが実行され、
前記第２のスキャンでは、前記第２の補正用ｋ空間の中心のデータを収集するための第２のシーケンスが実行され、
前記第１のシーケンスおよび前記第２のシーケンスは、位相エンコード勾配磁場および周波数エンコード勾配磁場を印加しないシーケンスである、請求項８に記載の磁気共鳴イメージング装置。
【請求項１０】
前記第１の補正用ｋ空間、前記第２の補正用ｋ空間、および前記画像用ｋ空間は、２次元ｋ空間又は３次元ｋ空間である、請求項１〜９のうちのいずれか一項に記載の磁気共鳴イメージング装置。
【請求項１１】
被検体の体動によって変位する第１の部位のｋ空間のデータを、３Ｄのパルスシーケンスを用いた複数のスキャンに分けて収集する磁気共鳴イメージング装置であって、
前記複数のスキャンのうちの第１のスキャンでは、第１の補正用ｋ空間の所定の位置を横切りｋｚ方向に延在する第１のデータと、画像用ｋ空間の第１の領域のデータとを収集し、
前記複数のスキャンのうちの第２のスキャンでは、第２の補正用ｋ空間の所定の位置を横切りｋｚ方向に延在する第２のデータと、前記画像用ｋ空間の第２の領域のデータとを収集し、
前記第１のデータ、前記第２のデータ、前記第１の領域のデータ、および前記第２の領域のデータを、ｋｚ方向に逆フーリエ変換する逆フーリエ変換手段と、
前記第１のデータを逆フーリエ変換することにより得られた第１の補正用データと、前記第２のデータを逆フーリエ変換することにより得られた第２の補正用データとに基づいて、前記第１の領域のデータを逆フーリエ変換することにより得られた第１の画像用データと、前記第２の領域のデータを逆フーリエ変換することにより得られた第２の画像用データとのうちの少なくともいずれか一方の画像用データを補正する補正手段と、
前記補正手段により補正された画像用データを用いて、前記第１の部位の画像を再構成する再構成手段と、
を有し、
前記補正手段は、前記少なくともいずれか一方の画像用データの強度を補正する、磁気共鳴イメージング装置。
【請求項１２】
前記補正手段は、
前記第１の補正用データの強度と前記第２の補正用データの強度とに基づいて、前記少なくともいずれか一方の画像用データの強度を補正する、請求項１１に記載の磁気共鳴イメージング装置。
【請求項１３】
前記補正手段は、
前記少なくともいずれか一方の画像用データの強度を補正した後、強度が補正された画像用データの位相を補正する、請求項１２に記載の磁気共鳴イメージング装置。
【請求項１４】
前記補正手段は、
前記第１の補正用データの位相と前記第２の補正用データの位相との差に基づいて、強度が補正された画像用データの位相を補正する、請求項１３に記載の磁気共鳴イメージング装置。
【請求項１５】
被検体の体動によって変位する第１の部位のｋ空間のデータを、３Ｄのパルスシーケンスを用いた複数のスキャンに分けて収集する磁気共鳴イメージング装置であって、
前記複数のスキャンのうちの第１のスキャンでは、第１の補正用ｋ空間の所定の位置を横切りｋｚ方向に延在する第１のデータと、画像用ｋ空間の第１の領域のデータとを収集し、
前記複数のスキャンのうちの第２のスキャンでは、第２の補正用ｋ空間の所定の位置を横切りｋｚ方向に延在する第２のデータと、前記画像用ｋ空間の第２の領域のデータとを収集し、
前記第１のデータ、前記第２のデータ、前記第１の領域のデータ、および前記第２の領域のデータを、ｋｚ方向に逆フーリエ変換する逆フーリエ変換手段と、
前記第１のデータを逆フーリエ変換することにより得られた第１の補正用データと、前記第２のデータを逆フーリエ変換することにより得られた第２の補正用データとに基づいて、前記第１の領域のデータを逆フーリエ変換することにより得られた第１の画像用データと、前記第２の領域のデータを逆フーリエ変換することにより得られた第２の画像用データとのうちの少なくともいずれか一方の画像用データを補正する補正手段と、
前記補正手段により補正された画像用データを用いて、前記第１の部位の画像を再構成する再構成手段と、
を有し、
前記補正手段は、前記少なくともいずれか一方の画像用データの位相を補正し、位相が補正された画像用データの強度を補正する、磁気共鳴イメージング装置。
【請求項１６】
前記補正手段は、
前記第１の補正用データの位相と前記第２の補正用データの位相との差に基づいて、前記少なくともいずれか一方の画像用データの位相を補正する、請求項１５に記載の磁気共鳴イメージング装置。
【請求項１７】
前記補正手段は、
前記第１の補正用データの強度と前記第２の補正用データの強度とに基づいて、位相が補正された画像用データの強度を補正する、請求項１５又は１６に記載の磁気共鳴イメージング装置。
【請求項１８】
前記第１の補正用ｋ空間の所定の位置は、前記第１の補正用ｋ空間の中心であり、
前記第２の補正用ｋ空間の所定の位置は、前記第２の補正用ｋ空間の中心である、請求項１１〜１７のうちのいずれか一項に記載の磁気共鳴イメージング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体の体動によって変位する所定の部位をスキャンする磁気共鳴イメージング装置に関する。
【背景技術】
【０００２】
横隔膜の移動に合わせてスライス位置を調整するスライストラッキング技術が知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開2004-305454号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
スライストラッキング技術では、撮影対象の臓器の動きに合わせて、スライス位置を変更し、磁気共鳴信号を収集する。しかし、臓器が動くことによって、臓器が感じる受信コイルの感度や静磁場が変化し、その結果、受信コイルの感度の変化がアーチファクトとなって画質を劣化させてしまうことがある。したがって、アーチファクトを低減することが望まれる。
【課題を解決するための手段】
【０００５】
第１の観点の発明は、被検体の体動によって変位する所定の部位のｋ空間のデータを、複数のスキャンに分けて収集する磁気共鳴イメージング装置であって、
前記複数のスキャンのうちの第１のスキャンでは、第１の補正用ｋ空間の中心のデータと、画像用ｋ空間の第１の領域のデータとが収集され、
前記複数のスキャンのうちの第２のスキャンでは、第２の補正用ｋ空間の中心のデータと、前記画像用ｋ空間の第２の領域のデータとが収集され、
前記第１の補正用ｋ空間の中心のデータと、前記第２の補正用ｋ空間の中心のデータとに基づいて、前記画像用ｋ空間の第１の領域および第２の領域のうちの少なくともいずれか一方の領域のデータを補正する補正手段と、
前記補正手段により補正された後の前記画像用ｋ空間のデータを用いて、前記所定の部位の画像を再構成する再構成手段と、を有する。
第２の観点の発明は、被検体の体動によって変位する所定の部位のｋ空間のデータを、３Ｄのパルスシーケンスを用いた複数のスキャンに分けて収集する磁気共鳴イメージング装置であって、
前記複数のスキャンのうちの第１のスキャンでは、第１の補正用ｋ空間の中心を横切りｋｚ方向に延在する第１の直線上のデータと、画像用ｋ空間の第１の領域のデータとを収集し、
前記複数のスキャンのうちの第２のスキャンでは、第２の補正用ｋ空間の中心を横切りｋｚ方向に延在する第２の直線上のデータと、前記画像用ｋ空間の第２の領域のデータとを収集し、
前記第１の直線上のデータ、前記第２の直線上のデータ、前記第１の領域のデータ、および前記第２の領域のデータを、ｋｚ方向に逆フーリエ変換する逆フーリエ変換手段と、
前記第１の直線上のデータを逆フーリエ変換することにより得られた第１の補正用データと、前記第２の直線上のデータを逆フーリエ変換することにより得られた第２の補正用データとに基づいて、前記第１の領域のデータを逆フーリエ変換することにより得られた第１の画像用データと、前記第２の領域のデータを逆フーリエ変換することにより得られた第２の画像用データとのうちの少なくともいずれか一方の画像用データを補正する補正手段と、
前記補正手段により補正された画像用データを用いて、前記所定の部位の画像を再構成する再構成手段と、を有する。
【発明の効果】
【０００６】
所定の部位の変位に伴って受信コイルの感度が変化しても、アーチファクトを低減することができる。
【図面の簡単な説明】
【０００７】
【図１】本発明の第１の実施形態の磁気共鳴イメージング装置の概略図である。
【図２】被検体１２の撮影領域を示す概略図である。
【図３】撮影部位を撮影するときに実行されるスキャンの説明図である。
【図４】１回目のスキャンおよび２回目のスキャンで使用されるイメージングシーケンスＩＳの一例を示す図である。
【図５】ＭＲＩ装置１００の処理フローを示す図である。
【図６】本スキャンによって収集されるｋ空間のデータを示す図である。
【図７】スカウト画像の一例を示す図である。
【図８】１回目のスキャンのナビゲータシーケンスＮＡＶによって検出された被検体１２の横隔膜１２ａの位置を概略的に示す図である。
【図９】補正用ｋ空間データ収集シーケンスＣａ１〜Ｃｎ１を実行するときに設定されるスライスを示す図である。
【図１０】補正用ｋ空間データ収集シーケンスＣａ１〜Ｃｎ１の一例を示す図である。
【図１１】２回目のスキャンのナビゲータシーケンスＮＡＶによって検出された被検体１２の横隔膜１２ａの位置を概略的に示す図である。
【図１２】１回目のスキャンにおけるスライスと、２回目のスキャンにおけるスライスとの違いを示す図である。
【図１３】ステップＳＴ１０およびＳＴ１１で実行される補正の説明図である。
【図１４】比較結果の説明図である。
【図１５】第２の実施形態で実行されるスキャンによって収集されるｋ空間のデータを示す図である。
【図１６】第２の実施形態における補正方法の説明図である。
【図１７】第３の実施形態のＭＲＩ装置３００を示す概略図である。
【図１８】第３の実施形態におけるスキャンの説明図である。
【図１９】第３の実施形態における補正方法の説明図である。
【発明を実施するための形態】
【０００８】
以下、図面を参照しながら、発明を実施するための形態を詳細に説明するが、本発明は、以下の形態に限定されることはない。
【０００９】
（１）第１の実施形態
図１は、本発明の第１の実施形態の磁気共鳴イメージング装置の概略図である。
磁気共鳴イメージング装置（以下、ＭＲＩ装置と呼ぶ。ＭＲＩ（Magnetic Resonance Imaging））１００は、磁場発生装置２と、テーブル３と、受信コイル４などを有している。
【００１０】
磁場発生装置２は、被検体１２が収容されるボア２１、超伝導コイル２２、勾配コイル２３、および送信コイル２４などを有している。超伝導コイル２２は静磁場Ｂ0を印加し、勾配コイル２３は勾配磁場を印加し、送信コイル２４はＲＦパルスを送信する。尚、超伝導コイル２２の代わりに、永久磁石を用いてもよい。
【００１１】
テーブル３は、クレードル３１を有している。クレードル３１は、ボア２１に移動できるように構成されている。クレードル３１がボア２１に移動することによって、被検体１２はボア２１に搬送される。
【００１２】
受信コイル４は、被検体１２の腹部から胸部に渡って取り付けられている。受信コイル４は、被検体１２からの磁気共鳴信号を受信する。
【００１３】
ＭＲＩ装置１００は、更に、シーケンサ５、送信器６、勾配磁場電源７、受信器８、中央処理装置９、入力装置１０、および表示装置１１を有している。
【００１４】
シーケンサ５は、中央処理装置９の制御を受けて、後述するスキャン（図３参照）を実行するための情報を送信器６および勾配磁場電源７に送る。具体的には、シーケンサ５は、中央処理装置９の制御を受けて、ＲＦパルスの情報（中心周波数、バンド幅など）を送信器６に送り、勾配磁場の情報（勾配磁場の強度など）を勾配磁場電源７に送る。
【００１５】
送信器６は、シーケンサ５から送られた情報に基づいて、ＲＦコイル２４を駆動する駆動信号を出力する。
【００１６】
勾配磁場電源７は、シーケンサ５から送られた情報に基づいて、勾配コイル２３を駆動する駆動信号を出力する。
【００１７】
受信器８は、受信コイル４で受信された磁気共鳴信号を信号処理し、中央処理装置９に伝送する。
【００１８】
中央処理装置９は、シーケンサ５および表示装置１１に必要な情報を伝送したり、受信器８から受け取った信号に基づいて画像を再構成するなど、ＭＲＩ装置１００の各種の動作を実現するように、ＭＲＩ装置１００の各部の動作を制御する。中央処理装置９は、例えばコンピュータ（computer）によって構成される。中央処理装置９は、補正手段９１および再構成手段９２などを有している。
【００１９】
補正手段９１は、強度補正手段９１ａおよび位相補正手段９１ｂを有している。強度補正手段９１ａは、ｋ空間の高周波領域のデータＤＡ２〜ＤＮ２（図１３参照）の強度を補正する。位相補正手段９１ｂは、ｋ空間の高周波領域の強度補正後のデータＤＡ２′〜ＤＮ２′（図１３参照）の位相を補正する。
【００２０】
再構成手段９２は、補正手段９１により補正されたｋ空間のデータを用いて、画像を再構成する。
【００２１】
中央処理装置９は、補正手段９１および再構成手段９２の一例であり、所定のプログラムを実行することにより、これらの手段として機能する。
【００２２】
入力装置１０は、オペレータ１３の操作に応じて、種々の命令を中央処理装置９に入力する。表示装置１１は種々の情報を表示する。
【００２３】
ＭＲＩ装置１００は、上記のように構成されている。次に、被検体１２を撮影するときに実行されるスキャンについて説明する。
【００２４】
図２は、被検体１２の撮影領域を示す概略図である。
第１の実施形態では、肝臓１２ｂを含む領域が、撮影領域Ｒimとして設定される。
【００２５】
図３は、撮影部位を撮影するときに実行されるスキャンの説明図である。
第１の実施形態では、スカウトスキャンを実行し、その後に、本スキャンを実行する。
【００２６】
スカウトスキャンは、スカウト画像のデータを収集するためのスキャンである。スカウト画像は、オペレータ１３が撮影領域Ｒim（図２参照）にスライスを設定するときに使用するための画像である。
【００２７】
本スキャンでは、ナビゲータ領域Ｒnavから横隔膜１２ａの位置Ｐを検出するためのデータを収集しながら、撮影領域Ｒimがスキャンされる。尚、横隔膜１２ａおよび肝臓１２ｂは、ｚ方向にのみ変位すると仮定する。本スキャンは、被検体１２に息止めをしてもらった状態で実行される２回のスキャンを有している。
【００２８】
１回目のスキャンでは、ナビゲータシーケンスＮＡＶと、補正用ｋ空間データ収集シーケンスＣａ１〜Ｃｎ１と、イメージングシーケンスＩＳとが実行される。ナビゲータシーケンスＮＡＶは、ナビゲータ領域Ｒnavからデータを収集し、横隔膜１２ａの位置Ｐを検出するためのシーケンスである。補正用ｋ空間データ収集シーケンスＣａ１〜Ｃｎ１は、補正用ｋ空間Ｓａ１〜Ｓｎ１の中心のデータＡ１〜Ｎ１を収集するためのシーケンスである（図６参照）。補正用ｋ空間データ収集シーケンスＣａ１〜Ｃｎ１については、後で詳細に説明する。イメージングシーケンスＩＳは、撮影領域Ｒimを画像化するためのシーケンスである。１回目のスキャンでは、イメージングシーケンスＩＳを複数回実行することによって、各スライスにおけるｋ空間ＳＡ（図６参照）の中央領域のデータが収集される。
【００２９】
２回目のスキャンでは、ナビゲータシーケンスＮＡＶと、補正用ｋ空間データ収集シーケンスＣａ２〜Ｃｎ２と、イメージングシーケンスＩＳとが実行される。補正用ｋ空間データ収集シーケンスＣａ２〜Ｃｎ２は、補正用ｋ空間Ｓａ２〜Ｓｎ２の中心のデータＡ２〜Ｎ２を収集するためのシーケンスである（図６参照）。補正用ｋ空間データ収集シーケンスＣａ２〜Ｃｎ２についても、後で詳細に説明する。補正用ｋ空間データ収集シーケンスＣａ２〜Ｃｎ２が実行された後、イメージングシーケンスＩＳを複数回実行することによって、各スライスにおけるｋ空間ＳＡ（図６参照）の高周波領域のデータが収集される。
【００３０】
図４は、１回目のスキャンおよび２回目のスキャンで使用されるイメージングシーケンスＩＳの一例を示す図である。
【００３１】
第１の実施形態では、イメージングシーケンスＩＳは、２Ｄ（Dimension）イメージング用のパルスシーケンスである。したがって、第１の実施形態では、ｋ空間は、２次元のｋ空間である。
【００３２】
次に、ＭＲＩ装置１００の動作について説明する。
図５は、ＭＲＩ装置１００の処理フローを示す図、図６は、本スキャンによって収集されるｋ空間のデータを示す図である。図５のフローの説明にあたっては、必要に応じて、図６などを参照しながら説明する。
【００３３】
ステップＳＴ１では、スカウト画像のデータを収集するためのスカウトスキャン（図３参照）が実行される。第１の実施形態では、肝臓１２ｂを撮影することを考えているので、肝臓１２ｂを含むスカウト画像のデータを収集するためのスカウトスキャンが行われる。
【００３４】
図７は、スカウト画像の一例を示す図である。
第１の実施形態では、スカウト画像は、コロナル画像である。オペレータ１３は、スカウト画像を参考にしながら、肝臓１２ｂの全体をスキャンすることができるように、スライス位置、スライス厚、およびスライス枚数等を決定する。スライス位置等を決定した後、ステップＳＴ２に進む。
【００３５】
ステップＳＴ２では、オペレータ１３が、被検体１２に対して、１回目の息止めの指示をする。この指示によって、被検体１２は息止めをする。被検体１２に対して息止めの指示をした後、ステップＳＴ３に進む。
【００３６】
ステップＳＴ３では、本スキャンの１回目のスキャンにおけるナビゲータシーケンスＮＡＶ（図６参照）が実行される。ナビゲータシーケンスＮＡＶを実行することによって、被検体１２が息止めをしているときの横隔膜１２ａの位置が検出される。
【００３７】
図８は、１回目のスキャンのナビゲータシーケンスＮＡＶによって検出された被検体１２の横隔膜１２ａの位置を概略的に示す図である。
【００３８】
図８では、検出された横隔膜１２ａの位置Ｐは「Ｐ１」である。横隔膜１２ａの位置を検出した後、ステップＳＴ４に進む。
【００３９】
ステップＳＴ４では、１回目のスキャンにおける補正用ｋ空間データ収集シーケンスＣａ１〜Ｃｎ１（図６参照）が実行される。
【００４０】
図９は、補正用ｋ空間データ収集シーケンスＣａ１〜Ｃｎ１を実行するときに設定されるスライスを示す図である。
【００４１】
図９に示すスライスＳＬａ〜ＳＬｎは、ステップＳＴ１においてオペレータ１３が決定したスライス位置、スライス厚、およびスライス枚数等に従って設定される。第１の実施形態では、スライスＳＬａ〜ＳＬｎの各々は、肝臓１２ｂの変位する方向（ｚ方向）に対して垂直に設定される。補正用ｋ空間データ収集シーケンスＣａ１〜Ｃｎ１は、図６に示すように、スライスＳＬａ〜ＳＬｎについて、補正用ｋ空間Ｓａ１〜Ｓｎ１の中心のデータＡ１〜Ｎ１を収集するためのシーケンスである。補正用ｋ空間Ｓａ１〜Ｓｎ１の中心のデータＡ１〜Ｎ１は、後述する２回目のスキャンのイメージングシーケンスＩＳにより収集されるｋ空間のデータＤＡ２〜ＤＮ２を補正するために使用されるデータである。
【００４２】
補正用ｋ空間データ収集シーケンスＣａ１〜Ｃｎ１は、例えば、図１０に示すように、ＲＦパルスの送信と、スライス選択勾配磁場Ｇｓの印加とを実行し、位相エンコード勾配磁場Ｇｐおよび周波数エンコード勾配磁場Ｇｆは印加しないパルスシーケンスによって実現することができる。
【００４３】
尚、肝臓１２ｂは被検体１２の呼吸に伴って変位するので、１回目の息止めをしたときの肝臓１２ｂの位置が、ステップＳＴ１で得られたスカウト画像における肝臓１２ｂの位置と同じであるとは限らない。したがって、スライスＳＬａ〜ＳＬｎのｚ方向の範囲Ｗ（図９参照）が狭すぎると、スカウト画像上では、肝臓１２ｂの全体を覆うようにスライスが設定されていても、１回目の息止めのときには、肝臓１２ｂの上端部や下端部が、スライスの範囲Ｗからはみ出てしまう恐れがある。そこで、オペレータ１３は、ステップＳＴ１でスライス位置等を決定する場合、１回目の息止めのときの肝臓１２ｂが、スライスＳＬａ〜ＳＬｎのｚ方向の範囲Ｗからはみ出ないように、スライスＳＬａ〜ＳＬｎをある程度広めに設定する必要がある。
【００４４】
補正用ｋ空間データ収集シーケンスＣａ１〜Ｃｎ１を実行した後、ステップＳＴ５に進む。
ステップＳＴ５では、１回目のスキャンにおけるイメージングシーケンスＩＳを複数回実行する。１回目のスキャンにおけるイメージングシーケンスＩＳは、スライスＳＬａ〜ＳＬｎについて、肝臓１２ｂを画像化するためのデータが配置されるｋ空間（以下、「画像用ｋ空間」と呼ぶ）ＳＡ〜ＳＮの中央領域のデータＤＡ１〜ＤＮ１を収集するためのシーケンスである。したがって、１回目のスキャンによって、画像用ｋ空間ＳＡ〜ＳＮの中央領域のデータＤＡ１〜ＤＮ１が収集される。
【００４５】
１回目のスキャンが終了したら、オペレータ１３は、被検体１２に対して、息止めの解除の指示をする。これによって、被検体１２は呼吸を再開することができる。被検体１２に呼吸を再開させた後、ステップＳＴ６に進む。
【００４６】
ステップＳＴ６では、オペレータ１３は、被検体１２に対して、２度目の息止めの指示をする。この指示によって、被検体１２は再び息止めをする。被検体１２に対して２回目の息止めの指示をした後、ステップＳＴ７に進む。
【００４７】
ステップＳＴ７では、２回目のスキャンのナビゲータシーケンスＮＡＶ（図６参照）が実行される。ナビゲータシーケンスＮＡＶを実行することによって、被検体１２が息止めをしているときの横隔膜１２ａの位置が検出される（図１１参照）。
【００４８】
図１１は、２回目のスキャンのナビゲータシーケンスＮＡＶによって検出された被検体１２の横隔膜１２ａの位置を概略的に示す図である。
【００４９】
図１１では、２回目の息止めのときの横隔膜１２ａおよび肝臓１２ｂは実線で示されており、１回目の息止めのときの横隔膜１２ａおよび肝臓１２ｂは破線で示されている。１回目の息止めでは、横隔膜１２ａの位置Ｐは、Ｐ＝Ｐ１であったが、２回目の息止めでは、横隔膜１２ａの位置Ｐは、Ｐ＝Ｐ２になっている。したがって、１回目の息止めと２回目の息止めとを比較すると、肝臓１２ｂの位置がずれていることがわかる。横隔膜１２ａの位置を検出した後、ステップＳＴ８に進む。
【００５０】
ステップＳＴ８では、２回目のスキャンにおける補正用ｋ空間データ収集シーケンスＣａ２〜Ｃｎ２を実行する。補正用ｋ空間データ収集シーケンスＣａ２〜Ｃｎ２は、図６に示すように、スライスＳＬａ〜ＳＬｎについて、補正用ｋ空間Ｓａ２〜Ｓｎ２の中心のデータＡ２〜Ｎ２を収集するためのシーケンスである。補正用ｋ空間Ｓａ２〜Ｓｎ２の中心のデータＡ２〜Ｎ２は、後述する２回目のスキャンのイメージングシーケンスＩＳにより収集されるｋ空間のデータＤＡ２〜ＤＮ２を補正するために使用されるデータである。補正用ｋ空間データ収集シーケンスＣａ２〜Ｃｎ２を実行することにより、スライスＳＬａ〜ＳＬｎについて、補正用ｋ空間Ｓａ２〜Ｓｎ２の中心のデータＡ２〜Ｎ２が収集される。しかし、１回目のスキャンと２回目のスキャンとを比較すると、横隔膜の位置にΔＰの差がある。そこで、２回目のスキャンにおける補正用ｋ空間データ収集シーケンスＣａ２〜Ｃｎ２を実行する場合、スライスＳＬａ〜ＳＬｎの位置を、１回目のスキャンにおけるスライスＳＬａ〜ＳＬｎに対して、差ΔＰだけ変位させる（図１２参照）。
【００５１】
図１２は、１回目のスキャンにおけるスライスと、２回目のスキャンにおけるスライスとの違いを示す図である。図１２では、２回目のスキャンにおける肝臓１２ｂとスライスが実線で示されており、１回目のスキャンにおける肝臓１２ｂとスライスとが破線で示されている。尚、図１２では、説明の便宜上、横隔膜１２ａは図示省略されている。
【００５２】
２回目のスキャンにおけるスライスＳＬａ〜ＳＬｎ（実線）の位置を、１回目のスキャンにおけるスライスＳＬａ〜ＳＬｎ（破線）に対して、差ΔＰだけ変位させることにより、スライスＳＬａ〜ＳＬｎによる肝臓１２ｂの励起位置を、１回目のスキャンと同じ励起位置にすることができる。
【００５３】
補正用ｋ空間Ｓａ２〜Ｓｎ２の中心のデータＡ２〜Ｎ２を収集した後、ステップＳＴ９に進む。
【００５４】
ステップＳＴ９では、２回目のスキャンにおけるイメージングシーケンスＩＳを複数回実行し、図６に示すように、１回目のスキャンで収集されなかった画像用ｋ空間ＳＡ〜ＳＮの残りのデータ（高周波領域のデータ）ＤＢ２〜ＤＮ２を収集する。
【００５５】
１回目のスキャンと２回目のスキャンによって、画像用ｋ空間ＳＡ〜ＳＮの全データを収集することができる。ただし、１回目のスキャンと２回目のスキャンとを比較すると、肝臓１２ｂの位置がｚ方向にΔＰずれているので、肝臓１２ｂと受信コイル４との相対的な位置関係も、ｚ方向にΔＰだけ変化している。したがって、１回目のスキャンと２回目のスキャンとを比較すると、肝臓１２ｂに対する受信コイル４の感度も変化する。このため、１回目のスキャンにより収集された画像用ｋ空間ＳＡ〜ＳＮの中央領域のデータＤＡ１〜ＤＮ１と、２回目のスキャンにより収集された画像用ｋ空間ＳＡ〜ＳＮの高周波領域のデータＤＡ２〜ＤＮ２とをそのまま使って各スライスＳＬａ〜ＳＬｎの画像を再構成すると、受信コイル４の感度の違いがアーチファクトとなって現れることがある。また、主磁場の不均一性もアーチファクトの原因となる。そこで、第１の実施形態では、アーチファクトを低減するために、２回目のスキャンで収集された画像用ｋ空間ＳＡ〜ＳＮの高周波領域のデータＤＡ２〜ＤＮ２を補正する。この補正は、ステップＳＴ１０およびＳＴ１１で行われる。
【００５６】
図１３は、ステップＳＴ１０およびＳＴ１１で実行される補正の説明図である。
ステップＳＴ１０では、強度補正が行われ、ステップＳＴ１１では位相補正が行われる。以下に、ステップＳＴ１０およびＳＴ１１の動作について順に説明する。
【００５７】
ステップＳＴ１０では、強度補正手段９１ａ（図１参照）が、２回目のスキャンで収集されたスライスＳＬａ〜ＳＬｎにおける画像用ｋ空間ＳＡ〜ＳＮの高周波領域のデータＤＡ２〜ＤＮ２の強度を補正する。例えばスライスＳＬａにおける画像用ｋ空間ＳＡの高周波領域のデータＤＡ２を補正する場合、１回目のスキャンで得られた補正用ｋ空間Ｓａ１の中心のデータＡ１と、２回目のスキャンで得られた補正用ｋ空間Ｓａ２の中心のデータＡ２とを使用する（図６参照）。画像用ｋ空間ＳＡの高周波領域のデータＤＡ２を強度補正した後のデータを「ＤＡ２′」とすると、強度補正後のデータＤＡ２′は、Ａ１およびＡ２を用いて、以下の式で表すことができる。
ＤＡ２′＝ＤＡ２×｜Ａ１｜／｜Ａ２｜・・・（１）
｜Ａ１｜：１回目のスキャンで得られた補正用ｋ空間Ｓａ１の中心のデータＡ１の強度
｜Ａ２｜：２回目のスキャンで得られた補正用ｋ空間Ｓａ２の中心のデータＡ２の強度
【００５８】
肝臓１２ｂが変位すると、受信コイルの感度が変化するので、ｋ空間のデータの強度が変動する。しかし、式（１）の｜Ａ１｜／｜Ａ２｜の値によって、ｋ空間のデータの強度の変動分を知ることができるので、｜Ａ１｜／｜Ａ２｜をＤＡ２に乗算することによって、ＤＡ２の強度を補正することができる。式（１）は、スライスＳＬａにおける画像用ｋ空間ＳＡの高周波領域のデータＤＡ２を強度補正する式を表している。しかし、他のスライスＳＬｂ〜ＳＬｎにおける画像用ｋ空間の高周波領域のデータＤＢ２〜ＤＮ２についても、１回目および２回目のスキャンで得られた補正用ｋ空間の中心のデータによって、データの強度の変動分を求めることができるので、強度補正をすることができる。例えば、スライスＳＬｂにおける画像用ｋ空間ＳＢの高周波領域のデータＤＢ２を強度補正する場合、強度補正後のデータＤＢ２′は、以下の式（２）で表すことができる。
ＤＢ２′＝ＤＢ２×｜Ｂ１｜／｜Ｂ２｜・・・（２）
｜Ｂ１｜：１回目のスキャンで得られた補正用ｋ空間Ｓｂ１の中心のデータＢ１の強度
｜Ｂ２｜：２回目のスキャンで得られた補正用ｋ空間Ｓｂ２の中心のデータＢ２の強度
【００５９】
また、スライスＳＬｎにおける画像用ｋ空間ＳＮの高周波領域のデータＤＮ２を強度補正する場合、強度補正後のデータＤＮ２′は、以下の式（３）で表すことができる。
ＤＮ２′＝ＤＮ２×｜Ｎ１｜／｜Ｎ２｜・・・（３）
｜Ｎ１｜：１回目のスキャンで得られた補正用ｋ空間Ｓｎ１の中心のデータＮ１の強度
｜Ｎ２｜：２回目のスキャンで得られた補正用ｋ空間Ｓｎ２の中心のデータＮ２の強度
【００６０】
したがって、図１３に示すように、ステップＳＴ１０では、スライスＳＬａ〜ＳＬｎにおける画像用ｋ空間ＳＡ〜ＳＮの高周波領域のデータＤＡ２〜ＤＮ２が、データＤＡ２′〜ＤＮ２′に強度補正される。強度補正を行った後、ステップＳＴ１１に進む。
【００６１】
ステップＳＴ１１では、位相補正手段９１ｂ（図１参照）が、画像用ｋ空間ＳＡ〜ＳＮの高周波領域の強度補正後のデータＤＡ２′〜ＤＮ２′の位相を補正する。例えば、スライスＳＬａにおける画像用ｋ空間ＳＡの高周波領域の強度補正後のデータＤＡ２′を位相補正する場合、強度補正と同様に、１回目のスキャンで得られた補正用ｋ空間Ｓａ１の中心のデータＡ１と、２回目のスキャンで得られた補正用ｋ空間Ｓａ２ｋ空間の中心のデータＡ２とを使用する（図６参照）。強度補正後のデータＤＡ２′を位相補正することにより得られるデータを、「ＤＡ２’’」とすると、位相補正後のデータＤＡ２’’は、以下の式で表すことができる。
ＤＡ２’’＝ＤＡ２′×ｅｘｐ（−ｉθａ）・・・（４）
θａ：データＡ１とＡ２との位相差
【００６２】
肝臓１２ｂが変位すると、受信コイルの感度が変化するので、ｋ空間のデータの位相が変動する。しかし、データＡ１とＡ２との位相差θａによって、ｋ空間のデータの位相の変動分を知ることができるので、ｅｘｐ（−ｉθａ）をＤＡ２′に乗算することによって、位相補正することができる。式（４）は、スライスＳＬａにおける強度補正後のデータＤＡ２′を位相補正する式を表している。しかし、他のスライスＳＬｂ〜ＳＬｎにおける強度補正後のデータＤＢ２′〜ＤＮ２′についても、１回目および２回目のスキャンで得られた補正用ｋ空間の中心のデータによって、位相補正をすることができる。例えば、スライスＳＬｂにおける画像用ｋ空間ＳＢの高周波領域の強度補正後のデータＤＢ２′を位相補正する場合、位相補正後のデータＤＢ２’’は、以下の式（５）で表すことができる。
ＤＢ２’’＝ＤＢ２′×ｅｘｐ（−ｉθｂ）・・・（５）
θｂ：データＢ１とＢ２との位相差
【００６３】
また、スライスＳＬｎにおける画像用ｋ空間ＳＮの高周波領域の強度補正後のデータＤＮ２′を位相補正する場合、位相補正後のデータＤＮ２’’は、以下の式（６）で表すことができる。
ＤＮ２’’＝ＤＮ２′×ｅｘｐ（−ｉθｎ）・・・（６）
θｎ：データＮ１とＮ２との位相差
【００６４】
したがって、図１３に示すように、ステップＳＴ１１では、スライスＳＬａ〜ＳＬｎにおける画像用ｋ空間ＳＡ〜ＳＮの高周波領域の強度補正後のデータＤＡ２′〜ＤＮ２′が、データＤＡ２’’〜ＤＮ２’’に位相補正される。位相補正を行った後、ステップＳＴ１２に進む。
【００６５】
ステップＳＴ１２では、再構成手段９２（図１参照）が、画像用ｋ空間ＳＡ〜ＳＮの中央領域のデータＤＡ１〜ＤＮ１と、強度補正および位相補正によって得られた高周波領域のデータＤＡ２’’〜ＤＮ２’’とを用いて画像を再構成する。このようにして、フローが終了する。
【００６６】
第１の実施形態では、１回目および２回目のスキャンにおいて、イメージングシーケンスＩＳを実行する前に、補正用ｋ空間の中心のデータを収集し、収集した補正用ｋ空間の中心のデータの強度および位相に基づいて、画像用ｋ空間ＳＡ〜ＳＮの高周波領域のデータＤＡ２〜ＤＮ２の強度および位相を補正している。したがって、肝臓１２ｂの変位によって受信コイルの感度が変化しても、アーチファクトの低減された画像を得ることができる。
【００６７】
尚、強度および位相の補正によってアーチファクトが低減されることを検証するため、強度および位相を補正した場合の画像と、強度および位相の補正をしない場合の画像とを取得し、画質を比較した。以下に、比較結果について説明する。
【００６８】
図１４は、比較結果の説明図である。
図１４（ａ）は、図６に示す本スキャンを実行することにより得られた画像用ｋ空間のデータを、強度および位相の補正をせずに再構成することにより得られたファントムの画像を示す図である。図１４（ｂ）は、図６に示す本スキャンを実行することにより得られた画像用ｋ空間のデータに対して強度および位相の補正を行い、強度および位相が補正されたデータを再構成することにより得られたファントムの画像を示す図である。尚、ファントムは、１回目のスキャンが終了した後、２回目のスキャンが開始される前に、１．７ｃｍだけ変位させている。
【００６９】
図１４（ａ）の画像を参照すると、矢印で示した部分に、アーチファクトが見られる。これに対して、図１４（ｂ）の画像では、アーチファクトが軽減されていることがわかる。
【００７０】
第１の実施形態では、画像用ｋ空間ＳＡ〜ＳＮの高周波領域のデータＤＡ２〜ＤＮ２の強度および位相を補正している。しかし、画像用ｋ空間ＳＡ〜ＳＮの高周波領域のデータＤＡ２〜ＤＮ２の強度および位相を補正する代わりに、画像用ｋ空間ＳＡ〜ＳＮの中央領域のデータＤＡ１〜ＤＮ１の強度および位相を補正してもよいし、画像用ｋ空間ＳＡ〜ＳＮの高周波領域のデータＤＡ２〜ＤＮ２と中央領域のデータＤＡ１〜ＤＮ１との両方について、強度および位相を補正してもよい。
【００７１】
第１の実施形態では、ナビゲータシーケンスＮＡＶによって、被検体１２の横隔膜の位置を検出している。しかし、ベローズで横隔膜の位置を検出してもよい。
【００７２】
第１の実施形態では、ｋ空間を中央領域と高周波領域とに分けているが、ｋ空間の分け方は、ｋ空間の正の領域と負の領域とに分けるなど、別の分け方でもよい。また、２回のスキャンに分けてｋ空間のデータを収集しているが、３回以上のスキャンに分けてｋ空間のデータを収集してもよいし、１回のスキャンでｋ空間のデータを収集してもよい。更に、スキャンは、息止めの状態で実行されているが、自由呼吸下でスキャンを実行してもよい。
【００７３】
第１の実施形態では、肝臓を撮影する場合について説明されているが、本発明は、被検体１２の体動により変位する部位であれば、肝臓以外の部位を撮影する場合にも適用できる。
【００７４】
第１の実施形態では、スライストラッキングが実行されているが、本発明は、必ずしもスライストラッキングを実行する必要はない。
【００７５】
（２）第２の実施形態
第１の実施形態では、イメージングシーケンスＩＳは、２Ｄ（Dimension）イメージング用のパルスシーケンスであるが、第２の実施形態では、イメージングシーケンスＩＳが、３Ｄイメージング用のパルスシーケンスであっても、補正用ｋ空間の中心のデータを用いて、強度補正および位相補正を行うことができる。第２の実施形態では、イメージングシーケンスＩＳが、３Ｄイメージング用のパルスシーケンスの場合に、補正用ｋ空間の中心のデータを用いて、強度補正および位相補正を行う例について説明する。イメージングシーケンスＩＳが、３Ｄイメージング用のパルスシーケンスであるので、ｋ空間は、３次元のｋ空間となる。
【００７６】
尚、第２の実施形態のＭＲＩ装置の構成は、第１の実施形態のＭＲＩ装置１００と同様に、図１で表すことができる。ただし、補正手段９１の動作については、第１の実施形態と相違する部分あるので、相違する部分については、後述する。
【００７７】
図１５は、第２の実施形態で実行されるスキャンによって収集されるｋ空間のデータを示す図である。
【００７８】
第２の実施形態でも、第１の実施形態と同様に、スカウトスキャンを実行し、その後に、本スキャンを実行する。スカウトスキャンは、第１の実施形態と同じであるので、スカウトスキャンの説明は省略し、本スキャンについて説明する。
【００７９】
第２の実施形態でも、第１の実施形態と同様に、本スキャンは、２回のスキャンを有している。
【００８０】
１回目のスキャンでは、ナビゲータシーケンスＮＡＶと、補正用ｋ空間データ収集シーケンスＣ１と、イメージングシーケンスＩＳとが実行される。ナビゲータシーケンスＮＡＶは、ナビゲータ領域Ｒnavからデータを収集し、横隔膜１２ａの位置Ｐを検出するためのシーケンスである。補正用ｋ空間データ収集シーケンスＣ１は、イメージングシーケンスＩＳの実行前に、補正用ｋ空間Ｓａ１の中心のデータＫc1を収集するためのシーケンスである。イメージングシーケンスＩＳは、撮影領域Ｒimを画像化するためのシーケンスであり、３Ｄイメージング用のパルスシーケンスである。イメージングシーケンスＩＳを複数回実行することによって、画像用ｋ空間ＳＡの中央領域のデータＤＫ１が収集される。
【００８１】
２回目のスキャンでは、ナビゲータシーケンスＮＡＶと、補正用ｋ空間データ収集シーケンスＣ２と、イメージングシーケンスＩＳとが実行される。２回目のスキャンでは、補正用ｋ空間データ収集シーケンスＣ２によって補正用ｋ空間Ｓａ２の中心のデータＫｃ２を収集した後、イメージングシーケンスＩＳを複数回を実行することによって、画像用ｋ空間ＳＡの高周波領域のデータＤＫ２が収集される。
【００８２】
次に、第２の実施形態における補正方法について説明する。
図１６は、第２の実施形態における補正方法の説明図である。
【００８３】
第２の実施形態では、先ず、強度補正手段９１ａ（図１参照）が、２回目のスキャンで収集された画像用ｋ空間ＳＡの高周波領域のデータＤＫ２の強度を補正する。画像用ｋ空間ＳＡの高周波領域のデータＤＫ２を強度補正した後のデータを「ＤＫ２′」とすると、強度補正後のデータＤＫ２′は、以下の式で表すことができる。
ＤＫ２′＝ＤＫ２×｜Ｋｃ１｜／｜Ｋｃ２｜・・・（７）
｜Ｋｃ１｜：１回目のスキャンで得られた補正用ｋ空間Ｓａ１の中心のデータＫｃ１の強度
｜Ｋｃ２｜：２回目のスキャンで得られた補正用ｋ空間Ｓａ２の中心のデータＫｃ２の強度
【００８４】
肝臓１２ｂが変位すると、受信コイルの感度が変化するので、ｋ空間のデータの強度が変動する。しかし、｜Ｋｃ１｜／｜Ｋｃ２｜の値によって、ｋ空間のデータの強度の変動分を知ることができるので、｜Ｋｃ１｜／｜Ｋｃ２｜をＤＫ２に乗算することによって、ＤＫ２の強度を補正することができる。
【００８５】
したがって、図１６に示すように、画像用ｋ空間ＳＡの高周波領域のデータＤＫ２が、データＤＫ２′に強度補正される。強度補正を行った後、位相補正手段９１ｂ（図１参照）が、画像用ｋ空間ＳＡの高周波領域の強度補正後のデータＤＫ２′の位相を補正する。強度補正後のデータＤＫ２′を位相補正することにより得られるデータを、「ＤＫ２’’」とすると、位相補正後のデータＤＫ２’’は、以下の式で表すことができる。
ＤＫ２’’＝ＤＫ２′×ｅｘｐ（−ｉθ）・・・（８）
θ：データＫｃ１とＫｃ２との位相差
【００８６】
肝臓１２ｂが変位すると、受信コイルの感度が変化するので、ｋ空間のデータの位相が変動する。しかし、データＫｃ１とＫｃ２との位相差θによって、ｋ空間のデータの位相の変動分を知ることができるので、ｅｘｐ（−ｉθ）をＤＫ２′に乗算することによって、位相補正することができる。
【００８７】
したがって、図１６に示すように、画像用ｋ空間ＳＡの高周波領域の強度補正後のデータＤＫ２′が、データＤＫ２’’に位相補正される。位相補正を行った後、画像用ｋ空間ＳＡの中央領域のデータＤＫ１と、強度補正および位相補正された高周波領域のデータＤＫ２’’とを用いて画像を再構成する。
【００８８】
上記のように、イメージングシーケンスＩＳが、３Ｄイメージング用のパルスシーケンスであっても、補正用ｋ空間の中心のデータＫｃ１およびＫｃ２を収集し、補正用ｋ空間の中心のデータＫｃ１およびＫｃ２を用いて、画像用ｋ空間ＳＡの高周波領域のデータＤＫ２の強度および位相を補正することにより、アーチファクトの低減された画像を得ることができる。
【００８９】
第２の実施形態では、先に強度補正し、次に位相補正を行っている。しかし、位相補正を先に行い、次に強度補正を行ってもよい。また、強度補正のみを行ってもよいし、位相補正のみを行ってもよい。
【００９０】
第２の実施形態では、画像用ｋ空間ＳＡの高周波領域のデータＤＫ２の強度および位相を補正している。しかし、画像用ｋ空間ＳＡの高周波領域のデータＤＫの強度および位相を補正する代わりに、画像用ｋ空間ＳＡの中央領域のデータＤＫ１の強度および位相を補正してもよいし、画像用ｋ空間ＳＡの高周波領域のデータＤＫと中央領域のデータＤＫ１との両方について、強度および位相を補正してもよい。
【００９１】
（３）第３の実施形態
第２の実施形態では、補正用ｋ空間の中心のデータを用いて、強度補正および位相補正を行う例について説明したが、第３の実施形態では、補正用ｋ空間の中心を横切り、ｋｚ方向に延在する直線上のデータを用いて、強度補正および位相補正を行う例について説明する。
【００９２】
図１７は、第３の実施形態のＭＲＩ装置３００を示す概略図である。
尚、第３の実施形態のＭＲＩ装置３００は、中央処理装置９が、第１の実施形態のＭＲＩ装置１００とは異なるが、その他は、第１の実施形態と同じである。したがって、第３の実施形態のＭＲＩ装置３００については、主に、中央処理装置９について説明する。
【００９３】
中央処理装置９は、逆フーリエ変換手段９０、補正手段９１、および再構成手段９２を有している。
【００９４】
逆フーリエ変換手段９０は、１回目および２回目のスキャンで収集されたｋ空間のデータＫｇ１、ＤＫ１およびＤＫ２、Ｋｇ２（後述する図１８参照）を、ｋｚ方向に１Ｄ逆フーリエ変換する。
【００９５】
補正手段９１は、強度補正手段９１ａおよび位相補正手段９１ｂを有している。強度補正手段９１ａは、１Ｄ逆フーリエ変換により得られた画像用データＤＫ１′およびＤＫ２′（後述する図１９参照）の強度を補正し、位相補正手段９１ｂは、画像用データＤＫ１′およびＤＫ２′の位相を補正する。
【００９６】
再構成手段９２は、強度補正および位相補正されたｋ空間のデータを用いて、画像を再構成する。
【００９７】
中央処理装置９は、逆フーリエ変換手段９０、補正手段９１、および再構成手段９２の一例であり、所定のプログラムを実行することにより、これらの手段として機能する。
【００９８】
次に、第３の実施形態において実行されるスキャンについて説明する。
図１８は、第３の実施形態におけるスキャンの説明図である。
【００９９】
第３の実施形態でも、第２の実施形態と同様に、スカウトスキャンを実行し、その後に、本スキャンを実行する。スカウトスキャンは、第２の実施形態と同じであるので、スカウトスキャンの説明は省略し、本スキャンについて説明する。
【０１００】
第３の実施形態でも、第２の実施形態と同様に、本スキャンは、２回のスキャンを有している。
【０１０１】
１回目のスキャンでは、ナビゲータシーケンスＮＡＶと、補正用ｋ空間データ収集シーケンスＧ１と、イメージングシーケンスＩＳとが実行される。ナビゲータシーケンスＮＡＶは、ナビゲータ領域Ｒnavからデータを収集し、横隔膜１２ａの位置Ｐを検出するためのシーケンスである。補正用ｋ空間データ収集シーケンスＧ１は、補正用ｋ空間Ｓａ１の中心を横切り、ｋｚ方向に延在する直線上のデータＫｇ１を収集するためのシーケンスである。イメージングシーケンスＩＳは、撮影領域Ｒimを画像化するためのシーケンスであり、３Ｄイメージング用のパルスシーケンスである。イメージングシーケンスＩＳを複数回実行することによって、画像用ｋ空間ＳＡの中央領域のデータＤＫ１が収集される。
【０１０２】
２回目のスキャンでは、ナビゲータシーケンスＮＡＶと、補正用ｋ空間データ収集シーケンスＧ２と、イメージングシーケンスＩＳとが実行される。２回目のスキャンでは、補正用ｋ空間データ収集シーケンスＧ２によって補正用ｋ空間Ｓａ２の直線上のデータＫｇ２を収集した後、イメージングシーケンスＩＳを複数回を実行することによって、画像用ｋ空間ＳＡの高周波領域のデータＤＫ２が収集される。
【０１０３】
第３の実施形態における補正方法について説明する。
図１９は、第３の実施形態における補正方法の説明図である。
【０１０４】
第３の実施形態では、先ず、逆フーリエ変換手段９０（図１７参照）が、１回目および２回目のスキャンで収集されたｋ空間のデータＫｇ１、ＤＫ１およびＤＫ２、Ｋｇ２を、ｋｚ方向に１Ｄ逆フーリエ変換する。ｋｚ方向に１Ｄ逆フーリエ変換することによって、補正用ｋ空間Ｓａ１のデータＫｇ１は、ｋｚ軸が実空間上のｚ軸に変換された補正用データＫｇ１′に変換される。また、画像用ｋ空間ＳＡのデータＤＫ１およびＤＫ２は、ｋｚ軸が実空間上のｚ軸に変換された画像用データＤＫ１′およびＤＫ２′に変換される。更に、補正用ｋ空間Ｓａ２のデータＫｇ２は、ｋｚ軸が実空間上のｚ軸に変換された補正用データＫｇ２′に変換される。したがって、ｋｚ方向に１Ｄ逆フーリエ変換することによって、第１の実施形態で説明した図６に示すように、スライスＳＬａ〜ＳＬｎごとのｋ空間のデータを得ることができる。ｋｚ方向に１Ｄ逆フーリエ変換した後は、第１の実施形態と同様の手順で、補正手段９１が強度補正および位相補正を実行し、再構成手段９２が画像を再構成することにより、アーチファクトの低減された画像を得ることができる。
【０１０５】
尚、第１〜第３の実施形態では、先に強度補正し、次に位相補正を行っている。しかし、位相補正を先に行い、次に強度補正を行ってもよい。また、強度補正のみを行ってもよいし、位相補正のみを行ってもよい。
【符号の説明】
【０１０６】
１００  ＭＲＩ装置
２  磁場発生装置
３  テーブル
４  受信コイル
５  シーケンサ
６  送信器
７  勾配磁場電源
８  受信器
９  中央処理装置
１０  入力装置
１１  表示装置
１２  被検体
１３  オペレータ
２１  ボア
２２  超伝導コイル
２３  勾配コイル
２４  送信コイル
３１  クレードル
９１  補正手段
９１ａ  強度補正手段
９１ｂ  位相補正手段
９２  再構成手段
【図１】
0005674356.tif 000002
【図２】
0005674356.tif 000003
【図３】
0005674356.tif 000004
【図４】
0005674356.tif 000005
【図５】
0005674356.tif 000006
【図６】
0005674356.tif 000007
【図８】
0005674356.tif 000008
【図９】
0005674356.tif 000009
【図１０】
0005674356.tif 000010
【図１１】
0005674356.tif 000011
【図１２】
0005674356.tif 000012
【図１３】
0005674356.tif 000013
【図１５】
0005674356.tif 000014
【図１６】
0005674356.tif 000015
【図１７】
0005674356.tif 000016
【図１８】
0005674356.tif 000017
【図１９】
0005674356.tif 000018
【図７】
0005674356.tif 000019
【図１４】
0005674356.tif 000020