JP6636461B2 - マルチエコーk空間収集を使用するMRイメージング - Google Patents

マルチエコーk空間収集を使用するMRイメージング Download PDF

Info

Publication number
JP6636461B2
JP6636461B2 JP2016569024A JP2016569024A JP6636461B2 JP 6636461 B2 JP6636461 B2 JP 6636461B2 JP 2016569024 A JP2016569024 A JP 2016569024A JP 2016569024 A JP2016569024 A JP 2016569024A JP 6636461 B2 JP6636461 B2 JP 6636461B2
Authority
JP
Japan
Prior art keywords
space
echo
signal
sampled
sequence
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
JP2016569024A
Other languages
English (en)
Other versions
JP2017516551A5 (ja
JP2017516551A (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.)
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 JP2017516551A publication Critical patent/JP2017516551A/ja
Publication of JP2017516551A5 publication Critical patent/JP2017516551A5/ja
Application granted granted Critical
Publication of JP6636461B2 publication Critical patent/JP6636461B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • 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/50NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo 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
    • 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]
    • 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/5616Echo 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 gradient refocusing, e.g. EPI
    • 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/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/5618Echo 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 both RF and gradient refocusing, e.g. GRASE
    • 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/56509Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling

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)

Description

本発明は、磁気共鳴(MR)イメージングの分野に関する。本発明は、物体のMRイメージング方法に関する。本発明は更に、MRデバイスと、MRデバイス上で実行されるコンピュータプログラムとに関する。
今日では、2次元又は3次元画像を形成するために磁場と核スピンとの相互作用を利用する画像形成MR方法は、軟組織のイメージングにおいて、他のイメージング方法よりも多くの側面において優れ、電離放射線を必要とせず、また、通常、侵襲的ではないため、特に医療診断分野において、広く使用されている。
MR方法では、一般に、例えば検査される患者の身体である物体は、強力な均一磁場内に配置される。当該磁場の方向は、同時に、測定の基準となる座標系の軸(通常は、z軸)を規定する。磁場は、所定の周波数(いわゆるラーモア(Larmor)周波数又はMR周波数)の交流電磁場(RF場)の印加によって励起(スピン共鳴)可能である磁場強度に依存して、個々の核スピンに対し様々なエネルギーレベルを生成する。巨視的観点では、個々の核スピンの分布は、適切な周波数の電磁パルス(RFパルス)の印加によって、平衡状態から外れる全般的な磁化を生成し、これにより、磁化は、z軸周りの歳差運動を行う。歳差運動は、円錐面を描き、その開口角は、フリップ角と呼ばれる。フリップ角の大きさは、印加される電磁パルスの強度及び持続時間に依存する。いわゆる90°パルスの場合、スピンは、z軸から横断面へと偏向される(フリップ角90°)。
RFパルスの終了後、磁化は、第1の時定数T(スピン格子又は縦緩和時間)で、z方向における磁化が再び蓄積される元の平衡状態に緩和して戻り、また、z方向に垂直な方向における磁化は、第2の時定数T(スピンスピン又は横緩和時間)で緩和する。磁化の変化は、当該磁化の変化がz軸に垂直な方向において測定されるように、MRデバイスの検査ボリューム内に配置され、方向付けられている受信RFコイルによって検出される。横方向磁化の減衰は、例えば90°パルスの印加後、同じ位相を有する秩序状態から、すべての位相角が均一に分布している(ディフェージング)状態への核スピンの遷移(局所的な磁場の不均一性によって誘起される)を伴う。ディフェージングは、リフォーカシングパルス(例えば180°パルス)によって相殺される。これは、受信コイルにおいてエコー信号(スピンエコー)を生成する。
身体内における空間分解能を実現するために、3つの主軸に沿って延在する定磁場勾配が、均一磁場に重畳され、スピン共鳴周波数の線形空間的依存性につながる。このとき、受信コイルにおいて捕捉される信号は、身体内の様々な場所に関連付けることができる様々な周波数成分を含む。受信コイルを介して得られる信号データは、空間周波数領域に対応し、k空間データと呼ばれる。k空間データは、通常、異なる位相エンコーディングで収集される複数のラインを含む。各ラインは、幾つかのサンプルを集めることによって、2値化される。k空間データのセットは、画像再構成アルゴリズムによって、MR画像に変換される。
いわゆる「ターボスピンエコー」(TSE)シーケンスは、よく知られたマルチエコーイメージングシーケンスである。TSEシーケンスの1つの「ショット」は、複数の高速印加される(通常180°の)リフォーカシングRFパルスが後に続く、磁気共鳴励起のための初期RFパルスを含み、異なって位相エンコードされた一連のスピンエコー信号が生成される。エコー信号が収集され、各エコー信号は、k空間ライン、即ち、k空間の1次元サンプルを表し、k空間ラインの位置は、シーケンスの印加された周波数エンコーディングと位相エンコーディングとが切替わった磁場勾配によって決定される。いわゆる「ターボ係数」(TF)は、各励起後に収集されたエコーの数である。通常、収集されたk空間プロファイルからMR画像を再構成することができるようにするために、マルチエコーシーケンスの複数のショットが、k空間を完全にサンプリングするために印加される。k空間は、通常、各セグメントが、シーケンスの各ショットによって生成される一連のエコーにおける1つのエコー指数に帰するようにセグメント化される。つまり、各ショットは、通常、k空間のすべてのセグメントからk空間プロファイルを収集し、k空間におけるセグメント数は、TFに等しい。
今日では、TSEシーケンスが、ほぼすべての応用において適用されている。その幅広い使用にも関わらず、セグメント化されたk空間の分散収集によって、動きに対して非常に敏感であることが知られている。例えば皮下脂肪の撮像物体の小さな動きが、関心部分に重なるゴーストを発生させ、走査の全体的な診断価値を妨げてしまう。
上記k空間セグメンテーションの原理及び問題は、スピンエコー(TSE)及び勾配エコーシーケンス(TFE)と、「エコープラーナ―イメージング」−EPI−シーケンスとの両方に当てはまる。
米国特許第8,155,389号から、位相エンコーディングラインのセット間の相関関係から動き情報を導出することが知られている。
米国特許出願公開第2008/061779号は、3次元MRイメージング用の磁気共鳴データ収集に関する。このデータ収集は、部分収集の分散である(3Dk空間における)kラインを採用する。つまり、各部分収集は、(周波数エンコーディング方向における)幾つかのkxプロファイルを含む。各kプロファイルは、異なる(k,k)2次元位相エンコーディングを有する。特に、(k,k)平面では、ラジアル走査アプローチが採用される。測定時間全体に亘って、kプロファイルの均等分布が、中心及び周辺k空間において(2次元位相エンコーディング平面に亘って)採用される。
上記から、改良されたMRイメージング技術が必要であることは、すぐに理解できる。したがって、本発明は、動きに対し、よりロバストである高速MRイメージングを可能にすることを目的とする。
本発明によれば、MRデバイスの検査ボリューム内に配置された物体のMRイメージング方法が開示される。当該方法は、
物体に、イメージングシーケンスを与えることによって、エコー信号を発生させるステップと、
エコー信号を収集するステップと、
収集したエコー信号からMR画像を再構成するステップとを含み、各エコー信号は、k空間ラインに帰し、k空間の一部において隣接して配置される幾つかのk空間ラインが、繰り返しサンプリングされ、当該幾つかのk空間ラインは、繰り返し毎に異なる順序でサンプリングされる。
k空間走査アプローチは、読み出しによってスパンする(周波数エンコーディング方向(k)及び位相エンコーディング方向(k))2次元k空間(k平面)に適用される。本発明によれば、k空間の一部からの幾つかのk空間ラインが、繰り返しサンプリングされる。これは、同じ幾つかのk空間ラインが、2回以上サンプリングされることを意味する。サンプリングは、例えば信号平均化のために繰り返される。様々な繰り返しにおいて、スピン条件又は磁化条件が、異なるように選択される。例えば平均化は、磁化の異なる、特に逆の位相において、繰り返し行われる。Dixon収集において、様々な繰り返しにおいて、異なるエコー時間が選択される。位相感応性反転収集では、様々な繰り返しが、それぞれ、スピン反転の適用あり及び適用なしで行われる。繰り返しサンプリングされる各k空間ラインは、動きアーチファクトを相殺するために、繰り返し毎に異なる順番でサンプリングされる。順番は、k空間ラインが連続的にサンプリングされる順序を規定する。順番は、収集したエコー信号の連続における各エコー信号の位置に対する各k空間ラインのk空間位置(即ち、位相エンコーディング)の帰属によって決定される。これは、つまり、k空間ラインがイメージングシーケンスによってサンプリングされる(例えばランダムな)順番は、繰り返し毎に変えられることを意味し、これは、上記幾つかのk空間ラインが収集される時間的な連続が、通常、k空間の対応する部分における対応する空間ラインの順次的な配列とは異なることを暗示する。このアプローチは、動きアーチファクトの特に効果的な平均化を達成する。
本発明にしたがって適用されるイメージングシーケンスは、シングルエコースピンエコー(SE)又は勾配エコー(GE)シーケンスである。代替実施形態では、イメージングシーケンスは、マルチエコーイメージングシーケンスであってよい。物体に、マルチエコーイメージングシーケンスの幾つかのショットが与えられ、各ショットによって、一連のエコー信号が発生する。この実施形態では、当該幾つかのk空間ラインが収集される順番は、各ショットによって発生される一連のエコー信号における各エコー信号の数を示すエコー指数に対する各k空間ラインのk空間位置の帰属によって決定される。
本発明の更なる好適な実施形態では、k空間は、中心k空間部と、1つ以上の周辺k空間部とに分割される。中心k空間部における信号平均値の数が、周辺k空間部における信号平均値の数よりも高い。この実施形態によれば、本発明は、k空間の中心部が、動きアーチファクトを相殺するために、信号平均化を目的とした多数の収集でサンプリングされる動きに反応しない部分平均化アプローチを提案する。異なる信号平均化が、中心k空間部と周辺k空間部とに、それぞれ適用される。周辺k空間部よりも中心k空間部においてより多い数の信号平均値を適用することによって、再構成されたMR画像において、動きアーチファクトは効率的に平均される。これは、中心k空間データが、画像全体のコントラストを主に決定するからである。
なお、信号平均化は、k空間において、本発明にしたがって、又は、画像空間において、例えば反復MR画像再構成技術を使用して行われる。
好適には、(中心及び/又は周辺k空間部の)少なくとも幾つかの繰り返しサンプリングされるk空間ラインは、順不同でサンプリングされる。ランダム化によって、収集されたk空間データの動きによって生じた不連続性が分布される。これにより、動きによって破損したk空間データの規則的なk空間間隔が回避され、したがって、再構成されたMR画像におけるゴースト発生が回避される。
本発明の更に別の好適な実施形態では、k空間の少なくとも一部が、可変密度でサンプリングされる。これは、つまり、収集されたk空間ラインは、位相エンコーディング方向において、不規則に分布されることを意味する。可変k空間密度は、k空間における動きによって破損したk空間サンプルの規則的分布を回避することによって、ゴースト発生を更に減少させる。
更に、k空間の少なくとも一部、好適には、周辺k空間部は、不完全にサンプリングされる。これは、各k空間領域におけるサンプリング密度が、イメージングされる視野に関して、ナイキスト制限を下回ることを意味する。
これは、例えば、動き最適化された可変密度の圧縮センシング(CS)による走査時間を短縮するために採用される。CS理論は、大幅な信号データ削減の大きな可能性を有することが知られている。CS理論では、変換領域においてスパーズ表現を有する信号データセットが、適切な正則化アルゴリズムの適用によって、アンダーサンプリングされた測定値から復元可能である。アンダーサンプリングの可能性は、収集時間を著しく短縮する。信号サンプリング及び再構成の数学的枠組みとして、CSは、信号データセットが正確に、又は、k空間サンプリング密度がナイキスト基準をはるかに下回る場合であっても、少なくとも高画質で再構成できる条件を規定し、また、そのような再構成のための方法も提供する。ほとんどの既存のCSベースのMR収集及び再構成スキームでは、MR信号データが変換領域において少ないという必須条件のみを利用する基本的なCS公式が使用される。
更に、本発明の方法は、それ自体は知られている「部分行列」技術と組み合わされる。「部分行列」は、k空間の複素共役性を利用する走査時間短縮方法である。負のk空間位置からのk空間プロファイルは、対応する正のk空間位置からの複素共役k空間プロファイルと実質的に同一であるため、本発明では、周辺の正又は負の部分のみがサンプリングされて、k空間全体のサンプリングが再現される。これは、最大分解能を維持しつつ、走査時間を大幅に短縮する。
本発明の好適な実施形態によれば、本発明のマルチエコーイメージングシーケンスは、スピンエコーシーケンスであり、そのうちの1つのショットは、
磁気共鳴の励起のための少なくとも1つのRFパルスと、
複数のリフォーカシングRFパルスと、
各エコー信号の周波数エンコーディングのための複数の切り替え磁場勾配と、
各エコー信号の位相コーディングのための複数の切り替え磁場勾配とを含む。つまり、本発明の方法は、TSE又はスピンエコータイプのEPIといったそれ自体は知られている「ターボ」スピンエコーシーケンスと組み合わせて適用される。或いは、マルチエコーイメージングシーケンスは、勾配エコーシーケンスであってもよく、そのうちの1つのショットは、
磁気共鳴の励起のための少なくとも1つのRFパルスと、
磁気共鳴をリフォーカスし、各エコー信号を周波数エンコーディングするための複数の切り替え磁場勾配と、
各エコー信号を位相エンコーディングするための複数の切り替え磁場勾配とを含む。つまり、本発明の方法は、TFE又はEPIといったそれ自体は知られている「ターボ」勾配エコーシーケンスと組み合わせて適用される。
本発明の方法は、デカルトk空間サンプリングスキームと組み合わされるのに特に適している。しかし、本発明の方法は、非デカルトサンプリングと組み合わせて適用されてもよい。あらゆるタイプのk空間サンプリングについて、本発明の意味において「中心k空間部」との用語は、k空間原点に近い又はk空間原点を取り囲むk空間プロファイルを指す。
更に別の好適な実施形態では、動き情報が、繰り返しサンプリングされたk空間ラインから収集されたエコー信号から導出される。例えば同じk空間位置から、異なる時点において収集されたエコー信号が関連付けられて、動き情報が抽出される。動き情報は、適切な後処理による収集された信号データの対応する動き補償のために利用される。或いは、検出された動きが、補償されるには重大過ぎることが分かると、動きによって破損したk空間プロファイルは廃棄するように自動的に決定されてもよい。したがって、動きアーチファクトは、本発明にしたがって平均されるか、破損データを補正又は廃棄することによって、完全に回避される。
MRイメージング分野の当業者には知られているように、多くの場合、水及び脂肪といった様々な化学種の信号全体への相対的寄与に関する情報を得て、これらの化学種の幾つかの寄与を抑制するか、これらの化学種の全ての寄与を別々に又は一緒に解析することが望まれる。これらの寄与は、異なるエコー時間において収集された2つ以上の対応するエコーからの情報が組み合わされると計算することができる。これは、化学シフトエンコーディングと考えられる。化学シフトエンコーディングでは、追加の次元、即ち、化学シフト次元が規定され、少し異なるエコー時間において2つ以上のエコー信号データセットを収集することによって符号化される。特に、水−脂肪分離のために、これらのタイプの実験は、しばしば、Dixon式測定と呼ばれる。Dixonイメージング又はDixon水/脂肪イメージングによって、水/脂肪分離は、異なるエコー時間において収集された2つ以上の対応するエコーから、水及び脂肪の寄与を計算することによって達成される。一般に、このような分離は、脂肪及び水において、水素の既知の歳差運動周波数差があることによって可能である。その最も単純な形において、水及び脂肪の画像は、「同位相(in phase)」及び「逆位相(out of phase)」データセットの加算か、減算かの何れかによって生成される。
好適な実施形態では、本発明の方法は、Dixonイメージングと組み合わせられる。したがって、つまり、イメージングシーケンスは、Dixonシーケンスであり、収集されたエコー信号への水及び脂肪からの信号寄与は分離される。本発明の特に動きロバストである収集スキームはDixonイメージングによく適し、再構成されたMR画像における動きアーチファクトが少ないことにより、診断画質が高くなることは、本発明の重要な見識である。Dixonイメージングの場合、同一k空間ラインの繰り返しサンプリングは、繰り返し毎に異なるエコー時間で行われることが有利である。
最も好適には、エコー信号は、Dixonシーケンスの各エコー時間において、異なるk空間ラインの順序で収集される。収集順序の異なるランダム化は、動きアーチファクトを最適に平均する。
これは、特に、いわゆる重み係数付加された(modulus-added weighted)脂肪及び水画像の再構成と組み合わせて適用され、水及び脂肪からの信号寄与は、MR画像を再構成するステップにおいて組み合わされる。重み係数付加された脂肪及び水画像を再構成するために、n乗された脂肪寄与の係数と、n乗された水寄与の係数との加重和のn乗根が、各画像ピクセル又はボクセルについて計算される。
より一般的に、本発明の方法は、k空間ラインが繰り返し収集され、繰り返し毎に、少なくとも1つのイメージングパラメータが変更される任意のタイプのイメージング技術と組み合わされてよい。したがって、本発明の方法は、反転回復技術(IR/PSIR)、拡散強調イメージング、磁化移動コントラスト(MTC)、化学交換飽和移動(CEST)等と有利に組み合わせることができる。
本発明の更なる態様では、位相感応性磁気共鳴画像が、収集された磁気共鳴信号の位相情報が再構成において保たれるという点において、再構成される。つまり、再構成された磁気共鳴画像のピクセルは、複素数の値である。
ここまで説明してきた本発明の方法は、検査ボリューム内に均一な静磁場を発生させる少なくとも1つの主磁石コイルと、検査ボリューム内の様々な空間方向において切り替え磁場勾配を発生させる幾つかの勾配コイルと、検査ボリューム内にRFパルスを発生させる及び/又は検査ボリューム内に配置される物体からのMR信号を受信する少なくとも1つのRFコイルと、RFパルス及び切り替え磁場勾配の時間的連続を制御する制御ユニットと、再構成ユニットとを含むMRデバイスによって実行される。本発明の方法は、例えばMRデバイスの再構成ユニット及び/又は制御ユニットの対応するプログラミングによって実施される。
本発明の方法は、今日、臨床使用されているほとんどのMRデバイスにおいて有利に実行可能である。このためには、MRデバイスが上記された本発明の方法のステップを行うように、当該MRデバイスが制御されるコンピュータプログラムを利用するだけでよい。コンピュータプログラムは、データ担体上にあるか、又は、MRデバイスの制御ユニット内へのインストールのために、ダウンロードされるようにデータネットワーク内にあってよい。
添付の図面は、本発明の好適な実施形態を開示する。しかし、当然ながら、図面は、例示目的のみにデザインされたものであり、本発明の制限の定義としてデザインされたものではない。
図1は、本発明の方法を実行するMRデバイスを示す。 図2は、本発明のサンプリングスキームの第1の実施形態を概略的に示すk空間の図を示す。 図3は、本発明のサンプリングスキームの第2の実施形態を概略的に示すk空間の図を示す。 図4は、本発明のサンプリングスキームの第3の実施形態を概略的に示すk空間の図を示す。 図5は、本発明のサンプリングスキームの第4の実施形態を概略的に示すk空間の図を示す。 図6は、本発明のサンプリングスキームの第5の実施形態を概略的に示すk空間の図を示す。 図7は、動きゴースト寄与を有する同位相画像、逆位相画像、水画像及び脂肪画像を概略的に示す。 図8は、Dixon測定における同位相及び逆位相のサンプリングを概略的に示すk空間の図を示す。
図1を参照するに、MRデバイス1が示される。デバイスは、実質的に均一で、時間的に一定の主磁場が、検査ボリュームを通るz軸に沿って作成されるように、超電導又は抵抗主磁石コイル2を含む。
磁気共鳴発生及び操作システムは、一連のRFパルス及び切り替え磁場勾配を印加して、核磁気スピンを反転又は励起し、磁気共鳴を誘導し、磁気共鳴をリフォーカスし、磁気共鳴を操作し、磁気共鳴を空間的及び他の方法でエンコードし、スピンを飽和させて、MRイメージングを行う。
より具体的には、勾配パルス増幅器3が、電流パルスを、検査ボリュームのx、y及びz軸に沿った全身勾配コイル4、5及び6のうちの選択された全身勾配コイルに印加する。デジタルRF周波数送信器7が、RFパルス又はパルスパケットを、送信/受信スイッチ8を介して、全身ボリュームRFコイル9に送信して、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)、エコーボリュームイメージング、勾配及びスピンエコーイメージング、高速スピンエコー(TSE)イメージング等といった複数のMRイメージングシーケンスのうちの何れかが生成される。選択されたシーケンスについて、受信器14は、各RF励起パルスの後に続いて、単一又は複数のMRデータラインを高速連続で受信する。データ収集システム16が、受信信号のアナログ‐デジタル変換を行い、各MRデータラインを、更なる処理に適したデジタル形式に変換する。最新式のMRデバイスでは、データ収集システム6は、ロー画像データの収集に特化した別個のコンピュータである。
最終的には、デジタルロー画像データは、再構成プロセッサ17によって、画像表現に再構成される。再構成プロセッサ17は、フーリエ変換、又は、他の適切な再構成アルゴリズムを適用する。MR画像は、患者の平面スライス、平行平面スライスのアレイ、3次元ボリューム等を表す。次に、画像は、画像メモリに記憶される。画像メモリにおいて、画像は、スライス、投影又は画像表現の他の部分を、例えばビデオモニタ18を介する可視化のための適切な形式に変換するためにアクセスされる。ビデオモニタ18は、結果として得られるMR画像の人間が読み取り可能な表示を提供する。
図1を引き続き参照し、更に図2乃至図5を参照するに、本発明の動きロバストな部分平均化アプローチの実施形態について説明する。
本発明によれば、身体10に対し、例えばTSEシーケンスといったマルチエコーイメージングシーケンスの幾つかのショットを行う。シーケンスの各ショットによって、一連のエコー信号が生成される。エコー信号は、収集され、各エコー信号は、k空間ラインに帰する。k空間ラインは、図2乃至図5において、k軸に平行に延在する細い水平ラインとして表される。信号平均化は、k空間ラインの少なくとも幾つかを繰り返しサンプリングすることによって行われる。繰り返しサンプリングするとは、同一のk空間ラインが、最終MR画像の再構成する信号収集処理中に、2回以上サンプリングされることを意味する。信号対ノイズ比の向上が、同じk空間位置からそれぞれ収集された幾つかのエコー信号の平均を取ることによって達成される。平均値の数は、収集の数(NSA)とも呼ばれる。信号対ノイズ比(SNR)の向上の近似量は、平均値の数の平方根として計算される。
図2に示されるように、k空間20は、中心k空間部21と、第1の周辺k空間部22と、第2の周辺k空間部23とに分割される。図示される実施形態では、中心k空間部21は、TSEシーケンスの各ショットによってサンプリングされ、第1の周辺k空間部22は、2番目のショット毎にサンプリングされ、第2の周辺k空間部23は、4番目のショット毎にサンプリングされる。提案される部分平均化アプローチによれば、k空間位置の可変の繰り返しサンプリングが行われ、中心k空間部21は、周辺部22、23に比べて、より頻繁に収集される。k空間ラインのその限定されたセットを有する中心k空間部21は、非常に頻繁にサンプリングされる。動き情報は、中心部21の繰り返しサンプリングされたk空間ラインから収集されたエコー信号から導出される。例えば、様々な平均化ステップ中に収集されたk空間データの相互相関が計算される。これに基づいて、自動動き補償が行われる。検出された動きが、重大過ぎることが分かると、対応するエコー信号データは拒絶される。撮像物体の位置が、動きによって不可逆的に変更されていることが分かると、走査パラメータが、位置シフトに従って自動的に更新される。k空間ラインのより大きいセットを有する中間周辺k空間部22は、それほど頻繁ではないが、走査されて、動きを平均し、アーチファイト及びT緩和といった信号減衰効果を相殺し、少なくともSNRの特定の向上を提供する。第2の周辺部23は、1回しか収集されない。第2の周辺部23は、例えばCS及び/又は部分行列技術が使用されるかどうかに依存して、1回未満で収集されてもよい。最後に、MR画像が、収集されたエコー信号から再構成される。この一般的なアプローチは、FFE、SE、TFE、TSE及びEPIシーケンスといったあらゆるタイプのシーケンスに適用可能である。
図3乃至図5に示される本発明の更なる態様は、中心k空間部21においては、密度の高いサンプリングで、周辺k空間部22、23においては、より密度の低いサンプリングで、k空間サンプリングが、可変密度で行われることである。図3乃至図5に示されるように、収集されたk空間ラインは、位相エンコーディング方向kにおいて、不規則に分布する。可変k空間密度は、k空間20における動きによって破損したk空間サンプルの規則的分布を回避することによって、ゴースト発生を減少させる。この可変密度収集は、画像コントラストが保たれ、基本的に、周辺k空間部22、23における収集の数の選択とは独立しているTSE又はTFEといったターボシーケンスにおいて特に有利である。図2に示される実施形態において見られるように、中心k空間部は、より頻繁に更新される一方で、周辺k空間部22、23は、低い頻度で更新される。したがって、中心k空間部21からの最も関連性のある信号支配的なk空間ラインが、短期間(例えばマルチエコーイメージングシーケンスの各ショット内)内に更新される。更に、k空間部21、22、23のサイズは、外側k空間まで増加し、これは、k空間20に亘る、動きによって生じたk空間の不連続性の更なる分布をもたらす。
図4は、それぞれ、ACQ1及びACQ2と表される第1及び第2の収集ステップ中に収集されたk空間ラインを示す。この実施形態では、周辺k空間部22におけるk空間ラインの位置が、第2の収集ステップACQ2において、k方向にシフトされている。第1の収集ステップACQ1におけるk空間ラインのk位置は、説明のために、図4の右図において、点線で示されている。対照的に、中心k空間21では、第1の収集ステップACQ1及び第2の収集ステップACQ2中に走査されるk空間ラインのk位置は、維持され、結果として、中心k空間部21におけるk空間ラインは、対応して繰り返しサンプリングされる一方で、周辺k空間部22における各k空間ラインは、1回しか走査されない。各個別の収集ステップACQ1、ACQ2中、周辺k空間部22は、不完全にサンプリングされる。これは、k方向におけるサンプリング密度が、ナイキスト閾値よりも低いことを意味する。両方の収集ステップACQ1、ACQ2においてサンプリングされたk空間ラインの組み合わせは、完全にサンプリングされたデータセットをもたらす。
本発明の主な特徴は、各繰り返しにおいて異なる順序で、各k空間ラインを繰り返しサンプリングすることであり、サンプリングの順序は、k空間ラインがサンプリングされる時間的順序を規定する。好適には、異なるランダムな順序が、収集ステップ毎に適用される。図2乃至図5に示される本発明のすべての実施形態と組み合わせて適用される本発明の技術は、動きアーチファクトを平均し、収集されたk空間データの、動きによって生じた不連続性を分布させ、結果として、動きによって破損したk空間データの不規則なk空間間隔がもたらされる。このようにすると、再構成されたMR画像におけるゴースト発生が回避される。例えば図4を参照するに、収集ステップACQ2中に、中心k空間部21からk空間ラインがサンプリングされる(好適にはランダムな)時間的順序は、収集ステップACQ1中に、中心k空間部21に適用される時間的サンプリング順序とは異なる。
図5に示される実施形態では、周辺k空間部22に、k空間サンプリンググリッドのランダム分布が適用されている。収集ステップACQ1中に収集し損ねた周辺k空間部22におけるk空間ラインが、収集ステップACQ2中にサンプリングされる。両方の収集ACQ1、ACQ2の組み合わせは、完全にサンプリングされたデータセットをもたらす。図4の実施形態と同様に、中心k空間部21は、繰り返しサンプリングが達成されるように、各収集ステップACQ1、ACQ2中に、完全にサンプリングされる。各収集ステップACQ1、ACQ2中の周辺k空間部22のランダムなアンダーサンプリングは、k空間に亘る動きによって破損したk空間サンプルの更に向上された分布、したがって、更に向上された動きロバスト性がもたらされる。
図6に示される実施形態では、再び、k空間20は、中心k空間部21と、周辺k空間部22とに分割される。中心k空間部21は、図示される実施形態では、シングルエコースピンエコー(SE)シーケンスの、対応して繰り返される適用によって、繰り返しサンプリングされる。中心k空間部21は、K1、K2、K3、K4、K5と表される隣接して配置される有限数のk空間ラインを含む。各k空間ラインK1、K2、K3、K4、K5は、SEシーケンスによって生成されたエコー信号に対応する。k空間ラインK1、K2、K3、K4、K5は、本発明によれば、繰り返し毎に、異なる順序で繰り返しサンプリングされる。順序は、例えば4回の繰り返しの間に、次の通りに変えられる。
1回目の繰り返し:K1、K2、K3、K4、K5
2回目の繰り返し:K2、K1、K4、K3、K5
3回目の繰り返し:K2、K4、K1、K5、K3
4回目の繰り返し:K3、K5、K2、K4、K1
信号平均化は、これに基づいて行われて、動きを平均し、アーチファクト及びT緩和といった信号減衰効果を相殺し、SNRの向上を提供する。周辺k空間部22は、1回しか走査されない。最後に、MR画像が、収集されたエコー信号から再構成される。
Dixon式測定では、動きゴーストは、通常、図7の単純な同位相逆位相Dixon例について示されるように、スペクトル源に依存して、水又は脂肪の画像においてエンコードされる(W=IP+OP、F=IP−OP)。ここでは、水の動きは、脂肪画像ではなく、水画像において、全振幅を有する動きゴーストにつながる。図7の左部分は、相関していない動きに関し、図7の右部分は、相関している動きに関する。
左下の単純な同位相逆位相Dixon例について示されるように、動きが相関されていない場合、動きゴーストは、水画像及び脂肪画像の両方に存在するが、アーチファクトは、2倍減少される。この平均化効果は、動きアーチファクトを効果的に減少させるように使用される。
DIXONのTSEでは、Dixon収集毎に、異なるプロファイル順序が、以下に示されるショットの異なる、例えばランダムにされた順序によって達成される。ショット毎に、プロファイルが、1つのk空間セグメントあたりに収集される(図8の例では、7個のk空間セグメントから7個のプロファイルがサンプリングされる)。ショットにおけるプロファイルの非等距離又はランダムサンプリング(青線)は、Dixon収集について、例えばIP及びOP収集について変えられる。動きがある場合、このDixon収集毎の異なるサンプリングは、IP及びOP収集において相関していない動きゴーストにつながる。Dixon再構成を適用すると、この動きゴーストは平均化される。
本発明の一態様によれば、k空間部は、Dixon収集毎に、異なる順序、好適には、ランダムな順序で、ランダムに順序付けられる。
TSE又はTFEといったターボシーケンスについて、ショットに亘り収集され、したがって、k空間セグメント内のk空間において近くに配置されるプロファイルは、ランダムに分布し、このランダムな順序は、Dixon収集と共に変化する。
FFEといった非ターボシーケンスについて、k空間は、ランダムセグメント順序でk空間セグメントに分割される。ランダムセグメント順序は、Dixon収集毎に変えられる。
平均化効果は、より多くのDixon収集が使用されるほど、高められる。したがって、典型的に6つのエコーが使用されるDixon量子化の場合、平均化効果は、最適化される。

Claims (15)

  1. MRデバイスの検査ボリューム内に配置される物体のMRイメージング方法であって、
    前記物体に、マルチエコーイメージングシーケンスを与えることによって、エコー信号を発生させるステップと、
    前記エコー信号を収集するステップと、
    収集された前記エコー信号からMR画像を再構成するステップと、
    を含み、
    前記物体に、前記マルチエコーイメージングシーケンスの幾つかのショットを与え、各ショットによって、一連のエコー信号が発生し、
    各エコー信号は、k空間ラインに帰し、
    (k,k)空間の一部において隣接して配置される幾つかのk空間ラインが、繰り返しサンプリングされ、前記幾つかのk空間ラインは、繰り返し毎に、異なる順序で、かつ、異なる選択されたスピン条件で、サンプリングされる、方法。
  2. 信号平均化が、前記幾つかのk空間ラインの前記繰り返しサンプリングによって行われ、k空間は、中心k空間部と、1つ以上の周辺k空間部とに分割され、前記中心k空間部における信号平均値の数が、前記1つ以上の周辺k空間部における信号平均値の数よりも高く、前記中心k空間部は、前記周辺k空間部よりも高いk空間密度でサンプリングされる、請求項1に記載の方法。
  3. k空間の少なくとも1つの部分は、可変密度でサンプリングされ、前記少なくとも1つの部分のサンプリング密度分布は、繰り返し間で異なる、請求項1又は2に記載の方法。
  4. 繰り返しサンプリングされるk空間ラインの少なくとも幾つかは、ランダムな順序でサンプリングされる、請求項1乃至3の何れか一項に記載の方法。
  5. 前記マルチエコーイメージングシーケンスは、スピンエコーシーケンスであり、そのうちの1つのショットは、
    磁気共鳴の励起のための少なくとも1つのRFパルスと、
    複数のリフォーカシングRFパルスと、
    各エコー信号の周波数エンコーディングのための複数の切り替え磁場勾配と、
    各エコー信号の位相コーディングのための複数の切り替え磁場勾配と、
    を含む、請求項1乃至4の何れか一項に記載の方法。
  6. 前記マルチエコーイメージングシーケンスは、勾配エコーシーケンスであり、そのうちの1つのショットは、
    磁気共鳴の励起のための少なくとも1つのRFパルスと、
    磁気共鳴をリフォーカスし、各エコー信号を周波数エンコーディングするための複数の切り替え磁場勾配と、
    各エコー信号を位相エンコーディングするための複数の切り替え磁場勾配と、
    を含む、請求項1乃至5の何れか一項に記載の方法。
  7. k空間の少なくとも一部は、不完全にサンプリングされる、請求項1乃至6の何れか一項に記載の方法。
  8. 動き情報は、前記繰り返しサンプリングされたkラインから収集される前記エコー信号から導出される、請求項1乃至7の何れか一項に記載の方法。
  9. 前記MR画像は、圧縮センシングを使用して、再構成される、請求項1乃至8の何れか一項に記載の方法。
  10. 前記マルチエコーイメージングシーケンスは、様々なエコー時間が各繰り返しにおいて選択されるDixonシーケンスであり、前記収集されたエコー信号への水及び脂肪からの信号寄与は分離される、請求項1乃至9の何れか一項に記載の方法。
  11. 前記エコー信号は、2つ以上の異なるエコー時間において収集され、前記エコー信号は、各エコー時間において、前記k空間ラインの異なる順序で収集される、請求項10に記載の方法。
  12. 水及び脂肪からの前記信号寄与は、前記MR画像を再構成する前記ステップにおいて組み合わされる、請求項10又は11に記載の方法。
  13. 前記マルチエコーイメージングシーケンスは、反転回復時間後の励起RFパルスが後に続く反転RFパルスを含む反転回復シーケンスである、請求項1乃至12の何れか一項に記載の方法。
  14. 検査ボリューム内に均一な静磁場を発生させる少なくとも1つの主磁石コイルと、前記検査ボリューム内の様々な空間方向において切り替え磁場勾配を発生させる幾つかの勾配コイルと、前記検査ボリューム内にRFパルスを発生させる及び/又は前記検査ボリューム内に配置される物体からのMR信号を受信する少なくとも1つのRFコイルと、RFパルス及び切り替え磁場勾配の時間的連続を制御する制御ユニットと、再構成ユニットとを含み、請求項1乃至12の何れか一項に記載の方法を行う、MRデバイス。
  15. 請求項1乃至12の何れか一項に記載の方法を行うための命令を含む、MRデバイス上で実行されるコンピュータプログラム。
JP2016569024A 2014-05-27 2015-05-22 マルチエコーk空間収集を使用するMRイメージング Active JP6636461B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP14170037 2014-05-27
EP14170037.7 2014-05-27
PCT/EP2015/061441 WO2015181091A2 (en) 2014-05-27 2015-05-22 Mr imaging using multi-echo k-space acquisition

Publications (3)

Publication Number Publication Date
JP2017516551A JP2017516551A (ja) 2017-06-22
JP2017516551A5 JP2017516551A5 (ja) 2018-06-28
JP6636461B2 true JP6636461B2 (ja) 2020-01-29

Family

ID=50846800

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016569024A Active JP6636461B2 (ja) 2014-05-27 2015-05-22 マルチエコーk空間収集を使用するMRイメージング

Country Status (5)

Country Link
US (1) US10466322B2 (ja)
EP (1) EP3149508B1 (ja)
JP (1) JP6636461B2 (ja)
CN (1) CN106461746B (ja)
WO (1) WO2015181091A2 (ja)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2893363B1 (en) * 2012-09-04 2020-08-19 Koninklijke Philips N.V. Propeller with dixon water fat separation
EP3039441A1 (en) * 2013-08-30 2016-07-06 Koninklijke Philips N.V. Dixon magnetic resonance imaging
EP3324844A4 (en) 2015-07-23 2018-08-22 Shanghai United Imaging Healthcare Co., Ltd. System and method for magnetic resonance imaging
CN108603921B (zh) * 2016-02-09 2021-08-10 皇家飞利浦有限公司 稳态磁共振指纹
EP3545325A1 (en) * 2016-11-24 2019-10-02 Koninklijke Philips N.V. Mr imaging with dixon-type water/fat separation
US10481232B2 (en) * 2017-04-13 2019-11-19 Canon Medical Systems Corporation Magnetic resonance imaging apparatus and magnetic resonance imaging method
CN108324276B (zh) * 2018-01-11 2021-07-30 上海东软医疗科技有限公司 磁共振成像方法和装置
CN108535674B (zh) * 2018-02-06 2020-11-20 苏州朗润医疗系统有限公司 用于快速自旋回波的多次平均减轻伪影的方法
EP3686621A1 (en) * 2019-01-23 2020-07-29 Siemens Healthcare GmbH Imaging method with variable density in k-space for multi-spin echo magnetic resonance imaging with different delay times after a preparation pulse and/or different echo times
DE102019201385B3 (de) 2019-02-04 2020-07-02 Siemens Healthcare Gmbh Verfahren zur Aufnahme von Magnetresonanzdaten, Magnetresonanzeinrichtung, Computerprogramm und elektronisch lesbarer Datenträger
CN110082696B (zh) * 2019-04-30 2021-06-08 上海东软医疗科技有限公司 一种多回波序列的成像方法和装置
JP7350568B2 (ja) * 2019-08-28 2023-09-26 キヤノンメディカルシステムズ株式会社 磁気共鳴データ収集方法及び磁気共鳴イメージング装置
CN110794351B (zh) * 2019-11-05 2021-11-16 国家纳米科学中心 一种双回波剪影磁共振成像方法及系统
DE102021203269A1 (de) * 2021-03-31 2022-10-06 Siemens Healthcare Gmbh Verfahren zur Reduzierung von Artefakten in mittels Magnetresonanz aufgenommenen Bilddatensätzen
US20230169704A1 (en) * 2021-11-29 2023-06-01 Canon Medical Systems Corporation Advanced signal combination for improved overlapping image correction
CN114236442B (zh) * 2021-12-14 2022-09-16 无锡鸣石峻致医疗科技有限公司 一种对核磁共振信号进行运动不敏感采集的方法、装置、计算机设备及核磁共振检测系统

Family Cites Families (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3332951B2 (ja) * 1992-05-29 2002-10-07 株式会社日立メディコ 磁気共鳴イメージング装置
JPH06169896A (ja) * 1992-12-10 1994-06-21 Hitachi Medical Corp 磁気共鳴イメージング装置
JPH06169869A (ja) * 1992-12-10 1994-06-21 Matsushita Electric Ind Co Ltd 電気掃除機
US6255820B1 (en) * 1999-03-08 2001-07-03 Picker International, Inc. Variable bandwidth MRI data collection
WO2005114246A1 (en) 2004-05-13 2005-12-01 Invivo Corporation Method of pseudopolar acquisition and reconstruction for dynamic mri
US7330602B2 (en) 2004-09-17 2008-02-12 The Board Of Trustees Of The Leland Stanford Junior University Variable density signal sampling
DE102006042998B4 (de) * 2006-09-13 2008-07-03 Siemens Ag Messsequenz für die dreidimensionale MR-Bildgebung sowie MR-Gerät
US7602179B2 (en) * 2006-11-22 2009-10-13 The General Hospital Corporation System and method for multi-echo bandwidth match imaging
US8155389B2 (en) 2007-10-02 2012-04-10 The University Of Utah Research Foundation Method and system for motion correction in imaging systems
CN101470178B (zh) 2007-12-29 2013-06-05 西门子(中国)有限公司 一种抑制残余运动伪影的方法及装置
US8030923B2 (en) * 2009-03-31 2011-10-04 General Electric Company Method and system to perform phase correction for species separation in magnetic resonance imaging
US8738113B2 (en) 2010-09-30 2014-05-27 University Of Utah Research Foundation Retrospectively correlated turbo spin echo imaging
US9500732B2 (en) 2011-01-28 2016-11-22 Hitachi, Ltd. Magnetic resonance imaging device and angiography using same
CN103502831B (zh) 2011-04-21 2016-08-31 皇家飞利浦有限公司 对运动中目标的磁共振成像
CN104395772B (zh) 2012-03-19 2017-05-10 皇家飞利浦有限公司 具有在对中心和外周k空间区的采样期间的呼吸运动检测的磁共振图像重建方法
US8768034B2 (en) 2012-03-20 2014-07-01 Siemens Medical Solutions Usa, Inc. Motion compensated MR imaging system
WO2014042026A1 (ja) 2012-09-11 2014-03-20 株式会社 東芝 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
DE102013213940B3 (de) * 2013-07-16 2015-01-08 Siemens Aktiengesellschaft Erfassung von MR-Daten mit unterschiedlicher Dichte

Also Published As

Publication number Publication date
WO2015181091A3 (en) 2016-01-21
CN106461746A (zh) 2017-02-22
US20170199258A1 (en) 2017-07-13
EP3149508A2 (en) 2017-04-05
EP3149508B1 (en) 2023-08-02
US10466322B2 (en) 2019-11-05
CN106461746B (zh) 2020-03-03
WO2015181091A2 (en) 2015-12-03
JP2017516551A (ja) 2017-06-22

Similar Documents

Publication Publication Date Title
JP6636461B2 (ja) マルチエコーk空間収集を使用するMRイメージング
US10191133B2 (en) MR imaging using multi-echo segmented k-space acquisition
JP6513398B2 (ja) 事前情報に制約される正則化を使用するmr画像再構成
JP6998218B2 (ja) 動き検出を用いるmr撮像
US9746539B2 (en) MR imaging with suppresion of flow artifacts
US10401456B2 (en) Parallel MR imaging with Nyquist ghost correction for EPI
US10203394B2 (en) Metal resistant MR imaging
US10234523B2 (en) MRI with dixon-type water/fat separation with estimation of the main magnetic field variations
JP2016512780A (ja) 圧縮センシングを使用したmr画像の再構成
CN106796274B (zh) 具有伪迹抑制的propeller-mr成像
JP2019522513A (ja) ディクソン法による水/脂肪分離を用いたmr撮像
JP6074126B1 (ja) k空間中心のサンプリングを用いるゼロエコー時間MR撮像
JP7446328B2 (ja) ソフト動きゲーティングによる3dラジアル又はスパイラル収集を用いたmr画像
JP2019535434A (ja) ディクソンタイプの水/脂肪分離を用いたmrイメージング
JP2020522344A (ja) 並列マルチスライスmr撮像
JP2019535435A (ja) プロペラmrイメージング
JP7128972B2 (ja) 歪み補正を伴うepi mr画像
EP3185029A1 (en) Mr imaging using propeller acquisition with t2 decay correction
EP3736593A1 (en) Dual echo steady state mr imaging using bipolar diffusion gradients
WO2018001759A1 (en) Diffusion weighted mr imaging using multi-shot epi with motion detection and modified sense reconstruction
Schilt Accelerating MRI with optimized sampling patterns in Pulseq
WO2023165865A1 (en) Diffusion mapping by mr fingerprinting

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180517

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180517

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20190128

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190305

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20190522

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190903

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: 20191119

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20191218

R150 Certificate of patent or registration of utility model

Ref document number: 6636461

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250