1
|
Shalom ES, Van Loo S, Khan A, Sourbron SP. Identifiability of spatiotemporal tissue perfusion models. Phys Med Biol 2024; 69:115034. [PMID: 38636525 DOI: 10.1088/1361-6560/ad4087] [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: 12/06/2023] [Accepted: 04/18/2024] [Indexed: 04/20/2024]
Abstract
Objective.Standard models for perfusion quantification in DCE-MRI produce a bias by treating voxels as isolated systems. Spatiotemporal models can remove this bias, but it is unknown whether they are fundamentally identifiable. The aim of this study is to investigate this question in silico using one-dimensional toy systems with a one-compartment blood flow model and a two-compartment perfusion model.Approach.For each of the two models, identifiability is explored theoretically and in-silico for three systems. Concentrations over space and time are simulated by forward propagation. Different levels of noise and temporal undersampling are added to investigate sensitivity to measurement error. Model parameters are fitted using a standard gradient descent algorithm, applied iteratively with a stepwise increasing time window. Model fitting is repeated with different initial values to probe uniqueness of the solution. Reconstruction accuracy is quantified for each parameter by comparison to the ground truth.Main results.Theoretical analysis shows that flows and volume fractions are only identifiable up to a constant, and that this degeneracy can be removed by proper choice of parameters. Simulations show that in all cases, the tissue concentrations can be reconstructed accurately. The one-compartment model shows accurate reconstruction of blood velocities and arterial input functions, independent of the initial values and robust to measurement error. The two-compartmental perfusion model was not fully identifiable, showing good reconstruction of arterial velocities and input functions, but multiple valid solutions for the perfusion parameters and venous velocities, and a strong sensitivity to measurement error in these parameters.Significance.These results support the use of one-compartment spatiotemporal flow models, but two-compartment perfusion models were not sufficiently identifiable. Future studies should investigate whether this degeneracy is resolved in more realistic 2D and 3D systems, by adding physically justified constraints, or by optimizing experimental parameters such as injection duration or temporal resolution.
Collapse
Affiliation(s)
- Eve S Shalom
- School of Physics and Astronomy, The University of Leeds, United Kingdom
| | - Sven Van Loo
- Department of Applied Physics, Ghent University, Belgium
| | - Amirul Khan
- School of Civil Engineering, The University of Leeds, United Kingdom
| | - Steven P Sourbron
- Division of Clinical Medicine, University of Sheffield, United Kingdom
| |
Collapse
|
2
|
Oh G, Moon Y, Moon WJ, Ye JC. Unpaired deep learning for pharmacokinetic parameter estimation from dynamic contrast-enhanced MRI without AIF measurements. Neuroimage 2024; 291:120571. [PMID: 38518829 DOI: 10.1016/j.neuroimage.2024.120571] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Revised: 02/28/2024] [Accepted: 03/08/2024] [Indexed: 03/24/2024] Open
Abstract
DCE-MRI provides information about vascular permeability and tissue perfusion through the acquisition of pharmacokinetic parameters. However, traditional methods for estimating these pharmacokinetic parameters involve fitting tracer kinetic models, which often suffer from computational complexity and low accuracy due to noisy arterial input function (AIF) measurements. Although some deep learning approaches have been proposed to tackle these challenges, most existing methods rely on supervised learning that requires paired input DCE-MRI and labeled pharmacokinetic parameter maps. This dependency on labeled data introduces significant time and resource constraints and potential noise in the labels, making supervised learning methods often impractical. To address these limitations, we present a novel unpaired deep learning method for estimating pharmacokinetic parameters and the AIF using a physics-driven CycleGAN approach. Our proposed CycleGAN framework is designed based on the underlying physics model, resulting in a simpler architecture with a single generator and discriminator pair. Crucially, our experimental results indicate that our method does not necessitate separate AIF measurements and produces more reliable pharmacokinetic parameters than other techniques.
Collapse
Affiliation(s)
- Gyutaek Oh
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), 291, Daehak-ro, Yuseong-gu, 34141, Daejeon, Republic of Korea
| | - Yeonsil Moon
- Department of Neurology, Konkuk University Medical Center, 120-1, Neungdong-ro, Gwangjin-gu, 05030, Seoul, Republic of Korea
| | - Won-Jin Moon
- Department of Radiology, Konkuk University Medical Center, 120-1, Neungdong-ro, Gwangjin-gu, 05030, Seoul, Republic of Korea.
| | - Jong Chul Ye
- Kim Jaechul Graduate School of AI, Korea Advanced Institute of Science and Technology (KAIST), 291, Daehak-ro, Yuseong-gu, 34141, Daejeon, Republic of Korea.
| |
Collapse
|
3
|
Rotkopf LT, Zhang KS, Tavakoli AA, Bonekamp D, Ziener CH, Schlemmer HP. Quantitative Analysis of DCE and DSC-MRI: From Kinetic Modeling to Deep Learning. ROFO-FORTSCHR RONTG 2022; 194:975-982. [PMID: 35211930 DOI: 10.1055/a-1762-5854] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Abstract
BACKGROUND Perfusion MRI is a well-established imaging modality with a multitude of applications in oncological and cardiovascular imaging. Clinically used processing methods, while stable and robust, have remained largely unchanged in recent years. Despite promising results from novel methods, their relatively minimal improvement compared to established methods did not generally warrant significant changes to clinical perfusion processing. RESULTS AND CONCLUSION Machine learning in general and deep learning in particular, which are currently revolutionizing computer-aided diagnosis, may carry the potential to change this situation and truly capture the potential of perfusion imaging. Recent advances in the training of recurrent neural networks make it possible to predict and classify time series data with high accuracy. Combining physics-based tissue models and deep learning, using either physics-informed neural networks or universal differential equations, simplifies the training process and increases the interpretability of the resulting models. Due to their versatility, these methods will potentially be useful in bridging the gap between microvascular architecture and perfusion parameters, akin to MR fingerprinting in structural MR imaging. Still, further research is urgently needed before these methods may be used in clinical practice. KEY POINTS · Machine learning offers promising methods for processing of perfusion data.. · Recurrent neural networks can classify time series with high accuracy.. · Data augmentation is essentially especially when using small datasets.. CITATION FORMAT · Rotkopf LT, Zhang KS, Tavakoli AA et al. Quantitative Analysis of DCE and DSC-MRI: From Kinetic Modeling to Deep Learning. Fortschr Röntgenstr 2022; DOI: 10.1055/a-1762-5854.
Collapse
Affiliation(s)
- Lukas T Rotkopf
- Department of Radiology, German Cancer Research Centre, Heidelberg, Germany
| | - Kevin Sun Zhang
- Department of Radiology, German Cancer Research Centre, Heidelberg, Germany
| | | | - David Bonekamp
- Department of Radiology, German Cancer Research Centre, Heidelberg, Germany
| | | | | |
Collapse
|
4
|
Flouri D, Lesnic D, Chrysochou C, Parikh J, Thelwall P, Sheerin N, Kalra PA, Buckley DL, Sourbron SP. Motion correction of free-breathing magnetic resonance renography using model-driven registration. MAGMA (NEW YORK, N.Y.) 2021; 34:805-822. [PMID: 34160718 PMCID: PMC8578117 DOI: 10.1007/s10334-021-00936-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/11/2021] [Revised: 05/24/2021] [Accepted: 06/08/2021] [Indexed: 11/29/2022]
Abstract
INTRODUCTION Model-driven registration (MDR) is a general approach to remove patient motion in quantitative imaging. In this study, we investigate whether MDR can effectively correct the motion in free-breathing MR renography (MRR). MATERIALS AND METHODS MDR was generalised to linear tracer-kinetic models and implemented using 2D or 3D free-form deformations (FFD) with multi-resolution and gradient descent optimization. MDR was evaluated using a kidney-mimicking digital reference object (DRO) and free-breathing patient data acquired at high temporal resolution in multi-slice 2D (5 patients) and 3D acquisitions (8 patients). Registration accuracy was assessed using comparison to ground truth DRO, calculating the Hausdorff distance (HD) between ground truth masks with segmentations and visual evaluation of dynamic images, signal-time courses and parametric maps (all data). RESULTS DRO data showed that the bias and precision of parameter maps after MDR are indistinguishable from motion-free data. MDR led to reduction in HD (HDunregistered = 9.98 ± 9.76, HDregistered = 1.63 ± 0.49). Visual inspection showed that MDR effectively removed motion effects in the dynamic data, leading to a clear improvement in anatomical delineation on parametric maps and a reduction in motion-induced oscillations on signal-time courses. DISCUSSION MDR provides effective motion correction of MRR in synthetic and patient data. Future work is needed to compare the performance against other more established methods.
Collapse
Affiliation(s)
- Dimitra Flouri
- Department of Applied Mathematics, University of Leeds, Leeds, UK.
- Department of Biomedical Imaging Sciences, University of Leeds, Leeds, UK.
- School of Biomedical Engineering and Imaging Sciences, Kings College London, London, UK.
- Department of Medical Physics and Biomedical Engineering, University College London, London, UK.
| | - Daniel Lesnic
- Department of Applied Mathematics, University of Leeds, Leeds, UK
| | - Constantina Chrysochou
- Department of Renal Medicine, Salford Royal National Health Service Foundation Trust, Salford, UK
| | - Jehill Parikh
- Translational and Clinical Research Institute, Newcastle University, Newcastle upon Tyne, UK
- Newcastle Magnetic Resonance Centre, Campus for Ageing and Vitality, University of Newcastle, Newcastle upon Tyne, UK
| | - Peter Thelwall
- Translational and Clinical Research Institute, Newcastle University, Newcastle upon Tyne, UK
- Newcastle Magnetic Resonance Centre, Campus for Ageing and Vitality, University of Newcastle, Newcastle upon Tyne, UK
| | - Neil Sheerin
- Translational and Clinical Research Institute, Newcastle University, Newcastle upon Tyne, UK
| | - Philip A Kalra
- Department of Renal Medicine, Salford Royal National Health Service Foundation Trust, Salford, UK
| | - David L Buckley
- Department of Biomedical Imaging Sciences, University of Leeds, Leeds, UK
| | - Steven P Sourbron
- Department of Biomedical Imaging Sciences, University of Leeds, Leeds, UK
- Department of Infection, Immunity and Cardiovascular Disease, University of Sheffield, Sheffield, UK
| |
Collapse
|
5
|
On the use of multicompartment models of diffusion and relaxation for placental imaging. Placenta 2021; 112:197-203. [PMID: 34392172 DOI: 10.1016/j.placenta.2021.07.302] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/11/2021] [Revised: 04/27/2021] [Accepted: 07/27/2021] [Indexed: 12/14/2022]
Abstract
Multi-compartment models of diffusion and relaxation are ubiquitous in magnetic resonance research especially applied to neuroimaging applications. These models are increasingly making their way into the world of placental imaging. This review provides a framework for their motivation and implementation and describes some of the outstanding questions that need to be answered before they can be routinely adopted.
Collapse
|
6
|
Chen J, Hagiwara M, Givi B, Schmidt B, Liu C, Chen Q, Logan J, Mikheev A, Rusinek H, Kim SG. Assessment of metastatic lymph nodes in head and neck squamous cell carcinomas using simultaneous 18F-FDG-PET and MRI. Sci Rep 2020; 10:20764. [PMID: 33247166 PMCID: PMC7695736 DOI: 10.1038/s41598-020-77740-5] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2020] [Accepted: 11/02/2020] [Indexed: 11/09/2022] Open
Abstract
In this study, we investigate the feasibility of using dynamic contrast enhanced magnetic resonance imaging (DCE-MRI), diffusion weighted imaging (DWI), and dynamic positron emission tomography (PET) for detection of metastatic lymph nodes in head and neck squamous cell carcinoma (HNSCC) cases. Twenty HNSCC patients scheduled for lymph node dissection underwent DCE-MRI, dynamic PET, and DWI using a PET-MR scanner within one week prior to their planned surgery. During surgery, resected nodes were labeled to identify their nodal levels and sent for routine clinical pathology evaluation. Quantitative parameters of metastatic and normal nodes were calculated from DCE-MRI (ve, vp, PS, Fp, Ktrans), DWI (ADC) and PET (Ki, K1, k2, k3) to assess if an individual or a combination of parameters can classify normal and metastatic lymph nodes accurately. There were 38 normal and 11 metastatic nodes covered by all three imaging methods and confirmed by pathology. 34% of all normal nodes had volumes greater than or equal to the smallest metastatic node while 4 normal nodes had SUV > 4.5. Among the MRI parameters, the median vp, Fp, PS, and Ktrans values of the metastatic lymph nodes were significantly lower (p = <0.05) than those of normal nodes. ve and ADC did not show any statistical significance. For the dynamic PET parameters, the metastatic nodes had significantly higher k3 (p value = 8.8 × 10-8) and Ki (p value = 5.3 × 10-8) than normal nodes. K1 and k2 did not show any statistically significant difference. Ki had the best separation with accuracy = 0.96 (sensitivity = 1, specificity = 0.95) using a cutoff of Ki = 5.3 × 10-3 mL/cm3/min, while k3 and volume had accuracy of 0.94 (sensitivity = 0.82, specificity = 0.97) and 0.90 (sensitivity = 0.64, specificity = 0.97) respectively. 100% accuracy can be achieved using a multivariate logistic regression model of MRI parameters after thresholding the data with Ki < 5.3 × 10-3 mL/cm3/min. The results of this preliminary study suggest that quantitative MRI may provide additional value in distinguishing metastatic nodes, particularly among small nodes, when used together with FDG-PET.
Collapse
Affiliation(s)
- Jenny Chen
- grid.137628.90000 0004 1936 8753Department of Radiology, Center for Biomedical Imaging (CBI), Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, 660 First Avenue, New York, NY 10016 USA
| | - Mari Hagiwara
- grid.137628.90000 0004 1936 8753Department of Radiology, Center for Biomedical Imaging (CBI), Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, 660 First Avenue, New York, NY 10016 USA
| | - Babak Givi
- grid.137628.90000 0004 1936 8753Department of Otolaryngology-Head and Neck Surgery, New York University School of Medicine, New York, NY USA
| | - Brian Schmidt
- grid.137628.90000 0004 1936 8753Department of Oral and Maxillofacial Surgery, Bluestone Center for Clinical Research, New York University College of Dentistry, New York, NY USA
| | - Cheng Liu
- grid.137628.90000 0004 1936 8753Department of Pathology, New York University School of Medicine, New York, NY USA
| | - Qi Chen
- grid.137628.90000 0004 1936 8753Department of Radiology, Center for Biomedical Imaging (CBI), Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, 660 First Avenue, New York, NY 10016 USA
| | - Jean Logan
- grid.137628.90000 0004 1936 8753Department of Radiology, Center for Biomedical Imaging (CBI), Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, 660 First Avenue, New York, NY 10016 USA
| | - Artem Mikheev
- grid.137628.90000 0004 1936 8753Department of Radiology, Center for Biomedical Imaging (CBI), Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, 660 First Avenue, New York, NY 10016 USA
| | - Henry Rusinek
- grid.137628.90000 0004 1936 8753Department of Radiology, Center for Biomedical Imaging (CBI), Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, 660 First Avenue, New York, NY 10016 USA
| | - Sungheon Gene Kim
- Department of Radiology, Center for Biomedical Imaging (CBI), Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, 660 First Avenue, New York, NY, 10016, USA. .,Department of Radiology, Weill Cornell Medical College, New York, NY, USA.
| |
Collapse
|
7
|
Giménez-Alventosa V, Segrelles JD, Moltó G, Roca-Sogorb M. APRICOT: Advanced Platform for Reproducible Infrastructures in the Cloud via Open Tools. Methods Inf Med 2020; 59:e33-e45. [PMID: 32777825 PMCID: PMC7746519 DOI: 10.1055/s-0040-1712460] [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] [Indexed: 12/02/2022]
Abstract
Background
Scientific publications are meant to exchange knowledge among researchers but the inability to properly reproduce computational experiments limits the quality of scientific research. Furthermore, bibliography shows that irreproducible preclinical research exceeds 50%, which produces a huge waste of resources on nonprofitable research at Life Sciences field. As a consequence, scientific reproducibility is being fostered to promote Open Science through open databases and software tools that are typically deployed on existing computational resources. However, some computational experiments require complex virtual infrastructures, such as elastic clusters of PCs, that can be dynamically provided from multiple clouds. Obtaining these infrastructures requires not only an infrastructure provider, but also advanced knowledge in the cloud computing field.
Objectives
The main aim of this paper is to improve reproducibility in life sciences to produce better and more cost-effective research. For that purpose, our intention is to simplify the infrastructure usage and deployment for researchers.
Methods
This paper introduces Advanced Platform for Reproducible Infrastructures in the Cloud via Open Tools (APRICOT), an open source extension for Jupyter to deploy deterministic virtual infrastructures across multiclouds for reproducible scientific computational experiments. To exemplify its utilization and how APRICOT can improve the reproduction of experiments with complex computation requirements, two examples in the field of life sciences are provided. All requirements to reproduce both experiments are disclosed within APRICOT and, therefore, can be reproduced by the users.
Results
To show the capabilities of APRICOT, we have processed a real magnetic resonance image to accurately characterize a prostate cancer using a Message Passing Interface cluster deployed automatically with APRICOT. In addition, the second example shows how APRICOT scales the deployed infrastructure, according to the workload, using a batch cluster. This example consists of a multiparametric study of a positron emission tomography image reconstruction.
Conclusion
APRICOT's benefits are the integration of specific infrastructure deployment, the management and usage for Open Science, making experiments that involve specific computational infrastructures reproducible. All the experiment steps and details can be documented at the same Jupyter notebook which includes infrastructure specifications, data storage, experimentation execution, results gathering, and infrastructure termination. Thus, distributing the experimentation notebook and needed data should be enough to reproduce the experiment.
Collapse
Affiliation(s)
- Vicent Giménez-Alventosa
- Instituto de Instrumentación para Imagen Molecular (I3M), Centro mixto CSIC-Universitat Politècnica de València, Valencia, Spain
| | - José Damián Segrelles
- Instituto de Instrumentación para Imagen Molecular (I3M), Centro mixto CSIC-Universitat Politècnica de València, Valencia, Spain
| | - Germán Moltó
- Instituto de Instrumentación para Imagen Molecular (I3M), Centro mixto CSIC-Universitat Politècnica de València, Valencia, Spain
| | - Mar Roca-Sogorb
- Quantitative Imaging Biomarkers in Medicine (QUIBIM), Valencia, Spain
| |
Collapse
|
8
|
Ghodasara S, Chen Y, Pahwa S, Griswold MA, Seiberlich N, Wright KL, Gulani V. Quantifying Perfusion Properties with DCE-MRI Using a Dictionary Matching Approach. Sci Rep 2020; 10:10210. [PMID: 32576843 PMCID: PMC7311534 DOI: 10.1038/s41598-020-66985-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2019] [Accepted: 05/25/2020] [Indexed: 12/18/2022] Open
Abstract
Perfusion properties can be estimated from pharmacokinetic models applied to DCE-MRI data using curve fitting algorithms; however, these suffer from drawbacks including the local minimum problem and substantial computational time. Here, a dictionary matching approach is proposed as an alternative. Curve fitting and dictionary matching were applied to simulated data using the dual-input single-compartment model with known perfusion property values and 5 in vivo DCE-MRI datasets. In simulation at SNR 60 dB, the dictionary estimate had a mean percent error of 0.4-1.0% for arterial fraction, 0.5-1.4% for distribution volume, and 0.0% for mean transit time. The curve fitting estimate had a mean percent error of 1.1-2.1% for arterial fraction, 0.5-1.3% for distribution volume, and 0.2-1.8% for mean transit time. In vivo, dictionary matching and curve fitting showed no statistically significant differences in any of the perfusion property measurements in any of the 10 ROIs between the methods. In vivo, the dictionary method performed over 140-fold faster than curve fitting, obtaining whole volume perfusion maps in just over 10 s. This study establishes the feasibility of using a dictionary matching approach as a new and faster way of estimating perfusion properties from pharmacokinetic models in DCE-MRI.
Collapse
Affiliation(s)
- Satyam Ghodasara
- School of Medicine, Case Western Reserve University, Cleveland, Ohio, USA
| | - Yong Chen
- Department of Radiology, University of North Carolina, Chapel Hill, North Carolina, USA
| | - Shivani Pahwa
- School of Medicine, Case Western Reserve University, Cleveland, Ohio, USA
- Department of Radiology, Case Western Reserve University, University Hospitals Cleveland Medical Center, Cleveland, Ohio, USA
| | - Mark A Griswold
- School of Medicine, Case Western Reserve University, Cleveland, Ohio, USA
- Department of Radiology, Case Western Reserve University, University Hospitals Cleveland Medical Center, Cleveland, Ohio, USA
- Department of Biomedical Engineering, Case Western Reserve University, Cleveland, Ohio, USA
| | - Nicole Seiberlich
- Department of Radiology, University of Michigan, Ann Arbor, Michigan, USA
| | - Katherine L Wright
- Department of Radiology, University of Michigan, Ann Arbor, Michigan, USA
| | - Vikas Gulani
- Department of Radiology, University of Michigan, Ann Arbor, Michigan, USA.
| |
Collapse
|
9
|
Dikaios N. Stochastic Gradient Langevin dynamics for joint parameterization of tracer kinetic models, input functions, and T1 relaxation-times from undersampled k-space DCE-MRI. Med Image Anal 2020; 62:101690. [DOI: 10.1016/j.media.2020.101690] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2019] [Revised: 02/16/2020] [Accepted: 03/13/2020] [Indexed: 02/04/2023]
|
10
|
Flouri D, Owen D, Aughwane R, Mufti N, Maksym K, Sokolska M, Kendall G, Bainbridge A, Atkinson D, Vercauteren T, Ourselin S, David AL, Melbourne A. Improved fetal blood oxygenation and placental estimated measurements of diffusion-weighted MRI using data-driven Bayesian modeling. Magn Reson Med 2019; 83:2160-2172. [PMID: 31742785 PMCID: PMC7064949 DOI: 10.1002/mrm.28075] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2019] [Revised: 10/22/2019] [Accepted: 10/22/2019] [Indexed: 12/24/2022]
Abstract
PURPOSE Motion correction in placental DW-MRI is challenging due to maternal breathing motion, maternal movements, and rapid intensity changes. Parameter estimates are usually obtained using least-squares methods for voxel-wise fitting; however, they typically give noisy estimates due to low signal-to-noise ratio. We introduce a model-driven registration (MDR) technique which incorporates a placenta-specific signal model into the registration process, and we present a Bayesian approach for Diffusion-rElaxation Combined Imaging for Detailed placental Evaluation model to obtain individual and population trends in estimated parameters. METHODS MDR exploits the fact that a placenta signal model is available and thus we incorporate it into the registration to generate a series of target images. The proposed registration method is compared to a pre-existing method used for DCE-MRI data making use of principal components analysis. The Bayesian shrinkage prior (BSP) method has no user-defined parameters and therefore measures of parameter variation in a region of interest are determined by the data alone. The MDR method and the Bayesian approach were evaluated on 10 control 4D DW-MRI singleton placental data. RESULTS MDR method improves the alignment of placenta data compared to the pre-existing method. It also shows a further reduction of the residual error between the data and the fit. BSP approach showed higher precision leading to more clearly apparent spatial features in the parameter maps. Placental fetal oxygen saturation (FO2 ) showed a negative linear correlation with gestational age. CONCLUSIONS The proposed pipeline provides a robust framework for registering DW-MRI data and analyzing longitudinal changes of placental function.
Collapse
Affiliation(s)
- Dimitra Flouri
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom.,Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| | - David Owen
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom.,Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| | - Rosalind Aughwane
- Institute for Women's Health, University College Hospital, London, United Kingdom
| | - Nada Mufti
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom.,Institute for Women's Health, University College Hospital, London, United Kingdom
| | - Kasia Maksym
- Institute for Women's Health, University College Hospital, London, United Kingdom
| | | | - Giles Kendall
- Institute for Women's Health, University College Hospital, London, United Kingdom
| | - Alan Bainbridge
- Medical Physics, University College Hospital, London, United Kingdom
| | - David Atkinson
- Centre for Medical Imaging, University College London, London, United Kingdom
| | - Tom Vercauteren
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom.,Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| | - Sebastien Ourselin
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Anna L David
- Institute for Women's Health, University College Hospital, London, United Kingdom.,University Hospital KU Leuven, Leuven, Belgium.,NIHR Biomedical Research Centre, University College London Hospitals, London, United Kingdom
| | - Andrew Melbourne
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom.,Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
| |
Collapse
|
11
|
Ahmed Z, Levesque IR. Pharmacokinetic modeling of dynamic contrast-enhanced MRI using a reference region and input function tail. Magn Reson Med 2019; 83:286-298. [PMID: 31393033 DOI: 10.1002/mrm.27913] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2019] [Revised: 06/18/2019] [Accepted: 06/18/2019] [Indexed: 12/20/2022]
Abstract
PURPOSE Quantitative analysis of dynamic contrast-enhanced MRI (DCE-MRI) requires an arterial input function (AIF) which is difficult to measure. We propose the reference region and input function tail (RRIFT) approach which uses a reference tissue and the washout portion of the AIF. METHODS RRIFT was evaluated in simulations with 100 parameter combinations at various temporal resolutions (5-30 s) and noise levels (σ = 0.01-0.05 mM). RRIFT was compared against the extended Tofts model (ETM) in 8 studies from patients with glioblastoma multiforme. Two versions of RRIFT were evaluated: one using measured patient-specific AIF tails, and another assuming a literature-based AIF tail. RESULTS RRIFT estimated the transfer constant K trans and interstitial volume v e with median errors within 20% across all simulations. RRIFT was more accurate and precise than the ETM at temporal resolutions slower than 10 s. The percentage error of K trans had a median and interquartile range of -9 ± 45% with the ETM and -2 ± 17% with RRIFT at a temporal resolution of 30 s under noiseless conditions. RRIFT was in excellent agreement with the ETM in vivo, with concordance correlation coefficients (CCC) of 0.95 for K trans , 0.96 for v e , and 0.73 for the plasma volume v p using a measured AIF tail. With the literature-based AIF tail, the CCC was 0.89 for K trans , 0.93 for v e and 0.78 for v p . CONCLUSIONS Quantitative DCE-MRI analysis using the input function tail and a reference tissue yields absolute kinetic parameters with the RRIFT method. This approach was viable in simulation and in vivo for temporal resolutions as low as 30 s.
Collapse
Affiliation(s)
- Zaki Ahmed
- Medical Physics Unit, McGill University, Montreal, Canada.,Department of Physics, McGill University, Montreal, Canada
| | - Ives R Levesque
- Medical Physics Unit, McGill University, Montreal, Canada.,Department of Physics, McGill University, Montreal, Canada.,Gerald Bronfman Department of Oncology, McGill University, Montreal, Canada.,Cancer Research Program, Research Institute of the McGill University Health Centre, Montreal, Canada
| |
Collapse
|
12
|
Ahmed Z, Levesque IR. An extended reference region model for DCE-MRI that accounts for plasma volume. NMR IN BIOMEDICINE 2018; 31:e3924. [PMID: 29745982 DOI: 10.1002/nbm.3924] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/30/2017] [Revised: 02/20/2018] [Accepted: 02/27/2018] [Indexed: 06/08/2023]
Abstract
The reference region model (RRM) for dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) provides pharmacokinetic parameters without requiring the arterial input function. A limitation of the RRM is that it assumes that the blood plasma volume in the tissue of interest is zero, but this is often not true in highly vascularized tissues, such as some tumours. This study proposes an extended reference region model (ERRM) to account for tissue plasma volume. Furthermore, ERRM was combined with a two-fit approach to reduce the number of fitting parameters, and this was named the constrained ERRM (CERRM). The accuracy and precision of RRM, ERRM and CERRM were evaluated in simulations covering a range of parameters, noise and temporal resolutions. These models were also compared with the extended Tofts model (ETM) on in vivo glioblastoma multiforme data. In simulations, RRM overestimated Ktrans by over 10% at vp = 0.01 under noiseless conditions. In comparison, ERRM and CERRM were both accurate, with CERRM showing better precision when noise was included. On in vivo data, CERRM provided maps that had the highest agreement with ETM, whilst also being robust at temporal resolutions as poor as 30 s. ERRM can provide pharmacokinetic parameters without an arterial input function in tissues with non-negligible vp where RRM provides inaccurate estimates. The two-fit approach, named CERRM, further improves on the accuracy and precision of ERRM.
Collapse
Affiliation(s)
- Zaki Ahmed
- Medical Physics Unit, McGill University, Montreal, QC, Canada
- Department of Physics, McGill University, Montreal, QC, Canada
| | - Ives R Levesque
- Medical Physics Unit, McGill University, Montreal, QC, Canada
- Department of Physics, McGill University, Montreal, QC, Canada
- Research Institute of the McGill University Health Centre, Montreal, QC, Canada
| |
Collapse
|
13
|
Kargar S, Borisch EA, Froemming AT, Kawashima A, Mynderse LA, Stinson EG, Trzasko JD, Riederer SJ. Robust and efficient pharmacokinetic parameter non-linear least squares estimation for dynamic contrast enhanced MRI of the prostate. Magn Reson Imaging 2017; 48:50-61. [PMID: 29278764 DOI: 10.1016/j.mri.2017.12.021] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2017] [Revised: 12/09/2017] [Accepted: 12/21/2017] [Indexed: 12/15/2022]
Abstract
PURPOSE To describe an efficient numerical optimization technique using non-linear least squares to estimate perfusion parameters for the Tofts and extended Tofts models from dynamic contrast enhanced (DCE) MRI data and apply the technique to prostate cancer. METHODS Parameters were estimated by fitting the two Tofts-based perfusion models to the acquired data via non-linear least squares. We apply Variable Projection (VP) to convert the fitting problem from a multi-dimensional to a one-dimensional line search to improve computational efficiency and robustness. Using simulation and DCE-MRI studies in twenty patients with suspected prostate cancer, the VP-based solver was compared against the traditional Levenberg-Marquardt (LM) strategy for accuracy, noise amplification, robustness to converge, and computation time. RESULTS The simulation demonstrated that VP and LM were both accurate in that the medians closely matched assumed values across typical signal to noise ratio (SNR) levels for both Tofts models. VP and LM showed similar noise sensitivity. Studies using the patient data showed that the VP method reliably converged and matched results from LM with approximate 3× and 2× reductions in computation time for the standard (two-parameter) and extended (three-parameter) Tofts models. While LM failed to converge in 14% of the patient data, VP converged in the ideal 100%. CONCLUSION The VP-based method for non-linear least squares estimation of perfusion parameters for prostate MRI is equivalent in accuracy and robustness to noise, while being more reliably (100%) convergent and computationally about 3× (TM) and 2× (ETM) faster than the LM-based method.
Collapse
Affiliation(s)
- Soudabeh Kargar
- Biomedical Engineering and Physiology Program, Mayo Graduate School, Rochester, MN, United States; Department of Radiology, Mayo Clinic, Rochester, MN, United States
| | - Eric A Borisch
- Department of Radiology, Mayo Clinic, Rochester, MN, United States
| | - Adam T Froemming
- Department of Radiology, Mayo Clinic, Rochester, MN, United States
| | - Akira Kawashima
- Department of Radiology, Mayo Clinic, Scottsdale, AZ, United States
| | - Lance A Mynderse
- Department of Urology, Mayo Clinic, Rochester, MN, United States
| | - Eric G Stinson
- Department of Radiology, Mayo Clinic, Rochester, MN, United States
| | - Joshua D Trzasko
- Department of Radiology, Mayo Clinic, Rochester, MN, United States
| | - Stephen J Riederer
- Biomedical Engineering and Physiology Program, Mayo Graduate School, Rochester, MN, United States; Department of Radiology, Mayo Clinic, Rochester, MN, United States.
| |
Collapse
|
14
|
Guo Y, Lingala SG, Bliesener Y, Lebel RM, Zhu Y, Nayak KS. Joint arterial input function and tracer kinetic parameter estimation from undersampled dynamic contrast-enhanced MRI using a model consistency constraint. Magn Reson Med 2017; 79:2804-2815. [PMID: 28905411 DOI: 10.1002/mrm.26904] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2017] [Revised: 08/11/2017] [Accepted: 08/16/2017] [Indexed: 12/13/2022]
Abstract
PURPOSE To develop and evaluate a model-based reconstruction framework for joint arterial input function (AIF) and kinetic parameter estimation from undersampled brain tumor dynamic contrast-enhanced MRI (DCE-MRI) data. METHODS The proposed method poses the tracer-kinetic (TK) model as a model consistency constraint, enabling the flexible inclusion of different TK models and TK solvers, and the joint estimation of the AIF. The proposed method is evaluated using an anatomic realistic digital reference object (DRO), and nine retrospectively down-sampled brain tumor DCE-MRI datasets. We also demonstrate application to 30-fold prospectively undersampled brain tumor DCE-MRI. RESULTS In DRO studies with up to 60-fold undersampling, the proposed method provided TK maps with low error that were comparable to fully sampled data and were demonstrated to be compatible with a third-party TK solver. In retrospective undersampling studies, this method provided patient-specific AIF with normalized root mean-squared-error (normalized by the 90th percentile value) less than 8% at up to 100-fold undersampling. In the 30-fold undersampled prospective study, the proposed method provided high-resolution whole-brain TK maps and patient-specific AIF. CONCLUSION The proposed model-based DCE-MRI reconstruction enables the use of different TK solvers with a model consistency constraint and enables joint estimation of patient-specific AIF. TK maps and patient-specific AIF with high fidelity can be reconstructed at up to 100-fold undersampling in k,t-space. Magn Reson Med 79:2804-2815, 2018. © 2017 International Society for Magnetic Resonance in Medicine.
Collapse
Affiliation(s)
- Yi Guo
- Ming Hsieh Department of Electrical Engineering, Viterbi School of Engineering, University of Southern California, Los Angeles, California, USA
| | - Sajan Goud Lingala
- Ming Hsieh Department of Electrical Engineering, Viterbi School of Engineering, University of Southern California, Los Angeles, California, USA
| | - Yannick Bliesener
- Ming Hsieh Department of Electrical Engineering, Viterbi School of Engineering, University of Southern California, Los Angeles, California, USA
| | | | - Yinghua Zhu
- Ming Hsieh Department of Electrical Engineering, Viterbi School of Engineering, University of Southern California, Los Angeles, California, USA
| | - Krishna S Nayak
- Ming Hsieh Department of Electrical Engineering, Viterbi School of Engineering, University of Southern California, Los Angeles, California, USA
| |
Collapse
|
15
|
Jafari R, Chhabra S, Prince MR, Wang Y, Spincemaille P. Vastly accelerated linear least-squares fitting with numerical optimization for dual-input delay-compensated quantitative liver perfusion mapping. Magn Reson Med 2017; 79:2415-2421. [PMID: 28833534 DOI: 10.1002/mrm.26888] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2017] [Revised: 08/02/2017] [Accepted: 08/03/2017] [Indexed: 12/12/2022]
Abstract
PURPOSE To propose an efficient algorithm to perform dual input compartment modeling for generating perfusion maps in the liver. METHODS We implemented whole field-of-view linear least squares (LLS) to fit a delay-compensated dual-input single-compartment model to very high temporal resolution (four frames per second) contrast-enhanced 3D liver data, to calculate kinetic parameter maps. Using simulated data and experimental data in healthy subjects and patients, whole-field LLS was compared with the conventional voxel-wise nonlinear least-squares (NLLS) approach in terms of accuracy, performance, and computation time. RESULTS Simulations showed good agreement between LLS and NLLS for a range of kinetic parameters. The whole-field LLS method allowed generating liver perfusion maps approximately 160-fold faster than voxel-wise NLLS, while obtaining similar perfusion parameters. CONCLUSIONS Delay-compensated dual-input liver perfusion analysis using whole-field LLS allows generating perfusion maps with a considerable speedup compared with conventional voxel-wise NLLS fitting. Magn Reson Med 79:2415-2421, 2018. © 2017 International Society for Magnetic Resonance in Medicine.
Collapse
Affiliation(s)
- Ramin Jafari
- Meinig School of Biomedical Engineering, Cornell University, Ithaca, New York, USA
| | - Shalini Chhabra
- Department of Radiology, Weill Cornell Medicine, New York, New York, USA
| | - Martin R Prince
- Department of Radiology, Weill Cornell Medicine, New York, New York, USA
| | - Yi Wang
- Meinig School of Biomedical Engineering, Cornell University, Ithaca, New York, USA.,Department of Radiology, Weill Cornell Medicine, New York, New York, USA
| | | |
Collapse
|
16
|
Garpebring A, Löfstedt T. Parameter estimation using weighted total least squares in the two-compartment exchange model. Magn Reson Med 2017; 79:561-567. [PMID: 28349618 DOI: 10.1002/mrm.26677] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2016] [Revised: 02/02/2017] [Accepted: 02/21/2017] [Indexed: 12/14/2022]
Abstract
PURPOSE The linear least squares (LLS) estimator provides a fast approach to parameter estimation in the linearized two-compartment exchange model. However, the LLS method may introduce a bias through correlated noise in the system matrix of the model. The purpose of this work is to present a new estimator for the linearized two-compartment exchange model that takes this noise into account. METHOD To account for the noise in the system matrix, we developed an estimator based on the weighted total least squares (WTLS) method. Using simulations, the proposed WTLS estimator was compared, in terms of accuracy and precision, to an LLS estimator and a nonlinear least squares (NLLS) estimator. RESULTS The WTLS method improved the accuracy compared to the LLS method to levels comparable to the NLLS method. This improvement was at the expense of increased computational time; however, the WTLS was still faster than the NLLS method. At high signal-to-noise ratio all methods provided similar precisions while inconclusive results were observed at low signal-to-noise ratio. CONCLUSION The proposed method provides improvements in accuracy compared to the LLS method, however, at an increased computational cost. Magn Reson Med 79:561-567, 2017. © 2017 International Society for Magnetic Resonance in Medicine.
Collapse
Affiliation(s)
| | - Tommy Löfstedt
- Department of Radiation Sciences, Umeå University, Umeå, Sweden
| |
Collapse
|
17
|
Kallehauge JF, Sourbron S, Irving B, Tanderup K, Schnabel JA, Chappell MA. Comparison of linear and nonlinear implementation of the compartmental tissue uptake model for dynamic contrast-enhanced MRI. Magn Reson Med 2016; 77:2414-2423. [PMID: 27605429 PMCID: PMC5484345 DOI: 10.1002/mrm.26324] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2016] [Revised: 05/10/2016] [Accepted: 06/08/2016] [Indexed: 12/14/2022]
Abstract
Purpose Fitting tracer kinetic models using linear methods is much faster than using their nonlinear counterparts, although this comes often at the expense of reduced accuracy and precision. The aim of this study was to derive and compare the performance of the linear compartmental tissue uptake (CTU) model with its nonlinear version with respect to their percentage error and precision. Theory and Methods The linear and nonlinear CTU models were initially compared using simulations with varying noise and temporal sampling. Subsequently, the clinical applicability of the linear model was demonstrated on 14 patients with locally advanced cervical cancer examined with dynamic contrast‐enhanced magnetic resonance imaging. Results Simulations revealed equal percentage error and precision when noise was within clinical achievable ranges (contrast‐to‐noise ratio >10). The linear method was significantly faster than the nonlinear method, with a minimum speedup of around 230 across all tested sampling rates. Clinical analysis revealed that parameters estimated using the linear and nonlinear CTU model were highly correlated (ρ ≥ 0.95). Conclusion The linear CTU model is computationally more efficient and more stable against temporal downsampling, whereas the nonlinear method is more robust to variations in noise. The two methods may be used interchangeably within clinical achievable ranges of temporal sampling and noise. Magn Reson Med 77:2414–2423, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.
Collapse
Affiliation(s)
- Jesper F Kallehauge
- Institute of Biomedical Engineering, Department of Engineering Science University of Oxford, Oxford, United Kingdom
| | - Steven Sourbron
- Division of Biomedical Imaging, University of Leeds, Leeds, United Kingdom
| | - Benjamin Irving
- Institute of Biomedical Engineering, Department of Engineering Science University of Oxford, Oxford, United Kingdom
| | - Kari Tanderup
- Department of Oncology, Aarhus University Hospital, Aarhus, Denmark
| | - Julia A Schnabel
- Division of Imaging Science and Biomedical Engineering, King's College London, London, United Kingdom
| | - Michael A Chappell
- Institute of Biomedical Engineering, Department of Engineering Science University of Oxford, Oxford, United Kingdom
| |
Collapse
|
18
|
Zöllner FG, Daab M, Sourbron SP, Schad LR, Schoenberg SO, Weisser G. An open source software for analysis of dynamic contrast enhanced magnetic resonance images: UMMPerfusion revisited. BMC Med Imaging 2016; 16:7. [PMID: 26767969 PMCID: PMC4712457 DOI: 10.1186/s12880-016-0109-0] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2015] [Accepted: 01/06/2016] [Indexed: 12/11/2022] Open
Abstract
Background Perfusion imaging has become an important image based tool to derive the physiological information in various applications, like tumor diagnostics and therapy, stroke, (cardio-) vascular diseases, or functional assessment of organs. However, even after 20 years of intense research in this field, perfusion imaging still remains a research tool without a broad clinical usage. One problem is the lack of standardization in technical aspects which have to be considered for successful quantitative evaluation; the second problem is a lack of tools that allow a direct integration into the diagnostic workflow in radiology. Results Five compartment models, namely, a one compartment model (1CP), a two compartment exchange (2CXM), a two compartment uptake model (2CUM), a two compartment filtration model (2FM) and eventually the extended Toft’s model (ETM) were implemented as plugin for the DICOM workstation OsiriX. Moreover, the plugin has a clean graphical user interface and provides means for quality management during the perfusion data analysis. Based on reference test data, the implementation was validated against a reference implementation. No differences were found in the calculated parameters. Conclusion We developed open source software to analyse DCE-MRI perfusion data. The software is designed as plugin for the DICOM Workstation OsiriX. It features a clean GUI and provides a simple workflow for data analysis while it could also be seen as a toolbox providing an implementation of several recent compartment models to be applied in research tasks. Integration into the infrastructure of a radiology department is given via OsiriX. Results can be saved automatically and reports generated automatically during data analysis ensure certain quality control.
Collapse
Affiliation(s)
- Frank G Zöllner
- Computer Assisted Clinical Medicine, Medical Faculty Mannheim, Heidelberg University, Theodor-Kutzer-Ufer 1-3, 68167, Mannheim, Germany.
| | - Markus Daab
- Computer Assisted Clinical Medicine, Medical Faculty Mannheim, Heidelberg University, Theodor-Kutzer-Ufer 1-3, 68167, Mannheim, Germany.
| | | | - Lothar R Schad
- Computer Assisted Clinical Medicine, Medical Faculty Mannheim, Heidelberg University, Theodor-Kutzer-Ufer 1-3, 68167, Mannheim, Germany.
| | - Stefan O Schoenberg
- Department of Clinical Radiology and Nuclear Medicine, University Medical Center Mannheim, Heidelberg University, Mannheim, Germany.
| | - Gerald Weisser
- Department of Clinical Radiology and Nuclear Medicine, University Medical Center Mannheim, Heidelberg University, Mannheim, Germany.
| |
Collapse
|