Development of whole-body representation and dose calculation in a commercial treatment planning system

IS2aR, a computational tool to transform voxelized reference phantoms into patient-specific whole-body virtual CTs for peripheral dose estimation

BACKGROUND

The risk of radiogenic cancer induction due to radiotherapy depends on the dose received by the patient's organs. Knowing the position of all organs is needed to assess this dose in a personalized way. However, radiotherapy planning computed tomography (pCT) scans contain truncated patient anatomy, limiting personalized dose evaluation.

PURPOSE

To develop a simple and freely available computational tool that adapts an ICRP reference computational phantom to generate a patient-specific whole-body CT for peripheral dose computations.

METHODS

Various bone-segmentation methods were explored onto fifteen pCTs, and the one with the highest Sørensen-Dice coefficient was implemented. The reference phantom is registered to the pCT, obtaining a registration transform matrix, which is then applied to create the whole-body virtual CT. Additional validation involved a comparison of absorbed doses to organs delineated on both the pCT and the virtual CT.

RESULTS

A dedicated graphical user interface was designed and implemented to house the developed functions for i) selecting a registration region on which automatic bone segmentation and rigid registration will occur, ii) displaying the results of these processes under selectable views, and iii) exporting the final patient-specific whole-body CT. This software was termed IS2aR. The tested whole-body virtual CT generated by IS2aR fulfilled our requirements.

CONCLUSIONS

IS2aR is a user-friendly computational software to create a personalized whole-body CT containing the original structures in the reference phantom. The personalized dose deposited in peripheral organs can be estimated further to assess second cancer induction risk in epidemiological studies.

Experimental Validation of an Analytical Program and a Monte Carlo Simulation for the Computation of the Far Out-of-Field Dose in External Beam Photon Therapy Applied to Pediatric Patients

Background

The out-of-the-field absorbed dose affects the probability of primary second radiation-induced cancers. This is particularly relevant in the case of pediatric treatments. There are currently no methods employed in the clinical routine for the computation of dose distributions from stray radiation in radiotherapy. To overcome this limitation in the framework of conventional teletherapy with photon beams, two computational tools have been developed—one based on an analytical approach and another depending on a fast Monte Carlo algorithm. The purpose of this work is to evaluate the accuracy of these approaches by comparison with experimental data obtained from anthropomorphic phantom irradiations.

Materials and Methods

An anthropomorphic phantom representing a 5-year-old child (ATOM, CIRS) was irradiated considering a brain tumor using a Varian TrueBeam linac. Two treatments for the same planned target volume (PTV) were considered, namely, intensity-modulated radiotherapy (IMRT) and volumetric modulated arc therapy (VMAT). In all cases, the irradiation was conducted with a 6-MV energy beam using the flattening filter for a prescribed dose of 3.6 Gy to the PTV. The phantom had natLiF : Mg, Cu, P (MCP-N) thermoluminescent dosimeters (TLDs) in its 180 holes. The uncertainty of the experimental data was around 20%, which was mostly attributed to the MCP-N energy dependence. To calculate the out-of-field dose, an analytical algorithm was implemented to be run from a Varian Eclipse TPS. This algorithm considers that all anatomical structures are filled with water, with the exception of the lungs which are made of air. The fast Monte Carlo code dose planning method was also used for computing the out-of-field dose. It was executed from the dose verification system PRIMO using a phase-space file containing 3x109 histories, reaching an average standard statistical uncertainty of less than 0.2% (coverage factor k = 1 ) on all voxels scoring more than 50% of the maximum dose. The standard statistical uncertainty of out-of-field voxels in the Monte Carlo simulation did not exceed 5%. For the Monte Carlo simulation the actual chemical composition of the materials used in ATOM, as provided by the manufacturer, was employed.

Results

In the out-of-the-field region, the absorbed dose was on average four orders of magnitude lower than the dose at the PTV. For the two modalities employed, the discrepancy between the central values of the TLDs located in the out-of-the-field region and the corresponding positions in the analytic model were in general less than 40%. The discrepancy in the lung doses was more pronounced for IMRT. The same comparison between the experimental and the Monte Carlo data yielded differences which are, in general, smaller than 20%. It was observed that the VMAT irradiation produces the smallest out-of-the-field dose when compared to IMRT.

Conclusions

The proposed computational methods for the routine calculation of the out-of-the-field dose produce results that are similar, in most cases, with the experimental data. It has been experimentally found that the VMAT irradiation produces the smallest out-of-the-field dose when compared to IMRT for a given PTV.