ja, so etwas kann man im prinzip "herausrechnen".
man geht davon aus, dass das beobachtete bild B aus dem urbild U und einer aufgefalteten störung S entstanden ist.
B = U x S
in deinem fall ist S die point spread funktion der abbildung oder besser die psf der gestörten abbildung. diese "unschärfe" wird nun auf jedes detail in U aufgefaltet und somit im ganzen bild B verbreitet.
um auf das urbild zu schliessen nutzt man die mathematische eigenheit, das die fouriertransformierte der faltung dem produkt der fouriertransformierten von U und S entspricht.
FT(B) = FT( U x S ) = FT(U) * FT(S)
diese FT's lassen sich für matrizen mittels gängiger fft techniken recht schnell berechnen. mit der gemessenen oder gut geschätzten funktion S und der inversion der FT matrizen kann man U errechnen (wenn ich keinen fehler gemacht habe [:)])
FT(B) * FT-1(S) = FT(U)
FT-1( FT(B) * FT-1(S) ) = U
gruss,
michael
ps: für weitere details kannst du dich direkt an mich wenden, das sollten wir nicht in formeltiefe hier austauschen...