The problem of image restoration of space variant blur is common and important. One of the most useful descriptions of this problem is in its algebraic form I=H*O, where O is the object represented as a column vector, I is the blur image represented as a column vector and H is the PSF matrix that represents the optical system. When inverting the problem to restore the geometric object from the blurred image and the known system matrix, restoration is limited in speed and quality by the system condition. Current optical design methods focus on image quality, therefore if additional image processing is needed the matrix condition is taken "as is". In this paper we would like to present a new optical approach which aims to improve the system condition by proper optical design. In this new method we use Singular Value Decomposition (SVD) to define the weak parts of the matrix condition. We design a second optical system based on those weak SVD parts and then we add the second system parallel to the first one. The original and second systems together work as an improved parallel optics system. Following that, we present a method for designing such a "parallel filter" for systems with a spread SVD pattern. Finally we present a study case in which by using our new method we improve a space variant image system with an initial condition number of 8.76e4, down to a condition number of 2.29e3. We use matrix inversion to simulate image restoration. Results show that the new parallel optics immunity to Additive White Gaussian Noise (AWGN) is much better then that of the original simple lens. Comparing the original and the parallel optics systems, the parallel optics system crosses the MSEIF=0 [db] limit in SNR value which is more than 50db lower then the SNR value in the case of the original simple lens. The new parallel optics system performance is also compared to another method based on the MTF approach.
[1]
J. Varah.
On the Numerical Solution of Ill-Conditioned Linear Systems with Applications to Ill-Posed Problems
,
1973
.
[2]
B. R. Hunt,et al.
Image restoration of space-variant blurs by sectioned methods
,
1978
.
[3]
H. Andrews,et al.
Singular value decompositions and digital image processing
,
1976
.
[4]
N George,et al.
Electronic imaging using a logarithmic asphere.
,
2001,
Optics letters.
[5]
Andrew R Harvey,et al.
Phase pupil functions for reduction of defocus and spherical aberrations.
,
2003,
Optics letters.
[6]
Xinmao Wang,et al.
Effect of small rank modification on the condition number of a matrix
,
2007,
Comput. Math. Appl..
[7]
Gonzalo Muyo,et al.
Circularly symmetric phase filters for control of primary third-order aberrations: coma and astigmatism.
,
2006,
Journal of the Optical Society of America. A, Optics, image science, and vision.
[8]
W. Cathey,et al.
Extended depth of field through wave-front coding.
,
1995,
Applied optics.
[9]
Alexander A. Sawchuk,et al.
Restoration of astigmatism and curvature of field
,
1975
.
[10]
J. Rolland,et al.
Model of wide-angle optical field propagation using scalar diffraction theory
,
2004
.
[11]
S Twomey.
Information content in remote sensing.
,
1974,
Applied optics.
[12]
W T Rhodes,et al.
Two-pupil synthesis of optical transfer functions.
,
1978,
Applied optics.
[13]
Thomas S. Huang,et al.
Inverse filtering for linear shift-variant imaging systems
,
1972
.