Octree indexing of DICOM images for voxel number reduction and improvement of Monte Carlo simulation computing efficiency.

The purpose of the present study is to introduce a compression algorithm for the CT (computed tomography) data used in Monte Carlo simulations. Performing simulations on the CT data implies large computational costs as well as large memory requirements since the number of voxels in such data reaches typically into hundreds of millions voxels. CT data, however, contain homogeneous regions which could be regrouped to form larger voxels without affecting the simulation's accuracy. Based on this property we propose a compression algorithm based on octrees: in homogeneous regions the algorithm replaces groups of voxels with a smaller number of larger voxels. This reduces the number of voxels while keeping the critical high-density gradient area. Results obtained using the present algorithm on both phantom and clinical data show that compression rates up to 75% are possible without losing the dosimetric accuracy of the simulation.

[1]  Facundo Ballester,et al.  A dosimetric study of Leipzig applicators. , 2005, International journal of radiation oncology, biology, physics.

[2]  D W O Rogers,et al.  Monte Carlo calculation of nine megavoltage photon beam spectra using the BEAM code. , 2002, Medical physics.

[3]  Irene Gargantini,et al.  Linear octtrees for fast processing of three-dimensional objects , 1982, Comput. Graph. Image Process..

[4]  J. Baró,et al.  PENELOPE: An algorithm for Monte Carlo simulation of the penetration and energy loss of electrons and positrons in matter , 1995 .

[5]  I. Kawrakow,et al.  3D electron dose calculation using a Voxel based Monte Carlo algorithm (VMC). , 1996, Medical physics.

[6]  Hans Lundqvist,et al.  Monte Carlo calculations of thermal neutron capture in gadolinium: a comparison of GEANT4 and MCNP with measurements. , 2006, Medical physics.

[7]  J. Dempsey,et al.  Evaluation of the gamma dose distribution comparison method. , 2003, Medical physics.

[8]  E. Poon,et al.  Accuracy of the photon and electron physics in GEANT4 for radiotherapy applications. , 2005, Medical physics.

[9]  S. Nehmeh,et al.  Validation of GATE Monte Carlo simulations of the GE Advance/Discovery LS PET scanners. , 2005, Medical physics.

[10]  Raj Acharya,et al.  Adjacency algorithms for linear octree nodes , 1990, Image Vis. Comput..

[11]  Anatoly B Rosenfeld,et al.  The role of nonelastic reactions in absorbed dose distributions from therapeutic proton beams in different medium. , 2004, Medical physics.

[12]  Tianfang Li,et al.  Density resolution of proton computed tomography. , 2005, Medical physics.

[13]  Facundo Ballester,et al.  Monte Carlo calculation of the TG-43 dosimetric parameters of a new BEBIG Ir-192 HDR source. , 2005, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[14]  D. Low,et al.  A technique for the quantitative evaluation of dose distributions. , 1998, Medical physics.

[15]  C Lartizien,et al.  GATE: a simulation toolkit for PET and SPECT. , 2004, Physics in medicine and biology.

[16]  J. Seuntjens,et al.  Consistency test of the electron transport algorithm in the GEANT4 Monte Carlo code , 2005, Physics in medicine and biology.

[17]  Donald Meagher,et al.  Geometric modeling using octree encoding , 1982, Comput. Graph. Image Process..

[18]  Jean-François Carrier,et al.  Impact of interseed attenuation and tissue composition for permanent prostate implants. , 2006, Medical physics.

[19]  F Ballester,et al.  Monte Carlo dosimetric study of best industries and Alpha Omega Ir-192 brachytherapy seeds. , 2004, Medical physics.

[20]  F Ballester,et al.  Monte Carlo and experimental derivation of TG43 dosimetric parameters for CSM-type Cs-137 sources. , 2004, Medical physics.

[21]  F Ballester,et al.  Monte Carlo dosimetric study of the BEBIG Co-60 HDR source , 2005, Physics in medicine and biology.

[22]  Hans Lundqvist,et al.  Gadolinium neutron capture brachytherapy (GdNCB), a new treatment method for intravascular brachytherapy. , 2005, Medical physics.

[23]  H Paganetti,et al.  Adaptation of GEANT4 to Monte Carlo dose calculations based on CT data. , 2004, Medical physics.

[24]  H Paganetti,et al.  Test of GEANT3 and GEANT4 nuclear models for 160 MeV protons stopping in CH2. , 2003, Medical physics.

[25]  J. Cygler,et al.  Commissioning and quality assurance of treatment planning computers. , 1993, International journal of radiation oncology, biology, physics.

[26]  Brian Wang,et al.  Simulation of organ-specific patient effective dose due to secondary neutrons in proton radiation treatment , 2005, Physics in medicine and biology.

[27]  C. L. Le Rest,et al.  Validation of a Monte Carlo simulation of the Philips Allegro/GEMINI PET systems using GATE , 2006, Physics in medicine and biology.

[28]  H Paganetti,et al.  Accurate Monte Carlo simulations for nozzle design, commissioning and quality assurance for a proton radiation therapy facility. , 2004, Medical physics.

[29]  J. Sempau,et al.  DPM, a fast, accurate Monte Carlo code optimized for photon and electron radiotherapy treatment planning dose calculations , 2000 .

[30]  Igor Mishustin,et al.  Neutrons from fragmentation of light nuclei in tissue-like media: a study with the GEANT4 toolkit. , 2005, Physics in medicine and biology.

[31]  Jean Pouliot,et al.  Relative biological effectiveness enhancement of a 125I brachytherapy seed with characteristic x rays from its constitutive materials. , 2002, Medical physics.

[32]  Sargur N. Srihari,et al.  Representation of Three-Dimensional Digital Images , 1981, CSUR.

[33]  Antonio M Lallena,et al.  Dosimetry characterization of 32P intravascular brachytherapy source wires using Monte Carlo codes PENELOPE and GEANT4. , 2004, Medical physics.

[34]  G. Starkschall,et al.  American Association of Physicists in Medicine Radiation Therapy Committee Task Group 53: quality assurance for clinical radiotherapy treatment planning. , 1998, Medical physics.

[35]  L. Beaulieu,et al.  Validation of GEANT4, an object-oriented Monte Carlo toolkit, for simulations in medical physics. , 2004, Medical physics.