Multilevel optimization algorithm for Inverse Problem in Electrocardiography
- The electric conductivity of cardiac tissue determines excitation propagation and is vital for quantifying ischemia and scar tissue and building personalized models. As scar tissue is generally characterized by different conduction of electrical excitation, we aim to estimate conductivity-related parameters in mathematical excitation models from endocardial mapping data, particularly the anisotropic conductivity tensor in the monodomain equation, which describes the cardiac excitation. Yet, estimating the distributed and anisotropic conductivity tensors reliably and efficiently from endocardial mapping data or electrocardiograms is a challenging inverse problem due to the computational complexity of the monodomain equation; Many expensive high-resolution computations for the monodomain equation on very fine space and time discretizations are involved. Thus, we aim at building an efficient multilevel method for accelerating the estimation procedure combining electrophysiology models of different complex- ity, which uses a computationally cheap eikonal model in addition to the more accurate monodomain model. Distributed parameter estimation, well-known as an ill-posed inverse problem, can be performed by minimizing the misfit between simulated and measured electrical activity on the endocardial surface subject to the monodomain model and some regularization, leading to a partial differential equation constrained optimization problem. We formulate this optimization problem, including scar tissue modeling and different regularizations, and design an efficient iterative solver. To this aim, we consider monodomain grid hi- erarchies, monodomain-eikonal model hierarchies, and the combination of both hierarchies in a recursive multilevel trust-region (RMTR) method. On the one hand, both the trust region method’s estimation quality and efficiency, independent of the data, are investigated from several numerical exam- ples. Endocardial mapping data of realistic density appears to be sufficient to provide quantitatively reasonable estimates of the location, size, and shape of scars close to the endocardial surface. In several situations, scar reconstruction based on eikonal and monodomain models differ significantly, suggesting the use of the more involved monodomain model for this purpose. Moreover, Eikonal models can accelerate the computations considerably, enabling the use of complex electrophysiology models for estimating myocardial scars from endocardial mapping data. In many situations, eikonal models approximate monodomain models well but are orders of magnitude faster to solve. Thus, eikonal models can utilize them to provide an RMTR acceleration with negligible overhead per iteration, resulting in a practical approach to estimating myocardial scars from endocardial mapping data. In addition, the multilevel solver is faster than a comparable single-level solver. On the other hand, we investigate different optimization approaches based on adjoint gradient computation for computing a maximum posterior estimate: steepest descent, limited memory BFGS, and recursive multilevel trust region methods using mesh hierarchies or heterogeneous model hierarchies. We compare overall performance, asymptotic convergence rate, and pre-asymptotic progress on selected examples in order to assess the benefit of our multifidelity acceleration.
Author: | Fatemeh Chegini |
---|---|
Document Type: | Doctoral Thesis |
Granting Institution: | Università della Svizzera italiana |
Advisor: | Rolf Krause, Martin Weiser |
Date of final exam: | 2022/03/17 |
Year of first publication: | 2022 |