Temporal Sparse Promoting Three Dimensional Imaging of Cardiac Activation

Toward Enabling Cardiac Digital Twins of Myocardial Infarction Using Deep Computational Models for Inverse Inference

Cardiac digital twins (CDTs) have the potential to offer individualized evaluation of cardiac function in a non-invasive manner, making them a promising approach for personalized diagnosis and treatment planning of myocardial infarction (MI). The inference of accurate myocardial tissue properties is crucial in creating a reliable CDT of MI. In this work, we investigate the feasibility of inferring myocardial tissue properties from the electrocardiogram (ECG) within a CDT platform. The platform integrates multi-modal data, such as cardiac MRI and ECG, to enhance the accuracy and reliability of the inferred tissue properties. We perform a sensitivity analysis based on computer simulations, systematically exploring the effects of infarct location, size, degree of transmurality, and electrical activity alteration on the simulated QRS complex of ECG, to establish the limits of the approach. We subsequently present a novel deep computational model, comprising a dual-branch variational autoencoder and an inference model, to infer infarct location and distribution from the simulated QRS. The proposed model achieves mean Dice scores of 0.457 ±0.317 and 0.302 ±0.273 for the inference of left ventricle scars and border zone, respectively. The sensitivity analysis enhances our understanding of the complex relationship between infarct characteristics and electrophysiological features. The in silico experimental results show that the model can effectively capture the relationship for the inverse inference, with promising potential for clinical application in the future. The code is available at https://github.com/lileitech/MI_inverse_inference.

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.

Reducing Line-of-block Artifacts in Cardiac Activation Maps Estimated Using ECG Imaging: A Comparison of Source Models and Estimation Methods

OBJECTIVE

To investigate cardiac activation maps estimated using electrocardiographic imaging and to find methods reducing line-of-block (LoB) artifacts, while preserving real LoBs.

METHODS

Body surface potentials were computed for 137 simulated ventricular excitations. Subsequently, the inverse problem was solved to obtain extracellular potentials (EP) and transmembrane voltages (TMV). From these, activation times (AT) were estimated using four methods and compared to the ground truth. This process was evaluated with two cardiac mesh resolutions. Factors contributing to LoB artifacts were identified by analyzing the impact of spatial and temporal smoothing on the morphology of source signals.

RESULTS

AT estimation using a spatiotemporal derivative performed better than using a temporal derivative. Compared to deflection-based AT estimation, correlation-based methods were less prone to LoB artifacts but performed worse in identifying real LoBs. Temporal smoothing could eliminate artifacts for TMVs but not for EPs, which could be linked to their temporal morphology. TMVs led to more accurate ATs on the septum than EPs. Mesh resolution had a negligible effect on inverse reconstructions, but small distances were important for cross-correlation-based estimation of AT delays.

CONCLUSION

LoB artifacts are mainly caused by the inherent spatial smoothing effect of the inverse reconstruction. Among the configurations evaluated, only deflection-based AT estimation in combination with TMVs and strong temporal smoothing can prevent LoB artifacts, while preserving real LoBs.

SIGNIFICANCE

Regions of slow conduction are of considerable clinical interest and LoB artifacts observed in non-invasive ATs can lead to misinterpretations. We addressed this problem by identifying factors causing such artifacts and methods to reduce them.

ECG Localization Method Based on Volume Conductor Model and Kalman Filtering

The 12-lead electrocardiogram was invented more than 100 years ago and is still used as an essential tool in the early detection of heart disease. By estimating the time-varying source of the electrical activity from the potential changes, several types of heart disease can be noninvasively identified. However, most previous studies are based on signal processing, and thus an approach that includes physics modeling would be helpful for source localization problems. This study proposes a localization method for cardiac sources by combining an electrical analysis with a volume conductor model of the human body as a forward problem and a sparse reconstruction method as an inverse problem. Our formulation estimates not only the current source location but also the current direction. For a 12-lead electrocardiogram system, a sensitivity analysis of the localization to cardiac volume, tilted angle, and model inhomogeneity was evaluated. Finally, the estimated source location is corrected by Kalman filter, considering the estimated electrocardiogram source as time-sequence data. For a high signal-to-noise ratio (greater than 20 dB), the dominant error sources were the model inhomogeneity, which is mainly attributable to the high conductivity of the blood in the heart. The average localization error of the electric dipole sources in the heart was 12.6 mm, which is comparable to that in previous studies, where a less detailed anatomical structure was considered. A time-series source localization with Kalman filtering indicated that source mislocalization could be compensated, suggesting the effectiveness of the source estimation using the current direction and location simultaneously. For the electrocardiogram R-wave, the mean distance error was reduced to less than 7.3 mm using the proposed method. Considering the physical properties of the human body with Kalman filtering enables highly accurate estimation of the cardiac electric signal source location and direction. This proposal is also applicable to electrode configuration, such as ECG sensing systems.

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.

Noninvasive Imaging of Epicardial and Endocardial Potentials With Low Rank and Sparsity Constraints

In this study, we explore the use of low rank and sparse constraints for the noninvasive estimation of epicardial and endocardial extracellular potentials from body-surface electrocardiographic data to locate the focus of premature ventricular contractions (PVCs). The proposed strategy formulates the dynamic spatiotemporal distribution of cardiac potentials by means of low rank and sparse decomposition, where the low rank term represents the smooth background and the anomalous potentials are extracted in the sparse matrix. Compared to the most previous potential-based approaches, the proposed low rank and sparse constraints are batch spatiotemporal constraints that capture the underlying relationship of dynamic potentials. The resulting optimization problem is solved using alternating direction method of multipliers. Three sets of simulation experiments with eight different ventricular pacing sites demonstrate that the proposed model outperforms the existing Tikhonov regularization (zero-order, second-order) and L1-norm based method at accurately reconstructing the potentials and locating the ventricular pacing sites. Experiments on a total of 39 cases of real PVC data also validate the ability of the proposed method to correctly locate ectopic pacing sites.

Localization of Activation Origin on Patient-Specific Endocardial Surface by the Equivalent Double Layer (EDL) Source Model With Sparse Bayesian Learning

OBJECTIVE

Ablation treatment of ventricular arrhythmias can be facilitated by pre-procedure planning aided by electrocardiographic inverse solution, which can help to localize the origin of arrhythmia. Our aim was to improve localization accuracy of the inverse solution for activation originating on the left-ventricular endocardial surface, by using a sparse Bayesian learning (SBL).

METHODS

The inverse problem of electrocardiography was solved by reconstructing endocardial potentials from time integrals of body-surface electrocardiograms and from patient-specific geometry of the heart and torso for three patients with structurally normal ventricular myocardium, who underwent endocardial catheter mapping that included pace mapping. Complementary simulations using dipole sources in patient-specific geometry were also performed. The proposed method is using sparse property of the equivalent-double-layer (EDL) model of cardiac sources; it employs the SBL and makes use of the spatio-temporal features of the cardiac action potentials.

RESULTS

The mean localization error of the proposed method for pooled pacing sites ( n=52) was significantly smaller ( p=0.0039) than that achieved for the same patients in the study of Erem et al. Simulation experiments localized the source dipoles ( n=48) from forward-simulated potentials with the error of 9.4 ± 4.5 mm (mean ± SD).

CONCLUSION

The results of our clinical and simulation experiments demonstrate that localization of left-ventricular endocardial activation by means of the Bayesian approach, based on sparse representation of sources by EDL, is feasible and accurate.

SIGNIFICANCE

The proposed approach to localizing endocardial sources may have important applications in pre-procedure assessment of arrhythmias and in guiding their ablation treatment.

Localization of Activation Origin on Patient-Specific Epicardial Surface by Empirical Bayesian Method

OBJECTIVE

Ablation treatment of ventricular arrhythmias can be facilitated by pre-procedure planning aided by electrocardiographic inverse solution, which can help to localize the origin of arrhythmia. Our aim was to improve localization accuracy of the inverse solution by using a novel Bayesian approach.

METHODS

The inverse problem of electrocardiography was solved by reconstructing epicardial potentials from 120 body-surface electrocardiograms and from patient-specific geometry of the heart and torso for four patients suffering from scar-related ventricular tachycardia who underwent epicardial catheter mapping, which included pace-mapping. Simulations using dipole sources in patient-specific geometry were also performed. The proposed method, using dynamic spatio-temporal a priori constraints of the solution, was compared with classical Tikhonov methods based on fixed constraints.

RESULTS

The mean localization error of the proposed method for all available pacing sites (n=78) was significantly smaller than that achieved by Tikhonov methods; specifically, the localization accuracy for pacing in the normal tissue (n=17) was [Formula: see text] mm (mean ± SD) versus [Formula: see text] mm reported in the previous study using the same clinical data and Tikhonov regularization. Simulation experiments further supported these clinical findings.

CONCLUSION

The promising results of in vivo and in silico experiments presented in this study provide a strong incentive to pursuing further investigation of data-driven Bayesian methods in solving the electrocardiographic inverse problem.

SIGNIFICANCE

The proposed approach to localizing origin of ventricular activation sequence may have important applications in pre-procedure assessment of arrhythmias and in guiding their ablation treatment.

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.

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.

Transfer Learning From Simulations on a Reference Anatomy for ECGI in Personalized Cardiac Resynchronization Therapy

GOAL

Noninvasive cardiac electrophysiology (EP) model personalisation has raised interest for instance in the scope of predicting EP cardiac resynchronization therapy (CRT) response. However, the restricted clinical applicability of current methods is due in particular to the limitation to simple situations and the important computational cost.

METHODS

We propose in this manuscript an approach to tackle these two issues. First, we analyze more complex propagation patterns (multiple onsets and scar tissue) using relevance vector regression and shape dimensionality reduction on a large simulated database. Second, this learning is performed offline on a reference anatomy and transferred onto patient-specific anatomies in order to achieve fast personalized predictions online.

RESULTS

We evaluated our method on a dataset composed of 20 dyssynchrony patients with a total of 120 different cardiac cycles. The comparison with a commercially available electrocardiographic imaging (ECGI) method shows a good identification of the cardiac activation pattern. From the cardiac parameters estimated in sinus rhythm, we predicted five different paced patterns for each patient. The comparison with the body surface potential mappings (BSPM) measured during pacing and the ECGI method indicates a good predictive power.

CONCLUSION

We showed that learning offline from a large simulated database on a reference anatomy was able to capture the main cardiac EP characteristics from noninvasive measurements for fast patient-specific predictions.

SIGNIFICANCE

The fast CRT pacing predictions are a step forward to a noninvasive CRT patient selection and therapy optimisation, to help clinicians in these difficult tasks.

Wavelet-promoted sparsity for non-invasive reconstruction of electrical activity of the heart

We investigated a novel sparsity-based regularization method in the wavelet domain of the inverse problem of electrocardiography that aims at preserving the spatiotemporal characteristics of heart-surface potentials. In three normal, anesthetized dogs, electrodes were implanted around the epicardium and body-surface electrodes were attached to the torso. Potential recordings were obtained simultaneously on the body surface and on the epicardium. A CT scan was used to digitize a homogeneous geometry which consisted of the body-surface electrodes and the epicardial surface. A novel multitask elastic-net-based method was introduced to regularize the ill-posed inverse problem. The method simultaneously pursues a sparse wavelet representation in time-frequency and exploits correlations in space. Performance was assessed in terms of quality of reconstructed epicardial potentials, estimated activation and recovery time, and estimated locations of pacing, and compared with performance of Tikhonov zeroth-order regularization. Results in the wavelet domain obtained higher sparsity than those in the time domain. Epicardial potentials were non-invasively reconstructed with higher accuracy than with Tikhonov zeroth-order regularization (p < 0.05), and recovery times were improved (p < 0.05). No significant improvement was found in terms of activation times and localization of origin of pacing. Next to improved estimation of recovery isochrones, which is important when assessing substrate for cardiac arrhythmias, this novel technique opens potentially powerful opportunities for clinical application, by allowing to choose wavelet bases that are optimized for specific clinical questions. Graphical Abstract The inverse problem of electrocardiography is to reconstruct heart-surface potentials from recorded bodysurface electrocardiograms (ECGs) and a torso-heart geometry. However, it is ill-posed and solving it requires additional constraints for regularization. We introduce a regularization method that simultaneously pursues a sparse wavelet representation in time-frequency and exploits correlations in space. Our approach reconstructs epicardial (heart-surface) potentials with higher accuracy than common methods. It also improves the reconstruction of recovery isochrones, which is important when assessing substrate for cardiac arrhythmias. This novel technique opens potentially powerful opportunities for clinical application, by allowing to choose wavelet bases that are optimized for specific clinical questions.

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.

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.

Non-invasive imaging of ventricular activation during pacing and arrhythmia: Methods and validation

Cardiovascular disease continued to be a leading killer world widely. Each year, about 400,000 cases of sudden cardiac arrest are reported in the U.S. alone. Clinically, radio-frequency ablative procedure has become widely applied in the treatment of ventricular arrhythmia. Non-invasive approaches have been demonstrated to be able to provide important information on the arrhythmogenesis and potentially assist in the clinical practice. In this work, we develop and validate a novel temporal sparse based imaging method, Cardiac Electrical Sparse Imaging (CESI). Computer simulation and animal validation results demonstrate that the CESI approach is capable of imaging with improved accuracy and robustness by exploiting the temporal sparse property underlying cellular electrophysiology. Overall, a CC of 0.8, RE of 0.2 and LE (localization error) of 7 mm has been achieved on human realistic simulation and good accuracy has been observed in canine simultaneous mapping studies. Also, the technique maintains full temporal resolution (RRE <; 0.04) in terms of the activation sequence under various disturbances and in various pathologies such as premature ventricular complex and ventricular tachycardia. Our promising results indicate the excellent performance of noninvasive imaging of cardiac activation under various arrhythmias, and its potential for aiding clinical management of lethal ventricular arrhythmia.

Influence of Modeling Errors on the Initial Estimate for Nonlinear Myocardial Activation Times Imaging Calculated With Fastest Route Algorithm

Noninvasive reconstruction of cardiac electrical activity has a great potential to support clinical decision making, planning, and treatment. Recently, significant progress has been made in the estimation of the cardiac activation from body surface potential maps (BSPMs) using boundary element method (BEM) with the equivalent double layer (EDL) as a source model. In this formulation, noninvasive assessment of activation times results in a nonlinear optimization problem with an initial estimate calculated with the fastest route algorithm (FRA). Each FRA-simulated activation sequence is converted into the ECG. The best initialization is determined by the sequence providing the highest correlation between predicted and measured potentials. We quantitatively assess the effects of the forward modeling errors on the FRA-based initialization. We present three simulation setups to investigate the effects of volume conductor model simplifications, neglecting the cardiac anisotropy and geometrical errors on the localization of ectopic beats starting on the ventricular surface. For the analysis, 12-lead ECG and 99 electrodes BSPM system were used. The areas in the heart exposing the largest localization errors were volume conductor model and electrode configuration specific with an average error <10 mm. The results show the robustness of the FRA-based initialization with respect to the considered modeling errors.