Reconstruction methods for sparse-data tomography
|
|
- Baldwin Snow
- 5 years ago
- Views:
Transcription
1 Reconstruction methods for sparse-data tomography Part B: filtered back-projection Samuli Siltanen Department of Mathematics and Statistics University of Helsinki, Finland INDAM intensive period Computational Methods for Inverse Problems in Imaging Como, Italy, June 26 28, 2018
2 Lotus root tomography YouTube search: lotus tomography Video: thanks to Tatiana Bubba, Andreas Hauptmann and Juho Rimpeläinen
3 Reconstruction of a function from its line integrals was first invented by Johann Radon in 1917 f (P) = 1 π 0 df p (q) q Johann Radon ( )
4 Outline X-ray attenuation as line integral The Radon transform and its inverse Illustration of Filtered Back-Projection (FBP)
5
6
7 Wilhelm Conrad Röntgen invented X-rays and was awarded the first Nobel Prize in Physics in 1901
8 X-ray intensity attenuates inside matter, here shown with a homogeneous block
9 Formula for X-ray attenuation along a line inside homogeneous matter An X-ray with intensity I 0 enters a homogeneous physical body. I 0 I 1 } {{ } s The intensity I 1 of the X-ray when it exits the material is I 1 = I 0 e µs, where s is the length of the path of the X-ray inside the body and µ > 0 is X-ray attenuation coefficient.
10 X-ray intensity attenuates inside matter, here shown with two homogeneous blocks
11 A digital X-ray detector counts how many photons arrive at each pixel photon count 1000 X-ray source 1000 Detector
12 Adding material between the source and detector reveals the exponential X-ray attenuation law photon count
13 We take logarithm of the photon counts to compensate for the exponential attenuation law 1000 photon count 1000 log
14 Final calibration step is to subtract the logarithms from the empty space value (here 6.9) photon count log line integral
15 Formula for X-ray attenuation along a line: Beer-Lambert law Let f : [a, b] R be a nonnegative function modelling X-ray attenuation along a line inside a physical body. Beer-Lambert law connects the initial and final intensities: We can also write it in the form I 1 = I 0 e b a f (x)dx. log(i 1 /I 0 ) = b a f (x)dx, where I 0 is known from calibration and I 1 from measurement.
16 Outline X-ray attenuation as line integral The Radon transform and its inverse Illustration of Filtered Back-Projection (FBP)
17
18 Let us start with a crucial calculation To reconstruct f at a point x, the most obvious data related to f (x) are the integrals over lines passing through x. Let us sum them all together, call the result Tf (x) and see what we get by introducing polar coordinates: Tf (x) = = = = π 0 2π 0 R 2 where stands for convolution. 0 f (x + tθ)dtdθ f (x + y) dy R 2 y f (y) x y dy = (f (y) 1 y )(x), f (x + tθ) tdtdθ t
19 The Calderón operator Λ is the inverse of T Recall that Fourier transform converts convolution to multiplication (ĝ h = ĝĥ) and that 1 y (ξ) = 1 ξ. Furthermore, define the Calderón operator Λ by ( ) Λf (x) := F 1 ξ ˆf (ξ) = 1 (2π) 2 e ix ξ ξ ˆf (ξ)dξ, R 2 where F 1 is the inverse Fourier transform. Note that Λ can be thought of as a high-pass filter. Now we see that Tf (ξ) = ˆf (ξ) ξ, and thus ΛTf = f.
20
21
22 Radon transform organizes sets of line integrals Rf (θ, s) = s f (x)dx x θ=s θ
23 Definition of the Radon transform Let f (x) = f (x 1, x 2 ) be the X-ray attenuation coefficient. The classical model for tomographic data is the Radon transform Rf (θ, s) = f (x)dx = x θ=s f (sθ + y)dy, y θ θ S 1, s R, where S 1 is the unit circle, θ is the orthogonal complement of the unit vector θ and x θ denotes vector inner product. Note that f is defined on R 2 and Rf is defined on S 1 R 1. Also: Rf (θ, s) = Rf ( θ, s). θ [0, 2π].
24 Radon transform as sinogram
25 Adjoint of the Radon transform R: the back-projection operator R Use Fubini s theorem and change of variables x = sθ + y with y θ and dx = dsdy and s = θ x to get Rf (θ, s) g(θ, s) dsdθ S 1 R 1 ( ) = f (sθ + y)dy g(θ, s) dsdθ S 1 R 1 θ ( ) = f (x) g(θ, x θ) dx dθ S 1 R 2 ( ) = f (x) g(θ, x θ)dθ dx R 2 S 1 =: f (x) R g(x) dx R 2 Both f and R g are defined on R 2. Note that Tf = R Rf.
26 We arrive at the modern way of writing the filtered back-projection reconstruction formula f = ΛR Rf ( f (P) = 1 π 0 ) df p (q) q Johann Radon ( )
27 Outline X-ray attenuation as line integral The Radon transform and its inverse Illustration of Filtered Back-Projection (FBP)
28 Simple example of tomographic imaging with a double-disc target
29 The inverse problem of tomography is to recover the unknown target from the measured X-ray data
30 Since we know the projection directions we can back-project the data into the image
31 Since we know the projection directions we can back-project the data into the image
32 Since we know the projection directions we can back-project the data into the image
33 Since we know the projection directions we can back-project the data into the image
34 Since we know the projection directions we can back-project the data into the image
35 Since we know the projection directions we can back-project the data into the image
36 Since we know the projection directions we can back-project the data into the image
37 Back-projection becomes more useful by summing up the images
38 Summing all the back-projections results in a blurred reconstruction
39 Summing all the back-projections results in a blurred reconstruction
40 Summing all the back-projections results in a blurred reconstruction
41 Summing all the back-projections results in a blurred reconstruction
42 Summing all the back-projections results in a blurred reconstruction
43 Summing all the back-projections results in a blurred reconstruction
44 Summing all the back-projections results in a blurred reconstruction
45 Here we use more directions, so the reconstruction quality is higher
46 Final reconstruction involves filtering on top of the back-projection ξ fb(ξ) fb(ξ) Multiplication with ice-cream cone FFT IFFT
47 These books are recommended for learning the mathematics of practical X-ray tomography 1983 Deans: The Radon Transform and Some of Its Applications 1986 Natterer: The mathematics of computerized tomography 1988 Kak & Slaney: Principles of computerized tomographic imaging 1996 Engl, Hanke & Neubauer: Regularization of inverse problems 1998 Hansen: Rank-deficient and discrete ill-posed problems 2001 Natterer & Wübbeling: Mathematical Methods in Image Reconstruction 2008 Buzug: Computed Tomography: From Photon Statistics to Modern Cone-Beam CT 2008 Epstein: Introduction to the mathematics of medical imaging 2010 Hansen: Discrete inverse problems 2012 Mueller & S: Linear and Nonlinear Inverse Problems with Practical Applications 2014 Kuchment: The Radon Transform and Medical Imaging
48 Thank you for your attention!
Cover 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 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 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 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 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 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 informationTomographic X-ray data of a walnut
Tomographic X-ray data of a walnut arxiv:1502.04064v1 [physics.data-an] 11 Feb 2015 K. Hämäläinen, L. Harhanen, A. Kallonen, A. Kujanpää, E. Niemi and S. Siltanen February 16, 2015 Abstract This is the
More informationPolar Coordinates. Chapter 10: Parametric Equations and Polar coordinates, Section 10.3: Polar coordinates 28 / 46
Polar Coordinates Polar Coordinates: Given any point P = (x, y) on the plane r stands for the distance from the origin (0, 0). θ stands for the angle from positive x-axis to OP. Polar coordinate: (r, θ)
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 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 informationPolar Coordinates. Chapter 10: Parametric Equations and Polar coordinates, Section 10.3: Polar coordinates 27 / 45
: Given any point P = (x, y) on the plane r stands for the distance from the origin (0, 0). θ stands for the angle from positive x-axis to OP. Polar coordinate: (r, θ) Chapter 10: Parametric Equations
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 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 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 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 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 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 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 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 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 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 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 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 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 informationTotal Variation and Tomographic Imaging from Projections
Total Variation and Tomographic Imaging from Projections Per Christian Hansen & Jakob Heide Jørgensen Technical University of Denmark Joint work with Dr. Emil Sidky & Prof. Xiaochuan Pan Univ. of Chicago
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 informationTranslational Computed Tomography: A New Data Acquisition Scheme
2nd International Symposium on NDT in Aerospace 2010 - We.1.A.3 Translational Computed Tomography: A New Data Acquisition Scheme Theobald FUCHS 1, Tobias SCHÖN 2, Randolf HANKE 3 1 Fraunhofer Development
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 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 ART of tomography
The ART of tomography Timothy G. Feeman Villanova University 2016 NCAA Men s Basketball National Champions In 1979, the Nobel Prize for Medicine and Physiology was awarded jointly to Allan McLeod Cormack
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 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 informationReconstruction of Tomographic Images From Limited Projections Using TVcim-p Algorithm
Reconstruction of Tomographic Images From Limited Projections Using TVcim-p Algorithm ABDESSALEM BENAMMAR, AICHA ALLAG and REDOUANE DRAI Research Center in Industrial Technologies (CRTI), P.O.Box 64, Cheraga
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 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 informationA discrete convolution kernel for No-DC MRI
Inverse Problems Inverse Problems 3 (5) 856 (pp) doi:.88/66-56/3/8/856 A discrete convolution kernel for No-DC MRI Gengsheng L Zeng, and Ya Li 3 Department of Engineering, Weber State University, Ogden,
More informationThe Role of Linear Algebra in Computed Tomography
The Role of Linear Algebra in Computed Tomography Katarina Gustafsson katarina.gustafsson96@gmail.com under the direction of Prof. Ozan Öktem and Emer. Jan Boman Department of Mathematics KTH - Royal Institution
More informationInterior Reconstruction Using the Truncated Hilbert Transform via Singular Value Decomposition
Interior Reconstruction Using the Truncated Hilbert Transform via Singular Value Decomposition Hengyong Yu 1, Yangbo Ye 2 and Ge Wang 1 1 CT Laboratory, Biomedical Imaging Division, VT-WFU School of Biomedical
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 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 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 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 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 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 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 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 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 informationMH2800/MAS183 - Linear Algebra and Multivariable Calculus
MH28/MAS83 - Linear Algebra and Multivariable Calculus SEMESTER II EXAMINATION 2-22 Solved by Tao Biaoshuai Email: taob@e.ntu.edu.sg QESTION Let A 2 2 2. Solve the homogeneous linear system Ax and write
More informationCover Page. The handle holds various files of this Leiden University dissertation
Cover Page The handle http://hdl.handle.net/188/48289 holds various files of this Leiden University dissertation Author: Plantagie, L. Title: Algebraic filters for filtered backprojection Issue Date: 201-04-13
More informationReconstruction from Projections
Reconstruction from Projections M.C. Villa Uriol Computational Imaging Lab email: cruz.villa@upf.edu web: http://www.cilab.upf.edu Based on SPECT reconstruction Martin Šámal Charles University Prague,
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 informationA new calibration-free beam hardening reduction method for industrial CT
A new calibration-free beam hardening reduction method for industrial CT ECC 2 for industrial CT Tobias Würfl 1, Nicole Maaß 2, Frank Dennerlein 2, Andreas K. Maier 1 1 Pattern Recognition Lab, FAU Erlangen-Nürnberg;
More informationAn FDK-like cone-beam SPECT reconstruction algorithm for non-uniform attenuated
Home Search Collections Journals About Contact us My IOPscience An FK-like cone-beam SPECT reconstruction algorithm for non-uniform attenuated projections acquired using a circular trajectory This content
More informationSpline-based Sparse Tomographic Reconstruction with Besov Priors
Spline-based Sparse Tomographic Reconstruction with Besov Priors Elham Sakhaee a and Alireza Entezari a a Department of CISE, University of Florida, Gainesville, FL ABSTRACT Tomographic reconstruction
More informationTwo Local FBP Algorithms for Helical Cone-beam Computed Tomography
Digital Industrial Radiology and Computed Tomography (DIR 215) 22-25 June 215, Belgium, Ghent - www.ndt.net/app.dir215 More Info at Open Access Database www.ndt.net/?id=187 Two Local FBP Algorithms for
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 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 informationFourier transforms and convolution
Fourier transforms and convolution (without the agonizing pain) CS/CME/BioE/Biophys/BMI 279 Oct. 26, 2017 Ron Dror 1 Why do we care? Fourier transforms Outline Writing functions as sums of sinusoids The
More informationComputed Tomography (CT) Scan Image Reconstruction on the SRC-7 David Pointer SRC Computers, Inc.
Computed Tomography (CT) Scan Image Reconstruction on the SRC-7 David Pointer SRC Computers, Inc. CT Image Reconstruction Herman Head Sinogram Herman Head Reconstruction CT Image Reconstruction for all
More informationSPECT reconstruction
Regional Training Workshop Advanced Image Processing of SPECT Studies Tygerberg Hospital, 19-23 April 2004 SPECT reconstruction Martin Šámal Charles University Prague, Czech Republic samal@cesnet.cz Tomography
More informationImage Reconstruction and Image Analysis in Tomography: Fan Beam and 3D Cone Beam
Pubblicazioni del Centro De Giorgi Proceedings of the Interdisciplinary Workshop on Mathematical Methods in Biomedical Imaging and Intensity-Modulated Radiation Therapy (IMRT), Pisa, Italy, October 2007.
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 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 informationA Survey of Algebraic Algorithms in Computerized Tomography
A Survey of Algebraic Algorithms in Computerized Tomography Martin A. Brooks A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree of Modelling and Computational Science (MSc) in
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 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 informationDEVELOPMENT OF SELECTIVE TOMOGRAPHY ALGORITHM FOR CT IMAGE RECONSTRUCTION
Abstract DEVEOPMENT OF SEECTIVE TOMOGRAPHY AGORITHM FOR CT IMAGE RECONSTRUCTION R Kishore 1, B R Karthikeyan 2, Preetham Shankpal 2 1-M.Sc. [Engg.] Student, 2-Assistant Professor, Department of EEE, M.S.Ramaiah
More informationScaling Calibration in the ATRACT Algorithm
Scaling Calibration in the ATRACT Algorithm Yan Xia 1, Andreas Maier 1, Frank Dennerlein 2, Hannes G. Hofmann 1, Joachim Hornegger 1,3 1 Pattern Recognition Lab (LME), Friedrich-Alexander-University Erlangen-Nuremberg,
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 informationApplication of the topological gradient method to tomography
Application of the topological gradient method to tomography Application of the topological gradient method to tomography D. Auroux * L. Jaafar-Belaid ** B. Rjaibi ** * Institut de Mathématiques de Toulouse
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 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 informationOptimisation of projection angle selection in Computed Tomography
Universiteit Leiden Master Thesis Mathematics Optimisation of projection angle selection in Computed Tomography Author: Francien G. Bossema, BSc. Supervisors: Dr. Floske Spieksma (UL) Prof. dr. Joost Batenburg
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 informationA Comparative Study of Two Methods for the Correction of Beam Hardening Artifacts in X-ray Computed Tomography
A Comparative Study of Two Methods for the Correction of Beam Hardening Artifacts in X-ray Computed Tomography by Maryam Khalid Alarfaj submitted to Oregon State University in partial fulfillment of the
More informationPolar Coordinates. Calculus 2 Lia Vas. If P = (x, y) is a point in the xy-plane and O denotes the origin, let
Calculus Lia Vas Polar Coordinates If P = (x, y) is a point in the xy-plane and O denotes the origin, let r denote the distance from the origin O to the point P = (x, y). Thus, x + y = r ; θ be the angle
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 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 information3D RECONSTRUCTION WITH X-RAY IMAGES
3D RECONSTRUCTION WITH X-RAY IMAGES 1 Hong-Yi Chen ( 陳弘毅 ), 2 Chiou-Shann Fuh ( 傅楸善 ), Cheng-Shih Wong ( 翁丞世 ), Tzu-Chia Tung ( 董子嘉 ), You-Hsien Lee ( 李祐賢 ), Ting-Chen Tsan ( 昝亭甄 ) Dept. of Computer Science
More informationART Exhibit. Per Christian Hansen. Tommy Elfving Touraj Nikazad Hans Henrik B. Sørensen. joint work with, among others
ART Exhibit Per Christian Hansen joint work with, among others Tommy Elfving Touraj Nikazad Hans Henrik B. Sørensen About Me Forward problem Interests: numerical methods for inverse problems and tomography,
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 informationComputed tomography - outline
Computed tomography - outline Computed Tomography Systems Jørgen Arendt Jensen and Mikael Jensen (DTU Nutech) October 6, 216 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering
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 informationExterior and limited-angle tomography in non-destructive evaluation
Inverse Problems 14 (1998) 339 353. Printed in the UK PII: S0266-5611(98)87954-9 Exterior and limited-angle tomography in non-destructive evaluation Eric Todd Quinto Department of Mathematics, Tufts University,
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 informationLearning Splines for Sparse Tomographic Reconstruction. Elham Sakhaee and Alireza Entezari University of Florida
Learning Splines for Sparse Tomographic Reconstruction Elham Sakhaee and Alireza Entezari University of Florida esakhaee@cise.ufl.edu 2 Tomographic Reconstruction Recover the image given X-ray measurements
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 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 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 informationCh. 4 Physical Principles of CT
Ch. 4 Physical Principles of CT CLRS 408: Intro to CT Department of Radiation Sciences Review: Why CT? Solution for radiography/tomography limitations Superimposition of structures Distinguishing between
More informationGradient and Directional Derivatives
Gradient and Directional Derivatives MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Fall 2011 Background Given z = f (x, y) we understand that f : gives the rate of change of z in
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 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 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 informationHenry Ford NERS/BIOE 481. Lecture 11 B Computed Tomography (CT)
NERS/BIOE 481 Lecture 11 B Computed Tomography (CT) Michael Flynn, Adjunct Prof Nuclear Engr & Rad. Science mikef@umich.edu mikef@rad.hfh.edu Henry Ford Health System RADIOLOGY RESEARCH VII Computed Tomography
More informationEM+TV Based Reconstruction for Cone-Beam CT with Reduced Radiation
Appears in ISVC 2011, Part I, LNCS 6938, pp. 1-10, 2011 EM+TV Based Reconstruction for Cone-Beam CT with Reduced Radiation Ming Yan 1, Jianwen Chen 2, Luminita A. Vese 1, John Villasenor 2, Ale Bui 3,
More informationMéthodes d imagerie pour les écoulements et le CND
Méthodes d imagerie pour les écoulements et le CND Journée scientifique FED3G CEA LIST/Lab Imagerie Tomographie et Traitement Samuel Legoupil 15 juin 2012 2D/3D imaging tomography Example Petrochemical
More information