Application of the single moving dipole inverse solution to the study of the Wolff-Parkinson-White syndrome in man

Noninvasive identification of two lesions with local repolarization changes using two dipoles in inverse solution simulation study

BACKGROUND

The method for inverse localization and identification of two distinct simultaneous lesions with changed repolarization in the ventricular myocardium (two-vessel disease) is proposed and its robustness to errors in input data is tested in this simulation study.

METHOD

The inverse solution was obtained from the difference between STT integral body surface potential map computed with repolarization changes and the STT integral map from normal activation. In a numerical model of ventricles 48 cases of two simultaneous lesions and 48 cases of a single lesion were modeled. The effect of the lesions was taken to be represented by two dipoles. The input data were disturbed by three types of added noise. Twenty three characteristics of every obtained inverse solution were defined and four of them were used as the features in discriminant analysis task distinguishing the correct inverse solutions identifying two lesions.

RESULTS

The mean localization error for identified two lesions was 1.1±0.7cm. The sensitivity and specificity of quadratic discriminant analysis with cross-validation and feature selection was higher than 90%.

CONCLUSIONS

The combination of the inverse solution with two dipoles and discriminant analysis allows the identification of two simultaneous lesions without a priori information about the number of lesions.

Noninvasive imaging of 3-dimensional myocardial infarction from the inverse solution of equivalent current density in pathological hearts

We propose a new approach to noninvasively image the 3-D myocardial infarction (MI) substrates based on equivalent current density (ECD) distribution that is estimated from the body surface potential maps (BSPMs) during S-T segment. The MI substrates were identified using a predefined threshold of ECD. Computer simulations were performed to assess the performance with respect to: 1) MI locations; 2) MI sizes; 3) measurement noise; 4) numbers of BSPM electrodes; and 5) volume conductor modeling errors. A total of 114 sites of transmural infarctions, 91 sites of epicardial infarctions, and 36 sites of endocardial infarctions were simulated. The simulation results show that: 1) Under 205 electrodes and 10-μV noise, the averaged accuracies of imaging transmural MI are 83.4% for sensitivity, 82.2% for specificity, 65.0% for Dice's coefficient, and 6.5 mm for distances between the centers of gravity (DCG). 2) For epicardial infarction, the averaged imaging accuracies are 81.6% for sensitivity, 75.8% for specificity, 45.3% for Dice's coefficient, and 7.5 mm for DCG; while for endocardial infarction, the imaging accuracies are 80.0% for sensitivity, 77.0% for specificity, 39.2% for Dice's coefficient, and 10.4 mm for DCG. 3) A reasonably good imaging performance was obtained under higher noise levels, fewer BSPM electrodes, and mild volume conductor modeling errors. The present results suggest that this method has the potential to aid in the clinical identification of the MI substrates.

Usefulness of ventricular endocardial electric reconstruction from body surface potential maps to noninvasively localize ventricular ectopic activity in patients

As radio frequency (RF) catheter ablation becomes increasingly prevalent in the management of ventricular arrhythmia in patients, an accurate and rapid determination of the arrhythmogenic site is of important clinical interest. The aim of this study was to test the hypothesis that the inversely reconstructed ventricular endocardial current density distribution from body surface potential maps (BSPMs) can localize the regions critical for maintenance of a ventricular ectopic activity. Patients with isolated and monomorphic premature ventricular contractions (PVCs) were investigated by noninvasive BSPMs and subsequent invasive catheter mapping and ablation. Equivalent current density (CD) reconstruction (CDR) during symptomatic PVCs was obtained on the endocardial ventricular surface in six patients (four men, two women, years 23-77), and the origin of the spontaneous ectopic activity was localized at the location of the maximum CD value. Compared with the last (successful) ablation site (LAS), the mean and standard deviation of localization error of the CDR approach were 13.8 and 1.3 mm, respectively. In comparison, the distance between the LASs and the estimated locations of an equivalent single moving dipole in the heart was 25.5 ± 5.5 mm. The obtained CD distribution of activated sources extending from the catheter ablation site also showed a high consistency with the invasively recorded electroanatomical maps. The noninvasively reconstructed endocardial CD distribution is suitable to predict a region of interest containing or close to arrhythmia source, which may have the potential to guide RF catheter ablation.

Localization of endocardial ectopic activity by means of noninvasive endocardial surface current density reconstruction

Localization of the source of cardiac ectopic activity has direct clinical benefits for determining the location of the corresponding ectopic focus. In this study, a recently developed current-density (CD)-based localization approach was experimentally evaluated in noninvasively localizing the origin of the cardiac ectopic activity from body-surface potential maps (BSPMs) in a well-controlled experimental setting. The cardiac ectopic activities were induced in four well-controlled intact pigs by single-site pacing at various sites within the left ventricle (LV). In each pacing study, the origin of the induced ectopic activity was localized by reconstructing the CD distribution on the endocardial surface of the LV from the measured BSPMs and compared with the estimated single moving dipole (SMD) solution and precise pacing site (PS). Over the 60 analyzed beats corresponding to ten pacing sites (six for each), the mean and standard deviation of the distance between the locations of maximum CD value and the corresponding PSs were 16.9 mm and 4.6 mm, respectively. In comparison, the averaged distance between the SMD locations and the corresponding PSs was slightly larger (18.4 ± 3.4 mm). The obtained CD distribution of activated sources extending from the stimulus site also showed high consistency with the endocardial potential maps estimated by a minimally invasive endocardial mapping system. The present experimental results suggest that the CD method is able to locate the approximate site of the origin of a cardiac ectopic activity, and that the distribution of the CD can portray the propagation of early activation of an ectopic beat.

Equivalent moving dipole localization of cardiac ectopic activity in a swine model during pacing

Localization of the initial site of cardiac ectopic activity has direct clinical benefits for treating focal cardiac arrhythmias. The aim of the present study is to experimentally evaluate the performance of the equivalent moving dipole technique on noninvasively localizing the origin of the cardiac ectopic activity from the recorded body surface potential mapping (BSPM) data in a well-controlled experimental setting. The cardiac ectopic activities were induced in four well-controlled intact pigs by either single-site pacing or dual-site pacing within the ventricles. In each pacing study, the initiation sites of cardiac ectopic activity were localized by estimating the locations of a single moving dipole (SMD) or two moving dipoles (TMDs) from the measured BSPM data and compared with the precise pacing sites (PSs). For the single-site pacing, the averaged SMD localization error was 18.6 ± 3.8 mm over 16 sites, while the averaged distance between the TMD locations and the two corresponding PSs was slightly larger (24.9 ± 6.2 mm over five pairs of sites), both occurring at the onset of the QRS complex (10-25 ms following the pacing spike). The obtained SMD trajectories originated near the stimulus site and then traversed across the heart during the ventricular depolarization. The present experimental results show that the initial location of the moving dipole can provide the approximate site of origin of a cardiac ectopic activity in vivo, and that the migration of the dipole can portray the passage of an ectopic beat across the heart.

Applicability of the Single Equivalent Moving Dipole Model in an Infinite Homogeneous Medium to Identify Cardiac Electrical Sources: A Computer Simulation Study in a Realistic Anatomic Geometry Torso Model

We have previously proposed an inverse algorithm for fitting potentials due to an arbitrary bio-electrical source to a single equivalent moving dipole (SEMD) model. The algorithm achieves fast identification of the SEMD parameters by employing a SEMD model embedded in an infinite homogeneous volume conductor. However, this may lead to systematic error in the identification of the SEMD parameters. In this paper, we investigate the accuracy of the algorithm in a realistic anatomic geometry torso model (forward problem). Specifically, we investigate the effect of measurement noise, dipole position and electrode configuration in the accuracy of the algorithm. The boundary element method was used to calculate the forward potential distribution at multiple electrode positions on the body surface due to a point dipole in the heart. We have found that the position and not the number of electrodes as well as the site of the origin of the arrhythmia in the heart have a significant effect on the accuracy of the inverse algorithm, while the measurement noise does not. Finally, we have shown that the inverse algorithm preserves the topology of the source distribution in the heart, thus potentially allowing the cardiac electrophysiologist to efficiently and accurately guide the tip of the catheter to the ablation site.

Statistical accuracy of a moving equivalent dipole method to identify sites of origin of cardiac electrical activation

While radio frequency (RF) catheter ablation (RCA) procedures for treating ventricular arrhythmias have evolved significantly over the past several years, the use of RCA has been limited to treating slow ventricular tachycardias (VTs). In this paper, we present preliminary results from computer and animal studies to evaluate the accuracy of an algorithm that uses the single equivalent moving dipole (SEMD) model in an infinite homogeneous volume conductor to guide the RF catheter to the site of origin of the arrhythmia. Our method involves measuring body surface electrocardiographic (ECG) signals generated by arrhythmic activity and by bipolar current pulses emanating from a catheter tip, and representing each of them by a SEMD model source at each instant of the cardiac cycle, thus enabling rapid repositioning of the catheter tip requiring only a few cycles of the arrhythmia. We found that the SEMD model accurately reproduced body surface ECG signals with a correlation coefficients > 0.95. We used a variety of methods to estimate the uncertainty of the SEMD parameters due to measurement noise and found that at the time when the arrhythmia is mostly localized during the cardiac cycle, the estimates of the uncertainty of the spatial SEMD parameters (from ECG signals) are between 1 and 3 mm. We used pacing data from spatially separated epicardial sites in a swine model as surrogates for focal ventricular arrhythmic sources and found that the spatial SEMD estimates of the two pacing sites agreed with both their physical separation and orientation with respect to each other. In conclusion, our algorithm to estimate the SEMD parameters from body surface ECG can potentially be a useful method for rapidly positioning the catheter tip to the arrhythmic focus during an RCA procedure.

Noninvasive study of ventricular preexcitation using multichannel magnetocardiography

In clinical practice, noninvasive classification of ventricular preexcitation (VPX) is usually done with ECG algorithms, which provide only a qualitative localization of accessory pathways. Since 1984, single or multichannel magnetocardiography (MMCG) has been used for three-dimensional localization of VPX sites, but a systematic study comparing the results of ECG and MMCG methods was lacking. This study evaluated the reliability of MMCG in an unshielded electrophysiological catheterization laboratory, and compared VPX classification as achieved with the five most recent ECG algorithms with that obtained by MMCG mapping and imaging techniques. A nine-channel direct current superconducting quantum interference device (DC-SQUID) MMCG system (sensitivity is 20 fT/Hz0.5) was used for sequential MMCG from 36 points on the anterior chest wall, within an area 20 x 20 cm. Twenty-eight patients with Wolff-Parkinson-White syndrome were examined at least twice, on the same day or after several months to test the reproducibility of the measurements. In eight patients, the reproducibility of MMCG was also evaluated using different MCG instrumentation during maximal VPX and/or atrioventricular reentrant tachycardia induced by transesophageal atrial pacing via a nonmagnetic catheter. The results of VPX localization with ECG algorithms and MMCG were compared. Equivalent current dipole, effective magnetic dipole, and distributed currents imaging models were used for the inverse solution. MMCG classification of VPX was found to be more accurate than ECG methods, and also provided additional information for the identification of paraseptal pathways. Furthermore, in patients with complex activation patterns during the delta wave, distributed currents imaging revealed two different activation patterns, suggesting the existence of multiple accessory pathways.

Applicability of the single equivalent point dipole model to represent a spatially distributed bio-electrical source

Although the single equivalent point dipole model has been used to represent well-localised bio-electrical sources, in realistic situations the source is distributed. Consequently, position estimates of point dipoles determined by inverse algorithms suffer from systematic error due to the non-exact applicability of the inverse model. In realistic situations, this systematic error cannot be avoided, a limitation that is independent of the complexity of the torso model used. This study quantitatively investigates the intrinsic limitations in the assignment of a location to the equivalent dipole due to distributed electrical source. To simulate arrhythmic activity in the heart, a model of a wave of depolarisation spreading from a focal source over the surface of a spherical shell is used. The activity is represented by a sequence of concentric belt sources (obtained by slicing the shell with a sequence of parallel plane pairs), with constant dipole moment per unit length (circumferentially) directed parallel to the propagation direction. The distributed source is represented by N dipoles at equal arc lengths along the belt. The sum of the dipole potentials is calculated at predefined electrode locations. The inverse problem involves finding a single equivalent point dipole that best reproduces the electrode potentials due to the distributed source. The inverse problem is implemented by minimising the chi2 per degree of freedom. It is found that the trajectory traced by the equivalent dipole is sensitive to the location of the spherical shell relative to the fixed electrodes. It is shown that this trajectory does not coincide with the sequence of geometrical centres of the consecutive belt sources. For distributed sources within a bounded spherical medium, displaced from the sphere's centre by 40% of the sphere's radius, it is found that the error in the equivalent dipole location varies from 3 to 20% for sources with size between 5 and 50% of the sphere's radius. Finally, a method is devised to obtain the size of the distributed source during the cardiac cycle.

Frequency series expansion of an explicit solution for a dipole inside a conducting sphere at low frequency

The electromagnetic field generated by a current dipole situated at an arbitrary position inside a conducting sphere is derived using the expansions of the spherical vector wave functions. The first few terms in a series expansion of this field with respect to the frequency are given for the normal magnetic field (used in magnetoencephalogram) and the tangential electric field (used in electroencephalogram) outside the conducting sphere at low frequency. It is shown that the first correction term to the static solution is linear in the frequency omega (the second correction term is proportional to omega 3/2) and, thus, the static solution can be used as a good approximation for the solution at a very low frequency.

Value of epicardial potential maps in localizing pre-excitation sites for radiofrequency ablation. A simulation study

Using computer simulations, we systematically investigated the limitations of an inverse solution that employs the potential distribution on the epicardial surface as an equivalent source model in localizing pre-excitation sites in Wolff-Parkinson-White syndrome. A model of the human ventricular myocardium that features an anatomically accurate geometry, an intramural rotating anisotropy and a computational implementation of the excitation process based on electrotonic interactions among cells, was used to simulate body surface potential maps (BSPMs) for 35 pre-excitation sites positioned along the atrioventricular ring. Two individualized torso models were used to account for variations in torso boundaries. Epicardial potential maps (EPMs) were computed using the L-curve inverse solution. The measure for accuracy of the localization was the distance between a position of the minimum in the inverse EPMs and the actual site of pre-excitation in the ventricular model. When the volume conductor properties and lead positions of the torso were precisely known and the measurement noise was added to the simulated BSPMs, the minimum in the inverse EPMs was at 12 ms after the onset on average within 0.65 +/- 0.26 cm of the pre-excitation site. When the standard torso model was used to localize the sites of onset of the pre-excitation sequence initiated in individualized male and female torso models, the mean distance between the minimum and the pre-excitation site was 0.67 +/- 0.31 cm for the male torso and 0.82 +/- 0.53 cm for the female torso. The findings of our study indicate that a location of the minimum in EPMs computed using the inverse solution can offer non-invasive means for pre-interventional planning of the ablative treatment.

Accuracy of single-dipole inverse solution when localising ventricular pre-excitation sites: simulation study

Different factors are investigated that may affect the accuracy of an inverse solution that uses a single-dipole equivalent generator, in a standardised inhomogeneous torso model, when localising the pre-excitation sites. An anatomical model of the human ventricular myocardium is used to simulate body surface potential maps (BSPMs) and magnetic field maps (MFMs) for 35 pre-excitation sites positioned on the epicardial surface along the atrioventricular ring. The sites of pre-excitation activity are estimated by the single-dipole method, and the measure for the accuracy of the localisation is the localisation error, defined as the distance between the location of the best-fitting single dipole and the actual site of pre-excitation in the ventricular model. The findings indicate that, when the electrical properties of the volume conductor and lead positions are precisely known and the 'measurement' noise is added to the simulated BSPMs and MFMs, the single-dipole method optimally localises the pre-excitation activity 20 ms after the onset of pre-excitation, within 0.71 +/- 0.28 cm and 0.65 +/- 0.30 cm using BSPMs and MFMs, respectively. When the standard torso model is used to localise the sites of onset of the pre-excitation sequence initiated in four individualised torso models, the maximum errors are as high as 2.6-3.0 cm (even though the average error, for both the BSPM and MFM localisations, remains within the 1.0-1.5 cm range). In spite of these shortcomings, it is thought that single-dipole localisations can be useful for non-invasive pre-interventional planning.

A bioelectric inverse imaging technique based on surface Laplacians

A new approach is proposed to solve bioelectric inverse problems by employing the surface Laplacian of the bioelectrical potential. A theoretical investigation was conducted to test the feasibility of epicardial inverse imaging of cardiac electrical activity. A two-sphere homogeneous volume conductor model, where the inner sphere represents the epicardium and the outer sphere the body surface, was used. Radial and tangential current dipoles were used to approximate localized wavefronts propagating from the endocardium to the epicardium, and ectopic myocardial activities. The epicardial potential distribution was reconstructed from the body surface Laplacians with the aid of the Tikhonov zero-order regularization technique, which then was compared with the results obtained from the body surface potentials using the same regularization scheme. The two inverse solutions were compared qualitatively via visual inspection of the reconstructed epicardial potential maps, and quantitatively by examining relative errors and correlation coefficients between the "true" and the reconstructed epicardial potentials. Both qualitative and quantitative results indicate that the surface Laplacians play a positive role in improving the ill-posed nature of the bioelectric inverse problem, which would enhance our capability of reconstructing important epicardial events such as extrema in the epicardial potential distribution. The present theoretical study suggests that the Laplacian-based inverse imaging technique may have important applications to epicardial inverse imaging and other bioelectric inverse imaging.

Magnetocardiographically-guided catheter ablation

After more than 30 years since the first magnetocardiographic (MCG) recording was carried out with induction coils, MCG is now approaching the threshold of clinical use. During the last 5 years, in fact, there has been a growing interest of clinicians in this new method which provides an unrivalled accuracy for noninvasive, three-dimensional localization of intracardiac source. An increasing number of laboratories are reporting data validating the use of MCG as an effective method for preoperative localization of arrhythmogenic substrates and for planning the best catheter ablation approach for different arrhythmogenic substrates. In this article, available data from literature have been reviewed. We consider the clinical use of MCG to localize arrhythmogenic substrates in patients with Wolff-Parkinson-White syndrome and in patients with ventricular tachycardia in order to assess the state-of-the-art of the method on a large number of patients. This article also addresses some suggestions for industrial development of more compact, medically oriented MCG equipments at reasonable cost.

An equivalent body surface charge model representing three-dimensional bioelectrical activity

A new surface-source model has been developed to account for the bioelectrical potential on the body surface. A single-layer surface-charge model on the body surface has been developed to equivalently represent bioelectrical sources inside the body. The boundary conditions on the body surface are discussed in relation to the surface-charge in a half-space conductive medium. The equivalent body surface-charge is shown to be proportional to the normal component of the electric field on the body surface just outside the body. The spatial resolution of the equivalent surface-charge distribution appears intermediate between those of the body surface potential distribution and the body surface Laplacian distribution. An analytic relationship between the equivalent surface-charge and the surface Laplacian of the potential was found for a half-space conductive medium. The effects of finite spatial sampling and noise on the reconstruction of the equivalent surface-charge were evaluated by computer simulations. It was found through computer simulations that the reconstruction of the equivalent body surface-charge from the body surface Laplacian distribution is very stable against noise and finite spatial sampling. The present results suggest that the equivalent body surface-charge model may provide an additional insight to our understanding of bioelectric phenomena.

Solvability of the electrocardiology inverse problem for a moving dipole

New formulations of the direct and inverse problems for the moving dipole are offered. It has been suggested to limit the study by a small area on the chest surface. This lowers the role of the medium inhomogeneity. When formulating the direct problem, irregular components are considered. The algorithm of simultaneous determination of the dipole and regular noise parameters has been described and analytically investigated. It is shown that temporal overdetermination of the equations offers a single solution of the inverse problem for the four leads.

Localization of accessory pathways in Wolff-Parkinson-White syndrome by high-resolution magnetocardiographic mapping

Fifteen patients with Wolff-Parkinson-White syndrome were studied with standard 12-lead electrocardiogram, invasive electrophysiologic study, and high-resolution magnetocardiographic (MCG) mapping. In addition, intraoperative epicardial mapping was performed in seven surgically treated patients. The MCG characteristics of ventricular preexcitation for different locations of the atrioventricular accessory pathways were described in terms of morphology and field patterns. Three mathematical source models in semi-infinite conducting space were used for localization computations: the current dipole model, the truncated current multipole model and the magnetic dipole model. Finally, the localization results of MCG and invasive mappings and electrocardiograms were compared. The mean three-dimensional distance between the localization results obtained from MCG maps and electrophysiologic study was 3.9 cm for the magnetic dipole model, 4.8 cm for the truncated current multipole model, and 7.3 cm for the current dipole model. The corresponding distances in the seven intraoperatively mapped cases were 2.3 cm for the magnetic dipole model, 5.2 cm for the truncated current multipole model, and 6.3 cm for the current dipole model. In conclusion, noninvasive MCG mapping may significantly contribute to the invasive catheter mapping for optimal preoperative localization of preexcitation site and atrioventricular accessory pathways in Wolff-Parkinson-White syndrome.

Magnetocardiographic functional localization using current multipole models

High-resolution magnetocardiographic mapping was applied to localize the ventricular preexcitation site in ten patients suffering form Wolff-Parkinson-White syndrome. Three different source models were tested, consisting of the dipole and quadrupole moments in a general multipole expansion. Noninvasive localizations were performed by computations based on measured magnetic maps without a priori assumptions of the source location and without imposing any constraints. In all cases, the computed results were compared with invasive localization results obtained by catheter mapping technique. Preoperative catheterization localizes the atrial end of the accessory pathway, while our method localizes the ventricular preexcitation site. Of the models used, the average three-dimensional difference between the invasive localization results and the HR-MCG results was smallest 2.9 cm for the source model consisting of the magnetic dipole. The preexcitation site was very deep in all cases. The current dipole alone was inaccurate in estimating the source depth, but inclusion of the quadrupole moments improved the results. Two of the patients underwent surgery to interrupt the accessory pathway, which provided further validation for the noninvasive localizations.

Sample size and dimensionality in multivariate classification: implications for body surface potential mapping

This paper presents empirically determined guidelines for specifying the number of features appropriate for multivariate classification studies for given sample sizes. Sample size was considered adequate if the mean distance between two sample sets, taken from the same continuous multivariate distribution and projected onto the best separating direction, remained below a prescribed level. To quantitate the sample size requirement, homogeneity of sample set pairs of equal size. N, taken from the same continuous multivariate distribution was studied as a function of dimensionality. M. Homogeneity was characterized by the maximum absolute distances (Dmax) between the corresponding pairs of empirical cumulative probability distributions on the best separating projection. Computer generated data sets were used to estimate the cumulative probability distribution, P(D)M.N, for sample sizes, N, ranging from 5 to 100 and the dimensionality, M, ranging from 1 to 4. An empirical relationship between the estimated step-polygons and the Kolmogorov type one dimensional limiting distribution L(z) has been established. Based on the sample size data of 34 key papers on clinical body surface potential mapping (BSPM) it is noted that in 30% of the cases only one, and in 6% of the cases only two parameters could be used for statistical group representation to ensure a reasonable reliability (Dmax less than 0.2). In 56% of the published cases the sample sizes could not guarantee this reliability even for one feature or parameter.

Computer model of excitation and recovery in the anisotropic myocardium. II. Excitation in the simplified left ventricle

A computer model of propagated excitation and recovery in anisotropic cardiac tissue was described in the first report of this series. The model consists of a large number of excitable elements whose subthreshold interactions are governed by the anisotropic bidomain theory but whose suprathreshold behavior (action potential) is largely preassigned. As described in the first report, the model's performance was tested in rectangular and cubic arrays of excitable elements. This second report deals with three-dimensional simulations in a simplified left ventricle with anisotropy; specifically, the activation process in the "normal" ventricle is described (exemplified by the activation sequences started from various endocardial, intramural, and epicardial sites). To further substantiate our model's validity, we compare simulated epicardial and body-surface potential distributions with experimental findings in isolated canine hearts and with clinical evidence provided by electrocardiographic body-surface mapping.

Multicategory classification of body surface potential maps

A statistical classification method is suggested for body surface potential maps (BSPM). The initial data reduction utilizes the Fourier expansion and time integration, resulting in physiological-oriented features. Based on Fischer's criterion, optimal discriminant vectors are used to map the features to an optimal subdomain. Experimental criteria determine the dimensionality of the subdomain and the number of features to be mapped into it. Classification is performed in two steps. In the first, a k-nearest neighbor (k-NN) rule is used for every two-category problem, the results of which are fed into a voting rule for final classification. The method is tested with 123 patients divided into four categories: normal (NR), ischemia (IS), myocardial infarction (MI), and left bundle branch block (LB) patients. The success is between 88% (for IS) and 100% (for LB) for QRS segment integration. Departure maps were used to explain the misclassified patterns.

Anatomic localization of a single electrical source within the boundary of the human torso

A closed prolate ellipsoid was used to approximate the surface of the Rush torso model to permit recovery of the site and orientation of known dipoles in 15 cardiac locations. Localization was found to be reasonably close, usually within 2 cm. When body surface potential maps of 37 subjects with right ventricular pacemakers were similarly treated, the discrepancy between known pacemaker site and the site of earliest activation was relatively large (mean, greater than 4 cm) and rapidly increased within the ensuing millisecond. The discrepancy not only emphasizes the wide variation in body shape and tissue distribution in living subjects, but also points to probable physical separation between stimulus site and earliest detectable activation site because of ischemia, infarction, or myocardial response to variation in current strength of the stimulus.

Application of dipole analysis for the diagnosis of myocardial infarction in the presence of left bundle branch block

The residue value on dipole analysis (the ratio of non-dipolar component to the measured body surface potentials) was estimated mathematically in 16 patients with left bundle branch block. Patients were classified into those with (group A, nine patients) and those without (group B, seven patients) a perfusion defect on thallium-201 myocardial scintigraphy. For the entire QRS complex the residue of group B was smaller than that of normal subjects (20.0 +/- 4.1% versus 24.6 +/- 3.5%, p less than 0.05). Group A showed a greater mean residue value than group B (27.4 +/- 4.4% versus 20.3 +/- 2.4%, p less than 0.01) only during the initial one-third of the QRS complex. All but one patient of group A and only one patient in group B showed a high peak on the residue curve during the initial stage of the QRS complex. The maximal residue value of group A during the initial QRS complex was significantly greater than that of group B (40.9 +/- 10.9% versus 23.4 +/- 5.4%, p less than 0.01). An arbitrarily selected criterion of the maximal residue value greater than or equal to 30% during the initial QRS complex showed a sensitivity of 89% with a specificity of 86% for the diagnosis of myocardial infarction in the presence of left bundle branch block. These results might be related to the complex ventricular activation around the infarcted area even in the presence of left bundle branch block in which intramyocardial conduction with a simple activation front predominates. Dipole analysis appeared to be a valuable method of diagnosing myocardial infarction in the presence of left bundle branch block.

Value of the resting 12 lead electrocardiogram and vectorcardiogram for locating the accessory pathway in patients with the Wolff-Parkinson-White syndrome

The resting 12 lead electrocardiogram and vectocardiogram were reviewed in 47 patients with the Wolff-Parkinson-White syndrome (a) who had pre-excitation on the resting 12 lead electrocardiogram, (b) who had a single anterograde conducting accessory pathway assessed and located during preoperative electrophysiological study and during epicardial mapping at operation, and (c) in whom surgical division of the accessory pathway resulted in loss of pre-excitation. The site of the accessory pathway established during operation was compared with that established by evaluating the polarity of the delta wave and QRS complex on the resting 12 lead electrocardiogram. The electrocardiogram was assessed by the Rosenbaum criteria (Wolff-Parkinson-White type A, left-sided pathway; or type B, right-sided pathway), the Gallagher criteria (atrial pacing resulting in maximal pre-excitation), and the World Health Organisation criteria (a composite of previous studies). The Gallagher and World Health Organisation criteria were derived from patients demonstrating maximal pre-excitation that often required atrial pacing. The present study was designed to determine whether these criteria could be accurately applied to the resting 12 lead electrocardiogram on which the degree of pre-excitation was variable. The Rosenbaum criteria correctly identified a left sided accessory pathway in 26 of 34 patients and a right-sided accessory pathway in nine of 13 patients. The Gallagher and World Health Organisation criteria correctly identified the location in only 15 (32%) of the 47 patients. The resting vectorcardiogram was inaccurate for locating the accessory pathway. Although published criteria are useful for identifying the site of the accessory pathway from an electrocardiogram obtained when rapid atrial pacing is being used to achieve maximal pre-excitation, they are not suitable for identifying the exact site of an accessory pathway from the resting 12 lead electrocardiogram.

Localization of cardiac ectopic activity in man by a single moving dipole. Comparison of different computation techniques

The accuracy of different computation techniques for the non-invasive localization of cardiac ectopic activity was evaluated. Body surface potentials were recorded from 63 leads in 14 patients with implanted pacemakers. The location, orientation and magnitude of a single moving dipole (SMD) were computed from the first eight terms of a truncated multipole expansion estimated from the body surface potentials. The SMD trajectories obtained during the QRS complex were plotted along with the heart outlines and pacing leads obtained independently from chest x-rays. The origin of the SMD trajectories was compared to the position of the pacing lead to evaluate the accuracy of the SMD. The optimum computation technique used a least-squares (LS) estimation of the multipole expansion truncated at 15 multipoles, in conjunction with a torso model that included regions of lower conductivity representing the lungs. With this method, the SMD trajectories originated near the pacing lead (25 +/- 12 mm) and adequately represented the progression of the ectopic wavefront across the entire heart silhouette. With the LS techniques using 8 or 24 multipoles, the spans of the trajectories were respectively too short, or too long to cover the heart, and the average distance between the SMD at QRS onset and the pacing lead was larger. With a surface integration technique, the SMD-pacing lead distances were similar, both for a finite homogeneous torso model with a fixed geometry, as well as for torso models adapted to the torso geometry of each patient. The SMD was found adequate to represent the progression of an ectopic wavefront, and to localize its origin in man.