JP5942271B2 - 磁気共鳴イメージング装置及び流体強調画像取得法 - Google Patents

磁気共鳴イメージング装置及び流体強調画像取得法 Download PDF

Info

Publication number
JP5942271B2
JP5942271B2 JP2012553658A JP2012553658A JP5942271B2 JP 5942271 B2 JP5942271 B2 JP 5942271B2 JP 2012553658 A JP2012553658 A JP 2012553658A JP 2012553658 A JP2012553658 A JP 2012553658A JP 5942271 B2 JP5942271 B2 JP 5942271B2
Authority
JP
Japan
Prior art keywords
fluid
region
image
pulse
prepulse
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2012553658A
Other languages
English (en)
Other versions
JPWO2012098955A1 (ja
Inventor
甲亮 平井
甲亮 平井
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Publication of JPWO2012098955A1 publication Critical patent/JPWO2012098955A1/ja
Application granted granted Critical
Publication of JP5942271B2 publication Critical patent/JP5942271B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/4828—Resolving the MR signals of different chemical species, e.g. water-fat imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/5635—Angiography, e.g. contrast-enhanced angiography [CE-MRA] or time-of-flight angiography [TOF-MRA]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
    • G01R33/4833—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
    • G01R33/4836—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices using an RF pulse being spatially selective in more than one spatial dimension, e.g. a 2D pencil-beam excitation pulse
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5602—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by filtering or weighting based on different relaxation times within the sample, e.g. T1 weighting using an inversion pulse
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5613—Generating steady state signals, e.g. low flip angle sequences [FLASH]

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Vascular Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Signal Processing (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Description

本発明は、核磁気共鳴(以下「NMR」と言う)現象を利用して断層像撮影(以下「MRI」と言う)を行なう際に、所望の組織とその他の組織との間でコントラストを強調した画像を取得する技術に関する。
NMR現象を利用して断層像撮影を行なうMRI装置は、被検体、特に人体の組織を構成する原子核スピンが発生するNMR信号を計測し、その頭部、腹部、四肢等の形態や機能を2次元的に或いは3次元的に画像化する装置である。撮影においては、NMR信号には、傾斜磁場によって異なる位相エンコードが付与されるとともに周波数エンコードされて、時系列データとして計測される。計測されたNMR信号は、2次元又は3次元フーリエ変換されることにより画像に再構成される。
上記MRI装置を用いた、血液や髄液などの被検体内を流れる流体を他の組織(静止組織等)と異なるコントラストにして画像を取得する撮像法の一つとして、IR(Inversion Recovery)パルスを先行パルス(RFプリパルス)として用いる方法が知られている(例えば、特許文献1)。具体的には、流体が通過する上流側の領域である第1の領域にIRパルスを照射して、該第1の領域における流体の縦磁化を180度フリップ(励起)した後、180度フリップされた流体の縦磁化がT1緩和によって回復しnull状態に到達した時点で、該IRパルスを経験した流体が流れ込んだ下流側である第2の領域においてエコー信号の計測を行なう。この上流側で流体の縦磁化を所定の角度にフリップする操作をラベリングという。これにより、流体からのエコー信号のみを抑制し、第2の領域の画像(絶対値画像)において、該流体と他の組織(静止部)との間に画素値(絶対値)の差を生じさせることでコントラストを強調する。
特開2009-10113号公報
Lauenstein TC et al;Evaluation of optimized inversion-recovery fat-suppression techniques for T2-weighted abdominal MR Imaging : J Magn Reson Imaging 2008:27:1448-1454 J. Pauly, D. Nishimura; A K-Space Analysis of Smaii-Tip-angle Excitation:J. Magn. Reson., 81, 43-56 (1989)
しかしながら、IRパルスをRFプリパルスとして用いて、エコー信号の計測前に流体の縦磁化をnullにするためには、所定時間(TI)待つ必要があり、撮像時間の短縮が困難となる課題が残っていた。また、IRパルスを用いて取得される画像においては、流体と他の組織との間のコントラスト強調が充分ではないという課題が残っていた。
そこで、本発明の目的は、上記課題を鑑みてなされたものであり、RFプリパルスとしてIRパルスを用いても撮像時間を延長することなく、流体部と静止部との間のコントラストを強調した画像を取得することが可能なMRI装置及び流体強調画像取得法を提供することである。
上記課題を解決するために、本発明は、撮像領域よりも上流側の領域にRFプリパルスを印加して流体部の縦磁化を負方向に励起し、撮像領域からエコー信号を計測し、該エコー信号を用いて再構成した画像の位相情報に基づいて、流体部と静止部との間のコントラストを強調した画像を取得する。
具体的には、本発明のMRI装置は、所定のパルスシーケンスに基づいて、流体部を含む被検体の撮像領域からエコー信号の計測を制御する計測制御部と、エコー信号を用いて、流体部と静止部との間でコントラストを強調した画像を取得する演算処理部と、を有して成り、パルスシーケンスは、流体部の縦磁化を負方向に励起して該流体部をラベリングするRFプリパルスを備えるRFプリパルス部と、該ラベリングされた流体部が流れ込んだ撮像領域からエコー信号を計測する計測シーケンス部と、を有して成り、演算処理部は、画像の位相情報に基づいて、流体部と静止部との間のコントラストを強調した流体強調画像を取得することを特徴とする。
また、本発明の流体強調画像取得法は、所定のパルスシーケンスに基づいて、流体部を含む被検体の撮像領域からエコー信号の計測を制御する計測ステップと、エコー信号を用いて、流体部と静止部との間でコントラストを強調した画像を取得する演算処理ステップと、を有して成り、パルスシーケンスは、流体部の縦磁化を負方向に励起して該流体部をラベリングするRFプリパルスを備えるRFプリパルス部と、該ラベリングされた流体部が流れ込んだ撮像領域からエコー信号を計測する計測シーケンス部と、を有して成り、計測ステップでは、撮像領域よりも上流側の領域にRFプリパルスが印加されて、流体部の縦磁化がゼロに回復する前に、計測シーケンス部が実行され、演算処理ステップでは、画像の位相情報に基づいて、流体部と静止部との間のコントラストを強調した流体強調画像を取得することを特徴とする。
本発明のMRI装置及び流体強調画像取得法によれば、RFプリパルスとしてIRパルスを用いても、撮像時間を延長することなく、流体部と静止部との間のコントラストを強調した画像を取得することが可能となる。
本発明に係るMRI装置の一実施例の全体構成を示すブロック図。 IRパルスにより180度フリップされた縦磁化がT1回復してnull状態を得て正方向の状態になるまでの長い時間を待ち時間(TI)とする場合を示し、(a)図は、RFパルス(RF)の印加タイミングとエコー信号(signal)の発生タイミングとを示し、このパルスシーケンスの各タイミングに合わせて第1の領域でラベリングされた流体部と第2の領域における静止部の磁化の挙動とをそれぞれ示す。(b)図は、静止部の中をU字管の中を左側から右側に向けて水が流れる流体部を備えて成るファントム例である。(c)図は、(a)図のパルスシーケンスを用いて、(b)図に示すファントムを撮像した場合に得られる絶対値画像と位相画像の例である。 待ち時間(TI)を短くして、IRパルスにより180度フリップされた縦磁化がnull(ゼロ)状態以上にT1回復する前に、エコー信号を計測する場合を示す。(a)図は、RFパルス(RF)の印加タイミングとエコー信号(signal)の発生タイミングとを示し、このパルスシーケンスの各タイミングに合わせて第1の領域における流体部と第2の領域における静止部の磁化の挙動とをそれぞれ示す。(b)図は、(a)図のパルスシーケンスを用いて、図2(b)に示すファントムを撮像した場合に得られる絶対値画像と位相画像の例を示す。 公知のGMN(Gradient Moment Nulling)法にリフェーズ傾斜磁場パルスの例を占めす。(a)図が一次のリフェーズ傾斜磁場パルス波形の一例を示し、(b)図が二次のリフェーズ傾斜磁場パルス波形の一例を示す。 実施例1における、パルスシーケンスの一例を表すシーケンスチャートを示す。 実施例1における、演算処理部114が有する各機能の機能ブロック図を示す。 実施例1における、処理フローを表すフローチャートを示す。 実施例1における、IRパルスを印加するラベリング領域804とその下流側の撮像領域(FOV)803の設定例を示す。 実施例1における、処理フローの各ステップを実施して得られる結果の一例を示す。 実施例2における、下肢部における各IRパルスによるラベリング領域の設定例を示す。 実施例2における、Main-Scanシーケンスのシーケンスチャートを示す。 実施例2における、処理フローの各ステップを実施して得られる結果の一例を示す。 実施例3における、xy平面上の形状のみが特定されたz方向に柱状領域を選択的に励起する2次元空間選択励起の例を示す。 実施例3における、下肢部における左右の動脈801の各々に対する空間選択ラベリングを行う領域の設定例を示す。 実施例3における、RFプリパルス部のシーケンスチャートを示す。 実施例3における、処理フローの各ステップを実施して得られる結果の一例を示す。 実施例4における、下肢部における各IRパルスによるラベリング領域の設定例を示す。 実施例4における、RFプリパルス部における各IRパルスの印加タイミングと縦磁化の挙動を示す。
以下、添付図面に従って本発明のMRI装置の好ましい実施例について詳説する。なお、発明の実施例を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。
最初に、本発明に係るMRI装置を図1に基づいて説明する。図1は、本発明に係るMRI装置の一実施例の全体構成を示すブロック図である。
このMRI装置は、NMR現象を利用して被検体101の断層画像を得るもので、図1に示すように、静磁場発生磁石102と、傾斜磁場コイル103及び傾斜磁場電源109と、RF送信コイル104及びRF送信部110と、RF受信コイル105及び信号検出部106と、信号処理部107と、計測制御部111と、全体制御部108と、表示・操作部113と、被検体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)が印加されて、NMR信号(エコー信号)にそれぞれの方向の位置情報がエンコードされる。
RF送信コイル104は、被検体101にRFパルスを照射するコイルであり、RF送信部110に接続され高周波パルス電流が供給される。これにより、被検体101の生体組織を構成する原子のスピンにNMR現象が誘起される。具体的には、RF送信部110が、後述の計測制御部111からの命令に従って駆動されて、高周波パルスが振幅変調され、増幅された後に被検体101に近接して配置されたRF送信コイル104に供給されることにより、RFパルスが被検体101に照射される。
RF受信コイル105は、被検体101の生体組織を構成するスピンのNMR現象により放出されるエコー信号を受信するコイルであり、信号検出部106に接続されて受信したエコー信号が信号検出部106に送られる。
信号検出部106は、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の撮像領域についての画像の再構成に必要なエコーデータの収集を制御する。繰り返しの際には、2次元撮像の場合には位相エンコード傾斜磁場の印加量を、3次元撮像の場合には更にスライスエンコード傾斜磁場の印加量も、変えて行なう。位相エンコードの数は通常1枚の画像あたり128、256、512等の値が選ばれ、スライスエンコードの数は、通常16、32、64等の値が選ばれる。これらの制御により信号処理部107からのエコーデータを全体制御部108に出力する。
全体制御部108は、計測制御部111の制御、及び、各種データ処理と処理結果の表示及び保存等の制御を行うものであって、CPU及びメモリを内部に有する演算処理部114と、光ディスク、磁気ディスク等の記憶部115とを有して成る。具体的には、計測制御部111を制御してエコーデータの収集を実行させ、計測制御部111からのエコーデータが入力されると、演算処理部114がそのエコーデータに印加されたエンコード情報に基づいて、メモリ内のk空間に相当する領域に記憶させる。以下、エコーデータをk空間に配置する旨の記載は、エコーデータをメモリ内のk空間に相当する領域に記憶させることを意味する。
また、メモリ内のk空間に相当する領域に記憶されたエコーデータ群をk空間データともいう。そして演算処理部114は、このk空間データに対して信号処理やフーリエ変換による画像再構成等の処理を実行し、その結果である被検体101の画像を、後述の表示・操作部113に表示させると共に記憶部115に記録させる。
表示・操作部113は、再構成された被検体101の画像を表示する表示部と、MRI装置の各種制御情報や上記全体制御部108で行う処理の制御情報を入力するトラックボール又はマウス及びキーボード等の操作部と、から成る。この操作部は表示部に近接して配置され、操作者が表示部を見ながら操作部を介してインタラクティブにMRI装置の各種処理を制御する。
現在MRI装置の撮像対象核種は、臨床で普及しているものとしては、被検体の主たる構成物質である水素原子核(プロトン)である。プロトン密度の空間分布や、励起状態の緩和時間の空間分布に関する情報を画像化することで、人体頭部、腹部、四肢等の形態または、機能を2次元もしくは3次元的に撮像する。
(本発明に係る磁化とその位相の説明)
次に、本発明の基礎となる、RFプリパルスとしてIRパルスを用いて、異なる組織の横磁化に位相差を設定する原理について説明する。本発明では、任意のフリップ角で縦磁化をフリップした後に任意のT1回復時間を設けるRFパルスを全てIRパルスと称する。ただし、本発明のRFプリパルスは、IRパルスに限定されず、縦磁化を所望の角度にフリップ(励起)できる全てのRFパルスが可能である。
なお、以降の説明では、縦磁化の方向に関して、フリップ(即ち励起)される前の縦磁化方向を正方向(静磁場方向)としその反対方向を負方向(反静磁場方向)とする。この方向設定では、フリップされる前の縦磁化は正の方向を向く最大状態であり、90°より大きくフリップされた直後は、負の方向を向く状態となる。そして、縦磁化がフリップされて生成される横磁化の方向は、この縦磁化方向に垂直な方向となる。
従来技術は、前述したように、流体の上流でラベリングを行い、流体の縦磁化がnullに回復する時点でエコー信号の計測を行うことで、流体と他の組織(静止部)との間にエコー信号の強度差を生じさせ、このエコー信号の強度差を画像に反映させて、流体と他の静止部との間のコントラストを強調する。
これに対して、本発明は、RFプリパルスを用いて流体と他の組織との間に、横磁化の位相差を生じさせ、この位相差を利用して、流体と他の組織との間のコントラストを強調する。なお、従来技術のように、流体と他の組織との間にエコー信号の強度差を生じさせた上で、さらに本発明の位相差を生じさせて、流体と他の組織との間のコントラストをさらに強調しても良い。
そこで、本発明は、被検体に、上流側の第1の領域の縦磁化を負方向にフリップ(励起)するRFプリパルスを備えたRFプリパルス部と、該RFプリパルスによって励起された縦磁化がnull(ゼロ)以上に回復する前に下流側の第2の領域でエコー信号を計測する計測シーケンス部とを有して成るパルスシーケンスを用いて、被検体からエコー信号を計測する。縦磁化を負方向に励起するためには、縦磁化を90度より大きく270度未満にフリップすればよいので、RFプリパルスは、縦磁化をα(90<α<270)度に励起するIRパルスとする。一方、計測シーケンス部は、公知のパルスシーケンスを用いて良い。計測シーケンス部における横磁化を生成するRFパルスのフリップ角をβ(0＜β≦90)度として、RFプリパルス部におけるIRパルスのフリップ角α度と計測シーケンス部における横磁化生成のためのRFパルスのフリップ角β度とを決定する。好ましくは、α＝180、β＝90とする。
最初に、RFプリパルス部のIRパルスの印加から計測シーケンス部の横磁化生成RFパルス(例えばフリップ角90度)までの待ち時間(TI)の長短により、横磁化の位相をπ異ならせることができる(つまり位相の極性を異ならせる)ことを、図2,3を用いて説明する。
比較のために、待ち時間(TI)を長くしてコントラスト強調する場合を、図2を用いて説明する。図2は、IRパルスにより180度フリップされた縦磁化がT1回復してnull(ゼロ)状態を得て正方向の状態になるまでの長い時間を待ち時間(TI)とする場合を示す。
図2(a)は、RFパルス(RF)の印加タイミングとエコー信号(signal)の発生タイミングとを示し、このパルスシーケンスの各タイミングに合わせて第1の領域でラベリングされた流体部と第2の領域における静止部の磁化の挙動とをそれぞれ示す図である。
IRパルス201は、上流側である第1の領域に印加されるため、その第1の領域に存在する流体の縦磁化は180度フリップされて負方向の最大状態となる。つまり、このIRパルス201により流体の縦磁化はラベリングされる。他方、第1の領域と異なる下流側の第2の領域の静止部は、IRパルス201が印加されないため、その縦磁化は変わらず正方向の最大状態を維持する。
そして、時間の経過と共にラベリングされた流体は、その流速に応じて第1の領域から第2の領域に向けて移動していくと共に、その縦磁化は負方向の最大状態から正方向の状態に向けて指数関数的にT1回復していく。そして、RFパルス201の印加からある時間経過後には、流体は、その縦磁化がnull状態となる。そして、さらに時間を空けてRFパルス201の印加から待ち時間(TI)後には、流体は、その位置が第2の領域に移動し、その縦磁化が正方向の状態となるが、この状態は正方向の最大状態よりは小さい。この時点で、計測シーケンス部における横磁化生成のための90度RFパルス202が第2の領域に印加されると、第1の領域でラベリングされ第2の領域に移動した流体は、その縦磁化が正方向の小さい状態のまま再度90度フリップされて、小さな横磁化状態となる。他方、第2の領域の静止部の縦磁化は、正方向の最大状態から90度フリップされて、最大の横磁化状態となる。この結果、横磁化は静止部の方が流体部より大きいが、励起直後の横磁化は共に同じ方向を向くので、横磁化の位相は同位相となる。
したがって、このような横磁化状態において計測シーケンス部で計測されるエコー信号及び該エコー信号から得られる再構成画像においては、絶対値は静止部が流体部より大きくなるが、位相は静止部と流体部とで同位相となる。つまり、IRパルスによりフリップされて負の状態となった流体の縦磁化が正の状態に回復する程度に長い待ち時間(TI)の場合には、再構成画像において、流体部と静止部の画素値に関して絶対値の差は生じるが位相の差は生じない。よって、静止部と流体部とでは、絶対値の大小のみでコントラストをつけざるを得ず、コントラストが不十分となる可能性がある。なお、静止部と流体部間で信号強度のみで最もコントラストが大きくなる待ち時間(TI)は、180度フリップされてラベリングされた流体部の縦磁化がnull状態になるまでの時間である。
図2(c)は、図2(a)のパルスシーケンスを用いて、図2(b)に示すファントムを撮像した場合に得られる絶対値画像と位相画像の例を示している。図2(b)に示すファントムは、静止部213の中のU字管の中を左側から右側に向けて水が流れる流体部214を備えて成り、上流部である第1の領域211にIRパルス201が照射され、下流部である第2の領域212にRFパルス202が照射される。第1の領域211でIRパルス201によりラベリングされた流体が、第2の領域212のU字部左側領域に移動するので、絶対値画像においては、U字部左側領域221の信号強度が低減しているが、以降のU字部領域では、流体部の縦磁化がT1回復により正方向の最大状態に戻るため、静止部213と同じ信号強度となる。一方、静止部213と流体部214の横磁化の位相が同位相となるので、位相画像においては、静止部213と流体部214の位相が同じ一様な画像となる。
一方、待ち時間(TI)を短くしてコントラスト強調する場合を、図3を用いて説明する。図3は、待ち時間(TI)を短くして、IRパルスにより180度フリップされた縦磁化がnull状態以上にT1回復する前に、エコー信号を計測する場合を示す。ここで、短い待ち時間(TI)とは、RFプリパルスにより90度より大きくフリップされた縦磁化が負方向の状態を維持してnull状態以上にT1回復しない程度に短い時間を意味する。より詳細には、RFプリパルス部における流体部の縦磁化を励起して該流体部をラベリングするIRパルスと、計測シーケンス部における横磁化生成RFパルスと、の間の待ち時間(TI)を、該IRパルスによりラベリングされた流体部の縦磁化が負方向の状態を維持するような時間とする。
図3(a)は、図2(a)と同様に、RFパルス(RF)の印加タイミングとエコー信号(signal)の発生タイミングとを示し、このパルスシーケンスの各タイミングに合わせて第1の領域における流体部と第2の領域における静止部の磁化の挙動とをそれぞれ示す図である。
IRパルス301は、上流側である第1の領域に印加されるため、その第1の領域に存在する流体の縦磁化は180度フリップされて負方向の最大状態となる。つまり、このIRパルス301により流体の縦磁化はラベリングされる。他方、第1の領域と異なる第2の領域の静止部は、IRパルス301が印加されないため、その縦磁化は変わらず正方向の最大状態を維持する。
そして、時間の経過と共にラベリングされた流体部は、その流速に応じて第1の領域から第2の領域に向けて移動していくと共に、その縦磁化は負方向の最大状態から指数関数的にT1回復していく。そして、流体部の縦磁化が負方向の状態を維持してnull以上にT1回復せず、かつ、第1の領域でラベリングされた流体が第2の領域に行き渡るような短い待ち時間(TI)後に、計測シーケンス部における横磁化生成のための90度RFパルス302が第2の領域に印加される。
その結果、第2の領域の静止部の縦磁化は、正方向の最大状態から90度フリップされて、正方向(ここでは静止部の横磁化方向を正方向とする)の最大の横磁化状態となる。他方、第1の領域でラベリングされ第2の領域に移動した流体部は、その縦磁化が負方向の大きい状態のまま再度90度フリップされて、大きな負方向の横磁化状態となる。つまり、短い待ち時間(TI)後の90度RFパルス302直後は、静止部と流体部とで、横磁化の大きさは大差無いが、横磁化の方向が相互に反対方向を向くので、横磁化の位相はπ異なる(あるいは、位相の極性が異なる)ことになる。したがって、このような横磁化状態から得られるエコー信号及び該エコー信号から得られる再構成画像においては、静止部と流体部とで、絶対値は大差無いが、位相はπ異なる(あるいは位相極性が異なる)。
図3(b)は、図3(a)のパルスシーケンスを用いて、図2(b)に示すファントムを撮像した場合に得られる絶対値画像と位相画像の例を示している。上流部である第1の領域211にIRパルス301が照射され、下流部である第2の領域212にRFパルス302が照射される。第1の領域211でIRパルス301によりラベリングされた流体が、第2の領域212のU字部左側領域に移動するが、待ち時間(TI)が短いために、絶対値画像321においては、U字部左側領域を含めて、静止部213と流体部214の信号強度差は少ない。他方、静止部213の横磁化の位相と第1の領域211でIRパルス301によりラベリングされた流体の横磁化の位相とが逆位相となる(つまりπ異なる)ので、位相画像322においては、静止部213の位相と流体部214のU字部左側領域325の位相が大きく異なる(つまりπ異なる)画像となる。
そこで、絶対値画像321の各画素値を位相画像322の対応する画素値(つまり位相値)に基づいて重み付けする演算処理を行うことによって、演算処理後の画像323においては、第1の領域211でIRパルス301によりラベリングされた流体が流れ込んだU字部左側領域326を他の静止部及び流体部に対してコントラストを強調した画像を得ることができる。
以上のことから、本発明は、IRパルスにより流体部をラベリングした後の待ち時間(TI)を短くして計測されたエコー信号を用いて再構成した被検体の画像に対して、該画像の位相情報に基づいて、ラベリングした流体部を他方の静止部に対して強調するコントラスト強調処理を施してコントラスト強調画像を取得する。具体的には、再構成した被検体の複素画像におけるラベリングした流体部と他の静止部との間の位相差を利用して、複素画像の絶対値をとった絶対値画像を重み付ける。これにより、待ち時間(TI)を長くして画素値の絶対値のみによるコントラスト強調の場合よりも、撮像時間を短くして、且つ、ラベリングした流体部と他の静止部との間のコントラストをさらに強調する。
(他の要因の位相誤差の除去について)
一般的に、複素画像には、RFプリパルス後の短い待ち時間(TI)によって付与された位相差π(反対位相極性)以外に撮像によって発生する位相誤差が混入するので、この位相誤差を除去する必要がある。
この位相誤差には、静磁場不均一やケミカルシフトなどの共鳴周波数ズレによって、画像用エコー信号の計測中に蓄積される位相誤差、A/Dに対する傾斜磁場印加タイミングの遅れ等のハードウェアの不完全性に起因した位相誤差、更には、被検体の動きに起因した位相誤差、が含まれる。
共鳴周波数のズレに起因して時間的に蓄積する位相誤差は、スピンエコー系のシーケンスでは無視できる。なぜなら、スピンエコー系のシーケンスでは、90度RFパルスによる励起からエコー時間(TE)までの間に180度再収束RFパルスを用いるので、位相誤差がキャンセルされるため、時間的に蓄積する位相誤差を無視できる。一方、グラジエントエコー系のシーケンスでは、180度再収束RFパルスが無いために、時間的に蓄積する位相誤差を無視できない。
そのため、RFプリパルスを加えない場合の位相画像(リファレンス位相画像)を予め事前計測(プリスキャン；Pre-Scan)にて撮像して求めておき、RFプリパルスを用いた場合の位相画像からリファレンス位相画像を差分処理して位相差画像を得ることで、時間的に蓄積する位相誤差を除去できる。また、プリスキャンにより得られるリファレンス位相画像にはハードウェアの不完全性に起因した位相誤差も含まれる。
即ち、リファレンス位相画像には、共鳴周波数のズレに起因して時間的に蓄積する位相誤差とハードウェアの不完全性に起因した位相誤差とが含まれるので、位相差画像から、これら二種の位相誤差が除去される。なお、これら二種の位相誤差は、空間的な位相変化が緩やかであるため、リファレンス位相画像は低空間分解能であっても十分な精度でこられ二種の位相誤差を表すことになる。そのため、リファレンス位相画像取得のためのプリスキャンは、撮像時間の短い低空間分解能(例えば32*32マトリクス程度)撮像で十分である。
また、エコー時間(TE)の異なる2つ以上のエコー信号を連続的に取得するマルチエコーシーケンスを用いることで、エコー信号間の時間差と位相差とから周波数ズレを算出し、その周波数ズレから目的のエコー時間(TE)における位相誤差を算出して除去することも可能である。
また、血流など、被検体自身の動き又はその内部の動き(等速運動や加速度運動)に起因した位相誤差に関しては、公知のGMN(Gradient Moment Nulling)法に基づく1次以上のリフェーズ傾斜磁場パルスをパルスシーケンスに加えることで、動きの影響を除去することが可能である。リフェーズ傾斜磁場パルスの一例を図4に示す。等速度運動(1次)による位相誤差を抑制するためには、図4(a)のような3つの傾斜磁場パルスの構成で、強度(絶対値)一定で面積比が1:-2:1の比率になるような傾斜磁場パルス波形を等速度運動方向に印加する。また、加速度運動(2次)による位相誤差を抑制するためには、図4(b)のような4つの傾斜磁場パルスの構成で、強度一定で面積比が1:-3:3:-1の比率になるような傾斜磁場パルス波形を加速度運動方向に印加する。
上記のようにプリスキャンによる位相計測、マルチエコー計測、1次以上のリフェーズ傾斜磁場パルスを組み合わせることで、位相差画像から各種位相誤差を除去することができるので、RFプリパルス後の短い待ち時間(TI)によって生じた、該RFプリパルスによってラベリングされた流体部と他の静止部との間の横磁化の位相差のみを抽出した位相差画像を取得することができる。そして、その位相差を用いて、ラベリングされた流体部と他の静止部との間のコントラストを強調した画像を取得することが可能となる。
次に、本発明のMRI装置及び流体強調画像取得法の実施例1を説明する。本実施例は、RFプリパルスにより流体部をラベリングした後に待ち時間(TI)を短くして画像を取得し、該画像の位相情報を用いて、ラベリングされた流体部と他の静止部との間のコントラストを強調したコントラスト強調画像を得る。具体的には、ラベリングされた流体部と静止部とで横磁化の位相をπ異ならせてエコー信号を計測し、計測されたエコー信号を用いて再構成された画像の位相画像に基づいて、該画像の画素毎の重み係数を決定し、該重み係数を絶対値画像に画素毎に掛け合わせて流体強調画像を取得する。以下、図5〜9に基づいて本実施例を詳細に説明する。
(実施例1のパルスシーケンス)
最初に、図5を用いて本実施例のパルスシーケンスを説明する。図5は、本実施例のパルスシーケンスの一例を表すシーケンスチャートであり、RF,Gs,Gp,Gr,A/D,及びSignalは、それぞれ、RFパルス、スライス傾斜磁場、位相エンコード傾斜磁場、読み出し傾斜磁場、エコー信号のサンプリング期間、及びエコー信号を示す(後述する他の実施例においても同様である)。図5(a)は、画像用エコー信号を計測するためのファーストスピンエコー(FSE;Fast-Spin Echo)シーケンスを用いる計測シーケンス部551の前に、RFプリパルスとしてのIRパルスを印加するRFプリパルス部550を加えたメインスキャン(Main-Scan)シーケンスの一例を示す。図5(b)は、図5(a)からRFプリパルス部550を除き、計測シーケンス部551におけるスライスと位相エンコードの各傾斜磁場パルスの変化量を大きくすることで、低空間分解能撮像に対応したプリスキャン(Pre-Scan)シーケンス560の一例を示す。なお、本実施例のベースとなるパルスシーケンスは、FSEシーケンスに限定されるものではなく、他のパルスシーケンスでもよい。これらMain-ScanシーケンスとPre-Scanシーケンスの各RFパルスと各傾斜磁場パルスの印加及びエコー信号の計測は計測制御部111により制御される。
最初に、RFプリパルス部550と計測シーケンス部551とを有して成るMain-Scanシーケンスの一例を図5(a)に基づいて説明する。このMain-Scanシーケンスにおいては、RFプリパルス部550のIRパルス501印加から計測シーケンス部551の横磁化生成用RFパルス504までの待ち時間(TI)を短くして、横磁化生成用RFパルス504によって生成される横磁化の位相が、IRパルス501によってラベリングされた流体部と他の静止部とでπ異なるようにする。
RFプリパルス部550は、IRパルス501とスライス傾斜磁場パルス502とを同時に印加し、その後にスポイル傾斜磁場パルス(503-1〜503-3)を印加する。IRパルス501とスライス傾斜磁場パルス502とで、所望の領域の縦磁化が選択的に180度反転される。本実施例では、所望の領域として、画像化する撮像領域(FOV)に流れ込む流体の上流側の領域とする。これにより、上流側の領域における流体部の縦磁化を180度フリップしてラベリングすることができる。そして、このIRパルス501の後に、スライス方向(Gs)、位相エンコード方向(Gp)、読み出し方向(Gr)の少なくとも1軸方向に、好ましくは3軸方向にスポイル傾斜磁場パルス(503-1〜503-3)を印加して、該IRパルス501で180度未満に励起されることで生じた横磁化を消失させる。
計測シーケンス部551は、FSEシーケンスに基づいてエコー信号の計測を行なう。RFプリパルス部550でラベリング対象となった上流側の領域よりも下流側の撮像領域(FOV)において、該撮像領域に流れ込んだラベリングされた流体部の縦磁化と他の静止部の縦磁化を共に90度フリップする90度パルス504と同時にスライス傾斜磁場パルス505を印加する。そして、動きによる影響を補正するために、印加強度の比率が1:-1:1、印加時間の比率が1:2:1となるような1次のリフェーズ傾斜磁場パルス(506、507)をスライス方向に印加する。次に180度再収束パルス511-1と同時にスライス傾斜磁場パルス512-1を印加し、その前後にも印加時間がスライス傾斜磁場パルス512-1の1/6で、強度(絶対値)が同じになるようなスライス方向のリフェーズ傾斜磁場パルス(509-1、513-1)を印加する。
次の180度再収束パルス511-2とスライス傾斜磁場(509-2,512-2,513-2)についても同様である。180度再収束パルス511-1の中心前後で、横磁化が感じる傾斜磁場極性は反転するため、509、512-1、513の傾斜磁場パルスの印加面積の比率は1:-3:3:-1となり、2次のリフェーズ傾斜磁場パルスとすることができる。読み出し方向(Gr)にも2次のリフェーズ用の傾斜磁場パルス(508、510、516-1)と読み出し傾斜磁場パルス517-1を印加する。
読み出し傾斜磁場パルス517-1の中心で、エコー信号519-1のピークを検出するため、508,510,516-1と517-1の中心までの傾斜磁場パルスを1つの傾斜磁場パルスのユニットとすると、印加時間を同一、傾斜磁場強度比が1:-3:-3:1となるように印加する。スライス傾斜磁場と同様、横磁化は180度再収束パルス511の影響で、強度比を1:-3:3:-1として感じるため、2次のリフェーズ効果を得ることができる。
また、読み出し方向(Gr)のリフェーズ傾斜磁場パルス516-1のタイミングで、スライス方向(Gs)にスライスエンコード傾斜磁場パルス514と、位相エンコード方向(Gp)に位相エンコード傾斜磁場パルス515を、それぞれ印加する。そして、読み出し傾斜磁場パルス517印加後は、スライス方向(Gs)と位相エンコード方向(Gp)にリワインド傾斜磁場パルス(520、521)を印加する。514、515、520、521の傾斜磁場パルスは180度再収束パルス毎に変化するように制御することで各種エンコードが実施される。
また、読み出し傾斜磁場パルス517-1印加時には、A/D518-1することでエコー信号519-1を計測する。読み出し方向(Gr)には、読み出し傾斜磁場パルス印加後に516-1と同じ形のリフェーズ傾斜磁場パルス522-1を印加し、次の180度再収束パルス511-2前の読み出し傾斜磁場パルス517-1の右半分とリフェーズ傾斜磁場パルス522-1、及び、後に再度繰り返されるリフェーズ傾斜磁場パルス516-2と読み出し傾斜磁場パルス517-2の左半分の傾斜磁場とを1つの傾斜磁場パルスのユニットとすると、横磁化が感じる傾斜磁場面積比は1:-3:3:-1となるため、2次のリフェーズが繰り返されることになる。
次に、計測シーケンス部551のみを有するPre-Scanシーケンスを図5(b)に基づいて説明する。図5(b)は、図5(a)からRFプリパルス部550を除き、計測シーケンス部551におけるスライス・位相エンコード傾斜磁場パルス(531,532,533,534)の変化量を大きくすることで、低空間分解能撮像に対応したPre-Scanシーケンスの一例である。他は、図5(a)のMain-Scanシーケンスと同じなので、詳細な説明を省略する。このPre-Scanシーケンスにより計測されたエコー信号を用いて位相画像を取得することにより、前述したように、RFプリパルスとしてのIRパルス501後の短い待ち時間(TI)によって生じた、ラベリングされた流体部と他の静止部の横磁化間の位相差以外の色々な位相誤差を纏めて取得することが可能となる。
(撮像条件)
次に、本実施例の好適な撮像条件について説明する。
最初に待ち時間(TI)の条件について説明する。待ち時間(TI)は、目的の流速(Vb)と流れ方向の撮像領域幅(FOVb)とに応じて設定する必要がある。具体的には、(1)式に基づいて設定する。
TI＝FOVb／Vb (1)
例えば、FOVb＝200mmで流速Vb＝500mm/s程度の血流をターゲットにした場合には、(1)式よりTI＝400msと求まる。目的の流速(Vb)は、例えば、予めPhase Contrastシーケンス等で流速を測定しておけば、TI設定は一意に決めることができる。
また、RFプリパルスによって反転させた流体部の縦磁化がnullになる前までに計測シーケンス部を実行する必要があるため、待ち時間(TI)はその上限値(limit TI)以下とする。limit TIは、流体のT1値(一定値)とIRパルスのフリップ角φによって一意的に定まり、φ＝180度の場合に最大となる。つまり、
TI＜＝limit TI(φ)＜＝limit TI(φ＝180度) (2)
とする必要がある。
次に、流れ方向の撮像領域幅(FOVb)の条件について説明する。(1)式と(2)式とから
FOVb＜limit TI×Vb (3)
となるので、流れ方向の撮像領域幅(FOVb)は、その上限値をlimit FOVbとすると
FOVb＜limit FOVb＝limit TI×Vb (4)
とする必要がある。
そこで、操作者により設定入力された流れ方向の撮像領域幅(FOVb)から(1)式に基づいて計算した待ち時間(TI)がlimit TIを越える場合には、操作者に対して、流れ方向の撮像領域幅(FOVb)をlimit FOVb以下とするように報知する。
なお、流れ方向のラベリング領域幅(つまりIRパルスの印加幅)は、流れ方向の撮像領域幅(FOVb)以上であればよい。
また、計測シーケンス部で計測するエコーデータの計測順序については、コントラストに寄与するｋ空間の低周波領域データを最初に取得するセントリックオーダーが好ましい。公知のセグメント計測法を適用する場合は、TI時間待って流入したラベリングされた磁化がFOV内から完全に流出する前に複数のセグメントデータを取りきる必要がある。そこで、ｋ空間の分割数であり、RFプリパルスを印加した後に、連続して収集するエコー数でもある、セグメント数(ファーストスピンエコーシーケンスの場合はエコートレイン数(ETL; Echo Train Length))に対して上限値(limit ETL)を設定するのがよい。ETLと180度RFパルスの間隔であるES(Echo Space)から算出される、計測シーケンス部の最短繰り返し時間minTRは式(5)で表され、流出する前にデータを取りきるための限界値limit max TRは、ラベリング幅ｄと流速Vbを使って式(6)のように表される。
min TR ＝ ETL * ES (5)
limit max TR ＝ d/Vb (6)
式(5)(6)を使うとlimit ETLは、式(7)のように算出できる。
Limit ETL ＝ d/Vb /ES (7)
例えば、ラベリングを行うIRパルスの印加幅d＝200mmでES＝10msecの撮像条件で、かつ流速Vb＝500mm/secの場合には、セグメント数の上限数はLimit ETL数＝40となり、設定できるエコートレイン数はこの上限数以下とするのが望ましい。設定画面上で撮像条件設定する際にはLimit ETL数を提示してもよい。
(本実施例の機能処理部の説明)
次に、本実施例の演算処理部114が有する、本実施例に係る各演算処理機能を図6に基づいて説明する。図6は、本実施例の演算処理部114が有する各機能の機能ブロック図である。本実施例に係る各演算処理機能は、シーケンス実行部601と、画像再構成部602と、位相画像演算部603と、位相差画像演算部604と、マスク処理部605と、位相アンラップ処理部606と、コントラスト強調処理部607と、リージョン・グローイング処理部608と、を有して成る。なお、これらの各演算処理機能の構成は、以降に説明する他の実施例おいても同様であるが、他の実施例では、一部の演算処理機能の処理内容が実施例に応じて一部異なり、異なる部分については、各実施例において説明する。
シーケンス実行部601は、操作者が設定入力した撮像条件の確認と修正処理を行って、撮像条件を確定し、確定した撮像条件に基づくPre-ScanシーケンスとMain-Scanシーケンスを計測制御部111に実行させる。
画像再構成部602は、Pre-ScanシーケンスとMain-Scanシーケンスとでそれぞれ計測されたエコー信号のデータ(エコーデータ)に対してフーリエ変換を施して、複素画像をそれぞれ再構成する。また、複素画像の各画素値の絶対値を演算して絶対値画像を得る。
位相画像演算部603は、複素画像の画素毎にその画素値である複素数の位相(偏角)を演算し、位相画像を得る。
位相差画像演算部604は、2つの位相画像を画素毎に差分演算して、位相差画像を得る。
マスク処理部605は、入力画像の画素毎にその画素値と所定閾値とを比較演算し、画素値を所定範囲の値(例えば、0〜1の値)に変換して、マスク画像を作成する。また、作成したマスク画像を他の画像に施して、つまり、画素毎に掛け合わせるマスク処理を行い、マスク処理後の画像を得る。
位相アンラップ処理部606は、入力された位相画像の各画素値において主値周りを除去する位相アンラップ処理を行、アンラップ処理後の位相画像を得る。
コントラスト強調処理部607は、位相差画像が有する位相情報に基づいて絶対値画像に重み付け演算を施してコントラスト強調処理を行う。具体的には、位相差画像の各画素の画素値(位相差)に基づいてその画素の重み係数を決定し、決定した重み係数を絶対値画像の対応する画素の画素値に掛け合わせてその画素値を重み付ける。この位相差画像に基づいた重み付け処理がコントラスト強調処理であり、コントラスト強調処理後の画像がコントラスト強調画像となる。
リージョン・グローイング処理部608は、位相画像上で、隣り合うピクセル間で位相の連続性が高い位相値は同一位相として処理する。例えば、隣り合うピクセル間で位相差が所定の閾値以内であれば、位相の連続性が高いと判断して同一位相とする。具体的には、最初に、所望の組織に対して操作者により手動で設定された開始点(Seed点)の位相値と該開始点に隣接するピクセルの位相値との間の位相差と閾値とを比較して、位相差が閾値以内であれば同一位相と判定し、開始点の位相値と同じ同一位相値領域とする。そして、同一位相値領域に隣接するピクセルの位相値と同一位相値領域の位相値との差が閾値以内であれば同一位相として、その隣接ピクセルも同一位相値領域とする。この処理を繰り返して、同一位相値領域を拡張していく。
以下、上記各機能部が連携して行なう本実施例の処理フローの具体的な説明を通して、これらの各機能部の具体的処理を説明する。
(本発明の処理フロー)
次に、図7を用いて本実施例の処理フローを説明する。図7は、本実施例の処理フローを表すフローチャートである。本処理フローは、予めプログラムとして記憶部115に記憶されており、演算処理部114が記憶部115からそのプログラムを読み込んで実行することにより実施される。また、本処理フローの説明では下肢部を想定し、図8に、上流領域であって、RFプリパルス部550でラベリングを行うためのIRパルスを印加するラベリング領域804とその下流側である撮像領域(FOV)803の設定例を示す。この2つの領域には動脈(太い実線)801に沿って静脈(太い点線)802も流れているが、本実施例では動脈801のみをラベリング対象とする。そして、図9に、図8に示す各領域に対して、図7に示す処理フローの各ステップを実施して得られる結果の一例を示す。以下、各ステップの処理の詳細を説明する。
ステップ701で、シーケンス実行部601は、図8に示すような位置決め画像を表示して、血流の上流側のラベリング領域804と、下流側の撮像領域803の設定入力を受け付ける。そして、シーケンス実行部601は、操作者が設定入力した他の撮像条件と合わせて、撮像可能か否かをチェックし、撮像不可であればその旨を操作者に報知すると共に、ラベリング領域804と撮像領域803の位置及び幅の修正入力や、他の撮像条件についての撮像可能な条件を操作者に提示して該他の撮像条件の修正入力を受け付ける。そして最終的に撮像可能な撮像条件を確定し、この確定した撮像条件に基づいて、図5(a)に示したMain-Scanシーケンスと図5(b)に示したPre-Scanシーケンスの実行に必要な各種制御データを具体的に計算する。
ステップ702で、シーケンス実行部601は、ステップ701で計算したPre-Scanシーケンスの各種制御データを計測制御部111に通知して該Pre-Scanシーケンスを実行させる。計測制御部111は、その指示を受けて、Pre-Scanシーケンスを実行してエコー信号の計測を制御し、計測したエコー信号のデータ(エコーデータ)を演算処理部114に通知する。画像再構成部602は、エコーデータをフーリエ変換して低空間分解能の複素画像を得る。そして、位相画像演算部603は、得られた複素画像からその低空間分解能の位相画像(第1の位相画像)901を求める。この第1の位相画像901は、前述したように、IRパルス501によってラベリングされた流体部と他の静止部との間の位相差以外の色々な位相誤差を纏めて含む。
ステップ703で、シーケンス実行部601は、ステップ701で計算したMain-Scanシーケンスの各種制御データを計測制御部111に通知して該Main-Scanシーケンスを実行させる。計測制御部111は、その指示を受けて、Main-Scanシーケンスを実行してエコー信号の計測を制御し、計測したエコー信号のデータ(エコーデータ)を演算処理部114に通知する。画像再構成部602は、そのエコーデータをフーリエ変換して複素画像及びその絶対値画像906を得る。そして、位相画像演算部603は、得られた複素画像からその位相画像(第2の位相画像)902を求める。
ステップ704で、位相差画像演算部604は、ステップ702で得られた第1の位相画像901を、ステップ703で得られた第2の位相画像902と同じ空間分解能の位相画像に変換した後に、第2の位相画像902との差分処理921を行い、位相差画像903を得る。この位相差画像903は、共鳴周波数のずれに起因した位相誤差とハードウェアの不完全性に起因した位相誤差とが除去されて、IRパルス501後の短い待ち時間(TI)によって生じた位相差のみが反映された位相画像となる。
ステップ705で、マスク処理部605は、ステップ703で得られた絶対値画像906の各画素の画素値(絶対値)に対して閾値(例えば、各画素値の内の最大値の20%)を設定して、その閾値より小さい画素値を持つ画素を背景(ノイズ領域)として除外することで、絶対値画像906における被検体領域のみを抽出するための第1のマスク画像908を作成する。具体的には、閾値より小さい画素値を持つ画素には0を、閾値より大きい画素値を持つ画素には1を、それぞれ割り当てて第1のマスク画像908を作成する。
ステップ706で、マスク処理部605はステップ704で得られた位相差画像903に、ステップ705で作成した第1のマスク画像908を施して、つまり、位相差画像903に第1のマスク画像908を画素毎に掛け合わせるマスク処理922を行い、位相差画像903から背景領域を除外し被検体領域のみを抽出した位相差画像904とする。除外した背景領域の画素値(位相値)には、所定の一定値(例えば0)を割り当てる。なお、第1のマスク画像908の背景領域の値は0なので、画素毎に掛け合わせれば結果画像の背景領域の値は必然的に0になる。
ステップ707で、位相アンラップ処理部606は、ステップ706でマスク処理された位相差画像904に対して、主値周りを除去する位相アンラップ処理を行う。さらに、静止部の位相値を基準位相θrefとして、全画素の位相値θから基準位相θrefとの差分(θ−θref)をとることで、つまり位相差画像の各画素値から基準位相を一様に引いた修正位相差画像を作成する。修正位相差画像は、静止部の位相値からの差分位相を表す画像であり、静止部の位相はゼロ、ラベリングされた血流部の位相はπとなる。
ステップ708で、リージョン・グローイング処理部608は、ステップ707で取得された修正位相差画像に対して、リージョン・グローイング処理を実施する。操作者により手動で、所望の組織に対して開始点(Seed点)を設定してもらい、この開始点を起点としてリージョン・グローイング処理を実施し、開始点の位相値と同一の位相値と判定できるピクセルを纏めて同一位相値領域とする。リージョン・グローイング処理の結果、修正位相差画像において同一位相値領域の各ピクセルの位相値は、開始点の位相値に置き換えられる。なお、本ステップ708の処理は省略しても良い。
ステップ709で、コントラスト強調処理部607は、ステップ707又は708で得られた修正位相差画像の各画素の画素値(位相差)に基づいてその画素の重み係数を決定し、決定した重み係数の分布を表す第2のマスク画像905を作成する。具体的には、ステップ707で得られた修正位相差画像の各画素の画素値に対して、所定の閾値(例えば±π/2)を設定して、画素値である位相値θの絶対値がその閾値未満の場合(つまり、−π/2＜θ＜＋π/2)に1、それ以外の場合(つまり、[θ≦−π/2] or [＋π/2≦θ])に[0〜1]の値に変換してその画素の重み係数とする。この変換により、ラベリングされた血流部の位相が[0〜1](例えば0.5)の重み係数に、静止部の位相が1の重み係数に、それぞれ変換されることになる。修正位相差画像の全画素に対して同様に重み係数を決定し、各画素の重み係数分布を表す第2のマスク画像905を作成する。この第2のマスク画像905がコントラスト強調用マスク画像となる。
ステップ710で、コントラスト強調処理部607は、ステップ709で得た第2のマスク画像(コントラスト強調用マスク画像)905をステップ703で得られた絶対値画像906に施す(923)。具体的には、絶対値画像906と第2のマスク画像905とを同一画素毎に画素値同士を掛け合わせる(923)ことで、絶対値画像906の各画素の画素値を、第2のマスク画像905の画素値で重み付け処理を行う。この第2のマスク画像905を用いた、即ち位相差画像903に基づいた重み付け処理(923)がコントラスト強調処理であり、このコントラスト強調処理によりコントラスト強調画像910を得る。
コントラスト強調画像910においては、静止領域に対してラベリングされた血流(動脈)領域が抑制された画像となる。即ち、絶対値画像906において、ラベリングされた血流(動脈)領域と他の静止領域との間のコントラストが強調された画像となる。図9に示すコントラスト強調画像910の例では、ラベリングされた血流(動脈)領域の信号が抑制されて、他の静止領域のみの輝度が強調された画像となっていることが理解される。
以上までが、本実施例のコントラスト強調画像取得方法の処理フローの説明である。
以上の構成により、本実施例のMRI装置及び流体強調画像取得法は、RFプリパルスの印加から計測シーケンス実行までの待ち時間(TI)を短く設定するので、撮像時間を短縮できる。さらに、ラベリングする血流部と他の静止部とでπの位相差を設定して位相差画像を得て、該位相差画像に基づいて、絶対値画像を重み付けすることで、待ち時間を長くして信号強度差のみでコントラストをつける手法と比較して、ラベリングした血流部と他の静止部間のコントラストが更に強調された画像を取得することができる。
また、前述のステップ709では、静止部の信号を血流部の信号に対して抑制するように重み係数を決定したが、逆に、静止部の信号を血流部の信号に対して抑制するように重み係数を決定してもよい。具体的には、修正位相差画像の各画素の画素値(位相値)θの絶対値が閾値未満の場合(つまり、−π/2＜θ＜＋π/2)に[0〜1]、それ以外の場合(つまり、[θ≦−π/2] or [＋π/2≦θ])に1の値に変換してその画素の重み係数としてもよい。
以上説明したように、本実施例のMRI装置及び流体強調画像取得法は、流体部の縦磁化を負方向に励起して該流体部をラベリングするRFプリパルス(IRパルス)を備えるRFプリパルス部と、該ラベリングされた流体部が流れ込んだ撮像領域からエコー信号を計測する計測シーケンス部と、を有して成るパルスシーケンスを用いて、撮像領域よりも上流側の領域にRFプリパルスを印加して流体部の縦磁化を負方向に励起してラベリングし、ラベリングされた流体部の縦磁化がnullに回復する前に、撮像領域からエコー信号を計測し、該エコー信号を用いて再構成した画像の位相情報に基づいて、流体部を静止部に対してコントラストを強調した画像を取得する。その結果、RFプリパルスとしてIRパルスを用いても、撮像時間を短縮しつつ、血流部と静止部との間のコントラストを強調した画像を取得することが可能となる。
次に、本発明のMRI装置及び流体強調画像取得法の実施例2について説明する。本実施例は、RFプリパルス部がフリップ角の異なる2つのRFプリパルスを有し、各RFプリパルスを、撮影領域(FOV)を間に挟んで流体部の上流側と下流側の領域に、それぞれ照射することで、流れ方向の異なる二つの流体部(例えば動脈と静脈)に対して異なるラベリングを行う。これにより、流体部と静止部間のみならず、2つの流体部間で異なるコントラストの画像を得る。2つのRFプリパルスは、例えばフリップ角の異なるIRパルスであって、フリップ角が異なることにより、流れ方向の異なる二つの流体部間でコントラスト強調を異ならせる。以下、RFプリパルスをIRパルスとして、流れ方向の異なる二つの流体部を動脈と静脈として、図10〜12を用いて本実施例を詳細に説明する。
最初に、本実施例2のRFプリパルス部について説明する。RFプリパルス部の2つのIRパルスの内の一方のIRパルス(以下IR1という)のフリップ角φ1を90度＜φ1＜270度(好ましくは、φ1＝180度)として、前述の実施例1と同様に、このIR1によりラベリングする血流部を位相差に基づいてコントラスト強調する。他方のIRパルス(以下IR2という)のフリップ角φ2を0度＜φ2≦90度(好ましくは、φ2＝90度)として、このIR2によりラベリングする血流部を信号強度差によりコントラスト強調する。
そこで、待ち時間(TI)に関する撮像条件として、特に、IR1に対応する待ち時間TI1はフリップ角φ1に励起された縦磁化がnullにT1回復する時間より短くする必要がある。即ち、
TI1＜limit TI(φ1)＜limit TI(φ1＝180度)
とする。なお、IR2に対応する待ち時間TI2は、IR2によりラベリングする血流部の縦磁化がT1回復して所望の信号強度となる時間とする。TI1<TI2であれば、IR2-IR1の順序で、TI1>TI2であれば、IR1-IR2の順序で、各IRパルスを印加する。
また、流れ方向の撮像領域幅(FOVb)に関する撮像条件として、動脈と静脈の流速をそれぞれVba,Vbvとすると、流れ方向の撮像領域幅(FOVb)の上限値limit FOVbは、IR1で動脈を、IR2で静脈をそれぞれラベリングする場合は、
FOVb＜limit FOVb＝MIN((Vba×TI1),(Vbv×TI2)) (8)
となり、逆に、IR2で動脈を、IR1で静脈をそれぞれラベリングする場合は、
FOVb＜limit FOVb＝MIN((Vba×TI2),(Vbv×TI1)) (9)
となる。ここで、MIN(a,b)をa,bのいずれか小さい方をあらわす関数とする。一般的には、Vba＞Vbvであることから、静脈を最初にラベリングして、次に動脈をラベリングする順序が、流れ方向に広い撮像領域幅(FOVb)を確保する点で好ましく、最初に静脈を、IR2でラベリングする場合(TI1<TI2となる)には(8)式で、IR1でラベリングする場合(TI1>TI2となる)は(9)式で、それぞれlimit FOVbが決定される。
以下、最初にIR1で静脈をラベリングし、次にIR2で動脈をラベリングする場合を説明する。この場合、待ち時間の関係は、limit TI＞TI1＞TI2となる。なお、動脈を最初にIR1又はIR2でラベリングすることも可能である。
次に、本実施例のラベリング領域の設定について図10に基づいて説明する。図10は、下肢部における動脈(実線)801及び静脈(点線)802に対する、各IRパルスによるラベリング領域の設定例を示す。IR1(φ1＝180度)で静脈をラベリングするために、IR1によるラベリング領域1002を静脈に関して上流側(動脈に関して下流側)に設定し、IR2によるラベリング領域1001を動脈に関して上流側(静脈に関して下流側に設定し、ラベリング領域1001と1002の間に撮像領域(FOV)1003を設定する。操作者は、このように位置決め画像上で各ラベリング領域及び撮像領域を設定する。さらに、図10に示すように、操作者が、位置決め画像上でラベリング領域毎に、そのラベリングが対象とする流速を入力できるようにして、(8)(9)式に基づいて、流れ方向の撮像領域幅(FOVb)が最適化されるようにしても良い。最適な流れ方向の撮像領域幅(FOVb)が求められて設定された後に、各ラベリング領域の位置と幅の調整が操作者による手動又は装置による自動で行われる。なお、逆に、IR1でラベリング領域1001を励起して動脈をラベリングし、IR2でラベリング領域1002を励起して静脈をラベリングしてもよい。
次に、本実施例のパルスシーケンスを説明する。本実施例のパルスシーケンスは、前述の実施例1のパルスシーケンスと同様に、Main-ScanシーケンスとPre-Scanシーケンスとからなり、Pre-Scanシーケンスは前述の実施例1と同様なので詳細な説明は省略する。
本実施例のMain-Scanシーケンスの一例を図11に示すシーケンスチャートを用いて説明する。図11に示すMain-Scanシーケンスは、フリップ角及び待ち時間の異なる2つのIRパルス(1101,1102)を有してなるRFプリパルス部1110と、計測シーケンス部551とからなる。計測シーケンス部551は前述の実施例1と同様なので詳細な説明を省略し、ここでは、RFプリパルス部1110について詳細に説明する。
RFプリパルス部1110では、最初に、IR1(例えばフリップ角φ1＝180度)に該当するIRパルス1101とスライス傾斜磁場1103とが同時に印加されて図10のラベリング領域1002が励起され、静脈の血流がラベリングされる。この際のIRパルス1101は、スライス傾斜磁場1103で定まるラベリング領域1002の共鳴周波数と周波数帯域を有するSinc波とされる。そして、IRパルス1101の後にスポイル傾斜磁場パルス(1104-1〜1104-3)が3軸に印加されて、IRパルス1101により180度未満に励起されて生成された横磁化が消失される。
次に、IR2(例えばフリップ角φ2＝90度)に該当するIRパルス1102とスライス傾斜磁場1106とが同時に印加されてラベリング領域1001が励起され、動脈の血流がラベリングされる。この際のIRパルス1102は、スライス傾斜磁場1106で定まるラベリング領域1001の共鳴周波数と周波数帯域を有するSinc波とされる。そして、IRパルス1102の後にスポイル傾斜磁場パルス(1107-1〜1107-3)が3軸に印加されて、IRパルス1102により生成された横磁化が消失される。
なお、2つのスポイル傾斜磁場パルス(1104,1107)は、軸毎に纏めて一つにして、1回で印加されてもよい。
そして、IRパルス1101からの待ち時間がTI1であり、IRパルス1102からの待ち時間がTI2のタイミングで計測シーケンス部551の横磁化生成用RFパルス504が印加されて該計測シーケンス部551が開始する。
次に、本実施例の処理フローを説明する。本実施例の処理フローは、前述の実施例1で説明した図7に示したフローチャートに基づく処理フローと同じであるが、一部のステップの処理内容が一部異なる。そこで、処理内容が異なる処理ステップのみを説明する。また、図12に、図7に示す処理フローの各ステップを実施して得られる結果の一例を示す。以下、図7の対応ステップ番号に「-2」をつけて本実施例の処理ステップであることを明確にする。
ステップ701-2で、シーケンス実行部601は、図10に示すような位置決め画像を表示して、動脈の上流側のラベリング領域1001と、静脈の上流側のラベリング領域1002と、間の撮像領域1003の設定入力を受け付ける。そして、シーケンス実行部601は、操作者が設定入力した他の撮像条件と合わせて、撮像可能か否かをチェックし、撮像不可であればその旨を操作者に報知すると共に、ラベリング領域1001及び1002と撮像領域1003の位置及び幅の修正入力や、他の撮像条件についての撮像可能な条件を操作者に提示して該他の撮像条件の修正入力を受け付ける。そして最終的に撮像可能な撮像条件を確定し、この確定した撮像条件に基づいて、図11に示したMain-Scanシーケンスと図5(b)に示したPre-Scanシーケンスの実行に必要な各種制御データを具体的に計算する。特に、IRパルス1101及びスライス傾斜磁場1103、と、IRパルス1102及びスライス傾斜磁場1106と、をそれぞれ静脈用のラベリング領域1002と動脈用のラベリング領域1001を独立に励起するように、共鳴周波数及び周波数帯域を設定したSinc波形と傾斜磁場強度とする。
ステップ702-2で、前述の実施例1のステップ702と同様の処理が行われる。この結果、低空間分解能の位相画像(第1の位相画像)1201が得られる。この第1の位相画像1201は、IRパルス1101とIRパルス1102によってそれぞれラベリングされた静脈と、動脈の血流部及び他の静止部と、の間の位相差以外の色々な位相誤差を纏めて含む。
ステップ703-2で、前述の実施例1のステップ703と同様の処理が行われる。この結果、複素画像及びその絶対値画像1206と、該複素画像からその位相画像(第2の位相画像)1202とが得られる。
ステップ704-2で、前述の実施例1のステップ704と同様の処理が行われる。即ち、位相差画像演算部604は、ステップ702-2で得られた第1の位相画像1201を、ステップ703-2で得られた第2の位相画像1202と同じ空間分解能の位相画像に変換した後に、第2の位相画像1202との差分処理1221を行い、位相差画像1203を得る。この位相差画像1203は、共鳴周波数のずれに起因した位相誤差とハードウェアの不完全性に起因した位相誤差とが除去されて、IRパルス1101後の短い待ち時間(TI1)によって生じた位相差のみが反映された位相画像となる。
ステップ705-2で、前述の実施例1のステップ705と同様の処理が行われる。この結果、第1のマスク画像1208が得られる。
ステップ706-2で、前述の実施例1のステップ706と同様の処理が行われる。即ち、マスク処理部605はステップ704-2で得られた位相差画像1203に、ステップ705-2で作成した第1のマスク画像1208を施して(1222)、位相差画像1203から背景領域(ノイズ領域)を除外し被検体領域のみを抽出した位相差画像1204とする。
ステップ707-2で、前述の実施例1のステップ707と同様の処理が行われる。この結果、修正位相差画像においては、静止部及びIR2でラベリングされた血流部(動脈)の位相はゼロ、IR1でラベリングされた血流部(静脈)の位相はπとなる。
ステップ708-2で、前述の実施例1のステップ708と同様の処理が行われる。
ステップ709-2で、前述の実施例1のステップ709と同様の処理が行われる。この結果、第2のマスク画像1205が得られる。例えば、静脈、動脈、静止部の強度比をそれぞれ0:1:1と重み付けする第2のマスク画像1205とすることができる。つまり、第2のマスク画像1205は、静脈と動脈及び静止部との間の位相差に基づいて、静脈と動脈及び静止部との間のコントラストを強調するものとなる。
ステップ710-2で、前述の実施例1のステップ710と同様の処理が行われる。即ち、コントラスト強調処理部607は、ステップ709-2で得た第2のマスク画像1205をステップ703-2で得られた絶対値画像1206に施す(1223)。これによりコントラスト強調画像1210が得られる。
コントラスト強調画像1210においては、静止部に対して、ラベリングされた血流部(動脈及び静脈)が抑制された画像となる。具体的には、静脈802は位相差に基づく第2のマスク画像1205によって動脈及び静止部との間で0：1の強度比で重み付けられ、動脈801は絶対値画像1206における信号強度差に基づいて静止部との間で例えば0.5：1の強度比で重み付けられる。そして、これら両方の重み付けの結果として静脈、動脈、静止部の強度比がそれぞれ0:0.5:1で強調処理されたコントラスト強調画像1210となる。即ち、ラベリングされた血流部(動脈および静脈)と静止部との間のコントラストのみならず、動脈801と静脈802との間のコントラストも強調された画像となる。
なお、絶対値画像1206と位相差画像(1203,1204)がそれぞれに動脈を強調した画像と静脈を強調した画像に対応しており、領域抽出が容易なため、動脈801と静脈802に任意の色付けをすることも可能である。
以上までが、本実施例の処理フローの説明である。なお、上記説明では、IR1で静脈802を位相差強調し、IR2で動脈801を信号強度差強調する場合を説明したが、IR1で動脈801をラベリングして位相差強調し、IR2で静脈802をラベリングして信号強度差強調しても良く、その場合も同様に処理することができる。
以上説明したように、本実施例のMRI装置及び流体強調画像取得法は、流体部が流れ方向の異なる第1の流体部(静脈)と第2の流体部(動脈)とを有し、RFプリパルス部は、フリップ角の異なる第1のRFプリパルス(IR1)と第2のRFプリパルス(IR2)とを有し、第1のRFプリパルスを第1の流体部の上流側の第1の領域(ラベリング領域1002)に印加して該第1の流体部をラベリングし、第2のRFプリパルスを第2の流体部の上流側の第2の領域(ラベリング領域1001)に印加して該第2の流体部をラベリングし、第1の領域と第2の領域の間の領域を撮像領域とし、第1の流体部と第2の流体部とを異なるコントラストにして、流体強調画像を取得する。これにより、撮像時間を延長することなく、流れ方向の異なる血流部(動脈と静脈)を静止部に対して異なるコントラストで強調できると共に、動脈と静脈との間も異なるコントラストで強調できるようになり、短時間で動脈と静脈、及び静止部を明瞭に識別できる画像を取得することが可能になる。つまり、動静脈を分離して描出することが可能となる。
次に、本発明のMRI装置及び流体強調画像取得法の実施例3について説明する。本実施例は、2次元空間選択励起を用いてラベリングする。前述の実施例2のラベリングは、IRパルスとスライス傾斜磁場パルスとを組み合わせて、流れ方向に垂直に所定に幅の領域を励起することにより、その領域を通過する流体部をラベリングする方法であったが、ラベリング領域を通過する流体部の他に必要の無い静止部も励起されてしまう。そこで、実施例3では、流体部を選択的に励起する2次元空間選択励起を行う。さらに、位置の異なる流体部を、フリップ角の異なる2次元空間選択励起して、それぞれラベリングを行う。具体的には、一つに流体部には、2次元空間選択的なIR1を用いてラベリングし、他の一つの流体部には、2次元空間選択的なIR2を用いてラベリングを行う。なお、一つの血流部のみを2次元空間選択励起してラベリングしても良い。以下、流体部として動脈を例にして、図13〜16を用いて本実施例を詳細に説明する。
最初に、図13を用いて2次元空間選択励起の概要を説明する。図13は、xy平面上の形状のみが特定されたz方向に柱状領域を選択的に励起する例を示す。ここでは、xy平面上で特定される形状を円とする。RF,Gx，Gy,及びGzは、それぞれ、RFパルス、x軸方向の傾斜磁場、y軸方向の傾斜磁場、及びz方向の傾斜磁場の印加タイミングを示す。図13に示すように、RFパルス(2DRF)1301をx軸方向の振動傾斜磁場(Gx)1302およびy軸方向の振動傾斜磁場(Gy)1303とともに印加する。これにより、z軸に平行なシリンダ形状の領域が選択的に励起される。このようにRFパルスと傾斜磁場パルスとを組み合わせて2次元空間選択励起を行う際に、その形状と位置及びフリップ角の調整は、RFパルスの周波数及び波形と、傾斜磁場パルスの波形とを制御して行う。以下、2次元空間選択励起を用いたラベリングを空間選択ラベリングという。なお、2次元空間選択励起の詳細は、(非特許文献2)に説明されているので、ここでの詳細な説明は省略する。
次に、本実施例の空間選択ラベリングを行う領域の設定について図14に基づいて説明する。図14は、下肢部における左右の動脈801の各々に対する空間選択ラベリングを行う領域の設定例を示す。右側の動脈1401の上流領域の縦磁化をフリップ角φ1(例えば180度)で空間選択ラベリング(以下、2D-IR1という)を行い、左側の動脈の上流領域1402の縦磁化をフリップ角φ2(例えば90度)で空間選択ラベリング(以下、2D-IR2という)を行う。各空間選択ラベリングを行う領域は、xy平面内では略円形で、撮像平面(紙面)に垂直なz方向に柱状の領域である。したがって、実際にラベリングされる血流部は、血管と円柱状領域との交差部分となる。これにより、右側の動脈の血流部に対しては位相差強調を行い、左側の動脈の血流部に対しては信号強度差強調を行うことになる。なお、ラベリング順序はどちらが先でも良い。また、左側の動脈の血流部に対して2D-IR1を印加して位相差強調を行い、右側の動脈の血流部に対して2D-IR2を印加して信号強度差強調を行っても良い。
また、待ち時間(TI)の条件については、フリップ角φ1で空間選択ラベリングを行う2D-IR1の待ち時間(TI)をlimit TIより短くする必要があり、(2)式の制限に従う。他方、流れ方向の撮像領域幅(FOVb)1403の条件は、2番目に印加される空間選択ラベリングに対応する待ち時間によって、以下の(10)式により定まる。
FOVb＜limit FOVb＝MIN((2D-IR1のTI1),(2D-IR2のTI2))×Vba (10)
さらに、図14に示すように、操作者が、位置決め画像上でラベリング領域毎に、そのラベリングが対象とする流速を入力できるようにして、(8)(9)式に基づいて、流れ方向の撮像領域幅(FOVb)が最適化されるようにしても良い。
次に、本実施例のパルスシーケンスを説明する。本実施例のパルスシーケンスは、前述の実施例1のパルスシーケンスと同様に、Main-ScanシーケンスとPre-Scanシーケンスとからなり、Pre-Scanシーケンスは前述の実施例1と同様なので詳細な説明は省略する。本実施例のMain-Scanシーケンスは、前述の実施例2と同様に、フリップ角及び待ち時間の異なる2つの空間選択ラベリング部(1501,1502)を有してなるRFプリパルス部1510と、計測シーケンス部551とからなる。計測シーケンス部551は前述の実施例1、2と同様なので詳細な説明を省略し、ここでは、RFプリパルス部1510について図15示すシーケンスチャートを用いて詳細に説明する。
RFプリパルス部1510では、最初に右側の動脈の上流領域1401を空間選択ラベリングするための2D-IR1(1501)を実行し、次に、左側の動脈の上流領域1402を空間選択ラベリングするための2D-IR2(1502)を実行する。各空間選択ラベリング(1501,1502)における、RF、Gs,Gp,Grは、それぞれ図13で示したRF,Gz,Gx,Gyの波形を有して、それぞれの対応する励起領域及びフリップ角に調整された波形を有する。そして、2D-IR1(1501)からの待ち時間がTI1であり、2D-IR2(1502)からの待ち時間がTI2のタイミングで計測シーケンス部551の横磁化生成用RFパルス504が印加されて該計測シーケンス部551が開始する。
さらに、前述の実施例2で説明した図11に示すRFプリパルス部1110と同様に、2D-IR1(1501)の後にスポイル傾斜磁場パルス(1104-1〜1104-3)が3軸に印加されて、2D-IR2(1502)の後にスポイル傾斜磁場パルス(1107-1〜1107-3)が3軸に印加されて、それぞれ生成された横磁化が消失される。なお、本実施例においても、2つのスポイル傾斜磁場(1104,1107)は、軸毎に纏めて一つにして、1回で印加されてもよい。
次に、本実施例の処理フローを説明する。本実施例の処理フローは、前述の実施例1で説明した図7に示したフローチャートに基づく処理フローと同じであるが、一部のステップの処理内容が一部異なる。そこで、処理内容が異なる処理ステップのみを説明する。また、図16に、図7に示す処理フローの各ステップを実施して得られる結果の一例を示す。以下、図7の対応ステップ番号に「-3」をつけて説明する。
ステップ701-3で、シーケンス実行部601は、図14に示すような位置決め画像を表示して、動脈の上流側の右側ラベリング領域1401と左側ラベリング領域1402及び撮像領域1403の設定入力を受け付ける。そして、シーケンス実行部601は、操作者が設定入力した他の撮像条件と合わせて、撮像可能か否かをチェックし、撮像不可であればその旨を操作者に報知すると共に、右側ラベリング領域1401と左側ラベリング領域1402及び撮像領域1403の位置や幅(直径)の修正入力や、他の撮像条件についての撮像可能な条件を操作者に提示して該他の撮像条件の修正入力を受け付ける。そして最終的に撮像可能な撮像条件を確定し、この確定した撮像条件に基づいて、図15に示したMain-Scanシーケンスと図5(b)に示したPre-Scanシーケンスの実行に必要な各種制御データを具体的に計算する。特に、各空間選択ラベリング(1501,1502)がそれぞれ右側ラベリング領域1401と左側ラベリング領域1402を励起するよう、各RF,Gs,Gp,Grを設定する。
ステップ702-3で、前述の実施例1のステップ702と同様の処理が行われる。この結果、低空間分解能の位相画像(第1の位相画像)1601が得られる。この第1の位相画像1601は、空間選択ラベリング(1501,1502)によってそれぞれラベリングされた動脈の左側と右側の血流部と他の静止部との間の位相差以外の色々な位相誤差を纏めて含む。
ステップ703-3で、前述の実施例1のステップ703と同様の処理が行われる。この結果、複素画像及びその絶対値画像1606と、該複素画像からその位相画像(第2の位相画像)1602とが得られる。
ステップ704-3で、前述の実施例1のステップ704と同様の処理が行われる。即ち、位相差画像演算部604は、ステップ702-3で得られた第1の位相画像1601を、ステップ703-3で得られた第2の位相画像1602と同じ空間分解能の位相画像に変換した後に、第2の位相画像1602との差分処理1621を行い、位相差画像1603を得る。この位相差画像1603は、共鳴周波数のずれに起因した位相誤差とハードウェアの不完全性に起因した位相誤差とが除去されて、空間選択ラベリング1501後の短い待ち時間(TI1)によって生じた位相差のみが反映された位相画像となる。
ステップ705-3で、前述の実施例1のステップ705と同様の処理が行われる。この結果、第1のマスク画像1608が得られる。
ステップ706-3で、前述の実施例1のステップ706と同様の処理が行われる。即ち、マスク処理部605はステップ704-3で得られた位相差画像1603に、ステップ705-3で作成した第1のマスク画像1608を施して(1622)、位相差画像1603から背景領域(ノイズ領域)を除外し被検体領域のみを抽出した位相差画像1604とする。
ステップ707-3で、前述の実施例1のステップ707と同様の処理が行われる。この結果、修正位相差画像においては、静止部及び2D-IR2でラベリングされた左側動脈の血流部の位相はゼロ、2D-IR1でラベリングされた右側動脈の血流部の位相はπとなる。
ステップ708-3で、前述の実施例1のステップ708と同様の処理が行われる。
ステップ709-3で、前述の実施例1のステップ709と同様の処理が行われる。この結果、第2のマスク画像1605が得られる。例えば、右側動脈、左側動脈、静止部の強度比をそれぞれ0.5:1:1と重み付けする第2のマスク画像1605とすることができる。
ステップ710-3で、前述の実施例1のステップ710と同様の処理が行われる。即ち、コントラスト強調処理部607は、ステップ709-3で得た第2のマスク画像1605をステップ703-3で得られた絶対値画像1606に施す(1623)。この第2のマスク画像1605を用いた、即ち位相差画像1603に基づいた重み付け処理(1623)がコントラスト強調処理であり、このコントラスト強調処理によりコントラスト強調画像1610が得られる。
コントラスト強調画像1610においては、静止部に対して、ラベリングされた左右の動脈が抑制された画像となる。具体的には、右側の動脈部は位相差に基づく第2のマスク画像1605によって静止部との間で0.5:1の強度比で重み付けられ、左側の動脈部は絶対値画像1606における信号強度差に基づいて、静止部との間で0:1の強度比で重み付けられる。そして、これら両方の重み付けの結果として右側の動脈部、左側の動脈部、静止部の強度比がそれぞれ0.5:0:1でそれぞれ強調処理されたコントラスト画像1610となる。即ち、ラベリングされた両動脈と他の静止領域との間のコントラストのみならず、左右の動脈801間のコントラストも強調された画像となる。
なお、絶対値画像1606と位相差画像(1603,1604)がそれぞれに左側動脈を強調した画像と、右側動脈を強調した画像に対応しており、領域抽出が容易なため、左右の動脈を任意の色付けをすることも可能である。つまり、左右の動脈を分離した画像を得ることが可能である。
以上までが、本実施例の処理フローの説明である。なお、上記説明では、空間選択ラベリング1501で右側動脈を位相差強調し、空間選択ラベリング1502で左側動脈を信号強度差強調する場合を説明したが、空間選択ラベリング1501で左側動脈を位相差強調し、空間選択ラベリング1502で右側動脈を信号強度差強調する場合も同様に処理することができる。また、左右の動脈の内のいずれか一方のみを2次元空間選択励起してラベリングし、該動脈のみのコントラスト強調画像を得ても良い。
以上説明したように、本実施例のMRI装置及び流体強調画像取得法は、RFプリパルスを、2次元空間選択励起を行うRFパルスとし、撮像領域の上流側の流体部の少なくとも一部と交差する領域に2次元空間選択励起を行うRFパルスを印加して、該流体部の少なくとも一部をラベリングする。流体部は第1の流体部(右側の動脈)と第2の流体部(左側の動脈)とを有して、これらを異なるコントラストで描出する場合には、RFプリパルスは、2次元空間選択励起を行うフリップ角の異なる第1のRFプリパルス(2D-IR1)と第2のRFプリパルス(2D-IR2)とを有し、第1のRFプリパルスを、撮像領域の上流側の第1の流体部の少なくとも一部と交差する領域に印加し、第2のRFプリパルスを、撮像領域の上流側の第2の流体部の少なくとも一部と交差する領域に印加する。これにより、撮像時間を延長することなく、左右の動脈を静止部に対して異なるコントラストで強調できると共に、左右の動脈間も異なるコントラストで強調できるようになり、短時間で左右の動脈、及び静止部を明瞭に識別できる画像を取得することが可能になる。
次に、本発明のMRI装置及び流体強調画像取得法の実施例4を説明する。本実施例は、RFプリパルス部が2つのRFプリパルスを有してなり、最初のRFプリパルスを流体部の上流側の第1の領域に印加し、負方向にフリップされた縦磁化がnullにT1回復するタイミングで、第2のRFプリパルスを第1の領域に隣接する下流側の第2の領域に印加する。そして、撮像領域は、第2の領域内に設定される。前述の各実施例は、血流部の縦磁化を負方向の状態にしてエコー信号を計測するものであったが、本実施例は、逆に撮像領域における静止部の縦磁化を負方向の状態にしてエコー信号を計測する。本実施例は、流速が遅く、撮像領域内にラベリングした流体が十分に流れ込まないことが想定される場合に特に有効となる。
以下、RFプリパルスをIRパルスとして、流体部として血流を例にして、図17,18を用いて、本実施例を詳細に説明する。図17は、下肢部における動脈(実線)801及び静脈(点線)802に対する、ラベリング領域の設定例を示す。図18は、本実施例の、RFプリパルス部(RF)の各IRパルスの印加タイミングと、この各タイミングに合わせて第1の領域における流体部と第2の領域における静止部の縦磁化の挙動とをそれぞれ示す図である。
最初に、図17を用いて各ラベリング領域の設定について説明する。最初のフリップ角φ1(例えば180度)で血流をラベリングするIRパルス1801による第1のラベリング領域1701が上流側に設定され、次のフリップ角φ2(例えば180度)で血流をラベリングするIRパルス1802による第2のラベリング領域1702が、第1のラベリング領域1701に隣接する下流側に設定される。そして、撮像領域(FOV)1703が、第2のラベリング領域1702内に設定され、該第2のラベリング領域1702の少なくとも一部を撮像領域1703とする。図17では、撮像領域1703がラベリング領域1702中に設定された例を示す。操作者は、このように位置決め画像上で各ラベリング領域及び撮像領域を設定する。さらに、位置決め画像上で操作者が、ラベリング領域毎に、そのラベリングが対象とする流速を入力できるようにして、(8)(9)式に基づいて、流れ方向の撮像領域幅(FOVb)が最適化できるようにしても良い。最適な流れ方向の撮像領域幅(FOVb)が求められて設定された後に、各ラベリング領域の位置と幅の調整が操作者による手動又は装置による自動で行われる。
次に、本実施例のパルスシーケンスを説明する。本実施例のパルスシーケンスは、前述の実施例1のパルスシーケンスと同様に、Main-ScanシーケンスとPre-Scanシーケンスとからなり、Pre-Scanシーケンスは前述の実施例1と同様なので詳細な説明は省略する。また、Main-Scanシーケンスは、RFプリパルス部のみが異なり、計測シーケンス部は前述の実施例と同様なので、以下、図18を用いてRFプリパルス部のみを詳細に説明する。本実施例のRFプリパルス部1810は、2つのIRパルス(1801,1802)を有して成り、第1のIRパルス1801で第1のラベリング領域1701を励起し、次いで第2のIRパルス1802で第2のラベリング領域1702を励起する。
ラベリング領域1701に第1のIRパルス1801を印加後、血流の縦磁化がnullになるまで時間TI1だけ待ち、その後、第2のラベリング領域1702に対して第2のIRパルス1802を印加する。その結果、第2のIRパルス1802の印加直後は、第2のラベリング領域及び撮像領域1703に流れ込んだ血流の縦磁化はnull状態のままとなるが、第2のラベリング領域の他の静止部の縦磁化は180度反転して負方向の状態となる。その後さらに、TI2時間待つことで、第1のラベリング領域1701でラベリングされた血流がさらに撮像領域1703内に流れ込むと共に、縦磁化がT1回復して正方向の状態となる。他方、待ち時間(TI2)を、第2のラベリング領域における他の静止部の縦磁化がnullにT1回復する前までの時間であって、負方向を維持する時間とする。従って、第2のIRパルス1802から待ち時間(TI2)後の、計測シーケンス部551の横磁化生成RFパルス504の直前では、撮像領域1703における他の静止部の縦磁化は負方向に反転したままの状態になる。このような縦磁化状態で計測シーケンス部551により撮像領域1703からエコー信号を計測し、該エコー信号から再構成された画像においては、血流部と他の静止部の横磁化の位相をπ異ならせることができる。さらに、待ち時間(TI2)は、撮像領域1703の静止部の縦磁化が負方向を維持している時間であれば任意でよく、そのため、撮像パラメータの設定自由度を高くすることができる。
計測シーケンス部551によりエコー信号を計測した後の処理は、前述の実施例1と同様なので、詳細な説明は省略する。ただし、本実施例では、撮像領域における静止部の縦磁化を負方向の状態にしてエコー信号を計測するので、前述のステップ704で得られる位相差画像においては静止部の位相がπとなるが、前述のステップ707の修正位相差画像においては、静止部の位相が基準に修正されるので、最終的には、前述の各実施例と同様に、修正位相差画像においては、静止部の位相はゼロ、ラベリングされた血流部の位相はπとなる。従って、ステップ708以降の処理は同様になる。
以上説明したように、本実施例のMRI装置及び流体強調画像取得法は、RFプリパルス部は、第1のRFプリパルス(IRパルス1801)と第2のRFプリパルス(IRパルス1802)とを有し、第1のRFプリパルスを流体部の上流側の第1の領域(第1のラベリング領域1701)に印加して流体部をラベリングし、第2のRFプリパルスを第1の領域の下流側に隣接する第2の領域(第2のラベリング領域1702)に印加し、第2の領域の少なくとも一部を撮像領域とする。そして、第1のRFプリパルスにより負方向にフリップされた流体部の縦磁化がnullにT1回復するタイミングで第2のRFプリパルスを印加し、第2のRFプリパルスにより負方向にフリップされた第2の領域の静止部の縦磁化が負方向を維持する時間内に、計測シーケンス部により撮像領域からのエコー信号を計測する。これにより、遅い流速であっても、撮像時間を延長することなく、血流部を他の静止部に対してコントラストを強調した流体強調画像を得ることが可能になる。
以上、本発明の各実施例を説明したがが、本発明はこれらの実施例に限定されるものではない。
前述の各実施例の説明では、Pre-Scanにより位相誤差を抽出する例を説明したが、高度に調整されたMRI装置においては、位相誤差が少ないので、Pre-Scanの必要が無い場合もあり得るので、Pre-Scanを省略して、Main-Scanのみの実施でも本発明は成立する。即ち、高度に調整されたMRI装置においては、Main-Scanデータから得た位相画像に直接第1のマスク画像を施して得た位相画像に基づいて第2のマスク画像を得ても良い。
101 被検体、102 静磁場発生磁石、103 傾斜磁場コイル、104 RF送信コイル、105 RF受信コイル、106 信号検出部、107 信号処理部、108 全体制御部、109 傾斜磁場電源、110 RF送信部、111 計測制御部、112 ベッド、113 表示・操作部、114 演算処理部、115 記憶部

Claims (17)

  1. 所定のパルスシーケンスに基づいて、流体部を含む被検体の撮像領域からエコー信号の計測を制御する計測制御部と、
    前記エコー信号を用いて、前記流体部と静止部との間でコントラストを強調した画像を取得する演算処理部と、
    を有して成り、
    前記パルスシーケンスは、前記流体部の縦磁化を負方向に励起して該流体部をラベリングするRFプリパルスを備えるRFプリパルス部と、該ラベリングされた流体部が流れ込んだ撮像領域からエコー信号を計測する計測シーケンス部と、を有し、
    前記演算処理部は、前記画像の位相情報に基づいて、前記流体部と前記静止部との間のコントラストを強調した流体強調画像を取得する磁気共鳴イメージング装置であって、
    前記計測制御部は、前記撮像領域よりも上流側の領域に前記RFプリパルスを印加して前記流体部をラベリングし、該ラベリングされた流体部の縦磁化がnull以上に回復する前に、前記計測シーケンス部を実行する
    ことを特徴とする磁気共鳴イメージング装置。
  2. 請求項1記載の磁気共鳴イメージング装置において、
    前記RFプリパルスは、縦磁化をα(90＜α＜270)度に励起するIRパルスを含み、
    前記計測シーケンス部は、縦磁化をβ(0＜β≦90)度に励起して横磁化を生成する横磁化生成RFパルスを含む
    ことを特徴とする磁気共鳴イメージング装置。
  3. 請求項2記載の磁気共鳴イメージング装置において、
    前記計測制御部は、前記RFプリパルス部における前記IRパルスと、前記計測シーケンス部における前記横磁化生成RFパルスと、の間の待ち時間(TI)を、該IRパルスにより前記ラベリングされた流体部の縦磁化が負方向の状態を維持する時間とする
    ことを特徴とする磁気共鳴イメージング装置。
  4. 請求項3記載の磁気共鳴イメージング装置において、
    前記計測制御部は、前記ラベリングされた流体部と前記静止部とで、前記横磁化生成RFパルスの印加直後の横磁化の位相をπ異ならせて、前記エコー信号の計測を制御する
    ことを特徴とする磁気共鳴イメージング装置。
  5. 請求項1乃至4のいずれか一項に記載の磁気共鳴イメージング装置において、
    前記演算処理部は、前記計測シーケンス部で計測されたエコー信号を用いて再構成された画像の位相画像に基づいて、該画像の画素毎の重み係数を決定し、該重み係数を用いて前記流体強調画像を取得する
    ことを特徴とする磁気共鳴イメージング装置。
  6. 請求項5記載の磁気共鳴イメージング装置において、
    前記演算処理部は、
    前記計測シーケンス部のみで計測されたエコー信号を用いて再構成された画像の位相画像と、前記RFプリパルス部と前記計測シーケンス部とで計測されたエコー信号を用いて再構成された画像の位相画像と、の位相差画像に基づいて前記画素毎の重み係数を決定し、
    前記重み係数の分布を表すマスク画像を作成し、
    前記マスク画像を前記再構成された画像の絶対値画像に画素毎に掛け合わせて、前記流体強調画像を取得する
    ことを特徴とする磁気共鳴イメージング装置。
  7. 請求項6記載の磁気共鳴イメージング装置において、
    前記演算処理部は、前記位相差画像において、前記流体部の位相を[0〜1]の値に、前記静止部の位相を[1]に変換して前記重み係数とする
    ことを特徴とする磁気共鳴イメージング装置。
  8. 請求項1乃至7のいずれか一項に記載の磁気共鳴イメージング装置において、
    前記流体部は流れ方向の異なる第1の流体部と第2の流体部とを有し、
    前記RFプリパルス部は、フリップ角の異なる第1のRFプリパルスと第2のRFプリパルスとを有し、
    前記計測制御部は、前記第1のRFプリパルスを前記第1の流体部の上流側の第1の領域に印加して該第1の流体部をラベリングし、前記第2のRFプリパルスを前記第2の流体部の上流側の第2の領域に印加して該第2の流体部をラベリングし、前記第1の領域と前記第2の領域の間の領域を前記撮像領域とし、
    前記演算処理部は、前記第1の流体部と前記第2の流体部とを異なるコントラストにして、前記流体強調画像を取得する
    ことを特徴とする磁気共鳴イメージング装置。
  9. 請求項8記載の磁気共鳴イメージング装置において、
    前記第1のRFプリパルスは、フリップ角φ 1 (90＜φ 1 ＜270)度のIRパルスであり、
    前記第2のRFプリパルスは、フリップ角φ 2 (0＜φ 2 ≦90)度のIRパルスであり、
    前記計測制御部は、前記第1の流体部の縦磁化がnull以上に回復する前に、前記計測シーケンス部を実行し、
    前記演算処理部は、前記第1の流体部を前記位相情報に基づいてコントラスト強調を行い、前記第2の流体部を前記静止部との信号強度の差に基づいてコントラスト強調を行う
    ことを特徴とする磁気共鳴イメージング装置。
  10. 請求項1乃至7のいずれか一項に記載の磁気共鳴イメージング装置において、
    前記RFプリパルスは、2次元空間選択励起を行うRFパルスであり、
    前記計測制御部は、前記撮像領域の上流側の流体部の少なくとも一部と交差する領域に2次元空間選択励起を行うRFパルスを印加して、該流体部の少なくとも一部をラベリングする
    ことを特徴とする磁気共鳴イメージング装置。
  11. 請求項10記載の磁気共鳴イメージング装置において、
    前記流体部は第1の流体部と第2の流体部とを有し、
    前記RFプリパルスは、前記2次元空間選択励起を行うフリップ角の異なる第1のRFプリパルスと第2のRFプリパルスとを有し、
    前記計測制御部は、前記第1のRFプリパルスを、前記撮像領域の上流側の前記第1の流体部の少なくとも一部と交差する領域に印加し、前記第2のRFプリパルスを、前記撮像領域の上流側の前記第2の流体部の少なくとも一部と交差する領域に印加する
    ことを特徴とする磁気共鳴イメージング装置。
  12. 請求項1乃至7のいずれか一項に記載の磁気共鳴イメージング装置において、
    前記RFプリパルス部は、第1のRFプリパルスと第2のRFプリパルスとを有し、
    前記計測制御部は、前記第1のRFプリパルスを前記流体部の上流側の第1の領域に印加して前記流体部をラベリングし、前記第2のRFプリパルスを前記第1の領域の下流側に隣接する第2の領域に印加し、前記第2の領域の少なくとも一部を前記撮像領域とする
    ことを特徴とする磁気共鳴イメージング装置。
  13. 請求項12記載の磁気共鳴イメージング装置において、
    前記計測制御部は、前記第1のRFプリパルスにより負方向にフリップされた流体部の縦磁化がnullにT1回復するタイミングで前記第2のRFプリパルスを印加し、前記第2のRFプリパルスにより負方向にフリップされた第2の領域の静止部の縦磁化が負方向を維持する時間内に、前記計測シーケンス部により前記撮像領域からのエコー信号の計測を制御する
    ことを特徴とする磁気共鳴イメージング装置。
  14. 所定のパルスシーケンスに基づいて、流体部を含む被検体の撮像領域からエコー信号の計測を制御する計測ステップと、
    前記エコー信号を用いて、前記流体部と静止部との間でコントラストを強調した画像を取得する演算処理ステップと、
    を有し、
    前記パルスシーケンスは、前記流体部の縦磁化を負方向に励起して該流体部をラベリングするRFプリパルスを備えるRFプリパルス部と、該ラベリングされた流体部が流れ込んだ撮像領域からエコー信号を計測する計測シーケンス部と、を有し、
    前記演算処理ステップでは、前記画像の位相情報に基づいて、前記流体部と前記静止部との間のコントラストを強調した流体強調画像を取得する磁気共鳴イメージング装置における流体強調画像取得法であって、
    前記計測ステップでは、前記撮像領域よりも上流側の領域に前記RFプリパルスが印加されて、前記流体部の縦磁化がnullに回復する前に、前記計測シーケンス部が実行される
    ことを特徴とする流体強調画像取得法。
  15. 請求項14記載の流体強調画像取得法において、
    前記流体部は方向の異なる第1の流体部と第2の流体部とを有し、
    前記RFプリパルス部は、フリップ角の異なる第1のRFプリパルスと第2のRFプリパルスとを有し、
    前記計測ステップでは、前記第1のRFプリパルスが前記第1の流体部の上流側の第1の領域に印加されて該第1の流体部がラベリングされ、前記第2のRFプリパルスが前記第2の流体部の上流側の第2の領域に印加されて該第2の流体部がラベリングされ、前記第1の領域と前記第2の領域の間の領域が前記撮像領域とされ、
    前記演算処理ステップでは、前記第1の流体部と前記第2の流体部とを異なるコントラストにして、前記流体強調画像が取得される
    ことを特徴とする流体強調画像取得法。
  16. 請求項14又は15記載の流体強調画像取得法において、
    前記RFプリパルスは、2次元空間選択励起を行うRFパルスであり、
    前記計測ステップでは、前記撮像領域の上流側の流体部の少なくとも一部と交差する領域に2次元選択励起を行うRFパルスが印加されて、該流体部の少なくとも一部がラベリングされる
    ことを特徴とする流体強調画像取得法。
  17. 請求項14乃至16のいずれか一項に記載の流体強調画像取得法において、
    前記RFプリパルス部は、第1のRFプリパルスと第2のRFプリパルスとを有し、
    前記計測ステップは、前記第1のRFプリパルスが前記流体部の上流側の第1の領域に印加されて前記流体部がラベリングされ、前記第2のRFプリパルスが前記第1の領域の下流側に隣接する第2の領域に印加され、前記第2の領域の少なくとも一部が前記撮像領域とされる
    ことを特徴とする流体強調画像取得法。
JP2012553658A 2011-01-17 2012-01-11 磁気共鳴イメージング装置及び流体強調画像取得法 Expired - Fee Related JP5942271B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2011006551 2011-01-17
JP2011006551 2011-01-17
PCT/JP2012/050295 WO2012098955A1 (ja) 2011-01-17 2012-01-11 磁気共鳴イメージング装置及び流体強調画像取得法

Publications (2)

Publication Number Publication Date
JPWO2012098955A1 JPWO2012098955A1 (ja) 2014-06-09
JP5942271B2 true JP5942271B2 (ja) 2016-06-29

Family

ID=46515586

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012553658A Expired - Fee Related JP5942271B2 (ja) 2011-01-17 2012-01-11 磁気共鳴イメージング装置及び流体強調画像取得法

Country Status (4)

Country Link
US (1) US9435871B2 (ja)
JP (1) JP5942271B2 (ja)
CN (1) CN103327891B (ja)
WO (1) WO2012098955A1 (ja)

Families Citing this family (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5611882B2 (ja) * 2010-05-31 2014-10-22 株式会社東芝 磁気共鳴イメージング装置
CN103260510B (zh) * 2010-12-07 2016-01-20 株式会社日立医疗器械 磁共振成像装置以及对比度增强图像取得方法
JP5925529B2 (ja) * 2011-03-31 2016-05-25 株式会社東芝 磁気共鳴イメージング装置
JP5974391B2 (ja) * 2013-01-16 2016-08-23 株式会社日立製作所 磁気共鳴イメージング装置及びそのタイミングずれ検出方法
KR101458557B1 (ko) * 2013-02-20 2014-11-07 삼성전자주식회사 상이한 숙임각을 갖는 자기 공명 영상 시스템에서 주자장 정보 및 라디오 펄스 관련 정보를 획득하기 위한 방법 및 장치
JP6398106B2 (ja) 2013-04-08 2018-10-03 国立大学法人京都大学 磁気共鳴イメージング装置
JP6208464B2 (ja) * 2013-04-26 2017-10-04 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置
WO2015005297A1 (ja) * 2013-07-12 2015-01-15 株式会社 日立メディコ 磁気共鳴イメージング装置および磁気共鳴イメージング方法
DE102013215703B3 (de) * 2013-08-08 2015-02-05 Siemens Aktiengesellschaft Bestimmung einer T1-Zeit von Wasser und einer T1-Zeit von Fett
JP6548204B2 (ja) * 2014-04-23 2019-07-24 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置
JP6433679B2 (ja) * 2014-05-02 2018-12-05 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
CN106659420A (zh) * 2014-08-08 2017-05-10 株式会社日立制作所 磁共振成像装置
JP6452994B2 (ja) * 2014-08-26 2019-01-16 キヤノンメディカルシステムズ株式会社 画像処理装置及び磁気共鳴イメージング装置
CN104523274B (zh) * 2014-12-25 2016-11-09 中国科学院深圳先进技术研究院 一种利用稳态自由进动序列的磁共振成像方法
JP6363028B2 (ja) 2015-01-30 2018-07-25 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置およびプログラム
US10921410B2 (en) * 2015-02-03 2021-02-16 Koninklijke Philips N.V. Method and system for susceptibility weighted magnetic resonance imaging
WO2016149888A1 (zh) * 2015-03-20 2016-09-29 中国科学院深圳先进技术研究院 一种磁共振成像方法及设备
CN104688228B (zh) * 2015-03-20 2016-11-23 中国科学院深圳先进技术研究院 一种磁共振成像方法及设备
WO2017038345A1 (ja) * 2015-08-31 2017-03-09 株式会社日立製作所 磁気共鳴イメージング装置および撮像シーケンス生成方法
DE102015223658B4 (de) * 2015-11-30 2017-08-17 Siemens Healthcare Gmbh Verfahren zum Erfassen von Magnetresonanz-Signalen eines Untersuchungsobjekts
CN108369266B (zh) * 2015-12-23 2021-08-06 伯拉考成像股份公司 比率计量的脉冲cest成像
EP3446141A1 (en) 2016-04-21 2019-02-27 Koninklijke Philips N.V. Magnetic resonance imaging of arterial structures
EP3446142A1 (en) * 2016-04-21 2019-02-27 Koninklijke Philips N.V. Magnetic resonance imaging of arterial structures
KR101775028B1 (ko) * 2016-09-26 2017-09-05 삼성전자주식회사 자기 공명 영상 장치 및 자기 공명 영상 획득 방법
US11500052B2 (en) * 2017-01-31 2022-11-15 Regents Of The University Of Minnesota System and method for producing temporally resolved images depicting late-gadolinium enhancement with magnetic resonance imaging
US10684343B2 (en) * 2017-05-31 2020-06-16 Canon Medical Systems Corporation Magnetic resonance imaging apparatus and magnetic resonance imaging method
CN111225609B (zh) * 2017-10-19 2024-02-02 博放医疗有限公司 用于热疗法的处理系统和动态校正方法
JP7455508B2 (ja) 2018-12-26 2024-03-26 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置および医用複素数画像処理装置
CN114325523B (zh) * 2020-09-27 2023-10-03 上海联影医疗科技股份有限公司 T1值确定方法、装置、电子设备和存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07308302A (ja) * 1993-11-26 1995-11-28 Siemens Medical Syst Inc 磁気共鳴を用いての、組織または流体の選択された領域の像の形成およびバックグランドの抑圧方法
JP2001252263A (ja) * 2000-03-14 2001-09-18 Toshiba Corp 選択励起インバージョンパルスを用いた磁気共鳴イメージングの方法及び装置
JP2004261619A (ja) * 2004-06-25 2004-09-24 Toshiba Corp 磁気共鳴イメージング装置
JP2009160122A (ja) * 2007-12-28 2009-07-23 Ge Medical Systems Global Technology Co Llc Ｍｒｉ装置
JP2010046473A (ja) * 2008-07-24 2010-03-04 Toshiba Corp 磁気共鳴イメージング装置
JP2011254905A (ja) * 2010-06-07 2011-12-22 Toshiba Corp 磁気共鳴イメージング装置

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5270653A (en) 1992-04-10 1993-12-14 The Board Of Trustees Of The Leland Stanford Junior University Selective three dimensional excitation for magnetic resonance applications
JP3576069B2 (ja) * 2000-04-27 2004-10-13 株式会社東芝 Ｍｒｉ装置
JP4342809B2 (ja) * 2003-02-14 2009-10-14 株式会社東芝 Ｍｒｉ装置
EP2309285B1 (en) 2008-07-24 2019-08-21 Toshiba Medical Systems Corporation Magnetic resonance imaging apparatus for contrast enhancement of flow images
JP5461962B2 (ja) * 2009-02-05 2014-04-02 株式会社東芝 磁気共鳴イメージング装置
DE102009036237B4 (de) * 2009-08-05 2013-10-31 Siemens Aktiengesellschaft Kontrastmittelfreie MR-Angiographie mit SSFP-Sequenzen
CN102018511A (zh) * 2009-09-18 2011-04-20 株式会社东芝 磁共振成像装置以及磁共振成像方法
US9123121B2 (en) * 2010-09-27 2015-09-01 Hitachi Medical Corporation Magnetic resonance imaging apparatus and fluid imaging method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07308302A (ja) * 1993-11-26 1995-11-28 Siemens Medical Syst Inc 磁気共鳴を用いての、組織または流体の選択された領域の像の形成およびバックグランドの抑圧方法
JP2001252263A (ja) * 2000-03-14 2001-09-18 Toshiba Corp 選択励起インバージョンパルスを用いた磁気共鳴イメージングの方法及び装置
JP2004261619A (ja) * 2004-06-25 2004-09-24 Toshiba Corp 磁気共鳴イメージング装置
JP2009160122A (ja) * 2007-12-28 2009-07-23 Ge Medical Systems Global Technology Co Llc Ｍｒｉ装置
JP2010046473A (ja) * 2008-07-24 2010-03-04 Toshiba Corp 磁気共鳴イメージング装置
JP2011254905A (ja) * 2010-06-07 2011-12-22 Toshiba Corp 磁気共鳴イメージング装置

Also Published As

Publication number Publication date
CN103327891B (zh) 2016-08-17
JPWO2012098955A1 (ja) 2014-06-09
US20130293231A1 (en) 2013-11-07
WO2012098955A1 (ja) 2012-07-26
US9435871B2 (en) 2016-09-06
CN103327891A (zh) 2013-09-25

Similar Documents

Publication Publication Date Title
JP5942271B2 (ja) 磁気共鳴イメージング装置及び流体強調画像取得法
JP6998218B2 (ja) 動き検出を用いるｍｒ撮像
US9939505B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
US8588890B2 (en) Magnetic resonance imaging apparatus
JP5719968B2 (ja) Ｍｒデータを収集する方法及び装置
US5928146A (en) Inspection apparatus using nuclear magnetic resonance
US7372269B2 (en) Magnetic resonance imaging method and apparatus
US7965079B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP5848713B2 (ja) 磁気共鳴イメージング装置及びコントラスト強調画像取得方法
JP2004526491A (ja) 拡散強調された磁気共鳴画像化データの取得方法および装置
JP5740307B2 (ja) 磁気共鳴イメージング装置及び傾斜磁場印加方法
JP2007117765A (ja) Ｅｐｉシーケンスにおける傾斜誘起された交差項磁場の測定および補正
US20160169999A1 (en) Apparatus and method for magnetic resonance imaging with high spatial temporal resolutions
JP2020522356A (ja) デュアルエコーディクソン型水／脂肪分離ｍｒ撮像
US20160349346A1 (en) Intrinsic navigation from velocity-encoding gradients in phase-contrast mri
JP4040745B2 (ja) Ｍｒ装置
JP2005111059A (ja) 磁気共鳴イメージング装置
JP3688795B2 (ja) 磁気共鳴イメージング装置
Munoz et al. Innovations in cardiovascular MR and PET-MR imaging
JP2005288026A (ja) 磁気共鳴イメージング装置、渦磁場分布推定方法、及び静磁場補正方法
JP6449374B2 (ja) 磁気共鳴イメージング装置
JP2016040024A (ja) 磁気共鳴イメージング装置
JPH09187436A (ja) 磁気共鳴診断装置における３次元血流速度画像投影方法
WO2012005137A1 (ja) 磁気共鳴イメージング装置及びｒｆパルス制御方法
JP2021520951A (ja) 磁化の反転状態の評価を伴う動脈スピンラベリング法

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20141222

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20141222

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160105

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160126

TRDD Decision of grant or rejection written
RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20160330

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20160405

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20160427

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160428

R150 Certificate of patent or registration of utility model

Ref document number: 5942271

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees