JP6656826B2 - 磁気共鳴イメージング装置 - Google Patents

磁気共鳴イメージング装置 Download PDF

Info

Publication number
JP6656826B2
JP6656826B2 JP2015125009A JP2015125009A JP6656826B2 JP 6656826 B2 JP6656826 B2 JP 6656826B2 JP 2015125009 A JP2015125009 A JP 2015125009A JP 2015125009 A JP2015125009 A JP 2015125009A JP 6656826 B2 JP6656826 B2 JP 6656826B2
Authority
JP
Japan
Prior art keywords
pulse
flip angle
slice
refocusing
center
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.)
Active
Application number
JP2015125009A
Other languages
English (en)
Other versions
JP2016007539A (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.)
Canon Medical Systems Corp
Original Assignee
Canon Medical Systems Corp
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 Canon Medical Systems Corp filed Critical Canon Medical Systems Corp
Publication of JP2016007539A publication Critical patent/JP2016007539A/ja
Application granted granted Critical
Publication of JP6656826B2 publication Critical patent/JP6656826B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image 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/5615Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5617Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using RF refocusing, e.g. RARE
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/288Provisions within MR facilities for enhancing safety during MR, e.g. reduction of the specific absorption rate [SAR], detection of ferromagnetic objects in the scanner room
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/4818MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/482MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory
    • G01R33/4822MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory in three dimensions
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563Image 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/5635Angiography, e.g. contrast-enhanced angiography [CE-MRA] or time-of-flight angiography [TOF-MRA]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/567Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution gated by physiological signals, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating
    • G01R33/5673Gating or triggering based on a physiological signal other than an MR signal, e.g. ECG gating or motion monitoring using optical systems for monitoring the motion of a fiducial marker

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Health & Medical Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • General Health & Medical Sciences (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Optics & Photonics (AREA)
  • Vascular Medicine (AREA)

Description

本発明の実施形態は、磁気共鳴イメージング装置に関する。
磁気共鳴イメージングは、静磁場中に置かれた被検体の原子核スピンを、そのラーモア(Larmor)周波数のRF(Radio Frequency)パルスで磁気的に励起し、励起に伴い発生する磁気共鳴信号のデータから画像を生成する撮像法である。磁気共鳴イメージングでは、被検体の生体組織に吸収されるRFパルスのエネルギーの指標である比吸収率(Specific Absorption Ratio:SAR)が、所定以下に抑えられることが求められている。
本発明が解決しようとする課題は、比吸収率の上昇を抑えつつ信号雑音比を効果的に高めることができる磁気共鳴イメージング装置を提供することである。
実施形態の磁気共鳴イメージング装置は、シーケンス制御部を備える。シーケンス制御部は、スライスエンコード方向において、3次元k空間の中心を含むハーフフーリエスライスで印加されるリフォーカシングパルスのフリップ角が他のハーフフーリエスライスで印加されるリフォーカシングパルスのフリップ角よりも大きくなるようにFASE(Fast Advanced Spin Echo)シーケンスを実行する。シーケンス制御部は、第1のスライスに印加するリフォーカシングパルスを第1のフリップ角で、第1のスライスよりk空間の中心に近い第2のスライスに印加するリフォーカシングパルスであって、実効TEパルス以外のリフォーカシングパルスを、第1のフリップ角より大きな第2のフリップ角で、第2のスライスよりk空間の中心に近い第3のスライスに印加するリフォーカシングパルスであって、実効TEパルス及び前記実効TEパルスの両側のパルス以外のリフォーカシングパルスを、第2のフリップ角より大きな第3のフリップ角で、第2のスライスに印加するリフォーカシングパルスであって、実効TEパルスであるリフォーカシングパルスを、第2のフリップ角より大きな第4のフリップ角で、第3のスライスに印加するリフォーカシングパルスであって、実効TEパルスの両隣のパルスを、第3のフリップ角より大きな第5のフリップ角で、第3のスライスに印加するリフォーカシングパルスであって、実効TEパルスであるリフォーカシングパルスを、第4のフリップ角及び前記第5のフリップ角より大きな第6のフリップ角で印加する。
図1は、実施形態に係るMRIシステムの一例を示すブロック図である。 図2は、従来技術に係るFASEイメージングにおけるパルスシーケンスの一例を示す図である。 図3は、従来技術に係るFASEイメージングにおけるk空間及びパルスシーケンスを示す図である。 図4は、実施形態に係るMRIシステムにおける処理を説明するためのフローチャートである。 図5は、実施形態に係るk空間及びパルスシーケンスを示す図(1)である。 図6は、実施形態に係るk空間及びパルスシーケンスを示す図(2)である。 図7は、実施形態に係るk空間及びパルスシーケンスを示す図(3)である。 図8は、実施形態に係る可変TR間隔の一例を示す図(1)である。 図9は、実施形態に係る可変TR間隔の一例を示す図(2)である。 図10は、その他の実施形態に係るk空間及びパルスシーケンスを示す図(1)である。 図11は、その他の実施形態に係るk空間及びパルスシーケンスを示す図(2)である。 図12は、その他の実施形態に係るk空間及びパルスシーケンスを示す図(3)である。
以下、図面を参照して、実施形態に係る磁気共鳴イメージング(Magnetic Resonance Imaging:MRI)装置を説明する。以下では、MRI装置の一例として、MRIシステムについて説明する。
(実施形態)
図1は、実施形態に係るMRIシステムの一例を示すブロック図である。図1に示すMRIシステムは、ガントリ10(概略断面で示す)と、これに接続された各種の関連システム構成要素20とを有する。少なくともガントリ10は、通常、シールドルーム内に配置される。図1に示すMRIシステムの構造は、実質的に同軸の円筒形に配置された静磁場B磁石12と、Gx、Gy、及びGzの傾斜磁場コイル14と、大型の全身用RF(Radio Frequency)コイル(Whole Body RF Coil:WBC)16とを有する。この円筒形に配置される要素の横軸に沿って、被検体用寝台11によって支持された被検体9の頭部を実質的に取り囲むように、イメージングボリューム18が示される。1又は複数のより小型のアレイRFコイル19を、イメージングボリューム18内で被検体の頭部に、より近接して結合してもよい。当業者には明らかなように、表面コイル等のように、WBCと比較して小さいコイルやアレイコイルは、特定の身体部分(例えば、腕、肩、肘、手首、膝、脚、胸、背骨等)に合わせて設計されることが多い。本明細書では、そのような小型RFコイルを、アレイコイル(Array Coil:AC)又はフェーズドアレイコイル(Phased Array Coil:PAC)と呼ぶ。これらは、RF信号をイメージングボリューム18内に送信するよう構成された少なくとも1つのコイルと、イメージングボリューム18において、上記の例における被検体の頭部等の被検体9からのRF信号を受信するよう構成された複数の受信コイルとを含んでもよい。心臓信号収集装置8は、例えば、被検体9に適宜取り付けられ、MRIシーケンスコントローラ30のトリガとなる抹消脈波ゲーティング信号や、心臓ゲーティング信号13を供給するよう使用される。
MRIシステムコントローラ22は、ディスプレイ24、キーボード26、及びプリンタ28に接続された入出力ポートを有する。当然のことながら、ディスプレイ24は、制御入力もできるようにタッチスクリーンタイプのものであってもよく、マウス等の入出力装置を設けてもよい。
MRIシステムコントローラ22はMRIシーケンスコントローラ30に接続され、MRIシーケンスコントローラ30は、Gx、Gy、及びGzの傾斜磁場コイルドライバ32、ならびにRF送信機34及び送受信スイッチ36(同じRFコイルが送信と受信の両方に使用される場合)を制御する。MRIシーケンスコントローラ30は、MRIイメージング(核磁気共鳴(Nuclear Magnetic Resonance:NMR)イメージングとしても知られている)技術を実装するための適切なプログラムコード構造38を含み、そのような技術にはパラレルイメージングも含まれ得る。以下に説明するように、診断MRI画像の元であるNMRエコーデータ(「エコーデータ」)を得るために、所定のパルスシーケンス(以下、単に「シーケンス」とも表記する)や設定パラメータに応じて生成されたパルスシーケンスを実行するようにMRIシーケンスコントローラ30を構成してもよい。実施形態によっては、パルスシーケンスはシングルショット高速スピンエコー(Single-Shot Fast Spin Echo:SS−FSE)法とも呼ばれる超高速スピンエコー(Fast Advanced Spin Echo:FASE)法、又はその同様のものであってよい。MRIシーケンスコントローラ30は、心臓信号収集装置8から受信した心臓ゲーティング信号により、パルスシーケンスを被検体の部位に送信し、心臓の収縮、拡張にそれぞれ対応するエコーデータを収集することで差分画像を得るように構成されていてもよい。また、MRIシーケンスコントローラ30は、EPIイメージングやパラレルイメージング用に構成されてもよい。さらに、MRIシーケンスコントローラ30により、1回以上の準備スキャン(プリスキャン)シーケンス、及びメインスキャンMR画像(診断画像と呼ぶこともある)を取得するためのスキャンシーケンスを容易に実行できる。
MRIシステムは、ディスプレイ24に送られる処理画像データを作成するために、入力をデータプロセッサ42に送るRF受信器40を有する。また、MRIデータプロセッサ42は、前に生成されたMRデータ、画像やマップ、システム構成パラメータ等を記憶するMRIメモリ46、画像再構成/減算プログラムコード構造44、及びプログラム記憶装置50にアクセスするように構成される。
また、図1に、MRIシステムのプログラム記憶装置50の一般的な説明を示す。MRIシステムのプログラム記憶装置50では、(例えば、画像再構成のため、以下に説明するような差分画像等の生成のため、選択されたMRI画像特性のシミュレーションのため、MRIの後処理のための)格納されたプログラムコード構造が、MRIシステムの各種データ処理構成要素へアクセス可能な非一時的コンピュータ可読記憶媒体に格納される。更に、オペレータの入力を受け付けるためのプログラムコードが記憶されていてもよい。オペレータが入力するものとして、例えば、超高速スピンエコー(Fast Advanced Spin Echo:FASE)等のパルスシーケンスにおけるリフォーカシングパルスのリフォーカシングフリップ角等のパルスシーケンスの設定や、TR(Repetition Time)間隔の設定等が挙げられる。当業者には明らかなように、プログラム記憶装置50をセグメント化して、少なくとも一部分を、システム構成要素20の処理コンピュータのうち、通常操作においてそのような格納されたプログラムコード構造を最優先で必要とする別のコンピュータに直接接続してもよい(すなわち、MRIシステムコントローラ22に普通に格納したり直接接続したりするのではなく)。
実際に、当業者には明らかなように、図1は、後述する例示的な実施形態を実現するために変更された典型的なMRIシステムの、非常に大まかな概略図を示したものである。システム構成要素20は様々な論理集合の「ボックス」に分割することができ、通常、多数のデジタル信号プロセッサ(Digital Signal Processor:DSP)と、マイクロプロセッサと、専用処理回路(例えば、高速AD(Analog to Digital)変換用、高速フーリエ変換用、アレイ処理等用)とを含む。これらのプロセッサの各々は、通常、クロック制御された「状態マシン」であり、物理データ処理回路は、クロックサイクル(又は、所定数のクロックサイクル)毎に、ある物理状態から別の物理状態に移る。
処理回路(例えば、CPU(Central Processing Unit)、レジスタ、バッファ、演算装置)の物理状態が、操作過程であるクロックサイクルから別のクロックサイクルに徐々に変化するだけでなく、関連データ記憶媒体の物理状態(例えば、磁気記憶媒体内のビット記憶場所)も、このようなシステムの操作過程において、ある状態から別の状態に変換される。例えば、画像再構成処理や、時として以下に説明するようなコントロール画像及びタグ画像からの差分画像生成の終了時に、物理的記憶媒体内のコンピュータ可読でアクセス可能なデータ値の記憶場所の配列は、ある先行状態(例えば、全て一様に「0」値、又は全て「1」値)から新しい状態に変換され、そのような配列における物理的場所の物理状態は、最小値と最大値との間で変化して、実世界の物理的事象及び物理的条件(例えば、イメージングボリューム空間内の被検体の内部物理構造)を表す。当業者には明らかなように、命令レジスタに順次読み込まれMRIシステムの1つ以上のCPUによって実行されたときに、MRIシステム内で特定のシーケンスの動作状態を引き起こし遷移させる特定構造のコンピュータ制御プログラムコードと同様に、そのような格納データ値の配列は物理的構造を表し構成する。
造影MRI技術は、臨床の現場において心血管撮影に日常的に利用され、通常短いスキャン時間で強調画像を形成するが、被検体への造影剤の注入を伴う。また、動静脈を明確に分離するためには、造影剤の量、注入速度について更なる研究の余地があると言わざるを得ない。これは、動脈相、静脈相が重なる人体の下肢部において特に顕著である。したがって、磁気共鳴血管撮影法(Magnetic Resonance Angiography:MRA)をはじめとする数多くのスキャン法において、造影剤を使用しない非造影技術が有用である。
多くの場合、非造影磁気共鳴血管撮影法(NC(non-contrast)−MRA)のシーケンスとして、3D FASE(Fast Advanced Spin Echo:FASE、Single-Shot FSE:SS−FSEとも呼ばれる)が利用される。例えば、3D NC−MRAを鮮血イメージング(Fresh Blood Imaging:FBI)法における心臓ゲーティングと組み合わせることで、血管が高信号強度を示すT2強調画像のような画像が得られる。
非特許文献1(参照によりその全体が本明細書に組み込まれる)において、Miyazakiらは、スライスエンコーディング毎にin−plane3DハーフフーリエFSEを心電図ゲーティングと同期して行うNC−MRA FBI法について記載している。非特許文献1には、ハーフフーリエFSEの様々な要素がNC−MRA FBI法を効果的にしていると述べている。具体的には、多くのFSE系シーケンスでは「ブラックブラッド(黒く映った状態)」の画像になることに対して、当該方法では「ブライトブラッド(明るく映った状態)」として描出されているのである。ブライトブラッド画像の描出に寄与する要素として、エコートレインスペーシング(Echo Train Spacing:ETS)長の短縮、血管の向きに合わせた位相エンコード(Phase Encoding:PE)方向、ハーフフーリエFSEにおけるk空間の中心近辺又は低周波部分からの信号の収集等が挙げられる。
ハーフフーリエFSEにおけるETS長の短縮により、シングルショット収集時間が短縮され、効果的に運動関連アーチファクトを低減し、磁化率効果の最小化が図られる。ETSが短いほど、ハーフフーリエFSEのシングルショット収集ウィンドウが小さくなり、さらにハーフフーリエによる画像再構成を行うことで、256×256マトリクスに対するシングルショット収集時間を1秒未満まで短縮可能となる。
位相エンコード方向を血管の向きに合わせることにより、隣接画素間における重複したT2信号のボケが起きた状態から信号強調が得られ、「ブライトブラッド」信号取得に寄与する。
さらに、信号がハーフフーリエFSEにおけるk空間の中心近辺又は低周波部分から収集されることで、読出し方向と比較して、位相エンコード方向では傾斜振幅が低くなり、フローディフェーズが生じにくくなる。
しかしながら、血管特異性は、収縮(例えば静脈が明るく動脈が暗い)や拡張(例えば動静脈ともに明るい)等によるトリガ遅延時間に多かれ少なかれ依存する。上記のFBI法を流速の早い血管で実施すると、拡張期画像では動静脈共にブライトブラッドとして映り、収縮期画像では動脈がブラックブラッドとして映り、静脈はブライトブラッドとして映る。通常、上述のFBI法は、位相エンコード方向が血管の向きに平行な状態で用いられる。したがって、流速の早い血管では、拡張期画像から収縮期画像を減算することで容易に動脈のみの画像が得られる。
しかし、上記のFBI法を末梢部又は流速の遅い血管で用いると結果が異なる。流速が遅いため、収縮期画像においても動脈がブライトブラッドとして映るのである。このため、動静脈の分離が困難となる。非特許文献2(参照によりその全体が本明細書に組み込まれる)において、Miyazakiらは、流速の遅い血管でも良好な画像を得るために、上述のFBI法にフロースポイルド傾斜手法を使用することについて記載している。位相エンコード方向ではなく、読出し方向を血管の向きと平行とすることで、自然にディフェージング効果が得られる。さらにフロースポイルド傾斜パルスを読出し方向にかけることで、フローディフェージング効果が促進される。
即ち、フロースポイルド傾斜パルスは、拡張期に比較的静止した状態の血液や、静止状態の背景組織の信号強度に大きな影響を及ぼさないため、これら二組の拡張期、収縮期画像の差分をとることで動脈画像が得られるのである。
フロースポイルド傾斜パルスは、静止状態の背景組織の信号強度に影響を及ぼすことはない。同様に、静脈は心周期を通じて、比較的遅く一定の流速に保たれるため、拡張期、収縮期ともにフロースポイルド傾斜パルスにより大きく影響されることはない。したがって、フロースポイルド傾斜パルスをかけることで、拡張期と収縮期との間における動脈信号強度差を増幅できるため、拡張期と収縮期の差分から動脈がより良好に描出できる。血管の流速に応じてフロースポイルド傾斜パルスの強度を変化させる場合がある。即ち、流速の遅い動脈では、拡張期と収縮期との間で信号強度に差が出るよう、より高強度のフロースポイルド傾斜パルスが必要となる。
非特許文献3において、Nakamuraらは、上述のフロースポイルドFBI法により、拡張期、収縮期データの信号差を利用して末梢動脈樹全体が描出可能であると述べている。この方法では、トリガ遅延時間の精度や、流量依存性の読出しスポイラー傾斜パルスの選択に応じて画質が決まる。Nakamuraらによると、フロースポイルドFBI法は1.5Tで高い性能を示す。
1.5Tで評価されるFBIだが、3Tでは比吸収率(Specific Absorption Rate:SAR)の問題が生じ、組織や血液のT1時間が長くなってしまう。より具体的には、3T(1.5Tを超える場合)では組織及び血液のT1時間がどうしても長くなってしまうため、信号雑音比(Signal to Noise Ratio:SNR)を適正にするには、TR間隔を延ばす必要が生じる。しかし、TR間隔を延ばすことで、スキャン時間(撮像時間)が延びてしまう。
したがって、以下のいくつかの実施形態では、NC−MRA用のFBI法においてSARを低減し、SNRを適正にすることを対象とする。より具体的には、従来のNC−MRA法では、固定されたリフォーカシングフリップ角(「リフォーカシングパルスのフリップ角」、「フロップ角」とも称する)でリフォーカシングパルスを印加する。これに対し、実施形態では、SARが所望の領域内にとどまるよう、TR間隔内でリフォーカシングフリップ角を可変とする。この際、リフォーカシングフリップ角は、k空間において重要な中心領域用のエコーデータを得るのに十分な大きさとすることができる。様々な大きさのリフォーカシングフリップ角が、k空間の中心に対応するリフォーカシングパルス(例えば、k空間の中心にマッピングされるエコー信号が生成されるリフォーカシングパルス)を中心に、対称的に配置されてよい。実施形態によっては、例えば、リフォーカシングフリップ角は、対応するエコー信号がk空間中心から少なくとも一方向に離れてマッピングされる程度に応じて変化する。さらに、実施形態のいくつかでは、リフォーカシングフリップ角のみでなく、TR間隔も変化させる。これにより、総スキャン時間を許容範囲内にした上で、十分に長いTRにより収集されたエコー信号に基づき、k空間の中心領域を設定できる。
言い換えると、シーケンス制御部としてのMRIシーケンスコントローラ30は、位相エンコード方向及びスライスエンコード方向の少なくとも一方において、k空間の中心を含む領域で印加されるリフォーカシングパルスのフリップ角が他の領域で印加されるリフォーカシングパルスのフリップ角よりも大きいシーケンスを実行する。
従来技術において、リフォーカシングフリップ角を変化させることは、TEを延ばすこと、エコートレインを一定振幅とすること、エコー振幅を平滑に変化させることのいずれかを目的とする。例えば、非特許文献4(参照によりその全体が本明細書に組み込まれる)において、Scottらは、固定されたリフォーカシングフリップ角の代わりに、可変リフォーカシングフリップ角によりT2強調を制御し、高信号でT2時間を短く保つことについて述べている。同文献中図1では、使用されるパルストレインにおいて開始点で大きなパルス振幅が緩やかに縮小し、60°程度の所定の低い数値に近接し(漸近的)、対象組織からほぼ一定の信号を取得し、信号のフリップ角は所定の信号展開法により最適化される。
可変リフォーカシングフリップ角によるT2制御を対象とした従来技術に対し、実施形態では、k空間の中心では大きく保たれたリフォーカシングフリップ角が、その他箇所では低くなっている。実施形態によっては、対応するエコーがk空間中心にマッピングされるリフォーカシングパルスを中心に、対応する可変エコー振幅が対称的(又は、有効TEエコーがパルストレインの中心点となるとは限らないため実質対称的)に配置されるよう可変リフォーカシングフリップ角が設定される。さらに、実施形態においては、NC−MRA FBI等のNC−MRAのような、例えば収縮期、拡張期エコー信号を収集するため、上述のSARやスキャン時間の問題がより深刻となる環境において、可変リフォーカシングフリップ角が採用される。
図2は、従来技術に係るFASEイメージングにおけるパルスシーケンスの一例を示す図である。図2には、上述のフロースポイルドFBI法に利用されるパルスシーケンス等の、FASE(又はFSE)用の従来のパルスシーケンスを示す。心電図のR波202は、心臓ゲーティングに利用され、収縮期、拡張期のデータは、心周期を通じてそれぞれ一定のタイミングで収集される。画像データ収集が開始されるトリガ点も、R波202と関連して図示されている。R波202から所定のトリガ遅延時間が経過した後に、データが収集される。従来では、このトリガ遅延時間は、固定された値である(例えば、2RR間隔毎に1回)。3DハーフフーリエFSEシーケンスは、スライスエンコーディング毎に心電同期されるため、全スライスパーティションで、心拍位相が同一となる。図2には、R波202に基づくスライスエンコーディング1とスライスエンコーディング2が示されている。RFパルスシーケンス及び対応するエコートレイン204、さらに対応する読出し方向の傾斜磁場パルス206も図示されている。
トリガリング(つまり、R波202からトリガ遅延時間が経過したタイミング)毎に、スライスエンコーディングが開始され、FASEパルスシーケンス等のパルスシーケンスがイメージングの対象に送信される。図示のように、スライスエンコーディング用のRFパルスシーケンス(例えば、スライスエンコーディング1)は、少なくとも1つのパルストレインを含む。各パルストレインは、励起パルスと、時間的に励起パルスの後に続く複数の連続したリフォーカシングパルスとを有する。励起パルスは、励起フリップ角90°であってよく、各リフォーカシングパルスはリフォーカシングフリップ角180°であってよい。リフォーカシングパルス毎に1つのエコー信号が得られる。例えば、リフォーカシングフリップ角を変更することでエコー信号の強度(例えば振幅)を変更することが可能であること、さらにRF励起パルス、リフォーカシングパルスの振幅、パルス幅、周波数がSARに寄与することが、当業者には理解されよう。
図2はさらに、エコートレイン204で示されたエコー信号を得るために、読出し方向に印加される傾斜磁場パルス206を示す。図示された傾斜磁場パルス206は、フロースポイルド傾斜を含み、フロースポイル領域を実際の読出し傾斜磁場パルスの両側に有する。非特許文献2はフロースポイルド傾斜磁場パルス、さらにその構成と効果についても言及している。
図3は、従来技術に係るFASEイメージングにおけるk空間及びパルスシーケンスを示す図である。図3には、イメージングスライス1〜5それぞれに対するFASE RFパルスシーケンスの、従来の適用例を概略的に示す。図3は、さらに、従来のハーフフーリエ法における、k空間内でのイメージングスライス1〜5に対応するNMRのエコー信号のマッピングを概略的に示す。例示的に図示されたk空間301は、スライス302等の、1〜5の複数のスライスを有する。読出し方向、位相エンコード方向、スライスエンコード(Slice Encoding:SE)方向(それぞれkx、ky、kzとする)は、k空間に対して定義される。3D FASEシーケンスは複数層のスライスに対して実行される。各スライス(例えばスライス1〜5)にFASEパルスシーケンス304を実行することで、エコー信号トレイン(図示省略)が得られる。各FASEパルスシーケンスは、励起パルス308及び後続の複数のリフォーカシングパルス310を有する。図3に概略的に示したPE、SE軸は時間軸座標も有することが当業者には理解されよう。各スライス1〜5内に、対応するRFパルスシーケンスにより、エコートレイン(図示省略)が生成される。あるRFパルスシーケンスにより生成されたエコートレインは、各リフォーカシングパルスにより生成された各エコー信号を有する。3D FASEシーケンスにおいて、シングルショットシーケンスにおける実効TE(Echo Time)がk空間の中心又は中心近傍(エコートレインのエコー内でk空間中心に最も近く)にあることが好ましい。なお、図3は、RFリフォーカシングパルスが全て同一のリフォーカシングフリップ角を有する従来技術を示す。SARに関連した問題として、リフォーカシングフリップ角が大きすぎる場合、SARが過剰にならないよう、TRの延長が必要となり得る。一方、リフォーカシングフリップ角が小さすぎる場合、SNRが十分に観測されない場合がある。
従来技術の3D FBIイメージングでは、特に3Tの場合、SARが高くなりすぎる(推奨/所望の安全性限界を超える)という問題がある。例えば「ブライトブラッド」を十分明確に表示するために、血管のSNRを高めるためには、リフォーカシングフリップ角を大きくする必要があるが、SAR限界により、リフォーカシングフリップ角の拡大及びTRの延長の少なくとも一方が規制される。したがって、スキャン時間が長くなってしまうのである。
3D FBIスキャンにおいて、血液信号を十分に高く保ちつつ、SARを所望範囲内にとどめるため、実施形態は、シングルショット(例えば、各々の又は多数のスライスエンコーディングショット)においてリフォーカシングフリップ角を可変とする。可変リフォーカシングフリップ角は、(例えば、各スライスに対応して)k空間の中心又は中心近傍に一又は複数の大きなリフォーカシングフリップ角が来るよう設定されてよい。これにより血液信号が十分に高く保たれる。k空間の中心又は中心近傍のスライスエンコーディングでは大きなリフォーカシングフリップ角が設定され、その他周辺でのスライスエンコーディングではより小さなリフォーカシングフリップ角が設定される。好適な実施形態では、3D FASEシーケンスの可変リフォーカシングフリップ角を、位相エンコード方向、スライスエンコード方向の両方に用いる(リフォーカシングフリップ角の二次元的変化)。
リフォーカシングパルスの可変リフォーカシングフリップ角の配置に関する、様々な候補のパターンを以下に例示する。
ショット毎にリフォーカシングフリップ角が可変のパターン。
スライスエンコーディングの中心では大きな可変リフォーカシングフリップ角、その他スライスエンコーディングではより小さなリフォーカシングフリップ角。
スライスエンコーディングの中心では大きな可変リフォーカシングフリップ角、その他スライスエンコーディングではやや大きい可変リフォーカシングフリップ角。
スライスエンコーディング範囲の中心では大きな可変リフォーカシングフリップ角、その他スライスエンコーディングではより小さなリフォーカシングフリップ角。
スライスエンコーディング範囲の中心では大きな可変リフォーカシングフリップ角、その他スライスエンコーディングではやや大きい可変リフォーカシングフリップ角。
いくつかの実施形態では、対応するエコー信号がk空間の中心にマッピングされるリフォーカシングパルスを中心に、異なる複数の大きいリフォーカシングパルスフリップ角が略対称に配置されるようなリフォーカシングパルスの構成を採る。言い換えると、シーケンスは、中心に対してフリップ角の大きさが略対称のリフォーカシングパルスを含む。
従来技術に対して、例示的な実施形態は、例えば、以下の効果を奏する。k空間ラインのうち、いくつか(例えば外側のもの)でリフォーカシングフリップ角を小さくすることでSARを低減し、中心側のk空間ラインではより大きなリフォーカシングフリップ角を設定することで血液信号のSNRを十分高く保つ。したがって、例えば外側のTRを短くし、スキャン時間を短縮する。
図4は、実施形態に係るMRIシステムにおける処理を説明するためのフローチャートである。図4には、1又は複数の実施形態における、FASEイメージング、例えばFBIイメージングにおけるSARを低減する処理400のフローチャートを示す。処理400は、例えば、図1に示すMRIシステム等のMRIシステムで実行されてもよい。処理400を実行する際、ステップS402〜ステップS414のうち1以上の操作を、図示した以外の順番で実行してもよく、又は実行しなくてもよく、又は1以上の他の操作と組み合わせてもよいことは当業者には理解されよう。
処理400が開始されると、ステップS402にて、被検体及びMRIシステムをスキャンのために準備する。ステップS402には、MRIシステムの送信RFコイルや受信RFコイルに対して撮像対象の被検体又は被検体の部位を位置決めする工程や、撮像を行うための一般的なパラメータや構成オプションを設定する工程が含まれてもよい。ここで説明される例示的な実施形態の技術は、システム構成及び被検体の位置を適宜設定した上で血管が検査される、頭部、頸部、心臓、肺、肝臓、脚、腸骨、大腿骨、腓腹、足、手等、被検体の多くの部位をイメージングするために利用される。イメージング対象は上記に限定されない。以下に説明するように、例えば、パルスシーケンス、スライスの箇所、大きさ等の要素は、対象画像の任意の特徴に合わせて調整可能である。例えば、撮像されている血管、又は身体もしくは臓器の特定部分の血流速度に応じて、構成を設定したり調整したりしてもよい。
実施形態によっては、準備段階として、例えば、被検体の位置合わせ及びコイル較正の少なくとも一方に供される1以上の低解像度MRI画像を得るよう、1以上のプリスキャンを収集する工程を含んでもよい。実施形態によっては、心電図の準備スキャンが1回以上行われても良い。これにより収集された画像から、対象となる少なくとも1つの血管(例えばイメージングされる少なくとも1つの血管/領域)に対して適切な心電図の遅延時間(トリガ遅延時間)を収縮期、拡張期に合わせて判定できる。このようにして選択された遅延時間が後に、スライスエンコーディング毎に心電図ゲーティングと同期された診断用3DハーフフーリエFSE収集に用いられる。実施形態によっては、さらにプリスキャンを実施し、例えば、取得されたプリスキャン画像に基づいて、フロースポイルド傾斜磁場パルスの強度を判断してもよい。上述のようなプリスキャンは、非特許文献2でも言及されており、それを例示的な実施形態で用いてもよい。
ステップS404において、可変リフォーカシングフリップ角及び関連する要素(例えば、TR間隔など)の設定パラメータが取得される。オペレータは、可変リフォーカシングフリップ角に関連する設定を、以下の手法の少なくとも1を通じて行うことができる。
・予め記憶されたパラメータ表から選択。
・各スライスや、単一スライスの各k空間のラインごとに適用される複数のパルスシーケンスの少なくとも1に対して、リフォーカシングフリップ角のパターンを指定。
実施形態によってはユーザインターフェース(UI)が設けられてよい。オペレータはUIを使用して、適用されるFASEパルスシーケンスの少なくとも1における少なくとも1のリフォーカシングフリップ角を個別に選択可能である。さらに、選択されたパルスに対して個別に、又はまとめてリフォーカシングフリップ角を設定可能である。実施形態によっては、ユーザは、実効TE(TE effective)に対応したリフォーカシングパルスのみに対してリフォーカシングフリップ角を設定してもよい。この実効TEに対応するエコー信号(以下、実効TEエコーとも表記する)が、k空間の中心に記録される。上述の構成はあくまで一例であり、その他様々なリフォーカシングフリップ角に関する設定が可能であることが当業者により理解されよう。
設定パラメータは、可変TR間隔に対しても任意に取得されてもよい。上述のように、TR間隔の長さもFBIスキャンに関連した総SARを増減させる重要な要素である。TR間隔は、当該分野で一般的に使用される用語であり、2つの励起パルス間の時間間隔である。実施形態においては、ある励起パルスに対する全てのリフォーカシングパルスはTR間隔内に収まる。
k空間中心スライス(例えばk空間中心に重なる、又は近傍のスライス)においてリフォーカシングフリップ角は大きいため、SARを所望の限界値以下に保つために長いTRが必要な場合等に、可変TR間隔は特に有効である。TR間隔の短縮はスキャン時間の短縮につながるが、スキャンに関連したSARが増加し、SNRが低下してしまう。一方、TR間隔が長いほど、スキャンタイムが長くなるが、SARが低減され、収集される信号の強度が増す。
オペレータは、異なる2つ以上のTR間隔長がスライスのエンコードのために設定されるように設定パラメータを定めてもよい。各TR間隔長においてエンコードされるスライスの数や割合も設定可能である。実施形態によっては、ユーザは、各スライスに対してTR間隔長を指定してもよい。
ステップS406において、取得した設定パラメータに応じて、収集パルスシーケンスが設定される。パルスシーケンスは、例えば、k空間の中心領域に寄与するリフォーカシングパルスが最大のリフォーカシングフリップ角を有するよう構成される。k空間の中心は、スライス(例えば2D平面)に関しても中心領域となり、複数スライス(例えば、スライス束又は3Dボリューム)に関しては3Dボリュームの中心となる。各パルスシーケンスにおけるリフォーカシングパルスは以下の条件を満たすよう配置される。
(1)各スライスにつき、k空間の中心に最も近いパルスに最大のリフォーカシングフリップ角が割り当てられ、k空間の中心から最も遠いパルスに最小のリフォーカシングフリップ角が割り当てられ、さらに、又は、それに代えて、リフォーカシングフリップ角がk空間中心から外側に向かうにつれて大きくなることがないよう、中間のリフォーカシングフリップ角が割り当てられる。
(2)スライス間では、3Dk空間の中心のスライスが最大のリフォーカシングフリップ角を有し、外側のスライスが最小のリフォーカシングフリップ角を有し、さらに、又は、それに代えて、スライスの2D中心におけるリフォーカシングフリップ角がk空間の中心から外側に向かうにつれて増大しない。
言い換えると、シーケンスは、3次元k空間の位相エンコード方向及びスライスエンコード方向において、中心を含む領域で印加されるリフォーカシングパルスのフリップ角が他の領域で印加されるリフォーカシングパルスのフリップ角よりも大きい。
これ以外の可変リフォーカシングフリップ角や可変TR間隔の配置も可能であり、それら設定も実施形態として挙げられる。例えば、(SARやスキャン時間が許容範囲を超えないという条件で)k空間の中心領域のみでなく、中心領域外の1以上の領域も大きいリフォーカシングフリップ角によるエコーデータ(エコー信号)を有する構成でもよい。
言い換えると、シーケンスは、位相エンコード方向において、パルストレインの開始点及び中心を含む領域で印加されるリフォーカシングパルスのフリップ角が他の領域で印加されるリフォーカシングパルスのフリップ角よりも大きい。
上記のように、設定パラメータに応じてパルスシーケンスを生成する際に、MRIシステムに連携したコンピュータのメモリに記憶された少なくとも1つの所定のルックアップテーブルを参照してもよい(例えば、パルスシーケンスにおける任意又は全てのリフォーカシングパルスに対する初期設定リフォーカシングフリップ角が記録された所定のルックアップテーブル)。ここで、最大リフォーカシングフリップ角、最小リフォーカシングフリップ角、及びリフォーカシングパルス数のうち少なくとも1が設定された場合、所定の演算機能により、各リフォーカシングパルスに対してリオフォーカシングフリップ角を決定する。又は、ユーザの指定に応じたパルスシーケンスが使用される(例えば、少なくとも1つの選択されたリフォーカシングパルスのリフォーカシングフリップ角がオペレータにより指定され、他は所定の初期設定値に設定されている場合)。
オペレータが可変TR間隔の設定も指示した場合、パルスシーケンスはさらに、k空間の3D中心近傍のスライスが最大のTR間隔を有するよう設定される。さらに、外縁のスライスが最短のTR間隔を有する設定や、中間のスライスはスライスTR間隔が中心から外縁に向かうにつれ増大しないような設定も可能である。
ステップS408では、上記の通り設定された設定済みのパルスシーケンスが実行される。この際、各パルスシーケンスをスライスに対し2度実行してもよい。この場合、1度目で収縮期に対応するエコー信号が取得され、2度目で拡張期に対応するエコー信号が取得される。さらに、ステップS406で設定された1つ又は複数の3D FASEパルスシーケンスに応じたRFパルスを、被検体の選択された部位に送信してもよい。
ステップS410において、実行された1つ又は複数のパルスシーケンスに対応するエコー信号が収集される。いかなる数のパルスシーケンスが実行されても、対応する応答エコー信号が得られるようにステップS408、ステップS410が実行されるということが当業者には理解されよう。実施形態によっては、ハーフフーリエ再構成に必要な全てのエコーデータは、収縮期、拡張期に対してそれぞれ1パルスシーケンスを実行することで収集される。その他の実施形態によっては、k空間のハーフフーリエ収集に必要な全エコーデータを取得するために2以上のパルスシーケンスを実行してもよい。
実施形態によっては、収集工程の結果、収縮期において取得された画像のデータセットと、同一の心周期における拡張期において取得された画像のデータセットが得られる。実施形態によっては、収縮期と拡張期のエコー信号は、例えば、未完成ではあるが完全な状態のMR画像のデータを提供可能な、k空間にそれぞれ対応する異なる記憶領域に収集される。
ステップS412において、画像(診断画像)が生成される。診断画像は、拡張期に収集されたエコー信号から得られた拡張期画像及び収縮期に収集されたエコー信号から得られた収縮期画像のうち一方から他方を減算することで生成される。
言い換えると、生成部としてのMRIシステムコントローラ22は、シーケンスによって収集されたエコー信号に基づいて、出力画像を生成する。例えば、シーケンス制御部としてのMRIシーケンスコントローラ30は、第1の心位相及び第2の心位相において、シーケンスをそれぞれ実行する。そして、生成部としてのMRIシステムコントローラ22は、第1の心位相において収集されたエコー信号から得られる第1の画像と、第2の心位相において収集されたエコー信号から得られる第2の画像との差分に基づいて、出力画像を生成する。
ステップS414において、生成された診断画像が出力される。生成された画像は、MRIシステムの表示部や、遠隔表示部を通じてオペレータに表示可能である。あるいは、生成された画像は、印刷されたり、ローカル、又は遠隔記憶装置に保存されたり、さらにMRI処理が施されたりしてもよい。
図5は、実施形態に係るk空間及びパルスシーケンスを示す図(1)である。図5には、ハーフフーリエ法により、NMRのエコーデータにより半分強が埋まった状態のk空間502を概略的に示す。図5に示すリフォーカシングパルス510の縦方向の大きさは、リフォーカシングフリップ角の大きさに対応する。k空間502の塗りつぶし領域は、単一のスライスに対して収集されたエコーデータを含む。実施形態において、イメージング対象物に対してFASEシーケンス504が実行される。FASEシーケンス504内の各リフォーカシングパルスに対応して生成されたエコーにより、エコートレイン(図示省略)が生成される。FASEシーケンス504は、励起パルス508と、後続のリフォーカシングパルス510を含む、複数のリフォーカシングパルス510を有する。生成されるエコートレインは、各リフォーカシングパルス510に対して生成されるエコーを有する。リフォーカシングパルス510の一つであるリフォーカシングパルス506は、k空間の中心にマッピングされるエコー信号(例えば、実効TEエコー)を生成する。リフォーカシングパルス506は、重要となるk空間の中心のエンコーディングを担うため、その他のリフォーカシングパルス510よりも大きなリフォーカシングフリップ角を有する。上述のように、リフォーカシングフリップ角が増大するにつれ、対応するエコー信号の強度(即ち振幅)が増し、SNRが上昇する。
図5に示すように、各2D FASEシングルショットは、実効TEエコーのリフォーカシングフリップ角がk空間の中心に来るようなk空間の可変リフォーカシングフリップ角のパターンを有してもよい。図示の例において、k空間の中心のリフォーカシングフリップ角は大きく設定されており、その他の部分のリフォーカシングフリップ角は、比較的小さく設定されている。k空間の中心の、リフォーカシングフリップ角が大きいリフォーカシングパルスの総数は、適宜設定可能である。
図6は、実施形態に係るk空間及びパルスシーケンスを示す図(2)である。図6には、位相エンコード(PE)方向、スライスエンコード(SE)方向の両方に適用された可変リフォーカシングフリップ角を概略的に示す。図6に示すリフォーカシングパルスの縦方向の大きさは、リフォーカシングフリップ角の大きさに対応する。血液信号のSNRを十分上げるため、PE、SEの両方向において、対応するエコーがk空間の中心にマッピングされるリフォーカシングパルスに大きなリフォーカシングフリップ角が割り当てられる。k空間602は5つのハーフフーリエスライスである、スライス1〜5を含む。FASEパルスシーケンス604により、スライス1〜5に対応するNMRのエコー信号が生成される(即ち、ハーフフーリエ法により、スライス1〜5は、FASEパルストレイン1〜5に対応するエコートレインで埋まる)。スライス3は、スライスエンコード方向においてk空間の中心に配置される。したがって、スライス3の中心点(即ち、スライス3平面において、位相エンコード方向、読出し(RO)方向それぞれの中心点から延びる直交軸の交点)は、3Dk空間の中心となる。したがって、スライス3の実効TEエコーは、3Dk空間の中心でエンコードされ、対応するリフォーカシングパルス606のリフォーカシングフリップ角は、スライス1〜5に割り当てられたリフォーカシングフリップ角のうち最大のものとなる。スライス2及び4はスライス3の両側に位置するため、SE方向においてk空間の中心に対し、スライス3の次に近いスライスとなる。したがって、図7に示す実施形態では、2番目に大きい(即ち、スライス3の実効TEエコーに対応するリフォーカシングパルス606に割り当てられたものの次に大きい)リフォーカシングフリップ角が、スライス3の実効TEエコーに隣接した、即ちスライス2及び4の実効TEエコーのリフォーカシングパルス606に割り当てられる。さらに、3Dk空間の中心近傍のスライス2〜4のリフォーカシングパルスのリフォーカシングフリップ角は、スライス1及び5のリフォーカシングパルスのリフォーカシングフリップ角よりも大きい。上述のように、リフォーカシングフリップ角の相対的な大きさは、対応するエコー信号の振幅により定まる。図7に示すリフォーカシングフリップ角の3Dパターンは、PE、SE「平面」における段階的な山と捉えることができる。多段階であるほど、「山」をなす要素がよりなだらかに変化することが理解されよう。
スライス1〜5のパルスシーケンス604は、例えば、以下の通りに設定される。スライス1及び5のリフォーカシングフリップ角はいずれも90°となる。スライス2及び4のリフォーカシングフリップ角(実効TEパルスのものは除く)は120°となる。スライス3のリフォーカシングフリップ角(実効TEパルス、及びその両側のパルスのものは除く)は120°となる。スライス2及び4の実効TEパルス、及びスライス3の実効TEパルスの両側のパルスのリフォーカシングフリップ角は150°となる。スライス3の実効TEパルスのリフォーカシングフリップ角は180°となる。例示的な実施形態において、k空間602は、256×256×50ボクセルのスキャン体積に対応し得る。当然、ハーフフーリエ法が使用されるため、実際にはエコー情報からデータの一部のみ収集すればよい。
図7は、実施形態に係るk空間及びパルスシーケンスを示す図(3)である。図7には、k空間702を概略的に示し、図6に示したk空間602を一側面から見た図に対応する。図7に示すリフォーカシングパルスの縦方向の大きさは、リフォーカシングフリップ角の大きさに対応する。k空間702において、複数のFASEパルスシーケンス708におけるリフォーカシングパルス710に対応する実効TEエコーが記録される。各FASEパルスシーケンス712,714,716,718,720からの実効TEに対応するエコー信号はk空間にマッピングされる。図示の通り、3D FASEパルスシーケンスが形成する可変リフォーカシングフリップ角のパターンにより、高いSNRを有する実効TEエコー710が得られ、さらに他部ではより小さなリフォーカシングフリップ角によりSARを低減させられる。図示の実施形態では、k空間の中心に最も近い3ラインに対応するリフォーカシングパルス716,718,720に対し、同一の大きなリフォーカシングフリップ角が割り当てられ、k空間の中心からより離れたラインは、より小さい角度のリフォーカシングパルス712,714を有する。大きな(又は最大の)リフォーカシングフリップ角がk空間の中心にあるパターンであれば、その他の設定も可能である。
図8は、実施形態に係る可変TR間隔の一例を示す図(1)である。図8には、可変TR間隔のパターン802を示す。可変TR間隔パターン802は、実施形態によっては、リフォーカシングパルスの可変リフォーカシングフリップ角と組み合わせて使用される。パターン802には、(例えばk空間702に示される)スライス1〜5が順次配置される3Dk空間において、k空間の中心から最も離れたスライス1及び5のそれぞれにおいて、FSEパルスシーケンスは2RRのTR間隔で収集され、スライス2〜4の場合、3RRのTR間隔で収集される。したがって、パターン802は、2段階可変TR間隔構造となっている。同パターンにおいて、k空間の中心に近いスライスは長いTR間隔が割り当てられ、その他スライスには短いTR間隔が割り当てられる。なお、RRの間隔は、例えば、トリガ(心電トリガ)の間隔に対応する。すなわち、本実施形態に係る繰り返し時間(TR間隔)は、シーケンスにおいて繰り返し印加される励起パルスの間隔であり、被検体から得られる心電信号の心電トリガの間隔の整数倍に設定される。
イメージングされる各スライスの収集(例えばエンコード)に要する時間の合計が総スキャン時間となる。したがって、全てのスライスエンコーディングでTRが同一の場合、1スライスの総TR間隔とスライス数との積に基づいて総スキャン時間が定まる。異なるTR間隔を持つスライスがある場合、総スキャン時間は全スライスのTR間隔の合計として表される。
実施形態によっては、使用されるTR間隔は、R波トリガからの遅延、3D FASE読出し、次のR波トリガまでの回復時間を含む。例示的な実施形態によっては、選択されたスライスのTR間隔は、FASE読出し取得に要する時間を変更することで変更される。具体的には、実施形態によっては、連続するリフォーカシングパルスの間隔(したがって読出し傾斜の間隔)が短いTR間隔ではエコー信号による縦磁化回復が低く、長いTR間隔ではエコー信号による縦磁化回復が高くなるように変更される。したがって、読出し間のTR間隔が3RRの場合、2RRの場合に比べて縦磁化(Mz)回復が上がるが、スキャン時間が延びてしまう。
SARはTRに関連しており、TRが短いほどFBIスキャンにおいてSARが高くなる。FBIにおいて、総スキャン時間はスライスエンコーディングの回数と1総スライスエンコーディングのTRに比例する。従来のFBIにおいて、各スライスエンコーディングは固定されたTR内で実行される。上述のように、TRはR波トリガからの遅延、3D FASE読出し、次のR波トリガまでの回復時間を含む。従来例において、TRは、例えば、2RR、3RR、4RR、5RRのうちいずれかの長さに固定されており、一旦TRの値が設定されると、スキャニングを通じて固定されたTRが使用される。通常、RR数の下限はSARにより定まり、心拍数の高い被検体にはRRを増やしてSARを下げる必要がある。
図9は、実施形態に係る可変TR間隔の一例を示す図(2)である。図9には、実施形態によっては使用可能な、別の可変TR間隔のパターン902を示す。図9に示すパターン902は、複数のスライスのエンコーディングに3以上の異なるTR間隔が使用される多段階パターンである。図9に示すように、パターン902では、SE方向にけるk空間中心に最も近い1又は複数のスライスに最長のTR間隔(例えば5RR)が割り当てられる。その他スライスは、k空間の中心から周囲又は外縁に向かうにつれTR間隔が増加しないように構成される。図9に示すように、当該パターン902は、(k空間の)最も中心にあるスライスが最長のTR間隔を有し、外縁のスライスが最短のTR間隔を持つように構成されているが、各TRが連続した所定数のスライスに割り当てられるようにしてもよい。例えば、図9に示す「繰り返し番号」列によると、番号dとなる(k空間の)最も中心のスライスには5RRのTRが割り当てられており、番号cとなる5RRスライスの両側にある次のスライスには4RRのTRが割り当てられており、番号bとなる4RRスライスの次のスライスには3RRが割り当てられており、両端の最後のスライス(番号a)には2RRが割り当てられている。位相エンコード方向やスライスエンコード方向において、TR間隔を変更してもよい。
言い換えると、シーケンスは、中心を含む領域の繰り返し時間である第1の時間間隔と、k空間の周辺を含む領域の繰り返し時間である第2の時間間隔と、中心を含む領域と周辺を含む領域との間の領域の繰り返し時間である第3の時間間隔とを有し、第3の時間間隔が第1の時間間隔より短く、第2の時間間隔より長い。
なお、図8及び図9は、一例に過ぎない。例えば、シーケンスは、位相エンコード方向及びスライスエンコード方向の少なくとも一方において、中心を含む領域の繰り返し時間が他の領域の繰り返し時間よりも長い場合であってもよい。また、例えば、シーケンスは、位相エンコード方向及びスライスエンコード方向において、中心を含む領域の繰り返し時間が他の領域の繰り返し時間よりも長い場合であってもよい。
(その他の実施形態)
上述した実施形態以外にも、種々の異なる形態にて実施されてもよい。
例えば、上記の説明では、実施形態がFSE(或いはFASE)に適用される場合を説明したが、これに限定されるものではない。例えば、実施形態は、可変リフォーカシングフリップ角を用いた撮像であれば、比吸収率の上昇を抑えつつ信号雑音比を効果的に高めるために、位相エンコード方向及びスライスエンコード方向の少なくとも一方において、k空間の中心を含む領域で印加されるリフォーカシングパルスのフリップ角が他の領域で印加されるリフォーカシングパルスのフリップ角よりも大きいシーケンスを実行してもよい。
(bSSFPへの適用)
また、上記の実施形態は、bSSFP(balanced Steady-State Free Precession)法に適用されても良い。例えば、シーケンス制御部としてのMRIシーケンスコントローラ30は、位相エンコード方向及びスライスエンコード方向の少なくとも一方において、k空間の中心を含む領域で印加される高周波パルスのα角が他の領域で印加される高周波パルスのα角よりも大きいシーケンスを実行する。
なお、以下に説明するbSSFPが適用されるアプリケーションの例としては、Time−SLIP(Time-Spatial Labeling Inversion Pulse)法やFSD(Flow-Spoiled Dephasing)法が挙げられる。ここで、Time−SLIP法とは、撮像領域とは独立に設定される標識化領域において組織の縦磁化成分を反転させる反転パルスを流体に印加し、所定時間が経過した後に関心領域のデータを収集することで、関心領域内に流入した流体を選択的に描出する撮像法である。また、FSD法とは、例えば90°−180°−90°のRFパルス等、0でない一次傾斜磁場モーメント(m1)を有する傾斜磁場プリパルスを用いて流れるスピンの位相を分散(ディフェーズ)させることで、DA(Dark Artery)画像(動脈が暗い画像)を得る撮像法である。すなわち、上記の実施形態がbSSFP法に適用される場合、Time−SLIP法による反転パルスやFSD法による傾斜磁場プリパルスが印加された後に、以下に説明するbSSFP法によってエコー信号の収集が行われる。
例えば、シーケンス制御部としてのMRIシーケンスコントローラ30は、k空間のサンプリングをセントリックに行う場合には、ランプアップ(Ramp up)を行った後にサンプリングを開始し、k空間のサンプリングをシーケンシャルに行う場合には、ランプアップを行わずにサンプリングを開始する。
図10は、その他の実施形態に係るk空間及びパルスシーケンスを示す図(1)である。図10には、2Dのk空間1001のサンプリングをセントリック(インターリーブ)に行う場合のパルスシーケンス1002を例示する。図10に示すRFパルスの縦方向の大きさは、α角の大きさに対応する。また、bSSFPにおいては、各RFパルスの間隔(τ)は、一定である。
図10に示すように、k空間1001のサンプリングをセントリックに行う場合には、k空間1001の中心付近から外縁に向かって順番にサンプリングが行われる。この場合、サンプリングが開始された時点で最も大きなα角のRFパルスが印加されるのが好ましい。このため、MRIシーケンスコントローラ30は、パルスシーケンス1002において、まず、RFパルスによる励起を繰り返し行うランプアップを行う(図10の「Ramp up」参照)。図10の例では、繰り返し印加されるRFパルスのα角が徐々に大きくなるように行われる。その後、MRIシーケンスコントローラ30は、bSSFPによるサンプリングを開始する。ここで、サンプリング1〜5は、k空間1001の位相エンコード方向における中心付近に対応している。このため、MRIシーケンスコントローラ30は、サンプリング1〜5において印加されるRFパルスのα角を最も大きし、サンプリング5より後のサンプリングにおいて印加されるRFパルスのα角を徐々に小さくする。
なお、図10は一例に過ぎない。例えば、MRIシーケンスコントローラ30は、必ずしもサンプリングが開始された時点で最も大きなα角のRFパルスを印加しなくてもよい。つまり、MRIシーケンスコントローラ30は、十分なSNRが確保できる範囲内の大きさのα角であれば、最大のα角でなくともサンプリングを開始してよい。また、図10では、k空間の全体に亘ってデータを収集する例を説明したが、実施形態はこれに限定されるものではない。例えば、k空間の対称性を利用してk空間の半分程度を充填するハーフフーリエ法によってデータが収集される場合にも、実施形態は同様に適用可能である。
図11は、その他の実施形態に係るk空間及びパルスシーケンスを示す図(2)である。図11には、2Dのk空間1101のサンプリングをシーケンシャルに行う場合のパルスシーケンス1102を例示する。図11に示すRFパルスの縦方向の大きさは、α角の大きさに対応する。
図11に示すように、k空間1101のサンプリングをシーケンシャルに行う場合には、k空間1101の一方の外縁からもう一方の外縁に向かって順番にサンプリングが行われる。この場合、中心付近のサンプリングが行われる時点で最も大きなα角のRFパルスが印加されるのが好ましい。このため、MRIシーケンスコントローラ30は、パルスシーケンス1102において、サンプリング1,2,3,4,5・・・と、k空間1101の外縁から順番にbSSFPによるサンプリングを開始する。そして、MRIシーケンスコントローラ30は、外縁から中心付近に向かうにつれてα角を大きくし、中心付近で一定以上のα角を保ち、中心付近から再び外縁に向かうにつれてα角を小さくする。
なお、図11は一例に過ぎない。例えば、図11に示す例では、ランプアップを行わずにサンプリングを開始する場合を説明したが、例えば、MRIシーケンスコントローラ30は、図10に示したように、ランプアップを行ってからサンプリングを開始してもよい。また、ハーフフーリエ法によってデータが収集される場合にも、同様に適用可能である。
図12は、その他の実施形態に係るk空間及びパルスシーケンスを示す図(3)である。図12には、3Dのk空間1201のサンプリングをシーケンシャルに行う場合のパルスシーケンス1202,1204,1206,1208,1210を例示する。ここで、各パルスシーケンス1202,1204,1206,1208,1210は、k空間1201の各スライス1,2,3,4,5で行われるパルスシーケンスに対応する。図12に示すRFパルスの縦方向の大きさは、α角の大きさに対応する。
図12に示すように、k空間1201のサンプリングをシーケンシャルに行う場合には、k空間1201の一方の外縁のスライスからもう一方の外縁のスライスに向かって順番にサンプリングが行われる。例えば、MRIシーケンスコントローラ30は、中心付近にあるスライス3において、最も大きいα角のRFパルスを印加するパルスシーケンス1206を行う。また、MRIシーケンスコントローラ30は、外縁にあるスライス1,5において、最も小さいα角のRFパルスを印加するパルスシーケンス1202,1210を行う。また、MRIシーケンスコントローラ30は、中心付近と外縁との間にあるスライス2,4において、スライス3より小さく、スライス1,5より大きいα角のRFパルスを印加するパルスシーケンス1204,1208を行う。
なお、図12は一例に過ぎない。例えば、MRIシーケンスコントローラ30は、k空間1201のサンプリングをシーケンシャルに行う場合においても、ランプアップを行ってからサンプリングを開始してもよい。
また、例えば、MRIシーケンスコントローラ30は、k空間1201のサンプリングをセントリックに行う場合においても、図12に示したパルスシーケンスと同様のパルスシーケンスでサンプリングを行ってもよい。つまり、MRIシーケンスコントローラ30は、SE方向にインターリーブにサンプリングを行う場合には、例えば、スライス3,2,4,1,5の順に、図12に示した大きさのα角のRFパルスを用いてサンプリングを行う。また、PE方向にインターリーブにサンプリングを行う場合には、例えば、図10に示したパルスシーケンスと同様のパルスシーケンスでサンプリングを行う。この場合、MRIシーケンスコントローラ30は、ランプアップを行ってからサンプリングを開始するのが好ましい。
これにより、その他の実施形態に係るMRIシステムは、bSSFPにおいてもサンプリングを行うことができる。
以上説明した少なくともひとつの実施形態によれば、比吸収率の上昇を抑えつつ信号雑音比を効果的に高めることができる。
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
22 MRIシステムコントローラ
30 MRIシーケンスコントローラ

Claims (9)

  1. スライスエンコード方向において、3次元k空間の中心を含むハーフフーリエスライスで印加されるリフォーカシングパルスのフリップ角が他のハーフフーリエスライスで印加されるリフォーカシングパルスのフリップ角よりも大きくなるようにFASE(Fast Advanced Spin Echo)シーケンスを実行するシーケンス制御部
    を備え
    前記シーケンス制御部は、
    第1のスライスに印加する前記リフォーカシングパルスを第1のフリップ角で、
    前記第1のスライスより前記k空間の中心に近い第2のスライスに印加する前記リフォーカシングパルスであって、実効TEパルス以外のリフォーカシングパルスを、前記第1のフリップ角より大きな第2のフリップ角で、
    前記第2のスライスより前記k空間の中心に近い第3のスライスに印加する前記リフォーカシングパルスであって、実効TEパルス及び前記実効TEパルスの両側のパルス以外のリフォーカシングパルスを、前記第2のフリップ角より大きな第3のフリップ角で、
    前記第2のスライスに印加する前記リフォーカシングパルスであって、前記実効TEパルスであるリフォーカシングパルスを、前記第2のフリップ角より大きな第4のフリップ角で、
    前記第3のスライスに印加する前記リフォーカシングパルスであって、前記実効TEパルスの両隣のパルスを、前記第3のフリップ角より大きな第5のフリップ角で、
    前記第3のスライスに印加する前記リフォーカシングパルスであって、前記実効TEパルスであるリフォーカシングパルスを、前記第4のフリップ角及び前記第5のフリップ角より大きな第6のフリップ角で印加する、磁気共鳴イメージング装置。
  2. 前記FASEシーケンスは、3次元k空間の位相エンコード方向及びスライスエンコード方向において、前記中心を含む領域で印加されるリフォーカシングパルスのフリップ角が他の領域で印加されるリフォーカシングパルスのフリップ角よりも大きい、請求項1に記載の磁気共鳴イメージング装置。
  3. 前記FASEシーケンスは、前記中心に対してフリップ角の大きさが略対称のリフォーカシングパルスを含む、請求項1又は2に記載の磁気共鳴イメージング装置。
  4. 前記FASEシーケンスは、位相エンコード方向において、パルストレインの開始時に対応する、前記3次元k空間の前記中心を含む領域外の領域及び前記中心を含む領域で印加されるリフォーカシングパルスのフリップ角が他の領域で印加されるリフォーカシングパルスのフリップ角よりも大きい、請求項1又は2に記載の磁気共鳴イメージング装置。
  5. 前記FASEシーケンスによって収集されたエコー信号に基づいて、出力画像を生成する生成部を更に備え、
    前記シーケンス制御部は、第1の心位相及び第2の心位相において、前記FASEシーケンスをそれぞれ実行し、
    前記生成部は、前記第1の心位相において収集されたエコー信号から得られる第1の画像と、前記第2の心位相において収集されたエコー信号から得られる第2の画像との差分に基づいて、前記出力画像を生成する、請求項1〜4のいずれか一つに記載の磁気共鳴イメージング装置。
  6. 前記FASEシーケンスは、位相エンコード方向及び前記スライスエンコード方向の少なくとも一方において、前記中心を含む領域の繰り返し時間が他の領域の繰り返し時間よりも長い、請求項1〜5のいずれか一つに記載の磁気共鳴イメージング装置。
  7. 前記FASEシーケンスは、前記位相エンコード方向及び前記スライスエンコード方向において、前記中心を含む領域の繰り返し時間が他の領域の繰り返し時間よりも長い、請求項6に記載の磁気共鳴イメージング装置。
  8. 前記繰り返し時間は、前記FASEシーケンスにおいて繰り返し印加される励起パルスの間隔であり、被検体から得られる心電信号の心電トリガの間隔の整数倍に設定される、請求項6又は7に記載の磁気共鳴イメージング装置。
  9. 前記FASEシーケンスは、前記中心を含む領域の繰り返し時間である第1の時間間隔と、前記k空間の周辺を含む領域の繰り返し時間である第2の時間間隔と、前記中心を含む領域と前記周辺を含む領域との間の領域の繰り返し時間である第3の時間間隔とを有し、前記第3の時間間隔が前記第1の時間間隔より短く、前記第2の時間間隔より長い、請求項6〜8のいずれか一つに記載の磁気共鳴イメージング装置。
JP2015125009A 2014-06-23 2015-06-22 磁気共鳴イメージング装置 Active JP6656826B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US14/312,121 US9983286B2 (en) 2014-06-23 2014-06-23 SAR reduction in fast advanced spin echo (FASE) or single-shot fast spin echo (SS-FSE) imaging
US14/312,121 2014-06-23

Publications (2)

Publication Number Publication Date
JP2016007539A JP2016007539A (ja) 2016-01-18
JP6656826B2 true JP6656826B2 (ja) 2020-03-04

Family

ID=54869428

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015125009A Active JP6656826B2 (ja) 2014-06-23 2015-06-22 磁気共鳴イメージング装置

Country Status (2)

Country Link
US (1) US9983286B2 (ja)
JP (1) JP6656826B2 (ja)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9720064B2 (en) * 2014-08-18 2017-08-01 Toshiba Medical Systems Corporation Variable TR (vTR) function in fresh blood imaging (FBI)
US10598750B2 (en) * 2015-04-02 2020-03-24 The General Hospital Corporation 3D balanced EPI magnetic resonance fingerprinting
US10145926B2 (en) * 2015-11-30 2018-12-04 Siemens Healthcare Gmbh Phase-cycled steady-state free precession imaging
CN105997075B (zh) * 2016-04-28 2018-05-04 上海联影医疗科技有限公司 一种磁共振数据采集方法及磁共振成像系统
WO2018085970A1 (zh) * 2016-11-08 2018-05-17 中国科学院深圳先进技术研究院 磁共振血管壁成像方法和设备
JP6956509B2 (ja) * 2017-04-21 2021-11-02 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
US11047942B2 (en) 2017-10-24 2021-06-29 University Of Cincinnati Magnetic resonance imaging method and system with optimal variable flip angles
EP3480616A1 (de) * 2017-11-06 2019-05-08 Siemens Healthcare GmbH Magnetresonanz-system und verfahren zum übertragen von steuerdaten an eine digitale steuereinrichtung eines magnetresonanz-systems
JP7223580B2 (ja) * 2018-01-26 2023-02-16 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
EP3524992B1 (de) * 2018-02-13 2021-04-07 Siemens Healthcare GmbH Verfahren zum betrieb einer magnetresonanzeinrichtung, magnetresonanzeinrichtung, computerprogramm und elektronisch lesbarer datenträger
US10746830B2 (en) 2018-08-28 2020-08-18 General Electric Company Systems and methods for hybrid slice encoding in three-dimensional magnetic resonance imaging
US11353527B2 (en) * 2019-07-19 2022-06-07 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for waveform determination in magnetic resonance imaging
CN110361682A (zh) * 2019-07-31 2019-10-22 张忠铭 快速确定磁共振射频激发翻转角度对应能量自动校准方法
CN111157932B (zh) * 2020-01-02 2022-08-30 华东师范大学 一种快速自旋回波脉冲序列中射频脉冲的优化方法
JP7374822B2 (ja) * 2020-03-09 2023-11-07 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング方法及び装置
US11686802B2 (en) * 2021-11-12 2023-06-27 Maier Stephan Method and magnetic resonance apparatus for diffusion image acquisition with motion offsetting and navigation-dependent segmentation

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3405813B2 (ja) * 1994-05-31 2003-05-12 株式会社東芝 磁気共鳴イメージング方法および磁気共鳴イメージング装置
JPH08117206A (ja) * 1994-10-28 1996-05-14 Toshiba Corp 磁気共鳴イメージング装置
JP2000041969A (ja) * 1998-07-28 2000-02-15 Toshiba Corp Mri装置およびmrイメージング方法
JP4225648B2 (ja) * 1999-09-08 2009-02-18 株式会社東芝 Mri装置
DE10219528A1 (de) * 2002-05-02 2003-11-13 Philips Intellectual Property Schnelles Kernresonanz-Bildgebungsverfahren mit Gradienten-Echos
JP4309632B2 (ja) * 2002-10-08 2009-08-05 株式会社東芝 磁気共鳴イメージング装置
US6885193B2 (en) * 2002-11-14 2005-04-26 Koninklijke Philips Electronics, N.V. Magnetization primer sequence for balanced steady state free precision imaging
US20120296199A1 (en) * 2011-03-21 2012-11-22 New York University Apparatus and Method of Non-Contrast Magnetic Resonance Angiography of Abdominal and Pelvic Arteries
JP5960134B2 (ja) * 2011-06-09 2016-08-02 株式会社日立製作所 磁気共鳴イメージング装置及びフリップ角決定方法
US9414765B2 (en) * 2012-03-08 2016-08-16 Toshiba Medical Systems Corporation Fresh blood imaging (FBI) with independently controlled MRI data acquisition parameters for diastolic and systolic MRI acquisition sub-sequences
JP6084430B2 (ja) * 2012-10-30 2017-02-22 株式会社日立製作所 再収束rfパルスのフリップ角制御法及び磁気共鳴イメージング装置
US10180481B2 (en) * 2013-03-15 2019-01-15 The Regents Of The University Of California System and method of magnetic resonance imaging using variable flip angle balanced steady-state free precession (VFA-BSSFP)
US9720064B2 (en) * 2014-08-18 2017-08-01 Toshiba Medical Systems Corporation Variable TR (vTR) function in fresh blood imaging (FBI)

Also Published As

Publication number Publication date
US9983286B2 (en) 2018-05-29
JP2016007539A (ja) 2016-01-18
US20150369891A1 (en) 2015-12-24

Similar Documents

Publication Publication Date Title
JP6656826B2 (ja) 磁気共鳴イメージング装置
Feng et al. 5D whole‐heart sparse MRI
US9839366B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
US9700220B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
US11199604B2 (en) MRI involving the acquisition of an angiography weighted image and of a perfusion weighted image
JP6625353B2 (ja) 磁気共鳴イメージング装置
JP5269342B2 (ja) 磁気共鳴イメージング装置および磁気共鳴イメージング装置における撮影条件設定方法
US20100249574A1 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
Pang et al. High efficiency coronary MR angiography with nonrigid cardiac motion correction
JP6077090B2 (ja) 磁気共鳴イメージング装置
JP6510273B2 (ja) 磁気共鳴イメージング装置、磁気共鳴イメージング方法及び磁気共鳴イメージングプログラム
JP6529774B2 (ja) 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
JP2015208679A (ja) 磁気共鳴イメージング装置、方法及びプログラム
JP6529807B2 (ja) 磁気共鳴イメージング装置、方法及びプログラム
Dietrich et al. Motion‐compensated fat‐water imaging for 3D cardiac MRI at ultra‐high fields
JP4230875B2 (ja) 磁気共鳴イメージング装置
JP6154161B2 (ja) 磁気共鳴イメージング装置
JP6348909B2 (ja) 磁気共鳴システム及び磁気共鳴方法
KR102342813B1 (ko) 볼륨 선택적 3차원 자기 공명 영상 생성 장치 및 방법
Kadbi et al. Assessment of flow and hemodynamics in the carotid artery using a reduced TE 4D flow spiral phase-contrast MRI
Schoormans et al. 3D black-blood DCE-MRI using radial stack-of-stars acquisition and CS reconstruction: application in carotid and femoral arteries
JP2023008743A (ja) 磁気共鳴イメージング装置
Yang et al. Single-shot spiral first-pass perfusion imaging: full heart coverage with high temporal resolution
Plein et al. Dealing with Cardiac Motion: How Do We Image the Beating Heart?

Legal Events

Date Code Title Description
RD01 Notification of change of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7421

Effective date: 20151102

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20160513

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20160929

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20161021

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180619

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20190124

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190219

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190419

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190903

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20191029

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20200107

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200205

R150 Certificate of patent or registration of utility model

Ref document number: 6656826

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150