Multi-level multi-leaf collimators: Optimization of layer thicknesses and a feasibility study.

PURPOSE: The function of multi-leaf collimators (MLC) is to modulate and shape the intensity of a radiotherapy beam by either blocking or unblocking beamlets. A variation on this functionality is tested in this work wherein the MLC is split into layers, with each layer attenuating the beam by a different amount. In this design, full blocking of a beamlet occurs only if all layers are blocked. This work suggests that such a device, a multi-layer MLC (MLMLC), can deliver dose distributions like a single layer MLC can deliver while requiring less time and monitor units (MU) METHODS: Optimal fluences were made for prostate plans using the Eclipse v13.6. An algorithm was developed to create step-and-shoot MLMLC patterns to match these optimal fluences when using up to six layers of MLC. Twelve MLMLC plans were made in total. These patterns were imported back into Eclipse as equivalent tungsten compensators and doses were calculated. Dose-volume histogram (DVH) values, total monitor units (MU), and total time to deliver were compared between arc-style MLMLC plans and nine-field step and shoot IMRT plans created completely in Eclipse using a single layer MLC. RESULTS: When using three or more layers, specified DVH values between the two sets agreed to within 5% while requiring roughly half as much time to deliver and about 20% fewer MU. CONCLUSIONS: Demonstrated that having multi-layer MLC can deliver dose distributions like a single layer MLC with less time and monitor units.

Dosimetric effects of the kV based image-guided radiation therapy of prone breast external beam radiation: Towards the optimized imaging frequency.

PURPOSE: For prone breast treatment, daily image-guided radiation therapy (IGRT) allows couch shifting to correct breast position relative to the treatment field. This work investigates the dosimetric effect of reducing kV imaging frequencies and the feasibility of optimizing the frequency using patient anatomy or their first 3-day shifts. METHOD: Thirty-seven prone breast patients who had been treated with skin marker alignment followed by daily kV were retrospectively analyzed. Three IGRT schemes (daily-kV, weekly-kV, no-kV) were simulated, assuming that fractions with kV imaging deliver a dose distribution equivalent to that in computed tomography (CT) planning, whereas other fractions yield a dose distribution as recreated by shifting the CT plan isocenter back to its position before the couch shift was applied. Treatment dose to targets (breast and lumpectomy cavity [LPC]) and organs at risks (OAR)s (heart, ipsilateral lung) in different schemes were calculated. Patient anatomy information on CT plans and first 3-day couch shift data were analyzed to investigate whether these factors could guide imaging scheme optimization. RESULTS: When kV imaging frequency was reduced, the percentage dose changes (δD) for breast and LPC objectives (average <1%) were smaller than those for heart and lung (average 28%-31% for Dmean ). In general, the δD of no-kV imaging was approximately that of weekly kV imaging × a factor of 1.2-1.4. Although most dose objectives were not affected, the potential higher heart dose may be of concern. No strong correlation was found between δD for different kV frequencies and patient anatomy size/distance or the first 3-day couch shift data. CONCLUSIONS: Despite resulting in lower imaging dose, time, cost, and similar target coverage, a reduction in kV imaging frequency may introduce higher heart complication risk. Daily kVs are needed more in left-sided breast patients. A less frequent imaging schedule, if considered, cannot be individually optimized using CT anatomic features or early shift data.

Evaluation of Computed Tomography Scanners for Feasibility of Using Averaged Hounsfield Unit-to-Stopping Power Ratio Calibration Curve.

PURPOSE: The purpose of this study was to quantify the variability of stoichiometric calibration curves for different computed tomography (CT) scanners and determine whether an averaged Hounsfield unit (HU)-to-stopping power ratio (SPR) calibration curve can be used across multiple CT scanners. MATERIALS AND METHODS: Five CT scanners were used to scan an electron density phantom to establish HU values of known material plugs. A stoichiometric calibration curve was calculated for CT scanners and for the average curve. Animal tissue surrogates were used to compare the water-equivalent thickness (WET) of the animal tissue surrogates calculated by the treatment planning system (TPS) and the WET values measured with a multilayered ionization chamber. The calibration curves were optimized to reduce the percentage of difference between measured and TPS-calculated WET values. A second set of tissue surrogates was then used to evaluate the overall range of uncertainty for the optimized CT-specific and average calibration curves. RESULTS: Overall, the average variation in HU for all 6 calibration curves before optimization was 8.3 HU. For both the averaged and CT-specific calibrations, the root mean square error (RMSE) of the percentage of difference between TPS-calculated and measured WET values before optimization was 4%. The RMSE of the percentage of difference for the TPS-calculated and multilayered ionization chamber measured WET values after the optimization for both averaged and CT-specific calibration curves was reduced to less than 1.5%. The overall RMSE of the TPS and the measured WET percentage of difference after optimization was 2.1% for both averaged and CT-specific calibration curves. CONCLUSION: Averaged CT calibration curves can be used to map the HU-to-SPR in TPSs, if the variations in HU values across all scanners is relatively small. Performing tissue surrogate optimization of the HU-to-SPR calibration curve has been shown to reduce the overall uncertainty of the calibration for averaged and CT-specific calibration curves and is recommended, especially if an averaged HU-to-SPR calibration curve is used.

Rectal dose to prostate cancer patients treated with proton therapy with or without rectal spacer.

The purpose of this study was to evaluate whether a spacer inserted in the prerectal space could reduce modeled rectal dose and toxicity rates for patients with prostate cancer treated in silico with pencil beam scanning (PBS) proton therapy. A total of 20 patients were included in this study who received photon therapy (12 with rectal spacer (DuraSeal™ gel) and 8 without). Two PBS treatment plans were retrospectively created for each patient using the following beam arrangements: (1) lateral-opposed (LAT) fields and (2) left and right anterior oblique (LAO/RAO) fields. Dose volume histograms (DVH) were generated for the prostate, rectum, bladder, and right and left femoral heads. The normal tissue complication probability (NTCP) for ≥grade 2 rectal toxicity was calculated using the Lyman-Kutcher-Burman model and compared between patients with and without the rectal spacer. A significantly lower mean rectal DVH was achieved in patients with rectal spacer compared to those without. For LAT plans, the mean rectal V70 with and without rectal spacer was 4.19 and 13.5%, respectively. For LAO/RAO plans, the mean rectal V70 with and without rectal spacer was 5.07 and 13.5%, respectively. No significant differences were found in any rectal dosimetric parameters between the LAT and the LAO/RAO plans generated with the rectal spacers. We found that ≥ 9 mm space resulted in a significant decrease in NTCP modeled for ≥grade 2 rectal toxicity. Rectal spacers can significantly decrease modeled rectal dose and predicted ≥grade 2 rectal toxicity in prostate cancer patients treated in silico with PBS. A minimum of 9 mm separation between the prostate and anterior rectal wall yields the largest benefit.

An analytical formalism to calculate phantom scatter factors for flattening filter free (FFF) mode photon beams.

Phantom Scatter Factors, Sp in the Khan formalism (Khan et al 1980 J. Radiat. Oncol. Biol. Phys. 6 745-51) describe medium-induced changes in photon-beam intensity as a function of size of the beam. According to the British Journal of Radiology, Supplement 25, megavoltage phantom scatter factors are invariant as a function of photon-beam energy. However, during the commissioning of an accelerator with flattening filter free (FFF) photon beams (Varian TrueBeam(TM) 6-MV FFF and 10-MV FFF), differences were noted in phantom scatter between the filtered beams and FFF-mode beams. The purpose of this work was to evaluate this difference and provide an analytical formalism to explain the phantom scatter differences between FFF-mode and the filtered mode. An analytical formalism was devised to demonstrate the source of phantom scatter differences between the filtered and the FFF-mode beams. The reason for the differences in the phantom scatter factors between the filtered and the FFF-mode beams is hypothesized to be the non-uniform beam profiles of the FFF-mode beams. The analytical formalism proposed here is based on this idea, taking the product of the filtered phantom scatter factors and the ratio of the off-axis ratio between the FFF-mode and the filtered beams. All measurements were performed using a Varian TrueBeam(TM) linear accelerator with photon energies of 6-MV and 10-MV in both filtered and FFF-modes. For all measurements, a PTW Farmer type chamber and a Scanditronix CC04 cylindrical ionization were used. The in-water measurements were made at depth dose maximum and 100 cm source-to-axis distance. The in-air measurements were done at 100 cm source-to-axis distance with appropriate build-up cap. From these measurements, the phantom scatter factors were derived for the filtered beams and the FFF-mode beams for both energies to be evaluated against the phantoms scatter factors calculated using the proposed algorithm. For 6-MV, the difference between the measured and the calculated FFF-mode phantom scatter factors ranged from -0.34% to 0.73%. The average per cent difference was -0.17% (1σ = 0.25%). For 10-MV, the difference ranged from -0.19% to 0.24%. The average per cent difference was -0.17% (1σ = 0.13%). An analytical formalism was presented to calculate the phantom scatter factors for FFF-mode beams using filtered phantom scatter factors as a basis. The overall differences between measurements and calculations were within ± 0.5% for 6-MV and ± 0.25% for 10-MV.

Evaluation of dose variation to normal and critical structures for lung hypofractionated stereotactic body radiation therapy.

PURPOSE: To quantify the dose received by normal and critical structures during lung stereotactic body radiation therapy (SBRT) when registered to tumor or bone. METHODS AND MATERIALS: Sixteen patients with lung cancer receiving a total dose of 50 Gy in 4fractions for lung SBRT were retrospectively studied. Cone-beam computed tomography (CT) was performed for all fractions, and the images obtained were registered with planning CT with respect tosoft tissue for target localization. Isocenter shifts were determined for each fraction from differences between the bony and tumor alignments; doses were then recalculated based on the new isocenters and summed over all 4 fractions to compare against the planned normal and critical tissue dose. The normal and critical structures evaluated were total and ipsilateral lung, spinal cord, and esophagus. The first data collected were isocenter coordinate shifts in all 3 Cartesian coordinates for both tumor andbony alignments. The second were the dose differences to the normal and critical structures fromthe planned and recalculated doses for alignment based on the tumor. RESULTS: The study showed that while the maximum isocenter coordinate shifts in any direction couldbe as much as 1.60 cm, the normal and critical structure dose variations between the original plans and the simulated plans showed almost no change. The mean volume of total lung that receivedat least 20Gy difference for total lung and ipsilateral lung were 0.01% and -0.04%, respectively. For the esophagus, spinal cord, and heart the maximum and mean dose differences were 0.25 Gy and -0.04 Gy, -0.08 Gy and -0.02 Gy, and 0.02 Gy and 0.05 Gy, respectively. CONCLUSIONS: Target localization using daily cone-beam CT with soft tissue registration was appropriate for minimizing the dose to the normal and critical structures without the need to re-plan due to the changes in the tumor position. For tumors located close to a critical structure, daily cone-beam CT is recommended to determine the appropriate isocenter shifts.

Feasibility of using two-dimensional array dosimeter for in vivo dose reconstruction via transit dosimetry.

Two-dimensional array dosimeters are commonly used to perform pretreatment quality assurance procedures, which makes them highly desirable for measuring transit fluences for in vivo dose reconstruction. The purpose of this study was to determine if an in vivo dose reconstruction via transit dosimetry using a 2D array dosimeter was possible. To test the accuracy of measuring transit dose distribution using a 2D array dosimeter, we evaluated it against the measurements made using ionization chamber and radiochromic film (RCF) profiles for various air gap distances (distance from the exit side of the solid water slabs to the detector distance; 0 cm, 30 cm, 40 cm, 50 cm, and 60 cm) and solid water slab thicknesses (10 cm and 20 cm). The backprojection dose reconstruction algorithm was described and evaluated. The agreement between the ionization chamber and RCF profiles for the transit dose distribution measurements ranged from -0.2% ~ 4.0% (average 1.79%). Using the backprojection dose reconstruction algorithm, we found that, of the six conformal fields, four had a 100% gamma index passing rate (3%/3 mm gamma index criteria), and two had gamma index passing rates of 99.4% and 99.6%. Of the five IMRT fields, three had a 100% gamma index passing rate, and two had gamma index passing rates of 99.6% and 98.8%. It was found that a 2D array dosimeter could be used for backprojection dose reconstruction for in vivo dosimetry.

High-precision GAFCHROMIC EBT film-based absolute clinical dosimetry using a standard flatbed scanner without the use of a scanner non-uniformity correction.

To report a study of the use of GAFCHROMIC EBT radiochromic film (RCF) digitized with a commercially available flatbed document scanner for accurate and reliable all-purpose two-dimensional (2D) absolute dosimetry within a clinical environment. We used a simplified methodology that yields high-precision dosimetry measurements without significant postirradiation correction. The Epson Expression 1680 Professional scanner and the Epson Expression 10000XL scanner were used to digitize the films. Both scanners were retrofitted with light-diffusing glass to minimize the effects of Newton rings. Known doses were delivered to calibration films. Flat and wedge fields were irradiated with variable depth of solid water and 5 cm back scatter solid water. No particular scanner nonuniformity effect corrections or significant post-scan image processing were carried out. The profiles were compared with CC04 ionization chamber profiles. The depth dose distribution was measured at a source-to-surface distance (SSD) of 100 cm with a field size of 10 x 10 cm2. Additionally, 22 IMRT fields were measured and evaluated using gamma index analysis. The overall accuracy of RCF with respect to CC04 was found to be 2%-4%. The overall accuracy of RCF was determined using the absolute mean of difference for all flat and wedge field profiles. For clinical IMRT fields, both scanners showed an overall gamma index passing rate greater than 90%. This work demonstrated that EBT films, in conjunction with a commercially available flatbed scanner, can be used as an accurate and precise absolute dosimeter. Both scanners showed that no significant scanner nonuniformity correction is necessary for accurate absolute dosimetry using the EBT films for field sizes smaller than or equal to 15 x 15 cm2.

Evaluation of similarity measures for use in the intensity-based rigid 2D-3D registration for patient positioning in radiotherapy.

PURPOSE: Rigid 2D-3D registration is an alternative to 3D-3D registration for cases where largely bony anatomy can be used for patient positioning in external beam radiation therapy. In this article, the authors evaluated seven similarity measures for use in the intensity-based rigid 2D-3D registration using a variation in Skerl's similarity measure evaluation protocol. METHODS: The seven similarity measures are partitioned intensity uniformity, normalized mutual information (NMI), normalized cross correlation (NCC), entropy of the difference image, pattern intensity (PI), gradient correlation (GC), and gradient difference (GD). In contrast to traditional evaluation methods that rely on visual inspection or registration outcomes, the similarity measure evaluation protocol probes the transform parameter space and computes a number of similarity measure properties, which is objective and optimization method independent. The variation in protocol offers an improved property in the quantification of the capture range. The authors used this protocol to investigate the effects of the downsampling ratio, the region of interest, and the method of the digitally reconstructed radiograph (DRR) calculation [i.e., the incremental ray-tracing method implemented on a central processing unit (CPU) or the 3D texture rendering method implemented on a graphics processing unit (GPU)] on the performance of the similarity measures. The studies were carried out using both the kilovoltage (kV) and the megavoltage (MV) images of an anthropomorphic cranial phantom and the MV images of a head-and-neck cancer patient. RESULTS: Both the phantom and the patient studies showed the 2D-3D registration using the GPU-based DRR calculation yielded better robustness, while providing similar accuracy compared to the CPU-based calculation. The phantom study using kV imaging suggested that NCC has the best accuracy and robustness, but its slow function value change near the global maximum requires a stricter termination condition for an optimization method. The phantom study using MV imaging indicated that PI, GD, and GC have the best accuracy, while NCC and NMI have the best robustness. The clinical study using MV imaging showed that NCC and NMI have the best robustness. CONCLUSIONS: The authors evaluated the performance of seven similarity measures for use in 2D-3D image registration using the variation in Skerl's similarity measure evaluation protocol. The generalized methodology can be used to select the best similarity measures, determine the optimal or near optimal choice of parameter, and choose the appropriate registration strategy for the end user in his specific registration applications in medical imaging.

Experimental investigation of the implementation of Fermi-Eyges-Hogstrom electron beam model of the Pinnacle3 system at extended SSDs.

A commercially available ADAC Pinnacle(3) radiation treatment planning system has been used to model electron beams from a Varian Clinac 2300C/D in the energy range of 6 to 22 MeV. Prior to clinical use, the dosimetric characteristics of the beams have to be modeled accurately. As a first step for beam modeling, a number of dose profile and depth dose measurements were taken at standard source-to-surface distance (SSD) of 100 cm. Dose profiles and depth dose measurements at extended SSDs up to 120 cm are important for ascertaining accuracy of the model, as well as their clinical usefulness in the treatment of some sites (e.g., head-and-neck tumors). Modeled and measured beam data were compared. Over 98% of comparison points (modeled vs. measured) at 100-cm SSD were within 2.5% or 2.5 mm. At 110 cm SSD, over 98% of compared points were within 4% or 4 mm, and at 120-cm SSD, over 98% of compared points were within 5% or 5 mm. Overall, more than 98% of compared points were within 4% or 4 mm. Better models were produced for lower energies (6 to 15 MeV) than higher energies (18 and 22 MeV). For 6, 9, 12, and 15 MeV, 89% of compared points were within 2% or 2 mm. For 18- and 22-MeV electron energies, 75% and 67%, respectively, were within 2% or 2 mm. These results are consistent with the recommendations of AAPM Task Group Report 53.

Dose variations with varying calculation grid size in head and neck IMRT.

Ever since the advent and development of treatment planning systems, the uncertainty associated with calculation grid size has been an issue. Even to this day, with highly sophisticated 3D conformal and intensity-modulated radiation therapy (IMRT) treatment planning systems (TPS), dose uncertainty due to grid size is still a concern. A phantom simulating head and neck treatment was prepared from two semi-cylindrical solid water slabs and a radiochromic film was inserted between the two slabs for measurement. Plans were generated for a 5,400 cGy prescribed dose using Philips Pinnacle(3) TPS for two targets, one shallow ( approximately 0.5 cm depth) and one deep ( approximately 6 cm depth). Calculation grid sizes of 1.5, 2, 3 and 4 mm were considered. Three clinical cases were also evaluated. The dose differences for the varying grid sizes (2 mm, 3 mm and 4 mm from 1.5 mm) in the phantom study were 126 cGy (2.3% of the 5,400 cGy dose prescription), 248.2 cGy (4.6% of the 5,400 cGy dose prescription) and 301.8 cGy (5.6% of the 5,400 cGy dose prescription), respectively for the shallow target case. It was found that the dose could be varied to about 100 cGy (1.9% of the 5,400 cGy dose prescription), 148.9 cGy (2.8% of the 5,400 cGy dose prescription) and 202.9 cGy (3.8% of the 5,400 cGy dose prescription) for 2 mm, 3 mm and 4 mm grid sizes, respectively, simply by shifting the calculation grid origin. Dose difference with a different range of the relative dose gradient was evaluated and we found that the relative dose difference increased with an increase in the range of the relative dose gradient. When comparing varying calculation grid sizes and measurements, the variation of the dose difference histogram was insignificant, but a local effect was observed in the dose difference map. Similar results were observed in the case of the deep target and the three clinical cases also showed results comparable to those from the phantom study.

Evaluation of surface and build-up region dose for intensity-modulated radiation therapy in head and neck cancer.

Despite much development, there remains dosimetric uncertainty in the surface and build-up regions in intensity-modulated radiation therapy treatment plans for head and neck cancers. Experiments were performed to determine the dosimetric discrepancies in the surface and build-up region between the treatment planning system (TPS) prediction and experimental measurement using radiochromic film. A head and neck compression film phantom was constructed from two semicylindrical solid water slabs. Treatment plans were generated using two commercial TPSs (PINNACLE3 and CORVUS) for two cases, one with a shallow (approximately 0.5 cm depth) target and another with a deep (approximately 6 cm depth) target. The plans were evaluated for a 54 Gy prescribed dose. For each case, two pieces of radiochromic film were used for dose measurement. A small piece of film strip was placed on the surface and another was inserted within the phantom. Overall, both TPSs showed good agreement with the measurement. For the shallow target case, the dose differences were within +/- 300 cGy (5.6% with respect to the prescribed dose) for PINNACLE3 and +/- 240 cGy (4.4%) for CORVUS in 90% of the region of interest. For the deep target case, the dose differences were +/- 350 (6.5%) for PINNACLE3 and +/- 260 cGy (4.8%) for CORVUS in 90% of the region of interest. However, it was found that there were significant discrepancies from the surface to about 0.2 cm in depth for both the shallow and deep target cases. It was concluded that both TPSs overestimated the surface dose for both shallow and deep target cases. The amount of overestimation ranges from 400 to 1000 cGy (approximately 7.4% to 18.5% with respect to the prescribed dose, 5400 cGy).