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

Lassi Roininen

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

ResearcherUnclaimed source record
math.STStatistics TheoryMethodologyApplicationsComputationComputational Engineering, Finance, and Scienceeess.IVeess.SPMachine Learningmath.PRphysics.comp-phphysics.geo-phq-fin.PR

Catalog footprint

What is connected

10works
13topics
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

10 published item(s)

preprint2022arXiv

Cauchy Markov Random Field Priors for Bayesian Inversion

The use of Cauchy Markov random field priors in statistical inverse problems can potentially lead to posterior distributions which are non-Gaussian, high-dimensional, multimodal and heavy-tailed. In order to use such priors successfully, sophisticated optimization and Markov chain Monte Carlo (MCMC) methods are usually required. In this paper, our focus is largely on reviewing recently developed Cauchy difference priors, while introducing interesting new variants, whilst providing a comparison. We firstly propose a one-dimensional second order Cauchy difference prior, and construct new first and second order two-dimensional isotropic Cauchy difference priors. Another new Cauchy prior is based on the stochastic partial differential equation approach, derived from Matérn type Gaussian presentation. The comparison also includes Cauchy sheets. Our numerical computations are based on both maximum a posteriori and conditional mean estimation.We exploit state-of-the-art MCMC methodologies such as Metropolis-within-Gibbs, Repelling-Attracting Metropolis, and No-U-Turn sampler variant of Hamiltonian Monte Carlo. We demonstrate the models and methods constructed for one-dimensional and two-dimensional deconvolution problems. Thorough MCMC statistics are provided for all test cases, including potential scale reduction factors.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bayesian quantification for coherent anti-Stokes Raman scattering spectroscopy

We propose a Bayesian statistical model for analyzing coherent anti-Stokes Raman scattering (CARS) spectra. Our quantitative analysis includes statistical estimation of constituent line-shape parameters, underlying Raman signal, error-corrected CARS spectrum, and the measured CARS spectrum. As such, this work enables extensive uncertainty quantification in the context of CARS spectroscopy. Furthermore, we present an unsupervised method for improving spectral resolution of Raman-like spectra requiring little to no \textit{a priori} information. Finally, the recently-proposed wavelet prism method for correcting the experimental artefacts in CARS is enhanced by using interpolation techniques for wavelets. The method is validated using CARS spectra of adenosine mono-, di-, and triphosphate in water, as well as, equimolar aqueous solutions of D-fructose, D-glucose, and their disaccharide combination sucrose.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Blind hierarchical deconvolution

Deconvolution is a fundamental inverse problem in signal processing and the prototypical model for recovering a signal from its noisy measurement. Nevertheless, the majority of model-based inversion techniques require knowledge on the convolution kernel to recover an accurate reconstruction and additionally prior assumptions on the regularity of the signal are needed. To overcome these limitations, we parametrise the convolution kernel and prior length-scales, which are then jointly estimated in the inversion procedure. The proposed framework of blind hierarchical deconvolution enables accurate reconstructions of functions with varying regularity and unknown kernel size and can be solved efficiently with an empirical Bayes two-step procedure, where hyperparameters are first estimated by optimisation and other unknowns then by an analytical formula.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Brexit Risk Implied by the SABR Martingale Defect in the EUR-GBP Smile

We construct a data-driven statistical indicator for quantifying the tail risk perceived by the EURGBP option market surrounding Brexit-related events. We show that under lognormal SABR dynamics this tail risk is closely related to the so-called martingale defect and provide a closed-form expression for this defect which can be computed by solving an inverse calibration problem. In order to cope with the the uncertainty which is inherent to this inverse problem, we adopt a Bayesian statistical parameter estimation perspective. We probe the resulting posterior densities with a combination of optimization and adaptive Markov chain Monte Carlo methods, thus providing a careful uncertainty estimation for all of the underlying parameters and the martingale defect indicator. Finally, to support the feasibility of the proposed method, we provide a Brexit "fever curve" for the year 2019.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Enhancing Industrial X-ray Tomography by Data-Centric Statistical Methods

X-ray tomography has applications in various industrial fields such as sawmill industry, oil and gas industry, chemical engineering, and geotechnical engineering. In this article, we study Bayesian methods for the X-ray tomography reconstruction. In Bayesian methods, the inverse problem of tomographic reconstruction is solved with help of a statistical prior distribution which encodes the possible internal structures by assigning probabilities for smoothness and edge distribution of the object. We compare Gaussian random field priors, that favour smoothness, to non-Gaussian total variation, Besov, and Cauchy priors which promote sharp edges and high-contrast and low-contrast areas in the object. We also present computational schemes for solving the resulting high-dimensional Bayesian inverse problem with 100,000-1,000,000 unknowns. In particular, we study the applicability of a no-U-turn variant of Hamiltonian Monte Carlo methods and of a more classical adaptive Metropolis-within-Gibbs algorithm for this purpose. These methods also enable full uncertainty quantification of the reconstructions. For faster computations, we use maximum a posteriori estimates with limited-memory BFGS optimisation algorithm. As the first industrial application, we consider sawmill industry X-ray log tomography. The logs have knots, rotten parts, and even possibly metallic pieces, making them good examples for non-Gaussian priors. Secondly, we study drill-core rock sample tomography, an example from oil and gas industry. We show that Cauchy priors produce smaller number of artefacts than other choices, especially with sparse high-noise measurements, and choosing Hamiltonian Monte Carlo enables systematic uncertainty quantification.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-Stationary Multi-layered Gaussian Priors for Bayesian Inversion

In this article, we study Bayesian inverse problems with multi-layered Gaussian priors. We first describe the conditionally Gaussian layers in terms of a system of stochastic partial differential equations. We build the computational inference method using a finite-dimensional Galerkin method. We show that the proposed approximation has a convergence-in-probability property to the solution of the original multi-layered model. We then carry out Bayesian inference using the preconditioned Crank--Nicolson algorithm which is modified to work with multi-layered Gaussian fields. We show via numerical experiments in signal deconvolution and computerized X-ray tomography problems that the proposed method can offer both smoothing and edge preservation at the same time.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A Bayesian-based approach to improving acoustic Born waveform inversion of seismic data for viscoelastic media

In seismic waveform inversion, the reconstruction of the subsurface properties is usually carried out using approximative wave propagation models to ensure computational efficiency. The viscoelastic nature of the subsurface is often unaccounted for, and two popular approximations--the acoustic and linearized Born inversion--are widely used. This leads to reconstruction errors since the approximations ignore realistic (physical) aspects of seismic wave propagation in the heterogeneous earth. In this study, we show that the Bayesian approximation error approach can be used to partially recover from errors, addressing elastic and viscous effects in acoustic Born inversion for viscoelastic media. The results of numerical examples indicate that neglecting the modelling errors induced by the approximations results in very poor recovery of the subsurface velocity fields.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cauchy difference priors for edge-preserving Bayesian inversion with an application to X-ray tomography

We study Cauchy-distributed difference priors for edge-preserving Bayesian statistical inverse problems. On the contrary to the well-known total variation priors, one-dimensional Cauchy priors are non-Gaussian priors also in the discretization limit. Cauchy priors have independent and identically distributed increments. One-dimensional Cauchy and Gaussian random walks are special cases of Lévy $α$-stable random walks with $α=1$ and $α=2$, respectively. Both random walks can be written in closed-form, and as priors, they provide smoothing and edge-preserving properties. We briefly discuss also continuous and discrete Lévy $α$-stable random walks, and generalize the methodology to two-dimensional priors. We apply the developed algorithm to one-dimensional deconvolution and two-dimensional X-ray tomography problems. We compute conditional mean estimates with single-component Metropolis-Hastings and maximum a posteriori estimates with Gauss-Newton-type optimization method. We compare the proposed tomography reconstruction method to filtered back-projection estimate and conditional mean estimates with Gaussian and total variation priors.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Elliptic boundary value problems with Gaussian white noise loads

Linear second order elliptic boundary value problems (BVP) on bounded Lipschitz domains are studied in the case of Gaussian white noise loads. Especially, Neumann and Robin BVPs are considered. The main obstacle for applying the usual variational approach is that the Gaussian white noise has irregular realizations. In particular, the corresponding Neumann boundary values are not well-defined in the ordinary sense. In this work, the BVP is reformulated by replacing the continuity of the boundary trace mapping with measurability. Instead of using variational methods alone, the reformulation of the BVP derives also from Cameron-Martin space techniques. The reformulation essentially returns the study of irregular white noise loads to study of $L^2$-loads. Admissibility of the reformulation is demonstrated by showing that usual finite element approximations of the BVP with discretized white noise loads converge to the solution of the reformulated problem. For Neumann and Robin BVPs, the finite dimensional approximations have been utilized before. However, here also the infinite-dimensional limit is considered.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sparse approximations of fractional Matérn fields

We consider a fast approximation method for a solution of a certain stochastic non-local pseudodifferential equation. This equation defines a Matérn class random field. The approximation method is based on the spectral compactness of the solution. We approximate the pseudodifferential operator with a Taylor expansion. By truncating the expansion, we can construct an approximation with Gaussian Markov random fields. We show that the solution of the truncated version can be constructed with an over-determined system of stochastic matrix equations with sparse matrices. We solve the system of equations with a sparse Cholesky decomposition. We consider the convergence of the discrete approximation of the solution to the continuous one. Finally numerical examples are given.

0 citations0 reviewsNo code linkedOpen access