A robust wavelet-based approach for dominant frequency analysis of atrial fibrillation in body surface signals.

OBJECTIVE: Atrial dominant frequency (DF) maps undergoing atrial fibrillation (AF) presented good spatial correlation with those obtained with the non-invasive body surface potential mapping (BSPM). In this study, a robust BSPM-DF calculation method based on wavelet analysis is proposed. APPROACH: Continuous wavelet transform along 40 scales in the pseudo-frequency range of 3-30 Hz is performed in each BSPM signal using a Gaussian mother wavelet. DFs are estimated from the intervals between the peaks, representing the activation times, in the maximum energy scale. The results are compared with the traditionally widely applied Welch periodogram and the robustness was tested on different protocols: increasing levels of white Gaussian noise, artificial DF harmonics presence and reduction in the number of leads. A total of 11 AF simulations and 12 AF patients are considered in the analysis. For each patient, intracardiac electrograms were acquired in 15 locations from both atria. The accuracy of both methods was assessed by calculating the absolute errors of the highest DF BSPM (HDF BSPM ) with respect to the atrial HDF, either simulated or intracardially measured, and assumed correct if ≤1 Hz. The spatial distribution of the errors between torso DFs and atrial HDFs were compared with atria driving mechanism locations. Torso HDF regions, defined as portions of the maps with [Formula: see text] Hz were identified and the percentage of the torso occuping these regions was compared between methods. The robustness of both methods to white Gaussian noise, ventricular influence and harmonics, and to lower spatial resolution BSPM lead layouts was analyzed: computer AF models (567 leads vs 256 leads down to 16 leads) and patient data (67 leads vs 32 and 16 leads). MAIN RESULTS: The proposed method allowed an improvement in non-invasive estimation of the atria HDF. For the models the median relative errors were 7.14% for the wavelet-based algorithm vs 60.00% for the Welch method; in patients, the errors were 10.03% vs 12.66%, respectively. The wavelet method outperformed the Welch approach in correct estimations of atrial HDFs in models (81.82% vs 45.45%, respectively) and patients (66.67% vs 41.67%). A low positive BSPM-DF map correlation was seen between the techniques (0.47 for models and 0.63 for patients), highlighting the overall differences in DF distributions. The wavelet-based algorithm was more robust to white Gaussian noise, residual ventricular activity and harmonics, and presented more consistent results in lead layouts with low spatial resolution. SIGNIFICANCE: Estimation of atrial HDFs using BSPM is improved by the proposed wavelet-based algorithm, helping to increase the non-invasive diagnostic ability in AF.

Adaptive step ODE algorithms for the 3D simulation of electric heart activity with graphics processing units.

In this paper we studied the implementation and performance of adaptive step methods for large systems of ordinary differential equations systems in graphics processing units, focusing on the simulation of three-dimensional electric cardiac activity. The Rush-Larsen method was applied in all the implemented solvers to improve efficiency. We compared the adaptive methods with the fixed step methods, and we found that the fixed step methods can be faster while the adaptive step methods are better in terms of accuracy and robustness.