We introduce multigrid priors to construct a Bayesian-inspired method to asses brain activity in functional magnetic resonance imaging (fMRI). A sequence of different scale grids is constructed over the image. Starting from the finest scale, coarse grain data variables are sequentially defined for each scale. Then we move back to finer scales, determining for each coarse scale a set of posterior probabilities. The posterior on a coarse scale is used as the prior for activity at the next finer scale. To test the method, we use a linear model with a given hemodynamic response function to construct the likelihood. We apply the method both to real and simulated data of a boxcar experiment. To measure the number of errors, we impose a decision to determine activity by setting a threshold on the posterior. Receiver operating characteristic (ROC) curves are used to study the dependence on threshold and on a few hyperparameters in the relation between specificity and sensitivity. We also study the deterioration of the results for real data, under information loss. This is done by decreasing the number of images in each period and also by decreasing the signal to noise ratio and compare the robustness to other methods.
Download full-text PDF |
Source |
---|---|
http://dx.doi.org/10.1016/j.neuroimage.2004.06.011 | DOI Listing |
IEEE Trans Vis Comput Graph
March 2024
We propose a new barrier-based box-constrained convex QP solver based on a primal-dual interior point method to efficiently solve large-scale pressure Poisson problems with non-negative pressure constraints, which commonly arise in liquid animation. The performance of prior active-set-based approaches is limited by the need to repeatedly update the active set. Our solver eliminates this issue by entirely avoiding the use of an active set, which in turn makes the inner problems of our Newton iteration process fully unconstrained.
View Article and Find Full Text PDFNumer Math (Heidelb)
November 2022
TU Wien, Institute of Analysis and Scientific Computing, Wiedner Hauptstr. 8-10/E101/4, Vienna, 1040 Austria.
We consider a linear symmetric and elliptic PDE and a linear goal functional. We design and analyze a goal-oriented adaptive finite element method, which steers the adaptive mesh-refinement as well as the approximate solution of the arising linear systems by means of a contractive iterative solver like the optimally preconditioned conjugate gradient method or geometric multigrid. We prove linear convergence of the proposed adaptive algorithm with optimal algebraic rates.
View Article and Find Full Text PDFBrain Topogr
January 2020
School of Physics and Astronomy, Sir Peter Mansfield Imaging Centre, University of Nottingham, Nottingham, NG7 2RD, UK.
A previously introduced Bayesian non-parametric multi-scale technique, called iterated Multigrid Priors (iMGP) method, is used to map the topographic organization of human primary somatosensory cortex (S1). We analyze high spatial resolution fMRI data acquired at ultra-high field (UHF, 7T) in individual subjects during vibrotactile stimulation applied to each distal phalange of the left hand digits using both a travelling-wave (TW) and event-related (ER) paradigm design. We compare the somatotopic digit representations generated in S1 using the iMGP method with those obtained using established fMRI paradigms and analysis techniques: Fourier-based analysis of travelling-wave data and General Linear Model (GLM) analysis of event-related data.
View Article and Find Full Text PDFBrain Topogr
March 2014
Sir Peter Mansfield Magnetic Resonance Centre, University of Nottingham, Nottingham, UK,
An important interest in event-related single trial fMRI is the possibility of studying cognitive processes that vary in time (e.g. learning or adaptation).
View Article and Find Full Text PDFIEEE Trans Image Process
November 2013
An image processing observational technique for the stereoscopic reconstruction of the waveform of oceanic sea states is developed. The technique incorporates the enforcement of any given statistical wave law modeling the quasi-Gaussianity of oceanic waves observed in nature. The problem is posed in a variational optimization framework, where the desired waveform is obtained as the minimizer of a cost functional that combines image observations, smoothness priors and a weak statistical constraint.
View Article and Find Full Text PDFEnter search terms and have AI summaries delivered each week - change queries or unsubscribe any time!