The broken-ray transform and its generalizations
|
|
- Angel Booker
- 5 years ago
- Views:
Transcription
1 The broken-ray transform and its generalizations Gaik Ambartsoumian University of Texas at Arlington 100 Years of the Radon Transform Linz, March 27-31, 2017
2 Acknowledgements The talk is based on results of collaborative work with Rim Gouia-Zarrad, American University of Sharjah, UAE Mohammad Latifi-Jebelli, University of Texas at Arlington Sunghwan Moon, UNIST, Korea Souvik Roy, Universität Würzburg, Germany Partially supported by NSF DMS Simons Foundation
3 Outline Some Motivating Imaging Modalities Prior Work and Terminology BRT in a Disc BRT in a Slab Geometry CRT in a Slab Geometry Partial Order and Positive Cones Cone Differentiation BRT in a Slab Geometry (revisited) Polyhedral CRT in a Slab Geometry Radon-Nikodym Theorem and Range Description
4 Single Scattering Optical Tomography (SSOT) Uses light, transmitted and scattered through an object, to determine the interior features of that object. If the object has moderate optical thickness it is reasonable to assume the majority of photons scatter once. Using collimated emitters/receivers one can measure the intensity of light scattered along various broken rays. Need to recover the spatially varying coefficients of light absorption and/or light scattering.
5 Florescu, Schotland and Markel (2009, 2010, 2011) So if the scattering coefficient is known, then the reconstruction of the absorption coefficient is reduced to inversion of a generalized Radon transform integrating along the broken rays.
6 Compton Scattering cos β = 1 mc2 E (E E)E R. Basko, G. Zeng, and G. Gullberg (1997, 1998) M. Nguyen, T. Truong, et al (2000 s)
7 Broken-Ray and Conical Transforms L. Florescu, V. Markel, J. Schotland, F. Zhao P. Grangeat, M. Morvidone, M. Nguyen, R. Régnier, T. Truong, H. Zaidi A. Katsevich, R. Krylov F. Terzioglu V. Palamodov R. Gouia-Zarrad M. Courdurier, F. Monard, A. Osses and F. Romero D. Finch, B. Sherson
8 Broken-Ray and Conical Transforms M. Cree, P. Bones and R. Basko, G. Zeng, G. Gullberg M. Allmaras, D. Darrow, Y. Hristova, G. Kanschat, P. Kuchment M. Haltmeier C. Jung, S. Moon V. Maxim D. Schiefeneder M. Lassas, M. Salo, G. Uhlmann M. Hubenthal, J. Ilmavirta
9 V-line Radon Transform (VRT) in 2D Definition The V-line Radon transform of function f (x, y) is the integral Rf (β, t) = f ds, (1) BR(β,t) of f along the broken ray BR(β, t) with respect to line measure ds. The problem of inversion is over-determined, so it is natural to consider a restriction of Rf to a two-dimensional set.
10 Geometry: Slab vs Disc Available directions Stability of reconstruction Hardware implementation (?)
11 Full Data Theorem If f (x, y) is a smooth function supported in the disc D(0, R sin θ), then f is uniquely determined by Rf (φ, d), φ [0, 2π], d [0, 2R].
12 Inversion Formula Rf (ψ φ, t d ) = Rf (φ, d) + Rf (φ + π, 2R d) Rf (φ, 2R), (2) for all values φ [0, 2π] and d [0, 2R]. f (x, y) = 1 4π 2π 0 H ( Rf t ) (ψ, x cos ψ + y sin ψ) dψ (3) where H is the Hilbert transform defined by Hh (t) = i sgn (r) ĥ(r) eirt dr. (4) 2π R and ĥ(r) is the Fourier transform of h(t), i.e. ĥ(r) = 1 h(t) e irt dt. (5) 2π R
13 Inversion Formula Issues with the support Interior problem Other methods without loss of information Rotation invariance
14 VRT in a Disc: Partial Data (G.A., S. Moon 2013) Theorem If f (x, y) is a smooth function supported in the disc D(0, R), then f is uniquely determined by Rf (φ, d), φ [0, 2π], d [0, R].
15 Fourier Expansions Denote g(β, t) := Rf (β, t). f (φ, ρ) = n= f n (ρ) e inφ, g(β, t) = where the Fourier coefficients are given by f n (ρ) = 1 2π 2π 0 n= f (φ, ρ) e inφ dφ, g n (t) = 1 2π 2π 0 g n (t) e inβ, g(β, t) e inβ dβ.
16 Inversion Formula Mf n (s) = Mg n (s 1), R(s) > 1 (6) 1/(s 1) + Mh n (s 1) where MF denotes the Mellin transform of function F MF (s) = 0 p s 1 F (p) dp, and h n is some fixed function. Hence for any t > 1 we have 1 t+ti f n (ρ) = lim ρ s Mg n (s 1) ds. (7) T 2πi t Ti 1/(s 1) + Mh n (s 1)
17 Definition of h n If 1 < t < 1 sin θ then 1 + t cos[ψ(t)] + t 2 sin[ψ(t)] h n (t) = ( 1) n e inψ(t) 1 + t 2 + 2t cos(ψ(t)) sin θ 1 t 2 sin 2 θ 1 t cos[2θ ψ(t)] + t 2 sin[2θ ψ(t)] e in[2θ ψ(t)] 1 + t 2 2t cos[2θ ψ(t)] sin θ 1 t 2 sin 2 θ, 1 + t cos[ψ(t)] + t 2 sin[ψ(t)] h n (t) = ( 1) n e inψ(t) 1 + t 2 + 2t cos[ψ(t)] sin θ 1 t 2 sin 2 θ, 0 < t 1 and h n (t) 0, for all t > 1 sin θ. Here ψ(t) = arcsin(t sin θ) + θ.
18 Numerical Reconstruction (G.A., S. Roy 2015)
19 Numerical Reconstruction (G.A., S.(b) rroy (a) r = M/8 = M/ ) Fig. 7: Reconstructions with 5% multiplicative Gaussian noise. (a) Phantom (b) Reconstructed (c) Reconstruction with 5% multiplicative Gaussian noise. Fig. 8: Reconstructions of a set of phantoms. In problems of invert that integrate a function standard results of microl of the object s (true) sin and which parts will be expect to recover stably be tangentially touched in the Radon data (e.g. However, if the generaliz image function along nonto do better than that. Fo produces images of excel basically two (angular) d the fact that the integratio corner at every point of reconstructions do not be The phantom edges do bl broken rays touch them away from the origin in mathematical study of th lack of it is not an easy ta by the authors and their c The artifacts of second of the (incomplete) Radon CT the integrals of the im lines with limited angular of full range [0, 2π]). Th streak artifacts along line equal to the endpoints of CT example above those
20 VRT in Slab Geometry (G.A., R. Gouia-Zarrad 2013) Theorem f C (R 2 ) in D = {(x, y) R 2 0 x x max, 0 y y max }. For (x v, y v ) R 2 and fixed β (0, π 2 ) consider the VRT g(x v, y v ). Then f (x, y) = cos β 2 ( ymax y g(x, y) + tan2 (β) y 2 ) g(x, t) dt. x 2
21 Numerical Implementation (2D)
22 Conical Surfaces of Various Flavors in 3D
23 CRT in Slab Geometry (G.A., R. Gouia-Zarrad 2013) Consider a function f C (R 3 ) supported in D = {(x, y, z) R 3 0 x x max, 0 y y max, 0 z z max }. For (x v, y v, z v ) R 3 we define the 3D conical Radon transform by g(x v, y v, z v ) = f (x, y, z) ds. C(x v,y v,z v )
24 3D Slab Geometry (G.A., R. Gouia-Zarrad 2013) Theorem An exact solution of the inversion problem for CRT is given by z [ ] d 2 2 x f λ,µ (z) = C(β) J 0 (u(z x)) z max dx 2 + u2 ĝ λ,µ (z v ) dz v dx z max where ĝ λ,µ (z v ) and f λ,µ (z) are the 2D Fourier transforms of the functions g(x v, y v, z v ) and f (x, y, z) with respect to the first two variables, C(β) = cos 2 β/(2π sin β) and u = tan β λ 2 + µ 2.
25 Partial Order in R n and Positive Cones A Partially Ordered Vector Space V is a vector space over R together with a partial order such that 1 if x y then x + z y + z for all z V 2 if x 0 then cx 0 for all c R + From the definition we have x y 0 y x and hence the order is completely determined by V + = {x V ; x 0} positive cone of V. Furthermore, for P V there is a partial order on V such that P = V + if and only if P ( P) = {0} P + P P c 0 cp P
26 Partial Order in R n and Positive Cones A Partially Ordered Vector Space V is a vector space over R together with a partial order such that 1 if x y then x + z y + z for all z V 2 if x 0 then cx 0 for all c R + From the definition we have x y 0 y x and hence the order is completely determined by V + = {x V ; x 0} positive cone of V. Furthermore, for P V there is a partial order on V such that P = V + if and only if P ( P) = {0} P + P P c 0 cp P
27 Partial Order in R n and Positive Cones A Partially Ordered Vector Space V is a vector space over R together with a partial order such that 1 if x y then x + z y + z for all z V 2 if x 0 then cx 0 for all c R + From the definition we have x y 0 y x and hence the order is completely determined by V + = {x V ; x 0} positive cone of V. Furthermore, for P V there is a partial order on V such that P = V + if and only if P ( P) = {0} P + P P c 0 cp P
28 Partial Order in R n and Negative Cones We consider partial orders in R n corresponding to negative cones (R n ) B generated by a set of fixed basis vectors B = {v 1,..., v n }, i.e. (R n ) B = { n i=1 c iv i ; c i 0}. In R 2 we will use linearly independent vectors u, v as a generating set for the negative cone. In this case the boundary of the negative cone is a broken line. For f L 1 (R n ) we define F on R n as F (x) = f (y)dµ y x where µ is the Lebesgue measure on R n and y x represents the negative cone at x with respect to partial order on R n.
29 Partial Order in R n and Negative Cones We consider partial orders in R n corresponding to negative cones (R n ) B generated by a set of fixed basis vectors B = {v 1,..., v n }, i.e. (R n ) B = { n i=1 c iv i ; c i 0}. In R 2 we will use linearly independent vectors u, v as a generating set for the negative cone. In this case the boundary of the negative cone is a broken line. For f L 1 (R n ) we define F on R n as F (x) = f (y)dµ y x where µ is the Lebesgue measure on R n and y x represents the negative cone at x with respect to partial order on R n.
30 Partial Order in R n and Negative Cones We consider partial orders in R n corresponding to negative cones (R n ) B generated by a set of fixed basis vectors B = {v 1,..., v n }, i.e. (R n ) B = { n i=1 c iv i ; c i 0}. In R 2 we will use linearly independent vectors u, v as a generating set for the negative cone. In this case the boundary of the negative cone is a broken line. For f L 1 (R n ) we define F on R n as F (x) = f (y)dµ y x where µ is the Lebesgue measure on R n and y x represents the negative cone at x with respect to partial order on R n.
31 Partial Order in R n and Negative Cones F (x) = y x f (y)dµ
32 Cone Differentiation Th. (G.A., M. Latifi-Jebelli 2016) If is the natural order on R, for an integrable function f and F (x) = f (y)dt y x we have F = f almost everywhere. Note that in this case F is absolutely continuous. Can we have a similar result in higher dimensions? We start from two dimensions. Let f be an integrable function on R 2 (with f < ) with respect to Lebesgue measure and define F (x) = y x f (y)dµ using the partial order made by u, v.
33 Cone Differentiation Th. (G.A., M. Latifi-Jebelli 2016) If is the natural order on R, for an integrable function f and F (x) = f (y)dt y x we have F = f almost everywhere. Note that in this case F is absolutely continuous. Can we have a similar result in higher dimensions? We start from two dimensions. Let f be an integrable function on R 2 (with f < ) with respect to Lebesgue measure and define F (x) = y x f (y)dµ using the partial order made by u, v.
34 Cone Differentiation Th. (G.A., M. Latifi-Jebelli 2016) If is the natural order on R, for an integrable function f and F (x) = f (y)dt y x we have F = f almost everywhere. Note that in this case F is absolutely continuous. Can we have a similar result in higher dimensions? We start from two dimensions. Let f be an integrable function on R 2 (with f < ) with respect to Lebesgue measure and define F (x) = y x f (y)dµ using the partial order made by u, v.
35 Cone Differentiation Th. (G.A., M. Latifi-Jebelli 2016) Define V t,s (x) as the average of f over the parallelogram centered at x, sides of length t, s and directions u, v. Note that the area of the parallelogram made with vectors tu, sv is equal to det (tu, sv) = ts det(u, v). Then V t,s (x) = 1 ts det(u, v) [F (x + t 2 u + s 2 v) F (x t 2 u + s 2 v) F (x + t 2 u s 2 v) + F (x t 2 u s 2 v)]
36 Cone Differentiation Th. (G.A., M. Latifi-Jebelli 2016) Define V t,s (x) as the average of f over the parallelogram centered at x, sides of length t, s and directions u, v. Note that the area of the parallelogram made with vectors tu, sv is equal to det (tu, sv) = ts det(u, v). Then V t,s (x) = 1 ts det(u, v) [F (x + t 2 u + s 2 v) F (x t 2 u + s 2 v) F (x + t 2 u s 2 v) + F (x t 2 u s 2 v)]
37 Cone Differentiation Th. (G.A., M. Latifi-Jebelli 2016) Define V t,s (x) as the average of f over the parallelogram centered at x, sides of length t, s and directions u, v. Note that the area of the parallelogram made with vectors tu, sv is equal to det (tu, sv) = ts det(u, v). Then V t,s (x) = 1 ts det(u, v) [F (x + t 2 u + s 2 v) F (x t 2 u + s 2 v) F (x + t 2 u s 2 v) + F (x t 2 u s 2 v)]
38 Cone Differentiation Th. (G.A., M. Latifi-Jebelli 2016) Likewise, for n dimensions by a geometric argument and induction over n we get the following averaging formula for f V t1,...,t n (x) = 1 t 1... t n det(v 1,..., v n ) α { 1 2, 1 2 }n sgn(α 1... α n )F (x + α 1 t 1 v α n t n v n ) In special case, to get a symmetric neighborhood of x we can let t 1 = = t n = t to get the average of f over P t, the parallelograms with sides of length t centered at x V t,...,t (x) = 1 t n det(v 1,..., v n ) P t f dµ = 1 µ(p t ) P t f dµ.
39 Cone Differentiation Th. (G.A., M. Latifi-Jebelli 2016) Likewise, for n dimensions by a geometric argument and induction over n we get the following averaging formula for f V t1,...,t n (x) = 1 t 1... t n det(v 1,..., v n ) α { 1 2, 1 2 }n sgn(α 1... α n )F (x + α 1 t 1 v α n t n v n ) In special case, to get a symmetric neighborhood of x we can let t 1 = = t n = t to get the average of f over P t, the parallelograms with sides of length t centered at x V t,...,t (x) = 1 t n det(v 1,..., v n ) P t f dµ = 1 µ(p t ) P t f dµ.
40 Cone Differentiation Th. (G.A., M. Latifi-Jebelli 2016) Now by averaging over this infinitesimal symmetric neighborhood of x and applying the Lebesgue Differentiation Theorem we have Theorem Let be an order in R n made from the positive cone of v 1,..., v n and for f L 1 (R n ) define F (x) = f (y)dµ. Then for almost every x we have y x f (x) = lim t 0 V t,...,t (x).
41 Cone Differentiation Th. (G.A., M. Latifi-Jebelli 2016) Theorem Let the hypothesis of the previous theorem be satisfied and f be continuous. Then 1 f (x) =... F (x), det(v 1,..., v n ) v 1 v n where v j denotes the directional derivative along vector v j. How does this help us with the inversion of the broken-ray or conical Radon transforms?
42 Cone Differentiation Th. (G.A., M. Latifi-Jebelli 2016) Theorem Let the hypothesis of the previous theorem be satisfied and f be continuous. Then 1 f (x) =... F (x), det(v 1,..., v n ) v 1 v n where v j denotes the directional derivative along vector v j. How does this help us with the inversion of the broken-ray or conical Radon transforms?
43 Inversion of BRT in 2D (G.A., M. Latifi-Jebelli 2016) Assume that L(x, y) is the unique broken ray with vertex at (x, y) and axis of symmetry α, where α = (α x, α y ) is a unit vector parallel to u+v 2. Also, let β be the angle between u and α. Theorem Let R be the broken ray transform on L 1 (R 2 ) defined by: (Rf )(x, y) = fdl. Then F (x, y) = 0 L(x,y) (Rf )(x + tα x, y + tα y ) sin β dt is the integral of f over the negative cone at (x, y). Hence f (x, y) = 1 det(u, v) u v 0 (Rf )(x + tα x, y + tα y ) sin β dt
44 Inversion of BRT in 2D (G.A., M. Latifi-Jebelli 2016) Assume that L(x, y) is the unique broken ray with vertex at (x, y) and axis of symmetry α, where α = (α x, α y ) is a unit vector parallel to u+v 2. Also, let β be the angle between u and α. Theorem Let R be the broken ray transform on L 1 (R 2 ) defined by: (Rf )(x, y) = fdl. Then F (x, y) = 0 L(x,y) (Rf )(x + tα x, y + tα y ) sin β dt is the integral of f over the negative cone at (x, y). Hence f (x, y) = 1 det(u, v) u v 0 (Rf )(x + tα x, y + tα y ) sin β dt
45 Inversion of BRT in 2D (G.A., M. Latifi-Jebelli 2016)
46 Polyhedral CRT in R n (G.A., M. Latifi-Jebelli 2016) For any x R n, we define (Rf )(x) to be integral over the boundary of polyhedral cone C generated by unit basis vectors u 1,..., u n starting from x, i.e. (Rf )(x) = f ds, where ds is n 1 dimensional Lebesgue measure on C. Assume that u i u j is constant for any i and j. Define w = u 1+ +u n u 1 + +u n. Let X i = span u 1,..., u i 1, u i+1,..., u n be the hyperplane containing a face of polyhedral cone and define y i to be a unit vector in X i. C
47 Polyhedral CRT in R n (G.A., M. Latifi-Jebelli 2016) For any x R n, we define (Rf )(x) to be integral over the boundary of polyhedral cone C generated by unit basis vectors u 1,..., u n starting from x, i.e. (Rf )(x) = f ds, where ds is n 1 dimensional Lebesgue measure on C. Assume that u i u j is constant for any i and j. Define w = u 1+ +u n u 1 + +u n. Let X i = span u 1,..., u i 1, u i+1,..., u n be the hyperplane containing a face of polyhedral cone and define y i to be a unit vector in X i. C
48 Polyhedral CRT in R n (G.A., M. Latifi-Jebelli 2016) For any x R n, we define (Rf )(x) to be integral over the boundary of polyhedral cone C generated by unit basis vectors u 1,..., u n starting from x, i.e. (Rf )(x) = f ds, where ds is n 1 dimensional Lebesgue measure on C. Assume that u i u j is constant for any i and j. Define w = u 1+ +u n u 1 + +u n. Let X i = span u 1,..., u i 1, u i+1,..., u n be the hyperplane containing a face of polyhedral cone and define y i to be a unit vector in X i. C
49 Polyhedral CRT in R n (G.A., M. Latifi-Jebelli 2016) For any x R n, we define (Rf )(x) to be integral over the boundary of polyhedral cone C generated by unit basis vectors u 1,..., u n starting from x, i.e. (Rf )(x) = f ds, where ds is n 1 dimensional Lebesgue measure on C. Assume that u i u j is constant for any i and j. Define w = u 1+ +u n u 1 + +u n. Let X i = span u 1,..., u i 1, u i+1,..., u n be the hyperplane containing a face of polyhedral cone and define y i to be a unit vector in X i. C
50 Polyhedral CRT in R n (G.A., M. Latifi-Jebelli 2016) Theorem Let R, w, y j be defined as above, then F (x) = 0 (Rf )(x + wt) w, y 1 dt is the integral of f over the polyhedral cone generated by u 1,... u n starting from x. Hence f (x) = 1... det(v 1,..., v n ) v 1 v n 0 (Rf )(x + wt) w, y 1 dt.
51 Range Description (G.A., M. Latifi-Jebelli 2016) Existence of f such that F (x) = y x f (y) dµ. What is the necessary and sufficient condition for a function F to be a cone integral of another function f 0 with respect to a given order structure in R n? In case of n = 1 the answer was provided by absolute continuity. We apply the Radon Nikodym Theorem to get the desired description of F. For a given F, we construct a corresponding measure ν that implies existence of f. We use the above conditions to obtain a range description for CRT.
52 Range Description (G.A., M. Latifi-Jebelli 2016) Existence of f such that F (x) = y x f (y) dµ. What is the necessary and sufficient condition for a function F to be a cone integral of another function f 0 with respect to a given order structure in R n? In case of n = 1 the answer was provided by absolute continuity. We apply the Radon Nikodym Theorem to get the desired description of F. For a given F, we construct a corresponding measure ν that implies existence of f. We use the above conditions to obtain a range description for CRT.
53 Range Description (G.A., M. Latifi-Jebelli 2016) Existence of f such that F (x) = y x f (y) dµ. What is the necessary and sufficient condition for a function F to be a cone integral of another function f 0 with respect to a given order structure in R n? In case of n = 1 the answer was provided by absolute continuity. We apply the Radon Nikodym Theorem to get the desired description of F. For a given F, we construct a corresponding measure ν that implies existence of f. We use the above conditions to obtain a range description for CRT.
54 Range Description (G.A., M. Latifi-Jebelli 2016) Existence of f such that F (x) = y x f (y) dµ. What is the necessary and sufficient condition for a function F to be a cone integral of another function f 0 with respect to a given order structure in R n? In case of n = 1 the answer was provided by absolute continuity. We apply the Radon Nikodym Theorem to get the desired description of F. For a given F, we construct a corresponding measure ν that implies existence of f. We use the above conditions to obtain a range description for CRT.
55 Range Description (G.A., M. Latifi-Jebelli 2016) Existence of f such that F (x) = y x f (y) dµ. What is the necessary and sufficient condition for a function F to be a cone integral of another function f 0 with respect to a given order structure in R n? In case of n = 1 the answer was provided by absolute continuity. We apply the Radon Nikodym Theorem to get the desired description of F. For a given F, we construct a corresponding measure ν that implies existence of f. We use the above conditions to obtain a range description for CRT.
56 Thanks for Your Attention!
Reconstructing a Function from its V-line Radon Transform in a D
Reconstructing a Function from its V-line Radon Transform in a Disc University of Texas at Arlington Inverse Problems and Application Linköping, Sweden April 5, 2013 Acknowledgements The talk is based
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 information3D Image Reconstruction from Compton camera data
3D Image Reconstruction from Compton camera data Peter Kuchment and Fatma Terzioglu Abstract In this paper, we address analytically and numerically the inversion of the integral transform (cone or Compton
More informationON A CLASS OF GENERALIZED RADON TRANSFORMS AND ITS APPLICATION IN IMAGING SCIENCE
ON A CLASS OF GENERALIZED RADON TRANSFORMS AND ITS APPLICATION IN IMAGING SCIENCE T.T. TRUONG 1 AND M.K. NGUYEN 2 1 University of Cergy-Pontoise, LPTM CNRS UMR 889, F-9532, France e-mail: truong@u-cergy.fr
More informationarxiv: v2 [math.na] 22 Dec 2016
INVERSION OF WEIGHTED DIVERGENT BEAM AND CONE TRANSFORMS arxiv:1612.6772v2 [math.na] 22 Dec 216 Peter Kuchment Department of Mathematics Texas A&M University College Station, TX 77843-3368, USA Fatma Terzioglu
More informationAppliedMathematics. Inversion of the attenuated V-line transform for SPECT with Compton cameras. Markus Haltmeier, Sunghwan Moon, Daniela Schiefeneder
Nr. 27 13. September 216 Leopold-Franzens-Universität Innsbruck Preprint-Series: Department of Mathematics - Applied Mathematics Inversion of the attenuated V-line transform for SPECT with Compton cameras
More informationInversion of a class of circular and elliptical Radon transforms
Inversion of a class of circular and elliptical Radon transforms Gaik Ambartsoumian and Venkateswaran P. Krishnan Abstract. The paper considers a class of elliptical and circular Radon transforms appearing
More informationBroken-Ray and Star Transforms and their Inversion
Broen-Ray and Star Transforms and their Inversion Vadim A Marel, Institut Fresnel Joint wor with Fan Zhao (f. @ UPenn) John Schotland (f. @ Upenn, now @UMich) Lucia Florescu (f. @ UPenn) Interest in using
More informationInversion of the Broken Ray Transform
University of Central Florida Electronic Theses and Dissertations Doctoral Dissertation (Open Access) Inversion of the Broken Ray Transform 2014 Roman Krylov University of Central Florida Find similar
More informationInversion of the broken ray transform in the case of energy-dependent attenuation
Physics in Medicine & Biology PAPER Inversion of the broken ray transform in the case of energy-dependent attenuation To cite this article: R Krylov and A Katsevich 0 Phys. Med. Biol. Manuscript version:
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 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 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 informationarxiv: v1 [math.dg] 24 Feb 2017
GEODESIC X-RAY TOMOGRAPHY FOR PIECEWISE CONSTANT FUNCTIONS ON NONTRAPPING MANIFOLDS JOONAS ILMAVIRTA, JERE LEHTONEN, AND MIKKO SALO arxiv:1702.07622v1 [math.dg] 24 Feb 2017 Abstract. We show that on a
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
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 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 informationcarries the circle w 1 onto the circle z R and sends w = 0 to z = a. The function u(s(w)) is harmonic in the unit circle w 1 and we obtain
4. Poisson formula In fact we can write down a formula for the values of u in the interior using only the values on the boundary, in the case when E is a closed disk. First note that (3.5) determines the
More informationDisturbance Decoupling Problem
DISC Systems and Control Theory of Nonlinear Systems, 21 1 Disturbance Decoupling Problem Lecture 4 Nonlinear Dynamical Control Systems, Chapter 7 The disturbance decoupling problem is a typical example
More informationQuadrature Formula for Computed Tomography
Quadrature Formula for Computed Tomography orislav ojanov, Guergana Petrova August 13, 009 Abstract We give a bivariate analog of the Micchelli-Rivlin quadrature for computing the integral of a function
More informationAn efficient numerical algorithm for the inversion of an integral transform arising in ultrasound imaging
Noname manuscript No. (will be inserted by the editor) An efficient numerical algorithm for the inversion of an integral transform arising in ultrasound imaging Souvik Roy Venkateswaran P. Krishnan Praveen
More informationCompton Camera. Compton Camera
Diagnostic Imaging II Student Project Compton Camera Ting-Tung Chang Introduction The Compton camera operates by exploiting the Compton Effect. It uses the kinematics of Compton scattering to contract
More information(b) If f L p (R), with 1 < p, then Mf L p (R) and. Mf L p (R) C(p) f L p (R) with C(p) depending only on p.
Lecture 3: Carleson Measures via Harmonic Analysis Much of the argument from this section is taken from the book by Garnett, []. The interested reader can also see variants of this argument in the book
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 RADON NIKODYM PROPERTY AND LIPSCHITZ EMBEDDINGS
THE RADON NIKODYM PROPERTY AND LIPSCHITZ EMBEDDINGS Motivation The idea here is simple. Suppose we have a Lipschitz homeomorphism f : X Y where X and Y are Banach spaces, namely c 1 x y f (x) f (y) c 2
More informationx 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7
Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)
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 information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationChange of Variables, Parametrizations, Surface Integrals
Chapter 8 Change of Variables, Parametrizations, Surface Integrals 8. he transformation formula In evaluating any integral, if the integral depends on an auxiliary function of the variables involved, it
More informationarxiv: v1 [math.dg] 7 Jan 2019
RECONSTRUCTION OF PIECEWISE CONSTANT FUNCTIONS FROM X-RAY DATA arxiv:1901.01909v1 [math.dg] 7 Jan 2019 VADIM LEBOVICI Abstract. We show that on a two-dimensional compact nontrapping Riemannian manifold
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationPart II Probability and Measure
Part II Probability and Measure Theorems Based on lectures by J. Miller Notes taken by Dexter Chua Michaelmas 2016 These notes are not endorsed by the lecturers, and I have modified them (often significantly)
More informationLecture 9 February 2, 2016
MATH 262/CME 372: Applied Fourier Analysis and Winter 26 Elements of Modern Signal Processing Lecture 9 February 2, 26 Prof. Emmanuel Candes Scribe: Carlos A. Sing-Long, Edited by E. Bates Outline Agenda:
More informationRestoration of Missing Data in Limited Angle Tomography Based on Helgason- Ludwig Consistency Conditions
Restoration of Missing Data in Limited Angle Tomography Based on Helgason- Ludwig Consistency Conditions Yixing Huang, Oliver Taubmann, Xiaolin Huang, Guenter Lauritsch, Andreas Maier 26.01. 2017 Pattern
More informationMultiple Choice. Compute the Jacobian, (u, v), of the coordinate transformation x = u2 v 4, y = uv. (a) 2u 2 + 4v 4 (b) xu yv (c) 3u 2 + 7v 6
.(5pts) y = uv. ompute the Jacobian, Multiple hoice (x, y) (u, v), of the coordinate transformation x = u v 4, (a) u + 4v 4 (b) xu yv (c) u + 7v 6 (d) u (e) u v uv 4 Solution. u v 4v u = u + 4v 4..(5pts)
More informationTomography of Highly Scattering Media with the Method of Rotated Reference Frames. Vadim Markel Manabu Machida George Panasyuk John Schotland
Tomography of Highly Scattering Media with the Method of Rotated Reference Frames Vadim Markel Manabu Machida George Panasyuk John Schotland Radiology/Bioengineering UPenn, Philadelphia MOTIVATION (The
More informationPROPERTIES OF SOME INTEGRAL TRANSFORMS ARISING IN TOMOGRAPHY. A Dissertation SUNGHWAN MOON
POPETIES OF SOME INTEGAL TANSFOMS AISING IN TOMOGAPHY A Dissertation by SUNGHWAN MOON Submitted to the Office of Graduate and Professional Studies of Texas A&M University in partial fulfillment of the
More informationSINGULAR CURVES OF AFFINE MAXIMAL MAPS
Fundamental Journal of Mathematics and Mathematical Sciences Vol. 1, Issue 1, 014, Pages 57-68 This paper is available online at http://www.frdint.com/ Published online November 9, 014 SINGULAR CURVES
More informationCIRCULAR arc Radon (CAR) transforms involve the integrals
50 IEEE TRANSACTIONS ON COMPUTATIONAL IMAGING, VOL. 2, NO., DECEMBER 2016 Numerical Inversion of Circular arc Radon Transform Tabish A. Syed, Venkateswaran P. Krishnan, and Jayanthi Sivaswamy Abstract
More information9 Radon-Nikodym theorem and conditioning
Tel Aviv University, 2015 Functions of real variables 93 9 Radon-Nikodym theorem and conditioning 9a Borel-Kolmogorov paradox............. 93 9b Radon-Nikodym theorem.............. 94 9c Conditioning.....................
More informationEfficient Solution Methods for Inverse Problems with Application to Tomography Radon Transform and Friends
Efficient Solution Methods for Inverse Problems with Application to Tomography Radon Transform and Friends Alfred K. Louis Institut für Angewandte Mathematik Universität des Saarlandes 66041 Saarbrücken
More informationAero III/IV Conformal Mapping
Aero III/IV Conformal Mapping View complex function as a mapping Unlike a real function, a complex function w = f(z) cannot be represented by a curve. Instead it is useful to view it as a mapping. Write
More informationSigned Measures. Chapter Basic Properties of Signed Measures. 4.2 Jordan and Hahn Decompositions
Chapter 4 Signed Measures Up until now our measures have always assumed values that were greater than or equal to 0. In this chapter we will extend our definition to allow for both positive negative values.
More informationQuantization of the open string on exact plane waves and non-commutative wave fronts
Quantization of the open string on exact plane waves and non-commutative wave fronts F. Ruiz Ruiz (UCM Madrid) Miami 2007, December 13-18 arxiv:0711.2991 [hep-th], with G. Horcajada Motivation On-going
More informationTomography and Reconstruction
Tomography and Reconstruction Lecture Overview Applications Background/history of tomography Radon Transform Fourier Slice Theorem Filtered Back Projection Algebraic techniques Measurement of Projection
More informationThe Mathematics of Maps Lecture 4. Dennis The The Mathematics of Maps Lecture 4 1/29
The Mathematics of Maps Lecture 4 Dennis The The Mathematics of Maps Lecture 4 1/29 Mercator projection Dennis The The Mathematics of Maps Lecture 4 2/29 The Mercator projection (1569) Dennis The The Mathematics
More informationSINGULARITIES OF LAGRANGIAN MEAN CURVATURE FLOW: MONOTONE CASE
SINGULARITIES OF LAGRANGIAN MEAN CURVATURE FLOW: MONOTONE CASE ANDRÉ NEVES Abstract. We study the formation of singularities for the mean curvature flow of monotone Lagrangians in C n. More precisely,
More informationContents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9
MATH 32B-2 (8W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Multiple Integrals 3 2 Vector Fields 9 3 Line and Surface Integrals 5 4 The Classical Integral Theorems 9 MATH 32B-2 (8W)
More informationMath 350 Solutions for Final Exam Page 1. Problem 1. (10 points) (a) Compute the line integral. F ds C. z dx + y dy + x dz C
Math 35 Solutions for Final Exam Page Problem. ( points) (a) ompute the line integral F ds for the path c(t) = (t 2, t 3, t) with t and the vector field F (x, y, z) = xi + zj + xk. (b) ompute the line
More informationCANONICAL EQUATIONS. Application to the study of the equilibrium of flexible filaments and brachistochrone curves. By A.
Équations canoniques. Application a la recherche de l équilibre des fils flexibles et des courbes brachystochrones, Mem. Acad. Sci de Toulouse (8) 7 (885), 545-570. CANONICAL EQUATIONS Application to the
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationInterpolation via weighted l 1 -minimization
Interpolation via weighted l 1 -minimization Holger Rauhut RWTH Aachen University Lehrstuhl C für Mathematik (Analysis) Matheon Workshop Compressive Sensing and Its Applications TU Berlin, December 11,
More information1. Introduction Boundary estimates for the second derivatives of the solution to the Dirichlet problem for the Monge-Ampere equation
POINTWISE C 2,α ESTIMATES AT THE BOUNDARY FOR THE MONGE-AMPERE EQUATION O. SAVIN Abstract. We prove a localization property of boundary sections for solutions to the Monge-Ampere equation. As a consequence
More informationRadon Transform An Introduction
Radon Transform An Introduction Yi-Hsuan Lin The Radon transform is widely applicable to tomography, the creation of an image from the scattering data associated to crosssectional scans of an object. If
More informationRich Tomography. Bill Lionheart, School of Mathematics, University of Manchester and DTU Compute. July 2014
Rich Tomography Bill Lionheart, School of Mathematics, University of Manchester and DTU Compute July 2014 What do we mean by Rich Tomography? Conventional tomography reconstructs one scalar image from
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 informationExercises for Multivariable Differential Calculus XM521
This document lists all the exercises for XM521. The Type I (True/False) exercises will be given, and should be answered, online immediately following each lecture. The Type III exercises are to be done
More informationMoment Measures. Bo az Klartag. Tel Aviv University. Talk at the asymptotic geometric analysis seminar. Tel Aviv, May 2013
Tel Aviv University Talk at the asymptotic geometric analysis seminar Tel Aviv, May 2013 Joint work with Dario Cordero-Erausquin. A bijection We present a correspondence between convex functions and Borel
More informationVectors in Special Relativity
Chapter 2 Vectors in Special Relativity 2.1 Four - vectors A four - vector is a quantity with four components which changes like spacetime coordinates under a coordinate transformation. We will write the
More informationMTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1.
MTH4101 CALCULUS II REVISION NOTES 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) 1.1 Introduction Types of numbers (natural, integers, rationals, reals) The need to solve quadratic equations:
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 informationAN INVERSE PROBLEM FOR THE WAVE EQUATION WITH A TIME DEPENDENT COEFFICIENT
AN INVERSE PROBLEM FOR THE WAVE EQUATION WITH A TIME DEPENDENT COEFFICIENT Rakesh Department of Mathematics University of Delaware Newark, DE 19716 A.G.Ramm Department of Mathematics Kansas State University
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationFeature Reconstruction in Tomography
Feature Reconstruction in Tomography Alfred K. Louis Institut für Angewandte Mathematik Universität des Saarlandes 66041 Saarbrücken http://www.num.uni-sb.de louis@num.uni-sb.de Wien, July 20, 2009 Images
More informationAverage theorem, Restriction theorem and Strichartz estimates
Average theorem, Restriction theorem and trichartz estimates 2 August 27 Abstract We provide the details of the proof of the average theorem and the restriction theorem. Emphasis has been placed on the
More informationHyperbolic Geometry on Geometric Surfaces
Mathematics Seminar, 15 September 2010 Outline Introduction Hyperbolic geometry Abstract surfaces The hemisphere model as a geometric surface The Poincaré disk model as a geometric surface Conclusion Introduction
More informationREPRESENTATION THEORY WEEK 7
REPRESENTATION THEORY WEEK 7 1. Characters of L k and S n A character of an irreducible representation of L k is a polynomial function constant on every conjugacy class. Since the set of diagonalizable
More informationHigh-Resolution. Transmission. Electron Microscopy
Part 4 High-Resolution Transmission Electron Microscopy 186 Significance high-resolution transmission electron microscopy (HRTEM): resolve object details smaller than 1nm (10 9 m) image the interior of
More information7a3 2. (c) πa 3 (d) πa 3 (e) πa3
1.(6pts) Find the integral x, y, z d S where H is the part of the upper hemisphere of H x 2 + y 2 + z 2 = a 2 above the plane z = a and the normal points up. ( 2 π ) Useful Facts: cos = 1 and ds = ±a sin
More informationSOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)
SOLUTIONS TO THE 18.02 FINAL EXAM BJORN POONEN December 14, 2010, 9:00am-12:00 (3 hours) 1) For each of (a)-(e) below: If the statement is true, write TRUE. If the statement is false, write FALSE. (Please
More informationLecture-05 Perturbation Theory and Feynman Diagrams
Lecture-5 Perturbation Theory and Feynman Diagrams U. Robkob, Physics-MUSC SCPY639/428 September 3, 218 From the previous lecture We end up at an expression of the 2-to-2 particle scattering S-matrix S
More informationMonte Carlo Radiation Transfer I
Monte Carlo Radiation Transfer I Monte Carlo Photons and interactions Sampling from probability distributions Optical depths, isotropic emission, scattering Monte Carlo Basics Emit energy packet, hereafter
More informationNIH Public Access Author Manuscript Inverse Probl. Author manuscript; available in PMC 2010 April 20.
NIH Public Access Author Manuscript Published in final edited form as: Inverse Probl. 2010 January 1; 26(3): 350131 3501329. doi:10.1088/0266-5611/26/3/035013. High Order Total Variation Minimization for
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 informationExact Solutions of the Einstein Equations
Notes from phz 6607, Special and General Relativity University of Florida, Fall 2004, Detweiler Exact Solutions of the Einstein Equations These notes are not a substitute in any manner for class lectures.
More informationMultivariable Calculus Notes. Faraad Armwood. Fall: Chapter 1: Vectors, Dot Product, Cross Product, Planes, Cylindrical & Spherical Coordinates
Multivariable Calculus Notes Faraad Armwood Fall: 2017 Chapter 1: Vectors, Dot Product, Cross Product, Planes, Cylindrical & Spherical Coordinates Chapter 2: Vector-Valued Functions, Tangent Vectors, Arc
More informationA Brief Introduction to Medical Imaging. Outline
A Brief Introduction to Medical Imaging Outline General Goals Linear Imaging Systems An Example, The Pin Hole Camera Radiations and Their Interactions with Matter Coherent vs. Incoherent Imaging Length
More informationMath 265H: Calculus III Practice Midterm II: Fall 2014
Name: Section #: Math 65H: alculus III Practice Midterm II: Fall 14 Instructions: This exam has 7 problems. The number of points awarded for each question is indicated in the problem. Answer each question
More informationHarmonic Functions and Brownian motion
Harmonic Functions and Brownian motion Steven P. Lalley April 25, 211 1 Dynkin s Formula Denote by W t = (W 1 t, W 2 t,..., W d t ) a standard d dimensional Wiener process on (Ω, F, P ), and let F = (F
More informationCompton Storage Rings
Compton Polarimetry @ Storage Rings Wolfgang Hillert ELectron Stretcher Accelerator Physics Institute of Bonn University Møller-Polarimeter Compton-Polarimeter Mott-Polarimeter Compton Scattering Differential
More informationI. Introduction. A New Method for Inverting Integrals Attenuated Radon Transform (SPECT) D to N map for Moving Boundary Value Problems
I. Introduction A New Method for Inverting Integrals Attenuated Radon Transform (SPECT) D to N map for Moving Boundary Value Problems F(k) = T 0 e k2 t+ikl(t) f (t)dt, k C. Integrable Nonlinear PDEs in
More informationLorentz Force. Acceleration of electrons due to the magnetic field gives rise to synchrotron radiation Lorentz force.
Set 10: Synchrotron Lorentz Force Acceleration of electrons due to the magnetic field gives rise to synchrotron radiation Lorentz force 0 E x E y E z dp µ dτ = e c F µ νu ν, F µ E x 0 B z B y ν = E y B
More informationMethod of Rotated Reference Frames
Method of Rotated Reference Frames Vadim Markel in collaboration with George Panasyuk Manabu Machida John Schotland Radiology/Bioengineering UPenn, Philadelphia Outline MOTIVATION (the inverse problems
More informationElectromagnetic. G. A. Krafft Jefferson Lab Jefferson Lab Professor of Physics Old Dominion University Physics 804 Electromagnetic Theory II
Physics 704/804 Electromagnetic Theory II G. A. Krafft Jefferson Lab Jefferson Lab Professor of Physics Old Dominion University 04-13-10 4-Vectors and Proper Time Any set of four quantities that transform
More information1 Vector fields, flows and Lie derivatives
Working title: Notes on Lie derivatives and Killing vector fields Author: T. Harmark Killingvectors.tex 5//2008, 22:55 Vector fields, flows and Lie derivatives Coordinates Consider an n-dimensional manifold
More informationIntroduction to Algebraic and Geometric Topology Week 14
Introduction to Algebraic and Geometric Topology Week 14 Domingo Toledo University of Utah Fall 2016 Computations in coordinates I Recall smooth surface S = {f (x, y, z) =0} R 3, I rf 6= 0 on S, I Chart
More informationSolutions to Tutorial 11 (Week 12)
THE UIVERSITY OF SYDEY SCHOOL OF MATHEMATICS AD STATISTICS Solutions to Tutorial 11 (Week 12) MATH3969: Measure Theory and Fourier Analysis (Advanced) Semester 2, 2017 Web Page: http://sydney.edu.au/science/maths/u/ug/sm/math3969/
More informationROI Reconstruction in CT
Reconstruction in CT From Tomo reconstruction in the 21st centery, IEEE Sig. Proc. Magazine (R.Clackdoyle M.Defrise) L. Desbat TIMC-IMAG September 10, 2013 L. Desbat Reconstruction in CT Outline 1 CT Radiology
More informationBessel Functions Michael Taylor. Lecture Notes for Math 524
Bessel Functions Michael Taylor Lecture Notes for Math 54 Contents 1. Introduction. Conversion to first order systems 3. The Bessel functions J ν 4. The Bessel functions Y ν 5. Relations between J ν and
More informationThe Jacobian. v y u y v (see exercise 46). The total derivative is also known as the Jacobian Matrix of the transformation T (u; v) :
The Jacobian The Jacobian of a Transformation In this section, we explore the concept of a "derivative" of a coordinate transformation, which is known as the Jacobian of the transformation. However, in
More informationCalculation of Momentum Distribution Function of a Non-Thermal Fermionic Dark Matter
Calculation of Momentum Distribution Function of a Non-Thermal Fermionic Dark Matter, March 8, 2017. arxiv:1612.02793, with Anirban Biswas. Aritra Gupta Why Non-Thermal? 1 / 31 The most widely studied
More informationLuminy Lecture 2: Spectral rigidity of the ellipse joint work with Hamid Hezari. Luminy Lecture April 12, 2015
Luminy Lecture 2: Spectral rigidity of the ellipse joint work with Hamid Hezari Luminy Lecture April 12, 2015 Isospectral rigidity of the ellipse The purpose of this lecture is to prove that ellipses are
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 informationRudin Real and Complex Analysis - Harmonic Functions
Rudin Real and Complex Analysis - Harmonic Functions Aaron Lou December 2018 1 Notes 1.1 The Cauchy-Riemann Equations 11.1: The Operators and Suppose f is a complex function defined in a plane open set
More informationElectromagnetic Waves Across Interfaces
Lecture 1: Foundations of Optics Outline 1 Electromagnetic Waves 2 Material Properties 3 Electromagnetic Waves Across Interfaces 4 Fresnel Equations 5 Brewster Angle 6 Total Internal Reflection Christoph
More informationGeneral Instructions
Math 240: Real Analysis Qualifying Exam May 28, 2015 Name: Student ID#: Problems/Page Numbers Total Points Your Score Problem 1 / Page 2-3 40 Points Problem 2 / Page 4 Problem 3 / Page 5 Problem 4 / Page
More informationFermionic Projectors and Hadamard States
Fermionic Projectors and Hadamard States Simone Murro Fakultät für Mathematik Universität Regensburg Foundations and Constructive Aspects of QFT Göttingen, 16th of January 2016 To the memory of Rudolf
More informationA Note on the Differential Equations with Distributional Coefficients
MATEMATIKA, 24, Jilid 2, Bil. 2, hlm. 115 124 c Jabatan Matematik, UTM. A Note on the Differential Equations with Distributional Coefficients Adem Kilicman Department of Mathematics, Institute for Mathematical
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 information