Commit ba0d9355 authored by mwaskom's avatar mwaskom

Improve resting example colormap and commentary

parent cf9ccbe3
"""
Display Volume Data
===================
Display Resting-State Correlations
==================================
PySurfer provides a function to sample a volume-encoded file
to the surface for improved visualization. At the moment, this
uses Freesurfer's mri_vol2surf routine.
In this example, we show how to build up a complex visualization of a
volume-based image showing resting-state correlations across the whole brain
from a seed in the angular gyrus. We'll plot several views of both hemispheres
in a single window and manipulate the colormap to best represent the nature of
the data.
"""
print __doc__
import os
from surfer import Brain, io
"""Bring up the visualization"""
......@@ -17,28 +20,47 @@ brain = Brain("fsaverage", "split", "inflated", views=['lat', 'med'],
"""Project the volume file and return as an array"""
mri_file = "example_data/resting_corr.nii.gz"
reg_file = "example_data/register.dat"
reg_file = os.path.join(os.environ["FREESURFER_HOME"],
"average/mni152.register.dat")
surf_data_lh = io.project_volume_data(mri_file, "lh", reg_file)
surf_data_rh = io.project_volume_data(mri_file, "rh", reg_file)
"""
You can pass this array to the add_overlay method for
a typical activation overlay (with thresholding, etc.)
You can pass this array to the add_overlay method for a typical activation
overlay (with thresholding, etc.).
"""
brain.add_overlay(surf_data_lh, min=.3, max=.7, name="ang_corr_lh", hemi='lh')
brain.add_overlay(surf_data_rh, min=.3, max=.7, name="ang_corr_rh", hemi='rh')
"""
You can also pass it to add_data for more control
over the visualzation. Here we'll plot the whole
over the visualization. Here we'll plot the whole
range of correlations
"""
for overlay in brain.overlays_dict["ang_corr_lh"]:
overlay.remove()
for overlay in brain.overlays_dict["ang_corr_rh"]:
overlay.remove()
brain.add_data(surf_data_lh, -.7, .7, colormap="jet", alpha=.7, hemi='lh')
brain.add_data(surf_data_rh, -.7, .7, colormap="jet", alpha=.7, hemi='rh')
"""
We want to use an appropriate color map for these data: a divergent map that
is centered on 0, which is a meaningful transition-point as it marks the change
from negative correlations to positive correlations.
We'll also plot the map with some transparency so that we can see through to the
underlying anatomy.
"""
brain.add_data(surf_data_lh, -.7, .7, colormap="RdBu", alpha=.75, hemi='lh')
brain.add_data(surf_data_rh, -.7, .7, colormap="RdBu", alpha=.75, hemi='rh')
"""
We'll next flip the colormap so negative correlations are blue and positive
correlations are red, which is better aligned with the conventions in fMRI
visualization.
"""
for hemi in ["lh", "rh"]:
for cbar in brain.data_dict[hemi]["colorbars"]:
cbar.reverse_lut = True
"""
This overlay represents resting-state correlations with a
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment