Use a PSF to find the best matching centres in a cube. Uses noise cube to clip bad pixels.
multi_cent_psf()