Currently we are simply building a nifti image with the data and our transformation matrix and let the image to do that decomposition for us. At some point we should compute that information without using the nifti image (we just don't know how to do it right now). Maybe we can use the nifti image for testing (compare our future implementation against it).