Software algorithm and hardware design for real-time implementation of new spectral estimator

BioMedical Engineering OnLine

Edward J Ciaccio²  Angelo B Biviano¹  Hasan Garan¹ 

[1] Department of Medicine–Division of Cardiology, Columbia University Medical Center, New York, USA

[2] Columbia University, Presbyterian Hospital 7 W-318, 630 West 168th Street, New York NY 10032, USA

关键词: Spectral analyzer;    Digital computer;    Circuit design;    Analog computer;    Algorithm;   

Others  :  793373
DOI  :  10.1186/1475-925X-13-61

 received in 2013-12-22, accepted in 2014-03-17,  发布年份 2014

【 摘 要 】

Background

Real-time spectral analyzers can be difficult to implement for PC computer-based systems because of the potential for high computational cost, and algorithm complexity. In this work a new spectral estimator (NSE) is developed for real-time analysis, and compared with the discrete Fourier transform (DFT).

Method

Clinical data in the form of 216 fractionated atrial electrogram sequences were used as inputs. The sample rate for acquisition was 977 Hz, or approximately 1 millisecond between digital samples. Real-time NSE power spectra were generated for 16,384 consecutive data points. The same data sequences were used for spectral calculation using a radix-2 implementation of the DFT. The NSE algorithm was also developed for implementation as a real-time spectral analyzer electronic circuit board.

Results

The average interval for a single real-time spectral calculation in software was 3.29 μs for NSE versus 504.5 μs for DFT. Thus for real-time spectral analysis, the NSE algorithm is approximately 150× faster than the DFT. Over a 1 millisecond sampling period, the NSE algorithm had the capability to spectrally analyze a maximum of 303 data channels, while the DFT algorithm could only analyze a single channel. Moreover, for the 8 second sequences, the NSE spectral resolution in the 3-12 Hz range was 0.037 Hz while the DFT spectral resolution was only 0.122 Hz. The NSE was also found to be implementable as a standalone spectral analyzer board using approximately 26 integrated circuits at a cost of approximately $500. The software files used for analysis are included as a supplement, please see the Additional files 1 and 2.

Conclusions

The NSE real-time algorithm has low computational cost and complexity, and is implementable in both software and hardware for 1 millisecond updates of multichannel spectra. The algorithm may be helpful to guide radiofrequency catheter ablation in real time.

【 授权许可】

   
2014 Ciaccio et al.; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

Figure 4.

Figure 5.

Figure 6.

【 参考文献 】

• [1]Bongiovanni G, Corsini P, Frosini G: Procedures for computing the discrete Fourier transform on staggered blocks. IEEE Transact Acoustics Speech Signal Process 1976, 24:132-137.
• [2]Lo PC, Lee YY: Real-time implementation of the moving FFT algorithm. Signal Process 1999, 79:251-259.
• [3]Nibouche O, Boussakta S, Darnell M, Benaissa M: Algorithms and pipeline architectures for 2-D FFT and FFT-like transforms. Digit Signal Process 2010, 20:1072-1086.
• [4]Huang HY, Lee YY, Lo PC: A novel algorithm for computing the 2D split-vector-radix FFT. Signal Process 2004, 84:561-570.
• [5]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.
• [6]Lin YJ, Tai CT, Kao T, Tso HW, Higa S, Tsao HM, Chang SL, Hsieh MH, Chen SA: Frequency analysis in different types of paroxysmal atrial fibrillation. J Am Coll Cardiol 2006, 47:1401-1407.
• [7]Sanders P, Berenfeld O, Hocini M, Jaïs P, Vaidyanathan R, Hsu LF, Garrigue S, Takahashi Y, Rotter M, Sacher F, Scavée C, Ploutz-Snyder R, Jalife J, Haïssaguerre M: Spectral analysis identifies sites of high-frequency activity maintaining atrial fibrillation in humans. Circulation 2005, 112:789-797.
• [8]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.
• [9]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. BioMed Central Full Text
• [10]Ciaccio EJ, Biviano AB, Garan H: Comparison of spectral estimators for characterizing fractionated atrial electrograms. BioMed Eng OnLine 2013, 12:72. BioMed Central Full Text
• [11]Ciaccio EJ, Biviano AB, Garan H: Computational method for high resolution spectral analysis of fractionated atrial electrograms. Comput Biol Med 2013, 43:1573-1582.
• [12]Holm M, Pehrson S, Ingemansson M, Sörnmo L, Johansson R, Sandhall L, Sunemark M, Smideberg B, Olsson C, Olsson SB: Non-invasive assessment of the atrial cycle length during atrial fibrillation in man: introducing, validating and illustrating a new ECG method. Cardiovasc Res 1998, 38:69-81.
• [13]Pehrson S, Holm M, Meurling C, Ingemansson M, Smideberg B, Sörnmo L, Olsson SB: Non-invasive assessment of magnitude and dispersion of atrial cycle length during chronic atrial fibrillation in man. Eur Heart J 1998, 19:1836-1844.
• [14]Jarman JW, Wong T, Kojodjojo P, Spohr H, Davies JE, Roughton M, Francis DP, Kanagaratnam P, Markides V, Davies DW, Peters NS: Spatiotemporal behavior of high dominant frequency during paroxysmal and persistent atrial fibrillation in the human left atrium. Circ Arrhythm Electrophysiol 2012, 5:650-658.
• [15]Salinet JL, Tuan JH, Sandilands AJ, Stafford PJ, Schlindwein FS, Ng GA: Distinctive patterns of dominant frequency trajectory behavior in drug-refractory persistent atrial fibrillation. J Cardiovasc Electrophysiol 2013. doi:10.1111/jce.12331
• [16]Press WH, Teukolsky SA, Vetterling WT, Flannery BP: Numerical Recipes in Fortran. NewYork: Cambridge University Press; 1992:501-502.
• [17]Ciaccio EJ, Biviano AB, Whang W, Gambhir A, Garan H: Spectral profiles of complex fractionated atrial electrograms are different in longstanding and acute onset atrial fibrillation atrial electrogram spectra. J Cardiovasc Electrophysiol 2012, 23:971-979.
• [18]Everett TH IV, Moorman JR, Kok LC, Akar JG, Haines DE: Assessment of global atrial fibrillation organization to optimize timing of atrial defibrillation. Circulation 2001, 103:2857-2861.
• [19]Everett TH IV, Verheule S, Wilson EE, Foreman S, Olgin JE: Left atrial dilatation resulting from chronic mitral regurgitation decreases spatiotemporal organization of atrial fibrillation in left atrium. Am J Physiol Heart Circ Physiol 2004, 286:H2452-H2460.
• [20]Mateo J, Rieta JJ: Radial basis function neural networks applied to efficient QRST cancellation in atrial fibrillation. Comput Biol Med 2013, 43:154-163.
• [21]benchFFT http://www.fftw.org/benchfft/ webcite
• [22]Ahmad A, Schlindwein FS, Ng GA: Comparison of computation time for estimation of dominant frequency of atrial electrograms: fast fourier transform, blackman tukey, autoregressive and multiple signal classification. J Biomed Sci Eng 2010, 3:843-847.
• [23]Nademanee K, McKenzie J, Kosar E, Schwab M, Sunsaneewitayakul B, Vasavakul T, Khunnawat C, Ngarmukos T: A new approach for catheter ablation of atrial fibrillation: mapping of the electrophysiologic substrate. J Am Coll Cardiol 2004, 43:2044-2053.
• [24]Nademanee K, Oketani N: The role of complex fractionated atrial electrograms in atrial fibrillation ablation. J Am Coll Cardiol 2009, 53:790-791.
• [25]Ciaccio EJ, Chow AW, Kaba RA, Davies DW, Segal OR, Peters NS: Detection of the diastolic pathway, circuit morphology, and inducibility of human postinfarction ventricular tachycardia from mapping in sinus rhythm. Heart Rhythm 2008, 5:981-991.
• [26]Narayan SM, Krummen DE, Rappel WJ: Clinical mapping approach to diagnose electrical rotors and focal impulse sources for human atrial fibrillation. J Cardiovasc Electrophysiol 2012, 23:447-454.
• [27]Atienza F, Almendral J, Jalife J, Zlochiver S, Ploutz-Snyder R, Torrecilla EG, Arenal A, Kalifa J, Fernández-Avilés F, Berenfeld O: Real-time dominant frequency mapping and ablation of dominant frequency sites in atrial fibrillation with left-to-right frequency gradients predicts long-term maintenance of sinus rhythm. Heart Rhythm 2009, 6:33-40.
• [28]Berenfeld O, Ennis S, Hwang E, Hooven B, Grzeda K, Mironov S, Yamazaki M, Kalifa J, Jalife J: Time-and frequency-domain analyses of atrial fibrillation activation rate: the optical mapping reference. Heart Rhythm 2011, 8:1758-1765.
• [29]Konings KT, Smeets JL, Penn OC, Wellens HJ, Allessie MA: Configuration of unipolar atrial electrograms during electrically induced atrial fibrillation in humans. Circulation 1997, 95:1231-1241.