Folie 1 - uni

Download Report

Transcript Folie 1 - uni

Medical Imaging
Dr. Mohammad Dawood
Department of Computer Science
University of Münster
Germany
Medical Imaging, SS-2011
Recap
Mohammad Dawood
2
Medical Imaging, SS-2011
α decay
Gamma scintigraphy
Mohammad Dawood
β- decay
PET
β+ decay
Sinogram
3
Medical Imaging, SS-2011
Ultrasound / Sonography
Mohammad Dawood
4
Medical Imaging, SS-2011
Sound spectrum
The frequency range of sound above 20kHz is known as ultrasound
Mohammad Dawood
5
Medical Imaging, SS-2011
Sound wave propagation
Mohammad Dawood
6
Medical Imaging, SS-2011
Ultrasound
Mohammad Dawood
7
Medical Imaging, SS-2011
Ultrasound
Mohammad Dawood
8
Medical Imaging, SS-2011
Ultrasound
- is produced through the conversion of electrical energy into mechanical energy
- is detected by the reverse process, by converting mechanical energy into
electrical energy.
- The transducer is a device that is both a transmitter and receiver of the
ultrasound signal and it serves a dual role in pulse echo imaging.
Mohammad Dawood
9
Medical Imaging, SS-2011
Inge Edler und Helmut Hertz
Mohammad Dawood
(1954)
10
Medical Imaging, SS-2011
Ultrasound
Mohammad Dawood
11
Medical Imaging, SS-2011
Ultrasound
Reflexion
Mohammad Dawood
Scatter
Transmission
12
Medical Imaging, SS-2011
Ultrasound
Depth
Resolution
Frequency
Mohammad Dawood
13
Medical Imaging, SS-2011
Ultrasound
Mohammad Dawood
14
Medical Imaging, SS-2011
Ultrasound
Mohammad Dawood
15
Medical Imaging, SS-2011
Mohammad Dawood
Medium
Speed of sound m/s
Air
331
Water
1483
Tissue
1540-1595
Liver
1549
Blood
1570
Glycerin
1923
Bones (Compact)
3600
16
Medical Imaging, SS-2011
Mohammad Dawood
17
Medical Imaging, SS-2011
Image Reconstruction
Mohammad Dawood
18
Medical Imaging, SS-2011
Reconstruction
Tomography
detector
x-ray source
collimator
object
Mohammad Dawood
19
Medical Imaging, SS-2011
Reconstruction
Law of Attenuation
Mohammad Dawood
20
Medical Imaging, SS-2011
Reconstruction
Parallel projections of a plane
Mohammad Dawood
21
Medical Imaging, SS-2011
Reconstruction
y
s
Radon Transformation
f
n
r
θ
x
Mohammad Dawood
22
Medical Imaging, SS-2011
Reconstruction
Radon Transformation (Line Integrals at different angles)
Mohammad Dawood
23
Medical Imaging, SS-2011
Reconstruction
Radon Transformation
Original
Mohammad Dawood
Sinogram (Radon Transform)
24
Medical Imaging, SS-2011
Reconstruction
Inverse Radon Transformation
H: Hilbert transform
Mohammad Dawood
25
Medical Imaging, SS-2011
Reconstruction
Inverse Radon Transformation
Problems with Missing data and Noise!
Mohammad Dawood
26
Medical Imaging, SS-2011
Reconstruction
Filtered Back Projection
Mohammad Dawood
27
Medical Imaging, SS-2011
Mohammad Dawood
28
Medical Imaging, SS-2011
Reconstruction
Filtered Back Projection
Mohammad Dawood
29
Medical Imaging, SS-2011
Reconstruction
FBP: Commonly used filters
1=Ram-Lak (ramp), 2=Shepp-Logan,
3=Cosine, and 4=Hamming
Mohammad Dawood
30
Medical Imaging, SS-2011
Reconstruction
Filtered Back Projection
2D/3D filtering is costly
Backproject
Filter 2D
Projections
Image
Filter 1D
Mohammad Dawood
Backproject
31
Medical Imaging, SS-2011
Reconstruction
Fourier slice theorem
Take a two-dimensional function f(r), project it onto a line, and do a
Fourier transform of that projection
Take that same function, but do a two-dimensional Fourier transform
first, and then slice it through its origin parallel to the projection line
Mohammad Dawood
32
Medical Imaging, SS-2011
Reconstruction
Fourier Slice Theorem
Mohammad Dawood
33
Medical Imaging, SS-2011
Mohammad Dawood
34
Medical Imaging, SS-2011
Mohammad Dawood
35
Medical Imaging, SS-2011
Reconstruction
Iterative Reconstruction
Kaczmarz Method (=ART: Algebraic Reconstruction Technique)
Mohammad Dawood
37
Medical Imaging, SS-2011
Reconstruction
Iterative Reconstruction
Kaczmarz Method (=ART: Algebraic Reconstruction Technique)
1. Start by setting x(0) = 0
2. Compute the forward projection from the n-th estimate, i.e. b(n) = A x(n)
3. Choose i and correct the current estimate x(n)
4. Iterate steps 2,3 until the difference between new forward projection b(n),
computed in 2, and the old one is below tolerance
Mohammad Dawood
38
Medical Imaging, SS-2011
Mohammad Dawood
39
Medical Imaging, SS-2011
Reconstruction
Iterative Reconstruction
EM (Expectation Maximization)
Mohammad Dawood
40