1
|
Milovic C, Tejos C, Silva J, Shmueli K, Irarrazaval P. XSIM: A structural similarity index measure optimized for MRI QSM. Magn Reson Med 2024. [PMID: 39176438 DOI: 10.1002/mrm.30271] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2024] [Revised: 07/03/2024] [Accepted: 08/06/2024] [Indexed: 08/24/2024]
Abstract
PURPOSE The structural similarity index measure (SSIM) has become a popular quality metric to evaluate QSM in a way that is closer to human perception than RMS error (RMSE). However, SSIM may overpenalize errors in diamagnetic tissues and underpenalize them in paramagnetic tissues, resulting in biasing. In addition, extreme artifacts may compress the dynamic range, resulting in unrealistically high SSIM scores (hacking). To overcome biasing and hacking, we propose XSIM: SSIM implemented in the native QSM range, and with internal parameters optimized for QSM. METHODS We used forward simulations from a COSMOS ground-truth brain susceptibility map included in the 2016 QSM Reconstruction Challenge to investigate the effect of QSM reconstruction errors on the SSIM, XSIM, and RMSE metrics. We also used these metrics to optimize QSM reconstructions of the in vivo challenge data set. We repeated this experiment with the QSM abdominal phantom. To validate the use of XSIM instead of SSIM for QSM quality assessment across a range of different reconstruction techniques/algorithms, we analyzed the reconstructions submitted to the 2019 QSM Reconstruction Challenge 2.0. RESULTS Our experiments confirmed the biasing and hacking effects on the SSIM metric applied to QSM. The XSIM metric was robust to those effects, penalizing the presence of streaking artifacts and reconstruction errors. Using XSIM to optimize QSM reconstruction regularization weights returned less overregularization than SSIM and RMSE. CONCLUSION XSIM is recommended over traditional SSIM to evaluate QSM reconstructions against a known ground truth, as it avoids biasing and hacking effects and provides a larger dynamic range of scores.
Collapse
Affiliation(s)
- Carlos Milovic
- School of Electrical Engineering, Pontificia Universidad Catolica de Valparaiso, Valparaiso, Chile
| | - Cristian Tejos
- Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Institute for Intelligent Healthcare Engineering (iHEALTH), Santiago, Chile
| | - Javier Silva
- Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
| | - Karin Shmueli
- Department of Medical Physics and Biomedical Engineering, University College London, London, UK
| | - Pablo Irarrazaval
- Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Institute for Intelligent Healthcare Engineering (iHEALTH), Santiago, Chile
- Institute for Biological and Medical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
| |
Collapse
|
2
|
Venkatesh V, Mathew RS, Yalavarthy PK. Spinet-QSM: model-based deep learning with schatten p-norm regularization for improved quantitative susceptibility mapping. MAGMA (NEW YORK, N.Y.) 2024; 37:411-427. [PMID: 38598165 DOI: 10.1007/s10334-024-01158-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/30/2023] [Revised: 03/16/2024] [Accepted: 03/19/2024] [Indexed: 04/11/2024]
Abstract
OBJECTIVE Quantitative susceptibility mapping (QSM) provides an estimate of the magnetic susceptibility of tissue using magnetic resonance (MR) phase measurements. The tissue magnetic susceptibility (source) from the measured magnetic field distribution/local tissue field (effect) inherent in the MR phase images is estimated by numerically solving the inverse source-effect problem. This study aims to develop an effective model-based deep-learning framework to solve the inverse problem of QSM. MATERIALS AND METHODS This work proposes a Schatten p -norm-driven model-based deep learning framework for QSM with a learnable norm parameter p to adapt to the data. In contrast to other model-based architectures that enforce the l2 -norm or l1 -norm for the denoiser, the proposed approach can enforce any p -norm ( 0 < p ≤ 2 ) on a trainable regulariser. RESULTS The proposed method was compared with deep learning-based approaches, such as QSMnet, and model-based deep learning approaches, such as learned proximal convolutional neural network (LPCNN). Reconstructions performed using 77 imaging volumes with different acquisition protocols and clinical conditions, such as hemorrhage and multiple sclerosis, showed that the proposed approach outperformed existing state-of-the-art methods by a significant margin in terms of quantitative merits. CONCLUSION The proposed SpiNet-QSM showed a consistent improvement of at least 5% in terms of the high-frequency error norm (HFEN) and normalized root mean squared error (NRMSE) over other QSM reconstruction methods with limited training data.
Collapse
Affiliation(s)
- Vaddadi Venkatesh
- Department of Computational and Data Sciences, Indian Institute of Science, Bangalore, Karnataka, 560012, India
| | - Raji Susan Mathew
- School of Data Science, Indian Institute of Science Education and Research, Thiruvananthapuram, Kerala, 695551, India
| | - Phaneendra K Yalavarthy
- Department of Computational and Data Sciences, Indian Institute of Science, Bangalore, Karnataka, 560012, India.
| |
Collapse
|
3
|
Gkotsoulias DG, Jäger C, Müller R, Gräßle T, Olofsson KM, Møller T, Unwin S, Crockford C, Wittig RM, Bilgic B, Möller HE. Chaos and COSMOS-Considerations on QSM methods with multiple and single orientations and effects from local anisotropy. Magn Reson Imaging 2024; 110:104-111. [PMID: 38631534 DOI: 10.1016/j.mri.2024.04.020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2024] [Revised: 04/07/2024] [Accepted: 04/14/2024] [Indexed: 04/19/2024]
Abstract
PURPOSE Field-to-susceptibility inversion in quantitative susceptibility mapping (QSM) is ill-posed and needs numerical stabilization through either regularization or oversampling by acquiring data at three or more object orientations. Calculation Of Susceptibility through Multiple Orientations Sampling (COSMOS) is an established oversampling approach and regarded as QSM gold standard. It achieves a well-conditioned inverse problem, requiring rotations by 0°, 60° and 120° in the yz-plane. However, this is impractical in vivo, where head rotations are typically restricted to a range of ±25°. Non-ideal sampling degrades the conditioning with residual streaking artifacts whose mitigation needs further regularization. Moreover, susceptibility anisotropy in white matter is not considered in the COSMOS model, which may introduce additional bias. The current work presents a thorough investigation of these effects in primate brain. METHODS Gradient-recalled echo (GRE) data of an entire fixed chimpanzee brain were acquired at 7 T (350 μm resolution, 10 orientations) including ideal COSMOS sampling and realistic rotations in vivo. Comparisons of the results included ideal COSMOS, in-vivo feasible acquisitions with 3-8 orientations and single-orientation iLSQR QSM. RESULTS In-vivo feasible and optimal COSMOS yielded high-quality susceptibility maps with increased SNR resulting from averaging multiple acquisitions. COSMOS reconstructions from non-ideal rotations about a single axis required additional L2-regularization to mitigate residual streaking artifacts. CONCLUSION In view of unconsidered anisotropy effects, added complexity of the reconstruction, and the general challenge of multi-orientation acquisitions, advantages of sub-optimal COSMOS schemes over regularized single-orientation QSM appear limited in in-vivo settings.
Collapse
Affiliation(s)
- Dimitrios G Gkotsoulias
- Nuclear Magnetic Resonance Methods & Development Group, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany.
| | - Carsten Jäger
- Department of Neurophysics, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
| | - Roland Müller
- Nuclear Magnetic Resonance Methods & Development Group, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
| | - Tobias Gräßle
- Epidemiology of Highly Pathogenic Microorganisms, Robert Koch-Institute, Berlin, Germany
| | | | | | - Steve Unwin
- Wildlife Health Australia, Canberra, Australia
| | - Catherine Crockford
- Department of Human Behavior, Ecology and Culture, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany; The Ape Social Mind Lab, Institut des Sciences Cognitives Marc Jeannerod, Bron, France; Taï Chimpanzee Project, Centre Suisse de Recherches Scientifiques, Abidjan, Côte d'Ivoire
| | - Roman M Wittig
- Department of Human Behavior, Ecology and Culture, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany; The Ape Social Mind Lab, Institut des Sciences Cognitives Marc Jeannerod, Bron, France; Taï Chimpanzee Project, Centre Suisse de Recherches Scientifiques, Abidjan, Côte d'Ivoire
| | - Berkin Bilgic
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Boston, MA, USA; Harvard Medical School, Boston, MA, USA; Harvard/MIT Health Sciences and Technology, Cambridge, MA, United States
| | - Harald E Möller
- Nuclear Magnetic Resonance Methods & Development Group, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
| |
Collapse
|
4
|
Bilgic B, Costagli M, Chan KS, Duyn J, Langkammer C, Lee J, Li X, Liu C, Marques JP, Milovic C, Robinson SD, Schweser F, Shmueli K, Spincemaille P, Straub S, van Zijl P, Wang Y. Recommended implementation of quantitative susceptibility mapping for clinical research in the brain: A consensus of the ISMRM electro-magnetic tissue properties study group. Magn Reson Med 2024; 91:1834-1862. [PMID: 38247051 PMCID: PMC10950544 DOI: 10.1002/mrm.30006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2023] [Revised: 10/31/2023] [Accepted: 12/14/2023] [Indexed: 01/23/2024]
Abstract
This article provides recommendations for implementing QSM for clinical brain research. It is a consensus of the International Society of Magnetic Resonance in Medicine, Electro-Magnetic Tissue Properties Study Group. While QSM technical development continues to advance rapidly, the current QSM methods have been demonstrated to be repeatable and reproducible for generating quantitative tissue magnetic susceptibility maps in the brain. However, the many QSM approaches available have generated a need in the neuroimaging community for guidelines on implementation. This article outlines considerations and implementation recommendations for QSM data acquisition, processing, analysis, and publication. We recommend that data be acquired using a monopolar 3D multi-echo gradient echo (GRE) sequence and that phase images be saved and exported in Digital Imaging and Communications in Medicine (DICOM) format and unwrapped using an exact unwrapping approach. Multi-echo images should be combined before background field removal, and a brain mask created using a brain extraction tool with the incorporation of phase-quality-based masking. Background fields within the brain mask should be removed using a technique based on SHARP or PDF, and the optimization approach to dipole inversion should be employed with a sparsity-based regularization. Susceptibility values should be measured relative to a specified reference, including the common reference region of the whole brain as a region of interest in the analysis. The minimum acquisition and processing details required when reporting QSM results are also provided. These recommendations should facilitate clinical QSM research and promote harmonized data acquisition, analysis, and reporting.
Collapse
Affiliation(s)
- Berkin Bilgic
- Martinos Center for Biomedical Imaging, Massachusetts General Hospital and Harvard Medical School, Charlestown, Massachusetts, USA
| | - Mauro Costagli
- Department of Neuroscience, Rehabilitation, Ophthalmology, Genetics, Maternal and Child Sciences (DINOGMI), University of Genoa, Genoa, Italy
- Laboratory of Medical Physics and Magnetic Resonance, IRCCS Stella Maris, Pisa, Italy
| | - Kwok-Shing Chan
- Martinos Center for Biomedical Imaging, Massachusetts General Hospital and Harvard Medical School, Charlestown, Massachusetts, USA
- Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands
| | - Jeff Duyn
- Advanced MRI Section, NINDS, National Institutes of Health, Bethesda, Maryland, USA
| | | | - Jongho Lee
- Department of Electrical and Computer Engineering, Seoul National University, Seoul, Republic of Korea
| | - Xu Li
- Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, USA
| | - Chunlei Liu
- Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, California, USA
- Helen Wills Neuroscience Institute, University of California, Berkeley, California, USA
| | - José P Marques
- Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands
| | - Carlos Milovic
- School of Electrical Engineering (EIE), Pontificia Universidad Catolica de Valparaiso, Valparaiso, Chile
| | - Simon Daniel Robinson
- High Field MR Centre, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
- Centre of Advanced Imaging, University of Queensland, Brisbane, Australia
| | - Ferdinand Schweser
- Buffalo Neuroimaging Analysis Center, Department of Neurology, Jacobs School of Medicine and Biomedical Sciences at the University at Buffalo, Buffalo, New York, USA
- Center for Biomedical Imaging, Clinical and Translational Science Institute at the University at Buffalo, Buffalo, New York, USA
| | - Karin Shmueli
- Department of Medical Physics and Biomedical Engineering, University College London, London, UK
| | - Pascal Spincemaille
- MRI Research Institute, Department of Radiology, Weill Cornell Medicine, New York, New York, USA
| | - Sina Straub
- Department of Radiology, Mayo Clinic, Jacksonville, Florida, USA
| | - Peter van Zijl
- Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, USA
| | - Yi Wang
- MRI Research Institute, Departments of Radiology and Biomedical Engineering, Cornell University, New York, New York, USA
| |
Collapse
|
5
|
Chen L, Shin HG, van Zijl PC, Li X. Exploiting gradient-echo frequency evolution: Probing white matter microstructure and extracting bulk susceptibility-induced frequency for quantitative susceptibility mapping. Magn Reson Med 2024; 91:1676-1693. [PMID: 38102838 PMCID: PMC10880384 DOI: 10.1002/mrm.29958] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2023] [Revised: 10/08/2023] [Accepted: 11/17/2023] [Indexed: 12/17/2023]
Abstract
PURPOSE This work is to investigate the microstructure-induced frequency shift in white matter (WM) with crossing fibers and to separate the microstructure-related frequency shift from the bulk susceptibility-induced frequency shift by model fitting the gradient-echo (GRE) frequency evolution for potentially more accurate quantitative susceptibility mapping (QSM). METHODS A hollow-cylinder fiber model (HCFM) with two fiber populations was developed to investigate GRE frequency evolutions in WM voxels with microstructural orientation dispersion. The simulated and experimentally measured TE-dependent local frequency shift was then fitted to a simplified frequency evolution model to obtain a microstructure-related frequency difference parameter (∆ f $$ \Delta f $$ ) and a TE-independent bulk susceptibility-induced frequency shift (C f $$ {C}_f $$ ). The obtainedC f $$ {C}_f $$ was then used for QSM reconstruction. Reconstruction performances were evaluated using a numerical head phantom and in vivo data and then compared to other multi-echo combination methods. RESULTS GRE frequency evolutions and∆ f $$ \Delta f $$ -based tissue parameters in both parallel and crossing fibers determined from our simulations were comparable to those observed in vivo. The TE-dependent frequency fitting method outperformed other multi-echo combination methods in estimatingC f $$ {C}_f $$ in simulations. The fitted∆ f $$ \Delta f $$ ,C f $$ {C}_f $$ , and QSM could be improved further by navigator-based B0 fluctuation correction. CONCLUSION A HCFM with two fiber populations can be used to characterize microstructure-induced frequency shifts in WM regions with crossing fibers. HCFM-based TE-dependent frequency fitting provides tissue contrast related to microstructure (∆ f $$ \Delta f $$ ) and in addition may help improve the quantification accuracy ofC f $$ {C}_f $$ and the corresponding QSM.
Collapse
Affiliation(s)
- Lin Chen
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, United States
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, Maryland, United States
| | - Hyeong-Geol Shin
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, United States
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, Maryland, United States
| | - Peter C.M. van Zijl
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, United States
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, Maryland, United States
| | - Xu Li
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, United States
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, Maryland, United States
| |
Collapse
|
6
|
García Saborit M, Jara A, Muñoz N, Milovic C, Tepper A, Alliende LM, Mena C, Iruretagoyena B, Ramirez-Mahaluf JP, Diaz C, Nachar R, Castañeda CP, González A, Undurraga J, Crossley N, Tejos C. Quantitative Susceptibility Mapping MRI in Deep-Brain Nuclei in First-Episode Psychosis. Schizophr Bull 2023; 49:1355-1363. [PMID: 37030007 PMCID: PMC10483330 DOI: 10.1093/schbul/sbad041] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 04/10/2023]
Abstract
BACKGROUND Psychosis is related to neurochemical changes in deep-brain nuclei, particularly suggesting dopamine dysfunctions. We used an magnetic resonance imaging-based technique called quantitative susceptibility mapping (QSM) to study these regions in psychosis. QSM quantifies magnetic susceptibility in the brain, which is associated with iron concentrations. Since iron is a cofactor in dopamine pathways and co-localizes with inhibitory neurons, differences in QSM could reflect changes in these processes. METHODS We scanned 83 patients with first-episode psychosis and 64 healthy subjects. We reassessed 22 patients and 21 control subjects after 3 months. Mean susceptibility was measured in 6 deep-brain nuclei. Using linear mixed models, we analyzed the effect of case-control differences, region, age, gender, volume, framewise displacement (FD), treatment duration, dose, laterality, session, and psychotic symptoms on QSM. RESULTS Patients showed a significant susceptibility reduction in the putamen and globus pallidus externa (GPe). Patients also showed a significant R2* reduction in GPe. Age, gender, FD, session, group, and region are significant predictor variables for QSM. Dose, treatment duration, and volume were not predictor variables of QSM. CONCLUSIONS Reduction in QSM and R2* suggests a decreased iron concentration in the GPe of patients. Susceptibility reduction in putamen cannot be associated with iron changes. Since changes observed in putamen and GPe were not associated with symptoms, dose, and treatment duration, we hypothesize that susceptibility may be a trait marker rather than a state marker, but this must be verified with long-term studies.
Collapse
Affiliation(s)
- Marisleydis García Saborit
- Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Institute for Intelligent Healthcare Engineering, Santiago, Chile
| | - Alejandro Jara
- Department of Statistics, Mathematics Faculty, Pontificia Universidad Catolica de Chile, Santiago, Chile
| | - Néstor Muñoz
- Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Institute for Intelligent Healthcare Engineering, Santiago, Chile
| | - Carlos Milovic
- School of Electrical Engineering, Pontificia Universidad Catolica de Valparaiso, Valparaiso, Chile
| | - Angeles Tepper
- Millennium Institute for Intelligent Healthcare Engineering, Santiago, Chile
- Department of Psychiatry, School of Medicine, Pontificia Universidad Catolica de Chile, Santiago, Chile
| | - Luz María Alliende
- Department of Psychiatry, School of Medicine, Pontificia Universidad Catolica de Chile, Santiago, Chile
| | - Carlos Mena
- Department of Psychiatry, School of Medicine, Pontificia Universidad Catolica de Chile, Santiago, Chile
| | - Bárbara Iruretagoyena
- Department of Psychiatry, School of Medicine, Pontificia Universidad Catolica de Chile, Santiago, Chile
| | | | - Camila Diaz
- Pharmacovigilance, Instituto Psiquiátrico Dr. J. Horwitz Barak, Santiago, Chile
| | - Ruben Nachar
- Pharmacovigilance, Instituto Psiquiátrico Dr. J. Horwitz Barak, Santiago, Chile
| | | | - Alfonso González
- Early Intervention Program, Instituto Psiquiátrico Dr J. Horwitz Barak, Santiago, Chile
- School of Medicine, Universidad Finis Terrae, Santiago, Chile
| | - Juan Undurraga
- Early Intervention Program, Instituto Psiquiátrico Dr J. Horwitz Barak, Santiago, Chile
- Department of Neurology and Psychiatry, Faculty of Medicine, Clínica Alemana Universidad del Desarrollo, Santiago, Chile
| | - Nicolas Crossley
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Institute for Intelligent Healthcare Engineering, Santiago, Chile
- Department of Psychiatry, School of Medicine, Pontificia Universidad Catolica de Chile, Santiago, Chile
| | - Cristian Tejos
- Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Institute for Intelligent Healthcare Engineering, Santiago, Chile
| |
Collapse
|
7
|
Gkotsoulias DG, Müller R, Jäger C, Schlumm T, Mildner T, Eichner C, Pampel A, Jaffe J, Gräßle T, Alsleben N, Chen J, Crockford C, Wittig R, Liu C, Möller HE. High angular resolution susceptibility imaging and estimation of fiber orientation distribution functions in primate brain. Neuroimage 2023; 276:120202. [PMID: 37247762 DOI: 10.1016/j.neuroimage.2023.120202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2022] [Revised: 05/21/2023] [Accepted: 05/27/2023] [Indexed: 05/31/2023] Open
Abstract
Uncovering brain-tissue microstructure including axonal characteristics is a major neuroimaging research focus. Within this scope, anisotropic properties of magnetic susceptibility in white matter have been successfully employed to estimate primary axonal trajectories using mono-tensorial models. However, anisotropic susceptibility has not yet been considered for modeling more complex fiber structures within a voxel, such as intersecting bundles, or an estimation of orientation distribution functions (ODFs). This information is routinely obtained by high angular resolution diffusion imaging (HARDI) techniques. In applications to fixed tissue, however, diffusion-weighted imaging suffers from an inherently low signal-to-noise ratio and limited spatial resolution, leading to high demands on the performance of the gradient system in order to mitigate these limitations. In the current work, high angular resolution susceptibility imaging (HARSI) is proposed as a novel, phase-based methodology to estimate ODFs. A multiple gradient-echo dataset was acquired in an entire fixed chimpanzee brain at 61 orientations by reorienting the specimen in the magnetic field. The constant solid angle method was adapted for estimating phase-based ODFs. HARDI data were also acquired for comparison. HARSI yielded information on whole-brain fiber architecture, including identification of peaks of multiple bundles that resembled features of the HARDI results. Distinct differences between both methods suggest that susceptibility properties may offer complementary microstructural information. These proof-of-concept results indicate a potential to study the axonal organization in post-mortem primate and human brain at high resolution.
Collapse
Affiliation(s)
- Dimitrios G Gkotsoulias
- Nuclear Magnetic Resonance Methods & Development Group, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany.
| | - Roland Müller
- Nuclear Magnetic Resonance Methods & Development Group, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
| | - Carsten Jäger
- Department of Neurophysics, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
| | - Torsten Schlumm
- Nuclear Magnetic Resonance Methods & Development Group, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
| | - Toralf Mildner
- Nuclear Magnetic Resonance Methods & Development Group, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
| | - Cornelius Eichner
- Department of Neuropsychology, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
| | - André Pampel
- Nuclear Magnetic Resonance Methods & Development Group, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
| | - Jennifer Jaffe
- Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany; Taï Chimpanzee Project, Centre Suisse de Recherches Scientifiques en Côte d'Ivoire, Côte d'Ivoire
| | - Tobias Gräßle
- Taï Chimpanzee Project, Centre Suisse de Recherches Scientifiques en Côte d'Ivoire, Côte d'Ivoire; Helmholtz Institute for One Health, Greifswald, Germany; Robert Koch Institute, Epidemiology of Highly Pathogenic Microorganisms, Berlin, Germany
| | - Niklas Alsleben
- Department of Neurophysics, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
| | - Jingjia Chen
- Electrical Engineering and Computer Sciences, University of California, Berkeley, CA, USA
| | - Catherine Crockford
- Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany; Taï Chimpanzee Project, Centre Suisse de Recherches Scientifiques en Côte d'Ivoire, Côte d'Ivoire; Institute of Cognitive Sciences, CNRS UMR5229 University of Lyon, Bron, France
| | - Roman Wittig
- Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany; Taï Chimpanzee Project, Centre Suisse de Recherches Scientifiques en Côte d'Ivoire, Côte d'Ivoire; Institute of Cognitive Sciences, CNRS UMR5229 University of Lyon, Bron, France
| | - Chunlei Liu
- Electrical Engineering and Computer Sciences, University of California, Berkeley, CA, USA; Helen Wills Neuroscience Institute, University of California, Berkeley, CA, USA
| | - Harald E Möller
- Nuclear Magnetic Resonance Methods & Development Group, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
| |
Collapse
|
8
|
Bilgic B, Costagli M, Chan KS, Duyn J, Langkammer C, Lee J, Li X, Liu C, Marques JP, Milovic C, Robinson S, Schweser F, Shmueli K, Spincemaille P, Straub S, van Zijl P, Wang Y. Recommended Implementation of Quantitative Susceptibility Mapping for Clinical Research in The Brain: A Consensus of the ISMRM Electro-Magnetic Tissue Properties Study Group. ARXIV 2023:arXiv:2307.02306v1. [PMID: 37461418 PMCID: PMC10350101] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Subscribe] [Scholar Register] [Indexed: 07/22/2023]
Abstract
This article provides recommendations for implementing quantitative susceptibility mapping (QSM) for clinical brain research. It is a consensus of the ISMRM Electro-Magnetic Tissue Properties Study Group. While QSM technical development continues to advance rapidly, the current QSM methods have been demonstrated to be repeatable and reproducible for generating quantitative tissue magnetic susceptibility maps in the brain. However, the many QSM approaches available give rise to the need in the neuroimaging community for guidelines on implementation. This article describes relevant considerations and provides specific implementation recommendations for all steps in QSM data acquisition, processing, analysis, and presentation in scientific publications. We recommend that data be acquired using a monopolar 3D multi-echo GRE sequence, that phase images be saved and exported in DICOM format and unwrapped using an exact unwrapping approach. Multi-echo images should be combined before background removal, and a brain mask created using a brain extraction tool with the incorporation of phase-quality-based masking. Background fields should be removed within the brain mask using a technique based on SHARP or PDF, and the optimization approach to dipole inversion should be employed with a sparsity-based regularization. Susceptibility values should be measured relative to a specified reference, including the common reference region of whole brain as a region of interest in the analysis, and QSM results should be reported with - as a minimum - the acquisition and processing specifications listed in the last section of the article. These recommendations should facilitate clinical QSM research and lead to increased harmonization in data acquisition, analysis, and reporting.
Collapse
Affiliation(s)
- Berkin Bilgic
- Martinos Center for Biomedical Imaging, Massachusetts General Hospital and Harvard Medical School, Charlestown, MA, United States
| | - Mauro Costagli
- Department of Neuroscience, Rehabilitation, Ophthalmology, Genetics, Maternal and Child Sciences (DINOGMI), University of Genoa, Genoa, Italy
- Laboratory of Medical Physics and Magnetic Resonance, IRCCS Stella Maris, Pisa, Italy
| | - Kwok-Shing Chan
- Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, the Netherlands
| | - Jeff Duyn
- Advanced MRI Section, NINDS, National Institutes of Health, Bethesda, MD, United States
| | | | - Jongho Lee
- Department of Electrical and Computer Engineering, Seoul National University, Seoul, Republic of Korea
| | - Xu Li
- Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, Baltimore, MD, United States
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, MD, United States
| | - Chunlei Liu
- Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, CA, USA
- Helen Wills Neuroscience Institute, University of California, Berkeley, CA, USA
| | - José P Marques
- Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, the Netherlands
| | - Carlos Milovic
- School of Electrical Engineering (EIE), Pontificia Universidad Catolica de Valparaiso, Valparaiso, Chile
| | - Simon Robinson
- High Field MR Centre, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Austria
| | - Ferdinand Schweser
- Buffalo Neuroimaging Analysis Center, Department of Neurology, Jacobs School of Medicine and Biomedical Sciences at the University at Buffalo, Buffalo, NY, USA
- Center for Biomedical Imaging, Clinical and Translational Science Institute at the University at Buffalo, Buffalo, NY, United States
| | - Karin Shmueli
- Department of Medical Physics and Biomedical Engineering, University College London, London, UK
| | - Pascal Spincemaille
- MRI Research Institute, Department of Radiology, Weill Cornell Medicine, New York, NY, United States
| | - Sina Straub
- Department of Radiology, Mayo Clinic, Jacksonville, FL, United States
| | - Peter van Zijl
- Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, Baltimore, MD, United States
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, MD, United States
| | - Yi Wang
- MRI Research Institute, Departments of Radiology and Biomedical Engineering, Cornell University, New York, NY, United States
| |
Collapse
|
9
|
Fang Z, Lai KW, van Zijl P, Li X, Sulam J. DeepSTI: Towards tensor reconstruction using fewer orientations in susceptibility tensor imaging. Med Image Anal 2023; 87:102829. [PMID: 37146440 PMCID: PMC10288385 DOI: 10.1016/j.media.2023.102829] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2022] [Revised: 03/11/2023] [Accepted: 04/18/2023] [Indexed: 05/07/2023]
Abstract
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.
Collapse
Affiliation(s)
- Zhenghan Fang
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21218, USA; Johns Hopkins Kavli Neuroscience Discovery Institute, Baltimore, MD 21218, USA
| | - Kuo-Wei Lai
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21218, USA
| | - Peter van Zijl
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, MD 21205, USA; Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, MD 21205, USA
| | - Xu Li
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, MD 21205, USA; Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, MD 21205, USA.
| | - Jeremias Sulam
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21218, USA; Johns Hopkins Kavli Neuroscience Discovery Institute, Baltimore, MD 21218, USA.
| |
Collapse
|
10
|
Silva J, Milovic C, Lambert M, Montalba C, Arrieta C, Irarrazaval P, Uribe S, Tejos C. Toward a realistic in silico abdominal phantom for QSM. Magn Reson Med 2023; 89:2402-2418. [PMID: 36695213 PMCID: PMC10952412 DOI: 10.1002/mrm.29597] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2022] [Revised: 12/18/2022] [Accepted: 01/06/2023] [Indexed: 01/26/2023]
Abstract
PURPOSE QSM outside the brain has recently gained interest, particularly in the abdominal region. However, the absence of reliable ground truths makes difficult to assess reconstruction algorithms, whose quality is already compromised by additional signal contributions from fat, gases, and different kinds of motion. This work presents a realistic in silico phantom for the development, evaluation and comparison of abdominal QSM reconstruction algorithms. METHODS Synthetic susceptibility andR 2 * $$ {R}_2^{\ast } $$ maps were generated by segmenting and postprocessing the abdominal 3T MRI data from a healthy volunteer. Susceptibility andR 2 * $$ {R}_2^{\ast } $$ values in different tissues/organs were assigned according to literature and experimental values and were also provided with realistic textures. The signal was simulated using as input the synthetic QSM andR 2 * $$ {R}_2^{\ast } $$ maps and fat contributions. Three susceptibility scenarios and two acquisition protocols were simulated to compare different reconstruction algorithms. RESULTS QSM reconstructions show that the phantom allows to identify the main strengths and limitations of the acquisition approaches and reconstruction algorithms, such as in-phase acquisitions, water-fat separation methods, and QSM dipole inversion algorithms. CONCLUSION The phantom showed its potential as a ground truth to evaluate and compare reconstruction pipelines and algorithms. The publicly available source code, designed in a modular framework, allows users to easily modify the susceptibility,R 2 * $$ {R}_2^{\ast } $$ and TEs, and thus creates different abdominal scenarios.
Collapse
Affiliation(s)
- Javier Silva
- Department of Electrical EngineeringPontificia Universidad Católica de Chile
SantiagoChile
- Biomedical Imaging CenterPontificia Universidad Católica de ChileSantiagoChile
- Millennium Institute for Intelligent Healthcare Engineering (iHEALTH)SantiagoChile
| | - Carlos Milovic
- Millennium Institute for Intelligent Healthcare Engineering (iHEALTH)SantiagoChile
- School of Electrical EngineeringPontificia Universidad Católica de ValparaísoValparaísoChile
| | - Mathias Lambert
- Department of Electrical EngineeringPontificia Universidad Católica de Chile
SantiagoChile
- Biomedical Imaging CenterPontificia Universidad Católica de ChileSantiagoChile
- Millennium Institute for Intelligent Healthcare Engineering (iHEALTH)SantiagoChile
| | - Cristian Montalba
- Biomedical Imaging CenterPontificia Universidad Católica de ChileSantiagoChile
- Millennium Institute for Intelligent Healthcare Engineering (iHEALTH)SantiagoChile
- Department of Radiology, School of MedicinePontificia Universidad Católica de ChileSantiagoChile
| | - Cristóbal Arrieta
- Biomedical Imaging CenterPontificia Universidad Católica de ChileSantiagoChile
- Millennium Institute for Intelligent Healthcare Engineering (iHEALTH)SantiagoChile
| | - Pablo Irarrazaval
- Department of Electrical EngineeringPontificia Universidad Católica de Chile
SantiagoChile
- Biomedical Imaging CenterPontificia Universidad Católica de ChileSantiagoChile
- Millennium Institute for Intelligent Healthcare Engineering (iHEALTH)SantiagoChile
- Institute for Biological and Medical Engineering, Pontificia Universidad Católica de ChileSantiagoChile
| | - Sergio Uribe
- Biomedical Imaging CenterPontificia Universidad Católica de ChileSantiagoChile
- Millennium Institute for Intelligent Healthcare Engineering (iHEALTH)SantiagoChile
- Department of Radiology, School of MedicinePontificia Universidad Católica de ChileSantiagoChile
| | - Cristian Tejos
- Department of Electrical EngineeringPontificia Universidad Católica de Chile
SantiagoChile
- Biomedical Imaging CenterPontificia Universidad Católica de ChileSantiagoChile
- Millennium Institute for Intelligent Healthcare Engineering (iHEALTH)SantiagoChile
| |
Collapse
|
11
|
Yao J, Morrison MA, Jakary A, Avadiappan S, Chen Y, Luitjens J, Glueck J, Driscoll T, Geschwind MD, Nelson AB, Villanueva-Meyer JE, Hess CP, Lupo JM. Comparison of quantitative susceptibility mapping methods for iron-sensitive susceptibility imaging at 7T: An evaluation in healthy subjects and patients with Huntington's disease. Neuroimage 2023; 265:119788. [PMID: 36476567 DOI: 10.1016/j.neuroimage.2022.119788] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Revised: 10/08/2022] [Accepted: 12/02/2022] [Indexed: 12/12/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Jingwen Yao
- Department of Radiology and Biomedical Imaging, UCSF, San Francisco, CA, USA
| | - Melanie A Morrison
- Department of Radiology and Biomedical Imaging, UCSF, San Francisco, CA, USA
| | - Angela Jakary
- Department of Radiology and Biomedical Imaging, UCSF, San Francisco, CA, USA
| | - Sivakami Avadiappan
- Department of Radiology and Biomedical Imaging, UCSF, San Francisco, CA, USA
| | - Yicheng Chen
- Department of Radiology and Biomedical Imaging, UCSF, San Francisco, CA, USA; UCSF/UC Berkeley Graduate Program in Bioengineering, San Francisco & Berkeley, CA, USA; Meta Platforms, Inc., Mountain View, CA, USA
| | - Johanna Luitjens
- Department of Radiology and Biomedical Imaging, UCSF, San Francisco, CA, USA; Department of Radiology, University Hospital, LMU Munich, Munich, Germany
| | - Julia Glueck
- Department of Neurology, University of California San Francisco, San Francisco, CA, USA
| | - Theresa Driscoll
- Department of Neurology, University of California San Francisco, San Francisco, CA, USA
| | - Michael D Geschwind
- Department of Neurology, University of California San Francisco, San Francisco, CA, USA
| | - Alexandra B Nelson
- Department of Neurology, University of California San Francisco, San Francisco, CA, USA
| | | | - Christopher P Hess
- Department of Radiology and Biomedical Imaging, UCSF, San Francisco, CA, USA; Department of Neurology, University of California San Francisco, San Francisco, CA, USA
| | - Janine M Lupo
- Department of Radiology and Biomedical Imaging, UCSF, San Francisco, CA, USA; UCSF/UC Berkeley Graduate Program in Bioengineering, San Francisco & Berkeley, CA, USA.
| |
Collapse
|
12
|
Kamesh Iyer S, Moon BF, Josselyn N, Kurtz RM, Song JW, Ware JB, Nabavizadeh SA, Witschey WR. Quantitative susceptibility mapping using plug-and-play alternating direction method of multipliers. Sci Rep 2022; 12:21679. [PMID: 36522372 PMCID: PMC9755132 DOI: 10.1038/s41598-022-22778-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2022] [Accepted: 10/19/2022] [Indexed: 12/23/2022] Open
Abstract
Quantitative susceptibility mapping employs regularization to reduce artifacts, yet many recent denoisers are unavailable for reconstruction. We developed a plug-and-play approach to QSM reconstruction (PnP QSM) and show its flexibility using several patch-based denoisers. We developed PnP QSM using alternating direction method of multiplier framework and applied collaborative filtering denoisers. We apply the technique to the 2016 QSM Challenge and in 10 glioblastoma multiforme datasets. We compared its performance with four published QSM techniques and a multi-orientation QSM method. We analyzed magnetic susceptibility accuracy using brain region-of-interest measurements, and image quality using global error metrics. Reconstructions on glioblastoma data were analyzed using ranked and semiquantitative image grading by three neuroradiologist observers to assess image quality (IQ) and sharpness (IS). PnP-BM4D QSM showed good correlation (β = 0.84, R2 = 0.98, p < 0.05) with COSMOS and no significant bias (bias = 0.007 ± 0.012). PnP-BM4D QSM achieved excellent quality when assessed using structural similarity index metric (SSIM = 0.860), high frequency error norm (HFEN = 58.5), cross correlation (CC = 0.804), and mutual information (MI = 0.475) and also maintained good conspicuity of fine features. In glioblastoma datasets, PnP-BM4D QSM showed higher performance (IQGrade = 2.4 ± 0.4, ISGrade = 2.7 ± 0.3, IQRank = 3.7 ± 0.3, ISRank = 3.9 ± 0.3) compared to MEDI (IQGrade = 2.1 ± 0.5, ISGrade = 2.1 ± 0.6, IQRank = 2.4 ± 0.6, ISRank = 2.9 ± 0.2) and FANSI-TGV (IQGrade = 2.2 ± 0.6, ISGrade = 2.1 ± 0.6, IQRank = 2.7 ± 0.3, ISRank = 2.2 ± 0.2). We illustrated the modularity of PnP QSM by interchanging two additional patch-based denoisers. PnP QSM reconstruction was feasible, and its flexibility was shown using several patch-based denoisers. This technique may allow rapid prototyping and validation of new denoisers for QSM reconstruction for an array of useful clinical applications.
Collapse
Affiliation(s)
- Srikant Kamesh Iyer
- Department of Radiology, University of Pennsylvania, Philadelphia, PA, USA.
- Perelman Center for Advanced Medicine, South Pavilion, Rm 11-155, Philadelphia, PA, USA.
| | - Brianna F Moon
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA, USA
| | - Nicholas Josselyn
- Department of Radiology, University of Pennsylvania, Philadelphia, PA, USA
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA, USA
| | - Robert M Kurtz
- Department of Radiology, University of Pennsylvania, Philadelphia, PA, USA
| | - Jae W Song
- Department of Radiology, University of Pennsylvania, Philadelphia, PA, USA
| | - Jeffrey B Ware
- Department of Radiology, University of Pennsylvania, Philadelphia, PA, USA
| | - S Ali Nabavizadeh
- Department of Radiology, University of Pennsylvania, Philadelphia, PA, USA
| | - Walter R Witschey
- Department of Radiology, University of Pennsylvania, Philadelphia, PA, USA
| |
Collapse
|
13
|
Boudreau M, Stikov N, Jezzard P. On the open-source landscape of Magnetic Resonance in Medicine. Magn Reson Med 2022; 88:1495-1497. [PMID: 35900923 DOI: 10.1002/mrm.29366] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2022] [Accepted: 05/30/2022] [Indexed: 11/10/2022]
Affiliation(s)
- Mathieu Boudreau
- Montreal Heart Institute, Montreal, Quebec, Canada.,NeuroPoly Lab, Polytechnique Montréal, Montreal, Quebec, Canada
| | - Nikola Stikov
- Montreal Heart Institute, Montreal, Quebec, Canada.,NeuroPoly Lab, Polytechnique Montréal, Montreal, Quebec, Canada
| | - Peter Jezzard
- Wellcome Centre for Integrative Neuroimaging, University of Oxford, Oxford, UK
| |
Collapse
|
14
|
Milovic C, Lambert M, Langkammer C, Bredies K, Irarrazaval P, Tejos C. Streaking artifact suppression of quantitative susceptibility mapping reconstructions via L1-norm data fidelity optimization (L1-QSM). Magn Reson Med 2022; 87:457-473. [PMID: 34350634 DOI: 10.1002/mrm.28957] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2021] [Revised: 07/15/2021] [Accepted: 07/17/2021] [Indexed: 01/05/2023]
Abstract
PURPOSE The presence of dipole-inconsistent data due to substantial noise or artifacts causes streaking artifacts in quantitative susceptibility mapping (QSM) reconstructions. Often used Bayesian approaches rely on regularizers, which in turn yield reduced sharpness. To overcome this problem, we present a novel L1-norm data fidelity approach that is robust with respect to outliers, and therefore prevents streaking artifacts. METHODS QSM functionals are solved with linear and nonlinear L1-norm data fidelity terms using functional augmentation, and are compared with equivalent L2-norm methods. Algorithms were tested on synthetic data, with phase inconsistencies added to mimic lesions, QSM Challenge 2.0 data, and in vivo brain images with hemorrhages. RESULTS The nonlinear L1-norm-based approach achieved the best overall error metric scores and better streaking artifact suppression. Notably, L1-norm methods could reconstruct QSM images without using a brain mask, with similar regularization weights for different data fidelity weighting or masking setups. CONCLUSION The proposed L1-approach provides a robust method to prevent streaking artifacts generated by dipole-inconsistent data, renders brain mask calculation unessential, and opens novel challenging clinical applications such asassessing brain hemorrhages and cortical layers.
Collapse
Affiliation(s)
- Carlos Milovic
- Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Nucleus for Cardiovascular Magnetic Resonance, Santiago, Chile
- Department of Medical Physics and Biomedical Engineering, University College London, London, UK
| | - Mathias Lambert
- Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Nucleus for Cardiovascular Magnetic Resonance, Santiago, Chile
| | - Christian Langkammer
- Department of Neurology, Medical University of Graz, Graz, Austria
- BioTechMed Graz, Graz, Austria
| | - Kristian Bredies
- BioTechMed Graz, Graz, Austria
- Institute of Mathematics and Scientific Computing, University of Graz, Graz, Austria
| | - Pablo Irarrazaval
- Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Nucleus for Cardiovascular Magnetic Resonance, Santiago, Chile
- Institute for Biological and Medical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
| | - Cristian Tejos
- Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Nucleus for Cardiovascular Magnetic Resonance, Santiago, Chile
| |
Collapse
|
15
|
Bilgic B, Langkammer C, Marques JP, Meineke J, Milovic C, Schweser F. QSM reconstruction challenge 2.0: Design and report of results. Magn Reson Med 2021; 86:1241-1255. [PMID: 33783037 DOI: 10.1002/mrm.28754] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2020] [Revised: 01/25/2021] [Accepted: 02/08/2021] [Indexed: 12/24/2022]
Abstract
PURPOSE The aim of the second quantitative susceptibility mapping (QSM) reconstruction challenge (Oct 2019, Seoul, Korea) was to test the accuracy of QSM dipole inversion algorithms in simulated brain data. METHODS A two-stage design was chosen for this challenge. The participants were provided with datasets of multi-echo gradient echo images synthesized from two realistic in silico head phantoms using an MR simulator. At the first stage, participants optimized QSM reconstructions without ground truth data available to mimic the clinical setting. At the second stage, ground truth data were provided for parameter optimization. Submissions were evaluated using eight numerical metrics and visual ratings. RESULTS A total of 98 reconstructions were submitted for stage 1 and 47 submissions for stage 2. Iterative methods had the best quantitative metric scores, followed by deep learning and direct inversion methods. Priors derived from magnitude data improved the metric scores. Algorithms based on iterative approaches and total variation (and its derivatives) produced the best overall results. The reported results and analysis pipelines have been made public to allow researchers to compare new methods to the current state of the art. CONCLUSION The synthetic data provide a consistent framework to test the accuracy and robustness of QSM algorithms in the presence of noise, calcifications and minor voxel dephasing effects. Total Variation-based algorithms produced the best results among all metrics. Future QSM challenges should assess whether this good performance with synthetic datasets translates to more realistic scenarios, where background fields and dipole-incompatible phase contributions are included.
Collapse
Affiliation(s)
- Berkin Bilgic
- Athinoula A. Martinos Center for Biomedical Imaging, Charlestown, Massachusetts, USA
- Department of Radiology, Harvard Medical School, Boston, Massachusetts, USA
- Harvard-MIT Health Sciences and Technology, MIT, Cambridge, Massachusetts, USA
| | | | - José P Marques
- Donders Centre for Cognitive Neuroimaging, Radboud University, Nijmegen, Netherlands
| | | | - Carlos Milovic
- Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Department of Medical Physics and Biomedical Engineering, University College London, London, UK
| | - Ferdinand Schweser
- Buffalo Neuroimaging Analysis Center, Department of Neurology, Jacobs School of Medicine and Biomedical Sciences, University at Buffalo, The State University of New York, Buffalo, New York, USA
- Center for Biomedical Imaging, Clinical and Translational Science Institute, University at Buffalo, The State University of New York, Buffalo, New York, USA
| |
Collapse
|
16
|
Marques JP, Meineke J, Milovic C, Bilgic B, Chan K, Hedouin R, van der Zwaag W, Langkammer C, Schweser F. QSM reconstruction challenge 2.0: A realistic in silico head phantom for MRI data simulation and evaluation of susceptibility mapping procedures. Magn Reson Med 2021; 86:526-542. [PMID: 33638241 PMCID: PMC8048665 DOI: 10.1002/mrm.28716] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2020] [Revised: 01/12/2021] [Accepted: 01/15/2021] [Indexed: 12/14/2022]
Abstract
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.
Collapse
Affiliation(s)
- José P. Marques
- Donders Institute for Brain, Cognition and BehaviorRadboud UniversityNijmegenthe Netherlands
| | | | - Carlos Milovic
- Department of Electrical EngineeringPontificia Universidad Catolica de ChileSantiagoChile
- Biomedical Imaging CenterPontificia Universidad Catolica de ChileSantiagoChile
- Department of Medical Physics and Biomedical EngineeringUniversity College LondonLondonUnited Kingdom
| | - Berkin Bilgic
- Athinoula A. Martinos Center for Biomedical ImagingCharlestownMassachusettsUSA
- Department of RadiologyHarvard Medical SchoolBostonMassachusettsUSA
- Harvard‐MIT Health Sciences and TechnologyMITCambridgeMassachusettsUSA
| | - Kwok‐Shing Chan
- Donders Institute for Brain, Cognition and BehaviorRadboud UniversityNijmegenthe Netherlands
| | - Renaud Hedouin
- Donders Institute for Brain, Cognition and BehaviorRadboud UniversityNijmegenthe Netherlands
- Centre Inria Rennes ‐ Bretagne AtlantiqueRennesFrance
| | | | | | - Ferdinand Schweser
- Buffalo Neuroimaging Analysis CenterDepartment of NeurologyJacobs School of Medicine and Biomedical SciencesUniversity at BuffaloThe State University of New YorkBuffaloNew YorkUSA
- Center for Biomedical Imaging, Clinical and Translational Science InstituteUniversity at BuffaloThe State University of New YorkBuffaloNew YorkUSA
| |
Collapse
|
17
|
Chen L, Cai S, van Zijl PC, Li X. Single-step calculation of susceptibility through multiple orientation sampling. NMR IN BIOMEDICINE 2021; 34:e4517. [PMID: 33822416 PMCID: PMC8184590 DOI: 10.1002/nbm.4517] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/13/2020] [Revised: 03/06/2021] [Accepted: 03/14/2021] [Indexed: 06/12/2023]
Abstract
Quantitative susceptibility mapping (QSM) was developed to estimate the spatial distribution of magnetic susceptibility from MR signal phase acquired using a gradient echo (GRE) sequence. The field-to-susceptibility inversion in QSM is known to be ill-posed and needs numerical stabilization through either regularization or data oversampling. The calculation of susceptibility through the multiple orientation sampling (COSMOS) method uses phase data acquired at three or more head orientations to achieve a well-conditioned field-to-susceptibility inversion and is often considered the gold standard for in vivo QSM. However, the conventional COSMOS approach, here named multistep COSMOS (MSCOSMOS), solves the dipole inversion from the local field derived from raw GRE phase through multiple steps of phase preprocessing. Error propagations between these consecutive phase processing steps can thus affect the final susceptibility quantification. On the other hand, recently proposed single-step QSM (SSQSM) methods aim to solve an integrated inversion from unprocessed or total phase to mitigate such error propagations but have been limited to single orientation QSM. This study therefore aimed to test the feasibility of using single-step COSMOS (SSCOSMOS) to jointly perform background field removal and dipole inversion with multiple orientation sampling, which could serve as a better standard for gauging SSQSM methods. We incorporated multiple spherical mean value (SMV) kernels of various radii with the dipole inversion in SSCOSMOS. QSM reconstructions with SSCOSMOS and MSCOSMOS were compared using both simulations with a numerical head phantom and in vivo human brain data. SSCOSMOS permitted integrated background removal and dipole inversion without the need to adjust any regularization parameters. In addition, with sufficiently large SMV kernels, SSCOSMOS performed consistently better than MSCOSMOS in all the tested error metrics in our simulations, giving better susceptibility quantification and smaller reconstruction error. Consistent tissue susceptibility values were obtained between SSCOSMOS and MSCOSMOS.
Collapse
Affiliation(s)
- Lin Chen
- Department of Electronic Science, Fujian Provincial Key Laboratory of Plasma and Magnetic Resonance, Xiamen University, Xiamen, China
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, United States
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, Maryland, United States
| | - Shuhui Cai
- Department of Electronic Science, Fujian Provincial Key Laboratory of Plasma and Magnetic Resonance, Xiamen University, Xiamen, China
| | - Peter C.M. van Zijl
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, United States
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, Maryland, United States
| | - Xu Li
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, United States
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, Maryland, United States
| |
Collapse
|
18
|
Milovic C, Prieto C, Bilgic B, Uribe S, Acosta-Cabronero J, Irarrazaval P, Tejos C. Comparison of parameter optimization methods for quantitative susceptibility mapping. Magn Reson Med 2021; 85:480-494. [PMID: 32738103 PMCID: PMC7722179 DOI: 10.1002/mrm.28435] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2020] [Revised: 06/10/2020] [Accepted: 06/26/2020] [Indexed: 12/28/2022]
Abstract
PURPOSE Quantitative Susceptibility Mapping (QSM) is usually performed by minimizing a functional with data fidelity and regularization terms. A weighting parameter controls the balance between these terms. There is a need for techniques to find the proper balance that avoids artifact propagation and loss of details. Finding the point of maximum curvature in the L-curve is a popular choice, although it is slow, often unreliable when using variational penalties, and has a tendency to yield overregularized results. METHODS We propose 2 alternative approaches to control the balance between the data fidelity and regularization terms: 1) searching for an inflection point in the log-log domain of the L-curve, and 2) comparing frequency components of QSM reconstructions. We compare these methods against the conventional L-curve and U-curve approaches. RESULTS Our methods achieve predicted parameters that are better correlated with RMS error, high-frequency error norm, and structural similarity metric-based parameter optimizations than those obtained with traditional methods. The inflection point yields less overregularization and lower errors than traditional alternatives. The frequency analysis yields more visually appealing results, although with larger RMS error. CONCLUSION Our methods provide a robust parameter optimization framework for variational penalties in QSM reconstruction. The L-curve-based zero-curvature search produced almost optimal results for typical QSM acquisition settings. The frequency analysis method may use a 1.5 to 2.0 correction factor to apply it as a stand-alone method for a wider range of signal-to-noise-ratio settings. This approach may also benefit from fast search algorithms such as the binary search to speed up the process.
Collapse
Affiliation(s)
- Carlos Milovic
- Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Nucleus for Cardiovascular Magnetic Resonance, Santiago, Chile
| | - Claudia Prieto
- Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, UK
| | - Berkin Bilgic
- Athinoula A. Martinos Center for Biomedical Imaging, Charlestown, MA, USA
- Department of Radiology, Harvard Medical School, Boston, MA, USA
- Harvard-MIT Health Sciences and Technology, MIT, Cambridge, MA, USA
| | - Sergio Uribe
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Nucleus for Cardiovascular Magnetic Resonance, Santiago, Chile
- Department of Radiology, Pontificia Universidad Catolica de Chile, Santiago, Chile
| | | | - Pablo Irarrazaval
- Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Nucleus for Cardiovascular Magnetic Resonance, Santiago, Chile
- Institute for Biological and Medical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
| | - Cristian Tejos
- Department of Electrical Engineering, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Nucleus for Cardiovascular Magnetic Resonance, Santiago, Chile
| |
Collapse
|