LECTURES ON MICROLOCAL CHARACTERIZATIONS IN LIMITED-ANGLE
|
|
- Emma Weaver
- 5 years ago
- Views:
Transcription
1 LECTURES ON MICROLOCAL CHARACTERIZATIONS IN LIMITED-ANGLE TOMOGRAPHY Jürgen Frikel
2 4 LECTURES 1 Nov. 11: Introduction to the mathematics of computerized tomography 2 Today: Introduction to the basic concepts of microlocal analysis 3 Nov. 25: Microlocal analysis of limited angle reconstructions in tomography I 4 Dec. 02: Microlocal analysis of limited angle reconstructions in tomography I 2 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
3 4 LECTURES 1 Nov. 11: Introduction to the mathematics of computerized tomography 2 Today: Introduction to the basic concepts of microlocal analysis 3 Nov. 25: Microlocal analysis of limited angle reconstructions in tomography I 4 Dec. 02: Microlocal analysis of limited angle reconstructions in tomography I References: L. Hörmander, The analysis of linear partial differential operators I: Distribution theory and Fourier analysis, vol Berlin: Springer-Verlag, JF and E. T. Quinto, Characterization and reduction of artifacts in limited angle tomography, Inverse Problems 29(12):125007, December Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
4 A REMINDER ON FBP TYPE INVERSION FORMULAS We consider the Radon transform R in 2D and filtered backprojection inversion formulas. Data: g = R f. Classical reconstruction: f (x) = 1 4π R Is 1 g = 1 4π R ( 2 s) 1/2 g = 1 I 1 4π S 1 s g(θ, x θ) dθ = 1 [ψ s g](θ, x θ) dθ, 4π S 1 where the filter ψ is defined in the Fourier domain via ψ(σ) = σ. 3 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
5 A REMINDER ON FBP TYPE INVERSION FORMULAS We consider the Radon transform R in 2D and filtered backprojection inversion formulas. Data: g = R f. Classical reconstruction: f (x) = 1 4π R Is 1 g = 1 4π R ( 2 s) 1/2 g = 1 I 1 4π S 1 s g(θ, x θ) dθ = 1 [ψ s g](θ, x θ) dθ, 4π S 1 where the filter ψ is defined in the Fourier domain via ψ(σ) = σ. Lambda reconstruction: Λ f (x) = 1 4π R ( 2 s)g = 1 [ψ Λ s g](θ, x θ) dθ 4π S 1 where the filter ψ is defined in the Fourier domain via ψ Λ (σ) = σ 2. 3 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
6 GENERAL FORM FBP TYPE INVERSION FORMULAS In Lambda reconstruction a differential operator ( 2 s) is used as filter, which can be written in terms of the Fourier transform as follows (recall F ( s g) = iσ ĝ): ( 2 s)g(s) = F ( 1 σ 2 ĝ(σ) ) = (2π) 1/2 σ 2 ĝ(σ)e isσ ds. R 4 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
7 GENERAL FORM FBP TYPE INVERSION FORMULAS In Lambda reconstruction a differential operator ( 2 s) is used as filter, which can be written in terms of the Fourier transform as follows (recall F ( s g) = iσ ĝ): ( 2 s)g(s) = F ( 1 σ 2 ĝ(σ) ) = (2π) 1/2 σ 2 ĝ(σ)e isσ ds. R In classical reconstructions a fractional power of a differential operator ( 2 s) 1/2 is used as filter. We have a similar representation in the Fourier domain ( 2 s) 1/2 g(s) = F 1 ( σ ĝ(σ)) = (2π) 1/2 σ ĝ(σ)e isσ ds. R However, the difference is huge as this operator is not local any more (involves Hilbert transforms). This is a so-called pseudodifferential operator. 4 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
8 GENERAL FORM FBP TYPE INVERSION FORMULAS In Lambda reconstruction a differential operator ( 2 s) is used as filter, which can be written in terms of the Fourier transform as follows (recall F ( s g) = iσ ĝ): ( 2 s)g(s) = F ( 1 σ 2 ĝ(σ) ) = (2π) 1/2 σ 2 ĝ(σ)e isσ ds. R In classical reconstructions a fractional power of a differential operator ( 2 s) 1/2 is used as filter. We have a similar representation in the Fourier domain ( 2 s) 1/2 g(s) = F 1 ( σ ĝ(σ)) = (2π) 1/2 σ ĝ(σ)e isσ ds. R However, the difference is huge as this operator is not local any more (involves Hilbert transforms). This is a so-called pseudodifferential operator. Let P be any of the above filters, then reconstruction means applying the operator to data g = R f. R := R P 4 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
9 GENERAL FORM FBP TYPE INVERSION FORMULAS In Lambda reconstruction a differential operator ( 2 s) is used as filter, which can be written in terms of the Fourier transform as follows (recall F ( s g) = iσ ĝ): ( 2 s)g(s) = F ( 1 σ 2 ĝ(σ) ) = (2π) 1/2 σ 2 ĝ(σ)e isσ ds. R In classical reconstructions a fractional power of a differential operator ( 2 s) 1/2 is used as filter. We have a similar representation in the Fourier domain ( 2 s) 1/2 g(s) = F 1 ( σ ĝ(σ)) = (2π) 1/2 σ ĝ(σ)e isσ ds. R However, the difference is huge as this operator is not local any more (involves Hilbert transforms). This is a so-called pseudodifferential operator. Let P be any of the above filters, then reconstruction means applying the operator to data g = R f. R := R P General FBP type reconstruction formulas use general pseudodifferential operators P as filters. 4 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
10 GENERAL FORM FBP TYPE INVERSION FORMULAS Since FBP type reconstructions are so easy to implement and since the yield very fast reconstructions, the are applied to all kinds of tomographic data, in particular to incomplete data 5 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
11 GENERAL FORM FBP TYPE INVERSION FORMULAS Since FBP type reconstructions are so easy to implement and since the yield very fast reconstructions, the are applied to all kinds of tomographic data, in particular to incomplete data Limited angle data Sparse angle data Interior data Exterior data... 5 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
12 GENERAL FORM FBP TYPE INVERSION FORMULAS Since FBP type reconstructions are so easy to implement and since the yield very fast reconstructions, the are applied to all kinds of tomographic data, in particular to incomplete data Limited angle data Sparse angle data Interior data Exterior data... In limited angle tomography, the projectios g θ = R θ f are known only for certain directions θ S 1 Φ S 1, for other directions θ the projections g θ = R θ f are unknown. FBP Reconstruction of limited angle data R ( R Φ f )(x) = 1 4π S 1 Φ [ψ s g](θ, x θ) dθ =??? 5 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
13 GENERAL FORM FBP TYPE INVERSION FORMULAS Since FBP type reconstructions are so easy to implement and since the yield very fast reconstructions, the are applied to all kinds of tomographic data, in particular to incomplete data Limited angle data Sparse angle data Interior data Exterior data... In limited angle tomography, the projectios g θ = R θ f are known only for certain directions θ S 1 Φ S 1, for other directions θ the projections g θ = R θ f are unknown. FBP Reconstruction of limited angle data R ( R Φ f )(x) = 1 4π S 1 Φ [ψ s g](θ, x θ) dθ =??? What do we reconstruct? 5 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
14 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 1 Original FBP Lambda [0, 20 ] 6 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
15 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 1 Original FBP Lambda [0, 40 ] 6 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
16 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 1 Original FBP Lambda [0, 60 ] 6 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
17 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 1 Original FBP Lambda [0, 80 ] 6 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
18 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 1 Original FBP Lambda [0, 100 ] 6 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
19 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 1 Original FBP Lambda [0, 120 ] 6 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
20 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 1 Original FBP Lambda [0, 140 ] 6 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
21 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 1 Original FBP Lambda [0, 160 ] 6 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
22 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 1 Original FBP Lambda [0, 180 ] 6 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
23 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 2 Original FBP Lambda [0, 90 ] 7 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
24 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 2 Original FBP Lambda [20, 110 ] 7 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
25 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 2 Original FBP Lambda [40, 130 ] 7 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
26 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 2 Original FBP Lambda [60, 150 ] 7 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
27 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 2 Original FBP Lambda [80, 170 ] 7 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
28 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 2 Original FBP Lambda [100, 190 ] 7 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
29 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 2 Original FBP Lambda [120, 210 ] 7 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
30 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 2 Original FBP Lambda [140, 230 ] 7 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
31 LIMITED ANGLE RECONSTRUCTIONS: EXPERIMENT 2 Original FBP Lambda [160, 250 ] 7 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
32 CHALLENGES IN LIMITED VIEW TOMOGRAPHY Observations at a first glance: Only certain features of the original object can be reconstructed Artifacts are generated 8 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
33 CHALLENGES IN LIMITED VIEW TOMOGRAPHY Observations at a first glance: Only certain features of the original object can be reconstructed Artifacts are generated Goal: Characterize image features that can be reliably reconstructed? Develop strategy to reduce artifacts? 8 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
34 CHALLENGES IN LIMITED VIEW TOMOGRAPHY Observations at a first glance: Only certain features of the original object can be reconstructed Artifacts are generated Goal: Characterize image features that can be reliably reconstructed? Develop strategy to reduce artifacts? Observations at a second glance: Common characteristic features for both types of reconstructions are edges Information about edges is given in terms of projection directions 8 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
35 WHAT ARE EDGES? Practically: Density jumps, boundaries between regions, etc. f 9 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
36 WHAT ARE EDGES? Practically: Density jumps, boundaries between regions, etc. f Singularities of f 9 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
37 WHAT ARE EDGES? Practically: Density jumps, boundaries between regions, etc. Mathematically: Where the function is not smooth singularities. f Singularities of f 9 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
38 WHAT ARE EDGES? Practically: Density jumps, boundaries between regions, etc. Mathematically: Where the function is not smooth singularities. Edges are local and oriented: What s an appropriate mathematical framework? f Singularities of f 9 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
39 WHAT ARE EDGES? Practically: Density jumps, boundaries between regions, etc. Mathematically: Where the function is not smooth singularities. Edges are local and oriented: What s an appropriate mathematical framework? Image processing: Location of singularities = locations of high gradient magnitude; Direction of singularities = gradient directions f Singularities of f 9 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
40 WHAT ARE EDGES? Practically: Density jumps, boundaries between regions, etc. Mathematically: Where the function is not smooth singularities. Edges are local and oriented: What s an appropriate mathematical framework? More powerful framework: Microlocal Analysis f Singularities of f 9 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
41 TODAY Global smoothness vs. decay of Fourier transform Singular locations Singular directions Wavefront sets: simultaneous description of singular locations and directions Examples 10 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
42 GLOBAL SMOOTHNESS VS. DECAY OF FOURIER TRANSFORM In what follows, we only consider real functions f : R 2 R. The smoothness of any function f can be charachterized by means of the decay of its Fourier transform f ˆ: f C k f ˆ(ξ) = O( ξ k ) for ξ. In particular, f C ˆ f (ξ) = O( ξ k ) for all k N. 11 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
43 GLOBAL SMOOTHNESS VS. DECAY OF FOURIER TRANSFORM In what follows, we only consider real functions f : R 2 R. The smoothness of any function f can be charachterized by means of the decay of its Fourier transform f ˆ: f C k f ˆ(ξ) = O( ξ k ) for ξ. In particular, f C ˆ f (ξ) = O( ξ k ) for all k N. Notation: A function which satisfies f (x) = O( x k ) for all k N is called rapidly decreasing. In all other cases slowly decreasing. Definition A function f which is not C is called singular. Equivalently, f is singular if the Fourier transform does not decay rapidly. 11 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
44 GLOBAL SMOOTHNESS VS. DECAY OF FOURIER TRANSFORM In what follows, we only consider real functions f : R 2 R. The smoothness of any function f can be charachterized by means of the decay of its Fourier transform f ˆ: f C k f ˆ(ξ) = O( ξ k ) for ξ. In particular, f C ˆ f (ξ) = O( ξ k ) for all k N. Notation: A function which satisfies f (x) = O( x k ) for all k N is called rapidly decreasing. In all other cases slowly decreasing. Definition A function f which is not C is called singular. Equivalently, f is singular if the Fourier transform does not decay rapidly. Above relations can be proven by using the property F (D α f ) = i α ξ α F f Above relations also hold for distributions in D and S If the Fourier transform of f has compact support (fastest decay one can imagine), then it can be shown that the function f is analytic (Paley-Wiener) One can also use the above relations to define smoothness of order α R (Sobolev-Spaces) 11 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
45 GLOBAL SMOOTHNESS VS. DECAY OF FOURIER TRANSFORM Limitation: No information about locations or directions is provided. Global smoothness is related to a non-directional decay of the Fourier transform. 12 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
46 GLOBAL SMOOTHNESS VS. DECAY OF FOURIER TRANSFORM Limitation: No information about locations or directions is provided. Global smoothness is related to a non-directional decay of the Fourier transform. Example: Let f be a function which is smooth except at the point 0 R 2, i.e. f C (R 2 \ 0). Let y R 2 be arbitrary and let f y (x) = f (x y). Then and hence fˆ y (ξ) = e iξy f ˆ(ξ) ˆ f y (ξ) = ˆ f (ξ). 12 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
47 GLOBAL SMOOTHNESS VS. DECAY OF FOURIER TRANSFORM Limitation: No information about locations or directions is provided. Global smoothness is related to a non-directional decay of the Fourier transform. Example: Let f be a function which is smooth except at the point 0 R 2, i.e. f C (R 2 \ 0). Let y R 2 be arbitrary and let f y (x) = f (x y). Then and hence fˆ y (ξ) = e iξy f ˆ(ξ) ˆ f y (ξ) = ˆ f (ξ). That is, we can generate functions f and f y that have equal decay properties but arbitrarily different locations of jums, discontinuities, etc.. Hence, the decay of the Fourier transform does not provide information about the location of discontinuities. 12 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
48 SINGULAR LOCATIONS How can we gain knowledge about locations of jump, discontinuieties, etc.? What is a proper concept of a singularity that can be used in a distributional setting? 13 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
49 SINGULAR LOCATIONS How can we gain knowledge about locations of jump, discontinuieties, etc.? What is a proper concept of a singularity that can be used in a distributional setting? Idea: Localize and study decay of the Fourier transforms! Definition Let f be a function or distribution. We say that x 0 R 2 is a regular point of f, if there exists a (cut-off) function ϕ C c (R 2 ) with ϕ(x 0 ) 0 such that the Fourier transform F (ϕ f ) decays rapidly (or in other words s.t. ϕ f C (R 2 )). 13 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
50 SINGULAR LOCATIONS How can we gain knowledge about locations of jump, discontinuieties, etc.? What is a proper concept of a singularity that can be used in a distributional setting? Idea: Localize and study decay of the Fourier transforms! Definition Let f be a function or distribution. We say that x 0 R 2 is a regular point of f, if there exists a (cut-off) function ϕ C c (R 2 ) with ϕ(x 0 ) 0 such that the Fourier transform F (ϕ f ) decays rapidly (or in other words s.t. ϕ f C (R 2 )). The complement of the (open) set of all regular points is called the singular support of f and is denoted by sing supp( f ). It can be explicitly stated as sing supp( f ) = { x R 2 : ϕ C c (R 2 ), ϕ(x) 0 : ϕ f C (R 2 ) }. sing supp( f ) supp( f ) sing supp( f ) = iff f C (R 2 ) 13 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
51 SINGULAR DIRECTIONS So far: Locations of singularities are given by the singular support and x sing supp f implies that F (ϕ f ) decays slowly for any cut-off function ϕ(x 0 ) Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
52 SINGULAR DIRECTIONS So far: Locations of singularities are given by the singular support and x sing supp f implies that F (ϕ f ) decays slowly for any cut-off function ϕ(x 0 ) 0. Can there exist directions along which F (ϕ f ) has rapid decay? 14 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
53 SINGULAR DIRECTIONS So far: Locations of singularities are given by the singular support and x sing supp f implies that F (ϕ f ) decays slowly for any cut-off function ϕ(x 0 ) 0. Can there exist directions along which F (ϕ f ) has rapid decay? (Example with a line singularity suggests that this might be the case) 14 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
54 SINGULAR DIRECTIONS So far: Locations of singularities are given by the singular support and x sing supp f implies that F (ϕ f ) decays slowly for any cut-off function ϕ(x 0 ) 0. Can there exist directions along which F (ϕ f ) has rapid decay? (Example with a line singularity suggests that this might be the case) Definition Let f be a be compactly supported distribution. A direction ξ 0 R 2 \ 0 is called regular direction if there exists a conic neighborhood N of ξ 0 such that fˆ decays fast in N. 14 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
55 SINGULAR DIRECTIONS So far: Locations of singularities are given by the singular support and x sing supp f implies that F (ϕ f ) decays slowly for any cut-off function ϕ(x 0 ) 0. Can there exist directions along which F (ϕ f ) has rapid decay? (Example with a line singularity suggests that this might be the case) Definition Let f be a be compactly supported distribution. A direction ξ 0 R 2 \ 0 is called regular direction if there exists a conic neighborhood N of ξ 0 such that fˆ decays fast in N. The complement of the (open) set of all regular directions is called the frequency set of f and is denoted by Σ( f ). It can be explicitly stated as Σ( f ) = { ξ R 2 \ 0 : conic neighb. N of ξ : ˆ f decays slowly in N }. Σ( f ) = iff f C (R 2 ) Singular directions only exist if singularities exist and vice versa. For any compactly supported distribution f and any ϕ Cc (R 2 ) we have Σ(ϕ f ) Σ( f ), i.e., multiplication by compactly supported smooth function does not add singular directions. 14 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
56 LOCATION AND DIRECTION OF A SINGULARITY How can we combine the concepts of location and direction of a singularity? 15 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
57 LOCATION AND DIRECTION OF A SINGULARITY How can we combine the concepts of location and direction of a singularity? The property Σ(ϕ f ) Σ( f ) is the key. 15 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
58 LOCATION AND DIRECTION OF A SINGULARITY How can we combine the concepts of location and direction of a singularity? The property Σ(ϕ f ) Σ( f ) is the key. Idea: ZOOM IN 15 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
59 LOCATION AND DIRECTION OF A SINGULARITY How can we combine the concepts of location and direction of a singularity? The property Σ(ϕ f ) Σ( f ) is the key. Idea: ZOOM IN Choose ϕ k C c (R 2 ) such that supp(ϕ k ) {x} as k, then the idea is that the limit lim k Σ(ϕ f ) will contain only singular directions that correspond to the singularity in x. 15 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
60 LOCATION AND DIRECTION OF A SINGULARITY How can we combine the concepts of location and direction of a singularity? The property Σ(ϕ f ) Σ( f ) is the key. Idea: ZOOM IN Choose ϕ k C c (R 2 ) such that supp(ϕ k ) {x} as k, then the idea is that the limit lim k Σ(ϕ f ) will contain only singular directions that correspond to the singularity in x. Definition For a distribution f and x 0 R 2 ( f ) we define the localized frequency set as Σ x0 ( f ) = Σ(ϕ f ). ϕ C c,ϕ(x 0 ) 0 15 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
61 LOCATION AND DIRECTION OF A SINGULARITY How can we combine the concepts of location and direction of a singularity? The property Σ(ϕ f ) Σ( f ) is the key. Idea: ZOOM IN Choose ϕ k C c (R 2 ) such that supp(ϕ k ) {x} as k, then the idea is that the limit lim k Σ(ϕ f ) will contain only singular directions that correspond to the singularity in x. Definition For a distribution f and x 0 R 2 ( f ) we define the localized frequency set as Σ x0 ( f ) = Σ(ϕ f ). ϕ C c,ϕ(x 0 ) 0 Σ x ( f ) Σ( f ) Σ x ( f ) iff x sing supp( f ) Σ x ( f ) is the set of singular directions of f at location x 15 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
62 THE WAVEFRONT SET Definition The wavefront set of a distribution is the set of all tuples (x, ξ) where x is a singular location of f and ξ is a singular direction of f at x. That is, WF( f ) = { (x, ξ) R 2 (R 2 \ 0) : ξ Σ x ( f ) }. 16 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
63 THE WAVEFRONT SET Definition The wavefront set of a distribution is the set of all tuples (x, ξ) where x is a singular location of f and ξ is a singular direction of f at x. That is, WF( f ) = { (x, ξ) R 2 (R 2 \ 0) : ξ Σ x ( f ) }. WF simultaneously describes locations and directions of singularities: where π 1 (x, ξ) = x and π 2 (x, ξ) = ξ. π 1 (WF( f )) = sing supp( f ), π 2 (WF( f )) = Σ( f ), 16 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
64 THE WAVEFRONT SET Definition The wavefront set of a distribution is the set of all tuples (x, ξ) where x is a singular location of f and ξ is a singular direction of f at x. That is, WF( f ) = { (x, ξ) R 2 (R 2 \ 0) : ξ Σ x ( f ) }. WF simultaneously describes locations and directions of singularities: where π 1 (x, ξ) = x and π 2 (x, ξ) = ξ. π 1 (WF( f )) = sing supp( f ), π 2 (WF( f )) = Σ( f ), Alternative representation: WF( f ) = {x} Σ x ( f ). x sing supp 16 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
65 THE WAVEFRONT SET Definition The wavefront set of a distribution is the set of all tuples (x, ξ) where x is a singular location of f and ξ is a singular direction of f at x. That is, WF( f ) = { (x, ξ) R 2 (R 2 \ 0) : ξ Σ x ( f ) }. WF simultaneously describes locations and directions of singularities: where π 1 (x, ξ) = x and π 2 (x, ξ) = ξ. π 1 (WF( f )) = sing supp( f ), π 2 (WF( f )) = Σ( f ), Alternative representation: WF( f ) = {x} Σ x ( f ). x sing supp A characterization of WF in terms of the complement: (x 0, ξ 0 ) WF( f ) ϕ C c (R 2 ), ϕ(x 0 ) 0 : ξ 0 Σ(ϕ f ) ϕ C c (R 2 ), ϕ(x 0 ) 0 : conic neighb. N(ξ 0 ) : k N : F (ϕ f ) = O( ξ k ) 16 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
66 EXAMPLES 17 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
67 EXAMPLES Example: Ω R 2 such that the boundary Ω is a smooth manifold: (x, ξ) WF(χ Ω ) x Ω, and ξ N x, where N x is the normal space to Ω at x Ω. 17 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
68 EXAMPLES Example: Ω R 2 such that the boundary Ω is a smooth manifold: (x, ξ) WF(χ Ω ) x Ω, and ξ N x, where N x is the normal space to Ω at x Ω. ξ x f distribution with supp( f ) D (x, ξdx) WF( f ) 17 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
69 EXAMPLES Example: Ω R 2 such that the boundary Ω is a smooth manifold: (x, ξ) WF(χ Ω ) x Ω, and ξ N x, where N x is the normal space to Ω at x Ω. x ξ (x, ξ) WF( f ) f Singularities of f 17 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
70 EXAMPLES Example: Ω R 2 such that the boundary Ω is a smooth manifold: (x, ξ) WF(χ Ω ) x Ω, and ξ N x, where N x is the normal space to Ω at x Ω. What happens at crossings and corners? x ξ (x, ξ) WF( f ) f Singularities of f 17 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
71 EXAMPLES Example: Ω R 2 such that the boundary Ω is a smooth manifold: (x, ξ) WF(χ Ω ) x Ω, and ξ N x, where N x is the normal space to Ω at x Ω. What happens at crossings and corners? If x is a corner or a crossing singularity (or of similar nature), Σ x ( f ) = R 2 \ 0. x ξ (x, ξ) WF( f ) f Singularities of f 17 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
72 PSEUDODIFFERENTIAL OPERATORS ΨDO S Theorem (Pseudolocal property of ΨDO s) For any pseudodifferential operator P, we have WF(P f ) WF( f ). If P is elliptic, we even have WF(P f ) = WF( f ). 18 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
73 PSEUDODIFFERENTIAL OPERATORS ΨDO S Theorem (Pseudolocal property of ΨDO s) For any pseudodifferential operator P, we have WF(P f ) WF( f ). If P is elliptic, we even have WF(P f ) = WF( f ). Pseudodifferential operators do not create new singularities! Example: The Riesz-Potentials are elliptic pseudodifferential operators. I α f = ( ) α/2 f = F 1 ( ξ α ˆ f (ξ)) 18 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
74 PSEUDODIFFERENTIAL OPERATORS ΨDO S Theorem (Pseudolocal property of ΨDO s) For any pseudodifferential operator P, we have WF(P f ) WF( f ). If P is elliptic, we even have WF(P f ) = WF( f ). Pseudodifferential operators do not create new singularities! The result is also valid for differential operators P as they are special cases of ΨDo s. The result can be also used in order to compute WF, e.g. for Heaviside function H we have H = δ so WF(H) = WF(δ) = {0}. Example: The Riesz-Potentials are elliptic pseudodifferential operators. I α f = ( ) α/2 f = F 1 ( ξ α ˆ f (ξ)) 18 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
75 CONSEQUENCES FOR TOMOGRAPHY In Λ-tomography we reconstruct Λ f = I 1 f. By ellipticity and the pseudolocal property we get WF(Λ f ) = WF( f ). 19 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
76 CONSEQUENCES FOR TOMOGRAPHY In Λ-tomography we reconstruct Λ f = I 1 f. By ellipticity and the pseudolocal property we get Λ-tomography reconstructs all singularities of f. WF(Λ f ) = WF( f ). 19 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
77 CONSEQUENCES FOR TOMOGRAPHY In Λ-tomography we reconstruct Λ f = I 1 f. By ellipticity and the pseudolocal property we get WF(Λ f ) = WF( f ). Λ-tomography reconstructs all singularities of f. From the inversion formula f = const I 1 n R R f we get similarly WF(R R f ) = WF( f ). 19 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
78 CONSEQUENCES FOR TOMOGRAPHY In Λ-tomography we reconstruct Λ f = I 1 f. By ellipticity and the pseudolocal property we get Λ-tomography reconstructs all singularities of f. WF(Λ f ) = WF( f ). From the inversion formula f = const I 1 n R R f we get similarly WF(R R f ) = WF( f ). Although singularities are differently pronounced in f, R R f and Λ f, we have WF( f ) = WF(R R f ) = WF(Λ f ). WF does not quantify the strength of singularities. To quantify the strength of singularities, a refined notion of the wavefront set can be defined (Sobolev wavefront set). 19 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
79 CONSEQUENCES FOR TOMOGRAPHY In Λ-tomography we reconstruct Λ f = I 1 f. By ellipticity and the pseudolocal property we get Λ-tomography reconstructs all singularities of f. WF(Λ f ) = WF( f ). From the inversion formula f = const I 1 n R R f we get similarly WF(R R f ) = WF( f ). Although singularities are differently pronounced in f, R R f and Λ f, we have WF( f ) = WF(R R f ) = WF(Λ f ). WF does not quantify the strength of singularities. To quantify the strength of singularities, a refined notion of the wavefront set can be defined (Sobolev wavefront set). For general FBP type reconstruction operators R P with general pseudodifferential operators P as filters, a minimum requirement on the filter can be formulated as In general, this is not even if P is elliptic. WF(R PR f )! = WF( f ). 19 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
80 Thanks! Next week: Microlocal analysis in limited angle tomography 20 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
LECTURES ON MICROLOCAL CHARACTERIZATIONS IN LIMITED-ANGLE
LECTURES ON MICROLOCAL CHARACTERIZATIONS IN LIMITED-ANGLE TOMOGRAPHY Jürgen Frikel 4 LECTURES 1 Today: Introduction to the mathematics of computerized tomography 2 Mathematics of computerized tomography
More informationA Walking Tour of Microlocal Analysis
A Walking Tour of Microlocal Analysis Jeff Schonert August 10, 2006 Abstract We summarize some of the basic principles of microlocal analysis and their applications. After reviewing distributions, we then
More informationANALYZING RECONSTRUCTION ARTIFACTS FROM ARBITRARY INCOMPLETE X-RAY CT DATA
ANALYZING RECONSTRUCTION ARTIFACTS FROM ARBITRARY INCOMPLETE X-RAY CT DATA LEISE BORG 1, JÜRGEN FRIKEL2, JAKOB SAUER JØRGENSEN 3, AND ERIC TODD QUINTO 4 Abstract. This article provides a mathematical analysis
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 informationMICROLOCAL ANALYSIS METHODS
MICROLOCAL ANALYSIS METHODS PLAMEN STEFANOV One of the fundamental ideas of classical analysis is a thorough study of functions near a point, i.e., locally. Microlocal analysis, loosely speaking, is analysis
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 informationMicrolocal Analysis : a short introduction
Microlocal Analysis : a short introduction Plamen Stefanov Purdue University Mini Course, Fields Institute, 2012 Plamen Stefanov (Purdue University ) Microlocal Analysis : a short introduction 1 / 25 Introduction
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 informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationarxiv: v4 [math.ap] 21 Jan 2015
Artifacts in incomplete data tomography with applications to photoacoustic tomography and sonar Jürgen Frikel Eric Todd Quinto arxiv:1407.3453v4 [math.ap] 21 Jan 2015 Abstract We develop a paradigm using
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 informationLOCAL SINGULARITY RECONSTRUCTION FROM INTEGRALS OVER CURVES IN R 3. Eric Todd Quinto. Hans Rullgård. (Communicated by Margaret Cheney)
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X Website: http://aimsciences.org pp. X XX LOCAL SINGULARITY RECONSTRUCTION FROM INTEGRALS OVER CURVES IN R 3 Eric Todd Quinto Tufts University
More informationSPREADING OF LAGRANGIAN REGULARITY ON RATIONAL INVARIANT TORI
SPREADING OF LAGRANGIAN REGULARITY ON RATIONAL INVARIANT TORI JARED WUNSCH Abstract. Let P h be a self-adjoint semiclassical pseudodifferential operator on a manifold M such that the bicharacteristic flow
More informationc 2015 Society for Industrial and Applied Mathematics
SIAM J. APPL. MATH. Vol. 75, No. 2, pp. 703 725 c 2015 Society for Industrial and Applied Mathematics ARTIFACTS IN INCOMPLETE DATA TOMOGRAPHY WITH APPLICATIONS TO PHOTOACOUSTIC TOMOGRAPHY AND SONAR JÜRGEN
More informationAnalyzing Reconstruction Artifacts from Arbitrary Incomplete X-ray CT Data \ast
SIAM J. IMAGING SCIENCES Vol. 11, No. 4, pp. 2786--2814 c\bigcirc 2018 Society for Industrial and Applied Mathematics Analyzing Reconstruction Artifacts from Arbitrary Incomplete X-ray CT Data \ast Leise
More informationBrief notes on Fourier Transform MSRI 18th 24th July 2009
Brief notes on Fourier Transform MSRI 18th 24th July 2009 P. Kuchment Contents 1 Some linear algebra 2 2 Fourier series 2 2.1 Fourier series expansions....................... 2 2.2 Properties of Fourier
More informationElliptic Regularity. Throughout we assume all vector bundles are smooth bundles with metrics over a Riemannian manifold X n.
Elliptic Regularity Throughout we assume all vector bundles are smooth bundles with metrics over a Riemannian manifold X n. 1 Review of Hodge Theory In this note I outline the proof of the following Fundamental
More informationIntroduction to Microlocal Analysis
Introduction to Microlocal Analysis First lecture: Basics Dorothea Bahns (Göttingen) Third Summer School on Dynamical Approaches in Spectral Geometry Microlocal Methods in Global Analysis University of
More informationBoundary problems for fractional Laplacians
Boundary problems for fractional Laplacians Gerd Grubb Copenhagen University Spectral Theory Workshop University of Kent April 14 17, 2014 Introduction The fractional Laplacian ( ) a, 0 < a < 1, has attracted
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 informationCHAPTER 5. Microlocalization
CHAPTER 5 Microlocalization 5.1. Calculus of supports Recall that we have already defined the support of a tempered distribution in the slightly round-about way: (5.1) if u S (R n ), supp(u) = {x R n ;
More informationInterpolation procedure in filtered backprojection algorithm for the limited-angle tomography
Interpolation procedure in filtered backprojection algorithm for the limited-angle tomography Aleksander Denisiuk 1 Abstract A new data completion procedure for the limited-angle tomography is proposed.
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 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 informationI. DISTRIBUTION SPACES
1.1 I. DISTIBUTION SPACES 1. Motivation and overview 1.1. Introduction. In the study of ordinary differential equations one can get very far by using just the classical concept of differentiability, working
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 informationLecture 2: Some basic principles of the b-calculus
Lecture 2: Some basic principles of the b-calculus Daniel Grieser (Carl von Ossietzky Universität Oldenburg) September 20, 2012 Summer School Singular Analysis Daniel Grieser (Oldenburg) Lecture 2: Some
More informationThe Schrödinger propagator for scattering metrics
The Schrödinger propagator for scattering metrics Andrew Hassell (Australian National University) joint work with Jared Wunsch (Northwestern) MSRI, May 5-9, 2003 http://arxiv.org/math.ap/0301341 1 Schrödinger
More informationPOINTWISE BOUNDS ON QUASIMODES OF SEMICLASSICAL SCHRÖDINGER OPERATORS IN DIMENSION TWO
POINTWISE BOUNDS ON QUASIMODES OF SEMICLASSICAL SCHRÖDINGER OPERATORS IN DIMENSION TWO HART F. SMITH AND MACIEJ ZWORSKI Abstract. We prove optimal pointwise bounds on quasimodes of semiclassical Schrödinger
More informationThe wave front set of a distribution
The wave front set of a distribution The Fourier transform of a smooth compactly supported function u(x) decays faster than any negative power of the dual variable ξ; that is for every number N there exists
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 informationSaturation Rates for Filtered Back Projection
Saturation Rates for Filtered Back Projection Matthias Beckmann, Armin Iske Fachbereich Mathematik, Universität Hamburg 3rd IM-Workshop on Applied Approximation, Signals and Images Bernried, 21st February
More informationApplication of local operators for numerical reconstruction of the singular support of a vector field by its known ray transforms
Application of local operators for numerical reconstruction of the singular support of a vector field by its known ray transforms E Yu Derevtsov 1, V V Pickalov 2 and T Schuster 3 1 Sobolev Institute of
More informationAn Introduction to X-ray tomography and Radon Transforms
Proceedings of Symposia in Applied Mathematics An Introduction to X-ray tomography and Radon Transforms Eric Todd Quinto Abstract. This article provides an introduction to the mathematics behind X-ray
More informationWavelets and multiresolution representations. Time meets frequency
Wavelets and multiresolution representations Time meets frequency Time-Frequency resolution Depends on the time-frequency spread of the wavelet atoms Assuming that ψ is centred in t=0 Signal domain + t
More informationMaximal Entropy for Reconstruction of Back Projection Images
Maximal Entropy for Reconstruction of Back Projection Images Tryphon Georgiou Department of Electrical and Computer Engineering University of Minnesota Minneapolis, MN 5545 Peter J Olver Department of
More informationIntroduction to analysis on manifolds with corners
Introduction to analysis on manifolds with corners Daniel Grieser (Carl von Ossietzky Universität Oldenburg) June 19, 20 and 21, 2017 Summer School and Workshop The Sen conjecture and beyond, UCL Daniel
More informationMORERA THEOREMS VIA MICROLOCAL ANALYSIS. Josip Globevnik and Eric Todd Quinto
MORERA THEOREMS VIA MICROLOCAL ANALYSIS Josip Globevnik and Eric Todd Quinto Abstract. We prove Morera theorems for curves in the plane using microlocal analysis. The key is that microlocal smoothness
More informationIntroduction. Christophe Prange. February 9, This set of lectures is motivated by the following kind of phenomena:
Christophe Prange February 9, 206 This set of lectures is motivated by the following kind of phenomena: sin(x/ε) 0, while sin 2 (x/ε) /2. Therefore the weak limit of the product is in general different
More informationNOTES FOR CARDIFF LECTURES ON MICROLOCAL ANALYSIS
NOTES FOR CARDIFF LECTURES ON MICROLOCAL ANALYSIS JARED WUNSCH Note that these lectures overlap with Alex s to a degree, to ensure a smooth handoff between lecturers! Our notation is mostly, but not completely,
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 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 informationPSEUDODIFFERENTIAL OPERATORS WITH GENERALIZED SYMBOLS AND REGULARITY THEORY
Electronic Journal of Differential Equations, Vol. 2005(2005), No. 116, pp. 1 43. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) PSEUDODIFFERENTIAL
More informationInverse Problem in Quantitative Susceptibility Mapping From Theory to Application. Jae Kyu Choi.
Inverse Problem in Quantitative Susceptibility Mapping From Theory to Application Jae Kyu Choi jaycjk@yonsei.ac.kr Department of CSE, Yonsei University 2015.08.07 Hangzhou Zhe Jiang University Jae Kyu
More informationA new class of pseudodifferential operators with mixed homogenities
A new class of pseudodifferential operators with mixed homogenities Po-Lam Yung University of Oxford Jan 20, 2014 Introduction Given a smooth distribution of hyperplanes on R N (or more generally on a
More informationMicrolocal analysis of. generalized functions: pseudodifferential techniques and propagation of singularities.
Loughborough University Institutional Repository Microlocal analysis of generalized functions: pseudodifferential techniques and propagation of singularities This item was submitted to Loughborough University's
More informationEQUIVARIANT AND FRACTIONAL INDEX OF PROJECTIVE ELLIPTIC OPERATORS. V. Mathai, R.B. Melrose & I.M. Singer. Abstract
j. differential geometry x78 (2008) 465-473 EQUIVARIANT AND FRACTIONAL INDEX OF PROJECTIVE ELLIPTIC OPERATORS V. Mathai, R.B. Melrose & I.M. Singer Abstract In this note the fractional analytic index,
More informationLECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI
LECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI Georgian Technical University Tbilisi, GEORGIA 0-0 1. Formulation of the corresponding
More informationThe Mathematics of Computerized Tomography
The Mathematics of Computerized Tomography The Mathematics of Computerized Tomography F. Natterer University of Münster Federal Republic of Germany B. G. TEUBNER Stuttgart @) JOHN WILEY & SONS Chichester.
More informationDiffraction by Edges. András Vasy (with Richard Melrose and Jared Wunsch)
Diffraction by Edges András Vasy (with Richard Melrose and Jared Wunsch) Cambridge, July 2006 Consider the wave equation Pu = 0, Pu = D 2 t u gu, on manifolds with corners M; here g 0 the Laplacian, D
More informationApplications of microlocal analysis to some hyperbolic inverse problems
Purdue University Purdue e-pubs Open Access Dissertations Theses and Dissertations Spring 2015 Applications of microlocal analysis to some hyperbolic inverse problems Andrew J Homan Purdue University Follow
More information数理解析研究所講究録別冊 = RIMS Kokyuroku Bessa (2008), B5:
Remarks on the Kernel Theorems in H TitleAnalysis and the Exact WKB Analysis Differential Equations) Author(s) LIESS, Otto; OKADA, Yasunori Citation 数理解析研究所講究録別冊 = RIMS Kokyuroku Bessa (2008), B5: 199-208
More informatione (x y)2 /4kt φ(y) dy, for t > 0. (4)
Math 24A October 26, 2 Viktor Grigoryan Heat equation: interpretation of the solution Last time we considered the IVP for the heat equation on the whole line { ut ku xx = ( < x
More informationSome uniqueness results for the determination of forces acting over Germain-Lagrange plates
Some uniqueness results for the determination of forces acting over Germain-Lagrange plates Three different techniques A. Kawano Escola Politecnica da Universidade de Sao Paulo A. Kawano (Poli-USP) Germain-Lagrange
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 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 informationMicrolocal analysis and inverse problems Lecture 4 : Uniqueness results in admissible geometries
Microlocal analysis and inverse problems Lecture 4 : Uniqueness results in admissible geometries David Dos Santos Ferreira LAGA Université de Paris 13 Wednesday May 18 Instituto de Ciencias Matemáticas,
More informationMicrolocal Analysis of Elliptical Radon Transforms with Foci on a Line
Microlocal Analysis of Elliptical Radon Transforms with Foci on a Line Venkateswaran P. Krishnan, Howard Levinson and Eric Todd Quinto We dedicate this article to the memory of Leon Ehrenpreis, a brilliant
More informationSEMICLASSICAL LAGRANGIAN DISTRIBUTIONS
SEMICLASSICAL LAGRANGIAN DISTRIBUTIONS SEMYON DYATLOV Abstract. These are (rather sloppy) notes for the talk given in UC Berkeley in March 2012, attempting to explain the global theory of semiclassical
More informationSharp estimates for a class of hyperbolic pseudo-differential equations
Results in Math., 41 (2002), 361-368. Sharp estimates for a class of hyperbolic pseudo-differential equations Michael Ruzhansky Abstract In this paper we consider the Cauchy problem for a class of hyperbolic
More informationLocal Tomography in 3-D SPECT
Pubblicazioni del Centro De Giorgi Proceedings of the Interdisciplinary Workshop on Mathematical Methods in Biomedical Imaging and Intensity-Modulated Radiation Therapy (IMRT), Pisa, Italy, October 2007.
More informationMesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain
Mesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain Stephen Edward Moore Johann Radon Institute for Computational and Applied Mathematics Austrian Academy of Sciences,
More informationEXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018
EXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018 While these notes are under construction, I expect there will be many typos. The main reference for this is volume 1 of Hörmander, The analysis of liner
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationWave equation on manifolds and finite speed of propagation
Wave equation on manifolds and finite speed of propagation Ethan Y. Jaffe Let M be a Riemannian manifold (without boundary), and let be the (negative of) the Laplace-Beltrami operator. In this note, we
More informationIndex theory on manifolds with corners: Generalized Gauss-Bonnet formulas
Index theory on singular manifolds I p. 1/4 Index theory on singular manifolds I Index theory on manifolds with corners: Generalized Gauss-Bonnet formulas Paul Loya Index theory on singular manifolds I
More informationTravel Time Tomography and Tensor Tomography, I
Travel Time Tomography and Tensor Tomography, I Plamen Stefanov Purdue University Mini Course, MSRI 2009 Plamen Stefanov (Purdue University ) Travel Time Tomography and Tensor Tomography, I 1 / 17 Alternative
More informationFourier Transform & Sobolev Spaces
Fourier Transform & Sobolev Spaces Michael Reiter, Arthur Schuster Summer Term 2008 Abstract We introduce the concept of weak derivative that allows us to define new interesting Hilbert spaces the Sobolev
More informationA local estimate from Radon transform and stability of Inverse EIT with partial data
A local estimate from Radon transform and stability of Inverse EIT with partial data Alberto Ruiz Universidad Autónoma de Madrid ge Joint work with P. Caro (U. Helsinki) and D. Dos Santos Ferreira (Paris
More informationMathematical Methods of Computed Tomography
Mathematical Methods of Computed Tomography Peter Kuchment Dedicated to the memory of Leon Ehrenpreis, Larry Shepp and Iosif Shneiberg friends, great mathematicians and human beings Young Researchers School
More information1 Radon Transform and X-Ray CT
Radon Transform and X-Ray CT In this section we look at an important class of imaging problems, the inverse Radon transform. We have seen that X-ray CT (Computerized Tomography) involves the reconstruction
More informationMATH 220: THE FOURIER TRANSFORM TEMPERED DISTRIBUTIONS. Beforehand we constructed distributions by taking the set Cc
MATH 220: THE FOURIER TRANSFORM TEMPERED DISTRIBUTIONS Beforehand we constructed distributions by taking the set Cc ( ) as the set of very nice functions, and defined distributions as continuous linear
More informationSELF-ADJOINTNESS OF SCHRÖDINGER-TYPE OPERATORS WITH SINGULAR POTENTIALS ON MANIFOLDS OF BOUNDED GEOMETRY
Electronic Journal of Differential Equations, Vol. 2003(2003), No.??, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) SELF-ADJOINTNESS
More informationMicrolocal Analysis of an Ultrasound Transform with Circular Source and Receiver Trajectories
Contemporary Mathematics Microlocal Analysis of an Ultrasound Transform with Circular Source and Receiver Trajectories G. Ambartsoumian, J. Boman, V. P. Krishnan, and E. T. Quinto Abstract. We consider
More informationHölder regularity estimation by Hart Smith and Curvelet transforms
Hölder regularity estimation by Hart Smith and Curvelet transforms Jouni Sampo Lappeenranta University Of Technology Department of Mathematics and Physics Finland 18th September 2007 This research is done
More informationMultiplication of Distributions
Multiplication of Distributions Christian Brouder Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie UPMC, Paris quantum field theory Feynman diagram x 2 x 3 x 1 x 4 x 5 x 6 G x 7 Feynman
More informationMicrolocal analysis and global solutions of some hyperbolic equations with discontinuous coefficients
Microlocal analysis and global solutions of some hyperbolic equations with discontinuous coefficients G. Hörmann (guenther@dix.mines.edu) and M. V. de Hoop (mdehoop@dix.mines.edu) Center for Wave Phenomena
More informationReminder Notes for the Course on Distribution Theory
Reminder Notes for the Course on Distribution Theory T. C. Dorlas Dublin Institute for Advanced Studies School of Theoretical Physics 10 Burlington Road, Dublin 4, Ireland. Email: dorlas@stp.dias.ie March
More informationEarlier this week we defined Boolean atom: A Boolean atom a is a nonzero element of a Boolean algebra, such that ax = a or ax = 0 for all x.
Friday, April 20 Today we will continue in Course Notes 3.3: Abstract Boolean Algebras. Earlier this week we defined Boolean atom: A Boolean atom a is a nonzero element of a Boolean algebra, such that
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 informationA SUPPORT THEOREM FOR THE GEODESIC RAY TRANSFORM OF SYMMETRIC TENSOR FIELDS. Venkateswaran P. Krishnan. Plamen Stefanov
A SUPPORT THEOREM FOR THE GEODESIC RAY TRANSFORM OF SYMMETRIC TENSOR FIELDS Venkateswaran P. Krishnan 110, 8 th Street, Rensselaer Polytechnic Institute Troy, NY 12180. Plamen Stefanov Department of Mathematics
More informationIntroduction to the Boundary Element Method
Introduction to the Boundary Element Method Salim Meddahi University of Oviedo, Spain University of Trento, Trento April 27 - May 15, 2015 1 Syllabus The Laplace problem Potential theory: the classical
More informationTHERMOACOUSTIC TOMOGRAPHY ARISING IN BRAIN IMAGING
THERMOACOUSTIC TOMOGRAPHY ARISING IN BRAIN IMAGING PLAMEN STEFANOV AND GUNTHER UHLMANN Abstract. We study the mathematical model of thermoacoustic and photoacoustic tomography when the sound speed has
More informationOn m-accretive Schrödinger operators in L p -spaces on manifolds of bounded geometry
On m-accretive Schrödinger operators in L p -spaces on manifolds of bounded geometry Ognjen Milatovic Department of Mathematics and Statistics University of North Florida Jacksonville, FL 32224 USA. Abstract
More informationMLISP: Machine Learning in Signal Processing Spring Lecture 8-9 May 4-7
MLISP: Machine Learning in Signal Processing Spring 2018 Prof. Veniamin Morgenshtern Lecture 8-9 May 4-7 Scribe: Mohamed Solomon Agenda 1. Wavelets: beyond smoothness 2. A problem with Fourier transform
More informationTHE FORM SUM AND THE FRIEDRICHS EXTENSION OF SCHRÖDINGER-TYPE OPERATORS ON RIEMANNIAN MANIFOLDS
THE FORM SUM AND THE FRIEDRICHS EXTENSION OF SCHRÖDINGER-TYPE OPERATORS ON RIEMANNIAN MANIFOLDS OGNJEN MILATOVIC Abstract. We consider H V = M +V, where (M, g) is a Riemannian manifold (not necessarily
More informationRIEMANN-HILBERT PROBLEM FOR THE MOISIL-TEODORESCU SYSTEM IN MULTIPLE CONNECTED DOMAINS
Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 310, pp. 1 5. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu RIEMANN-HILBERT PROBLEM FOR THE MOISIL-TEODORESCU
More informationOn support theorems for X-Ray transform with incomplete data
On for X-Ray transform with incomplete data Alexander Denisjuk Elblag University of Humanities and Economy Elblag, Poland denisjuk@euh-e.edu.pl November 9, 2009 1 / 40 2 / 40 Weighted X-ray Transform X
More informationLectures on. Sobolev Spaces. S. Kesavan The Institute of Mathematical Sciences, Chennai.
Lectures on Sobolev Spaces S. Kesavan The Institute of Mathematical Sciences, Chennai. e-mail: kesh@imsc.res.in 2 1 Distributions In this section we will, very briefly, recall concepts from the theory
More informationA PROPERTY OF SOBOLEV SPACES ON COMPLETE RIEMANNIAN MANIFOLDS
Electronic Journal of Differential Equations, Vol. 2005(2005), No.??, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) A PROPERTY
More informationSINGULARITY DETECTION AND PROCESSING WITH WAVELETS
SINGULARITY DETECTION AND PROCESSING WITH WAVELETS Stephane Mallat and Wen Liang Hwang Courant Institute of Mathematical Sciences New York University, New York, NY 10012 Technical Report March 1991 Abstract
More informationDefinition and basic properties of heat kernels I, An introduction
Definition and basic properties of heat kernels I, An introduction Zhiqin Lu, Department of Mathematics, UC Irvine, Irvine CA 92697 April 23, 2010 In this lecture, we will answer the following questions:
More informationSome Aspects of Solutions of Partial Differential Equations
Some Aspects of Solutions of Partial Differential Equations K. Sakthivel Department of Mathematics Indian Institute of Space Science & Technology(IIST) Trivandrum - 695 547, Kerala Sakthivel@iist.ac.in
More informationThe mathematics behind Computertomography
Radon transforms The mathematics behind Computertomography PD Dr. Swanhild Bernstein, Institute of Applied Analysis, Freiberg University of Mining and Technology, International Summer academic course 2008,
More informationCan There Be a General Theory of Fourier Integral Operators?
Can There Be a General Theory of Fourier Integral Operators? Allan Greenleaf University of Rochester Conference on Inverse Problems in Honor of Gunther Uhlmann UC, Irvine June 21, 2012 How I started working
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 informationBAND-LIMITED REFINABLE FUNCTIONS FOR WAVELETS AND FRAMELETS
BAND-LIMITED REFINABLE FUNCTIONS FOR WAVELETS AND FRAMELETS WEIQIANG CHEN AND SAY SONG GOH DEPARTMENT OF MATHEMATICS NATIONAL UNIVERSITY OF SINGAPORE 10 KENT RIDGE CRESCENT, SINGAPORE 119260 REPUBLIC OF
More informationarxiv:math/ v2 [math.ca] 18 Dec 1999
Illinois J. Math., to appear Non-symmetric convex domains have no basis of exponentials Mihail N. Kolountzakis 13 December 1998; revised October 1999 arxiv:math/9904064v [math.ca] 18 Dec 1999 Abstract
More informationComputational Harmonic Analysis (Wavelet Tutorial) Part II
Computational Harmonic Analysis (Wavelet Tutorial) Part II Understanding Many Particle Systems with Machine Learning Tutorials Matthew Hirn Michigan State University Department of Computational Mathematics,
More informationMicrolocal Analysis of an Ultrasound Transform with Circular Source and Receiver Trajectories
Contemporary Mathematics Microlocal Analysis of an Ultrasound Transform with Circular Source and Receiver Trajectories G. Ambartsoumian, J. Boman, V. P. Krishnan, and E. T. Quinto This article is dedicated
More information