Integrated Cardiac Electromechanics: Modeling and Personalization

Cardiac disease remains the leading cause of morbidity and mortality in the world. A variety of heart diagnosis techniques have been developed during the last century, and generally fall into two groups. The first group evaluates the electrical function of the heart using electrophysiological data such as electrocardiogram (ECG), while the second group aims to assess the mechanical function of the heart through medical imaging data. Nevertheless, the heart is an integrated electromechanical organ, where its cyclic pumping arises from the synergy of its electrical and mechanical function which requires first to be electrically excited in order to contract. At the same time, cardiac electrical function experiences feedback from mechanical contraction. This inter-dependent relationship determines that neither electrical function nor mechanical function alone can completely reflect the pathophysiological conditions of the heart. The aim of this thesis is working towards building an integrated framework for heart diagnosis through evaluation of electrical and mechanical functions simultaneously. The basic rational is to obtain quantitative interpretation of a subject-specific heart system by combining an electromechanical heart model and individual clinical measurements of the heart. To this end, we first develop a biologically-inspired mathematical model of the heart that provides a general, macroscopic description of cardiac electromechanics. The intrinsic electromechanical coupling arises from both excitation-induced contraction and deformation-induced mechano-electrical feedback. Then, as a first step towards a fully electromechanically integrated framework, we develop a model-based approach for investigating the effect of cardiac motion on noninvasive transmural imaging of cardiac electrophysiology. Specifically, we utilize the proposed heart model to obtain updated heart geometry through simulation, and further recover the electrical activities of the heart from body surface potential maps (BSPMs) by solving an optimization problem. Various simulations of the heart have been performed under healthy and abnormal conditions, which demonstrate the physiological plausibility of the proposed integrated electromechanical heart model. What\u27s more, this work presents the effect of cardiac motion to the solution of noninvasive estimation of cardiac electrophysiology and shows the importance of integrating cardiac electrical and mechanical functions for heart diagnosis. This thesis also paves the road for noninvasive evaluation of cardiac electromechanics

Computational modelling of the human heart and multiscale simulation of its electrophysiological activity aimed at the treatment of cardiac arrhythmias related to ischaemia and Infarction

[ES] Las enfermedades cardiovasculares constituyen la principal causa de morbilidad y mortalidad a nivel mundial, causando en torno a 18 millones de muertes cada año. De entre ellas, la más común es la enfermedad isquémica cardíaca, habitualmente denominada como infarto de miocardio (IM). Tras superar un IM, un considerable número de pacientes desarrollan taquicardias ventriculares (TV) potencialmente mortales durante la fase crónica del IM, es decir, semanas, meses o incluso años después la fase aguda inicial. Este tipo concreto de TV normalmente se origina por una reentrada a través de canales de conducción (CC), filamentos de miocardio superviviente que atraviesan la cicatriz del infarto fibrosa y no conductora. Cuando los fármacos anti-arrítmicos resultan incapaces de evitar episodios recurrentes de TV, la ablación por radiofrecuencia (ARF), un procedimiento mínimamente invasivo realizado mediante cateterismo en el laboratorio de electrofisiología (EF), se usa habitualmente para interrumpir de manera permanente la propagación eléctrica a través de los CCs responsables de la TV. Sin embargo, además de ser invasivo, arriesgado y requerir mucho tiempo, en casos de TVs relacionadas con IM crónico, hasta un 50% de los pacientes continúa padeciendo episodios recurrentes de TV tras el procedimiento de ARF. Por tanto, existe la necesidad de desarrollar nuevas estrategias pre-procedimiento para mejorar la planificación de la ARF y, de ese modo, aumentar esta tasa de éxito relativamente baja. En primer lugar, realizamos una revisión exhaustiva de la literatura referente a los modelos cardiacos 3D existentes, con el fin de obtener un profundo conocimiento de sus principales características y los métodos usados en su construcción, con especial atención sobre los modelos orientados a simulación de EF cardíaca. Luego, usando datos clínicos de un paciente con historial de TV relacionada con infarto, diseñamos e implementamos una serie de estrategias y metodologías para (1) generar modelos computacionales 3D específicos de paciente de ventrículos infartados que puedan usarse para realizar simulaciones de EF cardíaca a nivel de órgano, incluyendo la cicatriz del infarto y la región circundante conocida como zona de borde (ZB); (2) construir modelos 3D de torso que permitan la obtención del ECG simulado; y (3) llevar a cabo estudios in-silico de EF personalizados y pre-procedimiento, tratando de replicar los verdaderos estudios de EF realizados en el laboratorio de EF antes de la ablación. La finalidad de estas metodologías es la de localizar los CCs en el modelo ventricular 3D para ayudar a definir los objetivos de ablación óptimos para el procedimiento de ARF. Por último, realizamos el estudio retrospectivo por simulación de un caso, en el que logramos inducir la TV reentrante relacionada con el infarto usando diferentes configuraciones de modelado para la ZB. Validamos nuestros resultados mediante la reproducción, con una precisión razonable, del ECG del paciente en TV, así como en ritmo sinusal a partir de los mapas de activación endocárdica obtenidos invasivamente mediante sistemas de mapeado electroanatómico en este último caso. Esto permitió encontrar la ubicación y analizar las características del CC responsable de la TV clínica. Cabe destacar que dicho estudio in-silico de EF podría haberse efectuado antes del procedimiento de ARF, puesto que nuestro planteamiento está completamente basado en datos clínicos no invasivos adquiridos antes de la intervención real. Estos resultados confirman la viabilidad de la realización de estudios in-silico de EF personalizados y pre-procedimiento de utilidad, así como el potencial del abordaje propuesto para llegar a ser en un futuro una herramienta de apoyo para la planificación de la ARF en casos de TVs reentrantes relacionadas con infarto. No obstante, la metodología propuesta requiere de notables mejoras y validación por medio de es[CA] Les malalties cardiovasculars constitueixen la principal causa de morbiditat i mortalitat a nivell mundial, causant entorn a 18 milions de morts cada any. De elles, la més comuna és la malaltia isquèmica cardíaca, habitualment denominada infart de miocardi (IM). Després de superar un IM, un considerable nombre de pacients desenvolupen taquicàrdies ventriculars (TV) potencialment mortals durant la fase crònica de l'IM, és a dir, setmanes, mesos i fins i tot anys després de la fase aguda inicial. Aquest tipus concret de TV normalment s'origina per una reentrada a través dels canals de conducció (CC), filaments de miocardi supervivent que travessen la cicatriu de l'infart fibrosa i no conductora. Quan els fàrmacs anti-arítmics resulten incapaços d'evitar episodis recurrents de TV, l'ablació per radiofreqüència (ARF), un procediment mínimament invasiu realitzat mitjançant cateterisme en el laboratori de electrofisiologia (EF), s'usa habitualment per a interrompre de manera permanent la propagació elèctrica a través dels CCs responsables de la TV. No obstant això, a més de ser invasiu, arriscat i requerir molt de temps, en casos de TVs relacionades amb IM crònic fins a un 50% dels pacients continua patint episodis recurrents de TV després del procediment d'ARF. Per tant, existeix la necessitat de desenvolupar noves estratègies pre-procediment per a millorar la planificació de l'ARF i, d'aquesta manera, augmentar la taxa d'èxit, que es relativament baixa. En primer lloc, realitzem una revisió exhaustiva de la literatura referent als models cardíacs 3D existents, amb la finalitat d'obtindre un profund coneixement de les seues principals característiques i els mètodes usats en la seua construcció, amb especial atenció sobre els models orientats a simulació de EF cardíaca. Posteriorment, usant dades clíniques d'un pacient amb historial de TV relacionada amb infart, dissenyem i implementem una sèrie d'estratègies i metodologies per a (1) generar models computacionals 3D específics de pacient de ventricles infartats capaços de realitzar simulacions de EF cardíaca a nivell d'òrgan, incloent la cicatriu de l'infart i la regió circumdant coneguda com a zona de vora (ZV); (2) construir models 3D de tors que permeten l'obtenció del ECG simulat; i (3) dur a terme estudis in-silico de EF personalitzats i pre-procediment, tractant de replicar els vertaders estudis de EF realitzats en el laboratori de EF abans de l'ablació. La finalitat d'aquestes metodologies és la de localitzar els CCs en el model ventricular 3D per a ajudar a definir els objectius d'ablació òptims per al procediment d'ARF. Finalment, a manera de prova de concepte, realitzem l'estudi retrospectiu per simulació d'un cas, en el qual aconseguim induir la TV reentrant relacionada amb l'infart usant diferents configuracions de modelatge per a la ZV. Validem els nostres resultats mitjançant la reproducció, amb una precisió raonable, del ECG del pacient en TV, així com en ritme sinusal a partir dels mapes d'activació endocardíac obtinguts invasivament mitjançant sistemes de mapatge electro-anatòmic en aquest últim cas. Això va permetre trobar la ubicació i analitzar les característiques del CC responsable de la TV clínica. Cal destacar que aquest estudi in-silico de EF podria haver-se efectuat abans del procediment d'ARF, ja que el nostre plantejament està completament basat en dades clíniques no invasius adquirits abans de la intervenció real. Aquests resultats confirmen la viabilitat de la realització d'estudis in-silico de EF personalitzats i pre-procediment d'utilitat, així com el potencial de l'abordatge proposat per a arribar a ser en un futur una eina de suport per a la planificació de l'ARF en casos de TVs reentrants relacionades amb infart. No obstant això, la metodologia proposada requereix de notables millores i validació per mitjà d'estudis de simulació amb grans cohorts de pacients.[EN] Cardiovascular diseases represent the main cause of morbidity and mortality worldwide, causing around 18 million deaths every year. Among these diseases, the most common one is the ischaemic heart disease, usually referred to as myocardial infarction (MI). After surviving to a MI, a considerable number of patients develop life-threatening ventricular tachycardias (VT) during the chronic stage of the MI, that is, weeks, months or even years after the initial acute phase. This particular type of VT is typically sustained by reentry through slow conducting channels (CC), which are filaments of surviving myocardium that cross the non-conducting fibrotic infarct scar. When anti-arrhythmic drugs are unable to prevent recurrent VT episodes, radiofrequency ablation (RFA), a minimally invasive procedure performed by catheterization in the electrophysiology (EP) laboratory, is commonly used to interrupt the electrical conduction through the CCs responsible for the VT permanently. However, besides being invasive, risky and time-consuming, in the cases of VTs related to chronic MI, up to 50% of patients continue suffering from recurrent VT episodes after the RFA procedure. Therefore, there exists a need to develop novel pre-procedural strategies to improve RFA planning and, thereby, increase this relatively low success rate. First, we conducted an exhaustive review of the literature associated with the existing 3D cardiac models in order to gain a deep knowledge about their main features and the methods used for their construction, with special focus on those models oriented to simulation of cardiac EP. Later, using a clinical dataset of a chronically infarcted patient with a history of infarct-related VT, we designed and implemented a number of strategies and methodologies to (1) build patient-specific 3D computational models of infarcted ventricles that can be used to perform simulations of cardiac EP at the organ level, including the infarct scar and the surrounding region known as border zone (BZ); (2) construct 3D torso models that enable to compute the simulated ECG; and (3) carry out pre-procedural personalized in-silico EP studies, trying to replicate the actual EP studies conducted in the EP laboratory prior to the ablation. The goal of these methodologies is to allow locating the CCs into the 3D ventricular model in order to help in defining the optimal ablation targets for the RFA procedure. Lastly, as a proof-of-concept, we performed a retrospective simulation case study, in which we were able to induce an infarct-related reentrant VT using different modelling configurations for the BZ. We validated our results by reproducing with a reasonable accuracy the patient's ECG during VT, as well as in sinus rhythm from the endocardial activation maps invasively recorded via electroanatomical mapping systems in this latter case. This allowed us to find the location and analyse the features of the CC responsible for the clinical VT. Importantly, such in-silico EP study might have been conducted prior to the RFA procedure, since our approach is completely based on non-invasive clinical data acquired before the real intervention. These results confirm the feasibility of performing useful pre-procedural personalized in-silico EP studies, as well as the potential of the proposed approach to become a helpful tool for RFA planning in cases of infarct-related reentrant VTs in the future. Nevertheless, the developed methodology requires further improvements and validation by means of simulation studies including large cohorts of patients.During the carrying out of this doctoral thesis, the author Alejandro Daniel López Pérez was financially supported by the Ministerio de Economía, Industria y Competitividad of Spain through the program Ayudas para contratos predoctorales para la formación de doctores, with the grant number BES-2013-064089.López Pérez, AD. (2019). Computational modelling of the human heart and multiscale simulation of its electrophysiological activity aimed at the treatment of cardiac arrhythmias related to ischaemia and Infarction [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/124973TESI

A new algorithm to diagnose atrial ectopic origin from multi lead ECG systems - insights from 3D virtual human atria and torso

Rapid atrial arrhythmias such as atrial fibrillation (AF) predispose to ventricular arrhythmias, sudden cardiac death and stroke. Identifying the origin of atrial ectopic activity from the electrocardiogram (ECG) can help to diagnose the early onset of AF in a cost-effective manner. The complex and rapid atrial electrical activity during AF makes it difficult to obtain detailed information on atrial activation using the standard 12-lead ECG alone. Compared to conventional 12-lead ECG, more detailed ECG lead configurations may provide further information about spatio-temporal dynamics of the body surface potential (BSP) during atrial excitation. We apply a recently developed 3D human atrial model to simulate electrical activity during normal sinus rhythm and ectopic pacing. The atrial model is placed into a newly developed torso model which considers the presence of the lungs, liver and spinal cord. A boundary element method is used to compute the BSP resulting from atrial excitation. Elements of the torso mesh corresponding to the locations of the placement of the electrodes in the standard 12-lead and a more detailed 64-lead ECG configuration were selected. The ectopic focal activity was simulated at various origins across all the different regions of the atria. Simulated BSP maps during normal atrial excitation (i.e. sinoatrial node excitation) were compared to those observed experimentally (obtained from the 64-lead ECG system), showing a strong agreement between the evolution in time of the simulated and experimental data in the P-wave morphology of the ECG and dipole evolution. An algorithm to obtain the location of the stimulus from a 64-lead ECG system was developed. The algorithm presented had a success rate of 93%, meaning that it correctly identified the origin of atrial focus in 75/80 simulations, and involved a general approach relevant to any multi-lead ECG system. This represents a significant improvement over previously developed algorithms

Electrocardiographic Consequences Of Electrical And Anatomical Remodeling In Diabetic And Obese Humans

Background. Diabetes and obesity are two major risk factors for cardiovascular disease. Both can cause changes due to cardiac sources in body-surface potentials: BSPs), that is, in electrocardiograms: ECGs). By identifying the major effects of diabetes and obesity on BSPs, we hope to reveal the electrical phenotype of diabetes in body-surface ECGs in the presence of obesity. Methods. A Bidomain Platform was constructed to link heart-surface transmembrane potentials: TMPs) and BSPs. The Forward-Problem Module of the platform calculates BSPs from a bidomain-model of myocyte TMPs and torso anatomy. The platform also contains a Cardio-myocyte TMP Estimation Module in which an innovative method, named regularized waveform identification: RWI), was developed. It is a new approach used to reconstruct the TMPs from BSPs, that is, solving the electrocardiographic inverse bidomain problem. Using normal TMPs, BSPs were simulated on obese torsos and compared to BSPs on a normal torso to determine ECG changes that might accompany certain obese habitus. BSPs on a normal torso were also simulated with both normal TMPs and TMPs whose duration was increased in a manner expected to occur in the diabetic. In addition, BSPs were measured, heart and torso models were found on two adult male subjects: one normal and one obese diabetic. BSPs and estimated TMPs in these subjects, found by using the RWI method, were compared to identify ECG changes that might be found in the obese diabetic in a clinical setting. Results. Forward-problem solutions found for obese heart-torso models with normal TMPs compared to normal had relative errors: RE) of 12, 30, and 68\% for 20\% left-ventricular hypertrophy, 16\% abdomen extension, and displaced heart, respectively. These results suggest that standard 12-lead ECG measurement could be significantly affected by the anatomical changes associated with obesity. Simulation results also showed diabetic electrical remodeling may have a strong impact on BSPs. An RE of 125\% was observed between normal and diabetic BSPs due to prolongation of recovery that might accompany diabetes. Energy reduction of BSPs was found in both simulated and measured BSPs with obesity. Although QT interval prolongation found in simulated BSPs was not seen in the ECGs recorded from the obese diabetic subject, QT dispersion(QTd) was found increased in diabetic in both simulated and measure ECGs. Obviously, no statistical conclusions can be reached with our limited data set, but the suggestive results call for further clinical observations. TMPs were estimated in realistic, normal heart-torso model simulations using the RWI method. REs of about 15\% were found for up to 10\% noise added to BSPs; and for errors in heart size of 10\% and heart location of $\pm 1$ cm, which were significant improvements over conventional regularization methods alone. Conclusions. In this study, we characterized electrical changes with diabetes and anatomical changes with obesity; then independently evaluated their influences on body surface potentials: BSPs). These results suggest that standard 12-lead ECG measurements could be strongly influenced by the anatomical changes associated with obesity. Body-surface maps and inverse solutions to the heart-surface that minimize volume-conductor effects are likely to be more useful in investigating the influence of diabetic electrical remodeling among obese diabetic patients. Simulation results showed that the RWI inverse solution performed much better than traditional regularization methods alone and is robust in the presence of noise and geometric error. By incorporating temporal information, in the form of the basic TMP wave shape, estimation accuracy was enhanced while maintaining computational simplicity

Personalized noninvasive imaging of volumetric cardiac electrophysiology

Three-dimensionally distributed electrical functioning is the trigger of mechanical contraction of the heart. Disturbance of this electrical flow is known to predispose to mechanical catastrophe but, due to its amenability to certain intervention techniques, a detailed understanding of subject-specific cardiac electrophysiological conditions is of great medical interest. In current clinical practice, body surface potential recording is the standard tool for diagnosing cardiac electrical dysfunctions. However, successful treatments normally require invasive catheter mapping for a more detailed observation of these dysfunctions. In this dissertation, we take a system approach to pursue personalized noninvasive imaging of volumetric cardiac electrophysiology. Under the guidance of existing scientific knowledge of the cardiac electrophysiological system, we extract the subject specific cardiac electrical information from noninvasive body surface potential mapping and tomographic imaging data of individual subjects. In this way, a priori knowledge of system physiology leads the physiologically meaningful interpretation of personal data; at the same time, subject-specific information contained in the data identifies parameters in individual systems that differ from prior knowledge. Based on this perspective, we develop a physiological model-constrained statistical framework for the quantitative reconstruction of the electrical dynamics and inherent electrophysiological property of each individual cardiac system. To accomplish this, we first develop a coupled meshfree-BE (boundary element) modeling approach to represent existing physiological knowledge of the cardiac electrophysiological system on personalized heart-torso structures. Through a state space system approach and sequential data assimilation techniques, we then develop statistical model-data coupling algorithms for quantitative reconstruction of volumetric transmembrane potential dynamics and tissue property of 3D myocardium from body surface potential recoding of individual subjects. We also introduce a data integration component to build personalized cardiac electrophysiology by fusing tomographic image and BSP sequence of the same subject. In addition, we develop a computational reduction strategy that improves the efficiency and stability of the framework. Phantom experiments and real-data human studies are performed for validating each of the framework’s major components. These experiments demonstrate the potential of our framework in providing quantitative understanding of volumetric cardiac electrophysiology for individual subjects and in identifying latent threats in individual’s heart. This may aid in personalized diagnose, treatment planning, and fundamentally, prevention of fatal cardiac arrhythmia

Estimation of Atrial Electrical Complexity during Atrial Fibrillation by Solving the Inverse Problem of Electrocardiography

Tesis por compendio[ES] La fibrilación auricular (FA) es la arritmia más prevalente en el mundo y está asociada con una elevada morbilidad, mortalidad y costes sanitarios. A pesar de los avances en opciones de tratamiento farmacológico y terapia de ablación, el manejo de la FA todavía tiene margen de mejora. La imagen electrocardiográfica (ECGI) se ha destacado como un prometedor método no invasivo para evaluar la electrofisiología cardíaca y guiar las decisiones terapéuticas en casos de fibrilación auricular. No obstante, el ECGI se enfrenta a desafíos como la necesidad de resolver de manera precisa el denominado problema inverso de la electrocardiografía y de optimizar la calidad de las reconstrucciones de ECGI. Además, la integración del ECGI en los procesos clínicos rutinarios sigue siendo un reto, en gran medida debido a los costos que supone la necesidad de imágenes cardíacas. Por ello, los objetivos principales de esta tesis doctoral son impulsar la tecnología ECGI mediante la determinación de sus requisitos técnicos mínimos y la mejora de las metodologías existentes para obtener señales de ECGI precisas. Asimismo, buscamos evaluar la capacidad de ECGI para cuantificar de forma no invasiva la complejidad de la FA. Para lograr estos objetivos, se han llevado a cabo diversos estudios a lo largo de la tesis, desde el perfeccionamiento del ECGI hasta la evaluación de la FA utilizando esta tecnología. En primer lugar, se han estudiado los requisitos geométricos y de señal del problema inverso mediante el estudio de los efectos de la densidad de la malla del torso y la distribución de electrodos en la precisión del ECGI, lo que ha conducido a la identificación del número mínimo de nodos y su distribución en la malla del torso. Además, hemos identificado que para obtener señales de ECGI de alta calidad, es crucial la correcta disposición de los electrodos en la malla del torso reconstruido. Asimismo, se ha definido y evaluado una nueva metodología de ECGI sin necesidad de usar técnicas de imagen cardiaca. Para ello, hemos comparado métricas derivadas del ECGI calculadas con la geometría original del corazón de los pacientes con las métricas medidas en diferentes geometrías cardíacas. Nuestros resultados han mostrado que el ECGI sin necesidad de imágenes cardíacas es efectivo para la correcta cuantificación y localización de los patrones y zonas que mantienen la FA. En paralelo, hemos optimizado la regularización de Tikhonov de orden cero actual y la optimización de la curva L para el cálculo de las señales ECGI, investigando cómo el ruido eléctrico y las incertidumbres geométricas influyen en la regularización. A partir de ello, propusimos un nuevo criterio que realza la precisión de las soluciones de ECGI en escenarios con incertidumbre debido a condiciones de señal no ideales. En segundo lugar, en esta tesis doctoral, se han llevado a cabo múltiples análisis relativos a diferentes metodologías de procesado de señales y obtención métricas derivadas del ECGI con el fin de caracterizar mejor el sustrato cardíaco y la actividad reentrante en las señales de ECGI de pacientes con FA. Con el objetivo de obtener una comprensión más profunda de los mecanismos electrofisiológicos subyacentes a la FA, hemos establecido la estrategia de filtrado óptima para extraer patrones reentrantes específicos del paciente y métricas derivadas de señales ECGI. Además, hemos investigado la reproducibilidad de los mapas de reentradas derivados de las señales de ECGI y hemos encontrado su relación con el éxito de la ablación de venas pulmonares (PVI). Nuestros resultados han mostrado que una mayor reproducibilidad en los patrones reentrantes de FA detectados con ECGI está relacionada con el éxito de la PVI, creando una metodología para estratificar a los pacientes con FA antes de los procedimientos de ablación.[CA] La fibril·lació auricular (FA) és l'arrítmia més prevalent al món i està associada amb una elevada morbiditat, mortalitat i costos sanitaris. Malgrat els avanços en opcions de tractament farmacològic i teràpies d'ablació, el maneig de la FA encara té marge de millora. La imatge electrocardiogràfica (ECGI) s'ha destacat com un prometedor mètode no invasiu per a avaluar l'electrofisiologia cardíaca i guiar les decisions terapèutiques en casos de fibril·lació auricular. No obstant això, l'ECGI s'enfronta a desafiaments com la necessitat de resoldre de manera precisa el denominat problema invers de la electrocardiografia i d'optimitzar la qualitat de les reconstruccions de ECGI. A més, la integració del ECGI en els processos clínics rutinaris continua sent un repte, en gran manera a causa dels costos que suposa la necessitat d'imatges cardíaques. Per això, els objectius principals d'aquesta tesi doctoral són impulsar la tecnologia de l'ECGI mitjançant la determinació dels seus requisits tècnics mínims i la millora de les metodologies existents per obtenir senyals d'ECGI precises. A més, busquem avaluar la capacitat de l'ECGI per quantificar de forma no invasiva la complexitat de la FA. Per a aconseguir aquests objectius, s'han dut a terme diversos estudis al llarg de la tesi, des del perfeccionament de l'ECGI fins a l'avaluació de la FA utilitzant aquesta tecnologia. En primer lloc, hem estudiat els requisits geomètrics i de senyal del problema invers mitjançant l'estudi dels efectes de la densitat de la malla del tors i la distribució d'elèctrodes en la precisió de l'ECGI, el que ha conduït a la identificació del nombre mínim de nodes i la seva distribució en la malla del tors. A més, hem identificat que per obtindre senyals d'ECGI d'alta qualitat, és crucial la correcta disposició dels elèctrodes en la malla del tors reconstruïda. També s'ha definit i avaluat una nova metodologia d'ECGI sense necessitat d'utilitzar tècniques d'imatge cardíaca. Per a això, hem comparat mètriques derivades de l'ECGI calculades amb la geometria original del cor dels pacients amb les mètriques mesurades en diferents geometries cardíaques. Els nostres resultats han mostrat que l'ECGI sense necessitat d'imatges cardíaques és efectiu per a la correcta quantificació i localització dels patrons i zones que mantenen la FA. Paral·lelament, hem optimitzat la regularització de Tikhonov d'ordre zero actual i l'optimització de la corba L per al càlcul de les senyals d'ECGI, investigant com el soroll elèctric i les incerteses geomètriques influeixen en la regularització. Addicionalment, vam proposar un nou criteri que reforça la precisió de les solucions d'ECGI en escenaris amb incertesa degut a condicions de senyal no ideals. En segon lloc, en aquesta tesi doctoral, s'han dut a terme múltiples anàlisis relatius a diferents metodologies de processament de senyals i obtenció de mètriques derivades de l'ECGI amb l'objectiu de caracteritzar millor el substrat cardíac i l'activitat reentrant en les senyals d'ECGI de pacients amb FA. Amb l'objectiu d'obtindre una comprensió més profunda dels mecanismes electrofisiològics subjacents a la FA, hem establert l'estratègia de filtrat òptima per extreure patrons reentrants específics del pacient i mètriques derivades de senyals ECGI. A més, hem investigat la reproductibilitat dels mapes de reentrades derivats de les senyals d'ECGI i hem trobat la seva relació amb l'èxit de l'ablació de venes pulmonars (PVI). Els nostres resultats han mostrat que una major reproductibilitat en els patrons reentrants de FA detectats amb ECGI està relacionada amb l'èxit de la PVI, creant una metodologia per estratificar els pacients amb FA abans dels procediments d'ablació.[EN] Atrial fibrillation (AF) is the most prevalent arrhythmia in the world and is associated with significant morbidity, mortality, and healthcare costs. Despite advancements in pharmaceutical treatment alternatives and ablation therapy, AF management remains suboptimal. Electrocardiographic Imaging (ECGI) has emerged as a promising non-invasive method for assessing cardiac electrophysiology and guiding therapeutic decisions in atrial fibrillation. However, ECGI faces challenges in dealing with accurately resolving the ill-posed inverse problem of electrocardiography and optimizing the quality of ECGI reconstructions. Additionally, the integration of ECGI into clinical workflows is still a challenge that is hindered by the associated costs arising from the need for cardiac imaging. For this purpose, the main objectives of this PhD thesis are to advance ECGI technology by determining the minimal technical requirements and refining existing methodologies for acquiring accurate ECGI signals. In addition, we aim to assess the capacity of ECGI for noninvasively quantifying AF complexity. To fulfill these objectives, several studies were developed throughout the thesis, advancing from ECGI enhancement to AF evaluation using ECGI. Firstly, geometric and signal requirements of the inverse problem were addressed by studying the effects of torso mesh density and electrode distribution on ECGI accuracy, leading to the identification of the minimal number of nodes and their distribution on the torso mesh. Besides, we identified that the correct location of the electrodes on the reconstructed torso mesh is critical for the accurate ECGI signal obtention. Additionally, a new methodology of imageless ECGI was defined and assessed by comparing ECGI-derived drivers computed with the original heart geometry of the patients to the drivers measured in different heart geometries. Our results showed the ability of imageless ECGI to the correct quantification and location of atrial fibrillation drivers, validating the use of ECGI without the need for cardiac imaging. Also, the current state of-the-art zero-order Tikhonov regularization and L-curve optimization for computing ECGI signals were improved by investigating the impact of electrical noise and geometrical uncertainties on the regularization. We proposed a new criterion that enhances the accuracy and reliability of ECGI solutions in situations with uncertainty from unfavorable signal conditions. Secondly, in this PhD thesis, several analyses, signal processing methodologies, and ECGIderived metrics were investigated to better characterize the cardiac substrate and reentrant activity in ECGI signals from AF patients. With the objective of obtaining a deeper understanding of the electrophysiological mechanisms underlying AF, we established the optimal filtering strategy to extract patient-specific reentrant patterns and derived metrics in ECGI signals. Furthermore, we investigated the reproducibility of the obtained ECGI-reentrant maps and linked them to the success of PVI ablation. Our results showed that higher reproducibility on AF drivers detected with ECGI is linked with the success of PVI, creating a proof-of-concept mechanism for stratifying AF patients prior to ablation procedures.This work was supported by: Instituto de Salud Carlos III, and Ministerio de Ciencia e Innovación (supported by FEDER Fondo Europeo de Desarrollo Regional DIDIMO PLEC2021- 007614, ESSENCE PID2020-119364RB-I00, and RYC2018- 024346B-750), EIT Health (Activity code SAVE-COR 220385, EIT Health is supported by EIT, a body of the European Union) and Generalitat Valenciana Conselleria d’Educació, Investigació, Cultura i Esport (ACIF/2020/265 and BEFPI/2021/062).Molero Alabau, R. (2023). Estimation of Atrial Electrical Complexity during Atrial Fibrillation by Solving the Inverse Problem of Electrocardiography [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/199029Compendi

Identification of atrial fibrillation drivers by means of concentric ring electrodes

The prevalence of atrial fibrillation (AF) has tripled in the last 50 years due to population aging. High-frequency (DFdriver) activated atrial regions lead the activation of the rest of the atria, disrupting the propagation wavefront. Fourier based spectral analysis of body surface potential maps have been proposed for DFdriver identification, although these approaches present serious drawbacks due to their limited spectral resolution for short AF epochs and the blurring effect of the volume conductor. Laplacian signals (BC-ECG) from bipolar concentric ring electrodes (CRE) have been shown to outperform the spatial resolution achieved with conventional unipolar recordings. Our aimed was to determine the best DFdriver estimator in endocardial electrograms and to assess the BC-ECG capacity of CRE to quantify AF activity non-invasively

Bayesian Inference with Combined Dynamic and Sparsity Models: Application in 3D Electrophysiological Imaging

Data-driven inference is widely encountered in various scientific domains to convert the observed measurements into information that cannot be directly observed about a system. Despite the quickly-developing sensor and imaging technologies, in many domains, data collection remains an expensive endeavor due to financial and physical constraints. To overcome the limits in data and to reduce the demand on expensive data collection, it is important to incorporate prior information in order to place the data-driven inference in a domain-relevant context and to improve its accuracy. Two sources of assumptions have been used successfully in many inverse problem applications. One is the temporal dynamics of the system (dynamic structure). The other is the low-dimensional structure of a system (sparsity structure). In existing work, these two structures have often been explored separately, while in most high-dimensional dynamic system they are commonly co-existing and contain complementary information. In this work, our main focus is to build a robustness inference framework to combine dynamic and sparsity constraints. The driving application in this work is a biomedical inverse problem of electrophysiological (EP) imaging, which noninvasively and quantitatively reconstruct transmural action potentials from body-surface voltage data with the goal to improve cardiac disease prevention, diagnosis, and treatment. The general framework can be extended to a variety of applications that deal with the inference of high-dimensional dynamic systems

Influence of Myocardial Infarction on QRS Properties: A Simulation Study

The interplay between structural and electrical changes in the heart after myocardial infarction (MI) plays a key role in the initiation and maintenance of arrhythmia. The anatomical and electrophysiological properties of scar, border zone, and normal myocardium modify the electrocardiographic morphology, which is routinely analysed in clinical settings. However, the influence of various MI properties on the QRS is not intuitively predictable.In this work, we have systematically investigated the effects of 17 post-MI scenarios, varying the location, size, transmural extent, and conductive level of scarring and border zone area, on the forward-calculated QRS. Additionally, we have compared the contributions of different QRS score criteria for quantifying post-MI pathophysiology.The propagation of electrical activity in the ventricles is simulated via a Eikonal model on a unified coordinate system.The analysis has been performed on 49 subjects, and the results imply that the QRS is capable of identifying MI, suggesting the feasibility of inversely reconstructing infarct regions from QRS.There exist sensitivity variations of different QRS criteria for identifying 17 MI scenarios, which is informative for solving the inverse problem.Comment: 10 pages, accpeted by FIMH 202