Computationally Efficient Noninvasive Cardiac Activation Time Imaging

Summary

Objective: 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.

Methods: 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.

Results: 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.

Conclusions: As it was shown for data from a WPW-patient 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.

Keywords

Cardiac electrophysiology - computer modeling - inverse problem

• References

• 1 Schilling RJ. Which patient should be referred to an electrocardiologist: supraventricular tachycardia. Heart 2002; 36 (03) 299-304.
• 2 Cheng LK, Bodley JM, Pullan AJ. Comparison of potential- and activation-based formulations for the inverse problem of electrocardiology. IEEE Trans Biomed Eng 2003; 50 (01) 11-22.
• 3 Messnarz B, Tilg B, Modre R, Fischer G, Hanser F. A new spatiotemporal regularization approach for reconstruction of cardiac transmembrane potential patterns. IEEE Trans Biomed Eng 2004; 51 (02) 273-81.
• 4 Wu D, Ono K, Hosaka H, He B. Simulation of body surface Laplacian maps during ventricular pacing in a 3D inhomogeneous heart-torso model. Methods Inf Med 2000; 39 (02) 196-9.
• 5 Modre R, Tilg B, Fischer G, Hanser F, Messnarz B, Schocke M. et al Atrial noninvasive activation mapping of paced rhythm data. J Cardiovasc Electrophysiol 2003; 14 (07) 712-9.
• 6 Tilg B, Fischer G, Modre R, Hanser F, Messnarz B, Schocke M. et al Model-based imaging of cardiac electrical excitation in humans. IEEE Trans Med Imag 2002; 21 (09) 1031-9.
• 7 Ramanathan C, Ghanem RN, Jia P, Ryu K, Rudy Y. Noninvasive electrocardiographic imaging for cardiac electrophysiology and arrhythmia. Nat Med 2004; 10 (04) 422-8.
• 8 Greensite Fred. Cardiac electromagnetic imaging as an inverse problem. Electromagnetics 2001; 21 (07) 555-77.
• 9 Huiskamp G, van Oosterom A. The depolarization sequence of the human heart surface computed from measured body surface potentials. IEEE Trans Biomed Eng 1988; 35 (12) 1047-58.
• 10 Oostendorp T, Nenonen J, Korhonen P. Noninvasive determination of the activation sequence of the heart: application to patients with previous myocardial infarctions. J Electrocardiol 2002; 35 Supplement 75-80.
• 11 Pullan AJ, Cheng LK, Nash MP, Bradley CP, Paterson DJ. Noninvasive electrical imaging of the heart: theory and model development. Ann Biomed Eng 2001; 29 (10) 817-36.
• 12 Malmivuo J, Plonsey R. Bioelectromagnetism. Oxford University Press; 1995
• 13 Fischer G, Tilg B, Modre R, Hanser F, Messnarz B, Wach P. On modeling the Wilson terminal in the boundary and finite element method. IEEE Trans Biomed Eng 2002; 49 (03) 217-24.
• 14 Hoekema R, Uijen G, van Oosterom A. The number of independent signals in body surface maps. Methods Inf Med 1999; 38 (02) 119-24.
• 15 Huiskamp GJ. Difference formulas for the surface Laplacian on a triangulated surface. Journal of Computational Physics 1991; 91 (02) 477-96.
• 16 Fischer G, Tilg B, Wach P, Modre R, Leder U, Nowak H. Application of high-order boundary elements to the electrocardiographic inverse problem. Comput Methods Programs Biomed 1999; 58 (02) 119-31.
• 17 Press WH, Teukolsky SA, Vetterling WT, Flannery BP. Numerical Recipes in C++ – The Art of Scientific Computing. Oxford University Press; 2002
• 18 Fischer G, Seger M, Hanser F, Modre R, Tilg B. A Galerkin boundary element formulation for noninvasive cardiac activation time imaging – the forward problem. Ann Biomed Eng; submitted
• 19 Oppenheim AV, Willsky AS. Signals and Systems. Prentice Hall 1996
• 20 Weber H. Laplace-Transformation für Ingenieure der Elektrotechnik. Teubner 1987
• 21 Bailey JJ, Berson AS, Garson Jr A, Horan LG, Macfarlane PW, Mortara DW. et al Recommendations for standardization and specifications in automated electrocardiography: bandwidth and digital signal processing. A report for health professionals by an ad hoc writing group of the Committee on Electrocardiography and Cardiac Electrophysiology of the Council on Clinical Cardiology, American Heart Association. Circulation 1990; 81 (02) 730-9.
• 22 Huiskamp G, Greensite F. A new method for myocardial activation imaging. IEEE Trans Biomed Eng 1997; 44 (06) 433-46.
• 23 Modre RM, SegerFischer G, Hintermüller ChHayn D, Pfeifer B. et al Cardiac anisotropy: is it negligible for noninvasive activation time imaging. IEEE Trans Biomed Eng; in press
• 24 Modre R, Tilg B, Fischer G, Hanser F, Messnarz B, Seger M. et al Ventricular surface activation time imaging from electrocardiogram mapping data. Med Biol Eng Comput 2004; 42 (02) 146-50.
• 25 Ramlau R. TIGRA – an iterative algorithm for regularizing nonlinear ill-posed problems. Inverse Problems 2003; 19 (02) 433-65.
• 26 Seger M, Fischer G, Modre R, Messnarz B, Hanser F, Tilg B. Lead field computation for the electrocardiographic inverse problem – finite elements versus boundary elements. Comput Methods Programs Biomed 2005; 77 (03) 241-52.