Gaussian process manifold interpolation for probabilistic atrial activation maps and uncertain conduction velocity.
Ontology highlight
ABSTRACT:
SUBMITTER: Coveney S
PROVIDER: S-EPMC7287339 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Intensity-based image registration requires resampling images on a common grid to evaluate the similarity function. The uncertainty of interpolation varies across the image, depending on the location of resampled points relative to the base grid. We propose to perform Bayesian inference with Gaussian processes, where the covariance matrix of the Gaussian process posterior distribution estimates the uncertainty in interpolation. The Gaussian process replaces a single image with a distribution over images that we integrate into a generative model for registration. Marginalization over resampled images leads to a new similarity measure that includes the uncertainty of the interpolation. We demonstrate that our approach increases the registration accuracy and propose an efficient approximation scheme that enables seamless integration with existing registration methods.
| S-EPMC4219919 | biostudies-literature
Project description:We address the problem of interpolating randomly non-uniformly spatiotemporally scattered uncertain motion measurements, which arises in the context of soft tissue motion estimation. Soft tissue motion estimation is of great interest in the field of image-guided soft-tissue intervention and surgery navigation, because it enables the registration of pre-interventional/pre-operative navigation information on deformable soft-tissue organs. To formally define the measurements as spatiotemporally scattered motion signal samples, we propose a novel motion field representation. To perform the interpolation of the motion measurements in an uncertainty-aware optimal unbiased fashion, we devise a novel Gaussian process (GP) regression model with a non-constant-mean prior and an anisotropic covariance function and show through an extensive evaluation that it outperforms the state-of-the-art GP models that have been deployed previously for similar tasks. The employment of GP regression enables the quantification of uncertainty in the interpolation result, which would allow the amount of uncertainty present in the registered navigation information governing the decisions of the surgeon or intervention specialist to be conveyed.
| S-EPMC6433137 | biostudies-literature
Project description:AimsAccurate determination of intra-atrial conduction velocity (CV) is essential to identify arrhythmogenic areas. The most optimal, commonly used, estimation methodology to measure conduction heterogeneity, including finite differences (FiD), polynomial surface fitting (PSF), and a novel technique using discrete velocity vectors (DVV), has not been determined. We aim (i) to identify the most suitable methodology to unravel local areas of conduction heterogeneities using high-density CV estimation techniques, (ii) to quantify intra-atrial differences in CV, and (iii) to localize areas of CV slowing associated with paroxysmal atrial fibrillation (PAF).Methods and resultsIntra-operative epicardial mapping (>5000 sites, interelectrode distances 2 mm) of the right and left atrium and Bachmann's bundle (BB) was performed during sinus rhythm (SR) in 412 patients with or without PAF. The median atrial CV estimated using the DVV, PSF, and FiD techniques was 90.0 (62.4-116.8), 92.0 (70.6-123.2), and 89.4 (62.5-126.5) cm/s, respectively. The largest difference in CV estimates was found between PSF and DVV which was caused by smaller CV magnitudes detected only by the DVV technique. Using DVV, a lower CV at BB was found in PAF patients compared with those without atrial fibrillation (AF) [79.1 (72.2-91.2) vs. 88.3 (79.3-97.2) cm/s; P < 0.001].ConclusionsAreas of local conduction heterogeneities were most accurately identified using the DVV technique, whereas PSF and FiD techniques smoothen wavefront propagation thereby masking local areas of conduction slowing. Comparing patients with and without AF, slower wavefront propagation during SR was found at BB in PAF patients, indicating structural remodelling.
| S-EPMC8576284 | biostudies-literature
Project description:Previous studies have shown an association between elevated atrial NADPH-dependent oxidative stress and decreased plasma apelin in patients with atrial fibrillation (AF), though the basis for this relationship is unclear. In the current study, RT-PCR and immunofluorescence studies of human right atrial appendages (RAAs) showed expression of the apelin receptor, APJ, and reduced apelin content in the atria, but not in plasma, of patients with AF versus normal sinus rhythm. Disruption of the apelin gene in mice increased (2.4-fold) NADPH-stimulated superoxide levels and slowed atrial conduction velocities in optical mapping of a Langendorff-perfused isolated heart model, suggesting that apelin levels may influence AF vulnerability. Indeed, in mice with increased AF vulnerability (induced by chronic intense exercise), apelin administration reduced the incidence and duration of induced atrial arrhythmias in association with prolonged atrial refractory periods. Moreover, apelin decreased AF induction in isolated atria from exercised mice while accelerating conduction velocity and increasing action potential durations. At the cellular level, these changes were associated with increased atrial cardiomyocyte sodium currents. These findings support the conclusion that reduced atrial apelin is maladaptive in fibrillating human atrial myocardium and that increasing apelin bioavailability may be a worthwhile therapeutic strategy for treating and preventing AF.
| S-EPMC7526452 | biostudies-literature
Project description:In statistical modelling with Gaussian process regression, it has been shown that combining (few) high-fidelity data with (many) low-fidelity data can enhance prediction accuracy, compared to prediction based on the few high-fidelity data only. Such information fusion techniques for multi-fidelity data commonly approach the high-fidelity model f h(t) as a function of two variables (t, s), and then use f l(t) as the s data. More generally, the high-fidelity model can be written as a function of several variables (t, s 1, s 2….); the low-fidelity model f l and, say, some of its derivatives can then be substituted for these variables. In this paper, we will explore mathematical algorithms for multi-fidelity information fusion that use such an approach towards improving the representation of the high-fidelity function with only a few training data points. Given that f h may not be a simple function-and sometimes not even a function-of f l, we demonstrate that using additional functions of t, such as derivatives or shifts of f l, can drastically improve the approximation of f h through Gaussian processes. We also point out a connection with 'embedology' techniques from topology and dynamical systems. Our illustrative examples range from instructive caricatures to computational biology models, such as Hodgkin-Huxley neural oscillations.
| S-EPMC6501345 | biostudies-literature
Project description:Several lines of evidence have suggested that maintenance of atrial fibrillation (AF) depends on reentrant mechanisms. Maintenance of reentry necessitates a sufficiently short refractory period and/or delayed conduction, and AF has been associated with both alterations. Fibrosis, cellular dysfunction, and gap junction protein alterations occur in AF and cause conduction delay. We performed this study to test the hypothesis that gap junction protein overexpression would improve conduction and prevent AF.Thirty Yorkshire swine were randomized into 2 groups (sinus rhythm and AF), and each group into 3 subgroups: sham-operated control, gene therapy with adenovirus expressing connexin (Cx) 40, and gene therapy with adenovirus expressing Cx43 (n=5 per subgroup). All animals had epicardial gene painting; the AF group had burst atrial pacing. All animals underwent terminal study 7 days after gene transfer. Sinus rhythm animals had strong transgene expression but no atrial conduction changes. In AF animals, controls had reduced and lateralized Cx43 expression, and Cx43 gene transfer restored expression and cellular location to sinus rhythm control levels. In the AF group, both Cx40 and Cx43 gene transfer improved conduction and reduced AF relative to controls.Connexin gene therapy preserved atrial conduction and prevented AF.
| S-EPMC3260348 | biostudies-literature
Project description:Probability density functions (PDFs) are fundamental objects in mathematics with numerous applications in computer vision, machine learning and medical imaging. The feasibility of basic operations such as computing the distance between two PDFs and estimating a mean of a set of PDFs is a direct function of the representation we choose to work with. In this paper, we study the Gaussian mixture model (GMM) representation of the PDFs motivated by its numerous attractive features. (1) GMMs are arguably more interpretable than, say, square root parameterizations (2) the model complexity can be explicitly controlled by the number of components and (3) they are already widely used in many applications. The main contributions of this paper are numerical algorithms to enable basic operations on such objects that strictly respect their underlying geometry. For instance, when operating with a set of K component GMMs, a first order expectation is that the result of simple operations like interpolation and averaging should provide an object that is also a K component GMM. The literature provides very little guidance on enforcing such requirements systematically. It turns out that these tasks are important internal modules for analysis and processing of a field of ensemble average propagators (EAPs), common in diffusion weighted magnetic resonance imaging. We provide proof of principle experiments showing how the proposed algorithms for interpolation can facilitate statistical analysis of such data, essential to many neuroimaging studies. Separately, we also derive interesting connections of our algorithm with functional spaces of Gaussians, that may be of independent interest.
| S-EPMC4816648 | biostudies-literature
Project description:BackgroundConduction velocity (CV) heterogeneity and myocardial fibrosis both promote re-entry, but the relationship between fibrosis as determined by left atrial (LA) late-gadolinium enhanced cardiac magnetic resonance imaging (LGE-CMRI) and CV remains uncertain.ObjectiveAlthough average CV has been shown to correlate with regional LGE-CMRI in patients with persistent AF, we test the hypothesis that a localized relationship exists to underpin LGE-CMRI as a minimally invasive tool to map myocardial conduction properties for risk stratification and treatment guidance.Method3D LA electroanatomic maps during LA pacing were acquired from eight patients with persistent AF following electrical cardioversion. Local CVs were computed using triads of concurrently acquired electrograms and were co-registered to allow correlation with LA wall intensities obtained from LGE-CMRI, quantified using normalized intensity (NI) and image intensity ratio (IIR). Association was evaluated using multilevel linear regression.ResultsAn association between CV and LGE-CMRI intensity was observed at scales comparable to the size of a mapping electrode: -0.11 m/s per unit increase in NI (P < 0.001) and -0.96 m/s per unit increase in IIR (P < 0.001). The magnitude of this change decreased with larger measurement area. Reproducibility of the association was observed with NI, but not with IIR.ConclusionAt clinically relevant spatial scales, comparable to area of a mapping catheter electrode, LGE-CMRI correlates with CV. Measurement scale is important in accurately quantifying the association of CV and LGE-CMRI intensity. Importantly, NI, but not IIR, accounts for changes in the dynamic range of CMRI and enables quantitative reproducibility of the association.
| S-EPMC7693630 | biostudies-literature
Project description:ObjectiveLocal activation time (LAT) mapping of cardiac chambers is vital for targeted treatment of cardiac arrhythmias in catheter ablation procedures. Current methods require too many LAT observations for an accurate interpolation of the necessarily sparse LAT signal extracted from intracardiac electrograms (EGMs). Additionally, conventional performance metrics for LAT interpolation algorithms do not accurately measure the quality of interpolated maps. We propose, first, a novel method for spatial interpolation of the LAT signal which requires relatively few observations; second, a realistic sub-sampling protocol for LAT interpolation testing; and third, a new color-based metric for evaluation of interpolation quality that quantifies perceived differences in LAT maps.MethodsWe utilize a graph signal processing framework to reformulate the irregular spatial interpolation problem into a semi-supervised learning problem on the manifold with a closed-form solution. The metric proposed uses a color difference equation and color theory to quantify visual differences in generated LAT maps.ResultsWe evaluate our approach on a dataset consisting of seven LAT maps from four patients obtained by the CARTO electroanatomic mapping system during premature ventricular complex (PVC) ablation procedures. Random sub-sampling and re-interpolation of the LAT observations show excellent accuracy for relatively few observations, achieving on average 6% lower error than state-of-the-art techniques for only 100 observations.ConclusionOur study suggests that graph signal processing methods can improve LAT mapping for cardiac ablation procedures.SignificanceThe proposed method can reduce patient time in surgery by decreasing the number of LAT observations needed for an accurate LAT map.
| S-EPMC9549513 | biostudies-literature
Project description:The recent decline in energy, size and complexity scaling of traditional von Neumann architecture has resurrected considerable interest in brain-inspired computing. Artificial neural networks (ANNs) based on emerging devices, such as memristors, achieve brain-like computing but lack energy-efficiency. Furthermore, slow learning, incremental adaptation, and false convergence are unresolved challenges for ANNs. In this article we, therefore, introduce Gaussian synapses based on heterostructures of atomically thin two-dimensional (2D) layered materials, namely molybdenum disulfide and black phosphorus field effect transistors (FETs), as a class of analog and probabilistic computational primitives for hardware implementation of statistical neural networks. We also demonstrate complete tunability of amplitude, mean and standard deviation of the Gaussian synapse via threshold engineering in dual gated molybdenum disulfide and black phosphorus FETs. Finally, we show simulation results for classification of brainwaves using Gaussian synapse based probabilistic neural networks.
| S-EPMC6744503 | biostudies-literature