240 likes | 387 Views
Medical Imaging. Mohammad Dawood Department of Computer Science University of Münster Germany. Image Reconstruction. Reconstruction Law of Attenuation. Reconstruction Parallel projections of a plane. y. Reconstruction Radon Transformation. s. f. r. n. θ. x. Reconstruction
E N D
Medical Imaging Mohammad Dawood Department of Computer Science University of Münster Germany
Reconstruction Law of Attenuation
Reconstruction Parallel projections of a plane
y Reconstruction Radon Transformation s f r n θ x
Reconstruction Radon Transformation (Line Integrals at different angles)
Reconstruction Radon Transformation Original Sinogram (Radon Transform)
Reconstruction Inverse Radon Transformation H: Hilbert transform
Reconstruction Filtered Back Projection
Reconstruction Filtered Back Projection
Reconstruction Filtered Back Projection 2D/3D filtering is costly Backproject Filter 2D Projections Image Filter 1D Backproject
Reconstruction Fourier Slice Theorem
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
Reconstruction Fourier Slice Theorem
Reconstruction FBP: Commonly used filters 1=Ram-Lak (ramp), 2=Shepp-Logan, 3=Cosine, and 4=Hamming
Reconstruction Iterative Reconstruction b: measured values x: unknown attenuation coefficients aij: weights
Reconstruction Iterative Reconstruction Kaczmarz Method (=ART: Algebraic Reconstruction Technique)
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
Reconstruction Iterative Reconstruction EM (Expectation Maximization)
Reconstruction Iterative Reconstruction OSEM (Ordered Subset Expectation Maximization)