1
|
Kung HT, Cui SX, Kaplan JT, Joshi AA, Leahy RM, Nayak KS, Haldar JP. Diffusion tensor brain imaging at 0.55T: A feasibility study. Magn Reson Med 2024; 92:1649-1657. [PMID: 38725132 DOI: 10.1002/mrm.30156] [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: 11/06/2023] [Revised: 04/09/2024] [Accepted: 04/28/2024] [Indexed: 07/23/2024]
Abstract
PURPOSE To investigate the feasibility of diffusion tensor brain imaging at 0.55T with comparisons against 3T. METHODS Diffusion tensor imaging data with 2 mm isotropic resolution was acquired on a cohort of five healthy subjects using both 0.55T and 3T scanners. The signal-to-noise ratio (SNR) of the 0.55T data was improved using a previous SNR-enhancing joint reconstruction method that jointly reconstructs the entire set of diffusion weighted images from k-space using shared-edge constraints. Quantitative diffusion tensor parameters were estimated and compared across field strengths. We also performed a test-retest assessment of repeatability at each field strength. RESULTS After applying SNR-enhancing joint reconstruction, the diffusion tensor parameters obtained from 0.55T data were strongly correlated (R 2 ≥ 0 . 70 $$ {R}^2\ge 0.70 $$ ) with those obtained from 3T data. Test-retest analysis showed that SNR-enhancing reconstruction improved the repeatability of the 0.55T diffusion tensor parameters. CONCLUSION High-resolution in vivo diffusion MRI of the human brain is feasible at 0.55T when appropriate noise-mitigation strategies are applied.
Collapse
Affiliation(s)
- Hao-Ting Kung
- Signal and Image Processing Institute, Ming Hsieh Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, California, USA
| | - Sophia X Cui
- Siemens Medical Solutions USA, Los Angeles, California, USA
| | - Jonas T Kaplan
- Brain and Creativity Institute, University of Southern California, Los Angeles, California, USA
| | - Anand A Joshi
- Signal and Image Processing Institute, Ming Hsieh Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, California, USA
| | - Richard M Leahy
- Signal and Image Processing Institute, Ming Hsieh Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, California, USA
| | - Krishna S Nayak
- Signal and Image Processing Institute, Ming Hsieh Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, California, USA
| | - Justin P Haldar
- Signal and Image Processing Institute, Ming Hsieh Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, California, USA
- Brain and Creativity Institute, University of Southern California, Los Angeles, California, USA
| |
Collapse
|
2
|
Berglund J, Liljeblad M, Baron T. Unwrapping phase contrast MRI by iterative graph cuts. Magn Reson Med 2024; 92:1484-1495. [PMID: 38725423 DOI: 10.1002/mrm.30138] [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: 11/24/2023] [Revised: 04/10/2024] [Accepted: 04/15/2024] [Indexed: 07/23/2024]
Abstract
PURPOSE To develop and evaluate a phase unwrapping method for cine phase contrast MRI based on graph cuts. METHODS A proposed Iterative Graph Cuts method was evaluated in 10 cardiac patients with two-dimensional flow quantification which was repeated at low venc settings to provoke wrapping. The images were also unwrapped by a path-following method (ROMEO), and a Laplacian-based method (LP). Net flow was quantified using semi-automatic vessel segmentation. High venc images were also wrapped retrospectively to asses the residual amount of wrapped voxels. RESULTS The absolute net flow error after unwrapping at venc = 100 cm/s was 1.8 mL, which was 0.83 mL smaller than for LP. The repeatability error at high venc without unwrapping was 2.5 mL. The error at venc = 50 cm/s was 7.5 mL, which was 8.2 mL smaller than for ROMEO and 5.7 mL smaller than for LP. For retrospectively wrapped images with synthetic venc of 100/50/25 cm/s, the residual amount of wrapped voxels was 0.00/0.12/0.79%, which was 0.09/0.26/8.0 percentage points smaller than for LP. With synthetic venc of 25 cm/s, omitting magnitude information resulted in 3.2 percentage points more wrapped voxels, and only spatial/temporal unwrapping resulted in 4.6/21 percentage points more wrapped voxels compared to spatiotemporal unwrapping. CONCLUSION Iterative Graph Cuts enables unwrapping of cine phase contrast MRI with very small errors, except for at extreme blood velocities, with equal or better performance compared to ROMEO and LP. The use of magnitude information and spatiotemporal unwrapping is recommended.
Collapse
Affiliation(s)
- Johan Berglund
- Medical Physics, Uppsala University Hospital, Uppsala, Sweden
- Molecular Imaging and Medical Physics, Department of Surgical Sciences, Uppsala University, Uppsala, Sweden
| | - Mio Liljeblad
- Medical Physics, Uppsala University Hospital, Uppsala, Sweden
| | - Tomasz Baron
- Cardiology and Clinical Physiology, Department of Medical Sciences, Uppsala University, Uppsala, Sweden
- Uppsala Clinical Research Center, Uppsala University, Uppsala, Sweden
| |
Collapse
|
3
|
Ji S, Jang J, Kim M, Lee H, Kim W, Lee J, Shin HG. Comparison between R2'-based and R2*-based χ-separation methods: A clinical evaluation in individuals with multiple sclerosis. NMR IN BIOMEDICINE 2024; 37:e5167. [PMID: 38697612 DOI: 10.1002/nbm.5167] [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: 12/09/2023] [Revised: 03/28/2024] [Accepted: 04/02/2024] [Indexed: 05/05/2024]
Abstract
Susceptibility source separation, or χ-separation, estimates diamagnetic (χdia) and paramagnetic susceptibility (χpara) signals in the brain using local field and R2' (= R2* - R2) maps. Recently proposed R2*-based χ-separation methods allow for χ-separation using only multi-echo gradient echo (ME-GRE) data, eliminating the need for additional data acquisition for R2 mapping. Although this approach reduces scan time and enhances clinical utility, the impact of missing R2 information remains a subject of exploration. In this study, we evaluate the viability of two previously proposed R2*-based χ-separation methods as alternatives to their R2'-based counterparts: model-based R2*-χ-separation versus χ-separation and deep learning-based χ-sepnet-R2* versus χ-sepnet-R2'. Their performances are assessed in individuals with multiple sclerosis (MS), comparing them with their corresponding R2'-based counterparts (i.e., R2*-χ-separation vs. χ-separation and χ-sepnet-R2* vs. χ-sepnet-R2'). The evaluations encompass qualitative visual assessments by experienced neuroradiologists and quantitative analyses, including region of interest analyses and linear regression analyses. Qualitatively, R2*-χ-separation tends to report higher χpara and χdia values compared with χ-separation, leading to less distinct lesion contrasts, while χ-sepnet-R2* closely aligns with χ-sepnet-R2'. Quantitative analysis reveals a robust correlation between both R2*-based methods and their R2'-based counterparts (r ≥ 0.88). Specifically, in the whole-brain voxels, χ-sepnet-R2* exhibits higher correlation and better linearity than R2*-χ-separation (χdia/χpara from R2*-χ-separation: r = 0.88/0.90, slope = 0.79/0.86; χdia/χpara from χ-sepnet-R2*: r = 0.90/0.92, slope = 0.99/0.97). In MS lesions, both R2*-based methods display comparable correlation and linearity (χdia/χpara from R2*-χ-separation: r = 0.90/0.91, slope = 0.98/0.91; χdia/χpara from χ-sepnet-R2*: r = 0.88/0.88, slope = 0.91/0.95). Notably, χ-sepnet-R2* demonstrates negligible offsets, whereas R2*-χ-separation exhibits relatively large offsets (0.02 ppm in the whole brain and 0.01 ppm in the MS lesions), potentially indicating the false presence of myelin or iron in MS lesions. Overall, both R2*-based χ-separation methods demonstrated their viability as alternatives to their R2'-based counterparts. χ-sepnet-R2* showed better alignment with its R2'-based counterpart with minimal susceptibility offsets, compared with R2*-χ-separation that reported higher χpara and χdia values compared with R2'-based χ-separation.
Collapse
Affiliation(s)
- Sooyeon Ji
- Department of Electrical and Computer Engineering, Seoul National University, Seoul, South Korea
| | - Jinhee Jang
- Department of Radiology, Seoul St. Mary's Hospital, College of Medicine, The Catholic University of Korea, Seoul, South Korea
| | - Minjun Kim
- Department of Electrical and Computer Engineering, Seoul National University, Seoul, South Korea
| | - Hyebin Lee
- Department of Radiology, Seoul St. Mary's Hospital, College of Medicine, The Catholic University of Korea, Seoul, South Korea
| | - Woojun Kim
- Department of Neurology, Seoul St. Mary's Hospital, College of Medicine, The Catholic University of Korea, Seoul, South Korea
| | - Jongho Lee
- Department of Electrical and Computer Engineering, Seoul National University, Seoul, South Korea
| | - Hyeong-Geol Shin
- Department of Radiology, School of Medicine, Johns Hopkins University, Baltimore, Maryland, USA
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, USA
| |
Collapse
|
4
|
Rojas-López JA, Cabrera-Santiago A, Adragna C, Ibarra-Ortega BE, López-Luna JE, Contreras-Rodríguez JA, Martínez-Ortiz E. Commissioning of MRI-guided gynaecological brachytherapy using an MR-linac. Biomed Phys Eng Express 2024; 10:055032. [PMID: 39111326 DOI: 10.1088/2057-1976/ad6c54] [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: 05/30/2024] [Accepted: 08/07/2024] [Indexed: 08/28/2024]
Abstract
Purpose. To evaluate the feasibility of use of an 1.5 T magnetic resonance (MR)-linear accelerator MR-linac for imaging in gynaecologic high-dose-rate (HDR) brachytherapy.Method. Commissioning measurements for MR images quality control, geometric distortion, dwell position accuracy, applicator reconstruction and end-to-end test for a tandem-and-ring applicator were performed following the recommendations of American Brachytherapy Society, International Commission on Radiation Units and Measurements and Report of the Brachytherapy Working Group of the Spanish Society of Medical Physics. The values for MR-based IGABT were compared to the corresponding values with computed tomography (CT).Results. Measured distorsions for the MR images were less than 0.50 mm compared to the CT images. The differences between 3D displacements for all dwell positions were 0.66 mm and 0.62 mm for the tandem and ring, respectively. The maximum difference is 0.64 mm for the distances from the applicator tip obtained using the films. The CT and MR dose differences for the right and left 'A' points were 0.9% and -0.7%, respectively. Similar results were observed in terms of dose distribution for CT and Mr The gamma passing rate was 99.3% and 99.5%, respectively.Conclusion. The use of MR images from an MR-linac used in a radiotherapy service for gynaecological brachytherapy was proved to be feasible, safe and precise as the geometrical differences were less than 1 mm, and the dosimetric differences were less than 1% when comparing to the use of CT images for the same purpose.
Collapse
Affiliation(s)
- José Alejandro Rojas-López
- Hospital Almater SA de CV, Álvaro Obregón 1100, Segunda Sección, Mexicali, Baja California, Mexico
- Facultad de Astronomía, Matemáticas, Física y Computación, Universidad Nacional de Córdoba, Av. Medina Allende s/n, Ciudad Universitaria, CP:X5000HUA, Córdoba, Argentina
| | - Alexis Cabrera-Santiago
- Hospital Almater SA de CV, Álvaro Obregón 1100, Segunda Sección, Mexicali, Baja California, Mexico
- Unidad de Especialidades Médicas de Oncología, Av Claridad, Plutarco Elías Calles, 21376, Mexicali, Baja California, Mexico
| | - Celeste Adragna
- Instituto Oulton, Av Vélez Sarsfield 652, Córdoba, Argentina
| | | | - José Eleazar López-Luna
- Hospital Almater SA de CV, Álvaro Obregón 1100, Segunda Sección, Mexicali, Baja California, Mexico
| | | | - Efraín Martínez-Ortiz
- Hospital Almater SA de CV, Álvaro Obregón 1100, Segunda Sección, Mexicali, Baja California, Mexico
- Unidad de Especialidades Médicas de Oncología, Av Claridad, Plutarco Elías Calles, 21376, Mexicali, Baja California, Mexico
| |
Collapse
|
5
|
Naji N, Gee M, Jickling GC, Emery DJ, Saad F, McCreary CR, Smith EE, Camicioli R, Wilman AH. Quantifying cerebral microbleeds using quantitative susceptibility mapping from magnetization-prepared rapid gradient-echo. NMR IN BIOMEDICINE 2024; 37:e5139. [PMID: 38465729 DOI: 10.1002/nbm.5139] [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: 02/07/2024] [Accepted: 02/13/2024] [Indexed: 03/12/2024]
Abstract
T1-weighted magnetization-prepared rapid gradient-echo (MPRAGE) is commonly included in brain studies for structural imaging using magnitude images; however, its phase images can provide an opportunity to assess microbleed burden using quantitative susceptibility mapping (QSM). This potential application for MPRAGE-based QSM was evaluated using in vivo and simulated measurements. Possible factors affecting image quality were also explored. Detection sensitivity was evaluated against standard multiecho gradient echo (MEGE) QSM using 3-T in vivo data of 15 subjects with a combined total of 108 confirmed microbleeds. The two methods were compared based on the microbleed size and susceptibility measurements. In addition, simulations explored the detection sensitivity of MPRAGE-QSM at different representative magnetic field strengths and echo times using microbleeds of different size, susceptibility, and location. Results showed that in vivo microbleeds appeared to be smaller (× 0.54) and of higher mean susceptibility (× 1.9) on MPRAGE-QSM than on MEGE-QSM, but total susceptibility estimates were in closer agreement (slope: 0.97, r2: 0.94), and detection sensitivity was comparable. In simulations, QSM at 1.5 T had a low contrast-to-noise ratio that obscured the detection of many microbleeds. Signal-to-noise ratio (SNR) levels at 3 T and above resulted in better contrast and increased detection. The detection rates for microbleeds of minimum one-voxel diameter and 0.4-ppm susceptibility were 0.55, 0.80, and 0.88 at SNR levels of 1.5, 3, and 7 T, respectively. Size and total susceptibility estimates were more consistent than mean susceptibility estimates, which showed size-dependent underestimation. MPRAGE-QSM provides an opportunity to detect and quantify the size and susceptibility of microbleeds of at least one-voxel diameter at B0 of 3 T or higher with no additional time cost, when standard T2*-weighted images are not available or have inadequate spatial resolution. The total susceptibility measure is more robust against sequence variations and might allow combining data from different protocols.
Collapse
Affiliation(s)
- Nashwan Naji
- Department of Biomedical Engineering, University of Alberta, Edmonton, Alberta, Canada
| | - Myrlene Gee
- Division of Neurology, University of Alberta, Edmonton, Alberta, Canada
| | - Glen C Jickling
- Division of Neurology, University of Alberta, Edmonton, Alberta, Canada
| | - Derek J Emery
- Department of Radiology and Diagnostic Imaging, University of Alberta, Edmonton, Alberta, Canada
| | - Feryal Saad
- Radiology and Clinical Neurosciences, Hotchkiss Brain Institute, University of Calgary, Calgary, Alberta, Canada
| | - Cheryl R McCreary
- Radiology and Clinical Neurosciences, Hotchkiss Brain Institute, University of Calgary, Calgary, Alberta, Canada
- Seaman Family MR Research Centre, Foothills Medical Centre, Calgary, Alberta, Canada
| | - Eric E Smith
- Radiology and Clinical Neurosciences, Hotchkiss Brain Institute, University of Calgary, Calgary, Alberta, Canada
| | - Richard Camicioli
- Division of Neurology, University of Alberta, Edmonton, Alberta, Canada
| | - Alan H Wilman
- Department of Biomedical Engineering, University of Alberta, Edmonton, Alberta, Canada
- Department of Radiology and Diagnostic Imaging, University of Alberta, Edmonton, Alberta, Canada
| |
Collapse
|
6
|
Hervouin A, Bézy-Wendling J, Noury F. How to accurately quantify brain magnetic susceptibility in the context of Parkinson's disease: Validation on phantoms and healthy volunteers at 1.5 and 3 T. NMR IN BIOMEDICINE 2024:e5182. [PMID: 38993048 DOI: 10.1002/nbm.5182] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/12/2024] [Revised: 05/06/2024] [Accepted: 05/06/2024] [Indexed: 07/13/2024]
Abstract
Currently, brain iron content represents a new neuromarker for understanding the physiopathological mechanisms leading to Parkinson's disease (PD). In vivo quantification of biological iron is possible by reconstructing magnetic susceptibility maps obtained using quantitative susceptibility mapping (QSM). Applying QSM is challenging, as up to now, no standardization of acquisition protocols and phase image processing has emerged from referenced studies. Our objectives were to compare the accuracy and the sensitivity of 10 QSM pipelines built from algorithms from the literature, applied on phantoms data and on brain data. Two phantoms, with known magnetic susceptibility ranges, were created from several solutions of gadolinium chelate. Twenty healthy volunteers from two age groups were included. Phantoms and brain data were acquired at 1.5 and 3 T, respectively. Susceptibility-weighted images were obtained using a 3D multigradient-recalled-echo sequence. For brain data, 3D anatomical T1- and T2-weighted images were also acquired to segment the deep gray nuclei of interest. Concerning in vitro data, the linear dependence of magnetic susceptibility versus gadolinium concentration and deviations from the theoretically expected values were calculated. For brain data, the accuracy and sensitivity of the QSM pipelines were evaluated in comparison with results from the literature and regarding the expected magnetic susceptibility increase with age, respectively. A nonparametric Mann-Whitney U-test was used to compare the magnetic susceptibility quantification in deep gray nuclei between the two age groups. Our methodology enabled quantifying magnetic susceptibility in human brain and the results were consistent with those from the literature. Statistically significant differences were obtained between the two age groups in all cerebral regions of interest. Our results show the importance of optimizing QSM pipelines according to the application and the targeted magnetic susceptibility range, to achieve accurate quantification. We were able to define the optimal QSM pipeline for future applications on patients with PD.
Collapse
Affiliation(s)
| | | | - Fanny Noury
- Univ Rennes, Inserm, LTSI-UMR 1099, Rennes, France
| |
Collapse
|
7
|
Karakuzu A, Boudreau M, Stikov N. Reproducible Research Practices in Magnetic Resonance Neuroimaging: A Review Informed by Advanced Language Models. Magn Reson Med Sci 2024; 23:252-267. [PMID: 38897936 PMCID: PMC11234949 DOI: 10.2463/mrms.rev.2023-0174] [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] [Indexed: 06/21/2024] Open
Abstract
MRI has progressed significantly with the introduction of advanced computational methods and novel imaging techniques, but their wider adoption hinges on their reproducibility. This concise review synthesizes reproducible research insights from recent MRI articles to examine the current state of reproducibility in neuroimaging, highlighting key trends and challenges. It also provides a custom generative pretrained transformer (GPT) model, designed specifically for aiding in an automated analysis and synthesis of information pertaining to the reproducibility insights associated with the articles at the core of this review.
Collapse
Affiliation(s)
- Agah Karakuzu
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montréal, Montréal, Quebec, Canada
- Montréal Heart Institute, Montréal, Quebec, Canada
| | - Mathieu Boudreau
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montréal, Montréal, Quebec, Canada
| | - Nikola Stikov
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montréal, Montréal, Quebec, Canada
- Montréal Heart Institute, Montréal, Quebec, Canada
- Center for Advanced Interdisciplinary Research, Ss. Cyril and Methodius University, Skopje, North Macedonia
| |
Collapse
|
8
|
Jansen MG, Zwiers MP, Marques JP, Chan KS, Amelink JS, Altgassen M, Oosterman JM, Norris DG. The Advanced BRain Imaging on ageing and Memory (ABRIM) data collection: Study design, data processing, and rationale. PLoS One 2024; 19:e0306006. [PMID: 38905233 PMCID: PMC11192316 DOI: 10.1371/journal.pone.0306006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2024] [Accepted: 06/07/2024] [Indexed: 06/23/2024] Open
Abstract
To understand the neurocognitive mechanisms that underlie heterogeneity in cognitive ageing, recent scientific efforts have led to a growing public availability of imaging cohort data. The Advanced BRain Imaging on ageing and Memory (ABRIM) project aims to add to these existing datasets by taking an adult lifespan approach to provide a cross-sectional, normative database with a particular focus on connectivity, myelinization and iron content of the brain in concurrence with cognitive functioning, mechanisms of reserve, and sleep-wake rhythms. ABRIM freely shares MRI and behavioural data from 295 participants between 18-80 years, stratified by age decade and sex (median age 52, IQR 36-66, 53.20% females). The ABRIM MRI collection consists of both the raw and pre-processed structural and functional MRI data to facilitate data usage among both expert and non-expert users. The ABRIM behavioural collection includes measures of cognitive functioning (i.e., global cognition, processing speed, executive functions, and memory), proxy measures of cognitive reserve (e.g., educational attainment, verbal intelligence, and occupational complexity), and various self-reported questionnaires (e.g., on depressive symptoms, pain, and the use of memory strategies in daily life and during a memory task). In a sub-sample (n = 120), we recorded sleep-wake rhythms using an actigraphy device (Actiwatch 2, Philips Respironics) for a period of 7 consecutive days. Here, we provide an in-depth description of our study protocol, pre-processing pipelines, and data availability. ABRIM provides a cross-sectional database on healthy participants throughout the adult lifespan, including numerous parameters relevant to improve our understanding of cognitive ageing. Therefore, ABRIM enables researchers to model the advanced imaging parameters and cognitive topologies as a function of age, identify the normal range of values of such parameters, and to further investigate the diverse mechanisms of reserve and resilience.
Collapse
Affiliation(s)
- Michelle G. Jansen
- Donders Centre for Cognition, Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, the Netherlands
| | - Marcel P. Zwiers
- Donders Centre for Cognitive Neuroimaging, Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, the Netherlands
| | - Jose P. Marques
- Donders Centre for Cognitive Neuroimaging, Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, the Netherlands
| | - Kwok-Shing Chan
- Donders Centre for Cognitive Neuroimaging, Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, the Netherlands
| | - Jitse S. Amelink
- Donders Centre for Cognitive Neuroimaging, Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, the Netherlands
- Language and Genetics Department, Max Planck Institute for Psycholinguistics, Radboud University, Nijmegen, the Netherlands
| | - Mareike Altgassen
- Department of Psychology, Johannes Gutenberg-University Mainz, Mainz, Germany
| | - Joukje M. Oosterman
- Donders Centre for Cognition, Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, the Netherlands
| | - David G. Norris
- Donders Centre for Cognitive Neuroimaging, Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, the Netherlands
| |
Collapse
|
9
|
He Z, Lefebvre PM, Soullié P, Doguet M, Ambarki K, Chen B, Odille F. Phantom evaluation of electrical conductivity mapping by MRI: Comparison to vector network analyzer measurements and spatial resolution assessment. Magn Reson Med 2024; 91:2374-2390. [PMID: 38225861 DOI: 10.1002/mrm.30009] [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: 11/10/2023] [Revised: 12/21/2023] [Accepted: 12/22/2023] [Indexed: 01/17/2024]
Abstract
PURPOSE To evaluate the performance of various MR electrical properties tomography (MR-EPT) methods at 3 T in terms of absolute quantification and spatial resolution limit for electrical conductivity. METHODS Absolute quantification as well as spatial resolution performance were evaluated on homogeneous phantoms and a phantom with holes of different sizes, respectively. Ground-truth conductivities were measured with an open-ended coaxial probe connected to a vector network analyzer (VNA). Four widely used MR-EPT reconstruction methods were investigated: phase-based Helmholtz (PB), phase-based convection-reaction (PB-cr), image-based (IB), and generalized-image-based (GIB). These methods were compared using the same complex images from a 1 mm-isotropic UTE sequence. Alternative transceive phase acquisition sequences were also compared in PB and PB-cr. RESULTS In large homogeneous phantoms, all methods showed a strong correlation with ground truth conductivities (r > 0.99); however, GIB was the best in terms of accuracy, spatial uniformity, and robustness to boundary artifacts. In the resolution phantom, the normalized root-mean-squared error of all methods grew rapidly (>0.40) when the hole size was below 10 mm, with simplified methods (PB and IB), or below 5 mm, with generalized methods (PB-cr and GIB). CONCLUSION VNA measurements are essential to assess the accuracy of MR-EPT. In this study, all tested MR-EPT methods correlated strongly with the VNA measurements. The UTE sequence is recommended for MR-EPT, with the GIB method providing good accuracy for structures down to 5 mm. Structures below 5 mm may still be detected in the conductivity maps, but with significantly lower accuracy.
Collapse
Affiliation(s)
- Zhongzheng He
- IADI U1254, INSERM and Université de Lorraine, Nancy, France
| | | | - Paul Soullié
- IADI U1254, INSERM and Université de Lorraine, Nancy, France
| | - Martin Doguet
- IADI U1254, INSERM and Université de Lorraine, Nancy, France
- BioSerenity, Paris, France
| | | | - Bailiang Chen
- IADI U1254, INSERM and Université de Lorraine, Nancy, France
- CIC-IT 1433, INSERM, Université de Lorraine and CHRU Nancy, Nancy, France
| | - Freddy Odille
- IADI U1254, INSERM and Université de Lorraine, Nancy, France
- CIC-IT 1433, INSERM, Université de Lorraine and CHRU Nancy, Nancy, France
| |
Collapse
|
10
|
Motyka S, Weiser P, Bachrata B, Hingerl L, Strasser B, Hangel G, Niess E, Niess F, Zaitsev M, Robinson SD, Langs G, Trattnig S, Bogner W. Predicting dynamic, motion-related changes in B 0 field in the brain at a 7T MRI using a subject-specific fine-trained U-net. Magn Reson Med 2024; 91:2044-2056. [PMID: 38193276 DOI: 10.1002/mrm.29980] [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/04/2023] [Revised: 11/28/2023] [Accepted: 11/30/2023] [Indexed: 01/10/2024]
Abstract
PURPOSE Subject movement during the MR examination is inevitable and causes not only image artifacts but also deteriorates the homogeneity of the main magnetic field (B0 ), which is a prerequisite for high quality data. Thus, characterization of changes to B0 , for example induced by patient movement, is important for MR applications that are prone to B0 inhomogeneities. METHODS We propose a deep learning based method to predict such changes within the brain from the change of the head position to facilitate retrospective or even real-time correction. A 3D U-net was trained on in vivo gradient-echo brain 7T MRI data. The input consisted of B0 maps and anatomical images at an initial position, and anatomical images at a different head position (obtained by applying a rigid-body transformation on the initial anatomical image). The output consisted of B0 maps at the new head positions. We further fine-trained the network weights to each subject by measuring a limited number of head positions of the given subject, and trained the U-net with these data. RESULTS Our approach was compared to established dynamic B0 field mapping via interleaved navigators, which suffer from limited spatial resolution and the need for undesirable sequence modifications. Qualitative and quantitative comparison showed similar performance between an interleaved navigator-equivalent method and proposed method. CONCLUSION It is feasible to predict B0 maps from rigid subject movement and, when combined with external tracking hardware, this information could be used to improve the quality of MR acquisitions without the use of navigators.
Collapse
Affiliation(s)
- Stanislav Motyka
- High Field MR Center, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
- Christian Doppler Laboratory for Clinical Molecular MR Imaging, Vienna, Austria
| | - Paul Weiser
- Computational Imaging Research Lab, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Boston, Massachusetts, USA
| | - Beata Bachrata
- Department of Medical Engineering, Carinthia University of Applied Sciences, Klagenfurt, Austria
| | - Lukas Hingerl
- High Field MR Center, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
| | - Bernhard Strasser
- High Field MR Center, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
| | - Gilbert Hangel
- High Field MR Center, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
- Department of Neurosurgery, Medical University of Vienna, Vienna, Austria
| | - Eva Niess
- High Field MR Center, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
- Christian Doppler Laboratory for Clinical Molecular MR Imaging, Vienna, Austria
| | - Fabian Niess
- High Field MR Center, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
| | - Maxim Zaitsev
- Department of Radiology - Medical Physics, University of Freiburg, Freiburg, Germany
- Faculty of Medicine, University of Freiburg - Medical Centre, Freiburg, Germany
| | - Simon Daniel Robinson
- High Field MR Center, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
| | - Georg Langs
- Computational Imaging Research Lab, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
| | - Siegfried Trattnig
- High Field MR Center, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
| | - Wolfgang Bogner
- High Field MR Center, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
- Christian Doppler Laboratory for Clinical Molecular MR Imaging, Vienna, Austria
| |
Collapse
|
11
|
Yang H, Hong K, Baraboo JJ, Fan L, Larsen A, Markl M, Robinson JD, Rigsby CK, Kim D. GRASP reconstruction amplified with view-sharing and KWIC filtering reduces underestimation of peak velocity in highly-accelerated real-time phase-contrast MRI: A preliminary evaluation in pediatric patients with congenital heart disease. Magn Reson Med 2024; 91:1965-1977. [PMID: 38084397 PMCID: PMC10950531 DOI: 10.1002/mrm.29974] [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/10/2023] [Revised: 10/27/2023] [Accepted: 11/27/2023] [Indexed: 02/01/2024]
Abstract
PURPOSE To develop a highly-accelerated, real-time phase contrast (rtPC) MRI pulse sequence with 40 fps frame rate (25 ms effective temporal resolution). METHODS Highly-accelerated golden-angle radial sparse parallel (GRASP) with over regularization may result in temporal blurring, which in turn causes underestimation of peak velocity. Thus, we amplified GRASP performance by synergistically combining view-sharing (VS) and k-space weighted image contrast (KWIC) filtering. In 17 pediatric patients with congenital heart disease (CHD), the conventional GRASP and the proposed GRASP amplified by VS and KWIC (or GRASP + VS + KWIC) reconstruction for rtPC MRI were compared with respect to clinical standard PC MRI in measuring hemodynamic parameters (peak velocity, forward volume, backward volume, regurgitant fraction) at four locations (aortic valve, pulmonary valve, left and right pulmonary arteries). RESULTS The proposed reconstruction method (GRASP + VS + KWIC) achieved better effective spatial resolution (i.e., image sharpness) compared with conventional GRASP, ultimately reducing the underestimation of peak velocity from 17.4% to 6.4%. The hemodynamic metrics (peak velocity, volumes) were not significantly (p > 0.99) different between GRASP + VS + KWIC and clinical PC, whereas peak velocity was significantly (p < 0.007) lower for conventional GRASP. RtPC with GRASP + VS + KWIC also showed the ability to assess beat-to-beat variation and detect the highest peak among peaks. CONCLUSION The synergistic combination of GRASP, VS, and KWIC achieves 25 ms effective temporal resolution (40 fps frame rate), while minimizing the underestimation of peak velocity compared with conventional GRASP.
Collapse
Affiliation(s)
- Huili Yang
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
- Department of Biomedical Engineering, McCormick School of Engineering, Northwestern University, Evanston, Illinois, USA
| | - KyungPyo Hong
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
| | - Justin J Baraboo
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
- Department of Biomedical Engineering, McCormick School of Engineering, Northwestern University, Evanston, Illinois, USA
| | - Lexiaozi Fan
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
- Department of Biomedical Engineering, McCormick School of Engineering, Northwestern University, Evanston, Illinois, USA
| | - Andrine Larsen
- Department of Biomedical Engineering, Lehigh University, Bethlehem, Pennsylvania, USA
| | - Michael Markl
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
- Department of Biomedical Engineering, McCormick School of Engineering, Northwestern University, Evanston, Illinois, USA
| | - Joshua D Robinson
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
- Department of Pediatrics, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
- Division of Cardiology, Ann & Robert H. Lurie Children's Hospital, Chicago, Illinois, USA
| | - Cynthia K Rigsby
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
- Department of Pediatrics, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
- Department of Medical Imaging, Ann & Robert H. Lurie Children's Hospital, Chicago, Illinois, USA
| | - Daniel Kim
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
- Department of Biomedical Engineering, McCormick School of Engineering, Northwestern University, Evanston, Illinois, USA
| |
Collapse
|
12
|
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
|
13
|
Renton AI, Dao TT, Johnstone T, Civier O, Sullivan RP, White DJ, Lyons P, Slade BM, Abbott DF, Amos TJ, Bollmann S, Botting A, Campbell MEJ, Chang J, Close TG, Dörig M, Eckstein K, Egan GF, Evas S, Flandin G, Garner KG, Garrido MI, Ghosh SS, Grignard M, Halchenko YO, Hannan AJ, Heinsfeld AS, Huber L, Hughes ME, Kaczmarzyk JR, Kasper L, Kuhlmann L, Lou K, Mantilla-Ramos YJ, Mattingley JB, Meier ML, Morris J, Narayanan A, Pestilli F, Puce A, Ribeiro FL, Rogasch NC, Rorden C, Schira MM, Shaw TB, Sowman PF, Spitz G, Stewart AW, Ye X, Zhu JD, Narayanan A, Bollmann S. Neurodesk: an accessible, flexible and portable data analysis environment for reproducible neuroimaging. Nat Methods 2024; 21:804-808. [PMID: 38191935 PMCID: PMC11180540 DOI: 10.1038/s41592-023-02145-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2023] [Accepted: 11/27/2023] [Indexed: 01/10/2024]
Abstract
Neuroimaging research requires purpose-built analysis software, which is challenging to install and may produce different results across computing environments. The community-oriented, open-source Neurodesk platform ( https://www.neurodesk.org/ ) harnesses a comprehensive and growing suite of neuroimaging software containers. Neurodesk includes a browser-accessible virtual desktop, command-line interface and computational notebook compatibility, allowing for accessible, flexible, portable and fully reproducible neuroimaging analysis on personal workstations, high-performance computers and the cloud.
Collapse
Affiliation(s)
- Angela I Renton
- The University of Queensland, Queensland Brain Institute, St Lucia, Brisbane, Queensland, Australia.
- The University of Queensland, School of Electrical Engineering and Computer Science, St Lucia, Brisbane, Queensland, Australia.
| | - Thuy T Dao
- The University of Queensland, School of Electrical Engineering and Computer Science, St Lucia, Brisbane, Queensland, Australia
| | - Tom Johnstone
- Centre for Mental Health & Brain Sciences, Swinburne University of Technology, Hawthorn, Victoria, Australia
| | - Oren Civier
- Centre for Mental Health & Brain Sciences, Swinburne University of Technology, Hawthorn, Victoria, Australia
| | - Ryan P Sullivan
- The University of Sydney, School of Biomedical Engineering, Sydney, New South Wales, Australia
| | - David J White
- Centre for Mental Health & Brain Sciences, Swinburne University of Technology, Hawthorn, Victoria, Australia
| | - Paris Lyons
- Centre for Mental Health & Brain Sciences, Swinburne University of Technology, Hawthorn, Victoria, Australia
| | - Benjamin M Slade
- Centre for Mental Health & Brain Sciences, Swinburne University of Technology, Hawthorn, Victoria, Australia
| | - David F Abbott
- The Florey Institute of Neuroscience and Mental Health, The University of Melbourne, Melbourne, Victoria, Australia
| | - Toluwani J Amos
- School of Life Science and Technology, University of Electronic Science and Technology, Chengdu, China
| | - Saskia Bollmann
- The University of Queensland, School of Electrical Engineering and Computer Science, St Lucia, Brisbane, Queensland, Australia
| | - Andy Botting
- Australian Research Data Commons (ARDC), Sydney, New South Wales, Australia
| | - Megan E J Campbell
- School of Psychological Sciences, University of Newcastle, Newcastle, New South Wales, Australia
- Hunter Medical Research Institute Imaging Centre, Newcastle, New South Wales, Australia
| | - Jeryn Chang
- The University of Queensland, School of Biomedical Sciences, St Lucia, Brisbane, Queensland, Australia
| | - Thomas G Close
- The University of Sydney, School of Biomedical Engineering, Sydney, New South Wales, Australia
| | - Monika Dörig
- Integrative Spinal Research, Department of Chiropractic Medicine, Balgrist University Hospital, University of Zurich, Zurich, Switzerland
| | - Korbinian Eckstein
- The University of Queensland, School of Electrical Engineering and Computer Science, St Lucia, Brisbane, Queensland, Australia
| | - Gary F Egan
- The Turner Institute for Brain and Mental Health, School of Psychological Sciences, Monash University, Melbourne, Victoria, Australia
- Monash Biomedical Imaging, Monash University, Melbourne, Victoria, Australia
| | - Stefanie Evas
- School of Psychology, University of Adelaide, Adelaide, South Australia, Australia
- Human Health, Health & Biosecurity, CSIRO, Adelaide, South Australia, Australia
| | - Guillaume Flandin
- Wellcome Centre for Human Neuroimaging, University College London, London, UK
| | - Kelly G Garner
- School of Psychology, University of New South Wales, Sydney, New South Wales, Australia
- The University of Queensland, School of Psychology, St Lucia, Brisbane, Queensland, Australia
| | - Marta I Garrido
- Melbourne School of Psychological Sciences, he University of Melbourne, Melbourne, Victoria, Australia
- Graeme Clark Institute for Biomedical Engineering, The University of Melbourne, Melbourne, Victoria, Australia
| | - Satrajit S Ghosh
- McGovern Institute for Brain Research, Massachusetts Institute of Technology, Cambridge, MA, USA
- Department of Otolaryngology - Head and Neck Surgery, Harvard Medical School, Boston, MA, USA
| | - Martin Grignard
- GIGA CRC In-Vivo Imaging, University of Liège, Liège, Belgium
| | - Yaroslav O Halchenko
- Center for Open Neuroscience, Department of Psychological and Brain Sciences, Dartmouth College, Hanover, NH, USA
| | - Anthony J Hannan
- The Florey Institute of Neuroscience and Mental Health, The University of Melbourne, Melbourne, Victoria, Australia
| | - Anibal S Heinsfeld
- Department of Psychology, Center for Perceptual Systems, Institute for Neuroscience, Center For Learning and Memory, The University of Texas at Austin, Austin, TX, USA
| | - Laurentius Huber
- National Institute of Mental Health (NIMH), National Institutes Health, Bethesda, MD, USA
| | - Matthew E Hughes
- Centre for Mental Health & Brain Sciences, Swinburne University of Technology, Hawthorn, Victoria, Australia
| | - Jakub R Kaczmarzyk
- Department of Biomedical Informatics, Stony Brook University, New York, NY, USA
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, New York, NY, USA
| | - Lars Kasper
- BRAIN-TO Lab, Krembil Brain Institute, University Health Network, Toronto, Ontario, Canada
- Translational Neuromodeling Unit, Institute for Biomedical Engineering, University of Zurich and ETH Zurich, Zurich, Switzerland
| | - Levin Kuhlmann
- Department of Data Science and AI, Faculty of Information Technology, Monash University, Melbourne, Victoria, Australia
| | - Kexin Lou
- The University of Queensland, School of Electrical Engineering and Computer Science, St Lucia, Brisbane, Queensland, Australia
- Department of Biomedical Engineering, Southern University of Science and Technology, Shenzhen, China
| | - Yorguin-Jose Mantilla-Ramos
- Grupo Neuropsicología y Conducta (GRUNECO), Facultad de Medicina, Universidad de Antioquia, Medellín, Colombia
| | - Jason B Mattingley
- The University of Queensland, Queensland Brain Institute, St Lucia, Brisbane, Queensland, Australia
- The University of Queensland, School of Psychology, St Lucia, Brisbane, Queensland, Australia
| | - Michael L Meier
- Integrative Spinal Research, Department of Chiropractic Medicine, Balgrist University Hospital, University of Zurich, Zurich, Switzerland
| | - Jo Morris
- Australian Research Data Commons (ARDC), Sydney, New South Wales, Australia
| | - Akshaiy Narayanan
- School of Computer Science, The University of Auckland, Auckland, New Zealand
| | - Franco Pestilli
- Department of Psychology, Center for Perceptual Systems, Institute for Neuroscience, Center For Learning and Memory, The University of Texas at Austin, Austin, TX, USA
| | - Aina Puce
- Department of Psychological and Brain Sciences, Indiana University, Bloomington, IN, USA
| | - Fernanda L Ribeiro
- The University of Queensland, School of Electrical Engineering and Computer Science, St Lucia, Brisbane, Queensland, Australia
| | - Nigel C Rogasch
- The Turner Institute for Brain and Mental Health, School of Psychological Sciences, Monash University, Melbourne, Victoria, Australia
- Discipline of Psychiatry, Adelaide Medical School, University of Adelaide, Adelaide, South Australia, Australia
- Hopwood Centre for Neurobiology, Lifelong Health Theme, South Australian Health and Medical Research Institute (SAHMRI), Adelaide, South Australia, Australia
| | - Chris Rorden
- McCausland Center for Brain Imaging, Department of Psychology, University of South Carolina, Columbia, SC, USA
| | - Mark M Schira
- School of Psychology, University of Wollongong, Wollongong, New South Wales, Australia
| | - Thomas B Shaw
- The University of Queensland, School of Electrical Engineering and Computer Science, St Lucia, Brisbane, Queensland, Australia
- The University of Queensland, Centre for Advanced Imaging, St Lucia, Brisbane, Queensland, Australia
- Department of Neurology, Royal Brisbane and Women's Hospital, Brisbane, Queensland, Australia
| | - Paul F Sowman
- Macquarie University, School of Psychological Sciences, Sydney, New South Wales, Australia
| | - Gershon Spitz
- Department of Neuroscience, Central Clinical School, Faculty of Medicine, Nursing and Health Sciences, Monash University, Melbourne, Victoria, Australia
- Monash-Epworth Rehabilitation Research Centre, Turner Institute for Brain and Mental Health, School of Psychological Sciences, Monash University, Melbourne, Victoria, Australia
| | - Ashley W Stewart
- The University of Queensland, School of Electrical Engineering and Computer Science, St Lucia, Brisbane, Queensland, Australia
- ARC Training Centre for Innovation in Biomedical Imaging Technology, The University of Queensland, Brisbane, Queensland, Australia
| | - Xincheng Ye
- The University of Queensland, School of Electrical Engineering and Computer Science, St Lucia, Brisbane, Queensland, Australia
| | - Judy D Zhu
- Macquarie University, School of Psychological Sciences, Sydney, New South Wales, Australia
| | - Aswin Narayanan
- The University of Queensland, Centre for Advanced Imaging, St Lucia, Brisbane, Queensland, Australia
| | - Steffen Bollmann
- The University of Queensland, School of Electrical Engineering and Computer Science, St Lucia, Brisbane, Queensland, Australia.
- The University of Queensland, Centre for Advanced Imaging, St Lucia, Brisbane, Queensland, Australia.
- ARC Training Centre for Innovation in Biomedical Imaging Technology, The University of Queensland, Brisbane, Queensland, Australia.
- Queensland Digital Health Centre, The University of Queensland, Brisbane, Queensland, Australia.
| |
Collapse
|
14
|
Roberts AG, Romano DJ, Şişman M, Dimov AV, Spincemaille P, Nguyen TD, Kovanlikaya I, Gauthier SA, Wang Y. Maximum spherical mean value filtering for whole-brain QSM. Magn Reson Med 2024; 91:1586-1597. [PMID: 38169132 PMCID: PMC11416845 DOI: 10.1002/mrm.29963] [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: 04/21/2023] [Revised: 10/30/2023] [Accepted: 11/19/2023] [Indexed: 01/05/2024]
Abstract
PURPOSE To develop a tissue field-filtering algorithm, called maximum spherical mean value (mSMV), for reducing shadow artifacts in QSM of the brain without requiring brain-tissue erosion. THEORY AND METHODS Residual background field is a major source of shadow artifacts in QSM. The mSMV algorithm filters large field-magnitude values near the border, where the maximum value of the harmonic background field is located. The effectiveness of mSMV for artifact removal was evaluated by comparing existing QSM algorithms in numerical brain simulation as well as using in vivo human data acquired from 11 healthy volunteers and 93 patients. RESULTS Numerical simulation showed that mSMV reduces shadow artifacts and improves QSM accuracy. Better shadow reduction, as demonstrated by lower QSM variation in the gray matter and higher QSM image quality score, was also observed in healthy subjects and in patients with hemorrhages, stroke, and multiple sclerosis. CONCLUSION The mSMV algorithm allows QSM maps that are substantially equivalent to those obtained using SMV-filtered dipole inversion without eroding the volume of interest.
Collapse
Affiliation(s)
- Alexandra G. Roberts
- Department of Electrical and Computer Engineering, Cornell University, Ithaca NY, USA
- Department of Radiology, Weill Cornell Medicine, New York, NY, USA
| | - Dominick J. Romano
- Department of Radiology, Weill Cornell Medicine, New York, NY, USA
- Meinig School of Biomedical Engineering, Cornell University, Ithaca NY, USA
| | - Mert Şişman
- Department of Electrical and Computer Engineering, Cornell University, Ithaca NY, USA
- Department of Radiology, Weill Cornell Medicine, New York, NY, USA
| | - Alexey V. Dimov
- Department of Radiology, Weill Cornell Medicine, New York, NY, USA
| | | | - Thanh D. Nguyen
- Department of Radiology, Weill Cornell Medicine, New York, NY, USA
| | | | | | - Yi Wang
- Department of Electrical and Computer Engineering, Cornell University, Ithaca NY, USA
- Department of Radiology, Weill Cornell Medicine, New York, NY, USA
- Meinig School of Biomedical Engineering, Cornell University, Ithaca NY, USA
| |
Collapse
|
15
|
Kurian D, Hagberg GE, Scheffler K, Paul JS. A predictor-corrector phase unwrapping algorithm for temporally undersampled gradient-echo MRI. Magn Reson Med 2024; 91:1707-1722. [PMID: 38084410 DOI: 10.1002/mrm.29964] [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/02/2023] [Revised: 11/18/2023] [Accepted: 11/19/2023] [Indexed: 02/03/2024]
Abstract
PURPOSE To develop a method for unwrapping temporally undersampled and nonlinear gradient recalled echo (GRE) phase. THEORY AND METHODS Temporal unwrapping is performed as a sequential one step prediction of the echo phase, followed by a correction to the nearest integer wrap-count. A spatio-temporal extension of the 1D predictor corrector unwrapping (PCU) algorithm improves the prediction accuracy, and thereby maintains spatial continuity. The proposed method is evaluated using numerical phantom, physical phantom, and in vivo brain data at both 3 T and 9.4 T. The unwrapping performance is compared with the state-of-the-art temporal and spatial unwrapping algorithms, and the spatio-temporal iterative virtual-echo based Nyquist sampled (iVENyS) algorithm. RESULTS Simulation results showed significant reduction in unwrapping errors at higher echoes compared with the state-of-the-art algorithms. Similar to the iVENyS algorithm, the PCU algorithm was able to generate spatially smooth phase images for in vivo data acquired at 3 T and 9.4 T, bypassing the use of additional spatial unwrapping step. A key advantage over iVENyS algorithm is the superior performance of PCU algorithm at higher echoes. CONCLUSION PCU algorithm serves as a robust phase unwrapping method for temporally undersampled and nonlinear GRE phase, particularly in the presence of high field gradients.
Collapse
Affiliation(s)
- Deepu Kurian
- School of Electronic Systems & Automation, Digital University Kerala, Trivandrum, Kerala, India
| | - Gisela E Hagberg
- High Field Magnetic Resonance, Max-Planck-Institute for Biological Cybernetics, Tübingen, Germany
- Biomedical Magnetic Resonance, Department of Radiology, Eberhard Karl's University and University Hospital, Tübingen, Germany
| | - Klaus Scheffler
- High Field Magnetic Resonance, Max-Planck-Institute for Biological Cybernetics, Tübingen, Germany
- Biomedical Magnetic Resonance, Department of Radiology, Eberhard Karl's University and University Hospital, Tübingen, Germany
| | - Joseph Suresh Paul
- School of Electronic Systems & Automation, Digital University Kerala, Trivandrum, Kerala, India
- School of Informatics, Digital University Kerala, Trivandrum, Kerala, India
| |
Collapse
|
16
|
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
|
17
|
Guan X, Lancione M, Ayton S, Dusek P, Langkammer C, Zhang M. Neuroimaging of Parkinson's disease by quantitative susceptibility mapping. Neuroimage 2024; 289:120547. [PMID: 38373677 DOI: 10.1016/j.neuroimage.2024.120547] [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: 07/30/2023] [Revised: 02/02/2024] [Accepted: 02/17/2024] [Indexed: 02/21/2024] Open
Abstract
Parkinson's disease (PD) is a common neurodegenerative disease, and apart from a few rare genetic causes, its pathogenesis remains largely unclear. Recent scientific interest has been captured by the involvement of iron biochemistry and the disruption of iron homeostasis, particularly within the brain regions specifically affected in PD. The advent of Quantitative Susceptibility Mapping (QSM) has enabled non-invasive quantification of brain iron in vivo by MRI, which has contributed to the understanding of iron-associated pathogenesis and has the potential for the development of iron-based biomarkers in PD. This review elucidates the biochemical underpinnings of brain iron accumulation, details advancements in iron-sensitive MRI technologies, and discusses the role of QSM as a biomarker of iron deposition in PD. Despite considerable progress, several challenges impede its clinical application after a decade of QSM studies. The initiation of multi-site research is warranted for developing robust, interpretable, and disease-specific biomarkers for monitoring PD disease progression.
Collapse
Affiliation(s)
- Xiaojun Guan
- Department of Radiology, the Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China; Joint Laboratory of Clinical Radiology, the Second Affiliated Hospital, Zhejiang University School of Medicine, No.88 Jiefang Road, Shangcheng District, Hangzhou 31009, China
| | - Marta Lancione
- Laboratory of Medical Physics and Magnetic Resonance, IRCCS Stella Maris Foundation, Pisa, Italy
| | - Scott Ayton
- Florey Institute, The University of Melbourne, Australia
| | - Petr Dusek
- Department of Neurology and Center of Clinical Neuroscience, First Faculty of Medicine, Charles University and General University Hospital in Prague, Prague, Czechia; Department of Radiology, First Faculty of Medicine, Charles University and General University Hospital in Prague, Auenbruggerplatz 22, Prague 8036, Czechia
| | | | - Minming Zhang
- Department of Radiology, the Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China; Joint Laboratory of Clinical Radiology, the Second Affiliated Hospital, Zhejiang University School of Medicine, No.88 Jiefang Road, Shangcheng District, Hangzhou 31009, China.
| |
Collapse
|
18
|
Thomas GE, Hannaway N, Zarkali A, Shmueli K, Weil RS. Longitudinal Associations of Magnetic Susceptibility with Clinical Severity in Parkinson's Disease. Mov Disord 2024; 39:546-559. [PMID: 38173297 PMCID: PMC11141787 DOI: 10.1002/mds.29702] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Revised: 11/29/2023] [Accepted: 12/11/2023] [Indexed: 01/05/2024] Open
Abstract
BACKGROUND Dementia is common in Parkinson's disease (PD), but there is wide variation in its timing. A critical gap in PD research is the lack of quantifiable markers of progression, and methods to identify early stages of dementia. Atrophy-based magnetic resonance imaging (MRI) has limited sensitivity in detecting or tracking changes relating to PD dementia, but quantitative susceptibility mapping (QSM), sensitive to brain tissue iron, shows potential for these purposes. OBJECTIVE The objective of the paper is to study, for the first time, the longitudinal relationship between cognition and QSM in PD in detail. METHODS We present a longitudinal study of clinical severity in PD using QSM, including 59 PD patients (without dementia at study onset), and 22 controls over 3 years. RESULTS In PD, increased baseline susceptibility in the right temporal cortex, nucleus basalis of Meynert, and putamen was associated with greater cognitive severity after 3 years; and increased baseline susceptibility in basal ganglia, substantia nigra, red nucleus, insular cortex, and dentate nucleus was associated with greater motor severity after 3 years. Increased follow-up susceptibility in these regions was associated with increased follow-up cognitive and motor severity, with further involvement of hippocampus relating to cognitive severity. However, there were no consistent increases in susceptibility over 3 years. CONCLUSIONS Our study suggests that QSM may predict changes in cognitive severity many months prior to overt cognitive involvement in PD. However, we did not find robust longitudinal changes in QSM over the course of the study. Additional tissue metrics may be required together with QSM for it to monitor progression in clinical practice and therapeutic trials. © 2024 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.
Collapse
Affiliation(s)
| | - Naomi Hannaway
- Dementia Research CentreUCL Institute of NeurologyLondonUK
| | | | - Karin Shmueli
- Department of Medical Physics and Biomedical EngineeringUniversity College LondonLondonUK
| | - Rimona S. Weil
- Dementia Research CentreUCL Institute of NeurologyLondonUK
- Wellcome Centre for Human NeuroimagingUniversity College LondonLondonUK
- Movement Disorders ConsortiumUniversity College LondonLondonUK
| |
Collapse
|
19
|
Bachrata B, Bollmann S, Jin J, Tourell M, Dal-Bianco A, Trattnig S, Barth M, Ropele S, Enzinger C, Robinson SD. Super-resolution QSM in little or no additional time for imaging (NATIve) using 2D EPI imaging in 3 orthogonal planes. Neuroimage 2023; 283:120419. [PMID: 37871759 DOI: 10.1016/j.neuroimage.2023.120419] [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: 05/16/2023] [Revised: 09/22/2023] [Accepted: 10/20/2023] [Indexed: 10/25/2023] Open
Abstract
Quantitative Susceptibility Mapping has the potential to provide additional insights into neurological diseases but is typically based on a quite long (5-10 min) 3D gradient-echo scan which is highly sensitive to motion. We propose an ultra-fast acquisition based on three orthogonal (sagittal, coronal and axial) 2D simultaneous multi-slice EPI scans with 1 mm in-plane resolution and 3 mm thick slices. Images in each orientation are corrected for susceptibility-related distortions and co-registered with an iterative non-linear Minimum Deformation Averaging (Volgenmodel) approach to generate a high SNR, super-resolution data set with an isotropic resolution of close to 1 mm. The net acquisition time is 3 times the volume acquisition time of EPI or about 12 s, but the three volumes could also replace "dummy scans" in fMRI, making it feasible to acquire QSM in little or No Additional Time for Imaging (NATIve). NATIve QSM values agreed well with reference 3D GRE QSM in the basal ganglia in healthy subjects. In patients with multiple sclerosis, there was also a good agreement between the susceptibility values within lesions and control ROIs and all lesions which could be seen on 3D GRE QSMs could also be visualized on NATIve QSMs. The approach is faster than conventional 3D GRE by a factor of 25-50 and faster than 3D EPI by a factor of 3-5. As a 2D technique, NATIve QSM was shown to be much more robust to motion than the 3D GRE and 3D EPI, opening up the possibility of studying neurological diseases involving iron accumulation and demyelination in patients who find it difficult to lie still for long enough to acquire QSM data with conventional methods.
Collapse
Affiliation(s)
- Beata Bachrata
- High Field MR Centre, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Austria; Karl Landsteiner Institute for Clinical Molecular MR in Musculoskeletal Imaging, Vienna, Austria; Department of Medical Engineering, Carinthia University of Applied Sciences, Klagenfurt, Austria
| | - Steffen Bollmann
- Centre of Advanced Imaging, University of Queensland, Brisbane, Australia; ARC Training Centre for Innovation in Biomedical Imaging Technology, University of Queensland, Brisbane, Australia; School of Information Technology and Electrical Engineering, The University of Queensland, Brisbane, Australia
| | - Jin Jin
- Centre of Advanced Imaging, University of Queensland, Brisbane, Australia; Siemens Healthcare Pty Ltd, Australia
| | - Monique Tourell
- Centre of Advanced Imaging, University of Queensland, Brisbane, Australia; ARC Training Centre for Innovation in Biomedical Imaging Technology, University of Queensland, Brisbane, Australia
| | - Assunta Dal-Bianco
- Department of Neurology, Medical University of Vienna, Austria; Comprehensive Center for Clinical Neurosciences and Mental Health, Medical University of Vienna, Austria
| | - Siegfried Trattnig
- High Field MR Centre, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Austria; Karl Landsteiner Institute for Clinical Molecular MR in Musculoskeletal Imaging, Vienna, Austria
| | - Markus Barth
- Centre of Advanced Imaging, University of Queensland, Brisbane, Australia; ARC Training Centre for Innovation in Biomedical Imaging Technology, University of Queensland, Brisbane, Australia; School of Information Technology and Electrical Engineering, The University of Queensland, Brisbane, Australia
| | - Stefan Ropele
- Department of Neurology, Medical University of Graz, Austria
| | | | - Simon Daniel Robinson
- High Field MR Centre, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Austria; Centre of Advanced Imaging, University of Queensland, Brisbane, Australia; Department of Neurology, Medical University of Graz, Austria.
| |
Collapse
|
20
|
Naji N, Wilman A. Thin slab quantitative susceptibility mapping. Magn Reson Med 2023; 90:2290-2305. [PMID: 37526029 DOI: 10.1002/mrm.29800] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Revised: 06/27/2023] [Accepted: 06/30/2023] [Indexed: 08/02/2023]
Abstract
PURPOSE Susceptibility maps reconstructed from thin slabs may suffer underestimation due to background-field removal imperfections near slab boundaries and the increased difficulty of solving a 3D-inversion problem with reduced support, particularly in the direction of the main magnetic field. Reliable QSM reconstruction from thin slabs would enable focal acquisitions in a much-reduced scan time. METHODS This work proposes using additional rapid low-resolution data of extended spatial coverage to improve background-field estimation and regularize the inversion-to-susceptibility process for high resolution, thin slab data. The new method was tested using simulated and in-vivo brain data of high resolution (0.33 × 0.33 × 0.33 mm3 and 0.54 × 0.54 × 0.65 mm3 , respectively) at 3T, and compared to the standard large volume approach. RESULTS Using the proposed method, in-vivo high-resolution QSM at 3T was obtained from slabs of width as small as 10.4 mm, aided by a lower-resolution dataset of 24 times coarser voxels. Simulations showed that the proposed method produced more consistent measurements from slabs of at least eight slices. Reducing the mean ROI error to 5% required the low-resolution data to cover ˜60 mm in the direction of the main field, have at least 2-mm isotropic resolution that is not coarser than the high-resolution data by more than four-fold in any direction. CONCLUSION Applying the proposed method enabled focal QSM acquisitions at sub-millimeter resolution within reasonable acquisition time.
Collapse
Affiliation(s)
- Nashwan Naji
- Department of Biomedical Engineering, University of Alberta, Edmonton, Alberta, Canada
| | - Alan Wilman
- Department of Biomedical Engineering, University of Alberta, Edmonton, Alberta, Canada
| |
Collapse
|
21
|
Van AN, Montez DF, Laumann TO, Suljic V, Madison T, Baden NJ, Ramirez-Perez N, Scheidter KM, Monk JS, Whiting FI, Adeyemo B, Chauvin RJ, Krimmel SR, Metoki A, Rajesh A, Roland JL, Salo T, Wang A, Weldon KB, Sotiras A, Shimony JS, Kay BP, Nelson SM, Tervo-Clemmens B, Marek SA, Vizioli L, Yacoub E, Satterthwaite TD, Gordon EM, Fair DA, Tisdall MD, Dosenbach NU. Framewise multi-echo distortion correction for superior functional MRI. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.11.28.568744. [PMID: 38077010 PMCID: PMC10705259 DOI: 10.1101/2023.11.28.568744] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/19/2023]
Abstract
Functional MRI (fMRI) data are severely distorted by magnetic field (B0) inhomogeneities which currently must be corrected using separately acquired field map data. However, changes in the head position of a scanning participant across fMRI frames can cause changes in the B0 field, preventing accurate correction of geometric distortions. Additionally, field maps can be corrupted by movement during their acquisition, preventing distortion correction altogether. In this study, we use phase information from multi-echo (ME) fMRI data to dynamically sample distortion due to fluctuating B0 field inhomogeneity across frames by acquiring multiple echoes during a single EPI readout. Our distortion correction approach, MEDIC (Multi-Echo DIstortion Correction), accurately estimates B0 related distortions for each frame of multi-echo fMRI data. Here, we demonstrate that MEDIC's framewise distortion correction produces improved alignment to anatomy and decreases the impact of head motion on resting-state functional connectivity (RSFC) maps, in higher motion data, when compared to the prior gold standard approach (i.e., TOPUP). Enhanced framewise distortion correction with MEDIC, without the requirement for field map collection, furthers the advantage of multi-echo over single-echo fMRI.
Collapse
Affiliation(s)
- Andrew N. Van
- Department of Biomedical Engineering, Washington University in St. Louis, MO 63130
- Department of Neurology, Washington University School of Medicine, St. Louis, MO 63110
| | - David F. Montez
- Department of Neurology, Washington University School of Medicine, St. Louis, MO 63110
- Department of Psychiatry, Washington University School of Medicine, St. Louis, MO 63110
| | - Timothy O. Laumann
- Department of Psychiatry, Washington University School of Medicine, St. Louis, MO 63110
| | - Vahdeta Suljic
- Department of Neurology, Washington University School of Medicine, St. Louis, MO 63110
| | - Thomas Madison
- Institute of Child Development, University of Minnesota Medical School, Minneapolis, MN 55455
- Masonic Institute for the Developing Brain, University of Minnesota Medical School, Minneapolis, MN 55455
| | - Noah J. Baden
- Department of Neurology, Washington University School of Medicine, St. Louis, MO 63110
| | | | - Kristen M. Scheidter
- Department of Neurology, Washington University School of Medicine, St. Louis, MO 63110
| | - Julia S. Monk
- Department of Neurology, Washington University School of Medicine, St. Louis, MO 63110
| | - Forrest I. Whiting
- Department of Neurology, Washington University School of Medicine, St. Louis, MO 63110
| | - Babatunde Adeyemo
- Department of Neurology, Washington University School of Medicine, St. Louis, MO 63110
| | - Roselyne J. Chauvin
- Department of Neurology, Washington University School of Medicine, St. Louis, MO 63110
| | - Samuel R. Krimmel
- Department of Neurology, Washington University School of Medicine, St. Louis, MO 63110
| | - Athanasia Metoki
- Department of Neurology, Washington University School of Medicine, St. Louis, MO 63110
| | - Aishwarya Rajesh
- Department of Radiology, Washington University School of Medicine, St. Louis, MO 63110
| | - Jarod L. Roland
- Department of Neurosurgery, Washington University School of Medicine, St. Louis, MO 63110
| | - Taylor Salo
- Lifespan Informatics and Neuroimaging Center (PennLINC), Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, 19104
| | - Anxu Wang
- Department of Neurology, Washington University School of Medicine, St. Louis, MO 63110
- Division of Computation and Data Science, Washington University School of Medicine, St. Louis, MO 63110
| | - Kimberly B. Weldon
- Masonic Institute for the Developing Brain, University of Minnesota Medical School, Minneapolis, MN 55455
| | - Aristeidis Sotiras
- Department of Radiology, Washington University School of Medicine, St. Louis, MO 63110
- Institute for Informatics, Data Science & Biostatistics, Washington University School of Medicine, St. Louis, MO 63130
| | - Joshua S. Shimony
- Department of Radiology, Washington University School of Medicine, St. Louis, MO 63110
| | - Benjamin P. Kay
- Department of Neurology, Washington University School of Medicine, St. Louis, MO 63110
| | - Steven M. Nelson
- Masonic Institute for the Developing Brain, University of Minnesota Medical School, Minneapolis, MN 55455
- Department of Pediatrics, University of Minnesota Medical School, Minneapolis, MN 55455
| | - Brenden Tervo-Clemmens
- Masonic Institute for the Developing Brain, University of Minnesota Medical School, Minneapolis, MN 55455
- Department of Psychiatry & Behavioral Sciences, University of Minnesota Medical School, Minneapolis, MN 55455
| | - Scott A. Marek
- Department of Radiology, Washington University School of Medicine, St. Louis, MO 63110
| | - Luca Vizioli
- Center for Magnetic Resonance Research, University of Minnesota Medical School, Minneapolis, MN 55455
| | - Essa Yacoub
- Center for Magnetic Resonance Research, University of Minnesota Medical School, Minneapolis, MN 55455
| | - Theodore D. Satterthwaite
- Lifespan Informatics and Neuroimaging Center (PennLINC), Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, 19104
| | - Evan M. Gordon
- Department of Radiology, Washington University School of Medicine, St. Louis, MO 63110
| | - Damien A. Fair
- Institute of Child Development, University of Minnesota Medical School, Minneapolis, MN 55455
- Masonic Institute for the Developing Brain, University of Minnesota Medical School, Minneapolis, MN 55455
- Department of Pediatrics, University of Minnesota Medical School, Minneapolis, MN 55455
| | - M. Dylan Tisdall
- Department of Radiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104
| | - Nico U.F. Dosenbach
- Department of Biomedical Engineering, Washington University in St. Louis, MO 63130
- Department of Neurology, Washington University School of Medicine, St. Louis, MO 63110
- Department of Radiology, Washington University School of Medicine, St. Louis, MO 63110
- Department of Pediatrics, Washington University School of Medicine, St. Louis, MO 63110
| |
Collapse
|
22
|
Cheng J, Song M, Xu Z, Zheng Q, Zhu L, Chen W, Feng Y, Bao J, Cheng J. A new 3D phase unwrapping method by region partitioning and local polynomial modeling in abdominal quantitative susceptibility mapping. Front Neurosci 2023; 17:1287788. [PMID: 38033538 PMCID: PMC10684715 DOI: 10.3389/fnins.2023.1287788] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2023] [Accepted: 10/30/2023] [Indexed: 12/02/2023] Open
Abstract
Background Accurate phase unwrapping is a critical prerequisite for successful applications in phase-related MRI, including quantitative susceptibility mapping (QSM) and susceptibility weighted imaging. However, many existing 3D phase unwrapping algorithms face challenges in the presence of severe noise, rapidly changing phase, and open-end cutline. Methods In this study, we introduce a novel 3D phase unwrapping approach utilizing region partitioning and a local polynomial model. Initially, the method leverages phase partitioning to create initial regions. Noisy voxels connecting areas within these regions are excluded and grouped into residual voxels. The connected regions within the region of interest are then reidentified and categorized into blocks and residual voxels based on voxel count thresholds. Subsequently, the method sequentially performs inter-block and residual voxel phase unwrapping using the local polynomial model. The proposed method was evaluated on simulation and in vivo abdominal QSM data, and was compared with the classical Region-growing, Laplacian_based, Graph-cut, and PRELUDE methods. Results Simulation experiments, conducted under different signal-to-noise ratios and phase change levels, consistently demonstrate that the proposed method achieves accurate unwrapping results, with mean error ratios not exceeding 0.01%. In contrast, the error ratios of Region-growing (N/A, 84.47%), Laplacian_based (20.65%, N/A), Graph-cut (2.26%, 20.71%), and PRELUDE (4.28%, 10.33%) methods are all substantially higher than those of the proposed method. In vivo abdominal QSM experiments further confirm the effectiveness of the proposed method in unwrapping phase data and successfully reconstructing susceptibility maps, even in scenarios with significant noise, rapidly changing phase, and open-end cutline in a large field of view. Conclusion The proposed method demonstrates robust and accurate phase unwrapping capabilities, positioning it as a promising option for abdominal QSM applications.
Collapse
Affiliation(s)
- Junying Cheng
- Department of Magnetic Resonance Imaging, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
| | - Manli Song
- Department of Magnetic Resonance Imaging, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
| | - Zhongbiao Xu
- Department of Radiotherapy, Cancer Center, Guangdong Provincial People’s Hospital, Guangdong Academy of Medical Science, Guangzhou, China
| | - Qian Zheng
- College of Software Engineering, Zhengzhou University of Light Industry, Zhengzhou, China
| | - Li Zhu
- Department of Magnetic Resonance Imaging, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
| | - Wufan Chen
- School of Biomedical Engineering and Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, Guangzhou, China
| | - Yanqiu Feng
- School of Biomedical Engineering and Guangdong Provincial Key Laboratory of Medical Image Processing, Southern Medical University, Guangzhou, China
| | - Jianfeng Bao
- Department of Magnetic Resonance Imaging, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
| | - Jingliang Cheng
- Department of Magnetic Resonance Imaging, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
| |
Collapse
|
23
|
Robinson SD, Bachrata B, Eckstein K, Bollmann S, Bollmann S, Hodono S, Cloos M, Tourell M, Jin J, O'Brien K, Reutens DC, Trattnig S, Enzinger C, Barth M. Improved dynamic distortion correction for fMRI using single-echo EPI and a readout-reversed first image (REFILL). Hum Brain Mapp 2023; 44:5095-5112. [PMID: 37548414 PMCID: PMC10502646 DOI: 10.1002/hbm.26440] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2023] [Revised: 06/01/2023] [Accepted: 07/12/2023] [Indexed: 08/08/2023] Open
Abstract
The boundaries between tissues with different magnetic susceptibilities generate inhomogeneities in the main magnetic field which change over time due to motion, respiration and system instabilities. The dynamically changing field can be measured from the phase of the fMRI data and corrected. However, methods for doing so need multi-echo data, time-consuming reference scans and/or involve error-prone processing steps, such as phase unwrapping, which are difficult to implement robustly on the MRI host. The improved dynamic distortion correction method we propose is based on the phase of the single-echo EPI data acquired for fMRI, phase offsets calculated from a triple-echo, bipolar reference scan of circa 3-10 s duration using a method which avoids the need for phase unwrapping and an additional correction derived from one EPI volume in which the readout direction is reversed. This Reverse-Encoded First Image and Low resoLution reference scan (REFILL) approach is shown to accurately measure B0 as it changes due to shim, motion and respiration, even with large dynamic changes to the field at 7 T, where it led to a > 20% increase in time-series signal to noise ratio compared to data corrected with the classic static approach. fMRI results from REFILL-corrected data were free of stimulus-correlated distortion artefacts seen when data were corrected with static field mapping. The method is insensitive to shim changes and eddy current differences between the reference scan and the fMRI time series, and employs calculation steps that are simple and robust, allowing most data processing to be performed in real time on the scanner image reconstruction computer. These improvements make it feasible to routinely perform dynamic distortion correction in fMRI.
Collapse
Affiliation(s)
- Simon Daniel Robinson
- Centre of Advanced ImagingUniversity of QueenslandBrisbaneAustralia
- Department of NeurologyMedical University of GrazGrazAustria
- High Field MR Centre, Department of Biomedical Imaging and Image‐Guided TherapyMedical University of ViennaViennaAustria
- Karl Landsteiner Institute for Clinical Molecular MR in Musculoskeletal ImagingViennaAustria
| | - Beata Bachrata
- High Field MR Centre, Department of Biomedical Imaging and Image‐Guided TherapyMedical University of ViennaViennaAustria
- Karl Landsteiner Institute for Clinical Molecular MR in Musculoskeletal ImagingViennaAustria
- Department of Medical EngineeringCarinthia University of Applied SciencesKlagenfurtAustria
| | - Korbinian Eckstein
- High Field MR Centre, Department of Biomedical Imaging and Image‐Guided TherapyMedical University of ViennaViennaAustria
| | - Saskia Bollmann
- Centre of Advanced ImagingUniversity of QueenslandBrisbaneAustralia
| | - Steffen Bollmann
- School of Information Technology and Electrical EngineeringThe University of QueenslandBrisbaneAustralia
| | - Shota Hodono
- Centre of Advanced ImagingUniversity of QueenslandBrisbaneAustralia
- ARC Training Centre for Innovation in Biomedical Imaging Technology (CIBIT)The University of QueenslandBrisbaneAustralia
| | - Martijn Cloos
- Centre of Advanced ImagingUniversity of QueenslandBrisbaneAustralia
- ARC Training Centre for Innovation in Biomedical Imaging Technology (CIBIT)The University of QueenslandBrisbaneAustralia
| | - Monique Tourell
- Centre of Advanced ImagingUniversity of QueenslandBrisbaneAustralia
- Siemens Healthcare Pty Ltd.BrisbaneAustralia
| | - Jin Jin
- Siemens Healthcare Pty Ltd.BrisbaneAustralia
| | | | - David C. Reutens
- Centre of Advanced ImagingUniversity of QueenslandBrisbaneAustralia
- ARC Training Centre for Innovation in Biomedical Imaging Technology (CIBIT)The University of QueenslandBrisbaneAustralia
| | - Siegfried Trattnig
- High Field MR Centre, Department of Biomedical Imaging and Image‐Guided TherapyMedical University of ViennaViennaAustria
| | | | - Markus Barth
- Centre of Advanced ImagingUniversity of QueenslandBrisbaneAustralia
- School of Information Technology and Electrical EngineeringThe University of QueenslandBrisbaneAustralia
| |
Collapse
|
24
|
Lee JY, Mack AF, Mattheus U, Donato S, Longo R, Tromba G, Shiozawa T, Scheffler K, Hagberg GE. Distribution of corpora amylacea in the human midbrain: using synchrotron radiation phase-contrast microtomography, high-field magnetic resonance imaging, and histology. Front Neurosci 2023; 17:1236876. [PMID: 37869518 PMCID: PMC10586329 DOI: 10.3389/fnins.2023.1236876] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2023] [Accepted: 09/06/2023] [Indexed: 10/24/2023] Open
Abstract
Corpora amylacea (CA) are polyglucosan aggregated granules that accumulate in the human body throughout aging. In the cerebrum, CA have been found in proximity to ventricular walls, pial surfaces, and blood vessels. However, studies showing their three-dimensional spatial distribution are sparse. In this study, volumetric images of four human brain stems were obtained with MRI and phase-contrast X-ray microtomography, followed up by Periodic acid Schiff stain for validation. CA appeared as hyperintense spheroid structures with diameters up to 30 μm. An automatic pipeline was developed to segment the CA, and the spatial distribution of over 200,000 individual corpora amylacea could be investigated. A threefold-or higher-density of CA was detected in the dorsomedial column of the periaqueductal gray (860-4,200 CA count/mm3) than in the superior colliculus (150-340 CA count/mm3). We estimated that about 2% of the CA were located in the immediate vicinity of the vessels or in the peri-vascular space. While CA in the ependymal lining of the cerebral aqueduct was rare, the sub-pial tissue of the anterior and posterior midbrain contained several CA. In the sample with the highest CA density, quantitative maps obtained with MRI revealed high R2* values and a diamagnetic shift in a region which spatially coincided with the CA dense region.
Collapse
Affiliation(s)
- Ju Young Lee
- Graduate Training Centre of Neuroscience, Eberhard Karl's University of Tübingen, Tübingen, Germany
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Andreas F. Mack
- Institute of Clinical Anatomy and Cell Analysis, Eberhard Karl's University of Tübingen, Tübingen, Germany
| | - Ulrich Mattheus
- Institute of Clinical Anatomy and Cell Analysis, Eberhard Karl's University of Tübingen, Tübingen, Germany
| | - Sandro Donato
- Department of Physics and STAR-LAB, University of Calabria, Rende, Italy
- Division of Frascati, Istituto Nazionale di Fisica Nucleare (INFN), Frascati, Italy
| | - Renata Longo
- Department of Physics, University of Trieste, Trieste, Italy
- Division of Trieste, Istituto Nazionale di Fisica Nucleare (INFN), Trieste, Italy
| | | | - Thomas Shiozawa
- Institute of Clinical Anatomy and Cell Analysis, Eberhard Karl's University of Tübingen, Tübingen, Germany
| | - Klaus Scheffler
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
- Department of Biomedical Magnetic Resonance, University Hospital Tübingen, Tübingen, Germany
| | - Gisela E. Hagberg
- High Field Magnetic Resonance, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
- Department of Biomedical Magnetic Resonance, University Hospital Tübingen, Tübingen, Germany
| |
Collapse
|
25
|
Velikina JV, Zhao R, Buelo CJ, Samsonov AA, Reeder SB, Hernando D. Data adaptive regularization with reference tissue constraints for liver quantitative susceptibility mapping. Magn Reson Med 2023; 90:385-399. [PMID: 36929781 PMCID: PMC11057046 DOI: 10.1002/mrm.29644] [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/22/2022] [Revised: 02/24/2023] [Accepted: 03/05/2023] [Indexed: 03/18/2023]
Abstract
PURPOSE To improve repeatability and reproducibility across acquisition parameters and reduce bias in quantitative susceptibility mapping (QSM) of the liver, through development of an optimized regularized reconstruction algorithm for abdominal QSM. METHODS An optimized approach to estimation of magnetic susceptibility distribution is formulated as a constrained reconstruction problem that incorporates estimates of the input data reliability and anatomical priors available from chemical shift-encoded imaging. The proposed data-adaptive method was evaluated with respect to bias, repeatability, and reproducibility in a patient population with a wide range of liver iron concentration (LIC). The proposed method was compared to the previously proposed and validated approach in liver QSM for two multi-echo spoiled gradient-recalled echo protocols with different acquisition parameters at 3T. Linear regression was used for evaluation of QSM methods against a reference FDA-approvedR 2 $$ {R}_2 $$ -based LIC measure andR 2 ∗ $$ {R}_2^{\ast } $$ measurements; repeatability/reproducibility were assessed by Bland-Altman analysis. RESULTS The data-adaptive method produced susceptibility maps with higher subjective quality due to reduced shading artifacts. For both acquisition protocols, higher linear correlation with bothR 2 $$ {R}_2 $$ - andR 2 ∗ $$ {R}_2^{\ast } $$ -based measurements were observed for the data-adaptive method (r 2 = 0 . 74 / 0 . 69 $$ {r}^2=0.74/0.69 $$ forR 2 $$ {R}_2 $$ ,0 . 97 / 0 . 95 $$ 0.97/0.95 $$ forR 2 ∗ $$ {R}_2^{\ast } $$ ) than the standard method (r 2 = 0 . 60 / 0 . 66 $$ {r}^2=0.60/0.66 $$ and0 . 79 / 0 . 88 $$ 0.79/0.88 $$ ). For both protocols, the data-adaptive method enabled better test-retest repeatability (repeatability coefficients 0.19/0.30 ppm for the data-adaptive method, 0.38/0.47 ppm for the standard method) and reproducibility across protocols (reproducibility coefficient 0.28 vs. 0.53ppm) than the standard method. CONCLUSIONS The proposed data-adaptive QSM algorithm may enable quantification of LIC with improved repeatability/reproducibility across different acquisition parameters as 3T.
Collapse
Affiliation(s)
- Julia V Velikina
- Department of Radiology, University of Wisconsin, Madison, Wisconsin, USA
| | - Ruiyang Zhao
- Department of Radiology, University of Wisconsin, Madison, Wisconsin, USA
- Department of Medical Physics, University of Wisconsin, Madison, Wisconsin, USA
| | - Collin J Buelo
- Department of Radiology, University of Wisconsin, Madison, Wisconsin, USA
- Department of Medical Physics, University of Wisconsin, Madison, Wisconsin, USA
| | - Alexey A Samsonov
- Department of Radiology, University of Wisconsin, Madison, Wisconsin, USA
| | - Scott B Reeder
- Department of Radiology, University of Wisconsin, Madison, Wisconsin, USA
- Department of Medical Physics, University of Wisconsin, Madison, Wisconsin, USA
- Department of Biomedical Engineering, University of Wisconsin, Madison, Wisconsin, USA
- Department of Medicine, University of Wisconsin, Madison, WI, USA
- Department of Emergency Medicine, University of Wisconsin, Madison, Wisconsin, USA
| | - Diego Hernando
- Department of Radiology, University of Wisconsin, Madison, Wisconsin, USA
- Department of Medical Physics, University of Wisconsin, Madison, Wisconsin, USA
| |
Collapse
|
26
|
Bissell MM, Raimondi F, Ait Ali L, Allen BD, Barker AJ, Bolger A, Burris N, Carhäll CJ, Collins JD, Ebbers T, Francois CJ, Frydrychowicz A, Garg P, Geiger J, Ha H, Hennemuth A, Hope MD, Hsiao A, Johnson K, Kozerke S, Ma LE, Markl M, Martins D, Messina M, Oechtering TH, van Ooij P, Rigsby C, Rodriguez-Palomares J, Roest AAW, Roldán-Alzate A, Schnell S, Sotelo J, Stuber M, Syed AB, Töger J, van der Geest R, Westenberg J, Zhong L, Zhong Y, Wieben O, Dyverfeldt P. 4D Flow cardiovascular magnetic resonance consensus statement: 2023 update. J Cardiovasc Magn Reson 2023; 25:40. [PMID: 37474977 PMCID: PMC10357639 DOI: 10.1186/s12968-023-00942-z] [Citation(s) in RCA: 42] [Impact Index Per Article: 42.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2023] [Accepted: 05/30/2023] [Indexed: 07/22/2023] Open
Abstract
Hemodynamic assessment is an integral part of the diagnosis and management of cardiovascular disease. Four-dimensional cardiovascular magnetic resonance flow imaging (4D Flow CMR) allows comprehensive and accurate assessment of flow in a single acquisition. This consensus paper is an update from the 2015 '4D Flow CMR Consensus Statement'. We elaborate on 4D Flow CMR sequence options and imaging considerations. The document aims to assist centers starting out with 4D Flow CMR of the heart and great vessels with advice on acquisition parameters, post-processing workflows and integration into clinical practice. Furthermore, we define minimum quality assurance and validation standards for clinical centers. We also address the challenges faced in quality assurance and validation in the research setting. We also include a checklist for recommended publication standards, specifically for 4D Flow CMR. Finally, we discuss the current limitations and the future of 4D Flow CMR. This updated consensus paper will further facilitate widespread adoption of 4D Flow CMR in the clinical workflow across the globe and aid consistently high-quality publication standards.
Collapse
Affiliation(s)
- Malenka M Bissell
- Department of Biomedical Imaging Science, Leeds Institute of Cardiovascular and Metabolic Medicine (LICAMM), LIGHT Laboratories, Clarendon Way, University of Leeds, Leeds, LS2 9NL, UK.
| | | | - Lamia Ait Ali
- Institute of Clinical Physiology CNR, Massa, Italy
- Foundation CNR Tuscany Region G. Monasterio, Massa, Italy
| | - Bradley D Allen
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, IL, USA
| | - Alex J Barker
- Department of Radiology, Children's Hospital Colorado, University of Colorado Anschutz Medical Center, Aurora, USA
| | - Ann Bolger
- Department of Medicine, University of California, San Francisco, CA, USA
- Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
| | - Nicholas Burris
- Department of Radiology, University of Michigan, Ann Arbor, USA
| | - Carl-Johan Carhäll
- Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
| | | | - Tino Ebbers
- Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
| | | | - Alex Frydrychowicz
- Department of Radiology and Nuclear Medicine, University Hospital Schleswig-Holstein, Campus Lübeck and Universität Zu Lübeck, Lübeck, Germany
| | - Pankaj Garg
- Norwich Medical School, University of East Anglia, Norwich, UK
| | - Julia Geiger
- Department of Diagnostic Imaging, University Children's Hospital, Zurich, Switzerland
- Children's Research Center, University Children's Hospital Zurich, Zurich, Switzerland
| | - Hojin Ha
- Department of Mechanical and Biomedical Engineering, Kangwon National University, Chuncheon, South Korea
| | - Anja Hennemuth
- Institute of Computer-Assisted Cardiovascular Medicine, Charité - Universitätsmedizin, Berlin, Germany
- German Center for Cardiovascular Research (DZHK), Partner Site, Berlin, Germany
- Department of Diagnostic and Interventional Radiology and Nuclear Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Michael D Hope
- Department of Radiology and Biomedical Imaging, University of California, San Francisco, CA, USA
| | - Albert Hsiao
- Department of Radiology, University of California, San Diego, CA, USA
| | - Kevin Johnson
- Departments of Radiology and Medical Physics, University of Wisconsin, Madison, WI, USA
| | - Sebastian Kozerke
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Liliana E Ma
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, IL, USA
| | - Michael Markl
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, IL, USA
| | - Duarte Martins
- Department of Pediatric Cardiology, Hospital de Santa Cruz, Centro Hospitalar Lisboa Ocidental, Lisbon, Portugal
| | - Marci Messina
- Department of Radiology, Northwestern Medicine, Chicago, IL, USA
| | - Thekla H Oechtering
- Department of Radiology and Nuclear Medicine, University Hospital Schleswig-Holstein, Campus Lübeck and Universität Zu Lübeck, Lübeck, Germany
- Departments of Radiology and Medical Physics, University of Wisconsin, Madison, WI, USA
| | - Pim van Ooij
- Department of Radiology & Nuclear Medicine, Amsterdam Cardiovascular Sciences, Amsterdam Movement Sciences, Amsterdam University Medical Centers, Location AMC, Amsterdam, The Netherlands
- Department of Pediatric Cardiology, Division of Pediatrics, Wilhelmina Children's Hospital, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Cynthia Rigsby
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, IL, USA
- Department of Medical Imaging, Ann & Robert H Lurie Children's Hospital of Chicago, Chicago, IL, USA
| | - Jose Rodriguez-Palomares
- Department of Cardiology, Hospital Universitari Vall d´Hebron,Vall d'Hebron Institut de Recerca (VHIR), Universitat Autònoma de Barcelona, Barcelona, Spain
- Centro de Investigación Biomédica en Red-CV, CIBER CV, Madrid, Spain
| | - Arno A W Roest
- Department of Pediatric Cardiology, Willem-Alexander's Children Hospital, Leiden University Medical Center and Center for Congenital Heart Defects Amsterdam-Leiden, Leiden, The Netherlands
| | | | - Susanne Schnell
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, IL, USA
- Department of Medical Physics, Institute of Physics, University of Greifswald, Greifswald, Germany
| | - Julio Sotelo
- School of Biomedical Engineering, Universidad de Valparaíso, Valparaíso, Chile
- Biomedical Imaging Center, Pontificia Universidad Catolica de Chile, Santiago, Chile
- Millennium Institute for Intelligent Healthcare Engineering - iHEALTH, Santiago, Chile
| | - Matthias Stuber
- Département de Radiologie Médicale, Centre Hospitalier Universitaire Vaudois, Lausanne, Switzerland
| | - Ali B Syed
- Department of Radiology, Stanford University, Stanford, CA, USA
| | - Johannes Töger
- Clinical Physiology, Department of Clinical Sciences Lund, Lund University, Skåne University Hospital, Lund, Sweden
| | - Rob van der Geest
- Division of Image Processing, Department of Radiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Jos Westenberg
- CardioVascular Imaging Group (CVIG), Department of Radiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Liang Zhong
- National Heart Centre Singapore, Duke-NUS Medical School, National University of Singapore, Singapore, Singapore
| | - Yumin Zhong
- Department of Radiology, School of Medicine, Shanghai Children's Medical Center Affiliated With Shanghai Jiao Tong University, Shanghai, People's Republic of China
| | - Oliver Wieben
- Departments of Radiology and Medical Physics, University of Wisconsin, Madison, WI, USA
| | - Petter Dyverfeldt
- Department of Health, Medicine and Caring Sciences, Linköping University, Linköping, Sweden
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
| |
Collapse
|
27
|
Kraft M, Ryger S, Berman BP, Downs ME, Jordanova KV, Poorman ME, Oberdick SD, Ogier SE, Russek SE, Dagher J, Keenan KE. Towards a barrier-free anthropomorphic brain phantom for quantitative magnetic resonance imaging: Design, first construction attempt, and challenges. PLoS One 2023; 18:e0285432. [PMID: 37437022 DOI: 10.1371/journal.pone.0285432] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2023] [Accepted: 04/21/2023] [Indexed: 07/14/2023] Open
Abstract
Existing magnetic resonance imaging (MRI) reference objects, or phantoms, are typically constructed from simple liquid or gel solutions in containers with specific geometric configurations to enable multi-year stability. However, there is a need for phantoms that better mimic the human anatomy without barriers between the tissues. Barriers result in regions without MRI signal between the different tissue mimics, which is an artificial image artifact. We created an anatomically representative 3D structure of the brain that mimicked the T1 and T2 relaxation properties of white and gray matter at 3 T. While the goal was to avoid barriers between tissues, the 3D printed barrier between white and gray matter and other flaws in the construction were visible at 3 T. Stability measurements were made using a portable MRI system operating at 64 mT, and T2 relaxation time was stable from 0 to 22 weeks. The phantom T1 relaxation properties did change from 0 to 10 weeks; however, they did not substantially change between 10 weeks and 22 weeks. The anthropomorphic phantom used a dissolvable mold construction method to better mimic anatomy, which worked in small test objects. The construction process, though, had many challenges. We share this work with the hope that the community can build on our experience.
Collapse
Affiliation(s)
- Mikail Kraft
- National Institute of Standards and Technology, Physical Measurement Laboratory, Boulder, Colorado, United States of America
| | - Slavka Ryger
- National Institute of Standards and Technology, Physical Measurement Laboratory, Boulder, Colorado, United States of America
| | - Ben P Berman
- The MITRE Corporation, McLean, Virginia, United States of America
| | - Matthew E Downs
- The MITRE Corporation, McLean, Virginia, United States of America
| | - Kalina V Jordanova
- National Institute of Standards and Technology, Physical Measurement Laboratory, Boulder, Colorado, United States of America
| | - Megan E Poorman
- Hyperfine, Inc, Guilford, Connecticut, United States of America
| | - Samuel D Oberdick
- National Institute of Standards and Technology, Physical Measurement Laboratory, Boulder, Colorado, United States of America
- Department of Physics, University of Colorado, Boulder, Colorado, United States of America
| | - Stephen E Ogier
- National Institute of Standards and Technology, Physical Measurement Laboratory, Boulder, Colorado, United States of America
- Department of Physics, University of Colorado, Boulder, Colorado, United States of America
| | - Stephen E Russek
- National Institute of Standards and Technology, Physical Measurement Laboratory, Boulder, Colorado, United States of America
| | - Joseph Dagher
- The MITRE Corporation, McLean, Virginia, United States of America
| | - Kathryn E Keenan
- National Institute of Standards and Technology, Physical Measurement Laboratory, Boulder, Colorado, United States of America
| |
Collapse
|
28
|
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
|
29
|
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
|
30
|
Kiersnowski OC, Karsa A, Wastling SJ, Thornton JS, Shmueli K. Investigating the effect of oblique image acquisition on the accuracy of QSM and a robust tilt correction method. Magn Reson Med 2023; 89:1791-1808. [PMID: 36480002 PMCID: PMC10953050 DOI: 10.1002/mrm.29550] [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: 09/23/2022] [Revised: 10/28/2022] [Accepted: 11/16/2022] [Indexed: 12/13/2022]
Abstract
PURPOSE Quantitative susceptibility mapping (QSM) is used increasingly for clinical research where oblique image acquisition is commonplace, but its effects on QSM accuracy are not well understood. THEORY AND METHODS The QSM processing pipeline involves defining the unit magnetic dipole kernel, which requires knowledge of the direction of the main magnetic fieldB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ with respect to the acquired image volume axes. The direction ofB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ is dependent on the axis and angle of rotation in oblique acquisition. Using both a numerical brain phantom and in vivo acquisitions in 5 healthy volunteers, we analyzed the effects of oblique acquisition on magnetic susceptibility maps. We compared three tilt-correction schemes at each step in the QSM pipeline: phase unwrapping, background field removal and susceptibility calculation, using the RMS error and QSM-tuned structural similarity index. RESULTS Rotation of wrapped phase images gave severe artifacts. Background field removal with projection onto dipole fields gave the most accurate susceptibilities when the field map was first rotated into alignment withB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ . Laplacian boundary value and variable-kernel sophisticated harmonic artifact reduction for phase data background field removal methods gave accurate results without tilt correction. For susceptibility calculation, thresholded k-space division, iterative Tikhonov regularization, and weighted linear total variation regularization, all performed most accurately when local field maps were rotated into alignment withB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ before susceptibility calculation. CONCLUSION For accurate QSM, oblique acquisition must be taken into account. Rotation of images into alignment withB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ should be carried out after phase unwrapping and before background-field removal. We provide open-source tilt-correction code to incorporate easily into existing pipelines: https://github.com/o-snow/QSM_TiltCorrection.git.
Collapse
Affiliation(s)
- Oliver C. Kiersnowski
- Department of Medical Physics and Biomedical EngineeringUniversity College LondonLondonUnited Kingdom
| | - Anita Karsa
- Department of Medical Physics and Biomedical EngineeringUniversity College LondonLondonUnited Kingdom
| | - Stephen J. Wastling
- Neuroradiological Academic UnitUCL Queen Square Institute of NeurologyLondonUnited Kingdom
- Lysholm Department of NeuroradiologyNational Hospital for Neurology and NeurosurgeryLondonUnited Kingdom
| | - John S. Thornton
- Neuroradiological Academic UnitUCL Queen Square Institute of NeurologyLondonUnited Kingdom
- Lysholm Department of NeuroradiologyNational Hospital for Neurology and NeurosurgeryLondonUnited Kingdom
| | - Karin Shmueli
- Department of Medical Physics and Biomedical EngineeringUniversity College LondonLondonUnited Kingdom
| |
Collapse
|
31
|
D’Astous A, Cereza G, Papp D, Gilbert KM, Stockmann JP, Alonso-Ortiz E, Cohen-Adad J. Shimming toolbox: An open-source software toolbox for B0 and B1 shimming in MRI. Magn Reson Med 2023; 89:1401-1417. [PMID: 36441743 PMCID: PMC9910837 DOI: 10.1002/mrm.29528] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2022] [Revised: 10/21/2022] [Accepted: 10/25/2022] [Indexed: 11/29/2022]
Abstract
PURPOSE Introduce Shimming Toolbox ( https://shimming-toolbox.org), an open-source software package for prototyping new methods and performing static, dynamic, and real-time B0 shimming as well as B1 shimming experiments. METHODS Shimming Toolbox features various field mapping techniques, manual and automatic masking for the brain and spinal cord, B0 and B1 shimming capabilities accessible through a user-friendly graphical user interface. Validation of Shimming Toolbox was demonstrated in three scenarios: (i) B0 dynamic shimming in the brain at 7T using custom AC/DC coils, (ii) B0 real-time shimming in the spinal cord at 3T, and (iii) B1 static shimming in the spinal cord at 7T. RESULTS The B0 dynamic shimming of the brain at 7T took about 10 min to perform. It showed a 47% reduction in the standard deviation of the B0 field, associated with noticeable improvements in geometric distortions in EPI images. Real-time dynamic xyz-shimming in the spinal cord took about 5 min and showed a 30% reduction in the standard deviation of the signal distribution. B1 static shimming experiments in the spinal cord took about 10 min to perform and showed a 40% reduction in the coefficient of variation of the B1 field. CONCLUSION Shimming Toolbox provides an open-source platform where researchers can collaborate, prototype and conveniently test B0 and B1 shimming experiments. Future versions will include additional field map preprocessing techniques, optimization algorithms, and compatibility across multiple MRI manufacturers.
Collapse
Affiliation(s)
- Alexandre D’Astous
- NeuroPoly Lab, Institute of Biomedical Engineering,
Polytechnique Montréal, Montréal, QC, Canada
| | - Gaspard Cereza
- NeuroPoly Lab, Institute of Biomedical Engineering,
Polytechnique Montréal, Montréal, QC, Canada
| | - Daniel Papp
- NeuroPoly Lab, Institute of Biomedical Engineering,
Polytechnique Montréal, Montréal, QC, Canada
| | - Kyle M. Gilbert
- Centre for Functional and Metabolic Mapping, The
University of Western Ontario, London, Ontario, Canada
| | - Jason P. Stockmann
- Athinoula A. Martinos Center for Biomedical Imaging,
Massachusetts General Hospital, Charlestown, MA, United States
- Harvard Medical School, Boston, MA, United States
| | - Eva Alonso-Ortiz
- NeuroPoly Lab, Institute of Biomedical Engineering,
Polytechnique Montréal, Montréal, QC, Canada
| | - Julien Cohen-Adad
- NeuroPoly Lab, Institute of Biomedical Engineering,
Polytechnique Montréal, Montréal, QC, Canada
- Functional Neuroimaging Unit, CRIUGM, Université de
Montréal, Montréal, QC, Canada
- Mila - Quebec AI Institute, Montréal, QC,
Canada
- Centre de recherche du CHU Sainte-Justine,
Université de Montréal, Montréal, QC, Canada
| |
Collapse
|
32
|
Chen Z, Zhai X, Chen Z. Tilted quantitative susceptibility mapping at oblique MRI (tiltQSM). Comput Biol Med 2023; 157:106802. [PMID: 36965324 DOI: 10.1016/j.compbiomed.2023.106802] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2023] [Revised: 03/05/2023] [Accepted: 03/20/2023] [Indexed: 03/27/2023]
Abstract
OBJECTIVE If the phase image matrix was acquired from oblique MRI, it is needed to deal with the oblique effect for quantitative susceptibility mapping (QSM), as addressed in this paper. METHODS We proposed two methods for QSM reconstruction from slice-tilted MRI phase image (tiltQSM): 1) rotData per anti-tilting phase image rotation back into the B0-upright system, and 2) rotKernel per pro-tilting dipole kernel rotation into the same oblique setting as defined by the tilted phase image. Both matrix methods were implemented in an additional preprocessing subroutine to ensure that the phase image and the dipole kernel were represented in the same coordinate system (either in B0-upright system or in B0-tilted system); thereafter tiltQSM could be completed through a regular QSM procedure. Besides the oblique effect, tiltQSM also suffers from MRI anisotropy. We provided numeric simulations, phantom tests and in vivo brain experiments on tiltQSM with oblique MRI (axial slice tilting at 3T). RESULTS The tiltQSM reconstruction could attain a performance corr > 0.90 (spatial correlation conformance) for small tilting angles <10°. The tiltQSM performance could be further degraded by voxel anisotropy due to image matrix rotation (digital geometry error). CONCLUSIONS To seek inverse solutions of MRI phase images acquired at oblique MRI (e.g. in axial slice tilting), we proposed tiltQSM to deal with the oblique effect per matrix rotation (either rotData or rotKernel) in a preprocessing subroutine prior to a regular QSM procedure. In practice, it is always recommended to acquire MRI phase images in isotropic matrix at zero obliqueness (or limited to small tilting angles <10°) for maximal (optimal) QSM reconstruction.
Collapse
Affiliation(s)
- Zeyuan Chen
- Department of Computer Sciences, University of California-Davis, Davis, CA, USA; Microsoft Corporation, Seattle, WA, USA
| | | | - Zikuan Chen
- Zinv LLC, Albuquerque, NM, USA; Diagnostic Radiology, City of Hope National Medical Center, Duarte, CA, USA.
| |
Collapse
|
33
|
Renton AI, Dao TT, Johnstone T, Civier O, Sullivan RP, White DJ, Lyons P, Slade BM, Abbott DF, Amos TJ, Bollmann S, Botting A, Campbell MEJ, Chang J, Close TG, Eckstein K, Egan GF, Evas S, Flandin G, Garner KG, Garrido MI, Ghosh SS, Grignard M, Hannan AJ, Huber R, Kaczmarzyk JR, Kasper L, Kuhlmann L, Lou K, Mantilla-Ramos YJ, Mattingley JB, Morris J, Narayanan A, Pestilli F, Puce A, Ribeiro FL, Rogasch NC, Rorden C, Schira M, Shaw TB, Sowman PF, Spitz G, Stewart A, Ye X, Zhu JD, Hughes ME, Narayanan A, Bollmann S. Neurodesk: An accessible, flexible, and portable data analysis environment for reproducible neuroimaging. RESEARCH SQUARE 2023:rs.3.rs-2649734. [PMID: 36993557 PMCID: PMC10055538 DOI: 10.21203/rs.3.rs-2649734/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Neuroimaging data analysis often requires purpose-built software, which can be challenging to install and may produce different results across computing environments. Beyond being a roadblock to neuroscientists, these issues of accessibility and portability can hamper the reproducibility of neuroimaging data analysis pipelines. Here, we introduce the Neurodesk platform, which harnesses software containers to support a comprehensive and growing suite of neuroimaging software (https://www.neurodesk.org/). Neurodesk includes a browser-accessible virtual desktop environment and a command line interface, mediating access to containerized neuroimaging software libraries on various computing platforms, including personal and high-performance computers, cloud computing and Jupyter Notebooks. This community-oriented, open-source platform enables a paradigm shift for neuroimaging data analysis, allowing for accessible, flexible, fully reproducible, and portable data analysis pipelines.
Collapse
Affiliation(s)
- Angela I. Renton
- The University of Queensland, Queensland Brain Institute, St Lucia 4072, Australia
- The University of Queensland, School of Information Technology and Electrical Engineering, St Lucia 4072, Australia
| | - Thuy T. Dao
- The University of Queensland, School of Information Technology and Electrical Engineering, St Lucia 4072, Australia
| | - Tom Johnstone
- Centre for Mental Health & Brain Sciences, Swinburne University of Technology, Hawthorn 3122, Australia
| | - Oren Civier
- Centre for Mental Health & Brain Sciences, Swinburne University of Technology, Hawthorn 3122, Australia
| | - Ryan P. Sullivan
- The University of Sydney, School of Biomedical Engineering, Sydney, Australia
| | - David J. White
- Centre for Mental Health & Brain Sciences, Swinburne University of Technology, Hawthorn 3122, Australia
| | - Paris Lyons
- Centre for Mental Health & Brain Sciences, Swinburne University of Technology, Hawthorn 3122, Australia
| | - Benjamin M. Slade
- Centre for Mental Health & Brain Sciences, Swinburne University of Technology, Hawthorn 3122, Australia
| | - David F. Abbott
- The Florey Institute of Neuroscience and Mental Health, The University of Melbourne, Victoria, Australia
| | - Toluwani J. Amos
- School of Life Science and Technology, University of Electronic Science and Technology, China
| | - Saskia Bollmann
- The University of Queensland, School of Information Technology and Electrical Engineering, St Lucia 4072, Australia
| | - Andy Botting
- Australian Research Data Commons (ARDC), Australia
| | - Megan E. J. Campbell
- School of Psychological Sciences, University of Newcastle, Australia
- Hunter Medical Research Institute Imaging Centre, Newcastle, Australia
| | - Jeryn Chang
- The University of Queensland, School of Biomedical Sciences, St Lucia 4072, Australia
| | - Thomas G. Close
- The University of Sydney, School of Biomedical Engineering, Sydney, Australia
| | - Korbinian Eckstein
- The University of Queensland, School of Information Technology and Electrical Engineering, St Lucia 4072, Australia
| | - Gary F. Egan
- The Turner Institute for Brain and Mental Health, School of Psychological Sciences, Monash University, Victoria, Australia
- Monash Biomedical Imaging, Monash University, Victoria, Australia
| | - Stefanie Evas
- School of Psychology, University of Adelaide, Adelaide, 5000, Australia
- Human Health, Health & Biosecurity, CSIRO, Adelaide, 5000, Australia
| | - Guillaume Flandin
- Wellcome Centre for Human Neuroimaging, University College London, London, UK
| | - Kelly G. Garner
- The University of Queensland, Queensland Brain Institute, St Lucia 4072, Australia
- The University of Queensland, School of Psychology, St Lucia 4072, Australia
| | - Marta I. Garrido
- Melbourne School of Psychological Sciences, The University of Melbourne
- Graeme Clark Institute for Biomedical Engineering, The University of Melbourne
| | - Satrajit S. Ghosh
- McGovern Institute for Brain Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
- Department of Otolaryngology - Head and Neck Surgery, Harvard Medical School, Boston, MA, USA
| | - Martin Grignard
- GIGA CRC In-Vivo Imaging, University of Liege, Liege, Belgium
| | - Anthony J. Hannan
- The Florey Institute of Neuroscience and Mental Health, The University of Melbourne, Victoria, Australia
| | - Renzo Huber
- Functional Magnetic Resonance Imaging Core Facility (FMRIF), National Institute of Mental Health (NIMH), USA
| | - Jakub R. Kaczmarzyk
- Medical Scientist Training Program, Stony Brook University, Stony Brook, NY, United States of America
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, United States of America
| | - Lars Kasper
- Techna Institute, University Health Network, Toronto, Canada
| | - Levin Kuhlmann
- Department of Data Science and AI, Faculty of Information Technology, Monash University, Clayton VIC 3800, Australia
| | - Kexin Lou
- The University of Queensland, School of Information Technology and Electrical Engineering, St Lucia 4072, Australia
| | | | - Jason B. Mattingley
- The University of Queensland, Queensland Brain Institute, St Lucia 4072, Australia
- The University of Queensland, School of Psychology, St Lucia 4072, Australia
| | - Jo Morris
- Australian Research Data Commons (ARDC), Australia
| | | | - Franco Pestilli
- Department of Psychology, Center for Perceptual Systems, Center for Theoretical and Computational Neuroscience, Center on Aging and Population Sciences, Center for Learning and Memory, The University of Texas at Austin, 108 E Dean Keeton St, Austin, TX 78712, USA
| | - Aina Puce
- Department of Psychological and Brain Sciences, Indiana University, Bloomington, IN 47405, USA
| | - Fernanda L. Ribeiro
- The University of Queensland, School of Information Technology and Electrical Engineering, St Lucia 4072, Australia
| | - Nigel C. Rogasch
- Discipline of Psychiatry, Adelaide Medical School, University of Adelaide, Australia
- Hopwood Centre for Neurobiology, Lifelong Health Theme, South Australian Health and Medical Research Institute (SAHMRI), Adelaide, SA, Australia
- The Turner Institute for Brain and Mental Health, School of Psychological Sciences, Monash University, Victoria, Australia
| | - Chris Rorden
- McCausland Center for Brain Imaging, Department of Psychology, University of South Carolina, Columbia SC, 29208, USA
| | - Mark Schira
- School of Psychology, University of Wollongong, Wollongong, 2522, Australia
| | - Thomas B. Shaw
- The University of Queensland, School of Information Technology and Electrical Engineering, St Lucia 4072, Australia
- The University of Queensland, Centre for Advanced Imaging, St Lucia 4072, Australia
- Department of Neurology, Royal Brisbane and Women's Hospital, Brisbane, Australia
| | - Paul F. Sowman
- Macquarie University, School of Psychological Sciences, North Ryde 2112, Australia
| | - Gershon Spitz
- Department of Neuroscience, Central Clinical School, Faculty of Medicine, Nursing and Health Sciences, Monash University, Australia
- Monash-Epworth Rehabilitation Research Centre, Turner Institute for Brain and Mental Health, School of Psychological Sciences, Monash University, Clayton, 3168, Australia
| | - Ashley Stewart
- The University of Queensland, School of Information Technology and Electrical Engineering, St Lucia 4072, Australia
| | - Xincheng Ye
- The University of Queensland, School of Information Technology and Electrical Engineering, St Lucia 4072, Australia
| | - Judy D. Zhu
- Macquarie University, School of Psychological Sciences, North Ryde 2112, Australia
| | - Matthew E. Hughes
- Centre for Mental Health & Brain Sciences, Swinburne University of Technology, Hawthorn 3122, Australia
| | - Aswin Narayanan
- The University of Queensland, Centre for Advanced Imaging, St Lucia 4072, Australia
| | - Steffen Bollmann
- The University of Queensland, School of Information Technology and Electrical Engineering, St Lucia 4072, Australia
| |
Collapse
|
34
|
Barbieri M, Chaudhari AS, Moran CJ, Gold GE, Hargreaves BA, Kogan F. A method for measuring B 0 field inhomogeneity using quantitative double-echo in steady-state. Magn Reson Med 2023; 89:577-593. [PMID: 36161727 PMCID: PMC9712261 DOI: 10.1002/mrm.29465] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2022] [Revised: 08/31/2022] [Accepted: 09/01/2022] [Indexed: 12/13/2022]
Abstract
PURPOSE To develop and validate a method forB 0 $$ {B}_0 $$ mapping for knee imaging using the quantitative Double-Echo in Steady-State (qDESS) exploiting the phase difference (Δ θ $$ \Delta \theta $$ ) between the two echoes acquired. Contrary to a two-gradient-echo (2-GRE) method,Δ θ $$ \Delta \theta $$ depends only on the first echo time. METHODS Bloch simulations were applied to investigate robustness to noise of the proposed methodology and all imaging studies were validated with phantoms and in vivo simultaneous bilateral knee acquisitions. Two phantoms and five healthy subjects were scanned using qDESS, water saturation shift referencing (WASSR), and multi-GRE sequences.Δ B 0 $$ \Delta {B}_0 $$ maps were calculated with the qDESS and the 2-GRE methods and compared against those obtained with WASSR. The comparison was quantitatively assessed exploiting pixel-wise difference maps, Bland-Altman (BA) analysis, and Lin's concordance coefficient (ρ c $$ {\rho}_c $$ ). For in vivo subjects, the comparison was assessed in cartilage using average values in six subregions. RESULTS The proposed method for measuringΔ B 0 $$ \Delta {B}_0 $$ inhomogeneities from a qDESS acquisition providedΔ B 0 $$ \Delta {B}_0 $$ maps that were in good agreement with those obtained using WASSR.Δ B 0 $$ \Delta {B}_0 $$ ρ c $$ {\rho}_c $$ values were≥ $$ \ge $$ 0.98 and 0.90 in phantoms and in vivo, respectively. The agreement between qDESS and WASSR was comparable to that of a 2-GRE method. CONCLUSION The proposed method may allow B0 correction for qDESST 2 $$ {T}_2 $$ mapping using an inherently co-registeredΔ B 0 $$ \Delta {B}_0 $$ map without requiring an additional B0 measurement sequence. More generally, the method may help shorten knee imaging protocols that require an auxiliaryΔ B 0 $$ \Delta {B}_0 $$ map by exploiting a qDESS acquisition that also providesT 2 $$ {T}_2 $$ measurements and high-quality morphological imaging.
Collapse
Affiliation(s)
- Marco Barbieri
- Department of Radiology, Stanford University, Stanford, CA, U.S.A
| | - Akshay S. Chaudhari
- Department of Radiology, Stanford University, Stanford, CA, U.S.A
- Department of Biomedical Data Science, Stanford University, Stanford, CA, U.S.A
| | | | - Garry E. Gold
- Department of Radiology, Stanford University, Stanford, CA, U.S.A
- Department of Bioengineering, Stanford University, Stanford, CA, U.S.A
| | - Brian A. Hargreaves
- Department of Radiology, Stanford University, Stanford, CA, U.S.A
- Department of Bioengineering, Stanford University, Stanford, CA, U.S.A
- Department of Electrical Engineering, Stanford University, Stanford, CA, U.S.A
| | - Feliks Kogan
- Department of Radiology, Stanford University, Stanford, CA, U.S.A
| |
Collapse
|
35
|
Kames C, Doucette J, Rauscher A. Multi-echo dipole inversion for magnetic susceptibility mapping. Magn Reson Med 2023; 89:2391-2401. [PMID: 36695283 DOI: 10.1002/mrm.29588] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2022] [Revised: 12/08/2022] [Accepted: 12/31/2022] [Indexed: 01/26/2023]
Abstract
PURPOSE Reconstructing tissue magnetic susceptibility (QSM) from MRI phase data involves solving multiple consecutive ill-posed inverse problems such as phase unwrapping, background field removal, and field-to-source inversion. Multi-echo acquisitions present an additional challenge, as the magnetization field is typically computed from the multiple phase data prior to reconstructing the susceptibility map. Processing the multiple phase data introduces errors during the field estimation, violating assumptions of the subsequent inverse problems, manifesting as streaking artifacts in the susceptibility map. To address this challenge, we propose a multi-echo field-to-source forward model that forgoes the field estimation step. Moreover, we propose a fully general underestimation correction step to recover susceptibility sources that were regularized away during the field-to-source inversion. METHODS The multi-echo forward model and correction step were validated on the QSM Challenge 2.0 datasets and compared to the standard single field-to-source model in in vivo human brains using different types of deconvolution algorithms. RESULTS On the QSM Challenge 2.0 datasets the multi-echo forward model and correction step attain state-of-the-art results on all metrics by a wide margin. Experiments in in vivo brains show that the multi-echo model is in agreement with the single field-to-source model and that the proposed forward model and correction step can be used with any available dipole inversion method. CONCLUSION A multi-echo field-to-source forward model forgoes the need to fit multi-echo phase data and achieves state-of-the-art results on the QSM Challenge 2.0 data. Underestimated low-frequency susceptibility distributions can be partially recovered using a correction step.
Collapse
Affiliation(s)
- Christian Kames
- UBC MRI Research Centre, The University of British Columbia, Vancouver, British Columbia, Canada.,Department of Physics and Astronomy, The University of British Columbia, Vancouver, British Columbia, Canada
| | - Jonathan Doucette
- UBC MRI Research Centre, The University of British Columbia, Vancouver, British Columbia, Canada.,Department of Physics and Astronomy, The University of British Columbia, Vancouver, British Columbia, Canada
| | - Alexander Rauscher
- UBC MRI Research Centre, The University of British Columbia, Vancouver, British Columbia, Canada.,Department of Physics and Astronomy, The University of British Columbia, Vancouver, British Columbia, Canada.,Department of Pediatrics, The University of British Columbia, Vancouver, British Columbia, Canada
| |
Collapse
|
36
|
Bancelin D, Bachrata B, Bollmann S, de Lima Cardoso P, Szomolanyi P, Trattnig S, Robinson SD. Unsupervised physiological noise correction of functional magnetic resonance imaging data using phase and magnitude information (PREPAIR). Hum Brain Mapp 2022; 44:1209-1226. [PMID: 36401844 PMCID: PMC9875918 DOI: 10.1002/hbm.26152] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2022] [Revised: 09/29/2022] [Accepted: 10/23/2022] [Indexed: 11/21/2022] Open
Abstract
Of the sources of noise affecting blood oxygen level-dependent functional magnetic resonance imaging (fMRI), respiration and cardiac fluctuations are responsible for the largest part of the variance, particularly at high and ultrahigh field. Existing approaches to removing physiological noise either use external recordings, which can be unwieldy and unreliable, or attempt to identify physiological noise from the magnitude fMRI data. Data-driven approaches are limited by sensitivity, temporal aliasing, and the need for user interaction. In the light of the sensitivity of the phase of the MR signal to local changes in the field stemming from physiological processes, we have developed an unsupervised physiological noise correction method using the information carried in the phase and the magnitude of echo-planar imaging data. Our technique, Physiological Regressor Estimation from Phase and mAgnItude, sub-tR (PREPAIR) derives time series signals sampled at the slice TR from both phase and magnitude images. It allows physiological noise to be captured without aliasing, and efficiently removes other sources of signal fluctuations not related to physiology, prior to regressor estimation. We demonstrate that the physiological signal time courses identified with PREPAIR agree well with those from external devices and retrieve challenging cardiac dynamics. The removal of physiological noise was as effective as that achieved with the most used approach based on external recordings, RETROICOR. In comparison with widely used recording-free physiological noise correction tools-PESTICA and FIX, both performed in unsupervised mode-PREPAIR removed significantly more respiratory and cardiac noise than PESTICA, and achieved a larger increase in temporal signal-to-noise-ratio at both 3 and 7 T.
Collapse
Affiliation(s)
- David Bancelin
- High Field MR Centre, Department of Biomedical Imaging and Image‐guided TherapyMedical University of ViennaViennaAustria
| | - Beata Bachrata
- High Field MR Centre, Department of Biomedical Imaging and Image‐guided TherapyMedical University of ViennaViennaAustria,Karl Landsteiner Institute for Clinical Molecular MR in Musculoskeletal ImagingViennaAustria
| | - Saskia Bollmann
- Centre for Advanced ImagingThe University of QueenslandBrisbaneAustralia
| | - Pedro de Lima Cardoso
- High Field MR Centre, Department of Biomedical Imaging and Image‐guided TherapyMedical University of ViennaViennaAustria
| | - Pavol Szomolanyi
- High Field MR Centre, Department of Biomedical Imaging and Image‐guided TherapyMedical University of ViennaViennaAustria
| | - Siegfried Trattnig
- High Field MR Centre, Department of Biomedical Imaging and Image‐guided TherapyMedical University of ViennaViennaAustria,Karl Landsteiner Institute for Clinical Molecular MR in Musculoskeletal ImagingViennaAustria
| | - Simon Daniel Robinson
- High Field MR Centre, Department of Biomedical Imaging and Image‐guided TherapyMedical University of ViennaViennaAustria,Karl Landsteiner Institute for Clinical Molecular MR in Musculoskeletal ImagingViennaAustria,Centre for Advanced ImagingThe University of QueenslandBrisbaneAustralia,Department of NeurologyMedical University of GrazGrazAustria
| |
Collapse
|
37
|
Hagberg GE, Eckstein K, Tuzzi E, Zhou J, Robinson S, Scheffler K. Phase-based masking for quantitative susceptibility mapping of the human brain at 9.4T. Magn Reson Med 2022; 88:2267-2276. [PMID: 35754142 PMCID: PMC7613679 DOI: 10.1002/mrm.29368] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2021] [Revised: 05/05/2022] [Accepted: 05/31/2022] [Indexed: 12/15/2022]
Abstract
PURPOSE To develop improved tissue masks for QSM. METHODS Masks including voxels at the brain surface were automatically generated from the magnitude alone (MM) or combined with test functions from the first (PG) or second (PB) derivative of the sign of the wrapped phase. Phase images at 3T and 9.4T were simulated at different TEs and used to generate a mask, PItoh , with between-voxel phase differences less than π. MM, PG, and PB were compared with PItoh . QSM were generated from 3D multi-echo gradient-echo data acquired at 9.4T (21 subjects aged: 20-56y), and from the QSM2016 challenge 3T data using different masks, unwrapping, background removal, and dipole inversion algorithms. QSM contrast was quantified using age-based iron concentrations. RESULTS Close to air cavities, phase wraps became denser with increasing field and echo time, yielding increased values of the test functions. Compared with PItoh , PB had the highest Dice coefficient, while PG had the lowest and MM the highest percentage of voxels outside PItoh. Artifacts observed in QSM at 9.4T with MM were mitigated by stronger background filters but yielded a reduced QSM contrast. With PB, QSM contrast was greater and artifacts diminished. Similar results were obtained with challenge data, evidencing larger effects of mask close to air cavities. CONCLUSION Automatic, phase-based masking founded on the second derivative of the sign of the wrapped phase, including cortical voxels at the brain surface, was able to mitigate artifacts and restore QSM contrast across cortical and subcortical brain regions.
Collapse
Affiliation(s)
- Gisela E. Hagberg
- Department for Biomedical Magnetic Resonance, University Hospital Tübingen, Tübingen, Germany
- High Field Magnetic Resonance Center, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Korbinian Eckstein
- High Field MR Centre, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
- High Field Magnetic Resonance Center, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Elisa Tuzzi
- Department for Biomedical Magnetic Resonance, University Hospital Tübingen, Tübingen, Germany
- High Field Magnetic Resonance Center, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Jiazheng Zhou
- Department for Biomedical Magnetic Resonance, University Hospital Tübingen, Tübingen, Germany
- High Field Magnetic Resonance Center, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Simon Robinson
- Department of Neurology, Medical University of Graz, Graz, Austria
- High Field MR Centre, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
- Centre for Advanced Imaging, The University of Queensland, Brisbane, Australia
| | - Klaus Scheffler
- Department for Biomedical Magnetic Resonance, University Hospital Tübingen, Tübingen, Germany
- High Field Magnetic Resonance Center, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| |
Collapse
|
38
|
Daval-Frérot G, Massire A, Mailhe B, Nadar M, Vignaud A, Ciuciu P. Iterative static field map estimation for off-resonance correction in non-Cartesian susceptibility weighted imaging. Magn Reson Med 2022; 88:1592-1607. [PMID: 35735217 PMCID: PMC9545844 DOI: 10.1002/mrm.29297] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2021] [Revised: 04/01/2022] [Accepted: 04/19/2022] [Indexed: 11/21/2022]
Abstract
Purpose Patient‐induced inhomogeneities in the magnetic field cause distortions and blurring during acquisitions with long readouts such as in susceptibility‐weighted imaging (SWI). Most correction methods require collecting an additional ΔB0 field map to remove these artifacts. Theory The static ΔB0 field map can be approximated with an acceptable error directly from a single echo acquisition in SWI. The main component of the observed phase is linearly related to ΔB0 and the echo time (TE), and the relative impact of non‐ ΔB0 terms becomes insignificant with TE >20 ms at 3 T for a well‐tuned system. Methods The main step is to combine and unfold the multi‐channel phase maps wrapped many times, and several competing algorithms are compared for this purpose. Four in vivo brain data sets collected using the recently proposed 3D spreading projection algorithm for rapid k‐space sampling (SPARKLING) readouts are used to assess the proposed method. Results The estimated 3D field maps generated with a 0.6 mm isotropic spatial resolution provide overall similar off‐resonance corrections compared to reference corrections based on an external ΔB0 acquisitions, and even improved for 2 of 4 individuals. Although a small estimation error is expected, no aftermath was observed in the proposed corrections, whereas degradations were observed in the references. Conclusion A static ΔB0 field map estimation method was proposed to take advantage of acquisitions with long echo times, and outperformed the reference technique based on an external field map. The difference can be attributed to an inherent robustness to mismatches between volumes and external ΔB0 maps, and diverse other sources investigated. Click here for author‐reader discussions
Collapse
Affiliation(s)
- Guillaume Daval-Frérot
- Siemens Healthcare SAS, Saint-Denis, France.,CEA, NeuroSpin, CNRS, Université Paris-Saclay, Gif-sur-Yvette, France.,Inria, Palaiseau, France
| | | | - Boris Mailhe
- Siemens Healthineers, Digital Technology & Innovation, Princeton, New Jersey, USA
| | - Mariappan Nadar
- Siemens Healthineers, Digital Technology & Innovation, Princeton, New Jersey, USA
| | - Alexandre Vignaud
- CEA, NeuroSpin, CNRS, Université Paris-Saclay, Gif-sur-Yvette, France
| | - Philippe Ciuciu
- CEA, NeuroSpin, CNRS, Université Paris-Saclay, Gif-sur-Yvette, France.,Inria, Palaiseau, France
| |
Collapse
|
39
|
Bachrata B, Trattnig S, Robinson SD. Quantitative susceptibility mapping of the head-and-neck using SMURF fat-water imaging with chemical shift and relaxation rate corrections. Magn Reson Med 2022; 87:1461-1479. [PMID: 34850446 PMCID: PMC7612304 DOI: 10.1002/mrm.29069] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2021] [Revised: 09/23/2021] [Accepted: 10/15/2021] [Indexed: 12/19/2022]
Abstract
PURPOSE To address the challenges posed by fat-water chemical shift artifacts and relaxation rate discrepancies to quantitative susceptibility mapping (QSM) outside the brain, and to generate accurate susceptibility maps of the head-and-neck at 3 and 7 Tesla. METHODS Simultaneous Multiple Resonance Frequency (SMURF) imaging was extended to 7 Tesla and used to acquire head-and-neck gradient echo images at both 3 and 7 Tesla. Separated fat and water images were corrected for Type 1 (displacement) and Type 2 (phase discrepancy) chemical shift artefacts, and for the bias resulting from differences in T1 and T 2 ∗ relaxation rates, recombined and used as the basis for QSM. A novel phase signal-based masking approach was used to generate head-and-neck masks. RESULTS SMURF generated well-separated fat and water images of the head-and-neck. Corrections for chemical shift artefacts and relaxation rate differences removed overestimation of the susceptibility values, blurring in the susceptibility maps, and the disproportionate influence of fat in mixed voxels. The resulting susceptibility maps showed high correspondence between the paramagnetic areas and the locations of fatty tissues and the susceptibility estimates were similar to literature values. The proposed masking approach was shown to provide a simple means of generating head-and-neck masks. CONCLUSION Corrections for Type 1 and Type 2 chemical shift artefacts and for fat-water relaxation rate differences, mainly in T1 , were shown to be required for accurate susceptibility mapping of fatty-body regions. SMURF made it possible to apply these corrections and generate high-quality susceptibility maps of the entire head-and-neck at both 3 and 7 Tesla.
Collapse
Affiliation(s)
- Beata Bachrata
- High Field MR Centre, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
- Karl Landsteiner Institute for Clinical Molecular MR in Musculoskeletal Imaging, Vienna, Austria
| | - Siegfried Trattnig
- High Field MR Centre, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
- Karl Landsteiner Institute for Clinical Molecular MR in Musculoskeletal Imaging, Vienna, Austria
| | - Simon Daniel Robinson
- High Field MR Centre, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
- Karl Landsteiner Institute for Clinical Molecular MR in Musculoskeletal Imaging, Vienna, Austria
- Centre of Advanced Imaging, University of Queensland, Brisbane, Australia
- Department of Neurology, Medical University of Graz, Graz, Austria
| |
Collapse
|
40
|
Stewart AW, Robinson SD, O'Brien K, Jin J, Widhalm G, Hangel G, Walls A, Goodwin J, Eckstein K, Tourell M, Morgan C, Narayanan A, Barth M, Bollmann S. QSMxT: Robust masking and artifact reduction for quantitative susceptibility mapping. Magn Reson Med 2021; 87:1289-1300. [PMID: 34687073 DOI: 10.1002/mrm.29048] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2021] [Revised: 08/30/2021] [Accepted: 09/27/2021] [Indexed: 01/15/2023]
Abstract
PURPOSE Quantitative susceptibility mapping (QSM) estimates the spatial distribution of tissue magnetic susceptibilities from the phase of a gradient-echo signal. QSM algorithms require a signal mask to delineate regions with reliable phase for subsequent susceptibility estimation. Existing masking techniques used in QSM have limitations that introduce artifacts, exclude anatomical detail, and rely on parameter tuning and anatomical priors that narrow their application. Here, a robust masking and reconstruction procedure is presented to overcome these limitations and enable automated QSM processing. Moreover, this method is integrated within an open-source software framework: QSMxT. METHODS A robust masking technique that automatically separates reliable from less reliable phase regions was developed and combined with a two-pass reconstruction procedure that operates on the separated sources before combination, extracting more information and suppressing streaking artifacts. RESULTS Compared with standard masking and reconstruction procedures, the two-pass inversion reduces streaking artifacts caused by unreliable phase and high dynamic ranges of susceptibility sources. It is also robust across a range of acquisitions at 3 T in volunteers and phantoms, at 7 T in tumor patients, and in an in silico head phantom, with significant artifact and error reductions, greater anatomical detail, and minimal parameter tuning. CONCLUSION The two-pass masking and reconstruction procedure separates reliable from less reliable phase regions, enabling a more accurate QSM reconstruction that mitigates artifacts, operates without anatomical priors, and requires minimal parameter tuning. The technique and its integration within QSMxT makes QSM processing more accessible and robust to streaking artifacts.
Collapse
Affiliation(s)
- Ashley Wilton Stewart
- ARC Training Centre for Innovation in Biomedical Imaging Technology, University of Queensland, Brisbane, Queensland, Australia.,Centre for Advanced Imaging, University of Queensland, Brisbane, Australia
| | - Simon Daniel Robinson
- Centre for Advanced Imaging, University of Queensland, Brisbane, Australia.,Department of Neurology, Medical University of Graz, Graz, Austria.,Karl Landsteiner Institute for Clinical Molecular MR in Musculoskeletal Imaging, Vienna, Austria.,Department of Biomedical Imaging and Image-Guided Therapy, High Field MR Center, Medical University of Vienna, Vienna, Austria
| | - Kieran O'Brien
- ARC Training Centre for Innovation in Biomedical Imaging Technology, University of Queensland, Brisbane, Queensland, Australia.,Centre for Advanced Imaging, University of Queensland, Brisbane, Australia.,Siemens Healthcare Pty Ltd, Brisbane, Queensland, Australia
| | - Jin Jin
- ARC Training Centre for Innovation in Biomedical Imaging Technology, University of Queensland, Brisbane, Queensland, Australia.,Centre for Advanced Imaging, University of Queensland, Brisbane, Australia.,Siemens Healthcare Pty Ltd, Brisbane, Queensland, Australia
| | - Georg Widhalm
- Department of Neurosurgery, Medical University of Vienna, Vienna, Austria
| | - Gilbert Hangel
- Department of Biomedical Imaging and Image-Guided Therapy, High Field MR Center, Medical University of Vienna, Vienna, Austria.,Department of Neurosurgery, Medical University of Vienna, Vienna, Austria
| | - Angela Walls
- Clinical & Research Imaging Centre, South Australian Health and Medical Research Institute, Adelaide, South Australia, Australia
| | - Jonathan Goodwin
- Department of Radiation Oncology, Calvary Mater Hospital, Newcastle, New South Wales, Australia.,School of Mathematical and Physical Science, University of Newcastle, Newcastle, New South Wales, Australia
| | - Korbinian Eckstein
- Department of Biomedical Imaging and Image-Guided Therapy, High Field MR Center, Medical University of Vienna, Vienna, Austria
| | - Monique Tourell
- ARC Training Centre for Innovation in Biomedical Imaging Technology, University of Queensland, Brisbane, Queensland, Australia.,Centre for Advanced Imaging, University of Queensland, Brisbane, Australia
| | - Catherine Morgan
- School of Psychology and Centre for Brain Research, The University of Auckland, Auckland, New Zealand.,Centre of Research Excellence, Brain Research New Zealand-Rangahau Roro Aotearoa, Auckland, New Zealand.,Centre for Advanced MRI, The University of Auckland, Auckland, New Zealand
| | - Aswin Narayanan
- Centre for Advanced Imaging, University of Queensland, Brisbane, Australia
| | - Markus Barth
- ARC Training Centre for Innovation in Biomedical Imaging Technology, University of Queensland, Brisbane, Queensland, Australia.,Centre for Advanced Imaging, University of Queensland, Brisbane, Australia.,School of Information Technology and Electrical Engineering, The University of Queensland, Brisbane, Queensland, Australia
| | - Steffen Bollmann
- ARC Training Centre for Innovation in Biomedical Imaging Technology, University of Queensland, Brisbane, Queensland, Australia.,Centre for Advanced Imaging, University of Queensland, Brisbane, Australia.,School of Information Technology and Electrical Engineering, The University of Queensland, Brisbane, Queensland, Australia
| |
Collapse
|
41
|
Eckstein K, Bachrata B, Hangel G, Widhalm G, Enzinger C, Barth M, Trattnig S, Robinson SD. Improved susceptibility weighted imaging at ultra-high field using bipolar multi-echo acquisition and optimized image processing: CLEAR-SWI. Neuroimage 2021; 237:118175. [PMID: 34000407 PMCID: PMC7612087 DOI: 10.1016/j.neuroimage.2021.118175] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2021] [Revised: 04/28/2021] [Accepted: 05/13/2021] [Indexed: 02/07/2023] Open
Abstract
Purpose Susceptibility Weighted Imaging (SWI) has become established in the clinical investigation of stroke, microbleeds, tumor vascularization, calcification and iron deposition, but suffers from a number of shortcomings and artefacts. The goal of this study was to reduce the sensitivity of SWI to strong B1 and B0 inhomogeneities at ultra-high field to generate homogeneous images with increased contrast and free of common artefacts. All steps in SWI processing have been addressed −coil combination, phase unwrapping, image combination over echoes, phase filtering and homogeneity correction −and applied to an efficient bipolar multi-echo acquisition to substantially improve the quality of SWI. Principal results Our findings regarding the optimal individual processing steps lead us to propose a Contrast-weighted, Laplace-unwrapped, bipolar multi-Echo, ASPIRE-combined, homogeneous, improved Resolution SWI, or CLEAR-SWI. CLEAR-SWI was compared to two other multi-echo SWI methods and standard, single-echo SWI with the same acquisition time at 7 T in 10 healthy volunteers and with single-echo SWI in 13 patients with brain tumors. CLEAR-SWI had improved contrast-to-noise and homogeneity, reduced signal dropout and was not compromised by the artefacts which affected standard SWI in 10 out of 13 cases close to tumors (as assessed by expert raters), as well as generating T2* maps and phase images which can be used for Quantitative Susceptibility Mapping. In a comparison with other multi-echo SWI methods, CLEAR-SWI had the fewest artefacts, highest SNR and generally higher contrast-to-noise. Major conclusions CLEAR-SWI eliminates the artefacts common in standard, single-echo SWI, reduces signal dropouts and improves image homogeneity and contrast-to-noise. Applied clinically, in a study of brain tumor patients, CLEAR-SWI was free of the artefacts which affected standard, single-echo SWI.
Collapse
Affiliation(s)
- Korbinian Eckstein
- High Field MR Centre, Department of Biomedical Imaging and Image-guided Therapy, Medical University of Vienna, Vienna, Austria
| | - Beata Bachrata
- High Field MR Centre, Department of Biomedical Imaging and Image-guided Therapy, Medical University of Vienna, Vienna, Austria; Karl Landsteiner Institute for Clinical Molecular MR in Musculoskeletal Imaging, Vienna, Austria
| | - Gilbert Hangel
- High Field MR Centre, Department of Biomedical Imaging and Image-guided Therapy, Medical University of Vienna, Vienna, Austria; Department of Neurosurgery, Medical University of Vienna, Vienna, Austria
| | - Georg Widhalm
- Department of Neurosurgery, Medical University of Vienna, Vienna, Austria
| | | | - Markus Barth
- School of Information Technology and Electrical Engineering, Faculty of Engineering, Architecture and Information Technology, The University of Queensland, Brisbane, Australia
| | - Siegfried Trattnig
- High Field MR Centre, Department of Biomedical Imaging and Image-guided Therapy, Medical University of Vienna, Vienna, Austria; Karl Landsteiner Institute for Clinical Molecular MR in Musculoskeletal Imaging, Vienna, Austria
| | - Simon Daniel Robinson
- High Field MR Centre, Department of Biomedical Imaging and Image-guided Therapy, Medical University of Vienna, Vienna, Austria; Karl Landsteiner Institute for Clinical Molecular MR in Musculoskeletal Imaging, Vienna, Austria; Department of Neurology, Medical University of Graz, Graz, Austria; Centre for Advanced Imaging, The University of Queensland, Brisbane, Australia
| |
Collapse
|
42
|
Zhou H, Cheng C, Peng H, Liang D, Liu X, Zheng H, Zou C. The PHU-NET: A robust phase unwrapping method for MRI based on deep learning. Magn Reson Med 2021; 86:3321-3333. [PMID: 34272757 DOI: 10.1002/mrm.28927] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2020] [Revised: 06/22/2021] [Accepted: 06/24/2021] [Indexed: 12/19/2022]
Abstract
PURPOSE This work was aimed at designing a deep-learning-based approach for MR image phase unwrapping to improve the robustness and efficiency of traditional methods. METHODS A deep learning network called PHU-NET was designed for MR phase unwrapping. In this network, a novel training data generation method was proposed to simulate the wrapping patterns in MR phase images. The wrapping boundary and wrapping counts were explicitly estimated and used for network training. The proposed method was quantitatively evaluated and compared to other methods using a number of simulated datasets with varying signal-to-noise ratio (SNR) and MR phase images from various parts of the human body. RESULTS The results showed that our method performed better in the simulated data even under an extremely low SNR. The proposed method had less residual wrapping in the images from various parts of human body and worked well in the presence of severe anatomical discontinuity. Our method was also advantageous in terms of computational efficiency compared to the traditional methods. CONCLUSION This work proposed a robust and computationally efficient MR phase unwrapping method based on a deep learning network, which has promising performance in applications using MR phase information.
Collapse
Affiliation(s)
- Hongyu Zhou
- Paul C. Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China.,University of Chinese Academy of Sciences, Beijing, China
| | - Chuanli Cheng
- Paul C. Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China.,University of Chinese Academy of Sciences, Beijing, China
| | - Hao Peng
- Paul C. Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China.,School of Artificial Intelligence and Automation, Huazhong University of Science and Technology, Wuhan, China
| | - Dong Liang
- Paul C. Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China.,University of Chinese Academy of Sciences, Beijing, China.,Key Laboratory for Magnetic Resonance and Multimodality Imaging of Guangdong Province, Shenzhen, China
| | - Xin Liu
- Paul C. Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China.,University of Chinese Academy of Sciences, Beijing, China.,Key Laboratory for Magnetic Resonance and Multimodality Imaging of Guangdong Province, Shenzhen, China
| | - Hairong Zheng
- Paul C. Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China.,University of Chinese Academy of Sciences, Beijing, China.,Key Laboratory for Magnetic Resonance and Multimodality Imaging of Guangdong Province, Shenzhen, China
| | - Chao Zou
- Paul C. Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China.,University of Chinese Academy of Sciences, Beijing, China.,Key Laboratory for Magnetic Resonance and Multimodality Imaging of Guangdong Province, Shenzhen, China
| |
Collapse
|
43
|
Chen Z, Chen Z. Computed inverse MRI (CIMRI) for intrinsic brain magnetic susceptibility mapping. Comput Biol Med 2021; 134:104498. [PMID: 34051451 DOI: 10.1016/j.compbiomed.2021.104498] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2021] [Revised: 04/30/2021] [Accepted: 05/12/2021] [Indexed: 11/24/2022]
Abstract
In magnetic resonance imaging (MRI), tissue magnetization in the main field B0 is a necessary preparation for magnetic resonance signal formation that imposes an inherent dipole effect on MRI signals, which predisposes an artifact on tissue MRI. In the MRI principle, T2*-weighted MRI can be described by a cascade of data transformations: from the source of tissue magnetic susceptibility (denoted by χ) to the output of complex-valued T2* image (in a magnitude and phase pair). Under the linear approximation of the T2* phase MRI, we can computationally reconstruct the source χ by quantitative susceptibility mapping (QSM), which is an inverse solution that is modeled by computed inverse MRI (CIMRI). For a brain function study using MRI (fMRI), we can reconstruct a timeseries of brain χ images to represent the intrinsic brain function activity called functional QSM (fQSM). This intrinsic depiction is defined as the removal of the artifactual dipole effect and other MRI-introduced distortions from phase data through inverse mapping. With one high-resolution QSM experiment and one group (20 subjects) low-resolution fQSM experiment, we show that the dipole effect manifests as ripples around vessels and a spatial split at a local activation blob and that the dipole effect could be removed by CIMRI. In the context of inverse imaging or undoing MRI transformations (including dipole convolution), we computationally achieve brain intrinsic structural depiction by QSM and intrinsic functional depiction by fQSM.
Collapse
Affiliation(s)
- Zeyuan Chen
- Department of Computer Sciences, University of California-Davis, Davis, CA, 95616, USA.
| | - Zikuan Chen
- Department of Diagnostic Radiology, City of Hope National Medical Center, Duarte, CA, 91010, USA; Zinv LLC, Albuquerque, NM, 87108, USA.
| |
Collapse
|
44
|
Dymerska B, Eckstein K, Bachrata B, Siow B, Trattnig S, Shmueli K, Robinson SD. Phase unwrapping with a rapid opensource minimum spanning tree algorithm (ROMEO). Magn Reson Med 2020; 85:2294-2308. [PMID: 33104278 PMCID: PMC7821134 DOI: 10.1002/mrm.28563] [Citation(s) in RCA: 43] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2020] [Revised: 09/24/2020] [Accepted: 09/30/2020] [Indexed: 01/12/2023]
Abstract
PURPOSE To develop a rapid and accurate MRI phase-unwrapping technique for challenging phase topographies encountered at high magnetic fields, around metal implants, or postoperative cavities, which is sufficiently fast to be applied to large-group studies including Quantitative Susceptibility Mapping and functional MRI (with phase-based distortion correction). METHODS The proposed path-following phase-unwrapping algorithm, ROMEO, estimates the coherence of the signal both in space-using MRI magnitude and phase information-and over time, assuming approximately linear temporal phase evolution. This information is combined to form a quality map that guides the unwrapping along a 3D path through the object using a computationally efficient minimum spanning tree algorithm. ROMEO was tested against the two most commonly used exact phase-unwrapping methods, PRELUDE and BEST PATH, in simulated topographies and at several field strengths: in 3T and 7T in vivo human head images and 9.4T ex vivo rat head images. RESULTS ROMEO was more reliable than PRELUDE and BEST PATH, yielding unwrapping results with excellent temporal stability for multi-echo or multi-time-point data. It does not require image masking and delivers results within seconds, even in large, highly wrapped multi-echo data sets (eg, 9 seconds for a 7T head data set with 31 echoes and a 208 × 208 × 96 matrix size). CONCLUSION Overall, ROMEO was both faster and more accurate than PRELUDE and BEST PATH, delivering exact results within seconds, which is well below typical image acquisition times, enabling potential on-console application.
Collapse
Affiliation(s)
- Barbara Dymerska
- Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| | - Korbinian Eckstein
- High Field MR Center, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
| | - Beata Bachrata
- High Field MR Center, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria.,Christian Doppler Laboratory for Clinical Molecular MR Imaging, Medical University of Vienna, Vienna, Austria
| | - Bernard Siow
- Magnetic Resonance Imaging, The Francis Crick Institute, London, United Kingdom
| | - Siegfried Trattnig
- High Field MR Center, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria.,Christian Doppler Laboratory for Clinical Molecular MR Imaging, Medical University of Vienna, Vienna, Austria
| | - Karin Shmueli
- Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| | - Simon Daniel Robinson
- High Field MR Center, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria.,Centre for Advanced Imaging, University of Queensland, Australia.,Department of Neurology, Medical University of Graz, Graz, Austria
| |
Collapse
|