The Sylvester and Bézout resultant matrices for blind image deconvolution

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

Polynomial computations for blind image deconvolution

This paper considers the problem of blind image deconvolution (BID) when the blur arises from a spatially invariant point spread function (PSF) H, which implies that a blurred image G is formed by the convolution of H and the exact form F of G. Since the multiplication of two bivariate polynomials is performed by convolving their coefficient matrices, the equivalence of the formation of a blurred image and the product of two bivariate polynomials implies that BID can be performed by considering F, G and H to be bivariate polynomials on which polynomial operations are performed. These operations allow the PSF to be computed, which is then deconvolved from the blurred image G, thereby obtaining a deblurred image that is a good approximation of the exact image F. Computational results show that the deblurred image obtained using polynomial computations is better than the deblurred image obtained using other methods for blind image deconvolution