Assignment 2 Solutions
|
|
- Derick Poole
- 5 years ago
- Views:
Transcription
1 EE225E/BIOE265 Spring 2012 Principles of MRI Miki Lustig Assignment 2 Solutions 1. Read Nishimura Ch Non-Uniform Sampling. A student has an assignment to monitor the level of Hetch-Hetchi reservoir in Yosemite National Park. In order to meet the Nyquist rate, he has to drive and measure the water level at least once a day. Driving every day for 3 hours one way is no fun, so he came up with a new plan. In the new plan he visits the reservoir every 2 days. He arrives at yosemite early in the day and records the level of the reservoir. Then, he goes on a 6 hours hike and before heading back he records the level again. The samples he collects are g 1 [n] = g(nt ) and g 2 [n] = g(nt + 1/4), where T = 2days and g(t) is the continuous water level of the lake. Therefore, ĝ 1 (t) = g(t) 1 T X( t T ) and, ĝ 2 (t) = g(t) 1 T X( t T T ) where T = 1/4 We will show in this question that it is possible to reconstruct the continuous signal, even though the Nyquist rate appears to be violated. However, there are some costs associated with this. We will solve the problem in the frequency domain. Solution: The basic setup in this question is that instead of uniform sampling at the Nyquist rate we have non-uniform sampling, which on-average samples at the required rate (one sample too close, one too far). T ΔT T (a) g(t) is band limited. If T = 1 [days] is the Nyquist rate, what is the highest frequency fnyq? Solution: The highest frequency is (1/T)/2 which is 1/2[day 1 ]. Since we sample at a rate T = 2 [days], both Ĝ1(f) and Ĝ2(f) are periodic with a period 1/2 [day 1 ]. We would like to express recover G(f) from Ĝ1(f) and Ĝ2(f). Since G(f) is band limited between 1/2 < f < 1/2, we need only focus on that interval. 1
2 (b) For the interval 0 f 1/2[days 1 ] express Ĝ1(f) as a function of G(f). Solution: Computing the Fourier transform of ĝ 1 (t) yields: Ĝ 1 (f) = G(f) X(T f) So, G(f) is replicated every 0.5 interval. From the above figure we see that between 0 and 0.5 there are only two replicas. Therefore, Ĝ 1 (f) = G(f) + G(f 1/2) 0 f 1/2 (c) For the interval 0 f 1/2 [days 1 ] express Ĝ2(f) as a function of G(f). Solution: Computing the Fourier transform of ĝ 2 (t) yields: ( ) Ĝ 1 (f) = G(f) e i2πt T f X(T f) In a similar way we get a replication every 0.5 interval, but now each replica will have a different phase. For T = 2, T = 1/4, f = 1/2 we get that the phase of the replica centered at f = 1/2 is e i2πt T f = e iπ/2 = i. From the above figure we see that between 0 and 0.5 there are only two replicas. Therefore, Ĝ 1 (f) = G(f) ig(f 1/2) 0 f 1/2 (d) Use the result from (b)+(c) to find G(f) as a function of Ĝ 1 (f) and Ĝ2(f) for 0 f 1/2 [days 1 ]. 2
3 Solution: We can express the two equations in matrix form: [ ] [ 1 1 G(f) 1 i G(f 1/2) We then get that, [ and G(f) G(f 1/2) ] G(f) = iĝ1(f) + Ĝ2(f) i + 1 ] = = 1 [ i 1 i (e) Repeat (b)-(d) for the interval 1/2 f 0 [days 1 ]. Solution: In a similar way, G(f) = iĝ1(f) + Ĝ2(f) i + 1 [ Ĝ1 (f) Ĝ 2 (f) ] ] [ Ĝ1 (f) Ĝ 2 (f) ] 0 f 1/2 1/2 f 0 (f) In order to recover G(f) you had to invert a set of linear equations. Express the set of equations in (b)-(c) in matrix form. Plot the condition number (use the cond function in Matlab) for different shifts, e.g., 0 < T < 1/2. The condition number is the ratio between the highest and lowest eigen-values and indicates how much noise amplification occurs when inverting the matrix. From the plot, what can you make a conclusion on the cost of non-uniform sampling? Solutions: For T = 1/2 we have uniform sampling at the Nyquist rate and therefore a condition number equal to 1. However, as T is smaller and smaller our condition number worsens and the penalty is sensitivity to noise Condition number vs T Condition number T 3
4 3. Interpolation in k-space: In this problem we will see the effect of linear interpolation in k-space. This is a more involved version of the question that was given in the midterm in The following image was resampled in the vertical direction in k-space by linear interpolation at 3 4 the original rate. Interpolation in k-space corresponds to finer k-space sampling and therefor a larger FOV. So, the FOV of the interpolated image (right) is larger in the vertical dimension. There are two distinct artifacts that appear in the interpolated image and not in the original. The interpolated image is apodized (no longer has flat magnitude). In addition, there are some aliasing artifacts that appear as replication of the image in the interpolated dimension. These artifacts are the result of using a linear interpolation as opposed to a better interpolation scheme. In the following questions we will learn the source of these artifacts through a simple 1D example. You are given the following 1D signal that you want to sample in k-space. In this question, f(x) is the space signal and F (k) is its k-space. f(x) x, [cm] a) You decide to collect 256 k-space samples of F (k) with k = cycles/cm. Express ˆf(x), the result of sampling F (k) in k-space, as a function of f(x). Draw the signal ˆf(x). Emphasize any artifacts in the drawing and point out the differences with the original f(x). Solution: There are going to be two effects. Blurr due to the finite extent in k-space and replication due to the discrete sampling. Since k = and the number of samples is 256 the total extent in k-space is W kz = = 8 cycles/cm. The finite extent can therefore be written 4
5 as a multiplication with the rect function (k/8). The discrete sampling at rate of k = cycles/cm can be written as a multiplication with the Shah function 8 X( 8 k). So, we can write ˆF (k) as, ˆF (k) = F (k) X(256 8 k) (k 8 ). Taking the inverse Fourier transform yields, ˆf(x) = f(x) 8sinc(8x) X( 8 }{{} 256 x) blurr }{{} replication = 8 (f(x) 8sinc(8x)δ(x 32n) 256 n= ^ f(x) blurr ringing 1/ x, [cm] b) You prepare to show the results from part (a) to your adviser, however, you suddenly remember that he asked you to sample k-space at a finer interval of k = cycles/cm. Not wanting to repeat the experiment, you decide to interpolate ˆF (k) and resample it accordingly. Write an expression for F (k), the resampled signal, assuming linear interpolation. Express F (k) as a function of ˆF (k), HINT: { linear interpolation is a convolution with an appropriately scaled triangle function Λ(ak) = 1 ak ak < 1 0 otherwise Solution: 8 Linear interpolation is a convolution with a triangle function. Our samples are spaced at 256 cycles/cm interval, so the convolution is with Λ( k). The result Λ( 8 k) ˆF (k) is a continuous function. In order to resample, we need to multiply the continuous function with an impulse 6 train at 256 cycles/cm interval. Therefore we get that, ( F (k) = Λ( 256 ) 8 k) ˆF 256 (k) 6 X(256 6 k) c) Write an expression for f(x), the image domain representation of F (k) from part (b). First express f(x) as a function of ˆf(x) and then as a function of f(x). Draw f(x), pointing out any artifacts due to the interpolation and resampling. (HINT: First try to draw the effect of the interpolation, then draw the effect of the resampling.) Solution: For this question we need to know that the Fourier transform of a triangle function is, F 1 {Λ(k)} = sinc 2 (x). So, F 1 {Λ(ak)} = 1 a sinc2 (x/a). This can also be derived from the relation Λ(k) = (k) (k). 5
6 So, the convolution and resampling yields, ( 8 f(x) = ˆf(x) 256 sinc2 ( 8 ) 256 x) X( x). and, f(x) = Drawinf ˆf(x) sinc2 ( 8 256x) we get, ~ f(x) (( f(x) 8sinc(8x) X( 8 ) 256 x) sinc2 ( 8 ) 256 x) X( x) x, [cm] in which the replicas have not been fully suppressed and the signal heavily apodized. Convolving the result with X( 6 256x) we get, x, [cm] ~ f(x) aliasing x, [cm] d) Linear interpolation can be considered as interpolating with a poor low-pass filter. How would the result change if we use a filter with 0.1% ripple in the pass-band and stop band for the interpolation? Draw the result pointing out the differences. Solutions: A TBW=8 windowed sinc is a better filter and suppresses the replicas much better as illustrated below: 6
7 ~ f(x) slight apodization x, [cm] The replication of the signal will have negligible aliasing within the field of view! ~ f(x) almost no aliasing x, [cm] 4. Matlab assignment: 2D Fourier transforms: Use the function publish to make your report. In this assignment we will play with 2D Fourier transforms. As mentioned in class, in MRI the (0,0) coordinates both in k-space and of the image are located in the center. The functions fft2 and ifft2 in Matlab place the origin of the axis at the corner of the image. We will therefore need to write new functions that compute the centered DFT. (a) Program the functions F=fft2c(f) and f=ifft2c(f). Recall that >> F = ifftshift(fft2(fftshift(f))); >> f = fftshift(ifft2(ifftshift(f))); (b) Rectangles in k-space. Download the functions mysinc.m and myjinc.m. The functions implement the sinc and jinc function as described in class. They accepts a vector (or matrix) of inputs and return the results in the same format as the input. The Fourier transform of a rectangle with length 18[cm] and amplitude 1 is F (k) = 18 2 sinc(18k x )sinc(18k y ) We would like to sample the rectangle in k-space while achieving a resolution of 1mm and a 25.6cm Field of view. Recall from class that the resolution x = 1 W kx where W kx is the extent in k-space. Also recall that FOV x = 1 k x where k x is the spacing between samples in k-space. For 1mm, we get that W kx = 10 and k x = 10/256. 7
8 To sample the function we will create a grid using the function meshgrid. First, let s create the intervals: >> Wkx = 10; >> Wky = 10; >> Nx = 256; >> Ny = 256; >> kx_itvl = [ -Nx/2:Nx/2-1]/Nx * Wkx; >> ky_itvl = [ -Ny/2:Ny/2-1]/Ny * Wky; Now, we will create the grid and evaluate the function >> [kx, ky] = meshgrid(kx_itvl, ky_itvl); >> F = 18*18*mysinc(18*kx).*mysinc(18*ky); Since the sample lie on a grid, we can reconstruct the data using the inverse fft. In the inverse fft, the data is scaled by 1/(NxNy) = 1/(256 2 ). On the other hand, each point in k-space represents an area of k x k y = 10 2 / So, for the ifft to approximate the inverse Fourier transform we need to multiply the result by a 100. This is however not very important in MRI. The actual scaling of an image in MRI depends on many things and can vary significantly between scans. The contrast in the image is the important information. We will do it here just so the rectangle will have an amplitude of 1. But again, this is not so important. >> f = 100*ifft2c(F); >> figure(1), imshow(abs(f),[]), title( f ) >> figure(2), imshow(angle(f),[]), title( The phase of f ); Validate that the size of the rectangle is 180 pixels. In a similar way we can sample a circle in k-space. Here s an example of a circle with a diameter of 20cm: >> r = sqrt(kx.^2+ky.^2); >> F_circ = 20^2*myjinc(20*r); >> f_circ = 100*ifft2c(F_circ); >> figure(1), imshow(abs(f_circ),[]), title( f ) >> figure(2), imshow(angle(f_circ),[]), title( The phase of f ); To demonstrate the modulation property of the Fourier transform we will shift the k-space of the rectangle by 2 k x. We should expect 2 cycles of phase across the FOV in the x-axis. >> F_srect = 22*22*mysinc(22*(kx-10/256*2)).*mysinc(22*ky); >> f_srect = 100*ifft2c(F_srect); >> figure(1), imshow(abs(f_srect),[]), title( f ) >> figure(2), imshow(angle(f_srect),[]), title( The phase of f ); To shift a rectangle in the image domain, we need to modulate k-space accordingly. Here s an example of shifting a rectangle of length 5cm by 5cm in both axis. 8
9 >> F_shift = 10*10*mysinc(10*kx).*mysinc(10*ky).*exp(-i*2*pi*(5*kx+5*ky)); >> f_shift = 100*ifft2c(F_shift); >> figure(1), imshow(abs(f_shift),[]), title( f ) >> figure(2), imshow(angle(f_shift),[]), title( The phase of f ); (c) Write the function F = rectf(kx,ky,lx,ly,x0,y0). The function evaluates the Fourier transform of a rectangle. kx and ky are k-space coordinates matrices (can be generated by meshgrid for example). Lx and Ly are scalars corresponding to the length of the rectangle in the x and y axis. x0 and y0 are scalars corresponding to the coordinates of the center of the rectangle. Write the function F = circf(kx,ky,diam,x0,y0). The function evaluates the Fourier transform of a circle. kx and ky are k-space coordinates matrices (can be generated by meshgrid for example). diam is a scaler representing the diameter of the circle. x0 and y0 are scalars corresponding to the coordinates of the center of the circle. Using these functions create the function F = face(kx,ky). The function evaluates the Fourier transform of the following image: 1.5cm 3.5cm 1cm 1cm 20cm 4cm 1.5cm 2.5cm 10cm 3.5cm Intensity=0 Intensity=0.4 Intensity=0.8 (d) Sample and display the resulting image for the following parameters: i. FOV = cm, resolution = 4mm ii. FOV = 24 cm, resolution = 5mm iii. FOV = 25.6 cm, resolution = 1mm iv. FOV = 19.2 cm, resolution = 1mm What happens when the imaging FOV is smaller than the object? Solutions: The solutions are based on Frank Ong s hw solutions. Thanks Frank! 9
10 Table of Contents EE225e hw b... 1 Rectangle... 1 Circle... 3 Shifted k-space... 5 Modulated k-space d Sample and Display face d i) FOV = 102.4cm, res = 4mm d ii) FOV = 24cm, res = 5mm d iii) FOV = 25.6cm, res = 1mm d iv) FOV = 19.2cm, res = 1mm Aliasing happens when FOV is smaller than object EE225e hw2-3 1b close all clc clear Wkx = 10; Wky = 10; Nx = 256; Ny = 256; kx_itvl = [ -Nx/2:Nx/2-1]/Nx * Wkx; ky_itvl = [ -Ny/2:Ny/2-1]/Ny * Wky; [kx, ky] = meshgrid(kx_itvl, ky_itvl); F = 18*18*mysinc(18*kx).*mysinc(18*ky); Rectangle f = 100*ifft2c(F); figure(1), imshow(abs(f),[]), title(' f ') figure(2), imshow(angle(f),[]), title('the phase of f'); 1
11 2
12 Circle r = sqrt(kx.^2+ky.^2); F_circ = 20^2*myjinc(20*r); f_circ = 100*ifft2c(F_circ); figure(1), imshow(abs(f_circ),[]), title(' f ') figure(2), imshow(angle(f_circ),[]), title('the phase of f'); 3
13 4
14 Shifted k-space F_srect = 22*22*mysinc(22*(kx-10/256*2)).*mysinc(22*ky); f_srect = 100*ifft2c(F_srect); figure(1), imshow(abs(f_srect),[]), title(' f ') figure(2), imshow(angle(f_srect),[]), title('the phase of f'); 5
15 6
16 Modulated k-space F_shift = 10*10*mysinc(10*kx).*mysinc(10*ky).*exp(-i*2*pi*(5*kx+5*ky)); f_shift = 100*ifft2c(F_shift); figure(1), imshow(abs(f_shift),[]), title(' f ') figure(2), imshow(angle(f_shift),[]), title('the phase of f'); 7
17 8
18 1d Sample and Display face 1d i) FOV = 102.4cm, res = 4mm Wkx = 10/4; Wky = 10/4; Nx = 1024/4; Ny = 1024/4; kx_itvl = [ -Nx/2:Nx/2-1]/Nx * Wkx; ky_itvl = [ -Ny/2:Ny/2-1]/Ny * Wky; [kx, ky] = meshgrid(kx_itvl, ky_itvl); F = face(kx, ky); f = Wkx*Wky*ifft2c(F); figure(1), imshow(abs(f),[0 1]), title('face 1') 9
19 1d ii) FOV = 24cm, res = 5mm Wkx = 10/5; Wky = 10/5; Nx = 240/5; Ny = 240/5; kx_itvl = [ -Nx/2:Nx/2-1]/Nx * Wkx; ky_itvl = [ -Ny/2:Ny/2-1]/Ny * Wky; [kx, ky] = meshgrid(kx_itvl, ky_itvl); F = face(kx, ky); f = Wkx*Wky*ifft2c(F); figure(2), imshow(abs(f),[0 1]), title('face 2') 10
20 1d iii) FOV = 25.6cm, res = 1mm Wkx = 10; Wky = 10; Nx = 256; Ny = 256; kx_itvl = [ -Nx/2:Nx/2-1]/Nx * Wkx; ky_itvl = [ -Ny/2:Ny/2-1]/Ny * Wky; [kx, ky] = meshgrid(kx_itvl, ky_itvl); F = face(kx, ky); f = Wkx*Wky*ifft2c(F); figure(3), imshow(abs(f),[0 1]), title('face 3') 1d iv) FOV = 19.2cm, res = 1mm Wkx = 10; Wky = 10; Nx = 192; Ny = 192; kx_itvl = [ -Nx/2:Nx/2-1]/Nx * Wkx; ky_itvl = [ -Ny/2:Ny/2-1]/Ny * Wky; [kx, ky] = meshgrid(kx_itvl, ky_itvl); F = face(kx, ky); f = Wkx*Wky*ifft2c(F); figure(4), imshow(abs(f),[0 1]), title('face 4') 11
21 Aliasing happens when FOV is smaller than object. Published with MATLAB
Lecture 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 informationEE225E/BIOE265 Spring 2013 Principles of MRI. Assignment 9 Solutions. Due April 29th, 2013
EE5E/BIOE65 Spring 013 Principles of MRI Miki Lustig This is the last homework in class. Enjoy it. Assignment 9 Solutions Due April 9th, 013 1) In class when we presented the spin-echo saturation recovery
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 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 informationE2.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 informationTopic 7. Convolution, Filters, Correlation, Representation. Bryan Pardo, 2008, Northwestern University EECS 352: Machine Perception of Music and Audio
Topic 7 Convolution, Filters, Correlation, Representation Short time Fourier Transform Break signal into windows Calculate DFT of each window The Spectrogram spectrogram(y,1024,512,1024,fs,'yaxis'); A
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 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 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 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 informationName: INSERT YOUR NAME HERE. Due to dropbox by 6pm PDT, Wednesday, December 14, 2011
AMath 584 Name: INSERT YOUR NAME HERE Take-home Final UWNetID: INSERT YOUR NETID Due to dropbox by 6pm PDT, Wednesday, December 14, 2011 The main part of the assignment (Problems 1 3) is worth 80 points.
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 informationECS 332: Principles of Communications 2012/1. HW 4 Due: Sep 7
ECS 332: Principles of Communications 2012/1 HW 4 Due: Sep 7 Lecturer: Prapun Suksompong, Ph.D. Instructions (a) ONE part of a question will be graded (5 pt). Of course, you do not know which part will
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 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 informationTwo-Dimensional Signal Processing and Image De-noising
Two-Dimensional Signal Processing and Image De-noising Alec Koppel, Mark Eisen, Alejandro Ribeiro March 14, 2016 Before we considered (one-dimensional) discrete signals of the form x : [0, N 1] C where
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 informationNumerical Approximation Methods for Non-Uniform Fourier Data
Numerical Approximation Methods for Non-Uniform Fourier Data Aditya Viswanathan aditya@math.msu.edu 2014 Joint Mathematics Meetings January 18 2014 0 / 16 Joint work with Anne Gelb (Arizona State) Guohui
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 informationMEDE2500 Tutorial Nov-7
(updated 2016-Nov-4,7:40pm) MEDE2500 (2016-2017) Tutorial 3 MEDE2500 Tutorial 3 2016-Nov-7 Content 1. The Dirac Delta Function, singularity functions, even and odd functions 2. The sampling process and
More informationECE 301: Signals and Systems Homework Assignment #7
ECE 301: Signals and Systems Homework Assignment #7 Due on December 11, 2015 Professor: Aly El Gamal TA: Xianglun Mao 1 Aly El Gamal ECE 301: Signals and Systems Homework Assignment #7 Problem 1 Note:
More informationGrades will be determined by the correctness of your answers (explanations are not required).
6.00 (Fall 20) Final Examination December 9, 20 Name: Kerberos Username: Please circle your section number: Section Time 2 am pm 4 2 pm Grades will be determined by the correctness of your answers (explanations
More informationTwo-Dimensional Signal Processing and Image De-noising
Two-Dimensional Signal Processing and Image De-noising Alec Koppel, Mark Eisen, Alejandro Ribeiro March 12, 2018 Until now, we considered (one-dimensional) discrete signals of the form x : [0, N 1] C of
More informationRepresentation 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 informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Filtering in the Frequency Domain http://www.ee.unlv.edu/~b1morris/ecg782/ 2 Outline Background
More informationHonors Advanced Mathematics Determinants page 1
Determinants page 1 Determinants For every square matrix A, there is a number called the determinant of the matrix, denoted as det(a) or A. Sometimes the bars are written just around the numbers of the
More informationConstructing Approximation Kernels for Non-Harmonic Fourier Data
Constructing Approximation Kernels for Non-Harmonic Fourier Data Aditya Viswanathan aditya.v@caltech.edu California Institute of Technology SIAM Annual Meeting 2013 July 10 2013 0 / 19 Joint work with
More informationPART 1. Review of DSP. f (t)e iωt dt. F(ω) = f (t) = 1 2π. F(ω)e iωt dω. f (t) F (ω) The Fourier Transform. Fourier Transform.
PART 1 Review of DSP Mauricio Sacchi University of Alberta, Edmonton, AB, Canada The Fourier Transform F() = f (t) = 1 2π f (t)e it dt F()e it d Fourier Transform Inverse Transform f (t) F () Part 1 Review
More informationLAB 6: FIR Filter Design Summer 2011
University of Illinois at Urbana-Champaign Department of Electrical and Computer Engineering ECE 311: Digital Signal Processing Lab Chandra Radhakrishnan Peter Kairouz LAB 6: FIR Filter Design Summer 011
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 informationRepresentation 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 informationRAD229: Final Exam 2014/ SOLUTIONS You will have 3 hours to complete this Exam
RAD229: Final Exam 2014/2015 - SOLUTIONS You will have 3 hours to complete this Exam Solutions are given in Blue. In some cases, different interpretations may have led to different, but reasonable answers,
More informationECE 350 Signals and Systems Spring 2011 Final Exam - Solutions. Three 8 ½ x 11 sheets of notes, and a calculator are allowed during the exam.
ECE 35 Spring - Final Exam 9 May ECE 35 Signals and Systems Spring Final Exam - Solutions Three 8 ½ x sheets of notes, and a calculator are allowed during the exam Write all answers neatly and show your
More informationECE 301 Fall 2010 Division 2 Homework 10 Solutions. { 1, if 2n t < 2n + 1, for any integer n, x(t) = 0, if 2n 1 t < 2n, for any integer n.
ECE 3 Fall Division Homework Solutions Problem. Reconstruction of a continuous-time signal from its samples. Consider the following periodic signal, depicted below: {, if n t < n +, for any integer n,
More information2A1H Time-Frequency Analysis II Bugs/queries to HT 2011 For hints and answers visit dwm/courses/2tf
Time-Frequency Analysis II (HT 20) 2AH 2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 20 For hints and answers visit www.robots.ox.ac.uk/ dwm/courses/2tf David Murray. A periodic
More informationPost-Midterm Course Review
Post-Midterm Course Review EE 396B, Bloch & EPG, Gradient Echo Methods After Midterm: Spin-Echo Methods Sampling Radial, Spiral, EPI Measurement and Mapping Motion Diffusion 37 Spin Echo Sequences 2D Interleaved:
More informationEE 261 The Fourier Transform and its Applications Fall 2006 Final Exam Solutions. Notes: There are 7 questions for a total of 120 points
EE 6 The Fourier Transform and its Applications Fall 6 Final Exam Solutions Notes: There are 7 questions for a total of points Write all your answers in your exam booklets When there are several parts
More informationFall 2011, EE123 Digital Signal Processing
Lecture 6 Miki Lustig, UCB September 11, 2012 Miki Lustig, UCB DFT and Sampling the DTFT X (e jω ) = e j4ω sin2 (5ω/2) sin 2 (ω/2) 5 x[n] 25 X(e jω ) 4 20 3 15 2 1 0 10 5 1 0 5 10 15 n 0 0 2 4 6 ω 5 reconstructed
More informationDigital Image Processing. Filtering in the Frequency Domain
2D Linear Systems 2D Fourier Transform and its Properties The Basics of Filtering in Frequency Domain Image Smoothing Image Sharpening Selective Filtering Implementation Tips 1 General Definition: System
More informationEE422G Homework #9 Solution
EE4G Homework #9 Solution. (3 points) Sampling and Reconstruction a) Given the following discrete-time input x(n), sketch the outputs reconstructed continuous signal y(t) using Sample-And-Hold and Linear
More informationCh.11 The Discrete-Time Fourier Transform (DTFT)
EE2S11 Signals and Systems, part 2 Ch.11 The Discrete-Time Fourier Transform (DTFT Contents definition of the DTFT relation to the -transform, region of convergence, stability frequency plots convolution
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 informationCITS 4402 Computer Vision
CITS 4402 Computer Vision Prof Ajmal Mian Adj/A/Prof Mehdi Ravanbakhsh, CEO at Mapizy (www.mapizy.com) and InFarm (www.infarm.io) Lecture 04 Greyscale Image Analysis Lecture 03 Summary Images as 2-D signals
More informationEE123 Digital Signal Processing. M. Lustig, EECS UC Berkeley
EE123 Digital Signal Processing Today Last time: DTFT - Ch 2 Today: Continue DTFT Z-Transform Ch. 3 Properties of the DTFT cont. Time-Freq Shifting/modulation: M. Lustig, EE123 UCB M. Lustig, EE123 UCB
More informationMidterm 1 Review. Distance = (x 1 x 0 ) 2 + (y 1 y 0 ) 2.
Midterm 1 Review Comments about the midterm The midterm will consist of five questions and will test on material from the first seven lectures the material given below. No calculus either single variable
More informationEE591 Project Report December 2 nd, 2005
EE591 Project Report December 2 nd, 2005 Amrita Rajagopalan Department of Biomedical Engineering. arajagop@usc.edu Two Dimensional Spatially Selective RF Pulse Design Using Spiral Trajectory Abstract In
More information2A1H Time-Frequency Analysis II
2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 209 For any corrections see the course page DW Murray at www.robots.ox.ac.uk/ dwm/courses/2tf. (a) A signal g(t) with period
More information18.085: Summer 2016 Due: 3 August 2016 (in class) Problem Set 8
Problem Set 8 Unless otherwise specified, you may use MATLAB to assist with computations. provide a print-out of the code used and its output with your assignment. Please 1. More on relation between Fourier
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 informationFROM ANALOGUE TO DIGITAL
SIGNALS AND SYSTEMS: PAPER 3C1 HANDOUT 7. Dr David Corrigan 1. Electronic and Electrical Engineering Dept. corrigad@tcd.ie www.mee.tcd.ie/ corrigad FROM ANALOGUE TO DIGITAL To digitize signals it is necessary
More informationVID3: Sampling and Quantization
Video Transmission VID3: Sampling and Quantization By Prof. Gregory D. Durgin copyright 2009 all rights reserved Claude E. Shannon (1916-2001) Mathematician and Electrical Engineer Worked for Bell Labs
More informationLAB # 5 HANDOUT. »»» The N-point DFT is converted into two DFTs each of N/2 points. N = -W N Then, the following formulas must be used. = k=0,...
EEE4 Lab Handout. FAST FOURIER TRANSFORM LAB # 5 HANDOUT Data Sequence A = x, x, x, x3, x4, x5, x6, x7»»» The N-point DFT is converted into two DFTs each of N/ points. x, x, x4, x6 x, x3, x5, x7»»» N =e
More informationLecture 7: Statistics and the Central Limit Theorem. Philip Moriarty,
Lecture 7: Statistics and the Central Limit Theorem Philip Moriarty, philip.moriarty@nottingham.ac.uk NB Notes based heavily on lecture slides prepared by DE Rourke for the F32SMS module, 2006 7.1 Recap
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 informationCourse Review. Midterm Review: EE369B Concepts Simulations with Bloch Matrices, EPG SNR. B.Hargreaves - RAD 229. Section F1
Course Review Midterm Review: EE369B Concepts Simulations with Bloch Matrices, EPG SNR 1 Section F1 Bloch/Matrix Simulations M = [Mx My Mz] T RF and precession ~ 3x3 rotation matrices Relaxation ~ 3x3
More information(Refer Slide Time: 01:28 03:51 min)
Digital Signal Processing Prof. S. C. Dutta Roy Department of Electrical Engineering Indian Institute of Technology, Delhi Lecture 40 FIR Design by Windowing This is the 40 th lecture and our topic for
More informationLecture 04 Image Filtering
Institute of Informatics Institute of Neuroinformatics Lecture 04 Image Filtering Davide Scaramuzza 1 Lab Exercise 2 - Today afternoon Room ETH HG E 1.1 from 13:15 to 15:00 Work description: your first
More informationMachine vision, spring 2018 Summary 4
Machine vision Summary # 4 The mask for Laplacian is given L = 4 (6) Another Laplacian mask that gives more importance to the center element is given by L = 8 (7) Note that the sum of the elements in the
More informationELEN 4810 Midterm Exam
ELEN 4810 Midterm Exam Wednesday, October 26, 2016, 10:10-11:25 AM. One sheet of handwritten notes is allowed. No electronics of any kind are allowed. Please record your answers in the exam booklet. Raise
More informationComputational 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 information4. Sinusoidal solutions
16 4. Sinusoidal solutions Many things in nature are periodic, even sinusoidal. We will begin by reviewing terms surrounding periodic functions. If an LTI system is fed a periodic input signal, we have
More information2 Fourier Transforms and Sampling
2 Fourier ransforms and Sampling 2.1 he Fourier ransform he Fourier ransform is an integral operator that transforms a continuous function into a continuous function H(ω) =F t ω [h(t)] := h(t)e iωt dt
More informationL6: Short-time Fourier analysis and synthesis
L6: Short-time Fourier analysis and synthesis Overview Analysis: Fourier-transform view Analysis: filtering view Synthesis: filter bank summation (FBS) method Synthesis: overlap-add (OLA) method STFT magnitude
More informationE : Lecture 1 Introduction
E85.2607: Lecture 1 Introduction 1 Administrivia 2 DSP review 3 Fun with Matlab E85.2607: Lecture 1 Introduction 2010-01-21 1 / 24 Course overview Advanced Digital Signal Theory Design, analysis, and implementation
More informationMITOCW watch?v=rf5sefhttwo
MITOCW watch?v=rf5sefhttwo The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high-quality educational resources for free. To
More informationSound & Vibration Magazine March, Fundamentals of the Discrete Fourier Transform
Fundamentals of the Discrete Fourier Transform Mark H. Richardson Hewlett Packard Corporation Santa Clara, California The Fourier transform is a mathematical procedure that was discovered by a French mathematician
More informationChapter 4: Filtering in the Frequency Domain. Fourier Analysis R. C. Gonzalez & R. E. Woods
Fourier Analysis 1992 2008 R. C. Gonzalez & R. E. Woods Properties of δ (t) and (x) δ : f t) δ ( t t ) dt = f ( ) f x) δ ( x x ) = f ( ) ( 0 t0 x= ( 0 x0 1992 2008 R. C. Gonzalez & R. E. Woods Sampling
More informationFourier Transform 2D
Image Processing - Lesson 8 Fourier Transform 2D Discrete Fourier Transform - 2D Continues Fourier Transform - 2D Fourier Properties Convolution Theorem Eamples = + + + The 2D Discrete Fourier Transform
More informationToday s lecture. The Fourier transform. Sampling, aliasing, interpolation The Fast Fourier Transform (FFT) algorithm
Today s lecture The Fourier transform What is it? What is it useful for? What are its properties? Sampling, aliasing, interpolation The Fast Fourier Transform (FFT) algorithm Jean Baptiste Joseph Fourier
More informationDigital Signal Processing. Midterm 1 Solution
EE 123 University of California, Berkeley Anant Sahai February 15, 27 Digital Signal Processing Instructions Midterm 1 Solution Total time allowed for the exam is 8 minutes Some useful formulas: Discrete
More informationSimilar Shapes and Gnomons
Similar Shapes and Gnomons May 12, 2013 1. Similar Shapes For now, we will say two shapes are similar if one shape is a magnified version of another. 1. In the picture below, the square on the left is
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 informationContinuous Fourier transform of a Gaussian Function
Continuous Fourier transform of a Gaussian Function Gaussian function: e t2 /(2σ 2 ) The CFT of a Gaussian function is also a Gaussian function (i.e., time domain is Gaussian, then the frequency domain
More informationMath 56 Homework 5 Michael Downs
1. (a) Since f(x) = cos(6x) = ei6x 2 + e i6x 2, due to the orthogonality of each e inx, n Z, the only nonzero (complex) fourier coefficients are ˆf 6 and ˆf 6 and they re both 1 2 (which is also seen from
More informationVectors and Scalars. Scalar: A quantity specified by its magnitude only Vector: A quantity specified both by its magnitude and direction.
Vectors and Scalars Scalar: A quantity specified by its magnitude only Vector: A quantity specified both by its magnitude and direction. To distinguish a vector from a scalar quantity, it is usually written
More informationHow to manipulate Frequencies in Discrete-time Domain? Two Main Approaches
How to manipulate Frequencies in Discrete-time Domain? Two Main Approaches Difference Equations (an LTI system) x[n]: input, y[n]: output That is, building a system that maes use of the current and previous
More informationDOING PHYSICS WITH MATLAB ELECTRIC FIELD AND ELECTRIC POTENTIAL: POISSON S EQUATION
DOING PHYSICS WITH MATLAB ELECTRIC FIELD AND ELECTRIC POTENTIAL: POISSON S EQUATION Ian Cooper School of Physics, University of Sydney ian.cooper@sydney.edu.au DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS For
More informationLab #2 - Two Degrees-of-Freedom Oscillator
Lab #2 - Two Degrees-of-Freedom Oscillator Last Updated: March 0, 2007 INTRODUCTION The system illustrated in Figure has two degrees-of-freedom. This means that two is the minimum number of coordinates
More information6.869 Advances in Computer Vision. Bill Freeman, Antonio Torralba and Phillip Isola MIT Oct. 3, 2018
6.869 Advances in Computer Vision Bill Freeman, Antonio Torralba and Phillip Isola MIT Oct. 3, 2018 1 Sampling Sampling Pixels Continuous world 3 Sampling 4 Sampling 5 Continuous image f (x, y) Sampling
More informationCompensation. June 29, Arizona State University. Edge Detection from Nonuniform Data via. Convolutional Gridding with Density.
Arizona State University June 29, 2012 Overview gridding with density compensation is a common method for function reconstruction nonuniform Fourier data. It is computationally inexpensive compared with
More information1 1.27z z 2. 1 z H 2
E481 Digital Signal Processing Exam Date: Thursday -1-1 16:15 18:45 Final Exam - Solutions Dan Ellis 1. (a) In this direct-form II second-order-section filter, the first stage has
More informationSignal Processing COS 323
Signal Processing COS 323 Digital Signals D: functions of space or time e.g., sound 2D: often functions of 2 spatial dimensions e.g. images 3D: functions of 3 spatial dimensions CAT, MRI scans or 2 space,
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 informationExperimental Fourier Transforms
Chapter 5 Experimental Fourier Transforms 5.1 Sampling and Aliasing Given x(t), we observe only sampled data x s (t) = x(t)s(t; T s ) (Fig. 5.1), where s is called sampling or comb function and can be
More informationEE123 Digital Signal Processing
EE123 Digital Signal Processing Lecture 1 Time-Dependent FT Announcements! Midterm: 2/22/216 Open everything... but cheat sheet recommended instead 1am-12pm How s the lab going? Frequency Analysis with
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 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 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 informationDepartment of Electrical and Telecommunications Engineering Technology TEL (718) FAX: (718) Courses Description:
NEW YORK CITY COLLEGE OF TECHNOLOGY The City University of New York 300 Jay Street Brooklyn, NY 11201-2983 Department of Electrical and Telecommunications Engineering Technology TEL (718) 260-5300 - FAX:
More informationAcoustic holography. LMS Test.Lab. Rev 12A
Acoustic holography LMS Test.Lab Rev 12A Copyright LMS International 2012 Table of Contents Chapter 1 Introduction... 5 Chapter 2... 7 Section 2.1 Temporal and spatial frequency... 7 Section 2.2 Time
More informationLesson 8: Work and Energy
Name Period Lesson 8: Work and Energy 8.1 Experiment: What is Kinetic Energy? (a) Set up the cart, meter stick, pulley, hanging mass, and tape as you did in Lesson 5.1. You will examine the distance and
More informationHST.582J / 6.555J / J Biomedical Signal and Image Processing Spring 2007
MIT OpenCourseare http://ocw.mit.edu HST.58J / 6.555J / 16.56J Biomedical Signal and Image Processing Spring 7 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationEE123 Digital Signal Processing
EE123 Digital Signal Processing Discrete Time Fourier Transform M. Lustig, EECS UC Berkeley A couple of things Read Ch 2 2.0-2.9 It s OK to use 2nd edition Class webcast in bcourses.berkeley.edu or linked
More informationRF Pulse Design. Multi-dimensional Excitation I. M229 Advanced Topics in MRI Kyung Sung, Ph.D Class Business
RF Pulse Design Multi-dimensional Excitation I M229 Advanced Topics in MRI Kyung Sung, Ph.D. 2018.04.10 Class Business Office hours - Instructors: Fri 10-12pm TAs: Xinran Zhong and Zhaohuan Zhang (time:
More information2 Background: Fourier Series Analysis and Synthesis
Signal Processing First Lab 15: Fourier Series Pre-Lab and Warm-Up: You should read at least the Pre-Lab and Warm-up sections of this lab assignment and go over all exercises in the Pre-Lab section before
More informationFinal. Introduction to Artificial Intelligence. CS 188 Spring You have approximately 2 hours and 50 minutes.
CS 188 Spring 2014 Introduction to Artificial Intelligence Final You have approximately 2 hours and 50 minutes. The exam is closed book, closed notes except your two-page crib sheet. Mark your answers
More informationFiltering and Edge Detection
Filtering and Edge Detection Local Neighborhoods Hard to tell anything from a single pixel Example: you see a reddish pixel. Is this the object s color? Illumination? Noise? The next step in order of complexity
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 informationComputer Vision Lecture 3
Computer Vision Lecture 3 Linear Filters 03.11.2015 Bastian Leibe RWTH Aachen http://www.vision.rwth-aachen.de leibe@vision.rwth-aachen.de Demo Haribo Classification Code available on the class website...
More information