Not directly involved with this... I have a working TGVRestoration module for deconvolutions, that may be used to perform this sort of things with a minimum number of changes. Basically, all I need is a way to model the aberrations (field curvature, astigmatism, chromatic aberrations, etc.), in a global way, hopefully fast. Then, I need to calculate de adjoint operator of that degradation function. With both operations, we are set to feed the iterative algorithm, and restore the image.
If somebody gives me a hand with that, things may speed up. Right now I'm quite busy with another image processing problem, for my Ph. D. degree (magnetic susceptibility maps of brains with MRI).
Sorry for the slow advance in this.