(57)【要約】  短時間で高精度にRF送信コイルのB1分布を計測することが可能なMRI装置を提供する。そのために、MRI装置の撮像手段は、RF照射手段によるプリパルスの印加と、前記プリパルスの印加からの経過時間(TI)が異なる複数の信号取得シーケンスとを含みB1分布計測シーケンスを備える。信号取得シーケンスはRFパルスとしてフリプ角の小さいパルスを用いるとともに、プリパルス後の縦緩和が終わる前に実行される。演算手段は、複数の信号取得シーケンスにより、それぞれ得られたTIの異なる画像データを用いて、RF照射手段のB1分布を算出する。
2012060192.tif 000002
【特許請求の範囲】
【請求項１】
被検体に核磁気共鳴を起こさせるための高周波磁場(B1)を照射するRF照射部と、
高周波磁場プリパルスの印加から経過時間(TI)を設けて核磁気共鳴信号を取得する信号取得シーケンスを有して成るB1分布計測シーケンスを用いて前記被検体を撮像する撮像部と、
前記核磁気共鳴信号を用いて前記被検体の画像を再構成する演算部と、
を備え、
前記演算部は、前記画像を用いて前記RF照射部の照射磁場分布を求める磁気共鳴撮像装置であって、
前記演算部は、前記経過時間の異なる複数の画像を用いて前記照射磁場分布を求めることを特徴とする磁気共鳴撮像装置。
【請求項２】
請求項1記載の磁気共鳴撮像装置において、
前記B1分布計測シーケンスは、前記経過時間が異なる複数の信号取得シーケンスを含み、
前記演算部は、前記経過時間が異なる複数の信号取得シーケンスでそれぞれ取得された核磁気共鳴信号を用いて、前記経過時間の異なる複数の画像を再構成することを特徴とする磁気共鳴撮像装置。
【請求項３】
請求項1に記載の磁気共鳴撮像装置において、
前記撮像部は、前記高周波磁場プリパルスの印加条件が異なる複数のB1分布計測シーケンスを実行してそれぞれ画像を取得し、
前記演算部は、前記高周波磁場プリパルスの印加条件が異なる複数のB1分布計測シーケンスでそれぞれ得られた画像を用いて前記照射磁場分布を求めることを特徴とする磁気共鳴撮像装置。
【請求項４】
請求項3に記載の磁気共鳴撮像装置において、
前記複数のB1分布計測シーケンスの内で、第一のB1分布計測シーケンスは、前記高周波磁場プリパルスを用いずに実行される信号取得シーケンスであり、第二のB1分布計測シーケンスは、前記高周波磁場プリパルスの印加とそれに続く少なくとも2回の信号取得シーケンスを含むことを特徴とする磁気共鳴撮像装置。
【請求項５】
請求項4に記載の磁気共鳴撮像装置において、
前記演算部は、前記第一のB1分布計測シーケンスの信号取得シーケンスで得られた核磁気共鳴信号を再構成した画像の信号強度と、前記第二のB1分布計測シーケンスの複数の信号取得シーケンスでそれぞれ得られた核磁気共鳴信号を再構成した複数の画像の信号強度との比を用いて、行列式を解くことにより、前記照射磁場分布を算出することを特徴とする磁気共鳴撮像装置。
【請求項６】
請求項1に記載の磁気共鳴撮像装置において、
前記撮像部は、複数の断面について、それぞれ前記B1分布計測シーケンスを実行し、
前記演算部は、前記複数の断面のそれぞれについて、前記照射磁場分布を求めることを特徴とする磁気共鳴撮像装置。
【請求項７】
請求項6に記載の磁気共鳴撮像装置において、
前記信号取得シーケンスは、マルチスライス撮像シーケンス又は三次元撮像シーケンスであることを特徴とする磁気共鳴撮像装置。
【請求項８】
請求項1に記載の磁気共鳴撮像装置において、
前記撮像部は、前記高周波磁場プリパルスの印加条件が同一の複数の前記B1分布計測シーケンスを含み、該複数のB1分布計測シーケンスで前記画像の再構成に必要な全データを分割して取得することを特徴とする磁気共鳴撮像装置。
【請求項９】
請求項1に記載の磁気共鳴撮像装置において、
前記撮像部は、前記経過時間を異ならせて、前記B1分布計測シーケンスを複数回実行し、
前記演算部は、各B1分布計測シーケンスで得られた核磁気共鳴信号を用いて、前記経過時間の異なる複数の画像を再構成することを特徴とする磁気共鳴撮像装置。
【請求項１０】
請求項1に記載の磁気共鳴撮像装置において、
前記高周波磁場プリパルスは、フリップ角が90度以上の高周波磁場パルスであることを特徴とする磁気共鳴撮像装置。
【請求項１１】
請求項2に記載の磁気共鳴撮像装置において、
前記信号取得シーケンスは、フリップ角が10度以下の高周波磁場パルスを繰り返し印加して前記核磁気共鳴信号を取得することを特徴とする磁気共鳴撮像装置。
【請求項１２】
請求項11に記載の磁気共鳴撮像装置において、
前記撮像部は、前記複数の信号取得シーケンスは連続しており、前記高周波磁場パルスは一定の繰り返し時間で印加されることを特徴とする磁気共鳴撮像装置。
【請求項１３】
請求項12に記載の磁気共鳴撮像装置において、
前記信号取得シーケンスは、核磁気共鳴信号の取得を伴わない高周波磁場パルスの印加を含むことを特徴とする磁気共鳴撮像装置。
【請求項１４】
請求項1に記載の磁気共鳴撮像装置において、
前記演算部は、前記経過時間の異なる複数の画像の信号強度を用いて、連立方程式を解くことにより、前記照射磁場分布を求めることを特徴とする磁気共鳴撮像装置。
【請求項１５】
請求項1に記載の磁気共鳴撮像装置において、
前記信号取得シーケンスは、k空間を放射状に走査して前記核磁気共鳴信号を取得するパルスシーケンスであることを特徴とする磁気共鳴撮像装置。
【請求項１６】
請求項1に記載の磁気共鳴撮像装置において、
前記演算部は、算出した前記照射磁場分布に基づき、前記RF照射部に供給される高周波磁場の位相と振幅の少なくとも一方を調整することを特徴とする磁気共鳴撮像装置。
【請求項１７】
請求項1に記載の磁気共鳴撮像装置において、
前記RF照射部は、複数のチャンネルを有してなり、
前記撮像部は、前記チャンネル毎に前記B1分布計測シーケンスを実行して前記被検体を撮像し、
前記演算部は、前記チェンネル毎に前記照射磁場分布を求めることを特徴とする磁気共鳴撮像装置。
【請求項１８】
RF照射部を備えた磁気共鳴撮像装置における前記RF照射部の前記照射磁場分布を計測する方法であって、
前記RF照射部による高周波磁場プリパルスの印加から経過時間(TI)を設けて核磁気共鳴信号を取得する信号取得シーケンスを実行する計測ステップと、
前記核磁気共鳴信号を用いて前記被検体の画像を再構成する画像再構成ステップと、
前記画像を用いて前記RF照射部の照射磁場分布を求める照射磁場分布算出ステップと、
を備え、
前記照射磁場分布算出ステップは、前記経過時間の異なる複数の画像を用いて前記照射磁場分布を求めることを特徴とする照射磁場計測方法。
【請求項１９】
請求項18に記載の照射磁場計測方法であって、
前記計測ステップは、前記高周波磁場プリパルスの印加から前記経過時間が異なる複数の信号取得シーケンスを実行し、
前記画像再構成ステップは、前記経過時間が異なる複数の信号取得シーケンスでそれぞれ取得された核磁気共鳴信号を用いて、前記経過時間の異なる複数の画像を再構成することを特徴とする照射磁場計測方法。
【請求項２０】
請求項19に記載の照射磁場計測方法であって、
前記計測ステップは、前記高周波磁場プリパルスを用いずに信号取得シーケンスのみを実行するステップと、前記高周波磁場プリパルスの印加とそれに続く少なくとも2回の信号取得シーケンスを実行するステップと、を含み、
前記照射磁場分布算出ステップは、各信号取得シーケンスでそれぞれ得られた核磁気共鳴信号を再構成して得られた複数の画像の信号強度の比を用いて、行列式を解くことにより、前記照射磁場分布を算出することを特徴とする照射磁場計測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、磁気共鳴撮像装置(以下、MRI装置という)に係り、特に被検体に高周波磁場を照射する照射コイルの照射磁場分布を計測する機能を備えたMRI装置に関する。
【背景技術】
【０００２】
MRI装置は、均一な静磁場中に検査対象を配置した状態で、検査対象に高周波磁場パルスを印加することによって生じた核磁気共鳴(NMR)信号を計測し、NMR信号の演算により検査対象の画像を再構成する装置である。検査対象を配置する静磁場として高磁場の磁場発生装置を用いることにより、高いSNの画像を得ることができる。【０００３】
近年、超伝導磁石の開発に伴い、3T以上の高磁場を実現できる高磁場MRI装置が普及している。高磁場MRIでは、高いSNが得られるが、腹部撮像等で画像に輝度ムラが発生する問題がある。この輝度ムラの原因の一つとして、被検体の組織内の原子核スピンを励起する高周波磁場パルス(以下、RFパルスという)の磁場分布(照射磁場分布、B1分布)の空間的な不均一がある。一般に、励起のための高周波磁場の共鳴周波数は、静磁場強度に比例するため、高磁場MRIでは、従前の高周波磁場よりも高い周波数の磁場を照射する必要がある。その場合、生体内での高周波磁場の波長が生体(特に腹部)の大きさと同じスケールとなる。このため、高周波磁場の位相が生体内の位置によって変化し、画像ムラとなって現れる。【０００４】
この照射磁場分布(B1分布)の空間的不均一を解決するための技術としてRFシミングがある。RFシミングでは、複数のチャンネルを持つ送信用RFコイルを用い、それぞれのチャンネルに与えるRFパルスの強度と位相を独立に制御することでB1分布の不均一を低減する。各チャンネルに与えるRFパルスの強度と位相を決定するためには、被検体毎、撮像部位毎に各チャンネルのB1分布が必要であり、B1分布の計測方法が種々提案されている。【０００５】
B1分布を計測(B1計測)する一般的な方法は、Double Angle法(DAM)と呼ばれる手法であり、任意のフリップ角のRFパルスを用いて撮像した画像と、その2倍のフリップ角のRFパルスとを用いて撮像した画像との演算によりB1を計測する。また複数の異なるフリップ角のRFパルスを用いて複数の画像を取得し、パルスシーケンスにより決定される信号強度式により、複数の画像の信号強度(画素値)をフィッティングすることでB1分布を計算する手法も提案されている(非特許文献1)。また同様に取得した複数の画像に対し、フィッティングを用いないで、信号強度変化の周期からB1分布を計算する手法も提案されている(特許文献1)。【０００６】
さらに、高周波磁場プリパルス(以下、単にプリパルスと略記する)を印加して画像を取得するパルスシーケンスを、プリパルスの強度を変化させて繰り返し、取得した画像からB1分布を計算する手法(非特許文献2)や、プリパルスを印加直後に取得した画像と、プリパルスを印加せずに取得した画像との比を取ることでB1分布を計算する手法(非特許文献3)も提案されている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開2008−68830号公報
【非特許文献】
【０００８】
【非特許文献１】Hai-King Margaret Cheng,Graham A Wright 著、「Rapid High-resoluTIon T1 Mapping by Variable Flip Angles:Accurate and Precise Measurements in the Presence of Radiofrequency Field Inhomogeneity」、MagneTIc Resonance in Medicine 55:566-574
【非特許文献２】J.T.Vaughan,M.Garwood,C.M.Collins,W.Liu,L DelaBarre,G.Adriany,P.Andersen,H.Merkle,R.Goebel,M.B.Smith,K.Ugurbil 著、「7T vs 4T:RF Power,Homogeneity,and Signal-to-Noise Comparison in Head Images」、MagneTIc Resonance in Medicine 46:24-30(2001)
【非特許文献３】H-P.Fautz,M.Vogel,P.Gross,A.Kerr,and Y.Zur 著、「B1 mapping of coil arrays for parallel transmission」、Proc.Intl.Soc.Mag.Reson.Med.16(2008) 1247
【非特許文献４】R.Lattanzi,C.Glaser,A.V.Mikheev,C.Petchprapa,D.J.Mossa,S.Gyftopoulos,H.Rusinek,M.Recht,and D.Kim 著、「A B1-insensiTIve High ResoluTIon,2D T1 Mapping Pulse Sequence for Radial dGEMRIC of the Hip at 3T」,Prc.Intl.Soc.Mag.Reson.Med.19(2011) 504
【非特許文献５】John G.Sled,G.Bruce Pike:MagneTIc Resonance in Medicine 43:589-593 (2000)
【発明の概要】
【発明が解決しようとする課題】
【０００９】
しかしながら、上述した各手法には次のような問題がある。Double Angle法では、T1緩和の影響を除去するために、繰り返し時間TR(RFパルスの印加間隔)を一般的には5秒程度の時間に設定する。これにより、単純な計算式によりB1分布の計算が可能であるが、長い撮像時間を要し、撮像時間が約10分以上となることもある。特許文献1及び非特許文献1の技術では、T1値がフィッティングに使用する関数に含まれるため、その影響を除去する必要がなく、比較的短いTRでの撮像が可能であり、撮像時間はDouble Angle法より短くすることが可能である。しかしフィッティングや信号強度変化の周期を見出すために約20種類のフリップ角で画像を取得する必要があり、依然として撮像時間が長いという問題がある。【００１０】
また計算精度が撮像枚数やフィッティングの精度に依存するため、高精度にB1分布を計算するためにはさらに長い撮像時間が必要となる。【００１１】
非特許文献2の技術では、B1分布計算の基礎となる画像信号の強度変化の周期は、プリパルスのフリップ角のみに依存するため、画像信号を取得するためのパルスシーケンスによっては影響されない。従って、このパルスシーケンスにおけるTRを短く設定することが可能であり、Double Angle法に比べ撮像時間を短くすることが可能である。しかし、プリパルスの強度の変化による画像信号の強度の変化の周期を、フィッティングを用いて計算するため、高精度にB1分布を計算するためには、プリパルスのフリップ角を20回以上変化させて撮像する必要があり、依然として撮像時間が長いという問題がある。またB1分布計算には煩雑な計算が必要となる。【００１２】
非特許文献3の技術では、プリパルスの印加からエコー信号(NMR信号)取得までの時間(TI)を可能な限り短くすることで、高速にB1分布を計算することが可能になるが、計算式がTI＝0でしか使えない近似を用いているため、TIを長くした場合、誤差が大きくなるという問題がある。【００１３】
そこで本発明は、短時間で高精度にB1分布を計測することが可能なMRI装置を提供することを目的とする。
【課題を解決するための手段】
【００１４】
本発明は、B1分布計測に、プリパルスの印加と信号取得パルスシーケンスとを組み合わせるとともに、プリパルス印加からの経過時間(TI)が異なる複数の信号取得パルスシーケンスを実行し、これら複数の信号取得パルスシーケンスの実行によって取得したTIの異なる複数の画像(画像用データ)の演算によって、B1を求めるという新規な手法を提案するものである。【００１５】
即ち、本発明のMRI装置は、被検体に核磁気共鳴を起こさせるための高周波磁場(B1)を照射するRF照射部と、プリパルスの印加から経過時間(TI)を設けてエコー信号を取得する信号取得シーケンスを有して成るB1分布計測シーケンスを用いて被検体を撮像する撮像部と、エコー信号を用いて被検体の画像を再構成する演算部と、を備え、演算部は、画像を用いてRF照射部の照射磁場分布を求めるMRI装置であって、演算部は、経過時間の異なる複数の画像を用いて照射磁場分布を求めることを特徴とする。【００１６】
本発明の照射磁場計測方法は、RF照射部によるプリパルスの印加から経過時間(TI)を設けてエコー信号を取得する信号取得シーケンスを実行する計測ステップと、エコー信号を用いて被検体の画像を再構成する画像再構成ステップと、画像を用いてRF照射部の照射磁場分布を求める照射磁場分布算出ステップと、を備え、照射磁場分布算出ステップは、経過時間の異なる複数の画像を用いて照射磁場分布を求めることを特徴とする。
【発明の効果】
【００１７】
本発明によれば、プリパルスと、プリパルス印加からの経過時間が異なる複数の信号取得シーケンスとの組合せにより、極めて短時間で且つ精度の高いRFパルスの磁場分布計測(B1計測)を行うことができる。計測した磁場分布を用いて、RF照射手段を制御することにより、高精度なRFシミングが実現できる。【００１８】
また演算手段が、プリパルス無の画像の信号強度とプリパルス有の画像の信号強度との比を用いて、行列式を解くことにより、照射磁場分布を算出した場合には、B1分布と併せて、T1分布の算出を行うことができる。
【図面の簡単な説明】
【００１９】
【図１】本発明が適用されるMRI装置の一実施形態を示すブロック図
【図２】本発明によるB1分布計測手順の一実施形態を示すフロー図
【図３】B1分布計測シーケンスの第一の実施形態を示す図
【図４】信号取得シーケンスの一例を示すタイムチャート
【図５】信号取得シーケンスの他の例を示すタイムチャート
【図６】信号取得シーケンスのさらに別の例を示すタイムチャート
【図７】(a)は、B1分布計測シーケンスの第二の実施形態を示す図、(b)は、(a)の変更例を示す図
【図８】B1分布計測シーケンスの第二の実施形態を一般化した図
【図９】B1分布計測シーケンスの第三の実施形態を示す図
【図１０】第三実施形態の変更例を示す図
【図１１】B1分布計測シーケンスの第四の実施形態を示す図
【図１２】B1分布計測シーケンスの第五の実施形態を示す図
【図１３】(a)、(b)は、それぞれ、B1分布計測シーケンスの第六の実施形態を示す図
【図１４】本発明によるB1分布計測手順の他の実施形態を示すフロー図
【図１５】実施例1のB1分布算出結果を示す図
【図１６】実施例2のB1分布算出結果を示す図
【図１７】従来法(DAM)によるB1分布算出結果を示す図
【図１８】実施例3のB1分布算出結果を示す図
【図１９】本発明の手法と従来法(DAM)との計算精度を比較した結果を示す図
【発明を実施するための形態】
【００２０】
以下、本発明の実施の形態を説明する。まず本発明が適用されるMRI装置の全体構成について、図面を参照して説明する。
【００２１】
図1は、本発明が適用されるMRI装置の一実施形態を示すブロック図である。このMRI装置は、静磁場発生系2と、傾斜磁場発生系3と、送信系5と、受信系6と、信号処理系7と、シーケンサ4と、中央処理装置(CPU)8とを備えて構成される。
【００２２】
静磁場発生系2は、被検体1が置かれる空間に均一な静磁場を発生するものであり、永久磁石方式、常電導方式あるいは超電導方式の静磁場発生源(不図示)からなる。静磁場発生源は、垂直磁場方式であれば、被検体1の体軸と直交する方向に、水平磁場方式であれば、体軸方向に、均一な静磁場を発生させるように、配置されている。
【００２３】
傾斜磁場発生系3は、MRI装置の座標系(静止座標系)であるX、Y、Zの直交3軸方向に傾斜磁場を印加する傾斜磁場コイル9と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源10とから成る。後述のシ−ケンサ4からの命令に従ってそれぞれのコイルの傾斜磁場電源10を駆動することにより、X、Y、Zの3軸方向に所望の傾斜磁場Gx，Gy，Gzを印加することができる。傾斜磁場の印加の仕方によって、被検体の撮像スライスを選択的に励起し、また励起領域から発生するエコー信号(NMR信号)に位置情報を加えることができる。
【００２４】
シーケンサ4は、RFパルスと傾斜磁場パルスをある所定のパルスシーケンスで繰り返し印加する制御手段で、CPU8の制御で動作し、被検体1の断層画像のデータ収集に必要な種々の命令を送信系5、傾斜磁場発生系3、および受信系6に送る。
【００２５】
送信系5は、被検体1の生体組織を構成する原子の原子核スピンに核磁気共鳴を起こさせるために、被検体1にRFパルスを照射するもので、高周波発振器11と変調器12と高周波増幅器13と送信側の高周波コイル(送信コイル)14aとから成る。送信コイルは、本実施形態では、複数の給電点を有し、供給される高周波の強度と位相を調整できるように構成されている。高周波発振器11、変調器12及び高周波増幅器13は、各チャンネルに対応して複数備えられている。図では、2つの給電点がある場合を示しているが、給電点の数は2に限定されない。
【００２６】
高周波発振器11から出力されたRFパルスをシーケンサ4からの指令によるタイミングで変調器12により振幅変調し、この振幅変調されたRFパルスを高周波増幅器13で増幅した後に被検体1に近接して配置された高周波コイル14aに供給することにより、RFパルスが被検体1に照射される。シーケンサ4からのタイミングと変調器12による変調は、後述するB1分布の計測結果を反映して制御される。
【００２７】
受信系6は、被検体1の生体組織を構成する原子核スピンの核磁気共鳴により放出されるエコー信号を検出するもので、受信側の高周波コイル(受信コイル)14bと信号増幅器15と直交位相検波器16と、A/D変換器17とから成る。送信コイル14aから照射された電磁波によって誘起された被検体1の応答のNMR信号が被検体1に近接して配置された受信コイル14bで検出され、信号増幅器15で増幅された後、シーケンサ4からの指令によるタイミングで直交位相検波器16により直交する二系統の信号に分割され、それぞれがA/D変換器17でディジタル量に変換されて、信号処理系7に送られる。
【００２８】
なお図1では、送信用の高周波コイルと受信用の高周波コイルが、別個に設けられている構成を示しているが、一つの高周波コイル(マルチプルコイルを含む)が送信用及び受信用を兼ねる構成とすることも可能である。
【００２９】
信号処理系7は、CPU8と、各種データ処理と処理結果の表示及び保存等を行うもので、光ディスク19、磁気ディスク18等の外部記憶装置と、CRT等からなるディスプレイ20とを有する。受信系6からのデータがCPU8に入力されると、CPU8が信号処理、画像再構成等の処理を実行し、その結果である被検体1の断層画像をディスプレイ20に表示すると共に、外部記憶装置の磁気ディスク18等に記録する。
【００３０】
CPU8は、信号処理系7の演算部としての機能のほかに、装置の各要素を制御する制御部としての機能を有し、シーケンサ4を介して、種々のパルスシーケンスを実行させる。パルスシーケンスは、予めプログラムとして組み込まれている。本実施形態では、送信コイルによる照射磁場分布(B1分布)を計測するためのB1分布計測シーケンスを備えている。また信号処理系7は、このB1分布計測シーケンスの計測結果を用いて、B1分布の計算や送信コイルに与えられる高周波パルスの位相や振幅の計算を行い、この計算結果に基づき、送信コイルに与えられる高周波パルスの位相や振幅を制御する。
【００３１】
操作部25は、MRI装置の各種制御情報や上記信号処理系7で行う処理の制御情報を入力するもので、トラックボール又はマウス23、及び、キーボード24から成る。この操作部25はディスプレイ20に近接して配置され、操作者がディスプレイ20を見ながら操作部25を通してインタラクティブにMRI装置の各種処理を制御する。
【００３２】
なお、図1において、送信側の高周波コイル14aと傾斜磁場コイル9は、被検体1が挿入される静磁場発生系2の静磁場空間内に、垂直磁場方式であれば被検体1に対向して、水平磁場方式であれば被検体1を取り囲むようにして設置されている。また、受信側の高周波コイル14bは、被検体1に対向して、或いは取り囲むように設置されている。
【００３３】
次に上述したMRI装置において、B1分布計測を行う手順を説明する。図2に手順を示す。
【００３４】
まず静磁場空間内に被検体を配置し、目的とする撮像部位が静磁場空間のほぼ中央に位置するように位置決めを行う(ステップ201)。次いで、送信コイルのチャンネル毎に、B1分布計測シーケンスを実行し画像を取得する(ステップ202)。B1分布計測シーケンスについては、後述する。ステップ202で取得した画像を用いてチャンネル毎のB1分布を算出する(ステップ203)。算出したB1分布を用いて、そのチャンネルに印加する高周波パルスの振幅と位相を設定する(ステップ204、205)。その後、ステップ204で設定した条件で、所望の計測(撮像)を行う(ステップ206)。以下、B1分布計測の典型的な実施形態について、詳細に説明する。以下説明する各実施形態において、ステップ201とステップ204以降は共通であり、ステップ202とステップ203について、実施形態毎に説明する。
【００３５】
＜第一の実施形態＞
＜＜B1分布計測(ステップ202)＞＞
B1分布計測シーケンスは、比較的大きなフリップ角のRFパルスからなるプリパルスと、小さなフリップ角のRFパルスを用いた信号取得シーケンスとの組合せからなり、信号取得シーケンスでは、プリパルス印加後の経過時間(TI)が異なる複数の画像を得る。プリパルスと信号取得シーケンスの組み合わせの仕方によって、異なる態様が可能であり、その態様によってB1分布算出(ステップ203)が異なる。
【００３６】
本実施形態は、プリパルスの印加に引き続き、少なくとも3回の信号取得シーケンスを実行し、各信号取得シーケンスで得た画像から、B1分布を算出することが特徴である。
【００３７】
図3は、プリパルス301と信号取得シーケンス303、305、307との関係を示す図で、プリパルス301印加からの経過時間(TI)に依存した信号強度の変化を図中上側にグラフで示している。図示するグラフにおいて、横軸はプリパルス印加後の経過時間、縦軸は信号強度である。プリパルス301は、例えば非選択性のRFパルスで、フリップ角の大きい、例えば90度パルスである。このプリパルス301によって励起された原子核スピンが縦緩和している間に、少なくとも3回の信号取得シーケンス303、305、307を実行し、TIの異なる3つのk空間データ(又は画像データ)を取得する。なお、本発明では、プリパルスからの縦緩和が異なるデータを得ることが重要であり、プリパルスの影響が十分に残っている間にすべての信号取得シーケンスが完了するようにTIを設定する。これにより高精度にB1分布の計算を行うことができる。
【００３８】
各信号取得シーケンスは、短時間でk空間データを収集できるパルスシーケンスであれば特に限定されず、例えば、図4に示すようなグラディエントエコー(GrE)系のパルスシーケンスを採用できる。このGrE系シーケンスでは、スライス傾斜磁場パルス402とともに低フリップ角のRFパルス401を印加した後、位相エンコード傾斜磁場403を印加する。同時に読出し傾斜磁場404を印加し、極性の反転した読出し傾斜磁場の印加中にエコー信号405を計測する。最後に位相エンコード方向にリフェイズのための傾斜磁場406を印加する。RFパルス401として、縦磁化に対する影響を少なくするため、好ましくは10度以下、より好ましくは5度以下の低フリップ角パルスを用いる。
【００３９】
このパルスシーケンスは、RFパルス401として低フリップ角パルスを用いると共に、位相エンコード方向のリフェイズパルス406を用いているので、繰り返し時間TRを数ms(ミリ秒)程度にすることができる。RFパルス401からリフェイズ傾斜磁場406までを、位相エンコード傾斜磁場パルス403の強度を変えながら繰り返し、スライス傾斜磁場402によって選択されたスライスのデータ(k空間データ)を得る。
【００４０】
このk空間データは、後述するB1分布の計算に用いるものであり、マトリクスサイズは64×64程度でよい。これにより、極めて短時間、具体的には200ms程度の計測時間で全k空間データを取得することができる。
【００４１】
k空間データから形成される画像データのコントラストは、主としてk空間データの中央のデータにより支配されるので、各信号取得シーケンスにおいてk空間データの中心のエコーを計測するタイミングをプリパルス印加からの経過時間(TI)とする。図3に示す実施形態では、各信号取得シーケンスはk空間の中央のエコーから計測を開始する所謂セントリックオーダーのパルスシーケンスであり、信号取得シーケンスの開始時点が、それぞれTI1、TI2、TI3である。またこれら経過時間の関係は、本実施形態では、TI2＝2×TI1、TI3＝3×TI1の関係に設定されている。
【００４２】
＜＜B1分布算出(ステップ203)＞＞
次に、図3に示すB1分布計測シーケンスによって得られたk空間データから、B1分布を算出する手法を説明する。
【００４３】
最初の信号取得シーケンスで得られたk空間データを逆フーリエ変換することにより得た画像データの、ある注目画素の信号強度S(B1,TI)は式(1)で与えられる。
2012060192.tif 000003
式(1)中、Sseqは、プリパルスの後の信号取得シーケンスによって決まる信号強度を表し、αは設定したプリパルスのフリップ角を表し、TIはプリパルス印加からk空間中心の信号を収集するまでの時間を表し、T1は組織に依存する縦緩和時間を表す。
【００４４】
同様に、2番目、3番目の信号取得シーケンスから得られた画像データの注目画素の信号強度は、式(2)、式(3)で表わすことができる。
2012060192.tif 000004
ここで
2012060192.tif 000005
と定義すると、式(1)〜式(3)は、式(4)〜式(6)のように書くことができる。
2012060192.tif 000006
式(4)〜式(6)の連立方程式を解くことにより、式(7)、(8)よりX及びYが求められる。
2012060192.tif 000007
ここで、定義より、X＝1−cos(B1・α)であるので、B1は式(9)で計算できる。
2012060192.tif 000008
本実施形態では、TIの異なる少なくとも3回の信号取得シーケンスで得た画像から連立方程式を解くことによってB1分布を求めることができるので、数秒という極めて短い時間でB1分布計測を行うことができる。
【００４５】
上述したように、B1分布はTIの異なる複数回の計測結果(信号強度)から連立方程式を解くことによって正確に求めることができ、演算の容易性の点でも好適であるが、複数回の計測結果をフィッティングすることにより、B1分布を算出することも可能である。
【００４６】
また本実施形態では、上述した連立方程式を解くために、少なくとも3つの画像(3回の信号取得シーケンスの実行)が必要であるが、プリパルス印加後のスピン縦緩和時間の範囲であれば、3回以上であってもよい。また複数の信号取得シーケンスのTIは上述した例では整数比としたが、それぞれ異なっていれば良く、整数比以外でもよい。
【００４７】
信号取得シーケンスについても、短時間で画像データを取得するものであれば、図4に示すシーケンスに限らず、種々の変更が可能である。
【００４８】
例えば、図3に示すTIの異なる3つの信号取得シーケンスを、一定の間隔をあけて実行するのではなく、間隔を設けずに連続して複数の信号取得シーケンスを実行することも可能である。或いは信号取得シーケンスの前後で、エコーを計測することなく、信号取得シーケンスと同じTRでRFパルス(図3の301)を連続して印加することも可能である。これらの変更例では、RFパルスを連続して印加することにより、スピンを定常状態に維持し、信号取得シーケンス内でのコントラスト差が生じるのを抑制することができる。
【００４９】
また図4では、一方向の位相エンコードパルスを用いるGrE系のパルスシーケンスを示したが、二軸方向のエンコードパルスを組み合わせてk空間を放射状に走査する所謂ラジアル法のパルスシーケンスを採用することもできる。
【００５０】
ラジアル法のパルスシーケンスの一例を図5に示す。この信号取得シーケンスでは、RFパルス501を、スライス傾斜磁場パルス502と共に印加し、次いで、2軸方向(Gp方向およびGr方向)にそれぞれ傾斜磁場パルス504、506と、極性の反転した傾斜磁場パルス505、507を印加し、傾斜磁場パルス504、506の印加量と傾斜磁場パルス505、507の印加量の絶対値が同一になった時点でピークとなるグラディエントエコー信号508を計測する。励起からエコー信号を発生させるまでの時間すなわちエコー時間TＥは画像コントラストを決めるパラメータであり、撮影の目的に応じて任意に設定できる。
【００５１】
計測した1つのエコー信号508をk空間に配置したデータは、k空間の原点を通り、k空間の座標軸に対する角度が2軸の傾斜磁場パルスの印加量によって決まる一列のデータ(スポーク)となる。以上のステップを繰り返し時間TR毎に、2軸方向の傾斜磁場パルスの印加量とその比を異ならせながら回繰り返すことにより、k空間上を回転しながらデータを取得することができる。このパルスシーケンスを採用する場合も、RFパルス501のフリップ角を10度以下とし、TRを数ms程度とすることにより、繰り返し回数を200回(スポーク数を200)とした場合でも、0.1秒のオーダーで一つのk空間データを取得することができる。
【００５２】
ラジアル法のパルスシーケンスでは、各スポークのデータがk空間の中心を通るので、式(1)で表わす信号強度は、全データのうち最初のスポーク取得時から最後のスポーク取得時までの縦緩和の影響を平均した値となるが、上述したとおり、全データの取得時間が短いので、図4の信号取得シーケンスを用いたときと同様に、式(1)〜式(9)までの計算により、B1分布を算出することができる。
【００５３】
別の変更例として、図4に示す信号取得シーケンスのスライス軸(Gｓ)にスライスエンコード傾斜磁場が加えた3Dパルスシーケンス(図6)を採用することも可能である。
3Dパルスシーケンスを用いることにより、三次元のB1分布を計測することが可能になる。
【００５４】
＜第二の実施形態＞
＜＜B1分布計測(ステップ202)＞＞
本実施形態では、プリパルス印加からの経過時間が異なる少なくとも2回の画像用信号取得シーケンス(第一のB1分布計測シーケンス)と、第一のB1分布計測シーケンスとは異なるプリパルス印加の条件の信号取得シーケンス(第二のB1分布計測シーケンス)とを実行することによって取得した、少なくとも3つの画像から、B1分布を算出することが特徴である。
【００５５】
図7(a)は、第一のB1分布計測シーケンス710と、第二のB1分布計測シーケンス720との関係を示す図である。第一のB1分布計測シーケンス710は、図3に示す第一の実施形態と同様に、フリップ角が90度以上のRFパルスからなるプリパルス711と、それに続く少なくとも2回の信号取得シーケンス713、715とからなる。ただし、第一の実施形態では、プリパルス711印加後のスピンの縦緩和時間以内に少なくとも3回の信号取得シーケンスを実行したのに対し、本実施形態では、信号取得シーケンス713、715は2回以上であればよい。
【００５６】
第二のB1分布計測シーケンス720は、プリパルスを用いない単一の信号取得シーケンス721からなる。図7(a)では、第一のB1分布計測シーケンス710の前に、第二のB1分布計測シーケンス720が配置されているが、図7(b)に示すように、第一のB1分布計測シーケンス710の後の縦緩和時間TD経過後に第二のB1分布計測シーケンス720が配置されていてもよい。
【００５７】
信号取得シーケンス713、715及び721は、第一の実施形態と同様に、短時間でk空間データを収集できるパルスシーケンスであれば、種類を問わない。例えば、図4や図5に示すようなGrE系のパルスシーケンスを採用できる。
【００５８】
＜＜B1分布算出(ステップ203)＞＞
次に、図7に示すB1計測パルスシーケンスによって得られたk空間データから、B1を算出する手法を説明する。
【００５９】
第二のB1分布計測シーケンス720(信号取得シーケンス721)はフリップ角が0度のプリパルスを用いた信号取得シーケンスとみなすことができ、取得した画像の、ある注目画素の信号強度は、式(1)においてα＝0°とした、式(10)で与えられる。
2012060192.tif 000009
一方、プリパルス711のフリップ角をαとし、プリパルス711印加からの経過時間TI、TI×2において、k空間の中心データを取得した最初の信号取得シーケンス713、715では、注目画素の信号強度は式(11)、式(12)で与えられる。
2012060192.tif 000010
式(11)及び式(12)を、式(13)、式(14)のように変形し、
2012060192.tif 000011
式(13)を二乗し、式(14)で除算することにより、式(15)によりB1分布が求められる。
2012060192.tif 000012
なお図7に示す例では、第二のB1分布計測シーケンスはプリパルスを用いない信号取得シーケンスで構成されているが、本実施形態は、第一および第二のB1分布計測シーケンスにおけるプリパルスの印加条件(フリップ角)が異なっており、少なくとも3つの信号取得シーケンスが実行されればよく、例えば図8に示すように、第一のB1分布計測シーケンス810を、第1のフリップ角のプリパルス811の印加と2つの信号取得シーケンス813、815で、第二のB1分布計測シーケンス820を、第1のフリップ角とは異なる第2のフリップ角のプリパルス821の印加と1つの信号取得シーケンス823で構成することも可能である。第1のフリップ角及び第2のフリップ角との関係は特に限定されないが、例えば、90度と0度、90度と180度などの組み合わせとすることができる。
【００６０】
また図において、第一のB1分布計測シーケンス810と第二のB1分布計測シーケンス820とは、順序を入れ替えてもよい。いずれの場合にも、第一の計測と第二の計測との間(プリパルスの印加間隔TD)は、縦緩和時間以上に設定することが望ましく、それによりB1分布計測の精度を向上することができる。
【００６１】
本実施形態によれば、非特許文献2記載の技術に比べ、フリップ角の変更を1回で済ませられるため、撮像時間を［2÷(非特許文献2のフリップ角変化回数)］に短縮することができる。また第一実施形態と比較した場合、フリップ角のデータも用いているため、より高精度なB1分布を計測することができる。
【００６２】
なおB1分布計測シーケンスを繰り返す場合に、その間隔(プリパルス印加間隔TD)を短くすると、十分に緩和しない状態で次のプリパルスを印加することとなり、式(10)〜式(15)を用いた計算の誤差が大きくなる可能性がある。この問題は、上述したようにTDを縦緩和時間以上に設定することで解決できるが、B1分布計測シーケンスにおける最後の信号取得シーケンスの後に、強制的に縦緩和するRFパルスを追加し、プリパルス印加間隔を短縮することも可能である。強制的に縦緩和させるRFパルスとしては、プリパルスと強度が同じで位相が逆のRFパルスを用いることができ、このような強制縦緩和パルスを用いることで、撮像時間を短縮しても、精度の高いB1分布の計算を行うことができる。
【００６３】
＜第三の実施形態＞
＜＜B1分布計測(ステップ202)＞＞
本実施形態は、プリパルスを印加しないで信号を取得する基準画像取得シーケンスと、プリパルス印加後に、プリパルス印加からk空間中心の信号を取得するまでの時間TIが異なる複数の信号取得シーケンスを実行し、基準画像とTIの異なる複数の画像を用いた行列演算によってB1分布を算出することが特徴である。
【００６４】
図9に、本実施形態におけるプリパルスと、信号取得シーケンスとの関係を示す。図9は、プリパルスと信号取得シーケンスとの関係を示す図7に類似しているが、本実施形態では、信号取得シーケンスは、プリパルスの影響が大きい間に信号を取得する、短TIのシーケンスであり、第一及び第二実施形態では、最初の信号取得シーケンスのTIとそれに続く信号取得シーケンスのTIとの比が整数比であったのに対し、本実施形態では整数比に限らない。
【００６５】
信号取得シーケンス903、905、907・・・は、第一実施形態と同様に、図4に示すような短TRのグラディエントエコー系のパルスシーケンスとすることができ、セントリックオーダーであること、またフリップ角を小さく設定することが好ましい。プリパルス901は、非選択性のフリップ角の大きいRFパルスである。基準画像取得シーケンス900は、信号取得シーケンス903等と同じパルスシーケンスであり、撮像時間を短縮するために、プリパルス901を印加する直前に実行することが好ましい。
【００６６】
＜＜B1分布算出(ステップ203)＞＞
次に基準画像取得シーケンス900及び信号取得シーケンスで取得した画像を用いたB1分布の算出について説明する。
【００６７】
プリパルス印加後にk番目(k＝1,2,3・・・n)の信号取得シーケンスで取得した信号から再構成した画像の注目画素の信号強度は、k番目のTIをTIkとすると、第一実施形態で示した式(1)と同様の式(16)で与えられる。
2012060192.tif 000013
式中、式(1)と同じ符号は同じ意味を持つ。
【００６８】
一方、プリパルスの直前の基準画像取得シーケンスによって得られた画像の同じ注目画素に信号強度は、式(16)においてα＝0とした場合と同じであるため、式(10)と同一の式(17)で与えられる。
2012060192.tif 000014
式(16)を式(17)で除算し、自然対数を取ると、式(18)のように、log(1-cos(B1・α))と(−TIk/T1)の線形結合で表わすことができる。
2012060192.tif 000015
各信号取得シーケンスから得られたTIの異なる画像について同様の計算をすると、式(19)の連立方程式が得られる。
2012060192.tif 000016
ここで、Sは1×nの行列、Aは2×nの行列、Xは1×2の行列である。Wi(i＝1，2，3・・・n)はそれぞれのTIに対する重みを示し、任意に設定することができる。行列Aの擬似逆行列pinvAを左から掛けることで、式(19)を解くことができ、次式(20)、(21)のようにB1及びT1を求めることができる。
【００６９】
【数１】
2012060192.tif 000017
第一及び第二実施形態は、式(4)〜(6)或いは式(13)及び(14)の連立方程式を解くことにより、解(B1)を求めるので、解が発散しないためには、プリパルス後の信号取得シーケンスのTIを整数比としていたが、本実施形態では、そのようなTIの制限が要せずに解を求めることができるので、プリパルスの効果が高い、短TIで信号取得でき、計測時間を短縮できるとともに高精度なB1分布を得ることができる。また本実施形態では、B1分布と併せてT1分布を算出することが可能である。算出されたT1分布は、例えば、非特許文献4に記載されているように、骨関節炎の診断等に利用される。
【００７０】
なお図9では、プリパルス901の直前に、基準画像取得シーケンス900を実行し、基準画像を取得する場合を説明したが、プリパルス901印加後に十分の時間即ちT1緩和するに十分な時間経過後に、基準信号取得シーケンス900を実行してもよい。図10に、この変更例を示す。図示するように、プリパルス901印加後に異なるTIで信号取得シーケンス903、905、907・・・を実行してTIの異なる画像を取得し、最後にTIを十分長く設定して基準画像取得シーケンス900を実行して基準画像を取得する。この基準画像の信号強度S0は、式(17)で示したように、S0＝Sseqと表すことができるので、図9の場合と同様に、式(18)〜式(21)により、B1分布及びT1分布を算出することができる。
【００７１】
＜第四の実施形態＞
本実施形態は、複数の断面のB1分布計測をマルチスライスの手法で行うことを特徴とする。それ以外は、第一〜第三の実施形態と同様であり、異なる点を中心に説明する。
【００７２】
図11に本実施形態によるB1分布計測シーケンスの一例を示す。図示するように、B1分布計測シーケンス910は、プリパルス911とそれに続く信号取得シーケンス913、915、917からなる。このB1分布計測シーケンス910は、図3に示すB1分布計測シーケンスと類似しているが、各信号取得シーケンス913、915、917は、互いに直交する3つの断面のうちの一つを、それぞれ選択して行う。
【００７３】
即ち、例えば図4に示すシーケンスにおいて、RFパルス401と同時に印加されるスライス選択傾斜磁場402は、信号取得シーケンス毎に異なる軸のスライス傾斜磁場であり、異なる断面が励起される。これにより、各断面の画像は、プリパルス911印加からの経過時間TIが異なることになる。所定の間隔(TD)を置いて、同様のB1分布計測シーケンス920、930を繰り返す。その際、各B1分布計測シーケンスに含まれる信号取得シーケンスにおいて選択されるスライス(断面)の順序をサイクリックに変更する。例えば、最初のB1分布計測シーケンス910では、アキシャル面(AX)、サジタル面(SG)、コロナル面(COR)とし、次のB1分布計測シーケンス920では、SG→COR→AXの順、その次のB1分布計測シーケンス830では、COR→AX→SGの順とする。これによって、3回目のB1分布計測シーケンス930が終了した時点では、アキシャル面(AX)、サジタル面(SG)、コロナル面(COR)のそれぞれについて、TIの異なる3つの画像データ(TI1、TI2、TI3のk空間データ)を得ることができる。
【００７４】
図9は、複数の断面が互いに直交する3つの断面である場合を示したが、複数の断面が互いに平行な断面の場合にも、各信号取得シーケンスにおけるスライス選択条件を変えることにより、直交3断面の場合と同様に、複数の平行な断面について3組あるいはそれ以上の組の画像データを得ることができる。
【００７５】
その後、これら3組のデータを用いて、式(1)〜式(9)によりB1分布を算出することは第一の実施形態と同様である。
【００７６】
図9に示すB1分布計測シーケンスに代えて、図7(a)又は(b)に示すB1分布計測シーケンスとしてもよく、その場合にも、第二の実施形態と同様に、式(10)〜式(15)を用いて、各断面についてB1分布を算出することができる。
【００７７】
またプリパルスを印加せずに基準画像を得ることにより、第三の実施形態と同様に式(18)〜(21)を用いてB1分布を算出することができる。
【００７８】
本実施形態は、複数の断面のB1分布が必要となるRFシミングに有効である。本実施形態によれば、第一〜第三の実施形態と比較して、断面の数だけ計測時間が多くなるが、1回のB1分布計測シーケンスによる計測時間がわずか数秒程度であるので、10〜30秒程度の短い時間で3断面についてB1分布を計測することができる。
【００７９】
＜第五の実施形態＞
本実施形態は、B1分布の計算に用いるk空間データを複数回に分けて計測する。即ちマルチショットによる計測を行うことが特徴である。それ以外は、第一の実施形態〜第三の実施形態と同様であり、異なる点を中心に説明する。
【００８０】
本実施形態では、図12に示すように、プリパルス101と複数の信号取得シーケンス103、105、107とからなるB1分布計測シーケンス100を複数回行う。各回で同一フリップ角のプリパルス101を用いるが、1回の信号取得シーケンスではk空間データの一部を取得し、複数回で全k空間データを取得するようにする。k空間の分割の仕方は、特に限定されず、図示するように単純に複数の領域に分割してもよいし、それぞれの信号取得シーケンス110、120、130・・・で位相エンコードステップを粗くして計測し、k空間の中心に最も近いデータを常に所定のTIで計測するようにしてもよい。
【００８１】
第一〜第三の実施形態において、1回で全k空間データを取得する場合には、その中心データを取得するタイミングが、所定のTI(プリパルス印加からの経過時間)となるように設定されることを説明したが、その場合、k空間データにはTIが異なるデータが含まれることになる。本実施形態では、k空間データを分けて取得するので、k空間データに含まれる、TIが所定のTIと異なるデータの割合を少なくすることができる。その結果、信号取得時間はショット数倍延長することになるが、B1分布計測の精度を向上することができる。
【００８２】
本実施形態は、特に、信号取得シーケンスとして、図5に示すようなラジアルシーケンスや図6に示すような3Dパルスシーケンスを採用する場合に、有効である。
【００８３】
＜第六の実施形態＞
第一、第二の実施形態では、1回のプリパルスの後に続けてTIの異なる複数の信号取得シーケンスを実行する場合を説明したが、本実施形態は、1回のプリパルス後に1回の信号取得シーケンスを行い、この1回のプリパルスと1回の信号取得シーケンスからなるB1分布計測シーケンスを、信号取得シーケンスのTIを異ならせて繰り返し、TIの異なる複数の画像データを得ることが特徴である。
【００８４】
図13に、本実施形態のB1分布計測シーケンスを示す。図13(a)は、第一の実施形態のB1分布計測シーケンスを本実施形態に従い変更した場合、図13(b)は、第二の実施形態のB1分布計測シーケンスを本実施形態に従い変更した場合をそれぞれ示している。図13(a)において、3回のB1分布計測シーケンス1110、1120、1130で用いるプリパルス111及び信号取得シーケンス113、115、117は、3回とも同じであるが、TIのみが異なっている。また図13(b)において、3回のB1分布計測シーケンス1140、1150、1160で用いる信号取得シーケンス114、116、118は同じであるが、1回目のB1分布計測シーケンス1140は、プリパルスを用いてない(フリップ角0に相当)。また2回目のB1分布計測シーケンス1150と3回目のB1分布計測シーケンス1160は、TIのみが異なっている。
【００８５】
前述したように、信号取得シーケンスで用いるRFパルスは、フリップ角が10度以下の低フリップ角パルスであるが、連続して実行した場合には、後続の信号取得シーケンスでは先行する信号取得シーケンスで印加したRFパルスの影響を受ける。本実施形態では、TIの異なる信号取得シーケンスを別個に実行するので、互いのRFパルスによる影響を受けることなく、正確なTIを反映した画像データを得ることができる。またB1分布計測シーケンスの繰り返し回数は3回程度であるので、繰り返しによる計測時間の延長は少なく、図13(a)、(b)のいずれの場合も15秒程度でB1分布に必要なデータの収集を完了することができる。
【００８６】
以上、図2に示す本発明のMRI装置の動作のうち、コイルのチャンネル毎に行うB1計測ステップ202およびB1分布計算ステップ203の典型的な実施形態を説明したが、これら実施形態は適宜組み合わせることが可能であり、またパルスシーケンスやその繰り返し回数などは例示したものに限定されず種々の変更が可能である。
【００８７】
＜＜RFパルスの振幅および位相の設定(ステップ204、205)＞＞
次に、ステップ203で算出されたB1分布を用いたRFパルスの調整(図2のステップ204、205)を説明する。ここでは、RFコイルが複数の小型RFコイルからなるマルチプルアレイコイルであり、小型RFコイル毎に独立して制御できる給電点(チャンネル)を持つ場合を説明する。
【００８８】
RFコイルのチャンネル数がＮとして、各チャンネル毎に求めたB1分布をB1n(r)とし、各小型RFコイルに供給される高周波信号の振幅および位相をAn、φnとしたとき、全体としての磁場分布B1total(r)は式(22)で表わすことができる。
2012060192.tif 000018
式(22)の振幅及び位相の組(An, φn)を変化させて、磁場分布B1total(r)(rは実空間座標の位置)として、均一な磁場分布B1(r)を与える振幅及び位相の組を求める。この計算は、公知の非線形最適化アルゴリズムを用いて解くことができ、例えば、式(22)で求めたB1total(r)と目標とする磁場分布との平均2乗誤差の平方根を最小化する最適化アルゴリズムを用いて振幅及び位相の組(An,φn)を求めることができる。
【００８９】
求めた振幅及び位相の組を各小型RFコイルに設定する。具体的には、シーケンサ4および変調器12により、RFコイルの各チャンネルに供給される高周波パルスの振幅とタイミングとが調整される。
【００９０】
ステップ204、205で設定した振幅及び位相を用いて、所望の撮像を行う(ステップ206)。なお、ステップ202で計測したB1分布は、計測を行った被検体の部位に依存する。従って、被検体や撮像部位が変化した場合には、B1分布の再計測を行うことが好ましい。被検体や撮像部位の変更を含む手順を図14に示す。図14中、図2と同じステップは同じ符号で示している。このフローでは、撮像部位や被検体が変わったタイミングで(ステップ207)、ステップ201に戻り、B1分布の計測と、その結果を反映した各小型RFコイルの振幅及び位相の設定を行う。このように被検体や撮像部位が変わったときのみにB1分布計測を行うことにより、B1分布計測の回数を少なくすることができ、検査のスループットを向上させることができる。
【００９１】
なお、以上の実施形態では、RFコイルがマルチプルアレイコイルである場合を説明したが、RFコイルとして独立した複数のRFコイルを用いる場合にも同様に適用することができる。
【００９２】
また図2および図14に示すフローでは、複数のチャンネルからなる照射コイルについて、チャンネル毎にB1分布を計測する場合を示したが、2以上のチャンネルを組み合わせてB1分布を取得することも可能である。例えば、4チャンネルの照射コイルの場合、3チャンネルの組合せを異ならせて、4回の画像取得を行う。それにより、3チャンネルのB1分布の合成したB1分布が式(5)により算出される。例えば、4回のB1分布計測で得られるB1分布(B11〜B14)は次式(23)〜(26)のようになる。
【００９３】
B11＝B1a＋B1b＋B1c    (23)
B12＝B1b＋B1c＋B1d    (24)
B13＝B1a＋B1c＋B1d    (25)
B14＝B1a＋B1b＋B1d    (26)
よってこれら式から、次式(27)〜(30)により各チャンネルのB1a〜B1dを求めることができる。
【００９４】
B1a＝｛(B11＋B13＋B14)−2B12｝/3    (27)
B1b＝｛(B11＋B12＋B14)−2B13｝/3    (28)
B1c＝｛(B11＋B12＋B13)−2B14｝/3    (29)
B1d＝｛(B12＋B13＋B14)−2B11｝/3    (30)
各チャンネルのB1分布を求めた後のステップ204は、上述したとおりである。
【００９５】
また2チャンネルだけを用いてRFシミングを行う場合には、用いる2チャンネルを組み合わせて画像を取得してもよい。
【００９６】
また以上の実施形態では、B1分布計測の結果を用いて、RFコイルの位相と振幅を調整する場合を説明したが、本発明の特徴の一つは、B1分布計測の手法すなわち特定のB1分布計測シーケンスを備えたMRI装置であることに特徴があり、計測したB1分布を用いてRFコイルの位相と振幅を調整する代わりに、B1分布を用いて、撮像シーケンスにより収集したデータを補正することも可能である。その場合、単独のRFコイルについても本発明を適用することができる。
【００９７】
さらに計算したB1分布は、RFシミングに用いるのではなく、T2分布を取得するときの後処理として、B1不均一に起因するT2値の誤差を補正するのに用いることも可能である。B1分布を用いたT2値の補正については、例えば、非特許文献5に記載されている。
【００９８】
以上の各実施形態の説明で明らかになった本発明の特徴を纏めると以下のようになる。即ち、
本発明のMRI装置は、被検体に核磁気共鳴を起こさせるための高周波磁場(B1)を照射するRF照射部と、高周波磁場プリパルスの印加から経過時間(TI)を設けて核磁気共鳴信号を取得する信号取得シーケンスを有して成るB1分布計測シーケンスを用いて前記被検体を撮像する撮像部と、核磁気共鳴信号を用いて被検体の画像を再構成する演算部と、を備え、演算部は、画像を用いてRF照射部の照射磁場分布を求めるものであって、演算部は、経過時間の異なる複数の画像を用いて前記照射磁場分布を求めることを特徴とする。
【００９９】
好ましくは、B1分布計測シーケンスは、経過時間が異なる複数の信号取得シーケンスを含み、演算部は、経過時間が異なる複数の信号取得シーケンスでそれぞれ取得された核磁気共鳴信号を用いて、経過時間の異なる複数の画像を再構成する。
【０１００】
また好ましくは、撮像部は、高周波磁場プリパルスの印加条件が異なる複数のB1分布計測シーケンスを実行してそれぞれ画像を取得し、演算部は、高周波磁場プリパルスの印加条件が異なる複数のB1分布計測シーケンスでそれぞれ得られた画像を用いて前記照射磁場分布を求める。
【０１０１】
また好ましくは、複数のB1分布計測シーケンスの内で、第一のB1分布計測シーケンスは、高周波磁場プリパルスを用いずに実行される信号取得シーケンスであり、第二のB1分布計測シーケンスは、高周波磁場プリパルスの印加とそれに続く少なくとも2回の信号取得シーケンスを含む。
【０１０２】
また好ましくは、演算部は、第一のB1分布計測シーケンスの信号取得シーケンスで得られた核磁気共鳴信号を再構成した画像の信号強度と、第二のB1分布計測シーケンスの複数の信号取得シーケンスでそれぞれ得られた核磁気共鳴信号を再構成した複数の画像の信号強度との比を用いて、行列式を解くことにより、照射磁場分布を算出する。
【０１０３】
また好ましくは、撮像部は、複数の断面について、それぞれ前記B1分布計測シーケンスを実行し、演算部は、複数の断面のそれぞれについて、照射磁場分布を求める。
【０１０４】
また好ましくは、信号取得シーケンスは、マルチスライス撮像シーケンス又は三次元撮像シーケンスである。
【０１０５】
また好ましくは、撮像部は、高周波磁場プリパルスの印加条件が同一の複数の前記B1分布計測シーケンスを含み、該複数のB1分布計測シーケンスで画像の再構成に必要な全データを分割して取得する。
【０１０６】
また好ましくは、撮像部は、経過時間を異ならせて、B1分布計測シーケンスを複数回実行し、演算部は、各B1分布計測シーケンスで得られた核磁気共鳴信号を用いて、経過時間の異なる複数の画像を再構成する。
【０１０７】
また好ましくは、高周波磁場プリパルスは、フリップ角が90度以上の高周波磁場パルスである。
【０１０８】
また好ましくは、信号取得シーケンスは、フリップ角が10度以下の高周波磁場パルスを繰り返し印加して核磁気共鳴信号を取得する。
【０１０９】
また、好ましくは、撮像部は、複数の信号取得シーケンスは連続しており、高周波磁場パルスは一定の繰り返し時間で印加される。
【０１１０】
また好ましくは、信号取得シーケンスは、核磁気共鳴信号の取得を伴わない高周波磁場パルスの印加を含む。
【０１１１】
また好ましくは、演算部は、経過時間の異なる複数の画像の信号強度を用いて、連立方程式を解くことにより、照射磁場分布を求める。
【０１１２】
また好ましくは、信号取得シーケンスは、k空間を放射状に走査して核磁気共鳴信号を取得するパルスシーケンスである。
【０１１３】
また好ましくは、演算部は、算出した照射磁場分布に基づき、RF照射部に供給される高周波磁場の位相と振幅の少なくとも一方を調整する。
【０１１４】
また好ましくは、RF照射部は、複数のチャンネルを有してなり、撮像部は、チャンネル毎にB1分布計測シーケンスを実行して前記被検体を撮像し、演算部は、チェンネル毎に前記照射磁場分布を求める。
【０１１５】
また、本発明の照射磁場計測方法は、RF照射部を備えた磁気共鳴撮像装置における照射部の前記照射磁場分布を計測する方法であって、RF照射部による高周波磁場プリパルスの印加から経過時間(TI)を設けて核磁気共鳴信号を取得する信号取得シーケンスを実行する計測ステップと、核磁気共鳴信号を用いて被検体の画像を再構成する画像再構成ステップと、画像を用いてRF照射部の照射磁場分布を求める照射磁場分布算出ステップと、を備え、照射磁場分布算出ステップは、経過時間の異なる複数の画像を用いて照射磁場分布を求めることを特徴とする。
【０１１６】
好ましくは、計測ステップは、高周波磁場プリパルスの印加から経過時間が異なる複数の信号取得シーケンスを実行し、画像再構成ステップは、経過時間が異なる複数の信号取得シーケンスでそれぞれ取得された核磁気共鳴信号を用いて、経過時間の異なる複数の画像を再構成する。
【０１１７】
また好ましくは、計測ステップは、高周波磁場プリパルスを用いずに信号取得シーケンスのみを実行するステップと、高周波磁場プリパルスの印加とそれに続く少なくとも2回の信号取得シーケンスを実行するステップと、を含み、射磁場分布算出ステップは、各信号取得シーケンスでそれぞれ得られた核磁気共鳴信号を再構成して得られた複数の画像の信号強度の比を用いて、行列式を解くことにより、照射磁場分布を算出する。
【０１１８】
＜実施例1＞
図13(a)に示すB1分布計測シーケンスを用いて、以下のパラメータでTIの異なる3つの画像データを得た。撮像時間は15秒であった。
【０１１９】
プリパルス：フリップ角90度
TR：5000ms
位相エンコード数：64
TI：150ms、300ms、450ms
これら画像データを用いて、B1分布を算出した結果を図15に示す。
【０１２０】
＜実施例2＞
図13(b)に示すB1分布計測シーケンスを用いて、実施例1と同じパラメータでTIの異なる3つの画像データを得た。撮像時間は15秒であった。これら画像データを用いて、B1分布を算出した結果を図16に示す。
【０１２１】
これら実施例で得られた結果(図15、図16)は、従来法(Double  Angle法)で得られる結果(図17)と同等の精度であった。
【０１２２】
＜実施例3＞
図9に示すB1分布計測シーケンスを用いて、TR及び位相エンコード数は実施例1と同様にし、TI1＝50ms、TI2＝150msとして、一つの基準画像を含む3つの画像データを得た。撮像時間は、約500msであった。エンコード数を64から32に減らした場合には、さらに撮像時間を350msに短縮することができる。これら画像データを用いて、第三実施形態と同様の計算手法でB1分布を算出した結果を図18に示す。
【０１２３】
本発明の手法と従来法(DAM)との計算精度を比較した結果を図19に示す。図中、(a)はDAM、(b)は実施例1、(ｃ)は実施例3をそれぞれ示している。図19の横軸は、送信コイルに設定された実際のフリップ角を示し、縦軸はそれぞれの手法で算出しB1からも求めたフリップ角を示しており、計算の誤差が全くなければ縦軸と横軸は一致し、ｙ＝ｘのグラフとなる。実施例3の手法は最も精度よくB1分布を算出できることがわかる。
【０１２４】
また撮像時間については、同じ撮像パラメータで計測した場合、DAMでは10.8分かかるところ、実施例1、2は15秒と大幅に短縮でき、実施例3では500msとさらに短縮できた。
【産業上の利用可能性】
【０１２５】
本発明によれば、従来のRFパルスの磁場分布計測(B1計測)を極めて短時間で行うことができる。従って、撮像部位の変化に応じてリアルタイムでB1計測とその結果に基づくRFパルスの制御を行うことができるので、撮像時間延長による被検体への負担を軽減し、且つ被検体の体内の磁場の影響を受けやすい高磁場MRIにおいて、その影響を排除した診断能の高い画像を提供することができる。
【符号の説明】
【０１２６】
2  静磁場発生系、3  傾斜磁場発生系、4  シーケンサ、5  送信系、6  受信系、7  信号処理系、8  CPU、11  高周波発振器、12  変調器、13  増幅器、14a  高周波コイル(送信コイル)。
【図１】
2012060192.tif 000019
【図２】
2012060192.tif 000020
【図３】
2012060192.tif 000021
【図４】
2012060192.tif 000022
【図５】
2012060192.tif 000023
【図６】
2012060192.tif 000024
【図７】
2012060192.tif 000025
【図８】
2012060192.tif 000026
【図９】
2012060192.tif 000027
【図１０】
2012060192.tif 000028
【図１１】
2012060192.tif 000029
【図１２】
2012060192.tif 000030
【図１３】
2012060192.tif 000031
【図１４】
2012060192.tif 000032
【図１５】
2012060192.tif 000033
【図１６】
2012060192.tif 000034
【図１７】
2012060192.tif 000035
【図１８】
2012060192.tif 000036
【図１９】
2012060192.tif 000037
【手続補正書】
【提出日】平成25年8月5日(2013.8.5)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
被検体に核磁気共鳴を起こさせるための高周波磁場(B1)を照射するRF照射部と、
高周波磁場プリパルスの印加から経過時間(TI)を設けて核磁気共鳴信号を取得する信号取得シーケンスを有して成るB1分布計測シーケンスを用いて前記被検体を撮像する撮像部と、
前記核磁気共鳴信号を用いて前記被検体の画像を再構成する演算部と、
を備え、
前記演算部は、前記画像を用いて前記RF照射部の照射磁場分布を求める磁気共鳴撮像装置であって、
前記演算部は、前記経過時間の異なる複数の画像を用いて前記照射磁場分布を求めることを特徴とする磁気共鳴撮像装置。
【請求項２】
請求項1記載の磁気共鳴撮像装置において、
前記B1分布計測シーケンスは、前記経過時間が異なる複数の信号取得シーケンスを含み、
前記演算部は、前記経過時間が異なる複数の信号取得シーケンスでそれぞれ取得された核磁気共鳴信号を用いて、前記経過時間の異なる複数の画像を再構成することを特徴とする磁気共鳴撮像装置。
【請求項３】
請求項1に記載の磁気共鳴撮像装置において、
前記撮像部は、前記高周波磁場プリパルスの印加条件が異なる複数のB1分布計測シーケンスを実行してそれぞれ画像を取得し、
前記演算部は、前記高周波磁場プリパルスの印加条件が異なる複数のB1分布計測シーケンスでそれぞれ得られた画像を用いて前記照射磁場分布を求めることを特徴とする磁気共鳴撮像装置。
【請求項４】
請求項3に記載の磁気共鳴撮像装置において、
前記複数のB1分布計測シーケンスの内で、第一のB1分布計測シーケンスは、前記高周波磁場プリパルスを用いずに実行される信号取得シーケンスであり、第二のB1分布計測シーケンスは、前記高周波磁場プリパルスの印加とそれに続く少なくとも2回の信号取得シーケンスを含むことを特徴とする磁気共鳴撮像装置。
【請求項５】
請求項4に記載の磁気共鳴撮像装置において、
前記演算部は、前記第一のB1分布計測シーケンスの信号取得シーケンスで得られた核磁気共鳴信号を再構成した画像の信号強度と、前記第二のB1分布計測シーケンスの複数の信号取得シーケンスでそれぞれ得られた核磁気共鳴信号を再構成した複数の画像の信号強度との比を用いて、行列式を解くことにより、前記照射磁場分布を算出することを特徴とする磁気共鳴撮像装置。
【請求項６】
請求項1に記載の磁気共鳴撮像装置において、
前記撮像部は、複数の断面について、それぞれ前記B1分布計測シーケンスを実行し、
前記演算部は、前記複数の断面のそれぞれについて、前記照射磁場分布を求めることを特徴とする磁気共鳴撮像装置。
【請求項７】
請求項6に記載の磁気共鳴撮像装置において、
前記信号取得シーケンスは、マルチスライス撮像シーケンス又は三次元撮像シーケンスであることを特徴とする磁気共鳴撮像装置。
【請求項８】
請求項1に記載の磁気共鳴撮像装置において、
前記撮像部は、前記高周波磁場プリパルスの印加条件が同一の複数の前記B1分布計測シーケンスを含み、該複数のB1分布計測シーケンスで前記画像の再構成に必要な全データを分割して取得することを特徴とする磁気共鳴撮像装置。
【請求項９】
請求項1に記載の磁気共鳴撮像装置において、
前記撮像部は、前記経過時間を異ならせて、前記B1分布計測シーケンスを複数回実行し、
前記演算部は、各B1分布計測シーケンスで得られた核磁気共鳴信号を用いて、前記経過時間の異なる複数の画像を再構成することを特徴とする磁気共鳴撮像装置。
【請求項１０】
請求項1に記載の磁気共鳴撮像装置において、
前記高周波磁場プリパルスは、フリップ角が90度以上の高周波磁場パルスであることを特徴とする磁気共鳴撮像装置。
【請求項１１】
請求項2に記載の磁気共鳴撮像装置において、
前記信号取得シーケンスは、フリップ角が10度以下の高周波磁場パルスを繰り返し印加して前記核磁気共鳴信号を取得することを特徴とする磁気共鳴撮像装置。
【請求項１２】
請求項11に記載の磁気共鳴撮像装置において、
前記撮像部は、前記複数の信号取得シーケンスは連続しており、前記高周波磁場パルスは一定の繰り返し時間で印加されることを特徴とする磁気共鳴撮像装置。
【請求項１３】
請求項12に記載の磁気共鳴撮像装置において、
前記信号取得シーケンスは、核磁気共鳴信号の取得を伴わない高周波磁場パルスの印加を含むことを特徴とする磁気共鳴撮像装置。
【請求項１４】
請求項1に記載の磁気共鳴撮像装置において、
前記演算部は、前記経過時間の異なる複数の画像の信号強度を用いて、連立方程式を解くことにより、前記照射磁場分布を求めることを特徴とする磁気共鳴撮像装置。
【請求項１５】
請求項1に記載の磁気共鳴撮像装置において、
前記信号取得シーケンスは、k空間を放射状に走査して前記核磁気共鳴信号を取得するパルスシーケンスであることを特徴とする磁気共鳴撮像装置。
【請求項１６】
請求項1に記載の磁気共鳴撮像装置において、
前記演算部は、算出した前記照射磁場分布に基づき、前記RF照射部に供給される高周波磁場の位相と振幅の少なくとも一方を調整することを特徴とする磁気共鳴撮像装置。
【請求項１７】
請求項1に記載の磁気共鳴撮像装置において、
前記RF照射部は、複数のチャンネルを有してなり、
前記撮像部は、前記チャンネル毎に前記B1分布計測シーケンスを実行して前記被検体を撮像し、
前記演算部は、前記チェンネル毎に前記照射磁場分布を求めることを特徴とする磁気共鳴撮像装置。
【請求項１８】
RF照射部を備えた磁気共鳴撮像装置における前記RF照射部の照射磁場分布を計測する方法であって、
前記RF照射部による高周波磁場プリパルスの印加から経過時間(TI)を設けて核磁気共鳴信号を取得する信号取得シーケンスを実行する計測ステップと、
前記核磁気共鳴信号を用いて前記被検体の画像を再構成する画像再構成ステップと、
前記画像を用いて前記RF照射部の照射磁場分布を求める照射磁場分布算出ステップと、
を備え、
前記照射磁場分布算出ステップは、前記経過時間の異なる複数の画像を用いて前記照射磁場分布を求めることを特徴とする照射磁場計測方法。
【請求項１９】
請求項18に記載の照射磁場計測方法であって、
前記計測ステップは、前記高周波磁場プリパルスの印加から前記経過時間が異なる複数の信号取得シーケンスを実行し、
前記画像再構成ステップは、前記経過時間が異なる複数の信号取得シーケンスでそれぞれ取得された核磁気共鳴信号を用いて、前記経過時間の異なる複数の画像を再構成することを特徴とする照射磁場計測方法。
【請求項２０】
請求項19に記載の照射磁場計測方法であって、
前記計測ステップは、前記高周波磁場プリパルスを用いずに信号取得シーケンスのみを実行するステップと、前記高周波磁場プリパルスの印加とそれに続く少なくとも2回の信号取得シーケンスを実行するステップと、を含み、
前記照射磁場分布算出ステップは、各信号取得シーケンスでそれぞれ得られた核磁気共鳴信号を再構成して得られた複数の画像の信号強度の比を用いて、行列式を解くことにより、前記照射磁場分布を算出することを特徴とする照射磁場計測方法。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】０００３
【補正方法】変更
【補正の内容】
【０００３】
近年、超伝導磁石の開発に伴い、3T以上の高磁場を実現できる高磁場MRI装置が普及している。高磁場MRI装置では、高いSN画像が得られるが、腹部撮像等で画像に輝度ムラが発生する問題がある。この輝度ムラの原因の一つとして、被検体の組織内の原子核スピンを励起する高周波磁場パルス(以下、RFパルスという)の磁場分布(照射磁場分布、B1分布)の空間的な不均一がある。一般に、励起のための高周波磁場の共鳴周波数は、静磁場強度に比例するため、高磁場MRI装置では、従前の高周波磁場よりも高い周波数の磁場を照射する必要がある。その場合、生体内での高周波磁場の波長が生体(特に腹部)の大きさと同じスケールとなる。このため、高周波磁場の位相が生体内の位置によって変化し、画像ムラとなって現れる。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】０００５
【補正方法】変更
【補正の内容】
【０００５】
B1分布を計測(B1計測)する一般的な方法は、Double Angle法(DAM)と呼ばれる手法であり、任意のフリップ角のRFパルスを用いて撮像した画像と、その2倍のフリップ角のRFパルスとを用いて撮像した画像との演算によりB1分布を計測する。また複数の異なるフリップ角のRFパルスを用いて複数の画像を取得し、パルスシーケンスにより決定される信号強度式により、複数の画像の信号強度(画素値)をフィッティングすることでB1分布を計算する手法も提案されている(非特許文献1)。また同様に取得した複数の画像に対し、フィッティングを用いないで、信号強度変化の周期からB1分布を計算する手法も提案されている(特許文献1)。
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００１４
【補正方法】変更
【補正の内容】
【００１４】
本発明は、B1分布計測に、プリパルスの印加と信号取得パルスシーケンスとを組み合わせるとともに、プリパルス印加からの経過時間(TI)が異なる複数の信号取得パルスシーケンスを実行し、これら複数の信号取得パルスシーケンスの実行によって取得したTIの異なる複数の画像(画像用データ)の演算によって、B1分布を求めるという新規な手法を提案するものである。
【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】００４３
【補正方法】変更
【補正の内容】
【００４３】
最初の信号取得シーケンスで得られたk空間データを逆フーリエ変換することにより得た画像データの、ある注目画素の信号強度S(B1,TI)は式(1)で与えられる。
2012060192.tif 000038
式(1)中、Sseqは、プリパルスの後の信号取得シーケンスによって決まる信号強度を表し、αは設定したプリパルスのフリップ角を表し、TIはプリパルス印加からk空間中心の信号を収集するまでの時間を表し、T1は組織に依存する縦緩和時間を表す。
【手続補正６】
【補正対象書類名】明細書
【補正対象項目名】００５８
【補正方法】変更
【補正の内容】
【００５８】
＜＜B1分布算出(ステップ203)＞＞
次に、図7に示すB1計測パルスシーケンスによって得られたk空間データから、B1分布を算出する手法を説明する。
【手続補正７】
【補正対象書類名】明細書
【補正対象項目名】００６８
【補正方法】変更
【補正の内容】
【００６８】
一方、プリパルスの直前の基準画像取得シーケンスによって得られた画像の同じ注目画素の信号強度は、式(16)においてα＝0とした場合と同じであるため、式(10)と同一の式(17)で与えられる。
2012060192.tif 000039
式(16)を式(17)で除算し、自然対数を取ると、式(18)のように、log(1-cos(B1・α))と(−TIk/T1)の線形結合で表わすことができる。
2012060192.tif 000040
各信号取得シーケンスから得られたTIの異なる画像について同様の計算をすると、式(19)の連立方程式が得られる。
2012060192.tif 000041
ここで、Sは1×nの行列、Aは2×nの行列、Xは1×2の行列である。Wi(i＝1，2，3・・・n)はそれぞれのTIに対する重みを示し、任意に設定することができる。行列Aの擬似逆行列pinvAを左から掛けることで、式(19)を解くことができ、次式(20)、(21)のようにB1及びT1を求めることができる。
【手続補正８】
【補正対象書類名】明細書
【補正対象項目名】００６９
【補正方法】変更
【補正の内容】
【００６９】
2012060192.tif 000042
第一及び第二実施形態は、式(4)〜(6)或いは式(13)及び(14)の連立方程式を解くことにより、解(B1)を求めるので、解が発散しないためには、プリパルス後の信号取得シーケンスのTIを整数比としていたが、本実施形態では、そのようなTIの制限が要せずに解を求めることができるので、プリパルスの効果が高い、短TIで信号取得でき、計測時間を短縮できるとともに高精度なB1分布を得ることができる。また本実施形態では、B1分布と併せてT1分布を算出することが可能である。算出されたT1分布は、例えば、非特許文献4に記載されているように、骨関節炎の診断等に利用される。
【手続補正９】
【補正対象書類名】明細書
【補正対象項目名】００７０
【補正方法】変更
【補正の内容】
【００７０】
なお図9では、プリパルス901の直前に、基準画像取得シーケンス900を実行し、基準画像を取得する場合を説明したが、プリパルス901印加後に十分の時間即ちT1緩和するに十分な時間経過後に、基準信号取得シーケンス900を実行してもよい。図10に、この変更例を示す。図示するように、プリパルス901印加後に異なるTIで信号取得シーケンス903、905、907・・・を実行してTIの異なる画像を取得し、最後にTIを十分長く設定して基準画像取得シーケンス900を実行して基準画像を取得する。この基準画像の信号強度S0は、式(17)で示したように、S0＝Sseqと表すことができるので、図9の場合と同様に、式(18)〜式(21)により、B1分布及びT1分布を算出することができる。
【手続補正１０】
【補正対象書類名】明細書
【補正対象項目名】００７３
【補正方法】変更
【補正の内容】
【００７３】
即ち、例えば図4に示すシーケンスにおいて、RFパルス401と同時に印加されるスライス選択傾斜磁場402は、信号取得シーケンス毎に異なる軸のスライス傾斜磁場であり、異なる断面が励起される。これにより、各断面の画像は、プリパルス911印加からの経過時間TIが異なることになる。所定の間隔(TD)を置いて、同様のB1分布計測シーケンス920、930を繰り返す。その際、各B1分布計測シーケンスに含まれる信号取得シーケンスにおいて選択されるスライス(断面)の順序をサイクリックに変更する。例えば、最初のB1分布計測シーケンス910では、アキシャル面(AX)、サジタル面(SG)、コロナル面(COR)とし、次のB1分布計測シーケンス920では、SG→COR→AXの順、その次のB1分布計測シーケンス930では、COR→AX→SGの順とする。これによって、3回目のB1分布計測シーケンス930が終了した時点では、アキシャル面(AX)、サジタル面(SG)、コロナル面(COR)のそれぞれについて、TIの異なる3つの画像データ(TI1、TI2、TI3のk空間データ)を得ることができる。
【手続補正１１】
【補正対象書類名】明細書
【補正対象項目名】００８０
【補正方法】変更
【補正の内容】
【００８０】
本実施形態では、図12に示すように、プリパルス101と複数の信号取得シーケンス103、105、107とからなるB1分布計測シーケンス100を複数回行う。各回で同一フリップ角のプリパルス101を用いるが、1回の信号取得シーケンスではk空間データの一部を取得し、複数回で全k空間データを取得するようにする。k空間の分割の仕方は、特に限定されず、図示するように単純に複数の領域に分割してもよいし、それぞれの信号取得シーケンス100、120、130・・・で位相エンコードステップを粗くして計測し、k空間の中心に最も近いデータを常に所定のTIで計測するようにしてもよい。
【手続補正１２】
【補正対象書類名】図面
【補正対象項目名】図２
【補正方法】変更
【補正の内容】
【図２】
2012060192.tif 000043
【手続補正１３】
【補正対象書類名】図面
【補正対象項目名】図４
【補正方法】変更
【補正の内容】
【図４】
2012060192.tif 000044
【手続補正１４】
【補正対象書類名】図面
【補正対象項目名】図１１
【補正方法】変更
【補正の内容】
【図１１】
2012060192.tif 000045
【手続補正１５】
【補正対象書類名】図面
【補正対象項目名】図１２
【補正方法】変更
【補正の内容】
【図１２】
2012060192.tif 000046
【手続補正１６】
【補正対象書類名】図面
【補正対象項目名】図１４
【補正方法】変更
【補正の内容】
【図１４】
2012060192.tif 000047
【国際調査報告】
2012060192.tif 000048
2012060192.tif 000049
2012060192.tif 000050
2012060192.tif 000051