A comparison study between single- and dual-energy CT density extraction methods for neurological proton monte carlo treatment planning

Abstract

Monte Carlo proton dose calculations requires mass densities calculated from the patient CT image. This work investigates the impact of different single-energy CT (SECT) and dual-energy CT (DECT) to density conversion methods in proton dose distributions for brain tumours. Material and methods: Head CT scans for four patients were acquired in SECT and DECT acquisition modes. Commercial software was used to reconstruct DirectDensity™ images in Relative Electron Densities (RED, (Formula presented.)) and to obtain DECT-based pseudo-monoenergetic images (PMI). PMI and SECT images were converted to RED using piecewise linear interpolations calibrated on a head-sized phantom, these fits were referred to as “PMI2RED” and “CT2RED”. Two DECT-based calibration methods (“Hünemohr-15it” and “Saito-15it”) were also investigated. (Formula presented.) images were converted to mass-densities ((Formula presented.)) to investigate (Formula presented.) differences and one representative patient case was used to make a proton treatment plan. Using CT2RED as reference method, dose distribution differences in the target and in five organs-at-risk (OARs) were quantified. Results: In the phantom study, Saito-15it and Hünemohr-15it produced the lowest (Formula presented.) root-mean-square error (0.7%) and DirectDensity™ the highest error (2.7%). The proton plan evaluated in the Saito-15it and Hünemohr-15it datasets showed the largest relative differences compared to initial CT2RED plan down to −6% of the prescribed dose. Compared to CT2RED, average range differences were calculated: −0.1 ± 0.3 mm for PMI2RED; −0.8 ± 0.4 mm for Hünemohr-15it, and −1.2 ± 0.4 mm for Saito-15it. Conclusion: Given the wide choice of available conversion methods, studies investigating the density accuracy for proton dose calculations are necessary. However, there is still a gap between performing accuracy studies in reference (Formula presented.) phantoms and applying these methods in human CT images. For this treatment case, the PMI2RED method was equivalent to the conventional CT2RED method in terms of dose distribution, CTV coverage and OAR sparing, whereas Hünemohr-15it and Saito-15it presented the largest differences.