Computationally efficient noninvasive cardiac activation time imaging

The computer model-based computation of the cardiac activation sequence in humans has been recently subject of successful clinical validation. This method is of potential interest for guiding ablation therapy of arrhythmogenic substrates. However, computation times of almost an hour are unattractive in a clinical setting. Thus, the objective is the development of a method which performs the computation in a few minutes run time.

The computationally most expensive part is the product of the lead field matrix with a matrix containing the source pattern on the cardiac surface. The particular biophysical properties of both matrices are used for speeding up this operation by more than an order of magnitude. A conjugate gradient optimizer was developed using C++ for computing the activation map.

The software was tested on synthetic and clinical data. The increase in speed with respect to the previously used Fortran 77 implementation was a factor of 30 at a comparable quality of the results. As an additional finding the coupled regularization strategy, originally introduced for saving computation time, also reduced the sensitivity of the method to the choice of the regularization parameter.

As it was shown for data from a WPWpatient the developed software can deliver diagnostically valuable information at a much shorter span of time than current clinical routine methods. Its main application could be the localization of focal arrhythmogenic substrates.