Moscow-Bavarian Joint Advanced Student School 2006 / Medical Imaging Principles of Computerized Tomographic Imaging and Cone-Beam Reconstruction

Similar documents
BME I5000: Biomedical Imaging

MEDICAL IMAGING 2nd Part Computed Tomography

Ch. 4 Physical Principles of CT

Corso di laurea in Fisica A.A Fisica Medica 4 TC

Computer-Tomography I: Principles, History, Technology

MEDICAL IMAGING 2nd Part Computed Tomography

Biomedical Imaging. Computed Tomography. Patrícia Figueiredo IST

Introduction to Biomedical Imaging

Some reference material

Evaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT

Radon Transform and Filtered Backprojection

TEP Hounsfield units. Related topics Attenuation coefficient, Hounsfield units

X-ray tomography. X-ray tomography. Applications in Science. X-Rays. Notes. Notes. Notes. Notes

Introduction to Medical Imaging. Lecture 6: X-Ray Computed Tomography. CT number (in HU) = Overview. Klaus Mueller

Principles of Computerized Tomographic Imaging

Radiology. Marta Anguiano Millán. Departamento de Física Atómica, Molecular y Nuclear Facultad de Ciencias. Universidad de Granada

Emission Computed Tomography Notes

Digital Image Processing

Computed tomography - outline

Joint ICTP-TWAS Workshop on Portable X-ray Analytical Instruments for Cultural Heritage. 29 April - 3 May, 2013

NOVEL ALGORITHMS FOR X-RAY COMPUTED TOMOGRAPHY. Peter Lekeaka-Takunju

First CT Scanner. How it Works. Contemporary CT. Before and After CT. Computer Tomography: How It Works. Medical Imaging and Pattern Recognition

Computed Tomography January 2002 KTH A.K.

Index. aliasing artifacts and noise in CT images, 200 measurement of projection data, nondiffracting

ELEG404/604: Digital Imaging & Photography

Fundamentals of CT imaging

Image Acquisition Systems

Physical bases of X-ray diagnostics

ML reconstruction for CT

Chapter6 Image Reconstruction

Enhanced material contrast by dual-energy microct imaging

Shadow casting. What is the problem? Cone Beam Computed Tomography THE OBJECTIVES OF DIAGNOSTIC IMAGING IDEAL DIAGNOSTIC IMAGING STUDY LIMITATIONS

Diagnostic imaging techniques. Krasznai Zoltán. University of Debrecen Medical and Health Science Centre Department of Biophysics and Cell Biology

Measurement of Projection Data- The Nondiffracting Case

DUE to beam polychromacity in CT and the energy dependence

School on X-ray Imaging Techniques at the ESRF. Pierre BLEUET. ID22 Beamline

8/7/2017. Disclosures. MECT Systems Overview and Quantitative Opportunities. Overview. Computed Tomography (CT) CT Numbers. Polyenergetic Acquisition

CLASS HOURS: 4 CREDIT HOURS: 4 LABORATORY HOURS: 0

Light. Form of Electromagnetic Energy Only part of Electromagnetic Spectrum that we can really see

CT Systems and their standards

A Comparative Study of Two Methods for the Correction of Beam Hardening Artifacts in X-ray Computed Tomography

Introduction to XCT and its growth

Multi-slice CT Image Reconstruction Jiang Hsieh, Ph.D.

Medical Imaging BMEN Spring 2016

COMPARATIVE STUDIES OF DIFFERENT SYSTEM MODELS FOR ITERATIVE CT IMAGE RECONSTRUCTION

Reconstruction methods for sparse-data tomography

Spiral CT. Protocol Optimization & Quality Assurance. Ge Wang, Ph.D. Department of Radiology University of Iowa Iowa City, Iowa 52242, USA

Physics 1502: Lecture 28 Today s Agenda

Development and quantitative assessment of a beam hardening correction model for preclinical micro-ct

Computed tomography of simple objects. Related topics. Principle. Equipment TEP Beam hardening, artefacts, and algorithms

Midterm Review. Yao Wang Polytechnic University, Brooklyn, NY 11201

MEDICAL IMAGE ANALYSIS

Artifact Mitigation in High Energy CT via Monte Carlo Simulation

Tomographic Reconstruction

CT: Physics Principles & Equipment Design

Computed tomography (Item No.: P )

Image Reconstruction from Projection

Physics 210 Medical Physics Midterm Exam Fall 2012 October 12, 2012

arxiv: v2 [cond-mat.mtrl-sci] 5 Jan 2010

Empirical cupping correction: A first-order raw data precorrection for cone-beam computed tomography

What is Monte Carlo Modeling*?

Computed Tomography. Principles of Medical Imaging. Contents. Prof. Dr. Philippe Cattin. MIAC, University of Basel. Sep 26th/Oct 3rd, 2016

Tomography. Forward projectionsp θ (r) are known as a Radon transform. Objective: reverse this process to form the original image

GE s Revolution CT MATLAB III: CT. Kathleen Chen March 20, 2018

Purpose: To determine the index of refraction of glass, plastic and water.

DUAL energy CT (DECT) is a modality where one and. Empirical Dual Energy Calibration (EDEC) for Cone-Beam Computed Tomography

Computational Medical Imaging Analysis

Fast iterative beam hardening correction based on frequency splitting in computed tomography

Metal Streak Artifacts in X-ray Computed Tomography: A Simulation Study

Hot Sync. Materials Needed Today

Lecture 7 Notes: 07 / 11. Reflection and refraction

DICOM. Supplement 188 Multi-Energy CT Imaging. DICOM Working Group 21 Computed Tomography

MEDICAL EQUIPMENT: COMPUTED TOMOGRAPHY. Prof. Yasser Mostafa Kadah

Simulation of Mammograms & Tomosynthesis imaging with Cone Beam Breast CT images

The ART of tomography

Phase-Contrast Imaging and Tomography at 60 kev using a Conventional X-ray Tube

Development of a multi-axis X-ray CT for highly accurate inspection of electronic devices

Impact of X-ray Scatter When Using CT-based Attenuation Correction in PET: A Monte Carlo Investigation

Course Updates. Reminders: 1) Assignment #12 due today. 2) Polarization, dispersion. 3) Last HW (#13 posted) due Monday, May 3rd

Experiment 9. Law of reflection and refraction of light

UNIVERSITA DEGLI STUDI DI BOLOGNA

Effect of Scattering on the Image. Reducing Compton Scatter with a Grid

Corso di laurea in Fisica A.A Fisica Medica 5 SPECT, PET

Arion: a realistic projection simulator for optimizing laboratory and industrial micro-ct

Energy resolved X-ray diffraction Cl. J.Kosanetzky, G.Harding, U.Neitzel

Central Slice Theorem

RADIOLOGY AND DIAGNOSTIC IMAGING

Applying Hounsfield unit density calibration in SkyScan CT-analyser

Photon counting spectral CT versus conventional CT: comparative evaluation for breast imaging application

Medical Image Analysis, Assignment 3

Improved compressed sensing algorithm for sparse-view CT

Wavefronts and Rays. When light or other electromagnetic waves interact with systems much larger than the wavelength, it s a good approximation to

CALIBRATIONS FOR ANALYZING INDUSTRIAL SAMPLES ON MEDICAL CT SCANNERS

Design and performance characteristics of a Cone Beam CT system for Leksell Gamma Knife Icon

SYSTEM LINEARITY LAB MANUAL: 2 Modifications for P551 Fall 2013 Medical Physics Laboratory

Multilevel Optimization for Multi-Modal X-ray Data Analysis

Advanced Multi Material Decomposition of Dual Energy in Computed Tomography Image

Centre for Mathematical Sciences, Mathematics, LTH, December 2015

A New Approach for the Enhancement. of Dual-energy Computed Tomography Images. Kyung Kook Park

Digital Scatter Removal in Mammography to enable Patient Dose Reduction

Transcription:

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: P θ ( t) = ( θ, t) line ( x, A set of line integrals form a projection f y) ds

1 st generation tomographs Rotation-translation pencil beam Pencil beam, one detector For one angle a set of parallel line integrals is taken and form a projection. Then the angle is changed.

2 nd generation tomographs Rotation-translation fan beam Fan beam (10 ), ~30 detectors Multiple line integrals along a fan are taken simultaneously. Several parallel steps are done for different angles

3 rd generation tomographs Rotation-rotation single slices Fan beam (40-60 ), up to 1000 detectors A projection used only line integrals following a fan. All line integrals for one projection are taken simultaneously No translation is necessary

4 th generation tomographs Rotation-fix closed detector array Fan beam (40-60 ), up to 5000 fixed detectors Rotation only, a translation is not necessary

Purpose of tomography Conventional X-Ray provides only projections no spatial information averaging of all slices => low contrast Tomography tries to revert the projection

What is measured? The beam intensity and therewith the loss of beam intensity is measured Loss of intensity due to: Photoelectric absorption Compton effects Pair production (PET) Loss is energy-dependent

What is measured? μ is the attenuation coefficient of a material representing the photon loss rate due to the Compton effect and photoelectric absorption. μ is given by N 1 = µ N x

What is measured? The intensity after passing through x is given by:,whereas μ is assumed to be constant. As x N( x + x) = N( x) µ N( x) x goes to zero we obtain lim x 0 N( x + x) x N( x) = dn dx = µ N( x)

What is measured? Integration both sides,we obtain dn N( x) = µ dx ln N = µ x + C

What is measured? The number of photons as function of the position is given by: N( x) = Nin exp µ x [ ],where N in is the number of photons entering the object. This is the Lambert-Beer law This is only true for an x-ray beam consisting of monochromatic photons and a constant μ.

What is measured? μ(x,y) denotes the attenuation coefficient of a body The number of exiting photons is given by: N = d Nin exp µ( x, y) ds ray This is only true for an x-ray beam consisting of monochromatic photons.

Monochromatic and polychromatic x-ray Monochromatic: X-ray beam consists only of photons with the same energy (in practice this type of x-ray is not used) Polychromatic: X-ray consists of a spectrum of photons with different energy N [ ( x, y, E ds ] = S ( E)exp µ de d in ) Beam hardening => Artifacts

Detectors Xenon ionization detectors X-ray photons enter the detector chamber and ionize gas. The resulting current is measured. Used in 3rd generation tomographs

Detectors Scintillation detectors: Using a crystal the x-ray is transformed into photons with longer wavelength. These are measured using photo diodes. Used in 4th generation tomographs

Hounsfield scale Defined by Sir Godfrey Newbold Hounsfield 0 Hounsfield Units (HU) defined as radiodensity of distilled water -1000 HU defined as radiodensity of air Corresponds to the linear attenuation coefficient by: H = µ µ µ water water 1000 Substance Bone Calcification Congealed blood Grey matter White matter Water Fat Air Approx. value 80-1000 80-1000 56-76 36-46 22-32 0-100 -1000