Plot dropout regions on the flatmap

A dropout region is a region with very low EPI signal. In pycortex a crosshatch is used to display such dropout regions.

The crosshatches are created using the reference nifti image file with a threshold. Setting the with_dropout=True parameter in quickflat.make_figure takes the reference nifti image file and computes a thresholded version of this using the following formula:

FIXME: .. code-block:: python

rawdata[rawdata==0] = np.mean(rawdata[rawdata!=0]) normdata = (rawdata - rawdata.min()) / (rawdata.max() - rawdata.min()) normdata = (1 - normdata) ** power

Out:

Failed to get connection
** (inkscape:5219): CRITICAL **: 01:27:52.316: dbus_g_proxy_new_for_name: assertion 'connection != NULL' failed
** (inkscape:5219): CRITICAL **: 01:27:52.317: dbus_g_proxy_call: assertion 'DBUS_IS_G_PROXY (proxy)' failed
** (inkscape:5219): CRITICAL **: 01:27:52.317: dbus_g_connection_register_g_object: assertion 'connection != NULL' failed

import cortex
import numpy as np
np.random.seed(1234)

# Create a random pycortex Volume
volume = cortex.Volume.random(subject='S1', xfmname='fullhead')

# Plot a flatmap with the data projected onto the surface
# Highlight the curvature and dropout regions
_ = cortex.quickflat.make_figure(volume,
                                 with_curvature=True,
                                 with_dropout=True)