Add spins per voxel estimation function