"静磁場に配置され所定のスピン密度を有する物体の画像化のため磁気共鳴データを収集する方法であって、a少なくとも一の高周波パルスと、空間的符号化磁気共鳴信号を生成するために所定の空間的次元内の少なくとも一の視野における前記磁気共鳴信号の位相を異なる符号化法で符号化する少なくとも二つの設定の空間的符号化位相コントラスト傾斜を含む磁場傾斜と、を前記物体に印加するステップと、b少なくとも二つの磁気共鳴信号をそれぞれ前記少なくとも二つの設定の異なる空間的符号化位相コントラスト傾斜の一つを用いて取得するステップと、c前記信号間の位相差を計算することにより、前記空間的次元に沿った前記物体の少なくとも一の平均スピン密度位置を決定するステップと、を含んでなることを特徴とする方法。請求項1に記載の方法であって、前記所定の空間的次元内の前記視野に沿った前記設定の二つの異なる空間的符号化位相コントラスト傾斜は、強度G=0であるゼロ傾斜、及び強度G>0である単極傾斜の設定、又は強度G>0かつ反対の極性である二つの単極傾斜の設定、又は強度G1>0,G2>0かつG1≠G2での二つの単極傾斜の設定、又は強度G=0であるゼロ傾斜、及び空間的符号化定常スピンに対しては0次オーダー傾斜モーメント>0であり一定の速度で移動するスピンに対しては1次オーダー傾斜モーメント=0である非対称双極傾斜の設定を含んでなることを特徴とする方法。請求項1又は2に記載の方法であって、前記物体は、物体層に沿って分布するスピン密度を有しており、ステップbは、物体空間において投影面に前記物体の投影を表す二つの二次元磁気共鳴画像を、前記投影面に直交する空間的次元に沿った前記異なる設定の空間的符号化位相コントラスト傾斜を用いて、取得するステップを含み、ステップcは前記物体の三次元画像を決定するステップを含んでなることを特徴とする方法。請求項3に記載の方法であって、前記異なる空間的符号化位相コントラスト傾斜は、前記高周波パルス及び磁場傾斜の反復間隔の直ぐ前後の連続する反復間隔毎に、又は前記高周波パルス及び磁場傾斜を用いた連続する画像取得毎に連続的に交互に設定されてなることを特徴とする方法。請求項1~4のいずれか一に記載の方法であって、前記高周波パルス及び磁場傾斜がデカルト空間的周波数位相符号化、放射状符号化、螺旋状符号化、又は他の非デカルト空間的符号化に適合されてなることを特徴とする方法。請求項1~5のいずれか一に記載の方法であって、前記高周波パルス及び磁場傾斜は、低フリップ角傾斜エコーMRIシーケンス、スピンエコーMRIシーケンス、励起エコーMRIシーケンス又はエコープラナー画像化シーケンスを含んでなることを特徴とする方法。請求項6に記載の方法であって、前記低フリップ角傾斜エコーMRIFLASHシーケンスは、アンダーサンプリング放射状FLASHシーケンスを含んでなることを特徴とする方法。請求項1又は2に記載の方法であって、前記物体は、物体線に沿って分布するスピン密度を有する、又は所定の範囲に広がる物体スポットに集中するスピン密度を有し、前記磁場傾斜は周波数符号化傾斜を含み、ステップaは、前記周波数符号化傾斜の方向に直交する二つの所定の空間的次元の二つの視野に沿った少なくとも三つの異なる設定の空間的符号化位相コントラスト傾斜を含み、ステップbは、少なくとも三つの周波数符号化磁気共鳴信号をそれぞれ前記異なる設定の空間的符号化位相コントラスト傾斜の一つを印加した後に取得するステップを含み、ステップcは前記物体の三次元画像を決定するステップを含んでなることを特徴とする方法。請求項8に記載の方法であって、前記設定の少なくとも三つの異なる空間的符号化位相コントラスト傾斜は、強度G=0である一のゼロ傾斜の設定と、強度G>0で前記二つの所定の空間的次元の前記二つの視野に沿って生成される二つの傾斜の設定と、を含んでなることを特徴とする方法。請求項1又は2に記載の方法であって、前記物体は、物体線に沿って分布するスピン密度を有する、又は所定の範囲に広がる物体スポットに集中するスピン密度を有し、前記磁場傾斜は周波数符号化傾斜を含み、ステップaは、前記周波数符号化傾斜の方向に直交する所定の空間的次元の視野に沿った二つの異なる設定の空間的符号化位相コントラスト傾斜を含み、ステップbは、二つの周波数符号化磁気共鳴信号をそれぞれ前記異なる設定の空間的符号化位相コントラスト傾斜の一つを印加した後に取得するステップを含み、ステップcは前記物体の二次元画像を決定するステップを含んでなることを特徴とする方法。請求項1又は2に記載の方法であって、前記物体は、所定の範囲に広がる物体スポットに集中するスピン密度を有し、ステップbは周波数符号化傾斜がない状態において四つの磁気共鳴信号を取得するステップを含み、前記設定の異なる空間的符号化位相コントラスト傾斜は、強度G=0である一のゼロ傾斜と、強度G>0で互いに直交する空間的次元に沿った前記三つの視野に沿って生成される三つの傾斜と、を含んでなることを特徴とする方法。請求項1~11のいずれか一に記載の方法であって、前記高周波パルスはそれぞれ、非選択的パルス、周波数選択的パルス及び空間選択的な高周波パルスのうちの少なくとも一つを含んでなることを特徴とする方法。請求項1~12のいずれか一に記載の方法であって、ステップcは、前記設定の異なる空間的符号化位相コントラスト傾斜によって定義される少なくとも一の視野における前記物体の厚さを前記磁気共鳴信号の強度から推定するステップを含んでなることを特徴とする方法。請求項1~13のいずれか一に記載の方法であって、前記ステップa~cは、前記物体の動的画像化を可能にするような時間分解能の繰返速度で繰り返されてなることを特徴とする方法。請求項1~14のいずれか一に記載の方法であって、前記設定の異なる空間的符号化位相コントラスト傾斜が、一以下の閉じたスピン密度を有する前記物体の所定軸に沿って方向付けられてなることを特徴とする方法。請求項1~15のいずれか一に記載の方法を実行するプログラムコードを含むコンピュータ可読媒体に記憶されてなることを特徴とするコンピュータープログラム。静磁場に配置され所定のスピン密度を有する物体の画像化のため磁気共鳴データを収集するよう構成される制御装置であって、少なくとも一の高周波パルスと、空間的符号化磁気共鳴信号を生成するために磁場傾斜と、を前記物体に印加すると共に、所定の空間的次元内の少なくとも一の視野における前記磁気共鳴信号の位相を異なる符号化法で符号化する少なくとも二つの設定の空間的符号化位相コントラスト傾斜のため構成される位相コントラスト傾斜設定装置を有する、パルス傾斜制御装置と、少なくとも二つの磁気共鳴信号をそれぞれ前記少なくとも二つの設定の異なる空間的符号化位相コントラスト傾斜の一つを用いて取得するよう構成される信号取得装置と、前記信号間の位相差を計算することにより、前記空間的次元に沿った前記物体の少なくとも一の平均スピン密度位置を決定するよう構成される計算装置と、を備えることを特徴とする制御装置。請求項17に記載の制御装置であって、前記位相コントラスト傾斜設定装置は、強度G=0であるゼロ傾斜、及び強度G>0である単極傾斜の設定、又は強度G>0かつ反対の極性である二つの単極傾斜の設定、又は強度G1>0,G2>0かつG1≠G2での二つの単極傾斜の設定、又は強度G=0であるゼロ傾斜、及び空間的符号化定常スピンに対しては0次オーダー傾斜モーメント>0であり一定の速度で移動するスピンに対しては1次オーダー傾斜モーメント=0である非対称双極傾斜の設定のために構成されてなることを特徴とする制御装置。請求項17又は18に記載の制御装置であって、物体層に沿って分布するスピン密度を有する前記物体を画像化する磁気共鳴データを収集するよう構成される制御装置であり、前記信号取得装置は、物体空間において投影面に前記物体の投影を表す二つの二次元磁気共鳴画像を、前記投影面に直交する空間的次元に沿った前記異なる設定の空間的符号化位相コントラスト傾斜を用いて、取得するよう構成されており、前記計算装置は前記物体の三次元画像を決定するよう構成されてなることを特徴とする制御装置。請求項19に記載の制御装置であって、前記パルス傾斜制御装置は異なる空間的符号化位相コントラスト傾斜を前記高周波パルス及び磁場傾斜の反復間隔の直ぐ前後の連続する反復間隔毎に、又は前記高周波パルス及び磁場傾斜を用いた連続する画像取得毎に連続的に交互に設定するよう構成されてなることを特徴とする制御装置。請求項17~20のいずれか一に記載の制御装置であって、前記パルス傾斜制御装置がデカルト空間的周波数位相符号化、放射状符号化、螺旋状符号化、又は他の非デカルト空間的符号化に適合されてなることを特徴とする制御装置。請求項17~21のいずれか一に記載の制御装置であって、前記パルス傾斜制御装置は、低フリップ角傾斜エコーMRIシーケンス、スピンエコーMRIシーケンス、励起エコーMRIシーケンス又はエコープラナー画像化シーケンスを生成するよう構成されてなることを特徴とする制御装置。請求項22に記載の制御装置であって、前記低フリップ角傾斜エコーMRIFLASHシーケンスは、アンダーサンプリング放射状FLASHシーケンスを含んでなることを特徴とする制御装置。請求項17又は18に記載の制御装置であって、物体線に沿って分布するスピン密度を有する、又は所定の範囲に広がる物体スポットに集中するスピン密度を有する前記物体を画像化する磁気共鳴データを収集するよう構成される制御装置であり、前記パルス傾斜制御装置は周波数符号化傾斜を生成するよう構成されており、前記位相コントラスト傾斜設定装置は、前記周波数符号化傾斜の方向に直交する二つの所定の空間的次元の二つの視野に沿った少なくとも三つの異なる設定の空間的符号化位相コントラスト傾斜のために構成されており、前記信号取得装置は、少なくとも三つの周波数符号化磁気共鳴信号をそれぞれ前記異なる設定の空間的符号化位相コントラスト傾斜の一つを印加した後に取得するよう構成されており、前記計算装置は前記物体の三次元画像を決定するよう構成されてなることを特徴とする制御装置。請求項24に記載の制御装置であって、前記位相コントラスト傾斜設定装置は、強度G=0である一のゼロ傾斜の設定と、強度G>0で前記二つの所定の空間的次元の前記二つの視野に沿って生成される二つの傾斜の設定と、のために構成されてなることを特徴とする制御装置。請求項17又は18に記載の制御装置であって、物体線に沿って分布するスピン密度を有する、又は所定の範囲に広がる物体スポットに集中するスピン密度を有する前記物体を画像化する磁気共鳴データを収集するよう構成される制御装置であり、前記パルス傾斜制御装置は周波数符号化傾斜を生成するよう構成されており、前記位相コントラスト傾斜設定装置は、前記周波数符号化傾斜の方向に直交する所定の空間的次元の視野に沿った二つの異なる設定の空間的符号化位相コントラスト傾斜のために構成されており、前記信号取得装置は、二つの周波数符号化磁気共鳴信号をそれぞれ前記異なる設定の空間的符号化位相コントラスト傾斜の一つを印加した後に取得するよう構成されており、前記計算装置は前記物体の二次元画像を決定するよう構成されてなることを特徴とする制御装置。請求項17又は18に記載の制御装置であって、所定の範囲に広がる物体スポットに集中するスピン密度を有する前記物体を画像化する磁気共鳴データを収集するよう構成される制御装置であり、前記信号取得装置は周波数符号化傾斜がない状態において四つの磁気共鳴信号を取得するよう構成されており、前記設定の異なる空間的符号化位相コントラスト傾斜は、強度G=0である一のゼロ傾斜と、強度G>0で互いに直交する空間的次元に沿った前記三つの視野に沿って生成される三つの傾斜と、を含んでなることを特徴とする制御装置。請求項17~27のいずれか一に記載の制御装置であって、前記パルス傾斜制御装置は、それぞれが非選択的パルス、周波数選択的パルス及び空間選択的な高周波パルスのうちの少なくとも一つを含んでなる高周波パルスを生成するよう構成されてなることを特徴とする制御装置。請求項17~28のいずれか一に記載の制御装置であって、前記計算装置は、前記設定の異なる空間的符号化位相コントラスト傾斜によって定義される少なくとも一の視野における物体の厚さを前記磁気共鳴信号の強度から推定するよう構成されてなることを特徴とする制御装置。請求項17~29のいずれか一に記載の制御装置であって、前記パルス傾斜制御装置及び前記信号取得装置は、前記物体の動的画像化を可能にするような時間分解能の繰返速度で繰り返し動作するよう構成されてなることを特徴とする制御装置。請求項17~30のいずれか一に記載の制御装置であって、前記位相コントラスト傾斜設定装置は、一以下の閉じたスピン密度を有する前記物体の所定軸に沿った前記異なる空間的符号化位相コントラスト傾斜を設定するよう構成されてなることを特徴とする制御装置。磁気共鳴走査器と、前記磁気共鳴走査器と接続される請求項17~31のいずれか一に記載の制御装置とを備えてなることを特徴とする磁気共鳴画像化MRI装置。リアルタイムのMR画像に対する基本的な改良は、FLASHシーケンスと、連続未処理データ取得の放射状ラジアルデータサンプリング及びビューシェアリングviewsharingとの組み合わせによってなされたS.チャン他、非特許文献5参照。ラジアルデータサンプリングによって、まずまずのアンダーサンプリングファクタ約2が可能になり、1フレーム当たり約250ミリ秒で画像の未処理データを取得できる。一部のフレームにのみ対応する現在の画像未処理データを前の画像未処理データと共に用いる画像更新の再構成つまり、いわゆるスライディングウィンドウ法によれば、約50ミリ秒もの時間分解能が得られるようになり、1秒当たり20のMR画像のビデオフレームレートが得られる。しかしながら、S.チャン他の方法非特許文献5では、画像当たりの全取得時間は変わらず250ミリ秒のままであるという事実から、画質に関する欠点が生じ得る。非線形逆再構成プロセスを、ラジアルFLASH取得の更なる高アンダーサンプリング画像未処理データに適用することによって、画像化プロセスを更に改良できる非特許文献6~9参照。さらに、速度符号化位相コントラストMRIは、MRIの初期段階に開発された流速の動的なマッピングに対してはよく構築された方法であった例えば非特許文献1~4参照。その最も一般的な形式において、速度符号化位相コントラストMRIは異なる双極速度符号化傾斜によって取得した二つの断面画像間の位相差に基づく。これらの速度符号化傾斜は、自己補償波形及び結果として得られる0次オーダー傾斜モーメント=0によって特徴付けられる。0次オーダー傾斜モーメント=0では、定常スピンは影響を受けないが、1次オーダー傾斜モーメント>0に対応する流れるスピンの位相は影響を受ける。0次オーダー傾斜モーメントは傾斜の時間積分である。一方、1次オーダー傾斜モーメントは傾斜に時間を乗じた時間積分である。一定の速度で流れるスピンについては、正味位相差は速度に正比例する。速度符号化位相コントラストMRI測定の結果を、二つの空間的次元と一つの速度位相差次元からなる3Dデータセットで表現でき、各画素は固有の速度値で表現される。本発明の目的は、MRデータを収集する改良方法、特に所定のスピン密度を有する物体のMRIのため、例えば医療用画像化のため、従来技術の欠点を回避可能な改良方法を提供することにある。特に、本目的は、非常に例えば100倍短い取得時間での物体の3DのMRIのため、MRデータを収集する方法を提供することである。さらに、これによって短い取得時間で3DのMR画像を連続的に取得可能となる。この結果、実用的なフレームレートで連続画像を、特に画質を損なうことなく提供できる。さらに本発明の他の目的は、MRデータを収集する改良制御装置、特に短い取得時間で物体の3DのMRIのための改良制御装置を提供することにあり、またMRI装置、特にMRデータを収集する方法及び物体の画像化を行うよう構成されたMRI装置を提供することにある。本発明の好ましい実施形態を、位相コントラスト傾斜の設計及び本発明の技術によって得られた画像の例を特に参照して、以下に説明する。なおMRI装置並びにその動作、傾斜エコーシーケンスの構成並びに特定の画像化される物体へのその応用、入手可能なソフトウェアツールを用いる数式の数値的実行、画像再構成、及び任意選択可能な追加の画像処理工程等、従来のMRI技術において既知の事項については、詳細な説明を省略する。以下ではデカルト座標系x,y,zで3D空間に配置された物体を例示する。位相コントラスト傾斜は、x軸、y軸及び/又はz軸の一つと平行に例示する。なお、デカルト座標軸に関する説明は例示のためにのみ導入したものであることを付言しておく。特に、z軸は、MRI装置における静磁場の方向と必ずしも一致しない。位相コントラスト傾斜の方向は、空間における物体形状及び物体の向きに応じて選択される。層状、線状又はスポット形状を有する物体を例示して説明するが、本発明を、他の物体、特に位相コントラスト傾斜の方向に沿った寸法が有限である物体にも同様に適用できる。2D構造の3DのMRIMRデータ収集の第一実施形態本発明の第一実施形態では、物体層O1の画像化を、図1A及び図1Bを参照して、以下に説明する。物体層O1は、x-y面と基本的に平行に延設される曲面に沿って分布するスピン密度を有する。物体層O1の最小寸法はz方向である。したがって、z方向が位相コントラスト傾斜Gpcを印加する所定の空間的次元として選択され、x-y面と平行な面が投影面1として選択されている。物体層O1のMRIでは、物体空間において投影面1上に物体層O1の投影を表す二つの2DのMR画像が、z方向に沿った異なる設定の空間的符号化位相コントラスト傾斜で収集される。二つの異なる位相コントラスト傾斜で取得されたMR信号の間の位相差は、投影面1に対する物体層O1の平均スピン密度位置を表す。z方向と平行な視野FOV1の最大寸法を、以下の関係に基づいて計算できる。2π=γ×傾斜強度×傾斜継続時間×FOV1γは利用される核例えばプロトンの磁気回転比であり、空間的符号化位相コントラスト傾斜が矩形波である理論的な場合には、0次傾斜モーメントをその強度と継続時間との積から計算できる。実際上、この積は、傾斜波形の時間積分で置き換えられる速度符号化傾斜では常にゼロとなる。図1Bは、先に非特許文献5を参照して説明した放射状FLASHMRIシーケンスを利用した、物体層O1をマッピングするのためのRFパルス磁場傾斜シーケンスの実際の実施を示す。物体層O1に第一非選択的RF励起パルスを印加した後に、第一単極位相コントラスト傾斜Gpcを例えばz方向の正方向に印加する。次に、物体層O1に周波数符号化傾斜Gx及びGyを印加し、MR信号、例えば図の傾斜エコーが取得される。この第一反復間隔の直ぐ後に第二反復間隔が続く。再び、非選択的RF励起後、z方向の負方向の第二単極線位相コントラスト傾斜Gpcが印加される。第二MR信号は、第一反復間隔と同じ周波数符号化傾斜Gx及びGyを用いて取得される。投影面1上の物体層O1の二つの完全な投影画像のデータが収集されるまで、異なる位相コントラスト傾斜を用いた二つの反復間隔のシーケンスは、異なる設定の周波数符号化傾斜Gx及びGyを用いてn回繰り返される。投影面1上の各画素に対して、z方向における物体層O1の位置つまりそれぞれの位相差及びMR信号の関連する強度相対振幅が得られ、これにより、物体層O1の完全な3D画像を計算できる。所望により、非選択的RF励起を、空間選択的な、又は化学シフト選択的なRF励起で置き換えることもできる。さらに、図1Bに示す態様においては、典型的には小さい位相コントラスト傾斜を、間挿法で、つまり連続する反復間隔毎に交互に、印加する。しかしながら、異なる空間的位相コントラスト傾斜を用いて完全な画像を連続的に取得することも可能である。この態様は動的な応用では有用となる。フレームレートの二倍で位相差マップのスライディングウィンドウを再構成できるからである。また、従来のフーリエ法とは対照的に、離散的サンプリングパターンによって影響を受けることなく、又は点拡がり関数によって悪影響を受けることなく、位相差次元に沿った物体の位置を数値によって正確に定義されることを付言しておく。3D空間における曲面状構造を高速で動的に視覚化するために、空間的符号化位相コントラスト画像の放射状符号化を、アンダーサンプリングを大きくすることによって、実現することもできる。結果として、画像再構成を、放射状符号[非特許文献7を用いてパラレルMRIを改良するよう最近開発された正規化非線形インバージョンアルゴリズム非特許文献6を用いて、極めて高い時間分解能で達成できる非特許文献8及び9。これらの再構成は、反復的アルゴリズムの計算要求が高く、オフラインで行われるので、非特許文献5に記載されているように相補的放射状符号化を用いて複数の典型的に五つの連続的に取得されたデータセットを組み合わせることにより取得可能な、より完全なデータセットのグリッディングによって、マグニチュードイメージをリアルタイムにオンラインでモニタリングできる。空間的符号化位相コントラストMRIに対して非線形インバース再構成アルゴリズムを用いることにより、2Dコイル感度マップのみを用いて、3D空間におけるパラレルイメージングが正しく得られる。有効感度は、2D構造に沿った物理的3Dコイル感度の選択された画像化面投影面1上への射影である。3D空間において移動する面状構造に対しては、この射影は時間と共に急激に変化するので、選択された非線形インバースアルゴリズム非特許文献6~8によって保証されるようなコイル感度の連続的評価を必要とする。コイル感度の変化による位相誤差を回避するために、各再構成された画像にそれと同時に評価されたコイル感度をまず乗じて、チャネル毎に別々の画像を得る。それぞれの位相差を計算するとき、非共鳴効果及びコイル感度からの背景位相を除去する。最後に、すべてのチャネルからの位相コントラストマップは、各画素をそれぞれのチャネルのマグニチュードイメージにおけるエネルギーで重み付けして、組み合わせる。1D構造及び点源の2D/3DのMRIMRデータ収集の第二実施形態空間的符号化位相コントラストMRIの他の形態では、合計二又は三の反復間隔のみに対応する更に短い測定時間における1D物体つまり線状構造の2D及び3Dマッピングを説明する。これは、2D試験では、異なる位相符号化傾斜での単一周波数符号化MR信号の二回の取得だけしか必要としないからである。このようなデータは、各周波数符号化位置に対して位相差値を得るのに十分である。その結果としての2Dマップは一の周波数及び一の位相差次元を表す。その再構成では、二回の1Dフーリエ変換及びそれぞれの位相差の一回の計算を必要とするだけである。同様に、図2A~図2Cに示すように、3次元的曲線状1D物体の一の周波数符号化次元及び二つの空間的位相差次元を用いる3Dマッピングを、僅か三つの反復間隔の測定時間内に実行できる。図2Aは、3次元的曲線状1D物体物体線O2を3Dマッピングする本発明のMRデータ収集方法の応用を示す。図2Aでは、物体線O2は、y軸と平行なある決まった主軸を有する曲線状糸状体である。この主軸は、周波数符号化傾斜Gyが印加される投影軸2として選択されている。位相コントラスト傾斜Gpc1及びGpc2を印加するために選択された空間的次元は、投影軸2に直交するよう、つまりx方向及びz方向に方向付けられる。第一位相コントラスト傾斜Gpc1によって、第一FOV1がz方向に沿って広がり、第二位相コントラスト傾斜Gpc2によって、第二FOV2がx方向に沿って広がる。図2Bでは、物体線O2の3D画像化に対して、単一周波数符号化MR信号の三回の取得に、例えば二つの垂直単極位相符号化傾斜と空間的位相符号化を行わない一つの基準との取得が必要とされる。各周波数符号化位置に対して、二組のMR信号が得られる。一組は、二つの相符号化傾斜を用いてそれぞれ収集され、他方の組は空間的位相符号化を行わず収集される。各組の信号間の位相差は、それぞれの方向における平均スピン密度位置を、例えばz方向及びx方向におけるz1及びx1を与える。また他に、図2Cに示すように、空間的位相符号化の組み合わせることも可能である。空間的位相コントラストMRIにおける1D物体の通常の寸法厚さに関する情報は、個々の投影の振幅情報から得られる。同じことが2D物体の3DのMRI及び対応するマグニチュードイメージの振幅にも言える。この情報は緩和強調プロトン密度における差異を含んでいるが、物体内の相対的差異を保持する3D再構成に体素サイズを動的に割り当てるよう用いることもできる。点源の2D/3DのMRIMRデータ収集の第三実施形態本発明の第三実施形態では、点源の3Dマッピングを四つの励起及び更に短いTR間隔を用いて行うことができる。周波数符号化傾斜を用いることなく、それぞれの3Dマッピングを、RF励起パルスを用いるが傾斜を全く用いない一の反復間隔と、RF励起パルス及び空間的符号化位相コントラスト傾斜を用いる三つの反復間隔例えば三つの直交する方向と、からのデータ収集に基づいて、行うことができる。各反復間隔は、少なくとも一の複合データ点の取得を含んでいる。これを図3A及び図3Bに示す。ドット状の物体物体スポットO3の平均スピン密度位置を得るために、それぞれz方向、x方向及びy方向に方向付けられた三つの位相コントラスト傾斜Gpc1、Gpc2及びGpc3が印加される。したがって、3D空間は、直交するFOV1、FOV2及びFOV3に広がる。三組のMR信号が得られる。一組は、三つの相符号化傾斜を用いてそれぞれ収集され、他方の組は空間的位相符号化を行わず収集される。各組の信号間の位相差は、それぞれz方向、x方向及びy方向における平均スピン密度位置を与える。本発明の制御装置を有するMRI装置図4に、本発明の制御装置100とMRスキャナ200とを有するMRI装置300の実施形態を概略的に示す。制御装置100は、位相コントラスト傾斜設定装置111及びRF制御ユニット112を有するRFパルス傾斜制御装置110と、信号取得装置120と、計算装置130と、を有する。MRスキャナ200は、主磁界装置210と、磁気傾斜装置220と、励起/取得コイル装置230と、を有する。さらに、保持装置240が調査される物体O1を支持するよう配置される。構成要素210~240は、従来のMRスキャナにおいて既知の通り、構成できる。傾斜制御装置110の位相コントラスト傾斜設定装置111は、位相コントラスト傾斜を設定するよう構成された回路によって、又は磁気傾斜装置220の設定を制御するのに適したソフトウェアを実行するコンピュータユニットによって、実装することができる。RFパルス傾斜制御装置110は、励起/取得コイルデバイス230及び磁気傾斜装置220と接続され、上述の通り、励起RFパルス及び磁場傾斜を生成する。さらに、励起/取得コイルデバイス230は、信号取得装置120と接続され、物体O1において生成されたMR信号を検知する。信号取得装置120は、物体の少なくとも一の平均スピン密度位置を決定する又は上述した方法のうちの一つを用いて物体の完全な画像を再構成するのに適した計算装置130と接続される。模型及び人間のMRI研究本発明の方法の実験的検証を、市販のMRスキャナTIMトリオ、ジーメンスヘルスケア、エルランゲン、ドイツ及び標準的な32要素ヘッドコイルを用いて、3テスラの静磁場で行った。図1Bに示すように、空間的符号化位相コントラストMRIの好ましい実施では、空間的位相符号化するため空間非選択的RF励起パルス及び間挿法一反復間隔おき又は連続法一画像おきでRFスポイリング放射状FLASHシーケンスを用いた。不純物を添加した水を注入した管の異なる模型短縮されたT1緩和時間によって、空間的符号化位相コントラストMRIの原理を実証した。それぞれの画像を、全周波k空間サンプリングで間挿符号化法を用いて取得し、密度補償と、カイザーベッセル核を用いる畳み込みと、規則的グリッドに対する補間と、逆高速フーリエ変換と、k空間における補間を補償するロールオフ補正と、を含む従来のグリッディングアルゴリズムによって再構成した非特許文献5参照。図5は、全周波サンプリング放射状FLASHシーケンス及びグリッディング再構成を用いて得られた、3D螺旋として配置された水注入管の空間的符号化位相コントラストMR画像の例を示している。図5の上側は3D表現を示し、下側は対応する2D位相コントラストマップグレースケールで符号化された位相差のオーバーレイを用いたマグニチュードイメージを示す。画像は、非選択的RF励起パルスフリップ角12°を用いて1.0mmの線形分解能反復時間TR=2.35ミリ秒、エコー時間TE=1.61ミリ秒で得た。さらに三人の健常者が、2D構造の動的3DのMRIに対応する手と指の動きに関する研究に参加した。人体研究は、心臓血管のリアルタイムMRIのために既に開発されている二重エコー放射状FLASHMRIバージョンを用いて行った非特許文献9参照。適切に応用することにより、逆位相及び同位相条件で二つの2D位相コントラストマップを同時に記録でき、水/脂肪信号に貢献できる。これらの相補的位相情報をそれぞれの3D表現を改善するために利用した以下参照。高アンダーサンプリング放射状FLASH画像に対する取得パラメータTR=2.86ミリ秒、TE1/TE2=1.25/2.25ミリ秒、フリップ角3°、2.0mmの面分解能により、一画像当たり25ミリ秒9つのスポーク又は位相差マップ当たり50ミリ秒つまり2×9スポークの最小測定時間でリアルタイム3DのMRI動画を達成した。これらの応用では、異極性で二つの位相符号化画像の連続取得法を用いた。すべてのアンダーサンプリング画像は、正規化非線形インバージョン非特許文献6~8参照により、再構成した。位相コントラスト情報の空間的位置及び多次元視覚化への変換は、MATLABマスワークス、ナティック、マサチューセッツ州で作成されたルーチンを用いて行った。空間的位相情報が単一の数値として与えられるので、その3D視覚化が体素として、それぞれの大きさマグニチュード情報の信号振幅に対応する適当な厚さを通じて位置に任意の数を代入することによって、実現される。最小強度閾値の選択と共に、工程はノイズの影響を効果的に除去し、調査下の物体内の概略的な構造差を維持する。図6は、高アンダーサンプリング二重エコー放射状FLASHで記録された人間の手の空間的符号化位相コントラストMRI研究の基礎的な結果を示す。グレー値で符号化された位相オーバーレイを用いたマグニチュードイメージから、インビボ状態における位相情報は、逆位相画像及び同位相画像図6A及び図6Bの両方において、磁界不均質性及び水/脂肪化学シフト差によって特徴付けられることが実証される。これらの効果は対応する位相コントラストマップ図6Cにおいては打ち消し合うが、空間的符号化による位相差が保持される。また、逆位相状態により、例えば骨髄に強度がより低い領域が現れ図6上側、一方、同位相状態により、より連続した信号を得るのに役立ち、可能な限り多くの画素において信頼できる位相コントラスト情報が得られる図6下側ことが分かる。しかしながら、両方のマップから有用な相補的情報が得られるので、選択された視覚化法では、両方の位相コントラストマップにおいて一画素当たりの位相値と同位相マグニチュードイメージを組み合わせることにより、すべての情報を利用した。そして、結果として得られるマップを3D表現に変換した。図7に、空間的符号化位相コントラストMRI動画上部の連続する3D表現と、50ミリ秒の時間分解能で動いている人間の手の元の2D位相コントラストマップ底面図6と同じデータ及びパラメータと、を示す。動きを見易くするために、画像上側の左端から下側の右端までは、長い動画のうちの時間800ミリ秒の範囲の3番目のフレーム150ミリ秒ごとだけを表している。そして、実験的検証を、5.0ミリ秒の時間分解能での1D物体物体線の動的3DのMRIにも行った。後者の実験は三つの異なる空間的位相符号化図2を参照を用いて単一周波数符号化MR信号投影により行うので、それぞれの3D再構成では3回のフーリエ変換及び二回の対応する位相差計算しか必要としない。適切な計算は非常に高速であり、動的な線状物体のリアルタイムの3D視覚化には既存のMRIシステムで直ちに実行できる。図8に、1.0mmの線形分解能TR/TE=2.00/1.32ミリ秒で非選択的RF励起フリップ角12°を用いた周波数符号化傾斜エコーの二回又は三回の取得のみで収集された、水が注入された二つの異なる曲線状の管の、二つの空間的符号化位相コントラストMR画像を示す。図8Aに示す2D面内に延びる曲線状管の2D表現は、4.0ミリ秒の合計測定時間内の異なる位相符号化MR信号の二回の取得に基づく。図8Bは、合計測定時間6.0ミリ秒内のそれぞれ二つの直交する位相符号化傾斜を用いたMR信号の取得と空間的位相符号化を行わないMR信号の取得との三回の取得に基づく、3D体積内に延びる曲線状管の3D表現を示す。ここでは、グレー値が空間的な深さをよく視覚化している。捩れたゴムバンドに上端部が固定された状態で手を放した後の回転するバナナの選択した複数の図を表す、速く動いている線状物体の動的研究を図9に示す。データは、非選択的RF励起パルスフリップ角6°を用いて2.0mmの線形分解能TR/TE=1.68/1.15ミリ秒で収集した。対応する3DMRI動画は、5.0ミリ秒の時間分解能、つまり1秒当たり200フレームを達成した。画像上側の左端から下側の右端までは、長い動画のうちの時間505ミリ秒の範囲の20番目のフレーム毎100ミリ秒ごとを表している。上記説明、図面及び特許請求の範囲において開示した本発明の特徴は、本発明を異なる態様で個別に又は組み合わせて実現するものにも同様に有効である。"