BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Ronny Ramlau

Ronny Ramlau contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisastro-ph.IMmath.SP

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2022arXiv

LBT SOUL data as a science test bench for MICADO PSF-R tool

Current state-of-the-art adaptive optics (AO) provides ground-based, diffraction-limited observations with high Strehl ratios (SR). However, a detailed knowledge of the point spread function (PSF) is required to fully exploit the scientific potential of these data. This is even more crucial for the next generation AO instruments that will equip 30-meter class telescopes, as the characterization of the PSF will be mandatory to fulfill the planned scientific requirements. For this reason, there is a growing interest in developing tools that accurately reconstruct the observed PSF of AO systems, the so-called PSF reconstruction. In this context, a PSF-R service is a planned deliverable for the MICADO@ELT instrument and our group is in charge of its development. In the case of MICADO, a blind PSF-R approach is being pursued to have the widest applicability to science cases. This means that the PSF is reconstructed without extracting information from the science data, relying only on telemetry and calibrations. While our PSF-R algorithm is currently being developed, its implementation is mature enough to test performances with actual observations. In this presentation we will discuss the reliability of our reconstructed PSFs and the uncertainties introduced in the measurements of scientific quantities for bright, on-axis observations taken with the SOUL+LUCI instrument of the LBT. This is the first application of our algorithm to real data. It demonstrates its readiness level and paves the way to further testing. Our PSF-R algorithm is able to reconstruct the SR and full-width at half maximum of the observed PSFs with errors smaller than 2% and 4.5%, respectively. We carried out the scientific evaluation of the obtained reconstructed PSFs thanks to a dedicated set of simulated observations of an ideal science case.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Performance of an iterative wavelet reconstructor for the Multi-conjugate Adaptive Optics RelaY of ESO's ELT

The Multi-conjugate Adaptive Optics RelaY (MAORY) is one of the key Adaptive Optics (AO) systems on the European Southern Observatory's Extremely Large Telescope. MAORY aims to achieve a good wavefront correction over a large field of view, which involves a tomographic estimation of the 3D atmospheric wavefront disturbance. Mathematically, the reconstruction of turbulent layers in the atmosphere is severely ill-posed, hence, limits the achievable reconstruction accuracy. Moreover, the reconstruction has to be performed in real-time at a few hundred to one thousand Hertz frame rates. Huge amounts of data have to be processed and thousands of actuators of the deformable mirrors have to be controlled by elaborated algorithms. Even with extensive parallelization and pipelining, direct solvers, such as the Matrix Vector Multiplication (MVM) method, are extremely demanding. Thus, research in the last years shifted into the direction of iterative methods. In this paper we focus on the iterative Finite Element Wavelet Hybrid Algorithm (FEWHA). The key feature of FEWHA is a matrix-free representation of all operators involved, which makes the algorithm fast and enables on the fly system updates whenever parameters at the telescope or in the atmosphere change. We provide a performance analysis of the method regarding quality and run-time for the MAORY instrument using the AO software package COMPASS.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Conjugate-Gradient Approach to the Parameter Estimation Problem of Magnetic Resonance Advection Imaging

We consider the inverse problem of estimating the spatially varying pulse wave velocity in blood vessels in the brain from dynamic MRI data, as it appears in the recently proposed imaging technique of Magnetic Resonance Advection Imaging (MRAI). The problem is formulated as a linear operator equation with a noisy operator and solved using a conjugate gradient type approach. Numerical examples on experimental data show the usefulness and advantages of the developed algorithm in comparison to previously proposed methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A note on the minimization of a Tikhonov functional with $\ell^1$-penalty

In this paper, we consider the minimization of a Tikhonov functional with an $\ell_1$ penalty for solving linear inverse problems with sparsity constraints. One of the many approaches used to solve this problem uses the Nemskii operator to transform the Tikhonov functional into one with an $\ell_2$ penalty term but a nonlinear operator. The transformed problem can then be analyzed and minimized using standard methods. However, by the nature of this transform, the resulting functional is only once continuously differentiable, which prohibits the use of second order methods. Hence, in this paper, we propose a different transformation, which leads to a twice differentiable functional that can now be minimized using efficient second order methods like Newton's method. We provide a convergence analysis of our proposed scheme, as well as a number of numerical results showing the usefulness of our proposed approach.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Estimating solution smoothness and data noise with Tikhonov regularization

A main drawback of classical Tikhonov regularization is that often the parameters required to apply theoretical results, e.g., the smoothness of the sought-after solution and the noise level, are unknown in practice. In this paper we investigate in new detail the residuals in Tikhonov regularization viewed as functions of the regularization parameter. We show that the residual carries, with some restrictions, the information on both the unknown solution and the noise level. By calculating approximate solutions for a large range of regularization parameters, we can extract both parameters from the residual given only one set of noisy data and the forward operator. The smoothness in the residual allows to revisit parameter choice rules and relate a-priori, a-posteriori, and heuristic rules in a novel way that blurs the lines between the classical division of the parameter choice rules. All results are accompanied by numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Real-time implementation of an iterative solver for atmospheric tomography

The image quality of the new generation of earthbound Extremely Large Telescopes (ELTs) is heavily influenced by atmospheric turbulences. To compensate these optical distortions a technique called adaptive optics (AO) is used. Many AO systems require the reconstruction of the refractive index fluctuations in the atmosphere, called atmospheric tomography. The standard way of solving this problem is the Matrix Vector Multiplication, i.e., the direct application of a (regularized) generalized inverse of the system operator. However, over the last years the telescope sizes have increased significantly and the computational efficiency become an issue. Promising alternatives are iterative methods such as the Finite Element Wavelet Hybrid Algorithm (FEWHA), which is based on wavelets. Due to its efficient matrix-free representation of the underlying operators, the number of floating point operations and memory usage decreases significantly. In this paper, we focus on performance optimization techniques, such as parallel programming models, for the implementation of this iterative method on CPUs and GPUs. We evaluate the computational performance of our optimized, parallel version of FEWHA for ELT-sized test configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Nesterov's Accelerated Gradient Method for Nonlinear Ill-Posed Problems with a Locally Convex Residual Functional

In this paper, we consider Nesterov's Accelerated Gradient method for solving Nonlinear Inverse and Ill-Posed Problems. Known to be a fast gradient-based iterative method for solving well-posed convex optimization problems, this method also leads to promising results for ill-posed problems. Here, we provide a convergence analysis for ill-posed problems of this method based on the assumption of a locally convex residual functional. Furthermore, we demonstrate the usefulness of the method on a number of numerical examples based on a nonlinear diagonal operator and on an inverse problem in auto-convolution.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

On the singular value decomposition of n-fold integration operators

In theory and practice of inverse problems, linear operator equations $Tx=y$ with compact linear forward operators $T$ having a non-closed range $\mathcal{R}(T)$ and mapping between infinite dimensional Hilbert spaces plays some prominent role. As a consequence of the ill-posedness of such problems, regularization approaches are required, and due to its unlimited qualification spectral cut-off is an appropriate method for the stable approximate solution of corresponding inverse problems. For this method, however, the singular system $\{σ_i(T),u_i(T),v_i(T)\}_{i=1}^\infty$ of the compact operator $T$ is needed, at least for $i=1,2,...,N$, up to some stopping index $N$. In this note we consider $n$-fold integration operators $T=J^n\;(n=1,2,...)$ in $L^2([0,1])$ occurring in numerous applications, where the solution of the associated operator equation is characterized by the $n$-th generalized derivative $x=y^{(n)}$ of the Sobolev space function $y \in H^n([0,1])$. Almost all textbooks on linear inverse problems present the whole singular system $\{σ_i(J^1),u_i(J^1),v_i(J^1)\}_{i=1}^\infty$ in an explicit manner. However, they do not discuss the singular systems for $J^n,\;n \ge 2$. We will emphasize that this seems to be a consequence of the fact that for higher $n$ the eigenvalues $σ^2_i(J^n)$ of the associated ODE boundary value problems obey transcendental equations, the complexity of which is growing with $n$. We present the transcendental equations for $n=2,3,...$ and discuss and illustrate the associated eigenfunctions and some of their properties.

0 citations0 reviewsNo code linkedOpen access