Elsevier

NeuroImage

Volume 39, Issue 1, 1 January 2008, Pages 127-135
NeuroImage

Technical Note
Projection of fMRI data onto the cortical surface using anatomically-informed convolution kernels

https://doi.org/10.1016/j.neuroimage.2007.08.039Get rights and content

Abstract

As surface-based data analysis offer an attractive approach for intersubject matching and comparison, the projection of voxel-based 3D volumes onto the cortical surface is an essential problem. We present here a method that aims at producing representations of functional brain data on the cortical surface from functional MRI volumes. Such representations are for instance required for subsequent cortical-based functional analysis. We propose a projection technique based on the definition, around each node of the gray/white matter interface mesh, of convolution kernels whose shape and distribution rely on the geometry of the local anatomy. For one anatomy, a set of convolution kernels is computed that can be used to project any functional data registered with this anatomy. Therefore resulting in anatomically-informed projections of data onto the cortical surface, this kernel-based approach offers better sensitivity, specificity than other classical methods and robustness to misregistration errors. Influences of mesh and volumes spatial resolutions were also estimated for various projection techniques, using simulated functional maps.

Introduction

Functional magnetic resonance images (fMRI) are generally presented and analysed in their original 3D space of acquisition with little regard to the structural features of the human brain. In this context, surface-based analyses of cortical data have gained interest as they advocate for anatomical considerations of the object under study. First, as the main part of the brain activity is generated by cortical pyramidal neurons, descriptions of the cortex should respect its sheet-like structure, despite the intrinsic voxel-based nature of magnetic resonance images and the highly convoluted nature of the cortex. Moreover, distances computed in 3D voxel-based space do not take folds into account, and two points with low Euclidean distance can actually be located on the two opposite sides of a sulcus; with respect to this issue, geodesic distances, along the cortical surface, describe the anatomy better. Several papers in the literature (Clouchoux et al., 2005, Van Essen and Drury, 1997, Fischl et al., 1999, Toro and Burnod, 2003) advocate for this approach. As these works deal with visualisation, intersubjects registration or cortical localisation, little research has also been published on surface-based functional data analysis (Andrade et al., 2001, Flandin et al., 2002, Goebel and Singer, 1999, Kiebel et al., 2000). On this matter, one essential problem concerns the projection of functional data onto the anatomy, in other terms the mapping of functional voxel-based volumes onto cortical triangulated meshes. This is especially a crucial point to open the way to fMRI data analysis confined to the cortical ribbon.

Still, this problem has remained poorly studied and is rarely found in the literature. Some «geometric» simple methods do exist: some methods propose to average intensities along a normal direction or inside a sphere centered at each node of the mesh representing the cortical surface (e.g. when using free package BrainVISA), or to assign each node with the value of its containing voxel (Saad et al., 2004), or to compute trilinear interpolations (Andrade et al., 2001). Other methods try to embed some explicit anatomical information into the process: it can be done by defining each node’s influence scope through voxels distances to and along the cortical mesh (Warnking et al., 2002), or computing a geodesic Voronoï diagram so that each node is integrated within an associated set of voxels defined by the local anatomical geometry (Grova et al., 2006).

From these methods stems a common interest in delineating 3D areas, sometimes overlapping, onto which the signal emanating from the nodes would possibly be dispersed, taking local anatomical features into account. Nevertheless, the problems remains very difficult because of numerous factors such as the volumes resolution, the partial volume effects, the sensitivity to errors from previous processing steps like segmentation, mesh extraction and registration, all of this being added to the highly folded nature of the surface. Therefore, a good interpolation method should show robustness to these various parameters.

In this study, we propose an original method to produce surface-based representations of the cortical activity directly on the cortical mesh, relying on an expected distribution of the BOLD signal in EPI volumes, taking into account as much anatomical information as possible. The method itself is described in the next section. The third section presents experiments we ran to validate our process and evaluate its robustness to segmentation and registration errors. Results are presented in the fourth section and the method is discussed and compared to others in the last section.

Section snippets

Expected distribution of the signal

In order to represent the cortical activity from the whole 3D volume back to the surface where it was originated, we relied on an expected distribution of the BOLD signal around each node of the cortical mesh to address correlations between voxels and mesh nodes, with physiological and image-related motivations. Several distinct phenomena were considered: first, as a consequence of the columnar architecture of neurons inside the cortex (Mountcastle, 1978), we expect each cortical column to show

Simulated data sets

We led a series of experiments in order to assess the benefits of projecting functional data using convolution kernels, in terms of detection power, robustness to various kinds of errors and sensitivity to resolution changes. Most of these experiments were previously suggested by Grova et al. (2006) for the Voronoï-based approach. In this validation frame, several types of meshes were generated. We first extracted cortical meshes from the T1-weighted volume at three different resolution levels:

Results

We projected simulated functional data onto the cortical surfaces at three different resolutions. For each projection, we computed an AUCclose value. The results are presented in Table 1. They tend to illustrate kernel-based projection’s best accuracy in recovering, out from a simulated volume, the surface-based activation map initially taken as ground truth, in comparison to other methods, and for each resolution. With very-low resolution meshes, all the methods give results very close or

Choice of the orthogonality criterion

In order to compute geodesic and normal distances for any voxel in a node’s neighborhood, each voxel is assigned what was referred to as a normal node on the surface, from which the two distances, then the weights, were computed. These associations result from a Fast-Marching propagation algorithm. Relating voxels to surface nodes using geodesic distances is a relevant way to deal with the folded nature of cortex, rather than euclidean distances. However, the original orthogonality criterion we

Conclusion

In this paper, we propose a method allowing the projection of functional images onto the cortical surface. Relying on physiological and image-related hypotheses, it determines, for each mesh node, a specific interpolation area, strongly influenced by local anatomy, within which voxels intensities are averaged. Applying it to activation maps can serve visualisation purposes, as well as cortical localisation of activation foci through the use of a surface-based coordinate system (Clouchoux et

References (20)

There are more references available in the full text version of this article.

Cited by (23)

  • New acquisition techniques and their prospects for the achievable resolution of fMRI

    2021, Progress in Neurobiology
    Citation Excerpt :

    While not the focus of this review, it is important to note that not only the data acquisition, but also the data analysis greatly impacts the spatial specificity. The data acquired in fMRI studies are subject to complex analysis pipelines, in which they are interpolated (Glasser et al., 2013), projected (Operto et al., 2008), pooled (Polimeni et al., 2018), cleaned (Murphy and Bright, 2017) and smoothed (Blazejewska et al., 2019). This can drastically reduce the effective resolution of the data after the analysis.

  • Corresponding anatomical and coactivation architecture of the human precuneus showing similar connectivity patterns with macaques

    2019, NeuroImage
    Citation Excerpt :

    The primary limitations of the surface based analysis are missing the subcortical information and requiring accurate projection from fMRI to structure MRI (Tucholka et al., 2012). The accurate projection, in view of the typical resolution of fMRI data, from functional data onto the cortical surface is a real challenge in practice due to small cortical thickness, small distance between the width of a gyrus or the banks of a given sulci, and the artefacts due to geometric distortions that are not systematically and fully corrected during pre-processing (Hutton et al., 2002; Operto et al., 2008; Tucholka et al., 2012). Moreover, Tucholka et al. (2012) reported the incorrect projection of a computation task for surface-based analysis.

  • Analysis strategies for high-resolution UHF-fMRI data

    2018, NeuroImage
    Citation Excerpt :

    Averaging voxels along this straight-line path will cause sampling biases in small curved regions of cortex, where paths from multiple vertices tend to converge, and curvilinear paths that bend with the cortex may be beneficial (see Predicting geometry of columns and layers). Anatomically-informed interpolation has also been proposed, where the weights of the interpolation kernel are adapted to the local geometry of the cortex to avoid noise contamination from extracerebral CSF and signal dilution from subjacent white matter (Grova et al., 2006; Operto et al., 2008). It is also possible to perform anatomically informed analysis using cortical surface models without projecting the fMRI data onto the surfaces by performing the inverse operation of “projecting the surface vertices into voxels”.

  • Anatomically-adapted graph wavelets for improved group-level fMRI activation mapping

    2015, NeuroImage
    Citation Excerpt :

    Previous studies following CSM approaches (Andrade et al., 2001; Chung et al., 2005; Hagler et al., 2006; Qiu et al., 2006), which anatomically constrain the activation mapping, suffer from a few shortcomings. First, the interpolation effects related to the projection of 3-D volumes onto the cortical surface is a major challenge (Grova et al., 2006; Operto et al., 2008). Second, smoothing leads to spatial inaccuracy and artificial shift of activations, be it using conventional volume-based approaches (e.g., SPM) or CSM (Jo et al., 2007, 2008; Mikl et al., 2008).

  • An empirical comparison of surface-based and volume-based group studies in neuroimaging

    2012, NeuroImage
    Citation Excerpt :

    More generally, besides the fact that non-cortical regions are not analyzed by a surface-based procedure, this analysis suffers from an important limitation: it requires that functional data has been accurately projected on the cortical surface, i.e. that the correspondence of functional data with anatomical data is correct. This is a real challenge in practice, as the cortical thickness, and also the distance between the banks of a given sulci or the width of a gyrus are small, in view of the typical resolution of fMRI data (Operto et al., 2008). Moreover, fMRI data contain artefacts due to geometric distortions related to EPI acquisition.

View all citing articles on Scopus
View full text