anat *= np.sin(4*x)*np.sin(4*y)*np.sin(4*z) activation[activation < -.5] = -.6 activation -= activation.min() main_img = ni_api.Image( anat, ni_api.AffineTransform.from_start_step('ijk', xipy_ras, [0,0,0], [1,1,1]) ) over_img = ni_api.Image( activation, ni_api.AffineTransform.from_start_step('ijk', xipy_ras, [0,0,0], [1,1,1]) ) bi = BlendedImages(vtk_order=True) bi.main = main_img bi.over = over_img pd = tvtk.PointData() class PDataDB(object): def __init__(self, pd): self.pd = pd self.db = dict() def add_to_pdata(self, array, name): if len(array.shape) > 3: vtk_order = quick_convert_rgba_to_vtk(array) flat_shape = (np.prod(vtk_order.shape[:3]), vtk_order.shape[3]) else: