Abstract

This paper outlines a method for the determination of the unknown wave-front aberration function of an optical system from noisy measurements of the corresponding point-spread function. The problem is cast as a nonlinear least-squares estimation problem for the values of the wave-front aberration function at N points over the slit aperture, from measurements of the point-spread function at M points with $M⩾N$. Newton’s method is used to replace the nonlinear minimization problem with a sequence of linear problems. Each such problem requires the inversion of the Hessian matrix of the error metric that is shown to be both singular (with rank $⩽N-1$) and ill-conditioned. To overcome singularity, the pseudoinverse is used; to overcome ill-conditioning the pseudoinverse is calculated using singular value decomposition and the singular values then filtered. Attention is drawn to difficulties such as nonuniqueness, sensitivity of algorithms to initial guess, etc.; the ancillary mathematical details being set out in appendices. Some illustrative numerical results are presented and analyzed.

© 1980 Optical Society of America

Determination of the wave-front aberration function from measured values of the point-spread function: a two-dimensional phase retrieval problem

Richard Barakat and Barbara H. Sandler
J. Opt. Soc. Am. A 9(10) 1715-1723 (1992)

Wave-front estimation from wave-front slope measurements

W.H. Southwell
J. Opt. Soc. Am. 70(8) 998-1006 (1980)

Optimum balanced wave-front aberrations for radially symmetric amplitude distributions: Generalizations of Zernike polynomials

Richard Barakat
J. Opt. Soc. Am. 70(6) 739-742 (1980)

Retrieval of wave aberration from point spread function or optical transfer function data

Junji Maeda and Kazumi Murata
Appl. Opt. 20(2) 274-279 (1981)

Retrieval of wave aberration of human eyes from actual point-spread-function data

Pablo Artal, Javier Santamaría, and Julian Bescós
J. Opt. Soc. Am. A 5(8) 1201-1206 (1988)