Spatial resolution of body surface potential maps and magnetic field maps: a simulation study applied to the identification of ventricular pre-excitation sites

The magnetocardiogram

The magnetic field produced by the heart's electrical activity is called the magnetocardiogram (MCG). The first 20 years of MCG research established most of the concepts, instrumentation, and computational algorithms in the field. Additional insights into fundamental mechanisms of biomagnetism were gained by studying isolated hearts or even isolated pieces of cardiac tissue. Much effort has gone into calculating the MCG using computer models, including solving the inverse problem of deducing the bioelectric sources from biomagnetic measurements. Recently, most magnetocardiographic research has focused on clinical applications, driven in part by new technologies to measure weak biomagnetic fields.

Biomagnetism: The First Sixty Years

Biomagnetism is the measurement of the weak magnetic fields produced by nerves and muscle. The magnetic field of the heart-the magnetocardiogram (MCG)-is the largest biomagnetic signal generated by the body and was the first measured. Magnetic fields have been detected from isolated tissue, such as a peripheral nerve or cardiac muscle, and these studies have provided insights into the fundamental properties of biomagnetism. The magnetic field of the brain-the magnetoencephalogram (MEG)-has generated much interest and has potential clinical applications to epilepsy, migraine, and psychiatric disorders. The biomagnetic inverse problem, calculating the electrical sources inside the brain from magnetic field recordings made outside the head, is difficult, but several techniques have been introduced to solve it. Traditionally, biomagnetic fields are recorded using superconducting quantum interference device (SQUID) magnetometers, but recently, new sensors have been developed that allow magnetic measurements without the cryogenic technology required for SQUIDs.

Investigations of sensitivity and resolution of ECG and MCG in a realistically shaped thorax model

Solving the inverse problem of electrocardiography (ECG) and magnetocardiography (MCG) is often referred to as cardiac source imaging. Spatial properties of ECG and MCG as imaging systems are, however, not well known. In this modelling study, we investigate the sensitivity and point-spread function (PSF) of ECG, MCG, and combined ECG+MCG as a function of source position and orientation, globally around the ventricles: signal topographies are modelled using a realistically-shaped volume conductor model, and the inverse problem is solved using a distributed source model and linear source estimation with minimal use of prior information. The results show that the sensitivity depends not only on the modality but also on the location and orientation of the source and that the sensitivity distribution is clearly reflected in the PSF. MCG can better characterize tangential anterior sources (with respect to the heart surface), while ECG excels with normally-oriented and posterior sources. Compared to either modality used alone, the sensitivity of combined ECG+MCG is less dependent on source orientation per source location, leading to better source estimates. Thus, for maximal sensitivity and optimal source estimation, the electric and magnetic measurements should be combined.

From 3D to 4D imaging: is that useful for interventional cardiac electrophysiology?

Three-dimensional electroanatomical imaging is increasingly used in interventional cardiac electrophysiology, to guide catheter ablation of cardiac arrhythmias. At the same time, there is a growing interest for non-invasive methods, such as magnetocardiographic mapping (MCG), to localize the arrhythmogenic substrates, to test their reproducibility and to plan the most appropriate interventional approach. So far electroanatomical imaging has relayed on static mathematical modeling of the heart and more recently on direct merging with three-dimensional rendering of cardiac anatomy from multidetector computer tomography or magnetic resonance imaging. Merging electrophysiological information with static anatomical structures, can surely be a source of uncertainty for MCG-based pre-interventional localization of the arrhythmogenic substrate and causes mismatch between the real-time imaging of moving catheters and the static geometry of the cardiac chambers reconstructed with invasive electroanatomical imaging. The implementation of recent realistic numerical models of the beating heart in a breathing thorax can improve accuracy and fill the gap between non-invasive and interventional electroanatomical imaging.

Value and limitations of an inverse solution for two equivalent dipoles in localising dual accessory pathways

Investigations were carried out into whether an equivalent generator consisting of two dipoles could be used to detect dual sites of ventricular activity. A computer model of the human ventricular myocardium was used to simulate activation sequences initiated at eight different pairs of sites positioned on the epicardial surface of the atrio-ventricular ring. From these sequences, 117-lead body surface potentials (covering the anterior and posterior torso), 64-lead magnetic field maps (above the anterior chest) and 128-lead magnetic field maps (above the anterior and posterior chest) were simulated and were then used to localise dual accessory pathways employing pairs of equivalent dipoles. Average localisation errors were 12 mm, 12 mm and 9 mm, respectively, when body surface potentials, 64-lead and 128-lead magnetic fields were used. The results of the study suggest that solving the inverse problem for two dipoles could provide additional information on dual accessory pathways prior to electrophysiological study.

Classifying cases of fetal Wolff-Parkinson-White syndrome by estimating the accessory pathway from fetal magnetocardiograms

The paper presents an evaluation of the possibility of using fetal magnetocardiogram (FMCG) signals to estimate and classify the accessory pathway in fetal Wolff-Parkinson-White (WPW) syndrome. The FMCG signals of two fetuses with WPW syndrome (type A) were detected using a 64-channel superconducting quantum-interference device system. An average across the cycles of these signals was taken to obtain clear WPW signals. To determine the direction and position of the accessory pathway in a fetal heart accurately, the accessory pathway and activated pathway at the peak of the QRS complex thus obtained were estimated for each fetus, using a single-dipole model. The phase angle (about 90 degrees) between the equivalent current dipoles (ECDs) was the same for both fetuses. This angle suggested that the accessory pathway is in the left side of the heart, i.e. that the pathway exists in the left ventricle, which indicates type A WPW syndrome. Identification of the position of the accessory pathway in a fetus with WPW syndrome from the angle between the ECD of the accessory pathway and the ECD of the peak in the QRS complex was thus demonstrated.

Noninvasive characterisation of multiple ventricular events using electrocardiographic imaging

Distributions of epicardial potentials, calculated from body surface electrocardiograms (ECGs), were investigated to determine if they could enable detection of multiple sites of ventricular activity. An anatomical model of the human ventricular myocardium was used to simulate activation sequences initiated at nine different ventricular pairs of sites. From these sequences, body surface ECGs were simulated at 352 sites on the torso surface and then used to reconstruct epicardial potentials at 202 sites. The criterion for detection of dual ventricular events was the presence of two distinct primary potential minima in the reconstructed epicardial potentials. The shortest distance between the two events in the right ventricle that resulted in the reconstruction of epicardial potential patterns, featuring two minima, was 27 mm; the distance between the two events in the left ventricle was 23 mm. When Gaussian white noise in the simulated body surface potentials was increased from 3 microV to 15 microV and 50 microV, dual events became more difficult to distinguish. Findings indicate that calculated epicardial potentials provide useful visual information about the presence of multiple ventricular events that is not apparent in features of body surface ECGs, and could be particularly helpful in optimising mapping procedures during difficult or unsuccessful radiofrequency ablations of accessory pathways.

A method for detecting myocardial abnormality by using a total current-vector calculated from ST-segment deviation of a magnetocardiogram signal

A simple method to determine the state of ischaemia or fibrosis of myocardial cells has been developed. This method uses the ST wave of 64-channel magnetocardiogram (MCG) signals to calculate three parameters from the current-arrow map of the normal component signal of the MCG. One parameter is a total current vector that is obtained through summation of all current arrows. Another is a variance current vector calculated from the differential vector of two total current vectors at different times. The third is a flatness factor between the magnitude of the total current vector and the variance current vector. The three parameters are independent of the distance between the heart and the gradiometers. We measured the MCG signals of 29 healthy subjects, twenty patients with coronary artery disease (ten with previous myocardial infarction (MI) and ten with angina pectoris (AP)), and eight patients with cardiomyopathy (four with hypertrophic cardiomyopathy (HCM), three with dilated cardiomyopathy (DCM), and one with restrictive cardiomyopathy (RCM)). With our method, none of the healthy subjects tested positive for myocardial abnormalities, while 80% of the MI patients, 50% of the AP patients, and 100% of the cardiomyopathy patients tested positive. Although further testing is needed, we feel this simple technique enables easy diagnosis of myocardial damage.

The effect of geometric and topologic differences in boundary element models on magnetocardiographic localization accuracy

This study was performed to evaluate the changes in magnetocardiographic (MCG) source localization results when the geometry and the topology of the volume conductor model were altered. Boundary element volume conductor models of three patients were first constructed. These so-called reference torso models were then manipulated to mimic various sources of error in the measurement and analysis procedures. Next, equivalent current dipole localizations were calculated from simulated and measured multichannel MCG data. The localizations obtained with the reference models were regarded as the "gold standard." The effect of each modification was investigated by calculating three-dimensional distances from the gold standard localizations to the locations obtained with the modified model. The results show that the effect of the lungs and the intra-ventricular blood masses is significant for deep source locations and, therefore, the torso model should preferably contain internal inhomogeneities. However, superficial sources could be localized within a few millimeters even with nonindividual, so called standard torso models. In addition, the torso model should extend long enough in the pelvic region, and the positions of the lungs and the ventricles inside the model should be known in order to obtain accurate localizations.

Continuous localization of cardiac activation sites using a database of multichannel ECG recordings

Monomorphic ventricular tachycardia and ventricular extrasystoles have a specific exit site that can be localized using the multichannel surface electrocardiogram (ECG) and a database of paced ECG recordings. An algorithm is presented that improves on previous methods by providing a continuous estimate of the coordinates of the exit site instead of selecting one out of 25 predetermined segments. The accuracy improvement is greatest, and most useful, when adjacent pacing sites in individual patients are localized relative to each other. Important advantages of the new method are the objectivity and reproducibility of the localization results.

Localization of intramural necrotic regions using electrocardiographic imaging

Recent studies have demonstrated that electrocardiographic imaging (ECGI) is a novel noninvasive modality for exploring the spread of electrical activation within the ventricular wall. In this study, our goal was to explore the ability of ECGI in reconstructing epicardial potentials and electrograms in the ventricles damaged by localized necroses (<2 cm2). An anatomical model of the human ventricular myocardium was used to simulate activation sequences initiated at 428 epicardial and endocardial pacing sites distributed over the right ventricular and left ventricular free walls. From these realistic sequences, we simulated extracardiac potentials at epicardial (202 sites) and torso surfaces (352 sites) using boundary element model of the human torso. ECGI in terms of the L-curve was applied to compute epicardial potentials and unipolar electrograms (202 sites). Inversely computed electrograms correlated well with those simulated by an anatomical model (r > 0.9 at 68% of sites). Specifically, ECGI accurately reconstructed the following features that have been observed during measurements on the exposed canine hearts: (a) an epicardial potential pattern with a central minimum and two maxima, with the minimum positioned above the pacing site; (b) a complete transient loss of one of the positive areas in the epicardial potential pattern when the necrosis was located subepicardially; and (c) a transient gap in the expanding positive areas of the epicardial potential pattern when the necrosis was located intramurally or subendocardially. Findings of our study indicate that ECGI provides detailed reconstruction of patterns of myocardial activation in the presence of localized necroses and may be useful in the assessment of arrhythmogenic substrate in the clinical setting.

Bioelectromagnetic localization of a pacing catheter in the heart

The accuracy of localizing source currents within the human heart by non-invasive magneto- and electrocardiographic methods was investigated in 10 patients. A non-magnetic stimulation catheter inside the heart served as a reference current source. Biplane fluoroscopic imaging with lead ball markers was used to record the catheter position. Simultaneous multichannel magnetocardiographic (MCG) and body surface potential mapping (BSPM) recordings were performed during catheter pacing. Equivalent current dipole localizations were computed from MCG and BSPM data, employing standard and patient-specific boundary element torso models. Using individual models with the lungs included, the average MCG localization error was 7+/-3 mm, whereas the average BSPM localization error was 25+/-4 mm. In the simplified case of a single homogeneous standard torso model, an average error of 9+/-3 mm was obtained from MCG recordings. The MCG localization accuracies obtained in this study imply that the capability of multichannel MCG to locate dipolar sources is sufficient for clinical purposes, even without constructing individual torso models from x-ray or from magnetic resonance images.

Value of magnetocardiographic QRST integral maps in the identification of patients at risk of ventricular arrhythmias

It has been shown that regional ventricular repolarization properties can be reflected in body surface distributions of electrocardiographic QRST deflection areas (integrals). We hypothesize that these properties can be reflected also in the magnetocardiographic QRST areas and that this may be useful for predicting vulnerability to ventricular tachyarrhythmias. Magnetic field maps were obtained during sinus rhythm from 49 leads above the anterior chest in 22 healthy (asymptomatic) control subjects (group A) and in 29 patients with ventricular arrhythmias (group B). In each subject, the QRST deflection area was calculated for each lead and displayed as an integral map. The mean value of maximum was significantly larger in the control group A than in the patient group B (1,626+/-694 pTms vs. 582+/-547 pTms, P<0.0001). To quantitatively assess intragroup variability in the control group A and intergroup variability of the control and patient groups, we used the correlation coefficient r and covariance sigma. These indices showed significantly less intragroup than intergroup variation (e.g., in terms of sigma, 28.0x10(-6)+/-12.3x10(-6) vs. 3.4x10(-6)+/-12.5x10(-6), P<0.0001). Each QRST integral map was also represented as a weighted sum of 24 basis functions (eigenvectors) by means of Karhunen-Loeve transformation to calculate the contribution of the nondipolar eigenvectors (all eigenvectors beyond the third). This percentage nondipolar content of magnetocardiographic QRST integral maps was significantly higher in the patient group B than in the control group A (13.0%+/-9.1 % vs. 2.6%+/-2.0%, P<0.0001). Discriminations between control subjects and patients with ventricular arrhythmias based on magnitude of the maximum, covariance sigma, and nondipolar content were 90.2%, 90.2%, and 86.3% accurate, with a sensitivity of 89.7%, 93.1%, and 75.9%, and a specificity of 90.9%, 86.4%, and 100%. We have shown that magnitude of the maximum and indices of variability and nondipolarity of the magnetocardiographic QRST integral maps may predict arrhythmia vulnerability. This finding is in agreement with earlier studies that used body surface potential mapping and suggests that magneticfield mapping may also be a useful diagnostic tool for risk analysis.

Assessment of spatial resolution of pace mapping when using body surface potentials

Using computer simulations and statistical methods, the resolution of pace mapping when used in combination with body surface potentials was systematically investigated. In an anatomical model of the human ventricular myocardium, pre-excitation sequences were initiated at 69 sites positioned along the atrioventricular (AV) ring and corresponding body surface potential maps (BSPMs) were calculated at 32 leads placed on the anterior torso. For each time after the onset of pre-excitation (every 4 ms to 40 ms) and each root-mean-square (RMS) noise level (5, 10, 20 and 50 microV), BSPMs were cros-correlated and the spatial resolution defined as the largest pacing site separation at which the differences in correlation coefficients were not statistically significant (level p > or = 0.05). The findings indicate that when random RMS noise of 5 microV was added to the simulated BSPMs, average spatial resolution over all 60 sites was at 20 ms after the onset of pre-excitation within 3.5 +/- 0.9 mm. The results provide theoretical evidence that statistical analysis of BSPMs obtained during pace mapping can offer improved means for subcentimetre identification of accessory pathways located along the AV ring.

A comparison of simulated QRS isointegral maps resulting from pacing at adjacent sites: implications for the spatial resolution of pace mapping using body surface potentials

The precise localization of ventricular tachycardia (VT) foci is a prerequisite for the successful radiofrequency catheter ablation in patients. The purpose of this study was to systematically quantify over what distance adjacent sites in the right ventricular (RV) and left ventricular (LV) epicardium and LV endocardium could be distinguished by inspecting morphological features of QRS isointegral maps using statistical methods. We investigated the spatial resolution of QRS isointegral maps by means of an anatomically accurate computer model of the human ventricular myocardium that incorporates a bidomain model for simulating the realistic activation sequences and the oblique dipole model in combination with the boundary element method for calculating extracardiac potentials. In this model, we initiated activation sequences at a total of 183 epicardial and 75 LV endocardial pacing sites, positioned in three levels (basal, middle, and apical). For each of the 258 pacing sites, we calculated a set of 10 QRS isointegral maps with added Gaussian noise at 117 leads (covering the anterior and posterior torso) and at 32 leads (covering only the anterior torso), respectively. Sets of maps were then cross correlated and root-mean-square (RMS) values of difference maps were calculated for all possible pairs of pacing sites on the same level. We applied the nonparametric unpaired Kolmogorov-Smirnov test and defined the spatial resolution as the pacing site separation at which the differences in correlation coefficients and RMS differences were significant (level P < .05). We observed significant differences in maps when the distances between pacing sites were on average (+/- SD) greater than 4.3 +/- 1.0 mm. In more than 90% of pacing sites, the significant differences in maps were observed within 4 mm even when using a 32-lead mapping system. The findings of our study provide theoretical evidence that QRS isointegral maps may offer noninvasive means for preinterventional planning of the ablative treatment in localizing both endocardial and epicardial sites of origin of VT.

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.