On some analytic properties of the atmospheric tomography operator: Non-Uniqueness and reconstructability issues
Ronny Ramlau, Bernadett Stadler
TL;DR
The paper addresses the atmospheric tomography problem in adaptive optics, showing that data from a finite set of guide stars under a limited-angle geometry cannot uniquely determine the turbulent profile above the telescope, formalized via the operator $\mathbf{A}$. It develops a geometric overlap analysis using the overlap function $\omega_l$ and proves non-uniqueness by identifying regions with unit overlap that yield nontrivial nullspaces for $G>1$. The authors also investigate classical reconstruction methods, including $Tikhonov$ regularization, $Landweber$ iteration, and the $Kaczmarz$ method, demonstrating that they generally fail to reconstruct physically meaningful turbulence distributions in this setting. Numerical simulations corroborate the theory, highlighting fundamental ill-posedness and guiding the need for alternative modeling or regularization approaches in MCAO/LTAO/MOAO contexts.
Abstract
In this paper, we consider the atmospheric tomography operator, which describes the effect of turbulent atmospheric layers on incoming planar wavefronts. Given wavefronts from different guide stars, measured at a telescope, the inverse problem consists in the reconstruction of the turbulence above the telescope. We show that the available data is not sufficient to reconstruct the atmosphere uniquely. Additionally, we show that classical regularization methods as Tikhonov regularization or Landweber iteration will always fail to reconstruct a physically meaningful turbulence distribution.
