Fat navigators and Moiré phase tracking comparison for motion estimation and retrospective correction

Analysis of the effect of motion on highly accelerated 3D FatNavs in 3D brain images acquired at 3T

PURPOSE

3D FatNavs are rapid acquisitions of MRI fat-volumes within the head that can be used for retrospective motion correction for brain MRI. 3D FatNavs typically use very high acceleration factors and are reconstructed with the GRAPPA parallel imaging technique. However, the GRAPPA reconstruction is not expected to perform well on 3D FatNavs volumes in the presence of strong motion due to the mismatched calibration data acquired once at the start of the scan, leading to motion-parameter misestimation. This study aims to assess the accuracy and precision of 3D FatNav-derived motion-estimates in the presence of large changes in head position.

METHODS

Rigid motion parameters were simulated and applied retrospectively to the 3D FatNav volumes from MPRAGE datasets acquired at 3T. The transformed images were then re-reconstructed using GRAPPA to simulate real motion deterioration of the fat-navigator, and used to estimate the motion applied and evaluate the tracking inaccuracy. This information was then used to estimate the residual motion after 3D FatNav-based motion correction and applied to the original MPRAGE volumes. The effect of the misestimation was assessed using an image quality metric and the evaluation scores from two observers. Quality boundaries were then estimated to assess the motion tolerance when 3D FatNavs are used.

RESULTS

The GRAPPA reconstruction was shown to deteriorate for large changes in the head position, affecting the quality of 3D FatNav volumes and consequently degrading the accuracy of the motion-estimates. Based on our simulations, the estimated threshold of motion that led to a noticeable degradation in the MPRAGE image quality was up to RMS values of 3.7° and 3 mm for rotations and translations respectively.

CONCLUSIONS

3D FatNavs were shown to be able to correct for a wide range of motion levels and types. Boundaries of acceptable motion magnitudes for different levels of acceptable loss of image quality were determined.

Rapid and accurate navigators for motion and B0 tracking using QUEEN: Quantitatively enhanced parameter estimation from navigators

PURPOSE

To develop a framework that jointly estimates rigid motion and polarizing magnetic field (B0 ) perturbations (δ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ ) for brain MRI using a single navigator of a few milliseconds in duration, and to additionally allow for navigator acquisition at arbitrary timings within any type of sequence to obtain high-temporal resolution estimates.

THEORY AND METHODS

Methods exist that match navigator data to a low-resolution single-contrast image (scout) to estimate either motion orδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ . In this work, called QUEEN (QUantitatively Enhanced parameter Estimation from Navigators), we propose combined motion andδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ estimation from a fast, tailored trajectory with arbitrary-contrast navigator data. To this end, the concept of a quantitative scout (Q-Scout) acquisition is proposed from which contrast-matched scout data is predicted for each navigator. Finally, navigator trajectories, contrast-matched scout, andδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ are integrated into a motion-informed parallel-imaging framework.

RESULTS

Simulations and in vivo experiments show the need to modelδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ to obtain accurate motion parameters estimated in the presence of strongδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ . Simulations confirm that tailored navigator trajectories are needed to robustly estimate both motion andδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ . Furthermore, experiments show that a contrast-matched scout is needed for parameter estimation from multicontrast navigator data. A retrospective, in vivo reconstruction experiment shows improved image quality when using the proposed Q-Scout and QUEEN estimation.

CONCLUSIONS

We developed a framework to jointly estimate rigid motion parameters andδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ from navigators. Combing a contrast-matched scout with the proposed trajectory allows for navigator deployment in almost any sequence and/or timing, which allows for higher temporal-resolution motion andδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ estimates.

Robust retrospective motion correction of head motion using navigator-based and markerless motion tracking techniques

PURPOSE

This study investigated the artifacts arising from different types of head motion in brain MR images and how well these artifacts can be compensated using retrospective correction based on two different motion-tracking techniques.

METHODS

MPRAGE images were acquired using a 3 T MR scanner on a cohort of nine healthy participants. Subjects moved their head to generate circular motion (4 or 6 cycles/min), stepwise motion (small and large) and "simulated realistic" motion (nodding and slow diagonal motion), based on visual instructions. One MPRAGE scan without deliberate motion was always acquired as a "no motion" reference. Three dimensional fat-navigator (FatNavs) and a Tracoline markerless device (TracInnovations) were used to obtain motion estimates and images were separately reconstructed retrospectively from the raw data based on these different motion estimates.

RESULTS

Image quality was recovered from both motion tracking techniques in our stepwise and slow diagonal motion scenarios in almost all cases, with the apparent visual image quality comparable to the no-motion case. FatNav-based motion correction was further improved in the case of stepwise motion using a skull masking procedure to exclude non-rigid motion of the neck from the co-registration step. In the case of circular motion, both methods struggled to correct for all motion artifacts.

CONCLUSION

High image quality could be recovered in cases of stepwise and slow diagonal motion using both motion estimation techniques. The circular motion scenario led to more severe image artifacts that could not be fully compensated by the retrospective motion correction techniques used.

Mesoscopic in vivo human T2* dataset acquired using quantitative MRI at 7 Tesla

Mesoscopic (0.1-0.5 mm) interrogation of the living human brain is critical for advancing neuroscience and bridging the resolution gap with animal models. Despite the variety of MRI contrasts measured in recent years at the mesoscopic scale, in vivo quantitative imaging of T₂^* has not been performed. Here we provide a dataset containing empirical T₂^* measurements acquired at 0.35 × 0.35 × 0.35 mm³ voxel resolution using 7 Tesla MRI. To demonstrate unique features and high quality of this dataset, we generate flat map visualizations that reveal fine-scale cortical substructures such as layers and vessels, and we report quantitative depth-dependent T₂^* (as well as R₂^*) values in primary visual cortex and auditory cortex that are highly consistent across subjects. This dataset is freely available at https://doi.org/10.17605/OSF.IO/N5BJ7, and may prove useful for anatomical investigations of the human brain, as well as for improving our understanding of the basis of the T₂^*-weighted (f)MRI signal.

Accuracy investigations for volumetric head-motion navigators with and without EPI at 7 T

PURPOSE

Accuracy investigation of volumetric navigators for motion correction, with emphasis on geometric EPI distortions at ultrahigh field.

METHODS

High-resolution Dixon images were collected in different head positions and reconstructed to water, fat, T2 *, and B0 maps. Resolution reduction was performed, and the T2 * and B0 maps were used to apply effects of TE and EPI distortions to simulate various volumetric water and fat navigators. Registrations of the simulated navigators were compared with registrations of the original high-resolution images.

RESULTS

Increased accuracy was observed with increased spatial resolution for non-EPI navigators. When using EPI, the distortions had a negative effect on registration accuracy, which was most noticeable for high-resolution navigators. Parallel imaging helped to alleviate those caveats to a certain extent, and 5-fold acceleration gave close to similar accuracy to non-EPI in most cases. Shortening the TE by partial Fourier sampling was shown to be mostly beneficial, except for water navigators with long readout durations. The EPI blip direction had an influence on navigator accuracy, and positive blip gradient polarities (yielding mostly image stretching frontally) typically gave the best accuracy for water navigators, whereas no clear recommendation could be made for fat navigators. Generally, fat EPI navigators had lower accuracy than water EPI navigators with otherwise similar parameters.

CONCLUSIONS

Echo planar imaging has been widely used for MRI navigators, but the induced distortions reduce navigator accuracy at ultrahigh field. This study can help protocol optimization and guide the complex tradeoff between resolution and EPI acceleration in navigator parameter setup.

Motion-corrected 3D-EPTI with efficient 4D navigator acquisition for fast and robust whole-brain quantitative imaging

PURPOSE

To develop a motion estimation and correction method for motion-robust three-dimensional (3D) quantitative imaging with 3D-echo-planar time-resolved imaging.

THEORY AND METHODS

The 3D-echo-planar time-resolved imaging technique was designed with additional four-dimensional navigator acquisition (x-y-z-echoes) to achieve fast and motion-robust quantitative imaging of the human brain. The four-dimensional-navigator is inserted into the relaxation-recovery deadtime of the sequence in every pulse TR (∼2 s) to avoid extra scan time, and to provide continuous tracking of the 3D head motion and B0 -inhomogeneity changes. By using an optimized spatiotemporal encoding combined with a partial-Fourier scheme, the navigator acquires a large central k-t data block for accurate motion estimation using only four small-flip-angle excitations and readouts, resulting in negligible signal-recovery reduction to the 3D-echo-planar time-resolved imaging acquisition. By incorporating the estimated motion and B0 -inhomogeneity changes into the reconstruction, multi-contrast images can be recovered with reduced motion artifacts.

RESULTS

Simulation shows the cost to the SNR efficiency from the added navigator acquisitions is <1%. Both simulation and in vivo retrospective experiments were conducted, that demonstrate the four-dimensional navigator provided accurate estimation of the 3D motion and B0 -inhomogeneity changes, allowing effective reduction of image artifacts in quantitative maps. Finally, in vivo prospective undersampling acquisition was performed with and without head motion, in which the motion corrupted data after correction show close image quality and consistent quantifications to the motion-free scan, providing reliable quantitative measurements even with head motion.

CONCLUSION

The proposed four-dimensional navigator acquisition provides reliable tracking of the head motion and B0 change with negligible SNR cost, equips the 3D-echo-planar time-resolved imaging technique for motion-robust and efficient quantitative imaging.

Comparison of prospective and retrospective motion correction in 3D-encoded neuroanatomical MRI

PURPOSE

To compare prospective motion correction (PMC) and retrospective motion correction (RMC) in Cartesian 3D-encoded MPRAGE scans and to investigate the effects of correction frequency and parallel imaging on the performance of RMC.

METHODS

Head motion was estimated using a markerless tracking system and sent to a modified MPRAGE sequence, which can continuously update the imaging FOV to perform PMC. The prospective correction was applied either before each echo train (before-ET) or at every sixth readout within the ET (within-ET). RMC was applied during image reconstruction by adjusting k-space trajectories according to the measured motion. The motion correction frequency was retrospectively increased with RMC or decreased with reverse RMC. Phantom and in vivo experiments were used to compare PMC and RMC, as well as to compare within-ET and before-ET correction frequency during continuous motion. The correction quality was quantitatively evaluated using the structural similarity index measure with a reference image without motion correction and without intentional motion.

RESULTS

PMC resulted in superior image quality compared to RMC both visually and quantitatively. Increasing the correction frequency from before-ET to within-ET reduced the motion artifacts in RMC. A hybrid PMC and RMC correction, that is, retrospectively increasing the correction frequency of before-ET PMC to within-ET, also reduced motion artifacts. Inferior performance of RMC compared to PMC was shown with GRAPPA calibration data without intentional motion and without any GRAPPA acceleration.

CONCLUSION

Reductions in local Nyquist violations with PMC resulted in superior image quality compared to RMC. Increasing the motion correction frequency to within-ET reduced the motion artifacts in both RMC and PMC.

Quantitative evaluation of prospective motion correction in healthy subjects at 7T MRI

PURPOSE

Quantitative assessment of prospective motion correction (PMC) capability at 7T MRI for compliant healthy subjects to improve high-resolution images in the absence of intentional motion.

METHODS

Twenty-one healthy subjects were imaged at 7 T. They were asked not to move, to consider only unintentional motion. An in-bore optical tracking system was used to monitor head motion and consequently update the imaging volume. For all subjects, high-resolution T₁ (3D-MPRAGE), T₂ (2D turbo spin echo), proton density (2D turbo spin echo), and T2∗ (2D gradient echo) weighted images were acquired with and without PMC. The images were evaluated through subjective and objective analysis.

RESULTS

Subjective evaluation overall has shown a statistically significant improvement (5.5%) in terms of image quality with PMC ON. In a separate evaluation of every contrast, three of the four contrasts (T₁ , T₂ , and proton density) have shown a statistically significant improvement (9.62%, 9.85%, and 9.26%), whereas the fourth one ( T2∗ ) has shown improvement, although not statistically significant. In the evaluation with objective metrics, average edge strength has shown an overall improvement of 6% with PMC ON, which was statistically significant; and gradient entropy has shown an overall improvement of 2%, which did not reach statistical significance.

CONCLUSION

Based on subjective assessment, PMC improved image quality in high-resolution images of healthy compliant subjects in the absence of intentional motion for all contrasts except T2∗ , in which no significant differences were observed. Quantitative metrics showed an overall trend for an improvement with PMC, but not all differences were significant.

Comprehensive ultrahigh resolution whole brain in vivo MRI dataset as a human phantom

Here, we present an extension to our previously published structural ultrahigh resolution T1-weighted magnetic resonance imaging (MRI) dataset with an isotropic resolution of 250 µm, consisting of multiple additional ultrahigh resolution contrasts. Included are up to 150 µm Time-of-Flight angiography, an updated 250 µm structural T1-weighted reconstruction, 330 µm quantitative susceptibility mapping, up to 450 µm structural T2-weighted imaging, 700 µm T1-weighted back-to-back scans, 800 µm diffusion tensor imaging, one hour continuous resting-state functional MRI with an isotropic spatial resolution of 1.8 mm as well as more than 120 other structural T1-weighted volumes together with multiple corresponding proton density weighted acquisitions collected over ten years. All data are from the same participant and were acquired on the same 7 T scanner. The repository contains the unprocessed data as well as (pre-)processing results. The data were acquired in multiple studies with individual goals. This is a unique and comprehensive collection comprising a "human phantom" dataset. Therefore, we compiled, processed, and structured the data, making them publicly available for further investigation.

Combining prospective and retrospective motion correction based on a model for fast continuous motion

PURPOSE

Prospective motion correction (PMC) and retrospective motion correction (RMC) have different advantages and limitations. The present work aims to combine the advantages of both for rigid body motion, aiming at correcting for faster motions than was previously achievable. Additionally, it provides insights into the effects of motion on pulse sequences and MR signals with a goal of further improving motion correction in the future.

METHODS

The effective encoding trajectory and a global phase offset in a moving object are calculated based on complete gradient waveforms of an arbitrary sequence and a continuous motion model. These data are used to feed the forward signal model, which is then used in iterative image reconstruction to suppress the artifacts still present after the PMC.

RESULTS

Verification experiments with a rotation phantom and in vivo were performed. Predictions of simulated motion artifacts for PMC based on sequence waveforms are very accurate. The performance at combined PMC+RMC is limited by Nyquist violations in the sampled k-space and can be compensated by oversampling.

CONCLUSION

The combined correction results in better images than pure PMC in the presence of fast motion. The predictions of artifacts are very accurate, allowing for comparing sequences or protocols in simulations. The observed artifacts due to Nyquist violations are expected to be corrected by utilizing parallel imaging.

Free induction decay navigator motion metrics for prediction of diagnostic image quality in pediatric MRI

Purpose

To investigate the ability of free induction decay navigator (FIDnav)‐based motion monitoring to predict diagnostic utility and reduce the time and cost associated with acquiring diagnostically useful images in a pediatric patient cohort.

Methods

A study was carried out in 102 pediatric patients (aged 0‐18 years) at 3T using a 32‐channel head coil array. Subjects were scanned with an FID‐navigated MPRAGE sequence and images were graded by two radiologists using a five‐point scale to evaluate the impact of motion artifacts on diagnostic image quality. The correlation between image quality and four integrated FIDnav motion metrics was investigated, as well as the sensitivity and specificity of each FIDnav‐based metric to detect different levels of motion corruption in the images. Potential time and cost savings were also assessed by retrospectively applying an optimal detection threshold to FIDnav motion scores.

Results

A total of 12% of images were rated as non‐diagnostic, while a further 12% had compromised diagnostic value due to motion artifacts. FID‐navigated metrics exhibited a moderately strong correlation with image grade (Spearman's rho ≥ 0.56). Integrating the cross‐correlation between FIDnav signal vectors achieved the highest sensitivity and specificity for detecting non‐diagnostic images, yielding total time savings of 7% across all scans. This corresponded to a financial benefit of $2080 in this study.

Conclusions

Our results indicate that integrated motion metrics from FIDnavs embedded in structural MRI are a useful predictor of diagnostic image quality, which translates to substantial time and cost savings when applied to pediatric MRI examinations.

Most Small Cerebral Cortical Veins Demonstrate Significant Flow Pulsatility: A Human Phase Contrast MRI Study at 7T

Phase contrast MRI (pcMRI) has been used to investigate flow pulsatility in cerebral arteries, larger cerebral veins, and the cerebrospinal fluid (CSF). Such measurements of intracranial pulsatility and compliance are beginning to inform understanding of the pathophysiology of conditions including normal pressure hydrocephalus, multiple sclerosis, and dementias. We demonstrate the presence of flow pulsatility in small cerebral cortical veins, for the first time using pcMRI at 7 T, with the aim of improving our understanding of the hemodynamics of this little-studied vascular compartment. A method for establishing where venous flow is pulsatile is introduced, revealing significant pulsatility in 116 out of 146 veins, across eight healthy participants, assessed in parietal and frontal regions. Distributions of pulsatility index (PI) and pulse waveform delay were characterized, indicating a small, but statistically significant (p < 0.05), delay of 59 ± 41 ms in cortical veins with respect to the superior sagittal sinus, but no differences between veins draining different arterial supply territories. Measurements of pulsatility in smaller cortical veins, a hitherto unstudied compartment closer to the capillary bed, could lead to a better understanding of intracranial compliance and cerebrovascular (patho)physiology.

Improved susceptibility-weighted imaging for high contrast and resolution thalamic nuclei mapping at 7T

PURPOSE

The thalamus is an important brain structure and neurosurgical target, but its constituting nuclei are challenging to image non-invasively. Recently, susceptibility-weighted imaging (SWI) at ultra-high field has shown promising capabilities for thalamic nuclei mapping. In this work, several methodological improvements were explored to enhance SWI quality and contrast, and specifically its ability for thalamic imaging.

METHODS

High-resolution SWI was performed at 7T in healthy participants, and the following techniques were applied: (a) monitoring and retrospective correction of head motion and B₀ perturbations using integrated MR navigators, (b) segmentation and removal of venous vessels on the SWI data using vessel enhancement filtering, and (c) contrast enhancement by tuning the parameters of the SWI phase-magnitude combination. The resulting improvements were evaluated with quantitative metrics of image quality, and by comparison to anatomo-histological thalamic atlases.

RESULTS

Even with sub-millimeter motion and natural breathing, motion and field correction produced clear improvements in both magnitude and phase data quality (76% and 41%, respectively). The improvements were stronger in cases of larger motion/field deviations, mitigating the dependence of image quality on subject performance. Optimizing the SWI phase-magnitude combination yielded substantial improvements in image contrast, particularly in the thalamus, well beyond previously reported SWI results. The atlas comparisons provided compelling evidence of anatomical correspondence between SWI features and several thalamic nuclei, for example, the ventral intermediate nucleus. Vein detection performed favorably inside the thalamus, and vein removal further improved visualization.

CONCLUSION

Altogether, the proposed developments substantially improve high-resolution SWI, particularly for thalamic nuclei imaging.