DUAL energy X-ray radiography [1] can be used to separate
|
|
- Barnard Scott
- 5 years ago
- Views:
Transcription
1 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 53, NO. 1, FEBRUARY A Scatter Correction Using Thickness Iteration in Dual-Energy Radiography S. K. Ahn, G. Cho, and H. Jeon Abstract In dual-energy radiography with area detectors, a scattered signal causes dominant error in the separation of different materials. Several methods for scatter correction in dual energy radiography have been suggested, yielding improved results. Such methods, however, require additional lead blocks or detectors, and additional exposures to estimate the scatter fraction for every correction. In the present study we suggest a scatter correction method that uses a database of fractions and distributions of the scattered radiation. To verify the feasibility of this method we conducted a MCNP simulation for a two-material problem, aluminum and water. The generation of the scatter information for different thicknesses of an aluminum-water phantom has been simulated. Based on the uncorrected signals, the thickness of each material can be calculated by a conventional dual-energy algorithm. The scatter information of the corresponding thickness from the database, a look-up table, is then used to correct the original signals. The iteration of this scatter correction reduced relative-thickness error from 32% to 3.4% in aluminum, and from 41% to 2.8% in water. The proposed scatter correction method can be applied to two-material dual-energy radiography such as mammography, contrast imaging, and industrial inspections. I. INTRODUCTION DUAL energy X-ray radiography [1] can be used to separate soft and dense-material images in medical and industrial applications. It can be performed successfully with a line-scanning system because of its scatter-free nature. With area detectors, however, scattered radiation also contributes to the signal. This undesired behavior of scattered X-ray photons in radiography causes serious degradation of the contrast in the observed images. This results in poor separation of soft and dense-material images in dual-energy radiography. As the scatter fraction increases, the thickness is underestimated. As such, a dual-energy algorithm does not work well. According to previous studies [2], the percentage of scattered photons is typically 60% to 70% in the lungs and 80% to 95% in the mediastinum in a standard PA chest radiograph. Anti-scatter grids can reduce these scatter fractions to 20% to 30% in the lungs and 40% to 60% in the mediastinum. Several methods [3] [6] for scatter correction in dual-energy radiography have been suggested, yielding improved results. These software post-processing methods are mostly based on convolution filtering. However, these approaches require additional lead blocks or detectors, and also additional exposures to estimate the scatter fraction for each correction. Furthermore, Manuscript received November 15, This work was supported by the itrs, Innovative Technology Center for Radiation Safety, Korea. The authors are with the Department of Nuclear and Quantum Engineering, Korea Advanced Institute of Science and Technology, Daejon, Korea ( skahn76@kaist.ac.kr; gscho@kaist.ac.kr; shiloh1@kaist.ac.kr). Digital Object Identifier /TNS Fig. 1. Schematic of the aluminum-water phantom simulated with MCNP and the resultant low- and high-energy images. The phantom has four steps of aluminum (0.1/0.2/0.3/0.4 cm) overlapped at the top and four steps of water (1/2/3/4 cm) at the bottom. The incident X-ray has a top-down parallel direction covering the whole top surface, and the detection pixel array is located immediately under the phantom. The low-energy spectrum is 70 kvp filtered by 10 mm Al, and the high-energy spectrum is 150 kvp filtered by 20 mm Cu. previous studies used only one or two convolution kernels even though the scatter point spread functions have different shapes as a function of object thickness. The present study suggests a new scatter correction method for use in dual-energy radiography. Using a dual-energy algorithm, we can roughly estimate the object thickness, which can assign the unique scatter-spread function for a more precise correction. In this study, we assumed that the object is composed of aluminum and water only, and that this composition is known. All resultant image and scatter point spread functions are simulated with MCNPX. Image separation by the dual-energy algorithm and its scatter correction method are described in Section II, and quantitative results obtained by simulations are discussed in detail in Section III. II. MATERIALS AND METHODS Fig. 1 provides a description of the aluminum-water phantom and the simulation geometry for the dual-energy radiography, and the resultant low- and high-energy images. The phantom has four steps of aluminum (0.1/0.2/0.3/0.4 cm) overlapped at the top and four steps of water (1/2/3/4 cm) at the bottom. The incident X-ray has a top-down parallel direction over the whole top surface, and the detection pixel array is located immediately under the phantom. The simulation tool was MCNPX code version [7], and photoelectric, Rayleigh and Compton scattering were taken into consideration. The detection array has pixels and the pixel size is cm. To obtain statistically confident results with a reasonable simulation time, /$ IEEE
2 134 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 53, NO. 1, FEBRUARY 2006 Fig. 2. Separated images of (left) the aluminum step and (right) the water step by a dual-energy algorithm without any scatter correction. Due to scatter effects, the thickness maps are seriously distorted and the water thickness is especially underestimated. the pixel size is much larger than that utilized in a real situation. There is no influence of pixel size on the correction. X-ray spectra were generated by a SRS78, a diagnostic X-ray spectra generator [8]. The low-energy spectrum is 70 kvp filtered by 10 mm Al, and the high-energy spectrum is 150 kvp filtered by 20 mm Cu. These simulation conditions, including the phantom, X-ray spectra and the detection conditions, are the same for all simulations in this work. The gray level in the images is not an absolute value, but rather it is equalized by maximum and minimum values in each image. For these images, we applied the dual-energy algorithm to separate the two materials, and the suggested scatter correction method to improve dual-energy imaging. The detailed procedure is described below. A. Dual-Energy Algorithm With monoenergetic X-ray sources, there is an exact solution of thickness because of perfect linearity between logarithmic signals and attenuation coefficients. With broad spectra sources, however, an exact solution does not exist because of the nonlinear relationship arising from broad energy and attenuation coefficients. Therefore, the thickness should be estimated by an empirical model, and this is the situation in general dual-energy radiography. This situation is worsened if there is scattered radiation in the signals. Using a dual-energy algorithm, we can separate aluminum and water in the images. Fig. 2 shows the separated thickness maps and the aluminum-only and the water-only images. In the thickness maps, there are serious distortions due to scatter effects, and the water thickness is especially underestimated. These aluminum and water thickness images are used as the seed-thickness in the scatter correction, which assigns a unique scatter point spread function from the look-up table. B. Generation of the SPSFs and the Look-Up Table The characteristics of the scatter point-spread functions (SPSFs) have been studied by several research groups. Their results from Monte Carlo simulations and experiments showed that the spatial distribution of the SPSF is rotationally symmetric [9]. Therefore, we have simulated a quarter of SPSFs and built whole SPSFs with the assumption they are symmetric. The quarter-spsfs in a dimension of are expanded into whole-spsfs in 47 47, and this dimension is diminished into by cutting out the near-zero section. Therefore, the final SPSFs have pixels with cm pixel size.
3 AHN et al.: SCATTER CORRECTION USING THICKNESS ITERATION 135 Fig. 4. Pixel value of the central pixel (primary signal), scatter signal, and their ratio SPR as a function of pixel size. We can see that the SPR is decreasing with larger pixel size (detection pixel is the circle with various radiuses) and the relationship between the SPR and the pixel-size depends on the object thickness. and various combinations of thickness of aluminum and water, comprise the look-up table (LUT) (Fig. 5). Fig. 3. (a) Schematic of the MCNP simulation for generation of scatter-spread functions and (b) some of results. We assumed the center pixel value is the primary intensity, and the others are scattered. In the simulation, the aluminum and water thicknesses are varied from 0.0 to 0.4 cm and 0 to 5 cm, respectively. A schematic of the MCNP simulation for generation of SPSFs is shown in Fig. 3. The X-ray beam is pencil-shaped. Only the center pixel has primary X-ray intensity, and the others have scattered X-rays. In the simulation, the aluminum thickness is varied at 0.0/0.1/0.2/0.3/0.4 cm, and the water thickness at 0/1/2/3/4/5 cm. Therefore, 30 different SPSFs are generated for low-energy and 30 SPSFs for high-energy. 5 cm-thick water has been simulated because the estimated thickness for the correction is thicker than the real thickness in the early part of iteration, even though the thickest water step is 4 cm in the example image. The bottom graph in Fig. 3 shows one set of results, displaying six different shapes of SPSF for low-energy with 0.4 cm fixed aluminum thickness and various water thicknesses. Due to the assumption of a primary signal from the pixel value of the central pixel and finite size of the pixel, the primary signal includes the scattered photons in a simulation. These effects are shown in Fig. 4. In the graph, the simulated primary signal (pixel value of the central pixel) and scatter signal, and their ratio (scatter-to-primary ratio), are given as a function of the pixel size. Several other simulations show that SPR is also a function of object thickness. Therefore, a correction for pixel size effects is needed. We corrected SPSFs with the absolute primary intensities, which are estimated by fitting the SPRs as a function of pixel size with 2-order exponential decay functions. The corrected final SPSFs, which are for low- and high-energies C. Generation of Scatter Images Using the LUT and the estimated thickness, low- and highenergy scatter-only images are generated. These scatter-only images are an integration of the pixel values in SPSFs on a pixel by pixel basis (Fig. 6). SPSFs are selected by seed-thickness of aluminum and water. The scatter-corrected images are obtained by subtraction of these scatter-only images from the original low- and high-energy images. D. Scatter Correction Algorithm The procedure of the suggested scatter correction method (Fig. 7) is as follows. 1. Estimation of object thickness using a dual-energy algorithm without scatter correction: These thickness images, which have large error, are seed-thicknesses for scatter images (described in Section II-A). 2. For each pixel, generation of scatter images by integration of scatter point spread functions from the look-up-table: Thickness of each pixel is used for the appropriate scatter point spread functions (described in Sections II-B and II-C). 3. Scatter correction of original low- and high-energy images: Scatter-corrected images are obtained by subtraction of scatter from original images (described in Section II-C). 4. Estimation of new object thickness with corrected lowand high-energy images: These estimated thicknesses give more correct results than the first seed-thickness. The result is used as the new seed-thickness for the iterative correction. 5. Comparison of scatter-only images: After generation of new scatter-only images with the results from step 4, the difference between the previous scatter-only and the new scatter-only images can be obtained.
4 136 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 53, NO. 1, FEBRUARY 2006 Fig. 5. LUT for the scatter correction. The LUT consisted of low- and high-energy scatter point spread functions of various combination of thickness of aluminum and water. Fig. 6. Algorithm for the generation of scatter-only images. Using the LUT and estimated thickness (seed-thickness), low- and high-energy scatter-only images are generated. 6. Iterate process 2 through 5 : Continue steps 2 through 5 with newly estimated thickness images as new seed-thickness until the difference becomes less than the pre-set error. This scatter correction process converges into some finite error of scatter signals if the SPR is less than 1 [10]. In this Fig. 7. Algorithm of the scatter correction method in dual-energy radiography. The dual-energy processing and the generation of scatter images were described in Sections II-A and II-C. work, the object is not thick enough for the SPR to be greater than 1, and thus the results converged in 4 iterations in the correction process. Fig. 8 shows the final aluminum- and water-thickness maps and images after the scatter correction.
5 AHN et al.: SCATTER CORRECTION USING THICKNESS ITERATION 137 Fig. 8. Separated images of (left) aluminum and (right) water after scatter correction. The resultant separated images show much clear steps in the images and better results in thickness estimation on the maps. III. RESULTS AND DISCUSSION Longitudinal and transversal profiles of the results are plotted in Fig. 9 and the relative-thickness errors are listed in Table I. In each graph, the real values and the results without scatter correction are plotted as dashed lines and open dots. The results with the scatter correction are much closer to the real values than those without any scatter correction. Without scatter correction, the average relative-thickness errors are 32% for aluminum and 41% for water, and the maxima are 44% and 75% in the longitudinal direction. The scatter correction reduces these errors to 3.4% for aluminum and 2.8% for water on average, and 11% and 9% as a maximum. In the zero-thickness regions, which are located right immediately next to the first steps (0.1 cm-step in the aluminum object and 1 cm-step in the water object), the errors are greater than the average because the scatter point spread functions have quite different shapes at the edges of the objects. It should be noted that these effects are especially emphasized in this study because of the steps in the phantom. In the transversal direction, the average relative-thickness errors decrease from 14% without correction to 6.1% with correction for 3 mm of aluminum, and 44% to 1.9% for 3 cm of water. The errors are greater for the thin region and the maxima reach 21% and 27% for 1 mm of Al and for 1 cm of water after correction. Detailed error comparisons are listed in Table I. The thickness is slightly more overestimated in the left-hand side region because the other material, which is stacked as presented in Fig. 1, is thicker in that region. IV. CONCLUSION A new scatter correction method in dual-energy radiography was suggested. We performed scatter correction by iterative thickness estimation using unique SPSFs. Our method uses information from a dual-energy algorithm to correct the images. This is the most distinct aspect of the proposed method in comparison with other conventional correction methods, which are applied in the same manner in conventional radiography and dual-energy radiography. In order to verify the effectiveness of this method, we conducted a MCNP simulation. The scatter information for each combination of thickness of an aluminum-water phantom was generated. Based on the uncorrected signals, the thickness of each material could be calculated by a conventional dual-energy algorithm. The scatter information of corresponding thickness from the database is then used to correct the original signals. For the aluminum-water object, the suggested iterative correction reduced the rms error in the thickness estimation. The method described in this paper does not use additional blocks and therefore does not need additional exposure.
6 138 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 53, NO. 1, FEBRUARY 2006 TABLE I RELATIVE ERRORS WITH AND WITHOUT SCATTER CORRECTION One drawback of the proposed approach arises from a characteristic of the SPSF. SPSFs with several different orders of object materials or with a different air gap thickness may have different shapes and magnitudes. In contrast with a CT, dual-energy radiography does not specify the arrangement of the material, i.e., which material is located on the top and which is on the bottom. Also, SPSFs are not symmetric near the object edge. This may constitute a handicap for the application of perfect correction in a blind object situation. Inclusion of these effects can help to decrease errors in regions of abrupt material thickness variations. It is noted that these effects were not taken into account in the present study, because preliminary simulation results indicated that the differences are insignificant for the range of interest in this study. Nevertheless, the correction method improved the dual-energy imaging considerably, even though it is not perfect. Most applications are not in a blind situation. For example, abnormal dense tissue is always located inside soft tissue in mammography. This scatter correction method can be applied to two-material dual-energy radiography such as mammography, contrast imaging, and industrial inspections. Expansion of the correction algorithm to other materials will be the next step in developing the proposed scatter correction method. Fig. 9. (a), (b) Longitudinal and (c), (d) transversal profiles of central and iso-thickness lines in the aluminum and the water images. The average relative thickness errors are 3.4% for aluminum and 2.8% for water in the longitudinal, and % for 1 3 mm of aluminum and % for 1 3 cm of water in the transversal direction. REFERENCES [1] L. A. Lehmann et al., Generalized image combinations in dual KVP digital radiography, Med. Phys., vol. 8, no. 5, pp , [2] C. E. Floyd et al., Scatter compensation for digital chest radiography using maximum likelihood expectation maximization, Invest. Radiol., vol. 28, no. 5, pp , [3] C.-G. Shaw et al., Quantitative digital subtraction angiography: Two scanning techniques for correction of scattered radiation and veiling glare, Radiology, vol. 157, pp , [4] F. C. Wagner et al., Dual-energy X-ray projection imaging: Two sampling schemes for the correction of scattered radiation, Med. Phys., vol. 15, no. 5, pp , [5] D. G. Kruger et al., A regional convolution kernel algorithm for scatter correction in dual-energy images: Comparison to single-kernel algorithms, Med. Phys., vol. 21, no. 2, pp , [6] D. A. Hinshaw and J. T. Dobbins III, Plate scatter correction for improved performance in dual-energy imaging, Med. Phys., vol. 23, no. 6, pp , [7] MCNPX 2.4.0, Oak Ridge Nat. Lab., Oak Ridge, TN, [8] K. Cranley, B. J. Gilmore, G. W. A. Fogarty, and L. Desponds, Catalogue of Diagnostic X-Ray Spectra and Other Data, Inst. Phys. Eng. Med. Rep. no. 78. [9] M. Honda et al., A technique of scatter-glare correction using a digital filtration, Med. Phys., vol. 20, no. 1, pp , [10] D. E. G. Trotter et al., Thickness-dependent scatter correction algorithm for digital mammography, in Proc. SPIE, 2002, pp
Detailed analysis of scatter contribution from different simulated geometries of X-ray detectors
Detailed analysis of scatter contribution from different simulated geometries of X-ray detectors Elena Marimon 1,2*, Hammadi Nait-Charif 1, Asmar Khan 2, Philip A. Marsden 3, Oliver Diaz 4 1 Centre for
More informationEnergy resolved X-ray diffraction Cl. J.Kosanetzky, G.Harding, U.Neitzel
Proc. Of SPIE Vol 0626, Application of Optical Instrumentation in Medicine XIV and Picture Archiving and Communication Systems(PACS IV) for Medical Applications IV, ed. S J Dwyer/R H Schneider (Jan 1986)
More informationEvaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT
Evaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT Qiao Yang 1,4, Meng Wu 2, Andreas Maier 1,3,4, Joachim Hornegger 1,3,4, Rebecca Fahrig
More informationScatter Correction Methods in Dimensional CT
Scatter Correction Methods in Dimensional CT Matthias Baer 1,2, Michael Hammer 3, Michael Knaup 1, Ingomar Schmidt 3, Ralf Christoph 3, Marc Kachelrieß 2 1 Institute of Medical Physics, Friedrich-Alexander-University
More informationArtifact 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 informationDigital 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 informationSpiral CT. Protocol Optimization & Quality Assurance. Ge Wang, Ph.D. Department of Radiology University of Iowa Iowa City, Iowa 52242, USA
Spiral CT Protocol Optimization & Quality Assurance Ge Wang, Ph.D. Department of Radiology University of Iowa Iowa City, Iowa 52242, USA Spiral CT Protocol Optimization & Quality Assurance Protocol optimization
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 informationApplication 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 informationProton 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 informationContrast Enhancement with Dual Energy CT for the Assessment of Atherosclerosis
Contrast Enhancement with Dual Energy CT for the Assessment of Atherosclerosis Stefan C. Saur 1, Hatem Alkadhi 2, Luca Regazzoni 1, Simon Eugster 1, Gábor Székely 1, Philippe Cattin 1,3 1 Computer Vision
More informationBeam Attenuation Grid Based Scatter Correction Algorithm for. Cone Beam Volume CT
11th European Conference on Non-Destructive Testing (ECNDT 2014), October 6-10, 2014, Prague, Czech Republic Beam Attenuation Grid Based Scatter Correction Algorithm for More Info at Open Access Database
More informationAN X-ray system with a large area detector, such as is used
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 25, NO. 12, DECEMBER 2006 1573 Scatter Correction Method for X-Ray CT Using Primary Modulation: Theory and Preliminary Results Lei Zhu*, Member, IEEE, N. Robert
More informationComparison of internal and external dose conversion factors using ICRP adult male and MEET Man voxel model phantoms.
Comparison of internal and external dose conversion factors using ICRP adult male and MEET Man voxel model phantoms. D.Leone, A.Häußler Intitute for Nuclear Waste Disposal, Karlsruhe Institute for Technology,
More informationREMOVAL OF THE EFFECT OF COMPTON SCATTERING IN 3-D WHOLE BODY POSITRON EMISSION TOMOGRAPHY BY MONTE CARLO
REMOVAL OF THE EFFECT OF COMPTON SCATTERING IN 3-D WHOLE BODY POSITRON EMISSION TOMOGRAPHY BY MONTE CARLO Abstract C.S. Levin, Y-C Tai, E.J. Hoffman, M. Dahlbom, T.H. Farquhar UCLA School of Medicine Division
More informationScatter 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 informationMetal Artifact Reduction CT Techniques. Tobias Dietrich University Hospital Balgrist University of Zurich Switzerland
Metal Artifact Reduction CT Techniques R S S S Tobias Dietrich University Hospital Balgrist University of Zurich Switzerland N. 1 v o 4 1 0 2. Postoperative CT Metal Implants CT is accurate for assessment
More informationELECTRON 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 informationScatter Correction for Dual source Cone beam CT Using the Pre patient Grid. Yingxuan Chen. Graduate Program in Medical Physics Duke University
Scatter Correction for Dual source Cone beam CT Using the Pre patient Grid by Yingxuan Chen Graduate Program in Medical Physics Duke University Date: Approved: Lei Ren, Supervisor Fang Fang Yin, Chair
More informationBasics 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 informationValidation 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 informationSimulation of Mammograms & Tomosynthesis imaging with Cone Beam Breast CT images
Simulation of Mammograms & Tomosynthesis imaging with Cone Beam Breast CT images Tao Han, Chris C. Shaw, Lingyun Chen, Chao-jen Lai, Xinming Liu, Tianpeng Wang Digital Imaging Research Laboratory (DIRL),
More informationCh. 4 Physical Principles of CT
Ch. 4 Physical Principles of CT CLRS 408: Intro to CT Department of Radiation Sciences Review: Why CT? Solution for radiography/tomography limitations Superimposition of structures Distinguishing between
More informationSimulation 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 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 informationScatter reduction for grid-less mammography using the convolutionbased image post-processing technique
Scatter reduction for grid-less mammography using the convolutionbased image post-processing technique Elena Marimón *a,b, Hammadi Nait-Charif a, Asmar Khan b, Philip A. Marsden c, Oliver Diaz d a CDE
More informationLow-Dose Dual-Energy CT for PET Attenuation Correction with Statistical Sinogram Restoration
Low-Dose Dual-Energy CT for PET Attenuation Correction with Statistical Sinogram Restoration Joonki Noh, Jeffrey A. Fessler EECS Department, The University of Michigan Paul E. Kinahan Radiology Department,
More informationSimulation of Beam Hardening in Industrial CT with x-ray and Monoenergetic Source by Monte Carlo Code
2012, TextRoad Publication ISSN 2090-4304 Journal of Basic and Applied Scientific Research www.textroad.com Simulation of Beam Hardening in Industrial CT with x-ray and Monoenergetic Source by Monte Carlo
More informationInvestigation of tilted dose kernels for portal dose prediction in a-si electronic portal imagers
Investigation of tilted dose kernels for portal dose prediction in a-si electronic portal imagers Krista Chytyk MSc student Supervisor: Dr. Boyd McCurdy Introduction The objective of cancer radiotherapy
More informationMonte 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 informationCIVA Computed Tomography Modeling
CIVA Computed Tomography Modeling R. FERNANDEZ, EXTENDE, France S. LEGOUPIL, M. COSTIN, D. TISSEUR, A. LEVEQUE, CEA-LIST, France page 1 Summary Context From CIVA RT to CIVA CT Reconstruction Methods Applications
More information8/7/2017. Disclosures. MECT Systems Overview and Quantitative Opportunities. Overview. Computed Tomography (CT) CT Numbers. Polyenergetic Acquisition
Quantitative Multi-Energy Computed Tomography: Imaging and Therapy Advancements Disclosures MECT Systems Overview and Quantitative Opportunities The speaker receives research funding from GE Healthcare
More informationIndustrial Radiography Simulation by MCNPX for Pipeline Corrosion Detection
More Info at Open Access Database www.ndt.net/?id=18672 Industrial Radiography Simulation by MCNPX for Pipeline Corrosion Detection Mostafa KHODABANDELOU 1, Amir MOVAFEGHI 2 *, Peyman ROSTAMI 1 and Jahanghir
More informationCBCT Equivalent Source Generation Using HVL and Beam Profile Measurements. Johnny Little PSM - Medical Physics Graduate Student University of Arizona
CBCT Equivalent Source Generation Using HVL and Beam Profile Measurements. Johnny Little PSM - Medical Physics Graduate Student University of Arizona Introduction CBCT has become a routine procedure for
More informationLimitations 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 informationEnhanced material contrast by dual-energy microct imaging
Enhanced material contrast by dual-energy microct imaging Method note Page 1 of 12 2 Method note: Dual-energy microct analysis 1. Introduction 1.1. The basis for dual energy imaging Micro-computed tomography
More informationISOCS 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 informationRadiographic Simulator artist: Version 2
18th World Conference on Nondestructive Testing, 16-20 April 2012, Durban, South Africa Radiographic Simulator artist: Version 2 Carsten Bellon 1, Andreas Deresch 1, Christian Gollwitzer 1, Gerd-Rüdiger
More informationCT vs. VolumeScope: image quality and dose comparison
CT vs. VolumeScope: image quality and dose comparison V.N. Vasiliev *a, A.F. Gamaliy **b, M.Yu. Zaytsev b, K.V. Zaytseva ***b a Russian Sci. Center of Roentgenology & Radiology, 86, Profsoyuznaya, Moscow,
More informationDUE to beam polychromacity in CT and the energy dependence
1 Empirical Water Precorrection for Cone-Beam Computed Tomography Katia Sourbelle, Marc Kachelrieß, Member, IEEE, and Willi A. Kalender Abstract We propose an algorithm to correct for the cupping artifact
More information3/27/2012 WHY SPECT / CT? SPECT / CT Basic Principles. Advantages of SPECT. Advantages of CT. Dr John C. Dickson, Principal Physicist UCLH
3/27/212 Advantages of SPECT SPECT / CT Basic Principles Dr John C. Dickson, Principal Physicist UCLH Institute of Nuclear Medicine, University College London Hospitals and University College London john.dickson@uclh.nhs.uk
More informationCombining Analytical and Monte Carlo Modelling for Industrial Radiology
19 th World Conference on Non-Destructive Testing 2016 Combining Analytical and Monte Carlo Modelling for Industrial Radiology Carsten BELLON, Gerd-Rüdiger JAENISCH, Andreas DERESCH BAM Bundesanstalt für
More informationEmpirical cupping correction: A first-order raw data precorrection for cone-beam computed tomography
Empirical cupping correction: A first-order raw data precorrection for cone-beam computed tomography Marc Kachelrieß, a Katia Sourbelle, and Willi A. Kalender Institute of Medical Physics, University of
More informationDetermination 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 informationMedical 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 informationOptimization of scanner parameters for dual energy micro-ct
Optimization of scanner parameters for dual energy micro-ct E. PAUWELS* 1, J. DHAENE 1, A. DE MUYNCK 1 E., M. DIERICK 1, L. VAN HOOREBEKE 1 1 UGCT Dept. Physics and Astronomy, Ghent University, Proeftuinstraat
More informationEffects of the difference in tube voltage of the CT scanner on. dose calculation
Effects of the difference in tube voltage of the CT scanner on dose calculation Dong Joo Rhee, Sung-woo Kim, Dong Hyeok Jeong Medical and Radiological Physics Laboratory, Dongnam Institute of Radiological
More informationSpectral 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 informationImpact 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 informationBackground. Outline. Radiographic Tomosynthesis: Image Quality and Artifacts Reduction 1 / GE /
Radiographic Tomosynthesis: Image Quality and Artifacts Reduction Baojun Li, Ph.D Department of Radiology Boston University Medical Center 2012 AAPM Annual Meeting Background Linear Trajectory Tomosynthesis
More informationUNCOMPROMISING QUALITY
ION CHAMBERS UNCOMPROMISING QUALITY Designed with over 30 years of scientific integrity for a broad range of dosimetry measurements in diverse radiation beams Farmer-type Chambers For absolute dosimetry
More informationPhase-Contrast Imaging and Tomography at 60 kev using a Conventional X-ray Tube
Phase-Contrast Imaging and Tomography at 60 kev using a Conventional X-ray Tube T. Donath* a, F. Pfeiffer a,b, O. Bunk a, W. Groot a, M. Bednarzik a, C. Grünzweig a, E. Hempel c, S. Popescu c, M. Hoheisel
More information8/2/2016. Measures the degradation/distortion of the acquired image (relative to an ideal image) using a quantitative figure-of-merit
Ke Li Assistant Professor Department of Medical Physics and Department of Radiology School of Medicine and Public Health, University of Wisconsin-Madison This work is partially supported by an NIH Grant
More informationA fast scatter field estimator for Digital Breast Tomosynthesis
A fast scatter field estimator for Digital Breast Tomosynthesis Oliver Díaz a, David R. Dance b,c, Kenneth C. Young b,c, Premkumar Elangovan a, Predrag R. Bakic d and Kevin Wells a a Centre for Vision,
More information1. Deployment of a framework for drawing a correspondence between simple figure of merits (FOM) and quantitative imaging performance in CT.
Progress report: Development of assessment and predictive metrics for quantitative imaging in chest CT Subaward No: HHSN6801000050C (4a) PI: Ehsan Samei Reporting Period: month 1-18 Deliverables: 1. Deployment
More informationarxiv: 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 information18th 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 informationDevelopments in Dimensional Metrology in X-ray Computed Tomography at NPL
Developments in Dimensional Metrology in X-ray Computed Tomography at NPL Wenjuan Sun and Stephen Brown 10 th May 2016 1 Possible factors influencing XCT measurements Components Influencing variables Possible
More informationShadow 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 informationDosimetry Simulations with the UF-B Series Phantoms using the PENTRAN-MP Code System
Dosimetry Simulations with the UF-B Series Phantoms using the PENTRAN-MP Code System A. Al-Basheer, M. Ghita, G. Sjoden, W. Bolch, C. Lee, and the ALRADS Group Computational Medical Physics Team Nuclear
More informationDUAL-ENERGY CT IN PROTON THERAPY
10/31/17 DUAL-ENERGY CT IN PROTON THERAPY Isabel Almeida, MAASTRO Clinic 7th NCS Lustrum Symposium 1 10/31/17 http://zonptc.bouwwebcam.nl https://www.youtube.com/watch?v=3vvvf5bqn7g Range uncertainties
More informationElectron 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 informationRefinement of Imaging Processing of Scatter Correction and Beam Hardening Tools for Industrial Radiography and Cone Beam CT
Digital Industrial Radiology and Computed Tomography (DIR 2015) 22-25 June 2015, Belgium, Ghent - www.ndt.net/app.dir2015 More Info at Open Access Database www.ndt.net/?id=18073 Refinement of Imaging Processing
More informationSuperior Imaging on the Go
Superior Imaging on the Go About Samsung Samsung Electronics Co., Ltd. is a global leader in technology, opening new possibilities for people everywhere. Through relentless innovation and discovery, we
More informationComputer-Aided Diagnosis in Abdominal and Cardiac Radiology Using Neural Networks
Computer-Aided Diagnosis in Abdominal and Cardiac Radiology Using Neural Networks Du-Yih Tsai, Masaru Sekiya and Yongbum Lee Department of Radiological Technology, School of Health Sciences, Faculty of
More information8/1/2017. Current Technology: Energy Integrating Detectors. Principles, Pitfalls and Progress in Photon-Counting-Detector Technology.
Photon Counting Detectors and Their Applications in Medical Imaging Principles, Pitfalls and Progress in Photon-Counting-Detector Technology Taly Gilat Schmidt, PhD Associate Professor Department of Biomedical
More informationModelling I SPECIFIC REQUIREMENTS ON COMPUTER MODELING OF RADIOGRAPHIC INSPECTIONS OF NUCLEAR COMPONENTS
Modelling I Inspection of Complex Geometries using Radiographic Simulation in CIVA A. Schumm, EDF R&D-SINETICS, France ; R. Fernandez, CEA/LIST, France ; J. Tabary, CEA/LETI, France ABSTRACT The application
More informationAcknowledgments. High Performance Cone-Beam CT of Acute Traumatic Brain Injury
A. Sisniega et al. (presented at RSNA 214) High Performance Cone-Beam CT of Acute Traumatic Brain Injury A. Sisniega 1 W. Zbijewski 1, H. Dang 1, J. Xu 1 J. W. Stayman 1, J. Yorkston 2, N. Aygun 3 V. Koliatsos
More informationAssesing 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 informationAn Iterative Approach to the Beam Hardening Correction in Cone Beam CT (Proceedings)
Marquette University e-publications@marquette Biomedical Engineering Faculty Research and Publications Engineering, College of 1-1-1999 An Iterative Approach to the Beam Hardening Correction in Cone Beam
More informationPerformance Evaluation of the TINA Medical Image Segmentation Algorithm on Brainweb Simulated Images
Tina Memo No. 2008-003 Internal Memo Performance Evaluation of the TINA Medical Image Segmentation Algorithm on Brainweb Simulated Images P. A. Bromiley Last updated 20 / 12 / 2007 Imaging Science and
More informationEvaluation of RayXpert for shielding design of medical facilities
Evaluation of Raypert for shielding design of medical facilities Sylvie Derreumaux 1,*, Sophie Vecchiola 1, Thomas Geoffray 2, and Cécile Etard 1 1 Institut for radiation protection and nuclear safety,
More informationMonte Carlo modelling and applications to imaging
Monte Carlo modelling and applications to imaging The Monte Carlo method is a method to obtain a result through repeated random sampling of the outcome of a system. One of the earliest applications, in
More informationWorkshop on Quantitative SPECT and PET Brain Studies January, 2013 PUCRS, Porto Alegre, Brasil Corrections in SPECT and PET
Workshop on Quantitative SPECT and PET Brain Studies 14-16 January, 2013 PUCRS, Porto Alegre, Brasil Corrections in SPECT and PET Físico João Alfredo Borges, Me. Corrections in SPECT and PET SPECT and
More informationPhoton counting spectroscopic CT with dynamic beam attenuator
Photon counting spectroscopic CT with dynamic beam attenuator Haluk Atak 1) and Polad M. Shikhaliev a) 1) Department of Nuclear Engineering, Hacettepe University, Ankara, Turkey Abstract Purpose: Photon
More informationMonte Carlo simulation of scatter field for calculation of contrast of discs in synthetic CDMAM images
Monte Carlo simulation of scatter field for calculation of contrast of discs in synthetic CDMAM images Oliver Díaz 1, Mary Yip 1, Jorge Cabello 1, David R. Dance 2,3, Kenneth C. Young 2,3, and Kevin Wells
More informationACCURATE 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 informationMUSICA makes the difference
MUSICA makes the difference At busy Logan Hospital, MUSICA and the DR 600 are keeping the imaging workflow smooth and efficient, while meeting the quality and dose reduction needs for patients and staff
More informationA Radiometry Tolerant Method for Direct 3D/2D Registration of Computed Tomography Data to X-ray Images
A Radiometry Tolerant Method for Direct 3D/2D Registration of Computed Tomography Data to X-ray Images Transfer Function Independent Registration Boris Peter Selby 1, Georgios Sakas 2, Stefan Walter 1,
More informationDeep 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 information8/3/2016. Outline. The EPID Strikes Back: Future EPID Technology and Applications. Active Matrix Flat-Panel Imagers (AMFPIs)
8//6 The EPID Strikes Back: Future EPID Technology and Applications Larry E. Antonuk Department of Radiation Oncology University of Michigan, Ann Arbor Acknowledgements: Youcef El-Mohri, Qihua Zhao (U.
More informationA Software Tool for Measurement of the Modulation Transfer Function
A Software Tool for Measurement of the Modulation Transfer Function Yen-Hong Kao a, Michael Albert a, Ann-Katherine Carton a, Hilde Bosmans b, Andrew D. A. Maidment a a Department of Radiology, University
More informationImage-based Monte Carlo calculations for dosimetry
Image-based Monte Carlo calculations for dosimetry Irène Buvat Imagerie et Modélisation en Neurobiologie et Cancérologie UMR 8165 CNRS Universités Paris 7 et Paris 11 Orsay, France buvat@imnc.in2p3.fr
More informationDisclosure 7/24/2014. Validation of Monte Carlo Simulations For Medical Imaging Experimental validation and the AAPM Task Group 195 Report
Validation of Monte Carlo Simulations For Medical Imaging Experimental validation and the AAPM Task Group 195 Report Ioannis Sechopoulos, Ph.D., DABR Diagnostic Imaging Physics Lab Department of Radiology
More informationJapan Foundry Society, Inc. Application of Recent X-ray CT Technology to Investment Casting field. Kouichi Inagaki ICC / IHI Corporation
Japan Foundry Society, Inc. Application of Recent X-ray CT Technology to Investment Casting field Kouichi Inagaki ICC / IHI Corporation 13 th WORLD CONFERENCE ON INVESTMENT CASTING Paper: T3 Copyright
More informationConstruction of Voxel-type Phantom Based on Computed Tomographic Data of RANDO Phantom for the Monte Carlo Simulations
Construction of Voxel-type Phantom Based on Computed Tomographic Data of RANDO Phantom for the Monte Carlo Simulations K. Minami 1, K. Ejiri 1, M. Shimo 1, M. Kato, Y. Takeuchi, K. Yonemochi, H. Toyama
More informationCone-beam computed tomography with a flat-panel imager: Magnitude and effects of x-ray scatter
Cone-beam computed tomography with a flat-panel imager: Magnitude and effects of x-ray scatter Jeffrey H. Siewerdsen a) and David A. Jaffray Department of Radiation Oncology, William Beaumont Hospital,
More informationx' = c 1 x + c 2 y + c 3 xy + c 4 y' = c 5 x + c 6 y + c 7 xy + c 8
1. Explain about gray level interpolation. The distortion correction equations yield non integer values for x' and y'. Because the distorted image g is digital, its pixel values are defined only at integer
More informationMonte Carlo methods in proton beam radiation therapy. Harald Paganetti
Monte Carlo methods in proton beam radiation therapy Harald Paganetti Introduction: Proton Physics Electromagnetic energy loss of protons Distal distribution Dose [%] 120 100 80 60 40 p e p Ionization
More informationC E N T E R A T H O U S T O N S C H O O L of H E A L T H I N F O R M A T I O N S C I E N C E S. Image Operations I
T H E U N I V E R S I T Y of T E X A S H E A L T H S C I E N C E C E N T E R A T H O U S T O N S C H O O L of H E A L T H I N F O R M A T I O N S C I E N C E S Image Operations I For students of HI 5323
More informationABSTRACT 1. INTRODUCTION
Estimating breast thickness for dual-energy subtraction in contrastenhanced digital mammography using calibration phantoms Kristen C. Lau 1, Young Joon Kwon, Moez Karim Aziz, Raymond J. Acciavatti, Andrew
More informationHidenobu Tachibana The Cancer Institute Hospital of JFCR, Radiology Dept. The Cancer Institute of JFCR, Physics Dept.
2-D D Dose-CT Mapping in Geant4 Hidenobu Tachibana The Cancer Institute Hospital of JFCR, Radiology Dept. The Cancer Institute of JFCR, Physics Dept. Table of Contents Background & Purpose Materials Methods
More informationLarge 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 informationPartial Volume Simulation in Software Breast Phantoms
Partial Volume Simulation in Software Breast Phantoms Feiyu Chen, David Pokrajac, Xiquan Shi, Fengshan Liu, Andrew D.A. Maidment a, Predrag R. Bakic a Delaware State University, 1200 N DuPont Hwy, Dover
More informationOptimization of Beam Spectrum and Dose for Lower-Cost CT. Mary Esther Braswell. Graduate Program in Medical Physics Duke University.
Optimization of Beam Spectrum and Dose for Lower-Cost CT by Mary Esther Braswell Graduate Program in Medical Physics Duke University Date: Approved: James Dobbins, Supervisor Anuj Kapadia Robert Reiman
More informationComparison 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 informationCardiac Dual Energy CT: Technique
RSNA 2013, VSCA51-01, Chicago, Dec. 5, 2013 Cardiac Radiology Series Cardiac Dual Energy CT: Technique Willi A. Kalender, Ph.D. Institute of Medical Physics University of Erlangen www.imp.uni-erlangen.de
More informationRegistration concepts for the just-in-time artefact correction by means of virtual computed tomography
DIR 2007 - International Symposium on Digital industrial Radiology and Computed Tomography, June 25-27, 2007, Lyon, France Registration concepts for the just-in-time artefact correction by means of virtual
More informationSimulation of Internal Backscatter Effects on MTF and SNR of Pixelated Photon-counting Detectors
Simulation of Internal Backscatter Effects on MTF and SNR of Pixelated Photon-counting Detectors Alexander Korn, Juergen Giersch a and Martin Hoheisel b a Physikalisches Institut Universitaet Erlangen-Nuernberg,
More informationA study of densitometry comparison among three radiographic processing solutions
Iran. J. Radiat. Res., 2006; 4 (2): 81-86 A study of densitometry comparison among three radiographic processing solutions V. Changizi 1*, E. Jazayeri 1,A.Talaeepour 2 1 Department of Radiology Technology,
More information