Saturation Rates for Filtered Back Projection
|
|
- Agatha Miles
- 5 years ago
- Views:
Transcription
1 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 2018
2 Outline 1 Filtered Back Projection Basic Reconstruction Problem Approximate Reconstruction Formula 2 Analysis of the Reconstruction Error in the L 2 -Norm Error Estimate Convergence 3 L 2 -Error Analysis for C k -Window Functions Numerical Observations Error Estimate for C k -Windows 4 Asymptotic L 2 -Error Analysis Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 2 / 23
3 Basic Reconstruction Problem Problem formulation: Let Ω R 2 be bounded. Reconstruct a bivariate function f f x, y) with support suppf ) Ω from given Radon data {Rf t, θ) t R, θ [0, π)}, where the Radon transform Rf of f L 1 R 2 ) is defined as Rf t, θ) = f x, y) dx dy for t, θ) R [0, π). {x cosθ)+y sinθ)=t} Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 3 / 23
4 Basic Reconstruction Problem Problem formulation: Let Ω R 2 be bounded. Reconstruct a bivariate function f f x, y) with support suppf ) Ω from given Radon data {Rf t, θ) t R, θ [0, π)}, where the Radon transform Rf of f L 1 R 2 ) is defined as Rf t, θ) = f x, y) dx dy for t, θ) R [0, π). {x cosθ)+y sinθ)=t} y n θ = sinθ), cosθ) ) nθ = cosθ), sinθ) ) xl = t cosθ),tsinθ) ) lt,θ x Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 3 / 23
5 Basic Reconstruction Problem Problem formulation: Let Ω R 2 be bounded. Reconstruct a bivariate function f f x, y) with support suppf ) Ω from given Radon data {Rf t, θ) t R, θ [0, π)}, where the Radon transform Rf of f L 1 R 2 ) is defined as Rf t, θ) = f x, y) dx dy for t, θ) R [0, π). {x cosθ)+y sinθ)=t} Analytical solution: The inversion of R involves the back projection Bh of h L 1 R [0, π)), Bhx, y) = 1 π π 0 hx cosθ) + y sinθ), θ) dθ for x, y) R 2, and is given, for f L 1 R 2 ) CR 2 ), by the filtered back projection formula f x, y) = 1 2 B F 1 [ S FRf )S, θ)] ) x, y) x, y) R 2. Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 3 / 23
6 t The Shepp-Logan Phantom and its Radon Transform The Shepp-Logan phantom f θ The Radon transform Rf Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 4 / 23
7 The Back Projection of the Shepp-Logan Phantom Original Shepp Logan Phantom Unfiltered backprojection Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 5 / 23
8 The Filtered Back Projection Formula OBS! The filtered back projection formula f x, y) = 1 2 B F 1 [ S FRf )S, θ)] ) x, y) is ugly and unstable. Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 6 / 23
9 The Filtered Back Projection Formula OBS! The filtered back projection formula f x, y) = 1 2 B F 1 [ S FRf )S, θ)] ) x, y) is ugly and unstable. Stabilization: Replace S by a low-pass filter A L S) = S W S /L) with finite bandwidth L > 0 and even window W of compact suppw ) [ 1, 1]. Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 6 / 23
10 Approximate Reconstruction Stabilization: Replace the factor S by a low-pass filter A L : R R, A L S) = S W S /L) with finite bandwidth L > 0 and an even window function W : R R with compact support suppw ) [ 1, 1]. Approximate reconstruction formula: We can express the resulting approximate FBP reconstruction f L as f L = 1 2 B F 1 A L Rf ) = f K L, where we rely for f L 1 R 2 ) and g L 1 R [0, π)) on the standard relation Bg f = Bg Rf ) and define the convolution kernel K L : R 2 R as K L x, y) = 1 2 B F 1 A L ) x, y) for x, y) R 2. Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 7 / 23
11 Three Popular Examples for Low-Pass Filters The Ram-Lak filter. { S iff S L; A L S) = S L S) = 0 iff S > L. The Shepp-Logan filter. A L S) = S ) sinπs/2l)) L S) = πs/2l) { 2L π sinπs/2l)) iff S L; 0 iff S > L. The low-pass cosine filter. { S cosπs/2l)) iff S L; A L S) = S cosπs/2l)) L S) = 0 iff S > L. Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 8 / 23
12 Three Popular Examples for Low-Pass Filters 10 Ram Lak filter 7 Shepp Logan filter 4 low pass cosine filter Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 9 / 23
13 Analysis of the Reconstruction Error Aim Analyse the FBP reconstruction error e L = f f L depending on the window function W and the bandwidth L > 0. Previous results: Pointwise and L -error estimates by [Munshi et al., 1991, Munshi, 1992] L p -error estimates in terms of L p -moduli of continuity by [Madych, 1990] Our approach: L 2 -error estimates for target functions f from Sobolev spaces of fractional order, i.e., f H α R 2 ) = { g S R 2 ) g α < } for α > 0, where g 2 α = x 2 + y 2) α Fgx, y) 2 dx dy 2π R R L 2 -convergence rates L ) in terms of bandwidth L and smoothness α Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 10 / 23
14 A first L 2 Error Estimate Theorem A first L 2 -error estimate; Beckmann & I., 2015) Let f L 1 R 2 ) H α R 2 ), for some α > 0, W L R) and K L L 1 R 2 ). Then, the L 2 -norm of the FBP reconstruction error e L = f f L is bounded above by e L L2 R 2 ) 1 W,[ 1,1] f L2 R 2 ) + L α f α. Proof: For f L 1 R 2 ) L 2 R 2 ) and by the Rayleigh Plancherel theorem, we get e L 2 L 2 R 2 ) = f f K L 2 L 2 R 2 ) = 1 2π Ff Ff FK L 2 L 2 R 2 ;C) = 1 2π Ff W L Ff 2 L 2 R 2 ;C) = I 1 + I 2, where we used W L x, y) 2 ) = FK L x, y) for K L L 1 R 2 ), and where we let I 1 := 1 Ff W L Ff )x, y) 2 dx, y), 2π x,y) 2 L I 2 := 1 Ff x, y) 2 dx, y). 2π x,y) 2>L Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 11 / 23
15 A first L 2 Error Estimate Proof continued): For W L R), the first integral I 1 = 1 Ff W L Ff )x, y) 2 dx, y), 2π can be bounded above by x,y) 2 L I 1 1 2π 1 W L 2,[ L,L] Ff 2 L 2 R 2 ;C) = 1 W 2,[ 1,1] f 2 L 2 R 2 ). For f H α R 2 ), with α > 0, the second integral I 2 = 1 Ff x, y) 2 dx, y) 2π x,y) 2>L can be bounded above by I x 2 + y 2) α L 2α Ff x, y) 2 dx, y) L 2α f 2 2π α. x,y) 2>L Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 12 / 23
16 Refined L 2 -Error Analysis Theorem Refined L 2 -error estimate; Beckmann & I., 2016) Let f L 1 R 2 ) H α R 2 ) for some α > 0, let W L R) and K L L 1 R 2 ). Then, the L 2 -norm of the FBP reconstruction error e L = f f L is bounded above by e L L2 R 2 ) Φ 1 /2 α,w L) + L α) f α, where Φ α,w L) = 1 W S)) 2 sup S [ 1,1] 1 + L 2 S 2 ) α for L > 0. Proof: For f H α R 2 ), with α > 0, we bound integral I 2 as before and I 1 by I 1 = 1 2π x,y) 2 L sup S [ L,L] 1 W L x, y) 2 ) x 2 + y 2 ) α 1 + x 2 + y 2) α Ff x, y) 2 dx, y) ) x 2 + y 2) α Ff x, y) 2 dx dy. 2π 1 W L S)) S 2 ) α R R Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 13 / 23
17 Refined L 2 -Error Analysis Proof continued): Since 1 W L S)) 2 1 W sup S [ L,L] 1 + S 2 ) α = sup S /L)) 2 1 W S)) 2 S [ L,L] 1 + S 2 ) α = sup S [ 1,1] 1 + L 2 S 2 ) α, and with letting we find and so e L 2 L 2 R 2 ) I 1 Φ α,w L) = sup S [ 1,1] sup S [ 1,1] 1 W S)) 2 sup S [ 1,1] 1 + L 2 S 2 ) α for L > 0 ) 1 W S)) L 2 S 2 ) α f 2 α = Φ α,w L) f 2 α ) 1 W S)) L 2 S 2 ) α + L 2α f 2 α = Φ α,w L) + L 2α) f 2 α. Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 14 / 23
18 Refined L 2 -Error Analysis Theorem Convergence of Φ α,w ; Beckmann & I., 2016) Let W be continuous on [ 1, 1] and satisfy W 0) = 1. Then, for all α > 0, 1 W S)) 2 Φ α,w L) = max S [0,1] 1 + L 2 S 2 ) α 0 for L. Proof sketch): Let Sα,W,L [0, 1] be the smallest maximizer in [0, 1]of Φ α,w,l S) := 1 W S))2 1 + L 2 S 2 ) α for S [0, 1]. Case 1: S α,w,l uniformly bounded away from 0, S α,w,l c c α,w > 0. Then, 0 Φ α,w,l S α,w,l ) = 1 W S α,w,l ) ) L2 S α,w,l )2) α 1 W 2,[ 1,1] 1 + L 2 c 2 ) α L 0. Case 2: Sα,W,L 0 for L. Then, 0 Φ α,w,l S α,w,l ) = 1 W S α,w,l ) ) L2 S α,w,l )2) α 1 W S α,w,l) ) 2 L 0. Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 15 / 23
19 Convergence of the FBP Reconstruction Corollary L 2 convergence of FBP reconstruction; Beckmann & I., 16) Let f L 1 R 2 ) H α R 2 ) for some α > 0, K L L 1 R 2 ) and W C[ 1, 1]) with W 0) = 1. Then, the L 2 -norm of the FBP reconstruction error e L = f f L satisfies e L L2 R 2 ) = o1) for L. Basic Assumption Let S α,w,l be uniformly bounded away from 0, i.e., there is c α,w > 0 satisfying S α,w,l c α,w L > 0. OBS! Under the basic assumption we have Φ α,w L) c 2α α,w 1 W 2,[ 1,1] L 2α = OL 2α ) for L. Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 16 / 23
20 Rate of Convergence Theorem Rate of convergence; Beckmann & I., 2016) Let f L 1 R 2 ) H α R 2 ) for some α > 0, let K L L 1 R 2 ) and W C[ 1, 1]) with W 0) = 1. Further, let the basic assumption be satisfied for a constant c α,w. Then, the L 2 -norm of the FBP reconstruction error e L = f f L is bounded by ) e L L2 R 2 ) c α α,w 1 W,[ 1,1] + 1 L α f α. Therefore, e L L2 R 2 ) = OL α ) for L, i.e., the decay rate is determined by the smoothness α of the target function f. OBS! Let the window function W C[ 1, 1]), W χ [ 1,1], satisfy W S) = 1 S ε, ε) for some 0 < ε < 1. Then, the basic assumption is fulfilled with c α,w = ε. Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 17 / 23
21 Numerical Observations We investigated the behaviour of Sα,W,L numerically for standard low-pass filters: Shepp-Logan filter: W S) = sinc ) πs 2 χ[ 1,1] S), Cosine filter: W S) = cos ) πs 2 χ[ 1,1] S), Hamming filter for β [ 1 2, 1] ): W S) = β + 1 β) cosπs)) χ [ 1,1] S), Gaussian filter for β > 1): W S) = exp πs /β) 2) χ [ 1,1] S). For α < 2, we found that the basic assumption is fulfilled whereby c α,w > 0 L > 0 : S α,w,l c α,w Φ α,w L) = OL 2α ) for L. But for α 2, we observed Sα,W,L 0 for L. Moreover, in this case the convergence rate of Φ α,w stagnates at Φ α,w L) = OL 4 ) for L. Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 18 / 23
22 Numerical Observations Φα,W L 2α Φα,W Φα,W Φα,W L a) α = L b) α = L c) α = 2 Φα,W L Φα,W Φα,W Φα,W L d) α = L e) α = L f) α = 4 Fig.: Decay rate of Φ α,w for the Shepp-Logan filter Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 19 / 23
23 L 2 -Error Analysis for C k -Windows Theorem Convergence rate of Φ α,w for W C k ; Beckmann & I., 16) For k 2, let W C k [ 1, 1] satisfy Then, we have where the constant W 0) = 1, W j) 0) = 0 1 j k 1 Φ α,w L) { c 2 α,k for α > k k!) W k) 2 2,[ 1,1] L 2k 1 k!) W k) 2 2,[ 1,1] L 2α for α k, c α,k = k α k ) k/2 α k ) α/2 α is strictly monotonically decreasing in α > k. In particular, we have Φ α,w L) = O L 2 min{k,α}) for L. Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 20 / 23
24 L 2 -Error Analysis for C k -Windows Corollary L 2 -error estimate for C k -windows; Beckmann & I., 2016) Let f L 1 R 2 ) H α R 2 ), for some α > 0, let K L L 1 R 2 ) and W C k [ 1, 1]), for k 2, with W 0) = 1, W j) 0) = 0 1 j k 1. Then, the L 2 -norm of the FBP reconstruction error e L = f f L is bounded by cα,k k! W k),[ 1,1] L k + L ) f α α for α > k e L L 2 R 2 ) 1 k! W k),[ 1,1] L α + L ) f α α for α k. In particular, we have e L L 2 R 2 ) c W k),[ 1,1] L min{k,α} + L α) f α = O L min{k,α}). In conclusion, the rate of convergence saturates at OL k ). Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 21 / 23
25 Asymptotic L 2 -Error Analysis Theorem Asymptotic L 2 -error estimate; Beckmann & I., 2016) Let f L 1 R 2 ) H α R 2 ), for some α > 0, K L L 1 R 2 ) and let W L R) be k-times differentiable at the origin, k 2, with W 0) = 1, W j) 0) = 0 1 j k 1. Then, the L 2 -norm of the FBP reconstruction error e L = f f L is bounded by 2 k! c α,k W k) 0) L k + L ) f α α + ol k ) for α > k e L L2 R 2 ) 2 k! W k) 0) L α + L ) f α α + ol α ) for α k with the strictly monotonically decreasing constant k ) k/2 α k ) α/2 c α,k = for α > k. α k α In particular, we have e L L2 R 2 ) c W k) 0) L min{k,α} + L α) f α + o L min{k,α}). Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 22 / 23
26 References M. Beckmann, A. Iske: Error estimates and convergence rates for filtered back projection, 33 pages, to appear in Mathematics of Computation. M. Beckmann, A. Iske: Sobolev error estimates for filtered back projection reconstructions. IEEE 2017 International Conference Sampling Theory and Applications SampTA2017), 2017, M. Beckmann, A. Iske: Approximation of bivariate functions from fractional Sobolev spaces by filtered back projection. Hamburger Beiträge zur Angewandten Mathematik , 2017, M. Beckmann, A. Iske: Analysis of the inherent reconstruction error in filtered back projection. Hamburger Beiträge zur Angewandten Mathematik , 2016, M. Beckmann, A. Iske: Error estimates for filtered back projection. IEEE 2015 International Conference Sampling Theory and Applications SampTA2015), 2015, W.R. Madych: Summability and approximate reconstruction from Radon transform data. Contemporary Mathematics, Volume 113, AMS, Providence, 1990, P. Munshi: Error analysis of tomographic filters. I: theory. NDT & E International 254/5), 1992, P. Munshi, R.K.S. Rathore, K.S. Ram, M.S. Kalra: Error estimates for tomographic inversion. Inverse Problems 73), 1991, F. Natterer: The Mathematics of Computerized Tomography. Classics in Applied Mathematics 32, SIAM, Philadelphia, Matthias Beckmann University of Hamburg) Saturation Rates for Filtered Back Projection 23 / 23
Hamburger 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 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 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 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 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 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 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 informationOn Ridge Functions. Allan Pinkus. September 23, Technion. Allan Pinkus (Technion) Ridge Function September 23, / 27
On Ridge Functions Allan Pinkus Technion September 23, 2013 Allan Pinkus (Technion) Ridge Function September 23, 2013 1 / 27 Foreword In this lecture we will survey a few problems and properties associated
More informationReproducing Kernel Hilbert Spaces
Reproducing Kernel Hilbert Spaces Lorenzo Rosasco 9.520 Class 03 February 9, 2011 About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
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 information1.5 Approximate Identities
38 1 The Fourier Transform on L 1 (R) which are dense subspaces of L p (R). On these domains, P : D P L p (R) and M : D M L p (R). Show, however, that P and M are unbounded even when restricted to these
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 informationNIH Public Access Author Manuscript Proc Soc Photo Opt Instrum Eng. Author manuscript; available in PMC 2013 December 17.
NIH Public Access Author Manuscript Published in final edited form as: Proc Soc Photo Opt Instrum Eng. 2008 March 18; 6913:. doi:10.1117/12.769604. Tomographic Reconstruction of Band-limited Hermite Expansions
More informationLECTURES 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 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 informationA pointwise limit theorem for filtered backprojection in computed tomography
A pointwise it theorem for filtered backprojection in computed tomography Yangbo Ye a) Department of Mathematics University of Iowa Iowa City Iowa 52242 Jiehua Zhu b) Department of Mathematics University
More informationConvolution/Modulation Theorem. 2D Convolution/Multiplication. Application of Convolution Thm. Bioengineering 280A Principles of Biomedical Imaging
Bioengineering 8A Principles of Biomedical Imaging Fall Quarter 15 CT/Fourier Lecture 4 Convolution/Modulation Theorem [ ] e j πx F{ g(x h(x } = g(u h(x udu = g(u h(x u e j πx = g(uh( e j πu du = H( dxdu
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 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 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 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 informationApproximation by Conditionally Positive Definite Functions with Finitely Many Centers
Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal
More informationDesign of Tomography System for Industrial Plasmas
Design of Tomography System for Industrial Plasmas More info about this article: http://www.ndt.net/?id=22308 1 Kavita RATHORE, 2 Sudeep BHATTACHARJEE, 1 Prabhat MUNSHI 1 Nuclear Engineering and Technology
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 informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
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 informationReproducing Kernel Hilbert Spaces
Reproducing Kernel Hilbert Spaces Lorenzo Rosasco 9.520 Class 03 February 11, 2009 About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
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 informationHARMONIC ANALYSIS. Date:
HARMONIC ANALYSIS Contents. Introduction 2. Hardy-Littlewood maximal function 3. Approximation by convolution 4. Muckenhaupt weights 4.. Calderón-Zygmund decomposition 5. Fourier transform 6. BMO (bounded
More informationLECTURES 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 informationPointwise convergence rates and central limit theorems for kernel density estimators in linear processes
Pointwise convergence rates and central limit theorems for kernel density estimators in linear processes Anton Schick Binghamton University Wolfgang Wefelmeyer Universität zu Köln Abstract Convergence
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 informationReproducing Kernel Hilbert Spaces
Reproducing Kernel Hilbert Spaces Lorenzo Rosasco 9.520 Class 03 February 12, 2007 About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
More informationRemarks on Extremization Problems Related To Young s Inequality
Remarks on Extremization Problems Related To Young s Inequality Michael Christ University of California, Berkeley University of Wisconsin May 18, 2016 Part 1: Introduction Young s convolution inequality
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 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 informationStability of Kernel Based Interpolation
Stability of Kernel Based Interpolation Stefano De Marchi Department of Computer Science, University of Verona (Italy) Robert Schaback Institut für Numerische und Angewandte Mathematik, University of Göttingen
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 informationAP Calculus (BC) Chapter 10 Test No Calculator Section. Name: Date: Period:
AP Calculus (BC) Chapter 10 Test No Calculator Section Name: Date: Period: Part I. Multiple-Choice Questions (5 points each; please circle the correct answer.) 1. The graph in the xy-plane represented
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 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 informationFFTs in Graphics and Vision. Homogenous Polynomials and Irreducible Representations
FFTs in Graphics and Vision Homogenous Polynomials and Irreducible Representations 1 Outline The 2π Term in Assignment 1 Homogenous Polynomials Representations of Functions on the Unit-Circle Sub-Representations
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 informationMATH MEASURE THEORY AND FOURIER ANALYSIS. Contents
MATH 3969 - MEASURE THEORY AND FOURIER ANALYSIS ANDREW TULLOCH Contents 1. Measure Theory 2 1.1. Properties of Measures 3 1.2. Constructing σ-algebras and measures 3 1.3. Properties of the Lebesgue measure
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 informationRKHS, Mercer s theorem, Unbounded domains, Frames and Wavelets Class 22, 2004 Tomaso Poggio and Sayan Mukherjee
RKHS, Mercer s theorem, Unbounded domains, Frames and Wavelets 9.520 Class 22, 2004 Tomaso Poggio and Sayan Mukherjee About this class Goal To introduce an alternate perspective of RKHS via integral operators
More informationConvolution Theorem Modulation Transfer Function y(t)=g(t)*h(t)
MTF = Fourier Transform of PSF Bioengineering 28A Principles of Biomedical Imaging Fall Quarter 214 CT/Fourier Lecture 4 Bushberg et al 21 Convolution Theorem Modulation Transfer Function y(t=g(t*h(t g(t
More informationStability constants for kernel-based interpolation processes
Dipartimento di Informatica Università degli Studi di Verona Rapporto di ricerca Research report 59 Stability constants for kernel-based interpolation processes Stefano De Marchi Robert Schaback Dipartimento
More informationMultiscale Frame-based Kernels for Image Registration
Multiscale Frame-based Kernels for Image Registration Ming Zhen, Tan National University of Singapore 22 July, 16 Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image
More informationRegularity and compactness for the DiPerna Lions flow
Regularity and compactness for the DiPerna Lions flow Gianluca Crippa 1 and Camillo De Lellis 2 1 Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy. g.crippa@sns.it 2 Institut für Mathematik,
More informationBounded uniformly continuous functions
Bounded uniformly continuous functions Objectives. To study the basic properties of the C -algebra of the bounded uniformly continuous functions on some metric space. Requirements. Basic concepts of analysis:
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 informationOutline of Fourier Series: Math 201B
Outline of Fourier Series: Math 201B February 24, 2011 1 Functions and convolutions 1.1 Periodic functions Periodic functions. Let = R/(2πZ) denote the circle, or onedimensional torus. A function f : C
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 2: Radial Basis Function Interpolation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590
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 informationImage Reconstruction from Scattered Radon Data by Weighted Positive Definite Kernel Functions
Image econstruction from Scattered adon Data by Weighted Positive Definite Kernel Functions S. De Marchi A. Iske G. Santin November 3, 2017 Abstract We propose a novel kernel-based method for image reconstruction
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 informationNIH Public Access Author Manuscript Inverse Probl. Author manuscript; available in PMC 2010 April 20.
NIH Public Access Author Manuscript Published in final edited form as: Inverse Probl. 2010 January 1; 26(3): 350131 3501329. doi:10.1088/0266-5611/26/3/035013. High Order Total Variation Minimization for
More informationWhere is the bark, the tree and the forest, mathematical foundations of multi-scale analysis
Where is the bark, the tree and the forest, mathematical foundations of multi-scale analysis What we observe is not nature itself, but nature exposed to our method of questioning. -Werner Heisenberg M.
More informationMath 321 Final Examination April 1995 Notation used in this exam: N. (1) S N (f,x) = f(t)e int dt e inx.
Math 321 Final Examination April 1995 Notation used in this exam: N 1 π (1) S N (f,x) = f(t)e int dt e inx. 2π n= N π (2) C(X, R) is the space of bounded real-valued functions on the metric space X, equipped
More informationTowards Improved Sensitivity in Feature Extraction from Signals: one and two dimensional
Towards Improved Sensitivity in Feature Extraction from Signals: one and two dimensional Rosemary Renaut, Hongbin Guo, Jodi Mead and Wolfgang Stefan Supported by Arizona Alzheimer s Research Center and
More informationProblem Set 5: Solutions Math 201A: Fall 2016
Problem Set 5: s Math 21A: Fall 216 Problem 1. Define f : [1, ) [1, ) by f(x) = x + 1/x. Show that f(x) f(y) < x y for all x, y [1, ) with x y, but f has no fixed point. Why doesn t this example contradict
More information1.3.1 Definition and Basic Properties of Convolution
1.3 Convolution 15 1.3 Convolution Since L 1 (R) is a Banach space, we know that it has many useful properties. In particular the operations of addition and scalar multiplication are continuous. However,
More information1. Introduction. 2. Outlines
1. Introduction Graphs are beneficial because they summarize and display information in a manner that is easy for most people to comprehend. Graphs are used in many academic disciplines, including math,
More informationApproximating Infinity-Dimensional Stochastic Darcy s Equations without Uniform Ellipticity
Approximating Infinity-Dimensional Stochastic Darcy s Equations without Uniform Ellipticity Marcus Sarkis Jointly work with Juan Galvis USC09 Marcus Sarkis (WPI) SPDE without Uniform Ellipticity USC09
More informationRandom graphs: Random geometric graphs
Random graphs: Random geometric graphs Mathew Penrose (University of Bath) Oberwolfach workshop Stochastic analysis for Poisson point processes February 2013 athew Penrose (Bath), Oberwolfach February
More informationGenerated Covariates in Nonparametric Estimation: A Short Review.
Generated Covariates in Nonparametric Estimation: A Short Review. Enno Mammen, Christoph Rothe, and Melanie Schienle Abstract In many applications, covariates are not observed but have to be estimated
More informationA. Iosevich and I. Laba January 9, Introduction
K-DISTANCE SETS, FALCONER CONJECTURE, AND DISCRETE ANALOGS A. Iosevich and I. Laba January 9, 004 Abstract. In this paper we prove a series of results on the size of distance sets corresponding to sets
More information3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due ). Show that the open disk x 2 + y 2 < 1 is a countable union of planar elementary sets. Show that the closed disk x 2 + y 2 1 is a countable
More informationADAPTIVE ALGEBRAIC RECONSTRUCTION OF CT IMAGES IN COMPARISON TO FREQUENCY DOMAIN METHODS. Bernhard Bundschuh Universität-GH-Siegen, Germany
ADAPIVE ALGEBRAIC RECONSRUCION OF C IMAGES IN COMPARISON O FREQUENCY DOMAIN MEHODS Bernhard Bundschuh Universität-GH-Siegen German Abstract his paper presents a new algebraic algorithm for the reconstruction
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 informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
More informationSECTION A. f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes.
SECTION A 1. State the maximal domain and range of the function f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes. 2. By evaluating f(0),
More informationReminder Notes for the Course on Measures on Topological Spaces
Reminder Notes for the Course on Measures on Topological Spaces T. C. Dorlas Dublin Institute for Advanced Studies School of Theoretical Physics 10 Burlington Road, Dublin 4, Ireland. Email: dorlas@stp.dias.ie
More informationSemigroup Growth Bounds
Semigroup Growth Bounds First Meeting on Asymptotics of Operator Semigroups E.B. Davies King s College London Oxford, September 2009 E.B. Davies (KCL) Semigroup Growth Bounds Oxford, September 2009 1 /
More informationConfidence regions for images observed under the Radon transform
Confidence regions for images observed under the Radon transform Nicolai Bissantz, Hajo Holzmann and Katharina Proksch Fakultät für Mathematik Ruhr-Universität Bochum, Germany Fachbereich Mathematik und
More informationStrengthened Sobolev inequalities for a random subspace of functions
Strengthened Sobolev inequalities for a random subspace of functions Rachel Ward University of Texas at Austin April 2013 2 Discrete Sobolev inequalities Proposition (Sobolev inequality for discrete images)
More informationA Banach Gelfand Triple Framework for Regularization and App
A Banach Gelfand Triple Framework for Regularization and Hans G. Feichtinger 1 hans.feichtinger@univie.ac.at December 5, 2008 1 Work partially supported by EUCETIFA and MOHAWI Hans G. Feichtinger hans.feichtinger@univie.ac.at
More informationRepresentations of Gaussian measures that are equivalent to Wiener measure
Representations of Gaussian measures that are equivalent to Wiener measure Patrick Cheridito Departement für Mathematik, ETHZ, 89 Zürich, Switzerland. E-mail: dito@math.ethz.ch Summary. We summarize results
More informationConfidence regions for images observed under the Radon transform
Confidence regions for images observed under the adon transform Nicolai Bissantz, Hajo Holzmann and Katharina Proksch Fakultät für Mathematik uhr-universität Bochum, Germany Fachbereich Mathematik und
More informationSobolev spaces. May 18
Sobolev spaces May 18 2015 1 Weak derivatives The purpose of these notes is to give a very basic introduction to Sobolev spaces. More extensive treatments can e.g. be found in the classical references
More informationNonlinear error dynamics for cycled data assimilation methods
Nonlinear error dynamics for cycled data assimilation methods A J F Moodey 1, A S Lawless 1,2, P J van Leeuwen 2, R W E Potthast 1,3 1 Department of Mathematics and Statistics, University of Reading, UK.
More informationContinuous dependence estimates for the ergodic problem with an application to homogenization
Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,
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 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 informationInterpolation via weighted l 1 -minimization
Interpolation via weighted l 1 -minimization Holger Rauhut RWTH Aachen University Lehrstuhl C für Mathematik (Analysis) Matheon Workshop Compressive Sensing and Its Applications TU Berlin, December 11,
More informationMathematical modelling of collective behavior
Mathematical modelling of collective behavior Young-Pil Choi Fakultät für Mathematik Technische Universität München This talk is based on joint works with José A. Carrillo, Maxime Hauray, and Samir Salem
More informationFolland: Real Analysis, Chapter 8 Sébastien Picard
Folland: Real Analysis, Chapter 8 Sébastien Picard Problem 8.3 Let η(t) = e /t for t >, η(t) = for t. a. For k N and t >, η (k) (t) = P k (/t)e /t where P k is a polynomial of degree 2k. b. η (k) () exists
More informationReducing subspaces. Rowan Killip 1 and Christian Remling 2 January 16, (to appear in J. Funct. Anal.)
Reducing subspaces Rowan Killip 1 and Christian Remling 2 January 16, 2001 (to appear in J. Funct. Anal.) 1. University of Pennsylvania, 209 South 33rd Street, Philadelphia PA 19104-6395, USA. On leave
More informationConvergence rates in weighted L 1 spaces of kernel density estimators for linear processes
Alea 4, 117 129 (2008) Convergence rates in weighted L 1 spaces of kernel density estimators for linear processes Anton Schick and Wolfgang Wefelmeyer Anton Schick, Department of Mathematical Sciences,
More informationANALYSIS QUALIFYING EXAM FALL 2017: SOLUTIONS. 1 cos(nx) lim. n 2 x 2. g n (x) = 1 cos(nx) n 2 x 2. x 2.
ANALYSIS QUALIFYING EXAM FALL 27: SOLUTIONS Problem. Determine, with justification, the it cos(nx) n 2 x 2 dx. Solution. For an integer n >, define g n : (, ) R by Also define g : (, ) R by g(x) = g n
More informationON Y. NIEVERGELT S INVERSION FORMULA FOR THE RADON TRANSFORM. E. Ournycheva, B. Rubin. Abstract
ON Y. NIEVERGELT S INVERSION FORMULA FOR THE RADON TRANSFORM E. Ournycheva, B. Rubin Abstract In 986 Y. Nievergelt suggested a simple formula which allows to reconstruct a continuous compactly supported
More informationn 1 f = f m, m=0 n 1 k=0 Note that if n = 2, we essentially get example (i) (but for complex functions).
Chapter 4 Fourier Analysis 4.0.1 Intuition This discussion is borrowed from [T. Tao, Fourier Transform]. Fourier transform/series can be viewed as a way to decompose a function from some given space V
More information1 Continuity Classes C m (Ω)
0.1 Norms 0.1 Norms A norm on a linear space X is a function : X R with the properties: Positive Definite: x 0 x X (nonnegative) x = 0 x = 0 (strictly positive) λx = λ x x X, λ C(homogeneous) x + y x +
More informationFUNCTIONAL ANALYSIS LECTURE NOTES: WEAK AND WEAK* CONVERGENCE
FUNCTIONAL ANALYSIS LECTURE NOTES: WEAK AND WEAK* CONVERGENCE CHRISTOPHER HEIL 1. Weak and Weak* Convergence of Vectors Definition 1.1. Let X be a normed linear space, and let x n, x X. a. We say that
More informationMAS331: Metric Spaces Problems on Chapter 1
MAS331: Metric Spaces Problems on Chapter 1 1. In R 3, find d 1 ((3, 1, 4), (2, 7, 1)), d 2 ((3, 1, 4), (2, 7, 1)) and d ((3, 1, 4), (2, 7, 1)). 2. In R 4, show that d 1 ((4, 4, 4, 6), (0, 0, 0, 0)) =
More informationInterpolation by Harmonic Polynomials. on Radon Projections
Based on Radon Projections I Georgieva, C Hofreither joint work with C Koutschan, V Pillwein, T Thanatipanonda Sozopol September, 2, 2012 Collaborators Irina Georgieva Institute of Mathematics and Informatics,
More informationA RECONSTRUCTION FORMULA FOR BAND LIMITED FUNCTIONS IN L 2 (R d )
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 127, Number 12, Pages 3593 3600 S 0002-9939(99)04938-2 Article electronically published on May 6, 1999 A RECONSTRUCTION FORMULA FOR AND LIMITED FUNCTIONS
More informationAN IMPLEMENTATION OF THE BACK-PROJECTION ALGORITHM ACCORDING TO SANTOSA AND VOGELIUS
AN IMPLEMENTATION OF THE BACK-PROJECTION ALGORITHM ACCORDING TO SANTOSA AND VOGELIUS Pai Chi Nan Department of Mechanical Engineering - Polytechnic School of University of São Paulo Av. Prof. Mello de
More information