(57)【要約】      （修正有）【課題】ＡＤＣ（拡散係数）の算出精度を向上させることのできる磁気共鳴装置を提供する。【解決手段】ＭＰＧ（Ｍｏｔｉｏｎ  Ｐｒｏｂｉｎｇ  Ｇｒａｄｉｅｎｔ）を有するプリパレーションシーケンスと撮影部位からＭＲ信号を収集するためのイメージングシーケンスとを含む拡散強調スキャンを実行するＭＲ装置である。ＭＲ装置は、ＭＰＧのｂ値を変更し、ＭＰＧのｂ値が異なる拡散強調スキャンＤＷ１およびＤＷ２を実行するとともに、撮影部位のＴ１値を計測するためのＴ１値計測用スキャンＳＣ１、および撮影部位のＴ２値を計測するためのＴ２値計測用スキャンＳＣ２を実行するスキャン手段と、時間ＴＤＰ、Ｔｒｅｃ、ＭＲ信号の信号値Ｓ１、Ｓ２、Ｔ１値、およびＴ２値に基づいて、撮影部位のＡＤＣを算出する算出手段とを有する。【選択図】図６
2015123305.tif 000002
【特許請求の範囲】
【請求項１】
ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンを実行する磁気共鳴装置であって、
前記ＭＰＧの強さを表すｂ値を変更し、前記ＭＰＧのｂ値が異なる複数の拡散強調スキャンを実行するとともに、前記撮影部位のＴ１値を計測するためのＴ１値計測用スキャンを実行するスキャン手段と、
前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記複数の拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出手段と、
を有する、磁気共鳴装置。
【請求項２】
ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンと、ＭＰＧを有していない第３のシーケンスと前記撮影部位からＭＲ信号を収集するための第４のシーケンスとを含む他のスキャンとを実行する磁気共鳴装置であって、
前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、前記他のスキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出手段と、
を有する、磁気共鳴装置。
【請求項３】
前記スキャン手段は、
前記撮影部位のＴ２値を計測するためのＴ２値計測用スキャンを実行し、
前記算出手段は、
前記Ｔ２値計測用スキャンにより計測されたＴ２値も考慮して、前記見かけの拡散係数を算出する、請求項１又は２に記載の磁気共鳴装置。
【請求項４】
前記第１のシーケンスは、縦磁化をフリップさせるための第１のＲＦパルスと、横磁化を縦磁化に戻すための第２のＲＦパルスとを有し、
前記算出手段は、
前記第１のＲＦパルスと前記第２のＲＦパルスとの間の第２の時間も考慮して、前記見かけの拡散係数を算出する、請求項３に記載の磁気共鳴装置。
【請求項５】
前記第２のシーケンスは、前記撮影部位からＭＲ信号を収集するための複数のＲＦパルスを有し、
前記第１の時間は、前記第２のＲＦパルスと、前記複数のＲＦパルスのうちの一番最初に送信されるＲＦパルスとの間の時間である、請求項１〜４のうちのいずれか一項に記載の磁気共鳴装置。
【請求項６】
前記複数のＲＦパルスのフリップ角は次第に小さくなるように設定されている、請求項５に記載の磁気共鳴装置。
【請求項７】
前記第２のシーケンスは、前記複数のＲＦパルスの前に、ＭＲ信号を定常状態に移行させるためのＲＦパルスを有する、請求項５又は６に記載の磁気共鳴装置。
【請求項８】
前記第１のシーケンスは、横磁化を消去するためのキラーパルスを有する、請求項１〜７のうちのいずれか一項に記載の磁気共鳴装置。
【請求項９】
前記拡散強調スキャンは、第１のシーケンスと第２のシーケンスとの間に、脂肪を抑制するための第３のシーケンスを有する、請求項１〜８のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１０】
ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンと、前記撮影部位のＴ１値を計測するためのＴ１値計測用スキャンとを実行する磁気共鳴装置に適用されるプログラムであって、
前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記ＭＰＧの強さを表すｂ値が異なる複数の拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出処理を計算機に実行させるためのプログラム。
【請求項１１】
ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンと、ＭＰＧを有していない第３のシーケンスと撮影部位からＭＲ信号を収集するための第４のシーケンスとを含む他のスキャンとを実行する磁気共鳴装置に適用されるプログラムであって、
前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、前記他のスキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出処理を計算機に実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、拡散強調スキャンを実行する磁気共鳴装置、および磁気共鳴装置に適用されるプログラムに関する。
【背景技術】
【０００２】
拡散の情報を取得する方法としてＭＰＧ（Motion Probing Gradient）を用いた方法が知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１２−１５７６８７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ＭＰＧを用いて拡散の情報を得る方法の一例として、ＭＰＧをプリパルスとして印加したDiffusion Preparation（ＤＰ）と呼ばれる手法がある。この手法は、例えばＦＩＥＳＴＡ（Fast Imaging Employing Steady State Acquisition）でデータ収集する場合に適用されており、見かけの拡散係数ＡＤＣ（Apparent Diffusion Coefficient）などの拡散情報を比較的容易に得ることができる。しかし、ＤＰ法では、Ｔ１回復の影響を受けて、ＡＤＣの算出精度が悪くなることがある。したがって、ＡＤＣの算出精度を向上させることが望まれている。
【課題を解決するための手段】
【０００５】
本発明の第１の観点は、ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンを実行する磁気共鳴装置であって、
前記ＭＰＧの強さを表すｂ値を変更し、前記ＭＰＧのｂ値が異なる複数の拡散強調スキャンを実行するとともに、前記撮影部位のＴ１値を計測するためのＴ１値計測用スキャンを実行するスキャン手段と、
前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記複数の拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出手段と、
を有する、磁気共鳴装置である。
【０００６】
本発明の第２の観点は、ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンと、ＭＰＧを有していない第３のシーケンスと前記撮影部位からＭＲ信号を収集するための第４のシーケンスとを含む他のスキャンとを実行する磁気共鳴装置であって、
前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、前記他のスキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出手段と、
を有する、磁気共鳴装置である。
【０００７】
本発明の第３の観点は、ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンと、前記撮影部位のＴ１値を計測するためのＴ１値計測用スキャンとを実行する磁気共鳴装置に適用されるプログラムであって、
前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記ＭＰＧの強さを表すｂ値が異なる複数の拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出処理を計算機に実行させるためのプログラムである。
【０００８】
本発明の第４の観点は、ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンと、ＭＰＧを有していない第３のシーケンスと撮影部位からＭＲ信号を収集するための第４のシーケンスとを含む他のスキャンとを実行する磁気共鳴装置に適用されるプログラムであって、
前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、前記他のスキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出処理を計算機に実行させるためのプログラムである。
【発明の効果】
【０００９】
Ｔ１値を考慮してＡＤＣを算出するので、信頼性の高いＡＤＣが得られる。
【図面の簡単な説明】
【００１０】
【図１】本発明の一形態の磁気共鳴装置の概略図である。
【図２】ＡＤＣデータを取得するために実行されるスキャンの一例を示す図である。
【図３】拡散強調スキャンＤＷ１の説明図である。
【図４】拡散強調スキャンＤＷ２の説明図である。
【図５】本形態で実行されるスキャンの説明図である。
【図６】ＡＤＣデータを求めるときのフローを示す図である。
【図７】撮影部位のＴ１値を表すＴ１値データＥ１を概略的に示す図である。
【図８】撮影部位のＴ２値を表すＴ２値データＥ２を概略的に示す図である。
【図９】拡散強調スキャンＤＷ１により得られた撮影部位の信号値を表す信号値データＥ３を概略的に示す図である。
【図１０】拡散強調スキャンＤＷ２により得られた撮影部位の信号値を表す信号値データＥ４を概略的に示す図である。
【図１１】ＡＤＣの算出方法を概略的に示す図である。
【図１２】時間Ｔｒｅｃの他の例を示す図である。
【図１３】ＲＦパルスＸ１〜ＸＮのフリップ角α１〜αＮが次第に小さくなるように設定されたイメージングシーケンスＣを示す図である。
【図１４】式（１４）を用いてＡＤＣ値を計算する場合に実行されるスキャンの一例を示す図である。
【図１５】ｂ値が異なるｎ個の拡散強調スキャンＤＷ１〜ＤＷｎ（ｂ＝ｂ１、ｂ２、ｂ３、・・・ｂｎ）を実行する例を示す図である。
【発明を実施するための形態】
【００１１】
以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。
【００１２】
図１は、本発明の一形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ）１００は、マグネット２、テーブル３、受信コイル４などを有している。
【００１３】
マグネット２は、被検体１１が収容されるボア２１を有している。また、マグネット２には、超伝導コイル、勾配コイル、およびＲＦコイルなどが内蔵されている。
【００１４】
テーブル３は、被検体１１を支持するクレードル３ａを有している。クレードル３ａは、ボア２１内に移動できるように構成されている。クレードル３ａによって、被検体１１はボア２１に搬送される。
受信コイル４は、被検体１１からの磁気共鳴信号を受信する。
【００１５】
ＭＲ装置１００は、更に、送信器５、勾配磁場電源６、受信器７、制御部８、操作部９、および表示部１０などを有している。
【００１６】
送信器５はＲＦコイルに電流を供給し、勾配磁場電源６は勾配コイルに電流を供給する。受信器７は、受信コイル４から受け取った信号に対して、検波などの信号処理を実行する。尚、マグネット２、受信コイル４、送信器５、勾配磁場電源６，受信器７を合わせたものがスキャン手段に相当する。
【００１７】
制御部８は、表示部１０に必要な情報を伝送したり、受信器７から受け取ったデータに基づいて画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。制御部８は、データ処理手段８１およびＡＤＣ算出手段８２などを有している。
【００１８】
データ処理手段８１は、スキャンにより収集されたデータに基づいて、Ｔ１値、Ｔ２値、ＭＲ信号の信号値などを求める。
ＡＤＣ算出手段８２は、後述する式に基づいて、ＡＤＣを算出する。
【００１９】
制御部８は、データ処理手段８１およびＡＤＣ算出手段８２を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。
【００２０】
操作部９は、オペレータにより操作され、種々の情報を制御部８に入力する。表示部１０は種々の情報を表示する。
【００２１】
ＭＲ装置１００は、上記のように構成されている。本形態では、ＭＲ装置１００を用いてＡＤＣデータを取得する。以下に、ＡＤＣデータを取得するために実行されるスキャンについて説明する。
【００２２】
図２はＡＤＣデータを取得するために実行されるスキャンの一例を示す図である。
図２には、拡散強調データを取得するための２つの拡散強調スキャンＤＷ１およびＤＷ２が示されている。拡散強調スキャンＤＷ１およびＤＷ２は、ＭＰＧをプリパルスとして印加したDiffusion Preparation（ＤＰ）と呼ばれる手法で拡散情報を取得するスキャンである。以下、拡散強調スキャンＤＷ１およびＤＷ２について順に説明する。
【００２３】
図３は拡散強調スキャンＤＷ１の説明図である。
拡散強調スキャンＤＷ１では、シーケンス部ＳＥ１〜ＳＥｍが実行される。シーケンス部ＳＥ１は、プリパレーションシーケンスＡ、脂肪抑制シーケンスＢ、およびイメージングシーケンスＣを有している。
【００２４】
プリパレーションシーケンスＡは、３つのＲＦパルス（９０ｘパルス、１８０ｙパルス、９０−ｘパルス）を有している。９０ｘパルスは、縦磁化をフリップさせるためのＲＦパルスであり、９０−ｘパルスは、横磁化を縦磁化に戻すためのＲＦパルスである。９０ｘパルスと９０−ｘパルスとの間に、１８０ｙパルスが送信される。これらのＲＦパルスの間には、拡散強調を行うためのＭＰＧが印加される。ＭＰＧの強さを表すｂ値は、ｂ＝ｂ１に設定されている。９０−ｘパルスの直後に、横磁化を消去するためのキラーパルスが印加される。
【００２５】
プリパレーションシーケンスＡに続いて、脂肪抑制シーケンスＢが実行される。脂肪抑制シーケンスＢを実行した後に、イメージングシーケンスＣが実行される。
【００２６】
イメージングシーケンスＣは、撮影部位からＭＲ信号を収集するためのシーケンスである。イメージングシーケンスＣはＲＦパルスｘ０およびＲＦパルスｘ１〜ｘＮを有している。ＲＦパルスｘ０は、ＭＲ信号を定常状態に移行させるためのＲＦパルスである。尚、ＲＦパルスｘ０の代わりに、ランプアップパルスを用いてＭＲ信号を定常状態に移行させてもよい。
【００２７】
ＲＦパルスｘ０を送信した後に、ＭＲ信号を収集するためのＲＦパルスｘ１〜ｘＮが送信される。尚、説明の便宜上、イメージングシーケンスＣの勾配磁場は図示省略している。
【００２８】
図３には、シーケンス部ＳＥ１の構成が示されているが、他のシーケンス部ＳＥ２〜ＳＥｍも、シーケンス部ＳＥ１と同様に、プリパレーションシーケンスＡ、脂肪抑制シーケンスＢ、およびイメージングシーケンスＢを有している。
拡散強調スキャンＤＷ１を実行した後、次の拡散強調スキャンＤＷ２が実行される（図４参照）。
【００２９】
図４は拡散強調スキャンＤＷ２の説明図である。
拡散強調スキャンＤＷ２も、拡散強調スキャンＤＷ１と同様に、プリパレーションシーケンスＡ、脂肪抑制シーケンスＢ、およびイメージングシーケンスＣが繰り返し実行される。ただし、拡散強調スキャンＤＷ２で使用されるＭＰＧのｂ値は、拡散強調スキャンＤＷ１で使用されるＭＰＧのｂ値とは異なる値（ｂ＝ｂ２）に設定されている。
【００３０】
拡散強調スキャンＤＷ２を実行した後、拡散強調スキャンＤＷ１およびＤＷ２により収集されたデータに基づいて、ＡＤＣを算出することができる。
【００３１】
しかし、拡散強調スキャンＤＷ１およびＤＷ２では、プリパレーションシーケンスＡにおける９０−ｘパルスを送信してから、イメージングシーケンスＣのＲＦパルスｘ１を送信するまでの間に（時間Ｔｒｅｃ）、キラーパルスが印加され、更に脂肪抑制シーケンスＢが実行される。したがって、時間Ｔｒｅｃの間にＴ１回復が起こり、このＴ１回復が原因で、ＡＤＣの計算値とＡＤＣの実際の値との差が大きくなるという問題がある。そこで、本形態では、時間ＴｒｅｃにおけるＴ１回復の影響が補正されたＡＤＣを計算している。以下に、本形態において、Ｔ１回復の影響が補正されたＡＤＣを計算する方法について説明する。
【００３２】
先ず、拡散強調スキャンＤＷ１の時間Ｔｒｅｃ（図３参照）におけるＴ１回復について考える。時間ＴｒｅｃにおけるＴ１回復は、以下の式で表すことができる。
【数１】
2015123305.tif 000003
【００３３】
また、ＭDP1とＡＤＣとの関係は、以下の式で表すことができる。
【数２】
2015123305.tif 000004
【００３４】
式（２）を式（１）に代入し整理すると、以下の式が得られる。
【数３】
2015123305.tif 000005
【００３５】
ここで、拡散強調スキャンＤＷ１において、撮影部位内の位置（ｘ，ｙ，ｚ）から収集されるＭＲ信号の信号値を「Ｓ１」とする。Ｓ１がＭｚ１に比例すると仮定すると、Ｓ１は、以下の式で表すことができる。
【数４】
2015123305.tif 000006
式（４）において、Ｓ０は、Ｍ０の磁化を９０°倒したときに収集されるＭＲ信号の信号値を表している。
【００３６】
次に、拡散強調スキャンＤＷ２において収集されるＭＲ信号の信号値Ｓ２について考える。拡散強調スキャンＤＷ２ではｂ＝ｂ２（図４参照）である。したがって、拡散強調スキャンＤＷにおいて収集されるＭＲ信号の信号値Ｓ２は、式（４）のＳ１およびｂ１をそれぞれＳ２およびｂ２に置き換えることによって表すことができる。つまり、以下の式（５）で表すことができる。
【数５】
2015123305.tif 000007
【００３７】
式（４）および式（５）からＳ０を消去することにより、以下の式が得られる。
【数６】
2015123305.tif 000008
【００３８】
式（６）の中で、Ｓ１は、拡散強調スキャンＤＷ１により収集されたＭＲ信号の信号値であり、Ｓ２は、拡散強調スキャンＤＷ２により収集されたＭＲ信号の信号値である。したがって、Ｓ１およびＳ２は、それぞれ拡散強調スキャンＤＷ１およびＤＷ２を実行することにより知ることができる値である。また、ｂ１は拡散強調スキャンＤＷ１のｂ値であり、ｂ２は拡散強調スキャンＤＷ２のｂ値であるので、ｂ１およびｂ２は既知の値である。更に、ＴＤＰおよびＴｒｅｃは、拡散強調スキャンＤＷ１（ＤＷ２）によって決まる値であるので、ＴＤＰおよびＴｒｅｃも既知の値である。したがって、式（６）のＴ１およびＴ２を求めることができれば、Ｄを計算することができる。そこで、本形態では、被検体を撮影する場合、拡散強調スキャンＤＷ１およびＤＷ２の他に、Ｔ１およびＴ２を計測するためのスキャンも実行される（図５参照）。
【００３９】
図５は本形態で実行されるスキャンの説明図である。
本形態では、Ｔ１値計測用スキャンＳＣ１、Ｔ２値計測用スキャンＳＣ２、および拡散強調スキャンＤＷ１、拡散強調スキャンＤＷ２が実行される。
【００４０】
Ｔ１値計測スキャンＳＣ１は、撮影部位のＴ１値を計測するためのスキャンである。Ｔ１値を計測する手法としては、例えば、Look-Locker法を用いることができる。
【００４１】
Ｔ２計測スキャンＳＣ２は、撮影部位のＴ２値を計測するためのスキャンである。Ｔ２値計測スキャンでは、例えば、マルチエコースピンエコーシーケンスを用いることができる。
【００４２】
拡散強調スキャンＤＷ１およびＤＷ２は、図２〜図４に示す拡散強調スキャンＤＷ１およびＤＷ２と同じである。拡散強調スキャンＤＷ１はｂ＝ｂ１であり、一方、拡散強調スキャンＤＷ２はｂ＝ｂ２（≠ｂ１）である。
【００４３】
以下に、図５に示すスキャンを実行し、ＡＤＣデータを求めるフローについて説明する。
【００４４】
図６はＡＤＣデータを求めるときのフローを示す図である。
ステップＳＴ１では、Ｔ１値計測用スキャンＳＣ１を実行する。データ処理手段８１（図１参照）は、Ｔ１値計測用スキャンＳＣ１を実行することにより収集されたデータに基づいて、撮影部位のＴ１値を求める。図７に、撮影部位のＴ１値を表すＴ１値データＥ１を概略的に示す。Ｔ１値計測用スキャンＳＣ１を実行した後、ステップＳＴ２に進む。
【００４５】
ステップＳＴ２では、Ｔ２値計測用スキャンＳＣ２を実行する。データ処理手段８１は、Ｔ２値計測用スキャンＳＣ２を実行することにより収集されたデータに基づいて、撮影部位のＴ２値を求める。図８に、撮影部位のＴ２値を表すＴ２値データＥ２を概略的に示す。Ｔ２値計測用スキャンＳＣ２を実行した後、ステップＳＴ３に進む。
【００４６】
ステップＳＴ３では、拡散強調スキャンＤＷ１が実行される。データ処理手段８１は、拡散強調スキャンＤＷ１を実行することにより収集されたデータに基づいて、撮影部位の各位置（ｘ，ｙ，ｚ）におけるＭＲ信号の信号値を求める。図９に、拡散強調スキャンＤＷ１により得られた撮影部位の信号値を表す信号値データＥ３を概略的に示す。拡散強調スキャンＤＷ１はｂ＝ｂ１であるので、拡散強調スキャンＤＷ１を実行することにより、ｂ＝ｂ１における撮影部位の信号値データＥ３を取得することができる。拡散強調スキャンＤＷ１を実行した後、ステップＳＴ４に進む。
【００４７】
ステップＳＴ４では、拡散強調スキャンＤＷ２が実行される。データ処理手段８１は、拡散強調スキャンＤＷ２を実行することにより収集されたデータに基づいて、撮影部位の各位置（ｘ，ｙ，ｚ）におけるＭＲ信号の信号値を求める。図１０に、拡散強調スキャンＤＷ２により得られた撮影部位の信号値を表す信号値データＥ４を概略的に示す。拡散強調スキャンＤＷ２はｂ＝ｂ２であるので、拡散強調スキャンＤＷ２を実行することにより、ｂ＝ｂ２における撮影部位の信号値データＥ４を取得することができる。拡散強調スキャンＤＷ２を実行した後、ステップＳＴ５に進む。
【００４８】
ステップＳＴ５では、ＡＤＣ算出手段８２（図１参照）が、式（６）に基づいて撮影部位のＡＤＣを算出する。図１１にＡＤＣの算出方法を概略的に示す。ＡＤＣ算出手段８２は、先ず、データＥ１〜Ｅ４の各々から、同じボクセルのデータを取り出す。図１１では、データＥ１〜Ｅ４の各々からボクセルＶｉのデータが取り出された例が示されている。取り出されたデータは、式（６）に代入される。
【００４９】
式（６）のＴＤＰおよびＴｒｅｃは、拡散強調スキャンＤＷ１（ＤＷ２）によって決まる値である。また、ｂ１は拡散強調スキャンＤＷ１のｂ値である。したがって、ボクセルＶｉにおけるＡＤＣ値（Ｄ＝Ｄｉ）を算出することができる。
【００５０】
図１１では、ボクセルＶｉにおけるＡＤＣ値を求める場合について説明したが、他のボクセルのＡＤＣ値も、式（６）を用いて計算することができる。したがって、撮影部位内の各位置におけるＡＤＣ値を表すＡＤＣデータＥ５を取得することができる。
【００５１】
本形態では、Ｔ１値およびＴ２値が考慮された式（６）を用いてＡＤＣ値を計算している。したがって、時間ＴｒｅｃにおけるＴ１回復の影響が補正されたＡＤＣ値を求めることができる。
【００５２】
尚、上記の説明では、時間Ｔｒｅｃは、プリパレーションシーケンスＡにおける９０−ｘパルスを送信してから、イメージングシーケンスＣのＲＦパルスｘ１が送信されるまでの時間として定義されている。しかし、時間Ｔｒｅｃは、この定義に限定されることはない。図１２に、時間Ｔｒｅｃの他の例を示す。ＲＦパルスｘ０がＴ１回復にほとんど影響を与えないと考えられる場合は、Ｔｒｅｃ＝Ｔｒｅｃ１とすることができる。また、キラーパルスがＴ１回復にほとんど影響を与えないと考えられる場合は、Ｔｒｅｃ＝Ｔｒｅｃ２とすることができる。更に、ＲＦパルスｘ０およびキラーパルスがＴ１回復にほとんど影響を与えないと考えられる場合は、Ｔｒｅｃ＝Ｔｒｅｃ３とすることができる。
【００５３】
また、式（６）には、ＴＤＰが含まれている。しかし、ＴＤＰがＴ２よりも十分に短い場合には、式（６）は以下の式に簡略化することができる。
【数７】
2015123305.tif 000009
【００５４】
式（７）は、式（６）と比較すると、ＴＤＰとＴ２が消去されている。したがって、式（７）を用いてＡＤＣ値を計算する場合は、Ｔ２計測用スキャンＳＣ２が不要となるので、スキャン時間を短縮することができる。
【００５５】
また、本形態では、ＭＲ信号の信号値Ｓ１はＭｚ１に比例すると仮定して式（６）を導出した。しかし、実際には、イメージングシーケンスＣが実行されている間にもＴ１回復するので、イメージングシーケンスＣを実行している途中で、上記の仮定が成り立たなくなることがある。イメージングシーケンスＣを実行している途中で上記の仮定が成り立たなくなる場合、Ｔｒｅｃの間のＴ１回復を考慮するだけでは、信頼性の高いＡＤＣ値を求めることができない。したがって、より信頼性の高いＡＤＣ値を求めるためには、イメージングシーケンスＣが実行されている間にＴ１回復が進んでも、上記の仮定が成り立つようにイメージングシーケンスＣを構成することが望ましい。図１３に、上記の仮定が成り立つように構成されたイメージングシーケンスＣの一例を示す。図１３のイメージングシーケンスＣは、ＲＦパルスＸ１〜ＸＮのフリップ角α１〜αＮが次第に小さくなるように設定されている。このように、フリップ角を次第に小さくすることにより、イメージングシーケンスＣを実行している間にＴ１回復が進んでも、イメージングシーケンスＣを開始してから終了するまでの間に渡って、ＭＲ信号の信号値Ｓ１がＭｚ１に比例にする関係を保持することができるので、より信頼性の高いＡＤＣ値を求めることができる。
【００５６】
また、本形態では、式（６）を解くことによりＡＤＣ値を求めたが、ｂ１とｂ２との値によっては、式（６）を簡単に解くことができる。例えば、ｂ２＝２ｂ１となるように、ｂ１およびｂ２を設定した場合、式（６）は以下の式で表すことができる。
【数８】
2015123305.tif 000010
【００５７】
ここで、以下の式で表されるＸを定義する。
【数９】
2015123305.tif 000011
【００５８】
式（９）を用いると、式（８）はＸについての二次方程式で表すことができる。したがって、Ｘは、解の公式を用いて以下の式で表すことができる。
【数１０】
2015123305.tif 000012
【００５９】
式（９）を式（１０）に代入し整理すると、以下の式が得られる。
【数１１】
2015123305.tif 000013
【００６０】
式（１１）は、ルート（root）の前の符号が±である。したがって、Ｄの値は２つ得られることがわかる。しかし、実際には、Ｄの値は一つしかないので、ルートの前の符号は、−（マイナス）か＋（プラス）のどちらかである。そこで、ルートの前の符号が−（マイナス）であるべきか＋（プラス）であるべきかを決定するために、Ｔｒｅｃが十分に小さく、Ｔｒｅｃ＝０と見なせる場合について考える。Ｔｒｅｃ＝０の場合、ルートの前の符号を−（マイナス）にすると、Ｄ＝∞（無限大）と計算される。しかし、実際にはＤは有限の値であるので、ルートの前の符号を−（マイナス）にすると、Ｄの計算誤差が大きいと考えられる。したがって、ルートの前の符号は＋（プラス）
とすべきであると考えられるので、Ｄは以下の式で計算することができる。
【数１２】
2015123305.tif 000014
【００６１】
したがって、ｂ２＝２ｂ１に設定した場合は、式（１２）を用いてＡＤＣ値を算出することができる。
【００６２】
また、ｂ１又はｂ２のいずれか一方のｂ値をゼロとしてもよい。例えば、ｂ２＝０とすると、ｂ２＝０を式（６）に代入することにより以下の式が得られる。
【数１３】
2015123305.tif 000015
【００６３】
式（１３）を変形すると、以下の式が得られる。
【数１４】
2015123305.tif 000016
【００６４】
式（１４）を用いてＡＤＣ値を計算する場合、ｂ２＝０であるので、拡散強調スキャンＤＷ２はＭＰＧが不要となる。したがって、拡散強調スキャンＤＷ２の代わりに、ＭＰＧを備えていないスキャンを実行することにより、式（１４）を用いてＡＤＣを計算することができる。図１４に、式（１４）を用いてＡＤＣ値を計算する場合に実行されるスキャンの一例を示す。図１４のスキャンは、図５のスキャンと比較すると、Ｔ１値計測スキャンＳＣ１、Ｔ２値計測スキャンＳＣ２、および拡散強調スキャンＤＷ１が実行される点は同じである。しかし、図１４のスキャンでは、拡散強調スキャンＤＷ２は実行されず、その代わりに、ＭＰＧを備えていないスキャンＳＣ３が実行される。スキャンＳＣ３はＭＰＧを備えていないが、その他の点は、拡散強調スキャンＤＷ２と同じである。
【００６５】
図１４のスキャンを実行する場合は、式（１４）を用いることにより、ＡＤＣ値を計算することができる。
【００６６】
尚、式（１４）には、ＴＤＰが含まれている。しかし、ＴＤＰがＴ２よりも十分に短い場合には、式（１４）は以下の式に簡略化することができる。
【数１５】
2015123305.tif 000017
【００６７】
式（１５）は、式（１４）と比較すると、ＴＤＰとＴ２が含まれていないことがわかる。したがって、式（１５）を用いてＡＤＣ値を計算する場合は、Ｔ２計測用スキャンＳＣ２が不要となるので、スキャン時間を短縮することができる。
【００６８】
また、本形態では、ｂ値が異なる２つの式（４）および（５）を用いてＳ０を消去し、ＡＤＣ値を求めるための式（６）を導出している。しかし、以下に示すｂ値が異なるｎ（≧３）個の式（１６−１）〜（１６−ｎ）を考え、最小二乗法などを用いてＡＤＣを求めるための式を導出してもよい。
【数１６】
2015123305.tif 000018
【００６９】
この場合、図１５に示すように、ｂ値が異なるｎ個の拡散強調スキャンＤＷ１〜ＤＷｎ（ｂ＝ｂ１、ｂ２、ｂ３、・・・ｂｎ）を実行すればよい。尚、拡散強調スキャンＤＷ１〜ＤＷｎのうちの一つを、ＭＰＧを有していないスキャン（ｂ＝０）に置き換えてもよい。
【符号の説明】
【００７０】
２  マグネット
３  テーブル
３ａ  クレードル
４  受信コイル
５  送信器
６  勾配磁場電源
７  受信器
８  制御部
９  操作部
１０  表示部
１１  被検体
２１  ボア
８１  データ処理手段
８２  ＡＤＣ算出手段
１００  ＭＲ装置
【図１】
2015123305.tif 000019
【図２】
2015123305.tif 000020
【図３】
2015123305.tif 000021
【図４】
2015123305.tif 000022
【図５】
2015123305.tif 000023
【図６】
2015123305.tif 000024
【図７】
2015123305.tif 000025
【図８】
2015123305.tif 000026
【図９】
2015123305.tif 000027
【図１０】
2015123305.tif 000028
【図１１】
2015123305.tif 000029
【図１２】
2015123305.tif 000030
【図１３】
2015123305.tif 000031
【図１４】
2015123305.tif 000032
【図１５】
2015123305.tif 000033