(57)【特許請求の範囲】
【請求項１】
ｋ空間の第１の領域に含まれる格子点のデータと、前記第１の領域に隣接する第２の領域に含まれる格子点のデータとを収集するためのスキャンを実行する磁気共鳴装置であって、
ｋ空間の格子点に、格子点のデータが収集される順番を表す番号を付与するナンバリング手段と、
前記第１の領域と前記第２の領域との境界を設定する設定手段と、
前記第１の領域のうちの前記第２の領域側に位置する第１の複数の格子点のうちの一部の格子点に付与された番号と、前記第２の領域のうちの前記第１の領域側に位置する第２の複数の格子点のうちの一部の格子点に付与された番号とを交換する交換手段と、
前記スキャンを実行するスキャン手段であって、
ｋ空間の格子点のデータの収集が行われる第１のデータ収集期間において、前記第１の領域の中で、前記交換手段により番号の交換が行われていない格子点のデータと、前記第２の領域の中で、前記交換手段により番号の交換が行われた格子点のデータとを収集し、
ｋ空間の格子点のデータの収集が行われる第２のデータ収集期間において、前記第１の領域の中で、前記交換手段により番号の交換が行われた格子点のデータと、前記第２の領域の中で、前記交換手段により番号の交換が行われていない格子点のデータとを収集するスキャン手段と、
を有する磁気共鳴装置。
【請求項２】
前記第１の領域の中で番号の交換が行われていない格子点と前記第２の領域の中で番号の交換が行われた格子点とを含む第３の複数の格子点の番号を並び替えるソート手段を有する、請求項１に記載の磁気共鳴装置。
【請求項３】
前記スキャン手段は、
前記第１のデータ収集期間では、前記ソート手段により番号が並び替えられた後の前記第３の複数の格子点の番号に従って、前記第３の複数の格子点の各々のデータを収集する、請求項２に記載の磁気共鳴装置。
【請求項４】
前記ソート手段は、
格子点がｋ空間の中心から離れるに従って番号が大きくなるように、前記第３の複数の格子点の番号を並び替える、請求項２又は３に記載の磁気共鳴装置。
【請求項５】
前記交換手段は、
前記第１の領域の中から格子点の番号がα番づつ小さくなるように選択された各格子点の番号と、前記第２の領域の中から格子点の番号がβ番づつ大きくなるように選択された各格子点の番号とを交換する、請求項１〜４のうちのいずれか一項に記載の磁気共鳴装置。
【請求項６】
前記交換手段は、
前記第１の領域の各格子点と前記境界との距離に応じてαの値を変更し、前記第２の領域の各格子点と前記境界との距離に応じてβの値を変更する、請求項５に記載の磁気共鳴装置。
【請求項７】
αの値はβの値に等しい、請求項５又は６に記載の磁気共鳴装置。
【請求項８】
前記ナンバリング手段は、ｋ空間の格子点とｋ空間の中心との距離に基づいて、格子点に番号を付与する、請求項１〜７のうちのいずれか一項に記載の磁気共鳴装置。
【請求項９】
前記スキャンは、前記第１のデータ収集期間と同じ格子点のデータを収集するための他の第１のデータ収集期間を有する、請求項１〜８のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１０】
前記スキャン手段は、
前記第２のデータ収集期間は、複数のデータ収集期間に分けられており、
前記複数のデータ収集期間は、互いに異なる格子点のデータを収集する、請求項１〜９のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１１】
前記第１の領域はｋ空間の中心を含む低周波領域であり、
前記第２の領域は高周波領域である、請求項１〜１０のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１２】
前記スキャン手段は、
前記第１のデータ収集期間において、前記第１の複数の格子点のうちの一部の格子点のデータと、前記第２の複数の格子点のうちの一部の格子点のデータとを収集するための第１のシーケンスを実行し、
前記第２のデータ収集期間において、前記第１の複数の格子点のうちの、前記第１のデータ収集期間でデータが収集される格子点とは別の格子点のデータと、前記第２の複数の格子点のうちの、前記第１のデータ収集期間でデータが収集される格子点とは別の格子点のデータとを収集するための第２のシーケンスを実行する、請求項１〜１１のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１３】
前記スキャン手段は、
前記第１のデータ収集期間では、脂肪抑制パルスを印加した後で、前記第１のシーケンスを実行し、
前記第２のデータ収集期間では、脂肪抑制パルスを印加した後で、前記第２のシーケンスを実行する、請求項１２に記載の磁気共鳴装置。
【請求項１４】
前記第１のデータ収集期間又は前記第２のデータ収集期間では、複数の脂肪抑制パルスが印加される、請求項１３に記載の磁気共鳴装置。
【請求項１５】
前記スキャン手段は、
前記第１のデータ収集期間では、被検体の呼吸情報を含むナビゲータエコーを収集するための第１のナビゲータシーケンスを実行し、前記第１のナビゲータシーケンスを実行した後で、前記第１のシーケンスを実行し、  前記第２のデータ収集期間では、被検体の呼吸情報を含むナビゲータエコーを収集するための第２のナビゲータシーケンスを実行し、前記第２のナビゲータシーケンスを実行した後で、前記第２のシーケンスを実行する、請求項１２〜１４のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１６】
前記スキャン手段は、
前記第１のデータ収集期間では、血流信号を低下させるための第１の血流抑制シーケンスを実行し、前記第１の血流抑制シーケンスを実行した後で、前記第１のシーケンスを実行し、
前記第２のデータ収集期間では、血流信号を低下させるための第２の血流抑制シーケンスを実行し、前記第２の血流抑制シーケンスを実行した後で、前記第２のシーケンスを実行する、請求項１２〜１５のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１７】
前記スキャン手段は、造影剤が投与された被検体のスキャンを実行する、請求項１〜１６のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１８】
前記第１のデータ収集期間により収集されたｋ空間の格子点のデータと前記第２のデータ収集期間により収集されたデータとに基づいて画像を作成する画像作成手段を有する、請求項１〜１７のうちのいずれか一項に記載の磁気共鳴装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ｋ空間のデータを収集する磁気共鳴装置に関する。
【背景技術】
【０００２】
近年、造影剤を用いて各時相の画像を収集するダイナミックＭＲイメージングが普及している（例えば、特許文献１参照）。造影剤を用いたダイナミックＭＲイメージングでは、撮影部位における造影剤の濃度の時間変化を観察することが重要となるので、高い時間分解能で撮影を行うことが要求される。このような要求に応える手法として、ダイナミック撮像時には、ｋ空間の中心付近の低周波領域のデータのみを更新するキーホールイメージング法（例えば、非特許文献１参照）や、高周波領域をＮ個の領域に分割し、各領域のデータと低周波領域のデータとを交互に収集するＤＩＳＣＯと呼ばれる手法が知られている（例えば、非特許文献２参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１３−１７６６７３号公報
【非特許文献】
【０００４】
【非特許文献１】R.A.Jones et al., “Dynamic, contast enhanced, NMR perfusion imaging of regional cerebral ischaemia in rats using k space substitution,” SMRM 1992,p.1138.
【非特許文献２】Saranathan M, et al., “DIfferential Subsampling with Cartesian Ordering (DISCO): a high spatio-temporal resolution Dixon imaging sequence for multiphasic contrast enhanced abdominal imaging”, J Magn Reson Imging. 2012 35(6):1484-92.
【発明の概要】
【発明が解決しようとする課題】
【０００５】
造影剤を用いたダイナミックＭＲＩイメージングでは、脂肪が高信号で描出されてしまうと、画像診断の妨げになることがある。したがって、造影剤を用いたダイナミックＭＲＩイメージングでは、脂肪抑制を行うことも重要となる。そこで、脂肪抑制法を併用したイメージングが実行されている。脂肪抑制法としては、脂肪抑制パルスを用いた方法がある。しかし、脂肪抑制パルスを用いた場合、ｋ空間の低周波領域と高周波領域との境界付近で信号強度に大きな差が現れやすくなり、画像に、信号強度の差によるアーチファクトが現れるという問題がある。したがって、アーチファクトが低減された画像を取得するのに適したイメージング法が望まれている。
【課題を解決するための手段】
【０００６】
本発明の一観点は、ｋ空間の第１の領域に含まれる格子点のデータと、前記第１の領域に隣接する第２の領域に含まれる格子点のデータとを収集するためのスキャンを実行するスキャン手段を有する磁気共鳴装置であって、
前記スキャン手段は、
ｋ空間の格子点のデータの収集が行われる第１のデータ収集期間において、前記第１の領域のうちの前記第２の領域側に位置する第１の複数の格子点のうちの一部の格子点のデータと、前記第２の領域のうちの前記第１の領域側に位置する第２の複数の格子点のうちの一部の格子点のデータとを収集し、
ｋ空間の格子点のデータの収集が行われる第２のデータ収集期間において、前記第１の複数の格子点のうちの、前記第１のデータ収集期間でデータが収集される格子点とは別の格子点のデータと、前記第２の複数の格子点のうちの、前記第１のデータ収集期間でデータが収集される格子点とは別の格子点のデータとを収集する、磁気共鳴装置である。
【発明の効果】
【０００７】
第１の領域と第２の領域の境界付近において、第１のデータ収集期間で収集される格子点のデータと、第２のデータ収集期間で収集される格子点のデータを分散させることができる。したがって、第１の領域と第２の領域との境界付近では、信号強度の大きいデータと信号強度の小さいデータが分散するように配置されるので、信号強度の差に起因したアーチファクトを低減することができる。
【図面の簡単な説明】
【０００８】
【図１】本発明の第１の形態の磁気共鳴装置の概略図である。
【図２】プロセッサ９が実行する処理を示す図である。
【図３】キーホールイメージングの一例の説明図である。
【図４】ｋ空間を低周波領域と高周波領域に分けるときの一例を概略的に示す図である。
【図５】データ収集期間Ａ１、Ｂ、Ａ２およびＡ３で収集されるｋ空間のデータを示す図である。
【図６】脂肪抑制法を適用する例の説明図である。
【図７】第１の形態におけるダイナミックＭＲイメージングを示すフローを示す図である。
【図８】ステップＳＴ１の説明図である。
【図９】各格子点に付与された番号を示す図である。
【図１０】低周波領域と高周波領域との境界Ｌの一例を示す図である。
【図１１】番号が交換された後のｋ空間を概略的に示す図である。
【図１２】番号が交換された後のｋ空間を概略的に示す図である。
【図１３】番号が交換された後のｋ空間を概略的に示す図である。
【図１４】１≦ｎ≦１６の間で格子点の番号を交換したときの様子を概略的に示す図である。
【図１５】ステップＳＴ２において実行されるスキャンＳＣの説明図である。
【図１６】ｋ空間データＤ１の拡大図である。
【図１７】ｋ空間データＤ１の中央部分Ｒ１を示す図である。
【図１８】ｋ空間データＤ１の外縁部分Ｒ２を示す図である。
【図１９】境界Ｌの近傍部分Ｒ３を示す図である。
【図２０】データ収集期間Ｂに複数の脂肪抑制パルスを印加した例を示す図である。
【図２１】格子点４５〜格子点２２５のデータを複数のデータ収集期間に分けて収集する一例の説明図である。
【図２２】第２の形態におけるプロセッサ９が実行する処理の説明図である。
【図２３】第２の形態におけるステップＳＴ１のフローを示す図である。
【図２４】ソート手段９５による番号の並替えが終了した後の格子点１〜格子点４５の番号を示す図である。
【図２５】データ収集期間Ａ１に複数の脂肪抑制パルスを印加した例を示す図である。
【図２６】格子点１〜格子点４５のデータを奇数番目の格子点と偶数番目の格子点に分けて収集する一例を示す図である。
【発明を実施するための形態】
【０００９】
以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。
【００１０】
（１）第１の形態
図１は、本発明の第１の形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ）１００は、マグネット２、テーブル３、受信ＲＦコイル４などを有している。
【００１１】
マグネット２は、被検体１３が収容されるボア２１を有している。また、マグネット２には、超伝導コイル、勾配コイル、およびＲＦコイルなどが内蔵されている。
【００１２】
テーブル３は、被検体１３を支持するクレードル３ａを有している。クレードル３ａは、ボア２１内に移動できるように構成されている。クレードル３ａによって、被検体１３はボア２１に搬送される。
受信ＲＦコイル（以下、「受信コイル」と呼ぶ）４は、被検体１３からの磁気共鳴信号を受信する。
【００１３】
ＭＲ装置１００は、更に、送信器５、勾配磁場電源６、受信器７、コンピュータ８、操作部１１、および表示部１２などを有している。
【００１４】
送信器５はＲＦコイルに電流を供給し、勾配磁場電源６は勾配コイルに電流を供給する。受信器７は、受信コイル４から受け取った信号に対して、検波などの信号処理を実行する。尚、マグネット２、受信コイル４、送信器５、勾配磁場電源６、および受信器７を合わせたものがスキャン手段に相当する。
【００１５】
コンピュータ８は、表示部１２に必要な情報を伝送したり、画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。コンピュータ８は、プロセッサ９およびメモリ１０などを有している。
【００１６】
図２に、プロセッサ９が実行する処理の説明図である。メモリ１０には、プロセッサ９により実行されるプログラムが記憶されている。プロセッサ９は、メモリ１０に記憶されているプログラムを読み出し、プログラムに記述されている処理を実行する。プロセッサ９は、メモリ１０に記憶されているプログラムを読み出すことにより、ナンバリング手段９１〜画像作成手段９４などを構成する。
【００１７】
ナンバリング手段９１は、ｋ空間の格子点に番号を付与する。
設定手段９２は、高周波領域と低周波領域の境界を設定する。
交換手段９３は、高周波領域に含まれる一部の格子点の番号と、低周波領域に含まれる一部の格子点の番号とを交換する。
画像作成手段９４は、受信器７から受け取ったデータに基づいて、撮影部位の画像を作成する。
【００１８】
プロセッサ９は、ナンバリング手段９１〜画像作成手段９４を構成する一例であり、メモリ１０に記憶されたプログラムを実行することによりこの手段として機能する。尚、ナンバリング手段９１、設定手段９２、および交換手段９３を合わせたものが、決定手段に相当する。
【００１９】
操作部１１は、オペレータにより操作され、種々の情報をコンピュータ８に入力する。表示部１２は種々の情報を表示する。
ＭＲ装置１００は、上記のように構成されている。
【００２０】
また、図１には、被検体に造影剤を投与する造影剤投与装置１０１も示されている。
第１の形態では、造影剤投与装置１０１を用いて被検体１３に造影剤を投与し、ダイナミックＭＲイメージングを行う例について説明する。尚、以下の説明では、第１の形態におけるダイナミックＭＲイメージングの効果を明確にするために、第１の形態におけるダイナミックＭＲイメージングを説明する前に、ダイナミックＭＲイメージングで使用されている一手法であるキーホールイメージングについて説明する。
【００２１】
図３は、キーホールイメージングの一例の説明図である。
図３では、ｋ空間のデータを収集するための複数のデータ収集期間が示されている。ここでは、説明の便宜上、造影剤が投与された後の４つのデータ収集期間Ａ１、Ｂ、Ａ２およびＡ３のみを示している。各データ収集期間では、シーケンスＶが複数回実行される。図３には、シーケンスＶの一例として、３Ｄグラディエントエコー系のシーケンスが示されている。キーホールイメージングでは、ｋ空間を低周波領域と高周波領域に分けて、データ収集を行う。図４に、ｋ空間を低周波領域と高周波領域に分けるときの一例を概略的に示す。図４では、ｋ空間のｋｘ方向は図示省略し、ｋｙ−ｋｚ面が示されている。キーホールイメージングでは、ｋｙ―ｋｚ面を、低周波領域ＲＬと高周波領域ＲＨに分ける。説明の便宜上、低周波領域ＲＬには４９個の格子点（黒丸で示されている）が含まれており、高周波領域ＲＨには、１７６個の格子点（白丸で示されている）が含まれているとする。キーホールイメージングでは、図４に示すように、ｋ空間を低周波領域ＲＬと高周波領域ＲＨに分けて、データ収集期間Ａ１、Ｂ、Ａ２およびＡ３ごとにデータ収集を行う（図５参照）。
【００２２】
図５は、データ収集期間Ａ１、Ｂ、Ａ２およびＡ３で収集されるｋ空間のデータを示す図である。
データ収集期間Ａ１では、ｋ空間の低周波領域ＲＬ内の４９個の格子点（黒丸）の各々のデータが収集されるように、シーケンスＶの位相エンコード勾配磁場が調整される。図５では、データ収集期間Ａ１において収集された低周波領域のデータを、符号「ＤＬ１」で示してある。低周波領域のデータＤＬ１を収集した後、次のデータ収集期間Ｂに移行する。
【００２３】
データ収集期間Ｂでは、ｋ空間の高周波領域ＲＨ内の１７６個の格子点（白丸）の各々のデータが収集されるように、シーケンスＶの位相エンコード勾配磁場が調整される。図５では、データ収集期間Ｂにおいて収集された高周波領域のデータを、符号「ＤＨ」で示してある。したがって、データ収集期間Ａ１で得られたデータＤＬ１と、データ収集期間Ｂで得られたデータＤＨにより、全ての格子点のデータを含むｋ空間データＤ１を得ることができる。ｋ空間データＤ１を得た後、フーリエ変換を行うことにより、画像ＩＭ１を作成することができる。データ収集期間Ｂが終了した後、次のデータ収集期間Ａ２に移行する。
【００２４】
データ収集期間Ａ２では、データ収集期間Ａ１と同様に、ｋ空間の低周波領域ＲＬ内の４９個の格子点（黒丸）の各々のデータが収集されるように、シーケンスＶの位相エンコード勾配磁場が調整される。図５では、データ収集期間Ａ２において収集された低周波領域ＲＬのデータを、符号「ＤＬ２」で示してある。データ収集期間Ｂで得られたデータＤＨと、データ収集期間Ａ２で得られたデータＤＬ２により、２番目のｋ空間データＤ２を得ることができる。ｋ空間データＤ２を得た後、フーリエ変換を行うことにより、画像ＩＭ２を作成することができる。データ収集期間Ａ２が終了した後、次のデータ収集期間Ａ３に移行する。
【００２５】
データ収集期間Ａ３では、データ収集期間Ａ１と同様に、ｋ空間の低周波領域ＲＬ内の４９個の格子点（黒丸）の各々のデータが収集されるように、シーケンスＶの位相エンコード勾配磁場が調整される。図５では、データ収集期間Ａ３において収集された低周波領域ＲＬのデータを、符号「ＤＬ３」で示してある。データ収集期間Ｂで得られたデータＤＨと、データ収集期間Ａ３で得られたデータＤＬ３により、３番目のｋ空間のデータＤ３を得ることができる。ｋ空間のデータＤ３を得た後、フーリエ変換を行うことにより、画像ＩＭ３を作成することができる。
【００２６】
キーホールイメージングでは、データ収集期間Ｂで収集されるｋ空間の高周波領域ＲＨのデータＤＨが、画像ＩＭ１、ＩＭ２、およびＩＭ３を作成するために使用されるｋ空間の高周波領域のデータとして共通に使用される。したがって、データ収集期間Ａ１、Ａ２、およびＡ３の間は、低周波領域ＲＬのデータのみを収集すればよいので、短い時間間隔で画像ＩＭ１、ＩＭ２、およびＩＭ３を得ることができ、造影剤を用いたダイナミックＭＲＩイメージングの時間分解能を向上させることができる。
【００２７】
ただし、造影剤を用いたダイナミックＭＲＩイメージングでは、脂肪が高信号で描出されてしまうと、画像診断の妨げになることがある。したがって、造影剤を用いたダイナミックＭＲＩイメージングでは、脂肪抑制を行うことも重要となる。そこで、図５に示したスキャンＳＣに対して、脂肪抑制法を適用する例について考える（図６参照）。
【００２８】
図６は、脂肪抑制法を適用する例の説明図である。
図６では、脂肪抑制パルスＦを用いて脂肪抑制を行う例が示されている。以下に、脂肪抑制パルスＦの脂肪抑制効果について説明する。
データ収集期間Ａ１では、先ず、脂肪抑制パルスＦが印加される。脂肪抑制パルスＦが印加された後、シーケンスＶが実行され、低周波領域ＲＬのデータＤＬ１が収集される。脂肪抑制を行う場合、脂肪抑制効果を高めるために、一般的には、脂肪抑制パルスＦを印加したら、ｋ空間の中心Ｃに近い格子点からデータを収集する。したがって、データ収集期間Ａ１では、脂肪抑制パルスＦを印加した直後は、低周波領域の中心Ｃの格子点のデータが収集される。そして、中心Ｃに近い格子点から順にデータが収集され、データ収集期間Ａ１の終了時点に近づくに従い、低周波領域と高周波領域との境界Ｌの近傍に位置する格子点のデータが収集される。したがって、ｋ空間の中心Ｃの近傍では、脂肪抑制パルスＦの脂肪抑制効果が高いので、脂肪の信号強度は小さい（脂肪が十分に抑制されている）が、境界Ｌの近傍では、脂肪抑制パルスＦの脂肪抑制効果が低くなり、ｋ空間の中心Ｃの近傍よりも脂肪の信号強度が大きくなる。低周波領域ＲＬのデータを収集した後、次のデータ収集期間Ｂに移行する。
【００２９】
データ収集期間Ｂでは、データ収集期間Ａ１と同様に、先ず、脂肪抑制パルスＦが印加される。脂肪抑制パルスＦが印加された後、シーケンスＶが実行され、高周波領域ＲＨのデータＤＨが収集される。データ収集期間Ｂでも、脂肪抑制パルスＦの脂肪抑制効果を高めるために、低周波領域側の格子点（つまり、境界Ｌの近傍に位置する格子点）からデータが収集される。そして、データ収集期間Ｂの終了時点に近づくに従い、境界Ｌから離れた位置における格子点のデータが収集される。したがって、境界Ｌの近傍では、脂肪抑制パルスＦの脂肪抑制効果が高いので、脂肪の信号強度は小さい（脂肪が十分に抑制されている）が、境界Ｌから離れるに従って脂肪抑制パルスＦの脂肪抑制効果が低くなり、脂肪の信号強度が大きくなる。
【００３０】
データ収集期間Ａ１では低周波領域のデータＤＬ１が収集され、データ収集期間Ｂでは高周波領域のデータＤＨが収集されるので、全ての格子点のデータを含むｋ空間データＤ１を得ることができる。しかし、低周波領域のデータＤＬ１は、境界Ｌの近傍では脂肪の信号強度が大きく、一方、高周波領域のデータＤＨは、境界Ｌの近傍では脂肪の信号強度が小さいので、ｋ空間データＤ１の境界Ｌに、信号強度が急激に変化する段差が現れる。ｋ空間にこのような信号強度の段差が現れると、画像ＩＭ１にゴーストが発生してしまい、画質が劣化するという問題がある。
【００３１】
高周波領域ＲＨのデータを収集した後、データ収集期間Ａ２およびＡ３においてデータ収集が実行される。データ収集期間Ａ２およびＡ３でも、脂肪抑制パルスＦが印加され、シーケンスＶが実行される。したがって、ｋ空間データＤ２およびＤ３でも、ｋ空間のデータＤ１と同様に、境界Ｌにおいて信号強度の段差が現れる。このため、画像ＩＭ２および画像ＩＭ３にもゴーストが発生してしまい、画質が劣化するという問題がある。
【００３２】
そこで、第１の形態では、脂肪抑制パルスＦを用いたスキャンであってもゴーストが低減されるように、ダイナミックＭＲイメージングを実行している。以下に、第１の形態におけるダイナミックＭＲイメージングについて説明する。
【００３３】
図７は、第１の形態におけるダイナミックＭＲイメージングのフローを示す図である。
第１の形態のフローは、２つのステップＳＴ１およびＳＴ２に分けられる。ステップＳＴ１は、ｋ空間の格子点のデータの収集順序を決定するステップである。ステップＳＴ２は、ステップＳＴ１で決定したデータ収集順序に基づいてスキャンを実行するステップである。以下に、ステップＳＴ１およびＳＴ２について順に説明する。
【００３４】
図８は、ステップＳＴ１の説明図である。
ステップＳＴ１は、ステップＳＴ１１〜ＳＴ１３を有している。以下、ステップＳＴ１１〜ＳＴ１３について順に説明する。
【００３５】
ステップＳＴ１１では、ナンバリング手段９１（図２参照）が、ｋ空間の格子点に、ｋ空間の中心Ｃから近い順に１〜ｎの番号を付与する。図９に、各格子点に付与された番号を示す。第１の形態では、説明の便宜上、ｋ空間の格子点の数は２２５個であるとする。したがって、ｋ空間の格子点には、ｋ空間の中心Ｃから近い順に、１〜２２５の番号が付される。
【００３６】
図９では、一つの格子点を一つの四角形で示してある。四角形の中に記載されている数字は、各格子点に付与された番号を表している。番号「１」は、ｋ空間の中心Ｃ上の格子点に付与されている。そして、ｋ空間の各格子点とｋ空間の中心Ｃにおける格子点（番号「１」）との距離ｄを求め、距離ｄの短い順に、番号「２」、「３」、「４」、・・・「２２５」が付与されている。したがって、距離ｄが短いほど（格子点がｋ空間の中心Ｃに近いほど）、格子点に付与される番号は小さいが、距離ｄが長いほど（格子点が中心Ｃから離れているほど）、格子点に付与される番号は大きくなる。以下では、説明の便宜上、各格子点を、各格子点に付された番号で呼ぶことにする（例えば、番号「１」が付された格子点は、「格子点１」と呼ぶことにする）。尚、ｋ空間の格子点の中には、距離ｄが同じ値になる格子点がある（例えば、格子点２、格子点３、格子点４、および格子点５）。このような場合は、所定の条件（例えば、ｋｚ座標やｋｙ座標の座標値の違い）に基づいて、番号を付与する。したがって、２２５個の格子点に、異なる番号を付与することができる。各格子点に番号を付与した後、ステップＳＴ１２に進む。
【００３７】
ステップＳＴ１２では、設定手段９２（図２参照）が、ｋ空間上に、低周波領域と高周波領域との境界を設定する。図１０に、低周波領域と高周波領域との境界Ｌの一例を示す。境界Ｌは、４５個の格子点１〜格子点４５を囲むように設定されており、格子点４６〜格子点２２５は、境界Ｌの外側に位置している。図１０では、境界Ｌの内側の格子点と、境界Ｌの外側の格子点とを区別するために、境界Ｌの内側の格子点は、黒丸で示してある。境界Ｌを設定した後、ステップＳＴ１３に進む。
【００３８】
ステップＳＴ１３では、交換手段９３（図２参照）が、境界Ｌの内側の領域（低周波領域ＲＬ）に位置する４５個の格子点１〜格子点４５の一部の格子点の番号を、境界Ｌの外側の領域（高周波領域ＲＨ）に位置する格子点の番号と交換する。以下に、交換の方法について説明する。
【００３９】
第１の形態では、交換手段９３は、先ず、境界Ｌの内側の領域（低周波領域ＲＬ）の中から、以下の式で規定される番号Ｎｐが付与された格子点Ｎｐを特定する。
Ｎｐ＝４５−α（ｎ−１）  ・・・（１）
ここで、α：係数
ｎ：１≦ｎ≦ｎｕの範囲の整数
【００４０】
次に、交換手段９３は、境界Ｌの外側の領域（高周波領域ＲＨ）の中から、以下の式で規定される番号Ｎｑが付与された格子点Ｎｑを特定する。
Ｎｑ＝４５＋β・ｎ  ・・・（２）
ここで、β：係数
ｎ：１≦ｎ≦ｎｕの範囲の整数
【００４１】
尚、ｎの上限を規定するｎｕは、以下の式を満たすように設定される。
４５−α（ｎｕ−１）≧１  ・・・（３）
４５＋β・ｎｕ≦２２５    ・・・（４）
【００４２】
交換手段９３は、式（１）で決定される格子点Ｎｐの番号と、式（２）で決定される格子点Ｎｑの番号との交換を行う。尚、以下では、説明の便宜上、式（１）および（２）において、α＝β＝２とし、ｎｕ＝１６であるとする。したがって、式（１）のＮｐおよび式（２）のＮｑは、以下の式で表される。
Ｎｐ＝４５−２（ｎ−１）  ・・・（５）
ここで、ｎ：１≦ｎ≦１６の範囲の整数
Ｎｑ＝４５＋２・ｎ  ・・・（６）
ここで、ｎ：１≦ｎ≦１６の範囲の整数
【００４３】
交換手段９３は、式（５）および（６）のｎに、ｎ＝１、２、３、・・・１６を代入し、格子点Ｎｐの番号と格子点Ｎｑの番号との交換を行う。具体的には、以下の手順で格子点の番号を交換する。
【００４４】
交換手段９３は、先ず、式（５）および（６）のｎにｎ＝１を代入する。ｎ＝１では、Ｎｐ＝４５、Ｎｑ＝４７である。したがって、ｎ＝１では、格子点４５の番号「４５」と格子点４７の番号「４７」が選択される。交換手段９３は、選択された格子点４５の番号「４５」と、選択された格子点４７の番号「４７」とを交換する。図１１に、番号が交換された後のｋ空間を概略的に示す。この交換により、格子点４５は境界Ｌの外側の領域（高周波領域ＲＨ）に配置され、一方、格子点４７は境界Ｌの内側の領域（低周波領域ＲＬ）に配置される。
【００４５】
次に、交換手段９３は、式（５）および（６）のｎにｎ＝２を代入する。ｎ＝２では、Ｎｐ＝４３、Ｎｑ＝４９である。したがって、ｎ＝２では、格子点４３の番号「４３」と格子点４９の番号「４９」が選択される。交換手段９３は、選択された格子点４３の番号「４３」と、選択された格子点４９の番号「４９」とを交換する。図１２に、番号が交換された後のｋ空間を概略的に示す。この交換により、格子点４３は境界Ｌの外側の領域（高周波領域ＲＨ）に配置され、一方、格子点４９は境界Ｌの内側の領域（低周波領域ＲＬ）に配置される。
【００４６】
次に、交換手段９３は、式（５）および（６）のｎにｎ＝３を代入する。ｎ＝３では、Ｎｐ＝４１、Ｎｑ＝５１である。したがって、ｎ＝３では、格子点４１の番号「４１」と格子点５１の「番号「５１」が選択される。交換手段９３は、選択された格子点４１の番号「４１」と、選択された格子点５１の番号「５１」とを交換する。図１３に、番号が交換された後のｋ空間を概略的に示す。この交換により、格子点４１は境界Ｌの外側の領域（高周波領域ＲＨ）に配置され、一方、格子点５１は境界Ｌの内側の領域（低周波領域ＲＬ）に配置される。
【００４７】
以下同様に、交換手段９３は、ｎの値がｎ＝１６になるまで、ｎの値をインクリメントする。ｎ＝１６では、Ｎｐ＝１５、Ｎｑ＝７７である。したがって、交換手段９３は、格子点１５の番号「１５」と格子点７７の番号「７７」とを交換する。このようにして、格子点Ｎｐと格子点Ｎｑの番号が交換される。図１４に、１≦ｎ≦１６の間で格子点の番号を交換したときの様子を概略的に示す。
【００４８】
このように、ステップＳＴ１３では、低周波領域ＲＬに含まれる複数の格子点のうち、境界Ｌ側に位置する格子点に付与された番号と、高周波領域ＲＨに含まれる複数の格子点のうち、境界Ｌ側に位置する格子点に付与された番号とが交換される。したがって、格子点１〜格子点４５のうち番号が交換されない格子点（例えば、格子点１）は低周波領域ＲＬに残るが、格子点１〜格子点４５のうち番号が交換された格子点（例えば、格子点１５）は境界Ｌに対して高周波領域ＲＨ側に配置される。また、格子点４６〜格子点２２５のうち番号が交換されない格子点（例えば、格子点２２５）は高周波領域ＲＬに残るが、格子点４６〜格子点２２５のうち番号が交換された格子点（例えば、格子点４７）は境界Ｌに対して低周波領域ＲＬ側に配置される。
【００４９】
第１の形態では、図１４に示す各格子点に付与されている番号「１」〜「２２５」が、各格子点のデータが収集される順番を表している。このようにして、各格子点のデータの収集順序を決定することができる。また、図１４において、格子点１〜格子点４５は、後述するステップＳＴ２で実行されるスキャンのデータ収集期間Ａ１、Ａ２、およびＡ３（図１５参照）でデータ収集される格子点を表しており、一方、格子点４６〜格子点２２５は、後述するステップＳＴ２で実行されるスキャンのデータ収集期間Ｂ（図１５参照）でデータ収集される格子点を表している。したがって、ステップＳＴ１を実行することにより、ｋ空間に含まれる複数の格子点の中から、データ収集期間Ａ１、Ａ２、およびＡ３（図１５参照）でデータ収集される格子点１〜格子点４５と、データ収集期間Ｂ（図１５参照）でデータ収集される格子点４６〜格子点２２５を決定することができる。図１４に示す格子点１〜格子点２２５を決定したら、ステップＳＴ１のフローを終了する。ステップＳＴ１が終了した後、ステップＳＴ２（図７参照）に進む。
【００５０】
ステップＳＴ２では、格子点１〜格子点２５５のデータの収集するためのスキャンを実行する。
【００５１】
図１５は、ステップＳＴ２において実行されるスキャンＳＣの説明図である。
スキャンＳＣでは、先ず、データ収集期間Ａ１において、脂肪抑制パルスＦが印加される。脂肪抑制パルスＦが印加された直後に、シーケンスＶが実行される。データ収集期間Ａ１では、格子点に付与された番号順に、格子点１〜格子点４５のデータＤＡ１を収集する。図１５では、図面の制約上、データ収集期間Ａ１において収集された格子点１〜格子点４５のデータを「●」（黒丸）で示してある。格子点１〜格子点４５のデータを収集した後、次のデータ収集期間Ｂに移行する。
【００５２】
データ収集期間Ｂでも、データ収集期間Ａ１と同様に、脂肪抑制パルスＦが印加される。脂肪抑制パルスＦが印加された直後に、シーケンスＶが実行される。データ収集期間Ｂでは、格子点に付与された番号順に、格子点４６〜格子点２２５のデータを収集する。図１５では、図面の制約上、データ収集期間Ｂにおいて収集された格子点４６〜格子点２２５のデータを「○」（白丸）で示してある。
【００５３】
データ収集期間Ａ１で得られたデータＤＡ１と、データ収集期間Ｂで得られたデータＤＢにより、格子点１〜格子点２２５のデータを含むｋ空間データＤ１を得ることができる。
【００５４】
図１６は、ｋ空間データＤ１の拡大図である。
データ収集期間Ａ１では、以下の格子点のデータが収集される。
（ａ１）ｋ空間の中心の格子点１のデータ、および格子点１の近傍の格子点（例えば、格子点６）のデータ
（ａ２）低周波領域ＲＬのうちの高周波領域ＲＨ側に位置する複数の格子点のうちの一部の格子点（例えば、格子点３２、格子点３６）のデータ
（ａ３）高周波領域ＲＨのうちの低周波領域ＲＬ側に位置する複数の格子点のうちの一部の格子点（例えば、格子点３５、格子点３９）のデータ
【００５５】
一方、データ収集期間Ｂでは、以下の格子点のデータが収集される。
（ｂ１）高周波領域ＲＨの外縁に近い格子点（例えば、格子点２２５）のデータ
（ｂ２）高周波領域ＲＨのうちの低周波領域ＲＬ側に位置する複数の格子点のうちの一部の格子点（例えば、格子点４８、格子点６４）のデータ
（ｂ３）低周波領域ＲＬのうちの高周波領域ＲＨ側に位置する複数の格子点のうちの一部の格子点（例えば、格子点５１、格子点６３）のデータ
【００５６】
したがって、データ収集期間Ａ１では、低周波領域ＲＬの格子点のデータだけでなく、高周波領域ＲＨの格子点のデータも収集される。また、データ収集期間Ｂでは、高周波領域ＲＨの格子点のデータだけでなく、低周波領域ＲＬの格子点のデータも収集される。
【００５７】
図１７〜図１９は、ｋ空間データＤ１を３つの部分に分けて示した図である。
図１７はｋ空間データＤ１の中央部分Ｒ１を示す図、図１８はｋ空間データＤ１の外縁部分Ｒ２を示す図、図１９は境界Ｌの近傍部分Ｒ３を示す図である。以下、図１７〜図１９について説明する。
【００５８】
ｋ空間データＤ１の中央部分Ｒ１（図１７参照）では、データ収集期間Ａ１で収集されたデータのみが配置されており、ｋ空間データＤ１の外縁部分Ｒ２（図１８参照）では、データ収集期間Ｂで収集されたデータのみが配置されている。
【００５９】
しかし、境界Ｌの近傍部分Ｒ３（図１９参照）では、データ収集期間Ａ１で収集されたデータと、データ収集期間Ｂで収集されたデータとの両方が配置されている。例えば、境界Ｌの内側（低周波領域ＲＬ側）では、データ収集期間Ａ１で収集された格子点（例えば、格子点３２、格子点３６）のデータや、データ収集期間Ｂで収集された格子点（例えば、格子点５１、格子点６３）のデータが配置されている。また、境界Ｌの外側（高周波領域ＲＨ側）では、データ収集期間Ａ１で収集された格子点（例えば、格子点３５、格子点３９）のデータや、データ収集期間Ｂで収集された格子点（例えば、格子点４８、格子点６４）のデータが配置されている。したがって、境界Ｌの近傍部分Ｒ３では、データ収集期間Ａ１で収集された格子点のデータと、データ収集期間Ａ１で収集された格子点のデータとが、ｋｙ方向およびｋｚ方向に分散するように配置されていることがわかる。
図１５に戻って説明を続ける。
画像作成手段９４（図２参照）は、ｋ空間データＤ１をフーリエ変換する。これにより、画像ＩＭ１を作成することができる。
【００６０】
データ収集期間ＢにおいてデータＤＢを収集した後、データ収集期間Ａ２およびＡ３で更にデータが収集される。
【００６１】
データ収集期間Ａ２でも、脂肪抑制パルスＦが印加される。脂肪抑制パルスＦが印加された直後に、シーケンスＶが実行される。データ収集期間Ａ２では、格子点に付与された番号順に、格子点１〜格子点４５のデータＤＡ２を収集する。したがって、データ収集期間Ｂで得られたデータＤＢと、データ収集期間Ａ２で得られたデータＤＡ２により、格子点１〜格子点２２５のデータを含む２番目のｋ空間データＤ２を得ることができる。
【００６２】
ｋ空間データＤ２を得た後、画像作成手段９４は、ｋ空間データＤ２をフーリエ変換する。これにより、画像ＩＭ２を作成することができる。データ収集期間Ａ２においてデータＤＡ２を収集した後、データ収集期間Ａ３に移行する。
【００６３】
データ収集期間Ａ３でも、脂肪抑制パルスＦが印加される。脂肪抑制パルスＦが印加された直後に、シーケンスＶが実行される。データ収集期間Ａ３では、格子点に付与された番号順に、格子点１〜格子点４５のデータＤＡ３を収集する。したがって、データ収集期間Ｂで得られたデータＤＢと、データ収集期間Ａ３で得られたデータＤＡ３により、格子点１〜格子点２２５のデータを含む３番目のｋ空間データＤ３を得ることができる。
【００６４】
ｋ空間データＤ３を得た後、画像作成手段９４は、ｋ空間データＤ３をフーリエ変換する。これにより、画像ＩＭ３を作成することができる。
このようにしてスキャンＳＣが実行され、図７のフローが終了する。
【００６５】
第１の形態では、ステップＳＴ１３において、低周波領域に位置する一部の格子点の番号と、高周波領域に位置する一部の格子点の番号が交換される。したがって、図１４に示すように、格子点１〜格子点４５のうちの一部の格子点（例えば、格子点３５、格子点３９）は、境界Ｌに対して高周波領域ＲＨ側に位置する。また、格子点４６〜格子点２２５のうち一部の格子点（例えば、格子点５１、格子点６３）は境界Ｌに対して低周波領域ＲＬ側に位置する。そして、ステップＳＴ２において、データ収集期間Ａ１（Ａ２およびＡ３）では格子点１〜格子点４５のデータを収集し、データ収集期Ｂでは格子点４６〜格子点２２５のデータを収集する。したがって、境界Ｌの近傍部分Ｒ３（図１９参照）では、データ収集期間Ａ１（Ａ２およびＡ３）で収集された格子点のデータと、データ収集期間Ｂで収集された格子点のデータが、ｋｙ方向およびｋｚ方向に分散するように配置される。このため、境界Ｌの近傍部分Ｒ３において、信号強度の大きいデータと信号強度の小さいデータが分散するように配置されるので、信号強度の差が原因で生じるアーチファクト（例えば、ゴースト）を低減することができる。
【００６６】
尚、第１の形態では、データ収集期間Ａ１、Ｂ、Ａ２、およびＡ３では１個の脂肪抑制パルスＦが印加されている（図１５参照）。しかし、データ収集期間Ｂは、データ収集期間Ａ１、Ａ２、およびＡ３よりも収集するデータ数が多い。したがって、データ収集期間Ｂでは、データ収集期間Ｂの終了に近づくにつれて、脂肪抑制パルスＦの脂肪抑制効果が低下してしまい、データ収集期間Ａ１、Ａ２、およびＡ３と同等の脂肪抑制効果が得られないことが考えられる。そこで、データ収集期間Ｂでは、データ収集期間Ａ１、Ａ２、およびＡ３よりも収集するデータ数が多いことを考慮して、複数の脂肪抑制パルスを印加することが望ましい（図２０参照）。
【００６７】
図２０は、データ収集期間Ｂに複数の脂肪抑制パルスを印加した例を示す図である。
図２０には、データ収集期間Ｂに、４つの脂肪抑制パルスＦ１〜Ｆ４が印加される例が示されている。脂肪抑制パルスＦ１とＦ２との間に、格子点４６〜格子点９０のデータが収集される。脂肪抑制パルスＦ２とＦ３との間に、格子点９１〜格子点１３５のデータが収集される。脂肪抑制パルスＦ３とＦ４との間に、格子点１３６〜格子点１８０のデータが収集される。脂肪抑制パルスＦ４の後には、格子点１８１〜格子点２２５のデータが収集される。
【００６８】
データ収集期間Ｂでは、脂肪抑制パルスを印加した後に４５個の格子点のデータを収集するサイクルが４回繰り返されている。したがって、データ収集期間Ｂでも、データ収集期間Ａ１、Ａ２、およびＡ３と同様の脂肪抑制効果を得ることができる。
【００６９】
尚、第１の形態では、式（１）および（２）の係数αおよび係数βを、α＝β＝２に設定している。したがって、格子点Ｎｐの番号は２番づつ小さくなるように選択され、格子点Ｎｑの番号は２番づつ大きくなるように選択されている。しかし、α＝β＝２に限定されることはなく、別の値に設定してもよい。例えば、α＝β＝３に設定した場合は、格子点Ｎｐの番号は３番づつ小さくなるように選択され、格子点Ｎｑの番号は３番づつ大きくなるように選択されるので、交換される格子点の間隔を広げることができる。また、αおよびβは、必ずしも同じ値に設定する必要はなく、異なっていてもよい。更に、αおよびβを、格子点と境界との距離に応じて変更してもよい。例えば、格子点ＮｐおよびＮｑと境界Ｌとの距離が近いうちは、α＝β＝２で格子点Ｎｐの番号と格子点Ｎｑの番号との交換を行い、格子点ＮｐおよびＮｑと境界Ｌとの距離が遠くなってきたら、α＝β＝３に変更し格子点の番号の交換を行ってもよい。αおよびβを、格子点と境界との距離に応じて変更することにより、格子点のデータ収集順序に適した番号の交換を行うことができる。
【００７０】
第１の形態では、１つのデータ収集期間Ｂで格子点４５〜格子点２５５のデータを収集している。しかし、データ収集期間Ｂを複数のデータ収集期間に分け、格子点４５〜格子点２２５のデータを複数のデータ収集期間に分けて収集してもよい。以下に、格子点４５〜格子点２２５のデータを複数のデータ収集期間に分けて収集する一例について説明する。
【００７１】
図２１は、格子点４５〜格子点２２５のデータを複数のデータ収集期間に分けて収集する一例の説明図である。図２１では、データ収集期間Ｂを、２つのデータ収集期間Ｂ１およびＢ２に分けた例が示されている。
【００７２】
データ収集期間Ａ１、Ａ２、およびＡ３の各々は、格子点１〜格子点４５のデータを収集するための期間である。一方、データ収集期間Ｂ１およびＢ２は、格子点４５〜格子点２２５のデータを分割して収集するために設けられている。具体的には、データ収集期間Ｂ１は格子点４６〜格子点１３５のデータを収集するための期間であり、データ収集期間Ｂ２は格子点１３６〜格子点２２５のデータを収集するための期間である。
【００７３】
図２１では、データ収集期間Ａ１で収集されたデータＤＡ１と、データ収集期間Ｂ１で収集されたデータＤＢ１と、データ収集期間Ｂ２で収集されたデータＤＢ２により、１番目のｋ空間データＤ１が得られる。ｋ空間データＤ１をフーリエ変換することにより、画像ＩＭ１を作成することができる。
【００７４】
また、データ収集期間Ｂ１で収集されたデータＤＢ１と、データ収集期間Ａ２で収集されたデータＤＡ２と、データ収集期間Ｂ２で収集されたデータＤＢ２により、２番目のｋ空間データＤ２が得られる。ｋ空間データＤ２をフーリエ変換することにより、画像ＩＭ２を作成することができる。
【００７５】
そして、データ収集期間Ｂ１で収集されたデータＤＢ１と、データ収集期間Ｂ２で収集されたデータＤＢ２と、データ収集期間Ａ３で収集されたデータＤＡ３により、３番目のｋ空間データＤ３が得られる。ｋ空間データＤ３をフーリエ変換することにより、画像ＩＭ３を作成することができる。
【００７６】
尚、図２１では、格子点４６〜格子点１３５のデータをデータ収集期間Ｂ１で収集し、格子点１３６〜格子点２２５のデータをデータ収集期間Ｂ２で収集している。しかし、格子点４６〜格子点２２５のデータをデータ収集期間Ｂ１およびＢ２に分けて収集する方法は、図２１に限定されることはない。例えば、格子点４６〜格子点２２５のうちの奇数番号の格子点のデータをデータ収集期間Ｂ１で収集し、格子点４６〜格子点２２５のうちの偶数番号の格子点のデータをデータ収集期間Ｂ２で収集してもよい。
【００７７】
（２）第２の形態
第２の形態のＭＲ装置は、第１の形態のＭＲ装置と比較すると、プロセッサ９が実行する処理に相違点があるが、その他の点については第１の形態のＭＲ装置と同じである。したがって、第２の形態のＭＲ装置については、主に、プロセッサ９について説明する。
【００７８】
図２２は、第２の形態におけるプロセッサ９が実行する処理の説明図である。
メモリ１０には、プロセッサ９により実行されるプログラムが記憶されている。プロセッサ９は、メモリ１０に記憶されているプログラムを読み出し、プログラムに記述されている処理を実行する。プロセッサ９は、メモリ１０に記憶されているプログラムを読み出すことにより、ナンバリング手段９１〜ソート手段９５などを構成する。
【００７９】
ナンバリング手段９１〜画像作成手段９４は第１の形態と同じであるので説明は省略する。
ソート手段９５は、格子点の番号を交換した後に、格子点の番号を並び替える手段である。
第２の形態では、プロセッサ９は、ナンバリング手段９１〜ソート手段９５を構成する一例である。尚、ナンバリング手段９１、設定手段９２、交換手段９３、およびソート手段９５を合わせたものが、決定手段に相当する。
【００８０】
以下に、第２の形態におけるダイナミックＭＲイメージングのフローについて、図７を参照しながら説明する。
【００８１】
ステップＳＴ１では、格子点のデータ収集順序を決定する（図２３参照）。
図２３は、第２の形態におけるステップＳＴ１のフローを示す図である。
ステップＳＴ１１〜ステップＳＴ１３は、第１の形態と同じであるので、説明は省略する。ステップＳＴ１３において、番号を交換した後（図１４参照）、ステップＳＴ１４に進む。
【００８２】
ステップＳＴ１４では、ソート手段９５（図２２参照）が、図１４に示す格子点１〜格子点４５の番号を並べ替える。具体的には、格子点がｋ空間の中心（格子点１）から離れるに従って番号が大きくなるように、図１４に示す格子点１〜格子点４５の番号を並べ替える。図２４に、ソート手段９５による番号の並替えが終了した後の格子点１〜格子点４５の番号を示す。
【００８３】
第２の形態では、番号が並び替えられた後のｋ空間の各格子点に付与されている番号「１」〜「２２５」（図２４参照）が、各格子点のデータが収集される順番を表している。尚、第２の形態でも、第１の形態と同様に、格子点１〜格子点４５は、データ収集期間Ａ１、Ａ２、およびＡ３でデータ収集される格子点を表しており、一方、格子点４６〜格子点２２５は、データ収集期間Ｂでデータ収集される格子点を表している。
【００８４】
ステップＳＴ１４で番号を並び替えたら、ステップＳＴ１を終了する。ステップＳＴ１が終了した後、ステップＳＴ２（図７参照）に進む。
【００８５】
ステップＳＴ２ではスキャンが実行される。スキャンのデータ収集期間Ａ１、Ａ２、およびＡ３では、図２４に示す各格子点に付与された番号順に、格子点１〜格子点４５のデータが収集される。一方、データ収集期間Ｂでは、図２４に示す各格子点に付与された番号順に、格子点４６〜格子点２２５のデータを収集する。データを収集したら、フローを終了する。
【００８６】
第２の形態では、ソート手段９５により、格子点１〜格子点４５の番号は、格子点１に近い順に並べ替えられる。したがって、データ収集期間Ａ１、Ａ２、およびＡ３では、脂肪抑制パルスを印加した後、ｋ空間の中心に近い格子点のデータから先に収集されるので、脂肪抑制パルスの脂肪抑制効果を更に高めることができる。尚、図２４では、格子点１〜格子点４５の番号のみを並べ替えた例について示されているが、格子点４６〜格子点２２５の番号も、格子点がｋ空間の中心（格子点１）から離れるに従って番号が大きくなるように並べ替えてもよい。
【００８７】
第１および第２の形態では、データ収集期間Ａ１、Ａ２、およびＡ３において、１個の脂肪抑制パルスＦが印加されている。しかし、データ収集期間Ａ１、Ａ２、およびＡ３において、複数の脂肪抑制パルスが印加されるようにしてもよい。図２５は、データ収集期間Ａ１に複数の脂肪抑制パルスを印加した例を示す図である。
【００８８】
図２５には、データ収集期間Ａ１に２つの脂肪抑制パルスＦ１およびＦ２が印加される例が示されている。脂肪抑制パルスＦ１とＦ２との間には、格子点１〜格子点２２のデータが収集され、脂肪抑制パルスＦ２の後には、格子点２３〜格子４５のデータが収集される。
【００８９】
図２５では、脂肪抑制パルスＦ１の後に、２つめの脂肪抑制パルスＦ２が印加されている。したがって、データ収集期間Ａ１の途中で脂肪抑制パルスＦ１の脂肪抑制効果が低下しても、２つの脂肪抑制パルスＦ２により、データ収集期間Ａ１が終了するまで、十分な脂肪抑制効果を維持することができる。また、図２６に示すように、データ収集期間Ａ１の前半で格子点１〜格子点４５のうちの奇数番目の格子点のデータを先に収集し、データ収集期間Ａ１の後半で格子点１〜格子点４５のうちの偶数番目の格子点のデータを後に収集してもよい。
【００９０】
尚、第１および第２の形態では、ステップＳＴ１２において、低周波領域と高周波領域との境界Ｌを設定し、境界Ｌの付近で番号の交換を行っている。しかし、本発明は、低周波領域と高周波領域との境界Ｌを設定する例に限定されることはなく、境界Ｌとは別の境界を設定し、当該別の境界の付近で番号の交換を行ってもよい。更に、境界Ｌを設定せずに、格子点１からの距離に基づいて、図１４に示す格子点１〜格子点２２５を設定してもよい。
【００９１】
第１および第２の形態では、脂肪抑制パルスを用いた造影スキャンについて示されている。しかし、脂肪抑制パルスを用いない造影スキャンであっても、撮影部位に流入する造影剤の量に応じて、撮影部位から収集されるＭＲ信号の信号強度は変化するので、信号強度の差によるアーチファクトが現れることがある。したがって、本発明は、脂肪抑制パルスを用いない造影スキャンを実行する場合にも適用することができる。
【００９２】
また、第１および第２の形態では、造影剤を用いたスキャンについて説明しているが、本発明は、信号強度の差によるアーチファクトが現れるスキャンであれば、造影剤を使用するか否かに関わらず適用することができる。更に、本発明は、上記のスキャン以外にも、信号強度の差によるアーチファクトが現れる他のスキャンに適用することができる。他のスキャンの例としては、例えば、被検体の呼吸情報を含むナビゲータエコーを収集するためのナビゲータシーケンスを用いたスキャンや、血流信号を抑制させるための血流抑制シーケンス（例えば、MSDE（Motion Sensitized Driven Equilibrium）シーケンス）を用いたスキャンなどが考えられる。
【００９３】
第１および第２の形態では、シーケンスＶとして、３Ｄグラディエントエコー系のシーケンスＶを使用した例が示されている。しかし、本発明は、３Ｄグラディエントエコー系のシーケンスに限定されることはなく、３Ｄグラディエントエコー系のシーケンスの代わりに、３Ｄスピンエコー系のシーケンスを用いてもよい。また、３Ｄシーケンスの代わりに、各スライス面のデータを収集する２Ｄシーケンスを用いてもよい。
【００９４】
第１および第２の形態では、１回のシーケンスＶで１つの格子点に配置されるデータを収集している。しかし、ＥＰＩ（Echo Planar Imaging）法などを用いて、１回のシーケンスＶで複数の格子点のデータを収集してもよい。また、本発明は、位相エンコード数を少なくすることにより高速撮像を行うパラレルイメージング法でデータ収集を行う場合にも適用することができる。
【符号の説明】
【００９５】
２  マグネット
３  テーブル
３ａ  クレードル
４  受信コイル
５  送信器
６  勾配磁場電源
７  受信器
８  コンピュータ
９  プロセッサ
１０  メモリ
１１  操作部
１２  表示部
１３  被検体
２１  ボア
９１  ナンバリング手段
９２  設定手段
９３  交換手段
９４  画像作成手段
１００  ＭＲ装置
１０１  造影剤投与装置
【図１】
0006362950.tif 000002
【図２】
0006362950.tif 000003
【図３】
0006362950.tif 000004
【図４】
0006362950.tif 000005
【図５】
0006362950.tif 000006
【図６】
0006362950.tif 000007
【図７】
0006362950.tif 000008
【図８】
0006362950.tif 000009
【図９】
0006362950.tif 000010
【図１０】
0006362950.tif 000011
【図１１】
0006362950.tif 000012
【図１２】
0006362950.tif 000013
【図１３】
0006362950.tif 000014
【図１４】
0006362950.tif 000015
【図１５】
0006362950.tif 000016
【図１６】
0006362950.tif 000017
【図１７】
0006362950.tif 000018
【図１８】
0006362950.tif 000019
【図１９】
0006362950.tif 000020
【図２０】
0006362950.tif 000021
【図２１】
0006362950.tif 000022
【図２２】
0006362950.tif 000023
【図２３】
0006362950.tif 000024
【図２４】
0006362950.tif 000025
【図２５】
0006362950.tif 000026
【図２６】
0006362950.tif 000027