0000000000065791

AUTHOR

Matti Lassas

0000-0003-2043-3156

Uniqueness, reconstruction and stability for an inverse problem of a semi-linear wave equation

We consider the recovery of a potential associated with a semi-linear wave equation on Rn+1, n > 1. We show that an unknown potential a(x, t) of the wave equation ???u + aum = 0 can be recovered in a H & ouml;lder stable way from the map u|onnx[0,T] ???-> (11, avu|ac >= x[0,T])L2(oc >= x[0,T]). This data is equivalent to the inner product of the Dirichlet-to-Neumann map with a measurement function ???. We also prove similar stability result for the recovery of a when there is noise added to the boundary data. The method we use is constructive and it is based on the higher order linearization. As a consequence, we also get a uniqueness result. We also give a detailed presentation of the forw…

research product

The Poisson embedding approach to the Calderón problem

We introduce a new approach to the anisotropic Calder\'on problem, based on a map called Poisson embedding that identifies the points of a Riemannian manifold with distributions on its boundary. We give a new uniqueness result for a large class of Calder\'on type inverse problems for quasilinear equations in the real analytic case. The approach also leads to a new proof of the result by Lassas and Uhlmann (2001) solving the Calder\'on problem on real analytic Riemannian manifolds. The proof uses the Poisson embedding to determine the harmonic functions in the manifold up to a harmonic morphism. The method also involves various Runge approximation results for linear elliptic equations.

research product

Curvelet-based method for orientation estimation of particles from optical images

A method based on the curvelet transform is introduced to estimate the orientation distribution from two-dimensional images of small anisotropic particles. Orientation of fibers in paper is considered as a particular application of the method. Theoretical aspects of the suitability of this method are discussed and its efficiency is demonstrated with simulated and real images of fibrous systems. Comparison is made with two traditionally used methods of orientation analysis, and the new curvelet-based method is shown to perform better than these tradi- tional methods. peerReviewed

research product

Uniqueness and stability of an inverse problem for a semi-linear wave equation

We consider the recovery of a potential associated with a semi-linear wave equation on $\mathbb{R}^{n+1}$, $n\geq 1$. We show a H\"older stability estimate for the recovery of an unknown potential $a$ of the wave equation $\square u +a u^m=0$ from its Dirichlet-to-Neumann map. We show that an unknown potential $a(x,t)$, supported in $\Omega\times[t_1,t_2]$, of the wave equation $\square u +a u^m=0$ can be recovered in a H\"older stable way from the map $u|_{\partial \Omega\times [0,T]}\mapsto \langle\psi,\partial_\nu u|_{\partial \Omega\times [0,T]}\rangle_{L^2(\partial \Omega\times [0,T])}$. This data is equivalent to the inner product of the Dirichlet-to-Neumann map with a measurement fun…

research product

The Calderon problem in transversally anisotropic geometries

We consider the anisotropic Calderon problem of recovering a conductivity matrix or a Riemannian metric from electrical boundary measurements in three and higher dimensions. In the earlier work \cite{DKSaU}, it was shown that a metric in a fixed conformal class is uniquely determined by boundary measurements under two conditions: (1) the metric is conformally transversally anisotropic (CTA), and (2) the transversal manifold is simple. In this paper we will consider geometries satisfying (1) but not (2). The first main result states that the boundary measurements uniquely determine a mixed Fourier transform / attenuated geodesic ray transform (or integral against a more general semiclassical…

research product

Inverse problems and invisibility cloaking for FEM models and resistor networks

In this paper we consider inverse problems for resistor networks and for models obtained via the finite element method (FEM) for the conductivity equation. These correspond to discrete versions of the inverse conductivity problem of Calderón. We characterize FEM models corresponding to a given triangulation of the domain that are equivalent to certain resistor networks, and apply the results to study nonuniqueness of the discrete inverse problem. It turns out that the degree of nonuniqueness for the discrete problem is larger than the one for the partial differential equation. We also study invisibility cloaking for FEM models, and show how an arbitrary body can be surrounded with a layer …

research product

Curvelet-based method for orientation estimation of particles

A method based on the curvelet transform is introduced for estimating from two-dimensional images the orientation distribution of small anisotropic particles. Orientation of fibers in paper is considered as a particular application of the method. Theoretical aspects of the suitability of this method are discussed and its efficiency is demonstrated with simulated and real images of fibrous systems. Comparison is made with two traditionally used methods of orientation analysis, and the new curvelet-based method is shown to perform clearly better than these traditional methods.

research product

The Linearized Calderón Problem in Transversally Anisotropic Geometries

In this article we study the linearized anisotropic Calderon problem. In a compact manifold with boundary, this problem amounts to showing that products of harmonic functions form a complete set. Assuming that the manifold is transversally anisotropic, we show that the boundary measurements determine an FBI type transform at certain points in the transversal manifold. This leads to recovery of transversal singularities in the linearized problem. The method requires a geometric condition on the transversal manifold related to pairs of intersecting geodesics, but it does not involve the geodesic X-ray transform which has limited earlier results on this problem.

research product

"Towards a "fingerprint" of paper network; separating forgeries from genuine by the properties of fibre structure"

A novel method is introduced for distinguishing counterfeit banknotes from genuine samples. The method is based on analyzing differences in the networks of paper fibers. The main tool is a curvelet-based algorithm for measuring the distribution of overall fiber orientation and quantifying its anisotropy. The use of a couple or more appropriate parameters makes it possible to distinguish forgeries from genuine samples as concentrated point clouds in a two- or three-dimensional parameter space. Furthermore, the techniques of making watermarks is investigated by comparing genuine and counterfeit €50 banknotes. In addition, the so-called wire markings are shown to differ significantly from each…

research product

The Calder\'on problem for the conformal Laplacian

We consider a conformally invariant version of the Calder\'on problem, where the objective is to determine the conformal class of a Riemannian manifold with boundary from the Dirichlet-to-Neumann map for the conformal Laplacian. The main result states that a locally conformally real-analytic manifold in dimensions $\geq 3$ can be determined in this way, giving a positive answer to an earlier conjecture by Lassas and Uhlmann (2001). The proof proceeds as in the standard Calder\'on problem on a real-analytic Riemannian manifold, but new features appear due to the conformal structure. In particular, we introduce a new coordinate system that replaces harmonic coordinates when determining the co…

research product

Towards a "fingerprint" of paper network; separating forgeries from genuine by the properties of fibre structure

A novel method is introduced for distinguishing counterfeit banknotes from genuine samples. The method is based on analyzing differences in the networks of paper fibers. The main tool is a curvelet-based algorithm for measuring the distribution of overall fiber orientation and quantifying its anisotropy. The use of a couple or more appropriate parameters makes it possible to distinguish forgeries from genuine samples as concentrated point clouds in a two- or three-dimensional parameter space. Furthermore, the techniques of making watermarks is investigated by comparing genuine and counterfeit e50 banknotes. In addition, the so-called wire markings are shown to differ significantly from each…

research product

The Calderón problem for the conformal Laplacian

We consider a conformally invariant version of the Calderón problem, where the objective is to determine the conformal class of a Riemannian manifold with boundary from the Dirichlet-to-Neumann map for the conformal Laplacian. The main result states that a locally conformally real-analytic manifold in dimensions can be determined in this way, giving a positive answer to an earlier conjecture [LU02, Conjecture 6.3]. The proof proceeds as in the standard Calderón problem on a real-analytic Riemannian manifold, but new features appear due to the conformal structure. In particular, we introduce a new coordinate system that replaces harmonic coordinates when determining the conformal class in a …

research product

An Inverse Problem for the Relativistic Boltzmann Equation

We consider an inverse problem for the Boltzmann equation on a globally hyperbolic Lorentzian spacetime $(M,g)$ with an unknown metric $g$. We consider measurements done in a neighbourhood $V\subset M$ of a timelike path $\mu$ that connects a point $x^-$ to a point $x^+$. The measurements are modelled by a source-to-solution map, which maps a source supported in $V$ to the restriction of the solution to the Boltzmann equation to the set $V$. We show that the source-to-solution map uniquely determines the Lorentzian spacetime, up to an isometry, in the set $I^+(x^-)\cap I^-(x^+)\subset M$. The set $I^+(x^-)\cap I^-(x^+)$ is the intersection of the future of the point $x^-$ and the past of th…

research product

Partial data inverse problems and simultaneous recovery of boundary and coefficients for semilinear elliptic equations

We study various partial data inverse boundary value problems for the semilinear elliptic equation $\Delta u+ a(x,u)=0$ in a domain in $\mathbb R^n$ by using the higher order linearization technique introduced in [LLS 19, FO19]. We show that the Dirichlet-to-Neumann map of the above equation determines the Taylor series of $a(x,z)$ at $z=0$ under general assumptions on $a(x,z)$. The determination of the Taylor series can be done in parallel with the detection of an unknown cavity inside the domain or an unknown part of the boundary of the domain. The method relies on the solution of the linearized partial data Calder\'on problem [FKSU09], and implies the solution of partial data problems fo…

research product

Using the fibre structure of paper to determine authenticity of the documents: analysis of transmitted light images of stamps and banknotes.

A novel method is presented for distinguishing postal stamp forgeries and counterfeit banknotes from genuine samples. The method is based on analyzing differences in paper fibre networks. The main tool is a curvelet-based algorithm for measuring overall fibre orientation distribution and quantifying anisotropy. Using a couple of more appropriate parameters makes it possible to distinguish forgeries from genuine originals as concentrated point clouds in two- or three-dimensional parameter space.

research product

Curvelet-based method for orientation estimation of particles from optical images

A method based on the curvelet transform is introduced to estimate the orientation distribution from two-dimensional images of small anisotropic particles. Orientation of fibers in paper is considered as a particular application of the method. Theoretical aspects of the suitability of this method are discussed and its efficiency is demonstrated with simulated and real images of fibrous systems. Comparison is made with two traditionally used methods of orientation analysis, and the new curvelet-based method is shown to perform better than these tradi- tional methods. © The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of th…

research product

Evaluation of the orientation distribution of fibers from reflection images of fibrous samples

We consider illumination systems and mathematical algorithms for determination of the anisotropy and topographical features of an illuminated surface from its reflection images. As a particular example we study determination of the fiber orientation of paper surface. We also consider illumination systems with multiple light sources, and introduce optimization algorithms that exploit different spectral bands of these light sources. We show that a system of three light sources, e.g., a blue, green and red LED placed in a regular triangular form, efficiently prevents distortion of the above features. It is also easy to implement in applications, e.g., of the paper industry. We furthermore show…

research product

Evaluation of the areal material distribution of paper from its optical transmission image

International audience; The goal of this study was to evaluate the areal mass distribution (defined as the X-ray transmission image) of paper from its optical transmission image. A Bayesian inversion framework was used in the related deconvolution process so as to combine indirect optical information with a priori knowledge about the type of paper imaged. The a priori knowledge was expressed in the form of an empirical Besov space prior distribution constructed in a computationally effective way using the wavelet transform. The estimation process took the form of a large-scale optimization problem, which was in turn solved using the gradient descent method of Barzilai and Borwein. It was de…

research product

Inverse problems for elliptic equations with power type nonlinearities

We introduce a method for solving Calder\'on type inverse problems for semilinear equations with power type nonlinearities. The method is based on higher order linearizations, and it allows one to solve inverse problems for certain nonlinear equations in cases where the solution for a corresponding linear equation is not known. Assuming the knowledge of a nonlinear Dirichlet-to-Neumann map, we determine both a potential and a conformal manifold simultaneously in dimension $2$, and a potential on transversally anisotropic manifolds in dimensions $n \geq 3$. In the Euclidean case, we show that one can solve the Calder\'on problem for certain semilinear equations in a surprisingly simple way w…

research product

Modelling and analysing oriented fibrous structures

Abstract. A mathematical model for fibrous structures using a direction dependent scaling law is presented. The orientation of fibrous nets (e.g. paper) is analysed with a method based on the curvelet transform. The curvelet-based orientation analysis has been tested successfully on real data from paper samples: the major directions of fibrefibre orientation can apparently be recovered. Similar results are achieved in tests on data simulated by the new model, allowing a comparison with ground truth. peerReviewed

research product