afni.nimh.nih.gov

There is also program 3dmaskSVD.

And the related (slower) programs 3dLocalPV and 3dLocalSVD -- which give you the eigenvector(s) of the time series from a neighborhood around each voxel.