In radio astronomy it is very common to have bizarrely shaped PSF's (called beam patterns), and for several decades they have been using an algorithm called CLEAN to clean up radio survey maps to make them appear as they would had one used a perfect, more desirable, beam pattern.
When I blow up the images I process with PI, I see all the minor aberrations that remain, guiding errors, and so on. I'm wondering if there is a similar technique of subtracting out actual PSF's till the background residual image, much like MRT at large scale sizes, remains, then add back a better looking PSF?
The technique is iterative and only subtracts out some fraction of PSF from the brightest bumps in the image in every iteration. My guess is that this is to avoid overcorrecting.