Biomedical subjects

Olivier Faugeras

Publications and source records attributed to Olivier Faugeras.

10 recordsLinked to original sources

PubMed2006

Using nonlinear models in fMRI data analysis: model selection and activation detection.

There is an increasing interest in using physiologically plausible models in fMRI analysis. These models do raise new mathematical problems in terms of parameter estimation and interpretation of the measured data. In this paper, we show how to use physiological models to map and analyze brain activity from fMRI data. We describe a maximum likelihood parameter estimation algorithm and a statistical test that allow the following two actions: selecting the most statistically significant hemodynamic model for the measured data and deriving activation maps based on such model. Furthermore, as parameter estimation may leave much incertitude on the exact values of parameters, model identifiability characterization is a particular focus of our work. We applied these methods to different variations of the Balloon Model (Buxton, R.B., Wang, E.C., and Frank, L.R. 1998. Dynamics of blood flow and oxygenation changes during brain activation: the balloon model. Magn. Reson. Med. 39: 855-864; Buxton, R.B., Uludağ, K., Dubowitz, D.J., and Liu, T.T. 2004. Modelling the hemodynamic response to brain activation. NeuroImage 23: 220-233; Friston, K. J., Mechelli, A., Turner, R., and Price, C. J. 2000. Nonlinear responses in fMRI: the balloon model, volterra kernels, and other hemodynamics. NeuroImage 12: 466-477) in a visual perception checkerboard experiment. Our model selection proved that hemodynamic models better explain the BOLD response than linear convolution, in particular because they are able to capture some features like poststimulus undershoot or nonlinear effects. On the other hand, nonlinear and linear models are comparable when signals get noisier, which explains that activation maps obtained in both frameworks are comparable. The tools we have developed prove that statistical inference methods used in the framework of the General Linear Model might be generalized to nonlinear models.

Adult↗

PubMed2006

Generalized head models for MEG/EEG: boundary element method beyond nested volumes.

Accurate geometrical models of the head are necessary for solving the forward and inverse problems of magneto- and electro-encephalography (MEG/EEG). Boundary element methods (BEMs) require a geometrical model describing the interfaces between different tissue types. Classically, head models with a nested volume topology have been used. In this paper, we demonstrate how this constraint can be relaxed, allowing us to model more realistic head topologies. We describe the symmetric BEM for this new model. The symmetric BEM formulation uses both potentials and currents at the interfaces as unknowns and is in general more accurate than the alternative double-layer formulation.

Algorithms↗

PubMed2006

Bifurcation analysis of Jansen's neural mass model.

We present a mathematical model of a neural mass developed by a number of people, including Lopes da Silva and Jansen. This model features three interacting populations of cortical neurons and is described by a six-dimensional nonlinear dynamical system. We address some aspects of its behavior through a bifurcation analysis with respect to the input parameter of the system. This leads to a compact description of the oscillatory behaviors observed in Jansen and Rit (1995) (alpha activity) and Wendling, Bellanger, Bartolomei, and Chauvel (2000) (spike-like epileptic activity). In the case of small or slow variation of the input, the model can even be described as a binary unit. Again using the bifurcation framework, we discuss the influence of other parameters of the system on the behavior of the neural mass model.

Action Potentials↗

PubMed2005

Fast multipole acceleration of the MEG/EEG boundary element method.

The accurate solution of the forward electrostatic problem is an essential first step before solving the inverse problem of magneto- and electroencephalography (MEG/EEG). The symmetric Galerkin boundary element method is accurate but cannot be used for very large problems because of its computational complexity and memory requirements. We describe a fast multipole-based acceleration for the symmetric boundary element method (BEM). It creates a hierarchical structure of the elements and approximates far interactions using spherical harmonics expansions. The accelerated method is shown to be as accurate as the direct method, yet for large problems it is both faster and more economical in terms of memory consumption.

Algorithms↗

PubMed2005

A common formalism for the integral formulations of the forward EEG problem.

The forward electroencephalography (EEG) problem involves finding a potential V from the Poisson equation inverted Delta x (sigma inverted Delta V) f, in which f represents electrical sources in the brain, and sigma the conductivity of the head tissues. In the piecewise constant conductivity head model, this can be accomplished by the boundary element method (BEM) using a suitable integral formulation. Most previous work uses the same integral formulation, corresponding to a double-layer potential. In this paper we present a conceptual framework based on a well-known theorem (Theorem 1) that characterizes harmonic functions defined on the complement of a bounded smooth surface. This theorem says that such harmonic functions are completely defined by their values and those of their normal derivatives on this surface. It allows us to cast the previous BEM approaches in a unified setting and to develop two new approaches corresponding to different ways of exploiting the same theorem. Specifically, we first present a dual approach which involves a single-layer potential. Then, we propose a symmetric formulation, which combines single- and double-layer potentials, and which is new to the field of EEG, although it has been applied to other problems in electromagnetism. The three methods have been evaluated numerically using a spherical geometry with known analytical solution, and the symmetric formulation achieves a significantly higher accuracy than the alternative methods. Additionally, we present results with realistically shaped meshes. Beside providing a better understanding of the foundations of BEM methods, our approach appears to lead also to more efficient algorithms.

Algorithms↗

PubMed2004

Feature characterization in fMRI data: the Information Bottleneck approach.

Clustering is a well-known technique for the analysis of Functional Magnetic Resonance Imaging (fMRI) data, whose main advantage is certainly flexibility: given a metric on the dataset, it "summarizes" the main characteristics of the data. But intrinsic to this approach are also the problem of defining correctly the quantization accuracy, and the number of clusters necessary to describe the data. The Information Bottleneck (IB) approach to vector quantization, proposed by Bialek and Tishby, addresses these difficulties: (1) it deals with an explicit trade-off between quantization and data fidelity; (2) it does so during the clustering procedure and not post hoc; (3) it takes into account the full statistical distribution of the features within the feature space and not only their most likely value; last, it is principled through an information theoretic formulation, which is relevant in many situations. In this paper, we present how to benefit from this method to analyze fMRI data. Our application is the clustering of voxels according to the magnitude of their responses to several experimental conditions. The IB quantization provides a consistent representation of the data, allowing for an easy interpretation and comparison of datasets.

Algorithms↗

PubMed2004

Variational, geometric, and statistical methods for modeling brain anatomy and function.

We survey the recent activities of the Odyssée Laboratory in the area of the application of mathematics to the design of models for studying brain anatomy and function. We start with the problem of reconstructing sources in MEG and EEG, and discuss the variational approach we have developed for solving these inverse problems. This motivates the need for geometric models of the head. We present a method for automatically and accurately extracting surface meshes of several tissues of the head from anatomical magnetic resonance (MR) images. Anatomical connectivity can be extracted from diffusion tensor magnetic resonance images but, in the current state of the technology, it must be preceded by a robust estimation and regularization stage. We discuss our work based on variational principles and show how the results can be used to track fibers in the white matter (WM) as geodesics in some Riemannian space. We then go to the statistical modeling of functional magnetic resonance imaging (fMRI) signals from the viewpoint of their decomposition in a pseudo-deterministic and stochastic part that we then use to perform clustering of voxels in a way that is inspired by the theory of support vector machines and in a way that is grounded in information theory. Multimodal image matching is discussed next in the framework of image statistics and partial differential equations (PDEs) with an eye on registering fMRI to the anatomy. The paper ends with a discussion of a new theory of random shapes that may prove useful in building anatomical and functional atlases.

Algorithms↗

PubMed2003

The retinotopic organization of primate dorsal V4 and surrounding areas: A functional magnetic resonance imaging study in awake monkeys.

Using functional magnetic resonance imaging (fMRI), we mapped the retinotopic organization throughout the visual cortex of fixating monkeys. The retinotopy observed in areas V1, V2, and V3 was completely consistent with the classical view. V1 and V3 were bordered rostrally by a vertical meridian representation, and V2 was bordered by a horizontal meridian. More anterior in occipital cortex, both areas V3A and MT-V5 had lower and upper visual field representations split by a horizontal meridian. The rostral border of dorsal V4 was characterized by the gradual transition of a representation of the vertical meridian (dorsally) to a representation of the horizontal meridian (more ventrally). Central and ventral V4, on the other hand, were rostrally bordered by a representation of the horizontal meridian. The eccentricity lines ran perpendicular to the ventral V3-V4 border but were parallel to the dorsal V3-V4 border. These results indicate different retinotopic organizations within dorsal and ventral V4, suggesting that the latter regions may not be merely the lower and upper visual field representations of a single area. Moreover, because the present fMRI data are in agreement with previously published electrophysiological results, reported distinctions in the retinotopic organization of human and monkey dorsal V4 reflect genuine species differences that cannot be attributed to technical confounds. Finally, aside from dorsal V4, the retinotopic organization of macaque early visual cortex (V1, V2, V3, V3A, and ventral V4) is remarkably similar to that observed in human fMRI studies. This finding indicates that early visual cortex is mostly conserved throughout hominid evolution.

Anesthesia↗

PubMed2003

Symmetric BEM formulation for the M/EEG forward problem.

The forward M/EEG problem consists in simulating the electric potential and the magnetic field produced outside the head by currents in the brain related to neural activity. All previously proposed solutions using the Boundary Element Method (BEM) were based on a double-layer integral formulation. We have developed an alternative symmetric BEM formulation, achieving a significantly higher accuracy for sources close to tissue interfaces, namely in the cortex. Numerical experiments using a spherical semi-realistic multilayer head model with a known analytical solution are presented, showing that the new BEM performs better than the formulations used in our earlier comparisons, and in most cases outperforms the Finite Element Method (FEM) as far as accuracy is concerned, thus making the BEM a viable choice.

Algorithms↗

PubMed2003

Dynamical components analysis of fMRI data through kernel PCA.

In parallel with standard model-based methods for the analysis of fMRI data, exploratory methods--such as PCA, ICA, and clustering--have been developed to give an account of the dataset with minimal priors: no assumption is made on the data content itself, but the data structure is assumed to show some properties (decorrelation, independence) that allow for the detection of structures of interest. In this paper, we present an alternative that tries to take into account some relevant knowledge for the analysis of the dataset, e.g., the experimental paradigm, while keeping the flexibility of exploratory methods: we use a prior temporal modeling of the data that characterizes each voxel time course. Two implementations are proposed: one based on the General Linear Model, the other one on more flexible short-term predictors, whose complexity is controlled by a Minimum Description Length approach. However, our main concern here is the construction of a multivariate model; the latter is performed with the help of a kernel PCA method that builds a redundant representation of the data through the nonlinearity of the kernel. This allows for a refinement in the description of the (temporal) patterns of interest. In particular, this helps in the characterization of subtle variations in the response to different experimental conditions. We illustrate the usefulness of nonlinearity through the analysis of a synthetic dataset and show on a real dataset how it helps to interpret the experimental results.

Algorithms↗