LAUSR

Blind image deconvolution ($$\text {BID}$$BID) is one of the most important problems in image processing, and it requires the determination of an exact image $$\mathcal {F}$$F from a degraded form of it $$\mathcal {G}$$G when little or no information about $$\mathcal {F}$$F and the point spread function $$(\text {PSF}$$(PSF) $$\mathcal {H}$$H is known. Several methods have been developed for the solution of this problem, and one class of methods considers $$\mathcal {F}, \mathcal {G}$$F,G and $$\mathcal {H}$$H to be bivariate polynomials in which the polynomial computations are implemented by the Sylvester or Bézout resultant matrices. This paper compares these matrices for the solution of the problem of $$\text {BID}$$BID, and it is shown that it reduces to a comparison of their effectiveness for greatest common divisor ($$\text {GCD}$$GCD) computations. This is a difficult problem because the determination of the degree of the $$\text {GCD}$$GCD of two polynomials requires the calculation of the rank of a matrix, and this rank determines the size of the $$\text {PSF}$$PSF. It is shown that although the Bézout matrix is symmetric (unlike the Sylvester matrix) and it is smaller than the Sylvester matrix, which has computational advantages, it yields consistently worse results than the Sylvester matrix for the size and coefficients of the $$\text {PSF}$$PSF. Computational examples of blurred and deblurred images obtained with the Sylvester and Bézout matrices are shown, and the superior results obtained with the Sylvester matrix are evident.