Noninvasive Electroardiographic Imaging: Application of Hybrid Methods for Solving the Electrocardiography Inverse Problem

Computing the epicardial potentials from the body surface potentials constitutes one form of the ill-posed inverse problem of electrocardiography (ECG). In this paper, we employ hybrid methods combining the least square QR (LSQR) with truncated singular-value decomposition (TSVD) to solve the inverse problem of ECG. Hybrid methods are based on the Lanczos process, which yields a sequence of small bidiagonal systems approximating the original ill-posed problem, and on another additional direct regularization (the truncated SVD method is used in the present investigation), which is used to stabilize the iteration. The results show that determining of regularization parameters based on the final projected problem rather than on the original discretization one has firmer justification and it takes much less computational cost. The computation time could be reduced by several tenfolds typically, while the performance of the hybrid method is maintained well compared with TSVD, LSQR and GMRes methods. In addition, comparing with LSQR method, the hybrid method can obtain the inverse solutions without facing the "semi-convergence" problem.

[1]  P. Hansen Rank-Deficient and Discrete Ill-Posed Problems: Numerical Aspects of Linear Inversion , 1987 .

[2]  D. Calvetti,et al.  Noninvasive Electrocardiographic Imaging (ECGI): Application of the Generalized Minimal Residual (GMRes) Method , 2003, Annals of Biomedical Engineering.

[3]  Robert Modre,et al.  Lead field computation for the electrocardiographic inverse problem - finite elements versus boundary elements , 2005, Comput. Methods Programs Biomed..

[4]  Martin Hanke,et al.  On Lanczos Based Methods for the Regularization of Discrete Ill-Posed Problems , 2001 .

[5]  S Crozier,et al.  Electrodynamic Heart Model Construction and ECG Simulation , 2006, Methods of Information in Medicine.

[6]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[7]  R. Gulrajani,et al.  A new method for regularization parameter determination in the inverse problem of electrocardiography , 1997, IEEE Transactions on Biomedical Engineering.

[8]  L. Xia,et al.  Combination of the LSQR method and a genetic algorithm for solving the electrocardiography inverse problem. , 2007, Physics in medicine and biology.

[9]  D. O’Leary,et al.  A Bidiagonalization-Regularization Procedure for Large Scale Discretizations of Ill-Posed Problems , 1981 .

[10]  P. Favati,et al.  A framework for studying the regularizing properties of Krylov subspace methods , 2006 .

[11]  Y Rudy,et al.  The inverse problem in electrocardiography: solutions in terms of epicardial potentials. , 1988, Critical reviews in biomedical engineering.

[12]  Ling Xia,et al.  Simulation of Brugada syndrome using cellular and three-dimensional whole-heart modeling approaches , 2006, Physiological measurement.

[13]  Per Christian Hansen,et al.  Rank-Deficient and Discrete Ill-Posed Problems , 1996 .

[14]  G. Golub,et al.  Generalized cross-validation for large scale problems , 1997 .

[15]  Misha Elena Kilmer,et al.  Choosing Regularization Parameters in Iterative Methods for Ill-Posed Problems , 2000, SIAM J. Matrix Anal. Appl..

[16]  Ling Xia,et al.  Analysis of cardiac ventricular wall motion based on a three-dimensional electromechanical biventricular model , 2005, Physics in medicine and biology.