Radiography & Tomography

Radiography & Tomography

Radiography & Tomography Federica Marone Swiss Light Source, Paul Scherrer Institut, Villigen, Switzerland Monday, May 15, 2017 Federica Marone – HSC19 Contents § Mini-tutorial on X-ray tomography § Beer-Lambert law § Tomographic reconstruction algorithms § Analytical algorithms § Iterative algorithms § Examples § Practical issues § Flat field correction § Local tomography Monday, May 15, 2017 Federica Marone – HSC19 X-ray imaging yesterday…and today! First X-ray image, 1896 X-ray micro-radiography of a fruit fly Monday, May 15, 2017 Federica Marone – HSC19 X-ray sources of the 21st century Röntgen’s Lab, Late 19th century Swiss Light Source, Today Monday, May 15, 2017 Federica Marone – HSC19 Why a synchrotron for imaging ? dqz 1000-500 µm Hospital 50-10 µm Table Top micro CT quality2 Fµ res42× contr Bonse et al., Prog. Biophys. Molec. Biol., Vol. 65, No.1, pp. 133-169, 1996 1-0.1 µm SR - Tomographic Microscopy Monday, May 15, 2017 Federica Marone – HSC19 Mini-tutorial on X-ray tomography X-ray interaction with matter § Photoelectric effect § Elastic scattering Photoelectric Pb § Compton scattering effect Elastic scattering Compton scattering Source: NIST, Bertrand et al, 2012 Monday, May 15, 2017 Federica Marone – HSC19 X-ray interaction with matter § Photoelectric effect �� Cross section � ∝ �� �� Monday, May 15, 2017 Federica Marone – HSC19 Beer-Lambert law ./0 � = �+ , � � � = �� + = �� � � 1 � � = = �� + � � � • Monochromatic radiation μ = linear attenuation coefficient • Homogeneous object � ∝ � = �89 + �;< + �=< + … Monday, May 15, 2017 Federica Marone – HSC19 Slicing Imaging – Why ? § Tomography means imaging by sections or slices. (From the Greek word tomos, meaning "a section" or "a cut”) Radiographic projection Tomographic slice Monday, May 15, 2017 Federica Marone – HSC19 Beer-Lambert law ./0 � = �+ , � � � = �� + = �� � � 1 � � = = �� + � � � • Monochromatic radiation μ = linear attenuation coefficient • Homogeneous object � = �89 + �;< + �=< + … Monday, May 15, 2017 Federica Marone – HSC19 Beer-Lambert law ABCD I . ∫ / A,H 0H � = @ �+ � , � J �� + � � = ? I . ∫ / H 0H � = �+ , � J � 0 � = �� + = @ � � �� � + • Polychromatic radiation • Inhomogeneous object Radon Transform Monday, May 15, 2017 Federica Marone – HSC19 Beer-Lambert law T ./MNO ./PNO ./RNO . ∑ /RNO � = �+ , � , � ⋯ � = �+ , � RUM Z � ? � = �� + = V � Δ� � W W[\ Monday, May 15, 2017 Federica Marone – HSC19 Tomographic reconstruction algorithms § Radon Transform § Fourier Slice Theorem § Analytical algorithms § Direct Fourier Methods § Filtered Backprojection § Iterative algorithms § Algebraic § Statistical Monday, May 15, 2017 Federica Marone – HSC19 Radon transform § The Radon transform in 2D is the integral of a function over straight lines and therefore represents the projections data as obtained in a tomographic scan s R(t,φ) = ∫ µ (x,y)dl l y " I (t)% φ t R(t,φ) = Pφ (t) := −ln$ ' # I0 & +∞ +∞ R t,φ = µ x,y δ(x cosφ + ysinφ − t)dx dy x ( ) ∫ ∫ ( ) φ −∞ −∞ Wanted ! t μ(x,y) Monday, May 15, 2017 Federica Marone – HSC19 Radon transform § The Radon transform was introduced by Radon, who also provided a formula for the inverse transform § The inverse of the Radon transform can be used to reconstruct the original density from the projection data, and thus it forms the mathematical underpinning for tomographic reconstruction § Radon found the solution to the tomographic problem already in 1917 § but assumes an infinite number of projections and continuous projection functions § while we only have a finite number of projections and a finite number of detector points. +∞ +∞ R(t,φ) = ∫ ∫ µ (x,y)δ(x cosφ + ysinφ − t)dx dy −∞ −∞ Wanted ! Monday, May 15, 2017 Federica Marone – HSC19 Radon transform § The Radon transform is often called a sinogram because the Radon transform of a Dirac delta function is a distribution supported on the graph of a sine wave Monday, May 15, 2017 Federica Marone – HSC19 Fourier Slice Theorem The Fourier transform of a parallel projection of an object f(x,y) obtained at an angle θ equals a line of the 2D Fourier transform of f(x,y) taken at the same angle. Monday, May 15, 2017 Federica Marone – HSC19 Direct Fourier Methods § 1D FFT of the projections § Filtering § Resampling § 2D inverse FFT Monday, May 15, 2017 Federica Marone – HSC19 Direct Fourier Methods - resampling § Interpolation § Limited accuracy § Gridding § Most accurate Fourier reconstruction method § Mapping by convolution with FFT of w(x,y) § Gridrec § W(x): 1D Prolate Spheroidal Wave Functions (PSWF) of zeroth order § W(x,y): maximally concentrated in a square region of interest § FFT(w(x,y)): concentrated as much as possible around 0 § PSWF: calculated using rapidly converging expansion in terms of Legendre polynomials § PSWF and FFT (PSWF) can be efficiently computed and stored at run time Dowd et al., 1999; Marone et al., 2010 Monday, May 15, 2017 Federica Marone – HSC19 Filtered backprojection (mathematical background) +∞ +∞ " Iφ (r)% −ln$ ' := Pφ (r) = R(t,φ) = ∫ ∫ µ (x,y)δ(x cosφ + ysinφ − t)dx dy # I0 & −∞ −∞ ∞ ∞ Image function: µ(x, y) = ∫ ∫ M! (u,v)e j2π (ux+vy) du dv 2π−∞ −∞∞ Coord. transform: µ(x, y) = ∫ dθ ∫ M! (ω cosθ,ω sinθ)e j2πω ( xcosθ+ ysinθ ) g du��dv Cartesian to Polar 0 0 ⎪⎧ u = ω cosθ with ⎨ and ⎩⎪ v = ω sinθ 2π ∞ Fourier Slice Theorem: µ(x, y) = ∫ dθ ∫ P!(ω,θ)e j2πω ( xcosθ+ ysinθ )ω dω 0 0 Symmetry properties: P!(ω,θ + π ) = P!(−ω,θ) π ∞ j2 t Image function: µ(x, y) = ∫ dθ ∫ P!(ω,θ) ω e πω dω 0 −∞ A. C. Kak, M. Slaney, Principle of Computerized Tomographic Imaging”, SIAM Classics in Applied Mathematics 33, New York, 2001, ISBN 0-89871-494-X Monday, May 15, 2017 Federica Marone – HSC19 Reconstruction using backprojection 0-22.5° 0-45° § If no a priori information is known, the intensity of the object is assumed to be uniform along the beam path. § The projection intensity is evenly 0-67.5° 0-90° 0-112° distributed among all pixels along the ray path à Concept of backprojection! 0-135° 0-157° 0-180° Backprojection of a point … Monday, May 15, 2017 Federica Marone – HSC19 π ∞ Image reconstruction µ(x, y) = ∫ dθ ∫ P!(ω,θ) ω e j2πωt dω 0 −∞ without pre-filtering with filtering Monday, May 15, 2017 Federica Marone – HSC19 Filtered backprojection § Image characteristic can be influenced by the choice of a convolution kernel, whereby increasing spatial resolution or edge enhancement also means increasing image noise ! Monday, May 15, 2017 Federica Marone – HSC19 Analytical algorithms - Summary ü Simple and fast ✖ Need large number of projections § In the order of the number of detector rows ✖ Cannot include a-priori information Ø Not suitable for under-sampled datasets Monday, May 15, 2017 Federica Marone – HSC19 Direct inversion (historic) § Brute force approach § N equations with N unknowns § Need more than N2 equations to ensure linear independence § Sir Godfrey Hounsfield (in the late 60s) § Reconstructed the first human head § Solved 28000 equations simultaneously § Nobel Prize in 1979 together with A. Cormack Monday, May 15, 2017 Federica Marone – HSC19 Iterative algorithms Radon transform Guessed solution Guessed data ill-posed �� = � Data difference SMALL Guessed solution GOOD Monday, May 15, 2017 Federica Marone – HSC19 Iterative algorithms ? ? ? ? Algebraic ReconstructionTechnique (ART) Monday, May 15, 2017 Federica Marone – HSC19 Iterative algorithms ill-posed �� = � Optimization problem Cost function argmin � �, � + �(�) / Fidelity Regularization term term Monday, May 15, 2017 Federica Marone – HSC19 Iterative algorithms Optimization problem Cost function argmin � �, � + �(�) / Fidelity Regularization • Stability • A priori term term knowledge Algebraic methods 1 �(�, �) = �� − � n � � = 0 (ART,SIRT,…) 2 n n Tikhonov – small norm � � = � � n Lasso – sparsity � � = � � \ Total Variation (TV) – piecewise constant � � = � �� \ Preserve edges Dictionary, Nonlocal means, Nonlocal TV, … Monday, May 15, 2017 Federica Marone – HSC19 Iterative algorithms N – expected photon counts s � . �� � �� = �� , � i – each pixel in each view Optimization problem Cost function argmin � �, � + �(�) / Fidelity Regularization term term Monday, May 15, 2017 Federica Marone – HSC19 Iterative algorithms Optimization problem Cost function argmin � �, � + �(�) / Fidelity Regularization • Stability • A priori term term knowledge Statistical methods } Negative log-likelihood �(�|�) = V �z log �sz − �sz z Many cost functions (model geometries, artifacts, …) Many possibilities Monday, May 15, 2017 Federica Marone – HSC19 Iterative algorithms § FIRST STEP - Build the cost function § SECOND STEP – Find minimum/maximum § Optimization techniques § Linear/non-linear § Least-squares § Convex/non-convex § Constrained/unconstrained § Gradient methods § (Gauss-)Newton methods § Lagrangian methods § Expectation-maximization algorithms Monday, May 15, 2017 Federica Marone – HSC19 Iterative algorithms - Summary ü Can include a-priori information ü Flexible – can model almost anything ü Suitable for under-sampled datasets ✖ Computational intensity ✖ Highly dataset specific ✖ Parameters to tune Monday, May 15, 2017 Federica Marone – HSC19 Examples State-of-the-art SRXTM (1-50 um) 1 micron resolution routinely achieved at 10% MTF Monday, May 15, 2017 Federica Marone – HSC19 Gathering complex information Measure X-ray projections Reconstruct virtual slices Monday, May 15, 2017 Federica Marone – HSC19 Looking into a very small region of the lung Morphology of lung acini Haberthür et al., Journal of Synchrotron Radiation, 17(5), 2010 Schittny et al., American Journal of Physiolgy

View Full Text

Details

  • File Type
    pdf
  • Upload Time
    -
  • Content Languages
    English
  • Upload User
    Anonymous/Not logged-in
  • File Pages
    44 Page
  • File Size
    -

Download

Channel Download Status
Express Download Enable

Copyright

We respect the copyrights and intellectual property rights of all users. All uploaded documents are either original works of the uploader or authorized works of the rightful owners.

  • Not to be reproduced or distributed without explicit permission.
  • Not used for commercial purposes outside of approved use cases.
  • Not used to infringe on the rights of the original creators.
  • If you believe any content infringes your copyright, please contact us immediately.

Support

For help with questions, suggestions, or problems, please contact us