filtered backprojection. radon transformation radon transform in 2-d. named after the austrian...

23
Filtered Backprojection

Upload: christina-preston

Post on 31-Dec-2015

234 views

Category:

Documents


2 download

TRANSCRIPT

Page 1: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

Filtered Backprojection

Page 2: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

Radon Transformation

• Radon transform in 2-D.• Named after the Austrian mathematician

Johann Radon• RT is the integral transform consisting of

the integral of a function over straight lines.• The inverse of RT is used to reconstruct

images from medical computed tomography scans.

Page 3: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

Paralle-beam Projection

• A projection of a 2-D image f(x,y) is a set of line integrals.

• To represent an image, RT takes multiple, parallel-beam projections of the image from different angles by rotating the source around the center of the image.

Page 4: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

• For instance, the line integral of f(x,y) in the vertical direction is the projection of f(x,y) onto the x-axis.

Page 5: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

Math. and Geometry of theRadon Transform

Page 6: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

Fourier Slice Theorem• FT of the projection of a 2-D object is

equal to a slice through the origin of 2-D FT of the object.

Page 7: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

Collection of projections of an object at a number of angles

• In Fourier domain

Page 8: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform
Page 9: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform
Page 10: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

1D Fourier Transform

Page 11: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

1D Fourier Transform

Page 12: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

2D Fourier Transform

Page 13: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

2-D Fourier Basis

Page 14: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

2-D FT Example

Page 15: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform
Page 16: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

Backprojection

• A problem with backprojection is the blurring (star-like artifacts) that happens in the reconstructed image.

Page 17: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

Filtered Backprojection

• To remove the blurring, an optimal way is to apply a high-pass filter to eliminate these artifacts.

• Thus combine backprojection with high-pass filtering = filtered backproejction.

Page 18: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

Approaches to Backprojection

Projections

BackprojectUnblurr with a 2-D Filter

Unblurr with 1 2-D Filter

Backproject

Image

Image

Page 19: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

Filtered Projection

• Fourier Slice Theorem

• filtered back projection takes the Fourier Slice and applies a weighting

Page 20: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

Digital Filters

Page 21: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

RT Example

Page 22: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

Reconstruction Example

unfiltered filtered

Page 23: Filtered Backprojection. Radon Transformation Radon transform in 2-D. Named after the Austrian mathematician Johann Radon RT is the integral transform

Filtered Backprojection Algorithm

• In Matlab, implemented as iradon.m

• 1-D FFT

• Digital Filters

• Interpolation Functions

• 2-D Inverse FFT