prof. dr. philippecattin: computed tomography contents ... · prof. dr. philippecattin: computed...

33
Computed Tomography Principles of Medical Imaging Prof. Dr. Philippe Cattin MIAC, University of Basel Sep 26th/Oct 3rd, 2016 Sep 26th/Oct 3rd, 2016 Principles of Medical Imaging 1 of 65 26.09.2016 08:34 Contents 2 4 5 6 7 8 9 11 12 13 14 15 16 17 18 20 21 22 23 24 25 26 27 Prof. Dr. Philippe Cattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography Hounsfield's CT Prototype EMI-Scanner Detectors Important Terminology 2 Single Slice CT First Generation CT Scanner Design Second Generation CT Scanner Design Third Generation CT Scanner Design Fourth Generation CT Scanner Design Spiral Scanning CT Spiral Scanning CT (2) Spiral Scanning CT (3) Drawback of these Designs 3 Multi-Detector Row CT Multi-Detector Row CT Detector Design Detector Design (2) Detector Design (3) Detector Design (4) Detector Design (5) Detector Design (6) Dual Source CT Sep 26th/Oct 3rd, 2016 Principles of Medical Imaging 2 of 65 26.09.2016 08:34

Upload: others

Post on 22-Sep-2020

18 views

Category:

Documents


3 download

TRANSCRIPT

Page 1: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

ComputedTomography

Principles of MedicalImaging

Prof. Dr. Philippe Cattin

MIAC, University of Basel

Sep 26th/Oct 3rd, 2016

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

1 of 65 26.09.2016 08:34

Contents

2

4

5

6

7

8

9

11

12

13

14

15

16

17

18

20

21

22

23

24

25

26

27

Prof. Dr. Philippe Cattin: Computed Tomography

Contents

Abstract

1 Computed Tomography Basics

Introduction

Computed Tomography

Hounsfield's CT Prototype

EMI-Scanner

Detectors

Important Terminology

2 Single Slice CT

First Generation CT Scanner Design

Second Generation CT Scanner Design

Third Generation CT Scanner Design

Fourth Generation CT Scanner Design

Spiral Scanning CT

Spiral Scanning CT (2)

Spiral Scanning CT (3)

Drawback of these Designs

3 Multi-Detector Row CT

Multi-Detector Row CT

Detector Design

Detector Design (2)

Detector Design (3)

Detector Design (4)

Detector Design (5)

Detector Design (6)

Dual Source CT Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

2 of 65 26.09.2016 08:34

Page 2: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

28

29

31

32

33

36

37

39

40

41

42

43

44

46

47

48

49

51

52

53

54

56

Open Dual Source CT

Advantage of the DSCT

4 Electron Beam Tomography

Electron Beam Computed Tomography

Electron Beam Computed Tomography (2)

Dosage Comparison: EBCT vs. CT

5 Image Reconstruction

5.1 Introduction

Image Reconstruction

Image Reconstruction (2)

5.2 Radon Transform

Radon Transform

Parallel Projection

The Radon Transform

The Discrete Radon Transform

Radon Transform Examples

Radon Transform Examples (2)

5.3 Fourier Slice Theorem

Fourier Slice Theorem

Fourier Slice Theorem (2)

Reconstruction with the Fourier Slice Theorem

Reconstruction with the Fourier Slice Theorem(2)

5.4 Filtered Backprojection

Principle of Filtered Back-Projection

Numerical Back-Projection Example

Example Reconstructions

Example Reconstructions (2)

5.5 Helical Reconstruction

Helical Reconstruction Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

3 of 65 26.09.2016 08:34

57

58

60

62

63

65

66

67

68

69

71

360° Linear Interpolation

180° Linear Interpolation

5.6 Hounsfield Unit

Hounsfield Unit

6 Automatic Exposure Control

Automatic Exposure Control

Automatic Exposure Control (2)

7 Artefacts

Artefacts

Partial Volume Effect

High Density Artefacts

Gating in Cardio CT

CT and Medical Image Analysis

8 X-Ray Dosage Summary

X-Ray Dosage Summary

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

4 of 65 26.09.2016 08:34

Page 3: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

(2)

Prof. Dr. Philippe Cattin: Computed Tomography

Abstract

5 of 65 26.09.2016 08:34

ComputedTomography Basics

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

(4)Introduction

One of the major disadvantages associated with conventionalplanar radiography is its inability to produce sectional information.

The images produced on film represent the total attenuation of theX-ray beam as it passes through the patient. Depth information iscompletely lost!

Two general classes of tomography exist that solve this problem:

Linear tomography, which produces longitudinal sections

Computed axial tomography, which produces sectional or axial

slices

6 of 65 26.09.2016 08:34

Page 4: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Computed Tomography Basics

(5)

Prof. Dr. Philippe Cattin: Computed Tomography

Computed Tomography

→ Computed Tomography (CT)[http://en.wikipedia.org

/wiki/Computed_axial_tomography] originallyknown as Computed Axial Tomography(CAT) or Body Section Röntgenographyis a medical imaging modality used togenerate 3D images of the internals ofan object from a large series of 2DX-ray images taken around a single axisof rotation.

→ Godfrey Newbold Hounsfield[http://en.wikipedia.org

/wiki/Godfrey_Newbold_Hounsfield] conceivedthe CT scanner idea in 1967 andpublicly announced it in 1972. → AllanMcLeod Cormack [http://en.wikipedia.org

/wiki/Allan_McLeod_Cormack] independentlyinvented a similar process and theyshared the Nobel price in 1979.

Fig 3.1: CT Apparatus

It is claimed that the CT scanner was the greatest legacy ofthe Beatles; the massive profits from their record salesenabled EMI to fund scientific research

7 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Computed Tomography Basics

(6)

Prof. Dr. Philippe Cattin: Computed Tomography

Hounsfield's CT Prototype

The original 1971 prototype took parallel readings through angles, each apart, with each scan taking a little over fiveminutes. The images from these scans took hours to beprocessed by algebraic reconstruction techniques on a largecomputer.

Fig 3.2: Hounsfield's original CT prototype Fig 3.3: Principle of theprototype

8 of 65 26.09.2016 08:34

Page 5: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Computed Tomography Basics

(7)

Prof. Dr. Philippe Cattin: Computed Tomography

EMI-Scanner

The EMI-Scanner was the firstproduction X-ray CT machine. It waslimited to scan two adjacent slices ofthe brain, but acquired the image datain about . The computation timewas about per picture.

The scanner required the use of awater-filled Perspex tank with apre-shaped rubber head-cap at thefront. The water-tank was used toreduce the dynamic range of theradiation reaching the detectors(scanning outside the head vs. throughthe skull).

The images were relatively lowresolution, being composed of a matrixof only .

The CT scanner was a huge success: by1977 1130 machines were installedacross the world.

Fig 3.4: EMI brain scanner witha Data General Nova

minicomputer. The first scannerwas installed at Atkinson

Morley's Hospital, Wimbledon,England in 1971

9 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Computed Tomography Basics

(8)

Prof. Dr. Philippe Cattin: Computed Tomography

Detectors

→ Scintillator [http://en.wikipedia.org/wiki/Scintillator] Detectors

Low maximum count rate leads to longer scan times or more imagenoise

Xenon Gas Detectors

Pressurised Xe gas capable of higher count rates, but low detectionefficiency

Modern Ceramic → Scintillators [http://en.wikipedia.org/wiki/Scintillator]

Coupled with photodiodes these detectors offer the bestperformance

10 of 65 26.09.2016 08:34

Page 6: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Computed Tomography Basics

(9)

Prof. Dr. Philippe Cattin: Computed Tomography

Important Terminology

In-plane resolution:

acquisition resolution in the

-plane

Out-of-plane, through-plane

resolution: slice distance in

axis

Anisotropic scan: the

resolution in the axis is

generally less than in the

axis

Isotropic scan: the voxel

dimensions are equal in the

, and axis

Fig 3.5: Coordinate system generallyused

11 of 65 26.09.2016 08:34

Single Slice CT

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

(11)First Generation CT ScannerDesign

The generation of CT scanner usedthe translate-rotate geometry.

The EMI scanner, for instance, used apencil X-ray beam and a single detector.During translation of the gantry, theX-ray beam was sampled 160 times.After a rotation of a new profile wasacquired. This procedure was repeatedfor 180 different angles and tookroughly .

To minimise patient movement the headwas usually clamped.

Fig 3.6: First generation CTprinciple

12 of 65 26.09.2016 08:34

Page 7: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Single Slice CT

(12)

Prof. Dr. Philippe Cattin: Computed Tomography

Second Generation CTScanner Design

The generation scanner tried toreduce the excessive scan times byusing a small fan beam with multipledetectors (up to 30 in some designs).

Scan times of between werepossible with this design.

The introduction of multiple detectorswas an important development.

Fig 3.7: Second generation CTprinciple

13 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Single Slice CT

(13)

Prof. Dr. Philippe Cattin: Computed Tomography

Third Generation CT ScannerDesign

The generation brought down scantimes even further by using the rotate-rotate geometry.

As the large fan beam encompasses thepatient completely the translatorymotion of the previous designs can beavoided. The X-ray tube and thedetector array rotate as one about thepatient.

The number of detector elements istypically in the hundreds.

To avoid excessive variations in signalstrength various manufacturers use abow-tie shaped filter to suit the body orhead shape.

Fig 3.8: Third generation CTprinciple

14 of 65 26.09.2016 08:34

Page 8: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Single Slice CT

(14)

Prof. Dr. Philippe Cattin: Computed Tomography

Fourth Generation CTScanner Design

The generation CT uses arotate-fixed ring geometry where thering of detectors completely surroundsthe patient.

As the X-ray tube must be closer to thepatient than the detectors it has a poorradiographic geometry, i.e. largegeometric magnification.

Scan times as low as withinterscan delays of can beachieved with this type of geometry.

Using many thousand detectorelements a in-plane resolution of

can be obtained.

Fig 3.9: Fourth generation CTprinciple

15 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Single Slice CT

(15)

Prof. Dr. Philippe Cattin: Computed Tomography

Spiral Scanning CT

Advances in slip-ring technologyhave enabled the X-ray tube torotate continuously in the samedirection which overcomesproblems of interscan delays.

If the continuous motion of thegantry is combined with acontinuous advance of thepatient table along thelongitudinal axis we have aspiral/helical scanner.

The spiral scanning technologybrought about a significantreduction in scan times.

The gained speed came at aprice of increased complexity forreconstructing the helical data.

Fig 3.10: Illustration of helical scanning

Fig 3.11: Nice 3D rendering of helicalCT

16 of 65 26.09.2016 08:34

Page 9: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Single Slice CT

(16)

Prof. Dr. Philippe Cattin: Computed Tomography

Spiral Scanning CT (2)

The X-ray source iscollimated to a fan beamrotating around thepatient.

The X-ray tube and thedetectors are fixedtogether as a singlerotating unit.

Post patient collimationdefines the slicesensitivity profile.

Fig. 3.12: Basic design of a single slice CT usedin a spiral CT

17 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Single Slice CT

(17)

Prof. Dr. Philippe Cattin: Computed Tomography

Spiral Scanning CT (3)

In the context of helical scanning aparameter called Pitch is defined as the

Ratio of the distance that thepatient couch moves in onerotation to the collimationthickness (number of slices slice thickness)

(3.1)

In other words, for a couch advance of and a nominal collimation width

of , the pitch is 1. Pitch valuesare typically in the range of 1 to 2depending on the required spatialresolution in the direction of the couchmotion. Its a coverage indicator, inother words.

Fig 3.13: Pitch

18 of 65 26.09.2016 08:34

Page 10: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Single Slice CT

(18)

Prof. Dr. Philippe Cattin: Computed Tomography

Drawback of these Designs

Ideally, volume data are of high isotropic spatial resolution, haveminimal motion artefacts, and optimally utilise the contrast agentbolus.

To reduce motion artefacts CT examinations need to be completedwithin a certain time frame, e.g. on breath hold, forthe heart.

If, however, a large scan range such as the entire thorax has to becovered

a thick collimation (large inter slice distance)

must be used, leading to anisotropic voxel sizes (whilst the in-planeresolution only depends on the system geometry).

19 of 65 26.09.2016 08:34

Multi-Detector RowCT

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

(20)Multi-Detector Row CT

Strategies to achieve a betterlongitudinal resolution andfaster scans include thesimultaneous acquisition ofmultiple slices at a time, thustermed Multi-Detector Row CTor Multi-Slice CT (MSCT).

Interestingly, the very firstcommercial CT systems(EMI-Scanner and SiemensSiretom) were already two-slicesystems. Only the introduction ofthe helical scanning principleallowed to fully leverage theadvantages of multi-detector rowCT.

Fig 3.14: Multi-slice CT

Fig 3.15: SOMATOM Sensation 16Gantry (Siemens)

20 of 65 26.09.2016 08:34

Page 11: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Multi-Detector Row CT

(21)

Prof. Dr. Philippe Cattin: Computed Tomography

Detector Design

The figure shows, how different slice widths can be achieved byprepatient collimation for a single slice detector .

Fig 3.16: Prepatient collimation of the X-ray beam to obtain different slicethicknesses with a single detector row CT.

21 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Multi-Detector Row CT

(22)

Prof. Dr. Philippe Cattin: Computed Tomography

Detector Design (2)

The principle can be easily extended to slices if the sensor isseparated midway along the axis.

Fig 3.17: Collimation of the X-ray beam to obtain different slice thicknesses with atwo detector row CT.

For detectors a more elaborate detector design is required.

22 of 65 26.09.2016 08:34

Page 12: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Multi-Detector Row CT

(23)

Prof. Dr. Philippe Cattin: Computed Tomography

Detector Design (3)

The various manufacturers introduceddifferent detector designs in order toallow utmost flexibility in selecting slicewidths.

All designs combine several detectorrows electronically to a smaller numberof slices according to the selected slicewidth.

The total coverage of this detectordesign is (measured in theisocenter).

With prepatient collimation thefollowing slice widths can be realised:

, , , and .

Fig 3.18: Fixed array detector,16 rows, 4 slices

23 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Multi-Detector Row CT

(24)

Prof. Dr. Philippe Cattin: Computed Tomography

Detector Design (4)

A more efficient approach (needs lessdetector channels) uses the adaptivearray design.

This design allows the followingcollimated slice widths: two slices at

, four at , four at ,two at , and two at .

Fig 3.19: Adaptive arraydetector, 8 rows, 4 slices

24 of 65 26.09.2016 08:34

Page 13: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Multi-Detector Row CT

(25)

Prof. Dr. Philippe Cattin: Computed Tomography

Detector Design (5)

Sixteen-slice CT systems usually haveadaptive array detectors similar to theone depicted in Fig 3.20. It uses 24detector rows with a total coverage of

at the isocenter.

By properly combining the detectorrows, either 12 or 16 slices with

or can be acquiredsimultaneously.

Fig 3.20: Adaptive arraydetector, 24 rows, 16 slices

(Siemens)

25 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Multi-Detector Row CT

(26)

Prof. Dr. Philippe Cattin: Computed Tomography

Detector Design (6)

32, 40, and 64 slice systems are now available.

Fig 3.21: Toshiba detector mock-ups

26 of 65 26.09.2016 08:34

Page 14: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Multi-Detector Row CT

(27)

Prof. Dr. Philippe Cattin: Computed Tomography

Dual Source CT

A different approach to acquire moreslices in parallel was followed bySiemens with their → Dual SourceCT [http://www.siemens.com/dualsource]

(SOMATOM Definition).

Fig 3.22: Dual Source CT (SiemensSOMATOM Definition)

Fig 3.23: Comparison of LAD & Cxin diastole and systole

27 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Multi-Detector Row CT

(28)

Prof. Dr. Philippe Cattin: Computed Tomography

Open Dual Source CT

28 of 65 26.09.2016 08:34

Page 15: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Fig 3.24: Movie of the an open rotating dual source CT

29 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Multi-Detector Row CT

(29)

Prof. Dr. Philippe Cattin: Computed Tomography

Advantage of the DSCT

The scan is in cardiac-mode virtually independent of the heart

rate → no -blocker needed

If the two X-Ray tubes are operated with two different tube

voltages (other spectra) tissue types can be better

differentiated

Fig 3.25: HU values for different tissue types (theoretical simulation)

30 of 65 26.09.2016 08:34

Page 16: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Electron BeamTomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

31 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

(31)Electron Beam ComputedTomography

→ Electron beamcomputed tomography(EBCT or EBT)[http://en.wikipedia.org

/wiki/Electron_beam_tomography]

is a specific form of CTscanner in which theX-Ray tube is notmechanically spun inorder to rotate the sourceof X-Ray photons. Thisdifferent design wasexplicitly developed tobetter image heartstructures which neverstop moving.

As in conventional CTtechnology, the X-raysource still rotates aroundthe circle in spacecontaining an object to beimaged tomographically,but the X-Ray tube ismuch larger than theimaging circle and theelectron beam currentwithin the vacuum tube isswept electronically, in acircular path and focusedon a stationary tungstenanode target ring.

Fig 3.26: Patent illustration showing a cutawayview of an electron beam CT system.

Components are 22. electron gun, 23. electronbeam, 27. beam bending coil, 28-31. target

rings, 14. detector array. The electron beam isreflected by the target rings through the patient,to the detector on the opposite end of the scan

tube.

32 of 65 26.09.2016 08:34

Page 17: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Electron Beam Tomography

(32)

Prof. Dr. Philippe Cattin: Computed Tomography

Electron Beam ComputedTomography (2)

Design Advantages

The most advanced commercial EBCT can

perform image sweeps in as little as

compared to of the

mechanically swept X-Ray tube designs.

Design Disadvantages

Larger footprint and smaller scanning

volume

More than twice as expensive

Higher demands on room shielding from

electro-magnetical interferences

Sensitive to small vibrations

Fig 3.27: GE eSpeed300

Fig 3.28: Imatron C150

33 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Electron Beam Tomography

(33)

Prof. Dr. Philippe Cattin: Computed Tomography

Dosage Comparison: EBCTvs. CT

As a general rule, the electron beam CT delivers → onlyabout of the radiation [http://www.newportbodyscan.com

/radiationdosage.htm] than a conventional CT scanner would.

The primary explanation being that the EBCT is a fast shutteredcamera that only turns on the beam as needed to acquire theimage.

Conventional CT scanners have their X-ray emitter always on(often modulated) during the acquisition.

34 of 65 26.09.2016 08:34

Page 18: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

ImageReconstruction

Introduction

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

(36)Image Reconstruction

From the scanning process we have a set of image projections.Given these projections we want to determine the X-rayattenuation coefficients of the original image as accurate aspossible.

Fig 3.29: Image projections

Fig 3.30: A small section of the finalmatrix showing individual attenuation

values combined as a ray-sum

35 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Introduction

(37)

Prof. Dr. Philippe Cattin: Computed Tomography

Image Reconstruction (2)

Already in 1917 → Johann Radon [http://en.wikipedia.org/wiki/Johann_Radon]

published a paper with the mathematical theory, the → Radontransform [http://en.wikipedia.org/wiki/Radon_transform], useful toreconstruct a 2D image from multiple projections such as in CTsystems.

Hounsfield used, for the first CT scanner, an iterative technique toexactly solve the Radon transform. Its disadvantages are that it isslow and that all data must be collected before reconstruction canbegin.

Todays CT systems mainly use variants of the filteredback-projection approach that is computationally more efficient.

36 of 65 26.09.2016 08:34

Page 19: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Radon Transform

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

(39)Radon Transform

A straight line in Cartesiancoordinates can be either describedby its slope-intercept form

(3.2)

or by its normal representation

(3.3)

see Fig 3.31.

Fig. 3.31: Different linerepresentations

37 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Radon Transform

(40)

Prof. Dr. Philippe Cattin: Computed Tomography

Parallel Projection

An arbitrary point in the projection is given by the raysum alongthe line

(3.4)

in the continuous space the raysum is then given by

(3.5)

where is the impulse

function.

(3.6)

with .

The integrand is zerounless the argument inthe delta function is

zero. This is valid for allpoints on the line

.

Fig. 3.32: Projection geometry

38 of 65 26.09.2016 08:34

Page 20: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Radon Transform

(41)

Prof. Dr. Philippe Cattin: Computed Tomography

The Radon Transform

We can generalise this equation to arbitrary lines

(3.7)

This projection is called Radon transform. Often used notations forthe Radon transform of are

(3.8)

39 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Radon Transform

(42)

Prof. Dr. Philippe Cattin: Computed Tomography

The Discrete RadonTransform

In the discrete case the integrals in the Radon transform arereplaced by sums

(3.9)

The Radon transform forms the corner stone ofreconstruction from projections used e.g. in ComputedTomography, PET, SPECT.

40 of 65 26.09.2016 08:34

Page 21: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Radon Transform

(43)

Prof. Dr. Philippe Cattin: Computed Tomography

Radon Transform Examples

The figure below shows an example image with its Radontransform. The interpretation of the sinogram is still quite easy.

Fig. 3.33: Double box image

Fig. 3.34: Sinogram of the double boximage with projections over

41 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Radon Transform

(44)

Prof. Dr. Philippe Cattin: Computed Tomography

Radon Transform Examples(2)

The figure below shows an example phantom with its Radontransform. The interpretation of the sinogram is not possibleanymore, although the phantom's structure is quite simple.

Fig. 3.35: Shepp-Logan phantom

Fig. 3.36: Sinogram of the phantomwith projections over

42 of 65 26.09.2016 08:34

Page 22: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Fourier SliceTheorem

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

(46)Fourier Slice Theorem

In the following slide we will relate the Fourier transform of 1-Dprojection with the 2-D Fourier transform of the scanned object.Without loss of generality, we take the projection line to be the -axis in the derivation below. Given is the image and its

projection onto the -axis where

(3.10)

The Fourier transform of is

(3.11)

the slice at is then

(3.12)

which is the Fourier transform of .

43 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Fourier Slice Theorem

(47)

Prof. Dr. Philippe Cattin: Computed Tomography

Fourier Slice Theorem (2)

Fig. 3.37: Graphical representation of the Fourier slice theorem

The Fourier slice theorem states that, the 1-dimensionalFourier transform of a projection corresponds to the slice(line) - at the same angle - in the 2-dimensional Fouriertransform of the object

44 of 65 26.09.2016 08:34

Page 23: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Fourier Slice Theorem

(48)

Prof. Dr. Philippe Cattin: Computed Tomography

Reconstruction with theFourier Slice Theorem

In principle we could reconstruct the image by filling up the

Fourier space with the Fourier transforms of the individual

projections and then calculate the inverse Fourier transform. Thisapproach is, however, computationally very expensive.

Fig. 3.38: Image reconstructing using the Fourier slice theorem

45 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Fourier Slice Theorem

(49)

Prof. Dr. Philippe Cattin: Computed Tomography

Reconstruction with theFourier Slice Theorem (2)

If we just sum the spectra ofthe individual projectionbeams, the spectral densityfor low frequencies would betoo high as the beams arecloser to each other for smallradii → lower frequencies toostrong.

We therefore must correctthe spectrum with a suitableweighting factor. As thedensity of the projectionbeam goes with

(Frequency) the spectra mustbe multiplied with → ramp

filter.

Each projection directionthus has to be multiplied witha suitable weighting function

. As will be seen in the

next section, this can also beperformed as a convolutionwith the inverse Fouriertransform of in the

spatial domain → filteredback-projection.

Fig. 3.39: Spectral density must becorrected to get suitable results

46 of 65 26.09.2016 08:34

Page 24: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

FilteredBackprojection

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

(51)Principle of FilteredBack-Projection

(1) The measured projections are smeared back, i.e. combined as aray-sum, across the output matrix. (2) As the back-projected imageis heavily blurred and shows star artefacts it has to be filtered witha highpass yielding the final reconstructed image.

Fig 3.40: Filtered back-projection principle

47 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Filtered Backprojection

(52)

Prof. Dr. Philippe Cattin: Computed Tomography

Numerical Back-ProjectionExample

48 of 65 26.09.2016 08:34

Page 25: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Filtered Backprojection

(53)

Prof. Dr. Philippe Cattin: Computed Tomography

Example Reconstructions

Back-projectionwithout filteringusing 2, 4, 8,16, and 32projections.

Strong artefactscan be seen andthe images areheavily blurred.

Fig 3.41: Back-projection without filter

49 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Filtered Backprojection

(54)

Prof. Dr. Philippe Cattin: Computed Tomography

Example Reconstructions (2)

Back-projectionwith highpassfiltering usingthe same 2, 4, 8,16, and 32projections.

Strong artefactscan be seen andthe images areheavily blurred.

Fig 3.42: Back-projection with filter

50 of 65 26.09.2016 08:34

Page 26: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

HelicalReconstruction

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

(56)Helical Reconstruction

We would like to use the same filtered back-projection method asbefore:

Choose the

interpolation

position along

the z-axis

Only one

projection is

from the

reconstruction

position, others

are from

different

z-positions

Fig 3.43: Problem of the helical reconstruction

→ Introduces artefacts when the structures change along thez-axis.

51 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Helical Reconstruction

(57)

Prof. Dr. Philippe Cattin: Computed Tomography

360° Linear Interpolation

Idea: Use attenuation data from points apart on the helix forinterpolation

Fig 3.44: Interpolation

→ Interpolation makes the effective image width broader.

52 of 65 26.09.2016 08:34

Page 27: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Helical Reconstruction

(58)

Prof. Dr. Philippe Cattin: Computed Tomography

180° Linear Interpolation

Idea: Use attenuation data from complementary projections inaddition to points apart on the helix

Fig 3.45:Complementary

projections

Fig 3.46: Interpolation

→ Slice profile is narrower, as the z-axis distances are shorter thanin interpolation.

53 of 65 26.09.2016 08:34

Hounsfield Unit

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

(60)Hounsfield Unit

The tissue absorption coefficient depends on the tube voltage.

To make them comparable, theabsorption coefficients have tobe related to that of water a thesame tube voltage. This way anumber [Hounsfield unit = Hu]insensitive to tube voltage canbe obtained:

(3.13)

In practice CT values areproduced from for air,

for water, and between for bone.

Fig 3.47: Common CT numbers

54 of 65 26.09.2016 08:34

Page 28: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Automatic ExposureControl

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

55 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

(62)Automatic Exposure Control

To minimise radiation exposure and optimise radiation efficiencyautomatic exposure control (AEC) methods have been integratedinto todays CT systems. Several different fields to control the tubecurrent (mA) are explored:

Relative patient attenuationDepending on the patients body weight(diameter) tube current has to beadjusted.

Rotational attenuationAs the total X-ray attenuation is smallerin the AP-direction than in theLR-direction, tube current can bedynamically adapted.

z-Axis attenuationTube current is automatically adaptedto the patient geometry reducing theradiation exposure. In this example the6 year old boy was irradiated with onaverage instead of inthe standard protocol.

Heart cycle gatingFor heart and aorta scans heart-gatingto synchronise the acquisition with thepulse is required.

Tube power is generally reduced to in uninteresting parts of the heart

cycle.

56 of 65 26.09.2016 08:34

Page 29: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

57 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Automatic Exposure Control

(63)

Prof. Dr. Philippe Cattin: Computed Tomography

Automatic Exposure Control(2)

Allows the scanner to be controlled by setting a required image

quality level

Can avoid over- and under-exposure of patients and thus

minimises rescans

58 of 65 26.09.2016 08:34

Page 30: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Artefacts

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

(65)Artefacts

Several inherent CT artefacts have an important influence on theapplied Medical Image Analysis Methods and generally need to beaccounted for:

Partial Volume Effect

High density artefacts

Gating in Cardio CT

The Good News: CT data is geometrically very accurate. IfMR/US data is to be registered with CT, then CT should beused as the reference!

59 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Artefacts

(66)

Prof. Dr. Philippe Cattin: Computed Tomography

Partial Volume Effect

The partial volume effect,common to most medicalimaging modalities, poses animportant problem for manymedical image analysis methods.

The sampling of the imagingvolume renders itdifficult/impossible to exactlylocate the boundary of an object.

If not taken special care of, asimple shift of the object candrastically change the result,e.g. area measurement in Fig3.48(c)+(d). The measured areaof (d) is higher than thearea of (c).

Fig 3.48: (a) Original object, (b) objectsampled on a discrete grid, (c)

thresholded object, (d) thresholdedshifted object

60 of 65 26.09.2016 08:34

Page 31: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Artefacts

(67)

Prof. Dr. Philippe Cattin: Computed Tomography

High Density Artefacts

High density streak artefacts or Windmill artefacts result from thefinite width of the detector rows, which require interpolation. Theartefacts appear close to high contrast gradients.

Fig 3.49: High density artefact example

Fig 3.50: Windmill effect

61 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Artefacts

(68)

Prof. Dr. Philippe Cattin: Computed Tomography

Gating in Cardio CT

For many acquisitions of the heart and arterial system ECG gatingis used. To reduce exposure, the AEC reduces the tube current to

during systole (when no images are captured). Proper gating,however, depends on a regular sine rhythm not present in alldiseased patients.

Fig 3.51: Image of good quality withdecent SNR

Fig 3.52: Gating failed, image wasacquired with of the dose → very

noisy image

62 of 65 26.09.2016 08:34

Page 32: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

Artefacts

(69)

Prof. Dr. Philippe Cattin: Computed Tomography

CT and Medical ImageAnalysis

Medical image analysis is an indispensable tool in CT. Without theaid of advanced image analysis methods, radiologists would needsubstantially more time to find interesting location in the CTdatasets.

Fig 3.53: Plaque detection aid

Fig 3.54: Automatic dissection

63 of 65 26.09.2016 08:34

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

segmentation

64 of 65 26.09.2016 08:34

Page 33: Prof. Dr. PhilippeCattin: Computed Tomography Contents ... · Prof. Dr. PhilippeCattin: Computed Tomography Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography

X-Ray DosageSummary

Sep 26th/Oct 3rd, 2016Principles of Medical Imaging

(71)X-Ray Dosage Summary

ExaminationX-Ray

Dosage[mSv]

CT Dosage[mSv]

Ratio

Skull 0.03 2.6 87

Thorax 0.3 20.5 68

Abdomen 1.17 27.4 23

Spine 3.3 9.0 3

Extremities 0.06 1.0 17

The average X-ray exposition in central Europe isapproximately

65 of 65 26.09.2016 08:34