A new spatiotemporal regularization approach for reconstruction of cardiac transmembrane potential patterns

The single beat reconstruction of electrical cardiac sources from body-surface electrocardiogram data might become an important issue for clinical application. The feasibility and field of application of noninvasive imaging methods strongly depend on development of stable algorithms for solving the underlying ill-posed inverse problems. We propose a novel spatiotemporal regularization approach for the reconstruction of surface transmembrane potential (TMP) patterns. Regularization is achieved by imposing linearly formulated constraints on the solution in the spatial as well as in the temporal domain. In the spatial domain an operator similar to the surface Laplacian, weighted by a regularization parameter, is used. In the temporal domain monotonic nondecreasing behavior of the potential is presumed. This is formulated as side condition without the need of any regularization parameter. Compared to presuming template functions, the weaker temporal constraint widens the field of application because it enables the reconstruction of TMP patterns with ischemic and infarcted regions. Following the line of Tikhonov regularization, but considering all time points simultaneously, we obtain a linearly constrained sparse large-scale convex optimization problem solved by a fast interior point optimizer. We demonstrate the performance with simulations by comparing reconstructed TMP patterns with the underlying reference patterns.