A pointwise limit theorem for filtered backprojection in computed tomography
|
|
- Marshall Stafford
- 5 years ago
- Views:
Transcription
1 A pointwise it theorem for filtered backprojection in computed tomography Yangbo Ye a) Department of Mathematics University of Iowa Iowa City Iowa Jiehua Zhu b) Department of Mathematics University of Iowa Iowa City Iowa Ge ang c) Departments of Mathematics Biomedical Engineering and Radiology University of Iowa Iowa City Iowa Received 6 September 22; accepted for publication 4 January 23; published 22 April 23 Computed tomography CT is one of the most important areas in the modern science and technology. The most popular approach for image reconstruction is filtered backprojection. It is essential to understand the it behavior of the filtered backprojection algorithms. The classic results on the it of image reconstruction are typically done in the norm sense. In this paper we use the method of ited bandwidth to handle filtered backprojection-based image reconstruction when the spectrum of an underlying image is not absolutely integrable. Our main contribution is assuming the method of ited bandwidth to prove a pointwise it theorem for a class of functions practically relevant and quite general. Further work is underway to extend the theory and explore its practical applications. 23 American Association of Physicists in Medicine. DOI:.8/ Key words: computed tomography CT image reconstruction filtered backprojection error analysis ited bandwidth pointwise it I. INTRODUCTION Computed tomography CT is one of the most important areas in modern science and technology. The most popular approach for image reconstruction is filtered backprojection. 2 Let f (xy) be an absolutely integrable function on R 2 representing an object with a finite support such as defined on a unit disk. In this context the standard procedure for computed tomogrphy CT refers to the following steps: Step. Projection: P t f xyx y sin tdx dy Step 2. Filtration: Q t S wwe 2itw dw where S w Step 3. Backprojection: P te 2iwt dt f xy Q x y sin )d. The correctness of the above procedure is usually shown by citing the Fourier slice theorem that 2 3 S w fˆw w sin f xye 2ixw yw sin dx dy. Rigorously speaking the filtered backprojection algorithm defined by 2 and 3 is valid at almost every point (xy) for an absolutely integrable function f (x y) when its Fourier transform fˆ(u) is also absolutely integrable. hen fˆ(u) is not absolutely integrable it is essential to understand the it behavior of the filtered backprojection algorithms. There are various strategies to handle this divergence problem. Examples include the method of ited bandwidth the Abel method of summability and the Shepp Logan filter technique. 8 Generally speaking different strategies have their own advantages and disadvantages. As a matter of fact CT scanner manufacturers typically implement a number of reconstruction filters to balance spatial and contrast resolution for different applications. Typically the method of ited bandwidth produces higher spatial resolution while the Abel method of summability and the Shepp Logan filtering lead to better contrast resolution. The convolution kernel associated with the method of ited bandwidth is usually referred to as the standard kernel while the smooth kernels associated with the Shepp Logan method and other similar methods may be regarded as contrast kernels. The Abel method of summability is to introduce a convergence factor e w and then take : 4 86 Med. Phys. 3 5 May Õ23Õ3 5 Õ86Õ7Õ$2. 23 Am. Assoc. Phys. Med. 86
2 87 Ye Zhu and ang: A it theorem in CT 87 Q t S wwe w e 2itw dw f xy Q x y sin d for and f xy f xy. 7 For any absolutely integrable function f (xy) it is a classical result that f (xy) tends to f (xy) in the L norm as : f xy f xydx dy. Using the concept of the Lebesque set of f (xy) we can further prove that f (xy) tends to f (xy) as almost everywhere. 3 It is our goal in the present paper to study the it behavior of the method of ited bandwidth. This method is to assume that the bandwidth is ited and replace 2 and 3 by 9 and : Q t S wwe 2itw dw f xy Q x y sin d. Then we take the it to reconstruct the image f (xy): f xy f xy. As far as the method of ited bandwidth defined in 9 is concerned a corresponding L statement does not exist. In fact 9 can be written as Q t S ww/we 2itw dw 2 where is the characteristic function of the interval which can be interpreted as the characteristic function in polar coordinates on the unit disk x 2 y 2. Consequently f (xy) is the Fourier inversion of fˆ(uv)(u/v/). By the Fourier convolution formula f xy f xy * u/v/ where ((u/v/)) is expressed as follows: u/v/ xy Here in polar coordinates u/v/ e 2ixuyv du dv 2 ˆ xy. ˆ w w sin J 2w/w 3 4 is not absolutely integrable on the plane where J is the Bessel function of the first kind of order one see 2 below. Therefore f (xy) in3 is not of L in general. e need an L 2 theory to describe the it behavior of the ited bandwidth based reconstruction in. Assume f L (R 2 )L 2 (R 2 ) then fˆl 2 (R 2 ). The Fourier transform of an L 2 function g is defined as the L 2 it of ĝ xy guvu/v/ e 2iuxvy du dv. 5 By the Plancherel theorem 3 we immediately have f xy f xy 2 dx dy. 6 Similar to 6 the classic results on the it of image reconstruction are typically done in the norm sense. The most famous results are summarized in the book by Natterer. 2 Specifically they worked with C and assumed the sup norm defined on H ( n ) as the worst-case error. The only prior information they made use of is f H ( n ) with being about /2 and moderate see pp in Ref. 2 for details. In this paper we study the method of ited bandwidth to handle filtered backprojection based image reconstruction when the spectrum of an underlying image is not absolutely integrable. Particularly assuming the method of a ited bandwidth we prove a pointwise it theorem for a class of functions practically relevant and quite general. In the second section we present our main result the pointwise it theorem. In the third section we give two analytic examples while in the last section we discuss several issues and conclude the paper. II. POINTISE LIMIT THEOREM Because the pointwise it is stronger and much more specific than the L 2 it in this section we prove a pointwise it theorem to refine 6 for a class of quite general functions f (xy). Indeed we will consider functions f (xy) that are smooth but have finitely many jumps of finite magnitude along piecewise continuous curves. Theorem pointwise it theorem : Let f(xy) be an absolutely integrable function on R 2. Assume that other than finitely many jumps of finite magnitude of the function and/or its first and second partial derivatives along piecewise continuous curves this function has absolutely integrable second partial derivatives. Then the filtered backprojection using the method of ited bandwidth converges pointwise at any point (x y ) not on a jump i.e. f x y f x y. e remark that in the scenario of medical imaging there are distinctive components in the human body including solid organs soft tissues bones and so on. ithin each compo- Medical Physics Vol. 3 No. 5 May 23
3 88 Ye Zhu and ang: A it theorem in CT 88 nent the x-ray linear attenuation coefficient varies smoothly and can be well modeled by smooth functions with continuous second derivatives. On the other hand across the boundaries among the components jumps in the x-ray linear attenuation coefficient its first and second partial derivatives can be represented by a finite number of gate-like functions. The discontinuities across the boundaries in the first and second derivatives can be relevant in contrast dynamics such as bolus-driven perfusion for cancer studies. Before we prove Theorem let us formulate this class of functions. As we will use polar coordinates centered at (x y ) in our proof we express the curves of discontinuities in the following way. Note that this expression of our function f (xy) directly depends on a given point (x y ) not on a jump. For a different point (x y ) the polar coordinates will be different and the formulas for the function f (x r y r sin ) will also be different. ithout a loss of generality let rr j () j...m be piecewise continuous curves in. e may further assume that r ()r 2 () r m ()C for any. To capture various jumps of the function and its derivatives on these curves we use a polynomial of r in the region between two adjacent curves. More precisely define jump functions: j h j ra a j r a 5 j r 5 if r j rr j otherwise 7 for j...m where as coefficients a (k) j () are bounded a (k) j ()C 2 piecewise continuous functions of. Here we use a polynomial of degree 5 because we need to use six coefficients a (k) j () to match the function values its first and second derivatives in the direction of r on the outer and inner rims of the region between r j () and r j (). e remark that our jump function in 7 generalizes the gate function s j r a j if r j rr j otherwise. This generalization allows us to control jumps not only in value magnitude but also in derivatives. For our function f (xy) in Theorem we request that in polar coordinates centered at the point (x y ) it can be written as f x r y r sin gx r y r sin h r h m r 8 such that i g(xy) is absolutely integrable on R 2 ii g(xy) has all second partial derivatives everywhere in R 2 and iii its second partial derivatives are also absolutely integrable on R 2. Note that on the right side of 8 the second differentiability of g(x y) on jump curves is achieved by suitably choosing the coefficients a j (k) () of the polynomials h j (r). Actually we may replace conditions ii and iii by a weaker condition that 2gx r y r sin tdot 9 for some fixed ast. Here we used the Fourier ine transform ĝ t grrtdr. The assumption 9 essentially means that the Fourier ine transform of g(x y) along any half-line starting at (x y ) goes to in O(t ) as t. This is a version of the Riemann Lebesque theorem. 4 If as in Theorem g(x y) has absolutely integrable second partial derivatives the bound in 9 will be O(t 2 ). The magnitude of the jumps cannot be allowed to go to infinity as controlled by a j (k) ()C 2 for j...m and k...5. The jumps cannot occur at the point (x y ) either so that r j () cannot approach to. If we denote by cc(x y ) the shortest distance from (x y ) to the nearest discontinuous curve we then have r j ()c(x y ). e emphasize that the class of the above-defined functions f (xy) in8 should be regarded as a universal model for a variety of practical applications. Technically g(x r y r sin ) is smooth and describes continuous changes in the image domain while terms h j (r) capture the structural discontinuities. Combining the two kinds of bases we can satisfactorily approximate images that are regionwise smooth. In the proof below and in later sections we will use the Bessel function J n of order n of the first kind and the Gauss hypergeometric series 2 F defined as and J n z 2 2 e iz sin tnt dt 2F ab;c;z n aa anbb bn z n. cc cnn! See Ref. 6 for details. Proof of Theorem : By 3 we have f x y f xy * 2 ˆ xy x y 2 In polar coordinates it is f x xy y ˆ xydx dy. 2 2 Medical Physics Vol. 3 No. 5 May 23
4 89 Ye Zhu and ang: A it theorem in CT 89 f x y 2 2 d f x r y r sin J 2r rdr r 2 d f x r y r sin J 2rdr 2 2 d t f x 2 y t sin 2 J tdt. 22 Recall the Parseval formula for absolute integrable functions f and g on : f tĝ tdt fˆ tgtdt. 23 It is known that 5 t; J t t t 2 t. 24 Denote the function on the right side as g(t). Then by the inversion formula of the Fourier ine transform J t 2 g t. 25 Consequently f x y 2 2 r d f x 2 y r sin 2 r d 2 f x 2 y r sin 2 t 2 gtdt 2 tdt 2 r d 2 f x 2 y r sin tdt 2 t t 2 f x y 2 2 d f x r y r sin tdt 2t t If for a fixed 2 (g(x r y r sin ) (t)d O(/t ) then we have 2 2 d gx r y r sin 2t tdt t 2 dt t t 2 27 which converges and tends to as. Consequently the contribution to 26 of the term g in the expression of f (x r y r sin ) in8 will tend to. As for the terms a j (k) ()r k in h j (r) on the right side of 8 we have the following computation: a j u j r xa j sin r jx sin r j x x a j ru j r x a j d dx r j x r j x x a 2 j r 2 u j r x a 2 j d2 sin r j x sin r jx dx 2 x a 3 j r 3 u j r x a 3 j d3 r jx r j x dx 3 x a 4 j r 4 u j r x a 4 j d4 sin r jx sin r j x dx 4 x a 5 j r 5 u j r x a 5 j d5 r j x r jx dx 5 x where u j (r) for r j ()rr j () and vanishes elsewhere. Consequently Medical Physics Vol. 3 No. 5 May 23
5 82 Ye Zhu and ang: A it theorem in CT 82 a k j r k u j r x a k j r x j k sin r j x r j k sin r j x 28 for k...5 where represents terms with x 2 or higher powers of x in the denominators. These omitted terms will contribute expressions tending to when asin 27. The leading term in 28 in turn contributes 2 a k 2 j r j k sin 2r d j t dt t 2 2 a k 2 j r j k d sin 2r j t dt t a k j r j k J 2r j d 2 2 a k j r j k J 2r j d. 29 Recall that we assumed that a j () and r j () are piecewise continuous and bounded: a j ()C 2 c(x y )r j () C for 2. Consequently using the asymptotic expansion J z 2 z z/4oz3/2 3 we obtain a bound O(C k C 2 k /c(x y )) for 29 which tends to as. The theorem then follows immediately. III. EXAMPLES A. Disk function Theorem covers the characteristic function on the unit disk f xy x2 y 2 3 otherwise. Its Fourier transform fˆw w sin S w J 2w 32 w is not absolutely integrable over the plane. Now 9 becomes According to we have f xy Q x y sin d 2 d J 2w 2wx y sin dw. 34 hen (xy)() we may express x y sin as x 2 y 2 sin( ) for a certain phase shifting. Changing variables we have f xy2 d J 2w2wx 2 y 2 sin dw d 2 J wwx 2 y 2 sin dw. Switch the order of integration f xy 2 J wdw wx 2 y 2 sin d 2 J wj wx 2 y 2 dw hen (xy)() 36 is also true. If x 2 y 2 we have the discontinuous integral of eber Schafheitlin as the it f xy J wj wx 2 y 2 dw x2 y 2 x 2 y 2. Therefore 37 f xy f xy 38 for points (xy) not on the unit circle. If x 2 y 2 we have 2 J wj wdw 2 J Hence for x 2 y 2 f xy 2 J Q t S wwe 2itw dw w w J 2we 2iwt dw 2 J 2w2twdw. 33 B. Delta function If we apply the band-ited reconstruction algorithm 9 and to the Dirac delta function f (xy)(xy) the resultant (xy) exhibits interesting singular behaviors. In fact ˆ w w sin S w. 4 Medical Physics Vol. 3 No. 5 May 23
6 82 Ye Zhu and ang: A it theorem in CT 82 If t we have Q t we 2itw dw sin2t sin2 t. 42 t 2 t 2 hen (xy)() xy Q x y sin d sin2x2 y 2 sin x 2 y 2 sin sin2 x 2 y 2 sin d 43 2 x 2 y 2 sin 2 by changing variables. Denote Ax 2 y 2 and tsin. The integrand then equals sin2at At sin2 At A 2 t 2 Integrating it termwise we obtain xy A n 2 2 n n 2At 2n. 2n!2n2 44 n A 2n n!n! A J 2A x 2 y J 2x 2 y 2 2 if (xy)(). On the other hand it is easy to find that For a given point (xy)() we may use the asymptotic expansion J z 2 z z 3 4 O z 3/2 47 and derive an asymptotic expansion of 45 as : xy x 2 y 2 3/4 2x 2 y O 5/2. 48 /2 x 2 y 2 Consequently in addition to () according to 46 (xy) will have a larger and larger magnitude of vibration at any fixed location (xy)() as. Nevertheless (xy) indeed approaches (xy) as in the following sense that for a function f (xy) asin Theorem f x xy y xydx dy f x y 49 where the point (x y ) satisfies the requirement in Theorem. This is because the integral on the left-hand side of 49 is actually equal to f (x y ) and Theorem applies. In particular 2 xydx dy d r J 2rr dr 2 J 2rdr J zdz. 5 It is interesting to compare (xy) with (xy) obtained by the Abel method of summability: 2 xy x 2 y /2 It is well known that (xy) tends to as at any point (xy)(). IV. DISCUSSIONS AND CONCLUSION There are relative merits and drawbacks associated with different it methods for handling the case that fˆ(u) is not absolutely integrable. Given comparable parameters the method of ited bandwidth tends to produce finer spatial resolution higher image noise and stronger edge ringing while the Abel method of summability may result in smoother profiles at the t of poorer spatial resolution. A heuristic way to define the comparability of the parameters with the two methods may be to request that the area after hard/soft spectral truncation be identical. Mathematically it is requested that wdw which is equivalent to we w dw or. 53 Because f (xy) is a partial sum of a Fourier series it exhibits the Gibbs phenomenon around discontinuous points. An example is a reconstructed unit disk function for as shown in Fig.. On the other hand the Gibbs phenomenon does not occur with the Abel method of summability as demonstrated by f (xy) with.2.44 in Fig. 2. However the sharpness of the slope on the unit circle is better in Fig. than that in Fig. 2 which indicates better spatial resolution in the former than that in the latter. Quantitatively the magnitude of the gradient of f (xy) on the unit circle is computed as follows: Medical Physics Vol. 3 No. 5 May 23
7 822 Ye Zhu and ang: A it theorem in CT 822 FIG.. Radial profile of a unit disk reconstructed using the method of ited bandwidth. FIG. 2. Radial profile of a unit disk reconstructed using the Abel method of summability. grad f r 2 J wj ww dw 2 J 2 ww dw J 2 2J 2J and its counterpart is grad f r e w/2 J wj ww dw e w/2 J 2 ww dw F ;3; F ;3; It can be shown that for 3 we have grad f rgrad f r The current choice by CT manufacturers is based on the method of ited bandwidth and allows various degrees of modifications on the associated ramp filter while the Abel method of summability remains an interesting theoretical possibility. The popularity of the method of ited bandwidth is due to the discrete detection technology and the Shannon sampling theorem. However it is well known that the Fourier spectrum spreads over the whole space for any spatial function finitely supported. Suppose that we can record projection data continuously or extrapolate the Fourier spectrum well beyond the Nyquist bandwidth imposed by the sampling rate the Abel method of summability may have some practical merits. In addition to the methods we have discussed other methods are also possible based on various other ways to do the intermediate integrals which should be equivalent to using approximate delta function of different types. 7 In conclusion for the method of ited bandwidth we have proved a pointwise it theorem for a general class of functions. Our results have improved the understanding on filtered backprojection in CT. Although our work has been done in the parallel beam geometry the findings are instructive for the divergent beam geometry. Further work is underway to extend the theory and explore its practical applications. ACKNOLEDGMENTS The work is partially supported by NIH Grant No. RO DC359 and a grant from the Carver Scientific Research Initative Grant Program. The authors thank Seung ook Lee John Meinel and Xiang Li with the CT/Micro-CT Lab University of Iowa for discussions and are grateful to anonymous reviewers for constructive critiques. a Electronic mail: yangbo-ye@uiowa.edu b Electronic mail: jiehua-zhu@uiowa.edu c Electronic mail: ge-wang@uiowa.edu A. C. Kak and M. Slaney Principles of Computerized Tomographic Imaging Society for Industrial and Applied Mathematics Philadelphia 2. 2 F. Natterer The Mathematics of Computerized Tomography Society for Industrial and Applied Mathematics Philadelphia 2. 3 E. M. Stein and G. eiss Introduction to Fourier Analysis on Euclidean Spaces Princeton University Press Princeton M. J. Lighthill An Introduction to Fourier Analysis and Generalized Functions Cambridge University Press Great Britain H. Bateman Tables of Integral Transforms McGraw-Hill New York 954 Vol. I. 6 H. Bateman Higher Transcendental Functions McGraw-Hill New York 953 Vol. II. 7 G.. ei Discrete singular convolution for the solution of Fokker- Plank equation J. Chem. Phys F. Natterer Numerical methods in tomography Medical Physics Vol. 3 No. 5 May 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 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 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 informationThe Diffusion Equation with Piecewise Smooth Initial Conditions ABSTRACT INTRODUCTION
Malaysian Journal of Mathematical Sciences 5(1): 101-110 (011) The Diffusion Equation with Piecewise Smooth Initial Conditions Ravshan Ashurov and Almaz Butaev Institute of Advanced Technology (ITMA),Universiti
More informationMA 201: Differentiation and Integration of Fourier Series Applications of Fourier Series Lecture - 10
MA 201: Differentiation and Integration of Fourier Series Applications of Fourier Series ecture - 10 Fourier Series: Orthogonal Sets We begin our treatment with some observations: For m,n = 1,2,3,... cos
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 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 informationAdvanced Training Course on FPGA Design and VHDL for Hardware Simulation and Synthesis. 26 October - 20 November, 2009
2065-33 Advanced Training Course on FPGA Design and VHDL for Hardware Simulation and Synthesis 26 October - 20 November, 2009 Introduction to two-dimensional digital signal processing Fabio Mammano University
More informationJUHA KINNUNEN. Partial Differential Equations
JUHA KINNUNEN Partial Differential Equations Department of Mathematics and Systems Analysis, Aalto University 207 Contents INTRODUCTION 2 FOURIER SERIES AND PDES 5 2. Periodic functions*.............................
More informationDRAFT - Math 101 Lecture Note - Dr. Said Algarni
2 Limits 2.1 The Tangent Problems The word tangent is derived from the Latin word tangens, which means touching. A tangent line to a curve is a line that touches the curve and a secant line is a line that
More informationMathematics of Physics and Engineering II: Homework answers You are encouraged to disagree with everything that follows. P R and i j k 2 1 1
Mathematics of Physics and Engineering II: Homework answers You are encouraged to disagree with everything that follows Homework () P Q = OQ OP =,,,, =,,, P R =,,, P S = a,, a () The vertex of the angle
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 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 informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationHarmonic Analysis on the Cube and Parseval s Identity
Lecture 3 Harmonic Analysis on the Cube and Parseval s Identity Jan 28, 2005 Lecturer: Nati Linial Notes: Pete Couperus and Neva Cherniavsky 3. Where we can use this During the past weeks, we developed
More informationMIT 2.71/2.710 Optics 10/31/05 wk9-a-1. The spatial frequency domain
10/31/05 wk9-a-1 The spatial frequency domain Recall: plane wave propagation x path delay increases linearly with x λ z=0 θ E 0 x exp i2π sinθ + λ z i2π cosθ λ z plane of observation 10/31/05 wk9-a-2 Spatial
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More information1 Review of di erential calculus
Review of di erential calculus This chapter presents the main elements of di erential calculus needed in probability theory. Often, students taking a course on probability theory have problems with concepts
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 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 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 informationComplex Analysis Slide 9: Power Series
Complex Analysis Slide 9: Power Series MA201 Mathematics III Department of Mathematics IIT Guwahati August 2015 Complex Analysis Slide 9: Power Series 1 / 37 Learning Outcome of this Lecture We learn Sequence
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 informationTopics. Example. Modulation. [ ] = G(k x. ] = 1 2 G ( k % k x 0) G ( k + k x 0) ] = 1 2 j G ( k x
Topics Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2008 CT/Fourier Lecture 3 Modulation Modulation Transfer Function Convolution/Multiplication Revisit Projection-Slice Theorem Filtered
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationbe the set of complex valued 2π-periodic functions f on R such that
. Fourier series. Definition.. Given a real number P, we say a complex valued function f on R is P -periodic if f(x + P ) f(x) for all x R. We let be the set of complex valued -periodic functions f on
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 142, Final Exam. 12/7/10.
Math 4, Final Exam. /7/0. No notes, calculator, or text. There are 00 points total. Partial credit may be given. Write your full name in the upper right corner of page. Number the pages in the upper right
More information2.2 The derivative as a Function
2.2 The derivative as a Function Recall: The derivative of a function f at a fixed number a: f a f a+h f(a) = lim h 0 h Definition (Derivative of f) For any number x, the derivative of f is f x f x+h f(x)
More informationTopics. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2006 CT/Fourier Lecture 2
Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2006 CT/Fourier Lecture 2 Topics Modulation Modulation Transfer Function Convolution/Multiplication Revisit Projection-Slice Theorem Filtered
More informationf(s) e -i n π s/l d s
Pointwise convergence of complex Fourier series Let f(x) be a periodic function with period l defined on the interval [,l]. The complex Fourier coefficients of f( x) are This leads to a Fourier series
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 informationInfinite Series. 1 Introduction. 2 General discussion on convergence
Infinite Series 1 Introduction I will only cover a few topics in this lecture, choosing to discuss those which I have used over the years. The text covers substantially more material and is available for
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 information1.6 Computing and Existence
1.6 Computing and Existence 57 1.6 Computing and Existence The initial value problem (1) y = f(x,y), y(x 0 ) = y 0 is studied here from a computational viewpoint. Answered are some basic questions about
More informationLecture 4 Filtering in the Frequency Domain. Lin ZHANG, PhD School of Software Engineering Tongji University Spring 2016
Lecture 4 Filtering in the Frequency Domain Lin ZHANG, PhD School of Software Engineering Tongji University Spring 2016 Outline Background From Fourier series to Fourier transform Properties of the Fourier
More informationLet R be the line parameterized by x. Let f be a complex function on R that is integrable. The Fourier transform ˆf = F f is. e ikx f(x) dx. (1.
Chapter 1 Fourier transforms 1.1 Introduction Let R be the line parameterized by x. Let f be a complex function on R that is integrable. The Fourier transform ˆf = F f is ˆf(k) = e ikx f(x) dx. (1.1) It
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 informationChapter 7: Applications of Integration
Chapter 7: Applications of Integration Fall 214 Department of Mathematics Hong Kong Baptist University 1 / 21 7.1 Volumes by Slicing Solids of Revolution In this section, we show how volumes of certain
More informationResearch Article Exact Interior Reconstruction with Cone-Beam CT
Biomedical Imaging Volume 2007, Article ID 10693, 5 pages doi:10.1155/2007/10693 Research Article Exact Interior Reconstruction with Cone-Beam CT Yangbo Ye, 1 Hengyong Yu, 2 and Ge Wang 2 1 Department
More informationSample Final Questions: Solutions Math 21B, Winter y ( y 1)(1 + y)) = A y + B
Sample Final Questions: Solutions Math 2B, Winter 23. Evaluate the following integrals: tan a) y y dy; b) x dx; c) 3 x 2 + x dx. a) We use partial fractions: y y 3 = y y ) + y)) = A y + B y + C y +. Putting
More information2 The Fourier Transform
2 The Fourier Transform The Fourier transform decomposes a signal defined on an infinite time interval into a λ-frequency component, where λ can be any real(or even complex) number. 2.1 Informal Development
More informationPractice Final Exam Solutions
Important Notice: To prepare for the final exam, study past exams and practice exams, and homeworks, quizzes, and worksheets, not just this practice final. A topic not being on the practice final does
More informationAn idea how to solve some of the problems. diverges the same must hold for the original series. T 1 p T 1 p + 1 p 1 = 1. dt = lim
An idea how to solve some of the problems 5.2-2. (a) Does not converge: By multiplying across we get Hence 2k 2k 2 /2 k 2k2 k 2 /2 k 2 /2 2k 2k 2 /2 k. As the series diverges the same must hold for the
More informationFourier Series. Spectral Analysis of Periodic Signals
Fourier Series. Spectral Analysis of Periodic Signals he response of continuous-time linear invariant systems to the complex exponential with unitary magnitude response of a continuous-time LI system at
More informationImage Reconstruction from Projection
Image Reconstruction from Projection Reconstruct an image from a series of projections X-ray computed tomography (CT) Computed tomography is a medical imaging method employing tomography where digital
More informationContinuous Time Signal Analysis: the Fourier Transform. Lathi Chapter 4
Continuous Time Signal Analysis: the Fourier Transform Lathi Chapter 4 Topics Aperiodic signal representation by the Fourier integral (CTFT) Continuous-time Fourier transform Transforms of some useful
More informationSensors. Chapter Signal Conditioning
Chapter 2 Sensors his chapter, yet to be written, gives an overview of sensor technology with emphasis on how to model sensors. 2. Signal Conditioning Sensors convert physical measurements into data. Invariably,
More informationMath 210, Final Exam, Fall 2010 Problem 1 Solution. v cosθ = u. v Since the magnitudes of the vectors are positive, the sign of the dot product will
Math, Final Exam, Fall Problem Solution. Let u,, and v,,3. (a) Is the angle between u and v acute, obtuse, or right? (b) Find an equation for the plane through (,,) containing u and v. Solution: (a) The
More informationTuesday, September 29, Page 453. Problem 5
Tuesday, September 9, 15 Page 5 Problem 5 Problem. Set up and evaluate the integral that gives the volume of the solid formed by revolving the region bounded by y = x, y = x 5 about the x-axis. Solution.
More informationMath 489AB A Very Brief Intro to Fourier Series Fall 2008
Math 489AB A Very Brief Intro to Fourier Series Fall 8 Contents Fourier Series. The coefficients........................................ Convergence......................................... 4.3 Convergence
More informationLecture 6 January 21, 2016
MATH 6/CME 37: Applied Fourier Analysis and Winter 06 Elements of Modern Signal Processing Lecture 6 January, 06 Prof. Emmanuel Candes Scribe: Carlos A. Sing-Long, Edited by E. Bates Outline Agenda: Fourier
More informationTopics in Harmonic Analysis Lecture 1: The Fourier transform
Topics in Harmonic Analysis Lecture 1: The Fourier transform Po-Lam Yung The Chinese University of Hong Kong Outline Fourier series on T: L 2 theory Convolutions The Dirichlet and Fejer kernels Pointwise
More informationQualification Exam: Mathematical Methods
Qualification Exam: Mathematical Methods Name:, QEID#41534189: August, 218 Qualification Exam QEID#41534189 2 1 Mathematical Methods I Problem 1. ID:MM-1-2 Solve the differential equation dy + y = sin
More informationGrade: The remainder of this page has been left blank for your workings. VERSION D. Midterm D: Page 1 of 12
First Name: Student-No: Last Name: Section: Grade: The remainder of this page has been left blank for your workings. Midterm D: Page of 2 Indefinite Integrals. 9 marks Each part is worth marks. Please
More informationMATH 5640: Fourier Series
MATH 564: Fourier Series Hung Phan, UMass Lowell September, 8 Power Series A power series in the variable x is a series of the form a + a x + a x + = where the coefficients a, a,... are real or complex
More information12.7 Heat Equation: Modeling Very Long Bars.
568 CHAP. Partial Differential Equations (PDEs).7 Heat Equation: Modeling Very Long Bars. Solution by Fourier Integrals and Transforms Our discussion of the heat equation () u t c u x in the last section
More informationMathematics of Physics and Engineering II: Homework problems
Mathematics of Physics and Engineering II: Homework problems Homework. Problem. Consider four points in R 3 : P (,, ), Q(,, 2), R(,, ), S( + a,, 2a), where a is a real number. () Compute the coordinates
More informationCalculus I Sample Exam #01
Calculus I Sample Exam #01 1. Sketch the graph of the function and define the domain and range. 1 a) f( x) 3 b) g( x) x 1 x c) hx ( ) x x 1 5x6 d) jx ( ) x x x 3 6 . Evaluate the following. a) 5 sin 6
More informationINTRODUCTION TO REAL ANALYSIS II MATH 4332 BLECHER NOTES
INTRODUCTION TO REAL ANALYSIS II MATH 433 BLECHER NOTES. As in earlier classnotes. As in earlier classnotes (Fourier series) 3. Fourier series (continued) (NOTE: UNDERGRADS IN THE CLASS ARE NOT RESPONSIBLE
More informationThe Fourier series for a 2π-periodic function
The Fourier series for a 2π-periodic function Let f : ( π, π] R be a bounded piecewise continuous function which we continue to be a 2π-periodic function defined on R, i.e. f (x + 2π) = f (x), x R. The
More informationIntroduction to Fourier Analysis
Lecture Introduction to Fourier Analysis Jan 7, 2005 Lecturer: Nati Linial Notes: Atri Rudra & Ashish Sabharwal. ext he main text for the first part of this course would be. W. Körner, Fourier Analysis
More information17 The functional equation
18.785 Number theory I Fall 16 Lecture #17 11/8/16 17 The functional equation In the previous lecture we proved that the iemann zeta function ζ(s) has an Euler product and an analytic continuation to the
More informationMath 127: Course Summary
Math 27: Course Summary Rich Schwartz October 27, 2009 General Information: M27 is a course in functional analysis. Functional analysis deals with normed, infinite dimensional vector spaces. Usually, these
More informationReal Analysis Problems
Real Analysis Problems Cristian E. Gutiérrez September 14, 29 1 1 CONTINUITY 1 Continuity Problem 1.1 Let r n be the sequence of rational numbers and Prove that f(x) = 1. f is continuous on the irrationals.
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 informationSeries Solution of Linear Ordinary Differential Equations
Series Solution of Linear Ordinary Differential Equations Department of Mathematics IIT Guwahati Aim: To study methods for determining series expansions for solutions to linear ODE with variable coefficients.
More informationLecture 13 - Wednesday April 29th
Lecture 13 - Wednesday April 29th jacques@ucsdedu Key words: Systems of equations, Implicit differentiation Know how to do implicit differentiation, how to use implicit and inverse function theorems 131
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 informationFOURIER SERIES WITH THE CONTINUOUS PRIMITIVE INTEGRAL
Real Analysis Exchange Summer Symposium XXXVI, 2012, pp. 30 35 Erik Talvila, Department of Mathematics and Statistics, University of the Fraser Valley, Chilliwack, BC V2R 0N3, Canada. email: Erik.Talvila@ufv.ca
More informationChapter 7: Techniques of Integration
Chapter 7: Techniques of Integration MATH 206-01: Calculus II Department of Mathematics University of Louisville last corrected September 14, 2013 1 / 43 Chapter 7: Techniques of Integration 7.1. Integration
More informationReal Analysis Qualifying Exam May 14th 2016
Real Analysis Qualifying Exam May 4th 26 Solve 8 out of 2 problems. () Prove the Banach contraction principle: Let T be a mapping from a complete metric space X into itself such that d(tx,ty) apple qd(x,
More informationMAC2313 Final A. (5 pts) 1. How many of the following are necessarily true? i. The vector field F = 2x + 3y, 3x 5y is conservative.
MAC2313 Final A (5 pts) 1. How many of the following are necessarily true? i. The vector field F = 2x + 3y, 3x 5y is conservative. ii. The vector field F = 5(x 2 + y 2 ) 3/2 x, y is radial. iii. All constant
More informationMathematical Methods for Physics and Engineering
Mathematical Methods for Physics and Engineering Lecture notes for PDEs Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 The integration theory
More informationFourier Series. Fourier Transform
Math Methods I Lia Vas Fourier Series. Fourier ransform Fourier Series. Recall that a function differentiable any number of times at x = a can be represented as a power series n= a n (x a) n where the
More informationMath 234 Exam 3 Review Sheet
Math 234 Exam 3 Review Sheet Jim Brunner LIST OF TOPIS TO KNOW Vector Fields lairaut s Theorem & onservative Vector Fields url Divergence Area & Volume Integrals Using oordinate Transforms hanging the
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 informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationMATH 6B Spring 2017 Vector Calculus II Study Guide Final Exam Chapters 8, 9, and Sections 11.1, 11.2, 11.7, 12.2, 12.3.
MATH 6B pring 2017 Vector Calculus II tudy Guide Final Exam Chapters 8, 9, and ections 11.1, 11.2, 11.7, 12.2, 12.3. Before starting with the summary of the main concepts covered in the quarter, here there
More informationHOMEWORK SOLUTIONS MATH 1910 Sections 6.4, 6.5, 7.1 Fall 2016
HOMEWORK SOLUTIONS MATH 9 Sections 6.4, 6.5, 7. Fall 6 Problem 6.4. Sketch the region enclosed by x = 4 y +, x = 4y, and y =. Use the Shell Method to calculate the volume of rotation about the x-axis SOLUTION.
More informationSolutions to Homework 11
Solutions to Homework 11 Read the statement of Proposition 5.4 of Chapter 3, Section 5. Write a summary of the proof. Comment on the following details: Does the proof work if g is piecewise C 1? Or did
More informationPhysics 6303 Lecture 11 September 24, LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation
Physics 6303 Lecture September 24, 208 LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation, l l l l l l. Consider problems that are no axisymmetric; i.e., the potential depends
More informationLecture 7: Edge Detection
#1 Lecture 7: Edge Detection Saad J Bedros sbedros@umn.edu Review From Last Lecture Definition of an Edge First Order Derivative Approximation as Edge Detector #2 This Lecture Examples of Edge Detection
More informationFall 2013 Hour Exam 2 11/08/13 Time Limit: 50 Minutes
Math 8 Fall Hour Exam /8/ Time Limit: 5 Minutes Name (Print): This exam contains 9 pages (including this cover page) and 7 problems. Check to see if any pages are missing. Enter all requested information
More informationClass Meeting # 2: The Diffusion (aka Heat) Equation
MATH 8.52 COURSE NOTES - CLASS MEETING # 2 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 2: The Diffusion (aka Heat) Equation The heat equation for a function u(, x (.0.). Introduction
More informationThe Approximation Problem
EE 508 Lecture 3 The Approximation Problem Classical Approximating Functions - Thompson and Bessel Approximations Review from Last Time Elliptic Filters Can be thought of as an extension of the CC approach
More informationSome Fun with Divergent Series
Some Fun with Divergent Series 1. Preliminary Results We begin by examining the (divergent) infinite series S 1 = 1 + 2 + 3 + 4 + 5 + 6 + = k=1 k S 2 = 1 2 + 2 2 + 3 2 + 4 2 + 5 2 + 6 2 + = k=1 k 2 (i)
More informationChapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12
Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
More information1. If 1, ω, ω 2, -----, ω 9 are the 10 th roots of unity, then (1 + ω) (1 + ω 2 ) (1 + ω 9 ) is A) 1 B) 1 C) 10 D) 0
4 INUTES. If, ω, ω, -----, ω 9 are the th roots of unity, then ( + ω) ( + ω ) ----- ( + ω 9 ) is B) D) 5. i If - i = a + ib, then a =, b = B) a =, b = a =, b = D) a =, b= 3. Find the integral values for
More informationFourier Series. 1. Review of Linear Algebra
Fourier Series In this section we give a short introduction to Fourier Analysis. If you are interested in Fourier analysis and would like to know more detail, I highly recommend the following book: Fourier
More informationHigher Order Linear Equations
C H A P T E R 4 Higher Order Linear Equations 4.1 1. The differential equation is in standard form. Its coefficients, as well as the function g(t) = t, are continuous everywhere. Hence solutions are valid
More informationAdvanced Mathematics Unit 2 Limits and Continuity
Advanced Mathematics 3208 Unit 2 Limits and Continuity NEED TO KNOW Expanding Expanding Expand the following: A) (a + b) 2 B) (a + b) 3 C) (a + b)4 Pascals Triangle: D) (x + 2) 4 E) (2x -3) 5 Random Factoring
More informationAdvanced Mathematics Unit 2 Limits and Continuity
Advanced Mathematics 3208 Unit 2 Limits and Continuity NEED TO KNOW Expanding Expanding Expand the following: A) (a + b) 2 B) (a + b) 3 C) (a + b)4 Pascals Triangle: D) (x + 2) 4 E) (2x -3) 5 Random Factoring
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
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 information3 Applications of partial differentiation
Advanced Calculus Chapter 3 Applications of partial differentiation 37 3 Applications of partial differentiation 3.1 Stationary points Higher derivatives Let U R 2 and f : U R. The partial derivatives
More informationChapter 17. Finite Volume Method The partial differential equation
Chapter 7 Finite Volume Method. This chapter focusses on introducing finite volume method for the solution of partial differential equations. These methods have gained wide-spread acceptance in recent
More informationGreen s Functions and Distributions
CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where
More informationPurdue University Study Guide for MA Credit Exam
Purdue University Study Guide for MA 16010 Credit Exam Students who pass the credit exam will gain credit in MA16010. The credit exam is a two-hour long exam with multiple choice questions. No books or
More information