BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Tapio Helin

Tapio Helin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.APmath.PRastro-ph.IMmath.NAmath.STStatistics TheoryMachine LearningRobotics

Catalog footprint

What is connected

11works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

11 published item(s)

preprint2026arXiv

Recovering a (1+1)-dimensional wave equation from a single white noise boundary measurement

We consider the following inverse problem: Suppose a $(1+1)$-dimensional wave equation on $\mathbb{R}_+$ with zero initial conditions is excited with a Neumann boundary data modelled as a white noise process. Given also the Dirichlet data at the same point, determine the unknown first order coefficient function of the system. We first establish that direct problem is well-posed. The inverse problem is then solved by showing that correlations of the boundary data determine the Neumann-to-Dirichlet operator in the sense of distributions, which is known to uniquely identify the coefficient. This approach has applications in acoustic measurements of internal cross-sections of fluid pipes such as pressurised water supply pipes and vocal tract shape determination.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Laboratory Experiments of Model-based Reinforcement Learning for Adaptive Optics Control

Direct imaging of Earth-like exoplanets is one of the most prominent scientific drivers of the next generation of ground-based telescopes. Typically, Earth-like exoplanets are located at small angular separations from their host stars, making their detection difficult. Consequently, the adaptive optics (AO) system's control algorithm must be carefully designed to distinguish the exoplanet from the residual light produced by the host star. A new promising avenue of research to improve AO control builds on data-driven control methods such as Reinforcement Learning (RL). RL is an active branch of the machine learning research field, where control of a system is learned through interaction with the environment. Thus, RL can be seen as an automated approach to AO control, where its usage is entirely a turnkey operation. In particular, model-based reinforcement learning (MBRL) has been shown to cope with both temporal and misregistration errors. Similarly, it has been demonstrated to adapt to non-linear wavefront sensing while being efficient in training and execution. In this work, we implement and adapt an RL method called Policy Optimization for AO (PO4AO) to the GHOST test bench at ESO headquarters, where we demonstrate a strong performance of the method in a laboratory environment. Our implementation allows the training to be performed parallel to inference, which is crucial for on-sky operation. In particular, we study the predictive and self-calibrating aspects of the method. The new implementation on GHOST running PyTorch introduces only around 700 microseconds in addition to hardware, pipeline, and Python interface latency. We open-source well-documented code for the implementation and specify the requirements for the RTC pipeline. We also discuss the important hyperparameters of the method, the source of the latency, and the possible paths for a lower latency implementation.

0 citations0 reviewsNo code linkedOpen access
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
preprint2016arXiv

Correlation based passive imaging with a white noise source

Passive imaging refers to problems where waves generated by unknown sources are recorded and used to image the medium through which they travel. The sources are typically modelled as a random variable and it is assumed that some statistical information is available. In this paper we study the stochastic wave equation $\partial_t^2 u - Δ_g u = χW$, where $W$ is a random variable with the white noise statistics on ${\mathbb R}^{1+n}$, $n \ge 3$, $χ$ is a smooth function vanishing for negative times and outside a compact set in space, and $Δ_g$ is the Laplace-Beltrami operator associated to a smooth non-trapping Riemannian metric tensor $g$ on ${\mathbb R}^n$. The metric tensor $g$ models the medium to be imaged, and we assume that it coincides with the Euclidean metric outside a compact set. We consider the empirical correlations on an open set $\mathcal X \subset {\mathbb R}^n$, $$ C_T(t_1, x_1, t_2, x_2) = \frac 1 T \int_0^T u(t_1+s,x_1) u(t_2+s,x_2) ds, \quad t_1,t_2>0,\ x_1,x_2\in \mathcal X, $$ for $T>0$. Supposing that $χ$ is non-zero on $\mathcal X$ and constant in time after $t > 1$, we show that in the limit $T \to \infty$, the data $C_T$ becomes statistically stable, that is, independent of the realization of $W$. Our main result is that, with probability one, this limit determines the Riemannian manifold $({\mathbb R}^n,g)$ up to an isometry. To our knowledge, this is the first result showing that a medium can be determined in a passive imaging setting, without assuming a separation of scales.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Inverse scattering for a random potential

In this paper we consider an inverse problem for the $n$-dimensional random Schrödinger equation $(Δ-q+k^2)u = 0$. We study the scattering of plane waves in the presence of a potential $q$ which is assumed to be a Gaussian random function such that its covariance is described by a pseudodifferential operator. Our main result is as follows: given the backscattered far field, obtained from a single realization of the random potential $q$, we uniquely determine the principal symbol of the covariance operator of $q$. Especially, for $n=3$ this result is obtained for the full non-linear inverse backscattering problem. Finally, we present a physical scaling regime where the method is of practical importance.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Maximum a posteriori probability estimates in infinite-dimensional Bayesian inverse problems

A demanding challenge in Bayesian inversion is to efficiently characterize the posterior distribution. This task is problematic especially in high-dimensional non-Gaussian problems, where the structure of the posterior can be very chaotic and difficult to analyse. Current inverse problem literature often approaches the problem by considering suitable point estimators for the task. Typically the choice is made between the maximum a posteriori (MAP) or the conditional mean (CM) estimate. The benefits of either choice are not well-understood from the perspective of infinite-dimensional theory. Most importantly, there exists no general scheme regarding how to connect the topological description of a MAP estimate to a variational problem. The results by Dashti et. al. (2013) resolve this issue for non-linear inverse problems in Gaussian framework. In this work we improve the current understanding by introducing a novel concept called the weak MAP (wMAP) estimate. We show that any MAP estimate in the sense of Dashti et. al. (2013) is a wMAP estimate and, moreover, how in general infinite-dimensional non-Gaussian problems the wMAP estimate connects to a variational formulation. Such a formulation yields many properties of the estimate that were earlier impossible to study. In a recent work by Burger and Lucka (2014) the MAP estimator was studied in the context of Bayes cost method. Using Bregman distances, proper convex Bayes cost functions were introduced for which the MAP estimator is the Bayes estimator. Here, we generalize these results to the infinite-dimensional setting. Moreover, we discuss the implications of our results for some examples of prior models such as the Besov prior and hierarchical prior.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Towards analytical model optimization in atmospheric tomography

Modern ground-based telescopes rely on a technology called adaptive optics (AO) in order to compensate for the loss of image quality caused by atmospheric turbulence. Next-generation AO systems designed for a wide field of view require a stable and high-resolution reconstruction of the refractive index fluctuations in the atmosphere. By introducing a novel Bayesian method, we address the problem of estimating an atmospheric turbulence strength profile and reconstructing the refractive index fluctuations simultaneously, where we only use wavefront measurements of incoming light from guide stars. Most importantly, we demonstrate how this method can be used for model optimization as well. We propose two different algorithms for solving the maximum a posteriori estimate: the first approach is based on alternating minimization and has the advantage of integrability into existing atmospheric tomography methods. In the second approach, we formulate a convex non-differentiable optimization problem, which is solved by an iterative thresholding method. This approach clearly illustrates the underlying sparsity-enforcing mechanism for the strength profile. By introducing a tuning/regularization parameter, an automated model reduction of the layer structure of the atmosphere is achieved. Using numerical simulations, we demonstrate the performance of our method in practice.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Inverse acoustic scattering problem in half-space with anisotropic random impedance

We study an inverse acoustic scattering problem in half-space with a probabilistic impedance boundary value condition. The Robin coefficient (surface impedance) is assumed to be a Gaussian random function $λ= λ(x)$ with a pseudodifferential operator describing the covariance. We measure the amplitude of the backscattered field averaged over the frequency band and assume that the data is generated by a single realization of $λ$. Our main result is to show that under certain conditions the principal symbol of the covariance operator of $λ$ is uniquely determined. Most importantly, no approximations are needed and we can solve the full non-linear inverse problem. We concentrate on anisotropic models for the principal symbol, which leads to the analysis of a novel anisotropic spherical Radon transform and its invertibility.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Inverse problem for the wave equation with a white noise source

We consider a smooth Riemannian metric tensor $g$ on $\R^n$ and study the stochastic wave equation for the Laplace-Beltrami operator $\p_t^2 u - Δ_g u = F$. Here, $F=F(t,x,ω)$ is a random source that has white noise distribution supported on the boundary of some smooth compact domain $M \subset \R^n$. We study the following formally posed inverse problem with only one measurement. Suppose that $g$ is known only outside of a compact subset of $M^{int}$ and that a solution $u(t,x,ω_0)$ is produced by a single realization of the source $F(t,x,ω_0)$. We ask what information regarding $g$ can be recovered by measuring $u(t,x,ω_0)$ on $\R_+ \times \p M$? We prove that such measurement together with the realization of the source determine the scattering relation of the Riemannian manifold $(M, g)$ with probability one. That is, for all geodesics passing through $M$, the travel times together with the entering and exit points and directions are determined. In particular, if $(M,g)$ is a simple Riemannian manifold and $g$ is conformally Euclidian in $M$, the measurement determines the metric $g$ in $M$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Wavelet methods in multi-conjugate adaptive optics

The next generation ground-based telescopes rely heavily on adaptive optics for overcoming the limitation of atmospheric turbulence. In the future adaptive optics modalities, like multi-conjugate adaptive optics (MCAO), atmospheric tomography is the major mathematical and computational challenge. In this severely ill-posed problem a fast and stable reconstruction algorithm is needed that can take into account many real-life phenomena of telescope imaging. We introduce a novel reconstruction method for the atmospheric tomography problem and demonstrate its performance and flexibility in the context of MCAO. Our method is based on using locality properties of compactly supported wavelets, both in the spatial and frequency domain. The reconstruction in the atmospheric tomography problem is obtained by solving the Bayesian MAP estimator with a conjugate gradient based algorithm. An accelerated algorithm with preconditioning is also introduced. Numerical performance is demonstrated on the official end-to-end simulation tool OCTOPUS of European Southern Observatory.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An inverse problem for the wave equation with one measurement and the pseudorandom noise

We consider the wave equation $(\p_t^2-Δ_g)u(t,x)=f(t,x)$, in $\R^n$, $u|_{\R_-\times \R^n}=0$, where the metric $g=(g_{jk}(x))_{j,k=1}^n$ is known outside an open and bounded set $M\subset \R^n$ with smooth boundary $\p M$. We define a deterministic source $f(t,x)$ called the pseudorandom noise as a sum of point sources, $f(t,x)=\sum_{j=1}^\infty a_jδ_{x_j}(x)δ(t)$, where the points $x_j,\ j\in\Z_+$, form a dense set on $\p M$. We show that when the weights $a_j$ are chosen appropriately, $u|_{\R\times \p M}$ determines the scattering relation on $\p M$, that is, it determines for all geodesics which pass through $M$ the travel times together with the entering and exit points and directions. The wave $u(t,x)$ contains the singularities produced by all point sources, but when $a_j=λ^{-λ^{j}}$ for some $λ>1$, we can trace back the point source that produced a given singularity in the data. This gives us the distance in $(\R^n, g)$ between a source point $x_j$ and an arbitrary point $y \in \p M$. In particular, if $(\bar M,g)$ is a simple Riemannian manifold and $g$ is conformally Euclidian in $\bar M$, these distances are known to determine the metric $g$ in $M$. In the case when $(\bar M,g)$ is non-simple we present a more detailed analysis of the wave fronts yielding the scattering relation on $\p M$.

0 citations0 reviewsNo code linkedOpen access