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

Nuutti Hyvönen

Nuutti Hyvönen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.APmath.NANumerical Analysismath.OC

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
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

16 published item(s)

preprint2026arXiv

Spotlight inversion by orthogonal projections

Many computational problems involve solving a linear system of equations, although only a subset of the entries of the solution are needed. In inverse problems, where the goal is to estimate unknown parameters from indirect noisy observations, it is not uncommon that the forward model linking the observed variables to the unknowns depends on variables that are not of primary interest, often referred to as nuisance parameters. In this article, we consider linear problems, and propose a novel projection technique to eliminate, or at least mitigate, the contribution of the nuisance parameters in the model. We refer to this approach as spotlight inversion, as it allows to focus on only the portion of primary interest of the unknown parameter vector, leaving the uninteresting part in the shadow. The viability of the approach is illustrated with two computed examples, one where it works as model reduction for a finite element approximation of an elliptic PDE, the other amounting to local fanbeam X-ray tomography, spotlighting the region of interest that is part of the full target.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reconstruction of singular and degenerate inclusions in Calderón's problem

We consider the reconstruction of the support of an unknown perturbation to a known conductivity coefficient in Calderón's problem. In a previous result by the authors on monotonicity-based reconstruction, the perturbed coefficient is allowed to simultaneously take the values $0$ and $\infty$ in some parts of the domain and values bounded away from $0$ and $\infty$ elsewhere. We generalise this result by allowing the unknown coefficient to be the restriction of an $A_2$-Muckenhoupt weight in parts of the domain, thereby including singular and degenerate behaviour in the governing equation. In particular, the coefficient may tend to $0$ and $\infty$ in a controlled manner, which goes beyond the standard setting of Calderón's problem. Our main result constructively characterises the outer shape of the support of such a general perturbation, based on a local Neumann-to-Dirichlet map defined on an open subset of the domain boundary.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Series reversion in Calderón's problem

This work derives explicit series reversions for the solution of Calderón's problem. The governing elliptic partial differential equation is $\nabla\cdot(A\nabla u)=0$ in a bounded Lipschitz domain and with a matrix-valued coefficient. The corresponding forward map sends $A$ to a projected version of a local Neumann-to-Dirichlet operator, allowing for the use of partial boundary data and finitely many measurements. It is first shown that the forward map is analytic, and subsequently reversions of its Taylor series up to specified orders lead to a family of numerical methods for solving the inverse problem with increasing accuracy. The convergence of these methods is shown under conditions that ensure the invertibility of the Fréchet derivative of the forward map. The introduced numerical methods are of the same computational complexity as solving the linearised inverse problem. The analogous results are also presented for the smoothened complete electrode model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Mimicking relative continuum measurements by electrode data in two-dimensional electrical impedance tomography

This paper introduces a constructive method for approximating relative continuum measurements in two-dimensional electrical impedance tomography based on data originating from either the point electrode model or the complete electrode model. The upper bounds for the corresponding approximation errors explicitly depend on the number (and size) of the employed electrodes as well as on the regularity of the continuum current that is mimicked. In particular, if the input current and the object boundary are infinitely smooth, the discrepancy associated with the point electrode model converges to zero faster than any negative power of the number of electrodes. The results are first proven for the unit disk via trigonometric interpolation and quadrature rules, and they are subsequently extended to more general domains with the help of conformal mappings.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Monotonicity-based reconstruction of extreme inclusions in electrical impedance tomography

The monotonicity-based approach has become one of the fundamental methods for reconstructing inclusions in the inverse problem of electrical impedance tomography. Thus far the method has not been proven to be able to handle extreme inclusions that correspond to some parts of the studied domain becoming either perfectly conducting or perfectly insulating. The main obstacle has arguably been establishing suitable monotonicity principles for the corresponding Neumann-to-Dirichlet boundary maps. In this work, we tackle this shortcoming by first giving a convergence result in the operator norm for the Neumann-to-Dirichlet map when the conductivity coefficient decays to zero and/or grows to infinity in some given parts of the domain. This allows passing the necessary monotonicity principles to the limiting case. Subsequently, we show how the monotonicity method generalizes to the definite case of reconstructing either perfectly conducting or perfectly insulating inclusions, as well as to the indefinite case where the perturbed conductivity can take any values between, and including, zero and infinity.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On regularity of the logarithmic forward map of electrical impedance tomography

This work considers properties of the logarithm of the Neumann-to-Dirichlet boundary map for the conductivity equation in a Lipschitz domain. It is shown that the mapping from the (logarithm of) the conductivity, i.e. the (logarithm of) the coefficient in the divergence term of the studied elliptic partial differential equation, to the logarithm of the Neumann-to-Dirichlet map is continuously Fréchet differentiable between natural topologies. Moreover, for any essentially bounded perturbation of the conductivity, the Fréchet derivative defines a bounded linear operator on the space of square integrable functions living on the domain boundary, although the logarithm of the Neumann-to-Dirichlet map itself is unbounded in that topology. In particular, it follows from the fundamental theorem of calculus that the difference between the logarithms of any two Neumann-to-Dirichlet maps is always bounded on the space of square integrable functions. All aforementioned results also hold if the Neumann-to-Dirichlet boundary map is replaced by its inverse, i.e. the Dirichlet-to-Neumann map.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Optimal depth-dependent distinguishability bounds for electrical impedance tomography in arbitrary dimension

The inverse problem of electrical impedance tomography is severely ill-posed. In particular, the resolution of images produced by impedance tomography deteriorates as the distance from the measurement boundary increases. Such depth dependence can be quantified by the concept of distinguishability of inclusions. This paper considers the distinguishability of perfectly conducting ball inclusions inside a unit ball domain, extending and improving known two-dimensional results to an arbitrary dimension $d \geq 2$ with the help of Kelvin transformations. The obtained depth-dependent distinguishability bounds are also proven to be optimal.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Efficient inclusion of total variation type priors in quantitative photoacoustic tomography

Quantitative photoacoustic tomography is an emerging imaging technique aimed at estimating the distribution of optical parameters inside tissues from photoacoustic images, which are formed by combining optical information and ultrasonic propagation. This optical parameter estimation problem is ill-posed and needs to be approached within the framework of inverse problems. Photoacoustic images are three-dimensional and high-resolution. Furthermore, high-resolution reconstructions of the optical parameters are targeted. Therefore, in order to provide a practical method for quantitative photoacoustic tomography, the inversion algorithm needs to be able to perform successfully with problems of prominent size. In this work, an efficient approach for the inverse problem of quantitative photoacoustic tomography is proposed, assuming an edge-preferring prior for the optical parameters. The method is based on iteratively combining priorconditioned LSQR with a lagged diffusivity step and a linearisation of the measurement model, with the needed multiplications by Jacobians performed in a matrix-free manner. The algorithm is tested with three-dimensional numerical simulations. The results show that the approach can be used to produce accurate and high quality estimates of absorption and diffusion in complex three-dimensional geometries with moderate computation time and cost.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Edge-promoting reconstruction of absorption and diffusivity in optical tomography

In optical tomography a physical body is illuminated with near-infrared light and the resulting outward photon flux is measured at the object boundary. The goal is to reconstruct internal optical properties of the body, such as absorption and diffusivity. In this work, it is assumed that the imaged object is composed of an approximately homogeneous background with clearly distinguishable embedded inhomogeneities. An algorithm for finding the maximum a posteriori estimate for the absorption and diffusion coefficients is introduced assuming an edge-preferring prior and an additive Gaussian measurement noise model. The method is based on iteratively combining a lagged diffusivity step and a linearization of the measurement model of diffuse optical tomography with priorconditioned LSQR. The performance of the reconstruction technique is tested via three-dimensional numerical experiments with simulated measurement data.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stochastic Galerkin finite element method with local conductivity basis for electrical impedance tomography

The objective of electrical impedance tomography is to deduce information about the conductivity inside a physical body from electrode measurements of current and voltage at the object boundary. In this work, the unknown conductivity is modeled as a random field parametrized by its values at a set of pixels. The uncertainty in the pixel values is propagated to the electrode measurements by numerically solving the forward problem of impedance tomography by a stochastic Galerkin finite element method in the framework of the complete electrode model. For a given set of electrode measurements, the stochastic forward solution is employed in approximately parametrizing the posterior probability density of the conductivity and contact resistances. Subsequently, the conductivity is reconstructed by computing the maximum a posteriori and conditional mean estimates as well as the posterior covariance. The functionality of this approach is demonstrated with experimental water tank data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Construction of invisible conductivity perturbations for the point electrode model in electrical impedance tomography

We explain how to build invisible isotropic conductivity perturbations of the unit conductivity in the framework of the point electrode model for two-dimensional electrical impedance tomography. The theoretical approach, based on solving a fixed point problem, is constructive and allows the implementation of an algorithm for approximating the invisible perturbations. The functionality of the method is demonstrated via numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Edge-enhancing reconstruction algorithm for three-dimensional electrical impedance tomography

Electrical impedance tomography is an imaging modality for extracting information on the conductivity distribution inside a physical body from boundary measurements of current and voltage. In many practical applications, it is a priori known that the conductivity consists of embedded inhomogeneities in an approximately constant background. This work introduces an iterative reconstruction algorithm that aims at finding the maximum a posteriori estimate for the conductivity assuming an edge-preferring prior. The method is based on applying (a single step of) priorconditioned lagged diffusivity iteration to sequential linearizations of the forward model. The algorithm is capable of producing reconstructions on dense unstructured three-dimensional finite element meshes and with a high number of measurement electrodes. The functionality of the proposed technique is demonstrated with both simulated and experimental data in the framework of the complete electrode model, which is the most accurate model for practical impedance tomography.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Optimizing electrode positions in electrical impedance tomography

This work considers finding optimal positions for the electrodes within the Bayesian paradigm based on available prior information on the conductivity; the aim is to place the electrodes so that the posterior density of the (discretized) conductivity, i.e., the conditional density of the conductivity given the measurements, is as localized as possible. To make such an approach computationally feasible, the complete electrode forward model of impedance tomography is linearized around the prior expectation of the conductivity, allowing explicit representation for the (approximate) posterior covariance matrix. Two approaches are considered: minimizing the trace or the determinant of the posterior covariance. The introduced optimization algorithm is of the steepest descent type, with the needed gradients computed based on appropriate Fréchet derivatives of the complete electrode model. The functionality of the methodology is demonstrated via two-dimensional numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A note on analyticity properties of far field patterns

In scattering theory the far field pattern describes the directional dependence of a time-harmonic wave scattered by an obstacle or inhomogeneous medium, when observed sufficiently far away from these objects. Considering plane wave excitations, the far field pattern can be written as a function of two variables, namely the direction of propagation of the incident plane wave and the observation direction, and it is well-known to be separately real analytic with respect to each of them. We show that the far field pattern is in fact a jointly real analytic function of these two variables.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Point measurements for a Neumann-to-Dirichlet map and the Calderón problem in the plane

This work considers properties of the Neumann-to-Dirichlet map for the conductivity equation under the assumption that the conductivity is identically one close to the boundary of the examined smooth, bounded and simply connected domain. It is demonstrated that the so-called bisweep data, i.e., the (relative) potential differences between two boundary points when delta currents of opposite signs are applied at the very same points, uniquely determine the whole Neumann-to-Dirichlet map. In two dimensions, the bisweep data extend as a holomorphic function of two variables to some (interior) neighborhood of the product boundary. It follows that the whole Neumann-to-Dirichlet map is characterized by the derivatives of the bisweep data at an arbitrary point. On the diagonal of the product boundary, these derivatives can be given with the help of the derivatives of the (relative) boundary potentials at some fixed point caused by the distributional current densities supported at the same point, and thus such point measurements uniquely define the Neumann-to-Dirichlet map. This observation also leads to a new, truly local uniqueness result for the so-called Calderón inverse conductivity problem.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simultaneous reconstruction of outer boundary shape and admittivity distribution in electrical impedance tomography

The aim of electrical impedance tomography is to reconstruct the admittivity distribution inside a physical body from boundary measurements of current and voltage. Due to the severe ill-posedness of the underlying inverse problem, the functionality of impedance tomography relies heavily on accurate modelling of the measurement geometry. In particular, almost all reconstruction algorithms require the precise shape of the imaged body as an input. In this work, the need for prior geometric information is relaxed by introducing a Newton-type output least squares algorithm that reconstructs the admittivity distribution and the object shape simultaneously. The method is built in the framework of the complete electrode model and it is based on the Fréchet derivative of the corresponding current-to-voltage map with respect to the object boundary shape. The functionality of the technique is demonstrated via numerical experiments with simulated measurement data.

0 citations0 reviewsNo code linkedOpen access