MRI Tissue Magnetism Quantification through Total Field Inversion with Deep Neural Networks

04/11/2019 ∙ by Juan Liu, et al. ∙ Medical College of Wisconsin 4

Quantitative susceptibility mapping (QSM) utilizes MRI signal phase to infer estimates of local tissue magnetism (magnetic susceptibility), which has been shown useful to provide novel image contrast and as biomarkers of abnormal tissue. QSM requires addressing a challenging post-processing problem: filtering of image phase estimates and inversion of the phase to susceptibility relationship. A wide variety of quantification errors, robustness limitations, and artifacts plague QSM algorithms. To overcome these limitations, a robust deep-learning-based single-step QSM reconstruction approach is proposed and demonstrated. This neural network was trained using magnetostatic physics simulations based on in-vivo data sources. Random perturbations were added to the physics simulations to provide sufficient input-label pairs for the training purposes. The network was quantitatively tested using gold-standard in-silico labeled datasets against established QSM total field inversion approaches. In addition, the algorithm was applied to susceptibility-weighted imaging (SWI) data collected on a cohort of clinical subjects with brain hemmhorage. When quantitatively compared against gold-standard in-silico labels, the proposed algorithm outperformed the existing comparable approaches. High quality QSM were consistently estimated from clinical susceptibility-weighted data on 100 subjects without any noticeable inversion failures. The proposed approach was able to robustly generate high quality QSM with improved accuracy in in-silico gold-standard experiments. QSM produced by the proposed method can be generated in real-time on existing MRI scanner platforms and provide enhanced visualization and quantification of magnetism-based tissue contrasts.

READ FULL TEXT VIEW PDF
POST COMMENT

Comments

There are no comments yet.

Authors

page 12

page 14

page 15

page 16

page 17

page 18

page 19

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

I. Introduction

Magnetic susceptibility, which is a parametric measure of material magnetism 1, has long been utilized to generate diagnostic imaging contrast through magnetic resonance imaging (MRI) techniques such as T or enhanced susceptibility-weighted imaging (SWI) 2. Quantitative Susceptibility Mapping (QSM) is a MRI post-processing technique that expands upon these widely utilized susceptibility-weighted practices to provide estimates of underlying tissue magnetism 3.

SWI is widely-used for clinical assessments, including traumatic brain injury4, 5, vascular malformations 6, 7, brain tumors 8, 9, neurodegenerative diseases 10, 11, etc. Despite its widespread utility and usage, SWI suffers from blooming intensity artifacts and cannot differentiate calcifications from small microbleeds or iron deposits. Besides, SWI contrasts change with field strength, echo time settings, and depend on both the shape and the orientation of the blood vessels. By comparison, QSM quantitatively estimating the underlying tissue magnetic susceptibility is able to overcome the limitations of SWI. Existing studies have explored to usage of QSM to quantification of specific biomarkers including iron, calcium, gadolinium etc 12, 13, 14, 15. QSM has been explored to study brain tumors16, 17, neurodegenerative disorders and multiple sclerosis18, 19, 17, iron overload in the liver20, blood oxygenation assessment21, and mild traumatic brain injury22 as well.

QSM is performed by collecting phase-sensitive MRI, estimating local magnetic perturbations along the direction of the polarizing B magnetic field, and then using magnetic perturbations to estimate source magnetic susceptibility (i.e. magnetism) within tissue 3. The collection of phase-sensitive MRI and associated steps to develop local field perturbation estimates is well-established. On the contrary, the transformation of raw local field measurements to tissue susceptibility estimates is fraught with technical challenges. Historically, QSM transformation from raw perturbation fields usually consist of two steps: (i) background field removal to determine the local tissue field, and (ii) inversion from the local field to the tissue susceptibility. Both steps require solving ill-posed inverse computational problems.

Background field removal (BFR) methods typically rely on spatial filtering and/or dipole field modelling methods 23, 24, 25, 26. Existing BFR approaches can often perform well under ideal circumstances with perfect brain tissue masking. However, in the presence of these confounds, BFR is often compromised with residual background field and image artifacts in the local field, which then leads to artifacts and inaccuracies in final QSM estimates.

To avoid BFR error propagation into QSM estimation, several single-step QSM algorithms, which directly estimate the susceptibility distribution from the raw perturbation field, have been proposed. These approaches, hereby denoted as “total field inversion” (TFI) methods, utilize algorithmic adjustments such as Laplacian spatial conditioning operators 27, 28, 29, 30, preconditioning boosted iterative inversions 31, or total variation (TV) and/or Tikhonov regularized inversion 32, 33 to address the challenges imposed by TFI.

TFI methods have shown improved susceptibility estimation performance compared to conventional two-step QSM methods. However, all existing TFI methods require tuning of regularization parameters, which introduce estimation biases and artifacts. In addition, these methods require substantial computation times to solve their respected iterative optimization problems. Finally, our investigations have found that these methods still exhibit residual streaking artifacts and BFR failures in regions with large susceptibility variations, such as intracranial hemorrhages. Furthermore, we have found existing available methods to have compromised performance when applied to routinely-acquired clinical susceptibility-weighted imaging data, which is highly non-isotropic in resolution (due to examination time constraints).

To address the limitations of existing TFI methods, here we demonstrate a deep-learning-based TFI approach. Leveraging the well-defined forward susceptibility-to-field physics model 34, 35

, synthetic data are used for neural network training. This addresses the common machine-learning challenge of collecting and curating accurate ground-truth training data. The proposed TFI deep learning network is demonstrated on source MRI data collected under two different use cases with substantially different image resolutions. Quantitative in-silico gold-standard performance evaluations of the proposed method and several existing TFI and multi-step QSM algorithms are also presented. In addition, several qualitative demonstrations are presented to show superior performance of the proposed method relative to existing approaches.

Ii. Methods

Ii.A. Training Data

Multiple-Orientation QSM datasets, such as Calculation of Susceptibility through Multiple Orientation Sampling (COSMOS) 36

or Susceptibility Tensor Imaging (STI) 

37, which are not clinically viable acquisition approaches, are often treated as QSM golden-standard estimates. While these approaches do address concerns of artifacts and quantification concerns with the field inversion step, there are several concerns to consider when using such data for QSM deep learning training labels. First, the use of these approaches is assuming perfect background field removal of the training data. Second, these techniques are iterative regularized algorithms that are not guaranteed to provide exact estimates of the field-to-source magnetism relationship (i.e the “input” to “label” relationship is not a true “gold standard”). Rather than relying on such acquisitions for neural network training, we propose that the well-defined forward source-to-field model offers a more flexible and accurate gold-standard training label. In our approach, we utilize cohorts of in-vivo QSM estimates as an input to a computational training data generator, which allows for trivial data augmentation practices to provide large numbers of realistic training inputs with true gold-standard labels accurate to levels of machine precision. In this approach, the accuracy of the original in-vivo QSM estimate is largely irrelevant, as the network is being trained exclusively by the ensuing closed-form forward source-to-field computation.

200 in vivo QSM datasets were used as inputs to the simulated neural network training data generator. The resolution of this source data was 0.5x0.5x2.0 mm. From the four echo time images in each data set, QSM estimates were generated using the following existing tools: brain masking using SPM 38, BFR using the Regularization-enabled Sophisticated Harmonic Artifact Reduction on Phase data (RESHARP) 24 method, and susceptibility inversion was performed using a previously developed Approximated Susceptibility through Parcellated Encoder-decoder Networks (ASPEN) 39.

Besides, geometric shapes such as ellipsoid, sphere, cuboid and cylinder with random susceptibility values and random orientations are randomly placed on the susceptibility maps. Using the in-vivo QSM estimates, local tissue magnetic perturbations were calculated using well-known dipole convolution methods 35, 34. The background perturbation fields were then simulated from random magnetic susceptibility sources to mimic background field. The field perturbation for training total field inversion was then constructed from the superimposed local tissue perturbations and the background perturbation. Fig. 1 provides an example input (total field perturbations) and label (tissue susceptibility) from this training data construction process.

Ii.B. Neural Network Architecture and Training

A 3D convolutional neural network with encoder-decoder architecture was trained to perform total field QSM inversion with 3D unwrapped total field maps and brain masks as the inputs. In vanilla convolutional neural networks, the spatial invariant kernels is inappropriate to solve our problem due to large background fields close to brain boundary and invalid data outside brain mask. To overcome these limitations, gated convolution with LeakyReLU as activation function and Sigmoid for gating value was applied for automatically selecting adaptive features for each channel and each spatial location. Dilated gated convolution was applied in deeper convolutional layers to increase the receptive fields. Moreover, a non-local block was inserted in encoder-decoder network to capture long-range dependencies for non-local susceptibility estimation

40.

For performance evaluation using synthetic data and QSM challenge, the image resolution is 1.06 mm isotropic with neural network input size 160x160x160. For clinical data, the image resolution is 0.76x0.76x3.0 mm

with size 256x256x64. 5000 data were used for training. L1 loss between the label and output was utilized as a loss function. The RMSprop optimizer was used in the deep learning training. The initial learning rate was set as 0.0001, with exponential decay at every 200 steps. Two NVIDIA tesla k40 graphics processing units (GPUs) were used for training with batch size 2. The neural network was trained and evaluated using Keras with Tensorflow as backend.

Ii.C. Performance Evaluation

Ii.C.1. Synthetic Data

100 simulated data sets generated in similar fashion to the training data without containing randomly inserted geometric shapes were used to evaluate the performance of QSMTFINet compared with both single-step and two-steps QSM reconstruction approaches, were selected for comparison, including (i) single step total variation QSM (SS-TV-QSM) 32, (ii) least-norm QSM (LN-QSM) 33, (iii) RESHARP and Thresholded K-space Division (TKD), (iv) Projection onto Dipole Fields (PDF) 25 and Morphology Enabled Dipole Inversion (MEDI) 41.

Estimation errors from each technique were computed using root mean squared error (RMSE), high-frequency error norm (HFEN), and structural similarity (SSIM) index with the ground truth susceptibility maps.

Ii.C.2. 2016 QSM Reconstruction Challenge Dataset

QSM reconstruction challenge dataset was acquired from a healthy 30-year-old female subject. Data for COSMOS and susceptibility-tensor computation was computed using a heavily accelerated wave-CAIPI 42 acquisition at 1.06 mm isotropic resolution collected at 12 different head orientations at a single echo time 43.

Using this gold-standard evaluation dataset, the proposed method was compared to SS-TV-QSM, LN-QSM, TKD, and MEDI approaches. TKD and MEDI were performed using the provided local tissue field. TKD results were provided publically by the QSM challenge organizers, with threshold 0.19 which yields the best trade-off between quantification accuracy and artifacts. All methods were evaluated against the “gold standard” STI (3,3) component computational result provided with the challenge dataset 37.

Ii.C.3. Clinical Data

One hundred clinical QSM data were acquired using gradient echo T weighted angiography (SWAN, GE), a new method for SWI with short acquisition times, at a 3T MRI scanner (GE Healthcare MR750) with data acquisition parameters: in-plane data acquisition matrix 288x224, field of view 22 cm, slice thickness 3 mm, autocalibrated parallel imaging factors 2x1, number of slices 46-54, first echo time 12.6 ms, echo spacing 4.1 ms, number of echoes 7, flip angle 15, repetition time 39.7 ms, total scan time about 2 minutes.

The SWI images were processed by vendor reconstruction algorithms. The raw k-space data were saved for offline QSM processing. Multi-echo real and imaginary data were reconstructed from k-space data, with reconstruction matrix size 288x288, voxel size 0.76x0.76x3.0 mm. The field map was obtained by the fitting of multi-echo phases. Brain masks were obtained using the FSL brain extraction tool.

Using the total field map and brain mask, representative state-of-the-arts methods, both two-steps QSM reconstruction and single-step QSM reconstruction, were selected for comparison, including the five QSM reconstruction methods mentioned above. Three clinical cases were chosen to demonstrate TFINet performance.

For the purposes of performance evaluation, with RESHARP, spherical kernel radius was set as 6mm to trade off the background removal performance and brain erosion; for TKD, the threshold was set to 0.20; for MEDI, the regularization factor was set to 1000, which can yield a good trade-off between quantification accuracy and artifacts; the TV and Tikhonov regularization parameters for LN-QSM were set to 4x10 and 10 respectively. MEDI toolbox, SS-TV-QSM code, LN-QSM code publicly provided by the authors were used to calculate the QSM images44, 45, 46.

Iii. Results

Iii..1. Synthetic Data

Table 1 illustrates the RMSE, HFEN, and SSIM using five QSM reconstruction methods compared with ground truth from 100 simulated datasets. The proposed method achieved the best score in RMSE, HFEN, and SSIM. In Fig.3, QSM images reconstructed by five different reconstruction methods are compared with ground truth susceptibility maps. The residual maps compared with ground truth are shown in Fig.4. Based on visual assessment, TFICnn results have best susceptibility estimations and least residual errors.

Iii..2. 2016 QSM Reconstruction Challenge Dataset

Fig. 5 illustrates susceptibility maps from the 2016 ISMRM QSM Challenge dataset reconstructed using five methods, displayed in 3 reformatted planes(rows i, ii, iii). Streaking artifacts are clearly identified in the sagittal reformat for SS-TV-QSM maps. In addition, the SS-TV-QSM map shows clearly compromised spatial resolution relative to the other maps. These observations are amplified in the zoomed maps in the last two rows (iv, v), where clear performance improvements of the proposed method in reproducing the fine structure of the STI (3,3) map and no-visible streaking artifacts.

Iii..3. Clinical Data

In Fig. 6, QSM results on one 65-year-old patient with stroke is shown in three (axial/coronal/sagittal) views. In zoom-in axial (ii), a few hypointense regions of SWI image (black arrows) are hyperintense in QSM, indicating it iron deposition or hemorrhage. In SS-TV-QSM results, obvious image blurring is shown in axial plane and apparent streaking artifacts in coronal and saggital planes, as shown (c, i-iv). Though LN-QSM results show good image sharpness and non-visible streaking artifacts, it suffers image intensity variation across the brain volume and shading artifacts near the brain boundary, as shown (d, i-iv). The QSM images of PDF+MEDI show good image quality, yet with visible image blurring and streaking artifacts, as shown (e, i-iv). RESHARP+TKD results suffer from brain erosion, image blurring, and streaking artifacts, as shown (f, i-iv). Compared with other methods, TFINet results (f, i-iv) show fine details and negligible streaking artifacts.

In Fig. 7, QSM and SWI images of a 56-year-old subject with hemorrhagic intracranial metastases is illustrated. SS-TV-QSM, PDF+MEDI, and RESHARP+TKD suffers severe shading artifacts and streaking artifacts, as shown (b, d, e, i-iv). LN-QSM shows strong shading artifacts, as shown (c, i-iv). With comparison to other methods, TFINet results (f, i-iv) showed improved image sharpness, clear tissue structures, and no shading and streaking artifacts around bleeding regions. Besides, two small calcification is dark/hypointense on SWI images and diamagnetic on TFINet QSM images (white dash arrows).

In Fig. 8, QSM results and SWI images are shown from a 37-year-old subject with surgical planning in three (axial/coronal/sagittal) views. Close to the bleeding regions, the SS-TV-QSM images show severe image blurring and shading artifacts in axial plane and severe streaking artifacts in coronal/saggital planes, as shown (b, i-iv). LN-QSM results clearly show black shading artifacts close to the intracerebral hemorrhage region and brain boundaries, as shown (c, i-iv). PDF+MEDI and RESHARP+TKD results show blurring in axial view, and apparent streaking artifacts in coronal/sagittal views, as shown (d, e, i-iv). Compared with all other methods, TFINet can produce high quality QSM images with fine details and non-visible streaking artifacts, shown in (f, i-iv).

In Fig. 9, QSM results and SWI images are shown from from a 54-year-old subject with poststereotactic radiosurgery (SRS) brain metastasis. Close to the bleeding regions, SS-TV-QSM, LN-QSM, PDF+MEDI, and RESHARP+TKD results (b-e) shows severe shading and bleeding artifacts. Compared with all other methods, TFINet can produce high quality QSM images with fine details and non-visible streaking artifacts, shown in (f, i-iv).

Iv. Discussion

In this work, TFINet was evaluated on synthetic data, QSM challenge dataset, and clinical dataset. Quantitative metrics and visual assessment demonstrated that TFINet outperform existing single-step QSM inversion and two-steps QSM inversion approaches. TFINet results have improved image sharpness, well-preserved microstructures, and non-visible streaking artifacts.

With comparison with single-step QSM methods SS-TV-QSM and LN-QSM, TFINet can produce high quality QSM images with better image sharpness and invisible image artifacts. SS-TV-QSM and LN-QSM use regularization methods to solve the ill-posed inverse problem, which require careful regularization parameter tuning and long iterative computation times. Though they eliminate the background field removal step, the QSM images suffer from over-smoothing, streaking artifacts, shading artifacts, and large susceptibility quantification errors close to intracerebral hemorrhage regions. The results show that TFINet can overcome the limitations and greatly improve susceptibility quantification accuracy.

Compared with two-steps QSM methods, TFINet can not only speed up the QSM processing but also eliminate the background field removal error propagation into susceptibility estimation. Besides, the proposed methods overcome the constraints of existing background removal approaches, such as brain erosion, parameter tuning and image artifacts. In Fig. LABEL:fig_demo, the PDF+MEDI suffers from image blurring and streaking artifacts, especially large susceptibility errors and streaking artifacts near the brain boundary and brain hemorrhage regions, which is partly due to inaccurate background removal close to brain boundary in PDF.

The presented TFINet approach introduces several important innovations. First, it performs whole brain high-resolution QSM inversion using a neural network, which can avoid patch merging and tiling artifacts. Second, it utilizes a non-local block to increase the receptive fields and capture long-range information for non-local susceptibility estimation. Third, it uses gated convolutions to learn spatial information to spatially adaptive perform inner brain and close to brain boundary susceptibility estimation.

This feasibility study has also demonstrated the ability to use existing standard of SWI raw data to reconstruct QSM for clinical utility. This offers the possibility of QSM use in clinical operation without any additional scans beyond current standard of care protocols. Combining SWI magnitude and QSM estimation images may offer new diagnostic capabilities to assist radiological interpretation. In particular, it is well-known that SWI suffers from blooming artifacts and difficulties in differentiating calcifications and hemosiderin. QSM can overcome these limitations of SWI, which can expand the roles of SWI and QSM in neuroradiology clinical and research arenas. In Fig.7, the calcification is easily differentiated in QSM maps. From Fig.7 and Fig.8, TFINet results show no shading artifacts or streaking artifacts around the lesions, while also preserving the details of fine structures.

V. Conclusion

In summary, a deep-learning-based single-step QSM approach have been demonstrated. It can substantially improve brain susceptibility estimation using clinical QSM data. This capability opens up a wide array of QSM investigations using clinically acquired SWI data to derive QSM maps across a host of neuroimaging indications.

Vi. Tables

SS-TV-QSM LN-QSM PDF+MEDI RESHARP+TKD TFINet
RMSE ()
HFEN ()
SSIM (0-1)
Table 1: Numerical measures of QSM reconstruction quality on 100 synthetic data.

Vii. Figures

Figure 1: Illustration of training datasets.
Figure 2:

Network structure of SSQSMnet. A 3D encoder-decoder architecture was designed with 6 gated convolutional layers (kernel size 3x3x3, dilated rate 1x1x1), 3 gated convolutional layers (kernel size 3x3x3, dilated rate 2x2x2), 4 max pooling layer (pool size 2x2x2, stride size 2x2x2), 1 non-local block, 4 deconvolution layers (kernel size 3x3x3, stride size 2x2x2), 9 normalization layers, 5 feature concatenations, and 1 convolutional layer (kernel size 3x3x3, linear activation).

Figure 3: Total field (a), ground truth QSM images (b) and estimated QSM images (c-g) by five methods from one of synthetic test data in axial (i), coronal (iii), and sagittal (iv) views.
Figure 4: The residual error maps of Fig.3. TFINet results have the least errors compared with all other results.
Figure 5: STI (3,3), TKD, SS-TV-QSM, LN-QSM, and TFINet susceptibility maps from 2016 QSM reconstruction challenge in axial (i), zoom-in axial (ii), coronal (iii), and sagittal (iv) views.
Figure 6: Total fields (a), QSM images (b-f), SWI images (g), and FLAIR images (h) from a 65-year-old patient with stroke. In zoom-in axial (ii), a few hypointense regions of SWI image (black arrows) is hyperintense in QSM, indicating it iron deposition or hemorrhage. In saggital plane (iv), streaking and shading artifacts are clearly visible in SS-TV-QSM, LN-QSM, PDF+MEDI, and RESHARP+TKD images. TFINet can greatly suppress the streaking artifacts.
Figure 7: Total fields (a), QSM images (b-f), SWI images (g), and FLAIR images (h) from a 56-year-old subject with hemorrhagic intracranial metastases. In zoom-in axial (ii), image blurring and shading artifacts is clearly visible in SS-TV-QSM, PDF+MEDI, and RESHARP+TKD images. In TFINet, two small calcification (white dash arrows) is hypointense on SWI image and diamagnetic on QSM image. In coronal and saggital plane (iii, iv), streaking and shading artifacts are clearly visible in SS-TV-QSM, LN-QSM, PDF+MEDI, and RESHARP+TKD images. TFINet can greatly suppress the streaking artifacts.
Figure 8: Total fields (a), QSM images (b-f), SWI images (g), and FLAIR images (h) from a 37-year-old patient with surgical planning. In zoom-in axial (ii), image blurring and shading artifacts is clearly visible in SS-TV-QSM, PDF+MEDI, and RESHARP+TKD images. In coronal and saggital plane (iii, iv), streaking and shading artifacts are clearly visible in SS-TV-QSM, LN-QSM, PDF+MEDI, and RESHARP+TKD images. TFINet results show well-preserved microstrctures and non-visible artifacts.
Figure 9: Total fields (a), QSM images (b-f), SWI images (g), and FLAIR images (h) from from a 54-year-old subject with poststereotactic radiosurgery (SRS) brain metastasis. In zoom-in axial (ii), image blurring and shading artifacts is clearly visible in SS-TV-QSM, PDF+MEDI, and RESHARP+TKD images. In coronal and saggital plane (iii, iv), streaking and shading artifacts are clearly visible in SS-TV-QSM, LN-QSM, PDF+MEDI, and RESHARP+TKD images. TFINet results show well-preserved microstrctures and non-visible artifacts.

References

  • Schenck 1996 John F Schenck. The role of magnetic susceptibility in magnetic resonance imaging: Mri magnetic compatibility of the first and second kinds. Med. Phys., 23(6):815–850, 1996.
  • Haacke et al. 2004 E Mark Haacke, Yingbiao Xu, Yu-Chung N Cheng, and Jürgen R Reichenbach. Susceptibility weighted imaging (swi). Magn. Res. Med., 52(3):612–618, 2004.
  • Wang and Liu 2015 Yi Wang and Tian Liu. Quantitative susceptibility mapping (qsm): decoding mri data for a tissue magnetic biomarker. Magn. Res. Med., 73(1):82–101, 2015.
  • Babikian et al. 2005 Talin Babikian, M Catherin Freier, Karen A Tong, Joshua P Nickerson, Christopher J Wall, Barbara A Holshouser, Todd Burley, Matt L Riggs, and Stephen Ashwal. Susceptibility weighted imaging: neuropsychologic outcome and pediatric head injury. Pediatr. Neurol., 33(3):184–194, 2005.
  • Huang et al. 2015 Yen-Lin Huang, Ying-Sheng Kuo, Ying-Chi Tseng, David Yen-Ting Chen, Wen-Ta Chiu, and Chi-Jen Chen. Susceptibility-weighted mri in mild traumatic brain injury. Neurology, 84(6):580–585, 2015.
  • Thomas et al. 2008 Bejoy Thomas, Sivaraman Somasundaram, Krishnamoorthy Thamburaj, Chandrasekharan Kesavadas, Arun Kumar Gupta, Narendra K Bodhey, and Tirur Raman Kapilamoorthy. Clinical applications of susceptibility weighted mr imaging of the brain–a pictorial review. Neuroradiology, 50(2):105–116, 2008.
  • De Souza et al. 2008 JM De Souza, RC Domingues, LCH Cruz, FS Domingues, T Iasbeck, and EL Gasparetto. Susceptibility-weighted imaging for the evaluation of patients with familial cerebral cavernous malformations: a comparison with t2-weighted fast spin-echo and gradient-echo sequences. Am. J. Neuroradiol., 29(1):154–158, 2008.
  • Sehgal et al. 2006 Vivek Sehgal, Zachary Delproposto, Djamel Haddar, E Mark Haacke, Andrew E Sloan, Lucia J Zamorano, Geoffery Barger, Jiani Hu, Yingbiao Xu, Karthik Praveen Prabhakaran, et al. Susceptibility-weighted imaging to visualize blood products and improve tumor contrast in the study of brain masses. J. Mag. Reson. Im., 24(1):41–51, 2006.
  • Rauscher et al. 2005 Alexander Rauscher, Jan Sedlacik, Markus Barth, Hans-Joachim Mentzel, and Jürgen R Reichenbach. Magnetic susceptibility-weighted mr phase imaging of the human brain. Am. J. Neuroradiol., 26(4):736–742, 2005.
  • Haacke et al. 2009 E Mark Haacke, Malek Makki, Yulin Ge, Megha Maheshwari, Vivek Sehgal, Jiani Hu, Madeswaran Selvan, Zhen Wu, Zahid Latif, Yang Xuan, et al. Characterizing iron deposition in multiple sclerosis lesions using susceptibility weighted imaging. J. Mag. Reson. Im., 29(3):537–544, 2009.
  • Sehgal et al. 2005 Vivek Sehgal, Zachary Delproposto, E Mark Haacke, Karen A Tong, Nathaniel Wycliffe, Daniel K Kido, Yingbiao Xu, Jaladhar Neelavalli, Djamel Haddar, and Jürgen R Reichenbach. Clinical applications of neuroimaging with susceptibility-weighted imaging. J. Mag. Reson. Im., 22(4):439–450, 2005.
  • Haacke et al. 2005 E Mark Haacke, Norman YC Cheng, Michael J House, Qiang Liu, Jaladhar Neelavalli, Robert J Ogg, Asadullah Khan, Muhammad Ayaz, Wolff Kirsch, and Andre Obenaus. Imaging iron stores in the brain using magnetic resonance imaging. Magn. Reson. Imaging, 23(1):1–25, 2005.
  • Deistung et al. 2013a Andreas Deistung, Andreas Schäfer, Ferdinand Schweser, Uta Biedermann, Robert Turner, and Jürgen R Reichenbach. Toward in vivo histology: a comparison of quantitative susceptibility mapping (qsm) with magnitude-, phase-, and r2⁎-imaging at ultra-high magnetic field strength. NeuroImage, 65:299–314, 2013a.
  • Zhang et al. 2015 Jingwei Zhang, Tian Liu, Ajay Gupta, Pascal Spincemaille, Thanh D Nguyen, and Yi Wang. Quantitative mapping of cerebral metabolic rate of oxygen (cmro2) using quantitative susceptibility mapping (qsm). Magn. Res. Med., 74(4):945–952, 2015.
  • Zheng et al. 2013 Weili Zheng, Helen Nichol, Saifeng Liu, Yu-Chung N Cheng, and E Mark Haacke. Measuring iron in the brain using quantitative susceptibility mapping and x-ray fluorescence imaging. NeuroImage, 78:68–74, 2013.
  • Deistung et al. 2013b Andreas Deistung, Ferdinand Schweser, Benedikt Wiestler, Mario Abello, Matthias Roethke, Felix Sahm, Wolfgang Wick, Armin Michael Nagel, Sabine Heiland, Heinz-Peter Schlemmer, et al. Quantitative susceptibility mapping differentiates between blood depositions and calcifications in patients with glioblastoma. PloS One, 8(3):e57924, 2013b.
  • Wisnieff et al. 2015 Cynthia Wisnieff, Sriram Ramanan, John Olesik, Susan Gauthier, Yi Wang, and David Pitt. Quantitative susceptibility mapping (qsm) of white matter multiple sclerosis lesions: interpreting positive susceptibility and the presence of iron. Magn. Res. Med., 74(2):564–570, 2015.
  • Bilgic et al. 2012 Berkin Bilgic, Adolf Pfefferbaum, Torsten Rohlfing, Edith V Sullivan, and Elfar Adalsteinsson. Mri estimates of brain iron concentration in normal aging using quantitative susceptibility mapping. NeuroImage, 59(3):2625–2635, 2012.
  • Lotfipour et al. 2012 Ashley K Lotfipour, Samuel Wharton, Stefan T Schwarz, V Gontu, Andreas Schäfer, Andrew M Peters, Richard W Bowtell, Dorothee P Auer, Penny A Gowland, and Nin PS Bajaj. High resolution magnetic susceptibility mapping of the substantia nigra in parkinson’s disease. J. Magn. Reson Im., 35(1):48–55, 2012.
  • Sharma et al. 2015 Samir D Sharma, Diego Hernando, Debra E Horng, and Scott B Reeder. Quantitative susceptibility mapping in the abdomen as an imaging biomarker of hepatic iron overload. Magn. Res. Med., 74(3):673–683, 2015.
  • Fan et al. 2014 Audrey P Fan, Berkin Bilgic, Louis Gagnon, Thomas Witzel, Himanshu Bhat, Bruce R Rosen, and Elfar Adalsteinsson. Quantitative oxygenation venography from mri phase. Magn. Res. Med., 72(1):149–159, 2014.
  • Koch et al. 2018 Kevin M Koch, Timothy B Meier, Robin Karr, Andrew S Nencka, L Tugan Muftuler, and Michael McCrea. Quantitative susceptibility mapping after sports-related concussion. Am. J. Neuroradiol., 39(7):1215–1221, 2018.
  • Schweser et al. 2011 Ferdinand Schweser, Andreas Deistung, Berengar Wendel Lehr, and Jürgen Rainer Reichenbach. Quantitative imaging of intrinsic magnetic tissue properties using mri signal phase: an approach to in vivo brain iron metabolism? NeuroImage, 54(4):2789–2807, 2011.
  • Sun and Wilman 2014 Hongfu Sun and Alan H Wilman. Background field removal using spherical mean value filtering and tikhonov regularization. Magn. Res. Med., 71(3):1151–1157, 2014.
  • Liu et al. 2011 Tian Liu, Ildar Khalidov, Ludovic de Rochefort, Pascal Spincemaille, Jing Liu, A John Tsiouris, and Yi Wang. A novel background field removal method for mri using projection onto dipole fields. NMR Biomed., 24(9):1129–1136, 2011.
  • Zhou et al. 2014 Dong Zhou, Tian Liu, Pascal Spincemaille, and Yi Wang. Background field removal by solving the laplacian boundary value problem. NMR Biomed., 27(3):312–319, 2014.
  • Berkin et al. 2014 Bilgic Berkin, Langkammer Christian, L Wald Lawrence, and Setsompop Kawin. Single-step qsm with fast reconstruction. ISMRM, 3321, 2014.
  • Kristian et al. 2014 Bredies Kristian, Ropele Stefan, A Poser Benedikt, Barth Markus, and Langkammer Christian. Single-step quantitative susceptibility mapping using total generalized variation and 3d epi. ISMRM, 0604, 2014.
  • Liu et al. 2014 Tian Liu, Dong Zhou, Pascal Spincemaille, and Yi Wang. Differential approach to quantitative susceptibility mapping without background field removal. ISMRM, 0597, 2014.
  • Sharma et al. 2014 Samir D Sharma, Diego Hernando, Debra E Horng, and Scott B Reeder. A joint background field removal and dipole deconvolution approach for quantitative susceptibility mapping in the liver. ISMRM, 0606, 2014.
  • Liu et al. 2017 Zhe Liu, Youngwook Kee, Dong Zhou, Yi Wang, and Pascal Spincemaille. Preconditioned total field inversion (tfi) method for quantitative susceptibility mapping. Magn. Res. Med., 78(1):303–315, 2017.
  • Chatnuntawech et al. 2017 Itthi Chatnuntawech, Patrick McDaniel, Stephen F Cauley, Borjan A Gagoski, Christian Langkammer, Adrian Martin, P Ellen Grant, Lawrence L Wald, Kawin Setsompop, Elfar Adalsteinsson, et al. Single-step quantitative susceptibility mapping with variational penalties. NMR Biomed., 30(4):e3570, 2017.
  • Sun et al. 2018 Hongfu Sun, Yuhan Ma, M Ethan MacDonald, and G Bruce Pike. Whole head quantitative susceptibility mapping using a least-norm direct dipole inversion method. NeuroImage, 2018.
  • Marques and Bowtell 2005 JP Marques and R Bowtell. Application of a fourier-based method for rapid calculation of field inhomogeneity due to spatial variation of magnetic susceptibility. Concepts Magn. Reson. B., 25(1):65–78, 2005.
  • Salomir et al. 2003 Rares Salomir, Baudouin Denis de Senneville, and Chrit TW Moonen. A fast calculation method for magnetic field inhomogeneity due to an arbitrary distribution of bulk susceptibility. Concepts Magn. Reson. B., 19(1):26–34, 2003.
  • Liu et al. 2009 Tian Liu, Pascal Spincemaille, Ludovic De Rochefort, Bryan Kressler, and Yi Wang. Calculation of susceptibility through multiple orientation sampling (cosmos): a method for conditioning the inverse problem from measured magnetic field map to susceptibility source image in mri. Magn. Res. Med., 61(1):196–204, 2009.
  • Liu 2010 Chunlei Liu. Susceptibility tensor imaging. Magn. Res. Med., 63(6):1471–1477, 2010.
  • Ashburner and Friston 2005 John Ashburner and Karl J Friston. Unified segmentation. NeuroImage, 26(3):839–851, 2005.
  • Liu et al. 2019 Juan Liu, Andrew S Nencka, L Tugan Muftuler, Brad Swearingen, Robin Karr, and Kevin M Koch. Quantitative susceptibility inversion through parcellated multiresolution neural networks and k-space substitution. arXiv, 2019.
  • Wang et al. 2018 Xiaolong Wang, Ross Girshick, Abhinav Gupta, and Kaiming He. Non-local neural networks. In IEEE/CVF, pages 7794–7803, 2018.
  • Liu et al. 2012 Jing Liu, Tian Liu, Ludovic de Rochefort, James Ledoux, Ildar Khalidov, Weiwei Chen, A John Tsiouris, Cynthia Wisnieff, Pascal Spincemaille, Martin R Prince, et al. Morphology enabled dipole inversion for quantitative susceptibility mapping using structural consistency between the magnitude image and the susceptibility map. NeuroImage, 59(3):2560–2568, 2012.
  • Bilgic et al. 2015 Berkin Bilgic, Borjan A Gagoski, Stephen F Cauley, Audrey P Fan, Jonathan R Polimeni, P Ellen Grant, Lawrence L Wald, and Kawin Setsompop. Wave-caipi for highly accelerated 3d imaging. Magn. Res. Med., 73(6):2152–2162, 2015.
  • Langkammer et al. 2018 Christian Langkammer, Ferdinand Schweser, Karin Shmueli, Christian Kames, Xu Li, Li Guo, Carlos Milovic, Jinsuh Kim, Hongjiang Wei, Kristian Bredies, et al. Quantitative susceptibility mapping: report from the 2016 reconstruction challenge. Magn. Res. Med., 79(3):1661–1673, 2018.
  • Yi 2019 Wang Yi. Cornell MRI Research Lab, 2019. URL http://weill.cornell.edu/mri/pages/qsm.html.
  • Berkin 2019 Bilgic Berkin. Berkin Bilgic Personal Website, 2019. URL http://martinos.org/~berkin/software.html.
  • Hongfu 2019 Sun Hongfu. LN-QSM Github, 2019. URL https://github.com/sunhongfu/QSM.