Abstract

Wavefronts reconstructed from measured gradients are composed of a straightforward integration of the measured data, plus a correction term that disappears when there are no measurement errors. For regions of any shape, this term is a solution of Poisson’s equation with Dirichlet conditions ($V=0$ on the boundaries). We show that for rectangular regions, the correct solution is not a periodic one, but one expressed with Fourier cosine series. The correct solution has a lower variance than the periodic Fourier transform solution. Similar formulas exist for a circular region with obscuration. We present a near-optimal solution that is much faster than fast-Fourier-transform methods. By use of diagonal multigrid methods, a single iteration brings the correction term to within a standard deviation of 0.08, two iterations, to within 0.0064, etc.

© 2006 Optical Society of America

Modal wavefront reconstruction from its gradient

Iacopo Mochi and Kenneth A. Goldberg
Appl. Opt. 54(12) 3780-3785 (2015)

Jacobi circle and annular polynomials: modal wavefront reconstruction from wavefront gradient

Wenhan Sun, Shuai Wang, Xing He, and Bing Xu
J. Opt. Soc. Am. A 35(7) 1140-1148 (2018)

Noniterative boundary-artifact-free wavefront reconstruction from its derivatives

Pierre Bon, Serge Monneret, and Benoit Wattellier
Appl. Opt. 51(23) 5698-5704 (2012)

Multigrid preconditioned conjugate-gradient method for large-scale wave-front reconstruction

Luc Gilles, Curtis R. Vogel, and Brent L. Ellerbroek
J. Opt. Soc. Am. A 19(9) 1817-1822 (2002)

Multigrid algorithm for least-squares wavefront reconstruction

C. R. Vogel and Q. Yang
Appl. Opt. 45(4) 705-715 (2006)