磁気共鳴イメージング装置及び磁気共鳴イメージング方法
【課題】流体の流速を求めること。
【解決手段】実施形態に係る磁気共鳴イメージング装置は、収集部と、特定部と、取得部と、算出部とを備える。前記収集部は、被検体内を移動する流体の画像である流体画像を複数収集する。前記特定部は、前記複数の流体画像のうちの１つである基準画像と各流体画像との差分画像を用いて、前記流体の移動距離を特定する。前記取得部は、前記移動距離に対応する経過時間を、前記複数の流体画像の収集に用いられたパルスシーケンス情報から取得する。前記算出部は、前記移動距離を前記経過時間で除算することで前記流体の流速を算出する。
【解決手段】実施形態に係る磁気共鳴イメージング装置は、収集部と、特定部と、取得部と、算出部とを備える。前記収集部は、被検体内を移動する流体の画像である流体画像を複数収集する。前記特定部は、前記複数の流体画像のうちの１つである基準画像と各流体画像との差分画像を用いて、前記流体の移動距離を特定する。前記取得部は、前記移動距離に対応する経過時間を、前記複数の流体画像の収集に用いられたパルスシーケンス情報から取得する。前記算出部は、前記移動距離を前記経過時間で除算することで前記流体の流速を算出する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、磁気共鳴イメージング装置及び磁気共鳴イメージング方法に関する。
【背景技術】
【０００２】
従来、磁気共鳴イメージング装置（以下、ＭＲＩ（Magnetic Resonance Imaging）システム）による撮像法のひとつに、造影剤を用いずに被検体内を移動する流体を撮像する手法がある。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００８−２３３１７号公報
【特許文献２】米国特許出願公開第２００８／０００９７０５号明細書
【非特許文献】
【０００４】
【非特許文献１】Miyazaki, et al., “Non-contrast-enhanced MR angiography using 3D ECG-synchronized half-Fourier fast spin echo,” JMRI 12:776-783 (2000)
【非特許文献２】Suga, et al., Lung perfusion impairments in pulmonary embolic and airway obstruction with non contrast MR imaging,” J Appl Physiol 92:2439-2451 (2002)
【非特許文献３】Takahashi, et al., “Non-contrast-enhanced renal MRA using time-spatial labeling pulse(t-SLIP) with 3D balanced SSFP,” presented at the ISMRM 15th Annual Meeting, Berlin, Germany, page 179(2007)
【非特許文献４】Yamamoto, et al., “Non-contrast-enhanced MRDSM of peripheral arteries using continuous acquisitions of ECG-triggered 2D half-Fourier FSE within a cardiac cycle,” 12th Annual Meeting, Toronto, Canada, page 1709(2003)
【非特許文献５】Kanazawa, et al., “Time-spatial labeling inversion tag (t-SLIT) using a selective IR-tag on/off pulse in 2D and 3D half-Fourier FSE as arterial spin labeling,” presented at the ISMRM 10th Annual Meeting, Hawaii, page 140(2002)
【非特許文献６】Furudate, et al., “FBI-Navi for Easy Determination of Diastolic and Systolic Triggering Phase in Non-Contrast Fresh Blood Imaging(FBI),” ISMRM 16th Annual Meeting, Toronto, Canada, page 2902(2008)
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明が解決しようとする課題は、流体の流速を求めることである。
【課題を解決するための手段】
【０００６】
実施形態に係る磁気共鳴イメージング装置は、収集部と、特定部と、取得部と、算出部とを備える。前記収集部は、被検体内を移動する流体の画像である流体画像を複数収集する。前記特定部は、前記複数の流体画像のうちの１つである基準画像と各流体画像との差分画像を用いて、前記流体の移動距離を特定する。前記取得部は、前記移動距離に対応する経過時間を、前記複数の流体画像の収集に用いられたパルスシーケンス情報から取得する。前記算出部は、前記移動距離を前記経過時間で除算することで前記流体の流速を算出する。
【図面の簡単な説明】
【０００７】
【図１】図１は、実施形態に係るＭＲＩシステムの概略ブロック図である。
【図２】図２は、実施形態における心収縮期及び心拡張期に対する連続的な心電同期ＭＲスライス撮像シーケンスを示すタイミング図である。
【図３】図３は、図２に示すように、連続する心臓ゲーティング時間増分で「暗い（信号の輝度が低い）」心収縮期画像と「明るい（信号の輝度が高い）」心拡張期画像との間の差分により得られる連続する差分画像を示す概略図である。
【図４】図４は、図３の概略図と類似の概略図であるが、実施形態によって、標準又は平均の血流速度を計算できる方法を説明する注釈を含む。
【図５】図５は、図３〜図４とは異なる画像を示す類似の概略図であるが、実施形態によって、特定の速度を各連続する区間で計算できる方法を説明する注釈を含む。
【図６】図６は、図５と類似であるが、実施形態の一連の異なる画像の全体にわたって、全標準又は平均の速度を計算できる方法を具体的に示している。
【図７】図７は、Ｔｉｍｅ−ＳＬＩＰ（Time-Spatial Labeling Inversion Pulse）撮像法（流出及びタグ付け／タグ付けなし減算法）を使う別の実施形態による速度測定の概略を示す。
【図８】図８は、実施形態で利用可能なコンピュータプログラムコード構造の一例に関する概略フローチャートである。
【図９】図９は、２つの連続する撮像時間の間の蛇行状の経路に沿った血流に対する連続的な経路位置ポイントを操作者が定義する際に利用できる操作者画面の概略図である。
【図１０】図１０は、実施形態における時間及びスライス位置の両方又は一方に対する特定の速度及び平均速度の両方又は一方として考えられる出力表示グラフの概略図である。
【発明を実施するための形態】
【０００８】
図１に示すＭＲＩ（Magnetic Resonance Imaging）システムは、架台部１０（断面図にて示す）と、互いに接続される様々な関連のシステム構成要素２０とを含む。少なくとも架台部１０は、通常シールドルーム内に設置される。図１に示す１つのＭＲＩシステムは、静磁場Ｂ0磁石１２と、Ｇx、Ｇy、Ｇz傾斜磁場コイルセット１４と、ＲＦ（Radio Frequency）コイルアセンブリ１６との実質的に同軸円筒状の配置を含む。この円筒状に配置された要素の水平軸線に沿って、被検体テーブル１１によって支持された被検体９の頭部を取り囲むように示された撮像ボリューム１８がある。
【０００９】
ＭＲＩシステム制御部２２は、表示部２４、キーボード／マウス２６、及びプリンタ２８に接続された入力／出力ポートを備える。言うまでもなく、表示部２４は、制御入力もまた備えるような多様性のあるタッチスクリーンであってもよい。
【００１０】
ＭＲＩシステム制御部２２は、ＭＲＩシーケンス制御部３０とインタフェース接続する。ＭＲＩシーケンス制御部３０は、Ｇx、Ｇy、Ｇz傾斜磁場コイルドライバ３２、並びに、ＲＦ送信部３４及び送信／受信スイッチ３６（同じＲＦコイルが送信及び受信の両方に使用されている場合）を順に制御する。ＭＲＩシーケンス制御部３０は、既にＭＲＩシーケンス制御部３０にて可能な、ＭＲＩデータ取得シーケンス（血液等の流動する流体（flowing fluid）の撮像を含む）の実行に適切なプログラムコード構造３８を含む。心臓信号取得装置８（例えば、被検体９の身体に適切に装着された、ＥＣＧ（electrocardiogram）又は末梢血管パルストランスデューサ（単数又は複数））は、ＭＲＩシーケンス制御部３０のトリガー信号とする心電ゲート信号１３を出力することができる。
【００１１】
ＭＲＩシステムは、表示部２４に出力する処理された画像データを作成できるように、ＭＲＩデータ処理部４２に入力を供給するＲＦ受信部４０を含む。また、ＭＲＩデータ処理部４２を、画像再構成プログラムコード構造４４及びＭＲ画像記憶部４６にアクセスできるように構成してもよい（例えば、実施形態及び画像再構成プログラムコード構造４４に従った処理から得られたＭＲＩデータを格納するために）。
【００１２】
また、図１は、ＭＲＩシステムプログラム／データ格納部５０を一般化した描写を示す。ＭＲＩシステムプログラム／データ格納部５０に格納されるプログラムコード構造（例えば、画像再構成のため、及び、流速、距離／時間の測定のため、操作者入力等のための）は、ＭＲＩシステムの様々なデータ処理構成要素にアクセス可能なコンピュータ読み取り可能な記憶媒体に格納される。当業者には自明であるが、ＭＲＩシステムプログラム／データ格納部５０は、正常運転時にそのように格納されたプログラムコード構造に対して直近の必要性を有するＭＲＩシステムの処理コンピュータのうちの様々なコンピュータに分割し、且つ少なくとも一部を直結してもよい（すなわち、ＭＲＩシステム制御部２２に普通に格納したり、直結したりする代わりに）。
【００１３】
実際、当業者には自明であるが、図１の描写は、本明細書で後述する実施形態を実行できるように若干の変更を加えた一般的なＭＲＩシステムを非常に高度に簡素化した図である。システム構成要素は、様々な論理収集の「ボックス」に分割でき、通常、多数のデジタル信号処理装置（ＤＳＰ（Digital Signal Processors））、超小型演算処理装置、特殊用途向け処理回路（例えば、高速Ａ／Ｄ変換、高速フーリエ変換、アレイ処理用等）を含む。これら処理装置のそれぞれは、通常、各クロックサイクル（又は所定数のクロックサイクル）が発生すると、物理データ処理回路がある物理的状態から別の物理的状態へ進むクロック動作型の「状態機械」である。
【００１４】
動作中に、処理回路（例えば、ＣＰＵ（Central Processing Unit）、レジスタ、バッファ、計算ユニット等）の物理的状態が、あるクロックサイクルから別のクロックサイクルへ漸進的に変化するだけでなく、連結されているデータ格納媒体（例えば、磁気記憶媒体のビット格納部）の物理的状態も、そのようなシステムの動作中に、ある状態から別の状態へ変わる。例えば、ＭＲＩ再構成プロセスの終了時、物理的記憶媒体のコンピュータ読み取り可能なアクセス可能データ値格納場所のアレイは、いくつかの事前の状態（例えば、全部一律の「ゼロ」値又は全部「１」値）から新しい状態に変わる。その新しい状態では、そのようなアレイの物理的場所の物理的状態は、最小値と最大値との間で変動し、現実世界の物理的事象及び状況（例えば、撮像ボリューム空間内の被検体の動脈内を流れる血液）を表現する。当業者には自明であるが、格納されたデータ値のそのようなアレイは、物理的構造を表し且つ構成もする。つまり、命令レジスタの中に順次読み込まれてＭＲＩシステムの１つ以上のＣＰＵによって実行されたとき、動作状態の特定シーケンスが発生して、ＭＲＩシステム内中に移行されるコンピュータ制御プログラムコードの特定構造が構成される。
【００１５】
下記の実施形態は、流動する流体の速度（以下、適宜「流速」という）を含む、データ取得の処理及びＭＲ画像の生成表示の両方又は一方を行うための改良された方法を提供する。
【００１６】
心収縮期及び心拡張期のタイミングと同期して収集された画像間には、ＭＲ信号強度の変化がある。このＭＲ信号強度の変化は、例えば、信号値の高い心拡張期の画像から信号値の低い心収縮期の画像を減算すること、及び、信号値の低い心収縮期の画像から信号値の高い心拡張期の画像を減算することの両方あるいは一方を行うことによって、時間分解された流体（例えば、血液）の流れを非造影にて得るために用いることができる。これは、時間分解非造影ＭＲＡ（磁気共鳴血管造影法（Magnetic Resonance Angiography））と呼ばれることもある。また、流体の流れのための他の非造影撮像技術として、ＡＳＬ（Arterial Spin Labeling）と呼ばれることもあるＴｉｍｅ−ＳＬＩＰ（Time−Spatial Labeling Inversion Pulse）がある。
【００１７】
これらの非造影撮像技術においては、流体の流速測定を容易にするための操作者インタフェースが、望まれるであろう。
【００１８】
以下に示す実施形態によれば、非造影の流体画像における流体の流速測定は、少なくとも２つの手法で実現することができる。第１は、心拡張期と心収縮期との減算画像を用いる手法である。ＭＲＩシステムは、心拡張期と心収縮期との減算画像を用いて、流体の移動距離を特定することができる。例えば、動脈経路が比較的線形である場合、開始ポイントから最終ポイントまでの経路を単純な線であると考えることができる。あるいは、動脈経路が蛇行状の場合でも、経路に沿って効果的に置かれたポイントをつなぐことによって、開始ポイントから最終ポイントまでの経路を特定することができる。任意に形成された経路に沿って流動中の流体が移動した距離を自動的に追跡することも有効であるといえる。とにかく、流動中の流体が連続する撮像時間の間に通過した距離は、連続する撮像時間の間の特定の速度、並びに、一連のそのような撮像時間にわたる標準又は平均の速度を算出するために使用することができる。
【００１９】
第２は、非造影Ｔｉｍｅ−ＳＬＩＰ撮像技術を使用する手法である。非造影Ｔｉｍｅ−ＳＬＩＰ撮像技術を使用する場合、ｋ空間に対する位相エンコードの配列がセントリックオーダリングかシーケンシャルオーダリングかによって、ＦＡＳＥ（Fast Asymmetric Spin Echo）又はｂＳＳＦＰ（balanced Steady State Free Precession）などの若干異なる撮像シーケンスを含んでもよい。もちろん、流速の計算において分子は、やはり、信号を発生して流れる物質（例えば、血液）が移動した距離である。流速の計算において分母として使用される時間増分は、ＦＡＳＥ撮像法を用いたＴｉｍｅ−ＳＬＩＰの場合、ＢＢＴＩ（Black-Blood Time to Inversion）時間、プラス、実効エコー時間（ＴＥｅｆｆ）である。また、セントリックオーダリングでｂＳＳＦＰを用いたＴｉｍｅ−ＳＬＩＰの場合、ＢＢＴＩである。また、シーケンシャルオーダリングでｂＳＳＦＰを用いたＴｉｍｅ−ＳＬＩＰの場合、ＢＢＴＩ時間、プラス、１／２ＥＴＬ（エコートレイン長）（位相エンコード数の半数分に相当する時間）である。
【００２０】
速度測定は、異なるＥＣＧ（Electrocardiogram）（心電図ゲーティング）信号ポイント（例えば、心収縮期から心拡張期）間に流量変位を生成する、又は、ＦＢＩ（Fresh Blood Imaging）−Ｎａｖｉ信号差プロットを使用するＭＲ信号に対して行うことができる。流動するＭＲ信号変位又は移動したＭＲ信号の速度測定を計測できる。
【００２１】
上述のように、（ｋ空間配列に依存してＴＥｅｆｆ及びＥＴＬ／２の両方又は一方をおそらくプラスする）ＢＢＴＩ時間の間にＴｉｍｅ−ＳＬＩＰ撮像技術を使って流量を生成するＭＲ信号の変位もまた、平均速度の計算に使用できる。
【００２２】
一実施形態では、移動距離に沿った流速を選択的に計算するために（例えば、マウスボタンの押下によって）、簡単なＧＵＩ（Graphical User Interface）を使用できる。移動距離は、ユーザーによって（例えば、経路に沿って効果的に置かれたポイントを画定することによって）、又は、システム提供の自動追跡機能によって、輪郭を描出できる。とにかく、移動距離は、次に、（例えば、ＥＣＧ遅延撮像技術及びＴｉｍｅ−ＳＬＩＰ撮像技術の両方又は一方において減算された画像間の）その距離に関連する時間間隔によって除算される。
【００２３】
ＧＵＩは、平均速度及び特定速度の両方又は一方の流量移動の測定に使用できる。非造影技術では、速度測定は、少なくとも２つの方法で測定できる。第１の方法は、心拡張期−心収縮期減算撮像技術を使って、ＭＲ信号変位（最も近いポイント又は長く連結された線による）の距離に沿って変位線を描くことができる。ＭＲＩ表示システムは、自動的にこれを距離として記録して、次に、適切なボタンのクリック（又は、例えば、速度選択のための正しいマウスクリック）によって、（例えば、減算される心収縮期及び心拡張期の画像の有効経過時間の間の）適切な対応する時間差によってそのような距離を除算する計算を実現する。第２の方法（例えば、ＡＳＬ（arterial spin labeling）法）では、ＭＲＩ信号は、ＦＡＳＥを用いたＴｉｍｅ−ＳＬＩＰに対する（ＢＢＴＩ＋ＴＥｅｆｆ）、及びセントリックオーダリングでｂＳＳＦＰを用いたＴｉｍｅ−ＳＬＩＰに対する（ＢＢＴＩ）、及び（シーケンシャルオーダリングでｂＳＳＦＰを用いたＴｉｍｅ−ＳＬＩＰに対するＢＢＴＩ＋ＥＴＬ／２）によって除算された距離内で、最初の信号から最後の信号まで移動する流量を生成する。
【００２４】
さて、上述してきた実施形態に係るＭＲＩシステムを改めて説明すると、実施形態に係るＭＲＩシステムは、収集部と、特定部と、取得部と、算出部とを備える。収集部は、被検体内を移動する流体の画像である流体画像を複数収集する。特定部は、複数の流体画像のうちの１つである基準画像と各流体画像との差分画像を用いて、流体の移動距離を特定する。取得部は、移動距離に対応する経過時間を、複数の流体画像の収集に用いられたパルスシーケンス情報から取得する。算出部は、移動距離を経過時間で除算することで流体の流速を算出する。例えば、これらの部は、ＭＲＩシステム制御部２２内に備えられ（図示を省略）、ＭＲＩシステム制御部２２内に備えられたこれらの部が、ＭＲＩシーケンス制御部３０や架台部１０、その他関連する構成要素を制御する。以下、更にケースに分けて説明する。具体的には、心時相の異なる複数の流体画像を収集するケース１、及び、Ｔｉｍｅ−ＳＬＩＰ（Time-Spatial Labeling Inversion Pulse）撮像法によって時相の異なる複数の流体画像を収集するケース２を説明する。
【００２５】
なお、以下では、実施形態の例としていくつかのケースを説明するが、以下のケースに限定されるものではない。
【００２６】
［ケース１］
まず、ケース１の場合について説明する。
●ＦＡＳＥ（Fast Asymmetric Spin Echo）（又はＦＢＩ（Fresh Blood Imaging））シーケンスを使用する。位相エンコード（ＰＥ）方向が、末梢血管流量に対して垂直である場合、（１ＲＲ〜ｎＲＲｓなど）それらの間の小さい増分を伴う各ＦＡＳＥ画像のＭＲ信号輝度は、図２に示すように変動する。
●（例えば、心収縮期画像の）低輝度信号を（例えば、心拡張期起動期画像の）高輝度信号から減算すると、図３に示すように、移動する物体のように見える流動ボリュームの合成画像を得る。
●映像で表示した場合、非造影時間分解ＭＲＡを見ることができる。更に、蛇行状の血管に沿ってＭＲ流量信号を追跡することによって、画像間の時間差によって除算された移動血管長から速度が得られる。特定の（増分）流量速度ｖ１、ｖ２〜ｖｎの平均値から平均流量速度（もちろん、あらゆる中間の特定速度の計算を省略したい場合は、直接計算することもできる）が得られる。
【００２７】
ケース１に係る収集部は、被検体内を移動する流体の画像であって心時相の異なる複数の流体画像を収集する。ケース１に係る特定部は、複数の流体画像を用いて流体の移動距離を特定する。ケース１に係る取得部は、移動距離に対応する経過時間を、複数の流体画像の収集に用いられたパルスシーケンス情報から取得する。ケース１に係る算出部は、移動距離を経過時間で除算することで流体の流速を算出する。例えば、これらの部は、ＭＲＩシステム制御部２２内に備えられ（図示を省略）、ＭＲＩシステム制御部２２内に備えられたこれらの部が、ＭＲＩシーケンス制御部３０や架台部１０、その他関連する構成要素を制御する。
【００２８】
収集部は、例えばＦＢＩ撮像法を用いて、心収縮期と心拡張期との間にて心時相の異なる複数の流体画像を収集する。ＦＢＩ撮像法は、３Ｄ ＦＡＳＥによる血管撮像法であり、同期信号（例えばＲ波）からの適切な遅延時間を設定し、心電同期又は脈波同期を用いて収集を行うことで、心臓から拍出される新しい血液を描出する撮像法である。収集部によって収集されるＭＲ信号の輝度は、心収縮期と心拡張期との間で図２に示すように変動する。そこで、収集部は、心電同期信号からの遅延が徐々に増える遅延時間ｔ0msec、ｔ1msec、・・・を設定し、遅延時間ｔ0msecで同期してＭＲ信号Ｓ1を収集し、遅延時間ｔ1msecで同期してＭＲ信号Ｓ2を収集するというように、心時相の異なる複数の流体画像を収集する。
【００２９】
特定部は、例えば、所定の心時相にて収集された流体画像と、基準となる心時相にて収集された流体画像との差分画像を用いて、心時相の異なる流体画像毎に、流体の移動距離を特定する。上述したように、心時相の異なる流体画像間ではＭＲ信号の輝度が異なる。このため、例えば、基準となる心時相にて収集された流体画像から所定の心時相にて収集された流体画像を減算することで、その間に被検体内を移動した流体（例えば血液）の信号を描出することができる。例えば、図３において、「ｔn」は基準の心時相を示し、「Ｓ１（ｔn）」は心時相ｔnにて収集されたＭＲ信号を示す。図３に示すように、例えば、心収縮期にて心臓から拍出された血液のＭＲ信号は、その輝度が低いので（図３において例えば白色で表現）、遅延時間の増分が大きくなるほど、徐々に輝度の低い部分が増える。画像１〜画像６は、基準となる心時相の流体画像から各時相の流体画像を差し引いた画像であり、血液以外の情報が差し引かれ、血液のみが描出された画像となる。特定部は、例えば、この差分画像である画像１〜画像６を輝度に基づいて解析し、例えば、輝度の高い部分と輝度の低い部分とを区別することで、各心時相における流体の移動距離を特定する。例えば、特定部は、図４に示すように、移動距離Ｌ２〜Ｌ６を特定する。
【００３０】
取得部は、例えば、各流体画像の移動距離毎に、経過時間をパルスシーケンス情報から取得する。例えば、ケース１において、各心時相に対応する経過時間は、パルスシーケンス情報として設定された遅延時間に相当する。このため、取得部は、パルスシーケンス情報として設定された遅延時間を取得する。例えば、取得部は、図４に示すように、遅延時間として、ｔn+1msec、ｔn+2msec、ｔn+3msec、ｔn+4msec、ｔn+5msecを取得する。
【００３１】
そして、算出部は、例えば、各移動距離及び各経過時間を用いて流速を算出する。例えば、算出部は、各移動距離を累積した累積異動距離を、各経過時間を累積した累積経過時間で除算することで、流速を算出する。例えば、算出部は、図６に示すように、各心時相で累積した累積移動距離Ｌ６を累積経過時間（ｔn+5−ｔn）で除算することで速度を算出する。なお、速度の算出手法はこれに限られるものではない。算出部は、ある移動距離を、この移動距離に対応する経過時間で除算することで、その心時相に特定の特定速度を算出してもよい。例えば、図４に示すように、算出部は、移動距離Ｌ３を、経過時間（ｔn+2−ｔn）で除算することで、平均速度ＭＶ３を算出する。例えば、図５に示すように、算出部は、画像６と画像５との差分である移動距離ΔＬ６を、経過時間（ｔn+5−ｔn+4）で除算することで、特定速度ＳＶ（Specific Velocity）を算出してもよい。
【００３２】
［ケース２］
続いて、ケース２の場合について説明する。
●流出（非選択及び選択パルス）又は流入（選択パルスのみ）Ｔｉｍｅ−ＳＬＩＰ、並びに、別の取得及び減算（タグ付け及びタグ付けなし減算）もまた、平均速度の測定を可能にする。
○シングルショットＦＳＥ（ＦＡＳＥ）を使用する場合、速度は、流量移動時間としてのＢＢＴＩ時間（取得パルスにタグ付けする）及び有効ＴＥ（ＴＥｅｆｆ）内に移動した流量ボリュームを使って計算できる。したがって、速度は、（ＢＢＴＩ＋ＴＥｅｆｆ）によって除算された移動流体の長さ（Ｌ）として計算できる。
○ｂＳＳＦＰを使用する場合、速度は、流量移動時間としてのＢＢＴＩ時間（取得パルスにタグ付けする）内に移動した流量を使って計算できる。したがって、速度は、ＢＢＴＩで除算された移動流体の長さ（Ｌ）として計算できる。
【００３３】
ケース２に係る収集部は、被検体内を移動する流体のスピンに標識化を行い、所定時間経過後にスピンのエコー信号を受信する撮像を、所定時間を異ならせて複数回行い、流体の画像であって時相の異なる複数の流体画像を収集する。ケース２に係る特定部は、複数の流体画像を用いて、流体の移動距離を特定する。ケース２に係る取得部は、移動距離に対応する経過時間を、流体画像の収集に用いられたパルスシーケンス情報から取得する。ケース２に係る算出部は、移動距離を経過時間で除算することで流体の流速を算出する。例えば、これらの部は、ＭＲＩシステム制御部２２内に備えられ（図示を省略）、ＭＲＩシステム制御部２２内に備えられたこれらの部が、ＭＲＩシーケンス制御部３０や架台部１０、その他関連する構成要素を制御する。
【００３４】
ケース２の場合、収集部は、例えばＴｉｍｅ−ＳＬＩＰ撮像法を用いて、時相の異なる複数の流体画像を収集する。Ｔｉｍｅ−ＳＬＩＰ撮像法は、撮像領域に流入又は流出する流体を、この撮像領域とは独立した位置でラベリング（標識化）し、撮像領域に流入又は流出する流体の信号値を高く又は低くすることで、流体を描出する撮像法である。Ｔｉｍｅ−ＳＬＩＰ撮像法においては、同期信号（例えばＲ波）から一定の待ち時間後にＴｉｍｅ−ＳＬＩＰパルスを印加する。このＴｉｍｅ−ＳＬＩＰパルスは、領域非選択インバージョンパルス（図７において「非選択パルス」）及び領域選択インバージョンパルス（図７において「選択パルス」）を含み、領域非選択インバージョンパルスは、オン又はオフを設定することができる。領域選択インバージョンパルスによって撮像領域に流入（又は流出）する流体をラベリングすると、ＢＢＴＩ（Black-Blood Time to Inversion）時間後に流体が到達した部分の信号の輝度が高く（領域非選択インバージョンパルスがオフの場合は低く）なる。そこで、収集部は、ＢＢＴＩ時間を複数設定し、時相の異なる複数の流体画像を収集する。
【００３５】
なお、ラベリング位置が撮像領域外に設定された場合、ラベリングされた流体は撮像領域に流入するので、ここではこれを「フローイン」と呼ぶ。一方、ラベリング位置が撮像領域内に設定された場合、ラベリングされた流体は撮像領域に流出するので、ここではこれを「フローアウト」と呼ぶ。実施形態は、「フローイン」及び「フローアウト」のいずれにも適用することができる。
【００３６】
ここで、ケース２に係る収集部は、領域選択インバージョンパルスによってラベリングを行う収集と、領域選択インバージョンパルスによってラベリングを行わない収集とを、例えばスライスエンコード毎に交互に繰り返すことで２つの流体画像を収集する。なお、ラベリングを行わない収集については、スライスエンコード毎に毎回行わず、例えば１回行うものとしてもよい。更に、実施形態は、ラベリングを行わない収集を行わない場合にも、適用することができる。
【００３７】
そして、ケース２に係る特定部は、ラベリングを行わない撮像により収集された画像である基準画像と、ラベリングを行う撮像により収集された各流体画像との差分画像を用いて、時相の異なる流体画像毎に、流体の移動距離を特定する。例えば、特定部は、同一スライスエンコードについて収集された２つの流体画像間の差分をとることで、ラベリングされた部分のみを抽出し、背景信号を抑制する。あるいは、ラベリングを行わない収集が１回である場合は、特定部は、この１回に収集された基準画像と、ラベリングを行う撮像により収集された各流体画像との差分をとることで、ラベリングされた部分のみを抽出し、背景信号を抑制する。すると、図７に示すように、差分画像において、ＢＢＴＩ時間の増分が大きくなるほど、徐々に輝度の高い部分が増える（図７において例えば黒色で表現）。特定部は、例えば、各流体画像を輝度に基づいて解析し、例えば、輝度の高い部分と輝度の低い部分とを区別することで、各時相における流体の移動距離を特定する。
【００３８】
取得部は、各流体画像の移動距離毎に、経過時間をパルスシーケンス情報から取得する。例えば、ケース２において、収集部が、ＦＡＳＥ（Fast Asymmetric Spin Echo）法を用いたＴｉｍｅ−ＳＬＩＰ撮像法によって、時相の異なる複数の流体画像を収集した場合を想定する。この場合、取得部は、経過時間として、パルスシーケンス情報から、ＢＢＴＩ時間と実効エコー時間（TEeff（Effective Time to Echo））とを加算した値を取得する。
【００３９】
また、ケース２に係る取得部は、流体画像の生成に用いるｋ空間の充填方法に応じて、移動距離に対応する経過時間を取得してもよい。収集部が、ｂＳＳＦＰ（balanced Steady State Free Precession）法を用いたＴｉｍｅ−ＳＬＩＰ撮像法によって、時相の異なる複数の流体画像を収集した場合を想定する。この場合、取得部は、ｋ空間の中央から位相エンコードが配列されるセントリックオーダリングの場合には、経過時間として、パルスシーケンス情報からＢＢＴＩ時間を取得する。一方、取得部は、ｋ空間に位相エンコードが順次配列されるシーケンシャルオーダリングの場合には、経過時間として、パルスシーケンス情報から、ＢＢＴＩ時間と位相エンコード数の半数分に相当する時間とを加算した値を取得する。これらの経過時間は、ｋ空間の中央付近に充填されるＭＲ信号が収集された経過時間に対応することを意図するものである。
【００４０】
算出部は、ケース１と同様、例えば、各移動距離を累積した累積異動距離を、各経過時間を累積した累積経過時間で除算することで、流速を算出する。
【００４１】
なお、ケース２において、ラベリングの手法は、例えば、ラベリングのパルスを連続的に照射するｐＣＡＳＬ（Pulsed Continuous Arterial Spin Labeling）手法を用いてもよい。あるいは、ラベリングのパルスを連続的には照射しないものの、ラベリングの領域（標識化領域）の幅や大きさを厚く設定してもよい。
【００４２】
なお、ケース１及びケース２にかかわらず、望ましくは、全ての計算は、例えば、マウスクリックによって「速度」（「平均速度」及び「特定速度」の両方又は一方）のアイコンを選択することによって、ＭＲ画像コンソール上で行われる。
【００４３】
当業者には明白なように、非造影Ｔｉｍｅ−ＳＬＩＰ（ＡＳＬ）技術では、ＭＲＩシーケンスで複数の異なるＢＢＴＩ時間を使用することにより流量流体を観察する。
【００４４】
この実施形態では、ＭＲＩ表示システムから独立した計算ツールを使ってオフラインで計算するようなことは必要なく、速度測定（特定及び平均の両方又は一方）を容易にするための比較的単純なＧＵＩが提供される。このＧＵＩは、平均及び特定の両方又は一方の流量速度の容易かつ敏速な計算を可能にする。
【００４５】
図２に示すように、ＦＢＩ／Ｎａｖｉルーチンなどは、一般的なＰＱＲＳＴＵ心臓周期の心収縮期から心拡張期までの連続的な遅延時間における連続的スライス画像Ｓ１、Ｓ２、Ｓ３〜を取得するために利用できる。１つの完全な心臓周期は、「Ｒ〜Ｒ」サイクルと呼ばれることもある。もちろん、Ｒ心臓パルスの周波数は、被検体の脈拍数に匹敵する。
【００４６】
なお、「ＦＢＩ／Ｎａｖｉルーチン」とは、ＥＣＧ−Ｐｒｅｐ撮像によって収集された画像からＦＢＩ撮像に適した遅延時間を求めるための機能である。ＥＣＧ−Ｐｒｅｐ撮像は、ＦＢＩ撮像法における遅延時間を設定するために、ＦＢＩ撮像法による撮像に先行して行われる、２Ｄ ＦＡＳＥ撮像である。ＥＣＧ−Ｐｒｅｐ撮像によって、同期信号（例えばＲ波）からの遅延時間を変化させながら心時相の異なる複数の画像が収集され、収集された複数の画像、又は、複数の画像に基づいて解析された信号値の推移が表示部に表示される。ＦＢＩ−Ｎａｖｉは、ＥＣＧ−Ｐｒｅｐ撮像によって収集された複数の画像を解析して信号値の変化が大きい領域を抽出し、抽出した領域について、基準となる画像と各画像との信号値差を求め、求めた信号値差をグラフ表示する。
【００４７】
心臓（又は末梢血管脈拍）ゲーティングを適切に使用して、心臓周期の連続的な遅延ポイントで被検体にＭＲＩシーケンスをかけることができる。ｔ０における心収縮期から測定しながら、各連続するスライスＳ１、Ｓ２〜Ｓｎ（ｍａｘ）に対する遅延増分を、ｔｎ（ｍａｘ）における最大時間まで図２に描写する。
【００４８】
当業者には自明なように、ピクセル毎に、これらの画像のうちの１つを他の画像から減算することによって、異なる画像を生成できる。ＭＲＡの最大造影を得るためには、一般に、最少信号輝度を有する心収縮期画像を、最大信号輝度を有する心拡張期画像から減算し、それによって、単一のＭＲＡ画像を得るようになる。
【００４９】
しかし、連続する画像Ｓ２、Ｓ３などが、心収縮期画像Ｓ１から減算された場合、被検体の動脈に沿って進む「物体」として殆んど拍動性の流体流量の前縁を描出できる連続した異なる画像を生成することになる。
【００５０】
図３は、理想化された完全に線形の動脈の概略図である。該動脈は、心収縮期において時間ｔ０に撮像され、次に、ｔｎ＋１、ｔｎ＋２などにおいて心拡張期に向けて漸増的に増加する時間間隔で再び連続的に撮像されている。結果として生じる一連の画像は、（２枚の同一画像が減算されるので）ｎｕｌｌ画像１で始まり、その後進んで、撮像時間が心拡張期により近づくゆえに、流動中の流体が絶えず増加するＭＲ信号輝度を発生させるにつれて、上部から始まり下部に向かって進む撮像された動脈部分内を流れる流体が明らかになる（図３）。事実上、時間ｔｎから始まって、多様な連続時間間隔を経て、あらゆる長さの動脈部を通り抜ける流動中の流体の前縁を識別できる。
【００５１】
したがって、図４に示すように、各連続する区間の終わりで、流れる流体が移動した距離の全長を測定でき、更に、各撮像間隔の終わりで、その全長を経過時間で除算することにより、平均速度の測定が可能になる。
【００５２】
あるいは、図５に示すように、画像間の各漸増する時間間隔に対する特定速度もまた、画像間に移動された漸増距離の計算と、適切な関連した時間間隔による除算とによって計算できる。画像間の時間間隔が等しいと仮定すると、図５から分かるように、特定の速度は、画像２と比べて画像３で増加し、更に、画像３と比べて画像４で増加している。しかし、特定速度は、明らかに、画像４と比べて画像５で減少し、更にその後、再び、画像５に比べて画像６で増加している。
【００５３】
図６に示すように、同じ組になった連続する画像１〜６もまた、（画像６を画像１と比較して）全区間にわたる標準又は平均の速度のみを計算するように分析できる。
【００５４】
事実上、映像モードで示すように、心拡張期−心収縮期の画像を使用すれば、非造影時間分解ＭＲＡは、今や流速測定を含むことができる。流動中の信号源を追跡することによって、蛇行状の血管に沿ってでも、移動された血管部分の長さを、該移動に要した時間で除算することによって、速度（間隔又はより小さな間隔群の間で測定した場合、一連の間隔及び一連の特定速度の両方又は一方にわたる標準速度か、平均速度か）が得られる。すなわち、一区間から次の区間へ漸増する特定速度の標準値によって、全画像を包含するより長い時間にわたる平均又は標準の速度が得られる。
【００５５】
図７は、流体流量のＴｉｍｅ−ＳＬＩＰ非造影撮像を概略的に描いている。選択パルス及び非選択パルスを、異なる大きさの長方形で表している。当業者には自明なように、異なるＢＢＴＩ間隔を使用することによって、心拡張期−心収縮期の異なる画像に対して以前の図で説明した異なる時間増分における連続する画像に類似した、動脈、静脈などの中の流体流量の連続する画像を得ることができる。明らかなように、動脈、静脈などに沿って移動された距離の測定、及び、関連する移動時間による除算もまた、この非造影流量撮像技術を使って、速度測定を行うことができる。
【００５６】
図８は、コンピュータプログラムコード構造の速度測定モジュールを概略的に描いている。このモジュールは、主要な又は監視の操作システムなどに関連するあらゆる好適な手段（例えば、操作者による、好適に表示されたアイコンまたメニュー上でのマウスクリック、タッチスクリーン上でのフィンガー操作入力、キーボード入力などの使用）によって実行できる。
【００５７】
ステップＳ０１で、スライスカウンタｎを１の値に初期設定する。ステップＳ０２で、（例えば、自動的に機械が実行する距離追跡アルゴリズムを使用するために、基本設定ダイアログなどに制御パラメータを設定することによる）距離測定の自動追跡システムのためのオプションが設定されたか否かを見るテストを行う。設定されている場合、次に、自動追跡機能が、ステップＳ０３で作動し、特定の速度ＳＶが、ステップＳ０４で、ｔｎ−１とその後のｔｎとにおいて取得されたスライス間で計算される。
【００５８】
自動追跡機能を望まない場合、次に、ステップＳ０５で、経路カウンタＰを１の値に初期設定する。その後、ステップＳ０６で、操作者が第１の位置を定義するために、待機ループが実行される。例えば、図９に示すように、操作者は、カーソルを（撮像時間ｔｎ−１に関連する）第１のポイントＰ１上に移動させて、所定のマウスボタンをクリックすることによって、動脈などの内の流体流量の表示画像上の第１の位置を定義できる。操作者が、特別に定義したポイントを初期経路位置、中間経路位置、又は最終経路位置であることを識別するために選択できるポップアップメニュー９３もまた設けることができる。例えば、図９に示すように、操作者定義のポイントＰ１〜Ｐ５が、（撮像時間ｔｎ−１とｔｎとの間の）蛇行曲線状の血管に沿った効果的なポイントに定義されているので、流量経路に沿った明確なポイントの間の直線部分は、スライス画像Ｓｎ−１からスライス画像Ｓｎまでの流量の距離に厳密に近似している。一旦第１の位置を（図８の）ステップＳ０６で定義すれば、経路カウンタＰは、図８に示すように、ステップＳ０７で増分される。その後、経路に沿った更なる位置を操作者が定義するにつれて、待機ループが実行されて、経路カウンタは、最終位置が定義されるまで（ステップＳ０８）増分される。そこで、ステップＳ０９で、経路パラメータＰ（ｍａｘ）に対する最大値を、経路カウンタのその時点の最新値として定義する。次に、ステップＳ１０で、特定速度ＳＶが、連続するスライス画像ＳｎとＳｎ−１との間の流量の増分として計算される。
【００５９】
次に、最終の時間ｔｎに関連する特定速度ＳＶ（どのように計算されようとも）がステップＳ１１で表示され、格納され、又は（例えば、システムの好ましい部分に望ましく構成されるような遠隔場所に）出力される。その後、移動平均速度ＭＶはステップＳ１２で計算されて、同様にステップＳ１３で表示され、格納され、又は（例えば、システム内に以前構成されたように）望ましいように出力される。ステップＳ１４ではスライスカウンタが、最終値（すなわち、スライスＳｎ（ｍａｘ））に達したか否かを見るテストを行う。到達していない場合は、次に、スライスカウンタはステップＳ１５で増分されて、制御は図８に示すように判断ステップＳ０２か、又は、直接ブロックステップＳ０６の経路カウンタの再初期設定に戻される（すなわち、操作者に、自動追跡システムを作動させたくない場合）。すべての有効なスライスデータが処理されると、このモジュールは図８に示すように終了する。
【００６０】
図１０に示すように、特定速度ＳＶ及び移動平均速度ＭＶの両方又は一方は、時間（又は、スライス数など）の関数としてグラフ表示できる。
【００６１】
［表示］
ところで、実施形態に係るＭＲＩシステムは、流体画像、及び、流体画像との差分処理の対象とされる基準画像の少なくとも一つに、流速に関する情報を含めて表示部に表示する表示制御部を更に備えてもよい。例えば、表示制御部は、ＭＲＩシステム制御部２２内に備えられ（図示を省略）、ＭＲＩシステム制御部２２内に備えられた表示制御部が、関連する構成要素を制御する。実施形態において、収集部によって収集された流体画像は、例えばＦＢＩ撮像法によって収集された流体画像やＴｉｍｅ−ＳＬＩＰ撮像法によって収集された流体画像であり、分解能の高い画像であるといえる。実施形態に係るＭＲＩシステムは、この分解能の高い形態画像自体から、機能情報としての「流体の速度」を算出するので、形態画像と、機能情報の算出根拠とが一致していることになる。
【００６２】
そこで、例えば、表示制御部は、形態画像を表示しながら、その形態画像上に、流体の速度の情報を重畳して表示してもよい。例えば、表示制御部は、下肢の血管を描出した形態画像上において、流速が存在する箇所には流速が存在することを示す表示（例えば、下肢の血管の塗り潰し）を行い、流速が存在しない箇所には流速が存在しないことを示す表示（例えば、下肢の血管の斜線塗り潰し）を行ってもよい。また、例えば、表示制御部は、３Ｄの形態画像を表示しながら、その３Ｄの形態画像上に、２Ｄで、流体の映像データを重畳して表示（例えば、図３に例示した画像１〜画像６を連続再生して表示部に表示）してもよい。
【００６３】
以上述べた少なくとも一つの実施形態の磁気共鳴イメージング装置及び磁気共鳴イメージング方法によれば、流体の流速を求めることが可能になる。
【００６４】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【００６５】
８ 心臓信号取得装置
９ 被検体
１０ 架台部
１１ 被検体テーブル
１２ 静磁場Ｂ0磁石
１３ 心電ゲート信号
１４ Ｇｘ、Ｇｙ、Ｇｚ傾斜磁場コイルセット
１６ ＲＦコイルアセンブリ
１８ 撮像ボリューム
２０ システム構成要素
２２ ＭＲＩシステム制御部
２４ 表示部
２６ キーボード／マウス
２８ プリンタ
３０ ＭＲＩシーケンス制御部
３２ Ｇｘ、Ｇｙ、Ｇｚ傾斜磁場コイルドライバ
３４ ＲＦ送信部
３６ 送信／受信スイッチ
３８ プログラムコード構造
４０ ＲＦ受信部
４２ ＭＲＩデータ処理部
４４ 画像再構成プログラムコード構造
４６ ＭＲ画像記憶部
５０ ＭＲＩシステムプログラム／データ格納部
【技術分野】
【０００１】
本発明の実施形態は、磁気共鳴イメージング装置及び磁気共鳴イメージング方法に関する。
【背景技術】
【０００２】
従来、磁気共鳴イメージング装置（以下、ＭＲＩ（Magnetic Resonance Imaging）システム）による撮像法のひとつに、造影剤を用いずに被検体内を移動する流体を撮像する手法がある。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００８−２３３１７号公報
【特許文献２】米国特許出願公開第２００８／０００９７０５号明細書
【非特許文献】
【０００４】
【非特許文献１】Miyazaki, et al., “Non-contrast-enhanced MR angiography using 3D ECG-synchronized half-Fourier fast spin echo,” JMRI 12:776-783 (2000)
【非特許文献２】Suga, et al., Lung perfusion impairments in pulmonary embolic and airway obstruction with non contrast MR imaging,” J Appl Physiol 92:2439-2451 (2002)
【非特許文献３】Takahashi, et al., “Non-contrast-enhanced renal MRA using time-spatial labeling pulse(t-SLIP) with 3D balanced SSFP,” presented at the ISMRM 15th Annual Meeting, Berlin, Germany, page 179(2007)
【非特許文献４】Yamamoto, et al., “Non-contrast-enhanced MRDSM of peripheral arteries using continuous acquisitions of ECG-triggered 2D half-Fourier FSE within a cardiac cycle,” 12th Annual Meeting, Toronto, Canada, page 1709(2003)
【非特許文献５】Kanazawa, et al., “Time-spatial labeling inversion tag (t-SLIT) using a selective IR-tag on/off pulse in 2D and 3D half-Fourier FSE as arterial spin labeling,” presented at the ISMRM 10th Annual Meeting, Hawaii, page 140(2002)
【非特許文献６】Furudate, et al., “FBI-Navi for Easy Determination of Diastolic and Systolic Triggering Phase in Non-Contrast Fresh Blood Imaging(FBI),” ISMRM 16th Annual Meeting, Toronto, Canada, page 2902(2008)
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明が解決しようとする課題は、流体の流速を求めることである。
【課題を解決するための手段】
【０００６】
実施形態に係る磁気共鳴イメージング装置は、収集部と、特定部と、取得部と、算出部とを備える。前記収集部は、被検体内を移動する流体の画像である流体画像を複数収集する。前記特定部は、前記複数の流体画像のうちの１つである基準画像と各流体画像との差分画像を用いて、前記流体の移動距離を特定する。前記取得部は、前記移動距離に対応する経過時間を、前記複数の流体画像の収集に用いられたパルスシーケンス情報から取得する。前記算出部は、前記移動距離を前記経過時間で除算することで前記流体の流速を算出する。
【図面の簡単な説明】
【０００７】
【図１】図１は、実施形態に係るＭＲＩシステムの概略ブロック図である。
【図２】図２は、実施形態における心収縮期及び心拡張期に対する連続的な心電同期ＭＲスライス撮像シーケンスを示すタイミング図である。
【図３】図３は、図２に示すように、連続する心臓ゲーティング時間増分で「暗い（信号の輝度が低い）」心収縮期画像と「明るい（信号の輝度が高い）」心拡張期画像との間の差分により得られる連続する差分画像を示す概略図である。
【図４】図４は、図３の概略図と類似の概略図であるが、実施形態によって、標準又は平均の血流速度を計算できる方法を説明する注釈を含む。
【図５】図５は、図３〜図４とは異なる画像を示す類似の概略図であるが、実施形態によって、特定の速度を各連続する区間で計算できる方法を説明する注釈を含む。
【図６】図６は、図５と類似であるが、実施形態の一連の異なる画像の全体にわたって、全標準又は平均の速度を計算できる方法を具体的に示している。
【図７】図７は、Ｔｉｍｅ−ＳＬＩＰ（Time-Spatial Labeling Inversion Pulse）撮像法（流出及びタグ付け／タグ付けなし減算法）を使う別の実施形態による速度測定の概略を示す。
【図８】図８は、実施形態で利用可能なコンピュータプログラムコード構造の一例に関する概略フローチャートである。
【図９】図９は、２つの連続する撮像時間の間の蛇行状の経路に沿った血流に対する連続的な経路位置ポイントを操作者が定義する際に利用できる操作者画面の概略図である。
【図１０】図１０は、実施形態における時間及びスライス位置の両方又は一方に対する特定の速度及び平均速度の両方又は一方として考えられる出力表示グラフの概略図である。
【発明を実施するための形態】
【０００８】
図１に示すＭＲＩ（Magnetic Resonance Imaging）システムは、架台部１０（断面図にて示す）と、互いに接続される様々な関連のシステム構成要素２０とを含む。少なくとも架台部１０は、通常シールドルーム内に設置される。図１に示す１つのＭＲＩシステムは、静磁場Ｂ0磁石１２と、Ｇx、Ｇy、Ｇz傾斜磁場コイルセット１４と、ＲＦ（Radio Frequency）コイルアセンブリ１６との実質的に同軸円筒状の配置を含む。この円筒状に配置された要素の水平軸線に沿って、被検体テーブル１１によって支持された被検体９の頭部を取り囲むように示された撮像ボリューム１８がある。
【０００９】
ＭＲＩシステム制御部２２は、表示部２４、キーボード／マウス２６、及びプリンタ２８に接続された入力／出力ポートを備える。言うまでもなく、表示部２４は、制御入力もまた備えるような多様性のあるタッチスクリーンであってもよい。
【００１０】
ＭＲＩシステム制御部２２は、ＭＲＩシーケンス制御部３０とインタフェース接続する。ＭＲＩシーケンス制御部３０は、Ｇx、Ｇy、Ｇz傾斜磁場コイルドライバ３２、並びに、ＲＦ送信部３４及び送信／受信スイッチ３６（同じＲＦコイルが送信及び受信の両方に使用されている場合）を順に制御する。ＭＲＩシーケンス制御部３０は、既にＭＲＩシーケンス制御部３０にて可能な、ＭＲＩデータ取得シーケンス（血液等の流動する流体（flowing fluid）の撮像を含む）の実行に適切なプログラムコード構造３８を含む。心臓信号取得装置８（例えば、被検体９の身体に適切に装着された、ＥＣＧ（electrocardiogram）又は末梢血管パルストランスデューサ（単数又は複数））は、ＭＲＩシーケンス制御部３０のトリガー信号とする心電ゲート信号１３を出力することができる。
【００１１】
ＭＲＩシステムは、表示部２４に出力する処理された画像データを作成できるように、ＭＲＩデータ処理部４２に入力を供給するＲＦ受信部４０を含む。また、ＭＲＩデータ処理部４２を、画像再構成プログラムコード構造４４及びＭＲ画像記憶部４６にアクセスできるように構成してもよい（例えば、実施形態及び画像再構成プログラムコード構造４４に従った処理から得られたＭＲＩデータを格納するために）。
【００１２】
また、図１は、ＭＲＩシステムプログラム／データ格納部５０を一般化した描写を示す。ＭＲＩシステムプログラム／データ格納部５０に格納されるプログラムコード構造（例えば、画像再構成のため、及び、流速、距離／時間の測定のため、操作者入力等のための）は、ＭＲＩシステムの様々なデータ処理構成要素にアクセス可能なコンピュータ読み取り可能な記憶媒体に格納される。当業者には自明であるが、ＭＲＩシステムプログラム／データ格納部５０は、正常運転時にそのように格納されたプログラムコード構造に対して直近の必要性を有するＭＲＩシステムの処理コンピュータのうちの様々なコンピュータに分割し、且つ少なくとも一部を直結してもよい（すなわち、ＭＲＩシステム制御部２２に普通に格納したり、直結したりする代わりに）。
【００１３】
実際、当業者には自明であるが、図１の描写は、本明細書で後述する実施形態を実行できるように若干の変更を加えた一般的なＭＲＩシステムを非常に高度に簡素化した図である。システム構成要素は、様々な論理収集の「ボックス」に分割でき、通常、多数のデジタル信号処理装置（ＤＳＰ（Digital Signal Processors））、超小型演算処理装置、特殊用途向け処理回路（例えば、高速Ａ／Ｄ変換、高速フーリエ変換、アレイ処理用等）を含む。これら処理装置のそれぞれは、通常、各クロックサイクル（又は所定数のクロックサイクル）が発生すると、物理データ処理回路がある物理的状態から別の物理的状態へ進むクロック動作型の「状態機械」である。
【００１４】
動作中に、処理回路（例えば、ＣＰＵ（Central Processing Unit）、レジスタ、バッファ、計算ユニット等）の物理的状態が、あるクロックサイクルから別のクロックサイクルへ漸進的に変化するだけでなく、連結されているデータ格納媒体（例えば、磁気記憶媒体のビット格納部）の物理的状態も、そのようなシステムの動作中に、ある状態から別の状態へ変わる。例えば、ＭＲＩ再構成プロセスの終了時、物理的記憶媒体のコンピュータ読み取り可能なアクセス可能データ値格納場所のアレイは、いくつかの事前の状態（例えば、全部一律の「ゼロ」値又は全部「１」値）から新しい状態に変わる。その新しい状態では、そのようなアレイの物理的場所の物理的状態は、最小値と最大値との間で変動し、現実世界の物理的事象及び状況（例えば、撮像ボリューム空間内の被検体の動脈内を流れる血液）を表現する。当業者には自明であるが、格納されたデータ値のそのようなアレイは、物理的構造を表し且つ構成もする。つまり、命令レジスタの中に順次読み込まれてＭＲＩシステムの１つ以上のＣＰＵによって実行されたとき、動作状態の特定シーケンスが発生して、ＭＲＩシステム内中に移行されるコンピュータ制御プログラムコードの特定構造が構成される。
【００１５】
下記の実施形態は、流動する流体の速度（以下、適宜「流速」という）を含む、データ取得の処理及びＭＲ画像の生成表示の両方又は一方を行うための改良された方法を提供する。
【００１６】
心収縮期及び心拡張期のタイミングと同期して収集された画像間には、ＭＲ信号強度の変化がある。このＭＲ信号強度の変化は、例えば、信号値の高い心拡張期の画像から信号値の低い心収縮期の画像を減算すること、及び、信号値の低い心収縮期の画像から信号値の高い心拡張期の画像を減算することの両方あるいは一方を行うことによって、時間分解された流体（例えば、血液）の流れを非造影にて得るために用いることができる。これは、時間分解非造影ＭＲＡ（磁気共鳴血管造影法（Magnetic Resonance Angiography））と呼ばれることもある。また、流体の流れのための他の非造影撮像技術として、ＡＳＬ（Arterial Spin Labeling）と呼ばれることもあるＴｉｍｅ−ＳＬＩＰ（Time−Spatial Labeling Inversion Pulse）がある。
【００１７】
これらの非造影撮像技術においては、流体の流速測定を容易にするための操作者インタフェースが、望まれるであろう。
【００１８】
以下に示す実施形態によれば、非造影の流体画像における流体の流速測定は、少なくとも２つの手法で実現することができる。第１は、心拡張期と心収縮期との減算画像を用いる手法である。ＭＲＩシステムは、心拡張期と心収縮期との減算画像を用いて、流体の移動距離を特定することができる。例えば、動脈経路が比較的線形である場合、開始ポイントから最終ポイントまでの経路を単純な線であると考えることができる。あるいは、動脈経路が蛇行状の場合でも、経路に沿って効果的に置かれたポイントをつなぐことによって、開始ポイントから最終ポイントまでの経路を特定することができる。任意に形成された経路に沿って流動中の流体が移動した距離を自動的に追跡することも有効であるといえる。とにかく、流動中の流体が連続する撮像時間の間に通過した距離は、連続する撮像時間の間の特定の速度、並びに、一連のそのような撮像時間にわたる標準又は平均の速度を算出するために使用することができる。
【００１９】
第２は、非造影Ｔｉｍｅ−ＳＬＩＰ撮像技術を使用する手法である。非造影Ｔｉｍｅ−ＳＬＩＰ撮像技術を使用する場合、ｋ空間に対する位相エンコードの配列がセントリックオーダリングかシーケンシャルオーダリングかによって、ＦＡＳＥ（Fast Asymmetric Spin Echo）又はｂＳＳＦＰ（balanced Steady State Free Precession）などの若干異なる撮像シーケンスを含んでもよい。もちろん、流速の計算において分子は、やはり、信号を発生して流れる物質（例えば、血液）が移動した距離である。流速の計算において分母として使用される時間増分は、ＦＡＳＥ撮像法を用いたＴｉｍｅ−ＳＬＩＰの場合、ＢＢＴＩ（Black-Blood Time to Inversion）時間、プラス、実効エコー時間（ＴＥｅｆｆ）である。また、セントリックオーダリングでｂＳＳＦＰを用いたＴｉｍｅ−ＳＬＩＰの場合、ＢＢＴＩである。また、シーケンシャルオーダリングでｂＳＳＦＰを用いたＴｉｍｅ−ＳＬＩＰの場合、ＢＢＴＩ時間、プラス、１／２ＥＴＬ（エコートレイン長）（位相エンコード数の半数分に相当する時間）である。
【００２０】
速度測定は、異なるＥＣＧ（Electrocardiogram）（心電図ゲーティング）信号ポイント（例えば、心収縮期から心拡張期）間に流量変位を生成する、又は、ＦＢＩ（Fresh Blood Imaging）−Ｎａｖｉ信号差プロットを使用するＭＲ信号に対して行うことができる。流動するＭＲ信号変位又は移動したＭＲ信号の速度測定を計測できる。
【００２１】
上述のように、（ｋ空間配列に依存してＴＥｅｆｆ及びＥＴＬ／２の両方又は一方をおそらくプラスする）ＢＢＴＩ時間の間にＴｉｍｅ−ＳＬＩＰ撮像技術を使って流量を生成するＭＲ信号の変位もまた、平均速度の計算に使用できる。
【００２２】
一実施形態では、移動距離に沿った流速を選択的に計算するために（例えば、マウスボタンの押下によって）、簡単なＧＵＩ（Graphical User Interface）を使用できる。移動距離は、ユーザーによって（例えば、経路に沿って効果的に置かれたポイントを画定することによって）、又は、システム提供の自動追跡機能によって、輪郭を描出できる。とにかく、移動距離は、次に、（例えば、ＥＣＧ遅延撮像技術及びＴｉｍｅ−ＳＬＩＰ撮像技術の両方又は一方において減算された画像間の）その距離に関連する時間間隔によって除算される。
【００２３】
ＧＵＩは、平均速度及び特定速度の両方又は一方の流量移動の測定に使用できる。非造影技術では、速度測定は、少なくとも２つの方法で測定できる。第１の方法は、心拡張期−心収縮期減算撮像技術を使って、ＭＲ信号変位（最も近いポイント又は長く連結された線による）の距離に沿って変位線を描くことができる。ＭＲＩ表示システムは、自動的にこれを距離として記録して、次に、適切なボタンのクリック（又は、例えば、速度選択のための正しいマウスクリック）によって、（例えば、減算される心収縮期及び心拡張期の画像の有効経過時間の間の）適切な対応する時間差によってそのような距離を除算する計算を実現する。第２の方法（例えば、ＡＳＬ（arterial spin labeling）法）では、ＭＲＩ信号は、ＦＡＳＥを用いたＴｉｍｅ−ＳＬＩＰに対する（ＢＢＴＩ＋ＴＥｅｆｆ）、及びセントリックオーダリングでｂＳＳＦＰを用いたＴｉｍｅ−ＳＬＩＰに対する（ＢＢＴＩ）、及び（シーケンシャルオーダリングでｂＳＳＦＰを用いたＴｉｍｅ−ＳＬＩＰに対するＢＢＴＩ＋ＥＴＬ／２）によって除算された距離内で、最初の信号から最後の信号まで移動する流量を生成する。
【００２４】
さて、上述してきた実施形態に係るＭＲＩシステムを改めて説明すると、実施形態に係るＭＲＩシステムは、収集部と、特定部と、取得部と、算出部とを備える。収集部は、被検体内を移動する流体の画像である流体画像を複数収集する。特定部は、複数の流体画像のうちの１つである基準画像と各流体画像との差分画像を用いて、流体の移動距離を特定する。取得部は、移動距離に対応する経過時間を、複数の流体画像の収集に用いられたパルスシーケンス情報から取得する。算出部は、移動距離を経過時間で除算することで流体の流速を算出する。例えば、これらの部は、ＭＲＩシステム制御部２２内に備えられ（図示を省略）、ＭＲＩシステム制御部２２内に備えられたこれらの部が、ＭＲＩシーケンス制御部３０や架台部１０、その他関連する構成要素を制御する。以下、更にケースに分けて説明する。具体的には、心時相の異なる複数の流体画像を収集するケース１、及び、Ｔｉｍｅ−ＳＬＩＰ（Time-Spatial Labeling Inversion Pulse）撮像法によって時相の異なる複数の流体画像を収集するケース２を説明する。
【００２５】
なお、以下では、実施形態の例としていくつかのケースを説明するが、以下のケースに限定されるものではない。
【００２６】
［ケース１］
まず、ケース１の場合について説明する。
●ＦＡＳＥ（Fast Asymmetric Spin Echo）（又はＦＢＩ（Fresh Blood Imaging））シーケンスを使用する。位相エンコード（ＰＥ）方向が、末梢血管流量に対して垂直である場合、（１ＲＲ〜ｎＲＲｓなど）それらの間の小さい増分を伴う各ＦＡＳＥ画像のＭＲ信号輝度は、図２に示すように変動する。
●（例えば、心収縮期画像の）低輝度信号を（例えば、心拡張期起動期画像の）高輝度信号から減算すると、図３に示すように、移動する物体のように見える流動ボリュームの合成画像を得る。
●映像で表示した場合、非造影時間分解ＭＲＡを見ることができる。更に、蛇行状の血管に沿ってＭＲ流量信号を追跡することによって、画像間の時間差によって除算された移動血管長から速度が得られる。特定の（増分）流量速度ｖ１、ｖ２〜ｖｎの平均値から平均流量速度（もちろん、あらゆる中間の特定速度の計算を省略したい場合は、直接計算することもできる）が得られる。
【００２７】
ケース１に係る収集部は、被検体内を移動する流体の画像であって心時相の異なる複数の流体画像を収集する。ケース１に係る特定部は、複数の流体画像を用いて流体の移動距離を特定する。ケース１に係る取得部は、移動距離に対応する経過時間を、複数の流体画像の収集に用いられたパルスシーケンス情報から取得する。ケース１に係る算出部は、移動距離を経過時間で除算することで流体の流速を算出する。例えば、これらの部は、ＭＲＩシステム制御部２２内に備えられ（図示を省略）、ＭＲＩシステム制御部２２内に備えられたこれらの部が、ＭＲＩシーケンス制御部３０や架台部１０、その他関連する構成要素を制御する。
【００２８】
収集部は、例えばＦＢＩ撮像法を用いて、心収縮期と心拡張期との間にて心時相の異なる複数の流体画像を収集する。ＦＢＩ撮像法は、３Ｄ ＦＡＳＥによる血管撮像法であり、同期信号（例えばＲ波）からの適切な遅延時間を設定し、心電同期又は脈波同期を用いて収集を行うことで、心臓から拍出される新しい血液を描出する撮像法である。収集部によって収集されるＭＲ信号の輝度は、心収縮期と心拡張期との間で図２に示すように変動する。そこで、収集部は、心電同期信号からの遅延が徐々に増える遅延時間ｔ0msec、ｔ1msec、・・・を設定し、遅延時間ｔ0msecで同期してＭＲ信号Ｓ1を収集し、遅延時間ｔ1msecで同期してＭＲ信号Ｓ2を収集するというように、心時相の異なる複数の流体画像を収集する。
【００２９】
特定部は、例えば、所定の心時相にて収集された流体画像と、基準となる心時相にて収集された流体画像との差分画像を用いて、心時相の異なる流体画像毎に、流体の移動距離を特定する。上述したように、心時相の異なる流体画像間ではＭＲ信号の輝度が異なる。このため、例えば、基準となる心時相にて収集された流体画像から所定の心時相にて収集された流体画像を減算することで、その間に被検体内を移動した流体（例えば血液）の信号を描出することができる。例えば、図３において、「ｔn」は基準の心時相を示し、「Ｓ１（ｔn）」は心時相ｔnにて収集されたＭＲ信号を示す。図３に示すように、例えば、心収縮期にて心臓から拍出された血液のＭＲ信号は、その輝度が低いので（図３において例えば白色で表現）、遅延時間の増分が大きくなるほど、徐々に輝度の低い部分が増える。画像１〜画像６は、基準となる心時相の流体画像から各時相の流体画像を差し引いた画像であり、血液以外の情報が差し引かれ、血液のみが描出された画像となる。特定部は、例えば、この差分画像である画像１〜画像６を輝度に基づいて解析し、例えば、輝度の高い部分と輝度の低い部分とを区別することで、各心時相における流体の移動距離を特定する。例えば、特定部は、図４に示すように、移動距離Ｌ２〜Ｌ６を特定する。
【００３０】
取得部は、例えば、各流体画像の移動距離毎に、経過時間をパルスシーケンス情報から取得する。例えば、ケース１において、各心時相に対応する経過時間は、パルスシーケンス情報として設定された遅延時間に相当する。このため、取得部は、パルスシーケンス情報として設定された遅延時間を取得する。例えば、取得部は、図４に示すように、遅延時間として、ｔn+1msec、ｔn+2msec、ｔn+3msec、ｔn+4msec、ｔn+5msecを取得する。
【００３１】
そして、算出部は、例えば、各移動距離及び各経過時間を用いて流速を算出する。例えば、算出部は、各移動距離を累積した累積異動距離を、各経過時間を累積した累積経過時間で除算することで、流速を算出する。例えば、算出部は、図６に示すように、各心時相で累積した累積移動距離Ｌ６を累積経過時間（ｔn+5−ｔn）で除算することで速度を算出する。なお、速度の算出手法はこれに限られるものではない。算出部は、ある移動距離を、この移動距離に対応する経過時間で除算することで、その心時相に特定の特定速度を算出してもよい。例えば、図４に示すように、算出部は、移動距離Ｌ３を、経過時間（ｔn+2−ｔn）で除算することで、平均速度ＭＶ３を算出する。例えば、図５に示すように、算出部は、画像６と画像５との差分である移動距離ΔＬ６を、経過時間（ｔn+5−ｔn+4）で除算することで、特定速度ＳＶ（Specific Velocity）を算出してもよい。
【００３２】
［ケース２］
続いて、ケース２の場合について説明する。
●流出（非選択及び選択パルス）又は流入（選択パルスのみ）Ｔｉｍｅ−ＳＬＩＰ、並びに、別の取得及び減算（タグ付け及びタグ付けなし減算）もまた、平均速度の測定を可能にする。
○シングルショットＦＳＥ（ＦＡＳＥ）を使用する場合、速度は、流量移動時間としてのＢＢＴＩ時間（取得パルスにタグ付けする）及び有効ＴＥ（ＴＥｅｆｆ）内に移動した流量ボリュームを使って計算できる。したがって、速度は、（ＢＢＴＩ＋ＴＥｅｆｆ）によって除算された移動流体の長さ（Ｌ）として計算できる。
○ｂＳＳＦＰを使用する場合、速度は、流量移動時間としてのＢＢＴＩ時間（取得パルスにタグ付けする）内に移動した流量を使って計算できる。したがって、速度は、ＢＢＴＩで除算された移動流体の長さ（Ｌ）として計算できる。
【００３３】
ケース２に係る収集部は、被検体内を移動する流体のスピンに標識化を行い、所定時間経過後にスピンのエコー信号を受信する撮像を、所定時間を異ならせて複数回行い、流体の画像であって時相の異なる複数の流体画像を収集する。ケース２に係る特定部は、複数の流体画像を用いて、流体の移動距離を特定する。ケース２に係る取得部は、移動距離に対応する経過時間を、流体画像の収集に用いられたパルスシーケンス情報から取得する。ケース２に係る算出部は、移動距離を経過時間で除算することで流体の流速を算出する。例えば、これらの部は、ＭＲＩシステム制御部２２内に備えられ（図示を省略）、ＭＲＩシステム制御部２２内に備えられたこれらの部が、ＭＲＩシーケンス制御部３０や架台部１０、その他関連する構成要素を制御する。
【００３４】
ケース２の場合、収集部は、例えばＴｉｍｅ−ＳＬＩＰ撮像法を用いて、時相の異なる複数の流体画像を収集する。Ｔｉｍｅ−ＳＬＩＰ撮像法は、撮像領域に流入又は流出する流体を、この撮像領域とは独立した位置でラベリング（標識化）し、撮像領域に流入又は流出する流体の信号値を高く又は低くすることで、流体を描出する撮像法である。Ｔｉｍｅ−ＳＬＩＰ撮像法においては、同期信号（例えばＲ波）から一定の待ち時間後にＴｉｍｅ−ＳＬＩＰパルスを印加する。このＴｉｍｅ−ＳＬＩＰパルスは、領域非選択インバージョンパルス（図７において「非選択パルス」）及び領域選択インバージョンパルス（図７において「選択パルス」）を含み、領域非選択インバージョンパルスは、オン又はオフを設定することができる。領域選択インバージョンパルスによって撮像領域に流入（又は流出）する流体をラベリングすると、ＢＢＴＩ（Black-Blood Time to Inversion）時間後に流体が到達した部分の信号の輝度が高く（領域非選択インバージョンパルスがオフの場合は低く）なる。そこで、収集部は、ＢＢＴＩ時間を複数設定し、時相の異なる複数の流体画像を収集する。
【００３５】
なお、ラベリング位置が撮像領域外に設定された場合、ラベリングされた流体は撮像領域に流入するので、ここではこれを「フローイン」と呼ぶ。一方、ラベリング位置が撮像領域内に設定された場合、ラベリングされた流体は撮像領域に流出するので、ここではこれを「フローアウト」と呼ぶ。実施形態は、「フローイン」及び「フローアウト」のいずれにも適用することができる。
【００３６】
ここで、ケース２に係る収集部は、領域選択インバージョンパルスによってラベリングを行う収集と、領域選択インバージョンパルスによってラベリングを行わない収集とを、例えばスライスエンコード毎に交互に繰り返すことで２つの流体画像を収集する。なお、ラベリングを行わない収集については、スライスエンコード毎に毎回行わず、例えば１回行うものとしてもよい。更に、実施形態は、ラベリングを行わない収集を行わない場合にも、適用することができる。
【００３７】
そして、ケース２に係る特定部は、ラベリングを行わない撮像により収集された画像である基準画像と、ラベリングを行う撮像により収集された各流体画像との差分画像を用いて、時相の異なる流体画像毎に、流体の移動距離を特定する。例えば、特定部は、同一スライスエンコードについて収集された２つの流体画像間の差分をとることで、ラベリングされた部分のみを抽出し、背景信号を抑制する。あるいは、ラベリングを行わない収集が１回である場合は、特定部は、この１回に収集された基準画像と、ラベリングを行う撮像により収集された各流体画像との差分をとることで、ラベリングされた部分のみを抽出し、背景信号を抑制する。すると、図７に示すように、差分画像において、ＢＢＴＩ時間の増分が大きくなるほど、徐々に輝度の高い部分が増える（図７において例えば黒色で表現）。特定部は、例えば、各流体画像を輝度に基づいて解析し、例えば、輝度の高い部分と輝度の低い部分とを区別することで、各時相における流体の移動距離を特定する。
【００３８】
取得部は、各流体画像の移動距離毎に、経過時間をパルスシーケンス情報から取得する。例えば、ケース２において、収集部が、ＦＡＳＥ（Fast Asymmetric Spin Echo）法を用いたＴｉｍｅ−ＳＬＩＰ撮像法によって、時相の異なる複数の流体画像を収集した場合を想定する。この場合、取得部は、経過時間として、パルスシーケンス情報から、ＢＢＴＩ時間と実効エコー時間（TEeff（Effective Time to Echo））とを加算した値を取得する。
【００３９】
また、ケース２に係る取得部は、流体画像の生成に用いるｋ空間の充填方法に応じて、移動距離に対応する経過時間を取得してもよい。収集部が、ｂＳＳＦＰ（balanced Steady State Free Precession）法を用いたＴｉｍｅ−ＳＬＩＰ撮像法によって、時相の異なる複数の流体画像を収集した場合を想定する。この場合、取得部は、ｋ空間の中央から位相エンコードが配列されるセントリックオーダリングの場合には、経過時間として、パルスシーケンス情報からＢＢＴＩ時間を取得する。一方、取得部は、ｋ空間に位相エンコードが順次配列されるシーケンシャルオーダリングの場合には、経過時間として、パルスシーケンス情報から、ＢＢＴＩ時間と位相エンコード数の半数分に相当する時間とを加算した値を取得する。これらの経過時間は、ｋ空間の中央付近に充填されるＭＲ信号が収集された経過時間に対応することを意図するものである。
【００４０】
算出部は、ケース１と同様、例えば、各移動距離を累積した累積異動距離を、各経過時間を累積した累積経過時間で除算することで、流速を算出する。
【００４１】
なお、ケース２において、ラベリングの手法は、例えば、ラベリングのパルスを連続的に照射するｐＣＡＳＬ（Pulsed Continuous Arterial Spin Labeling）手法を用いてもよい。あるいは、ラベリングのパルスを連続的には照射しないものの、ラベリングの領域（標識化領域）の幅や大きさを厚く設定してもよい。
【００４２】
なお、ケース１及びケース２にかかわらず、望ましくは、全ての計算は、例えば、マウスクリックによって「速度」（「平均速度」及び「特定速度」の両方又は一方）のアイコンを選択することによって、ＭＲ画像コンソール上で行われる。
【００４３】
当業者には明白なように、非造影Ｔｉｍｅ−ＳＬＩＰ（ＡＳＬ）技術では、ＭＲＩシーケンスで複数の異なるＢＢＴＩ時間を使用することにより流量流体を観察する。
【００４４】
この実施形態では、ＭＲＩ表示システムから独立した計算ツールを使ってオフラインで計算するようなことは必要なく、速度測定（特定及び平均の両方又は一方）を容易にするための比較的単純なＧＵＩが提供される。このＧＵＩは、平均及び特定の両方又は一方の流量速度の容易かつ敏速な計算を可能にする。
【００４５】
図２に示すように、ＦＢＩ／Ｎａｖｉルーチンなどは、一般的なＰＱＲＳＴＵ心臓周期の心収縮期から心拡張期までの連続的な遅延時間における連続的スライス画像Ｓ１、Ｓ２、Ｓ３〜を取得するために利用できる。１つの完全な心臓周期は、「Ｒ〜Ｒ」サイクルと呼ばれることもある。もちろん、Ｒ心臓パルスの周波数は、被検体の脈拍数に匹敵する。
【００４６】
なお、「ＦＢＩ／Ｎａｖｉルーチン」とは、ＥＣＧ−Ｐｒｅｐ撮像によって収集された画像からＦＢＩ撮像に適した遅延時間を求めるための機能である。ＥＣＧ−Ｐｒｅｐ撮像は、ＦＢＩ撮像法における遅延時間を設定するために、ＦＢＩ撮像法による撮像に先行して行われる、２Ｄ ＦＡＳＥ撮像である。ＥＣＧ−Ｐｒｅｐ撮像によって、同期信号（例えばＲ波）からの遅延時間を変化させながら心時相の異なる複数の画像が収集され、収集された複数の画像、又は、複数の画像に基づいて解析された信号値の推移が表示部に表示される。ＦＢＩ−Ｎａｖｉは、ＥＣＧ−Ｐｒｅｐ撮像によって収集された複数の画像を解析して信号値の変化が大きい領域を抽出し、抽出した領域について、基準となる画像と各画像との信号値差を求め、求めた信号値差をグラフ表示する。
【００４７】
心臓（又は末梢血管脈拍）ゲーティングを適切に使用して、心臓周期の連続的な遅延ポイントで被検体にＭＲＩシーケンスをかけることができる。ｔ０における心収縮期から測定しながら、各連続するスライスＳ１、Ｓ２〜Ｓｎ（ｍａｘ）に対する遅延増分を、ｔｎ（ｍａｘ）における最大時間まで図２に描写する。
【００４８】
当業者には自明なように、ピクセル毎に、これらの画像のうちの１つを他の画像から減算することによって、異なる画像を生成できる。ＭＲＡの最大造影を得るためには、一般に、最少信号輝度を有する心収縮期画像を、最大信号輝度を有する心拡張期画像から減算し、それによって、単一のＭＲＡ画像を得るようになる。
【００４９】
しかし、連続する画像Ｓ２、Ｓ３などが、心収縮期画像Ｓ１から減算された場合、被検体の動脈に沿って進む「物体」として殆んど拍動性の流体流量の前縁を描出できる連続した異なる画像を生成することになる。
【００５０】
図３は、理想化された完全に線形の動脈の概略図である。該動脈は、心収縮期において時間ｔ０に撮像され、次に、ｔｎ＋１、ｔｎ＋２などにおいて心拡張期に向けて漸増的に増加する時間間隔で再び連続的に撮像されている。結果として生じる一連の画像は、（２枚の同一画像が減算されるので）ｎｕｌｌ画像１で始まり、その後進んで、撮像時間が心拡張期により近づくゆえに、流動中の流体が絶えず増加するＭＲ信号輝度を発生させるにつれて、上部から始まり下部に向かって進む撮像された動脈部分内を流れる流体が明らかになる（図３）。事実上、時間ｔｎから始まって、多様な連続時間間隔を経て、あらゆる長さの動脈部を通り抜ける流動中の流体の前縁を識別できる。
【００５１】
したがって、図４に示すように、各連続する区間の終わりで、流れる流体が移動した距離の全長を測定でき、更に、各撮像間隔の終わりで、その全長を経過時間で除算することにより、平均速度の測定が可能になる。
【００５２】
あるいは、図５に示すように、画像間の各漸増する時間間隔に対する特定速度もまた、画像間に移動された漸増距離の計算と、適切な関連した時間間隔による除算とによって計算できる。画像間の時間間隔が等しいと仮定すると、図５から分かるように、特定の速度は、画像２と比べて画像３で増加し、更に、画像３と比べて画像４で増加している。しかし、特定速度は、明らかに、画像４と比べて画像５で減少し、更にその後、再び、画像５に比べて画像６で増加している。
【００５３】
図６に示すように、同じ組になった連続する画像１〜６もまた、（画像６を画像１と比較して）全区間にわたる標準又は平均の速度のみを計算するように分析できる。
【００５４】
事実上、映像モードで示すように、心拡張期−心収縮期の画像を使用すれば、非造影時間分解ＭＲＡは、今や流速測定を含むことができる。流動中の信号源を追跡することによって、蛇行状の血管に沿ってでも、移動された血管部分の長さを、該移動に要した時間で除算することによって、速度（間隔又はより小さな間隔群の間で測定した場合、一連の間隔及び一連の特定速度の両方又は一方にわたる標準速度か、平均速度か）が得られる。すなわち、一区間から次の区間へ漸増する特定速度の標準値によって、全画像を包含するより長い時間にわたる平均又は標準の速度が得られる。
【００５５】
図７は、流体流量のＴｉｍｅ−ＳＬＩＰ非造影撮像を概略的に描いている。選択パルス及び非選択パルスを、異なる大きさの長方形で表している。当業者には自明なように、異なるＢＢＴＩ間隔を使用することによって、心拡張期−心収縮期の異なる画像に対して以前の図で説明した異なる時間増分における連続する画像に類似した、動脈、静脈などの中の流体流量の連続する画像を得ることができる。明らかなように、動脈、静脈などに沿って移動された距離の測定、及び、関連する移動時間による除算もまた、この非造影流量撮像技術を使って、速度測定を行うことができる。
【００５６】
図８は、コンピュータプログラムコード構造の速度測定モジュールを概略的に描いている。このモジュールは、主要な又は監視の操作システムなどに関連するあらゆる好適な手段（例えば、操作者による、好適に表示されたアイコンまたメニュー上でのマウスクリック、タッチスクリーン上でのフィンガー操作入力、キーボード入力などの使用）によって実行できる。
【００５７】
ステップＳ０１で、スライスカウンタｎを１の値に初期設定する。ステップＳ０２で、（例えば、自動的に機械が実行する距離追跡アルゴリズムを使用するために、基本設定ダイアログなどに制御パラメータを設定することによる）距離測定の自動追跡システムのためのオプションが設定されたか否かを見るテストを行う。設定されている場合、次に、自動追跡機能が、ステップＳ０３で作動し、特定の速度ＳＶが、ステップＳ０４で、ｔｎ−１とその後のｔｎとにおいて取得されたスライス間で計算される。
【００５８】
自動追跡機能を望まない場合、次に、ステップＳ０５で、経路カウンタＰを１の値に初期設定する。その後、ステップＳ０６で、操作者が第１の位置を定義するために、待機ループが実行される。例えば、図９に示すように、操作者は、カーソルを（撮像時間ｔｎ−１に関連する）第１のポイントＰ１上に移動させて、所定のマウスボタンをクリックすることによって、動脈などの内の流体流量の表示画像上の第１の位置を定義できる。操作者が、特別に定義したポイントを初期経路位置、中間経路位置、又は最終経路位置であることを識別するために選択できるポップアップメニュー９３もまた設けることができる。例えば、図９に示すように、操作者定義のポイントＰ１〜Ｐ５が、（撮像時間ｔｎ−１とｔｎとの間の）蛇行曲線状の血管に沿った効果的なポイントに定義されているので、流量経路に沿った明確なポイントの間の直線部分は、スライス画像Ｓｎ−１からスライス画像Ｓｎまでの流量の距離に厳密に近似している。一旦第１の位置を（図８の）ステップＳ０６で定義すれば、経路カウンタＰは、図８に示すように、ステップＳ０７で増分される。その後、経路に沿った更なる位置を操作者が定義するにつれて、待機ループが実行されて、経路カウンタは、最終位置が定義されるまで（ステップＳ０８）増分される。そこで、ステップＳ０９で、経路パラメータＰ（ｍａｘ）に対する最大値を、経路カウンタのその時点の最新値として定義する。次に、ステップＳ１０で、特定速度ＳＶが、連続するスライス画像ＳｎとＳｎ−１との間の流量の増分として計算される。
【００５９】
次に、最終の時間ｔｎに関連する特定速度ＳＶ（どのように計算されようとも）がステップＳ１１で表示され、格納され、又は（例えば、システムの好ましい部分に望ましく構成されるような遠隔場所に）出力される。その後、移動平均速度ＭＶはステップＳ１２で計算されて、同様にステップＳ１３で表示され、格納され、又は（例えば、システム内に以前構成されたように）望ましいように出力される。ステップＳ１４ではスライスカウンタが、最終値（すなわち、スライスＳｎ（ｍａｘ））に達したか否かを見るテストを行う。到達していない場合は、次に、スライスカウンタはステップＳ１５で増分されて、制御は図８に示すように判断ステップＳ０２か、又は、直接ブロックステップＳ０６の経路カウンタの再初期設定に戻される（すなわち、操作者に、自動追跡システムを作動させたくない場合）。すべての有効なスライスデータが処理されると、このモジュールは図８に示すように終了する。
【００６０】
図１０に示すように、特定速度ＳＶ及び移動平均速度ＭＶの両方又は一方は、時間（又は、スライス数など）の関数としてグラフ表示できる。
【００６１】
［表示］
ところで、実施形態に係るＭＲＩシステムは、流体画像、及び、流体画像との差分処理の対象とされる基準画像の少なくとも一つに、流速に関する情報を含めて表示部に表示する表示制御部を更に備えてもよい。例えば、表示制御部は、ＭＲＩシステム制御部２２内に備えられ（図示を省略）、ＭＲＩシステム制御部２２内に備えられた表示制御部が、関連する構成要素を制御する。実施形態において、収集部によって収集された流体画像は、例えばＦＢＩ撮像法によって収集された流体画像やＴｉｍｅ−ＳＬＩＰ撮像法によって収集された流体画像であり、分解能の高い画像であるといえる。実施形態に係るＭＲＩシステムは、この分解能の高い形態画像自体から、機能情報としての「流体の速度」を算出するので、形態画像と、機能情報の算出根拠とが一致していることになる。
【００６２】
そこで、例えば、表示制御部は、形態画像を表示しながら、その形態画像上に、流体の速度の情報を重畳して表示してもよい。例えば、表示制御部は、下肢の血管を描出した形態画像上において、流速が存在する箇所には流速が存在することを示す表示（例えば、下肢の血管の塗り潰し）を行い、流速が存在しない箇所には流速が存在しないことを示す表示（例えば、下肢の血管の斜線塗り潰し）を行ってもよい。また、例えば、表示制御部は、３Ｄの形態画像を表示しながら、その３Ｄの形態画像上に、２Ｄで、流体の映像データを重畳して表示（例えば、図３に例示した画像１〜画像６を連続再生して表示部に表示）してもよい。
【００６３】
以上述べた少なくとも一つの実施形態の磁気共鳴イメージング装置及び磁気共鳴イメージング方法によれば、流体の流速を求めることが可能になる。
【００６４】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【００６５】
８ 心臓信号取得装置
９ 被検体
１０ 架台部
１１ 被検体テーブル
１２ 静磁場Ｂ0磁石
１３ 心電ゲート信号
１４ Ｇｘ、Ｇｙ、Ｇｚ傾斜磁場コイルセット
１６ ＲＦコイルアセンブリ
１８ 撮像ボリューム
２０ システム構成要素
２２ ＭＲＩシステム制御部
２４ 表示部
２６ キーボード／マウス
２８ プリンタ
３０ ＭＲＩシーケンス制御部
３２ Ｇｘ、Ｇｙ、Ｇｚ傾斜磁場コイルドライバ
３４ ＲＦ送信部
３６ 送信／受信スイッチ
３８ プログラムコード構造
４０ ＲＦ受信部
４２ ＭＲＩデータ処理部
４４ 画像再構成プログラムコード構造
４６ ＭＲ画像記憶部
５０ ＭＲＩシステムプログラム／データ格納部
【特許請求の範囲】
【請求項１】
被検体内を移動する流体の画像である流体画像を複数収集する収集部と、
前記複数の流体画像のうちの１つである基準画像と各流体画像との差分画像を用いて、前記流体の移動距離を特定する特定部と、
前記移動距離に対応する経過時間を、前記複数の流体画像の収集に用いられたパルスシーケンス情報から取得する取得部と、
前記移動距離を前記経過時間で除算することで前記流体の流速を算出する算出部と
を備えたことを特徴とする磁気共鳴イメージング装置。
【請求項２】
前記取得部は、前記流体画像の生成に用いるｋ空間の充填方法に応じて、前記移動距離に対応する経過時間を取得することを特徴とする請求項１に記載の磁気共鳴イメージング装置。
【請求項３】
前記取得部は、ｋ空間の中央から位相エンコードが配列されるセントリックオーダリングの場合には、前記経過時間として、前記パルスシーケンス情報からＢＢＴＩ（Black-Blood Time to Inversion）時間を取得し、ｋ空間に位相エンコードが順次配列されるシーケンシャルオーダリングの場合には、前記経過時間として、前記パルスシーケンス情報から、ＢＢＴＩ時間と位相エンコード数の半数分に相当する時間とを加算した値を取得することを特徴とする請求項２に記載の磁気共鳴イメージング装置。
【請求項４】
前記取得部は、各流体画像の移動距離毎に、前記経過時間を取得し、
前記算出部は、各移動距離を累積した累積異動距離を、各経過時間を累積した累積経過時間で除算することで、前記流速を算出することを特徴とする請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項５】
前記収集部は、心時相の異なる複数の流体画像を収集し、
前記特定部は、基準となる心時相にて収集された流体画像である前記基準画像と、所定の心時相にて収集された各流体画像との差分画像を用いて、心時相の異なる流体画像毎に前記流体の移動距離を特定し、
前記取得部は、各流体画像の移動距離毎に、前記経過時間を取得し、
前記算出部は、各移動距離を累積した累積異動距離を、各経過時間を累積した累積経過時間で除算することで、前記流速を算出することを特徴とする請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項６】
前記収集部は、被検体内を移動する流体のスピンに標識化を行い、所定時間経過後に前記スピンのエコー信号を受信する撮像を、前記所定時間を異ならせて複数回行い、該流体の画像であって時相の異なる複数の流体画像を収集し、
前記特定部は、前記標識化を行わない撮像により収集された画像である前記基準画像と、前記標識化を行う撮像により収集された各流体画像との差分画像を用いて、時相の異なる流体画像毎に前記流体の移動距離を特定し、
前記取得部は、各流体画像の移動距離毎に、前記経過時間を取得し、
前記算出部は、各移動距離を累積した累積異動距離を、各経過時間を累積した累積経過時間で除算することで、前記流速を算出することを特徴とする請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項７】
前記収集部は、Ｔｉｍｅ−ＳＬＩＰ（Time-Spatial Labeling Inversion Pulse）撮像法によって、時相の異なる複数の流体画像を収集することを特徴とする請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項８】
前記収集部は、ＦＡＳＥ（Fast Asymmetric Spin Echo）法を用いたＴｉｍｅ−ＳＬＩＰ撮像法によって、時相の異なる複数の流体画像を収集し、
前記取得部は、前記経過時間として、前記パルスシーケンス情報から、ＢＢＴＩ時間と実効エコー時間（TEeff（Effective Time to Echo））とを加算した値を取得することを特徴とする請求項１に記載の磁気共鳴イメージング装置。
【請求項９】
前記流体画像、及び、該流体画像との差分処理の対象とされる基準画像の少なくとも一つに、前記流速に関する情報を含めて表示部に表示する表示制御部を更に備えたことを特徴とする請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項１０】
前記特定部は、前記流体画像を表示部に表示し、前記流体が到達した位置の指定を受け付けることで前記移動距離を特定することを特徴とする請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項１１】
前記特定部は、前記画像を解析し、前記流体の移動経路を追跡することで前記移動距離を特定することを特徴とする請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項１２】
磁気共鳴イメージング装置で実行される磁気共鳴イメージング方法であって、
被検体内を移動する流体の画像である流体画像を複数収集する収集工程と、
前記複数の流体画像のうちの１つである基準画像と各流体画像との差分画像を用いて、前記流体の移動距離を特定する特定工程と、
前記移動距離に対応する経過時間を、前記複数の流体画像の収集に用いられたパルスシーケンス情報から取得する取得工程と、
前記移動距離を前記経過時間で除算することで前記流体の流速を算出する算出工程と
を含んだことを特徴とする磁気共鳴イメージング方法。
【請求項１】
被検体内を移動する流体の画像である流体画像を複数収集する収集部と、
前記複数の流体画像のうちの１つである基準画像と各流体画像との差分画像を用いて、前記流体の移動距離を特定する特定部と、
前記移動距離に対応する経過時間を、前記複数の流体画像の収集に用いられたパルスシーケンス情報から取得する取得部と、
前記移動距離を前記経過時間で除算することで前記流体の流速を算出する算出部と
を備えたことを特徴とする磁気共鳴イメージング装置。
【請求項２】
前記取得部は、前記流体画像の生成に用いるｋ空間の充填方法に応じて、前記移動距離に対応する経過時間を取得することを特徴とする請求項１に記載の磁気共鳴イメージング装置。
【請求項３】
前記取得部は、ｋ空間の中央から位相エンコードが配列されるセントリックオーダリングの場合には、前記経過時間として、前記パルスシーケンス情報からＢＢＴＩ（Black-Blood Time to Inversion）時間を取得し、ｋ空間に位相エンコードが順次配列されるシーケンシャルオーダリングの場合には、前記経過時間として、前記パルスシーケンス情報から、ＢＢＴＩ時間と位相エンコード数の半数分に相当する時間とを加算した値を取得することを特徴とする請求項２に記載の磁気共鳴イメージング装置。
【請求項４】
前記取得部は、各流体画像の移動距離毎に、前記経過時間を取得し、
前記算出部は、各移動距離を累積した累積異動距離を、各経過時間を累積した累積経過時間で除算することで、前記流速を算出することを特徴とする請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項５】
前記収集部は、心時相の異なる複数の流体画像を収集し、
前記特定部は、基準となる心時相にて収集された流体画像である前記基準画像と、所定の心時相にて収集された各流体画像との差分画像を用いて、心時相の異なる流体画像毎に前記流体の移動距離を特定し、
前記取得部は、各流体画像の移動距離毎に、前記経過時間を取得し、
前記算出部は、各移動距離を累積した累積異動距離を、各経過時間を累積した累積経過時間で除算することで、前記流速を算出することを特徴とする請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項６】
前記収集部は、被検体内を移動する流体のスピンに標識化を行い、所定時間経過後に前記スピンのエコー信号を受信する撮像を、前記所定時間を異ならせて複数回行い、該流体の画像であって時相の異なる複数の流体画像を収集し、
前記特定部は、前記標識化を行わない撮像により収集された画像である前記基準画像と、前記標識化を行う撮像により収集された各流体画像との差分画像を用いて、時相の異なる流体画像毎に前記流体の移動距離を特定し、
前記取得部は、各流体画像の移動距離毎に、前記経過時間を取得し、
前記算出部は、各移動距離を累積した累積異動距離を、各経過時間を累積した累積経過時間で除算することで、前記流速を算出することを特徴とする請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項７】
前記収集部は、Ｔｉｍｅ−ＳＬＩＰ（Time-Spatial Labeling Inversion Pulse）撮像法によって、時相の異なる複数の流体画像を収集することを特徴とする請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項８】
前記収集部は、ＦＡＳＥ（Fast Asymmetric Spin Echo）法を用いたＴｉｍｅ−ＳＬＩＰ撮像法によって、時相の異なる複数の流体画像を収集し、
前記取得部は、前記経過時間として、前記パルスシーケンス情報から、ＢＢＴＩ時間と実効エコー時間（TEeff（Effective Time to Echo））とを加算した値を取得することを特徴とする請求項１に記載の磁気共鳴イメージング装置。
【請求項９】
前記流体画像、及び、該流体画像との差分処理の対象とされる基準画像の少なくとも一つに、前記流速に関する情報を含めて表示部に表示する表示制御部を更に備えたことを特徴とする請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項１０】
前記特定部は、前記流体画像を表示部に表示し、前記流体が到達した位置の指定を受け付けることで前記移動距離を特定することを特徴とする請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項１１】
前記特定部は、前記画像を解析し、前記流体の移動経路を追跡することで前記移動距離を特定することを特徴とする請求項１又は２に記載の磁気共鳴イメージング装置。
【請求項１２】
磁気共鳴イメージング装置で実行される磁気共鳴イメージング方法であって、
被検体内を移動する流体の画像である流体画像を複数収集する収集工程と、
前記複数の流体画像のうちの１つである基準画像と各流体画像との差分画像を用いて、前記流体の移動距離を特定する特定工程と、
前記移動距離に対応する経過時間を、前記複数の流体画像の収集に用いられたパルスシーケンス情報から取得する取得工程と、
前記移動距離を前記経過時間で除算することで前記流体の流速を算出する算出工程と
を含んだことを特徴とする磁気共鳴イメージング方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【公開番号】特開２０１２−８１２７６（Ｐ２０１２−８１２７６Ａ）
【公開日】平成２４年４月２６日（２０１２．４．２６）
【国際特許分類】
【出願番号】特願２０１１−２２６１５９（Ｐ２０１１−２２６１５９）
【出願日】平成２３年１０月１３日（２０１１．１０．１３）
【出願人】（０００００３０７８）株式会社東芝 (54,554)
【出願人】（５９４１６４５４２）東芝メディカルシステムズ株式会社 (4,066)
【Ｆターム（参考）】
【公開日】平成２４年４月２６日（２０１２．４．２６）
【国際特許分類】
【出願日】平成２３年１０月１３日（２０１１．１０．１３）
【出願人】（０００００３０７８）株式会社東芝 (54,554)
【出願人】（５９４１６４５４２）東芝メディカルシステムズ株式会社 (4,066)
【Ｆターム（参考）】
[ Back to top ]