JP2022521277A - ウェーブ符号化を用いたパラレルｍｒイメージング - Google Patents
ウェーブ符号化を用いたパラレルｍｒイメージング Download PDFInfo
- Publication number
- JP2022521277A JP2022521277A JP2021549205A JP2021549205A JP2022521277A JP 2022521277 A JP2022521277 A JP 2022521277A JP 2021549205 A JP2021549205 A JP 2021549205A JP 2021549205 A JP2021549205 A JP 2021549205A JP 2022521277 A JP2022521277 A JP 2022521277A
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- field gradient
- signal
- imaging
- signal profile
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 32
- 238000007689 inspection Methods 0.000 claims abstract description 10
- 238000004590 computer program Methods 0.000 claims abstract description 6
- 230000035945 sensitivity Effects 0.000 claims description 8
- 238000000264 spin echo pulse sequence Methods 0.000 claims description 5
- 230000007423 decrease Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 2
- 238000005070 sampling Methods 0.000 abstract description 6
- 230000005415 magnetization Effects 0.000 description 9
- 238000013459 approach Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 210000004556 brain Anatomy 0.000 description 3
- 229920013655 poly(bisphenol-A sulfone) Polymers 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000011503 in vivo imaging Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
- G01R33/4833—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
- G01R33/4835—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices of multiple slices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image 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/5611—Parallel 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56563—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image 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/5615—Echo 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/5616—Echo 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image 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/5615—Echo 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/5617—Echo 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
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (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)
- Spectroscopy & Molecular Physics (AREA)
- Optics & Photonics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本発明は、ＭＲ装置１の検査ボリューム内に配置される対象１０のＭＲ撮像方法に関する。この方法は、対象を撮像シーケンスに曝すことによりＭＲ信号を生成するステップと、デカルトｋ空間サンプリングスキームでＭＲ信号プロファイルを取得するステップであって、各ＭＲ信号プロファイルが、読み出し方向に沿った時間的に一定の磁場勾配及び位相符号化方向に沿った正弦波変調される磁場勾配の存在下で取得される、ステップと、磁場勾配の変調スキームを考慮して、取得されたＭＲ信号プロファイルからＭＲ画像を再構成するステップと、を有する。本発明は、磁場勾配の正弦波変調の周波数が各ＭＲ信号プロファイルの取得中に変化されることを提案する。更に、本発明は、この方法を実行するＭＲ装置、並びにＭＲ装置上で実行されるコンピュータプログラムに関する。
Description
本発明は、磁気共鳴（ＭＲ）イメージングの分野に関する。本発明は、対象のＭＲイメージング方法に関する。本発明はまた、ＭＲ装置、及びＭＲ装置上で実行されるコンピュータプログラムに関する。
２次元又は３次元画像を形成するために磁場と原子核スピンとの間の相互作用を利用する画像形成ＭＲ方法は、軟組織の撮像に関して、多くの点で他のイメージング方法よりも優れており、電離放射線を必要とせず、通常侵襲性ではないので、今日、特に医療診断の分野で広く使用されている。一般的なＭＲ方法によれば、対象、例えば検査されるべき患者の身体は、強力で均一な磁場内に配置され、かかる磁場の方向は同時に、測定値の基礎となる座標系の軸（通常はｚ軸）を規定する。磁場は、定義された周波数（いわゆるラーモア周波数又はＭＲ周波数）の交番電磁場（ＲＦ場）の印加によって励起され得る磁場強度（スピン共鳴）に依存して、個々の核スピンについて異なるエネルギーレベルを生成する。巨視的観点から、個々の核スピンの分布は、適当な周波数の電磁パルス（ＲＦパルス）の印加により平衡状態から逸らされて偏向されることができる全体的な磁化を生成し、その結果、磁化は、ｚ軸を中心に歳差運動を行う。歳差運動は、開口角をフリップ角と呼ぶ円錐の表面を描写する。フリップ角の大きさは、印加される電磁パルスの強度と持続時間に依存する。いわゆる９０°パルスの場合、スピンは、ｚ軸から横断方向の平面（フリップ角９０°）に偏向される。
ＲＦパルスの終了後、磁化は元の平衡状態に戻り、この場合、ｚ方向の磁化は、第１の時定数Ｔ１（スピン格子又は縦緩和時間）で再び構築され、ｚ方向に直交する方向の磁化は、第２の時定数Ｔ２（スピン－スピン又は横緩和時間）で緩和される。磁化の変動は、磁化の変動がｚ軸に垂直な方向において測定されるように、ＭＲ装置の検査ボリューム内に配置され方向付けられた受信ＲＦコイルことによって検出されることができる。横方向磁化の減衰は例えば、９０°パルスの印加の後、同じ位相を有する秩序状態から全ての位相角が均一に分布する状態（ディフェージング）への原子核スピンの遷移（局所的な磁場の不均一性によって引き起こされる）を伴う。デフェージングは、リフォーカシングパルス（例えば、１８０°パルス）によって補償されることができる。これにより、受信コイルにエコー信号（スピンエコー）を生成する。
体内の空間分解能を実現するために、３つの主軸に沿って伸びる磁場勾配が均一な磁場に重畳され、スピン共鳴周波数の空間依存性をもたらす。次いで、受信コイル内においてピックアップされる信号は、体内の異なるロケーションに関連付けられることができる異なる周波数成分を含む。受信コイルを介して得られる信号データは、空間周波数ドメインに対応し、ｋ空間データと呼ばれる。ｋ空間データは、通常、異なる位相符号化で取得された複数のライン（プロファイルとも呼ばれる）を含む。各ラインは、多数のサンプルを収集することによってデジタル化される。ｋ空間データの組は、画像再構成アルゴリズムによってＭＲ画像に変換される。
近年、パラレル取得と呼ばれるＭＲ取得を加速する技術が開発されている。このカテゴリの方法は、SENSE（Sensitivity Encoding)、SMASH(Simultaneous Acquisition of Spatial Harmonics）、GRAPPA（Generalized Auto-calibrating Partially Parallel Acquisition）である。SENSE、SMASH、及びGRAPPA及び他のパラレル取得技術は、複数のＲＦ受信コイルから並行にサブサンプリングｋ空間データ取得を使用する。これらの方法では、複数のコイルからの（複素）信号データは、最終的に再構成されたＭＲ画像におけるサブサンプリングアーチファクト（エイリアシング）を抑圧するようなやり方で、複雑な重み付けと組み合わせる。このタイプの複雑なアレイ信号の組み合わせは、空間フィルタリングと呼ばれることもあり、ｋ空間ドメイン（SMASH及びGRAPPAにおけるように）又は画像ドメイン（検知におけるように）において実行される組み合わせ、並びにハイブリッドである方法を含む。
Larkman他（Journal of Magnetic Resonance Imaging, 13, 313-317, 2001）は、マルチスライスイメージングの場合、感度符号化をスライス方向にも適用して、スキャン効率を高めることを提案した。Breuer他（Magnetic Resonance in Medicine, 53, 684-691, 2005）は、"controlled aliasing in parallel imaging results in higher acceleration"（ＣＡＩＰＩＲＩＮＨＡ又はＣＡＩＰＩ）というアプローチを提案してこの基本的な考え方を改良した。この技術は、マルチスライス取得中の各個々のスライスにおけるエイリアシングアーチファクトの出現を修正し、後続のパラレル画像再構成プロシージャを改善する。従って、ＣＡＩＰＩは、純粋な後処理アプローチのみを使用する他のマルチスライスパラレルイメージング概念と比較して、より効率的なパラレルマルチスライスイメージング技術である。ＣＡＩＰＩでは、位相変調マルチスライスＲＦパルス（既知のＨａｄａｍａｒｄパルスに類似）の使用と同時に任意の厚さと距離の複数スライスが励起される。取得されたＭＲ信号データはサブサンプリングされ、互いにシフトされたように見える重ね合わせスライス画像が得られる。エイリアシングされたスライス画像のシフトは、フーリエシフト定理に従ってＲＦパルスの位相変調方式によって制御される。位相符号化ステップから位相符号化ステップまで、マルチスライスＲＦパルスは、個々の位相シフトを各スライスのＭＲ信号に適用する。関係するスライスの個々の信号寄与を分離する逆再構成問題の数値条件付けが、このシフトを用いて改善される。ＣＡＩＰＩはまた、使用されるＲＦ受信コイルのコイル感度が、撮像される個々のスライスにおいて劇的に異なることがないように、重ね合わせられたスライス画像の分離を改善する可能性を有する。しかしながら、ＣＡＩＰＩには限界がある。
既知のＣＡＩＰＩ方法を更に改良した、いわゆるウェーブＣＡＩＰＩ（Ｗａｖｅ－ＣＡＩＰＩ）方法が提案されており、かかる方法は、無視できる幾何学的ファクタと低いアーチファクトレベルで、高度に加速される３Ｄイメージングを実現する（例えば、Berkin Bilgic et al., Magn. Reson. Med. 2015, 73(6), 2152-2162）。ウェーブＣＡＩＰＩ法はそれぞれのｋｘ符号化プロファイルの読み出し中に、位相及びスライス符号化方向に沿って正弦波変調された磁場勾配（ＧＹ、ＧＺ）を生成することを含む。その結果得られた取得は、エイリアシングを全ての空間方向に均等に拡げ、それにより、高レベルのサブサンプリングによる加速パラレルイメージングのための３次元コイル感度プロファイルを完全に利用する。勾配変調のボクセル拡散効果は、データグリッドを必要としない効率的な再構成スキームが使用されることができるように、画像空間における（変換された）点拡がり関数による畳み込みに対応することが示されている。
本発明によれば、ＭＲ装置の検査ボリューム内に配置された対象のＭＲ撮像方法が開示される。この方法は、対象を撮像シーケンスに曝すことによりＭＲ信号を生成するステップと、デカルトｋ空間サンプリングスキームでＭＲ信号プロファイルを取得するステップであって、各ＭＲ信号プロファイルが、読み出し方向に沿った時間的に一定の磁場勾配及び位相符号化方向に沿った正弦波変調された磁場勾配の存在下で取得される、ステップと、磁場勾配の変調スキームを考慮してＭＲ画像を再構成するステップと、を有する。本発明は、磁場勾配の正弦変調の周波数が各ＭＲ信号プロファイルの取得中に変化されることを提案する。
本発明の技術は一般に、既知のウェーブＣＡＩＰＩアプローチに対応する。しかしながら、本発明によれば、読み出し方向に沿ったエイリアシング伝播を更に改善し、潜在的な渦電流関連の悪影響を低減し、（変換された）点拡がり関数（ＰＳＦ）のサイドローブを低減するために、ｋ空間軌道を符号化する可変周波数変調波が提案される。ウェーブＣＡＩＰＩ技術に関する既知の最適化戦略は、ウェーブ符号化磁場勾配のより大きな振幅が、勾配変調のサイクル数を調整することと比較して、エイリアシング広がり効果を改善し、幾何学的因子を低減するために重要な役割を果たすことを実証した（Polak et al., ISMRM 2018, p. 0937を参照）。しかしながら、従来の一定周波数ウェーブ符号化のＰＳＦは、依然として、読み出し方向に沿った大きな振動と大きな振幅のサイドローブを持つ分布に従う。更に、磁場勾配振幅を増加させることは、渦電流の問題を導入し、特に、各ＭＲ信号プロファイルの取得間隔の終了時に、ウェーブ符号化勾配が大きな振幅で停止する方向に沿って（例えば、余弦波の形の勾配波形の場合）、問題が生じる。これらの残差渦電流はＲＦスライス選択プロファイルに更に悪影響を及ぼし、その結果、局所的な画像アーチファクト及び大域的な信号強度損失が生じる。このようなアーチファクト及び一定振幅及び／又は一定周波数のウェーブ符号化に関連する他の問題は、本発明によって克服される。
本発明に従って適用される撮像シーケンスは、任意の従来の２次元又は３次元撮像シーケンス、又は空間―時間（時間分解）撮像シーケンスなどの多次元拡張とすることができる。撮像シーケンスは、２次元、３次元又はより高次元のスピンエコーシーケンス、好ましくはターボスピンエコーシーケンスとすることができる。本発明の方法はまた、２次元、３次元又はより高次元の勾配エコーシーケンス、好ましくはターボフィールドエコーシーケンスと組み合わせて適用されることができる。より高い次元は、スペクトル次元、又は上で言及されるように、時間次元を含むことができる。このシーケンスは、化学的シフト（例えば、水分／脂肪分離）、緩和時間（Ｔ１、Ｔ２、Ｔ２
＊）、呼吸状態、心周期、ｂ値及び方向、拡散などを調査するために設計されることができる。２次元の場合、磁場勾配変調による本発明の実現は、既知のBunch符号化（ＢＰＥ）技術と等価である。
本発明の好ましい実施形態では、磁場勾配の正弦波変調の瞬時周波数が、各ＭＲ信号プロファイルの取得時間隔の前半の間に増加し、取得時間隔の後半の間に減少してその初期値に戻る。この（ミラー状の）周波数変動スキームにおいて、コルクスクリュー形状の読み出し軌道は、中央ｋ空間の近傍でより頻繁に回転し、周辺ｋ空間に近づくと反対の挙動をする。このようにして、ｋ空間サンプリング間隔は、可変密度サンプリングを模倣する読み出し方向に沿って変化し、サブサンプリング誘導エイリアシングをより均等に拡散する一方で、ＰＳＦの最大サイドローブを低減する。
更なる好ましい実施形態では、磁場勾配変調の振幅が、ＭＲ信号の取得中に変化される。振幅は例えば、各ＭＲ信号プロファイルの取得時間隔の開始時及び終了時に最小となるように変調することができる。ウェーブ符号化勾配波形に対する追加の振幅変調を採用することで、渦電流を低減し、従って、スライス選択プロファイルを改善することができる。
本発明の好ましい実施形態によれば、ＭＲ信号は、画像スライスの面内方向におけるサブサンプリングを伴って取得される。画像スライスのＭＲ画像は、この場合、SENSE、SMASH又はGRAPPAのようなそれ自体公知のパラレル画像再構成アルゴリズムを使用することによって再構成されることができる。
これまでに説明した本発明の方法は、検査ボリューム内に均一な静磁場を生成する少なくとも１つの主磁石コイルと、検査ボリューム内の異なる空間方向に切り替え磁場勾配を生成する複数の勾配コイル、身体からのＭＲ信号を受信するＲＦコイルの組、時間的に連続するＲＦパルス及び切り替え磁場勾配を制御する制御ユニットと、再構成ユニットと、を有するＭＲ装置によって実行されることができる。本発明の方法は例えば、ＭＲ装置の再構成ユニット及び／又は制御ユニットの対応するプログラムによって実現されることができる。
本発明の方法は、現在臨床上使用されているほとんどのＭＲ装置において有利に実行されることができる。この目的のために、本発明の上述の方法ステップを実行するようにＭＲ装置が制御されるコンピュータプログラムを利用することのみが必要である。コンピュータプログラムは、ＭＲ装置の制御ユニットにインストールするためにダウンロードされるように、データ担体上に存在してもよいし、データネットワーク内に存在してもよい。
添付の図面は、本発明の好ましい実施形態を開示する。しかしながら、図面は説明のみを目的として設計されたものであり、本発明の限定の定義として設計されたものではないことを理解されたい。図面において：
図１を参照すると、ＭＲ装置１が示されている。この装置は検査ボリュームを通ってｚ軸に沿って実質的に均一で時間的に一定の主磁場が生成されるように、超電導又は抵抗性の主磁石コイル２を備えている。
磁気共鳴生成及び操作システムは一連のＲＦパルス及び切り替え磁場勾配を適用して、核磁気スピンを反転又は励起し、磁気共鳴を誘導し、磁気共鳴をリフォーカスし、磁気共鳴を操作し、磁気共鳴を空間的に及び他の方法で符号化し、スピンを飽和させ、及び同様のものを行ってＭＲ撮像を行う。
より具体的には、勾配パルス増幅器３が、検査ボリュームのｘ、ｙ及びｚ軸に沿って、全身勾配コイル４、５及び６の選択されたものに電流パルスを印加する。デジタルＲＦ周波数送信器７は、送信／受信スイッチ８を介して、ＲＦパルス又はパルスパケットを全身ボリュームＲＦコイル９に送信して、ＲＦパルスを検査ボリュームに送信する。典型的なＭＲ撮像シーケンスは、短い持続時間のＲＦパルスセグメントのパケットから構成されており、これらのセグメントは、互いに一緒になって、任意の印加された磁場勾配と共に、核磁気共鳴の選択された操作を達成する。ＲＦパルスは、共鳴を飽和し、共鳴を励起し、磁化を反転し、共鳴をリフォーカスし、又は共鳴を操作するとともに、検査ボリューム内に配置された身体１０の一部を選択するために使用される。
パラレルイメージングによる身体１０の限られた領域のＭＲ画像の生成のために、一組のローカルＲＦコイル１１、１２、１３が、撮像のために選択された領域に隣接して配置される。
得られたＭＲ信号は、ＲＦコイル１１、１２、１３によってピックアップされ、好ましくは前置増幅器（図示せず）を有する受信器１４によって復調される。受信器１４は、送受信スイッチ８を介してＲＦコイル９、１１、１２及び１３に接続される。
ホストコンピュータ１５は、勾配パルス増幅器３及び送信器７を制御して、エコープラナー画像、エコーボリューム画像、勾配及びスピンエコー画像、高速スピンエコー画像等のような複数のＭＲ画像シーケンスのいずれかを生成する。選択されたシーケンスに関して、受信器１４は、各ＲＦ励起パルスに続いて、単一又は複数のＭＲデータラインを迅速に連続して受信する。データ収集システム１６は受信信号のアナログ－デジタル変換を実行し、各ＭＲデータラインを更なる処理に適したデジタル形式に変換する。今日のＭＲ装置では、データ収集システム１６は、生画像データの収集に特化した別個のコンピュータである。最終的に、デジタルの生画像データは、フーリエ変換又は他の適切な再構成アルゴリズムを適用する再構成プロセッサ１７によって画像表現に再構成される。ＭＲ画像は、患者を横切る平面スライス、平行な平面スライスのアレイ、３次元ボリュームなどを表すことができる。次いで、画像は画像メモリに記憶され、画像メモリは、画像表現のスライス、投影、又は他の部分を、例えば、結果として得られるＭＲ画像の人間可読の表示を提供するビデオモニタ１８を介して、視覚化のための適切なフォーマットに変換するためにアクセスされることができる。
引き続き図１を参照し、図２～図４を更に参照して、本発明の撮像アプローチを説明する。
本発明によれば、撮像シーケンス（例えば、既知のＴＳＥシーケンスのようなスピンエコーシーケンス）が適用される。ＭＲ信号は、異なる空間感度を有するＲＦコイル１１、１２、１３を介してパラレルに受信される。ＭＲ信号プロファイルは、位相符号化方向にサブサンプリングを持つデカルトｋ空間スキームで取得される。ウェーブＣＡＩＰＩ技術によれば、各ＭＲ信号プロファイルは、読み出し方向に沿った時間的に一定の磁場勾配と、位相符号化方向に沿った正弦波変調された磁場勾配の存在下で取得される。言い換えれば、ＭＲ信号取得は、位相符号化方向Ｇｙ及びＧｚに沿って正弦的に変化する磁場勾配と、線形読み出し磁場勾配Ｇｘを使用するデカルト読み出しの重ね合わせである。
本発明によれば、磁場勾配の正弦波変調の周波数は、各ＭＲ信号プロファイルの取得中に変化される。磁場勾配変調の瞬時周波数は一般に、ｆ（ｔ）＝ｆｃ＋ｈ（ｔ），ｔ∈［０，Ｔａｃｑ］として定義することができ、ここで、ｆｃ及びｈ（ｔ）は、それぞれ、一定周波数成分及び時間変化する周波数成分に対応する。ｈ（ｔ）は、対称、非対称、区分線形、非線形関数の任意の種類の関数でよい。本発明の実施形態では、対称及び非線形関数が、瞬時周波数、すなわち、ｆ（ｔ）＝ｆｃ－２πｆｍβｃｏｓ（２πｆｍｔ）を変調するために使用される。ここで、パラメータｆｍ及びβは、周波数変動パターン及び波形形状をそれぞれ制御する。典型的には、ｆｍ＝１／Ｔａｃｑが、収集時間隔Ｔａｃｑの前半又は後半の間に単調増加又は減少特性を達成するために選択され、一方、βは、パフォーマンスを最適化し、最大勾配強度及びスルーレートを含むハードウェア制約に適応するための調整可能なパラメータである。これにより、ウェーブ符号化ｋ空間軌跡（Ｃ１（ｔ）及びＣ２（ｔ））と、対応する勾配波形（ｇｙ（ｔ）及びｇｚ（ｔ））は、以下のように定式化されることができる：
ここで、Ｇｙ及びＧｚは、それぞれｙ及びｚ方向に沿った最大勾配振幅を示す。上記の理想的な勾配波形は、勾配変調伝達関数（すなわち、渦電流補償）の効果のために補正され、次いで、ＭＲスキャンにおいて適用され得る。また、３Ｄ ＴＳＥスキャンのためのCarr-Purcell-Meiboom-Gill（ＣＰＭＧ）条件を維持するために、追加の事前位相勾配及びリワインディング勾配が適用されることもできる。
ＭＲ画像再構成のために、正弦波勾配ＧｙとＧｚによって生成される符号化が、分離可能な点拡がり関数（ＰＳＦ）を通じて補足されることができることが利用される。基礎となる３次元画像内の各ＭＲ信号プロファイルは、ｙ方向及びｚ方向の空間的位置に依存するＰＳＦと畳み込まれ、取得されたＭＲ画像が得られる。ハイブリッド（ｋｘ，ｙ，ｚ）空間では，ＭＲ信号データの各ｋｘ線が所定の（ｙ，ｚ）位置にあるとき，ＰＳＦに含まれる固有の位相変調を得ることができる。
図２は、提案されている可変周波数のウェーブ符号化磁界勾配波形の実施例（図２のアルファベットＶで示す）を、一定周波数のウェーブ符号化波形（アルファベットＣで示す）と比較して示したものである。勾配波形は、図２（ａ）に示されているが、対応するｋ空間の軌跡は、図２（ｂ）に示されており、各時間は一定周波数及び可変周波数のウェーブ符号化スキームである。図２（ｂ）のウェーブ符号化ｋ空間軌跡は、インビボイメージングから較正される。図２（ａ）の対応する勾配波形は、数値微分とスケーリングによってシミュレートされる。
較正したウェーブ符号化ｋ空間軌道と採用されるサブサンプリングマスクにより、（変換した）ＰＳＦが合成されることができる。図３は、一定周波数（Ｃ）及び可変周波数（Ｖ）のウェーブ符号化スキームに関するＰＳＦの比較を示す。図３ａは、ｙ方向及びｚ方向の両方における減少係数２と、完全にサンプリングされる中央領域とを有する、使用されたＣＡＩＰＩサブサンプリングマスクを示す。図３ｂ及び図３ｃは、２つの異なる（ｙ，ｚ）位置について得られるＰＳＦを示す。図３から、本発明の可変周波数のウェーブ符号化スキームは、読み出し方向に沿ったエイリアシング伝播を改善し、ＰＳＦのサイドローブの振幅を著しく減少させることが分かる。
図４ａと図４ｂは、それぞれ脳スキャンの中央と横方向のスライスである。（「Full（フル）」で示される）図４ａ及び図４ｂの上段には、完全にサンプリングされたＭＲ信号データから再構成されたスライス画像が示されている。従来のデカルトサンプリングスキームは、左列の画像（「Cartesian（デカルト）」によって示される）に使用されている。中央の列（「Wave C（ウェーブＣ）」で示される）の画像には一定周波数のウェーブ符号化が用いられ、右の列（「Wave V（ウェーブＶ）」で示される）の画像には可変周波数のウェーブ符号化が用いられた。図４ａ及び図４ｂの最下段（「Sub（サブ）」によって示される）の画像は、位相符号化方向における低減係数３及びスライス符号化方向における低減係数２によりサブサンプリングされたＭＲ信号データから再構成される。図４ａ及び図４ｂから分かるように、３ｘ２ウェーブＣＡＩＰＩインビボ脳加速実験において、本発明の可変周波数のウェーブ符号化アプローチは、中央及び横方向の画像スライスにおいて、より優れたエイリアシング抑制結果を提供することができる。残差エイリアシングアーチファクトが可変周波数のウェーブ符号化により低減されることは、図４ａの最下段において明瞭に見ることができる。更に、可変周波数ウェーブ符号化は、横方向スライスにおける渦電流誘起スライスプロファイル劣化による信号損失を著しく低減することができる。明確な信号低下は、右側の列には存在しない図４ｂの中央の列に明確に見ることができる。
Claims (9)
- ＭＲ装置の検査ボリューム内に配置される対象のＭＲ撮像方法であって、
前記対象を撮像シーケンスに曝すことによってＭＲ信号を生成するステップと、
デカルトｋ空間スキームでＭＲ信号プロファイルを取得するステップであって、各ＭＲ信号プロファイルが、読み出し方向に沿った時間的に一定の磁場勾配及び位相符号化方向に沿った正弦波変調される磁場勾配の存在下で取得される、ステップと、
前記磁場勾配の変調スキームを考慮してＭＲ画像を再構成するステップと、
を有し、前記磁場勾配の正弦波変調の周波数が、各ＭＲ信号プロファイルの取得中に変化される、方法。 - 前記撮像シーケンスが、２次元、３次元、又はより高い次元のスピンエコーシーケンスであり、又はターボスピンエコーシーケンスである、請求項１に記載の方法。
- 前記撮像シーケンスは、２次元、３次元、又はより高い次元の勾配エコーシーケンスであり、又はターボフィールドエコーシーケンスである、請求項１に記載の方法。
- 前記磁場勾配の正弦波変調の瞬時周波数は、各ＭＲ信号プロファイルの前記取得時間隔の前半の間に増加し、前記取得時間隔の後半の間に減少してその初期値に戻る、請求項１乃至３のいずれか１項に記載の方法。
- 前記磁場勾配変調の振幅が、前記ＭＲ信号の取得中に変化される、請求項１乃至４のいずれか１項に記載の方法。
- 前記ＭＲ信号は、異なる空間感度プロファイルを有する少なくとも２つのＲＦコイルを通じたサブサンプリングにより取得され、前記ＭＲ画像は、SENSE、SMASH、GRAPPA又はパラレルイメージング再構成アルゴリズムを使用して再構成される、請求項１乃至５のいずれか１項に記載の方法。
- 前記撮像シーケンスは、２又はそれより多くの空間的に分離した画像スライスを同時に励起するためのマルチスライスＲＦパルスを有し、前記異なる画像スライスからのＭＲ信号の寄与が、前記少なくとも２つのＲＦコイルの空間感度プロファイルに基づいて分離される、請求項１乃至６のいずれか１項に記載の方法。
- 検査ボリューム内に均一な静磁場を生成する少なくとも１つの主磁石コイルと、前記検査ボリューム内の異なる空間方向に切り替え磁場勾配を生成する複数の勾配コイルと、少なくとも１又は複数のＲＦコイルと、時間的に連続するＲＦパルス及び切り替え磁場勾配を制御する制御ユニットと、再構成ユニットとを有するＭＲ装置であって、
前記対象を撮像シーケンスに曝すことによってＭＲ信号を生成するステップと、
デカルトｋ空間スキームでＭＲ信号プロファイルを取得するステップであって、各ＭＲ信号プロファイルが、読み出し方向に沿った時間的に一定の磁場勾配及び位相符号化方向に沿った正弦波変調された磁場勾配の存在下で取得される、ステップと、
前記磁場勾配の変調スキームを考慮して、前記取得されたＭＲ信号プロファイルからＭＲ画像を再構成するステップと、
を実行するように構成され、前記磁場勾配の正弦波変調の周波数が、各ＭＲ信号プロファイルの取得中に変化される。請求項１乃至７のいずれか１項に記載の方法を実行するＭＲ装置。 - ＭＲ装置において実行されるコンピュータプログラムであって、
撮像シーケンスを生成するステップと、
デカルトｋ空間スキームでＭＲ信号プロファイルを取得するステップであって、各ＭＲ信号プロファイルが、読み出し方向に沿った時間的に一定の磁場勾配とスライス及び位相符号化方向に沿った正弦波変調される磁場勾配との存在下で取得され、前記ＭＲ信号が、検査ボリューム内で異なる空間感度プロファイルを有する少なくとも２つのＲＦコイルの組を介してパラレルに受信される、ステップと、
前記磁場勾配の変調スキームを考慮して、前記取得されたＭＲ信号プロファイルからＭＲ画像を再構成するステップと、
を実行する命令を有し、前記磁場勾配の正弦波変調の周波数が各ＭＲ信号プロファイルの取得中に変化される、コンピュータプログラム。
Applications Claiming Priority (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US201962809779P | 2019-02-25 | 2019-02-25 | |
| US62/809,779 | 2019-02-25 | ||
| EP19170549.0A EP3730962A1 (en) | 2019-04-23 | 2019-04-23 | Parallel mr imaging using wave-encoding |
| EP19170549.0 | 2019-04-23 | ||
| PCT/EP2020/054159 WO2020173749A1 (en) | 2019-02-25 | 2020-02-18 | Parallel mr imaging using wave-encoding |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2022521277A true JP2022521277A (ja) | 2022-04-06 |
| JPWO2020173749A5 JPWO2020173749A5 (ja) | 2023-02-17 |
Family
ID=69526280
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2021549205A Pending JP2022521277A (ja) | 2019-02-25 | 2020-02-18 | ウェーブ符号化を用いたパラレルｍｒイメージング |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US11815577B2 (ja) |
| EP (1) | EP3931586B1 (ja) |
| JP (1) | JP2022521277A (ja) |
| WO (1) | WO2020173749A1 (ja) |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7301342B2 (en) | 2005-04-15 | 2007-11-27 | Case Western Reserve University | Bunched phase encoding (BPE) |
| DE102007035176B4 (de) * | 2007-07-27 | 2010-03-18 | Siemens Ag | Verfahren zur Aufzeichnung und Verarbeitung einer Folge von zeitlich aufeinander folgenden Bilddatensätzen sowie Magnet-Resonanz-Gerät |
| DE102013100349B4 (de) * | 2013-01-14 | 2016-05-12 | Siemens Aktiengesellschaft | Echoplanare MR-Bildgebung mit zickzack-artigen k-Raum-Trajektorien |
| US10132889B2 (en) * | 2013-05-22 | 2018-11-20 | General Electric Company | System and method for reducing acoustic noise level in MR imaging |
| WO2015036340A1 (en) * | 2013-09-10 | 2015-03-19 | Koninklijke Philips N.V. | Metal resistant mr imaging |
| CN104714199B (zh) * | 2013-12-17 | 2018-04-24 | 西门子(深圳)磁共振有限公司 | 一种磁共振成像方法和装置 |
-
2020
- 2020-02-18 EP EP20704325.8A patent/EP3931586B1/en active Active
- 2020-02-18 JP JP2021549205A patent/JP2022521277A/ja active Pending
- 2020-02-18 US US17/432,921 patent/US11815577B2/en active Active
- 2020-02-18 WO PCT/EP2020/054159 patent/WO2020173749A1/en unknown
Non-Patent Citations (2)
| Title |
|---|
| BERKIN BILGIC, ET AL.: "Wave-CAIPI for Highly Accelerated 3D Imaging", MAGNETIC RESONANCE IN MEDICINE, vol. 73, JPN6023015134, 2015, pages 2152 - 2162, ISSN: 0005039382 * |
| RICHARD D. HOGE, ET AL.: "Density Compensation Functions for Spiral MRI", MRM, vol. 38, JPN6023015133, 1997, pages 117 - 128, XP055025777, ISSN: 0005039383, DOI: 10.1002/mrm.1910380117 * |
Also Published As
| Publication number | Publication date |
|---|---|
| EP3931586A1 (en) | 2022-01-05 |
| US11815577B2 (en) | 2023-11-14 |
| EP3931586B1 (en) | 2022-06-22 |
| US20220155396A1 (en) | 2022-05-19 |
| WO2020173749A1 (en) | 2020-09-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3635425B1 (en) | Mr imaging using a stack-of-stars acquisition with variable contrast | |
| EP3191862B1 (en) | Zero echo time mr imaging | |
| US9159145B2 (en) | Fast dual contrast MR imaging | |
| JP2023109791A (ja) | 並列マルチスライスｍｒ撮像 | |
| EP3447517A1 (en) | Dixon-type water/fat separation mr imaging | |
| CN107810425B (zh) | 消除非t2加权信号贡献的t2加权mr成像 | |
| EP3517988A1 (en) | Mr imaging using a stack-of-stars acquisition with intrinsic motion correction | |
| WO2018114554A1 (en) | Dixon-type water/fat separation mr imaging | |
| JP2019535435A (ja) | プロペラｍｒイメージング | |
| US11959986B2 (en) | MR imaging with spiral acquisition | |
| EP2581756A1 (en) | MR imaging using parallel signal acquisition | |
| EP3736593A1 (en) | Dual echo steady state mr imaging using bipolar diffusion gradients | |
| EP3748385A1 (en) | Optimized k-space profile ordering for 3d radial or spiral mr imaging | |
| EP3931588A1 (en) | Epi mr imaging with distortion correction | |
| JP2022521277A (ja) | ウェーブ符号化を用いたパラレルｍｒイメージング | |
| EP3432019A1 (en) | Parallel multi-slice mr imaging using signal averaging | |
| EP4012434A1 (en) | Dixon-type water/fat separation mr imaging | |
| EP3730962A1 (en) | Parallel mr imaging using wave-encoding | |
| US20220229139A1 (en) | Multi-echo mr imaging with spiral acquisition | |
| EP3792647A1 (en) | Dixon-type water/fat separation mr imaging | |
| WO2022096545A1 (en) | Spin echo mr imaging with spiral acquisition |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230209 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230209 |
|
| A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20230209 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20230418 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20230714 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20231016 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20231212 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20240305 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240611 |