Introduction to Tomography - Engineering School Class Web Sites
Download
Report
Transcript Introduction to Tomography - Engineering School Class Web Sites
Introduction to Tomography
Presented by Scott Lichtor
Overview
Problem Statement
Tomographic Applications
The Mathematics
Necessary Math
Fourier Slice Theorem
Filtered Backprojection
Matlab Example
Problem
Can’t see inside of people to diagnose problems
Can’t see inside of machinery to diagnose problems
How do take a picture of a place where you can’t fit a
camera?
Solution
Tomography
Reconstructs a function using line integrals
Goal: recover the interior structure of a body using exterior
measurements
Routine for medicine, earth sciences
Image taken from http://media-2.web.britannica.com
Tomography Applications
Single photon emission computed tomography (SPECT) is used
for gamma imaging
Gamma-emitting radio-isotope is injected into the body
Gamma camera returns a 2-D image of the object
Reconstruction then returns a 3-D image of the object
Used for medical imaging (tumor imaging, functional brain
imaging)
Image taken from http://www.biocompresearch.org
Tomography Applications
Positron emission tomography (PET) acquires data from electron
positron annihilation
Positron-emitting tracer is injected into the body
System detects gamma rays produced by tracer
Uses PET to reconstruct 3-D image
Used for oncology, neurology, cardiology, etc.
Image taken from http://www.ibfm.cnr.it
Tomography Applications
Computed tomography (CT) is used for X-ray imaging
X-rays are produced and sent through the body
Record the line integrals
Calculate the shape of the imaged object
Used extensively for medical imaging
Also used for non-destructive materials testing
Image taken from http://www.csmc.edu
Tomography
I’ll focus on X-ray tomography
Get interior structure of body by X-raying the object from
many different directions
When an X-ray goes through an object, it is attenuated by the
object
Very dense objects will weaken the strength of the ray
considerable
Less dense objects will affect the strength of the ray less
History of Computed Tomography
Alessandro Vallebona proposed representing a slice of the
body on radiographic film in the early 1900s
First commercially viable CT scanner invented by Sir
Godfrey Hounsfield at EMI Laboratories in 1972
Originally, water tanks were needed for imaging on humans
Necessary Mathematics
Line integrals are integrals along a line
Coordinate system: (x,y)->(Ѳ,t)
Ѳ: Angle, t: distance along source
Fourier Transform: F(w) = ∫f(t)e-j2πwtdt
F(u,v)=∫∫f(x,y)e-j2π(ux+vy)dxdy
Image taken from http://www.mindef.gov.sg
Tomography
A projection is composed of a bunch of line integrals
Easiest example: line integrals with the same Ѳ but different
t’s (parallel line integrals).
The value of a line integral:
P Ѳ(t) = ∫(Ѳ,t)line f(x,y)ds
P Ѳ(t) = ∫∫f(x,y)δ(x cos (Ѳ)+y sin (Ѳ)-t)dxdy
Radon transform
Fourier Slice Theorem
Object function is f
Fourier transform of f is F
Projection P
Fourier transform of P is S
F(u,v)=∫∫f(x,y)e-j2π(ux+vy)dxdy
S Ѳ (w) = ∫P Ѳ (t)e-j2πwtdt
To demonstrate the Fourier Slice Theorem, let Ѳ=0
Fourier Slice Theorem
Suppose v=0
F(u,0) = ∫∫f(x,y)e-j2πuxdxdy
•
= ∫(∫f(x,y)dy)e-j2πuxdx
P Ѳ=0(x) = ∫ f(x,y)dy
So F(u,0) = ∫ P Ѳ=0(x) e-j2πuxdx Image taken from http://www.eng.warwick.ac.uk
There’s a relationship between the projection data and the object
image
Specifically, each projection gives a slice of the Fourier transform
of the overall image
Filtered Backprojection
Filtered backprojection is the algorithm used to reconstruct the
object image
Idea: use the projection data to get slices of the Fourier transform
of the object image. Then, calculate the object image
Image taken from http://www.eng.warwick.ac.uk
Filtered Backprojection
Procedure:
For all angles K
1.
Get projections P Ѳ
2.
Apply Fourier transform and get S Ѳ (w)
3.
Place the inverse Fourier transforms of the
projections on the approximation of the original image
In this way an approximation of the original image can be
obtained (this is only the algorithm for parallel projections)
Example
Matlab illustration
Typical image to reconstruct:
Example
Create projection data
Use the radon function
The radon function applies the radon transform to an image
Example
18 projections
Example
36 projections
Example
90 projections
Example
The imradon function reconstructs images from projection
data
Example
Reconstruction with 36 projections
Example
Reconstruction with 36 projections
Example
Reconstruction with 90 projections
Sources
An Introduction to X-ray tomography and Radon
Transforms, by Eric Todd Quinto
Principles of Computerized Tomographic Imaging, by Avinash
C. Kak and Malcolm Slaney
Wikipedia