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