(57)【特許請求の範囲】
【請求項１】
励起パルスと、前記励起パルスの位相に対して４５°の位相差を有する複数のリフォーカスパルスとを含むシーケンスを実行する磁気共鳴装置であって、
前記シーケンスにより得られる第１エコー信号に含まれるスピンエコーと前記第１エコー信号の次に得られる第２エコー信号に含まれる刺激エコーとに基づいてＢ１マップを作成するＢ１マップ作成手段を有する、磁気共鳴装置。
【請求項２】
前記Ｂ１マップに基づいてＴ２マップを作成するＴ２マップ作成手段を有する、請求項１に記載の磁気共鳴装置。
【請求項３】
前記シーケンスを繰り返し実行することにより取得されたエコー信号に基づいて、複数の画像データを作成する画像データ作成手段を有し、
前記Ｔ２マップ作成手段は、
前記複数の画像データの各ピクセルの信号強度の実測値と、前記Ｂ１マップのＢ１値とに基づいて、前記Ｔ２マップを作成する、請求項２に記載の磁気共鳴装置。
【請求項４】
前記画像データ作成手段は、
同じ番号のエコー信号を用いて各画像データを作成する、請求項３に記載の磁気共鳴装置。
【請求項５】
前記リフォーカスパルスのフリップ角は１７０°以下である、請求項２〜４のうちのいずれか一項に記載の磁気共鳴装置。
【請求項６】
前記Ｂ１マップ作成手段は、
前記スピンエコーの大きさと、前記刺激エコーの大きさとに基づいて、前記Ｂ１マップを作成する、請求項１〜５のうちのいずれか一項に記載の磁気共鳴装置。
【請求項７】
励起パルスと、前記励起パルスの位相に対して４５°の位相差を有する複数のリフォーカスパルスとを含むシーケンスを実行する磁気共鳴装置であって、
前記シーケンスにより得られる第１エコー信号に含まれるスピンエコーと前記第１エコー信号の次に得られる第２エコー信号に含まれる刺激エコーとに基づいてＴ２マップを作成するＴ２マップ作成手段を有する、磁気共鳴装置。
【請求項８】
前記シーケンスは２個のリフォーカスパルスを含んでいる、請求項１〜７のうちのいずれか一項に記載の磁気共鳴装置。
【請求項９】
励起パルスと、前記励起パルスの位相に対して４５°の位相差を有する複数のリフォーカスパルスとを含むシーケンスを実行する磁気共鳴装置のプログラムであって、
前記シーケンスにより得られる第１エコー信号に含まれるスピンエコーと前記第１エコー信号の次に得られる第２エコー信号に含まれる刺激エコーとに基づいてＢ１マップを作成するＢ１マップ作成処理を計算機に実行させるためのプログラム。
【請求項１０】
励起パルスと、前記励起パルスの位相に対して４５°の位相差を有する複数のリフォーカスパルスとを含むシーケンスを実行する磁気共鳴装置のプログラムであって、
前記シーケンスにより得られる第１エコー信号に含まれるスピンエコーと前記第１エコー信号の次に得られる第２エコー信号に含まれる刺激エコーとに基づいてＴ２マップを作成するＴ２マップ作成処理を計算機に実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体のスキャンを行う磁気共鳴装置、およびその磁気共鳴装置に適用されるプログラムに関する。
【背景技術】
【０００２】
近年、人体の組織のＴ２値を画像化したＴ２マップが、医学的な診断に使用されている。Ｔ２マップを作成する方法として、ＣＰＭＧ（Carr-Purcell-Meiboom-Gill）法を用いた技術が知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１２−０１０７２８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
特許文献１では、スピンエコーを収集し、収集したスピンエコーに基づいて、Ｔ２マップを作成している。しかし、Ｂ１不均一により、リフォーカスパルスのフリップ角が１８０°から外れると、エコー信号に含まれる刺激エコー（Stimulated Echo）の影響が強くなり、Ｔ２の推定値の誤差が大きくなるという問題がある。そこで、Ｔ２の推定値の誤差を小さくする方法として、スピンエコーと刺激エコーとを分離することが考えられる。スピンエコーと刺激エコーとを分離することができれば、刺激エコーの影響を低減することができるので、Ｔ２の推定値の誤差を小さくすることが可能となる。しかし、スピンエコーと刺激エコーとを分離することは容易ではない。したがって、スピンエコーと刺激エコーとを容易に分離することが望まれている。
【課題を解決するための手段】
【０００５】
本発明の第１の態様は、励起パルスと、前記励起パルスの位相に対して４５°の位相差を有する複数のリフォーカスパルスとを含むシーケンスを実行する磁気共鳴装置である。
【０００６】
本発明の第２の態様は、励起パルスと、前記励起パルスの位相に対して４５°の位相差を有する複数のリフォーカスパルスとを含むシーケンスを実行する磁気共鳴装置のプログラムであって、前記シーケンスにより得られる第１エコー信号に含まれるスピンエコーと、前記第１エコー信号の次に得られる第２エコー信号に含まれる刺激エコーとに基づいて、Ｂ１マップを作成するＢ１マップ作成処理を計算機に実行させるためのプログラムである。
【発明の効果】
【０００７】
励起パルスと、前記励起パルスの位相に対して４５°の位相差を有する複数のリフォーカスパルスとを含むシーケンスにより、スピンエコーと刺激エコーとを容易に分離することができる。
【図面の簡単な説明】
【０００８】
【図１】本発明の一形態の磁気共鳴装置の概略図である。
【図２】本形態で実行されるシーケンスを示す図である。
【図３】エコー信号の位相プロファイルの説明図である。
【図４】Ｂ１マップおよびＴ２マップを作成するときのフローを示す図である。
【図５】シーケンスＳ１〜Ｓｎを実行することにより得られるエコー信号を概略的に示す図である。
【図６】画像データを作成するときの説明図である。
【図７】Ｂ１マップの作成方法の説明図である。
【図８】画像データＤ１〜ＤｍのピクセルＰｉにおける信号強度Ｃ１ｉ〜ＣｍｉとＴ２との関係を求めるときの説明図である。
【図９】Ｂ１＝ｂをＳ１ｉ（Ｂ１，Ｔ２）〜Ｓｍｉ（Ｂ１，Ｔ２）に代入した様子を示す図である。
【図１０】画像データＤ１〜ＤｍのピクセルＰｉの信号強度の実測値を概略的に示す図である。
【図１１】Ｔ２の値を算出するときの説明図である。
【図１２】シミュレーション結果を示す図である。
【発明を実施するための形態】
【０００９】
以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。
【００１０】
図１は、本発明の一形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ。ＭＲ：Magnetic Resonance）１００は、マグネット２、テーブル３、受信コイル４などを有している。
【００１１】
マグネット２は、被検体１１が収容されるボア２１、超伝導コイル２２、勾配コイル２３、送信コイル２４などを有している。超伝導コイル２２は静磁場を印加し、勾配コイル２３は勾配磁場を印加し、送信コイル２４はＲＦパルスを送信する。尚、超伝導コイル２２の代わりに、永久磁石を用いてもよい。
【００１２】
テーブル３は、被検体１１を支持するためのクレードル３ａを有している。クレードル３ａがボア２１に移動することによって、被検体１１がボアに搬入される。
受信コイル４は被検体１１に取り付けられている。
【００１３】
ＭＲ装置１００は、更に、送信器５、勾配磁場電源６、受信器７、制御部８、操作部９、および表示部１０を有している。
【００１４】
送信器５は、ＲＦコイル２４に電流を供給する。
勾配磁場電源６は、勾配コイル２３に電流を供給する。
受信器７は、受信コイル４で受信された磁気共鳴信号を信号処理し、制御部８に出力する。
【００１５】
制御部８は、表示部１０に必要な情報を伝送したり、受信器７から受け取ったデータに基づいて画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。制御部８は、例えばコンピュータ（computer）によって構成される。制御部８は、画像データ作成手段８１、Ｂ１マップ作成手段８２、Ｔ２マップ作成手段８３などを有している。
【００１６】
画像データ作成手段８１は、シーケンスＳ１〜Ｓｎ（図２参照）より得られたエコー信号に基づいて、画像データを作成する。Ｂ１マップ作成手段８２はＢ１マップ（送信磁場マップ）を作成し、Ｔ２マップ作成手段８３はＴ２マップを作成する。
【００１７】
制御部８は、画像データ作成手段８１、Ｂ１マップ作成手段８２、Ｔ２マップ作成手段８３を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。
【００１８】
操作部９は、オペレータにより操作され、種々の情報を制御部８に入力する。表示部１０は種々の情報を表示する。
ＭＲ装置１００は、上記のように構成されている。
【００１９】
本形態では、ＭＲ装置１００を用いて、Ｂ１マップおよびＴ２マップを取得する。以下に、Ｂ１マップおよびＴ２マップを取得する方法について説明する。
【００２０】
図２は、本形態で実行されるシーケンスを示す図である。
シーケンスＳ１〜Ｓｎは、ＦＳＥ法を用いたシーケンスである。シーケンスＳ１は、励起パルスＲＦ０と、複数のリフォーカスパルスＲＦ１〜ＲＦｍとを有している。本形態では、励起パルスＲＦ０は９０°パルスであり、リフォーカスパルスＲＦ１〜ＲＦｍは１８０°パルスであるが、別の角度のＲＦパルスを用いてもよい。励起パルスＲＦ０とリフォーカスパルスＲＦ１〜ＲＦｍは、位相差が４５°になるように設定されている。図２では、励起パルスＲＦ０の位相φをφ＝０°とし、リフォーカスパルスＲＦ１〜ＲＦｍの位相φをφ＝４５°とすることにより、位相差を４５°に設定しているが、位相φはφ＝０°、４５°に限定されることはない。
【００２１】
シーケンスＳ１では、励起パルスＲＦ０と、リフォーカスパルスＲＦ１〜ＲＦｍを送信することにより、第１エコー信号Ｓｅｉ〜第ｍエコー信号Ｓｅｍが得られる。
【００２２】
他のシーケンスＳ２〜Ｓｎも、シーケンスＳ１と同様に、励起パルスＲＦ０とリフォーカスパルスＲＦ１〜ＲＦｍとを有している。したがって、各シーケンスＳ１〜Ｓｎを実行するたびに、第１エコー信号Ｓｅ１〜第ｍエコー信号Ｓｅｍが得られる。尚、図２では説明の便宜上、勾配磁場は図示省略している。
【００２３】
次に、第１エコー信号Ｓｅ１に含まれるエコー成分と、第２エコー信号Ｓｅ２に含まれるエコー成分ついて、図３を参照しながら説明する。
【００２４】
（１）第１エコー信号Ｓｅ１のエコー成分について
時点ｔ０において励起パルスＲＦ０が送信される。図３には、励起パルスＲＦ０が送信された後の位相変化を表す２つの位相プロファイルＡおよびＢが示されている。位相プロファイルＡは、時点ｔ０において位相が９０°であるが、時点ｔ１において位相が１８０°に変化した例を示している。一方、位相プロファイルＢは、時点ｔ０においては位相が９０°であるが、時点ｔ１において位相が１３５°に変化した例を示している。
【００２５】
次に、時点ｔ１において、１番目のリフォーカスパルスＲＦ１が送信される。リフォーカスパルスＲＦ１の位相は励起パルスＲＦ０の位相に対して４５°ずれているので、時点ｔ１におけるエコー信号の位相は、位相４５°を中心として反転する。その結果、位相プロファイルＡの位相は、１８０°から−９０°になり、位相プロファイルＢの位相は、１３５°から−４５°になる。位相が反転した後、位相プロファイルＡの位相は、時点１１において０°になり、位相プロファイルＢの位相も、時点ｔ１１において０°になる。このように、時点ｔ１１において位相が一致するので、スピンエコーＳＥ（１）が得られる。したがって、第１エコー信号Ｓｅ１の主なエコー成分は、スピンエコーＳＥ（１）となる。
【００２６】
（２）第２エコー信号Ｓｅ２のエコー成分について
スピンエコーＳＥ（１）が得られた後、位相プロファイルＡの位相は次第に大きくなり、時点ｔ２において９０°になる。一方、位相プロファイルＢは、位相プロファイルＡよりも位相変化が緩やかであるので、時点ｔ２において位相はまだ４５°である。
【００２７】
次に、時点ｔ２において、２番目のリフォーカスパルスＲＦ２が送信される。リフォーカスパルスＲＦ２の位相は励起パルスＲＦ０の位相に対して４５°ずれているので、時点ｔ２におけるエコー信号の位相は、位相４５°を中心として反転する。その結果、位相プロファイルＡの位相は、９０°から０°になる。一方、位相プロファイルＢの位相は、時点ｔ２において４５°であるので、リフォーカスパルスＲＦ２が送信されても、位相プロファイルＢの位相は４５°のままである。リフォーカスパルスＲＦ２が送信された後、位相プロファイルＡおよびＢの位相は、時点ｔ２１において９０°になる。このように、時点ｔ２１において位相が一致するので、スピンエコーＳＥ（２）が得られる。したがって、第２エコー信号Ｓｅ２には、スピンエコーＳＥ（２）が含まれている。
【００２８】
ただし、第２エコー信号Ｓｅ２は、スピンエコーＳＥ（２）の他に、刺激エコーＳＴＥ（２）も含んでいる。次に、刺激エコーＳＴＥ（２）がどのようにして生成されるかについて説明する。
【００２９】
時点ｔ１において送信されるリフォーカスパルスＲＦ１は１８０°パルスであるので、フリップ角は理想的には１８０°である。しかし、全てのスピンのフリップ角が１８０°になるわけではなく、フリップ角が１８０°からずれるスピンも存在する。このようなスピンによって、リフォーカスパルスＲＦ１が送信されても位相プロファイルに実質的な変化が見られない場合がある。この例は、位相プロファイルＣおよびＤで示されている。位相プロファイルＣでは、時点ｔ１において位相は１８０°であり、時点ｔ２でも位相は１８０°である。また、位相プロファイルＤでは、時点ｔ１において位相は１３５°であり、時点ｔ２でも位相は１３５°である。
【００３０】
そして、時点ｔ２において、２番目のリフォーカスパルスＲＦ２が送信される。リフォーカスパルスＲＦ２の位相は励起パルスＲＦ０の位相に対して４５°ずれているので、位相プロファイルＣの位相は、１８０°から−９０°になり、位相プロファイルＤの位相は、１３５°から−４５°になる。位相が反転した後、位相プロファイルＣの位相は、時点２１において０°になり、位相プロファイルＤの位相も、時点ｔ２１において０°になる。このように、時点ｔ２１において位相が一致するので、刺激エコーＳＴＥ（２）が得られる。したがって、第２エコー信号Ｓｅ２には、スピンエコーＳＥ（２）だけでなく、刺激エコーＳＴＥ（２）が含まれる。
【００３１】
本形態では、励起パルスＲＦ０とリフォーカスパルスＲＦ１およびＲＦ２は、位相差が４５°になるように設定されている。このように位相差が４５°に設定されたパルスを用いることにより、第２エコー信号Ｓｅ２に含まれるスピンエコーＳＥ（２）と刺激エコーＳＴＥ（２）との位相を９０°ずらすことができる。したがって、スピンエコーＳＥ（２）と刺激エコーＳＴＥ（２）との分離を容易に行うことができる。
【００３２】
また、本形態では、スピンエコーＳＥ（２）の位相と刺激エコーＳＴＥ（２）の位相を９０°ずらすことができるシーケンスＳ１〜Ｓｎを用いて、Ｂ１マップおよびＴ２マップを作成する。以下に、Ｂ１マップおよびＴ２マップを作成する手順について説明する。
【００３３】
図４は、Ｂ１マップおよびＴ２マップを作成するときのフローを示す図である。
ステップＳＴ１では、シーケンスＳ１〜Ｓｎを実行する。図５は、シーケンスＳ１〜Ｓｎを実行することにより得られるエコー信号を概略的に示す図である。１回目のシーケンスＳ１を実行することにより、第１エコー信号Ｓｅ１〜第ｍエコー信号Ｓｅｍが得られる。１回目のシーケンスＳ１が終了した後、２回目〜ｎ回目のシーケンスＳ２〜Ｓｎが順に実行される。２回目〜ｎ回目のシーケンスＳ２〜Ｓｎについても、１回目のシーケンスＳ１と同様に、第１エコー信号Ｓｅ１〜第ｍエコー信号Ｓｅｍが得られる。シーケンスＳ１〜Ｓｎを実行した後、ステップＳＴ２に進む。
【００３４】
ステップＳＴ２では、画像データ作成手段８１（図１参照）が、ステップＳＴ１で取得されたエコー信号に基づいて、画像データを作成する。
【００３５】
図６は、画像データを作成するときの説明図である。
画像データＤ１〜Ｄｍは、１回目〜ｎ回目のシーケンスＳ１〜Ｓｎで得られた同じ番号のエコー信号を用いて作成される。本形態では、各シーケンスＳ１〜Ｓｎでｍ個のエコー信号が得られているので、ｍ個の画像データＤ１〜Ｄｍが作成される。
画像データＤ１〜Ｄｍを作成した後、ステップＳＴ３に進む。
【００３６】
ステップＳＴ３では、Ｂ１マップ作成手段８２（図１参照）が、画像データＤ１およびＤ２に基づいて、Ｂ１マップを作成する。以下に、Ｂ１マップの作成方法について説明する。
【００３７】
図７は、Ｂ１マップの作成方法の説明図である。
画像データＤ１の各ピクセルＰｉ（ｉ＝１〜ｚ）において、第１エコー信号Ｓｅ１に含まれるスピンエコーＳＥ（１）の大きさをＭｘ１とする。また、画像データＤ２のピクセルＰｉにおいて、第２エコー信号Ｓｅ２に含まれる刺激エコーＳＴＥ（２）の大きさをＭｘ２とする。この場合、Ｍｘ１およびＭｘ２は、以下の式で表される。
【数１】
0006033572.tif 000002
【００３８】
尚、式（１）および（２）の右辺には、「ａ」および「ｂ」が含まれている。ａは励起パルスＲＦ０による実際のフリップ角であり、ｂはリフォーカスパルスによる実際のフリップ角である。
【００３９】
一般的に、Ｔ１値は、エコー間隔ＥＳＰよりも十分に大きいので、式（３）のＥ１は、Ｅ１＝１に近似できる。Ｅ１＝１とすると、式（１）および（２）から、「Ｅ２」および「ａ」を消去することができ、以下の式（５）が得られる。
【数２】
0006033572.tif 000003
【００４０】
式（５）を参照すると、ｂの値は、第１エコー信号Ｓｅ１に含まれるスピンエコーＳＥ（１）の大きさＭｘ１と、第２エコー信号Ｓｅ２に含まれる刺激エコーＳＴＥ（２）の大きさＭｘ２とによって表されることがわかる。したがって、ピクセルＰｉごとにＭｘ１の値とＭｘ２の値を求め、求めた値を式（５）に代入すれば、ピクセルＰｉごとにｂの値を算出することができる。
【００４１】
尚、第２エコー信号Ｓｅ２は、刺激エコーＳＴＥ（２）だけでなく、スピンエコーＳＥ（２）も含んでいる（図３参照）。したがって、式（５）に刺激エコーＳＴＥ（２）の大きさＭｘ２を代入するためには、第２エコー信号Ｓｅ２の中から、スピンエコーＳＥ（２）と刺激エコーＳＴＥ（２）とを分離し、刺激エコーＳＴＥ（２）を抽出する必要がある。本形態では、図３を参照しながら説明したように、スピンエコーＳＥ（２）の位相と刺激エコーＳＴＥ（２）の位相を９０°ずらすことができるので、第２エコー信号Ｓｅ２の中から、スピンエコーＳＥ（２）と刺激エコーＳＴＥ（２）とを容易に分離することができる。したがって、式（５）を用いて、各ピクセルごとにｂの値を算出することができるので、Ｂ１マップを作成することができる。Ｂ１マップを作成した後、ステップＳＴ４に進む。
【００４２】
ステップＳＴ４では、Ｔ２マップ作成手段８３（図１参照）がＴ２マップを作成する。
【００４３】
図８〜図１１は、Ｔ２マップを作成するときの説明図である。
本形態では、Ｔ２マップを作成するにあたり、先ず、画像データＤ１〜ＤｍのピクセルＰｉ（ｉ＝１〜ｚ）における信号強度Ｃ１ｉ〜Ｃｍｉと、Ｔ２との関係を求める（図８参照）。
【００４４】
図８は、画像データＤ１〜ＤｍのピクセルＰｉにおける信号強度Ｃ１ｉ〜ＣｍｉとＴ２との関係を求めるときの説明図である。
【００４５】
信号強度Ｃ１ｉ〜Ｃｍｉは、Ｂ１の値とＴ２の値とを変数とした関数Ｓ１ｉ（Ｂ１，Ｔ２）〜Ｓｍ１（Ｂ１，Ｔ２）で表すことができる。したがって、Ｂ１の値がわかれば、信号強度Ｃ１ｉ〜ＣｍｉとＴ２との関係を求めることができる。本形態では、ステップＳＴ３においてＢ１マップを作成しているので、ピクセルＰｉにおけるＢ１の値（Ｂ１＝ｂ）を知ることができる。したがって、Ｂ１＝ｂをＳ１ｉ（Ｂ１，Ｔ２）〜Ｓｍｉ（Ｂ１，Ｔ２）に代入することにより、信号強度Ｃ１ｉ〜ＣｍｉとＴ２との関係を求めることができる。図９に、Ｂ１＝ｂをＳ１ｉ（Ｂ１，Ｔ２）〜Ｓｍｉ（Ｂ１，Ｔ２）に代入した様子を示す。したがって、信号強度Ｃ１ｉ〜ＣｍｉとＴ２との関係を求めることができる。
【００４６】
次に、画像データＤ１〜ＤｍのピクセルＰｉの信号強度の実測値について考える。本形態では、ステップＳＴ２において、画像データＤ１〜Ｄｍを作成しているので、画像データＤ１〜ＤｍのピクセルＰｉの信号強度の実測値は既に求められている。図１０に、画像データＤ１〜ＤｍのピクセルＰｉの信号強度の実測値を概略的に示す。図１０では、画像データＤ１〜ＤｍのピクセルＰｉの信号強度の実測値を、Ｓ１ｉ〜Ｓｍｉで示してある。
【００４７】
そして、ピクセルＰｉの信号強度Ｓ１ｉ（ｂ，Ｔ２）〜Ｓｍｉ（ｂ，Ｔ２）と、ピクセルＰｉの信号強度の実測値Ｓ１ｉ〜Ｓｍｉとに基づいて、Ｔ２の値を算出する。具体的には、以下のようにしてＴ２の値を算出する（図１１参照）。
【００４８】
図１１は、Ｔ２の値を算出するときの説明図である。
Ｔ２マップ作成手段８３（図１参照）は、ピクセルＰｉの信号強度Ｓ１ｉ（ｂ，Ｔ２）〜Ｓｍｉ（ｂ，Ｔ２）と信号強度の実測値Ｓ１ｉ〜Ｓｍｉとの間のばらつきを最小にするＴ２の値を算出する。このようなＴ２値は、例えば以下の式（６）を用いて求めることができる。
【数３】
0006033572.tif 000004
【００４９】
Ｔ２マップ作成手段８３は、式（６）で表される関数Ｆの値が最小になるときのＴ２を算出することにより、信号強度Ｓ１ｉ（ｂ，Ｔ２）〜Ｓｍｉ（ｂ，Ｔ２）と信号強度の実測値Ｓ１ｉ〜Ｓｍｉとの間のばらつきを最小にするＴ２の値を算出することができる。したがって、推定誤差の小さいＴ２の値が得られる。このようにして各ピクセルごとにＴ２の値を算出することができるので、Ｔ２マップを作成することができる。Ｔ２マップを作成した後、フローを終了する。
【００５０】
本形態では、励起パルスＲＦ０とリフォーカスパルスＲＦ１およびＲＦ２は、位相差が４５°になるように設定されている。したがって、スピンエコーＳＥ（２）の位相と刺激エコーＳＴＥ（２）の位相を９０°ずらすことができるので、第２エコー信号Ｓｅ２の中から、スピンエコーＳＥ（２）と刺激エコーＳＴＥ（２）とを容易に分離することができる。また、式（５）を用いて、各ピクセルごとにｂの値を算出することができるので、Ｂ１マップを作成することができる。
【００５１】
本形態では、Ｂ１マップに基づいてＴ２マップを作成しているので、Ｂ１不均一によるＴ２値の誤差を小さくすることができ、信頼性の高いＴ２マップを作成することができる。また、Ｔ２マップを取得するためのシーケンスＳ１〜Ｓｎを実行することにより得られたデータに基づいて、Ｂ１マップも作成することができるので、Ｂ１マップのデータを取得するための専用のシーケンスを実行する必要がなく、スキャン時間を短縮することができる。
【００５２】
本形態では、各シーケンスはｍ個のリフォーカスパルスＲＦ１〜ＲＦｍを含んでいる。ｍは、ｍ≧２であればよいが、ｍの値を大きくしすぎると、スキャン時間が長くなる。したがって、スキャン時間をできるだけ短くしたい場合は、リフォーカスパルスの数を２個とすればよい。２個のリフォーカスパルスを送信すれば、第１エコー信号と第２エコー信号が得られるので、Ｂ１マップおよびＴ２マップを作成することができる。
【００５３】
本形態では、ステップＳＴ３でＢ１マップを作成し、ステップＳＴ４でＴ２マップを作成しているが、Ｂ１マップのみを作成してもよいし、Ｔ２マップのみを作成してもよい。
【００５４】
また、本形態では、リフォーカスパルスＲＦ１〜ＲＦｍの位相は、励起パルスＲＦ０の位相に対して４５°ずれている（図３参照）。このように、位相を４５°ずらすことによって、Ｔ２値の推定誤差を小さくすることができる。このことを検証するため、２つのシミュレーションＱ１およびＱ２を実行した。シミュレーションＱ１は、リフォーカスパルスＲＦ１〜ＲＦｍの位相を、励起パルスＲＦ０の位相に対して９０°ずらしたときのシミュレーションである。一方、シミュレーションＱ２は、本形態のように、リフォーカスパルスＲＦ１〜ＲＦｍの位相を、励起パルスＲＦ０の位相に対して４５°ずらしたときのシミュレーションである。
【００５５】
図１２はシミュレーション結果を示す図である。
図１２（ａ１）、（ｂ１）、および（ｃ１）は、シミュレーションＱ１、即ち位相を９０°ずらしたときのシミュレーション結果であり、図１２（ａ２）、（ｂ２）、および（ｃ２）は、シミュレーションＱ２、即ち位相を４５°ずらしたときのシミュレーション結果である。
【００５６】
尚、図１２（ａ１）および（ａ２）は、平衡状態における縦磁化に対するノイズの標準偏差とＴ２値の推定誤差との関係を示すグラフ、図１２（ｂ１）および（ｂ２）は、Ｔ２値とＴ２値の推定誤差との関係を示すグラフ、図１２（ｃ１）および（ｃ２）は、リフォーカスパルスのフリップ角とＴ２値の推定誤差との関係を示すグラフである。
【００５７】
臨床的な観点から考えると、Ｔ２の推定誤差は３ｍｓ以下にすることが望ましい。そこで、各グラフのＴ２の推定誤差を比較すると、シミュレーションＱ２（位相を４５°ずらした場合）の方が、シミュレーションＱ１（位相を９０°ずらした場合）よりも、広範囲にわたってＴ２値の推定誤差３ｍｓを下回っていることがわかる。
【００５８】
また、図１２（ａ２）、（ｂ２）、および（ｃ２）の各グラフには、２本の折れ線Ｌ１およびＬ２が示されている。折れ線Ｌ１は、本形態の方法で作成したＢ１マップを用いてＴ２値を求めたときのシミュレーション結果であり、一方、折れ線Ｌ２は、本形態の方法で作成したＢ１マップを用いずに、2−parameter fittingによりＴ２値を求めたときのシミュレーション結果である。２本の折れ線Ｌ１およびＬ２を比較すると、本形態の方法で作成したＢ１マップを用いてＴ２値を求めた方が、2−parameter fittingによりＴ２値を求めるよりも、Ｔ２値の推定誤差を更に小さくできることがわかる。
【００５９】
更に、図１２（ｃ２）を参照すると、フリップ角が１７０°以下であれば、Ｔ２の推定誤差は３ｍｓ以下になることがわかる。したがって、リフォーカスパルスのフリップ角は１７０°以下に設定することが望ましい。
【００６０】
本形態では、Ｂ１マップとＴ２マップを作成する例について説明しているが、本発明はＢ１マップとＴ２マップとを作成する例に限定されることはなく、スピンエコーと刺激エコーとを分離することが必要な場合に適用することができる。例えば、スピンエコーと刺激エコーとを分離することにより、エコーの減衰曲線を滑らかにすることができるので、ＦＳＥ（Fast Spin Echo）で撮影する場合に生じるアーチファクトを軽減するのにも有効である。
【符号の説明】
【００６１】
２  マグネット
３  テーブル
３ａ  クレードル
４  受信コイル
５  送信器
６  勾配磁場電源
７  受信器
８  制御部
９  操作部
１０  表示部
１１  被検体
２１  ボア
２２  超伝導コイル
２３  勾配コイル
２４  送信コイル
８１  画像データ作成手段
８２  Ｂ１マップ作成手段
８３  Ｔ２マップ作成手段
１００  ＭＲ装置
【図１】
0006033572.tif 000005
【図２】
0006033572.tif 000006
【図３】
0006033572.tif 000007
【図４】
0006033572.tif 000008
【図５】
0006033572.tif 000009
【図６】
0006033572.tif 000010
【図７】
0006033572.tif 000011
【図８】
0006033572.tif 000012
【図９】
0006033572.tif 000013
【図１０】
0006033572.tif 000014
【図１１】
0006033572.tif 000015
【図１２】
0006033572.tif 000016