def write_mask_image(self):
"""Show the overlap of each run's mask."""
mask_imgs = [nib.load(f) for f in self.inputs.masks]
mask_data = [img.get_data()[..., np.newaxis] for img in mask_imgs]
each_mask = np.concatenate(mask_data, axis=-1)
any_mask = each_mask.max(axis=-1)
self.every_mask = each_mask.min(axis=-1)
m = Mosaic(self.inputs.anatomy, each_mask, any_mask,
step=self.step, show_mask=False)
m.plot_mask_edges()
out_fname = op.abspath("mask_overlap.png")
self.summary_files.append(out_fname)
m.savefig(out_fname)
m.close()
def write_mask_image(self):
"""Show the overlap of each run's mask."""
mask_imgs = [nib.load(f) for f in self.inputs.masks]
mask_data = [img.get_data()[..., np.newaxis] for img in mask_imgs]
each_mask = np.concatenate(mask_data, axis=-1)
any_mask = each_mask.max(axis=-1)
self.every_mask = each_mask.min(axis=-1)
m = Mosaic(self.inputs.anatomy, each_mask, any_mask,
stat_interp="nearest", step=self.step, show_mask=False)
m.plot_mask_edges()
out_fname = op.abspath("mask_overlap.png")
self.summary_files.append(out_fname)
m.savefig(out_fname)
m.close()
