OSA's Digital Library

Applied Optics

Applied Optics


  • Editor: Joseph N. Mait
  • Vol. 50, Iss. 21 — Jul. 20, 2011
  • pp: 4124–4133

Multiscale stochastic approach for phase screens synthesis

Alessandro Beghi, Angelo Cenedese, and Andrea Masiero  »View Author Affiliations

Applied Optics, Vol. 50, Issue 21, pp. 4124-4133 (2011)

View Full Text Article

Enhanced HTML    Acrobat PDF (439 KB)

Browse Journals / Lookup Meetings

Browse by Journal and Year


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools



Simulating the turbulence effect on ground telescope observations is of fundamental importance for the design and test of suitable control algorithms for adaptive optics systems. In this paper we propose a multiscale approach for efficiently synthesizing turbulent phases at very high resolution. First, the turbulence is simulated at low resolution, taking advantage of a previously developed method for generating phase screens [ J. Opt. Soc. Am. A 25, 515 (2008)]. Then, high-resolution phase screens are obtained as the output of a multiscale linear stochastic system. The multiscale approach significantly improves the computational efficiency of turbulence simulation with respect to recently developed methods [ Opt. Express 14, 988 (2006)] [ J. Opt. Soc. Am. A 25, 515 (2008)] [ J. Opt. Soc. Am. A 25, 463 (2008)]. Furthermore, the proposed procedure ensures good accuracy in reproducing the statistical characteristics of the turbulent phase.

© 2011 Optical Society of America

OCIS Codes
(010.1330) Atmospheric and oceanic optics : Atmospheric turbulence
(350.5030) Other areas of optics : Phase

ToC Category:
Atmospheric and Oceanic Optics: Wave-front Sensing

Original Manuscript: February 1, 2011
Revised Manuscript: June 5, 2011
Manuscript Accepted: June 6, 2011
Published: July 19, 2011

Alessandro Beghi, Angelo Cenedese, and Andrea Masiero, "Multiscale stochastic approach for phase screens synthesis," Appl. Opt. 50, 4124-4133 (2011)

Sort:  Author  |  Year  |  Journal  |  Reset  


  1. F. Roddier, Adaptive Optics in Astronomy (Cambridge University, 1999). [CrossRef]
  2. F. Assemat, R. Wilson, and E. Gendron, “Method for simulating infinitely long and non stationary phase screens with optimized memory storage,” Opt. Express 14, 988–999 (2006). [CrossRef] [PubMed]
  3. A. Beghi, A. Cenedese, and A. Masiero, “Stochastic realization approach to the efficient simulation of phase screens,” J. Opt. Soc. Am. A 25, 515–525 (2008). [CrossRef]
  4. D. Fried and T. Clark, “Extruding Kolmogorov-type phase screen ribbons,” J. Opt. Soc. Am. A 25, 463–468 (2008). [CrossRef]
  5. P. S. Maybeck, Stochastic Models, Estimation, and Control (Academic, 1979), Vol.  1.
  6. R. G. Lane and A. Glindemann and J. C. Dainty, “Simulation of a Kolmogorov phase screen,” Waves Random Media 2, 209–224 (1992). [CrossRef]
  7. E. Thiebaut and M. Tallon, “Fast minimum variance wavefront reconstruction for extremely large telescope,” J. Opt. Soc. Am. A 27, 1046–1059 (2010). [CrossRef]
  8. M. Luettgen, W. Karl, A. Willsky, and R. Tenney, “Multiscale representations of Markov random fields,” IEEE Trans. Signal Process. 41, 3377–3396 (1993). [CrossRef]
  9. A. Benveniste, R. Nikoukhah, and A. Willsky, “Multiscale system theory,” IEEE Trans. Circuits Syst. I 41, 2–15 (1994). [CrossRef]
  10. K. Daoudi, A. Frakt, and A. Willsky, “Multiscale autoregressive models and wavelets,” IEEE Trans. Inf. Theory 45, 828–845 (1999). [CrossRef]
  11. W. Irving and A. Willsky, “A canonical correlations approach to multiscale stochastic realization,” IEEE Trans. Autom. Control 46, 1514–1528 (2001). [CrossRef]
  12. A. Frakt and A. Willsky, “Computationally efficient stochastic realization for internal multiscale autoregressive models,” Multidimens. Syst. Signal Process. 12, 109–142 (2001). [CrossRef]
  13. I. Daubechies, Ten Lectures on Wavelets, CBMS-NSF Regional Conference Series in Applied Mathematics 61 (SIAM, 1992). [CrossRef]
  14. S. Mallat, A Wavelet Tour of Signal Processing (Academic, 1999).
  15. R. Conan, “Modelisation des effets de l’echelle externe de coherence spatiale du front d’onde pour l’observation a haute resolution angulaire en astronomie,” Ph.D. thesis (Université Nice Sophia Antipolis, 2000).
  16. F. Roddier, “The effects of atmospheric turbulence in optical astronomy,” Prog. Opt. 19, 281–376 (1981). [CrossRef]
  17. From Eq. , wi(u) and wi′(u) are practically uncorrelated for large |u−u′|. Furthermore, intuitively, the conditioned value of wi(u), given the local values of xi (i.e., the process at lower resolution), is much more uncorrelated with wi(u′), u≠u′, with respect to wi(u).
  18. J. Woods, “Two-dimensional discrete Markovian fields,” IEEE Trans. Inf. Theory 18, 232–240 (1972). [CrossRef]
  19. M. B. Priestley, Spectral Analysis and Time Series (Academic, 1982), Vol.  1.
  20. A. Steinhardt, “Correlation matching by finite length sequences,” IEEE Trans. Acoust. Speech Signal Process. 36, 545–559 (1988). [CrossRef]
  21. M. Le Ravalec, B. Noetinger, and L. Hu, “The FFT moving average (FFT-MA) generator: an efficient numerical method for generating and conditioning Gaussian simulations,” Math. Geol. 32, 701–723 (2000). [CrossRef]
  22. A. Oppenheim and R. Schafer, Digital Signal Processing(Prentice-Hall, 1975).
  23. Imposing F(θi)(h)=F(ce,i)(h), the solution is restricted to be symmetric. However, there might exist a shorter nonsymmetric sequence of coefficients leading to the same covariances {ce,i(−di),…,ce,i(di)}.
  24. To be more precise, let m¯ and n¯ be the dimensions of a single phase screen (if the wind velocity is parallel to u, then typically m¯=m and n¯≪n) and assume m0≪m and n0≪n; then the number of operations computed by the algorithm to generate an m¯×n¯ phase screen are approximately 64d¯2m¯n¯. When iteratively generating an m×n screen, with n≫n¯, dividing it in approximately n/n¯ phase screens, there is an extra computational load due to the extra computations to ensure the continuity between successive phase screens. Anyway, such extra computational load is usually a minor term in the overall complexity (lower than 10% in our simulations). Similar considerations can be repeated for the memory requirements.
  25. Wavelets different from the Haar transform have better multiscale prediction ability but worse complexity and memory requirements.

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited