def __init__(self,
model,
contrast,
corrected=True,
stat_thresh=2.3,
stat_range=(2, 5),
stat_cmap="OrRd_r",
stat_alpha=.85,
sharp=False,
label_slices=True):
"""Initialize the object but do not plot anything yet."""
anat_img = nib.load(op.join(data_dir, "average_anat.nii.gz"))
contrast_dir = op.join(analysis_dir, exp + "-" + model, "group", "mni",
contrast)
if corrected:
stat_file = "zstat1_threshold.nii.gz"
else:
stat_file = "zstat1.nii.gz"
stat_img = nib.load(op.join(contrast_dir, stat_file))
mask_img = nib.load(op.join(contrast_dir, "mask.nii.gz"))
self.anat_img = VolumeImg(anat_img.get_data(), anat_img.get_affine(),
"mni")
self.stat_img = VolumeImg(stat_img.get_data(), stat_img.get_affine(),
"mni")
self.mask_img = VolumeImg(mask_img.get_data(),
mask_img.get_affine(),
"mni",
interpolation="nearest")
self.stat_thresh = stat_thresh
self.stat_cmap = stat_cmap
self.stat_alpha = stat_alpha
self.stat_vmin, self.stat_vmax = stat_range
self.label_slices = label_slices
self.sharp = sharp
