LECTURES ON MICROLOCAL CHARACTERIZATIONS IN LIMITED-ANGLE
|
|
- Asher George
- 5 years ago
- Views:
Transcription
1 LECTURES ON MICROLOCAL CHARACTERIZATIONS IN LIMITED-ANGLE TOMOGRAPHY Jürgen Frikel
2 4 LECTURES 1 Today: Introduction to the mathematics of computerized tomography 2 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
3 4 LECTURES 1 Today: Introduction to the mathematics of computerized tomography 2 Nov. 18: Introduction to the basic concepts of microlocal analysis 2 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
4 4 LECTURES 1 Today: Introduction to the mathematics of computerized tomography 2 Nov. 18: Introduction to the basic concepts of microlocal analysis 3 Nov. 25: Microlocal analysis of limited angle reconstructions in tomography I 2 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
5 4 LECTURES 1 Today: Introduction to the mathematics of computerized tomography 2 Nov. 18: 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)
6 4 LECTURES 1 Today: Introduction to the mathematics of computerized tomography 2 Nov. 18: 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: F. Natterer, The mathematics of computerized tomography. Stuttgart: B. G. Teubner, 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)
7 ORIGINS OF TOMOGRAPHY Image taken from Johan Radon, Über die Bestimmung von Funktionen durch ihre Integralwerte Längs gewisser Manningsfaltigkeiten, Berichte Sächsische Akadamie der Wissenschaften, Leipzig, Math.-Phys. Kl., 69, pp , Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
8 ORIGINS OF TOMOGRAPHY Allan Cormack Godfrey Hounsfield The Nobel Prize in Physiology or Medicine 1979 was awarded jointly to Allan M. Cormack and Godfrey N. Hounsfield for the development of computer assisted tomography 4 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
9 ORIGINS OF TOMOGRAPHY First scanner, $300 Modern scanner, > $1 million 4 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
10 ORIGINS OF TOMOGRAPHY First clinical scan 1971 Modern scan Case courtesy of Dr Maxime St-Amant, Radiopaedia.org 4 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
11 TODAY Analytical approach to computerized tomography Principle of tomography Radon transform - a mathematical model of tomography Reconstruction via backprojection Fourier slice theorem Inversion formulas & Filtered backprojection Ill-posedness & regularization 5 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
12 PRINCIPLE OF COMPUTED TOMOGRAPHY Beer Lambert law X-rays are attenuated when traveling through object according to di = f (γ(t)) I(t) dt for t R I(0) = I 0 Source f I 0 f = attenuation coefficient, γ = x-ray path γ(t) = s θ + t θ = line with direction θ starting at x detector = sθ Detector I(θ, p) 6 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
13 PRINCIPLE OF COMPUTED TOMOGRAPHY Beer Lambert law X-rays are attenuated when traveling through object according to di = f (γ(t)) I(t) dt for t R I(0) = I 0 Source f I 0 f = attenuation coefficient, γ = x-ray path γ(t) = s θ + t θ = line with direction θ starting at x detector = sθ Detector I(θ, p) Solution of the initial value problem is given by { tdetector } I(θ, s) = I 0 exp f (s θ + t θ ) dt 0 6 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
14 PRINCIPLE OF COMPUTED TOMOGRAPHY Beer Lambert law X-rays are attenuated when traveling through object according to di = f (γ(t)) I(t) dt for t R I(0) = I 0 Source f I 0 f = attenuation coefficient, γ = x-ray path γ(t) = s θ + t θ = line with direction θ starting at x detector = sθ Detector I(θ, p) Solution of the initial value problem is given by { tdetector } I(θ, s) = I 0 exp f (s θ + t θ ) dt 0 { } = I 0 exp f (x) dx. L(θ,s) 6 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
15 PRINCIPLE OF COMPUTED TOMOGRAPHY Beer Lambert law X-rays are attenuated when traveling through object according to di = f (γ(t)) I(t) dt for t R I(0) = I 0 Source f I 0 f = attenuation coefficient, γ = x-ray path γ(t) = s θ + t θ = line with direction θ starting at x detector = sθ Detector I(θ, p) Solution of the initial value problem is given by { tdetector } I(θ, s) = I 0 exp f (s θ + t θ ) dt 0 { } = I 0 exp f (x) dx. L(θ,s) Mathematical model of the measurement process ( ) I 0 R f (θ, s) = f (x) dx = ln L(θ,s) I(θ, s) 6 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
16 RECONSTRUCTION PROBLEM??? Data / Sinogram Sought image Mathematical problem Solve the integral equation y = R f 7 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
17 RECONSTRUCTION PROBLEM??? Data / Sinogram Sought image Mathematical problem Solve the integral equation y = R f I. Algebraic reconstruction Fully discretize formulation of the problem Linear system of equations Rx = y This is an algebraic problem! Examples: ART, SART, SIRT, statistical reconstruction methods such as ML-EM, variational methods, etc. 7 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
18 RECONSTRUCTION PROBLEM??? Data / Sinogram Sought image Mathematical problem Solve the integral equation y = R f I. Algebraic reconstruction Fully discretize formulation of the problem Linear system of equations Rx = y This is an algebraic problem! Examples: ART, SART, SIRT, statistical reconstruction methods such as ML-EM, variational methods, etc. II. Analytical reconstruction Study of the continuous problem (above) Derivation of inversion formulas Discretization of analytic reconstruction formulas Examples: Filtered backprojection (FBP), Fourier inversion, etc. 7 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
19 RECONSTRUCTION PROBLEM??? Data / Sinogram Sought image Mathematical problem Solve the integral equation y = R f I. Algebraic reconstruction Fully discretize formulation of the problem Linear system of equations Rx = y This is an algebraic problem! Examples: ART, SART, SIRT, statistical reconstruction methods such as ML-EM, variational methods, etc. II. Analytical reconstruction Study of the continuous problem (above) Derivation of inversion formulas Discretization of analytic reconstruction formulas Examples: Filtered backprojection (FBP), Fourier inversion, etc. 7 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
20 RADON TRANSFORM IN R n Definition Let f : R n R be a suitably chosen function. The Radon transform of f, denoted by R f, is defined as (1) R f (θ, s) = f (x) dσ(x), (θ, s) S n 1 R, H(θ,s) where H(θ, s) = { x R n : x, θ = s } is the hyperplane with normal vector θ and the signed distance from the origin s. 8 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
21 RADON TRANSFORM IN R n Definition Let f : R n R be a suitably chosen function. The Radon transform of f, denoted by R f, is defined as (1) R f (θ, s) = f (x) dσ(x), (θ, s) S n 1 R, H(θ,s) where H(θ, s) = { x R n : x, θ = s } is the hyperplane with normal vector θ and the signed distance from the origin s. Notation: R θ f (s) = R f (θ, s). This function called projections along direction θ. 8 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
22 RADON TRANSFORM IN R n Definition Let f : R n R be a suitably chosen function. The Radon transform of f, denoted by R f, is defined as (1) R f (θ, s) = f (x) dσ(x), (θ, s) S n 1 R, H(θ,s) where H(θ, s) = { x R n : x, θ = s } is the hyperplane with normal vector θ and the signed distance from the origin s. Notation: R θ f (s) = R f (θ, s). This function called projections along direction θ. Radon transform maps functions to its hyperplane itegrals, 8 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
23 RADON TRANSFORM IN R n Definition Let f : R n R be a suitably chosen function. The Radon transform of f, denoted by R f, is defined as (1) R f (θ, s) = f (x) dσ(x), (θ, s) S n 1 R, H(θ,s) where H(θ, s) = { x R n : x, θ = s } is the hyperplane with normal vector θ and the signed distance from the origin s. Notation: R θ f (s) = R f (θ, s). This function called projections along direction θ. Radon transform maps functions to its hyperplane itegrals, Radon transform is an even function R f ( θ, s) = R f (θ, s) Measurment only w.r.t. half-sphere 8 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
24 RADON TRANSFORM IN R n Definition Let f : R n R be a suitably chosen function. The Radon transform of f, denoted by R f, is defined as (1) R f (θ, s) = f (x) dσ(x), (θ, s) S n 1 R, H(θ,s) where H(θ, s) = { x R n : x, θ = s } is the hyperplane with normal vector θ and the signed distance from the origin s. Notation: R θ f (s) = R f (θ, s). This function called projections along direction θ. Radon transform maps functions to its hyperplane itegrals, Radon transform is an even function R f ( θ, s) = R f (θ, s) Measurment only w.r.t. half-sphere Continuity of the Radon transform = Stability of the measurement process R is continuous on many standard function spaces, such as L 1 (R n ), L 2 (Ω), S(R n ) and many distributional spaces. 8 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
25 RADON TRANSFORM OF RADIAL FUNCTIONS f : R 2 R is radial if the value f (x) only depends on x, i.e., if there is a function ϕ : [0, ) R such that f (x) = ϕ( x ). 9 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
26 RADON TRANSFORM OF RADIAL FUNCTIONS f : R 2 R is radial if the value f (x) only depends on x, i.e., if there is a function ϕ : [0, ) R such that f (x) = ϕ( x ). In this case we have R f (θ, s) = R = f (sθ + tθ ) dt ϕ( s 2 + t 2 ) dt 9 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
27 RADON TRANSFORM OF RADIAL FUNCTIONS f : R 2 R is radial if the value f (x) only depends on x, i.e., if there is a function ϕ : [0, ) R such that f (x) = ϕ( x ). In this case we have R f (θ, s) = R = = 2 0 f (sθ + tθ ) dt ϕ( s 2 + t 2 ) dt ϕ( s 2 + t 2 ) dt 9 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
28 RADON TRANSFORM OF RADIAL FUNCTIONS f : R 2 R is radial if the value f (x) only depends on x, i.e., if there is a function ϕ : [0, ) R such that f (x) = ϕ( x ). In this case we have R f (θ, s) = R = = 2 0 f (sθ + tθ ) dt ϕ( s 2 + t 2 ) dt ϕ( s 2 + t 2 ) dt substitution r 2 = s 2 + t 2 gives 2r dr = 2t dt, and hence ϕ(r)r = 2 s r 2 s dr 2 9 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
29 RADON TRANSFORM OF RADIAL FUNCTIONS f : R 2 R is radial if the value f (x) only depends on x, i.e., if there is a function ϕ : [0, ) R such that f (x) = ϕ( x ). In this case we have R f (θ, s) = R = = 2 0 f (sθ + tθ ) dt ϕ( s 2 + t 2 ) dt ϕ( s 2 + t 2 ) dt substitution r 2 = s 2 + t 2 gives 2r dr = 2t dt, and hence ϕ(r)r = 2 s r 2 s dr 2 The Radon transform of radial functions is independent of θ! 1 Projection enough to reconstruct a radial function 9 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
30 RADON TRANSFORM OF RADIAL FUNCTIONS Example: f (x) = χ B(0,1) (x) = χ [0,1] ( x ): R f (θ, s) = 2 s χ [0,1] (r)r r 2 s 2 dr 10 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
31 RADON TRANSFORM OF RADIAL FUNCTIONS Example: f (x) = χ B(0,1) (x) = χ [0,1] ( x ): R f (θ, s) = 2 s χ [0,1] (r)r r 2 s 2 dr s > 1 : R f (θ, s) = 0 10 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
32 RADON TRANSFORM OF RADIAL FUNCTIONS Example: f (x) = χ B(0,1) (x) = χ [0,1] ( x ): R f (θ, s) = 2 s χ [0,1] (r)r r 2 s 2 dr s > 1 : R f (θ, s) = 0 1 s < 1 : R f (θ, s) = 2 s r r 2 s dr = 2[ r 1/2] 1 s 2 = 2 1 s Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
33 RADON TRANSFORM OF RADIAL FUNCTIONS Example: f (x) = χ B(0,1) (x) = χ [0,1] ( x ): R f (θ, s) = 2 s χ [0,1] (r)r r 2 s 2 dr s > 1 : R f (θ, s) = 0 1 s < 1 : R f (θ, s) = 2 s r r 2 s dr = 2[ r 1/2] 1 s 2 = 2 1 s So 2 1 s 2, s < 1 R f (θ, s) = 0, otherwise R θ f (s) s 10 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
34 BACKPROJECTION OPERATOR Definition Let g be a (sinogram) function on S n 1 R. Given a projection along the direction θ, we define the backprojection operators along direction θ via Rθ g(x) = g(θ, x θ). 11 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
35 BACKPROJECTION OPERATOR Definition Let g be a (sinogram) function on S n 1 R. Given a projection along the direction θ, we define the backprojection operators along direction θ via Rθ g(x) = g(θ, x θ). The backprojection operator R as follows: R g(x) = R S n 1 θg(x) dσ(θ) = g(θ, x θ) dσ(θ). S n 1 11 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
36 BACKPROJECTION OPERATOR Definition Let g be a (sinogram) function on S n 1 R. Given a projection along the direction θ, we define the backprojection operators along direction θ via Rθ g(x) = g(θ, x θ). The backprojection operator R as follows: R g(x) = R S n 1 θg(x) dσ(θ) = g(θ, x θ) dσ(θ). S n 1 The value R g(x) is an average of all measurements g = R f (θ, s) which correspond to hyperplanes passing through the point x R n. 11 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
37 BACKPROJECTION OPERATOR Definition Let g be a (sinogram) function on S n 1 R. Given a projection along the direction θ, we define the backprojection operators along direction θ via Rθ g(x) = g(θ, x θ). The backprojection operator R as follows: R g(x) = R S n 1 θg(x) dσ(θ) = g(θ, x θ) dσ(θ). S n 1 The value R g(x) is an average of all measurements g = R f (θ, s) which correspond to hyperplanes passing through the point x R n. The operator R is the L 2 Hilbert space adjoint of the Radon transform R, i.e., for f L 2 (Ω) and g L 2 (S n 1 R) R f (θ, s)g(θ, s) dθ ds = f (x)r g(x) dx. R R n S n 1 11 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
38 BACKPROJECTION OPERATOR Definition Let g be a (sinogram) function on S n 1 R. Given a projection along the direction θ, we define the backprojection operators along direction θ via Rθ g(x) = g(θ, x θ). The backprojection operator R as follows: R g(x) = R S n 1 θg(x) dσ(θ) = g(θ, x θ) dσ(θ). S n 1 The value R g(x) is an average of all measurements g = R f (θ, s) which correspond to hyperplanes passing through the point x R n. The operator R is the L 2 Hilbert space adjoint of the Radon transform R, i.e., for f L 2 (Ω) and g L 2 (S n 1 R) R f (θ, s)g(θ, s) dθ ds = f (x)r g(x) dx. R R n S n 1 As an adjoint operator, R is continuous whenever R is 11 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
39 RECONSTRUCTION VIA BACKPROJECTION Theorem For f S(R n ) we have R R f (x) = S n 2 ( f 1 ) (x) 12 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
40 RECONSTRUCTION VIA BACKPROJECTION Theorem For f S(R n ) we have R R f (x) = S n 2 ( f 1 ) (x) Original Backprojection reconstruction 12 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
41 RECONSTRUCTION VIA BACKPROJECTION Theorem For f S(R n ) we have R R f (x) = S n 2 ( f 1 ) (x) Backprojection reconstruction produces a function that is smoother than the original function (by 1 order): high frequencies are damped 12 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
42 RECONSTRUCTION VIA BACKPROJECTION Theorem For f S(R n ) we have R R f (x) = S n 2 ( f 1 ) (x) Backprojection reconstruction produces a function that is smoother than the original function (by 1 order): high frequencies are damped Since R is a continuous operator, the backprojection reconstruction is stable 12 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
43 RECONSTRUCTION VIA BACKPROJECTION Theorem For f S(R n ) we have R R f (x) = S n 2 ( f 1 ) (x) Backprojection reconstruction produces a function that is smoother than the original function (by 1 order): high frequencies are damped Since R is a continuous operator, the backprojection reconstruction is stable To reconstruct the original function f, we need to restore (amplify) high frequencies filtering (sharpening) step 12 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
44 RECONSTRUCTION VIA BACKPROJECTION Theorem For f S(R n ) we have R R f (x) = S n 2 ( f 1 ) (x) Backprojection reconstruction produces a function that is smoother than the original function (by 1 order): high frequencies are damped Since R is a continuous operator, the backprojection reconstruction is stable To reconstruct the original function f, we need to restore (amplify) high frequencies filtering (sharpening) step Using the above theorem, the normal equation R R f = R g (up to a constant) reads f 1 = R g, for data g = R f tomographic reconstruction can be interpreted as a deconvolution problem 12 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
45 FOURIER TRANSFORM Fourier transform turns out to be a very useful tool for studying the Radon transform. Definition (Fourier transform and its inverse) Let f L 1 (R n ). The Fourier transform of f is defined via F f (ξ) ˆ f (ξ) = (2π) n/2 R n f (x)e ix ξ dx. Let g L 1 (R n ). The inverse Fourier transform of f is defined via F 1 g(x) ǧ(x) = (2π) n/2 R n f (ξ)e ix ξ dξ. 13 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
46 FOURIER TRANSFORM Fourier transform turns out to be a very useful tool for studying the Radon transform. Definition (Fourier transform and its inverse) Let f L 1 (R n ). The Fourier transform of f is defined via F f (ξ) ˆ f (ξ) = (2π) n/2 R n f (x)e ix ξ dx. Let g L 1 (R n ). The inverse Fourier transform of f is defined via F 1 g(x) ǧ(x) = (2π) n/2 R n f (ξ)e ix ξ dξ. Sometimes it s useful to calculate the 1D Fourier transform of the projection function g(θ, s) = R f (θ, s) with respect to the second variable s. To make that clear, we will write F s g(θ, σ) = (2π) 1/2 g(θ, s)e is σ ds R for the Fourier transform of g(θ, s) with respect to the variable s (here θ is considered to be a fixed parameter). Whenever we write R f (θ, σ) the Fourier transform of R f has to be understood in that sense. Same holds for the inverse Fourier transform. 13 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
47 RELATION TO THE FOURIER TRANSFORM Theorem (Fourier slice theorem) Let f S(R n ). Then, for σ R, F s R f (θ, σ) = (2π) (n 1)/2 f (σθ). 14 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
48 RELATION TO THE FOURIER TRANSFORM Theorem (Fourier slice theorem) Let f S(R n ). Then, for σ R, F s R f (θ, σ) = (2π) (n 1)/2 f (σθ). f 1D Fourier transform fˆ θ R θ f 14 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
49 RELATION TO THE FOURIER TRANSFORM Theorem (Fourier slice theorem) Let f S(R n ). Then, for σ R, F s R f (θ, σ) = (2π) (n 1)/2 f (σθ). f 1D Fourier transform fˆ θ R θ f This theorem can be used derive a reconstruction procedure Fourier reconstructions 14 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
50 INJECTIVITY OF THE RADON TRANSFORM Many proprities of the Radon transform can be derived from the properties of the Fourier transform. Theorem The Radon transform R : L 1 (R n ) L 1 (S n 1 R) is an injective operator. 15 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
51 INJECTIVITY OF THE RADON TRANSFORM Many proprities of the Radon transform can be derived from the properties of the Fourier transform. Theorem The Radon transform R : L 1 (R n ) L 1 (S n 1 R) is an injective operator. Proof. Suppose R f 0 for f L 1 (R n ). Then, the Fourier slice theorem implies f (σθ) = (2π) (1 n)/2 F s R f (θ, σ) = 0 for all (θ, σ) S n 1 R. Hence, f 0 and the injectivity of the Fourier transform implies that f Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
52 INJECTIVITY OF THE RADON TRANSFORM Many proprities of the Radon transform can be derived from the properties of the Fourier transform. Theorem The Radon transform R : L 1 (R n ) L 1 (S n 1 R) is an injective operator. So far: There is a unique solution to the reconstruction problem R f = y, formally given as f = R 1 y. 15 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
53 INJECTIVITY OF THE RADON TRANSFORM Many proprities of the Radon transform can be derived from the properties of the Fourier transform. Theorem The Radon transform R : L 1 (R n ) L 1 (S n 1 R) is an injective operator. So far: There is a unique solution to the reconstruction problem R f = y, formally given as f = R 1 y. Are there explicit expressions for R 1? 15 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
54 INJECTIVITY OF THE RADON TRANSFORM Many proprities of the Radon transform can be derived from the properties of the Fourier transform. Theorem The Radon transform R : L 1 (R n ) L 1 (S n 1 R) is an injective operator. So far: There is a unique solution to the reconstruction problem R f = y, formally given as f = R 1 y. Are there explicit expressions for R 1? YES! We will derive them in a minitute. 15 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
55 INJECTIVITY OF THE RADON TRANSFORM Many proprities of the Radon transform can be derived from the properties of the Fourier transform. Theorem The Radon transform R : L 1 (R n ) L 1 (S n 1 R) is an injective operator. So far: There is a unique solution to the reconstruction problem R f = y, formally given as f = R 1 y. Are there explicit expressions for R 1? YES! We will derive them in a minitute. Is R 1 continuous? 15 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
56 INJECTIVITY OF THE RADON TRANSFORM Many proprities of the Radon transform can be derived from the properties of the Fourier transform. Theorem The Radon transform R : L 1 (R n ) L 1 (S n 1 R) is an injective operator. So far: There is a unique solution to the reconstruction problem R f = y, formally given as f = R 1 y. Are there explicit expressions for R 1? YES! We will derive them in a minitute. Is R 1 continuous? Unfortuntely, NO. It can be shown that R 1 is not a continuous operator and, hence, that the reconstruction problem R f = y is ill-posed. However, the ill-posedness is mild (maybe later). 15 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
57 INVERSION FORMULAS For f S(R n ) and α < n we define the Riesz potential (which is a linear operator) via I α f = ( ) α/2 f = F 1 ( ξ α f (ξ) ). For g (S n 1 R) we analogously define the Riesz-potential with respect to the second variable Is α g(θ, s) = ( ) 2 α/2 ( s R f (θ, s) = F 1 s σ α F s g(θ, σ) ). 16 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
58 INVERSION FORMULAS For f S(R n ) and α < n we define the Riesz potential (which is a linear operator) via I α f = ( ) α/2 f = F 1 ( ξ α f (ξ) ). For g (S n 1 R) we analogously define the Riesz-potential with respect to the second variable Is α g(θ, s) = ( ) 2 α/2 ( s R f (θ, s) = F 1 s σ α F s g(θ, σ) ). Theorem Lef f S(R n ). Then, for any α < n, the following inversion formulas hold: f = 1 2 (2π)1 n I α R I α n+1 s R f. 16 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
59 VARIATIONS OF INVERSION FORMULAS n = 3 and α = 2: f = 1 8π 2 R R f 17 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
60 VARIATIONS OF INVERSION FORMULAS n = 3 and α = 2: f = 1 8π 2 R R f n = 3 and α = 0: f = 1 8π 2 R I 2 s R f = 1 8π 2 R ( 2 s)r f 17 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
61 VARIATIONS OF INVERSION FORMULAS n = 3 and α = 2: f = 1 8π 2 R R f n = 3 and α = 0: f = 1 8π 2 R I 2 s R f = 1 8π 2 R ( 2 s)r f n = 2 and α = 0: f = 1 4π R I 1 s R f = 1 4π R ( 2 s) 1/2 R f 17 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
62 VARIATIONS OF INVERSION FORMULAS n = 3 and α = 2: f = 1 8π 2 R R f n = 3 and α = 0: f = 1 8π 2 R I 2 s R f = 1 8π 2 R ( 2 s)r f n = 2 and α = 0: f = 1 4π R I 1 s R f = 1 4π R ( 2 s) 1/2 R f Structure: Filter & Backproject 17 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
63 VARIATIONS OF INVERSION FORMULAS n = 3 and α = 2: f = 1 8π 2 R R f n = 3 and α = 0: f = 1 8π 2 R I 2 s R f = 1 8π 2 R ( 2 s)r f n = 2 and α = 0: f = 1 4π R I 1 s R f = 1 4π R ( 2 s) 1/2 R f Structure: Filter & Backproject Sinogram Filtered sinogram 17 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
64 VARIATIONS OF INVERSION FORMULAS n = 3 and α = 2: f = 1 8π 2 R R f n = 3 and α = 0: f = 1 8π 2 R I 2 s R f = 1 8π 2 R ( 2 s)r f n = 2 and α = 0: f = 1 4π R I 1 s R f = 1 4π R ( 2 s) 1/2 R f Structure: Filter & Backproject Backprojection Filtered backprojection 17 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
65 INVERSION FORMULAS IN EVEN AND ODD DIMENSIONS For α = 0 we obtain a classical filtered backprojection form f = 1 2 (2π)1 n R I 1 n s R f, where the operator I 1 n s acts as a filter on the data g = R f. 18 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
66 INVERSION FORMULAS IN EVEN AND ODD DIMENSIONS For α = 0 we obtain a classical filtered backprojection form f = 1 2 (2π)1 n R I 1 n s R f, where the operator I 1 n s acts as a filter on the data g = R f. F s [Is 1 n R θ f ](σ) = σ n 1 F s R θ (σ) = (sgn(σ) σ) n 1 F s R θ (σ) 18 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
67 INVERSION FORMULAS IN EVEN AND ODD DIMENSIONS For α = 0 we obtain a classical filtered backprojection form f = 1 2 (2π)1 n R I 1 n s R f, where the operator I 1 n s acts as a filter on the data g = R f. F s [Is 1 n R θ f ](σ) = σ n 1 F s R θ (σ) = (sgn(σ) σ) n 1 F s R θ (σ) = (sgn(σ)) n 1 σ n 1 F s R θ (σ) 18 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
68 INVERSION FORMULAS IN EVEN AND ODD DIMENSIONS For α = 0 we obtain a classical filtered backprojection form f = 1 2 (2π)1 n R I 1 n s R f, where the operator I 1 n s acts as a filter on the data g = R f. F s [Is 1 n R θ f ](σ) = σ n 1 F s R θ (σ) = (sgn(σ) σ) n 1 F s R θ (σ) = (sgn(σ)) n 1 σ n 1 F s R θ (σ) = (sgn(σ)) n 1 i (n 1) i n 1 σ n 1 F s R θ (σ) 18 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
69 INVERSION FORMULAS IN EVEN AND ODD DIMENSIONS For α = 0 we obtain a classical filtered backprojection form f = 1 2 (2π)1 n R I 1 n s R f, where the operator I 1 n s acts as a filter on the data g = R f. F s [Is 1 n R θ f ](σ) = σ n 1 F s R θ (σ) = (sgn(σ) σ) n 1 F s R θ (σ) = (sgn(σ)) n 1 σ n 1 F s R θ (σ) = (sgn(σ)) n 1 i (n 1) i n 1 σ n 1 F s R θ (σ) = ( i sgn(σ)) n 1 (iσ) n 1 F s R θ (σ) 18 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
70 INVERSION FORMULAS IN EVEN AND ODD DIMENSIONS For α = 0 we obtain a classical filtered backprojection form f = 1 2 (2π)1 n R I 1 n s R f, where the operator I 1 n s acts as a filter on the data g = R f. F s [Is 1 n R θ f ](σ) = σ n 1 F s R θ (σ) because F s [ n sg](σ) = (iσ) n ĝ(σ) we have = (sgn(σ) σ) n 1 F s R θ (σ) = (sgn(σ)) n 1 σ n 1 F s R θ (σ) = (sgn(σ)) n 1 i (n 1) i n 1 σ n 1 F s R θ (σ) = ( i sgn(σ)) n 1 (iσ) n 1 F s R θ (σ) = ( i sgn(σ)) n 1 F s [ n 1 s R θ ](σ) 18 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
71 INVERSION FORMULAS IN EVEN AND ODD DIMENSIONS For α = 0 we obtain a classical filtered backprojection form f = 1 2 (2π)1 n R I 1 n s R f, where the operator I 1 n s acts as a filter on the data g = R f. F s [Is 1 n R θ f ](σ) = σ n 1 F s R θ (σ) because F s [ n sg](σ) = (iσ) n ĝ(σ) we have where H is defined in the Fourier domain via = (sgn(σ) σ) n 1 F s R θ (σ) = (sgn(σ)) n 1 σ n 1 F s R θ (σ) = (sgn(σ)) n 1 i (n 1) i n 1 σ n 1 F s R θ (σ) = ( i sgn(σ)) n 1 (iσ) n 1 F s R θ (σ) = ( i sgn(σ)) n 1 F s [ n 1 s R θ ](σ) = F s [H n 1 n 1 s R θ ](σ), Ĥg(σ) = i sgn(σ) ĝ(σ). 18 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
72 INVERSION FORMULAS IN EVEN AND ODD DIMENSIONS Now observe that for n 2 we have ( i sgn(σ)) n 1 ( 1) (n 2)/2 ( i sgn(σ)), = ( 1) (n 1)/2, for n even for n odd 19 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
73 INVERSION FORMULAS IN EVEN AND ODD DIMENSIONS Now observe that for n 2 we have ( i sgn(σ)) n 1 ( 1) (n 2)/2 ( i sgn(σ)), = ( 1) (n 1)/2, for n even for n odd Therefore f = 1 2 (2π)1 n R ( 1) (n 2)/2 H n 1 s R f (θ, x θ), for n even ( 1) (n 1)/2 n 1 s R f (θ, x θ), for n odd 19 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
74 INVERSION FORMULAS IN EVEN AND ODD DIMENSIONS Now observe that for n 2 we have ( i sgn(σ)) n 1 ( 1) (n 2)/2 ( i sgn(σ)), = ( 1) (n 1)/2, for n even for n odd Therefore f = 1 2 (2π)1 n R ( 1) (n 2)/2 H n 1 s R f (θ, x θ), for n even ( 1) (n 1)/2 n 1 s R f (θ, x θ), for n odd If n is odd, the filter is a differential operator and, hence, a local operator. As a result, the inversion formula is local, i.e., the inversion formula can be evaluated at point x if the data R f (θ, y θ) is known for all θ S n 1 and for y in a neighbourhood of x. 19 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
75 INVERSION FORMULAS IN EVEN AND ODD DIMENSIONS Now observe that for n 2 we have ( i sgn(σ)) n 1 ( 1) (n 2)/2 ( i sgn(σ)), = ( 1) (n 1)/2, for n even for n odd Therefore f = 1 2 (2π)1 n R ( 1) (n 2)/2 H n 1 s R f (θ, x θ), for n even ( 1) (n 1)/2 n 1 s R f (θ, x θ), for n odd If n is odd, the filter is a differential operator and, hence, a local operator. As a result, the inversion formula is local, i.e., the inversion formula can be evaluated at point x if the data R f (θ, y θ) is known for all θ S n 1 and for y in a neighbourhood of x. If n is even, the inversion formula is not local, since H is an integral operator Hg(s) = 1 g(t) π R s t dt H is the so-called Hilbert transform. 19 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
76 LAMBDA RECONSTRUCTION To make the reconstruction formula local for n = 2, the strategy is to have filter which is local (differential operator). 20 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
77 LAMBDA RECONSTRUCTION To make the reconstruction formula local for n = 2, the strategy is to have filter which is local (differential operator). Choosing n = 2 and α = 1 gives I 1 f = 1 4π R Is 2 R f = 1 4π R ( 2 s)r f Instead of reconstructing f we reconstruct Λ f := I 1 f. This formula is local and of filtered backprojection type. 20 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
78 LAMBDA RECONSTRUCTION To make the reconstruction formula local for n = 2, the strategy is to have filter which is local (differential operator). Choosing n = 2 and α = 1 gives I 1 f = 1 4π R Is 2 R f = 1 4π R ( 2 s)r f Instead of reconstructing f we reconstruct Λ f := I 1 f. This formula is local and of filtered backprojection type. In the implementation, the standard FBP filter σ has to be replaces by σ 2 high frequecies are amplified even more, low frequencies are damped. 20 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
79 LAMBDA RECONSTRUCTION To make the reconstruction formula local for n = 2, the strategy is to have filter which is local (differential operator). Choosing n = 2 and α = 1 gives I 1 f = 1 4π R Is 2 R f = 1 4π R ( 2 s)r f Instead of reconstructing f we reconstruct Λ f := I 1 f. This formula is local and of filtered backprojection type. In the implementation, the standard FBP filter σ has to be replaces by σ 2 high frequecies are amplified even more, low frequencies are damped. Instead of values f (x), it reconstructs the singularities of f. 20 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
80 LAMBDA RECONSTRUCTION To make the reconstruction formula local for n = 2, the strategy is to have filter which is local (differential operator). Choosing n = 2 and α = 1 gives I 1 f = 1 4π R Is 2 R f = 1 4π R ( 2 s)r f Instead of reconstructing f we reconstruct Λ f := I 1 f. This formula is local and of filtered backprojection type. In the implementation, the standard FBP filter σ has to be replaces by σ 2 high frequecies are amplified even more, low frequencies are damped. Instead of values f (x), it reconstructs the singularities of f. Original FBP Lambda 20 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
81 ILL-POSEDNESS OF COMPUTED TOMOGRAPHY Reconstruction is mildly ill-posed (of order (n-1)/2) c f H α 0 R f H α+(n 1)/2 C f H α 0 21 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
82 ILL-POSEDNESS OF COMPUTED TOMOGRAPHY Reconstruction is mildly ill-posed (of order (n-1)/2) c f H α 0 R f H α+(n 1)/2 C f H α 0 What causes instabilities in the reconstruction??? 21 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
83 ILL-POSEDNESS OF COMPUTED TOMOGRAPHY Reconstruction is mildly ill-posed (of order (n-1)/2) c f H α 0 R f H α+(n 1)/2 C f H α 0 What causes instabilities in the reconstruction??? We know R and R continuous (stable) operations. 21 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
84 ILL-POSEDNESS OF COMPUTED TOMOGRAPHY Reconstruction is mildly ill-posed (of order (n-1)/2) c f H α 0 R f H α+(n 1)/2 C f H α 0 What causes instabilities in the reconstruction??? We know R and R continuous (stable) operations. Given the inversion formula f = 1/2(2π) 1 n R Is 1 n R f the instabilities must come from filtering: 21 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
85 ILL-POSEDNESS OF COMPUTED TOMOGRAPHY Reconstruction is mildly ill-posed (of order (n-1)/2) c f H α 0 R f H α+(n 1)/2 C f H α 0 What causes instabilities in the reconstruction??? We know R and R continuous (stable) operations. Given the inversion formula f = 1/2(2π) 1 n R Is 1 n R f the instabilities must come from filtering: F s [Is 1 n R f ](θ, σ) = σ n 1 F s R f (θ, σ) = (2π) (n 1)/2 σ n 1 f ˆ(σθ) 21 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
86 ILL-POSEDNESS OF COMPUTED TOMOGRAPHY Reconstruction is mildly ill-posed (of order (n-1)/2) c f H α 0 R f H α+(n 1)/2 C f H α 0 What causes instabilities in the reconstruction??? We know R and R continuous (stable) operations. Given the inversion formula f = 1/2(2π) 1 n R Is 1 n R f the instabilities must come from filtering: F s [Is 1 n R f ](θ, σ) = σ n 1 F s R f (θ, σ) = (2π) (n 1)/2 σ n 1 f ˆ(σθ) High frequencies are amplified! 21 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
87 ILL-POSEDNESS OF COMPUTED TOMOGRAPHY Reconstruction is mildly ill-posed (of order (n-1)/2) c f H α 0 R f H α+(n 1)/2 C f H α 0 What causes instabilities in the reconstruction??? We know R and R continuous (stable) operations. Given the inversion formula f = 1/2(2π) 1 n R Is 1 n R f the instabilities must come from filtering: F s [Is 1 n R f ](θ, σ) = σ n 1 F s R f (θ, σ) = (2π) (n 1)/2 σ n 1 f ˆ(σθ) High frequencies are amplified! This causes instabilities since noise is a high frequency phenomenon. 21 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
88 ILL-POSEDNESS OF COMPUTED TOMOGRAPHY Reconstruction is mildly ill-posed (of order (n-1)/2) c f H α 0 R f H α+(n 1)/2 C f H α 0 What causes instabilities in the reconstruction??? We know R and R continuous (stable) operations. Given the inversion formula f = 1/2(2π) 1 n R Is 1 n R f the instabilities must come from filtering: F s [Is 1 n R f ](θ, σ) = σ n 1 F s R f (θ, σ) = (2π) (n 1)/2 σ n 1 f ˆ(σθ) High frequencies are amplified! This causes instabilities since noise is a high frequency phenomenon. Regularizationation by replacing the filtered with a band-limited version: σ n 1 w(σ) σ n 1 21 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
89 REGULARIZED FBP For α > 0, let ω α : ( 1/α, 1/α) [0, ) be smooth such that ω α (σ) σ as α 0 ( σ), and set ψ α (σ) := F 1 (ω α (σ) σ ). Stabilized FBP inversion R α (g)(x) := 1 (g θ ψ α )(x θ) dσ(θ) 4π S 1 Remark f α = R α (g) is a low-pass filtered version of f Plot of ˆψ α = ω α(σ) σ 22 Mathematics of computerized tomography Jürgen Frikel (DTU Compute)
90 See you next week! 23 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 Nov. 11: Introduction to the mathematics of computerized tomography 2 Today: Introduction to the basic concepts
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 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 informationThe Radon transform. Chris Stolk. December 18, 2014
The Radon transform Chris Stolk December 18, 2014 1 Introduction In two dimensions the Radon transform is an integral transform that maps a function to its integrals over lines. Let θ S 1 and s R then
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 informationWeighted Radon transforms for which the Chang approximate inversion formula is precise. R.G. Novikov
ECOLE POLYTECHNIQUE CENTE DE MATHÉMATIQUES APPLIQUÉES UM CNS 764 928 PALAISEAU CEDEX (FANCE). Tel.: 0 69 33 46 00. Fax: 0 69 33 46 46 http : //www.cmap.polytechnique.fr/ Weighted adon transforms for which
More informationError Estimates and Convergence Rates for Filtered Back Projection Reconstructions
Universität Hamburg Dissertation Error Estimates and Convergence ates for Filtered Back Projection econstructions Dissertation zur Erlangung des Doktorgrades an der Fakultät für Mathematik, Informatik
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 informationAn Application of Abel s Method to the Inverse Radon Transform
J. Inverse Ill-Posed Probl.? (????), 1 14 DOI 1.1515/jip-????-??? de Gruyter???? An Application of Abel s Method to the Inverse Radon Transform Shavkat Alimov, Joseph David, Alexander Nolte and Julie Sherman
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 informationAn analog of Chang inversion formula for weighted Radon transforms in multidimensions
An analog of Chang inversion formula for weighted adon transforms in multidimensions Fedor Goncharov, oman Novikov To cite this version: Fedor Goncharov, oman Novikov. An analog of Chang inversion formula
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 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 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 informationInverse problems and medical imaging
Inverse problems and medical imaging Bastian von Harrach harrach@math.uni-stuttgart.de Chair of Optimization and Inverse Problems, University of Stuttgart, Germany Rhein-Main Arbeitskreis Mathematics of
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 informationWeighted Radon transforms for which the Chang approximate inversion formula is precise
Weighted adon transforms for which the Chang approximate inversion formula is precise oman Novikov To cite this version: oman Novikov. Weighted adon transforms for which the Chang approximate inversion
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 informationInverse problems and medical imaging
Inverse problems and medical imaging Bastian von Harrach harrach@math.uni-frankfurt.de Institute of Mathematics, Goethe University Frankfurt, Germany Colloquium of the Department of Mathematics Saarland
More informationKernel-based Image Reconstruction from scattered Radon data by (anisotropic) positive definite functions
Kernel-based Image Reconstruction from scattered Radon data by (anisotropic) positive definite functions Stefano De Marchi 1 February 5, 2016 1 Joint work with A. Iske (Hamburg, D), A. Sironi (Lousanne,
More informationImproved Radon Based Imaging using the Shearlet Transform
Improved Radon Based Imaging using the Shearlet Transform Glenn R. Easley a, Flavia Colonna b, Demetrio Labate c a System Planning Corporation, Arlington, Virginia b George Mason University, Fairfax, Virginia
More informationComputed Tomography Lab phyc30021 Laboratory and Computational Physics 3
Computed Tomography Lab phyc30021 Laboratory and Computational Physics 3 By Joshua P. Ellis School of Physics Last compiled 8th August 2017 Contents 1. Introduction 1 1.1. Aim.............................................
More informationInverse problems in statistics
Inverse problems in statistics Laurent Cavalier (Université Aix-Marseille 1, France) YES, Eurandom, 10 October 2011 p. 1/27 Table of contents YES, Eurandom, 10 October 2011 p. 2/27 Table of contents 1)
More informationInverse problems and medical imaging
Inverse problems and medical imaging Bastian von Harrach harrach@math.uni-frankfurt.de Institute of Mathematics, Goethe University Frankfurt, Germany Seminario di Calcolo delle Variazioni ed Equazioni
More informationTHE RADON TRANSFORM JOSHUA BENJAMIN III. CLASS OF 2020
THE RADON TRANSFORM JOSHUA BENJAMIN III CLASS OF 2020 jbenjamin@college.harvard.edu Acknowledgments I would like to thank: Fabian Haiden for allowing me to create an exposition and for offering advice
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 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 informationKernel-based Image Reconstruction from Scattered Radon Data
Volume 9 2016 Pages 1 12 Kernel-based Image Reconstruction from Scattered Radon Data Stefano De Marchi a Armin Iske b Amos Sironi c Abstract Computerized tomography requires suitable numerical methods
More informationKernel-based Image Reconstruction from Scattered Radon Data
Proceedings of Kernel-based Methods and Function Approximation 2016, Volume 9 2016 Pages 19 31 Kernel-based Image Reconstruction from Scattered Radon Data Stefano De Marchi a Armin Iske b Amos Sironi c
More informationEE 4372 Tomography. Carlos E. Davila, Dept. of Electrical Engineering Southern Methodist University
EE 4372 Tomography Carlos E. Davila, Dept. of Electrical Engineering Southern Methodist University EE 4372, SMU Department of Electrical Engineering 86 Tomography: Background 1-D Fourier Transform: F(
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 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 informationThe Mathematics of CT-Scans
The Mathematics of CT-Scans Tomography has become one of the most important applications of mathematics to the problems of keeping us alive. Modern medicine relies heavily on imaging methods, beginning
More informationInversions of ray transforms on simple surfaces
Inversions of ray transforms on simple surfaces François Monard Department of Mathematics, University of Washington. June 09, 2015 Institut Henri Poincaré - Program on Inverse problems 1 / 42 Outline 1
More informationKey-Words: - fast algorithms, boundary value problems, partial differential equations, Radon transform, MATLAB software
Fast Algorithms and MATLAB Software for Solution of the Dirichlet Boundary Value Problems for Elliptic Partial Differential Equations in Domains with Complicated Geometry ALEXANDRE GREBENNIKOV Faculty
More informationWhat are inverse problems?
What are inverse problems? 6th Forum of Polish Mathematicians Warsaw University of Technology Joonas Ilmavirta University of Jyväskylä 8 September 2015 Joonas Ilmavirta (Jyväskylä) Inverse problems 8 September
More 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 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 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 informationOn single-photon emission computed tomography imaging based on an exact formula for the nonuniform attenuation correction
INSTITUTE OF PHYSICSPUBLISHING Inverse Problems 18 (2002) L11 L19 INVERSE PROBLEMS PII: S0266-5611(02)39995-7 LETTER TO THE EDITOR On single-photon emission computed tomography imaging based on an exact
More informationMB-JASS D Reconstruction. Hannes Hofmann March 2006
MB-JASS 2006 2-D Reconstruction Hannes Hofmann 19 29 March 2006 Outline Projections Radon Transform Fourier-Slice-Theorem Filtered Backprojection Ramp Filter 19 29 Mar 2006 Hannes Hofmann 2 Parallel Projections
More informationTopics. EM spectrum. X-Rays Computed Tomography Direct Inverse and Iterative Inverse Backprojection Projection Theorem Filtered Backprojection
Bioengineering 28A Principles of Biomedical Imaging Fall Quarter 24 X-Rays/CT Lecture Topics X-Rays Computed Tomography Direct Inverse and Iterative Inverse Backprojection Projection Theorem Filtered Backprojection
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 informationChap. 15 Radiation Imaging
Chap. 15 Radiation Imaging 15.1 INTRODUCTION Modern Medical Imaging Devices Incorporating fundamental concepts in physical science and innovations in computer technology Nobel prize (physics) : 1895 Wilhelm
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 informationVisualization of sound field and sound source vibration using laser measurement method
Visualization of sound field and sound source vibration using laser measurement method PACS: 43.58.Fm ABSTRACT Yasuhiro Oikawa (1), Tomomi Hasegawa (1), Yasuhiro Ouchi (2), Yoshio Yamasaki (1) and Yusuke
More information2D X-Ray Tomographic Reconstruction From Few Projections
2D X-Ray Tomographic Reconstruction From Few Projections Application of Compressed Sensing Theory CEA-LID, Thalès, UJF 6 octobre 2009 Introduction Plan 1 Introduction 2 Overview of Compressed Sensing Theory
More informationTopics. EM spectrum. X-Rays Computed Tomography Direct Inverse and Iterative Inverse Backprojection Projection Theorem Filtered Backprojection
Bioengineering 28A Principles of Biomedical Imaging Fall Quarter 25 X-Rays/CT Lecture Topics X-Rays Computed Tomography Direct Inverse and Iterative Inverse Backprojection Projection Theorem Filtered Backprojection
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 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 informationThe spherical mean value operator with centers on a sphere
The spherical mean value operator with centers on a sphere David Finch Department of Mathematics Oregon State University Corvallis, OR 9733-465 Email: finch@math.oregonstate.edu Rakesh Department of Mathematics
More informationBioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2005 X-Rays/CT Lecture 1. Topics
Bioengineering 28A Principles of Biomedical Imaging Fall Quarter 25 X-Rays/CT Lecture Topics X-Rays Computed Tomography Direct Inverse and Iterative Inverse Backprojection Projection Theorem Filtered Backprojection
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 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 informationThe Radon Transform and the Mathematics of Medical Imaging
Colby College Digital Commons @ Colby Honors Theses Student Research 2012 The Radon Transform and the Mathematics of Medical Imaging Jen Beatty Colby College Follow this and additional works at: http://digitalcommons.colby.edu/honorstheses
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 informationMATHEMATICAL SOLUTION OF SOME INDUSTRIAL PROBLEM IN TECHNICAL ELECTRODYNAMICS. Irina DMITRIEVA 1. Abstract
Bulletin of the Transilvania University of Braşov Series III: Mathematics, Informatics, Physics, Vol 5(54) 2012, Special Issue: Proceedings of the Seventh Congress of Romanian Mathematicians, 103-110,
More informationOptimized analytic reconstruction for SPECT
Jean-Pol Guillement, Roman Novikov To cite this version: Jean-Pol Guillement, Roman Novikov. Optimized analytic reconstruction for SPECT. Journal of Inverse ill-posed Problems, De Gruter, 2012, 20 (4),
More informationComputerized Tomography and its Applications: a Guided Tour
Computerized Tomography and its Applications: a Guided Tour J.B.T.M. Roerdink Dept. of Computing Science University of Groningen P.O. Box 800, 9700 AV Groningen The Netherlands Abstract We present a review
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 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 informationInverse problems in statistics
Inverse problems in statistics Laurent Cavalier (Université Aix-Marseille 1, France) Yale, May 2 2011 p. 1/35 Introduction There exist many fields where inverse problems appear Astronomy (Hubble satellite).
More informationIntroduction to the Mathematics of Medical Imaging
Introduction to the Mathematics of Medical Imaging Second Edition Charles L. Epstein University of Pennsylvania Philadelphia, Pennsylvania EiaJTL Society for Industrial and Applied Mathematics Philadelphia
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 informationWave equation techniques for attenuating multiple reflections
Wave equation techniques for attenuating multiple reflections Fons ten Kroode a.tenkroode@shell.com Shell Research, Rijswijk, The Netherlands Wave equation techniques for attenuating multiple reflections
More informationSeismic inverse scattering by reverse time migration
Seismic inverse scattering by reverse time migration Chris Stolk 1, Tim Op t Root 2, Maarten de Hoop 3 1 University of Amsterdam 2 University of Twente 3 Purdue University MSRI Inverse Problems Seminar,
More informationMedical Imaging. Transmission Measurement
Medical Imaging Image Recontruction from Projection Prof Ed X. Wu Tranmiion Meaurement ( meaurable unknown attenuation, aborption cro-ection hadowgram ) I(ϕ,ξ) µ(x,) I = I o exp A(x, )dl L X-ra ource X-Ra
More informationSampling in tomography
in tomography Echantillonnage en tomographie L. Desbat TIMC-IMAG, UMR 5525, UJF-Grenoble 1 15 June 2012 L. Desbat in tomography Outline CT, Radon transform and its inversion 1 CT, Radon transform and its
More informationRecovering a function from its spherical mean values in two and three dimensions
Recovering a function from its spherical mean values in two and three dimensions David Finch Rakesh 1 Introduction The problem of determining a function from a subset of its spherical means has a rich
More informationBackprojection. Projections. Projections " $ & cosθ & Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2014 CT/Fourier Lecture 2
Backprojection Bioengineering 280A Principles of Biomedical Imaging 3 0 0 0 0 3 0 0 0 0 0 3 0 0 3 0 3 Fall Quarter 2014 CT/Fourier Lecture 2 0 0 0 1 1 1 0 0 0 1 0 0 1 2 1 0 0 1 1 1 0 1 3 1 0 1 1 1 1 1
More informationAppliedMathematics. Iterative Methods for Photoacoustic Tomography with Variable Sound Speed. Markus Haltmeier, Linh V. Nguyen
Nr. 29 7. December 26 Leopold-Franzens-Universität Innsbruck Preprint-Series: Department of Mathematics - Applied Mathematics Iterative Methods for Photoacoustic Tomography with Variable Sound Speed Markus
More 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 informationInverse Transport Problems and Applications. II. Optical Tomography and Clear Layers. Guillaume Bal
Inverse Transport Problems and Applications II. Optical Tomography and Clear Layers Guillaume Bal Department of Applied Physics & Applied Mathematics Columbia University http://www.columbia.edu/ gb23 gb23@columbia.edu
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Image econstruction from Scattered adon Data by Weighted Positive Definite Kernel Functions S. De Marchi, A. Iske, G. Santin Nr. 2017-02 January 2017 Image
More informationNon-abelian Radon transform and its applications
Non-abelian Radon transform and its applications Roman Novikov CNRS, Centre de Mathématiques Appliquées, Ecole Polytechnique March 23, 2017 1/25 Consider the equation θ x ψ +A(x,θ)ψ = 0, x R d, θ S d 1,
More informationPolynomial Approximation: The Fourier System
Polynomial Approximation: The Fourier System Charles B. I. Chilaka CASA Seminar 17th October, 2007 Outline 1 Introduction and problem formulation 2 The continuous Fourier expansion 3 The discrete Fourier
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 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 informationFitting functions to data
1 Fitting functions to data 1.1 Exact fitting 1.1.1 Introduction Suppose we have a set of real-number data pairs x i, y i, i = 1, 2,, N. These can be considered to be a set of points in the xy-plane. They
More informationAnalysis and X-ray tomography
Joonas Ilmavirta joonas.ilmavirta@jyu.fi November 2017 arxiv:1711.06557v1 [math.fa] 15 Nov 2017 These are lecture notes for the course MATS4300 Analysis and X-ray tomography given at the University of
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 informationLinear Inverse Problems
Linear Inverse Problems Ajinkya Kadu Utrecht University, The Netherlands February 26, 2018 Outline Introduction Least-squares Reconstruction Methods Examples Summary Introduction 2 What are inverse problems?
More informationFourier Series. (Com S 477/577 Notes) Yan-Bin Jia. Nov 29, 2016
Fourier Series (Com S 477/577 otes) Yan-Bin Jia ov 9, 016 1 Introduction Many functions in nature are periodic, that is, f(x+τ) = f(x), for some fixed τ, which is called the period of f. Though function
More informationApproximation numbers of Sobolev embeddings - Sharp constants and tractability
Approximation numbers of Sobolev embeddings - Sharp constants and tractability Thomas Kühn Universität Leipzig, Germany Workshop Uniform Distribution Theory and Applications Oberwolfach, 29 September -
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 informationNew relationship between the divergent beam projection and the Radon transform
Journal of X-Ray Science and Technology 9 ( 385 4 385 DOI.333/XST--3 IOS Press New relationship between the divergent beam projection and the Radon transform Yuchuan Wei a, and Ge Wang a,b a Biomedical
More informationPrinciples of Computed Tomography (CT)
Page 298 Princiles of Comuted Tomograhy (CT) The theoretical foundation of CT dates back to Johann Radon, a mathematician from Vienna who derived a method in 1907 for rojecting a 2-D object along arallel
More informationarxiv: v2 [math.ap] 27 Jul 2016
THE BROKEN RAY TRANSFORM IN n DIMENSIONS WITH FLAT REFLECTING BOUNDARY MARK HUBENTHAL arxiv:1310.7156v2 [math.ap] 27 Jul 2016 Abstract. We study the broken ray transform on n-dimensional Euclidean domains
More informationGaussian Measure of Sections of convex bodies
Gaussian Measure of Sections of convex bodies A. Zvavitch Department of Mathematics, University of Missouri, Columbia, MO 652, USA Abstract In this paper we study properties of sections of convex bodies
More informationYERNAT M. ASSYLBEKOV AND PLAMEN STEFANOV
SHARP STABILITY ESTIMATE FOR THE GEODESIC RAY TRANSFORM YERNAT M. ASSYLBEKOV AND PLAMEN STEFANOV Abstract. We prove a sharp L 2 H stability estimate for the geodesic X-ray transform of tensor fields of
More informationWhat is a Hilbert C -module? arxiv:math/ v1 [math.oa] 29 Dec 2002
What is a Hilbert C -module? arxiv:math/0212368v1 [math.oa] 29 Dec 2002 Mohammad Sal Moslehian Dept. of maths., Ferdowsi Univ., P.O.Box 1159, Mashhad 91775, Iran E-mail: msalm@math.um.ac.ir abstract.in
More informationOptimal mass transport as a distance measure between images
Optimal mass transport as a distance measure between images Axel Ringh 1 1 Department of Mathematics, KTH Royal Institute of Technology, Stockholm, Sweden. 21 st of June 2018 INRIA, Sophia-Antipolis, France
More information1D Convolution. Convolution g[m] = g[0]δ[m] + g[1]δ[m 1] + g[2]δ[m 2] h[m',k] = L[δ[m k]] = h[ m $ k]
cm 5 cm Bioengineering 28A Principles of Biomeical Imaging Assume =/2 Fall Quarter 24 X-Rays Lecture 2 cm 5 cm cm cm Assume =/2 Assume =/2 Convolution g[m] = g[]δ[m] + g[]δ[m ] + g[2]δ[m 2] h[m,k] = L[δ[m
More informationIterative Correction Methods for the Beam Hardening Artifacts in Computed Tomography. Chelsea Hall Department of Mathematics Oregon State University
Iterative Correction Methods for the Beam Hardening Artifacts in Computed Tomography Chelsea Hall Department of Mathematics Oregon State University December 2013 Abstract In this paper, we discuss three
More informationSuperiorized Inversion of the Radon Transform
Superiorized Inversion of the Radon Transform Gabor T. Herman Graduate Center, City University of New York March 28, 2017 The Radon Transform in 2D For a function f of two real variables, a real number
More informationNonlinear Flows for Displacement Correction and Applications in Tomography
Nonlinear Flows for Displacement Correction and Applications in Tomography Guozhi Dong 1 and Otmar Scherzer 1,2 1 Computational Science Center, University of Vienna, Oskar-Morgenstern-Platz 1, 1090 Wien,
More informationBiomedical Engineering Image Formation
Biomedical Engineering Image Formation PD Dr. Frank G. Zöllner Computer Assisted Clinical Medicine Medical Faculty Mannheim Learning objectives! Understanding the process of image formation! Point spread
More informationInverse Scattering Problems: To Overcome the Ill-posedness
Inverse Scattering Problems: To Overcome the Ill-posedness Gang Bao Zhejiang University IMS, NUS, Sept 25, 2018 1 Acknowledgements Collaborators: Peijun Li, Hai Zhang, Junshan Lin, Faouzi Triki, Tao Yin
More informationarxiv: v1 [math.ca] 23 Oct 2018
A REMARK ON THE ARCSINE DISTRIBUTION AND THE HILBERT TRANSFORM arxiv:80.08v [math.ca] 3 Oct 08 RONALD R. COIFMAN AND STEFAN STEINERBERGER Abstract. We prove that if fx) ) /4 L,) and its Hilbert transform
More informationImage amplitudes in reverse time migration/inversion
Image amplitudes in reverse time migration/inversion Chris Stolk 1, Tim Op t Root 2, Maarten de Hoop 3 1 University of Amsterdam 2 University of Twente 3 Purdue University TRIP Seminar, October 6, 2010
More information