Image Reconstruction from Projection
|
|
- Easter Nichols
- 5 years ago
- Views:
Transcription
1 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 geometry processing is used to generate a three-dimensional image of the internals of an object from a large series of two-dimensional X-ray images taken around a single axis of rotation. 3/7/18 1
2 Backprojection In computed tomography or other imaging techniques requiring reconstruction from multiple projections, an algorithm for calculating the contribution of each voxel of the structure to the measured ray data, to generate an image; the oldest and simplest method of image reconstruction. 3/7/18 2
3 Medical Imaging (Before 1895) Only way to see is to cut!
4 Engineering Advances 1 st x-ray (1895) x-ray today
5 Projection X-Ray Projection Format Small Dose Fast Inexpensive Pneumonia Projection
6 Computed Tomography (CT) Head CT abdominal CT Many Projections
7 Soft, uniform tissue Image Reconstruction: Introduction Intensity is proportional to absorption Uniform with higher absorption Tumor 3/7/18 7
8 Image Reconstruction: Introduction 3/7/18 8
9 3/7/18 9
10 3/7/18 10
11 Other CTs Electron beam CT (Fifth-generation CT) Electron beam tomography (EBCT) was introduced in the early 1980s, by medical physicist Andrew Castagnini, as a method of improving the temporal resolution of CT scanners. High cost of EBCT equipment, and poor flexibility Helical (or spiral) cone beam computed tomography (sixth-generation) A type of three dimensional computed tomography (CT) in which the source (usually of x-rays) describes a helical trajectory relative to the object while a two dimensional array of detectors measures the transmitted radiation on part of a cone of rays emitting from the source 3/7/18 11
12 Other CTs Multislice CT (seventh-generation) The major benefit of multi-slice CT Significant increase in detail Utilizes X-ray tubes more economically Reducing cost and potentially reducing dosage 3/7/18 12
13 Impulse lines and line-integrals In 1D Z 1 In 2D Z f(x) (x)dx = f(0) f(x, y) (x)dx = f(0,y) A 1D Line Integral Z 1 1 Z 1 1 f(x, y) (x)dy dx= Z 1 f(0,y)dy 1 Integral of the 1D
14 Impulse lines and line-integrals y f(x,y) 1 2 x 1 dy Z f(x, y) (x)dx = f(0,y)=u(y) Z 1 dy =1 y 1 dy Z f(x, y) (x 1)dx = f(1,y)=u(2y) what about line integral with δ(x-u)? Z 1 dy y =0.5
15 Line Integral and Projection y f(x,y) 1 2 x p(u) 1 2 u p(u) = Z 1 1 Z 1 1 f(x, y) (x u)dxdy = Z 1 1 f(u,y)dy
16 Projections and the Radon Transform 3/7/18 16
17 Projections and the Radon Transform g( ρ, θ ) = f ( x, y) δ ( x cosθ + y sin θ ρ ) dxdy j k k k j 3/7/18 17
18 General Projections p(, ) = Z 1 1 Z 1 1 f(x, y) ( x cos y sin )dxdy Line integral projection x-ray f(x, y) p ( )
19 Many Projections - Tomography
20 Radon Transform p(, ) = Z 1 1 Z 1 1 f(x, y) ( x cos y sin )dxdy 0.5 y f(x,y) 0.5 y f(x,y) x x p(ρ,0) 1 p(ρ,0.4π) ρ 2 ρ
21 Radon Transform: Sinogram Also called Sinogram Impulse Sinusoid y x p(, ) 2
22 g( ρ) = 0 otherwise A r ρ ρ r 3/7/18 22
23 Computed Tomography Sinogram cross-section FBP x-ray source
24 Image Reconstruction f ( x, y) = g( x cosθ + y sin θ, θ ) θ π f ( x, y) = f ( x, y) dθ θ π f ( x, y) = f ( x, y) θ 0 θ = 0 A back-projected image formed is referred to as a laminogram 3/7/18 24
25 Examples: Laminogram 3/7/18 25
26 Projection Slice Theorem (Bracewell) f(x, y) F 1D {p(, )} = F ( cos, sine cos ) y F (f x,f y ) x p(, 0) = p(x) 1D FT F (f x, 0)
27 Projection Slice Theorem (Bracewell) f(x, y) F 1D {p(, )} = F ( cos, cos ) y F (f x,f y ) sine x p(, ) 1D FT F ( cos, cos )
28 Projection Slice Theorem (Bracewell) f(x, y) F 1D {p(, )} = F ( cos, cos ) y F (f x,f y ) sine x p(, ) 1D FT F ( cos, cos )
29 Projection Slice Theorem (Bracewell)
30 Reconstruction Using Parallel-Beam Filtered Backprojections j2 ( ux vy) f ( x, y) = π + F( u, v) e dudv Let u = wcos θ, v = w sin θ, then dudv = wdwdθ, 2π j2 π w( xcosθ + ysin θ ) f ( x, y) = F( wcos, wsin ) e wdwd 0 0 = θ θ θ 2π j2 π w( xcosθ + ysin θ ) G( w, θ ) e wdwd 0 0! G( w, θ ) = G( w, θ ) π 2 ( cos sin ) (, ) (, ) j π f x y w G w e w x θ + y θ = θ dwd 0 θ θ 3/7/18 30
31 Reconstruction Using Parallel-Beam Filtered Backprojections π j2 π w( xcosθ + ysin θ ) f ( x, y) = w G( w, θ ) e dwd 0 It s not integrable = π 0 (, ) j2π wρ w G w θ e dw d ρ = xcosθ + ysinθ θ θ Approach: Window the ramp so it becomes zero outside of a defined frequency interval. That is, a window band-limits the ramp filter. 3/7/18 31
32 Hamming / Hann Widow 2π w c + ( c 1) cos 0 w ( M 1) h( w) = M 1 0 otherwise c c = = 0.54, the function is called the Hamming window 0.5, the function is called the Han window 3/7/18 32
33 The Plot of Hamming Widow 3/7/18 33
34 Filtered Backprojection The complete, filtered backprojection (to obtain the reconstructed image f(x,y) ) is described as follows: 1. Compute the 1-D Fourier transform of each projection 2. Multiply each Fourier transform by the filter function w which has been multiplied by a suitable (e.g., Hamming) window 3. Obtain the inverse 1-D Fourier transform of each resulting filtered transform 4. Integrate (sum) all the 1-D inverse transforms from step 3 3/7/18 34
35 Implementation of Filtered Backprojection in Spatial Domain Fourier transform of the product of two frequency domain functions is equal to the convolution of the spatial representation Let s(p) denote the inverse Fourier transform of w π j2π wρ f ( x, y) = w G( w, θ ) e dw d 0 ρ = xcosθ + ysinθ = π [ ( ) (, )] s ρ g ρ θ dθ 0 ρ = xcosθ + ysinθ π = g( ρ, θ ) s( x cosθ y sin θ ρ) d ρ + dθ 0 θ 3/7/18 35
36 Partly Discrete Reconstruction Let s assume continuous angle ϴ, discrete ρ f(x, y) y F (apple x, apple y ) apple =! 2 x DTFT p [n] ) P [apple]
37 Partly Discrete Reconstruction Let s assume continuous angle ϴ, discrete ρ f(x, y) y F (apple x, apple y ) apple =! 2 x DTFT p [n] ) P [apple]
38 Reconstruction From Polar Coordinates f[n, m] = = Z 0.5 Z Z Z F (apple x, apple y )e 2 j(apple xn+apple y m) dapple x dapple y F (, )e 2 j( cos( )n+ sin( )m) d d Polar frequency data must be multiplied by ρ Also called a rho filter
39 Discrete Reconstruction Let s assume discrete angle ϴm, discrete ρ f(x, y) y F (apple apple =! x, apple y ) 2 x DFT p m [n] ) P m [l]
40 Discrete Reconstruction Let s assume discrete angle ϴm, discrete ρ f(x, y) y F (apple apple =! x, apple y ) 2 x DFT p m [n] ) P m [l]
41 Filtered Back Projection Replace integrals with sums. Sum over radius and angle Define a (filtered) backprojection: C m [n x,n y ]= So, (N/2) 1 X l= N/2 F [l, m ]e 2 j(l/n cos( m)n x +l/n sin( m )n y ) l/n f[n x,n y ]= X m C m [n x,n y ] =
42 Example
43 Example Convolution Back Projection For ϴ=0 C 0 [n x,n y ]= (N/2) 1 X l= N/2 F [l, 0] l/n e 2 j(l/nn x) FFT Rho back proj. F [l, 0] l/n IFFT
44 Example Convolution Back Projection For ϴ=π/2 C /2 [n x,n y ]= (N/2) 1 X l= N/2 F [l, /2] l/n e 2 j(l/nn y) FFT Rho back proj. F [l, /2] l/n IFFT
45 Convolution Back Projection
46 Filtered Back Projection Back projection Filtered Back projection
47 How Many Projections? FOV F (apple x, apple y ) 1/FOV
48 How Many Projections? 256 Proj. 128 Proj 64 Proj
49 Fan Beam CT Single Source Many detectors How to reconstruct?
50 Fan Beam CT Single Source Many detectors How to reconstruct? Re-binning!
51 Fan Beam CT Single Source Many detectors How to reconstruct? Re-binning!
Digital 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationFirst CT Scanner. How it Works. Contemporary CT. Before and After CT. Computer Tomography: How It Works. Medical Imaging and Pattern Recognition
Computer Tomography: How t Works Medical maging and Pattern Recognition Lecture 7 Computed Tomography Oleh Tretiak Only one plane is illuminated. Source-subject motion provides added information. 2 How
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 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 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 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 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 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 informationComputed Tomography. Principles, Design, Artifacts, and Recent Advances. Jiang Hsieh THIRD EDITION. SPIE PRESS Bellingham, Washington USA
Computed Tomography Principles, Design, Artifacts, and Recent Advances THIRD EDITION Jiang Hsieh SPIE PRESS Bellingham, Washington USA Table of Contents Preface Nomenclature and Abbreviations xi xv 1 Introduction
More informationTopic 8: Tomographic Imaging
Topic 8: Tomographic maging Aim This section covers tomographic imaging and some of its applications mainly related to medical imaging. Contents: ntroduction The adon Transform ourier nversion iltered
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 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 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 informationMath 11 Fall 2016 Section 1 Monday, October 17, 2016
Math 11 Fall 16 Section 1 Monday, October 17, 16 First, some important points from the last class: f(x, y, z) dv, the integral (with respect to volume) of f over the three-dimensional region, is a triple
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 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 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 informationEEE 512 ADVANCED DIGITAL SIGNAL AND IMAGE PROCESSING
UNIVERSITI SAINS MALAYSIA Semester I Examination Academic Session 27/28 October/November 27 EEE 52 ADVANCED DIGITAL SIGNAL AND IMAGE PROCESSING Time : 3 hours INSTRUCTION TO CANDIDATE: Please ensure that
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 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 informationComputer-Tomography II: Image reconstruction and applications
Computer-Tomography II: Image reconstruction and applications Prof. Dr. U. Oelfke DKFZ Heidelberg Department of Medical Physics (E040) Im Neuenheimer Feld 280 69120 Heidelberg, Germany u.oelfke@dkfz.de
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 informationImprovements of FDK-type Helical Cone-beam CT Imaging
Improvements of FDK-type Helical Cone-beam CT Imaging Liang Hong Zhu School of Electrical & Electronic Engineering A thesis submitted to the in fulfillment of the requirements for the degree of Doctor
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 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 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 informationDiscrete Estimation of Data Completeness for 3D Scan Trajectories with Detector Offset
Discrete Estimation of Data Completeness for 3D Scan Trajectories with Detector Offset Andreas Maier 1, Patrick Kugler 2, Günter Lauritsch 2, Joachim Hornegger 1 1 Pattern Recognition Lab and SAOT Erlangen,
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 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 informationCorso di laurea in Fisica A.A Fisica Medica 4 TC
Corso di laurea in Fisica A.A. 2007-2008 Fisica Medica 4 TC Computed Tomography Principles 1. Projection measurement 2. Scanner systems 3. Scanning modes Basic Tomographic Principle The internal structure
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 informationCHAPTER 4 RAY COMPUTATION. 4.1 Normal Computation
CHAPTER 4 RAY COMPUTATION Ray computation is the second stage of the ray tracing procedure and is composed of two steps. First, the normal to the current wavefront is computed. Then the intersection of
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 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 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 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 informationTopic 8: Tomographic Imaging
8.1 Introduction Topic 8: Tomographic Imaging Most generally tomography is the reconstruction of an image from a series of projections, and is most commonly used to image the internal structure of a semi-transparent
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 informationDigital Image Processing. Image Enhancement in the Frequency Domain
Digital Image Processing Image Enhancement in the Frequency Domain Topics Frequency Domain Enhancements Fourier Transform Convolution High Pass Filtering in Frequency Domain Low Pass Filtering in Frequency
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 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 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 informationThe Pennsylvania State University. The Graduate School. College of Engineering DEVELOPMENT OF A NEUTRON COMPUTED TOMOGRAPHY SYSTEM
The Pennsylvania State University The Graduate School College of Engineering DEVELOPMENT OF A NEUTRON COMPUTED TOMOGRAPHY SYSTEM AT THE PENNSYLVANIA STATE UNIVERISTY A Thesis in Nuclear Engineering by
More informationImprovement of Efficiency and Flexibility in Multi-slice Helical CT
J. Shanghai Jiaotong Univ. (Sci.), 2008, 13(4): 408 412 DOI: 10.1007/s12204-008-0408-x Improvement of Efficiency and Flexibility in Multi-slice Helical CT SUN Wen-wu 1 ( ), CHEN Si-ping 2 ( ), ZHUANG Tian-ge
More informationMULTI-PURPOSE 3D COMPUTED TOMOGRAPHY SYSTEM
MULTI-PURPOSE 3D COMPUTED TOMOGRAPHY SYSTEM M. Simon, C. Sauerwein, I. Tiseanu, S. Burdairon Hans Wälischmiller GmbH Klingleweg 8, D-88709 Meersburg, Germany e-mail: ms@hwm.com ABSTRACT A new flexible
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 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 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 informationX-ray Tomography. A superficial introduction, but sufficient enough to get us started in surgical navigation.
X-ray Tomography A superficial introduction, but sufficient enough to get us started in surgical navigation. X-ray absorption in homogeneous tissue I o I o / I d m = density I I=I o e -kdm k= constant
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 informationComputed Tomography January 2002 KTH A.K.
CT A.K. Computed Tomography January KTH 1 Introduction X-ray was discovered (accidentally) by a German physicist, Wilhelm Konrad Röntgen in 1895. A few years later, in 191, Röntgen was awarded the first
More informationTheoretically-exact CT-reconstruction from experimental data
Theoretically-exact CT-reconstruction from experimental data T Varslot, A Kingston, G Myers, A Sheppard Dept. Applied Mathematics Research School of Physics and Engineering Australian National University
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 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 informationCoordinate Transformations in Advanced Calculus
Coordinate Transformations in Advanced Calculus by Sacha Nandlall T.A. for MATH 264, McGill University Email: sacha.nandlall@mail.mcgill.ca Website: http://www.resanova.com/teaching/calculus/ Fall 2006,
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 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 informationComputer-Tomography I: Principles, History, Technology
Computer-Tomography I: Principles, History, Technology Prof. Dr. U. Oelfke DKFZ Heidelberg Department of Medical Physics (E040) Im Neuenheimer Feld 280 69120 Heidelberg, Germany u.oelfke@dkfz.de History
More informationParametric Surfaces. Substitution
Calculus Lia Vas Parametric Surfaces. Substitution Recall that a curve in space is given by parametric equations as a function of single parameter t x = x(t) y = y(t) z = z(t). A curve is a one-dimensional
More informationComputed Tomography. Principles of Medical Imaging. Contents. Prof. Dr. Philippe Cattin. MIAC, University of Basel. Sep 26th/Oct 3rd, 2016
Computed Tomography Principles of Medical Imaging Prof. Dr. Philippe Cattin MIAC, University of Basel Contents Abstract 1 Computed Tomography Basics Introduction Computed Tomography Hounsfield's CT Prototype
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 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 informationTo Do. Advanced Computer Graphics. Sampling and Reconstruction. Outline. Sign up for Piazza
Advanced Computer Graphics CSE 63 [Spring 207], Lecture 3 Ravi Ramamoorthi http://www.cs.ucsd.edu/~ravir Sign up for Piazza To Do Assignment, Due Apr 28. Anyone need help finding partners? Any issues with
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 informationMath Boot Camp: Coordinate Systems
Math Boot Camp: Coordinate Systems You can skip this boot camp if you can answer the following question: Staying on a sphere of radius R, what is the shortest distance between the point (0, 0, R) on the
More informationVisual Pathways to the Brain
Visual Pathways to the Brain 1 Left half of visual field which is imaged on the right half of each retina is transmitted to right half of brain. Vice versa for right half of visual field. From each eye
More informationRadiology. Marta Anguiano Millán. Departamento de Física Atómica, Molecular y Nuclear Facultad de Ciencias. Universidad de Granada
Departamento de Física Atómica, Molecular y Nuclear Facultad de Ciencias. Universidad de Granada Overview Introduction Overview Introduction Tecniques of imaging in Overview Introduction Tecniques of imaging
More information