JP2019503750A - 動きに起因する拡散勾配不整合の補正を備えるdti - Google Patents

動きに起因する拡散勾配不整合の補正を備えるdti Download PDF

Info

Publication number
JP2019503750A
JP2019503750A JP2018532370A JP2018532370A JP2019503750A JP 2019503750 A JP2019503750 A JP 2019503750A JP 2018532370 A JP2018532370 A JP 2018532370A JP 2018532370 A JP2018532370 A JP 2018532370A JP 2019503750 A JP2019503750 A JP 2019503750A
Authority
JP
Japan
Prior art keywords
echo
diffusion
gradient
motion
shot
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.)
Granted
Application number
JP2018532370A
Other languages
English (en)
Other versions
JP2019503750A5 (ja
JP6873137B2 (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 JP2019503750A publication Critical patent/JP2019503750A/ja
Publication of JP2019503750A5 publication Critical patent/JP2019503750A5/ja
Application granted granted Critical
Publication of JP6873137B2 publication Critical patent/JP6873137B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/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/56341Diffusion imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • 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/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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion

Landscapes

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

Abstract

本発明は、MRデバイス1の検査ボリューム内に配置される対象物10のMR撮像の方法に関する。本発明の目的は、動きに起因する拡散符号化の変化によって生じるアーチファクトを最小限にした拡散テンソル撮像(DTI)を可能にすることである。本発明の方法は、前記対象物10に、マルチエコー撮像シーケンスの少なくとも1つのショットを受けさせるステップであって、前記ショットによってエコー信号の列が生成され、前記マルチエコー撮像シーケンスは、拡散勾配を含む、ステップと、様々な拡散勾配のエコー信号を取得するステップと、前記エコー信号の取得中に前記対象物10の動きを検出するステップと、検出された前記動きに基づいて、前記対象物10の複数の動きの状態を決定するステップであって、前記対象物10は、検出された動きの各状態においては静止しているステップと、検出された前記対象物10の前記動きの状態についての前記対象物10の基準フレームに対する拡散勾配の方向をそれぞれ決定するステップと、決定された前記拡散勾配の相対的方向について取得されるエコー信号から拡散係数を導出するステップと、を含む。さらには、本発明は、この方法を実行するためのMRデバイス並びにMRデバイス上で実行されるコンピュータプログラムに関する。

Description

本発明は、磁気共鳴(MR)撮像の分野に関する。本発明は、対象物のMR撮像の方法に関する。本発明はまた、MRデバイス及びMRデバイス上で実行されるコンピュータプログラムに関する。
2次元又は3次元画像を形成するために磁場と核スピンとの間の相互作用を利用する画像形成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画像に変換される。
MR撮像は、拡散への感度が高い。よく知られる拡散強調撮像(DWI)技術は、通常、拡散勾配を含む撮像シーケンスを使用することによって実行され、拡散勾配の方向に沿った(水分子の)プロトンの拡散は、取得されたMR信号の振幅を減少させる。拡散テンソル撮像(DTI)は、DWIのより洗練された形式であり、拡散の大きさと方向との両方を決定することを可能にする。例えば、DTIは、MR脳撮像において白質線維の視覚化を可能にし、脳梗塞、多発性硬化症、てんかんなどの疾患に関連付けられる白質の微細な変化をマッピングすることができる。いわゆる異方性度(FA)は、MR画像の各ボクセル位置での拡散テンソルの形についての情報を提供する。FAは、拡散テンソルの固有値の分散から決定される。したがって、FAは、所与の画像位置における等方性拡散と線形拡散との間の差異を反映する。DTIの変形として、昨今、拡散テンソルトラクトグラフィ(DTT)と呼ばれる技術が開発されてきた。この技術は、脳内の神経線維の非侵襲的な追跡を可能にする。白質繊維の軌道は、繊維の長手軸に対応すると仮定される最も速い拡散の方向を追跡することによって再構成される。
脳のDWI技術は、動きから生じる信号減衰が関心対象の測定を混乱させる可能性があるため、巨視的な頭部の動きに対して特に脆弱である。MR検査中の被験者の動きは、子供、高齢者、又はパーキンソン病のような静止状態になることを妨げる病状を有する患者などの人々において、特に問題となる可能性がある。動きは、撮像される脳組織をシフトする(再構成されるMR画像にゴーストアーチファクトをもたらす)、及び不正確な拡散符号化にさらすという主に2つの点でデータに影響を与える。
拡散テンソルを決定する前の遡及的な動き補正方法が広く使用されている。動きの時間を遡及的に補正する基本的かつ一般的なやり方は、拡散強調MR画像と基準MR画像(強調されていない)との相互位置合わせと、各画像位置での動きを考慮に入れた拡散勾配方向の連続的な再配向とを利用する。このような工程は空間補間を含み、補間は部分容積効果及び拡散テンソル計算における伝播を伴うDWIの分散特性に影響を及ぼし得る。
動きに起因する重大なアーチファクトを回避するために、DWIデータは、一般的に、シングルショットエコープラナー撮像(EPI)などのシングルショット撮像シーケンスを使用して取得される。しかしながら、シングルショットDWIでは画質が低く、空間分解能が制限される可能性がある。重大な幾何学的歪み及び制限された空間分解能は、拡散特性を高精度で測定することを困難にする。昨今、シングルショットDWIの限界に対処するための試みがなされている。
米国特許出願公開第2014/0002078A1号には、並列取得を使用し、動きによるショット間位相変動を固有の方式で補正し、これによりゴーストアーチファクトを回避するマルチショットDWI技術(MUSE−多重化センシティビティ符号化と呼ばれる)が記載されている。
Joeng等による「Magnetic Resonance in Medicine」volume 69(3), pages 793-802, 2013は、高速並行画像取得に一般的に使用される標準SENSEアルゴリズムの調整を用いたマルチショットDWI技術を提案している。この調整は、ショット間の動きに起因する位相誤差に対処する。この知られている技術は、画像空間サンプリング機能を用いた画像再構成(IRIS)と呼ばれる。
しかしながら、これらの技術は、動きに起因する誤った拡散符号化を考慮していない。
MURAT AKSOYによる「Effects of motion and b-matrix correction for high resolution DTI with short-axis PROPELLER-EPI」では、スキャナの基準フレームから患者の基準フレームへの回転行列、並進行列R及びΔrを含む方程式を解くことによって、複雑なk空間データから直接推定される拡散テンソルが提示される。方程式は、非線形共役勾配アルゴリズムを使用して解かれ、この詳細は、「Single-step nonlinear diffusion tensor estimation in the presence of microscopic and macroscopic motion」に記載されている。
前述の記載から、改善されたDWI技術が必要であることは容易に理解される。したがって、本発明の目的は、拡散符号化の動きに起因する変化によって引き起こされるアーチファクトを最小限にしたDWI(及びDTI/DTT)を可能にすることである。
本発明によれば、MRデバイスの検査ボリューム内に配置される対象物のMR撮像の方法が開示される。この方法は、
対象物に、拡散勾配を含むマルチエコー撮像シーケンスの少なくとも1つのショットを受けさせるステップであって、前記ショットによってエコー信号の列が生成されるステップと、
様々な拡散勾配のエコー信号を取得するステップと、
エコー信号の取得中に対象物(10)の動きを検出するステップと、
検出された動きに基づいて、対象物(10)の複数の動きの状態を決定するステップであって、対象物(10)は、検出された動きの各状態においては静止しているステップと、
検出された対象物(10)の動きの状態についての対象物(10)の基準フレームに対する拡散勾配の相対的方向をそれぞれ決定するステップと、
決定された拡散勾配の相対的方向について取得されるエコー信号から拡散係数を導出するステップと、を含む。
言い換えれば、本発明は、拡散勾配の相対的方向における動きに起因する変化を付加的な拡散勾配方向として扱い、これらを拡散係数の計算に直接使用することを提案する。
本発明における基準フレームの意味は、例えば、患者の頭部のような検査される対象物に対して固定された空間的関係にある座標系である。
本発明の技術では、対象物がある特定の動きの状態にある間、対象物の動きがなく、したがって対象物の基準フレームにおける拡散勾配方向の変動が生じないと仮定する。動きが検出されると、撮像された対象物の位置及び向きが変化する。したがって、検出された動きは、本発明の意味の範囲内で、一つの動きの状態を他方の動きの状態から時間的に分離する。取得されたエコー信号は、様々な動きの状態に属性付けられる。最終的に拡散テンソルが対象物の基準フレーム内で取得されるように、同一の動きの状態に属性付けられるエコー信号及び決定された拡散勾配方向から、各拡散係数が本発明に従って計算される。
このアプローチでは、拡散テンソルの要素を計算する前に、拡散符号化の変動を補正する必要はない。
好ましくは、本発明に従って取得される拡散係数から画像が生成される。画像は、臨床研究で一般的に使用される、異方性度(FA)マップ、平均拡散率(MD)マップ、放射状拡散率(RD)マップ、若しくは軸方向拡散率(AD)マップ、又は拡散係数から導出される他の任意のスカラー計量のマップである。
本発明によるマルチエコー撮像シーケンスの1つの「ショット」は、少なくとも1つの(典型的には180°の)再収束RFパルスが後に続く磁気共鳴の励起のための初期RFパルスを含み、拡散勾配は、再収束RFパルスの前後に印加される。このパルスのシーケンスは、別々に位相符号化及び周波数符号化された勾配リコールのエコー信号の列として測定される拡散符号化スピンエコーを生成する。これらのエコー信号が取得され、各エコー信号は、k空間プロファイルを表す。取得された信号データから完全なMR画像を再構成できるように、マルチエコーシーケンスの1つの単一ショット又は複数のショットが、k空間の完全なサンプリングのために適用され得る。高画質を達成し、高精度で拡散特性を測定するために、マルチショットマルチエコー撮像シーケンスが好適であり得る。
本発明の一実施形態では、動きは、マルチエコー撮像シーケンスの異なるショット間でのみ検出され得る。この場合、各ショットが対象物の1つの動きの状態に属性付けられ得る(複数のショットが同じ動きの状態に属性付けられ得る)ように、シーケンスの各ショットの間には動きが生じないと仮定される。k空間がシーケンスの単一のショット又はショットのグループの中で完全にサンプリングされる場合、このショット内/ショットのグループ内で取得されたエコーデータから完全な画像を再構成し、このショットに/ショットのグループに属性付けられる相対的勾配方向に対する拡散計数を計算することが可能である。次いで、拡散テンソルの残りのテンソル要素は、撮像シーケンスの他のショットに属性付けられるエコー信号データ及び相対的勾配方向から計算することができる。
同様に、信号対雑音比(SNR)を改善するために信号平均化を使用して、エコー信号を取得することができる。この場合、異なる信号平均化ステップ間の動きが検出され、各信号平均化ステップは、対象物の1つの動きの状態に属性付けられる。1つの信号平均化ステップで取得されるエコーデータからMR画像が再構成され、この平均化ステップに属性付けられる相対的勾配方向について拡散係数が計算される。次いで、拡散テンソルの残りのテンソル要素は、撮像シーケンスの他の信号平均化ステップに属性付けられるエコー信号データ及び相対的勾配方向から計算することができる。
DTIでは、2つ以上の異なる拡散勾配方向のセットの信号データを取得し、異なる拡散勾配方向のそれぞれについての拡散係数を計算することが必要である。勾配方向は、いくつかの連続したステップの中で変化する。拡散勾配方向の変化の間で動きが検出され、次いで、各拡散勾配方向が対象物の1つの動きの状態に属性付けられ得る。拡散勾配方向ごとに取得されたエコーデータからMR画像が再構成され、対象物の基準フレーム内の対応する勾配方向について拡散係数が計算される。次いで、拡散テンソルの残りのテンソル要素が、他の拡散勾配変動ステップに属性付けられるエコー信号データ及び相対勾配方向から計算され得る。
1つの可能な実施形態において、本発明の方法において使用される撮像シーケンスは、マルチショットEPIシーケンスである。
動きは、例えば、既知のナビゲータ技術によって検出されてもよい。この目的のために、ナビゲータエコーの生成が、本発明によって使用される撮像シーケンスに組み込まれ得る。あるいは、動きは、k空間整合性基準の使用により、取得されたエコー信号データの一部から再構成された低解像度MR画像のレジストレーションにより、又は内部及び/若しくは外部の動きセンサにより検出されてもよい。動きが検出された場合、対象物の基準フレームに対する拡散勾配の変更された方向が決定され、対象物の新しい動きの状態に属性付けられる。
これまでに説明した本発明の方法は、検査ボリューム内に均一な静磁場を発生させるための少なくとも1つの主磁石コイルと、検査ボリューム内の異なる空間方向において切替磁場勾配を発生させるためのいくつかの勾配コイルと、検査ボリューム内にRFパルスを生成するため、及び/又は検査ボリューム内に配置される対象物からMR信号を受信するための少なくとも1つのRFコイルと、RFパルス及び切替磁場勾配の時間的連続を制御する制御ユニットと、再構成ユニットと、を備えるMRデバイスによって実行され得る。本発明の方法は、例えば、再構成ユニット及び/又はMRデバイスの制御ユニットに対応するプログラミングによって実施することができる。
本発明の方法は、現在の臨床用途のほとんどのMRデバイスにおいて有利に実行することができる。この目的のためには、本発明の上述の方法ステップを実行するようにMRデバイスを制御するコンピュータプログラムを利用するだけでよい。コンピュータプログラムは、MRデバイスの制御ユニットにインストールされるためにダウンロードされるように、データキャリア上に存在しても、データネットワーク内に存在してもよい。
添付の図面は、本発明の好ましい実施形態を開示する。しかしながら、図面は説明の目的のためにのみデザインされており、本発明の範囲の限定の定義ではないことを理解されたい。
図1は、本発明の方法を実行するためのMRデバイスを示す。 図2は、本発明の実施形態で使用される撮像シーケンスの図を示す。 図3は、本発明の勾配不整合補償を用いて及び用いないで取得されるFAマップを示す。
図1を参照すると、MRデバイス1が示されている。このデバイスは、検査ボリューム全体に、z軸に沿った、実質的に均一で、時間的に一定である主磁場が形成されるように、超伝導又は抵抗性の主磁石コイル2を含む。
磁気共鳴生成及び操作システムは、MR撮像を実行するために、一連のRFパルス及び切替磁場勾配を印加して、核磁気スピンを反転又は励起し、磁気共鳴を誘発し、磁気共鳴を再収束し、磁気共鳴を操作し、磁気共鳴を空間的に又はそれ以外のやり方で符号化し、スピンを飽和させる。
より具体的には、勾配パルス増幅器3は、検査ボリュームのx軸、y軸及びz軸に沿った全身勾配コイル4、5、及び6のうちから選択される1つに電流パルスを印加する。デジタル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を制御して、拡散強調エコープラナー撮像(DWI−EPI)などのような複数のMR撮像シーケンスのいずれかを生成する。選択されたシーケンスについて、受信器14は、各RF励起パルスに続いて、単一又は複数のMRデータラインを続々と受信する。データ収集システム16は、受信信号のアナログ−デジタル変換を行い、各MRデータラインをさらなる処理に適したデジタルフォーマットに変換する。最近のMRデバイスでは、データ収集システム16は、生画像データの取得に特化した独立したコンピュータである。
最終的に、デジタル生画像データは、フーリエ変換又はSENSE、SMASH、若しくはGRAPPAのような他の適切な再構成アルゴリズムを適用する再構成プロセッサ17によって画像表現に再構成される。MR画像は、患者を通る平面スライス、平行な平面スライスのアレイ、又は三次元ボリュームなどを表すことができる。次いで、画像は、画像メモリの中に記憶され、そこで、画像表現のスライス、投影、又は他の部分を視覚化に適切な形式に変換するために処理され得、これは、例えば、結果としてできたMR画像のための、人間が解読可能なディスプレイを提供するビデオモニタ18を介して行われる。
引き続き図1を、さらに図2及び図3も参照して、本発明の方法の実施形態を下記に説明する。
図2に示されるように、身体10は、マルチエコー撮像シーケンスの複数のショットを受ける。撮像シーケンスは、スピンエコー拡散強調マルチショットEPIシーケンスである。第2の180°再収束パルスの後に、ナビゲータエコー取得NAVが続く。実際の拡散撮像及びナビゲータエコー取得NAVに使用されるエコー信号取得IMGは、SENSEによって加速されてもよい。異なる拡散勾配方向は、撮像シーケンスの複数のショットを使用して適用される。図2の点線は、最初の180°RF再収束パルスの前後に適用された拡散勾配を表している。エコー信号は、様々な拡散勾配に対して取得される。
対象物の動きは、ナビゲータエコーを使用して検出される。検出された動きから、身体10に対する拡散勾配方向が決定され、これは、拡散勾配の方向は、身体10の各検出された動きの状態についての身体10の基準フレームに対して決定されることを意味する。このように決定された相対的勾配方向のセットは、身体10の各動きの状態に属性付けられる。
例えば、エコー信号は、図2に示される撮像シーケンスを使用して拡散勾配の3つの異なる方向について、取得される。完全なk空間エコー信号データセットが、撮像シーケンスの4つのショットの中で取得される。3つの信号平均化ステップが使用される(平均信号数NSA=3)。これは、3×4×3=36ショットが実行されることを意味する。12ショットの各グループは、同じ拡散勾配方向を有する。動きが検出され、例えば、身体10の5つの異なる動きの状態が識別されることになる。身体10の基準フレームに対する拡散勾配の方向は、ナビゲータエコーからの各動きの状態ごとに決定される。これら5つの相対的勾配方向は、この例では、撮像シーケンスの8、10、4、6、8ショットのグループにそれぞれ属性付けられる。これらの5つのグループのそれぞれの取得中に動きは生じない。上述のIRIS又はMUSE技術は、ショットの各グループの画像の再構成のために使用され、これによって動きに起因する位相不整合を補正することができる。次に、これらの5つの画像は、画像空間に位置合わせされ、5つの相対的勾配方向の拡散係数が導出される。
図3は、本発明の拡散勾配不整合補償を用いて、及び用いずに得られたFAマップを示す。図3aは、いかなる動きもないFAマップを示す。図3b及び図3cでは、撮像される頭部に取得中40°の回転が発生した。図3bは、位相及び大きさの不整合補正(MUSE又はIRIS方式を使用して)のみで得られたFAマップを示す。図3cは、位相及び大きさの不整合補正(MUSE/IRIS)、並びに本発明による勾配方向の考慮と共に得られるFAマップを示す。図3bの白い矢印は、動きに起因する不整合の拡散符号化によるアーチファクトの位置を示している。このアーチファクトは、図3cでは見られない。
本発明の効果をさらに説明するために、図3の画像データにおいて、脳梁の膝節及び膨大部(GCC及びSCC)のそれぞれで2つの関心領域が選択される。以下の表は、本発明の方法を使用しないで(図3a及び図3bに基づくデータを比較することによって)決定された、及び本発明の方法を使用して(図3a及び図3cに基づくデータを比較することによって)決定された拡散テンソルの主要な固有ベクトルの角度偏差(AD)が度単位で示され、また、同様の態様で、FAの平均平方誤差(RMSE)が示される。
Figure 2019503750

Claims (12)

  1. MRデバイスの検査ボリューム内に配置される対象物のMR撮像の方法であって、前記方法は、
    前記対象物に、拡散勾配を含むマルチエコー撮像シーケンスの少なくとも1つのショットを受けさせるステップであって、前記ショットによってエコー信号の列が生成される、ステップと、
    様々な印加の方向を有する拡散勾配のエコー信号を取得するステップと、
    前記エコー信号の取得中に前記対象物の動きを検出するステップと、
    検出された前記動きに基づいて、前記対象物の複数の動きの状態を決定するステップであって、前記対象物は、検出された前記動きの各状態においては静止している、ステップと、
    検出された前記対象物の前記動きの状態についての当該対象物の基準フレームに対する拡散勾配の相対的方向をそれぞれ決定するステップと、
    決定された前記拡散勾配の相対的方向について取得されるエコー信号から拡散係数を導出するステップと、を含む、方法。
  2. 前記対象物に、前記マルチエコー撮像シーケンスのいくつかのショットを受けさせる、請求項1に記載の方法。
  3. 前記マルチエコー撮像シーケンスは、ナビゲータエコーを含む、請求項1又は2に記載の方法。
  4. 前記動きは、前記ナビゲータエコーを使用して検出される、請求項3に記載の方法。
  5. 前記ナビゲータエコーからの各動きの状態について、前記対象物の前記基準フレームに対する拡散勾配の相対的方向が決定される、請求項3又は4に記載の方法。
  6. 前記エコー信号は、信号平均化ステップを使用して取得される、請求項1乃至5の何れか一項に記載の方法。
  7. 様々な前記信号平均化ステップの間の動きが検出され、前記信号平均化ステップの各々は、前記対象物の1つの動きの状態に属性付けられる、請求項6に記載の方法。
  8. 前記マルチエコー撮像シーケンスは、マルチショットEPIシーケンスである、請求項1乃至7の何れか一項に記載の方法。
  9. 取得された前記エコー信号を様々な動きの状態に属性付けるステップと、
    前記エコー信号と、同じ動きの状態に属性付けられる決定された前記拡散勾配の相対的方向の各々とによって各拡散係数を計算するステップと、をさらに含む、請求項1乃至8の何れか一項に記載の方法。
  10. 動きに起因する位相不整合について、再構成された画像を補正するステップをさらに含む、請求項1乃至9の何れか一項に記載の方法。
  11. 検査ボリューム内に均一な静磁場を発生させるための少なくとも1つの主磁石コイルと、前記検査ボリューム内の様々な空間方向において切替磁場勾配を発生させるためのいくつかの勾配コイルと、前記検査ボリューム内にRFパルスを生成するため、及び/又は前記検査ボリューム内に配置される対象物からMR信号を受信するための少なくとも1つのRFコイルと、前記RFパルス及び切替磁場勾配の時間的連続を制御する制御ユニットと、再構成ユニットとを備える、MRデバイスであって、前記MRデバイスは、
    前記対象物に、拡散勾配を含むマルチエコー撮像シーケンスの少なくとも1つのショットを受けさせることであって、前記ショットによってエコー信号の列が生成されることと、
    様々な印加の方向を有する拡散勾配に対するエコー信号を取得することと、
    前記エコー信号の取得中に前記対象物の動きを検出することと、
    検出された前記動きに基づいて、前記対象物の複数の動きの状態を決定することであって、前記対象物は、検出された前記動きの各状態においては静止していることと、
    検出された前記対象物の前記動きの状態についての前記対象物の基準フレームに対する拡散勾配の相対的方向をそれぞれ決定することと、
    決定された前記拡散勾配の相対的方向について取得されるエコー信号から拡散係数を導出することとを実行する、MRデバイス。
  12. MRデバイス上で実行されるコンピュータプログラムであって、
    拡散勾配を含むマルチエコー撮像シーケンスの少なくとも1つのショットを生成するステップであって、前記ショットによってエコー信号の列が生成されるステップ、
    様々な印加の方向を有する拡散勾配のエコー信号を取得するステップ、
    前記エコー信号の取得中に前記対象物の動きを検出するステップ、
    検出された前記動きに基づいて、前記対象物の複数の動きの状態を決定するステップであって、前記対象物は、検出された前記動きの各状態においては静止しているステップ、
    検出された前記対象物の前記動きの状態についての前記対象物の基準フレームに対する拡散勾配の相対的方向をそれぞれ決定するステップ、及び
    決定された前記拡散勾配の相対的方向について取得されるエコー信号から、拡散係数を導出するステップ、
    のための命令を含む、コンピュータプログラム。
JP2018532370A 2015-12-22 2016-12-15 動きに起因する拡散勾配不整合の補正を備えるdti Active JP6873137B2 (ja)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
CNPCT/CN2015/098240 2015-12-22
CN2015098240 2015-12-22
EP16166963 2016-04-26
EP16166963.5 2016-04-26
PCT/EP2016/081088 WO2017108543A1 (en) 2015-12-22 2016-12-15 Dti with correction of motion-induced diffusion gradient inconsistency

Publications (3)

Publication Number Publication Date
JP2019503750A true JP2019503750A (ja) 2019-02-14
JP2019503750A5 JP2019503750A5 (ja) 2020-07-27
JP6873137B2 JP6873137B2 (ja) 2021-05-19

Family

ID=57543053

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018532370A Active JP6873137B2 (ja) 2015-12-22 2016-12-15 動きに起因する拡散勾配不整合の補正を備えるdti

Country Status (5)

Country Link
US (1) US10746832B2 (ja)
EP (1) EP3394628A1 (ja)
JP (1) JP6873137B2 (ja)
CN (1) CN108431625B (ja)
WO (1) WO2017108543A1 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3633401A1 (en) * 2018-10-04 2020-04-08 Siemens Healthcare GmbH Prevention of compensating a wrongly detected motion in mri
CN109752681B (zh) * 2019-01-22 2021-10-29 奥泰医疗系统有限责任公司 倾斜平面回波成像方法及核磁共振成像系统
EP3699624A1 (en) * 2019-02-25 2020-08-26 Koninklijke Philips N.V. Calculation of a b0 image using multiple diffusion weighted mr images

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5833609A (en) 1996-11-26 1998-11-10 Picker International, Inc. Rotating diffusion MR imaging reduced motion artifacts
DE102004017852B4 (de) * 2004-04-13 2008-11-27 Siemens Ag Bewegungskorrigiertes Multi-Shot-Verfahren zur diffusionsgewichteten Bildgebung in der Magnetresonanztomographie
US7902825B2 (en) * 2008-05-19 2011-03-08 The Board Of Trustees Of The Leland Stanford Junior University Motion corrected tensor magnetic resonance imaging
DE102011006851B4 (de) * 2011-04-06 2012-12-06 Siemens Aktiengesellschaft Gruppierung von diffusionsgewichteten MR-Bildern nach den verwendeten Diffusionsgradienten zur Erstellung einer Diffusionsinformation
US9687172B2 (en) * 2012-01-24 2017-06-27 National Institute of Health (NIH), The United States of America, U.S. Dept. of Health and Human Services (DHHS) System for motion corrected MR diffusion imaging
SE537065C2 (sv) * 2012-05-04 2014-12-23 Cr Dev Ab Pulssekvensförfarande för MRI
EP2858559B1 (en) * 2012-06-28 2021-01-20 Duke University Multi-shot scan protocols for high-resolution mri incorporating multiplexed sensitivity-encoding (muse)
JP6371776B2 (ja) 2012-12-12 2018-08-08 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 磁気共鳴の拡散強調画像(dwi)のための動作検出および修正方法

Also Published As

Publication number Publication date
US10746832B2 (en) 2020-08-18
CN108431625B (zh) 2021-08-24
US20180372829A1 (en) 2018-12-27
WO2017108543A1 (en) 2017-06-29
CN108431625A (zh) 2018-08-21
EP3394628A1 (en) 2018-10-31
JP6873137B2 (ja) 2021-05-19

Similar Documents

Publication Publication Date Title
US10444315B2 (en) MRI with motion correction using navigators acquired using a dixon technique
JP6998218B2 (ja) 動き検出を用いるmr撮像
JP6243522B2 (ja) 正則化された検出再構成を使用するマルチエコーディクソン水−脂肪分離及びb0歪補正による並列mri
CN105556326B (zh) 具有dixon类型的水/脂肪分离的mr成像
CN111758041A (zh) 使用多梯度回波序列进行的Dixon MR成像
CN110720047B (zh) 双回波Dixon型水/脂肪分离MR成像
JP6873137B2 (ja) 動きに起因する拡散勾配不整合の補正を備えるdti
CN109983358B (zh) Propeller mr成像
CN109716155B (zh) 具有迪克逊型水/脂肪分离的mr成像
US11815582B2 (en) Dual echo steady state MR imaging using bipolar diffusion gradients
JP7128972B2 (ja) 歪み補正を伴うepi mr画像
O'Brien et al. 3-d residual eddy current field characterisation: applied to diffusion weighted magnetic resonance imaging
CN111164444B (zh) 具有经改进的脂肪位移校正的Dixon型水/脂肪分离MR成像
EP3480617A1 (en) Diffusion weighted turbo spin echo mr imaging with motion compensation
WO2018001759A1 (en) Diffusion weighted mr imaging using multi-shot epi with motion detection and modified sense reconstruction
WO2024061852A1 (en) Mr imaging with water/fat/b0 mapping
CN114402214A (zh) 狄克逊型水/脂肪分离mr成像

Legal Events

Date Code Title Description
A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20191212

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20191212

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200612

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20200612

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20200821

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200901

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20201130

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210209

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210420

R150 Certificate of patent or registration of utility model

Ref document number: 6873137

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150