1
|
Endres J, Weinmüller S, Dang HN, Zaiss M. Phase distribution graphs for fast, differentiable, and spatially encoded Bloch simulations of arbitrary MRI sequences. Magn Reson Med 2024; 92:1189-1204. [PMID: 38576164 DOI: 10.1002/mrm.30055] [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/26/2023] [Revised: 02/02/2024] [Accepted: 02/02/2024] [Indexed: 04/06/2024]
Abstract
PURPOSE An analytical approach to Bloch simulations for MRI sequences is introduced that enables time efficient calculations of signals free of Monte-Carlo noise, while providing full flexibility and differentiability in RF flip angles, RF phases, magnetic field gradients and time, as well as insights into image formation. THEORY AND METHODS We present an implementation of the extended phase graph (EPG) concept implemented in PyTorch, including an efficient state selection algorithm. This simulation is compared with an isochromat-based Bloch simulation with random isochromat distribution as well as with in vivo measurements using the Pulseq standard. Additionally, different sequences are tested and analyzed using this novel simulation approach. RESULTS Our simulation outperforms isochromat-based simulations in terms of computation time as well as signal quality, without exhibiting any kind of Monte-Carlo noise. The novel approach allows extracting useful information about the magnetization evolution not available otherwise. Our approach extends the common state-tensor-based EPG simulation approach for the contribution of dephased states including spatial encoding andT 2 ' $$ {T}_2^{\prime } $$ effects, and arbitrary timing. This allows calculation of echo shapes in addition to echo amplitudes only. Our implementation provides full differentiability in all input parameters allowing gradient descent optimization. Simulation of non-instantaneous pulses via hard-pulse approximation is left for future work, as the performance and accuracy characteristics are not yet analyzed. CONCLUSIONS Phase distribution graphs provide fast, differentiable, and spatially encoded Bloch simulations for most MRI sequences. It allows efficient simulation and optimization of arbitrary MRI sequences, which was previously only possible via high isochromat counts.
Collapse
Affiliation(s)
- Jonathan Endres
- Institute of Neuroradiology, University Hospital Erlangen, Friedrich-Alexander University Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Simon Weinmüller
- Institute of Neuroradiology, University Hospital Erlangen, Friedrich-Alexander University Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Hoai Nam Dang
- Institute of Neuroradiology, University Hospital Erlangen, Friedrich-Alexander University Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Moritz Zaiss
- Institute of Neuroradiology, University Hospital Erlangen, Friedrich-Alexander University Erlangen-Nürnberg (FAU), Erlangen, Germany
- Magnetic Resonance Center, Max-Planck-Institute, Tübingen, Germany
- Department Artificial Intelligence in Biomedical Engineering, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany
| |
Collapse
|
2
|
Mahmud SZ, Singh M, van Zijl P, Heo HY. Fast and motion-robust saturation transfer MRI with inherent B 0 correction using rosette trajectories and compressed sensing. Magn Reson Med 2024. [PMID: 39129199 DOI: 10.1002/mrm.30249] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2024] [Revised: 07/25/2024] [Accepted: 07/26/2024] [Indexed: 08/13/2024]
Abstract
PURPOSE To implement rosette readout trajectories with compressed sensing reconstruction for fast and motion-robust CEST and magnetization transfer contrast imaging with inherent correction of B0 inhomogeneity. METHODS A pulse sequence was developed for fast saturation transfer imaging using a stack of rosette trajectories with a higher sampling density near the k-space center. Each rosette lobe was segmented into two halves to generate dual-echo images. B0 inhomogeneities were estimated using the phase difference between the images and corrected subsequently. The rosette-based imaging was evaluated in comparison to a fully sampled Cartesian trajectory and demonstrated on CEST phantoms (creatine solutions and egg white) and healthy volunteers at 3 T. RESULTS Compared with the conventional Cartesian acquisition, compressed sensing reconstructed rosette images provided image quality with overall higher contrast-to-noise ratio and significantly faster readout time. Accurate B0 map estimation was achieved from the rosette acquisition with a negligible bias of 0.01 Hz between the rosette and dual-echo Cartesian gradient echo B0 maps, using the latter as ground truth. The water-saturation spectra (Z-spectra) and amide proton transfer weighted signals obtained from the rosette-based sequence were well preserved compared with the fully sampled data, both in the phantom and human studies. CONCLUSIONS Fast, motion-robust, and inherent B0-corrected CEST and magnetization transfer contrast imaging using rosette trajectories could improve subject comfort and compliance, contrast-to-noise ratio, and provide inherent B0 homogeneity information. This work is expected to significantly accelerate the translation of CEST-MRI into a robust, clinically viable approach.
Collapse
Affiliation(s)
- Sultan Z Mahmud
- Department of Radiology, Johns Hopkins University, Baltimore, Maryland, USA
| | - Munendra Singh
- Department of Radiology, Johns Hopkins University, Baltimore, Maryland, USA
| | - Peter van Zijl
- Department of Radiology, Johns Hopkins University, Baltimore, Maryland, USA
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, USA
| | - Hye-Young Heo
- Department of Radiology, Johns Hopkins University, Baltimore, Maryland, USA
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, USA
| |
Collapse
|
3
|
Iyyakkunnel S, Weigel M, Bieri O. Fast bias-corrected conductivity mapping using stimulated echoes. MAGMA (NEW YORK, N.Y.) 2024:10.1007/s10334-024-01194-3. [PMID: 39105952 DOI: 10.1007/s10334-024-01194-3] [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/23/2024] [Revised: 07/15/2024] [Accepted: 07/22/2024] [Indexed: 08/07/2024]
Abstract
OBJECTIVE To demonstrate the potential of a double angle stimulated echo (DA-STE) method for fast and accurate "full" homogeneous Helmholtz-based electrical properties tomography using a simultaneous B 1 + magnitude and transceive phase measurement. METHODS The combination of a spin and stimulated echo can be used to yield an estimate of both B 1 + magnitude and transceive phase and thus provides the means for "full" EPT reconstruction. An interleaved 2D acquisition scheme is used for rapid acquisition. The method was validated in a saline phantom and compared to a double angle method based on two single gradient echo acquisitions (GRE-DAM). The method was evaluated in the brain of a healthy volunteer. RESULTS The B 1 + magnitude obtained with DA-STE showed excellent agreement with the GRE-DAM method. Conductivity values based on the "full" EPT reconstruction also agreed well with the expectations in the saline phantom. In the brain, the method delivered conductivity values close to literature values. DISCUSSION The method allows the use of the "full" Helmholtz-based EPT reconstruction without the need for additional measurements. As a result, quantitative conductivity values are improved compared to phase-based EPT reconstructions. DA-STE is a fast complex- B 1 + mapping technique that could render EPT clinically relevant at 3 T.
Collapse
Affiliation(s)
- Santhosh Iyyakkunnel
- Division of Radiological Physics, Department of Radiology , University Hospital Basel, Basel, Switzerland.
- Department of Biomedical Engineering, University of Basel, Basel, Switzerland.
| | - Matthias Weigel
- Division of Radiological Physics, Department of Radiology , University Hospital Basel, Basel, Switzerland
- Department of Biomedical Engineering, University of Basel, Basel, Switzerland
- Translational Imaging in Neurology (ThINk) Basel, Department of Biomedical Engineering, Faculty of Medicine, University Hospital Basel and University of Basel, Basel, Switzerland
- Neurologic Clinic and Policlinic, MS Center and Research Center for Clinical Neuroimmunology and Neuroscience Basel (RC2NB), University Hospital Basel and University of Basel, Basel, Switzerland
| | - Oliver Bieri
- Division of Radiological Physics, Department of Radiology , University Hospital Basel, Basel, Switzerland
- Department of Biomedical Engineering, University of Basel, Basel, Switzerland
| |
Collapse
|
4
|
Bhattacharya S, Price AN, Uus A, Sousa HS, Marenzana M, Colford K, Murkin P, Lee M, Cordero-Grande L, Teixeira RPAG, Malik SJ, Deprez M. In vivo T2 measurements of the fetal brain using single-shot fast spin echo sequences. Magn Reson Med 2024; 92:715-729. [PMID: 38623934 DOI: 10.1002/mrm.30094] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2023] [Revised: 02/18/2024] [Accepted: 03/08/2024] [Indexed: 04/17/2024]
Abstract
PURPOSE We propose a quantitative framework for motion-corrected T2 fetal brain measurements in vivo and validate the single-shot fast spin echo (SS-FSE) sequence to perform these measurements. METHODS Stacks of two-dimensional SS-FSE slices are acquired with different echo times (TE) and motion-corrected with slice-to-volume reconstruction (SVR). The quantitative T2 maps are obtained by a fit to a dictionary of simulated signals. The sequence is selected using simulated experiments on a numerical phantom and validated on a physical phantom scanned on a 1.5T system. In vivo quantitative T2 maps are obtained for five fetuses with gestational ages (GA) 21-35 weeks on the same 1.5T system. RESULTS The simulated experiments suggested that a TE of 400 ms combined with the clinically utilized TEs of 80 and 180 ms were most suitable for T2 measurements in the fetal brain. The validation on the physical phantom confirmed that the SS-FSE T2 measurements match the gold standard multi-echo spin echo measurements. We measured average T2s of around 200 and 280 ms in the fetal brain grey and white matter, respectively. This was slightly higher than fetal T2* and the neonatal T2 obtained from previous studies. CONCLUSION The motion-corrected SS-FSE acquisitions with varying TEs offer a promising practical framework for quantitative T2 measurements of the moving fetus.
Collapse
Affiliation(s)
- Suryava Bhattacharya
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Anthony N Price
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- Guy's and St Thomas' NHS Foundation Trust, London, UK
- Centre for the Developing Brain, King's College London, London, UK
| | - Alena Uus
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- Centre for the Developing Brain, King's College London, London, UK
| | - Helena S Sousa
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | | | - Kathleen Colford
- Centre for the Developing Brain, King's College London, London, UK
| | - Peter Murkin
- Guy's and St Thomas' NHS Foundation Trust, London, UK
- Centre for the Developing Brain, King's College London, London, UK
| | - Maggie Lee
- Guy's and St Thomas' NHS Foundation Trust, London, UK
- Centre for the Developing Brain, King's College London, London, UK
| | - Lucilio Cordero-Grande
- Biomedical Image Technologies, ETSI Telecomunicración, Universidad Politécnica de Madrid and CIBER-BBN, Madrid, Spain
| | - Rui Pedro A G Teixeira
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- Centre for the Developing Brain, King's College London, London, UK
| | - Shaihan J Malik
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- Centre for the Developing Brain, King's College London, London, UK
| | - Maria Deprez
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- Centre for the Developing Brain, King's College London, London, UK
| |
Collapse
|
5
|
van der Heide O, van den Berg CAT, Sbrizzi A. GPU-accelerated Bloch simulations and MR-STAT reconstructions using the Julia programming language. Magn Reson Med 2024; 92:618-630. [PMID: 38441315 DOI: 10.1002/mrm.30074] [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/30/2023] [Revised: 02/15/2024] [Accepted: 02/16/2024] [Indexed: 06/02/2024]
Abstract
PURPOSE MR-STAT is a relatively new multiparametric quantitative MRI technique in which quantitative paramater maps are obtained by solving a large-scale nonlinear optimization problem. Managing reconstruction times is one of the main challenges of MR-STAT. In this work we leverage GPU hardware to reduce MR-STAT reconstruction times. A highly optimized, GPU-compatible Bloch simulation toolbox is developed as part of this work that can be utilized for other quantitative MRI techniques as well. METHODS The Julia programming language was used to develop a flexible yet highly performant and GPU-compatible Bloch simulation toolbox called BlochSimulators.jl. The runtime performance of the toolbox is benchmarked against other Bloch simulation toolboxes. Furthermore, a (partially matrix-free) modification of a previously presented (matrix-free) MR-STAT reconstruction algorithm is proposed and implemented using the Julia language on GPU hardware. The proposed algorithm is combined with BlochSimulators.jl and the resulting MR-STAT reconstruction times on GPU hardware are compared to previously presented MR-STAT reconstruction times. RESULTS The BlochSimulators.jl package demonstrates superior runtime performance on both CPU and GPU hardware when compared to other existing Bloch simulation toolboxes. The GPU-accelerated partially matrix-free MR-STAT reconstruction algorithm, which relies on BlochSimulators.jl, allows for reconstructions of 68 seconds per two-dimensional (2D slice). CONCLUSION By combining the proposed Bloch simulation toolbox and the partially matrix-free reconstruction algorithm, 2D MR-STAT reconstructions can be performed in the order of one minute on a modern GPU card. The Bloch simulation toolbox can be utilized for other quantitative MRI techniques as well, for example for online dictionary generation for MR Fingerprinting.
Collapse
Affiliation(s)
- Oscar van der Heide
- Computational Imaging Group for MR Diagnostics and Therapy, Center for Image Sciences, University Medical Center Utrecht, Utrecht, The Netherlands
- Department of Radiotherapy, Division of Imaging and Oncology, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Cornelis A T van den Berg
- Computational Imaging Group for MR Diagnostics and Therapy, Center for Image Sciences, University Medical Center Utrecht, Utrecht, The Netherlands
- Department of Radiotherapy, Division of Imaging and Oncology, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Alessandro Sbrizzi
- Computational Imaging Group for MR Diagnostics and Therapy, Center for Image Sciences, University Medical Center Utrecht, Utrecht, The Netherlands
- Department of Radiotherapy, Division of Imaging and Oncology, University Medical Center Utrecht, Utrecht, The Netherlands
| |
Collapse
|
6
|
Barzaghi L, Paoletti M, Monforte M, Bortolani S, Bonizzoni C, Thorsten F, Bergsland N, Santini F, Deligianni X, Tasca G, Ballante E, Figini S, Ricci E, Pichiecchio A. Muscle diffusion tensor imaging in facioscapulohumeral muscular dystrophy. Muscle Nerve 2024; 70:248-256. [PMID: 38873946 DOI: 10.1002/mus.28179] [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/18/2023] [Revised: 05/16/2024] [Accepted: 05/26/2024] [Indexed: 06/15/2024]
Abstract
INTRODUCTION/AIMS Muscle diffusion tensor imaging has not yet been explored in facioscapulohumeral muscular dystrophy (FSHD). We assessed diffusivity parameters in FSHD subjects compared with healthy controls (HCs), with regard to their ability to precede any fat replacement or edema. METHODS Fat fraction (FF), water T2 (wT2), mean, radial, axial diffusivity (MD, RD, AD), and fractional anisotropy (FA) of thigh muscles were calculated in 10 FSHD subjects and 15 HCs. All parameters were compared between FSHD and controls, also exploring their gradient along the main axis of the muscle. Diffusivity parameters were tested in a subgroup analysis as predictors of disease involvement in muscle compartments with different degrees of FF and wT2 and were also correlated with clinical severity scores. RESULTS We found that MD, RD, and AD were significantly lower in FSHD subjects than in controls, whereas we failed to find a difference for FA. In contrast, we found a significant positive correlation between FF and FA and a negative correlation between MD, RD, and AD and FF. No correlation was found with wT2. In our subgroup analysis we found that muscle compartments with no significant fat replacement or edema (FF < 10% and wT2 < 41 ms) showed a reduced AD and FA compared with controls. Less involved compartments showed different diffusivity parameters than more involved compartments. DISCUSSION Our exploratory study was able to demonstrate diffusivity parameter abnormalities even in muscles with no significant fat replacement or edema. Larger cohorts are needed to confirm these preliminary findings.
Collapse
Affiliation(s)
- Leonardo Barzaghi
- Department of Mathematics, University of Pavia, Pavia, Italy
- Advanced Imaging and Artificial Intelligence Center, Neuroradiology Department, IRCCS Mondino Foundation, Pavia, Italy
- INFN, Group of Pavia, Pavia, Italy
| | - Matteo Paoletti
- Advanced Imaging and Artificial Intelligence Center, Neuroradiology Department, IRCCS Mondino Foundation, Pavia, Italy
| | - Mauro Monforte
- UOC di Neurologia, Fondazione Policlinico Universitario A. Gemelli IRCCS, Rome, Italy
| | - Sara Bortolani
- UOC di Neurologia, Fondazione Policlinico Universitario A. Gemelli IRCCS, Rome, Italy
| | - Chiara Bonizzoni
- Advanced Imaging and Artificial Intelligence Center, Neuroradiology Department, IRCCS Mondino Foundation, Pavia, Italy
| | | | - Niels Bergsland
- Department of Neurology, Jacobs School of Medicine and Biomedical Sciences, Buffalo Neuroimaging Analysis Center, University of Buffalo, The State University of New York, Buffalo, New York, USA
- IRCCS, Fondazione Don Carlo Gnocchi ONLUS, Milan, Italy
| | - Francesco Santini
- Department of Radiology, University Hospital Basel, Basel, Switzerland
- Basel Muscle MRI, Department of Biomedical Engineering, University of Basel, Basel, Switzerland
| | - Xeni Deligianni
- Department of Radiology, University Hospital Basel, Basel, Switzerland
- Basel Muscle MRI, Department of Biomedical Engineering, University of Basel, Basel, Switzerland
| | - Giorgio Tasca
- UOC di Neurologia, Fondazione Policlinico Universitario A. Gemelli IRCCS, Rome, Italy
- John Walton Muscular Dystrophy Research Centre, Newcastle University and Newcastle Hospitals NHS Foundation Trusts, Newcastle upon Tyne, UK
| | - Elena Ballante
- Department of Political and Social Sciences, University of Pavia, Pavia, Italy
- BioData Science Center, IRCCS Mondino Foundation, Pavia, Italy
| | - Silvia Figini
- Department of Political and Social Sciences, University of Pavia, Pavia, Italy
- BioData Science Center, IRCCS Mondino Foundation, Pavia, Italy
| | - Enzo Ricci
- UOC di Neurologia, Fondazione Policlinico Universitario A. Gemelli IRCCS, Rome, Italy
| | - Anna Pichiecchio
- Advanced Imaging and Artificial Intelligence Center, Neuroradiology Department, IRCCS Mondino Foundation, Pavia, Italy
- Department of Brain and Behavioral Sciences, University of Pavia, Pavia, Italy
| |
Collapse
|
7
|
Lowen D, Pracht ED, Gras V, Massire A, Mauconduit F, Stoecker T, Boulant N. Design of calibration-free RF pulses for T 2 $$ {}_2 $$ -weighted single-slab 3D turbo-spin-echo sequences at 7T utilizing parallel transmission. Magn Reson Med 2024. [PMID: 39054786 DOI: 10.1002/mrm.30212] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2024] [Revised: 06/03/2024] [Accepted: 06/19/2024] [Indexed: 07/27/2024]
Abstract
PURPOSE T 2 $$ {}_2 $$ -weighted turbo-spin-echo (TSE) sequences are a fundamental technique in brain imaging but suffer from field inhomogeneities at ultra-high fields. Several methods have been proposed to mitigate the problem, but were limited so far to nonselective three-dimensional (3D) measurements, making short acquisitions difficult to achieve when targeting very high resolution images, or needed additional calibration procedures, thus complicating their application. METHODS Slab-selective excitation pulses were designed for flexible placement utilizing the concept of k T $$ {}_T $$ -spokes. Phase-coherent refocusing universal pulses were subsequently optimized with the Gradient Ascent Pulse Engineering algorithm and tested in vivo for improved signal homogeneity. RESULTS Implemented within a 3D variable flip angle TSE sequence, these pulses led to a signal-to-noise ratio (SNR) improvement ranging from 10% to 30% compared to a two-dimensional (2D) T2w TSE sequence employingB 1 + $$ {\mathrm{B}}_1^{+} $$ -shimmed pulses.B 1 + $$ {\mathrm{B}}_1^{+} $$ field inhomogeneities could be mitigated and artifacts fromB 0 $$ {\mathrm{B}}_0 $$ deviations reduced. The concept of universal pulses was successfully applied. CONCLUSION We present a pulse design method which provides a set of calibration-free universal pulses (UPs) for slab-selective excitation and phase-coherent refocusing in slab-selective TSE sequences.
Collapse
Affiliation(s)
- Daniel Lowen
- MR Physics, German Center for Neurodegenerative Diseases (DZNE), Bonn, Germany
| | - Eberhard D Pracht
- MR Physics, German Center for Neurodegenerative Diseases (DZNE), Bonn, Germany
| | - Vincent Gras
- Commissariat à l'Energie Atomique, CNRS, NeuroSpin, BAOBAB, Université Paris-Saclay, Gif sur Yvette, France
| | | | - Franck Mauconduit
- Commissariat à l'Energie Atomique, CNRS, NeuroSpin, BAOBAB, Université Paris-Saclay, Gif sur Yvette, France
| | - Tony Stoecker
- MR Physics, German Center for Neurodegenerative Diseases (DZNE), Bonn, Germany
- Department of Physics & Astronomy, University of Bonn, Bonn, Germany
| | - Nicolas Boulant
- Commissariat à l'Energie Atomique, CNRS, NeuroSpin, BAOBAB, Université Paris-Saclay, Gif sur Yvette, France
| |
Collapse
|
8
|
Li S, Wang L, Priest AN, Horvat-Menih I, Mendichovszky IA, Gallagher FA, Wang H, Li H. Highly accelerated parameter mapping using model-based alternating reconstruction coupling fitting. Phys Med Biol 2024; 69:145014. [PMID: 38917824 DOI: 10.1088/1361-6560/ad5bb8] [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: 03/01/2024] [Accepted: 06/25/2024] [Indexed: 06/27/2024]
Abstract
Objective.A model-based alternating reconstruction coupling fitting, termed Model-based Alternating Reconstruction COupling fitting (MARCO), is proposed for accurate and fast magnetic resonance parameter mapping.Approach.MARCO utilizes the signal model as a regularization by minimizing the bias between the image series and the signal produced by the suitable signal model based on iteratively updated parameter maps when reconstructing. The technique can incorporate prior knowledge of both image series and parameters by adding sparsity constraints. The optimization problem is decomposed into three subproblems and solved through three alternating steps involving reconstruction and nonlinear least-square fitting, which can produce both contrast-weighted images and parameter maps simultaneously.Main results.The algorithm is applied toT2mapping with extended phase graph algorithm integrated and validated on undersampled multi-echo spin-echo data from both phantom and in vivo sources. Compared with traditional compressed sensing and model-based methods, the proposed approach yields more accurateT2maps with more details at high acceleration factors.Significance.The proposed method provides a basic framework for quantitative MR relaxometry, theoretically applicable to all quantitative MR relaxometry. It has the potential to improve the diagnostic utility of quantitative imaging techniques.
Collapse
Affiliation(s)
- Shaohang Li
- Institute of Science and Technology for Brain-inspired Intelligence, Fudan University, Shanghai 200433, People's Republic of China
| | - Lili Wang
- Institute of Science and Technology for Brain-inspired Intelligence, Fudan University, Shanghai 200433, People's Republic of China
| | - Andrew N Priest
- Department of Radiology, University of Cambridge, Cambridge CB2 0QQ, United Kingdom
- Department of Radiology, Addenbrooke's Hospital, Cambridge University Hospitals NHS Foundation Trust, Cambridge CB2 0QQ, United Kingdom
| | - Ines Horvat-Menih
- Department of Radiology, University of Cambridge, Cambridge CB2 0QQ, United Kingdom
| | - Iosif A Mendichovszky
- Department of Radiology, University of Cambridge, Cambridge CB2 0QQ, United Kingdom
- Department of Radiology, Addenbrooke's Hospital, Cambridge University Hospitals NHS Foundation Trust, Cambridge CB2 0QQ, United Kingdom
| | - Ferdia A Gallagher
- Department of Radiology, University of Cambridge, Cambridge CB2 0QQ, United Kingdom
- Department of Radiology, Addenbrooke's Hospital, Cambridge University Hospitals NHS Foundation Trust, Cambridge CB2 0QQ, United Kingdom
| | - He Wang
- Institute of Science and Technology for Brain-inspired Intelligence, Fudan University, Shanghai 200433, People's Republic of China
| | - Hao Li
- Institute of Science and Technology for Brain-inspired Intelligence, Fudan University, Shanghai 200433, People's Republic of China
- Department of Radiology, University of Cambridge, Cambridge CB2 0QQ, United Kingdom
| |
Collapse
|
9
|
von Deuster C, Nanz D. Enhancing fluid signal in driven-equilibrium short-TI inversion-recovery imaging with short TR times: A feasibility study. Magn Reson Med 2024. [PMID: 38987979 DOI: 10.1002/mrm.30215] [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: 11/29/2023] [Revised: 06/07/2024] [Accepted: 06/22/2024] [Indexed: 07/12/2024]
Abstract
PURPOSE Fluid-sensitive turbo spin echo (TSE) MRI with short-TI inversion-recovery preparation for fat suppression (STIR) plays a critical role in the diagnostics of the musculoskeletal system (e.g., close to metal implants). Potential advantages of 3D acquisitions, however, are difficult to exploit due to long acquisition times. Shortening the TR incurs a signal loss, and a driven-equilibrium (DE) extension reduces fluid signal even further. METHODS The phase of the flip-back pulse was changed by 180° relative to the conventional implementation (i.e., 90° along the positive x-axis (90°x) instead of -90°x). After signal modeling and numerical simulations, the modification was implemented in STIR-TSE sequences and tested on a clinical 3T system. Imaging was performed in the lumbar spine, and long-TR images without DE were acquired as reference. CSF SNR and fluid-muscle contrast were measured and compared between the sequences. Imaging was repeated in a metal implant phantom. RESULTS A shortening of TR by 43%-57% reduced the CSF SNR by 39%-59%. A conventional DE module further reduced SNR to 26%-40%, whereas the modified DE recovered SNR to 59%-108% compared with the long-TR acquisitions. Fluid-tissue contrast was increased by about 340% with the modified DE module compared with the conventional extension. Similar results were obtained in implant measurements. CONCLUSIONS The proposed DE element for TSE-STIR sequences has the potential to accelerate the acquisition of fluid-sensitive images. DE-STIR may work most efficiently for 3D acquisitions, in which no temporo-spatial interleaving of inversion and imaging pulses is possible.
Collapse
Affiliation(s)
- Constantin von Deuster
- Advanced Clinical Imaging Technology, Siemens Healthineers International AG, Zurich, Switzerland
- Swiss Center for Musculoskeletal Imaging, Balgrist Campus AG, Zurich, Switzerland
| | - Daniel Nanz
- Swiss Center for Musculoskeletal Imaging, Balgrist Campus AG, Zurich, Switzerland
- Medical Faculty, University of Zurich, Zurich, Switzerland
| |
Collapse
|
10
|
de Buck MHS, Hess AT, Jezzard P. Simulation-based optimization and experimental comparison of intracranial T2-weighted DANTE-SPACE vessel wall imaging at 3T and 7T. Magn Reson Med 2024. [PMID: 38970460 DOI: 10.1002/mrm.30203] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Revised: 04/30/2024] [Accepted: 06/12/2024] [Indexed: 07/08/2024]
Abstract
PURPOSE T2-weighted DANTE-SPACE (Delay Alternating with Nutation for Tailored Excitation - Sampling Perfection with Application optimized Contrasts using different flip angle Evolution) sequences facilitate non-invasive intracranial vessel wall imaging at 7T through simultaneous suppression of blood and CSF. However, the achieved vessel wall delineation depends closely on the selected sequence parameters, and little information is available about the performance of the sequence using more widely available 3T MRI. Therefore, in this paper a comprehensive DANTE-SPACE simulation framework is used for the optimization and quantitative comparison of T2-weighted DANTE-SPACE at both 7T and 3T. METHODS Simulations are used to propose optimized sequence parameters at both 3T and 7T. At 7T, an additional protocol which uses a parallel transmission (pTx) shim during the DANTE preparation for improved suppression of inflowing blood is also proposed. Data at both field strengths using optimized and literature protocols are acquired and quantitatively compared in six healthy volunteers. RESULTS At 7T, more vessel wall signal can be retained while still achieving sufficient CSF suppression by using fewer DANTE pulses than described in previous implementations. The use of a pTx shim during DANTE at 7T provides a modest further improvement to the inner vessel wall delineation. At 3T, aggressive DANTE preparation is required to achieve CSF suppression, resulting in reduced vessel wall signal. As a result, the achievable vessel wall definition at 3T is around half that of 7T. CONCLUSION Simulation-based optimization of DANTE parameters facilitates improved T2-weighted DANTE-SPACE contrasts at 7T. The improved vessel definition of T2-weighted DANTE-SPACE at 7T makes DANTE preparation more suitable for T2-weighted VWI at 7T than at 3T.
Collapse
Affiliation(s)
- Matthijs H S de Buck
- Wellcome Centre for Integrative Neuroimaging, FMRIB Division, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, UK
- Spinoza Centre for Neuroimaging, Netherlands Institute for Neuroscience, Royal Netherlands Academy for Arts and Sciences (KNAW), Amsterdam, The Netherlands
| | - Aaron T Hess
- Wellcome Centre for Integrative Neuroimaging, FMRIB Division, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, UK
| | - Peter Jezzard
- Wellcome Centre for Integrative Neuroimaging, FMRIB Division, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, UK
| |
Collapse
|
11
|
Lasič S, Chakwizira A, Lundell H, Westin CF, Nilsson M. Tuned exchange imaging: Can the filter exchange imaging pulse sequence be adapted for applications with thin slices and restricted diffusion? NMR IN BIOMEDICINE 2024:e5208. [PMID: 38961745 DOI: 10.1002/nbm.5208] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/05/2024] [Revised: 06/04/2024] [Accepted: 06/06/2024] [Indexed: 07/05/2024]
Abstract
Filter exchange imaging (FEXI) is a double diffusion-encoding (DDE) sequence that is specifically sensitive to exchange between sites with different apparent diffusivities. FEXI uses a diffusion-encoding filtering block followed by a detection block at varying mixing times to map the exchange rate. Long mixing times enhance the sensitivity to exchange, but they pose challenges for imaging applications that require a stimulated echo sequence with crusher gradients. Thin imaging slices require strong crushers, which can introduce significant diffusion weighting and bias exchange rate estimates. Here, we treat the crushers as an additional encoding block and consider FEXI as a triple diffusion-encoding sequence. This allows the bias to be corrected in the case of multi-Gaussian diffusion, but not easily in the presence of restricted diffusion. Our approach addresses challenges in the presence of restricted diffusion and relies on the ability to independently gauge sensitivities to exchange and restricted diffusion for arbitrary gradient waveforms. It follows two principles: (i) the effects of crushers are included in the forward model using signal cumulant expansion; and (ii) timing parameters of diffusion gradients in filter and detection blocks are adjusted to maintain the same level of restriction encoding regardless of the mixing time. This results in the tuned exchange imaging (TEXI) protocol. The accuracy of exchange mapping with TEXI was assessed through Monte Carlo simulations in spheres of identical sizes and gamma-distributed sizes, and in parallel hexagonally packed cylinders. The simulations demonstrate that TEXI provides consistent exchange rates regardless of slice thickness and restriction size, even with strong crushers. However, the accuracy depends on b-values, mixing times, and restriction geometry. The constraints and limitations of TEXI are discussed, including suggestions for protocol adaptations. Further studies are needed to optimize the precision of TEXI and assess the approach experimentally in realistic, heterogeneous substrates.
Collapse
Affiliation(s)
- Samo Lasič
- Department of Diagnostic Radiology, Lund University, Lund, Sweden
- Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital - Amager and Hvidovre, Copenhagen, Denmark
| | - Arthur Chakwizira
- Department of Medical Radiation Physics, Lund, Lund University, Lund, Sweden
| | - Henrik Lundell
- Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital - Amager and Hvidovre, Copenhagen, Denmark
- MR Section, DTU Health Tech, Technical University of Denmark, Lyngby, Denmark
| | - Carl-Fredrik Westin
- Department of Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts, USA
| | - Markus Nilsson
- Department of Clinical Sciences Lund, Radiology, Lund University, Lund, Sweden
| |
Collapse
|
12
|
Wong YL, Li T, Liu C, Lee HFV, Cheung LYA, Hui ESK, Cao P, Cai J. Reconstruction of multi-phase parametric maps in 4D-magnetic resonance fingerprinting (4D-MRF) by optimization of local T1 and T2 sensitivities. Med Phys 2024; 51:4721-4735. [PMID: 38386904 DOI: 10.1002/mp.17001] [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: 03/21/2023] [Revised: 01/18/2024] [Accepted: 01/23/2024] [Indexed: 02/24/2024] Open
Abstract
BACKGROUND Time-resolved magnetic resonance fingerprinting (MRF), or 4D-MRF, has been demonstrated its feasibility in motion management in radiotherapy (RT). However, the prohibitive long acquisition time is one of challenges of the clinical implementation of 4D-MRF. The shortening of acquisition time causes data insufficiency in each respiratory phase, leading to poor accuracies and consistencies of the predicted tissues' properties of each phase. PURPOSE To develop a technique for the reconstruction of multi-phase parametric maps in four-dimensional magnetic resonance fingerprinting (4D-MRF) through the optimization of local T1 and T2 sensitivities. METHODS The proposed technique employed an iterative optimization to tailor the data arrangement of each phase by manipulation of inter-phase frames, such that the T1 and T2 sensitivities, which were quantified by the modified Minkowski distance, of the truncated signal evolution curve was maximized. The multi-phase signal evolution curves were modified by sliding window reconstruction and inter-phase frame sharing (SWIFS). Motion correction (MC) and dot product matching were sequentially performed on the modified signal evolution and dictionary to reconstruct the multi-parametric maps. The proposed technique was evaluated by numerical simulations using the extended cardiac-torso (XCAT) phantom with regular and irregular breathing patterns, and by in vivo MRF data of three health volunteers and six liver cancer patients acquired at a 3.0 T scanner. RESULTS In simulation study, the proposed SWIFS approach achieved the overall mean absolute percentage error (MAPE) of 8.62% ± 1.59% and 16.2% ± 3.88% for the eight-phases T1 and T2 maps, respectively, in the sagittal view with irregular breathing patterns. In contrast, the overall MAPE of T1 and T2 maps generated by the conventional approach with multiple MRF repetitions were 22.1% ± 11.0% and 30.8% ± 14.9%, respectively. For in-vivo study, the predicted mean T1 and T2 of liver by the proposed SWIFS approach were 795 ms ± 38.9 ms and 58.3 ms ± 11.7 ms, respectively. CONCLUSIONS Both simulation and in vivo results showed that the approach empowered by T1 and T2 sensitivities optimization and sliding window under the shortened acquisition of MRF had superior performance in the estimation of multi-phase T1 and T2 maps as compared to the conventional approach with oversampling of MRF data.
Collapse
Affiliation(s)
- Yat Lam Wong
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University, Hong Kong, China
- Department of Clinical Oncology, Queen Mary Hospital, Hong Kong, China
| | - Tian Li
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University, Hong Kong, China
| | - Chenyang Liu
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University, Hong Kong, China
| | - Ho-Fun Victor Lee
- Department of Clinical Oncology, The University of Hong Kong, Hong Kong, China
| | - Lai-Yin Andy Cheung
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University, Hong Kong, China
- Department of Clinical Oncology, Oncology Center, St. Paul's Hospital, Hong Kong, China
| | - Edward Sai Kam Hui
- Department of Imaging and Interventional Radiology, The Chinese University of Hong Kong, Hong Kong, China
| | - Peng Cao
- Department of Diagnostic Radiology, The University of Hong Kong, Hong Kong, China
| | - Jing Cai
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University, Hong Kong, China
| |
Collapse
|
13
|
de Buck MHS, Jezzard P, Hess AT. An extended phase graph-based framework for DANTE-SPACE simulations including physiological, temporal, and spatial variations. Magn Reson Med 2024; 92:332-345. [PMID: 38469983 DOI: 10.1002/mrm.30071] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2023] [Revised: 01/18/2024] [Accepted: 02/09/2024] [Indexed: 03/13/2024]
Abstract
PURPOSE The delay alternating with nutation for tailored excitation (DANTE)-sampling perfection with application-optimized contrasts (SPACE) sequence facilitates 3D intracranial vessel wall imaging with simultaneous suppression of blood and CSF. However, the achieved image contrast depends closely on the selected sequence parameters, and the clinical use of the sequence is limited in vivo by observed signal variations in the vessel wall, CSF, and blood. This paper introduces a comprehensive DANTE-SPACE simulation framework, with the aim of providing a better understanding of the underlying contrast mechanisms and facilitating improved parameter selection and contrast optimization. METHODS An extended phase graph formalism was developed for efficient spin ensemble simulation of the DANTE-SPACE sequence. Physiological processes such as pulsatile flow velocity variation, varying flow directions, intravoxel velocity variation, diffusion, andB 1 + $$ {\mathrm{B}}_1^{+} $$ effects were included in the framework to represent the mechanisms behind the achieved signal levels accurately. RESULTS Intravoxel velocity variation improved temporal stability and robustness against small velocity changes. Time-varying pulsatile velocity variation affected CSF simulations, introducing periods of near-zero velocity and partial rephasing. Inclusion of diffusion effects was found to substantially reduce the CSF signal. Blood flow trajectory variations had minor effects, butB 1 + $$ {\mathrm{B}}_1^{+} $$ differences along the trajectory reduced DANTE efficiency in low-B 1 + $$ {\mathrm{B}}_1^{+} $$ areas. Introducing low-velocity pulsatility of both CSF and vessel wall helped explain the in vivo observed signal heterogeneity in both tissue types. CONCLUSION The presented simulation framework facilitates a more comprehensive optimization of DANTE-SPACE sequence parameters. Furthermore, the simulation framework helps to explain observed contrasts in acquired data.
Collapse
Affiliation(s)
- Matthijs H S de Buck
- Wellcome Centre for Integrative Neuroimaging, FMRIB Division, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, UK
- Spinoza Centre for Neuroimaging, Amsterdam, The Netherlands
| | - Peter Jezzard
- Wellcome Centre for Integrative Neuroimaging, FMRIB Division, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, UK
| | - Aaron T Hess
- Wellcome Centre for Integrative Neuroimaging, FMRIB Division, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, UK
| |
Collapse
|
14
|
Yablonski M, Zhou Z, Cao X, Schauman S, Liao C, Setsompop K, Yeatman JD. Fast and reliable quantitative measures of white matter development with magnetic resonance fingerprinting. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.06.26.600735. [PMID: 38979185 PMCID: PMC11230456 DOI: 10.1101/2024.06.26.600735] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
Developmental cognitive neuroscience aims to shed light on evolving relationships between brain structure and cognitive development. To this end, quantitative methods that reliably measure individual differences in brain tissue properties are fundamental. Standard qualitative MRI sequences are influenced by scan parameters and hardware-related biases, and also lack physical units, making the analysis of individual differences problematic. In contrast, quantitative MRI can measure physical properties of the tissue but with the cost of long scan durations and sensitivity to motion. This poses a critical limitation for studying young children. Here, we examine the reliability and validity of an efficient quantitative multiparameter mapping method - Magnetic Resonance Fingerprinting (MRF) - in children scanned longitudinally. We focus on T1 values in white matter, since quantitative T1 values are known to primarily reflect myelin content, a key factor in brain development. Forty-nine children aged 8-13y (mean 10.3y ±1.4) completed two scanning sessions 2-4 months apart. In each session, two 2-minute 3D-MRF scans at 1mm isotropic resolution were collected to evaluate the effect of scan duration on image quality and scan-rescan reliability. A separate calibration scan was used to measure B0 inhomogeneity and correct for bias. We examined the impact of scan time and B0 inhomogeneity correction on scan-rescan reliability of values in white matter, by comparing single 2-min and combined two 2-min scans, with and without B0-correction. Whole-brain voxel-based reliability analysis showed that combining two 2-min MRF scans improved reliability (pearson's r=0.87) compared with a single 2-min scan (r=0.84), while B0-correction had no effect on reliability in white matter (r=0.86 and 0.83 4-min vs 2-min). Using diffusion tractography, we delineated MRF-derived T1 profiles along major white matter fiber tracts and found similar or higher reliability for T1 from MRF compared to diffusion parameters (based on a 10-minute dMRI scan). Lastly, we found that T1 values in multiple white matter tracts were significantly correlated with age. In sum, MRF-derived T1 values were highly reliable in a longitudinal sample of children and replicated known age effects. Reliability in white matter was improved by longer scan duration but was not affected by B0-correction, making it a quick and straightforward scan to collect. We propose that MRF provides a promising avenue for acquiring quantitative brain metrics in children and patient populations where scan time and motion are of particular concern.
Collapse
|
15
|
Kobayashi N. Optimization of flip angle and radiofrequency pulse phase to maximize steady-state magnetization in three-dimensional missing pulse steady-state free precession. NMR IN BIOMEDICINE 2024; 37:e5112. [PMID: 38299770 PMCID: PMC11078623 DOI: 10.1002/nbm.5112] [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: 07/18/2023] [Revised: 12/07/2023] [Accepted: 01/09/2024] [Indexed: 02/02/2024]
Abstract
Missing pulse (MP) steady-state free precession (SSFP) is a magnetic resonance imaging (MRI) pulse sequence that is highly tolerant to the magnetic field inhomogeneity. In this study, optimal flip angle and radiofrequency (RF) phase scheduling in three-dimensional (3D) MP-SSFP is introduced to maximize the steady-state magnetization while keeping broadband excitation to cover widely distributed frequencies generated by inhomogeneous magnetic fields. Numerical optimization based on extended phase graph (EPG) simulation was performed to maximize the MP-SSFP steady-state magnetization. To limit the specific absorption rate (SAR) associated with the broadband excitation in 3D MP-SSFP, SAR constraint was introduced in the numerical optimization. Optimized flip angle and RF phase settings were experimentally tested by introducing a linear inhomogeneous magnetic field in a range of 10-20 mT/m and using a phantom with known T1/T2 relaxation and diffusion parameters at 3 T. The experimental results were validated through comparisons with EPG simulation. Image contrasts and molecular diffusion effects were investigated in in vivo human brain imaging with 3D MP-SSFP with the optimal flip angle and RF phase settings. In the phantom measurements, the optimal flip angle and RF phase settings improved the MP-SSFP steady-state magnetization/signal-to-noise ratio by up to 41% under the fixed SAR conditions, which matched well with EPG simulation results. In vivo brain imaging with the optimal RF pulse settings provided T2-like image contrasts. Diffusion effects were relatively minor with the linear inhomogeneous field of 10-20 mT/m for white and gray matter, but cerebrospinal fluid showed conspicuous signal intensity attenuation as the linear inhomogeneous field increased. Numerical optimization achieved significant improvement in the steady-state magnetization in MP-SSFP compared with the RF pulse settings used in previous studies. The proposed flip angle and RF phase optimization is promising to improve 3D MP-SSFP image quality for MRI in inhomogeneous magnetic fields.
Collapse
Affiliation(s)
- Naoharu Kobayashi
- Center for Magnetic Resonance Research, Department of Radiology, University of Minnesota, Minneapolis, Minnesota, USA
| |
Collapse
|
16
|
Liao C, Cao X, Iyer SS, Schauman S, Zhou Z, Yan X, Chen Q, Li Z, Wang N, Gong T, Wu Z, He H, Zhong J, Yang Y, Kerr A, Grill-Spector K, Setsompop K. High-resolution myelin-water fraction and quantitative relaxation mapping using 3D ViSTa-MR fingerprinting. Magn Reson Med 2024; 91:2278-2293. [PMID: 38156945 PMCID: PMC10997479 DOI: 10.1002/mrm.29990] [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/11/2023] [Revised: 12/11/2023] [Accepted: 12/11/2023] [Indexed: 01/03/2024]
Abstract
PURPOSE This study aims to develop a high-resolution whole-brain multi-parametric quantitative MRI approach for simultaneous mapping of myelin-water fraction (MWF), T1, T2, and proton-density (PD), all within a clinically feasible scan time. METHODS We developed 3D visualization of short transverse relaxation time component (ViSTa)-MRF, which combined ViSTa technique with MR fingerprinting (MRF), to achieve high-fidelity whole-brain MWF and T1/T2/PD mapping on a clinical 3T scanner. To achieve fast acquisition and memory-efficient reconstruction, the ViSTa-MRF sequence leverages an optimized 3D tiny-golden-angle-shuffling spiral-projection acquisition and joint spatial-temporal subspace reconstruction with optimized preconditioning algorithm. With the proposed ViSTa-MRF approach, high-fidelity direct MWF mapping was achieved without a need for multicompartment fitting that could introduce bias and/or noise from additional assumptions or priors. RESULTS The in vivo results demonstrate the effectiveness of the proposed acquisition and reconstruction framework to provide fast multi-parametric mapping with high SNR and good quality. The in vivo results of 1 mm- and 0.66 mm-isotropic resolution datasets indicate that the MWF values measured by the proposed method are consistent with standard ViSTa results that are 30× slower with lower SNR. Furthermore, we applied the proposed method to enable 5-min whole-brain 1 mm-iso assessment of MWF and T1/T2/PD mappings for infant brain development and for post-mortem brain samples. CONCLUSIONS In this work, we have developed a 3D ViSTa-MRF technique that enables the acquisition of whole-brain MWF, quantitative T1, T2, and PD maps at 1 and 0.66 mm isotropic resolution in 5 and 15 min, respectively. This advancement allows for quantitative investigations of myelination changes in the brain.
Collapse
Affiliation(s)
- Congyu Liao
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Xiaozhi Cao
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Siddharth Srinivasan Iyer
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Sophie Schauman
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Zihan Zhou
- Department of Radiology, Stanford University, Stanford, CA, USA
- Center for Brain Imaging Science and Technology, College of Biomedical Engineering & Instrument Science, Zhejiang University, Hangzhou, Zhejiang, China
| | - Xiaoqian Yan
- Department of Psychology, Stanford University, Stanford, CA, USA
- Institute of Science and Technology for Brain-Inspired Intelligence, Fudan University, Shanghai, China
| | - Quan Chen
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Zhitao Li
- Department of Radiology, Stanford University, Stanford, CA, USA
| | - Nan Wang
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Ting Gong
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital and Harvard Medical School, Charlestown, MA, USA
| | - Zhe Wu
- Techna Institute, University Health Network, Toronto, ON, Canada
| | - Hongjian He
- Center for Brain Imaging Science and Technology, College of Biomedical Engineering & Instrument Science, Zhejiang University, Hangzhou, Zhejiang, China
- School of Physics, Zhejiang University, Hangzhou, Zhejiang, China
| | - Jianhui Zhong
- Center for Brain Imaging Science and Technology, College of Biomedical Engineering & Instrument Science, Zhejiang University, Hangzhou, Zhejiang, China
- Department of Imaging Sciences, University of Rochester, Rochester, NY, USA
| | - Yang Yang
- Department of Radiology and Biomedical Imaging, University of California, San Francisco, CA, USA
| | - Adam Kerr
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
- Stanford Center for Cognitive and Neurobiological Imaging, Stanford University, Stanford, CA, USA
| | | | - Kawin Setsompop
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| |
Collapse
|
17
|
Martinez JA, Yu VY, Tringale KR, Otazo R, Cohen O. Phase-sensitive deep reconstruction method for rapid multiparametric MR fingerprinting and quantitative susceptibility mapping in the brain. Magn Reson Imaging 2024; 109:147-157. [PMID: 38513790 PMCID: PMC11042874 DOI: 10.1016/j.mri.2024.03.023] [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: 12/01/2023] [Revised: 03/15/2024] [Accepted: 03/16/2024] [Indexed: 03/23/2024]
Abstract
INTRODUCTION This study explores the potential of Magnetic Resonance Fingerprinting (MRF) with a novel Phase-Sensitivity Deep Reconstruction Network (PS-DRONE) for simultaneous quantification of T1, T2, Proton Density, B1+, phase and quantitative susceptibility mapping (QSM). METHODS Data were acquired at 3 T in vitro and in vivo using an optimized EPI-based MRF sequence. Phantom experiments were conducted using a standardized phantom for T1 and T2 maps and a custom-made agar-based gadolinium phantom for B1 and QSM maps. In vivo experiments included five healthy volunteers and one patient diagnosed with brain metastasis. PSDRONE maps were compared to reference maps obtained through standard imaging sequences. RESULTS Total scan time was 2 min for 32 slices and a resolution of [1 mm, 1 mm, 4.5 mm]. The reconstruction of T1, T2, Proton Density, B1+ and phase maps were reconstructed within 1 s. In the phantoms, PS-DRONE analysis presented accurate and strongly correlated T1 and T2 maps (r = 0.99) compared to the reference maps. B1 maps from PS-DRONE showed slightly higher values, though still correlated (r = 0.6) with the reference. QSM values showed a small bias but were strongly correlated (r = 0.99) with reference data. In the in vivo analysis, PS-DRONE-derived T1 and T2 values for gray and white matter matched reference values in healthy volunteers. PS-DRONE B1 and QSM maps showed strong correlations with reference values. CONCLUSION The PS-DRONE network enables concurrent acquisition of T1, T2, PD, B1+, phase and QSM maps, within 2 min of acquisition time and 1 s of reconstruction time.
Collapse
Affiliation(s)
- Jessica A Martinez
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York 10065, NY, USA.
| | - Victoria Y Yu
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York 10065, NY, USA
| | - Kathryn R Tringale
- Department of Radiation Oncology, Memorial Sloan Kettering Cancer Center, New York 10065, NY, USA
| | - Ricardo Otazo
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York 10065, NY, USA
| | - Ouri Cohen
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York 10065, NY, USA
| |
Collapse
|
18
|
Lancione M, Cencini M, Scaffei E, Cipriano E, Buonincontri G, Schulte RF, Pirkl CM, Buchignani B, Pasquariello R, Canapicchi R, Battini R, Biagi L, Tosetti M. Magnetic resonance fingerprinting-based myelin water fraction mapping for the assessment of white matter maturation and integrity in typical development and leukodystrophies. NMR IN BIOMEDICINE 2024; 37:e5114. [PMID: 38390667 DOI: 10.1002/nbm.5114] [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: 08/29/2023] [Revised: 01/10/2024] [Accepted: 01/11/2024] [Indexed: 02/24/2024]
Abstract
A quantitative biomarker for myelination, such as myelin water fraction (MWF), would boost the understanding of normative and pathological neurodevelopment, improving patients' diagnosis and follow-up. We quantified the fraction of a rapidly relaxing pool identified as MW using multicomponent three-dimensional (3D) magnetic resonance fingerprinting (MRF) to evaluate white matter (WM) maturation in typically developing (TD) children and alterations in leukodystrophies (LDs). We acquired DTI and 3D MRF-based R1, R2 and MWF data of 15 TD children and 17 LD patients (9 months-12.5 years old) at 1.5 T. We computed normative maturation curves in corpus callosum and corona radiata and performed WM tract profile analysis, comparing MWF with R1, R2 and fractional anisotropy (FA). Normative maturation curves demonstrated a steep increase for all tissue parameters in the first 3 years of age, followed by slower growth for MWF while R1, R2R2 and FA reached a plateau. Unlike FA, MWF values were similar for regions of interest (ROIs) with different degrees of axonal packing, suggesting independence from fiber bundle macro-organization and higher myelin specificity. Tract profile analysis indicated a specific spatial pattern of myelination in the major fiber bundles, consistent across subjects. LD were better distinguished from TD by MWF rather than FA, showing reduced MWF with respect to age-matched controls in both ROI-based and tract analysis. In conclusion, MRF-based MWF provides myelin-specific WM maturation curves and is sensitive to alteration due to LDs, suggesting its potential as a biomarker for WM disorders. As MRF allows fast simultaneous acquisition of relaxometry and MWF, it can represent a valuable diagnostic tool to study and follow up developmental WM disorders in children.
Collapse
Affiliation(s)
| | - Matteo Cencini
- Pisa Division, National Institute for Nuclear Physics (INFN), Pisa, Italy
| | | | - Emilio Cipriano
- IRCCS Stella Maris Foundation, Pisa, Italy
- Department of Physics, University of Pisa, Pisa, Italy
| | | | | | | | | | | | | | - Roberta Battini
- IRCCS Stella Maris Foundation, Pisa, Italy
- Department of Clinical and Experimental Medicine, Università di Pisa, Pisa, Italy
| | | | | |
Collapse
|
19
|
Jun Y, Arefeen Y, Cho J, Fujita S, Wang X, Ellen Grant P, Gagoski B, Jaimes C, Gee MS, Bilgic B. Zero-DeepSub: Zero-shot deep subspace reconstruction for rapid multiparametric quantitative MRI using 3D-QALAS. Magn Reson Med 2024; 91:2459-2482. [PMID: 38282270 PMCID: PMC11005062 DOI: 10.1002/mrm.30018] [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/03/2023] [Revised: 12/15/2023] [Accepted: 01/06/2024] [Indexed: 01/30/2024]
Abstract
PURPOSE To develop and evaluate methods for (1) reconstructing 3D-quantification using an interleaved Look-Locker acquisition sequence with T2 preparation pulse (3D-QALAS) time-series images using a low-rank subspace method, which enables accurate and rapid T1 and T2 mapping, and (2) improving the fidelity of subspace QALAS by combining scan-specific deep-learning-based reconstruction and subspace modeling. THEORY AND METHODS A low-rank subspace method for 3D-QALAS (i.e., subspace QALAS) and zero-shot deep-learning subspace method (i.e., Zero-DeepSub) were proposed for rapid and high fidelity T1 and T2 mapping and time-resolved imaging using 3D-QALAS. Using an ISMRM/NIST system phantom, the accuracy and reproducibility of the T1 and T2 maps estimated using the proposed methods were evaluated by comparing them with reference techniques. The reconstruction performance of the proposed subspace QALAS using Zero-DeepSub was evaluated in vivo and compared with conventional QALAS at high reduction factors of up to nine-fold. RESULTS Phantom experiments showed that subspace QALAS had good linearity with respect to the reference methods while reducing biases and improving precision compared to conventional QALAS, especially for T2 maps. Moreover, in vivo results demonstrated that subspace QALAS had better g-factor maps and could reduce voxel blurring, noise, and artifacts compared to conventional QALAS and showed robust performance at up to nine-fold acceleration with Zero-DeepSub, which enabled whole-brain T1, T2, and PD mapping at 1 mm isotropic resolution within 2 min of scan time. CONCLUSION The proposed subspace QALAS along with Zero-DeepSub enabled high fidelity and rapid whole-brain multiparametric quantification and time-resolved imaging.
Collapse
Affiliation(s)
- Yohan Jun
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Boston, MA, United States
- Department of Radiology, Harvard Medical School, Boston, MA, United States
| | - Yamin Arefeen
- Chandra Family Department of Electrical and Computer Engineering, The University of Texas, Austin, TX, United States
- Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, Houston, TX, United States
| | - Jaejin Cho
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Boston, MA, United States
- Department of Radiology, Harvard Medical School, Boston, MA, United States
| | - Shohei Fujita
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Boston, MA, United States
- Department of Radiology, Harvard Medical School, Boston, MA, United States
| | - Xiaoqing Wang
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Boston, MA, United States
- Department of Radiology, Harvard Medical School, Boston, MA, United States
| | - P. Ellen Grant
- Department of Radiology, Harvard Medical School, Boston, MA, United States
- Fetal-Neonatal Neuroimaging & Developmental Science Center, Boston Children’s Hospital, Boston, MA, United States
| | - Borjan Gagoski
- Department of Radiology, Harvard Medical School, Boston, MA, United States
- Fetal-Neonatal Neuroimaging & Developmental Science Center, Boston Children’s Hospital, Boston, MA, United States
| | - Camilo Jaimes
- Department of Radiology, Harvard Medical School, Boston, MA, United States
- Department of Radiology, Massachusetts General Hospital, Boston, MA, United States
| | - Michael S. Gee
- Department of Radiology, Harvard Medical School, Boston, MA, United States
- Department of Radiology, Massachusetts General Hospital, Boston, MA, United States
| | - Berkin Bilgic
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Boston, MA, United States
- Department of Radiology, Harvard Medical School, Boston, MA, United States
- Harvard/MIT Health Sciences and Technology, Cambridge, MA, United States
| |
Collapse
|
20
|
Brackenier Y, Wang N, Liao C, Cao X, Schauman S, Yurt M, Cordero-Grande L, Malik SJ, Kerr A, Hajnal JV, Setsompop K. Rapid and accurate navigators for motion and B 0 tracking using QUEEN: Quantitatively enhanced parameter estimation from navigators. Magn Reson Med 2024; 91:2028-2043. [PMID: 38173304 DOI: 10.1002/mrm.29976] [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: 06/17/2023] [Revised: 11/28/2023] [Accepted: 11/29/2023] [Indexed: 01/05/2024]
Abstract
PURPOSE To develop a framework that jointly estimates rigid motion and polarizing magnetic field (B0 ) perturbations (δ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ ) for brain MRI using a single navigator of a few milliseconds in duration, and to additionally allow for navigator acquisition at arbitrary timings within any type of sequence to obtain high-temporal resolution estimates. THEORY AND METHODS Methods exist that match navigator data to a low-resolution single-contrast image (scout) to estimate either motion orδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ . In this work, called QUEEN (QUantitatively Enhanced parameter Estimation from Navigators), we propose combined motion andδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ estimation from a fast, tailored trajectory with arbitrary-contrast navigator data. To this end, the concept of a quantitative scout (Q-Scout) acquisition is proposed from which contrast-matched scout data is predicted for each navigator. Finally, navigator trajectories, contrast-matched scout, andδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ are integrated into a motion-informed parallel-imaging framework. RESULTS Simulations and in vivo experiments show the need to modelδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ to obtain accurate motion parameters estimated in the presence of strongδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ . Simulations confirm that tailored navigator trajectories are needed to robustly estimate both motion andδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ . Furthermore, experiments show that a contrast-matched scout is needed for parameter estimation from multicontrast navigator data. A retrospective, in vivo reconstruction experiment shows improved image quality when using the proposed Q-Scout and QUEEN estimation. CONCLUSIONS We developed a framework to jointly estimate rigid motion parameters andδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ from navigators. Combing a contrast-matched scout with the proposed trajectory allows for navigator deployment in almost any sequence and/or timing, which allows for higher temporal-resolution motion andδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ estimates.
Collapse
Affiliation(s)
| | - Nan Wang
- Department of Radiology, Stanford University, Stanford, California, USA
| | - Congyu Liao
- Department of Radiology, Stanford University, Stanford, California, USA
| | - Xiaozhi Cao
- Department of Radiology, Stanford University, Stanford, California, USA
| | - Sophie Schauman
- Department of Radiology, Stanford University, Stanford, California, USA
| | - Mahmut Yurt
- Department of Radiology, Stanford University, Stanford, California, USA
- Department of Electrical Engineering, Stanford University, Stanford, California, USA
| | - Lucilio Cordero-Grande
- Biomedical Image Technologies, ETSI Telecomunicación, Universidad Politécnica de Madrid and CIBER-BNN, Madrid, Spain
| | - Shaihan J Malik
- Biomedical Engineering Department, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- Center for the Developing Brain, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Adam Kerr
- Department of Electrical Engineering, Stanford University, Stanford, California, USA
- Cognitive and Neurobiological Imaging, Stanford University, Stanford, California, USA
| | - Joseph V Hajnal
- Biomedical Engineering Department, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- Center for the Developing Brain, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Kawin Setsompop
- Department of Radiology, Stanford University, Stanford, California, USA
- Department of Electrical Engineering, Stanford University, Stanford, California, USA
| |
Collapse
|
21
|
Ayaz A, Al Khalil Y, Amirrajab S, Lorenz C, Weese J, Pluim J, Breeuwer M. Brain MR image simulation for deep learning based medical image analysis networks. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2024; 248:108115. [PMID: 38503072 DOI: 10.1016/j.cmpb.2024.108115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/15/2023] [Revised: 02/02/2024] [Accepted: 03/02/2024] [Indexed: 03/21/2024]
Abstract
BACKGROUND AND OBJECTIVE As large sets of annotated MRI data are needed for training and validating deep learning based medical image analysis algorithms, the lack of sufficient annotated data is a critical problem. A possible solution is the generation of artificial data by means of physics-based simulations. Existing brain simulation data is limited in terms of anatomical models, tissue classes, fixed tissue characteristics, MR sequences and overall realism. METHODS We propose a realistic simulation framework by incorporating patient-specific phantoms and Bloch equations-based analytical solutions for fast and accurate MRI simulations. A large number of labels are derived from open-source high-resolution T1w MRI data using a fully automated brain classification tool. The brain labels are taken as ground truth (GT) on which MR images are simulated using our framework. Moreover, we demonstrate that the T1w MR images generated from our framework along with GT annotations can be utilized directly to train a 3D brain segmentation network. To evaluate our model further on larger set of real multi-source MRI data without GT, we compared our model to existing brain segmentation tools, FSL-FAST and SynthSeg. RESULTS Our framework generates 3D brain MRI for variable anatomy, sequence, contrast, SNR and resolution. The brain segmentation network for WM/GM/CSF trained only on T1w simulated data shows promising results on real MRI data from MRBrainS18 challenge dataset with a Dice scores of 0.818/0.832/0.828. On OASIS data, our model exhibits a close performance to FSL, both qualitatively and quantitatively with a Dice scores of 0.901/0.939/0.937. CONCLUSIONS Our proposed simulation framework is the initial step towards achieving truly physics-based MRI image generation, providing flexibility to generate large sets of variable MRI data for desired anatomy, sequence, contrast, SNR, and resolution. Furthermore, the generated images can effectively train 3D brain segmentation networks, mitigating the reliance on real 3D annotated data.
Collapse
Affiliation(s)
- Aymen Ayaz
- Biomedical Engineering Department, Eindhoven University of Technology, Eindhoven, the Netherlands.
| | - Yasmina Al Khalil
- Biomedical Engineering Department, Eindhoven University of Technology, Eindhoven, the Netherlands.
| | - Sina Amirrajab
- Biomedical Engineering Department, Eindhoven University of Technology, Eindhoven, the Netherlands.
| | | | - Jürgen Weese
- Philips Research Laboratories, Hamburg, Germany.
| | - Josien Pluim
- Biomedical Engineering Department, Eindhoven University of Technology, Eindhoven, the Netherlands.
| | - Marcel Breeuwer
- Biomedical Engineering Department, Eindhoven University of Technology, Eindhoven, the Netherlands; MR R&D - Clinical Science, Philips Healthcare, Best, the Netherlands.
| |
Collapse
|
22
|
Schneider A, Munoz C, Hua A, Ellis S, Jeljeli S, Kunze KP, Neji R, Reader AJ, Reyes E, Ismail TF, Botnar RM, Prieto C. Non-rigid motion-compensated 3D whole-heart T 2 mapping in a hybrid 3T PET-MR system. Magn Reson Med 2024; 91:1951-1964. [PMID: 38181169 DOI: 10.1002/mrm.29973] [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: 06/07/2023] [Revised: 11/21/2023] [Accepted: 11/26/2023] [Indexed: 01/07/2024]
Abstract
PURPOSE Simultaneous PET-MRI improves inflammatory cardiac disease diagnosis. However, challenges persist in respiratory motion and mis-registration between free-breathing 3D PET and 2D breath-held MR images. We propose a free-breathing non-rigid motion-compensated 3D T2 -mapping sequence enabling whole-heart myocardial tissue characterization in a hybrid 3T PET-MR system and provides non-rigid respiratory motion fields to correct also simultaneously acquired PET data. METHODS Free-breathing 3D whole-heart T2 -mapping was implemented on a hybrid 3T PET-MRI system. Three datasets were acquired with different T2 -preparation modules (0, 28, 55 ms) using 3-fold undersampled variable-density Cartesian trajectory. Respiratory motion was estimated via virtual 3D image navigators, enabling multi-contrast non-rigid motion-corrected MR reconstruction. T2 -maps were computed using dictionary-matching. Approach was tested in phantom, 8 healthy subjects, 14 MR only and 2 PET-MR patients with suspected cardiac disease and compared with spin echo reference (phantom) and clinical 2D T2 -mapping (in-vivo). RESULTS Phantom results show a high correlation (R2 = 0.996) between proposed approach and gold standard 2D T2 mapping. In-vivo 3D T2 -mapping average values in healthy subjects (39.0 ± 1.4 ms) and patients (healthy tissue) (39.1 ± 1.4 ms) agree with conventional 2D T2 -mapping (healthy = 38.6 ± 1.2 ms, patients = 40.3 ± 1.7 ms). Bland-Altman analysis reveals bias of 1.8 ms and 95% limits of agreement (LOA) of -2.4-6 ms for healthy subjects, and bias of 1.3 ms and 95% LOA of -1.9 to 4.6 ms for patients. CONCLUSION Validated efficient 3D whole-heart T2 -mapping at hybrid 3T PET-MRI provides myocardial inflammation characterization and non-rigid respiratory motion fields for simultaneous PET data correction. Comparable T2 values were achieved with both 3D and 2D methods. Improved image quality was observed in the PET images after MR-based motion correction.
Collapse
Affiliation(s)
- Alina Schneider
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Camila Munoz
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Alina Hua
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Sam Ellis
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Sami Jeljeli
- PET Centre, St Thomas' Hospital, King's College London & Guys and St Thomas' NHS Foundation Trust, London, UK
| | - Karl P Kunze
- MR Research Collaborations, Siemens Healthcare Limited, Camberley, UK
| | - Radhouene Neji
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Andrew J Reader
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Eliana Reyes
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Tevfik F Ismail
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - René M Botnar
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- Escuela de Ingeniería, Pontificia Universidad Católica de Chile, Santiago, Chile
- Millenium Institute for Intelligent Healthcare Engineering iHEALTH, Santiago, Chile
| | - Claudia Prieto
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
- Escuela de Ingeniería, Pontificia Universidad Católica de Chile, Santiago, Chile
- Millenium Institute for Intelligent Healthcare Engineering iHEALTH, Santiago, Chile
| |
Collapse
|
23
|
Huaroc Moquillaza E, Weiss K, Stelter J, Steinhelfer L, Lee YJ, Amthor T, Koken P, Makowski MR, Braren R, Doneva M, Karampinos DC. Accelerated liver water T 1 mapping using single-shot continuous inversion-recovery spiral imaging. NMR IN BIOMEDICINE 2024; 37:e5097. [PMID: 38269568 DOI: 10.1002/nbm.5097] [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: 08/09/2023] [Revised: 11/21/2023] [Accepted: 12/06/2023] [Indexed: 01/26/2024]
Abstract
PURPOSE Liver T1 mapping techniques typically require long breath holds or long scan time in free-breathing, need correction for B 1 + inhomogeneities and process composite (water and fat) signals. The purpose of this work is to accelerate the multi-slice acquisition of liver water selective T1 (wT1) mapping in a single breath hold, improving the k-space sampling efficiency. METHODS The proposed continuous inversion-recovery (IR) Look-Locker methodology combines a single-shot gradient echo spiral readout, Dixon processing and a dictionary-based analysis for liver wT1 mapping at 3 T. The sequence parameters were adapted to obtain short scan times. The influence of fat, B 1 + inhomogeneities and TE on the estimation of T1 was first assessed using simulations. The proposed method was then validated in a phantom and in 10 volunteers, comparing it with MRS and the modified Look-Locker inversion-recovery (MOLLI) method. Finally, the clinical feasibility was investigated by comparing wT1 maps with clinical scans in nine patients. RESULTS The phantom results are in good agreement with MRS. The proposed method encodes the IR-curve for the liver wT1 estimation, is minimally sensitive to B 1 + inhomogeneities and acquires one slice in 1.2 s. The volunteer results confirmed the multi-slice capability of the proposed method, acquiring nine slices in a breath hold of 11 s. The present work shows robustness to B 1 + inhomogeneities (wT 1 , No B 1 + = 1.07 wT 1 , B 1 + - 45.63 , R 2 = 0.99 ) , good repeatability (wT 1 , 2 ° = 1 . 0 wT 1 , 1 ° - 2.14 , R 2 = 0.96 ) and is in better agreement with MRS (wT 1 = 0.92 wT 1 MRS + 103.28 , R 2 = 0.38 ) than is MOLLI (wT 1 MOLLI = 0.76 wT 1 MRS + 254.43 , R 2 = 0.44 ) . The wT1 maps in patients captured diverse lesions, thus showing their clinical feasibility. CONCLUSION A single-shot spiral acquisition can be combined with a continuous IR Look-Locker method to perform rapid repeatable multi-slice liver water T1 mapping at a rate of 1.2 s per slice without a B 1 + map. The proposed method is suitable for nine-slice liver clinical applications acquired in a single breath hold of 11 s.
Collapse
Affiliation(s)
- Elizabeth Huaroc Moquillaza
- Department of Diagnostic and Interventional Radiology, School of Medicine, Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
| | | | - Jonathan Stelter
- Department of Diagnostic and Interventional Radiology, School of Medicine, Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
| | - Lisa Steinhelfer
- Department of Diagnostic and Interventional Radiology, School of Medicine, Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
| | | | | | | | - Marcus R Makowski
- Department of Diagnostic and Interventional Radiology, School of Medicine, Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
| | - Rickmer Braren
- Department of Diagnostic and Interventional Radiology, School of Medicine, Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
| | | | - Dimitrios C Karampinos
- Department of Diagnostic and Interventional Radiology, School of Medicine, Klinikum rechts der Isar, Technical University of Munich, Munich, Germany
| |
Collapse
|
24
|
Menon RG, Monga A, Kijowski R, Regatte RR. Characterization of Age-Related and Sex-Related Differences of Relaxation Parameters in the Intervertebral Disc Using MR-Fingerprinting. J Magn Reson Imaging 2024; 59:1312-1324. [PMID: 37610269 PMCID: PMC10935608 DOI: 10.1002/jmri.28925] [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/14/2023] [Revised: 07/13/2023] [Accepted: 07/14/2023] [Indexed: 08/24/2023] Open
Abstract
BACKGROUND Multiparameter characterization using MR fingerprinting (MRF) can quantify multiple relaxation parameters of intervertebral disc (IVD) simultaneously. These parameters may vary by age and sex. PURPOSE To investigate age- and sex-related differences in the relaxation parameters of the IVD of the lumbar spine using a multiparameter MRF technique. STUDY TYPE Prospective. SUBJECTS 17 healthy subjects (8 male; mean age = 34 ± 10 years, range 20-60 years). FIELD STRENGTH/SEQUENCE 3D-MRF sequence for simultaneous acquisition of proton density, T1 , T2 , and T1ρ maps at 3.0T. ASSESSMENT Global mean T1 , T2 , and T1ρ of all lumbar IVDs and mean T1 , T2 , and T1ρ of each individual IVD (L1-L5) were measured. Gray level co-occurrence matrix was used to quantify textural features (median, contrast, correlation, energy, and homogeneity) from T1 , T2 , and T1ρ maps. STATISTICAL TESTS Spearman rank correlations (R) evaluated the association between age and T1 , T2 , and T1ρ of IVD. Mann-Whitney U-tests evaluated differences between males and females in T1 , T2 , and T1ρ of IVD. Statistical significance was defined as P-value <0.05. RESULTS There was a significant negative correlation between age and global mean values of all IVDs for T1 (R = -0.637), T2 (R = -0.509), and T1ρ (R = -0.726). For individual IVDs, there was a significant negative correlation between age and mean T1 at all IVD segments (R range = -0.530 to -0.708), between age and mean T2 at L2-L3, L3-L4, and L4-L5 (R range = -0.493 to 0.640), and between age and mean T1ρ at all segments except L1-L2 (R range = -0.632 to -0.763). There were no significant differences between sexes in global mean T1 , T2, and T1ρ (P-value = 0.23-0.76) The texture features with the highest significant correlations with age for all IVDs were global T1ρ mean (R = -0.726), T1 energy (R = -0.681), and T1 contrast (R = 0.709). CONCLUSION This study showed that the 3D-MRF technique has potential to characterize age-related differences in T1 , T2, or T1ρ of IVD in healthy subjects. LEVEL OF EVIDENCE 2 TECHNICAL EFFICACY: Stage 1.
Collapse
Affiliation(s)
- Rajiv G. Menon
- Center for Biomedical Imaging, Department of Radiology, New York University Grossman School of Medicine
| | - Anmol Monga
- Center for Biomedical Imaging, Department of Radiology, New York University Grossman School of Medicine
| | - Richard Kijowski
- Center for Biomedical Imaging, Department of Radiology, New York University Grossman School of Medicine
| | - Ravinder R. Regatte
- Center for Biomedical Imaging, Department of Radiology, New York University Grossman School of Medicine
| |
Collapse
|
25
|
Marriott A, Rioux J, Brewer K. Nonuniform sliding-window reconstruction for accelerated dual contrast agent quantification with MR fingerprinting. MAGMA (NEW YORK, N.Y.) 2024; 37:273-282. [PMID: 38217784 PMCID: PMC10994993 DOI: 10.1007/s10334-023-01140-9] [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: 10/17/2023] [Revised: 12/07/2023] [Accepted: 12/11/2023] [Indexed: 01/15/2024]
Abstract
OBJECTIVE MR fingerprinting (MRF) can enable preclinical studies of cell tracking by quantifying multiple contrast agents simultaneously, but faster scan times are required for in vivo applications. Sliding window (SW)-MRF is one option for accelerating MRF, but standard implementations are not sufficient to preserve the accuracy of T2*, which is critical for tracking iron-labelled cells in vivo. PURPOSE To develop a SW approach to MRF which preserves the T2* accuracy required for accelerated concentration mapping of iron-labelled cells on single-channel preclinical systems. METHODS A nonuniform SW was applied to the MRF sequence and dictionary. Segments of the sequence most sensitive to T2* were subject to a shorter window length, preserving the T2* sensitivity. Phantoms containing iron-labelled CD8+ T cells and gadolinium were used to compare 24× undersampled uniform and nonuniform SW-MRF parameter maps. Dual concentration maps were generated for both uniform and nonuniform MRF and compared. RESULTS Lin's concordance correlation coefficient, compared to gold standard parameter values, was much greater for nonuniform SW-MRF than for uniform SW-MRF. A Wilcoxon signed-rank test showed no significant difference between nonuniform SW-MRF and gold standards. Nonuniform SW-MRF outperformed the uniform SW-MRF concentration maps for all parameters, providing a balance between T2* sensitivity of short window lengths, and SNR of longer window lengths. CONCLUSIONS Nonuniform SW-MRF improves the accuracy of matching compared to uniform SW-MRF, allowing higher accelerated concentration mapping for preclinical systems.
Collapse
Affiliation(s)
- Anna Marriott
- Biomedical MRI Research Laboratory (BMRL), IWK Health Centre, 5850/5980 University Avenue, Halifax, NS, B3K 6R8, Canada
- Department of Physics and Atmospheric Science, Dalhousie University, Halifax, NS, Canada
| | - James Rioux
- Department of Physics and Atmospheric Science, Dalhousie University, Halifax, NS, Canada
- Biomedical Translational Imaging Centre (BIOTIC), NS Health, Halifax, NS, Canada
- Department of Diagnostic Radiology, Dalhousie University, Halifax, NS, Canada
| | - Kimberly Brewer
- Biomedical MRI Research Laboratory (BMRL), IWK Health Centre, 5850/5980 University Avenue, Halifax, NS, B3K 6R8, Canada.
- Department of Physics and Atmospheric Science, Dalhousie University, Halifax, NS, Canada.
- Department of Diagnostic Radiology, Dalhousie University, Halifax, NS, Canada.
- School of Biomedical Engineering, Dalhousie University, Halifax, NS, Canada.
- Department of Microbiology and Immunology, Dalhousie University, Halifax, NS, Canada.
| |
Collapse
|
26
|
Assländer J, Gultekin C, Mao A, Zhang X, Duchemin Q, Liu K, Charlson RW, Shepherd TM, Fernandez-Granda C, Flassbeck S. Rapid quantitative magnetization transfer imaging: Utilizing the hybrid state and the generalized Bloch model. Magn Reson Med 2024; 91:1478-1497. [PMID: 38073093 DOI: 10.1002/mrm.29951] [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: 06/15/2023] [Revised: 10/30/2023] [Accepted: 11/14/2023] [Indexed: 02/03/2024]
Abstract
PURPOSE To explore efficient encoding schemes for quantitative magnetization transfer (qMT) imaging with few constraints on model parameters. THEORY AND METHODS We combine two recently proposed models in a Bloch-McConnell equation: the dynamics of the free spin pool are confined to the hybrid state, and the dynamics of the semi-solid spin pool are described by the generalized Bloch model. We numerically optimize the flip angles and durations of a train of radio frequency pulses to enhance the encoding of three qMT parameters while accounting for all eight parameters of the two-pool model. We sparsely sample each time frame along this spin dynamics with a three-dimensional radial koosh-ball trajectory, reconstruct the data with subspace modeling, and fit the qMT model with a neural network for computational efficiency. RESULTS We extracted qMT parameter maps of the whole brain with an effective resolution of 1.24 mm from a 12.6-min scan. In lesions of multiple sclerosis subjects, we observe a decreased size of the semi-solid spin pool and longer relaxation times, consistent with previous reports. CONCLUSION The encoding power of the hybrid state, combined with regularized image reconstruction, and the accuracy of the generalized Bloch model provide an excellent basis for efficient quantitative magnetization transfer imaging with few constraints on model parameters.
Collapse
Affiliation(s)
- Jakob Assländer
- Center for Biomedical Imaging, Department of Radiology, NYU School of Medicine, New York, New York, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), Department of Radiology, NYU School of Medicine, New York, New York, USA
| | - Cem Gultekin
- Courant Institute of Mathematical Sciences, New York University, New York, New York, USA
| | - Andrew Mao
- Center for Biomedical Imaging, Department of Radiology, NYU School of Medicine, New York, New York, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), Department of Radiology, NYU School of Medicine, New York, New York, USA
- Vilcek Institute of Graduate Biomedical Sciences, NYU School of Medicine, New York, New York, USA
| | - Xiaoxia Zhang
- Center for Biomedical Imaging, Department of Radiology, NYU School of Medicine, New York, New York, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), Department of Radiology, NYU School of Medicine, New York, New York, USA
| | - Quentin Duchemin
- Laboratoire d'analyse et de mathématiques appliquées, Université Gustave Eiffel, Champs-sur-Marne, France
| | - Kangning Liu
- Center for Data Science, New York University, New York, New York, USA
| | - Robert W Charlson
- Department of Neurology, NYU School of Medicine, New York, New York, USA
| | - Timothy M Shepherd
- Center for Biomedical Imaging, Department of Radiology, NYU School of Medicine, New York, New York, USA
| | - Carlos Fernandez-Granda
- Courant Institute of Mathematical Sciences, New York University, New York, New York, USA
- Center for Data Science, New York University, New York, New York, USA
| | - Sebastian Flassbeck
- Center for Biomedical Imaging, Department of Radiology, NYU School of Medicine, New York, New York, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), Department of Radiology, NYU School of Medicine, New York, New York, USA
| |
Collapse
|
27
|
Cao X, Liao C, Zhou Z, Zhong Z, Li Z, Dai E, Iyer SS, Hannum AJ, Yurt M, Schauman S, Chen Q, Wang N, Wei J, Yan Y, He H, Skare S, Zhong J, Kerr A, Setsompop K. DTI-MR fingerprinting for rapid high-resolution whole-brain T 1 , T 2 , proton density, ADC, and fractional anisotropy mapping. Magn Reson Med 2024; 91:987-1001. [PMID: 37936313 PMCID: PMC11068310 DOI: 10.1002/mrm.29916] [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/14/2023] [Revised: 10/17/2023] [Accepted: 10/18/2023] [Indexed: 11/09/2023]
Abstract
PURPOSE This study aims to develop a high-efficiency and high-resolution 3D imaging approach for simultaneous mapping of multiple key tissue parameters for routine brain imaging, including T1 , T2 , proton density (PD), ADC, and fractional anisotropy (FA). The proposed method is intended for pushing routine clinical brain imaging from weighted imaging to quantitative imaging and can also be particularly useful for diffusion-relaxometry studies, which typically suffer from lengthy acquisition time. METHODS To address challenges associated with diffusion weighting, such as shot-to-shot phase variation and low SNR, we integrated several innovative data acquisition and reconstruction techniques. Specifically, we used M1-compensated diffusion gradients, cardiac gating, and navigators to mitigate phase variations caused by cardiac motion. We also introduced a data-driven pre-pulse gradient to cancel out eddy currents induced by diffusion gradients. Additionally, to enhance image quality within a limited acquisition time, we proposed a data-sharing joint reconstruction approach coupled with a corresponding sequence design. RESULTS The phantom and in vivo studies indicated that the T1 and T2 values measured by the proposed method are consistent with a conventional MR fingerprinting sequence and the diffusion results (including diffusivity, ADC, and FA) are consistent with the spin-echo EPI DWI sequence. CONCLUSION The proposed method can achieve whole-brain T1 , T2 , diffusivity, ADC, and FA maps at 1-mm isotropic resolution within 10 min, providing a powerful tool for investigating the microstructural properties of brain tissue, with potential applications in clinical and research settings.
Collapse
Affiliation(s)
- Xiaozhi Cao
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Congyu Liao
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Zihan Zhou
- Department of Radiology, Stanford University, Stanford, CA, USA
| | - Zheng Zhong
- Department of Radiology, Stanford University, Stanford, CA, USA
| | - Zhitao Li
- Department of Radiology, Stanford University, Stanford, CA, USA
| | - Erpeng Dai
- Department of Radiology, Stanford University, Stanford, CA, USA
| | - Siddharth Srinivasan Iyer
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering and Computer Science, MIT, Cambridge, Massachusetts, USA
| | - Ariel J Hannum
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Bioengineering, Stanford University, Stanford, CA, USA
| | - Mahmut Yurt
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Sophie Schauman
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Quan Chen
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Nan Wang
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Jintao Wei
- Center for Brain Imaging Science and Technology, Zhejiang University, Hangzhou, Zhejiang, China
| | - Yifan Yan
- School of Public Health and the Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
| | - Hongjian He
- Center for Brain Imaging Science and Technology, Zhejiang University, Hangzhou, Zhejiang, China
- School of Physics, Zhejiang University, Hangzhou, Zhejiang, China
| | - Stefan Skare
- Department of Clinical Neuroscience, Karolinska Institute, Stockholm, Sweden
| | - Jianhui Zhong
- Department of Imaging Sciences, University of Rochester, NY, USA
| | - Adam Kerr
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Kawin Setsompop
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| |
Collapse
|
28
|
Slioussarenko C, Baudin PY, Reyngoudt H, Marty B. Bi-component dictionary matching for MR fingerprinting for efficient quantification of fat fraction and water T 1 in skeletal muscle. Magn Reson Med 2024; 91:1179-1189. [PMID: 37867467 DOI: 10.1002/mrm.29901] [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/28/2023] [Revised: 09/15/2023] [Accepted: 10/06/2023] [Indexed: 10/24/2023]
Abstract
PURPOSE To propose an efficient bi-component MR fingerprinting (MRF) fitting method using a Variable Projection (VARPRO) strategy, applied to the quantification of fat fraction (FF) and water T1 (T 1 H 2 0 $$ \mathrm{T}{1}_{{\mathrm{H}}_20} $$ ) in skeletal muscle tissues. METHODS The MRF signals were analyzed in a two-step process by comparing them to the elements of separate water and fat dictionaries (bi-component dictionary matching). First, each pair of water and fat dictionary elements was fitted to the acquired signal to determine an optimal FF that was used to merge the fingerprints in a combined water/fat dictionary. Second, standard dictionary matching was applied to the combined dictionary for determining the remaining parameters. A clustering method was implemented to further accelerate the fitting. Accuracy, precision, and matching time of this approach were evaluated on both numerical and in vivo datasets, and compared to the reference dictionary-matching approach that includes FF as a dictionary parameter. RESULTS In numerical phantoms, all MRF parameters showed high correlation with ground truth for the reference and the bi-component method (R2 > 0.98). In vivo, the estimated parameters from the proposed method were highly correlated with those from the reference approach (R2 > 0.997). The bi-component method achieved an acceleration factor of up to 360 compared to the reference dictionary matching. CONCLUSION The proposed bi-component fitting approach enables a significant acceleration of the reconstruction of MRF parameter maps for fat-water imaging, while maintaining comparable precision and accuracy to the reference on FF andT 1 H 2 0 $$ \mathrm{T}{1}_{{\mathrm{H}}_20} $$ estimation.
Collapse
Affiliation(s)
| | - Pierre-Yves Baudin
- Institute of Myology, Neuromuscular Investigation Center, NMR Laboratory, Paris, France
| | - Harmen Reyngoudt
- Institute of Myology, Neuromuscular Investigation Center, NMR Laboratory, Paris, France
| | - Benjamin Marty
- Institute of Myology, Neuromuscular Investigation Center, NMR Laboratory, Paris, France
| |
Collapse
|
29
|
Monga A, Singh D, de Moura HL, Zhang X, Zibetti MVW, Regatte RR. Emerging Trends in Magnetic Resonance Fingerprinting for Quantitative Biomedical Imaging Applications: A Review. Bioengineering (Basel) 2024; 11:236. [PMID: 38534511 DOI: 10.3390/bioengineering11030236] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2024] [Revised: 02/21/2024] [Accepted: 02/22/2024] [Indexed: 03/28/2024] Open
Abstract
Magnetic resonance imaging (MRI) stands as a vital medical imaging technique, renowned for its ability to offer high-resolution images of the human body with remarkable soft-tissue contrast. This enables healthcare professionals to gain valuable insights into various aspects of the human body, including morphology, structural integrity, and physiological processes. Quantitative imaging provides compositional measurements of the human body, but, currently, either it takes a long scan time or is limited to low spatial resolutions. Undersampled k-space data acquisitions have significantly helped to reduce MRI scan time, while compressed sensing (CS) and deep learning (DL) reconstructions have mitigated the associated undersampling artifacts. Alternatively, magnetic resonance fingerprinting (MRF) provides an efficient and versatile framework to acquire and quantify multiple tissue properties simultaneously from a single fast MRI scan. The MRF framework involves four key aspects: (1) pulse sequence design; (2) rapid (undersampled) data acquisition; (3) encoding of tissue properties in MR signal evolutions or fingerprints; and (4) simultaneous recovery of multiple quantitative spatial maps. This paper provides an extensive literature review of the MRF framework, addressing the trends associated with these four key aspects. There are specific challenges in MRF for all ranges of magnetic field strengths and all body parts, which can present opportunities for further investigation. We aim to review the best practices in each key aspect of MRF, as well as for different applications, such as cardiac, brain, and musculoskeletal imaging, among others. A comprehensive review of these applications will enable us to assess future trends and their implications for the translation of MRF into these biomedical imaging applications.
Collapse
Affiliation(s)
- Anmol Monga
- Center of Biomedical Imaging, Department of Radiology, New York University Grossman School of Medicine, New York, NY 10016, USA
| | - Dilbag Singh
- Center of Biomedical Imaging, Department of Radiology, New York University Grossman School of Medicine, New York, NY 10016, USA
| | - Hector L de Moura
- Center of Biomedical Imaging, Department of Radiology, New York University Grossman School of Medicine, New York, NY 10016, USA
| | - Xiaoxia Zhang
- Center of Biomedical Imaging, Department of Radiology, New York University Grossman School of Medicine, New York, NY 10016, USA
| | - Marcelo V W Zibetti
- Center of Biomedical Imaging, Department of Radiology, New York University Grossman School of Medicine, New York, NY 10016, USA
| | - Ravinder R Regatte
- Center of Biomedical Imaging, Department of Radiology, New York University Grossman School of Medicine, New York, NY 10016, USA
| |
Collapse
|
30
|
Kumar S, Saber H, Charron O, Freeman L, Tamir JI. Correcting synthetic MRI contrast-weighted images using deep learning. Magn Reson Imaging 2024; 106:43-54. [PMID: 38092082 DOI: 10.1016/j.mri.2023.11.015] [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: 08/30/2023] [Revised: 11/30/2023] [Accepted: 11/30/2023] [Indexed: 12/17/2023]
Abstract
Synthetic magnetic resonance imaging (MRI) offers a scanning paradigm where a fast multi-contrast sequence can be used to estimate underlying quantitative tissue parameter maps, which are then used to synthesize any desirable clinical contrast by retrospectively changing scan parameters in silico. Two benefits of this approach are the reduced exam time and the ability to generate arbitrary contrasts offline. However, synthetically generated contrasts are known to deviate from the contrast of experimental scans. The reason for contrast mismatch is the necessary exclusion of some unmodeled physical effects such as partial voluming, diffusion, flow, susceptibility, magnetization transfer, and more. The inclusion of these effects in signal encoding would improve the synthetic images, but would make the quantitative imaging protocol impractical due to long scan times. Therefore, in this work, we propose a novel deep learning approach that generates a multiplicative correction term to capture unmodeled effects and correct the synthetic contrast images to better match experimental contrasts for arbitrary scan parameters. The physics inspired deep learning model implicitly accounts for some unmodeled physical effects occurring during the scan. As a proof of principle, we validate our approach on synthesizing arbitrary inversion recovery fast spin-echo scans using a commercially available 2D multi-contrast sequence. We observe that the proposed correction visually and numerically reduces the mismatch with experimentally collected contrasts compared to conventional synthetic MRI. Finally, we show results of a preliminary reader study and find that the proposed method statistically significantly improves in contrast and SNR as compared to synthetic MR images.
Collapse
Affiliation(s)
- Sidharth Kumar
- Chandra Family Department of Electrical and Computer Engineering, The University of Texas at Austin, Austin 78712, TX, USA.
| | - Hamidreza Saber
- Dell Medical School Department of Neurology, The University of Texas at Austin, Austin 78712, TX, USA; Dell Medical School Department of Neurosurgery, The University of Texas at Austin, Austin 78712, TX, USA
| | - Odelin Charron
- Dell Medical School Department of Neurology, The University of Texas at Austin, Austin 78712, TX, USA
| | - Leorah Freeman
- Dell Medical School Department of Neurology, The University of Texas at Austin, Austin 78712, TX, USA; Dell Medical School Department of Diagnostic Medicine, The University of Texas at Austin, Austin 78712, TX, USA
| | - Jonathan I Tamir
- Chandra Family Department of Electrical and Computer Engineering, The University of Texas at Austin, Austin 78712, TX, USA; Dell Medical School Department of Diagnostic Medicine, The University of Texas at Austin, Austin 78712, TX, USA; Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin 78712, TX, USA
| |
Collapse
|
31
|
Kadalie E, Trotier AJ, Corbin N, Miraux S, Ribot EJ. Rapid whole brain 3D T 2 mapping respiratory-resolved Double-Echo Steady State (DESS) sequence with improved repeatability. Magn Reson Med 2024; 91:221-236. [PMID: 37794821 DOI: 10.1002/mrm.29847] [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: 03/03/2023] [Revised: 08/11/2023] [Accepted: 08/11/2023] [Indexed: 10/06/2023]
Abstract
PURPOSE To propose a quantitative 3D double-echo steady-state (DESS) sequence that offers rapid and repeatable T2 mapping of the human brain using different encoding schemes that account for respiratory B0 variation. METHODS A retrospective self-gating module was firstly implemented into the standard DESS sequence in order to suppress the respiratory artifact via data binning. A compressed-sensing trajectory (CS-DESS) was then optimized to accelerate the acquisition. Finally, a spiral Cartesian encoding (SPICCS-DESS) was incorporated to further disrupt the coherent respiratory artifact. These different versions were compared to a standard DESS sequence (fully DESS) by assessing the T2 distribution and repeatability in different brain regions of eight volunteers at 3 T. RESULTS The respiratory artifact correction was determined to be optimal when the data was binned into seven respiratory phases. Compared to the fully DESS, T2 distribution was improved for the CS-DESS and SPICCS-DESS with interquartile ranges reduced significantly by a factor ranging from 2 to 12 in the caudate, putamen, and thalamus regions. In the gray and white matter areas, average absolute test-retest T2 differences across all volunteers were respectively 3.5 ± 2% and 3.1 ± 2.1% for the SPICCS-DESS, 4.6 ± 4.6% and 4.9 ± 5.1% for the CS-DESS, and 15% ± 13% and 7.3 ± 5.6% for the fully DESS. The SPICCS-DESS sequence's acquisition time could be reduced by half (<4 min) while maintaining its efficient T2 mapping. CONCLUSION The respiratory-resolved SPICCS-DESS sequence offers rapid, robust, and repeatable 3D T2 mapping of the human brain, which can be especially effective for longitudinal monitoring of cerebral pathologies.
Collapse
Affiliation(s)
- Emile Kadalie
- Univ. Bordeaux, CNRS, Centre de Résonance Magnétique des Systèmes Biologiques (CRMSB), UMR 5536, F-33000, Bordeaux, France
| | - Aurélien J Trotier
- Univ. Bordeaux, CNRS, Centre de Résonance Magnétique des Systèmes Biologiques (CRMSB), UMR 5536, F-33000, Bordeaux, France
| | - Nadège Corbin
- Univ. Bordeaux, CNRS, Centre de Résonance Magnétique des Systèmes Biologiques (CRMSB), UMR 5536, F-33000, Bordeaux, France
| | - Sylvain Miraux
- Univ. Bordeaux, CNRS, Centre de Résonance Magnétique des Systèmes Biologiques (CRMSB), UMR 5536, F-33000, Bordeaux, France
| | - Emeline J Ribot
- Univ. Bordeaux, CNRS, Centre de Résonance Magnétique des Systèmes Biologiques (CRMSB), UMR 5536, F-33000, Bordeaux, France
| |
Collapse
|
32
|
Mickevicius NJ. Magnetic resonance coherence pathway unraveling. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2024; 358:107613. [PMID: 38134509 DOI: 10.1016/j.jmr.2023.107613] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/11/2023] [Revised: 12/08/2023] [Accepted: 12/15/2023] [Indexed: 12/24/2023]
Abstract
Efficiently acquiring multi-contrast magnetic resonance imaging data is crucial for patient comfort and clinical throughput. Developing scan acceleration methods tailored for specific applications drastically improves the value of an MRI examination. Here, we propose a novel method to control the aliasing of simultaneously acquired images of multiple spin echo coherence pathways with the goal of producing high quality multi-contrast images from a single acquisition. Modulating the radiofrequency phase of several pulses applied in brief succession also uniquely modulates the phase of spin echo coherence pathways. A method, termed magnetic resonance coherence pathway unraveling (MR-CPU), to control the aliasing of simultaneously acquired coherence pathway images is developed here along with parallel imaging-based reconstruction methods to separate them. MR-CPU was validated in phantom experiments and tested in vivo. High levels of correlation between reference pathway images and MR-CPU-derived coherence pathway images were found from the phantom experiments. Minimal artifacts arising from the separation of the overlapped coherence pathway images were observed in vivo. MR-CPU provides a novel mechanism through which to acquire and separate multiple overlapped coherence pathway images, thus adding to the diagnostic potential of an MRI exam without the penalty of additional scan time.
Collapse
|
33
|
Schäper J, Bauman G, Bieri O. Improved gray-white matter contrast using magnetization prepared fast imaging with steady-state free precession (MP-FISP) brain imaging at 0.55 T. Magn Reson Med 2024; 91:162-173. [PMID: 37598421 DOI: 10.1002/mrm.29838] [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: 03/16/2023] [Revised: 07/10/2023] [Accepted: 08/02/2023] [Indexed: 08/22/2023]
Abstract
PURPOSE To improve the gray/white matter contrast of magnetization prepared rapid gradient echo (MP-RAGE) MRI at 0.55 T by optimizing the acquisition and sequence kernel parameters. METHODS A segmented magnetization prepared rapid gradient echo prototype sequence was implemented with (MP-RAGE*) and without (MP-FISP*) radiofrequency spoiling. Optimized parameters were derived with the assistance of an extended phase graph signal simulation as a function of the relaxation times, the flip angle, the delay times, and the effective inversion time using segmentation. The resulting protocols were compared to the MP-RAGE product sequence offered by the vendor in terms of signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR). A tissue segmentation reproducibility study was performed on three volunteers for the product MP-RAGE and the MP-FISP*. RESULTS The MP-RAGE simulation reproduced the parameters already used in the product MP-RAGE on the scanner. An average CNR improvement of 15% for the custom MP-RAGE* over the product MP-RAGE and additional 22% for the MP-FISP* over the MP-RAGE* were observed, which is in accordance with the simulation results. The total improvement, averaged over all volunteers and regions, was 41%. The reproducibility study did not yield a significant difference between MP-RAGE and MP-FISP*. CONCLUSION We presented some easy-to-implement adjustments to the MP-RAGE sequence at 0.55 T, which can lead to an overall average improvement of 41% in CNR.
Collapse
Affiliation(s)
- Jessica Schäper
- Department of Biomedical Engineering, University of Basel, Basel, Switzerland
- Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland
| | - Grzegorz Bauman
- Department of Biomedical Engineering, University of Basel, Basel, Switzerland
- Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland
| | - Oliver Bieri
- Department of Biomedical Engineering, University of Basel, Basel, Switzerland
- Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland
| |
Collapse
|
34
|
Cabini RF, Barzaghi L, Cicolari D, Arosio P, Carrazza S, Figini S, Filibian M, Gazzano A, Krause R, Mariani M, Peviani M, Pichiecchio A, Pizzagalli DU, Lascialfari A. Fast deep learning reconstruction techniques for preclinical magnetic resonance fingerprinting. NMR IN BIOMEDICINE 2024; 37:e5028. [PMID: 37669779 DOI: 10.1002/nbm.5028] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/18/2023] [Revised: 07/05/2023] [Accepted: 07/27/2023] [Indexed: 09/07/2023]
Abstract
We propose a deep learning (DL) model and a hyperparameter optimization strategy to reconstruct T1 and T2 maps acquired with the magnetic resonance fingerprinting (MRF) methodology. We applied two different MRF sequence routines to acquire images of ex vivo rat brain phantoms using a 7-T preclinical scanner. Subsequently, the DL model was trained using experimental data, completely excluding the use of any theoretical MRI signal simulator. The best combination of the DL parameters was implemented by an automatic hyperparameter optimization strategy, whose key aspect is to include all the parameters to the fit, allowing the simultaneous optimization of the neural network architecture, the structure of the DL model, and the supervised learning algorithm. By comparing the reconstruction performances of the DL technique with those achieved from the traditional dictionary-based method on an independent dataset, the DL approach was shown to reduce the mean percentage relative error by a factor of 3 for T1 and by a factor of 2 for T2 , and to improve the computational time by at least a factor of 37. Furthermore, the proposed DL method enables maintaining comparable reconstruction performance, even with a lower number of MRF images and a reduced k-space sampling percentage, with respect to the dictionary-based method. Our results suggest that the proposed DL methodology may offer an improvement in reconstruction accuracy, as well as speeding up MRF for preclinical, and in prospective clinical, investigations.
Collapse
Affiliation(s)
- Raffaella Fiamma Cabini
- Department of Mathematics, University of Pavia, Pavia, Italy
- INFN, Istituto Nazionale di Fisica Nucleare, Pavia, Italy
| | - Leonardo Barzaghi
- Department of Mathematics, University of Pavia, Pavia, Italy
- INFN, Istituto Nazionale di Fisica Nucleare, Pavia, Italy
- Advanced Imaging and Artificial Intelligence, Department of Neuroradiology, IRCCS Mondino Foundation, Pavia, Italy
| | - Davide Cicolari
- Department of Physics, University of Pavia, Pavia, Italy
- Department of Physics, University of Milan, Milan, Italy
- INFN, Istituto Nazionale di Fisica Nucleare, Milan, Italy
- Department of Medical Physics, ASST GOM Niguarda, Milan, Italy
| | - Paolo Arosio
- Department of Physics, University of Milan, Milan, Italy
- INFN, Istituto Nazionale di Fisica Nucleare, Milan, Italy
| | - Stefano Carrazza
- Department of Physics, University of Milan, Milan, Italy
- INFN, Istituto Nazionale di Fisica Nucleare, Milan, Italy
| | - Silvia Figini
- INFN, Istituto Nazionale di Fisica Nucleare, Pavia, Italy
- Department of Social and Political Science, University of Pavia, Pavia, Italy
| | - Marta Filibian
- INFN, Istituto Nazionale di Fisica Nucleare, Pavia, Italy
- Centro Grandi Strumenti, University of Pavia, Pavia, Italy
| | - Andrea Gazzano
- Laboratory of Cellular and Molecular Neuropharmacology, Department of Biology and Biotechnology "L. Spallanzani", University of Pavia, Pavia, Italy
| | | | - Manuel Mariani
- Department of Physics, University of Pavia, Pavia, Italy
| | - Marco Peviani
- Laboratory of Cellular and Molecular Neuropharmacology, Department of Biology and Biotechnology "L. Spallanzani", University of Pavia, Pavia, Italy
| | - Anna Pichiecchio
- Advanced Imaging and Artificial Intelligence, Department of Neuroradiology, IRCCS Mondino Foundation, Pavia, Italy
- Department of Brain and Behavioural Sciences, University of Pavia, Pavia, Italy
| | | | - Alessandro Lascialfari
- INFN, Istituto Nazionale di Fisica Nucleare, Pavia, Italy
- Department of Physics, University of Pavia, Pavia, Italy
| |
Collapse
|
35
|
Zhang M, Arango N, Arefeen Y, Guryev G, Stockmann JP, White J, Adalsteinsson E. Stochastic-offset-enhanced restricted slice excitation and 180° refocusing designs with spatially non-linear ΔB 0 shim array fields. Magn Reson Med 2023; 90:2572-2591. [PMID: 37667645 PMCID: PMC10699120 DOI: 10.1002/mrm.29827] [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: 01/19/2023] [Revised: 06/30/2023] [Accepted: 07/26/2023] [Indexed: 09/06/2023]
Abstract
PURPOSE Developing a general framework with a novel stochastic offset strategy for the design of optimized RF pulses and time-varying spatially non-linear ΔB0 shim array fields for restricted slice excitation and refocusing with refined magnetization profiles within the intervals of the fixed voxels. METHODS Our framework uses the decomposition property of the Bloch equations to enable joint design of RF-pulses and shim array fields for restricted slice excitation and refocusing with auto-differentiation optimization. Bloch simulations are performed independently on orthogonal basis vectors, Mx, My, and Mz, which enables designs for arbitrary initial magnetizations. Requirements for refocusing pulse designs are derived from the extended phase graph formalism obviating time-consuming sub-voxel isochromatic simulations to model the effects of crusher gradients. To refine resultant slice-profiles because of voxelwise optimization functions, we propose an algorithm that stochastically offsets spatial points at which loss is computed during optimization. RESULTS We first applied our proposed design framework to standard slice-selective excitation and refocusing pulses in the absence of non-linear ΔB0 shim array fields and compared them against pulses designed with Shinnar-Le Roux algorithm. Next, we demonstrated our technique in a simulated setup of fetal brain imaging in pregnancy for restricted-slice excitation and refocusing of the fetal brain. CONCLUSIONS Our proposed framework for optimizing RF pulse and time-varying spatially non-linear ΔB0 shim array fields achieve high fidelity restricted-slice excitation and refocusing for fetal MRI, which could enable zoomed fast-spin-echo-MRI and other applications.
Collapse
Affiliation(s)
- Molin Zhang
- Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Nicolas Arango
- Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Yamin Arefeen
- Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Georgy Guryev
- Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Jason P. Stockmann
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, USA
- Department of Radiology, Harvard Medical School, Charlestown, MA, USA
| | - Jacob White
- Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Elfar Adalsteinsson
- Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, MA, USA
- Harvard-MIT Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, USA
- Institute for Medical Engineering and Science, Massachusetts Institute of Technology, Cambridge, MA, USA
| |
Collapse
|
36
|
Wang K, Doneva M, Meineke J, Amthor T, Karasan E, Tan F, Tamir JI, Yu SX, Lustig M. High-fidelity direct contrast synthesis from magnetic resonance fingerprinting. Magn Reson Med 2023; 90:2116-2129. [PMID: 37332200 DOI: 10.1002/mrm.29766] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Revised: 05/03/2023] [Accepted: 05/31/2023] [Indexed: 06/20/2023]
Abstract
PURPOSE This work was aimed at proposing a supervised learning-based method that directly synthesizes contrast-weighted images from the Magnetic Resonance Fingerprinting (MRF) data without performing quantitative mapping and spin-dynamics simulations. METHODS To implement our direct contrast synthesis (DCS) method, we deploy a conditional generative adversarial network (GAN) framework with a multi-branch U-Net as the generator and a multilayer CNN (PatchGAN) as the discriminator. We refer to our proposed approach as N-DCSNet. The input MRF data are used to directly synthesize T1-weighted, T2-weighted, and fluid-attenuated inversion recovery (FLAIR) images through supervised training on paired MRF and target spin echo-based contrast-weighted scans. The performance of our proposed method is demonstrated on in vivo MRF scans from healthy volunteers. Quantitative metrics, including normalized root mean square error (nRMSE), peak signal-to-noise ratio (PSNR), structural similarity (SSIM), learned perceptual image patch similarity (LPIPS), and Fréchet inception distance (FID), were used to evaluate the performance of the proposed method and compare it with others. RESULTS In-vivo experiments demonstrated excellent image quality with respect to that of simulation-based contrast synthesis and previous DCS methods, both visually and according to quantitative metrics. We also demonstrate cases in which our trained model is able to mitigate the in-flow and spiral off-resonance artifacts typically seen in MRF reconstructions, and thus more faithfully represent conventional spin echo-based contrast-weighted images. CONCLUSION We present N-DCSNet to directly synthesize high-fidelity multicontrast MR images from a single MRF acquisition. This method can significantly decrease examination time. By directly training a network to generate contrast-weighted images, our method does not require any model-based simulation and therefore can avoid reconstruction errors due to dictionary matching and contrast simulation (code available at:https://github.com/mikgroup/DCSNet).
Collapse
Affiliation(s)
- Ke Wang
- Electrical Engineering and Computer Sciences, University of California at Berkeley, Berkeley, California, USA
- International Computer Science Institute, University of California at Berkeley, Berkeley, California, USA
| | | | | | | | - Ekin Karasan
- Electrical Engineering and Computer Sciences, University of California at Berkeley, Berkeley, California, USA
| | - Fei Tan
- Bioengineering, UC Berkeley-UCSF, San Francisco, California, USA
| | - Jonathan I Tamir
- Chandra Family Department of Electrical and Computer Engineering, The University of Texas at Austin, Austin, Texas, USA
| | - Stella X Yu
- Electrical Engineering and Computer Sciences, University of California at Berkeley, Berkeley, California, USA
- International Computer Science Institute, University of California at Berkeley, Berkeley, California, USA
- Computer Science and Engineering, University of Michigan, Ann Arbor, Michigan, USA
| | | |
Collapse
|
37
|
Barbieri M, Watkins LE, Mazzoli V, Desai AD, Rubin E, Schmidt A, Gold GE, Hargreaves BA, Chaudhari AS, Kogan F. [Formula: see text] Field inhomogeneity correction for qDESS [Formula: see text] mapping: application to rapid bilateral knee imaging. MAGMA (NEW YORK, N.Y.) 2023; 36:711-724. [PMID: 37142852 PMCID: PMC10524110 DOI: 10.1007/s10334-023-01094-y] [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/23/2022] [Revised: 04/17/2023] [Accepted: 04/18/2023] [Indexed: 05/06/2023]
Abstract
PURPOSE [Formula: see text] mapping is a powerful tool for studying osteoarthritis (OA) changes and bilateral imaging may be useful in investigating the role of between-knee asymmetry in OA onset and progression. The quantitative double-echo in steady-state (qDESS) can provide fast simultaneous bilateral knee [Formula: see text] and high-resolution morphometry for cartilage and meniscus. The qDESS uses an analytical signal model to compute [Formula: see text] relaxometry maps, which require knowledge of the flip angle (FA). In the presence of [Formula: see text] inhomogeneities, inconsistencies between the nominal and actual FA can affect the accuracy of [Formula: see text] measurements. We propose a pixel-wise [Formula: see text] correction method for qDESS [Formula: see text] mapping exploiting an auxiliary [Formula: see text] map to compute the actual FA used in the model. METHODS The technique was validated in a phantom and in vivo with simultaneous bilateral knee imaging. [Formula: see text] measurements of femoral cartilage (FC) of both knees of six healthy participants were repeated longitudinally to investigate the association between [Formula: see text] variation and [Formula: see text]. RESULTS The results showed that applying the [Formula: see text] correction mitigated [Formula: see text] variations that were driven by [Formula: see text] inhomogeneities. Specifically, [Formula: see text] left-right symmetry increased following the [Formula: see text] correction ([Formula: see text] = 0.74 > [Formula: see text] = 0.69). Without the [Formula: see text] correction, [Formula: see text] values showed a linear dependence with [Formula: see text]. The linear coefficient decreased using the [Formula: see text] correction (from 24.3 ± 1.6 ms to 4.1 ± 1.8) and the correlation was not statistically significant after the application of the Bonferroni correction (p value > 0.01). CONCLUSION The study showed that [Formula: see text] correction could mitigate variations driven by the sensitivity of the qDESS [Formula: see text] mapping method to [Formula: see text], therefore, increasing the sensitivity to detect real biological changes. The proposed method may improve the robustness of bilateral qDESS [Formula: see text] mapping, allowing for an accurate and more efficient evaluation of OA pathways and pathophysiology through longitudinal and cross-sectional studies.
Collapse
Affiliation(s)
- Marco Barbieri
- Department of Radiology, Stanford University, Stanford, CA, USA
| | - Lauren E. Watkins
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Bioengineering, Stanford University, Stanford, CA, USA
| | | | - Arjun D. Desai
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Elka Rubin
- Department of Radiology, Stanford University, Stanford, CA, USA
| | - Andrew Schmidt
- Department of Radiology, Stanford University, Stanford, CA, USA
| | - Garry Evan Gold
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Bioengineering, Stanford University, Stanford, CA, USA
| | - Brian Andrew Hargreaves
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Bioengineering, Stanford University, Stanford, CA, USA
- Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Akshay Sanjay Chaudhari
- Department of Radiology, Stanford University, Stanford, CA, USA
- Department of Biomedical Data Science, Stanford University, Stanford, CA, USA
| | - Feliks Kogan
- Department of Radiology, Stanford University, Stanford, CA, USA
| |
Collapse
|
38
|
Qiu S, Ma S, Wang L, Chen Y, Fan Z, Moser FG, Maya M, Sati P, Sicotte NL, Christodoulou AG, Xie Y, Li D. Direct synthesis of multi-contrast brain MR images from MR multitasking spatial factors using deep learning. Magn Reson Med 2023; 90:1672-1681. [PMID: 37246485 PMCID: PMC10524469 DOI: 10.1002/mrm.29715] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2022] [Revised: 04/27/2023] [Accepted: 05/03/2023] [Indexed: 05/30/2023]
Abstract
PURPOSE To develop a deep learning method to synthesize conventional contrast-weighted images in the brain from MR multitasking spatial factors. METHODS Eighteen subjects were imaged using a whole-brain quantitative T1 -T2 -T1ρ MR multitasking sequence. Conventional contrast-weighted images consisting of T1 MPRAGE, T1 gradient echo, and T2 fluid-attenuated inversion recovery were acquired as target images. A 2D U-Net-based neural network was trained to synthesize conventional weighted images from MR multitasking spatial factors. Quantitative assessment and image quality rating by two radiologists were performed to evaluate the quality of deep-learning-based synthesis, in comparison with Bloch-equation-based synthesis from MR multitasking quantitative maps. RESULTS The deep-learning synthetic images showed comparable contrasts of brain tissues with the reference images from true acquisitions and were substantially better than the Bloch-equation-based synthesis results. Averaging on the three contrasts, the deep learning synthesis achieved normalized root mean square error = 0.184 ± 0.075, peak SNR = 28.14 ± 2.51, and structural-similarity index = 0.918 ± 0.034, which were significantly better than Bloch-equation-based synthesis (p < 0.05). Radiologists' rating results show that compared with true acquisitions, deep learning synthesis had no notable quality degradation and was better than Bloch-equation-based synthesis. CONCLUSION A deep learning technique was developed to synthesize conventional weighted images from MR multitasking spatial factors in the brain, enabling the simultaneous acquisition of multiparametric quantitative maps and clinical contrast-weighted images in a single scan.
Collapse
Affiliation(s)
- Shihan Qiu
- Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, California, USA
- Department of Bioengineering, UCLA, Los Angeles, California, USA
| | - Sen Ma
- Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, California, USA
| | - Lixia Wang
- Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, California, USA
| | - Yuhua Chen
- Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, California, USA
- Department of Bioengineering, UCLA, Los Angeles, California, USA
| | - Zhaoyang Fan
- Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, California, USA
- Departments of Radiology and Radiation Oncology, University of Southern California, Los Angeles, California, USA
| | - Franklin G. Moser
- Department of Imaging, Cedars-Sinai Medical Center, Los Angeles, California, USA
| | - Marcel Maya
- Department of Imaging, Cedars-Sinai Medical Center, Los Angeles, California, USA
| | - Pascal Sati
- Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, California, USA
- Department of Neurology, Cedars-Sinai Medical Center, Los Angeles, California, USA
| | - Nancy L. Sicotte
- Department of Neurology, Cedars-Sinai Medical Center, Los Angeles, California, USA
| | - Anthony G. Christodoulou
- Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, California, USA
- Department of Bioengineering, UCLA, Los Angeles, California, USA
| | - Yibin Xie
- Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, California, USA
| | - Debiao Li
- Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, California, USA
- Department of Bioengineering, UCLA, Los Angeles, California, USA
| |
Collapse
|
39
|
Cohen O, Otazo R. Global deep learning optimization of chemical exchange saturation transfer magnetic resonance fingerprinting acquisition schedule. NMR IN BIOMEDICINE 2023; 36:e4954. [PMID: 37070221 PMCID: PMC10896067 DOI: 10.1002/nbm.4954] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/05/2023] [Revised: 04/10/2023] [Accepted: 04/11/2023] [Indexed: 05/06/2023]
Abstract
Chemical exchange saturation transfer (CEST) MRI is a promising molecular imaging technique but suffers from long scan times and complicated processing. CEST was recently combined with magnetic resonance fingerprinting (MRF) to address these shortcomings. However, the CEST-MRF signal depends on multiple acquisition and tissue parameters so selecting an optimal acquisition schedule is challenging. In this work, we propose a novel dual-network deep learning framework to optimize the CEST-MRF acquisition schedule. The quality of the optimized schedule was assessed in a digital brain phantom and compared with alternate deep learning optimization approaches. The effect of schedule length on the reconstruction error was also investigated. A healthy subject was scanned with optimized and random schedules and with a conventional CEST sequence for comparison. The optimized schedule was also tested in a subject with metastatic renal cell carcinoma. Reproducibility was assessed via test-retest experiments and the concordance correlation coefficient calculated for white matter (WM) and grey matter (GM). The optimized schedule was 12% shorter but yielded equal or lower normalized root mean square error for all parameters. The proposed optimization also provided a lower error compared with alternate methodologies. Longer schedules generally yielded lower error. In vivo maps obtained with the optimized schedule showed reduced noise and improved delineation of GM and WM. CEST curves synthesized from the optimized parameters were highly correlated (r = 0.99) with measured conventional CEST. The mean concordance correlation coefficient in WM/GM for all tissue parameters was 0.990/0.978 for the optimized schedule but only 0.979/0.975 for the random schedule. The proposed schedule optimization is widely applicable to MRF pulse sequences and provides accurate and reproducible tissue maps with reduced noise at a shorter scan time than a randomly generated schedule.
Collapse
Affiliation(s)
- Ouri Cohen
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, New York, USA
| | - Ricardo Otazo
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, New York, USA
- Department of Radiology, Memorial Sloan Kettering Cancer Center, New York, New York, USA
| |
Collapse
|
40
|
Liu C, Li T, Cao P, Hui ES, Wong YL, Wang Z, Xiao H, Zhi S, Zhou T, Li W, Lam SK, Cheung ALY, Lee VHF, Ying M, Cai J. Respiratory-Correlated 4-Dimensional Magnetic Resonance Fingerprinting for Liver Cancer Radiation Therapy Motion Management. Int J Radiat Oncol Biol Phys 2023; 117:493-504. [PMID: 37116591 DOI: 10.1016/j.ijrobp.2023.04.015] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2022] [Revised: 04/04/2023] [Accepted: 04/18/2023] [Indexed: 04/30/2023]
Abstract
PURPOSE The objective of this study was to develop a respiratory-correlated (RC) 4-dimensional (4D) imaging technique based on magnetic resonance fingerprinting (MRF) (RC-4DMRF) for liver tumor motion management in radiation therapy. METHODS AND MATERIALS Thirteen patients with liver cancer were prospectively enrolled in this study. k-space MRF signals of the liver were acquired during free-breathing using the fast acquisition with steady-state precession sequence on a 3T scanner. The signals were binned into 8 respiratory phases based on respiratory surrogates, and interphase displacement vector fields were estimated using a phase-specific low-rank optimization method. Hereafter, the tissue property maps, including T1 and T2 relaxation times, and proton density, were reconstructed using a pyramid motion-compensated method that alternatively optimized interphase displacement vector fields and subspace images. To evaluate the efficacy of RC-4DMRF, amplitude motion differences and Pearson correlation coefficients were determined to assess measurement agreement in tumor motion between RC-4DMRF and cine magnetic resonance imaging (MRI); mean absolute percentage errors of the RC-4DMRF-derived tissue maps were calculated to reveal tissue quantification accuracy using digital human phantom; and tumor-to-liver contrast-to-noise ratio of RC-4DMRF images was compared with that of planning CT and contrast-enhanced MRI (CE-MRI) images. A paired Student t test was used for statistical significance analysis with a P value threshold of .05. RESULTS RC-4DMRF achieved excellent agreement in motion measurement with cine MRI, yielding the mean (± standard deviation) Pearson correlation coefficients of 0.95 ± 0.05 and 0.93 ± 0.09 and amplitude motion differences of 1.48 ± 1.06 mm and 0.81 ± 0.64 mm in the superior-inferior and anterior-posterior directions, respectively. Moreover, RC-4DMRF achieved high accuracy in tissue property quantification, with mean absolute percentage errors of 8.8%, 9.6%, and 5.0% for T1, T2, and proton density, respectively. Notably, the tumor contrast-to-noise ratio in RC-4DMRI-derived T1 maps (6.41 ± 3.37) was found to be the highest among all tissue property maps, approximately equal to that of CE-MRI (6.96 ± 1.01, P = .862), and substantially higher than that of planning CT (2.91 ± 1.97, P = .048). CONCLUSIONS RC-4DMRF demonstrated high accuracy in respiratory motion measurement and tissue properties quantification, potentially facilitating tumor motion management in liver radiation therapy.
Collapse
Affiliation(s)
- Chenyang Liu
- Department of Health Technology and Informatics, Hong Kong Polytechnic University, Hong Kong SAR, China
| | - Tian Li
- Department of Health Technology and Informatics, Hong Kong Polytechnic University, Hong Kong SAR, China
| | - Peng Cao
- Department of Diagnostic Radiology, University of Hong Kong, Hong Kong SAR, China
| | - Edward S Hui
- Department of Imaging and Interventional Radiology, Chinese University of Hong Kong, Hong Kong SAR, China; Department of Psychiatry, Chinese University of Hong Kong, Hong Kong SAR, China
| | - Yat-Lam Wong
- Department of Health Technology and Informatics, Hong Kong Polytechnic University, Hong Kong SAR, China
| | - Zuojun Wang
- Department of Diagnostic Radiology, University of Hong Kong, Hong Kong SAR, China
| | - Haonan Xiao
- Department of Health Technology and Informatics, Hong Kong Polytechnic University, Hong Kong SAR, China
| | - Shaohua Zhi
- Department of Health Technology and Informatics, Hong Kong Polytechnic University, Hong Kong SAR, China
| | - Ta Zhou
- Department of Health Technology and Informatics, Hong Kong Polytechnic University, Hong Kong SAR, China
| | - Wen Li
- Department of Health Technology and Informatics, Hong Kong Polytechnic University, Hong Kong SAR, China
| | - Sai Kit Lam
- Department of Biomedical Engineering, Hong Kong Polytechnic University, Hong Kong SAR, China; Research Institute for Smart Ageing, Hong Kong Polytechnic University, Hong Kong SAR, China
| | | | - Victor Ho-Fun Lee
- Department of Clinical Oncology, University of Hong Kong, Hong Kong SAR, China
| | - Michael Ying
- Department of Health Technology and Informatics, Hong Kong Polytechnic University, Hong Kong SAR, China.
| | - Jing Cai
- Department of Health Technology and Informatics, Hong Kong Polytechnic University, Hong Kong SAR, China; Research Institute for Smart Ageing, Hong Kong Polytechnic University, Hong Kong SAR, China.
| |
Collapse
|
41
|
Zhao C, Shao X, Shou Q, Ma SJ, Gokyar S, Graf C, Stollberger R, Wang DJ. Whole-Cerebrum distortion-free three-dimensional pseudo-continuous arterial spin labeling at 7T. Neuroimage 2023; 277:120251. [PMID: 37364741 PMCID: PMC10528743 DOI: 10.1016/j.neuroimage.2023.120251] [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/24/2023] [Revised: 06/21/2023] [Accepted: 06/23/2023] [Indexed: 06/28/2023] Open
Abstract
Fulfilling potentials of ultrahigh field for pseudo-Continuous Arterial Spin Labeling (pCASL) has been hampered by B1/B0 inhomogeneities that affect pCASL labeling, background suppression (BS), and the readout sequence. This study aimed to present a whole-cerebrum distortion-free three-dimensional (3D) pCASL sequence at 7T by optimizing pCASL labeling parameters, BS pulses, and an accelerated Turbo-FLASH (TFL) readout. A new set of pCASL labeling parameters (Gave = 0.4 mT/m, Gratio = 14.67) was proposed to avoid interferences in bottom slices while achieving robust labeling efficiency (LE). An OPTIM BS pulse was designed based on the range of B1/B0 inhomogeneities at 7T. A 3D TFL readout with 2D-CAIPIRINHA undersampling (R = 2 × 2) and centric ordering was developed, and the number of segments (Nseg) and flip angle (FA) were varied in simulation to achieve the optimal trade-off between SNR and spatial blurring. In-vivo experiments were performed on 19 subjects. The results showed that the new set of labeling parameters effectively achieved whole-cerebrum coverage by eliminating interferences in bottom slices while maintaining a high LE. The OPTIM BS pulse achieved 33.3% higher perfusion signal in gray matter (GM) than the original BS pulse with a cost of 4.8-fold SAR. Incorporating a moderate FA (8°) and Nseg (2), whole-cerebrum 3D TFL-pCASL imaging was achieved with a 2 × 2 × 4 mm3 resolution without distortion and susceptibility artifacts compared to 3D GRASE-pCASL. In addition, 3D TFL-pCASL showed a good to excellent test-retest repeatability and potential of higher resolution (2 mm isotropic). The proposed technique also significantly improved SNR when compared to the same sequence at 3T and simultaneous multislice TFL-pCASL at 7T. By combining a new set of labeling parameters, OPTIM BS pulse, and accelerated 3D TFL readout, we achieved high resolution pCASL at 7T with whole-cerebrum coverage, detailed perfusion and anatomical information without distortion, and sufficient SNR.
Collapse
Affiliation(s)
- Chenyang Zhao
- Laboratory of FMRI Technology (LOFT), Mark & Mary Stevens Neuroimaging and Informatics Institute, Keck School of Medicine, University of Southern California, United States
| | - Xingfeng Shao
- Laboratory of FMRI Technology (LOFT), Mark & Mary Stevens Neuroimaging and Informatics Institute, Keck School of Medicine, University of Southern California, United States
| | - Qinyang Shou
- Laboratory of FMRI Technology (LOFT), Mark & Mary Stevens Neuroimaging and Informatics Institute, Keck School of Medicine, University of Southern California, United States
| | - Samantha J Ma
- Siemens Medical Solutions USA, Los Angeles, CA, United States
| | - Sayim Gokyar
- Laboratory of FMRI Technology (LOFT), Mark & Mary Stevens Neuroimaging and Informatics Institute, Keck School of Medicine, University of Southern California, United States
| | - Christina Graf
- Institute of Biomedical Imaging, Graz University of Technology, Austria
| | | | - Danny Jj Wang
- Laboratory of FMRI Technology (LOFT), Mark & Mary Stevens Neuroimaging and Informatics Institute, Keck School of Medicine, University of Southern California, United States; Department of Neurology, Keck School of Medicine, University of Southern California, United States.
| |
Collapse
|
42
|
Gaur S, Panda A, Fajardo JE, Hamilton J, Jiang Y, Gulani V. Magnetic Resonance Fingerprinting: A Review of Clinical Applications. Invest Radiol 2023; 58:561-577. [PMID: 37026802 PMCID: PMC10330487 DOI: 10.1097/rli.0000000000000975] [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] [Indexed: 04/08/2023]
Abstract
ABSTRACT Magnetic resonance fingerprinting (MRF) is an approach to quantitative magnetic resonance imaging that allows for efficient simultaneous measurements of multiple tissue properties, which are then used to create accurate and reproducible quantitative maps of these properties. As the technique has gained popularity, the extent of preclinical and clinical applications has vastly increased. The goal of this review is to provide an overview of currently investigated preclinical and clinical applications of MRF, as well as future directions. Topics covered include MRF in neuroimaging, neurovascular, prostate, liver, kidney, breast, abdominal quantitative imaging, cardiac, and musculoskeletal applications.
Collapse
Affiliation(s)
- Sonia Gaur
- Department of Radiology, Michigan Medicine, Ann Arbor, MI
| | - Ananya Panda
- All India Institute of Medical Sciences, Jodhpur, Rajasthan, India
| | | | - Jesse Hamilton
- Department of Radiology, Michigan Medicine, Ann Arbor, MI
| | - Yun Jiang
- Department of Radiology, Michigan Medicine, Ann Arbor, MI
| | - Vikas Gulani
- Department of Radiology, Michigan Medicine, Ann Arbor, MI
| |
Collapse
|
43
|
Wang Z, Ramasawmy R, Feng X, Campbell-Washburn AE, Mugler JP, Meyer CH. Concomitant magnetic-field compensation for 2D spiral-ring turbo spin-echo imaging at 0.55T and 1.5T. Magn Reson Med 2023; 90:552-568. [PMID: 37036033 PMCID: PMC10578525 DOI: 10.1002/mrm.29663] [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: 12/23/2022] [Revised: 03/08/2023] [Accepted: 03/21/2023] [Indexed: 04/11/2023]
Abstract
PURPOSE To develop 2D turbo spin-echo (TSE) imaging using annular spiral rings (abbreviated "SPRING-RIO TSE") with compensation of concomitant gradient fields and B0 inhomogeneity at both 0.55T and 1.5T for fast T2 -weighted imaging. METHODS Strategies of gradient waveform modifications were implemented in SPRING-RIO TSE for compensation of self-squared concomitant gradient terms at the TE and across echo spacings, along with reconstruction-based corrections to simultaneously compensate for the residual concomitant gradient and B0 field induced phase accruals along the readout. The signal pathway disturbance caused by time-varying and spatially dependent concomitant fields was simulated, and echo-to-echo phase variations before and after sequence-based compensation were compared. Images from SPRING-RIO TSE with no compensation, with compensation, and Cartesian TSE were also compared via phantom and in vivo acquisitions. RESULTS Simulation showed how concomitant fields affected the signal evolution with no compensation, and both simulation and phantom studies demonstrated the performance of the proposed sequence modifications, as well as the readout off-resonance corrections. Volunteer data showed that after full correction, the SPRING-RIO TSE sequence achieved high image quality with improved SNR efficiency (15%-20% increase), and reduced RF SAR (˜50% reduction), compared to the standard Cartesian TSE, presenting potential benefits, especially in regaining SNR at low-field (0.55T). CONCLUSION Implementation of SPRING-RIO TSE with concomitant field compensation was tested at 0.55T and 1.5T. The compensation principles can be extended to correct for other trajectory types that are time-varying along the echo train and temporally asymmetric in TSE-based imaging.
Collapse
Affiliation(s)
- Zhixing Wang
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA
| | - Rajiv Ramasawmy
- Cardiovascular Branch, Division of Intramural Research, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, Maryland, USA
| | - Xue Feng
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA
| | - Adrienne E. Campbell-Washburn
- Cardiovascular Branch, Division of Intramural Research, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, Maryland, USA
| | - John P. Mugler
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA
- Department of Radiology & Medical Imaging, University of Virginia, Charlottesville, Virginia, USA
| | - Craig H. Meyer
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA
- Department of Radiology & Medical Imaging, University of Virginia, Charlottesville, Virginia, USA
| |
Collapse
|
44
|
Castillo‐Passi C, Coronado R, Varela‐Mattatall G, Alberola‐López C, Botnar R, Irarrazaval P. KomaMRI.jl: An open-source framework for general MRI simulations with GPU acceleration. Magn Reson Med 2023; 90:329-342. [PMID: 36877139 PMCID: PMC10952765 DOI: 10.1002/mrm.29635] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Revised: 02/20/2023] [Accepted: 02/22/2023] [Indexed: 03/07/2023]
Abstract
PURPOSE To develop an open-source, high-performance, easy-to-use, extensible, cross-platform, and general MRI simulation framework (Koma). METHODS Koma was developed using the Julia programming language. Like other MRI simulators, it solves the Bloch equations with CPU and GPU parallelization. The inputs are the scanner parameters, the phantom, and the pulse sequence that is Pulseq-compatible. The raw data is stored in the ISMRMRD format. For the reconstruction, MRIReco.jl is used. A graphical user interface utilizing web technologies was also designed. Two types of experiments were performed: one to compare the quality of the results and the execution speed, and the second to compare its usability. Finally, the use of Koma in quantitative imaging was demonstrated by simulating Magnetic Resonance Fingerprinting (MRF) acquisitions. RESULTS Koma was compared to two well-known open-source MRI simulators, JEMRIS and MRiLab. Highly accurate results (with mean absolute differences below 0.1% compared to JEMRIS) and better GPU performance than MRiLab were demonstrated. In an experiment with students, Koma was proved to be easy to use, eight times faster on personal computers than JEMRIS, and 65% of test subjects recommended it. The potential for designing acquisition and reconstruction techniques was also shown through the simulation of MRF acquisitions, with conclusions that agree with the literature. CONCLUSIONS Koma's speed and flexibility have the potential to make simulations more accessible for education and research. Koma is expected to be used for designing and testing novel pulse sequences before implementing them in the scanner with Pulseq files, and for creating synthetic data to train machine learning models.
Collapse
Affiliation(s)
- Carlos Castillo‐Passi
- School of Biomedical Engineering and Imaging SciencesKing's College LondonLondonUK
- Institute for Biological and Medical EngineeringPontificia Universidad Católica de ChileSantiagoChile
- Millennium Institute for Intelligent Healthcare Engineering (iHEALTH)Pontificia Universidad Católica de ChileSantiagoChile
| | - Ronal Coronado
- Institute for Biological and Medical EngineeringPontificia Universidad Católica de ChileSantiagoChile
- Millennium Institute for Intelligent Healthcare Engineering (iHEALTH)Pontificia Universidad Católica de ChileSantiagoChile
- Electrical EngineeringPontificia Universidad Católica de ChileSantiagoChile
| | - Gabriel Varela‐Mattatall
- Centre for Functional and Metabolic Mapping (CFMM), Robarts Research InstituteWestern UniversityLondonOntarioCanada
- Department of Medical Biophysics, Schulich School of Medicine and DentistryWestern UniversityLondonOntarioCanada
| | | | - René Botnar
- School of Biomedical Engineering and Imaging SciencesKing's College LondonLondonUK
- Institute for Biological and Medical EngineeringPontificia Universidad Católica de ChileSantiagoChile
- Millennium Institute for Intelligent Healthcare Engineering (iHEALTH)Pontificia Universidad Católica de ChileSantiagoChile
| | - Pablo Irarrazaval
- Institute for Biological and Medical EngineeringPontificia Universidad Católica de ChileSantiagoChile
- Millennium Institute for Intelligent Healthcare Engineering (iHEALTH)Pontificia Universidad Católica de ChileSantiagoChile
- Electrical EngineeringPontificia Universidad Católica de ChileSantiagoChile
- Laboratorio de Procesado de ImagenUniversidad de ValladolidValladolidSpain
| |
Collapse
|
45
|
Cartesian vs radial MR-STAT: An efficiency and robustness study. Magn Reson Imaging 2023; 99:7-19. [PMID: 36709010 DOI: 10.1016/j.mri.2023.01.017] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2022] [Revised: 10/21/2022] [Accepted: 01/14/2023] [Indexed: 01/27/2023]
Abstract
MR Spin TomogrAphy in Time-domain ("MR-STAT") is quantitative MR technique in which multiple quantitative parameters are estimated from a single short scan by solving a large-scale non-linear optimization problem. In this work we extended the MR-STAT framework to non-Cartesian gradient trajectories. Cartesian MR-STAT and radial MR-STAT were compared in terms of time-efficiency and robustness in simulations, gel phantom measurements and in vivo measurements. In simulations, we observed that both Cartesian and radial MR-STAT are highly robust against undersampling. Radial MR-STAT does have a lower spatial encoding power because the outer corners of k-space are never sampled. However, especially in T2, this is compensated by a higher dynamic encoding power that comes from sampling the k-space center with each readout. In gel phantom measurements, Cartesian MR-STAT was observed to be robust against overfitting whereas radial MR-STAT suffered from high-frequency artefacts in the parameter maps at later iterations. These artefacts are hypothesized to be related to hardware imperfections and were (partially) suppressed with image filters. The time-efficiencies were higher for Cartesian MR-STAT in all vials. In-vivo, the radial reconstruction again suffered from overfitting artefacts. The robustness of Cartesian MR-STAT over the entire range of experiments may make it preferable in a clinical setting, despite radial MR-STAT resulting in a higher T1 time-efficiency in white matter.
Collapse
|
46
|
Ostenson J, Robison RK, Brittain EL, Damon BM. Feasibility of joint mapping of triglyceride saturation and water longitudinal relaxation in a single breath hold applied to high fat-fraction adipose depots in the periclavicular anatomy. Magn Reson Imaging 2023; 99:58-66. [PMID: 36764629 PMCID: PMC10088071 DOI: 10.1016/j.mri.2023.02.001] [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/18/2022] [Revised: 01/27/2023] [Accepted: 02/06/2023] [Indexed: 02/11/2023]
Abstract
INTRODUCTION Simultaneous mapping of triglyceride (TAG) saturation and tissue water relaxation may improve the characterization of the structure and function of anatomies with significant adipose tissue. While several groups have demonstrated in vivo TAG saturation imaging using MRI, joint mapping of relaxation and TAG saturation is understudied. Such mappings may avoid bias from physiological motion, if they can be done within a single breath-hold, and also account for static and applied magnetic field heterogeneity. METHODS We propose a transient-state/MR fingerprinting single breath-hold sequence at 3 T, a low-rank reconstruction, and a parameter estimation pipeline that jointly estimates the number of double bonds (NDB), number of methylene interrupted double bonds (NMIDB), and tissue water T1, while accounting for non-ideal radiofrequency transmit scaling and off-resonance effects. We test the proposed method in simulations, in phantom against MR spectroscopy (MRS), and in vivo regions in and around high fat fraction (FF) periclavicular adipose tissue. Partial volume and multi-peak transverse relaxation effects are explored. RESULTS The simulation results demonstrate accurate NDB, NMIDB, and water T1 estimates across a range of NDB, NMIDB, and T1 values. In phantoms, the proposed method's estimates of NDB and NMIDB correlate with those from MR spectroscopy (Pearson correlation ≥0.98), while the water T1 estimates are concordant with a standard phantom. The NDB and NMIDB are sensitive to partial volumes of water, showing increasing bias at FF < 40%. This bias is found to be due to noise and transverse relaxation effects. The in vivo periclavicular adipose tissue has high FF (>90%). The adipose tissue NDB and NMIDB, and muscle T1 estimates are comparable to those reported in the literature. CONCLUSION Robust estimation of NDB, NMIDB at high FF and water T1 across a broad range of FFs are feasible using the proposed methods. Further reduction of noise and model bias are needed to employ the proposed technique in low FF anatomies and pathologies.
Collapse
Affiliation(s)
- Jason Ostenson
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, United States of America; Dept. of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, United States of America.
| | - Ryan K Robison
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, United States of America; Philips, Gainesville, FL, United States of America
| | - Evan L Brittain
- Division of Cardiovascular Medicine, Vanderbilt University Medical Center, Nashville, TN, United States of America
| | - Bruce M Damon
- Dept. of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, United States of America; Carle Clinical Imaging Research Program, Urbana, IL, United States of America; Stephens Family Clinical Research Institute, Carle Health, Urbana, IL, United States of America; Department of Bioengineering and Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, IL, United States of America; Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, United States of America
| |
Collapse
|
47
|
Poojar P, Qian E, Fernandes TT, Nunes RG, Fung M, Quarterman P, Jambawalikar SR, Lignelli A, Geethanath S. Tailored magnetic resonance fingerprinting. Magn Reson Imaging 2023; 99:81-90. [PMID: 36764630 DOI: 10.1016/j.mri.2023.02.002] [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/2021] [Revised: 01/27/2023] [Accepted: 02/06/2023] [Indexed: 02/11/2023]
Abstract
Neuroimaging of certain pathologies requires both multi-parametric qualitative and quantitative imaging. The role of the quantitative MRI (qMRI) is well accepted but suffers from long acquisition times leading to patient discomfort, especially in geriatric and pediatric patients. Previous studies show that synthetic MRI can be used in order to reduce the scan time and provide qMRI as well as multi-contrast data. However, this approach suffers from artifacts such as partial volume and flow. In order to increase the scan efficiency (the number of contrasts and quantitative maps acquired per unit time), we designed, simulated, and demonstrated rapid, simultaneous, multi-contrast qualitative (T1 weighted, T1 fluid attenuated inversion recovery (FLAIR), T2 weighted, water, and fat), and quantitative imaging (T1 and T2 maps) through the approach of tailored MR fingerprinting (TMRF) to cover whole-brain in approximately four minutes. We performed TMRF on in vivo four healthy human brains and in vitro ISMRM/NIST phantom and compared with vendor supplied gold standard (GS) and MRF sequences. All scans were performed on a 3 T GE Premier system and images were reconstructed offline using MATLAB. The reconstructed qualitative images were then subjected to custom DL denoising and gradient anisotropic diffusion denoising. The quantitative tissue parametric maps were reconstructed using a dense neural network to gain computational speed compared to dictionary matching. The grey matter and white matter tissues in qualitative and quantitative data for the in vivo datasets were segmented semi-automatically. The SNR and mean contrasts were plotted and compared across all three methods. The GS images show better SNR in all four subjects compared to MRF and TMRF (GS > TMRF>MRF). The T1 and T2 values of MRF are relatively overestimated as compared to GS and TMRF. The scan efficiency for TMRF is 1.72 min-1 which is higher compared to GS (0.32 min-1) and MRF (0.90 min-1).
Collapse
Affiliation(s)
- Pavan Poojar
- Icahn School of Medicine at Mt. Sinai, New York, NY, USA; Columbia Magnetic Resonance Research Center, Columbia University in the city of New York, NY, USA
| | - Enlin Qian
- Columbia Magnetic Resonance Research Center, Columbia University in the city of New York, NY, USA
| | - Tiago T Fernandes
- Institute for Systems and Robotics and Department of Bioengineering, Instituto Superior Técnico, Universidade de Lisboa, Lisbon, Portugal
| | - Rita G Nunes
- Institute for Systems and Robotics and Department of Bioengineering, Instituto Superior Técnico, Universidade de Lisboa, Lisbon, Portugal
| | - Maggie Fung
- GE Healthcare Applied Sciences Laboratory East, New York, NY, USA
| | | | - Sachin R Jambawalikar
- Department of Radiology, Columbia University Irving Medical Center, Columbia University in the city of New York, NY, USA
| | - Angela Lignelli
- Department of Radiology, Columbia University Irving Medical Center, Columbia University in the city of New York, NY, USA
| | - Sairam Geethanath
- Icahn School of Medicine at Mt. Sinai, New York, NY, USA; Columbia Magnetic Resonance Research Center, Columbia University in the city of New York, NY, USA.
| |
Collapse
|
48
|
O'Reilly T, Börnert P, Liu H, Webb A, Koolstra K. 3D magnetic resonance fingerprinting on a low-field 50 mT point-of-care system prototype: evaluation of muscle and lipid relaxation time mapping and comparison with standard techniques. MAGMA (NEW YORK, N.Y.) 2023:10.1007/s10334-023-01092-0. [PMID: 37202655 PMCID: PMC10386962 DOI: 10.1007/s10334-023-01092-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/09/2022] [Revised: 04/11/2023] [Accepted: 04/17/2023] [Indexed: 05/20/2023]
Abstract
OBJECTIVE To implement magnetic resonance fingerprinting (MRF) on a permanent magnet 50 mT low-field system deployable as a future point-of-care (POC) unit and explore the quality of the parameter maps. MATERIALS AND METHODS 3D MRF was implemented on a custom-built Halbach array using a slab-selective spoiled steady-state free precession sequence with 3D Cartesian readout. Undersampled scans were acquired with different MRF flip angle patterns and reconstructed using matrix completion and matched to the simulated dictionary, taking excitation profile and coil ringing into account. MRF relaxation times were compared to that of inversion recovery (IR) and multi-echo spin echo (MESE) experiments in phantom and in vivo. Furthermore, B0 inhomogeneities were encoded in the MRF sequence using an alternating TE pattern, and the estimated map was used to correct for image distortions in the MRF images using a model-based reconstruction. RESULTS Phantom relaxation times measured with an optimized MRF sequence for low field were in better agreement with reference techniques than for a standard MRF sequence. In vivo muscle relaxation times measured with MRF were longer than those obtained with an IR sequence (T1: 182 ± 21.5 vs 168 ± 9.89 ms) and with an MESE sequence (T2: 69.8 ± 19.7 vs 46.1 ± 9.65 ms). In vivo lipid MRF relaxation times were also longer compared with IR (T1: 165 ± 15.1 ms vs 127 ± 8.28 ms) and with MESE (T2: 160 ± 15.0 ms vs 124 ± 4.27 ms). Integrated ΔB0 estimation and correction resulted in parameter maps with reduced distortions. DISCUSSION It is possible to measure volumetric relaxation times with MRF at 2.5 × 2.5 × 3.0 mm3 resolution in a 13 min scan time on a 50 mT permanent magnet system. The measured MRF relaxation times are longer compared to those measured with reference techniques, especially for T2. This discrepancy can potentially be addressed by hardware, reconstruction and sequence design, but long-term reproducibility needs to be further improved.
Collapse
Affiliation(s)
- Thomas O'Reilly
- Radiology, C.J. Gorter Center for MRI, Leiden University Medical Center, Albinusdreef 2, 2333 ZA, Leiden, The Netherlands
| | - Peter Börnert
- Radiology, C.J. Gorter Center for MRI, Leiden University Medical Center, Albinusdreef 2, 2333 ZA, Leiden, The Netherlands
- Philips Research, Röntgenstraβe 24-26, 22335, Hamburg, Germany
| | - Hongyan Liu
- Computational Imaging Group for MR Diagnostics & Therapy, Center for Imaging Sciences, University Medical Center Utrecht, Heidelberglaan 100, 3584 CX, Utrecht, The Netherlands
| | - Andrew Webb
- Radiology, C.J. Gorter Center for MRI, Leiden University Medical Center, Albinusdreef 2, 2333 ZA, Leiden, The Netherlands
| | - Kirsten Koolstra
- Radiology, Division of Image Processing, Leiden University Medical Center, Albinusdreef 2, 2333 ZA, Leiden, The Netherlands.
| |
Collapse
|
49
|
Guenthner C, Peereboom SM, Dillinger H, McGrath C, Albannay MM, Vishnevskiy V, Fuetterer M, Luechinger R, Jenneskens T, Sturzenegger U, Overweg J, Koken P, Börnert P, Kozerke S. Ramping down a clinical 3 T scanner: a journey into MRI and MRS at 0.75 T. MAGMA (NEW YORK, N.Y.) 2023:10.1007/s10334-023-01089-9. [PMID: 37171689 PMCID: PMC10386956 DOI: 10.1007/s10334-023-01089-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/30/2022] [Revised: 03/06/2023] [Accepted: 04/01/2023] [Indexed: 05/13/2023]
Abstract
OBJECT Lower-field MR is reemerging as a viable, potentially cost-effective alternative to high-field MR, thanks to advances in hardware, sequence design, and reconstruction over the past decades. Evaluation of lower field strengths, however, is limited by the availability of lower-field systems on the market and their considerable procurement costs. In this work, we demonstrate a low-cost, temporary alternative to purchasing a dedicated lower-field MR system. MATERIALS AND METHODS By ramping down an existing clinical 3 T MRI system to 0.75 T, proton signals can be acquired using repurposed 13C transmit/receive hardware and the multi-nuclei spectrometer interface. We describe the ramp-down procedure and necessary software and hardware changes to the system. RESULTS Apart from presenting system characterization results, we show in vivo examples of cardiac cine imaging, abdominal two- and three-point Dixon-type water/fat separation, water/fat-separated MR Fingerprinting, and point-resolved spectroscopy. In addition, the ramp-down approach allows unique comparisons of, e.g., gradient fidelity of the same MR system operated at different field strengths using the same receive chain, gradient coils, and amplifiers. DISCUSSION Ramping down an existing MR system may be seen as a viable alternative for lower-field MR research in groups that already own multi-nuclei hardware and can also serve as a testing platform for custom-made multi-nuclei transmit/receive coils.
Collapse
Affiliation(s)
- Christian Guenthner
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland.
| | | | - Hannes Dillinger
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Charles McGrath
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | | | - Valery Vishnevskiy
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Max Fuetterer
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Roger Luechinger
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | | | | | | | | | | | - Sebastian Kozerke
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| |
Collapse
|
50
|
Soustelle L, Troalen T, Hertanu A, Ranjeva JP, Guye M, Varma G, Alsop DC, Duhamel G, Girard OM. Quantitative magnetization transfer MRI unbiased by on-resonance saturation and dipolar order contributions. Magn Reson Med 2023. [PMID: 37154400 DOI: 10.1002/mrm.29678] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Revised: 03/22/2023] [Accepted: 04/01/2023] [Indexed: 05/10/2023]
Abstract
PURPOSE To demonstrate the bias in quantitative MT (qMT) measures introduced by the presence of dipolar order and on-resonance saturation (ONRS) effects using magnetization transfer (MT) spoiled gradient-recalled (SPGR) acquisitions, and propose changes to the acquisition and analysis strategies to remove these biases. METHODS The proposed framework consists of SPGR sequences prepared with simultaneous dual-offset frequency-saturation pulses to cancel out dipolar order and associated relaxation (T1D ) effects in Z-spectrum acquisitions, and a matched quantitative MT (qMT) mathematical model that includes ONRS effects of readout pulses. Variable flip angle and MT data were fitted jointly to simultaneously estimate qMT parameters (macromolecular proton fraction [MPF], T2,f , T2,b , R, and free pool T1 ). This framework is compared with standard qMT and investigated in terms of reproducibility, and then further developed to follow a joint single-point qMT methodology for combined estimation of MPF and T1 . RESULTS Bland-Altman analyses demonstrated a systematic underestimation of MPF (-2.5% and -1.3%, on average, in white and gray matter, respectively) and overestimation of T1 (47.1 ms and 38.6 ms, on average, in white and gray matter, respectively) if both ONRS and dipolar order effects are ignored. Reproducibility of the proposed framework is excellent (ΔMPF = -0.03% and ΔT1 = -19.0 ms). The single-point methodology yielded consistent MPF and T1 values with respective maximum relative average bias of -0.15% and -3.5 ms found in white matter. CONCLUSION The influence of acquisition strategy and matched mathematical model with regard to ONRS and dipolar order effects in qMT-SPGR frameworks has been investigated. The proposed framework holds promise for improved accuracy with reproducibility.
Collapse
Affiliation(s)
- Lucas Soustelle
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France
- APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| | | | - Andreea Hertanu
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France
- APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| | - Jean-Philippe Ranjeva
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France
- APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| | - Maxime Guye
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France
- APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| | - Gopal Varma
- Division of MR Research, Department of Radiology, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, Massachusetts, USA
| | - David C Alsop
- Division of MR Research, Department of Radiology, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, Massachusetts, USA
| | - Guillaume Duhamel
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France
- APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| | - Olivier M Girard
- Aix Marseille Univ, CNRS, CRMBM, Marseille, France
- APHM, Hôpital Universitaire Timone, CEMEREM, Marseille, France
| |
Collapse
|