Reconstructing conductivities with boundary corrected D-bar method
|
|
- Hillary York
- 5 years ago
- Views:
Transcription
1 Reconstructing conductivities with boundary corrected D-bar method Janne Tamminen June 24, 2011
2 Short introduction to EIT The Boundary correction procedure The D-bar method Simulation of measurement data, numerical D-bar
3 Short introduction to EIT The aim of electrical impedance tomography (eit) is to form an image of the conductivity distribution inside an unknown body using electric boundary measurements. Applications in medical imaging, nondestructive testing, subsurface monitoring: monitoring heart and lungs of unconscious patients, detecting pulmonary edema, breast cancer detection, detecting cracks in concrete structures, environmental applications...
4 Short introduction to EIT The inverse conductivity problem introduced by Calderón is to find the conductivity σ when the Dirichlet-to-Neumann -map (dn-map) Λ σ is known. { (σ u) = 0 in Ω1, (1) u = f on Ω 1. Λ σ f = σ u ν Ω 1 (2) The problem of reconstructing σ from Λ σ is nonlinear and ill-posed.
5 Short introduction to EIT In the two-dimensional case, unique reconstruction can be obtained: For σ W 2,p (Ω 1 ),p > 1, Nachman 1996 (the D-bar method), Knudsen-Lassas-Mueller-Siltanen 2009 For σ W 1,p (Ω 1 ),p > 2, Brown-Uhlmann 1997 and for σ L (Ω 1 ) Astala-Päivärinta Nachman s method reduces the reconstruction problem to the case σ = 1 near the boundary. This procedure is tested numerically in this joint work with S. Siltanen.
6 The Boundary correction procedure Short introduction to EIT The Boundary correction procedure The D-bar method Simulation of measurement data, numerical D-bar
7 The Boundary correction procedure We transform the original conductivity equation to Schrödinger equation: by writing u = σ 1/2 ũ we get ( q)ũ = 0 in Ω 1, (3) where q = σ 1/2 (σ 1/2 ). This means we have to have σ c 0 > 0. The dn-map becomes Λ q = σ 1/2 (Λ σ + 1 σ 2 ν )σ 1/2. (4) In order for this transformation to be useful we need to have σ 1 near Ω 1 so that Λ q = Λ σ.
8 The Boundary correction procedure Ω 2 Ω 1 Λ σ Λ γ? We extend σ: { σ(x), x Ω 1, γ(x) = σ(x), x Ω 2 \Ω 1, σ σ W 2,p σ 1 σ Ω1 = σ Ω1 (5) σ ν Ω 1 = σ ν Ω 1 (6) giving us γ W 2,p (Ω 2 ) and Λ q = Λ γ.
9 The Boundary correction procedure In this work, Ω 1 = D(0,r 1 ), r 1 = 1.0 Ω 2 = D(0,r 2 ), r 2 = 1.2 σ L (Ω 1 ) σ W 2,p (Ω 2 \Ω 1 ) γ L (Ω 2 ),
10 The Boundary correction procedure Define two Dirichlet problems, for j = 1,2, ( σ u j ) = 0 in Ω 2 \Ω 1 u j = f j on Ω j u j = 0 on Ω i, i = 1,2, i j. (7) Four new dn maps in Ω 2 \Ω 1 can be characterized by so λ ij : H 1/2 ( Ω j ) H 1/2 ( Ω i ). Λ ij f j = σ u j ν Ω i, i,j = 1,2 (8)
11 The Boundary correction procedure Proposition (6.1 in Nachman 1996) Let γ W 2,p (Ω 2 ),p > 1, then Λ γ = Λ 22 +Λ 21 (Λ σ Λ 11 ) 1 Λ 12. (9) The proof consists of showing that the operator Λ σ Λ 11 is invertible, and the identities (Λ γ Λ 22 )f 2 = Λ 21 (u Ω1 ) (10) (Λ σ Λ 11 )(u Ω1 ) = Λ 12 f 2 (11) for any u that solves the conductivity equation in Ω 2 with u Ω2 = f 2.
12 The D-bar method Short introduction to EIT The Boundary correction procedure The D-bar method Simulation of measurement data, numerical D-bar
13 The D-bar method The D-bar method: based on Nachman 1996 robust algorithm was given by Siltanen-Mueller-Isaacson The method has been successfully tested on a chest phantom and on in vivo human chest data (Isaacson-Mueller-Newell-Siltanen). Λ γ t(k) γ
14 The D-bar method The CGO-solution (Complex Geometric Optics) on Ω 2 can be solved from ψ(,k) Ω2 = e ikx S k (Λ γ Λ 1 )ψ(,k) Ω2, (12) in the Sobolev space H 1/2 ( Ω 2 ) for all k C\{0}. Here S k is a single-layer operator (S k φ)(x) := G k (x y)φ(y)ds, Ω 2 and G k is Faddeev s Green function defined by G k (x) := e ikx g k (x), g k (x) := 1 R2 e ix ξ (2π) 2 ξ 2 +2k(ξ 1 +iξ 2 ) dξ.
15 The D-bar method Using the CGO-solution we can define the scattering transform: t(k) = e i k x (Λ γ Λ 1 )ψ(,k)ds. (13) Ω 2 For each fixed x Ω 2, we would solve the following integral formulation of the D-bar equation: µ(x,k) = 1+ 1 t(k ) x+k x) (2π) 2 µ(x,k )dk R 2 (k k ) k ei(k 1 dk 2, (14) where µ(x,k) = exp( ik(x 1 +ix 2 ))ψ(x,k). Then the conductivity would be perfectly reconstructed as γ(x) = µ(x,0) 2.
16 Simulation of measurement data, numerical D-bar Short introduction to EIT The Boundary correction procedure The D-bar method Simulation of measurement data, numerical D-bar
17 Simulation of measurement data, numerical D-bar In this work we omit the requirement of continuity of the derivative γ ν and use the method of Nakamura-Tanuma-Siltanen-Wang to calculate g(θ) σ Ω1. We extend σ to γ using the following extension in polar coordinates: σ(ρ,θ), ρ r 1, γ(ρ,θ) = (g(θ) 1)f m (ρ)+1, r 1 < ρ r e, 1, r e < ρ r 2, (15) where r e = and f m (ρ) 0 is a suitable third-degree polynomial satisfying f m (r 1 ) = 1 and f m (r e ) = 0.
18 Simulation of measurement data, numerical D-bar Any linear operator Λ : H 1/2 ( Ω i ) H 1/2 ( Ω j ) can be represented by a matrix in the following way. Define a truncated orthonormal basis at the boundary Ω k : φ (n) k (θ) = 1 2πrk e inθ, n = N,...,N, k = 1,2. (16) Write any function f : Ω i C as a vector f = [ˆf( N),ˆf( N+1),...,ˆf(N 1),ˆf(N)] T, ˆf(n) = Ω i fφ (n) i ds. Then the operator Λ is approximated by the matrix L = [û(l,n)], where û(l,n) = (Λφ (n) i )φ (l) j ds. (17) Ω j
19 Simulation of measurement data, numerical D-bar In practice in eit we measure the Neumann-to-Dirichlet -map approximated by the matrix R σ. We simulate measurement noise by using the matrix R ε σ := R σ +ce, (18) where E is a matrix with random entries independently distributed according to the Gaussian normal density. Then, the noisy DN-matrix L ε σ is roughly speaking the inverse of Rε σ, and the boundary correction procedure gives us L ε γ = L22 +L 21 (L ε σ L11 ) 1 L 12, (19) provided that the matrix L ε σ L 11 is invertible.
20 Simulation of measurement data, numerical D-bar The D-bar method numerically: expand e ikx Ω2 as a vector g in our finite trigonometric basis, solve the CGO-solution ψ k := [I +S k (L ε γ L 1 )] 1 g. (20) for k ranging in a grid inside the disc k < R ( the truncation radius R > 0). Define the truncated scattering transform by { t R (k) = Ω 2 e i k x F 1 ((L ε γ L 1) ψ k )(x)ds for k < R, (21) 0, otherwise, Finally solve the equation (14) with the numerical algorithm of Knudsen-Mueller-Siltanen using t R and denote the solution by µ R (x,k). Then γ(x) µ R (x,0) 2.
21 Short introduction to EIT The Boundary correction procedure The D-bar method Simulation of measurement data, numerical D-bar
22 Solutions to the Dirichlet problems are calculated using FEM, with triangles in the disk Ω 1 and triangles in the annulus Ω 2 \Ω 1. R ǫ 1 R 1 / R 1 = (the ACT3 impedance tomography imager of Rensselaer Polytechnic Institute has SNR of 95.5 db noise of %) ǫ fem = R th 1 R 1 / R th 1 = , the error R ǫ σ R σ / R σ ranges between and , The condition number of the matrix L ε σ L 11 was less than 27 in all our test cases. The error L ǫ γ L 2 γ / L 2 γ, where L 2 γ is the dn map calculated directly on the boundary Ω 2, was less than 2.2% in all cases.
23 % L ǫ γ=1 L2 γ=1 / L2 γ= Lǫ γ=1 L γ=1 / L γ= r 2
24 Example 1 σ Ω π 2π θ over boundary
25 Example 2 6 σ Ω π 2π θ over boundary
26 Example 3 7 σ Ω π 2π θ over boundary
27 Example 4 σ Ω π 2π θ over boundary
28 Re(t(k)) Im(t(k)) non-noisy L σ 10 noisy L ε σ Scattering transform, real and imaginary parts, of Example noisy and corrected L ε γ
29 R = Original Uncorrected Corrected
30 % 50 Example 1 50 Example % Example 3 Example R R
31 R = 3.0 Uncorrected R = 5.0 Uncorrected R = 6.0 Uncorrected Example 1 Original conductivity 30% 25% 37% R = 3.0 Corrected R = 5.0 Corrected R = 6.0 Corrected 26% 18% 49%
32 R = 3.0 Uncorrected R = 5.4 Uncorrected R = 6.0 Uncorrected Example 2 Original conductivity 42% 29% 35% R = 3.0 Corrected R = 5.4 Corrected R = 6.0 Corrected 35% 15% 39%
33 R = 3.0 Uncorrected R = 5.0 Uncorrected R = 6.0 Uncorrected Example 3 Original conductivity 66% 60% 67% R = 3.0 Corrected R = 5.0 Corrected R = 6.0 Corrected 59% 39% 75%
34 R = 3.0 Uncorrected R = 4.8 Uncorrected R = 6.0 Uncorrected Example 4 Original conductivity 25% 22% 49% R = 3.0 Corrected R = 4.8 Corrected R = 6.0 Corrected 25% 21% 63%
35 Thank you!
36 K. Astala and L. Päivärinta, Calderón s inverse conductivity problem in the plane, Ann. of Math., 163 (2006), R. M. Brown and G. Uhlmann, Uniqueness in the inverse conductivity problem for nonsmooth conductivities in two dimensions, Comm. Partial Differential Equations, 22 (1997), A. P. Calderón, On an inverse boundary value problem Seminar on Numerical Analysis and its Applications to Continuum Physics, Soc. Brasileira de Matemàtica, 1980, R. D. Cook, G. J. Saulnier and J. C. Goble, A phase sensitive voltmeter for a high-speed, high-precision electrical impedance tomograph, Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Vol.13, No.1 (1991), D. Isaacson, J. L. Mueller, J. C. Newell and S. Siltanen, Reconstructions of chest phantoms by the D-bar method for electrical impedance tomography, IEEE Trans. Med. Im., 23 (2004),
37 D. Isaacson, J. L. Mueller, J. C. Newell and S. Siltanen, Imaging cardiac activity by the D-bar method for electrical impedance tomography, Physiol. Meas., 27 (2006), S43 S50. K. Knudsen, On the Inverse Conductivity Problem, Ph.D. thesis, Department of Mathematical Sciences, Aalborg University, Denmark (2002) K. Knudsen, J. L. Mueller and S. Siltanen, Numerical solution method for the D-bar-equation in the plane, J. Comp. Phys., 198 (2004), K. Knudsen, M. J. Lassas, J. L. Mueller and S. Siltanen, D-bar method for electrical impedance tomography with discontinuous conductivities SIAM J. Appl. Math., 67 (2007), K. Knudsen, M. Lassas, J. L. Mueller and S. Siltanen, Reconstructions of piecewise constant conductivities by the D-bar method for Electrical Impedance Tomography, Journal of Physics: Conference Series, 124 (2008),. K. Knudsen, M. Lassas, J. L. Mueller and S. Siltanen, Regularized D-bar method for the inverse conductivity problem, Inverse Problems and Imaging, 3 (2009),
38 J. L. Mueller and S. Siltanen, Direct reconstructions of conductivities from boundary measurements, SIAM J. Sci. Comp., 24 (2003), A. I. Nachman, Reconstructions from boundary measurements, Ann. of Math., 128 (1988), A. I. Nachman, Global uniqueness for a two-dimensional inverse boundary value problem, Ann. of Math., 143 (1996), G. Nakamura, K. Tanuma, S. Siltanen and S. Wang, Numerical recovery of conductivity at the boundary from the localized Dirichlet to Neumann map, Computing, 75 (2005), S. Siltanen, J. Mueller and D. Isaacson, An implementation of the reconstruction algorithm of A. Nachman for the 2-D inverse conductivity problem, Inverse Problems, 16 (2000), J. Sylvester and G. Uhlmann, A global uniqueness theorem for an inverse boundary value problem, Ann. of Math, 125 (1987),
Reconstructing conductivities with boundary corrected D-bar method
J. Inverse Ill-Posed Probl., Ahead of Print DOI 10.1515/jip-2013-0042 de Gruyter 2014 Reconstructing conductivities with boundary corrected D-bar method Samuli Siltanen and Janne P. Tamminen Abstract.
More informationComputational inversion based on complex geometrical optics solutions
Computational inversion based on complex geometrical optics solutions Samuli Siltanen Department of Mathematics and Statistics University of Helsinki, Finland samuli.siltanen@helsinki.fi http://www.siltanen-research.net
More informationNumerical computation of complex geometrical optics solutions to the conductivity equation
Numerical computation of complex geometrical optics solutions to the conductivity equation Kari Astala a, Jennifer L Mueller b, Lassi Päivärinta a, Samuli Siltanen a a Department of Mathematics and Statistics,
More informationA CALDERÓN PROBLEM WITH FREQUENCY-DIFFERENTIAL DATA IN DISPERSIVE MEDIA. has a unique solution. The Dirichlet-to-Neumann map
A CALDERÓN PROBLEM WITH FREQUENCY-DIFFERENTIAL DATA IN DISPERSIVE MEDIA SUNGWHAN KIM AND ALEXANDRU TAMASAN ABSTRACT. We consider the problem of identifying a complex valued coefficient γ(x, ω) in the conductivity
More informationComplex geometrical optics solutions for Lipschitz conductivities
Rev. Mat. Iberoamericana 19 (2003), 57 72 Complex geometrical optics solutions for Lipschitz conductivities Lassi Päivärinta, Alexander Panchenko and Gunther Uhlmann Abstract We prove the existence of
More informationCalderón s inverse problem in 2D Electrical Impedance Tomography
Calderón s inverse problem in 2D Electrical Impedance Tomography Kari Astala (University of Helsinki) Joint work with: Matti Lassas, Lassi Päivärinta, Samuli Siltanen, Jennifer Mueller and Alan Perämäki
More informationDIRECT ELECTRICAL IMPEDANCE TOMOGRAPHY FOR NONSMOOTH CONDUCTIVITIES
DIRECT ELECTRICAL IMPEDANCE TOMOGRAPHY FOR NONSMOOTH CONDUCTIVITIES K ASTALA, J L MUELLER, A PERÄMÄKI, L PÄIVÄRINTA AND S SILTANEN Abstract. A new reconstruction algorithm is presented for eit in dimension
More informationContributors and Resources
for Introduction to Numerical Methods for d-bar Problems Jennifer Mueller, presented by Andreas Stahel Colorado State University Bern University of Applied Sciences, Switzerland Lexington, Kentucky, May
More informationAALBORG UNIVERSITY. A new direct method for reconstructing isotropic conductivities in the plane. Department of Mathematical Sciences.
AALBORG UNIVERSITY A new direct method for reconstructing isotropic conductivities in the plane by Kim Knudsen January 23 R-23-1 Department of Mathematical Sciences Aalborg University Fredrik Bajers Vej
More informationNONLINEAR FOURIER ANALYSIS FOR DISCONTINUOUS CONDUCTIVITIES: COMPUTATIONAL RESULTS
NONLINEAR FOURIER ANALYSIS FOR DISCONTINUOUS CONDUCTIVITIES: COMPUTATIONAL RESULTS K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN Contents. Introduction. The nonlinear Fourier transform 6 3. Inverse
More informationElectrical Impedance Tomography: 3D Reconstructions using Scattering Transforms
Downloaded from orbit.dtu.dk on: Nov 11, 2018 Electrical Impedance Tomography: 3D Reconstructions using Scattering Transforms Delbary, Fabrice; Hansen, Per Christian; Knudsen, Kim Published in: Applicable
More informationInverse Electromagnetic Problems
Title: Name: Affil./Addr. 1: Inverse Electromagnetic Problems Gunther Uhlmann, Ting Zhou University of Washington and University of California Irvine /gunther@math.washington.edu Affil./Addr. 2: 340 Rowland
More informationMapping Properties of the Nonlinear Fourier Transform in Dimension Two
Communications in Partial Differential Equations, 3: 59 60, 007 Copyright Taylor & Francis Group, LLC ISSN 0360-530 print/53-433 online DOI: 0.080/036053005005304 Mapping Properties of the Nonlinear Fourier
More informationComputational inversion using complex geometric optics solutions and nonlinear Fourier transforms
Computational inversion using complex geometric optics solutions and nonlinear Fourier transforms Samuli Siltanen Department of Mathematics and Statistics University of Helsinki, Finland samuli.siltanen@helsinki.fi
More informationarxiv: v2 [math.ap] 17 Jun 2014
NONLINEAR FOURIER ANALYSIS FOR DISCONTINUOUS CONDUCTIVITIES: COMPUTATIONAL RESULTS K. ASTALA, L. PÄIVÄRINTA, J. M. REYES AND S. SILTANEN arxiv:30.705v [math.ap] 7 Jun 04 Contents. Introduction. The nonlinear
More informationFast shape-reconstruction in electrical impedance tomography
Fast shape-reconstruction in electrical impedance tomography Bastian von Harrach bastian.harrach@uni-wuerzburg.de (joint work with Marcel Ullrich) Institut für Mathematik - IX, Universität Würzburg The
More informationAALBORG UNIVERSITY. The Calderón problem with partial data for less smooth conductivities. Kim Knudsen. Department of Mathematical Sciences
AALBORG UNIVERSITY The Calderón problem with partial data for less smooth conductivities by Kim Knudsen R-2005-07 February 2005 Department of Mathematical Sciences Aalborg University Fredrik Bajers Vej
More informationis measured (n is the unit outward normal) at the boundary, to define the Dirichlet-
SIAM J MATH ANAL Vol 45, No 5, pp 700 709 c 013 Society for Industrial and Applied Mathematics ON A CALDERÓN PROBLEM IN FREQUENCY DIFFERENTIAL ELECTRICAL IMPEDANCE TOMOGRAPHY SUNGWHAN KIM AND ALEXANDRU
More informationSome issues on Electrical Impedance Tomography with complex coefficient
Some issues on Electrical Impedance Tomography with complex coefficient Elisa Francini (Università di Firenze) in collaboration with Elena Beretta and Sergio Vessella E. Francini (Università di Firenze)
More informationStatistical and Computational Inverse Problems with Applications Part 2: Introduction to inverse problems and example applications
Statistical and Computational Inverse Problems with Applications Part 2: Introduction to inverse problems and example applications Aku Seppänen Inverse Problems Group Department of Applied Physics University
More informationDownloaded 11/20/13 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 33, No. 1, pp. 328 341 c 2011 Society for Industrial and Applied Mathematics A SIMPLE NUMERICAL METHOD FOR COMPLEX GEOMETRICAL OPTICS SOLUTIONS TO THE CONDUCTIVITY EQUATION KUI
More informationA Direct Method for reconstructing inclusions from Electrostatic Data
A Direct Method for reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with:
More informationarxiv: v1 [math.ap] 25 Sep 2018
arxiv:1809.09272v1 [math.ap] 25 Sep 2018 NACHMAN S RECONSTRUCTION METHO FOR THE CALERÓN PROBLEM WITH ISCONTINUOUS CONUCTIVITIES GEORGE LYTLE, PETER PERRY, AN SAMULI SILTANEN Abstract. We show that Nachman
More informationWhat are inverse problems?
What are inverse problems? 6th Forum of Polish Mathematicians Warsaw University of Technology Joonas Ilmavirta University of Jyväskylä 8 September 2015 Joonas Ilmavirta (Jyväskylä) Inverse problems 8 September
More informationarxiv: v1 [math.na] 19 Sep 2018
SPECTRAL APPROACHES TO D-BAR PROBLEMS ARISING IN ELECTRICAL IMPEDANCE TOMOGRAPHY arxiv:89.8294v [math.na] 9 Sep 28 CHRISTIAN KLEIN AND NIKOLA STOILOV Abstract. In this work we present spectral approaches
More informationProceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 SOME INVERSE SCATTERING PROBLEMS FOR TWO-DIMENSIONAL SCHRÖDINGER
More informationCOMPLEX SPHERICAL WAVES AND INVERSE PROBLEMS IN UNBOUNDED DOMAINS
COMPLEX SPHERICAL WAVES AND INVERSE PROBLEMS IN UNBOUNDED DOMAINS MIKKO SALO AND JENN-NAN WANG Abstract. This work is motivated by the inverse conductivity problem of identifying an embedded object in
More informationarxiv: v2 [math.na] 2 Feb 2016 J. P. TAMMINEN, T. TARVAINEN AND S. SILTANEN
THE D-BAR METHOD FOR DIFFUSE OPTICAL TOMOGRAPHY: A COMPUTATIONAL STUDY arxiv:6.4636v2 [math.na] 2 Feb 26 J. P. TAMMINEN, T. TARVAINEN AND S. SILTANEN Abstract. The D-bar method at negative energy is numerically
More informationarxiv: v1 [math.ap] 6 Mar 2017
A Fidelity-embedded Regularization Method for Robust Electrical Impedance Tomography Kyounghun Lee, Eung Je Woo, and Jin Keun Seo arxiv:1703.01821v1 [math.ap] 6 Mar 2017 March 7, 2017 Abstract Electrical
More informationOn uniqueness in the inverse conductivity problem with local data
On uniqueness in the inverse conductivity problem with local data Victor Isakov June 21, 2006 1 Introduction The inverse condictivity problem with many boundary measurements consists of recovery of conductivity
More informationPut Paper Number Here
Proceedings of Proceedings of icipe : rd International Conference on Inverse Problems in Engineering June,, Port Ludlow, WA, USA Put Paper Number Here NEW RECONSTRUCTION ALGORITHMS IN ELECTRICAL IMPEDANCE
More informationDETERMINING A FIRST ORDER PERTURBATION OF THE BIHARMONIC OPERATOR BY PARTIAL BOUNDARY MEASUREMENTS
DETERMINING A FIRST ORDER PERTURBATION OF THE BIHARMONIC OPERATOR BY PARTIAL BOUNDARY MEASUREMENTS KATSIARYNA KRUPCHYK, MATTI LASSAS, AND GUNTHER UHLMANN Abstract. We consider an operator 2 + A(x) D +
More informationStability and instability in inverse problems
Stability and instability in inverse problems Mikhail I. Isaev supervisor: Roman G. Novikov Centre de Mathématiques Appliquées, École Polytechnique November 27, 2013. Plan of the presentation The Gel fand
More informationEstimation of transmission eigenvalues and the index of refraction from Cauchy data
Estimation of transmission eigenvalues and the index of refraction from Cauchy data Jiguang Sun Abstract Recently the transmission eigenvalue problem has come to play an important role and received a lot
More informationElectrostatic Imaging via Conformal Mapping. R. Kress. joint work with I. Akduman, Istanbul and H. Haddar, Paris
Electrostatic Imaging via Conformal Mapping R. Kress Göttingen joint work with I. Akduman, Istanbul and H. Haddar, Paris Or: A new solution method for inverse boundary value problems for the Laplace equation
More informationA uniqueness result and image reconstruction of the orthotropic conductivity in magnetic resonance electrical impedance tomography
c de Gruyter 28 J. Inv. Ill-Posed Problems 16 28), 381 396 DOI 1.1515 / JIIP.28.21 A uniqueness result and image reconstruction of the orthotropic conductivity in magnetic resonance electrical impedance
More informationarxiv: v2 [math.na] 26 Mar 2018
APPROXIMATION OF FULL-BOUNDARY DATA FROM PARTIAL-BOUNDARY ELECTRODE MEASUREMENTS A. HAUPTMANN arxiv:1703.05550v2 [math.na] 26 Mar 2018 Abstract. Measurements on a subset of the boundary are common in electrical
More informationInverse conductivity problem for inaccurate measurements
Inverse Problems 12 (1996) 869 883. Printed in the UK Inverse conductivity problem for inaccurate measurements Elena Cherkaeva and Alan C Tripp The University of Utah, Department of Geophysics, 717 Browning
More informationInverse problems and medical imaging
Inverse problems and medical imaging Bastian von Harrach harrach@math.uni-frankfurt.de Institute of Mathematics, Goethe University Frankfurt, Germany Colloquium of the Department of Mathematics Saarland
More informationINVERSE BOUNDARY VALUE PROBLEMS FOR THE PERTURBED POLYHARMONIC OPERATOR
INVERSE BOUNDARY VALUE PROBLEMS FOR THE PERTURBED POLYHARMONIC OPERATOR KATSIARYNA KRUPCHYK, MATTI LASSAS, AND GUNTHER UHLMANN Abstract. We show that a first order perturbation A(x) D + q(x) of the polyharmonic
More informationDetecting stochastic inclusions in electrical impedance tomography
Detecting stochastic inclusions in electrical impedance tomography Bastian von Harrach harrach@math.uni-frankfurt.de (joint work with A. Barth, N. Hyvönen and L. Mustonen) Institute of Mathematics, Goethe
More informationUniqueness and stability of Lamé parameters in elastography
Uniqueness stability of Lamé parameters in elastography Ru-Yu Lai Abstract This paper concerns an hybrid inverse problem involving elastic measurements called Transient Elastography TE which enables detection
More informationarxiv: v1 [math.ap] 11 Jun 2007
Inverse Conductivity Problem for a Parabolic Equation using a Carleman Estimate with One Observation arxiv:0706.1422v1 [math.ap 11 Jun 2007 November 15, 2018 Patricia Gaitan Laboratoire d Analyse, Topologie,
More informationOn Transformation Optics based Cloaking
On Transformation Optics based Cloaking TING ZHOU Northeastern University CIRM Outline Inverse Problem and Invisibility 1 Inverse Problem and Invisibility 2 Singular cloaking medium 3 Regularization and
More informationA note on the MUSIC algorithm for impedance tomography
A note on the MUSIC algorithm for impedance tomography Martin Hanke Institut für Mathematik, Johannes Gutenberg-Universität, 55099 Mainz, Germany E-mail: hanke@math.uni-mainz.de Abstract. We investigate
More informationAllan Greenleaf, Matti Lassas, and Gunther Uhlmann
Mathematical Research Letters 10, 685 693 (2003) ON NONUNIQUENESS FOR CALDERÓN S INVERSE PROBLEM Allan Greenleaf, Matti Lassas, and Gunther Uhlmann Abstract. We construct anisotropic conductivities with
More informationINVERSE BOUNDARY VALUE PROBLEMS FOR THE MAGNETIC SCHRÖDINGER EQUATION
INVERSE BOUNDARY VALUE PROBLEMS FOR THE MAGNETIC SCHRÖDINGER EQUATION MIKKO SALO Abstract. We survey recent results on inverse boundary value problems for the magnetic Schrödinger equation. 1. Introduction
More informationINCREASING STABILITY IN THE INVERSE PROBLEM FOR THE SCHRÖDINGER EQUATION. A Dissertation by. Li Liang
INCREASING STABILITY IN THE INVERSE PROBLEM FOR THE SCHRÖDINGER EQUATION A Dissertation by Li Liang Bachelor of Science, Wichita State University, 20 Master of Science, Wichita State University, 203 Submitted
More informationInverse problems and medical imaging
Inverse problems and medical imaging Bastian von Harrach harrach@math.uni-stuttgart.de Chair of Optimization and Inverse Problems, University of Stuttgart, Germany Rhein-Main Arbeitskreis Mathematics of
More informationRegularizing the optimization-based solution of inverse coefficient problems with monotonicity constraints
Regularizing the optimization-based solution of inverse coefficient problems with monotonicity constraints Bastian von Harrach (joint work with Mach Nguyet Minh) http://numerical.solutions Institute of
More informationResistor Networks and Optimal Grids for Electrical Impedance Tomography with Partial Boundary Measurements
Resistor Networks and Optimal Grids for Electrical Impedance Tomography with Partial Boundary Measurements Alexander Mamonov 1, Liliana Borcea 2, Vladimir Druskin 3, Fernando Guevara Vasquez 4 1 Institute
More informationarxiv: v1 [math.ap] 21 Dec 2014
CONSTRUCTION OF INVISIBLE CONUCTIVITY PERTURBATIONS FOR THE POINT ELECTROE MOEL IN ELECTRICAL IMPEANCE TOMOGRAPHY LUCAS CHESNEL, NUUTTI HYVÖNEN, AN STRATOS STABOULIS Abstract. We explain how to build invisible
More informationMonotonicity arguments in electrical impedance tomography
Monotonicity arguments in electrical impedance tomography Bastian Gebauer gebauer@math.uni-mainz.de Institut für Mathematik, Joh. Gutenberg-Universität Mainz, Germany NAM-Kolloquium, Georg-August-Universität
More informationThe Mathematics of Invisibility: Cloaking Devices, Electromagnetic Wormholes, and Inverse Problems. Lectures 1-2
CSU Graduate Workshop on Inverse Problems July 30, 2007 The Mathematics of Invisibility: Cloaking Devices, Electromagnetic Wormholes, and Inverse Problems Lectures 1-2 Allan Greenleaf Department of Mathematics,
More informationIncreasing stability in an inverse problem for the acoustic equation
Increasing stability in an inverse problem for the acoustic equation Sei Nagayasu Gunther Uhlmann Jenn-Nan Wang Abstract In this work we study the inverse boundary value problem of determining the refractive
More informationUniversity of North Carolina at Charlotte Charlotte, USA Norwegian Institute of Science and Technology Trondheim, Norway
1 A globally convergent numerical method for some coefficient inverse problems Michael V. Klibanov, Larisa Beilina University of North Carolina at Charlotte Charlotte, USA Norwegian Institute of Science
More informationNumerical Harmonic Analysis on the Hyperbolic Plane
Numerical Harmonic Analysis on the Hyperbolic Plane Buma Fridman, Peter Kuchment, Kirk Lancaster, Serguei Lissianoi, Mila Mogilevsky, Daowei Ma, Igor Ponomarev, and Vassilis Papanicolaou Mathematics and
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes Fachrichtung 6.1 Mathematik Preprint Nr. 362 An iterative method for EIT involving only solutions of Poisson equations. I: Mesh-free forward solver Thorsten Hohage and Sergej
More informationThe Factorization Method for the Reconstruction of Inclusions
The Factorization Method for the Reconstruction of Inclusions Martin Hanke Institut für Mathematik Johannes Gutenberg-Universität Mainz hanke@math.uni-mainz.de January 2007 Overview Electrical Impedance
More informationUse of the magnetic field generated by the internal distribution of injected currents for Electrical Impedance Tomography (MR-EIT)
ELEKTRİK, VOL.6, NO.3, 1998, c TÜBİTAK Use of the magnetic field generated by the internal distribution of injected currents for Electrical Impedance Tomography (MR-EIT) Y. Ziya İder & Özlem Birgül Middle
More informationElectric potentials with localized divergence properties
Electric potentials with localized divergence properties Bastian Gebauer bastian.gebauer@oeaw.ac.at Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy of Sciences,
More informationDiscretization-invariant Bayesian inversion and Besov space priors. Samuli Siltanen RICAM Tampere University of Technology
Discretization-invariant Bayesian inversion and Besov space priors Samuli Siltanen RICAM 28.10.2008 Tampere University of Technology http://math.tkk.fi/inverse-coe/ This is a joint work with Matti Lassas
More informationHybrid Inverse Problems and Internal Functionals. Guillaume Bal
Hybrid Inverse Problems and Internal Functionals Department of Applied Physics & Applied Mathematics Columbia University Joint w.: Cédric Bellis, Eric Bonnetier, Matias Courdurier, Chenxi Guo, Sébastien
More informationInverse parameter identification problems
Inverse parameter identification problems Bastian von Harrach harrach@math.uni-stuttgart.de Chair of Optimization and Inverse Problems, University of Stuttgart, Germany ICP - Institute for Computational
More informationQuadrature Formula for Computed Tomography
Quadrature Formula for Computed Tomography orislav ojanov, Guergana Petrova August 13, 009 Abstract We give a bivariate analog of the Micchelli-Rivlin quadrature for computing the integral of a function
More informationA Direct Sampling Method for Electrical Impedance Tomography
A Direct Sampling Method for Electrical Impedance Tomography Yat Tin Chow Kazufumi Ito Jun Zou Abstract This work investigates the electrical impedance tomography (EIT) problem in the case when only one
More informationCONSTRUCTION OF THE HALF-LINE POTENTIAL FROM THE JOST FUNCTION
CONSTRUCTION OF THE HALF-LINE POTENTIAL FROM THE JOST FUNCTION Tuncay Aktosun Department of Mathematics and Statistics Mississippi State University Mississippi State, MS 39762 Abstract: For the one-dimensional
More informationReconstructing inclusions from Electrostatic Data
Reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with: W. Rundell Purdue
More informationA mathematical model and inversion procedure for Magneto-Acousto-Electric Tomography (MAET)
A mathematical model and inversion procedure for Magneto-Acousto-Electric Tomography (MAET) Leonid Kunyansky University of Arizona, Tucson, AZ Suppported in part by NSF grant DMS-0908243 NSF grant "we
More informationConstructing Approximation Kernels for Non-Harmonic Fourier Data
Constructing Approximation Kernels for Non-Harmonic Fourier Data Aditya Viswanathan aditya.v@caltech.edu California Institute of Technology SIAM Annual Meeting 2013 July 10 2013 0 / 19 Joint work with
More informationSolving approximate cloaking problems using finite element methods
Solving approximate cloaking problems using finite element methods Q. T. Le Gia 1 H. Gimperlein M. Maischak 3 E. P. Stephan 4 June, 017 Abstract Motivating from the approximate cloaking problem, we consider
More informationSurvey of Inverse Problems For Hyperbolic PDEs
Survey of Inverse Problems For Hyperbolic PDEs Rakesh Department of Mathematical Sciences University of Delaware Newark, DE 19716, USA Email: rakesh@math.udel.edu January 14, 2011 1 Problem Formulation
More informationCOMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY
COMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY BIN HAN AND HUI JI Abstract. In this paper, we provide a family of compactly supported orthonormal complex wavelets with dilation
More informationReconstruction of inclusions in an elastic body
J. Math. Pures Appl. 91 009 569 58 www.elsevier.com/locate/matpur Reconstruction of inclusions in an elastic body Gunther Uhlmann a,,1, Jenn-Nan Wang b,, Chin-Tien Wu c,3 a Department of Mathematics, University
More informationarxiv: v3 [math.ap] 4 Jan 2017
Recovery of an embedded obstacle and its surrounding medium by formally-determined scattering data Hongyu Liu 1 and Xiaodong Liu arxiv:1610.05836v3 [math.ap] 4 Jan 017 1 Department of Mathematics, Hong
More informationProblems that led me to Gunther Uhlmann. David Isaacson RPI
Problems that led me to Gunther Uhlmann David Isaacson RPI 1. Inverse problem in electrocardiography. 2. Inverse boundary value problem for conductivity. GU Can we improve the diagnosis and treatment of
More informationRiesz bases of Floquet modes in semi-infinite periodic waveguides and implications
Riesz bases of Floquet modes in semi-infinite periodic waveguides and implications Thorsten Hohage joint work with Sofiane Soussi Institut für Numerische und Angewandte Mathematik Georg-August-Universität
More informationEIT Reconstruction Notes
EIT Reconstruction Notes WRB Lionheart September 23, 2010 This introduction is taken mainly from the book chapter Andy Adler, Romina Gaburro, William Lionheart, EIT, chapter in Handbook of Mathematical
More informationInverse problems and medical imaging
Inverse problems and medical imaging Bastian von Harrach harrach@math.uni-frankfurt.de Institute of Mathematics, Goethe University Frankfurt, Germany Seminario di Calcolo delle Variazioni ed Equazioni
More informationCONVEXIFICATION OF COEFFICIENT INVERSE PROBLEMS
CONVEXIFICATION OF COEFFICIENT INVERSE PROBLEMS Michael V. Klibanov Department of Mathematics and Statistics University of North Carolina at Charlotte Charlotte, NC 28223, USA THE MOST RECENT RESULT PDE-based
More informationA new inversion method for dissipating electromagnetic problem
A new inversion method for dissipating electromagnetic problem Elena Cherkaev Abstract The paper suggests a new technique for solution of the inverse problem for Maxwell s equations in a dissipating medium.
More informationApproximation of inverse boundary value problems by phase-field methods
Approximation of inverse boundary value problems by phase-field methods Luca RONDI Università degli Studi di Trieste Dipartimento di Matematica e Informatica MSRI Berkeley, 3 December 2010 Luca RONDI (Università
More informationOn Approximate Cloaking by Nonsingular Transformation Media
On Approximate Cloaking by Nonsingular Transformation Media TING ZHOU MIT Geometric Analysis on Euclidean and Homogeneous Spaces January 9, 2012 Invisibility and Cloaking From Pendry et al s paper J. B.
More informationM.P. Ramirez T., IAENG Member, M. C. Robles G., R. A. Hernandez-Becerril
, July 6-8, 2011, London, U.K. On the Construction of Separable-Variables Conductivity Functions, and Their Application for Approaching Numerical Solutions of the Two-Dimensional Electrical Impedance Equation
More informationPROBING FOR INCLUSIONS IN HEAT CONDUCTIVE BODIES. Patricia Gaitan. Hiroshi Isozaki. Olivier Poisson. Samuli Siltanen.
Inverse Problems and Imaging Volume 6, No. 3, 01, 43 446 doi:10.3934/ipi.01.6.43 PROBING FOR INCLUSIONS IN HEAT CONDUCTIVE BODIES Patricia Gaitan Aix-Marseille University, France Hiroshi Isozaki University
More informationRecent Progress in Electrical Impedance Tomography
Recent Progress in Electrical Impedance Tomography Martin Hanke and Martin Brühl Fachbereich Mathematik und Informatik, Johannes Gutenberg-Universität, 55099 Mainz, Germany E-mail: hanke@math.uni-mainz.de,
More informationIBVPs for linear and integrable nonlinear evolution PDEs
IBVPs for linear and integrable nonlinear evolution PDEs Dionyssis Mantzavinos Department of Applied Mathematics and Theoretical Physics, University of Cambridge. Edinburgh, May 31, 212. Dionyssis Mantzavinos
More informationRESISTOR NETWORKS AND FINITE ELEMENT MODELS
RESISTOR NETWORKS AND FINITE ELEMENT MODELS A thesis submitted to the University of Manchester for the degree of Doctor of Philosophy in the Faculty of Engineering and Physical Sciences 2011 Abdulaziz
More informationNon-uniqueness result for a hybrid inverse problem
Non-uniqueness result for a hybrid inverse problem Guillaume Bal and Kui Ren Abstract. Hybrid inverse problems aim to combine two imaging modalities, one that displays large contrast and one that gives
More informationNumerical differentiation by means of Legendre polynomials in the presence of square summable noise
www.oeaw.ac.at Numerical differentiation by means of Legendre polynomials in the presence of square summable noise S. Lu, V. Naumova, S. Pereverzyev RICAM-Report 2012-15 www.ricam.oeaw.ac.at Numerical
More informationNovel tomography techniques and parameter identification problems
Novel tomography techniques and parameter identification problems Bastian von Harrach harrach@ma.tum.de Department of Mathematics - M1, Technische Universität München Colloquium of the Institute of Biomathematics
More informationLiouville-type theorem for the Lamé system with singular coefficients
Liouville-type theorem for the Lamé system with singular coefficients Blair Davey Ching-Lung Lin Jenn-Nan Wang Abstract In this paper, we study a Liouville-type theorem for the Lamé system with rough coefficients
More informationIterative regularization of nonlinear ill-posed problems in Banach space
Iterative regularization of nonlinear ill-posed problems in Banach space Barbara Kaltenbacher, University of Klagenfurt joint work with Bernd Hofmann, Technical University of Chemnitz, Frank Schöpfer and
More informationAn inverse elliptic problem of medical optics with experimental data
An inverse elliptic problem of medical optics with experimental data Jianzhong Su, Michael V. Klibanov, Yueming Liu, Zhijin Lin Natee Pantong and Hanli Liu Abstract A numerical method for an inverse problem
More informationNon-uniqueness result for a hybrid inverse problem
Non-uniqueness result for a hybrid inverse problem Guillaume Bal and Kui Ren Abstract. Hybrid inverse problems aim to combine two imaging modalities, one that displays large contrast and one that gives
More informationAn eigenvalue method using multiple frequency data for inverse scattering problems
An eigenvalue method using multiple frequency data for inverse scattering problems Jiguang Sun Abstract Dirichlet and transmission eigenvalues have important applications in qualitative methods in inverse
More informationElectrical Impedance Tomography
Electrical Impedance Tomography David Isaacson Jonathan Newell Gary Saulnier RPI Part 1. Adaptive Current Tomography and Eigenvalues With help from D.G.Gisser, M.Cheney, E. Somersalo, J.Mueller, S.Siltanen
More informationThe Imaging of Anisotropic Media in Inverse Electromagnetic Scattering
The Imaging of Anisotropic Media in Inverse Electromagnetic Scattering Fioralba Cakoni Department of Mathematical Sciences University of Delaware Newark, DE 19716, USA email: cakoni@math.udel.edu Research
More informationFast and accurate methods for the discretization of singular integral operators given on surfaces
Fast and accurate methods for the discretization of singular integral operators given on surfaces James Bremer University of California, Davis March 15, 2018 This is joint work with Zydrunas Gimbutas (NIST
More informationCopyright c1999 by Academia Scientiarum Fennica ISSN??????? ISBN??????? Received 4 October Mathematics Subject Classication: 35J5,35R3,78A7,
ANNALES ACADEMISCIENTIARUM FENNIC MATHEMATICA DISSERTATIONES??? ELECTRICAL IMPEDANCE TOMOGRAPHY AND FADDEEV GREEN'S FUNCTIONS SAMULI SILTANEN Dissertation for the Degree of Doctor of Technology to be presented
More information