Inverse electrocardiographic source localization of ischemia: An optimization framework and finite element solution

Authors:
Dafang Wang;Robert M. Kirby;Rob S. Macleod;Chris R. Johnson
Affiliations:
-;-;-;-
Venue:
Journal of Computational Physics
Year:
2013

Citing 7
Cited 0

Iterative methods for total variation denoising

SIAM Journal on Scientific Computing - Special issue on iterative methods in numerical linear algebra; selected papers from the Colorado conference
A Nonlinear Primal-Dual Method for Total Variation-Based Image Restoration

SIAM Journal on Scientific Computing
Convex Optimization

Convex Optimization
Parallel Lagrange--Newton--Krylov--Schur Methods for PDE-Constrained Optimization. Part I: The Krylov--Schur Solver

SIAM Journal on Scientific Computing
On the use of the resting potential and level set methods for identifying ischemic heart disease: An inverse problem

Journal of Computational Physics
Statistical inverse problems: discretization, model reduction and inverse crimes

Journal of Computational and Applied Mathematics - Special issue: Applied computational inverse problems
Duality-based algorithms for total-variation-regularized image restoration

Computational Optimization and Applications

Quantified Score

Hi-index	31.45

Visualization

Abstract

With the goal of non-invasively localizing cardiac ischemic disease using body-surface potential recordings, we attempted to reconstruct the transmembrane potential (TMP) throughout the myocardium with the bidomain heart model. The task is an inverse source problem governed by partial differential equations (PDE). Our main contribution is solving the inverse problem within a PDE-constrained optimization framework that enables various physically-based constraints in both equality and inequality forms. We formulated the optimality conditions rigorously in the continuum before deriving finite element discretization, thereby making the optimization independent of discretization choice. Such a formulation was derived for the L"2-norm Tikhonov regularization and the total variation minimization. The subsequent numerical optimization was fulfilled by a primal-dual interior-point method tailored to our problem's specific structure. Our simulations used realistic, fiber-included heart models consisting of up to 18,000 nodes, much finer than any inverse models previously reported. With synthetic ischemia data we localized ischemic regions with roughly a 10% false-negative rate or a 20% false-positive rate under conditions up to 5% input noise. With ischemia data measured from animal experiments, we reconstructed TMPs with roughly 0.9 correlation with the ground truth. While precisely estimating the TMP in general cases remains an open problem, our study shows the feasibility of reconstructing TMP during the ST interval as a means of ischemia localization.