Performance study of a fan beam collimator for a multi- modality small animal imaging device

Size: px
Start display at page:

Download "Performance study of a fan beam collimator for a multi- modality small animal imaging device"

Transcription

1 Performance study of a fan beam collimator for a multi- modality small animal imain device ASM SABBIR AHMED 1 Wolfrad Semmler 2 Jor Peter 2 1 University of Saskatchewan Saskatoon Canada & 2 German Cancer Research Centre University of Heidelber Germany Trabzon Turkey May

2 Study Objectives Desin of a fan beam collimator for a small animal imain device to perform rapid dynamic study with improved efficiency and spatial resolution. A Monte Carlo model to optimize the collimator s parameters to conduct a performance study by analyzin the scatterin and noise indexes for a 3D eometry. A S Ahmed Trabzon Turkey May Slide:*

3 Study Importance Pharmacokinetic data distribution revels important functional information of the livin bodies. Hih sensitive collimator is essentially required for rapid dynamic acquisition where hiher count rates are to be recorded in short time frame. A S Ahmed Trabzon Turkey May Slide: 3

4 Introduction Collimator s role in scintillation imain * The desired collimator material has hih linear attenuation coefficient; e.. tunsten lead tantalum etc. LLNL * Most desirable eometrically collimated rays pass the collimator without interactin the septal wall JAERI * Undesirable penetratin amma rays pass throuh one or more septal walls and scattered photons deflected by scatterin. Scintillation imain without proper collimation spoils the ultimate objective * Penetratin and scattered photons produce noises and substantially derade imae quality A S Ahmed Trabzon Turkey May Slide: 4

5 Introduction The small animal imain device Patent: EP26/ ; J. Peter et.al. IEEE 25 A S Ahmed Trabzon Turkey May Slide: 5

6 Introduction Theoretical formulation for modelin efficiency and resolution [Scatterin/attenuation/penetration=] A photon from source s z will pass throuh the collimator if the photon passes within front and back projections of the hole on imae plane. s As a collimator is composed of repeated lattice structure with respect to a source so the collimator response function is periodic with respect to source position and it depends on front and back aperture function: i.e. c r s z j k zbl hs r r z ' jk zb hs r r z ' jk A S Ahmed Trabzon Turkey May Slide: 6

7 Introduction Theory: efficiency and resolution If the collimator response function is decomposed in a Fourier series the above equation becomes N M MN m MN s i z r c z s r c ]. exp[ 2 where M and N are referred to as the hole-pattern harmonics. For a parallel hole collimator the point source response function is independent of the source position s; In that case all the harmonics except c vanish. Collimator eometric transfer function can be written in terms of the photon fluence at the collimator face and the two collimator-aperture functions where is the effective lenth of the collimator =d+t 2 is the area of a lattice. Hence hole pattern harmonics except M=N= can be nelected. N M MN m MN s i z c z z s r c z r z s 2 ]. exp[2 4 1 Takin the Fourier transform N M MN MN MN e m s i m l B z H m l l B z H l 2 ]. exp[ l e A S Ahmed Trabzon Turkey May Slide: 7

8 Introduction Theory: efficiency and resolution The Averae efficiency for the parallel-hole collimator p can be written as p H H 2 4l e Hole _ area 2 4l Lattice _ area e Total lattice area >FWHM The resolution of a parallel-hole collimator measures the width of the point source response function PSRF as a function of r. The width of the PSRF is stronly dependent on the source position in front of the collimator. The eometric resolution of the parallel-hole collimator R is p defined as the FWHM of the PSRF by the relation R p zd l e d l e B l b Henkin RE et.al. Nucl Med Mosby Ins 1996 A S Ahmed Trabzon Turkey May Slide: 8

9 Introduction Theory: : Fan beam efficiency and resolution For fan beam eometry individual hole-axis alins differently at every position on the collimator face and hence the efficiency of the fan beam collimator depends stronly on source position. The eometric efficiency for the fan beam collimator is iven as f b p [cos 2 ] f f b 2 wherez B l b is object distance p is efficiency of the parallel collimator q is the anle between collimator axis and the object position. The distance b varies from to f; i.e. <b<f. The resolution of the fan beam collimator R f R p 1 cos le / 2 b f b [Tsui BMW IEEE 1996] A S Ahmed Trabzon Turkey May Slide: 9

10 Introduction Theory: : Equation derivation for line and point source If a line source of lenth L is placed in air the eometric efficiencies for a parallel hole and a fan beam collimator are : p l p L b b f l f dz L b L 1 cos 2 2 L b f b f f f L l p No scatterin or attenuation included For attenuatin medium: exp h h p a p exp h b h b f a f L p L b b p al p e L dz b z L L 1 ] exp[ 1 [Point source] [Line source] B e B BE A e A AE e f L h L B i A i A p al f cos 2 2 No scatterin included A S Ahmed Trabzon Turkey May Slide: 1

11 Methodoloy The small animal imain device: Collimator s Specifications The scintillation crystal consists of 66x66 array.2 cm separation of square shaped individual NaITl crystal elements.13 x.13 x.6 cm 3. The proposed collimator is desined to match with the crystal dimension. Tunsten material attenuation coefficient m=36.34 cm -1 is proposed for this collimator. Thickness t=.2 cm Hole size d=.13 cm focal lenth f=18 cm. A S Ahmed Trabzon Turkey May Slide: 11

12 Methodoloy Collimator Desin consideration The proposed fan beam collimator was desined to study the whole body dynamic study. The focal lenth f was calculated as 18 cm. A hole lenth l =1.8 cm was calculated for the new converin collimator. A fan beam is a combination of parallel and converin collimator. Collimator s hole axis s are convered toward the object short axis and parallel to object lon axis. A S Ahmed Trabzon Turkey May Slide: 12

13 Methodoloy Collimator Desin A S Ahmed Trabzon Turkey May Slide: 13

14 Methodoloy Calculation of Scatterin index SI Monte Carlo method MCNP was used to calculate the fluence for individual detector-element for point and line sources. This fluence included scattered and unscattered photons emitted from the source. F4 tally calculated the number of photons over a cell that was normalized per incident particle. F4 tally r E t vn r E t Same number of particles were tallied for all cases. r E t dedt V E T dv V wvt V w V T L Scatterin index SI= Total fluence unscattered fluence Particle density N = Particle weiht/volume w/v T L = vt = Track lenth in the medium 17x17 detector blocks were included in the model study. The sources were biased within the solid anle. A S Ahmed Trabzon Turkey May Slide: 14

15 Methodoloy Anular dependency of 3D object A water filled cylindrical object r=3 cm that compares to standard mouse was placed at a distance b from collimator face. Because of symmetry the upper half A B C was only considered where A and B exist in the most and the least depth points inside attenuatin medium C in the middle. For every q anle at the object center a correspondin q anle was calculated at the focal point. q varied from f to +f -f<q<f. A S Ahmed Trabzon Turkey May Slide: 15

16 Methodoloy Calculation of Sinal noise index SNI A rin source radius 1.5 cm consists of four point sources were placed inside water filled cylinder at different depths. All sources make same manitude of anle with the collimator axis. The shortest distance between two consecutive peaks was lon enouh not to interrupt each other. The sinal at the rin-center can be taken as the averae of the sinals at j k m n with four sources at these locations. The noise at the center is calculated as the noise calculated at j k m n with a point source at the center c. s c SNI = N c s s s s S j k m n c ; sources j k. m. n 4 N N N N N j k m n c ; source c 4 A S Ahmed Trabzon Turkey May Slide: 16

17 Results & Discussion Point source profile for the collimator A point source in air at 1 cm away from the collimator face. Focal lenth f=18 cm b=1 cm A S Ahmed Trabzon Turkey May Slide: 17

18 Results & Discussion Point and line source in air For point source the detector efficiency increases 1.2E+3 times near the focal point. For line source it s increases 3.3E+1 times A S Ahmed Trabzon Turkey May Slide: 18

19 Results & Discussion Point and line source in water At.5 cm depth scatterin 14 % for parallel 19 % for fan beam At.5 cm depth scatterin 19 % for parallel 3 % for fan beam A S Ahmed Trabzon Turkey May Slide: 19

20 Results & Discussion Relative efficiency and Scatterin index for point source in water The attenuatin medium was fixed at 3 cm away from the collimator face. For point source relative efficiency increases about 6 times than that of Parallel The fan beam scatterin index increases about 11 times than that of parallel collimator A S Ahmed Trabzon Turkey May Slide: 2

21 Results & Discussion Relative efficiency and Scatterin index for line source in water The relative efficiency is linear for line source The scatterin index is linear for line source A S Ahmed Trabzon Turkey May Slide: 21

22 Results & Discussion Anular dependency of a 3D object A S Ahmed Trabzon Turkey May Slide: 22

23 Results & Discussion Anular dependency of a 3D object A S Ahmed Trabzon Turkey May Slide: 23

24 Results & Discussion Sinal noise index SNI With fan beam the SNI increases about 2.3 times than the that of the parallel collimator A S Ahmed Trabzon Turkey May Slide: 26

25 Conclusion * The performance characteristic of the proposed fan beam collimator has been studied for different source confiurations based on analytic formulation and Monte Carlo simulation studies. * The results shown in this study revealed that a fan beam collimator is useful for scintillation imain. But the advantae of fan beam is not predictable from mathematical derivation with a point source in air alone. Details particle transport calculation is essential. * Its performance stronly dependent on the source distribution and differs sinificantly for different source confiuration. * Inside attenuatin medium the increased attenuatin effect outweihs the fan beam increased eometric efficiency when a line source moves toward the focal point Supported by: Tsui BMW IEEE 1996; brain imain. * The performance of a fan beam collimator should be evaluated for specific application and purpose. A S Ahmed Trabzon Turkey May Slide: 25

26 A S Ahmed Trabzon Turkey May QUESTIONS?

Chapter 5 THE MODULE FOR DETERMINING AN OBJECT S TRUE GRAY LEVELS

Chapter 5 THE MODULE FOR DETERMINING AN OBJECT S TRUE GRAY LEVELS Qian u Chapter 5. Determinin an Object s True Gray evels 3 Chapter 5 THE MODUE OR DETERMNNG AN OJECT S TRUE GRAY EVES This chapter discusses the module for determinin an object s true ray levels. To compute

More information

I. INTRODUCTION. Figure 1. Radiation room model at Dongnai General Hospital

I. INTRODUCTION. Figure 1. Radiation room model at Dongnai General Hospital International Journal of Computational Engineering Research Vol, 04 Issue, 4 Simulation of Photon and Electron dose distributions 5 code for the treatment area using the linear electron accelerator (LINAC)

More information

Application of MCNP Code in Shielding Design for Radioactive Sources

Application of MCNP Code in Shielding Design for Radioactive Sources Application of MCNP Code in Shielding Design for Radioactive Sources Ibrahim A. Alrammah Abstract This paper presents three tasks: Task 1 explores: the detected number of as a function of polythene moderator

More information

Determination of Three-Dimensional Voxel Sensitivity for Two- and Three-Headed Coincidence Imaging

Determination of Three-Dimensional Voxel Sensitivity for Two- and Three-Headed Coincidence Imaging IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 50, NO. 3, JUNE 2003 405 Determination of Three-Dimensional Voxel Sensitivity for Two- and Three-Headed Coincidence Imaging Edward J. Soares, Kevin W. Germino,

More information

Basics of treatment planning II

Basics of treatment planning II Basics of treatment planning II Sastry Vedam PhD DABR Introduction to Medical Physics III: Therapy Spring 2015 Dose calculation algorithms! Correction based! Model based 1 Dose calculation algorithms!

More information

Breaking Through the Barriers to GPU Accelerated Monte Carlo Particle Transport

Breaking Through the Barriers to GPU Accelerated Monte Carlo Particle Transport Breaking Through the Barriers to GPU Accelerated Monte Carlo Particle Transport GTC 2018 Jeremy Sweezy Scientist Monte Carlo Methods, Codes and Applications Group 3/28/2018 Operated by Los Alamos National

More information

Medical Physics Research Center, Mashhad University of Medical Sciences, Mashhad, Iran.

Medical Physics Research Center, Mashhad University of Medical Sciences, Mashhad, Iran. DXRaySMCS First User Friendly Interface Developed for Prediction of Diagnostic Radiology X-Ray Spectra Produced by Monte Carlo (MCNP-4C) Simulation in Iran M.T. Bahreyni Toosi a*, H. Moradi b, H. Zare

More information

Artifact Mitigation in High Energy CT via Monte Carlo Simulation

Artifact Mitigation in High Energy CT via Monte Carlo Simulation PIERS ONLINE, VOL. 7, NO. 8, 11 791 Artifact Mitigation in High Energy CT via Monte Carlo Simulation Xuemin Jin and Robert Y. Levine Spectral Sciences, Inc., USA Abstract The high energy (< 15 MeV) incident

More information

Monte Carlo simulations. Lesson FYSKJM4710 Eirik Malinen

Monte Carlo simulations. Lesson FYSKJM4710 Eirik Malinen Monte Carlo simulations Lesson FYSKJM4710 Eirik Malinen MC simulations 1 Simulations of stochastic processes Interactions are stochastic: the path of a single ionizing particle may not be predicted Interactions

More information

Iterative SPECT reconstruction with 3D detector response

Iterative SPECT reconstruction with 3D detector response Iterative SPECT reconstruction with 3D detector response Jeffrey A. Fessler and Anastasia Yendiki COMMUNICATIONS & SIGNAL PROCESSING LABORATORY Department of Electrical Engineering and Computer Science

More information

Monte Carlo simulations

Monte Carlo simulations MC simulations Monte Carlo simulations Eirik Malinen Simulations of stochastic processes Interactions are stochastic: the path of a single ioniing particle may not be predicted Interactions are quantified

More information

NKS GammaUser

NKS GammaUser Uncertainties Marie-Christine Lépy Laboratoire National Henri Becquerel CEA Saclay, F-91191 Gif-sur-Yvette cedex, France E-mail: marie-christine.lepy@cea.fr NKS GammaUser 2014 Helsinki 6-8 October 2014

More information

MCNP4C3-BASED SIMULATION OF A MEDICAL LINEAR ACCELERATOR

MCNP4C3-BASED SIMULATION OF A MEDICAL LINEAR ACCELERATOR Computational Medical Physics Working Group Workshop II, Sep 3 Oct 3, 7 University of Florida (UF), Gainesville, Florida USA on CD-ROM, American Nuclear Society, LaGrange Park, IL (7) MCNP4C3-BASED SIMULATION

More information

Ray Optics I. Last time, finished EM theory Looked at complex boundary problems TIR: Snell s law complex Metal mirrors: index complex

Ray Optics I. Last time, finished EM theory Looked at complex boundary problems TIR: Snell s law complex Metal mirrors: index complex Phys 531 Lecture 8 20 September 2005 Ray Optics I Last time, finished EM theory Looked at complex boundary problems TIR: Snell s law complex Metal mirrors: index complex Today shift gears, start applying

More information

Supplementary Figure 1 Optimum transmissive mask design for shaping an incident light to a desired

Supplementary Figure 1 Optimum transmissive mask design for shaping an incident light to a desired Supplementary Figure 1 Optimum transmissive mask design for shaping an incident light to a desired tangential form. (a) The light from the sources and scatterers in the half space (1) passes through the

More information

A SYSTEM OF DOSIMETRIC CALCULATIONS

A SYSTEM OF DOSIMETRIC CALCULATIONS A SYSTEM OF DOSIMETRIC CALCULATIONS INTRODUCTION Dose calculation based on PDD and TAR have Limitations The dependence of PDD on SSD Not suitable for isocentric techniques TAR and SAR does not depend on

More information

Deep Scatter Estimation (DSE): Feasibility of using a Deep Convolutional Neural Network for Real-Time X-Ray Scatter Prediction in Cone-Beam CT

Deep Scatter Estimation (DSE): Feasibility of using a Deep Convolutional Neural Network for Real-Time X-Ray Scatter Prediction in Cone-Beam CT Deep Scatter Estimation (DSE): Feasibility of using a Deep Convolutional Neural Network for Real-Time X-Ray Scatter Prediction in Cone-Beam CT Joscha Maier 1,2, Yannick Berker 1, Stefan Sawall 1,2 and

More information

What is Monte Carlo Modeling*?

What is Monte Carlo Modeling*? What is Monte Carlo Modeling*? Monte Carlo Modeling is a statisitcal method used here to simulate radiative transfer by simulating photon (or more exactly light rays/beams) interaction with a medium. MC

More information

CHAPTER 11 NUCLEAR MEDICINE IMAGING DEVICES

CHAPTER 11 NUCLEAR MEDICINE IMAGING DEVICES CHAPTER 11 M.A. LODGE, E.C. FREY Russell H. Morgan Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, Maryland, United States of America 11.1. INTRODUCTION Imaging

More information

Monte Carlo Simulation for Neptun 10 PC Medical Linear Accelerator and Calculations of Electron Beam Parameters

Monte Carlo Simulation for Neptun 10 PC Medical Linear Accelerator and Calculations of Electron Beam Parameters Monte Carlo Simulation for Neptun 1 PC Medical Linear Accelerator and Calculations of Electron Beam Parameters M.T. Bahreyni Toossi a, M. Momen Nezhad b, S.M. Hashemi a a Medical Physics Research Center,

More information

ISOCS Characterization of Sodium Iodide Detectors for Gamma-Ray Spectrometry

ISOCS Characterization of Sodium Iodide Detectors for Gamma-Ray Spectrometry ISOCS Characterization of Sodium Iodide Detectors for Gamma-Ray Spectrometry Sasha A. Philips, Frazier Bronson, Ram Venkataraman, Brian M. Young Abstract--Activity measurements require knowledge of the

More information

The Monte Carlo simulation of a Package formed by the combination of three scintillators: Brillance380, Brillance350, and Prelude420.

The Monte Carlo simulation of a Package formed by the combination of three scintillators: Brillance380, Brillance350, and Prelude420. EURONS I3 506065 JRA9 RHIB Report made during stay IEM-CSIC Madrid december 2006 MINISTERIO DE ASUNTOS EXTERIORES Y DE COOPERACIÓN AECI VICESECRETARÍA GENERAL The Monte Carlo simulation of a Package formed

More information

Assesing multileaf collimator effect on the build-up region using Monte Carlo method

Assesing multileaf collimator effect on the build-up region using Monte Carlo method Pol J Med Phys Eng. 2008;14(3):163-182. PL ISSN 1425-4689 doi: 10.2478/v10013-008-0014-0 website: http://www.pjmpe.waw.pl M. Zarza Moreno 1, 2, N. Teixeira 3, 4, A. P. Jesus 1, 2, G. Mora 1 Assesing multileaf

More information

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

arxiv: v2 [cond-mat.mtrl-sci] 5 Jan 2010 Gamma scattering scanning of concrete block for detection of voids. Shivaramu 1, Arijit Bose 2 and M. Margret 1 1 Radiological Safety Division, Safety Group, IGCAR, Kalpakaam - 63 12 (India) 2 Chennai

More information

Development and Verification of an SP 3 Code Using Semi-Analytic Nodal Method for Pin-by-Pin Calculation

Development and Verification of an SP 3 Code Using Semi-Analytic Nodal Method for Pin-by-Pin Calculation Journal of Physical Science and Application 7 () (07) 0-7 doi: 0.765/59-5348/07.0.00 D DAVID PUBLISHIN Development and Verification of an SP 3 Code Usin Semi-Analytic Chuntao Tan Shanhai Nuclear Enineerin

More information

Spectral analysis of non-stationary CT noise

Spectral analysis of non-stationary CT noise Spectral analysis of non-stationary CT noise Kenneth M. Hanson Los Alamos Scientific Laboratory Int. Symposium and Course on Computed Tomography, Las Vegas, April 7-11, 1980 This presentation available

More information

Limitations in the PHOTON Monte Carlo gamma transport code

Limitations in the PHOTON Monte Carlo gamma transport code Nuclear Instruments and Methods in Physics Research A 480 (2002) 729 733 Limitations in the PHOTON Monte Carlo gamma transport code I. Orion a, L. Wielopolski b, * a St. Luke s/roosevelt Hospital, Columbia

More information

Engineered Diffusers Intensity vs Irradiance

Engineered Diffusers Intensity vs Irradiance Engineered Diffusers Intensity vs Irradiance Engineered Diffusers are specified by their divergence angle and intensity profile. The divergence angle usually is given as the width of the intensity distribution

More information

Iterative and analytical reconstruction algorithms for varying-focal-length cone-beam

Iterative and analytical reconstruction algorithms for varying-focal-length cone-beam Home Search Collections Journals About Contact us My IOPscience Iterative and analytical reconstruction algorithms for varying-focal-length cone-beam projections This content has been downloaded from IOPscience.

More information

BME I5000: Biomedical Imaging

BME 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 information

Outline. Monte Carlo Radiation Transport Modeling Overview (MCNP5/6) Monte Carlo technique: Example. Monte Carlo technique: Introduction

Outline. Monte Carlo Radiation Transport Modeling Overview (MCNP5/6) Monte Carlo technique: Example. Monte Carlo technique: Introduction Monte Carlo Radiation Transport Modeling Overview () Lecture 7 Special Topics: Device Modeling Outline Principles of Monte Carlo modeling Radiation transport modeling with Utilizing Visual Editor (VisEd)

More information

4. A bulb has a luminous flux of 2400 lm. What is the luminous intensity of the bulb?

4. A bulb has a luminous flux of 2400 lm. What is the luminous intensity of the bulb? 1. Match the physical quantities (first column) with the units (second column). 4. A bulb has a luminous flux of 2400 lm. What is the luminous intensity of the bulb? (π=3.) Luminous flux A. candela Radiant

More information

Crystal Quality Analysis Group

Crystal Quality Analysis Group Crystal Quality Analysis Group Contents Contents 1. Overview...1 2. Measurement principles...3 2.1 Considerations related to orientation and diffraction conditions... 3 2.2 Rocking curve measurement...

More information

Transitioning from pencil beam to Monte Carlo for electron dose calculations

Transitioning from pencil beam to Monte Carlo for electron dose calculations Transitioning from pencil beam to Monte Carlo for electron dose calculations Jessie Huang-Vredevoogd (jyhuang4@wisc.edu) University of Wisconsin NCC AAPM October 12, 2019 1 Topics to cover Background RayStation

More information

Rectangular Lenslet Array

Rectangular Lenslet Array Rectangular Lenslet Array INTRODUCTION Lenslet arrays are used in a variety of applications that include beam homogenization. This knowledge base article demonstrates the setup of an imaging lenslet array

More information

Particle track plotting in Visual MCNP6 Randy Schwarz 1,*

Particle track plotting in Visual MCNP6 Randy Schwarz 1,* Particle track plotting in Visual MCNP6 Randy Schwarz 1,* 1 Visual Editor Consultants, PO Box 1308, Richland, WA 99352, USA Abstract. A visual interface for MCNP6 has been created to allow the plotting

More information

Electron Dose Kernels (EDK) for Secondary Particle Transport in Deterministic Simulations

Electron Dose Kernels (EDK) for Secondary Particle Transport in Deterministic Simulations Electron Dose Kernels (EDK) for Secondary Particle Transport in Deterministic Simulations A. Al-Basheer, G. Sjoden, M. Ghita Computational Medical Physics Team Nuclear & Radiological Engineering University

More information

ACCURATE TEXTURE MEASUREMENTS ON THIN FILMS USING A POWDER X-RAY DIFFRACTOMETER

ACCURATE TEXTURE MEASUREMENTS ON THIN FILMS USING A POWDER X-RAY DIFFRACTOMETER ACCURATE TEXTURE MEASUREMENTS ON THIN FILMS USING A POWDER X-RAY DIFFRACTOMETER MARK D. VAUDIN NIST, Gaithersburg, MD, USA. Abstract A fast and accurate method that uses a conventional powder x-ray diffractometer

More information

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

Corso di laurea in Fisica A.A Fisica Medica 5 SPECT, PET Corso di laurea in Fisica A.A. 2007-2008 Fisica Medica 5 SPECT, PET Step 1: Inject Patient with Radioactive Drug Drug is labeled with positron (β + ) emitting radionuclide. Drug localizes

More information

Simulation of Radiographic Testing for POD Assessment

Simulation of Radiographic Testing for POD Assessment 4th European-American Workshop on Reliability of NDE - Th.5.B.2 Simulation of Radiographic Testing for POD Assessment Gerd-Rüdiger JAENISCH, Carsten BELLON, Uwe EWERT, BAM, Berlin, Germany Abstract. Efficient

More information

Recognition and Measurement of Small Defects in ICT Testing

Recognition and Measurement of Small Defects in ICT Testing 19 th World Conference on Non-Destructive Testing 2016 Recognition and Measurement of Small Defects in ICT Testing Guo ZHIMIN, Ni PEIJUN, Zhang WEIGUO, Qi ZICHENG Inner Mongolia Metallic Materials Research

More information

Figure 1: Derivation of Bragg s Law

Figure 1: Derivation of Bragg s Law What is Bragg s Law and why is it Important? Bragg s law refers to a simple equation derived by English physicists Sir W. H. Bragg and his son Sir W. L. Bragg in 1913. This equation explains why the faces

More information

Deep Scatter Estimation (DSE): Accurate Real-Time Scatter Estimation for X-Ray CT using a Deep Convolutional Neural Network

Deep Scatter Estimation (DSE): Accurate Real-Time Scatter Estimation for X-Ray CT using a Deep Convolutional Neural Network Deep Scatter Estimation (DSE): Accurate Real-Time Scatter Estimation for X-Ray CT using a Deep Convolutional Neural Network Joscha Maier 1,2, Stefan Sawall 1,2, and Marc Kachelrieß 1,2 1 German Cancer

More information

Basic Radiation Oncology Physics

Basic Radiation Oncology Physics Basic Radiation Oncology Physics T. Ganesh, Ph.D., DABR Chief Medical Physicist Fortis Memorial Research Institute Gurgaon Acknowledgment: I gratefully acknowledge the IAEA resources of teaching slides

More information

EGS5 Monte Carlo Code: Installation

EGS5 Monte Carlo Code: Installation EGS5 Monte Carlo Code: Installation What do you need to install? EGS5 Cygwin CGView gvim (or any editor) EGS5 Monte Carlo Code: Installation Installing CYGWIN Installing EGS5 Running an EGS5 TEST Any trouble

More information

18th World Conference on Nondestructive Testing, April 2012, Durban, South Africa

18th World Conference on Nondestructive Testing, April 2012, Durban, South Africa 18th World Conference on Nondestructive Testing, 16-20 April 2012, Durban, South Africa TOWARDS ESTABLISHMENT OF STANDARDIZED PRACTICE FOR ASSESSMENT OF SPATIAL RESOLUTION AND CONTRAST OF INTERNATIONAL

More information

NON-COLLIMATED SCATTERED RADIATION TOMOGRAPHY

NON-COLLIMATED SCATTERED RADIATION TOMOGRAPHY NON-COLLIMATED SCATTERED RADIATION TOMOGRAPHY Gorshkov V.A., Space Research Institute, Moscow, Russia Yumashev V.M., State corporation "Rosatom", Centre "Atom-innovation", Moscow, Russia Kirilenko K.V.,

More information

Improved Detector Response Characterization Method in ISOCS and LabSOCS

Improved Detector Response Characterization Method in ISOCS and LabSOCS P Improved Detector Response Characterization Method in ISOCS and LabSOCS *1 1 1 1 1 R. VenkataramanP P, F. BronsonP P, V. AtrashkevichP P, M. FieldP P, and B.M. YoungP P. 1 PCanberra Industries, 800 Research

More information

PHYSICS. Chapter 34 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT

PHYSICS. Chapter 34 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT PHYSICS FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E Chapter 34 Lecture RANDALL D. KNIGHT Chapter 34 Ray Optics IN THIS CHAPTER, you will learn about and apply the ray model of light Slide 34-2

More information

P H Y L A B 1 : G E O M E T R I C O P T I C S

P H Y L A B 1 : G E O M E T R I C O P T I C S P H Y 1 4 3 L A B 1 : G E O M E T R I C O P T I C S Introduction Optics is the study of the way light interacts with other objects. This behavior can be extremely complicated. However, if the objects in

More information

Emission Computed Tomography Notes

Emission 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 information

Summary. Introduction

Summary. Introduction Imae amplitude compensations: propaator amplitude recovery and acquisition aperture correction Jun Cao and Ru-Shan Wu Department of Earth and Planetary Sciences/IGPP, University of California, Santa Cruz

More information

Implementation and evaluation of a fully 3D OS-MLEM reconstruction algorithm accounting for the PSF of the PET imaging system

Implementation and evaluation of a fully 3D OS-MLEM reconstruction algorithm accounting for the PSF of the PET imaging system Implementation and evaluation of a fully 3D OS-MLEM reconstruction algorithm accounting for the PSF of the PET imaging system 3 rd October 2008 11 th Topical Seminar on Innovative Particle and Radiation

More information

Medical Imaging BMEN Spring 2016

Medical Imaging BMEN Spring 2016 Name Medical Imaging BMEN 420-501 Spring 2016 Homework #4 and Nuclear Medicine Notes All questions are from the introductory Powerpoint (based on Chapter 7) and text Medical Imaging Signals and Systems,

More information

Image Acquisition Systems

Image 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 information

ELECTRON DOSE KERNELS TO ACCOUNT FOR SECONDARY PARTICLE TRANSPORT IN DETERMINISTIC SIMULATIONS

ELECTRON DOSE KERNELS TO ACCOUNT FOR SECONDARY PARTICLE TRANSPORT IN DETERMINISTIC SIMULATIONS Computational Medical Physics Working Group Workshop II, Sep 30 Oct 3, 2007 University of Florida (UF), Gainesville, Florida USA on CD-ROM, American Nuclear Society, LaGrange Park, IL (2007) ELECTRON DOSE

More information

Monitor Unit (MU) Calculation

Monitor Unit (MU) Calculation Monitor Unit (MU) Calculation Timothy C. Zhu 1, Haibo Lin 1, and JiaJian Shen 2 1 University of Pennsylvania, Philadelphia, PA 2 Mayo Clinic, Phoenix, AZ Introduction Pencil-beam based dose/mu algorithms

More information

MEDICAL IMAGING 2nd Part Computed Tomography

MEDICAL 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 information

Graphical User Interface for Simplified Neutron Transport Calculations

Graphical User Interface for Simplified Neutron Transport Calculations Graphical User Interface for Simplified Neutron Transport Calculations Phase 1 Final Report Instrument No: DE-SC0002321 July 20, 2009, through April 19, 2010 Recipient: Randolph Schwarz, Visual Editor

More information

Validation of GEANT4 for Accurate Modeling of 111 In SPECT Acquisition

Validation of GEANT4 for Accurate Modeling of 111 In SPECT Acquisition Validation of GEANT4 for Accurate Modeling of 111 In SPECT Acquisition Bernd Schweizer, Andreas Goedicke Philips Technology Research Laboratories, Aachen, Germany bernd.schweizer@philips.com Abstract.

More information

Surface Area and Volume

Surface Area and Volume Surface Area and Volume Level 1 2 1. Calculate the surface area and volume of each shape. Use metres for all lengths. Write your answers to 4 decimal places: a) 0.8 m Surface Area: Volume: b) 1 m 0.2 m

More information

Chapter 7: Geometrical Optics

Chapter 7: Geometrical Optics Chapter 7: Geometrical Optics 7. Reflection at a Spherical Surface L.O 7.. State laws of reflection Laws of reflection state: L.O The incident ray, the reflected ray and the normal all lie in the same

More information

Comparison of Predictions by MCNP and EGSnrc of Radiation Dose

Comparison of Predictions by MCNP and EGSnrc of Radiation Dose Comparison of Predictions by MCNP and EGSnrc of Radiation Dose Computational Medical Physics Working Group Workshop II, Sep 30 Oct 3, 2007 at UF imparted to various Material Targets by Beams and small

More information

Transmission Electron Microscopy 2. Scattering and Diffraction

Transmission Electron Microscopy 2. Scattering and Diffraction Transmission Electron Microscopy 2. Scattering and Diffraction EMA 6518 Spring 2007 01/07 Outline Why are we interested in electron scattering? Terminology of scattering The characteristics of electron

More information

REFLECTION & REFRACTION

REFLECTION & REFRACTION REFLECTION & REFRACTION OBJECTIVE: To study and verify the laws of reflection and refraction using a plane mirror and a glass block. To see the virtual images that can be formed by the reflection and refraction

More information

Identification of Shielding Material Configurations Using NMIS Imaging

Identification of Shielding Material Configurations Using NMIS Imaging Identification of Shielding Material Configurations Using NMIS Imaging B. R. Grogan, J. T. Mihalczo, S. M. McConchie, and J. A. Mullens Oak Ridge National Laboratory, P.O. Box 2008, MS-6010, Oak Ridge,

More information

Chapter 7: Geometrical Optics. The branch of physics which studies the properties of light using the ray model of light.

Chapter 7: Geometrical Optics. The branch of physics which studies the properties of light using the ray model of light. Chapter 7: Geometrical Optics The branch of physics which studies the properties of light using the ray model of light. Overview Geometrical Optics Spherical Mirror Refraction Thin Lens f u v r and f 2

More information

Large Plastic Scintillation Detectors for the Nuclear Materials Identification System

Large Plastic Scintillation Detectors for the Nuclear Materials Identification System Large Plastic Scintillation Detectors for the Nuclear Materials Identification System J.S. Neal, J.T. Mihalczo, M. T. Hiatt, J. D. Edwards Oak Ridge National Laboratory P. O. Box 2008, Oak Ridge, Tennessee

More information

Proton dose calculation algorithms and configuration data

Proton dose calculation algorithms and configuration data Proton dose calculation algorithms and configuration data Barbara Schaffner PTCOG 46 Educational workshop in Wanjie, 20. May 2007 VARIAN Medical Systems Agenda Broad beam algorithms Concept of pencil beam

More information

Optics Course (Phys 311) Geometrical Optics Refraction through Lenses

Optics Course (Phys 311) Geometrical Optics Refraction through Lenses Optics Course (Phys ) Geometrical Optics Refraction through Lenses Lecturer: Dr Zeina Hashim Slide 1 Objectives covered in this lesson : 1. Refraction through single spherical refracting surfaces. 2. Lenses:

More information

Dose Calculations: Where and How to Calculate Dose. Allen Holder Trinity University.

Dose Calculations: Where and How to Calculate Dose. Allen Holder Trinity University. Dose Calculations: Where and How to Calculate Dose Trinity University www.trinity.edu/aholder R. Acosta, W. Brick, A. Hanna, D. Lara, G. McQuilen, D. Nevin, P. Uhlig and B. Slater Dose Calculations - Why

More information

Deliverable D10.2. WP10 JRA04 INDESYS Innovative solutions for nuclear physics detectors

Deliverable D10.2. WP10 JRA04 INDESYS Innovative solutions for nuclear physics detectors MS116 Characterization of light production, propagation and collection for both organic and inorganic scintillators D10.2 R&D on new and existing scintillation materials: Report on the light production,

More information

2017 Summer Course on Optical Oceanography and Ocean Color Remote Sensing. Monte Carlo Simulation

2017 Summer Course on Optical Oceanography and Ocean Color Remote Sensing. Monte Carlo Simulation 2017 Summer Course on Optical Oceanography and Ocean Color Remote Sensing Curtis Mobley Monte Carlo Simulation Delivered at the Darling Marine Center, University of Maine July 2017 Copyright 2017 by Curtis

More information

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

Design and performance characteristics of a Cone Beam CT system for Leksell Gamma Knife Icon Design and performance characteristics of a Cone Beam CT system for Leksell Gamma Knife Icon WHITE PAPER Introduction Introducing an image guidance system based on Cone Beam CT (CBCT) and a mask immobilization

More information

Monte Carlo Method for Solving Inverse Problems of Radiation Transfer

Monte Carlo Method for Solving Inverse Problems of Radiation Transfer INVERSE AND ILL-POSED PROBLEMS SERIES Monte Carlo Method for Solving Inverse Problems of Radiation Transfer V.S.Antyufeev. ///VSP/// UTRECHT BOSTON KÖLN TOKYO 2000 Contents Chapter 1. Monte Carlo modifications

More information

Practical training Inga Kindermann; Dominik Haeussler

Practical training Inga Kindermann; Dominik Haeussler liht interference Materials: Laserpointer mountin material diffraction ratin with 40 slits per centimetre, CD-ROM (and ), cardboard security advise: Never let the direct beam of a laser hit your eyes!

More information

dq dt I = Irradiance or Light Intensity is Flux Φ per area A (W/m 2 ) Φ =

dq dt I = Irradiance or Light Intensity is Flux Φ per area A (W/m 2 ) Φ = Radiometry (From Intro to Optics, Pedrotti -4) Radiometry is measurement of Emag radiation (light) Consider a small spherical source Total energy radiating from the body over some time is Q total Radiant

More information

Digital Scatter Removal in Mammography to enable Patient Dose Reduction

Digital Scatter Removal in Mammography to enable Patient Dose Reduction Digital Scatter Removal in Mammography to enable Patient Dose Reduction Mary Cocker Radiation Physics and Protection Oxford University Hospitals NHS Trust Chris Tromans, Mike Brady University of Oxford

More information

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

Shadow casting. What is the problem? Cone Beam Computed Tomography THE OBJECTIVES OF DIAGNOSTIC IMAGING IDEAL DIAGNOSTIC IMAGING STUDY LIMITATIONS Cone Beam Computed Tomography THE OBJECTIVES OF DIAGNOSTIC IMAGING Reveal pathology Reveal the anatomic truth Steven R. Singer, DDS srs2@columbia.edu IDEAL DIAGNOSTIC IMAGING STUDY Provides desired diagnostic

More information

An Investigation of a Model of Percentage Depth Dose for Irregularly Shaped Fields

An Investigation of a Model of Percentage Depth Dose for Irregularly Shaped Fields Int. J. Cancer (Radiat. Oncol. Invest): 96, 140 145 (2001) 2001 Wiley-Liss, Inc. Publication of the International Union Against Cancer An Investigation of a Model of Percentage Depth Dose for Irregularly

More information

Optics Course (Phys 311) Geometrical Optics Refraction through Lenses

Optics Course (Phys 311) Geometrical Optics Refraction through Lenses Optics Course (Phys ) Geometrical Optics Refraction through Lenses Lecturer: Dr Zeina Hashim Slide 1 Objectives covered in this lesson : 1. The refracting power of a thin lens. 2. Thin lens combinations.

More information

2/26/2016. Chapter 23 Ray Optics. Chapter 23 Preview. Chapter 23 Preview

2/26/2016. Chapter 23 Ray Optics. Chapter 23 Preview. Chapter 23 Preview Chapter 23 Ray Optics Chapter Goal: To understand and apply the ray model of light. Slide 23-2 Chapter 23 Preview Slide 23-3 Chapter 23 Preview Slide 23-4 1 Chapter 23 Preview Slide 23-5 Chapter 23 Preview

More information

Chapter 36. Image Formation

Chapter 36. Image Formation Chapter 36 Image Formation Apr 22, 2012 Light from distant things We learn about a distant thing from the light it generates or redirects. The lenses in our eyes create images of objects our brains can

More information

2.5D Forward Solver To Model Scattering Of Long Dielectric Cylinders In An Active Millimeter Wave Imaging System

2.5D Forward Solver To Model Scattering Of Long Dielectric Cylinders In An Active Millimeter Wave Imaging System 2.5D Forward Solver To Model Scattering Of Long Dielectric Cylinders In An Active Millimeter Wave Imaging System S. Van den Bulcke and A. Franchois Department of Information Technology Ghent University,

More information

Mathematical methods and simulations tools useful in medical radiation physics

Mathematical methods and simulations tools useful in medical radiation physics Mathematical methods and simulations tools useful in medical radiation physics Michael Ljungberg, professor Department of Medical Radiation Physics Lund University SE-221 85 Lund, Sweden Major topic 1:

More information

Gamma spectroscopic measurements using the PID350 pixelated CdTe radiation detector

Gamma spectroscopic measurements using the PID350 pixelated CdTe radiation detector Gamma spectroscopic measurements using the PID350 pixelated CdTe radiation detector K. Karafasoulis, K. Zachariadou, S. Seferlis, I. Papadakis, D. Loukas, C. Lambropoulos, C. Potiriadis Abstract Spectroscopic

More information

Scatter in an uncollimated x-ray CT machine based on a Geant4 Monte Carlo simulation. Wadeson, Nicola and Morton, Edward and Lionheart, William

Scatter in an uncollimated x-ray CT machine based on a Geant4 Monte Carlo simulation. Wadeson, Nicola and Morton, Edward and Lionheart, William Scatter in an uncollimated x-ray CT machine based on a Geant4 Monte Carlo simulation Wadeson, Nicola and Morton, Edward and Lionheart, William 2010 MIMS EPrint: 2010.66 Manchester Institute for Mathematical

More information

DOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS

DOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS DOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS RAYLEIGH-SOMMERFELD DIFFRACTION RECTANGULAR APERTURES Ian Cooper School of Physics, University of Sydney ian.cooper@sydney.edu.au DOWNLOAD DIRECTORY FOR MATLAB

More information

Physics Midterm I

Physics Midterm I Phys121 - February 6, 2009 1 Physics 121 - Midterm I Last Name First Name Student Number Signature Tutorial T.A. (circle one): Ricky Chu Firuz Demir Maysam Emadi Alireza Jojjati Answer ALL 10 questions.

More information

A fast method for estimation of light flux in fluorescence image guided surgery

A fast method for estimation of light flux in fluorescence image guided surgery A fast method for estimation of light flux in fluorescence image guided surgery 1. Introduction In this document, we present a theoretical method to estimate the light flux in near-infrared fluorescence

More information

DETERMINISTIC 3D RADIATION TRANSPORT SIMULATION FOR DOSE DISTRIBUTION AND ORGAN DOSE EVALUATION IN DIAGNOSTIC CT

DETERMINISTIC 3D RADIATION TRANSPORT SIMULATION FOR DOSE DISTRIBUTION AND ORGAN DOSE EVALUATION IN DIAGNOSTIC CT DETERMINISTIC 3D RADIATION TRANSPORT SIMULATION FOR DOSE DISTRIBUTION AND ORGAN DOSE EVALUATION IN DIAGNOSTIC CT Monica Ghita,, Glenn Sjoden, Manuel Arreola, Ahmad Al-Basheer Basheer, Choonsik Lee, Wesley

More information

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

Impact of X-ray Scatter When Using CT-based Attenuation Correction in PET: A Monte Carlo Investigation 26 IEEE Nuclear Science Symposium Conference Record M6-349 Impact of X-ray Scatter When Using CT-based Attenuation Correction in PET: A Monte Carlo Investigation Habib Zaidi, Senior Member, IEEE and Mohammad

More information

Institute of Mechatronics and Information Systems

Institute of Mechatronics and Information Systems EXERCISE 4 Free vibrations of an electrical machine model Target Getting familiar with the fundamental issues of free vibrations analysis of a simplified model of an electrical machine, with the use of

More information

Algebra Based Physics

Algebra Based Physics Slide 1 / 66 Slide 2 / 66 Algebra Based Physics Geometric Optics 2015-12-01 www.njctl.org Table of ontents Slide 3 / 66 lick on the topic to go to that section Reflection Spherical Mirror Refraction and

More information

CHAPTER 10: TALLYING IN MCNP

CHAPTER 10: TALLYING IN MCNP _or_.e_sa_m_h_us_se_in 6:..:...7 ---"M.=o:.c.;;nte.:-C-"-=arlo Particle Transport with MCNP CHAPTER 10: TALLYING IN MCNP Tallying is the process of scoring the parameters of interest, Le. providing the

More information

GPU Based Convolution/Superposition Dose Calculation

GPU Based Convolution/Superposition Dose Calculation GPU Based Convolution/Superposition Dose Calculation Todd McNutt 1, Robert Jacques 1,2, Stefan Pencea 2, Sharon Dye 2, Michel Moreau 2 1) Johns Hopkins University 2) Elekta Research Funded by and Licensed

More information

Monte-Carlo modeling used to simulate propagation of photons in a medium

Monte-Carlo modeling used to simulate propagation of photons in a medium Monte-Carlo modeling used to simulate propagation of photons in a medium Nils Haëntjens Ocean Optics Class 2017 based on lectures from Emmanuel Boss and Edouard Leymarie What is Monte Carlo Modeling? Monte

More information

Modeling the ORTEC EX-100 Detector using MCNP

Modeling the ORTEC EX-100 Detector using MCNP Modeling the ORTEC EX-100 Detector using MCNP MCNP is a general-purpose Monte Carlo radiation transport code for modeling the interaction of radiation with materials based on composition and density. MCNP

More information

NDT OF SPECIMEN OF COMPLEX GEOMETRY USING ULTRASONIC ADAPTIVE

NDT OF SPECIMEN OF COMPLEX GEOMETRY USING ULTRASONIC ADAPTIVE NDT OF SPECIMEN OF COMPLEX GEOMETRY USING ULTRASONIC ADAPTIVE TECHNIQUES - THE F.A.U.S.T. SYSTEM INTRODUCTION O. Roy, S. Mahaut, M. Serre Commissariat a I'Energie Atomique CEAlCEREM, CE Saclay France Phased

More information

ACOMPTON-SCATTERING gamma camera records two

ACOMPTON-SCATTERING gamma camera records two IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 53, NO. 5, OCTOBER 2006 2787 Filtered Back-Projection in 4 Compton Imaging With a Single 3D Position Sensitive CdZnTe Detector Dan Xu, Student Member, IEEE, and

More information