Exploring linearity of deep neural network trained QSM: QSMnet+

Comparison between R2'-based and R2*-based χ-separation methods: A clinical evaluation in individuals with multiple sclerosis

Susceptibility source separation, or χ-separation, estimates diamagnetic (χdia) and paramagnetic susceptibility (χpara) signals in the brain using local field and R2' (= R2* - R2) maps. Recently proposed R2*-based χ-separation methods allow for χ-separation using only multi-echo gradient echo (ME-GRE) data, eliminating the need for additional data acquisition for R2 mapping. Although this approach reduces scan time and enhances clinical utility, the impact of missing R2 information remains a subject of exploration. In this study, we evaluate the viability of two previously proposed R2*-based χ-separation methods as alternatives to their R2'-based counterparts: model-based R2*-χ-separation versus χ-separation and deep learning-based χ-sepnet-R2* versus χ-sepnet-R2'. Their performances are assessed in individuals with multiple sclerosis (MS), comparing them with their corresponding R2'-based counterparts (i.e., R2*-χ-separation vs. χ-separation and χ-sepnet-R2* vs. χ-sepnet-R2'). The evaluations encompass qualitative visual assessments by experienced neuroradiologists and quantitative analyses, including region of interest analyses and linear regression analyses. Qualitatively, R2*-χ-separation tends to report higher χpara and χdia values compared with χ-separation, leading to less distinct lesion contrasts, while χ-sepnet-R2* closely aligns with χ-sepnet-R2'. Quantitative analysis reveals a robust correlation between both R2*-based methods and their R2'-based counterparts (r ≥ 0.88). Specifically, in the whole-brain voxels, χ-sepnet-R2* exhibits higher correlation and better linearity than R2*-χ-separation (χdia/χpara from R2*-χ-separation: r = 0.88/0.90, slope = 0.79/0.86; χdia/χpara from χ-sepnet-R2*: r = 0.90/0.92, slope = 0.99/0.97). In MS lesions, both R2*-based methods display comparable correlation and linearity (χdia/χpara from R2*-χ-separation: r = 0.90/0.91, slope = 0.98/0.91; χdia/χpara from χ-sepnet-R2*: r = 0.88/0.88, slope = 0.91/0.95). Notably, χ-sepnet-R2* demonstrates negligible offsets, whereas R2*-χ-separation exhibits relatively large offsets (0.02 ppm in the whole brain and 0.01 ppm in the MS lesions), potentially indicating the false presence of myelin or iron in MS lesions. Overall, both R2*-based χ-separation methods demonstrated their viability as alternatives to their R2'-based counterparts. χ-sepnet-R2* showed better alignment with its R2'-based counterpart with minimal susceptibility offsets, compared with R2*-χ-separation that reported higher χpara and χdia values compared with R2'-based χ-separation.

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.

Machine learning-based amide proton transfer imaging using partially synthetic training data

PURPOSE

Machine learning (ML) has been increasingly used to quantify CEST effect. ML models are typically trained using either measured data or fully simulated data. However, training with measured data often lacks sufficient training data, whereas training with fully simulated data may introduce bias because of limited simulations pools. This study introduces a new platform that combines simulated and measured components to generate partially synthetic CEST data, and to evaluate its feasibility for training ML models to predict amide proton transfer (APT) effect.

METHODS

Partially synthetic CEST signals were created using an inverse summation of APT effects from simulations and the other components from measurements. Training data were generated by varying APT simulation parameters and applying scaling factors to adjust the measured components, achieving a balance between simulation flexibility and fidelity. First, tissue-mimicking CEST signals along with ground truth information were created using multiple-pool model simulations to validate this method. Second, an ML model was trained individually on partially synthetic data, in vivo data, and fully simulated data, to predict APT effect in rat brains bearing 9 L tumors.

RESULTS

Experiments on tissue-mimicking data suggest that the ML method using the partially synthetic data is accurate in predicting APT. In vivo experiments suggest that our method provides more accurate and robust prediction than the training using in vivo data and fully synthetic data.

CONCLUSION

Partially synthetic CEST data can address the challenges in conventional ML methods.

Quantitative susceptibility mapping through model-based deep image prior (MoDIP)

The data-driven approach of supervised learning methods has limited applicability in solving dipole inversion in Quantitative Susceptibility Mapping (QSM) with varying scan parameters across different objects. To address this generalization issue in supervised QSM methods, we propose a novel training-free model-based unsupervised method called MoDIP (Model-based Deep Image Prior). MoDIP comprises a small, untrained network and a Data Fidelity Optimization (DFO) module. The network converges to an interim state, acting as an implicit prior for image regularization, while the optimization process enforces the physical model of QSM dipole inversion. Experimental results demonstrate MoDIP's excellent generalizability in solving QSM dipole inversion across different scan parameters. It exhibits robustness against pathological brain QSM, achieving over 32 % accuracy improvement than supervised deep learning methods. It is also 33 % more computationally efficient and runs 4 times faster than conventional DIP-based approaches, enabling 3D high-resolution image reconstruction in under 4.5 min.

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.

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.

A spatially adaptive regularization based three-dimensional reconstruction network for quantitative susceptibility mapping

Objective.Quantitative susceptibility mapping (QSM) is a new imaging technique for non-invasive characterization of the composition and microstructure ofin vivotissues, and it can be reconstructed from local field measurements by solving an ill-posed inverse problem. Even for deep learning networks, it is not an easy task to establish an accurate quantitative mapping between two physical quantities of different units, i.e. field shift in Hz and susceptibility value in ppm for QSM.Approach. In this paper, we propose a spatially adaptive regularization based three-dimensional reconstruction network SAQSM. A spatially adaptive module is specially designed and a set of them at different resolutions are inserted into the network decoder, playing a role of cross-modality based regularization constraint. Therefore, the exact information of both field and magnitude data is exploited to adjust the scale and shift of feature maps, and thus any information loss or deviation occurred in previous layers could be effectively corrected. The network encoding has a dynamic perceptual initialization, which enables the network to overcome receptive field intervals and also strengthens its ability to detect features of various sizes.Main results. Experimental results on the brain data of healthy volunteers, clinical hemorrhage and simulated phantom with calcification demonstrate that SAQSM can achieve more accurate reconstruction with less susceptibility artifacts, while perform well on the stability and generalization even for severe lesion areas.Significance. This proposed framework may provide a valuable paradigm to quantitative mapping or multimodal reconstruction.

Model resolution-based deconvolution for improved quantitative susceptibility mapping

Quantitative susceptibility mapping (QSM) utilizes the relationship between the measured local field and the unknown susceptibility map to perform dipole deconvolution. The aim of this work is to introduce and systematically evaluate the model resolution-based deconvolution for improved estimation of the susceptibility map obtained using the thresholded k-space division (TKD). A two-step approach has been proposed, wherein the first step involves the TKD susceptibility map computation and the second step involves the correction of this susceptibility map using the model-resolution matrix. The TKD-estimated susceptibility map can be expressed as the weighted average of the true susceptibility map, where the weights are determined by the rows of the model-resolution matrix, and hence a deconvolution of the TKD susceptibility map using the model-resolution matrix yields a better approximation to the true susceptibility map. The model resolution-based deconvolution is realized using closed-form, iterative, and sparsity-regularized implementations. The proposed approach was compared with L2 regularization, TKD, rescaled TKD in superfast dipole inversion, the modulated closed-form method, and iterative dipole inversion, as well as sparsity-regularized dipole inversion. It was observed that the proposed approach showed a substantial reduction in the streaking artifacts across 94 test volumes considered in this study. The proposed approach also showed better error reduction and edge preservation compared with other approaches. The proposed model resolution-based deconvolution compensates for the truncation of zero coefficients in the dipole kernel at the magic angle and hence provides a closer approximation to the true susceptibility map compared with other direct methods.

Cortical Iron Accumulation as an Imaging Marker for Neurodegeneration in Clinical Cognitive Impairment Spectrum: A Quantitative Susceptibility Mapping Study

OBJECTIVE

Cortical iron deposition has recently been shown to occur in Alzheimer's disease (AD). In this study, we aimed to evaluate how cortical gray matter iron, measured using quantitative susceptibility mapping (QSM), differs in the clinical cognitive impairment spectrum.

MATERIALS AND METHODS

This retrospective study evaluated 73 participants (mean age ± standard deviation, 66.7 ± 7.6 years; 52 females and 21 males) with normal cognition (NC), 158 patients with mild cognitive impairment (MCI), and 48 patients with AD dementia. The participants underwent brain magnetic resonance imaging using a three-dimensional multi-dynamic multi-echo sequence on a 3-T scanner. We employed a deep neural network (QSMnet+) and used automatic segmentation software based on FreeSurfer v6.0 to extract anatomical labels and volumes of interest in the cortex. We used analysis of covariance to investigate the differences in susceptibility among the clinical diagnostic groups in each brain region. Multivariable linear regression analysis was performed to study the association between susceptibility values and cognitive scores including the Mini-Mental State Examination (MMSE).

RESULTS

Among the three groups, the frontal (P < 0.001), temporal (P = 0.004), parietal (P = 0.001), occipital (P < 0.001), and cingulate cortices (P < 0.001) showed a higher mean susceptibility in patients with MCI and AD than in NC subjects. In the combined MCI and AD group, the mean susceptibility in the cingulate cortex (β = -216.21, P = 0.019) and insular cortex (β = -276.65, P = 0.001) were significant independent predictors of MMSE scores after correcting for age, sex, education, regional volume, and APOE4 carrier status.

CONCLUSION

Iron deposition in the cortex, as measured by QSMnet+, was higher in patients with AD and MCI than in NC participants. Iron deposition in the cingulate and insular cortices may be an early imaging marker of cognitive impairment related neurodegeneration.

DeepSTI: Towards tensor reconstruction using fewer orientations in susceptibility tensor imaging

Susceptibility tensor imaging (STI) is an emerging magnetic resonance imaging technique that characterizes the anisotropic tissue magnetic susceptibility with a second-order tensor model. STI has the potential to provide information for both the reconstruction of white matter fiber pathways and detection of myelin changes in the brain at mm resolution or less, which would be of great value for understanding brain structure and function in healthy and diseased brain. However, the application of STI in vivo has been hindered by its cumbersome and time-consuming acquisition requirement of measuring susceptibility induced MR phase changes at multiple head orientations. Usually, sampling at more than six orientations is required to obtain sufficient information for the ill-posed STI dipole inversion. This complexity is enhanced by the limitation in head rotation angles due to physical constraints of the head coil. As a result, STI has not yet been widely applied in human studies in vivo. In this work, we tackle these issues by proposing an image reconstruction algorithm for STI that leverages data-driven priors. Our method, called DeepSTI, learns the data prior implicitly via a deep neural network that approximates the proximal operator of a regularizer function for STI. The dipole inversion problem is then solved iteratively using the learned proximal network. Experimental results using both simulation and in vivo human data demonstrate great improvement over state-of-the-art algorithms in terms of the reconstructed tensor image, principal eigenvector maps and tractography results, while allowing for tensor reconstruction with MR phase measured at much less than six different orientations. Notably, promising reconstruction results are achieved by our method from only one orientation in human in vivo, and we demonstrate a potential application of this technique for estimating lesion susceptibility anisotropy in patients with multiple sclerosis.

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.

Depth-wise profiles of iron and myelin in the cortex and white matter using χ-separation: A preliminary study

The in-vivo profiling of iron and myelin across cortical depths and underlying white matter has important implications for advancing knowledge about their roles in brain development and degeneration. Here, we utilize χ-separation, a recently-proposed advanced susceptibility mapping that creates positive (χ_pos) and negative (χ_neg) susceptibility maps, to generate the depth-wise profiles of χ_pos and χ_neg as surrogate biomarkers for iron and myelin, respectively. Two regional sulcal fundi of precentral and middle frontal areas are profiled and compared with findings from previous studies. The results show that the χ_pos profiles peak at superificial white matter (SWM), which is an area beneath cortical gray matter known to have the highest accumulation of iron within the cortex and white matter. On the other hand, the χ_neg profiles increase in SWM toward deeper white matter. These characteristics in the two profiles are in agreement with histological findings of iron and myelin. Furthermore, the χ_neg profiles report regional differences that agree with well-known distributions of myelin concentration. When the two profiles are compared with those of QSM and R₂*, different shapes and peak locations are observed. This preliminary study offers an insight into one of the possible applications of χ-separation for exploring microstructural information of the human brain, as well as clinical applications in monitoring changes of iron and myelin in related diseases.

Affine transformation edited and refined deep neural network for quantitative susceptibility mapping

Deep neural networks have demonstrated great potential in solving dipole inversion for Quantitative Susceptibility Mapping (QSM). However, the performances of most existing deep learning methods drastically degrade with mismatched sequence parameters such as acquisition orientation and spatial resolution. We propose an end-to-end AFfine Transformation Edited and Refined (AFTER) deep neural network for QSM, which is robust against arbitrary acquisition orientation and spatial resolution up to 0.6 mm isotropic at the finest. The AFTER-QSM neural network starts with a forward affine transformation layer, followed by a Unet for dipole inversion, then an inverse affine transformation layer, followed by a Residual Dense Network (RDN) for QSM refinement. Simulation and in-vivo experiments demonstrated that the proposed AFTER-QSM network architecture had excellent generalizability. It can successfully reconstruct susceptibility maps from highly oblique and anisotropic scans, leading to the best image quality assessments in simulation tests and suppressed streaking artifacts and noise levels for in-vivo experiments compared with other methods. Furthermore, ablation studies showed that the RDN refinement network significantly reduced image blurring and susceptibility underestimation due to affine transformations. In addition, the AFTER-QSM network substantially shortened the reconstruction time from minutes using conventional methods to only a few seconds.

Rapid high-fidelity T 2 * mapping using single-shot overlapping-echo acquisition and deep learning reconstruction

PURPOSE

To develop and evaluate a single-shot quantitative MRI technique called GRE-MOLED (gradient-echo multiple overlapping-echo detachment) for rapid T 2 * $$ {T}_2^{\ast } $$ mapping.

METHODS

In GRE-MOLED, multiple echoes with different TEs are generated and captured in a single shot of the k-space through MOLED encoding and EPI readout. A deep neural network, trained by synthetic data, was employed for end-to-end parametric mapping from overlapping-echo signals. GRE-MOLED uses pure GRE acquisition with a single echo train to deliver T 2 * $$ {T}_2^{\ast } $$ maps less than 90 ms per slice. The self-registered B₀ information modulated in image phase was utilized for distortion-corrected parametric mapping. The proposed method was evaluated in phantoms, healthy volunteers, and task-based FMRI experiments.

RESULTS

The quantitative results of GRE-MOLED T 2 * $$ {T}_2^{\ast } $$ mapping demonstrated good agreement with those obtained from the multi-echo GRE method (Pearson's correlation coefficient = 0.991 and 0.973 for phantom and in vivo brains, respectively). High intrasubject repeatability (coefficient of variation <1.0%) were also achieved in scan-rescan test. Enabled by deep learning reconstruction, GRE-MOLED showed excellent robustness to geometric distortion, noise, and random subject motion. Compared to the conventional FMRI approach, GRE-MOLED also achieved a higher temporal SNR and BOLD sensitivity in task-based FMRI.

CONCLUSION

GRE-MOLED is a new real-time technique for T 2 * $$ {T}_2^{\ast } $$ quantification with high efficiency and quality, and it has the potential to be a better quantitative BOLD detection method.

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.

Incomplete spectrum QSM using support information

Introduction

Reconstructing a bounded object from incomplete k-space data is a well posed problem, and it was recently shown that this incomplete spectrum approach can be used to reconstruct undersampled MRI images with similar quality to compressed sensing approaches. Here, we apply this incomplete spectrum approach to the field-to-source inverse problem encountered in quantitative magnetic susceptibility mapping (QSM). The field-to-source problem is an ill-posed problem because of conical regions in frequency space where the dipole kernel is zero or very small, which leads to the kernel's inverse being ill-defined. These "ill-posed" regions typically lead to streaking artifacts in QSM reconstructions. In contrast to compressed sensing, our approach relies on knowledge of the image-space support, more commonly referred to as the mask, of our object as well as the region in k-space with ill-defined values. In the QSM case, this mask is usually available, as it is required for most QSM background field removal and reconstruction methods.

Methods

We tuned the incomplete spectrum method (mask and band-limit) for QSM on a simulated dataset from the most recent QSM challenge and validated the QSM reconstruction results on brain images acquired in five healthy volunteers, comparing incomplete spectrum QSM to current state-of-the art-methods: FANSI, nonlinear dipole inversion, and conventional thresholded k-space division.

Results

Without additional regularization, incomplete spectrum QSM performs slightly better than direct QSM reconstruction methods such as thresholded k-space division (PSNR of 39.9 vs. 39.4 of TKD on a simulated dataset) and provides susceptibility values in key iron-rich regions similar or slightly lower than state-of-the-art algorithms, but did not improve the PSNR in comparison to FANSI or nonlinear dipole inversion. With added (ℓ1-wavelet based) regularization the new approach produces results similar to compressed sensing based reconstructions (at sufficiently high levels of regularization).

Discussion

Incomplete spectrum QSM provides a new approach to handle the "ill-posed" regions in the frequency-space data input to QSM.

Quantitative susceptibility mapping (QSM) of the cardiovascular system: challenges and perspectives

Quantitative susceptibility mapping (QSM) is a powerful, non-invasive, magnetic resonance imaging (MRI) technique that relies on measurement of magnetic susceptibility. So far, QSM has been employed mostly to study neurological disorders characterized by iron accumulation, such as Parkinson's and Alzheimer's diseases. Nonetheless, QSM allows mapping key indicators of cardiac disease such as blood oxygenation and myocardial iron content. For this reason, the application of QSM offers an unprecedented opportunity to gain a better understanding of the pathophysiological changes associated with cardiovascular disease and to monitor their evolution and response to treatment. Recent studies on cardiovascular QSM have shown the feasibility of a non-invasive assessment of blood oxygenation, myocardial iron content and myocardial fibre orientation, as well as carotid plaque composition. Significant technical challenges remain, the most evident of which are related to cardiac and respiratory motion, blood flow, chemical shift effects and susceptibility artefacts. Significant work is ongoing to overcome these challenges and integrate the QSM technique into clinical practice in the cardiovascular field.

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.

In this article, A review of the results of articles on the use of QSM technique to identify nuclei exposed to iron deposition in various neurodegenerative diseases was performed.

Brain nuclei with the highest changes in iron deposition were identified as a biomarker for the identification of specific neurological diseases

By studying recent articles, The best toolbox for each step of the QSM processing algorithm was introduced.

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) χ₃₃ 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 χ₃₃-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 were also available at https://osf.io/y6rc3/.

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.

Deep learning-based quantitative susceptibility mapping (QSM) in the presence of fat using synthetically generated multi-echo phase training data

PURPOSE

To enable a fast and automatic deep learning-based QSM reconstruction of tissues with diverse chemical shifts, relevant to most regions outside the brain.

METHODS

A UNET was trained to reconstruct susceptibility maps using synthetically generated, unwrapped, multi-echo phase data as input. The RMS error with respect to synthetic validation data was computed. The method was tested on two in vivo knee and two pelvis data sets. Comparisons were made to a conventional fat-water separation pipeline by applying a commonly used graph-cut algorithm, both without and with an extended mask for background field removal (FWS-CONV-QSM and FWS-MASK-CONV-QSM, respectively). Several regions of interest were segmented and compared. Furthermore, the approach was tested on a prostate cancer patient receiving low-dose-rate brachytherapy, to detect and localize the seeds by MRI.

RESULTS

The RMS error was 0.292 ppm with FWS-CONV-QSM and 0.123 ppm for the UNET approach. Susceptibility maps were reconstructed much faster (< 10 s) and completely automatically (no background masking needed) by the UNET compared with the other applied techniques (5 min 51 s and 22 min 44 s for CONV-QSM and FWS-MASK-CONV-QSM, respectively. Background artifacts, fat-water swaps, and hypointense artifacts between I-125 seeds of a patient receiving low-dose brachytherapy in the prostate were largely reduced in the UNET approach.

CONCLUSIONS

Deep learning-based QSM reconstruction, trained solely with synthetic data, is well-suited to rapidly reconstructing high-quality susceptibility maps in the presence of fat without needing masking for background field removal.

Choroid Plexus Volume and Permeability at Brain MRI within the Alzheimer Disease Clinical Spectrum

Background Mounting evidence suggests that the choroid plexus (CP) plays an important role in the pathophysiology of Alzheimer disease (AD), but its imaging profile in cognitive impairment remains unclear. Purpose To evaluate CP volume, permeability, and susceptibility by using MRI in patients at various stages of cognitive impairment. Materials and Methods This retrospective study evaluated patients with cognitive symptoms who underwent 3.0-T MRI of the brain, including dynamic contrast-enhanced (DCE) imaging and quantitative susceptibility mapping (QSM), between January 2013 and May 2020. CP volume was automatically segmented using three-dimensional T1-weighted sequences; the volume transfer constant (ie, K^trans) and fractional plasma volume (ie, V_p) were determined using DCE MRI, and susceptibility was assessed using QSM. The effects of CP volume, expressed as the ratio to intracranial volume, on cognition were evaluated using multivariable linear regression adjusted for age, sex, education, apolipoprotein E ε4 allele status, and volumetric measures. Results A total of 532 patients with cognitive symptoms (mean age, 72 years ± 9 [SD]; 388 women) were included: 78 with subjective cognitive impairment (SCI), 158 with early mild cognitive impairment (MCI), 149 with late MCI, and 147 with AD. Among these, 132 patients underwent DCE MRI and QSM. CP volume was greater in patients at more severe stages (ratio of intracranial volume × 10³: 0.9 ± 0.3 for SCI, 1.0 ± 0.3 for early MCI, 1.1 ± 0.3 for late MCI, and 1.3 ± 0.4 for AD; P < .001). Lower K^trans (r = -0.19; P = .03) and V_p (r = -0.20; P = .02) were negatively associated with CP volume; susceptibility was not (r = 0.15; P = .10). CP volume was negatively associated with memory (B = -0.67; standard error of the mean [SEM], 0.21; P = .01), executive function (B = -0.90; SEM, 0.31; P = .01), and global cognition (B = -0.82; SEM, 0.32; P = .01). Conclusion Among patients with cognitive symptoms, larger choroid plexus volume was associated with severity of cognitive impairment in the Alzheimer disease spectrum. Published under a CC BY 4.0 license. Online supplemental material is available for this article. See also the editorial by Chiang in this issue.

[Brain Iron Imaging in Aging and Cognitive Disorders: MRI Approaches]

Iron has a vital role in the human body, including the central nervous system. Increased deposition of iron in the brain has been reported in aging and important neurodegenerative diseases. Owing to the unique magnetic resonance properties of iron, MRI has great potential for in vivo assessment of iron deposition, distribution, and non-invasive quantification. In this paper, we will review the MRI methods for iron assessment and their changes in aging and neurodegenerative diseases, focusing on Alzheimer's disease. In addition, we will summarize the limitations of current approaches and introduce new areas and MRI methods for iron imaging that are expected in the future.

Method to Minimize the Errors of AI: Quantifying and Exploiting Uncertainty of Deep Learning in Brain Tumor Segmentation

Despite the unprecedented success of deep learning in various fields, it has been recognized that clinical diagnosis requires extra caution when applying recent deep learning techniques because false prediction can result in severe consequences. In this study, we proposed a reliable deep learning framework that could minimize incorrect segmentation by quantifying and exploiting uncertainty measures. The proposed framework demonstrated the effectiveness of a public dataset: Multimodal Brain Tumor Segmentation Challenge 2018. By using this framework, segmentation performances, particularly for small lesions, were improved. Since the segmentation of small lesions is difficult but also clinically significant, this framework could be effectively applied to the medical imaging field.

Basal Ganglia Iron Content Increases with Glioma Severity Using Quantitative Susceptibility Mapping: A Potential Biomarker of Tumor Severity

BACKGROUND AND PURPOSE

Gliomas have been found to alter iron metabolism and transport in ways that result in an expansion of their intracellular iron compartments to support aggressive tumor growth. This study used deep neural network trained quantitative susceptibility mapping to assess basal ganglia iron concentrations in glioma patients.

MATERIALS AND METHODS

Ninety-two patients with brain lesions were initially enrolled in this study and fifty-nine met the inclusion criteria. Susceptibility-weighted images were collected at 3.0 T and used to construct quantitative susceptibility maps via a deep neural network-based method. The regions of interest were manually drawn within basal ganglia structures and the mean voxel intensities were extracted and averaged across multiple slices. One-way ANCOVA tests were conducted to compare the susceptibility values of groups of patients based on tumor grade while controlling for age, sex, and tumor type.

RESULTS

The mean basal ganglia susceptibility for patients with grade IV tumors was higher than that for patients with grade II tumors (p = 0.00153) and was also higher for patients with grade III tumors compared to patients with grade II tumors (p = 0.020), after controlling for age, sex, and tumor type. Patient age influenced susceptibility values (p = 0.00356), while sex (p = 0.69) and tumor type (p = 0.11) did not.

CONCLUSIONS

The basal ganglia iron content increased with glioma severity. Basal ganglia iron levels may thus be a useful biomarker in glioma prognosis and treatment, especially with regard to iron-based cancer therapies.

Learn Less, Infer More: Learning in the Fourier Domain for Quantitative Susceptibility Mapping

Quantitative susceptibility mapping (QSM) aims to evaluate the distribution of magnetic susceptibility from magnetic resonance phase measurements by solving the ill-conditioned dipole inversion problem. Removing the artifacts and preserving the anisotropy of tissue susceptibility simultaneously is still a challenge in QSM. To deal with this issue, a novel k-QSM network is proposed to resolve dipole inversion issues in QSM reconstruction. The k-QSM network converts the results obtained by truncated k-space division (TKD) into the Fourier domain as inputs. After passing through several convolutional and residual blocks, the ill-posed signals of TKD are corrected by making the network output close to the calculation of susceptibility through multiple orientation sampling (COSMOS)-labeled QSM. To evaluate the superiority of k-QSM, comparisons with several state-of-the-art methods are performed in terms of QSM artifacts removing, anisotropy preserving, generalization ability, and clinical applications. Compared to existing methods, the k-QSM achieves a 22.31% lower normalized root mean square error, 10.30% higher peak signal-to-noise ratio (PSNR), 33.10% lower high-frequency error norm, and 1.06% higher structural similarity. In addition, the orientation-dependent susceptibility variation obtained by k-QSM is significant, verifying that k-QSM has the ability to preserve susceptibility anisotropy. When the trained models are tested on the dataset from different centers, our k-QSM shows a strong generalization ability with the highest PSNR. Moreover, by comparing the susceptibility maps between healthy controls and drug addicts with different methods, we found the proposed k-QSM is more sensitive to the susceptibility abnormality in the patients. The proposed k-QSM method learns less—only to fix the ill-posed signals of TKD, but infers more—both COSMOS-like and anisotropy-preserving QSM results. Its generalization ability and great sensitivity to susceptibility changes can make it a potential method for distinguishing some diseases.

DIFFnet: Diffusion Parameter Mapping Network Generalized for Input Diffusion Gradient Schemes and b-Value

In MRI, deep neural networks have been proposed to reconstruct diffusion model parameters. However, the inputs of the networks were designed for a specific diffusion gradient scheme (i.e., diffusion gradient directions and numbers) and a specific b-value that are the same as the training data. In this study, a new deep neural network, referred to as DIFFnet, is developed to function as a generalized reconstruction tool of the diffusion-weighted signals for various gradient schemes and b-values. For generalization, diffusion signals are normalized in a q-space and then projected and quantized, producing a matrix (Qmatrix) as an input for the network. To demonstrate the validity of this approach, DIFFnet is evaluated for diffusion tensor imaging (DIFFnet_DTI) and for neurite orientation dispersion and density imaging (DIFFnet_NODDI). In each model, two datasets with different gradient schemes and b-values are tested. The results demonstrate accurate reconstruction of the diffusion parameters at substantially reduced processing time (approximately 8.7 times and 2240 times faster processing time than conventional methods in DTI and NODDI, respectively; less than 4% mean normalized root-mean-square errors (NRMSE) in DTI and less than 8% in NODDI). The generalization capability of the networks was further validated using reduced numbers of diffusion signals from the datasets and a public dataset from Human Connection Project. Different from previously proposed deep neural networks, DIFFnet does not require any specific gradient scheme and b-value for its input. As a result, it can be adopted as an online reconstruction tool for various complex diffusion imaging.

Accelerating quantitative susceptibility and R2* mapping using incoherent undersampling and deep neural network reconstruction

Quantitative susceptibility mapping (QSM) and R2* mapping are MRI post-processing methods that quantify tissue magnetic susceptibility and transverse relaxation rate distributions. However, QSM and R2* acquisitions are relatively slow, even with parallel imaging. Incoherent undersampling and compressed sensing reconstruction techniques have been used to accelerate traditional magnitude-based MRI acquisitions; however, most do not recover the full phase signal, as required by QSM, due to its non-convex nature. In this study, a learning-based Deep Complex Residual Network (DCRNet) is proposed to recover both the magnitude and phase images from incoherently undersampled data, enabling high acceleration of QSM and R2* acquisition. Magnitude, phase, R2*, and QSM results from DCRNet were compared with two iterative and one deep learning methods on retrospectively undersampled acquisitions from six healthy volunteers, one intracranial hemorrhage and one multiple sclerosis patients, as well as one prospectively undersampled healthy subject using a 7T scanner. Peak signal to noise ratio (PSNR), structural similarity (SSIM), root-mean-squared error (RMSE), and region-of-interest susceptibility and R2* measurements are reported for numerical comparisons. The proposed DCRNet method substantially reduced artifacts and blurring compared to the other methods and resulted in the highest PSNR, SSIM, and RMSE on the magnitude, R2*, local field, and susceptibility maps. Compared to two iterative and one deep learning methods, the DCRNet method demonstrated a 3.2% to 9.1% accuracy improvement in deep grey matter susceptibility when accelerated by a factor of four. The DCRNet also dramatically shortened the reconstruction time of single 2D brain images from 36-140 seconds using conventional approaches to only 15-70 milliseconds.

MoDL-QSM: Model-based deep learning for quantitative susceptibility mapping

Quantitative susceptibility mapping (QSM) has demonstrated great potential in quantifying tissue susceptibility in various brain diseases. However, the intrinsic ill-posed inverse problem relating the tissue phase to the underlying susceptibility distribution affects the accuracy for quantifying tissue susceptibility. Recently, deep learning has shown promising results to improve accuracy by reducing the streaking artifacts. However, there exists a mismatch between the observed phase and the theoretical forward phase estimated by the susceptibility label. In this study, we proposed a model-based deep learning architecture that followed the STI (susceptibility tensor imaging) physical model, referred to as MoDL-QSM. Specifically, MoDL-QSM accounts for the relationship between STI-derived phase contrast induced by the susceptibility tensor terms (χ₁₃, χ₂₃ and χ₃₃) and the acquired single-orientation phase. The convolutional neural networks are embedded into the physical model to learn a regularization term containing prior information. χ₃₃ and phase induced by χ₁₃ and χ₂₃ terms were used as the labels for network training. Quantitative evaluation metrics were compared with recently developed deep learning QSM methods. The results showed that MoDL-QSM achieved superior performance, demonstrating its potential for future applications.

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.

The Calcium Versus Hemorrhage Trial: Developing Diagnostic Criteria for Chronic Intracranial Susceptibility Lesions Using Single-Energy Computed Tomography, Dual-Energy Computed Tomography, and Quantitative Susceptibility Mapping

PURPOSE

Chronic susceptibility lesions in the brain can be either hemorrhagic (potentially dangerous) or calcific (usually not dangerous) but are difficult to discriminate on routine imaging. We proposed to develop quantitative diagnostic criteria for single-energy computed tomography (SECT), dual-energy computed tomography (DECT), and quantitative susceptibility mapping (QSM) to distinguish hemorrhage from calcium.

MATERIALS AND METHODS

Patients with positive susceptibility lesions on routine T2*-weighted magnetic resonance of the brain were recruited into this prospective imaging clinical trial, under institutional review board approval and with informed consent. The SECT, DECT, and QSM images were obtained, the lesions were identified, and the regions of interest were defined, with the mean values recorded. Criteria for quantitative interpretation were developed on the first 50 patients, and then applied to the next 45 patients. Contingency tables, scatter plots, and McNemar test were applied to compare classifiers.

RESULTS

There were 95 evaluable patients, divided into a training set of 50 patients (328 lesions) and a validation set of 45 patients (281 lesions). We found the following classifiers to best differentiate hemorrhagic from calcific lesions: less than 68 Hounsfield units for SECT, calcium level of less than 15 mg/mL (material decomposition value) for DECT, and greater than 38 ppb for QSM. There was general mutual agreement among the proposed criteria. The proposed criteria outperformed the current published criteria.

CONCLUSIONS

We provide the updated criteria for the classification of chronic positive susceptibility brain lesions as hemorrhagic versus calcific for each major clinically available imaging modality. These proposed criteria have greater internal consistency than the current criteria and should likely replace it as gold standard.

Brain/MINDS beyond human brain MRI project: A protocol for multi-level harmonization across brain disorders throughout the lifespan

Psychiatric and neurological disorders are afflictions of the brain that can affect individuals throughout their lifespan. Many brain magnetic resonance imaging (MRI) studies have been conducted; however, imaging-based biomarkers are not yet well established for diagnostic and therapeutic use. This article describes an outline of the planned study, the Brain/MINDS Beyond human brain MRI project (BMB-HBM, FY2018 ~ FY2023), which aims to establish clinically-relevant imaging biomarkers with multi-site harmonization by collecting data from healthy traveling subjects (TS) at 13 research sites. Collection of data in psychiatric and neurological disorders across the lifespan is also scheduled at 13 sites, whereas designing measurement procedures, developing and analyzing neuroimaging protocols, and databasing are done at three research sites. A high-quality scanning protocol, Harmonization Protocol (HARP), was established for five high-quality 3 T scanners to obtain multimodal brain images including T1 and T2-weighted, resting-state and task functional and diffusion-weighted MRI. Data are preprocessed and analyzed using approaches developed by the Human Connectome Project. Preliminary results in 30 TS demonstrated cortical thickness, myelin, functional connectivity measures are comparable across 5 scanners, suggesting sensitivity to subject-specific connectome. A total of 75 TS and more than two thousand patients with various psychiatric and neurological disorders are scheduled to participate in the project, allowing a mixed model statistical harmonization. The HARP protocols are publicly available online, and all the imaging, demographic and clinical information, harmonizing database will also be made available by 2024. To the best of our knowledge, this is the first project to implement a prospective, multi-level harmonization protocol with multi-site TS data. It explores intractable brain disorders across the lifespan and may help to identify the disease-specific pathophysiology and imaging biomarkers for clinical practice.

xQSM: quantitative susceptibility mapping with octave convolutional and noise-regularized neural networks

Quantitative susceptibility mapping (QSM) provides a valuable MRI contrast mechanism that has demonstrated broad clinical applications. However, the image reconstruction of QSM is challenging due to its ill-posed dipole inversion process. In this study, a new deep learning method for QSM reconstruction, namely xQSM, was designed by introducing noise regularization and modified octave convolutional layers into a U-net backbone and trained with synthetic and in vivo datasets, respectively. The xQSM method was compared with two recent deep learning (QSMnet⁺ and DeepQSM) and two conventional dipole inversion (MEDI and iLSQR) methods, using both digital simulations and in vivo experiments. Reconstruction error metrics, including peak signal-to-noise ratio, structural similarity, normalized root mean squared error and deep gray matter susceptibility measurements, were evaluated for comparison of the different methods. The results showed that the proposed xQSM network trained with in vivo datasets achieved the best reconstructions of all the deep learning methods. In particular, it led to, on average, 32.3%, 25.4% and 11.7% improvement in the accuracy of globus pallidus susceptibility estimation for digital simulations and 39.3%, 21.8% and 6.3% improvements for in vivo acquisitions compared with DeepQSM, QSMnet⁺ and iLSQR, respectively. It also exhibited the highest linearity against different susceptibility intensity scales and demonstrated the most robust generalization capability to various spatial resolutions of all the deep learning methods. In addition, the xQSM method also substantially shortened the reconstruction time from minutes using MEDI to only a few seconds.

Deep model-based magnetic resonance parameter mapping network (DOPAMINE) for fast T1 mapping using variable flip angle method

Quantitative tissue characteristics, which provide valuable diagnostic information, can be represented by magnetic resonance (MR) parameter maps using magnetic resonance imaging (MRI); however, a long scan time is necessary to acquire them, which prevents the application of quantitative MR parameter mapping to real clinical protocols. For fast MR parameter mapping, we propose a deep model-based MR parameter mapping network called DOPAMINE that combines a deep learning network with a model-based method to reconstruct MR parameter maps from undersampled multi-channel k-space data. DOPAMINE consists of two networks: 1) an MR parameter mapping network that uses a deep convolutional neural network (CNN) that estimates initial parameter maps from undersampled k-space data (CNN-based mapping), and 2) a reconstruction network that removes aliasing artifacts in the parameter maps with a deep CNN (CNN-based reconstruction) and an interleaved data consistency layer by an embedded MR model-based optimization procedure. We demonstrated the performance of DOPAMINE in brain T₁ map reconstruction with a variable flip angle (VFA) model. To evaluate the performance of DOPAMINE, we compared it with conventional parallel imaging, low-rank based reconstruction, model-based reconstruction, and state-of-the-art deep-learning-based mapping methods for three different reduction factors (R = 3, 5, and 7) and two different sampling patterns (1D Cartesian and 2D Poisson-disk). Quantitative metrics indicated that DOPAMINE outperformed other methods in reconstructing T₁ maps for all sampling patterns and reduction factors. DOPAMINE exhibited quantitatively and qualitatively superior performance to that of conventional methods in reconstructing MR parameter maps from undersampled multi-channel k-space data. The proposed method can thus reduce the scan time of quantitative MR parameter mapping that uses a VFA model.

SEPIA-Susceptibility mapping pipeline tool for phase images

Quantitative susceptibility mapping (QSM) is a physics-driven computational technique that has a high sensitivity in quantifying iron deposition based on MRI phase images. Furthermore, it has a unique ability to distinguish paramagnetic and diamagnetic contributions such as haemorrhage and calcification based on image contrast. These properties have contributed to a growing interest to use QSM not only in research but also in clinical applications. However, it is challenging to obtain high quality susceptibility map because of its ill-posed nature, especially for researchers who have less experience with QSM and the optimisation of its pipeline. In this paper, we present an open-source processing pipeline tool called SuscEptibility mapping PIpeline tool for phAse images (SEPIA) dedicated to the post-processing of MRI phase images and QSM. SEPIA connects various QSM toolboxes freely available in the field to offer greater flexibility in QSM processing. It also provides an interactive graphical user interface to construct and execute a QSM processing pipeline, simplifying the workflow in QSM research. The extendable design of SEPIA also allows developers to deploy their methods in the framework, providing a platform for developers and researchers to share and utilise the state-of-the-art methods in QSM.

DeepResp: Deep learning solution for respiration-induced B0 fluctuation artifacts in multi-slice GRE

Respiration-induced B₀ fluctuation corrupts MRI images by inducing phase errors in k-space. A few approaches such as navigator have been proposed to correct for the artifacts at the expense of sequence modification. In this study, a new deep learning method, which is referred to as DeepResp, is proposed for reducing the respiration-artifacts in multi-slice gradient echo (GRE) images. DeepResp is designed to extract the respiration-induced phase errors from a complex image using deep neural networks. Then, the network-generated phase errors are applied to the k-space data, creating an artifact-corrected image. For network training, the computer-simulated images were generated using artifact-free images and respiration data. When evaluated, both simulated images and in-vivo images of two different breathing conditions (deep breathing and natural breathing) show improvements (simulation: normalized root-mean-square error (NRMSE) from 7.8 ± 5.2% to 1.3 ± 0.6%; structural similarity (SSIM) from 0.88 ± 0.08 to 0.99 ± 0.01; ghost-to-signal-ratio (GSR) from 7.9 ± 7.2% to 0.6 ± 0.6%; deep breathing: NRMSE from 13.9 ± 4.6% to 5.8 ± 1.4%; SSIM from 0.86 ± 0.03 to 0.95 ± 0.01; GSR 20.2 ± 10.2% to 5.7 ± 2.3%; natural breathing: NRMSE from 5.2 ± 3.3% to 4.0 ± 2.5%; SSIM from 0.94 ± 0.04 to 0.97 ± 0.02; GSR 5.7 ± 5.0% to 2.8 ± 1.1%). Our approach does not require any modification of the sequence or additional hardware, and may therefore find useful applications. Furthermore, the deep neural networks extract respiration-induced phase errors, which is more interpretable and reliable than results of end-to-end trained networks.

Learned Proximal Networks for Quantitative Susceptibility Mapping

Quantitative Susceptibility Mapping (QSM) estimates tissue magnetic susceptibility distributions from Magnetic Resonance (MR) phase measurements by solving an ill-posed dipole inversion problem. Conventional single orientation QSM methods usually employ regularization strategies to stabilize such inversion, but may suffer from streaking artifacts or over-smoothing. Multiple orientation QSM such as calculation of susceptibility through multiple orientation sampling (COSMOS) can give well-conditioned inversion and an artifact free solution but has expensive acquisition costs. On the other hand, Convolutional Neural Networks (CNN) show great potential for medical image reconstruction, albeit often with limited interpretability. Here, we present a Learned Proximal Convolutional Neural Network (LP-CNN) for solving the ill-posed QSM dipole inversion problem in an iterative proximal gradient descent fashion. This approach combines the strengths of data-driven restoration priors and the clear interpretability of iterative solvers that can take into account the physical model of dipole convolution. During training, our LP-CNN learns an implicit regularizer via its proximal, enabling the decoupling between the forward operator and the data-driven parameters in the reconstruction algorithm. More importantly, this framework is believed to be the first deep learning QSM approach that can naturally handle an arbitrary number of phase input measurements without the need for any ad-hoc rotation or re-training. We demonstrate that the LP-CNN provides state-of-the-art reconstruction results compared to both traditional and deep learning methods while allowing for more flexibility in the reconstruction process.