4.6 Article

Simultaneous Reconstruction of Outer Boundary Shape and Admittivity Distribution in Electrical Impedance Tomography

Journal

SIAM JOURNAL ON IMAGING SCIENCES
Volume 6, Issue 1, Pages 176-198

Publisher

SIAM PUBLICATIONS
DOI: 10.1137/120877301

Keywords

electrical impedance tomography; shape derivative; model inaccuracies; output least squares; complete electrode model; unknown boundary shape

Funding

  1. Finnish Funding Agency for Technology and Innovation TEKES [40370/08]
  2. Academy of Finland [135979, 141044]
  3. Academy of Finland (the Centre of Excellence in Inverse Problems Research and decision) [140280]
  4. Academy of Finland (AKA) [141044, 140280, 135979, 140280, 141044, 135979] Funding Source: Academy of Finland (AKA)

Ask authors/readers for more resources

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 is based on the Frechet 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.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.6
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

Article Mathematics, Applied

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

Henrik Garde, Nuutti Hyvonen

Summary: 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.

NUMERISCHE MATHEMATIK (2021)

Article Mathematics, Applied

Sequentially optimized projections in x-ray imaging *

M. Burger, A. Hauptmann, T. Helin, N. Hyvonen, J. P. Puska

Summary: This work applies Bayesian experimental design to select optimal projection geometries in (discretized) parallel beam x-ray tomography with Gaussian prior and additive noise. The introduced greedy exhaustive optimization algorithm proceeds sequentially, allowing redefining the region of interest after each projection and considering both A and D-optimality. Two-dimensional numerical experiments demonstrate the functionality of the approach.

INVERSE PROBLEMS (2021)

Article Automation & Control Systems

Uniform null controllability for parabolic equations with discontinuous diffusion coefficients*

Jeremi Darde, Sylvain Ervedoza, Roberto Morales

Summary: This article studies the null-controllability of a heat equation in a domain consisting of two media with different constant conductivities. The behavior of the system when the conductivity of the non-actuated medium approaches infinity is of particular interest, leading to a uniform null-controllability result under suitable geometric conditions. The strategy is based on analyzing the controllability of wave operators and using the transmutation technique to explain the geometric conditions.

ESAIM-CONTROL OPTIMISATION AND CALCULUS OF VARIATIONS (2021)

Article Mathematics, Applied

Approximation error method for imaging the human head by electrical impedance tomography*

V Candiani, N. Hyvonen, J. P. Kaipio, V Kolehmainen

Summary: This study focuses on electrical impedance tomography imaging of the human head for locating and classifying strokes in emergency care. By modeling natural variations in head geometry and skull, as well as the impact of misplaced electrodes, it enables reliable reconstruction of stroke-induced conductivity perturbations in an average head model. The functionality of a specific edge-preferring reconstruction algorithm for locating strokes is demonstrated through numerical experiments using simulated three-dimensional data.

INVERSE PROBLEMS (2021)

Article Mathematics, Applied

Mumford-Shah regularization in electrical impedance tomography with complete electrode model

Jyrki Jauhiainen, Aku Seppanen, Tuomo Valkonen

Summary: Electrical impedance tomography (EIT) is a method that aims to determine the conductivity within a target body by measuring electrical signals on its surface. The inverse conductivity problem is challenging due to limited measurement data, hence requiring regularization. Traditional regularization methods focus on promoting smooth features, but for targets consisting of multiple distinct objects or materials, the Mumford-Shah (M-S) regularization familiar in image segmentation is more suitable. However, it poses numerical challenges. The study shows through theoretical analysis that a modification of the Ambrosio-Tortorelli approximation of the M-S regularizer is applicable to EIT, specifically in the complete electrode model of boundary measurements. Numerical and experimental studies confirm the practicality of this approach, producing higher quality results compared to typical regularizations employed in EIT when the conductivity of the target consists of distinct smoothly-varying regions.

INVERSE PROBLEMS (2022)

Article Mathematics, Applied

SERIES REVERSION IN CALDER ′ ON'S PROBLEM

Henrik Garde, Nuutti Hyvonen

Summary: This work presents numerical methods for solving Calderon's problem using series reversions. By reversing the series of the forward map, a family of methods for solving the inverse problem is obtained, and the convergence of these methods is proven. The introduced numerical methods have the same computational complexity as solving the linearized inverse problem.

MATHEMATICS OF COMPUTATION (2022)

Article Mathematics, Applied

CONTACT ADAPTING ELECTRODE MODEL FOR ELECTRICAL IMPEDANCE TOMOGRAPHY

J. Darde, N. Hyvonen, T. Kuutela, T. Valkonens

Summary: This study proposes a new robust modeling method for contact electrodes in electrical impedance tomography. By assuming approximate knowledge about the electrodes' whereabouts and using a boundary admittivity function to determine their actual locations, the proposed method enables simultaneous reconstruction of the positions and strengths of the contacts.

SIAM JOURNAL ON APPLIED MATHEMATICS (2022)

Article Mathematics, Applied

EDGE-PROMOTING ADAPTIVE BAYESIAN EXPERIMENTAL DESIGN FOR X-RAY IMAGING

Tapio Helin, Nuutti Hyvonen, Juha-Pekka Puska

Summary: This work focuses on sequential edge-promoting Bayesian experimental design for linear inverse problems, specifically X-ray tomography. It interprets the computation process of total variation-type absorption reconstruction inside the imaged body using lagged diffusivity iteration in the Bayesian framework. By assuming a Gaussian additive noise model, an approximate Gaussian posterior with a covariance structure containing information about the location of edges in the posterior mean is obtained. The next projection geometry is then determined using A- or D-optimal Bayesian design, which minimizes the trace or determinant of the updated posterior covariance matrix that accounts for the new projection. Two- and three-dimensional numerical examples based on simulated data demonstrate the effectiveness of this approach.

SIAM JOURNAL ON SCIENTIFIC COMPUTING (2022)

Article Mathematics, Applied

The Morozov's principle applied to data assimilation problems

Laurent Bourgeois, Jeremi Darde

Summary: This paper focuses on the application of Morozov's principle to an abstract data assimilation framework and illustrates how to regularize ill-posed problems through three examples. It extends known results related to Morozov's choice of the regularization parameter to the case where the operator does not have a dense range, and computes the solution satisfying Morozov's principle using optimization duality.

ESAIM-MATHEMATICAL MODELLING AND NUMERICAL ANALYSIS (2022)

Article Mathematics, Applied

Series reversion for electrical impedance tomography with modeling errors *

H. Garde, N. Hyvonen, T. Kuutela

Summary: This work extends the series reversion method for Calderon's problem to realistic electrode measurements, treating both the internal admittivity of the investigated body and the contact admittivity at the electrode-object interfaces as unknowns. A family of numerical methods for solving the inverse problem of electrical impedance tomography is obtained through Taylor series reversion, with different parametrizations for the unknown internal and boundary admittivities. The functionality and convergence of the methods are established under certain conditions, and the methods are heuristically extended to more general settings using regularization motivated by Bayesian inversion. The performance of the regularized approach is tested through three-dimensional numerical examples.

INVERSE PROBLEMS (2023)

Article Engineering, Biomedical

Effects of atlas-based anatomy on modelled light transport in the neonatal head

Pauliina Hirvi, Topi Kuutela, Qianqian Fang, Antti Hannukainen, Nuutti Hyvonen, Ilkka Nissila

Summary: This study investigates the approximation error when using an atlas instead of the neonate's own anatomical model in diffuse optical tomography (DOT). It is found that there is a considerable approximation error, but the method using atlas-based imaging can still detect small changes in brain absorption. Individual-level atlas models can be used in DOT difference imaging when an exactly age-appropriate atlas is not available.

PHYSICS IN MEDICINE AND BIOLOGY (2023)

Article Mathematics, Applied

SERIES REVERSION IN CALDERON'S PROBLEM

Henrik Garde, Nuutti Hyvonen

Summary: This study presents explicit series reversions for the solution of Calderon's problem and derives a family of numerical methods to improve the accuracy of the solution by reversing the Taylor series of the forward map. The convergence of these numerical methods is shown under conditions that ensure the invertibility of the Frechet derivative of the forward map. The introduced numerical methods have the same computational complexity as solving the linearized inverse problem.

MATHEMATICS OF COMPUTATION (2022)

Article Mathematics, Applied

RECONSTRUCTION OF SINGULAR AND DEGENERATE INCLUSIONS IN CALDERON'S PROBLEM

Henrik Garde, Nuutti Hyvonen

Summary: We investigate the reconstruction of the support of an unknown perturbation to a known conductivity coefficient in Calderon's problem. We generalize a previous result by allowing the unknown coefficient to be the restriction of an A(2)-Muckenhoupt weight in parts of the domain, incorporating singular and degenerate behavior. Our main constructive result characterizes the outer shape of the support based on a local Neumann-to-Dirichlet map defined on an open subset of the domain boundary.

INVERSE PROBLEMS AND IMAGING (2022)

Article Mathematics, Applied

INVERSE HEAT SOURCE PROBLEM AND EXPERIMENTAL DESIGN FOR DETERMINING IRON LOSS DISTRIBUTION

Antti Hannukainen, Nuutti Hyvonen, Lauri Perkkio

Summary: Iron loss determination in the magnetic core of electrical machines is treated as an inverse heat source problem, with sensor positions optimized for minimizing uncertainty in reconstruction. The study focuses on formulating the problem and efficiently solving the discretized sensor optimization and source reconstruction problems. A semi-realistic linear model is discretized using finite elements and numerically studied.

SIAM JOURNAL ON SCIENTIFIC COMPUTING (2021)

No Data Available