For the latest deconvolution results, click here.
Blind Deconvolution
Blind deconvolution is has a lofty goal: To recover both the deconvolved image and the PSF using just the blurred image and some constraints as input. Results from two algorithms are presented here: 1. the Iterative Blind Deconvolvution method of Ayers & Dainty (IBD Method) and 2. The BiggsLucy Iterative methods (BL Method). In both algorithms the required constraint information is quite minimal (i.e. nonnegativity and support). The PSF is not known (often just a plain or random pixel field is used as prior estimate).
In the Biaram program DeconIB I have implemented an IBD algorithm based on the Ayers & Dainty method: G. R. Ayers and J. C. Dainty, Iterative blind deconvolution method and its applications. Optics Letters, (1988); 13:547549. This method iterates between the Fourier and spatial domain imposing these constraints at each step to refine the solution. The results of applying this method to a 64 x 64 pixel synthetic image blurred in the computer by a PSF (a 9x9 pixel symmetrical geometric intensity progression) are presented here. The images are shown at twice actual size and are contrast stretched.
In the Biaram program DeconBMI I have implemented the algorithm (see: D. S. C. Biggs and M. Andrews, Asymmetric iterative blind deconvolution of multiframe images. Proc. SPIE, (1998); 346:133.
This method iterates between a number of imagerefineing cycles and a (usually different) number of PSF refining cycles at each major iteration to refine the solution. The 'cycles' are just iterations of a usualtype LucyRichardson method which uses the currest PSF estimate when doing image refinement and then switches the roles of image and PSF to use the current image estimate (as if it were a PSF) to refine the current PSF estimate The results of applying this method to a 64 x 64 pixel synthetic image blurred in the computer by a PSF (a 15x15 pixel asymmetrical PSF in the shape of a Greek 'Psi') are presented here. The images are shown at twice actual size and the PSF images are shown contrast stretched.
23.02.2008
Experiment 1: Ayers & Dainty Method
Experiment 2: BiggsLucy Method
Input Images: a) The blurred test image, b) The estimate of the PSF (a plain pixel field). For this experiment I also use the known support area of the PSF (i.e. 15x15 pixels). Knowledge of the maximum and minimum grey levels in the image solution was used for constraints but no information about the spatial support of the image is used. Ground truth image and PSF are not shown but are very similar to the results of the deconvolution shown below. 
a) b) 
Result: a) The estimate of the deconvolved image, b) The estimate of the PSF. These results were obtained after 1048576 major iterations. The image and PSF cycels were equally matched in this experiment at 1 cycle each per major iteration. 
a) b) 
Convergence Measure: The value for the 'residual' is calculated as the sum of the absolute differences from 1.0 of the ratio: (convolved estimate) / (input image) The result is divided by the number of pixels in the solution to give residual per pixel and is then multiplied by 100 to give the residual per 100 pixels. Iterations are shown on a log scale. Note the sudden 'bend' in the graph at about 100 iterations. This indicates a sudden 'snap' of the solution as it catches onto the 'right' path. It is clear that the solution velocity is complicated and not simply linear or logarithmic. 

