Dobutamine stress magnetic resonance imaging suffices for the demonstration of myocardial ischaemia and viability.

We report three patients in whom dobutamine stress magnetic imaging (DS-MRI) was essential in assessing myocardial ischaemia. Two patients were referred to the cardiologist because of chest pain. Patient A had typical exertional angina and a normal resting electrocardiogram (ECG). Patient B had typical exercise-induced angina and had recently experienced an attack of severe chest pain at rest for 15 minutes. The ECG showed a complete left bundle branch block (LBBB). Patient C was referred for heart failure of unknown origin. There were no symptoms of chest pain during rest or exercise. Echocardiography in this patient demonstrated global left ventricular (LV) dilatation, systolic dysfunction and a small dyskinetic segment in the inferior wall. In all these patients exercise stress testing had failed to demonstrate myocardial ischaemia. Patients A and C produced normal findings whereas in patient B the abnormal repolarisation due to pre-existent LBBB precluded a diagnosis of ischaemia. Breath-hold DS-MRI was performed to study LV wall motion and wall thickening at rest through increasing doses of dobutamine. A test was considered positive for myocardial ischaemia if wall motion abnormalities developed at high-dose levels of the drug (20 µg/kg/min or more with a maximum of 40 µg/kg/min) in previously normal vascular territories or worsened in a segment that was normal at baseline. Recovery of wall thickening in a previously hypokinetic or akinetic segment at a low dose of dobutamine (5-10 µg/kg/min) was taken as proof of viability. Patients A and B developed hypokinesia progressing into akinesia at high-dose dobutamine in the anteroseptal area of the LV indicative of ischaemia. These findings were corroborated by coronary angiography demonstrating severe coronary artery disease which led to coronary artery bypass grafting (CABG) in patient A and balloon angioplasty in patient B. In patient C global recovery of LV contractions during low-dose dobutamine was followed by hypokinesia in the inferoseptal area during high-dose dobutamine. This biphasic response indicates myocardial viability as well as ischaemia. CABG was carried out because of multiple stenoses in the left coronary artery. Post-operatively LV function normalised. DS-MRI is a valuable method for detecting myocardial ischaemia and viability in patients with suspected coronary artery, and can be applied in every hospital with MRI equipment at its disposal.

[Dobutamine stress magnetic resonance imaging (DS-MRI), a valuable tool for the diagnosis of ischemic heart disease]. / Dobutaminestress-MRI van het hart: een waardevolle onderzoekstechniek voor de diagnostiek van ischemische hartziekten.

OBJECTIVE: Assessment of the clinical applicability of DS-MRI for the detection of myocardial ischemia and myocardial viability. DESIGN: Prospective. METHOD: In the period from 1 November 1999 to 31 October 2000, patients with suspected coronary artery disease who could not be studied by means of conventional bicycle ergometry underwent breath-hold DS-MRI (1 Tesla) 4 days after cessation of anti-ischemic medication. Three left ventricular short-axis planes were examined for the occurrence of disorders in wall movement during infusion of increasing doses of dobutamine (10, 20, 30 and 40 micrograms/kg/min). Temporary recovery of wall thickening in a previously diminished or non-contracting segment under 5 micrograms/kg/min of dobutamine was considered proof of viability. Development of hypo-, a- or dyskinesia at higher doses of dobutamine was taken to indicate ischemia. If the DS-MRI test was positive for ischemia, coronary angiography was performed. If indicated, this was followed by revascularisation. If DS-MRI did not reveal ischemia, the patient was seen at the outpatient department. RESULTS: Of the 100 patients (62 men and 38 women with an average age of 62 years, SD = 12) subjected to DS-MRI, 95 yielded results that were suitable for diagnosis. Of the 42 patients with DS-MRI scans that were considered positive for ischemia and in whom coronary angiography was subsequently performed, 41 had such coronary abnormalities that revascularisation was indicated. One patient was false-positive. All 53 patients with non-ischemic DS-MRI scans were followed-up for 11-23 months (mean 17 months). One patient died suddenly 2 weeks after the MRI-test. The other 52 patients did not experience any coronary events nor sudden cardiac death. The predictive value of a positive DS-MRI scan for ischemia was 98% and the predictive value of a negative DS-MRI scan was also 98%. CONCLUSION: DS-MRI is a safe diagnostic method for the detection or exclusion of myocardial ischemia and viability in patients with suspected coronary artery disease.

Employment of intra-individual variability to improve computerized ECG interpretation.

One of the reasons for the limited practical utility of computer programs for interpretation of electrocardiograms (ECGs) is their susceptibility to intra-individual variability. Two of the most prominent sources of intra-individual variability in ECGs, electrode placement variations and respiration, were studied for their effects on computerized ECG interpretation. Previous research has shown that the effects of intra-individual variability on computerized ECG interpretation depend largely on the individual ECG. To enable the assessment of chest electrode position variations for individual standard 12- lead ECGs, ECGs resulting from simulations of such position variations were interpreted. Variability due to respiration was assessed by interpreting all individual ECG beats instead of an averaged beat. In this paper two methods are presented that employ information about the intra-individual variability in individual ECGs. The first method provides an estimate of the reliability of the interpretation, the second attempts to improve the interpretation itself. In the first method we quantified the variation in interpretation caused by the two sources of intra-individual variability with the use of a stability index, a high index value indicating a low variation in interpretation. This index was subsequently studied using two sets of ECGs. For the first set a â clinical' reference interpretation was obtained from discharge letters. For the second set three cardiologists provided a â cardiologists' reference. The performance of subgroups of ECGs having stability indices higher than a particular value was computed. It appeared that for the â cardiologists' reference, the interpretations of ECGs with a high stability index were more often correct. No effect was found for the â clinical' reference. In the second method we attempted to improve the original interpretation by combining the alternative interpretations into a new interpretation. This was done by taking the median or the average of the quantified alternatives. These combined interpretations proved to perform better than the original interpretation when a cardiologist's interpretation was taken as a reference. This paper shows that intra-individual ECG variability can be used to improve original interpretations. This can be done without having to record multiple ECGs, provided that a model is available to simulate intra-individual variability. The presented methods do not depend on the classification algorithm that is used. They can be used both during classifier design to correct imperfections, and in routine use of the classifier to produce more representative classifications.

The coming of age of computerized ECG processing: can it replace the cardiologist in epidemiological studies and clinical trials?

In spite of decades of research and widespread use of computer programs for the analysis of electrocardiograms (ECGs), the accuracy and usefulness of computerized ECG processing has been questioned. To determine whether ECG computer programs can replace cardiologists in epidemiological studies and clinical trials, we reviewed the literature for evidence, concentrating on one influential ECG measurement, viz. QT interval duration, and one classification method, the Minnesota Code, which is the de facto standard for ECG coding. We compared interobserver variabilities of cardiologists with differences between computer programs and cardiologists, in order not to prejudice against the computer. Studies that contain this type of information indicate that interobserver variabilities are at least as large as differences between computer and cardiologist. This suggests that ECG computer programs perform at least equally well as human observers in ECG measurement and coding, and can replace the cardiologist in epidemiological studies and clinical trials.

Accurate automatic detection of electrode interchange in the electrocardiogram.

In recording an electrocardiogram (ECG), an interchange of electrodes may easily go unnoticed. Automatic detection would be desirable, but current algorithms, when dealing with more than left arm-right arm reversal, have moderate sensitivity. We propose a novel approach that uses the redundancy of information in the standard 12-lead ECG. We assume that each of the 8 independent electrocardiographic leads can be reconstructed from the 7 others in reasonable approximation. The correlation between any electrocardiographic lead and its reconstruction should be higher if the electrodes are correctly placed than when some interchange were present. The difference in correlation should have discriminative power. This was verified on a set of 3,305 ECGs for 14 common electrode interchange errors. The material was split in a learning and test set, and general reconstruction coefficients were computed from the learning set. For each interchange, electrode-error ECGs were derived by rearranging leads of the unaltered ECGs. Correlations between the actual leads and their reconstructions were computed for all ECGs. From the differences in lead correlation, decision rules were derived for each kind of interchange. All 14 rules had specificities of > or =99.5% in the test set. Sensitivities were > or =93% for 11 rules, and left arm-left leg electrode reversal scored low.

Computerized ST depression analysis improves prediction of all-cause and cardiovascular mortality: the strong heart study.

BACKGROUND: Nonspecific ST depression assessed by standard visual Minnesota coding (MC) has been demonstrated to predict risk. Although computer analysis has been applied to digital ECGs for MC, the prognostic value of computerized MC and computerized ST depression analyses have not been examined in relation to standard visual MC. METHODS: The predictive value of nonspecific ST depression as determined by visual and computerized MC codes 4.2 or 4.3 was compared with computer-measured ST depression >or= 50 microV in 2,127 American Indian participants in the first Strong Heart Study examination. Computerized MC and ST depression were determined using separate computerized-ECG analysis programs and visual MC was performed by an experienced ECG core laboratory. RESULTS: The prevalence of MC 4.2 or 4.3 by computer was higher than by visual analysis (6.4 vs 4.4%, P < 0.001). After mean follow-up of 3.7 +/- 0.9 years, there were 73 cardiovascular deaths and 227 deaths from all causes. In univariate Cox analyses, visual MC (relative risk [RR] 4.8, 95% confidence interval [CI] 2.6-9.1), computerized MC (RR 6.0, 95% CI 3.5-10.3), and computer-measured ST depression (RR 7.6, 95% CI 4.5-12.9) were all significant predictors of cardiovascular death. In separate multivariate Cox regression analyses that included age, sex, diabetes, HDL and LDL cholesterol, body mass index, systolic and diastolic blood pressure, microalbuminuria, smoking, and the presence of coronary heart disease, computerized MC (RR 3.0, 95% CI 1.6-5.6) and computer-measured ST depression (RR 3.1, 95% CI 1.7-5.7), but not visual MC, remained significant predictors of cardiovascular mortality. When both computerized MC and computer-measured ST depression were entered into the multivariate Cox regression, each variable provided independent risk stratification (RR 2.1, 95% CI 1.0-4.4, and RR 2.1, 95% CI 1.0-4.4, respectively). Similarly, computerized MC and computer-measured ST depression, but not visual MC, were independent predictors of all-cause mortality after controlling for standard risk factors. CONCLUSIONS: Computer analysis of the ECG, using computerized MC and computer-measured ST depression, provides independent and additive risk stratification for cardiovascular and all-cause mortality, and improves risk stratification compared with visual MC. These findings support the use of routine computer analysis of ST depression on the rest ECG for assessment of risk and suggest that computerized MC can replace visual MC for this purpose.

Within-subject electrocardiographic differences at equal heart rates: role of the autonomic nervous system.

Various combinations of sympathetic and vagal tone can yield the same heart rate, while ventricular electrophysiology differs. To demonstrate this in humans, we studied healthy volunteers in the sitting position with horizontal legs. First, heart rate was increased by lowering the legs to 60 degrees and back. Thereafter, heart rate was increased by handgrip. In each subject, a leg-lowering angle was selected at which heart rate matched best with heart rate in the third handgrip minute. Thirteen subjects had a heart rate match better than 1%. Heart rate (control: 65.2+/-9.0 bpm) increased to 72.1+/-8.7 (leg lowering) and to 72.1+/-8.8 (handgrip) bpm. QRS azimuth, QRS duration, maximal T vector, T azimuth, T elevation, ST duration, QRS-T angle and QT interval differed significantly (P<0.05) between leg lowering and handgrip (QT interval 418+/-15 versus 435+/-21 ms). Also, septal dispersion of repolarization, assessed as the time difference between the apex and the end of the T wave in the V2 and V3 leads, differed significantly (V2: 96.7+/-19.3 versus 110.0+/-23.3 ms, P<0.01; V3: 88.7+/-19.3 versus 97.3+/-23.3 ms; P<0.01). Hence, leg lowering and handgrip cause different ventricular depolarization and repolarization. The hypertensive handgrip manoeuvre entails a longer QT interval and probably an increased septal dispersion of repolarization.

QT dispersion as an attribute of T-loop morphology.

BACKGROUND: The suggestion that increased QT dispersion (QTD) is due to increased differences in local action potential durations within the myocardium is wanting. An alternative explanation was sought by relating QTD to vectorcardiographic T-loop morphology. METHODS AND RESULTS: The T loop is characterized by its amplitude and width (defined as the spatial angle between the mean vectors of the first and second halves of the loop). We reasoned that small, wide ("pathological") T loops produce larger QTD than large, narrow ("normal") loops. To quantify the relationship between QTD and T-loop morphology, we used a program for automated analysis of ECGs and a database of 1220 standard simultaneous 12-lead ECGs. For each ECG, QT durations, QTD, and T-loop parameters were computed. T-loop amplitude and width were dichotomized, with 250 microV (small versus large amplitudes) and 30 degrees (narrow versus wide loops) taken as thresholds. Over all 1220 ECGs, QTDs were smallest for large, narrow T loops (54.2+/-27.1 ms) and largest for small, wide loops (69. 5+/-33.5 ms; P<0.001). CONCLUSIONS: QTD is an attribute of T-loop morphology, as expressed by T-loop amplitude and width.

T axis as an indicator of risk of cardiac events in elderly people.

BACKGROUND: The T axis was postulated to be a general marker of repolarisation abnormality, indicative of subclinical myocardial damage. The aim of this investigation was to assess the prognostic importance of the T axis for fatal and non-fatal cardiac events, in a prospective cohort study of men and women aged 55 years and older. METHODS: 2352 men and 3429 women from the population-based Rotterdam Study took part in the study. Electrocardiograms were done, and T axes were categorised as normal, borderline, or abnormal. Data were analysed with Cox's proportional-hazards models; adjustment for age and sex was done where appropriate. FINDINGS: During 3-6 (mean 4) years of follow-up of the 5781 participants, 165 (2.9%) fatal and 192 (3.3%) non-fatal cardiac events occurred. Participants with an abnormal T axis (n=609) had an increased risk of cardiac death (hazard ratio 3.9 [95% CI 2.8-5.6]), sudden cardiac death (4.4 [2.6-7.4]), non-fatal cardiac events (2.7 [1.9-3.9]), and combined fatal or non-fatal cardiac events (3.2 [2.5-4.1]); p<0.001 for each. Additional adjustment for established cardiovascular risk factors resulted in lower, but still significant risk for all endpoints. The risk associated with an abnormal T axis was higher than those for any other cardiovascular risk factor. Additional subgroup analyses indicated that the risk of cardiac death was not substantially modified by age, sex, or history of myocardial infarction. INTERPRETATION: The T axis is a strong and independent risk indicator of fatal and non-fatal cardiac events in the elderly.

Reproducibility of computerized ECG measurements and coding in a nonhospitalized elderly population.

The standard 12-lead electrocardiogram (ECG) is used in many epidemiologic studies to diagnose and predict cardiovascular disease. In view of this, knowledge about the reproducibility of ECG measurements and coding is essential. Minute-to-minute, day-to-day, and year-to-year variability of ECG measurements, composite scores, and Minnesota Code classification were assessed by use of a computer program, in 101 nonhospitalized elderly men and women. Interval ECG measurements were more reproducible than amplitude measurements. The best reproducibility was found for the overall QTc interval (coefficient of variation 3.1%, 4.0%, and 5.2% for the minute-to-minute, day-to-day, and year-to-year groups, respectively) and the poorest was found for the Cardiac Infarction Injury Score (coefficient of variation 67.1%, 78.5%, and 94.3%, respectively). Minnesota Code discrepancies occurred in 16%, 19%, and 22% of the ECGs in the minute-to-minute, day-to-day, and year-to-year groups, respectively. Reproducibility within specific code categories was much better. Overall, variability tended to increase with time. In the routine setting, electrode positioning had relatively little effect on total variability.

Measurement error as a source of QT dispersion: a computerised analysis.

OBJECTIVE: To establish a general method to estimate the measuring error in QT dispersion (QTD) determination, and to assess this error using a computer program for automated measurement of QTD. SUBJECTS: Measurements were done on 1220 standard simultaneous 12 lead electrocardiograms. DESIGN: The computer program was validated against two observers on a random subset of 100 electrocardiograms. Simple laws of physics require that at least five of the six extremity leads have the same QT duration. This allows the direct assessment of the error in measuring QTD derived from five extremity leads (QTD5). It also enables ST-T amplitude dependent distributions of measurement error in determining QT duration to be established. These QT error distributions were then used to estimate the error in measuring QTD from all 12 leads (QTD12). MAIN OUTCOME MEASURES: Mean and standard deviation of error in measuring QTD duration, QTD5, and QTD12. RESULTS: Performance of the program was comparable to that of observers. Errors in measuring QT duration (measured QT minus reference QT) fell from a mean (SD) of 6.9 (17.1) ms for ST-T amplitudes < 50 microV to -1.4 (6.3) ms for amplitudes > 350 microV. Measurement errors of QTD5 and QTD12 were 20.4 (11.5) ms and 29.4 (14.9) ms. CONCLUSIONS: The fact that no QTD can exist between five of the six extremity leads provides a means of estimating QTD measurement error. Measuring error of QT duration is dependent on ST-T amplitude. QTD measurement error is large compared with typical QTD values reported.

T-loop morphology as a marker of cardiac events in the elderly.

ST-T wave changes of electrocardiographic (ECG) leads have long been recognized as predictors of future cardiac events, but they only imperfectly characterize T-loop morphology. Using vectorcardiographic (VCG) parameters, we investigated the predictive value of T-loop abnormality for fatal and nonfatal cardiac events in a prospective cohort study among 5,815 elderly. Separately, the predictive value of an easily obtainable T-loop parameter, the T axis, was also assessed. Measurements were determined by a computer program, using VCGs reconstructed from the standard 12-lead ECGs. During the 3 to 6 (mean 4) years of follow-up, 166 fatal and 193 nonfatal cardiac events occurred. Subjects with an abnormal T-loop morphology had increased risks for fatal cardiac events (hazard ratio 4.3; 95% CI 3.0-6.4) and nonfatal cardiac events (3.0; 1.9-4.8). Risks associated with an abnormal T axis alone were only slightly lower. Additional adjustment for established cardiovascular risk indicators resulted in lower, but still highly significant risks. Both T-loop and T-axis abnormalities appear to be strong, independent risk indicators of cardiac events in the elderly.

Diagnostic interpretation of electrocardiograms in population-based research: computer program research physicians, or cardiologists?

We assessed the performance of diagnostic electrocardiogram (ECG) interpretation by the computer program MEANS and by research physicians, compared to cardiologists, in a physician-based study. To establish a strategy for ECG interpretation in health surveys, we also studied the diagnostic capacity of three scenarios: use of the computer program alone (A), computer program and cardiologist (B), and computer program, research physician, and cardiologist (C). A stratified random sample of 381 ECGs was drawn from ECGs collected in the Rotterdam Study (n = 3057), which were interpreted both by a trained research physician using a form for structured clinical evaluation and by MEANS. All ECGs were interpreted independently by two cardiologists; if they disagreed (n = 175) the ECG was judged by a third cardiologist. Five ECG diagnoses were considered: anterior and inferior myocardial infarction (MI), left and right bundle branch block (LBBB and RBBB), and left ventricular hypertrophy (LVH). Overall, sensitivities and specificities of MEANS and the research physicians were high. The sensitivity of MEANS ranged from 73.8% to 92.9% and of the research physician ranged from 71.8% to 96.9%. The specificity of MEANS ranged from 97.5% to 99.8% and of the research physician from 96.3% to 99.6%. To diagnose LVH, LBBB, and RBBB, use of the computer program alone gives satisfactory results. Preferably, all positive findings of anterior and inferior MI by the program should be verified by a cardiologist. We conclude that diagnostic ECG interpretation by computer can be very helpful in population-based research, being at least as good as ECG interpretation by a trained research physician, but much more efficient and therefore less expensive.

Effect of electrode positioning on ECG interpretation by computer.

The aim of this study was to assess the variability in automated electrocardiogram (ECG) interpretation due to electrode positioning variations. Such variations were simulated by using a set of 746 body surface potential mappings from apparently healthy individuals and patients with myocardial infarction or left ventricular hypertrophy. Four types of electrode position changes were simulated, and the effect on ECG measurements and diagnostic classifications was determined by a computer program. At most 6% of the cases showed important changes in classification for longitudinal shifts. Transversal shifts causes less than 1.5% of important changes. An expert cardiologist, who analyzed a subset of 80 cases, agreed with the computer in 38 of 40 cases in which it made no change. In the 40 cases with large diagnostic changes, the cardiologist made no change in 18 cases. The effect of electrode position changes on ECG classification by an expert cardiologist was about half of the effect determined by computerized ECG classification. The effects on classification are significant; therefore, correct placement of chest electrodes remains mandatory.

Use of the standard 12-lead ECG to simulate electrode displacements.

Placement of the precordial electrodes for recording a 12-lead electrocardiogram (ECG) is subject to variation. Previous research has shown that displacement, especially in the longitudinal direction, can lead to changes in diagnosis. In practice, both the displacement and the effects of displacement on an individual ECG are unknown. To assess this effect for a given ECG, the authors developed a method to simulate ECGs at different displacements using only the recorded ECG. The material consisted of 746 body surface potential maps (BSPMs) containing 232 cases without abnormalities, 277 with myocardial infarction (MI), and 237 with left ventricular hypertrophy. By interpolating BSPMs, ECGs from closely spaced electrode positions could be derived. Taking electrode positioning errors that may be encountered in practice, 40 ECGs at different electrode displacements (displaced ECGs) for each BSPM were derived. Using half of the BSPMs, for each displacement, a transformation matrix that transforms the ECG at the standard 12-lead electrode positions (standard ECG) to the displaced ECG was determined. Using the other half of the BSPMs, each displaced ECG was compared with the ECG yielded by the corresponding transformation matrix (transformed ECG). For each comparison, the differences were assessed between the two sets of ECG signals and between the diagnostic computer classifications of the two sets. Signal differences were expressed as mean absolute amplitude differences over the QRS. Computer interpretation of MI and left ventricular hypertrophy was graded in five levels of certainty (no, consider, possible, probable, definite). For instance, for the largest longitudinal displacement studied of about one intercostal space, the 96th percentile mean absolute amplitude difference over the test set was 204 microV. The percentage of cases showing a change in MI classification of more than two certainty levels was 2.7% for this displacement. When comparing the standard ECG with the displaced ECG, these figures were 434 microV and 8.3%, respectively. It is concluded that ECGs from displaced electrodes can be well simulated by transforming the standard ECG, both for the ECG signal and diagnostic classifications.

Validation of a new computer program for Minnesota coding.

The Minnesota code (MC) is a classification system for electrocardiograms (ECGs) that is used for ECG coding in epidemiologic studies. As the MC measurement procedures and rules are complex, visual coding is time-consuming and error-prone. Automation should reduce measurement and coding errors. The authors developed an MC program, closely adhering to the MC regulations. To validate the program, a test set of 300 ECGs containing a wide variety of codable patterns was collected. The ECGs were coded independently by the program and by an experienced human reader. A reference code ("truth") was established by resolving disagreements through a consensus procedure. If the computer and human agreed, they were considered to be correct. Sensitivity and specificity were computed for each of the nine main code categories of the MC, both for the computer and for visual coding. The results show that the program is as good as or better than the human reader for sensitivity and specificity of all MC categories. Particularly noteworthy is the good program performance for arrhythmia coding. Most coding differences between the program and truth arise from small, borderline measurement differences in combination with the all-or-none character of the coding criteria. In conclusion, computerized Minnesota coding is a valuable alternative or supplement to visual coding.

A method to reduce the effect of electrode position variations on automated ECG interpretation.

To reduce the effect of electrode position variations on the diagnostic interpretation of an ECG, ECG and VCG interpretations were combined. The reduction was assessed by generating ECGs with displaced electrodes for a group of subjects using Body Surface Potential Maps (BSPMs). VCGs were reconstructed from the ECGs. The group consisted of normals, cases with myocardial infarction (MI), and with left ventricular hypertrophy (LVH). The effects of four types of electrode position changes were assessed for the diagnostic categories MI and LVH. The combined interpretation proved to be less sensitive to large changes than either the ECG or the VCG interpretation alone. The number of small changes increased for the combined interpretation. The combined interpretation showed higher agreement with a human expert than the ECG interpretation alone.

Interpolation of body surface potential maps.

The performance of four methods for interpolation of body surface potential maps (BSPMs) for different electrode grid densities was assessed. This study is part of a research project on the influence of the variability of 12-lead electrocardiograms on computer interpretation due to small electrode position changes. Interpolated BSPMs can be used to simulate this variability. The set of BSPMs studied, derived from a 117-electrode grid with relatively many electrodes on the left precordial part of the thorax, consisted of 232 cases without abnormalities, 277 with infarction, and 237 with left ventricular hypertrophy. The interpolation methods used were fast Fourier transforms, Chebyshev polynomials, linear functions, and cubic splines (CS). In the horizontal plane, a reference signal was first interpolated and, thereafter, resampled using 11 different sets of electrodes with the number of electrodes ranging from 18 down to 8. In the vertical direction, five grids with electrodes only on the front of the thorax and nine grids with electrodes on the front and back were examined. As a performance measure for interpolation, mean absolute error (MAE) was used: the absolute differences between the reference signal and the interpolated signal, averaged over the QRS on all maps. All methods showed deteriorating performance for decreasing grid density. In the horizontal direction, CS proved to be slightly superior to other methods for the left precordial electrodes for all but the densest grid (e.g., MAE = 22.8 microV vs MAE > 24.8 microV for a 12-electrode grid). For electrodes not in that area, CS performed the best as well (MAE = 16.1 microV for the same grid), with differences with the other methods being small (MAE > 16.4 microV). In the vertical direction, CS showed the best results on the front, both for the dense nonperiodic (MAE = 19.1 microV vs MAE > 26.6 microV for a 6-electrode grid) and periodic grids (MAE = 25.1 microV vs MAE > 26.6 microV for a 12-electrode grid). Linear functions performed best for sparse nonperiodic grids and sparse periodic grids for electrodes on the back, with the difference with CS for the last case being small. The method CS performed best overall, and is recommended for interpolating BSPMs.