COMP Report: An updated algorithm to estimate medical physics staffing levels for radiation oncology.

PURPOSE: Medical physics staffing models require periodic review due to the rapid evolution of technology and clinical techniques in radiation oncology. We present an update to a grid-based physics staffing algorithm for radiation oncology (originally published in 2012) that has been widely used in Canada over the last decade. MATERIALS AND METHODS: The physics staffing algorithm structure was modified to improve the clarity and consistency of input data. We collected information on clinical procedures, equipment inventory, and teaching activities from 15 radiation treatment centers in the province of Ontario from April 1, 2018, to March 31, 2019. Using these data sets, the algorithm's weighting parameters were adjusted to align the prediction of full-time equivalent (FTE) personnel with actual staffing levels in Ontario. The algorithm computes FTE estimates for medical physicists, physics assistants, engineering (electrical and mechanical), and information technology (IT) support. The performance of the algorithm was also tested in eight Canadian cancer centers outside of Ontario. RESULTS: The mean difference between the algorithm and actual staffing for the 23 Canadian cancer centers did not exceed 0.5 FTE for any staffing group. The results were slightly better in Ontario than in other provinces, as expected since the algorithm was optimized using Ontario data. There was a linear correlation between the algorithm predictions and the number of annual-treated cases for physicists, and physicists plus physics assistants. For other staff categories, the algorithm weighting parameters were not significantly altered, except for a reduction in mechanical engineering staff. Comparison with other published models suggests that the updated algorithm should be considered as a minimum recommended staffing level for the clinical support of radiation oncology programs. CONCLUSIONS: We support the use of grid-based physics staffing algorithms that account for clinical workload with flexibility to adapt to local conditions with variable academic and research demands.

Technical Note: A fast inverse direct aperture optimization algorithm for volumetric-modulated arc therapy.

PURPOSE: In a recent article, our group proposed a fast direct aperture optimization (DAO) algorithm for fixed-gantry intensity-modulated radiation therapy (IMRT) called fast inverse direct aperture optimization (FIDAO). When tested on fixed-gantry IMRT plans, we observed up to a 200-fold increase in the optimization speed. Compared to IMRT, rotational volumetric-modulated arc therapy (VMAT) is a much larger optimization problem and has many more delivery constraints. The purpose of this work is to extend and evaluate FIDAO for inverse planning of VMAT plans. METHODS: A prototype FIDAO algorithm for VMAT treatment planning was developed in MATLAB using the open-source treatment planning toolkit matRad (v2.2 dev_VMAT build). VMAT treatment plans using one 3600 arc were generated on the AAPM TG-119 phantom, as well as sample clinical liver and prostate cases. The plans were created by first performing fluence map optimization on 28° equispaced beams, followed by aperture sequencing and arc sequencing with a gantry angular sampling rate of 4°. After arc sequencing, a copy of the plan underwent DAO using the prototype FIDAO algorithm, while another copy of the plan underwent DAO using matRad's DAO method, which served as the conventional algorithm. RESULTS: Both algorithms achieved similar plan quality, although the FIDAO plans had considerably fewer hot spots in the unspecified normal tissue. The optimization time (number of iterations) for FIDAO and the conventional DAO algorithm, respectively, were: 65 s (245) vs 602 s (275) in the TG-119 phantom case; 25 s (85) vs 803 s (159) in the liver case; and 99 s (174) vs 754 s (149) in the prostate case. CONCLUSIONS: This study demonstrated promising speed enhancements in using FIDAO for the direct aperture optimization of VMAT plans.

[18F]FDG cardiac PET imaging in a canine model of radiation-induced cardiovascular disease associated with breast cancer radiotherapy.

Radiotherapy for the treatment of left-sided breast cancer increases the long-term risk of cardiovascular disease. The purpose of the present study was to noninvasively image the progression of radiation-induced cardiac inflammation in a large animal model using a hybrid PET and MRI system. Five canines were imaged using [18F]fluorodeoxyglucose PET to assess changes in myocardial inflammation. All animals were imaged at baseline, 1 wk, and 1, 3, 6, and 12 mo after focused cardiac external beam irradiation with image guidance. Radiation was delivered in a single fraction. The linear quadratic model was used to convert a typical multifractionated heart dose to a corrected single-fraction biologically equivalent dose. Immunohistochemistry was performed on excised left ventricular tissue samples from all five irradiated canines and one nonirradiated control canine to confirm the presence of inflammation. The mean doses delivered to the entire heart, left ventricle, left anterior descending artery, and left circumflex artery were 1.7 ± 0.2, 2.7 ± 0.2, 5.5 ± 0.9, and 1.1 ± 0.4 Gy, respectively. FDG standard uptake values remained persistently elevated compared with baseline (1.1 ± 0.03 vs. 2.6 ± 0.19, P < 0.05). The presence of myocardial inflammation was confirmed histologically and correlated with myocardial dose. This study suggests a global inflammatory response that is persistent up to 12 mo postirradiation. Inflammation PET imaging should be considered in future clinical studies to monitor the early changes in cardiac function that may play a role in the ultimate development of radiation-induced cardiac toxicity. NEW & NOTEWORTHY Using advanced cardiac PET imaging, we have shown the spatial and quantitative relationship between radiation dose deposition and temporal changes in inflammation. We have shown that the progression of radiation-induced cardiac inflammation is immediate and does not subside for up to 1 yr after radiation. Results are presented in a large animal model that closely resembles the size and vessel architecture of humans. The proposed imaging protocol can be easily replicated for clinical use.

A fast inverse direct aperture optimization algorithm for intensity-modulated radiation therapy.

PURPOSE: The goal of this work was to develop and evaluate a fast inverse direct aperture optimization (FIDAO) algorithm for IMRT treatment planning and plan adaptation. METHODS: A previously proposed fluence map optimization algorithm called fast inverse dose optimization (FIDO) was extended to optimize the aperture shapes and weights of IMRT beams. FIDO is a very fast fluence map optimization algorithm for IMRT that finds the global minimum using direct matrix inversion without unphysical negative beam weights. In this study, an equivalent second-order Taylor series expansion of the FIDO objective function was used, which allowed for the objective function value and gradient vector to be computed very efficiently during direct aperture optimization, resulting in faster optimization. To evaluate the speed gained with FIDAO, a proof-of-concept algorithm was developed in MATLAB using an interior-point optimization method to solve the reformulated aperture-based FIDO problem. The FIDAO algorithm was used to optimize four step-and-shoot IMRT cases: on the AAPM TG-119 phantom as well as a liver, prostate, and head-and-neck clinical cases. Results were compared with a conventional DAO algorithm that uses the same interior-point method but using the standard formulation of the objective function and its gradient vector. RESULTS: A substantial gain in optimization speed was obtained with the prototype FIDAO algorithm compared to the conventional DAO algorithm while producing plans of similar quality. The optimization time (number of iterations) for the prototype FIDAO algorithm vs the conventional DAO algorithm was 0.3 s (17) vs 56.7 s (50); 2.0 s (28) vs 134.1 s (57); 2.5 s (26) vs 180.6 s (107); and 6.7 s (20) vs 469.4 s (482) in the TG-119 phantom, liver, prostate, and head-and-neck examples, respectively. CONCLUSIONS: A new direct aperture optimization algorithm based on FIDO was developed. For the four IMRT test cases examined, this algorithm executed approximately 70-200 times faster without compromising the IMRT plan quality.

Online daily assessment of dose change in head and neck radiotherapy without dose-recalculation.

BACKGROUND: Head and neck cancers are commonly treated with radiation therapy, but due to possible volume changes, plan adaptation may be required during the course of treatment. Currently, plan adaptations consume significant clinical resources. Existing methods to evaluate the need for plan adaptation requires deformable image registration (DIR) to a new CT simulation or daily cone beam CT (CBCT) images and the recalculation of the dose distribution. In this study, we explore a tool to assist the decision for plan adaptation using a CBCT without re-computation of dose, allowing for rapid online assessment. METHODS: This study involved 18 head and neck cancer patients treated with CBCT image guidance who had their treatment plan modified based on a new CT simulation (ReCT). Dose changes were estimated using different methods and compared to the current gold standard of using DIR between the planning CT scan (PCT) and ReCT with recomputed dose. The first and second methods used DIR between the PCT and daily CBCT with the planned dose or recalculated dose from the ReCT respectively, with the dose transferred to the CBCT using rigid registration. The necessity of plan adaptation was assessed by the change in dose to 95% of the planning target volume (D95) and mean dose to the parotids. RESULTS: The treatment plans were adapted clinically for all 18 patients but only 7 actually needed an adaptation yielding 11 unnecessary adaptations. Applying a method using the daily CBCT with the planned dose distribution would have yielded only four unnecessary adaptations and no missed adaptations: a significant improvement from that done clinically. CONCLUSION: Using the DIR between the planning CT and daily CBCT can flag cases for plan adaptation before every fraction while not requiring a new re-planning CT scan and dose recalculation.

Patient-specific calibration of cone-beam computed tomography data sets for radiotherapy dose calculations and treatment plan assessment.

PURPOSE: In this work, we propose a new method of calibrating cone beam computed tomography (CBCT) data sets for radiotherapy dose calculation and plan assessment. The motivation for this patient-specific calibration (PSC) method is to develop an efficient, robust, and accurate CBCT calibration process that is less susceptible to deformable image registration (DIR) errors. METHODS: Instead of mapping the CT numbers voxel-by-voxel with traditional DIR calibration methods, the PSC methods generates correlation plots between deformably registered planning CT and CBCT voxel values, for each image slice. A linear calibration curve specific to each slice is then obtained by least-squares fitting, and applied to the CBCT slice's voxel values. This allows each CBCT slice to be corrected using DIR without altering the patient geometry through regional DIR errors. A retrospective study was performed on 15 head-and-neck cancer patients, each having routine CBCTs and a middle-of-treatment re-planning CT (reCT). The original treatment plan was re-calculated on the patient's reCT image set (serving as the gold standard) as well as the image sets produced by voxel-to-voxel DIR, density-overriding, and the new PSC calibration methods. Dose accuracy of each calibration method was compared to the reference reCT data set using common dose-volume metrics and 3D gamma analysis. A phantom study was also performed to assess the accuracy of the DIR and PSC CBCT calibration methods compared with planning CT. RESULTS: Compared with the gold standard using reCT, the average dose metric differences were ≤ 1.1% for all three methods (PSC: -0.3%; DIR: -0.7%; density-override: -1.1%). The average gamma pass rates with thresholds 3%, 3 mm were also similar among the three techniques (PSC: 95.0%; DIR: 96.1%; density-override: 94.4%). CONCLUSIONS: An automated patient-specific calibration method was developed which yielded strong dosimetric agreement with the results obtained using a re-planning CT for head-and-neck patients.

Technical Note: Evaluation of an iterative reconstruction algorithm for optical CT radiation dosimetry.

PURPOSE: Iterative CT reconstruction algorithms are gaining popularity as GPU-based computation becomes more accessible. These algorithms are desirable in x-ray CT for their ability to achieve similar image quality at a fraction of the dose required for standard filtered backprojection reconstructions. In optical CT dosimetry, the noise reduction capability of such algorithms is similarly desirable because noise has a detrimental effect on the precision of dosimetric analysis, and can create misleading test results. In this study, we evaluate an iterative CT reconstruction algorithm for gel dosimetry, with special attention to the challenging dosimetry of small fields. METHODS: An existing ordered subsets convex algorithm using total variation minimization regularization (OSC-TV) was implemented. Three datasets, which represent the extreme cases of gel dosimetry, were examined: a large, 15 cm diameter uniform phantom, a 1.35 cm diameter finger phantom, and a 15 cm gel dosimeter irradiated with 3 × 3, 2 × 2, 1 × 1, and 0.6 × 0.6 cm fields. These were scanned on an in-house scanning laser system, and reconstructed with both filtered backprojection and OSC-TV with a range of regularization constants. The contrast to artifact + noise ratio (CANR) and penumbra width measurements (80% to 20% and 95% to 5% distances) were used to compare reconstructions. RESULTS: Our results showed that OSC-TV can achieve 3-5× improvement in contrast to artifact + noise ratio compared to filtered backprojection, while preserving the shape of steep dose gradients. For very small objects (≤ 0.6 × 0.6 cm fields in a 16 × 16 cm field of view), the mean value in the center of the object can be suppressed if the regularization constant is improperly set, which must be avoided. CONCLUSIONS: Overall, the results indicate that OSC-TV is a suitable reconstruction algorithm for gel dosimetry, provided care is taken in setting the regularization parameter when reconstructing objects that are small compared to the scanner field of view.

Using gamma index to flag changes in anatomy during image-guided radiation therapy of head and neck cancer.

During radiation therapy of head and neck cancer, the decision to consider replanning a treatment because of anatomical changes has significant resource implications. We developed an algorithm that compares cone-beam computed tomography (CBCT) image pairs and provides an automatic alert as to when remedial action may be required. Retrospective CBCT data from ten head and neck cancer patients that were replanned during their treatment was used to train the algorithm on when to recommend a repeat CT simulation (re-CT). An additional 20 patients (replanned and not replanned) were used to validate the predictive power of the algorithm. CBCT images were compared in 3D using the gamma index, combining Hounsfield Unit (HU) difference with distance-to-agreement (DTA), where the CBCT study acquired on the first fraction is used as the reference. We defined the match quality parameter (MQPx ) as a difference between the xth percentiles of the failed-pixel histograms calculated from the reference gamma comparison and subsequent comparisons, where the reference gamma comparison is taken from the first two CBCT images acquired during treatment. The decision to consider re-CT was based on three consecutive MQP values being less than or equal to a threshold value, such that re-CT recommendations were within ±3 fractions of the actual re-CT order date for the training cases. Receiver-operator characteristic analysis showed that the best trade-off in sensitivity and specificity was achieved using gamma criteria of 3 mm DTA and 30 HU difference, and the 80th percentile of the failed-pixel histogram. A sensitivity of 82% and 100% was achieved in the training and validation cases, respectively, with a false positive rate of ~30%. We have demonstrated that gamma analysis of CBCT-acquired anatomy can be used to flag patients for possible replanning in a manner consistent with local clinical practice guidelines.

Representing the dosimetric impact of deformable image registration errors.

Deformable image registration (DIR) is emerging as a tool in radiation therapy for calculating the cumulative dose distribution across multiple fractions of treatment. Unfortunately, due to the variable nature of DIR algorithms and dependence of performance on image quality, registration errors can result in dose accumulation errors. In this study, landmarked images were used to characterize the DIR error throughout an image space and determine its impact on dosimetric analysis. Ten thoracic 4DCT images with 300 landmarks per image study matching the end-inspiration and end-expiration phases were obtained from 'dir-labs'. DIR was performed using commercial software MIM Maestro. The range of dose uncertainty (RDU) was calculated at each landmark pair as the maximum and minimum of the doses within a sphere around the landmark in the end-expiration phase. The radius of the sphere was defined by a measure of DIR error which included either the actual DIR error, mean DIR error per study, constant errors of 2 or 5 mm, inverse consistency error, transitivity error or the distance discordance metric (DDM). The RDUs were evaluated using the magnitude of dose uncertainty (MDU) and inclusion rate (IR) of actual error lying within the predicted RDU. The RDU was calculated for 300 landmark pairs on each 4DCT study for all measures of DIR error. The most representative RDU was determined using the actual DIR error with a MDU of 2.5 Gy and IR of 97%. Across all other measures of DIR error, the DDM was most predictive with a MDU of 2.5 Gy and IR of 86%, closest to the actual DIR error. The proposed method represents the range of dosimetric uncertainty of DIR error using either landmarks at specific voxels or measures of registration accuracy throughout the volume.

Scanning laser optical computed tomography system for large volume 3D dosimetry.

Stray light causes artifacts in optical computed tomography (CT) that negatively affect the accuracy of radiation dosimetry in gels or solids. Scatter effects are exacerbated by a large dosimeter volume, which is desirable for direct verification of modern radiotherapy treatment plans such as multiple-isocenter radiosurgery. The goal in this study was to design and characterize an optical CT system that achieves high accuracy primary transmission measurements through effective stray light rejection, while maintaining sufficient scan speed for practical application. We present an optical imaging platform that uses a galvanometer mirror for horizontal scanning, and a translation stage for vertical movement of a laser beam and small area detector for minimal stray light production and acceptance. This is coupled with a custom lens-shaped optical CT aquarium for parallel ray sampling of projections. The scanner images 15 cm diameter, 12 cm height cylindrical volumes at 0.33 mm resolution in approximately 30 min. Attenuation coefficients reconstructed from CT scans agreed with independent cuvette measurements within 2% for both absorbing and scattering solutions as well as small 1.25 cm diameter absorbing phantoms placed within a large, scattering medium that mimics gel. Excellent linearity between the optical CT scanner and the independent measurement was observed for solutions with between 90% and 2% transmission. These results indicate that the scanner should achieve highly accurate dosimetry of large volume dosimeters in a reasonable timeframe for clinical application to radiotherapy dose verification procedures.

Optical CT imaging of solid radiochromic dosimeters in mismatched refractive index solutions using a scanning laser and large area detector.

PURPOSE: The practical use of the PRESAGE® solid plastic dosimeter is limited by the inconvenience of immersing it in high-viscosity oils to achieve refractive index matching for optical computed tomography (CT) scanning. The oils are slow to mix and difficult to clean from surfaces, and the dosimeter rotation can generate dynamic Schlieren inhomogeneity patterns in the reference liquid, limiting the rotational and overall scan speed. Therefore, it would be beneficial if lower-viscosity, water-based solutions with slightly unmatched refractive index could be used instead. The purpose of this work is to demonstrate the feasibility of allowing mismatched conditions when using a scanning laser system with a large acceptance angle detector. A fiducial-based ray path measurement technique is combined with an iterative CT reconstruction algorithm to reconstruct images. METHODS: A water based surrounding liquid with a low viscosity was selected for imaging PRESAGE® solid dosimeters. Liquid selection was optimized to achieve as high a refractive index as possible while avoiding rotation-induced Schlieren effects. This led to a refractive index mismatch of 6% between liquid and dosimeters. Optical CT scans were performed with a fan-beam scanning-laser optical CT system with a large area detector to capture most of the refracted rays. A fiducial marker placed on the wall of a cylindrical sample occludes a given light ray twice. With knowledge of the rotation angle and the radius of the cylindrical object, the actual internal path of each ray through the dosimeter can be calculated. Scans were performed with 1024 projections of 512 data samples each, and rays were rebinned to form 512 parallel-beam projections. Reconstructions were performed on a 512 × 512 grid using 100 iterations of the SIRT iterative CT algorithm. Proof of concept was demonstrated with a uniformly attenuating solution phantom. PRESAGE® dosimeters (11 cm diameter) were irradiated with Cobalt-60 irradiator to achieve either a uniform dose or a 2-level "step-dose" pattern. RESULTS: With 6% refractive index mismatching, a circular field of view of 85% of the diameter of a cylindrical sample can be reconstructed accurately. Reconstructed images of the test solution phantom were uniform (within 3%) inside this radius. However, the dose responses of the PRESAGE® samples were not spatially uniform, with variations of at least 5% in sensitivity. The variation appears as a "cupping" artifact with less sensitivity in the middle than at the periphery of the PRESAGE® cylinder. Polarization effects were also detected for these samples. CONCLUSIONS: The fiducial-based ray path measurement scheme, coupled with an iterative reconstruction algorithm, enabled optical CT scanning of PRESAGE® dosimeters immersed in mismatched refractive index solutions. However, improvements to PRESAGE® dose response uniformity are required.

Evaluation of unified intensity-modulated arc therapy for the radiotherapy of head-and-neck cancer.

PURPOSE: Recently our group developed a unified intensity-modulated arc therapy (UIMAT) technique which allows for the simultaneous inverse-optimization and the combined delivery of volume-modulated arc therapy (VMAT) and intensity-modulated radiation therapy (IMRT). The aim of this study was to evaluate the dosimetric benefits of UIMAT plans for radiation treatment of complex head-and-neck cancer cases. METHODS AND MATERIALS: A retrospective treatment planning study was performed on 30 head-and-neck cases, 15 of which were treated clinically with VMAT while the other 15 were treated with step-and-shoot IMRT. These cases were re-planned using our UIMAT technique and the results were compared with the clinically delivered plans. Plans were assessed in terms of clinically relevant metrics describing target volume coverage, dose conformity, and the sparing of organs at risk. RESULTS: When compared to stand-alone VMAT or IMRT, UIMAT plans offered slightly better tumor volume coverage (Median D95: 98.1% vs. 97.5%, p=0.01) and similar dose conformity (Median CI: 0.69 vs. 0.69, p=0.09). More significantly, UIMAT plans had substantially lower doses to all organs at risk, including the spinal cord (Median D2%: 29.9Gy vs. 35.6Gy, p<0.01), brainstem (Median D2%: 21.2Gy vs. 25.6Gy, p<0.01), left parotid (Median DMean: 26.1Gy vs. 28.0Gy, p<0.01), and right parotid (Median DMean: 23.6Gy vs. 27.2Gy, p<0.01). The reduction in OAR doses did not result from the redistribution of dose to unspecified tissue. Furthermore, UIMAT plans can be delivered with comparable delivery times to VMAT (Median time: 135s vs. 168s, p=0.394) but with fewer monitor units (Median MU: 486 vs. 635, p<0.01). CONCLUSIONS: Compared to stand-alone IMRT or VMAT, UIMAT was demonstrated to have a dosimetric advantage for the radiation treatment of head-and-neck cancer.

Stray light in cone beam optical computed tomography: I. Measurement and reduction strategies with planar diffuse source.

Optical cone-beam computed tomographic (CBCT) scanning of 3D radiochromic dosimeters may provide a practical method for 3D dose verification in radiation therapy. However, in cone-beam geometry stray light contaminates the projection images, degrading the accuracy of reconstructed linear attenuation coefficients. Stray light was measured using a beam pass aperture array (BPA) and structured illumination methods. The stray-to-primary ray ratio (SPR) along the central axis was found to be 0.24 for a 5% gelatin hydrogel, representative of radiochromic hydrogels. The scanner was modified by moving the spectral filter from the detector to the source, changing the light's spatial fluence pattern and lowering the acceptance angle by extending distance between the source and object. These modifications reduced the SPR significantly from 0.24 to 0.06. The accuracy of the reconstructed linear attenuation coefficients for uniform carbon black liquids was compared to independent spectrometer measurements. Reducing the stray light increased the range of accurate transmission readings. In order to evaluate scanner performance for the more challenging application to small field dosimetry, a carbon black finger gel phantom was prepared. Reconstructions of the phantom from CBCT and fan-beam CT scans were compared. The modified source resulted in improved agreement. Subtraction of residual stray light, measured with BPA or structured illumination from each projection further improved agreement. Structured illumination was superior to BPA for measuring stray light for the smaller 1.2 and 0.5 cm diameter phantom fingers. At the costs of doubling the scanner size and tripling the number of scans, CBCT reconstructions of low-scattering hydrogel dosimeters agreed with those of fan-beam CT scans.

Stray light in cone beam optical computed tomography: II. Reduction using a convergent light source.

Optical cone beam computed tomography (CBCT) using a broad beam and CCD camera is a fast method for densitometry of 3D optical gel dosimeters. However, diffuse light sources introduce considerable stray light into the imaging system, leading to underestimation of attenuation coefficients and non-uniformities in CT images unless corrections are applied to each projection image. In this study, the light source of a commercial optical CT scanner is replaced with a convergent cone beam source consisting of almost exclusively image forming primary rays. The convergent source is achieved using a small isotropic source and a Fresnel lens. To characterize stray light effects, full-field cone beam CT imaging is compared to fan beam CT (FBCT) using a 1 cm high fan beam aperture centered on the optic axis of the system. Attenuating liquids are scanned within a large 96 mm diameter uniform phantom and in a small 13.5 mm diameter finger phantom. For the uniform phantom, cone and fan beam CT attenuation coefficients agree within a maximum deviation of (1 ± 2)% between mean values over a wide range from 0.036 to 0.43 cm(-1). For the finger phantom, agreement is found with a maximum deviation of (4 ± 2)% between mean values over a range of 0.1-0.47 cm(-1). With the convergent source, artifacts associated with refractive index mismatch and vessel optical features are more pronounced. Further optimization of the source size to achieve a balance between quantitative accuracy and artifact reduction should enable practical, accurate 3D dosimetry, avoiding time consuming 3D scatter measurements.

Feasibility of a unified approach to intensity-modulated radiation therapy and volume-modulated arc therapy optimization and delivery.

PURPOSE: To study the feasibility of unified intensity-modulated arc therapy (UIMAT) which combines intensity-modulated radiotherapy (IMRT) and volumetric-modulated arc therapy (VMAT) optimization and delivery to produce superior radiation treatment plans, both in terms of dose distribution and efficiency of beam delivery when compared with either VMAT or IMRT alone. METHODS: An inverse planning algorithm for UIMAT was prototyped within the pinnacle treatment planning system (Philips Healthcare). The IMRT and VMAT deliveries are unified within the same arc, with IMRT being delivered at specific gantry angles within the arc. Optimized gantry angles for the IMRT and VMAT phases are assigned automatically by the inverse optimization algorithm. Optimization of the IMRT and VMAT phases is done simultaneously using a direct aperture optimization algorithm. Five treatment plans each for prostate, head and neck, and lung were generated using a unified optimization technique and compared with clinical IMRT or VMAT plans. Delivery verification was performed with an ArcCheck phantom (Sun Nuclear) on a Varian TrueBeam linear accelerator (Varian Medical Systems). RESULTS: In this prototype implementation, the UIMAT plans offered the same target dose coverage while reducing mean doses to organs at risk by 8.4% for head-and-neck cases, 5.7% for lung cases, and 3.5% for prostate cases, compared with the VMAT or IMRT plans. In addition, UIMAT can be delivered with similar efficiency as VMAT. CONCLUSIONS: In this proof-of-concept work, a novel radiation therapy optimization and delivery technique that interlaces VMAT or IMRT delivery within the same arc has been demonstrated. Initial results show that unified VMAT/IMRT has the potential to be superior to either standard IMRT or VMAT.

Forcing lateral electron disequilibrium to spare lung tissue: a novel technique for stereotactic body radiation therapy of lung cancer.

Stereotactic body radiation therapy (SBRT) has quickly become a preferred treatment option for early-stage lung cancer patients who are ineligible for surgery. This technique uses tightly conformed megavoltage (MV) x-ray beams to irradiate a tumour with ablative doses in only a few treatment fractions. Small high energy x-ray fields can cause lateral electron disequilibrium (LED) to occur within low density media, which can reduce tumour dose. These dose effects may be challenging to predict using analytic dose calculation algorithms, especially at higher beam energies. As a result, previous authors have suggested using low energy photons (<10 MV) and larger fields (>5 × 5 cm(2)) for lung cancer patients to avoid the negative dosimetric effects of LED. In this work, we propose a new form of SBRT, described as LED-optimized SBRT (LED-SBRT), which utilizes radiotherapy (RT) parameters designed to cause LED to advantage. It will be shown that LED-SBRT creates enhanced dose gradients at the tumour/lung interface, which can be used to manipulate tumour dose, and/or normal lung dose. To demonstrate the potential benefits of LED-SBRT, the DOSXYZnrc (National Research Council of Canada, Ottawa, ON) Monte Carlo (MC) software was used to calculate dose within a cylindrical phantom and a typical lung patient. 6 MV or 18 MV x-ray fields were focused onto a small tumour volume (diameter ∼1 cm). For the phantom, square fields of 1 × 1 cm(2), 3 × 3 cm(2), or 5 × 5 cm(2) were applied. However, in the patient, 3 × 1 cm(2), 3 × 2 cm(2), 3 × 2.5 cm(2), or 3 × 3 cm(2) field sizes were used in simulations to assure target coverage in the superior-inferior direction. To mimic a 180° SBRT arc in the (symmetric) phantom, a single beam profile was calculated, rotated, and beams were summed at 1° segments to accumulate an arc dose distribution. For the patient, a 360° arc was modelled with 36 equally weighted (and spaced) fields focused on the tumour centre. A planning target volume (PTV) was generated by considering the extent of tumour motion over the patient's breathing cycle and set-up uncertainties. All patient dose results were normalized such that at least 95% of the PTV received at least 54 Gy (i.e. D95 = 54 Gy). Further, we introduce 'LED maps' as a novel clinical tool to compare the magnitude of LED resulting from the various SBRT arc plans. Results from the phantom simulation suggest that the best lung sparing occurred for RT parameters that cause severe LED. For equal tumour dose coverage, normal lung dose (2 cm outside the target region) was reduced from 92% to 23%, comparing results between the 18 MV (5 × 5 cm(2)) and 18 MV (1 × 1 cm(2)) arc simulations. In addition to reduced lung dose for the 18 MV (1 × 1 cm(2)) arc, maximal tumour dose increased beyond 125%. Thus, LED can create steep dose gradients to spare normal lung, while increasing tumour dose levels (if desired). In the patient simulation, a LED-optimized arc plan was designed using either 18 MV (3 × 1 cm(2)) or 6 MV (3 × 3cm(2)) beams. Both plans met the D95 dose coverage requirement for the target. However, the LED-optimized plan increased the maximum, mean, and minimum dose within the PTV by as much as 80 Gy, 11 Gy, and 3 Gy, respectively. Despite increased tumour dose levels, the 18 MV (3 × 1 cm(2)) arc plan improved or maintained the V20, V5, and mean lung dose metrics compared to the 6 MV (3 × 3 cm(2)) simulation. We conclude that LED-SBRT has the potential to increase dose gradients, and dose levels within a small lung tumour. The magnitude of tumour dose increase or lung sparing can be optimized through manipulation of RT parameters (e.g. beam energy and field size).

Dosimetric and radiobiological consequences of computed tomography-guided adaptive strategies for intensity modulated radiation therapy of the prostate.

PURPOSE: To examine a range of scenarios for image-guided adaptive radiation therapy of prostate cancer, including different schedules for megavoltage CT imaging, patient repositioning, and dose replanning. METHODS AND MATERIALS: We simulated multifraction dose distributions with deformable registration using 35 sets of megavoltage CT scans of 13 patients. We computed cumulative dose-volume histograms, from which tumor control probabilities and normal tissue complication probabilities (NTCPs) for rectum were calculated. Five-field intensity modulated radiation therapy (IMRT) with 18-MV x-rays was planned to achieve an isocentric dose of 76 Gy to the clinical target volume (CTV). The differences between D95, tumor control probability, V70Gy, and NTCP for rectum, for accumulated versus planned dose distributions, were compared for different target volume sizes, margins, and adaptive strategies. RESULTS: The CTV D95 for IMRT treatment plans, averaged over 13 patients, was 75.2 Gy. Using the largest CTV margins (10/7 mm), the D95 values accumulated over 35 fractions were within 2% of the planned value, regardless of the adaptive strategy used. For tighter margins (5 mm), the average D95 values dropped to approximately 73.0 Gy even with frequent repositioning, and daily replanning was necessary to correct this deficit. When personalized margins were applied to an adaptive CTV derived from the first 6 treatment fractions using the STAPLE (Simultaneous Truth and Performance Level Estimation) algorithm, target coverage could be maintained using a single replan 1 week into therapy. For all approaches, normal tissue parameters (rectum V(70Gy) and NTCP) remained within acceptable limits. CONCLUSIONS: The frequency of adaptive interventions depends on the size of the CTV combined with target margins used during IMRT optimization. The application of adaptive target margins (<5 mm) to an adaptive CTV determined 1 week into therapy minimizes the need for subsequent dose replanning.

Correction for 'artificial' electron disequilibrium due to cone-beam CT density errors: implications for on-line adaptive stereotactic body radiation therapy of lung.

Cone-beam computed tomography (CBCT) has rapidly become a clinically useful imaging modality for image-guided radiation therapy. Unfortunately, CBCT images of the thorax are susceptible to artefacts due to scattered photons, beam hardening, lag in data acquisition, and respiratory motion during a slow scan. These limitations cause dose errors when CBCT image data are used directly in dose computations for on-line, dose adaptive radiation therapy (DART). The purpose of this work is to assess the magnitude of errors in CBCT numbers (HU), and determine the resultant effects on derived tissue density and computed dose accuracy for stereotactic body radiation therapy (SBRT) of lung cancer. Planning CT (PCT) images of three lung patients were acquired using a Philips multi-slice helical CT simulator, while CBCT images were obtained with a Varian On-Board Imaging system. To account for erroneous CBCT data, three practical correction techniques were tested: (1) conversion of CBCT numbers to electron density using phantoms, (2) replacement of individual CBCT pixel values with bulk CT numbers, averaged from PCT images for tissue regions, and (3) limited replacement of CBCT lung pixels values (LCT) likely to produce artificial lateral electron disequilibrium. For each corrected CBCT data set, lung SBRT dose distributions were computed for a 6 MV volume modulated arc therapy (VMAT) technique within the Philips Pinnacle treatment planning system. The reference prescription dose was set such that 95% of the planning target volume (PTV) received at least 54 Gy (i.e. D95). Further, we used the relative depth dose factor as an a priori index to predict the effects of incorrect low tissue density on computed lung dose in regions of severe electron disequilibrium. CT number profiles from co-registered CBCT and PCT patient lung images revealed many reduced lung pixel values in CBCT data, with some pixels corresponding to vacuum (-1000 HU). Similarly, CBCT data in a plastic lung phantom were reduced by 200 HU compared with known CT number values. For the three patients, dose results using the CBCT number data registered with PCT showed a prescription dose reduction ranging from 4 to 13% (D95 = 47 Gy). Therefore, accurate determination of lung density, especially for very low lung density (<0.2 g cm(-3)) is essential, but difficult to achieve using the CBCT data. Applying corrective techniques (1) and (2) to CBCT patient data produced unacceptable dose differences. For one typical VMAT SBRT patient, the D95 for the corrected CBCT and BCT image-based plans differed by -4% (D95 = 52 Gy) and 9% (D95 = 59 Gy) compared to the co-registered PCT image-based plan. However, corrective technique (3) produced negligible dose differences comparing LCT and PCT image-based plans. With regard to implementing on-line DART, dose errors must be minimized because they affect re-optimization decisions, and prevent accurate accumulation of the dose distribution.

Categorizing segmentation quality using a quantitative quality assurance algorithm.

INTRODUCTION: Obtaining high levels of contouring consistency is a major limiting step in optimizing the radiotherapeutic ratio. We describe a novel quantitative methodology for the quality assurance (QA) of contour compliance referenced against a community set of contouring experts. METHODS: Two clinical tumour site scenarios (10 lung cases and one prostate case) were used with QA algorithm. For each case, multiple physicians (lung: n = 6, prostate: n = 25) segmented various target/organ at risk (OAR) structures to define a set of community reference contours. For each set of community contours, a consensus contour (Simultaneous Truth and Performance Level Estimation (STAPLE)) was created. Differences between each individual community contour versus the group consensus contour were quantified by consensus-based contouring penalty metric (PM) scores. New observers segmented these same cases to calculate individual PM scores (for each unique target/OAR) for each new observer-STAPLE pair for comparison against the community and consensus contours. RESULTS: Four physicians contoured the 10 lung cases for a total of 72 contours for quality assurance evaluation against the previously derived community consensus contours. A total of 16 outlier contours were identified by the QA system of which 11 outliers were due to over-contouring discrepancies, three were due to over-/under-contouring discrepancies, and two were due to missing/incorrect nodal contours. In the prostate scenario involving six physicians, the QA system detected a missing penile bulb contour, systematic inner-bladder contouring, and under-contouring of the upper/anterior rectum. CONCLUSION: A practical methodology for QA has been demonstrated with future clinical trial credentialing, medical education and auto-contouring assessment applications.