Computed Tomography (Part 2)
|
|
- Bertina Morton
- 5 years ago
- Views:
Transcription
1 EL-GY 6813 / BE-GY 603 / G Medical Imaging Computed Tomography (Part ) Yao Wang Polytechnic School of Engineering New York University, Brooklyn, NY 1101 Based on Prince and Links, Medical Imaging Signals and Systems and Lecture Notes by Prince. Figures are from the book.
2 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 Last Lecture F16 Yao NYU
3 This Lecture Review of Parallel Ray Projection and Reconstruction Practical implementation with samples Fan Beam Reconstruction Signal to Noise in CT F16 Yao NYU 3
4 Review: Projection Slice Theorem Projection Slice theorem The Fourier Transform of a projection at angle θ is a line in the Fourier transform of the image at the same angle. The image cannot be displayed. Your computer may not have enough memory to open the image, or the image may have been corrupted. Restart your computer, and then open the file again. If the red x still appears, you may have to delete the image and then insert it again. F16 Yao NYU 4
5 Reconstruction Algorithm for Parallel Projections Backprojection Sum: Backprojection of each projection Sum Filtered backprojection: FT of each projection Filtering each projection in frequency domain Inverse FT Backprojection Sum Convolution backprojection Convolve each projection with the ramp filter Backprojection Sum Which one gives more accurate result? Which one takes less computation? F16 Yao NYU 5
6 Practical Implementation Projections g(l, q) are only measured at finite intervals l=nt; t chosen based on maximum frequency in G(ρ, θ), W 1/t >=W or t <=1/W (Nyquist Sampling Theorem) W can be estimated by the number of cycles/cm in the projection direction in the most detailed area in the slice to be scanned For filtered backprojection: Fourier transform G(ρ, θ) is obtained via FFT using samples g(nt, q) If N sample are taken, ideally N point FFT should be taken by zero padding g(nt, q) Recall convolving two signals of length N leads to a single signal of length N-1 For convolution backprojection The ramp-filter is sampled at l=nt Sampled Ram-Lak Filter (limiting the rho-filter to within (-1/t, -1/t ) c( n) = 1/ 4τ ; 1/ ( nπτ ) 0; ; n = 0 n = odd n = even See paper at F16 Yao NYU 6
7 Practical Implementation of Filtered P θ j (l):!projection!data!along!θ j ; R(m,n)= π K! Backprojection s θ j (k):!dft!of!p θ j (l); B θ j (m,n):!backprojected!image!from!filtered!projection R(m,n):Sum!of!backprojections. L:!number!of!samples!along!each!projection,!assuming!L!is!odd K :!number!of!projection!angles. s θ j (k)= L 1 P θ (l)exp iπ lk j L, k = L 1,...,L 1,L 1 = floor L l= L 1 B θ j (m,n)= L 1 k= L 1 K j=1 s θ (k) k j L exp iπ mcosθ + nθ j j k L B θ j (m,n) From paper at Note this implementation perform L point DFT assuming the projection data has L points. L is odd. The reconstructed values near the boundary may have ringing artifacts. Note that backprojection is implemented within the evaluation of inverse FT. F16 Yao NYU 7
8 Practical Implementation of Convolution Backprojection c(l):!predesigned!filter B θ j (m,n)= L 1 l= L 1 ( ) P θ (l) c mcosθ j j + nsinθ j l R(m,n)= π! K K j=1 B θ j (m,n) Here again: Back projection is implemented within convolution! F16 Yao NYU 8
9 The Ram-Lak Filter (from [Kak&Slaney]) F16 Yao NYU 9
10 1 st Generation CT: Parallel Projections F16 Yao NYU 10
11 3G: Fan Beam Much faster than G F16 Yao NYU 11
12 Fan Beam: Equiangular Ray We will focus on the equiangular detector setting on the right in this lecture. F16 Yao NYU 1
13 Equiangular Ray Projection Source location is described by ( β, D) D is typically fixed, β varies to provide a largeview angle. To provide complete view, β (0,π ) a For a given source with angle β, γ specifies the detector position or the projection line. For each β, γ varies over a range( γ, γ m m ) (D, β, γ ) completely specifies the line of projection : θ = β + γ, l = D sin( γ ) Instead of g(l, θ ), we can use p( γ, β ) to represent a projection β + α + γ = π / θ + α = π / β + γ = θ F16 Yao NYU 13
14 Equiangular Ray Reconstruction g D r φ β φ Reconstructed image is represented in the polar coordinate using (r,ϕ). The relative position of a pixel at (r,ϕ) to the source at (D, β) is specified by (D',γ) : ( ) + ( r cos(β ϕ) ) ( r cos(β ϕ) ) ( D + rsin(β ϕ) ) D' (r,ϕ) = D + rsin(β ϕ) tanγ(r,ϕ) = Reconstruction formular: f (r,ϕ) = π 0 1 q(γ, β)dβ ( D' ) q(γ, β) = p'(γ, β)*c f (γ) p'(γ, β) = p(γ, β)cos(γ) c f (γ) = 1 D # γ & % ( $ sinγ ' Convolution weighted backprojection c(γ) c(γ) is the ramp filter used in parallel projection. c f (γ) is weighted based on the angle γ. Derivation not required for this class. Detail can be found at [Kak&Slaney]. Note typos in [Prince&Links], 1 st ed. F16 Yao NYU 14
15 Typos in [Prince&Links, 1 st Ed] P. 07, Eq. (6.38), change to c( D'sin γ ) = γ D'sin γ c( γ ) Eq. (6.39) change to c f (γ) = 1 D! γ $ # & " sinγ % c(γ) Eq. (6.40),(6.41) p(γ, β) p'(γ, β) p'(γ, β) = p(γ, β)cos(γ) q(γ, β) = p'(γ, β)*c f (γ) F16 Yao NYU 15
16 Practical Implementation Projections P(γ,β) are only measured at finite intervals γ=nα; α chosen based on maximum frequency in γ direction, W 1/α>=W or α<=1/w For convolution backprojection The filter c f (γ) is sampled at γ=nα Sampled Filter g(nα) q(nα,β) = p (nα,β) *g(nα) For backprojection For given (r,φ), for a given β, determine (D,γ) Use interpolation to determine q(γ,β) from known values at γ=nα Modify the projection based on D (weighted backprojection) Back projection and sum D' -1/ ( r, φ) = tanγ ( r, φ) = ( D + r sin( β φ) ) + ( r cos( β φ) ) ( r cos( β φ) ) ( D + r sin( β φ )) f (r,ϕ) = π 0 1 q(γ, β)dβ ( D' ) q(γ, β) = p'(γ, β)*c f (γ) p'(γ, β) = p(γ, β)cos(γ) F16 Yao NYU 16
17 Matlab Functions for Fan Beam CT Relevant functions: fanbeam(), ifanbeam() F16 Yao NYU 17
18 Fan Beam Reconstruction Through Rebinning [Smith&Webb] F16 Yao NYU 18
19 Constrained Optimization Formulation (Algebraic Reconstruction Techniques or ART) The!measurement!g(l,θ)!in!general!can!be!related!to!the!unknown!f (m,n)!as g(l,θ) = m,n f (m,n)a(m,n; l,θ) 1 (m,n)!!is!on!the!line!(l,θ) A(m,n;l,θ) = 0 otherwise More!generally,!we!can!order!all!possible!(l,θ)!pairs!into!1D!indices,!and!represent g(l,θ)!!as!a!vector!g,!and!order!all!pixel!positions!(m,n)also!into!!1d!indices,!so!that! all!possible!f (m,n)!can!be!put!into!a!vector!f,!and!all!possible!a(m,n;!l,θ)!can!be!put! into!a!matrix!a.!!then!the!above!equations!can!generally!be!written!as! Af = g So!we!can!reconstruct!the!image!f!by!solving!the!above!matrix!equation. If!there!are!more!measurements!than!the!number!of!pixels,!then!the!problem!can!be!solved using!the!least!squares!formulation,!i.e., min!! Af g,!with!a!unique!solution!f = A T A ( ) '1 A T g When!there!are!fewer!measurements!than!the!number!of!pixels,!we!have!an!under!determined! problem,!and!we!need!to!put!additional!constraints!on!f.!!for!example,!we!may!require!that!the!image!f is!sparse!in!some!transform!domain,!i.e.!!tf!has!only!few!nonizero!coefficients,!or!l0!norm!of!tf!is!small,! and!solve!the!following!problem min!! Af g + λ Tf 0 This!is!the!classifical!compressive!sensing!problem.! Instead!of!using!the!L0!norm,!we!can!relax!to!the!L1!norm,!and!solve!the!following!convex!optimization!probem min!! Af g + λ Tf 1! There!are!many!fast!algorithms!to!solve!the!above!problem. F16 Yao NYU 19
20 Blurring Effect SNR CT Quality Evaluation F16 Yao NYU 0
21 Effect of Area Detector Practical detector integrates the detected photons over an area Mathematically, the detector can be characterized by an indicator function s(l) (aka impulse response) The measured projection g (l,q) is related to real projection g(l,q) by g (l,q)= g(l,q) * s(l) G (r,q)= G(r,q) S(r) F16 Yao NYU 1
22 Windowing Function Recall that the ideal filter c(r) is typically modified by a window function W(r) Overall Effect fˆ( x, y) can be thought of as the reconstructedimagefromthe projection gˆ( l, θ), whose Fourier transform is Gˆ( ρ, θ) = G( ρ, θ) S( ρ) W ( ρ) gˆ( l, θ) = g( l, θ)* s( l)* w( l) F16 Yao NYU
23 Blurred Projection h(x,y): PSF of the blurring F16 Yao NYU 3
24 Circular Symmetry of Blurring F16 Yao NYU 4
25 PSF given by Hankel Transform r = x + y F16 Yao NYU 5
26 Circularly Symmetric Functions and Hankel Transform Circularly symmetric: f(x,y) = f(r), only depends on the distance to the origin, not angle Fourier transform of circularly symmetric function is also circularly symmetric F(u,v)=F(ρ) Let x = r cosφ, y = r sin φ; u = ρ cosθ, v = ρ sinθ F( ρ, θ ) = = f ( r, φ)exp{ jπ ( rρ cosφ cosθ + rρ sin φ sinθ ) f ( r, φ)exp{ jπrρ cos( φ θ )} rdrdφ } rdrdφ If f ( x, y) = F( ρ, θ ) = { exp{ jπrρ cos( φ θ )} dφ} f ( r) f ( r) rdr = π 0 f ( r) J 0 (πρr) rdr = F( ρ) Hankel Transform F( ρ) = π 0 f ( r) J 0 (πρr) rdr; J 0 1 ( r) = π π 0 cos( r sin φ) dφ F16 Yao NYU 6
27 Common Transform pairs See Table.3 Fourier => Hankel π ( x + y ) { e } πr { e } = = e e π ( u πρ + v ) Hankel { sin c( r) } = rect( q) π 1 4q Scaling property 1 Hankel a a { f(ar) } = F( q / ) Duality: If h(r) <-> H(ρ), then H(r)<->h(ρ) Derivation of Hankel transform pairs are not required. But you should be able to use given transform pairs, to determine the blur function. F16 Yao NYU 7
28 Example Example 6.5 in [Prince&Links] Detector: rectangular detector with width d S(l)=rect(l/d) Rectangular window function Solution W(ρ)=rect(ρ/ρ o ); ρ o >>1/d S(l)=rect(l/d) <-> S(ρ)=d sinc(dρ) ρo >>1/d -> H(ρ)=S(ρ) W(ρ)~= S(ρ) =dsinc(dρ) (Hankel transform of h(r)) Illustrate and explain h(r) h( r) = inverse Hankel Hankel h( r) = { sinc( r) } rect ( r / d) 4r { H ( ρ) } rect ( ρ) 1 Hankel{ f(ar) } = F( ρ / a) a Using duality property π d = π 1 4ρ F16 Yao NYU 8
29 Noise in CT Measurement Pr{ N ij = k} = k a k! e a ; k = 0,1,... E{ N ij = k} = a Var{ N ij = k} = a F16 Yao NYU 9
30 What about the measured projection F16 Yao NYU 30
31 CBP Approximation F16 Yao NYU 31
32 Definitions and Assumptions g ij are independent because N ij are independent Deriving mean and variance of µ(x,y) based on the independence assumption See [Prince&Links] for derivation F16 Yao NYU 3
33 Variance increases with ρ 0 (cut-off freq. of filter), and T (detector spacing), decreases with M (number of angles),\bar N (or N 0 ) (x-ray intensity) F16 Yao NYU 33
34 SNR of the Reconstructed Image C: fractional change of m from \bar m F16 Yao NYU 34
35 SNR in a good design F16 Yao NYU 35
36 SNR in Fan Beam N=Nf/D w=l/d SNR decreases as D increases. Reason: Convolution of the projection reading with the ramp filter couples the noise between detectors, and effectively increases the noise as the number of detector increases But larger D is desired to obtain a good resolution. F16 Yao NYU 36
37 Rule of Thumb F16 Yao NYU 37
38 Aliasing Artifacts Nyquist Sampling theorem: If the maximum freq of a signal is fmax, it should be sampled with a freq fs>=max, or sampling interval T<=1/fmax If sampled at a lower freq. without pre-filtering, aliasing will occur High freq. content fold over to low freq Prefilter to lower fmax, and then sample If the number of samples in each projection (D) or the number of projection angles (M) are not sufficiently dense, the reconstructed image will have streak artifacts Caused by aliasing Practical detectors are area detectors and perform pre-filtering implicitly F16 Yao NYU 38
39 From [Kak&Slaney] Fig. 5.1 F16 Yao NYU 39
40 Summary Parallel projection reconstruction Backprojection summation Fourier method (projection slice theorem) Filtered backprojection Convolution backprojection Practical implementation: using finite samples Fan beam projection and reconstruction Weighted backprojection Blurring due to non-ideal filters and detectors Approximate the overall effect by a filter: h(l)=w(l)*s(l); H(r)=W(r) S(r) Circularly symmetric functions and Hankel transform Equivalent spatial domain filter h(r)=inverse Hankel {H(q)} Noise in measurement and reconstructed image Factors influencing the SNR of reconstructed image Average X-ray intensity, Number of angles (M), number of samples per angle (D), filter cut-off r o Impact of number of projection angles and samples on reconstruction image quality Nyquist sampling theorem Streak artifacts F16 Yao NYU 40
41 Reference Prince and Links, Medical Imaging Signals and Systems, Chap 6. A. C. Kak and M. Slaney, Principles of Computerized Tomographic Imaging. Originally published by IEEE, E-copy available at Chap 3 Contain detailed derivation of reconstruction algorithms both for parallel and fan beam projections. Have discussions both in continuous domain and implementation with sampled discrete signals. Chap 5 discusses noise in measurement and reconstructed image. Chap 5 also covers aliasing effect with more mathematical interpretations A useful lecture note from Prof. Fessler Lecture note by Prof. Parra For more discussion on aliasing due to under sampling F16 Yao NYU 41
42 Reading: Homework Prince and Links, Medical Imaging Signals and Systems, Chap 6, Sec Note down all the corrections for Ch. 6 on your copy of the textbook based on the provided errata. Problems for Chap 6 of the text book: P.6.9 P.6.10 (part e is not required) P Hint: solution for part (a) should be P.6.19 P.6.0 3µ g( l,60) = 3µ 0 ( a / + l) ( a / l) a / l 0 0 l a / otherwise F16 Yao NYU 4
43 Computer Assignment 1. Learn how do fanbeam. ifanbeam work; summarize their functionalities. Type demos on the command line, then select toolbox - > image processing -> transform -> reconstructing an image from projection data. Alternatively, you can use help for each particular function.. Write a MATLAB program that 1) generate a phantom image (you can use a standard phantom provided by MATLAB or construct your own), ) produce equiangular fan beam projections; 3) reconstruct the phantom using filtered backprojection algorithm; Your program should allow the user to specify the number of fan beams, and the number of projections per fan beam, the angular spacing between the projections. Run your program with different number of projections for the same view angle, and with different view angles, and compare the quality. Use the same filter and interpolation algorithm for all the comparisons. Compare the reconstructed image quality obtained with different number of view angles and number of projections per view angle. Also, compare the image quality with those obtained with parallel projections for the same phantom image when the same total number of measurements are used (from your last assignment). You can use the fanbeam() and ifanbeam() functions in MATLAB. F16 Yao NYU 43
Midterm Review. Yao Wang Polytechnic University, Brooklyn, NY 11201
Midterm Review Yao Wang Polytechnic University, Brooklyn, NY 11201 Based on J. L. Prince and J. M. Links, Medical maging Signals and Systems, and lecture notes by Prince. Figures are from the textbook.
More informationModern CT system generations Measurement of attenuation
CT reconstruction repetition & hints Reconstruction in CT and hints to the assignments Jørgen Arendt Jensen October 4, 16 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering
More informationReconstruction in CT and hints to the assignments
Reconstruction in CT and hints to the assignments Jørgen Arendt Jensen October 24, 2016 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound Imaging
More informationMedical Image Reconstruction Term II 2012 Topic 6: Tomography
Medical Image Reconstruction Term II 2012 Topic 6: Tomography Professor Yasser Mostafa Kadah Tomography The Greek word tomos means a section, a slice, or a cut. Tomography is the process of imaging a cross
More informationReconstruction in CT and relation to other imaging modalities
Reconstruction in CT and relation to other imaging modalities Jørgen Arendt Jensen November 16, 2015 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound
More informationCentral Slice Theorem
Central Slice Theorem Incident X-rays y f(x,y) R x r x Detected p(, x ) The thick line is described by xcos +ysin =R Properties of Fourier Transform F [ f ( x a)] F [ f ( x)] e j 2 a Spatial Domain Spatial
More informationTopics. Projections. Review Filtered Backprojection Fan Beam Spiral CT Applications. Bioengineering 280A Principles of Biomedical Imaging
Bioengineering 28A Principles of Biomedical Imaging Fall Quarter 24 X-Rays/CT Lecture 2 Topics Review Filtered Backprojection Fan Beam Spiral CT Applications s I θ (r) = I exp µ(x, y)ds Lr,θ = I exp µ(rcosθ
More informationRadon Transform and Filtered Backprojection
Radon Transform and Filtered Backprojection Jørgen Arendt Jensen October 13, 2016 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound Imaging Department
More informationTopics. Review Filtered Backprojection Fan Beam Spiral CT Applications. Bioengineering 280A Principles of Biomedical Imaging
Bioengineering 28A Principles of Biomedical Imaging Fall Quarter 24 X-Rays/CT Lecture 2 Topics Review Filtered Backprojection Fan Beam Spiral CT Applications Projections I θ (r) = I exp Lr,θ µ(x,y)ds =
More informationGE s Revolution CT MATLAB III: CT. Kathleen Chen March 20, 2018
GE s Revolution CT MATLAB III: CT Kathleen Chen chens18@rpi.edu March 20, 2018 https://www.zmescience.com/medicine/inside-human-body-real-time-gifs-demo-power-ct-scan/ Reminders Make sure you have MATLAB
More informationMulti-slice CT Image Reconstruction Jiang Hsieh, Ph.D.
Multi-slice CT Image Reconstruction Jiang Hsieh, Ph.D. Applied Science Laboratory, GE Healthcare Technologies 1 Image Generation Reconstruction of images from projections. textbook reconstruction advanced
More informationImage Reconstruction from Projection
Image Reconstruction from Projection Reconstruct an image from a series of projections X-ray computed tomography (CT) Computed tomography is a medical imaging method employing tomography where digital
More informationReconstruction in CT and relation to other imaging modalities
Reconstruction in CT and relation to other imaging modalities Jørgen Arendt Jensen November 1, 2017 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound
More information2D Fan Beam Reconstruction 3D Cone Beam Reconstruction. Mario Koerner
2D Fan Beam Reconstruction 3D Cone Beam Reconstruction Mario Koerner Moscow-Bavarian Joint Advanced Student School 2006 March 19 2006 to March 29 2006 Overview 2D Fan Beam Reconstruction Shortscan Reconstruction
More informationELEG404/604: Digital Imaging & Photography
: Digital Imaging & Photography Gonzalo R. Arce Department of Electrical and Computer Engineering University of Delaware Chapter VIII History X-Ray discovery In 1895 Wilhelm Rontgen discovered the X-rays,
More informationIndex. aliasing artifacts and noise in CT images, 200 measurement of projection data, nondiffracting
Index Algebraic equations solution by Kaczmarz method, 278 Algebraic reconstruction techniques, 283-84 sequential, 289, 293 simultaneous, 285-92 Algebraic techniques reconstruction algorithms, 275-96 Algorithms
More informationX-ray tomography. X-ray tomography. Applications in Science. X-Rays. Notes. Notes. Notes. Notes
X-ray tomography Important application of the Fast Fourier transform: X-ray tomography. Also referred to as CAT scan (Computerized Axial Tomography) X-ray tomography This has revolutionized medical diagnosis.
More informationDigital Image Processing
Digital Image Processing Image Restoration and Reconstruction (Image Reconstruction from Projections) Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science and Engineering
More informationBiomedical Image Analysis. Spatial Filtering
Biomedical Image Analysis Contents: Spatial Filtering The mechanics of Spatial Filtering Smoothing and sharpening filters BMIA 15 V. Roth & P. Cattin 1 The Mechanics of Spatial Filtering Spatial filter:
More informationImage Restoration. Yao Wang Polytechnic Institute of NYU, Brooklyn, NY 11201
Image Restoration Yao Wang Polytechnic Institute of NYU, Brooklyn, NY 11201 Partly based on A. K. Jain, Fundamentals of Digital Image Processing, and Gonzalez/Woods, Digital Image Processing Figures from
More information2D Fan Beam Reconstruction 3D Cone Beam Reconstruction
2D Fan Beam Reconstruction 3D Cone Beam Reconstruction Mario Koerner March 17, 2006 1 2D Fan Beam Reconstruction Two-dimensional objects can be reconstructed from projections that were acquired using parallel
More informationFeldkamp-type image reconstruction from equiangular data
Journal of X-Ray Science and Technology 9 (2001) 113 120 113 IOS Press Feldkamp-type image reconstruction from equiangular data Ben Wang a, Hong Liu b, Shiying Zhao c and Ge Wang d a Department of Elec.
More informationEECS490: Digital Image Processing. Lecture #16
Lecture #16 Wiener Filters Constrained Least Squares Filter Computed Tomography Basics Reconstruction and the Radon Transform Fourier Slice Theorem Filtered Backprojections Fan Beams Motion Blurring Model
More informationEnhancement Image Quality of CT Using Single Slice Spiral Technique
Enhancement Image Quality of CT Using Single Slice Spiral Technique Doaa. N. Al Sheack 1 and Dr.Mohammed H. Ali Al Hayani 2 1 2 Electronic and Communications Engineering Department College of Engineering,
More informationTomography. Forward projectionsp θ (r) are known as a Radon transform. Objective: reverse this process to form the original image
C. A. Bouman: Digital Image Processing - January 9, 217 1 Tomography Many medical imaging systems can only measure projections through an object with density f(x,y). Projections must be collected at every
More informationImage Reconstruction 3 Fully 3D Reconstruction
Image Reconstruction 3 Fully 3D Reconstruction Thomas Bortfeld Massachusetts General Hospital, Radiation Oncology, HMS HST.S14, February 25, 2013 Thomas Bortfeld (MGH, HMS, Rad. Onc.) Image Reconstruction
More informationMaterial for Chapter 6: Basic Principles of Tomography M I A Integral Equations in Visual Computing Material
Material for Chapter : Integral Equations in Visual Computing Material Basic Principles of Tomography c 00 Bernhard Burgeth 0 Source: Images Figure : Radon Transform: ttenuation http://en.wikimedia.org/wiki/image:radon_transform.png
More informationA Novel Two-step Method for CT Reconstruction
A Novel Two-step Method for CT Reconstruction Michael Felsberg Computer Vision Laboratory, Dept. EE, Linköping University, Sweden mfe@isy.liu.se Abstract. In this paper we address the parallel beam 2D
More informationBiomedical Imaging. Computed Tomography. Patrícia Figueiredo IST
Biomedical Imaging Computed Tomography Patrícia Figueiredo IST 2013-2014 Overview Basic principles X ray attenuation projection Slice selection and line projections Projection reconstruction Instrumentation
More informationThe Non-uniform Fast Fourier Transform. in Computed Tomography (No ) Supervisor: Prof. Mike Davies
The Non-uniform Fast Fourier Transform in Computed Tomography (No. 2.4.4) Masters Report Supervisor: Prof. Mike Davies Billy Junqi Tang s1408760 AUGEST 12, 2015 UNIVERSITY F EDINBURGH Abstract This project
More informationA practical software platform of integrated reconstruction algorithms and forward algorithms for 2D industrial CT
Journal of X-Ray Science and Technology 13 (2005) 9 21 9 IOS Press A practical software platform of integrated reconstruction algorithms and forward algorithms for 2D industrial CT Hui Li a, Zhaotian Zhang
More informationOptimal threshold selection for tomogram segmentation by reprojection of the reconstructed image
Optimal threshold selection for tomogram segmentation by reprojection of the reconstructed image K.J. Batenburg 1 and J. Sijbers 1 University of Antwerp, Vision Lab, Universiteitsplein 1, B-2610 Wilrijk,
More informationPhys. 428, Lecture 6
Phys. 428, Lecture 6 Mid-Term & Weekly questions Since we are behind in the material, there will be no midterm Instead the final project will be subdivided into graded stages Class Project Pick: An imaging
More informationChapter6 Image Reconstruction
Chapter6 Image Reconstruction Preview 61I 6.1 Introduction 6.2 Reconstruction by Fourier Inversion 6.3 Reconstruction by convolution and backprojection 6.4 Finite series-expansion 1 Preview Reconstruction
More informationCoE4TN4 Image Processing. Chapter 5 Image Restoration and Reconstruction
CoE4TN4 Image Processing Chapter 5 Image Restoration and Reconstruction Image Restoration Similar to image enhancement, the ultimate goal of restoration techniques is to improve an image Restoration: a
More informationReconstruction methods for sparse-data tomography
Reconstruction methods for sparse-data tomography Part B: filtered back-projection Samuli Siltanen Department of Mathematics and Statistics University of Helsinki, Finland samuli.siltanen@helsinki.fi www.siltanen-research.net
More informationTOMOGRAPHIC reconstruction problems are found in
4750 IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 23, NO. 11, NOVEMBER 2014 Improving Filtered Backprojection Reconstruction by Data-Dependent Filtering Abstract Filtered backprojection, one of the most
More informationNon-Stationary CT Image Noise Spectrum Analysis
Non-Stationary CT Image Noise Spectrum Analysis Michael Balda, Björn J. Heismann,, Joachim Hornegger Pattern Recognition Lab, Friedrich-Alexander-Universität Erlangen Siemens Healthcare, Erlangen michael.balda@informatik.uni-erlangen.de
More informationEmission Computed Tomography Notes
Noll (24) ECT Notes: Page 1 Emission Computed Tomography Notes Introduction Emission computed tomography (ECT) is the CT applied to nuclear medicine. There are two varieties of ECT: 1. SPECT single-photon
More informationIntroduction to Medical Imaging. Lecture 6: X-Ray Computed Tomography. CT number (in HU) = Overview. Klaus Mueller
Overview Introduction to Medical Imaging Lecture 6: X-Ray Computed Tomography Scanning: rotate source-detector pair around the patient Klaus Mueller data Computer Science Department Stony Brook University
More informationTopics. View Aliasing. CT Sampling Requirements. Sampling Requirements in CT Sampling Theory Aliasing
Topics Bioengineering 280A Principles of Biomedical Imaging Sampling Requirements in CT Sampling Theory Aliasing Fall Quarter 2010 CT/Fourier Lecture 5 CT Sampling Requirements View Aliasing What should
More informationComputer Vision and Graphics (ee2031) Digital Image Processing I
Computer Vision and Graphics (ee203) Digital Image Processing I Dr John Collomosse J.Collomosse@surrey.ac.uk Centre for Vision, Speech and Signal Processing University of Surrey Learning Outcomes After
More informationAPPLICATION OF RADON TRANSFORM IN CT IMAGE MATCHING Yufang Cai, Kuan Shen, Jue Wang ICT Research Center of Chongqing University, Chongqing, P.R.
APPLICATION OF RADON TRANSFORM IN CT IMAGE MATCHING Yufang Cai, Kuan Shen, Jue Wang ICT Research Center of Chongqing University, Chongqing, P.R.China Abstract: When Industrial Computerized Tomography (CT)
More information9. Computed Tomography (CT) 9.1 Absorption of X-ray
Absorption of X-ray 9. 9.1 Absorption of X-ray X-ray radiography Object I κ Incident Transmitted d I(y) = I e κ(y)d κ : attenuation coefficient In the case of X-ray, it depends on the atomic number. (Heavy
More informationTHE FAN-BEAM scan for rapid data acquisition has
190 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 26, NO. 2, FEBRUARY 2007 Hilbert Transform Based FBP Algorithm for Fan-Beam CT Full Partial Scans Jiangsheng You*, Member, IEEE, Gengsheng L. Zeng, Senior
More informationComputer Vision I. Announcements. Fourier Tansform. Efficient Implementation. Edge and Corner Detection. CSE252A Lecture 13.
Announcements Edge and Corner Detection HW3 assigned CSE252A Lecture 13 Efficient Implementation Both, the Box filter and the Gaussian filter are separable: First convolve each row of input image I with
More informationDEVELOPMENT OF CONE BEAM TOMOGRAPHIC RECONSTRUCTION SOFTWARE MODULE
Rajesh et al. : Proceedings of the National Seminar & Exhibition on Non-Destructive Evaluation DEVELOPMENT OF CONE BEAM TOMOGRAPHIC RECONSTRUCTION SOFTWARE MODULE Rajesh V Acharya, Umesh Kumar, Gursharan
More informationconvolution shift invariant linear system Fourier Transform Aliasing and sampling scale representation edge detection corner detection
COS 429: COMPUTER VISON Linear Filters and Edge Detection convolution shift invariant linear system Fourier Transform Aliasing and sampling scale representation edge detection corner detection Reading:
More informationTomographic Reconstruction
Tomographic Reconstruction 3D Image Processing Torsten Möller Reading Gonzales + Woods, Chapter 5.11 2 Overview Physics History Reconstruction basic idea Radon transform Fourier-Slice theorem (Parallel-beam)
More informationJoint ICTP-TWAS Workshop on Portable X-ray Analytical Instruments for Cultural Heritage. 29 April - 3 May, 2013
2455-5 Joint ICTP-TWAS Workshop on Portable X-ray Analytical Instruments for Cultural Heritage 29 April - 3 May, 2013 Lecture NoteBasic principles of X-ray Computed Tomography Diego Dreossi Elettra, Trieste
More informationAliasing. Can t draw smooth lines on discrete raster device get staircased lines ( jaggies ):
(Anti)Aliasing and Image Manipulation for (y = 0; y < Size; y++) { for (x = 0; x < Size; x++) { Image[x][y] = 7 + 8 * sin((sqr(x Size) + SQR(y Size)) / 3.0); } } // Size = Size / ; Aliasing Can t draw
More informationCover Page. The handle holds various files of this Leiden University dissertation
Cover Page The handle http://hdl.handle.net/1887/8289 holds various files of this Leiden University dissertation Author: Plantagie, L. Title: Algebraic filters for filtered backprojection Issue Date: 2017-0-13
More informationMEDICAL IMAGING 2nd Part Computed Tomography
MEDICAL IMAGING 2nd Part Computed Tomography Introduction 2 In the last 30 years X-ray Computed Tomography development produced a great change in the role of diagnostic imaging in medicine. In convetional
More informationAlgebraic Iterative Methods for Computed Tomography
Algebraic Iterative Methods for Computed Tomography Per Christian Hansen DTU Compute Department of Applied Mathematics and Computer Science Technical University of Denmark Per Christian Hansen Algebraic
More informationSingle Particle Reconstruction Techniques
T H E U N I V E R S I T Y of T E X A S S C H O O L O F H E A L T H I N F O R M A T I O N S C I E N C E S A T H O U S T O N Single Particle Reconstruction Techniques For students of HI 6001-125 Computational
More informationImage Restoration Chapter 5. Prof. Vidya Manian Dept. of Electrical and Computer Engineering INEL 5327 ECE, UPRM
Image Processing Image Restoration Chapter 5 Prof. Vidya Manian Dept. of Electrical and Computer Engineering g 1 Overview A model of the Image Degradation/Restoration Process Noise Models Restoration in
More informationLecture 5: Frequency Domain Transformations
#1 Lecture 5: Frequency Domain Transformations Saad J Bedros sbedros@umn.edu From Last Lecture Spatial Domain Transformation Point Processing for Enhancement Area/Mask Processing Transformations Image
More informationRevisit of the Ramp Filter
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 62, NO. 1, FEBRUARY 2015 131 Revisit of the Ramp Filter Gengsheng L. Zeng, Fellow, IEEE Abstract An important part of the filtered backprojection (FBP) algorithm
More informationTomographic Image Reconstruction in Noisy and Limited Data Settings.
Tomographic Image Reconstruction in Noisy and Limited Data Settings. Syed Tabish Abbas International Institute of Information Technology, Hyderabad syed.abbas@research.iiit.ac.in July 1, 2016 Tabish (IIIT-H)
More informationDigital Image Processing
Digital Image Processing Third Edition Rafael C. Gonzalez University of Tennessee Richard E. Woods MedData Interactive PEARSON Prentice Hall Pearson Education International Contents Preface xv Acknowledgments
More informationPrinciples of Computerized Tomographic Imaging
Principles of Computerized Tomographic Imaging Parallel CT, Fanbeam CT, Helical CT and Multislice CT Marjolein van der Glas August 29, 2000 Abstract The total attenuation suffered by one beam of x-rays
More informationAn Introduc+on to Mathema+cal Image Processing IAS, Park City Mathema2cs Ins2tute, Utah Undergraduate Summer School 2010
An Introduc+on to Mathema+cal Image Processing IAS, Park City Mathema2cs Ins2tute, Utah Undergraduate Summer School 2010 Luminita Vese Todd WiCman Department of Mathema2cs, UCLA lvese@math.ucla.edu wicman@math.ucla.edu
More informationPhase problem and the Radon transform
Phase problem and the Radon transform Andrei V. Bronnikov Bronnikov Algorithms The Netherlands The Radon transform and applications Inverse problem of phase-contrast CT Fundamental theorem Image reconstruction
More informationSolutions For Homework #7
Solutions For Homework #7 Problem :[ pts] Let f(r) = r = x2 + y 2 () We compute the Hankel Transform of f(r) by first computing its Abel Transform and then calculating the D Fourier Transform of the result.
More informationEE123 Digital Signal Processing
Multi-Dimensional Signals EE23 Digital Signal Processing Lecture Our world is more complex than D Images: f(x,y) Videos: f(x,y,t) Dynamic 3F scenes: f(x,y,z,t) Medical Imaging 3D Video Computer Graphics
More informationApproximating Algebraic Tomography Methods by Filtered Backprojection: A Local Filter Approach
Fundamenta Informaticae 135 (2014) 1 19 1 DOI 10.3233/FI-2014-1109 IOS Press Approximating Algebraic Tomography Methods by Filtered Backprojection: A Local Filter Approach Linda Plantagie Centrum Wiskunde
More informationsurface Image reconstruction: 2D Fourier Transform
2/1/217 Chapter 2-3 K-space Intro to k-space sampling (chap 3) Frequenc encoding and Discrete sampling (chap 2) Point Spread Function K-space properties K-space sampling principles (chap 3) Basic Contrast
More informationCover Page. The handle holds various files of this Leiden University dissertation
Cover Page The handle http://hdl.handle.net/887/48289 holds various files of this Leiden University dissertation Author: Plantagie, L. Title: Algebraic filters for filtered backprojection Issue Date: 207-04-3
More informationBME I5000: Biomedical Imaging
1 Lucas Parra, CCNY BME I5000: Biomedical Imaging Lecture 4 Computed Tomography Lucas C. Parra, parra@ccny.cuny.edu some slides inspired by lecture notes of Andreas H. Hilscher at Columbia University.
More informationCone-beam and fan-beam image reconstruction algorithms based on spherical and circular harmonics
INSTITUTE OF PHYSICS PUBLISHING Phys. Med. Biol. 49 (24) 2239 2256 PHYSICS IN MEDICINE AND BIOLOGY PII: S31-9155(4)7992-6 Cone-beam and fan-beam image reconstruction algorithms based on spherical and circular
More informationIntroduction to Medical Imaging. Cone-Beam CT. Klaus Mueller. Computer Science Department Stony Brook University
Introduction to Medical Imaging Cone-Beam CT Klaus Mueller Computer Science Department Stony Brook University Introduction Available cone-beam reconstruction methods: exact approximate algebraic Our discussion:
More informationNIH Public Access Author Manuscript Med Phys. Author manuscript; available in PMC 2009 March 13.
NIH Public Access Author Manuscript Published in final edited form as: Med Phys. 2008 February ; 35(2): 660 663. Prior image constrained compressed sensing (PICCS): A method to accurately reconstruct dynamic
More informationCompressed Sensing for Electron Tomography
University of Maryland, College Park Department of Mathematics February 10, 2015 1/33 Outline I Introduction 1 Introduction 2 3 4 2/33 1 Introduction 2 3 4 3/33 Tomography Introduction Tomography - Producing
More informationBME I5000: Biomedical Imaging
BME I5000: Biomedical Imaging Lecture 1 Introduction Lucas C. Parra, parra@ccny.cuny.edu 1 Content Topics: Physics of medial imaging modalities (blue) Digital Image Processing (black) Schedule: 1. Introduction,
More informationAnalytical Reconstructions of Regions of Interest in Medical Imaging
Analytical s of Regions of Interest in Medical Imaging Kévin Polisano Supervised by : Laurent Desbat 09/18/2012 Kévin Polisano Analytical s of Regions of Interest in Medical Imaging 09/18/2012 1 / 45 1
More informationA study on the section sensitivity profile in multi-row-detector spiral CT
Journal of X-Ray Science and Technology (3) IOS Press A study on the section sensitivity profile in multi-row-detector spiral CT Ge Wang a, Mark Madsen a, Katherine Redford a, Shiying Zhao b and Michael
More informationMoscow-Bavarian Joint Advanced Student School 2006 / Medical Imaging Principles of Computerized Tomographic Imaging and Cone-Beam Reconstruction
Line Integrals Line integrals represent the integral of some parameter of the object along the line (e.g. attenuation of x-rays) Object: f(x,y) Line: x cosθ + y sinθ = t Line integral / Radon transform:
More informationBME I5000: Biomedical Imaging
1 Lucas Parra, CCNY BME I5000: Biomedical Imaging Lecture 11 Point Spread Function, Inverse Filtering, Wiener Filtering, Sharpening,... Lucas C. Parra, parra@ccny.cuny.edu Blackboard: http://cityonline.ccny.cuny.edu/
More informationGRAPHICAL USER INTERFACE (GUI) TO STUDY DIFFERENT RECONSTRUCTION ALGORITHMS IN COMPUTED TOMOGRAPHY
GRAPHICAL USER INTERFACE (GUI) TO STUDY DIFFERENT RECONSTRUCTION ALGORITHMS IN COMPUTED TOMOGRAPHY A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Engineering
More informationContinuous and Discrete Image Reconstruction
25 th SSIP Summer School on Image Processing 17 July 2017, Novi Sad, Serbia Continuous and Discrete Image Reconstruction Péter Balázs Department of Image Processing and Computer Graphics University of
More information2-D Reconstruction Hannes Hofmann. 2-D Reconstruction. MB-JASS 2006, March 2006
2-D Reconstruction MB-JASS 2006, 19 29 March 2006 Computer Tomography systems use X-rays to acquire images from the inside of the human body. Out of the projection images the original object is reconstructed.
More informationImage Acquisition Systems
Image Acquisition Systems Goals and Terminology Conventional Radiography Axial Tomography Computer Axial Tomography (CAT) Magnetic Resonance Imaging (MRI) PET, SPECT Ultrasound Microscopy Imaging ITCS
More informationAdaptive algebraic reconstruction technique
Adaptive algebraic reconstruction technique Wenkai Lua) Department of Automation, Key State Lab of Intelligent Technology and System, Tsinghua University, Beijing 10084, People s Republic of China Fang-Fang
More informationDigital Image Processing. Lecture 6
Digital Image Processing Lecture 6 (Enhancement in the Frequency domain) Bu-Ali Sina University Computer Engineering Dep. Fall 2016 Image Enhancement In The Frequency Domain Outline Jean Baptiste Joseph
More informationImplementing Fast Hierarchical Back Projection
Jason ECE 558 Final Project Paper May 9, 2007 Abstract Filtered back projection (FBP) is a commonly used image reconstruction scheme for parallel beam tomography. The algorithm implemented in this project,
More informationOutline. Foundations of Computer Graphics (Spring 2012)
Foundations of Computer Graphics (Spring 2012) CS 184, Lectures 19: Sampling and Reconstruction http://inst.eecs.berkeley.edu/~cs184 Basic ideas of sampling, reconstruction, aliasing Signal processing
More informationFourier Transform in Image Processing. CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012)
Fourier Transform in Image Processing CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012) 1D: Common Transform Pairs Summary source FT Properties: Convolution See book DIP 4.2.5:
More informationEvaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT
Evaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT Qiao Yang 1,4, Meng Wu 2, Andreas Maier 1,3,4, Joachim Hornegger 1,3,4, Rebecca Fahrig
More informationImage Alignment and Application of 2D Fast Rotational Matching in Single Particle cryo-electron Microscopy. Yao Cong. Baylor College of Medicine
Image Alignment and Application of 2D Fast Rotational Matching in Single Particle cryo-electron Microscopy Yao Cong Baylor College of Medicine cryo-em & single particle analysis Single-particle electron
More informationConsistency in Tomographic Reconstruction by Iterative Methods
Consistency in Tomographic Reconstruction by Iterative Methods by M. Reha Civanlar and H.J. Trussell Center for Communications and Signal Processing Department of Electrical and Computer Engineering North
More informationBiophysical Techniques (BPHS 4090/PHYS 5800)
Biophysical Techniques (BPHS 4090/PHYS 5800) Instructors: Prof. Christopher Bergevin (cberge@yorku.ca) Schedule: MWF 1:30-2:30 (CB 122) Website: http://www.yorku.ca/cberge/4090w2017.html York University
More informationBiomedical Image Analysis. Point, Edge and Line Detection
Biomedical Image Analysis Point, Edge and Line Detection Contents: Point and line detection Advanced edge detection: Canny Local/regional edge processing Global processing: Hough transform BMIA 15 V. Roth
More informationComputational Aspects of MRI
David Atkinson Philip Batchelor David Larkman Programme 09:30 11:00 Fourier, sampling, gridding, interpolation. Matrices and Linear Algebra 11:30 13:00 MRI Lunch (not provided) 14:00 15:30 SVD, eigenvalues.
More informationTomographic Algorithm for Industrial Plasmas
Tomographic Algorithm for Industrial Plasmas More info about this article: http://www.ndt.net/?id=22342 1 Sudhir K. Chaudhary, 1 Kavita Rathore, 2 Sudeep Bhattacharjee, 1 Prabhat Munshi 1 Nuclear Engineering
More informationAn approximate cone beam reconstruction algorithm for gantry-tilted CT
An approximate cone beam reconstruction algorithm for gantry-tilted CT Ming Yan a, Cishen Zhang ab, Hongzhu Liang a a School of Electrical & Electronic Engineering, Nanyang Technological University, Singapore;
More informationFourier analysis and sampling theory
Reading Required: Shirley, Ch. 9 Recommended: Fourier analysis and sampling theory Ron Bracewell, The Fourier Transform and Its Applications, McGraw-Hill. Don P. Mitchell and Arun N. Netravali, Reconstruction
More informationIntroduction to Biomedical Imaging
Alejandro Frangi, PhD Computational Imaging Lab Department of Information & Communication Technology Pompeu Fabra University www.cilab.upf.edu X-ray Projection Imaging Computed Tomography Digital X-ray
More informationMedical Imaging BMEN Spring 2016
Name Medical Imaging BMEN 420-501 Spring 2016 Homework #4 and Nuclear Medicine Notes All questions are from the introductory Powerpoint (based on Chapter 7) and text Medical Imaging Signals and Systems,
More informationImage Sampling and Quantisation
Image Sampling and Quantisation Introduction to Signal and Image Processing Prof. Dr. Philippe Cattin MIAC, University of Basel 1 of 46 22.02.2016 09:17 Contents Contents 1 Motivation 2 Sampling Introduction
More informationTESTING OF THE CIRCLE AND LINE ALGORITHM IN THE SETTING OF MICRO-CT
SCA2016-080 1/7 TESTING OF THE CIRCLE AND LINE ALGORITHM IN THE SETTING OF MICRO-CT Alexander Katsevich 1, 2 and Michael Frenkel 1 1 itomography Corp., 2 University of Central Florida (UCF) This paper
More information