Optimized flow compensation for contrast-enhanced T1-weighted Wave-CAIPI 3D MPRAGE imaging of the brain

Flow-related artifacts have been observed in highly accelerated T1-weighted contrast-enhanced wave-controlled aliasing in parallel imaging (CAIPI) magnetization-prepared rapid gradient-echo (MPRAGE) imaging and can lead to diagnostic uncertainty. We developed an optimized flow-mitigated Wave-CAIPI MPRAGE acquisition protocol to reduce these artifacts through testing in a custom-built flow phantom. In the phantom experiment, maximal flow artifact reduction was achieved with the combination of flow compensation gradients and radial reordered k-space acquisition and was included in the optimized sequence. Clinical evaluation of the optimized MPRAGE sequence was performed in 64 adult patients, who all underwent contrast-enhanced Wave-CAIPI MPRAGE imaging without flow-compensation and with optimized flow-compensation parameters. All images were evaluated for the presence of flow-related artifacts, signal-to-noise ratio (SNR), gray-white matter contrast, enhancing lesion contrast, and image sharpness on a 3-point Likert scale. In the 64 cases, the optimized flow mitigation protocol reduced flow-related artifacts in 89% and 94% of the cases for raters 1 and 2, respectively. SNR, gray-white matter contrast, enhancing lesion contrast, and image sharpness were rated as equivalent for standard and flow-mitigated Wave-CAIPI MPRAGE in all subjects. The optimized flow mitigation protocol successfully reduced the presence of flow-related artifacts in the majority of cases.Relevance statementAs accelerated MRI using novel encoding schemes become increasingly adopted in clinical practice, our work highlights the need to recognize and develop strategies to minimize the presence of unexpected artifacts and reduction in image quality as potential compromises to achieving short scan times.Key points• Flow-mitigation technique led to an 89-94% decrease in flow-related artifacts.• Image quality, signal-to-noise ratio, enhancing lesion conspicuity, and image sharpness were preserved with the flow mitigation technique.• Flow mitigation reduced diagnostic uncertainty in cases where flow-related artifacts mimicked enhancing lesions.

Motion guidance lines for robust data consistency-based retrospective motion correction in 2D and 3D MRI

PURPOSE

To develop a robust retrospective motion-correction technique based on repeating k-space guidance lines for improving motion correction in Cartesian 2D and 3D brain MRI.

METHODS

The motion guidance lines are inserted into the standard sequence orderings for 2D turbo spin echo and 3D MPRAGE to inform a data consistency-based motion estimation and reconstruction, which can be guided by a low-resolution scout. The extremely limited number of required guidance lines are repeated during each echo train and discarded in the final image reconstruction. Thus, integration within a standard k-space acquisition ordering ensures the expected image quality/contrast and motion sensitivity of that sequence.

RESULTS

Through simulation and in vivo 2D multislice and 3D motion experiments, we demonstrate that respectively 2 or 4 optimized motion guidance lines per shot enables accurate motion estimation and correction. Clinically acceptable reconstruction times are achieved through fully separable on-the-fly motion optimizations (˜1 s/shot) using standard scanner GPU hardware.

CONCLUSION

The addition of guidance lines to scout accelerated motion estimation facilitates robust retrospective motion correction that can be effectively introduced without perturbing standard clinical protocols and workflows.

Clinical Evaluation of Scout Accelerated Motion Estimation and Reduction Technique for 3D MR Imaging in the Inpatient and Emergency Department Settings

BACKGROUND AND PURPOSE

A scout accelerated motion estimation and reduction (SAMER) framework has been developed for efficient retrospective motion correction. The goal of this study was to perform an initial evaluation of SAMER in a series of clinical brain MR imaging examinations.

MATERIALS AND METHODS

Ninety-seven patients who underwent MR imaging in the inpatient and emergency department settings were included in the study. SAMER motion correction was retrospectively applied to an accelerated T1-weighted MPRAGE sequence that was included in brain MR imaging examinations performed with and without contrast. Two blinded neuroradiologists graded images with and without SAMER motion correction on a 5-tier motion severity scale (none = 1, minimal = 2, mild = 3, moderate = 4, severe = 5).

RESULTS

The median SAMER reconstruction time was 1 minute 47 seconds. SAMER motion correction significantly improved overall motion grades across all examinations (P < .005). Motion artifacts were reduced in 28% of cases, unchanged in 64% of cases, and increased in 8% of cases. SAMER improved motion grades in 100% of moderate motion cases and 75% of severe motion cases. Sixty-nine percent of nondiagnostic motion cases (grades 4 and 5) were considered diagnostic after SAMER motion correction. For cases with minimal or no motion, SAMER had negligible impact on the overall motion grade. For cases with mild, moderate, and severe motion, SAMER improved the motion grade by an average of 0.3 (SD, 0.5), 1.1 (SD, 0.3), and 1.1 (SD, 0.8) grades, respectively.

CONCLUSIONS

SAMER improved the diagnostic image quality of clinical brain MR imaging examinations with motion artifacts. The improvement was most pronounced for cases with moderate or severe motion.

Deep learning-based motion quantification from k-space for fast model-based magnetic resonance imaging motion correction

BACKGROUND

Intra-scan rigid-body motion is a costly and ubiquitous problem in clinical magnetic resonance imaging (MRI) of the head.

PURPOSE

State-of-the-art methods for retrospective motion correction in MRI are often computationally expensive or in the case of image-to-image deep learning (DL) based methods can be prone to undesired alterations of the image (hallucinations'). In this work we introduce a novel rigid-body motion correction method which combines the advantages of classical model-driven and data-consistency (DC) preserving approaches with a novel DL algorithm, to provide fast and robust retrospective motion correction.

METHODS

The proposed Motion Parameter Estimating Densenet (MoPED) retrospectively estimates subject head motion during MRI acquisitions using a DL network with DenseBlocks and multitask learning. It quantifies the 2D rigid in-plane motion parameters slice-wise for each echo train (ET) of a Cartesian T2-weighted 2D Turbo-Spin-Echo sequence. The network receives a center patch of the motion corrupted k-space as well as an additional motion-free low-resolution reference scan to provide the ground truth orientation. The supervised training utilizes motion simulations based on 28 acquisitions with subject-wise training, validation, and test data splits of 70%, 23%, and 7%. During inference, MoPED is embedded in an iterative DC-driven motion correction algorithm which alternatingly updates estimates of the motion parameters and motion-corrected low-resolution k-space data. The estimated motion parameters are then used to reconstruct the final motion corrected image. The mean absolute/squared error and the Pearson correlation coefficient were used to analyze the motion parameter estimation quality on in-silico data in a quantitative evaluation. Structural similarity (SSIM), DC error and root mean squared error (RMSE) were used as metrics of image quality improvement. Furthermore, the generalization capability of the network was analyzed on two in-vivo motion volumes with 28 slices each and on one simulated T1-weighted volume.

RESULTS

The motion estimation achieves a Pearson correlation of 0.968 to the simulated ground-truth of the 2433 test data slices used. In-silico results indicate that MoPED decreases the time for the optimization by a factor of around 27 compared to a conventional method and is able to reduce the RMSE of the reconstructions and average DC error by more than a factor of two compared to uncorrected images. In-vivo experiments show a decrease in computation time by a factor of around 20, a RMSE decrease from 0.055 to 0.033 and an SSIM increase from 0.795 to 0.862. Furthermore, contrast independence is demonstrated as MoPED is also able to correct T1-weighted images in simulations without retraining. Due to the model-based correction, no hallucinations were observed.

CONCLUSIONS

Incorporating DL in a model-based motion correction algorithm shows great benefit on the optimization and computation time. The k-space-based estimation also allows a data consistent correction and therefore avoids the risk of hallucinations of image-to-image approaches.

Scout accelerated motion estimation and reduction (SAMER)

PURPOSE

To demonstrate a navigator/tracking-free retrospective motion estimation technique that facilitates clinically acceptable reconstruction time.

METHODS

Scout accelerated motion estimation and reduction (SAMER) uses a single 3-5 s, low-resolution scout scan and a novel sequence reordering to independently determine motion states by minimizing the data-consistency error in a SENSE plus motion forward model. This eliminates time-consuming alternating optimization as no updates to the imaging volume are required during the motion estimation. The SAMER approach was assessed quantitatively through extensive simulation and was evaluated in vivo across multiple motion scenarios and clinical imaging contrasts. Finally, SAMER was synergistically combined with advanced encoding (Wave-CAIPI) to facilitate rapid motion-free imaging.

RESULTS

The highly accelerated scout provided sufficient information to achieve accurate motion trajectory estimation (accuracy ~0.2 mm or degrees). The novel sequence reordering improved the stability of the motion parameter estimation and image reconstruction while preserving the clinical imaging contrast. Clinically acceptable computation times for the motion estimation (~4 s/shot) are demonstrated through a fully separable (non-alternating) motion search across the shots. Substantial artifact reduction was demonstrated in vivo as well as corresponding improvement in the quantitative error metric. Finally, the extension of SAMER to Wave-encoding enabled rapid high-quality imaging at up to R = 9-fold acceleration.

CONCLUSION

SAMER significantly improved the computational scalability for retrospective motion estimation and correction.

Direct magnetic field estimation based on echo planar raw data

Gradient recalled echo echo planar imaging is widely used in functional magnetic resonance imaging. The fast data acquisition is, however, very sensitive to field inhomogeneities which manifest themselves as artifacts in the images. Typically used correction methods have the common deficit that the data for the correction are acquired only once at the beginning of the experiment, assuming the field inhomogeneity distribution B(0) does not change over the course of the experiment. In this paper, methods to extract the magnetic field distribution from the acquired k-space data or from the reconstructed phase image of a gradient echo planar sequence are compared and extended. A common derivation for the presented approaches provides a solid theoretical basis, enables a fair comparison and demonstrates the equivalence of the k-space and the image phase based approaches. The image phase analysis is extended here to calculate the local gradient in the readout direction and improvements are introduced to the echo shift analysis, referred to here as "k-space filtering analysis." The described methods are compared to experimentally acquired B(0) maps in phantoms and in vivo. The k-space filtering analysis presented in this work demonstrated to be the most sensitive method to detect field inhomogeneities.

SENSE shimming (SSH): A fast approach for determining B(0) field inhomogeneities using sensitivity coding

The pursuit of ever higher field strengths and faster data acquisitions has led to the construction of coil arrays with high numbers of elements. With the sensitivity encoding (SENSE) technique, it has been shown that the sensitivity of those elements can be used for spatial image encoding. Here, a proof-of-principle is presented of a method that can be considered an extreme case of the SENSE approach, completely abstaining from using encoding gradients. The resulting sensitivity encoded free-induction decay (FID) data are then not used for imaging, but for determining B(0) field inhomogeneity distribution. The method has therefore been termed "SENSE shimming" (SSH). In phantom experiments the method's ability to detect inhomogeneities of up to the second order is demonstrated.