// Now we inverting the masking operation, to go back to a full 3D// representation
(x, y, z) = mean_img.shape
components = np.zeros((x, y, z, n_components))
components[mask] = components_masked
// Using a masked array is important to have transparency in the figures