JP2016512780A - 圧縮センシングを使用したmr画像の再構成 - Google Patents

圧縮センシングを使用したmr画像の再構成 Download PDF

Info

Publication number
JP2016512780A
JP2016512780A JP2016503743A JP2016503743A JP2016512780A JP 2016512780 A JP2016512780 A JP 2016512780A JP 2016503743 A JP2016503743 A JP 2016503743A JP 2016503743 A JP2016503743 A JP 2016503743A JP 2016512780 A JP2016512780 A JP 2016512780A
Authority
JP
Japan
Prior art keywords
space
image data
signal data
data
signal
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.)
Withdrawn
Application number
JP2016503743A
Other languages
English (en)
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.)
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 JP2016512780A publication Critical patent/JP2016512780A/ja
Withdrawn 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
    • 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
    • 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/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56545Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by finite or discrete sampling, e.g. Gibbs ringing, truncation artefacts, phase aliasing artefacts
    • 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/4828Resolving the MR signals of different chemical species, e.g. water-fat imaging

Landscapes

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

Abstract

本発明は、MR機器の検査ボリュームの中に配置された対象物のMR画像化方法に関する。本発明は、MR信号データを取得するために画像化シーケンスに対象物をかけるステップであり、MR信号データはk空間のサブサンプリングを用いた不規則k空間サンプリングパターンを使用することによりk空間位置と時間の関数として取得されるステップと、MR信号データからMR画像データを再構成するステップであり、MR画像データは空間次元と周波数次元を含み、MR画像データにおけるサブサンプリングアーチファクトを抑圧するために変換領域におけるMR画像データの希薄性が利用されるステップと、を含む。さらに、本発明は、MR機器およびコンピュータプログラムに関する。

Description

本発明は、磁気共鳴(MR)画像化の分野に関する。より特定的には、MR機器の検査ボリュームの中に配置された患者の体の少なくとも一部分のMR画像化方法に関する。本発明は、また、MR機器およびMR機器上で実行されるコンピュータプログラムにも関する。
2次元または3次元画像を形成するために、磁場(magnetic field)と核スピン(nuclear spin)との間の相互作用を利用する画像形成MR方法が、特に医療診断の分野において、現在幅広く使用されている。軟部組織(soft tissue)の画像化のために、多くの観点で他の画像化方法よりも優れているからである。電離放射線(ionizing radiation)を必要としないし、たいていは侵襲的でない。
一般的にMR方法に従えば、検査されるべき患者の体は、強力で均一な磁場の中に配置される。磁場の方向は、測定が基礎とする座標系の軸(通常はz軸)を、同時に定義する。磁場は、磁場強度に従って、個々の核スピンに対して異なるエネルギーレベルを生成する。定められた周波数(いわゆるラーモア(Larmor)周波数、またはMR周波数)の電磁気交番(alternating)磁界(RF field)の適用によって励起され得るもの(スピン共鳴)である。マクロ的な観点から、個々の核スピンの分布は、全体の磁化を生成し、適切な周波数の電磁気パルス(RFパルス)の適用によって平衡状態の外に偏向され得る。一方で、RFパルスの磁場は、z軸に対して垂直に延びている。磁化が、z軸に関して歳差運動(precession)を行うようにである。磁化のこの動作は、開口の角度がフリップ角(flip angle)として参照される円錐形の表面を記述する。フリップ角の大きさは、適用される電磁気パルスの強度と持続時間に依存している。いわゆる90°パルスの場合、スピンは、z軸から横断面(フリップ角90°)へ偏向される。RFパルスは、MR機器のRFコイル構成を介して患者の体に向けて放射される。RFコイル構成は、典型的に、患者の体が配置される検査ボリュームを取り囲んでいる。
RFパルスの終了後に、磁化は緩和して最初の平衡状態へ戻る。そこでは、第1時定数T(スピン格子または縦緩和時間)を伴うz方向の磁化が再び確立される。そして、z方向に垂直な磁化が第2時定数T(スピン−スピンまたは横緩和時間)を伴って緩和する。磁化の変化は、MR機器の検査ボリュームの中に配置され方向付された受信RFアンテナまたはコイルによって検出され得る。磁化の変化がz軸に対して垂直な方向で測定されるような方法においてである。例えば、90°パルスの適用の後には、横磁化の崩壊(decay)を同時に伴う。核スピン(ローカル磁場の不均質性によって引き起こされるもの)の、同一の位相を伴う秩序状態から、全ての位相拡張領域が不均一に分布されている状態への移行(ディフェイジング(dephasing))によるものである。ディフェイジングは、パルスの再フォーカスによって(例えば180°パルス)補償され得る。このことは、受信コイルにおいてエコー信号(スピンエコー)を生成する。
体における空間的解像度を実現するために、均一な磁場の上に、3つの主軸に沿って延びる線形磁場勾配が重ねられて(superposed)、スピン共鳴周波数の線形空間依存性を導いている。受信コイルにおいてピックアップされた信号は、すると、体の中の異なる場所に関連し得る異なる周波数のコンポーネントを含んでいる。受信RFアンテナまたはコイルを介して取得された信号データは、空間周波数領域に対応しており、k空間データと呼ばれる。k空間データは、たいてい、異なる位相エンコードを用いて取得された複数のラインを含んでいる。それぞれのラインは、数多くのサンプルを収集することによってデジタル化される。一式のk空間データは、フーリエ変換、または、他の本来知られた再構成技術によって、MR画像へ転換される。
多くのMR画像化アプリケーションにおいては、画像化スピードが重要である。しかしながら、MR画像の再構成のために要求されるMR信号が収集され得るところのスピードは、基本的に、物理的および生理学制約によって制限されている。従って、MR画像化領域における多くの最近の開発は、再構成されたMR画像の品質を劣化することなく、取得される信号データの量を削減することを目指している。多くのそうした開発の中で、圧縮センシング(Compressed Sensing(CS))の理論は、顕著な信号データ削減のための多大なポテンシャルを有している。CS理論においては、変換領域において希薄な表現を有する信号データセットが、適切な正規化(regularisation)アルゴリズムの適用によって、サブサンプルされた測定から回復され得る。サブサンプリングの可能性は、著しく削減された取得時間を導く。信号サンプリングと再構成のための数学的フレームワークとして、k空間サンプリング密度がナイキスト基準(Nyquist criterion)よりはるかに低い場合でさえ、信号データセットが正確に又は少なくとも高画像品質を伴って再構成され得る条件を規定する。そして、また、そうした再構成の方法も提供する。大多数の既存のCS基盤のMR取得および再構成スキームにおいては、基本的CS定式が使用される。MR信号データは変換領域において希薄であるという必須条件だけを利用するものである。例えば、M.Lustigらは、迅速なMR画像化のためのCS適用を提案してきた(M.Lustig et al.:”Sparse MRI:The Application of Compressed Sensing for Rapid MR Imaging”、Magnetic Resonances in Medicine、58、1182−1195、2007)。所定のアプリケーションにおいては知られていないMR画像に関する追加の従来情報が利用可能であり得るので、この従来情報をCS再構成の中に組み入れることが有利であることも、また、知られている。
MR取得を加速するためのさらなる技術は、並行獲得(parallel acquisition)技術として知られている。このカテゴリーにおける方法は、SENSE(Pruessmann et al.,”SENSE:Sensitivity Encoding for Fast MRI”、Magnetic Resonances in Medicine 1999、42(5)、1952−1962)およびSMASH(Sodickson et al.,”Simultaneous acqusition of spatial harmonics(SMASH):Fast imaging with radiofrequency coil arrays”、Magnetic Resonances in Medicine 1997、38、591−603)である。SENSEおよびSMASHは、複数のRF受信アンテナ(コイル)から並行して獲得されたサブサンプルk空間データ取得を使用する。これらの方法においては、複数のコイルからの(複合)信号データが、複合重み付けを用いて結合される。最終的な再構成されたMR画像において、サブサンプリングアーチファクトを抑圧するようにである。このタイプの複合アレイ結合は、ときどき、空間フィルタリングとして参照され、かつ、k空間領域(SMASHにおけるように)において、または、画像領域(SENSEにおけるように)において実行される結合を含んでいる。ハイブリッドである方法も同様である。SENSEまたはSMASHのいずれでも、高い正確性を伴って、適切な重み付け又はコイルの感度を知ることが非常に重要である。コイルの感度、つまり、信号検出のために使用されるアレイRFコイルの空間感度(spatial sensitivity)プロフィール、を取得するために、実際の画像を取得する以前及び/又は後で、較正プリスキャン(pre−scan)が典型的には実行される。プリスキャン、ときどき参照スキャンとしても参照されるもの、において、MR信号は、たいてい、最終の診断MR画像のために要求される解像度よりも著しく低い解像度において取得される。
MR画像の取得および再構成のための技術が、最近知られるようになった。そこでは、圧縮センシング再構成が並行MR信号取得と組み合わされている。
本発明の目的は、希薄にサンプルされたMR信号データから高品質MR画像を再構成できるようにする改善した方法を提供することである。
本発明に従って、MR機器の検査ボリュームの中に置かれた患者の少なくとも一部分に係るMR画像化の方法が開示される。本方法は、
−MR信号データを取得するために画像化シーケンスに前記対象物をかけるステップであり、前記MR信号データは、k空間のサブサンプリングを用いた不規則k空間サンプリングパターンを使用することにより、k空間位置と時間の関数として取得されるステップと、
−前記MR信号データからMR画像データを再構成するステップであり、前記MR画像データは空間次元と周波数次元を含み、前記MR画像データにおけるサブサンプリングアーチファクトを抑圧するために変換領域における前記MR画像データの希薄性が利用されるステップと、を含んでいる。
本発明の非常に重要な態様は、サブサンプルされたMR信号の時間次元を用いた取得である。時間次元は、再構成されたMR画像データの周波数次元に対応するものである。本発明の方法は、従って、空間−スペクトラムMR画像化(spatial−spectral MR imaging)の方法であると言われ得る。この目的のために、磁気共鳴の励起と実際のMR信号取得の瞬間との間の異なる時間間隔に対して本発明に従ってMR信号が取得される。MR信号は、例えば、スピンエコー又は勾配エコーとして生成されてよい。そこで、MR信号は、本発明に従って、異なるエコー時間値に対して取得される。異なるエコー時間値は、この場合におけるMR信号データの時間次元に及ぶ。本発明に係る方法によって、空間座標に従属され得るスペクトラム情報が提供される。異なるタイプの画像コントラストを提供するために、スペクトラム情報が使用され得る。スペクトラム情報は、さらに、脂肪プロトン(fat proton)から発しているMR信号を、水プロトン(water proton)から発しているMR信号と分離するために使用され得る。効果的な水/脂肪の分離が、このようにして達成され得る。
規則的サンプリングと反対に、本発明に従った不規則k空間サンプリングは、MR画像の再構成をより難しくしている。それ自身で、この複雑性は不利である。しかしながら、k空間の規則的サブサンプルは、再構成されたMR画像の中で少数の分離した位置においてサブサンプリングアーチファクトを生じさせることを考慮に入れるべきである。それに反して、本発明に従って不規則k空間サンプリングを用いると、不完全なMR画像再構成が同様な(または、さらにより表された)アーチファクトを結果として生じる。しかし、これらのアーチファクトは、再構成されたMR画像において「スクランブル(”scramble”)」される。別の言葉で言えば、不規則k空間サンプリングは、インコヒーレント(incoherent)なエイリアシング(aliasing)アーチファクトを結果として生じる。そうしたインコヒーレントなアーチファクトは、k空間の規則的サンプリングによって取得されたコヒーレントなアーチファクトよりも容易に抑圧され得ることが、本発明に従って、利用され得る。用語「エイリアシングアーチファクト」は、このコンテクストにおいては、検査下にある患者の体部分のサイズに関するナイキスト基準によって要求される密度と比較して、k空間のサブサンプリングから現れる全てのアーチファクトを参照する。これらのインコヒーレントなアーチファクトは、本発明に従って、変換領域におけるMR画像データの希薄性を利用することによって取り除かれる。例えば、本来知られた圧縮センシング技術によるものである。変換領域においては、画像データが希薄であるが、アーチファクトは有用な情報から分離されている。例えば、非線形フィルタリングオペレーションを使用することによるものである。
変換領域、そこではMR画像データが希薄であるが、への変換は、「希薄化(”sparsifying”)」変換として参照され得る。本発明の望ましい実施例に従えば、周波数次元におけるMR画像データに係るこの希薄化変換の特性は、空間次元における特性とは異なるものである。このことは、次の原則に基づいてよい。つまり、MR画像データのエッジ(つまり、隣接する画像特徴間の境界)は、空間次元において希薄であると想定され、一方、MR画像データ自身は、周波数次元において希薄であると想定される、ことである。周波数領域においては、各空間位置について、非常に限定された一式の周波数値だけに対して、MR画像データが非ゼロであると期待される。MR画像データが非ゼロであると期待される周波数位置は、典型的に、時たまの突然な変化を除いて、空間座標における連続関数である。そうした変化は、おそらく、既定の値と一致する。例えば、3.4ppm(つまり、脂肪プロトンと水プロトンとの間の周波数における差異)である。さらに、これらの変化は、おそらく、空間次元におけるエッジと一致する。他の基準が想像可能である。
k空間次元および(エコー)時間次元の両方におけるMR信号取得のための不規則サンプリングパターンの利点は、特には、切替え磁場勾配を正確にコントロールする必要が無い、短いスキャン時間である。一方で、正確なk空間座標および時間値が、もちろん、知られる必要がある。それぞれのMR信号サンプルに対して、再構成ができるようにするためである。
本発明の一つの実施例においては、連続したランダムk空間軌道を含むように、k空間サンプリングスキームが選択される。例えば、MR信号取得の最中に、ランダムまたは擬似ランダムに磁場勾配を変化させることが可能である。時間値と同様にk空間座標値も、再構成できるようにするために、取得プロシージャの最中にそれぞれのMR信号サンプルに対して割り当てられる。本発明のこの実施例では、使用されるMR装置の勾配システムに対する技術要求は低い。特別に速い勾配切替えは必要とされない。勾配切替えによって生じるノイズレベルは、検査される患者について許容できる範囲内にうまく維持され得る。
本発明の方法は、望ましくは、非デカルトサンプリングスキームと組み合わされる。スパイラルスキャニングが適用されてよい。ラジアルk空間サンプリングも、また、可能な変形である。
本発明に従った、k空間スキャニングは不規則であるが、k空間のスキャニングの密度への程度のコントロールが、k空間サンプリングのタイミングスキームと同様に、本発明に係る再構成問題を最適に条件付けるために確立されるべきである。再構成問題の条件付けに関しては、MR信号データの取得のためのk空間領域を再訪問(re−visit)する連続的なk空間軌道を含むように、k空間サンプリングパターンを選択することが有利である、ということがわかる。再訪問とは、異なる(エコー)時間値においてk空間の同一領域からMR信号が取得されることを意味する。k空間領域を再訪問することによって、k空間次元および(エコー)時間次元によって拡張された空間が、MR画像の再構成に関して適切にカバーされることが保証される。別の言葉で言えば、不規則(ランダム)k空間サンプリングスキームが、本発明に従って、制限された空間−時間的(spatio−temporal)サンプリング密度を用いて有利に適用され得る
本発明の一つの実施例においては、スキャンされた領域にわたるほとんど均一な密度が目的とされ得る。別の実施例においては、k空間サンプリングパターンは、k空間の中心領域が周辺領域よりも高い密度でサンプルされるところで好まれ得る。
本発明に係る方法は、例えば、以下の記事に開示されたk空間サンプリングと有利に組み合わされてよい。”Fast Three−Dimensional k−Space Trajectory Design Using Missile Guidance Ideas”、Magnetic Resonances in Medicine 52:329−336(2004)である。一式のそうした軌道、k空間中心から異なる方向でそれぞれに動き始めるもの、が使用され得る。このようにして、k空間のそれぞれの領域は、それぞれの軌道に対する異なるエコー時間を用いて、軌道ごとに一度訪問される。
k空間において外側に、次に内側に再び進行しているスパイラルまたはスパイラルのような軌道も、また、本発明に従って、k空間サンプリングに対して適用され得る。k空間における各スパイラルの延長は、ある程度ランダムに行われてよい。k空間の外側領域が、内側領域よりも頻度が低く訪問されるようにである。望ましくは、k空間サンプリングの軌道は、正確なスパイラルでなく、ランダムな摂動を伴うスパイラルである。サンプリング密度の制限が定められてよく、k空間中心からの距離と時間の関数として、サンプルされるべきk空間位置の平均数量を決定する。
さらなる変形においては、本発明に従って、ラジアルタイプの軌道が適用され得る。ここで、k空間軌道は、k空間において外側に進み、次に中心に戻り、そして、別の角度に対して継続する。連続したラジアル軌道間の角度は、望ましくは黄金角(”Golden Angle”)に等しくされ得る。または、そうでなければ、以前の軌道でなく又は以前の軌道によって希薄にしか訪問されていない方向に、新たな軌道それぞれをガイドするように選択され得る。3次元取得において、新たなラジアル軌道それぞれの方向は、2次元の黄金手段に従って選択され得る2つの角度によって定められ得る。Magnetic Resonances in Medicine 61:354−363(2009)における”Temporal Stability of Adaptive 3D Radial MRI Using Multidimensional Golden Means”に記述された方法と同様である。さらに、「軌道長(”trajectory length”)」(つまり、軌道がk空間中心に戻る以前にk空間中心から周辺へ進行している軌道によってカバーされる距離)は、ラジアル軌道間でランダムに変化してよい。
本発明の望ましい実施例に従って、MR信号データは、異なる空間感度プロフィールを有する2つまたはそれ以上のRF受信アンテナを介する並行画像化(parallel imaging)手段によって取得される。並行画像化は、取得プロシージャのさらなるスピードアップのために使用され得る。MR画像データの再構成は、MR信号データの規則的k空間グリッドへの再グリッド化を含み得る。例えば、従来のオペレーションの適用によるものであり、ここでは、コンボリューションカーネル(convolution kernel)が適切に選択される。このアプローチの基本的なアイデアは、不規則k空間位置から規則的グリッドへ直接的コンボリューションを適用することである。コンボリューションカーネルは、(事前に取得された)空間感度プロフィールに基づいて、それぞれのk空間位置に対して最適化され得る。カーネルは、既知のGRAPPA再構成スキームと同様に、密集してスキャンされたk空間の領域を使用して取得され得る。コンボリューションオペレーションは、一式のN個RFアンテナデータセット(Nは、MR信号取得のために使用されるRF受信アンテナの数量)から一つの単一な再グリッドk空間データセットへの変換を実行する。次のステップとして、再グリッドMR信号データは、フーリエ変換され、かつ、非線形ノイズフィルタリングが適用され得る(例えば、”Nonlinear Anisotropic Filtering of MR Data”、Guido Gerig et al、IEEE Transaction on Medical Imaging、11(2)、221−232、1992、または、”Noise Reduction for Magnetic Resonance Images via Adaptive Multiscale Products Thresholding”、Paul Bao and Lei Zhang、IEEE Transaction on Medical Imaging、22(9)、1089−1099、2003、もしくは、国際特許出願公開第WO97/23844号公報、において記述されるように)。そうしたフィルタリングアルゴリズムは、本発明に従って、希薄な変換領域への変換を暗黙のうちに含んでいる。MR画像の再構成は、不可欠なステップを伴う繰返しアルゴリズムとして実施され得る。ステップは、1.再グリッド化(特定的に計算されたカーネルを使用する)、2.フーリエ変換、3.非線形フィルタリング、および4.ステップ1−3の繰返し、である。
これまで説明された本発明に係る方法は、以下のMR機器によって実行され得る。MR機器は、検査ボリュームの中で均一で安定した磁場を生成するための少なくとも一つのメイン磁石コイルと、検査ボリュームの中で異なる空間的方向における切替え磁場勾配を生成するための数多くの勾配コイルと、検査ボリュームの中でRFパルスを生成、及び/又は、検査ボリュームの中に配置された対象物からMR信号を受信するための少なくとも一つのRFコイルと、時間的に一連のRFパルスおよび切替え磁場勾配をコントロールするためのコントロールユニットと、再構成ユニットと、および、可視化ユニットと、を含むものである。本発明に係る方法は、望ましくは、MR機器の再構成ユニット、可視化ユニット、及び/又は、コントロールユニットの対応するプログラミングによって実施される。
本発明に係る方法は、現時点で医療用途での大部分のMR機器において有利に実行され得る。この目的のために、コンピュータプログラムを利用することが、単に必要である。それにより、MR機器は、上記に説明された本発明に係る方法のステップを実施するようにコントロールされる。コンピュータプログラムは、データキャリア上に在るか、または、MR機器のコントロールユニットの中にインストールするためにダウンロードされるようにデータネットワーク内に在るか、いずれでもよい。
付属の図面は、本発明の望ましい実施例を開示する。しかしながら、図面は、説明目的だけのためにデザインされたものであり、本発明に係る限定を定めるものではないことが理解されるべきである。
図1は、本発明を実行するためのMR機器を示している。 図2は、本発明に従った、ブロックダイヤグラムとしてのMR画像の再構成のための繰返し規則(iteration scheme)である。 図3は、本発明に従って、不規則k空間サンプリング(irregular k−space sampling)を示している。 図4は、本発明に従って、空間次元(spatial dimension)と周波数次元(frequency dimension)を有するMR画像データの実施例に係るダイヤグラムを示している。 図5は、本発明に従って、適用可能な非デカルトk空間サンプリングスキーム(non−Cartesian k−space sampling scheme)を模式的に示している。
図1を参照すると、MR機器1が示されている。本機器は、超伝導または抵抗性のメイン磁石コイル2を含む。検査ボリュームを通じてz軸に沿って、実質的に均一で、時間的に一定なメイン磁場が形成されるようにである。
磁気共鳴生成及び操作システムは、一連のRFパルスと切替え(switched)磁場勾配を適用する。MR画像化を実行するように、核磁気スピンを転化または励起し、磁気共鳴を誘起し、磁気共鳴に再フォーカスし、磁気共鳴を操作し、空間的及び別の方法で磁気共鳴をエンコードし、スピンを飽和し、等をなすためである。
より特定的には、勾配パルス増幅器3は、勾配コイル4、5、および6のうち選択された一つに対して、検査ボリュームのx、y、およびz軸に沿って、電流パルスを適用する。デジタルRF周波数送信器7は、RFパルスまたはパルスパケットを、送信/受信スイッチ8を介して、全身ボリュームRFコイル9へ送信する。RFパルスを検査ボリュームの中へ送信するためである。典型的なMR画像化シーケンスは、お互いに一緒に取得される短期(short duration)のRFパルスセグメントのパケットで構成されており、そして、あらゆる適用された磁場勾配は、核磁気共鳴の選択された操作を達成する。RFパルスは、飽和、共鳴励起、磁化反転、共鳴再フォーカス、のために、または、共鳴を操作して検査ボリュームの中に配置された体10の一部分を選択するために使用される。MR信号は、また、全身ボリュームRFコイル9によってもピックアップされる。
体10の限定された領域のMR画像を生成するために、一式のローカルアレイRFコイル11、12、13が、画像化のために選択された領域に隣接して配置されている。アレイコイル11、12、13は、並行画像化のために使用され得るものであり、体−コイルRF送信によって誘発されたMR信号を受信する。
結果として生じるMR信号は、全身ボリュームRFコイル9によって、及び/又は、アレイコイル11、12、13によってピックアップされ、望ましくは増幅器(図示なし)を含んでいる受信器14によって復調される。受信器14は、送信/受信スイッチ8を介して、RFコイル9、11、12、および13に接続されている。
ホストコンピュータ15は、本発明に係る画像化シーケンスを生成するように、勾配パルス増幅器3と送信器7をコントロールする。受信器14は、各励起パルスに続く異なるエコー時間値におけるMR信号を受信する。データ取得システム16は、受信した信号のアナログ−デジタル変換を実行して、各MRデータラインをさらなる処理に好適なデジタルフォーマットに変換する。現在のMR機器において、データ取得システム16は、生の画像データの取得に特化された別個のコンピュータである。最後に、デジタルの生画像データは、好適な再構成アルゴリズムを適用する再構成プロセッサ17によって、画像表現へと再構成される。MR画像は、患者を貫く平面スライス、平行な平面スライスのアレイ、3次元ボリューム、等を表現することができる。画像は、次に、画像メモリの中に保管される。画像メモリは、スライス、プロジェクション、または、画像表現の他の部分を視覚化のために好適なフォーマットへと変換するためにアクセスされ得る。例えば、結果として生じるMR画像の人間が読取り可能な表示を提供するビデオモニタリング18を介するものである。
本発明に係るMR画像取得および再構成技術が、図2において説明される。
本方法は、MR信号データを取得するために画像化シーケンスに対して体の部分を投入することによって、ステップ41で開始する。ここで、MR信号データは、k空間のサブサンプリングを用いた不規則k空間サンプリングパターンを使用することにより、(k、k、k、t)空間におけるデータとして、つまり、k空間位置とエコー時間の関数として、取得される。MR信号では、RF受信アンテナ11、12、13を介して取得される。次のステップ42として、(k、k、k、t)空間における規則的グリッド(regular grid)へのMR信号データの再グリッド化(re−gridding)が実行される。その後で、ステップ43において、(x、y、z、ω)空間におけるMR画像データが取得されるように、再グリッド化されたMR信号データのフーリエ変換が実行される。ステップ44においては、非線形フィルタ(例えば、L1−ノルムといったような、従来の圧縮センシング(Compressed Sensing)におけるもの)がフーリエ変換されたデータに対して適用される。非線形フィルタリングは、圧縮センシングによって要求されるので、暗黙のうちに、希薄(sparse)領域への変換を実行する。画像の再構成は、ステップ44の後で終了してよい。しかしながら、k−t−空間の取得スキームおよび取得されたMR信号データの品質に応じて、繰返しが有益であり得る。この目的のために、ステップ45において、フィルタされたデータが、逆フーリエ変換によって(k、k、k、t)空間に戻すように変換される。一式の中間MR信号データを取得するためである。これらの中間MR信号データは、最初に取得されたMR信号データセットの(k、k、k、t)空間におけるサンプリング位置について補間される。最初に取得されたMR信号データから補間された中間MR信号データが引き算されて、一式の異なるMR信号データが得られる。そして、次の繰返しが開始され、再び異なるMR信号データが再グリッド化される。フーリエ変換43および非線形フィルタリング44に係るステップが、異なるMR信号データを再グリッド化するために繰り返される。ここで、ステップ44においてノイズ低減フィルタが適用される前に、現在および以前の繰返しに係るフーリエ変換されたデータが足し算される。数多くの繰返しの後で、最終的なMR画像データ46が(x、y、z、ω)空間で得られる。
本発明に従って、取得されたMR信号データの正確なk空間位置がランダム化される。サブサンプリングによるあらゆるエイリアシング(aliasing)アーチファクトは、このようにしてインコヒーレント(incoherent)にされる。インコヒーレントなアーチファクト(または少なくともその部分)の除去が、画像データが希薄な領域において、本発明に従って実行される。その領域においては、既に上述のように、非線形フィルタリングステップを使用して、アーチファクトが有益な情報から分離される。異なるタイプのそうした「希薄化(”sparsifying”)」変換が適用可能である。簡潔な変換は、例えば、MR画像データからエッジ画像を本質的に作成する勾配オペレーションである。ここでは、すべてのピクセル値が、それぞれのピクセル値と一つまたはそれ以上の隣接するピクセルの値との間の差異によって置き換えられる。このコンテクストにおいては、MR画像データの「希薄化」変換が、周波数次元においては、空間次元とは異なって動作することが重要である。MR画像データのエッジ(つまり、隣接する画像特徴の境界)は、空間次元において希薄であると推測され、一方で、MR画像データ自身は、周波数次元において希薄であると推測される。各空間位置において、MR画像データは、非常に限定された一式の周波数値に対してだけ非ゼロであることが期待される(図4に見られるように)。MR画像データが非ゼロである周波数座標は、時々の突然の変化を除いて、空間座標において典型的には連続関数である(例えば、手術機器といった、金属パーツに近いもの)。そうした変化は、おそらく既定の値と一致している。例えば、3.4ppmである(つまり、脂肪プロトンと水プロトンとの周波数における差異である)。さらに、これらの変化は、おそらく空間次元におけるエッジと一致している。このように、MR画像をエッジ画像へと変換する変換(例えば、各ボクセル位置において画像勾配を計算することによるもの)が、空間次元において使用され得るし、一方で、単一(unity)変換が、周波数次元において使用される。さらに精巧な「希薄化」変換が、もちろん、考えられる。
図3と図5は、本発明に従って使用され得る、不規則な非デカルトk空間サンプリングパターンの実施例を概説する。不規則(ランダム)k空間サンプリングスキームは、制限された空間−時間的サンプリング密度を用いて有利に適用され得る。再構成問題の調整に関して、k空間サンプリングパターンを選択することが有利である。パターンが、MR信号データを取得するためのk空間領域を再訪問(re−visit)するk空間の軌道を含むようにである(ここで、それぞれのk空間領域の「訪問(”visit”)」は時間において実質的に異なるものである)。サンプリング密度の制限が定められてよく、k空間中心からの距離と時間の関数としてサンプルされるべきk空間位置の平均数量を決定する。図5aは、k空間において外側に進み、次に内側に再び進行しているスパイラルまたはスパイラルのような軌道を示している。スパイラルは、k空間の外側領域が内側よりも少なく訪問されるように決定される。k空間サンプリングの軌道は、正確なスパイラルではなく、ランダムな摂動を伴うスパイラルである。図5bは、ラジアルタイプの軌道を示している。ここで、それぞれのk空間軌道(「スポーク(”spoke”)」)は、k空間において外側に進み、次に中心に戻り、そして、別の角度に対して継続する。連続したスポーク間の角度は、黄金角(”Golden Angle”)に等しい。各新たなスポークが、以前のスポークによっては訪問されなかった方向においてガイドされるようにである。図5Bにおける数字は、スポークの時間的な取得順序を示している。
図4は、仮想事例を示している。そこで、k空間サンプリングパターンは(例えば、図3に示されるように)、無限に密集している。この場合に、MR信号データは、(k、k、k、t)空間における規則的なグリッドに対してグリッド化され得る(tは、例えば、MR信号がスピンエコーまたは勾配エコーとして生成される場合にエコー時間であってよい)。そして、(x、y、z、ω)空間へとフーリエ変換する。結果として生じるMRデータは、従って、空間次元と周波数次元(ω)を含んでいる。図6は、明確のために、(x、ω)次元だけを示している。実際には、もちろん、k空間スキャン軌道は、無限の密集程度について、(x、y、z、ω)空間をサンプルしない。このように、データの普通のグリッド化およびフーリエ変換は、非常に強力なアーチファクトを有するMR画像データを結果として生じる。こうしたアーチファクトを抑圧するために、例えば、圧縮センシングに従って、上述のように、変換領域におけるMR画像データの希薄性が利用される。
多くのMR画像化アプリケーションにおいては、画像化スピードが重要である。しかしながら、MR画像の再構成のために要求されるMR信号が収集され得るところのスピードは、基本的に、物理的および生理学制約によって制限されている。従って、MR画像化領域における多くの最近の開発は、再構成されたMR画像の品質を劣化することなく、取得される信号データの量を削減することを目指している。多くのそうした開発の中で、圧縮センシング(Compressed Sensing(CS))の理論は、顕著な信号データ削減のための多大なポテンシャルを有している。CS理論においては、変換領域において希薄な表現を有する信号データセットが、適切な正規化(regularisation)アルゴリズムの適用によって、サブサンプルされた測定から回復され得る。サブサンプリングの可能性は、著しく削減された取得時間を導く。信号サンプリングと再構成のための数学的フレームワークとして、k空間サンプリング密度がナイキスト基準(Nyquist criterion)よりはるかに低い場合でさえ、信号データセットが正確に又は少なくとも高画像品質を伴って再構成され得る条件を規定する。そして、また、そうした再構成の方法も提供する。大多数の既存のCS基盤のMR取得および再構成スキームにおいては、基本的CS定式が使用される。MR信号データは変換領域において希薄であるという必須条件だけを利用するものである。例えば、M.Lustigらは、迅速なMR画像化のためのCS適用を提案してきた(M.Lustig et al.:”Sparse MRI:The Application of Compressed Sensing for Rapid MR Imaging”、Magnetic Resonances in Medicine、58、1182−1195、2007)。所定のアプリケーションにおいては知られていないMR画像に関する追加の従来情報が利用可能であり得るので、この従来情報をCS再構成の中に組み入れることが有利であることも、また、知られている。この論文は、圧縮センシングが(医療診断)MR信号データに対して適用され得ることを述べている。特に、ランダムに選択された位相エンコードラインを用いたインコヒーレントなサンプリング軌道が論じられている。さらに、圧縮サンプリングと並行画像化の組み合せに言及している。

Claims (10)

  1. MR機器の検査ボリュームの中に配置された対象物のMR画像化方法であって、
    MR信号データを取得するために画像化シーケンスに前記対象物をかけるステップであり、前記MR信号データは、k空間のサブサンプリングを用いた不規則k空間サンプリングパターンを使用することにより、k空間位置と時間の関数として取得される、ステップと、
    前記MR信号データからMR画像データを再構成するステップであり、前記MR画像データは空間次元と周波数次元を含み、前記MR画像データにおけるサブサンプリングアーチファクトを抑圧するために、変換領域における前記MR画像データの希薄性が利用される、ステップと、
    を含む、方法。
  2. 前記MR信号データは、異なる空間感度プロフィールを有する2つまたはそれ以上のRF受信アンテナを介する並行画像化によって取得される、
    請求項1に記載の方法。
  3. 前記再構成は、前記MR信号データの規則的k空間グリッドへの再グリッド化を含む、
    請求項2に記載の方法。
  4. 前記k空間サンプリングパターンは、MR信号取得の最中に、ランダムまたは擬似ランダムに磁場勾配を変化させることより生成される、連続したランダムk空間軌道を含む、
    請求項1乃至3いずれか一項に記載の方法。
  5. 前記MR信号データの取得のために、k区間の非デカルトサンプリングが適用される、
    請求項1乃至4いずれか一項に記載の方法。
  6. 前記k空間サンプリングパターンは、前記MR信号データを獲得するためにk空間領域を再訪問するk空間軌道を含む、
    請求項1乃至5いずれか一項に記載の方法。
  7. 前記MR画像データは、圧縮センシングによって再構成される、
    請求項1乃至6いずれか一項に記載の方法。
  8. 前記MR画像データの再構成は、前記MR画像データが希薄である前記変換領域への変換を含み、前記周波数次元における前記変換は、前記空間次元における前記変換とは異なって動作する、
    請求項1乃至7いずれか一項に記載の方法。
  9. 請求項1乃至8いずれか一項に記載の方法を実行するためのMR機器であって、
    前記MR機器は、
    検査ボリュームの中で均一で安定した磁場を生成するための少なくとも一つのメイン磁石コイルと、
    前記検査ボリュームの中で異なる空間的方向における切替え磁場勾配を生成するための複数の勾配コイルと、
    前記検査ボリュームの中でRFパルスを生成、及び/又は、前記検査ボリュームの中に配置された対象物からMR信号を受信するための少なくとも一つのRFコイルと、
    時間的に一連のRFパルスおよび切替え磁場勾配をコントロールするためのコントロールユニットと、
    再構成ユニットと、を含み、
    前記MR機器は、
    MR信号データを取得するために画像化シーケンスに前記対象物をかけるステップであり、前記MR信号データは、k空間のサブサンプリングを用いた不規則k空間サンプリングパターンを使用することにより、k空間位置と時間の関数として取得される、ステップと、
    前記MR信号データからMR画像データを再構成するステップであり、前記MR画像データは空間次元と周波数次元を含み、前記MR画像データにおけるサブサンプリングアーチファクトを抑圧するために、変換領域における前記MR画像データの希薄性が利用される、ステップと、
    を実行するように構成されている、MR機器。
  10. MR機器で実行されるコンピュータプログラムであって、
    MR信号データを取得するために画像化シーケンスを生成するステップであり、前記MR信号データは、k空間のサブサンプリングを用いた不規則k空間サンプリングパターンを使用することにより、k空間位置と時間の関数として取得される、ステップと、
    前記MR信号データからMR画像データを再構成するステップであり、前記MR画像データは空間次元と周波数次元を含み、前記MR画像データにおけるサブサンプリングアーチファクトを抑圧するために、変換領域における前記MR画像データの希薄性が利用される、ステップと、
    を実行するためのインストラクションを含む、
    コンピュータプログラム。
JP2016503743A 2013-03-21 2014-03-12 圧縮センシングを使用したmr画像の再構成 Withdrawn JP2016512780A (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201361803929P 2013-03-21 2013-03-21
US61/803,929 2013-03-21
PCT/IB2014/059660 WO2014147518A2 (en) 2013-03-21 2014-03-12 Mr image reconstruction using compressed sensing

Publications (1)

Publication Number Publication Date
JP2016512780A true JP2016512780A (ja) 2016-05-09

Family

ID=50624888

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016503743A Withdrawn JP2016512780A (ja) 2013-03-21 2014-03-12 圧縮センシングを使用したmr画像の再構成

Country Status (5)

Country Link
US (1) US9964615B2 (ja)
EP (1) EP2976654A2 (ja)
JP (1) JP2016512780A (ja)
CN (1) CN105074491B (ja)
WO (1) WO2014147518A2 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018051649A1 (ja) * 2016-09-13 2018-03-22 株式会社日立製作所 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
JP2020085494A (ja) * 2018-11-16 2020-06-04 ヤフー株式会社 情報処理装置、情報処理方法、及び情報処理プログラム

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2824472A1 (en) * 2013-07-10 2015-01-14 Siemens Aktiengesellschaft MRI with self-navigated motion compensation
CN105259525B (zh) * 2015-10-28 2019-08-13 北京大学 基于邻域共享压缩感知的动态对比增强磁共振快速成像方法
US10908235B2 (en) * 2016-04-08 2021-02-02 The Johns Hopkins University Method of fast imaging of NMR parameters with variably-accelerated sensitivity encoding
US11041926B2 (en) * 2016-06-02 2021-06-22 Koninklijke Philips N.V. Dixon-type water/fat separation MR imaging
CN106093814A (zh) * 2016-06-02 2016-11-09 浙江理工大学 一种基于多尺度低秩模型的心脏磁共振成像方法
CN106569158B (zh) * 2016-10-24 2019-02-05 浙江大学 一种降低特殊吸收率的压缩型螺旋梯度射频脉冲设计方法
US10145924B2 (en) 2017-01-13 2018-12-04 Siemens Healthcare Gmbh Magnetic resonance imaging with undersampling
CN110226100B (zh) * 2017-01-25 2022-03-15 上海联影医疗科技股份有限公司 用于磁共振成像的系统和方法
US10775464B2 (en) * 2017-01-31 2020-09-15 Regents Of The University Of Minnesota System and method for dynamic, cardiac phase-resolved quantitative longitudinal relaxation parameter mapping
JP2021532875A (ja) 2018-07-30 2021-12-02 ハイパーファイン,インコーポレイテッド 磁気共鳴画像再構成のための深層学習技術
JP2021534852A (ja) 2018-08-15 2021-12-16 ハイパーファイン,インコーポレイテッド 磁気共鳴画像においてアーティファクトを抑制するためのディープラーニング技術
US10672151B1 (en) * 2019-01-07 2020-06-02 Uih America, Inc. Systems and methods for magnetic resonance image reconstruction
JP2022526718A (ja) 2019-03-14 2022-05-26 ハイパーファイン,インコーポレイテッド 空間周波数データから磁気共鳴画像を生成するための深層学習技術
EP3748385A1 (en) 2019-06-04 2020-12-09 Koninklijke Philips N.V. Optimized k-space profile ordering for 3d radial or spiral mr imaging
CN110544254A (zh) * 2019-09-12 2019-12-06 杭州海睿博研科技有限公司 用于核磁同步自动化成像与分割的系统和方法
CN113254214B (zh) * 2021-06-08 2023-03-10 西安科技大学 基于OpenMP的倒装焊芯片声时频域和时域成像方法

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11501140A (ja) 1995-12-21 1999-01-26 フィリップス エレクトロニクス エヌ ベー 方向適合ノイズ減少
US6437568B1 (en) * 2000-10-02 2002-08-20 General Electric Company Low noise MRI scanner
WO2005078470A2 (en) * 2004-02-10 2005-08-25 Koninklijke Philips Electronics N.V. Magnetic resonance imaging method
DE102006042998B4 (de) * 2006-09-13 2008-07-03 Siemens Ag Messsequenz für die dreidimensionale MR-Bildgebung sowie MR-Gerät
US20100145185A1 (en) * 2008-12-04 2010-06-10 Xiaole Hong Extending the resolution of mri data by combining subsets from plural image acquisitions
CN102264291B (zh) * 2008-12-26 2013-10-23 株式会社日立医疗器械 磁共振成像装置及读出梯度磁场误差修正方法
US8310233B2 (en) 2009-02-18 2012-11-13 Mayo Foundation For Medical Education And Research Method for image reconstruction from undersampled medical imaging data
US8897515B2 (en) 2009-09-08 2014-11-25 Mayo Foundation For Medical Education And Research Method for compressed sensing image reconstruction using a priori knowledge of spatial support
EP2402780A1 (en) * 2010-06-23 2012-01-04 Koninklijke Philips Electronics N.V. Method of reconstructing a magnetic resonance image of an object considering higher-order dynamic fields
US20130099786A1 (en) * 2010-07-02 2013-04-25 Koninklijke Philips Electronics N.V. Parallel magnetic resonance imaging using undersampled coil data for coil sensitivity estimation
CN101975936A (zh) * 2010-09-03 2011-02-16 杭州电子科技大学 一种基于cs压缩感知技术的快速磁共振成像方法
CN102018514B (zh) * 2010-12-30 2015-07-08 中国科学院深圳先进技术研究院 磁共振弥散张量成像方法及系统
EP2503349A1 (en) * 2011-03-24 2012-09-26 Koninklijke Philips Electronics N.V. MR image reconstruction using prior information-constrained regularization
CN102540116A (zh) * 2011-12-08 2012-07-04 中国科学院深圳先进技术研究院 磁共振成像方法和系统
CN102565737A (zh) * 2011-12-12 2012-07-11 中国科学院深圳先进技术研究院 快速磁共振成像方法及系统

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018051649A1 (ja) * 2016-09-13 2018-03-22 株式会社日立製作所 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
JP2020085494A (ja) * 2018-11-16 2020-06-04 ヤフー株式会社 情報処理装置、情報処理方法、及び情報処理プログラム
JP7036704B2 (ja) 2018-11-16 2022-03-15 ヤフー株式会社 情報処理装置、情報処理方法、及び情報処理プログラム

Also Published As

Publication number Publication date
EP2976654A2 (en) 2016-01-27
CN105074491B (zh) 2018-07-17
WO2014147518A3 (en) 2014-11-27
US9964615B2 (en) 2018-05-08
US20160291106A1 (en) 2016-10-06
WO2014147518A2 (en) 2014-09-25
CN105074491A (zh) 2015-11-18

Similar Documents

Publication Publication Date Title
US9964615B2 (en) MR image reconstruction using compressed sensing
JP6513398B2 (ja) 事前情報に制約される正則化を使用するmr画像再構成
RU2557334C2 (ru) Динамическая контрастная улучшенная мр визуализация с реконструкцией сжатого измерения
JP6243522B2 (ja) 正則化された検出再構成を使用するマルチエコーディクソン水−脂肪分離及びb0歪補正による並列mri
CN106574954B (zh) 针对epi的具有奈奎斯特伪影校正的并行mr成像
JP6270709B2 (ja) 磁気共鳴生データの再構成方法および装置
US10203394B2 (en) Metal resistant MR imaging
US10823806B2 (en) Magnetic resonance imaging using dataset undersampling
CN107209238B (zh) 具有对边带伪迹的抑制的并行多切片mr成像
EP3080634B1 (en) Zero echo time mr imaging with water/fat separation
JP6684781B2 (ja) ゼロエコー時間mrイメージング
CN106796274B (zh) 具有伪迹抑制的propeller-mr成像
JP2016509935A (ja) パラレルのマルチスライスmr画像法
CN113544526B (zh) 使用波编码的并行mr成像
EP4239357A1 (en) Diffusion mapping by mr fingerprinting
US11815577B2 (en) Parallel MR imaging using wave-encoding

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150925

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20170310

A761 Written withdrawal of application

Free format text: JAPANESE INTERMEDIATE CODE: A761

Effective date: 20170718