Iterative algorithms for deblurring and deconvolution with constraints

Abstract

The deblurring problem is that of recovering the function from (noisy) values . The discrete finite version of the problem is to solve the system of linear equations Hc=d for c, where H is a matrix and d and c are vectors. When the kernel is a function of the difference , the deblurring problem becomes a deconvolution problem. The use of iterative algorithms to effect deblurring subject to non-negativity constraints on c has been presented by Snyder et al for the case of non-negative kernel function h. In this paper we extend these algorithms to include upper and lower bounds on the entries of the desired solution. We show that any linear deblurring problem involving a real kernel h can be transformed into a linear deblurring problem involving a non-negative kernel. Therefore our algorithms apply to general deblurring and deconvolution problems. These algorithms converge to a solution of the system of equations y = Px, with , for , , satisfying the vector inequalities , whenever such a solution exists. When there is no solution satisfying the constraints the simultaneous versions converge to an approximate solution that minimizes a cost function related to the Kullback-Leibler cross-entropy and the Fermi-Dirac generalized entropy.