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.

InVAErt networks for amortized inference and identifiability analysis of lumped parameter hemodynamic models

InVAErt networks: A data-driven framework for model synthesis and identifiability analysis

Lnista-Net: Learning Nonlocal Regularized Iterative Shrinkage-Thresholding Network for Noninvasive Cardiac Transmembrane Potential Imaging.

Towards Enabling Cardiac Digital Twins of Myocardial Infarction Using Deep Computational Models for Inverse Inference

Inference of ventricular activation properties from non-invasive electrocardiography

Simulation-based Inference for Cardiovascular Models

Real-time whole-heart electromechanical simulations using Latent Neural Ordinary Differential Equations

Influence of Resting Venous Blood Volume Fraction on Dynamic Causal Modeling and System Identifiability

Machine learning in cardiovascular flows modeling: Predicting arterial blood pressure from non-invasive 4D flow MRI data using physics-informed neural networks

Physics-informed graph neural networks for flow field estimation in carotid arteries

Cardiovascular system identification: Simulation study using arterial and central venous pressures

Machine Learning Identification Framework of Hemodynamics of Blood Flow in Patient-Specific Coronary Arteries with Abnormality

HyPer-EP: Meta-Learning Hybrid Personalized Models for Cardiac Electrophysiology

Rapid Estimation of Left Ventricular Contractility with a Physics-Informed Neural Network Inverse Modeling Approach

Multimodal Variational Autoencoder for Low-cost Cardiac Hemodynamics Instability Detection

Identifiability of spatiotemporal tissue perfusion models

Variational Inverting Network for Statistical Inverse Problems of Partial Differential Equations

A machine learning method for real-time numerical simulations of cardiac electromechanics

Physics-informed neural networks (PINNs) for 4D hemodynamics prediction: An investigation of optimal framework based on vascular morphology