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

Multi-modal and multi-subject modular organization of human brain networks

The human brain is a complex network of anatomically interconnected brain areas. Spontaneous neural activity is constrained by this architecture, giving rise to patterns of statistical dependencies between the activity of remote neural elements. The non-trivial relationship between structural and functional connectivity poses many unsolved challenges about cognition, disease, development, learning and aging. While numerous studies have focused on statistical relationships between edge weights in anatomical and functional networks, less is known about dependencies between their modules and communities. In this work, we investigate and characterize the relationship between anatomical and functional modular organization of the human brain, developing a novel multi-layer framework that expands the classical concept of multi-layer modularity. By simultaneously mapping anatomical and functional networks estimated from different subjects into communities, this approach allows us to carry out a multi-subject and multi-modal analysis of the brain's modular organization. Here, we investigate the relationship between anatomical and functional modules during resting state, finding unique and shared structures. The proposed framework constitutes a methodological advance in the context of multi-layer network analysis and paves the way to further investigate the relationship between structural and functional network organization in clinical cohorts, during cognitively demanding tasks, and in developmental or lifespan studies

Computer- und Videospieler

In der vorliegenden Arbeit wurden mögliche Zusammenhänge zwischen Kompetenz- und Kontrollüberzeugungen, Motiven von Computer- und Videospielern und Gewaltspielpräferenz untersucht. Zur Sammlung der Daten von Computer- und Videospielern wurde eine Onlinebefragung durchgeführt. Das Erhebungsinstrument vereinte selbst konstruierte und aus der Literatur übernommene Items. Erhoben wurden soziodemographische Daten, Charakteristika der Nutzung von Bildschirmspielen, die für das Erkenntnisinteresse relevanten Spielmotive sowie die Kompetenz- und Kontrollüberzeugungen der Befragten. Die Ergebnisse zeigen unter anderem, dass Vielspieler mit stark ausgeprägten externalen Kontrollüberzeugungen eine höhere Motivation aufweisen aggressive Kontrolle im Spiel auszuüben, als Vielspieler mit geringen Ausprägungen in Externalität. Die quantitativen Datenauswertungen belegen außerdem, dass es Spielern mit einem hohen Selbstkonzept eigener Fähigkeiten wichtiger ist Leistungen im Spiel zu realisieren, als Spielern mit durchschnittlich ausgeprägten Kompetenzüberzeugungen. Zur Interpretation dieser Befunde kann der Uses-and-Gratifications Approach herangezogen werden. Computer- und Videospieler nutzen virtuelle Welten demnach, um Kontrolldefizite zu kompensieren und handlungsbezogenen Ideenreichtum auszuleben. Die Studie konnte auch zeigen, dass die Motivation von Spielern, auf aggressive Weise Kontrolle im Spiel auszuüben, mit der Präferenz für Gewaltspiele zusammenhängt. Theoretische Erkenntnisse zum Aggressionsbegriff untermauern die Interpretation, dass Gewalt in Spielen besonders geeignet ist, um Gefühle von Macht und Kontrolle zu erzeugen. Daher greifen Spieler mit höherer aggressiver Kontrollmotivation eher zu violenten Spielen

A geothermal aquifer in the dilation zones on the southern margin of the Dublin Basin

This is the author accepted manuscript. the final version is available from Oxford University Press via the DOI in this recordWe present modelling of the geophysical data from the Newcastle area, west of Dublin, Ireland within the framework of the IRETHERM project. IRETHERM's overarching objective was to facilitate a more thorough strategic understanding of Ireland's geothermal energy potential through integrated modelling of new and existing geophysical, geochemical and geological data. The Newcastle area, one of the target localities, is situated at the southern margin of the Dublin Basin, close to the largest conurbation on the island of Ireland in the City of Dublin and surrounds. As part of IRETHERM, magnetotelluric (MT) soundings were carried out in the highly urbanized Dublin suburb in 2011 and 2012, and a description of MT data acquisition, processing methods, multi-dimensional geoelectrical models and porosity modelling with other geophysical data are presented. The MT time series were heavily noise-contaminated and distorted due to electromagnetic noise from nearby industry and Dublin City tram/railway systems. Time series processing was performed using several modern robust codes to obtain reasonably reliable and interpretable MT impedance and geomagnetic transfer function ‘tipper’ estimates at most of the survey locations. The most ‘quiet’ 3-hour subsets of data during the night time, when the DC ‘LUAS’ tram system was not operating, were used in multi-site and multivariate processing. The final 2-D models underwent examination using a stability technique, and the final two 2-D profiles, with reliability estimations expressed through conductance and resistivity, were derived. In the final stage of this study, 3-D modelling of all magnetotelluric data in the Newcastle area was also undertaken. Comparison of the MT models and their interpretation with existing seismic profiles in the area reveals that the Blackrock to Newcastle Fault (BNF) zone is visible in the models as a conductive feature down to depths of 4 km. The investigated area below Newcastle can be divided into two domains of different depths, formed as depth zones. The first zone, from the surface down to 1–2 km, is dominated by NE-SW oriented conductors connected with shallow faults or folds probably filled with less saline waters. The conductors are also crossing the surface trace of the BNF. The second depth domain can be identified from depths of 2 km to 4 km, where structures are oriented along the BNF and the observed conductivity is lower. The deeper conductive layers are interpreted as geothermal-fluid-bearing rocks. Porosity and permeability estimations from the lithological borehole logs indicate the geothermal potential of the bedrock, to deliver warm water to the surface. The fluid permeability estimation, based on Archie's law for porous structures and synthetic studies of fractured zones, suggests a permeability in the range 100 mD–100 D in the study area, which is prospective for geothermal energy exploitation.Science Foundation Ireland (SFI)Slovak Academy of Sciences (SAS)European Union FP7APVVSlovak Grant Agency VEG

Exploiting cache locality at run-time

With the increasing gap between the speeds of the processor and memory system, memory access has become a major performance bottleneck in modern computer systems. Recently, Symmetric Multi-Processor (SMP) systems have emerged as a major class of high-performance platforms. Improving the memory performance of Parallel applications with dynamic memory-access patterns on Symmetric Multi-Processors (SMP) is a hard problem. The solution to this problem is critical to the successful use of the SMP systems because dynamic memory-access patterns occur in many real-world applications. This dissertation is aimed at solving this problem.;Based on a rigorous analysis of cache-locality optimization, we propose a memory-layout oriented run-time technique to exploit the cache locality of parallel loops. Our technique have been implemented in a run-time system. Using simulation and measurement, we have shown our run-time approach can achieve comparable performance with compiler optimizations for those regular applications, whose load balance and cache locality can be well optimized by tiling and other program transformations. However, our approach was shown to improve significantly the memory performance for applications with dynamic memory-access patterns. Such applications are usually hard to optimize with static compiler optimizations.;Several contributions are made in this dissertation. We present models to characterize the complexity and present a solution framework for optimizing cache locality. We present an effective estimation technique for memory-access patterns to support efficient locality optimizations and information integration. We present a memory-layout oriented run-time technique for locality optimization. We present efficient scheduling algorithms to trade off locality and load imbalance. We provide a detailed performance evaluation of the run-time technique

Domain decomposition methods for coupled Stokes-Darcy flows

This thesis studies the numerical methods for coupled Stokes-Darcy problem. It consists of three major parts: First, a non-overlapping domain decomposition method is presented for Stokes-Darcy problem by partitioning the computational domain into multiple subdomains, upon which families of coupled local problems of lower complexity are formulated. The coupling is based on appropriate interface matching conditions. The global problem is reduced to an interface problem by eliminating the interior subdomain variables, which can be solved by an iterative procedure. FETI approach is used for floating Stokes subdomains. The condition number of the resulting algebraic system is analyzed and numerical tests on matching grids verifying the theoretical estimates are provided. Second, a multiscale flux basis algorithm is developed based on the domain decomposition with multiscale mortar mixed finite element method. The algorithm involves precomputing a multiscale flux basis, which consists of the flux (or velocity trace) response from each mortar degree of freedom. It is computed by each subdomain independently before the interface iteration begins. The subdomain solves required at each iteration are substituted by a linear combination of the multiscale basis. This may lead to a significant reduction in computational cost since the number of subdomain solves is fixed, depending only on the number of mortar degrees of freedom associated with a subdomain. Several numerical examples are carried out to demonstrate the efficiency of the multiscale flux basis implementation. Third, a multiscale flux basis implementation is presented for coupled Stokes-Darcy flows with stochastic permeability, with its log represented as a sum of local Karhunen-Lo\`{e}ve expansions. The problem is approximated by stochastic collocation on either a tensor product or a sparse grid, coupled with multiscale mortar mixed finite element method using non-overlapping domain decomposition for the spatial discretization. Two algorithms based on deterministic or stochastic multiscale flux basis are introduced. Some numerical tests are presented to illustrate the performances of these algorithms, with the stochastic multiscale flux basis showing a great advantage in computational cost among all

A Stochastic Modeling Approach to Region-and Edge-Based Image Segmentation

The purpose of image segmentation is to isolate objects in a scene from the background. This is a very important step in any computer vision system since various tasks, such as shape analysis and object recognition, require accurate image segmentation. Image segmentation can also produce tremendous data reduction. Edge-based and region-based segmentation have been examined and two new algorithms based on recent results in random field theory have been developed. The edge-based segmentation algorithm uses the pixel gray level intensity information to allocate object boundaries in two stages: edge enhancement, followed by edge linking. Edge enhancement is accomplished by maximum energy filters used in one-dimensional bandlimited signal analysis. The issue of optimum filter spatial support is analyzed for ideal edge models. Edge linking is performed by quantitative sequential search using the Stack algorithm. Two probabilistic search metrics are introduced and their optimality is proven and demonstrated on test as well as real scenes. Compared to other methods, this algorithm is shown to produce more accurate allocation of object boundaries. Region-based segmentation was modeled as a MAP estimation problem in which the actual (unknown) objects were estimated from the observed (known) image by a recursive classification algorithms. The observed image was modeled by an Autoregressive (AR) model whose parameters were estimated locally, and a Gibbs-Markov random field (GMRF) model was used to model the unknown scene. A computational study was conducted on images having various types of texture images. The issues of parameter estimation, neighborhood selection, and model orders were examined. It is concluded that the MAP approach for region segmentation generally works well on images having a large content of microtextures which can be properly modeled by both AR and GMRF models. On these texture images, second order AR and GMRF models were shown to be adequate

Document preprocessing and fuzzy unsupervised character classification

This dissertation presents document preprocessing and fuzzy unsupervised character classification for automatically reading daily-received office documents that have complex layout structures, such as multiple columns and mixed-mode contents of texts, graphics and half-tone pictures. First, the block segmentation algorithm is performed based on a simple two-step run-length smoothing to decompose a document into single-mode blocks. Next, the block classification is performed based on the clustering rules to classify each block into one of the types such as text, horizontal or vertical lines, graphics, and pictures. The mean white-to-black transition is shown as an invariance for textual blocks, and is useful for block discrimination. A fuzzy model for unsupervised character classification is designed to improve the robustness, correctness, and speed of the character recognition system. The classification procedures are divided into two stages. The first stage separates the characters into seven typographical categories based on word structures of a text line. The second stage uses pattern matching to classify the characters in each category into a set of fuzzy prototypes based on a nonlinear weighted similarity function. A fuzzy model of unsupervised character classification, which is more natural in the representation of prototypes for character matching, is defined and the weighted fuzzy similarity measure is explored. The characteristics of the fuzzy model are discussed and used in speeding up the classification process. After classification, the character recognition procedure is simply applied on the limited versions of the fuzzy prototypes. To avoid information loss and extra distortion, an topography-based approach is proposed to apply directly on the fuzzy prototypes to extract the skeletons. First, a convolution by a bell-shaped function is performed to obtain a smooth surface. Second, the ridge points are extracted by rule-based topographic analysis of the structure. Third, a membership function is assigned to ridge points with values indicating the degrees of membership with respect to the skeleton of an object. Finally, the significant ridge points are linked to form strokes of skeleton, and the clues of eigenvalue variation are used to deal with degradation and preserve connectivity. Experimental results show that our algorithm can reduce the deformation of junction points and correctly extract the whole skeleton although a character is broken into pieces. For some characters merged together, the breaking candidates can be easily located by searching for the saddle points. A pruning algorithm is then applied on each breaking position. At last, a multiple context confirmation can be applied to increase the reliability of breaking hypotheses

Bayesian jackknife tests with a small number of subsets: application to HERA 21 cm power spectrum upper limits

We present a Bayesian jackknife test for assessing the probability that a data set contains biased subsets, and, if so, which of the subsets are likely to be biased. The test can be used to assess the presence and likely source of statistical tension between different measurements of the same quantities in an automated manner. Under certain broadly applicable assumptions, the test is analytically tractable. We also provide an open-source code, CHIBORG, that performs both analytic and numerical computations of the test on general Gaussian-distributed data. After exploring the information theoretical aspects of the test and its performance with an array of simulations, we apply it to data from the Hydrogen Epoch of Reionization Array (HERA) to assess whether different sub-seasons of observing can justifiably be combined to produce a deeper 21 cm power spectrum upper limit. We find that, with a handful of exceptions, the HERA data in question are statistically consistent and this decision is justified. We conclude by pointing out the wide applicability of this test, including to CMB experiments and the H0 tension