(57)【特許請求の範囲】
【請求項１】
被検体の周期的な体動についての体動情報を検出する体動情報検出部と、
撮像シーケンスに基づいて、前記周期的な体動情報に同期させた同期撮像を実行して、3次元K空間データの計測を制御する計測制御部と、
前記3次元K空間データを用いて前記被検体の血管画像を再構成する演算処理部と、
を有する磁気共鳴イメージング装置であって、
前記撮像シーケンスは、3次元K空間(各軸をkr,k1,k2とする)におけるリードアウト方向krに垂直な2方向k1,k2に非平行な計測軌跡に沿ってエコーデータを計測するものであり、
前記計測制御部は、前記撮像シーケンスの繰り返し時間が前記体動情報の複数周期となるように、前記同期撮像を制御し、
前記撮像シーケンスは、前記3次元K空間のリードアウト方向に垂直な2方向に対応する実空間の方向に、それぞれ、前記被検体の血流の核磁化の位相を変調するためのリフェーズ傾斜磁場又はディフェーズ傾斜磁場を印加するものであることを特徴とする磁気共鳴イメージング装置。
【請求項２】
請求項1記載の磁気共鳴イメージング装置において、
前記計測軌跡は、一つの直線状軌跡を、前記3次元K空間のリードアウト方向を回転軸として回転させて得られる複数の直線状計測軌跡であり、
前記計測制御部は、前記撮像シーケンスを繰り返して、異なる直線状軌跡に沿うエコーデータをそれぞれ計測することを特徴とする磁気共鳴イメージング装置。
【請求項３】
請求項1又は2記載の磁気共鳴イメージング装置において、
前記周期的な体動情報は、心電図であり、
前記同期撮像は、心電図のR波に同期する撮像であり、
前記計測制御部は、前記R波から所定の遅延時間(DT)を空けて前記撮像シーケンスを開始することを特徴とする磁気共鳴イメージング装置。
【請求項４】
請求項3記載の磁気共鳴イメージング装置において、
前記遅延時間は、前記撮像シーケンスの実行が、拡張期又は収縮期となる時間であることを特徴とする磁気共鳴イメージング装置。
【請求項５】
請求項1又は2記載の磁気共鳴イメージング装置において、
前記演算処理部は、リファレンススキャンにより前記被検体を事前に撮像して取得されたデータに基づいて撮像パラメータ値を決定し、該決定した撮像パラメータ値に基づいて前記撮像シーケンスを設定し、
前記計測制御部は、前記設定された撮像シーケンスを実行することを特徴とする磁気共鳴イメージング装置。
【請求項６】
請求項5記載の磁気共鳴イメージング装置において、
前記撮像パラメータ値は、前記撮像シーケンスの実行が、前記被検体の心電図において収縮期又は拡張期となるような、心電図R波からの遅延時間(DT)を含むことを特徴とする磁気共鳴イメージング装置。
【請求項７】
請求項5又は6記載の磁気共鳴イメージング装置において、
前記計測制御部は、PC法パルスシーケンスを用いて前記リファレンススキャンを行い、
前記演算処理部は、前記リファレンススキャンで取得されたデータに基づいて、注目する血流部分の流速変化情報を取得することを特徴する磁気共鳴イメージング装置。
【請求項８】
請求項7記載の磁気共鳴イメージング装置において、
前記流速変化情報を流速変化グラフとして表示する表示部と、
前記表示部に表示された前記流速変化グラフ上で、遅延時間(DT)の設定入力を受け付ける入力部と、
を備えたことを特徴とする磁気共鳴イメージング装置。
【請求項９】
請求項1乃至8のいずれか一項に記載の磁気共鳴イメージング装置において、
前記計測軌跡は、複数の平行直線状計測軌跡を含んで成る単位軌跡群を、前記3次元K空間のリードアウト方向を回転軸として回転させて得られる複数の単位軌跡群であることを特徴とする磁気共鳴イメージング装置。
【請求項１０】
被検体の周期的な体動についての体動情報を検出する体動情報検出部と、
撮像シーケンスに基づいて、前記周期的な体動情報に同期させた同期撮像を実行して、
3次元K空間データの計測を制御する計測制御部と、
前記3次元K空間データを用いて前記被検体の血管画像を再構成する演算処理部と、
を有する磁気共鳴イメージング装置であって、
前記撮像シーケンスは、3次元K空間(各軸をkr,k1,k2とする)におけるリードアウト方向krに垂直な2方向k1,k2に非平行な計測軌跡に沿ってエコーデータを計測するものであり、
前記計測制御部は、前記撮像シーケンスの繰り返し時間が前記体動情報の複数周期となるように、前記同期撮像を制御し、
前記計測軌跡は、2つの線分軌跡が接続して成る折れ線計測軌跡であって、前記2つの線分軌跡間の角度が異なる複数の折れ線軌跡であることを特徴とする磁気共鳴イメージング装置。
【請求項１１】
請求項3乃至8のいずれか一項又は10記載の磁気共鳴イメージング装置において、
前記計測軌跡は、スパイラル計測軌跡であることを特徴とする磁気共鳴イメージング装置。
【請求項１２】
被検体の心電図に撮像シーケンスを同期させる同期撮像によって、3次元K空間(各軸をkr,k1,k2とする)内の計測軌跡に沿ってエコーデータを計測する計測ステップと、
前記計測されたエコーデータを用いて、前記被検体の血管画像を取得するステップと、
を有して成る血管画像撮像方法であって、
前記計測軌跡は、前記3次元K空間におけるリードアウト方向krに垂直な2方向k1,k2に非平行な計測軌跡であり、
前記撮像シーケンスは、繰り返し時間が前記心電図の複数周期であり、前記3次元K空間のリードアウト方向に垂直な2方向に対応する実空間の方向に、それぞれ、前記被検体の血流の核磁化の位相を変調するためのリフェーズ傾斜磁場又はディフェーズ傾斜磁場を印加するものであることを特徴とする血管画像撮像方法。
【請求項１３】
請求項12記載の血管画像撮像方法において、
前記計測軌跡は、一つの直線状軌跡を、前記3次元K空間のリードアウト方向を回転軸として回転させて得られる複数の直線状計測軌跡であり、
前記計測ステップでは、前記撮像シーケンスを繰り返して、異なる直線状軌跡に沿うエコーデータがそれぞれ計測されることを特徴とする血管画像撮像方法。
【請求項１４】
請求項13記載の血管画像撮像方法において、
リファレンススキャンにより前記被検体を事前に撮像してデータを取得するステップと、
前記リファレンススキャンにより取得されたデータに基づいて、前記撮像シーケンスを用いて前記血管画像を取得するための撮像パラメータ値を決定するステップと、
前記決定された撮像パラメータ値に基づいて前記撮像シーケンスを設定するステップと、
を更に備えることを特徴とする血管画像撮像方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、磁気共鳴イメージング(以下、「MRI」という)における、被検体の体動情報に同期して撮像を行い、取得する血管画像の画質を向上させる技術に関する。
【背景技術】
【０００２】
MRI装置は、核磁気共鳴(NMR)現象を利用して被検体の画像を得る計測装置であって、被検体に高周波磁場(以下、RFという)パルスを照射し、その応答として、被検体の組織を構成する原子核スピンが発生するNMR信号を計測する。そして計測したNMR信号に基づいて、被検体の頭部、腹部、四肢等の形態や機能を2次元的に或いは3次元的に画像化する。撮像の際には、傾斜磁場によって、NMR信号は、異なる位相エンコードやスライスエンコードが付与されると共に、周波数エンコードが付与されて、時系列データとして計測される。計測されたNMR信号は、2次元又は3次元フーリエ変換されることにより画像に再構成される。【０００３】
上記MRI装置を用いて、非造影で、即ち造影剤が被検体に投与されることなく、被検体の血管画像(MRA画像)を取得することが行われている。【０００４】
その方法の一つとして、特許文献1のように、時相検出手段により収集された被検体の心時相を表す信号に同期して、該信号から設定された遅延時間(Delay Time、DT)を空けて所定のスライスエンコード量分のエコー信号を高速スピンエコー(FSE)シーケンスで収集する動作を、複数心拍毎に繰り返す方法がある。例えば、遅延時間を収縮期に設定してエコー信号を収集すれば静脈が主に描出された静脈画像が得られ、拡張期に設定すれば動脈と静脈が共に描出された動静脈画像が得られる。また、これらの2つの画像データを差分することで動脈が主に描出された動脈画像が得られる。【０００５】
また、特許文献2では、位相エンコード方向に走行する血管の描出能を向上するために、位相エンコード方向にディフェーズ(Dephase)又はリフェーズ(Rephase)傾斜磁場パルスを印加する。リードアウト(Readout)方向にディフェーズ又はリフェーズ傾斜磁場パルスを印加することも可能であり、この場合はリードアウト方向に走行する血管の描出能が向上する。さらに、リードアウト方向、位相エンコード方向の両方にディフェーズ又はリフェーズ傾斜磁場パルスを印加することも可能である。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第4090619号公報
【特許文献２】特許第4309632号公報
【特許文献３】特開平7-284485号公報
【非特許文献】
【０００７】
【非特許文献１】Ｊ. Ｉ. Jackson  et. al., Selection of a Convolution Function for Fourier Inversion Using Gridding， IEEE Trans. Med. Imaging, vol. 10, pp. 473-478, 1991
【発明の概要】
【発明が解決しようとする課題】
【０００８】
特許文献1は、位相エンコード方向又はスライスエンコード方向、のどちらか1方向にT2減衰の影響が入ったエコーデータが取得され、そのエコーデータをフーリエ変換して画像化することで非造影MRA画像を取得する。そのため、どちらか1方向にT2減衰の影響が入ったエコーデータを画像化すると、位相エンコード方向又はスライスエンコード方向にボケた非造影MRA画像となり、診断に悪影響を与える可能性がある。【０００９】
また、特許文献2では、スライスエンコード方向にはディフェーズ又はリフェーズ傾斜磁場パルスを印加しないため、スライスエンコード方向に走行する血管の描出能については考慮されていない。【００１０】
そこで、本発明は、エコーデータにおけるT2減衰の影響があっても、血管のボケが抑制され描出能が向上された非造影MRA像を取得することを目的とする。
【課題を解決するための手段】
【００１１】
上記目的を達成するために、本発明は、3次元K空間(各軸をkr,k1,k2とする)におけるリードアウト方向krに垂直な2方向k1,k2に非平行な計測軌跡に沿ってエコーデータを計測する撮像シーケンスを、被検体の周期的体動情報に同期させて実行し、その際、撮像シーケンスの繰り返し時間を周期的体動情報の複数周期となるようにする。
【００１２】
具体的には、本発明のMRI装置は、被検体の周期的な体動についての体動情報を検出する体動情報検出部と、撮像シーケンスに基づいて、周期的な体動情報に同期させた同期撮像を実行して、3次元K空間データの計測を制御する計測制御部と、3次元K空間データを用いて被検体の血管画像を再構成する演算処理部と、を有し、撮像シーケンスは、3次元K空間(各軸をkr,k1,k2とする)におけるリードアウト方向krに垂直な2方向k1,k2に非平行な計測軌跡に沿ってエコーデータを計測するものであり、計測制御部は、撮像シーケンスの繰り返し時間が体動情報の複数周期となるように、同期撮像を制御することを特徴とする。
【００１３】
また、本発明の血管画像撮像方法は、被検体の心電図に撮像シーケンスを同期させる同期撮像によって、3次元K空間内の所定の計測軌跡に沿ってエコーデータを計測する計測ステップと、計測されたエコーデータを用いて、被検体の血管画像を取得するステップと、を有し、計測軌跡は、3次元K空間(各軸をkr,k1,k2とする)におけるリードアウト方向krに垂直な2方向k1,k2に非平行な計測軌跡であり、撮像シーケンスの繰り返し時間が心電図の複数周期であることを特徴とする。
【００１４】
また、本発明の血管画像撮像方法は、被検体の心電図に撮像シーケンスを同期させる同期撮像によって、3次元K空間内の所定の計測軌跡に沿ってエコーデータを計測する計測ステップと、計測されたエコーデータを用いて、被検体の血管画像を取得するステップと、を有し、計測軌跡は、3次元K空間におけるリードアウト方向に垂直な2方向に非平行な計測軌跡であり、撮像シーケンスの繰り返し時間(TR)が心電図の複数周期であることを特徴とする。
【発明の効果】
【００１５】
本発明のMRI装置及び血管画像撮像方法によれば、エコーデータにおけるT2減衰の影響があっても、血管のボケが抑制され描出能が向上された非造影MRA像を取得することができる。特に、リードアウト方向以外の方向(例えばリードアウト方向がH-F(Head-Foot)の場合は、R-L(Right-Left)とA-P(Anterior-Posterior))に走行する血管に対する描出能が向上された非造影MRA画像を取得できることができる。
【図面の簡単な説明】
【００１６】
【図１】本発明に係るMRI装置の一実施例の全体構成を示すブロック図。
【図２】実施例1における、(k1-k2)空間を非直交系サンプリングするための計測軌跡を示す図。
【図３】実施例1における、図2に示した計測軌跡に沿ってエコーデータの計測を行なう撮像シーケンスを表すシーケンスチャート。
【図４】実施例1の動作フローを表すフローチャート。
【図５】実施例1における、リファレンススキャンで用いる、PC(Phase Contrast)法パルスシーケンスの一例を示すシーケンスチャート。
【図６】撮像領域内の動脈と静脈の血流速度の変化と心電図との関係を表す流速変化グラフの一例を示す図。
【図７】心電図に同期する同期撮像であって、撮像シーケンスの繰り返し時間(TR)を複数の心拍(R-R)毎とする一例を示す図。(a)図は遅延時間(DT)が収縮期に設定されて静脈画像を取得する例を、(b)図は遅延時間(DT)が拡張期にされて動静脈画像を取得する例を、それぞれ示す。
【図８】実施例2の撮像シーケンスを示すシーケンスチャートであって、図3に示した実施例1の撮像シーケンスに、ディフェーズ傾斜磁場パルス又はリフェーズ傾斜磁場パルスを3方向に追加したもの。
【図９】実施例3における、(k1-k2)空間を非直交系サンプリングするための計測軌跡を示す図。
【図１０】実施例3における、図9に示した計測軌跡に沿ってエコーデータの計測を行なう撮像シーケンスを示すシーケンスチャート。
【図１１】実施例4における、(k1-k2)空間を非直交系サンプリングするための計測軌跡を示す図。
【図１２】実施例4における、図9に示した計測軌跡に沿ってエコーデータの計測を行なう撮像シーケンスを示すシーケンスチャート。
【図１３】実施例5における、(k1-k2)空間を非直交系サンプリングするための計測軌跡を示す図。
【図１４】実施例5における、図9に示した計測軌跡に沿ってエコーデータの計測を行なう撮像シーケンスを示すシーケンスチャート。
【図１５】実施例6における、(k1-k2)空間の格子点のエコーデータをジグザグに非直交系サンプリングするための一つの基本ジグザグ計測軌跡を示す図。(a)図は、(b)図よりk1(k2)方向の幅の狭い基本ジグザグ計測軌跡を示す。
【図１６】実施例6における、(k1-k2)空間の格子点のエコーデータをランダムに非直交系サンプリングを行なう計測軌跡を示す図。
【発明を実施するための形態】
【００１７】
以下、添付図面に従って本発明のMRI装置の好ましい実施形態について詳説する。なお、発明の実施形態を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。
【００１８】
最初に、本発明に係るMRI装置の一例の全体概要を図1に基づいて説明する。図1は、本発明に係るMRI装置の一実施例の全体構成を示すブロック図である。このMRI装置は、NMR現象を利用して被検体101の断層画像を得るもので、図1に示すように、静磁場発生磁石102と、傾斜磁場コイル103及び傾斜磁場電源109と、送信RFコイル104及びRF送信部110と、受信RFコイル105及び信号検出部106と、信号処理部107と、計測制御部111と、全体制御部108と、表示・操作部113と、被検体101を搭載してその被検体101を静磁場発生磁石102の内部に出し入れするベッド112と、を備えて構成される。
【００１９】
静磁場発生磁石102は、垂直磁場方式であれば被検体101の体軸と直交する方向に、水平磁場方式であれば体軸方向に、それぞれ均一な静磁場を発生させるもので、被検体101の周りに永久磁石方式、常電導方式あるいは超電導方式の静磁場発生源が配置されている。
【００２０】
傾斜磁場コイル103は、MRI装置の実空間座標系(静止座標系)であるX，Y，Zの3軸方向に巻かれたコイルであり、それぞれの傾斜磁場コイルは、それを駆動する傾斜磁場電源109に接続され電流が供給される。具体的には、各傾斜磁場コイルの傾斜磁場電源109は、それぞれ後述の計測制御部111からの命令に従って駆動されて、それぞれの傾斜磁場コイルに電流を供給する。これにより、X，Y，Zの3軸方向に傾斜磁場Gx，Gy，Gzが発生する。
【００２１】
2次元スライス面の撮像時には、スライス面(撮像断面)に直交する方向にスライス傾斜磁場パルス(Gs)が印加されて被検体101に対するスライス面が設定され、そのスライス面に直交して且つ互いに直交する残りの2つの方向に位相エンコード傾斜磁場パルス(Gp)と周波数エンコード(リードアウト)傾斜磁場パルス(Gf)が印加されて、エコー信号にそれぞれの方向の位置情報がエンコードされる。本発明に係る3次元領域の撮像時の制御については後述する。
【００２２】
送信RFコイル104は、被検体101にRFパルスを照射するコイルであり、RF送信部110に接続され高周波パルス電流が供給される。これにより、被検体101の生体組織を構成する原子の原子核スピンにNMR現象が誘起される。具体的には、RF送信部110が、後述の計測制御部111からの命令に従って駆動されて、高周波パルスが振幅変調され、増幅された後に被検体101に近接して配置された送信RFコイル104に供給されることにより、RFパルスが被検体101に照射される。
【００２３】
受信RFコイル105は、被検体101の生体組織を構成する原子核スピンのNMR現象により放出されるNMR信号(エコー信号)を受信するコイルであり、信号検出部106に接続されて受信したエコー信号が信号検出部106に送られる。信号検出部106は、受信RFコイル105で受信されたエコー信号の検出処理を行う。具体的には、RF送信コイル104から照射されたRFパルスによって誘起された被検体101の応答のエコー信号が被検体101に近接して配置された受信RFコイル105で受信され、後述の計測制御部111からの命令に従って、信号検出部106が、受信されたエコー信号を増幅し、直交位相検波により直交する二系統の信号に分割し、それぞれを所定数(例えば128,256,512等)サンプリングし、各サンプリング信号をA/D変換してディジタル量に変換し、後述の信号処理部107に送る。 従って、エコー信号は所定数のサンプリングデータからなる時系列のデジタルデータ(以下、エコーデータという)として得られる。
【００２４】
信号処理部107は、エコーデータに対して各種処理を行い、処理されたエコーデータを計測制御部111に送る。
【００２５】
計測制御部111は、被検体101の断層画像の再構成に必要なデータ収集のための種々の命令を、主に、傾斜磁場電源109と、RF送信部110と、信号検出部106に送信してこれらを制御する制御部である。具体的には、計測制御部111は、後述する全体制御部108の制御で動作し、ある所定のパルスシーケンスに基づいて、傾斜磁場電源109、RF送信部110及び信号検出部106を制御して、被検体101へのRFパルスと傾斜磁場パルスの印加及び被検体101からのエコー信号の検出を繰り返し実行し、被検体101の撮像領域についての画像の再構成に必要なエコーデータを収集する。
【００２６】
全体制御部108は、計測制御部111の制御、及び、各種データ処理と処理結果の表示及び保存等の制御を行うものであって、CPU及びメモリを内部に有する演算処理部114と、光ディスク、磁気ディスク等の記憶部115とを有して成る。具体的には、計測制御部111を制御してエコーデータの収集を実行させ、計測制御部111からのエコーデータが入力されると、演算処理部114がそのエコーデータに印加されたエンコード情報に基づいて、メモリのK空間に相当する領域に記憶させる。メモリのK空間に相当する領域に記憶されたエコーデータ群をK空間データともいう。そして演算処理部114はこのK空間データに対して信号処理やフーリエ変換による画像再構成等の処理を実行し、その結果である被検体101の画像を、後述の表示・操作部113に表示させると共に記憶部115に記録する。
【００２７】
表示・操作部113は、再構成された被検体101の画像を表示する表示部と、MRI装置の各種制御情報や上記全体制御部108で行う処理の制御情報を入力するトラックボール又はマウス及びキーボード等の操作部と、から成る。この操作部は表示部に近接して配置され、操作者が表示部を見ながら操作部を通してインタラクティブにMRI装置の各種処理を制御する。
【００２８】
また、本発明に係るMRI装置は、被検体の体動情報を検出する体動情報検出部を備える。この体動情報検出部は、被検体101に装着されて被検体の体動情報を検出するセンサー部116と、センサー部116からの信号を処理して、その処理した体動情報を計測制御部111に送る体動情報処理部117とを有してなる。体動情報検出部が被検体の心電図(心電波形)を検出するものであれば、センサー部116は心電図を検出する電極であり、体動情報処理部117は電極からのアナログ信号を処理する。計測制御部111は、体動情報検出部で検出された被検体の体動情報に同期させて、パルスシーケンスの実行による撮像を行なう同期撮像を制御する。
【００２９】
なお、図1において、送信側のRF送信コイル104と傾斜磁場コイル103は、被検体101が挿入される静磁場発生磁石102の静磁場空間内に、垂直磁場方式であれば被検体101に対向して、水平磁場方式であれば被検体101を取り囲むようにして設置されている。また、受信側の受信RFコイル105は、被検体101に対向して、或いは取り囲むように設置されている。
【００３０】
現在のMRI装置の撮像対象核種は、臨床で普及しているものとしては、被検体の主たる構成物質である水素原子核(プロトン)である。プロトン密度の空間分布や、励起状態の緩和時間の空間分布に関する情報を画像化することで、人体頭部、腹部、四肢等の形態または、機能を2次元もしくは3次元的に撮像する。
【００３１】
(本発明の計測軌跡の特徴について)
本発明は、3次元K空間(kr、k1、k2)において、リードアウト方向に対応するkr方向に直交する方向であるk1方向及びk2方向からなる(k1-k2)空間内で、(k1-k2)空間の座標軸(k1、k2)に非平行(非直交的)な計測軌跡に沿ってエコーデータを計測する(以下このような計測を非直交系(Non-Cartesian)サンプリングともいう)。計測軌跡に沿うエコーデータの計測は、該計測軌跡を等間隔又は不等間隔のいずれでも良い。その結果、K空間の格子点から外れたエコーデータが殆どとなる。これに対して、従来はK空間座標軸の一つに平行(直交系的)な計測軌跡に沿ってK空間格子点上のエコーデータの計測を行なう(以下このような計測を直交系(Cartesian)サンプリングともいう)。
【００３２】
(k1-k2)空間の座標軸(k1、k2)に非平行な計測軌跡の例としては、基本計測軌跡を、3次元K空間のリードアウト(kr)方向を回転軸として回転させた計測軌跡とすることができる。換言すれば、(k1-k2)空間内の任意の基準点(例えば原点又は原点近傍の点或いはこれら以外の任意の点)の周りで基本計測軌跡を回転させた計測軌跡とすることができる。或いは、(k1-k2)空間内をランダムに計測しても良い。
【００３３】
以上の様に本発明は(k1-k2)空間を非直交系サンプリングするので、従来の直交系サンプリングにおける位相エンコード方向とスライスエンコード方向の区別が全く無くなり、これらの両エンコード方向を区別して定義することができない。それ故に、本発明では、リードアウト方向に垂直な2方向をk1方向とk2方向とし、これらの2方向が張る空間を(k1-k2)空間とした。当業者には、直交系サンプリングにおいては、3次元K空間における位相エンコード方向とスライスエンコード方向を明確に意識して、それぞれkｐ方向及びkｓ方向と記載する習慣が有るが、本発明は非直交系サンプリングに基づき、位相エンコード方向とスライスエンコード方向の意識が全く無い。このような当業者の習慣に従うことは誤解を招くので、本明細書では、敢えて異なる記載を行なう。
【００３４】
なお、K空間の各座標軸と、実空間における各傾斜磁場の印加方向とはそれぞれ対応する。即ち、実空間におけるリードアウト傾斜磁場の印加方向がK空間のリードアウト方向に対応し、実空間におけるリードアウト傾斜磁場の印加方向に垂直な2方向が、K空間のリードアウト方向に垂直な2方向にそれぞれ対応する。以下の説明では、この対応に基づいて、適宜、2つの座標空間を区別することなく同じ表現で方向を指定する。
【００３５】
(本発明の撮像シーケンスの特徴について)
本発明の撮像シーケンスは、(k1-k2)空間を非直交系サンプリングするようにエンコード傾斜磁場の印加量を制御する。具体的には、(k1-k2)空間における任意の計測点(k1、k2)に対応するエンコード傾斜磁場の印加強度(G1,G2)を、矩形波で傾斜磁場パルスを印加する場合を想定すると、以下の(1)式で表すことができる。
【００３６】
G1  ＝  k1/(γ・FOV1・T)        (1)
G2  ＝  k2/(γ・FOV2・T)
ここでTはエンコード傾斜磁場の印加時間、FOV1はk1方向の視野サイズ、FOV2はk2方向の視野サイズ、γは磁気回転比を表す。即ち、本発明では、非直交サンプリングする(k1-k2)空間における任意の計測点のエコーデータを計測する際には、その計測点の座標に応じて(1)式で定まるエンコード傾斜磁場を印加して、該エコーデータを計測する。
【００３７】
以上説明したように、本発明では、(k1-k2)空間を非直交サンプリングするため、取得されたエコーデータにおけるT2減衰の影響が、リードアウト(kr)方向のみならず、k1方向及びk2方向にも分散することになる。つまり、T2減衰の影響が3次元K空間において3次元的に分散することになる。これにより、フーリエ変換して得られるMRA画像のボケを低減することができる。従来は、(k1-k2)空間を直交系サンプリングするため、T2減衰の影響が特定の方向に集中し、その方向にMRA画像がぼけてしまう課題があった。本発明ではこの課題が解消される。
【００３８】
また好ましくは、本発明は、MRA画像の画質を向上させるために、少なくともk1方向とk2方向にディフェーズ傾斜磁場又はリフェーズ傾斜磁場を印加する。好ましくはkr方向も含めた3方向にディフェーズ傾斜磁場又はリフェーズ傾斜磁場を印加する。これにより、MRA画像における血管の描出能を向上させることができる。
【００３９】
また好ましくは、本発明は、上記非直交系サンプリング行なう撮像シーケンスを、被検体の体動情報に同期させた同期撮像を行なうことで、MRA画像を取得する。そこで、撮像シーケンスの繰り返し時間(TR)を複数心拍とするが、これはT2強調(T2W)画像を取得するためである。血液のT1値を考慮すると、一般的な心拍数(60前後)の被検体では、2〜3心拍程度の繰り返し時間とすることが望ましい。
【００４０】
また、同期タイミング(例えば心電図のR波)から撮像シーケンスが起動されるまでの時間或いは実効TEで指定されるエコー信号のピーク位置までの時間(つまり実効TEまでの時間)である遅延時間(Delay Time;DT)を制御する。具体的には、遅延時間(DT)を収縮期に設定して静脈画像を取得する。また、遅延時間(DT)を拡張期に設定して動静脈画像を取得する。同期方法は心電同期又は脈波同期のいずれでも良い。以下、心電同期を例にして、遅延時間(DT)をR波から実効TEまでの時間として、本発明に係る同期非直交系サンプリングについての複数の実施例を説明する。
【００４１】
(実施例1)
次に、本発明のMRI装置及び血管画像撮像方法についての実施例1を説明する。本実施例は、一つの直線状軌跡を、3次元K空間のリードアウト(kr)方向を回転軸として回転させて得られる複数の直線状計測軌跡に沿ってエコーデータを計測する。即ち、本実施例は、直線状計測軌跡を基本計測軌跡とし、(k1-k2)空間内でリードアウト(kr)方向を回転軸としてこの直線状計測軌跡を回転させた計測軌跡を用い、このような回転対称な直線状計測軌跡に沿ってラディアル的な非直交系サンプリングを行なう。本実施例に係る撮像シーケンスは、これらの直線状計測軌跡に沿ってエコーデータを計測する。これにより、エコーデータにおけるT2減衰の影響をリードアウト方向以外の2方向にも分散させることができるため、MRA画像において血管のボケを抑制して血管の描出能を向上させることができる。
【００４２】
図2に、本実施例に係り、(k1-k2)空間を非直交系サンプリングするための計測軌跡を示す。図2は(k1-k2)空間のみを示しており、リードアウト(kr)方向はこの(k1-k2)空間に垂直な方向であって紙面に垂直な方向である。以降の説明においても同様である。
【００４３】
本実施例の計測軌跡は、(k1-k2)空間の原点又は原点近傍の点(任意の基準点)を通過する直線状の計測軌跡である。各直線状計測軌跡は、リードアウト(kr)方向を回転軸とする回転対称な関係にあり、任意の一つの直線状計測軌跡を(k1-k2)空間内の原点又は原点近傍の点の周りに所定角度回転して他の直線状計測軌跡が得られる。原点周りの回転であれば回転軸はリードアウト軸となり、任意の基準点周りの回転であれば回転軸は該基準点を通るリードアウト軸に平行な直線となる。別の表現をすれば、本実施例は、(k1-k2)空間内でラディアルサンプリングを行なう計測軌跡とする。そして、各直線状計測軌跡に沿って所定の間隔でエコーデータが計測される。各直線状計測軌跡に沿う複数のエコーデータが、撮像シーケンスの1繰り返し時間(TR)内で計測されるか、又は、複数の繰り返し時間に分割して計測される。
【００４４】
図2は、(k1-k2)空間の原点周りに直線状計測軌跡を回転させる例であって、原点を通る直線状計測軌跡201に沿って、等間隔又は不等間隔の7つの計測点(202-1〜202-7)のエコーデータが計測される例を示す。他の直線状計測軌跡に沿うエコーデータの計測も同様である。
【００４５】
また、本実施例は、上記(k1-k2)空間の原点を通過する直線状計測軌跡に沿うエコーデータの計測を、被検体の体動情報に同期させて行う。具体的には、撮像シーケンスの繰り返し時間(TR)を複数心拍にして、複数の繰り返し時間(TR)間で、直線状計測軌跡の回転角を変える。計測軌跡の回転のさせ方は、撮像シーケンスの繰り返し時間毎に所定角度ずつ回転させても良いし、ランダムに回転させても良い。撮像が終了した時点で、画像再構成に必要な3次元K空間データが取得されていれば良い。
【００４６】
(実施例1に係る撮像シーケンス)
次に、本実施例に係り、図2に示した(k1-k2)空間の原点を通過する直線状計測軌跡に沿うエコーデータの計測を行なう撮像シーケンスを図3に基づいて説明する。図3は、本実施例の撮像シーケンスのシーケンスチャート(タイミング図)であり、エコーファクター＝7の例、即ち、1繰り返し時間(TR)内或いは1回の励起で7つのエコー信号を計測する例を示す。ECG、RF/Echo、G1(G2)、G2(G1)、及びGrは、それぞれ心電図(心電波形)、RFパルス/エコー信号、k1(k2)方向に印加する傾斜磁場パルス波形、k2(k1)方向に印加する傾斜磁場波形、及びリードアウト方向に印加する傾斜磁場波形を意味する。以下、後述する各シーケンスチャートにおいて同様とする。なお、k1方向とk2方向は、特に区別する必要が無く、いずれがk1方向又はk2方向で良いので、k1(k2)方向とk2(k1)方向の様に記載している。
【００４７】
計測制御部111は、図3に示すシーケンスチャートに基づいて、傾斜磁場電源109、RF送信部110及び信号検出部106を制御してエコー信号の計測を行なう。具体的には、90°RFパルス301がスライス選択傾斜磁場310と共に印加されて所望の撮像領域(後述のステップ401で設定される撮像領域)が励起される。なお、図3では、心電図R波を簡略表示してある。その後、所定の時間間隔で180°反転RFパルス302が複数回印加される。図3の例はエコーファクター＝7なので、7回の180°反転RFパルス(302-1〜302-7)が印加される。そして180°反転RFパルス302の後でエコー信号303が発生し、7つのエコー信号(303-1〜303-7)がそれぞれ計測される。この際、心電図でR波からの遅延時間(DT)は、R波から実効TEまでの時間である(つまり、R波から、K空間中心に配置されるエコー信号303-4のピーク位置までの時間)。また、180°反転RFパルスは、スライス選択又はスライス非選択のどちらでもよい。スライス非選択にすることでプロファイルの改善、血管描出能の改善が見込める可能性がある。図3はスライス非選択の例を示す。図3以降のシーケンスチャートも同様である。
【００４８】
各エコー信号303には、(k1-k2)空間内の直線軌跡上のエコーデータとするために、k1方向とk2方向にエンコード傾斜磁場304と306がそれぞれ印加される。各エコー信号303の計測の後には、印加されたエンコード傾斜磁場304,306の印加量をキャンセルして横磁化の位相を元のゼロに戻すために、エンコード傾斜磁場と逆極性で絶対値が同じ印加量(＝印加波形と時間軸との囲む面積)を有するリワインド傾斜磁場305と307がk1方向とk2方向にそれぞれ印加される。これらのエンコード傾斜磁場304,306及びリワインド傾斜磁場305,307の印加量は、(k1-k2)空間内の計測点の位置に応じて(1)式に基づいて変更される。図3の例では、エンコード傾斜磁場304,306は、エコー信号の計測毎に負極性大振幅→振幅0(ゼロ)→正極性大振幅と変更され、反対にリワインド傾斜磁場305,307は、エコー信号の計測毎に正極性大振幅→振幅0(ゼロ)→負極性大振幅と変更される。これにより、図2に示す一つの直線状計測軌跡201上において、一方の端部から原点を通って他方の端部迄の各計測点のエコーデータがそれぞれ計測される。図2の例では、直線状計測軌跡201上の7つの計測点202-1〜202-7のエコーデータが計測されることになる。
【００４９】
k1方向とk2方向のエンコード傾斜磁場304,306の振幅比は、(k1-k2)空間内での直線状計測軌跡の回転角度に応じて変わる。即ち、回転角ゼロで定まる最大振幅のエンコード傾斜磁場を、直線状計測軌跡の回転角度に応じて、k1方向とk2方向のエンコード傾斜磁場にそれぞれ分配される。したがって、直線状計測軌跡の回転角度を変えて計測する際には、これらのk1方向とk2方向に分配されるエンコード傾斜磁場の印加量が回転角度に応じて変更される。つまり、図3のように、計測される各エコー信号に付与されるk1、k2方向のエンコード量を回転角度に対応して共に変化させることで、図2のような直線状計測軌跡を(k1-k2)空間の原点の周りで回転させながら、ラディアル的なk空間走査を行なうことが可能になる。
【００５０】
なお、各エコー信号303の計測の際には、リードアウト傾斜磁場309が印加されて、リードアウト方向の空間位置情報が各エコー信号303にエンコードされる。また、各エコー信号をそのサンプリング期間の略中央でピークとなるように計測するために、各リードアウト傾斜磁場309の半分の印加量となるディフェーズ傾斜磁場308がリードアウト傾斜磁場309の前であってスライス選択傾斜磁場310の後に印加される。
【００５１】
計測制御部111は、以上の一繰り返し時間(TR)分のパルスシーケンスを、複数心拍ごとに繰り返す。その際、複数の繰り返し時間(TR)の間で、直線状計測軌跡の回転角を変える。即ち、上述したように、エンコード傾斜磁場及びリワインド傾斜磁場の印加量を変えて、直線状計測軌跡の回転角度を異ならせる。好ましくは、繰り返し時間(TR)毎に直線状計測軌跡の回転角度を異ならせる。異ならせ方は、前述したとおり、所定の一定角度毎でも良いし、ランダムであっても良い。或いは、一つの直線状計測軌跡に沿うエコーデータの計測を複数の繰り返しで分割して行っても良い。
【００５２】
以上の本実施例の撮像シーケンスによれば、リードアウト方向に垂直な(k1-k2)空間内でエコーデータをラディアル的に計測することで、計測されたエコーデータにおけるT2減衰の影響を3次元的に(少なくともk1方向とk2方向に)分散できるようになり、フーリエ変換して得られるMRA画像における血管のボケを低減することができる。
【００５３】
(リファレンススキャン)
また、本実施例は、上記撮像シーケンスを用いて、被検体に造影剤を投与することなく、即ち非造影で、MRA画像を取得する。そこで、計測制御部111は、リファレンススキャンを行い、上記撮像シーケンスを用いて所望の非造影MRA画像の取得に好適な撮像パラメータの値を決定するためのデータを取得する。リファレンススキャンは撮像シーケンスの実施前に実施される。撮像シーケンスは、このリファレンススキャンで取得されたデータを用いて決定された好適な撮像パラメータ値に基づいて実施される。
【００５４】
撮像シーケンスにより取得される非造影MRA画像の画質を左右する撮像パラメータは、心周期における拡張期と収縮期の時相を特定するためのＲ波からの遅延時間(DT)と、撮像領域における血管内の血流速度と、を含むので、決定すべき撮像パラメータ値としては、その遅延時間(DT)血流速度の値である。
【００５５】
例えば、リファレンススキャンとして、図5に示すPC(Phase Contrast)法パルスシーケンスを用いる場合は、所望の(後述するステップ401で設定される)撮像領域を時系列で複数撮像し、該撮像領域の位相画像を時系列で複数取得する。そして、得られた時系列位相画像を用いて、注目する血流部分の位相の時間変化を求めて、該血流の流速変化グラフを作成する。位相画像の位相と血流速度との関係は以下の(2)式に基づいて求めることができる。
0005858791.tif 000002
ここで、γはラーモア周波数、Gは傾斜磁場強度、vは血流速度、tは傾斜磁場印加時間である。
【００５６】
図6に流速変化グラフの一例を示す。この流速変化グラフは、心電図と同じ時間スケールで並列に表示したものである。実線が動脈血流の流速変化グラフであり、点線が静脈の流速変化グラフである。これにより、心電時相と血流速度との関係が一目瞭然となり、拡張期又は収縮期を特定するための遅延時間(DT)と、撮像領域における血管内の血流速度の設定が容易となる。設定は、例えば、操作者が表示・操作部113を介して、この流速変化グラフ上で所望の点又は期間を指定することにより行なわれる。
【００５７】
一般的には、図6に示す様に、心電図の特徴箇所にP、Q、R、S、Tの符号が付されており、R-T期間が収縮期であり、T-R期間が拡張期である。そこで、R-T期間のいずれかを計測期間として設定すると、その期間のエコーデータを用いて収縮期画像(静脈画像)が取得され、T-R期間のいずれかを計測期間として設定すると、その期間のエコーデータを用いて拡張期画像(動静脈画像)が取得さる。
【００５８】
また、リファレンススキャンとして、本撮像シーケンスに準ずるパルスシーケンスでリファレンススキャンを行う場合は、操作者は、得られた画像における動脈、静脈の描出能を視覚的に評価して、遅延時間(DT)や血流速度を決定し、表示・操作部113を介して、それらの値を入力設定する。
【００５９】
(実施例1に係る処理フロー)
次に、本実施例の、(k1-k2)空間内計測軌跡に沿うエコーデータの計測を行なう撮像シーケンスを用いた非造影MRA画像の取得を実現する処理フローを図4に基づいて説明する。図4は、本実施例の処理フローを表すフローチャートである。この動作フローの全体フロー及び各ステップにおける個別処理はプログラムとして予め磁気ディスク等の記憶部115に記憶されており、CPUが必要に応じてメモリに読み込んで実行することにより実施される。以下、各ステップを詳細に説明する。なお、非造影MRA画像を取得するので、被検体に造影剤を投与するステップは無い。
【００６０】
ステップ401で、操作者は、表示・操作部113を介して、撮像シーケンスの撮像条件(撮像領域、FOV，リードアウト方向、画像のマトリックス数等)を設定する。特にリードアウト方向の設定は、(k1-k2)空間をラディアル的にサンプリングするため、位相エンコード方向とスライスエンコード方向の区別が無くなるので、リードアウト方向のみを設定する。リードアウト方向は、H-F(Head-Foot)、R-L(Right-Left)、A-P(Anterior-Posterior)のいずれか1方向に実質的に一致させることが好ましい。さらに、リードアウト方向を血管の走行方向に合わせることが望ましい。例えば、下肢の非造影MRA画像を取得したい場合には、下肢の血管の走行方向は主にH-F方向となるので、リードアウト方向をH-F方向に合わせることが望ましい。
【００６１】
また、操作者は、取得するMRA画像種として動静脈画像(拡張期画像)か静脈画像(収縮期画像)かを決定する。この決定に基づいて、後述するステップにおける画像演算の手法が設定される。
【００６２】
ステップ402で、計測制御部111は、ステップ401で設定された撮像領域に対し、リファレンススキャンを実行する。リファレンススキャンにより計測されたエコーデータ又は画像は、後述のステップにおいて、撮像シーケンスを用いた所望の非造影MRA画像の取得に好適な撮像パラメータ値を決定するために用いられる。
【００６３】
リファレンススキャンに用いるパルスシーケンスとしては、前述したとおり、図5に示すような、ベロシティーエンコード(VENC)パルスを用いる公知のPC法に基づくパルスシーケンスでもよし、撮像シーケンスに準ずるパルスシーケンスでもよい。詳細は前述したとおりである。
【００６４】
ステップ403で、演算処理部114は、ステップ402のリファレンススキャンで計測されたデータ(エコーデータ又は画像データ)に基づいて、撮像シーケンスを用いた所望の非造影MRA画像の取得に好適な撮像パラメータ値を決定する。導出すべき撮像パラメータ値については前述したとおりである。そして、演算処理部114は、ステップ401で設定された撮像条件とステップ403で決定された撮像パラメータ値に基づいて、前述した図3の撮像シーケンスを具体的に設定する。
【００６５】
ステップ404で、計測制御部111は、ステップ403で具体的に設定された撮像シーケンスを用いて、本実施例の非直交系サンプリングを例えば被検体から検出した心電図に同期させて、同期撮像(本撮像)を開始する。
【００６６】
同期撮像の際には、計測制御部111は、撮像シーケンスの繰り返し時間(TR)を2以上の複数心拍とし、複数の繰り返し時間(TR)の間で直線状計測軌跡の回転角を変えるために、直線状計測軌跡の回転角に応じたエンコード傾斜磁場の印加強度や印加量を変更して印加する。
【００６７】
また、ステップ401で設定されたMRA画像種とステップ403で決定された遅延時間(DT)に基づいて、計測制御部111は、R波からの遅延時間を設定する。具体的には、静脈画像を取得する場合には心電図R波からの遅延時間(DT)を収縮期に設定し、動静脈画像を取得する場合には心電図R波からの遅延時間(DT)を拡張期に設定する。
【００６８】
計測制御部111は、以上のような、繰り返し時間(TR)を複数心拍とする心電同期非直交系サンプリングを行なう撮像シーケンスを設定して開始する。
【００６９】
図7は、心電図に同期する同期撮像であって、撮像シーケンスの繰り返し時間(TR)を複数の心拍(R-R)毎とする一例として、TR＝3心拍(3R-R)の場合を示す。また、図7は、遅延時間(DT)の後の黒枠の期間に撮像シーケンスが実行される例を示す。図7(a)は、遅延時間(DT)が収縮期に設定されて静脈画像を取得する例を、図7(b)は遅延時間(DT)が拡張期にされて動静脈画像を取得する例を、それぞれ示す。
【００７０】
ステップ405で、計測制御部111は、ステップ404で設定し開始した、繰り返し時間(TR)を複数心拍とする心電同期非直交系サンプリングを行なう撮像シーケンスを繰り返し実行して、所望の非造影MRA画像を取得するためのエコーデータを計測する。その際、計測制御部111は、直線状計測軌跡が、3次元K空間のリードアウト(kr)方向を回転軸として、(k1-k2)空間内で回転するように各エンコード傾斜磁場の出力(印加強度、印加量)を制御し、各直線状計測軌跡に沿うエコーデータを計測する。好ましくは、計測制御部111は、繰り返し時間(TR)毎に計測軌跡の角度を変える。そして、各直線状計測軌跡に沿うエコーデータの計測の際に、計測制御部111は、撮像シーケンスを心電図Ｒ波に同期させて、ステップ404で設定された遅延時間(DT)を空けてこのような非直交系サンプリングを行なう。
【００７１】
ステップ406で、計測制御部111は、ステップ401で設定された撮像条件に基づくエコーデータ量、即ち、画像再構成に必要なエコーデータ量の計測が完了したか否かを判定し、否(No)であれば、ステップ405に戻って、本実施例の非直交系サンプリングを継続する。その際、計測制御部111は、直線状計測軌跡の回転軸の周りの回転角を変えて、異なる直線状計測軌跡に沿ってエコーデータの計測を継続する。完了(Yes)であればステップ407に移行する。
【００７２】
ステップ407で、演算処理部114は、ステップ405で計測されたエコーデータを、3次元K空間の各格子点に再配置(グリッディング)する。グリッディング処理は、例えばSinc関数やKaiser-Bessel関数の補間用関数を用いて行う(非特許文献1)。
【００７３】
ステップ408で、演算処理部114は、ステップ407でグリッディングされた3次元K空間データをフーリエ変換して3次元画像を再構成する。そして、ステップ401で設定された画像種に応じて、演算処理部114は、静脈画像(収縮期画像)と動静脈画像(拡張期画像)間で各種演算を行う。例えば、ステップ401で動脈画像の取得が設定されていたならば、収縮期画像と拡張期画像との間で差分演算を施して、その差分演算の結果取得された3次元画像を3次元動脈画像とする。
【００７４】
ステップ409で、演算処理部114は、ステップ408で演算の結果取得した3次元画像データを用いて、所望の方向の投影画像を作成し、最終的な非造影MRA画像とする。投影画像を作成する処理としては、例えば公知のMIP(Maximum Intensity Projection)法やボリュームレンダリング法を用いることができる。
【００７５】
以上までが本実施例の処理フローの説明である。なお、リファレンススキャンを実行することなく、所定の、或いは、事前に決定された撮像パラメータ値を用いてもよい。この場合は、これらの事前の撮像パラメータ値を用いて撮像シーケンスを生成しておくことになり、前述のステップ402と403を実行する必要はない。
【００７６】
以上説明したように、本実施例のMRI装置及び血管画像撮像方法は、3次元K空間のリードアウト方向を回転軸とする複数の直線状計測軌跡に沿ってエコーデータを計測する撮像シーケンスを、被検体の体動情報に同期させて行なう。その際、撮像シーケンスの繰り返し時間(TR)を複数心拍にして、複数の繰り返し時間(TR)の間で、直線状計測軌跡の(k1-k2)空間内での回転角を変える。これにより、計測されたエコーデータにおけるT2減衰の影響が3次元的に分散するため、このように計測されたエコーデータをフーリエ変換して得られる非造影MRA画像のボケを低減することができ、画質を向上させることができる。
【００７７】
(実施例2)
本発明のMRI装置及び血管画像撮像方法についての実施例2を説明する。本実施例は、3次元K空間のリードアウト方向に垂直な2方向に対応する実空間の方向に、それぞれ、前記被検体の血流の核磁化の位相を変調するためのリフェーズ傾斜磁場又はディフェーズ傾斜磁場を印加する。具体的には、前述の実施例1の撮像シーケンスにおいて、3軸方向の少なくともk1方向とk2方向に、リフェーズ傾斜磁場又はディフェーズ傾斜磁場を印加する。好ましくはリードアウト(kr)方向を含めた全方向に印加する。例えば、H-F方向をリードアウト方向とすると、H-F方向、R-L方向、及びA-P方向にリフェーズ傾斜磁場又はディフェーズ傾斜磁場を印加する。これにより、MRA画像において、リフェーズ傾斜磁場又はディフェーズ傾斜磁場を印加した方向に走行する血管に対する描出能を向上させることができる。他については、前述の実施例1と同じなので説明を省略する。以下、図8に基づいて、本実施例を詳細に説明する。
【００７８】
図8は、本実施例の撮像シーケンスのシーケンスチャートを示し、前述の実施例1で説明した図3に示した撮像シーケンスに、ディフェーズ傾斜磁場パルス又はリフェーズ傾斜磁場パルスを3方向に追加したものである。
【００７９】
具体的には、各方向の傾斜磁場に追加された黒塗矩形の傾斜磁場パルス(801,802,,810,811)が、ディフェーズ傾斜磁場パルス又はリフェーズ傾斜磁場パルスに相当する(以下まとめて位相制御傾斜磁場という)。具体的には、リードアウト傾斜磁場(Gr)の印加方向に関しては、傾斜磁場パルス810,309,811の印加量比が略1:2:1の関係となるようにすれば、血流の核磁化の位相をリフェーズすることができる。逆に傾斜磁場パルス810,309,811の印加量比が略1:2:1の関係からずれるようにすれば、血流の核磁化の位相をディフェーズすることができる。一方、G1(G2)、及びG2(G1)傾斜磁場の印加方向に関しては、傾斜磁場パルス801,802はディフェーズ傾斜磁場パルスとなり血流の核磁化の位相をディフェーズすることができる。計測制御部111は、図8の撮像シーケンスを実施するよう、各方向の傾斜磁場印加を制御する。
【００８０】
上述の通り、k1(k2)方向のエンコード傾斜磁場G1(G2)と、k2(k1)方向のエンコード傾斜磁場G2(G1)には、90°パルス後に、位相制御傾斜磁場801,802がそれぞれ印加される。しかし、これらの位相制御傾斜磁場801,802によって、エコー信号に印加されるべき所望のエンコード量が影響を受ける。そこで、k1(k2)方向とk2(k1)方向では、位相制御傾斜磁場801,802の印加量に相当する分を、元々の各エンコード傾斜磁場304,306とリワインド傾斜磁場305,307の印加量に対して加算又は減算して得た印加量を、本実施例のエンコード傾斜磁場804,806とリワインド傾斜磁場805,807とする。図8の点線枠で示した部分が加算又は減算分に相当する。具体的には、エンコード傾斜磁場G1(G2)では、元々のエンコード傾斜磁場304の印加量から位相制御傾斜磁場801の印加量相当分を減算したエンコード傾斜磁場804が印加される。また、元々のリワインド傾斜磁場305の印加量に位相制御傾斜磁場801の印加量相当分を加算したリワインド傾斜磁場805が印加される。同様に、エンコード傾斜磁場G2(G1)では、元々のエンコード傾斜磁場306の印加量から位相制御傾斜磁場802の印加量相当分を減算したエンコード傾斜磁場806が印加される。また、元々のリワインド傾斜磁場307の印加量に位相制御傾斜磁場802の印加量相当分を加算したリワインド傾斜磁場807が印加される。
【００８１】
また、リードアウト(kr)方向では、エコー信号の計測に際に印加するリードアウト傾斜磁場309の前後に位相制御傾斜磁場810,811がそれぞれ印加される。これにより、本来印加すべきディフェーズ傾斜磁場308の印加量が影響を受けるので、追加された位相制御傾斜磁場810,811の印加量に相当する印加量(点線枠で示した部分)を元々のディフェーズ傾斜磁場308から減算したディフェーズ傾斜磁場808が印加される。図8の例では、元々のディフェーズ傾斜磁場308の印加量と減算分とが相殺して、ディフェーズ傾斜磁場808の印加量が0(ゼロ)となった場合を示す。
【００８２】
以上の血流の核磁化の位相を変調するためのリフェーズ傾斜磁場パルス又はディフェーズ傾斜磁場パルスは、3軸方向の少なくともk1方向とk2方向に印加する。好ましくはリードアウト(kr)方向を含めた全方向に印加する。これにより、リフェーズ傾斜磁場パルス又はディフェーズ傾斜磁場パルスが印加された方向に走行する血管の描出能を向上させることができる。例えば、リードアウト方向をH-F方向としてK1方向とK2方向とをそれぞれR-L方向とA-P方向とした場合には、血管画像において、R-L方向とA-P方向に走行する血管の描出能を向上することができる。
【００８３】
また、ディフェーズ傾斜磁場パルス又はリフェーズ傾斜磁場パルスの印加量を、前述の実施例1で説明したラディアル的な非直交系サンプリングの場合に、k空間上で直線状計測軌跡の方向により変化させてもよい。例えば、直線状計測軌跡の方向がk1方向又はk2方向に一致している場合にのみ、ディフェーズ傾斜磁場パルス又はリフェーズ傾斜磁場パルスを印加してもよい。
【００８４】
以上の本実施例の撮像シーケンスを用いて、(k1-k2)空間を非直交系サンプリングすることや、被検体の体動情報に同期させた同期撮像を行なうことは前述の実施例1と同じなので詳細な説明は省略する。なお、本実施例は非造影MRA画像の取得に好適であり、血管描出能の向上効果が非造影MRA画像において顕著となる。
【００８５】
以上説明したように、本実施例のMRI装置及び血管画像撮像方法は、少なくともk1方向とk2方向にディフェーズ傾斜磁場パルス又はリフェーズ傾斜磁場パルスを印加してMRA画像を取得するので、そのMRA画像において、ディフェーズ傾斜磁場又はリフェーズ傾斜磁場を印加した方向に走行する血管の描出能力を向上させることができる。特に、リードアウト(kr)方向をH-F方向とした場合には、MRA画像において、リードアウト(kr)方向ではない方向(R-L方向とA-P方向)に走行する血管の描出能を向上させることができる。特に非造影MRA画像の取得において、本実施例の効果が顕著となる。
【００８６】
(実施例3)
本発明のMRI装置及び血管画像撮像方法についての実施例3を説明する。本実施例は、複数の平行直線状計測軌跡を含んで成る単位軌跡群(ブレード)を、リードアウト(kr)方向を回転軸として回転させて得られる複数の単位軌跡群に沿ってエコーデータを計測する非直交系サンプリングを行う。即ち、本実施例は、ブレードを構成する複数の平行直線状計測軌跡を基本計測軌跡とし、ブレードを(k1-k2)空間内の任意の基準点の回りで回転させた計測軌跡とする。したがって、前述の実施例1と異なる箇所は撮像シーケンス及び計測軌跡の形状である。他については前述の実施例1と同じなので詳細な説明を省略する。以下、本実施例の撮像シーケンス及び計測軌跡の形状について詳細に説明する。
【００８７】
図9は、本実施例の(k1-k2)空間における非直交系計測軌跡の一例を示し、図10は、図9の計測軌跡に沿ってエコーデータの計測を行なう撮像シーケンスを表すシーケンスチャートを示す。以下、このようなパルスシーケンスをハイブリッドラディアルシーケンスという。
【００８８】
図9に示す(k1-k2)空間内の計測軌跡は、複数の平行直線状計測軌跡を含んで成るブレードを、リードアウト(kr)方向を回転軸として回転させたものであり、回転後の各ブレードを構成する複数の平行直線状計測軌跡に沿ったエコーデータがそれぞれ計測される。なお、図9は、(k1-k2)空間の原点でブレードを回転させた例を示すが、回転中心は原点以外の原点近傍の点又は任意の基準点でよい。
【００８９】
具体的には、前述の実施例1で説明した図2に示した一つの直線状計測軌跡と、これに平行な複数の直線状計測軌跡と、を含んで成る平行直線状計測軌跡群をブレードとする。このブレードを(k1-k2)空間内でリードアウト(kr)方向を回転軸として回転させて複数のブレードを生成する。各ブレードを構成する直線状計測軌跡の数、ブレードの個数及び回転角度は、所望の画像が取得されるように決定することができ、例えば図4のステップ401で操作者が設定可能である。図9はブレードを構成する平行直線状計測軌跡の数が3つ(901,902,903)の例を示している。これらのブレードを構成する複数の平行直線状計測軌跡に沿ったエコーデータが計測される。
【００９０】
なお、本実施例の非直交系計測軌跡では、グリッディング処理が必要となるため、本実施例のフローチャートは実施例1のフローチャートである図2と同様となる。
【００９１】
図10に示すハイブリッドラディアルシーケンスは、一つのブレードにおいて複数の平行直線状計測軌跡に沿ってエコー信号を計測するために、エンコード傾斜磁場G1(G2)とG2(G1)には、90°RFパルス後に、オフセット傾斜磁場1001、1002が印加される。このオフセット傾斜磁場1001、1002は、平行直線状計測軌跡の数とそれらの間隔に応じて、印加量及び印加ステップ数が異なり、(3)式に応じて決定することができる。
【００９２】
G  ＝    k/(γ・FOV・T)        (3)
ここで、kはオフセット傾斜磁場のステップ番号、Tはオフセット傾斜磁場の印加時間、FOVはオフセット傾斜磁場印加方向の視野サイズ、γは磁気回転比を表す。
【００９３】
計測制御部111は、この(3)式で回転角度毎に求められるオフセット傾斜磁場1001、1002を、当該回転角のブレードを構成する複数の平行直線状計測軌跡に沿ったエコーデータの計測を行なう際に、印加する。他の傾斜磁場については、前述の実施例1で説明した図3の撮像シーケンスと同様なので詳細な説明は省略する。
【００９４】
上記ハイブリッドラディアルシーケンスは、(k1-k2)空間における回転の基準点(図9の場合は原点)近傍のデータを特に密に、或いは、重複させて計測することができるので、被検体の体動にロバストといわれている。
【００９５】
なお、本実施例のハイブリッドラディアルシーケンスにおいても、実施例1で説明した同期撮像や、実施例2で説明したディフェーズ傾斜磁場とリフェーズ傾斜磁場の印加を行って、MRA画像を取得することは同様であり、詳細な説明は省略する。そして、これらを組み合わせて非造影MRA画像の取得を行えば、実施例1と同様の効果が得られる。
【００９６】
以上説明したように、本実施例のMRI装置及び血管画像撮像方法は、複数の平行直線状計測軌跡を含んで成るブレードを、リードアウト(kr)方向を回転軸として回転させて得られる複数の平行直線状計測軌跡に沿ってエコー信号を計測する非直交系サンプリングを行う。これにより、前述の実施例1と同様の効果を得ることができる。さらに、このような計測軌跡を用いるので被検体の体動にロバストなMRA画像を取得することが可能になる。
【００９７】
(実施例4)
次に、本発明のMRI装置及び血管画像撮像方法についての実施例4を説明する。本実施例は、2つの線分軌跡が(k1-k2)空間の原点又は原点近傍の点で接続して成る計測軌跡(以下、折れ線計測軌跡という)とし、2つの線分間の角度が異なる複数の折れ線計測軌跡に沿ったエコーデータを計測する非直交系サンプリングとする。即ち、本実施例は、折れ線計測軌跡を単位計測軌跡とし、折れ線計測軌跡を折れ曲がり角度を変えて、(k1-k2)空間内の任意の基準点の回りで回転させた計測軌跡とする。したがって、前述の実施例1と異なる箇所は撮像シーケンス及び計測軌跡の形状である。他については前述の実施例1と同じなので詳細な説明を省略する。以下、本実施例の撮像シーケンス及び計測軌跡の形状について詳細に説明する。
【００９８】
図11は、本実施例の、(k1-k2)空間における折れ線計測軌跡の一例を示し、図12は、図11の折れ線計測軌跡に沿ってエコーデータの計測を行なう撮像シーケンスを表すシーケンスチャートを示す。
【００９９】
図11に示す(k1-k2)空間内の任意の一つの計測軌跡は、(k1-k2)空間の原点又は原点近傍の点で、角度θで折れ曲がり、結果として2つの線分が原点又は原点近傍の点で接続して成る折れ線計測軌跡である。つまり、(k1-k2)空間の端部(高域側)から中央部(低域)に向かう線分軌跡と、中央部(低域)から(k1-k2)空間の別の端部(高域側)に向かう線分軌跡と、が角度θを成して原点又は原点近傍の点で接続する折れ線計測軌跡である。そして、接続点が同一であって2つの線分間の角度が異なる複数の折れ線計測軌跡で(k1-k2)空間を充填する。折れ線計測軌跡の番号をkとすると、折れ曲がり角θとの関係を、例えば(4)式のようにすることができる。
【０１００】
θ＝(π/12)Ｘ(2k-1)      (k＝1，2，3・・11)      (4)
k＝1の計測軌跡の前半部を図11の線分軌跡1101-1とすると、その後半部の線分軌跡1102-1との折れ曲がり角θは、θ＝π/12となる。同様に、k＝2の計測軌跡の前半部を図11の線分軌跡1101-2とすると、その後半部の線分軌跡1102-2との折れ曲がり角θは、θ＝π/4となる。以下同様に、計測軌跡の前半部の線分軌跡1101を反時計回りに順次選択していくと、各計測軌跡の後半部の線分軌跡1102は時計回りとなり、これら2つの線分軌跡間の角度が順次開いていくことになる。したがって、これら全ての折れ線計測軌跡で、(k1-k2)空間を充填することができる。図11及び(4)式は、11本の折れ線軌跡で(k1-k2)空間を充填する例を示す。なお、図11の折れ線計測軌跡は、全て(k1-k2)空間原点で折れ曲がっているが、複数の折れ線計測軌跡の折れ曲がり点を異ならせても良い。
【０１０１】
このような折れ線計測軌跡の前半線分と後半線分上のエコーデータは、(k1-k2)空間原点に関して互いに複素共役でない関係となる。一般的にK空間原点に関して複素共役関係にあるエコーデータは実質的に同じ情報量を持つ。つまり、本実施例の折れ線計測軌跡に沿うエコーデータ群は、実施例1の直線状計測軌跡に沿うエコーデータ群よりもより多くの情報量を持っていることになる。したがって、本実施例の折れ線計測軌跡は多くの情報量を短時間で取得でき、K空間の非対称計測・再構成に好適である。本実施例の折れ曲がり計測軌跡を用いて非対称計測・再構成を行なう場合は、複数の折れ線計測軌跡の折れ曲がり点を異ならせて(k1-k2)空間原点近傍の低域エコーデータを対称に取得すると良い。
【０１０２】
なお、本実施例の非直交系計測軌跡では、グリッディング処理が必要となるため、本実施例のフローチャートは実施例1のフローチャートである図2と同様となる。
【０１０３】
図12に示す、上記計測軌跡に沿ったエコーデータの計測を行なう撮像シーケンスは、前述の実施例1で説明した図3に示した撮像シーケンスと比較して、点線枠部分のエンコード傾斜磁場パルス1206及びリワインド傾斜磁場パルス1207の極性及び振幅が異なる。これは、前述の図3の撮像シーケンスは、直線状計測軌跡に沿ってエコーデータを計測するものであるため、計測軌跡が(k1-k2)空間の第3象限から第1象限に伸びる直線である場合に、エンコード傾斜磁場及びリワインド傾斜磁場も、それらの振幅が単調に増加又は減少することになる。これに対して、図12の撮像シーケンスは、折れ線計測軌跡に沿ってエコーデータを計測するものであるため、折れ線計測軌跡の前半部に沿うエコーデータの計測に際には、図3の撮像シーケンスと同様に、エンコード傾斜磁場パルスとリワインド傾斜磁場パルスは、それらの振幅が単調に増加又は減少することになる。しかし、折れ線計測軌跡の後半部に沿うエコーデータの計測の際には、計測軌跡が折れ曲がっているために、エンコード傾斜磁場1206とリワインド傾斜磁場1207は、前半部に引き続いて単調に増加又は減少するものとはならず、増加又は減少の仕方が反転する。即ち、図12のシーケンスチャートの点線枠部分に着目すると、図3では、エンコード傾斜磁場パルスは正極性で振幅が単調増加しているのに対して、図12では、エンコード傾斜磁場パルス1206が負極性で単調減少している。振幅の変化の仕方は、折れ線の折れ曲がり角度に対応して変化する。また、リワインド傾斜磁場パルス1207は、エンコード傾斜磁場パルス1206と逆極性の変化となる。エンコード傾斜磁場及びリワインド傾斜磁場のこのような変化により、折れ線軌跡に沿ったエコーデータの計測が可能になる。
【０１０４】
なお、本実施例の撮像シーケンスにおいても、実施例1で説明した同期撮像や、実施例2で説明したディフェーズ傾斜磁場とリフェーズ傾斜磁場の印加を行って、MRA画像を取得することは同様であり、詳細な説明は省略する。そして、これらを組み合わせて非造影MRA画像の取得を行えば、実施例1と同様の効果が得られる。
【０１０５】
以上説明したように本実施例のMRI装置及び血管画像撮像方法は、(k1-k2)空間内の原点又は原点近傍の点で折れ曲がり、かつ、折れ線がなす角度が計測軌跡毎に異なるような折れ線計測軌跡に沿ってエコーデータを計測する非直交系サンプリングを行なうので、前述の実施例1と同様の効果を得ることができる。さらに、折れ線計測軌跡とすることで、多くの情報量を短時間で取得できる効果を得ることができる。
【０１０６】
(実施例5)
次に、本発明のMRI装置及び血管画像撮像方法についての実施例5を説明する。本実施例は、(k1-k2)空間の原点又は原点近傍の点を通過する渦巻状(スパイラル状)の計測軌跡とする。そして、一つのスパイラル軌跡を、リードアウト(kr)方向を回転軸として回転させて得た複数のスパイラル軌跡に沿ってエコーデータを計測する非直交系サンプリングを行なう。即ち、本実施例は、スパイラル計測軌跡を単位計測軌跡とし、このスパイラル計測軌跡を(k1-k2)空間内の原点又は原点近傍の点の回りで回転させた計測軌跡とする。したがって、前述の実施例1と異なる箇所は撮像シーケンス及び計測軌跡の形状である。他については前述の実施例1と同じなので説明を省略する。以下、本実施例の撮像シーケンス及び計測軌跡の形状について詳細に説明する。
【０１０７】
図13は、本実施例の(k1-k2)空間におけるスパイラル計測軌跡の一例を示し、図14は、図13のスパイラル計測軌跡に沿ってエコーデータの計測を行なう撮像シーケンスを表すシーケンスチャートを示す。
【０１０８】
図13に示す(k1-k2)空間内の任意の一つの計測軌跡は、(k1-k2)空間の原点又は原点近傍の点を通過する渦巻状(スパイラル状)の軌跡である。そして、一つのスパイラル計測軌跡を、リードアウト(kr)方向を回転軸として所定角度回転させて複数のスパイラル計測軌跡を得る。図13は、実線で表したスパイラル計測軌跡と、これを原点周りに回転させて得た点線のスパイラル計測軌跡を示す。各スパイラル計測軌跡の原点又は原点近傍の点回りの回転数は任意に設定可能であり、(k1-k2)空間の原点の周りを1周以下の回転数とする計測軌跡でも良いし、複数周回転する計測軌跡でも良い。
【０１０９】
これらのスパイラル計測軌跡上を中心(低域)から端部(高域)に向けて、或いは逆に端部(高域)から中心(低域)に向けて順にエコーデータを計測しても良いし、スパイラル計測軌跡上をランダムに移動してエコーデータを計測してもよい。図13は、1301-1〜1301-7の順で、つまり一つのスパイラル軌跡を端部側から中心に向けて、さらに同じ軌跡上を中心から端部側に戻る順序で、エコーデータを計測する例を示す。
【０１１０】
このようなスパイラル計測軌跡は、一つの計測軌跡で(k1-k2)空間を略均等に走査できるので、計測されたエコーデータにおけるT2減衰の影響や、体動の影響をより3次元的に分散できるようになるので、MRA画像におけるボケや体動アーチファクトの影響を低減でき、画質を向上させることができる。
【０１１１】
なお、本実施例の非直交系計測軌跡では、グリッディング処理が必要となるため、本実施例のフローチャートは実施例1のフローチャートである図2と同様となる。
【０１１２】
図14に示す、上記スパイラル計測軌跡に沿ったエコーデータの計測を行なう撮像シーケンスは、前述の実施例1で説明した図3に示した撮像シーケンスと比較して、エンコード傾斜磁場G1(G2)とG2(G1)が異なる。図3に示した撮像シーケンスにおけるエンコード傾斜磁場G1(G2)とG2(G1)は、直線状計測軌跡に沿ったエコーデータを計測するためには、各エンコード傾斜磁場パルスの振幅又は印加量が単調に増加又は減少するものであった。これに対して、本実施例のエンコード傾斜磁場G1(G2)とG2(G1)は、スパイラル計測軌跡に沿ってエコーデータを計測するために印加される波形であるため、前述の(1)式に基づいて複雑に変化する。例えば、図13のスパイラル計測軌跡上の計測点1301-1〜1301-7のエコーデータが図14のシーケンスチャートにおけるエコー信号303-1〜303-7にそれぞれ対応するので、これらのエコー信号303-1〜303-7を計測するために印加されるエンコード傾斜磁場パルス1404,1406及びリワインド傾斜磁場パルス1405,1407は、図14に示すとおりとなる。なお、スパイラル計測軌跡の詳細は、例えば特許文献3に説明されており、詳細な説明は省略する。
【０１１３】
なお、本実施例の撮像シーケンスにおいても、実施例1で説明した同期撮像や、実施例2で説明したディフェーズ傾斜磁場とリフェーズ傾斜磁場の印加を行って、MRA画像を取得することは同様であり、詳細な説明は省略する。そして、これらを組み合わせて非造影MRA画像の取得を行えば、実施例1と同様の効果が得られる。
【０１１４】
以上説明したように、本実施例のMRI装置及び血管画像撮像方法は、(k1-k2)空間内でその原点又は原点近傍の点を通過するスパイラル測軌跡に沿ってエコーデータを計測する非直交系サンプリングを行なうので、前述の実施例1と同様の効果を得ることができる。さらに、スパイラル計測軌跡とすることで、MRA画像におけるボケや体動アーチファクトの影響を低減でき、画質を向上させることができる。
【０１１５】
(実施例6)
次に、本発明のMRI装置及び血管画像撮像方法についての実施例6を説明する。本実施例は、3次元K空間におけるリードアウト方向に垂直な2方向において、複数の格子点をジグザグ又はランダムに通過する計測軌跡とする。即ち、本実施例は、(k1-k2)空間内で複数の格子点をジグザグ又はランダムに通過する計測軌跡に沿ってエコーデータを計測する非直交系サンプリングを行なう。したがって、前述の実施例1と異なる箇所は撮像シーケンス及び計測軌跡の形状とグリッディングが不要となる点である。他については前述の実施例1と同じなので説明を省略する。以下、本実施例の撮像シーケンス及び計測軌跡の形状について詳細に説明する。
【０１１６】
図15，16は、本実施例の、(k1-k2)空間の格子点をジグザグ又はランダムに移動しながらそのエコーデータを計測するための計測軌跡の一例を示す。
【０１１７】
図15に示す計測軌跡は、(k1-k2)空間の格子点のエコーデータをジグザグに非直交系サンプリングするための一つの基本ジグザグ計測軌跡の一例である。この基本ジグザグ計測軌跡を原点又は任意の基準点の周りに所定の角度で回転させながら、最寄りの格子点のエコーデータを計測する。好ましくは、撮像シーケンスの繰り返し時間(TR)毎に基本ジグザグ計測軌跡を所定の角度で回転させながら最寄りの格子点のエコーデータを計測する。図15(a)は、k1(k2)方向の幅の狭い基本ジグザグ計測軌跡の一例であって、図の点線矢印方向に各格子点のエコーデータを計測する例を示す。図15(b)は、図15(a)の基本ジグザグ計測軌跡よりもk1(k2)方向の幅が広い基本ジグザグ計測軌跡の例を示す。(k1-k2)空間内で回転させることは図15 (a)と同じである。或いは、前述の実施例1〜実施例4で説明した各非直交系計測軌跡に最寄りの格子点のエコーデータを計測するようにしてもよい。
【０１１８】
図16は、(k1-k2)空間の格子点のエコーデータをランダムに非直交系サンプリングを行なう例であり、点線矢印の順序で各格子点のエコーデータを計測する例を示す。各格子点のエコーデータをランダムに計測するために、計測制御部111は、擬似乱数を発生させて格子点を決定し、決定した格子点位置に応じて、撮像シーケンスのエンコード傾斜磁場G1、G2の印加量を(1)式に基づいて決定する。そして、計測制御部111は、撮像シーケンスの繰り返し時間毎に、計測する格子点の位置を変える。好ましくは、計測格子点の位置が重複しないように、繰り返し時間毎の計測格子点の位置を異ならせる。また、好ましくは、実効TEのエコーデータだけはK空間中心データとなるようにエンコード傾斜磁場の印加を制御する。
【０１１９】
上記のように、(k1-k2)空間の格子点をジグザグ又はランダムに移動しながらそのエコーデータを計測するための、本実施例に係る撮像シーケンスは、前述の実施例1で説明した図3の撮像シーケンスと比較してエンコード傾斜磁場G1(G2)とG2(G1)が異なる。具体的には、各エンコード傾斜磁場パルス及び各リワインド傾斜磁場パルスの振幅又は印加量を、ジグザグ又はランダムに計測する格子点毎に、該格子点の座標値に応じて(1)式に基づいて設定すればよい。その結果、各エンコード傾斜磁場パルス及び各リワインド傾斜磁場パルスの振幅又は印加量は不規則になる。詳細な図示及び説明を省略する。
【０１２０】
なお、本実施例の撮像シーケンスにおいても、実施例1で説明した同期撮像や、実施例2で説明したディフェーズ傾斜磁場とリフェーズ傾斜磁場の印加を行って、MRA画像を取得することは同様であり、詳細な説明は省略する。そして、これらを組み合わせて非造影MRA画像の取得を行えば、実施例1と同様の効果が得られる。
【０１２１】
以上のように、ジグザグ又はランダムいずれの場合もK空間の格子点位置のエコーデータが計測されるので、画像再構成する際にグリッディング処理を行う必要が無くなる。そのため、演算処理部114は、実施例1で説明した図4の処理フローにおいて、ステップ407のグリッディング処理を行わずにスキップして、ステップ408で計測されたK空間データをそのままフーリエ変換する。
【０１２２】
以上説明したように、本実施例のMRI装置及び血管画像撮像方法は、(k1-k2)空間内で、その格子点のエコーデータをジグザグ又はランダムとなるような計測軌跡に沿ってエコー信号を計測する非直交系サンプリングを行なうので、前述の実施例1と同様の効果を得ることができる。さらに、格子点のエコーデータを直接計測するのでグリッディング処理が必要ない。そのため、グリッディング処理の必要が無い本実施例においては、画像再構成処理を簡略・短時間化することが可能になる。
【０１２３】
以上が、本発明を適用する具体的な実施例である。しかし、本発明は、以上の各実施例で開示された内容にとどまらず、本発明の趣旨を踏まえた上で各種形態を取りえる。例えば、各実施例で説明した計測軌跡を混合した非直交系サンプリングを行なうことができる。必要なことは、計測されたエコーデータにおけるT2減衰の影響が3次元的に分散さるような非直交系計測軌跡とすることである。
【０１２４】
また、本発明の各実施例の説明では、被検体に造影剤を投与しないで、非造影MRA画像を取得する例を説明したが、造影剤を投与して造影MRA画像を取得する場合にも、本発明の各実施例を適用することにより、画質を向上させることができる。
【符号の説明】
【０１２５】
101  被検体、102  静磁場発生磁石、103  傾斜磁場コイル、104  送信RFコイル、105  受信RFコイル、106  信号検出部、107  信号処理部、108  全体制御部、109  傾斜磁場電源、110  RF送信部、111  計測制御部、112  ベッド、113  表示・操作部、114  演算処理部、115  記憶部、116  センサー部、117  体動情報処理部
【図１】
0005858791.tif 000003
【図２】
0005858791.tif 000004
【図３】
0005858791.tif 000005
【図４】
0005858791.tif 000006
【図５】
0005858791.tif 000007
【図６】
0005858791.tif 000008
【図７】
0005858791.tif 000009
【図８】
0005858791.tif 000010
【図９】
0005858791.tif 000011
【図１０】
0005858791.tif 000012
【図１１】
0005858791.tif 000013
【図１２】
0005858791.tif 000014
【図１３】
0005858791.tif 000015
【図１４】
0005858791.tif 000016
【図１５】
0005858791.tif 000017
【図１６】
0005858791.tif 000018