Abstract-This contribution addresses the extraction of atrial activity (AA) from real electrocardiogram (ECG) recordings of atrial fibrillation (AF). We show the appropriateness of independent component analysis (ICA) to tackle this biomedical challenge when regarded as a blind source separation (BSS) problem. ICA is a statistical tool able to reconstruct the unobservable independent sources of bioelectric activity which generate, through instantaneous linear mixing, a measurable set of signals. The three key hypothesis that make ICA applicable in the present scenario are discussed and validated: 1) AA and ventricular activity (VA) are generated by sources of independent bioelectric activity; 2) AA and VA present non-Gaussian distributions; and 3) the generation of the surface ECG potentials from the cardioelectric sources can be regarded as a narrow-band linear propagation process. To empirically endorse these claims, an ICA algorithm is applied to recordings from seven patients with persistent AF. We demonstrate that the AA source can be identified using a kurtosis-based reordering of the separated signals followed by spectral analysis of the sub-Gaussian sources. In contrast to traditional methods, the proposed BSS-based approach is able to obtain a unified AA signal by exploiting the atrial information present in every ECG lead, which results in an increased robustness with respect to electrode selection and placement.

BACKGROUND Ablation is an effective therapy in patients with atrial fibrillation (AF) in which an electrical driver can be identified. OBJECTIVE The aim of this study was to present and discuss a novel and strictly noninvasive approach to map and identify atrial regions responsible for AF perpetuation. METHODS Surface potential recordings of 14 patients with AF were recorded using a 67-lead recording system. Singularity points (SPs) were identified in surface phase maps after band-pass filtering at the highest dominant frequency (HDF). Mathematical models of combined atria and torso were constructed and used to investigate the ability of surface phase maps to estimate rotor activity in the atrial wall. RESULTS The simulations show that surface SPs originate at atrial SPs, but not all atrial SPs are reflected at the surface. Stable SPs were found in AF signals during 8.3% ± 5.7% vs 73.1% ± 16.8% of the time in unfiltered vs HDF-filtered patient data, respectively (P < .01). The average duration of each rotational pattern was also lower in unfiltered than in HDF-filtered AF signals (160 ± 43 ms vs 342 ± 138 ms; P < .01), resulting in 2.8 ± 0.7 rotations per rotor. Band-pass filtering reduced the apparent meandering of surface HDF rotors by reducing the effect of the atrial electrical activity occurring at different frequencies. Torso surface SPs representing HDF rotors during AF were reflected at specific areas corresponding to the fastest atrial location. CONCLUSION Phase analysis of surface potential signals after HDF filtering during AF shows reentrant drivers localized to either the left atrium or the right atrium, helping in localizing ablation targets.

The analysis and characterization of atrial tachyarrhythmias requires, in a previous step, the extraction of the atrial activity (AA) free from ventricular activity and other artefacts. This contribution adopts the blind source separation (BSS) approach to AA estimation from multilead electrocardiograms (ECGs). Previously proposed BSS methods for AA extraction--e.g., independent component analysis (ICA)--exploit only the spatial diversity introduced by the multiple spatially-separated electrodes. However, AA typically shows certain degree of temporal correlation, with a narrowband spectrum featuring a main frequency peak around 3.5-9 Hz. Taking advantage of this observation, we put forward a novel two-step BSS-based technique which exploits both spatial and temporal information contained in the recorded ECG signals. The spatiotemporal BSS algorithm is validated on simulated and real ECGs from a significant number of atrial fibrillation (AF) and atrial flutter (AFL) episodes, and proves consistently superior to a spatial-only ICA method. In simulated ECGs, a new methodology for the synthetic generation of realistic AF episodes is proposed, which includes a judicious comparison between the known AA content and the estimated AA sources. Using this methodology, the ICA technique obtains correlation indexes of 0.751, whereas the proposed approach obtains a correlation of 0.830 and an error in the estimated signal reduced by a factor of 40%. In real ECG recordings, we propose to measure performance by the spectral concentration (SC) around the main frequency peak. The spatiotemporal algorithm outperforms the ICA method, obtaining a SC of 58.8% and 44.7%, respectively.

Atrial fibrillation (AF) is the most common arrhythmia encountered in clinical practice. Neither the natural history of AF nor its response to therapy is sufficiently predictable by clinical and echocardiographic parameters. The purpose of this article is to describe technical aspects of novel electrocardiogram (ECG) analysis techniques and to present research and clinical applications of these methods for characterization of both the fibrillatory process and the ventricular response during AF. Atrial fibrillatory frequency (or rate) can reliably be assessed from the surface ECG using digital signal processing (extraction of atrial signals and spectral analysis). This measurement shows large inter-individual variability and correlates well with intra-atrial cycle length, a parameter which appears to have primary importance in AF maintenance and response to therapy. AF with a low fibrillatory rate is more likely to terminate spontaneously and responds better to antiarrhythmic drugs or cardioversion, whereas high-rate AF is more often persistent and refractory to therapy. Ventricular responses during AF can be characterized by a variety of methods, which include analysis of heart rate variability, RR-interval histograms, Lorenz plots, and non-linear dynamics. These methods have all shown a certain degree of usefulness, either in scientific explorations of atrioventricular (AV) nodal function or in selected clinical questions such as predicting response to drugs, cardioversion, or AV nodal modification. The role of the autonomic nervous system for AF sustenance and termination, as well as for ventricular rate responses, can be explored by different ECG analysis methods. In conclusion, non-invasive characterization of atrial fibrillatory activity and ventricular response can be performed from the surface ECG in AF patients. Different signal processing techniques have been suggested for identification of underlying AF pathomechanisms and prediction of therapy efficacy.

Background Ablation of high-frequency sources in patients with atrial fibrillation (AF) is an effective therapy to restore sinus rhythm. However, this strategy may be ineffective in patients without a significant dominant frequency (DF) gradient. The aim of this study was to investigate whether sites with high-frequency activity in human AF can be identified noninvasively, which should help intervention planning and therapy. Methods and Results In 14 patients with a history of AF, 67-lead body surface recordings were simultaneously registered with 15 endocardial electrograms from both atria including the highest DF site, which was predetermined by atrial-wide real-time frequency electroanatomical mapping. Power spectra of surface leads and the body surface location of the highest DF site were compared with intracardiac information. Highest DFs found on specific sites of the torso showed a significant correlation with DFs found in the nearest atrium (ρ=0.96 for right atrium and ρ=0.92 for left atrium) and the DF gradient between them (ρ=0.93). The spatial distribution of power on the surface showed an inverse relationship between the frequencies versus the power spread area, consistent with localized fast sources as the AF mechanism with fibrillatory conduction elsewhere. Conclusions Spectral analysis of body surface recordings during AF allows a noninvasive characterization of the global distribution of the atrial DFs and the identification of the atrium with the highest frequency, opening the possibility for improved noninvasive personalized diagnosis and treatment.

Introduction Ablation of high dominant frequency (DF) sources in patients with atrial fibrillation (AF) is an effective treatment option for paroxysmal AF. The aim of this study was to evaluate the accuracy of noninvasive estimation of DF and electrical patterns determination by solving the inverse problem of the electrocardiography. Methods Four representative AF patients with left-to-right and right-to-left atrial DF patterns were included in the study. For each patient, intracardiac electrograms from both atria were recorded simultaneously together with 67-lead body surface recordings. In addition to clinical recordings, realistic mathematical models of atria and torso anatomy with different DF patterns of AF were used. For both mathematical models and clinical recordings, inverse-computed electrograms were compared to intracardiac electrograms in terms of voltage, phase, and frequency spectrum relative errors. Results Comparison between intracardiac and inverse computed electrograms for AF patients showed 8.8 ± 4.4% errors for DF, 32 ± 4% for voltage, and 65 ± 4% for phase determination. These results were corroborated by mathematical simulations showing that the inverse problem solution was able to reconstruct the frequency spectrum and the DF maps with relative errors of 5.5 ± 4.1%, whereas the reconstruction of the electrograms or the instantaneous phase presented larger relative errors (i.e., 38 ± 15% and 48 ± 14 % respectively, P < 0.01). Conclusions Noninvasive reconstruction of atrial frequency maps can be achieved by solving the inverse problem of electrocardiography with a higher accuracy than temporal distribution patterns.

Local stretching increases the electrophysiological heterogeneity of myocardium and accelerates and increases the complexity of VF in the stretched area, without significantly modifying the occurrences of the types of VF activation patterns in the nonstretched zone.

Analysis of atrial rhythm is important in the treatment and management of patients with atrial fibrillation. Several algorithms exist for extracting the atrial signal from the electrocardiogram (ECG) in atrial fibrillation, but there are few reports on how well these techniques are able to recover the atrial signal. We assessed and compared three algorithms for extracting the atrial signal from the 12-lead ECG. The 12-lead ECGs of 30 patients in atrial fibrillation were analyzed. Atrial activity was extracted by three algorithms, Spatiotemporal QRST cancellation (STC), principal component analysis (PCA), and independent component analysis (ICA). The amplitude and frequency characteristics of the extracted atrial signals were compared between algorithms and against reference data. Mean (standard deviation) amplitude of QRST segments of V1 was 0.99 (0.54) mV, compared to 0.18 (0.11) mV (STC), 0.19 (0.13) mV (PCA), and 0.29 (0.22) mV (ICA). Hence, for all algorithms there were significant reductions in the amplitude of the ventricular activity compared with that in V1. Reference atrial signal amplitude in V1 was 0.18 (0.11) mV, compared to 0.17 (0.10) mV (STC), 0.12 (0.09) mV (PCA), and 0.18 (0.13) mV (ICA) in the extracted atrial signals. PCA tended to attenuate the atrial signal in these segments. There were no significant differences for any of the algorithms when comparing the amplitude of the reference atrial signal with that of the extracted atrial signals in segments in which ventricular activity had been removed. There were no significant differences between algorithms in the frequency characteristics of the extracted atrial signals. There were discrepancies in amplitude and frequency characteristics of the atrial signal in only a few cases resulting from notable residual ventricular activity for PCA and ICA algorithms. In conclusion, the extracted atrial signals from these algorithms exhibit very similar amplitude and frequency characteristics. Users of these algorithms should be observant of residual ventricular activities which can affect the analysis of the fibrillatory waveform in clinical practice.