Assessment of bias in scoring of AI-based radiotherapy segmentation and planning studies using modified TRIPOD and PROBAST guidelines as an example.

BACKGROUND AND PURPOSE: Studies investigating the application of Artificial Intelligence (AI) in the field of radiotherapy exhibit substantial variations in terms of quality. The goal of this study was to assess the amount of transparency and bias in scoring articles with a specific focus on AI based segmentation and treatment planning, using modified PROBAST and TRIPOD checklists, in order to provide recommendations for future guideline developers and reviewers. MATERIALS AND METHODS: The TRIPOD and PROBAST checklist items were discussed and modified using a Delphi process. After consensus was reached, 2 groups of 3 co-authors scored 2 articles to evaluate usability and further optimize the adapted checklists. Finally, 10 articles were scored by all co-authors. Fleiss' kappa was calculated to assess the reliability of agreement between observers. RESULTS: Three of the 37 TRIPOD items and 5 of the 32 PROBAST items were deemed irrelevant. General terminology in the items (e.g., multivariable prediction model, predictors) was modified to align with AI-specific terms. After the first scoring round, further improvements of the items were formulated, e.g., by preventing the use of sub-questions or subjective words and adding clarifications on how to score an item. Using the final consensus list to score the 10 articles, only 2 out of the 61 items resulted in a statistically significant kappa of 0.4 or more demonstrating substantial agreement. For 41 items no statistically significant kappa was obtained indicating that the level of agreement among multiple observers is due to chance alone. CONCLUSION: Our study showed low reliability scores with the adapted TRIPOD and PROBAST checklists. Although such checklists have shown great value during development and reporting, this raises concerns about the applicability of such checklists to objectively score scientific articles for AI applications. When developing or revising guidelines, it is essential to consider their applicability to score articles without introducing bias.

Stereotactic Arrhythmia Radioablation (STAR): Assessment of cardiac and respiratory heart motion in ventricular tachycardia patients - A STOPSTORM.eu consortium review.

AIM: To identify the optimal STereotactic Arrhythmia Radioablation (STAR) strategy for individual patients, cardiorespiratory motion of the target volume in combination with different treatment methodologies needs to be evaluated. However, an authoritative overview of the amount of cardiorespiratory motion in ventricular tachycardia (VT) patients is missing. METHODS: In this STOPSTORM consortium study, we performed a literature review to gain insight into cardiorespiratory motion of target volumes for STAR. Motion data and target volumes were extracted and summarized. RESULTS: Out of the 232 studies screened, 56 provided data on cardiorespiratory motion, of which 8 provided motion amplitudes in VT patients (n = 94) and 10 described (cardiac/cardiorespiratory) internal target volumes (ITVs) obtained in VT patients (n = 59). Average cardiac motion of target volumes was < 5 mm in all directions, with maximum values of 8.0, 5.2 and 6.5 mm in Superior-Inferior (SI), Left-Right (LR), Anterior-Posterior (AP) direction, respectively. Cardiorespiratory motion of cardiac (sub)structures showed average motion between 5-8 mm in the SI direction, whereas, LR and AP motions were comparable to the cardiac motion of the target volumes. Cardiorespiratory ITVs were on average 120-284% of the gross target volume. Healthy subjects showed average cardiorespiratory motion of 10-17 mm in SI and 2.4-7 mm in the AP direction. CONCLUSION: This review suggests that despite growing numbers of patients being treated, detailed data on cardiorespiratory motion for STAR is still limited. Moreover, data comparison between studies is difficult due to inconsistency in parameters reported. Cardiorespiratory motion is highly patient-specific even under motion-compensation techniques. Therefore, individual motion management strategies during imaging, planning, and treatment for STAR are highly recommended.

Treatment plan quality assessment for radiotherapy of rectal cancer patients using prediction of organ-at-risk dose metrics.

BACKGROUND AND PURPOSE: Radiotherapy centers frequently lack simple tools for periodic treatment plan verification and feedback on current plan quality. It is difficult to measure treatment quality over different years or during the planning process. Here, we implemented plan quality assurance (QA) by developing a database of dose-volume histogram (DVH) metrics and a prediction model. These tools were used to assess automatically optimized treatment plans for rectal cancer patients, based on cohort analysis. MATERIAL AND METHODS: A treatment plan QA framework was established and an overlap volume histogram based model was used to predict DVH parameters for cohorts of patients treated in 2018 and 2019 and grouped according to planning technique. A training cohort of 22 re-optimized treatment plans was used to make the prediction model. The prediction model was validated on 95 automatically generated treatment plans (automatically optimized cohort) and 93 manually optimized plans (manually optimized cohort). RESULTS: For the manually optimized cohort, on average the prediction deviated less than 0.3 ± 1.4 Gy and -4.3 ± 5.5 Gy, for the mean doses to the bowel bag and bladder, respectively; for the automatically optimized cohort a smaller deviation was observed: -0.1 ± 1.1 Gy and -0.2 ± 2.5 Gy, respectively. The interquartile range of DVH parameters was on average smaller for the automatically optimized cohort, indicating less variation within each parameter compared to manual planning. CONCLUSION: An automated framework to monitor treatment quality with a DVH prediction model was successfully implemented clinically and revealed less variation in DVH parameters for automated in comparison to manually optimized plans. The framework also allowed for individual feedback and DVH estimation.

Challenges and caveats of a multi-center retrospective radiomics study: an example of early treatment response assessment for NSCLC patients using FDG-PET/CT radiomics.

BACKGROUND: Prognostic models based on individual patient characteristics can improve treatment decisions and outcome in the future. In many (radiomic) studies, small size and heterogeneity of datasets is a challenge that often limits performance and potential clinical applicability of these models. The current study is example of a retrospective multi-centric study with challenges and caveats. To highlight common issues and emphasize potential pitfalls, we aimed for an extensive analysis of these multi-center pre-treatment datasets, with an additional 18F-fluorodeoxyglucose (FDG) positron emission tomography/computed tomography (PET/CT) scan acquired during treatment. METHODS: The dataset consisted of 138 stage II-IV non-small cell lung cancer (NSCLC) patients from four different cohorts acquired from three different institutes. The differences between the cohorts were compared in terms of clinical characteristics and using the so-called 'cohort differences model' approach. Moreover, the potential prognostic performances for overall survival of radiomic features extracted from CT or FDG-PET, or relative or absolute differences between the scans at the two time points, were assessed using the LASSO regression method. Furthermore, the performances of five different classifiers were evaluated for all image sets. RESULTS: The individual cohorts substantially differed in terms of patient characteristics. Moreover, the cohort differences model indicated statistically significant differences between the cohorts. Neither LASSO nor any of the tested classifiers resulted in a clinical relevant prognostic model that could be validated on the available datasets. CONCLUSION: The results imply that the study might have been influenced by a limited sample size, heterogeneous patient characteristics, and inconsistent imaging parameters. No prognostic performance of FDG-PET or CT based radiomics models can be reported. This study highlights the necessity of extensive evaluations of cohorts and of validation datasets, especially in retrospective multi-centric datasets.

Development of a virtual spacer to support the decision for the placement of an implantable rectum spacer for prostate cancer radiotherapy: Comparison of dose, toxicity and cost-effectiveness.

INTRODUCTION: Previous studies have shown that the implantable rectum spacer (IRS) is not beneficial for all patients. A virtual IRS (V-IRS) was constructed to help identify the patients for whom it is cost-effective to implant an IRS, and its viability as a tool to tailor the decision of an IRS implantation to be beneficial for the specified patient was assessed. Please watch animation: (https://www.youtube.com/watch?v=tDlagSXMKqw) MATERIALS AND METHODS: The V-IRS was tested on 16 patients: 8 with a rectal balloon implant (RBI) and 8 with a hydrogel spacer. A V-IRS was developed using 7 computed tomography (CT) scans of patients with a RBI. To examine the V-IRS, CT scans before and after the implantation of an IRS were used. IMRT plans were made based on CT scans before the IRS, after IRS and with the V-IRS, prescribing 70 Gray (Gy) to the planning target volume. Toxicity was accessed using externally validated normal tissue complication probability (NTCP) models, and the Cost-effectiveness was analyzed using a published Markov model. RESULTS: The rectum volume receiving 75Gy (V75) were improved by both the IRS and the V-IRS with on average 4.2% and 4.3% respectively. The largest NTCP reduction resulting from the IRS and the V-IRS was 4.0% and 3.9% respectively. The RBI was cost-effective for 1 out of 8 patients, and the hydrogel was effective for 2 out of 8 patients, and close to effective for a third patient. The classification accuracy of the model, regarding cost-effectiveness, was 100%. CONCLUSION: The V-IRS approach in combination with a toxicity prediction model and a cost-effectiveness analyses is a promising basis for a decision support tool for the implantation of either a hydrogel spacer or a rectum balloon implant.

Quality assessment of positron emission tomography scans: recommendations for future multicentre trials.

BACKGROUND: Standardization protocols and guidelines for positron emission tomography (PET) in multicenter trials are available, despite a large variability in image acquisition and reconstruction parameters exist. In this study, we investigated the compliance of PET scans to the guidelines of the European Association of Nuclear Medicine (EANM). From these results, we provide recommendations for future multicenter studies using PET. MATERIAL AND METHODS: Patients included in a multicenter randomized phase II study had repeated PET scans for early response assessment. Relevant acquisition and reconstruction parameters were extracted from the digital imaging and communications in medicine (DICOM) header of the images. The PET image parameters were compared to the guidelines of the EANM for tumor imaging version 1.0 recommended parameters. RESULTS: From the 223 included patients, 167 baseline scans and 118 response scans were available from 15 hospitals. Scans of 19% of the patients had an uptake time that fulfilled the Uniform Protocols for Imaging in Clinical Trials response assessment criteria. The average quality score over all hospitals was 69%. Scans with a non-compliant uptake time had a larger standard deviation of the mean standardized uptake value (SUVmean) of the liver than scans with compliant uptake times. CONCLUSIONS: Although a standardization protocol was agreed on, there was a large variability in imaging parameters. For future, multicenter studies including PET imaging a prospective central quality review during patient inclusion is needed to improve compliance with image standardization protocols as defined by EANM.

Quantitative assessment of Zirconium-89 labeled cetuximab using PET/CT imaging in patients with advanced head and neck cancer: a theragnostic approach.

Biomarkers predicting treatment response to the monoclonal antibody cetuximab in locally advanced head and neck squamous cell carcinomas (LAHNSCC) are lacking. We hypothesize that tumor accessibility is an important factor in treatment success of the EGFR targeting drug. We quantified uptake of cetuximab labeled with Zirconium-89 (89Zr) using PET/CT imaging.Seventeen patients with stage III-IV LAHNSCC received a loading dose unlabeled cetuximab, followed by 10 mg 54.5±9.6 MBq 89Zr-cetuximab. PET/CT images were acquired either 3 and 6 or 4 and 7 days post-injection. 89Zr-cetuximab uptake was quantified using standardized uptake value (SUV) and tumor-to-background ratio (TBR), and correlated to EGFR immunohistochemistry. TBR was compared between scan days to determine optimal timing.Uptake of 89Zr-cetuximab varied between patients (day 6-7: SUVpeak range 2.5-6.2). TBR increased significantly (49±28%, p < 0.01) between first (1.1±0.3) and second scan (1.7±0.6). Between groups with a low and high EGFR expression a significant difference in SUVmean (2.1 versus 3.0) and SUVpeak (3.2 versus 4.7) was found, however, not in TBR. Data is available at www.cancerdata.org (DOI: 10.17195/candat.2016.11.1).In conclusion, 89Zr-cetuximab PET imaging shows large inter-patient variety in LAHNSCC and provides additional information over FDG-PET and EGFR expression. Validation of the predictive value is recommended with scans acquired 6-7 days post-injection.

[18F]FDG PET/CT-based response assessment of stage IV non-small cell lung cancer treated with paclitaxel-carboplatin-bevacizumab with or without nitroglycerin patches.

PURPOSE: Nitroglycerin (NTG) is a vasodilating drug, which increases tumor blood flow and consequently decreases hypoxia. Therefore, changes in [18F] fluorodeoxyglucose positron emission tomography ([18F]FDG PET) uptake pattern may occur. In this analysis, we investigated the feasibility of [18F]FDG PET for response assessment to paclitaxel-carboplatin-bevacizumab (PCB) treatment with and without NTG patches. And we compared the [18F]FDG PET response assessment to RECIST response assessment and survival. METHODS: A total of 223 stage IV non-small cell lung cancer (NSCLC) patients were included in a phase II study (NCT01171170) randomizing between PCB treatment with or without NTG patches. For 60 participating patients, a baseline and a second [18F]FDG PET/computed tomography (CT) scan, performed between day 22 and 24 after the start of treatment, were available. Tumor response was defined as a 30 % decrease in CT and PET parameters, and was compared to RECIST response at week 6. The predictive value of these assessments for progression free survival (PFS) and overall survival (OS) was assessed with and without NTG. RESULTS: A 30 % decrease in SUVpeak assessment identified more patients as responders compared to a 30 % decrease in CT diameter assessment (73 % vs. 18 %), however, this was not correlated to OS (SUVpeak30 p = 0.833; CTdiameter30 p = 0.557). Changes in PET parameters between the baseline and the second scan were not significantly different for the NTG group compared to the control group (p value range 0.159-0.634). The CT-based (part of the [18F]FDG PET/CT) parameters showed a significant difference between the baseline and the second scan for the NTG group compared to the control group (CT diameter decrease of 7 ± 23 % vs. 19 ± 14 %, p = 0.016, respectively). CONCLUSIONS: The decrease in tumoral FDG uptake in advanced NSCLC patients treated with chemotherapy with and without NTG did not differ between both treatment arms. Early PET-based response assessment showed more tumor responders than CT-based response assessment (part of the [18F]FDG PET/CT); this was not correlated to survival. This might be due to timing of the [18F]FDG PET shortly after the bevacizumab infusion.

Imaging for Target Volume Definition and Response Assessment in Lung Cancer.

Target volume definition is of obvious importance in successful radiotherapy. Single-energy CT scans remain the standard, but FDG-PET-CT scans aid the determination of which lymph nodes should be included in the gross tumor volume and to fine-tune areas of cancer involvement. FDG-PET-CT imaging remains the gold standard in clinical practice. Hypoxia and proliferation tracers are still investigational, as is PET-guided redistribution of the radiation dose within the tumor. Contrast-enhanced CT as well as 4D CT scans contain information such as the characteristics of the lungs that are related to individual radiosensitivity, ventilation, and perfusion. Dual-energy CT imaging holds promise for the future for characterization of both tumor and normal tissues. The assessment of response after radiotherapy on the basis of CT scans remains difficult because of inflammatory and fibrotic changes. RECIST is still the standard. FDG avidity suffers from too high rates of false positive and false negative signals and is therefore not recommended, except on clinical indication.

Assessment of tumour size in PET/CT lung cancer studies: PET- and CT-based methods compared to pathology.

BACKGROUND: Positron emission tomography (PET) may be useful for defining the gross tumour volume for radiation treatment planning and for response monitoring of non-small cell lung cancer (NSCLC) patients. The purpose of this study was to compare tumour sizes obtained from CT- and various more commonly available PET-based tumour delineation methods to pathology findings. METHODS: Retrospective non-respiratory gated whole body [18F]-fluoro-2-deoxy-D-glucose PET/CT studies from 19 NSCLC patients were used. Several (semi-)automatic PET-based tumour delineation methods and manual CT-based delineation were used to assess the maximum tumour diameter. RESULTS: 50%, adaptive 41% threshold-based and contrast-oriented delineation methods showed good agreement with pathology after removing two outliers (R2=0.82). An absolute SUV threshold of 2.5 also showed a good agreement with pathology after the removal of 5 outliers (R2: 0.79), but showed a significant overestimation in the maximum diameter (19.8 mm, p<0.05). Adaptive 50%, relative threshold level and gradient-based methods did not show any outliers, provided only small, non-significant differences in maximum tumour diameter (<4.7 mm, p>0.10), and showed fair correlation (R2>0.62) with pathology. Although adaptive 70% threshold-based methods showed underestimation compared to pathology (36%), it provided the best precision (SD: 14%) together with good correlation (R2=0.81). Good correlation between CT delineation and pathology was observed (R2=0.77). However, CT delineation showed a significant overestimation compared with pathology (3.8 mm, p<0.05). CONCLUSIONS: PET-based tumour delineation methods provided tumour sizes in agreement with pathology and may therefore be useful to define the (metabolically most) active part of the tumour for radiotherapy and response monitoring purposes.

Response assessment using 18F-FDG PET early in the course of radiotherapy correlates with survival in advanced-stage non-small cell lung cancer.

UNLABELLED: This study investigated the possibility of early response assessment based on (18)F-FDG uptake during radiotherapy with respect to overall survival in patients with non-small cell lung cancer. METHODS: (18)F-FDG PET/CT was performed before radiotherapy and was repeated in the second week of radiotherapy for 34 consecutive lung cancer patients. The CT volume and standardized uptake value (SUV) parameters of the primary tumor were quantified at both time points. Changes in volume and SUV parameters correlated with 2-y overall survival. RESULTS: The average change in mean SUV in the primary tumor of patients with a 2-y survival was a decrease by 20% ± 21%-significantly different (P < 0.007) from nonsurvivors, who had an increase by 2% ± 22%. A sensitivity and specificity of 63% and 93%, respectively, to separate the 2 groups was reached for a decrease in mean SUV of 15%. Survival curves were significantly different using this cutoff (P = 0.001). The hazard ratio for a 1% decrease in mean SUV was 1.032 (95% confidence interval, 1.010-1.055). Changes in tumor volume defined on CT did not correlate with overall survival. CONCLUSION: The use of repeated (18)F-FDG PET to assess treatment response early during radiotherapy is possible in patients undergoing radiotherapy or sequential or concurrent chemoradiotherapy. A decrease in (18)F-FDG uptake by the primary tumor correlates with higher long-term overall survival.

3D dose delivery verification using repeated cone-beam imaging and EPID dosimetry for stereotactic body radiotherapy of non-small cell lung cancer.

PURPOSE: To implement a 3D dose verification procedure, based on in-room cone-beam CT imaging and portal dosimetry, for lung cancer patients treated with stereotactic body radiotherapy (SBRT). MATERIALS AND METHODS: MV cone-beam CT scans were made for patient positioning and calibrated for dose calculation purposes. Prior to treatment, the treatment fields were captured using a calibrated electronic portal imaging device (EPID). A Monte Carlo dose reconstruction model was used to estimate the 3D dose delivered to the patient inside the cone-beam CT images. The planned and delivered dose distributions were compared for 4 patients and 10 treatment fractions using dose-volume histograms and gamma analysis. RESULTS: The gamma analysis showed a good agreement between the planned and delivered dose distributions for patients without changes in anatomy. The delivered mean dose per fraction inside the target volume deviated on average 1.1+/-1.4% from the planned dose. For the critical organs, only minor differences were observed between the reconstructed and planned dose. CONCLUSIONS: A method was presented that allows verification of the dose delivered in 3D for lung cancer patients treated with SBRT. The procedure is independent of the treatment planning system and uses in-room MV cone-beam CT imaging and portal dosimetry.

3D in vivo dosimetry using megavoltage cone-beam CT and EPID dosimetry.

PURPOSE: To develop a method that reconstructs, independently of previous (planning) information, the dose delivered to patients by combining in-room imaging with transit dose measurements during treatment. METHODS AND MATERIALS: A megavoltage cone-beam CT scan of the patient anatomy was acquired with the patient in treatment position. During treatment, delivered fields were measured behind the patient with an electronic portal imaging device. The dose information in these images was back-projected through the cone-beam CT scan and used for Monte Carlo simulation of the dose distribution inside the cone-beam CT scan. Validation was performed using various phantoms for conformal and IMRT plans. Clinical applicability is shown for a head-and-neck cancer patient treated with IMRT. RESULTS: For single IMRT beams and a seven-field IMRT step-and-shoot plan, the dose distribution was reconstructed within 3%/3mm compared with the measured or planned dose. A three-dimensional conformal plan, verified using eight point-dose measurements, resulted in a difference of 1.3 +/- 3.3% (1 SD) compared with the reconstructed dose. For the patient case, planned and reconstructed dose distribution was within 3%/3mm for about 95% of the points within the 20% isodose line. Reconstructed mean dose values, obtained from dose-volume histograms, were within 3% of prescribed values for target volumes and normal tissues. CONCLUSIONS: We present a new method that verifies the dose delivered to a patient by combining in-room imaging with the transit dose measured during treatment. This verification procedure opens possibilities for offline adaptive radiotherapy and dose-guided radiotherapy strategies taking into account the dose distribution delivered during treatment sessions.

The next step in patient-specific QA: 3D dose verification of conformal and intensity-modulated RT based on EPID dosimetry and Monte Carlo dose calculations.

BACKGROUND AND PURPOSE: A method was evaluated to reconstruct the 3D dose distribution in patients using their planning CT-scan in combination with a Monte Carlo calculation, and the energy fluence of the actual treatment beams measured pre-treatment with an EPID without the patient or a phantom in the beam. MATERIALS AND METHODS: Nine plans of lung cancer patients treated with a 3D conformal technique, calculated using a simple convolution algorithm (CA), as well as five IMRT treatments of head-and-neck cancer patients, calculated with a more advanced superposition algorithm (SA), were verified. Differences between planned and reconstructed dose distributions were quantified in terms of DVH parameters. RESULTS: For the lung cancer group, differences between the reconstructed mean PTV dose and the values calculated with the TPS were 5.0+/-4.2% (1SD) and -1.4+/-1.5% for the CA and SA algorithm, respectively. No large differences in the lung and spinal cord DVH parameters were found. For the IMRT treatments, the average dose differences in the PTV were generally below 3%. The reconstructed mean parotid gland dose was 3.2+/-1.2% lower, while the maximum spinal cord dose was on average 3.1+/-1.9% higher. CONCLUSIONS: EPID dosimetry combined with 3D dose reconstruction is a useful procedure for patient-specific QA of complex treatments. DVH parameters can be used to interpret the dose distribution delivered to the patient in the same way as during standard treatment plan evaluation.

Treatment verification in the presence of inhomogeneities using EPID-based three-dimensional dose reconstruction.

Treatment verification is a prerequisite for the verification of complex treatments, checking both the treatment planning process and the actual beam delivery. Pretreatment verification can detect errors introduced by the treatment planning system (TPS) or differences between planned and delivered dose distributions. In a previous paper we described the reconstruction of three-dimensional (3-D) dose distributions in homogeneous phantoms using an in-house developed model based on the beams delivered by the linear accelerator measured with an amorphous silicon electronic portal imaging device (EPID), and a dose calculation engine using the Monte Carlo code XVMC. The aim of the present study is to extend the method to situations in which tissue inhomogeneities are present and to make a comparison with the dose distributions calculated by the TPS. Dose distributions in inhomogeneous phantoms, calculated using the fast-Fourier transform convolution (FFTC) and multigrid superposition (MGS) algorithms present in the TPS, were verified using the EPID-based dose reconstruction method and compared to film and ionization chamber measurements. Differences between dose distributions were evaluated using the gamma-evaluation method (3%/3 mm) and expressed as a mean gamma and the percentage of points with gamma> 1 (P(gamma>1)). For rectangular inhomogeneous phantoms containing a low-density region, the differences between film and reconstructed dose distributions were smaller than 3%. In low-density regions there was an overestimation of the planned dose using the FFTC and MGS algorithms of the TPS up to 20% and 8%, respectively, for a 10 MV photon beam and a 3 x 3 cm2 field. For lower energies and larger fields (6 MV, 5 x 5 cm2), these differences reduced to 6% and 3%, respectively. Dose reconstruction performed in an anthropomorphic thoracic phantom for a 3-D conformal and an IMRT plan, showed good agreement between film data and reconstructed dose values (P(gamma>1) <6%). The algorithms of the TPS underestimated the dose in the low-dose regions outside the treatment field, due to an implementation error of the jaws and multileaf collimator of the linac in the TPS. The FFTC algorithm of the TPS showed differences up to 6% or 6 mm at the interface between lung and breast. Two intensity-modulated radiation therapy head and neck plans, reconstructed in a commercial phantom having a bone-equivalent insert and an air cavity, showed good agreement between film measurement, reconstructed and planned dose distributions using the FFTC and MGS algorithm, except in the bone-equivalent regions where both TPS algorithms underestimated the dose with 4%. Absolute dose verification was performed at the isocenter where both planned and reconstructed dose were within 2% of the measured dose. Reproducibility for the EPID measurements was assessed and found to be of negligible influence on the reconstructed dose distribution. Our 3-D dose verification approach is based on the actual dose measured with an EPID in combination with a Monte Carlo dose engine, and therefore independent of a TPS. Because dose values are reconstructed in 3-D, isodose surfaces and dose-volume histograms can be used to detect dose differences in target volume and normal tissues. Using our method, the combined planning and treatment delivery process is verified, offering an easy to use tool for the verification of complex treatments.

A Monte Carlo based three-dimensional dose reconstruction method derived from portal dose images.

The verification of intensity-modulated radiation therapy (IMRT) is necessary for adequate quality control of the treatment. Pretreatment verification may trace the possible differences between the planned dose and the actual dose delivered to the patient. To estimate the impact of differences between planned and delivered photon beams, a three-dimensional (3-D) dose verification method has been developed that reconstructs the dose inside a phantom. The pretreatment procedure is based on portal dose images measured with an electronic portal imaging device (EPID) of the separate beams, without the phantom in the beam and a 3-D dose calculation engine based on the Monte Carlo calculation. Measured gray scale portal images are converted into portal dose images. From these images the lateral scattered dose in the EPID is subtracted and the image is converted into energy fluence. Subsequently, a phase-space distribution is sampled from the energy fluence and a 3-D dose calculation in a phantom is started based on a Monte Carlo dose engine. The reconstruction model is compared to film and ionization chamber measurements for various field sizes. The reconstruction algorithm is also tested for an IMRT plan using 10 MV photons delivered to a phantom and measured using films at several depths in the phantom. Depth dose curves for both 6 and 10 MV photons are reconstructed with a maximum error generally smaller than 1% at depths larger than the buildup region, and smaller than 2% for the off-axis profiles, excluding the penumbra region. The absolute dose values are reconstructed to within 1.5% for square field sizes ranging from 5 to 20 cm width. For the IMRT plan, the dose was reconstructed and compared to the dose distribution with film using the gamma evaluation, with a 3% and 3 mm criterion. 99% of the pixels inside the irradiated field had a gamma value smaller than one. The absolute dose at the isocenter agreed to within 1% with the dose measured with an ionization chamber. It can be concluded that our new dose reconstruction algorithm is able to reconstruct the 3-D dose distribution in phantoms with a high accuracy. This result is obtained by combining portal dose images measured prior to treatment with an accurate dose calculation engine.