Simulation study on potential accuracy gains from dual energy CT tissue segmentation for low-energy brachytherapy Monte Carlo dose calculations

This work compares Monte Carlo (MC) dose calculations for (125)I and (103)Pd low-dose rate (LDR) brachytherapy sources performed in virtual phantoms containing a series of human soft tissues of interest for brachytherapy. The geometries are segmented (tissue type and density assignment) based on simulated single energy computed tomography (SECT) and dual energy (DECT) images, as well as the all-water TG-43 approach. Accuracy is evaluated by comparison to a reference MC dose calculation performed in the same phantoms, where each voxel's material properties are assigned with exactly known values. The objective is to assess potential dose calculation accuracy gains from DECT. A CT imaging simulation package, ImaSim, is used to generate CT images of calibration and dose calculation phantoms at 80, 120, and 140 kVp. From the high and low energy images electron density ρ(e) and atomic number Z are obtained using a DECT algorithm. Following a correction derived from scans of the calibration phantom, accuracy on Z and ρ(e) of ±1% is obtained for all soft tissues with atomic number Z ∊ [6,8] except lung. GEANT4 MC dose calculations based on DECT segmentation agreed with the reference within ±4% for (103)Pd, the most sensitive source to tissue misassignments. SECT segmentation with three tissue bins as well as the TG-43 approach showed inferior accuracy with errors of up to 20%. Using seven tissue bins in our SECT segmentation brought errors within ±10% for (103)Pd. In general (125)I dose calculations showed higher accuracy than (103)Pd. Simulated image noise was found to decrease DECT accuracy by 3-4%. Our findings suggest that DECT-based segmentation yields improved accuracy when compared to SECT segmentation with seven tissue bins in LDR brachytherapy dose calculation for the specific case of our non-anthropomorphic phantom. The validity of our conclusions for clinical geometry as well as the importance of image noise in the tissue segmentation procedure deserves further experimental investigation.

[1]  F Verhaegen,et al.  SpekCalc: a program to calculate photon spectra from tungsten anode x-ray tubes , 2009, Physics in medicine and biology.

[2]  D E Raeside,et al.  The accuracy of single-seed dose superposition for I-125 implants. , 1989, Medical physics.

[3]  L. Anderson,et al.  Dosimetry of interstitial brachytherapy sources: Recommendations of the AAPM Radiation Therapy Committee Task Group No. 43 , 1995 .

[4]  J. Williamson,et al.  Update of AAPM Task Group No. 43 Report: A revised AAPM protocol for brachytherapy dose calculations. , 2004 .

[5]  Frank Verhaegen,et al.  Sensitivity of low energy brachytherapy Monte Carlo dose calculations to uncertainties in human tissue composition. , 2010, Medical physics.

[6]  E. Pedroni,et al.  The calibration of CT Hounsfield units for radiotherapy treatment planning. , 1996, Physics in medicine and biology.

[7]  Mitchell M Goodsitt,et al.  Accuracies of the synthesized monochromatic CT numbers and effective atomic numbers obtained with a rapid kVp switching dual energy CT scanner. , 2011, Medical physics.

[8]  R M Thomson,et al.  Monte Carlo dosimetry for I125 and Pd103 eye plaque brachytherapy. , 2008, Medical physics.

[9]  T. Bortfeld,et al.  Correlation between CT numbers and tissue parameters needed for Monte Carlo simulations of clinical dose distributions. , 2000, Physics in medicine and biology.

[10]  Frank Verhaegen,et al.  The difference of scoring dose to water or tissues in Monte Carlo dose calculations for low energy brachytherapy photon sources. , 2011, Medical physics.

[11]  Andreas H Mahnken,et al.  Spectral rhoZ-projection method for characterization of body fluids in computed tomography: ex vivo experiments. , 2009, Academic radiology.

[12]  Frank Verhaegen,et al.  Postimplant dosimetry using a Monte Carlo dose calculation engine: a new clinical standard. , 2007, International journal of radiation oncology, biology, physics.

[13]  Jean-François Carrier,et al.  SU‐FF‐T‐406: Toward a More Accurate Dose Calculation Technique Using a Semiautomatic Organ Contouring in Monte Carlo Post‐Implant Assessment of Breast LDR Brachytherapy , 2009 .

[14]  James A. Scott.  Photon, Electron, Proton and Neutron Interaction Data for Body Tissues ICRU Report 46. International Commission on Radiation Units and Measurements, Bethesda, 1992, $40.00 , 1993 .

[15]  K. Stierstorfer,et al.  Density and atomic number measurements with spectral x-ray attenuation method , 2003 .

[16]  D. R. White,et al.  The composition of body tissues. , 1986, The British journal of radiology.

[17]  F Verhaegen,et al.  Monte Carlo study of LDR seed dosimetry with an application in a clinical brachytherapy breast implant. , 2009, Medical physics.

[18]  R. Mohan,et al.  Theoretical variance analysis of single- and dual-energy computed tomography methods for calculating proton stopping power ratios of biological tissues , 2010, Physics in medicine and biology.

[19]  Philip M Evans,et al.  Calculation of x-ray spectra emerging from an x-ray tube. Part I. electron penetration characteristics in x-ray targets. , 2007, Medical physics.

[20]  J. Valentin.  Basic anatomical and physiological data for use in radiological protection: reference values , 2002, Annals of the ICRP.

[21]  Xiaochuan Pan,et al.  Optimal noise control in and fast reconstruction of fan-beam computed tomography image. , 1999, Medical physics.

[22]  Frank Verhaegen,et al.  Tissue segmentation in Monte Carlo treatment planning: a simulation study using dual-energy CT images. , 2008, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[23]  Frank Verhaegen,et al.  MO‐D‐303A‐06: ImaSim, An Animated Tool for Teaching Imaging , 2009 .

[24]  I. Kawrakow,et al.  History by history statistical estimators in the BEAM code system. , 2002, Medical physics.

[25]  F. Verhaegen,et al.  Dual-energy CT-based material extraction for tissue segmentation in Monte Carlo dose calculations , 2008, Physics in medicine and biology.

[26]  A. Boyer,et al.  Tissue Inhomogeneity Corrections for Megavoltage Photon Beams , 2004 .

[27]  M. Endo,et al.  Electron density measurement with dual-energy x-ray CT using synchrotron radiation. , 2003, Physics in medicine and biology.

[28]  B. Heismann,et al.  Quantitative image-based spectral reconstruction for computed tomography. , 2009, Medical physics.

[29]  F. Verhaegen,et al.  Sensitivity study for CT image use in Monte Carlo treatment planning , 2005, Physics in medicine and biology.

[30]  Omar Chibani,et al.  MCPI©: A sub-minute Monte Carlo dose calculation engine for prostate implants. , 2005, Medical physics.

[31]  B. Logan,et al.  Reconstructing Interior Head Tissue from X-Ray Transmissions , 1974 .

[32]  G. Poludniowski.  Calculation of x-ray spectra emerging from an x-ray tube. Part II. X-ray production and filtration in x-ray targets. , 2007, Medical physics.

[33]  H. E. Johns,et al.  Physics of Radiology , 1983 .

[34]  B. Whiting,et al.  On two-parameter models of photon cross sections: application to dual-energy CT imaging. , 2006, Medical physics.

[35]  Frank Verhaegen,et al.  Influence of breast composition and interseed attenuation in dose calculations for post-implant assessment of permanent breast 103Pd seed implant. , 2010, Physics in medicine and biology.