Xavier Tricoche

Influence of tissue conductivity anisotropy on EEG/MEG field and return current computation in a realistic head model: a simulation and visualization study using high-resolution finite element modeling.

To achieve a deeper understanding of the brain, scientists, and clinicians use electroencephalography (EEG) and magnetoencephalography (MEG) inverse methods to reconstruct sources in the cortical sheet of the human brain. The influence of structural and electrical anisotropy in both the skull and the white matter on the EEG and MEG source reconstruction is not well understood. In this paper, we report on a study of the sensitivity to tissue anisotropy of the EEG/MEG forward problem for deep and superficial neocortical sources with differing orientation components in an anatomically accurate model of the human head. The goal of the study was to gain insight into the effect of anisotropy of skull and white matter conductivity through the visualization of field distributions, isopotential surfaces, and return current flow and through statistical error measures. One implicit premise of the study is that factors that affect the accuracy of the forward solution will have at least as strong an influence over solutions to the associated inverse problem. Major findings of the study include (1) anisotropic white matter conductivity causes return currents to flow in directions parallel to the white matter fiber tracts; (2) skull anisotropy has a smearing effect on the forward potential computation; and (3) the deeper a source lies and the more it is surrounded by anisotropic tissue, the larger the influence of this anisotropy on the resulting electric and magnetic fields. Therefore, for the EEG, the presence of tissue anisotropy both for the skull and white matter compartment substantially compromises the forward potential computation and as a consequence, the inverse source reconstruction. In contrast, for the MEG, only the anisotropy of the white matter compartment has a significant effect. Finally, return currents with high amplitudes were found in the highly conducting cerebrospinal fluid compartment, underscoring the need for accurate modeling of this space.

Efficient Computation and Visualization of Coherent Structures in Fluid Flow Applications

The recently introduced notion of Finite-Time Lyapunov Exponent to characterize Coherent Lagrangian Structures provides a powerful framework for the visualization and analysis of complex technical flows. Its definition is simple and intuitive, and it has a deep theoretical foundation. While the application of this approach seems straightforward in theory, the associated computational cost is essentially prohibitive. Due to the Lagrangian nature of this technique, a huge number of particle paths must be computed to fill the space-time flow domain. In this paper, we propose a novel scheme for the adaptive computation of FTLE fields in two and three dimensions that significantly reduces the number of required particle paths. Furthermore, for three-dimensional flows, we show on several examples that meaningful results can be obtained by restricting the analysis to a well-chosen plane intersecting the flow domain. Finally, we examine some of the visualization aspects of FTLE-based methods and introduce several new variations that help in the analysis of specific aspects of a flow.

A topology simplification method for 2D vector fields

Topology analysis of plane, turbulent vector fields results in visual clutter caused by critical points indicating vortices of finer and finer scales. A simplification can be achieved by merging critical points within a prescribed radius into higher order critical points. After building clusters containing the singularities to merge, the method generates a piecewise linear representation of the vector field in each cluster containing only one (higher order) singularity. Any visualization method can be applied to the result after this process. Using different maximal distances for the critical points to be merged results in a hierarchy of simplified vector fields that can be used for analysis on different scales.

Continuous topology simplification of planar vector fields

Vector fields can present complex structural behavior, especially in turbulent computational fluid dynamics. The topological analysis of these data sets reduces the information, but one is usually still left with too many details for interpretation. In this paper, we present a simplification approach that removes pairs of critical points from the data set, based on relevance measures. In contrast to earlier methods, no grid changes are necessary, since the whole method uses small local changes of the vector values defining the vector field. An interpretation in terms of bifurcations underlines the continuous, natural flavor of the algorithm.

Surface techniques for vortex visualization

This paper presents powerful surface based techniques for the analysis of complex flow fields resulting from CFD simulations. Emphasis is put on the examination of vortical structures. An improved method for stream surface computation that delivers accurate results in regions of intricate flow is presented, along with a novel method to determine boundary surfaces of vortex cores. A number of surface techniques are presented that aid in understanding the flow behavior displayed by these surfaces. Furthermore, a scheme for phenomenological extraction of vortex core lines using stream surfaces is discussed and its accuracy is compared to one of the most established standard techniques.

Topology Tracking for the Visualization of Time-Dependent Two-Dimensional Flows

The paper presents a topology-based visualization method for time-dependent two-dimensional vector elds. A time interpolation enables the accurate tracking of critical points and closed orbits as well as the detection and identication of structural changes. This completely characterizes the topology of the unsteady ow. Bifurcation theory provides the theoretical framework. The results are conveyed by surfaces that separate subvolumes of uniform ow behavior in a three-dimensional space-time domain.

Generation of Accurate Integral Surfaces in Time-Dependent Vector Fields

We present a novel approach for the direct computation of integral surfaces in time-dependent vector fields. As opposed to previous work, which we analyze in detail, our approach is based on a separation of integral surface computation into two stages: surface approximation and generation of a graphical representation. This allows us to overcome several limitations of existing techniques. We first describe an algorithm for surface integration that approximates a series of time lines using iterative refinement and computes a skeleton of the integral surface. In a second step, we generate a well-conditioned triangulation. Our approach allows a highly accurate treatment of very large time-varying vector fields in an efficient, streaming fashion. We examine the properties of the presented methods on several example datasets and perform a numerical study of its correctness and accuracy. Finally, we investigate some visualization aspects of integral surfaces.

Tracking of Vector Field Singularities in Unstructured 3D Time-Dependent Datasets

We present an approach for monitoring the positions of vector field singularities and related structural changes in time-dependent datasets. The concept of singularity index is discussed and extended from the well-understood planar case to the more intricate three-dimensional setting. Assuming a tetrahedral grid with linear interpolation in space and time, vector field singularities obey rules imposed by fundamental invariants (Poincare index), which we use as a basis for an efficient tracking algorithm. We apply the presented algorithm to CFD datasets to illustrate its purpose. We examine structures that exhibit topological variations with time and describe some of the insight gained with our method. Examples are given that show a correlation in the evolution of physical quantities that play a role in vortex breakdown.

Delineating white matter structure in diffusion tensor MRI with anisotropy creases.

Geometric models of white matter architecture play an increasing role in neuroscientific applications of diffusion tensor imaging, and the most popular method for building them is fiber tractography. For some analysis tasks, however, a compelling alternative may be found in the first and second derivatives of diffusion anisotropy. We extend to tensor fields the notion from classical computer vision of ridges and valleys, and define anisotropy creases as features of locally extremal tensor anisotropy. Mathematically, these are the loci where the gradient of anisotropy is orthogonal to one or more eigenvectors of its Hessian. We propose that anisotropy creases provide a basis for extracting a skeleton of the major white matter pathways, in that ridges of anisotropy coincide with interiors of fiber tracts, and valleys of anisotropy coincide with the interfaces between adjacent but distinctly oriented tracts. The crease extraction algorithm we present generates high-quality polygonal models of crease surfaces, which are further simplified by connected-component analysis. We demonstrate anisotropy creases on measured diffusion MRI data, and visualize them in combination with tractography to confirm their anatomic relevance.

Visualization of Intricate Flow Structures for Vortex Breakdown Analysis

Vortex breakdowns and flow recirculation are essential phenomena in aeronautics where they appear as a limiting factor in the design of modern aircrafts. Because of the inherent intricacy of these features, standard flow visualization techniques typically yield cluttered depictions. The paper addresses the challenges raised by the visual exploration and validation of two CFD simulations involving vortex breakdown. To permit accurate and insightful visualization we propose a new approach that unfolds the geometry of the breakdown region by letting a plane travel through the structure along a curve. We track the continuous evolution of the associated projected vector field using the theoretical framework of parametric topology. To improve the understanding of the spatial relationship between the resulting curves and lines we use direct volume rendering and multidimensional transfer functions for the display of flow-derived scalar quantities. This enriches the visualization and provides an intuitive context for the extracted topological information. Our results offer clear, synthetic depictions that permit new insight into the structural properties of vortex breakdowns.