An iterative linearized optimization technique for non‐linear ill‐posed problems applied to cardiac activation time imaging

Abstract

A promising approach for the solution of the electrocardiographic inverse problem is the calculation of the cardiac activation sequence from body surface potential (BSP) mapping data. Here, a two‐fold regularization scheme is applied in order to stabilize the inverse solution of this intrinsically ill‐posed problem. The solution of the inverse problem is defined by the minimum of a non‐linear cost function. The L‐curve method can be applied for regularization parameter determination. Solving the optimization problem by a Newton‐like method, the L‐curve may be of pronged shape. Then a numerically unique determination of the optimal regularization parameter will become difficult. This problem can be avoided applying an iterative linearized algorithm. It is shown that activation time imaging due to temporal and spatial regularization is stable with respect to large model errors. Even neglecting cardiac anisotropy in activation time imaging results in an acceptable inverse solution.