Efficient Estimation of Cardiac Conductivities via POD-DEIM Model Order Reduction

Clinical oriented applications of computational electrocardiology require efficient and reliable identification of patient-specific parameters of mathematical models based on available measures. In particular, the estimation of cardiac conductivities in models of potential propagation is crucial, since they have major quantitative impact on the solution. Available estimates of cardiac conductivities are significantly diverse in the literature and the definition of experimental/mathematical estimation techniques is an open problem with important practical implications in clinics. We have recently proposed a methodology based on a variational procedure, where the reliability is confirmed by numerical experiments. In this paper we explore model-order-reduction techniques to fit the estimation procedure into timelines of clinical interest. Specifically we consider the Monodomain model and resort to Proper Orthogonal Decomposition (POD) techniques to take advantage of an off-line step when solving iteratively the electrocardiological forward model online. In addition, we perform the Discrete Empirical Interpolation Method (DEIM) to tackle the nonlinearity of the model. While standard POD techniques usually fail in this kind of problems, due to the wave-front propagation dynamics, an educated novel sampling of the parameter space based on the concept of Domain of Effectiveness introduced here dramatically reduces the computational cost of the inverse solver by at least 95%.

[1]  L. Clerc.  Directional differences of impulse spread in trabecular muscle from mammalian heart. , 1976, The Journal of physiology.

[2]  O. C. Zienkiewicz,et al.  The Finite Element Method: Its Basis and Fundamentals , 2005 .

[3]  Luca Gerardo-Giorda,et al.  A model-based block-triangular preconditioner for the Bidomain system in electrocardiology , 2009, J. Comput. Phys..

[4]  Kawal S. Rhode,et al.  Personalization of Atrial Anatomy and Electrophysiology as a Basis for Clinical Modeling of Radio-Frequency Ablation of Atrial Fibrillation , 2013, IEEE Transactions on Medical Imaging.

[5]  A. Patera,et al.  A PRIORI CONVERGENCE OF THE GREEDY ALGORITHM FOR THE PARAMETRIZED REDUCED BASIS METHOD , 2012 .

[6]  Charles Pierre,et al.  Comparing the bidomain and monodomain models in electro-cardiology through convergence analysis , 2010 .

[7]  Marcelo Lobosco,et al.  Simulation of Ectopic Pacemakers in the Heart: Multiple Ectopic Beats Generated by Reentry inside Fibrotic Regions , 2015, BioMed research international.

[8]  Jean-Frédéric Gerbeau,et al.  Parameter Identification in Cardiac Electrophysiology Using Proper Orthogonal Decomposition Method , 2011, FIMH.

[9]  Aslak Tveito,et al.  Optimal monodomain approximations of the bidomain equations , 2007, Appl. Math. Comput..

[10]  Michal Rewienski,et al.  A trajectory piecewise-linear approach to model order reduction of nonlinear dynamical systems , 2003 .

[11]  G. Plank,et al.  A Novel Rule-Based Algorithm for Assigning Myocardial Fiber Orientation to Computational Heart Models , 2012, Annals of Biomedical Engineering.

[12]  Stefan Volkwein,et al.  A posteriori error estimation for semilinear parabolic optimal control problems with application to model reduction by POD , 2013 .

[13]  A. McCulloch,et al.  Patient-specific modeling of ventricular activation pattern using surface ECG-derived vectorcardiogram in bundle branch block. , 2014, Progress in biophysics and molecular biology.

[14]  F. Fenton,et al.  Visualization of spiral and scroll waves in simulated and experimental cardiac tissue , 2008 .

[15]  R. Ragno,et al.  Effects of Mentha suaveolens Essential Oil on Chlamydia trachomatis , 2015, BioMed research international.

[16]  Gianluigi Rozza,et al.  Certified reduced basis approximation for parametrized partial differential equations and applications , 2011 .

[17]  P. C. Franzone,et al.  A PARALLEL SOLVER FOR REACTION-DIFFUSION SYSTEMS IN COMPUTATIONAL ELECTROCARDIOLOGY , 2004 .

[18]  Andrew J. Pullan,et al.  Mathematically Modelling the Electrical Activity of the Heart: From Cell to Body Surface and Back Again , 2005 .

[19]  N. Nguyen,et al.  REDUCED BASIS APPROXIMATION AND A POSTERIORI ERROR ESTIMATION FOR THE PARAMETRIZED UNSTEADY BOUSSINESQ EQUATIONS , 2011 .

[20]  A. M. Scher,et al.  Effect of Tissue Anisotropy on Extracellular Potential Fields in Canine Myocardium in Situ , 1982, Circulation research.

[21]  A. McCulloch,et al.  A collocation-Galerkin finite element model of cardiac action potential propagation , 1994, IEEE Transactions on Biomedical Engineering.

[22]  Karen Willcox,et al.  Surrogate and reduced-order modeling: a comparison of approaches for large-scale statistical inverse problems [Chapter 7] , 2010 .

[23]  J. Peraire,et al.  A ‘best points’ interpolation method for efficient approximation of parametrized functions , 2008 .

[24]  A. Veneziani,et al.  Estimation of cardiac conductivities in ventricular tissue by a variational approach , 2015 .

[25]  RewieÅ ski,et al.  A trajectory piecewise-linear approach to model order reduction of nonlinear dynamical systems , 2003 .

[26]  Charbel Farhat,et al.  Progressive construction of a parametric reduced‐order model for PDE‐constrained optimization , 2014, ArXiv.

[27]  Joël M. H. Karel,et al.  Inverse reconstruction of epicardial potentials improved by vectorcardiography and realistic potentials , 2013, Computing in Cardiology 2013.

[28]  Karen Willcox,et al.  A Survey of Projection-Based Model Reduction Methods for Parametric Dynamical Systems , 2015, SIAM Rev..

[29]  A. Veneziani,et al.  A model reduction approach for the variational estimation of vascular compliance by solving an inverse fluid–structure interaction problem , 2014 .

[30]  Natalia Trayanova,et al.  Defibrillation of the heart: insights into mechanisms from modelling studies , 2006, Experimental physiology.

[31]  Qiang Du,et al.  Centroidal Voronoi Tessellations: Applications and Algorithms , 1999, SIAM Rev..

[32]  E. Sachs,et al.  Trust-region proper orthogonal decomposition for flow control , 2000 .

[33]  Stefan Volkwein,et al.  Galerkin Proper Orthogonal Decomposition Methods for a General Equation in Fluid Dynamics , 2002, SIAM J. Numer. Anal..

[34]  M Boulakia,et al.  Reduced-order modeling for cardiac electrophysiology. Application to parameter identification. , 2012, International journal for numerical methods in biomedical engineering.

[35]  David Kilpatrick,et al.  Estimation of the Bidomain Conductivity Parameters of Cardiac Tissue From Extracellular Potential Distributions Initiated by Point Stimulation , 2010, Annals of Biomedical Engineering.

[36]  Cesare Corrado,et al.  Identification of weakly coupled multiphysics problems. Application to the inverse problem of electrocardiography , 2015, J. Comput. Phys..

[37]  Peter R Johnston,et al.  A sensitivity study of conductivity values in the passive bidomain equation. , 2011, Mathematical biosciences.

[38]  David Galbally,et al.  Non‐linear model reduction for uncertainty quantification in large‐scale inverse problems , 2009 .

[39]  Richard J. Beckman,et al.  A Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output From a Computer Code , 2000, Technometrics.

[40]  Stefan Volkwein,et al.  A-Posteriori Error Estimation of Discrete POD Models for PDE-Constrained Optimal Control , 2017 .

[41]  N. Nguyen,et al.  An ‘empirical interpolation’ method: application to efficient reduced-basis discretization of partial differential equations , 2004 .

[42]  K. Willcox,et al.  Aerodynamic Data Reconstruction and Inverse Design Using Proper Orthogonal Decomposition , 2004 .

[43]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[44]  Jeroen J. Bax,et al.  Assessment of left ventricular mechanical dyssynchrony by phase analysis of ECG-gated SPECT myocardial perfusion imaging , 2008, Journal of nuclear cardiology : official publication of the American Society of Nuclear Cardiology.

[45]  Huanhuan Yang.  Parameter Estimation and Reduced-order Modeling in Electrocardiology , 2015 .

[46]  Charbel Farhat,et al.  A Compact Proper Orthogonal Decomposition Basis for Optimization-Oriented Reduced-Order Models , 2008 .

[47]  Danny C. Sorensen,et al.  Nonlinear Model Reduction via Discrete Empirical Interpolation , 2010, SIAM J. Sci. Comput..