Application of robust Generalised Cross-Validation to the inverse problem of electrocardiology

Robust Generalised Cross-Validation was proposed recently as a method for determining near optimal regularisation parameters in inverse problems. It was introduced to overcome a problem with the regular Generalised Cross-Validation method in which the function that is minimised to obtain the regularisation parameter often has a broad, flat minimum, resulting in a poor estimate for the parameter. The robust method defines a new function to be minimised which has a narrower minimum, but at the expense of introducing a new parameter called the robustness parameter. In this study, the Robust Generalised Cross-Validation method is applied to the inverse problem of electrocardiology. It is demonstrated that, for realistic situations, the robustness parameter can be set to zero. With this choice of robustness parameter, it is shown that the robust method is able to obtain estimates of the regularisation parameter in the inverse problem of electrocardiology that are comparable to, or better than, many of the standard methods that are applied to this inverse problem.

[1]  Leo K. Cheng,et al.  The inverse problem of electrocardiography , 2011 .

[2]  Mark A. Lukas,et al.  Comparing parameter choice methods for regularization of ill-posed problems , 2011, Math. Comput. Simul..

[3]  F. J. Claydon,et al.  A volume conductor model of the thorax for the study of defibrillation fields , 1988, IEEE Transactions on Biomedical Engineering.

[4]  Per Christian Hansen,et al.  Truncated Singular Value Decomposition Solutions to Discrete Ill-Posed Problems with Ill-Determined Numerical Rank , 1990, SIAM J. Sci. Comput..

[5]  R. Macleod,et al.  The Forward Problem of Electrocardiography , 2010 .

[6]  A. N. Tikhonov,et al.  Solutions of ill-posed problems , 1977 .

[7]  Robert Michael Kirby,et al.  Inverse electrocardiographic source localization of ischemia: An optimization framework and finite element solution , 2013, J. Comput. Phys..

[8]  Olaf Dössel,et al.  Use of a Tolerance Operator in Wavefront-Based ECG Imaging of Transmembrane Voltages , 2013 .

[9]  Robert Michael Kirby,et al.  Resolution Strategies for the Finite-Element-Based Solution of the ECG Inverse Problem , 2010, IEEE Transactions on Biomedical Engineering.

[10]  Alireza Mazloumi Gavgani,et al.  Noise reduction using anisotropic diffusion filter in inverse electrocardiology , 2012, 2012 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[11]  G. Huiskamp,et al.  An improved method for estimating epicardial potentials from the body surface , 1998, IEEE Transactions on Biomedical Engineering.

[12]  C. A. Brebbia.  Boundary elements ten, an introduction , 1988 .

[13]  R. Gulrajani.  Bioelectricity and biomagnetism , 1998 .

[14]  William H. Press,et al.  Numerical recipes , 1990 .

[15]  Johnston,et al.  au An Analysis of the Zero-Crossing Method for Choosing Regularisation Parameters , 2017 .

[16]  Alireza Mazloumi Gavgani,et al.  Genetic algorithm-based regularization parameter estimation for the inverse electrocardiography problem using multiple constraints , 2012, Medical & Biological Engineering & Computing.

[17]  Andrew J. Pullan,et al.  Effects of experimental and modeling errors on electrocardiographic inverse formulations , 2003, IEEE Transactions on Biomedical Engineering.

[18]  M. A. Lukas.  Robust generalized cross-validation for choosing the regularization parameter , 2006 .

[19]  Y. Rudy,et al.  Application of the Method of Fundamental Solutions to Potential-based Inverse Electrocardiography , 2006, Annals of Biomedical Engineering.

[20]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[21]  Linwei Wang.  Computational reduction for noninvasive transmural electrophysiological imaging , 2013, Comput. Biol. Medicine.

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

[23]  G. Huiskamp,et al.  A new method for myocardial activation imaging , 1997, IEEE Transactions on Biomedical Engineering.

[24]  Y. Rudy,et al.  The use of temporal information in the regularization of the inverse problem of electrocardiography , 1990, IEEE Transactions on Biomedical Engineering.

[25]  G. Wahba.  Practical Approximate Solutions to Linear Operator Equations When the Data are Noisy , 1977 .

[26]  Quan Ni,et al.  Three-Dimensional Activation Mapping in Ventricular Muscle: Interpolation and Approximation of Activation Times , 1999, Annals of Biomedical Engineering.

[27]  R M Gulrajani,et al.  A computer heart model incorporating anisotropic propagation. III. Simulation of ectopic beats. , 1993, Journal of electrocardiology.

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

[29]  B. M. Horacek,et al.  Numerical Model of an Inhomogeneous Human Torso1 , 1974 .

[30]  Carlos Alberto Brebbia,et al.  Boundary Elements: An Introductory Course , 1989 .

[31]  Alireza Mazloumi Gavgani,et al.  Use of genetic algorithm for selection of regularization parameters in multiple constraint inverse ECG problem , 2011, 2011 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[32]  Per Christian Hansen,et al.  Analysis of Discrete Ill-Posed Problems by Means of the L-Curve , 1992, SIAM Rev..

[33]  Rob S. MacLeod,et al.  Inverse electrocardiography by simultaneous imposition of multiple constraints , 1999, IEEE Transactions on Biomedical Engineering.

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

[35]  R M Gulrajani,et al.  A computer heart model incorporating anisotropic propagation. I. Model construction and simulation of normal activation. , 1993, Journal of electrocardiology.