JP2021090760A - スタックオブスター獲得を使用したmr像形成 - Google Patents

スタックオブスター獲得を使用したmr像形成 Download PDF

Info

Publication number
JP2021090760A
JP2021090760A JP2021001435A JP2021001435A JP2021090760A JP 2021090760 A JP2021090760 A JP 2021090760A JP 2021001435 A JP2021001435 A JP 2021001435A JP 2021001435 A JP2021001435 A JP 2021001435A JP 2021090760 A JP2021090760 A JP 2021090760A
Authority
JP
Japan
Prior art keywords
space
radial
image
signal
density
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.)
Pending
Application number
JP2021001435A
Other languages
English (en)
Inventor
ガブリエル マリアンヌ べック
Marianne Beck Gabriele
ガブリエル マリアンヌ べック
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2021090760A publication Critical patent/JP2021090760A/ja
Pending legal-status Critical Current

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/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/4824MR 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 non-Cartesian trajectory
    • G01R33/4826MR 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 non-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/5601Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution involving use of a contrast agent for contrast manipulation, e.g. a paramagnetic, super-paramagnetic, ferromagnetic or hyperpolarised contrast agent
    • 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/5611Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
    • 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/56308Characterization of motion or flow; Dynamic imaging

Landscapes

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

Abstract

【課題】低いレベルの筋状アーチファクトを伴うスタックオブスター獲得スキームを使用して高速MR像形成を可能にする方法を提供する。【解決手段】少なくとも1つのRFパルスとスイッチングされる磁場勾配との像形成シーケンスを物体に施すステップと、スタックオブスタースキームに従ってMR信号を獲得するステップであって、MR信号が、スライス方向に沿った異なる位置に配置された複数の平行スライスから、放射状k空間プロファイルS1〜S12として獲得されk空間プロファイルの放射状密度が、スライス位置の関数として変化し、放射状密度が、より中心のk空間位置においてより高く、より周辺のk空間位置においてより低く、k空間プロファイルが、より周辺のk空間位置におけるスライスからより、より中心の位置におけるスライスから、より高い時間密度で獲得される、ステップと、MR信号からMR像を再構成するステップとを有する、方法。【選択図】図4

Description

本発明は、磁気共鳴(MR:magnetic resonance)像形成の分野に関する。本発明は、物体のMR像形成の方法に関係する。本発明は、さらに、MRデバイス上で動作するMRデバイス及びコンピュータプログラムに関する。
二次元又は三次元像を形成するために磁場と核スピンとの間における相互作用を使用する像形成MR方法は、軟組織の像形成の場合において、多くの点で他の像形成方法より優れており、電離放射を必要とせず、通常、非侵襲性であるので、特に医療診断の分野において現在、広く使用されている。
概して、MR方法によると、物体、例えば、検査される患者の体は、強く一様な磁場内に配置され、同時に磁場の方向が、測定の基準となる座標系の軸(通常はz軸)を規定する。磁場は、規定された周波数(いわゆる、ラーモア周波数、又はMR周波数)の交番電磁場(RF場)の印加により、励起(スピン共鳴)され得る磁場強度に応じて、個々の核スピンに対して異なるエネルギーレベルを生成する。巨視的な視点によると、個々の核スピンの分布は、適切な周波数の電磁パルス(RFパルス)の印加により平衡状態から外れることにより磁化がz軸周りの歳差運動をもたらし得る全体の磁化を生成する。歳差運動は、コーンの表面を規定し、コーンの開口の角度は、フリップ角と呼ばれる。フリップ角の大きさは、印加される電磁パルスの強度及び持続期間に依存する。いわゆる90°パルスの場合、スピンはz軸から、横断面(フリップ角90°)に外れる。
RFパルスの終了後、磁化が緩和して元の平衡状態に戻り、その際、z方向の磁化が第1の時定数T(スピン格子又は縦緩和時間)で再度構築され、z方向に直交する方向の磁化が第2の時定数T(スピン−スピン又は横緩和時間)で緩和する。磁化の変化は、z軸に直交する方向において磁化の変化が測定されるような手法で、MRデバイスの検査ボリューム内に配置及び配向された受信RFコイルにより検出され得る。例えば、90°パルスの印加後、横磁化の減衰の後に、同位相の秩序状態から、すべての位相角が一様に分布(ディフェーズ)した状態への(局所磁場の不均一さにより誘導される)核スピンの遷移が追随する。ディフェーズは、再合焦パルス(例えば、180°パルス)により補償され得る。これが、受信コイル内にエコー信号(スピンエコー)を生成する。
体内における空間分解能を実現するために、3つの主軸に沿って延びる一定の磁場勾配が、一様な磁場に重ね合わされて、スピン共鳴周波数の線形空間依存性をもたらす。受信コイル内において拾われた信号は、したがって、体内の異なる位置に関係付けられ得る異なる周波数の成分を含む。受信コイルを介して取得された信号データは、空間周波数領域に対応し、k空間データと呼ばれる。k空間データは、通常、異なる位相符号化により獲得された複数のラインを含む。各ラインは、多くのサンプルを収集することによりデジタル化される。k空間データのセットが、像再構成アルゴリズムによりMR像に変換される。
知られたいわゆる三次元(3D)スタックオブスター(stack-of-stars)獲得スキーム(例えば、WO2013/159044A1を参照)において、複数の空間的に非選択的又はスラブ選択的なRF励起が適用され、各々の後に、1つ又は複数のMR信号(例えば、勾配エコー信号)の獲得が続き、各MR信号は、k空間プロファイルを表す。MR信号は、複数の平行スライスから放射状k空間プロファイルとして獲得される。スライスは、スライス方向に沿った異なる位置に配置される。中心(k=k=0)の周りで回転する半径方向「スポーク」に沿って各単一スライス内においてMR信号が獲得されながら、スライス方向(例えば、k方向)において、標準的なデカルト位相符号化が実施される。これは、積み重ねられたディスク(「スタックオブスター」)から構成される円柱k空間カバレッジをもたらす。技術的に、これは、スライスの面内方向における磁場勾配を同時に生成することにより、及びそれらの振幅を変調することにより実現される。k空間プロファイル獲得ステップの時間的な順序を選択するために、異なるスキームが使用され得る。例えば、異なる角度位置でk空間プロファイルが獲得される前に、スライス方向に沿ったすべての位相符号化ステップが連続的に獲得され得る。このことが、デカルトサンプリングの期間が短く維持されることを確実にして、スライスの積み重ね内における高いデータ一貫性をもたらし、スタックオブスターアプローチのための半径方向サンプリングの全体的な動きロバスト性を維持する。デカルト位相符号化ステップは、中心スライスからk空間周辺(中心外)に、又は、−kz,maxから+kz,maxに線形順序で実施される。角度的な順序付けの場合、像形成シーケンスは、複数のインターリーブを伴う等角距離サンプリング、又は、いわゆる黄金角度スキームを使用し得る。等距離スキームにおいて、角距離は、ΔΦ=180°/ntotalに従って計算され、ntotalは、スポークの総数である。インターリーブがk空間内における時間的なコヒーレンスを低減するので、複数のインターリーブ(又は、「回転」)を使用してスポークを獲得することが有益である。したがって、動きの非一貫性がk空間内に広げられ、アーチファクトが減衰される。黄金角度スキームにおいて、k空間プロファイルの角度は、各時点において、180°に黄金比を乗算した値に対応するΔΦ=111.25°ずつインクリメントされる。したがって、続いてサンプリングされたスポークは、以前にサンプリングされたスポークのセット内における最大ギャップを埋めながら、常に補足情報を追加する。その結果、任意の連続的な獲得されたスポークのセットが、ほぼ一様にk空間をカバーし、このことが、時間的なサブフレームの再構成を可能にし、黄金角度スキームを動的像形成研究に十分に適したものにする。
上述の3D半径方向スタックオブスタースキームは、高い動きロバスト性、良性エイリアシングアーチファクト、及びk空間中心の連続更新など、臨床MR像形成のためのいくつかの有望な利点を提案する。しかし、これらの利点にもかかわらず、筋状アーチファクトが多くの場合に問題となる。より高く要求される半径方向サンプリング密度と、その結果としての、より長い走査期間が導かれる。
D.C.PetersらによるISMRM−アブストラクト2007 p.207「Radial undersampling that is variable in k」は、k空間内におけるスタックオブスターサンプリングアプローチを開示する。k次元において変化する半径方向アンダーサンプリングパターンが使用されて、中心のk空間に対してより多くの射影を獲得する。
上述の事項により、改善されたMR像形成技術が必要とされることが容易に理解される。したがって、筋状アーチファクトのレベルを低減したスタックオブスター獲得スキームを使用した高速MR像形成を可能にすることが、本発明の目的である。
本発明によると、MRデバイスの検査ボリューム内に位置する物体のMR像形成の方法が開示される。本方法は、
少なくとも1つのRFパルスとスイッチングされる磁場勾配との像形成シーケンスを物体に施すステップと、
スタックオブスタースキームに従ってMR信号を獲得するステップであって、MR信号が、スライス方向に沿った異なる位置に配置された複数の平行スライスから複数の異なる角度位置で放射状k空間プロファイルとして獲得され、k空間プロファイルの放射状密度がスライス位置の関数として変化し、放射状密度が、より中心のk空間位置においてより高く、より周辺のk空間位置においてより低く、k空間プロファイルが、より周辺のk空間位置におけるスライスからより、より中心のk空間位置におけるスライスから、より高い時間密度で、すなわち、単位時間当たりに獲得されたk空間プロファイルの数で獲得される、ステップと、
MR信号からMR像を再構成するステップと、
を有する。
本発明は、修正されたスタックオブスター獲得技法を適用することを提案する。修正は、放射状k空間サンプリング密度(すなわち、スライス当たりの異なるように配向された放射状k空間プロファイルの数)がスライスごとに変更されるというものである。k空間プロファイルの放射状密度は、k空間中心(k=0)により近く位置するスライス内においてより高く、k空間プロファイルの放射状密度は、k空間中心からより遠くに位置するスライス内においてより低い。
k空間中心付近において放射状k空間サンプリングの密度が高められることは、筋状アーチファクトのレベルを低減するとともに、総走査期間が最小に維持され得る。
k空間の中心領域の定期的な更新に起因して、放射状k空間サンプリングが目的の領域内における造影剤の取り込みを追跡することを可能にするので、放射状k空間サンプリングは造影MR像形成の場合において興味深いことが一般的に知られている。それにもかかわらず、コントラスト操作のためにk空間で重み付けされた像コントラスト法(Songら、Magnetic Resonance in Medicine、vol.44、2004、825〜832ページを参照)と組み合わせた場合でさえ、従来の放射状スタックオブスター獲得から再構成された時系列に沿ったMR像は、不鮮明な造影プロファイルを示すので、造影剤が特定の像ボクセルに達する時点が明確に識別不能である。これは、本発明の方法により改善され得る。
好ましい実施形態において、本発明は、より周辺のk空間位置におけるスライスからより、より中心のk空間位置におけるスライスから、放射状k空間プロファイルがより高い時間密度で獲得される獲得順序を提案する。言い換えると、中心のナイキスト楕円体内の放射状k空間プロファイルが、近い時点において獲得される一方で、周辺の放射状k空間プロファイル(すなわち、ナイキスト楕円体の外部)は、より遠い時点においてサンプリングされる。k空間からより多くのサンプルが獲得されるので、内部におけるサンプリングされたk空間プロファイルの(半径方向)密度が等しいか又はより高い中心のナイキスト楕円体の寸法に対し、視界のためのナイキスト基準がk空間内において半径方向外向きに徐々に上昇し、周辺のk空間内におけるサンプリングされたk空間プロファイルの(放射状)密度は視界のためのナイキスト基準以下となることに留意されたい。ナイキスト楕円体の実際の寸法は、再構成されたMR像の空間分解能を決定する。
この手法は、造影MR像形成を改善する。中心のk空間楕円体において高められた時間密度は、造影剤が特定の像ボクセルに達するコントラスト時点が捕捉されることを確実にする。
各MR像は、次に、ナイキスト楕円体内において、(より遠い時点においてサンプリングされた)周辺のk空間領域から獲得されたk空間プロファイルと組み合わせて、それぞれ最も近い時点において獲得されたk空間プロファイルのみを使用して再構成される。これは、互いに時間的に最も近くで獲得されたMR信号部分が、MR像を再構成するために中心の−完全にサンプリングされた、又はオーバーサンプリングされた−楕円体k空間領域内において選択されることを意味する。周辺のk空間MR楕円体に対しては、外向きの放射状k空間位置から延びた、k空間内における部分的な放射状スポークのみが選択される。これらの部分的なスポークは、周辺のk空間MR楕円体へと外向きになるにつれ、これらの部分的なスポークは、より高い角密度をもつ。
このアプローチによれば、時間的に近くで獲得されたMR信号のみがk空間中心内におけるMR像再構成のために使用されるので、造影剤が特定の像ボクセルに達する時点がより明確に識別され得る。したがって、測定された造影プロファイルの不鮮明さが低減する。
楕円体の外部のk空間位置におけるMR信号は、MR像を再構成する前に除去される。これは、k空間内における楕円体「シャッター」を構成する。これは、筋状アーチファクトを効果的に除去する。さらに有益な実施態様において、k空間の中心領域は、(ナイキスト基準より高いサンプリング密度で)オーバーサンプリングされ、このことが、MR信号の獲得を、動きによるアーチファクトに対してよりロバストなものとする。
MR像の再構成に対するより改良された実施態様において、周辺のk空間からのMR信号は、外向きの放射状k空間位置から延びたk空間内の部分的な放射状スポークに沿って選択され、これらの部分的なスポークは、それらの外向きの放射状k空間位置がk空間内において、より半径方向外向きであるので、k空間内においてより高い角密度をもつ。部分的なスポークは、スタックオブスターのk軸から半径方向に分離されたk平面の各々内にある、それらの外向きの放射状k空間位置から延びる。これは、例えば、まず、黄金角度だけk空間内において角度的にオフセットされる(3つの)放射状スポークに沿ってサンプリングすることにより、個々のk平面内において、中心領域が、完全にサンプリングされるか、又は、さらには、オーバーサンプリングされるという点で、実現される。k空間内における中心領域の周囲における1つ又はいくつかの中間の環状領域内で、再構成のためのMR信号の選択は、k軸から生じて中心のk空間領域をさらにサンプリングする放射状スポークに沿った選択に加えて、k軸からオフセットされた外向きの半径方向位置から延びた、いくつかの部分的な放射状スポークに沿ったものである。k空間の周辺において、1つ又は複数の中間の環状領域を半径方向に越えて、選択は、中間の環状領域を越えて外向きの半径方向位置から延びた、k空間内におけるさらなる部分的なスポークに沿って追加される。k空間内における部分的な放射状スポークの角密度は、半径方向外向きに増加する。
k空間の中心領域のサンプリングは、外向きの(k空間内の軸方向及び半径方向の)領域内でサンプルを追加する場合より高い時間レートでサンプルが追加されるように行われる。信号獲得において、スタックオブスターのそのスライスの中心から半径方向外向きに延びた半径方向プロファイルに沿ってk空間プロファイルを獲得することが好ましい。再構成のために、部分的なスポークからのMR信号のみが選択され得る。知られたPROPELLER像形成技術において獲得されたk空間プロファイルは、本発明の意味における放射状k空間プロファイルともみなされることに注意が必要である。PROPELLER概念において、MR信号は、デカルトベースのk空間サンプリングスキームにおけるL個の最低周波数位相符号化ラインに対応した、各々がL個の平行なk空間プロファイルからなる、N個のストリップ内のk空間内において獲得される。k空間ブレードとも呼ばれる各ストリップが、k空間内で、例えば180°/Nの角度だけ回転されることにより、MR信号セットの総セットがk空間内の円にわたって広がる。したがって、本発明の意味において、スタックオブスターは、「スタックオブブレード」と同じである。
本発明の好ましい実施形態によると、k空間プロファイルの放射状密度は、所与の視界(FOV:field−of−view)に従ったナイキスト基準がk空間の中心の周囲の楕円体内で満たされるように変化する。楕円体の寸法は、個々の像形成タスクの獲得速度と像品質との間で必要とされるトレードオフに従って選択され得る。本発明は、像エネルギーの支配的な中心のk空間内でオーバーサンプリング(すなわち、ナイキスト基準により要求されるものより高い半径方向サンプリング密度)を提供しながら、全走査期間を最小化する。半径方向サンプリング密度は、周辺のk空間領域内のより低い半径方向サンプリング密度まで徐々に低下する。周辺の領域内、すなわち、楕円体の外部において、放射状k空間サンプリング密度は、像品質に重大な影響を与えることなく、ナイキスト閾値未満にさえなる。k空間のオーバーサンプリングは、上述のように、ナイキスト楕円体内において本発明の方法に従って発生する。これは、像品質について妥協することなくMR像の各々を再構成するために、オーバーサンプリングされた中心のk空間領域内において互いに時間的に最も近くで獲得されたMR信号のみを選択することを可能にする。
本発明の好ましい実施形態によると、連続的に獲得されたk空間プロファイルの角度インクリメントは、黄金角度である。これは、上述の黄金角度スキームに対応し、黄金角度スキームにおいて、k空間プロファイルの角度は、各時点においてΔΦ=111.25°ずつインクリメントされる。これは、動的像形成の場合に特に有益である。
好ましくは、k空間プロファイルは、異なる角度位置でさらなるk空間プロファイルを獲得する前に、同じ角度位置で少なくとも2つの異なるスライスから獲得される。言い換えると、スライス方向に沿った位相符号化ステップは、異なる黄金角度位置でk空間プロファイルが獲得される前に、連続的に獲得される。これは、デカルトサンプリングの期間が短く維持されることを確実にし、このことが、スライスの積み重ね内における高いデータ一貫性をもたらし、全体的な動きロバスト性を維持する。
本発明に従って、k空間が可変密度でサンプリングされる。これは、(動き最適化された)可変密度の圧縮された検出により、走査期間を短くするために使用される。k空間プロファイルは、さらには、周辺のk空間領域内、すなわち、中でk空間が完全にサンプリングされる楕円体の外部において不規則に分布する。CSの理論は、大幅な信号データの低減のための大きな可能性を含むことが知られている。CS理論において、変換領域内のスパース表現を含む信号データセットは、適切な正規化アルゴリズムの適用により、アンダーサンプリングされた測定値から復元され得る。アンダーサンプリングの可能性は、大幅に短くなった獲得期間をもたらす。信号のサンプリング及び再構成ための数学的枠組みとして、CSが条件を規定し、この条件の下で、k空間サンプリング密度がナイキスト基準にはるかに満たない場合でも、信号データセットが厳密に、又は少なくとも高い像品質で再構成され得、CSはこのような再構成のための方法をさらに提供する。多くの既存のCSベースのMR獲得及び再構成スキームにおいて、MR信号データが変換領域内に散在するという条件のみを利用する基本CS式が使用される。未知のMR像に関する追加的な事前情報が、特定の用途において利用可能であるので、CS再構成内にこの事前情報を組み込むことが有益であることがさらに知られている。
本発明の技術は、完全にサンプリングされた小さな中心のk空間楕円体と、周辺のk空間領域内におけるアンダーサンプリングとを伴う、CSに十分に適した3D可変密度k空間サンプリングスキームを可能にする。
本発明の方法は、直接的な手法で従来のビューシェアリング、キーホール、又はk−tサンプリング技術と組み合わされる。半走査などの加速技術、又は、SENSE又は反復SENSEのような平行像形成技術を使用して、k空間プロファイルが獲得され、MR像へと再構成される。
例えば、時系列に沿ったMR像は、連続的に獲得されたMR信号から再構成され、周辺のk空間位置から獲得されたMR信号は、時間分解能を高めるために、系列の2つ以上の連続したMR像の間で共有(ビューシェアリング)される。
これまで説明された本発明の方法は、今のところ、検査ボリューム内に一様な静磁場を生成するための少なくとも1つの主電磁コイルと、検査ボリューム内において異なる空間的方向にスイッチングされる磁場勾配を生成するための複数の勾配コイルと、検査ボリューム内にRFパルスを生成するための、及び/又は、検査ボリューム内に位置する物体からMR信号を受信するための少なくとも1つのRFコイルと、RFパルスの時間的な連続性とスイッチングされる磁場勾配とを制御するための制御ユニットと、受信されたMR信号からMR像を再構成するための再構成ユニットとを含むMRデバイスにより実行され得る。本発明の方法は、例えば、MRデバイスの再構成ユニット及び/又は制御ユニットの対応するプログラムにより実施され得る。
本発明の方法は、現在の臨床用の多くのMRデバイスにおいて有益に実行され得る。この目的を達成するために、MRデバイスが本発明の上述の方法のステップを実施するようにMRデバイスが制御されるコンピュータプログラムを使用することが必要であるにすぎない。コンピュータプログラムは、MRデバイスの制御ユニットへのインストールのためにダウンロードするために、データ担体上に存在するか、又は、データネットワーク内に存在する。
添付の図面は、本発明の好ましい実施形態を開示する。しかし、図面は単に例示を目的として描かれることと、本発明の限定の定義を示すのではないこととが理解されなければならない。
本発明の方法を実施するためのMRデバイスを示す図である。 本発明のサンプリングスキームの一実施形態を概略的に示すk空間の図である。 ナイキスト限界に近い、より低い放射状密度で従来のように獲得されたMR像(図3a)と比較して、本発明に従った、中心のk空間スライス内においてより高い放射状密度で獲得及び再構成されたMR像(図3b)を示す図である。 動的像形成のための本発明の適用を概略的に示すk空間の図である。
図1を参照すると、MRデバイス1が示される。本デバイスは、実質的に一様で時間的に一定の主磁場が、検査ボリュームを通してz軸に沿って生成されるように、超伝導又は抵抗主電磁コイル2を備える。
磁気共鳴生成及び操作システムは、核磁気スピンを反転又は励起し、磁気共鳴を誘導し、磁気共鳴を再合焦し、磁気共鳴を操作し、磁気共鳴を空間的に及び別のやり方で符号化し、スピンを飽和させるなどのために、一連のRFパルスとスイッチングされる磁場勾配とを印加して、MR像形成を実施する。
より具体的には、勾配パルス増幅器3は、検査ボリュームのx軸、y軸、及びz軸に沿って全身勾配コイル4、5、及び6のうちの選択されたものに電流パルスを印加する。デジタルRF周波数送信器7は、送信/受信スイッチ8を介して全身ボリュームRFコイル9にRFパルス又はパルスパケットを送信して、検査ボリューム内にRFパルスを送信する。典型的なMR像形成シーケンスは、互いに一緒に取得された短い持続期間のRFパルスセグメントのパケットから構成され、任意の印加された磁場勾配が核磁気共鳴の選択された操作を達成する。RFパルスは、飽和させるために、共鳴を励起するために、磁化を反転させるために、共鳴を再合焦するために、又は共鳴を操作するために、及び検査ボリューム内に位置する物体10の一部を選択するために使用される。MR信号は、全身ボリュームRFコイル9によっても拾われる。
物体10の制限された領域のMR像の生成のために、ローカルアレイRFコイル11、12、13のセットが、像形成のために選択された領域に連続して位置する。アレイコイル11、12、13は、全身コイルRF送信により誘導されたMR信号を受信するために使用され得る。
結果として得られるMR信号は、全身ボリュームRFコイル9により、及び/又は、アレイRFコイル11、12、13により拾われ、好ましくは、前置増幅器(図示せず)を含む受信器14により復調される。受信器14は、送信/受信スイッチ8を介してRFコイル9、11、12、及び13に接続される。
ホストコンピュータ15は、勾配パルス増幅器3と送信器7とを制御して、例えば、エコープラナー像形成(EPI:echo planar imaging)、エコーボリューム像形成、勾配及びスピンエコー像形成、高速スピンエコー(TSE)像形成など、複数のMR像形成シーケンスのうちの任意のものを生成する。選択されたシーケンスに対して、受信器14は、各RF励起パルスに続いて迅速に連続して、単一又は複数のMRデータラインを受信する。データ獲得システム16は、受信された信号のアナログ・デジタル変換を実施して、各MRデータラインをさらなる処理に適したデジタル形式に変換する。現代のMRデバイスにおいて、データ獲得システム16は、生の像データの獲得に特化した独立したコンピュータである。
最終的に、デジタルの生の像データは、フーリエ変換又は他の適切な再構成アルゴリズムを適用する再構成プロセッサ17により像表現へと再構成される。MR像は、患者を通した平面スライス、平行な平面スライスのアレイ、三次元ボリュームなどを表す。次に、像が像メモリに記憶され、スライス、射影、又は像表現の他の部分を、例えば、結果として得られるMR像の、人間により可読な表示を提供するビデオ表示器18を介した可視化のための適切な形式に変換するために、像メモリにおいて像がアクセスされる。
続けて図1を参照し、図2〜図4をさらに参照すると、本発明の像形成アプローチの実施形態が説明される。
図2aに示されるように、3Dスタックオブスター獲得スキームに従って、(例えば、従来の勾配エコー像形成シーケンスを使用して)MR信号が獲得される。複数の空間的に非選択的又はスラブ選択的なRF励起の各々の後、1つ又は複数のMR信号が獲得され、各MR信号は、k空間プロファイルを表す。MR信号は、複数の(図2aに示す実施形態における7つの)平行スライス21、22、23、24、25、26、27から、放射状k空間プロファイルとして獲得される。スライスは、スライス方向kに沿って異なる位置に配置される。中心(k=k=0)の周りで回転する半径方向「スポーク」S1、S2、S3、S4、S5、S6、S7、S8に沿って各単一スライス内においてMR信号が獲得されながら、k方向においてデカルト位相符号化が実施される。これは、積み重ねられたディスクから構成された円柱k空間カバレッジをもたらす。図示された実施形態において、異なる角度位置でk空間プロファイルが獲得される前に、スライス方向kに沿った位相符号化ステップが連続的に獲得される。スポークS1、S2、S3、S4、S5、S6、S7、S8の角度的な順序付けに対して、上述の黄金角度スキームが使用される。スポークの角度は、各時点においてΔΦ=111.25°ずつインクリメントされる。本発明によると、k空間プロファイルの放射状密度、すなわち、スライス当たりの獲得されたスポークS1、S2、S3、S4、S5、S6、S7、S8の数は、スライス位置の関数として変化し、放射状密度は、より中心のk空間位置においてより高く、より周辺のk空間位置においてより低い。これは、例えば、次のように、図2aに示す実施形態において達成される。すなわち、第1のステップにおいて、スポークS1は、内側の3つのスライス23、24、及び25から獲得される。これは、スポークS2及びS3に対して連続的に行われる。高いデータ一貫性と全体的な動きロバスト性とを確実なものとするために必須の異なる黄金角度位置でk空間プロファイルをサンプリングする前に、スライス方向に沿った位相符号化ステップの連続した獲得が実施される。第2のステップにおいて、内側の5つのスライス22、23、24、25、及び26からのみ、スポークS4及びS5が獲得される。第3のステップにおいて、黒矢印により示されるように、すべての7つのスライス21、22、23、24、25、26、及び27から、スポークS6、S7、及びS8が獲得される。このように、k空間中心(k=0の周囲)は、k空間周辺より密にサンプリングされる。所与のFOVに従ったナイキスト基準が、k空間の中心の周囲(k=k=k=0)で楕円体29内において満たされる(図2b参照)ように、k空間プロファイル(スポーク)の放射状密度は、この手法により変えられる。図2aに示されるように、楕円体は楕円29に沿ってスライス21と交差する。楕円29内において、ナイキスト基準が満たされるとともに、k空間は、楕円29の外部においてアンダーサンプリングされる。スライス22、23、24内におけるそれぞれの楕円の半径は、これに対応してk=0に向けて大きくなる。任意選択的に、楕円体28が十分に大きなk空間の中心領域をカバーするので、そのk空間内において十分に大きな楕円体が楕円体「シャッター」として機能することが、獲得されたMR信号からMR像を再構成する前に、楕円体29の外部のk空間位置におけるMR信号を除去することにより使用される。
本発明は、ナイキスト基準に適合するk空間サンプリングを使用して全走査期間を最小化する。しかし、それはさらに、像エネルギーの支配的な中心のk空間内にオーバーサンプリングを提供し得る。半径方向サンプリング密度は、中心のスライス23、24、25から、周辺のスライス21、27内におけるより低い半径方向サンプリング密度まで徐々に低下する。楕円体29の外部において、放射状k空間密度は、像品質に重大な影響を与えることなく、ナイキスト閾値未満にさえなる。
k空間中心の周囲におけるより高い半径方向サンプリング密度に起因して、最小走査期間内で、本発明に従って筋状アーチファクトが効率良く低減される。
図3は、より高い放射状密度の効果を示す。より低い放射状密度で獲得された図3aに示すMR像と比較して、本発明に従ってより高い放射状密度でサンプリングする図3bに示すMR像内において、筋状アーチファクトが効率良く低減される。
k空間の中心領域の定期的な更新に起因して、本発明のk空間サンプリング技術は動的造影MR像形成の場合に興味深い。本発明に従って連続的に獲得されたMR信号から、時系列に沿ったMR像が再構成され、放射状k空間プロファイルは、より周辺のk空間位置におけるスライスからより、より中心のk空間位置におけるスライスから、より高い時間密度で獲得される。
これは、時間tの関数としてのk空間サンプリングの進行を示す図4aに示される。獲得は、k=0における中心のk空間スライス(左図)からの、3つのk空間スポークの獲得(黄金角度獲得)から始まる。ナイキスト基準は、描かれた塗りつぶした円内において満たされる。次のステップにおいて、さらに3つのスポークが3つの内側スライスの各々から獲得される。後続のステップにおいて、k空間プロファイルが、すべての5つのスライスから別の3つの半径方向に沿って獲得される。中でナイキスト基準が満たされる領域は、すべての3つの図内に塗りつぶした円として示される。このプロファイル順序スキームは、中でMR信号が近い時点において獲得される中心の3D k空間ナイキスト楕円体28をもたらす一方で、k空間の楕円体k空間シャッター29の外部の周辺は、図4bに示されるように、より遠い時点においてサンプリングされる。
k空間のオーバーサンプリングは、ナイキスト楕円体28内において発生する。図4cは、連続的に獲得されたスポークS1〜S12を含む、図4aに示すスライスの積み重ねのうちのk=0のスライスを示す。k空間の最も中心の円形領域内において、スポークS1、S2、S3のMR信号のみが、像再構成のために選択される。このMR信号データは、k空間の中間環状領域内におけるスポークS1、S2、S3、S4、S5、S6からのMR信号により補完される。k空間の残りの周辺の領域内におけるすべてのスポークS1〜S12からのMR信号によりMR信号データが完全なものとなる。このように、より遠い時点においてサンプリングされた、周辺のk空間領域から獲得されたk空間プロファイルと組み合わせて、それぞれ最も近い時点において獲得されたナイキスト楕円体内のk空間プロファイルのみを使用して、時系列の各MR像が再構成される。図4aに示される例において、時間が進行するにつれて、まず、より中心のk空間位置がサンプリングされ、次に、より周辺のk空間位置がサンプリングされる。この順序は、反転されてよい。代替的に、サンプリング位置のk値は、時間が進行するにつれて、例えば、線形に増加又は減少する。
この技法は、造影MR像形成の改善を可能にする。目的の領域内における造影剤の取り込みは、従来技術と比べて、より高い精度で追跡され得る。時間的に近くで獲得されたMR信号のみが、k空間中心内におけるMR像再構成のために使用されるので、造影剤が特定の像ボクセルに達する時点が、より明確に識別され得る。したがって、測定された造影プロファイルの不鮮明さが低減される。

Claims (13)

  1. MRデバイスの検査ボリューム内に位置する物体のMR像形成の方法であって、前記方法が、
    少なくとも1つのRFパルスとスイッチングされる磁場勾配との像形成シーケンスを前記物体に施すステップと、
    スタックオブスター(stack-of-stars)スキームに従って、MR信号を獲得するステップであって、前記MR信号が、複数の異なる角度位置で半径方向にk空間プロファイルとして、及び、スライス方向に沿った異なる位置に配置された複数の平行スライスから獲得され、前記k空間プロファイルの放射状密度が、スライス位置の関数として変化し、前記放射状密度が、より中心のk空間位置においてより高く、より周辺のk空間領域内の位置においてより低く、前記k空間プロファイルが、より周辺のk空間位置におけるスライスからより、より中心のk空間位置におけるスライスから、より高い時間密度で獲得されるステップと、
    前記MR信号からMR像を再構成するステップと、
    を有する、方法。
  2. 前記MR像を再構成するために、k空間内における放射状スポークに沿って、互いに時間的に最も近くで獲得されたMR信号が、k空間の中心の周りの楕円体内において選択され、前記k空間プロファイルの前記放射状密度は、所与の視界に従ったナイキスト基準が前記k空間の中心の周りの前記楕円体内において満たされるように、放射状密度変化する、
    請求項1に記載の方法。
  3. 前記MR像を再構成するために、中心の楕円体を越えた周辺のk空間領域内にあるk空間内の部分的な放射状スポークからのMR信号のみが選択され、前記k空間内のる前記部分的な放射状スポークは、外向きの放射状k空間位置から延び、前記部分的な放射状スポークは外向きに周辺のk空間MR楕円体となるにつれて、前記部分的な放射状スポークがより高い角密度をもつ、
    請求項1又は2に記載の方法。
  4. 前記楕円体の外部のk空間位置におけるMR信号が、前記MR像を再構成する前に除去される、
    請求項1乃至3のいずれか一項に記載の方法。
  5. 前記k空間プロファイルが、異なる角度位置でさらなるk空間プロファイルを獲得する前に、同じ角度位置で少なくとも2つの異なるスライスから獲得される、
    請求項1乃至4のいずれか一項に記載の方法。
  6. 前記楕円体の外部よりも前記楕円体内において、所与の視界に従ったナイキスト基準より高い放射状密度が使用され、前記楕円体の外部における前記k空間プロファイルの前記放射状密度が、前記ナイキスト基準以下であるように、前記k空間プロファイルの前記放射状密度が変化する、
    請求項1乃至5のいずれか一項に記載の方法。
  7. 連続的に獲得されたk空間プロファイルの角度インクリメントが、黄金角度である、
    請求項1乃至6のいずれか一項に記載の方法。
  8. 前記k空間プロファイルが、k空間の周辺の領域内において不規則に分布する、
    請求項1乃至7のいずれか一項に記載の方法。
  9. 前記MR像が、圧縮された検出を使用して再構成される、
    請求項1から請求項8のいずれか一項に記載の方法。
  10. 時系列のMR像が、連続的に獲得されたMR信号から再構成され、前記楕円体を越えた周辺のk空間位置から獲得されたMR信号が、前記時系列の2つ以上の連続したMR像の間で共有される、
    請求項1乃至9のいずれか一項に記載の方法。
  11. 前記MR信号が、複数のRF受信アンテナを介して並列に獲得され、前記MR像が、SENSEを使用して再構成される、
    請求項1乃至10のいずれか一項に記載の方法。
  12. 検査ボリューム内に一様な静磁場を生成するための少なくとも1つの主電磁コイルと、前記検査ボリューム内において異なる空間的方向にスイッチングされる磁場勾配を生成するための複数の勾配コイルと、前記検査ボリューム内にRFパルスを生成するための、及び/又は、前記検査ボリューム内に位置する物体からMR信号を受信するための少なくとも1つのRFコイルと、前記RFパルスの時間的な連続性とスイッチングされる磁場勾配とを制御するための制御ユニットと、受信された前記MR信号からMR像を再構成するための再構成ユニットと、を含むMRデバイスであって、前記MRデバイスが、
    少なくとも1つのRFパルスとスイッチングされる磁場勾配との像形成シーケンスを前記物体に施すことと、
    スタックオブスター(stack-of-stars)スキームに従ってMR信号を獲得することであって、前記MR信号が、スライス方向に沿った異なる位置に配置された複数の平行スライスから、半径方向のk空間プロファイルとして獲得され、前記k空間プロファイルの放射状密度が、スライス位置の関数として変化し、前記放射状密度が、より中心のk空間位置においてより高く、より周辺のk空間領域内においてより低く、前記k空間プロファイルが、より周辺のk空間位置におけるスライスからより、より中心のk空間位置におけるスライスから、より高い時間密度で獲得される、前記MR信号を獲得することと、
    前記MR信号からMR像を再構成することと、
    を実施する、MRデバイス。
  13. MRデバイス上で動作するコンピュータプログラムであって、前記コンピュータプログラムが、
    少なくとも1つのRFパルスとスイッチングされる磁場勾配との像形成シーケンスを生成することと、
    スタックオブスター(stack-of-stars)スキームに従ってMR信号を獲得することであって、前記MR信号が、スライス方向に沿った異なる位置に配置された複数の平行スライスから、半径方向のk空間プロファイルとして獲得され、前記k空間プロファイルの前記放射状密度が、スライス位置の関数として変化し、前記放射状密度が、より中心のk空間位置においてより高く、より周辺のk空間領域内においてより低く、
    前記k空間プロファイルが、より周辺のk空間位置におけるスライスからより、より中心のk空間位置におけるスライスから、より高い時間密度で獲得される、前記MR信号を獲得することと、
    前記MR信号からMR像を再構成することと、
    を行うための命令を含む、コンピュータプログラム。
JP2021001435A 2015-06-15 2021-01-07 スタックオブスター獲得を使用したmr像形成 Pending JP2021090760A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP15172087.7 2015-06-15
EP15172087 2015-06-15

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2017564875A Division JP6865697B2 (ja) 2015-06-15 2016-06-12 スタックオブスター獲得を使用したmr像形成

Publications (1)

Publication Number Publication Date
JP2021090760A true JP2021090760A (ja) 2021-06-17

Family

ID=53397943

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2017564875A Active JP6865697B2 (ja) 2015-06-15 2016-06-12 スタックオブスター獲得を使用したmr像形成
JP2021001435A Pending JP2021090760A (ja) 2015-06-15 2021-01-07 スタックオブスター獲得を使用したmr像形成

Family Applications Before (1)

Application Number Title Priority Date Filing Date
JP2017564875A Active JP6865697B2 (ja) 2015-06-15 2016-06-12 スタックオブスター獲得を使用したmr像形成

Country Status (5)

Country Link
US (1) US10670678B2 (ja)
EP (1) EP3308186B1 (ja)
JP (2) JP6865697B2 (ja)
CN (1) CN107750338B (ja)
WO (1) WO2016202707A1 (ja)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10353034B2 (en) * 2016-05-05 2019-07-16 Siemens Healthcare Gmbh Multi-echo pseudo-golden angle stack of stars thermometry with high spatial and temporal resolution using k-space weighted image contrast
EP3413071A1 (en) 2017-06-09 2018-12-12 Koninklijke Philips N.V. Mr imaging using a stack-of-stars acquisition with variable contrast
US20190054194A1 (en) * 2017-08-17 2019-02-21 Case Western Reserve University System and Method For Evaluation of Subjects Using Magnetic Resonance Imaging and Oxygen-17
EP3454070A1 (en) * 2017-09-06 2019-03-13 Commissariat à l'Energie Atomique et aux Energies Alternatives Method and apparatus for accelerated magnetic resonance imaging
EP3467533A1 (en) * 2017-10-09 2019-04-10 ETH Zurich Magnetic resonance imaging method with hybrid filling of k-space
EP3470869A1 (en) 2017-10-10 2019-04-17 Koninklijke Philips N.V. Mr imaging using motion-dependent radial or spiral k-space sampling
EP3517988A1 (en) 2018-01-29 2019-07-31 Koninklijke Philips N.V. Mr imaging using a stack-of-stars acquisition with intrinsic motion correction
EP3579009A1 (en) * 2018-06-05 2019-12-11 Koninklijke Philips N.V. Zero echo time mr imaging with water-fat separation
EP3598162B1 (en) * 2018-07-19 2024-10-02 Siemens Healthineers AG Magnetic resonance method and apparatus for suppressing metal artifacts
US10746830B2 (en) * 2018-08-28 2020-08-18 General Electric Company Systems and methods for hybrid slice encoding in three-dimensional magnetic resonance imaging
CN109696647B (zh) * 2019-02-21 2021-05-28 奥泰医疗系统有限责任公司 三维多次激发扩散加权成像的k空间采集方法及重建方法
EP3709040A1 (de) * 2019-03-13 2020-09-16 Siemens Healthcare GmbH Passive magnetfeldkamera und verfahren zum betrieb der passiven magnetfeldkamera
EP3709042A1 (en) * 2019-03-14 2020-09-16 Koninklijke Philips N.V. Mr imaging using a 3d radial or spiral acquisition with soft motion gating
JP7292930B2 (ja) * 2019-04-03 2023-06-19 キヤノンメディカルシステムズ株式会社 画像再構成方法及び再構成装置
EP3748385A1 (en) 2019-06-04 2020-12-09 Koninklijke Philips N.V. Optimized k-space profile ordering for 3d radial or spiral mr imaging
EP3842819A1 (en) 2019-12-23 2021-06-30 Koninklijke Philips N.V. 3d mr imaging with intrinsic motion detection
CN111257809B (zh) * 2020-01-23 2022-05-17 上海东软医疗科技有限公司 磁共振成像方法、装置、存储介质及医疗设备
JP7509546B2 (ja) * 2020-02-12 2024-07-02 キヤノンメディカルシステムズ株式会社 画像再構成装置
US11802926B2 (en) * 2020-11-12 2023-10-31 Oxford University Innovation Limited Methods and systems for optimizing sampling efficiency of dynamically sampled magnetic resonance imaging (MRI)
CN112557978B (zh) * 2020-11-26 2022-03-25 上海交通大学 一种基于三维非笛卡尔轨迹的多回波采样方法
GB2608404B (en) * 2021-06-30 2023-07-05 Siemens Healthcare Gmbh Method of reconstructing a dynamic series of motion-compensated magnetic resonance images
US11940516B2 (en) 2022-01-29 2024-03-26 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for magnetic resonance imaging
CN115430828B (zh) * 2022-09-22 2024-07-23 济南海圣机电科技有限公司 一种浇注机铁水定量定速浇注控制方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11276456A (ja) * 1998-03-30 1999-10-12 Shimadzu Corp Mrイメージング装置
JP2008529643A (ja) * 2005-02-11 2008-08-07 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 短いエコー時間での磁気共鳴イメージング
JP2011507603A (ja) * 2007-12-20 2011-03-10 ウイスコンシン アラムナイ リサーチ ファウンデーシヨン 先験的画像制限画像再構成法
JP2012505709A (ja) * 2008-10-14 2012-03-08 メイヨ フォンデーシヨン フォー メディカル エジュケーション アンド リサーチ 移動テーブルmri装置及び方法
JP2012196537A (ja) * 2006-09-06 2012-10-18 Toshiba Corp 磁気共鳴イメージング装置
JP2014508622A (ja) * 2011-03-24 2014-04-10 コーニンクレッカ フィリップス エヌ ヴェ 事前情報に制約される正則化を使用するmr画像再構成
US20150077112A1 (en) * 2012-04-19 2015-03-19 New York University System, method and computer-accessible medium for highly-accelerated dynamic magnetic resonance imaging using golden-angle radial sampling and compressed sensing
JP2015510804A (ja) * 2012-03-19 2015-04-13 コーニンクレッカ フィリップス エヌ ヴェ K空間領域の中央および周辺のサンプリング最中での呼吸運動の検出を用いた磁気共鳴画像の再構成

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6404194B1 (en) 1994-08-05 2002-06-11 The Board Of Trustees Of The Leland Stanford Junior University Three dimensional magnetic resonance imaging using spiral trajectories in k-space
US5713358A (en) 1996-03-26 1998-02-03 Wisconsin Alumni Research Foundation Method for producing a time-resolved series of 3D magnetic resonance angiograms during the first passage of contrast agent
US5881728A (en) 1996-07-26 1999-03-16 Wisconsin Alumni Research Foundation Digital subtraction magnetic resonance angiography with image artifact suppression
US5873825A (en) 1997-04-11 1999-02-23 Wisconsin Alumni Research Foundation Three dimensional digital subtraction magnetic resonance angiography with limited k-space mask
DE19816961A1 (de) * 1998-04-17 1999-06-17 Siemens Ag Verfahren zur Gewinnung eines dreidimensionalen Rohdatensatzes für die MR-Bildgebung und Vorrichtung zur Durchführung des Verfahrens
DE102006017470B4 (de) * 2006-04-13 2010-07-22 Siemens Ag Verfahren zur Erzeugung von Magnetresonanzbildern und Magnetresonanzanlage hierfür
US8120358B2 (en) 2006-04-13 2012-02-21 The Regents Of The University Of California Magnetic resonance imaging with high spatial and temporal resolution
WO2007124151A2 (en) 2006-04-21 2007-11-01 The Trustees Of The University Of Pennsylvania Rapid 3-dimensional bilateral breast mr imaging
EP2461175A1 (en) * 2010-12-02 2012-06-06 Koninklijke Philips Electronics N.V. MR imaging using a multi-point Dixon technique
WO2012102338A1 (ja) 2011-01-28 2012-08-02 株式会社 日立メディコ 磁気共鳴イメージング装置及びそれを用いた血管撮像方法
DE102011078273B4 (de) * 2011-06-29 2014-11-20 Siemens Aktiengesellschaft MR-Angiographie mit nicht-kartesischer Signalakquisition
JP2013192578A (ja) * 2012-03-16 2013-09-30 Ge Medical Systems Global Technology Co Llc 磁気共鳴イメージング装置及びプログラム
US8744551B2 (en) 2012-07-12 2014-06-03 Northshore University Healthsystem Method for non-contrast enhanced magnetic resonance angiography
WO2014085288A1 (en) * 2012-11-29 2014-06-05 The Regents Of The University Of California Noninvasive 4-d time-resolved dynamic magnetic resonance angiography
DE102013205208B4 (de) * 2013-03-25 2015-02-12 Siemens Aktiengesellschaft Multiecho-Magnetresonanz-Messsequenz mit erhöhter Auflösung
US10052033B2 (en) * 2015-03-12 2018-08-21 Siemens Healthcare Gmbh System and method for motion-robust 3D magnetic resonance imaging of vessel walls

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11276456A (ja) * 1998-03-30 1999-10-12 Shimadzu Corp Mrイメージング装置
JP2008529643A (ja) * 2005-02-11 2008-08-07 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 短いエコー時間での磁気共鳴イメージング
JP2012196537A (ja) * 2006-09-06 2012-10-18 Toshiba Corp 磁気共鳴イメージング装置
JP2011507603A (ja) * 2007-12-20 2011-03-10 ウイスコンシン アラムナイ リサーチ ファウンデーシヨン 先験的画像制限画像再構成法
JP2012505709A (ja) * 2008-10-14 2012-03-08 メイヨ フォンデーシヨン フォー メディカル エジュケーション アンド リサーチ 移動テーブルmri装置及び方法
JP2014508622A (ja) * 2011-03-24 2014-04-10 コーニンクレッカ フィリップス エヌ ヴェ 事前情報に制約される正則化を使用するmr画像再構成
JP2015510804A (ja) * 2012-03-19 2015-04-13 コーニンクレッカ フィリップス エヌ ヴェ K空間領域の中央および周辺のサンプリング最中での呼吸運動の検出を用いた磁気共鳴画像の再構成
US20150077112A1 (en) * 2012-04-19 2015-03-19 New York University System, method and computer-accessible medium for highly-accelerated dynamic magnetic resonance imaging using golden-angle radial sampling and compressed sensing

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
LI, ZHITAO: ""Highly Accelerated 3D Dynamic Imaging with Variable Density Golden Angle Stack-of-Stars Sampling"", PROC.INTL.SOC.MAG.RESON.MED., vol. 21, JPN6020013817, 26 April 2013 (2013-04-26), US, pages 3797, ISSN: 0004930508 *

Also Published As

Publication number Publication date
JP6865697B2 (ja) 2021-04-28
WO2016202707A1 (en) 2016-12-22
US20180149721A1 (en) 2018-05-31
CN107750338A (zh) 2018-03-02
EP3308186B1 (en) 2022-02-02
JP2018519050A (ja) 2018-07-19
CN107750338B (zh) 2020-12-22
EP3308186A1 (en) 2018-04-18
US10670678B2 (en) 2020-06-02

Similar Documents

Publication Publication Date Title
JP6865697B2 (ja) スタックオブスター獲得を使用したmr像形成
JP6998218B2 (ja) 動き検出を用いるmr撮像
US9964615B2 (en) MR image reconstruction using compressed sensing
JP7075420B2 (ja) 可変コントラストのスタック・オブ・スター収集を使用したmrイメージング
JP6445470B2 (ja) パラレルのマルチスライスmr画像法
JP6074126B1 (ja) k空間中心のサンプリングを用いるゼロエコー時間MR撮像
JP7446328B2 (ja) ソフト動きゲーティングによる3dラジアル又はスパイラル収集を用いたmr画像
US20170307716A1 (en) Propeller mr imaging with artefact suppression
CN107810425B (zh) 消除非t2加权信号贡献的t2加权mr成像
JP7023954B2 (ja) プロペラmrイメージング
JP7060685B2 (ja) 体動依存型ラジアル又はスパイラルk空間サンプリングを使用するMRイメージング
JP2020522344A (ja) 並列マルチスライスmr撮像
JP2021528143A (ja) スパイラル取得のためのmrファントム
EP3118643A1 (en) Dynamic propeller mr imaging
JP2022539993A (ja) 強調3d放射状mr撮像

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210204

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210204

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220125

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220304

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20220601

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20221128