Representation of 1D Function
|
|
- Ilene Mason
- 5 years ago
- Views:
Transcription
1 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ξ. To gain intuition, consider the approximation g(x) n= g(nδx) 1 Δx Π x nδx Δx. Δx g(x) 1
2 Representation of 2D Function Similarly, we can write a 2D function as g(x, y) = g(ξ,η)δ(x ξ,y η)dξdη. To gain intuition, consider the approximation g(x, y) g(nδx,mδy) 1 Δx Π x nδx 1 n= Δx Δy Π y mδy ΔxΔy. m= Δy Impulse Response Intuition: the impulse response is the response of a system to an input of infinitesimal width and unit area. Original Image Blurred Image Since any input can be thought of as the weighted sum of impulses, a linear system is characterized by its impulse response(s). 2
3 Impulse Response The impulse response characterizes the response of a system over all space to a Dirac delta impulse function at a certain location. h(x 2 ;ξ) = L[ δ( x 1 ξ) ] 1D Impulse Response [ ( )] 2D Impulse Response h(x 2, y 2 ;ξ,η) = L δ x 1 ξ, y 1 η y 1 Impulse at ξ,η y 2 h(x2,y 2 ;ξ,η) x 1 x 2 Superposition Integral What is the response to an arbitrary function g(x 1,y 1 )? Write g(x 1,y 1 ) = g(ξ,η)δ(x 1 ξ, y 1 η)dξdη. - - The response is given by [ ] I(x 2,y 2 ) = L g 1 (x 1,y 1 ) [ - - ] = L g(ξ,η)δ(x 1 ξ,y 1 η)dξdη [ ] = g(ξ,η)l δ(x 1 ξ, y 1 η) dξdη - - = g(ξ,η)h(x 2,y 2 ;ξ,η) dξdη - - 3
4 Space Invariance If a system is space invariant, the impulse response depends only on the difference between the output coordinates and the position of ( ) the impulse and is given by h(x 2, y 2 ;ξ,η) = h x 2 ξ,y 2 η I(x 2,y 2 ) = 2D Convolution For a space invariant linear system, the superposition integral becomes a convolution integral. - - g(ξ,η)h(x 2,y 2 ;ξ,η) dξdη = g(ξ,η)h(x 2 ξ, y 2 η) dξdη - - = g(x 2,y 2 ) **h(x 2,y 2 ) where ** denotes 2D convolution. This will sometimes be abbreviated as *, e.g. I(x 2, y 2 )= g(x 2, y 2 )*h(x 2, y 2 ). 4
5 I(x) = 1D Convolution For completeness, here is the 1D version. - g(ξ)h(x;ξ)dξ = g(ξ)h(x ξ) dξ - = g(x) h(x) Convolution with Dirac delta function g(x) δ(x) = g(ξ)δ(x ξ) dξ - = g(x) g(x) δ(x Δ) = g(ξ)δ(x Δ ξ) dξ - = g(x Δ) Convolution of g(x) with a shifted Dirac delta function just shifts g(x) g(x,y) δ(x,y) = g(ξ,η)δ(x ξ,y η) dξdη - - = g(x,y) g(x,y) δ(x x 0,y y 0 ) = g(ξ,η)δ(x x 0 ξ,y y 0 η) dξdη - - = g(x x 0,y y 0 ) 5
6 2D Convolution Example g(x)= δ(x+1/2,y) + δ(x,y) y h(x)=rect(x,y) y x -1/2 12 x I(x,y)=g(x)**h(x,y) x 2D Convolution Example 6
7 1D Convolution Review g(x) h(x) = - g(ξ)h(x ξ)dξ Basic Rule: Flip one function, slide it past the other function, and integrate as you go. g(x)=rect(x) h(x)=rect(x-1/2) -1/2 1/2 1 1D Convolution Review h(-1/2-ξ) g(ξ) I(x) h(-ξ) h(1/2-ξ) -1/2 1/2 3/2 h(3/2-ξ) 7
8 Convolution/Modulation Theorem { } = g(u) h(x u)du F g(x) h(x) [ ] e j 2πk xx = g(u) h(x u) e j 2πk xx = g(u)h(k x )e j 2πk xu du = G(k x )H(k x ) dxdu dx Convolution in the spatial domain transforms into multiplication in the frequency domain. Dual is modulation F{ g(x)h(x) } = G( k x ) H(k x ) Convolution 2D Convolution/Multiplication [ ] = G(k x,k y )H(k x,k y ) F g(x, y) h(x,y) Multiplication [ ] = G(k x,k y ) H(k x,k y ) F g(x, y)h(x, y) 8
9 Application of Convolution Thm. 1 x x <1 Λ(x) = 0 otherwise F(Λ(x)) = 1 ( 1 x ) e j 2πkxx dx =?? Application of Convolution Thm. Λ(x) = Π(x) Π(x) F(Λ(x)) = sinc 2 ( ) k x -1 1 = * 9
10 Eigenfunctions The fundamental nature of the convolution theorem may be better understood by observing that the complex exponentials are eigenfunctions of the convolution operator. e j 2πk xx g(x) z(x) z(x) = g(x) e j 2πk x x = g(u) e j 2πk x (x u) du = G(k x )e j 2πk x x The response of a linear shift invariant system to a complex exponential is simply the exponential multiplied by the FT of the system s impulse response. Eigenfunctions Now consider an arbitrary input h(x). h(x) g(x) z(x) Recall that we can express h(x) as the integral of weighted complex exponentials. h(x) = H(k x )e j 2πk x x dk x Each of these exponentials is weighted by G(k x ) so that the response may be written as z(x) = G(k x )H(k x )e j 2πk x x dk x 10
11 Convolution Example Analog vs. Digital The Analog World: Continuous time/space, continuous valued signals or images, e.g. vinyl records, photographs, x-ray films. The Digital World: Discrete time/space, discrete-valued signals or images, e.g. CD-Roms, DVDs, digital photos, digital x-rays, CT, MRI, ultrasound. 11
12 The Process of Sampling g(x) Δx sample g[n]=g(n Δx) Questions How finely do we need to sample? What happens if we don t sample finely enough? Can we reconstruct the original signal or image from its samples? 12
13 Sampling in the Time Domain Sampling in Image Space 13
14 Sampling in k-space Sampling in k-space 14
15 Comb Function comb(x) = δ(x n) n= x Other names: Impulse train, bed of nails, shah function. Scaled Comb Function comb x = Δx n= δ( x Δx n) = δ( x nδx ) Δx n= = Δx δ(x nδx) n= x Δx 15
16 1D spatial sampling g S (x) = g(x) 1 Δx comb x Δx Recall the sifting property But we can also write = g(x) δ(x nδx) n= = g(nδx)δ(x nδx) n= So, g(x)δ(x a) = g(a)δ(x a) g(x)δ(x a) = g(a) g(a)δ(x a) = g(a) δ(x a) = g(a) 1D spatial sampling g(x) comb(x/δx)/ Δx x Δx g S (x) x 16
17 F comb(x) Fourier Transform of comb(x) [ ] = comb(k x ) = δ(k x n) n= F 1 Δx comb( x Δx ) = 1 Δx Δxcomb(k xδx) = δ(k x Δx n) n= = 1 Δx n= δ(k x n Δx ) Fourier Transform of comb(x/ Δx) comb(x/ Δx)/ Δx Δx F x 1/Δx comb(k x Δx) 1/Δx k x 17
18 F g S (x) Fourier Transform of g S (x) Δx [ ] = F g(x) 1 Δx comb x = G(k x ) F 1 Δx comb x Δx = G(k x ) 1 Δx = 1 Δx = 1 Δx n= n= n= δ k x n Δx G(k x ) δ k x n Δx G k x n Δx Fourier Transform of g S (x) G(k x ) k x 1/Δx G S (k x ) k x 1/Δx 18
19 Nyquist Condition G(k x ) K S -B B k x G S (k x ) k x K S =1/Δx To avoid overlap, we require that 1/Δx>2B or K S > 2B where K S =1/ Δx is the sampling frequency Example Assume that the highest spatial frequency in an object is B = 2 cm -1. Thus, smallest spatial period is 0.5 cm.. Nyquist theorem says we need to sample with Δx < 1/2B = 0.25 cm This corresponds to 2 samples per spatial period. 19
20 Reconstruction from Samples K S G S (k x ) K S =1/Δx Multiply by (1/K S )rect(k x /K S ) (1/K S ) G S (k x )rect(k x /K S ) =G(k x ) Example Cosine Reconstruction cos(2πk 0 x) -k 0 k 0 K S K S >2k 0 -k 0 k 0 KS K S =2k 0 -k 0 k 0 KS 20
21 Reconstruction from Samples If the Nyquist condition is met, then G ˆ S (k x ) = 1 G S (k x )rect(k x /K S ) = G(k x ) K S And the signal can be reconstructed by convolving the sample with a sinc function g ˆ S (x) = g S (x) sinc(k s x) = g(nδx)δ( x nδx) sinc(k s x) n= = g(nδx) sinc(k s (x nδx)) n= g(x) Reconstruction from Samples g S (x) Sample at Δx x ˆ g S (x) sinc(k s x) = sinc(x / Δx) 21
22 Cosine Example with K S =2k 0 Example with K s =4k 0 22
23 Example with K s =8k 0 Aliasing Example 23
24 Aliasing -B B G(k x ) k x K S Aliasing occurs when the Nyquist condition is not satisfied. This occurs for K S 2B Aliasing Example cos(2πk 0 x) -k 0 k 0 K S =k 0 -k 0 k 0 K S 24
25 Aliasing Example cos(2πk 0 x) -k 0 k 0 2k 0 >K S >k 0 -k 0 k 0 K S Practical Considerations Why sample higher than the Nyquist frequency? -- true sinc interpolation is not practical since the sinc function goes from -infinity to infinity -- the requirements on the low-pass filter are reduced. Hard -k 0 k 0 Easier -k 0 k 0 K S 25
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 informationTopics. 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 informationFourier 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 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 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 informationSpin-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 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 informationBioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2004 Lecture 2 Linear Systems. Topics
Bioegieerig 280A Priciples of Biomedical Imagig Fall Quarter 2004 Lecture 2 Liear Systems Topics 1. Liearity 2. Impulse Respose ad Delta fuctios 3. Superpositio Itegral 4. Shift Ivariace 5. 1D ad 2D covolutio
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 informationBioengineering 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 informationTopics. 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 informationAssignment 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 informationFourier 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 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 informationGBS765 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 informationLinear Operators and Fourier Transform
Linear Operators and 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
More informationDiscrete 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 informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Spring 2014 TTh 14:30-15:45 CBC C313 Lecture 05 Image Processing Basics 13/02/04 http://www.ee.unlv.edu/~b1morris/ecg782/
More informationBME 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 informationConvolution and Linear Systems
CS 450: Introduction to Digital Signal and Image Processing Bryan Morse BYU Computer Science Introduction Analyzing Systems Goal: analyze a device that turns one signal into another. Notation: f (t) g(t)
More informationDiscrete 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 informationG52IVG, 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 informationImages 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 informationFourier Transform in Image Processing. CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012)
Fourier Transform in Image Processing CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012) Basis Decomposition Write a function as a weighted sum of basis functions f ( x) wibi(
More informationFourier transform. Stefano Ferrari. Università degli Studi di Milano Methods for Image Processing. academic year
Fourier transform Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Methods for Image Processing academic year 27 28 Function transforms Sometimes, operating on a class of functions
More informationBiophysics 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 informationImage preprocessing in spatial domain
Image preprocessing in spatial domain Sampling theorem, aliasing, interpolation, geometrical transformations Revision: 1.4, dated: May 25, 2006 Tomáš Svoboda Czech Technical University, Faculty of Electrical
More informationImage preprocessing in spatial domain
Image preprocessing in spatial domain Sampling theorem, aliasing, interpolation, geometrical transformations Revision: 1.3, dated: December 7, 2005 Tomáš Svoboda Czech Technical University, Faculty of
More informationModule 3 : Sampling and Reconstruction Lecture 22 : Sampling and Reconstruction of Band-Limited Signals
Module 3 : Sampling and Reconstruction Lecture 22 : Sampling and Reconstruction of Band-Limited Signals Objectives Scope of this lecture: If a Continuous Time (C.T.) signal is to be uniquely represented
More informationIB Paper 6: Signal and Data Analysis
IB Paper 6: Signal and Data Analysis Handout 5: Sampling Theory S Godsill Signal Processing and Communications Group, Engineering Department, Cambridge, UK Lent 2015 1 / 85 Sampling and Aliasing All of
More informationDISCRETE 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 informationFrequency2: Sampling and Aliasing
CS 4495 Computer Vision Frequency2: Sampling and Aliasing Aaron Bobick School of Interactive Computing Administrivia Project 1 is due tonight. Submit what you have at the deadline. Next problem set stereo
More informationBackprojection. Projections. Projections " $ & cosθ & Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2014 CT/Fourier Lecture 2
Backprojection Bioengineering 280A Principles of Biomedical Imaging 3 0 0 0 0 3 0 0 0 0 0 3 0 0 3 0 3 Fall Quarter 2014 CT/Fourier Lecture 2 0 0 0 1 1 1 0 0 0 1 0 0 1 2 1 0 0 1 1 1 0 1 3 1 0 1 1 1 1 1
More informationReview: Continuous Fourier Transform
Review: Continuous Fourier Transform Review: convolution x t h t = x τ h(t τ)dτ Convolution in time domain Derivation Convolution Property Interchange the order of integrals Let Convolution Property By
More informationESS Finite Impulse Response Filters and the Z-transform
9. Finite Impulse Response Filters and the Z-transform We are going to have two lectures on filters you can find much more material in Bob Crosson s notes. In the first lecture we will focus on some of
More informationQuality Improves with More Rays
Recap Quality Improves with More Rays Area Area 1 shadow ray 16 shadow rays CS348b Lecture 8 Pat Hanrahan / Matt Pharr, Spring 2018 pixelsamples = 1 jaggies pixelsamples = 16 anti-aliased Sampling and
More informationFourier 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 informationEE 224 Signals and Systems I Review 1/10
EE 224 Signals and Systems I Review 1/10 Class Contents Signals and Systems Continuous-Time and Discrete-Time Time-Domain and Frequency Domain (all these dimensions are tightly coupled) SIGNALS SYSTEMS
More informationComputer Vision & Digital Image Processing
Computer Vision & Digital Image Processing Image Restoration and Reconstruction I Dr. D. J. Jackson Lecture 11-1 Image restoration Restoration is an objective process that attempts to recover an image
More informationImage Processing and Computer Vision
Image Processing and Computer Vision Processing of continuous images Image Processing and Computer Vision linear filtering Fourier transformation Wiener filtering Nonlinear diffusion Visual Computing:
More information1-D MATH REVIEW CONTINUOUS 1-D FUNCTIONS. Kronecker delta function and its relatives. x 0 = 0
-D MATH REVIEW CONTINUOUS -D FUNCTIONS Kronecker delta function and its relatives delta function δ ( 0 ) 0 = 0 NOTE: The delta function s amplitude is infinite and its area is. The amplitude is shown as
More informationChap 4. Sampling of Continuous-Time Signals
Digital Signal Processing Chap 4. Sampling of Continuous-Time Signals Chang-Su Kim Digital Processing of Continuous-Time Signals Digital processing of a CT signal involves three basic steps 1. Conversion
More informationAdditional Pointers. Introduction to Computer Vision. Convolution. Area operations: Linear filtering
Additional Pointers Introduction to Computer Vision CS / ECE 181B andout #4 : Available this afternoon Midterm: May 6, 2004 W #2 due tomorrow Ack: Prof. Matthew Turk for the lecture slides. See my ECE
More informationOptical Transfer Function
Optical Transfer Function y x y' x' electro-optical system input I(x,y) output image U(x,y) The OTF describes resolution of an image-system. It parallels the concept of electrical transfer function, about
More information6 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 informationThe Poisson summation formula, the sampling theorem, and Dirac combs
The Poisson summation ormula, the sampling theorem, and Dirac combs Jordan Bell jordanbell@gmailcom Department o Mathematics, University o Toronto April 3, 24 Poisson summation ormula et S be the set o
More informationBiomedical 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 informationTherefore 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 informationFourier Series Example
Fourier Series Example Let us compute the Fourier series for the function on the interval [ π,π]. f(x) = x f is an odd function, so the a n are zero, and thus the Fourier series will be of the form f(x)
More informationSignal Processing and Linear Systems1 Lecture 4: Characterizing Systems
Signal Processing and Linear Systems Lecture : Characterizing Systems Nicholas Dwork www.stanford.edu/~ndwork Our goal will be to develop a way to learn how the system behaves. In general, this is a very
More informationESE 531: Digital Signal Processing
ESE 531: Digital Signal Processing Lec 9: February 13th, 2018 Downsampling/Upsampling and Practical Interpolation Lecture Outline! CT processing of DT signals! Downsampling! Upsampling 2 Continuous-Time
More informationImage Processing and Computer Vision. Visual Computing: Joachim M. Buhmann 1/66
Image Processing and Computer Vision Visual Computing: Joachim M. Buhmann 1/66 Image Processing and Computer Vision Processing of continuous images linear filtering Fourier transformation Wiener filtering
More informationThe 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 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 informationModule 1: Signals & System
Module 1: Signals & System Lecture 6: Basic Signals in Detail Basic Signals in detail We now introduce formally some of the basic signals namely 1) The Unit Impulse function. 2) The Unit Step function
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 information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 2.161 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Massachusetts
More information! Circular Convolution. " Linear convolution with circular convolution. ! Discrete Fourier Transform. " Linear convolution through circular
Previously ESE 531: Digital Signal Processing Lec 22: April 18, 2017 Fast Fourier Transform (con t)! Circular Convolution " Linear convolution with circular convolution! Discrete Fourier Transform " Linear
More informationDiscrete Fourier transform (DFT)
Discrete Fourier transform (DFT) Signal Processing 2008/9 LEA Instituto Superior Técnico Signal Processing LEA (IST) Discrete Fourier transform 1 / 34 Periodic signals Consider a periodic signal x[n] with
More informationSolutions of differential equations using transforms
Solutions of differential equations using transforms Process: Take transform of equation and boundary/initial conditions in one variable. Derivatives are turned into multiplication operators. Solve (hopefully
More informationNov : Lecture 18: The Fourier Transform and its Interpretations
3 Nov. 04 2005: Lecture 8: The Fourier Transform and its Interpretations Reading: Kreyszig Sections: 0.5 (pp:547 49), 0.8 (pp:557 63), 0.9 (pp:564 68), 0.0 (pp:569 75) Fourier Transforms Expansion of a
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
More informationSEISMIC WAVE PROPAGATION. Lecture 2: Fourier Analysis
SEISMIC WAVE PROPAGATION Lecture 2: Fourier Analysis Fourier Series & Fourier Transforms Fourier Series Review of trigonometric identities Analysing the square wave Fourier Transform Transforms of some
More informationECE 425. Image Science and Engineering
ECE 425 Image Science and Engineering Spring Semester 2000 Course Notes Robert A. Schowengerdt schowengerdt@ece.arizona.edu (520) 62-2706 (voice), (520) 62-8076 (fa) ECE402 DEFINITIONS 2 Image science
More informationFourier Transform Chapter 10 Sampling and Series
Fourier Transform Chapter 0 Sampling and Series Sampling Theorem Sampling Theorem states that, under a certain condition, it is in fact possible to recover with full accuracy the values intervening between
More informationFourier Transform for Continuous Functions
Fourier Transform for Continuous Functions Central goal: representing a signal by a set of orthogonal bases that are corresponding to frequencies or spectrum. Fourier series allows to find the spectrum
More informationSignals & Systems. Chapter 7: Sampling. Adapted from: Lecture notes from MIT, Binghamton University, and Purdue. Dr. Hamid R.
Signals & Systems Chapter 7: Sampling Adapted from: Lecture notes from MIT, Binghamton University, and Purdue Dr. Hamid R. Rabiee Fall 2013 Outline 1. The Concept and Representation of Periodic Sampling
More informationCONTINUOUS IMAGE MATHEMATICAL CHARACTERIZATION
c01.fm Page 3 Friday, December 8, 2006 10:08 AM 1 CONTINUOUS IMAGE MATHEMATICAL CHARACTERIZATION In the design and analysis of image processing systems, it is convenient and often necessary mathematically
More informationCS-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 informationAnalog vs. discrete signals
Analog vs. discrete signals Continuous-time signals are also known as analog signals because their amplitude is analogous (i.e., proportional) to the physical quantity they represent. Discrete-time signals
More informationLecture Schedule: Week Date Lecture Title
http://elec34.org Sampling and CONVOLUTION 24 School of Information Technology and Electrical Engineering at The University of Queensland Lecture Schedule: Week Date Lecture Title 2-Mar Introduction 3-Mar
More informationChap 2. Discrete-Time Signals and Systems
Digital Signal Processing Chap 2. Discrete-Time Signals and Systems Chang-Su Kim Discrete-Time Signals CT Signal DT Signal Representation 0 4 1 1 1 2 3 Functional representation 1, n 1,3 x[ n] 4, n 2 0,
More informationIntroduction to Linear Systems
cfl David J Fleet, 998 Introduction to Linear Systems David Fleet For operator T, input I, and response R = T [I], T satisfies: ffl homogeniety: iff T [ai] = at[i] 8a 2 C ffl additivity: iff T [I + I 2
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 information3. Lecture. Fourier Transformation Sampling
3. Lecture Fourier Transformation Sampling Some slides taken from Digital Image Processing: An Algorithmic Introduction using Java, Wilhelm Burger and Mark James Burge Separability ² The 2D DFT can be
More informationComputer Vision & Digital Image Processing. Periodicity of the Fourier transform
Computer Vision & Digital Image Processing Fourier Transform Properties, the Laplacian, Convolution and Correlation Dr. D. J. Jackson Lecture 9- Periodicity of the Fourier transform The discrete Fourier
More informationSpatial-Domain Convolution Filters
Spatial-Domain Filtering 9 Spatial-Domain Convolution Filters Consider a linear space-invariant (LSI) system as shown: The two separate inputs to the LSI system, x 1 (m) and x 2 (m), and their corresponding
More informationMultidimensional digital signal processing
PSfrag replacements Two-dimensional discrete signals N 1 A 2-D discrete signal (also N called a sequence or array) is a function 2 defined over thex(n set 1 of, n 2 ordered ) pairs of integers: y(nx 1,
More informationLecture 3: Linear Filters
Lecture 3: Linear Filters Professor Fei Fei Li Stanford Vision Lab 1 What we will learn today? Images as functions Linear systems (filters) Convolution and correlation Discrete Fourier Transform (DFT)
More informationSampling in 1D ( ) Continuous time signal f(t) Discrete time signal. f(t) comb
Sampling in 2D 1 Sampling in 1D Continuous time signal f(t) Discrete time signal t ( ) f [ k] = f( kt ) = f( t) δ t kt s k s f(t) comb k 2 Nyquist theorem (1D) At least 2 sample/period are needed to represent
More informationThe Convolution Sum for Discrete-Time LTI Systems
The Convolution Sum for Discrete-Time LTI Systems Andrew W. H. House 01 June 004 1 The Basics of the Convolution Sum Consider a DT LTI system, L. x(n) L y(n) DT convolution is based on an earlier result
More informationRepresentation of Signals and Systems. Lecturer: David Shiung
Representation of Signals and Systems Lecturer: David Shiung 1 Abstract (1/2) Fourier analysis Properties of the Fourier transform Dirac delta function Fourier transform of periodic signals Fourier-transform
More informationSummary of Lecture 5
Summary of Lecture 5 In lecture 5 we learnt how to pick the reproduction levels for the given thresholds. We learnt how to design MSQE optimal (Lloyd-Max) quantizers. We reviewed linear systems, linear
More informationContents. Signals as functions (1D, 2D)
Fourier Transform The idea A signal can be interpreted as en electromagnetic wave. This consists of lights of different color, or frequency, that can be split apart usign an optic prism. Each component
More informationModule 3. Convolution. Aim
Module Convolution Digital Signal Processing. Slide 4. Aim How to perform convolution in real-time systems efficiently? Is convolution in time domain equivalent to multiplication of the transformed sequence?
More informationLecture 3: Linear Filters
Lecture 3: Linear Filters Professor Fei Fei Li Stanford Vision Lab 1 What we will learn today? Images as functions Linear systems (filters) Convolution and correlation Discrete Fourier Transform (DFT)
More informationSignal 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 informationFourier Transform 4: z-transform (part 2) & Introduction to 2D Fourier Analysis
052600 VU Signal and Image Processing Fourier Transform 4: z-transform (part 2) & Introduction to 2D Fourier Analysis Torsten Möller + Hrvoje Bogunović + Raphael Sahann torsten.moeller@univie.ac.at hrvoje.bogunovic@meduniwien.ac.at
More informationThe Fourier Transform
The Fourier Transorm Fourier Series Fourier Transorm The Basic Theorems and Applications Sampling Bracewell, R. The Fourier Transorm and Its Applications, 3rd ed. New York: McGraw-Hill, 2. Eric W. Weisstein.
More informationTutorial Sheet #2 discrete vs. continuous functions, periodicity, sampling
2.39 utorial Sheet #2 discrete vs. continuous functions, periodicity, sampling We will encounter two classes of signals in this class, continuous-signals and discrete-signals. he distinct mathematical
More informationAssignment 2 Solutions
EE225E/BIOE265 Spring 2012 Principles of MRI Miki Lustig Assignment 2 Solutions 1. Read Nishimura Ch. 3 2. Non-Uniform Sampling. A student has an assignment to monitor the level of Hetch-Hetchi reservoir
More informationBioE 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 informationInformation and Communications Security: Encryption and Information Hiding
Short Course on Information and Communications Security: Encryption and Information Hiding Tuesday, 10 March Friday, 13 March, 2015 Lecture 5: Signal Analysis Contents The complex exponential The complex
More informationWave Phenomena Physics 15c. Lecture 11 Dispersion
Wave Phenomena Physics 15c Lecture 11 Dispersion What We Did Last Time Defined Fourier transform f (t) = F(ω)e iωt dω F(ω) = 1 2π f(t) and F(w) represent a function in time and frequency domains Analyzed
More information1 otherwise. Note that the area of the pulse is one. The Dirac delta function (a.k.a. the impulse) can be defined using the pulse as follows:
The Dirac delta function There is a function called the pulse: { if t > Π(t) = 2 otherwise. Note that the area of the pulse is one. The Dirac delta function (a.k.a. the impulse) can be defined using the
More informationInterchange of Filtering and Downsampling/Upsampling
Interchange of Filtering and Downsampling/Upsampling Downsampling and upsampling are linear systems, but not LTI systems. They cannot be implemented by difference equations, and so we cannot apply z-transform
More informationFrequency-Domain C/S of LTI Systems
Frequency-Domain C/S of LTI Systems x(n) LTI y(n) LTI: Linear Time-Invariant system h(n), the impulse response of an LTI systems describes the time domain c/s. H(ω), the frequency response describes the
More informationESE 531: Digital Signal Processing
ESE 531: Digital Signal Processing Lec 8: February 7th, 2017 Sampling and Reconstruction Lecture Outline! Review " Ideal sampling " Frequency response of sampled signal " Reconstruction " Anti-aliasing
More informationECE 301 Fall 2011 Division 1 Homework 10 Solutions. { 1, for 0.5 t 0.5 x(t) = 0, for 0.5 < t 1
ECE 3 Fall Division Homework Solutions Problem. Reconstruction of a continuous-time signal from its samples. Let x be a periodic continuous-time signal with period, such that {, for.5 t.5 x(t) =, for.5
More informationReview of Analog Signal Analysis
Review of Analog Signal Analysis Chapter Intended Learning Outcomes: (i) Review of Fourier series which is used to analyze continuous-time periodic signals (ii) Review of Fourier transform which is used
More information