Made it into a function
Working with real data
selecting specific voxels and doing economic SVD