regularizepsf is a Python package (with Cython speed improvements) for determining and correcting
point spread functions in astronomical images.
It was originally developed for the PUNCH mission and is documented in a forthcoming
Astrophysical Journal paper (see arXiv version until publication).
Below is an example of correcting model data using the package. An initial image of a simplified starfield (a) is synthetically observed with a slowly varying PSF (b), then regularized with this technique (c). The final image visually matches a direct convolution of the initial image with the target PSF (d). The panels are gamma-corrected to highlight the periphery of the model PSFs.