EE123 Digital Signal Processing
|
|
- Pauline Anderson
- 6 years ago
- Views:
Transcription
1 Multi-Dimensional Signals EE23 Digital Signal Processing Lecture Our world is more complex than D Images: f(x,y) Videos: f(x,y,t) Dynamic 3F scenes: f(x,y,z,t) Medical Imaging 3D Video Computer Graphics We will focus on 2D Continuous-Time 2D functions δ(x,y): Impulse at x=0, y=0 δ(x) : Impulse line (vertical or horizontal?) (x,y) : 2D rect function cos(2π(fxx + fyy) - Spatial harmonic Circularly Symmetic: δ(r): Impulse ring (x,y): Pillbox u(r) (r r 0 ) 2 r o Spatial Frequency What is a spatial frequency? Complex Harmonic: e j( xx+ y y) = e j2 (f xx+f y y) Units (for example): x,y - cm fx,fy - /cm Ωx, Ωy - rad/cm
2 Spatial Frequency Vinyl Record Transforms a temporal signal to a spatial signal Spatial Frequency CD ROM encodes digital temporal signals to spatial signals What is the frequency? What is the frequency? (-,) sin(2 (f x x + f y y)) (-,) sin(2 (f x x + f y y)) (-,-) a) fx=2, fy=2 b) fx=, fy=0 (,-) c) fx = 4, fy=0 d) none of the above (-,-) (,-) 2 cycles for 2 cm fx= cm -
3 What is the frequency? What is the frequency? sin(2 (f x x + f y y)) (-,) or cos(2 (f x x + f y y)) (-,) cos(2 (f x x + f y y)) (-,-) (,-) (-,-) (,-) a) sin, fx=0, fy=2 c) cos, fx = 0, fy=2 a) fx=, fy=2 c) fx = 2, fy= b) cos, fx=0, fy=4 d) none of the above b) fx=4, fy=2 d) fx=2, fy=4 What is the frequency? What is the Temporal Frequency? cos(2 (f x x + f y y)) Vinyl rotates at Hz (-,) (-,) (-,-) (,-) (-,-) (,-) a) fx=, fy=2 c) fx = 2, fy= a) cos(2π8t) c) cos(2π4t) b) fx=4, fy=2 d) fx=2, fy=4 b) cos(2π8t 2 ) d) cos(2π4t 2 )
4 What is the Temporal Frequency? (-,) Vinyl rotates at Hz Challenge: What is the password? Aliasing! 2sec (-,-) a) cos(2π00t) b) cos(2π00t 2 ) (,-) c) cos(2π40t) d) none of the answers 2D Fourier Transform Z Z F (f x,f y )= f(x, y)e j2 (f xx+f y y) dxdy
5 sinc(ax)sinc(by) 2D DTFT F (! x,! y )= X X f[n x,n y ]e j(! xn x +! y n y ) n x = n y = apple! x,! y apple sinc(ax)sinc(by) F (apple x, apple y )= I prefer 2nd X X n x = n y = f[n x,n y ]e j2 (apple xn x +apple y n y ) 0.5 apple apple x, apple y apple 0.5 Massaging the DTFT leads to separable transforms in each axis
6 Separability of 2D DTFT white board... 2D - DFT Similarly to D: Forward: F [k x,k y ]= Inverse: NX MX n x =0 n y =0 f[n x,n y ]e j2 (n xk x /N +n y k y /M ) apple x = k x /N, apple y = k y /M f[n x,n y ]= NM NX MX k x =0 k y =0 F [k x,k y ]e +j2 (n xk x /N +n y k y /M ) 2D - DFT F [k x,k y ]= M NX MX n x =0 n y =0 f[n x,n y ]e j2 (n xk x /N +n y k y /M ) f[n x,n y ] F [k x,k y ] (0,) Properties of 2D DFT Circular Convolution f[n x,n y ] h[n x,n y ]=F [k x,k y ]H[k x,k y ] ny ky Circular shift (0,0) N nx (0,0) kx (,0) Need to fftshift in 2D to get it to look like DTFT. f[(n x m x ) N, (n y m y ) M ]=e j2 (k xm x /N +k y m y /M ) F [k x,k y ]
7
8 . Projection w Many Projections - Tomography (A fr:x,j) J(f-xco,e-j>l'ne)rfx:dj / \:. J projection x-ray
9 / \:. J. Radon Transform w Radon Transform: Sinogram (A fr:x,j) J(f-xco,e-j>l'ne)rfx:dj / \:. J Computed Tomography Projection Slice Theorem (Bracewell) Sinogram cross-section,f F0<. (), ey'" fy{) FBP M. lrd ft spatial frequency x-ray source central-section/ projection slice theorem (Bracewell 56)?voo P: o-tu : Ar fcx,j)dj F( /.;_,o) =- JJ G "' o ='l k, or6f 9 /J-fr{ t,:x-.. f; jj d-xjj =-
10 Projection Slice Theorem (Bracewell) Discrete Reconstruction?voo P : o-tu : A r G "' o ='l k, or6f 9.. fcx,j)dj F( /.;_,o) =- J J /J-fr{ t,:x-.. f; jj d-xjj =- -\A - (.) ' Discretization /)L{c(2EJ7.7:.[.'
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 informationCentral Slice Theorem
Central Slice Theorem Incident X-rays y f(x,y) R x r x Detected p(, x ) The thick line is described by xcos +ysin =R Properties of Fourier Transform F [ f ( x a)] F [ f ( x)] e j 2 a Spatial Domain Spatial
More informationModern CT system generations Measurement of attenuation
CT reconstruction repetition & hints Reconstruction in CT and hints to the assignments Jørgen Arendt Jensen October 4, 16 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering
More informationReconstruction in CT and hints to the assignments
Reconstruction in CT and hints to the assignments Jørgen Arendt Jensen October 24, 2016 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound Imaging
More informationDigital Image Processing
Digital Image Processing Image Restoration and Reconstruction (Image Reconstruction from Projections) Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science and Engineering
More informationBiomedical Image Analysis. Spatial Filtering
Biomedical Image Analysis Contents: Spatial Filtering The mechanics of Spatial Filtering Smoothing and sharpening filters BMIA 15 V. Roth & P. Cattin 1 The Mechanics of Spatial Filtering Spatial filter:
More information9. Computed Tomography (CT) 9.1 Absorption of X-ray
Absorption of X-ray 9. 9.1 Absorption of X-ray X-ray radiography Object I κ Incident Transmitted d I(y) = I e κ(y)d κ : attenuation coefficient In the case of X-ray, it depends on the atomic number. (Heavy
More informationRadon Transform and Filtered Backprojection
Radon Transform and Filtered Backprojection Jørgen Arendt Jensen October 13, 2016 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound Imaging Department
More informationContinuous Space Fourier Transform (CSFT)
C. A. Bouman: Digital Image Processing - January 8, 8 Continuous Space Fourier Transform (CSFT) Forward CSFT: F(u,v) = Inverse CSFT: f(x,y) = f(x,y)e jπ(ux+vy) dxdy F(u,v)e jπ(ux+vy) dudv Space coordinates:.
More informationLecture 5: Frequency Domain Transformations
#1 Lecture 5: Frequency Domain Transformations Saad J Bedros sbedros@umn.edu From Last Lecture Spatial Domain Transformation Point Processing for Enhancement Area/Mask Processing Transformations Image
More informationLecture 2: 2D Fourier transforms and applications
Lecture 2: 2D Fourier transforms and applications B14 Image Analysis Michaelmas 2017 Dr. M. Fallon Fourier transforms and spatial frequencies in 2D Definition and meaning The Convolution Theorem Applications
More informationEECS490: Digital Image Processing. Lecture #16
Lecture #16 Wiener Filters Constrained Least Squares Filter Computed Tomography Basics Reconstruction and the Radon Transform Fourier Slice Theorem Filtered Backprojections Fan Beams Motion Blurring Model
More informationFourier analysis and sampling theory
Reading Required: Shirley, Ch. 9 Recommended: Fourier analysis and sampling theory Ron Bracewell, The Fourier Transform and Its Applications, McGraw-Hill. Don P. Mitchell and Arun N. Netravali, Reconstruction
More informationComputer Vision and Graphics (ee2031) Digital Image Processing I
Computer Vision and Graphics (ee203) Digital Image Processing I Dr John Collomosse J.Collomosse@surrey.ac.uk Centre for Vision, Speech and Signal Processing University of Surrey Learning Outcomes After
More informationSampling: Application to 2D Transformations
Sampling: Application to 2D Transformations University of the Philippines - Diliman August Diane Lingrand lingrand@polytech.unice.fr http://www.essi.fr/~lingrand Sampling Computer images are manipulated
More informationDigital Image Processing. Lecture 6
Digital Image Processing Lecture 6 (Enhancement in the Frequency domain) Bu-Ali Sina University Computer Engineering Dep. Fall 2016 Image Enhancement In The Frequency Domain Outline Jean Baptiste Joseph
More informationBME I5000: Biomedical Imaging
1 Lucas Parra, CCNY BME I5000: Biomedical Imaging Lecture 11 Point Spread Function, Inverse Filtering, Wiener Filtering, Sharpening,... Lucas C. Parra, parra@ccny.cuny.edu Blackboard: http://cityonline.ccny.cuny.edu/
More informationComputer Vision 2. SS 18 Dr. Benjamin Guthier Professur für Bildverarbeitung. Computer Vision 2 Dr. Benjamin Guthier
Computer Vision 2 SS 18 Dr. Benjamin Guthier Professur für Bildverarbeitung Computer Vision 2 Dr. Benjamin Guthier 1. IMAGE PROCESSING Computer Vision 2 Dr. Benjamin Guthier Content of this Chapter Non-linear
More informationDigital Signal Processing Lecture Notes 22 November 2010
Digital Signal Processing Lecture otes 22 ovember 2 Topics: Discrete Cosine Transform FFT Linear and Circular Convolution Rate Conversion Includes review of Fourier transforms, properties of Fourier transforms,
More informationTomographic Reconstruction
Tomographic Reconstruction 3D Image Processing Torsten Möller Reading Gonzales + Woods, Chapter 5.11 2 Overview Physics History Reconstruction basic idea Radon transform Fourier-Slice theorem (Parallel-beam)
More informationTopics. Projections. Review Filtered Backprojection Fan Beam Spiral CT Applications. Bioengineering 280A Principles of Biomedical Imaging
Bioengineering 28A Principles of Biomedical Imaging Fall Quarter 24 X-Rays/CT Lecture 2 Topics Review Filtered Backprojection Fan Beam Spiral CT Applications s I θ (r) = I exp µ(x, y)ds Lr,θ = I exp µ(rcosθ
More informationTopics. Review Filtered Backprojection Fan Beam Spiral CT Applications. Bioengineering 280A Principles of Biomedical Imaging
Bioengineering 28A Principles of Biomedical Imaging Fall Quarter 24 X-Rays/CT Lecture 2 Topics Review Filtered Backprojection Fan Beam Spiral CT Applications Projections I θ (r) = I exp Lr,θ µ(x,y)ds =
More informationChapter6 Image Reconstruction
Chapter6 Image Reconstruction Preview 61I 6.1 Introduction 6.2 Reconstruction by Fourier Inversion 6.3 Reconstruction by convolution and backprojection 6.4 Finite series-expansion 1 Preview Reconstruction
More informationGE s Revolution CT MATLAB III: CT. Kathleen Chen March 20, 2018
GE s Revolution CT MATLAB III: CT Kathleen Chen chens18@rpi.edu March 20, 2018 https://www.zmescience.com/medicine/inside-human-body-real-time-gifs-demo-power-ct-scan/ Reminders Make sure you have MATLAB
More informationThe main goal of Computer Graphics is to generate 2D images 2D images are continuous 2D functions (or signals)
Motivation The main goal of Computer Graphics is to generate 2D images 2D images are continuous 2D functions (or signals) monochrome f(x,y) or color r(x,y), g(x,y), b(x,y) These functions are represented
More informationJoint ICTP-TWAS Workshop on Portable X-ray Analytical Instruments for Cultural Heritage. 29 April - 3 May, 2013
2455-5 Joint ICTP-TWAS Workshop on Portable X-ray Analytical Instruments for Cultural Heritage 29 April - 3 May, 2013 Lecture NoteBasic principles of X-ray Computed Tomography Diego Dreossi Elettra, Trieste
More informationABSTRACT INTRODUCTION
3-D f-k filtering 3-D f-k filtering Wai-Kin Chan and Robert R. Stewart ABSTRACT Stewart and Schieck (1993) show that a conventional 2-pass 2-D f-k filter fails to provide an axially-symmetric response.
More informationCOORDINATE GEOMETRY CHANGE OF AXES EXERCISE 1. The point to which the origin should be shifted in order to eliminate x and y terms in the equation 4x 9y 8x 36y 4 0 is (1, ) ( 4, ( 1,) 4) (1, ). In order
More informationMEDICAL IMAGING 2nd Part Computed Tomography
MEDICAL IMAGING 2nd Part Computed Tomography Introduction 2 In the last 30 years X-ray Computed Tomography development produced a great change in the role of diagnostic imaging in medicine. In convetional
More informationMotivation. The main goal of Computer Graphics is to generate 2D images. 2D images are continuous 2D functions (or signals)
Motivation The main goal of Computer Graphics is to generate 2D images 2D images are continuous 2D functions (or signals) monochrome f(x,y) or color r(x,y), g(x,y), b(x,y) These functions are represented
More informationImage preprocessing in spatial domain
Image preprocessing in spatial domain Sampling theorem, aliasing, interpolation, geometrical transformations Revision:.4, dated: May 25, 26 Tomáš Svoboda Czech Technical University, Faculty of Electrical
More informationComputed Tomography (Part 2)
EL-GY 6813 / BE-GY 603 / G16.446 Medical Imaging Computed Tomography (Part ) Yao Wang Polytechnic School of Engineering New York University, Brooklyn, NY 1101 Based on Prince and Links, Medical Imaging
More informationImage preprocessing in spatial domain
Image preprocessing in spatial domain Sampling theorem, aliasing, interpolation, geometrical transformations Revision:.3, dated: December 7, 25 Tomáš Svoboda Czech Technical University, Faculty of Electrical
More informationLOWELL WEEKLY JOURNAL.
N $ N N N 5 * N F N F *» ) N F N * )» N F 5 N F ** *» FXN* X N F N N * N q» N q $ $ $ q* 5 5 * q 5 * 5 q N * * F N / N N» N N X )» * * F F N * * * N * / * N F F F N F X F * N F N ** N N F * NF N N F X
More informationDigital Image Processing
Digital Image Processing Third Edition Rafael C. Gonzalez University of Tennessee Richard E. Woods MedData Interactive PEARSON Prentice Hall Pearson Education International Contents Preface xv Acknowledgments
More informationIntroduction to Medical Imaging. Cone-Beam CT. Klaus Mueller. Computer Science Department Stony Brook University
Introduction to Medical Imaging Cone-Beam CT Klaus Mueller Computer Science Department Stony Brook University Introduction Available cone-beam reconstruction methods: exact approximate algebraic Our discussion:
More informationLab 1: Basic operations on images in Matlab
Lab 1: Basic operations on images in Matlab Maria Magnusson with contributions by Michael Felsberg, 2017, Computer Vision Laboratory, Department of Electrical Engineering, Linköping University 1 Introduction
More informationTopics. 2D Signal. a b. c d. 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. 2D Signal Representations 2. 2D Fourier Transform 3. Transform Pairs 4. FT Properties 2D Signal
More informationBiophysical Techniques (BPHS 4090/PHYS 5800)
Biophysical Techniques (BPHS 4090/PHYS 5800) Instructors: Prof. Christopher Bergevin (cberge@yorku.ca) Schedule: MWF 1:30-2:30 (CB 122) Website: http://www.yorku.ca/cberge/4090w2017.html York University
More information3. Image formation, Fourier analysis and CTF theory. Paula da Fonseca
3. Image formation, Fourier analysis and CTF theory Paula da Fonseca EM course 2017 - Agenda - Overview of: Introduction to Fourier analysis o o o o Sine waves Fourier transform (simple examples of 1D
More informationImage restoration. Lecture 14. Milan Gavrilovic Centre for Image Analysis Uppsala University
Image restoration Lecture 14 Milan Gavrilovic milan@cb.uu.se Centre for Image Analysis Uppsala University Computer Assisted Image Analysis 2009-05-08 M. Gavrilovic (Uppsala University) L14 Image restoration
More informationComputational Aspects of MRI
David Atkinson Philip Batchelor David Larkman Programme 09:30 11:00 Fourier, sampling, gridding, interpolation. Matrices and Linear Algebra 11:30 13:00 MRI Lunch (not provided) 14:00 15:30 SVD, eigenvalues.
More informationReading. 2. Fourier analysis and sampling theory. Required: Watt, Section 14.1 Recommended:
Reading Required: Watt, Section 14.1 Recommended: 2. Fourier analysis and sampling theory Ron Bracewell, The Fourier Transform and Its Applications, McGraw-Hill. Don P. Mitchell and Arun N. Netravali,
More informationReconstruction in CT and relation to other imaging modalities
Reconstruction in CT and relation to other imaging modalities Jørgen Arendt Jensen November 16, 2015 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound
More informationLECTURE 3-1 AREA OF A REGION BOUNDED BY CURVES
7 CALCULUS II DR. YOU 98 LECTURE 3- AREA OF A REGION BOUNDED BY CURVES If y = f(x) and y = g(x) are continuous on an interval [a, b] and f(x) g(x) for all x in [a, b], then the area of the region between
More informationELEC Dr Reji Mathew Electrical Engineering UNSW
ELEC 4622 Dr Reji Mathew Electrical Engineering UNSW Dynamic Range and Weber s Law HVS is capable of operating over an enormous dynamic range, However, sensitivity is far from uniform over this range Example:
More informationOptics Quiz #2 April 30, 2014
.71 - Optics Quiz # April 3, 14 Problem 1. Billet s Split Lens Setup I the ield L(x) is placed against a lens with ocal length and pupil unction P(x), the ield (X) on the X-axis placed a distance behind
More informationImage Restoration. Diffusion Denoising Deconvolution Super-resolution Tomographic Reconstruction
Image Restoration Image Restoration Diffusion Denoising Deconvolution Super-resolution Tomographic Reconstruction Diffusion Term Consider only the regularization term E-L equation: (Laplace equation) Steepest
More informationIndex. aliasing artifacts and noise in CT images, 200 measurement of projection data, nondiffracting
Index Algebraic equations solution by Kaczmarz method, 278 Algebraic reconstruction techniques, 283-84 sequential, 289, 293 simultaneous, 285-92 Algebraic techniques reconstruction algorithms, 275-96 Algorithms
More informationMedical Image Reconstruction Term II 2012 Topic 6: Tomography
Medical Image Reconstruction Term II 2012 Topic 6: Tomography Professor Yasser Mostafa Kadah Tomography The Greek word tomos means a section, a slice, or a cut. Tomography is the process of imaging a cross
More informationX-ray tomography. X-ray tomography. Applications in Science. X-Rays. Notes. Notes. Notes. Notes
X-ray tomography Important application of the Fast Fourier transform: X-ray tomography. Also referred to as CAT scan (Computerized Axial Tomography) X-ray tomography This has revolutionized medical diagnosis.
More informationTutorial 5. Jun Xu, Teaching Asistant March 2, COMP4134 Biometrics Authentication
Tutorial 5 Jun Xu, Teaching Asistant nankaimathxujun@gmail.com COMP4134 Biometrics Authentication March 2, 2017 Table of Contents Problems Problem 1: Answer The Questions Problem 2: Indeterminate Region
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 2010 CT/Fourier Lecture 5 CT Sampling Requirements View Aliasing What should
More informationContents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Homework 1 - Solutions 3. 2 Homework 2 - Solutions 13
MATH 32B-2 (8) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Homework - Solutions 3 2 Homework 2 - Solutions 3 3 Homework 3 - Solutions 9 MATH 32B-2 (8) (L) G. Liu / (TA) A. Zhou Calculus
More informationCoE4TN4 Image Processing. Chapter 5 Image Restoration and Reconstruction
CoE4TN4 Image Processing Chapter 5 Image Restoration and Reconstruction Image Restoration Similar to image enhancement, the ultimate goal of restoration techniques is to improve an image Restoration: a
More informationAn Introduc+on to Mathema+cal Image Processing IAS, Park City Mathema2cs Ins2tute, Utah Undergraduate Summer School 2010
An Introduc+on to Mathema+cal Image Processing IAS, Park City Mathema2cs Ins2tute, Utah Undergraduate Summer School 2010 Luminita Vese Todd WiCman Department of Mathema2cs, UCLA lvese@math.ucla.edu wicman@math.ucla.edu
More informationA Novel Two-step Method for CT Reconstruction
A Novel Two-step Method for CT Reconstruction Michael Felsberg Computer Vision Laboratory, Dept. EE, Linköping University, Sweden mfe@isy.liu.se Abstract. In this paper we address the parallel beam 2D
More informationRefractive Index Map Reconstruction in Optical Deflectometry Using Total-Variation Regularization
Refractive Index Map Reconstruction in Optical Deflectometry Using Total-Variation Regularization L. Jacques 1, A. González 2, E. Foumouo 2 and P. Antoine 2 1 Institute of Information and Communication
More informationInvestigation on reconstruction methods applied to 3D terahertz computed Tomography
Investigation on reconstruction methods applied to 3D terahertz computed Tomography B. Recur, 3 A. Younus, 1, P. Mounaix 1, S. Salort, 2 B. Chassagne, 2 P. Desbarats, 3 J-P. Caumes, 2 and E. Abraham 1
More informationange Data Reconstruction Using Fourier Slice Theorern*
ange Data Reconstruction Using Fourier Slice Theorern* Shih-Schon Lin and Chiou-Shann Fuh Department of nformation Engineering and Computer Science National Taiwan University E-mail: r3506011@image.csie.ntu.edu.tw,
More informationComputer Vision I. Announcements. Fourier Tansform. Efficient Implementation. Edge and Corner Detection. CSE252A Lecture 13.
Announcements Edge and Corner Detection HW3 assigned CSE252A Lecture 13 Efficient Implementation Both, the Box filter and the Gaussian filter are separable: First convolve each row of input image I with
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) 1D: Common Transform Pairs Summary source FT Properties: Convolution See book DIP 4.2.5:
More informationLecture Image Enhancement and Spatial Filtering
Lecture Image Enhancement and Spatial Filtering Harvey Rhody Chester F. Carlson Center for Imaging Science Rochester Institute of Technology rhody@cis.rit.edu September 29, 2005 Abstract Applications of
More informationAliasing. Can t draw smooth lines on discrete raster device get staircased lines ( jaggies ):
(Anti)Aliasing and Image Manipulation for (y = 0; y < Size; y++) { for (x = 0; x < Size; x++) { Image[x][y] = 7 + 8 * sin((sqr(x Size) + SQR(y Size)) / 3.0); } } // Size = Size / ; Aliasing Can t draw
More informationEEE 512 ADVANCED DIGITAL SIGNAL AND IMAGE PROCESSING
UNIVERSITI SAINS MALAYSIA Semester I Examination Academic Session 27/28 October/November 27 EEE 52 ADVANCED DIGITAL SIGNAL AND IMAGE PROCESSING Time : 3 hours INSTRUCTION TO CANDIDATE: Please ensure that
More informationTomography. Forward projectionsp θ (r) are known as a Radon transform. Objective: reverse this process to form the original image
C. A. Bouman: Digital Image Processing - January 9, 217 1 Tomography Many medical imaging systems can only measure projections through an object with density f(x,y). Projections must be collected at every
More informationThis is a good time to refresh your memory on double-integration. We will be using this skill in the upcoming lectures.
Chapter 5: JOINT PROBABILITY DISTRIBUTIONS Part 1: Sections 5-1.1 to 5-1.4 For both discrete and continuous random variables we will discuss the following... Joint Distributions (for two or more r.v. s)
More informationSteen Moeller Center for Magnetic Resonance research University of Minnesota
Steen Moeller Center for Magnetic Resonance research University of Minnesota moeller@cmrr.umn.edu Lot of material is from a talk by Douglas C. Noll Department of Biomedical Engineering Functional MRI Laboratory
More informationCorso di laurea in Fisica A.A Fisica Medica 4 TC
Corso di laurea in Fisica A.A. 2007-2008 Fisica Medica 4 TC Computed Tomography Principles 1. Projection measurement 2. Scanner systems 3. Scanning modes Basic Tomographic Principle The internal structure
More informationExam 2 is Tue Nov 21. Bring a pencil and a calculator. Discuss similarity to exam1. HW3 is due Tue Dec 5.
Stat 100a: Introduction to Probability. Outline for the day 1. Bivariate and marginal density. 2. CLT. 3. CIs. 4. Sample size calculations. 5. Review for exam 2. Exam 2 is Tue Nov 21. Bring a pencil and
More informationTheoretically Perfect Sensor
Sampling 1/67 Sampling The ray tracer samples the geometry, only gathering information from the parts of the world that interact with a finite number of rays In contrast, a scanline renderer can push all
More informationReconstruction in CT and relation to other imaging modalities
Reconstruction in CT and relation to other imaging modalities Jørgen Arendt Jensen November 1, 2017 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound
More informationSampling, Ordering, Interleaving
Sampling, Ordering, Interleaving Sampling patterns and PSFs View ordering Modulation due to transients Temporal modulations Slice interleaving Sequential, Odd/even, bit-reversed Arbitrary Other considerations:
More informationTheoretically Perfect Sensor
Sampling 1/60 Sampling The ray tracer samples the geometry, only gathering information from the parts of the world that interact with a finite number of rays In contrast, a scanline renderer can push all
More informationVisual Pathways to the Brain
Visual Pathways to the Brain 1 Left half of visual field which is imaged on the right half of each retina is transmitted to right half of brain. Vice versa for right half of visual field. From each eye
More informationImage Restoration. Yao Wang Polytechnic Institute of NYU, Brooklyn, NY 11201
Image Restoration Yao Wang Polytechnic Institute of NYU, Brooklyn, NY 11201 Partly based on A. K. Jain, Fundamentals of Digital Image Processing, and Gonzalez/Woods, Digital Image Processing Figures from
More informationLaser sensors. Transmitter. Receiver. Basilio Bona ROBOTICA 03CFIOR
Mobile & Service Robotics Sensors for Robotics 3 Laser sensors Rays are transmitted and received coaxially The target is illuminated by collimated rays The receiver measures the time of flight (back and
More informationImage Enhancement in Spatial Domain. By Dr. Rajeev Srivastava
Image Enhancement in Spatial Domain By Dr. Rajeev Srivastava CONTENTS Image Enhancement in Spatial Domain Spatial Domain Methods 1. Point Processing Functions A. Gray Level Transformation functions for
More informationExamination in Digital Image Processing, TSBB08
Examination in Digital Image Processing, TSBB8 Time: 215-1-23, 14.-18. Location: TER4 Examiner: Maria Magnusson will visit the classroom at 15. and 16.45, tel. 281336, 73-84 38 67 Permitted equipment:
More informationMaterial for Chapter 6: Basic Principles of Tomography M I A Integral Equations in Visual Computing Material
Material for Chapter : Integral Equations in Visual Computing Material Basic Principles of Tomography c 00 Bernhard Burgeth 0 Source: Images Figure : Radon Transform: ttenuation http://en.wikimedia.org/wiki/image:radon_transform.png
More informationFilters and Fourier Transforms
Filters and Fourier Transforms NOTE: Before reading these notes, see 15-462 Basic Raster notes: http://www.cs.cmu.edu/afs/cs/academic/class/15462/web/notes/notes.html OUTLINE: The Cost of Filtering Fourier
More informationLecture note on Senior Laboratory in Phyhsics. Fraunhofer diffraction and double-split experiment
Lecture note on Senior Laboratory in Phyhsics Fraunhofer diffraction and double-split experiment Masatsugu Suzuki and I.S. Suzuki Department of Physics State University of New York at Binghamton, Binghamton,
More informationTemperature Distribution Measurement Based on ML-EM Method Using Enclosed Acoustic CT System
Sensors & Transducers 2013 by IFSA http://www.sensorsportal.com Temperature Distribution Measurement Based on ML-EM Method Using Enclosed Acoustic CT System Shinji Ohyama, Masato Mukouyama Graduate School
More informationFirst CT Scanner. How it Works. Contemporary CT. Before and After CT. Computer Tomography: How It Works. Medical Imaging and Pattern Recognition
Computer Tomography: How t Works Medical maging and Pattern Recognition Lecture 7 Computed Tomography Oleh Tretiak Only one plane is illuminated. Source-subject motion provides added information. 2 How
More informationCamera Model and Calibration
Camera Model and Calibration Lecture-10 Camera Calibration Determine extrinsic and intrinsic parameters of camera Extrinsic 3D location and orientation of camera Intrinsic Focal length The size of the
More informationFourier transform of images
Fourier transform of images Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Methods for Image Processing academic year 2014 2015 Extension to bidimensional domain The concepts
More informationNon-Stationary CT Image Noise Spectrum Analysis
Non-Stationary CT Image Noise Spectrum Analysis Michael Balda, Björn J. Heismann,, Joachim Hornegger Pattern Recognition Lab, Friedrich-Alexander-Universität Erlangen Siemens Healthcare, Erlangen michael.balda@informatik.uni-erlangen.de
More informationDrawing a Triangle (and an introduction to sampling)
Lecture 4: Drawing a Triangle (and an introduction to sampling) Computer Graphics CMU 15-462/15-662, Spring 2017 Assignment 1 is out! https://15462-s17.github.io/asst1_drawsvg/ Let s draw some triangles
More informationPhys. 428, Lecture 6
Phys. 428, Lecture 6 Mid-Term & Weekly questions Since we are behind in the material, there will be no midterm Instead the final project will be subdivided into graded stages Class Project Pick: An imaging
More informationMRI Physics II: Gradients, Imaging
MRI Physics II: Gradients, Imaging Douglas C., Ph.D. Dept. of Biomedical Engineering University of Michigan, Ann Arbor Magnetic Fields in MRI B 0 The main magnetic field. Always on (0.5-7 T) Magnetizes
More informationBlind Sparse Motion MRI with Linear Subpixel Interpolation
Blind Sparse Motion MRI with Linear Subpixel Interpolation Anita Möller, Marco Maaß, Alfred Mertins Institute for Signal Processing, University of Lübeck moeller@isip.uni-luebeck.de Abstract. Vital and
More informationDigital Image Fundamentals. Prof. George Wolberg Dept. of Computer Science City College of New York
Digital Image Fundamentals Prof. George Wolberg Dept. of Computer Science City College of New York Objectives In this lecture we discuss: - Image acquisition - Sampling and quantization - Spatial and graylevel
More informationImage Acquisition + Histograms
Image Processing - Lesson 1 Image Acquisition + Histograms Image Characteristics Image Acquisition Image Digitization Sampling Quantization Histograms Histogram Equalization What is an Image? An image
More information2D Fan Beam Reconstruction 3D Cone Beam Reconstruction
2D Fan Beam Reconstruction 3D Cone Beam Reconstruction Mario Koerner March 17, 2006 1 2D Fan Beam Reconstruction Two-dimensional objects can be reconstructed from projections that were acquired using parallel
More informationTo Do. Advanced Computer Graphics. Sampling and Reconstruction. Outline. Sign up for Piazza
Advanced Computer Graphics CSE 63 [Spring 207], Lecture 3 Ravi Ramamoorthi http://www.cs.ucsd.edu/~ravir Sign up for Piazza To Do Assignment, Due Apr 28. Anyone need help finding partners? Any issues with
More informationLecture 6 Basic Signal Processing
Lecture 6 Basic Signal Processing Copyright c1996, 1997 by Pat Hanrahan Motivation Many aspects of computer graphics and computer imagery differ from aspects of conventional graphics and imagery because
More informationLecture 2 Image Processing and Filtering
Lecture 2 Image Processing and Filtering UW CSE vision faculty What s on our plate today? Image formation Image sampling and quantization Image interpolation Domain transformations Affine image transformations
More informationSolutions For Homework #7
Solutions For Homework #7 Problem :[ pts] Let f(r) = r = x2 + y 2 () We compute the Hankel Transform of f(r) by first computing its Abel Transform and then calculating the D Fourier Transform of the result.
More informationImportant!!! First homework is due on Monday, September 26 at 8:00 am.
Important!!! First homework is due on Monday, September 26 at 8:00 am. You can solve and submit the homework on line using webwork: http://webwork.dartmouth.edu/webwork2/m3cod/. If you do not have a user
More informationImage Reconstruction 3 Fully 3D Reconstruction
Image Reconstruction 3 Fully 3D Reconstruction Thomas Bortfeld Massachusetts General Hospital, Radiation Oncology, HMS HST.S14, February 25, 2013 Thomas Bortfeld (MGH, HMS, Rad. Onc.) Image Reconstruction
More information