- Home
- Documents
*Computed Tomography (Part 1)eeweb.poly.edu/~yao/EL6813/Lect4_CT_ch6_part1_F18_JM.pdf Computed...*

prev

next

out of 58

View

8Download

0

Embed Size (px)

Computed Tomography (Part 1)

Jonathan Mamou and Yao Wang Tandon School of Engineering

New York University, Brooklyn, NY 11201

Based on Prince and Links, Medical Imaging Signals and Systems and Lecture Notes by Prince. Figures are from the book.

EL-GY 6813 / BE-GY 6203 / G16.4426 Medical Imaging

F16 Yao Wang @ NYU

Lecture Outline • Instrumentation

– CT Generations – X-ray source and collimation – CT detectors

• Image Formation – Line integrals – Parallel Ray Reconstruction

• Radon transform • Back projection • Filtered backprojection • Convolution backprojection • Implementation issues

2

F16 Yao Wang @ NYU

Limitation of Projection Radiography • Projection radiography

– Projection of a 2D slice along one direction only – Can only see the “shadow” of the 3D body

• CT: generating many 1D projections in different angles – When angle spacing is sufficiently small and angular coverage is

complete, 2D slices can be reconstructed successfully

3

F16 Yao Wang @ NYU

1st Generation CT: Parallel Projections

4

F16 Yao Wang @ NYU

2nd Generation

5

F16 Yao Wang @ NYU

3G: Fan Beam

Much faster than 2G

6

F16 Yao Wang @ NYU

4G

Fast Cannot use collimator at detector, hence affected by scattering

7

F16 Yao Wang @ NYU

5G: Electron Beam CT (EBCT)

Stationary source and detector. Used for fast (cine) whole heart imaging Source of x-ray moves around by steering an electron beam around X-ray tube anode.

8

F16 Yao Wang @ NYU

6G: Helical (or spiral) CT

Entire abdomen or chest can be completed in 30 sec.

9

F16 Yao Wang @ NYU

7G: Multislice

From http://www.kau.edu.sa/Files/0008512/Files/19500_2nd_presentation_final.pdf

10

F16 Yao Wang @ NYU

Reduced scan time and increased Z-resolution (thin slices) Most modern MSCT systems generates 64 slices per rotation, can image whole body (1.5 m) in 30 sec.

From http://www.kau.edu.sa/Files/0008512/Files/19500_2nd_presentation_final.pdf

11

F16 Yao Wang @ NYU

From http://www.kau.edu.sa/Files/0008512/Files/19500_2nd_presentation_final.pdf

12

F16 Yao Wang @ NYU

From http://www.kau.edu.sa/Files/0008512/Files/19500_2nd_presentation_final.pdf

13

F16 Yao Wang @ NYU

X-ray Source

14

Dual Energy CT • Recall that linear attenuation coefficient depends on the

X-ray energy • Dual energy CT enables simultanious acquisition to two

images of the organ, reflecting their attenuation profiles corresponding to different energies

• The scanner has two x-ray tubes using two different voltage levels to generate the x-rays.

F16 Yao Wang @ NYU 15

F16 Yao Wang @ NYU

X-ray Detectors

Convert detected photons to lights

Convert light to electric current

16

Solid-state detector Xenon gas detector Multiple solid-state detector array

F16 Yao Wang @ NYU

CT Measurement Model

17

F16 Yao Wang @ NYU 18

I0 must be calibrated for each detector, measured in the absence of attenuating objects

F16 Yao Wang @ NYU

CT Number

Requires at least 12 bits for representation (e.g., h=-1000 for air)

Godfrey Hounsfield, together with Alan Cormack, invented X-ray CT!

19

CT numbers of Different Tissues at 70KeV [Smith&Webb]

F16 Yao Wang @ NYU 20

F16 Yao Wang @ NYU

Parameterization of a Line

s

x

y

Option 1 (parameterized by s):

Option 2:

Each projection line L is defined by (l,θ)

A point on this line (x,y) can be specified with two options

21

F16 Yao Wang @ NYU

Line Integral: parametric form

22

F16 Yao Wang @ NYU

Line Integral: set form

23

F16 Yao Wang @ NYU

Physical meaning of “f” & “g”

24

F16 Yao Wang @ NYU

What is g(l,θ)?

25

F16 Yao Wang @ NYU

Example • Example 1: Consider an image slice which contains a

single square in the center. What is its projections along 0, 45, 90, 135 degrees?

• Example 2: Instead of a square, we have a rectangle. Repeat.

26

F16 Yao Wang @ NYU

Sinogram

27

Bottom row: θ=0 (vertical proj), middle θ =π/2 (horizontal proj), top θ = π -Δθ

F16 Yao Wang @ NYU

Backprojection • The simplest method for reconstructing an image from a

projection along an angle is by backprojection – Assigning every point in the image along the line defined by (l,θ) the

projected value g(l, θ), repeat for all l for the given θ

s

x

y

28

F16 Yao Wang @ NYU 29

F16 Yao Wang @ NYU

Example • Continue with the example of the image with a square in

the center. Determine the backprojected image from each projection and the reconstruction by summing different number of backprojections

30

F16 Yao Wang @ NYU

Two Ways of Performing Backprojection • Option 1: assigning value of g(l,θ) to all points on the line (l, θ)

– g(l, θ) is only measured at certain l: ln=n ∆l – If l is coarsely sampled (∆l is large), many points in an image will not be

assigned a value – Many points on the line may not be a sample point in a digital image

• Option 2: For each θ, go through all sampling points (x,y) in an image, find its corresponding “l=x cos θ +y sin θ”, take the g value for (l, θ)

– g(l, θ) is only measured at certain l: ln=n Δl – must interpolate g(l,θ) for any l from given g(ln,θ)

• Option 2 is better, as it makes sure all sample points in an image are assigned a value

• For more accurate results, the backprojected value at each point should be divided by the length of the underlying image in the projection direction (if known)

31

F16 Yao Wang @ NYU

Backprojection Summation

Replaced by a sum in practice

32

F16 Yao Wang @ NYU

Implementation Issues

From L. Parra at CUNY, http://bme.ccny.cuny.edu/faculty/parra/teaching/med-imaging/lecture4.pdf 33

g(:,phi) stores the projection data at angle phi corresponding to excentricities stored in s_n

F16 Yao Wang @ NYU

Implementation: Projection • Creating projection data using computers will yield similar problems.

Possible l and θ are both quantized. If you first specify (l,θ), then find (x,y) that are on this line. It is not easy. Instead, for given θ, you can go through all (x,y) and determine corresponding l, quantize l to one of those you want to collect data.

• Sample matlab code (for illustration purpose only) – f(x,y) stores the original image data – G(l,phi) stores projection data, ql is the desired quantization stepsize for l.

N=ceil(sqrt(I*I+J*J))+1; %(assume image size IxJ), N0 is maximum lateral distance N0= floor((N-1)/2); ql=1; G=zeros(N,180); for phi=0:179

for (x=-J/2:J/2-1; y=-I/2:I/2-1) l=x*cos(phi*pi/180)+y*sin(phi*pi/180); l=round(l/ql)+N0+1; If (l>=1) && (l

F16 Yao Wang @ NYU

Problems with Backprojection Blurring

35

F16 Yao Wang @ NYU

Projection Slice Theorem

The 1D Fourier Transform of a projection at angle θ is a line in the 2D Fourier transform of the image at the same angle passing through the origin.

If (l, θ) are sampled sufficiently dense, then from g (l,θ) we essentially know F(u,v) (on the polar coordinate), and by inverse transform we can obtain f(x,y)!

dlljlgG }2exp{),(),( πρθθρ −= ∫ ∞

∞−

36

F16 Yao Wang @ NYU

Illustration of the Projection Slice Theorem

37

F16 Yao Wang @ NYU

Proof • Go through on the board • Using the set form of the line integral • See Prince&Links, P. 203

dlljlgG }2exp{),(),( πρθθρ −= ∫ ∞

∞−

38

F16 Yao Wang @ NYU

The Fourier Method • The projection slice theorem leads to the following

conceptually simple reconstruction method – Take 1D FT of each projection g(l,θ) to obtain G(ρ, θ) for all θ – Convert G(ρ, θ) to Cartesian grid F(u,v) – Take inverse 2D FT to obtain f(x,y)

• Not used because – Difficult to interpolate polar data onto a Cartesian grid – Inverse 2D FT is time consuming

• But is important for conceptual understanding – Take inverse 2D FT of G(ρ, θ) in the polar coordinate grid leads

to the widely used Filtered Backprojection algorithm

39

F16 Yao Wang @ NYU

Filtered Backprojection • Inverse 2D FT in Cartesian coordinate:

• Inverse 2D FT in Pol