Convolution Theorem Modulation Transfer Function y(t)=g(t)*h(t)

Size: px
Start display at page:

Download "Convolution Theorem Modulation Transfer Function y(t)=g(t)*h(t)"

Transcription

1 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 + = + h(t + + H(f Y(f=G(fH(f G(f f Bushberg et al 21 f 1

2 Convolution/Modulation Theorem [ ] e j 2πk xx F{ g(x h(x } = g(u h(x udu = g(u h(x u e j 2πk xx = g(uh(k x e j 2πk xu du = G(k x H(k x dxdu dx 2D Convolution/Multiplication Convolution [ ] = G(k x H(k x F g(x, y h(x,y Multiplication [ ] = G(k x H(k x F g(x, yh(x, y Convolution in the spatial domain transforms into multiplication in the frequency domain. Dual is modulation F{ g(xh(x } = G( k x H(k x Application of Convolution Thm. $ 1 x x <1 Λ(x = % & otherwise 1 F(Λ(x = 1 x e j 2πkxx dx =?? 1 ( Application of Convolution Thm. Λ(x = Π(x Π(x F(Λ(x = sinc 2 ( k x = * 2

3 Convolution Response of an Imaging System g(x,y G(k x h 1 (x,y h 2 (x,y h 3 (x,y MODULE 1 MODULE 2 MODULE 3 H 1 (k x H 2 (k x H 3 (k x z(x,y Z(k x z(x,y=g(x,y**h 1 (x,y**h 2 (x,y**h 3 (x,y Z(k x =G(k x H 1 (k x H 2 (k x H 3 (k x System MTF = Product of MTFs of Components Useful Approximation FWHM System = FWHM 12 + FWHM 22 + FWHM N 2 FWHM 1 =1mm FWHM 2 = 2mm FWHM system = 5 = 2.24mm Bushberg et al 21 PollEv.com/be28a 3

4 Modulation F g(xe j 2πk x [ ] = G(k x δ(k x k = G( k x k F[ g(xcos( 2πk x ] = 1 2 G ( k k x G ( k + k x F[ g(xsin( 2πk x ] = 1 2 j G ( k k x 1 2 j G ( k + k x PollEv.com/be28a Amplitude Modulation (e.g. AM Radio g(t 2g(t cos(2πf t Modulation x = 2cos(2πf t G(f -f f G(f-f + G(f+f * = 4

5 Summary of Basic Properties Linearity F [ ag(x, y+ bh(x, y ] = ag(k x + bh(k x Scaling F [ g(ax,by ] = 1 ab G! k x a, k $ x # & " b % Duality F{ G(x } = g( k x Shift F [ g(x a, y b ] = G(k x e j2π (k xa+k yb Convolution F [ g(x, y**h(x, y ] = G(k x H(k x Multiplication F [ g(x, yh(x, y ] = G(k x **H(k x Modulation F( g(x, ye j2π (xa+yb * + = G(k x a b In-class Exercise µ(x, y = rect(x, ycos( 2π (x + y Sketch the object Calculate and sketch the Fourier Transform of the Object G(ρ,θ = Projection-Slice Theorem g(l,θe j 2πρl dl = f (x,y δ(x cosθ + y sinθ le j 2πρl dx dy dl = f (x,y e j 2πρ x cosθ +y sinθ [ ] u= ρ cosθ,v= ρ sinθ = F 2D f (x,y ( dx dy g(x, y = Π(xΠ(y G(k x = sinc(k x sinc(k y y 1/2 (sinc/rect -1/2 1/2 x -1/2 Prince&Links 26 5

6 Projection-Slice Theorem Projection-Slice Theorem g(l, l l U(k x, = µ(x, ye j 2π (k x x +k y y dxdy [ ] = µ(x, ydy e j 2πk x x dx = g(x, e j 2πk x x dx = g(l, e j 2πkl dl In-Class : µ(x, y = cos2πx l g(l,θ [ ] U(k x = µ(x, ye j 2π (k x x +k y y dxdy = F 2D µ(x, y U(k x = G(k,θ k x = k cosθ k y = k sinθ k = k 2 2 x + k y F G(k,θ = l g(l,θe j 2πkl dl G(ρ,θ = Projection-Slice Theorem g(l,θe j 2πρl dl = f (x,y δ(x cosθ + y sinθ le j 2πρl dx dy dl = f (x,y e j 2πρ x cosθ +y sinθ [ ] u= ρ cosθ,v= ρ sinθ = F 2D f (x,y ( dx dy In-class Exercise µ(x, y = rect(x, ycos( 2π (x + y Sketch this object. What are the projections at theta = and 9 degrees? For what angle is the projection maximized? PollEv.com/be28a Prince&Links 26 6

7 Fourier Reconstruction Polar Version of Inverse FT F Interpolate onto Cartesian grid then take inverse transform µ(x, y = 2π G(k x e j 2π (k xx +k y y dk x dk y = G(k,θe j 2π (xk cosθ +yk sinθ kdkdθ π = G(k,θe j 2πk(x cosθ +y sinθ k dkdθ Note : g(l,θ + π = g( l,θ So G(k,θ + π = G( k,θ µ(x, y = Filtered Backprojection π π G(k,θe j 2π (xk cosθ +yk sinθ k dkdθ = k G(k,θe j 2πkl dkdθ π = g (l,θdθ where l = x cosθ + y sinθ g (l,θ = k G(k,θe j 2πkl dk = g(l,θ F 1 k = g(l,θ q(l [ ] Backproject a filtered projection Fourier Interpretation Density N circumference N 2π k Low frequencies are oversampled. So to compensate for this, multiply the k-space data by k before inverse transforming. Kak and Slaney; 7

8 Ram-Lak Filter Projection Reconstruction Path F k max =1/Δs x F -1 Filtered Projection Back- Project Projection Reconstruction Path * Filtered Projection Back- Project Kak and Slaney 8

9 Prince and Links 25 Prince and Links 25 9

Convolution/Modulation Theorem. 2D Convolution/Multiplication. Application of Convolution Thm. Bioengineering 280A Principles of Biomedical Imaging

Convolution/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 information

Topics. Example. Modulation. [ ] = G(k x. ] = 1 2 G ( k % k x 0) G ( k + k x 0) ] = 1 2 j G ( k x

Topics. 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 information

Topics. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2006 CT/Fourier Lecture 2

Topics. 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 information

Fourier Reconstruction. Polar Version of Inverse FT. Filtered Backprojection [ ] Bioengineering 280A Principles of Biomedical Imaging

Fourier Reconstruction. Polar Version of Inverse FT. Filtered Backprojection [ ] Bioengineering 280A Principles of Biomedical Imaging Fourier Reconstruction Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2014 CT/Fourier Lecture 5 F Interpolate onto Cartesian grid with appropriate density correction then take inverse

More information

Backprojection. Projections. Projections " $ & cosθ & Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2014 CT/Fourier Lecture 2

Backprojection. 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 information

Representation of 1D Function

Representation of 1D Function Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2005 Linear Systems Lecture 2 Representation of 1D Function From the sifting property, we can write a 1D function as g(x) = g(ξ)δ(x ξ)dξ.

More information

Topics. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2004 Lecture 4 2D Fourier Transforms

Topics. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2004 Lecture 4 2D Fourier Transforms Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2004 Lecture 4 2D Fourier Transforms Topics 1. 2D Signal Representations 2. 2D Fourier Transform 3. Transform Pairs 4. FT Properties 1

More information

Spin-Warp Pulse Sequence

Spin-Warp Pulse Sequence Bioengineering 28A Principles of Biomedical Imaging Fall Quarter 25 Linear Systems Lecture Spin-Warp Pulse Sequence RF G x (t G y (t k y k x K-space trajectories EPI Spiral k y k y k x k x Credit: Larry

More information

Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2005 X-Rays/CT Lecture 1. Topics

Bioengineering 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 information

Topics. EM spectrum. X-Rays Computed Tomography Direct Inverse and Iterative Inverse Backprojection Projection Theorem Filtered Backprojection

Topics. 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 information

Topics. EM spectrum. X-Rays Computed Tomography Direct Inverse and Iterative Inverse Backprojection Projection Theorem Filtered Backprojection

Topics. 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 information

BME I5000: Biomedical Imaging

BME I5000: Biomedical Imaging 1 BME I5000: Biomedical Imaging Lecture FT Fourier Transform Review Lucas C. Parra, parra@ccny.cuny.edu Blackboard: http://cityonline.ccny.cuny.edu/ Fourier Transform (1D) The Fourier Transform (FT) is

More information

Representation of 1D Function

Representation of 1D Function Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2005 Linear Systems Lecture 2 Representation of 1D Function From the sifting property, we can write a 1D function as g( = g(ξδ( ξdξ. To

More information

MIT 2.71/2.710 Optics 10/31/05 wk9-a-1. The spatial frequency domain

MIT 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 information

Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2006 X-Rays Lecture 3. Topics

Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2006 X-Rays Lecture 3. Topics Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2006 X-Rays Lecture 3 Topics Address questions from last lecture Review image equation Linearity Delta Functions Superposition Convolution

More information

EE 438 Essential Definitions and Relations

EE 438 Essential Definitions and Relations May 2004 EE 438 Essential Definitions and Relations CT Metrics. Energy E x = x(t) 2 dt 2. Power P x = lim T 2T T / 2 T / 2 x(t) 2 dt 3. root mean squared value x rms = P x 4. Area A x = x(t) dt 5. Average

More information

Solutions For Homework #5

Solutions For Homework #5 Solutions For Homework #5 Problem :[5 pts] ( ) The function f(x, y) = e π x a +y has a Fourier Transform as follows f(x, y) a e π(a u +v ) () which is obtained by applying the Similarity Theorem to the

More information

Discrete Fourier Transform

Discrete Fourier Transform Last lecture I introduced the idea that any function defined on x 0,..., N 1 could be written a sum of sines and cosines. There are two different reasons why this is useful. The first is a general one,

More information

Chain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics

Chain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics 3.33pt Chain Rule MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Spring 2019 Single Variable Chain Rule Suppose y = g(x) and z = f (y) then dz dx = d (f (g(x))) dx = f (g(x))g (x)

More information

HERMITE MULTIPLIERS AND PSEUDO-MULTIPLIERS

HERMITE MULTIPLIERS AND PSEUDO-MULTIPLIERS PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 124, Number 7, July 1996 HERMITE MULTIPLIERS AND PSEUDO-MULTIPLIERS JAY EPPERSON Communicated by J. Marshall Ash) Abstract. We prove a multiplier

More information

6 The Fourier transform

6 The Fourier transform 6 The Fourier transform In this presentation we assume that the reader is already familiar with the Fourier transform. This means that we will not make a complete overview of its properties and applications.

More information

Topics. View Aliasing. CT Sampling Requirements. Sampling Requirements in CT Sampling Theory Aliasing

Topics. View Aliasing. CT Sampling Requirements. Sampling Requirements in CT Sampling Theory Aliasing Topics Bioengineering 280A Principles of Biomedical Imaging Sampling Requirements in CT Sampling Theory Aliasing Fall Quarter 2007 CT/Fourier Lecture 4 CT Sampling Requirements View Aliasing What should

More information

Assignment for next week

Assignment for next week Assignment for next week Due: 6:30 PM Monday April 11 th (2016) -> PDF by email only 5 points: Pose two questions based on this lecture 5 points: what (and why) is sinc(0)? 5 points: ketch (by hand okay)

More information

MB-JASS D Reconstruction. Hannes Hofmann March 2006

MB-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 information

Biophysics 230: Nuclear Magnetic Resonance Haacke Chapter 11

Biophysics 230: Nuclear Magnetic Resonance Haacke Chapter 11 Biophysics 230: Nuclear Magnetic Resonance Haacke Chapter 11 Daniel B. Rowe, Ph.D. daniel.rowe@marquette.edu Department of Math, Stat, Comp Sci Marquette University dbrowe@mcw.edu Department of Biophysics

More information

Intro to harmonic analysis on groups Risi Kondor

Intro to harmonic analysis on groups Risi Kondor Risi Kondor Any (sufficiently smooth) function f on the unit circle (equivalently, any 2π periodic f ) can be decomposed into a sum of sinusoidal waves f(x) = k= c n e ikx c n = 1 2π f(x) e ikx dx 2π 0

More information

28. Fourier transforms in optics, part 3

28. Fourier transforms in optics, part 3 28. Fourier transforms in optics, part 3 Magnitude and phase some examples amplitude and phase of light waves what is the spectral phase, anyway? The Scale Theorem Defining the duration of a pulse the

More information

Fourier Transform. sin(n# x)), where! = 2" / L and

Fourier Transform. sin(n# x)), where! = 2 / L and Fourier Transform Henning Stahlberg Introduction The tools provided by the Fourier transform are helpful for the analysis of 1D signals (time and frequency (or Fourier) domains), as well as 2D/3D signals

More information

GBS765 Electron microscopy

GBS765 Electron microscopy GBS765 Electron microscopy Lecture 1 Waves and Fourier transforms 10/14/14 9:05 AM Some fundamental concepts: Periodicity! If there is some a, for a function f(x), such that f(x) = f(x + na) then function

More information

Midterm Review. Yao Wang Polytechnic University, Brooklyn, NY 11201

Midterm Review. Yao Wang Polytechnic University, Brooklyn, NY 11201 Midterm Review Yao Wang Polytechnic University, Brooklyn, NY 11201 Based on J. L. Prince and J. M. Links, Medical maging Signals and Systems, and lecture notes by Prince. Figures are from the textbook.

More information

Optics for Engineers Chapter 11

Optics for Engineers Chapter 11 Optics for Engineers Chapter 11 Charles A. DiMarzio Northeastern University Nov. 212 Fourier Optics Terminology Field Plane Fourier Plane C Field Amplitude, E(x, y) Ẽ(f x, f y ) Amplitude Point Spread

More information

Non-Uniform Fast Fourier Transformation (NUFFT) and Magnetic Resonace Imaging 2008 년 11 월 27 일. Lee, June-Yub (Ewha Womans University)

Non-Uniform Fast Fourier Transformation (NUFFT) and Magnetic Resonace Imaging 2008 년 11 월 27 일. Lee, June-Yub (Ewha Womans University) Non-Uniform Fast Fourier Transformation (NUFFT) and Magnetic Resonace Imaging 2008 년 11 월 27 일 Lee, June-Yub (Ewha Womans University) 1 NUFFT and MRI Nonuniform Fast Fourier Transformation (NUFFT) Basics

More information

Introduction to Cryptography

Introduction to Cryptography B504 / I538: Introduction to Cryptography Spring 2017 Lecture 15 Assignment 3 is due! Assignment 4 is out and is due in three weeks! 1 Recall: One-way functions (OWFs) Intuitively, a one-way function (OWF)

More information

Optics for Engineers Chapter 11

Optics for Engineers Chapter 11 Optics for Engineers Chapter 11 Charles A. DiMarzio Northeastern University Apr. 214 Fourier Optics Terminology Apr. 214 c C. DiMarzio (Based on Optics for Engineers, CRC Press) slides11r1 1 Fourier Optics

More information

Biomedical Engineering Image Formation II

Biomedical Engineering Image Formation II Biomedical Engineering Image Formation II PD Dr. Frank G. Zöllner Computer Assisted Clinical Medicine Medical Faculty Mannheim Fourier Series - A Fourier series decomposes periodic functions or periodic

More information

MAT 211 Final Exam. Spring Jennings. Show your work!

MAT 211 Final Exam. Spring Jennings. Show your work! MAT 211 Final Exam. pring 215. Jennings. how your work! Hessian D = f xx f yy (f xy ) 2 (for optimization). Polar coordinates x = r cos(θ), y = r sin(θ), da = r dr dθ. ylindrical coordinates x = r cos(θ),

More information

3 Applications of partial differentiation

3 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 information

E2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2)

E2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2) E.5 Signals & Linear Systems Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & ) 1. Sketch each of the following continuous-time signals, specify if the signal is periodic/non-periodic,

More information

1 Assignment 1: Nonlinear dynamics (due September

1 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 information

Image Reconstruction from Projection

Image 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 information

The Fourier transform, indicated by the operator F, constructs a spectrum A(k x,k y ) = F {E(x,y)} from a spatial distribution E(x,y):

The Fourier transform, indicated by the operator F, constructs a spectrum A(k x,k y ) = F {E(x,y)} from a spatial distribution E(x,y): 1 4. FOURIER TRANSFORMS (vs. 3.0) DEFINITIONS GENERAL PROPERTIES linearity conjugate functions multiplication and convolution multiplication and correlation Parseval's relation time shifting frequency

More information

conditional cdf, conditional pdf, total probability theorem?

conditional cdf, conditional pdf, total probability theorem? 6 Multiple Random Variables 6.0 INTRODUCTION scalar vs. random variable cdf, pdf transformation of a random variable conditional cdf, conditional pdf, total probability theorem expectation of a random

More information

Spin. Nuclear Spin Rules

Spin. Nuclear Spin Rules Spin Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 203 MRI Lecture Intrinsic angular momentum of elementary particles -- electrons, protons, neutrons. Spin is quantized. Key concept

More information

ECON 186 Class Notes: Optimization Part 2

ECON 186 Class Notes: Optimization Part 2 ECON 186 Class Notes: Optimization Part 2 Jijian Fan Jijian Fan ECON 186 1 / 26 Hessians The Hessian matrix is a matrix of all partial derivatives of a function. Given the function f (x 1,x 2,...,x n ),

More information

G52IVG, School of Computer Science, University of Nottingham

G52IVG, School of Computer Science, University of Nottingham Image Transforms Fourier Transform Basic idea 1 Image Transforms Fourier transform theory Let f(x) be a continuous function of a real variable x. The Fourier transform of f(x) is F ( u) f ( x)exp[ j2πux]

More information

1,3. f x x f x x. Lim. Lim. Lim. Lim Lim. y 13x b b 10 b So the equation of the tangent line is y 13x

1,3. f x x f x x. Lim. Lim. Lim. Lim Lim. y 13x b b 10 b So the equation of the tangent line is y 13x 1.5 Topics: The Derivative lutions 1. Use the limit definition of derivative (the one with x in it) to find f x given f x 4x 5x 6 4 x x 5 x x 6 4x 5x 6 f x x f x f x x0 x x0 x xx x x x x x 4 5 6 4 5 6

More information

Math 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 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 information

worked out from first principles by parameterizing the path, etc. If however C is a A path C is a simple closed path if and only if the starting point

worked out from first principles by parameterizing the path, etc. If however C is a A path C is a simple closed path if and only if the starting point III.c Green s Theorem As mentioned repeatedly, if F is not a gradient field then F dr must be worked out from first principles by parameterizing the path, etc. If however is a simple closed path in the

More information

Fourier Sampling. Fourier Sampling. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2005 Linear Systems Lecture 3.

Fourier Sampling. Fourier Sampling. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2005 Linear Systems Lecture 3. Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2005 Linear Sstems Lecture 3 Fourier Sampling F Instead of sampling the signal, we sample its Fourier Transform Sample??? F -1 Fourier

More information

Final Exam SOLUTIONS MAT 131 Fall 2011

Final Exam SOLUTIONS MAT 131 Fall 2011 1. Compute the following its. (a) Final Exam SOLUTIONS MAT 131 Fall 11 x + 1 x 1 x 1 The numerator is always positive, whereas the denominator is negative for numbers slightly smaller than 1. Also, as

More information

Final Exam Review Quesitons

Final Exam Review Quesitons Final Exam Review Quesitons. Compute the following integrals. (a) x x 4 (x ) (x + 4) dx. The appropriate partial fraction form is which simplifies to x x 4 (x ) (x + 4) = A x + B (x ) + C x + 4 + Dx x

More information

F(u) = f e(t)cos2πutdt + = f e(t)cos2πutdt i f o(t)sin2πutdt = F e(u) + if o(u).

F(u) = f e(t)cos2πutdt + = f e(t)cos2πutdt i f o(t)sin2πutdt = F e(u) + if o(u). The (Continuous) Fourier Transform Matrix Computations and Virtual Spaces Applications in Signal/Image Processing Niclas Börlin niclas.borlin@cs.umu.se Department of Computing Science, Umeå University,

More information

Mathematics for Physical Sciences III

Mathematics for Physical Sciences III Mathematics for Physical Sciences III Change of lecturer: First 4 weeks: myself again! Remaining 8 weeks: Dr Stephen O Sullivan Continuous Assessment Test Date to be announced (probably Week 7 or 8) -

More information

CS-9645 Introduction to Computer Vision Techniques Winter 2018

CS-9645 Introduction to Computer Vision Techniques Winter 2018 Table of Contents Spectral Analysis...1 Special Functions... 1 Properties of Dirac-delta Functions...1 Derivatives of the Dirac-delta Function... 2 General Dirac-delta Functions...2 Harmonic Analysis...

More information

Random Signals and Systems. Chapter 3. Jitendra K Tugnait. Department of Electrical & Computer Engineering. Auburn University.

Random Signals and Systems. Chapter 3. Jitendra K Tugnait. Department of Electrical & Computer Engineering. Auburn University. Random Signals and Systems Chapter 3 Jitendra K Tugnait Professor Department of Electrical & Computer Engineering Auburn University Two Random Variables Previously, we only dealt with one random variable

More information

Integration Techniques

Integration Techniques Review for the Final Exam - Part - Solution Math Name Quiz Section The following problems should help you review for the final exam. Don t hesitate to ask for hints if you get stuck. Integration Techniques.

More information

2017 HSC Mathematics Extension 1 Marking Guidelines

2017 HSC Mathematics Extension 1 Marking Guidelines 07 HSC Mathematics Extension Marking Guidelines Section I Multiple-choice Answer Key Question Answer A B 3 B 4 C 5 D 6 D 7 A 8 C 9 C 0 B NESA 07 HSC Mathematics Extension Marking Guidelines Section II

More information

NIH Public Access Author Manuscript Proc Soc Photo Opt Instrum Eng. Author manuscript; available in PMC 2013 December 17.

NIH 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 information

DISCRETE FOURIER TRANSFORM

DISCRETE FOURIER TRANSFORM DD2423 Image Processing and Computer Vision DISCRETE FOURIER TRANSFORM Mårten Björkman Computer Vision and Active Perception School of Computer Science and Communication November 1, 2012 1 Terminology:

More information

1D Convolution. Convolution g[m] = g[0]δ[m] + g[1]δ[m 1] + g[2]δ[m 2] h[m',k] = L[δ[m k]] = h[ m $ k]

1D 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 information

Advanced Training Course on FPGA Design and VHDL for Hardware Simulation and Synthesis. 26 October - 20 November, 2009

Advanced 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 information

Saturation Rates for Filtered Back Projection

Saturation 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 information

51. General Surface Integrals

51. General Surface Integrals 51. General urface Integrals The area of a surface in defined parametrically by r(u, v) = x(u, v), y(u, v), z(u, v) over a region of integration in the input-variable plane is given by d = r u r v da.

More information

Starting from Heat Equation

Starting from Heat Equation Department of Applied Mathematics National Chiao Tung University Hsin-Chu 30010, TAIWAN 20th August 2009 Analytical Theory of Heat The differential equations of the propagation of heat express the most

More information

Created by T. Madas LINE INTEGRALS. Created by T. Madas

Created by T. Madas LINE INTEGRALS. Created by T. Madas LINE INTEGRALS LINE INTEGRALS IN 2 DIMENSIONAL CARTESIAN COORDINATES Question 1 Evaluate the integral ( x + 2y) dx, C where C is the path along the curve with equation y 2 = x + 1, from ( ) 0,1 to ( )

More information

Name: Instructor: Multiple Choice. x 3. = lim x 3 x 3 x (x 2 + 7) 16 = lim. (x 3)( x ) x 3 (x 3)( x ) = lim.

Name: Instructor: Multiple Choice. x 3. = lim x 3 x 3 x (x 2 + 7) 16 = lim. (x 3)( x ) x 3 (x 3)( x ) = lim. Multiple Choice 1.(6 pts.) Evaluate the following limit: x + 7 4 lim. x 3 x 3 lim x 3 x + 7 4 x 3 x + 7 4 x + 7 + 4 x 3 x 3 x + 7 + 4 (x + 7) 16 x 3 (x 3)( x + 7 + 4) x 9 x 3 (x 3)( x + 7 + 4) x 3 (x 3)(x

More information

FRACTIONAL DUAL SOLUTIONS AND CORRESPONDING SOURCES

FRACTIONAL DUAL SOLUTIONS AND CORRESPONDING SOURCES Progress In Electromagnetics Research, PIER 5, 3 38, 000 FRACTIONAL DUAL SOLUTIONS AND CORRESPONDING SOURCES Q. A. Naqvi and A. A. Rizvi Communications Lab. Department of Electronics Quaidi-i-Azam University

More information

e x2 dxdy, e x2 da, e x2 x 3 dx = e

e x2 dxdy, e x2 da, e x2 x 3 dx = e STS26-4 Calculus II: The fourth exam Dec 15, 214 Please show all your work! Answers without supporting work will be not given credit. Write answers in spaces provided. You have 1 hour and 2minutes to complete

More information

Images have structure at various scales

Images have structure at various scales Images have structure at various scales Frequency Frequency of a signal is how fast it changes Reflects scale of structure A combination of frequencies 0.1 X + 0.3 X + 0.5 X = Fourier transform Can we

More information

Solution. This is a routine application of the chain rule.

Solution. This is a routine application of the chain rule. EXAM 2 SOLUTIONS 1. If z = e r cos θ, r = st, θ = s 2 + t 2, find the partial derivatives dz ds chain rule. Write your answers entirely in terms of s and t. dz and dt using the Solution. This is a routine

More information

Derivatives and Integrals

Derivatives and Integrals Derivatives and Integrals Definition 1: Derivative Formulas d dx (c) = 0 d dx (f ± g) = f ± g d dx (kx) = k d dx (xn ) = nx n 1 (f g) = f g + fg ( ) f = f g fg g g 2 (f(g(x))) = f (g(x)) g (x) d dx (ax

More information

Solutions of Math 53 Midterm Exam I

Solutions of Math 53 Midterm Exam I Solutions of Math 53 Midterm Exam I Problem 1: (1) [8 points] Draw a direction field for the given differential equation y 0 = t + y. (2) [8 points] Based on the direction field, determine the behavior

More information

Hamburger Beiträge zur Angewandten Mathematik

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 information

Computational Methods for Astrophysics: Fourier Transforms

Computational Methods for Astrophysics: Fourier Transforms Computational Methods for Astrophysics: Fourier Transforms John T. Whelan (filling in for Joshua Faber) April 27, 2011 John T. Whelan April 27, 2011 Fourier Transforms 1/13 Fourier Analysis Outline: Fourier

More information

BioE Exam 1 10/9/2018 Answer Sheet - Correct answer is A for all questions. 1. The sagittal plane

BioE Exam 1 10/9/2018 Answer Sheet - Correct answer is A for all questions. 1. The sagittal plane BioE 1330 - Exam 1 10/9/2018 Answer Sheet - Correct answer is A for all questions 1. The sagittal plane A. is perpendicular to the coronal plane. B. is parallel to the top of the head. C. represents a

More information

Math 425 Notes 9. We state a rough version of the fundamental theorem of calculus. Almost all calculations involving integrals lead back to this.

Math 425 Notes 9. We state a rough version of the fundamental theorem of calculus. Almost all calculations involving integrals lead back to this. Multiple Integrals: efintion Math 425 Notes 9 We state a rough version of the fundamental theorem of calculus. Almost all calculations involving integrals lead back to this. efinition 1. Let f : R R and

More information

Volumes of Solids of Revolution Lecture #6 a

Volumes of Solids of Revolution Lecture #6 a Volumes of Solids of Revolution Lecture #6 a Sphereoid Parabaloid Hyperboloid Whateveroid Volumes Calculating 3-D Space an Object Occupies Take a cross-sectional slice. Compute the area of the slice. Multiply

More information

Discrete Fourier Transform

Discrete Fourier Transform Discrete Fourier Transform DD2423 Image Analysis and Computer Vision Mårten Björkman Computational Vision and Active Perception School of Computer Science and Communication November 13, 2013 Mårten Björkman

More information

Implicit Functions, Curves and Surfaces

Implicit Functions, Curves and Surfaces Chapter 11 Implicit Functions, Curves and Surfaces 11.1 Implicit Function Theorem Motivation. In many problems, objects or quantities of interest can only be described indirectly or implicitly. It is then

More information

Dimensions = xyz dv. xyz dv as an iterated integral in rectangular coordinates.

Dimensions = xyz dv. xyz dv as an iterated integral in rectangular coordinates. Math Show Your Work! Page of 8. () A rectangular box is to hold 6 cubic meters. The material used for the top and bottom of the box is twice as expensive per square meter than the material used for the

More information

Books. B. Borden, Radar Imaging of Airborne Targets, Institute of Physics, 1999.

Books. B. Borden, Radar Imaging of Airborne Targets, Institute of Physics, 1999. Books B. Borden, Radar Imaging of Airborne Targets, Institute of Physics, 1999. C. Elachi, Spaceborne Radar Remote Sensing: Applications and Techniques, IEEE Press, New York, 1987. W. C. Carrara, R. G.

More information

Higher order derivative

Higher order derivative 2 Î 3 á Higher order derivative Î 1 Å Iterated partial derivative Iterated partial derivative Suppose f has f/ x, f/ y and ( f ) = 2 f x x x 2 ( f ) = 2 f x y x y ( f ) = 2 f y x y x ( f ) = 2 f y y y

More information

LESSON 25: LAGRANGE MULTIPLIERS OCTOBER 30, 2017

LESSON 25: LAGRANGE MULTIPLIERS OCTOBER 30, 2017 LESSON 5: LAGRANGE MULTIPLIERS OCTOBER 30, 017 Lagrange multipliers is another method of finding minima and maxima of functions of more than one variable. In fact, many of the problems from the last homework

More information

Signal Analysis. Filter Banks and. One application for filter banks is to decompose the input signal into different bands or channels

Signal Analysis. Filter Banks and. One application for filter banks is to decompose the input signal into different bands or channels Filter banks Multi dimensional Signal Analysis A common type of processing unit for discrete signals is a filter bank, where some input signal is filtered by n filters, producing n channels Channel 1 Lecture

More information

Therefore the new Fourier coefficients are. Module 2 : Signals in Frequency Domain Problem Set 2. Problem 1

Therefore the new Fourier coefficients are. Module 2 : Signals in Frequency Domain Problem Set 2. Problem 1 Module 2 : Signals in Frequency Domain Problem Set 2 Problem 1 Let be a periodic signal with fundamental period T and Fourier series coefficients. Derive the Fourier series coefficients of each of the

More information

2.2 Separable Equations

2.2 Separable Equations 2.2 Separable Equations Definition A first-order differential equation that can be written in the form Is said to be separable. Note: the variables of a separable equation can be written as Examples Solve

More information

Limit Laws- Theorem 2.2.1

Limit Laws- Theorem 2.2.1 Limit Laws- Theorem 2.2.1 If L, M, c, and k are real numbers and lim x c f (x) =L and lim g(x) =M, then x c 1 Sum Rule: lim x c (f (x)+g(x)) = L + M 2 Difference Rule: lim x c (f (x) g(x)) = L M 3 Constant

More information

Math 210, Final Exam, Spring 2012 Problem 1 Solution. (a) Find an equation of the plane passing through the tips of u, v, and w.

Math 210, Final Exam, Spring 2012 Problem 1 Solution. (a) Find an equation of the plane passing through the tips of u, v, and w. Math, Final Exam, Spring Problem Solution. Consider three position vectors (tails are the origin): u,, v 4,, w,, (a) Find an equation of the plane passing through the tips of u, v, and w. (b) Find an equation

More information

3F1 Random Processes Examples Paper (for all 6 lectures)

3F1 Random Processes Examples Paper (for all 6 lectures) 3F Random Processes Examples Paper (for all 6 lectures). Three factories make the same electrical component. Factory A supplies half of the total number of components to the central depot, while factories

More information

Math 121A: Homework 6 solutions

Math 121A: Homework 6 solutions Math A: Homework 6 solutions. (a) The coefficients of the Fourier sine series are given by b n = π f (x) sin nx dx = x(π x) sin nx dx π = (π x) cos nx dx nπ nπ [x(π x) cos nx]π = n ( )(sin nx) dx + π n

More information

Continuous-Time Fourier Transform

Continuous-Time Fourier Transform Signals and Systems Continuous-Time Fourier Transform Chang-Su Kim continuous time discrete time periodic (series) CTFS DTFS aperiodic (transform) CTFT DTFT Lowpass Filtering Blurring or Smoothing Original

More information

E X A M. Probability Theory and Stochastic Processes Date: December 13, 2016 Duration: 4 hours. Number of pages incl.

E X A M. Probability Theory and Stochastic Processes Date: December 13, 2016 Duration: 4 hours. Number of pages incl. E X A M Course code: Course name: Number of pages incl. front page: 6 MA430-G Probability Theory and Stochastic Processes Date: December 13, 2016 Duration: 4 hours Resources allowed: Notes: Pocket calculator,

More information

MATH 18.01, FALL PROBLEM SET # 6 SOLUTIONS

MATH 18.01, FALL PROBLEM SET # 6 SOLUTIONS MATH 181, FALL 17 - PROBLEM SET # 6 SOLUTIONS Part II (5 points) 1 (Thurs, Oct 6; Second Fundamental Theorem; + + + + + = 16 points) Let sinc(x) denote the sinc function { 1 if x =, sinc(x) = sin x if

More information

Math 350 Solutions for Final Exam Page 1. Problem 1. (10 points) (a) Compute the line integral. F ds C. z dx + y dy + x dz C

Math 350 Solutions for Final Exam Page 1. Problem 1. (10 points) (a) Compute the line integral. F ds C. z dx + y dy + x dz C Math 35 Solutions for Final Exam Page Problem. ( points) (a) ompute the line integral F ds for the path c(t) = (t 2, t 3, t) with t and the vector field F (x, y, z) = xi + zj + xk. (b) ompute the line

More information

Math 165 Final Exam worksheet solutions

Math 165 Final Exam worksheet solutions C Roettger, Fall 17 Math 165 Final Exam worksheet solutions Problem 1 Use the Fundamental Theorem of Calculus to compute f(4), where x f(t) dt = x cos(πx). Solution. From the FTC, the derivative of the

More information

Lesson 3: Linear differential equations of the first order Solve each of the following differential equations by two methods.

Lesson 3: Linear differential equations of the first order Solve each of the following differential equations by two methods. Lesson 3: Linear differential equations of the first der Solve each of the following differential equations by two methods. Exercise 3.1. Solution. Method 1. It is clear that y + y = 3 e dx = e x is an

More information

Find the indicated derivative. 1) Find y(4) if y = 3 sin x. A) y(4) = 3 cos x B) y(4) = 3 sin x C) y(4) = - 3 cos x D) y(4) = - 3 sin x

Find the indicated derivative. 1) Find y(4) if y = 3 sin x. A) y(4) = 3 cos x B) y(4) = 3 sin x C) y(4) = - 3 cos x D) y(4) = - 3 sin x Assignment 5 Name Find the indicated derivative. ) Find y(4) if y = sin x. ) A) y(4) = cos x B) y(4) = sin x y(4) = - cos x y(4) = - sin x ) y = (csc x + cot x)(csc x - cot x) ) A) y = 0 B) y = y = - csc

More information

Partial Derivatives Formulas. KristaKingMath.com

Partial Derivatives Formulas. KristaKingMath.com Partial Derivatives Formulas KristaKingMath.com Domain and range of a multivariable function A function f of two variables is a rule that assigns to each ordered pair of real numbers (x, y) in a set D

More information

Chapter 4 Notes, Calculus I with Precalculus 3e Larson/Edwards

Chapter 4 Notes, Calculus I with Precalculus 3e Larson/Edwards 4.1 The Derivative Recall: For the slope of a line we need two points (x 1,y 1 ) and (x 2,y 2 ). Then the slope is given by the formula: m = y x = y 2 y 1 x 2 x 1 On a curve we can find the slope of a

More information

Math 201 Solutions to Assignment 1. 2ydy = x 2 dx. y = C 1 3 x3

Math 201 Solutions to Assignment 1. 2ydy = x 2 dx. y = C 1 3 x3 Math 201 Solutions to Assignment 1 1. Solve the initial value problem: x 2 dx + 2y = 0, y(0) = 2. x 2 dx + 2y = 0, y(0) = 2 2y = x 2 dx y 2 = 1 3 x3 + C y = C 1 3 x3 Notice that y is not defined for some

More information