The Multigaussian method: a new approach to mitigating spatial heterogeneities with multichannel radiochromic film dosimetry

An open-source development based on photogrammetry for a real-time IORT treatment planning system

PURPOSE

This study presents a treatment planning system for intraoperative low-energy photon radiotherapy based on photogrammetry from real images of the surgical site taken in the operating room.

MATERIAL AND METHODS

The study population comprised 15 patients with soft-tissue sarcoma. The system obtains the images of the area to be irradiated with a smartphone or tablet, so that the absorbed doses in the tissue can be calculated from the reconstruction without the need for computed tomography. The system was commissioned using 3D printing of the reconstructions of the tumor beds. The absorbed doses at various points were verified using radiochromic films that were suitably calibrated for the corresponding energy and beam quality.

RESULTS

The average reconstruction time of the 3D model from the video sequence in the 15 patients was 229,6±7,0 s. The entire procedure, including video capture, reconstruction, planning, and dose calculation was 520,6±39,9 s. Absorbed doses were measured on the 3D printed model with radiochromic film, the differences between these measurements and those calculated by the treatment planning system were 1.4% at the applicator surface, 2.6% at 1 cm, 3.9% at 2 cm and 6.2% at 3 cm.

CONCLUSIONS

The study shows a photogrammetry-based low-energy photon IORT planning system, capable of obtaining real-time images inside the operating room, immediately after removal of the tumor and immediately before irradiation. The system was commissioned with radiochromic films measurements in 3D-printed model.

Monte Carlo-based independent dose verification of radiosurgery HyperArc plans

PURPOSE

To investigate the feasibility of using the free PRIMO Monte Carlo software for independent dose check of cranial SRS plans designed with the Varian HyperArc (HA) technique.

MATERIALS AND METHODS

In this study, the PRIMO Monte Carlo software v. 0.3.64.1800 was used with the phase-space files (v. 2, Feb. 27, 2013) provided by Varian for 6 MV flattening-filter-free (FFF) photon beams from a Varian TrueBeam linear accelerator (linac), equipped with a Millennium 120 multileaf collimator (MLC). This configuration was validated by comparing the percentage depth doses (PDDs), lateral profiles and relative output factors (OFs) simulated in a water phantom against measurements for field sizes from 1 × 1 to 40 × 40 cm². The agreement between simulated and experimental relative dose curves was evaluated using a global (G) gamma index analysis. In addition, the accuracy of PRIMO to model the MLC was investigated (dosimetric leaf gap, tongue and groove, leaf transmission and interleaf leakage). Thirty-five HA SRS plans computed in the Eclipse treatment planning system (TPS) were simulated in PRIMO. The Acuros XB algorithm v. 16.10 (dose to medium) was used in Eclipse. Sixty targets with diameters ranging from 6 to 33 mm were included. Agreement between the dose distributions given by Eclipse and PRIMO was evaluated in terms of 3D global gamma passing rates (GPRs) for the 2 %/2 mm criteria.

RESULTS

Average GPR greater than 95 % with the 2 %(G)/1 mm criteria were obtained over the PDD and profiles of each field size. Differences between PRIMO calculated and measured OFs were within 0.5 % in all fields, except for the 1 × 1 cm² with a discrepancy of 1.5 %. Regarding the MLC modeling in PRIMO, an agreement within 3 % was achieved between calculated and experimental doses. Excellent agreement between PRIMO and Eclipse was found for the 35 HA plans. The 3D global GPRs (2 %/2 mm) for the targets and external patient contour were 99.6 % ± 1.1 % and 99.8 % ± 0.5 %, respectively.

CONCLUSIONS

According to the results described in this study, the PRIMO Monte Carlo software, in conjunction with the 6X FFF Varian phase-space files, can be used as secondary dose calculation software to check stereotactic radiosurgery plans from Eclipse using the HyperArc technique.

A protocol for accurate radiochromic film dosimetry using Radiochromic.com

BACKGROUND

Radiochromic films have many applications in radiology and radiation therapy. Generally, the dosimetry system for radiochromic film dosimetry is composed of radiochromic films, flatbed scanner, and film analysis software. The purpose of this work is to present the effectiveness of a protocol for accurate radiochromic film dosimetry using Radiochromic.com as software for film analysis.

MATERIALS AND METHODS

Procedures for image acquisition, lot calibration, and dose calculation are explained and analyzed. Radiochromic.com enables state-of-the-art models and corrections for radiochromic film dosimetry, such as the Multigaussian model for multichannel film dosimetry, and lateral, inter-scan, and re-calibration corrections of the response.

RESULTS

The protocol presented here provides accurate dose results by mitigating the sources of uncertainty that affect radiochromic film dosimetry.

CONCLUSIONS

Appropriate procedures for film and scanner handling in combination with Radiochromic.com as software for film analysis make easy and accurate radiochromic film dosimetry feasible.

Validation of virtual water phantom software for pre-treatment verification of single-isocenter multiple-target stereotactic radiosurgery

The aim of this study was to benchmark the accuracy of the VIrtual Phantom Epid dose Reconstruction (VIPER) software for pre-treatment dosimetric verification of multiple-target stereotactic radiosurgery (SRS). VIPER is an EPID-based method to reconstruct a 3D dose distribution in a virtual phantom from in-air portal images. Validation of the VIPER dose calculation was assessed using several MLC-defined fields for a 6 MV photon beam. Central axis percent depth doses (PDDs) and output factors were measured with an ionization chamber in a water tank, while dose planes at a depth of 10 cm in a solid flat phantom were acquired with radiochromic films. The accuracy of VIPER for multiple-target SRS plan verification was benchmarked against Monte Carlo simulations. Eighteen multiple-target SRS plans designed with the Eclipse treatment planning system were mapped to a cylindrical water phantom. For each plan, the 3D dose distribution reconstructed by VIPER within the phantom was compared with the Monte Carlo simulation, using a 3D gamma analysis. Dose differences (VIPER vs. measurements) generally within 2% were found for the MLC-defined fields, while film dosimetry revealed gamma passing rates (GPRs) ≥95% for a 3%/1 mm criteria. For the 18 multiple-target SRS plans, average 3D GPRs greater than 93% and 98% for the 3%/2 mm and 5%/2 mm criteria, respectively. Our results validate the use of VIPER as a dosimetric verification tool for pre-treatment QA of single-isocenter multiple-target SRS plans. The method requires no setup time on the linac and results in an accurate 3D characterization of the delivered dose.

Eigencolor radiochromic film dosimetry

PURPOSE

The goal of this work is to propose a new multichannel method correcting for systematic thickness disturbances and to evaluate its precision in relevant radiation dosimetry applications.

METHODS

The eigencolor ratio technique is introduced and theoretically developed to provide a method correcting for thickness disturbances. The method is applied to EBT3 Gafchromic^TM film irradiated with cobalt-60 and 6 MV photon beams and digitized with an Epson 10000XL photo scanner. Dose profiles and output factors of different field sizes are measured and analyzed. Variance analysis of the previous method of Bouchard et al. ["On the characterization and uncertainty analysis of radiochromic film dosimetry" Med Phys. 2009;36:1931-1946] is adapted to the new approach. Uncertainties are predicted for relevant applications.

RESULTS

Results show that systematic disturbances attributed to thickness variations are efficiently corrected. The method is shown efficient to identify and correct for dark spots which cause systematic errors in single-channel distributions. Applications of the method in the context of relative dosimetry yields standard uncertainties ranging between 0.8% and 1.9%, depending on the region of interest (ROI) size and the film irradiation. Variance analysis predicts that uncertainty levels between 0.3% and 0.6% are achievable with repeated measurements. Uncertainties are found to vary with absorbed dose and ROI size.

CONCLUSIONS

The proposed multichannel method is efficient for accurate dosimetry, reaching uncertainty levels comparable to previous publications with EBT film. The method is also promising for applications beyond clinical QA, such as machine characterization and other advanced dosimetry applications.

A protocol for absolute dose verification of SBRT/SRS treatment plans using Gafchromic™ EBT-XD films

PURPOSE

To provide a practical protocol for absolute dose verification of stereotactic body radiotherapy (SBRT) and stereotactic radiosurgery (SRS) treatment plans, based on our clinical experience. It aims to be a concise summary of the main aspects to be considered when establishing an accurate film dosimetry system.

METHODS

Procedures for film calibration and conversion to dose are described for a dosimetry system composed of Gafchromic™ EBT-XD films and a flatbed document scanner. Factors that affect the film-scanner response are also reviewed and accounted for. The accuracy of the proposed methodology was assessed by taking a set of strips irradiated to known doses and its applicability is illustrated for ten SBRT/SRS treatment plans. The film response was converted to dose using red and triple channel dosimetry. The agreement between the planned and measured dose distributions was evaluated using global gamma analysis with criteria of 3%/2mm 10% threshold (TH), 2%/2mm 10% TH, and 2%/2mm 20% TH.

RESULTS

The differences between the expected and determined doses from the strips analysis were 0.9 ± 0.6% for the red channel and 1.1 ± 0.7% for the triple channel method. Regarding the SBRT/SRS plans verification, the mean gamma passing rates were 99.5 ± 1.0% vs 99.6 ± 1.0% (3%/2mm 10% TH), 96.9 ± 3.5% vs 99.1 ± 1.3% (2%/2mm 10% TH) and 98.4 ± 1.8% vs 98.8 ± 1.5% (2%/2mm 20% TH) for red and triple channel dosimetry, respectively.

CONCLUSIONS

The proposed protocol allows for accurate absolute dose verification of SBRT/SRS treatment plans, applying both single and triple channel methods. It may work as a guide for users that intend to implement a film dosimetry system.

In vivo dosimetry in low-voltage IORT breast treatments with XR-RV3 radiochromic film

PURPOSE

The objectives of the study were to establish a procedure for in vivo film-based dosimetry for intraoperative radiotherapy (IORT), evaluate the typical doses delivered to organs at risk, and verify the dose prescription.

MATERIALS AND METHODS

In vivo dose measurements were studied using XR-RV3 radiochromic films in 30 patients with breast cancer undergoing IORT using the Axxent® device (Xoft Inc.). The stability of the radiochromic films in the energy ranges used was verified by taking measurements at different depths. The stability of the scanner response was tested, and 5 different calibration curves were constructed for different beam qualities. Six pieces of film were placed in each of the 30 patients. All the pieces were correctly sterilized and checked to ensure that the process did not affect the outcome. All calibration and dose measurements were analyzed using the Radiochromic.com software application.

RESULTS

The doses were measured for 30 patients. The doses in contact with the applicator (prescription zone) were 19.8 ± 0.9 Gy. In the skin areas, the doses were as follows: 1-2 cm from the applicator, 1.86 ± 0.77 Gy; 2-5 cm, 0.73 ± 0.14 Gy; and greater than 5 cm, 0.28 ± 0.17 Gy. The dose delivered to the pectoral muscle (tungsten shielding disc) was 0.51 ± 0.27 Gy.

CONCLUSIONS

The study demonstrated the viability of XR-RV3 films for in vivo dose measurement in the dose and energy ranges applied in a complex procedure, such as breast IORT. The doses in organs at risk were far below the tolerances for cases such as those studied.

A new approach to radiochromic film dosimetry based on non-local means

Radiochromic film in conjunction with flatbed scanners are frequently employed as dosimeters for advanced techniques in radiotherapy. Their strengths are as follows: light element composition, low energy dependence, near biological tissue equivalence and high spatial resolution. However, they have some weaknesses as well: non-uniformities, read out noise, and scanning artifacts. Several processing protocols have been proposed intending to correct the perturbations these weaknesses produce. The aim of this paper is to present a new processing protocol for radiochromic film dosimetry based on a non-local means denoising algorithm. Three dose distributions of open square fields and a spatial combination of these fields using different angles of incidence and monitor units have been employed to validate the protocol. The dose distributions are traceable to ionization chamber measurements. Additionally, a real dose distribution of a treatment was used to simulate scanning data with noise and scanning lateral artifact, and to study how the protocol behaves under these perturbations. The same measured raw data have been processed by means of an implementation of the multichannel protocol (multigaussian method). It has been found that the proposed protocol reduces dose uncertainty even though it uses fewer scans than the multichannel protocol.

Comparison of three film analysis softwares using EBT2 and EBT3 films in radiotherapy

Introduction The purpose of the study was to compare the results of gamma value based film analysis according to the used type of self-developer film and software product. Material and methods The films were irradiated with different treatment techniques such as 3D conformal and intensity modulated radiotherapy with static and rotational delivery. Stereotactic plans with conformal and intensity modulated arc techniques, using coplanar and non-coplanar beam setup were also evaluated. The data of irradiated film were compared with the planned planar dose distribution exported from the treatment planning system. Three film analysis software programs were evaluated: PTW Mephysto (PTW), FilmQA Pro (FQP) and radiohromic.com(RC). Both EBT2 and EBT3 types of films were examined. The comparisons of dose distributions were performed with gamma analysis using 10% cut-off level. Results The results of the gamma analysis for larger fields were between 78.3% and 98.3%, 75.7% and 100%, 80.2% and 98.8% with PTW, FQP and RC, respectively. The results of evaluation in case of stereotactic measurements were 76.8%-99.2% for PTW, 95.7%-100% for FQP and 91.2%-99.9% for RC. Conclusions All the three software programs are suitable for calibrating and evaluating films, performing gamma analysis, and can be used for patient specific quality assurance measurements. There is no direct connection between gamma passing rate and absolute accuracy or software quality, it is just a feature of the software. The interpretation of own results has to be defined on an institutional level according to given workflow and preliminary results.

Output correction factors for small static fields in megavoltage photon beams for seven ionization chambers in two orientations - perpendicular and parallel

Purpose

The goal of the present work was to provide a large set of detector‐specific output correction factors for seven small volume ionization chambers on two linear accelerators in four megavoltage photon beams utilizing perpendicular and parallel orientation of ionization chambers in the beam for nominal field sizes ranging from 0.5 cm² × 0.5 cm² to 10 cm² × 10 cm². The present study is the second part of an extensive research conducted by our group.

Methods

Output correction factors kQclin,Qreffclin,fref were experimentally determined on two linacs, Elekta Versa HD and Varian TrueBeam for 6 and 10 MV beams with and without flattening filter for nine square fields ranging from 0.5 cm² × 0.5 cm² to 10 cm² × 10 cm², for seven mini and micro ionization chambers, IBA CC04, IBA Razor, PTW 31016 3D PinPoint, PTW 31021 3D Semiflex, PTW 31022 3D PinPoint, PTW 31023 PinPoint, and SI Exradin A16. An Exradin W1 plastic scintillator and EBT3 radiochromic films were used as the reference detectors.

Results

For all ionization chambers, values of output correction factors kQclin,Qreffclin,fref were lower for parallel orientation compared to those obtained in the perpendicular orientation. Five ionization chambers from our study set, IBA Razor, PTW 31016 3D PinPoint, PTW 31022 3D PinPoint, PTW 31023 PinPoint, and SI Exradin A16, fulfill the requirement recommended in the TRS‐483 Code of Practice, that is, 0.95<kQclin,Qreffclin,fref<1.05, down to the field size 0.8 cm² × 0.8 cm², when they are positioned in parallel orientation; two of the ionization chambers, IBA Razor and PTW 31023 PinPoint, satisfy this condition down to the field size of 0.5 cm² × 0.5 cm².

Conclusions

The present paper provides experimental results of detector‐specific output correction factors for seven small volume ionization chambers. Output correction factors were determined in 6 and 10 MV photon beams with and without flattening filter down to the square field size of 0.5 cm² × 0.5 cm² for two orientations of ionization chambers — perpendicular and parallel. Our main finding is that output correction factors are smaller if they are determined in a parallel orientation compared to those obtained in a perpendicular orientation for all ionization chambers regardless of the photon beam energy, filtration, or linear accelerator being used. Based on our findings, we recommend using ionization chambers in parallel orientation, to minimize corrections in the experimental determination of field output factors. Latter holds even for field sizes below 1.0 cm² × 1.0 cm², whenever necessary corrections remain within 5%, which was the case for several ionization chambers from our set.

TRS‐483 recommended perpendicular orientation of ionization chambers for the determination of field output factors. The present study presents results for both perpendicular and parallel orientation of ionization chambers. When validated by other researchers, the present results for parallel orientation can be considered as a complementary dataset to those given in TRS‐483.

A novel method for the determination of field output factors and output correction factors for small static fields for six diodes and a microdiamond detector in megavoltage photon beams

Purpose

The goal of this work is to provide a large and consistent set of data for detector‐specific output correction factors, kQclin,Qreffclin,fref, for small static fields for seven solid‐state detectors and to determine field output factors, ΩQclin,Qreffclin,fref, using EBT3 radiochromic films and W1 plastic scintillator as reference detectors on two different linear accelerators and four megavoltage photon beams. Consistent measurement conditions and recommendations given in the International Code of Practice TRS‐483 for small‐field dosimetry were followed throughout the study.

Methods

ΩQclin,Qreffclin,fref were determined on two linacs, Elekta Versa HD and Varian TrueBeam, for 6 and 10 MV beams with and without flattening filter and for nine fields ranging from 0.5 × 0.5 cm² to 10 × 10 cm². Signal readings obtained with EBT3 radiochromic films and W1 plastic scintillator were fitted by an analytical function. Volume averaging correction factors, determined from two‐dimensional (2D) dose matrices obtained with EBT3 films and fitted to bivariate Gaussian function, were used to correct measured signals. kQclin,Qreffclin,fref were determined empirically for six diodes, IBA SFD, IBA Razor, PTW 60008 P, PTW 60012 E, PTW 60018 SRS, and SN EDGE, and a PTW 60019 microDiamond detector.

Results

Field output factors and detector‐specific kQclin,Qreffclin,fref are presented in the form of analytical functions as well as in the form of discrete values. It is found that in general, for a given linac, small‐field output factors need to be determined for every combination of beam energy and filtration (WFF or FFF) and field size as the differences between them can be statistically significant (P < 0.05). For different beam energies, the present data for kQclin,Qreffclin,fref are found to differ significantly (P < 0.05) from the corresponding data published in TRS‐483 mostly for the smallest fields (<1.5 cm). For the PTW microDiamond detector, statistically significant differences (P < 0.05) between kQclin,Qreffclin,fref values were found for all investigated beams on an Elekta Versa HD linac for field sizes 0.5 × 0.5 cm² and 0.8 × 0.8 cm². Significant differences in kQclin,Qreffclin,fref between beams of a given energy but with and without flattening filters are found for measurements made in small fields (<1.5 cm) at a given linac. Differences in kQclin,Qreffclin,fref are also found when measurements are made at different linacs using the same beam energy filtration combination; for the PTW microDiamond detector, these differences were found to be around 6% and were considered as significant.

Conclusions

Selection of two reference detectors, EBT3 films and W1 plastic scintillator, and use of an analytical function, is a novel approach for the determination of ΩQclin,Qreffclin,fref for small static fields in megavoltage photon beams. Large set of kQclin,Qreffclin,fref data for seven solid‐state detectors and four beam energies determined on two linacs by a single group of researchers can be considered a valuable supplement to the literature and the TRS‐483 dataset.