Noninvasive mapping of transmural potentials during activation in swine hearts from body surface electrocardiograms

Systematic review of computational techniques, dataset utilization, and feature extraction in electrocardiographic imaging

This study aimed to analyze computational techniques in ECG imaging (ECGI) reconstruction, focusing on dataset identification, problem-solving, and feature extraction. We employed a PRISMA approach to review studies from Scopus and Web of Science, applying Cochrane principles to assess risk of bias. The selection was limited to English peer-reviewed papers published from 2010 to 2023, excluding studies that lacked computational technique descriptions. From 99 reviewed papers, trends show a preference for traditional methods like the boundary element and Tikhonov methods, alongside a rising use of advanced technologies including hybrid techniques and deep learning. These advancements have enhanced cardiac diagnosis and treatment precision. Our findings underscore the need for robust data utilization and innovative computational integration in ECGI, highlighting promising areas for future research and advances. This shift toward tailored cardiac care suggests significant progress in diagnostic and treatment methods.

Nonlocal based FISTA network for noninvasive cardiac transmembrane potential imaging

Objective. The primary aim of our study is to advance our understanding and diagnosis of cardiac diseases. We focus on the reconstruction of myocardial transmembrane potential (TMP) from body surface potential mapping.Approach. We introduce a novel methodology for the reconstruction of the dynamic distribution of TMP. This is achieved through the integration of convolutional neural networks with conventional optimization algorithms. Specifically, we utilize the subject-specific transfer matrix to describe the dynamic changes in TMP distribution and ECG observations at the body surface. To estimate the TMP distribution, we employ LNFISTA-Net, a learnable non-local regularized iterative shrinkage-thresholding network. The coupled estimation processes are iteratively repeated until convergence.Main results. Our experiments demonstrate the capabilities and benefits of this strategy. The results highlight the effectiveness of our approach in accurately estimating the TMP distribution, thereby providing a reliable method for the diagnosis of cardiac diseases.Significance. Our approach demonstrates promising results, highlighting its potential utility for a range of applications in the medical field. By providing a more accurate and dynamic reconstruction of TMP, our methodology could significantly improve the diagnosis and treatment of cardiac diseases, thereby contributing to advancements in healthcare.

Electrocardiographic imaging for cardiac arrhythmias and resynchronization therapy

Use of the 12-lead electrocardiogram (ECG) is fundamental for the assessment of heart disease, including arrhythmias, but cannot always reveal the underlying mechanism or the location of the arrhythmia origin. Electrocardiographic imaging (ECGi) is a non-invasive multi-lead ECG-type imaging tool that enhances conventional 12-lead ECG. Although it is an established technology, its continuous development has been shown to assist in arrhythmic activation mapping and provide insights into the mechanism of cardiac resynchronization therapy (CRT). This review addresses the validity, reliability, and overall feasibility of ECGi for use in a diverse range of arrhythmias. A systematic search limited to full-text human studies published in peer-reviewed journals was performed through Medline via PubMed, using various combinations of three key concepts: ECGi, arrhythmia, and CRT. A total of 456 studies were screened through titles and abstracts. Ultimately, 42 studies were included for literature review. Evidence to date suggests that ECGi can be used to provide diagnostic insights regarding the mechanistic basis of arrhythmias and the location of arrhythmia origin. Furthermore, ECGi can yield valuable information to guide therapeutic decision-making, including during CRT. Several studies have used ECGi as a diagnostic tool for atrial and ventricular arrhythmias. More recently, studies have tested the value of this technique in predicting outcomes of CRT. As a non-invasive method for assessing cardiovascular disease, particularly arrhythmias, ECGi represents a significant advancement over standard procedures in contemporary cardiology. Its full potential has yet to be fully explored.

Prediction of atrial fibrillation using the recurrence complex network of body surface potential mapping signals

OBJECTIVE:

Atrial fibrillation (AF) is the most common type of persistent arrhythmia. Early diagnosis and intervention of AF is essential to avert the further fatality. The technique of noninvasive electrical mapping, especially the body surface potential mapping (BSPM), has a more practical application in the study of predicting AF, when compared with the invasive electrical mapping methods such as the epicardial mapping and interventional catheter mapping. However, the prediction of AF with noninvasive signals has been inadequately studied. Thus, the aim of this paper was to analyze the properties of atrial dynamic system based on the noninvasive BSPM signals (BSPMs), using the recurrence complex network, and consequently to evaluate its role in predicting the recurrence of AF in clinical aspect.

METHOD:

Twelve patients with persistent AF were included in this study. Their preoperative and postoperative BSPMs were recorded. Initially, the preoperative BSPMs were transformed into the recurrence complex network to characterize the complexity property of the atria. Subsequently, the parameters of recurrence ratio (REC), determinism (DET), entropy of the diagonal structure distribution (ENTR), and laminarity (LAM) were calculated. Furthermore, the difference in the parameters in the four regions of the body and the difference obtained from the dominant frequency (DF) method were compared. Finally, the results obtained for the atrial dynamic system complexity from a 12-lead electrocardiogram (ECG) from the BSPMs were discussed.

RESULTS:

Our study revealed that the patients whose REC is greater than an average threshold, and with a lower LAM presented a much higher possibility of AF recurrence, after the AF surgery.

CONCLUSIONS:

The recurrence complex network is a useful and convenient way to evaluate the nonlinear properties of the BSPMs in patients with AF. It has good immunity to the lead position and has a potential role in the understanding of predicting the recurrence of AF.

Noninvasive Mapping and Electrocardiographic Imaging in Atrial and Ventricular Arrhythmias (CardioInsight)

Electrocardiographic imaging is a mapping technique aiming to noninvasively characterize cardiac electrical activity using signals collected from the torso to reconstruct epicardial potentials. Its efficacy has been demonstrated clinically, from mapping premature ventricular complexes and accessory pathways to of complex arrhythmias. Electrocardiographic imaging uses a standardized workflow. Signals should be checked manually to avoid automatic processing errors. Reentry is confirmed in the presence of local activation covering the arrhythmia cycle length. Focal breakthroughs demonstrate a QS pattern associated with centrifugal activation. Electrocardiographic imaging offers a unique opportunity to better understand the mechanism of cardiac arrhythmias and guide ablation.

Noninvasive Activation Imaging of Ventricular Arrhythmias by Spatial Gradient Sparse in Frequency Domain-Application to Mapping Reentrant Ventricular Tachycardia

The aim of this paper is to develop and evaluate a novel imaging method [spatial gradient sparse in frequency domain (SSF)] for the reconstruction of activation sequences of ventricular arrhythmia from noninvasive body surface potential map (BSPM) measurements. We formulated and solved the electrocardiographic inverse problem in the frequency domain, and the activation time was encoded in the phase information of the imaging solution. A cellular automaton heart model was used to generate focal ventricular tachycardia (VT). Different levels of Gaussian white noise were added to simulate noise-contaminated BSPM. The performance of SSF was compared with that of weighted minimum norm inverse solution. We also evaluated the method in a swine model with simultaneous intracardiac and body surface recordings. Four reentrant VTs were observed in pigs with myocardial infarction generated by left anterior descending artery occlusion. The imaged activation sequences of reentrant VTs were compared with those obtained from intracardiac electrograms. In focal VT simulation, SSF has increased the correlation coefficient (CC) by 5% and decreased localization errors (LEs) by 2.7 mm on average under different noise levels. In the animal validation with reentrant VT, SSF has achieved an average CC of 88% and an average LE of 6.3 mm in localizing the earliest and latest activation site in the reentry circuit. Our promising results suggest that the SSF provides noninvasive imaging capability of detecting and mapping macro-reentrant circuits in 3-D ventricular space. The SSF may become a useful imaging tool of identifying and localizing the potential targets for ablation of focal and reentrant VT.

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.

Localization of Origins of Premature Ventricular Contraction by Means of Convolutional Neural Network From 12-Lead ECG

OBJECTIVE

This paper proposes a novel method to localize origins of premature ventricular contractions (PVCs) from 12-lead electrocardiography (ECG) using convolutional neural network (CNN) and a realistic computer heart model.

METHODS

The proposed method consists of two CNNs (Segment CNN and Epi-Endo CNN) to classify among ventricular sources from 25 segments and from epicardium (Epi) or endocardium (Endo). The inputs are the full time courses and the first half of QRS complexes of 12-lead ECG, respectively. After registering the ventricle computer model with an individual patient's heart, the training datasets were generated by multiplying ventricular current dipoles derived from single pacing at various locations with patient-specific lead field. The origins of PVC are localized by calculating the weighted center of gravity of classification returned by the CNNs. A number of computer simulations were conducted to evaluate the proposed method under a variety of noise levels and heart registration errors. Furthermore, the proposed method was evaluated on 90 PVC beats from nine human patients with PVCs and compared against ablation outcome in the same patients.

RESULTS

The computer simulation evaluation returned relatively high accuracies for Segment CNN (∼78%) and Epi-Endo CNN (∼90%). Clinical testing in nine PVC patients resulted an averaged localization error of 11 mm.

CONCLUSION

Our simulation and clinical evaluation results demonstrate the capability and merits of the proposed CNN-based method for localization of PVC.

SIGNIFICANCE

This paper suggests a new approach for cardiac source localization of origin of arrhythmias using only the 12-lead ECG by means of CNN, and may have important applications for future real-time monitoring and localizing origins of cardiac arrhythmias guiding ablation treatment.

Activation recovery interval imaging of premature ventricular contraction

Dispersion of ventricular repolarization due to abnormal activation contributes to the susceptibility to cardiac arrhythmias. However, the global pattern of repolarization is difficult to assess clinically. Activation recovery interval (ARI) has been used to understand the properties of ventricular repolarization. In this study, we developed an ARI imaging technique to noninvasively reconstruct three-dimensional (3D) ARI maps in 10 premature ventricular contraction (PVC) patients and evaluated the results with the endocardial ARI maps recorded by a clinical navigation system (CARTO). From the analysis results of a total of 100 PVC beats in 10 patients, the average correlation coefficient is 0.86±0.05 and the average relative error is 0.06±0.03. The average localization error is 4.5±2.3 mm between the longest ARI sites in 3D ARI maps and those in CARTO endocardial ARI maps. The present results suggest that ARI imaging could serve as an alternative of evaluating global pattern of ventricular repolarization noninvasively and could assist in the future investigation of the relationship between global repolarization dispersion and the susceptibility to cardiac arrhythmias.

[Rhythm analysis of body surface potential mapping recordings from atrial fibrillation patients based on autocorrelation function]

The study of atrial fibrillation (AF) has been known as a hot topic of clinical concern. Body surface potential mapping (BSPM), a noninvasive electrical mapping technology, has been widely used in the study of AF. This study adopted 10 AF patients' preoperative and postoperative BSPM data (each patient's data contained 128 channels), and applied the autocorrelation function method to obtain the activation interval of the BSPM signals. The activation interval results were compared with that of manual counting method and the applicability of the autocorrelation function method was verified. Furthermore, we compared the autocorrelation function method with the commonly used fast Fourier transform (FFT) method. It was found that the autocorrelation function method was more accurate. Finally, to find a simple rule to predict the recurrence of atrial fibrillation, the autocorrelation function method was used to analyze the preoperative BSPM signals of 10 patients with persistent AF. Consequently, we found that if the patient's proportion of channels with dominant frequency larger than 2.5 Hz in the anterior left region is greater than the other three regions (the anterior right region, the posterior left region, and the posterior right region), he or she might have a higher possibility of AF recurrence. This study verified the rationality of the autocorrelation function method for rhythm analysis and concluded a simple rule of AF recurrence prediction based on this method.

Model-Based Feature Augmentation for Cardiac Ablation Target Learning From Images

GOAL

We present a model-based feature augmentation scheme to improve the performance of a learning algorithm for the detection of cardiac radio-frequency ablation (RFA) targets with respect to learning from images alone.

METHODS

Initially, we compute image features from delayed-enhanced magnetic resonance imaging (DE-MRI) to describe local tissue heterogeneities and feed them into a machine learning framework with uncertainty assessment for the identification of potential ablation targets. Next, we introduce the use of a patient-specific image-based model derived from DE-MRI coupled with the Mitchell-Schaeffer electrophysiology model and a dipole formulation for the simulation of intracardiac electrograms. Relevant features are extracted from these simulated signals which serve as a feature augmentation scheme for the learning algorithm. We assess the classifier's performance when using only image features and with model-based feature augmentation.

RESULTS

We obtained average classification scores of 97.2 % accuracy, 82.4 % sensitivity, and 95.0 % positive predictive value by using a model-based feature augmentation scheme. Preliminary results also show that training the algorithm on the closest patient from the database, instead of using all the patients, improves the classification results.

CONCLUSION

We presented a feature augmentation scheme based on biophysical cardiac electrophysiology modeling to increase the prediction scores of a machine learning framework for the RFA target prediction.

SIGNIFICANCE

The results derived from this study are a proof of concept that the use of model-based feature augmentation strengthens the performance of a purely image driven learning scheme for the prediction of cardiac ablation targets.

Noninvasive Personalization of a Cardiac Electrophysiology Model From Body Surface Potential Mapping

GOAL

We use noninvasive data (body surface potential mapping, BSPM) to personalize the main parameters of a cardiac electrophysiological (EP) model for predicting the response to different pacing conditions.

METHODS

First, an efficient forward model is proposed, coupling the Mitchell-Schaeffer transmembrane potential model with a current dipole formulation. Then, we estimate the main parameters of the cardiac model: activation onset location and tissue conductivity. A large patient-specific database of simulated BSPM is generated, from which specific features are extracted to train a machine learning algorithm. The activation onset location is computed from a Kernel Ridge Regression and a second regression calibrates the global ventricular conductivity.

RESULTS

The evaluation of the results is done both on a benchmark dataset of a patient with premature ventricular contraction (PVC) and on five nonischaemic implanted cardiac resynchonization therapy (CRT) patients with a total of 21 different pacing conditions. Good personalization results were found in terms of the activation onset location for the PVC (mean distance error, MDE = 20.3 mm), for the pacing sites (MDE = 21.7 mm) and for the CRT patients (MDE = 24.6 mm). We tested the predictive power of the personalized model for biventricular pacing and showed that we could predict the new electrical activity patterns with a good accuracy in terms of BSPM signals.

CONCLUSION

We have personalized the cardiac EP model and predicted new patient-specific pacing conditions.

SIGNIFICANCE

This is an encouraging first step towards a noninvasive preoperative prediction of the response to different pacing conditions to assist clinicians for CRT patient selection and therapy planning.

ECG imaging of ventricular tachycardia: evaluation against simultaneous non-contact mapping and CMR-derived grey zone

ECG imaging is an emerging technology for the reconstruction of cardiac electric activity from non-invasively measured body surface potential maps. In this case report, we present the first evaluation of transmurally imaged activation times against endocardially reconstructed isochrones for a case of sustained monomorphic ventricular tachycardia (VT). Computer models of the thorax and whole heart were produced from MR images. A recently published approach was applied to facilitate electrode localization in the catheter laboratory, which allows for the acquisition of body surface potential maps while performing non-contact mapping for the reconstruction of local activation times. ECG imaging was then realized using Tikhonov regularization with spatio-temporal smoothing as proposed by Huiskamp and Greensite and further with the spline-based approach by Erem et al. Activation times were computed from transmurally reconstructed transmembrane voltages. The results showed good qualitative agreement between the non-invasively and invasively reconstructed activation times. Also, low amplitudes in the imaged transmembrane voltages were found to correlate with volumes of scar and grey zone in delayed gadolinium enhancement cardiac MR. The study underlines the ability of ECG imaging to produce activation times of ventricular electric activity-and to represent effects of scar tissue in the imaged transmembrane voltages.

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.

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.

Data-driven estimation of cardiac electrical diffusivity from 12-lead ECG signals

Diagnosis and treatment of dilated cardiomyopathy (DCM) is challenging due to a large variety of causes and disease stages. Computational models of cardiac electrophysiology (EP) can be used to improve the assessment and prognosis of DCM, plan therapies and predict their outcome, but require personalization. In this work, we present a data-driven approach to estimate the electrical diffusivity parameter of an EP model from standard 12-lead electrocardiograms (ECG). An efficient forward model based on a mono-domain, phenomenological Lattice-Boltzmann model of cardiac EP, and a boundary element-based mapping of potentials to the body surface is employed. The electrical diffusivity of myocardium, left ventricle and right ventricle endocardium is then estimated using polynomial regression which takes as input the QRS duration and electrical axis. After validating the forward model, we computed 9500 EP simulations on 19 different DCM patients in just under three seconds each to learn the regression model. Using this database, we quantify the intrinsic uncertainty of electrical diffusion for given ECG features and show in a leave-one-patient-out cross-validation that the regression method is able to predict myocardium diffusion within the uncertainty range. Finally, our approach is tested on the 19 cases using their clinical ECG. 84% of them could be personalized using our method, yielding mean prediction errors of 18.7ms for the QRS duration and 6.5° for the electrical axis, both values being within clinical acceptability. By providing an estimate of diffusion parameters from readily available clinical data, our data-driven approach could therefore constitute a first calibration step toward a more complete personalization of cardiac EP.

Identifying model inaccuracies and solution uncertainties in noninvasive activation-based imaging of cardiac excitation using convex relaxation

Noninvasive imaging of cardiac electrical function has begun to move towards clinical adoption. Here, we consider one common formulation of the problem, in which the goal is to estimate the spatial distribution of electrical activation times during a cardiac cycle. We address the challenge of understanding the robustness and uncertainty of solutions to this formulation. This formulation poses a nonconvex, nonlinear least squares optimization problem. We show that it can be relaxed to be convex, at the cost of some degree of physiological realism of the solution set, and that this relaxation can be used as a framework to study model inaccuracy and solution uncertainty. We present two examples, one using data from a healthy human subject and the other synthesized with the ECGSIM software package. In the first case, we consider uncertainty in the initial guess and regularization parameter. In the second case, we mimic the presence of an ischemic zone in the heart in a way which violates a model assumption. We show that the convex relaxation allows understanding of spatial distribution of parameter sensitivity in the first case, and identification of model violation in the second.

Using transmural regularization and dynamic modeling for noninvasive cardiac potential imaging of endocardial pacing with imprecise thoracic geometry

Cardiac electrical imaging from body surface potential measurements is increasingly being seen as a technology with the potential for use in the clinic, for example for pre-procedure planning or during-treatment guidance for ventricular arrhythmia ablation procedures. However several important impediments to widespread adoption of this technology remain to be effectively overcome. Here we address two of these impediments: the difficulty of reconstructing electric potentials on the inner (endocardial) as well as outer (epicardial) surfaces of the ventricles, and the need for full anatomical imaging of the subject's thorax to build an accurate subject-specific geometry. We introduce two new features in our reconstruction algorithm: a nonlinear low-order dynamic parameterization derived from the measured body surface signals, and a technique to jointly regularize both surfaces. With these methodological innovations in combination, it is possible to reconstruct endocardial activation from clinically acquired measurements with an imprecise thorax geometry. In particular we test the method using body surface potentials acquired from three subjects during clinical procedures where the subjects' hearts were paced on their endocardia using a catheter device. Our geometric models were constructed using a set of CT scans limited in axial extent to the immediate region near the heart. The catheter system provides a reference location to which we compare our results. We compare our estimates of pacing site localization, in terms of both accuracy and stability, to those reported in a recent clinical publication , where a full set of CT scans were available and only epicardial potentials were reconstructed.

Computing ischemic regions in the heart with the bidomain model--first steps towards validation

We investigate whether it is possible to use the bidomain model and body surface potential maps (BSPMs) to compute the size and position of ischemic regions in the human heart. This leads to a severely ill posed inverse problem for a potential equation. We do not use the classical inverse problems of electrocardiography, in which the unknown sources are the epicardial potential distribution or the activation sequence. Instead we employ the bidomain theory to obtain a model that also enables identification of ischemic regions transmurally. This approach makes it possible to distinguish between subendocardial and transmural cases, only using the BSPM data. The main focus is on testing a previously published algorithm on clinical data, and the results are compared with images taken with perfusion scintigraphy. For the four patients involved in this study, the two modalities produce results that are rather similar: The relative differences between the center of mass and the size of the ischemic regions, suggested by the two modalities, are 10.8% ± 4.4% and 7.1% ± 4.6%, respectively. We also present some simulations which indicate that the methodology is robust with respect to uncertainties in important model parameters. However, in contrast to what has been observed in investigations only involving synthetic data, inequality constraints are needed to obtain sound results.