(57)【要約】      （修正有）【課題】信号消失などのアーチファクトを低減することができる磁気共鳴装置を提供する。【解決手段】スキャンＳＣを実行し、エコートレインＥＴ１〜ＥＴｍを収集する。そして、エコートレインＥＴ１〜ＥＴｍの第１エコーＥ１１〜Ｅｍ１のピーク値Ｐ１１〜Ｐｍ１を求める。次に、第１エコーＥ１１〜Ｅｍ１の中から、ピーク値が閾値Ｔより小さい第１エコーＥ２１およびＥｍ１を特定する。決定手段は、ピーク値が閾値Ｔより小さい第１エコーＥ２１含むエコートレインＥＴ２と、ピーク値が閾値Ｔより小さい第１エコーＥｍ１含むエコートレインＥＴｍとを再収集すると決定する。【選択図】図１０
2014147466.tif 000002
【特許請求の範囲】
【請求項１】
画像再構成に使用されるエコーと画像再構成に使用されないエコーとを含むエコートレインを収集するためのシーケンスを複数回実行する磁気共鳴装置であって、
画像再構成に使用されないエコーに基づいて、被検体の体動がエコーにどの程度の影響を与えているかを表す指標を求める手段と、
前記指標に基づいて、複数回のシーケンスにより収集された複数のエコートレインの中から、再収集するエコートレインを決定する決定手段と、
を有する磁気共鳴装置。
【請求項２】
エコートレインを再収集するための追加のスキャンを実行する、請求項１に記載の磁気共鳴装置。
【請求項３】
前記指標は、画像再構成に使用されないエコーのピーク値である、請求項２に記載の磁気共鳴装置。
【請求項４】
前記決定手段は、
画像再構成に使用されないエコーのピーク値が大きいか小さいかを判断するための閾値を求め、
前記ピーク値と前記閾値との比較結果に基づいて、前記複数のエコートレインの中から、再収集するエコートレインを決定する、請求項３に記載の磁気共鳴装置。
【請求項５】
操作者は、前記追加のスキャンが実行される前、又は前記追加のスキャンが実行されている間に、前記閾値を変更することができる、請求項４に記載の磁気共鳴装置。
【請求項６】
前記指標は、画像再構成に使用されないエコーの面積である、請求項２に記載の磁気共鳴装置。
【請求項７】
前記決定手段は、
画像再構成に使用されないエコーの面積が大きいか小さいかを判断するための閾値を求め、
前記面積と前記閾値との比較結果に基づいて、前記複数のエコートレインの中から、再収集するエコートレインを決定する、請求項６に記載の磁気共鳴装置。
【請求項８】
操作者は、前記追加のスキャンが実行される前、又は前記追加のスキャンが実行されている間に、前記閾値を変更することができる、請求項７に記載の磁気共鳴装置。
【請求項９】
前記体動は心拍動である、請求項１〜８のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１０】
前記指標を求める手段は、
画像再構成に使用されないエコーを実空間のプロファイルに変換し、前記プロファイルの所定の位置における信号値を、前記指標として求める、請求項２に記載の磁気共鳴装置。
【請求項１１】
前記決定手段は、
前記プロファイルの前記所定の位置における信号値が大きいか小さいかを判断するための閾値を求め、
前記所定の位置における信号値と前記閾値との比較結果に基づいて、前記複数のエコートレインの中から、再収集するエコートレインを決定する、請求項１０に記載の磁気共鳴装置。
【請求項１２】
操作者は、前記追加のスキャンが実行される前、又は前記追加のスキャンが実行されている間に、前記閾値を変更することができる、請求項１１に記載の磁気共鳴装置。
【請求項１３】
前記シーケンスは、肝臓を含む撮影領域の画像データを収集するために実行され、
前記所定の位置は、肝臓のうちの心臓に近い部位の位置である、請求項１０〜１２のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１４】
画像再構成に使用されるエコーと画像再構成に使用されないエコーとを含むエコートレインを収集するためのシーケンスを複数回実行する磁気共鳴装置のプログラムであって、
画像再構成に使用されないエコーに基づいて、被検体の体動がエコーにどの程度の影響を与えているかを表す指標を求める処理と、
前記指標に基づいて、複数回のシーケンスにより収集された複数のエコートレインの中から、再収集するエコートレインを決定する決定処理と、
を計算機に実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エコートレインを収集するためのシーケンスを複数回実行する磁気共鳴装置、およびこの磁気共鳴装置に適用可能なプログラムに関する。
【背景技術】
【０００２】
拡散の情報を取得する方法として、スピンの動きを検出するためのＭＰＧ（Motion Probing Gradient）を用いた方法が知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１２−１５７６８７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ＭＰＧを用いて拡散の情報を得る方法の一例として、ＭＰＧをプリパルスとして印加したdiffusion preparationと呼ばれる手法がある。この手法は、例えば、ＦＳＥ（Fast Spin Echo）を収集するパルスシーケンスを実行する場合に適用されており、比較的容易に拡散の情報を得ることができる。
【０００５】
しかし、上記の手法では、心拍動などの体動の影響を受けて、画像データに信号消失などのアーチファクトが現れることがある。特に、ＭＰＧのｂ値が高い場合、信号消失は現れやすくなる。したがって、信号消失などのアーチファクトを低減することができる技術が望まれている。
【課題を解決するための手段】
【０００６】
本発明の第１の観点は、画像再構成に使用されるエコーと画像再構成に使用されないエコーとを含むエコートレインを収集するためのシーケンスを複数回実行する磁気共鳴装置であって、
画像再構成に使用されないエコーに基づいて、被検体の体動がエコーにどの程度の影響を与えているかを表す指標を求める手段と、
前記指標に基づいて、複数回のシーケンスにより収集された複数のエコートレインの中から、再収集するエコートレインを決定する決定手段と、
を有する磁気共鳴装置である。
【０００７】
本発明の第２の観点は、画像再構成に使用されるエコーと画像再構成に使用されないエコーとを含むエコートレインを収集するためのシーケンスを複数回実行する磁気共鳴装置のプログラムであって、
画像再構成に使用されないエコーに基づいて、被検体の体動がエコーにどの程度の影響を与えているかを表す指標を求める処理と、
前記指標に基づいて、複数回のシーケンスにより収集された複数のエコートレインの中から、再収集するエコートレインを決定する決定処理と、
を計算機に実行させるためのプログラムである。
【発明の効果】
【０００８】
画像再構成に使用されないエコーに基づいて、被検体の体動がエコーにどの程度の影響を与えているかを表す指標を求め、この指標を用いて、複数のエコートレインの中から、再収集するエコートレインを決定している。したがって、アーチファクトが低減された画像データを得ることができる。
【図面の簡単な説明】
【０００９】
【図１】本発明の第１の形態の磁気共鳴装置の概略図である。
【図２】第１の形態における撮影領域を概略的に示す図である。
【図３】撮影領域Ｒからデータを収集するために実行されるスキャンＳＣの説明図である。
【図４】スキャンＳＣにおいて実行されるシーケンスの説明図である。
【図５】シーケンスＳＥ１〜ＳＥｍにより収集されたエコーをｋ空間に埋めるときの様子を示す図である。
【図６】第１の形態におけるスキャンを実行するときのフローを示す図である。
【図７】シーケンスＳＥ１により収集されたエコートレインＥＴ１を概略的に示す図である。
【図８】シーケンスＳＥ２により収集されたエコートレインＥＴ２を概略的に示す図である。
【図９】シーケンスＳＥ１〜ＳＥｍにより収集されたエコートレインＥＴ１〜ＥＴｍと、各エコートレインの第１エコーのピーク値Ｐ１１〜Ｐｍ１を示す図である。
【図１０】エコートレインＥＴ１〜ＥＴｍの中から再収集する必要があるエコートレインを決定する方法の一例の説明図である。
【図１１】第２エコーＥ１２〜Ｅｍ２のピーク値に基づいて、エコートレインＥＴ１〜ＥＴｍの中から、再収集するエコートレインを決定するときの説明図である。
【図１２】追加のスキャンの説明図である。
【図１３】エコーの面積を求める方法の一例の説明図である。
【図１４】第２の形態のＭＲ装置の概略図である。
【図１５】スキャンＳＣを概略的に示す図である。
【図１６】シーケンスＳＥ２により収集されたエコートレインＥＴ２を概略的に示す図である。
【図１７】シーケンスＳＥ１〜ＳＥｍにより収集されたエコートレインＥＴ１〜ＥＴｍと、第１エコーＥ１１〜Ｅｍ１を逆フーリエ変換することにより得られたプロファイルＦ１〜Ｆｍを概略的に示す図である。
【図１８】エコートレインＥＴ１〜ＥＴｍの中から再収集する必要があるエコートレインを決定する方法の一例の説明図である。
【図１９】追加のスキャンの説明図である。
【発明を実施するための形態】
【００１０】
以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。
【００１１】
（１）第１の形態
図１は、本発明の第１の形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ。ＭＲ：Magnetic Resonance）１００は、マグネット２、テーブル３、受信コイル４などを有している。
【００１２】
マグネット２は、被検体１０が収容されるボア２１を有している。また、マグネット２は、超伝導コイル、勾配コイル、およびＲＦコイルなどが内蔵されている。
【００１３】
テーブル３は、被検体１０を支持するクレードル３ａを有している。クレードル３ａは、ボア２１内に移動できるように構成されている。クレードル３ａによって、被検体１０はボア２１に搬送される。
【００１４】
受信コイル４は、被検体１０の胴部に取り付けられている。受信コイル４は、被検体１０からの磁気共鳴信号を受信する。
【００１５】
ＭＲ装置１００は、更に、送信器５、勾配磁場電源６、制御部７、操作部８、および表示部９などを有している。
送信器５はＲＦコイルに電流を供給し、勾配磁場電源６は勾配コイルに電流を供給する。
【００１６】
制御部７は、表示部９に必要な情報を伝送したり、受信コイル４から受け取ったデータに基づいて画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。制御部７は、ピーク値検出手段７１および決定手段７２などを有している。
【００１７】
ピーク値検出手段７１は、エコートレインに含まれる複数のエコーのうちの第１エコーのピーク値を検出する。
決定手段７２は、第１エコーのピーク値に基づいて、再収集するエコートレインを決定する。
【００１８】
制御部７は、ピーク値検出手段７１および決定手段７２を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。尚、ピーク値検出手段７１は、指標を求める手段の一例に相当する。
【００１９】
操作部８は、オペレータにより操作され、種々の情報を制御部７に入力する。表示部９は種々の情報を表示する。
ＭＲ装置１００は、上記のように構成されている。
【００２０】
図２は第１の形態における撮影領域を概略的に示す図である。第１の形態では、肝臓を含む領域が、撮影領域Ｒとして設定されている。
【００２１】
図３は、撮影領域Ｒからデータを収集するために実行されるスキャンＳＣの説明図、図４は、スキャンＳＣにおいて実行されるシーケンスの説明図である。
スキャンＳＣでは、シーケンスＳＥｉ（ｉ＝１〜ｍ）が実行される。図４に示すように、シーケンスＳＥｉは、プリパレーション部Ａとデータ収集部Ｂとを有している。
【００２２】
プリパレーション部Ａは、スピンの動きを検出するための傾斜磁場ＭＰＧを有している。データ収集部Ｂは、撮影領域ＲからエコートレインＥＴｉを収集するためのシーケンス部である。第１の形態では、シーケンスＳＥ１〜ＳＥｍが実行されるので、ｍ個のエコートレインＥＴ１〜ＥＴｍが得られる。
【００２３】
エコートレインＥＴｉはｎ個のエコーＥｉ１〜Ｅｉｎを有している。第１の形態では、エコーＥｉ１〜Ｅｉｎのうち第１エコーＥｉ１および第２エコーＥｉ２は、画像再構成のデータとしては使用されず、破棄される。そして、第３エコーＥｉ３〜第ｎエコーＥｉｎが画像再構成のデータとして使用される。画像再構成のデータとして使用される第３エコーＥｉ３〜第ｎエコーＥｉｎの期間は、傾斜磁場ＧｓｓおよびＧｐｅのエンコード量は増減する。しかし、画像再構成のデータとして使用されない第１エコーＥｉ１〜第２エコーＥｉ２の期間は、傾斜磁場ＧｓｓおよびＧｐｅのエンコード量はゼロに設定されている。
図３には、シーケンスＳＥ１、ＳＥ２、ＳＥｊ、およびＳＥｍにより収集されるエコートレインＥＴ１、ＥＴ２、ＥＴｊ、およびＥＴｍが概略的に示されている。
【００２４】
図５は、シーケンスＳＥ１〜ＳＥｍにより収集されたエコーをｋ空間に埋めるときの様子を示す図である。
シーケンスＳＥ１を実行することによりエコートレインＥＴ１（第１エコーＥ１１〜第ｎエコーＥ１ｎ）が収集される。上記のように、第１エコーＥ１１および第２エコーＥ１２は画像再構成のデータとしては使用されず、第３エコーＥ１３〜第ｎエコーＥ１ｎが画像再構成のデータとして使用される。エコートレインＥＴ１の第３エコーＥ１３〜第ｎエコーＥ１ｎは、ｋｚ＝１のｋ空間のデータとして使用される。シーケンスＳＥ１を実行した後、次のシーケンスＳＥ２が実行される。
【００２５】
シーケンスＳＥ２を実行することによりエコートレインＥＴ２（第１エコーＥ２１〜第ｎエコーＥ２ｎ）が収集される。第１エコーＥ２１および第２エコーＥ２２は画像再構成のデータとしては使用されず、第３エコーＥ２３〜第ｎエコーＥ２ｎが画像再構成のデータとして使用される。エコートレインＥＴ２の第３エコーＥ２３〜第ｎエコーＥ２ｎは、ｋｚ＝２のｋ空間のデータとして使用される。
【００２６】
以下同様に、各シーケンスが実行され、第３エコー〜第ｎエコーが画像再構成のデータとして使用される。例えば、シーケンスＳＥｊを実行することによりエコートレインＥＴｊ（第１エコーＥｊ１〜第ｎエコーＥｊｎ）が収集され、第３エコーＥｊ３〜第ｎエコーＥｊｎが、ｋｚ＝ｊのｋ空間のデータとして使用される。そして、最後にシーケンスＳＥｍが実行される。
【００２７】
シーケンスＳＥｍを実行することによりエコートレインＥＴｍ（第１エコーＥｍ１〜第ｎエコーＥｍｎ）が収集される。第１エコーＥｍ１および第２エコーＥｍ２は画像再構成のデータとしては使用されず、第３エコーＥｍ３〜第ｎエコーＥｍｎが画像再構成のデータとして使用される。エコートレインＥＴｍの第３エコーＥｍ３〜第ｎエコーＥｍｎは、ｋｚ＝ｍのｋ空間のデータとして使用される。
【００２８】
このようにして、ｋｚ＝１〜ｍのｋ空間のデータが収集される。ｋｚ＝１〜ｍのｋ空間のデータを実空間のデータに変換することにより、撮影領域Ｒの拡散の情報を含む画像データを得ることができる。しかし、シーケンスＳＥｉの傾斜磁場ＭＰＧ（図４参照）のｂ値が大きい場合、エコートレインのエコーが心拍動などの体動の影響を受けて、画像データに、信号消失などのアーチファクトが現れやすくなる。このようなアーチファクトは画質を劣化させるので、アーチファクトはできるだけ低減する必要がある。そこで、第１の形態では、アーチファクトを低減するために、以下のようなフローに従ってスキャンを実行する（図６参照）。
【００２９】
図６は、第１の形態におけるスキャンを実行するときのフローを示す図である。
ステップＳＴ１では、上記のスキャンＳＣ（図３参照）を実行する。スキャンＳＣでは、最初に、シーケンスＳＥ１が実行される。シーケンスＳＥ１が実行されることにより、エコートレインＥＴ１が収集される。図７に、シーケンスＳＥ１により収集されたエコートレインＥＴ１を概略的に示す。エコートレインＥＴ１の第３エコーＥ１３〜第ｎエコーＥ１ｎは、ｋｚ＝１に配置される。また、第１の形態では、ピーク値検出手段７１（図１参照）が、画像再構成に使用されない第１エコーＥ１１のピーク値Ｐ１１を検出する。
シーケンスＳＥ１を実行した後、次のシーケンスＳＥ２が実行される。
【００３０】
シーケンスＳＥ２が実行されることにより、エコートレインＥＴ２が収集される。図８に、シーケンスＳＥ２により収集されたエコートレインＥＴ２を概略的に示す。エコートレインＥＴ２の第３エコーＥ２３〜第ｎエコーＥ２ｎは、ｋｚ＝２のｋ空間のデータとして使用される。また、ピーク値検出手段７１は、画像再構成に使用されない第１エコーＥ２１のピーク値Ｐ２１を検出する。
【００３１】
以下同様に、各シーケンスを順に実行し、各シーケンスで得られた第１エコーのピーク値を求める。図９に、シーケンスＳＥ１〜ＳＥｍにより収集されたエコートレインＥＴ１〜ＥＴｍと、各エコートレインの第１エコーのピーク値Ｐ１１〜Ｐｍ１を示す。また、図９の左下には、第１エコーＥ１１〜Ｅｍ１のピーク値Ｐ１１〜Ｐｍ１の大きさの違いを概略的に表すグラフも示されている。
スキャンＳＣを実行した後、ステップＳＴ２に進む。
【００３２】
ステップＳＴ２では、決定手段７２（図１参照）が、エコートレインＥＴ１〜ＥＴｍの中から、再収集するエコートレインを決定する。
心拍などの体動はエコートレインのエコーに影響を与えることがあり、体動の影響を受けているエコーを画像再構成のエコーとして使用すると、画像の信号消失などのアーチファクトの原因となる。そこで、第１の形態では、エコートレインＥＴ１〜ＥＴｍの中に、体動の影響を大きく受けている可能性の高いエコートレインが含まれている場合は、エコートレインの再収集を行う。決定手段７２は、エコートレインＥＴ１〜ＥＴｍの中から、再収集する必要があるエコートレインを決定する。以下に、この決定方法について説明する。
【００３３】
図１０は、エコートレインＥＴ１〜ＥＴｍの中から再収集する必要があるエコートレインを決定する方法の一例の説明図である。
決定手段７２は、先ず、エコートレインＥＴ１〜ＥＴｍの第１エコーＥ１１〜Ｅｍ１のピーク値Ｐ１１〜Ｐｍ１に基づいて、体動の影響を大きく受けている第１エコーを特定する。エコーのピーク値は体動の影響により小さくなる傾向があるので、エコーのピーク値を求めることによって、体動がエコーにどの程度の影響を与えているかを知ることができる。上記のように、エコーのピーク値は体動の影響により小さくなる傾向があるので、ピーク値の小さい第１エコーを特定することにより、体動の影響を大きく受けている第１エコーを特定することができる。第１の形態では、決定手段７２は、第１エコーＥ１１〜Ｅｍ１のピーク値が大きいか小さいかを判断するための閾値Ｔを求め、第１エコーＥ１１〜Ｅｍ１の中で、ピーク値が閾値Ｔより小さい第１エコーを、体動の影響を大きく受けている第１エコーと判断する。閾値Ｔを求める方法の一例としては、ｍ個のピーク値Ｐ１１〜Ｐｍ１の中から、ピーク値の大きい順にｋ（＜ｍ）個のピーク値を抽出し、抽出されたピーク値の平均値ＭＳを閾値Ｔとする方法がある。また、以下の式を用いて閾値Ｔを求めてもよい。
Ｔ＝ＭＳ＊Ｋ    ・・・（１）
ここで、ＭＳ：抽出されたピーク値の平均値
Ｋ：係数（例えば、Ｋ＝０．９）
【００３４】
図１０には、ピーク値が閾値Ｔより小さい第１エコーとして、第１エコーＥ２１およびＥｍ１が示されている。尚、第１エコーＥ２１およびＥｍ１の他にも、ピーク値が閾値Ｔより小さい第１エコーは存在するが、以下では、説明の便宜上、第１のエコーＥ２１のピーク値Ｐ２１と、第１のエコーＥｍ１のピーク値Ｐｍ１のみが、閾値Ｔより小さいとする。したがって、第１エコーＥ２１およびＥｍ１が、体動の影響を大きく受けている第１エコーとして特定される。
【００３５】
エコートレインの中の第１エコーが体動の影響を大きく受けている場合、第１エコーに続いて収集される第２エコー〜第ｎエコーも体動の影響を大きく受けている可能性が高いと考えられる。体動の影響を大きく受けているエコーを画像再構成のエコーとして使用すると、画像の信号消失などのアーチファクトとなって現れるので、画像劣化の原因となる。そこで、第１エコーのピーク値が閾値Ｔより小さい場合は、エコートレインを再収集すると決定する。第１の形態では、第１エコーＥ２１およびＥｍ１のピーク値が閾値Ｔより小さいので、決定手段７２は、第１エコーＥ２１含むエコートレインＥＴ２と、第１エコーＥｍ１含むエコートレインＥＴｍとを再収集すると決定する。
【００３６】
尚、図１０では、エコートレインＥＴ１〜ＥＴｍの第１エコーＥ１１〜Ｅｍ１のピーク値Ｐ１１〜Ｐｍ１に基づいて、エコートレインＥＴ１〜ＥＴｍの中から、再収集するエコートレインを決定している。しかし、第１エコーＥ１１〜Ｅｍ１のピーク値の代わりに第２エコーＥ１２〜Ｅｍ２のピーク値に基づいて、エコートレインＥＴ１〜ＥＴｍの中から、再収集するエコートレインを決定してもよい（図１１参照）。
【００３７】
図１１は、第２エコーＥ１２〜Ｅｍ２のピーク値に基づいて、エコートレインＥＴ１〜ＥＴｍの中から、再収集するエコートレインを決定するときの説明図である。
第２エコーＥ１２〜Ｅｍ２は、第１エコーＥ１１〜Ｅｍ１と同様に、体動の影響を受けるとピーク値が小さくなる傾向があるので、第２エコーＥ１２〜Ｅｍ２のピーク値Ｐ１２〜Ｐｍ２を比較しても、第３エコー〜第ｎエコーが体動の影響を大きく受けているか否かを判断することができる。したがって、第２エコーＥ１２〜Ｅｍ２のピーク値Ｐ１２〜Ｐｍ２に対して閾値Ｔを設定し、第２エコーＥ１２〜Ｅｍ２の中から、ピーク値が閾値Ｔより小さい第２エコーを求めることによって、エコートレインＥＴ１〜ＥＴｍの中から、再収集するエコートレインを決定してもよい。
【００３８】
尚、エコートレインＥＴ１〜ＥＴｍには第１エコーおよび第２エコーの他に、画像再構成に使用される第３エコー〜第ｎエコーも含まれている。したがって、第１エコーや第２エコーの代わりに、画像再構成に使用される第３エコー〜第ｎエコーを用いて、エコートレインＥＴ１〜ＥＴｍの中から、再収集するエコートレインを決定することも考えられる。しかし、第３エコー〜第ｎエコーを収集する場合、傾斜磁場Ｇｓｓのエンコード量や傾斜磁場Ｇｐｅのエンコード量を増減させる必要があるので（図４参照）、エコーのピーク値の低下が、体動が原因で生じているのか、エンコード量の違いが原因で生じているのか判断できない。これに対し、第１エコーおよび第２エコーを収集する場合、傾斜磁場Ｇｓｓのエンコード量や傾斜磁場Ｇｐｅのエンコード量は増減させないので、エンコード量は同じ値である（図４では、第１エコーおよび第２エコーを収集するときのエンコード量はゼロに設定されている）。したがって、第１エコーおよび第２エコーの場合は、エコーのピーク値の低下は体動が原因で生じていると考えることができる。このような理由から、本形態では、第１のエコー又は第２のエコーを用いて、再収集するエコートレインを決定している。尚、第１の形態では、第１エコーおよび第２エコーを収集するときの傾斜磁場Ｇｓｓのエンコード量はゼロであるが、エンコード量を同じ値にすることができるのであれば、エンコード量はゼロより大きい値であってもよい。同様に、傾斜磁場Ｇｐｅのエンコード量もゼロにより大きい値に設定してもよい。
エコートレインＥＴ１〜ＥＴｍの中から、再収集するエコートレインＥＴ２およびＥＴｍを決定した後、ステップＳＴ３に進む。
【００３９】
ステップＳＴ３では、エコートレインＥＴ２およびＥＴｍを再収集するための追加のスキャンを実行する（図１２参照）。
図１２は追加のスキャンの説明図である。
追加のスキャンＳＣｐでは、エコートレインＥＴ２を再収集するためのシーケンスＳＥ２と、エコートレインＥＴｍを再収集するためのシーケンスＳＥｍが実行される。
【００４０】
最初に、１回目のシーケンスＳＥ２が実行される。ピーク値検出手段７１は、１回目のシーケンスＳＥ２により収集された第１エコーＥ２１のピーク値ａ１を検出する。そして、決定手段７２はピーク値ａ１と閾値Ｔとを比較する。図１２では、ピーク値ａ１は、閾値Ｔよりも小さい値になっている。上述したように、第１エコーのピーク値が閾値Ｔよりも小さい場合は、体動の影響を大きく受けていると考えられる。したがって、１回目のシーケンスＳＥ２により収集された第３エコーＥ２３〜第ｎエコーＥ２ｎは、ｋｚ＝２におけるｋ空間のデータとして採用できない。そこで、追加のスキャンＳＣｐでは、第１エコーＥ２１のピーク値が閾値Ｔより大きくなるまで、シーケンスＳＥ２を繰り返し実行する。図１２では、１回目〜ｖ−１回目のシーケンスＳＥ２で収集された第１エコーＥ２１のピーク値ａ１〜ａｖ−１は閾値Ｔより小さいが、ｖ回目のシーケンスＳＥ２で収集された第１エコーＥ２１のピーク値ａｖは閾値Ｔより大きい。したがって、ｖ回目のシーケンスＳＥ２で収集されたエコートレインＥＴ２の第３エコーＥ２３〜第ｎエコーＥ２ｎが、ｋｚ＝２におけるｋ空間のデータとして採用される。
【００４１】
ｖ回目のシーケンスＳＥ２を実行した後、エコートレインＥＴｍを再収集するためのシーケンスＳＥｍが実行される。シーケンスＳＥｍも、第１エコーＥｍ１のピーク値が閾値Ｔより大きくなるまで繰り返し実行される。シーケンスＳＥｍについては、ｗ回目で第１エコーＥｍ１のピーク値ｂｗが閾値Ｔより大きくなったとする。したがって、ｗ回目のシーケンスＳＥｍで収集されたエコートレインＥＴｍの第３エコーＥｍ３〜第ｎエコーＥｍｎが、ｋｚ＝ｍにおけるｋ空間のデータとして採用される。
【００４２】
したがって、ステップＳＴ１のスキャンＳＣと、ステップＳＴ３の追加のスキャンＳＣｐによって、画像再構成に必要なｋ空間のデータＤ′を取得することができる。
【００４３】
尚、閾値Ｔの値が大きすぎる場合、追加のスキャンＳＣｐで得られた第１エコーのピーク値が閾値Ｔを超えにくくなるので、追加のスキャンＳＣｐのスキャン時間が長くなることがある。そこで、追加のスキャンＳＣｐの直前、あるいは追加のスキャンＳＣｐの最中に、オペレータが閾値Ｔの値を変更できるようにしてもよい。これにより、オペレータは、追加のスキャンＳＣｐの直前、あるいは追加のスキャンＳＣｐの最中に、閾値Ｔの値を小さくすることができる。閾値Ｔの値を小さくすることにより、追加のスキャンＳＣｐで得られた第１エコーのピーク値が閾値Ｔを超えやすくなるので、追加のスキャンＳＣｐのスキャン時間を短くすることができる。
追加のスキャンＳＣｐを実行した後、ステップＳＴ４に進む。
【００４４】
ステップＳＴ４では、図１２に示すｋ空間のデータＤ′を用いて画像再構成が行われ、フローが終了する。
【００４５】
第１の形態では、スキャンＳＣを実行することにより、エコートレインＥＴ１〜ＥＴｍを収集し、エコートレインＥＴ１〜ＥＴｍの第１エコーＥ１１〜Ｅｍ１のピーク値Ｐ１１〜Ｐｍ１を求める（図９参照）。そして、第１エコーＥ１１〜Ｅｍ１のピーク値Ｐ１１〜Ｐｍ１と閾値Ｔとを比較する（図１０参照）。第１エコーのピーク値が閾値Ｔより低い場合は、第１エコーに続いて収集される第２エコー〜第ｎエコーも体動の影響を大きく受けている可能性が高いと考えられる。そこで、第１の形態では、第１エコーのピーク値が閾値Ｔより低い場合は、第１エコーのピーク値が閾値Ｔより大きくなるまでエコートレインを再収集するための追加のスキャンＳＣｐを実行する（図１２参照）。したがって、体動の影響をあまり受けていないエコーを用いて画像再構成が行われるので、信号消失などのアーチファクトが低減された画像データを得ることができる。
【００４６】
尚、第１の形態では、体動がエコーにどの程度の影響を与えているかを知るために、エコーのピーク値を求めている。しかし、エコーのピーク値の代わりに、エコーの面積を用いてもよい。体動の影響を大きく受けたエコーは信号値が小さくなる傾向があるので、それに伴いエコーの面積も小さくなる傾向がある。したがって、エコーのピーク値の代わりにエコーの面積を求めても、体動がエコーにどの程度の影響を与えているかを知ることができる。図１３は、エコーの面積を求める方法の一例の説明図である。例えば、エコーの包絡線により囲まれた面積Ａをエコーの面積としてもよいし、エコーの絶対値の波形の面積Ｂをエコーの面積としてもよい。したがって、閾値などを用いて、第１エコーＥ１１〜Ｅｍ１の中から面積Ａ（又は面積Ｂ）が小さい第１エコーを特定することにより、再収集するエコートレインを決定することができる。
【００４７】
（２）第２の形態
図１４は、第２の形態のＭＲ装置の概略図である。
第２の形態のＭＲ装置２００は、第１の形態のＭＲ装置１００（図１参照）と比較すると、制御部７が異なっているが、その他の構成は同じである。したがって、第２の形態のＭＲ装置２００については、主に、制御部７について説明する。
【００４８】
第２の形態では、制御部７は、変換手段７１１、検出手段７１２、および決定手段７１３などを有している。
変換手段７１１は、エコーを実空間のプロファイルに変換する。
検出手段７１２は、変換手段７１１により得られたプロファイルの所定の位置における信号値を検出する。
決定手段７１３は、検出手段７１２により検出された信号値に基づいて、再収集するエコートレインを決定する。
【００４９】
制御部７は、変換手段７１１、検出手段７１２、および決定手段７１３を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。尚、変換手段７１１および検出手段７１２を合わせたものが、指標を求める手段の一例に相当する。
【００５０】
次に、第２の形態のＭＲ装置の動作について、図１の形態と同様に、図６のフローを参照しながら説明する。
【００５１】
ステップＳＴ１では、スキャンＳＣを実行する。図１５にスキャンＳＣを概略的に示す。スキャンＳＣでは、最初に、シーケンスＳＥ１が実行される。シーケンスＳＥ１が実行されることにより、エコートレインＥＴ１が収集される。エコートレインＥＴ１の第３エコーＥ１３〜第ｎエコーＥ１ｎは、ｋｚ＝１に配置される。また、変換手段７１１（図１４参照）は、画像再構成に使用されない第１エコーＥ１１をｘ方向（周波数エンコード方向）に逆フーリエ変換し、実空間上のプロファイルＦ１を求める。プロファイルＦ１は、ｘ方向（周波数エンコード方向）の信号値の強度変化を表している。
シーケンスＳＥ１を実行した後、次のシーケンスＳＥ２が実行される。
【００５２】
シーケンスＳＥ２が実行されることにより、エコートレインＥＴ２が収集される。図１６に、シーケンスＳＥ２により収集されたエコートレインＥＴ２を概略的に示す。エコートレインＥＴ２の第３エコーＥ２３〜第ｎエコーＥ２ｎは、ｋｚ＝２のｋ空間のデータとして使用される。また、変換手段７１１は、画像再構成に使用されない第１エコーＥ２１をｘ方向に逆フーリエ変換し、実空間上のプロファイルＦ２を求める。
【００５３】
以下同様に、各シーケンスを順に実行し、各シーケンスで得られた第１エコーを逆フーリエ変換し、実空間上のプロファイルを求める。図１７に、シーケンスＳＥ１〜ＳＥｍにより収集されたエコートレインＥＴ１〜ＥＴｍと、第１エコーＥ１１〜Ｅｍ１を逆フーリエ変換することにより得られたプロファイルＦ１〜Ｆｍを概略的に示す。
スキャンＳＣを実行した後、ステップＳＴ２に進む。
【００５４】
ステップＳＴ２では、決定手段７１３（図１４参照）が、エコートレインＥＴ１〜ＥＴｍの中から、再収集するエコートレインを決定する。
図１８は、エコートレインＥＴ１〜ＥＴｍの中から再収集する必要があるエコートレインを決定する方法の一例の説明図である。
【００５５】
先ず、肝臓のうちの心臓の近くに位置する部位Ａを考える。ここで、部位Ａのｘ方向の位置をｘ＝ｘ１とする。検出手段７１２（図１４参照）は、プロファイルＦ１〜Ｆｍの各々について、ｘ＝ｘ１における信号値Ｓ（＝Ｓ１〜Ｓｍ）を検出する。エコーが心拍動の影響を受けた場合、心臓の近くに位置する部位Ａの信号が消失する現象が生じることがある。部位Ａの信号が消失した場合、ｘ＝ｘ１における信号値Ｓは小さくなる。したがって、プロファイルＦ１〜Ｆｍの中で、信号値Ｓの小さいプロファイルが、部位Ａの信号が消失している可能性が高いプロファイルと考えられる。第２の形態では、決定手段７１３は、信号値Ｓ１〜Ｓｍが大きいか小さいかを判断するための閾値Ｔを求め、プロファイルＦ１〜Ｆｍの中で、ｘ＝ｘ１における信号値Ｓが閾値Ｔより小さいプロファイルを、部位Ａの信号が消失している可能性が高いプロファイルと判断する。閾値Ｔを求める方法の一例としては、ｍ個の信号値Ｓ１〜Ｓｍの中から、信号値の大きい順にｋ（＜ｍ）個の信号値を抽出し、抽出された信号値の平均値ＭＳを閾値Ｔとする方法がある。また、以下の式を用いて閾値Ｔを求めてもよい。
Ｔ＝ＭＳ＊Ｋ    ・・・（２）
ここで、ＭＳ：抽出された信号値の平均値
Ｋ：係数（例えば、Ｋ＝０．９）
【００５６】
図１８には、閾値Ｔより小さい信号値として、信号値Ｓ２およびＳｍが示されている。尚、信号値Ｓ２およびＳｍの他にも、閾値Ｔより小さい信号値は存在するが、以下では、説明の便宜上、信号値Ｓ２およびＳｍのみが、閾値Ｔより小さいとする。したがって、２つのプロファイルＦ２およびＦｍが、部位Ａの信号が消失している可能性が高いプロファイルとして特定される。
【００５７】
プロファイルＦ２は、エコートレインＥＴ２の第１エコーＥ２１を逆フーリエ変換することにより得られており、プロファイルＦｍは、エコートレインＥＴｍの第１エコーＥｍ１を逆フーリエ変換することにより得られている。したがって、第１エコーＥ２１およびＥｍ１が、心拍動の影響を大きく受けている第１エコーとして特定される。
【００５８】
エコートレインの中の第１エコーが心拍動の影響を大きく受けている場合、第１エコーに続いて収集される第２エコー〜第ｎエコーも心拍動の影響を大きく受けている可能性が高いと考えられる。心拍動の影響を受けているエコーを画像再構成のエコーとして使用すると、画像の信号消失などのアーチファクトとなって現れるので、画像劣化の原因となる。そこで、決定手段７１３は、第１エコーＥ２１含むエコートレインＥＴ２と、第１エコーＥｍ１含むエコートレインＥＴｍとを再収集すると決定する。エコートレインＥＴ２およびＥＴｍを再収集すると決定した後、ステップＳＴ３に進む。
【００５９】
ステップＳＴ３では、エコートレインＥＴ２およびＥＴｍを再収集するための追加のスキャンを実行する（図１９参照）。
【００６０】
図１９は追加のスキャンの説明図である。
追加のスキャンＳＣｐでは、エコートレインＥＴ２を再収集するためのシーケンスＳＥ２と、エコートレインＥＴｍを再収集するためのシーケンスＳＥｍが実行される。
最初に、１回目のシーケンスＳＥ２が実行される。変換手段７１１は、１回目のシーケンスＳＥ２により収集された第１エコーＥ２１を逆フーリエ変換し、プロファイルＦ２＿１を求める。そして、検出手段７１２は、プロファイルＦ２＿１のｘ＝ｘ１における信号値Ｓ２＿１を検出する。次に、決定手段７１３は、信号値Ｓ２＿１と閾値Ｔとを比較する。図１９では、信号値Ｓ２＿１は、閾値Ｔよりも小さい値になっている。上述したように、プロファイルのｘ＝ｘ１における信号値が閾値Ｔよりも小さい場合は、エコーは心拍動の影響を大きく受けていると考えられる。したがって、１回目のシーケンスＳＥ２により収集された第３エコーＥ２３〜第ｎエコーＥ２ｎは、ｋｚ＝２におけるｋ空間のデータとして採用できない。そこで、シーケンスＳＥ２を実行する場合、プロファイルのｘ＝ｘ１における信号値Ｓが閾値Ｔより大きくなるまで、シーケンスＳＥ２を繰り返し実行する。図１９の場合、１回目〜ｖ−１回目のシーケンスＳＥ２では、プロファイルのｘ＝ｘ１における信号値Ｓは閾値Ｔより小さいが、ｖ回目のシーケンスＳＥ２では、プロファイルのｘ＝ｘ１における信号値Ｓ（＝Ｓ２＿ｖ）が閾値Ｔより大きくなったとする。したがって、ｖ回目のシーケンスＳＥ２で収集されたエコートレインＥＴ２の第３エコーＥ２３〜第ｎエコーＥ２ｎが、ｋｚ＝２におけるｋ空間のデータとして採用される。
【００６１】
ｖ回目のシーケンスＳＥ２を実行した後、エコートレインＥＴｍを再収集するためのシーケンスＳＥｍが実行される。シーケンスＳＥｍも、プロファイルのｘ＝ｘ１における信号値Ｓが閾値Ｔより大きくなるまで繰り返し実行される。シーケンスＳＥｍについては、ｗ回目でプロファイルのｘ＝ｘ１における信号値Ｓ（＝Ｓｍ＿ｗ）が閾値Ｔより大きくなったとする。したがって、ｗ回目のシーケンスＳＥｍで収集されたエコートレインＥＴｍの第３エコーＥｍ３〜第ｎエコーＥｍｎが、ｋｚ＝ｍにおけるｋ空間のデータとして採用される。
【００６２】
したがって、ステップＳＴ１のスキャンＳＣと、ステップＳＴ３の追加のスキャンＳＣｐによって、画像再構成に必要なｋ空間のデータＤ′を取得することができる。
【００６３】
尚、閾値Ｔの値が大きすぎる場合、追加のスキャンＳＣｐで得られたプロファイルのｘ＝ｘ１における信号値Ｓが閾値Ｔを超えにくくなるので、追加のスキャンＳＣｐのスキャン時間が長くなることがある。そこで、追加のスキャンＳＣｐの直前、あるいは追加のスキャンＳＣｐの最中に、オペレータが閾値Ｔの値を変更できるようにしてもよい。これにより、オペレータは、追加のスキャンＳＣｐの直前、あるいは追加のスキャンＳＣｐの最中に、閾値Ｔの値を小さくすることができる。閾値Ｔの値を小さくすることにより、追加のスキャンＳＣｐで得られたプロファイルのｘ＝ｘ１における信号値Ｓが閾値Ｔを超えやすくなるので、追加のスキャンＳＣｐのスキャン時間を短くすることができる。
追加のスキャンＳＣｐを実行した後、ステップＳＴ４に進む。
【００６４】
ステップＳＴ４では、図１９に示すｋ空間のデータＤ′を用いて画像再構成が行われ、フローが終了する。
【００６５】
第２の形態では、スキャンＳＣを実行することにより、エコートレインＥＴ１〜ＥＴｍを収集し、エコートレインＥＴ１〜ＥＴｍの第１エコーＥ１１〜Ｅｍ１を実空間のプロファイルＦ１〜Ｆｍに変換する（図１７参照）。そして、肝臓の位置ｘ＝ｘ１における信号値Ｓ１〜Ｓｍと閾値Ｔとを比較する（図１８参照）。ｘ＝ｘ１における信号値が閾値Ｔより低い場合は、第１エコーに続いて収集される第２エコー〜第ｎエコーも心拍動の影響を大きく受けている可能性が高いと考えられる。そこで、第２の形態では、ｘ＝ｘ１における信号値が閾値Ｔより低い場合は、ｘ＝ｘ１における信号値が閾値Ｔより大きくなるまでエコートレインを再収集するための追加のスキャンＳＣｐを実行する（図１９参照）。したがって、心拍動の影響をあまり受けていないエコーを用いて画像再構成が行われるので、信号消失などのアーチファクトが低減された画像データを得ることができる。
【００６６】
尚、第１および第２の形態では、第１エコーおよび第２エコーを、画像再構成に使用されないエコーとしている。しかし、画像再構成に使用されないエコーは、第１エコーや第２エコーに限定されることはなく、第１エコー〜第ｎエコーのうちのどのエコーであってもよい。
【符号の説明】
【００６７】
２  マグネット
３  テーブル
３ａ  クレードル
４  受信コイル
５  送信器
６  勾配磁場電源
７  制御部
８  操作部
９  表示部
１０  被検体
２１  ボア
７１  ピーク値検出手段
７２、７１３  決定手段
１００、２００  ＭＲ装置
７１１  変換手段
７１２  検出手段
【図１】
2014147466.tif 000003
【図２】
2014147466.tif 000004
【図３】
2014147466.tif 000005
【図４】
2014147466.tif 000006
【図５】
2014147466.tif 000007
【図６】
2014147466.tif 000008
【図７】
2014147466.tif 000009
【図８】
2014147466.tif 000010
【図９】
2014147466.tif 000011
【図１０】
2014147466.tif 000012
【図１１】
2014147466.tif 000013
【図１２】
2014147466.tif 000014
【図１３】
2014147466.tif 000015
【図１４】
2014147466.tif 000016
【図１５】
2014147466.tif 000017
【図１６】
2014147466.tif 000018
【図１７】
2014147466.tif 000019
【図１８】
2014147466.tif 000020
【図１９】
2014147466.tif 000021