The stationary wavelet transform as an efficient reductor of powerline interference for atrial bipolar electrograms in cardiac electrophysiology

[EN] Objective :The most relevant source of signal contamination in the cardiac electrophysiology (EP) laboratory is the ubiquitous powerline interference (PLI). To reduce this perturbation, algorithms including common fixed-bandwidth and adaptive-notch filters have been proposed. Although such methods have proven to add artificial fractionation to intra-atrial electrograms (EGMs), they are still frequently used. However, such morphological alteration can conceal the accurate interpretation of EGMs, specially to evaluate the mechanisms supporting atrial fibrillation (AF), which is the most common cardiac arrhythmia. Given the clinical relevance of AF, a novel algorithm aimed at reducing PLI on highly contaminated bipolar EGMs and, simultaneously, preserving their morphology is proposed. Approach: The method is based on the wavelet shrinkage and has been validated through customized indices on a set of synthesized EGMs to accurately quantify the achieved level of PLI reduction and signal morphology alteration. Visual validation of the algorithm¿s performance has also been included for some real EGM excerpts. Main results: The method has outperformed common filtering-based and wavelet-based strategies in the analyzed scenario. Moreover, it possesses advantages such as insensitivity to amplitude and frequency variations in the PLI, and the capability of joint removal of several interferences. Significance: The use of this algorithm in routine cardiac EP studies may enable improved and truthful evaluation of AF mechanisms.Research supported by grants DPI2017-83952-C3 MINECO/AEI/FEDER, UE and SBPLY/17/180501/000411 from Junta de Comunidades de Castilla-La Mancha.Martinez-Iniesta, M.; Ródenas, J.; Rieta, JJ.; Alcaraz, R. (2019). The stationary wavelet transform as an efficient reductor of powerline interference for atrial bipolar electrograms in cardiac electrophysiology. Physiological Measurement. 40(7):1-17. https://doi.org/10.1088/1361-6579/ab2cb8S117407Addison, P. S. (2017). The Illustrated Wavelet Transform Handbook. doi:10.1201/9781315372556Allen, D. P. (2009). A frequency domain Hampel filter for blind rejection of sinusoidal interference from electromyograms. Journal of Neuroscience Methods, 177(2), 303-310. doi:10.1016/j.jneumeth.2008.10.019Atienza, F., Almendral, J., Jalife, J., Zlochiver, S., Ploutz-Snyder, R., Torrecilla, E. G., … Berenfeld, O. (2009). 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, 6(1), 33-40. doi:10.1016/j.hrthm.2008.10.024Atienza, F., Almendral, J., Moreno, J., Vaidyanathan, R., Talkachou, A., Kalifa, J., … Berenfeld, O. (2006). Activation of Inward Rectifier Potassium Channels Accelerates Atrial Fibrillation in Humans. Circulation, 114(23), 2434-2442. doi:10.1161/circulationaha.106.633735Bahoura, M., & Ezzaidi, H. (2011). FPGA-Implementation of Discrete Wavelet Transform with Application to Signal Denoising. Circuits, Systems, and Signal Processing, 31(3), 987-1015. doi:10.1007/s00034-011-9355-0Baumert, M., Sanders, P., & Ganesan, A. (2016). Quantitative-Electrogram-Based Methods for Guiding Catheter Ablation in Atrial Fibrillation. Proceedings of the IEEE, 104(2), 416-431. doi:10.1109/jproc.2015.2505318Castillo, E., Morales, D. P., García, A., Martínez-Martí, F., Parrilla, L., & Palma, A. J. (2013). Noise Suppression in ECG Signals through Efficient One-Step Wavelet Processing Techniques. Journal of Applied Mathematics, 2013, 1-13. doi:10.1155/2013/763903Chen, S.-W., & Chen, Y.-H. (2015). Hardware Design and Implementation of a Wavelet De-Noising Procedure for Medical Signal Preprocessing. Sensors, 15(10), 26396-26414. doi:10.3390/s151026396El B’charri, O., Latif, R., Elmansouri, K., Abenaou, A., & Jenkal, W. (2017). ECG signal performance de-noising assessment based on threshold tuning of dual-tree wavelet transform. BioMedical Engineering OnLine, 16(1). doi:10.1186/s12938-017-0315-1Everett, T. H., Lai-Chow Kok, Vaughn, R. H., Moorman, R., & Haines, D. E. (2001). Frequency domain algorithm for quantifying atrial fibrillation organization to increase defibrillation efficacy. IEEE Transactions on Biomedical Engineering, 48(9), 969-978. doi:10.1109/10.942586Ganesan, A. N., Kuklik, P., Lau, D. H., Brooks, A. G., Baumert, M., Lim, W. W., … Sanders, P. (2013). Bipolar Electrogram Shannon Entropy at Sites of Rotational Activation. Circulation: Arrhythmia and Electrophysiology, 6(1), 48-57. doi:10.1161/circep.112.976654Goldberger, A. L., Amaral, L. A. N., Glass, L., Hausdorff, J. M., Ivanov, P. C., Mark, R. G., … Stanley, H. E. (2000). PhysioBank, PhysioToolkit, and PhysioNet. Circulation, 101(23). doi:10.1161/01.cir.101.23.e215GOLDBERGER, J. J. (2017). Substrate Ablation for Treatment of Atrial Fibrillation: Back to Basics. Journal of Cardiovascular Electrophysiology, 28(2), 156-158. doi:10.1111/jce.13149Gutiérrez-Gnecchi, J. A., Morfin-Magaña, R., Lorias-Espinoza, D., Tellez-Anguiano, A. del C., Reyes-Archundia, E., Méndez-Patiño, A., & Castañeda-Miranda, R. (2017). DSP-based arrhythmia classification using wavelet transform and probabilistic neural network. Biomedical Signal Processing and Control, 32, 44-56. doi:10.1016/j.bspc.2016.10.005Hamilton, P. S. (1996). A comparison of adaptive and nonadaptive filters for reduction of power line interference in the ECG. IEEE Transactions on Biomedical Engineering, 43(1), 105-109. doi:10.1109/10.477707Heijman, J., Algalarrondo, V., Voigt, N., Melka, J., Wehrens, X. H. T., Dobrev, D., & Nattel, S. (2015). The value of basic research insights into atrial fibrillation mechanisms as a guide to therapeutic innovation: a critical analysis. Cardiovascular Research, 109(4), 467-479. doi:10.1093/cvr/cvv275Honarbakhsh, S., Schilling, R. J., Orini, M., Providencia, R., Keating, E., Finlay, M., … Hunter, R. J. (2019). Structural remodeling and conduction velocity dynamics in the human left atrium: Relationship with reentrant mechanisms sustaining atrial fibrillation. Heart Rhythm, 16(1), 18-25. doi:10.1016/j.hrthm.2018.07.019Jadidi, A. S., Lehrmann, H., Weber, R., Park, C.-I., & Arentz, T. (2013). Optimizing Signal Acquisition and Recording in an Electrophysiology Laboratory. Cardiac Electrophysiology Clinics, 5(2), 137-142. doi:10.1016/j.ccep.2013.01.005JARMAN, J. W. E., WONG, T., KOJODJOJO, P., SPOHR, H., DAVIES, J. E. R., ROUGHTON, M., … PETERS, N. S. (2014). Organizational Index Mapping to Identify Focal Sources During Persistent Atrial Fibrillation. Journal of Cardiovascular Electrophysiology, 25(4), 355-363. doi:10.1111/jce.12352KOTTKAMP, H., BERG, J., BENDER, R., RIEGER, A., & SCHREIBER, D. (2015). Box Isolation of Fibrotic Areas (BIFA): A Patient-Tailored Substrate Modification Approach for Ablation of Atrial Fibrillation. Journal of Cardiovascular Electrophysiology, 27(1), 22-30. doi:10.1111/jce.12870Krijthe, B. P., Kunst, A., Benjamin, E. J., Lip, G. Y. H., Franco, O. H., Hofman, A., … Heeringa, J. (2013). Projections on the number of individuals with atrial fibrillation in the European Union, from 2000 to 2060. European Heart Journal, 34(35), 2746-2751. doi:10.1093/eurheartj/eht280Lau, D. H., Schotten, U., Mahajan, R., Antic, N. A., Hatem, S. N., Pathak, R. K., … Sanders, P. (2015). Novel mechanisms in the pathogenesis of atrial fibrillation: practical applications. European Heart Journal, 37(20), 1573-1581. doi:10.1093/eurheartj/ehv375Lenis, G., Pilia, N., Loewe, A., Schulze, W. H. W., & Dössel, O. (2017). Comparison of Baseline Wander Removal Techniques considering the Preservation of ST Changes in the Ischemic ECG: A Simulation Study. Computational and Mathematical Methods in Medicine, 2017, 1-13. doi:10.1155/2017/9295029Levkov, C., Mihov, G., Ivanov, R., Daskalov, I., Christov, I., & Dotsinsky, I. (2005). Removal of power-line interference from the ECG: a review of the subtraction procedure. BioMedical Engineering OnLine, 4(1). doi:10.1186/1475-925x-4-50Lian, J., Garner, G., Muessig, D., & Lang, V. (2010). A simple method to quantify the morphological similarity between signals. Signal Processing, 90(2), 684-688. doi:10.1016/j.sigpro.2009.07.010Liu, H., Li, Y., Zhou, Y., Jing, X., & Truong, T.-K. (2018). Joint power line interference suppression and ECG signal recovery in transform domains. Biomedical Signal Processing and Control, 44, 58-66. doi:10.1016/j.bspc.2018.04.001Mamun, M., Al-Kadi, M., & Marufuzzaman, M. (2013). Effectiveness of Wavelet Denoising on Electroencephalogram Signals. Journal of Applied Research and Technology, 11(1), 156-160. doi:10.1016/s1665-6423(13)71524-4Martens, S. M. M., Mischi, M., Oei, S. G., & Bergmans, J. W. M. (2006). An Improved Adaptive Power Line Interference Canceller for Electrocardiography. IEEE Transactions on Biomedical Engineering, 53(11), 2220-2231. doi:10.1109/tbme.2006.883631Mart�nez, M., R�denas, J., Alcaraz, R., & Rieta, J. J. (2017). Application of the Stationary Wavelet Transform to Reduce Power-line Interference in Atrial Electrograms. 2017 Computing in Cardiology Conference (CinC). doi:10.22489/cinc.2017.112-049Martínez-Iniesta, M., Ródenas, J., Alcaraz, R., & Rieta, J. J. (2017). Waveform Integrity in Atrial Fibrillation: The Forgotten Issue of Cardiac Electrophysiology. Annals of Biomedical Engineering, 45(8), 1890-1907. doi:10.1007/s10439-017-1832-6Metting van Rijn, A. C., Peper, A., & Grimbergen, C. A. (1990). High-quality recording of bioelectric events. Medical & Biological Engineering & Computing, 28(5), 389-397. doi:10.1007/bf02441961Nademanee, K., McKenzie, J., Kosar, E., Schwab, M., Sunsaneewitayakul, B., Vasavakul, T., … Ngarmukos, T. (2004). A new approach for catheter ablation of atrial fibrillation: mapping of the electrophysiologic substrate. Journal of the American College of Cardiology, 43(11), 2044-2053. doi:10.1016/j.jacc.2003.12.054Naseri, H., & Homaeinezhad, M. R. (2012). Computerized quality assessment of phonocardiogram signal measurement-acquisition parameters. Journal of Medical Engineering & Technology, 36(6), 308-318. doi:10.3109/03091902.2012.684832NG, J., BORODYANSKIY, A. I., CHANG, E. T., VILLUENDAS, R., DIBS, S., KADISH, A. H., & GOLDBERGER, J. J. (2010). Measuring the Complexity of Atrial Fibrillation Electrograms. Journal of Cardiovascular Electrophysiology, 21(6), 649-655. doi:10.1111/j.1540-8167.2009.01695.xNg, J., Sehgal, V., Ng, J. K., Gordon, D., & Goldberger, J. J. (2014). Iterative Method to Detect Atrial Activations and Measure Cycle Length From Electrograms During Atrial Fibrillation. IEEE Transactions on Biomedical Engineering, 61(2), 273-278. doi:10.1109/tbme.2013.2290003Oesterlein, T. G., Lenis, G., Rudolph, D.-T., Luik, A., Verma, B., Schmitt, C., & Dössel, O. (2015). Removing ventricular far-field signals in intracardiac electrograms during stable atrial tachycardia using the periodic component analysis. Journal of Electrocardiology, 48(2), 171-180. doi:10.1016/j.jelectrocard.2014.12.004Platonov, P. G., Corino, V. D. A., Seifert, M., Holmqvist, F., & Sornmo, L. (2014). Atrial fibrillatory rate in the clinical context: natural course and prediction of intervention outcome. Europace, 16(suppl 4), iv110-iv119. doi:10.1093/europace/euu249Poornachandra, S., & Kumaravel, N. (2008). A novel method for the elimination of power line frequency in ECG signal using hyper shrinkage function. Digital Signal Processing, 18(2), 116-126. doi:10.1016/j.dsp.2007.03.011Quintanilla, J. G., Pérez-Villacastín, J., Pérez-Castellano, N., Pandit, S. V., Berenfeld, O., Jalife, J., & Filgueiras-Rama, D. (2016). Mechanistic Approaches to Detect, Target, and Ablate the Drivers of Atrial Fibrillation. Circulation: Arrhythmia and Electrophysiology, 9(1). doi:10.1161/circep.115.002481Ravelli, F., Masè, M., Cristoforetti, A., Marini, M., & Disertori, M. (2014). The logical operator map identifies novel candidate markers for critical sites in patients with atrial fibrillation. Progress in Biophysics and Molecular Biology, 115(2-3), 186-197. doi:10.1016/j.pbiomolbio.2014.07.006Schanze, T. (2017). Removing noise in biomedical signal recordings by singular value decomposition. Current Directions in Biomedical Engineering, 3(2), 253-256. doi:10.1515/cdbme-2017-0052Schotten, U., Dobrev, D., Platonov, P. G., Kottkamp, H., & Hindricks, G. (2016). Current controversies in determining the main mechanisms of atrial fibrillation. Journal of Internal Medicine, 279(5), 428-438. doi:10.1111/joim.12492Singh, B. N., & Tiwari, A. K. (2006). Optimal selection of wavelet basis function applied to ECG signal denoising. Digital Signal Processing, 16(3), 275-287. doi:10.1016/j.dsp.2005.12.003Van Wagoner, D. R., Piccini, J. P., Albert, C. M., Anderson, M. E., Benjamin, E. J., Brundel, B., … Wehrens, X. H. T. (2015). Progress toward the prevention and treatment of atrial fibrillation: A summary of the Heart Rhythm Society Research Forum on the Treatment and Prevention of Atrial Fibrillation, Washington, DC, December 9–10, 2013. Heart Rhythm, 12(1), e5-e29. doi:10.1016/j.hrthm.2014.11.011Venkatachalam, K. L., Herbrandson, J. E., & Asirvatham, S. J. (2011). Signals and Signal Processing for the Electrophysiologist. Circulation: Arrhythmia and Electrophysiology, 4(6), 965-973. doi:10.1161/circep.111.96430

Analysis of Atrial Electrograms

This work provides methods to measure and analyze features of atrial electrograms - especially complex fractionated atrial electrograms (CFAEs) - mathematically. Automated classification of CFAEs into clinical meaningful classes is applied and the newly gained electrogram information is visualized on patient specific 3D models of the atria. Clinical applications of the presented methods showed that quantitative measures of CFAEs reveal beneficial information about the underlying arrhythmia

Signal Processing Methods for the Analysis of the Electrocardiogram

Das Elektrokardiogramm (EKG) zeichnet die elektrische Aktivität des Herzens auf der Brust- oberfläche auf. Dieses Signal kann einfach und kostengünstig aufgenommen werden und wird daher in einer Vielzahl von mobilen und stationären Anwendungen genutzt. Es ist über die letzten 100 Jahre zum Goldstandard bei der Diagnose vieler kardiologischer Krankheiten geworden. Herzerkrankungen bleiben ein relevantes Thema in unserer Gesellschaft, da sie zu 30 % aller Todesfälle weltweit führen. Allein die koronare Herzkrankheit ist die häufigste Todesursache überhaupt. Weiterhin sind 2 bis 3 % der Europäer von Herzrhythmusstörungen wie Vorhofflimmern und Vorhofflattern betroffen. Die damit verbundenen geschätzten Kosten in der Europäischen Union belaufen sich auf 26 Milliarden Euro pro Jahr. In allen diesen Fällen ist die Aufzeichnung des EKGs der erste unumgängliche Schritt für eine verlässliche Diagnose und erfolgreiche Therapie. Im Rahmen dieser Dissertation wurden eine Reihe von Algorithmen zur Signalver- arbeitung des EKG entwickelt, die automatisch die rhythmischen und morphologischen Eigenschaften aus dem EKG extrahieren und dadurch den diagnostischen Prozess und die Entscheidungsfindung des Arztes unterstützen. In einem ersten Projekt wurde das Phänomen der postextrasystolischen T-Wellen-Änderung (PEST) untersucht. Die aus der PEST ex- trahierten Biomarker haben wir als Prädiktoren für Herzversagen postuliert. Ein zweites Projekt handelte vom Entwurf eines akkuraten Algorithmus zur Detektion und Annotation der P-Welle im EKG. Als Referenz während der Entwicklung wurden intrakardial gemessene Signale verwendet. Eine dritte Untersuchg hatte das Ziel, das physiologische Phänomen der respiratorischen Sinusarrhythmie (RSA) besser zu verstehen. In diesem Projekt wurde ein Algorithmus zur Trennung der Herzratenvariabilität (HRV) in ihre atmungsabhängige und ihre atmungsunabhn ̈gige Komponente untersucht. Letzterer Anteil der HRV könnte neue Erkenntnisse über die Regulationsmechanismen des kardiovaskulären Systems liefern. In der vierten und letzten Studie wurde der Einfluss mentaler Belastung auf das EKG während der Autofahrt untersucht. Eine Vielzahl von Deskriptoren wurden gefunden, die eine gefährliche mentale Beanspruchung detektieren und somit den Fahrer vor einem möglichen Unfall schützen können. Wir schließen aus diesen Untersuchungen, dass gut entwickelte Methoden der Signalver- arbeitung des EKG das Potential haben, die Belastung der Patienten, die an Herzerkrankungen leiden, und die Anzahl der Verkehrsunfälle zu reduzieren

A new transform for the analysis of complex fractionated atrial electrograms

Representation of independent biophysical sources using Fourier analysis can be inefficient because the basis is sinusoidal and general. When complex fractionated atrial electrograms (CFAE) are acquired during atrial fibrillation (AF), the electrogram morphology depends on the mix of distinct nonsinusoidal generators. Identification of these generators using efficient methods of representation and comparison would be useful for targeting catheter ablation sites to prevent arrhythmia reinduction. A data-driven basis and transform is described which utilizes the ensemble average of signal segments to identify and distinguish CFAE morphologic components and frequencies. Calculation of the dominant frequency (DF) of actual CFAE, and identification of simulated independent generator frequencies and morphologies embedded in CFAE, is done using a total of 216 recordings from 10 paroxysmal and 10 persistent AF patients. The transform is tested versus Fourier analysis to detect spectral components in the presence of phase noise and interference. Correspondence is shown between ensemble basis vectors of highest power and corresponding synthetic drivers embedded in CFAE. Results: The ensemble basis is orthogonal, and efficient for representation of CFAE components as compared with Fourier analysis (p ≤ 0.002). When three synthetic drivers with additive phase noise and interference were decomposed, the top three peaks in the ensemble power spectrum corresponded to the driver frequencies more closely as compared with top Fourier power spectrum peaks (p ≤ 0.005). The synthesized drivers with phase noise and interference were extractable from their corresponding ensemble basis with a mean error of less than 10%. Conclusions: The new transform is able to efficiently identify CFAE features using DF calculation and by discerning morphologic differences. Unlike the Fourier transform method, it does not distort CFAE signals prior to analysis, and is relatively robust to jitter in periodic events. Thus the ensemble method can provide a useful alternative for quantitative characterization of CFAE during clinical study

High-Density Mapping Analysis of Electrical Spatiotemporal Behaviour in Atrial Fibrillation

Tese de mestrado integrado, Engenharia Biomédica e Biofísica (Sinais e Imagens Médicas), 2022, Universidade de Lisboa, Faculdade de CiênciasDoenças cardiovasculares, tais como arritmias, são a principal causa de morte no mundo, especialmente no Sul e no Este da Ásia, e nos Estados Unidos da América [1]. As arritmas são caracterizadas pela alteração no ritmo sinusal normal do coração. Em particular, a fibrilhação auricular (FA) é a arritmia cardíaca mais comum na prática clínica, contribuindo para mais de 200 mil mortes globalmente em 2017 [2]. Caracteriza-se pela contração rápida e dessincronizada das aurículas, e está associada ao aumento da mortalidade e afecta de forma negativa a qualidade de vida dos pacientes. A FA é geralmente tratada através de medicação, porém quando esta falha, a ablação por cateter é indicada, sendo um tratamento de referência para combater esta patologia. A ablação apresenta uma taxa de sucesso de aproximadamente 50% no primeiro procedimento, sendo necessário efectuar vários procedimentos para aumentar a eficácia do tratamento [3]. A detecção desta patologia envolve, numa primeira fase, a realização de um electrocardiograma (ECG) e, posteriormente um estudo electrofisiológico para saber com precisão onde se localiza e o mecanismo subjacente à mesma. Este último implica o registo da actividade eléctrica através de electrogramas (EGM) locais em diferentes pontos das aurículas e dos ventrículos, com o auxílio de sistemas de mapeamento tridimensionais (3D) electroanatómicos, sendo um procedimento invasivo. Existem diversos métodos lineares e não lineares que permitem a análise dos EGMs nos domínios do tempo, frequência, fase, entre outros, com a finalidade de melhor compreender os mecanismos subjacentes à FA e, consequentemente aumentar a taxa de sucesso do processo de ablação e melhorar a sua eficiência. Esta área de estudo progrediu significativamente, tanto a nível de hardware, como de software. Apesar disso, os métodos desenvolvidos não têm nem acrescentado benefícios adicionais, nem melhorado significativamente a taxa de sucesso do processo de ablação. Existem várias razões para tal, e grande parte deve-se ao facto destes métodos de análise estarem incorporados nos sistemas de mapeamento e o seu software ser exclusivo. Isto leva a que não consigamos perceber como é que os algoritmos funcionam nos diferentes sistemas de mapeamento para comparar as suas diferenças e semelhanças. Devido a estes constrangimentos, os investigadores são compelidos a desenvolver os seus próprios métodos de análise e técnicas de mapeamento, o que leva à existência de uma multitude de métodos e técnicas de mapeamento que parecem ser diferentes entre si, resultando em informação ambígua e conflituosa no que diz respeito aos mecanismos da FA, e a conclusões distintas entre estudos. O sucesso do tratamento poderia aumentar se tivéssemos uma melhor compreensão dos métodos de análise e da sua aplicação no contexto da FA; perceber se os métodos apontam para o mesmo fenómeno de fibrilhação, se existe alguma correlação entre os métodos, e se a informação fornecida pelos mesmos é complementar ou redundante. Assim, o objectivo deste trabalho consistiu em implementar diferentes métodos para analisar os EGMs e a estrutura 3D da aurícula esquerda (AE) de doentes com FA, numa tentativa de responder às questões que motivaram a realização deste projecto. Em última análise, ao observar os mapas 3D da AE tendo uma melhor compreensão dos métodos, poderemos identificar com precisão as regiões na AE responsáveis por iniciar a FA, e ter mais conhecimento sobre os mecanismos responsáveis pela mesma. Desta forma, o processo de ablação poderá alcançar o seu potencial. Para este projecto, foram incluídos os mapas 3D electroanatómicos da AE de dez doentes com FA paroxística ou persistente do hospital de Santa Marta, recolhidos com o sistema de mapeamento CARTO 3. Cada ponto electroanatómico dos mapas inclui as 12 derivações do ECG, e os EGMs unipolares e bipolares registados com o cateter de mapeamento Pentaray de 20 pólos. Porém, apenas os EGMs bipolares foram incluídos na análise. Processaram-se os sinais bipolares e, devido a algumas limitações, foi possível apenas a implementação de dois métodos diferentes para os analisar: um no domínio da frequência – Frequência Dominante (FD) –, e outro no domínio da Teoria da Informação – a entropia de Shannon. De seguida, criaram-se três tipos de mapas 3D electroanatómicos da AE para cada doente: um de voltagem, cuja informação foi adquirida com o sistema de mapeamento, um de FD, e outro de entropia. A informação de cada mapa estava organizada segundo um padrão de cores. Observando os diferentes tipos de mapas da AE paralelamente, foi possível comparar os métodos, e perceber que tipo de informação cada um deles fornecia, numa tentativa de melhor compreender os mecanismos da FA. Foi possível observar em algumas regiões da AE, principalmente nos mapas de voltagem e de FD, a presença de “centros de activação” ou “centros de fibrilhação”, que poderão ser os gatilhos responsáveis por desencadear ou manter o mecanismo de fibrilhação. Para confirmar se de facto aquelas regiões eram os gatilhos de fibrilhação, seria necessário submeter os doentes ao processo de ablação e queimar essas zonas; e posteriormente acompanhar os doentes para observar os efeitos do procedimento e confirmar a hipótese. Contudo, dadas as limitações do trabalho e o facto desta área de investigação ser pouco explorada, é fulcral obter um maior número de estudo comparativos entre mais métodos de diferentes domínios e confirmar se apontam ou não para o mesmo fenómeno de fibrilhação. Apesar de terem sido implementados apenas dois métodos de análise dos EGMs, o projecto permitiu a comparação entre os mesmos, uma área de estudo por onde ainda há muito para investigar. Com mais conhecimento sobre os diferentes métodos, a sua aplicação, inter-relação e adequação no estudo dos mecanismos da FA e das propriedades electrofisiológicas desta patologia, é possível desenvolver procedimentos de ablação mais eficientes e selectivos, de forma a diminuir os riscos e aumentar a taxa de sucesso do tratamento.Atrial fibrillation (AF) is the most frequent cardiac arrhythmia in clinical practice and is described by rapid and irregular contractions of the atria. Despite catheter ablation (CA) being a well-established treatment for AF, it is sub-optimal, with a success rate of approximately 50 % after a single procedure, with some patients requiring multiple procedures to achieve long-term freedom from this pathology. This prompted the proposal and development of various quantitative electrogram (EGM)-based methods along with different mapping systems with their respective mapping techniques, to better understand the mechanisms responsible for initiating and maintaining AF, thus improving ablation outcomes. However, this diversification of methods and tools resulted in disperse and inconsistent data regarding the mechanisms of AF. This work consisted of employing two different methods to analyse the electrograms (EGM): dominant frequency (DF) and Shannon entropy (ShEn). From these EGMs, metrics were then extracted and displayed in colour-coded fashion on a 3D mesh of the left atrium (LA) from patients with paroxysmal or persistent AF. The two methods were compared to understand whether or not these indicated different phenomena/mechanisms, and if these could locate sites suspected of triggering and maintaining AF. The results, while not fully conforming to the literature, allowed the comparison between different EGM analysis methods, a field of study that requires further research. Overall, this project highlighted the limited data available within the topic, hindering our understanding of AF mechanisms and development of more effective and selective ablation procedures to avoid unnecessary complications, and ultimately improve the effects of the treatment's outcomes

Study on the non-linear metrics contribution to estimate atrial fibrillation organization from the surface electrocardiogram

[EN] Atrial fibrillation (AF) is the most frequently diagnosed arrhythmia, characterized by an uncoordinated atrial electrical activation, thus causing the atria to be unable to pump blood effectively. The prevalence of AF is expected to increase significantly in the next decades as the population ages. However, both the knowledge and the treatment of this arrhythmia still have to experiment a significant progress. Previous studies have reported that AF organization, which can be defined as the repetitiveness degree of the atrial activity pattern, correlates with the arrhythmia status as well as with the therapy outcome. Thus, estimating AF organization from surface electrocardiographic (ECG) recordings constitutes a very interesting approach because ECG recordings are easy and cheap to obtain. The objective of this doctoral thesis is to assess the use of a variety of nonlinear indices in the estimation of AF organization from single-lead noninvasive ECG recordings. Apart from the most common noninvasive AF organization estimators, such as Sample Entropy (SampEn) and the dominant atrial frequency (DAF), the following nonlinear indices have been studied: Fuzzy Entropy, Spectral Entropy, Lempel-Ziv Complexity and Hurst Exponents. Moreover, since the presence of noise and ventricular residuals affects the performance of nonlinear methods, the application of a strategy aimed at reducing these nuisances has been evaluated. Therefore, the application of these metrics over the atrial activity fundamental waveform, named the main atrial wave (MAW), has been proposed. In this doctoral thesis, the following scenarios involving AF organization have been considered: the prediction of paroxysmal AF spontaneous termination, the study of the earlier signs anticipating AF termination and the classification between paroxysmal and persistent AF from short ECG recordings. Firstly, the performance of the studied metrics discriminating events related to AF organization was tested making use of a reference database aimed at predicting AF spontaneous termination. In this study, most of the proposed indices provided higher accuracy than traditional AF organization estimators. Accuracy values higher than 90% were obtained with several indices. In particular, the generalized Hurst exponents of order 1 and 2, H(1) and H(2), achieved outstanding results, thus being selected for later studies in this thesis. Furthermore, the computation of H(2) depends on two critical parameters, namely, the analyzed interval length (L) and the maximum search window for self-similarities (tau). Hence, a study with 660 combinations on these two parameters was performed, together with the sampling frequency (fs) of the recording, in order to obtain their optimal combination in computing AF organization. On the other hand, previous works analyzing the spontaneous termination of AF have been only focused on the last 2 minutes preceding the termination. In contrast, a different scenario considering longer recordings to detect the earlier signs anticipating paroxysmal AF termination has been analyzed for the first time in this thesis. H(2) was selected for the study because of its highest accuracy in AF termination prediction. Additionally, the DAF and SampEn were also computed as references. Through this study it has been corroborated that AF organization only varies significantly within the last 3 minutes before spontaneous termination. As a consequence, the early prediction of paroxysmal AF spontaneous termination does not seem feasible through the current signal analysis tools. Finally, H(2) was applied in the classification between paroxysmal and persistent AF from short ECG recordings, achieving a higher diagnostic accuracy than DAF and SampEn. This result suggests that the analysis of ambulatory ECG recordings through H(2) could be a future alternative to the use of Holter ECG recordings in the classification between paroxysmal and persistent AF.[ES] La fibrilación auricular (FA) es la arritmia más frecuente y se caracteriza por una actividad auricular descoordinada, que impide que las aurículas bombeen sangre de manera eficaz. Se espera que la prevalencia de la FA aumente significativamente en las próximas décadas debido al envejecimiento de la población. Sin embargo, tanto el conocimiento relativo a esta arritmia como su tratamiento son todavía mejorables. Estudios previos han relacionado la organización de la FA, que se puede definir como el grado de repetitividad de la actividad auricular, con el estado de la arritmia o su respuesta al tratamiento. Además, la estimación de la organización de la FA a partir de registros electrocardiográficos (ECG) de superficie resulta especialmente interesante porque su obtención es sencilla y barata. El objetivo de esta tesis doctoral es evaluar el uso de distintos índices no lineales para estimar la organización de la FA a partir del ECG. Además de los estimadores no invasivos de organización más comunes, como la entropía muestral (SampEn) y la frecuencia auricular dominante (DAF), se han estudiado los siguientes métodos no lineales: la entropía borrosa, la entropía espectral, la complejidad Lempel-Ziv y los exponentes de Hurst. Además, se ha estudiado el uso de una estrategia destinada a la reducción del ruido y los residuos de actividad ventricular para mejorar el desempeño de métodos no lineales. Así, los índices estudiados también se han aplicado sobre la forma de onda fundamental de la actividad auricular, conocida como la onda auricular principal (MAW). Se han considerado los siguientes escenarios relacionados con la organización de la FA: la predicción de la terminación espontánea de la FA paroxística, el estudio de los primeros indicios de terminación espontánea de la FA y la clasificación entre FA paroxística y FA persistente a partir de registros ECG de corta duración. Primero, se estudió la capacidad de los índices estudiados para distinguir eventos relacionados con la organización de la FA mediante el análisis de una base de datos de referencia para la predicción de su terminación espontánea. La mayoría de los índices propuestos consiguieron una mayor precisión que los estimadores tradicionales de organización. Así, varios de los índices obtuvieron una precisión superior al 90% en la predicción de la terminación espontánea de la FA. En particular, los exponentes de Hurst generalizados de orden 1 y 2, H(1) y H(2), lograron los mejores resultados de clasificación. Puesto que el cálculo de H(2) depende de dos parámetros críticos, la longitud del intervalo analizado (L) y el tamaño máximo de la ventana donde buscar similitudes (tau), se llevó a cabo un estudio con 660 combinaciones de esos dos parámetros junto con la frecuencia de muestreo (fs) del registro para determinar el uso óptimo de este índice. Por otra parte, los trabajos previos que han estudiado la terminación espontánea de la FA se han centrado en los últimos 2 minutos antes de la terminación. Por contra, en esta tesis doctoral se han estudiado por primera vez registros de mayor duración para detectar los primeros indicios de la terminación de la FA. Para ello, se eligió el uso de H(2) por su alta precisión en la predicción de la terminación de la FA. Además, la DAF y SampEn se calcularon como referencias. En este estudio se ha comprobado que la organización de la FA solamente presenta variaciones significativas en los últimos 3 minutos antes de su terminación espontánea. Por ello, la predicción temprana de la terminación no parece posible con los medios actuales de análisis de la señal. Por último, se aplicó H(2) para clasificar entre FA paroxística y FA persistente a partir de ECGs de corta duración, obteniendo una mayor precisión diagnóstica que la DAF y SampEn. Este resultado sugiere que el análisis de ECGs ambulatorios por medio de H(2) puede ser en el futuro una alte[CA] La fibril·lació auricular (FA) és l'arítmia més freqüent i es caracteritza per una activitat auricular descoordinada, que impedix que les aurícules bomben sang de manera eficaç. S'espera que la prevalença de la FA augmente significativament en les pròximes dècades a causa de l'envelliment de la població. No obstant això, tant el coneixement relatiu a esta arítmia com el seu tractament són encara millorables. Estudis previs han relacionat l'organització de la FA, que es pot definir com el grau de repetitivitat de l'activitat auricular, amb l'estat de l'arítmia o la seua resposta al tractament. A més, l'estimació de l'organització de la FA a partir de registres electrocardiogràfics (ECG) de superfície resulta especialment interessant perquè la seua obtenció és senzilla i barata. L'objectiu d'esta tesi doctoral és avaluar l'ús de distints índexs no lineals en l'estimació de l'organització de la FA a partir de l'ECG de superfície. A més dels estimadors no invasius d'organització més comuns, com l'entropia mostral (SampEn) i la freqüència auricular dominant (DAF), s'han estudiat els següents mètodes no lineals: l'entropia borrosa, l'entropia espectral, la complexitat Lempel-Ziv i els exponents de Hurst. A més, s'ha estudiat l'ús d'una estratègia destinada a la reducció del soroll i els residus d'activitat ventricular per a millorar la seua capacitat d'estimar l'organització. Així, doncs, els índexs estudiats també s'han aplicat sobre la forma d'onda fonamental de l'activitat auricular, coneguda com l'onda auricular principal (MAW). S'han considerat els següents escenaris relacionats amb l'organització de la FA: la predicció de la terminació espontània de la FA paroxística, l'estudi dels primers indicis de terminació espontània de la FA i la classificació entre FA paroxística i FA persistent a partir de registres ECG de curta duració. Primer, es va estudiar la capacitat dels índexs estudiats per a distingir esdeveniments relacionats amb l'organització de la FA per mitjà de l'anàlisi d'una base de dades de referència per a la predicció de la seua terminació espontània. La majoria dels índexs proposats van aconseguir una major precisió que els estimadors tradicionals d'organització de la FA. Així, alguns dels índexs van obtindre una precisió superior al 90% en la predicció de la terminació espontània de la FA. En particular, els exponents de Hurst generalitzats d'orde 1 i 2, H(1) i H(2), van aconseguir els millors resultats de classificació. Com el càlcul de H(2) depén de dos paràmetres crítics, la longitud de l'interval analitzat (L) i la grandària màxima de la finestra on buscar similituds (tau), es va dur a terme un estudi amb 660 combinacions d'eixos dos paràmetres junt amb la freqüència de mostratge (fs) del registre per a determinar la combinació òptima de valors per a estimar l'organització de la FA. D'altra banda, els treballs previs que han estudiat la terminació espontània de la FA s'han centrat en els últims 2 minuts abans de la terminació. Per contra, en esta tesi doctoral s'han estudiat per primera vegada registres de major duració amb l'objectiu de detectar els primers indicis de la terminació de la FA. Es va triar l'ús de H(2) per a este estudi per la seua alta precisió en la predicció de la terminació de la FA. A més, la DAF i SampEn es van calcular com a referències. En este estudi s'ha comprovat que l'organització de la FA només presenta variacions significatives en els últims 3 minuts abans de la seua terminació espontània. Per això, la predicció primerenca de la terminació no pareix possible amb els mitjans actuals d'anàlisi del senyal. Finalment, es va aplicar H(2) per a classificar entre FA paroxística i FA persistent a partir d'ECGs de curta duració, obtenint una millor precisió diagnòstica que amb la DAF i SampEn. Este resultat suggerix que l'anàlisi d'ECGs ambulatoris per mitjà de H(2) pot ser en eJulián Seguí, M. (2015). Study on the non-linear metrics contribution to estimate atrial fibrillation organization from the surface electrocardiogram [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/56150TESI

Classification of De novo post-operative and persistent atrial fibrillation using multi-channel ECG recordings

Atrial fibrillation (AF) is the most sustained arrhythmia in the heart and also the most common complication developed after cardiac surgery. Due to its progressive nature, timely detection of AF is important. Currently, physicians use a surface electrocardiogram (ECG) for AF diagnosis. However, when the patient develops AF, its various development stages are not distinguishable for cardiologists based on visual inspection of the surface ECG signals. Therefore, severity detection of AF could start from differentiating between short-lasting AF and long-lasting AF. Here, de novo post-operative AF (POAF) is a good model for short-lasting AF while long-lasting AF can be represented by persistent AF. Therefore, we address in this paper a binary severity detection of AF for two specific types of AF. We focus on the differentiation of these two types as de novo POAF is the first time that a patient develops AF. Hence, comparing its development to a more severe stage of AF (e.g., persistent AF) could be beneficial in unveiling the electrical changes in the atrium. To the best of our knowledge, this is the first paper that aims to differentiate these different AF stages. We propose a method that consists of three sets of discriminative features based on fundamentally different aspects of the multi-channel ECG data, namely based on the analysis of RR intervals, a greyscale image representation of the vectorcardiogram, and the frequency domain representation of the ECG. Due to the nature of AF, these features are able to capture both morphological and rhythmic changes in the ECGs. Our classification system consists of a random forest classifier, after a feature selection stage using the ReliefF method. The detection efficiency is tested on 151 patients using 5-fold cross-validation. We achieved 89.07% accuracy in the classification of de novo POAF and persistent AF. The results show that the features are discriminative to reveal the severity of AF. Moreover, inspection of the most important features sheds light on the different characteristics of de novo post-operative and persistent AF.</p

Statistical and Graph-Based Signal Processing: Fundamental Results and Application to Cardiac Electrophysiology

The goal of cardiac electrophysiology is to obtain information about the mechanism, function, and performance of the electrical activities of the heart, the identification of deviation from normal pattern and the design of treatments. Offering a better insight into cardiac arrhythmias comprehension and management, signal processing can help the physician to enhance the treatment strategies, in particular in case of atrial fibrillation (AF), a very common atrial arrhythmia which is associated to significant morbidities, such as increased risk of mortality, heart failure, and thromboembolic events. Catheter ablation of AF is a therapeutic technique which uses radiofrequency energy to destroy atrial tissue involved in the arrhythmia sustenance, typically aiming at the electrical disconnection of the of the pulmonary veins triggers. However, recurrence rate is still very high, showing that the very complex and heterogeneous nature of AF still represents a challenging problem. Leveraging the tools of non-stationary and statistical signal processing, the first part of our work has a twofold focus: firstly, we compare the performance of two different ablation technologies, based on contact force sensing or remote magnetic controlled, using signal-based criteria as surrogates for lesion assessment. Furthermore, we investigate the role of ablation parameters in lesion formation using the late-gadolinium enhanced magnetic resonance imaging. Secondly, we hypothesized that in human atria the frequency content of the bipolar signal is directly related to the local conduction velocity (CV), a key parameter characterizing the substrate abnormality and influencing atrial arrhythmias. Comparing the degree of spectral compression among signals recorded at different points of the endocardial surface in response to decreasing pacing rate, our experimental data demonstrate a significant correlation between CV and the corresponding spectral centroids. However, complex spatio-temporal propagation pattern characterizing AF spurred the need for new signals acquisition and processing methods. Multi-electrode catheters allow whole-chamber panoramic mapping of electrical activity but produce an amount of data which need to be preprocessed and analyzed to provide clinically relevant support to the physician. Graph signal processing has shown its potential on a variety of applications involving high-dimensional data on irregular domains and complex network. Nevertheless, though state-of-the-art graph-based methods have been successful for many tasks, so far they predominantly ignore the time-dimension of data. To address this shortcoming, in the second part of this dissertation, we put forth a Time-Vertex Signal Processing Framework, as a particular case of the multi-dimensional graph signal processing. Linking together the time-domain signal processing techniques with the tools of GSP, the Time-Vertex Signal Processing facilitates the analysis of graph structured data which also evolve in time. We motivate our framework leveraging the notion of partial differential equations on graphs. We introduce joint operators, such as time-vertex localization and we present a novel approach to significantly improve the accuracy of fast joint filtering. We also illustrate how to build time-vertex dictionaries, providing conditions for efficient invertibility and examples of constructions. The experimental results on a variety of datasets suggest that the proposed tools can bring significant benefits in various signal processing and learning tasks involving time-series on graphs. We close the gap between the two parts illustrating the application of graph and time-vertex signal processing to the challenging case of multi-channels intracardiac signals

Novel approaches for quantitative electrogram analysis for rotor identification: Implications for ablation in patients with atrial fibrillation

University of Minnesota Ph.D. dissertation. May 2017. Major: Biomedical Engineering. Advisor: Elena Tolkacheva. 1 computer file (PDF); xxviii, 349 pages + 4 audio/video filesAtrial fibrillation (AF) is the most common sustained cardiac arrhythmia that causes stroke affecting more than 2.3 million people in the US. Catheter ablation with pulmonary vein isolation (PVI) to terminate AF is successful for paroxysmal AF but suffers limitations with persistent AF patients as current mapping methods cannot identify AF active substrates outside of PVI region. Recent evidences in the mechanistic understating of AF pathophysiology suggest that ectopic activity, localized re-entrant circuit with fibrillatory propagation and multiple circuit re-entries may all be involved in human AF. Accordingly, the hypothesis that rotor is an underlying AF mechanism is compatible with both the presence of focal discharges and multiple wavelets. Rotors are stable electrical sources which have characteristic spiral waves like appearance with a pivot point surrounded by peripheral region. Targeted ablation at the rotor pivot points in several animal studies have demonstrated efficacy in terminating AF. The objective of this dissertation was to develop robust spatiotemporal mapping techniques that can fully capture the intrinsic dynamics of the non-stationary time series intracardiac electrogram signal to accurately identify the rotor pivot zones that may cause and maintain AF. In this thesis, four time domain approaches namely multiscale entropy (MSE) recurrence period density entropy (RPDE), kurtosis and intrinsic mode function (IMF) complexity index and one frequency domain approach namely multiscale frequency (MSF) was proposed and developed for accurate identification of rotor pivot points. The novel approaches were validated using optical mapping data with induced ventricular arrhythmia in ex-vivo isolated rabbit heart with single, double and meandering rotors (including numerically simulated data). The results demonstrated the efficacy of the novel approaches in accurate identification of rotor pivot point. The chaotic nature of rotor pivot point resulted in higher complexity measured by MSE, RPDE, kurtosis, IMF and MSF compared to the stable rotor periphery that enabled its accurate identification. Additionally, the feasibility of using conventional catheter mapping system to generate patient specific 3D maps for intraprocedural guidance for catheter ablation using these novel approaches was demonstrated with 1055 intracardiac electrograms obtained from both atria’s in a persistent AF patient. Notably, the 3D maps did not provide any clinically significant information on rotor pivot point identification or the presence of rotors themselves. Validation of these novel approaches is required in large datasets with paroxysmal and persistent AF patients to evaluate their clinical utility in rotor identification as potential targets for AF ablation

Spatial Characterization and Estimation of Intracardiac Propagation Patterns During Atrial Fibrillation

This doctoral thesis is in the field of biomedical signal processing with focus on methods for the analysis of atrial fibrillation (AF). Paper I of the present thesis addresses the challenge of extracting spatial properties of AF from body surface signals. Different parameters are extracted to estimate the preferred direction of atrial activation and the complexity of the atrial activation pattern. In addition, the relation of the spatial properties to AF organization, which is quantified by AF frequency, is evaluated. While no significant correlation between the preferred direction of atrial activation and AF frequency could be observed, the complexity of the atrial activation pattern was found to increase with AF frequency. The remaining three papers deal with the analysis of the propagation of the electrical activity in the atria during AF based on intracardiac signals. In Paper II, a time-domain method to quantify propagation patterns along a linear catheter based on the detected atrial activation times is developed. Taking aspects on intra-atrial signal organization into account, the detected activation times are combined into wavefronts, and parameters related to the consistency of the wavefronts over time and the activation order along the catheter are extracted. Furthermore, the potential relationship of the extracted parameters to established measures from body surface signals is investigated. While the degree of wavefront consistency was not reflected by the applied body surface measures, AF frequency could distinguish between recordings with different degrees of intra-atrial signal organization. This supports the role of AF frequency as an organization measure of AF. In Paper III, a novel method to analyze intracardiac propagation patterns based on causality analysis in the frequency domain is introduced. In particular, the approach is based on the partial directed coherence (PDC), which evaluates directional coupling between multiple signals in the frequency domain. The potential of the method is illustrated with simulation scenarios based on a detailed ionic model of the human atrial cell as well as with real data recordings, selected to present typical propagation mechanisms and recording situations in atrial tachyarrhythmias. For simulated data, the PDC is correctly reflecting the direction of coupling and thus the propagation between all recording sites. For real data, clear propagation patterns are identified which agree with previous clinical observations. Thus, the results illustrate the ability of the novel approach to identify propagation patterns from intracardiac signals during AF which can provide important information about the underlying AF mechanisms, potentially improving the planning and outcome of ablation. However, spurious couplings over long distances can be observed when analyzing real data comprised by a large number of simultaneously recorded signals, which gives room for further improvement of the method. The derivation of the PDC is entirely based on the fit of a multivariate autoregressive (MVAR) model, commonly estimated by the least-squares (LS) method. In Paper IV, the adaptive group least absolute selection and shrinkage operator (LASSO) is introduced in order to avoid overfitting of the MVAR model and to incorporate prior information such as sparsity of the solution. The sparsity can be motivated by the observation that direct couplings over longer distances are likely to be zero during AF; an information which has been further incorporated by proposing distance-adaptive group LASSO. In simulations, adaptive and distance-adaptive group LASSO are found to be superior to LS estimation in terms of both detection and estimation accuracy. In addition, the results of both simulations and real data analysis indicate that further improvements can be achieved when the distance between the recording sites is known or can be estimated. This further promotes the PDC as a method for analysis of AF propagation patterns, which may contribute to a better understanding of AF mechanisms as well as improved AF treatment