Visual maps from anatomy

Overview

Benson and his colleagues (including Aguirre and Brainard) developed methods to identify the positions of the early visual field maps from anatomical (T1-weighted) MRI data. They can achieve this because the detailed properties of the cortical folds are well-correlated with the positions of these functional maps. The general principle - correlation between structural folds and functional regions - is an emerging theme from work in occipital and temporal cortex (e.g., see the papers by Weiner on face and Witthoft on hV4 in the (Grill-Spector lab)[http://vpnl.stanford.edu/].

The docker container (also a Flywheel Gear) that we describe here implements Noah Benson's algorithm for identifying V1/V2/V3 from an anatomical T1.

It also returns the assignments of all the visual field maps represented in the Wang et al. atlas onto the T1-weighted image.

##Input

T1 anatomy nifti file.

##Outputs

More detailed information for the docker outputs can be found here. Below we list the files the most relevant files:

scanner.template_areas.nii.gz - A nifti file at the same resolution and in the same space as the T1 anatomical file. The value at each voxel indicates the visual area. For example, the visual area template will consist of values ranging from 0-3 (1: V1, 2: V2, 3: V3, 0: none of the above).

scanner.wang2015_atlas.nii.gz - A nifti file at the same resolution and in the same space as the T1 anatomical file. The value at each voxel indicates the visual area from the Wang et al. atlas [http://cercor.oxfordjournals.org/content/early/2016/05/25/cercor.bhw157.full]. Values range between 0 to 25. 0 indicates voxels that are not labeled as a visual area.

##TODO

A few comments on things we should change.

• The word 'scanner' may be misleading. The output files are in the same coordinate frame as the T1-weighted input. Those data may or may not be in scanner (image) coordinates. Indeed, we often run the container on T1-weighted data that have been registered into AC-PC space.
• We should probably regularize 'wang2015_atlas' and 'template_areas' to both be, say, 'wang2015_atlas' and benson2014_atlas'. We think that would be helpful.
• The algorithm also returns gifti files. These can be downloaded and viewed with a gifti viewer (e.g. [http://www.artefact.tk/software/matlab/gifti/].

The conversion from MGZ to GIFTI seems to be broken in Flywheel. Let's see if we can't figure out what the deal is.

##Visualizing

In Flywheel simply select the visualization icon and choose Papaya or Slice Drop. The color map "Spectrum" works well in Papaya.

##References

The retinotopic organization of striate cortex is well predicted by surface topology. Benson NC, Butt OH, Datta R, Radoeva PD, Brainard DH, Aguirre GK. Current Biology. 2012;22(21):2081–2085. doi: 10.1016/j.cub.2012.09.014. pmid:23041195

Correction of Distortion in Flattened Representations of the Cortical Surface Allows Prediction of V1-V3 Functional Organization from Anatomy Noah C. Benson, Omar H. Butt, David H. Brainard, Geoffrey K. Aguirre

Wang L, Mruczek REB, Arcaro MJ, Kastner S. Probabilistic Maps of Visual Topography in Human Cortex. Cerebral Cortex (New York, NY). 2015;25(10):3911-3931. doi:10.1093/cercor/bhu277.

Witthoft N, Nguyen ML, Golarai G, et al. Where Is Human V4? Predicting the Location of hV4 and VO1 from Cortical Folding. Cerebral Cortex (New York, NY). 2014;24(9):2401-2408. doi:10.1093/cercor/bht092.