Identification of weakly coupled multiphysics problems. Application to the inverse problem of electrocardiography

This work addresses the inverse problem of electrocardiography from a new perspective, by combining electrical and mechanical measurements. Our strategy relies on the defini-tion of a model of the electromechanical contraction which is registered on ECG data but also on measured mechanical displacements of the heart tissue typically extracted from medical images. In this respect, we establish in this work the convergence of a sequential estimator which combines for such coupled problems various state of the art sequential data assimilation methods in a unified consistent and efficient framework. Indeed we ag-gregate a Luenberger observer for the mechanical state and a Reduced Order Unscented Kalman Filter applied on the parameters to be identified and a POD projection of the electrical state. Then using synthetic data we show the benefits of our approach for the estimation of the electrical state of the ventricles along the heart beat compared with more classical strategies which only consider an electrophysiological model with ECG measurements. Our numerical results actually show that the mechanical measurements improve the identifiability of the electrical problem allowing to reconstruct the electrical state of the coupled system more precisely. Therefore, this work is intended to be a first proof of concept, with theoretical justifications and numerical investigations, of the ad-vantage of using available multi-modal observations for the estimation and identification of an electromechanical model of the heart

Debatable issues in automated ECG reporting

Although automated ECG analysis has been available for many years, there are some aspects which require to be re-assessed with respect to their value while newer techniques which are worthy of review are beginning to find their way into routine use. At the annual International Society of Computerized Electrocardiology conference held in April 2017, four areas in particular were debated. These were a) automated 12 lead resting ECG analysis; b) real time out of hospital ECG monitoring; c) ECG imaging; and d) single channel ECG rhythm interpretation. One speaker presented the positive aspects of each technique and another outlined the more negative aspects. Debate ensued. There were many positives set out for each technique but equally, more negative features were not in short supply, particularly for out of hospital ECG monitoring

Non-invasive identification of atrial fibrillation drivers

Atrial fibrillation (AF) is one of the most common cardiac arrhythmias. Nowadays the fibrillatory process is known to be provoked by the high-frequency reentrant activity of certain atrial regions that propagates the fibrillatory activity to the rest of the atrial tissue, and the electrical isolation of these key regions has demonstrated its effectiveness in terminating the fibrillatory process. The location of the dominant regions represents a major challenge in the diagnosis and treatment of this arrhythmia. With the aim to detect and locate the fibrillatory sources prior to surgical procedure, non-invasive methods have been developed such as body surface electrical mapping (BSPM) which allows to record with high spatial resolution the electrical activity on the torso surface or the electrocardiographic imaging (ECGI) which allows to non-invasively reconstruct the electrical activity in the atrial surface. Given the novelty of these systems, both technologies suffer from a lack of scientific knowledge about the physical and technical mechanisms that support their operation. Therefore, the aim of this thesis is to increase that knowledge, as well as studying the effectiveness of these technologies for the localization of dominant regions in patients with AF. First, it has been shown that BSPM systems are able to noninvasively identify atrial rotors by recognizing surface rotors after band-pass filtering. Furthermore, the position of such surface rotors is related to the atrial rotor location, allowing the distinction between left or right atrial rotors. Moreover, it has been found that the surface electrical maps in AF suffer a spatial smoothing effect by the torso conductor volume, so the surface electrical activity can be studied with a relatively small number of electrodes. Specifically, it has been seen that 12 uniformly distributed electrodes are sufficient for the correct identification of atrial dominant frequencies, while at least 32 leads are needed for non-invasive identification of atrial rotors. Secondly, the effect of narrowband filtering on the effectiveness of the location of reentrant patterns was studied. It has been found that this procedure allows isolating the reentrant electrical activity caused by the rotor, increasing the detection rate for both invasive and surface maps. However, the spatial smoothing caused by the regularization of the ECGI added to the temporal filtering causes a large increase in the spurious reentrant activity, making it difficult to detect real reentrant patterns. However, it has been found that maps provided by the ECGI without temporal filtering allow the correct detection of reentrant activity, so narrowband filtering should be applied for intracavitary or surface signal only. Finally, we studied the stability of the markers used to detect dominant regions in ECGI, such as frequency maps or the rotor presence. It has been found that in the presence of alterations in the conditions of the inverse problem, such as electrical or geometrical noise, these markers are significantly more stable than the ECGI signal morphology from which they are extracted. In addition, a new methodology for error reduction in the atrial spatial location based on the curvature of the curve L has been proposed. The results presented in this thesis showed that BSPM and ECGI systems allows to non-invasively locate the presence of high-frequency rotors, responsible for the maintenance of AF. This detection has been proven to be unambiguous and robust, and the physical and technical mechanisms that support this behavior have been studied. These results indicate that both non-invasive systems provide information of great clinical value in the treatment of AF, so their use can be helpful for selecting and planning atrial ablation procedures.La fibrilación auricular (FA) es una de las arritmias cardiacas más frecuentes. Hoy en día se sabe que el proceso fibrilatorio está provocado por la actividad reentrante a alta frecuencia de ciertas regiones auriculares que propagan la actividad fibrilatoria en el resto del tejido auricular, y se ha demostrado que el aislamiento eléctrico de estas regiones dominantes permite detener el proceso fibrilatorio. La localización de las regiones dominantes supone un gran reto en el diagnóstico y tratamiento de la FA. Con el objetivo de poder localizar las fuentes fibrilatorias con anterioridad al procedimiento quirúrgico, se han desarrollado métodos no invasivos como la cartografía eléctrica de superficie (CES) que registra con gran resolución espacial la actividad eléctrica en la superficie del torso o la electrocardiografía por imagen (ECGI) que permite reconstruir la actividad eléctrica en la superficie auricular. Dada la novedad de estos sistemas, existe una falta de conocimiento científico sobre los mecanismos físicos y técnicos que sustentan su funcionamiento. Por lo tanto, el objetivo de esta tesis es aumentar dicho conocimiento, así como estudiar la eficacia de ambas tecnologías para la localización de regiones dominantes en pacientes con FA. En primer lugar, ha visto que los sistemas CES permiten identificar rotores auriculares mediante el reconocimiento de rotores superficiales tras el filtrado en banda estrecha. Además, la posición de los rotores superficiales está relacionada con la localización de dichos rotores, permitiendo la distinción entre rotores de aurícula derecha o izquierda. Por otra parte, se ha visto que los mapas eléctricos superficiales durante FA sufren una gran suavizado espacial por el efecto del volumen conductor del torso, lo que permite que la actividad eléctrica superficial pueda ser estudiada con un número relativamente reducido de electrodos. Concretamente, se ha visto que 12 electrodos uniformemente distribuidos son suficientes para una correcta identificación de frecuencias dominantes, mientras que son necesarios al menos 32 para una correcta identificación de rotores auriculares. Por otra parte, también se ha estudiado el efecto del filtrado en banda estrecha sobre la eficacia de la localización de patrones reentrantes. Así, se ha visto que este procedimiento permite aislar la actividad eléctrica reentrante provocada por el rotor, aumentando la tasa de detección tanto para señal obtenida de manera invasiva como para los mapas superficiales. No obstante, este filtrado temporal sobre la señal de ECGI provoca un gran aumento de la actividad reentrante espúrea que dificulta la detección de patrones reentrantes reales. Sin embargo, los mapas ECGI sin filtrado temporal permiten la detección correcta de la actividad reentrante, por lo el filtrado debería ser aplicado únicamente para señal intracavitaria o superficial. Por último, se ha estudiado la estabilidad de los marcadores utilizados en ECGI para detectar regiones dominantes, como son los mapas de frecuencia o la presencia de rotores. Se ha visto que en presencia de alteraciones en las condiciones del problema inverso, como ruido eléctrico o geométrico, estos marcadores son significativamente más estables que la morfología de la propia señal ECGI. Además, se ha propuesto una nueva metodología para la reducción del error en la localización espacial de la aurícula basado en la curvatura de la curva L. Los resultados presentados en esta tesis revelan que los sistemas de CES y ECGI permiten localizar de manera no invasiva la presencia de rotores de alta frecuencia. Esta detección es univoca y robusta, y se han estudiado los mecanismos físicos y técnicos que sustentan dicho comportamiento. Estos resultados indican que ambos sistemas no invasivos proporcionan información de gran valor clínico en el tratamiento de la FA, por lo que su uso puede ser de gran ayuda para la selección y planificaciLa fibril·lació auricular (FA) és una de les arítmies cardíaques més freqüents. Hui en dia es sabut que el procés fibrilatori està provocat per l'activitat reentrant de certes regions auriculars que propaguen l'activitat fibril·latoria a la resta del teixit auricular, i s'ha demostrat que l'aïllament elèctric d'aquestes regions dominants permet aturar el procés fibrilatori. La localització de les regions dominants suposa un gran repte en el diagnòstic i tractament d'aquesta arítmia. Amb l'objectiu de poder localitzar fonts fibril·latories amb anterioritat al procediment quirúrgic s'han desenvolupat mètodes no invasius com la cartografia elèctrica de superfície (CES) que registra amb gran resolució espacial l'activitat elèctrica en la superfície del tors o l'electrocardiografia per imatge (ECGI) que permet obtenir de manera no invasiva l'activitat elèctrica en la superfície auricular. Donada la relativa novetat d'aquests sistemes, existeix una manca de coneixement científic sobre els mecanismes físics i tècnics que sustenten el seu funcionament. Per tant, l'objectiu d'aquesta tesi és augmentar aquest coneixement, així com estudiar l'eficàcia d'aquestes tecnologies per a la localització de regions dominants en pacients amb FA. En primer lloc, s'ha vist que els sistemes CES permeten identificar rotors auriculars mitjançant el reconeixement de rotors superficials després del filtrat en banda estreta. A més, la posició dels rotors superficials està relacionada amb la localització d'aquests rotors, permetent la distinció entre rotors de aurícula dreta o esquerra. També s'ha vist que els mapes elèctrics superficials durant FA pateixen un gran suavitzat espacial per l'efecte del volum conductor del tors, el que permet que l'activitat elèctrica superficial pugui ser estudiada amb un nombre relativament reduït d'elèctrodes. Concretament, s'ha vist que 12 elèctrodes uniformement distribuïts són suficients per a una correcta identificació de freqüències dominants auriculars, mentre que són necessaris almenys 32 per a una correcta identificació de rotors auriculars. D'altra banda, també s'ha estudiat l'efecte del filtrat en banda estreta sobre l'eficàcia de la localització de patrons reentrants. Així, s'ha vist que aquest procediment permet aïllar l'activitat elèctrica reentrant provocada pel rotor, augmentant la taxa de detecció tant pel senyal obtingut de manera invasiva com per als mapes superficials. No obstant això, aquest filtrat temporal sobre el senyal de ECGI provoca un gran augment de l'activitat reentrant espúria que dificulta la detecció de patrons reentrants reals. A més, els mapes proporcionats per la ECGI sense filtrat temporal permeten la detecció correcta de l'activitat reentrant, per la qual cosa el filtrat hauria de ser aplicat únicament per a senyal intracavitària o superficial. Per últim, s'ha estudiat l'estabilitat dels marcadors utilitzats en ECGI per a detectar regions auriculars dominants, com són els mapes de freqüència o la presència de rotors. S'ha vist que en presència d'alteracions en les condicions del problema invers, com soroll elèctric o geomètric, aquests marcadors són significativament més estables que la morfologia del mateix senyal ECGI. A més, s'ha proposat una nova metodologia per a la reducció de l'error en la localització espacial de l'aurícula basat en la curvatura de la corba L. Els resultats presentats en aquesta tesi revelen que els sistemes de CES i ECGI permeten localitzar de manera no invasiva la presència de rotors d'alta freqüència. Aquesta detecció és unívoca i robusta, i s'han estudiat els mecanismes físics i tècnics que sustenten aquest comportament. Aquests resultats indiquen que els dos sistemes no invasius proporcionen informació de gran valor clínic en el tractament de la FA, pel que el seu ús pot ser de gran ajuda per a la selecció i planificació de procediments d'ablació auricular.Rodrigo Bort, M. (2016). Non-invasive identification of atrial fibrillation drivers [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/75346TESISPremios Extraordinarios de tesis doctorale

Physiology-based regularization of the electrocardiographic inverse problem

The inverse problem of electrocardiography aims at noninvasively reconstructing electrical activity of the heart from recorded body-surface electrocardiograms. A crucial step is regularization, which deals with ill-posedness of the problem by imposing constraints on the possible solutions. We developed a regularization method that includes electrophysiological input. Body-surface potentials are recorded and a computed tomography scan is performed to obtain the torso-heart geometry. Propagating waveforms originating from several positions at the heart are simulated and used to generate a set of basis vectors representing spatial distributions of potentials on the heart surface. The real heart-surface potentials are then reconstructed from the recorded body-surface potentials by finding a sparse representation in terms of this basis. This method, which we named 'physiology-based regularization' (PBR), was compared to traditional Tikhonov regularization and validated using in vivo recordings in dogs. PBR recovered details of heart-surface electrograms that were lost with traditional regularization, attained higher correlation coefficients and led to improved estimation of recovery times. The best results were obtained by including approximate knowledge about the beat origin in the PBR basis

Numerical simulation of electrocardiograms for full cardiac cycles in healthy and pathological conditions

This work is dedicated to the simulation of full cycles of the electrical activity of the heart and the corresponding body surface potential. The model is based on a realistic torso and heart anatomy, including ventricles and atria. One of the specificities of our approach is to model the atria as a surface, which is the kind of data typically provided by medical imaging for thin volumes. The bidomain equations are considered in their usual formulation in the ventricles, and in a surface formulation on the atria. Two ionic models are used: the Courtemanche-Ramirez-Nattel model on the atria, and the "Minimal model for human Ventricular action potentials" (MV) by Bueno-Orovio, Cherry and Fenton in the ventricles. The heart is weakly coupled to the torso by a Robin boundary condition based on a resistor- capacitor transmission condition. Various ECGs are simulated in healthy and pathological conditions (left and right bundle branch blocks, Bachmann's bundle block, Wolff-Parkinson-White syndrome). To assess the numerical ECGs, we use several qualitative and quantitative criteria found in the medical literature. Our simulator can also be used to generate the signals measured by a vest of electrodes. This capability is illustrated at the end of the article

Reconstruction of atrial ectopic focal and Re-entrant excitations from body surface potentials. Insights from 3D virtual human atria and torso

Non-invasive electrocardiographic imaging has been seen as a painless and economic method to map the electrical functions of the heart. However, it is still a great challenge to obtain accurate reconstruction of cardiac electrical activity from body surface potentials (BSP) due to the ill-posed behaviour of the cardiac inverse-problem. Though some advances have been made in solving the inverse-problem, few studies have been conducted for the atria, which have dramatic differences to the ventricles in their anatomical structures and electrophysiological properties. It is unclear either how the spatial resolution of electrodes on the BSP and rapid excitation rates of atrial activation during atrial fibrillation affect the accuracy of the inverse-problem. In this study, we used a biophysically detailed model of the human atria and torso to investigate effects of multi-lead ECG on the accuracy of reconstructed atrial excitation pattern on the epicardiac surface during the time courses of atrial fibrillation induced by electrical remodelling. It was shown that the solution of the atrial inverse-problem was dependent on the spatial resolution of electrodes on the body surface. The solution was also influenced by the morphology of the AP, rate and types of atrial excitation as well as the implantation of variant orders of the Tikhonov regularization method

Doctor of Philosophy

dissertationInverse Electrocardiography (ECG) aims to noninvasively estimate the electrophysiological activity of the heart from the voltages measured at the body surface, with promising clinical applications in diagnosis and therapy. The main challenge of this emerging technique lies in its mathematical foundation: an inverse source problem governed by partial differential equations (PDEs) which is severely ill-conditioned. Essential to the success of inverse ECG are computational methods that reliably achieve accurate inverse solutions while harnessing the ever-growing complexity and realism of the bioelectric simulation. This dissertation focuses on the formulation, optimization, and solution of the inverse ECG problem based on finite element methods, consisting of two research thrusts. The first thrust explores the optimal finite element discretization specifically oriented towards the inverse ECG problem. In contrast, most existing discretization strategies are designed for forward problems and may become inappropriate for the corresponding inverse problems. Based on a Fourier analysis of how discretization relates to ill-conditioning, this work proposes refinement strategies that optimize approximation accuracy o f the inverse ECG problem while mitigating its ill-conditioning. To fulfill these strategies, two refinement techniques are developed: one uses hybrid-shaped finite elements whereas the other adapts high-order finite elements. The second research thrust involves a new methodology for inverse ECG solutions called PDE-constrained optimization, an optimization framework that flexibly allows convex objectives and various physically-based constraints. This work features three contributions: (1) fulfilling optimization in the continuous space, (2) formulating rigorous finite element solutions, and (3) fulfilling subsequent numerical optimization by a primal-dual interiorpoint method tailored to the given optimization problem's specific algebraic structure. The efficacy o f this new method is shown by its application to localization o f cardiac ischemic disease, in which the method, under realistic settings, achieves promising solutions to a previously intractable inverse ECG problem involving the bidomain heart model. In summary, this dissertation advances the computational research of inverse ECG, making it evolve toward an image-based, patient-specific modality for biomedical research