PM 6 4 Filtered Backprojection
PM 6 4 Filtered Backprojection
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.
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.
For instance, the line integral of f(x,y) in
the vertical direction is the projection of
f(x,y) onto the x-axis.
Math. and Geometry of the
Radon 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.
Collection of projections of an
object at a number of angles
In Fourier domain
1D Fourier Transform
1D Fourier Transform
2D Fourier Transform
2-D Fourier Basis
2-D FT Example
Backprojection
A problem with backprojection is the
blurring (star-like artifacts) that happens in
the reconstructed image.
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.
Approaches to Backprojection
Projections
Backproject
Unblurr with
a 2-D Filter
Unblurr with
1 2-D Filter
Backproject
Image
Image
Filtered Projection
Fourier Slice Theorem
filtered back projection takes the Fourier
Slice and applies a weighting
Digital Filters
RT Example
Reconstruction Example
unfiltered filtered
Filtered Backprojection Algorithm
In Matlab, implemented as iradon.m
1-D FFT
Digital Filters
Interpolation Functions
2-D Inverse FFT