for j in range(3):
            input_vol3d[:, j, :] = peppers(size=512)[0, :, :]
            mask3d[:, j, :] = mask2d