Imaging and visualization of 3-D cardiac electric activity

Biological current source imaging method based on acoustoelectric effect: A systematic review

Neuroimaging can help reveal the spatial and temporal diversity of neural activity, which is of utmost importance for understanding the brain. However, conventional non-invasive neuroimaging methods do not have the advantage of high temporal and spatial resolution, which greatly hinders clinical and basic research. The acoustoelectric (AE) effect is a fundamental physical phenomenon based on the change of dielectric conductivity that has recently received much attention in the field of biomedical imaging. Based on the AE effect, a new imaging method for the biological current source has been proposed, combining the advantages of high temporal resolution of electrical measurements and high spatial resolution of focused ultrasound. This paper first describes the mechanism of the AE effect and the principle of the current source imaging method based on the AE effect. The second part summarizes the research progress of this current source imaging method in brain neurons, guided brain therapy, and heart. Finally, we discuss the problems and future directions of this biological current source imaging method. This review explores the relevant research literature and provides an informative reference for this potential non-invasive neuroimaging method.

Electrocardiographic Imaging for Atrial Fibrillation: A Perspective From Computer Models and Animal Experiments to Clinical Value

Electrocardiographic imaging (ECGI) is a technique to reconstruct non-invasively the electrical activity on the heart surface from body-surface potential recordings and geometric information of the torso and the heart. ECGI has shown scientific and clinical value when used to characterize and treat both atrial and ventricular arrhythmias. Regarding atrial fibrillation (AF), the characterization of the electrical propagation and the underlying substrate favoring AF is inherently more challenging than for ventricular arrhythmias, due to the progressive and heterogeneous nature of the disease and its manifestation, the small volume and wall thickness of the atria, and the relatively large role of microstructural abnormalities in AF. At the same time, ECGI has the advantage over other mapping technologies of allowing a global characterization of atrial electrical activity at every atrial beat and non-invasively. However, since ECGI is time-consuming and costly and the use of electrical mapping to guide AF ablation is still not fully established, the clinical value of ECGI for AF is still under assessment. Nonetheless, AF is known to be the manifestation of a complex interaction between electrical and structural abnormalities and therefore, true electro-anatomical-structural imaging may elucidate important key factors of AF development, progression, and treatment. Therefore, it is paramount to identify which clinical questions could be successfully addressed by ECGI when it comes to AF characterization and treatment, and which questions may be beyond its technical limitations. In this manuscript we review the questions that researchers have tried to address on the use of ECGI for AF characterization and treatment guidance (for example, localization of AF triggers and sustaining mechanisms), and we discuss the technological requirements and validation. We address experimental and clinical results, limitations, and future challenges for fruitful application of ECGI for AF understanding and management. We pay attention to existing techniques and clinical application, to computer models and (animal or human) experiments, to challenges of methodological and clinical validation. The overall objective of the study is to provide a consensus on valuable directions that ECGI research may take to provide future improvements in AF characterization and treatment guidance.

Validation and Opportunities of Electrocardiographic Imaging: From Technical Achievements to Clinical Applications

Electrocardiographic imaging (ECGI) reconstructs the electrical activity of the heart from a dense array of body-surface electrocardiograms and a patient-specific heart-torso geometry. Depending on how it is formulated, ECGI allows the reconstruction of the activation and recovery sequence of the heart, the origin of premature beats or tachycardia, the anchors/hotspots of re-entrant arrhythmias and other electrophysiological quantities of interest. Importantly, these quantities are directly and non-invasively reconstructed in a digitized model of the patient's three-dimensional heart, which has led to clinical interest in ECGI's ability to personalize diagnosis and guide therapy. Despite considerable development over the last decades, validation of ECGI is challenging. Firstly, results depend considerably on implementation choices, which are necessary to deal with ECGI's ill-posed character. Secondly, it is challenging to obtain (invasive) ground truth data of high quality. In this review, we discuss the current status of ECGI validation as well as the major challenges remaining for complete adoption of ECGI in clinical practice. Specifically, showing clinical benefit is essential for the adoption of ECGI. Such benefit may lie in patient outcome improvement, workflow improvement, or cost reduction. Future studies should focus on these aspects to achieve broad adoption of ECGI, but only after the technical challenges have been solved for that specific application/pathology. We propose 'best' practices for technical validation and highlight collaborative efforts recently organized in this field. Continued interaction between engineers, basic scientists, and physicians remains essential to find a hybrid between technical achievements, pathological mechanisms insights, and clinical benefit, to evolve this powerful technique toward a useful role in clinical practice.

Three-Dimensional Noninvasive Imaging of Ventricular Arrhythmias in Patients With Premature Ventricular Contractions

OBJECTIVE

Noninvasive imaging of cardiac electrical activity promises to provide important information regarding the underlying arrhythmic substrates for successful ablation intervention and further understanding of the mechanism of such lethal disease. The aim of this study is to evaluate the performance of a novel 3-D cardiac activation imaging technique to noninvasively localize and image origins of focal ventricular arrhythmias in patients undergoing radio frequency ablation.

METHODS

Preprocedural ECG gated contrast enhanced cardiac CT images and body surface potential maps were collected from 13 patients within a week prior to the ablation. The electrical activation images were estimated over the 3-D myocardium using a cardiac electric sparse imaging technique, and compared with CARTO activation maps and the ablation sites in the same patients.

RESULTS

Noninvasively-imaged activation sequences were consistent with the CARTO mapping results with an average correlation coefficient of 0.79, average relative error of 0.19, and average relative resolution error of 0.017. The imaged initiation sites of premature ventricular contractions (PVCs) were, on average, within 8 mm of the last successful ablation site and within 3 mm of the nearest ablation site.

CONCLUSION

The present results demonstrate the excellent performance of the 3-D cardiac activation imaging technique in imaging the activation sequence associated with PVC, and localizing the initial sites of focal ventricular arrhythmias in patients. These promising results suggest that the 3-D cardiac activation imaging technique may become a useful tool for aiding clinical diagnosis and management of ventricular arrhythmias.

Noninvasive Imaging of Human Atrial Activation during Atrial Flutter and Normal Rhythm from Body Surface Potential Maps

Background

Knowledge of atrial electrophysiological properties is crucial for clinical intervention of atrial arrhythmias and the investigation of the underlying mechanism. This study aims to evaluate the feasibility of a novel noninvasive cardiac electrical imaging technique in imaging bi-atrial activation sequences from body surface potential maps (BSPMs).

Methods

The study includes 7 subjects, with 3 atrial flutter patients, and 4 healthy subjects with normal atrial activations. The subject-specific heart-torso geometries were obtained from MRI/CT images. The equivalent current densities were reconstructed from 208-channel BSPMs by solving the inverse problem using individual heart-torso geometry models. The activation times were estimated from the time instant corresponding to the highest peak in the time course of the equivalent current densities. To evaluate the performance, a total of 32 cycles of atrial flutter were analyzed. The imaged activation maps obtained from single beats were compared with the average maps and the activation maps measured from CARTO, by using correlation coefficient (CC) and relative error (RE).

Results

The cardiac electrical imaging technique is capable of imaging both focal and reentrant activations. The imaged activation maps for normal atrial activations are consistent with findings from isolated human hearts. Activation maps for isthmus-dependent counterclockwise reentry were reconstructed on three patients with typical atrial flutter. The method was capable of imaging macro counterclockwise reentrant loop in the right atrium and showed inter-atria electrical conduction through coronary sinus. The imaged activation sequences obtained from single beats showed good correlation with both the average activation maps (CC = 0.91±0.03, RE = 0.29±0.05) and the clinical endocardial findings using CARTO (CC = 0.70±0.04, RE = 0.42±0.05).

Conclusions

The noninvasive cardiac electrical imaging technique is able to reconstruct complex atrial reentrant activations and focal activation patterns in good consistency with clinical electrophysiological mapping. It offers the potential to assist in radio-frequency ablation of atrial arrhythmia and help defining the underlying arrhythmic mechanism.

Noninvasive Imaging of High-Frequency Drivers and Reconstruction of Global Dominant Frequency Maps in Patients With Paroxysmal and Persistent Atrial Fibrillation

OBJECTIVE

Highest dominant-frequency (DF) drivers maintaining atrial fibrillation (AF) activities are effective ablation targets for restoring sinus rhythms in patients. This study aims to investigate whether AF drivers with highest activation rate can be noninvasively localized by means of a frequency-based cardiac electrical imaging (CEI) technique, which may aid in the planning of ablation strategy and the investigation of the underlying mechanisms of AF.

METHOD

A total of seven out of 13 patients were recorded with spontaneous paroxysmal or persistent AF and analyzed. The biatrial DF maps were reconstructed by coupling 5-s BSPM with CT-determined patient geometry. The CEI results were compared with ablation sites and DFs found from BSPMs.

RESULTS

CEI imaged left-to-right maximal frequency gradient (7.42 ± 0.66 Hz versus 5.85 ± 1.2 Hz, LA versus RA, p < 0.05) in paroxysmal AF patients. Patients with persistent AF were imaged with a loss of the intrachamber frequency gradient and a dispersion of the fast sources in both chambers. CEI was able to capture the AF behaviors, which were characterized by short-term stability, dynamic transition, and spatial repetition of the highest DF sites. The imaged highest DF sites were consistent with ablation sites in patients studied.

CONCLUSIONS

The frequency-based CEI allows localization of AF drivers with highest DF and characterization of the spatiotemporal frequency behaviors, suggesting the possibility for individualizing treatment strategy and advancing understanding of the underlying AF mechanisms.

SIGNIFICANCE

The establishment of noninvasive imaging techniques localizing AF drivers would facilitate management of this significant cardiac arrhythmia.

Examining the Impact of Prior Models in Transmural Electrophysiological Imaging: A Hierarchical Multiple-Model Bayesian Approach

Noninvasive cardiac electrophysiological (EP) imaging aims to mathematically reconstruct the spatiotemporal dynamics of cardiac sources from body-surface electrocardiographic (ECG) data. This ill-posed problem is often regularized by a fixed constraining model. However, a fixed-model approach enforces the source distribution to follow a pre-assumed structure that does not always match the varying spatiotemporal distribution of actual sources. To understand the model-data relation and examine the impact of prior models, we present a multiple-model approach for volumetric cardiac EP imaging where multiple prior models are included and automatically picked by the available ECG data. Multiple models are incorporated as an Lp-norm prior for sources, where p is an unknown hyperparameter with a prior uniform distribution. To examine how different combinations of models may be favored by different measurement data, the posterior distribution of cardiac sources and hyperparameter p is calculated using a Markov Chain Monte Carlo (MCMC) technique. The importance of multiple-model prior was assessed in two sets of synthetic and real-data experiments, compared to fixed-model priors (using Laplace and Gaussian priors). The results showed that the posterior combination of models (the posterior distribution of p) as determined by the ECG data differed substantially when reconstructing sources with different sizes and structures. While the use of fixed models is best suited in situations where the prior assumption fits the actual source structures, the use of an automatically adaptive set of models may have the ability to better address model-data mismatch and to provide consistent performance in reconstructing sources with different properties.

Temporal Sparse Promoting Three Dimensional Imaging of Cardiac Activation

A new Cardiac Electrical Sparse Imaging (CESI) technique is proposed to image cardiac activation throughout the three-dimensional myocardium from body surface electrocardiogram (ECG) with the aid of individualized heart-torso geometry. The sparse property of cardiac electrical activity in the time domain is utilized in the temporal sparse promoting inverse solution, one formulated to achieve higher spatial-temporal resolution, stronger robustness and thus enhanced capability in imaging cardiac electrical activity. Computer simulations were carried out to evaluate the performance of this imaging method under various circumstances. A total of 12 single site pacing and 7 dual sites pacing simulations with artificial and the hospital recorded sensor noise were used to evaluate the accuracy and stability of the proposed method. Simulations with modeling error on heart-torso geometry and electrode-torso registration were also performed to evaluate the robustness of the technique. In addition to the computer simulations, the CESI algorithm was further evaluated using experimental data in an animal model where the noninvasively imaged activation sequences were compared with those measured with simultaneous intracardiac mapping. All of the CESI results were compared with conventional weighted minimum norm solutions. The present results show that CESI can image with better accuracy, stability and stronger robustness in both simulated and experimental circumstances. In sum, we have proposed a novel method for cardiac activation imaging, and our results suggest that the CESI has enhanced performance, and offers the potential to image the cardiac activation and to assist in the clinical management of ventricular arrhythmias.

Sensitivity of Noninvasive Cardiac Electrophysiological Imaging to Variations in Personalized Anatomical Modeling

OBJECTIVE

Noninvasive cardiac electrophysiological (EP) imaging techniques rely on anatomically-detailed heart-torso models derived from high-quality tomographic images of individual subjects. However, anatomical modeling involves variations that lead to unresolved uncertainties in the outcome of EP imaging, bringing questions to the robustness of these methods in clinical practice. In this study, we design a systematic statistical approach to assess the sensitivity of EP imaging methods to the variations in personalized anatomical modeling.

METHODS

We first quantify the variations in personalized anatomical models by a novel application of statistical shape modeling. Given the statistical distribution of the variation in personalized anatomical models, we then employ unscented transform to determine the sensitivity of EP imaging outputs to the variation in input personalized anatomical modeling.

RESULTS

We test the feasibility of our proposed approach using two of the existing EP imaging methods: epicardial-based electrocardiographic imaging and transmural electrophysiological imaging. Both phantom and real-data experiments show that variations in personalized anatomical models have negligible impact on the outcome of EP imaging.

CONCLUSION

This study verifies the robustness of EP imaging methods to the errors in personalized anatomical modeling and suggests the possibility to simplify the process of anatomical modeling in future clinical practice.

SIGNIFICANCE

This study proposes a systematic statistical approach to quantify anatomical modeling variations and assess their impact on EP imaging, which can be extended to find a balance between the quality of personalized anatomical models and the accuracy of EP imaging that may improve the clinical feasibility of EP imaging.

Imaging cardiac activation sequence during ventricular tachycardia in a canine model of nonischemic heart failure

Noninvasive cardiac activation imaging of ventricular tachycardia (VT) is important in the clinical diagnosis and treatment of arrhythmias in heart failure (HF) patients. This study investigated the ability of the three-dimensional cardiac electrical imaging (3DCEI) technique for characterizing the activation patterns of spontaneously occurring and norepinephrine (NE)-induced VTs in a newly developed arrhythmogenic canine model of nonischemic HF. HF was induced by aortic insufficiency followed by aortic constriction in three canines. Up to 128 body-surface ECGs were measured simultaneously with bipolar recordings from up to 232 intramural sites in a closed-chest condition. Data analysis was performed on the spontaneously occurring VTs (n=4) and the NE-induced nonsustained VTs (n=8) in HF canines. Both spontaneously occurring and NE-induced nonsustained VTs initiated by a focal mechanism primarily from the subendocardium, but occasionally from the subepicardium of left ventricle. Most focal initiation sites were located at apex, right ventricular outflow tract, and left lateral wall. The NE-induced VTs were longer, more rapid, and had more focal sites than the spontaneously occurring VTs. Good correlation was obtained between imaged activation sequence and direct measurements (averaged correlation coefficient of ∼0.70 over 135 VT beats). The reconstructed initiation sites were ∼10 mm from measured initiation sites, suggesting good localization in such a large animal model with cardiac size similar to a human. Both spontaneously occurring and NE-induced nonsustained VTs had focal initiation in this canine model of nonischemic HF. 3DCEI is feasible to image the activation sequence and help define arrhythmia mechanism of nonischemic HF-associated VTs.

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.

Noninvasive transmural electrophysiological imaging based on minimization of total-variation functional

While tomographic imaging of cardiac structure and kinetics has improved substantially, electrophysiological mapping of the heart is still restricted to the surface with little or no depth information beneath. The progress in reconstructing 3-D action potential from surface voltage data has been hindered by the intrinsic ill-posedness of the problem and the lack of a unique solution in the absence of prior assumptions. In this work, we propose a novel adaption of the total-variation (TV) prior to exploit the unique spatial property of transmural action potential of being piecewise smooth with a steep boundary (gradient) separating depolarized and repolarized regions. We present a variational TV-prior instead of a common discrete TV-prior for improved robustness to mesh resolution, and solve the TV-minimization by a sequence of weighted, first-order L2-norm minimization. In a large set of phantom experiments, the proposed method is shown to outperform existing quadratic methods in preserving the steep gradient of action potential along the border of infarcts, as well as in capturing the disruption to the normal path of electrical wavefronts. Real-data experiments also further demonstrate the potential of the proposed method in revealing the location and shape of infarcts when quadratic methods fail to do so.

Hierarchical multiple-model Bayesian approach to transmural electrophysiological imaging

Noninvasive electrophysiological (EP) imaging of the heart aims to mathematically reconstruct the spatiotemporal dynamics of cardiac current sources from body-surface electrocardiographic (ECG) data. This ill-posed problem is often regularized by a fixed constraining model. However, this approach enforces the source distribution to follow a preassumed spatial structure that does not always match the varying spatiotemporal distribution of current sources. We propose a hierarchical Bayesian approach to transmural EP imaging that employs a continuous combination of multiple models, each reflecting a specific spatial property for current sources. Multiple models are incorporated as an Lp-norm prior for current sources, where p is an unknown hyperparameter with a prior probabilistic distribution. The current source estimation is obtained as an optimally weighted combination of solutions across all models, the weight being determined from the posterior distribution of p inferred from ECG data. The accuracy of our approach is assessed in a set of synthetic and real-data experiments on human heart-torso models. While the use of fixed models (L1- and L2-norm) only properly recovers sources with specific structures, our method delivers consistent performance in reconstructing sources with various extents and structures.

Noninvasive cardiac activation imaging of ventricular arrhythmias during drug-induced QT prolongation in the rabbit heart

BACKGROUND

Imaging myocardial activation from noninvasive body surface potentials promises to aid in both cardiovascular research and clinical medicine.

OBJECTIVE

To investigate the ability of a noninvasive 3-dimensional cardiac electrical imaging technique for characterizing the activation patterns of dynamically changing ventricular arrhythmias during drug-induced QT prolongation in rabbits.

METHODS

Simultaneous body surface potential mapping and 3-dimensional intracardiac mapping were performed in a closed-chest condition in 8 rabbits. Data analysis was performed on premature ventricular complexes, couplets, and torsades de pointes (TdP) induced during intravenous administration of clofilium and phenylephrine with combinations of various infusion rates.

RESULTS

The drug infusion led to a significant increase in the QT interval (from 175 ± 7 to 274 ± 31 ms) and rate-corrected QT interval (from 183 ± 5 to 262 ± 21 ms) during the first dose cycle. All the ectopic beats initiated by a focal activation pattern. The initial beat of TdPs arose at the focal site, whereas the subsequent beats were due to focal activity from different sites or 2 competing focal sites. The imaged results captured the dynamic shift of activation patterns and were in good correlation with the simultaneous measurements, with a correlation coefficient of 0.65 ± 0.02 averaged over 111 ectopic beats. Sites of initial activation were localized to be ~5 mm from the directly measured initiation sites.

CONCLUSIONS

The 3-dimensional cardiac electrical imaging technique could localize the origin of activation and image activation sequence of TdP during QT prolongation induced by clofilium and phenylephrine in rabbits. It offers the potential to noninvasively investigate the proarrhythmic effects of drug infusion and assess the mechanisms of arrhythmias on a beat-to-beat basis.

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.

Noninvasive reconstruction of the three-dimensional ventricular activation sequence during pacing and ventricular tachycardia in the rabbit heart

Ventricular arrhythmias represent one of leading causes for sudden cardiac death, a significant problem in public health. Noninvasive imaging of cardiac electric activities associated with ventricular arrhythmias plays an important role in better our understanding of the mechanisms and optimizing the treatment options. The present study aims to rigorously validate a novel three-dimensional (3-D) cardiac electrical imaging (3-DCEI) technique with the aid of 3-D intra-cardiac mapping during paced rhythm and ventricular tachycardia (VT) in the rabbit heart. Body surface potentials and intramural bipolar electrical recordings were simultaneously measured in a closed-chest condition in thirteen healthy rabbits. Single-site pacing and dual-site pacing were performed from ventricular walls and septum. VTs and premature ventricular complexes (PVCs) were induced by intravenous norepinephrine (NE). The non-invasively imaged activation sequence correlated well with invasively measured counterparts, with a correlation coefficient of 0.72 and a relative error of 0.30 averaged over all paced beats and NE-induced PVCs and VT beats. The averaged distance from imaged site of initial activation to measured site determined from intra-cardiac mapping was ∼5mm. These promising results suggest that 3-DCEI is feasible to non-invasively localize the origins and image activation sequence of focal ventricular arrhythmias.

Noninvasive reconstruction of the three-dimensional ventricular activation sequence during pacing and ventricular tachycardia in the canine heart

Single-beat imaging of myocardial activation promises to aid in both cardiovascular research and clinical medicine. In the present study we validate a three-dimensional (3D) cardiac electrical imaging (3DCEI) technique with the aid of simultaneous 3D intracardiac mapping to assess its capability to localize endocardial and epicardial initiation sites and image global activation sequences during pacing and ventricular tachycardia (VT) in the canine heart. Body surface potentials were measured simultaneously with bipolar electrical recordings in a closed-chest condition in healthy canines. Computed tomography images were obtained after the mapping study to construct realistic geometry models. Data analysis was performed on paced rhythms and VTs induced by norepinephrine (NE). The noninvasively reconstructed activation sequence was in good agreement with the simultaneous measurements from 3D cardiac mapping with a correlation coefficient of 0.74 ± 0.06, a relative error of 0.29 ± 0.05, and a root mean square error of 9 ± 3 ms averaged over 460 paced beats and 96 ectopic beats including premature ventricular complexes, couplets, and nonsustained monomorphic VTs and polymorphic VTs. Endocardial and epicardial origins of paced beats were successfully predicted in 72% and 86% of cases, respectively, during left ventricular pacing. The NE-induced ectopic beats initiated in the subendocardium by a focal mechanism. Sites of initial activation were estimated to be ∼7 mm from the measured initiation sites for both the paced beats and ectopic beats. For the polymorphic VTs, beat-to-beat dynamic shifts of initiation site and activation pattern were characterized by the reconstruction. The present results suggest that 3DCEI can noninvasively image the 3D activation sequence and localize the origin of activation of paced beats and NE-induced VTs in the canine heart with good accuracy. This 3DCEI technique offers the potential to aid interventional therapeutic procedures for treating ventricular arrhythmias arising from epicardial or endocardial sites and to noninvasively assess the mechanisms of these arrhythmias.

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.

A Kalman filter-based approach to reduce the effects of geometric errors and the measurement noise in the inverse ECG problem

In this article, we aimed to reduce the effects of geometric errors and measurement noise on the inverse problem of Electrocardiography (ECG) solutions. We used the Kalman filter to solve the inverse problem in terms of epicardial potential distributions. The geometric errors were introduced into the problem via wrong determination of the size and location of the heart in simulations. An error model, which is called the enhanced error model (EEM), was modified to be used in inverse problem of ECG to compensate for the geometric errors. In this model, the geometric errors are modeled as additive Gaussian noise and their noise variance is added to the measurement noise variance. The Kalman filter method includes a process noise component, whose variance should also be estimated along with the measurement noise. To estimate these two noise variances, two different algorithms were used: (1) an algorithm based on residuals, (2) expectation maximization algorithm. The results showed that it is important to use the correct noise variances to obtain accurate results. The geometric errors, if ignored in the inverse solution procedure, yielded incorrect epicardial potential distributions. However, even with a noise model as simple as the EEM, the solutions could be significantly improved.

Noninvasive imaging of three-dimensional cardiac activation sequence during pacing and ventricular tachycardia

BACKGROUND

Imaging cardiac excitation within ventricular myocardium is important in the treatment of cardiac arrhythmias and might help improve our understanding of arrhythmia mechanisms.

OBJECTIVE

This study sought to rigorously assess the imaging performance of a 3-dimensional (3D) cardiac electrical imaging (3DCEI) technique with the aid of 3D intracardiac mapping from up to 216 intramural sites during paced rhythm and norepinephrine (NE)-induced ventricular tachycardia (VT) in the rabbit heart.

METHODS

Body surface potentials and intramural bipolar electrical recordings were simultaneously measured in a closed-chest condition in 13 healthy rabbits. Single-site pacing and dual-site pacing were performed from ventricular walls and septum. VTs and premature ventricular complexes (PVCs) were induced by intravenous NE. Computed tomography images were obtained to construct geometry models.

RESULTS

The noninvasively imaged activation sequence correlated well with invasively measured counterpart, with a correlation coefficient of 0.72 ± 0.04, and a relative error of 0.30 ± 0.02 averaged over 520 paced beats as well as 73 NE-induced PVCs and VT beats. All PVCs and VT beats initiated in the subendocardium by a nonreentrant mechanism. The averaged distance from the imaged site of initial activation to the pacing site or site of arrhythmias determined from intracardiac mapping was ∼5 mm. For dual-site pacing, the double origins were identified when they were located at contralateral sides of ventricles or at the lateral wall and the apex.

CONCLUSION

3DCEI can noninvasively delineate important features of focal or multifocal ventricular excitation. It offers the potential to aid in localizing the origins and imaging activation sequences of ventricular arrhythmias, and to provide noninvasive assessment of the underlying arrhythmia mechanisms.

Electrocardiologic and related methods of non-invasive detection and risk stratification in myocardial ischemia: state of the art and perspectives

Background: Electrocardiographic methods still provide the bulk of cardiovascular diagnostics. Cardiac ischemia is associated with typical alterations in cardiac biosignals that have to be measured, analyzed by mathematical algorithms and allegorized for further clinical diagnostics. The fast growing fields of biomedical engineering and applied sciences are intensely focused on generating new approaches to cardiac biosignal analysis for diagnosis and risk stratification in myocardial ischemia.

Objectives: To present and review the state of the art in and new approaches to electrocardiologic methods for non-invasive detection and risk stratification in coronary artery disease (CAD) and myocardial ischemia; secondarily, to explore the future perspectives of these methods.

Methods: In follow-up to the Expert Discussion at the 2008 Workshop on "Biosignal Analysis" of the German Society of Biomedical Engineering in Potsdam, Germany, we comprehensively searched the pertinent literature and databases and compiled the results into this review. Then, we categorized the state-of-the-art methods and selected new approaches based on their applications in detection and risk stratification of myocardial ischemia. Finally, we compared the pros and cons of the methods and explored their future potentials for cardiology.

Results: Resting ECG, particularly suited for detecting ST-elevation myocardial infarctions, and exercise ECG, for the diagnosis of stable CAD, are state-of-the-art methods. New exercise-free methods for detecting stable CAD include cardiogoniometry (CGM); methods for detecting acute coronary syndrome without ST elevation are Body Surface Potential Mapping, functional imaging and CGM. Heart rate variability and blood pressure variability analyses, microvolt T-wave alternans and signal-averaged ECG mainly serve in detecting and stratifying the risk for lethal arrythmias in patients with myocardial ischemia or previous myocardial infarctions. Telemedicine and ambient-assisted living support the electrocardiological monitoring of at-risk patients.

Conclusions: There are many promising methods for the exercise-free, non-invasive detection of CAD and myocardial ischemia in the stable and acute phases. In the coming years, these new methods will help enhance state-of-the-art procedures in routine diagnostics. The future can expect that equally novel methods for risk stratification and telemedicine will transition into clinical routine.

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.

Noninvasive three-dimensional cardiac activation imaging on a rabbit model

Three-dimensional cardiac activation imaging (3-DCAI) aims at imaging the activation sequence throughout the 3-D myocardium. In the present study, the performance of 3-DCAI was validated through both in vivo animal experiments and computer simulations under a pacing protocol. The non-invasively imaged activation sequence from body surface potential maps (BSPMs) was quantitatively compared with the measured activation sequence obtained from the simultaneous intramural recording using a 3-D intra-cardiac mapping technique in a rabbit model. In addition, computer simulations were conducted to provide further assessment of the performance of the 3-DCAI algorithm in a realistic-geometry rabbit heart-torso model. The encouraging results suggest that 3-DCAI can non-invasively image the activation sequence and localize the origin of activation with good accuracy.

Cardiac source localization by means of a single moving dipole solution during endocardial pacing in an animal model

The accuracy of localizing the initiation site of cardiac activation by noninvasively estimating a single moving dipole (SMD) was investigated in a swine model. Body surface potential mapping (BSPM) and intracavitary noncontact mapping (NCM) were performed simultaneously during acute left ventricular (LV) endocardial pacing. For each animal, the boundary element model was constructed from preoperative magnetic resonance images (MRI). In each pacing study, the initiation site was localized by inversely estimating the location of an SMD from BSPM data. The results were compared with the precise pacing sites recorded by the NCM system. In total, four pacing sites from two pigs were analyzed, and the averaged source localization error was 16.8 +/- 2.3 mm. The present results indicate the potential of localizing focal cardiac events by estimating single moving dipole.

Noninvasive three-dimensional cardiac activation imaging from body surface potential maps: a computational and experimental study on a rabbit model

Three-dimensional (3-D) cardiac activation imaging (3-DCAI) is a recently developed technique that aims at imaging the activation sequence throughout the the ventricular myocardium. 3-DCAI entails the modeling and estimation of the cardiac equivalent current density (ECD) distribution from which the activation time at any myocardial site is determined as the time point with the peak amplitude of local ECD estimates. In this paper, we report, for the first time, an in vivo validation study assessing the feasibility of 3-DCAI in comparison with the 3-D intracardiac mapping, for a group of four healthy rabbits undergoing the ventricular pacing from various locations. During the experiments, the body surface potentials and the intramural bipolar electrical recordings were simultaneously measured in a closed-chest condition. The ventricular activation sequence noninvasively imaged from the body surface measurements by using 3-DCAI was generally in agreement with that obtained from the invasive intramural recordings. The quantitative comparison between them showed a root mean square (rms) error of 7.42 +/-0.61 ms, a relative error (RE) of 0.24 +/-0.03, and a localization error (LE) of 5.47 +/-1.57 mm. The experimental results were also consistent with our computer simulations conducted in well-controlled and realistic conditions. The present study suggest that 3-DCAI can noninvasively capture some important features of ventricular excitation (e.g., the activation origin and the activation sequence), and has the potential of becoming a useful imaging tool aiding cardiovascular research and clinical diagnosis of cardiac diseases.

Ultrasound current source density imaging

Surgery to correct severe heart arrhythmias usually requires detailed maps of the cardiac activation wave prior to ablation. The pinpoint electrical mapping procedure is laborious and limited by its spatial resolution (5-10 mm). We propose ultrasound current source density imaging (UCSDI), a direct 3-D imaging technique that potentially facilitates existing mapping procedures with superior spatial resolution. The technique is based on a pressure-induced change in resistivity known as the acoustoelectric (AE) effect, which is spatially confined to the ultrasound focus. AE-modulated voltage recordings are used to map and reconstruct current densities. In this preliminary study, we tested UCSDI under controlled conditions and compared it with conventional electrical mapping techniques. A 2-D dipole field was produced by a pair of electrodes in a bath of 0.9% NaCl solution. Boundary electrodes detected the AE signal while a 7.5-MHz focused ultrasound transducer was scanned across the bath. UCSDI located the current source and sink to within 1 mm of their actual positions. A future UCSDI system potentially provides real-time 3-D images of the cardiac activation wave coregistered with anatomical ultrasound and would greatly facilitate corrective procedures for heart abnormalities.

Noninvasive estimation of three-dimensional cardiac electrical activities from body surface potential maps

A noninvasive three-dimensional (3D) cardiac electrical imaging (3DCEI) approach, which can estimate the location of the initiation site (IS) of activation and the resultant 3D activation sequence (AS) from body surface potential maps (BSPMs), was validated in an intact large mammalian model (swine) during acute ventricular pacing. Body surface potential mapping and intracavitary noncontact mapping (NCM) were performed simultaneously during pacing from both right ventricular (RV) sites (intramural) and left ventricular (LV) sites (endocardial). Subsequent 3DCEI analyses were performed on the measured BSPMs. In total, 5 RV and 5 LV sites from control and heart failure animals were paced. The averaged localization error of the RV and LV sites were 7.0+/-1.1 mm and 6.6+/-1.9 mm, respectively. The endocardial ASs as a subset of the estimated 3D ASs by 3DCEI were consistent with those reconstructed from the NCM system. The present experimental results demonstrate that the noninvasive 3DCEI approach can localize the initiation site and estimate cardiac activation sequence with good accuracy in an in vivo setting, under control, paced and/or diseased conditions.

Estimation of global ventricular activation sequences by noninvasive three-dimensional electrical imaging: validation studies in a Swine model during pacing

BACKGROUND

A novel noninvasive imaging technique, the heart-model-based three-dimensional cardiac electrical imaging (3DCEI) approach was previously developed and validated to estimate the initiation site (IS) of cardiac activity and the activation sequence (AS) from body surface potential maps (BSPMs) in a rabbit model. The aim of the present study was to validate the 3DCEI in an intact large mammalian model (swine) during acute ventricular pacing.

METHODS AND RESULTS

The heart-torso geometries were constructed from preoperative magnetic resonance (MR) images acquired from each animal. Body surface potential mapping and intracavitary noncontact mapping (NCM) were performed simultaneously during pacing from both right ventricular (RV) (intramural) and left ventricular (LV) sites (endocardial). Subsequent 3DCEI analyses were performed from the measured BSPMs. The estimated ISs were compared with the precise pacing locations and estimated ASs were compared with those recorded by the NCM system. In total, five RV and five LV sites from control and heart failure (HF) animals were paced and sequences of 100 paced beats were analyzed (10 for each site). The averaged localization error (LE) of the RV and LV sites were 7.3 +/- 1.8 mm (n = 50) and 7.0 +/- 2.2 mm (n = 50), respectively. The global 3D ASs throughout the ventricular myocardium were also derived. The endocardial ASs as a subset of the estimated 3D ASs were consistent with those reconstructed from the NCM system.

CONCLUSION

The present experimental results demonstrate that the noninvasive 3DCEI approach can localize the IS and estimate AS with good accuracy in an in vivo setting under control, paced, and/or diseased conditions.

Three-Dimensional Cardiac Electrical Imaging From Intracavity Recordings

A novel approach is proposed to image 3-D cardiac electrical activity from intracavity electrical recordings with the aid of a catheter. The feasibility and performance were evaluated by computer simulation studies, where a 3-D cellular-automaton heart model and a finite-element thorax volume conductor model were utilized. The finite-element method (FEM) was used to simulate the intracavity recordings induced by a single-site and dual-site pacing protocol. The 3-D ventricular activation sequences as well as the locations of the initial activation sites were inversely estimated by minimizing the dissimilarity between the intracavity potential "measurements" and the model-generated intracavity potentials. Under single-site pacing, the relative error (RE) between the true and estimated activation sequences was 0.03 +/- 0.01 and the localization error (LE) (of the initiation site) was 1.88 +/- 0.92 mm, as averaged over 12 pacing trials when considering 25 microV additive measurement noise using 64 catheter electrodes. Under dual-site pacing, the RE was 0.04 +/- 0.01 over 12 pacing trials and the LE over 24 initial pacing sites was 2.28 +/- 1.15 mm, when considering 25 microV additive measurement noise using 64 catheter electrodes. The proposed 3-D cardiac electrical imaging approach using intracavity electrical recordings was also tested under various simulated conditions and robust inverse solutions obtained. The present promising simulation results suggest the feasibility of obtaining 3-D information of cardiac electrical activity from intracavity recordings. The application of this inverse method has the potential of enhancing electrocardiographic mapping by catheters in electrophysiology laboratories, aiding cardiac resynchronization therapy, and other clinical applications.

Effects of cardiac anisotropy on three-dimensional inverse activation time imaging: a cellular-automaton-based model study

The aim of the present study was to examine the effect of cardiac anisotropy on the ECG activation time imaging solutions using a cellular-automaton-heart-model-based approach. The anisotropic heart model was used in forward pacing simulation while anisotropic and isotropic heart models were used in inverse procedures, respectively. Twenty-four sites throughout the ventricles were paced and the corresponding cardiac activation sequences were inversely estimated. The inverse solutions of activation time using anisotropic or isotropic heart model were compared. The computer simulation results suggest that the cardiac anisotropy does not seem to have significant influence on the inverse estimation of activation time using the present 3D heart-model- based inverse imaging approach during single paced activation.

Three-dimensional electrocardiographic imaging

We review our recent work on the development and evaluation of three-dimensional electrocardiographic imaging technology (3DEIT). Cardiac electrophysiological properties, including activation time and transmembrane potentials, are estimated from body surface ECG signals with the aid of a realistic geometry heart model in which electrophysiological a priori information is incorporated. We have conducted computer simulation studies to demonstrate the feasibility of imaging activation sequence throughout the three-dimensional myocardium, and localizing sites of origin of activation and arrhythmias using our 3DEIT approach. We will also review the pilot experimental studies in evaluating the 3DEIT approach in a patient with pacemaker and experimental animals with intracardiac recordings. Our promising results to date suggest the feasibility of the 3DEIT technology that we are developing and that it merits further investigation.

Three-dimensional activation sequence imaging in a rabbit model

This paper evaluates a biophysical-model based three-dimensional (3-D) activation sequence imaging approach in a rabbit model. In this approach, cardiac electrical sources within the myocardial volume are represented by distributed equivalent current densities; a realistic heart-torso volume conductor model is built from the CT scans of the rabbit's torso; spatial-temporal regularization is applied when solving the inverse problem of current density estimation; and the activation time at every myocardial location is determined as the time point when the estimated local current density reaches its maximum amplitude. Computer simulations have been conducted to image the activation sequence initiated by pacing 11 sites throughout the ventricular myocardium. Under 20muV Gaussian white noise, the average correlation coefficient (CC) between the imaged and the simulated activation sequences is 0.92, the average relative error (RE) is 0.19, and the average localization error (LE) is 4.99mm averaged over 11 pacing sites. Even under 60muV Gaussian white noise, reasonable results can still be achieved by the present approach with CC = 0.89, RE = 0.22, and LE = 6.85mm. The simulation results demonstrate that the present 3-D imaging approach has reasonable accuracy and robustness against recording noises.

Noninvasive reconstruction of three-dimensional ventricular activation sequence from the inverse solution of distributed equivalent current density

We propose a new electrocardiographic (ECG) inverse approach for imaging the three-dimensional (3-D) ventricular activation sequence based on the modeling and estimation of the equivalent current density throughout the entire volume of the ventricular myocardium. The spatio-temporal coherence of the ventricular excitation process has been utilized to derive the activation time from the estimated time course of the equivalent current density. In the present study, we explored four different linear inverse algorithms (the minimum norm and weighted minimum norm estimates in combination with two regularization schemes: the instant-by-instant regularization and the isotropy method) to estimate the current density at each time instant during the ventricular depolarization. The activation time at any given location within the ventricular myocardium was determined as the time point with the occurrence of the maximum local current density estimate. Computer simulations were performed to evaluate this approach using single- and dual-site pacing protocols in a physiologically realistic cellular automaton heart model. The performance and stability of the proposed approach was evaluated with respect to the various levels of measurement noise (0, 5, 10, 20, 40, and 60 microV), the various numbers of ECG electrodes and the modeling errors on the torso geometry and heart position. The simulation results demonstrate that: 1) the single-site paced 3-D activation sequence can be well reconstructed from 200-channel body surface potential maps with additive Gaussian white noise of 20 microV (correlation coefficient = 0.90, relative error = 0.19, and localization error = 5.49 mm); 2) a higher imaging accuracy can be obtained when the activation is initiated from the left/right ventricle (LV/RV) compared to from the septum; 3) the isotropy method gives rise to a better performance than the conventional instant-by-instant regularization; 4) a decreased imaging accuracy results from a larger noise level, a fewer number of electrodes, or the volume conductor modeling errors; however, a reasonable imaging accuracy can still be obtained with a 60 microV noise level, 64 electrodes, or mild errors on both the torso geometry and heart position, respectively; 5) the dual-site paced 3-D activation sequence can be imaged when the two sites are paced either simultaneously or with a time delay of 20 ms; 6) two pacing sites can be resolved and localized in the imaged 3-D activation sequence when they are located at the contralateral sides of ventricles or at the ventricular lateral wall and the apex, respectively.

Three-dimensional myocardial activation imaging in a rabbit model

This study evaluated a recently developed three-dimensional (3-D) electrocardiographic imaging (3DECI) approach in a closed-chest rabbit model. First, the performance and sensitivity of parameters of 3DECI were evaluated using a geometrically realistic rabbit heart-torso model. Second, a 3-D intracardiac mapping procedure was evaluated based on the heart-torso rabbit model. Third, comparisons were made among the forward-simulated ventricular activation sequence, the activation sequence derived by the intracardiac mapping, and the 3DECI inverse solution. Finally, the present procedures were tested in vivo in a rabbit, in which the relative error between the measured and imaged activation sequence was 0.20 and the localization error was 5.1 mm. The present simulation and experimental results suggest the merits of the 3DECI imaging approach, and the validity of intracardiac mapping as a tool to evaluate the 3DECI.

Noninvasive three-dimensional electrocardiographic imaging of ventricular activation sequence

Imaging the myocardial activation sequence is critical for improved diagnosis and treatment of life-threatening cardiac arrhythmias. It is desirable to reveal the underlying cardiac electrical activity throughout the three-dimensional (3-D) myocardium (rather than just the endocardial or epicardial surface) from noninvasive body surface potential measurements. A new 3-D electrocardiographic imaging technique (3-DEIT) based on the boundary element method (BEM) and multiobjective nonlinear optimization has been applied to reconstruct the cardiac activation sequences from body surface potential maps. Ultrafast computerized tomography scanning was performed for subsequent construction of the torso and heart models. Experimental studies were then conducted, during left and right ventricular pacing, in which noninvasive assessment of ventricular activation sequence by means of 3-DEIT was performed simultaneously with 3-D intracardiac mapping (up to 200 intramural sites) using specially designed plunge-needle electrodes in closed-chest rabbits. Estimated activation sequences from 3-DEIT were in good agreement with those constructed from simultaneously recorded intracardiac electrograms in the same animals. Averaged over 100 paced beats (from a total of 10 pacing sites), total activation times were comparable (53.3 +/- 8.1 vs. 49.8 +/- 5.2 ms), the localization error of site of initiation of activation was 5.73 +/- 1.77 mm, and the relative error between the estimated and measured activation sequences was 0.32 +/- 0.06. The present experimental results demonstrate that the 3-D paced ventricular activation sequence can be reconstructed by using noninvasive multisite body surface electrocardiographic measurements and imaging of heart-torso geometry. This new 3-D electrocardiographic imaging modality has the potential to guide catheter-based ablative interventions for the treatment of life-threatening cardiac arrhythmias.

Localization of the site of origin of reentrant arrhythmia from body surface potential maps: a model study

We have developed a model-based imaging approach to estimate the site of origin of reentrant arrhythmia from body surface potential maps (BSPMs), with the aid of a cardiac arrhythmia model. The reentry was successfully simulated and maintained in the cardiac model, and the simulated ECG waveforms over the body surface corresponding to a maintained reentry have evident characteristics of ventricular tachycardia. The performance of the inverse imaging approach was evaluated by computer simulations. The present simulation results show that an averaged localization error of about 1.5 mm, when 5% Gaussian white noise was added to the BSPMs, was detected. The effects of the heart-torso geometry uncertainty on the localization were also initially assessed and the simulation results suggest that no significant influence was observed when 10% torso geometry uncertainty or 10 mm heart position shifting was considered. The present simulation study suggests the feasibility of localizing the site of origin of reentrant arrhythmia from non-invasive BSPMs, with the aid of a cardiac arrhythmia model.

Non-invasive estimation of myocardial infarction by means of a heart-model-based imaging approach: A simulation study

In the study, a new myocardial infarction (MI) estimation method was developed for estimating MI in the three-dimensional myocardium by means of a heart-model-based inverse approach. The site and size of MI are estimated from body surface electrocardiograms by minimising multiple objective functions of the measured body surface potential maps (BSPMs) and the heart-model-generated BSPMs. Computer simulations were conducted to evaluate the performance of the developed method, using a single-site MI and dual-site MI protocols. The simulation results show that, for the single-site MI, the averaged spatial distance (SD) between the weighting centres of the 'true' and estimated MIs, and the averaged relative error (RE) between the numbers of the 'true' and estimated infarcted units are 3.0 +/- 0.6/3.6 +/- 0.6 mm and 0.11 +/- 0.02/0.14 +/- 0.02, respectively, when 5 microV/10 microV Gaussian white noise was added to the body surface potentials. For the dual-site MI, the averaged SD between the weighting centres of the 'true' and estimated MIs, and the averaged RE between the numbers of the 'true' and estimated infarcted units are 3.8 +/- 0.7/3.9 +/- 0.7mm and 0.12 +/- 0.02/0.14 +/- 0.03, respectively, when 5 microV/10 microV Gaussian white noise was added to the body surface potentials. The simulation results suggest the feasibility of applying the heart-model-based imaging approach to the estimation of myocardial infarction from body surface potentials.

Noninvasive imaging of cardiac transmembrane potentials within three-dimensional myocardium by means of a realistic geometry anisotropic heart model

We have developed a new approach for imaging cardiac transmembrane potentials (TMPs) within the three-dimensional (3-D) myocardium by means of an anisotropic heart model. The cardiac TMP distribution is estimated from body surface electrocardiograms by minimizing objective functions of the "measured" body surface potential maps (BSPMs) and the heart-model-generated BSPMs. Computer simulation studies have been conducted to evaluate the present 3-D TMP imaging approach using pacing protocols. Simulations of single-site pacing at 24 sites throughout the ventricles, as well as dual-site pacing at 12 pairs of sites in the vicinity of atrio-ventricular ring were performed. The present simulation results show that the correlation coefficient (CC) and relative error (RE) between the "true" and inversely estimated TMP distributions were 0.9915 +/- 0.0041 and 0.1266 +/- 0.0326, for single-site pacing, and 0.9889 +/- 0.0034 and 0.1473 +/- 0.0237 for dual-site pacing, respectively, when 10 microV Gaussian white noise (GWN) was added to the BSPMs. The effects of heart and torso geometry uncertainty were also evaluated by shifting the heart position by 10 mm and altering the torso size by 10%. The CC between the "true" and inversely estimated TMP distributions was above 0.97 when these geometry uncertainties were considered. The present simulation results demonstrate the feasibility of noninvasive estimation of TMP distribution throughout the ventricles from body surface electrocardiographic measurements, and suggest that the present method may become a useful alternative in noninvasive imaging of distributed cardiac electrophysiological processes within the 3-D myocardium.

Noninvasive localization of the site of origin of paced cardiac activation in human by means of a 3-D heart model

A recently developed heart-model-based localization approach is experimentally evaluated in noninvasively localizing the site of origin of cardiac activation in a patient with a pacemaker. The heart-torso model of the patient was constructed from the contrast ultrafast computed tomography images. The site of initial paced activation in the patient was quantitatively localized and compared with the tip position of the pacemaker lead. The localization error of the inverse estimation was found to be 5.2 mm with respect to the true lead tip position. The promising result of this pilot experimental study suggests the feasibility of localizing the site of origin of cardiac activation in an experimental setting. The heart-model-based localization approach may become an alternative tool in localizing the site of origin of cardiac activation.

Body surface Laplacian electrocardiogram of ventricular depolarization in normal human subjects

INTRODUCTION

The body surface Laplacian electrocardiogram (ECG) mapping provides a noninvasive means for spatiotemporal mapping of cardiac electrical events. The aim of the present study was to explore the relationship between the Laplacian ECG and the underlying cardiac activities during ventricular depolarization in healthy human subjects.

METHODS AND RESULTS

A 95-channel body surface potential ECG was recorded over the anterolateral chest from 11 healthy male subjects. The surface Laplacian (SL) ECG was estimated from the recorded potentials during QRS complex by means of a novel spline SL estimator, as well as by the conventional 5-point SL estimator for comparison purpose. A simulation study was also conducted using a realistic geometry heart-torso model in an attempt to qualitatively interpret the experimental results. For all subjects, more spatial details were observed in the SL ECG maps compared with the potential ECG maps, with spline SL more robust against noise than the 5-point SL. In total, three positive activities (denoted as P1, P2, P3) and four negative activities (denoted as N1, N2, N3, N4) in the spline SL ECG maps were observed during ventricular depolarization. Initial localized P1 and N1 activities were observed in 11 and 8 subjects, respectively. Then, the initial P1 was divided into three positive activities (P1, P2, P3) in 9 subjects. After the appearance of multiple positive activities, three negative activities (N2, N3, N4) appeared in 11, 8, and 9 subjects, respectively. Similar findings were obtained in the computer simulation study.

CONCLUSION

The present study demonstrates that the SL ECG provides more spatial details than the potential ECG, and multiple simultaneously active ventricular activities could be revealed in the SL ECG maps. The results suggest that the SL ECG may provide an alternative for noninvasive mapping of cardiac electrical activity.

Noninvasive three-dimensional activation time imaging of ventricular excitation by means of a heart-excitation model

We propose a new method for imaging activation time within three-dimensional (3D) myocardium by means of a heart-excitation model. The activation time is estimated from body surface electrocardiograms by minimizing multiple objective functions of the measured body surface potential maps (BSPMs) and the heart-model-generated BSPMs. Computer simulation studies have been conducted to evaluate the proposed 3D myocardial activation time imaging approach. Single-site pacing at 24 sites throughout the ventricles, as well as dual-site pacing at 12 pairs of sites in the vicinity of atrioventricular ring, was performed. The present simulation results show that the average correlation coefficient (CC) and relative error (RE) for single-site pacing were 0.9992+/-0.0008/0.9989+/-0.0008 and 0.05+/-0.02/0.07+/-0.03, respectively, when 5 microV/10 microV Gaussian white noise (GWN) was added to the body surface potentials. The average CC and RE for dual-site pacing were 0.9975+/-0.0037 and 0.08+/-0.04, respectively, when 10 microV GWN was added to the body surface potentials. The present simulation results suggest the feasibility of noninvasive estimation of activation time throughout the ventricles from body surface potential measurement, and suggest that the proposed method may become an important alternative in imaging cardiac electrical activity noninvasively.