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

Solving the Inverse Problem of Electrocardiography for Cardiac Digital Twins: A Survey

Cardiac digital twins (CDTs) are personalized virtual representations used to understand complex cardiac mechanisms. A critical component of CDT development is solving the ECG inverse problem, which enables the reconstruction of cardiac sources and the estimation of patient-specific electrophysiology (EP) parameters from surface ECG data. Despite challenges from complex cardiac anatomy, noisy ECG data, and the ill-posed nature of the inverse problem, recent advances in computational methods have greatly improved the accuracy and efficiency of ECG inverse inference, strengthening the fidelity of CDTs. This paper aims to provide a comprehensive review of the methods for solving ECG inverse problems, their validation strategies, their clinical applications, and their future perspectives. For the methodologies, we broadly classify state-of-the-art approaches into two categories: deterministic and probabilistic methods, including both conventional and deep learning-based techniques. Integrating physics laws with deep learning models holds promise, but challenges such as capturing dynamic electrophysiology accurately, accessing accurate domain knowledge, and quantifying prediction uncertainty persist. Integrating models into clinical workflows while ensuring interpretability and usability for healthcare professionals is essential. Overcoming these challenges will drive further research in CDTs.

Cardiac Activation Maps Reconstruction: A Comparative Study Between Data-Driven and Physics-Based Methods

One of the essential diagnostic tools of cardiac arrhythmia is activation mapping. Noninvasive current mapping procedures include electrocardiographic imaging. It allows reconstructing heart surface potentials from measured body surface potentials. Then, activation maps are generated using the heart surface potentials. Recently, a study suggests to deploy artificial neural networks to estimate activation maps directly from body surface potential measurements. Here we carry out a comparative study between the data-driven approach DirectMap and noninvasive classic technique based on reconstructed heart surface potentials using both Finite element method combined with L1-norm regularization (FEM-L1) and the spatial adaptation of Time-delay neural networks (SATDNN-AT). In this work, we assess the performance of the three approaches using a synthetic single paced-rhythm dataset generated on the atria surface. The results show that data-driven approach DirectMap quantitatively outperforms the two other methods. In fact, we observe an absolute activation time error and a correlation coefficient, respectively, equal to 7.20 ms, 93.2% using DirectMap, 14.60 ms, 76.2% using FEM-L1 and 13.58 ms, 79.6% using SATDNN-AT. In addition, results show that data-driven approaches (DirectMap and SATDNN-AT) are strongly robust against additive gaussian noise compared to FEM-L1.

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.

Localization of origins of premature ventricular contraction in the whole ventricle based on machine learning and automatic beat recognition from 12-lead ECG

OBJECTIVE

The localization of origins of premature ventricular contraction (PVC) is the key factor for the success of ablation of ventricular arrhythmias. Existing methods rely heavily on manual extraction of PVC beats, which limits their application to the automatic PVC recognition from long-term data recorded by ECG monitors before and during operation. In addition, research identifying PVC sources in the whole ventricle have not been reported. The purpose of this study was to validate the feasibility of localization of origins of PVC in the whole ventricle and to explore an automatic algorithm for recognition of PVC beats based on long-term 12-lead ECG.

APPROACH

This study included 249 patients with spontaneous PVCs or pacing-induced PVCs. A novel algorithm was used to automatically extract PVC beats from a massive amount of original ECG data, which was collected by different acquisition devices. After clustering and labelling, 374 sample groups, each containing dozens to hundreds of PVC beats, formed the entire dataset of 11 categories corresponding to 11 regions of PVC origins in the whole ventricle. To choose the best classification model for the current task, four machine learning methods, support vector machine (SVM), random forest (RF), gradient-boosting decision tree (GBDT) and Gaussian naïve Bayes (GNB), were compared by randomly selecting 70% of the entire dataset (sample groups = 257) for training and the remaining 30% (sample groups = 117) for testing. The average performance of each model was estimated by the bootstrap method using 1000 resampling trials.

MAIN RESULTS

For PVC beat recognition, the achieved testing accuracy, sensitivity and specificity is 97.6%, 98.3% and 96.7%, respectively. For localization purpose, the achieved testing accuracy varies slightly from 70.7% to 74.1% among four classifiers, and when neighboring regions were combined, the testing rank accuracy is improved to a range of 91.5% to 93.2%.

SIGNIFICANCE

The proposed algorithm can automatically recognize PVC beats and map them to one of the 11 regions in the whole ventricle. Owing to the high accuracy of PVC beat recognition and the capability to target the potential PVC origins in multi regions, it is expected to be a predominant technique being used in clinical settings to automatically analyze huge ECG data before and during operation so as to replace the tedious manual identification.

Non-invasive reconstruction of dynamic myocardial transmembrane potential with graph-based total variation constraints

Non-invasive reconstruction of electrophysiological activity in the heart is of great significance for clinical disease prevention and surgical treatment. The distribution of transmembrane potential (TMP) in three-dimensional myocardium can help us diagnose heart diseases such as myocardial ischemia and ectopic pacing. However, the problem of solving TMP is ill-posed, and appropriate constraints need to be added. The existing state-of-art method total variation minimisation only takes advantage of the local similarity in space, which has the problem of over-smoothing, and fails to take into account the relationship among frames in the dynamic TMP sequence. In this work, the authors introduce a novel regularisation method called graph-based total variation to make up for the above shortcomings. The graph structure takes the TMP value of a time sequence on each heart node as the criterion to establish the similarity relationship among the heart. Two sets of phantom experiments were set to verify the superiority of the proposed method over the traditional constraints: infarct scar reconstruction and activation wavefront reconstruction. In addition, experiments with ten real premature ventricular contractions patient data were used to demonstrate the accuracy of the authors’ method in clinical applications.

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.