1 / 25

Introduction to Tomography

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

cate
Télécharger la présentation

Introduction to Tomography

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Introduction to Tomography Presented by Scott Lichtor

  2. Overview • Problem Statement • Tomographic Applications • The Mathematics • Necessary Math • Fourier Slice Theorem • Filtered Backprojection • Matlab Example

  3. 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?

  4. 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

  5. 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

  6. 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

  7. 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

  8. 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

  9. 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

  10. 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

  11. 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

  12. 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

  13. 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 • 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 Image taken from http://www.eng.warwick.ac.uk

  14. 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

  15. 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)

  16. Example • Matlab illustration • Typical image to reconstruct:

  17. Example • Create projection data • Use the radon function • The radon function applies the radon transform to an image

  18. Example • 18 projections

  19. Example • 36 projections

  20. Example • 90 projections

  21. Example • The imradon function reconstructs images from projection data

  22. Example • Reconstruction with 36 projections

  23. Example • Reconstruction with 36 projections

  24. Example • Reconstruction with 90 projections

  25. 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

More Related