Constrained least squares estimation is a technique for solution of integral equations of the first kind. The problem of image restoration requires the solution of an integral equation of the first kind. However, application of constrained least squares estimation to image restoration requires the solution of extremely large linear systems of equations. In this paper we demonstrate that, for convolution-type models of image restoration, special properties of the linear system of equations can be used to reduce the computational requirements. The necessary computations can be carried out by the fast Fourier transform, and the constrained least squares estimate can be constructed in the discrete frequency domain. A practical procedure for constrained least squares estimation is presented, and two examples are shown as output from a program for the CDC 7600 computer which performs the constrained least squares restoration of digital images.
[1]
David L. Phillips,et al.
A Technique for the Numerical Solution of Certain Integral Equations of the First Kind
,
1962,
JACM.
[2]
Thomas G. Stockham,et al.
High-speed convolution and correlation
,
1966,
AFIPS '66 (Spring).
[3]
B. Hunt.
A matrix theory proof of the discrete convolution theorem
,
1971
.
[4]
B. Hunt.
A theorem on the difficulty of numerical deconvolution
,
1972
.
[5]
M. M. Sondhi,et al.
Image restoration: The removal of spatially invariant degradations
,
1972
.
[6]
B. R. Hunt.
Data structures and computational organization in digital image enhancement
,
1972
.