Computational inversion based on complex geometrical optics solutions
|
|
- Kristopher Wiggins
- 5 years ago
- Views:
Transcription
1 Computational inversion based on complex geometrical optics solutions Samuli Siltanen Department of Mathematics and Statistics University of Helsinki, Finland International Conference on Inverse Problems and Applications in Honor of Gunther Uhlmann s 60th Birthday Hangzhou, China September 19, 2012
2
3 This is a joint work with Takanori Ide, Tokyo Metropolitan University, Japan Masaru Ikehata, Gunma University, Japan David Isaacson, Rensselaer Polytechnic Institute, USA Hiroshi Isozaki, Tsukuba University, Japan Kim Knudsen, Technical University of Denmark Matti Lassas, University of Helsinki, Finland Susumu Nakata, Ritsumeikan University, Japan Jon Newell, Rensselaer Polytechnic Institute, USA Jennifer Mueller, Colorado State University, USA Gunther Uhlmann, University of Washington, USA
4 Outline Electrical impedance tomography The D-bar method for EIT Electrical inclusion detection
5 Outline Electrical impedance tomography The D-bar method for EIT Electrical inclusion detection
6 Electrical impedance tomography (EIT) is an emerging medical imaging technique Feed electric currents through electrodes. Measure the resulting voltages. Repeat the measurement for several current patterns. Reconstruct distribution of electric conductivity inside the patient. Different tissues have different conductivities, so EIT gives an image of the patient s inner structure. EIT is a harmless and painless imaging method suitable for long-term monitoring.
7 This talk concentrates on applications of EIT to chest imaging Applications: monitoring cardiac activity, lung function, and pulmonary perfusion. Also, electrocardiography (ECG) can be enhanced using knowledge about conductivity distribution.
8 The mathematical model of EIT is the inverse conductivity problem introduced by Calderón Let Ω R 2 be the unit disc and let conductivity σ : Ω R satisfy 0 < M 1 σ(z) M. Ω Applying voltage f at the boundary Ω leads to the elliptic PDE { σ u = 0 in Ω, u Ω = f. Boundary measurements are modelled by the Dirichlet-to-Neumann map Λ σ : f σ u n Ω. Calderón s problem is to recover σ from the knowledge of Λ σ. It is a nonlinear and ill-posed inverse problem.
9 Why is Calderón s problem nonlinear? Define a quadratic form P σ for functions f : Ω R by P σ (f ) = σ u 2 dz, (1) where u is the solution of the Dirichlet problem { σ u = 0 in Ω, u Ω = f. Ω Now the map σ P σ is nonlinear because u depends on σ in (1). Physically, P σ (f ) is the power needed for maintaining the voltage potential f on the boundary Ω. Integrate by parts in (1): P σ (f ) = Ω f (σ u n ) ds = Ω Thus the map σ Λ σ cannot be linear in σ. f (Λ σ f ) ds.
10 We illustrate the ill-posedness of Calderón s problem using a simulated example σ 1 σ 2
11 We apply the voltage distribution f (θ) = cos θ at the boundary of the two different phantoms σ 1 u 1 θ σ 2 u 2
12 The measurement is the distribution of current through the boundary σ 1 u 1 θ σ 1 u 1 n θ σ 2 u 2 σ 2 u 2 n θ
13 The measurements are very similar, although the conductivities are quite different σ 1 σ 1 u 1 n σ 2 u 2 n σ 2
14 Let us apply the more oscillatory distribution f (θ) = cos 2θ of voltage at the boundary σ 1 u 1 σ 2 u 2
15 The measurement is again the distribution of current through the boundary σ 1 u 1 σ 1 u 1 n θ σ 2 u 2 σ 2 u 2 n θ
16 The current distribution measurements are again very similar σ 1 σ 1 u 1 n σ 2 u 2 n σ 2
17 EIT is an ill-posed problem: big differences in conductivity cause only small effect in data σ 1 cos θ cos 4θ cos 2θ cos 5θ σ 2 cos 3θ cos 6θ
18 EIT is an ill-posed problem: noise in data causes serious difficulties in interpreting the data σ 1 cos θ cos 4θ cos 2θ cos 5θ σ 2 cos 3θ cos 6θ
19 Many different types of reconstruction methods have been suggested for EIT in the literature Linearization: Barber, Bikowski, Brown, Calderón, Cheney, Isaacson, Mueller, Newell Iterative regularization: Dobson, Gehre, Hua, Jin, Kaipio, Kindermann, Kluth, Leitão, Lechleiter, Lipponen, Maass, Neubauer, Rieder, Rondi, Santosa, Seppänen, Tompkins, Webster, Woo Bayesian inversion: Fox, Kaipio, Kolehmainen, Nicholls, Pikkarainen, Ronkanen, Somersalo, Vauhkonen, Voutilainen Resistor network methods: Borcea, Druskin, Mamonov, Vasquez Layer stripping: Cheney, Isaacson, Isaacson, Somersalo D-bar methods: Astala, Bikowski, Bowerman, Delbary, Hansen, Isaacson, Kao, Knudsen, Lassas, Mueller, Murphy, Nachman, Newell, Päivärinta, Perämäki, Saulnier, S, Tamasan Teichmüller space methods: Kolehmainen, Lassas, Ola, S Methods for partial information: Alessandrini, Ammari, Bilotta, Brühl, Eckel, Erhard, Gebauer, Hanke, Harrach, Hyvönen, Ide, Ikehata, Isozaki, Kang, Kim, Kress, Kwon, Lechleiter, Lim, Morassi, Nakamura, Nakata, Potthast, Rossetand, Seo, Sheen, S, Turco, Uhlmann, Wang, and others
20 The forward map F : X D(F ) Y of an ill-posed problem does not have a continuous inverse Model space X = L (Ω) Data space Y σ F Λ σ δ Λ σ δ D(F ) F (D(F ))
21 Regularization means constructing a continuous map Γ α : Y X that inverts F approximately Model space X = L (Ω) Data space Y F σ Λ σ Λδ δ σ Γ α D(F ) F (D(F ))
22 A regularization strategy needs to be constructed so that the assumptions below are satisfied A set Γ α : Y X of continuous maps parameterized by α > 0 is a regularization strategy for F if for each fixed σ D(F ) we have lim Γ α(λ σ ) σ X = 0. α 0 Further, a regularization strategy with a choice α = α(δ) of regularization parameter is called admissible if α(δ) 0 as δ 0, and for any fixed σ D(F ) the following holds: { sup Γα(δ) (Λ δ σ) σ X : Λ δ σ Λ σ Y δ } 0 as δ 0. Λ δ σ
23 Outline Electrical impedance tomography The D-bar method for EIT Electrical inclusion detection
24 History of CGO-based methods for real 2D EIT Infinite-precision data Practical data 1980 Calderón 2008 Bikowski-Mueller 1987 Sylvester-Uhlmann (d 3) 2008 Boverman-Isaacson-Kao-Saulnier-Newell 1988 Nachman 2010 Bikowski-Knudsen-Mueller 1988 R. G. Novikov 2011 Delbary-Hansen-Knudsen 1996 Nachman (σ C 2 (Ω)) 2000 S-Mueller-Isaacson 1997 Liu 2003 Mueller-S 2004 Isaacson-Mueller-Newell-S 2006 Isaacson-Mueller-Newell-S 2007 Murphy-Mueller 2008 Knudsen-Lassas-Mueller-S 2009 Knudsen-Lassas-Mueller-S 2012 S-Tamminen 1997 Brown-Uhlmann (σ C 1 (Ω)) 2001 Knudsen-Tamasan 2001 Barceló-Barceló-Ruiz 2003 Knudsen 2003 Astala-Päivärinta (σ L (Ω)) 2009 Astala-Mueller-Päivärinta-S 2005 Astala-Lassas-Päivärinta 2011 Astala-Mueller-Päivärinta Barceló-Faraco-Ruiz Perämäki-S 2008 Clop-Faraco-Ruiz
25 Nachman s 1996 uniqueness proof in 2D uses complex geometric optics (CGO) solutions Define a potential q by setting q(z) 0 for z outside Ω and q(z) = σ(z) σ(z) for z Ω. Then q C 0 (Ω). We look for solutions of the Schrödinger equation ( + q)ψ(, k) = 0 in R 2 parametrized by k C \ 0 and satisfying the asymptotic condition e ikz ψ(z, k) 1 W 1, p (R 2 ), where p > 2 and ikz = i(k 1 + ik 2 )(x + iy).
26 The CGO solutions are constructed using a generalized Lippmann-Schwinger equation Define µ(z, k) = e ikz ψ(z, k). Then ( + q)ψ = 0 implies ( 4ik z + q)µ(, k) = 0, (2) where the D-bar operator is defined by z = 1 2 ( x + i y ). A solution of (2) satisfying µ(z, k) 1 W 1, p (R 2 ) can be constructed using the Lippmann-Schwinger type equation µ = 1 g k (qµ), where g k satisfies ( 4ik z )g k = δ and is defined by g k (z) = 1 4π 2 R2 e iz ξ ξ 2 + 2k(ξ 1 + iξ 2 ) dξ 1dξ 2.
27 The Faddeev fundamental solution g 1 (z) has a logarithmic singularity at z = 0 Real part of g 1 (z) Imaginary part of g 1 (z) It is enough to know g 1 (z) because of the relation g k (z) = g 1 (kz).
28 The conductivity σ can be recovered from the functions µ(z, k) at k = 0 Recall that with the asymptotics Substituting k = 0 gives ( 4ik z + q)µ(, k) = 0 µ(z, k) 1 W 1, p (R 2 ). ( + σ σ )µ(, 0) = 0, (3) and setting µ(z, 0) = σ(z) gives the unique solution of (3) satisfying µ(z, 0) 1 W 1, p (R 2 ).
29 The crucial intermediate object in the proof is the non-physical scattering transform t(k) We denote z = x + iy C or z = (x, y) R 2 whenever needed. The scattering transform t : C C is defined by t(k) := e ikz q(z)ψ(z, k) dxdy. (4) R 2 Sometimes (4) is called the nonlinear Fourier transform of q. This is because asymptotically ψ(z, k) e ikz as z, and substituting e ikz in place of ψ(z, k) into (4) results in e i(kz+kz) q(z)dxdy R 2 = e i( 2k 1,2k 2 ) (x,y) q(z)dxdy R 2 = q( 2k 1, 2k 2 ).
30 Infinite-precision data: Solve boundary integral equation ψ(, k) Ω = e ikz S k (Λ σ Λ 1 )ψ for every complex number k C \ 0. Practical data: Solve boundary integral equation ψ δ (, k) Ω = e ikz S k (Λ δ σ Λ 1 )ψ δ for all 0 < k < R = 1 10 log δ. Evaluate the scattering transform: t(k) = e ikz (Λ σ Λ 1 )ψ(, k) ds. Ω Fix z Ω. Solve D-bar equation k µ(z, k) = t(k) 4πk e i(kz+kz) µ(z, k) with µ(z, ) 1 L r L (C). For k R set t δ R (k) = 0. For k < R t δ R(k) = e ikz (Λ δ σ Λ 1 )ψ δ (, k) ds. Ω Fix z Ω. Solve D-bar equation k µδ R(z, k) = tδ R (k) 4πk e i(kz+kz) µ δ R (z, k) with µ δ R (z, ) 1 Lr L (C). Reconstruct: σ(z) = (µ(z, 0)) 2. Set Γ 1/R(δ) (Λ δ σ) := (µ δ R (z, 0))2.
31 We define spaces for our regularization strategy = L (Ω) Let M > 0 and 0 < ρ < 1. The domain D(F ) consists of functions σ : Ω R with σ C 2 (Ω) M, σ(z) M 1, σ(z) 1 for ρ < z < 1. Bounded linear operators A : H 1/2 ( Ω) H 1/2 ( Ω) satisfying A(1) = 0, A(f ) ds = 0. Ω
32 Main result: nonlinear low-pass filtering yields a regularization strategy with convergence speed Theorem (Knudsen, Lassas, Mueller & S 2009) There exists a constant 0 < δ 0 < 1, depending only on M and ρ, with the following properties. Let σ D(F ) be arbitrary and assume given noisy data Λ δ σ satisfying Then Γ α with the choice R(δ) = 1 10 Λ δ σ Λ σ Y δ < δ 0. log δ, α(δ) = 1 R(δ), is well-defined, admissible and satisfies the estimate Γ α(δ) (Λ δ σ) σ L (Ω) C( log δ) 1/14.
33 Let us analyze how the regularization works using a simulated heart-and-lungs phantom
34 This is how the actual scattering transform looks like in the disc k < 10, computed by knowing σ Real part of t(k) Imaginary part
35 Scattering transform in the disc k < 10, here computed from noisy measurement Λ δ σ Real part of t(k) Imaginary part
36 Regularized reconstructions from simulated data with noise amplitude δ = Λ δ σ Λ σ Y δ 10 6 δ 10 5 δ 10 4 δ 10 3 δ 10 2 The percentages are the relative square norm errors in the reconstructions.
37 The method works for real data as well, including laboratory phantoms and in vivo human data Saline and agar phantom Reconstruction (R = 4) [Isaacson, Mueller, Newell & S 2006] [Montoya & Mueller 2012]
38 Unknown boundary shape can be estimated from EIT data using Teichmüller space methods [Kolehmainen, Lassas, Ola & S 2012]
39 Outline Electrical impedance tomography The D-bar method for EIT Electrical inclusion detection
40 This is a brief history of the enclosure method 2000 Ikehata Enclosure method for EIT } 2000 Ikehata & S Implementation for 2D EIT 2000 Brühl & Hanke 2002 Ikehata Regularization analysis 2002 Ikehata & Ohe Polygonal cavities in 2D EIT 2004 Ikehata & S Regularized cone probing in 2D EIT 2007 Ide, Isozaki, Nakata, S & Uhlmann Spherical probing in 2D EIT 2008 Ikehata & Ohe Locating cracks in 2D EIT 2008 Uhlmann & Wang General 2D probing 2010 Ide, Isozaki, Nakata & S Spherical probing in 3D EIT 2010 Zhou Maxwell s equations 2012 Ikehata, Niemi & S Inverse obstacle scattering in 2D
41 Different probing strategies Half-space Cone Spherical
42 Main theorem of spherical probing, part 1 (Ide, Isozaki, Nakata, S and Uhlmann 2007) We denote the inclusion set by Ω 1. The background conductivity is σ 0 1. Assume that σ(x) = σ 0 (x) for x Ω 1 and σ(x) > σ 0 (x) on Ω 1. Ω σ σ 0 Take x 0 from the outside of the convex hull of Ω, and let R > 0. Then there exists a smooth solution u τ (x) of (σ 0 (x) u τ (x)) = 0 on Ω with the following special properties. Ω 1 σ > σ 0 x 0 R
43 Main theorem of spherical probing, part 2 Hyperbolic transformation: y 1 = x x 2 2 r 2 (x 1 + r) 2 + x2 2, y 2 = 2x 2 r (x 1 + r) 2 + x 2 2 (x 1, x 2 )-plane: u τ (x) (y 1, y 2 )-plane: exp( τy 1 + iτy 2 ) Decay Growth Decay Growth
44 Main theorem of spherical probing, part 3 Let f τ = u τ Ω. There is a 0 < δ = δ(ω 1, x 0, R) such that (Λ σ Λ σ0 )f τ, f τ < Ce δτ (Λ σ Λ σ0 )f τ, f τ > Ce δτ
45 Practical scheme for detecting inclusions Take 0 < τ 1 < τ 2. (Λ σ Λ σ0 )f τ, f τ < Ce δτ (Λ σ Λ σ0 )f τ, f τ > Ce δτ (Λ σ Λ σ0 )f τ1, f τ1 > (Λ σ Λ σ0 )f τ2, f τ2 (Λ σ Λ σ0 )f τ1, f τ1 < (Λ σ Λ σ0 )f τ2, f τ2
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194 Spherical probing from ideal and noisy EIT data Conductivity Maximal region Reconstruction (ideal data) Reconstruction (noisy data)
195 Spherical probing from ideal and noisy EIT data in 3D (Ide, Isozaki, Nakata & S 2010)
196
197 T. Ide, H. Isozaki, S. Nakata, and S. Siltanen, Local detection of three-dimensional inclusions in electrical impedance tomography, Inverse problems, 26 (2010), p T. Ide, H. Isozaki, S. Nakata, S. Siltanen, and G. Uhlmann, Probing for electrical inclusions with complex spherical waves, Communications on pure and applied mathematics, 60 (2007), pp M. Ikehata, E. Niemi, and S. Siltanen, Inverse obstacle scattering with limited-aperture data, Inverse Problems and Imaging, 6 (2012), pp M. Ikehata and S. Siltanen, Numerical method for finding the convex hull of an inclusion in conductivity from boundary measurements, Inverse Problems, 16 (2000), pp M. Ikehata and S. Siltanen, Electrical impedance tomography and Mittag-Leffler s function, Inverse Problems, 20 (2004), pp D. Isaacson, J. Mueller, J. Newell, and S. Siltanen, Imaging cardiac activity by the D-bar method for electrical impedance tomography, Physiological Measurement, 27 (2006), pp. S43 S50. 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 Transactions on Medical Imaging, 23 (2004), pp K. Knudsen, M. Lassas, J. Mueller, and S. Siltanen, Regularized D-bar method for the inverse conductivity problem, Inverse Problems and Imaging, 3 (2009), pp J. Mueller and S. Siltanen, Direct reconstructions of conductivities from boundary measurements, SIAM Journal on Scientific Computing, 24 (2003), pp A. I. Nachman, Global uniqueness for a two-dimensional inverse boundary value problem, Annals of Mathematics, 143 (1996), pp 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), pp
Computational 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 informationReconstructing conductivities with boundary corrected D-bar method
Reconstructing conductivities with boundary corrected D-bar method Janne Tamminen June 24, 2011 Short introduction to EIT The Boundary correction procedure The D-bar method Simulation of measurement data,
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 informationReconstructing 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationRecent progress on the factorization method for electrical impedance tomography
Recent progress on the factorization method for electrical impedance tomography Bastian von Harrach harrach@math.uni-stuttgart.de Chair of Optimization and Inverse Problems, University of Stuttgart, Germany
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationSampling methods for low-frequency electromagnetic imaging
Sampling methods for low-frequency electromagnetic imaging Bastian Gebauer gebauer@math.uni-mainz.de Institut für Mathematik, Joh. Gutenberg-Universität Mainz, Germany Joint work with Jin Keun Seo, Yonsei
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 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 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 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 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 informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Mathematics 236 (2012) 4645 4659 Contents lists available at SciVerse ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationarxiv: v1 [math.ap] 13 Jun 2017
ETECTING STOCHSTIC INCLUSIONS IN ELECTRICL IMPENCE TOMOGRPHY NRE BRTH, BSTIN HRRCH, NUUTTI HYVÖNEN, N LURI MUSTONEN arxiv:706.0396v [math.p] 3 Jun 07 bstract. This work considers the inclusion detection
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 informationOn hp-adaptive Solution of Complete Electrode Model Forward Problems of Electrical Impedance Tomography
On hp-adaptive Solution of Complete Electrode Model Forward Problems of Electrical Impedance Tomography Harri Hakula, Nuutti Hyvönen, Tomi Tuominen Aalto University School of Science Department of Mathematics
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 informationResolving the White Noise Paradox in the Regularisation of Inverse Problems
1 / 32 Resolving the White Noise Paradox in the Regularisation of Inverse Problems Hanne Kekkonen joint work with Matti Lassas and Samuli Siltanen Department of Mathematics and Statistics University of
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 informationSimultaneous reconstruction of outer boundary shape and conductivity distribution in electrical impedance tomography
Simultaneous reconstruction of outer boundary shape and conductivity distribution in electrical impedance tomography Nuutti Hyvönen Aalto University nuutti.hyvonen@aalto.fi joint work with J. Dardé, A.
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 informationTrends in hybrid data tomography Workshop at DTU Compute Wednesday January 24, 2018 Room 324/050
Trends in hybrid data tomography Workshop at DTU Compute Wednesday January 24, 2018 Room 324/050 Program 09:00 09:50 Stochastic Gradient Descent for Inverse Problems Bangti Jin, University College London
More informationA local estimate from Radon transform and stability of Inverse EIT with partial data
A local estimate from Radon transform and stability of Inverse EIT with partial data Alberto Ruiz Universidad Autónoma de Madrid ge Joint work with P. Caro (U. Helsinki) and D. Dos Santos Ferreira (Paris
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 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 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 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 informationModel-aware Newton-type regularization in electrical impedance tomography
Model-aware Newton-type regularization in electrical impedance tomography Andreas Rieder Robert Winkler FAKULTÄT FÜR MATHEMATIK INSTITUT FÜR ANGEWANDTE UND NUMERISCHE MATHEMATIK SFB 1173 KIT University
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 informationc 2009 Society for Industrial and Applied Mathematics
SIAM J. APPL. MATH. Vol. 69, No. 6, pp. 1662 1681 c 2009 Society for Industrial and Applied Mathematics DETECTING INCLUSIONS IN ELECTICAL IMPEDANCE TOMOGAPHY WITHOUT EFEENCE MEASUEMENTS ASTIAN HAACH AND
More informationKey words. Electrical impedance tomography, stochastic conductivity, inclusion detection, factorization method, monotonicity method
ETECTING STOCHSTIC INCLUSIONS IN ELECTRICL IMPENCE TOMOGRPHY NRE BRTH, BSTIN HRRCH, NUUTTI HYVÖNEN, N LURI MUSTONEN bstract. This work considers the inclusion detection problem of electrical impedance
More informationRegularization in Banach Space
Regularization in Banach Space Barbara Kaltenbacher, Alpen-Adria-Universität Klagenfurt joint work with Uno Hämarik, University of Tartu Bernd Hofmann, Technical University of Chemnitz Urve Kangro, University
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 informationInverse Transport Problems and Applications. II. Optical Tomography and Clear Layers. Guillaume Bal
Inverse Transport Problems and Applications II. Optical Tomography and Clear Layers Guillaume Bal Department of Applied Physics & Applied Mathematics Columbia University http://www.columbia.edu/ gb23 gb23@columbia.edu
More informationBayesian Methods and Uncertainty Quantification for Nonlinear Inverse Problems
Bayesian Methods and Uncertainty Quantification for Nonlinear Inverse Problems John Bardsley, University of Montana Collaborators: H. Haario, J. Kaipio, M. Laine, Y. Marzouk, A. Seppänen, A. Solonen, Z.
More informationA local estimate from Radon transform and stability of Inverse EIT with partial data
A local estimate from Radon transform and stability of Inverse EIT with partial data Alberto Ruiz Universidad Autónoma de Madrid U. California, Irvine.June 2012 w/ P. Caro (U. Helsinki) and D. Dos Santos
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 informationAn inverse scattering problem in random media
An inverse scattering problem in random media Pedro Caro Joint work with: Tapio Helin & Matti Lassas Computational and Analytic Problems in Spectral Theory June 8, 2016 Outline Introduction and motivation
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 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 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 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 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 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 informationFactorization method in inverse
Title: Name: Affil./Addr.: Factorization method in inverse scattering Armin Lechleiter University of Bremen Zentrum für Technomathematik Bibliothekstr. 1 28359 Bremen Germany Phone: +49 (421) 218-63891
More informationEnclosure methods for Helmholtz-type equations
Inside Out II MSRI Publications Volume 60, 2012 Enclosure methods for Helmholtz-type equations JENN-NAN WANG AN TING HOU The inverse problem under consideration is to reconstruct the shape information
More informationThe Factorization Method for Inverse Scattering Problems Part I
The Factorization Method for Inverse Scattering Problems Part I Andreas Kirsch Madrid 2011 Department of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research Center
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 informationJUSTIFICATION OF POINT ELECTRODE MODELS IN ELECTRICAL IMPEDANCE TOMOGRAPHY
JUSTIFICATION OF POINT ELECTRODE MODELS IN ELECTRICAL IMPEDANCE TOMOGRAPHY MARTIN HANKE, BASTIAN HARRACH, AND NUUTTI HYVÖNEN Abstract. The most accurate model for real-life electrical impedance tomography
More informationToward imaging modalities with high(er) resolution
Toward imaging modalities with high(er) resolution François Monard Department of Mathematics, University of Michigan February 9, 2016 University of North Carolina at Charlotte Introduction Topic du jour
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 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 informationThe inverse conductivity problem with power densities in dimension n 2
The inverse conductivity problem with power densities in dimension n 2 François Monard Guillaume Bal Dept. of Applied Physics and Applied Mathematics, Columbia University. June 19th, 2012 UC Irvine Conference
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 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 informationThe Linear Sampling Method and the MUSIC Algorithm
CODEN:LUTEDX/(TEAT-7089)/1-6/(2000) The Linear Sampling Method and the MUSIC Algorithm Margaret Cheney Department of Electroscience Electromagnetic Theory Lund Institute of Technology Sweden Margaret Cheney
More informationMonotonicity-based inverse scattering
Monotonicity-based inverse scattering Bastian von Harrach http://numerical.solutions Institute of Mathematics, Goethe University Frankfurt, Germany (joint work with M. Salo and V. Pohjola, University of
More informationAn inverse source problem in optical molecular imaging
An inverse source problem in optical molecular imaging Plamen Stefanov 1 Gunther Uhlmann 2 1 2 University of Washington Formulation Direct Problem Singular Operators Inverse Problem Proof Conclusion Figure:
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 informationSparse tomography. Samuli Siltanen. Department of Mathematics and Statistics University of Helsinki, Finland
Sparse tomography Samuli Siltanen Department of Mathematics and Statistics University of Helsinki, Finland Minisymposium: Fourier analytic methods in tomographic image reconstruction Applied Inverse Problems
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 informationImproved Rotational Invariance for Statistical Inverse in Electrical Impedance Tomography
Improved Rotational Invariance for Statistical Inverse in Electrical Impedance Tomography Jani Lahtinen, Tomas Martinsen and Jouko Lampinen Laboratory of Computational Engineering Helsinki University of
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 informationReconstruction Scheme for Active Thermography
Reconstruction Scheme for Active Thermography Gen Nakamura gnaka@math.sci.hokudai.ac.jp Department of Mathematics, Hokkaido University, Japan Newton Institute, Cambridge, Sept. 20, 2011 Contents.1.. Important
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 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 Problems and Applications
Inverse Problems and Applications June 3 9, 2002 Cortona, Italy Organizers Giovanni Alessandrini (Università degli Studi di Trieste, Italy) Gunther Uhlmann (University of Washington, USA) Inverse Problems
More information. LIMITED VIEW RESOLVING POWER OF CONDUCTIVITY IMAGING FROM BOUNDARY MEASUREMENTS
. LIMITED VIEW RESOLVING POWER OF CONDUCTIVITY IMAGING FROM BOUNDARY MEASUREMENTS HABIB AMMARI JOSSELIN GARNIER AND KNUT SØLNA Abstract. In this paper we consider resolution estimates for a linearized
More informationA CONVERGENCE ANALYSIS OF THE NEWTON-TYPE REGULARIZATION CG-REGINN WITH APPLICATION TO IMPEDANCE TOMOGRAPHY
A CONVERGENCE ANALYSIS OF THE NEWTON-TYPE REGULARIZATION CG-REGINN WITH APPLICATION TO IMPEDANCE TOMOGRAPHY ARMIN LECHLEITER AND ANDREAS RIEDER January 8, 2007 Abstract. The Newton-type regularization
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationEnhanced Newton-Raphson Method for Inverse Problems Solution of Electrical Impedance Tomography
2015, TextRoad Publication ISSN: 2090-4274 Journal of Applied Environmental and Biological Sciences www.textroad.com Enhanced Newton-Raphson Method for Inverse Problems Solution of Electrical Impedance
More information