(57)【特許請求の範囲】
【請求項１】
被検体に核磁気共鳴を起させる手段と、前記核磁気共鳴信号に位相エンコードを付与する手段と、前記核磁気共鳴信号を測定し、前記核磁気共鳴信号に基づき画像を再構成する手段と、前記各手段を制御し、位相エンコードで決まるｋ空間の計測を制御する制御手段とを備え、
前記制御手段は、前記ｋ空間上に、ｋ空間中央領域を含む必須計測領域と、前記必須計測領域と重複した領域を持たない周辺計測領域とを設定し、前記ｋ空間にデータとして配置される核磁気共鳴信号の計測を繰り返し行い、各回の計測が、前記必須計測領域のデータと、前記周辺計測領域から選択した一部のデータの計測を含み、各回毎に選択する部分を順次変更するように制御する磁気共鳴イメージング装置において、
前記制御手段は、前記周辺計測領域に存在する全サンプル点を、その空間的分布がほぼ均等な複数の群に分け、前記計測の各回毎に選択する群を順次変更するように制御することを特徴とする磁気共鳴イメージング装置。
【請求項２】
請求項１に記載の磁気共鳴イメージング装置において、
前記画像再構成手段は、今回の計測ステップで計測しなかった群のデータを、過去の計測ステップで計測された群のデータから生成し、周辺計測領域に補充し、この補充されたｋ空間のデータに基づいて画像を生成することを特徴とする磁気共鳴イメージング装置。
【請求項３】
請求項１または２に記載の磁気共鳴イメージング装置において、
前記ｋ空間は３次元のｋ空間であって、前記必須計測領域と周辺計測領域の設定は、当該３次元のｋ空間全体を、前記３次元のｋ空間を規定する３つの座標軸のうちの少なくとも２つの座標軸で決まる２次元空間について前記必須領域と周辺領域を設定することを特徴とする磁気共鳴イメージング装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、核磁気共鳴現象を利用して被検体の断層画像を得る磁気共鳴イメージング装置(ＭＲＩ)に関し、特に、磁気共鳴イメージング装置おけるダイナミック磁気共鳴アンジオグラフィ撮像に好適な計測方法に関するものである。
【０００２】
【従来の技術】
ＭＲＩ（磁気共鳴イメージング装置）は、ＮＭＲ現象を利用して、被検体中の所望の検査部位における原子核スピンの密度分布、緩和時間分布等を計測して、その計測データから被検体の任意の断面を画像表示するものである。
【０００３】
また、このような磁気共鳴イメージング装置を利用して血流画像を得る技術は、ＭＲＡ（磁気共鳴アンジオグラフィ）と呼ばれる。このＭＲＡは、造影剤を使用しない手法と造影剤を使用する手法に大別することができる。
【０００４】
造影剤を用いるＭＲＡの手法としては、Gd-DTPA等のT1短縮型の造影剤とグラジエントエコー系のTRの短いパルスシーケンスを組み合わせる手法が一般的である。この手法では、数ms〜数10msの短時間間隔の高周波磁場による同一領域の励起の繰り返しによって、血流以外の部分のスピンを飽和させ、T1の短い造影剤を含むためにスピンの飽和が起こりにくい血流部分からの信号を高信号で得る。造影剤が血管内にとどまっている間に血管を含む体積の計測を行うことにより血管を他組織に対し高コントラストで描出することができる。さらに、造影前後の画像差分をとるによって、血流部分のコントラストを上げた画像を得る手法DSA（Digital Substraction Angiography）も行われている。
【０００５】
このような造影剤を用いるＭＲＡの手法では、通常肘静脈から造影剤を注入する。注入された造影剤は、心臓、動脈系、毛細血管、静脈系と順に循環する。そこで、ダイナミックＭＲＡと呼ばれる技術では、造影剤の循環の各段階において繰り返し計測を行い、各部の血流の時系列画像が得られるようにしている。このような造影剤を用いるＭＲＡの手法やダイナミックＭＲＡについては、例えば「3D Contrast ＭＲＡngiography 2nd edition. Prince MR, Grist TM and Debatin JF, Springer, pp3-39,1988」に詳細に記載されている。
【０００６】
このようなダイナミックＭＲＡにおける1回の撮像時間は、２次元計測の場合は繰り返し時間TRと位相エンコードステップ数を乗じた時間となり、３次元計測の場合は繰り返し時間TRと位相エンコードステップ数にスライスエンコードステップ数を乗じた時間となる。したがって、空間分解能を向上するためには位相エンコードステップ数やスライスエンコードステップ数を大きくとることが望ましい一方、そのようにすると時間分解能は低下する。すなわち基本的に空間分解能と時間分解能はトレードオフの関係にある。
【０００７】
そこで、このような問題を解決するために、計測の繰り返しにおいて前の回で計測したデータを共有するＭＲフルオロスコピーと呼ばれる手法や、全エンコードステップについては計測したデータをレファレンスとして、以後は一部の領域（低周波領域）のみを計測し更新するキーホール法と呼ばれる手法がある。さらに米国特許5,713,358号、米国特許5,830,143号には、次のような計測方法を提案している。
【０００８】
すなわち、これらの技術では、一回の計測においてｋ空間全てを計測するのではなく、ｋy-kz空間をky方向に複数の領域に分割し、分割した領域を単位として、ｋ空間の中心領域（低周波数領域）をより高い頻度で、周辺領域をより低い頻度で計測するように、各回で計測するｋ空間の部分を制御する。そして、今回計測しなかったｋ空間の部分については、それ以前に計測した回の計測結果を流用する。または、それ以前に計測した複数の回の計測結果から補間によって生成する。
【０００９】
この技術によれば、各回の計測において位相エンコードステップの全てを実行しないために、各回の計測時間を短縮し時間分解能を向上することができると共に、比較的、診断に重要なコントラストを決定するｋ空間の中心領域（低周波領域）については経時的な変化を確実に捉えることができる。
【００１０】
【発明が解決しようとする課題】
しかしながら、前記米国特許に記載される技術では、高周波領域のデータについては低い頻度でしか計測されないため、高周波数領域の時間分解能の犠牲が著しく大きかった。また、ｋ空間の分割の仕方がky方向のみと一次元方向のみであるために、計測結果に方向依存性が生じ、これによって生成される画像の品質を劣化する場合があった。
【００１１】
そこで、本発明は、ＭＲＩにおいて繰り返し撮像を行う際に、高周波数領域の時間分解能を大きく劣化させることなく、各回の計測時間を短縮化する計測方法を提供することを課題とする。また、このような計測方法において、さらに、計測結果の方向依存性を排除することを課題とする。
【００１２】
【課題を解決するための手段】
上記課題を達成するために、本発明は、ＭＲＩ装置を用いて被検体の計測を繰り返し行う計測方法であって、ｋ空間上に、ｋ空間中央領域を含む必須計測領域と、各々前記必須計測領域と重複した領域を持たない周辺計測領域とを設定する初期設定ステップと、周辺計測領域に存在する全サンプル点を、その空間的分布がほぼ均等な複数の群に分けるステップと、前記ｋ空間に配置されるデータとして前記被検体からの核磁気共鳴信号を計測する計測ステップと、前記計測ステップを繰り返し行う繰り返しステップとを有し、前記繰り返しステップによって繰り返し行われる各回の計測ステップは、前記必須計測領域のデータと、前記周辺計測領域の複数の群から選択した一つの群のデータの計測を含み、各回毎に選択する群を順次変更することを特徴とするＭＲＩ装置における計測方法を提供する。
【００１３】
ここで「サンプル点の空間的分布がほぼ均等な群」とは、各群のサンプル点が、それぞれ空間的に偏りがない（ｋ空間の中心に対し、ほぼ等方的に分布する）状態で存在し、かつ各群が互いに同等であることを意味する。
【００１４】
本発明のＭＲＩにおける計測方法は、さらに、計測ステップで計測したデータを用いて画像再構成するステップを含み、当該画像再構成するステップは、今回の計測ステップで計測しなかった群のデータを、他の回（過去の回）の計測ステップで計測された群のデータを用いて生成し、周辺計測領域に補充し、この補充されたｋ空間のデータに基づいて画像を生成することを特徴とする。
【００１５】
このような計測方法によれば、各回の計測ステップにおける測定されるデータに、高周波領域の成分として、必ず周辺計測領域分の一つの群のデータが含まれることになるため、前記従来のｋ空間分割計測手法（米国特許5,713,358号、米国特許5,830,143号等）に比べ、高周波数領域の時間分解能の劣化は小さく、生成画質の劣化も小さい。しかも各回の周辺計測領域のデータは、空間的に均等であるので、計測結果に方向依存性が生じることがない。
【００１６】
本発明の計測方法は、好適には３次元のｋ空間の計測に適用される。この場合、必須領域と周辺領域を設定するステップにおいて、当該３次元のｋ空間全体を、前記３次元のｋ空間を規定する３つの座標軸のうちの少なくとも２つの座標軸で決まる２次元空間について必須計測領域と周辺計測領域を設定する。
一般に３次元ｋ空間の計測法においては、スライスエンコード方向（ｋｚ軸）のサンプル数および／または位相エンコード方向（ｋｙ軸）のエンコードステップを間引いたときに計測時間の低減効果を得ることができる。従って、本発明の計測方法においてもこれらに２軸で決まる空間について、周波数領域の群の分割を行うことにより、計測時間低減による時間分解能の向上を高めることができる。
【００１７】
また、このような計測方法において、ｋ空間の全領域を計測するのではなく、その一部を非計測とすることができる。即ち、必須計測領域および周辺計測領域の和集合空間に外接する直方体状のｋ空間は、和集合空間と重複しない、非計測領域を有するものとすることができる。この場合、非計測領域についてはゼロデータを埋めて画像再構成する。このように、直方体状のｋ空間のうちの、計測目的上重要でない空間周波数部分を非計測領域とすることにより、計測目的への適合度を大きく損なうことなく、時間分解能を向上することができる。
【００１８】
本発明の計測方法では、さらに、周辺測定領域を低中周波領域と高周波領域に分割し、低中周波領域と高周波領域で、分割する群の数を異ならせることも可能である。例えば低中周波領域の分割数を少なくすることにより、低中周波領域で計測するサンプル点の密度を高くし、この領域が画像コントラストに寄与する部分の画質の劣化を防止することができる。
【００１９】
本発明のＭＲＩ装置は、被検体に核磁気共鳴を起させる手段と、前記核磁気共鳴信号に位相エンコードを付与する手段と、前記核磁気共鳴信号を測定し、前記核磁気共鳴信号に基づき画像を再構成する手段とを備えた磁気共鳴イメージング装置であって、前記各手段を制御し、位相エンコードで決るｋ空間の計測を制御する制御手段と、ｋ空間上の、ｋ空間中央領域を含む必須計測領域と、前記必須計測領域と重複した領域を持たない周辺計測領域との設定および前記周辺計測領域のサンプル点を空間的分布が均等である複数の群に分割した群の設定を記憶する記憶手段とを有し、制御手段は、被検体の計測を繰り返し行う制御を行うと共に、繰り返し行う各回の計測において、前記必須計測領域のデータと、計測の各回毎に順次選択した一つの周辺計測領域の群のデータとを計測し、今回の計測で計測しなかった群のデータを、他の回の計測で計測された群のデータを用いて生成し、周辺計測領域に補充することを特徴とする。
【００２０】
【発明の実施の形態】
以下、本発明の実施形態について、ダイナミックＭＲＡへの適用を例にとり説明する。
【００２１】
図1に、本実施形態に係るＭＲＩ装置の構成を示す。
図示するように、このＭＲＩ装置は、ＮＭＲ現象を利用して被検体の断層像を得るもので、静磁場発生磁石2と、傾斜磁場発生系3と、シーケンサ4と、送信系5と、受信系6と、信号処理系7と、中央処理装置(CPU)8とを備えて成る。
【００２２】
静磁場発生磁石2は、被検体1の周りにその体軸方向または体軸と直交する方向に均一な静磁場を発生させるもので、被検体1の周りのある広がりをもった空間に永久磁石方式または常電導方式あるいは超電導方式の磁場発生手段が配置されている。
【００２３】
傾斜磁場発生系3は、X,Y,Zの三軸方向に巻かれた傾斜磁場コイル9と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源10とから成り、後述のシーケンサ4からの命令に従ってそれぞれのコイルの傾斜磁場電源10を駆動することにより、X,Y,Zの三軸方向の傾斜磁場Gx,Gy,Gzを被検体1に印加するようになっている。この傾斜磁場の加え方により被検体1の特定のスライス又はスラブを選択的に励起することができ、また計測空間（ｋ空間）における計測点（サンプル点）の位置、計測順序を規定することができる。
【００２４】
シーケンサ4は、CPU8の制御で動作し、被検体1の断層像のデータ収集に必要な種々の命令を、傾斜磁場発生系3、送信系5及び受信系6に送るようになっている。シーケンサ4が制御する傾斜磁場発生系3、送信系5及び受信系6の動作タイミングはパルスシーケンスと呼ばれ、ここではパルスシーケンスの一つとして三次元血流撮像のためのシーケンスが採用される。このパルスシーケンスはCPU8に備えられたメモリに予めプログラムとして組み込まれており、他のパルスシーケンスと同様、使用者が撮影の目的に応じて適宜選択することにより実行することができる。シーケンサ4の制御については後に詳述する。
【００２５】
送信系5は、シーケンサ4から送り出される高周波パルスにより被検体1の生体組織を構成する原子の原子核に核磁気共鳴を起こさせるために高周波磁場を照射するもので、高周波発振器11と変調器12と高周波増幅器13と送信側の高周波コイル14aとから成る。高周波発振器11から出力された高周波パルスをシーケンサ4の命令にしたがって変調器12で振幅変調し、この振幅変調された高周波パルスを高周波増幅器13で増幅した後に被検体1に近接して配置された高周波コイル14aに供給することにより、電磁波が被検体1に照射されるようになっている。
【００２６】
受信系6は、被検体1の生体組織の原子核の核磁気共鳴により放出されるエコー信号(ＮＭＲ信号)を検出するもので、受信側の高周波コイル14bと増幅器15と直交位相検波器16とA/D変換器17とから成る。上記送信側の高周波コイル14aから照射された電磁波による被検体1の応答の電磁波(ＮＭＲ信号)は被検体1に近接して配置された高周波コイル14bで検出される。検出されたエコー信号は、増幅器15及び直交位相検波器16を介してA/D変換器17に入力されディジタル量に変換され、さらにシーケンサ4からの命令によるタイミングで直交位相検波器16によりサンプリングされた二系列の収集データとされ、信号処理系7に送られる。
【００２７】
信号処理系7は、CPU8と、磁気ディスク18及び磁気テープ19等の記録装置と、CRT等のディスプレイ20とから成り、CPU8で３次元フーリエ変換、補正係数計算像再構成等の処理を行い、任意断面の信号強度分布あるいは複数の信号に適当な演算を行って得られた分布を画像化してディスプレイ20に断層像として表示するようになっている。なお、図1において、送信側及び受信側の高周波コイル14a、14bと傾斜磁場コイル9は、被検体1の周りの空間に配置された静磁場発生磁石2の磁場空間内に設置されている。
【００２８】
次に、このようなＭＲＩ装置における、ダイナミックＭＲＡ撮像動作について説明する。
まず、CPU8の入力装置を介して造影剤を用いたＭＲＡが選択されると、シーケンス4はCPU8によって制御され、三次元ＭＲＡシーケンスを実行する。
【００２９】
このパルスシーケンスは、例えば図２に示すように、グラディエントエコー法を基本とするシーケンスで、三次元ＭＲＡシーケンスに一般的なものである。ただし、Gradient Moment Nulling 等のフローコンペンセーションなどが付加される場合もある。
【００３０】
図示した例では、高周波磁場パルスRFを印加して目的血管を含む領域（スラブ）を励起した後、Z方向の位相エンコード傾斜磁場GzおよびY方向の位相エンコード傾斜磁場Gyを印加し、次いで読み出し/周波数エンコード傾斜磁場Gxを印加するとともにその極性を反転させてエコー信号を計測している。３次元イメージングでは、このようなパルスシーケンスを、Z方向の位相エンコード傾斜磁場Gz、Y方向の位相エンコード傾斜磁場Gyの磁場強度の組み合わせを変えながら、所定の繰り返し時間TRで繰り返すことによって、三次元データを得る。
【００３１】
ここで、Z、Y方向のエンコードステップ数はZ、Y方向の領域の分割数を決めるもので、例えばそれぞれ128、256などに設定されている。また、このZ、 Y方向のエンコードステップ数によってｋ空間のky−kz平面のマトリクスサイズが規定される。即ち、図２のシーケンスにおいて、傾斜磁場強度Gzのある値と傾斜磁場強度Gyのある値の組み合わせのときに計測された信号は、Gz、Gyの値に対応するｋ空間の座標（ky,kz）に配置されることになる。逆に言えば、ｋ空間のky−kz平面のある座標（ky1、kz1）をサンプルするということは、ky1に対応するGyとkz1に対応するGzの組み合わせにおいて、図２のパルスシーケンスを実行することであり、座標（ky1、kz1）をサンプルしないということは、ky1に対応するGyとkz1に対応するGzの組み合わせにおいて、図２のパルスシーケンスの実行を省略するということである。
【００３２】
従来の基本的なダイナミックＭＲＡでは、各回の計測におけるZ、Y方向のエンコードステップ数、ステップ幅を固定として、毎回の計測において、全てのGyとGzの組み合わせについて図２のパルスシーケンスを実行し、全ての座標にデータが充足された３次元のｋ空間を得る。
【００３３】
これに対し、本実施形態では、ダイナミックＭＲＡ撮像において、各回の計測におけるZ、Y方向のエンコードステップ数、ステップ幅を固定せずに、各回の計測を以下のように制御する。
【００３４】
すなわち、本実施形態では、X方向の周波数エンコード数、Z, Y方向の位相エンコードステップ数に応じて大きさを定めたkx, ky, kzの座標軸を有する３次元ｋ空間のky-kz平面上に、原点を中心とする略円形の中央領域である必須計測領域と、その周辺に周辺計測領域を設定する。そして周辺計測領域については、周辺計測領域に存在する全てのサンプル点を複数の群に分割する。その際、分割された各群において、サンプル点は空間的な偏りがないよう分布し、かつ各群が互いに同等であるように分割する。このようなサンプル点の設定分割は、後述する計測の順序とともに、CPU８のメモリに記憶される。
【００３５】
各回の計測において、必須計測領域は必ず計測し、周辺計測領域の各群については所望の順序で一つずつ計測する。従って各群の数と同じ数の計測を終了した時点で、周辺計測領域の全てのサンプル点についてのデータを得ることになる。
【００３６】
図３および図４に、第1の実施形態として周辺計測領域を３つの群に分割した場合を示す。ここでは、説明を簡単にするために、マトリックスサイズが５×９であるｋｙ−ｋｚ空間について説明する。
【００３７】
本例では、図３に示すように、黒丸で示す、座標原点を含む１１のサンプル点が必須計測領域であり、必須計測領域を除く他の領域が周辺計測領域である。周辺計測領域の３４のサンプル点は、図３(a)、(b)、(c)に斜線が付された丸で示すように１２、１１、１１に３分割されており、それぞれ空間的にほぼ均等に分散している。
【００３８】
最初の計測では、まず図４に示すように、必須計測領域の１１のサンプル点と(a)に示す群(以下、群aとする)の計測を行い、必須計測領域と周辺計測領域の１２のサンプル点のデータを収集する。次に必須計測領域の１１のサンプル点と(b)に示す群bの計測を行い、必須計測領域と周辺計測領域の群aとは異なる１１のサンプル点のデータを収集する。さらに３回目の計測では、必須計測領域必須計測領域の１１のサンプル点と周辺計測領域のうち群a、b以外の１１のサンプル点のデータを収集する。
【００３９】
この３回の計測で周辺計測領域について全てのデータが揃うので、ここでｋ空間の全領域のデータを用いて画像を再構成する。すなわち、３次元のｋ空間に対して３次元フーリエ変換を施し、３次元データを得、これを投影した２次元画像を得る。
【００４０】
４回目の計測では、1回目の計測と同様に必須計測領域と群aについての計測を行い、必須計測領域のデータと周辺計測領域の一部のデータを得る。この回で得られなかった周辺計測領域のデータについては、既に２回目および３回目の計測で得た周辺計測領域のデータを用い、必須計測領域と全ての周辺計測領域のデータが充填された３次元データを作成し、この３次元データを用いて画像を再構成する。
【００４１】
５回目以降は、同様に、その回で得た必須計測領域のデータと周辺計測領域の一部のデータと、その回に最も近い過去２回の計測で得た周辺計測領域のデータを用いて、３次元データを作成し、画像再構成する。
【００４２】
このようにすることで、1回の計測時間を、ｋ空間の全領域のデータを計測する場合に比べ、大幅に(前掲の例で約1／2程度)短縮することができ、その1回の計測時間と同じ間隔で画像を再構成することができるので、ダイナミックＭＲＡで表示される画像の時間分解能を向上することができる。また形成される画像（すなわち更新される画像）は、ｋ空間の全領域のデータを用いているので、低周波領域のみならず、高周波領域についても画像の劣化がない。さらに低周波領域のデータについては常にその領域内の全サンプル点のデータを含んでいるので、コントラストの高い血流画像を得ることができる。
【００４３】
次にサンプル点の計測順序について説明する。ｋ空間を計測する順序としては、セントリックオーダー、シーケンシャルオーダー等が知られている。本発明の計測方法では、これらのいずれも採用できるが、図４に示す例では、サンプル点の原点からの距離によって順序を決めている。すなわち、本例では、まずｋ空間をｋｙ軸またはｋｚ軸に沿って２分割し、一方の領域(例えば、図中下側のＥ―Ｃ領域)では、サンプル点を原点からの距離が遠い順に計測順序をソートし、他方の領域(図中、上側の領域Ｃ−Ｅ領域)では、サンプル点を原点からの距離が近い順に計測順序をソートする。このように計測順序をソートした後、Ｅ−Ｃ領域の原点からの距離が最も遠いサンプル点から計測を開始し、その後はソートした計測順序に従って計測を進め、Ｃ−Ｅ領域の、原点からの距離が最も遠いサンプル点で計測を終了する。あるサンプル点を計測するということは、既に説明したように、そのサンプル点の座標（ｋｙ、ｋｚ）に対応するスライスエンコード傾斜磁場強度および位相エンコード傾斜磁場強度で、例えば図２のパルスシーケンスを実行することを意味する。
【００４４】
このような計測順序は、群が異なっていても同様とする。この計測順序によれば、計測開始から終了までの中央の時点を含む比較的広い時間範囲で低周波領域のデータを計測することになるので、目的とする血管の造影剤濃度が最も高い時点で、低周波領域を計測するようにタイミングを合わせることが容易となるという利点がある。これにより、目的血管についての画像を高コントラストで描画することができる。
【００４５】
但し、本発明の計測方法は、このような計測順序に限定されず、例えば図５に示すように、図中、左下端のサンプル点を開始点として、順番に右端まで進み、右端においてｋｚ軸を一つ上に進み、今度は左端に進むというようなシーケンシャルオーダーで行うことも可能である。
【００４６】
次に本発明の第２の実施形態として、周辺計測領域をさらに低中周波領域と高周波領域に分割し、それぞれについて領域内のサンプル点を空間的に均等な群に分ける場合を説明する。
【００４７】
本例では、図６に示すように、ｋ空間を原点を含む低周波領域の必須計測領域６１と、その外側の低中周波領域６２と、さらにその外側の高周波領域６３に分割する。各領域は互いに排他的な領域であって重ならない。必須計測領域６１は、第１の実施形態と同様に、各回の計測でデータ計測を行う領域であり、低中周波領域６２と高周波領域６３は、サンプル点が空間的に均等である複数の群に分割されており、各回の計測では群の一つずつを順次計測する。この場合、低中周波領域６２と高周波領域６３とでは、群の分割数を異ならせ、例えば低中周波領域６２は２つの群a、bからなり、高周波領域は４つの群a〜ｄからなる。これは、１回の計測で計測されるサンプル点の密度が、必須計測領域６１を１とするとき、低中周波領域６２ではその1/2、高周波領域６３では1/4になることと同義である。
【００４８】
この場合、１回目の計測では、必須計測領域６１と低中周波領域の群aと高周波領域の群aを計測し、次の計測では必須計測領域６１と低中周波領域の群bと高周波領域の群bを計測し、３回目の計測では必須計測領域６１と低中周波領域の群aと高周波領域の群cを計測する、というように、低中周波領域および高周波領域は、それぞれの群を順次、サイクリックに計測する。
【００４９】
高周波領域の計測が一巡する４回目で、３次元ｋ空間を埋めるデータが得られ、この３次元データを用いて画像を再構成する。以降は、各回の計測毎にその回で計測したデータと、その回で計測されず、過去の回で計測されたデータを用いて３次元データを作成し画像を再構成する。
【００５０】
この実施形態においても、サンプル点の計測順序は、図４に示したようにＥ−Ｃ領域およびＣ−Ｅ領域について、それぞれサンプル点を原点からの距離によって計測順序にソートし計測してもよいし、図５に示したようにシーケンシャルオーダーで計測してもよい。
【００５１】
第２の実施形態においても、第１の実施形態と同様に、１回の計測時間を短縮し、画像の時間分解能を向上することができ、その際、低周波領域のみならず高周波領域についても画像の劣化をなくすることができる。またこの実施形態では、周辺計測領域を２つに分割したことにより、目的とする血管の太さなどに応じて高コントラストで血管描出することができる。なお、周辺計測領域の分割数や、周辺計測領域をさらに分割した群のサンプル密度は、撮影の目的に併せて任意に変更することができる。
【００５２】
さらに第３の実施形態を説明する。この実施形態では、図７に示すように、ｋ空間を必須計測領域７１、周辺計測領域７２、非計測領域７３に３分割する。必須計測領域７１は、図３および図６に示す実施例と同様に、各回の計測で計測する領域である。周辺計測領域７２は、この必須計測領域７１の外側の領域で、第１の実施形態と同様に、サンプル点が空間的に均等である複数の群に分割されている。この分割数も第１の実施形態と同様に任意である。各計測において、群は一つずつ順次計測される。非計測領域７３は、計測を行わない領域である。
【００５３】
なお、図では必須計測領域７１を原点を中心とする円形(或いは楕円形)の領域とし、周辺計測領域をそれと同心円状の領域とする場合を示したが、これらは必須計測領域をｋ空間座標の辺と平行な或いは４５度傾斜した四角形の組み合わせや、四角形と円形の組み合わせとしてもよい。これは周辺計測領域を２つに分割した図６に示す第２の実施形態についても同様である。
【００５４】
この第３の実施形態についても、各回における必須計測領域７１および周辺計測領域７２の計測は、第１の実施形態と同じであるが、診断等に比較的重要でない３次元のｋ空間上の角の領域を非計測領域とすることにより、より短い時間で計測を繰り返せるようになるので、比較的小さな画質の劣化において時間分解能を向上することができる。
【００５５】
この場合にも画像再構成は、各回毎に行われ、その回で得られた必須計測領域のデータおよび周辺計測領域の一つの群のデータと、その回で得られず最近の過去の回で得られた他の群のデータとを用いると共に、非計測領域についてはゼロフィリングしたデータを用いて３次元データとし、これを３次元フーリエ変換することにより３次元画像データを得る。この３次元画像データから２次元投影画像を作成することは前述の二つの実施形態と同様である。
【００５６】
３次元画像データからの画像表示法および２次元投影画像の作成は、公知の手法を用いて行うことができる。例えば、再構成した３次元データと、造影剤投入前に計測しておいて３次元データとの差分をとり、これを最終的な３次元データとしてもよい。この差分は、複素差分であることが好ましいが、強度絶対値の差分であっても良い。また、今回のk空間データと造影剤投入前に計測しておいたk空間データとのk空間データ同士で複素差分を取った後に再構成を行うようにしてもよい。
【００５７】
さらに３次元データを、ＭＩＰ或いはＭｉｎＩＰ等の光軸軌跡法を用いて、冠状断、矢状断、軸横断等の任意の方向に投影し２次元画像を生成し、表示する。この場合、ある軸を中心として投影方向を回転させながら順次二次元画像を生成し、これらをつなげて動画像を生成するようにしてもよい。
【００５８】
さらにSurface rendering、Volume renderingなどの周知のレンダリング手法を用いてもよい。
以上、本発明の実施形態として、ky-kz方向について必須計測領域および周辺計測領域の境界を設定する場合について説明したが、これはkx,ky,kzの３方向のうちの任意の方向、または、任意の2つの方向の組み合わせ、または、３方向の全てについて、必須計測領域および周辺計測領域の境界を設定するようにしてもよい。
【００５９】
また、以上では、Z方向とY方向について位相エンコーディングを行うことにより３次元計測を行う場合について説明したが、Z方向についてはスライス選択傾斜磁場とRFパルスの組み合わせによって、その位置をエンコードする３次元計測を行う場合についても、以上の必須計測領域および周辺計測領域を用いた計測は同様に適用することができる。
【００６０】
さらに以上では３次元計測を行う場合について説明したが、２次元計測（1スライス毎の計測）を行う場合についてもkx-ky方向について必須計測領域および周辺計測領域の境界を設定することにより同様に適用することができる、
【００６１】
本発明のＭＲＩ装置およびそれを用いた計測方法は、造影剤と用いたダイナミックＭＲＡのみならず、ＭＲＩにおいて任意の繰り返し計測を行う場合について同様に適用することができる。
【００６２】
【発明の効果】
以上のように、本発明によれば、ＭＲＩにおいて繰り返し撮像を行う際に、高周波数領域の時間分解能を大きく劣化させることなく、各回の計測時間を短縮化する計測方法を提供することができる。また、このような計測方法において、さらに、計測結果の方向依存性を排除することができる。
【図面の簡単な説明】
【図１】本発明の実施形態に係る磁気共鳴イメージング装置の構成を示すブロック図。
【図２】本発明の実施形態に係る計測に用いることのできるパルスシーケンスの例を示すタイミングチャート。
【図３】本発明の第１の実施形態によるｋｙ−ｋｚ空間の分割を示す図。
【図４】本発明の第１の実施形態による計測法を示す図。
【図５】本発明の第1の実施形態による計測法を示す図。
【図６】本発明の第２の実施形態による計測法を示す図。
【図７】本発明の第３の実施形態による計測法を示す図。
【符号の説明】
1…被検体
2…磁場発生装置
3…傾斜磁場発生系
4…シーケンサ
5…送信系
6…受信系
7…信号処理系
8…CPU
9…傾斜磁場コイル
14a…送信側の高周波コイル
14b…受信側の高周波コイル
【図１】
0003847554.tif 000002
【図２】
0003847554.tif 000003
【図３】
0003847554.tif 000004
【図４】
0003847554.tif 000005
【図５】
0003847554.tif 000006
【図６】
0003847554.tif 000007
【図７】
0003847554.tif 000008