An illustrated comparison of processing methods for MR phase imaging and QSM: combining array coil signals and phase unwrapping

Quantification of tissue stiffness with magnetic resonance elastography and finite difference time domain (FDTD) simulation-based spatiotemporal neural network

Quantification of tissue stiffness with magnetic resonance elastography (MRE) is an inverse problem that is sensitive to noise. Conventional methods for the purpose include direct inversion (DI) and local frequency estimation (LFE). In this study, we propose to train a spatiotemporal neural network using MRE data simulated by the Finite Difference Time Domain method (FDTDNet), and to use the trained network to estimate tissue stiffness from MRE data. The proposed method showed significantly better robustness to noise than DI or LFE. For simulated data with signal-to-noise ratio (SNR) of 15 dB, tissue stiffness by FDTDNet had mean absolute error of 0.41 kPa or 7 %, 77.8 % and 84.4 % lower than those by DI and LFE respectively (P < 0.0001). For a homogeneous phantom with driver power decreasing from 30 % to 5 %, FDTDNet, DI and LFE provided stiffness estimates with deviation of 6.9 % (0.21 kPa), 9.2 % (0.28 kPa) and 45.8 % (1.20 kPa) of the respective stiffness level at driver power of 30 %. Detectability of small inclusions in estimated stiffness maps is also critical. For simulated data with inclusions of radius of 0.31 cm, FDTDNet achieved contrast-to-noise ratio (CNR) of 4.20, 6900 % and 347 % higher than DI and LFE respectively (P < 0.0001), and structural similarity index (SSIM) of 0.61, 27 % and 177 % higher than DI and LFE respectively (P < 0.0001). For phantom with inclusion of radius 0.39 cm, CNR of FDTDNet was 2.98, 90 % and 80 % higher than DI and LFE respectively (P < 0.0001) and SSIM was 0.80, 89 % and 28 % higher than DI and LFE respectively (P < 0.0001). We also demonstrated the feasibility of FDTDNet in MRE data acquired from calf muscles of human subjects. In conclusion, by using a spatiotemporal neural network trained with simulated data, FDTDNet estimated tissue stiffness from MRE with superior noise robustness and detectability of focal inclusions, therefore showed potential in precisely quantifying MRE of human subjects.

Determinants of long-term paramagnetic rim lesion evolution in people with multiple sclerosis

OBJECTIVE

Baseline paramagnetic rim lesion (PRL) load predicts disease progression in people with multiple sclerosis (pwMS). Understanding how PRLs relate to other known MS-related factors, and the practical utility of PRLs in clinical trials, is crucial for informing clinical decision-making and guiding development of novel disease-modifying treatments (DMTs).

METHODS

This study included 152 pwMS enrolled in a larger prospective, longitudinal cohort study who had 3T MRI scans and clinical assessments at baseline and 5- or 10-year follow-ups. PRLs were identified on baseline 3T quantitative susceptibility maps and classified as persisting, disappearing, or newly appearing at follow-up. The relationships between PRL evolution and clinical, radiological, environmental, and genetic characteristics were assessed, and clinical trial sample sizes were estimated using PRL appearance or disappearance as outcome measures.

RESULTS

DMT use was associated with lower odds of new PRL appearance (for high-efficacy DMTs: odds ratio = 0.088, p = 0.024), but not disappearance. Current smoking status was associated with greater baseline PRL number (B = 0.527 additional PRLs, p = 0.013). A 24-month clinical trial in people with progressive MS for a DMT that doubles the rate of PRL rim disappearance would require an estimated 118 people with progressive MS per group at 80% statistical power.

INTERPRETATION

Early MS diagnosis and subsequent DMT initiation may reduce new chronic active inflammation. However, the utility of PRL disappearance or new PRL appearance as outcome measures in clinical trials is limited by potentially large sample sizes that are needed for moderate efficacy drugs.

Thalamic Magnetic Susceptibility (χ) Alterations in Neurodegenerative Diseases: A Systematic Review and Meta-Analysis of Quantitative Susceptibility Mapping Studies

BACKGROUND

Quantitative Susceptibility Mapping (QSM) provides a non-invasive post-processing method to investigate alterations in magnetic susceptibility (χ), reflecting iron content within brain regions implicated in neurodegenerative diseases (NDDs).

PURPOSE

To investigate alterations in thalamic χ in patients with NDDs using QSM.

STUDY TYPE

Systematic review and meta-analysis.

POPULATION

A total of 696 patients with NDDs and 760 healthy controls (HCs) were included in 27 studies.

FIELD STRENGTH/SEQUENCE

Three-dimensional multi-echo gradient echo sequence for QSM at mostly 3 Tesla.

ASSESSMENT

Studies reporting QSM values in the thalamus of patients with NDDs were included. Following PRISMA 2020, we searched the four major databases including PubMed, Scopus, Web of Science, and Embase for peer-reviewed studies published until October 2024.

STATISTICAL TESTS

Meta-analysis was conducted using a random-effects model to calculate the standardized mean difference (SMD) between patients and HCs.

RESULTS

The pooled SMD indicated a significant increase in thalamic χ in NDDs compared to HCs (SMD = 0.42, 95% CI: 0.05-0.79; k = 27). Notably, amyotrophic lateral sclerosis patients showed a significant increase in thalamic χ (1.09, 95% CI: 0.65-1.53, k = 2) compared to HCs. Subgroup analyses revealed significant χ alterations in younger patients (mean age ≤ 62 years; 0.56, 95% CI: 0.10-1.02, k = 11) and studies using greater coil channels (coil channels > 16; 0.64, 95% CI: 0.28-1.00, k = 9). Publication bias was not detected and quality assessment indicated that studies with a lower risk of bias presented more reliable findings (0.75, 95% CI: 0.32-1.18, k = 9). Disease type was the primary driver of heterogeneity, while other factors, such as coil type and geographic location, also contributed to variability.

DATA CONCLUSION

Our findings support the potential of QSM for investigating thalamic involvement in NDDs. Future research should focus on disease-specific patterns, thalamic-specific nucleus analysis, and temporal evolution.

PLAIN LANGUAGE SUMMARY

Our research investigated changes in iron levels within the thalamus, a brain region crucial for motor and cognitive functions, in patients with various neurodegenerative diseases (NDDs). The study utilized a specific magnetic resonance imaging technique called Quantitative Susceptibility Mapping (QSM) to measure iron content. It identified a significant increase in thalamic iron levels in NDD patients compared to healthy individuals. This increase was particularly prominent in patients with Amyotrophic Lateral Sclerosis, younger individuals, and studies employing advanced imaging equipment.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 2.

Association between paramagnetic rim lesions and pulvinar iron depletion in persons with multiple sclerosis

BACKGROUND

The deep gray matter (DGM), especially the pulvinar, and the white matter surrounding chronic active lesions have demonstrated depleted iron levels, indicating a possible mechanistic link. However, no studies have investigated the potential relationship between these phenomena.

OBJECTIVES

The study aimed to determine whether PRLs were associated with pulvinar iron depletion and, if so, whether this relationship was spatially mediated.

METHODS

This retrospective analysis included 139 people with MS (pwMS) and 43 healthy controls (HCs) scanned at 3T MRI at baseline and after 5.4 ± 0.6 years. Pulvinar iron concentrations (cFe) and iron masses (mFe) were estimated from quantitative susceptibility maps and tested for associations with PRLs. A separate cohort of 96 pwMS with PRLs and propensity-matched HCs was included to evaluate peri‑plaque normal-appearing white matter (NAWM) abnormalities.

RESULTS

PRL number was associated with greater decline in pulvinar cFe (β = -0.265, p = 0.005) and mFe (β = -0.256, p = 0.006). Peri-plaque NAWM susceptibility was increased 11 mm surrounding PRLs, outside which shorter PRL-to-pulvinar distance was associated with greater decline in pulvinar cFe (β = 0.380, p = 0.005) and mFe (β = 0.348, p = 0.022).

CONCLUSIONS

Our findings support a spatially-mediated relationship between PRLs and chronic pulvinar iron depletion.

Associations Between Paramagnetic Rim Lesion Evolution and Clinical and Radiologic Disease Progression in Persons With Multiple Sclerosis

BACKGROUND AND OBJECTIVES

Recent technological advances have enabled visualizing in vivo a subset of chronic active brain lesions in persons with multiple sclerosis (pwMS), referred to as "paramagnetic rim lesions" (PRLs), with iron-sensitive MRI. PRLs predict future clinical disease progression, making them a promising clinical and translational imaging marker. However, it is unknown how disease progression is modified by PRL evolution (PRL disappearance, new PRL appearance). This is key to understanding MS pathophysiology and may help inform selection of sensitive endpoints for clinical trials targeting chronic active inflammation. To this end, we assessed the longitudinal associations between PRL disappearance and new PRL appearance and clinical disability progression and brain atrophy.

METHODS

PwMS and healthy controls (HCs) were included from a larger prospective, longitudinal cohort study at the University at Buffalo if they had available 3T MRI and clinical visits at baseline and follow-up timepoints. PwMS with sufficient clinical data for confirmed disability progression (CDP) analysis were included in a Disability Progression Cohort, and pwMS and HCs with brain volumetry data at baseline and follow-up were included in MS and HC Brain Atrophy cohorts. PRLs were assessed at baseline and follow-up and assigned as disappearing, newly appearing, or persisting at follow-up. Linear models were fit to compare annualized PRL disappearance rates or new PRL appearance (yes/no) with annualized rates of CDP and progression independent of relapse activity (PIRA) or with annualized rates of brain atrophy, adjusting for covariates including baseline PRL number and follow-up time. Statistical analyses were corrected for false discovery rate (FDR; i.e., q-value).

RESULTS

In total, 160 pwMS (73.8% female; mean baseline age 46.6 ± 11.4 years; mean baseline disease duration 13.8 ± 10.6 years; median follow-up time 5.6 [interquartile range 5.2-7.8] years; 26.9% progressive MS) and 27 HCs (74.1% female; mean baseline age 43.9 ± 13.6 years; median follow-up time 5.4 [5.2-5.6] years) were enrolled. Greater PRL disappearance rates were associated with reduced rates of CDP (β mean = -0.262, 95% CI -0.475 to -0.049, q = 0.028) and PIRA (β mean = -0.283, 95% CI -0.492 to -0.073, q = 0.036), and new PRL appearance was associated with increased rates of PIRA (β mean = 0.223, 95% CI 0.049-0.398, q = 0.036). By contrast, no associations between new PRL appearance or PRL disappearance and brain volume changes survived FDR correction (q > 0.05).

DISCUSSION

Our results show that resolution of existing PRLs and lack of new PRLs are associated with improved clinical outcomes. These findings further motivate the need for novel therapies targeting microglia-mediated brain inflammation and adoption of clinical strategies to prevent appearance of new PRL.

Combining Quantitative Susceptibility Mapping With the Gray Matter Volume to Predict Neurological Deficits in Patients With Small Artery Occlusion

BACKGROUND

Currently, there is still a lack of valuable neuroimaging markers to assess the clinical severity of stroke patients with small artery occlusion (SAO). Quantitative susceptibility mapping (QSM) is a quantitative processing method for neuroradiological diagnostics. Gray matter (GM) volume changes in stroke patients are also proved to be associated with neurological deficits. This study aims to explore the predictive value of QSM and GM volume in neurological deficits of patients with SAO.

METHODS

As neurological deficits, the National Institutes of Health Stroke Scale (NIHSS) was used. Sixty-six SAO participants within 24 h of first onset were enrolled and divided into mild and moderate groups based on NIHSS. QSM values of infarct area and GM volume were calculated from magnetic resonance imaging (MRI) data. Two-sample t-tests were used to compare differences in QSM value and GM volume between the two groups, and the diagnostic efficacy of the combination of QSM value and GM volume was evaluated.

RESULTS

The results revealed both the QSM value and GM volume within the infarct area of the moderate group were lower compared to the mild group. Moderate group exhibited lower GM volume in some specific gyrus compared with mild group in the case of voxel-wise GM volume on whole-brain voxel level. The support vector machine (SVM) classifier's analysis showed a high power for the combination of QSM value, GM volume within the infarct area, and voxel-wise GM volume.

CONCLUSION

Our research first reported the combination of QSM value, GM volume within the infarct area, and voxel-wise GM volume could be used to predict neurological impairment of patients with SAO, which provides new insights for further understanding the SAO stroke.

Unwrapping phase contrast MRI by iterative graph cuts

PURPOSE

To develop and evaluate a phase unwrapping method for cine phase contrast MRI based on graph cuts.

METHODS

A proposed Iterative Graph Cuts method was evaluated in 10 cardiac patients with two-dimensional flow quantification which was repeated at low venc settings to provoke wrapping. The images were also unwrapped by a path-following method (ROMEO), and a Laplacian-based method (LP). Net flow was quantified using semi-automatic vessel segmentation. High venc images were also wrapped retrospectively to asses the residual amount of wrapped voxels.

RESULTS

The absolute net flow error after unwrapping at venc = 100 cm/s was 1.8 mL, which was 0.83 mL smaller than for LP. The repeatability error at high venc without unwrapping was 2.5 mL. The error at venc = 50 cm/s was 7.5 mL, which was 8.2 mL smaller than for ROMEO and 5.7 mL smaller than for LP. For retrospectively wrapped images with synthetic venc of 100/50/25 cm/s, the residual amount of wrapped voxels was 0.00/0.12/0.79%, which was 0.09/0.26/8.0 percentage points smaller than for LP. With synthetic venc of 25 cm/s, omitting magnitude information resulted in 3.2 percentage points more wrapped voxels, and only spatial/temporal unwrapping resulted in 4.6/21 percentage points more wrapped voxels compared to spatiotemporal unwrapping.

CONCLUSION

Iterative Graph Cuts enables unwrapping of cine phase contrast MRI with very small errors, except for at extreme blood velocities, with equal or better performance compared to ROMEO and LP. The use of magnitude information and spatiotemporal unwrapping is recommended.

Differently Distributed Iron Deposition in Mild and Moderate Stroke with Small Artery Occlusion: A Pilot Voxel-Based Susceptibility Mapping Study

BACKGROUND

Small artery occlusion (SAO) is a common ischemic stroke subtype. However, its clinical outcome can be more severe than commonly understood. The severity of SAO can vary, ranging from mild to moderate. Iron deposition has been associated with the development and progression of stroke. However, its specific distribution and relationship with stroke severity in SAO remain unclear. The study's purpose is to investigate the differences in iron deposition between mild stroke with SAO (SAO-MiS) and moderate stroke with SAO (SAO-MoS) through quantitative susceptibility mapping (QSM) and its association with neurological deficits.

METHODS

Sixty-eight SAO participants within 24 hours of first onset were enrolled and separated into SAO-MiS and SAO-MoS according to the National Institutes of Health Stroke Scale (NIHSS) scores. QSM helped calculate the susceptibility maps, reflecting the iron content within the brain. The susceptibility maps were analyzed using voxel-wise statistical analysis to compare the iron deposition between SAO-MiS and SAO-MoS. Then, differentially distributed iron deposition helped differentiate between mild and moderate stroke using support vector machine (SVM) methods.

RESULTS

Compared with SAO-MiS, SAO-MoS depicted elevated iron deposition in the left pallidum, parahippocampal gyrus, and superior frontal gyrus medial region, and is lower in the right superior/middle frontal gyrus and bilateral supplementary motor area. Based on iron deposition, the SVM classifier's analysis revealed a high power to discriminate SAO-MoS from SAO-MiS. In addition, fibrinogen, triglyceride (TG), and total cholesterol (TC) were linked with QSM values in specific brain regions.

CONCLUSIONS

Our study first revealed the brain iron distribution after SAO and differently distributed iron deposition in SAO-MiS and SAO-MoS. The results indicate that iron deposition could play a role in the pathophysiology of SAO and its correlation with stroke severity.

A fully automatic parenchyma extraction method for MRI T2* relaxometry of iron loaded liver in transfusion-dependent patients

PURPOSE

To develop a fully automatic parenchyma extraction method for the T2* relaxometry of iron overload liver.

METHODS

A retrospective multicenter collection of liver MR examinations from 177 transfusion-dependent patients was conducted. The proposed method extended a semiautomatic parenchyma extraction algorithm to a fully automatic approach by introducing a modified TransUNet on the R2* (1/T2*) map for liver segmentation. Axial liver slices from 129 patients at 1.5 T were allocated to training (85%) and internal test (15%) sets. Two external test sets separately included 1.5 T data from 20 patients and 3.0 T data from 28 patients. The final T2* measurement was obtained by fitting the average signal of the extracted liver parenchyma. The agreement between T2* measurements using fully and semiautomatic parenchyma extraction methods was assessed using coefficient of variation (CoV) and Bland-Altman plots.

RESULTS

Dice of the deep network-based liver segmentation was 0.970 ± 0.019 on the internal dataset, 0.960 ± 0.035 on the external 1.5 T dataset, and 0.958 ± 0.014 on the external 3.0 T dataset. The mean difference bias between T2* measurements of the fully and semiautomatic methods were separately 0.12 (95% CI: -0.37, 0.61) ms, 0.04 (95% CI: -1.0, 1.1) ms, and 0.01 (95% CI: -0.25, 0.23) ms on the three test datasets. The CoVs between the two methods were 4.2%, 4.8% and 2.0% on the internal test set and two external test sets.

CONCLUSIONS

The developed fully automatic parenchyma extraction approach provides an efficient and operator-independent T2* measurement for assessing hepatic iron content in clinical practice.

Phase-sensitive deep reconstruction method for rapid multiparametric MR fingerprinting and quantitative susceptibility mapping in the brain

INTRODUCTION

This study explores the potential of Magnetic Resonance Fingerprinting (MRF) with a novel Phase-Sensitivity Deep Reconstruction Network (PS-DRONE) for simultaneous quantification of T1, T2, Proton Density, B1+, phase and quantitative susceptibility mapping (QSM).

METHODS

Data were acquired at 3 T in vitro and in vivo using an optimized EPI-based MRF sequence. Phantom experiments were conducted using a standardized phantom for T1 and T2 maps and a custom-made agar-based gadolinium phantom for B1 and QSM maps. In vivo experiments included five healthy volunteers and one patient diagnosed with brain metastasis. PSDRONE maps were compared to reference maps obtained through standard imaging sequences.

RESULTS

Total scan time was 2 min for 32 slices and a resolution of [1 mm, 1 mm, 4.5 mm]. The reconstruction of T1, T2, Proton Density, B1+ and phase maps were reconstructed within 1 s. In the phantoms, PS-DRONE analysis presented accurate and strongly correlated T1 and T2 maps (r = 0.99) compared to the reference maps. B1 maps from PS-DRONE showed slightly higher values, though still correlated (r = 0.6) with the reference. QSM values showed a small bias but were strongly correlated (r = 0.99) with reference data. In the in vivo analysis, PS-DRONE-derived T1 and T2 values for gray and white matter matched reference values in healthy volunteers. PS-DRONE B1 and QSM maps showed strong correlations with reference values.

CONCLUSION

The PS-DRONE network enables concurrent acquisition of T1, T2, PD, B1+, phase and QSM maps, within 2 min of acquisition time and 1 s of reconstruction time.

Recommended implementation of quantitative susceptibility mapping for clinical research in the brain: A consensus of the ISMRM electro-magnetic tissue properties study group

This article provides recommendations for implementing QSM for clinical brain research. It is a consensus of the International Society of Magnetic Resonance in Medicine, Electro-Magnetic Tissue Properties Study Group. While QSM technical development continues to advance rapidly, the current QSM methods have been demonstrated to be repeatable and reproducible for generating quantitative tissue magnetic susceptibility maps in the brain. However, the many QSM approaches available have generated a need in the neuroimaging community for guidelines on implementation. This article outlines considerations and implementation recommendations for QSM data acquisition, processing, analysis, and publication. We recommend that data be acquired using a monopolar 3D multi-echo gradient echo (GRE) sequence and that phase images be saved and exported in Digital Imaging and Communications in Medicine (DICOM) format and unwrapped using an exact unwrapping approach. Multi-echo images should be combined before background field removal, and a brain mask created using a brain extraction tool with the incorporation of phase-quality-based masking. Background fields within the brain mask should be removed using a technique based on SHARP or PDF, and the optimization approach to dipole inversion should be employed with a sparsity-based regularization. Susceptibility values should be measured relative to a specified reference, including the common reference region of the whole brain as a region of interest in the analysis. The minimum acquisition and processing details required when reporting QSM results are also provided. These recommendations should facilitate clinical QSM research and promote harmonized data acquisition, analysis, and reporting.

Incorporating a-priori information in deep learning models for quantitative susceptibility mapping via adaptive convolution

Quantitative susceptibility mapping (QSM) has attracted considerable interest for tissue characterization (e.g., iron and calcium accumulation, myelination, venous vasculature) in the human brain and relies on extensive data processing of gradient-echo MRI phase images. While deep learning-based field-to-susceptibility inversion has shown great potential, the acquisition parameters applied in clinical settings such as image resolution or image orientation with respect to the magnetic field have not been fully accounted for. Furthermore, the lack of comprehensive training data covering a wide range of acquisition parameters further limits the current QSM deep learning approaches. Here, we propose the integration of a priori information of imaging parameters into convolutional neural networks with our approach, adaptive convolution, that learns the mapping between the additional presented information (acquisition parameters) and the changes in the phase images associated with these varying acquisition parameters. By associating a-priori information with the network parameters itself, the optimal set of convolution weights is selected based on data-specific attributes, leading to generalizability towards changes in acquisition parameters. Moreover, we demonstrate the feasibility of pre-training on synthetic data and transfer learning to clinical brain data to achieve substantial improvements in the computation of susceptibility maps. The adaptive convolution 3D U-Net demonstrated generalizability in acquisition parameters on synthetic and in-vivo data and outperformed models lacking adaptive convolution or transfer learning. Further experiments demonstrate the impact of the side information on the adaptive model and assessed susceptibility map computation on simulated pathologic data sets and measured phase data.

Renal MRI: From Nephron to NMR Signal

Renal diseases pose a significant socio-economic burden on healthcare systems. The development of better diagnostics and prognostics is well-recognized as a key strategy to resolve these challenges. Central to these developments are MRI biomarkers, due to their potential for monitoring of early pathophysiological changes, renal disease progression or treatment effects. The surge in renal MRI involves major cross-domain initiatives, large clinical studies, and educational programs. In parallel with these translational efforts, the need for greater (patho)physiological specificity remains, to enable engagement with clinical nephrologists and increase the associated health impact. The ISMRM 2022 Member Initiated Symposium (MIS) on renal MRI spotlighted this issue with the goal of inspiring more solutions from the ISMRM community. This work is a summary of the MIS presentations devoted to: 1) educating imaging scientists and clinicians on renal (patho)physiology and demands from clinical nephrologists, 2) elucidating the connection of MRI parameters with renal physiology, 3) presenting the current state of leading MR surrogates in assessing renal structure and functions as well as their next generation of innovation, and 4) describing the potential of these imaging markers for providing clinically meaningful renal characterization to guide or supplement clinical decision making. We hope to continue momentum of recent years and introduce new entrants to the development process, connecting (patho)physiology with (bio)physics, and conceiving new clinical applications. We envision this process to benefit from cross-disciplinary collaboration and analogous efforts in other body organs, but also to maximally leverage the unique opportunities of renal physiology. LEVEL OF EVIDENCE: 1 TECHNICAL EFFICACY STAGE: 2.

A new 3D phase unwrapping method by region partitioning and local polynomial modeling in abdominal quantitative susceptibility mapping

Background

Accurate phase unwrapping is a critical prerequisite for successful applications in phase-related MRI, including quantitative susceptibility mapping (QSM) and susceptibility weighted imaging. However, many existing 3D phase unwrapping algorithms face challenges in the presence of severe noise, rapidly changing phase, and open-end cutline.

Methods

In this study, we introduce a novel 3D phase unwrapping approach utilizing region partitioning and a local polynomial model. Initially, the method leverages phase partitioning to create initial regions. Noisy voxels connecting areas within these regions are excluded and grouped into residual voxels. The connected regions within the region of interest are then reidentified and categorized into blocks and residual voxels based on voxel count thresholds. Subsequently, the method sequentially performs inter-block and residual voxel phase unwrapping using the local polynomial model. The proposed method was evaluated on simulation and in vivo abdominal QSM data, and was compared with the classical Region-growing, Laplacian_based, Graph-cut, and PRELUDE methods.

Results

Simulation experiments, conducted under different signal-to-noise ratios and phase change levels, consistently demonstrate that the proposed method achieves accurate unwrapping results, with mean error ratios not exceeding 0.01%. In contrast, the error ratios of Region-growing (N/A, 84.47%), Laplacian_based (20.65%, N/A), Graph-cut (2.26%, 20.71%), and PRELUDE (4.28%, 10.33%) methods are all substantially higher than those of the proposed method. In vivo abdominal QSM experiments further confirm the effectiveness of the proposed method in unwrapping phase data and successfully reconstructing susceptibility maps, even in scenarios with significant noise, rapidly changing phase, and open-end cutline in a large field of view.

Conclusion

The proposed method demonstrates robust and accurate phase unwrapping capabilities, positioning it as a promising option for abdominal QSM applications.

BFRnet: A deep learning-based MR background field removal method for QSM of the brain containing significant pathological susceptibility sources

INTRODUCTION

Background field removal (BFR) is a critical step required for successful quantitative susceptibility mapping (QSM). However, eliminating the background field in brains containing significant susceptibility sources, such as intracranial hemorrhages, is challenging due to the relatively large scale of the field induced by these pathological susceptibility sources.

METHOD

This study proposes a new deep learning-based method, BFRnet, to remove the background field in healthy and hemorrhagic subjects. The network is built with the dual-frequency octave convolutions on the U-net architecture, trained with synthetic field maps containing significant susceptibility sources. The BFRnet method is compared with three conventional BFR methods and one previous deep learning method using simulated and in vivo brains from 4 healthy and 2 hemorrhagic subjects. Robustness against acquisition field-of-view (FOV) orientation and brain masking are also investigated.

RESULTS

For both simulation and in vivo experiments, BFRnet led to the best visually appealing results in the local field and QSM results with the minimum contrast loss and the most accurate hemorrhage susceptibility measurements among all five methods. In addition, BFRnet produced the most consistent local field and susceptibility maps between different sizes of brain masks, while conventional methods depend drastically on precise brain extraction and further brain edge erosions. It is also observed that BFRnet performed the best among all BFR methods for acquisition FOVs oblique to the main magnetic field.

CONCLUSION

The proposed BFRnet improved the accuracy of local field reconstruction in the hemorrhagic subjects compared with conventional BFR algorithms. The BFRnet method was effective for acquisitions of tilted orientations and retained whole brains without edge erosion as often required by traditional BFR methods.

Improved dynamic distortion correction for fMRI using single-echo EPI and a readout-reversed first image (REFILL)

The boundaries between tissues with different magnetic susceptibilities generate inhomogeneities in the main magnetic field which change over time due to motion, respiration and system instabilities. The dynamically changing field can be measured from the phase of the fMRI data and corrected. However, methods for doing so need multi-echo data, time-consuming reference scans and/or involve error-prone processing steps, such as phase unwrapping, which are difficult to implement robustly on the MRI host. The improved dynamic distortion correction method we propose is based on the phase of the single-echo EPI data acquired for fMRI, phase offsets calculated from a triple-echo, bipolar reference scan of circa 3-10 s duration using a method which avoids the need for phase unwrapping and an additional correction derived from one EPI volume in which the readout direction is reversed. This Reverse-Encoded First Image and Low resoLution reference scan (REFILL) approach is shown to accurately measure B0 as it changes due to shim, motion and respiration, even with large dynamic changes to the field at 7 T, where it led to a > 20% increase in time-series signal to noise ratio compared to data corrected with the classic static approach. fMRI results from REFILL-corrected data were free of stimulus-correlated distortion artefacts seen when data were corrected with static field mapping. The method is insensitive to shim changes and eddy current differences between the reference scan and the fMRI time series, and employs calculation steps that are simple and robust, allowing most data processing to be performed in real time on the scanner image reconstruction computer. These improvements make it feasible to routinely perform dynamic distortion correction in fMRI.

Ferric carboxymaltose effects on restless legs syndrome and on brain iron in patients with iron deficiency anemia

OBJECTIVE

Brain iron status is fundamental in RLS pathogenesis. The aim of this study was to determine the clinical efficacy and brain iron concentration improvement in RLS patients with IDA, using 1500 mg FCM.

METHODS

This is a randomized, double-blinded, placebo-controlled study. RLS patients with IDA were grouped into either 1500 mg FCM or placebo. The primary outcomes were the change from baseline on the International Restless Legs Syndrome Study Group scale (IRLS) and brain iron measured by QSM and R2∗.

RESULTS

A total of 18 RLS patients with IDA were enrolled, 10 in the FCM group and 8 in the placebo. At the week 6 endpoint, the FCM group showed significant improvement in both IRLS (-13.60 ± 9.47 vs. -3.63 ± 5.40, p = 0.011) and VAS (-40.50 ± 28.81 vs. -0.63 ± 28.28, p = 0.004) from baseline. Change from baseline with R2∗ techniques showed a treatment effect for the thalamus and QSM technique for both the substantia nigra and pulvinar. A correlation was proved between the IRLS difference and the difference of QSM in thalamus (p = 0.028).

CONCLUSION

This study demonstrates that 1500 mg FCM effectively treats RLS symptoms in IDA patients over six weeks, with MRI measurements of improved brain iron content serving as a potential biomarker for RLS patients.

High angular resolution susceptibility imaging and estimation of fiber orientation distribution functions in primate brain

Uncovering brain-tissue microstructure including axonal characteristics is a major neuroimaging research focus. Within this scope, anisotropic properties of magnetic susceptibility in white matter have been successfully employed to estimate primary axonal trajectories using mono-tensorial models. However, anisotropic susceptibility has not yet been considered for modeling more complex fiber structures within a voxel, such as intersecting bundles, or an estimation of orientation distribution functions (ODFs). This information is routinely obtained by high angular resolution diffusion imaging (HARDI) techniques. In applications to fixed tissue, however, diffusion-weighted imaging suffers from an inherently low signal-to-noise ratio and limited spatial resolution, leading to high demands on the performance of the gradient system in order to mitigate these limitations. In the current work, high angular resolution susceptibility imaging (HARSI) is proposed as a novel, phase-based methodology to estimate ODFs. A multiple gradient-echo dataset was acquired in an entire fixed chimpanzee brain at 61 orientations by reorienting the specimen in the magnetic field. The constant solid angle method was adapted for estimating phase-based ODFs. HARDI data were also acquired for comparison. HARSI yielded information on whole-brain fiber architecture, including identification of peaks of multiple bundles that resembled features of the HARDI results. Distinct differences between both methods suggest that susceptibility properties may offer complementary microstructural information. These proof-of-concept results indicate a potential to study the axonal organization in post-mortem primate and human brain at high resolution.

Recommended Implementation of Quantitative Susceptibility Mapping for Clinical Research in The Brain: A Consensus of the ISMRM Electro-Magnetic Tissue Properties Study Group

This article provides recommendations for implementing quantitative susceptibility mapping (QSM) for clinical brain research. It is a consensus of the ISMRM Electro-Magnetic Tissue Properties Study Group. While QSM technical development continues to advance rapidly, the current QSM methods have been demonstrated to be repeatable and reproducible for generating quantitative tissue magnetic susceptibility maps in the brain. However, the many QSM approaches available give rise to the need in the neuroimaging community for guidelines on implementation. This article describes relevant considerations and provides specific implementation recommendations for all steps in QSM data acquisition, processing, analysis, and presentation in scientific publications. We recommend that data be acquired using a monopolar 3D multi-echo GRE sequence, that phase images be saved and exported in DICOM format and unwrapped using an exact unwrapping approach. Multi-echo images should be combined before background removal, and a brain mask created using a brain extraction tool with the incorporation of phase-quality-based masking. Background fields should be removed within the brain mask using a technique based on SHARP or PDF, and the optimization approach to dipole inversion should be employed with a sparsity-based regularization. Susceptibility values should be measured relative to a specified reference, including the common reference region of whole brain as a region of interest in the analysis, and QSM results should be reported with - as a minimum - the acquisition and processing specifications listed in the last section of the article. These recommendations should facilitate clinical QSM research and lead to increased harmonization in data acquisition, analysis, and reporting.

Evaluation of multi-channel phase reconstruction methods for quantitative susceptibility mapping on postmortem human brain

Quantitative Susceptibility Mapping (QSM) is an established Magnetic Resonance Imaging (MRI) technique with high potential in brain iron studies associated to several neurodegenerative diseases. Unlike other MRI techniques, QSM relies on phase images to estimate tissue's relative susceptibility, therefore requiring a reliable phase data. Phase images from a multi-channel acquisition should be reconstructed in a proper way. On this work it was compared the performance of combination of phase matching algorithms (MCPC3D-S and VRC) and phase combination methods based on a complex weighted sum of phases, considering the magnitude at different powers (k = 0 to 4) as the weighting factor. These reconstruction methods were applied in two datasets: a simulated brain dataset for a 4-coil array and data of 22 postmortem subjects acquired at a 7T scanner using a 32 channels coil. For the simulated dataset, differences between the ground truth and the Root Mean Squared Error (RMSE) were evaluated. For both simulated and postmortem data, the mean (MS) and standard deviation (SD) of susceptibility values of five deep gray matter regions were calculated. For the postmortem subjects, MS and SD were statistically compared across all subjects. A qualitative analysis indicated no differences between methods, except for the Adaptive approach on postmortem data, which showed intense artifacts. In the 20% noise level case, the simulated data showed increased noise in central regions. Quantitative analysis showed that both MS and SD were not statistically different when comparing k = 1 and k = 2 on postmortem brain images, however visual inspection showed some boundaries artifacts on k = 2. Furthermore, the RMSE decreased (on regions near the coils) and increased (on central regions and on overall QSM) with increasing k. In conclusion, for reconstruction of phase images from multiple coils with no reference available, alternative methods are needed. In this study it was found that overall, the phase combination with k = 1 is preferred over other powers of k.

Tilted quantitative susceptibility mapping at oblique MRI (tiltQSM)

OBJECTIVE

If the phase image matrix was acquired from oblique MRI, it is needed to deal with the oblique effect for quantitative susceptibility mapping (QSM), as addressed in this paper.

METHODS

We proposed two methods for QSM reconstruction from slice-tilted MRI phase image (tiltQSM): 1) rotData per anti-tilting phase image rotation back into the B₀-upright system, and 2) rotKernel per pro-tilting dipole kernel rotation into the same oblique setting as defined by the tilted phase image. Both matrix methods were implemented in an additional preprocessing subroutine to ensure that the phase image and the dipole kernel were represented in the same coordinate system (either in B₀-upright system or in B₀-tilted system); thereafter tiltQSM could be completed through a regular QSM procedure. Besides the oblique effect, tiltQSM also suffers from MRI anisotropy. We provided numeric simulations, phantom tests and in vivo brain experiments on tiltQSM with oblique MRI (axial slice tilting at 3T).

RESULTS

The tiltQSM reconstruction could attain a performance corr > 0.90 (spatial correlation conformance) for small tilting angles <10°. The tiltQSM performance could be further degraded by voxel anisotropy due to image matrix rotation (digital geometry error).

CONCLUSIONS

To seek inverse solutions of MRI phase images acquired at oblique MRI (e.g. in axial slice tilting), we proposed tiltQSM to deal with the oblique effect per matrix rotation (either rotData or rotKernel) in a preprocessing subroutine prior to a regular QSM procedure. In practice, it is always recommended to acquire MRI phase images in isotropic matrix at zero obliqueness (or limited to small tilting angles <10°) for maximal (optimal) QSM reconstruction.

A method for measuring B0 field inhomogeneity using quantitative double-echo in steady-state

PURPOSE

To develop and validate a method forB 0 $$ {B}_0 $$ mapping for knee imaging using the quantitative Double-Echo in Steady-State (qDESS) exploiting the phase difference (Δ θ $$ \Delta \theta $$ ) between the two echoes acquired. Contrary to a two-gradient-echo (2-GRE) method,Δ θ $$ \Delta \theta $$ depends only on the first echo time.

METHODS

Bloch simulations were applied to investigate robustness to noise of the proposed methodology and all imaging studies were validated with phantoms and in vivo simultaneous bilateral knee acquisitions. Two phantoms and five healthy subjects were scanned using qDESS, water saturation shift referencing (WASSR), and multi-GRE sequences.Δ B 0 $$ \Delta {B}_0 $$ maps were calculated with the qDESS and the 2-GRE methods and compared against those obtained with WASSR. The comparison was quantitatively assessed exploiting pixel-wise difference maps, Bland-Altman (BA) analysis, and Lin's concordance coefficient (ρ c $$ {\rho}_c $$ ). For in vivo subjects, the comparison was assessed in cartilage using average values in six subregions.

RESULTS

The proposed method for measuringΔ B 0 $$ \Delta {B}_0 $$ inhomogeneities from a qDESS acquisition providedΔ B 0 $$ \Delta {B}_0 $$ maps that were in good agreement with those obtained using WASSR.Δ B 0 $$ \Delta {B}_0 $$ ρ c $$ {\rho}_c $$ values were≥ $$ \ge $$ 0.98 and 0.90 in phantoms and in vivo, respectively. The agreement between qDESS and WASSR was comparable to that of a 2-GRE method.

CONCLUSION

The proposed method may allow B0 correction for qDESST 2 $$ {T}_2 $$ mapping using an inherently co-registeredΔ B 0 $$ \Delta {B}_0 $$ map without requiring an additional B0 measurement sequence. More generally, the method may help shorten knee imaging protocols that require an auxiliaryΔ B 0 $$ \Delta {B}_0 $$ map by exploiting a qDESS acquisition that also providesT 2 $$ {T}_2 $$ measurements and high-quality morphological imaging.

Comparison of quantitative susceptibility mapping methods for iron-sensitive susceptibility imaging at 7T: An evaluation in healthy subjects and patients with Huntington's disease

Quantitative susceptibility mapping (QSM) is a promising tool for investigating iron dysregulation in neurodegenerative diseases, including Huntington's disease (HD). Many diverse methods have been proposed to generate accurate and robust QSM images. In this study, we evaluated the performance of different dipole inversion algorithms for iron-sensitive susceptibility imaging at 7T on healthy subjects of a large age range and patients with HD. We compared an iterative least-squares-based method (iLSQR), iterative methods that use regularization, single-step approaches, and deep learning-based techniques. Their performance was evaluated by comparing: (1) deviations from a multiple-orientation QSM reference; (2) visual appearance of QSM maps and the presence of artifacts; (3) susceptibility in subcortical brain regions with age; (4) regional brain susceptibility with published postmortem brain iron quantification; and (5) susceptibility in HD-affected basal ganglia regions between HD subjects and healthy controls. We found that single-step QSM methods with either total variation or total generalized variation constraints (SSTV/SSTGV) and the single-step deep learning method iQSM generally provided the best performance in terms of correlation with iron deposition and were better at differentiating between healthy controls and premanifest HD individuals, while deep learning QSM methods trained with multiple-orientation susceptibility data created QSM maps that were most similar to the multiple orientation reference and with the best visual scores.

Unsupervised physiological noise correction of functional magnetic resonance imaging data using phase and magnitude information (PREPAIR)

Of the sources of noise affecting blood oxygen level-dependent functional magnetic resonance imaging (fMRI), respiration and cardiac fluctuations are responsible for the largest part of the variance, particularly at high and ultrahigh field. Existing approaches to removing physiological noise either use external recordings, which can be unwieldy and unreliable, or attempt to identify physiological noise from the magnitude fMRI data. Data-driven approaches are limited by sensitivity, temporal aliasing, and the need for user interaction. In the light of the sensitivity of the phase of the MR signal to local changes in the field stemming from physiological processes, we have developed an unsupervised physiological noise correction method using the information carried in the phase and the magnitude of echo-planar imaging data. Our technique, Physiological Regressor Estimation from Phase and mAgnItude, sub-tR (PREPAIR) derives time series signals sampled at the slice TR from both phase and magnitude images. It allows physiological noise to be captured without aliasing, and efficiently removes other sources of signal fluctuations not related to physiology, prior to regressor estimation. We demonstrate that the physiological signal time courses identified with PREPAIR agree well with those from external devices and retrieve challenging cardiac dynamics. The removal of physiological noise was as effective as that achieved with the most used approach based on external recordings, RETROICOR. In comparison with widely used recording-free physiological noise correction tools-PESTICA and FIX, both performed in unsupervised mode-PREPAIR removed significantly more respiratory and cardiac noise than PESTICA, and achieved a larger increase in temporal signal-to-noise-ratio at both 3 and 7 T.

Phase-based masking for quantitative susceptibility mapping of the human brain at 9.4T

PURPOSE

To develop improved tissue masks for QSM.

METHODS

Masks including voxels at the brain surface were automatically generated from the magnitude alone (MM) or combined with test functions from the first (PG) or second (PB) derivative of the sign of the wrapped phase. Phase images at 3T and 9.4T were simulated at different TEs and used to generate a mask, PItoh , with between-voxel phase differences less than π. MM, PG, and PB were compared with PItoh . QSM were generated from 3D multi-echo gradient-echo data acquired at 9.4T (21 subjects aged: 20-56y), and from the QSM2016 challenge 3T data using different masks, unwrapping, background removal, and dipole inversion algorithms. QSM contrast was quantified using age-based iron concentrations.

RESULTS

Close to air cavities, phase wraps became denser with increasing field and echo time, yielding increased values of the test functions. Compared with PItoh , PB had the highest Dice coefficient, while PG had the lowest and MM the highest percentage of voxels outside PItoh. Artifacts observed in QSM at 9.4T with MM were mitigated by stronger background filters but yielded a reduced QSM contrast. With PB, QSM contrast was greater and artifacts diminished. Similar results were obtained with challenge data, evidencing larger effects of mask close to air cavities.

CONCLUSION

Automatic, phase-based masking founded on the second derivative of the sign of the wrapped phase, including cortical voxels at the brain surface, was able to mitigate artifacts and restore QSM contrast across cortical and subcortical brain regions.

Quantitative susceptibility mapping (QSM) and R2* of silent cerebral infarcts in sickle cell anemia

Silent cerebral infarction (SCI) is the most commonly reported radiological abnormality in patients with sickle cell anemia (SCA) and is associated with future clinical stroke risk. To date, there have been few histological and quantitative MRI studies of SCI and multiple radiological definitions exist. As a result, the tissue characteristics and composition of SCI remain elusive. The objective of this work was therefore to investigate the composition of segmented SCI lesions using quantitative MRI for R2 * and quantitative magnetic susceptibility mapping (QSM). 211 SCI lesions were segmented from 32 participants with SCA and 6 controls. SCI were segmented according to two definitions (FLAIR+/-T1w-based threshold) using a semi-automated pipeline. Magnetic susceptibility (χ) and R2 * maps were calculated from a multi-echo gradient echo sequence and mean SCI values were compared to an equivalent region of interest in normal appearing white matter (NAWM). SCI χ and R2 * were investigated as a function of SCI definition, patient demographics, anatomical location, and cognition. Compared to NAWM, SCI were significantly less diamagnetic (χ = -0.0067 ppm vs. -0.0153 ppm, p < 0.001) and had significantly lower R2 * (16.7 s-1 vs. 19.2 s-1, p < 0.001). SCI definition had a significant effect on the mean SCI χ and R2 * , with lesions becoming significantly less diamagnetic and having significantly lower R2 * after the application of a more stringent T1w-based threshold. SCI-NAWM R2 * decrease was significantly greater in patients with SCA compared with controls (-2.84 s-1 vs. -0.64 s-1, p < 0.0001). No significant association was observed between mean SCI-NAWM χ or R2* differences and subject age, lesion anatomical location, or cognition. The increased χ and decreased R2 * in SCI relative to NAWM observed in both patients and controls is indicative of lower myelin or increased water content within the segmented lesions. The significant SCI-NAWM R2 * differences observed between SCI in patients with SCA and controls suggests there may be differences in tissue composition relative to NAWM in SCI in the two populations. Quantitative MRI techniques such as QSM and R2 * mapping can be used to enhance our understanding of the pathophysiology and composition of SCI in patients with SCA as well as controls.

Quantitative susceptibility mapping as an imaging biomarker for Alzheimer’s disease: The expectations and limitations

Alzheimer’s disease (AD) is the most common type of dementia and a distressing diagnosis for individuals and caregivers. Researchers and clinical trials have mainly focused on β-amyloid plaques, which are hypothesized to be one of the most important factors for neurodegeneration in AD. Meanwhile, recent clinicopathological and radiological studies have shown closer associations of tau pathology rather than β-amyloid pathology with the onset and progression of Alzheimer’s symptoms. Toward a biological definition of biomarker-based research framework for AD, the 2018 National Institute on Aging–Alzheimer’s Association working group has updated the ATN classification system for stratifying disease status in accordance with relevant pathological biomarker profiles, such as cerebral β-amyloid deposition, hyperphosphorylated tau, and neurodegeneration. In addition, altered iron metabolism has been considered to interact with abnormal proteins related to AD pathology thorough generating oxidative stress, as some prior histochemical and histopathological studies supported this iron-mediated pathomechanism. Quantitative susceptibility mapping (QSM) has recently become more popular as a non-invasive magnetic resonance technique to quantify local tissue susceptibility with high spatial resolution, which is sensitive to the presence of iron. The association of cerebral susceptibility values with other pathological biomarkers for AD has been investigated using various QSM techniques; however, direct evidence of these associations remains elusive. In this review, we first briefly describe the principles of QSM. Second, we focus on a large variety of QSM applications, ranging from common applications, such as cerebral iron deposition, to more recent applications, such as the assessment of impaired myelination, quantification of venous oxygen saturation, and measurement of blood– brain barrier function in clinical settings for AD. Third, we mention the relationships among QSM, established biomarkers, and cognitive performance in AD. Finally, we discuss the role of QSM as an imaging biomarker as well as the expectations and limitations of clinically useful diagnostic and therapeutic implications for AD.

Multi-echo quantitative susceptibility mapping: how to combine echoes for accuracy and precision at 3 Tesla

Purpose

To compare different multi‐echo combination methods for MRI QSM. Given the current lack of consensus, we aimed to elucidate how to optimally combine multi‐echo gradient‐recalled echo signal phase information, either before or after applying Laplacian‐base methods (LBMs) for phase unwrapping or background field removal.

Methods

Multi‐echo gradient‐recalled echo data were simulated in a numerical head phantom, and multi‐echo gradient‐recalled echo images were acquired at 3 Tesla in 10 healthy volunteers. To enable image‐based estimation of gradient‐recalled echo signal noise, 5 volunteers were scanned twice in the same session without repositioning. Five QSM processing pipelines were designed: 1 applied nonlinear phase fitting over TEs before LBMs; 2 applied LBMs to the TE‐dependent phase and then combined multiple TEs via either TE‐weighted or SNR‐weighted averaging; and 2 calculated TE‐dependent susceptibility maps via either multi‐step or single‐step QSM and then combined multiple TEs via magnitude‐weighted averaging. Results from different pipelines were compared using visual inspection; summary statistics of susceptibility in deep gray matter, white matter, and venous regions; phase noise maps (error propagation theory); and, in the healthy volunteers, regional fixed bias analysis (Bland–Altman) and regional differences between the means (nonparametric tests).

Results

Nonlinearly fitting the multi‐echo phase over TEs before applying LBMs provided the highest regional accuracy of χ and the lowest phase noise propagation compared to averaging the LBM‐processed TE‐dependent phase. This result was especially pertinent in high‐susceptibility venous regions.

Conclusion

For multi‐echo QSM, we recommend combining the signal phase by nonlinear fitting before applying LBMs.

Click here for author‐reader discussions

Iterative static field map estimation for off-resonance correction in non-Cartesian susceptibility weighted imaging

Purpose

Patient‐induced inhomogeneities in the magnetic field cause distortions and blurring during acquisitions with long readouts such as in susceptibility‐weighted imaging (SWI). Most correction methods require collecting an additional ΔB0 field map to remove these artifacts.

Theory

The static ΔB0 field map can be approximated with an acceptable error directly from a single echo acquisition in SWI. The main component of the observed phase is linearly related to ΔB0 and the echo time (TE), and the relative impact of non‐ ΔB0 terms becomes insignificant with TE >20 ms at 3 T for a well‐tuned system.

Methods

The main step is to combine and unfold the multi‐channel phase maps wrapped many times, and several competing algorithms are compared for this purpose. Four in vivo brain data sets collected using the recently proposed 3D spreading projection algorithm for rapid k‐space sampling (SPARKLING) readouts are used to assess the proposed method.

Results

The estimated 3D field maps generated with a 0.6 mm isotropic spatial resolution provide overall similar off‐resonance corrections compared to reference corrections based on an external ΔB0 acquisitions, and even improved for 2 of 4 individuals. Although a small estimation error is expected, no aftermath was observed in the proposed corrections, whereas degradations were observed in the references.

Conclusion

A static ΔB0 field map estimation method was proposed to take advantage of acquisitions with long echo times, and outperformed the reference technique based on an external field map. The difference can be attributed to an inherent robustness to mismatches between volumes and external ΔB0 maps, and diverse other sources investigated.

Click here for author‐reader discussions

Instant tissue field and magnetic susceptibility mapping from MRI raw phase using Laplacian enhanced deep neural networks

Quantitative susceptibility mapping (QSM) is an MRI post-processing technique that produces spatially resolved magnetic susceptibility maps from phase data. However, the traditional QSM reconstruction pipeline involves multiple non-trivial steps, including phase unwrapping, background field removal, and dipole inversion. These intermediate steps not only increase the reconstruction time but accumulates errors. This study aims to overcome existing limitations by developing a Laplacian-of-Trigonometric-functions (LoT) enhanced deep neural network for near-instant quantitative field and susceptibility mapping (i.e., iQFM and iQSM) from raw MRI phase data. The proposed iQFM and iQSM methods were compared with established reconstruction pipelines on simulated and in vivo datasets. In addition, experiments on patients with intracranial hemorrhage and multiple sclerosis were also performed to test the generalization of the proposed neural networks. The proposed iQFM and iQSM methods in healthy subjects yielded comparable results to those involving the intermediate steps while dramatically improving reconstruction accuracies on intracranial hemorrhages with large susceptibilities. High susceptibility contrast between multiple sclerosis lesions and healthy tissue was also achieved using the proposed methods. Comparative studies indicated that the most significant contributor to iQFM and iQSM over conventional multi-step methods was the elimination of traditional Laplacian unwrapping. The reconstruction time on the order of minutes for traditional approaches was shortened to around 0.1 seconds using the trained iQFM and iQSM neural networks.

Hybrid data fidelity term approach for quantitative susceptibility mapping

Purpose

Methods

Results

Conclusions

Overview of quantitative susceptibility mapping using deep learning: Current status, challenges and opportunities

Quantitative susceptibility mapping (QSM) has gained broad interest in the field by extracting bulk tissue magnetic susceptibility, predominantly determined by myelin, iron and calcium from magnetic resonance imaging (MRI) phase measurements in vivo. Thereby, QSM can reveal pathological changes of these key components in a variety of diseases. QSM requires multiple processing steps such as phase unwrapping, background field removal and field-to-source inversion. Current state-of-the-art techniques utilize iterative optimization procedures to solve the inversion and background field correction, which are computationally expensive and require a careful choice of regularization parameters. With the recent success of deep learning using convolutional neural networks for solving ill-posed reconstruction problems, the QSM community also adapted these techniques and demonstrated that the QSM processing steps can be solved by efficient feed forward multiplications not requiring either iterative optimization or the choice of regularization parameters. Here, we review the current status of deep learning-based approaches for processing QSM, highlighting limitations and potential pitfalls, and discuss the future directions the field may take to exploit the latest advances in deep learning for QSM.

Imaging white matter microstructure with gradient-echo phase imaging: Is ex vivo imaging with formalin-fixed tissue a good approximation of the in vivo brain?

Purpose

Ex vivo imaging is a commonly used approach to investigate the biophysical mechanism of orientation‐dependent signal phase evolution in white matter. Yet, how phase measurements are influenced by the structural alteration in the tissue after formalin fixation is not fully understood. Here, we study the effects on magnetic susceptibility, microstructural compartmentalization, and chemical exchange measurement with a postmortem formalin‐fixed whole‐brain human tissue.

Methods

A formalin‐fixed, postmortem human brain specimen was scanned with multiple orientations to the main magnetic field direction for robust bulk magnetic susceptibility measurement with conventional quantitative susceptibility imaging models. White matter samples were subsequently excised from the whole‐brain specimen and scanned in multiple rotations on an MRI scanner to measure the anisotropic magnetic susceptibility and microstructure‐related contributions in the signal phase and to validate the findings of the whole‐brain data.

Results

The bulk isotropic magnetic susceptibility of ex vivo whole‐brain imaging is comparable to in vivo imaging, with noticeable enhanced nonsusceptibility contributions. The excised specimen experiment reveals that anisotropic magnetic susceptibility and compartmentalization phase effect were considerably reduced in the formalin‐fixed white matter specimens.

Conclusions

Formalin‐fixed postmortem white matter exhibits comparable isotropic magnetic susceptibility to previous in vivo imaging findings. However, the measured phase and magnitude data of the fixed white matter tissue shows a significantly weaker orientation dependency and compartmentalization effect. Alternatives to formalin fixation are needed to better reproduce the in vivo microstructural effects in postmortem samples.

Comparison of Phase Estimation Methods for Quantitative Susceptibility Mapping Using a Rotating-Tube Phantom

Quantitative Susceptibility Mapping (QSM) is an MRI tool with the potential to reveal pathological changes from magnetic susceptibility measurements. Before phase data can be used to recover susceptibility (Δχ), the QSM process begins with two steps: data acquisition and phase estimation. We assess the performance of these steps, when applied without user intervention, on several variations of a phantom imaging task. We used a rotating-tube phantom with five tubes ranging from Δχ=0.05 ppm to Δχ=0.336 ppm. MRI data was acquired at nine angles of rotation for four different pulse sequences. The images were processed by 10 phase estimation algorithms including Laplacian, region-growing, branch-cut, temporal unwrapping, and maximum-likelihood methods, resulting in approximately 90 different combinations of data acquisition and phase estimation methods. We analyzed errors between measured and expected phases using the probability mass function and Cumulative Distribution Function. Repeatable acquisition and estimation methods were identified based on the probability of relative phase errors. For single-echo GRE and segmented EPI sequences, a region-growing method was most reliable with Pr (relative error <0.1) = 0.95 and 0.90, respectively. For multiecho sequences, a maximum-likelihood method was most reliable with Pr (relative error <0.1) = 0.97. The most repeatable multiecho methods outperformed the most repeatable single-echo methods. We found a wide range of repeatability and reproducibility for off-the-shelf MRI acquisition and phase estimation approaches, and this variability may prevent the techniques from being widely integrated in clinical workflows. The error was dominated in many cases by spatially discontinuous phase unwrapping errors. Any postprocessing applied on erroneous phase estimates, such as QSM's background field removal and dipole inversion, would suffer from error propagation. Our paradigm identifies methods that yield consistent and accurate phase estimates that would ultimately yield consistent and accurate Δχ estimates.

QSMxT: Robust masking and artifact reduction for quantitative susceptibility mapping

PURPOSE

Quantitative susceptibility mapping (QSM) estimates the spatial distribution of tissue magnetic susceptibilities from the phase of a gradient-echo signal. QSM algorithms require a signal mask to delineate regions with reliable phase for subsequent susceptibility estimation. Existing masking techniques used in QSM have limitations that introduce artifacts, exclude anatomical detail, and rely on parameter tuning and anatomical priors that narrow their application. Here, a robust masking and reconstruction procedure is presented to overcome these limitations and enable automated QSM processing. Moreover, this method is integrated within an open-source software framework: QSMxT.

METHODS

A robust masking technique that automatically separates reliable from less reliable phase regions was developed and combined with a two-pass reconstruction procedure that operates on the separated sources before combination, extracting more information and suppressing streaking artifacts.

RESULTS

Compared with standard masking and reconstruction procedures, the two-pass inversion reduces streaking artifacts caused by unreliable phase and high dynamic ranges of susceptibility sources. It is also robust across a range of acquisitions at 3 T in volunteers and phantoms, at 7 T in tumor patients, and in an in silico head phantom, with significant artifact and error reductions, greater anatomical detail, and minimal parameter tuning.

CONCLUSION

The two-pass masking and reconstruction procedure separates reliable from less reliable phase regions, enabling a more accurate QSM reconstruction that mitigates artifacts, operates without anatomical priors, and requires minimal parameter tuning. The technique and its integration within QSMxT makes QSM processing more accessible and robust to streaking artifacts.

QSM reconstruction challenge 2.0: Design and report of results

PURPOSE

The aim of the second quantitative susceptibility mapping (QSM) reconstruction challenge (Oct 2019, Seoul, Korea) was to test the accuracy of QSM dipole inversion algorithms in simulated brain data.

METHODS

A two-stage design was chosen for this challenge. The participants were provided with datasets of multi-echo gradient echo images synthesized from two realistic in silico head phantoms using an MR simulator. At the first stage, participants optimized QSM reconstructions without ground truth data available to mimic the clinical setting. At the second stage, ground truth data were provided for parameter optimization. Submissions were evaluated using eight numerical metrics and visual ratings.

RESULTS

A total of 98 reconstructions were submitted for stage 1 and 47 submissions for stage 2. Iterative methods had the best quantitative metric scores, followed by deep learning and direct inversion methods. Priors derived from magnitude data improved the metric scores. Algorithms based on iterative approaches and total variation (and its derivatives) produced the best overall results. The reported results and analysis pipelines have been made public to allow researchers to compare new methods to the current state of the art.

CONCLUSION

The synthetic data provide a consistent framework to test the accuracy and robustness of QSM algorithms in the presence of noise, calcifications and minor voxel dephasing effects. Total Variation-based algorithms produced the best results among all metrics. Future QSM challenges should assess whether this good performance with synthetic datasets translates to more realistic scenarios, where background fields and dipole-incompatible phase contributions are included.

Receiver phase alignment using fitted SVD derived sensitivities from routine prescans

Magnetic resonance imaging radio frequency arrays are composed of multiple receive coils that have their signals combined to form an image. Combination requires an estimate of the radio frequency coil sensitivities to align signal phases and prevent destructive interference. At lower fields this can be accomplished using a uniform physical reference coil. However, at higher fields, uniform volume coils are lacking and, when available, suffer from regions of low receive sensitivity that result in poor sensitivity estimation and combination. Several approaches exist that do not require a physical reference coil but require manual intervention, specific prescans, or must be completed post-acquisition. This makes these methods impractical for large multi-volume datasets such as those collected for novel types of functional MRI or quantitative susceptibility mapping, where magnitude and phase are important. This pilot study proposes a fitted SVD method which utilizes existing combination methods to create a phase sensitive combination method targeted at large multi-volume datasets. This method uses any multi-image prescan to calculate the relative receive sensitivities using voxel-wise singular value decomposition. These relative sensitivities are fitted to the solid harmonics using an iterative least squares fitting algorithm. Fits of the relative sensitivities are used to align the phases of the receive coils and improve combination in subsequent acquisitions during the imaging session. This method is compared against existing approaches in the human brain at 7 Tesla by examining the combined data for the presence of singularities and changes in phase signal-to-noise ratio. Two additional applications of the method are also explored, using the fitted SVD method in an asymmetrical coil and in a case with subject motion. The fitted SVD method produces singularity-free images and recovers between 95-100% of the phase signal-to-noise ratio depending on the prescan data resolution. Using solid harmonic fitting to interpolate singular value decomposition derived receive sensitivities from existing prescans allows the fitted SVD method to be used on all acquisitions within a session without increasing exam duration. Our fitted SVD method is able to combine imaging datasets accurately without supervision during online reconstruction.

Improved susceptibility weighted imaging at ultra-high field using bipolar multi-echo acquisition and optimized image processing: CLEAR-SWI

Purpose

Susceptibility Weighted Imaging (SWI) has become established in the clinical investigation of stroke, microbleeds, tumor vascularization, calcification and iron deposition, but suffers from a number of shortcomings and artefacts. The goal of this study was to reduce the sensitivity of SWI to strong B₁ and B₀ inhomogeneities at ultra-high field to generate homogeneous images with increased contrast and free of common artefacts. All steps in SWI processing have been addressed −coil combination, phase unwrapping, image combination over echoes, phase filtering and homogeneity correction −and applied to an efficient bipolar multi-echo acquisition to substantially improve the quality of SWI.

Principal results

Our findings regarding the optimal individual processing steps lead us to propose a Contrast-weighted, Laplace-unwrapped, bipolar multi-Echo, ASPIRE-combined, homogeneous, improved Resolution SWI, or CLEAR-SWI. CLEAR-SWI was compared to two other multi-echo SWI methods and standard, single-echo SWI with the same acquisition time at 7 T in 10 healthy volunteers and with single-echo SWI in 13 patients with brain tumors. CLEAR-SWI had improved contrast-to-noise and homogeneity, reduced signal dropout and was not compromised by the artefacts which affected standard SWI in 10 out of 13 cases close to tumors (as assessed by expert raters), as well as generating T2* maps and phase images which can be used for Quantitative Susceptibility Mapping. In a comparison with other multi-echo SWI methods, CLEAR-SWI had the fewest artefacts, highest SNR and generally higher contrast-to-noise.

Major conclusions

CLEAR-SWI eliminates the artefacts common in standard, single-echo SWI, reduces signal dropouts and improves image homogeneity and contrast-to-noise. Applied clinically, in a study of brain tumor patients, CLEAR-SWI was free of the artefacts which affected standard, single-echo SWI.

A novel phantom with dia- and paramagnetic substructure for quantitative susceptibility mapping and relaxometry

PURPOSE

A phantom is presented in this study that allows for an experimental evaluation of QSM reconstruction algorithms. The phantom contains susceptibility producing particles with dia- and paramagnetic properties embedded in an MRI visible medium and is suitable to assess the performance of algorithms that attempt to separate isotropic dia- and paramagnetic susceptibility at the sub-voxel level.

METHODS

The phantom was built from calcium carbonate (diamagnetic) and tungsten carbide particles (paramagnetic) embedded in gelatin and surrounded by agarose gel. Different mass fractions and mixing ratios of both susceptibility sources were used. Gradient echo data were acquired at 1.5 T, 3 T and 7 T. Susceptibility maps were calculated using the MEDI toolbox and relaxation rates ΔR₂^∗ were determined using exponential fitting.

RESULTS

Relaxation rates as well as susceptibility values generally coincide with the theoretical values for particles fulfilling the assumptions of the the static dephasing regime with stronger deviations for relaxation rates at higher field strength and for high susceptibility values. MRI raw data are available for free academic use as supplementary material.

CONCLUSIONS

In this study, a susceptibility phantom is presented that might find its application in the development and quantitative validation of current and future QSM reconstruction algorithms which aim to separate the influence of isotropic dia- and paramagnetic substructure in quantitative susceptibility mapping.

Deep grey matter quantitative susceptibility mapping from small spatial coverages using deep learning

INTRODUCTION

Quantitative Susceptibility Mapping (QSM) is generally acquired with full brain coverage, even though many QSM brain-iron studies focus on the deep grey matter (DGM) region only. Reducing the spatial coverage to the DGM vicinity can substantially shorten the scan time or enhance the spatial resolution without increasing scan time; however, this may lead to significant DGM susceptibility underestimation.

METHOD

A recently proposed deep learning-based QSM method, namely xQSM, is investigated to assess the accuracy of dipole inversion on reduced brain coverages. The xQSM method is compared with two conventional dipole inversion methods using simulated and in vivo experiments from 4 healthy subjects at 3T. Pre-processed magnetic field maps are extended symmetrically from the centre of globus pallidus in the coronal plane to simulate QSM acquisitions of difference spatial coverages, ranging from 100% (∼32mm) to 400% (∼128mm) of the actual DGM physical size.

RESULTS

The proposed xQSM network led to the lowest DGM contrast loss in both simulated and in vivo subjects, with the smallest susceptibility variation range across all spatial coverages. For the digital brain phantom simulation, xQSM improved the DGM susceptibility underestimation more than 20% in small spatial coverages, as compared to conventional methods. For the in vivo acquisition, less than 5% DGM susceptibility error was achieved in 48mm axial slabs using the xQSM network, while a minimum of 112mm coverage was required for conventional methods. It is also shown that the background field removal process performed worse in reduced brain coverages, which further deteriorated the subsequent dipole inversion.

CONCLUSION

The recently proposed deep learning-based xQSM method significantly improves the accuracy of DGM QSM from small spatial coverages as compared with conventional QSM algorithms, which can shorten DGM QSM acquisition time substantially.

The PHU-NET: A robust phase unwrapping method for MRI based on deep learning

PURPOSE

This work was aimed at designing a deep-learning-based approach for MR image phase unwrapping to improve the robustness and efficiency of traditional methods.

METHODS

A deep learning network called PHU-NET was designed for MR phase unwrapping. In this network, a novel training data generation method was proposed to simulate the wrapping patterns in MR phase images. The wrapping boundary and wrapping counts were explicitly estimated and used for network training. The proposed method was quantitatively evaluated and compared to other methods using a number of simulated datasets with varying signal-to-noise ratio (SNR) and MR phase images from various parts of the human body.

RESULTS

The results showed that our method performed better in the simulated data even under an extremely low SNR. The proposed method had less residual wrapping in the images from various parts of human body and worked well in the presence of severe anatomical discontinuity. Our method was also advantageous in terms of computational efficiency compared to the traditional methods.

CONCLUSION

This work proposed a robust and computationally efficient MR phase unwrapping method based on a deep learning network, which has promising performance in applications using MR phase information.

Investigating the effect of flow compensation and quantitative susceptibility mapping method on the accuracy of venous susceptibility measurement

Quantitative susceptibility mapping (QSM) is a promising non-invasive method for obtaining information relating to oxygen metabolism. However, the optimal acquisition sequence and QSM reconstruction method for reliable venous susceptibility measurements are unknown. Full flow compensation is generally recommended to correct for the influence of venous blood flow, although the effect of flow compensation on the accuracy of venous susceptibility values has not been systematically evaluated. In this study, we investigated the effect of different acquisition sequences, including different flow compensation schemes, and different QSM reconstruction methods on venous susceptibilities. Ten healthy subjects were scanned with five or six distinct QSM sequence designs using monopolar readout gradients and different flow compensation schemes. All data sets were processed using six different QSM pipelines and venous blood susceptibility was evaluated in whole-brain segmentations of the venous vasculature and single veins. The quality of vein segmentations and the accuracy of venous susceptibility values were analyzed and compared between all combinations of sequences and reconstruction methods. The influence of the QSM reconstruction method on average venous susceptibility values was found to be 2.7-11.6 times greater than the influence of the acquisition sequence, including flow compensation. The majority of the investigated QSM reconstruction methods tended to underestimate venous susceptibility values in the vein segmentations that were obtained. In summary, we found that multi-echo gradient-echo acquisition sequences without full flow compensation yielded venous susceptibility values comparable to sequences with full flow compensation. However, the QSM reconstruction method had a great influence on susceptibility values and thus needs to be selected carefully for accurate venous QSM.

QSM reconstruction challenge 2.0: A realistic in silico head phantom for MRI data simulation and evaluation of susceptibility mapping procedures

PURPOSE

To create a realistic in silico head phantom for the second QSM reconstruction challenge and for future evaluations of processing algorithms for QSM.

METHODS

We created a digital whole-head tissue property phantom by segmenting and postprocessing high-resolution (0.64 mm isotropic), multiparametric MRI data acquired at 7 T from a healthy volunteer. We simulated the steady-state magnetization at 7 T using a Bloch simulator and mimicked a Cartesian sampling scheme through Fourier-based processing. Computer code for generating the phantom and performing the MR simulation was designed to facilitate flexible modifications of the phantom in the future, such as the inclusion of pathologies as well as the simulation of a wide range of acquisition protocols. Specifically, the following parameters and effects were implemented: TR and TE, voxel size, background fields, and RF phase biases. Diffusion-weighted imaging phantom data are provided, allowing future investigations of tissue-microstructure effects in phase and QSM algorithms.

RESULTS

The brain part of the phantom featured realistic morphology with spatial variations in relaxation and susceptibility values similar to the in vivo setting. We demonstrated some of the phantom's properties, including the possibility of generating phase data with nonlinear evolution over TE due to partial-volume effects or complex distributions of frequency shifts within the voxel.

CONCLUSION

The presented phantom and computer programs are publicly available and may serve as a ground truth in future assessments of the faithfulness of quantitative susceptibility reconstruction algorithms.

Emerging Applications for Quantitative Susceptibility Mapping in the Detection of Traumatic Brain Injury Pathology

Traumatic brain injury (TBI) is a common but heterogeneous injury underpinned by numerous complex and interrelated pathophysiological mechanisms. An essential trace element, iron is abundant within the brain and involved in many fundamental neurobiological processes, including oxygen transportation, oxidative phosphorylation, myelin production and maintenance, as well as neurotransmitter synthesis and metabolism. Excessive levels of iron are neurotoxic and thus iron homeostasis is tightly regulated in the brain, however, many details about the mechanisms by which this is achieved are yet to be elucidated. A key mediator of oxidative stress, mitochondrial dysfunction and neuroinflammatory response, iron dysregulation is an important contributor to secondary injury in TBI. Advances in neuroimaging that leverage magnetic susceptibility properties have enabled increasingly comprehensive investigations into the distribution and behaviour of iron in the brain amongst healthy individuals as well as disease states such as TBI. Quantitative Susceptibility Mapping (QSM) is an advanced neuroimaging technique that promises quantitative estimation of local magnetic susceptibility at the voxel level. In this review, we provide an overview of brain iron and its homeostasis, describe recent advances enabling applications of QSM within the context of TBI and summarise the current state of the literature. Although limited, the emergent research suggests that QSM is a promising neuroimaging technique that can be used to investigate a host of pathophysiological changes that are associated with TBI.

Multiecho complex total field inversion method (mcTFI) for improved signal modeling in quantitative susceptibility mapping

PURPOSE

Typical quantitative susceptibility mapping (QSM) reconstruction steps consist of first estimating the magnetization field from the gradient-echo images, and then reconstructing the susceptibility map from the estimated field. The errors from the field-estimation steps may propagate into the final QSM map, and the noise in the estimated field map may no longer be zero-mean Gaussian noise, thus, causing streaking artifacts in the resulting QSM. A multiecho complex total field inversion (mcTFI) method was developed to compute the susceptibility map directly from the multiecho gradient echo images using an improved signal model that retains the Gaussian noise property in the complex domain. It showed improvements in QSM reconstruction over the conventional field-to-source inversion.

METHODS

The proposed mcTFI method was compared with the nonlinear total field inversion (nTFI) method in a numerical brain with hemorrhage and calcification, the numerical brains provided by the QSM Challenge 2.0, 18 brains with intracerebral hemorrhage scanned at 3T, and 6 healthy brains scanned at 7T.

RESULTS

Compared with nTFI, the proposed mcTFI showed more accurate QSM reconstruction around the lesions in the numerical simulations. The mcTFI reconstructed QSM also showed the best image quality with the least artifacts in the brains with intracerebral hemorrhage scanned at 3T and healthy brains scanned at 7T.

CONCLUSION

The proposed multiecho complex total field inversion improved QSM reconstruction over traditional field-to-source inversion through better signal modeling.

Multi-echo gradient-recalled-echo phase unwrapping using a Nyquist sampled virtual echo train in the presence of high-field gradients

PURPOSE

To develop a spatio-temporal approach to accurately unwrap multi-echo gradient-recalled echo phase in the presence of high-field gradients.

THEORY AND METHODS

Using the virtual echo-based Nyquist sampled (VENyS) algorithm, the temporal unwrapping procedure is modified by introduction of one or more virtual echoes between the first lower and the immediate higher echo, so as to reinstate the Nyquist condition at locations with high-field gradients. An iterative extension of the VENyS algorithm maintains spatial continuity by adjusting the phase rotations to make the neighborhood phase differences less than π. The algorithm is evaluated using simulated data, Gadolinium contrast-doped phantom, and in vivo brain, abdomen, and chest data sets acquired at 3 T and 9.4 T. The unwrapping performance is compared with the standard temporal unwrapping algorithm used in the morphology-enabled dipole inversion-QSM pipeline as a benchmark for validation.

RESULTS

Quantitative evaluation using numerical phantom showed significant reduction in unwrapping errors in regions of large field gradients, and the unwrapped phase revealed an exact match with the linear concentration profile of vials in a gadolinium contrast-doped phantom data acquired at 9.4 T. Without the need for additional spatial unwrapping, the iterative VENyS algorithm was able to generate spatially continuous phase images. Application to in vivo data resulted in better unwrapping performance, especially in regions with large susceptibility changes such as the air/tissue interface.

CONCLUSION

The iterative VENyS algorithm serves as a robust unwrapping method for multi-echo gradient-recalled echo phase in the presence of high-field gradients.

Automatic determination of the regularization weighting for wavelet-based compressed sensing MRI reconstructions

PURPOSE

To present a method that automatically, rapidly, and in a noniterative manner determines the regularization weighting for wavelet-based compressed sensing reconstructions. This method determines level-specific regularization weighting factors from the wavelet transform of the image obtained from zero-filling in k-space.

METHODS

We compare reconstruction results obtained by our method, λ auto , to the ones obtained by the L-curve, λ Lcurve , and the minimum NMSE, λ NMSE . The comparisons are done using in vivo data; then, simulations are used to analyze the impact of undersampling and noise. We use NMSE, Pearson's correlation coefficient, high-frequency error norm, and structural similarity as reconstruction quality indices.

RESULTS

Our method, λ auto , provides improved reconstructed image quality to that obtained by λ Lcurve regardless of undersampling or SNR and comparable quality to λ NMSE at high SNR. The method determines the regularization weighting prospectively with negligible computational time.

CONCLUSION

Our main finding is an automatic, fast, noniterative, and robust procedure to determine the regularization weighting. The impact of this method is to enable prospective and tuning-free wavelet-based compressed sensing reconstructions.

Phase unwrapping with a rapid opensource minimum spanning tree algorithm (ROMEO)

PURPOSE

To develop a rapid and accurate MRI phase-unwrapping technique for challenging phase topographies encountered at high magnetic fields, around metal implants, or postoperative cavities, which is sufficiently fast to be applied to large-group studies including Quantitative Susceptibility Mapping and functional MRI (with phase-based distortion correction).

METHODS

The proposed path-following phase-unwrapping algorithm, ROMEO, estimates the coherence of the signal both in space-using MRI magnitude and phase information-and over time, assuming approximately linear temporal phase evolution. This information is combined to form a quality map that guides the unwrapping along a 3D path through the object using a computationally efficient minimum spanning tree algorithm. ROMEO was tested against the two most commonly used exact phase-unwrapping methods, PRELUDE and BEST PATH, in simulated topographies and at several field strengths: in 3T and 7T in vivo human head images and 9.4T ex vivo rat head images.

RESULTS

ROMEO was more reliable than PRELUDE and BEST PATH, yielding unwrapping results with excellent temporal stability for multi-echo or multi-time-point data. It does not require image masking and delivers results within seconds, even in large, highly wrapped multi-echo data sets (eg, 9 seconds for a 7T head data set with 31 echoes and a 208 × 208 × 96 matrix size).

CONCLUSION

Overall, ROMEO was both faster and more accurate than PRELUDE and BEST PATH, delivering exact results within seconds, which is well below typical image acquisition times, enabling potential on-console application.

QSMART: Quantitative susceptibility mapping artifact reduction technique

PURPOSE

Quantitative susceptibility mapping (QSM) is a novel MR technique that allows mapping of tissue susceptibility values from MR phase images. QSM is an ill-conditioned inverse problem, and although several methods have been proposed in the field, in the presence of a wide range of susceptibility sources, streaking artifacts appear around high susceptibility regions and contaminate the whole QSM map. QSMART is a post-processing pipeline that uses two-stage parallel inversion to reduce the streaking artifacts and remove banding artifact at the cortical surface and around the vasculature.

METHOD

Tissue and vein susceptibility values were separately estimated by generating a mask of vasculature driven from the magnitude data using a Frangi filter. Spatially dependent filtering was used for the background field removal step and the two susceptibility estimates were combined in the final QSM map. QSMART was compared to RESHARP/iLSQR and V-SHARP/iLSQR inversion in a numerical phantom, 7T in vivo single and multiple-orientation scans, 9.4T ex vivo mouse data, and 4.7T in vivo rat brain with induced focal ischemia.

RESULTS

Spatially dependent filtering showed better suppression of phase artifacts near cortex compared to RESHARP and V-SHARP, while preserving voxels located within regions of interest without brain edge erosion. QSMART showed successful reduction of streaking artifacts as well as improved contrast between different brain tissues compared to the QSM maps obtained by RESHARP/iLSQR and V-SHARP/iLSQR.

CONCLUSION

QSMART can reduce QSM artifacts to enable more robust estimation of susceptibility values in vivo and ex vivo.