(57)【要約】【課題】体動による空間位相分布の位置ずれの影響を回避し、安定した温度変化の情報を提供すること。【解決手段】温度計測用パルスシーケンスを実行して検査対象内の温度情報を含むエコー信号を計測し、時系列データを得て、その空間位相分布から温度変化分布を求める際に、呼吸動に同期してエコー信号に付与する位相ンコードを制御し、同じ変位範囲にあるタイミングで計測したエコー信号でｋ空間の低周波領域を充填する。このような時系列データを用いて温度分布画像を求めることにより、信頼性の高い温度分布画像の情報を提供できる。【選択図】図４
2004049478.tif 000002
【特許請求の範囲】
【請求項１】
検査対象に静磁場を与える静磁場発生手段と、前記検査対象に傾斜磁揚を印加する傾斜磁場発生手段と、前記検査対象を構成する原子の原子核に磁気共鳴を起こさせる高周波パルスを発生する高周波パルス発生手段と、前記検査対象が発生する核磁気共鳴信号をｋ空間データとして検出する検出手段と、時系列的に計測されたｋ空間データを用いて前記検査対象の温度変化分布を演算する演算手段と、前記各手段の制御を行う制御手段と、前記検査対象の周期動を計測する手段を備え、前記制御手段は、前記周期動に応じて前記傾斜磁場発生手段を制御し、前記ｋ空間の少なくとも低周波領域データを、前記周期動における同じ変位範囲で計測した核磁気共鳴信号で充填することを特徴とする磁気共鳴イメージング装置。
【請求項２】
前記制御手段は、前記周期動の少なくとも２つの異なる変位範囲に対応して２以上のｋ空間データを収集することを特徴とする請求項１記載の磁気共鳴イメージング装置。
【請求項３】
前記検査対象の周期動を計測する手段は、核磁気共鳴信号を用いて前記検査対象の周期動を算出することを特徴とする請求項１記載の磁気共鳴イメージング装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、磁気共鳴イメージング（以下、ＭＲＩと称す）装置に関し、特にＭＲＩ装置による生体内温度分布画像の計測に好適な技術に関する。
【０００２】
【従来の技術】
近年、ＭＲＩ装置を術中モニタとして使用するインターベンショナルＭＲＩ（Ｉｎｔｅｒｖｅｎｔｉｏｎａｌ　ＭＲＩ：以下、ＩＶ−ＭＲＩと称す）が注目されている。ＩＶ−ＭＲＩで行われる治療法には、レーザ治療、マイクロ波凝固術、エタノールなどの薬物注入、ＲＦ照射切除、低温治療などがある。これらの治療において、ＭＲＩは、患部に穿刺針や細管を到達させるためのリアルタイムイメージングによるガイド及び治療中の組織変化の可視化、加熱・冷却治療中の局所温度のモニタなどの役割を果たす。ＩＶ−ＭＲＩの典型的な応用例としては、マイクロ波凝固術中における体内の温度分布の画像化が挙げられる。
【０００３】
温度分布の画像化手法には、信号強度から求める方法、拡散係数から求める方法、プロトンの位相シフトから求める方法（ＰＰＳ法：Ｐｒｏｔｏｎ　Ｐｈａｓｅ　Ｓｈｉｆｔ法）等があるが、この中でＰＰＳ法が最も測定精度に優れている。この機能を活用して、生体内の温度をモニタし、レーザ照射治療のモニタや、マイクロ波凝固術中の温度モニタができるようになりつつある。
【０００４】
ＰＰＳ法は、エコー信号の位相成分に含まれる（共鳴周波数）×（静磁場）成分の温度依存性を利用するものであり、グラディエントエコー系のパルスシーケンスで計測したエコー信号の位相情報から温度分布を求める。図１は、このような位相情報を計測する方法の一例である。以下、図１に示す計測方法と位相分布・温度分布の計算方法を説明する。
【０００５】
まず、スライス選択のための傾斜磁場ＧｓｌＯ２の印加と共に、９０°高周波パルス１０１（以下、９０°パルスと称す）で試料のスピンを励起し、スピンの位相を変化させるために傾斜磁場Ｇｐ１０３、読み出し傾斜磁場Ｇｒ１０４を順次印加し、グラディエントエコー１０６を発生させる。位相エンコードの異なる複数のエコー信号をフーリエ変換して得られる複素画像の実部と虚部から位相分布を、例えば、式（１）より求める。
【０００６】
【数１】
2004049478.tif 000003
【０００７】
そして、得られた位相分布、エコー信号が最大となる時点と９０°パルスとの間隔ＴＥ（１０６）、共鳴周波数ｆ、水の温度系数から、例えば、式（２）より温度Ｔを求める。
【０００８】
【数２】
2004049478.tif 000004
【０００９】
上記手法を用いて、異なる時刻ｔｌ、ｔｎ（ｎ：撮影回数）で取得した信号からそれぞれ計算した温度分布の差分をとり、ある時間における被検体の温度変化の分布を取得することができる。
【００１０】
【発明が解決しようとしている課題】
前述のようにＭＲＩによる温度モニタリングを行う場合、連続した時系列データを取得し、異なる時刻で得られた空間位相分布の差分を温度変化に換算するため、常に同一の温度変化領域を撮像する必要がある。
【００１１】
しかし、生体への適用においては呼吸動等の体動を抑制することは困難であり、温度変化領域の空間位置が移動することが多々あり、安定して同一の温度変化領域を計測することは困難である。このような温度変化領域の移動は、温度モニタリングの信頼性を大幅に損なうことになる。例えば、加熱治療を例にとると、呼吸動の抑制が不十分な状態で計測した時系列データには、加熱部位の温度上昇の情報を含むデータと、含まないデータが混在することになり、後者の場合は加熱による温度上昇の情報が得られないことになる。このような温度画像を加熱中にリアルタイムで計測・表示した場合は、温度が上昇したり、しなかったり、場合によっては突然加熱領域が広がったり、消えたりして、安定した温度モニタリングを行うことができない。
【００１２】
一般に形態画像の取得を目的としたＭＲＩにおいて、呼吸動によるアーチファクトを抑制する手法は種々提案されている。例えば、Ｒｅｓｐｉｒａｔｏｒｙ　ｏｒｄｅｒｅｄ　ｐｈａｓｅ　ｅｎｃｏｄｉｎｇ（ＲＯＰＥ）：Ａ　ｍｅｔｈｏｄ　ｆｏｒ　ｒｅｄｕｃｉｎｇ　ｒｅｓｐｉｒａｔｏｒｙ　ｍｏｔｉｏｎ　ａｒｔｉｆａｃｔｓ　ｉｎ　ＭＲ　ｉｍａｇｉｎｇ；Ｊｏｕｒｎａｌ　ｏｆ　Ｃｏｍｐｕｔｅｒ　　Ａｓｓｉｓｔｅｄ　Ｔｏｍｏｇｒａｐｈｙ　９（４）：８３５−８３８（１９８５）では、位相エンコードを制御することにより体動アーチファクトを抑制する方法が提案されている。
【００１３】
しかし温度計測においては、単に画像における体動アーチファクトを抑えるだけでなく、臓器の位置変動自体を抑える必要がある。
そこで本発明は、温度変化分布を計測する機能を備えたＭＲＩ装置において、体動による影響を回避して、温度モニタリングの正確性、信頼性を向上させ、温度変化分布画像を安定して提供することを目的とする。
【００１４】
【課題を解決するための手段】
上記目的を達成するために、本発明のＭＲＩ装置は、検査対象に静磁場を与える静磁場発生手段と、前記検査対象に傾斜磁揚を印加する傾斜磁場発生手段と、前記検査対象を構成する原子の原子核に磁気共鳴を起こさせる高周波パルスを発生する高周波パルス発生手段と、前記検査対象が発生する核磁気共鳴信号をｋ空間データとして検出する検出手段と、時系列的に計測されたｋ空間データを用いて前記検査対象の温度変化分布を演算する演算手段と、前記各手段の制御を行う制御手段と、前記検査対象の周期動を計測する手段を備え、前記制御手段は、前記周期動に応じて前記傾斜磁場発生手段を制御し、前記ｋ空間の少なくとも低周波領域データを、前記周期動における同じ変位範囲で計測した核磁気共鳴信号で充填することを特徴とする。
【００１５】
本発明のＭＲＩ装置において、制御手段は、周期動の少なくとも２つの異なる変位範囲に対応して２以上のｋ空間データを収集してもよい。
本発明のＭＲＩ装置によれば、温度変化分布を演算するためのｋ空間データのうち、少なくとも低周波領域のデータを、検査対象の変位がほぼ同一のときに計測した核磁気共鳴信号とするので、体動の影響を抑制した計測結果を得ることができ、これにより温度モニタリングの正確性、信頼性を向上させることができる。
【００１６】
【発明の実施の形態】
以下、本発明の実施形態を図面を参照して詳細に説明する。まず、本実施形態におけるＭＲＩ装置の構成を図２により説明する。
【００１７】
ＭＲＩ装置は、被検体２０１内部に一様な静磁場Ｈ０を発生させるための電磁石または永久磁石より構成された静磁場発生磁気回路２０２、被検体２０１に直交するｘ、ｙおよびｚの３方向に強度が線形に変化する傾斜磁場Ｇｘ，Ｇｙ，Ｇｚを発生する傾斜磁場コイル２０９、被検体２０１に高周波磁場を発生する送信コイル２１４ａ、被検体２０１から生じる核磁気共鳴信号を検出するための検出コイル２１４ｂ、傾斜磁場、高周波パルスを所定のタイミングで発生させるためのシーケンサ２０７、シーケンサの制御や画像処理等の各種処理を行うコンピュータ２０８、コンピュータ２０８が行なう信号処理や制御のためのプログラムや必要なデータ及び処理結果を格納するとともに処理結果などを画像として表示する信号処理系２０６、各種パラメータの設定等の操作を行う操作部２２１等を有して構成される。なお、本実施例では、送信コイル２１４ａと検出コイル　２１４ｂは別々に設けているが、送受信両用のものでもよい。
【００１８】
次に動作の概要を説明する。シンセサイザ２１１により発生させた高周波を変調器２１２で変調し電力増幅器２１３で増幅し、送信コイル２１４ａに供給することにより被検体２０１の内部に高周波磁場を発生させ、核スピンを励起させる。通常は１Ｈを対象とするが、３１Ｐ、１３Ｃ等、核スピンを有する他の原子核を対象とすることもある。この際、エコー信号に位置情報を付与するために、スライス方向、位相エンコード方向および読み出し方向の傾斜磁場を印加する。
【００１９】
被検体２０１から放出される核磁気共鳴信号（エコー信号）は検出コイル２１４ｂにより受信され、増幅器２１５を通った後、検波器２１６で直交位相検波され、Ａ／Ｄ変換器２１７を介してコンピュータ２０８へ入力される。サンプリング後のエコー信号は、位相エンコード（ｋｙ）を縦軸とするｋ空間に配置される。
【００２０】
コンピュータ２０８はｋ空間データにフーリエ変換等の演算を行なった後、核スピンの密度分布、緩和時間分布、スペクトル分布等に対応する画像を作成し、ＣＲＴ等のディスプレイ２２８に表示する。またコンピュータ２０８は、シーケンサ２０７に指令を送り、傾斜磁場発生系２０３、送信系２０４及び検出系２０５を制御する。この制御のシーケンスは、撮像方法によって種々のものが用意されており、予め信号処理系２０６の記憶装置（ＲＯＭ２２４、ＲＡＭ２２５）に格納されている。本発明のＭＲＩ装置では、温度計測のためのパルスシーケンス、例えば図１に示すパルスシーケンスを備えると共に、このパルスシーケンスの実行によって取得したエコー信号について温度を求める演算を行う。
【００２１】
さらにこのＭＲＩ装置は、被検体２０１の体動をモニタするための外部位置センサ２３０が備えられている。外部位置センサ２３０は、例えば呼吸動をモニタする場合には、圧力センサを採用することができる。圧力センサを被検体の腹部に固定し、圧力の変動を計測することにより、呼吸動により位置変位を計測することができる。尚、以下の説明では、呼吸動の検出として、腹壁の動きを呼吸動とみなすこととするが、横隔膜の動きをモニタしても同様に応用できる。
【００２２】
圧力センサ２３０から出力される腹壁の変位量に対応する信号は、コンピュータ２０８に入力される。コンピュータ２０８は、この変位量に基き、シーケンサ２０７に指令を送り、温度計測用パルスシーケンスにおける位相エンコードを制御する。位相エンコードの制御は、位相エンコード方向の傾斜磁場の強度及びその変更順序について行う。通常、１枚の断面の情報を得るために、６４、１２８、２５６の位相エンコードが付与される。
【００２３】
次に、温度計測及び撮影手順について説明する。図３は、その手順を示すフロー図である。まず、圧力センサにより腹壁の動きをモニタし、被検体の呼吸動を検出する（ステップ３０１）。図４（ａ）は、呼吸動（腹部の変位）をグラフとして表示したものである。このように呼吸動は周期的なパターンを示す。この変位パターンをもとに、ｋ空間の低周波領域、高周波領域を充填する仕方を決定する（ステップ３０２）。
【００２４】
よく知られているように、ｋ空間の中央（低周波領域）に配置されるデータが画質に最も影響を与える。そこで、このようなｋ空間の性質を利用して、同一変位であって動きの少ないときに、画像のコントラストを決定する領域のデータを計測するように位相エンコード量を決定する。
【００２５】
例えば、図４において、被検体の変位が範囲４１１にあるときに計測するエコー信号に、低周波域の位相エンコード４１０を割り当てる。変位が範囲４３１では、被検体の動きが大きいので、このとき計測するエコー信号には高周波領域の位相エンコード４３０ａ、４３０ｂを割り当てる。変位の範囲４２１では、被検体の動きは少ないが、低周波域４１０を計測したときと腹部の変位が異なり、目的とする臓器の位置が異なることが予想されるため、低周波領域４１０の計測には用いずに、高周波領域４３０ａ、４３０ｂとの間の領域４２０ａ、４２０ｂの計測を行う。
【００２６】
予め計測した呼吸動パターンにおいて最大近傍の範囲４１１よりも最小値近傍の範囲４２１の方が変動が少ない時間が長い場合には、範囲４２１にあるときにｋ空間の低周波領域４１０を充填し、範囲４１１でその外側の領域４２０を充填するものとする。この場合にも、ｋ空間の高周波領域４３０は、呼吸動の大きいタイミング４３１において計測する。このようなｋ空間を充填する領域の割合は、予め計測した呼吸パターンから計算する。
このように変位の範囲毎にｋ空間領域を割り当てたものは、例えばテーブルとして、コンピュータ２０８の記憶装置に格納される。
【００２７】
次に、図１のパルスシーケンスを実行し、目的部位の空間分布の計測を行う（ステップ３０３）。このパルスシーケンスの繰り返し毎に、コンピュータ２０８は圧力センサ２３０で計測した被検体（腹壁）の変位をもとに、位相エンコード傾斜磁場１０３を制御し、例えば腹部の変位が小さいタイミング３１１では低周波域の位相エンコード量となる傾斜磁場を印加してｋ空間の低周波領域４１０を計測する。
【００２８】
このようにｋ空間の低周波領域４１０を常に同じ変位のときに計測することにより、各データにおいて呼吸動による臓器の移動の影響は解消され、どの画像においても、臓器の位置は同様の位置となる。
こうして得られた各データの複素差分画像から式（１）により位相分布を求め、さらに式（２）を用いて温度変化分布を計算する（ステップ３０４）。
【００２９】
【数３】
2004049478.tif 000005
【００３０】
こうして得られる温度変化分布は、温度変化領域の空間位置が常に一定であり正確な温度モニタリングを可能にする。
なおｋ空間充填法は、変位が同じであって移動が少ないときに、ｋ空間の低周波領域を計測するものであれば、図４に示すものに限定これに限るものではなく、種々の変更が可能である。
【００３１】
図５〜図７に他の実施例を示す。図５に示す実施例は、被検体の腹部の位置に応じて、呼期と吸期においてそれぞれの温度画像を計測する方法である。この方法でも、温度情報を含む空間位相分布を計測する前に、被検体の呼吸動を予めモニタしておく。そして、その呼吸パターンから、腹部の変位の最大と最小を予測し、その中間値より大きい場合と小さい場合の２つのパターン（範囲）５１１、５２１に分け、それぞれに異なるｋ空間を割り当てる。即ち、腹部の変位が、最大値を含む範囲５１１となるタイミングではｋ空間１（５１０）を充填し、最小値を含む範囲５２１となるタイミングではｋ空間２（５２０）をそれぞれ充填する。
【００３２】
こうして計測された時系列データにおいても、図４の実施例と同様に、呼吸動による臓器の移動の影響は減少する。これら時系列データをそれぞれフーリエ変換して得られた複素差分画像から位相分布、温度変化分布を計算することにより、呼期と吸期のそれぞれについて、呼吸動によるずれが緩和された温度変化領域のモニタリングができる。なお、呼期と吸期の２パターンだけでなく、さらに細かい時相に分割し、それぞれ異なるｋ空間を充填すれば、より呼吸動の影響の少ない温度変化分布を得ることができる。
【００３３】
図６に示す実施例は、図４及び図５の実施例を組み合わせた手法である。本実施例でも、吸期６１１、呼期６２１の２つのパターンに分け、それぞれのタイミングで計測したデータで異なるｋ空間１、２を充填する。但し、本実施例では、吸期６１１及び呼期６２１を、更に腹部の変動の小さいタイミング６１１ａ、６２１ａと大きいタイミング６１１ｂ、６２１ｂに分ける。そして、吸期６１１で計測したデータのうち、変動の小さいタイミング計測したデータでｋ空間１の低周波領域を充填し、変動の大きいタイミング計測したデータで高周波領域を充填する。同様に、呼期６２１で計測したデータのうち、変動の小さいタイミング計測したデータでｋ空間２の低周波領域を充填し、変動の大きいタイミング計測したデータで高周波領域を充填する。これによって、より安定した温度変化分布を得ることができる。
【００３４】
さらに図７に示す実施例では、体動パターンを、吸期６１１及び呼期６２１の腹部変動の小さいタイミング７１１、７１３とその中間７１２とに分け、吸期７１１で計測したデータでｋ空間１を充填し、呼期７１３で計測したデータでｋ空間２を充填するとともに、中間７１２で計測したデータは、ｋ空間１の高周波領域（７１０ａ）とｋ空間２の高周波領域（７２０ａ）で共有する。これにより計測の時間短縮を図ることができる。このように本実施形態によれば、予め計測した体動パターンをもとにｋ空間充填法を設定し、その際、変位が同じときに計測したデータでｋ空間の低周波領域を充填するようにしたので、呼吸動によるずれが緩和された温度変化領域のモニタリングができる。
【００３５】
尚、上記実施形態では、呼吸動の検出に圧力センサなどの外部位置センサを用いた場合を説明したが、このような外部位置センサの代わりに、ナビゲーションエコーを用いた変位量のモニタも可能である。ナビゲーションエコーを用いた体動モニタリングについては、形態画像において体動アーチファクトを抑制する手法として、例えば、「Ｐｒｏｓｐｅｃｔｉｖｅ　ａｄａｐｔｉｖｅ　ｎａｖｉｇａｔｏｒ　ｃｏｒｒｅｃｔｉｏｎ　ｆｏｒ　ｂｒｅａｔｈ−ｈｏｌｄ　ＭＲ　ｃｏｒｏｎａｒｙ　ａｎｇｉｏｇｒａｐｈｙ；Ｍａｇｎｅｔｉｃ　Ｒｅｓｏｎａｎｃｅ　ｏｆ　Ｍｅｄｉｃｉｎｅ　３７：１４８−　１５２（１９９７）」に記載された手法などを採用することができる。
【００３６】
この場合には、例えば、図１のパルスシーケンスの繰り返し時間内で、位相エンコードを付与しないエコー信号をナビゲーションエコーとして計測し、このナビゲーションエコーを読み出し方向にフーリエ変換したデータから変位を求めることが可能である。
また上記実施形態では、温度計測のための撮影方法として、図１のパルスシーケンスを採用したが、パルスシーケンスとしては、グラディエントエコー系のシーケンスであれば、他のシーケンスであっても良い。シーケンスの具体例としては、ＳＡＲＧＥ、ＴＲＳＡＲＧＥ、ＲＦＳＡＲＧＥなどの高速ＧｒＥシーケンス、ＳＳＦＰ（Ｓｔｅａｄｙ　Ｓｔａｔｅ　Ｆｒｅｅ　Ｐｒｅｃｅｓｓｉｏｎ）のシーケンス、ＧｒＥ型のＥＰＩ（Ｅｃｈｏ　Ｐｌａｎａｒ　Ｉｍａｇｉｎｇ）シーケンスなど公知のパルスシーケンスが利用できる。
【００３７】
【発明の効果】
本発明によれば、被検体の呼吸動を検出し、呼吸に同期して温度変化情報を含む空間位相分布を連続して計測し、これを用いて温度分布画像を求めることにより、信頼性の高い温度分布画像の情報を提供できる。
【図面の簡単な説明】
【図１】従来のＰＰＳ法及び本発明に用いる撮影シーケンスを示す図。
【図２】本発明が適用される装置の概要を示すブロック図。
【図３】本発明による温度計測の手順を示すフロー図。
【図４】本発明による位相エンコード制御の一実施例を示す図。
【図５】本発明による位相エンコード制御の他の実施例を示す図。
【図６】本発明による位相エンコード制御の他の実施例を示す図。
【図７】本発明による位相エンコード制御の他の実施例を示す図。
【符号の説明】
２０１…被検体、２０２…静磁場発生磁気回路、２０３…傾斜磁場発生系、２０４…送信系、２０６…信号処理系、２０７…シーケンサ、２０８…コンピュータ、２３０…圧力センサ
【図１】
2004049478.tif 000006
【図２】
2004049478.tif 000007
【図３】
2004049478.tif 000008
【図４】
2004049478.tif 000009
【図５】
2004049478.tif 000010
【図６】
2004049478.tif 000011
【図７】
2004049478.tif 000012