9. Computed Tomography (CT) 9.1 Absorption of X-ray
|
|
- Sharlene Mosley
- 5 years ago
- Views:
Transcription
1 Absorption of X-ray 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 atom large κ.) I If κ is depth dependent, κd d κ(x,y)dx κ(x,y)dx (κ(x,y) = Not Object) Only the integral of κ along optical path can be obtained from X-ray radiography. 149/197
2 Projection from several directions 9.2 Projection from several directions I I(ξ,θ = ) I ξ I(ξ,θ = π 2 ) The distribution of κ including depth distribution, which is called tomography, can be obtained from several projected data with different directions. Computed Tomography ξ 15/197
3 Projection from several directions Number of projections and Num. of internal nodes 3x3 unknowns 3 projs. (observed) 3 projs. (observed) Num. of unknowns > Num. of Obs. Cannot solve. Num. of unknowns = Num. of Obs. Cannot distinguish. To obtain more projection, other projections with different directions are needed. θ [,π] 151/197
4 Schematic of forward and backward-projection 9.3 Schematic of forward and backward-projection Forward projection (measuring process) Forward projection Integral along beam path p(ξ,θ) = η κ(x,y)dη θ (x x(ξ,θ), y y(ξ,θ)) Accumulate along path deg 45 deg = 9 deg 135 deg /197
5 Schematic of forward and backward-projection Backward projection (1) Simple backprojection 1 Map averaged value of the projection data along path for each angle. = 2 Take an average of mapped data for each pixel Blurrier than original. 8 Original = deg 45 deg deg 135 deg /197
6 Schematic of forward and backward-projection Backward projection (2) Filtered backward projection In the simple backward projection, the reconstracted result is blurred. To reduce the blurring, edge enhancement filter is applied to the projection data. Edge enhancement filter g n = f n kf n = f n k(f n 1 2f n +f n+1 ) = kf n 1 +(2k +1)f n kf n+1 (k = 1) g n = +1 i= 1 w i f n i (w 1,w,w +1 ) = ( 1,+3, 1) 154/197
7 Schematic of forward and backward-projection Filtered backward projection 1 Edge enhancement of projection data: p (e.g.) p j = +1 w i p j i i= 1 (w 1,w,w +1 ) = ( 1,3, 1) p 1 p Apply simple backward projection using p. 1 = In this example, the reconstructed field is identical to original. deg 45 deg deg 135 deg /197
8 Radon Transform 9.4 Radon Transform Coordinate system η I (ξ,θ) y r θ ξ x I(ξ, θ) ξ Inner product of e x : xe x e x +ye x e y = ξe x e ξ +ηe x e η η y x = ξcosθ ηsinθ θ ξ x Expression of point r r = xe x +ye y { x = +ξcosθ ηsinθ y = +ξsinθ +ηcosθ { ξ = +xcosθ +ysinθ η = xsinθ+ycosθ = ξe ξ +ηe η 156/197
9 Radon Transform Radon Transform Projected data (known) I(ξ,θ) = I e L κ(r)dl L {r(ξ,η;θ) ξ = ξ (const.)} Sinogram (known) p(ξ,θ) log I(ξ,θ) I Radon Transform Integral over straight line in 2-D space. p(ξ,θ) = κ(r(ξ, θ)) dl L Extend from line integral to 2-D area integral ξ = xcosθ+ysinθ L [ ]dl = [ ]dη ξ=ξ = [ ]δ(ξ ξ )dξ dη η y dl(ξ = ξ ) = [ ]δ(ξ ξ (x,y))dxdy = [ ]δ(ξ (xcosθ+ysinθ))dxdy x ξ ξ 157/197
10 Projection slice theorem 9.5 Projection slice theorem Projection (Sinogram) p(ξ,θ) = κ(x,y)δ(ξ (xcosθ+ysinθ))dxdy (1) FT with respect to ξ P(k ξ,θ) = κ(x,y)δ(ξ (xcosθ+ysinθ))e jkξξ dxdydξ = κ(x,y)e jkξ(xcosθ+ysinθ) dxdy (2) 158/197
11 Projection slice theorem 2-D Fourier Transform in polar coordinate system Forward transform F(k x,k y ) = f(x,y)e j(kxx+kyy) dxdy (k x = kcosθ, k y = ksinθ) = f(x,y)e jk(xcosθ+ysinθ) dxdy F (k,θ) (3) Inverse transform f(x,y) = 1 4π 2 F(k x,k y )e +j(kxx+kyy) dk x dk y = 1 4π 2 2π ( dk x dk y = 2π kdθdk) F (k,θ)e +jk(xcosθ+ysinθ) k dθdk (4) 159/197
12 Projection slice theorem Eq. (2) and Eq. (3) are same. Projection Slice Theorem P(k ξ,θ)isexpressed byfourier transformofκ(x,y) inpolar coordinate system. Since P(k ξ,θ) is known, κ(x,y) is obtained by the inverse FT using (4). κ(x,y) = 1 4π 2 2π P(k ξ,θ)e +jk ξ(xcosθ+ysinθ) k ξ dθdk ξ (5) 16/197
13 Projection slice theorem Projection from opposite direction κ(x,y) = 1 4π 2 = 1 4π 2 2π π P(k ξ,θ)e +jk ξ(xcosθ+ysinθ) k ξ dθdk ξ P(k ξ,θ)e +jk ξ(xcosθ+ysinθ) k ξ dθdk ξ 2π π 2π k ξ dθdk ξ = k ξ dθdk ξ + k ξ dθdk ξ π Projection from opposite direction ( ) p(ξ,θ±π) = p( ξ,θ) θ θ±π P(k ξ ξ ξ,θ ±π) = P( k ξ,θ) e +jkξ(xcos(θ±π)+ysin(θ±π)) = e jk ξ(xcosθ+ysinθ) 2π 2nd Term = P(k ξ,θ)e +jkξ(xcosθ+ysinθ) ( k ξ dθdk ξ θ = θ π ) π π = P( k ξ,θ )e jkξ(xcosθ +ysinθ ) k ξ dθ ( ) dk ξ k ξ = k ξ π = P(k ξ,θ )e +jk ξ (xcosθ +ysinθ ) ( k ξ )dθ ( dk ξ ) π = P(k ξ,θ )e +jk ξ (xcosθ +ysinθ ) k ξ dθ dk ξ 161/197
14 Reconstruction by using Fourier transform 9.6 Reconstruction by using Fourier transform κ(x,y) = 1 4π 2 = 1 2π π π q(ξ(x,y),θ) = ξ {}}{ P(k ξ,θ)e +jk ξ(xcosθ+ysinθ) kξ dθdk ξ [P(k ξ,θ) k ξ ]e +jkξξ dk ξ dθ = π 1 2π } {{ } =F 1 k {P((k ξ,θ) k ξ } q(ξ,θ) ξ [P(k ξ,θ) k ξ ]e +jk ξ(xcosθ+ysinθ) dk ξ π q(ξ(x,y),θ)dθ }{{} Average with θ. κ(x,y) = 1 2 } F 1 k ξ {F ξ {p(ξ,θ)} ξ H(k ξ ) k ξ θ (H(k ξ ) = k ξ in the case of Ramp function) 162/197
15 Reconstruction by using Fourier transform 1 Sinogram : p(ξ, θ) 2 FWD FT with ξ : P(k ξ,θ) = 3 Filtering (Weight with k ξ ) : k ξ P(k ξ,θ) p(ξ,θ)e jk ξξ dξ 4 INV FT with k ξ : q(ξ,θ) = 1 P(k ξ,θ) k ξ e +jkξξ dk ξ 2π 5 Backward projection (Coordinate transform and integrate with θ) : κ(x,y) = 1 π q(xcosθ+ysinθ,θ)dθ 2π Two FT (FWD and INV) are needed for a certain θ. 163/197
16 Reconstruction by using Fourier transform Filtered Back-Projection Before IFT, k ξ is multiplied. Since this factor k ξ is considered as a filter in the spectral domain, the method based on FT is called Filtered Back-projection (FBP). In the actual computation, k ξ [, ] [ k max,+k max ]. k max is Nyquist frequency determined by sampling interval. To avoid ringing artifact caused by high frequency component, another filter can be applied. (e.g. Shepp-Logan filter) H(k ξ ) = 2k max sin πk ξ π k max Ramp filter Shepp-Logan filter 2k max k max +k max 164/197
17 Reconstruction by using Fourier transform Handling of discrete data p(ξ n,θ m ) = p(n ξ,m θ) (n,m) : Integer e.g. (Interpolation for ξ) q(xcosθ+ysinθ,θ) q(ξ n+1,θ) q(ξ n,θ) κ(x i,y j ) = κ(i x,j y) (i,j) : Integer ξ ξ n ξ n+1 xcosθ+ysinθ In order to evaluate q(xcosθ+ysinθ,θ) from q(ξ,θ) interpolation are needed. 165/197
18 Reconstruction by using Fourier transform Reconstruction using convolution { c(x) = a(x )b(x x )dx C(k) = A(k)B(k) q(ξ,θ) = 1 2π kmax } c(x) = 1 A(k)B(k)e +jkx dk 2π k max P(k ξ,θ)h(k ξ )e +jk ξξ dk ξ = ξmax ξ max p(ξ,θ)h(ξ ξ )dξ h(ξ) If H(k ξ ) = k ξ, h(ξ) = 1 2π = kmax k ξ e +jkξξ dk ξ k max 1 2π k2 max (ξ = ) 1 k max sin(k max ξ)+ 1 π ξ ξ 2 (cos(k maxξ) 1) (ξ ) 5 5 ξ/ ξ Subtract by neighbors Edge Enhancement 166/197
19 Reconstruction by using Fourier transform 1 Sinogram : p(ξ, θ) 2 Convolution : q(ξ,θ) = p(ξ,θ)h(ξ ξ )dξ 3 Back-projection : κ(x,y) = 1 π q(xcosθ +ysinθ,θ)dθ 2π Only one convolution for each θ. No FT. κ(x, y) p (x,y;) q (x,y;) ( p,q : contrast is enhanced to display (erf(i/σ)) ) Fourier transform Convolutinal integral Number of multiplifications for each θ DFT 2 2N 2 FFT 2 2N logn All points N 2 Neighboring M pts. (M N) MN Faster computation than DFT, if the convolution is applied to neighboring points only. 167/197
20 Reconstruction by using Fourier transform Filtered Back-projection and Simple BP Filtered Back-Projection κ(x,y) = 1 } F 1 2 k ξ {F ξ {p(ξ,θ)} ξ H(k ξ ) Simple Back-projection H(k) = 1 κ(x,y) = 1 2 = 1 2π F 1 { k ξ π F ξ {p(ξ,θ)} ξ } p(ξ,θ)dθ No need to FT Fast The reconstructed distribution is blurred. k ξ θ k ξ θ = 1 2 p(ξ,θ) θ Iterate two procedures of projection and back-projection. 168/197
21 Iterative Reconstruction 9.7 Iterative Reconstruction Applying the forward projection (FP) to the reconstructed field obtained by the backward projection (BP), we can evaluate the error. The BP of the error is added to the field obtained in the previous step. The simple BP is used for the BP algorithm, since the simple BP is fast. κ(x, y) p(ξ, θ) κ True p Obs κ 2 F B (1) (3) p 1 (5) κ 1 (4) p 1 κ 1 (2) (1) BP for the proj. : κ 1 = B { p Obs} (2) FP for the field : p 1 = F {κ 1 } (3) Under-estimation : p 1 = p Obs p 1 (4) BP for the under-est. : κ 1 = B{ p 1 } (5) Update the field : κ 2 = κ 1 +α κ 1 ( α : Relaxation factor for stable reconstruction < α 1 ) 169/197
22 Iterative Reconstruction Simple Back-projection Sinogram p(ξ n,θ m ) = κ(x,y)dl κ(x,y L nm ) l nm L nm κ(x,y L nm ) = p(ξ ( ) n,θ m ) l nm = dl l nm L nm Fraction of projection is mapped onto the internal distribution. κ(x i,y j ) = 1 p(ξ n,θ m ) a i,j,n,m N n l nm m n a i,j,n,m : Overlap area fraction between the pixel (x i,y j ) and the beam L nm with width ξ. 17/197
23 Example of reconstruction by simulation 9.8 Example of reconstruction by simulation True κ(x,y) Filtered Back-Projection (Filter : Ramp function) κ(x, y) [,2] κ(x, y) [-.2,+.2] Sinogram p(ξ, θ) κ(x,y) [,2], N x = N y = 1, N ξ = 1, N θ = 45( θ = 4deg) κ 2 κ 2 =.15 Line artifact 171/197
24 Example of reconstruction by simulation True κ(x,y) Iterative reconstruction N = 1 N = 2 N = 1 N = 2 Sinogram p(ξ, θ) κ(x,y) [,2], N x = N y = 1, N ξ = 1, N θ = 45( θ = 4deg) κ 2 =.42 κ 2 =.25 κ 2 =.17 Reduction of edge blurring κ 2 = /197
Image 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationTheoretically Perfect Sensor
Sampling 1/67 Sampling The ray tracer samples the geometry, only gathering information from the parts of the world that interact with a finite number of rays In contrast, a scanline renderer can push all
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 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 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 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 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 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 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 informationRefractive Index Map Reconstruction in Optical Deflectometry Using Total-Variation Regularization
Refractive Index Map Reconstruction in Optical Deflectometry Using Total-Variation Regularization L. Jacques 1, A. González 2, E. Foumouo 2 and P. Antoine 2 1 Institute of Information and Communication
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 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 informationTheoretically Perfect Sensor
Sampling 1/60 Sampling The ray tracer samples the geometry, only gathering information from the parts of the world that interact with a finite number of rays In contrast, a scanline renderer can push all
More informationComputed Tomography (Part 2)
EL-GY 6813 / BE-GY 603 / G16.446 Medical Imaging Computed Tomography (Part ) Yao Wang Polytechnic School of Engineering New York University, Brooklyn, NY 1101 Based on Prince and Links, Medical Imaging
More informationA method and algorithm for Tomographic Imaging of highly porous specimen using Low Frequency Acoustic/Ultrasonic signals
More Info at Open Access Database www.ndt.net/?id=15210 A method and algorithm for Tomographic Imaging of highly porous specimen using Low Frequency Acoustic/Ultrasonic signals Subodh P S 1,a, Reghunathan
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 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 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 informationApplication of optimal sampling lattices on CT image reconstruction and segmentation or three dimensional printing
Application of optimal sampling lattices on CT image reconstruction and segmentation or three dimensional printing XIQIANG ZHENG Division of Health and Natural Sciences, Voorhees College, Denmark, SC 29042
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 informationMicro and Nano- Tomography
Micro and Nano- Tomography Wilhelm Conrad Roentgen first Nobel prize in physics (1901) Use of X-rays for imaging Compressing 3-d information into a 2-d image Development of tomography Sir Godfrey N. Hounsfield
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 informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi:10.1038/nature10934 Supplementary Methods Mathematical implementation of the EST method. The EST method begins with padding each projection with zeros (that is, embedding
More informationChapter 11 Image Processing
Chapter Image Processing Low-level Image Processing Operates directly on a stored image to improve or enhance it. Stored image consists of a two-dimensional array of pixels (picture elements): Origin (0,
More informationGeneralized Filtered Backprojection for Digital Breast Tomosynthesis Reconstruction
Generalized Filtered Backprojection for Digital Breast Tomosynthesis Reconstruction Klaus Erhard a, Michael Grass a, Sebastian Hitziger b, Armin Iske b and Tim Nielsen a a Philips Research Europe Hamburg,
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 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 informationCS 179: GPU Programming. Lecture 12 / Homework 4
CS 179: GPU Programming Lecture 12 / Homework 4 Admin Lab 4 is out Due Wednesday, April 27 @3pm Come to OH this week, this set is more difficult than before. Breadth-First Search Given source vertex S:
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 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 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 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 informationMotivation. The main goal of Computer Graphics is to generate 2D images. 2D images are continuous 2D functions (or signals)
Motivation The main goal of Computer Graphics is to generate 2D images 2D images are continuous 2D functions (or signals) monochrome f(x,y) or color r(x,y), g(x,y), b(x,y) These functions are represented
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 informationMidterm 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 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 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 informationReading. 2. Fourier analysis and sampling theory. Required: Watt, Section 14.1 Recommended:
Reading Required: Watt, Section 14.1 Recommended: 2. Fourier analysis and sampling theory Ron Bracewell, The Fourier Transform and Its Applications, McGraw-Hill. Don P. Mitchell and Arun N. Netravali,
More informationSteen Moeller Center for Magnetic Resonance research University of Minnesota
Steen Moeller Center for Magnetic Resonance research University of Minnesota moeller@cmrr.umn.edu Lot of material is from a talk by Douglas C. Noll Department of Biomedical Engineering Functional MRI Laboratory
More informationCS 179: GPU Programming. Lecture 11 / Homework 4
CS 179: GPU Programming Lecture 11 / Homework 4 Breadth-First Search Given source vertex S: Find min. #edges to reach every vertex from S (Assume source is vertex 0) 0 1 1 2 2 3 Sequential pseudocode:
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 informationComputer Graphics. Texture Filtering & Sampling Theory. Hendrik Lensch. Computer Graphics WS07/08 Texturing
Computer Graphics Texture Filtering & Sampling Theory Hendrik Lensch Overview Last time Texture Parameterization Procedural Shading Today Texturing Filtering 2D Texture Mapping Forward mapping Object surface
More informationImage Processing. Application area chosen because it has very good parallelism and interesting output.
Chapter 11 Slide 517 Image Processing Application area chosen because it has very good parallelism and interesting output. Low-level Image Processing Operates directly on stored image to improve/enhance
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 informationThe main goal of Computer Graphics is to generate 2D images 2D images are continuous 2D functions (or signals)
Motivation The main goal of Computer Graphics is to generate 2D images 2D images are continuous 2D functions (or signals) monochrome f(x,y) or color r(x,y), g(x,y), b(x,y) These functions are represented
More informationCompressive Sensing based on Optical Coding
Compressive Sensing based on Optical Coding Marko Stamenkovic Supervisor: Dipl.-Phys. Markus Kamm January 25, 2017 Abstract In this work we analyze the properties and specificities of the dynamical, electrically
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 informationEdge detection. Gradient-based edge operators
Edge detection Gradient-based edge operators Prewitt Sobel Roberts Laplacian zero-crossings Canny edge detector Hough transform for detection of straight lines Circle Hough Transform Digital Image Processing:
More informationAttenuation map reconstruction from TOF PET data
Attenuation map reconstruction from TOF PET data Qingsong Yang, Wenxiang Cong, Ge Wang* Department of Biomedical Engineering, Rensselaer Polytechnic Institute, Troy, NY 80, USA *Ge Wang (ge-wang@ieee.org)
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 information3. Image formation, Fourier analysis and CTF theory. Paula da Fonseca
3. Image formation, Fourier analysis and CTF theory Paula da Fonseca EM course 2017 - Agenda - Overview of: Introduction to Fourier analysis o o o o Sine waves Fourier transform (simple examples of 1D
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 informationXRDUG Seminar III Edward Laitila 3/1/2009
XRDUG Seminar III Edward Laitila 3/1/2009 XRDUG Seminar III Computer Algorithms Used for XRD Data Smoothing, Background Correction, and Generating Peak Files: Some Features of Interest in X-ray Diffraction
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 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 informationMapping Between Digital and Continuous Projections via the Discrete Radon Transform in Fourier Space
Mapping Between Digital and Continuous Projections via the Discrete Radon Transform in Fourier Space Andrew Kingston and Imants Svalbe Center for X-ray Physics and Imaging, School of Physics and Materials
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 informationDigital Signal Processing Lecture Notes 22 November 2010
Digital Signal Processing Lecture otes 22 ovember 2 Topics: Discrete Cosine Transform FFT Linear and Circular Convolution Rate Conversion Includes review of Fourier transforms, properties of Fourier transforms,
More informationTopics. 2D Signal. a b. c d. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2004 Lecture 4 2D Fourier Transforms
Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2004 Lecture 4 2D Fourier Transforms Topics. 2D Signal Representations 2. 2D Fourier Transform 3. Transform Pairs 4. FT Properties 2D Signal
More informationx' = c 1 x + c 2 y + c 3 xy + c 4 y' = c 5 x + c 6 y + c 7 xy + c 8
1. Explain about gray level interpolation. The distortion correction equations yield non integer values for x' and y'. Because the distorted image g is digital, its pixel values are defined only at integer
More informationTemperature Distribution Measurement Based on ML-EM Method Using Enclosed Acoustic CT System
Sensors & Transducers 2013 by IFSA http://www.sensorsportal.com Temperature Distribution Measurement Based on ML-EM Method Using Enclosed Acoustic CT System Shinji Ohyama, Masato Mukouyama Graduate School
More informationEstimation of Missing Fan-Beam Projections using Frequency Consistency Conditions
Estimation of Missing Fan-Beam Projections using Frequency Consistency Conditions Marcel Pohlmann, Martin Berger, Andreas Maier, Joachim Hornegger and Rebecca Fahrig Abstract Reducing radiation dose is
More informationPlanar tomosynthesis reconstruction in a parallel-beam framework via virtual object reconstruction
Planar tomosynthesis reconstruction in a parallel-beam framework via virtual object reconstruction Brian E. Nett a,shuaileng a and Guang-Hong Chen a,b a Department of Medical Physics, University of Wisconsin-Madison,
More informationExpectation Maximization and Total Variation Based Model for Computed Tomography Reconstruction from Undersampled Data
Expectation Maximization and Total Variation Based Model for Computed Tomography Reconstruction from Undersampled Data Ming Yan and Luminita A. Vese Department of Mathematics, University of California,
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 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 informationMichael Moody School of Pharmacy University of London 29/39 Brunswick Square London WC1N 1AX, U.K.
This material is provided for educational use only. The information in these slides including all data, images and related materials are the property of : Michael Moody School of Pharmacy University of
More informationR θ (t) 2. (θ) s. (t) θ 2. f(x,y) f(x,y) s 1. Fourier transform v. object
Parallelisation of Tomographic Reconstruction Methods M.A. Westenberg and J.B.T.M. Roerdink Department of Mathematics and Computing Science, University of Groningen P.O. Bo 800, 9700 AV Groningen, The
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 informationAcknowledgments and financial disclosure
AAPM 2012 Annual Meeting Digital breast tomosynthesis: basic understanding of physics principles James T. Dobbins III, Ph.D., FAAPM Director, Medical Physics Graduate Program Ravin Advanced Imaging Laboratories
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 informationInvestigation on reconstruction methods applied to 3D terahertz computed Tomography
Investigation on reconstruction methods applied to 3D terahertz computed Tomography B. Recur, 3 A. Younus, 1, P. Mounaix 1, S. Salort, 2 B. Chassagne, 2 P. Desbarats, 3 J-P. Caumes, 2 and E. Abraham 1
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 informationMultilevel Optimization for Multi-Modal X-ray Data Analysis
Multilevel Optimization for Multi-Modal X-ray Data Analysis Zichao (Wendy) Di Mathematics & Computer Science Division Argonne National Laboratory May 25, 2016 2 / 35 Outline Multi-Modality Imaging Example:
More informationBeing edited by Prof. Sumana Gupta 1
Being edited by Prof. Sumana Gupta 1 Introduction Digital Image Processing: refers to processing of 2-D picture by a digital Computer- in a broader context refers to digital proc of any 2-d data. A digital
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 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 informationSurface Roughness Scattering in Grazing-incidence X-ray Telescopes. Outline
Surface Roughness Scattering in Grazing-incidence X-ray Telescopes Ping Zhao Chandra X-ray Center / SAO, Cambridge, MA 0138 Outline 1. Reflecting Surface: Smooth or Rough?. Random Rough Surface and Power
More informationAdaptive Reconstruction Methods for Low-Dose Computed Tomography
Adaptive Reconstruction Methods for Low-Dose Computed Tomography Joseph Shtok Ph.D. supervisors: Prof. Michael Elad, Dr. Michael Zibulevsky. Technion IIT, Israel, 011 Ph.D. Talk, Apr. 01 Contents of this
More informationTomographic reconstruction of shock layer flows
Tomographic reconstruction of shock layer flows A thesis submitted for the degree of Doctor of Philosophy at The Australian National University May 2005 Rado Faletič Department of Physics Faculty of Science
More informationVariations on Image Reconstruction: Weighted Back Projection and Fourier Expectation Maximization
EG UK Computer Graphics & Visual Computing (2014) Rita Borgo, Wen Tang (Editors) Variations on Image Reconstruction: Weighted Back Projection and Fourier Expectation Maximization A. Ryan 1 and B. Mora
More informationDigital Image Fundamentals II
Digital Image Fundamentals II 1. Image modeling and representations 2. Pixels and Pixel relations 3. Arithmetic operations of images 4. Image geometry operation 5. Image processing with Matlab - Image
More informationSliced Ridgelet Transform for Image Denoising
IOSR Journal of Computer Engineering (IOSR-JCE) e-issn: 2278-0661, p- ISSN: 2278-8727Volume 14, Issue 1 (Sep. - Oct. 2013), PP 17-21 Sliced Ridgelet Transform for Image Denoising V.Krishna Naik 1, Prof.Dr.G.Manoj
More informationDigital Image Processing
Digital Image Processing SPECIAL TOPICS CT IMAGES Hamid R. Rabiee Fall 2015 What is an image? 2 Are images only about visual concepts? We ve already seen that there are other kinds of image. In this lecture
More information