特許 6461464 磁気共鳴装置およびプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6461464

(24)【登録日】2019年1月11日

(45)【発行日】2019年1月30日

(54)【発明の名称】磁気共鳴装置およびプログラム

(51)【国際特許分類】

   A61B 5/055 20060101AFI20190121BHJP

【ＦＩ】

   A61B5/055 370

   A61B5/055 376

【請求項の数】12

【全頁数】31

(21)【出願番号】特願2013-227103(P2013-227103)

(22)【出願日】2013年10月31日

(65)【公開番号】特開2015-85002(P2015-85002A)

(43)【公開日】2015年5月7日

【審査請求日】2016年10月21日

(73)【特許権者】

【識別番号】300019238

【氏名又は名称】ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー

(74)【代理人】

【識別番号】100137545

【弁理士】

【氏名又は名称】荒川 聡志

(72)【発明者】

【氏名】岩舘 雄治

【審査官】 松本 隆彦

(56)【参考文献】

【文献】 特開２００３−２２５２２３（ＪＰ，Ａ）

【文献】 特開２０１３−２０８２０８（ＪＰ，Ａ）

【文献】 特開２００８−１４８９１８（ＪＰ，Ａ）

【文献】 米国特許第０６２７５７２０（ＵＳ，Ｂ１）

【文献】 特開平１１−１２８２０２（ＪＰ，Ａ）

【文献】 特開平１１−１１３８７７（ＪＰ，Ａ）

【文献】 特開２０１４−１４７４６６（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ａ６１Ｂ５／０５５

Ｇ０１Ｒ３３／５６

(57)【特許請求の範囲】

【請求項1】

撮影部位から複数のエコートレインを収集するためのスキャンを実行するスキャン手段と、
各エコートレインに含まれる第ｉエコーを、前記撮影部位における所定方向の信号値の変化を表すプロファイルに変換する変換手段と、
前記変換手段によりに得られた複数のプロファイルの中から基準プロファイルとして使用されるプロファイルを設定する設定手段と、
前記基準プロファイルと前記複数のプロファイルの各々との相関係数の度数分布を表す度数分布データを作成し、前記度数分布データに基づいて、再収集するエコートレインを決定する決定手段と、
を有する磁気共鳴装置。

【請求項2】

前記度数分布データは、相関係数の度数が最大となる第１のビンを有しており、
前記決定手段は、
相関係数が前記第１のビンに含まれていない場合、エコートレインを再収集すると決定する、請求項１に記載の磁気共鳴装置。

【請求項3】

撮影部位から複数のエコートレインを収集するためのスキャンを実行するスキャン手段と、
各エコートレインに含まれる第ｉエコーを、前記撮影部位における所定方向の信号値の変化を表すプロファイルに変換する変換手段と、
前記変換手段によりに得られた複数のプロファイルの中から、２つ以上のプロファイルを基準プロファイルの候補として選択する候補選択手段と、
前記候補選択手段により選択された基準プロファイルの各候補に対して、前記基準プロファイルの候補と複数のプロファイルの各々との相関係数の度数分布を表す度数分布データを作成する作成手段と、
前記度数分布データに基づいて、前記候補選択手段により選択された２つ以上の基準プロファイルの候補の中から、１つのプロファイルを、基準プロファイルとして選択するプロファイル選択手段と、
前記基準プロファイルと前記複数のプロファイルの各々との相関係数に基づいて、再収集するエコートレインを決定する決定手段と、
を有する磁気共鳴装置。

【請求項4】

前記決定手段は、
前記作成手段により作成された複数の度数分布データの中から、再収集するエコートレインを決定するときに使用する第１の度数分布データを特定する、請求項３に記載の磁気共鳴装置。

【請求項5】

前記第１の度数分布データは、相関係数の度数が最大となる第１のビンを有しており、
前記決定手段は、
相関係数が前記第１のビンに含まれていない場合、エコートレインを再収集すると決定する、請求項４に記載の磁気共鳴装置。

【請求項6】

前記変換手段は、
再収集されたエコートレインに含まれる前記第ｉのエコーを、前記撮影部位における所定方向の信号値の変化を表す第１のプロファイルに変換し、
前記スキャン手段は、
前記基準プロファイルと前記第１のプロファイルとの相関係数が前記第１のビンに含まれていない場合、相関係数が前記第１のビンに含まれるまで、エコートレインを再収集するためのシーケンスを実行する、請求項２又は５に記載の磁気共鳴装置。

【請求項7】

前記スキャンは、
ｋ空間の第１の領域に配置されるエコーを含むエコートレインを収集するためのシーケンスが複数回実行される第１の期間と、
前記第１の期間で収集された複数のエコートレインの一部のエコートレインを再収集するためのシーケンスが実行される第２の期間と、
ｋ空間の第２の領域に配置されるエコーを含むエコートレインを収集するためのシーケンスが複数回実行される第３の期間と、
を有し、
前記決定手段は、前記第１の期間のシーケンスが実行されたときに、前記度数分布データを作成する、請求項２又は５に記載の磁気共鳴装置。

【請求項8】

前記変換手段は、
前記第１の期間で収集されたエコートレインに含まれる前記第ｉのエコーを、前記撮影部位における所定方向の信号値の変化を表す第１のプロファイルに変換し、
前記スキャン手段は、
前記基準プロファイルと前記第１のプロファイルとの相関係数が前記第１のビンに含まれていない場合、相関係数が前記第１のビンに含まれるまで、エコートレインを再収集するためのシーケンスを実行する、請求項７に記載の磁気共鳴装置。

【請求項9】

前記変換手段は、
前記第３の期間で収集されたエコートレインに含まれる前記第ｉのエコーを、前記撮影部位における所定方向の信号値の変化を表す第２のプロファイルに変換し、
前記スキャン手段は、
前記基準プロファイルと前記第２のプロファイルとの相関係数が前記第１のビンに含まれていない場合、相関係数が前記第１のビンに含まれるまで、エコートレインを再収集するためのシーケンスを実行する、請求項７又は８に記載の磁気共鳴装置。

【請求項10】

前記第ｉエコーは第１エコーである、請求項１〜９のうちのいずれか一項に記載の磁気共鳴装置。

【請求項11】

撮影部位から複数のエコートレインを収集するためのスキャンを実行する磁気共鳴装置に適用されるプログラムであって、
各エコートレインに含まれる第ｉエコーを、前記撮影部位における所定方向の信号値の変化を表すプロファイルに変換する変換処理と、
前記変換処理によりに得られた複数のプロファイルの中から基準プロファイルとして使用されるプロファイルを設定する設定処理と、
前記基準プロファイルと前記複数のプロファイルの各々との相関係数の度数分布を表す度数分布データを作成し、前記度数分布データに基づいて、再収集するエコートレインを決定する決定処理と、
を計算機に実行させるためのプログラム。

【請求項12】

撮影部位から複数のエコートレインを収集するためのスキャンを実行する磁気共鳴装置に適用されるプログラムであって、
各エコートレインに含まれる第ｉエコーを、前記撮影部位における所定方向の信号値の変化を表すプロファイルに変換する変換処理と、
前記変換処理によりに得られた複数のプロファイルの中から、２つ以上のプロファイルを基準プロファイルの候補として選択する候補選択処理と、
前記候補選択処理により選択された基準プロファイルの各候補に対して、前記基準プロファイルの候補と複数のプロファイルの各々との相関係数の度数分布を表す度数分布データを作成する作成処理と、
前記度数分布データに基づいて、前記候補選択処理により選択された２つ以上の基準プロファイルの候補の中から、１つのプロファイルを、基準プロファイルとして選択するプロファイル選択処理と、
前記基準プロファイルと前記複数のプロファイルの各々との相関係数に基づいて、再収集するエコートレインを決定する決定処理と、
を計算機に実行させるためのプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、エコートレインを収集するためのシーケンスを複数回実行する磁気共鳴装置、およびこの磁気共鳴装置に適用可能なプログラムに関する。

【背景技術】

【0002】

拡散の情報を取得する方法として、スピンの動きを検出するためのＭＰＧ（Motion Probing Gradient）を用いた方法が知られている（特許文献１参照）。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特開２０１２−１５７６８７号公報

【発明の概要】

【発明が解決しようとする課題】

【0004】

ＭＰＧを用いて拡散の情報を得る方法の一例として、ＭＰＧをプリパルスとして印加したdiffusion preparationと呼ばれる手法がある。この手法は、例えば、ＦＳＥ（Fast Spin Echo）を収集するパルスシーケンスを実行する場合に適用されており、比較的容易に拡散の情報を得ることができる。

【0005】

しかし、上記の手法では、心拍動などの体動の影響を受けて、ゴーストなどのアーチファクトが現れることがある。アーチファクトは画質を劣化させる原因であるので、アーチファクトを低減することができる技術が望まれている。

【課題を解決するための手段】

【0006】

本発明の第１の観点は、撮影部位から複数のエコートレインを収集するためのスキャンを実行するスキャン手段と、
各エコートレインに含まれる第ｉエコーを、前記撮影部位における所定方向の信号値の変化を表すプロファイルに変換する変換手段と、
前記変換手段によりに得られた複数のプロファイルの中から基準プロファイルとして使用されるプロファイルを設定する設定手段と、
前記基準プロファイルと前記複数のプロファイルの各々との相関係数に基づいて、再収集するエコートレインを決定する決定手段と、
を有する磁気共鳴装置である。

【0007】

本発明の第２の観点は、撮影部位から複数のエコートレインを収集するためのスキャンを実行する磁気共鳴装置に適用されるプログラムであって、
各エコートレインに含まれる第ｉエコーを、前記撮影部位における所定方向の信号値の変化を表すプロファイルに変換する変換処理と、
前記変換処理によりに得られた複数のプロファイルの中から基準プロファイルとして使用されるプロファイルを設定する設定処理と、
前記基準プロファイルと前記複数のプロファイルの各々との相関係数に基づいて、再収集するエコートレインを決定する決定処理と、
を計算機に実行させるためのプログラムである。

【発明の効果】

【0008】

本発明では、基準プロファイルと各プロファイルとの相関係数に基づいて、体動の影響を大きく受けているエコートレインを特定することができる。したがって、体動の影響を大きく受けているエコートレインを再収集することができ、アーチファクトが低減された画像を得ることができる。

【図面の簡単な説明】

【0009】

【図1】本発明の第１の形態の磁気共鳴装置の概略図である。

【図2】第１の形態における撮影部位を概略的に示す図である。

【図3】撮影部位Ｒからデータを収集するために実行されるスキャンＳＣの一例を示す図である。

【図4】シーケンスＳＥ_１〜ＳＥ_２１を実行することにより収集されたエコートレインを概略的に示す図である。

【図5】シーケンスＳＥ_１〜ＳＥ_２１により収集されたエコーをｋ空間に埋めるときの様子を示す図である。

【図6】第１の形態で実行されるスキャンＳＣの説明図である。

【図7】第１の形態におけるスキャンＳＣを実行するときのフローを示す図である。

【図8】シーケンスＳＥ_１が実行されるときの様子を示す図である。

【図9】シーケンスＳＥ_２により収集されたエコートレインＥＴ_２を概略的に示す図である。

【図10】シーケンスＳＥ_２１が実行されたときの様子を示す図である。

【図11】基準プロファイルＦ_ｒｅｆを示す図である。

【図12】算出された相関係数Ｃを示す図である。

【図13】第１エコーＥ_１１〜Ｅ_２１，１と相関係数Ｃとの対応関係を表す図である。

【図14】相関係数の度数分布データＤ０を概略的に示す図である。

【図15】エコートレインＥＴ_１を再収集するための１回目のシーケンスＳＥ_１を実行するときの説明図である。

【図16】期間α２においてシーケンスＳＥ_１が３回実行されたときの様子を示す図である。

【図17】シーケンスＳＥ_２を実行したときの様子を示す図である。

【図18】第２の形態のＭＲ装置の概略図である。

【図19】第２の形態のフローの説明図である。

【図20】基準プロファイルＦ_ｒｅｆの候補として選択されたプロファイルを示す図である。

【図21】算出された相関係数Ｃを示す図である。

【図22】相関係数の度数分布データＤ１を概略的に示す図である。

【図23】算出された相関係数Ｃを示す図である。

【図24】相関係数の度数分布データＤ２を概略的に示す図である。

【図25】ノイズの影響が小さい方のプロファイルを選択するときの説明図である。

【図26】第３の形態におけるスキャンＳＣの説明図である。

【図27】第３の形態におけるスキャンＳＣを実行するときのフローを示す図である。

【図28】シーケンスＳＥ_１が実行されるときの様子を示す図である。

【図29】シーケンスＳＥ_２により収集されたエコートレインＥＴ_２を概略的に示す図である。

【図30】シーケンスＳＥ_１０が実行されたときの様子を示す図である。

【図31】基準プロファイルＦ_ｒｅｆを示す図である。

【図32】算出された相関係数Ｃを示す図である。

【図33】第１エコーＥ_１１〜Ｅ_１０１と相関係数Ｃとの対応関係を示す図である。

【図34】相関係数の度数分布データＤ３を概略的に示す図である。

【図35】エコートレインＥＴ_１を再収集するための１回目のシーケンスＳＥ_１を実行するときの説明図である。

【図36】期間β２においてシーケンスＳＥ_１が３回実行されたときの様子を示す図である。

【図37】シーケンスＳＥ_２を実行したときの様子を示す図である。

【図38】シーケンスＳＥ_１１を実行するときの説明図である。

【図39】期間β３においてシーケンスＳＥ_１１が３回実行されたときの様子を示す図である。

【図40】シーケンスＳＥ_１２〜ＳＥ_２１を実行したときの様子を概略的に示す図である。

【発明を実施するための形態】

【0010】

以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。

【0011】

（１）第１の形態
図１は、本発明の第１の形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ。ＭＲ：Magnetic Resonance）１００は、マグネット２、テーブル３、受信コイル４などを有している。

【0012】

マグネット２は、被検体１１が収容されるボア２１を有している。また、マグネット２は、超伝導コイル、勾配コイル、およびＲＦコイルなどが内蔵されている。

【0013】

テーブル３は、被検体１１を支持するクレードル３ａを有している。クレードル３ａは、ボア２１内に移動できるように構成されている。クレードル３ａによって、被検体１１はボア２１に搬送される。

【0014】

受信コイル４は、被検体１１の頭部に取り付けられている。受信コイル４は、被検体１１からの磁気共鳴信号を受信する。

【0015】

ＭＲ装置１００は、更に、送信器５、勾配磁場電源６、受信器７、制御部８、操作部９、および表示部１０などを有している。

【0016】

送信器５はＲＦコイルに電流を供給し、勾配磁場電源６は勾配コイルに電流を供給する。受信器７は、受信コイル４から受け取った信号に対して、検波などの信号処理を実行する。尚、マグネット２、受信コイル４、送信器５、勾配磁場電源６、受信器７を合わせたものが、スキャン手段に相当する。

【0017】

制御部８は、表示部１０に必要な情報を伝送したり、受信コイル４から受け取ったデータに基づいて画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。制御部８は、変換手段８１〜決定手段８４などを有している。

【0018】

変換手段８１は、第１エコーを逆フーリエ変換し、後述するプロファイルを求める。
設定手段８２は、基準プロファイルとして使用されるプロファイルを設定する。
算出手段８３は、基準プロファイルと各プロファイルとの相関係数を算出する。
決定手段８４は、相関係数に基づいて、再収集するエコートレインを決定する。

【0019】

制御部８は、変換手段８１〜決定手段８４を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。

【0020】

操作部９は、オペレータにより操作され、種々の情報を制御部８に入力する。表示部１０は種々の情報を表示する。
ＭＲ装置１００は、上記のように構成されている。

【0021】

図２は第１の形態における撮影部位を概略的に示す図である。第１の形態では、頭部が撮影部位Ｒとして設定されている。

【0022】

図３は、撮影部位Ｒからデータを収集するために実行されるスキャンＳＣの一例を示す図である。

【0023】

図３では、シーケンスＳＥ_ｉ（ｉ＝１〜２１）が繰り返し実行される例が示されている。シーケンスＳＥ_ｉは、プリパレーション部Ａとデータ収集部Ｂとを有している。

【0024】

プリパレーション部Ａは、拡散強調を行うための傾斜磁場ＭＰＧを有している。データ収集部Ｂは、撮影部位ＲからエコートレインＥＴ_ｉを収集するためのシーケンス部である。

【0025】

エコートレインＥＴ_ｉはｎ個のエコーＥ_ｉ１〜Ｅ_ｉｎを有している。第１の形態では、エコーＥ_ｉ１〜Ｅ_ｉｎのうち第２エコーＥ_ｉ２〜第ｎエコーＥ_ｉｎは画像再構成のデータとして使用される。しかし、第１エコーＥ_ｉ１は、画像再構成のデータとしては使用されず、エコートレインを再収集するか否かを判断するために使用される。第１エコーＥ_ｉ１を用いてエコートレインを再収集するか否かを判断する方法については後述する。

【0026】

画像再構成のデータとして使用される第２エコーＥ_ｉ２〜第ｎエコーＥ_ｉｎを収集する場合、傾斜磁場Ｇ_ｓｓおよびＧ_ｐｅのエンコード量は増減する。しかし、画像再構成のデータとして使用されない第１エコーＥ_ｉ１を収集する場合、傾斜磁場Ｇ_ｓｓおよびＧ_ｐｅのエンコード量はゼロに設定されている。

【0027】

シーケンスＳＥ_１〜ＳＥ_２１を実行することにより、シーケンスごとに、エコートレインが収集される。図４に、シーケンスＳＥ_１〜ＳＥ_２１を実行することにより収集されたエコートレインを概略的に示す。図４では、シーケンスＳＥ_１、ＳＥ_２、ＳＥ_３、およびＳＥ_２１により収集されるエコートレインＥＴ_１、ＥＴ_２、ＥＴ_３、およびＥＴ_２１が示されている。

【0028】

図５は、シーケンスＳＥ_１〜ＳＥ_２１により収集されたエコーをｋ空間に埋めるときの様子を示す図である。図５には、ｋｙ−ｋｚ面が示されている。
シーケンスＳＥ_１を実行することによりエコートレインＥＴ_１（第１エコーＥ_１１〜第ｎエコーＥ_１ｎ）が収集される。上記のように、第１エコーＥ_１１は画像再構成のデータとしては使用されず、第２エコーＥ_１２〜第ｎエコーＥ_１ｎが画像再構成のデータとして使用される。エコートレインＥＴ_１の第２エコーＥ_１２〜第ｎエコーＥ_１ｎは、ｋｚ＝−１０のｋ空間のデータとして使用される。シーケンスＳＥ_１を実行した後、次のシーケンスＳＥ_２が実行される。

【0029】

シーケンスＳＥ_２を実行することによりエコートレインＥＴ_２（第１エコーＥ_２１〜第ｎエコーＥ_２ｎ）が収集される。第１エコーＥ_２１は画像再構成のデータとしては使用されず、第２エコーＥ_２２〜第ｎエコーＥ_２ｎが画像再構成のデータとして使用される。エコートレインＥＴ_２の第２エコーＥ_２２〜第ｎエコーＥ_２ｎは、ｋｚ＝−９のｋ空間のデータとして使用される。

【0030】

以下同様に、各シーケンスが実行され、第２エコー〜第ｎエコーが画像再構成のデータとして使用される。最後にシーケンスＳＥ_２１が実行される。

【0031】

シーケンスＳＥ_２１を実行することによりエコートレインＥＴ_２１（第１エコーＥ_２１，１〜第ｎエコーＥ_２１，ｎ）が収集される。第１エコーＥ_２１，１は画像再構成のデータとしては使用されず、第２エコーＥ_２１，２〜第ｎエコーＥ_２１，ｎが画像再構成のデータとして使用される。エコートレインＥＴ_２１の第２エコーＥ_２１，２〜第ｎエコーＥ_２１，ｎは、ｋｚ＝１０のｋ空間のデータとして使用される。

【0032】

このようにしてｋ空間のデータが収集される。尚、図５では、シーケンスの実行回数は２１回であるが、画像再構成に必要なデータを収集することができるのであれば、シーケンスの実行回数は２１回に限定されることはなく、２１回より少なくてもよいし、２１回より多くてもよい。

【0033】

ｋ空間のデータを実空間のデータに変換することにより、撮影部位Ｒの拡散の情報を含む画像を得ることができる。しかし、心拍などの体動はエコートレインのエコーに影響を与えることがあり、体動の影響を受けているエコーを画像再構成のエコーとして使用すると、ゴーストなどのアーチファクトの原因となる。そこで、第１の形態では、アーチファクトの原因となる可能性の高いエコートレインが含まれている場合、そのエコートレインの再収集が行われる。以下に、第１の形態のスキャンについて説明する。

【0034】

図６は、第１の形態で実行されるスキャンＳＣの説明図である。
第１の形態では、スキャンＳＣは、２つの期間α１およびα２に分けられる。期間α１は、シーケンスＳＥ_１〜ＳＥ_２１を実行するための期間である。一方、期間α２は、エコートレインを再収集するためのシーケンスが実行される期間である。以下、第１の形態のスキャンＳＣを実行するときのフローについて説明する。

【0035】

図７は、第１の形態におけるスキャンＳＣを実行するときのフローを示す図である。
ステップＳＴ１では、スキャンＳＣの期間α１におけるシーケンスが実行される（図８参照）。

【0036】

図８〜図１０は、スキャンＳＣの期間α１において実行されるシーケンスの説明図である。
期間α１では、先ず、シーケンスＳＥ_１が実行される。図８に、シーケンスＳＥ_１が実行されるときの様子を示す。シーケンスＳＥ_１が実行されることにより、エコートレインＥＴ_１が収集される。エコートレインＥＴ_１の第２エコーＥ_１２〜第ｎエコーＥ_１ｎは、ｋｚ＝−１０に配置される。また、変換手段８１（図１参照）は、画像再構成に使用されない第１エコーＥ_１１をｘ方向（周波数エンコード方向）に逆フーリエ変換し、これにより、プロファイルＦ_１が得られる。プロファイルＦ_１は、ｘ方向（周波数エンコード方向）の信号値の変化を表している。
シーケンスＳＥ_１を実行した後、次のシーケンスＳＥ_２が実行される。

【0037】

シーケンスＳＥ_２が実行されることにより、エコートレインＥＴ_２が収集される。図９に、シーケンスＳＥ_２により収集されたエコートレインＥＴ_２を概略的に示す。エコートレインＥＴ_２の第２エコーＥ_２２〜第ｎエコーＥ_２ｎは、ｋｚ＝−９に配置される。また、変換手段８１は、画像再構成に使用されない第１エコーＥ_２１をｘ方向（周波数エンコード方向）に逆フーリエ変換し、これにより、プロファイルＦ_２が得られる。プロファイルＦ_２は、ｘ方向（周波数エンコード方向）の信号値の変化を表している。

【0038】

以下同様に、各シーケンスを順に実行し、各シーケンスで得られた第１エコーを逆フーリエ変換し、プロファイルを求める。図１０に、最後のシーケンスＳＥ_２１が実行されたときの様子を示す。エコートレインＥＴ_２１の第２エコーＥ_２１，２〜第ｎエコーＥ_２１，ｎは、ｋｚ＝１０に配置される。また、変換手段８１は、画像再構成に使用されない第１エコーＥ_２１，１をｘ方向（周波数エンコード方向）に逆フーリエ変換し、プロファイルＦ_２１を求める。プロファイルＦ_２１は、ｘ方向（周波数エンコード方向）の信号値の変化を表している。

【0039】

このようにして、シーケンスＳＥ_１〜シーケンスＳＥ_２１が実行される。シーケンスＳＥ_１〜シーケンスＳＥ_２１を実行した後、ステップＳＴ２に進む。

【0040】

ステップＳＴ２では、エコートレインＥＴ_１〜ＥＴ_２１の中から、再収集する必要があるエコートレインを特定する。

【0041】

心拍などの体動はエコートレインのエコーに影響を与えることがあり、体動の影響を受けているエコーを画像再構成のエコーとして使用すると、ゴーストなどのアーチファクトの原因となる。そこで、第１の形態では、エコートレインＥＴ_１〜ＥＴ_２１の中に、アーチファクトの原因となる可能性の高いエコートレインが含まれている場合、そのエコートレインを特定する。以下に、この特定方法について説明する。

【0042】

先ず、ステップＳＴ２１において、設定手段８２（図１参照）が、プロファイルＦ_１〜Ｆ_２１の中から基準プロファイルＦ_ｒｅｆとして使用されるプロファイルを設定する。基準プロファイルＦ_ｒｅｆは、次のステップＳＴ２２においてプロファイルの相関係数を求めるときの基準となるプロファイルである。ここでは、プロファイルＦ_１を基準プロファイルＦ_ｒｅｆとして設定する。図１１に、基準プロファイルＦ_ｒｅｆを示す。
基準プロファイルＦ_ｒｅｆを設定した後、ステップＳＴ２２に進む。

【0043】

ステップＳＴ２２では、算出手段８３（図１参照）が、基準プロファイルＦ_ｒｅｆと、プロファイルＦ_１〜Ｆ_２１の各々との相関係数Ｃを求める。図１２に、算出された相関係数Ｃを示す。図１２では、図面のスペースの制約上、代表して、以下の３つの相関係数が示されている。
（１）基準プロファイルＦ_ｒｅｆとプロファイルＦ_１との相関係数Ｃ
（２）基準プロファイルＦ_ｒｅｆとプロファイルＦ_２との相関係数Ｃ
（３）基準プロファイルＦ_ｒｅｆとプロファイルＦ_２１との相関係数Ｃ

【0044】

本形態では、プロファイルＦ_１が基準プロファイルＦ_ｒｅｆとして使用されているので、基準プロファイルＦ_ｒｅｆとプロファイルＦ_１との相関係数Ｃは、Ｃ＝１となる。また、第１の形態では、基準プロファイルＦ_ｒｅｆとプロファイルＦ_２との相関係数Ｃは、Ｃ＝０．９７であり、基準プロファイルＦ_ｒｅｆとプロファイルＦ_２１との相関係数Ｃは、Ｃ＝０．９２であるとする。

【0045】

上記のように相関係数Ｃを求めることにより、第１エコーＥ_１１〜Ｅ_２１，１と相関係数Ｃとの対応関係を知ることができる。図１３に、第１エコーＥ_１１〜Ｅ_２１，１と相関係数Ｃとの対応関係を表す散布図が示されている。散布図の横軸は第１エコーＥ_１１〜Ｅ_２１，１を表しており、縦軸は相関係数Ｃを表す。散布図では、相関係数ＣはＣ＝０．７〜１．０の範囲が示されている。第１のエコーＥ_１１〜Ｅ_２１，１の各々における相関係数Ｃの値は、記号「○」で示されている。相関係数Ｃを求めた後、ステップＳＴ２３に進む。

【0046】

ステップＳＴ２３では、決定手段８４（図１参照）が、相関係数Ｃに基づいて、再収集するエコートレインを決定する。以下に、相関係数Ｃに基づいて、再収集するエコートレインを決定する方法について説明する。
先ず、決定手段８４は、相関係数の度数分布を表す度数分布データを作成する（図１４参照）。

【0047】

図１４は相関係数の度数分布データＤ０を概略的に示す図である。
決定手段８４は、先ず、散布図の相関係数の範囲０．７〜１．０を複数の区間（ビン）に分ける。図１４では、散布図の相関係数の範囲０．７〜１．０を６個の区間、即ち、ビン１〜ビン６に分けたときの様子が示されている。

【0048】

散布図の相関係数の範囲０．７〜１．０をビン１〜ビン６に分けた後、決定手段８４は、各ビンに含まれている相関係数の個数（度数）を求める。したがって、相関係数の度数分布データＤ０を作成することができる。

【0049】

度数分布データＤ０を作成した後、決定手段８４は、ビン１〜ビン６の中から、相関係数の度数が最大となるビンを特定する。図１４では、ビン１の度数は２個、ビン２の度数は１９個、ビン３〜ビン６の度数は０個である。したがって、ビン２が、相関係数の度数が最大となるビンとして特定される。

【0050】

一般的に、被検体から収集した複数のエコーが体動の影響をあまり受けていない場合、これらの複数のエコーの各々のプロファイルを比較すると、プロファイルは類似する傾向がある。つまり、被検体から収集した複数のエコーが体動の影響をあまり受けていない場合、相関係数のばらつきは小さくなる傾向がある。一方、被検体から取集した複数のエコーが体動の影響を大きく受けている場合、これらの複数のエコーのプロファイルはあまり類似せず、相関係数のばらつきは大きくなる傾向がある。したがって、第１のエコーが体動の影響をあまり受けていない場合、相関係数Ｃは同一のビンに集中して現れる。図１４の度数分布データＤ０を参照すると、ビン２の度数が最大であるので、相関係数はビン２に集中して現れていることが分かる。ビン２には、第１エコーＥ_３１〜Ｅ_２１，１の相関係数が含まれている。したがって、第１エコーＥ_３１〜Ｅ_２１，１は体動の影響をそれほど受けていないと考えられる。

【0051】

一方、ビン２とは別のビンに含まれている第１エコーは、体動の影響を大きく受けていると考えられる。第１エコーが体動の影響を大きく受けている場合、その直後に収集される第２エコー〜第ｎエコーも体動の影響を大きく受けていると考えられる。したがって、これらの第２エコー〜第ｎエコーを画像再構成のエコーとして使用すると、ゴーストなどのアーチファクトが発生する原因となる。そこで、第１エコーの相関係数がビン２に含まれていない場合、エコートレインを再収集すると決定する。図１４を参照すると、第１エコーＥ_１１およびＥ_２１がビン２に含まれていないことが分かる。したがって、決定手段８４は、第１エコーＥ_１１を含むエコートレインＥＴ_１と、第１エコーＥ_２１を含むエコートレインＥＴ_２とを再収集すると決定する。
再収集するエコートレインＥＴ_１およびＥＴ_２を決定した後、ステップＳＴ３に進む。

【0052】

ステップＳＴ３では、エコートレインＥＴ_１およびＥＴ_２の再収集が行われる。第１の形態では、スキャンＳＣの期間α１の直後に、エコートレインＥＴ_１およびＥＴ_２を再収集するための期間α２が設けられる（図１５〜図１７参照）。

【0053】

図１５〜図１７は、期間α２にエコートレインＥＴ_１およびＥＴ_２を再収集するときの説明図である。
先ず、エコートレインＥＴ_１を再収集するための１回目のシーケンスＳＥ_１が実行される（図１５参照）。

【0054】

図１５は、エコートレインＥＴ_１を再収集するための１回目のシーケンスＳＥ_１を実行するときの説明図である。
変換手段８１は、１回目のシーケンスＳＥ_１により収集された第１エコーＥ_１１を逆フーリエ変換し、プロファイルＦ_１を求める。そして、算出手段８３は基準プロファイルＦ_ｒｅｆとプロファイルＦ_１との相関係数Ｃを計算する。ここでは、Ｃ＝０．８２であるとする。次に、決定手段８４は、相関係数０．８２がビン２（相関係数０．９０〜０．９５）に含まれているか否かを判断する。図１５を参照すると、相関係数０．８２はビン２には含まれていないことが分かる。上記のように、相関係数がビン２に含まれていない場合、第１エコーは体動の影響を大きく受けていると考えられる。したがって、１回目のシーケンスＳＥ_１により再収集された第２エコーＥ_１２〜第ｎエコーＥ_１ｎは、ｋ空間のデータとして採用することができないと判断される。そこで、期間α２では、相関係数がビン２に含まれるまで、エコートレインＥＴ_１を再収集するためのシーケンスＳＥ_１が繰り返し実行される（図１６参照）。

【0055】

図１６は、期間α２においてシーケンスＳＥ_１が３回実行されたときの様子を示す図である。

【0056】

図１６では、１回目および２回目のシーケンスＳＥ_１で再収集された第１エコーＥ_１１の相関係数はビン２には含まれていないが、３回目のシーケンスＳＥ_１で収集された第１エコーＥ_１１の相関係数Ｃがビン２に含まれたとする。したがって、３回目のシーケンスＳＥ_１で再収集されたエコートレインＥＴ_１の第２エコーＥ_１２〜第ｎエコーＥ_１ｎが、ｋ空間のデータとして採用される。

【0057】

ｋ空間のデータとして採用されるエコートレインＥＴ_１を再収集することができたら、エコートレインＥＴ_２を再収集するためのシーケンスＳＥ_２を実行する（図１７参照）。

【0058】

図１７は、シーケンスＳＥ_２を実行したときの様子を示す図である。
シーケンスＳＥ_２は、第１エコーＥ_２１の相関係数がビン２に含まれるまで繰り返し実行される。シーケンスＳＥ_２については、２回目の再収集で第１エコーＥ_２１の相関係数がビン２に含まれたとする。したがって、２回目の再収集で得られたエコートレインＥＴ_２の第２エコーＥ_２２〜第ｎエコーＥ_２ｎが、ｋ空間のデータとして採用される。

【0059】

このようにして、体動の影響をあまり受けていないエコートレインＥＴ_１およびＥＴ_２を再収集することができる。尚、エコートレインＥＴ_３〜ＥＴ_２１は、ステップＳＴ１において、体動の影響をあまり受けていないエコートレインとして収集されているので、エコートレインＥＴ_３〜ＥＴ_２１の再収集は行われない。

【0060】

したがって、エコートレインＥＴ_１およびＥＴ_２が再収集されたｋ空間のデータＤ′を得ることができる。エコートレインＥＴ_１およびＥＴ_２を再収集した後、ステップＳＴ４に進む。

【0061】

ステップＳＴ４では、図１７に示すｋ空間のデータＤ′を用いて画像再構成が行われ、フローが終了する。

【0062】

第１の形態では、ステップＳＴ１でエコートレインＥＴ_１〜ＥＴ_２１を収集し、エコートレインＥＴ_１〜ＥＴ_２１の第１エコーＥ_１１〜Ｅ_２１，１を逆フーリエ変換し、プロファイルＦ_１〜Ｆ_２１を求める（図１０参照）。そして、基準プロファイルＦ_ｒｅｆを設定し（図１１参照）、基準プロファイルＦ_ｒｅｆとプロファイルＦ_１〜Ｆ_２１の各々との相関係数Ｃを計算する（図１２参照）。第１のエコーが体動の影響をあまり受けていない場合は相関係数Ｃのばらつきは小さいが、第１のエコーが体動の影響を大きく受けている場合は相関係数Ｃのばらつきは大きくなる傾向がある。したがって、度数が最も大きいビンを特定することにより、体動の影響をあまり受けていないエコートレインと、体動の影響を大きく受けているエコートレインとを区別することができる。本形態では、体動の影響を大きく受けているエコートレインはステップＳＴ３において再収集されるので、ゴーストなどのアーチファクトが低減された画像を得ることができる。

【0063】

尚、第１の形態では、プロファイルＦ_１を基準プロファイルＦ_ｒｅｆとして使用している。しかし、プロファイルＦ_１とは別のプロファイルを基準プロファイルＦ_ｒｅｆとして使用してもよい。

【0064】

また、第１の形態では、相関係数Ｃの範囲の下限値を０．７に設定している（例えば、図１７参照）。しかし、ステップＳＴ２２で算出した相関係数に、０．７を下回る値が含まれている場合は、その値が含まれるように、相関係数の範囲の下限値を設定すればよい。例えば、ステップＳＴ２２で算出した相関係数の最小値が、Ｃ＝０．６５の場合、相関係数Ｃの範囲の下限値は０．６５以下の値（例えば、０．６）に設定すればよい。

【0065】

（２）第２の形態
第１の形態では、プロファイルＦ_１が基準プロファイルＦ_ｒｅｆとして使用されている。しかし、プロファイルＦ_１がノイズの影響を強く受けてしまっている場合、本来は再収集しなければならないエコートレインが、再収集する必要がないと判断されてしまい、ゴーストなどのアーチファクトを十分に低減することができないことがある。したがって、基準プロファイルＦ_ｒｅｆは、ノイズの影響をあまり受けていないことが望ましい。そこで、第２の形態では、ノイズの影響をあまり受けていないプロファイルを基準プロファイルとして設定する例について説明する。

【0066】

図１８は、第２の形態のＭＲ装置の概略図である。
第２の形態のＭＲ装置２００は、第１の形態のＭＲ装置１００（図１参照）と比較すると、制御部８が異なっているが、その他の構成は同じである。したがって、第２の形態のＭＲ装置２００については、主に、制御部８について説明する。制御部８は、変換手段８０１〜決定手段８０５などを有している。

【0067】

変換手段８０１は、第１エコーを逆フーリエ変換し、後述するプロファイルを求める。
候補選択手段８０２は、複数のプロファイルの中から、基準プロファイルの候補となる２つ以上のプロファイルを選択する。
算出手段８０３は、基準プロファイルの候補となるプロファイルと各プロファイルとの相関係数を算出する。
設定手段８０４は、相関係数に基づいて、基準プロファイルの候補の中から、基準プロファイルとして使用されるプロファイルを設定する。設定手段８０４は、度数分布データ作成手段８０４ａと、プファイル選択手段８０４ｂとを有している。度数分布データ作成手段８０４ａは、基準プロファイルの候補ごとに、相関係数の度数分布を表す度数分布データを作成する。プロファイル選択手段８０４ｂは、基準プロファイルの候補の中から、基準プロファイルとして使用されるプロファイルを選択する。
決定手段８０５は、相関係数に基づいて、再収集するエコートレインを決定する。

【0068】

制御部８は、変換手段８０１〜決定手段８０５を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。

【0069】

第２の形態では、制御部８は上記のように構成されている。次に、第２の形態においてスキャンを実行するときのフローについて説明する。

【0070】

図１９は、第２の形態のフローの説明図である。
ステップＳＴ１は、第１の形態と同様に、シーケンスＳＥ_１〜シーケンスＳＥ_２１が実行される。したがって、第１の形態において図８〜図１０を参照しながら説明したように、シーケンスＳＥ_１〜シーケンスＳＥ_２１が実行されることにより、エコートレインＥＴ_１〜ＥＴ_２１が取得される。変換手段８０１（図１８参照）は、第１エコーＥ_１１〜Ｅ_２１，１を逆フーリエ変換する。これにより、プロファイルＦ_１〜Ｆ_２１が得られる。シーケンスＳＥ_１〜シーケンスＳＥ_２１を実行した後、ステップＳＴ２に進む。

【0071】

ステップＳＴ２では、エコートレインＥＴ_１〜ＥＴ_２１の中から、再収集する必要があるエコートレインを特定する。ステップＳＴ２はステップＳＴ２０１〜ＳＴ２０３を有している。以下、各ステップＳＴ２０１〜ＳＴ２０３について説明する。

【0072】

ステップＳＴ２０１では、候補選択手段８０２（図１８参照）が、プロファイルＦ_１〜Ｆ_２０の中から基準プロファイルＦ_ｒｅｆの候補となるプロファイルを選択する。図２０に、基準プロファイルＦ_ｒｅｆの候補として選択されたプロファイルを示す。第２の形態では、プロファイルＦ_１およびＦ_２が、基準プロファイルＦ_ｒｅｆの候補として選択されるが、３つ以上のプロファイルを、基準プロファイルＦ_ｒｅｆの候補として選択してもよい。基準プロファイルＦ_ｒｅｆの候補Ｆ_１およびＦ_２を選択した後、ステップＳＴ２０２に進む。

【0073】

ステップＳＴ２０２では、基準プロファイルＦ_ｒｅｆの候補として選択されたプロファイルＦ_１およびＦ_２のうち、どちらのプロファイルがノイズの影響が小さいかを判断する。ステップＳＴ２０２は、ステップＳＴ２０２ａ〜ＳＴ２０２ｅを有しているので、各ステップＳＴ２０２ａ〜ＳＴ２０２ｅについて説明する。

【0074】

ステップＳＴ２０２ａでは、算出手段８０３（図１８参照）が、基準プロファイルＦ_ｒｅｆの候補として選択されたプロファイルＦ_１と、プロファイルＦ_１〜Ｆ_２１の各々との相関係数を求める。図２１に、算出された相関係数Ｃを示す。また、図２１の左下には、第１のエコーＥ_１１〜Ｅ_２１，１と相関係数Ｃとの対応関係を表す散布図が示されている。散布図の横軸は第１のエコーＥ_１１〜Ｅ_２１，１を表しており、縦軸は相関係数Ｃを表す。散布図では、相関係数ＣはＣ＝０．７〜１．０の範囲が示されている。第１のエコーＥ_１１〜Ｅ_２１，１の各々における相関係数Ｃの値は、記号「○」で示されている。相関係数Ｃを求めた後、ステップＳＴ２０２ｂに進む。

【0075】

ステップＳＴ２０２ｂでは、度数分布データ作成手段８０４ａ（図１８参照）が、相関係数の度数分布を表す度数分布データを作成する（図２２参照）。

【0076】

図２２は、相関係数の度数分布データＤ１を概略的に示す図である。
度数分布データ作成手段８０４ａは、先ず、散布図の相関係数の範囲０．７〜１．０を複数の区間（ビン）に分ける。図２２では、散布図の相関係数の範囲０．７〜１．０を６個の区間、即ち、ビン１〜ビン６に分けたときの様子が示されている。

【0077】

散布図の相関係数の範囲０．７〜１．０をビン１〜ビン６に分けた後、度数分布データ作成手段８０４ａは、各ビンに含まれている相関係数の個数（度数）を求める。このようにして、度数分布データＤ１を作成することができる。図２２では、ビン１の度数は２個、ビン２の度数は１９個、ビン３〜ビン６の度数は０個であるので、ビン２の度数が最大になっている。度数分布データＤ１を作成した後、ステップＳＴ２０２ｃに進む。

【0078】

ステップＳＴ２０２ｃでは、算出手段８０３が、基準プロファイルＦ_ｒｅｆの候補として選択されたプロファイルＦ_２と、プロファイルＦ_１〜Ｆ_２１の各々との相関係数を求める。図２３に、算出された相関係数Ｃを示す。また、図２３の左下には、第１のエコーＥ_１１〜Ｅ_２１，１と相関係数Ｃとの対応関係を表す散布図が示されている。散布図の横軸は第１のエコーＥ_１１〜Ｅ_２１，１を表しており、縦軸は相関係数Ｃを表す。散布図では、相関係数ＣはＣ＝０．７〜１．０の範囲が示されている。第１のエコーＥ_１１〜Ｅ_２１，１の各々における相関係数Ｃの値は、記号「○」で示されている。相関係数Ｃを求めた後、ステップＳＴ２０２ｄに進む。

【0079】

ステップＳＴ２０２ｄでは、度数分布データ作成手段８０４ａが、相関係数の度数分布を表す度数分布データを作成する（図２４参照）。

【0080】

図２４は、相関係数の度数分布データＤ２を概略的に示す図である。
度数分布データ作成手段８０４ａは、先ず、散布図の相関係数の範囲０．７〜１．０を複数の区間（ビン）に分ける。図２４では、散布図の相関係数の範囲０．７〜１．０を６個の区間、即ち、ビン１〜ビン６に分けたときの様子が示されている。

【0081】

散布図の相関係数の範囲０．７〜１．０をビン１〜ビン６に分けた後、度数分布データ作成手段８０４ａは、各ビンに含まれている相関係数の個数（度数）を求める。このようにして、度数分布データＤ２を作成することができる。図２４では、ビン１の度数は１９個、ビン２の度数は２個、ビン３〜ビン６の度数は０個であるので、ビン１の度数が最大になっている。度数分布データＤ２を作成した後、ステップＳＴ２０２ｅに進む。

【0082】

ステップＳＴ２０２ｅでは、プロファイル選択手段８０４ｂが、度数分布データＤ１（図２２参照）と度数分布データＤ２（図２４参照）とに基づいて、プロファイルＦ_１およびＦ_２の中から、ノイズの影響が小さい方のプロファイルを選択する（図２５参照）。

【0083】

図２５は、ノイズの影響が小さい方のプロファイルを選択するときの説明図である。
図２５には、図２２に示す度数分布データＤ１（基準プロファイルＦ_ｒｅｆの候補：プロファイルＦ_１）と、図２４に示す度数分布データＤ２（基準プロファイルＦ_ｒｅｆの候補：プロファイルＦ_２）とが示されている。

【0084】

プロファイル選択手段８０４ｂは、先ず、度数分布データＤ１の中から、度数が最大となるビン２に含まれる相関係数Ｃを抽出する。そして、ビン２に含まれる相関係数の平均値を求める。ここでは、ビン２に含まれる相関係数の平均値を０．９３とする。

【0085】

次に、プロファイル選択手段８０４ｂは、度数分布データＤ２の中から、度数が最大となるビン１に含まれる相関係数Ｃを抽出する。そして、ビン１に含まれる相関係数の平均値を求める。ここでは、ビン１に含まれる相関係数の平均値を０．９７とする。

【0086】

相関係数の平均値を求めた後、プロファイル選択手段８０４ｂは、これらの相関係数の平均値を比較する。一般的に、プロファイルがノイズの影響をあまり受けていない場合、相関係数は大きい値（Ｃ＝１に近い値）になる傾向がある。したがって、相関係数の平均値が大きいほど、プロファイルはノイズの影響をあまり受けていないと考えることができる。度数分布データＤ１と度数分布データＤ２とを比較すると、度数分布データＤ２の相関係数の平均値（０．９７）は、度数分布データＤ１の相関係数の平均値（０．９３）よりも大きい。したがって、プロファイルＦ_１とプロファイルＦ_２とを比較すると、プロファイルＦ_２は、プロファイルＦ_１よりも、ノイズの影響を受けていないと考えられる。そこで、プロファイル選択手段８０４ｂは、プロファイルＦ_２を、ノイズの影響をあまり受けていないプロファイルとして選択する。このようにして選択されたプロファイルＦ_２が基準プロファイルＦ_ｒｅｆとして設定される。基準プロファイルＦ_ｒｅｆを設定した後、ステップＳＴ２０３に進む。

【0087】

ステップＳＴ２０３では、決定手段８０５（図１８参照）が、再収集するエコートレインを決定する。再収集するエコートレインは、以下の手順で決定する。

【0088】

決定手段８０５は、先ず、２つの度数分布データＤ１およびＤ２のうちの一方の度数分布データを、再収集するエコートレインを決定するときに使用する度数分布データとして選択する。第２の形態では、プロファイルＦ_２が基準プロファイルＦ_ｒｅｆとして設定されたので（ステップＳＴ２０２ｅ参照）、決定手段８０５は、度数分布データＤ２を、再収集するエコートレインを決定するときに使用する度数分布データとして選択する。

【0089】

度数分布データＤ２を選択した後、決定手段８０５は、度数分布データＤ２のビン１〜ビン６の中で、度数が最大となるビンを特定する。図２４を参照すると、度数分布データＤ２では、ビン１の度数が最大となっている。第１の形態で説明したように、第１のエコーが体動の影響をあまり受けていない場合、相関係数Ｃは同一のビンに集中して現れるので、ビン１に含まれる第１エコーは体動の影響をそれほど受けていないと考えられる。

【0090】

一方、ビン１とは別のビンに含まれている第１エコーは、体動の影響を大きく受けていると考えられる。図２４の度数分布データＤ２を参照すると、第１エコ−Ｅ_１１およびＥ_２１，１がビン１に含まれていないことが分かる。したがって、決定手段８０５は、第１エコーＥ_１１を含むエコートレインＥＴ_１と、第１エコーＥ_２１，１を含むエコートレインＥＴ_２１とを再収集すると決定する。
再収集するエコートレインを決定した後、ステップＳＴ３に進む。

【0091】

ステップＳＴ３では、相関係数が度数分布データＤ２のビン１に含まれるまで、エコートレインＥＴ_１およびＥＴ_２１が再収集される。エコートレインを再収集した後、ステップＳＴ４に進み、画像再構成を行いフローを終了する。

【0092】

第２の形態では、ノイズの影響をあまり受けていないプロファイルを基準プロファイルとして用いることができるので、エコートレインを再収集するか否かの判断を、より的確に行うことができる。

【0093】

（３）第３の形態
以下に、第３の形態について説明する。尚、第３の形態のＭＲ装置の基本的な構成は、第１の形態のＭＲ装置１００と同じである。したがって、必要に応じて図１を参照しながら、第３の形態について説明する。

【0094】

図２６は、第３の形態におけるスキャンＳＣの説明図である。
第３の形態では、スキャンＳＣは、３つの期間β１〜β３に分けられる。

【0095】

期間β１は、シーケンスＳＥ_１〜ＳＥ_１０が実行される期間である。
期間β２は、エコートレインを再収集するためのシーケンスが実行される期間である。
期間β３は、シーケンスＳＥ_１１〜ＳＥ_２１が実行され、必要に応じてエコートレインを再収集するためのシーケンスが実行される期間である。以下、第３の形態で実行されるスキャンＳＣについて説明する。

【0096】

図２７は、第３の形態におけるスキャンＳＣを実行するときのフローを示す図である。
ステップＳＴで１は、スキャンＳＣの期間β１におけるシーケンスが実行される（図２８参照）。

【0097】

図２８〜図３０は、スキャンＳＣの期間β１において実行されるシーケンスの説明図である。

【0098】

期間β１では、先ず、シーケンスＳＥ_１が実行される。図２８に、シーケンスＳＥ_１が実行されるときの様子を示す。シーケンスＳＥ_１が実行されることにより、エコートレインＥＴ_１が収集される。エコートレインＥＴ_１の第２エコーＥ_１２〜第ｎエコーＥ_１ｎは、ｋｚ＝−１０に配置される。また、変換手段８１（図１参照）は、第１エコーＥ_１１をｘ方向（周波数エンコード方向）に逆フーリエ変換する。これにより、プロファイルＦ_１が得られる。プロファイルＦ_１は、ｘ方向（周波数エンコード方向）の信号値の変化を表している。
シーケンスＳＥ_１を実行した後、次のシーケンスＳＥ_２が実行される。

【0099】

シーケンスＳＥ_２が実行されることにより、エコートレインＥＴ_２が収集される。図２９に、シーケンスＳＥ_２により収集されたエコートレインＥＴ_２を概略的に示す。エコートレインＥＴ_２の第２エコーＥ_２２〜第ｎエコーＥ_２ｎは、ｋｚ＝−９に配置される。また、第１エコーＥ_２１をｘ方向（周波数エンコード方向）に逆フーリエ変換することにより、プロファイルＦ_２が得られる。プロファイルＦ_２は、ｘ方向（周波数エンコード方向）の信号値の変化を表している。

【0100】

以下同様に、各シーケンスを順に実行し、各シーケンスで得られた第１エコーを逆フーリエ変換し、プロファイルを求める。図３０に、シーケンスＳＥ_１０が実行されたときの様子を示す。エコートレインＥＴ_１０の第２エコーＥ_１０，２〜第ｎエコーＥ_１０，ｎは、ｋｚ＝−１に配置される。また、第１エコーＥ_１０，１をｘ方向（周波数エンコード方向）に逆フーリエ変換することにより、プロファイルＦ_１０が得られる。プロファイルＦ_１０は、ｘ方向（周波数エンコード方向）の信号値の変化を表している。

【0101】

このようにして、シーケンスＳＥ_１〜シーケンスＳＥ_１０が実行される。シーケンスＳＥ_１〜シーケンスＳＥ_１０を実行した後、ステップＳＴ２に進む。

【0102】

ステップＳＴ２では、エコートレインＥＴ_１〜ＥＴ_１０の中から、再収集する必要があるエコートレインを特定する。以下、ステップＳＴ２に含まれるステップＳＴ２１〜ＳＴ２３について順に説明する。

【0103】

先ず、ステップＳＴ２１において、設定手段８２（図１参照）が、プロファイルＦ_１〜Ｆ_１０の中から、基準プロファイルＦ_ｒｅｆとして使用されるプロファイルを設定する。第３の形態では、第１の形態と同様に、プロファイルＦ_１が基準プロファイルＦ_ｒｅｆとして設定されるとする。図３１に、基準プロファイルＦ_ｒｅｆを示す。

【0104】

基準プロファイルＦ_ｒｅｆを設定した後、ステップＳＴ２２に進む。

【0105】

ステップＳＴ２２では、算出手段８３（図１参照）が、基準プロファイルＦ_ｒｅｆと、プロファイルＦ_１〜Ｆ_１０の各々との相関係数Ｃを求める。図３２に、算出された相関係数Ｃを示す。図３２では、図面のスペースの制約上、代表して、以下の３つの相関係数が示されている。
（１）基準プロファイルＦ_ｒｅｆとプロファイルＦ_１との相関係数Ｃ
（２）基準プロファイルＦ_ｒｅｆとプロファイルＦ_２との相関係数Ｃ
（３）基準プロファイルＦ_ｒｅｆとプロファイルＦ_１０との相関係数Ｃ

【0106】

本形態では、プロファイルＦ_１が基準プロファイルＦ_ｒｅｆとして使用されているので、基準プロファイルＦ_ｒｅｆとプロファイルＦ_１との相関係数Ｃは、Ｃ＝１となる。また、第３の形態では、基準プロファイルＦ_ｒｅｆとプロファイルＦ_２との相関係数Ｃは、Ｃ＝０．９７であり、基準プロファイルＦ_ｒｅｆとプロファイルＦ_１０との相関係数Ｃは、Ｃ＝０．９２であるとする。

【0107】

上記のように相関係数Ｃを求めることにより、第１エコーＥ_１１〜Ｅ_１０，１と相関係数Ｃとの対応関係を知ることができる。図３３に、第１エコーＥ_１１〜Ｅ_１０，１と相関係数Ｃとの対応関係を表す散布図が示されている。散布図の横軸は第１エコーＥ_１１〜Ｅ_１０，１を表しており、縦軸は相関係数Ｃを表す。散布図では、相関係数ＣはＣ＝０．７〜１．０の範囲が示されている。第１のエコーＥ_１１〜Ｅ_１０，１の各々における相関係数Ｃの値は、記号「○」で示されている。相関係数Ｃを求めた後、ステップＳＴ２３に進む。

【0108】

ステップＳＴ２３では、決定手段８４（図１参照）が、相関係数に基づいて、再収集するエコートレインを決定する。以下に、相関係数に基づいて、再収集するエコートレインを決定する方法について説明する。

【0109】

先ず、決定手段８４は、相関係数の度数分布を表す度数分布データを作成する（図３４参照）。

【0110】

図３４は相関係数の度数分布データＤ３を概略的に示す図である。
決定手段８４は、先ず、散布図の相関係数の範囲０．７〜１．０を複数の区間（ビン）に分ける。図３４では、散布図の相関係数の範囲０．７〜１．０を６個の区間、即ち、ビン１〜ビン６に分けたときの様子が示されている。

【0111】

散布図の相関係数の範囲０．７〜１．０をビン１〜ビン６に分けた後、決定手段８４は、各ビンに含まれている相関係数の個数（度数）を求める。図３４では、ビン１の度数は２個、ビン２の度数は８個、ビン３〜ビン６の度数は０個である。したがって、図３４では、ビン２の度数が最大になっている。

【0112】

第１の形態で説明したように、第１のエコーが体動の影響をあまり受けていない場合、相関係数Ｃは同一のビンに集中して現れる。図３４の度数分布データＤ３を参照すると、ビン２の度数が最大であるので、相関係数はビン２に集中して現れていることが分かる。ビン２には、第１エコーＥ_３１〜Ｅ_１０，１の相関係数が含まれている。したがって、第１エコーＥ_３１〜Ｅ_１０，１は体動の影響をそれほど受けていないと考えられる。

【0113】

一方、第１エコーＥ_１１およびＥ_２１はビン２に含まれていないので、体動の影響を受けていると考えられる。したがって、決定手段８４は、第１エコーＥ_１１を含むエコートレインＥＴ_１と、第１エコーＥ_２１を含むエコートレインＥＴ_２とを再収集すると決定する。
再収集するエコートレインＥＴ_１およびＥＴ_２を決定した後、ステップＳＴ３に進む。

【0114】

ステップＳＴ３では、エコートレインＥＴ_１およびＥＴ_２の再収集が行われる。第３の形態では、スキャンＳＣの期間β１の直後に、エコートレインＥＴ_１およびＥＴ_２を再収集するための期間β２が設けられる（図３５〜図３７参照）。

【0115】

図３５〜図３７は、期間β２にエコートレインＥＴ_１およびＥＴ_２を再収集するときの説明図である。
先ず、エコートレインＥＴ_１を再収集するための１回目のシーケンスＳＥ_１が実行される（図３５参照）。

【0116】

図３５は、エコートレインＥＴ_１を再収集するための１回目のシーケンスＳＥ_１を実行するときの説明図である。

【0117】

変換手段８１は、１回目のシーケンスＳＥ_１により収集された第１エコーＥ_１１を逆フーリエ変換し、プロファイルＦ_１を求める。そして、算出手段８３は基準プロファイルＦ_ｒｅｆとプロファイルＦ_１との相関係数Ｃを計算する。ここでは、Ｃ＝０．８２であるとする。次に、決定手段８４は、相関係数０．８２がビン２（相関係数０．９０〜０．９５）に含まれているか否かを判断する。図３５を参照すると、相関係数０．８２はビン２には含まれていないことが分かる。上記のように、相関係数がビン２に含まれていない場合、第１エコーは体動の影響を大きく受けていると考えられる。したがって、１回目のシーケンスＳＥ_１により再収集された第２エコーＥ_１２〜第ｎエコーＥ_１ｎは、ｋ空間のデータとして採用することができない。そこで、期間β２では、相関係数がビン２に含まれるまで、エコートレインＥＴ_１を再収集するためのシーケンスＳＥ_１が繰り返し実行される（図３６参照）。

【0118】

図３６は、期間β２においてシーケンスＳＥ_１が３回実行されたときの様子を示す図である。
図３６では、１回目および２回目のシーケンスＳＥ_１で再収集された第１エコーＥ_１１の相関係数はビン２には含まれていないが、３回目のシーケンスＳＥ_１で収集された第１エコーＥ_１１の相関係数Ｃがビン２に含まれたとする。したがって、３回目のシーケンスＳＥ_１で再収集されたエコートレインＥＴ_１の第２エコーＥ_１２〜第ｎエコーＥ_１ｎが、ｋ空間のデータとして採用される。

【0119】

ｋ空間のデータとして採用されるエコートレインＥＴ_１を再収集することができたら、エコートレインＥＴ_２を再収集するためのシーケンスＳＥ_２を実行する（図３７参照）。

【0120】

図３７は、シーケンスＳＥ_２を実行したときの様子を示す図である。
シーケンスＳＥ_２は、第１エコーＥ_２１の相関係数がビン２に含まれるまで繰り返し実行される。シーケンスＳＥ_２については、２回目の再収集で第１エコーＥ_２１の相関係数がビン２に含まれたとする。したがって、２回目の再収集で得られたエコートレインＥＴ_２の第２エコーＥ_２２〜第ｎエコーＥ_２ｎが、ｋ空間のデータとして採用される。

【0121】

このようにして、体動の影響をあまり受けていないエコートレインＥＴ_１およびＥＴ_２を再収集することができる。尚、エコートレインＥＴ_３〜ＥＴ_１０は、ステップＳＴ１において、体動の影響をあまり受けていないエコートレインとして収集されているので、エコートレインＥＴ_３〜ＥＴ_１０の再収集は行われない。

【0122】

エコートレインＥＴ_１およびＥＴ_２を再収集した後、ステップＳＴ４に進む。

【0123】

ステップＳＴ３では、ｋｚ＝０〜１０に配置されるデータを収集するためのシーケンスＳＥ_１１〜ＳＥ_２１が実行される。第３の形態では、スキャンＳＣの期間β２の直後に、シーケンスＳＥ_１１〜ＳＥ_２１を実行するための期間β３が設けられる（図３８〜図４０参照）。

【0124】

図３８〜図４０は、期間β３で実行されるシーケンスの説明図である。
期間β３では、先ず、ｋｚ＝０のデータを収集するためのシーケンスＳＥ_１１が実行される（図３８参照）。

【0125】

図３８はシーケンスＳＥ_１１を実行するときの説明図である。
シーケンスＳＥ_１１が実行されることにより、エコートレインＥＴ_１１が収集される。変換手段８１は、シーケンスＳＥ_１１により収集された第１エコーＥ_１１，１を逆フーリエ変換し、プロファイルＦ_１１を求める。そして、算出手段８３は基準プロファイルＦ_ｒｅｆとプロファイルＦ_１１との相関係数Ｃを計算する。ここでは、Ｃ＝０．９７であるとする。相関係数０．９７はビン２（相関係数０．９０〜０．９５）に含まれていなので、第１エコーＥ_１１，１は体動の影響を大きく受けていると考えられる。したがって、シーケンスＳＥ_１１により収集された第２エコーＥ_１１，２〜第ｎエコーＥ_１１，ｎは、ｋ空間のデータとして採用することができない。そこで、期間β３では、相関係数がビン２に含まれるまで、エコートレインを再収集するためのシーケンスＳＥ_１１が実行される（図３９参照）。

【0126】

図３９は、期間β３においてシーケンスＳＥ_１１が３回実行されたときの様子を示す図である。
図３９では、１回目および２回目のシーケンスＳＥ_１１で再収集された第１エコーＥ_１１，１の相関係数はビン２には含まれていないが、３回目のシーケンスＳＥ_１１で収集された第１エコーＥ_１１１の相関係数Ｃがビン２に含まれたとする。したがって、３回目のシーケンスＳＥ_１１で再収集されたエコートレインＥＴ_１１の第２エコーＥ_１１，２〜第ｎエコーＥ_１１，ｎが、ｋｚ＝０のデータとして採用される。

【0127】

ｋｚ＝０のデータを取得した後、残りのｋｚ＝１〜１０のデータを収集するためのシーケンスＳＥ_１２〜ＳＥ_２１を実行する。図４０に、シーケンスＳＥ_１２〜ＳＥ_２１を実行したときの様子を概略的に示す。シーケンスＳＥ_１２〜ＳＥ_２１を実行する場合も、シーケンスＳＥ_１１と同様に、相関係数Ｃがビン２に含まれるまでシーケンスを実行する。したがって、体動の影響をあまり受けていないエコートレインがｋ空間のデータとして採用される。

【0128】

上記のようにして、ｋ空間のデータを取得することができる。ｋ空間のデータを取得した後、ステップＳＴ５に進み、画像再構成が行われ、フローが終了する。

【0129】

第３の形態では、ｋ空間の一部の領域（ｋｚ＝−９〜−１）のデータが収集された時点で、度数分布データＤ３（図３４）を作成している。このように、ｋ空間の全領域のデータを収集する前に度数分布データＤ３を作成しても、エコートレインが体動の影響を大きく受けているか否かを判断することができる。したがって、体動の影響を大きく受けているエコートレインは再収集されるので、ゴーストなどのアーチファクトが低減された画像を取得することができる。

【0130】

尚、第３の形態においても、第２の形態のように、基準プロファイルの候補となる２つ以上のプロファイルを選択し、選択された２つ以上のプロファイルの中から、基準プロファイルとして使用されるプロファイルを設定してもよい。

【0131】

また、第３の形態は、第１の形態と比較すると、度数分布データを作成するために必要な相関係数の個数が少なくて済むので、より短時間で度数分布データを作成することができる。

【0132】

尚、第１〜第３の形態では、第１エコーを、画像再構成に使用されないエコーとしている。しかし、画像再構成に使用されないエコーは、第１エコーに限定されることはなく、第１エコー〜第ｎエコーのうちのどのエコーであってもよい。

【0133】

また、第１〜第３の形態では、第１エコーと相関係数との関係を表す度数分布データを作成し、再収集するエコートレインを決定している。しかし、画像再構成に使用される第２エコー〜第ｎエコーの中から選択された第ｉエコー（２≦ｉ≦ｎ）のプロファイルを求め、第ｉエコーと相関係数との関係を表す度数分布データに基づいて、再収集するエコートレインを決定してもよい。ただし、第２エコー〜第ｎエコーは位相エンコード勾配が印加されたときに収集されるエコーであるので、位相分散などの影響で、エコーのピーク値は小さくなると考えられる。エコーのピーク値が小さい場合、相関係数の誤差が大きくなることが考えられるので、エコーのピーク値はできるだけ大きいことが望ましい。したがって、位相分散などの影響をできるだけ受けないように、位相エンコード勾配が印加されていないときに収集された第１エコーを用いて度数分布データを作成することが望ましい。

【符号の説明】

【0134】

２ マグネット
３ テーブル
３ａ クレードル
４ 受信コイル
５ 送信器
６ 勾配磁場電源
７ 制御部
８ 操作部
９ 表示部
１０ 被検体
２１ ボア
７１ ピーク値検出手段
７２、７１３ 決定手段
１００、２００ ＭＲ装置
７１１ 変換手段
７１２ 検出手段

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】

【図18】

【図19】

【図20】

【図21】

【図22】

【図23】

【図24】

【図25】

【図26】

【図27】

【図28】

【図29】

【図30】

【図31】

【図32】

【図33】

【図34】

【図35】

【図36】

【図37】

【図38】

【図39】

【図40】