Premium
Fast uncertainty quantification of activation sequences in patient‐specific cardiac electrophysiology meeting clinical time constraints
Author(s) -
Quaglino A.,
Pezzuto S.,
Koutsourelakis P. S.,
Auricchio A.,
Krause R.
Publication year - 2018
Publication title -
international journal for numerical methods in biomedical engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.741
H-Index - 63
eISSN - 2040-7947
pISSN - 2040-7939
DOI - 10.1002/cnm.2985
Subject(s) - solver , computer science , uncertainty quantification , graphics , graphics processing unit , bayesian optimization , eikonal equation , field (mathematics) , algorithm , artificial intelligence , machine learning , mathematics , computer graphics (images) , programming language , operating system , mathematical analysis , pure mathematics
We present a fast, patient‐specific methodology for uncertainty quantification in electrophysiology, aimed at meeting the time constraints of clinical practitioners. We focus on computing the statistics of the activation map, given the uncertainties associated with the conductivity tensor modeling the fiber orientation in the heart. We use a fast parallel solution method implemented on a graphics processing unit for the eikonal approximation, in order to compute the activation map and to sample the random fiber field with correlation on the basis of geodesic distances. While this enables to perform uncertainty quantification studies with a manageable computational effort, the required time frame still exceeds clinically suitable time expectations. In order to reduce it further by 2 orders of magnitude, we rely on Bayesian multifidelity methods. In particular, we propose a low‐fidelity model that is patient‐specific and free from the additional training cost associated with reduced models. This is achieved by a sound physics–based simplification of the full eikonal model. The low‐fidelity output is then corrected by the standard multifidelity framework.[1][Koutsourelakis PS, 2009] In practice, the complete procedure only requires approximately 100 new runs of our eikonal graphics processing unit solver for producing the sought estimates and their associated credible intervals, enabling a full online analysis in less than 5 minutes.