Search

The body surface electrocardiogram (ECG) is a direct result of electrical activity generated by the myocardium. Using the body surface ECGs to reconstruct cardiac electrical activity is called the inverse problem of electrocardiography. The method to solve the inverse problem depends on the chosen cardiac source model to describe cardiac electrical activity. In this paper, we describe the theoretical basis of two inverse methods based on the most commonly used cardiac source models: the epicardial potential model and the equivalent dipole layer model. We discuss similarities and differences in applicability, strengths and weaknesses and sketch a road towards improved inverse solutions by targeted use, sequential application or a combination of the two methods.

Introduction: Premature ventricular contractions (PVCs) are one of the most commonly targeted pathologies for ECGI validation, often through ventricular stimulation to mimic the ectopic beat. However, it remains unclear if such stimulated beats faithfully reproduce spontaneously occurring PVCs, particularly in the case of the R-on-T phenomenon. The objective of this study was to determine the differences in ECGI accuracy when reconstructing spontaneous PVCs as compared to ventricular-stimulated beats and to explore the impact of pathophysiological perturbation on this reconstruction accuracy.Methods: Langendorff-perfused pig hearts (n = 3) were suspended in a human torso-shaped tank, and local hyperkalemia was induced through perfusion of a high-K+ solution (8 mM) into the LAD. Recordings were taken simultaneously from the heart and tank surfaces during ventricular pacing and during spontaneous PVCs (including R-on-T), both at baseline and high K+. Epicardial potentials were reconstructed from torso potentials using ECGI.Results: Spontaneously occurring PVCs were better reconstructed than stimulated beats at baseline in terms of electrogram morphology [correlation coefficient (CC) = 0.74 ± 0.05 vs. CC = 0.60 ± 0.10], potential maps (CC = 0.61 ± 0.06 vs. CC = 0.51 ± 0.12), and activation time maps (CC = 0.86 ± 0.07 vs. 0.76 ± 0.10), though there was no difference in the localization error (LE) of epicardial origin (LE = 14 ± 6 vs. 15 ± 11 mm). High K+ perfusion reduced the accuracy of ECGI reconstructions in terms of electrogram morphology (CC = 0.68 ± 0.10) and AT maps (CC = 0.70 ± 0.12 and 0.59 ± 0.23) for isolated PVCs and paced beats, respectively. LE trended worse, but the change was not significant (LE = 17 ± 9 and 20 ± 12 mm). Spontaneous PVCs were less well when the R-on-T phenomenon occurred and the activation wavefronts encountered a line of block.Conclusion: This study demonstrates the differences in ECGI accuracy between spontaneous PVCs and ventricular-paced beats. We also observed a reduction in this accuracy near regions of electrically inactive tissue. These results highlight the need for more physiologically realistic experimental models when evaluating the accuracy of ECGI methods. In particular, reconstruction accuracy needs to be further evaluated in the presence of R-on-T or isolated PVCs, particularly when encountering obstacles (functional or anatomical) which cause line of block and re-entry.

Background: The optimal way to determine repolarization time (RT) from the intracardiac unipolar electrogram (UEG) has been a topic of debate for decades. RT is typically determined by either the Wyatt method or the “alternative method,” which both consider UEG T-wave slope, but differently.Objective: To determine the optimal method to measure RT on the UEG.Methods: Seven pig hearts surrounded by an epicardial sock with 100 electrodes were Langendorff-perfused with selective cannulation of the left anterior descending (LAD) coronary artery and submersed in a torso-shaped tank containing 256 electrodes on the torso surface. Repolarization was prolonged in the non-LAD-regions by infusing dofetilide and shortened in the LAD-region using pinacidil. RT was determined by the Wyatt (tWyatt) and alternative (tAlt) methods, in both invasive (recorded with epicardial electrodes) and in non-invasive UEGs (reconstructed with electrocardiographic imaging). tWyatt and tAlt were compared to local effective refractory period (ERP).Results: With contact mapping, mean absolute error (MAE) of tWyatt and tAlt vs. ERP were 21 ms and 71 ms, respectively. Positive T-waves typically had an earlier ERP than negative T-waves, in line with theory. tWyatt -but not tAlt-shortened by local infusion of pinacidil. Similar results were found for the non-invasive UEGs (MAE of tWyatt and tAlt vs. ERP were 30 ms and 92 ms, respectively).Conclusion: The Wyatt method is the most accurate to determine RT from (non) invasive UEGs, based on novel and historical analyses. Using it to determine RT could unify and facilitate repolarization assessment and amplify its role in cardiac electrophysiology.

Pulmonary vein isolation (PVI) is the most successful treatment for atrial fibrillation (AF) nowadays. However, not all AF patients benefit from PVI. In this study, we evaluate the use of ECGI to identify reentries and relate rotor density in the pulmonary vein (PV) area as an indicator of PVI outcome. Rotor maps were computed in a set of 29 AF patients using a new rotor detection algorithm. The relationship between the distribution of reentrant activity and the clinical outcome after PVI was studied. The number of rotors and proportion of PSs in different atrial regions were computed and compared retrospectively in two groups of patients: patients that remained in sinus rhythm 6 months after PVI and patients with arrhythmia recurrence. The total number of rotors obtained was higher in patients returning to arrhythmia after the ablation (4.31 ± 2.77 vs. 3.58 ± 2.67%, p = 0.018). However, a significantly higher concentration of PSs in the pulmonary veins was found in patients that remained in sinus rhythm (10.20 ± 12.40% vs. 5.19 ± 9.13%, p = 0.011) 6 months after PVI. The results obtained show a direct relationship between the expected AF mechanism and the electrophysiological parameters provided by ECGI, suggesting that this technology offers relevant information to predict the clinical outcome after PVI in AF patients.

Cardiac disease is a leading cause of morbidity and mortality in developed countries. Currently, non-invasive techniques that can identify patients at risk and provide accurate diagnosis and ablation guidance therapy are under development. One of these is electrocardiographic imaging (ECGI). In ECGI, the first step is to formulate a forward problem that relates the unknown potential sources on the cardiac surface to the measured body surface potentials. Then, the unknown potential sources on the cardiac surface are reconstructed through the solution of an inverse problem. Unfortunately, ECGI still lacks accuracy due to the underlying inverse problem being ill-posed, and this consequently imposes limitations on the understanding and treatment of many cardiac diseases. Therefore, it is necessary to improve the solution of the inverse problem. In this work, we transfer and adapt four inverse problem methods to the ECGI setting: algebraic reconstruction technique (ART), random ART, ART Split Bregman (ART-SB) and range restricted generalized minimal residual (RRGMRES) method. We test all these methods with data from the Experimental Data and Geometric Analysis Repository (EDGAR) and compare their solution with the recorded epicardial potentials provided by EDGAR and a generalized minimal residual (GMRES) iterative method computed solution. Activation maps are also computed and compared. The results show that ART achieved the most stable solutions and, for some datasets, returned the best reconstruction. Differences between the solutions derived from ART and random ART are almost negligible, and the accuracy of their solutions is followed by RRGMRES, ART-SB and finally the GMRES (which returned the worst reconstructions). The RRGMRES method provided the best reconstruction for some datasets but appeared to be less stable than ART when comparing different datasets. In conclusion, we show that the proposed methods (ART, random ART, and RRGMRES) improve the GMRES solution, which has been suggested as inverse problem solution for ECGI.

Background: Non-invasive cardiac mapping—also known as Electrocardiographic imaging (ECGi)—is a novel, painless and relatively economic method to map the electrical activation and repolarization patterns of the heart, providing a valuable tool for early identification and diagnosis of conduction abnormalities and arrhythmias. Moreover, the ability to obtain information on cardiac electrical activity non-invasively using ECGi provides the potential for a priori information to guide invasive surgical procedures, improving success rates, and reducing procedure time.Previous studies have shown the influence of clinical variables, such as heart rate, heart size, endocardial wall, and body composition on surface electrocardiogram (ECG) measurements. The influence of clinical variables on the ECG variability has provided information on cardiovascular control and its abnormalities in various pathologies. However, the effects of such clinical variables on the Body Surface Potential (BSP) and ECGi maps have yet to be systematically investigated.Methods: In this study we investigated the effects of heart size, intracardiac thickness, and heart rate on BSP and ECGi maps using a previously-developed 3D electrophysiologically-detailed ventricles-torso model. The inverse solution was solved using the three different Tikhonov regularization methods.Results: Through comparison of multiple measures of error/accuracy on the ECGi reconstructions, our results showed that using different heart geometries to solve the forward and inverse problems produced a larger estimated focal excitation location. An increase of ~2 mm in the Euclidean distance error was observed for an increase in the heart size. However, the estimation of the location of focal activity was still able to be obtained. Similarly, a Euclidean distance increase was observed when the order of regularization was reduced.For the case of activation maps reconstructed at the same ectopic focus location but different heart rates, an increase in the errors and Euclidean distance was observed when the heart rate was increased.Conclusions: Non-invasive cardiac mapping can still provide useful information about cardiac activation patterns for the cases when a different geometry is used for the inverse problem compared to the one used for the forward solution; rapid pacing rates can induce order-dependent errors in the accuracy of reconstruction.

Background: Non-invasive electrocardiographic imaging (ECGI) is a promising tool to provide high-resolution panoramic imaging of cardiac electrical activity noninvasively from body surface potential measurements. Current experimental methods for ECGI validation are limited to comparison with unipolar electrograms and the relatively low spatial resolution of cardiac mapping arrays. We aim to develop a novel experimental set up combining a human shaped torso tank with high-resolution optical mapping allowing the validation of ECGI reconstructions.Methods: Langendorff-perfused pig hearts (n = 3) were suspended in a human torso-shaped tank, with the left anterior descending artery (LAD) cannulated on a separate perfusion. Electrical signals were recorded from an 108-electrode epicardial sock and 128 electrodes embedded in the tank surface. Simultaneously, optical mapping of the heart was performed through the anterior surface of the tank. Recordings were made in sinus rhythm and ventricular pacing (n = 55), with activation and repolarization heterogeneities induced by perfusion of hot and cold solutions as well as Sotalol through the LAD. Fluoroscopy provided 3D cardiac and electrode geometries in the tank that were transformed to the 2D optical mapping window using an optimization algorithm. Epicardial unipolar electrograms were reconstructed from torso potentials using ECGI and validated using optical activation and repolarization maps.Results: The transformation and alignment of the 3D geometries onto the 2D optical mapping window was good with an average correlation of 0.87 ± 0.10 and error of 7.7 ± 3.1 ms with activation derived from the sock. The difference in repolarization times were more substantial (error = 17.4 ± 3.7 ms) although the sock and optical repolarization patterns themselves were very similar (correlation = 0.83 ± 0.13). Validation of ECGI reconstructions revealed ECGI accurately captures the pattern of activation (correlation = 0.79 ± 0.11) and identified regions of late and/or early repolarization during different perfusions through LAD. ECGI also correctly demonstrated gradients in both activation and repolarization, although in some cases these were under or over-estimated or shifted slightly in space.Conclusion: A novel experimental setup has been developed, combining a human-shaped torso tank with optical mapping, which can be effectively used in the validation of ECGI techniques; including the reconstruction of activation and repolarization patterns and gradients.

Electrocardiography provides some information useful for ischemic diagnosis. However, more recently there has been substantial growth in the area of ECG imaging, which by solving the inverse problem of electrocardiography aims to produce high-resolution mapping of the electrical and magnetic dynamics of the heart. Most inverse studies use the full resolution of the body surface potential (BSP) to reconstruct the epicardial potentials, however using a limited number of torso electrodes to interpolate the BSP is more clinically relevant and has an important effect on the reconstruction which must be quantified. A circular ischemic lesion on the right ventricle lateral wall 27 mm in radius is reconstructed using three Tikhonov methods along with 6 different electrode configurations ranging from 32 leads to 1,024 leads. The 2nd order Tikhonov solution performed the most accurately (~80% lesion identified) followed by the 1st (~50% lesion identified) and then the 0 order Tikhonov solution performed the worst with a maximum of ~30% lesion identified regardless of how many leads were used. With an increasing number of leads the solution produces less error, and the error becomes more localised around the lesion for all three regularisation methods. In noisy conditions, the relative performance gap of the 1st and 2nd order Tikhonov solutions was reduced, and determining an accurate regularisation parameter became relatively more difficult. Lesions located on the left ventricle walls were also able to be identified but comparatively to the right ventricle lateral wall performed marginally worse with lesions located on the interventricular septum being able to be indicated by the reconstructions but not successfully identified against the error. The quality of reconstruction was found to decrease as the lesion radius decreased, with a lesion radius of 512 torso electrodes.

Electrocardiographic imaging aims at reconstructing cardiac electrical events from electrical signals measured on the body surface. The most common approach relies on the inverse solution of the Laplace equation in the torso to reconstruct epicardial potential maps from body surface potential maps. Here we apply a method based on a parameter identification problem to reconstruct both activation and repolarization times. From an ansatz of action potential, based on the Mitchell-Schaeffer ionic model, we compute body surface potential signals. The inverse problem is reduced to the identification of the parameters of the Mitchell-Schaeffer model. We investigate whether solving the inverse problem with the endocardium improves the results or not. We solved the parameter identification problem on two different meshes: one with only the epicardium, and one with both the epicardium and the endocardium. We compared the results on both the heart (activation and repolarization times) and the torso. The comparison was done on validation data of sinus rhythm and ventricular pacing. We found similar results with both meshes in 6 cases out of 7: the presence of the endocardium slightly improved the activation times. This was the most visible on a sinus beat, leading to the conclusion that inclusion of the endocardium would be useful in situations where endo-epicardial gradients in activation or repolarization times play an important role.

The accurate generation of forward models is an important element in general research in electrocardiography, and in particular for the techniques for ElectroCardioGraphic Imaging (ECGI). Recent research efforts have been devoted to the reliable and fast generation of forward models. However, these model can suffer from several sources of inaccuracy, which in turn can lead to considerable error in both the forward simulation of body surface potentials and even more so for ECGI solutions. In particular, the accurate localization of the heart within the torso is sensitive to movements due to respiration and changes in position of the subject, a problem that cannot be resolved with better imaging and segmentation alone. Here, we propose an algorithm to localize the position of the heart using electrocardiographic recordings on both the heart and torso surface over a sequence of cardiac cycles. We leverage the dependency of electrocardiographic forward models on the underlying geometry to parameterize the forward model with respect to the position (translation) and orientation of the heart, and then estimate these parameters from heart and body surface potentials in a numerical inverse problem. We show that this approach is capable of localizing the position of the heart in synthetic experiments and that it reduces the modeling error in the forward models and resulting inverse solutions in canine experiments. Our results show a consistent decrease in error of both simulated body surface potentials and inverse reconstructed heart surface potentials after re-localizing the heart based on our estimated geometric correction. These results suggest that this method is capable of improving electrocardiographic models used in research settings and suggest the basis for the extension of the model presented here to its application in a purely inverse setting, where the heart potentials are unknown.