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