- Research
- Open Access

# Identification of recurring patterns in fractionated atrial electrograms using new transform coefficients

- Edward J Ciaccio
^{1}Email author, - Angelo B Biviano
^{1}, - William Whang
^{1}and - Hasan Garan
^{1}

**11**:4

https://doi.org/10.1186/1475-925X-11-4

© Ciaccio et al.; licensee BioMed Central Ltd 2012

**Received:**20 November 2011**Accepted:**11 January 2012**Published:**19 January 2012

## Abstract

### Background

Identification of recurrent patterns in complex fractionated atrial electrograms (CFAE) has been used to differentiate paroxysmal from persistent atrial fibrillation (AF). Detection of the atrial CFAE patterns might therefore be assistive in guiding radiofrequency catheter ablation to drivers of the arrhythmia. In this study a technique for robust detection and classification of recurrent CFAE patterns is described.

### Method

CFAE were obtained from the four pulmonary vein ostia, and from the anterior and posterior left atrium, in 10 patients with paroxysmal AF and 10 patients with longstanding persistent AF (216 recordings in total). Sequences 8.4 s in length were analyzed (8,192 sample points, 977 Hz sampling). Among the 216 sequences, two recurrent patterns A and B were substituted for 4 and 5 of the sequences, respectively. To this data, random interference, and random interference + noise were separately added. Basis vectors were constructed using a new transform that is derived from ensemble averaging. Patterns A and B were then detected and classified using a threshold level of Euclidean distance between spectral signatures as constructed with transform coefficients.

### Results

In the presence of interference, sensitivity to detect and distinguish two patterns A and B was 96.2%, while specificity to exclude nonpatterns was 98.0%. In the presence of interference + noise, sensitivity was 89.1% while specificity was 97.0%.

### Conclusions

Transform coefficients computed from ensemble averages can be used to succinctly quantify synchronized patterns present in AF data. The technique is useful to automatically detect recurrent patterns in CFAE that are embedded in interference without user bias. This quantitation can be implemented in real-time to map the AF substrate prior to and during catheter ablation.

## Keywords

- atrial fibrillation
- catheter ablation
- ensemble averaging
- pattern recognition
- transform

## Background

Radiofrequency catheter ablation is often used for successful treatment of atrial fibrillation (AF), and is guided in part by the morphology of electrograms recorded from the catheter tip. Of particular interest are complex fractionated atrial electrograms (CFAE), which are composed of multiple deflections with varying baseline, or continuous deflections with low voltage [1]. The CFAE may represent the arrhythmogenic substrate for AF. Ablating CFAE can increase the cycle length of the arrhythmia, suggesting the importance of some of these regions as drivers to maintain AF [2]. Ablation of CFAE may improve outcome of the catheter ablation procedure [3]. Progress in the development of signal processing algorithms to identify CFAE can improve the efficacy of ablation strategies [4, 5]. Although AF may originate in the pulmonary veins (PV), in many patients, particularly those with the longstanding persistent type, other areas of the left atrium must be ablated to successfully stop the arrhythmia. However, the precise characteristics of atrial electrograms that would suggest that a particular area should be ablated is currently the subject of debate. Moreover, differing techniques for CFAE quantitation do not necessarily identify the same areas for ablation [6, 7]. Furthermore, CFAE identified outside the PVs often represent a large surface area of tissue that when ablated in its entirety, can increase procedure time and the possibility of patient morbidity. Thus the need to recognize CFAE with special quantitative characteristics that when ablated can improve outcome. The presence of morphologic differences and repetitive patterns in CFAE have been observed and quantified as a way to distinguish paroxysmal from longstanding persistent AF patients [8, 9]. If such patterns could be distinguished from one recorded sequence to the next, then the substrate could be mapped based on pattern recurrence, which is likely related to degree of arrhythmogenicity [8, 9].

In previous work we described a method of spectral estimation and transformation for analysis of atrial fibrillation data [10–12]. Since this transform is data-driven, the orthogonal basis vectors are unique to the particular data set being analyzed. Thus it is possible to extract the original patterns from which the signals were generated from any additive noise and interference that may be present. As is shown in this study, if a recurrent pattern is present in CFAE, it can be detected by generating and then comparing transform coefficients, and is robust to presence of additive random noise and interference. Pattern recognition techniques are then used to distinguish two recurring patterns from nonpatterns present in the CFAE [13, 14].

## Methods

### A Clinical data acquisition

Electrograms were recorded in a series of twenty patients referred to the Columbia University Medical Center cardiac electrophysiology (EP) laboratory for catheter ablation of AF. These recordings were obtained prospectively as approved by the Internal Review Board at Columbia University Medical Center, but analyzed retrospectively after the catheter ablation procedures were completed using our standard clinical protocols. Ten patients had documented clinical paroxysmal (acute) AF, with a normal sinus rhythm as their baseline rhythm in the electrophysiology laboratory. Atrial fibrillation was induced by burst pacing from the coronary sinus or the lateral right atrial wall, and the arrhythmia persisted for at least 10 minutes for those signals to be included in the retrospective analysis. Ten other patients had persistent (longstanding) AF, and had been in AF without interruption for 1-6 years prior to the catheter mapping and ablation procedure. Bipolar electrograms of at least 10 seconds in duration, recorded from the distal ablation electrode during arrhythmia, were bandpass filtered by the system at acquisition to remove baseline drift and high frequency noise (30-500 Hz), sampled at 977 Hz, and stored. Although the bandpass high corner was slightly greater than the Nyquist frequency, negligible signal energy resides in the region [11].

Only digitized signals identified as CFAE by two cardiac electrophysiologists were included in the retrospective analysis. The CFAE recordings were obtained from two sites outside the ostia of each of the four PVs. Similar recordings were obtained at two sites on the endocardial surface of the left atrial free wall, one in the mid-posterior wall, and another on the anterior ridge at the base of the left atrial appendage. From each of these recordings, 8.4-second sequences (8,192 sample points) were extracted and analyzed. A total of 240 such sequences were acquired - 120 from paroxysmal and 120 from longstanding AF patients. Subsequently, only 216 of the recordings were confirmed as CFAE, and only these were used for subsequent analysis. All CFAE signals were normalized to mean zero and unity variance prior to further processing [12].

### B Transform coefficients and the spectral signature

*w*is then calculated by averaging

*n*successive mean zero segments of signal $\underline{x}$ having a length

*N*, with each segment being of integer length

*w*:

*w*that are extracted from

*x*and summed. Thus ${\underline{e}}_{w}$ is computed by summing segments of the signal having period

*w*sample points. The summing matrix ${\mathbf{U}}_{w}$ is padded at right and bottom edge if $N/w$ is not an integer, as described elsewhere [12]. The relationship between frequency

*f*and period

*w*is given by:

*w*is given by:

*w*summations used to calculate the power, and the

*n*summations used to form each ensemble average. The transformation matrix is:

*x*of length

*N*can then be decomposed using the linear transformation:

*n*is the number of summations for ensemble averaging, and ${\underline{a}}_{w}$ and $\underline{x}$ are $N\times 1$ in dimension. The basis vectors ${\underline{a}}_{w}$ for all

*w*characterize the periodic behaviors present in signal vector $\underline{x}$, and they are linearly independent except when small integer relationships exist between periods ${w}_{i}$ and ${w}_{j}$ of any particular pair [12]. Columnwise, each identity submatrix in Equation 5b extracts and sums one segment of

*w*sample points in $\underline{x}$ (Equation 6), with the sum total being projected onto the canonical basis. Rowwise the identity matrices serve to repeat the ensemble average of length

*w*over a total length

*N*during construction of ${\underline{a}}_{w}$. Thus the transformation matrix of Equations 5a and 5b decomposes the signals into periodic ensemble averages. These orthogonal basis vectors can be used to project signal $\underline{x}$ into ensemble space:

*w*is correlated with the ensemble average at

*w*(LHS), the resulting transform coefficient is the ensemble average power at

*w*(RHS). The power spectrum can either be plotted versus period

*w*, or versus frequency

*f*, and to level the noise floor, which depends on the number of summations for averaging, it is scaled by $\sqrt{n}$ when graphed [12]. The transformation process can be analyzed by expanding Equation 7:

In the middle part of Equation 8, when the mathematical operations are done from right to left, then starting with signal $\underline{x}$, the ensemble average is generated by ${\mathbf{U}}_{w}\phantom{\rule{0.25em}{0ex}}\underline{x}$, orthogonal basis vectors are formed by ${\mathbf{U}}_{w}^{T}\phantom{\rule{0.25em}{0ex}}{\mathbf{U}}_{w}\phantom{\rule{0.25em}{0ex}}\underline{x}$, and transform coefficients by ${\underline{x}}^{T}\phantom{\rule{0.25em}{0ex}}{\mathbf{U}}_{w}^{T}\phantom{\rule{0.25em}{0ex}}{\mathbf{U}}_{w}\phantom{\rule{0.25em}{0ex}}\underline{x}$.

*w*are not power spectra. Rather, they represent the power in each basis vector derived from $\underline{z}$ that is correlated to $\underline{x}$ and to $\underline{y}$, respectively. Hence they can be defined as spectral signatures, having similarities to the power spectra of $\underline{x}$ and $\underline{y}$ respectively, but using, for reconstruction, basis vectors that are only partially correlated with each signal. Considering the relationship between $\underline{x}$, $\underline{y}$, and $\underline{z}$, several inequalities should be noted:

Equation 15 states that the sum of the spectral signatures of $\underline{x}$ and $\underline{y}$ with respect to $\underline{z}$ equals the power spectrum of $\underline{z}$.

Now suppose that $\underline{z}$ is an average of many signals, some of which contain a particular pattern A or a different pattern B. These two repeating patterns A and B will be reinforced, and the random content will be reduced, by the summation that forms $\underline{z}$. The basis vectors formed from this average will therefore mostly be constructed from a combination of the features from pattern A and those of pattern B. Hence we would expect that if a particular input signal happened to contain pattern A or pattern B, then the resulting transform coefficients ${C}_{w}$ would be similar to ${P}_{w}$. In contrast, if the input signal is a random vector unrelated to pattern A or B, the resulting transform coefficients ${C}_{w}$ for all *w* would approach zero and contain both positive and negative values, due to the lack of correlation of the signal with the basis vectors. These properties can be exploited for detection of two recurring patterns A and B, as described in the next section.

### C Paradigm for recognition of recurrent patterns in atrial fibrillation signals

*m*random signals, a few of which contain pattern A or pattern B embedded in additive random noise and interference, as represented at top left in Figure 1. If the

*m*signals are averaged, only recurring patterns will be reinforced in the resulting mean signal (top middle, Figure 1). When the power spectrum

*P*is generated from this mean signal (middle portion of Figure 1), it will contain elements of any recurring patterns. If the transform coefficients of each of the

*m*individual signals are obtained using the basis vectors constructed from the mean signal, the resulting spectral signatures ${S}_{i}$ will be similar to

*P*only for those signals which contain a recurring pattern. As suggested previously [12], the Euclidean distance (ED) between power spectrum

*P*and a particular spectral signature ${S}_{i}$ will be large if ${S}_{i}$ lacks a recurring pattern, and small if it contains a recurring pattern. A first threshold level (Th1) of Euclidean distance ED1 can therefore be used to detect the presence of a candidate pattern in ${S}_{i}$, for all

*i*(lower right, Figure 1). Once the candidate patterns in the series are detected, their spectral signatures can be compared one to another using a second Euclidean distance and threshold (ED2 and Th2). The Euclidean distance will be shorter if ${S}_{i}$ and ${S}_{j}$ are generated from two signals

*i*and

*j*containing the same pattern. The Euclidean distance will be longer if ${S}_{i}$ and ${S}_{j}$ are generated from signals containing different patterns, or if one or both of the signals contain no pattern. The latter can occur if the first threshold is set to a longer value of Euclidean distance, so that some random signals lacking patterns are initially identified as candidate patterns. Therefore, threshold selection is a tradeoff between excluding actual patterns (shorter Euclidean distance) versus including nonpatterns (longer Euclidean distance).

The method described above can be automated to detect recurring patterns without manual intervention, except to set threshold level Th1 for pattern detection, and threshold level Th2 to cluster and classify the detected patterns. For the procedure to work, the patterns contained in the signals must be synchronous so that they reinforce upon averaging. This can often be achieved either by simultaneously recording from many atrial sites, or by successive recording using a suitable trigger such as an *F* wave peak when present in the electrocardiogram.

Of the final series of $m-2=214$ signals with interferences, 27 thus contained one of two recurring patterns A and B due to the method of combination described in Equation 16. This final set was used for transformation, with patterns A and B detected according to the flow diagram in Figure 1. Since the values of *w* ranged from 50 to 1,000, each spectrum and spectral signature was 951-dimensional. The process was repeated with different patterns selected at random for a total of 10 trials. As an additional test, all of the steps above were repeated with random noise added to each CFAE as well as interference. The random noise vectors with Gaussian distribution were approximately mean zero, and truncated with a standard deviation of ±2 millivolts, twice the normalized CFAE standard deviation.

## Results

**Euclidean distance between candidate patterns**

Pattern | A | n | A | B | A | A | A | B | B | n | B |
---|---|---|---|---|---|---|---|---|---|---|---|

A | 0.000 | 0.142 | 0.056 | 0.133 | 0.092 | 0.074 | 0.065 | 0.128 | 0.135 | 0.221 | 0.135 |

n | 0.142 | 0.000 | 0.151 | 0.149 | 0.166 | 0.123 | 0.131 | 0.188 | 0.156 | 0.143 | 0.122 |

A | 0.056 | 0.151 | 0.000 | 0.136 | 0.097 | 0.068 | 0.068 | 0.143 | 0.145 | 0.214 | 0.138 |

B | 0.133 | 0.149 | 0.136 | 0.000 | 0.167 | 0.127 | 0.144 | 0.101 | 0.096 | 0.185 | 0.063 |

A | 0.092 | 0.166 | 0.097 | 0.167 | 0.000 | 0.095 | 0.101 | 0.196 | 0.206 | 0.271 | 0.151 |

A | 0.074 | 0.123 | 0.068 | 0.127 | 0.095 | 0.000 | 0.082 | 0.156 | 0.124 | 0.191 | 0.116 |

A | 0.065 | 0.131 | 0.068 | 0.144 | 0.101 | 0.082 | 0.000 | 0.156 | 0.151 | 0.212 | 0.152 |

B | 0.128 | 0.188 | 0.143 | 0.101 | 0.196 | 0.156 | 0.156 | 0.000 | 0.105 | 0.241 | 0.102 |

B | 0.135 | 0.156 | 0.145 | 0.096 | 0.206 | 0.124 | 0.151 | 0.105 | 0.000 | 0.161 | 0.104 |

n | 0.221 | 0.143 | 0.214 | 0.185 | 0.271 | 0.191 | 0.212 | 0.241 | 0.161 | 0.000 | 0.168 |

B | 0.135 | 0.122 | 0.138 | 0.063 | 0.151 | 0.116 | 0.152 | 0.102 | 0.104 | 0.168 | 0.000 |

**Statistics for pattern classification**

Trial number | sen: int | spe: int | sen: int+n | sen: int+n |
---|---|---|---|---|

1 | 100.0 | 100.0 | 91.1 | 100.0 |

2 | 97.8 | 100.0 | 97.8 | 95.0 |

3 | 93.3 | 100.0 | 82.2 | 100.0 |

4 | 95.6 | 100.0 | 88.9 | 100.0 |

5 | 93.3 | 100.0 | 88.9 | 100.0 |

6 | 91.1 | 80.0 | 88.9 | 75.0 |

7 | 95.6 | 100.0 | 88.9 | 100.0 |

8 | 97.8 | 100.0 | 91.1 | 100.0 |

9 | 97.8 | 100.0 | 84.4 | 100.0 |

10 | 100.0 | 100.0 | 88.9 | 100.0 |

mean | 96.2 ± 3.0 | 98.0 ± 6.3 | 89.1 ± 4.1 | 97.0 ± 7.9 |

## Discussion

### Summary

In this study a new transform was used to characterize recurring patterns in CFAE. First, ensemble averages were computed from signal segments of length *w*, repeated for all *w* in the frequency range of interest as given by Equation 3. From each ensemble average an orthogonal basis vector is constructed by repeating the ensemble average of length *w* for the entire signal length *N* (Equation 6). The inner product between basis vector and original signal produces a transform coefficient, which is the signal power at that frequency (Equation 7). The power spectrum is a plot of the entire series of transform coefficients versus frequency. Transform coefficients resulting from the inner product of one signal with the basis vectors of another signal can take on negative as well as positive values, and will have an average level near zero if the signals are uncorrelated (Equations 11a and 11b). The correlation coefficients formed from correlated signal $\underline{x}$ with the basis vectors of $\underline{z}$, as described in the Methods, can be similar to the spectrum of $\underline{x}$ and is termed the spectral signature. Transform coefficients were used to detect two recurring patterns in a sequence of CFAE, embedded in interference and random noise, and to distinguish them from each other and from nonpatterns. The method was implemented and repeated for 10 trials. No manual intervention was used except to set initial threshold levels of Euclidean distance for identification of correlated content, i.e., for pattern extraction, and to distinguish the extracted patterns (Figure 1).

### Prior work in pattern recognition

Our study made use of correlation in the frequency domain to discern repetitive patterns from nonpatterns in CFAE. Electrogram pattern recognition using correlation waveform analysis has also been found useful to discern electrograms arising from arrhythmia in another study [15]. In this prior work, the electrogram shape at any particular recording site was shown to remain stable during electrophysiologic study, even during such interventions as overdrive pacing and infusion of pharmacologic agents, which included epinephrine and isoproterenol [16]. Thus we expect that our method for detecting and discerning patterns in CFAE will likely be robust to the typical interventions that are done during the course of a clinical electrophysiologic study. We also found the transformation to the frequency domain to be useful to characterize each signal and its relationship to the mean signal. Similarly, in a prior study the wavelet transform was found useful to discern atrial electrogram patterns by categorizing them into one of four classes of fractionation based on frequency and phasic relationships [4, 5]. Although the patterns used in our study were synchronized artificially for recognition, multielectrode recordings will likely be useful to simultaneously acquire data from many recording sites in which any patterns present will be synchronized. This would simplify the mapping procedure using the paradigm of Figure 1. Elsewhere, it has been shown that multiple simultaneously obtained recordings are indeed useful for rapid and accurate classification of CFAE patterns [16, 17]. Although we have not yet proposed a paradigm to relate the detected patterns to the electrophysiologic substrate, it is planned to develop such a technique for a future prospective study. Adopting a standardized description of CFAE morphology and use of reproducible methodology would enable ease of comparison between clinical trials [17].

Regarding the possible relationship of observed pattern to electrophysiologic properties of the substrate, the complexity of CFAE, as determined by pattern type, is believed to be related to the degree of organization of electrical activity [18, 19]. To reduce complexity, adaptive template matching can be used to normalize the signals with respect to atrial cycle length (*x*-scale), as well as amplitude differences (*y*-scale) [20]. This would be useful to compare patterns in persistent versus paroxysmal AF independent of the dominant frequency, which tends to be higher in the persistent type, and to compare regions of atrial fractionation having low electrogram amplitude to areas with higher electrogram amplitude [11].

### Potential advantages of the new method

The spectral signature is a graph of the correlated content between two signals in frequency space, which can be exploited for pattern recognition (see the Methods). If a series of signals is averaged and basis vectors of the mean are used to obtain the spectral signature of each individual signal, then there will be correlation between the spectrum of the mean, and the spectral signature of the individual signal, when the individual signal contains a synchronous pattern that recurs within the series. By measuring the Euclidean distance between all individual signals having spectral signatures similar to the power spectrum of the mean signal, patterns contained in the sequence can be identified, distinguished from one another, and distinguished from nonpatterns when the nonpatterns are mostly uncorrelated with respect to the mean signal. Thus the new technique is potentially useful to automatically identify and distinguish repetitive patterns present in a series of signals, once threshold levels for the Euclidean distance estimate to detect candidate patterns, and to discern patterns, are established. In the latter step, if multiple patterns are present they can also be discerned using a single threshold level, since the Euclidean distance will be short only with respect to members of the same class.

We found that best threshold levels are sensitive to the degree of additive random noise and interference, and differed when only interference was added versus interference + random noise (0.105 versus 0.132, respectively). Once the threshold levels were established based on a test trial, they were much less sensitive to particular patterns and patients i.e., the results for each trial were similar (Table 2).

### Clinical correlates

When considering the potential for recurrent patterns to appear in CFAE signals, fractionated electrogram deflection morphology is probably mechanism-dependent. Possible sources of fractionation include areas of slow electrical activation, wavefront collision, anchor points at regions driven by a reentrant circuit, and presence of multiple activation wavelets as triggered, for example, from ganglionated plexi [1]. These mechanisms are potentially distinguished both by spatiotemporal occurrence, and as described in this study by the frequency characteristics of the actual patterns. When using an automated technique for pattern detection, user bias is eliminated as a variable for defining ablation targets [1, 13, 14]. A goal of CFAE pattern recognition software is to visualize the spatial distribution of CFAE for catheter ablation [21]. Our technique can be used for spatial mapping of CFAE by type of pattern present at each recording site. CFAE are observed at both PVs and the left atrial free wall, although their occurrence at the free wall is more common in persistent AF [21]. Acute and longstanding persistent AF were shown by our group to have CFAE that differ in morphologic characteristics and degree of repetitiveness [8, 9]. Therefore, patterns in CFAE and their frequency of occurrence are likely to differ by AF type, and to present differently during mapping procedures. Although in previous work, recurring features in the form of individual deflections were measured on the order of tens of milliseconds in duration, in the current study, patterns in the data were measured over the entire sequence length used, 8,192 sample points (8.4 s). Thus the new technique assumes statistical stationarity over this time interval. Indeed, it has been shown in a previous study that 8-second lengths are ideal for characterizing the dominant frequency of CFAE [22]. Use of a noncontact catheter during clinical electrophysiologic study would enable simultaneous recording of CFAE from multiple sites. From the simultaneously recorded CFAE, far-field electrical activation patterns originating from distant drivers can be detected and localized [23]. Analysis of atrial electrogram patterns using a noncontact catheter is planned for future work. Recurring patterns in data obtained from sequential recording sites using a standard catheter may also be detectable when synchronized by a trigger such as the peak of the *F* wave in the electrocardiogram, also planned as a future study.

### Conclusions

A method was developed to recognize recurring patterns in AF data using a data-driven transform. It was shown to be efficacious for detection of recurring patterns with additive interferences and noise, simulated from actual data obtained from paroxysmal and persistent AF patients. Although the additive interference and noise rendered the original patterns almost unrecognizable by visual observation (Figures 2 6 7 8 and 9), the spectral signatures of each signal with additive interference and random noise, as related to the correlated content in the mean signal, were useful for detecting patterns, for distinguishing between two different patterns that were present, and for discerning patterns from nonpatterns using two threshold Euclidean distances. This paradigm may also be useful to develop a mechanistic understanding of paroxysmal and persistent AF, because presence of recurrent patterns can be compared and contrasted between the two AF types. It may also be useful to apply to other types of biomedical data such as ventricular tachyarrhythmias, and to videocapsule images of the small intestine, where spectral estimation from signal averaging has been described previously [24, 25].

### Limitations

Although the method described in this study was accurate for detecting simulated recurring patterns embedded in interference, real patterns may differ from one instance to the next, for example as caused by spatial and temporal jitter, which can reduce sensitivity and specificity [11]. Furthermore, the technique was shown to be useful for synchronous pattern data only. If the patterns are highly out of phase, they will not be reinforced when the mean signal is formed (Figure 1), and thus will not appear as correlated content. Therefore, synchronized pattern data is necessary for accuracy, as would more likely be effected by acquiring individual signal data from many recording sites simultaneously, or by synchronizing successively recorded data, for example by triggering the onset of the atrial recording to the presence of an *F* wave in the electrocardiogram.

## Author contributions

EJC developed the mathematical transform, conducted the data analysis, and wrote the manuscript. ABB, WW, and HG made helpful suggestions, provided the clinical data, and determined which recordings were complex fractionated atrial electrograms. All authors have read and approved the final manuscript.

## Declarations

### Acknowledgements

No outside funding was received for this study.

## Authors’ Affiliations

## References

- Monir G, Pollak SJ:
**Consistency of the CFAE phenomena using custom software for automated detection of complex fractionated atrial electrograms (CFAEs) in the left atrium during atrial fibrillation.***J Cardiovasc Electrophysiol*2008,**19:**915–919. 10.1111/j.1540-8167.2008.01153.xView ArticleGoogle Scholar - Hunter RJ, Diab I, Tayebjee M, Richmond L, Sporton S, Earley MJ, Schilling RJ:
**Characterization of fractionated atrial electrograms critical for maintenance of atrial fibrillation: a randomized, controlled trial of ablation strategies (the CFAE AF trial).***Circ Arrhythm Electrophysiol*2011,**4:**622–629. 10.1161/CIRCEP.111.962928View ArticleGoogle Scholar - Hunter RJ, Berriman TJ, Diab I, Baker V, Finlay M, Richmond L, Duncan E, Kamdar R, Thomas G, Abrams D, Dhinoja M, Sporton S, Earley MJ, Schilling RJ:
**Long-term efficacy of catheter ablation for atrial fibrillation: impact of additional targeting of fractionated electrograms.***Heart*2010,**96:**1372–1378. 10.1136/hrt.2009.188128View ArticleGoogle Scholar - Kremen V, Lhotská L, Macas M, Cihák R, Vancura V, Kautzner J, Wichterle D:
**A new approach to automated assessment of fractionation of endocardial electrograms during atrial fibrillation.***Physiol Meas*2008,**29:**1371–1381. 10.1088/0967-3334/29/12/002View ArticleGoogle Scholar - Kremen V, Kordík P, Lhotská L:
**Comparison of several classifiers to evaluate endocardial electrograms fractionation in human.***Conf Proc IEEE Eng Med Biol Soc*2009,**2009:**2502–2505.Google Scholar - Lee G, Roberts-Thomson K, Madry A, Spence S, Teh A, Heck PM, Kumar S, Kistler PM, Morton JB, Sanders P, Kalman JM:
**Relationship among complex signals, short cycle length activity, and dominant frequency in patients with long-lasting persistent AF: a high-density epicardial mapping study in humans.***Heart Rhythm*2011,**8:**1714–1719. 10.1016/j.hrthm.2011.05.021View ArticleGoogle Scholar - Orlov MV:
**A farewell to arms: Are complex fractionated atrial electrograms doomed as a target for AF ablation?***Heart Rhythm*2011,**8:**1720–1721. 10.1016/j.hrthm.2011.06.013View ArticleGoogle Scholar - Ciaccio EJ, Biviano AB, Whang W, Gambhir A, Garan H:
**Different characteristics of complex fractionated atrial electrograms in acute paroxysmal versus long-standing persistent atrial fibrillation.***Heart Rhythm*2010,**7:**1207–1215. 10.1016/j.hrthm.2010.06.018View ArticleGoogle Scholar - Ciaccio EJ, Biviano AB, Whang W, Gambhir A, Einstein AJ, Garan H:
**Differences in repeating patterns of complex fractionated left atrial electrograms in longstanding persistent as compared with paroxysmal atrial fibrillation.***Circ Arrhythm Electrophysiol*2011,**4:**470–477. 10.1161/CIRCEP.110.960153View ArticleGoogle Scholar - Ciaccio EJ, Biviano AB, Whang W, Wit AL, Garan H, Coromilas J:
**New methods for estimating local electrical activation rate during atrial fibrillation.***Heart Rhythm*2009,**6:**21–32. 10.1016/j.hrthm.2008.10.016View ArticleGoogle Scholar - Ciaccio EJ, Biviano AB, Whang W, Wit AL, Coromilas J, Garan H:
**Optimized measurement of activation rate at left atrial sites with complex fractionated electrograms during atrial fibrillation.***J Cardiovasc Electrophysiol*2010,**21:**133–143. 10.1111/j.1540-8167.2009.01595.xView ArticleGoogle Scholar - Ciaccio EJ, Biviano AB, Whang W, Coromilas J, Garan H:
**A new transform for the analysis of complex fractionated atrial electrograms.***Biomed Eng Online*2011,**10:**35. 10.1186/1475-925X-10-35View ArticleGoogle Scholar - Ciaccio EJ, Dunn SM, Akay M:
**Biosignal pattern-recognition and interpretation systems. 1. Fundamental concepts.***IEEE Eng Med Biol Mag*1993,**12:**89–98.View ArticleGoogle Scholar - Ciaccio EJ, Dunn SM, Akay M:
**Biosignal pattern-recognition and interpretation systems. 3. Clustering and classification.***IEEE Eng Med Biol Mag*1994,**13:**129–135.View ArticleGoogle Scholar - Finelli CJ, DiCarlo LA, Jenkins JM, Winston SA, Li PC:
**Effects of increased heart rate and sympathetic tone on intraventricular electrogram morphology.***Am J Cardiol*1991,**68:**1321–1328. 10.1016/0002-9149(91)90239-HView ArticleGoogle Scholar - Andrikopoulos G, Tzeis S, Vardas PE:
**Invasive therapy for atrial fibrillation: recent developments in ablation, navigation and mapping technology.***Heart*2011,**97:**237–243. 10.1136/hrt.2009.190017View ArticleGoogle Scholar - Hunter RJ, Diab I, Thomas G, Duncan E, Abrams D, Dhinoja M, Sporton S, Earley MJ, Schilling RJ:
**Validation of a classification system to grade fractionation in atrial fibrillation and correlation with automated detection systems.***Europace*2009,**11:**1587–1596. 10.1093/europace/eup351View ArticleGoogle Scholar - Nollo G, Marconcini M, Faes L, Bovolo F, Ravelli F, Bruzzone L:
**An automatic system for the analysis and classification of human atrial fibrillation patterns from intracardiac electrograms.***IEEE Trans Biomed Eng*2008,**55:**2275–2285.View ArticleGoogle Scholar - Faes L, Nollo G, Antolini R, Gaita F, Ravelli F:
**A method for quantifying atrial fibrillation organization based on wave-morphology similarity.***IEEE Trans Biomed Eng*2002,**49**(12 Pt 2):1504–1513.View ArticleGoogle Scholar - Theuns DA, Rivero-Ayerza M, Goedhart DM, Miltenburg M, Jordaens LJ:
**Morphology discrimination in implantable cardioverter-defibrillators: consistency of template match percentage during atrial tachyarrhythmias at different heart rates.***Europace*2008,**10:**1060–1066. 10.1093/europace/eun194View ArticleGoogle Scholar - Wu J, Estner H, Luik A, Ucer E, Reents T, Pflaumer A, Zrenner B, Hessling G, Deisenhofer I:
**Automatic 3D mapping of complex fractionated atrial electrograms (CFAE) in patients with paroxysmal and persistent atrial fibrillation.***J Cardiovasc Electrophysiol*2008,**19:**897–903. 10.1111/j.1540-8167.2008.01145.xView ArticleGoogle Scholar - Biviano AB, Coromilas J, Ciaccio EJ, Whang W, Hickey K, Garan H:
**Frequency domain and time complex analyses manifest low correlation and temporal variability when calculating activation rates in atrial fibrillation patients.***Pacing Clin Electrophysiol*2011,**34:**540–548. 10.1111/j.1540-8159.2010.02993.xView ArticleGoogle Scholar - Ciaccio EJ:
**Localization of the slow conduction zone during reentrant ventricular tachycardia.***Circulation*2000,**102:**464–469. 10.1161/01.CIR.102.4.464View ArticleGoogle Scholar - Ciaccio EJ, Coromilas J, Wit AL, Garan H:
**Onset dynamics of ventricular tachyarrhythmias as measured by dominant frequency.***Heart Rhythm*2011,**8:**615–623. 10.1016/j.hrthm.2010.11.028View ArticleGoogle Scholar - Ciaccio EJ, Tennyson CA, Bhagat G, Lewis SK, Green PH:
**Robust spectral analysis of videocapsule images acquired from celiac disease patients.***Biomed Eng Online*2011,**10:**78. 10.1186/1475-925X-10-78View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.