Efficient probabilistic model personalization integrating uncertainty on data and parameters: Application to eikonal-diffusion models in cardiac electrophysiology.

Biophysical models are increasingly used for medical applications at the organ scale. However, model predictions are rarely associated with a confidence measure although there are important sources of uncertainty in computational physiology methods. For instance, the sparsity and noise of the clinical data used to adjust the model parameters (personalization), and the difficulty in modeling accurately soft tissue physiology. The recent theoretical progresses in stochastic models make their use computationally tractable, but there is still a challenge in estimating patient-specific parameters with such models. In this work we propose an efficient Bayesian inference method for model personalization using polynomial chaos and compressed sensing. This method makes Bayesian inference feasible in real 3D modeling problems. We demonstrate our method on cardiac electrophysiology. We first present validation results on synthetic data, then we apply the proposed method to clinical data. We demonstrate how this can help in quantifying the impact of the data characteristics on the personalization (and thus prediction) results. Described method can be beneficial for the clinical use of personalized models as it explicitly takes into account the uncertainties on the data and the model parameters while still enabling simulations that can be used to optimize treatment. Such uncertainty handling can be pivotal for the proper use of modeling as a clinical tool, because there is a crucial requirement to know the confidence one can have in personalized models.

[1]  Yuan He,et al.  Reconstructing parameters of the FitzHugh–Nagumo system from boundary potential measurements , 2007, Journal of Computational Neuroscience.

[2]  Heng Tao Shen,et al.  Dimensionality Reduction , 2009, Encyclopedia of Database Systems.

[3]  Y. Marzouk,et al.  A stochastic collocation approach to Bayesian inference in inverse problems , 2009 .

[4]  Alejandro F. Frangi,et al.  Personalization of Fast Conduction Purkinje System in Eikonal-Based Electrophysiological Models with Optical Mapping Data , 2010, STACOM/CESC.

[5]  M. Lustig,et al.  Compressed Sensing MRI , 2008, IEEE Signal Processing Magazine.

[6]  Mats G. Larson,et al.  The Finite Element Method: Theory, Implementation, and Applications , 2013 .

[7]  Daniel Steven,et al.  Epicardial ventricular tachycardia ablation a multicenter safety study. , 2010, Journal of the American College of Cardiology.

[8]  Roger G. Ghanem,et al.  On the construction and analysis of stochastic models: Characterization and propagation of the errors associated with limited data , 2006, J. Comput. Phys..

[9]  Nicholas Zabaras,et al.  An efficient Bayesian inference approach to inverse problems based on an adaptive sparse grid collocation method , 2009 .

[10]  Andrew J. Pullan,et al.  A Finite Element Method for an Eikonal Equation Model of Myocardial Excitation Wavefront Propagation , 2002, SIAM J. Appl. Math..

[11]  D. Xiu Efficient collocational approach for parametric uncertainty analysis , 2007 .

[12]  Y. Rudy,et al.  Noninvasive electrocardiographic imaging for cardiac electrophysiology and arrhythmia , 2004, Nature Medicine.

[13]  Habib N. Najm,et al.  Dimensionality reduction and polynomial chaos acceleration of Bayesian inference in inverse problems , 2008, J. Comput. Phys..

[14]  S Niederer,et al.  Inter-model consistency and complementarity: learning from ex-vivo imaging and electrophysiological data towards an integrated understanding of cardiac physiology. , 2011, Progress in biophysics and molecular biology.

[15]  B M Steinhaus,et al.  Estimating Cardiac Transmembrane Activation and Recovery Times From Unipolar and Bipolar Extracellular Electrograms: A Simulation Study , 1989, Circulation research.

[16]  D. Xiu Fast numerical methods for stochastic computations: A review , 2009 .

[17]  Hervé Delingette,et al.  A Recursive Anisotropic Fast Marching Approach to Reaction Diffusion Equation: Application to Tumor Growth Modeling , 2007, IPMI.

[18]  Huafeng Liu,et al.  Noninvasive Computational Imaging of Cardiac Electrophysiology for 3-D Infarct , 2011, IEEE Transactions on Biomedical Engineering.

[19]  D. Xiu,et al.  Modeling Uncertainty in Steady State Diffusion Problems via Generalized Polynomial Chaos , 2002 .

[20]  Roy C. P. Kerckhoffs,et al.  Patient-specific modeling of dyssynchronous heart failure: a case study. , 2011, Progress in biophysics and molecular biology.

[21]  H Zhang,et al.  Models of cardiac tissue electrophysiology: progress, challenges and open questions. , 2011, Progress in biophysics and molecular biology.

[22]  Robert Michael Kirby,et al.  Application of Stochastic Finite Element Methods to Study the Sensitivity of ECG Forward Modeling to Organ Conductivity , 2008, IEEE Transactions on Biomedical Engineering.

[23]  Habib N. Najm,et al.  Stochastic spectral methods for efficient Bayesian solution of inverse problems , 2005, J. Comput. Phys..

[24]  P. C. Franzone,et al.  Wavefront propagation in an activation model of the anisotropic cardiac tissue: asymptotic analysis and numerical simulations , 1990, Journal of mathematical biology.

[25]  Henggui Zhang,et al.  Cardiac cell modelling: observations from the heart of the cardiac physiome project. , 2011, Progress in biophysics and molecular biology.

[26]  Hervé Delingette,et al.  Coupled personalization of cardiac electrophysiology models for prediction of ischaemic ventricular tachycardia , 2011, Interface Focus.

[27]  H. Halperin,et al.  Integrated Electroanatomic Mapping With Three-Dimensional Computed Tomographic Images for Real-Time Guided Ablations , 2006, Circulation.

[28]  Habib N. Najm,et al.  Stochastic spectral methods for efficient Bayesian solution of inverse problems , 2005 .

[29]  E. Somersalo,et al.  Statistical and computational inverse problems , 2004 .

[30]  Simon R. Arridge,et al.  Model-Based Imaging of Cardiac Apparent Conductivity and Local Conduction Velocity for Diagnosis and Planning of Therapy , 2008, IEEE Transactions on Medical Imaging.

[31]  Eric Voth,et al.  Validation of Computed Tomography Image Integration into the EnSite NavX Mapping System to Perform Catheter Ablation of Atrial Fibrillation , 2008, Journal of cardiovascular electrophysiology.

[32]  A. Morelli Inverse Problem Theory , 2010 .

[33]  Houman Owhadi,et al.  A non-adapted sparse approximation of PDEs with stochastic inputs , 2010, J. Comput. Phys..

[34]  Fabio Nobile,et al.  A Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data , 2008, SIAM J. Numer. Anal..

[35]  M. Girolami,et al.  Riemann manifold Langevin and Hamiltonian Monte Carlo methods , 2011, Journal of the Royal Statistical Society: Series B (Statistical Methodology).

[36]  N. Rashevsky,et al.  Mathematical biology , 1961, Connecticut medicine.

[37]  Jürgen Kurths,et al.  Nonlinear Dynamical System Identification from Uncertain and Indirect Measurements , 2004, Int. J. Bifurc. Chaos.

[38]  E. Somersalo,et al.  Approximation errors and model reduction with an application in optical diffusion tomography , 2006 .

[39]  Masashi Sugiyama,et al.  Dual-Augmented Lagrangian Method for Efficient Sparse Reconstruction , 2009, IEEE Signal Processing Letters.

[40]  Hervé Delingette,et al.  A multi-front eikonal model of cardiac electrophysiology for interactive simulation of radio-frequency ablation , 2011, Comput. Graph..