Sensitivity- and effort-gain analysis: multilead ECG electrode array selection for activation time imaging.

Methods for noninvasive imaging of electric function of the heart might become clinical standard procedure the next years. Thus, the overall procedure has to meet clinical requirements as an easy and fast application. In this paper, we propose a new electrode array which improves the resolution of methods for activation time imaging considering clinical constraints such as easy to apply and compatibility with routine leads. For identifying the body-surface regions where the body surface potential (BSP) is most sensitive to changes in transmembrane potential (TMP), a virtual array method was used to compute local linear dependency (LLD) maps. The virtual array method computes a measure for the LLD in every point on the body surface. The most suitable number and position of the electrodes within the sensitive body surface regions was selected by constructing effort gain (EG) plots. Such a plot depicts the relative attainable rank of the leadfield matrix in relation to the increase in number of electrodes required to build the electrode array. The attainable rank itself was computed by a detector criterion. Such a criterion estimates the maximum number of source space eigenvectors not covered by noise when being mapped to the electrode space by the leadfield matrix and recorded by a detector. From the sensitivity maps, we found that the BSP is most sensitive to changes in TMP on the upper left frontal and dorsal body surface. These sensitive regions are covered best by an electrode array consisting of two L-shaped parts of approximately 30 cm x 30 cm and approximately 20 cm x 20 cm. The EG analysis revealed that the array meeting clinical requirements best and improving the resolution of activation time imaging consists of 125 electrodes with a regular horizontal and vertical spacing of 2-3 cm.

Cardiac anisotropy: is it negligible regarding noninvasive activation time imaging?

The aim of this study was to quantify the effect of cardiac anisotropy in the activation-based inverse problem of electrocardiography. Differences of the patterns of simulated body surface potential maps for isotropic and anisotropic conditions were investigated with regard to activation time (AT) imaging of ventricular depolarization. AT maps were estimated by solving the nonlinear inverse ill-posed problem employing spatio-temporal regularization. Four different reference AT maps (sinus rhythm, right-ventricular and septal pacing, accessory pathway) were calculated with a bidomain theory based anisotropic finite-element heart model in combination with a cellular automaton. In this heart model a realistic fiber architecture and conduction system was implemented. Although the anisotropy has some effects on forward solutions, effects on inverse solutions are small indicating that cardiac anisotropy might be negligible for some clinical applications (e.g., imaging of focal events) of our AT imaging approach. The main characteristic events of the AT maps were estimated despite neglected electrical anisotropy in the inverse formulation. The worst correlation coefficient of the estimated AT maps was 0.810 in case of sinus rhythm. However, all characteristic events of the activation pattern were found. The results of this study confirm our clinical validation studies of noninvasive AT imaging in which cardiac anisotropy was neglected.

Model-based imaging of cardiac electrical excitation in humans.

Activation time (AT) imaging from electrocardiographic (ECG) mapping data has been developing for several years. By coupling ECG mapping and three-dimensional (3-D) + time anatomical data, the electrical excitation sequence can be imaged completely noninvasively in the human heart. In this paper, a bidomain theory-based surface heart model AT imaging approach was applied to single-beat data of atrial and ventricular depolarization in two patients with structurally normal hearts. In both patients, the AT map was reconstructed from sinus and paced rhythm data. Pacing sites were the apex of the right ventricle and the coronary sinus (CS) ostium. For CS pacing, the reconstructed AT pattern on the endocardium of the right atrium was compared with the CARTO map in both patients. The localization errors of the origins of the initial endocardial breakthroughs were determined to be 6 and 12 mm. The sites of early activation and the areas with late activation were estimated with sufficient accuracy. The reconstructed sinus rhythm sequence was in good qualitative agreement with the pattern previously published for the isolated Langendorff-perfused human heart.

Noninvasive myocardial activation time imaging: a novel inverse algorithm applied to clinical ECG mapping data.

Linear approaches like the minimum-norm least-square algorithm show insufficient performance when it comes to estimating the activation time map on the surface of the heart from electrocardiographic (ECG) mapping data. Additional regularization has to be considered leading to a nonlinear problem formulation. The Gauss-Newton approach is one of the standard mathematical tools capable of solving this kind of problem. To our experience, this algorithm has specific drawbacks which are caused by the applied regularization procedure. In particular, under clinical conditions the amount of regularization cannot be determined clearly. For this reason, we have developed an iterative algorithm solving this nonlinear problem by a sequence of regularized linear problems. At each step of iteration, an individual L-curve is computed. Subsequent iteration steps are performed with the individual optimal regularization parameter. This novel approach is compared with the standard Gauss-Newton approach. Both methods are applied to simulated ECG mapping data as well as to single beat sinus rhythm data from two patients recorded in the catheter laboratory. The proposed approach shows excellent numerical and computational performance, even under clinical conditions at which the Gauss-Newton approach begins to break down.

A new spatiotemporal regularization approach for reconstruction of cardiac transmembrane potential patterns.

The single-beat reconstruction of electrical cardiac sources from body-surface electrocardiogram data might become an important issue for clinical application. The feasibility and field of application of noninvasive imaging methods strongly depend on development of stable algorithms for solving the underlying ill-posed inverse problems. We propose a novel spatiotemporal regularization approach for the reconstruction of surface transmembrane potential (TMP) patterns. Regularization is achieved by imposing linearly formulated constraints on the solution in the spatial as well as in the temporal domain. In the spatial domain an operator similar to the surface Laplacian, weighted by a regularization parameter, is used. In the temporal domain monotonic nondecreasing behavior of the potential is presumed. This is formulated as side condition without the need of any regularization parameter. Compared to presuming template functions, the weaker temporal constraint widens the field of application because it enables the reconstruction of TMP patterns with ischemic and infarcted regions. Following the line of Tikhonov regularization, but considering all time points simultaneously, we obtain a linearly constrained sparse large-scale convex optimization problem solved by a fast interior point optimizer. We demonstrate the performance with simulations by comparing reconstructed TMP patterns with the underlying reference patterns.

A comparison of noninvasive reconstruction of epicardial versus transmembrane potentials in consideration of the null space.

We compare two source formulations for the electrocardiographic forward problem in consideration of their implications for regularizing the ill-posed inverse problem. The established epicardial potential source model is compared with a bidomain-theory-based transmembrane potential source formulation. The epicardial source approach is extended to the whole heart surface including the endocardial surfaces. We introduce the concept of the numerical null and signal space to draw attention to the problems associated with the nonuniqueness of the inverse solution and show that reconstruction of null-space components is an important issue for physiologically meaningful inverse solutions. Both formulations were tested with simulated data generated with an anisotropic heart model and with clinically measured data of two patients. A linear and a recently proposed quasi-linear inverse algorithm were applied for reconstructions of the epicardial and transmembrane potential, respectively. A direct comparison of both formulations was performed in terms of computed activation times. We found the transmembrane potential-based formulation is a more promising source formulation as stronger regularization by incorporation of biophysical a priori information is permitted.

On modeling the Wilson terminal in the boundary and finite element method.

In clinical electrocardiography, the zero-potential is commonly defined by the Wilson central terminal. In the electrocardiographic forward and inverse problem, the zero-potential is often defined in a different way, e.g., by the sum of all node potentials yielding zero. This study presents relatively simple to implement techniques, which enable the incorporation of the Wilson Terminal in the boundary element method (BEM) and finite element method (FEM). For the BEM, good results are obtained when properly adopting matrix deflation for modeling the Wilson terminal. Applying other zero-potential-definitions, the obtained solutions contained a remarkable offset with respect to the reference defined by the Wilson terminal. In the inverse problem (nonlinear dipole fit), errors introduced by an erroneous zero-potential-definition can lead to displacements of more than 5 mm in the computed dipole location. For the FEM, a method similar to matrix deflation is proposed in order to properly consider the Wilson central terminal. The matrix obtained from this manipulation is symmetric, sparse and positive definite enabling the application of standard FEM-solvers.

Electrocardiographic imaging of atrial and ventricular electrical activation.

Inverse electrocardiography has been developing for several years. By combining measurements obtained by electrocardiographic body surface mapping with three-dimensional anatomical data, one can non-invasively image the electrical activation sequence in the human heart. In this study, an imaging approach that uses a bidomain theory-based surface heart model was applied to single-beat data of atrial and ventricular activation. We found that for sinus and paced rhythms, the sites of early activation and the areas with late activation were estimated with sufficient accuracy. In particular, for focal arrhythmias, this model-based imaging approach might allow the guidance and evaluation of antiarrhythmic interventions, for instance, in case of catheter ablation or drug therapy.

Semiautomatic volume conductor modeling pipeline for imaging the cardiac electrophysiology noninvasively.

In this paper we present an approach for extracting patient individual volume conductor models (VCM) using volume data acquired from Magnetic Resonance Imaging (MRI) for computational biology of electrical excitation in the patient's heart. The VCM consists of the compartments chest surface, lung surfaces, the atrial and ventricular myocardium, and the blood masses. For each compartment a segmentation approach with no or little necessity of user interaction was implemented and integrated into a VCM segmentation pipeline to enable the inverse problem of electrocardiography to become clinical applicable. The segmentation pipeline was tested using volume data from ten patients with structurally normal hearts.

Single-beat noninvasive imaging of cardiac electrophysiology of ventricular pre-excitation.

OBJECTIVES: The aim of this study was to determine whether noninvasive imaging of cardiac electrophysiology (NICE) is feasible in patients with Wolff-Parkinson-White (WPW) syndrome in the clinical setting of a catheter laboratory and to test the accuracy of the noninvasively obtained ventricular activation sequences as compared with that of standard invasive electroanatomic mapping. BACKGROUND: NICE of ventricular activation could serve as a useful tool in the treatment of cardiac arrhythmias and might help improve our understanding of arrhythmia mechanisms. METHODS: NICE works by fusing the data from high-resolution electrocardiographic mapping and a model of the patient's cardiac anatomy obtained by magnetic resonance imaging. The ventricular activation sequence was computed with a bidomain theory-based heart model to solve this inverse problem. Noninvasive imaging of cardiac electrophysiology was performed in 7 patients with WPW syndrome undergoing catheter ablation of the accessory pathway. The position error of NICE was defined as the distance between the site of earliest activation computed by NICE and the successful ablation site identified by electroanatomic mapping (CARTO; Biosense Webster, Diamond Bar, California) for normal atrioventricular (AV) conduction as well as for adenosine-induced AV block. RESULTS: The error introduced by geometric coupling of the CARTO data and the NICE model was 5 +/- 3 mm (model discretization 10 mm). All ventricular accessory pathway insertion sites were identified with an accuracy of 18.7 +/- 5.8 mm (baseline) and 18.7 +/- 6.4 mm (adenosine). CONCLUSIONS: The individual cardiac anatomy model obtained for each patient enables accurate noninvasive electrocardiographic imaging of ventricular pre-excitation in patients with WPW syndrome. Noninvasive imaging of cardiac electrophysiology might be used as a complementary noninvasive approach to localize the origin and help identify and understand the underlying mechanisms of cardiac arrhythmias.

Rank Based Selection of Electrode Positions for a Multi-Lead ECG Electrode Array.

Methods for noninvasive imaging of electrical function of the heart seem to become a clinical standard procedure the next years. Thus, the overall procedure has to meet clinical requirements as easy and fast application. In this study we propose a new electrode array meeting clinical requirements such as easy to apply and compatibility with routine leads. Within body surface regions of high sensitivity, identified in a prior, information content based study, the number of required electrodes was optimized using effort-gain plots. These plots were generated by applying a so called type one detector criterion. The optimal array was selected from a set of 12 electrode arrays. Each of them consists of two L-shaped regular spaced parts. The optimal array was found by comparing several layouts and electrode densities to the electrode array we use for clinical studies. It consists of 125 electrodes with a regular spacing between 2cm and 3cm.

Effects of cardiac resynchronization therapy on ventricular repolarization in patients with congestive heart failure.

INTRODUCTION: Biventricular pacing has been shown to improve the clinical status of patients with congestive heart failure, but little is known about its influence on ventricular repolarization. The aim of our study was to evaluate the effect of biventricular pacing on ECG markers of ventricular repolarization in patients with congestive heart failure. METHODS AND RESULTS: Twenty-five patients with congestive heart failure, sinus rhythm (SR), and complete LBBB (6 females; age 61 +/- 8 years; NYHA class II-III; echocardiographic ejection fraction 21 +/- 5%; QRS > or = 130 ms) underwent permanent biventricular DDDR pacemaker implantation. A high-resolution 65-lead body-surface ECG recording was performed at baseline and during right-, left-, and biventricular pacing, and the total 65-lead root mean square curve of the QRST complex and the interlead QT dispersion were assessed. The QRS duration was increased during right (RV)- and left ventricular (LV) pacing (127 +/- 26% and 117 +/- 40%; P < 0.05), as compared to SR (100%) and biventricular pacing (93 +/- 16%; ns). The QTc interval was increased during RV and LV pacing (112 +/- 12% and 114 +/- 14%; P < 0.05) as compared to SR (100%) or biventricular pacing (99 +/- 12%). There was no effect on JT interval during all pacing modes. The T(peak-end) interval was increased during right (120 +/- 34%; P < 0.01) and LV pacing (113 +/- 29%; P < 0.05) but decreased during biventricular pacing (81 +/- 19%; P < 0.01). A similar effect was found for the T(peak-end) integral and the T(peak) amplitude. QT dispersion was increased during right ventricular (129 +/- 16 ms; P < 0.05) and decreased during biventricular pacing (90 +/- 12 ms; P < 0.01), as compared to SR (114 +/- 22 ms). CONCLUSIONS: Using a high-resolution surface ECG, biventricular pacing resulted in a significant reduction of ECG markers of ventricular dispersion of repolarization.

Multi-lead ECG electrode array for clinical application of electrocardiographic inverse problem.

Methods for noninvasive imaging of electric function of the heart might become clinical standard procedure the next years. Thus, the overall procedure has to meet clinical requirements as easy and fast application. In this study we propose a new electrode array which improves the information content in the ECG map, considering clinical constraints such as easy to apply and compatibility with routine leads. A major challenge is the development of an electrode array which yields a high information content even for a large interindividual variation in torso shape. For identifying regions of high information content we introduce the concept of a locally applied virtual electrode array. As a result of our analysis we constructed a new electrode array consisting of two L-shaped regular spaced parts and compared it to the electrode array we use for clinical studies upon activation time imaging. We assume that one side effect caused by the regular shape and spacing of the new array be that the reconstruction of electrodes placed on the patients back is simplified. It may be sufficient to record a few characteristic electrode positions and merge them with a model of the posterior array.

Simulation of atrial electrophysiology and body surface potentials for normal and abnormal rhythm.

The effect of different atrial electrical activation sequences (sinus rhythm and atrial flutter circling in the right atrium) on the body surface potentials is investigated in this study. A realistic volume conductor model consisting of atria, lungs, chest and blood masses is generated from image stacks recorded by magnetic resonance imaging. The electrical sources-the transmembrane potentials-within the atrial volumetric model are simulated for different atrial rhythms employing a cellular automaton capable of considering different parameters depending on the specific properties of the tissues. The potentials on the torso surface are computed applying the finite element method for solving the differential equations derived from the bidomain theory. Both the simulated atrial activation patterns and the computed torso potentials for atrial sinus rhythm and atrial flutter are in qualitatively and quantitatively good agreement with data measured in humans. The simulation of body surface potentials generated by different electrical activation sequences in the atria or ventricles allows testing and assessing noninvasive imaging of cardiac electrophysiology, as both the potentials on the body surface and the reference activation in the heart are available.

Atrial noninvasive activation mapping of paced rhythm data.

INTRODUCTION: Atrial arrhythmias have emerged as a topic of great interest for clinical electrophysiologists. Noninvasive imaging of electrical function in humans may be useful for computer-aided diagnosis and treatment of cardiac arrhythmias, which can be accomplished by the fusion of data from ECG mapping and magnetic resonance imaging (MRI). METHODS AND RESULTS: In this study, a bidomain-theory-based surface heart model activation time (AT) imaging approach was applied to paced rhythm data from four patients. Pacing sites were the right superior pulmonary vein, left inferior pulmonary vein, left superior pulmonary vein, coronary sinus, posterior wall of right atrium, and high right atrium. For coronary sinus pacing, the AT pattern of the right atrium was compared with a CARTO map. The root mean square error between CARTO geometry (85 nodal points) and the surface model of the right atrium was 8.6 mm. The correlation coefficient of the noninvasively obtained AT map of the right atrium and the CARTO map was 0.76. All pulmonary vein pacing sites were identified. The reconstructed pacing site of right posterior atrial pacing correlates with the invasively determined pacing catheter position with a localization distance of 4 mm. CONCLUSION: The individual anatomic model of the atria of each patient enables accurate noninvasive AT imaging within the atria, resulting in a localization error for the pacing sites within 10 mm. Our findings may have implications for imaging of atrial activity in patients with focal arrhythmias or focal triggers.