Quantitative Susceptibility Mapping Reconstruction Research Articles

Quantitative susceptibility mapping in amyotrophic lateral sclerosis: automatic quantification of the magnetic susceptibility in the subcortical nuclei

Objective: Previous studies have suggested a link between dysregulation of cortical iron levels and neuronal loss in amyotrophic lateral sclerosis (ALS) patients. However, few studies have reported differences in quantitative susceptibility mapping (QSM) values in subcortical nuclei between patients with ALS and healthy controls (HCs). Methods: MRI was performed using a 3 Tesla Prisma scanner (64-channel head coil), including 3D T1-MPRAGE and multi-echo 3D GRE for QSM reconstruction. Automated QSM segmentation was used to measure susceptibility values in the subcortical nuclei, which were compared between the groups. Correlations with clinical scales were analyzed. Group comparisons were performed using independent t-tests, with p < 0.05 considered significant. Correlations were assessed using Pearson’s correlation, with p < 0.05 considered significant. Cohen’s d was reported to compare the standardized mean difference (SMD) of QSM. Results: Twelve patients with limb-onset ALS (mean age 48.7 years, 75% male) and 13 age-, sex-, and handedness-matched HCs (mean age 44.6 years, 69% male) were included. Compared to HCs, ALS patients demonstrated significantly lower susceptibility in the left caudate nucleus (CN) (SMD = −0.845), right CN (SMD = −0.851), whole CN (SMD = −1.016), and left subthalamic nucleus (STN) (SMD = −1.000). Susceptibility in the left putamen (SMD = −0.857), left thalamus (SMD = −1.081), and whole thalamus (SMD = −0.968) was significantly higher in the patients. The susceptibility of the substantia nigra (SN), CN, and pulvinar was positively correlated with disease duration. Conclusions: QSM detects abnormal iron accumulation patterns in the subcortical gray matter of ALS patients, which correlates with disease characteristics, supporting its potential as a neuroimaging biomarker.

DF-QSM: Data Fidelity based Hybrid Approach for Improved Quantitative Susceptibility Mapping of the Brain.

Quantitative Susceptibility Mapping (QSM) is an advanced magnetic resonance imaging (MRI) technique to quantify the magnetic susceptibility of the tissue under investigation. Deep learning methods have shown promising results in deconvolving the susceptibility distribution from the measured local field obtained from the MR phase. Although existing deep learning based QSM methods can produce high-quality reconstruction, they are highly biased toward training data distribution with less scope for generalizability. This work proposes a hybrid two-step reconstruction approach to improve deep learning based QSM reconstruction. The susceptibility map prediction obtained from the deep learning methods has been refined in the framework developed in this work to ensure consistency with the measured local field. The developed method was validated on existing deep learning and model-based deep learning methods for susceptibility mapping of the brain. The developed method resulted in improved reconstruction for MRI volumes obtained with different acquisition settings, including deep learning models trained on constrained (limited) data settings.

Spinet-QSM: model-based deep learning with schatten p-norm regularization for improved quantitative susceptibility mapping.

Quantitative susceptibility mapping (QSM) provides an estimate of the magnetic susceptibility of tissue using magnetic resonance (MR) phase measurements. The tissue magnetic susceptibility (source) from the measured magnetic field distribution/local tissue field (effect) inherent in the MR phase images is estimated by numerically solving the inverse source-effect problem. This study aims to develop an effective model-based deep-learning framework to solve the inverse problem of QSM. This work proposes a Schatten -norm-driven model-based deep learning framework for QSM with a learnable norm parameter to adapt to the data. In contrast to other model-based architectures that enforce the l -norm or l -norm for the denoiser, the proposed approach can enforce any -norm ( ) on a trainable regulariser. The proposed method was compared with deep learning-based approaches, such as QSMnet, and model-based deep learning approaches, such as learned proximal convolutional neural network (LPCNN). Reconstructions performed using 77 imaging volumes with different acquisition protocols and clinical conditions, such as hemorrhage and multiple sclerosis, showed that the proposed approach outperformed existing state-of-the-art methods by a significant margin in terms of quantitative merits. The proposed SpiNet-QSM showed a consistent improvement of at least 5% in terms of the high-frequency error norm (HFEN) and normalized root mean squared error (NRMSE) over other QSM reconstruction methods with limited training data.

A subject-specific unsupervised deep learning method for quantitative susceptibility mapping using implicit neural representation

Quantitative susceptibility mapping (QSM) is an MRI-based technique that estimates the underlying tissue magnetic susceptibility based on phase signal. Deep learning (DL)-based methods have shown promise in handling the challenging ill-posed inverse problem for QSM reconstruction. However, they require extensive paired training data that are typically unavailable and suffer from generalization problems. Recent model-incorporated DL approaches also overlook the non-local effect of the tissue phase in applying the source-to-field forward model due to patch-based training constraint, resulting in a discrepancy between the prediction and measurement and subsequently suboptimal QSM reconstruction. This study proposes an unsupervised and subject-specific DL method for QSM reconstruction based on implicit neural representation (INR), referred to as INR-QSM. INR has emerged as a powerful framework for learning a high-quality continuous representation of the signal (image) by exploiting its internal information without training labels. In INR-QSM, the desired susceptibility map is represented as a continuous function of the spatial coordinates, parameterized by a fully-connected neural network. The weights are learned by minimizing a loss function that includes a data fidelity term incorporated by the physical model and regularization terms. Additionally, a novel phase compensation strategy is proposed for the first time to account for the non-local effect of tissue phase in data consistency calculation to make the physical model more accurate. Our experiments show that INR-QSM outperforms traditional established QSM reconstruction methods and the compared unsupervised DL method both qualitatively and quantitatively, and is competitive against supervised DL methods under data perturbations.

Plug-and-Play latent feature editing for orientation-adaptive quantitative susceptibility mapping neural networks

Quantitative susceptibility mapping (QSM) is a post-processing technique for deriving tissue magnetic susceptibility distribution from MRI phase measurements. Deep learning (DL) algorithms hold great potential for solving the ill-posed QSM reconstruction problem. However, a significant challenge facing current DL-QSM approaches is their limited adaptability to magnetic dipole field orientation variations during training and testing. In this work, we propose a novel Orientation-Adaptive Latent Feature Editing (OA-LFE) module to learn the encoding of acquisition orientation vectors and seamlessly integrate them into the latent features of deep networks. Importantly, it can be directly Plug-and-Play (PnP) into various existing DL-QSM architectures, enabling reconstructions of QSM from arbitrary magnetic dipole orientations. Its effectiveness is demonstrated by combining the OA-LFE module into our previously proposed phase-to-susceptibility single-step instant QSM (iQSM) network, which was initially tailored for pure-axial acquisitions. The proposed OA-LFE-empowered iQSM, which we refer to as iQSM+, is trained in a simulated-supervised manner on a specially-designed simulation brain dataset. Comprehensive experiments are conducted on simulated and in vivo human brain datasets, encompassing subjects ranging from healthy individuals to those with pathological conditions. These experiments involve various MRI platforms (3T and 7T) and aim to compare our proposed iQSM+ against several established QSM reconstruction frameworks, including the original iQSM. The iQSM+ yields QSM images with significantly improved accuracies and mitigates artifacts, surpassing other state-of-the-art DL-QSM algorithms. The PnP OA-LFE module’s versatility was further demonstrated by its successful application to xQSM, a distinct DL-QSM network for dipole inversion. In conclusion, this work introduces a new DL paradigm, allowing researchers to develop innovative QSM methods without requiring a complete overhaul of their existing architectures.

Harmonic field extension for QSM with reduced spatial coverage using physics-informed generative adversarial network

Quantitative susceptibility mapping (QSM) is frequently employed in investigating brain iron related to brain development and diseases within deep gray matter (DGM). Nonetheless, the acquisition of whole-brain QSM data is time-intensive. An alternative approach, focusing the QSM specifically on areas of interest such as the DGM by reducing the field-of-view (FOV), can significantly decrease scan times. However, severe susceptibility value underestimations have been reported during QSM reconstruction with a limited FOV, largely attributable to artifacts from incorrect background field removal in the boundary region. This presents a considerable barrier to the clinical use of QSM with small spatial coverages using conventional methods alone.To mitigate the propagation of these errors, we proposed a harmonic field extension method based on a physics-informed generative adversarial network. Both quantitative and qualitative results demonstrate that our method outperforms conventional methods and delivers results comparable to those obtained with full FOV. Furthermore, we demonstrate the versatility of our method by applying it to data acquired prospectively with limited FOV and to data from patients with Parkinson's disease. The method has shown significant improvements in local field results, with QSM outcomes. In a clear illustration of its feasibility and effectiveness in real clinical environments, our proposed method addresses the prevalent issue of susceptibility underestimation in QSM with small spatial coverage.

NoQSM-net: Combining Convolutional Neural Network With Numerical Optimization Algorithm for Quantitative Susceptibility Mapping Reconstruction

NoQSM-net: Combining Convolutional Neural Network With Numerical Optimization Algorithm for Quantitative Susceptibility Mapping Reconstruction

A latent code based multi-variable modulation network for susceptibility mapping.

Quantitative susceptibility mapping (QSM) is a technique for obtaining quantitative information on tissue susceptibility and has shown promising potential for clinical applications, in which the magnetic susceptibility is calculated by solving an ill-posed inverse problem. Recently, deep learning-based methods are proposed to address this issue, but the diversity of data distribution was not well considered, and thus the model generalization is limited in clinical applications. In this paper, we propose a Latent Code based Multi-Variable modulation network for QSM reconstruction (LCMnet). Particularly, a specific modulation module is exploited to incorporate three variables, i.e., field map, magnitude image, and initial susceptibility. The latent code in the modulated convolution is learned from feature maps of the field data using the encoder-decoder framework. The susceptibility map pre-estimated from simple thresholding is the constant input of the module, thereby enhancing the network stability and accelerating training convergence. As another input, multi-level features generated by a cross-fusion block integrate the information of field and magnitude data effectively. Experimental results on in vivo human brain data, challenge data, clinical data and synthetic data demonstrate that the proposed method LCMnet can achieve outstanding performance on accurate susceptibility measurement and also excellent generalization.

Exploiting gradient-echo frequency evolution: Probing white matter microstructure and extracting bulk susceptibility-induced frequency for quantitative susceptibility mapping.

This work is to investigate the microstructure-induced frequency shift in white matter (WM) with crossing fibers and to separate the microstructure-related frequency shift from the bulk susceptibility-induced frequency shift by model fitting the gradient-echo (GRE) frequency evolution for potentially more accurate quantitative susceptibility mapping (QSM). A hollow-cylinder fiber model (HCFM) with two fiber populations was developed to investigate GRE frequency evolutions in WM voxels with microstructural orientation dispersion. The simulated and experimentally measured TE-dependent local frequency shift was then fitted to a simplified frequency evolution model to obtain a microstructure-related frequency difference parameter ( ) and a TE-independent bulk susceptibility-induced frequency shift ( ). The obtained was then used for QSM reconstruction. Reconstruction performances were evaluated using a numerical head phantom and in vivo data and then compared to other multi-echo combination methods. GRE frequency evolutions and -based tissue parameters in both parallel and crossing fibers determined from our simulations were comparable to those observed in vivo. The TE-dependent frequency fitting method outperformed other multi-echo combination methods in estimating in simulations. The fitted , , and QSM could be improved further by navigator-based B0 fluctuation correction. A HCFM with two fiber populations can be used to characterize microstructure-induced frequency shifts in WM regions with crossing fibers. HCFM-based TE-dependent frequency fitting provides tissue contrast related to microstructure ( ) and in addition may help improve the quantification accuracy of and the corresponding QSM.

Free-breathing high isotropic resolution quantitative susceptibility mapping (QSM) of liver using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction.

To enable free-breathing and high isotropic resolution liver quantitative susceptibility mapping (QSM) using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction. Using 3D multi-echo UTE cones MRI, a respiratory motion was estimated from the k-space center of the imaging data. After sorting the k-space data with estimated motion, respiratory motion state-resolved reconstruction was performed for multi-echo data followed by nonlinear least-squares fitting for proton density fat fraction (PDFF), , and fat-corrected B0 field maps. PDFF and B0 field maps were subsequently used for QSM reconstruction. The proposed method was compared with motion-averaged (gridding) reconstruction and conventional 3D multi-echo Cartesian MRI in moving gadolinium phantom and in vivo studies. Region of interest (ROI)-based linear regression analysis was performed on these methods to investigate correlations between gadolinium concentration and QSM in the phantom study and between and QSM in in vivo study. Cones with motion-resolved reconstruction showed sharper image quality compared to motion-averaged reconstruction with a substantial reduction of motion artifacts in both moving phantom and in vivo studies. For ROI-based linear regression analysis of the phantom study, susceptibility values from cones with motion-resolved reconstruction ( = 0.31 0.05, = 0.999) and Cartesian without motion ( = 0.32 0.04, = 1.000) showed linear relationships with gadolinium concentrations and showed good agreement with each other. For in vivo, motion-resolved reconstruction showed higher goodness of fit ( = 0.00261 0.524, = 0.977) compared to motion-averaged reconstruction ( = 0.0021 0.572, = 0.723) in ROI-based linear regression analysis between and QSM. Feasibility of free-breathing liver QSM was demonstrated with motion-resolved 3D multi-echo UTE cones MRI, achieving high isotropic resolution currently unachievable in conventional Cartesian MRI.

Quantitative susceptibility mapping using multi-channel convolutional neural networks with dipole-adaptive multi-frequency inputs.

Quantitative susceptibility mapping (QSM) quantifies the distribution of magnetic susceptibility and shows great potential in assessing tissue contents such as iron, myelin, and calcium in numerous brain diseases. The accuracy of QSM reconstruction was challenged by an ill-posed field-to-susceptibility inversion problem, which is related to the impaired information near the zero-frequency response of the dipole kernel. Recently, deep learning methods demonstrated great capability in improving the accuracy and efficiency of QSM reconstruction. However, the construction of neural networks in most deep learning-based QSM methods did not take the intrinsic nature of the dipole kernel into account. In this study, we propose a dipole kernel-adaptive multi-channel convolutional neural network (DIAM-CNN) method for the dipole inversion problem in QSM. DIAM-CNN first divided the original tissue field into high-fidelity and low-fidelity components by thresholding the dipole kernel in the frequency domain, and it then inputs the two components as additional channels into a multichannel 3D Unet. QSM maps from the calculation of susceptibility through multiple orientation sampling (COSMOS) were used as training labels and evaluation reference. DIAM-CNN was compared with two conventional model-based methods [morphology enabled dipole inversion (MEDI) and improved sparse linear equation and least squares (iLSQR) and one deep learning method (QSMnet)]. High-frequency error norm (HFEN), peak signal-to-noise-ratio (PSNR), normalized root mean squared error (NRMSE), and the structural similarity index (SSIM) were reported for quantitative comparisons. Experiments on healthy volunteers demonstrated that the DIAM-CNN results had superior image quality to those of the MEDI, iLSQR, or QSMnet results. Experiments on data with simulated hemorrhagic lesions demonstrated that DIAM-CNN produced fewer shadow artifacts around the bleeding lesion than the compared methods. This study demonstrates that the incorporation of dipole-related knowledge into the network construction has a potential to improve deep learning-based QSM reconstruction.

MsQSM: Morphology-based self-supervised deep learning for quantitative susceptibility mapping

Quantitative susceptibility mapping (QSM) has been applied to the measurement of iron deposition and the auxiliary diagnosis of neurodegenerative disease. There still exists a dipole inversion problem in QSM reconstruction. Recently, deep learning approaches have been proposed to resolve this problem. However, most of these approaches are supervised methods that need pairs of the input phase and ground-truth. It remains a challenge to train a model for all resolutions without using the ground-truth and only using one resolution data. To address this, we proposed a self-supervised QSM deep learning method based on morphology. It consists of a morphological QSM builder to decouple the dependency of the QSM on acquisition resolution, and a morphological loss to reduce artifacts effectively and save training time efficiently. The proposed method can reconstruct arbitrary resolution QSM on both human data and animal data, regardless of whether the resolution is higher or lower than that of the training set. Our method outperforms the previous best unsupervised method with a 3.6% higher peak signal-to-noise ratio, 16.2% lower normalized root mean square error, and 22.1% lower high-frequency error norm. The morphological loss reduces training time by 22.1% with respect to the cycle gradient loss used in the previous unsupervised methods. Experimental results show that the proposed method accurately measures QSM with arbitrary resolutions, and achieves state-of-the-art results among unsupervised deep learning methods. Research on applications in neurodegenerative diseases found that our method is robust enough to measure significant increase in striatal magnetic susceptibility in patients during Alzheimer’s disease progression, as well as significant increase in substantia nigra susceptibility in Parkinson’s disease patients, and can be used as an auxiliary differential diagnosis tool for Alzheimer’s disease and Parkinson’s disease.

Tilted quantitative susceptibility mapping at oblique MRI (tiltQSM)

ObjectiveIf 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. MethodsWe proposed two methods for QSM reconstruction from slice-tilted MRI phase image (tiltQSM): 1) rotData per anti-tilting phase image rotation back into the B0-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 B0-upright system or in B0-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). ResultsThe 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). ConclusionsTo 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.

Test-retest analysis of cerebral oxygen extraction estimates in healthy volunteers: comparison of methods based on quantitative susceptibility mapping and dynamic susceptibility contrast magnetic resonance imaging

BackgroundEstimation of the oxygen extraction fraction (OEF) by quantitative susceptibility mapping (QSM) magnetic resonance imaging (MRI) is promising but requires systematic evaluation. Extraction of OEF-related information from the tissue residue function in dynamic susceptibility contrast MRI (DSC-MRI) has also been proposed. In this study, whole-brain OEF repeatability was investigated, as well as the relationships between QSM-based OEF and DSC-MRI-based parameters, i.e., mean transit time (MTT) and an oxygen extraction index, referred to as apparent OEF (AOEF). MethodTest-retest data were obtained from 20 healthy volunteers at 3 T. QSM maps were reconstructed from 3D gradient-echo MRI phase data, using morphology-enabled dipole inversion. DSC-MRI was accomplished using gradient-echo MRI at a temporal resolution of 1.24 s. ResultsThe whole-brain QSM-based OEF was (40.4±4.8) % and, in combination with a previously published cerebral blood flow (CBF) estimate, this corresponds to a cerebral metabolic rate of oxygen level of CMRO2 = 3.36 ml O2/min/100 g. The intra-class correlation coefficient [ICC(2,1)] for OEF test-retest data was 0.73. The MTT-versus-OEF and AOEF-versus-OEF relationships showed correlation coefficients of 0.61 (p = 0.004) and 0.52 (p = 0.019), respectively. DiscussionQSM-based OEF showed a convincing absolute level and good test-retest results in terms of the ICC. Moderate to good correlations between QSM-based OEF and DSC-MRI-based parameters were observed. The present results constitute an indicator of the level of robustness that can be achieved without applying extraordinary resources in terms of MRI equipment, imaging protocol, QSM reconstruction, and OEF analysis.

Deep learning-regularized, single-step quantitative susceptibility mapping quantification.

The purpose of the current study was to develop deep learning-regularized, single-step quantitative susceptibility mapping (QSM) quantification, directly generating QSM from the total phase map. A deep learning-regularized, single-step QSM quantification model, named SS-POCSnet, was trained with datasets created using the QSM synthesis approach in QSM reconstruction challenge 2.0. In SS-POCSnet, a data fidelity term based on a single-step model was iteratively applied that combined the spherical mean value kernel and dipole model. Meanwhile, SS-POCSnet regularized susceptibility maps, avoiding underestimating susceptibility values. We evaluated the SS-POCSnet on 10 synthetic datasets, 24 clinical datasets with lesions of cerebral microbleed (CMB) and calcification, and 10 datasets with multiple sclerosis (MS).On synthetic datasets, SS-POCSnet showed the best performance among the methods evaluated, with a normalized root mean squared error of 37.3% ± 4.2%, susceptibility-tuned structured similarity index measure of 0.823 ± 0.02, high-frequency error norm of 37.0 ± 5.7, and peak signal-to-noise ratio of 42.8 ± 1.1. SS-POCSnet also reduced the underestimations of susceptibility values in deep brain nuclei compared with those from the other models evaluated. Furthermore, SS-POCSnet was sensitive to CMB/calcification and MS lesions, demonstrating its clinical applicability. Our method also supported variable imaging parameters, including matrix size and resolution. It was concluded that deep learning-regularized, single-step QSM quantification can mitigate underestimating susceptibility values in deep brain nuclei.

S2Q-Net: Mining the High-Pass Filtered Phase Data in Susceptibility Weighted Imaging for Quantitative Susceptibility Mapping.

Susceptibility weighted imaging (SWI) is a routine magnetic resonance imaging (MRI) sequence that combines the magnitude and high-pass filtered phase images to qualitatively enhance the image contrasts related to tissue susceptibility. Tremendous amounts of the high-pass filtered phase data with low signal to noise ratio and incomplete background field removal have thus been collected under default clinical settings. Since SWI cannot quantitatively estimate the susceptibility, it is thus non-trivial to derive quantitative susceptibility mapping (QSM) directly from these redundant phase data, which effectively promotes the mining of the SWI data collected previously. To this end, a novel deep learning based SWI-to-QSM-Net (S2Q-Net) is proposed for QSM reconstruction from SWI high-pass filtered phase data. S2Q-Net firstly estimates the edge maps of QSM to integrate edge prior into features, which benefits the network to reconstruct QSM with realistic and clear tissue boundaries. Furthermore, a novel Second-order Cross Dense Block is proposed in S2Q-Net, which can capture rich inter-region interactions to provide more non-local phase information related to local tissue susceptibility. Experimental results on both simulated and in-vivo data indicate its superiority over all the compared deep learning based QSM reconstruction methods.

Early differentiation of neurodegenerative diseases using the novel QSM technique: what is the biomarker of each disorder?

During neurodegenerative diseases, the brain undergoes morphological and pathological changes; Iron deposits are one of the causes of pathological changes in the brain. The Quantitative susceptibility mapping (QSM) technique, a type of magnetic resonance (MR) image reconstruction, is one of the newest diagnostic methods for iron deposits to detect changes in magnetic susceptibility. Numerous research projects have been conducted in this field. The purpose of writing this review article is to identify the first deep brain nuclei that undergo magnetic susceptibility changes during neurodegenerative diseases such as Alzheimer's or Parkinson's disease. The purpose of this article is to identify the brain nuclei that are prone to iron deposition in any specific disorder. In addition to the mentioned purpose, this paper proposes the optimal scan parameters and appropriate algorithms of each QSM reconstruction step by reviewing the results of different articles. As a result, The QSM technique can identify nuclei exposed to iron deposition in various neurodegenerative diseases. Also, the selection of scan parameters is different based on the sequence and purpose; an example of the parameters is placed in the tables. The BET toolbox in FSL, Laplacian-based phase-unwrapping process, the V_SHARP algorithm, and morphology-enabled dipole inversion (MEDI) method are the most widely used algorithms in various stages of QSM reconstruction.

Towards in vivo ground truth susceptibility for single-orientation deep learning QSM: A multi-orientation gradient-echo MRI dataset

Recently, deep neural networks have shown great potential for solving dipole inversion of quantitative susceptibility mapping (QSM) with improved results. However, these studies utilized their limited dataset for network training and inference, which may lead to untrustworthy conclusions. Thus, a common dataset is needed for a fair comparison between different QSM reconstruction networks. Additionally, finding an in vivo reference susceptibility map that matches acquired single-orientation phase data remains an open problem. Susceptibility tensor imaging (STI) χ33 and Calculation of Susceptibility through Multiple Orientation Sampling (COSMOS) are considered reference susceptibility candidates. However, a large number of multi-orientation GRE data for both STI and COSMOS reconstruction are now unavailable for training supervised neural networks for QSM. In this study, we reported the largest multi-orientation dataset, to the best of our knowledge in the QSM research field, with a total of 144 scans from 8 healthy subjects collected using a 3D GRE sequence from the same MR scanner. In addition, the parcellation of deep gray matter is also provided for automatically extracting susceptibility values. Five recently developed deep neural networks, i.e., xQSM, QSMnet, autoQSM, LPCNN, and MoDL-QSM were performed on this dataset. This potential data source could provide a common framework and labels to test the accuracy and robustness of deep neural networks for QSM reconstruction. This dataset has the potential to provide a benchmark of reference susceptibility for the deep learning-based QSM methods. Additionally, the trained COSMOS-labeled and χ33-labeled networks were tested on the pathological data to explore their potential applications. The data together with deep gray matter parcellation maps are now publicly available via an open repository at https://osf.io/yfms7/, and the raw multi-orientation GRE data are also available at https://osf.io/y6rc3/.

Unsupervised resolution-agnostic quantitative susceptibility mapping using adaptive instance normalization.

Quantitative susceptibility mapping (QSM) is a useful magnetic resonance imaging (MRI) technique that provides the spatial distribution of magnetic susceptibility values of tissues. QSMs can be obtained by deconvolving the dipole kernel from phase images, but the spectral nulls in the dipole kernel make the inversion ill-posed. In recent years, deep learning approaches have shown a comparable QSM reconstruction performance to the classic approaches, in addition to the fast reconstruction time. Most of the existing deep learning methods are, however, based on supervised learning, so matched pairs of input phase images and ground-truth maps are needed. Moreover, it was reported that the deep learning-based methods fail to reconstruct QSM when the resolution of test data is different from the trained resolution. To address this, here we propose an unsupervised resolution-agnostic QSM deep learning method. The proposed method does not require QSM labels for training and reconstructs QSM with various resolutions by using adaptive instance normalization. Experimental results and clinical validation confirm that the proposed method provides accurate QSM with various resolutions compared to other deep learning approaches, and shows competitive performance to the best classical approaches in addition to the ultra-fast reconstruction.

Notes on quantitative susceptibility mapping reconstruction accuracy under challenging conditions: Phantom measurements and simulations

Magnetic susceptibility can be assessed by quantitative susceptibility mapping (QSM), based on measured magnetic resonance imaging (MRI) phase data. The QSM reconstruction process is, however, mathematically challenging and still not fully robust. A signal-generating holmium [Ho(III)] aqueous solution with air-equivalent magnetic susceptibility was prepared, and used as a surrounding medium in a water phantom with tubes filled with a solution of gadolinium contrast agent at various concentrations. Extended analyses under controlled conditions were accomplished by simulations of the phantom construction. Without surrounding holmium solution, a gadolinium tube positioned centrally, parallel with B0, showed a susceptibility difference that agreed well with theoretical values, whereas a peripheral parallel tube position showed larger deviation. Orientation perpendicular to B0 resulted in less variation between the internal tube positions. Air-equivalent magnetic susceptibility corresponded to 16.5 mM Ho(III) solution. With surrounding holmium solution, several post-processing steps became challenging. Simulations indicated higher degree of underestimation when the theoretical susceptibility difference increased. Details in the mathematical implementation, for example, background field removal can strongly influence the result. Simulated results were, in part, unexpected, and provided awareness of limitations in the reconstruction technique, mainly related to conditions with large susceptibility differences between compartments.