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

Otmar Scherzer

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

ResearcherUnclaimed source record
math.APmath.NAmath.OCmath.FAComputer Visionmath-phmath.MPNumerical Analysismath.DGmath.SPmath.STphysics.comp-phphysics.flu-dynphysics.med-phphysics.opticsStatistics TheorySystems and Control

Catalog footprint

What is connected

44works
17topics
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

44 published item(s)

preprint2022arXiv

Universal approximation properties of shallow quadratic neural networks

In this paper we study shallow neural network functions which are linear combinations of compositions of activation and quadratic functions, replacing standard affine linear functions, often called neurons. We show the universality of this approximation and prove convergence rates results based on the theory of wavelets and statistical learning. We show for simple test cases that this ansatz requires a smaller numbers of neurons than standard affine linear neural networks. Moreover, we investigate the efficiency of this approach for clustering tasks with the MNIST data set. Similar observations are made when comparing deep (multi-layer) networks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A data-driven iteratively regularized Landweber iteration

We derive and analyse a new variant of the iteratively regularized Landweber iteration, for solving linear and nonlinear ill-posed inverse problems. The method takes into account training data, which are used to estimate the interior of a black box, which is used to define the iteration process. We prove convergence and stability for the scheme in infinite dimensional Hilbert spaces. These theoretical results are complemented by several numerical experiments for solving linear inverse problems for the Radon transform and a nonlinear inverse problem for Schlieren tomography.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Adjoint-state method for Hybridizable Discontinuous Galerkin discretization, application to the inverse acoustic wave problem

In this paper, we perform non-linear minimization using the Hybridizable Discontinuous Galerkin method (HDG) for the discretization of the forward problem, and implement the adjoint-state method for the computation of the functional derivatives. Compared to continuous and discontinuous Galerkin discretizations, HDG reduces the computational cost by working with the numerical traces, hence removing the degrees of freedom that are inside the cells. It is particularly attractive for large-scale time-harmonic quantitative inverse problems which make repeated use of the forward discretization as they rely on an iterative minimization procedure. HDG is based upon two levels of linear problems: a global system to find the numerical traces, followed by local systems to construct the volume solution. This technicality requires a careful derivation of the adjoint-state method, that we address in this paper. We work with the acoustic wave equations in the frequency domain and illustrate with a three-dimensional experiment using partial reflection-data, where we further employ the features of DG-like methods to efficiently handle the topography with p-adaptivity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inverse Problems of Single Molecule Localization Microscopy

Single molecule localization microscopy is a recently developed superresolution imaging technique to visualize structural properties of single cells. The basic principle consists in chemically attaching fluorescent dyes to the molecules, which after excitation with a strong laser may emit light. To achieve superresolution, signals of individual fluorophores are separated in time. In this paper we follow the physical and chemical literature and derive mathematical models describing the propagation of light emitted from dyes in single molecule localization microscopy experiments via Maxwell's equations. This forms the basis of formulating inverse problems related to single molecule localization microscopy. We also show that the current status of reconstruction methods is a simplification of more general inverse problems for Maxwell's equations as discussed here.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Asymptotic Expansions for Higher Order Elliptic Equations with an Application to Quantitative Photoacoustic Tomography

In this paper, we derive new asymptotic expansions for the solutions of higher order elliptic equations in the presence of small inclusions. As a byproduct, we derive a topological derivative based algorithm for the reconstruction of piecewise smooth functions. This algorithm can be used for edge detection in imaging, topological optimization, and for inverse problems, such as Quantitative Photoacoustic Tomography, for which we demonstrate the effectiveness of our asymptotic expansion method numerically.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Eigenvector Model Descriptors for Solving an Inverse Problem of Helmholtz Equation: Extended Materials

We study the seismic inverse problem for the recovery of subsurface properties in acoustic media. In order to reduce the ill-posedness of the problem, the heterogeneous wave speed parameter to be recovered is represented using a limited number of coefficients associated with a basis of eigenvectors of a diffusion equation, following the regularization by discretization approach. We compare several choices for the diffusion coefficient in the partial differential equations, which are extracted from the field of image processing. We first investigate their efficiency for image decomposition (accuracy of the representation with respect to the number of variables and denoising). Next, we implement the method in the quantitative reconstruction procedure for seismic imaging, following the Full Waveform Inversion method, where the difficulty resides in that the basis is defined from an initial model where none of the actual structures is known. In particular, we demonstrate that the method is efficient for the challenging reconstruction of media with salt-domes. We employ the method in two and three-dimensional experiments and show that the eigenvector representation compensates for the lack of low frequency information, it eventually serves us to extract guidelines for the implementation of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Preconditioning Inverse Problems for Hyperbolic Equations with Applications to Photoacoustic Tomography

This paper is concerned with robust preconditioning of wave equations constrained linear inverse problems from boundary observation data. The main result of this paper is a concept for regularization parameter robust preconditioning. Analogous concepts have been developed for control problems based on elliptic partial equations before.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convergence of Tikhonov regularization for solving ill-posed operator equations with solutions defined on surfaces

We study Tikhonov regularization for solving ill--posed operator equations where the solutions are functions defined on surfaces. One contribution of this paper is an error analysis of Tikhonov regularization which takes into account perturbations of the surfaces, in particular when the surfaces are approximated by spline surfaces. Another contribution is that we highlight the analysis of regularization for functions with range in vector bundles over surfaces. We also present some practical applications, such as an inverse problem of gravimetry and an imaging problem for denoising vector fields on surfaces, and show the numerical verification.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamical optical flow of saliency maps for predicting visual attention

Saliency maps are used to understand human attention and visual fixation. However, while very well established for static images, there is no general agreement on how to compute a saliency map of dynamic scenes. In this paper we propose a mathematically rigorous approach to this prob- lem, including static saliency maps of each video frame for the calculation of the optical flow. Taking into account static saliency maps for calculating the optical flow allows for overcoming the aperture problem. Our ap- proach is able to explain human fixation behavior in situations which pose challenges to standard approaches, such as when a fixated object disappears behind an occlusion and reappears after several frames. In addition, we quantitatively compare our model against alternative solutions using a large eye tracking data set. Together, our results suggest that assessing optical flow information across a series of saliency maps gives a highly accurate and useful account of human overt attention in dynamic scenes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Inverse Problems of Combined Photoacoustic and Optical Coherence Tomography

Optical coherence tomography (OCT) and photoacoustic tomography (PAT) are emerging non-invasive biological and medical imaging techniques. It is a recent trend in experimental science to design experiments that perform PAT and OCT imaging at once. In this paper we present a mathematical model describing the dual experiment. Since OCT is mathematically modelled by Maxwell's equations or some simplifications of it, whereas the light propagation in quantitative photoacoustics is modelled by (simplifications of) the radiative transfer equation, the first step in the derivation of a mathematical model of the dual experiment is to obtain a unified mathematical description, which in our case are Maxwell's equations. As a by-product we therefore derive a new mathematical model of photoacoustic tomography based on Maxwell's equations. It is well known by now, that without additional assumptions on the medium, it is not possible to uniquely reconstruct all optical parameters from either one of these modalities alone. We show that in the combined approach one has additional information, compared to a single modality, and the inverse problem of reconstruction of the optical parameters becomes feasible.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Photoacoustic Tomography With Spatially Varying Compressibility and Density

This paper investigates photoacoustic tomography with two spatially varying acoustic parameters, the compressibility and the density. We consider the reconstruction of the absorption density parameter (imaging parameter of Photoacoustics) with complete and partial measurement data. We investigate and analyze three different numerical methods for solving the imaging problem and compare the results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Singular Values of the Attenuated Photoacoustic Imaging Operator

We analyse the ill-posedness of the photoacoustic imaging problem in the case of an attenuating medium. To this end, we introduce an attenuated photoacoustic operator and determine the asymptotic behaviour of its singular values. Dividing the known attenuation models into strong and weak attenuation classes, we show that for strong attenuation, the singular values of the attenuated photoacoustic operator decay exponentially, and in the weak attenuation case the singular values of the attenuated photoacoustic operator decay with the same rate as the singular values of the non-attenuated photoacoustic operator.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Direct Method for Photoacoustic Tomography with Inhomogeneous Sound Speed

The standard approach for photoacoustic imaging with variable speed of sound is time reversal, which consists in solving a well-posed final-boundary value problem for the wave equation backwards in time. This paper investigates the iterative Landweber regularization algorithm, where convergence is guaranteed by standard regularization theory, notably also in cases of trapping sound speed or for short measurement times. We formulate and solve the direct and inverse problem on the whole Euclidean space, what is common in standard photoacoustic imaging, but not for time-reversal algorithms, where the problems are considered on a domain enclosed by the measurement devices. We formulate both the direct and adjoint photoacoustic operator as the solution of an interior and an exterior differential equation which are coupled by transmission conditions. The prior is solved numerically using a Galerkin scheme in space and finite difference discretization in time, while the latter consists in solving a boundary integral equation. We therefore use a BEM-FEM approach for numerical solution of the forward operators. We analyze this method, prove convergence, and provide numerical tests. Moreover, we compare the approach to time reversal.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A variational method for quantitative photoacoustic tomography with piecewise constant coefficients

In this article, we consider the inverse problem of determining spatially heterogeneous absorption and diffusion coefficients from a single measurement of the absorbed energy (in the steady-state diffusion approximation of light transfer). This problem, which is central in quantitative photoacoustic tomography, is in general ill-posed since it admits an infinite number of solution pairs. We show that when the coefficients are known to be piecewise constant functions, a unique solution can be obtained. For the numerical determination of the coefficients, we suggest a variational method based based on an Ambrosio-Tortorelli-approximation of a Mumford-Shah-like functional, which we implemented numerically and tested on simulated two-dimensional data.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Computational Realization of a Non-Equidistant Grid Sampling in Photoacoustics with a Non-Uniform FFT

To obtain the initial pressure from the collected data on a planar sensor arrangement in Photoacoustic tomography, there exists an exact analytic frequency domain reconstruction formula. An efficient realization of this formula needs to cope with the evaluation of the datas Fourier transform on a non-equispaced mesh. In this paper, we use the non-uniform fast Fourier transform to handle this issue and show its feasibility in 3D experiments. This is done in comparison to the standard approach that uses polynomial interpolation. Moreover, we investigate the effect and the utility of flexible sensor location on the quality of photoacoustic image reconstruction. The computational realization is accomplished by the use of a multi-dimensional non-uniform fast Fourier algorithm, where non-uniform data sampling is performed both in frequency and spatial domain. We show that with appropriate sampling the imaging quality can be significantly improved. Reconstructions with synthetic and real data show the superiority of this method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Non-Equispaced Grid Sampling in Photoacoustics with a Non-Uniform FFT

To obtain the initial pressure from the collected data on a planar sensor arrangement in photoacoustic tomography, there exists an exact analytic frequency domain reconstruction formula. An efficient realization of this formula needs to cope with the evaluation of the data's Fourier transform on a non-equispaced mesh. In this paper, we use the non-uniform fast Fourier transform to handle this issue and show its feasibility in 3D experiments with real and synthetic data. This is done in comparison to the standard approach that uses linear, polynomial or nearest neighbor interpolation. Moreover, we investigate the effect and the utility of flexible sensor location to make optimal use of a limited number of sensor points. The computational realization is accomplished by the use of a multi-dimensional non-uniform fast Fourier algorithm, where non-uniform data sampling is performed both in frequency and spatial domain. Examples with synthetic and real data show that both approaches improve image quality.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optical Flow on Evolving Sphere-Like Surfaces

In this work we consider optical flow on evolving Riemannian 2-manifolds which can be parametrised from the 2-sphere. Our main motivation is to estimate cell motion in time-lapse volumetric microscopy images depicting fluorescently labelled cells of a live zebrafish embryo. We exploit the fact that the recorded cells float on the surface of the embryo and allow for the extraction of an image sequence together with a sphere-like surface. We solve the resulting variational problem by means of a Galerkin method based on vector spherical harmonics and present numerical results computed from the aforementioned microscopy data.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimal Convergence Rates Results for Linear Inverse Problems in Hilbert Spaces

In this paper, we prove optimal convergence rates results for regularisation methods for solving linear ill-posed operator equations in Hilbert spaces. The result generalises existing convergence rates results on optimality to general source conditions, such as logarithmic source conditions. Moreover, we also provide optimality results under variational source conditions and show the connection to approximative source conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantitative photoacoustic tomography with piecewise constant material parameters

The goal of quantitative photoacoustic tomography is to determine optical and acoustical material properties from initial pressure maps as obtained, for instance, from photoacoustic imaging. The most relevant parameters are absorption, diffusion and Grueneisen coefficients, all of which can be heterogeneous. Recent work by Bal and Ren shows that in general, unique reconstruction of all three parameters is impossible, even if multiple measurements of the initial pressure (corresponding to different laser excitation directions at a single wavelength) are available. Here, we propose a restriction to piecewise constant material parameters. We show that in the diffusion approximation of light transfer, piecewise constant absorption, diffusion and Grüneisen coefficients can be recovered uniquely from photoacoustic measurements at a single wavelength. In addition, we implemented our ideas numerically and tested them on simulated three-dimensional data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Decomposition of Optical Flow on the Sphere

We propose a number of variational regularisation methods for the estimation and decomposition of motion fields on the $2$-sphere. While motion estimation is based on the optical flow equation, the presented decomposition models are motivated by recent trends in image analysis. In particular we treat $u+v$ decomposition as well as hierarchical decomposition. Helmholtz decomposition of motion fields is obtained as a natural by-product of the chosen numerical method based on vector spherical harmonics. All models are tested on time-lapse microscopy data depicting fluorescently labelled endodermal cells of a zebrafish embryo.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Generalized Convergence Rates Results for Linear Inverse Problems in Hilbert Spaces

In recent years, a series of convergence rates conditions for regularization methods has been developed. Mainly, the motivations for developing novel conditions came from the desire to carry over convergence rates results from the Hilbert space setting to generalized Tikhonov regularization in Banach spaces. For instance, variational source conditions have been developed and they were expected to be equivalent to standard source conditions for linear inverse problems in a Hilbert space setting. We show that this expectation does not hold. However, in the standard Hilbert space setting these novel conditions are optimal, which we prove by using some deep results from Neubauer, and generalize existing convergence rates results. The key tool in our analysis is a novel source condition, which we put into relation to the existing source conditions from the literature. As a positive by-product, convergence rates results can be proven without spectral theory, which is the standard technique for proving convergence rates for linear inverse problems in Hilbert spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Inverse boundary value problem for the Helmholtz equation: Multi-level approach and iterative reconstruction

We study the inverse boundary value problem for the Helmholtz equation using the Dirichlet-to-Neumann map at selected frequency as the data. We develop an explicit reconstruction of the wavespeed using a multi-level nonlinear projected steepest descent iterative scheme in Banach spaces. We consider wavespeeds containing discontinuities. A conditional Lipschitz stability estimate for the inverse problem holds for wavespeeds of the form of a linear combination of piecewise constant functions with an underlying domain partitioning, and gives a framework in which the scheme converges. The stability constant grows exponentially as the number of subdomains in the domain partitioning increases. To mitigate this growth of the stability constant, we introduce hierarchical compressive approximations of the solution to the inverse problem with piecewise constant functions. We establish an optimal bound of the stability constant, which leads to a condition on the compression rate pertaining to these approximations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mathematical Modelling of Optical Coherence Tomography

In this chapter a general mathematical model of Optical Coherence Tomography (OCT) is presented on the basis of the electromagnetic theory. OCT produces high resolution images of the inner structure of biological tissues. Images are obtained by measuring the time delay and the intensity of the backscattered light from the sample considering also the coherence properties of light. The scattering problem is considered for a weakly scattering medium located far enough from the detector. The inverse problem is to reconstruct the susceptibility of the medium given the measurements for different positions of the mirror. Different approaches are addressed depending on the different assumptions made about the optical properties of the sample. This procedure is applied to a full field OCT system and an extension to standard (time and frequency domain) OCT is briefly presented.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Optical Flow on Evolving Surfaces with Space and Time Regularisation

We extend the concept of optical flow with spatiotemporal regularisation to a dynamic non-Euclidean setting. Optical flow is traditionally computed from a sequence of flat images. The purpose of this paper is to introduce variational motion estimation for images that are defined on an evolving surface. Volumetric microscopy images depicting a live zebrafish embryo serve as both biological motivation and test data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stability in the linearized problem of quantitative elastography

The goal of quantitative elastography is to identify biomechanical parameters from interior displacement data, which are provided by other modalities, such as ultrasound or magnetic resonance imaging. In this paper, we analyze the stability of several linearized problems in quantitative elastography. Our method is based on the theory of redundant systems of linear partial differential equations. We analyze the ellipticity properties of the corresponding PDE systems augmented with the interior displacement data; we explicitly characterize the kernel of the forward operators and show injectivity for particular linearizations. Stability criteria can then be deduced. Our results show stability of shear modulus, pressure and density; they indicate that singular strain fields should be avoided, and show how additional measurements can help in ensuring stability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Texture Generation for Photoacoustic Elastography

Elastographic imaging is a widely used technique which can in principle be implemented on top of every imaging modality. In elastography, the specimen is exposed to a force causing local displacements in the probe, and imaging is performed before and during the displacement experiment. From the computed displacements material parameters can be deduced, which in turn can be used for clinical diagnosis. Photoacoustic imaging is an emerging image modality, which exhibits functional and morphological contrast. However, opposed to ultrasound imaging, for instance, it is considered a modality which is not suited for elastography, because it does not reveal speckle patterns. However, this is somehow counter-intuitive, because photoacoustic imaging makes available the whole frequency spectrum as opposed to single frequency standard ultrasound imaging. In this work, we show that in fact artificial speckle patterns can be introduced by using only a band-limited part of the measurement data. We also show that after introduction of artificial speckle patterns, deformation estimation can be implemented more reliably in photoacoustic imaging.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An Approach to the Minimization of the Mumford-Shah Functional using Γ-convergence and Topological Asymptotic Expansion

In this paper, we present a method for the numerical minimization of the Mumford-Shah functional that is based on the idea of topological asymptotic expansions. The basic idea is to cover the expected edge set with balls of radius ε> 0 and use the number of balls, multiplied with 2ε, as an estimate for the length of the edge set. We introduce a functional based on this idea and prove that it converges in the sense of Γ-limits to the Mumford-Shah functional. Moreover, we show that ideas from topological asymptotic analysis can be used for determining where to position the balls covering the edge set. The results of the proposed method are presented by means of two numerical examples and compared with the results of the classical approximation due to Ambrosio and Tortorelli.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exact Solutions of One-Dimensional TGV

L2-TGV-regularization has been introduced by Bredies, Kunisch, and Pock. This regularization method requires careful tuning of two regularization parameters. The focus of this paper is to derive analytical results, which allow for characterizing parameter settings, which make this method in fact different from L2-TV (the ROF-model) and L2-TV2 regularization, respectively. In this paper we also provide explicit solutions of TGV-denoising for particular one-dimensional function data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Optical Flow on Evolving Surfaces with an Application to the Analysis of 4D Microscopy Data

We extend the concept of optical flow to a dynamic non-Euclidean setting. Optical flow is traditionally computed from a sequence of flat images. It is the purpose of this paper to introduce variational motion estimation for images that are defined on an evolving surface. Volumetric microscopy images depicting a live zebrafish embryo serve as both biological motivation and test data.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A convergence analysis of a multi-level projected steepest descent iteration for nonlinear inverse problems in Banach spaces subject to stability constraints

We consider nonlinear inverse problems described by operator equations in Banach spaces. Assuming conditional stability of the inverse problem, that is, assuming that stability holds on a closed, convex subset of the domain of the operator, we introduce a novel nonlinear projected steepest descent iteration and analyze its convergence to an approximate solution given limited accuracy data. We proceed with developing a multi-level algorithm based on a nested family of closed, convex subsets on which stability holds and the stability constants are ordered. Growth of the stability constants is coupled to the increase in accuracy of approximation between neighboring levels to ensure that the algorithm can continue from level to level until the iterate satisfies a desired discrepancy criterion, after a finite number of steps.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Local Analysis of Inverse Problems: Hölder Stability and Iterative Reconstruction

We consider a class of inverse problems defined by a nonlinear map from parameter or model functions to the data. We assume that solutions exist. The space of model functions is a Banach space which is smooth and uniformly convex; however, the data space can be an arbitrary Banach space. We study sequences of parameter functions generated by a nonlinear Landweber iteration and conditions under which these strongly converge, locally, to the solutions within an appropriate distance. We express the conditions for convergence in terms of Hölder stability of the inverse maps, which ties naturally to the analysis of inverse problems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonparametric instrumental regression with non-convex constraints

This paper considers the nonparametric regression model with an additive error that is dependent on the explanatory variables. As is common in empirical studies in epidemiology and economics, it also supposes that valid instrumental variables are observed. A classical example in microeconomics considers the consumer demand function as a function of the price of goods and the income, both variables often considered as endogenous. In this framework, the economic theory also imposes shape restrictions on the demand function, like integrability conditions. Motivated by this illustration in microeconomics, we study an estimator of a nonparametric constrained regression function using instrumental variables by means of Tikhonov regularization. We derive rates of convergence for the regularized model both in a deterministic and stochastic setting under the assumption that the true regression function satisfies a projected source condition including, because of the non-convexity of the imposed constraints, an additional smallness condition.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Photoacoustic imaging in attenuating acoustic media based on strongly causal models

In this paper we derive time reversal imaging functionals for two strongly causal acoustic attenuation models, which have been proposed recently. The time reversal techniques are based on recently proposed ideas of Ammari et al for the thermo-viscous wave equation. Here and there an asymptotic analysis provides reconstruction functionals from first order corrections for the attenuating effect. In addition, we present a novel approach for higher order corrections.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Reconstruction Formulas for a Single Scattering Model in Photoacoustic Imaging and Applications to Sectional Imaging

There has been devoted significant mathematical research to model the light propagation in tissue and to recover the absorption and scattering coefficients after photoacoustic inversion. Typically, the basic light propagation models considered there are the diffusion limits of the radiative heat transfer model. These equations are well suited for models where the elastic scattering is the dominant effect. If, however, the scattering is less pronounced, a single scattering approximation model is practicable. As we show in this paper, this approach is practically relevant in focused/sectional photoacoustic imaging. In this paper, we study analytical reconstruction formulas for the single scattering case. To realise the single scattering approach, we propose concrete physical experiments based on photoacoustical sectional imaging.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Levenberg-Marquardt Iteration for Numerical Inversion of the Power Density Operator

In this paper we develop a convergence analysis in an infinite dimensional setting of the Levenberg-Marquardt iteration for the solution of a hybrid conductivity imaging problem. The problem consists in determining the spatially varying conductivity $σ$ from a series of measurements of power densities for various voltage inductions. Although this problem has been very well studied in the literature, convergence and regularizing properties of iterative algorithms in an infinite dimensional setting are still rudimentary. We provide a partial result under the assumptions that the derivative of the operator, mapping conductivities to power densities, is injective and the data is noise-free. Moreover, we implemented the Levenberg-Marquardt algorithm and tested it on simulated data.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Convergence of Variational Regularization Methods for Imaging on Riemannian Manifolds

We consider abstract operator equations $Fu=y$, where $F$ is a compact linear operator between Hilbert spaces $U$ and $V$, which are function spaces on \emph{closed, finite dimensional Riemannian manifolds}, respectively. This setting is of interest in numerous applications such as Computer Vision and non-destructive evaluation. In this work, we study the approximation of the solution of the ill-posed operator equation with Tikhonov type regularization methods. We prove well-posedness, stability, convergence, and convergence rates of the regularization methods. Moreover, we study in detail the numerical analysis and the numerical implementation. Finally, we provide for three different inverse problems numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Reconstruction Formulas for Photoacoustic Sectional Imaging

The literature on reconstruction formulas for photoacoustic tomography (PAT) is vast. The various reconstruction formulas differ by used measurement devices and geometry on which the data are sampled. In standard photoacoustic imaging (PAI), the object under investigation is illuminated uniformly. Recently, sectional photoacoustic imaging techniques, using focusing techniques for initializing and measuring the pressure along a plane, appeared in the literature. This paper surveys existing and provides novel exact reconstruction formulas for sectional photoacoustic imaging.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Simultaneous Reconstructions of Absorption Density and Wave Speed with Photoacoustic Measurements

In this paper we propose an approach for \emph{simultaneous} identification of the \emph{absorption density} and the \emph{speed of sound} by photoacoustic measurements. Experimentally our approach can be realized with sliced photoacoustic experiments. The mathematical model for such an experiment is developed and exact reconstruction formulas for both parameters are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Residual Method for Regularizing Ill-Posed Problems

Although the \emph{residual method}, or \emph{constrained regularization}, is frequently used in applications, a detailed study of its properties is still missing. This sharply contrasts the progress of the theory of Tikhonov regularization, where a series of new results for regularization in Banach spaces has been published in the recent years. The present paper intends to bridge the gap between the existing theories as far as possible. We develop a stability and convergence theory for the residual method in general topological spaces. In addition, we prove convergence rates in terms of (generalized) Bregman distances, which can also be applied to non-convex regularization functionals. We provide three examples that show the applicability of our theory. The first example is the regularized solution of linear operator equations on $L^p$-spaces, where we show that the results of Tikhonov regularization generalize unchanged to the residual method. As a second example, we consider the problem of density estimation from a finite number of sampling points, using the Wasserstein distance as a fidelity term and an entropy measure as regularization term. It is shown that the densities obtained in this way depend continuously on the location of the sampled points and that the underlying density can be recovered as the number of sampling points tends to infinity. Finally, we apply our theory to compressed sensing. Here, we show the well-posedness of the method and derive convergence rates both for convex and non-convex regularization under rather weak conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Photoacoustic Imaging Taking into Account Attenuation

First, we review existing attenuation models and discuss their causality properties, which we believe to be essential for algorithms for inversion with attenuated data. Then, we survey causality properties of common attenuation models. We also derive integro-differential equations which the attenuated waves are satisfying. In addition we discuss the ill--conditionness of the inverse problem for calculating the unattenuated wave from the attenuated one.

0 citations0 reviewsNo code linkedOpen access