![Page 1: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/1.jpg)
Electrical Impedance Tomography
and Hybrid Imaging
Allan Greenleaf
University of Rochester, USA
LMS–EPSRC Durham Symposium
Mathematical and Computational Aspects of Maxwell’s Equations
July 15, 2016
Partially supported by DMS-1362271 and a Simons Foundation Fellowship
![Page 2: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/2.jpg)
Electrical impedance tomography (EIT)
Calderon’s inverse conductivity problem:Imaging an electrical conductivity σ(x) vianoninvasive voltage/current measurements atthe surface of an object.
![Page 3: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/3.jpg)
Electrical impedance tomography (EIT)
Calderon’s inverse conductivity problem:Imaging an electrical conductivity σ(x) vianoninvasive voltage/current measurements atthe surface of an object.
+ Major theoretical and numerical advancesover last 35 years.
– Plain EIT has seen limited application inclinical/industrial settings.
![Page 4: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/4.jpg)
Electrical impedance tomography (EIT)
Calderon’s inverse conductivity problem:Imaging an electrical conductivity σ(x) vianoninvasive voltage/current measurements atthe surface of an object.
+ Major theoretical and numerical advancesover last 35 years.
– Plain EIT has seen limited application inclinical/industrial settings.
→ Hybrid imaging developed to overcomedisadvantages of EIT and other modalities.
![Page 5: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/5.jpg)
• Interior of a region Ω ⊂ Rn, n = 2, 3,filled with matter having conductivity σ(x).(Ω = human body, industrial part,...)
![Page 6: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/6.jpg)
• Interior of a region Ω ⊂ Rn, n = 2, 3,filled with matter having conductivity σ(x).(Ω = human body, industrial part,...)
• Place electrodes on the boundary, ∂Ω.
Connect to DC sources to create aprescribed voltage distribution, f , on ∂Ω.
f induces a electric potential u(x) in Ω.
![Page 7: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/7.jpg)
• Interior of a region Ω ⊂ Rn, n = 2, 3,filled with matter having conductivity σ(x).(Ω = human body, industrial part,...)
• Place electrodes on the boundary, ∂Ω.
Connect to DC sources to create aprescribed voltage distribution, f , on ∂Ω.
f induces a electric potential u(x) in Ω.
•Measure resulting current flow I across ∂Ω.Ohm’s Law =⇒
I = σ · ∂u∂ν
![Page 8: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/8.jpg)
Quasi-static regime: Electric potential u(x)satisfies conductivity equation,
∇ · (σ∇u)(x) = 0 on Ω,
with Dirichlet boundary condition
u|∂Ω = f = prescribed voltage on ∂Ω.
![Page 9: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/9.jpg)
Quasi-static regime: Electric potential u(x)satisfies conductivity equation,
∇ · (σ∇u)(x) = 0 on Ω,
with Dirichlet boundary condition
u|∂Ω = f = prescribed voltage on ∂Ω.
Dirichlet-to-Neumann operator
f −→ σ · ∂u∂ν
=: Λσ(f ) on ∂Ω.
Λσ : H12(∂Ω)→ H−
12(∂Ω) bounded lin. oper.
![Page 10: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/10.jpg)
Calderon’s Inverse Problem (Isotropic)
(i) Uniqueness: Does Λσ1 = Λσ2 =⇒ σ1 = σ2?
![Page 11: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/11.jpg)
Calderon’s Inverse Problem (Isotropic)
(i) Uniqueness: Does Λσ1 = Λσ2 =⇒ σ1 = σ2?
(ii) Reconstruction: Can we find σ(x) from Λσ?
![Page 12: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/12.jpg)
Calderon’s Inverse Problem (Isotropic)
(i) Uniqueness: Does Λσ1 = Λσ2 =⇒ σ1 = σ2?
(ii) Reconstruction: Can we find σ(x) from Λσ?
(iii, . . . ) Stability of Λσ → σ, numerics, . . .
![Page 13: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/13.jpg)
Calderon’s Inverse Problem (Isotropic)
(i) Uniqueness: Does Λσ1 = Λσ2 =⇒ σ1 = σ2?
(ii) Reconstruction: Can we find σ(x) from Λσ?
(iii, . . . ) Stability of Λσ → σ, numerics, . . .
A: Yes to (i), (ii), but poor stability.
![Page 14: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/14.jpg)
Progress on isotropic Calderon problem
1980, Calderon: linearization around σ ≡ 1
1984, Kohn and Vogelius: uniqueness forpiecewise -Cω conductivities
![Page 15: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/15.jpg)
Progress on isotropic Calderon problem
1980, Calderon: linearization around σ ≡ 1
1984, Kohn and Vogelius: uniqueness forpiecewise -Cω conductivities
1986, Sylvester and Uhlmann: uniquenessfor σ ∈ C2, n ≥ 3. Introduced CGO solutions.
![Page 16: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/16.jpg)
Progress on isotropic Calderon problem
1980, Calderon: linearization around σ ≡ 1
1984, Kohn and Vogelius: uniqueness forpiecewise -Cω conductivities
1986, Sylvester and Uhlmann: uniquenessfor σ ∈ C2, n ≥ 3. Introduced CGO solutions.
1988, Nachman: reconstruction, n ≥ 3
1996, Nachman: uniqueness+reconstr., n = 2
![Page 17: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/17.jpg)
2006, Astala and Pavarinta: uniqueness andreconstruction for σ ∈ L∞, n = 2
![Page 18: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/18.jpg)
2006, Astala and Pavarinta: uniqueness andreconstruction for σ ∈ L∞, n = 2
For n = 3:
2013, Haberman and Tataru: uniqueness forσ ∈ C1 or Lipschitz close to constant.
2015, Caro and Rogers: ! for Lipschitz σ.
2015, Haberman: ! for σ ∈ W 1,3+ε.
![Page 19: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/19.jpg)
2006, Astala and Pavarinta: uniqueness andreconstruction for σ ∈ L∞, n = 2
For n = 3:
2013, Haberman and Tataru: uniqueness forσ ∈ C1 or Lipschitz close to constant.
2015, Caro and Rogers: ! for Lipschitz σ.
2015, Haberman: ! for σ ∈ W 1,3+ε.
Q.: Does uniqueness hold for σ ∈ L∞ ?
![Page 20: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/20.jpg)
Problem: EIT has high contrast sensitivity,but low spatial resolution.
Figure 1: EIT tank and measurements. Source: Kaipio lab, Univ. of Kuopio, Finland
![Page 21: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/21.jpg)
Hybrid inverse problems
Q: Can one improve imaging by using datafrom more than one type of wave?
(i) Image registration, e.g., CT+MRI
![Page 22: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/22.jpg)
Hybrid inverse problems
Q: Can one improve imaging by using datafrom more than one type of wave?
(i) Image registration, e.g., CT+MRI
(ii) Stabilization: collect data for two typesof waves, X and Y , simultaneously. Eitheruse
• Y data to provide a priori informationthat stabilizes reconstruction from X data;
or
• an algorithm using both X and Y data.
![Page 23: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/23.jpg)
Ex.: Current Density Impedance Imaging(CDI) - Magnetic Resonance EIT (MREIT):
Measure both voltage/current at ∂Ω andcurrent density σ|∇u| in the interior (via MRI).
→ J-substitution algo. of Kwon, Woo, et al.
![Page 24: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/24.jpg)
Ex.: Current Density Impedance Imaging(CDI) - Magnetic Resonance EIT (MREIT):
Measure both voltage/current at ∂Ω andcurrent density σ|∇u| in the interior (via MRI).
→ J-substitution algo. of Kwon, Woo, et al.
However, want to discuss
(iii) ‘Multi-physics’ hybrid methods in whichtwo different kinds of waves are physicallycoupled.
![Page 25: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/25.jpg)
Multi-physics methods often combine twoillumination and detection modalities,one with
• high contrast sensitivity but low resolution,
![Page 26: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/26.jpg)
Multi-physics methods often combine twoillumination and detection modalities,one with
• high contrast sensitivity but low resolution,
and the other one
• low contrast but high resolution,
linked by a physical interaction.
![Page 27: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/27.jpg)
Multi-physics methods often combine twoillumination and detection modalities,one with
• high contrast sensitivity but low resolution,
and the other one
• low contrast but high resolution,
linked by a physical interaction.
Mathematically: couple an elliptic PDEwith a hyperbolic PDE.
![Page 28: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/28.jpg)
Thermo-acoustic tomography (TAT)
Illuminate object with short microwave pulse.EM energy is absorbed preferentially bysubregions of interest, e.g., tumors.
![Page 29: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/29.jpg)
Thermo-acoustic tomography (TAT)
Illuminate object with short microwave pulse.EM energy is absorbed preferentially bysubregions of interest, e.g., tumors.
Photo-acoustic effect: Thermal expansionproduces acoustic waves (often ultrasound)with sources at loci of high EM absorption.
![Page 30: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/30.jpg)
Thermo-acoustic tomography (TAT)
Illuminate object with short microwave pulse.EM energy is absorbed preferentially bysubregions of interest, e.g., tumors.
Photo-acoustic effect: Thermal expansionproduces acoustic waves (often ultrasound)with sources at loci of high EM absorption.
Acoustic waves then propagate out to ∂Ω,where measured.
EM governed by diffusion eqn. (elliptic),US by acoustic wave eqn. (hyperbolic).
![Page 31: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/31.jpg)
• Solve hyperbolic inverse problem for US.
Reconstructs with good spatial resolutionan internal measurement: a functional F (x, u,∇u)of the solution u(x) of the elliptic problem forthe EM field.
• Then solve the elliptic inverse problem offinding absorption coefficient in Ω from
• u on ∂Ω
• F on Ω
• Other a priori information/assumptions
![Page 32: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/32.jpg)
Photo-acoustic Tomography (PAT):
Illumination by infrared EM,detection by ultrasound.
![Page 33: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/33.jpg)
Photo-acoustic Tomography (PAT):
Illumination by infrared EM,detection by ultrasound.
Ultrasound Modulated Optical Tomography(UMOT):
Illumination by ultrasound,detection by infrared.
![Page 34: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/34.jpg)
Photo-acoustic Tomography (PAT):
Illumination by infrared EM,detection by ultrasound.
Ultrasound Modulated Optical Tomography(UMOT):
Illumination by ultrasound,detection by infrared.
Acousto-Electric Tomography (AET/UMEIT):
Illumination by ultrasound,detection by EIT.
![Page 35: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/35.jpg)
Model of PAT
• Illuminate with short pulse.Scalar EM field in Ω satisfies
−∇ · (σ(x)∇u(x)) + a(x)u(x) = 0,
u|∂Ω (known)
a(x) = absorption coeff. (desired)
σ(x) = diffusion coeff.
![Page 36: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/36.jpg)
• Resulting pressure p(x, t) satisfies
(∂2t − c(x)2∆
)p(x, t) = 0 on Ω× [0,∞)
p(x, 0) = F (x, u(x)), ∂tp(x, 0) = 0
F = Γ(x)a(x)u(x), where Γ(x) = Gruneisen coeff.
Then: (1) solve hyperbolic IP and find F (x, u(x))from p|∂Ω×[0,T0]
(2) solve problem finding a(x) from F, u|∂Ω, Γ
![Page 37: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/37.jpg)
Real principal type (RPT) operators
P (x,D) ∈ Ψm(Rn), n ≥ 2, is of RPT if
(i) principal symbol pm(x, ξ) is R-valued
(ii) dpm(x, ξ) 6= (0, 0) at
ΣP = (x, ξ) ∈ T ∗Rn, ξ 6= 0 : pm(x, ξ) = 0
![Page 38: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/38.jpg)
Real principal type (RPT) operators
P (x,D) ∈ Ψm(Rn), n ≥ 2, is of RPT if
(i) principal symbol pm(x, ξ) is R-valued
(ii) dpm(x, ξ) 6= (0, 0) at
ΣP = (x, ξ) ∈ T ∗Rn, ξ 6= 0 : pm(x, ξ) = 0Thus, ΣP is foliated by bicharacteristics :=integral curves of
Hpm :=∑ ∂pm
∂ξj
∂
∂xj− ∂pm∂xj
∂
∂ξj
(iii) No bichar is trapped over a compactset K ⊂ Rn.
![Page 39: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/39.jpg)
Duistermaat and Hormander (FIO II):constructed parametrices for RPT ops,showed they are locally solvable,and singularities of Pu = f propagate alongthe bicharacteristics.
![Page 40: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/40.jpg)
Duistermaat and Hormander (FIO II):constructed parametrices for RPT ops,showed they are locally solvable,and singularities of Pu = f propagate alongthe bicharacteristics.
Thm. For all f ∈ E ′(X), Pu = f is solvable,and if (x, ξ) ∈ WF (u) \ WF (f ), then WF (u)contains the bicharacteristic through (x, ξ).
![Page 41: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/41.jpg)
Duistermaat and Hormander (FIO II):constructed parametrices for RPT ops,showed they are locally solvable,and singularities of Pu = f propagate alongthe bicharacteristics.
Thm. For all f ∈ E ′(X), Pu = f is solvable,and if (x, ξ) ∈ WF (u) \ WF (f ), then WF (u)contains the bicharacteristic through (x, ξ).
Did this by conjugating P (x,D) to model,
Q1(x,D) =∂
∂x1+ Q−∞(x,D)
whose Green’s function ∂∂x1
is H(x1)·δ(x′)+ . . .
![Page 42: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/42.jpg)
Current work: virtual ‘hybrid’ imaging in 2D
Only one kind of wave: electrostatic.
Good propagation of singularities is obtainednot via coupling with another physics,but by mathematical analysis.
![Page 43: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/43.jpg)
Current work: virtual ‘hybrid’ imaging in 2D
Only one kind of wave: electrostatic.
Good propagation of singularities is obtainednot via coupling with another physics,but by mathematical analysis.
After a transformation, singularities of D2Ndata propagate interior details efficientlyfrom any x0 ∈ Ω to any y0 ∈ ∂Ω.
Q.: What kind of PDE have this kind ofpropagation of singularities?
![Page 44: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/44.jpg)
Complex principal type (CPT) operators
pm(x, ξ) = pRm(x, ξ) + i pIm(x, ξ) with
(i) ∇x,ξ pR, ∇x,ξ pI linearly indep. at
Σ = (x, ξ) : pm(x, ξ) = 0 (codim 2)
(ii) Poisson bracket pRm, pIm :=
(∇ξ pRm) · (∇x pIm)− (∇x pRm) · (∇ξ pIm) ≡ 0 on Σ
(i), (ii) ⇐⇒ Σ is a codimension 2 coisotropicsubmanifold of T ∗X. =⇒
![Page 45: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/45.jpg)
Σ is foliated by 2-dim bicharacteristic leaves,which project to characteristic surfaces in X.
(iii) a nontrapping assumption.
Thm. (D.-H.) P (x,D) is locally solvable andif Pu = f , then
WF (u) \WF (f )
is a union of bicharacteristic leaves.
![Page 46: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/46.jpg)
Virtual hybrid edge detection: Exploit CPToperator structure underlying EIT to extractinformation about interior singularities of theconductivity.
Singularities propagate efficientlyalong 2D characteristics to ∂Ω.
Recently available CGO solutions have madenumerics doable.
![Page 47: Electrical Impedance Tomography and Hybrid Imaging · PDF fileElectrical Impedance Tomography and Hybrid Imaging ... July 15, 2016 Partially supported ... (RPT) operators P(x;D)](https://reader031.vdocument.in/reader031/viewer/2022020120/5a82131d7f8b9a682c8db441/html5/thumbnails/47.jpg)
Virtual hybrid edge detection: Exploit CPToperator structure underlying EIT to extractinformation about interior singularities of theconductivity.
Singularities propagate efficientlyalong 2D characteristics to ∂Ω.
Recently available CGO solutions have madenumerics doable.
Thank you!