Validation of a method for “dose of the day” calculation in head-neck tomotherapy by using planning ct-to-MVCT deformable image registration

PurposeThe aim of this study was to test the feasibility and dosimetric accuracy of a method that employs planning CT-to-MVCT deformable image registration (DIR) for calculation of the daily dose for head and neck (HN) patients treated with Helical Tomotherapy (HT).MethodsFor each patient, the planning kVCT (CTplan) was deformably registered to the MVCT acquired at the 15th therapy session (MV15) with a B-Spline Free Form algorithm using Mattes mutual information (open-source software 3D Slicer), resulting in a deformed CT (CTdef). On the same day as MVCT15, a kVCT was acquired with the patient in the same treatment position (CT15). The original HT plans were recalculated both on CTdef and CT15, and the corresponding dose distributions were compared; local dose differences <2% of the prescribed dose (DD2%) and 2D/3D gamma-index values (2%-2 mm) were assessed respectively with Mapcheck SNC Patient software (Sun Nuclear) and with 3D-Slicer.ResultsOn average, 87.9% ± 1.2% of voxels were found for DD2% (on average 27 slices available for each patient) and 94.6% ± 0.8% of points passed the 2D gamma analysis test while the 3D gamma test was satisfied in 94.8% ± 0.8% of body’s voxels.ConclusionsThis study represents the first demonstration of the dosimetric accuracy of kVCT-to-MVCT DIR for dose of the day computations. The suggested method is sufficiently fast and reliable to be used for daily delivered dose evaluations in clinical strategies for adaptive Tomotherapy of HN cancer.     

Evaluation criteria for film based intensity modulated radiation therapy quality assurance

The aim of this study was to use different gamma histogram criteria for the comparison of planned dose with irradiated dose distribution and find that what percent of pixels passing a certain criteria imitate a good quality plan. The dose was calculated for 156 patients by inverse planning optimization using the Corvus treatment planning system. Gafchromic films in combination with 2571 0.6 cm³ Farmer type ionization chamber and Farmer 2570/1 electrometer from NE Technology were used to measure the delivered dose in solid water phantom. All the measurements were performed on Varian CL21EX linear accelerator (Varian Medical Systems, Palo Alto, CA) fitted with a Millennium 120 leaf collimator. In this study the mean value of the percent of passing pixels within the region of interest under the criterion of 3% DD and 3 mm DTA is 90.2 ± 7.1% for head and neck cases and 92.2 ± 5.8% for non-head and neck cases. If we choose the criteria of 3% DD and 3 mm DTA then 96.3% head and neck plans have the percent of passing pixels ≥ 75% and 95.1% non-head and neck plans have the percent of passing pixels ≥ 80%. It is evident from the results of this study that the criterion of 5% DD and 3 mm DTA with the percent of passing pixels ≥ 90 for non-head and neck cases while the percent of passing pixels ≥ 85 for head and neck cases endorse that a plan is good. The results of this study may be useful for other institutions which use verification software and EBT films for patient specific IMRT QA.     

A patient immobilization device for prone breast radiotherapy: Dosimetric effects and inclusion in the treatment planning system

PurposeTo assess the dosimetric impact of a patient positioning device for prone breast radiotherapy and assess the accuracy of a treatment planning system (TPS) in predicting this impact.MethodsBeam attenuation and build-up dose perturbations, quantified by ionization chamber and radiochromic film dosimetry, were evaluated for 3 components of the patient positioning device: the carbon fiber baseplate, the support cushions and the support wedge for the contralateral breast. Dose calculations were performed using the XVMC dose engine implemented in the Monaco TPS. All components were included during planning CT acquisition.ResultsBeam attenuation amounted to 7.57% (6 MV) and 5.33% (15 MV) for beams obliquely intersecting the couchtop–baseplate combination. Beams traversing large sections of the support wedge were attenuated by 12.28% (6 MV) and 9.37% (15 MV). For the support cushion foam, beam attenuation remained limited to 0.11% (6 MV) and 0.08% (15 MV) per centimeter thickness. A substantial loss of dose build-up was detected when irradiating through any of the investigated components. TPS dose calculations accurately predicted beam attenuation by the baseplate and support wedge. A manual density overwrite was needed to model attenuation by the support cushion foam. TPS dose calculations in build-up regions differed considerably from measurements for both open beams and beams traversing the device components.ConclusionsIrradiating through the components of the positioning device resulted in a considerable degradation of skin sparing. Inclusion of the device components in the treatment planning CT allowed to accurately model the most important attenuation effect, but failed to accurately predict build-up doses.     

3D-printed applicators for high dose rate brachytherapy: Dosimetric assessment at different infill percentage

PurposeDosimetric assessment of high dose rate (HDR) brachytherapy applicators, printed in 3D with acrylonitrile butadiene styrene (ABS) at different infill percentage.Materials and methodsA low-cost, desktop, 3D printer (Hamlet 3DX100, Hamlet, Dublin, IE) was used for manufacturing simple HDR applicators, reproducing typical geometries in brachytherapy: cylindrical (common in vaginal treatment) and flat configurations (generally used to treat superficial lesions). Printer accuracy was investigated through physical measurements. The dosimetric consequences of varying the applicator’s density by tuning the printing infill percentage were analysed experimentally by measuring depth dose profiles and superficial dose distribution with Gafchromic EBT3 films (International Specialty Products, Wayne, NJ). Dose distributions were compared to those obtained with a commercial superficial applicator.ResultsMeasured printing accuracy was within 0.5 mm. Dose attenuation was not sensitive to the density of the material. Surface dose distribution comparison of the 3D printed flat applicators with respect to the commercial superficial applicator showed an overall passing rate greater than 94% for gamma analysis with 3% dose difference criteria, 3 mm distance-to-agreement criteria and 10% dose threshold.ConclusionLow-cost 3D printers are a promising solution for the customization of the HDR brachytherapy applicators. However, further assessment of 3D printing techniques and regulatory materials approval are required for clinical application.     

Virtual patient 3D dose reconstruction using in air EPID measurements and a back-projection algorithm for IMRT and VMAT treatments

PurposeAt our institute, a transit back-projection algorithm is used clinically to reconstruct in vivo patient and in phantom 3D dose distributions using EPID measurements behind a patient or a polystyrene slab phantom, respectively. In this study, an extension to this algorithm is presented whereby in air EPID measurements are used in combination with CT data to reconstruct ‘virtual’ 3D dose distributions. By combining virtual and in vivo patient verification data for the same treatment, patient-related errors can be separated from machine, planning and model errors.Methods and materialsThe virtual back-projection algorithm is described and verified against the transit algorithm with measurements made behind a slab phantom, against dose measurements made with an ionization chamber and with the OCTAVIUS 4D system, as well as against TPS patient data. Virtual and in vivo patient dose verification results are also compared.ResultsVirtual dose reconstructions agree within 1% with ionization chamber measurements. The average γ-pass rate values (3% global dose/3 mm) in the 3D dose comparison with the OCTAVIUS 4D system and the TPS patient data are 98.5 ± 1.9%(1SD) and 97.1 ± 2.9%(1SD), respectively. For virtual patient dose reconstructions, the differences with the TPS in median dose to the PTV remain within 4%.ConclusionsVirtual patient dose reconstruction makes pre-treatment verification based on deviations of DVH parameters feasible and eliminates the need for phantom positioning and re-planning. Virtual patient dose reconstructions have additional value in the inspection of in vivo deviations, particularly in situations where CBCT data is not available (or not conclusive).     

Measurement of the influence of titanium hip prosthesis on therapeutic electron beam dose distributions in a novel pelvic phantom

PurposeTo investigate the degree of 18 and 22 MeV electron beam dose perturbations caused by unilateral hip titanium (Ti) prosthesis.MethodsMeasurements were acquired using Gafchromic EBT2 film in a novel pelvic phantom made out of Nylon-12 slices in which a Ti-prosthesis is embedded. Dose perturbations were measured and compared using depth doses for 8 × 8, 10 × 10 and 11 × 11 cm² applicator-defined field sizes at 95 cm source-surface-distance (SSD). Comparisons were also made between film data at 100 cm SSD for a 10 × 10 cm² field and dose calculations made on CMS XiO treatment planning system utilizing the pencil beam algorithm. The extent of dose deviations caused by the Ti prosthesis based on film data was quantified through the dose enhancement factor (DEF), defined as the ratio of the dose influenced by the prosthesis and the unchanged beam.ResultsAt the interface between Nylon-12 and the Ti implant on the prosthesis entrance side, the dose increased to values of 21 ± 1% and 23 ± 1% for 18 and 22 MeV electron beams, respectively. DEFs increased with increasing electron energy and field size, and were found to fall off quickly with distance from the nylon-prosthesis interface. A comparison of film and XiO depth dose data for 18 and 22 MeV gave relative errors of 20% and 25%, respectively.ConclusionThis study outlines the lack of accuracy of the XiO TPS for electron planning in highly heterogeneous media. So a dosimetric error of 20–25% could influence clinical outcome.     

Modeling of couch transmission in the RayStation treatment planning system

PurposeTo present our methods and results regarding the modeling of a carbon fiber couch (Varian Exact IGRT) in the RayStation treatment planning system (TPS).MethodsThree geometrical-models (GMs) were implemented in the TPS to represent the three different regions of the couch (thick, medium and thin). The materials and densities of each GM component were tuned to maximize the agreement between measured and calculated attenuations. Moreover, a couch computed-tomography (CT) scan was acquired and dosimetrically compared with the GMs. For validation, plan-specific quality assurance (QA) of VMAT plans (TG-119 cases, 5 prostate and 5 H&N clinical cases) was performed by comparing measured dose distributions with doses computed with and without including the GMs in the TPS.ResultsCouch attenuations up to 4.3% were measured (energy: 6MV). Compared to couch CT, GMs could be modified to optimize the agreement with measurements and reduce dependence on the dose grid resolution. For both couch CT and GM, absolute deviations between measured and calculated attenuations were within 1.0%. When including the GMs in plan-specific QA, global 2%/2 mm γ-pass rates showed an average improvement of 4.8% (p-value < 0.001, max +18.6%). The couch reduced the mean dose to targets by up to 2.4% of the prescribed dose for prostate cases and up to 1.4% for H&N cases.ConclusionsRayStation accurately considers the implemented couch GMs replicating measured attenuations within an uncertainty of 1.0%. Materials and densities are proposed for the Varian Exact IGRT couch. The results obtained justify introducing couch GMs in clinical routine.     

Absolute dose verification of static intensity modulated radiation therapy (IMRT) with ion chambers of various volumes and TLD detectors

Aim

This study aims at examining absolute dose verification of step-and-shoot intensity modulated radiation treatment (IMRT) of prostate and brain patients by use of ion chambers of two different volumes and thermoluminescent detectors (TLD).

Surface dose measurements with GafChromic EBT film for 6 and 18 MV photon beams

The aim of this study was to determine the surface doses using GafChromic EBT films and compare them with plane-parallel ionization chamber measurements for 6 and 18 MV high energy photon beams. The measurements were made in a water equivalent solid phantom in the build-up region of the 6 and 18 MV photon beams at 100 cm SSD for various field sizes. Markus type plane-parallel ion chamber with fixed-separation between collecting electrodes was used to measure the percent depth doses. GafChromic EBT film measurements were performed both on the phantom surface and maximum dose depth at the same geometry with ion chamber measurements. The surface doses found using GafChromic EBT film were 15%, 20%, 29%and 39% ± 2% (1SD) for 6 MV photons, 6%, 11%, 23% and 32% ± 2% (1SD) for 18 MV photons at 5, 10, 20 and 30 cm² field sizes, respectively. GafChromic EBT film provides precise measurements for surface dose in the high energy photons. Agreement between film and plane-parallel chamber measurements was found to be within ±3% for 18 MV photon beams. There was 5% overestimate on the surface doses when compared with the plane-parallel chamber measurements for all field sizes in the 6 MV photon beams.     

Are regions of the lumbar multifidus differentially activated during walking at varied speed and inclination?

PurposeLumbar multifidus is a complex muscle with multi-fascicular morphology shown to be differentially controlled in healthy individuals during sagittal-plane motion. The normal behaviour of multifidus muscle regions during walking has only received modest attention in the literature. This study aimed to determine activation patterns for deep and superficial multifidus in young adults during walking at different speeds and inclination.MethodsThis observational cohort study evaluated ten healthy volunteers in their twenties (three women, seven men) as they walked on a treadmill in eight conditions; at 2 km/h and 4 km/h, each at 0, 1, 5, and 10% inclination. Intramuscular EMG was recorded from the deep and superficial multifidus unilaterally at L5. Activity was characterized by: amplitude of the peak of activation, position of peak within the gait cycle (0–100%), and duration relative to the full gait cycle.ResultsAcross all conditions superficial multifidus showed higher normalised EMG amplitude (p < 0.01); superficial multifidus peak amplitude was 232 ± 115% higher when walking at 4 km/h/10%, versus only 172 ± 77% higher for deeper region (p < 0.01). The percentage of the gait cycle where peak EMG amplitude was detected did not differ between regions (49 ± 13%). Deep multifidus duration of activation was longer when walking at the faster vs slower speed at all inclinations (p < 0.01), which was not evident for superficial multifidus (p < 0.05). Thus, a significantly longer activation of deep multifidus was observed compared to superficial multifidus when walking at 4 km/h (p < 0.05).ConclusionsDifferential activation within lumbar multifidus was shown in young adults during walking. The prolonged, more tonic activation of deep relative to superficial regions of multifidus during gait supports a postural function of deeper fibres.     

Application of Gafchromic EBT2 film for intraoperative radiation therapy quality assurance

PurposeIntraoperative radiation therapy (IORT) using electron beam is commonly done by mobile dedicated linacs that have a variable range of electron energies. This paper focuses on the evaluation of the EBT2 film response in the green and red colour channels for IORT quality assurance (QA).MethodsThe calibration of the EBT2 films was done in two ranges; 0–8 Gy for machine QA by red channel and 8–24 Gy for patient-specific QA by green channel analysis. Irradiation of calibration films and relative dosimetries were performed in a water phantom. To evaluate the accuracy of the film response in relative dosimetry, gamma analysis was used to compare the results of the Monte Carlo simulation and ionometric dosimetry. Ten patients with early stage breast cancer were selected for in-vivo dosimetry using the green channel of the EBT2 film.ResultsThe calibration curves were obtained by linear fitting of the green channel and a third-order polynomial function in the red channel (R2 = 0.99). The total dose uncertainty was up to 4.2% and 4.7% for the red and green channels, respectively. There was a good agreement between the relative dosimetries of films by the red channel, Monte Carlo simulations and ionometric values. The mean dose difference of the in-vivo dosimetry by green channel of this film and the expected values was about 1.98% ± 0.75.ConclusionThe results of this study showed that EBT2 film can be considered as an appropriate tool for machine and patient-specific QA in IORT.     

Evaluation of the target dose coverage of stereotactic body radiotherapy for lung cancer using helical tomotherapy: A dynamic phantom study

AimTo evaluate the target dose coverage for lung stereotactic body radiotherapy (SBRT) using helical tomotherapy (HT) with the internal tumor volume (ITV) margin settings adjusted according to the degree of tumor motion.BackgroundLung SBRT with HT may cause a dosimetric error when the target motion is large.Materials and methodsTwo lung SBRT plans were created using a tomotherapy planning station. Using these original plans, five plans with different ITV margins (4.0–20.0 mm for superior-inferior [SI] dimension) were generated. To evaluate the effects of respiratory motion on HT, an original dynamic motion phantom was developed. The respiratory wave of a healthy volunteer was used for dynamic motion as the typical tumor respiratory motion. Five patterns of motion amplitude that corresponded to five ITV margin sizes and three breathing cycles of 7, 14, and 28 breaths per minute were used. We evaluated the target dose change between a static delivery and a dynamic delivery with each motion pattern.ResultsThe target dose difference increased as the tumor size decreased and as the tumor motion increased. Although a target dose difference of <5 % was observed at ≤10 mm of tumor motion for each condition, a maximum difference of -9.94 % ± 7.10 % was observed in cases of small tumors with 20 mm of tumor motion under slow respiration.ConclusionsMinimizing respiratory movement is recommended as much as possible for lung SBRT with HT, especially for cases involving small tumors.     

NTCP modeling analysis of acute hematologic toxicity in whole pelvic radiation therapy for gynecologic malignancies – A dosimetric comparison of IMRT and spot-scanning proton therapy (SSPT)

PurposeThis treatment planning study was conducted to determine whether spot scanning proton beam therapy (SSPT) reduces the risk of grade ⩾3 hematologic toxicity (HT3+) compared with intensity modulated radiation therapy (IMRT) for postoperative whole pelvic radiation therapy (WPRT).Methods and materialsThe normal tissue complication probability (NTCP) of the risk of HT3+ was used as an in silico surrogate marker in this analysis. IMRT and SSPT plans were created for 13 gynecologic malignancy patients who had received hysterectomies. The IMRT plans were generated using the 7-fields step and shoot technique. The SSPT plans were generated using anterior-posterior field with single field optimization. Using the relative biological effectives (RBE) value of 1.0 for IMRT and 1.1 for SSPT, the prescribed dose was 45 Gy(RBE) in 1.8 Gy(RBE) per fractions for 95% of the planning target volume (PTV). The homogeneity index (HI) and the conformity index (CI) of the PTV were also compared.ResultsThe bone marrow (BM) and femoral head doses using SSPT were significantly lower than with IMRT. The NTCP modeling analysis showed that the risk of HT3+ using SSPT was significantly lower than with IMRT (NTCP = 0.04 ± 0.01 and 0.19 ± 0.03, p = 0.0002, respectively). There were no significant differences in the CI and HI of the PTV between IMRT and SSPT (CI = 0.97 ± 0.01 and 0.96 ± 0.02, p = 0.3177, and HI = 1.24 ± 0.11 and 1.27 ± 0.05, p = 0.8473, respectively).ConclusionThe SSPT achieves significant reductions in the dose to BM without compromising target coverage, compared with IMRT. The NTCP value for HT3+ in SSPT was significantly lower than in IMRT.     

In vivo skin dose measurement using MOSkin detectors in tangential breast radiotherapy

The purpose of this study is to measure patient skin dose in tangential breast radiotherapy. Treatment planning dose calculation algorithm such as Pencil Beam Convolution (PBC) and in vivo dosimetry techniques such as radiochromic film can be used to accurately monitor radiation doses at tissue depths, but they are inaccurate for skin dose measurement. A MOSFET-based (MOSkin) detector was used to measure skin dose in this study. Tangential breast radiotherapies (“bolus” and “no bolus”) were simulated on an anthropomorphic phantom and the skin doses were measured. Skin doses were also measured in 13 patients undergoing each of the techniques. In the patient study, the EBT2 measurements and PBC calculation tended to over-estimate the skin dose compared with the MOSkin detector (p < 0.05) in the “no bolus radiotherapy”. No significant differences were observed in the “bolus radiotherapy” (p > 0.05). The results from patients were similar to that of the phantom study. This shows that the EBT2 measurement and PBC calculation, while able to predict accurate doses at tissue depths, are inaccurate in predicting doses at build-up regions. The clinical application of the MOSkin detectors showed that the average total skin doses received by patients were 1662 ± 129 cGy (medial) and 1893 ± 199 cGy (lateral) during “no bolus radiotherapy”. The average total skin doses were 4030 ± 72 cGy (medial) and 4004 ± 91 cGy (lateral) for “bolus radiotherapy”. In some cases, patient skin doses were shown to exceed the dose toxicity level for skin erythema. Hence, a suitable device for in vivo dosimetry is necessary to accurately determine skin dose.     

Application of Monte Carlo calculations for validation of a treatment planning system in high dose rate brachytherapy

AimThe accuracy of treatment planning systems is of vital importance in treatment outcomes in brachytherapy. In the current study the accuracy of dose calculations of a high dose rate (HDR) brachytherapy treatment planning system (TPS) was validated using the Monte Carlo method.Materials and methodsThree ⁶⁰Co sources of the GZP6 afterloading brachytherapy system were modelled using MCNP4C Monte Carlo (MC) code. The dose distribution around all the sources was calculated by MC and a dedicated treatment planning system. The results of both methods were compared.ResultsThere was good agreement (<2%) between TPS and MC calculated dose distributions except at a point near the sources (<1 cm) and beyond the tip of the sources.ConclusionsOur study confirmed the accuracy of TPS calculated dose distributions for clinical use in HDR brachytherapy.     

Investigating the spatial accuracy of CBCT-guided cranial radiosurgery: A phantom end-to-end test study

PurposeTo evaluate the spatial accuracy of a frameless cone-beam computed tomography (CBCT)-guided cranial radiosurgery (SRS) using an end-to-end (E2E) phantom test methodology.Methods and materialsFive clinical SRS plans were mapped to an acrylic phantom containing a radiochromic film. The resulting phantom-based plans (E2E plans) were delivered four times. The phantom was setup on the treatment table with intentional misalignments, and CBCT-imaging was used to align it prior to E2E plan delivery. Comparisons (global gamma analysis) of the planned and delivered dose to the film were performed using a commercial triple-channel film dosimetry software. The necessary distance-to-agreement to achieve a 95% (DTA95) gamma passing rate for a fixed 3% dose difference provided an estimate of the spatial accuracy of CBCT-guided SRS. Systematic (∑) and random (σ) error components, as well as 95% confidence levels were derived for the DTA95 metric.ResultsThe overall systematic spatial accuracy averaged over all tests was 1.4 mm (SD: 0.2 mm), with a corresponding 95% confidence level of 1.8 mm. The systematic (Σ) and random (σ) spatial components of the accuracy derived from the E2E tests were 0.2 mm and 0.8 mm, respectively.ConclusionsThe E2E methodology used in this study allowed an estimation of the spatial accuracy of our CBCT-guided SRS procedure. Subsequently, a PTV margin of 2.0 mm is currently used in our department.     

Full Monte Carlo and measurement-based overall performance assessment of improved clinical implementation of eMC algorithm with emphasis on lower energy range

New version 13.6.23 of the electron Monte Carlo (eMC) algorithm in Varian Eclipse™ treatment planning system has a model for 4 MeV electron beam and some general improvements for dose calculation. This study provides the first overall accuracy assessment of this algorithm against full Monte Carlo (MC) simulations for electron beams from 4 MeV to 16 MeV with most emphasis on the lower energy range. Beams in a homogeneous water phantom and clinical treatment plans were investigated including measurements in the water phantom. Two different material sets were used with full MC: (1) the one applied in the eMC algorithm and (2) the one included in the Eclipse™ for other algorithms. The results of clinical treatment plans were also compared to those of the older eMC version 11.0.31. In the water phantom the dose differences against the full MC were mostly less than 3% with distance-to-agreement (DTA) values within 2 mm. Larger discrepancies were obtained in build-up regions, at depths near the maximum electron ranges and with small apertures. For the clinical treatment plans the overall dose differences were mostly within 3% or 2 mm with the first material set. Larger differences were observed for a large 4 MeV beam entering curved patient surface with extended SSD and also in regions of large dose gradients. Still the DTA values were within 3 mm. The discrepancies between the eMC and the full MC were generally larger for the second material set. The version 11.0.31 performed always inferiorly, when compared to the 13.6.23.     

“Edge-on” MOSkin detector for stereotactic beam measurement and verification

Dosimetry in small radiation field is challenging and complicated because of dose volume averaging and beam perturbations in a detector. We evaluated the suitability of the “Edge-on” MOSkin (MOSFET) detector in small radiation field measurement. We also tested the feasibility for dosimetric verification in stereotactic radiosurgery (SRS) and stereotactic radiotherapy (SRT). “Edge-on” MOSkin detector was calibrated and the reproducibility and linearity were determined. Lateral dose profiles and output factors were measured using the “Edge-on” MOSkin detector, ionization chamber, SRS diode and EBT2 film. Dosimetric verification was carried out on two SRS and five SRT plans. In dose profile measurements, the “Edge-on” MOSkin measurements concurred with EBT2 film measurements. It showed full width at half maximum of the dose profile with average difference of 0.11 mm and penumbral width with difference of ±0.2 mm for all SRS cones as compared to EBT2 film measurement. For output factor measurements, a 1.1% difference was observed between the “Edge-on” MOSkin detector and EBT2 film for 4 mm SRS cone. The “Edge-on” MOSkin detector provided reproducible measurements for dose verification in real-time. The measured doses concurred with the calculated dose for SRS (within 1%) and SRT (within 3%). A set of output correction factors for the “Edge-on” MOSkin detector for small radiation fields were derived from EBT2 film measurement and presented. This study showed that the “Edge-on” MOSkin detector is a suitable tool for dose verification in small radiation field.