I think this would be better-named ``subtract_psf`` - "removing" is a bit more ambiguous, and `prf`->`psf` because it's called "psf_photometry" (and the module is named `psf.py`, not `prf.py`).