JP2015123305A - 磁気共鳴装置およびプログラム - Google Patents

Abstract

【課題】ＡＤＣ（拡散係数）の算出精度を向上させることのできる磁気共鳴装置を提供する。【解決手段】ＭＰＧ（Ｍｏｔｉｏｎ Ｐｒｏｂｉｎｇ Ｇｒａｄｉｅｎｔ）を有するプリパレーションシーケンスと撮影部位からＭＲ信号を収集するためのイメージングシーケンスとを含む拡散強調スキャンを実行するＭＲ装置である。ＭＲ装置は、ＭＰＧのｂ値を変更し、ＭＰＧのｂ値が異なる拡散強調スキャンＤＷ１およびＤＷ２を実行するとともに、撮影部位のＴ１値を計測するためのＴ１値計測用スキャンＳＣ１、および撮影部位のＴ２値を計測するためのＴ２値計測用スキャンＳＣ２を実行するスキャン手段と、時間ＴＤＰ、Ｔｒｅｃ、ＭＲ信号の信号値Ｓ１、Ｓ２、Ｔ１値、およびＴ２値に基づいて、撮影部位のＡＤＣを算出する算出手段とを有する。【選択図】図６

Description

本発明は、拡散強調スキャンを実行する磁気共鳴装置、および磁気共鳴装置に適用されるプログラムに関する。

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

特開２０１２−１５７６８７号公報

ＭＰＧを用いて拡散の情報を得る方法の一例として、ＭＰＧをプリパルスとして印加したDiffusion Preparation（ＤＰ）と呼ばれる手法がある。この手法は、例えばＦＩＥＳＴＡ（Fast Imaging Employing Steady State Acquisition）でデータ収集する場合に適用されており、見かけの拡散係数ＡＤＣ（Apparent Diffusion Coefficient）などの拡散情報を比較的容易に得ることができる。しかし、ＤＰ法では、Ｔ１回復の影響を受けて、ＡＤＣの算出精度が悪くなることがある。したがって、ＡＤＣの算出精度を向上させることが望まれている。

本発明の第１の観点は、ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンを実行する磁気共鳴装置であって、
前記ＭＰＧの強さを表すｂ値を変更し、前記ＭＰＧのｂ値が異なる複数の拡散強調スキャンを実行するとともに、前記撮影部位のＴ１値を計測するためのＴ１値計測用スキャンを実行するスキャン手段と、
前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記複数の拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出手段と、
を有する、磁気共鳴装置である。

本発明の第２の観点は、ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンと、ＭＰＧを有していない第３のシーケンスと前記撮影部位からＭＲ信号を収集するための第４のシーケンスとを含む他のスキャンとを実行する磁気共鳴装置であって、
前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、前記他のスキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出手段と、
を有する、磁気共鳴装置である。

本発明の第３の観点は、ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンと、前記撮影部位のＴ１値を計測するためのＴ１値計測用スキャンとを実行する磁気共鳴装置に適用されるプログラムであって、
前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記ＭＰＧの強さを表すｂ値が異なる複数の拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出処理を計算機に実行させるためのプログラムである。

本発明の第４の観点は、ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンと、ＭＰＧを有していない第３のシーケンスと撮影部位からＭＲ信号を収集するための第４のシーケンスとを含む他のスキャンとを実行する磁気共鳴装置に適用されるプログラムであって、
前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、前記他のスキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出処理を計算機に実行させるためのプログラムである。

Ｔ１値を考慮してＡＤＣを算出するので、信頼性の高いＡＤＣが得られる。

本発明の一形態の磁気共鳴装置の概略図である。 ＡＤＣデータを取得するために実行されるスキャンの一例を示す図である。 拡散強調スキャンＤＷ１の説明図である。 拡散強調スキャンＤＷ２の説明図である。 本形態で実行されるスキャンの説明図である。 ＡＤＣデータを求めるときのフローを示す図である。 撮影部位のＴ１値を表すＴ１値データＥ１を概略的に示す図である。 撮影部位のＴ２値を表すＴ２値データＥ２を概略的に示す図である。 拡散強調スキャンＤＷ１により得られた撮影部位の信号値を表す信号値データＥ３を概略的に示す図である。 拡散強調スキャンＤＷ２により得られた撮影部位の信号値を表す信号値データＥ４を概略的に示す図である。 ＡＤＣの算出方法を概略的に示す図である。 時間Ｔ_ｒｅｃの他の例を示す図である。 ＲＦパルスＸ_１〜Ｘ_Ｎのフリップ角α_１〜α_Ｎが次第に小さくなるように設定されたイメージングシーケンスＣを示す図である。 式（１４）を用いてＡＤＣ値を計算する場合に実行されるスキャンの一例を示す図である。 ｂ値が異なるｎ個の拡散強調スキャンＤＷ１〜ＤＷｎ（ｂ＝ｂ１、ｂ２、ｂ３、・・・ｂｎ）を実行する例を示す図である。

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

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

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

テーブル３は、被検体１１を支持するクレードル３ａを有している。クレードル３ａは、ボア２１内に移動できるように構成されている。クレードル３ａによって、被検体１１はボア２１に搬送される。
受信コイル４は、被検体１１からの磁気共鳴信号を受信する。

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

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

制御部８は、表示部１０に必要な情報を伝送したり、受信器７から受け取ったデータに基づいて画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。制御部８は、データ処理手段８１およびＡＤＣ算出手段８２などを有している。

データ処理手段８１は、スキャンにより収集されたデータに基づいて、Ｔ１値、Ｔ２値、ＭＲ信号の信号値などを求める。
ＡＤＣ算出手段８２は、後述する式に基づいて、ＡＤＣを算出する。

制御部８は、データ処理手段８１およびＡＤＣ算出手段８２を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。

操作部９は、オペレータにより操作され、種々の情報を制御部８に入力する。表示部１０は種々の情報を表示する。

ＭＲ装置１００は、上記のように構成されている。本形態では、ＭＲ装置１００を用いてＡＤＣデータを取得する。以下に、ＡＤＣデータを取得するために実行されるスキャンについて説明する。

図２はＡＤＣデータを取得するために実行されるスキャンの一例を示す図である。
図２には、拡散強調データを取得するための２つの拡散強調スキャンＤＷ１およびＤＷ２が示されている。拡散強調スキャンＤＷ１およびＤＷ２は、ＭＰＧをプリパルスとして印加したDiffusion Preparation（ＤＰ）と呼ばれる手法で拡散情報を取得するスキャンである。以下、拡散強調スキャンＤＷ１およびＤＷ２について順に説明する。

図３は拡散強調スキャンＤＷ１の説明図である。
拡散強調スキャンＤＷ１では、シーケンス部ＳＥ_１〜ＳＥ_ｍが実行される。シーケンス部ＳＥ_１は、プリパレーションシーケンスＡ、脂肪抑制シーケンスＢ、およびイメージングシーケンスＣを有している。

プリパレーションシーケンスＡは、３つのＲＦパルス（９０_ｘパルス、１８０_ｙパルス、９０_−ｘパルス）を有している。９０_ｘパルスは、縦磁化をフリップさせるためのＲＦパルスであり、９０_−ｘパルスは、横磁化を縦磁化に戻すためのＲＦパルスである。９０_ｘパルスと９０_−ｘパルスとの間に、１８０_ｙパルスが送信される。これらのＲＦパルスの間には、拡散強調を行うためのＭＰＧが印加される。ＭＰＧの強さを表すｂ値は、ｂ＝ｂ_１に設定されている。９０_−ｘパルスの直後に、横磁化を消去するためのキラーパルスが印加される。

プリパレーションシーケンスＡに続いて、脂肪抑制シーケンスＢが実行される。脂肪抑制シーケンスＢを実行した後に、イメージングシーケンスＣが実行される。

イメージングシーケンスＣは、撮影部位からＭＲ信号を収集するためのシーケンスである。イメージングシーケンスＣはＲＦパルスｘ_０およびＲＦパルスｘ_１〜ｘ_Ｎを有している。ＲＦパルスｘ_０は、ＭＲ信号を定常状態に移行させるためのＲＦパルスである。尚、ＲＦパルスｘ_０の代わりに、ランプアップパルスを用いてＭＲ信号を定常状態に移行させてもよい。

ＲＦパルスｘ_０を送信した後に、ＭＲ信号を収集するためのＲＦパルスｘ_１〜ｘ_Ｎが送信される。尚、説明の便宜上、イメージングシーケンスＣの勾配磁場は図示省略している。

図３には、シーケンス部ＳＥ_１の構成が示されているが、他のシーケンス部ＳＥ_２〜ＳＥ_ｍも、シーケンス部ＳＥ_１と同様に、プリパレーションシーケンスＡ、脂肪抑制シーケンスＢ、およびイメージングシーケンスＢを有している。
拡散強調スキャンＤＷ１を実行した後、次の拡散強調スキャンＤＷ２が実行される（図４参照）。

図４は拡散強調スキャンＤＷ２の説明図である。
拡散強調スキャンＤＷ２も、拡散強調スキャンＤＷ１と同様に、プリパレーションシーケンスＡ、脂肪抑制シーケンスＢ、およびイメージングシーケンスＣが繰り返し実行される。ただし、拡散強調スキャンＤＷ２で使用されるＭＰＧのｂ値は、拡散強調スキャンＤＷ１で使用されるＭＰＧのｂ値とは異なる値（ｂ＝ｂ_２）に設定されている。

拡散強調スキャンＤＷ２を実行した後、拡散強調スキャンＤＷ１およびＤＷ２により収集されたデータに基づいて、ＡＤＣを算出することができる。

しかし、拡散強調スキャンＤＷ１およびＤＷ２では、プリパレーションシーケンスＡにおける９０_−ｘパルスを送信してから、イメージングシーケンスＣのＲＦパルスｘ_１を送信するまでの間に（時間Ｔ_ｒｅｃ）、キラーパルスが印加され、更に脂肪抑制シーケンスＢが実行される。したがって、時間Ｔ_ｒｅｃの間にＴ１回復が起こり、このＴ１回復が原因で、ＡＤＣの計算値とＡＤＣの実際の値との差が大きくなるという問題がある。そこで、本形態では、時間Ｔ_ｒｅｃにおけるＴ１回復の影響が補正されたＡＤＣを計算している。以下に、本形態において、Ｔ１回復の影響が補正されたＡＤＣを計算する方法について説明する。

先ず、拡散強調スキャンＤＷ１の時間Ｔ_ｒｅｃ（図３参照）におけるＴ１回復について考える。時間Ｔ_ｒｅｃにおけるＴ１回復は、以下の式で表すことができる。

また、Ｍ_DP1とＡＤＣとの関係は、以下の式で表すことができる。

式（２）を式（１）に代入し整理すると、以下の式が得られる。

ここで、拡散強調スキャンＤＷ１において、撮影部位内の位置（ｘ，ｙ，ｚ）から収集されるＭＲ信号の信号値を「Ｓ_１」とする。Ｓ_１がＭ_ｚ１に比例すると仮定すると、Ｓ_１は、以下の式で表すことができる。
式（４）において、Ｓ_０は、Ｍ０の磁化を９０°倒したときに収集されるＭＲ信号の信号値を表している。

次に、拡散強調スキャンＤＷ２において収集されるＭＲ信号の信号値Ｓ_２について考える。拡散強調スキャンＤＷ２ではｂ＝ｂ_２（図４参照）である。したがって、拡散強調スキャンＤＷにおいて収集されるＭＲ信号の信号値Ｓ_２は、式（４）のＳ_１およびｂ_１をそれぞれＳ_２およびｂ_２に置き換えることによって表すことができる。つまり、以下の式（５）で表すことができる。

式（４）および式（５）からＳ_０を消去することにより、以下の式が得られる。

式（６）の中で、Ｓ_１は、拡散強調スキャンＤＷ１により収集されたＭＲ信号の信号値であり、Ｓ_２は、拡散強調スキャンＤＷ２により収集されたＭＲ信号の信号値である。したがって、Ｓ_１およびＳ_２は、それぞれ拡散強調スキャンＤＷ１およびＤＷ２を実行することにより知ることができる値である。また、ｂ_１は拡散強調スキャンＤＷ１のｂ値であり、ｂ_２は拡散強調スキャンＤＷ２のｂ値であるので、ｂ_１およびｂ_２は既知の値である。更に、Ｔ_ＤＰおよびＴ_ｒｅｃは、拡散強調スキャンＤＷ１（ＤＷ２）によって決まる値であるので、Ｔ_ＤＰおよびＴ_ｒｅｃも既知の値である。したがって、式（６）のＴ１およびＴ２を求めることができれば、Ｄを計算することができる。そこで、本形態では、被検体を撮影する場合、拡散強調スキャンＤＷ１およびＤＷ２の他に、Ｔ１およびＴ２を計測するためのスキャンも実行される（図５参照）。

図５は本形態で実行されるスキャンの説明図である。
本形態では、Ｔ１値計測用スキャンＳＣ１、Ｔ２値計測用スキャンＳＣ２、および拡散強調スキャンＤＷ１、拡散強調スキャンＤＷ２が実行される。

Ｔ１値計測スキャンＳＣ１は、撮影部位のＴ１値を計測するためのスキャンである。Ｔ１値を計測する手法としては、例えば、Look-Locker法を用いることができる。

Ｔ２計測スキャンＳＣ２は、撮影部位のＴ２値を計測するためのスキャンである。Ｔ２値計測スキャンでは、例えば、マルチエコースピンエコーシーケンスを用いることができる。

拡散強調スキャンＤＷ１およびＤＷ２は、図２〜図４に示す拡散強調スキャンＤＷ１およびＤＷ２と同じである。拡散強調スキャンＤＷ１はｂ＝ｂ_１であり、一方、拡散強調スキャンＤＷ２はｂ＝ｂ_２（≠ｂ_１）である。

以下に、図５に示すスキャンを実行し、ＡＤＣデータを求めるフローについて説明する。

図６はＡＤＣデータを求めるときのフローを示す図である。
ステップＳＴ１では、Ｔ１値計測用スキャンＳＣ１を実行する。データ処理手段８１（図１参照）は、Ｔ１値計測用スキャンＳＣ１を実行することにより収集されたデータに基づいて、撮影部位のＴ１値を求める。図７に、撮影部位のＴ１値を表すＴ１値データＥ１を概略的に示す。Ｔ１値計測用スキャンＳＣ１を実行した後、ステップＳＴ２に進む。

ステップＳＴ２では、Ｔ２値計測用スキャンＳＣ２を実行する。データ処理手段８１は、Ｔ２値計測用スキャンＳＣ２を実行することにより収集されたデータに基づいて、撮影部位のＴ２値を求める。図８に、撮影部位のＴ２値を表すＴ２値データＥ２を概略的に示す。Ｔ２値計測用スキャンＳＣ２を実行した後、ステップＳＴ３に進む。

ステップＳＴ３では、拡散強調スキャンＤＷ１が実行される。データ処理手段８１は、拡散強調スキャンＤＷ１を実行することにより収集されたデータに基づいて、撮影部位の各位置（ｘ，ｙ，ｚ）におけるＭＲ信号の信号値を求める。図９に、拡散強調スキャンＤＷ１により得られた撮影部位の信号値を表す信号値データＥ３を概略的に示す。拡散強調スキャンＤＷ１はｂ＝ｂ_１であるので、拡散強調スキャンＤＷ１を実行することにより、ｂ＝ｂ_１における撮影部位の信号値データＥ３を取得することができる。拡散強調スキャンＤＷ１を実行した後、ステップＳＴ４に進む。

ステップＳＴ４では、拡散強調スキャンＤＷ２が実行される。データ処理手段８１は、拡散強調スキャンＤＷ２を実行することにより収集されたデータに基づいて、撮影部位の各位置（ｘ，ｙ，ｚ）におけるＭＲ信号の信号値を求める。図１０に、拡散強調スキャンＤＷ２により得られた撮影部位の信号値を表す信号値データＥ４を概略的に示す。拡散強調スキャンＤＷ２はｂ＝ｂ_２であるので、拡散強調スキャンＤＷ２を実行することにより、ｂ＝ｂ_２における撮影部位の信号値データＥ４を取得することができる。拡散強調スキャンＤＷ２を実行した後、ステップＳＴ５に進む。

ステップＳＴ５では、ＡＤＣ算出手段８２（図１参照）が、式（６）に基づいて撮影部位のＡＤＣを算出する。図１１にＡＤＣの算出方法を概略的に示す。ＡＤＣ算出手段８２は、先ず、データＥ１〜Ｅ４の各々から、同じボクセルのデータを取り出す。図１１では、データＥ１〜Ｅ４の各々からボクセルＶｉのデータが取り出された例が示されている。取り出されたデータは、式（６）に代入される。

式（６）のＴ_ＤＰおよびＴ_ｒｅｃは、拡散強調スキャンＤＷ１（ＤＷ２）によって決まる値である。また、ｂ_１は拡散強調スキャンＤＷ１のｂ値である。したがって、ボクセルＶｉにおけるＡＤＣ値（Ｄ＝Ｄｉ）を算出することができる。

図１１では、ボクセルＶｉにおけるＡＤＣ値を求める場合について説明したが、他のボクセルのＡＤＣ値も、式（６）を用いて計算することができる。したがって、撮影部位内の各位置におけるＡＤＣ値を表すＡＤＣデータＥ５を取得することができる。

本形態では、Ｔ１値およびＴ２値が考慮された式（６）を用いてＡＤＣ値を計算している。したがって、時間Ｔ_ｒｅｃにおけるＴ１回復の影響が補正されたＡＤＣ値を求めることができる。

尚、上記の説明では、時間Ｔ_ｒｅｃは、プリパレーションシーケンスＡにおける９０_−ｘパルスを送信してから、イメージングシーケンスＣのＲＦパルスｘ_１が送信されるまでの時間として定義されている。しかし、時間Ｔ_ｒｅｃは、この定義に限定されることはない。図１２に、時間Ｔ_ｒｅｃの他の例を示す。ＲＦパルスｘ_０がＴ１回復にほとんど影響を与えないと考えられる場合は、Ｔ_ｒｅｃ＝Ｔ_ｒｅｃ１とすることができる。また、キラーパルスがＴ１回復にほとんど影響を与えないと考えられる場合は、Ｔ_ｒｅｃ＝Ｔ_ｒｅｃ２とすることができる。更に、ＲＦパルスｘ_０およびキラーパルスがＴ１回復にほとんど影響を与えないと考えられる場合は、Ｔ_ｒｅｃ＝Ｔ_ｒｅｃ３とすることができる。

また、式（６）には、Ｔ_ＤＰが含まれている。しかし、Ｔ_ＤＰがＴ２よりも十分に短い場合には、式（６）は以下の式に簡略化することができる。

式（７）は、式（６）と比較すると、Ｔ_ＤＰとＴ２が消去されている。したがって、式（７）を用いてＡＤＣ値を計算する場合は、Ｔ２計測用スキャンＳＣ２が不要となるので、スキャン時間を短縮することができる。

また、本形態では、ＭＲ信号の信号値Ｓ_１はＭ_ｚ１に比例すると仮定して式（６）を導出した。しかし、実際には、イメージングシーケンスＣが実行されている間にもＴ１回復するので、イメージングシーケンスＣを実行している途中で、上記の仮定が成り立たなくなることがある。イメージングシーケンスＣを実行している途中で上記の仮定が成り立たなくなる場合、Ｔ_ｒｅｃの間のＴ１回復を考慮するだけでは、信頼性の高いＡＤＣ値を求めることができない。したがって、より信頼性の高いＡＤＣ値を求めるためには、イメージングシーケンスＣが実行されている間にＴ１回復が進んでも、上記の仮定が成り立つようにイメージングシーケンスＣを構成することが望ましい。図１３に、上記の仮定が成り立つように構成されたイメージングシーケンスＣの一例を示す。図１３のイメージングシーケンスＣは、ＲＦパルスＸ_１〜Ｘ_Ｎのフリップ角α_１〜α_Ｎが次第に小さくなるように設定されている。このように、フリップ角を次第に小さくすることにより、イメージングシーケンスＣを実行している間にＴ１回復が進んでも、イメージングシーケンスＣを開始してから終了するまでの間に渡って、ＭＲ信号の信号値Ｓ_１がＭ_ｚ１に比例にする関係を保持することができるので、より信頼性の高いＡＤＣ値を求めることができる。

また、本形態では、式（６）を解くことによりＡＤＣ値を求めたが、ｂ_１とｂ_２との値によっては、式（６）を簡単に解くことができる。例えば、ｂ_２＝２ｂ_１となるように、ｂ_１およびｂ_２を設定した場合、式（６）は以下の式で表すことができる。

ここで、以下の式で表されるＸを定義する。

式（９）を用いると、式（８）はＸについての二次方程式で表すことができる。したがって、Ｘは、解の公式を用いて以下の式で表すことができる。

式（９）を式（１０）に代入し整理すると、以下の式が得られる。

式（１１）は、ルート（root）の前の符号が±である。したがって、Ｄの値は２つ得られることがわかる。しかし、実際には、Ｄの値は一つしかないので、ルートの前の符号は、−（マイナス）か＋（プラス）のどちらかである。そこで、ルートの前の符号が−（マイナス）であるべきか＋（プラス）であるべきかを決定するために、Ｔ_ｒｅｃが十分に小さく、Ｔ_ｒｅｃ＝０と見なせる場合について考える。Ｔ_ｒｅｃ＝０の場合、ルートの前の符号を−（マイナス）にすると、Ｄ＝∞（無限大）と計算される。しかし、実際にはＤは有限の値であるので、ルートの前の符号を−（マイナス）にすると、Ｄの計算誤差が大きいと考えられる。したがって、ルートの前の符号は＋（プラス）
とすべきであると考えられるので、Ｄは以下の式で計算することができる。

したがって、ｂ_２＝２ｂ_１に設定した場合は、式（１２）を用いてＡＤＣ値を算出することができる。

また、ｂ_１又はｂ_２のいずれか一方のｂ値をゼロとしてもよい。例えば、ｂ_２＝０とすると、ｂ_２＝０を式（６）に代入することにより以下の式が得られる。

式（１３）を変形すると、以下の式が得られる。

式（１４）を用いてＡＤＣ値を計算する場合、ｂ_２＝０であるので、拡散強調スキャンＤＷ２はＭＰＧが不要となる。したがって、拡散強調スキャンＤＷ２の代わりに、ＭＰＧを備えていないスキャンを実行することにより、式（１４）を用いてＡＤＣを計算することができる。図１４に、式（１４）を用いてＡＤＣ値を計算する場合に実行されるスキャンの一例を示す。図１４のスキャンは、図５のスキャンと比較すると、Ｔ１値計測スキャンＳＣ１、Ｔ２値計測スキャンＳＣ２、および拡散強調スキャンＤＷ１が実行される点は同じである。しかし、図１４のスキャンでは、拡散強調スキャンＤＷ２は実行されず、その代わりに、ＭＰＧを備えていないスキャンＳＣ３が実行される。スキャンＳＣ３はＭＰＧを備えていないが、その他の点は、拡散強調スキャンＤＷ２と同じである。

図１４のスキャンを実行する場合は、式（１４）を用いることにより、ＡＤＣ値を計算することができる。

尚、式（１４）には、Ｔ_ＤＰが含まれている。しかし、Ｔ_ＤＰがＴ２よりも十分に短い場合には、式（１４）は以下の式に簡略化することができる。

式（１５）は、式（１４）と比較すると、Ｔ_ＤＰとＴ２が含まれていないことがわかる。したがって、式（１５）を用いてＡＤＣ値を計算する場合は、Ｔ２計測用スキャンＳＣ２が不要となるので、スキャン時間を短縮することができる。

また、本形態では、ｂ値が異なる２つの式（４）および（５）を用いてＳ_０を消去し、ＡＤＣ値を求めるための式（６）を導出している。しかし、以下に示すｂ値が異なるｎ（≧３）個の式（１６−１）〜（１６−ｎ）を考え、最小二乗法などを用いてＡＤＣを求めるための式を導出してもよい。

この場合、図１５に示すように、ｂ値が異なるｎ個の拡散強調スキャンＤＷ１〜ＤＷｎ（ｂ＝ｂ１、ｂ２、ｂ３、・・・ｂｎ）を実行すればよい。尚、拡散強調スキャンＤＷ１〜ＤＷｎのうちの一つを、ＭＰＧを有していないスキャン（ｂ＝０）に置き換えてもよい。

２ マグネット
３ テーブル
３ａ クレードル
４ 受信コイル
５ 送信器
６ 勾配磁場電源
７ 受信器
８ 制御部
９ 操作部
１０ 表示部
１１ 被検体
２１ ボア
８１ データ処理手段
８２ ＡＤＣ算出手段
１００ ＭＲ装置

Claims (11)

1. ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンを実行する磁気共鳴装置であって、
   前記ＭＰＧの強さを表すｂ値を変更し、前記ＭＰＧのｂ値が異なる複数の拡散強調スキャンを実行するとともに、前記撮影部位のＴ１値を計測するためのＴ１値計測用スキャンを実行するスキャン手段と、
   前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記複数の拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出手段と、
   を有する、磁気共鳴装置。
2. ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンと、ＭＰＧを有していない第３のシーケンスと前記撮影部位からＭＲ信号を収集するための第４のシーケンスとを含む他のスキャンとを実行する磁気共鳴装置であって、
   前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、前記他のスキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出手段と、
   を有する、磁気共鳴装置。
3. 前記スキャン手段は、
   前記撮影部位のＴ２値を計測するためのＴ２値計測用スキャンを実行し、
   前記算出手段は、
   前記Ｔ２値計測用スキャンにより計測されたＴ２値も考慮して、前記見かけの拡散係数を算出する、請求項１又は２に記載の磁気共鳴装置。
4. 前記第１のシーケンスは、縦磁化をフリップさせるための第１のＲＦパルスと、横磁化を縦磁化に戻すための第２のＲＦパルスとを有し、
   前記算出手段は、
   前記第１のＲＦパルスと前記第２のＲＦパルスとの間の第２の時間も考慮して、前記見かけの拡散係数を算出する、請求項３に記載の磁気共鳴装置。
5. 前記第２のシーケンスは、前記撮影部位からＭＲ信号を収集するための複数のＲＦパルスを有し、
   前記第１の時間は、前記第２のＲＦパルスと、前記複数のＲＦパルスのうちの一番最初に送信されるＲＦパルスとの間の時間である、請求項１〜４のうちのいずれか一項に記載の磁気共鳴装置。
6. 前記複数のＲＦパルスのフリップ角は次第に小さくなるように設定されている、請求項５に記載の磁気共鳴装置。
7. 前記第２のシーケンスは、前記複数のＲＦパルスの前に、ＭＲ信号を定常状態に移行させるためのＲＦパルスを有する、請求項５又は６に記載の磁気共鳴装置。
8. 前記第１のシーケンスは、横磁化を消去するためのキラーパルスを有する、請求項１〜７のうちのいずれか一項に記載の磁気共鳴装置。
9. 前記拡散強調スキャンは、第１のシーケンスと第２のシーケンスとの間に、脂肪を抑制するための第３のシーケンスを有する、請求項１〜８のうちのいずれか一項に記載の磁気共鳴装置。
10. ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンと、前記撮影部位のＴ１値を計測するためのＴ１値計測用スキャンとを実行する磁気共鳴装置に適用されるプログラムであって、
    前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記ＭＰＧの強さを表すｂ値が異なる複数の拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出処理を計算機に実行させるためのプログラム。
11. ＭＰＧを有する第１のシーケンスと撮影部位からＭＲ信号を収集するための第２のシーケンスとを含む拡散強調スキャンと、ＭＰＧを有していない第３のシーケンスと撮影部位からＭＲ信号を収集するための第４のシーケンスとを含む他のスキャンとを実行する磁気共鳴装置に適用されるプログラムであって、
    前記第１のシーケンスと前記第２のシーケンスとの間の第１の時間、前記拡散強調スキャンを実行することにより収集されたＭＲ信号の信号値、前記他のスキャンを実行することにより収集されたＭＲ信号の信号値、および前記Ｔ１値計測用スキャンにより計測されたＴ１値に基づいて、前記撮影部位の見かけの拡散係数を算出する算出処理を計算機に実行させるためのプログラム。