(57)【要約】      （修正有）【課題】ゴーストなどのアーチファクトを低減する磁気共鳴装置を提供する。【解決手段】撮影部位から複数のエコートレインを収集するためのスキャンを実行するスキャン手段と、各エコートレインに含まれる第１エコーを、撮影部位における所定方向の信号値の変化を表すプロファイルに変換する変換手段と、変換手段によりに得られた複数のプロファイルの中から基準プロファイルとして使用されるプロファイルを設定する設定手段と、基準プロファイルと複数のプロファイルの各々との相関係数に基づいて、再収集するエコートレインを決定する決定手段とを有する磁気共鳴装置。【選択図】図７
2015085002.tif 000002
【特許請求の範囲】
【請求項１】
撮影部位から複数のエコートレインを収集するためのスキャンを実行するスキャン手段と、
各エコートレインに含まれる第ｉエコーを、前記撮影部位における所定方向の信号値の変化を表すプロファイルに変換する変換手段と、
前記変換手段によりに得られた複数のプロファイルの中から基準プロファイルとして使用されるプロファイルを設定する設定手段と、
前記基準プロファイルと前記複数のプロファイルの各々との相関係数に基づいて、再収集するエコートレインを決定する決定手段と、
を有する磁気共鳴装置。
【請求項２】
前記決定手段は、
前記相関係数の度数分布を表す度数分布データを作成し、前記度数分布データに基づいて、再収集するエコートレインを決定する、請求項１に記載の磁気共鳴装置。
【請求項３】
前記度数分布データは、相関係数の度数が最大となる第１のビンを有しており、
前記決定手段は、
相関係数が前記第１のビンに含まれていない場合、エコートレインを再収集すると決定する、請求項２に記載の磁気共鳴装置。
【請求項４】
前記複数のプロファイルの中から、前記基準プロファイルの候補となる２つ以上のプロファイルを選択する候補選択手段を有し、
前記設定手段は、前記基準プロファイルの候補ごとに算出された相関係数に基づいて、前記基準プロファイルの候補の中から前記基準プロファイルとして使用されるプロファイルを設定する、請求項１に記載の磁気共鳴装置。
【請求項５】
前記設定手段は、
前記基準プロファイルの候補ごとに、前記相関係数の度数分布を表す度数分布データを作成する作成手段と、
前記度数分布データに基づいて、前記基準プロファイルの候補の中から前記基準プロファイルとして使用されるプロファイルを選択するプロファイル選択手段と、
を有する、請求項４に記載の磁気共鳴装置。
【請求項６】
前記決定手段は、
前記作成手段により作成された複数の度数分布データの中から、再収集するエコートレインを決定するときに使用する第１の度数分布データを特定する、請求項５に記載の磁気共鳴装置。
【請求項７】
前記第１の度数分布データは、相関係数の度数が最大となる第１のビンを有しており、
前記決定手段は、
相関係数が前記第１のビンに含まれていない場合、エコートレインを再収集すると決定する、請求項６に記載の磁気共鳴装置。
【請求項８】
前記変換手段は、
再収集されたエコートレインに含まれる前記第ｉのエコーを、前記撮影部位における所定方向の信号値の変化を表す第１のプロファイルに変換し、
前記スキャン手段は、
前記基準プロファイルと前記第１のプロファイルとの相関係数が前記第１のビンに含まれていない場合、相関係数が前記第１のビンに含まれるまで、エコートレインを再収集するためのシーケンスを実行する、請求項３又は７に記載の磁気共鳴装置。
【請求項９】
前記スキャンは、
ｋ空間の第１の領域に配置されるエコーを含むエコートレインを収集するためのシーケンスが複数回実行される第１の期間と、
前記第１の期間で収集された複数のエコートレインの一部のエコートレインを再収集するためのシーケンスが実行される第２の期間と、
ｋ空間の第２の領域に配置されるエコーを含むエコートレインを収集するためのシーケンスが複数回実行される第３の期間と、
を有し、
前記決定手段は、前記第１の期間のシーケンスが実行されたときに、前記度数分布データを作成する、請求項３又は７に記載の磁気共鳴装置。
【請求項１０】
前記変換手段は、
前記第１の期間で収集されたエコートレインに含まれる前記第ｉのエコーを、前記撮影部位における所定方向の信号値の変化を表す第１のプロファイルに変換し、
前記スキャン手段は、
前記基準プロファイルと前記第１のプロファイルとの相関係数が前記第１のビンに含まれていない場合、相関係数が前記第１のビンに含まれるまで、エコートレインを再収集するためのシーケンスを実行する、請求項９に記載の磁気共鳴装置。
【請求項１１】
前記変換手段は、
前記第３の期間で収集されたエコートレインに含まれる前記第ｉのエコーを、前記撮影部位における所定方向の信号値の変化を表す第２のプロファイルに変換し、
前記スキャン手段は、
前記基準プロファイルと前記第２のプロファイルとの相関係数が前記第１のビンに含まれていない場合、相関係数が前記第１のビンに含まれるまで、エコートレインを再収集するためのシーケンスを実行する、請求項９又は１０に記載の磁気共鳴装置。
【請求項１２】
前記第ｉエコーは第１エコーである、請求項１〜１１のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１３】
撮影部位から複数のエコートレインを収集するためのスキャンを実行する磁気共鳴装置に適用されるプログラムであって、
各エコートレインに含まれる第ｉエコーを、前記撮影部位における所定方向の信号値の変化を表すプロファイルに変換する変換処理と、
前記変換処理によりに得られた複数のプロファイルの中から基準プロファイルとして使用されるプロファイルを設定する設定処理と、
前記基準プロファイルと前記複数のプロファイルの各々との相関係数に基づいて、再収集するエコートレインを決定する決定処理と、
を計算機に実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エコートレインを収集するためのシーケンスを複数回実行する磁気共鳴装置、およびこの磁気共鳴装置に適用可能なプログラムに関する。
【背景技術】
【０００２】
拡散の情報を取得する方法として、スピンの動きを検出するためのＭＰＧ（Motion Probing Gradient）を用いた方法が知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１２−１５７６８７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ＭＰＧを用いて拡散の情報を得る方法の一例として、ＭＰＧをプリパルスとして印加したdiffusion preparationと呼ばれる手法がある。この手法は、例えば、ＦＳＥ（Fast Spin Echo）を収集するパルスシーケンスを実行する場合に適用されており、比較的容易に拡散の情報を得ることができる。
【０００５】
しかし、上記の手法では、心拍動などの体動の影響を受けて、ゴーストなどのアーチファクトが現れることがある。アーチファクトは画質を劣化させる原因であるので、アーチファクトを低減することができる技術が望まれている。
【課題を解決するための手段】
【０００６】
本発明の第１の観点は、撮影部位から複数のエコートレインを収集するためのスキャンを実行するスキャン手段と、
各エコートレインに含まれる第ｉエコーを、前記撮影部位における所定方向の信号値の変化を表すプロファイルに変換する変換手段と、
前記変換手段によりに得られた複数のプロファイルの中から基準プロファイルとして使用されるプロファイルを設定する設定手段と、
前記基準プロファイルと前記複数のプロファイルの各々との相関係数に基づいて、再収集するエコートレインを決定する決定手段と、
を有する磁気共鳴装置である。
【０００７】
本発明の第２の観点は、撮影部位から複数のエコートレインを収集するためのスキャンを実行する磁気共鳴装置に適用されるプログラムであって、
各エコートレインに含まれる第ｉエコーを、前記撮影部位における所定方向の信号値の変化を表すプロファイルに変換する変換処理と、
前記変換処理によりに得られた複数のプロファイルの中から基準プロファイルとして使用されるプロファイルを設定する設定処理と、
前記基準プロファイルと前記複数のプロファイルの各々との相関係数に基づいて、再収集するエコートレインを決定する決定処理と、
を計算機に実行させるためのプログラムである。
【発明の効果】
【０００８】
本発明では、基準プロファイルと各プロファイルとの相関係数に基づいて、体動の影響を大きく受けているエコートレインを特定することができる。したがって、体動の影響を大きく受けているエコートレインを再収集することができ、アーチファクトが低減された画像を得ることができる。
【図面の簡単な説明】
【０００９】
【図１】本発明の第１の形態の磁気共鳴装置の概略図である。
【図２】第１の形態における撮影部位を概略的に示す図である。
【図３】撮影部位Ｒからデータを収集するために実行されるスキャンＳＣの一例を示す図である。
【図４】シーケンスＳＥ１〜ＳＥ２１を実行することにより収集されたエコートレインを概略的に示す図である。
【図５】シーケンスＳＥ１〜ＳＥ２１により収集されたエコーをｋ空間に埋めるときの様子を示す図である。
【図６】第１の形態で実行されるスキャンＳＣの説明図である。
【図７】第１の形態におけるスキャンＳＣを実行するときのフローを示す図である。
【図８】シーケンスＳＥ１が実行されるときの様子を示す図である。
【図９】シーケンスＳＥ２により収集されたエコートレインＥＴ２を概略的に示す図である。
【図１０】シーケンスＳＥ２１が実行されたときの様子を示す図である。
【図１１】基準プロファイルＦｒｅｆを示す図である。
【図１２】算出された相関係数Ｃを示す図である。
【図１３】第１エコーＥ１１〜Ｅ２１，１と相関係数Ｃとの対応関係を表す図である。
【図１４】相関係数の度数分布データＤ０を概略的に示す図である。
【図１５】エコートレインＥＴ１を再収集するための１回目のシーケンスＳＥ１を実行するときの説明図である。
【図１６】期間α２においてシーケンスＳＥ１が３回実行されたときの様子を示す図である。
【図１７】シーケンスＳＥ２を実行したときの様子を示す図である。
【図１８】第２の形態のＭＲ装置の概略図である。
【図１９】第２の形態のフローの説明図である。
【図２０】基準プロファイルＦｒｅｆの候補として選択されたプロファイルを示す図である。
【図２１】算出された相関係数Ｃを示す図である。
【図２２】相関係数の度数分布データＤ１を概略的に示す図である。
【図２３】算出された相関係数Ｃを示す図である。
【図２４】相関係数の度数分布データＤ２を概略的に示す図である。
【図２５】ノイズの影響が小さい方のプロファイルを選択するときの説明図である。
【図２６】第３の形態におけるスキャンＳＣの説明図である。
【図２７】第３の形態におけるスキャンＳＣを実行するときのフローを示す図である。
【図２８】シーケンスＳＥ１が実行されるときの様子を示す図である。
【図２９】シーケンスＳＥ２により収集されたエコートレインＥＴ２を概略的に示す図である。
【図３０】シーケンスＳＥ１０が実行されたときの様子を示す図である。
【図３１】基準プロファイルＦｒｅｆを示す図である。
【図３２】算出された相関係数Ｃを示す図である。
【図３３】第１エコーＥ１１〜Ｅ１０１と相関係数Ｃとの対応関係を示す図である。
【図３４】相関係数の度数分布データＤ３を概略的に示す図である。
【図３５】エコートレインＥＴ１を再収集するための１回目のシーケンスＳＥ１を実行するときの説明図である。
【図３６】期間β２においてシーケンスＳＥ１が３回実行されたときの様子を示す図である。
【図３７】シーケンスＳＥ２を実行したときの様子を示す図である。
【図３８】シーケンスＳＥ１１を実行するときの説明図である。
【図３９】期間β３においてシーケンスＳＥ１１が３回実行されたときの様子を示す図である。
【図４０】シーケンスＳＥ１２〜ＳＥ２１を実行したときの様子を概略的に示す図である。
【発明を実施するための形態】
【００１０】
以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。
【００１１】
（１）第１の形態
図１は、本発明の第１の形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ。ＭＲ：Magnetic Resonance）１００は、マグネット２、テーブル３、受信コイル４などを有している。
【００１２】
マグネット２は、被検体１１が収容されるボア２１を有している。また、マグネット２は、超伝導コイル、勾配コイル、およびＲＦコイルなどが内蔵されている。
【００１３】
テーブル３は、被検体１１を支持するクレードル３ａを有している。クレードル３ａは、ボア２１内に移動できるように構成されている。クレードル３ａによって、被検体１１はボア２１に搬送される。
【００１４】
受信コイル４は、被検体１１の頭部に取り付けられている。受信コイル４は、被検体１１からの磁気共鳴信号を受信する。
【００１５】
ＭＲ装置１００は、更に、送信器５、勾配磁場電源６、受信器７、制御部８、操作部９、および表示部１０などを有している。
【００１６】
送信器５はＲＦコイルに電流を供給し、勾配磁場電源６は勾配コイルに電流を供給する。受信器７は、受信コイル４から受け取った信号に対して、検波などの信号処理を実行する。尚、マグネット２、受信コイル４、送信器５、勾配磁場電源６、受信器７を合わせたものが、スキャン手段に相当する。
【００１７】
制御部８は、表示部１０に必要な情報を伝送したり、受信コイル４から受け取ったデータに基づいて画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。制御部８は、変換手段８１〜決定手段８４などを有している。
【００１８】
変換手段８１は、第１エコーを逆フーリエ変換し、後述するプロファイルを求める。
設定手段８２は、基準プロファイルとして使用されるプロファイルを設定する。
算出手段８３は、基準プロファイルと各プロファイルとの相関係数を算出する。
決定手段８４は、相関係数に基づいて、再収集するエコートレインを決定する。
【００１９】
制御部８は、変換手段８１〜決定手段８４を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。
【００２０】
操作部９は、オペレータにより操作され、種々の情報を制御部８に入力する。表示部１０は種々の情報を表示する。
ＭＲ装置１００は、上記のように構成されている。
【００２１】
図２は第１の形態における撮影部位を概略的に示す図である。第１の形態では、頭部が撮影部位Ｒとして設定されている。
【００２２】
図３は、撮影部位Ｒからデータを収集するために実行されるスキャンＳＣの一例を示す図である。
【００２３】
図３では、シーケンスＳＥｉ（ｉ＝１〜２１）が繰り返し実行される例が示されている。シーケンスＳＥｉは、プリパレーション部Ａとデータ収集部Ｂとを有している。
【００２４】
プリパレーション部Ａは、拡散強調を行うための傾斜磁場ＭＰＧを有している。データ収集部Ｂは、撮影部位ＲからエコートレインＥＴｉを収集するためのシーケンス部である。
【００２５】
エコートレインＥＴｉはｎ個のエコーＥｉ１〜Ｅｉｎを有している。第１の形態では、エコーＥｉ１〜Ｅｉｎのうち第２エコーＥｉ２〜第ｎエコーＥｉｎは画像再構成のデータとして使用される。しかし、第１エコーＥｉ１は、画像再構成のデータとしては使用されず、エコートレインを再収集するか否かを判断するために使用される。第１エコーＥｉ１を用いてエコートレインを再収集するか否かを判断する方法については後述する。
【００２６】
画像再構成のデータとして使用される第２エコーＥｉ２〜第ｎエコーＥｉｎを収集する場合、傾斜磁場ＧｓｓおよびＧｐｅのエンコード量は増減する。しかし、画像再構成のデータとして使用されない第１エコーＥｉ１を収集する場合、傾斜磁場ＧｓｓおよびＧｐｅのエンコード量はゼロに設定されている。
【００２７】
シーケンスＳＥ１〜ＳＥ２１を実行することにより、シーケンスごとに、エコートレインが収集される。図４に、シーケンスＳＥ１〜ＳＥ２１を実行することにより収集されたエコートレインを概略的に示す。図４では、シーケンスＳＥ１、ＳＥ２、ＳＥ３、およびＳＥ２１により収集されるエコートレインＥＴ１、ＥＴ２、ＥＴ３、およびＥＴ２１が示されている。
【００２８】
図５は、シーケンスＳＥ１〜ＳＥ２１により収集されたエコーをｋ空間に埋めるときの様子を示す図である。図５には、ｋｙ−ｋｚ面が示されている。
シーケンスＳＥ１を実行することによりエコートレインＥＴ１（第１エコーＥ１１〜第ｎエコーＥ１ｎ）が収集される。上記のように、第１エコーＥ１１は画像再構成のデータとしては使用されず、第２エコーＥ１２〜第ｎエコーＥ１ｎが画像再構成のデータとして使用される。エコートレインＥＴ１の第２エコーＥ１２〜第ｎエコーＥ１ｎは、ｋｚ＝−１０のｋ空間のデータとして使用される。シーケンスＳＥ１を実行した後、次のシーケンスＳＥ２が実行される。
【００２９】
シーケンスＳＥ２を実行することによりエコートレインＥＴ２（第１エコーＥ２１〜第ｎエコーＥ２ｎ）が収集される。第１エコーＥ２１は画像再構成のデータとしては使用されず、第２エコーＥ２２〜第ｎエコーＥ２ｎが画像再構成のデータとして使用される。エコートレインＥＴ２の第２エコーＥ２２〜第ｎエコーＥ２ｎは、ｋｚ＝−９のｋ空間のデータとして使用される。
【００３０】
以下同様に、各シーケンスが実行され、第２エコー〜第ｎエコーが画像再構成のデータとして使用される。最後にシーケンスＳＥ２１が実行される。
【００３１】
シーケンスＳＥ２１を実行することによりエコートレインＥＴ２１（第１エコーＥ２１，１〜第ｎエコーＥ２１，ｎ）が収集される。第１エコーＥ２１，１は画像再構成のデータとしては使用されず、第２エコーＥ２１，２〜第ｎエコーＥ２１，ｎが画像再構成のデータとして使用される。エコートレインＥＴ２１の第２エコーＥ２１，２〜第ｎエコーＥ２１，ｎは、ｋｚ＝１０のｋ空間のデータとして使用される。
【００３２】
このようにしてｋ空間のデータが収集される。尚、図５では、シーケンスの実行回数は２１回であるが、画像再構成に必要なデータを収集することができるのであれば、シーケンスの実行回数は２１回に限定されることはなく、２１回より少なくてもよいし、２１回より多くてもよい。
【００３３】
ｋ空間のデータを実空間のデータに変換することにより、撮影部位Ｒの拡散の情報を含む画像を得ることができる。しかし、心拍などの体動はエコートレインのエコーに影響を与えることがあり、体動の影響を受けているエコーを画像再構成のエコーとして使用すると、ゴーストなどのアーチファクトの原因となる。そこで、第１の形態では、アーチファクトの原因となる可能性の高いエコートレインが含まれている場合、そのエコートレインの再収集が行われる。以下に、第１の形態のスキャンについて説明する。
【００３４】
図６は、第１の形態で実行されるスキャンＳＣの説明図である。
第１の形態では、スキャンＳＣは、２つの期間α１およびα２に分けられる。期間α１は、シーケンスＳＥ１〜ＳＥ２１を実行するための期間である。一方、期間α２は、エコートレインを再収集するためのシーケンスが実行される期間である。以下、第１の形態のスキャンＳＣを実行するときのフローについて説明する。
【００３５】
図７は、第１の形態におけるスキャンＳＣを実行するときのフローを示す図である。
ステップＳＴ１では、スキャンＳＣの期間α１におけるシーケンスが実行される（図８参照）。
【００３６】
図８〜図１０は、スキャンＳＣの期間α１において実行されるシーケンスの説明図である。
期間α１では、先ず、シーケンスＳＥ１が実行される。図８に、シーケンスＳＥ１が実行されるときの様子を示す。シーケンスＳＥ１が実行されることにより、エコートレインＥＴ１が収集される。エコートレインＥＴ１の第２エコーＥ１２〜第ｎエコーＥ１ｎは、ｋｚ＝−１０に配置される。また、変換手段８１（図１参照）は、画像再構成に使用されない第１エコーＥ１１をｘ方向（周波数エンコード方向）に逆フーリエ変換し、これにより、プロファイルＦ１が得られる。プロファイルＦ１は、ｘ方向（周波数エンコード方向）の信号値の変化を表している。
シーケンスＳＥ１を実行した後、次のシーケンスＳＥ２が実行される。
【００３７】
シーケンスＳＥ２が実行されることにより、エコートレインＥＴ２が収集される。図９に、シーケンスＳＥ２により収集されたエコートレインＥＴ２を概略的に示す。エコートレインＥＴ２の第２エコーＥ２２〜第ｎエコーＥ２ｎは、ｋｚ＝−９に配置される。また、変換手段８１は、画像再構成に使用されない第１エコーＥ２１をｘ方向（周波数エンコード方向）に逆フーリエ変換し、これにより、プロファイルＦ２が得られる。プロファイルＦ２は、ｘ方向（周波数エンコード方向）の信号値の変化を表している。
【００３８】
以下同様に、各シーケンスを順に実行し、各シーケンスで得られた第１エコーを逆フーリエ変換し、プロファイルを求める。図１０に、最後のシーケンスＳＥ２１が実行されたときの様子を示す。エコートレインＥＴ２１の第２エコーＥ２１，２〜第ｎエコーＥ２１，ｎは、ｋｚ＝１０に配置される。また、変換手段８１は、画像再構成に使用されない第１エコーＥ２１，１をｘ方向（周波数エンコード方向）に逆フーリエ変換し、プロファイルＦ２１を求める。プロファイルＦ２１は、ｘ方向（周波数エンコード方向）の信号値の変化を表している。
【００３９】
このようにして、シーケンスＳＥ１〜シーケンスＳＥ２１が実行される。シーケンスＳＥ１〜シーケンスＳＥ２１を実行した後、ステップＳＴ２に進む。
【００４０】
ステップＳＴ２では、エコートレインＥＴ１〜ＥＴ２１の中から、再収集する必要があるエコートレインを特定する。
【００４１】
心拍などの体動はエコートレインのエコーに影響を与えることがあり、体動の影響を受けているエコーを画像再構成のエコーとして使用すると、ゴーストなどのアーチファクトの原因となる。そこで、第１の形態では、エコートレインＥＴ１〜ＥＴ２１の中に、アーチファクトの原因となる可能性の高いエコートレインが含まれている場合、そのエコートレインを特定する。以下に、この特定方法について説明する。
【００４２】
先ず、ステップＳＴ２１において、設定手段８２（図１参照）が、プロファイルＦ１〜Ｆ２１の中から基準プロファイルＦｒｅｆとして使用されるプロファイルを設定する。基準プロファイルＦｒｅｆは、次のステップＳＴ２２においてプロファイルの相関係数を求めるときの基準となるプロファイルである。ここでは、プロファイルＦ１を基準プロファイルＦｒｅｆとして設定する。図１１に、基準プロファイルＦｒｅｆを示す。
基準プロファイルＦｒｅｆを設定した後、ステップＳＴ２２に進む。
【００４３】
ステップＳＴ２２では、算出手段８３（図１参照）が、基準プロファイルＦｒｅｆと、プロファイルＦ１〜Ｆ２１の各々との相関係数Ｃを求める。図１２に、算出された相関係数Ｃを示す。図１２では、図面のスペースの制約上、代表して、以下の３つの相関係数が示されている。
（１）基準プロファイルＦｒｅｆとプロファイルＦ１との相関係数Ｃ
（２）基準プロファイルＦｒｅｆとプロファイルＦ２との相関係数Ｃ
（３）基準プロファイルＦｒｅｆとプロファイルＦ２１との相関係数Ｃ
【００４４】
本形態では、プロファイルＦ１が基準プロファイルＦｒｅｆとして使用されているので、基準プロファイルＦｒｅｆとプロファイルＦ１との相関係数Ｃは、Ｃ＝１となる。また、第１の形態では、基準プロファイルＦｒｅｆとプロファイルＦ２との相関係数Ｃは、Ｃ＝０．９７であり、基準プロファイルＦｒｅｆとプロファイルＦ２１との相関係数Ｃは、Ｃ＝０．９２であるとする。
【００４５】
上記のように相関係数Ｃを求めることにより、第１エコーＥ１１〜Ｅ２１，１と相関係数Ｃとの対応関係を知ることができる。図１３に、第１エコーＥ１１〜Ｅ２１，１と相関係数Ｃとの対応関係を表す散布図が示されている。散布図の横軸は第１エコーＥ１１〜Ｅ２１，１を表しており、縦軸は相関係数Ｃを表す。散布図では、相関係数ＣはＣ＝０．７〜１．０の範囲が示されている。第１のエコーＥ１１〜Ｅ２１，１の各々における相関係数Ｃの値は、記号「○」で示されている。相関係数Ｃを求めた後、ステップＳＴ２３に進む。
【００４６】
ステップＳＴ２３では、決定手段８４（図１参照）が、相関係数Ｃに基づいて、再収集するエコートレインを決定する。以下に、相関係数Ｃに基づいて、再収集するエコートレインを決定する方法について説明する。
先ず、決定手段８４は、相関係数の度数分布を表す度数分布データを作成する（図１４参照）。
【００４７】
図１４は相関係数の度数分布データＤ０を概略的に示す図である。
決定手段８４は、先ず、散布図の相関係数の範囲０．７〜１．０を複数の区間（ビン）に分ける。図１４では、散布図の相関係数の範囲０．７〜１．０を６個の区間、即ち、ビン１〜ビン６に分けたときの様子が示されている。
【００４８】
散布図の相関係数の範囲０．７〜１．０をビン１〜ビン６に分けた後、決定手段８４は、各ビンに含まれている相関係数の個数（度数）を求める。したがって、相関係数の度数分布データＤ０を作成することができる。
【００４９】
度数分布データＤ０を作成した後、決定手段８４は、ビン１〜ビン６の中から、相関係数の度数が最大となるビンを特定する。図１４では、ビン１の度数は２個、ビン２の度数は１９個、ビン３〜ビン６の度数は０個である。したがって、ビン２が、相関係数の度数が最大となるビンとして特定される。
【００５０】
一般的に、被検体から収集した複数のエコーが体動の影響をあまり受けていない場合、これらの複数のエコーの各々のプロファイルを比較すると、プロファイルは類似する傾向がある。つまり、被検体から収集した複数のエコーが体動の影響をあまり受けていない場合、相関係数のばらつきは小さくなる傾向がある。一方、被検体から取集した複数のエコーが体動の影響を大きく受けている場合、これらの複数のエコーのプロファイルはあまり類似せず、相関係数のばらつきは大きくなる傾向がある。したがって、第１のエコーが体動の影響をあまり受けていない場合、相関係数Ｃは同一のビンに集中して現れる。図１４の度数分布データＤ０を参照すると、ビン２の度数が最大であるので、相関係数はビン２に集中して現れていることが分かる。ビン２には、第１エコーＥ３１〜Ｅ２１，１の相関係数が含まれている。したがって、第１エコーＥ３１〜Ｅ２１，１は体動の影響をそれほど受けていないと考えられる。
【００５１】
一方、ビン２とは別のビンに含まれている第１エコーは、体動の影響を大きく受けていると考えられる。第１エコーが体動の影響を大きく受けている場合、その直後に収集される第２エコー〜第ｎエコーも体動の影響を大きく受けていると考えられる。したがって、これらの第２エコー〜第ｎエコーを画像再構成のエコーとして使用すると、ゴーストなどのアーチファクトが発生する原因となる。そこで、第１エコーの相関係数がビン２に含まれていない場合、エコートレインを再収集すると決定する。図１４を参照すると、第１エコーＥ１１およびＥ２１がビン２に含まれていないことが分かる。したがって、決定手段８４は、第１エコーＥ１１を含むエコートレインＥＴ１と、第１エコーＥ２１を含むエコートレインＥＴ２とを再収集すると決定する。
再収集するエコートレインＥＴ１およびＥＴ２を決定した後、ステップＳＴ３に進む。
【００５２】
ステップＳＴ３では、エコートレインＥＴ１およびＥＴ２の再収集が行われる。第１の形態では、スキャンＳＣの期間α１の直後に、エコートレインＥＴ１およびＥＴ２を再収集するための期間α２が設けられる（図１５〜図１７参照）。
【００５３】
図１５〜図１７は、期間α２にエコートレインＥＴ１およびＥＴ２を再収集するときの説明図である。
先ず、エコートレインＥＴ１を再収集するための１回目のシーケンスＳＥ１が実行される（図１５参照）。
【００５４】
図１５は、エコートレインＥＴ１を再収集するための１回目のシーケンスＳＥ１を実行するときの説明図である。
変換手段８１は、１回目のシーケンスＳＥ１により収集された第１エコーＥ１１を逆フーリエ変換し、プロファイルＦ１を求める。そして、算出手段８３は基準プロファイルＦｒｅｆとプロファイルＦ１との相関係数Ｃを計算する。ここでは、Ｃ＝０．８２であるとする。次に、決定手段８４は、相関係数０．８２がビン２（相関係数０．９０〜０．９５）に含まれているか否かを判断する。図１５を参照すると、相関係数０．８２はビン２には含まれていないことが分かる。上記のように、相関係数がビン２に含まれていない場合、第１エコーは体動の影響を大きく受けていると考えられる。したがって、１回目のシーケンスＳＥ１により再収集された第２エコーＥ１２〜第ｎエコーＥ１ｎは、ｋ空間のデータとして採用することができないと判断される。そこで、期間α２では、相関係数がビン２に含まれるまで、エコートレインＥＴ１を再収集するためのシーケンスＳＥ１が繰り返し実行される（図１６参照）。
【００５５】
図１６は、期間α２においてシーケンスＳＥ１が３回実行されたときの様子を示す図である。
【００５６】
図１６では、１回目および２回目のシーケンスＳＥ１で再収集された第１エコーＥ１１の相関係数はビン２には含まれていないが、３回目のシーケンスＳＥ１で収集された第１エコーＥ１１の相関係数Ｃがビン２に含まれたとする。したがって、３回目のシーケンスＳＥ１で再収集されたエコートレインＥＴ１の第２エコーＥ１２〜第ｎエコーＥ１ｎが、ｋ空間のデータとして採用される。
【００５７】
ｋ空間のデータとして採用されるエコートレインＥＴ１を再収集することができたら、エコートレインＥＴ２を再収集するためのシーケンスＳＥ２を実行する（図１７参照）。
【００５８】
図１７は、シーケンスＳＥ２を実行したときの様子を示す図である。
シーケンスＳＥ２は、第１エコーＥ２１の相関係数がビン２に含まれるまで繰り返し実行される。シーケンスＳＥ２については、２回目の再収集で第１エコーＥ２１の相関係数がビン２に含まれたとする。したがって、２回目の再収集で得られたエコートレインＥＴ２の第２エコーＥ２２〜第ｎエコーＥ２ｎが、ｋ空間のデータとして採用される。
【００５９】
このようにして、体動の影響をあまり受けていないエコートレインＥＴ１およびＥＴ２を再収集することができる。尚、エコートレインＥＴ３〜ＥＴ２１は、ステップＳＴ１において、体動の影響をあまり受けていないエコートレインとして収集されているので、エコートレインＥＴ３〜ＥＴ２１の再収集は行われない。
【００６０】
したがって、エコートレインＥＴ１およびＥＴ２が再収集されたｋ空間のデータＤ′を得ることができる。エコートレインＥＴ１およびＥＴ２を再収集した後、ステップＳＴ４に進む。
【００６１】
ステップＳＴ４では、図１７に示すｋ空間のデータＤ′を用いて画像再構成が行われ、フローが終了する。
【００６２】
第１の形態では、ステップＳＴ１でエコートレインＥＴ１〜ＥＴ２１を収集し、エコートレインＥＴ１〜ＥＴ２１の第１エコーＥ１１〜Ｅ２１，１を逆フーリエ変換し、プロファイルＦ１〜Ｆ２１を求める（図１０参照）。そして、基準プロファイルＦｒｅｆを設定し（図１１参照）、基準プロファイルＦｒｅｆとプロファイルＦ１〜Ｆ２１の各々との相関係数Ｃを計算する（図１２参照）。第１のエコーが体動の影響をあまり受けていない場合は相関係数Ｃのばらつきは小さいが、第１のエコーが体動の影響を大きく受けている場合は相関係数Ｃのばらつきは大きくなる傾向がある。したがって、度数が最も大きいビンを特定することにより、体動の影響をあまり受けていないエコートレインと、体動の影響を大きく受けているエコートレインとを区別することができる。本形態では、体動の影響を大きく受けているエコートレインはステップＳＴ３において再収集されるので、ゴーストなどのアーチファクトが低減された画像を得ることができる。
【００６３】
尚、第１の形態では、プロファイルＦ１を基準プロファイルＦｒｅｆとして使用している。しかし、プロファイルＦ１とは別のプロファイルを基準プロファイルＦｒｅｆとして使用してもよい。
【００６４】
また、第１の形態では、相関係数Ｃの範囲の下限値を０．７に設定している（例えば、図１７参照）。しかし、ステップＳＴ２２で算出した相関係数に、０．７を下回る値が含まれている場合は、その値が含まれるように、相関係数の範囲の下限値を設定すればよい。例えば、ステップＳＴ２２で算出した相関係数の最小値が、Ｃ＝０．６５の場合、相関係数Ｃの範囲の下限値は０．６５以下の値（例えば、０．６）に設定すればよい。
【００６５】
（２）第２の形態
第１の形態では、プロファイルＦ１が基準プロファイルＦｒｅｆとして使用されている。しかし、プロファイルＦ１がノイズの影響を強く受けてしまっている場合、本来は再収集しなければならないエコートレインが、再収集する必要がないと判断されてしまい、ゴーストなどのアーチファクトを十分に低減することができないことがある。したがって、基準プロファイルＦｒｅｆは、ノイズの影響をあまり受けていないことが望ましい。そこで、第２の形態では、ノイズの影響をあまり受けていないプロファイルを基準プロファイルとして設定する例について説明する。
【００６６】
図１８は、第２の形態のＭＲ装置の概略図である。
第２の形態のＭＲ装置２００は、第１の形態のＭＲ装置１００（図１参照）と比較すると、制御部８が異なっているが、その他の構成は同じである。したがって、第２の形態のＭＲ装置２００については、主に、制御部８について説明する。制御部８は、変換手段８０１〜決定手段８０５などを有している。
【００６７】
変換手段８０１は、第１エコーを逆フーリエ変換し、後述するプロファイルを求める。
候補選択手段８０２は、複数のプロファイルの中から、基準プロファイルの候補となる２つ以上のプロファイルを選択する。
算出手段８０３は、基準プロファイルの候補となるプロファイルと各プロファイルとの相関係数を算出する。
設定手段８０４は、相関係数に基づいて、基準プロファイルの候補の中から、基準プロファイルとして使用されるプロファイルを設定する。設定手段８０４は、度数分布データ作成手段８０４ａと、プファイル選択手段８０４ｂとを有している。度数分布データ作成手段８０４ａは、基準プロファイルの候補ごとに、相関係数の度数分布を表す度数分布データを作成する。プロファイル選択手段８０４ｂは、基準プロファイルの候補の中から、基準プロファイルとして使用されるプロファイルを選択する。
決定手段８０５は、相関係数に基づいて、再収集するエコートレインを決定する。
【００６８】
制御部８は、変換手段８０１〜決定手段８０５を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。
【００６９】
第２の形態では、制御部８は上記のように構成されている。次に、第２の形態においてスキャンを実行するときのフローについて説明する。
【００７０】
図１９は、第２の形態のフローの説明図である。
ステップＳＴ１は、第１の形態と同様に、シーケンスＳＥ１〜シーケンスＳＥ２１が実行される。したがって、第１の形態において図８〜図１０を参照しながら説明したように、シーケンスＳＥ１〜シーケンスＳＥ２１が実行されることにより、エコートレインＥＴ１〜ＥＴ２１が取得される。変換手段８０１（図１８参照）は、第１エコーＥ１１〜Ｅ２１，１を逆フーリエ変換する。これにより、プロファイルＦ１〜Ｆ２１が得られる。シーケンスＳＥ１〜シーケンスＳＥ２１を実行した後、ステップＳＴ２に進む。
【００７１】
ステップＳＴ２では、エコートレインＥＴ１〜ＥＴ２１の中から、再収集する必要があるエコートレインを特定する。ステップＳＴ２はステップＳＴ２０１〜ＳＴ２０３を有している。以下、各ステップＳＴ２０１〜ＳＴ２０３について説明する。
【００７２】
ステップＳＴ２０１では、候補選択手段８０２（図１８参照）が、プロファイルＦ１〜Ｆ２０の中から基準プロファイルＦｒｅｆの候補となるプロファイルを選択する。図２０に、基準プロファイルＦｒｅｆの候補として選択されたプロファイルを示す。第２の形態では、プロファイルＦ１およびＦ２が、基準プロファイルＦｒｅｆの候補として選択されるが、３つ以上のプロファイルを、基準プロファイルＦｒｅｆの候補として選択してもよい。基準プロファイルＦｒｅｆの候補Ｆ１およびＦ２を選択した後、ステップＳＴ２０２に進む。
【００７３】
ステップＳＴ２０２では、基準プロファイルＦｒｅｆの候補として選択されたプロファイルＦ１およびＦ２のうち、どちらのプロファイルがノイズの影響が小さいかを判断する。ステップＳＴ２０２は、ステップＳＴ２０２ａ〜ＳＴ２０２ｅを有しているので、各ステップＳＴ２０２ａ〜ＳＴ２０２ｅについて説明する。
【００７４】
ステップＳＴ２０２ａでは、算出手段８０３（図１８参照）が、基準プロファイルＦｒｅｆの候補として選択されたプロファイルＦ１と、プロファイルＦ１〜Ｆ２１の各々との相関係数を求める。図２１に、算出された相関係数Ｃを示す。また、図２１の左下には、第１のエコーＥ１１〜Ｅ２１，１と相関係数Ｃとの対応関係を表す散布図が示されている。散布図の横軸は第１のエコーＥ１１〜Ｅ２１，１を表しており、縦軸は相関係数Ｃを表す。散布図では、相関係数ＣはＣ＝０．７〜１．０の範囲が示されている。第１のエコーＥ１１〜Ｅ２１，１の各々における相関係数Ｃの値は、記号「○」で示されている。相関係数Ｃを求めた後、ステップＳＴ２０２ｂに進む。
【００７５】
ステップＳＴ２０２ｂでは、度数分布データ作成手段８０４ａ（図１８参照）が、相関係数の度数分布を表す度数分布データを作成する（図２２参照）。
【００７６】
図２２は、相関係数の度数分布データＤ１を概略的に示す図である。
度数分布データ作成手段８０４ａは、先ず、散布図の相関係数の範囲０．７〜１．０を複数の区間（ビン）に分ける。図２２では、散布図の相関係数の範囲０．７〜１．０を６個の区間、即ち、ビン１〜ビン６に分けたときの様子が示されている。
【００７７】
散布図の相関係数の範囲０．７〜１．０をビン１〜ビン６に分けた後、度数分布データ作成手段８０４ａは、各ビンに含まれている相関係数の個数（度数）を求める。このようにして、度数分布データＤ１を作成することができる。図２２では、ビン１の度数は２個、ビン２の度数は１９個、ビン３〜ビン６の度数は０個であるので、ビン２の度数が最大になっている。度数分布データＤ１を作成した後、ステップＳＴ２０２ｃに進む。
【００７８】
ステップＳＴ２０２ｃでは、算出手段８０３が、基準プロファイルＦｒｅｆの候補として選択されたプロファイルＦ２と、プロファイルＦ１〜Ｆ２１の各々との相関係数を求める。図２３に、算出された相関係数Ｃを示す。また、図２３の左下には、第１のエコーＥ１１〜Ｅ２１，１と相関係数Ｃとの対応関係を表す散布図が示されている。散布図の横軸は第１のエコーＥ１１〜Ｅ２１，１を表しており、縦軸は相関係数Ｃを表す。散布図では、相関係数ＣはＣ＝０．７〜１．０の範囲が示されている。第１のエコーＥ１１〜Ｅ２１，１の各々における相関係数Ｃの値は、記号「○」で示されている。相関係数Ｃを求めた後、ステップＳＴ２０２ｄに進む。
【００７９】
ステップＳＴ２０２ｄでは、度数分布データ作成手段８０４ａが、相関係数の度数分布を表す度数分布データを作成する（図２４参照）。
【００８０】
図２４は、相関係数の度数分布データＤ２を概略的に示す図である。
度数分布データ作成手段８０４ａは、先ず、散布図の相関係数の範囲０．７〜１．０を複数の区間（ビン）に分ける。図２４では、散布図の相関係数の範囲０．７〜１．０を６個の区間、即ち、ビン１〜ビン６に分けたときの様子が示されている。
【００８１】
散布図の相関係数の範囲０．７〜１．０をビン１〜ビン６に分けた後、度数分布データ作成手段８０４ａは、各ビンに含まれている相関係数の個数（度数）を求める。このようにして、度数分布データＤ２を作成することができる。図２４では、ビン１の度数は１９個、ビン２の度数は２個、ビン３〜ビン６の度数は０個であるので、ビン１の度数が最大になっている。度数分布データＤ２を作成した後、ステップＳＴ２０２ｅに進む。
【００８２】
ステップＳＴ２０２ｅでは、プロファイル選択手段８０４ｂが、度数分布データＤ１（図２２参照）と度数分布データＤ２（図２４参照）とに基づいて、プロファイルＦ１およびＦ２の中から、ノイズの影響が小さい方のプロファイルを選択する（図２５参照）。
【００８３】
図２５は、ノイズの影響が小さい方のプロファイルを選択するときの説明図である。
図２５には、図２２に示す度数分布データＤ１（基準プロファイルＦｒｅｆの候補：プロファイルＦ１）と、図２４に示す度数分布データＤ２（基準プロファイルＦｒｅｆの候補：プロファイルＦ２）とが示されている。
【００８４】
プロファイル選択手段８０４ｂは、先ず、度数分布データＤ１の中から、度数が最大となるビン２に含まれる相関係数Ｃを抽出する。そして、ビン２に含まれる相関係数の平均値を求める。ここでは、ビン２に含まれる相関係数の平均値を０．９３とする。
【００８５】
次に、プロファイル選択手段８０４ｂは、度数分布データＤ２の中から、度数が最大となるビン１に含まれる相関係数Ｃを抽出する。そして、ビン１に含まれる相関係数の平均値を求める。ここでは、ビン１に含まれる相関係数の平均値を０．９７とする。
【００８６】
相関係数の平均値を求めた後、プロファイル選択手段８０４ｂは、これらの相関係数の平均値を比較する。一般的に、プロファイルがノイズの影響をあまり受けていない場合、相関係数は大きい値（Ｃ＝１に近い値）になる傾向がある。したがって、相関係数の平均値が大きいほど、プロファイルはノイズの影響をあまり受けていないと考えることができる。度数分布データＤ１と度数分布データＤ２とを比較すると、度数分布データＤ２の相関係数の平均値（０．９７）は、度数分布データＤ１の相関係数の平均値（０．９３）よりも大きい。したがって、プロファイルＦ１とプロファイルＦ２とを比較すると、プロファイルＦ２は、プロファイルＦ１よりも、ノイズの影響を受けていないと考えられる。そこで、プロファイル選択手段８０４ｂは、プロファイルＦ２を、ノイズの影響をあまり受けていないプロファイルとして選択する。このようにして選択されたプロファイルＦ２が基準プロファイルＦｒｅｆとして設定される。基準プロファイルＦｒｅｆを設定した後、ステップＳＴ２０３に進む。
【００８７】
ステップＳＴ２０３では、決定手段８０５（図１８参照）が、再収集するエコートレインを決定する。再収集するエコートレインは、以下の手順で決定する。
【００８８】
決定手段８０５は、先ず、２つの度数分布データＤ１およびＤ２のうちの一方の度数分布データを、再収集するエコートレインを決定するときに使用する度数分布データとして選択する。第２の形態では、プロファイルＦ２が基準プロファイルＦｒｅｆとして設定されたので（ステップＳＴ２０２ｅ参照）、決定手段８０５は、度数分布データＤ２を、再収集するエコートレインを決定するときに使用する度数分布データとして選択する。
【００８９】
度数分布データＤ２を選択した後、決定手段８０５は、度数分布データＤ２のビン１〜ビン６の中で、度数が最大となるビンを特定する。図２４を参照すると、度数分布データＤ２では、ビン１の度数が最大となっている。第１の形態で説明したように、第１のエコーが体動の影響をあまり受けていない場合、相関係数Ｃは同一のビンに集中して現れるので、ビン１に含まれる第１エコーは体動の影響をそれほど受けていないと考えられる。
【００９０】
一方、ビン１とは別のビンに含まれている第１エコーは、体動の影響を大きく受けていると考えられる。図２４の度数分布データＤ２を参照すると、第１エコ−Ｅ１１およびＥ２１，１がビン１に含まれていないことが分かる。したがって、決定手段８０５は、第１エコーＥ１１を含むエコートレインＥＴ１と、第１エコーＥ２１，１を含むエコートレインＥＴ２１とを再収集すると決定する。
再収集するエコートレインを決定した後、ステップＳＴ３に進む。
【００９１】
ステップＳＴ３では、相関係数が度数分布データＤ２のビン１に含まれるまで、エコートレインＥＴ１およびＥＴ２１が再収集される。エコートレインを再収集した後、ステップＳＴ４に進み、画像再構成を行いフローを終了する。
【００９２】
第２の形態では、ノイズの影響をあまり受けていないプロファイルを基準プロファイルとして用いることができるので、エコートレインを再収集するか否かの判断を、より的確に行うことができる。
【００９３】
（３）第３の形態
以下に、第３の形態について説明する。尚、第３の形態のＭＲ装置の基本的な構成は、第１の形態のＭＲ装置１００と同じである。したがって、必要に応じて図１を参照しながら、第３の形態について説明する。
【００９４】
図２６は、第３の形態におけるスキャンＳＣの説明図である。
第３の形態では、スキャンＳＣは、３つの期間β１〜β３に分けられる。
【００９５】
期間β１は、シーケンスＳＥ１〜ＳＥ１０が実行される期間である。
期間β２は、エコートレインを再収集するためのシーケンスが実行される期間である。
期間β３は、シーケンスＳＥ１１〜ＳＥ２１が実行され、必要に応じてエコートレインを再収集するためのシーケンスが実行される期間である。以下、第３の形態で実行されるスキャンＳＣについて説明する。
【００９６】
図２７は、第３の形態におけるスキャンＳＣを実行するときのフローを示す図である。
ステップＳＴで１は、スキャンＳＣの期間β１におけるシーケンスが実行される（図２８参照）。
【００９７】
図２８〜図３０は、スキャンＳＣの期間β１において実行されるシーケンスの説明図である。
【００９８】
期間β１では、先ず、シーケンスＳＥ１が実行される。図２８に、シーケンスＳＥ１が実行されるときの様子を示す。シーケンスＳＥ１が実行されることにより、エコートレインＥＴ１が収集される。エコートレインＥＴ１の第２エコーＥ１２〜第ｎエコーＥ１ｎは、ｋｚ＝−１０に配置される。また、変換手段８１（図１参照）は、第１エコーＥ１１をｘ方向（周波数エンコード方向）に逆フーリエ変換する。これにより、プロファイルＦ１が得られる。プロファイルＦ１は、ｘ方向（周波数エンコード方向）の信号値の変化を表している。
シーケンスＳＥ１を実行した後、次のシーケンスＳＥ２が実行される。
【００９９】
シーケンスＳＥ２が実行されることにより、エコートレインＥＴ２が収集される。図２９に、シーケンスＳＥ２により収集されたエコートレインＥＴ２を概略的に示す。エコートレインＥＴ２の第２エコーＥ２２〜第ｎエコーＥ２ｎは、ｋｚ＝−９に配置される。また、第１エコーＥ２１をｘ方向（周波数エンコード方向）に逆フーリエ変換することにより、プロファイルＦ２が得られる。プロファイルＦ２は、ｘ方向（周波数エンコード方向）の信号値の変化を表している。
【０１００】
以下同様に、各シーケンスを順に実行し、各シーケンスで得られた第１エコーを逆フーリエ変換し、プロファイルを求める。図３０に、シーケンスＳＥ１０が実行されたときの様子を示す。エコートレインＥＴ１０の第２エコーＥ１０，２〜第ｎエコーＥ１０，ｎは、ｋｚ＝−１に配置される。また、第１エコーＥ１０，１をｘ方向（周波数エンコード方向）に逆フーリエ変換することにより、プロファイルＦ１０が得られる。プロファイルＦ１０は、ｘ方向（周波数エンコード方向）の信号値の変化を表している。
【０１０１】
このようにして、シーケンスＳＥ１〜シーケンスＳＥ１０が実行される。シーケンスＳＥ１〜シーケンスＳＥ１０を実行した後、ステップＳＴ２に進む。
【０１０２】
ステップＳＴ２では、エコートレインＥＴ１〜ＥＴ１０の中から、再収集する必要があるエコートレインを特定する。以下、ステップＳＴ２に含まれるステップＳＴ２１〜ＳＴ２３について順に説明する。
【０１０３】
先ず、ステップＳＴ２１において、設定手段８２（図１参照）が、プロファイルＦ１〜Ｆ１０の中から、基準プロファイルＦｒｅｆとして使用されるプロファイルを設定する。第３の形態では、第１の形態と同様に、プロファイルＦ１が基準プロファイルＦｒｅｆとして設定されるとする。図３１に、基準プロファイルＦｒｅｆを示す。
【０１０４】
基準プロファイルＦｒｅｆを設定した後、ステップＳＴ２２に進む。
【０１０５】
ステップＳＴ２２では、算出手段８３（図１参照）が、基準プロファイルＦｒｅｆと、プロファイルＦ１〜Ｆ１０の各々との相関係数Ｃを求める。図３２に、算出された相関係数Ｃを示す。図３２では、図面のスペースの制約上、代表して、以下の３つの相関係数が示されている。
（１）基準プロファイルＦｒｅｆとプロファイルＦ１との相関係数Ｃ
（２）基準プロファイルＦｒｅｆとプロファイルＦ２との相関係数Ｃ
（３）基準プロファイルＦｒｅｆとプロファイルＦ１０との相関係数Ｃ
【０１０６】
本形態では、プロファイルＦ１が基準プロファイルＦｒｅｆとして使用されているので、基準プロファイルＦｒｅｆとプロファイルＦ１との相関係数Ｃは、Ｃ＝１となる。また、第３の形態では、基準プロファイルＦｒｅｆとプロファイルＦ２との相関係数Ｃは、Ｃ＝０．９７であり、基準プロファイルＦｒｅｆとプロファイルＦ１０との相関係数Ｃは、Ｃ＝０．９２であるとする。
【０１０７】
上記のように相関係数Ｃを求めることにより、第１エコーＥ１１〜Ｅ１０，１と相関係数Ｃとの対応関係を知ることができる。図３３に、第１エコーＥ１１〜Ｅ１０，１と相関係数Ｃとの対応関係を表す散布図が示されている。散布図の横軸は第１エコーＥ１１〜Ｅ１０，１を表しており、縦軸は相関係数Ｃを表す。散布図では、相関係数ＣはＣ＝０．７〜１．０の範囲が示されている。第１のエコーＥ１１〜Ｅ１０，１の各々における相関係数Ｃの値は、記号「○」で示されている。相関係数Ｃを求めた後、ステップＳＴ２３に進む。
【０１０８】
ステップＳＴ２３では、決定手段８４（図１参照）が、相関係数に基づいて、再収集するエコートレインを決定する。以下に、相関係数に基づいて、再収集するエコートレインを決定する方法について説明する。
【０１０９】
先ず、決定手段８４は、相関係数の度数分布を表す度数分布データを作成する（図３４参照）。
【０１１０】
図３４は相関係数の度数分布データＤ３を概略的に示す図である。
決定手段８４は、先ず、散布図の相関係数の範囲０．７〜１．０を複数の区間（ビン）に分ける。図３４では、散布図の相関係数の範囲０．７〜１．０を６個の区間、即ち、ビン１〜ビン６に分けたときの様子が示されている。
【０１１１】
散布図の相関係数の範囲０．７〜１．０をビン１〜ビン６に分けた後、決定手段８４は、各ビンに含まれている相関係数の個数（度数）を求める。図３４では、ビン１の度数は２個、ビン２の度数は８個、ビン３〜ビン６の度数は０個である。したがって、図３４では、ビン２の度数が最大になっている。
【０１１２】
第１の形態で説明したように、第１のエコーが体動の影響をあまり受けていない場合、相関係数Ｃは同一のビンに集中して現れる。図３４の度数分布データＤ３を参照すると、ビン２の度数が最大であるので、相関係数はビン２に集中して現れていることが分かる。ビン２には、第１エコーＥ３１〜Ｅ１０，１の相関係数が含まれている。したがって、第１エコーＥ３１〜Ｅ１０，１は体動の影響をそれほど受けていないと考えられる。
【０１１３】
一方、第１エコーＥ１１およびＥ２１はビン２に含まれていないので、体動の影響を受けていると考えられる。したがって、決定手段８４は、第１エコーＥ１１を含むエコートレインＥＴ１と、第１エコーＥ２１を含むエコートレインＥＴ２とを再収集すると決定する。
再収集するエコートレインＥＴ１およびＥＴ２を決定した後、ステップＳＴ３に進む。
【０１１４】
ステップＳＴ３では、エコートレインＥＴ１およびＥＴ２の再収集が行われる。第３の形態では、スキャンＳＣの期間β１の直後に、エコートレインＥＴ１およびＥＴ２を再収集するための期間β２が設けられる（図３５〜図３７参照）。
【０１１５】
図３５〜図３７は、期間β２にエコートレインＥＴ１およびＥＴ２を再収集するときの説明図である。
先ず、エコートレインＥＴ１を再収集するための１回目のシーケンスＳＥ１が実行される（図３５参照）。
【０１１６】
図３５は、エコートレインＥＴ１を再収集するための１回目のシーケンスＳＥ１を実行するときの説明図である。
【０１１７】
変換手段８１は、１回目のシーケンスＳＥ１により収集された第１エコーＥ１１を逆フーリエ変換し、プロファイルＦ１を求める。そして、算出手段８３は基準プロファイルＦｒｅｆとプロファイルＦ１との相関係数Ｃを計算する。ここでは、Ｃ＝０．８２であるとする。次に、決定手段８４は、相関係数０．８２がビン２（相関係数０．９０〜０．９５）に含まれているか否かを判断する。図３５を参照すると、相関係数０．８２はビン２には含まれていないことが分かる。上記のように、相関係数がビン２に含まれていない場合、第１エコーは体動の影響を大きく受けていると考えられる。したがって、１回目のシーケンスＳＥ１により再収集された第２エコーＥ１２〜第ｎエコーＥ１ｎは、ｋ空間のデータとして採用することができない。そこで、期間β２では、相関係数がビン２に含まれるまで、エコートレインＥＴ１を再収集するためのシーケンスＳＥ１が繰り返し実行される（図３６参照）。
【０１１８】
図３６は、期間β２においてシーケンスＳＥ１が３回実行されたときの様子を示す図である。
図３６では、１回目および２回目のシーケンスＳＥ１で再収集された第１エコーＥ１１の相関係数はビン２には含まれていないが、３回目のシーケンスＳＥ１で収集された第１エコーＥ１１の相関係数Ｃがビン２に含まれたとする。したがって、３回目のシーケンスＳＥ１で再収集されたエコートレインＥＴ１の第２エコーＥ１２〜第ｎエコーＥ１ｎが、ｋ空間のデータとして採用される。
【０１１９】
ｋ空間のデータとして採用されるエコートレインＥＴ１を再収集することができたら、エコートレインＥＴ２を再収集するためのシーケンスＳＥ２を実行する（図３７参照）。
【０１２０】
図３７は、シーケンスＳＥ２を実行したときの様子を示す図である。
シーケンスＳＥ２は、第１エコーＥ２１の相関係数がビン２に含まれるまで繰り返し実行される。シーケンスＳＥ２については、２回目の再収集で第１エコーＥ２１の相関係数がビン２に含まれたとする。したがって、２回目の再収集で得られたエコートレインＥＴ２の第２エコーＥ２２〜第ｎエコーＥ２ｎが、ｋ空間のデータとして採用される。
【０１２１】
このようにして、体動の影響をあまり受けていないエコートレインＥＴ１およびＥＴ２を再収集することができる。尚、エコートレインＥＴ３〜ＥＴ１０は、ステップＳＴ１において、体動の影響をあまり受けていないエコートレインとして収集されているので、エコートレインＥＴ３〜ＥＴ１０の再収集は行われない。
【０１２２】
エコートレインＥＴ１およびＥＴ２を再収集した後、ステップＳＴ４に進む。
【０１２３】
ステップＳＴ３では、ｋｚ＝０〜１０に配置されるデータを収集するためのシーケンスＳＥ１１〜ＳＥ２１が実行される。第３の形態では、スキャンＳＣの期間β２の直後に、シーケンスＳＥ１１〜ＳＥ２１を実行するための期間β３が設けられる（図３８〜図４０参照）。
【０１２４】
図３８〜図４０は、期間β３で実行されるシーケンスの説明図である。
期間β３では、先ず、ｋｚ＝０のデータを収集するためのシーケンスＳＥ１１が実行される（図３８参照）。
【０１２５】
図３８はシーケンスＳＥ１１を実行するときの説明図である。
シーケンスＳＥ１１が実行されることにより、エコートレインＥＴ１１が収集される。変換手段８１は、シーケンスＳＥ１１により収集された第１エコーＥ１１，１を逆フーリエ変換し、プロファイルＦ１１を求める。そして、算出手段８３は基準プロファイルＦｒｅｆとプロファイルＦ１１との相関係数Ｃを計算する。ここでは、Ｃ＝０．９７であるとする。相関係数０．９７はビン２（相関係数０．９０〜０．９５）に含まれていなので、第１エコーＥ１１，１は体動の影響を大きく受けていると考えられる。したがって、シーケンスＳＥ１１により収集された第２エコーＥ１１，２〜第ｎエコーＥ１１，ｎは、ｋ空間のデータとして採用することができない。そこで、期間β３では、相関係数がビン２に含まれるまで、エコートレインを再収集するためのシーケンスＳＥ１１が実行される（図３９参照）。
【０１２６】
図３９は、期間β３においてシーケンスＳＥ１１が３回実行されたときの様子を示す図である。
図３９では、１回目および２回目のシーケンスＳＥ１１で再収集された第１エコーＥ１１，１の相関係数はビン２には含まれていないが、３回目のシーケンスＳＥ１１で収集された第１エコーＥ１１１の相関係数Ｃがビン２に含まれたとする。したがって、３回目のシーケンスＳＥ１１で再収集されたエコートレインＥＴ１１の第２エコーＥ１１，２〜第ｎエコーＥ１１，ｎが、ｋｚ＝０のデータとして採用される。
【０１２７】
ｋｚ＝０のデータを取得した後、残りのｋｚ＝１〜１０のデータを収集するためのシーケンスＳＥ１２〜ＳＥ２１を実行する。図４０に、シーケンスＳＥ１２〜ＳＥ２１を実行したときの様子を概略的に示す。シーケンスＳＥ１２〜ＳＥ２１を実行する場合も、シーケンスＳＥ１１と同様に、相関係数Ｃがビン２に含まれるまでシーケンスを実行する。したがって、体動の影響をあまり受けていないエコートレインがｋ空間のデータとして採用される。
【０１２８】
上記のようにして、ｋ空間のデータを取得することができる。ｋ空間のデータを取得した後、ステップＳＴ５に進み、画像再構成が行われ、フローが終了する。
【０１２９】
第３の形態では、ｋ空間の一部の領域（ｋｚ＝−９〜−１）のデータが収集された時点で、度数分布データＤ３（図３４）を作成している。このように、ｋ空間の全領域のデータを収集する前に度数分布データＤ３を作成しても、エコートレインが体動の影響を大きく受けているか否かを判断することができる。したがって、体動の影響を大きく受けているエコートレインは再収集されるので、ゴーストなどのアーチファクトが低減された画像を取得することができる。
【０１３０】
尚、第３の形態においても、第２の形態のように、基準プロファイルの候補となる２つ以上のプロファイルを選択し、選択された２つ以上のプロファイルの中から、基準プロファイルとして使用されるプロファイルを設定してもよい。
【０１３１】
また、第３の形態は、第１の形態と比較すると、度数分布データを作成するために必要な相関係数の個数が少なくて済むので、より短時間で度数分布データを作成することができる。
【０１３２】
尚、第１〜第３の形態では、第１エコーを、画像再構成に使用されないエコーとしている。しかし、画像再構成に使用されないエコーは、第１エコーに限定されることはなく、第１エコー〜第ｎエコーのうちのどのエコーであってもよい。
【０１３３】
また、第１〜第３の形態では、第１エコーと相関係数との関係を表す度数分布データを作成し、再収集するエコートレインを決定している。しかし、画像再構成に使用される第２エコー〜第ｎエコーの中から選択された第ｉエコー（２≦ｉ≦ｎ）のプロファイルを求め、第ｉエコーと相関係数との関係を表す度数分布データに基づいて、再収集するエコートレインを決定してもよい。ただし、第２エコー〜第ｎエコーは位相エンコード勾配が印加されたときに収集されるエコーであるので、位相分散などの影響で、エコーのピーク値は小さくなると考えられる。エコーのピーク値が小さい場合、相関係数の誤差が大きくなることが考えられるので、エコーのピーク値はできるだけ大きいことが望ましい。したがって、位相分散などの影響をできるだけ受けないように、位相エンコード勾配が印加されていないときに収集された第１エコーを用いて度数分布データを作成することが望ましい。
【符号の説明】
【０１３４】
２  マグネット
３  テーブル
３ａ  クレードル
４  受信コイル
５  送信器
６  勾配磁場電源
７  制御部
８  操作部
９  表示部
１０  被検体
２１  ボア
７１  ピーク値検出手段
７２、７１３  決定手段
１００、２００  ＭＲ装置
７１１  変換手段
７１２  検出手段
【図１】
2015085002.tif 000003
【図２】
2015085002.tif 000004
【図３】
2015085002.tif 000005
【図４】
2015085002.tif 000006
【図５】
2015085002.tif 000007
【図６】
2015085002.tif 000008
【図７】
2015085002.tif 000009
【図８】
2015085002.tif 000010
【図９】
2015085002.tif 000011
【図１０】
2015085002.tif 000012
【図１１】
2015085002.tif 000013
【図１２】
2015085002.tif 000014
【図１３】
2015085002.tif 000015
【図１４】
2015085002.tif 000016
【図１５】
2015085002.tif 000017
【図１６】
2015085002.tif 000018
【図１７】
2015085002.tif 000019
【図１８】
2015085002.tif 000020
【図１９】
2015085002.tif 000021
【図２０】
2015085002.tif 000022
【図２１】
2015085002.tif 000023
【図２２】
2015085002.tif 000024
【図２３】
2015085002.tif 000025
【図２４】
2015085002.tif 000026
【図２５】
2015085002.tif 000027
【図２６】
2015085002.tif 000028
【図２７】
2015085002.tif 000029
【図２８】
2015085002.tif 000030
【図２９】
2015085002.tif 000031
【図３０】
2015085002.tif 000032
【図３１】
2015085002.tif 000033
【図３２】
2015085002.tif 000034
【図３３】
2015085002.tif 000035
【図３４】
2015085002.tif 000036
【図３５】
2015085002.tif 000037
【図３６】
2015085002.tif 000038
【図３７】
2015085002.tif 000039
【図３８】
2015085002.tif 000040
【図３９】
2015085002.tif 000041
【図４０】
2015085002.tif 000042