Modelling of the focal spot intensity distribution and the off-focal spot radiation in kilovoltage x-ray tubes for imaging.

X-ray tubes for medical applications typically generate x-rays by accelerating electrons, emitted from a cathode, with an interelectrode electric field, towards an anode target. X-rays are not emitted from one point, but from an irregularly shaped area on the anode, the focal spot. Focal spot intensity distributions and off-focal radiation negatively affect the imaging spatial resolution and broadens the beam penumbra. In this study, a Monte Carlo simulation model of an x-ray tube was developed to evaluate the spectral and spatial characteristics of off-focal radiation for multiple photon energies. Slit camera measurements were used to determine the horizontal and vertical intensity profiles of the small and the large focal spot of a diagnostic x-ray tube. First, electron beamlet weighting factors were obtained via an iterative optimization method to represent both focal spot sizes. These weighting factors were then used to extract off-focal spot radiation characteristics for the small and large focal spot sizes at 80, 100, and 120 kV. Finally, 120 kV simulations of a steel sphere (d = 4 mm) were performed to investigate image blurring with a point source, the small focal spot, and the large focal spot. The magnitude of off-focal radiation strongly depends on the anode size and the electric field coverage, and only minimally on the tube potential and the primary focal spot size. In conclusion, an x-ray tube Monte Carlo simulation model was developed to simulate focal spot intensity distributions and to evaluate off-focal radiation characteristics at several energies. This model can be further employed to investigate focal spot correction methods and to improve cone-beam CT image quality.

The effect of different image reconstruction techniques on pre-clinical quantitative imaging and dual-energy CT.

OBJECTIVE:: To analyse the effect of different image reconstruction techniques on image quality and dual energy CT (DECT) imaging metrics. METHODS:: A software platform for pre-clinical cone beam CT X-ray image reconstruction was built using the open-source reconstruction toolkit. Pre-processed projections were reconstructed with filtered back-projection and iterative algorithms, namely Feldkamp, Davis, and Kress (FDK), Iterative FDK, simultaneous algebraic reconstruction technique (SART), simultaneous iterative reconstruction technique and conjugate gradient. Imaging metrics were quantitatively assessed, using a quality assurance phantom, and DECT analysis was performed to determine the influence of each reconstruction technique on the relative electron density (ρe) and effective atomic number (Zeff) values. RESULTS:: Iterative reconstruction had favourable results for the DECT analysis: a significantly smaller spread for each material in the ρe-Zeff space and lower Zeff and ρe residuals (on average 24 and 25% lower, respectively). In terms of image quality assurance, the techniques FDK, Iterative FDK and SART provided acceptable results. The three reconstruction methods showed similar geometric accuracy, uniformity and CT number results. The technique SART had a contrast-to-noise ratio up to 76% higher for solid water and twice as high for Teflon, but resolution was up to 28% lower when compared to the other two techniques. CONCLUSIONS:: Advanced image reconstruction can be beneficial, but the benefit is small, and calculation times may be unacceptable with current technology. The use of targeted and downscaled reconstruction grids, larger, yet practicable, pixel sizes and GPU are recommended. ADVANCES IN KNOWLEDGE:: An iterative CBCT reconstruction platform was build using RTK.

Automatic multiatlas based organ at risk segmentation in mice.

OBJECTIVE:: During the treatment planning of a preclinical small animal irradiation, which has time limitations for reasons of animal wellbeing and workflow efficiency, the time consuming organ at risk (OAR) delineation is performed manually. This work aimed to develop, demonstrate, and quantitatively evaluate an automated contouring method for six OARs in a preclinical irritation treatment workflow. METHODS:: Microcone beam CT images of nine healthy mice were contoured with an in-house developed multiatlas-based image segmentation (MABIS) algorithm for six OARs: kidneys, eyes, heart, and brain. The automatic contouring was compared with the manual delineation using three quantitative metrics: the Dice Similarity Coefficient (DSC), 95th percentile Hausdorff Distance, and the centre of mass displacement. RESULTS:: A good agreement between manual and automatic contouring was found for OARs with sharp organ boundaries. For the brain and the heart, the median DSC was larger than 0.94, the median 95th Hausdorff Distance smaller than 0.44 mm, and the median centre of mass displacement smaller than 0.20 mm. Lower DSC values were obtained for the other OARs, but the median DSC was still larger than 0.74 for the left eye, 0.69 for the right eye, 0.89 for the left kidney and 0.80 for the right kidney. CONCLUSION:: The MABIS algorithm was able to delineate six OARs with a relatively high accuracy. Segmenting OARs with sharp organ boundaries performed better than low contrast OARs. ADVANCES IN KNOWLEDGE:: A MABIS algorithm is developed, evaluated, and demonstrated in a preclinical small animal irradiation research workflow.

A novel approach to EPID-based 3D volumetric dosimetry for IMRT and VMAT QA.

Intensity modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT) are relatively complex treatment delivery techniques and require quality assurance (QA) procedures. Pre-treatment dosimetric verification represents a fundamental QA procedure in daily clinical routine in radiation therapy. The purpose of this study is to develop an EPID-based approach to reconstruct a 3D dose distribution as imparted to a virtual cylindrical water phantom to be used for plan-specific pre-treatment dosimetric verification for IMRT and VMAT plans. For each depth, the planar 2D dose distributions acquired in air were back-projected and convolved by depth-specific scatter and attenuation kernels. The kernels were obtained by making use of scatter and attenuation models to iteratively estimate the parameters from a set of reference measurements. The derived parameters served as a look-up table for reconstruction of arbitrary measurements. The summation of the reconstructed 3D dose distributions resulted in the integrated 3D dose distribution of the treatment delivery. The accuracy of the proposed approach was validated in clinical IMRT and VMAT plans by means of gamma evaluation, comparing the reconstructed 3D dose distributions with Octavius measurement. The comparison was carried out using (3%, 3 mm) criteria scoring 99% and 96% passing rates for IMRT and VMAT, respectively. An accuracy comparable to the one of the commercial device for 3D volumetric dosimetry was demonstrated. In addition, five IMRT and five VMAT were validated against the 3D dose calculation performed by the TPS in a water phantom using the same passing rate criteria. The median passing rates within the ten treatment plans was 97.3%, whereas the lowest was 95%. Besides, the reconstructed 3D distribution is obtained without predictions relying on forward dose calculation and without external phantom or dosimetric devices. Thus, the approach provides a fully automated, fast and easy QA procedure for plan-specific pre-treatment dosimetric verification.

VOXSI: A voxelized single- and dual-energy CT scenario generator for quantitative imaging.

BACKGROUND AND PURPOSE: Dedicated CT simulation models have the potential to investigate several acquisition, reconstruction, or post-processing parameters without giving any radiation dose to patients. A software program was developed for the simulation and the analysis of single-energy and dual-energy CT images. Simulation and analysis functionalities of the software are described. MATERIALS AND METHODS: In the software, named VOXSI (VOXelized CT SImulator), the X-ray source, user specified simulation geometry, CT setup and the detector energy response can be varied. CT image reconstructions can be performed with an implementation of the ASTRA toolbox. In the DECT post processing toolkit, GUI tools are provided to calculate effective atomic number, relative electron density, pseudo-monoenergetic images, and material map images. Quantitative CT number validation, based on a RMI 467 tissue characterization phantom model, was performed between experimental and simulated CT scans at three different X-ray tube potentials (80, 120, and 140 kVp) with a third generation CT scanner. RESULTS: Overall, a good agreement was found for the mean CT numbers of the RMI 467 inserts. For all energies, the maximum difference in CT numbers between experimental and simulated data was below 17 HU for the soft tissues and below 48 HU for the osseous tissues. CONCLUSION: The software's simulation algorithm showed a good agreement between the CT measurements and CT simulations of the RMI 467 phantom at different energies. The capabilities of the software are demonstrated by an elaborated dual-energy CT research example.

A novel system for commissioning brachytherapy applicators: example of a ring applicator.

A novel system was developed to improve commissioning and quality assurance of brachytherapy applicators used in high dose rate (HDR). It employs an imaging panel to create reference images and to measure dwell times and dwell positions. As an example: two ring applicators of the same model were evaluated. An applicator was placed on the surface of an imaging panel and a HDR 192Ir source was positioned in an imaging channel above the panel to generate an image of the applicator, using the gamma photons of the brachytherapy source. The applicator projection image was overlaid with the images acquired by capturing the gamma photons emitted by the source dwelling inside the applicator. We verified 0.1, 0.2, 0.5 and 1.0 cm interdwell distances for different offsets, applicator inclinations and transfer tube curvatures. The data analysis was performed using in-house developed software capable of processing the data in real time, defining catheters and creating movies recording the irradiation procedure. One applicator showed up to 0.3 cm difference from the expected position for a specific dwell position. The problem appeared intermittently. The standard deviations of the remaining dwell positions (40 measurements) were less than 0.05 cm. The second ring applicator had a similar reproducibility with absolute coordinate differences from expected values ranging from -0.10 up to 0.18 cm. The curvature of the transfer tube can lead to differences larger than 0.1 cm whilst the inclination of the applicator showed a negligible effect. The proposed method allows the verification of all steps of the irradiation, providing accurate information about dwell positions and dwell times. It allows the verification of small interdwell positions (⩽0.1 cm) and reduces measurement time. In addition, no additional radiation source is necessary since the HDR 192Ir source is used to generate an image of the applicator.

Online pretreatment verification of high-dose rate brachytherapy using an imaging panel.

Brachytherapy is employed to treat a wide variety of cancers. However, an accurate treatment verification method is currently not available. This study describes a pre-treatment verification system that uses an imaging panel (IP) to verify important aspects of the treatment plan. A detailed modelling of the IP was only possible with an extensive calibration performed using a robotic arm. Irradiations were performed with a high dose rate (HDR) 192Ir source within a water phantom. An empirical fit was applied to measure the distance between the source and the detector so 3D Cartesian coordinates of the dwell positions can be obtained using a single panel. The IP acquires 7.14 fps to verify the dwell times, dwell positions and air kerma strength (Sk). A gynecological applicator was used to create a treatment plan that was registered with a CT image of the water phantom used during the experiments for verification purposes. Errors (shifts, exchanged connections and wrong dwell times) were simulated to verify the proposed verification system. Cartesian source positions (panel measurement plane) have a standard deviation of about 0.02 cm. The measured distance between the source and the panel (z-coordinate) have a standard deviation up to 0.16 cm and maximum absolute error of ≈0.6 cm if the signal is close to sensitive limit of the panel. The average response of the panel is very linear with Sk. Therefore, Sk measurements can be performed with relatively small errors. The measured dwell times show a maximum error of 0.2 s which is consistent with the acquisition rate of the panel. All simulated errors were clearly identified by the proposed system. The use of IPs is not common in brachytherapy, however, it provides considerable advantages. It was demonstrated that the IP can accurately measure Sk, dwell times and dwell positions.

Detection of anatomical changes in lung cancer patients with 2D time-integrated, 2D time-resolved and 3D time-integrated portal dosimetry: a simulation study.

The aim of this work is to assess the performance of 2D time-integrated (2D-TI), 2D time-resolved (2D-TR) and 3D time-integrated (3D-TI) portal dosimetry in detecting dose discrepancies between the planned and (simulated) delivered dose caused by simulated changes in the anatomy of lung cancer patients. For six lung cancer patients, tumor shift, tumor regression and pleural effusion are simulated by modifying their CT images. Based on the modified CT images, time-integrated (TI) and time-resolved (TR) portal dose images (PDIs) are simulated and 3D-TI doses are calculated. The modified and original PDIs and 3D doses are compared by a gamma analysis with various gamma criteria. Furthermore, the difference in the D 95% (ΔD 95%) of the GTV is calculated and used as a gold standard. The correlation between the gamma fail rate and the ΔD 95% is investigated, as well the sensitivity and specificity of all combinations of portal dosimetry method, gamma criteria and gamma fail rate threshold. On the individual patient level, there is a correlation between the gamma fail rate and the ΔD 95%, which cannot be found at the group level. The sensitivity and specificity analysis showed that there is not one combination of portal dosimetry method, gamma criteria and gamma fail rate threshold that can detect all simulated anatomical changes. This work shows that it will be more beneficial to relate portal dosimetry and DVH analysis on the patient level, rather than trying to quantify a relationship for a group of patients. With regards to optimizing sensitivity and specificity, different combinations of portal dosimetry method, gamma criteria and gamma fail rate should be used to optimally detect certain types of anatomical changes.

Investigating deformable image registration and scatter correction for CBCT-based dose calculation in adaptive IMPT.

PURPOSE: This work aims at investigating intensity corrected cone-beam x-ray computed tomography (CBCT) images for accurate dose calculation in adaptive intensity modulated proton therapy (IMPT) for prostate and head and neck (H&N) cancer. A deformable image registration (DIR)-based method and a scatter correction approach using the image data obtained from DIR as prior are characterized and compared on the basis of the same clinical patient cohort for the first time. METHODS: Planning CT (pCT) and daily CBCT data (reconstructed images and measured projections) of four H&N and four prostate cancer patients have been considered in this study. A previously validated Morphons algorithm was used for DIR of the planning CT to the current CBCT image, yielding a so-called virtual CT (vCT). For the first time, this approach was translated from H&N to prostate cancer cases in the scope of proton therapy. The warped pCT images were also used as prior for scatter correction of the CBCT projections for both tumor sites. Single field uniform dose and IMPT (only for H&N cases) treatment plans have been generated with a research version of a commercial planning system. Dose calculations on vCT and scatter corrected CBCT (CBCTcor) were compared by means of the proton range and a gamma-index analysis. For the H&N cases, an additional diagnostic replanning CT (rpCT) acquired within three days of the CBCT served as additional reference. For the prostate patients, a comprehensive contour comparison of CBCT and vCT, using a trained physician's delineation, was performed. RESULTS: A high agreement of vCT and CBCTcor was found in terms of the proton range and gamma-index analysis. For all patients and indications between 95% and 100% of the proton dose profiles in beam's eye view showed a range agreement of better than 3 mm. The pass rate in a (2%,2 mm) gamma-comparison was between 96% and 100%. For H&N patients, an equivalent agreement of vCT and CBCTcor to the reference rpCT was observed. However, for the prostate cases, an insufficient accuracy of the vCT contours retrieved from DIR was found, while the CBCTcor contours showed very high agreement to the contours delineated on the raw CBCT. CONCLUSIONS: For H&N patients, no considerable differences of vCT and CBCTcor were found. For prostate cases, despite the high dosimetric agreement, the DIR yields incorrect contours, probably due to the more pronounced anatomical changes in the abdomen and the reduced soft-tissue contrast in the CBCT. Using the vCT as prior, these inaccuracies can be overcome and images suitable for accurate delineation and dose calculation in CBCT-based adaptive IMPT can be retrieved from scatter correction of the CBCT projections.

Time-resolved versus time-integrated portal dosimetry: the role of an object's position with respect to the isocenter in volumetric modulated arc therapy.

The aim of this work is to compare time-resolved (TR) and time-integrated (TI) portal dosimetry, focussing on the role of an object's position with respect to the isocenter in volumetric modulated arc therapy (VMAT). Portal dose images (PDIs) are simulated and measured for different cases: a sphere (1), a bovine bone (2) and a patient geometry (3). For the simulated case (1) and the experimental case (2), several transformations are applied at different off-axis positions. In the patient case (3), three simple plans with different isocenters are created and pleural effusion is simulated in the patient. The PDIs before and after the sphere transformations, as well as the PDIs with and without simulated pleural effusion, are compared using a TI and TR gamma analysis. In addition, the performance of the TI and TR gamma analyses for the detection of real geometric changes in patients treated with clinical plans is investigated and a correlation analysis is performed between gamma fail rates and differences in dose volume histogram (DVH) metrics. The TI gamma analysis can show large differences in gamma fail rates for the same transformation at different off-axis positions (or for different plan isocenters). The TR gamma analysis, however, shows consistent gamma fail rates. For the detection of real geometric changes in patients treated with clinical plans, the TR gamma analysis has a higher sensitivity than the TI gamma analysis. However, the specificity for the TR gamma analysis is lower than for the TI gamma analysis. Both the TI and TR gamma fail rates show no correlation with changes in DVH metrics. This work shows that TR portal dosimetry is fundamentally superior to TI portal dosimetry, because it removes the strong dependence of the gamma fail rate on the off-axis position/plan isocenter. However, for 2D TR portal dosimetry, it is still difficult to interpret gamma fail rates in terms of changes in DVH metrics for patients treated with VMAT.

Dose rate mapping of VMAT treatments.

Human tissues exhibit a varying response to radiation dose depending on the dose rate and fractionation scheme used. Dose rate effects have been reported for different radiations, and tissue types. The literature indicates that there is not a significant difference in response for low-LET radiation when using dose rates between 1 Gy min(-1) and 12 Gy min(-1) but lower dose rates have an observable sparing effect on tissues and a differential effect between tissues. In intensity-modulated radiotherapy such as volumetric modulated arc therapy (VMAT) the dose can be delivered with a wide range of dose rates. In this work we developed a method based on time-resolved Monte Carlo simulations to quantify the dose rate frequency distribution for clinical VMAT treatments for three cancer sites, head and neck, lung, and pelvis within both planning target volumes (PTV) and normal tissues. The results show a wide range of dose rates are used to deliver dose in VMAT and up to 75% of the PTV can have its dose delivered with dose rates <1 Gy min(-1). Pelvic plans on average have a lower mean dose rate within the PTV than lung or head and neck plans but a comparable mean dose rate within the organs at risk. Two VMAT plans that fulfil the same dose objectives and constraints may be delivered with different dose rate distributions, particularly when comparing single arcs to multiple arc plans. It is concluded that for dynamic plans, the dose rate range used varies to a larger degree than previously assumed. The effect of the dose rate range in VMAT on clinical outcome is unknown.

Simulation of pseudo-CT images based on deformable image registration of ultrasound images: A proof of concept for transabdominal ultrasound imaging of the prostate during radiotherapy.

PURPOSE: Imaging of patient anatomy during treatment is a necessity for position verification and for adaptive radiotherapy based on daily dose recalculation. Ultrasound (US) image guided radiotherapy systems are currently available to collect US images at the simulation stage (USsim), coregistered with the simulation computed tomography (CT), and during all treatment fractions. The authors hypothesize that a deformation field derived from US-based deformable image registration can be used to create a daily pseudo-CT (CTps) image that is more representative of the patients' geometry during treatment than the CT acquired at simulation stage (CTsim). METHODS: The three prostate patients, considered to evaluate this hypothesis, had coregistered CT and US scans on various days. In particular, two patients had two US-CT datasets each and the third one had five US-CT datasets. Deformation fields were computed between pairs of US images of the same patient and then applied to the corresponding USsim scan to yield a new deformed CTps scan. The original treatment plans were used to recalculate dose distributions in the simulation, deformed and ground truth CT (CTgt) images to compare dice similarity coefficients, maximum absolute distance, and mean absolute distance on CT delineations and gamma index (γ) evaluations on both the Hounsfield units (HUs) and the dose. RESULTS: In the majority, deformation did improve the results for all three evaluation methods. The change in gamma failure for dose (γDose, 3%, 3 mm) ranged from an improvement of 11.2% in the prostate volume to a deterioration of 1.3% in the prostate and bladder. The change in gamma failure for the CT images (γCT, 50 HU, 3 mm) ranged from an improvement of 20.5% in the anus and rectum to a deterioration of 3.2% in the prostate. CONCLUSIONS: This new technique may generate CTps images that are more representative of the actual patient anatomy than the CTsim scan.

Is integrated transit planar portal dosimetry able to detect geometric changes in lung cancer patients treated with volumetric modulated arc therapy?

BACKGROUND: Geometric changes are frequent during the course of treatment of lung cancer patients. This may potentially result in deviations between the planned and actual delivered dose. Electronic portal imaging device (EPID)-based integrated transit planar portal dosimetry (ITPD) is a fast method for absolute in-treatment dose verification. The aim of this study was to investigate if ITPD could detect geometric changes in lung cancer patients. MATERIALS AND METHODS: A total of 460 patients treated with volumetric modulated arc therapy (VMAT) following daily cone beam computed tomography (CT)-based setup were visually inspected for geometrical changes on a daily basis. Forty-six patients were subject to changes and had a re-CT and an adaptive treatment plan. The reasons for adaptation were: change in atelectasis (n = 18), tumor regression (n = 9), change in pleural effusion (n = 8) or other causes (n = 11). The ITPDs were calculated on both the initial planning CT and the re-CT and compared with a global gamma (γ) evaluation (criteria: 3%\3mm). A treatment fraction failed when the percentage of pixels failing in the radiation fields exceeded 10%. Dose-volume histograms (DVHs) were compared between the initial plan versus the plan re-calculated on the re-CT. RESULTS: The ITPD threshold method detected 76% of the changes in atelectasis, while only 50% of the tumor regression cases and 42% of the pleural effusion cases were detected. Only 10% of the cases adapted for other reasons were detected with ITPD. The method has a 17% false-positive rate. No significant correlations were found between changes in DVH metrics and γ fail-rates. CONCLUSIONS: This study showed that most cases with geometric changes caused by atelectasis could be captured by ITPD, however for other causes ITPD is not sensitive enough to detect the clinically relevant changes and no predictive power of ITPD was found.

Weekly kilovoltage cone-beam computed tomography for detection of dose discrepancies during (chemo)radiotherapy for head and neck cancer.

BACKGROUND: Use of highly conformal radiotherapy in patients with head and neck carcinoma may lead to under-/overdosage of gross target volume (GTV) and organs at risk (OAR) due to changes in patients' anatomy. A method to achieve more effective radiation treatment combined with less toxicity is dose-guided radiotherapy (DGRT). The aim of this study was to evaluate discrepancies between planned and actually delivered radiation dose in head and neck patients and to identify predictive factors. METHODS: In this retrospective analysis, 20 patients with cT2-4 N0-3 M0 carcinoma originating from oropharynx, oral cavity, larynx and hypopharynx (Cohort 1), and seven patients with cT1-4 N0-3 M0 nasopharyngeal carcinoma (Cohort 2) treated with primary (chemo)radiotherapy and undergoing weekly kV-CBCT scans were included. Radiation dose was recalculated on 184 kV-CBCT images, which was quantified by D95% (GTV), Dmean (parotid and submandibular glands) and D2% (spinal cord). Predictive factors investigated for changes in these dose metrics were: gender, age, cT/N-stage, tumor grade, HPV-status, systemic therapy, body mass index at start of treatment, weight loss and volume change over the duration of the radiotherapy. RESULTS: There was no significant difference between the planned and delivered dose for GTV and OARs of Week 1 to subsequent weeks for Cohort 1. In Cohort 2, actually delivered Dmean to parotid glands was significant higher than planned dose (1.1 Gy, p = 0.002). No clinically relevant correlations between dose changes and predictive factors were found. CONCLUSION: Weekly dose calculations do not seem to improve dose delivery for patients with tumors of the oral cavity, oropharynx, larynx and hypopharynx. In patients with nasopharyngeal carcinoma, however, mid-treatment imaging may facilitate DGRT.

High dose rate and flattening filter free irradiation can be safely implemented in clinical practice.

PURPOSE: We hypothesize that flattening filter free (FFF) high dose rate irradiation will decrease cell survival in normal and cancer cells with more pronounced effects in DNA repair deficient cells. Additionally, we hypothesize that removal of the flattening filter will result in an enhanced relative biological effectiveness independent of the dose rate. MATERIALS AND METHODS: Clonogenic survival was assessed after exposure to dose rates of 4 or 24 Gy/min (FFF 10 megavolt [MV] photon beam) using a Varian TrueBeam accelerator. Additionally, cells were exposed to 4 Gy/min with or without flattening filter. Relative biological effectiveness estimations were performed comparing the different beam photon spectra. RESULTS: Cell survival in tumor and normal cell lines was not influenced by high dose rate irradiation. The intrinsic radiation sensitivity of DNA repair deficient cells was not affected by high dose rate compared to normal dose rate. Furthermore, the relative biological effectiveness was not significantly different from unity in any of the cell lines for both FFF and conventional flattened beam exposures. CONCLUSIONS: High dose rate irradiation did not affect long-term survival and DNA repair for cell lines of different tissues. This suggests that high dose rate does not influence treatment outcome or treatment toxicity and could be safely implemented in clinical routine.

What level of accuracy is achievable for preclinical dose painting studies on a clinical irradiation platform?

Advancements made over the past decades in both molecular imaging and radiotherapy planning and delivery have enabled studies that explore the efficacy of heterogeneous radiation treatment ("dose painting") of solid cancers based on biological information provided by different imaging modalities. In addition to clinical trials, preclinical studies may help contribute to identifying promising dose painting strategies. The goal of this current study was twofold: to develop a reproducible positioning and set-up verification protocol for a rat tumor model to be imaged and treated on a clinical platform, and to assess the dosimetric accuracy of dose planning and delivery for both uniform and positron emission tomography-computed tomography (PET-CT) based heterogeneous dose distributions. We employed a syngeneic rat rhabdomyosarcoma model, which was irradiated by volumetric modulated arc therapy (VMAT) with uniform or heterogeneous 6 MV photon dose distributions. Mean dose to the gross tumor volume (GTV) as a whole was kept at 12 Gy for all treatment arms. For the nonuniform plans, the dose was redistributed to treat the 30% of the GTV representing the biological target volume (BTV) with a dose 40% higher than the rest of the GTV (GTV - BTV) (~15 Gy was delivered to the BTV vs. ~10.7 Gy was delivered to the GTV - BTV). Cone beam computed tomography (CBCT) images acquired for each rat prior to irradiation were used to correctly reposition the tumor and calculate the delivered 3D dose. Film quality assurance was performed using a water-equivalent rat phantom. A comparison between CT or CBCT doses and film measurements resulted in passing rates >98% with a gamma criterion of 3%/2 mm using 2D dose images. Moreover, between the CT and CBCT calculated doses for both uniform and heterogeneous plans, we observed maximum differences of <2% for mean dose to the tumor and mean dose to the biological target volumes. In conclusion, we have developed a robust method for dose painting in a rat tumor model on a clinical platform, with a high accuracy achieved in the delivery of complex dose distributions. Our work demonstrates the technical feasibility of this approach and enables future investigations on the therapeutic effect of preclinical dose painting strategies using a state-of-the-art clinical platform.

HDR 192Ir source speed measurements using a high speed video camera.

PURPOSE: The dose delivered with a HDR (192)Ir afterloader can be separated into a dwell component, and a transit component resulting from the source movement. The transit component is directly dependent on the source speed profile and it is the goal of this study to measure accurate source speed profiles. METHODS: A high speed video camera was used to record the movement of a (192)Ir source (Nucletron, an Elekta company, Stockholm, Sweden) for interdwell distances of 0.25-5 cm with dwell times of 0.1, 1, and 2 s. Transit dose distributions were calculated using a Monte Carlo code simulating the source movement. RESULTS: The source stops at each dwell position oscillating around the desired position for a duration up to (0.026 ± 0.005) s. The source speed profile shows variations between 0 and 81 cm/s with average speed of ∼ 33 cm/s for most of the interdwell distances. The source stops for up to (0.005 ± 0.001) s at nonprogrammed positions in between two programmed dwell positions. The dwell time correction applied by the manufacturer compensates the transit dose between the dwell positions leading to a maximum overdose of 41 mGy for the considered cases and assuming an air-kerma strength of 48 000 U. The transit dose component is not uniformly distributed leading to over and underdoses, which is within 1.4% for commonly prescribed doses (3-10 Gy). CONCLUSIONS: The source maintains its speed even for the short interdwell distances. Dose variations due to the transit dose component are much lower than the prescribed treatment doses for brachytherapy, although transit dose component should be evaluated individually for clinical cases.

A novel time dependent gamma evaluation function for dynamic 2D and 3D dose distributions.

Modern external beam radiotherapy requires detailed verification and quality assurance so that confidence can be placed on both the delivery of a single treatment fraction and on the consistency of delivery throughout the treatment course. To verify dose distributions, a comparison between prediction and measurement must be made. Comparisons between two dose distributions are commonly performed using a Gamma evaluation which is a calculation of two quantities on a pixel by pixel basis; the dose difference, and the distance to agreement. By providing acceptance criteria (e.g. 3%, 3 mm), the function will find the most appropriate match within its two degrees of freedom. For complex dynamic treatments such as IMRT or VMAT it is important to verify the dose delivery in a time dependent manner and so a gamma evaluation that includes a degree of freedom in the time domain via a third parameter, time to agreement, is presented here. A C++ (mex) based gamma function was created that could be run on either CPU and GPU computing platforms that would allow a degree of freedom in the time domain. Simple test cases were created in both 2D and 3D comprising of simple geometrical shapes with well-defined boundaries varying over time. Changes of varying magnitude in either space or time were introduced and repeated gamma analyses were performed varying the criteria. A clinical VMAT case was also included, artificial air bubbles of varying size were introduced to a patient geometry, along with shifts of varying magnitude in treatment time. For all test cases where errors in distance, dose or time were introduced, the time dependent gamma evaluation could accurately highlight the errors.The time dependent gamma function presented here allows time to be included as a degree of freedom in gamma evaluations. The function allows for 2D and 3D data sets which are varying over time to be compared using appropriate criteria without penalising minor offsets of subsequent radiation fields in time.

Time dependent pre-treatment EPID dosimetry for standard and FFF VMAT.

Methods to calibrate Megavoltage electronic portal imaging devices (EPIDs) for dosimetry have been previously documented for dynamic treatments such as intensity modulated radiotherapy (IMRT) using flattened beams and typically using integrated fields. While these methods verify the accumulated field shape and dose, the dose rate and differential fields remain unverified. The aim of this work is to provide an accurate calibration model for time dependent pre-treatment dose verification using amorphous silicon (a-Si) EPIDs in volumetric modulated arc therapy (VMAT) for both flattened and flattening filter free (FFF) beams. A general calibration model was created using a Varian TrueBeam accelerator, equipped with an aS1000 EPID, for each photon spectrum 6 MV, 10 MV, 6 MV-FFF, 10 MV-FFF. As planned VMAT treatments use control points (CPs) for optimization, measured images are separated into corresponding time intervals for direct comparison with predictions. The accuracy of the calibration model was determined for a range of treatment conditions. Measured and predicted CP dose images were compared using a time dependent gamma evaluation using criteria (3%, 3 mm, 0.5 sec). Time dependent pre-treatment dose verification is possible without an additional measurement device or phantom, using the on-board EPID. Sufficient data is present in trajectory log files and EPID frame headers to reliably synchronize and resample portal images. For the VMAT plans tested, significantly more deviation is observed when analysed in a time dependent manner for FFF and non-FFF plans than when analysed using only the integrated field. We show EPID-based pre-treatment dose verification can be performed on a CP basis for VMAT plans. This model can measure pre-treatment doses for both flattened and unflattened beams in a time dependent manner which highlights deviations that are missed in integrated field verifications.

Measured vs simulated portal images for low MU fields on three accelerator types: possible consequences for 2D portal dosimetry.

PURPOSE: As external beam treatment plans become more dynamic and the dose to normal tissue is further constrained, treatments may consist of a larger number of beams, each delivering smaller doses (or monitor units, MU), in, e.g., volumetric modulated arc therapy (VMAT). Electronic portal imaging devices (EPID) may be used to verify external beam treatments on integrated fractions as well as in a more time dependent manner such as field by field. For treatment verification performed during a fraction (e.g., individual fields or VMAT control points), the lower limit of EPID measurement capability becomes important. The authors quantified the signal and timing accuracy of EPID images for low MU intensity modulated radiotherapy (IMRT) and conformal fields. METHODS: EPID images were collected from three different vendor's accelerators for low MU fields and compared to expected images. Simulations were performed to replicate the EPID acquisition pattern and to enhance the understanding of EPID readout schemes. RESULTS: Large discrepancies between observed and predicted images were noted due to an under-response to single low MU fields. It is shown that a variability of up to 37% can be observed for low MU fields in clinically used EPID acquisition modes and that the majority of this variability can be accounted for by the readout scheme, integration, and timing of EPID acquisitions. Simulations have confirmed the causes of the discrepancies. The occurrence and extent of the variation has been estimated for clinical settings. CONCLUSIONS: Incorrect absolute EPID signals collected for low MU fields in external beam treatments will negatively affect quantitative applications such as individual field based EPID dosimetry, typically appearing as an underdose, unless corrections to currently employed EPID readout schemes are made.