1
|
Weine J, McGrath C, Dirix P, Buoso S, Kozerke S. CMRsim-A python package for cardiovascular MR simulations incorporating complex motion and flow. Magn Reson Med 2024; 91:2621-2637. [PMID: 38234037 DOI: 10.1002/mrm.30010] [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: 12/15/2023] [Accepted: 12/22/2023] [Indexed: 01/19/2024]
Abstract
PURPOSE To present an open-source MR simulation framework that facilitates the incorporation of complex motion and flow for studying cardiovascular MR (CMR) acquisition and reconstruction. METHODS CMRsim is a Python package that allows simulation of CMR images using dynamic digital phantoms with complex motion as input. Two simulation paradigms are available, namely, numerical and analytical solutions to the Bloch equations, using a common motion representation. Competitive simulation speeds are achieved using TensorFlow for GPU acceleration. To demonstrate the capability of the package, one introductory and two advanced CMR simulation experiments are presented. The latter showcase phase-contrast imaging of turbulent flow downstream of a stenotic section and cardiac diffusion tensor imaging on a contracting left ventricle. Additionally, extensive documentation and example resources are provided. RESULTS The Bloch simulation with turbulent flow using approximately 1.5 million particles and a sequence duration of 710 ms for each of the seven different velocity encodings took a total of 29 min on a NVIDIA Titan RTX GPU. The results show characteristic phase contrast and magnitude modulation present in real data. The analytical simulation of cardiac diffusion tensor imaging with bulk-motion phase sensitivity took approximately 10 s per diffusion-weighted image, including preparation and loading steps. The results exhibit the expected alteration of diffusion metrics due to strain. CONCLUSION CMRsim is the first simulation framework that allows one to feasibly incorporate complex motion, including turbulent flow, to systematically study advanced CMR acquisition and reconstruction approaches. The open-source package features modularity and transparency, facilitating maintainability and extensibility in support of reproducible research.
Collapse
Affiliation(s)
- Jonathan Weine
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Charles McGrath
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Pietro Dirix
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Stefano Buoso
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Sebastian Kozerke
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| |
Collapse
|
2
|
Alemany I, Rose JN, Ferreira PF, Pennell DJ, Nielles‐Vallespin S, Scott AD, Doorly DJ. Realistic numerical simulations of diffusion tensor cardiovascular magnetic resonance: The effects of perfusion and membrane permeability. Magn Reson Med 2023; 90:1641-1656. [PMID: 37415339 PMCID: PMC10952789 DOI: 10.1002/mrm.29737] [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: 01/13/2023] [Revised: 04/19/2023] [Accepted: 05/16/2023] [Indexed: 07/08/2023]
Abstract
PURPOSE To study the sensitivity of diffusion tensor cardiovascular magnetic resonance (DT-CMR) to microvascular perfusion and changes in cell permeability. METHODS Monte Carlo (MC) random walk simulations in the myocardium have been performed to simulate self-diffusion of water molecules in histology-based media with varying extracellular volume fraction (ECV) and permeable membranes. The effect of microvascular perfusion on simulations of the DT-CMR signal has been incorporated by adding the contribution of particles traveling through an anisotropic capillary network to the diffusion signal. The simulations have been performed considering three pulse sequences with clinical gradient strengths: monopolar stimulated echo acquisition mode (STEAM), monopolar pulsed-gradient spin echo (PGSE), and second-order motion-compensated spin echo (MCSE). RESULTS Reducing ECV intensifies the diffusion restriction and incorporating membrane permeability reduces the anisotropy of the diffusion tensor. Widening the intercapillary velocity distribution results in increased measured diffusion along the cardiomyocytes long axis when the capillary networks are anisotropic. Perfusion amplifies the mean diffusivity for STEAM while the opposite is observed for short diffusion encoding time sequences (PGSE and MCSE). CONCLUSION The effect of perfusion on the measured diffusion tensor is reduced using an increased reference b-value. Our results pave the way for characterization of the response of DT-CMR to microstructural changes underlying cardiac pathology and highlight the higher sensitivity of STEAM to permeability and microvascular circulation due to its longer diffusion encoding time.
Collapse
Affiliation(s)
- Ignasi Alemany
- Department of AeronauticsImperial College LondonLondonUK
- Cardiovascular Magnetic Resonance Unit, Royal Brompton HospitalGuy's and St Thomas' NHS Foundation TrustLondonUK
| | - Jan N. Rose
- Department of AeronauticsImperial College LondonLondonUK
| | - Pedro F. Ferreira
- Cardiovascular Magnetic Resonance Unit, Royal Brompton HospitalGuy's and St Thomas' NHS Foundation TrustLondonUK
- National Heart and Lung InstituteImperial College LondonLondonUK
| | - Dudley J. Pennell
- Cardiovascular Magnetic Resonance Unit, Royal Brompton HospitalGuy's and St Thomas' NHS Foundation TrustLondonUK
- National Heart and Lung InstituteImperial College LondonLondonUK
| | - Sonia Nielles‐Vallespin
- Cardiovascular Magnetic Resonance Unit, Royal Brompton HospitalGuy's and St Thomas' NHS Foundation TrustLondonUK
- National Heart and Lung InstituteImperial College LondonLondonUK
| | - Andrew D. Scott
- Cardiovascular Magnetic Resonance Unit, Royal Brompton HospitalGuy's and St Thomas' NHS Foundation TrustLondonUK
- National Heart and Lung InstituteImperial College LondonLondonUK
| | | |
Collapse
|
3
|
Jing Y, Magnin IE, Frindel C. Monte Carlo simulation of water diffusion through cardiac tissue models. Med Eng Phys 2023; 120:104013. [PMID: 37673779 DOI: 10.1016/j.medengphy.2023.104013] [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: 10/22/2022] [Revised: 05/13/2023] [Accepted: 06/22/2023] [Indexed: 09/08/2023]
Abstract
Monte Carlo diffusion simulations are commonly used to establish a reliable ground truth of tissue microstructure, including for the validation of diffusion-weighted MRI. However, selecting simulation parameters is challenging and affects validity and reproducibility. We conducted experiments to investigate critical conditions in Monte Carlo simulations, such as tissue representation complexity, simulated molecules, update duration, and compartment size. Results show significant changes in microstructure characteristics when parameters are altered, emphasizing the importance of careful control for a reliable ground truth.
Collapse
Affiliation(s)
- Yuhan Jing
- Univ Lyon, INSA-Lyon, Université Claude Bernard Lyon 1, UJM-Saint Etienne, CNRS, Inserm, CREATIS UMR 5220, U1294, 21 Avenue Jean Capelle, Lyon, 69621, France
| | - Isabelle E Magnin
- Univ Lyon, INSA-Lyon, Université Claude Bernard Lyon 1, UJM-Saint Etienne, CNRS, Inserm, CREATIS UMR 5220, U1294, 21 Avenue Jean Capelle, Lyon, 69621, France
| | - Carole Frindel
- Univ Lyon, INSA-Lyon, Université Claude Bernard Lyon 1, UJM-Saint Etienne, CNRS, Inserm, CREATIS UMR 5220, U1294, 21 Avenue Jean Capelle, Lyon, 69621, France.
| |
Collapse
|
4
|
Lashgari M, Ravikumar N, Teh I, Li JR, Buckley DL, Schneider JE, Frangi AF. Three-dimensional micro-structurally informed in silico myocardium-Towards virtual imaging trials in cardiac diffusion weighted MRI. Med Image Anal 2022; 82:102592. [PMID: 36095906 DOI: 10.1016/j.media.2022.102592] [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/27/2021] [Revised: 08/14/2022] [Accepted: 08/18/2022] [Indexed: 11/24/2022]
Abstract
In silico tissue models (viz. numerical phantoms) provide a mechanism for evaluating quantitative models of magnetic resonance imaging. This includes the validation and sensitivity analysis of imaging biomarkers and tissue microstructure parameters. This study proposes a novel method to generate a realistic numerical phantom of myocardial microstructure. The proposed method extends previous studies by accounting for the variability of the cardiomyocyte shape, water exchange between the cardiomyocytes (intercalated discs), disorder class of myocardial microstructure, and four sheetlet orientations. In the first stage of the method, cardiomyocytes and sheetlets are generated by considering the shape variability and intercalated discs in cardiomyocyte-cardiomyocyte connections. Sheetlets are then aggregated and oriented in the directions of interest. The morphometric study demonstrates no significant difference (p>0.01) between the distribution of volume, length, and primary and secondary axes of the numerical and real (literature) cardiomyocyte data. Moreover, structural correlation analysis validates that the in-silico tissue is in the same class of disorderliness as the real tissue. Additionally, the absolute angle differences between the simulated helical angle (HA) and input HA (reference value) of the cardiomyocytes (4.3°±3.1°) demonstrate a good agreement with the absolute angle difference between the measured HA using experimental cardiac diffusion tensor imaging (cDTI) and histology (reference value) reported by (Holmes et al., 2000) (3.7°±6.4°) and (Scollan et al. 1998) (4.9°±14.6°). Furthermore, the angular distance between eigenvectors and sheetlet angles of the input and simulated cDTI is much smaller than those between measured angles using structural tensor imaging (as a gold standard) and experimental cDTI. Combined with the qualitative results, these results confirm that the proposed method can generate richer numerical phantoms for the myocardium than previous studies.
Collapse
Affiliation(s)
- Mojtaba Lashgari
- Centre for Computational Imaging and Simulation Technologies in Biomedicine (CISTIB), School of Computing, University of Leeds, Leeds, UK; Biomedical Imaging Science Department, Leeds Institute for Cardiovascular and Metabolic Medicine (LICAMM), School of Medicine, University of Leeds, Leeds, UK.
| | - Nishant Ravikumar
- Centre for Computational Imaging and Simulation Technologies in Biomedicine (CISTIB), School of Computing, University of Leeds, Leeds, UK; Biomedical Imaging Science Department, Leeds Institute for Cardiovascular and Metabolic Medicine (LICAMM), School of Medicine, University of Leeds, Leeds, UK
| | - Irvin Teh
- Biomedical Imaging Science Department, Leeds Institute for Cardiovascular and Metabolic Medicine (LICAMM), School of Medicine, University of Leeds, Leeds, UK
| | - Jing-Rebecca Li
- INRIA Saclay, Equipe DEFI, CMAP, Ecole Polytechnique, Route de Saclay, 91128 Palaiseau Cedex, France
| | - David L Buckley
- Biomedical Imaging Science Department, Leeds Institute for Cardiovascular and Metabolic Medicine (LICAMM), School of Medicine, University of Leeds, Leeds, UK
| | - Jurgen E Schneider
- Biomedical Imaging Science Department, Leeds Institute for Cardiovascular and Metabolic Medicine (LICAMM), School of Medicine, University of Leeds, Leeds, UK
| | - Alejandro F Frangi
- Centre for Computational Imaging and Simulation Technologies in Biomedicine (CISTIB), School of Computing, University of Leeds, Leeds, UK; Biomedical Imaging Science Department, Leeds Institute for Cardiovascular and Metabolic Medicine (LICAMM), School of Medicine, University of Leeds, Leeds, UK; INRIA Saclay, Equipe DEFI, CMAP, Ecole Polytechnique, Route de Saclay, 91128 Palaiseau Cedex, France; Medical Imaging Research Center (MIRC), Department of Cardiovascular Sciences, KU Leuven, Leuven, Belgium; Medical Imaging Research Center (MIRC), Department of Electrical Engineering, KU Leuven, Leuven, Belgium; Alan Turing Institute, London, UK.
| |
Collapse
|
5
|
Scott AD, Jackson T, Khalique Z, Gorodezky M, Pardoe B, Begum L, Bruno VD, Chowdhury RA, Ferreira PF, Nielles‐Vallespin S, Roehl M, McCarthy KP, Sarathchandra P, Rose JN, Doorly DJ, Pennell DJ, Ascione R, de Silva R, Firmin DN. Development of a cardiovascular magnetic resonance-compatible large animal isolated heart model for direct comparison of beating and arrested hearts. NMR IN BIOMEDICINE 2022; 35:e4692. [PMID: 35040195 PMCID: PMC9286060 DOI: 10.1002/nbm.4692] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/21/2021] [Revised: 12/21/2021] [Accepted: 01/07/2022] [Indexed: 06/02/2023]
Abstract
Cardiac motion results in image artefacts and quantification errors in many cardiovascular magnetic resonance (CMR) techniques, including microstructural assessment using diffusion tensor cardiovascular magnetic resonance (DT-CMR). Here, we develop a CMR-compatible isolated perfused porcine heart model that allows comparison of data obtained in beating and arrested states. Ten porcine hearts (8/10 for protocol optimisation) were harvested using a donor heart retrieval protocol and transported to the remote CMR facility. Langendorff perfusion in a 3D-printed chamber and perfusion circuit re-established contraction. Hearts were imaged using cine, parametric mapping and STEAM DT-CMR at cardiac phases with the minimum and maximum wall thickness. High potassium and lithium perfusates were then used to arrest the heart in a slack and contracted state, respectively. Imaging was repeated in both arrested states. After imaging, tissue was removed for subsequent histology in a location matched to the DT-CMR data using fiducial markers. Regular sustained contraction was successfully established in six out of 10 hearts, including the final five hearts. Imaging was performed in four hearts and one underwent the full protocol, including colocalised histology. The image quality was good and there was good agreement between DT-CMR data in equivalent beating and arrested states. Despite the use of autologous blood and dextran within the perfusate, T2 mapping results, DT-CMR measures and an increase in mass were consistent with development of myocardial oedema, resulting in failure to achieve a true diastolic-like state. A contiguous stack of 313 5-μm histological sections at and a 100-μm thick section showing cell morphology on 3D fluorescent confocal microscopy colocalised to DT-CMR data were obtained. A CMR-compatible isolated perfused beating heart setup for large animal hearts allows direct comparisons of beating and arrested heart data with subsequent colocalised histology, without the need for onsite preclinical facilities.
Collapse
Affiliation(s)
- Andrew D. Scott
- Cardiovascular Magnetic Resonance UnitRoyal Brompton HospitalLondonUK
- National Heart and Lung InstituteImperial CollegeLondonUK
| | - Tim Jackson
- Department of PerfusionRoyal Brompton HospitalLondonUK
| | - Zohya Khalique
- Cardiovascular Magnetic Resonance UnitRoyal Brompton HospitalLondonUK
- National Heart and Lung InstituteImperial CollegeLondonUK
| | - Margarita Gorodezky
- Cardiovascular Magnetic Resonance UnitRoyal Brompton HospitalLondonUK
- National Heart and Lung InstituteImperial CollegeLondonUK
| | - Ben Pardoe
- Department of PerfusionRoyal Brompton HospitalLondonUK
| | - Lale Begum
- Department of PerfusionRoyal Brompton HospitalLondonUK
| | - V. Domenico Bruno
- Translational Biomedical Research CentreUniversity of BristolBristolUK
- Bristol Heart InstituteUniversity Hospital Bristol NHS Foundation TrustBristolUK
| | - Rasheda A. Chowdhury
- National Heart and Lung InstituteImperial CollegeLondonUK
- Imperial Centre for Cardiac EngineeringImperial CollegeLondonUK
| | - Pedro F. Ferreira
- Cardiovascular Magnetic Resonance UnitRoyal Brompton HospitalLondonUK
- National Heart and Lung InstituteImperial CollegeLondonUK
| | - Sonia Nielles‐Vallespin
- Cardiovascular Magnetic Resonance UnitRoyal Brompton HospitalLondonUK
- National Heart and Lung InstituteImperial CollegeLondonUK
| | - Malte Roehl
- Cardiovascular Magnetic Resonance UnitRoyal Brompton HospitalLondonUK
- National Heart and Lung InstituteImperial CollegeLondonUK
| | | | - Padmini Sarathchandra
- National Heart and Lung InstituteImperial CollegeLondonUK
- Magdi Yacoub Institute, National Heart and Lung InstituteImperial CollegeLondonUK
| | - Jan N. Rose
- Department of AeronauticsImperial CollegeLondonUK
| | | | - Dudley J. Pennell
- Cardiovascular Magnetic Resonance UnitRoyal Brompton HospitalLondonUK
- National Heart and Lung InstituteImperial CollegeLondonUK
| | - Raimondo Ascione
- Translational Biomedical Research CentreUniversity of BristolBristolUK
- Bristol Heart InstituteUniversity Hospital Bristol NHS Foundation TrustBristolUK
| | - Ranil de Silva
- Cardiovascular Magnetic Resonance UnitRoyal Brompton HospitalLondonUK
- National Heart and Lung InstituteImperial CollegeLondonUK
| | - David N. Firmin
- Cardiovascular Magnetic Resonance UnitRoyal Brompton HospitalLondonUK
- National Heart and Lung InstituteImperial CollegeLondonUK
| |
Collapse
|
6
|
Random walk diffusion simulations in semi-permeable layered media with varying diffusivity. Sci Rep 2022; 12:10759. [PMID: 35750717 PMCID: PMC9232609 DOI: 10.1038/s41598-022-14541-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Accepted: 06/08/2022] [Indexed: 11/09/2022] Open
Abstract
In this paper we present random walk based solutions to diffusion in semi-permeable layered media with varying diffusivity. We propose a novel transit model for solving the interaction of random walkers with a membrane. This hybrid model is based on treating the membrane permeability and the step change in diffusion coefficient as two interactions separated by an infinitesimally small layer. By conducting an extensive analytical flux analysis, the performance of our hybrid model is compared with a commonly used membrane transit model (reference model). Numerical simulations demonstrate the limitations of the reference model in dealing with step changes in diffusivity and show the capability of the hybrid model to overcome this limitation and to offer substantial gains in computational efficiency. The suitability of both random walk transit models for the application to simulations of diffusion tensor cardiovascular magnetic resonance (DT-CMR) imaging is assessed in a histology-based domain relevant to DT-CMR. In order to demonstrate the usefulness of the new hybrid model for other possible applications, we also consider a larger range of permeabilities beyond those commonly found in biological tissues.
Collapse
|
7
|
Weine J, van Gorkum RJH, Stoeck CT, Vishnevskiy V, Kozerke S. Synthetically Trained Convolutional Neural Networks for Improved Tensor Estimation from Free-Breathing Cardiac DTI. Comput Med Imaging Graph 2022; 99:102075. [DOI: 10.1016/j.compmedimag.2022.102075] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2021] [Revised: 03/15/2022] [Accepted: 05/05/2022] [Indexed: 10/18/2022]
|
8
|
Rahman T, Moulin K, Perotti LE. Cardiac Diffusion Tensor Biomarkers of Chronic Infarction Based on In Vivo Data. APPLIED SCIENCES-BASEL 2022; 12. [PMID: 36032414 PMCID: PMC9408809 DOI: 10.3390/app12073512] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/10/2022]
Abstract
In vivo cardiac diffusion tensor imaging (cDTI) data were acquired in
swine subjects six to ten weeks post-myocardial infarction (MI) to identify
microstructural-based biomarkers of MI. Diffusion tensor invariants, diffusion
tensor eigenvalues, and radial diffusivity (RD) are evaluated in the infarct,
border, and remote myocardium, and compared with extracellular volume fraction
(ECV) and native T1 values. Additionally, to aid the interpretation of the
experimental results, the diffusion of water molecules was numerically simulated
as a function of ECV. Finally, findings based on in vivo measures were confirmed
using higher-resolution and higher signal-to-noise data acquired ex vivo in the
same subjects. Mean diffusivity, diffusion tensor eigenvalues, and RD increased
in the infarct and border regions compared to remote myocardium, while
fractional anisotropy decreased. Secondary (e2) and tertiary
(e3) eigenvalues increased more significantly than the primary
eigenvalue in the infarct and border regions. These findings were confirmed by
the diffusion simulations. Although ECV presented the largest increase in
infarct and border regions, e2, e3, and RD increased the
most among non-contrast-based biomarkers. RD is of special interest as it
summarizes the changes occurring in the radial direction and may be more robust
than e2 or e3 alone.
Collapse
Affiliation(s)
- Tanjib Rahman
- Department of Mechanical and Aerospace Engineering,
University of Central Florida, Orlando, FL 32816, USA
| | - Kévin Moulin
- CREATIS Laboratory, Univ. Lyon, UJM-Saint-Etienne, INSA,
CNRS UMR 5520, INSERM, 69100 Villeurbanne, France
- Department of Radiology, University Hospital Saint-Etienne,
42270 Saint-Priest-en-Jarez, France
| | - Luigi E. Perotti
- Department of Mechanical and Aerospace Engineering,
University of Central Florida, Orlando, FL 32816, USA
- Correspondence:
| |
Collapse
|
9
|
Connecting macroscopic diffusion metrics of cardiac diffusion tensor imaging and microscopic myocardial structures based on simulation. Med Image Anal 2022; 77:102325. [DOI: 10.1016/j.media.2021.102325] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2019] [Revised: 11/26/2021] [Accepted: 11/29/2021] [Indexed: 11/20/2022]
|
10
|
Dillinger H, McGrath C, Guenthner C, Kozerke S. Fundamentals of turbulent flow spectrum imaging. Magn Reson Med 2021; 87:1231-1249. [PMID: 34786764 PMCID: PMC9299145 DOI: 10.1002/mrm.29001] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2021] [Revised: 08/12/2021] [Accepted: 08/18/2021] [Indexed: 12/13/2022]
Abstract
PURPOSE To introduce a mathematical framework and in-silico validation of turbulent flow spectrum imaging (TFSI) of stenotic flow using phase-contrast MRI, evaluate systematic errors in quantitative turbulence parameter estimation, and propose a novel method for probing the Lagrangian velocity spectra of turbulent flows. THEORY AND METHODS The spectral response of velocity-encoding gradients is derived theoretically and linked to turbulence parameter estimation including the velocity autocorrelation function spectrum. Using a phase-contrast MRI simulation framework, the encoding properties of bipolar gradient waveforms with identical first gradient moments but different duration are investigated on turbulent flow data of defined characteristics as derived from computational fluid dynamics. Based on theoretical insights, an approach using velocity-compensated gradient waveforms is proposed to specifically probe desired ranges of the velocity autocorrelation function spectrum with increased accuracy. RESULTS Practical velocity-encoding gradients exhibit limited encoding power of typical turbulent flow spectra, resulting in up to 50% systematic underestimation of intravoxel SD values. Depending on the turbulence level in fluids, the error due to a single encoding gradient spectral response can vary by 20%. When using tailored velocity-compensated gradients, improved quantification of the Lagrangian velocity spectrum on a voxel-by-voxel basis is achieved and used for quantitative correction of intravoxel SD values estimated with velocity-encoding gradients. CONCLUSION To address systematic underestimation of turbulence parameters using bipolar velocity-encoding gradients in phase-contrast MRI of stenotic flows with short correlation times, tailored velocity-compensated gradients are proposed to improve quantitative mapping of turbulent blood flow characteristics.
Collapse
Affiliation(s)
- Hannes Dillinger
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Charles McGrath
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Christian Guenthner
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Sebastian Kozerke
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| |
Collapse
|
11
|
Abstract
Advances in technology have made it possible to image the microstructure of the heart with diffusion-weighted magnetic resonance. The technique provides unique insights into the cellular architecture of the myocardium and how this is perturbed in a range of disease contexts. In this review, the physical basis of diffusion MRI and the challenges of implementing it in the beating heart are discussed. Cutting edge acquisition and analysis techniques, as well as the results of initial clinical studies, are reported.
Collapse
Affiliation(s)
- David E Sosnovik
- Cardiology Division, Cardiovascular Research Center, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA; Department of Radiology, Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA; Division of Health Sciences and Technology, Harvard Medical School and Massachusetts Institute of Technology, Cambridge, MA, USA.
| |
Collapse
|
12
|
Naughton NM, Tennyson CG, Georgiadis JG. Lattice Boltzmann method for simulation of diffusion magnetic resonance imaging physics in multiphase tissue models. Phys Rev E 2020; 102:043305. [PMID: 33212689 DOI: 10.1103/physreve.102.043305] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2019] [Accepted: 08/31/2020] [Indexed: 06/11/2023]
Abstract
We report an implementation of the lattice Boltzmann method (LBM) to integrate the Bloch-Torrey equation, which describes the evolution of the transverse magnetization vector and the fate of the signal of diffusion magnetic resonance imaging (dMRI). Motivated by the need to interpret dMRI experiments in biological tissues, and to offset the small time-step limitation of classical LBM, a hybrid LBM scheme is introduced and implemented to solve the Bloch-Torrey equation. A membrane boundary condition is presented which is able to accurately represent the effects of thin curvilinear membranes typically found in biological tissues. As implemented, the hybrid LBM scheme accommodates piece-wise uniform transport, dMRI parameters, periodic and mirroring outer boundary conditions, and finite membrane permeabilities on non-boundary-conforming inner boundaries. By comparing with analytical solutions of limiting cases, we demonstrate that the hybrid LBM scheme is more accurate than the classical LBM scheme. The proposed explicit LBM scheme maintains second-order spatial accuracy, stability, and first-order temporal accuracy for a wide range of parameters. The parallel implementation of the hybrid LBM code in a multi-CPU computer system, as well as on GPUs, is straightforward and efficient. Along with offering certain advantages over finite element or Monte Carlo schemes, the proposed hybrid LBM constitutes a flexible scheme that can by easily adapted to model more complex interfacial conditions and physics in heterogeneous multiphase tissue models and to accommodate sophisticated dMRI sequences.
Collapse
Affiliation(s)
- Noel M Naughton
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | | | - John G Georgiadis
- Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
- Department of Biomedical Engineering, Illinois Institute of Technology, Chicago, Illinois 60616, USA
| |
Collapse
|
13
|
Stoeck CT, von Deuster C, van Gorkum RJH, Kozerke S. Motion and eddy current-induced signal dephasing in in vivo cardiac DTI. Magn Reson Med 2019; 84:277-288. [PMID: 31868257 DOI: 10.1002/mrm.28132] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2019] [Revised: 11/11/2019] [Accepted: 11/25/2019] [Indexed: 12/14/2022]
Abstract
PURPOSE To address motion in cardiac DWI, stimulated-echo acquisition mode (STEAM) and second-order motion-compensated spin-echo (SE) sequences have been proposed. Despite applying motion-compensation strategies, residual motion can cause misleading signal attenuation. The purpose of this study is to estimate the motion-induced error in both sequences by analysis of image phase. METHODS Diffusion-weighted motion-compensated SE sequences and STEAM imaging was applied in vivo with diffusion encoding along 3 orthogonal directions. A b-value range of 100 to 600 s/mm2 and trigger delays of 25%, 50%, and 75% of end systole and middiastole were used. Eddy-current contributions were obtained from phantom measurements. After computation of motion-induced phase maps, the amount of signal dephasing was computed from phase gradients, and the resulting errors in diffusion tensor parameters were calculated. RESULTS Motion-induced dephasing from the STEAM sequence showed less dependency on the b-value and no dependency on the heart phase, whereas SE imaging performed best at 75% end systole followed by 50% end systole and middiastole. For a typical experimental setting, errors of 3.3%/3.0% mean diffusivity, 4.9%/4.8% fractional anisotropy, 2.9º/3.2º helix angulation, 0.8º/0.7º transverse angulation, and 9.9º/10.0º sheet angulation (SE/STEAM) were calculated. CONCLUSION Image phase contains valuable information regarding uncompensated motion and eddy currents in cardiac DTI. Although the trigger delay window for SE is narrower compared with the STEAM-based approach, imaging in both systole and diastole is feasible and both sequences perform similarly if the trigger delays are selected carefully with SE.
Collapse
Affiliation(s)
- Christian T Stoeck
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | | | | | - Sebastian Kozerke
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| |
Collapse
|
14
|
Diffusion tensor cardiovascular magnetic resonance in hypertrophic cardiomyopathy: a comparison of motion-compensated spin echo and stimulated echo techniques. MAGNETIC RESONANCE MATERIALS IN PHYSICS BIOLOGY AND MEDICINE 2019; 33:331-342. [PMID: 31758419 PMCID: PMC7230046 DOI: 10.1007/s10334-019-00799-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/12/2019] [Revised: 10/15/2019] [Accepted: 11/07/2019] [Indexed: 11/18/2022]
Abstract
Objectives Diffusion tensor cardiovascular magnetic resonance (DT-CMR) interrogates myocardial microstructure. Two frequently used in vivo DT-CMR techniques are motion-compensated spin echo (M2-SE) and stimulated echo acquisition mode (STEAM). Whilst M2-SE is strain-insensitive and signal to noise ratio efficient, STEAM has a longer diffusion time and motion compensation is unnecessary. Here we compare STEAM and M2-SE DT-CMR in patients. Materials and methods Biphasic DT-CMR using STEAM and M2-SE, late gadolinium imaging and pre/post gadolinium T1-mapping were performed in a mid-ventricular short-axis slice, in ten hypertrophic cardiomyopathy (HCM) patients at 3 T. Results Adequate quality data were obtained from all STEAM, but only 7/10 (systole) and 4/10 (diastole) M2-SE acquisitions. Compared with STEAM, M2-SE yielded higher systolic mean diffusivity (MD) (p = 0.02) and lower fractional anisotropy (FA) (p = 0.02, systole). Compared with segments with neither hypertrophy nor late gadolinium, segments with both had lower systolic FA using M2-SE (p = 0.02) and trend toward higher MD (p = 0.1). The negative correlation between FA and extracellular volume fraction was stronger with STEAM than M2-SE (r2 = 0.29, p < 0.001 STEAM vs. r2 = 0.10, p = 0.003 M2-SE). Discussion In HCM, only STEAM reliably assesses biphasic myocardial microstructure. Higher MD and lower FA from M2-SE reflect the shorter diffusion times. Further work will relate DT-CMR parameters and microstructural changes in disease. Electronic supplementary material The online version of this article (10.1007/s10334-019-00799-3) contains supplementary material, which is available to authorized users.
Collapse
|
15
|
Naughton NM, Georgiadis JG. Global sensitivity analysis of skeletal muscle dMRI metrics: Effects of microstructural and pulse parameters. Magn Reson Med 2019; 83:1458-1470. [DOI: 10.1002/mrm.28014] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2019] [Revised: 08/26/2019] [Accepted: 09/05/2019] [Indexed: 12/13/2022]
Affiliation(s)
- Noel M. Naughton
- Department of Mechanical Science and Engineering University of Illinois at Urbana‐Champaign Urbana Illinois
| | - John G. Georgiadis
- Department of Mechanical Science and Engineering University of Illinois at Urbana‐Champaign Urbana Illinois
- Department of Biomedical Engineering Illinois Institute of Technology Chicago Illinois
| |
Collapse
|
16
|
Spinner GR, Stoeck CT, Mathez L, von Deuster C, Federau C, Kozerke S. On probing intravoxel incoherent motion in the heart‐spin‐echo versus stimulated‐echo DWI. Magn Reson Med 2019; 82:1150-1163. [DOI: 10.1002/mrm.27777] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2018] [Revised: 03/06/2019] [Accepted: 03/27/2019] [Indexed: 02/06/2023]
Affiliation(s)
- Georg R. Spinner
- Institute for Biomedical Engineering University and ETH Zurich Zurich Switzerland
| | - Christian T. Stoeck
- Institute for Biomedical Engineering University and ETH Zurich Zurich Switzerland
| | - Linda Mathez
- Institute for Biomedical Engineering University and ETH Zurich Zurich Switzerland
| | | | - Christian Federau
- Institute for Biomedical Engineering University and ETH Zurich Zurich Switzerland
| | - Sebastian Kozerke
- Institute for Biomedical Engineering University and ETH Zurich Zurich Switzerland
| |
Collapse
|
17
|
Rose JN, Nielles-Vallespin S, Ferreira PF, Firmin DN, Scott AD, Doorly DJ. Novel insights into in-vivo diffusion tensor cardiovascular magnetic resonance using computational modeling and a histology-based virtual microstructure. Magn Reson Med 2018; 81:2759-2773. [PMID: 30350880 PMCID: PMC6637383 DOI: 10.1002/mrm.27561] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2018] [Revised: 09/12/2018] [Accepted: 09/14/2018] [Indexed: 12/13/2022]
Abstract
Purpose To develop histology‐informed simulations of diffusion tensor cardiovascular magnetic resonance (DT‐CMR) for typical in‐vivo pulse sequences and determine their sensitivity to changes in extra‐cellular space (ECS) and other microstructural parameters. Methods We synthesised the DT‐CMR signal from Monte Carlo random walk simulations. The virtual tissue was based on porcine histology. The cells were thickened and then shrunk to modify ECS. We also created idealised geometries using cuboids in regular arrangement, matching the extra‐cellular volume fraction (ECV) of 16–40%. The simulated voxel size was 2.8 × 2.8 × 8.0 mm3 for pulse sequences covering short and long diffusion times: Stejskal–Tanner pulsed‐gradient spin echo, second‐order motion‐compensated spin echo, and stimulated echo acquisition mode (STEAM), with clinically available gradient strengths. Results The primary diffusion tensor eigenvalue increases linearly with ECV at a similar rate for all simulated geometries. Mean diffusivity (MD) varies linearly, too, but is higher for the substrates with more uniformly distributed ECS. Fractional anisotropy (FA) for the histology‐based geometry is higher than the idealised geometry with low sensitivity to ECV, except for the long mixing time of the STEAM sequence. Varying the intra‐cellular diffusivity (DIC) results in large changes of MD and FA. Varying extra‐cellular diffusivity or using stronger gradients has minor effects on FA. Uncertainties of the primary eigenvector orientation are reduced using STEAM. Conclusions We found that the distribution of ECS has a measurable impact on DT‐CMR parameters. The observed sensitivity of MD and FA to ECV and DIC has potentially interesting applications for interpreting in‐vivo DT‐CMR parameters.
Collapse
Affiliation(s)
- Jan N Rose
- Department of Aeronautics, Imperial College London, London, United Kingdom
| | - Sonia Nielles-Vallespin
- Cardiovascular Magnetic Resonance Unit, The Royal Brompton Hospital, London, United Kingdom.,National Heart and Lung Institute, Imperial College London, London, United Kingdom
| | - Pedro F Ferreira
- Cardiovascular Magnetic Resonance Unit, The Royal Brompton Hospital, London, United Kingdom.,National Heart and Lung Institute, Imperial College London, London, United Kingdom
| | - David N Firmin
- Cardiovascular Magnetic Resonance Unit, The Royal Brompton Hospital, London, United Kingdom.,National Heart and Lung Institute, Imperial College London, London, United Kingdom
| | - Andrew D Scott
- Cardiovascular Magnetic Resonance Unit, The Royal Brompton Hospital, London, United Kingdom.,National Heart and Lung Institute, Imperial College London, London, United Kingdom
| | - Denis J Doorly
- Department of Aeronautics, Imperial College London, London, United Kingdom
| |
Collapse
|