Physics-Informed Neural Networks for Cardiac Activation Mapping

A critical procedure in diagnosing atrial fibrillation is the creation of electro-anatomic activation maps. Current methods generate these mappings from interpolation using a few sparse data points recorded inside the atria; they neither include prior knowledge of the underlying physics nor uncertainty of these recordings. Here we propose a physics-informed neural network for cardiac activation mapping that accounts for the underlying wave propagation dynamics and we quantify the epistemic uncertainty associated with these predictions. These uncertainty estimates not only allow us to quantify the predictive error of the neural network, but also help to reduce it by judiciously selecting new informative measurement locations via active learning. We illustrate the potential of our approach using a synthetic benchmark problem and a personalized electrophysiology model of the left atrium. We show that our new method outperforms linear interpolation and Gaussian process regression for the benchmark problem and linear interpolation at clinical densities for the left atrium. In both cases, the active learning algorithm achieves lower error levels than random allocation. Our findings open the door toward physics-based electro-anatomic mapping with the ultimate goals to reduce procedural time and improve diagnostic predictability for patients affected by atrial fibrillation. Open source code is available at https://github.com/fsahli/EikonalNet.

[1]  Jeremy E. Oakley,et al.  Probabilistic Interpolation of Uncertain Local Activation Times on Human Atrial Manifolds , 2020, IEEE Transactions on Biomedical Engineering.

[2]  S. De,et al.  Integrating machine learning and multiscale modeling—perspectives, challenges, and opportunities in the biological, biomedical, and behavioral sciences , 2019, npj Digital Medicine.

[3]  K. Kunisch,et al.  Inverse localization of earliest cardiac activation sites from activation maps based on the viscous Eikonal equation , 2019, Journal of Mathematical Biology.

[4]  Mark Potse,et al.  On Sampling Spatially-Correlated Random Fields for Complex Geometries , 2019, FIMH.

[5]  Paris Perdikaris,et al.  Multi-fidelity classification using Gaussian processes: accelerating the prediction of large-scale computational models , 2019, Computer Methods in Applied Mechanics and Engineering.

[6]  F. S. Costabal,et al.  Machine learning in drug development: Characterizing the effect of 30 drugs on the QT interval using Gaussian process regression, sensitivity analysis, and uncertainty quantification. , 2019, Computer methods in applied mechanics and engineering.

[7]  Paris Perdikaris,et al.  Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations , 2019, J. Comput. Phys..

[8]  Paris Perdikaris,et al.  Physics-Constrained Deep Learning for High-dimensional Surrogate Modeling and Uncertainty Quantification without Labeled Data , 2019, J. Comput. Phys..

[9]  Steven Niederer,et al.  A technique for measuring anisotropy in atrial conduction to estimate conduction velocity and atrial fibre direction , 2019, Comput. Biol. Medicine.

[10]  Wolfgang A. Wall,et al.  Automatic mapping of atrial fiber orientations for patient‐specific modeling of cardiac electromechanics using image registration , 2018, International journal for numerical methods in biomedical engineering.

[11]  Paris Perdikaris,et al.  Adversarial Uncertainty Quantification in Physics-Informed Neural Networks , 2018, J. Comput. Phys..

[12]  Marianna Meo,et al.  Universal atrial coordinates applied to visualisation, registration and construction of patient specific meshes , 2018, Medical Image Anal..

[13]  Yibo Yang,et al.  Physics-informed deep generative models , 2018, ArXiv.

[14]  Paris Perdikaris,et al.  Learning Parameters and Constitutive Relationships with Physics Informed Deep Neural Networks , 2018, 1808.03398.

[15]  Albin Cassirer,et al.  Randomized Prior Functions for Deep Reinforcement Learning , 2018, NeurIPS.

[16]  A. Quaglino,et al.  Fast uncertainty quantification of activation sequences in patient‐specific cardiac electrophysiology meeting clinical time constraints , 2018, International journal for numerical methods in biomedical engineering.

[17]  Lippincott Williams Wilkins,et al.  Correction to: Heart Disease and Stroke Statistics-2018 Update: A Report From the American Heart Association. , 2018, Circulation.

[18]  Henry Chubb,et al.  Local activation time sampling density for atrial tachycardia contact mapping: how much is enough? , 2017, Europace : European pacing, arrhythmias, and cardiac electrophysiology : journal of the working groups on cardiac pacing, arrhythmias, and cardiac cellular electrophysiology of the European Society of Cardiology.

[19]  Ellen Kuhl,et al.  Interpreting Activation Mapping of Atrial Fibrillation: A Hybrid Computational/Physiological Study , 2018, Annals of Biomedical Engineering.

[20]  Ellen Kuhl,et al.  The importance of mechano-electrical feedback and inertia in cardiac electromechanics. , 2017, Computer methods in applied mechanics and engineering.

[21]  Marie E. Rognes,et al.  Cbcbeat: an Adjoint-enabled Framework for Computational Cardiac Electrophysiology , 2017, J. Open Source Softw..

[22]  Mark Potse,et al.  Evaluation of a Rapid Anisotropic Model for ECG Simulation , 2017, Front. Physiol..

[23]  Houman Owhadi,et al.  Gaussian processes A hands-on tutorial , 2017 .

[24]  Ellen Kuhl,et al.  Generating Purkinje networks in the human heart. , 2016, Journal of biomechanics.

[25]  Yuan Yu,et al.  TensorFlow: A system for large-scale machine learning , 2016, OSDI.

[26]  Mathias Niepert,et al.  Learning Convolutional Neural Networks for Graphs , 2016, ICML.

[27]  Maria A. Zuluaga,et al.  Left Atrial Segmentation Challenge 2013: MRI training , 2015 .

[28]  Jimmy Ba,et al.  Adam: A Method for Stochastic Optimization , 2014, ICLR.

[29]  J. Wong,et al.  Generating fibre orientation maps in human heart models using Poisson interpolation , 2014, Computer methods in biomechanics and biomedical engineering.

[30]  G. V. Steeg Non-parametric Entropy Estimation Toolbox (NPEET) , 2014 .

[31]  Vincent Jacquemet,et al.  An eikonal-diffusion solver and its application to the interpolation and the simulation of reentrant cardiac activations , 2012, Comput. Methods Programs Biomed..

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

[33]  F. Ravelli,et al.  Automatic reconstruction of activation and velocity maps from electro-anatomic data by radial basis functions , 2010, 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology.

[34]  Yoshua Bengio,et al.  Understanding the difficulty of training deep feedforward neural networks , 2010, AISTATS.

[35]  Burr Settles,et al.  Active Learning Literature Survey , 2009 .

[36]  A. Kraskov,et al.  Estimating mutual information. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[37]  F. Fenton,et al.  Vortex dynamics in three-dimensional continuous myocardium with fiber rotation: Filament instability and fibrillation. , 1998, Chaos.

[38]  Dimitrios I. Fotiadis,et al.  Artificial neural networks for solving ordinary and partial differential equations , 1997, IEEE Trans. Neural Networks.

[39]  David A. Cohn,et al.  Active Learning with Statistical Models , 1996, NIPS.

[40]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

[41]  M. Stein Large sample properties of simulations using latin hypercube sampling , 1987 .

[42]  H. B. Mann,et al.  On a Test of Whether one of Two Random Variables is Stochastically Larger than the Other , 1947 .

[43]  F. Wilcoxon Individual Comparisons by Ranking Methods , 1945 .