Singular Values of the Attenuated Photoacoustic Imaging Operator
|
|
- Anabel Cole
- 5 years ago
- Views:
Transcription
1 Singular Values of the Attenuated Photoacoustic Imaging Operator Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 1 Computational Science Center, University Vienna, Austria & 2 Radon Institute of Computational and Applied Mathematics (RICAM), Linz, Austria Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 2 17, Radon 2017Institute 1 / of48com
2 Outline 1 Photoacoustics Applications 2 Modeling aspects Attenuation Examples of attenuation models 3 Inversion of the integrated photoacoustic operator Strongly attenuating media Weakly attenuating media 4 If time allows: Time reversal with variable sound speed (joint with Z. Belhachmi (Metz), T. Glatz (UVIE) Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 2 17, Radon 2017Institute 2 / of48com
3 Photoacoustics Photoacoustic imaging Lightning and Thunder (L.H. Wang) Specimen is uniformly illuminated by a short/pulsed electromagnetic pulse (visible or near infrared light - Photoacoustics, microwaves - Thermoacoustics) eter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 2 17, Radon 2017Institute 3 / of48com
4 Photoacoustics Photoacoustic imaging Lightning and Thunder (L.H. Wang) Specimen is uniformly illuminated by a short/pulsed electromagnetic pulse (visible or near infrared light - Photoacoustics, microwaves - Thermoacoustics) Two-step conversion process: Absorbed EM energy is converted into heat Material reacts with expansion Expansion produces pressure waves eter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 2 17, Radon 2017Institute 3 / of48com
5 Photoacoustics Photoacoustic imaging Lightning and Thunder (L.H. Wang) Specimen is uniformly illuminated by a short/pulsed electromagnetic pulse (visible or near infrared light - Photoacoustics, microwaves - Thermoacoustics) Two-step conversion process: Absorbed EM energy is converted into heat Material reacts with expansion Expansion produces pressure waves Imaging: Pressure waves are detected at the boundary of the object (over time) and are used for reconstruction of conversion parameter (EM energy into expansion/ waves) eter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 2 17, Radon 2017Institute 3 / of48com
6 Photoacoustics The method in between DOT (diffuse optical tomography), OCT (optical coherence tomography), SPIM (single plane imaging). MAD (mandibular advancement device), 2P (2 photon),... Fig. 1 Citation Steven Y. Leigh, Ye Chen, Jonathan T.C. Liu, "Modulated-alignment dual-axis (MAD) confocal microscopy for deep optical sectioning in tissues," Biomed. Opt. Express 5, (2014); Image 2014 Optical Society of America and may be used for noncommercial purposes only. Report a copyright concern regarding this image. Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 2 17, Radon 2017Institute 4 / of48com
7 Photoacoustics Schematic representation Acoustic Part Photoacoustic PAI Equations for wave propagation Thermodynamics Light Propagation QPAT Equations for photon density Figure: Three kind of problems Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 2 17, Radon 2017Institute 5 / of48com
8 Photoacoustics Governing equation of photoacoustics: Standard Wave equation for the pressure 1 c p dj (x, t) p(x, t) = t2 dt (t)h(x) Parameters and functions: H absorbed electromagnetic energy c 0 speed of sound - in most mathematical studies normalized to constant 1 j(t) can be considered to approximate a δ-impulse (Lightning) Alternative if j = δ: Initial value problem for homogeneous wave equation p(x, 0) = H(x), p t (x, 0) = 0 Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 2 17, Radon 2017Institute 6 / of48com
9 Photoacoustics Inverse problem of photoacoustics measurement data: p(x, t) for x M, t > 0 measurement region Reconstruction of H(x) for x Ω Exact Reconstruction Formulae for H depending on measurement geometry Sphere, Cylinder, Plane [Xu, Wang, 2002], [Finch, Patch, Rakesh, 2004] Circle [Finch, Haltmeier, Rakesh, 2005] Universal Backprojection [Wang et al, 2005, Natterer 2012] Palamodov Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 2 17, Radon 2017Institute 7 / of48com
10 Photoacoustics Applications Setup Figure: Photoacoustic setup: Robert Nuster and Günther Paltauf (University Graz) Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 2 17, Radon 2017Institute 8 / of48com
11 Photoacoustics Applications Applications: Microscopy Figure: Chicken embryo (5 days old)(1.2 mm) and mouse heart. Data by Berooz Zebian and Wolfgang Drexler, Robert Nuster et al. Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 2 17, Radon 2017Institute 9 / of48com
12 Photoacoustics Applications Medical applications Figure: Human studies of combined OCT (optical coherence tomography) and PAI: data by Boris Hermann and Wolfgang Drexler Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 10 / of48com
13 Modeling aspects Modeling aspects Standard Photoacoustics does not model variable sound speed, attenuation and variable illumination and does not recover physical parameters 1 Quantitative Photoacoustics (separation of parameters in H) [Arridge, Bal, Ken, Scotland, Uhlmann,...] 2 Sound speed variations: [Agranovsky, Hristova, Kuchment, Stefanov, Uhlmann,...] 3 Attenuation and dispersion [Anastasio, Patch, Riviere, Burgholzer, Kowar, S., Ammari, Wahab,...] 4 Variable Illumination [Wang, Bal,...] 5 Finite bandwidth detectors [Haltmeier et al] Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 11 / of48com
14 Modeling aspects Attenuation Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 12 / of48com
15 Modeling aspects Attenuation Attenuation: Model equation A κ p(x, t) p(x, t) = dδ dt (t)h(x) p(x, t) = 0 for t < 0 L: pseudo-differential operator w.r.t. t, with symbol κ 2 (ω) Iκ(ω) attenuation law Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 13 / of48com
16 Modeling aspects Attenuation Solution of forward problem Inverse Fourier-transform with respect to t: ˇf (ω) = 1 2π Solution of attenuated equation: ˇp(ω; x) = R 3 iω 4π 2π f (t)e iωt dt e iκ(ω) x y H(y)dy x y Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 14 / of48com
17 Modeling aspects Attenuation Solution of forward problem Inverse Fourier-transform with respect to t: ˇf (ω) = 1 2π Solution of attenuated equation: ˇp(ω; x) = R 3 iω 4π 2π f (t)e iωt dt e iκ(ω) x y H(y)dy x y Object of interest: Time-integrated photoacoustic operator in frequency domain: 1 R3 ˇP κ H(ω; x) = 4π e iκ(ω) x y H(y)dy 2π x y Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 14 / of48com
18 Modeling aspects Attenuation Formal definition of attenuation operator A κ : S (R R 3 ) S (R R 3 ) is defined by the action on the tensor products φ ψ S(R R 3 ), given by (φ ψ)(t, x) = φ(t)ψ(x): A κ u, φ ψ S,S = u, (F 1 κ 2 Fφ) ψ S,S, where F : S(R) S(R), Fφ(ω) = 1 2π φ(t)e iωt dt denotes the Fourier transform Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 15 / of48com
19 Modeling aspects Attenuation Assumptions on κ κ C (R; H) with H = {z C : Iz > 0} is called attenuation coefficient if 1 all derivatives of κ are polynomially bounded: l N0 κ1 >0,N N κ (l) (ω) κ 1 (1 + ω ) N 2 there exists a continuous continuation κ : H H of κ such that κ : H H is holomorphic. Moreover, κ1 >0,N N κ(z) κ 1 (1 + z )Ñ z H 3 κ( ω) = κ(ω) for all ω R Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 16 / of48com
20 Modeling aspects Attenuation Consequences of these assumptions 1 A κ is well-defined: κ 2 u S (R) u S (R) 2 p S (R R 3 ) is causal: that is supp(p) [0, ) R 3 3 A κ maps real-valued distributions onto itself Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 17 / of48com
21 Modeling aspects Attenuation Finite speed of propagation p S (R R 3 ) propagates with finite speed c > 0 if supp p {(t, x) R R 3 : x ct + R} h supp h B R (0) Lemma Let κ be an attenuation coefficient with the holomorphic extension κ : H H. Then, p propagates with finite speed if and only if κ(iω) lim ω iω > 0 In this case, it propagates with the speed c = lim ω iω κ(iω) Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 18 / of48com
22 Modeling aspects Examples of attenuation models Examples of attenuation models κ C (R; H) is called a strong attenuation coefficient if κ0,β>0,ω 0 0 Iκ(ω) κ 0 ω β, ω R, ω ω 0 weak attenuation coefficient if it is of the form c>0,κ 0 κ C (R) L 2 (R) κ(ω) = ω c + iκ + κ (ω), ω R Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 19 / of48com
23 Modeling aspects Examples of attenuation models Thermo-viscous model Name: Thermo-viscous model ω Coefficient: κ : R C, κ(ω) =, 1 iτω Upper bound: Speed: Type: Range of κ: κ(z) z for all z H iω c = lim ω κ(iω) = lim 1 + τω = ω Strong attenuation coefficient τ > 0, κ(z) = κ(z) i i κ i i Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 20 / of48com
24 Modeling aspects Examples of attenuation models Kowar-S-Bonnefond model Attenuation: Speed: Type: Range of κ: ( ) α κ : R C, κ(ω) = ω 1 +, 1+( iτω) γ κ(z) = κ(z) iω c = lim ω κ(iω) = lim ω 1 + Strong attenuation coefficient 1 α 1+(τω) γ = 1 i i κ i i Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 21 / of48com
25 Modeling aspects Examples of attenuation models Power law Szabo Model: Power law Coefficient: κ : R C, κ(ω) = ω + iα( iω) γ, κ(z) = κ(z) Parameters: γ (0, 1), α > 0 Upper bound: κ(z) z + α z γ α(1 γ) + (1 + αγ) z Propagation speed: iω c = lim ω κ(iω) = lim 1 ω 1 + αω γ 1 = 1 Attenuation type: Strong attenuation coefficient Range of κ: i i κ i i Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 22 / of48com
26 Modeling aspects Examples of attenuation models Modified Szabo s model Coefficient: κ : R C, κ(ω) = ω 1 + α( iω) γ 1 Parameters: γ (0, 1), α > 0 Upper bound: κ(z) 1 2 α(1 γ) + (1 + α 2 (1 + γ)) z iω Speed: c = lim ω κ(iω) = lim 1 = 1 ω 1 + αω γ 1 Type: Strong attenuation coefficient Range of κ: i i κ i i Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 23 / of48com
27 Modeling aspects Examples of attenuation models Nachman-Smith-Waag Coefficient: κ : R C, κ(ω) = ω 1 i τω c 0 1 iτω Parameters: c 0 > 0, τ > 0, τ (0, τ) z 1 i τz Extension: κ : H C, c 0 1 iτz, κ(z) 1 c 0 z for all z H iω 1 + τω τ Speed: c = lim ω κ(iω) = lim c 0 ω 1 + τω = c 0 τ Type: Weak attenuation coefficient Range of κ: i i κ i i Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 24 / of48com
28 Inversion of the integrated photoacoustic operator Photoacoustic equation in F-domain: recalled The solution p of the attenuated wave equation in Fourier domain: ˇp(ω, x) = G κ (ω, x y)h(y)dy R 3 = iω 4π 2π R3 e iκ(ω) x y x y Measurement data of photoacoustics H(y)dy, ω R, x R 3 ˇm(ω, ξ) = ˇp(ω, ξ) ω R, ξ Ω Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 25 / of48com
29 Inversion of the integrated photoacoustic operator Alternative equation: integrated photoacoustic operator Inverse problem of attenuated photoacoustics: Find H such that 1 1 ˇm(ω, ξ) = ˇp(ω, ξ), iω iω ω R, ξ Ω Meaning that the data are integrated in time-domain before inversion! Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 26 / of48com
30 Inversion of the integrated photoacoustic operator Integrated photoacoustic operator Let Ω R 3 be a bounded Lipschitz domain and Ω ε = {x Ω dist(x, Ω) > ε} κ denotes either a strong or a weak attenuation coefficient Definition ˇP κ : L 2 (Ω ε ) L 2 (R Ω), H 1 4π e 2π Ωε iκ(ω) ξ y H(y)dy ξ y is called integrated photoacoustic operator of the attenuation coefficient κ in frequency domain Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 27 / of48com
31 Inversion of the integrated photoacoustic operator Properties of integrated photoacoustic operator Let ˇP κ : L 2 (Ω ε ) L 2 (R Ω) be the integrated photoacoustic operator of a weak or a strong attenuation coefficient κ: Theorem Then, ˇP κ ˇP κ : L 2 (Ω ε ) L 2 (Ω ε ) is a self-adjoint integral operator with kernel F L 2 (Ω ε Ω ε ) given by that is F κ (x, y) = 1 32π 3 Ω e iκ(ω) ξ y iκ(ω) ξ x ds(ξ)dω, ξ y ξ x ˇP κ ˇP κ h(x) = F κ (x, y)h(y)dy Ω ε In particular, ˇP κ ˇP κ is a Hilbert Schmidt operator and thus compact Proof: lengthy Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 28 / of48com
32 Inversion of the integrated photoacoustic operator Strongly attenuating media Singular values of IPO for strongly attenuating media Kernel estimates 1 Lemma For strongly attenuation κ the kernel F κ of ˇP κ ˇP κ satisfies 1 B,b>0 sup sup j j! x,y Ω ε s j F κ (x, y + sv) Bbj j ( N β 1)j j N 0 v S 2 s=0 1 N N denotes the exponent for l = 0 (polynomial condition) and β (0, N] is the exponent in the condition for the strong attenuation coefficient κ Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 29 / of48com
33 Inversion of the integrated photoacoustic operator Strongly attenuating media Singular values of IPO for strongly attenuating media Let Ω be a bounded Lipschitz domain in R 3 and ε > 0 Corollary Let ˇP κ : L 2 (Ω ε ) L 2 (R Ω) be the integrated photoacoustic operator of a strong attenuation coefficient κ. Then, there exist constants C, c > 0 such that the eigenvalues (λ n ( ˇP κ ˇP κ )) n N (in decreasing order) satisfy λ n ( ˇP κ ˇP κ ) Cn m ) n exp ( cn β Nm n N Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 30 / of48com
34 Inversion of the integrated photoacoustic operator Weakly attenuating media Singular values of IPO for weakly attenuating media Split ˇP (0) κ h(ω, ξ) = ˇP (1) κ h(ω, ξ) = ˇP κ = ˇP (0) κ + ˇP (1) κ, 1 4π 2π 1 4π 2π Ωε e i ω c ξ y ξ y e κ ξ y h(y)dy Ωε e i ω c ξ y ξ y e κ ξ y (e iκ (ω) ξ y 1)h(y)dy ( ˇP κ (0), ˇP κ (1) ) photoacoustic operator with constant attenuation and the perturbation Lemma Let κ be a weak attenuation coefficient. Then, there exist constants C 1, C 2 > 0 such that we have C 1 n 2 3 λn ( ˇP κ ˇP κ ) C 2 n 2 3 n N Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 31 / of48com
35 Inversion of the integrated photoacoustic operator Weakly attenuating media Summary Medium type Examples Singular values decay No attenuation κ(ω) = ω c 0 n 2/3 (Palamodov) Weak attenuation Nachman-Smith-Waag κ(ω) = z 1 i τω c 0 1 iτω n 2/3 Strong attenuation Thermo-viscous model κ(ω) = ω 1 iτω e cnc Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 32 / of48com
36 Inversion of the integrated photoacoustic operator Weakly attenuating media References P. Elbau, O. Scherzer, and C. Shi. Singular values of the attenuated photoacoustic imaging operator. Preprint on ArXiv arxiv: , University of Vienna, Austria, R. Kowar and O. Scherzer. Attenuation models in photoacoustics. In H. Ammari, editor, Mathematical Modeling in Biomedical Imaging II: Optical, Ultrasound, and Opto-Acoustic Tomographies. Springer Verlag, Berlin Heidelberg, R. Kowar, O. Scherzer, and X. Bonnefond. Causality analysis of frequency-dependent wave attenuation. Math. Methods Appl. Sci., 34: , A. I. Nachman, J. F. Smith, III, and R. C. Waag. An equation for acoustic propagation in inhomogeneous media with relaxation losses. J. Acoust. Soc. Amer., 88(3): , V. P. Palamodov. Remarks on the general Funk transform and thermoacoustic tomography. Inverse Probl. Imaging, 4(4): , T.L. Szabo. Time domain wave equations for lossy media obeying a frequency power law. J. Acoust. Soc. Amer., 96: , Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 33 / of48com
37 Inversion of the integrated photoacoustic operator Weakly attenuating media Sound Speed Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 34 / of48com
38 If time allows: Time reversal with variable sound speed (joint with Z. Belhachmi (Metz), T. Glatz (UVIE) Photoacoustic with variable sound speed Reconstruct H in from measurements of 1 c 2 (x) p (x, t) p(x, t) = 0 in R n (0, ) p(x, 0) = H(x) in R n p (x, 0) = 0 in R n m = p Γ (0, ) Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 35 / of48com
39 If time allows: Time reversal with variable sound speed (joint with Z. Belhachmi (Metz), T. Glatz (UVIE) Literature Exact reconstruction formulas Convergence results require non-trapping conditions M. Agranovsky and P. Kuchment. Uniqueness of reconstruction and an inversion procedure for thermoacoustic and photoacoustic tomography with variable sound speed. Inverse Probl., 23(5): , P. Stefanov and G. Uhlmann. Thermoacoustic tomography with variable sound speed. Inverse Probl., 25(7):075011, 16, 2009 Y. Hristova, P. Kuchment, and L. Nguyen. Reconstruction and time reversal in thermoacoustic tomography in acoustically homogeneous and inhomogeneous media. Inverse Probl., 24(5):055006, Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 36 / of48com
40 If time allows: Time reversal with variable sound speed (joint with Z. Belhachmi (Metz), T. Glatz (UVIE) Our approach: Operator equation formulation Let L : H 1 0(Ω) L 2 (Γ (0, )) H p Γ (0, ) We consider the solution of L[H] = m with Landweber iteration H k+1 = H k λl [L[H k ] m] Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 37 / of48com
41 If time allows: Time reversal with variable sound speed (joint with Z. Belhachmi (Metz), T. Glatz (UVIE) Adjoint problem Let h C (Γ (0, )). Then the L 2 -adjoint is given z := z(h) = L 2 h, where z solves 1 c 2 z z = 0 in R n \ Γ (0, T ) z(t ) = z (T ) = 0 in R n [z] = 0 on Γ (0, T ) [ ν z] = h on Γ (0, T ) Moreover, [ z L h = 1 ] (0) c 2 Z. Belhachmi, T. Glatz, O.Scherzer A Direct Method for Photoacoustic Tomography with Inhomogeneous Sound Speed. Inverse Problems 4, , 2016 Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 38 / of48com
42 If time allows: Time reversal with variable sound speed (joint with Z. Belhachmi (Metz), T. Glatz (UVIE) Classical regularization theory Let H be the minimum-norm-solution H { } H1 = min H H1 : L[H] = m Then 1 If m = p(h ). Then H k H. 2 If m δ m δ, then H δ k (δ) H, where k (δ) is the largest index such that p(hk δ ) (δ) mδ > δ C. W. Groetsch. The Theory of Tikhonov Regularization for Fredholm Equations of the First Kind. Pitman, Boston, 1984 Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 39 / of48com
43 If time allows: Time reversal with variable sound speed (joint with Z. Belhachmi (Metz), T. Glatz (UVIE) Time reversal Let h L 2 (Γ (0, )) and Γ = Ω. Moreover, let ˆL[h] be defined z c 2 z = 0 in Ω (0, T ) z(t ) = 0, z (T ) = 0 in R n z = h on Γ (0, T ) Standard time reversal: Ĥ = ˆL[m] M. Fink. Time reversal of ultrasonic fields. i. basic principles. IEEE Trans. Ultrason., Ferroeletr., Freq. Control, 39(5): , H. Grün, G. Paltauf, M. Haltmeier, and P. Burgholzer. Photoacoustic tomography using a fiber based Fabry-Perot interferometer as an integrating line detector and image reconstruction by model-based time reversal method. Volume 6631 of Proceedings of SPIE, page SPIE, 2007 Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 40 / of48com
44 If time allows: Time reversal with variable sound speed (joint with Z. Belhachmi (Metz), T. Glatz (UVIE) Time reversal: Stefanov & Uhlmann Let h L 2 (Γ (0, )) and Γ = Ω. Moreover, let L[h] = z(0) be defined z c 2 z = 0 in Ω (0, T ) z(t ) = φ, z (T ) = 0 in R n z = h on Γ (0, T ) where φ = 0 in Ω and φ Γ = h(t ) on Γ Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 41 / of48com
45 If time allows: Time reversal with variable sound speed (joint with Z. Belhachmi (Metz), T. Glatz (UVIE) Time reversal of Stefanov & Uhlmann First iteration (time reversal) H = L[m] Moreover, LL = I K (K compact) and H = K m Lm m=0 Probably not minimum norm solution. P. Stefanov and G. Uhlmann. Thermoacoustic tomography with variable sound speed. Inverse Probl., 25(7):075011, 16, 2009 Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 42 / of48com
46 If time allows: Time reversal with variable sound speed (joint with Z. Belhachmi (Metz), T. Glatz (UVIE) Non-trapping speed, T = 4T 0 Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 43 / of48com
47 If time allows: Time reversal with variable sound speed (joint with Z. Belhachmi (Metz), T. Glatz (UVIE) Non-trapping speed, T = 1.2T 0 Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 44 / of48com
48 If time allows: Time reversal with variable sound speed (joint with Z. Belhachmi (Metz), T. Glatz (UVIE) Landweber: Non-trapping speed, T = 1.2T 0 Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 45 / of48com
49 If time allows: Time reversal with variable sound speed (joint with Z. Belhachmi (Metz), T. Glatz (UVIE) Trapping speed, T = 2T 0. Time reversal; 5th Neumann sum; Landweber iterate 20 Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 46 / of48com
50 If time allows: Time reversal with variable sound speed (joint with Z. Belhachmi (Metz), T. Glatz (UVIE) Generalizations κ(x)y (x, t) (ρ(x) 1 y(x, t) ) = 0 in R n (0, ) y(x, 0) = H(x) in R n y (x, 0) = 0 in R n. κ compressibility ρ density c(x) = 1 κρ Z. Belhachmi, T. Glatz, O.Scherzer Photoacoustic Tomography With Spatially Varying Compressibility and Density. JIIPP, 2016 Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 47 / of48com
51 If time allows: Time reversal with variable sound speed (joint with Z. Belhachmi (Metz), T. Glatz (UVIE) Thank you for your attention Peter Elbau 1, Otmar Scherzer 1,2, Cong Shi 1 Singular ( 1 Computational Values of the Science Attenuated Center, Photoacoustic University Vienna, Imaging Austria February Operator & 17, 2 Radon 2017 Institute 48 / of48com
Bayesian approach to image reconstruction in photoacoustic tomography
Bayesian approach to image reconstruction in photoacoustic tomography Jenni Tick a, Aki Pulkkinen a, and Tanja Tarvainen a,b a Department of Applied Physics, University of Eastern Finland, P.O. Box 167,
More informationTHERMO AND PHOTOACOUSTIC TOMOGRAPHY WITH VARIABLE SPEED AND PLANAR DETECTORS
THERMO AND PHOTOACOUSTIC TOMOGRAPHY WITH VARIABLE SPEED AND PLANAR DETECTORS PLAMEN STEFANOV AND YANG YANG Abstract. We analyze the mathematical model of multiwave tomography with a variable speed with
More informationAppliedMathematics. Iterative Methods for Photoacoustic Tomography in Attenuating Acoustic Media. Markus Haltmeier, Richard Kowar and Linh V.
Nr. 37. May 7 Leopold-Franzens-Universität Innsbruck Preprint-Series: Department of Mathematics - Applied Mathematics Iterative Methods for Photoacoustic Tomography in Attenuating Acoustic Media Markus
More informationarxiv: v2 [math.na] 9 Aug 2017
Iterative Methods for Photoacoustic Tomography in Attenuating Acoustic Media Markus Haltmeier Richard Kowar Linh V. Nguyen arxiv:74.74v [math.na] 9 Aug 7 Department of Mathematics, University of Innsbruck
More informationA one-step reconstruction algorithm for quantitative photoacoustic imaging
A one-step reconstruction algorithm for quantitative photoacoustic imaging Tian Ding Kui Ren Sarah Vallélian Abstract Quantitative photoacoustic tomography (QPAT) is a recent hybrid imaging modality that
More informationQuantitative thermo-acoustics and related problems
Quantitative thermo-acoustics and related problems TING ZHOU MIT joint work with G. Bal, K. Ren and G. Uhlmann Conference on Inverse Problems Dedicated to Gunther Uhlmann s 60th Birthday UC Irvine June
More informationQuantitative Thermo-acoustics and related problems
Quantitative Thermo-acoustics and related problems Guillaume Bal Department of Applied Physics & Applied Mathematics, Columbia University, New York, NY 10027 E-mail: gb2030@columbia.edu Kui Ren Department
More informationOn reconstruction and time reversal in thermoacoustic tomography in acoustically homogeneous and inhomogeneous media
On reconstruction and time reversal in thermoacoustic tomography in acoustically homogeneous and inhomogeneous media Yulia Hristova, Peter Kuchment, and Linh Nguyen Mathematics Department, Texas A& M University,
More informationPhoto-Acoustic imaging in layered media
Photo-Acoustic imaging in layered media Faouzi TRIKI Université Grenoble-Alpes, France (joint works with Kui Ren, Texas at Austin, USA.) 4 juillet 2016 Photo-Acoustic effect Photo-Acoustic effect (Bell,
More informationAppliedMathematics. Operator learning approach for the limited view problem in photoacoustic tomography
Nr. 38 09. May 2017 Leopold-Franzens-Universität Innsbruck Preprint-Series: Department of Mathematics - Applied Mathematics Operator learning approach for the limited view problem in photoacoustic tomography
More informationOverview I. Approaches to Reconstruction in Photoacoustic Tomography. Bell and the Photoacoustic Effect I. Overview II
I Introduction Photoacoustic Tomography Approaches to Reconstruction in Photoacoustic Tomography March 9, 2011 Photoacoustic tomography (PAT) and the related thermoacoustic tomography (TAT) are important
More informationA Neumann series based method for photoacoustic tomography on irregular domains
A Neumann series based method for photoacoustic tomography on irregular domains Eric Chung, Chi Yeung Lam, and Jianliang Qian Abstract. Recently, a Neumann series based numerical method is developed for
More informationA one-step reconstruction algorithm for quantitative photoacoustic imaging
Inverse Problems PAPER A one-step reconstruction algorithm for quantitative photoacoustic imaging To cite this article: Tian Ding et al 0 Inverse Problems 000 Manuscript version: Accepted Manuscript Accepted
More informationRecent progress on the explicit inversion of geodesic X-ray transforms
Recent progress on the explicit inversion of geodesic X-ray transforms François Monard Department of Mathematics, University of Washington. Geometric Analysis and PDE seminar University of Cambridge, May
More informationInjectivity and Exact Inversion of Ultrasound Operators in the Sph
Injectivity and Exact Inversion of Ultrasound Operators in the Spherical Geometry University of Texas at Arlington Department of Mathematics Geometric Analysis on Euclidean and Homogeneous Spaces January
More informationA Second-order Finite Difference Scheme For The Wave Equation on a Reduced Polar Grid DRAFT
A Second-order Finite Difference Scheme For The Wave Equation on a Reduced Polar Grid Abstract. This paper presents a second-order numerical scheme, based on finite differences, for solving the wave equation
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 informationAppliedMathematics. Analysis of the Linearized Problem of Quantitative Photoacoustic Tomography
Nr. 3 6. February 27 Leopold-Franzens-Universität Innsbruck Preprint-Series: Department of Mathematics - Applied Mathematics Analysis of the Linearized Problem of Quantitative Photoacoustic Tomography
More informationReference: Cheney M. and Borden B., Fundamentals of Radar Imaging, SIAM, 2009
Inverse Problems Inverse Problems are problems where causes for a desired or an observed effect are to be determined. They lie at the heart of scientific inquiry and technological development. Applications
More informationRECONSTRUCTION OF MICROWAVE ABSORPTION PROPERTIES IN HETEROGENEOUS TISSUE FOR MICROWAVE-INDUCED THERMO-ACOUSTIC TOMOG- RAPHY
Progress In Electromagnetics Research, Vol. 130, 225 240, 2012 RECONSTRUCTION OF MICROWAVE ABSORPTION PROPERTIES IN HETEROGENEOUS TISSUE FOR MICROWAVE-INDUCED THERMO-ACOUSTIC TOMOG- RAPHY J.-G. Wang 1,
More informationMultiwave imaging in an enclosure with variable wave speed
Multiwave imaging in an enclosure with variable wave speed Sebastián Acosta 1 and Carlos Montalto 2 1 Department of Pediatrics Cardiology, Baylor College of Medicine, TX, USA 2 Department of Mathematics,
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 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 informationA global stability estimate for the photo-acoustic inverse problem in layered media
A global stability estimate for the photo-acoustic inverse problem in layered media Kui Ren Faouzi Triki May 11, 217 Abstract This paper is concerned with the stability issue in determining absorption
More informationA hybrid reconstruction method for quantitative PAT
A hybrid reconstruction method for quantitative PAT Kui Ren Hao Gao Hongkai Zhao January 20, 2012 Abstract The objective of quantitative photoacoustic tomography (qpat) is to reconstruct the diffusion
More informationHyperbolic inverse problems and exact controllability
Hyperbolic inverse problems and exact controllability Lauri Oksanen University College London An inverse initial source problem Let M R n be a compact domain with smooth strictly convex boundary, and let
More informationINSTABILITY OF THE LINEARIZED PROBLEM IN MULTIWAVE TOMOGRAPHY OF RECOVERY BOTH THE SOURCE AND THE SPEED
INSTABILITY OF THE LINEARIZED PROBLEM IN MULTIWAVE TOMOGRAPHY OF RECOVERY BOTH THE SOURCE AND THE SPEED PLAMEN STEFANOV AND GUNTHER UHLMANN ABSTRACT. In this paper we consider the linearized problem of
More informationShort note on compact operators - Monday 24 th March, Sylvester Eriksson-Bique
Short note on compact operators - Monday 24 th March, 2014 Sylvester Eriksson-Bique 1 Introduction In this note I will give a short outline about the structure theory of compact operators. I restrict attention
More informationThe Levenberg-Marquardt Iteration for Numerical Inversion of the Power Density Operator
The Levenberg-Marquardt Iteration for Numerical Inversion of the Power Density Operator G. Bal (gb2030@columbia.edu) 1 W. Naetar (wolf.naetar@univie.ac.at) 2 O. Scherzer (otmar.scherzer@univie.ac.at) 2,3
More informationInversion of the circular Radon transform on an annulus
Inversion of the circular Radon transform on an annulus Gaik Ambartsoumian 1, Rim Gouia-Zarrad 1 and Matthew A Lewis 2 1 Department of Mathematics, University of Texas at Arlington, Arlington, TX 7617-48,
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 Methods for geodesic X-ray transforms and applications to open theoretical questions
Numerical Methods for geodesic X-ray transforms and applications to open theoretical questions François Monard Department of Mathematics, University of Washington. Nov. 13, 2014 UW Numerical Analysis Research
More informationRADAR Imaging A Mathematical Perspective
RADAR Imaging A Mathematical Perspective P. ELBAU (UNIVERSITY OF VIENNA) AND O. SCHERZER (UNIVERSITY OF VIENNA AND RICAM, AUSTRIAN ACADEMY OF SCIENCES, LINZ) In this paper we explain with simple geometric
More informationThe spherical mean value operator with centers on a sphere
The spherical mean value operator with centers on a sphere David Finch Department of Mathematics Oregon State University Corvallis, OR 9733-465 Email: finch@math.oregonstate.edu Rakesh Department of Mathematics
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 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 informationA new method for the solution of scattering problems
A new method for the solution of scattering problems Thorsten Hohage, Frank Schmidt and Lin Zschiedrich Konrad-Zuse-Zentrum Berlin, hohage@zibde * after February 22: University of Göttingen Abstract We
More informationA method of the forward problem for magneto-acousto-electrical tomography
Technology and Health Care 24 (216) S733 S738 DOI 1.3233/THC-16122 IOS Press S733 A method of the forward problem for magneto-acousto-electrical tomography Jingxiang Lv a,b,c, Guoqiang Liu a,, Xinli Wang
More informationA hybrid reconstruction method for quantitative PAT
A hybrid reconstruction method for quantitative PAT Kui Ren Hao Gao Hongkai Zhao October 24, 2012 Abstract The objective of quantitative photoacoustic tomography (qpat) is to reconstruct the diffusion
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 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 informationDetecting Interfaces in a Parabolic-Elliptic Problem
Detecting Interfaces in a Parabolic-Elliptic Problem Bastian Gebauer bastian.gebauer@oeaw.ac.at Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy of Sciences, Linz,
More informationNonlinear error dynamics for cycled data assimilation methods
Nonlinear error dynamics for cycled data assimilation methods A J F Moodey 1, A S Lawless 1,2, P J van Leeuwen 2, R W E Potthast 1,3 1 Department of Mathematics and Statistics, University of Reading, UK.
More informationMicrowave-induced thermoacoustic tomography using multi-sector scanning
Microwave-induced thermoacoustic tomography using multi-sector scanning Minghua Xu, Geng Ku, and Lihong V. Wang a) Optical Imaging Laboratory, Biomedical Engineering Program, Texas A&M University, 3120
More information19 Mathematics of Photoacoustic and Thermoacoustic Tomography
19 Mathematics of Photoacoustic and Thermoacoustic Tomography Peter Kuchment Leonid Kunyansky 19.1 Introduction...819 19.2 MathematicalModels oftat...820 19.2.1 Point Detectors andthe WaveEquation Model...820
More informationReconstruction of the refractive index by near- eld phaseless imaging
Reconstruction of the refractive index by near- eld phaseless imaging Victor Palamodov November 0, 207 Abstract. An explicit method is described for reconstruction of the complex refractive index by the
More informationInverse Transport Problems (Generalized) Stability Estimates. Guillaume Bal
Inverse Transport Problems (Generalized) Stability Estimates Guillaume Bal Department of Applied Physics & Applied Mathematics Columbia University Joint with Alexandre Jollivet Data Acquisition in CT-scan
More informationPart 2 Introduction to Microlocal Analysis
Part 2 Introduction to Microlocal Analysis Birsen Yazıcı & Venky Krishnan Rensselaer Polytechnic Institute Electrical, Computer and Systems Engineering August 2 nd, 2010 Outline PART II Pseudodifferential
More informationExact Reconstruction in Photoacoustic Tomography with Circular Integrating Detectors II: Spherical Geometry
Exact Reconstruction in Photoacoustic Tomography with Circular Integrating Detectors II: Spherical Geometry Gerhard Zangerl 1 Otmar Scherzer 1,2 1 Department of Mathematics 2 Radon Institute of Computational
More informationAnalysis of the anomalous localized resonance
Hyundae Lee(Inha University, Korea) Joint work with Habib Ammari, Giulio Ciraolo, Hyeonbae Kang, Graeme Milton. UCI, June 22, 2012 A conference on inverse problems in honor of Gunther Uhlmann Outline Introduction
More informationMathematics of Photoacoustic and Thermoacoustic Tomography
Mathematics of Photoacoustic and Thermoacoustic Tomography arxiv:0912.2022v1 [math.ap] 10 Dec 2009 P. Kuchment and L. Kunyansky Contact information: Peter Kuchment Mathematics Department Texas A & M University
More informationRegularizations of general singular integral operators
Regularizations of general singular integral operators Texas A&M University March 19th, 2011 This talk is based on joint work with Sergei Treil. This work is accepted by Revista Matematica Iberoamericano
More informationNonlinear Flows for Displacement Correction and Applications in Tomography
Nonlinear Flows for Displacement Correction and Applications in Tomography Guozhi Dong 1 and Otmar Scherzer 1,2 1 Computational Science Center, University of Vienna, Oskar-Morgenstern-Platz 1, 1090 Wien,
More informationA time domain probe method for threedimensional rough surface reconstructions
A time domain probe method for threedimensional rough surface reconstructions Article Accepted Version Burkard, C. and Potthast, R. (2009) A time domain probe method for three dimensional rough surface
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Error Estimates for Filtered Back Projection Matthias Beckmann and Armin Iske Nr. 2015-03 January 2015 Error Estimates for Filtered Back Projection Matthias
More informationA series solution and a fast algorithm for the inversion of the spherical mean Radon transform
IOP PUBLISHING Inverse Problems 23 (2007) S11 S20 INVERSE PROBLEMS doi:10.1088/0266-5611/23/6/s02 A series solution and a fast algorithm for the inversion of the spherical mean Radon transform Leonid A
More informationAmbiguity of optical coherence tomography measurements due to rough surface scattering
Ambiguity of optical coherence tomography measurements due to rough surface scattering Y. Ashtamker, 1 V Freilikher, 1,* and J C Dainty 2 1 Department of Physics, Bar-Ilan University, Ramat Gan 52900,
More informationInversions of ray transforms on simple surfaces
Inversions of ray transforms on simple surfaces François Monard Department of Mathematics, University of Washington. June 09, 2015 Institut Henri Poincaré - Program on Inverse problems 1 / 42 Outline 1
More informationA SHARP STABILITY ESTIMATE IN TENSOR TOMOGRAPHY
A SHARP STABILITY ESTIMATE IN TENSOR TOMOGRAPHY PLAMEN STEFANOV 1. Introduction Let (M, g) be a compact Riemannian manifold with boundary. The geodesic ray transform I of symmetric 2-tensor fields f is
More informationExplicit Reconstructions in QPAT, QTAT, TE, and MRE
Explicit Reconstructions in QPAT, QTAT, TE, and MRE Guillaume Bal February 14, 2012 Abstract Photo-acoustic Tomography (PAT) and Thermo-acoustic Tomography (TAT) are medical imaging modalities that combine
More informationarxiv: v3 [physics.class-ph] 27 Jan 2014
Comparison of fractional wave equations for power law attenuation in ultrasound and elastography Sverre Holm a, Sven Peter Näsholm b a Department of Informatics, University of Oslo, P. O. Box 18, NO 316
More informationOn stable inversion of the attenuated Radon transform with half data Jan Boman. We shall consider weighted Radon transforms of the form
On stable inversion of the attenuated Radon transform with half data Jan Boman We shall consider weighted Radon transforms of the form R ρ f(l) = f(x)ρ(x, L)ds, L where ρ is a given smooth, positive weight
More informationContinuous Frames and Sampling
NuHAG University of Vienna, Faculty for Mathematics Marie Curie Fellow within the European network HASSIP HPRN-CT-2002-285 SampTA05, Samsun July 2005 Joint work with Massimo Fornasier Overview 1 Continuous
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 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 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 informationRecovering a function from its spherical mean values in two and three dimensions
Recovering a function from its spherical mean values in two and three dimensions David Finch Rakesh 1 Introduction The problem of determining a function from a subset of its spherical means has a rich
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 informationStatistical Inverse Problems and Instrumental Variables
Statistical Inverse Problems and Instrumental Variables Thorsten Hohage Institut für Numerische und Angewandte Mathematik University of Göttingen Workshop on Inverse and Partial Information Problems: Methodology
More informationAN IMAGE CORRECTION METHOD BASED ON ELEC- TROMAGNETIC SIMULATION FOR MICROWAVE IN- DUCED THERMO-ACOUSTIC TOMOGRAPHY SYSTEM
Progress In Electromagnetics Research B, Vol. 43, 19 33, 2012 AN IMAGE CORRECTION METHOD BASED ON ELEC- TROMAGNETIC SIMULATION FOR MICROWAVE IN- DUCED THERMO-ACOUSTIC TOMOGRAPHY SYSTEM J. Song 1, Z. Q.
More informationMathematical analysis of ultrafast ultrasound imaging
Mathematical analysis of ultrafast ultrasound imaging Giovanni S Alberti Department of Mathematics, University of Genoa Joint work with H. Ammari (ETH), F. Romero (ETH) and T. Wintz (Sony). IPMS 18, May
More informationMicrolocal Methods in X-ray Tomography
Microlocal Methods in X-ray Tomography Plamen Stefanov Purdue University Lecture I: Euclidean X-ray tomography Mini Course, Fields Institute, 2012 Plamen Stefanov (Purdue University ) Microlocal Methods
More informationThe structure of laser pulses
1 The structure of laser pulses 2 The structure of laser pulses Pulse characteristics Temporal and spectral representation Fourier transforms Temporal and spectral widths Instantaneous frequency Chirped
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 informationThe Method of Small Volume Expansions and Applications to Emerging Medical Imaging
The Method of Small Volume Expansions and Applications to Emerging Medical Imaging Hyeonbae Kang hbkang@inha.ac.kr people.math.inha.ac.kr/ hbkang Inha University Oct. 14, 2009, KAIST p. 1/5 Contents Method
More informationTwo-parameter regularization method for determining the heat source
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 8 (017), pp. 3937-3950 Research India Publications http://www.ripublication.com Two-parameter regularization method for
More informationParameter Identification in Partial Differential Equations
Parameter Identification in Partial Differential Equations Differentiation of data Not strictly a parameter identification problem, but good motivation. Appears often as a subproblem. Given noisy observation
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
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 informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationPart 2 Introduction to Microlocal Analysis
Part 2 Introduction to Microlocal Analysis Birsen Yazıcı& Venky Krishnan Rensselaer Polytechnic Institute Electrical, Computer and Systems Engineering March 15 th, 2010 Outline PART II Pseudodifferential(ψDOs)
More information10. OPTICAL COHERENCE TOMOGRAPHY
1. OPTICAL COHERENCE TOMOGRAPHY Optical coherence tomography (OCT) is a label-free (intrinsic contrast) technique that enables 3D imaging of tissues. The principle of its operation relies on low-coherence
More informationModeling power law absorption and dispersion for acoustic propagation using the fractional Laplacian
Modeling power law absorption and dispersion for acoustic propagation using the fractional Laplacian Bradley E. Treeby a and B. T. Cox Department of Medical Physics and Bioengineering, University College
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 informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationCONTINUITY OF CP-SEMIGROUPS IN THE POINT-STRONG OPERATOR TOPOLOGY
J. OPERATOR THEORY 64:1(21), 149 154 Copyright by THETA, 21 CONTINUITY OF CP-SEMIGROUPS IN THE POINT-STRONG OPERATOR TOPOLOGY DANIEL MARKIEWICZ and ORR MOSHE SHALIT Communicated by William Arveson ABSTRACT.
More informationConvergence of the Ensemble Kalman Filter in Hilbert Space
Convergence of the Ensemble Kalman Filter in Hilbert Space Jan Mandel Center for Computational Mathematics Department of Mathematical and Statistical Sciences University of Colorado Denver Parts based
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More informationDunkl operators and Clifford algebras II
Translation operator for the Clifford Research Group Department of Mathematical Analysis Ghent University Hong Kong, March, 2011 Translation operator for the Hermite polynomials Translation operator for
More informationVelocity averaging a general framework
Outline Velocity averaging a general framework Martin Lazar BCAM ERC-NUMERIWAVES Seminar May 15, 2013 Joint work with D. Mitrović, University of Montenegro, Montenegro Outline Outline 1 2 L p, p >= 2 setting
More informationOn the hydrodynamic diffusion of rigid particles
On the hydrodynamic diffusion of rigid particles O. Gonzalez Introduction Basic problem. Characterize how the diffusion and sedimentation properties of particles depend on their shape. Diffusion: Sedimentation:
More informationQuantum Optical Coherence Tomography
Quantum Optical Coherence Tomography Bahaa Saleh Alexander Sergienko Malvin Teich Quantum Imaging Lab Department of Electrical & Computer Engineering & Photonics Center QuickTime and a TIFF (Uncompressed)
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 informationOn the split equality common fixed point problem for quasi-nonexpansive multi-valued mappings in Banach spaces
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 (06), 5536 5543 Research Article On the split equality common fixed point problem for quasi-nonexpansive multi-valued mappings in Banach spaces
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 informationViscosity Iterative Approximating the Common Fixed Points of Non-expansive Semigroups in Banach Spaces
Viscosity Iterative Approximating the Common Fixed Points of Non-expansive Semigroups in Banach Spaces YUAN-HENG WANG Zhejiang Normal University Department of Mathematics Yingbing Road 688, 321004 Jinhua
More informationAppliedMathematics. Iterative Methods for Photoacoustic Tomography with Variable Sound Speed. Markus Haltmeier, Linh V. Nguyen
Nr. 29 7. December 26 Leopold-Franzens-Universität Innsbruck Preprint-Series: Department of Mathematics - Applied Mathematics Iterative Methods for Photoacoustic Tomography with Variable Sound Speed Markus
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 information