A FREQUENCY DEPENDENT PRECONDITIONED WAVELET METHOD FOR ATMOSPHERIC TOMOGRAPHY

Atmospheric tomography, i.e. the reconstruction of the turbulence in the atmosphere, is a main task for the adaptive optics systems of the next generation telescopes. For extremely large telescopes, such as the European Extremely Large Telescope, this problem becomes overly complex and an efficient algorithm is needed to reduce numerical costs. Recently, a conjugate gradient method based on wavelet parametrization of turbulence layers was introduced [5]. An iterative algorithm can only be numerically efficient when the number of iterations required for a sufficient reconstruction is low. A way to achieve this is to design an efficient preconditioner. In this paper we propose a new frequency– dependent preconditioner for the wavelet method. In the context of a multi conjugate adaptive optics (MCAO) system simulated on the official end-to-end simulation tool OCTOPUS of the European Southern Observatory we demonstrate robustness and speed of the preconditioned algorithm. We show that three iterations are sufficient for a good reconstruction.

[1]  Steven A. Orszag,et al.  CBMS-NSF REGIONAL CONFERENCE SERIES IN APPLIED MATHEMATICS , 1978 .

[2]  Ingrid Daubechies,et al.  Ten Lectures on Wavelets , 1992 .

[3]  B. Welsh,et al.  Imaging Through Turbulence , 1996 .

[4]  L M Mugnier,et al.  Optimal wave-front reconstruction strategies for multiconjugate adaptive optics. , 2001, Journal of the Optical Society of America. A, Optics, image science, and vision.

[5]  Z. Strakos,et al.  On error estimation in the conjugate gradient method and why it works in finite precision computations. , 2002 .

[6]  Roberto Ragazzoni,et al.  MAD the ESO multi-conjugate adaptive optics demonstrator , 2003, SPIE Astronomical Telescopes + Instrumentation.

[7]  Brent L. Ellerbroek,et al.  Simulations of closed-loop wavefront reconstruction for multiconjugate adaptive optics on giant telescopes , 2003, SPIE Optics + Photonics.

[8]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[9]  E. Somersalo,et al.  Statistical and computational inverse problems , 2004 .

[10]  Luc Gilles,et al.  Robustness study of the pseudo open-loop controller for multiconjugate adaptive optics. , 2005, Applied optics.

[11]  Qiang Yang,et al.  Fourier Domain Preconditioned Conjugate Gradient Algorithm for Atmospheric Tomography , 2005 .

[12]  Curtis R Vogel,et al.  Fast optimal wavefront reconstruction for multi-conjugate adaptive optics using the Fourier domain preconditioned conjugate gradient algorithm. , 2006, Optics express.

[13]  Norbert Hubin,et al.  Adaptive optics simulations for the European Extremely Large Telescope , 2006, SPIE Astronomical Telescopes + Instrumentation.

[14]  Richard M. Clare,et al.  Closed-loop ground layer adaptive optics simulations with elongated spots : impact of modeling noise correlations , 2010 .

[15]  Michel Tallon,et al.  Fast minimum variance wavefront reconstruction for extremely large telescopes. , 2010, Journal of the Optical Society of America. A, Optics, image science, and vision.

[16]  Fabien Momey,et al.  Fractal iterative method for fast atmospheric tomography on extremely large telescopes , 2010, Astronomical Telescopes + Instrumentation.

[17]  Ronny Ramlau,et al.  An efficient solution to the atmospheric turbulence tomography problem using Kaczmarz iteration , 2012 .

[18]  Matthias Rosensteiner,et al.  Wavefront reconstruction for extremely large telescopes via CuRe with domain decomposition. , 2012, Journal of the Optical Society of America. A, Optics, image science, and vision.

[19]  Tapio Helin,et al.  Wavelet methods in multi-conjugate adaptive optics , 2013, 1302.3734.