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.

[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.

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.

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.

Medical physics staffing for radiation oncology: a decade of experience in Ontario, Canada.

The January 2010 articles in The New York Times generated intense focus on patient safety in radiation treatment, with physics staffing identified frequently as a critical factor for consistent quality assurance. The purpose of this work is to review our experience with medical physics staffing, and to propose a transparent and flexible staffing algorithm for general use. Guided by documented times required per routine procedure, we have developed a robust algorithm to estimate physics staffing needs according to center-specific workload for medical physicists and associated support staff, in a manner we believe is adaptable to an evolving radiotherapy practice. We calculate requirements for each staffing type based on caseload, equipment inventory, quality assurance, educational programs, and administration. Average per-case staffing ratios were also determined for larger-scale human resource planning and used to model staffing needs for Ontario, Canada over the next 10 years. The workload specific algorithm was tested through a survey of Canadian cancer centers. For center-specific human resource planning, we propose a grid of coefficients addressing specific workload factors for each staff group. For larger scale forecasting of human resource requirements, values of 260, 700, 300, 600, 1200, and 2000 treated cases per full-time equivalent (FTE) were determined for medical physicists, physics assistants, dosimetrists, electronics technologists, mechanical technologists, and information technology specialists, respectively.

The effect of an inconsistent breathing amplitude on the relationship between an external marker and internal lung deformation in a porcine model.

PURPOSE: Investigate the relationship between the motion of the Varian Real-time Position Management (RPM) device and the internal motion of a pig during induced inconsistencies in the amplitude of breathing. METHODS: Twelve studies were performed on four ventilated female Landrace cross pigs using a GE Healthcare, Discovery CT 750 HD scanner. In each study, a 4.0 cm section (64 slices) of the pig's lungs was repeatedly scanned 20 times using cine mode, each time lasting more than one breathing cycle. During these cine scans, a Varian RPM device was used to collect respiratory amplitudes and the ventilator air return tube was periodically crimped to induce inconsistent breathing amplitudes. Each breathing cycle and its associated cine scan were categorized as either consistent or inconsistent, based on thresholds of the minimum expiration and maximum inspiration amplitudes. From the group of consistent amplitude cine scans in a study, a reference scan was chosen. The effect of inconsistent breathing amplitudes on the relationship between the motion of the RPM marker and the motion within three regions of interest (in each lung and the chest wall) was investigated with two methods: (1) A 4D-CT sorting algorithm based on RPM amplitude was used to sort volumes into 4D-CT phase bins. Within each phase bin, the nonlinear deformation of volumes collected during consistent and inconsistent breathing amplitude was calculated with respect to the reference volume. The magnitude of the deformations (in mm) were compared to determine if inconsistent breathing amplitude caused greater deformations. (2) Nonlinear deformations between each CT volume from a cine scan and the maximum expiration volume of the reference scan were calculated. Regression analyses between the nonlinear deformations within three regions of interest (in each lung and the chest wall) and the RPM amplitudes were performed to test the effect of inconsistent breathing amplitudes on the linearity of the relationship between the 3D motion of internal anatomy and the 1D motion of the RPM external marker. RESULTS: (1) Inconsistent versus consistent breathing amplitudes caused a significant increase in deformation relative to the reference scan within the left lung (1.40 +/- 0.42 versus 1.29 +/- 0.36 mm, p < 0.05). (2) One-to-one correspondences between motions of internal anatomies and motion of the RPM external marker did not exist. The regression lines between the two types of motions did not yield an identity relationship (unity slope and zero intercept). Inconsistent breathing produced significantly different regression lines than consistent breathing in ten of the 12 studies within a left lung region of interest. CONCLUSIONS: The results of these two studies indicate that inconsistency in the amplitude of breathing disrupted the correspondence between the motion of the external marker and internal anatomies. As a consequence, radiation therapy of tumors embedded in lung tissue may be prone to significant errors if inconsistent breathing amplitudes occur during treatment.

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.

Three-dimensional dosimetry of small megavoltage radiation fields using radiochromic gels and optical CT scanning.

The dosimetry of small fields as used in stereotactic radiotherapy, radiosurgery and intensity-modulated radiation therapy can be challenging and inaccurate due to partial volume averaging effects and possible disruption of charged particle equilibrium. Consequently, there exists a need for an integrating, tissue equivalent dosimeter with high spatial resolution to avoid perturbing the radiation beam and artificially broadening the measured beam penumbra. In this work, radiochromic ferrous xylenol-orange (FX) and leuco crystal violet (LCV) micelle gels were used to measure relative dose factors (RDFs), percent depth dose profiles and relative lateral beam profiles of 6 MV x-ray pencil beams of diameter 28.1, 9.8 and 4.9 mm. The pencil beams were produced via stereotactic collimators mounted on a Varian 2100 EX linear accelerator. The gels were read using optical computed tomography (CT). Data sets were compared quantitatively with dosimetric measurements made with radiographic (Kodak EDR2) and radiochromic (GAFChromic EBT) film, respectively. Using a fast cone-beam optical CT scanner (Vista), corrections for diffusion in the FX gel data yielded RDFs that were comparable to those obtained by minimally diffusing LCV gels. Considering EBT film-measured RDF data as reference, cone-beam CT-scanned LCV gel data, corrected for scattered stray light, were found to be in agreement within 0.5% and -0.6% for the 9.8 and 4.9 mm diameter fields, respectively. The validity of the scattered stray light correction was confirmed by general agreement with RDF data obtained from the same LCV gel read out with a laser CT scanner that is less prone to the acceptance of scattered stray light. Percent depth dose profiles and lateral beam profiles were found to agree within experimental error for the FX gel (corrected for diffusion), LCV gel (corrected for scattered stray light), and EBT and EDR2 films. The results from this study reveal that a three-dimensional dosimetry method utilizing optical CT-scanned radiochromic gels allows for the acquisition of a self-consistent volumetric data set in a single exposure, with sufficient spatial resolution to accurately characterize small fields.

The use of CT density changes at internal tissue interfaces to correlate internal organ motion with an external surrogate.

The purpose of this paper is to describe a non-invasive method to monitor the motion of internal organs affected by respiration without using external markers or spirometry, to test the correlation with external markers, and to calculate any time shift between the datasets. Ten lung cancer patients were CT scanned with a GE LightSpeed Plus 4-Slice CT scanner operating in a ciné mode. We retrospectively reconstructed the raw CT data to obtain consecutive 0.5 s reconstructions at 0.1 s intervals to increase image sampling. We defined regions of interest containing tissue interfaces, including tumour/lung interfaces that move due to breathing on multiple axial slices and measured the mean CT number versus respiratory phase. Tumour motion was directly correlated with external marker motion, acquired simultaneously, using the sample coefficient of determination, r(2). Only three of the ten patients showed correlation higher than r(2) = 0.80 between tumour motion and external marker position. However, after taking into account time shifts (ranging between 0 s and 0.4 s) between the two data sets, all ten patients showed correlation better than r(2) = 0.8. This non-invasive method for monitoring the motion of internal organs is an effective tool that can assess the use of external markers for 4D-CT imaging and respiratory-gated radiotherapy on a patient-specific basis.

A fully automated non-external marker 4D-CT sorting algorithm using a serial cine scanning protocol.

Current 4D-CT methods require external marker data to retrospectively sort image data and generate CT volumes. In this work we develop an automated 4D-CT sorting algorithm that performs without the aid of data collected from an external respiratory surrogate. The sorting algorithm requires an overlapping cine scan protocol. The overlapping protocol provides a spatial link between couch positions. Beginning with a starting scan position, images from the adjacent scan position (which spatial match the starting scan position) are selected by maximizing the normalized cross correlation (NCC) of the images at the overlapping slice position. The process was continued by 'daisy chaining' all couch positions using the selected images until an entire 3D volume was produced. The algorithm produced 16 phase volumes to complete a 4D-CT dataset. Additional 4D-CT datasets were also produced using external marker amplitude and phase angle sorting methods. The image quality of the volumes produced by the different methods was quantified by calculating the mean difference of the sorted overlapping slices from adjacent couch positions. The NCC sorted images showed a significant decrease in the mean difference (p < 0.01) for the five patients.

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.

Three-dimensional dose verification for intensity-modulated radiation therapy in the radiological physics centre head-and-neck phantom using optical computed tomography scans of ferrous xylenol-orange gel dosimeters.

PURPOSE: To extend the Radiological Physics Centre (RPC) intensity-modulated radiation therapy dose verification protocol to three dimensions using optical computed tomography (CT) scans of ferrous xylenol-orange (FX) gels. METHODS AND MATERIALS: The dosimetry insert in the RPC head-and-neck phantom was replaced with an FX cylindrical gel dosimeter. Two gels were calibrated, independently irradiated with 6-MV X-ray beams and scanned using laser and cone-beam (Vista) optical CT, respectively. For matching dose slices, measured dose distributions were compared with Pinnacle3 computed distributions. RESULTS: Within high-dose regions and low gradients, doses measured using laser CT were 2% to 3% less than the computed dose, whereas with cone-beam CT they were 4% to 5% less. Inside the central 90% of the gel cylinder diameter, the fraction of voxels satisfying the two-dimensional gamma analysis (5% dose difference, 3-mm distance to agreement) with laser-beam- and cone-beam-measured dose distributions were 98.4% and 99.0%, respectively. A three-dimensional gamma analysis with cone-beam data revealed that 96.7% of voxels within the central 90% gel volume satisfied the above criteria. Within the axial and sagittal planes through the primary planning target volume (PTV), computed and measured doses using GAFChromicEBT film (RPC measured) and cone-beam scanned FX gel generally agreed. At equivalent points in the planning target volumes, computed, thermoluminescent dosimeter (RPC-measured), and gel point doses agreed to within 5.1% in absolute dose. CONCLUSIONS: Laser and cone-beam CT yield comparable dose distributions in high-dose regions. The RPC head phantom and optical CT-scanned FX gels can be used for accurate intensity-modulated radiation therapy dose verification in three dimensions.

Experience-driven dose-volume histogram maps of NTCP risk as an aid for radiation treatment plan selection and optimization.

Commonly, the quality of treatment plans is judged by a dose-volume histogram (DVH) in regards to satisfying a series of dose-volume constraints. This paper presents a novel technique for mapping normal tissue complication probabilities (NTCP) onto regions of dose-volume space with statistical considerations of risk. Mapping is done for DVHs specific to one treatment technique for one disease site. In this study, the method is illustrated for simplified intensity modulated arc therapy of the prostate, and the resulting NTCP values apply to complications in the rectum. The method consists of implementing a Monte Carlo algorithm that creates a large set of DVH curves by simulating random walks through dose-volume space. The walks are guided by a base set of clinical DVHs. Grid points in the dose-volume space have an associated NTCP spectrum for curves passing above right of the grid point of interest. After a DVH is simulated and the NTCP estimate calculated using the Lyman model, dose-volume points located to the bottom left of the curve are scored with this NTCP value and contributed to the spectrum of each point. A NTCP tolerance value is then selected and the risk of violating this tolerance is identified by a gray-scale map in regions of dose-volume space. The generated maps distinguish technique-specific, high-risk regions, a feature which is advantageous over fixed single-point dose-volume constraints commonly used. The maps also provide a visualization tool to help select safe and robust treatment plans and open the possibility for improving the efficiency of biologically based plan optimization by focusing on the more critical sections of DVH curves.

An apparent threshold dose response in ferrous xylenol-orange gel dosimeters when scanned with a yellow light source.

Freshly prepared radiochromic ferrous xylenol-orange (FX) gels optically scanned with a light source exhibit a threshold dose response that is thermally and wavelength dependent. Correction for this threshold dose leads to accurate dose calibration and better reproducibility in multiple fraction radiation exposures. The objective of this study was to determine the cause of the threshold dose effect and to control it through improved dose calibration procedures. The results of a systematic investigation into the chemical cause revealed that impurities within the various FX gel constituents (i.e. xylenol-orange, gelatin, sulfuric acid and ferrous ammonium sulfate) were not directly responsible for the threshold dose. Rather, it was determined that the threshold dose response stems from a spectral sensitivity to different chemical complexes that are formed at different dose levels in FX gels between ferric (Fe(III)) ions and xylenol-orange (XO), i.e. Fe(III)i:XOj. A double Fe(III)2:XO1 complex preferentially absorbs at longer wavelengths (i.e. yellow), while at shorter wavelengths (i.e. green) the sensitivity is biased toward the single Fe(III)1:XO1 complex. As a result, when scanning with yellow light, freshly prepared FX gels require a minimum concentration of Fe(III) ions to shift the equilibrium concentration to favor the predominant production of the double Fe(III)2:XO1 complex at low doses. This can be accomplished via pre-irradiation of freshly prepared gels to a priming dose of approximately 0.5 Gy or allowing auto-oxidation to generate the startup concentration of Fe(III) ions required to negate the apparent threshold dose response.

Ionization chamber volume determination and quality assurance using micro-CT imaging.

Ion chamber dosimetry requires a high degree of precision, at all steps within the dosimetric process, in order to ensure accurate dose measurements. This work presents a novel technique for ion chamber volume determination and quality assurance, using micro-computed tomography (micro-CT). Four nominally identical Exradin A1SL chambers (0.056 cm(3)) (Standard Imaging, WI, USA) were imaged using a micro-CT system (GE Locus, GE Healthcare, London, Ontario) and irradiated in a 6 MV x-ray reference field. Air volumes were calculated from the CT datasets using 3D analysis software (Microview 2.1.1, General Electric Healthcare, London, Ontario). Differences in the volumes of each chamber determined using micro-CT images agreed with differences in the ionization response within 1% for each chamber. Calibration coefficients were then compared through cross-calibration with a calibrated ion chamber and from the CT-measured volumes. The average ratio of these values was found to be 0.958 +/- 0.009 indicating good correlation. The results demonstrate the promise of using micro-CT imaging for the absolute volumetric characterization of ion chambers. The images have the potential to be an important clinical tool for quality assurance of ion chamber construction and integrity after routine clinical usage.

A prospective evaluation of helical tomotherapy.

PURPOSE: To report results from two clinical trials evaluating helical tomotherapy (HT). METHODS AND MATERIALS: Patients were enrolled in one of two prospective trials of HT (one for palliative and one for radical treatment). Both an HT plan and a companion three-dimensional conformal radiotherapy (3D-CRT) plan were generated. Pretreatment megavoltage computed tomography was used for daily image guidance. RESULTS: From September 2004 to January 2006, a total of 61 sites in 60 patients were treated. In all but one case, a clinically acceptable tomotherapy plan for treatment was generated. Helical tomotherapy plans were subjectively equivalent or superior to 3D-CRT in 95% of plans. Helical tomotherapy was deemed equivalent or superior in two thirds of dose-volume point comparisons. In cases of inferiority, differences were either clinically insignificant and/or reflected deliberate tradeoffs to optimize the HT plan. Overall imaging and treatment time (median) was 27 min (range, 16-91 min). According to a patient questionnaire, 78% of patients were satisfied to very satisfied with the treatment process. CONCLUSIONS: Helical tomotherapy demonstrated clear advantages over conventional 3D-CRT in this diverse patient group. The prospective trials were helpful in deploying this technology in a busy clinical setting.

Dosimetric evaluation of daily rigid and nonrigid geometric correction strategies during on-line image-guided radiation therapy (IGRT) of prostate cancer.

The purpose of this study is to evaluate a geometric image guidance strategy that simultaneously correct for various inter-fractional rigid and nonrigid geometric uncertainties in an on-line environment, using field shape corrections (called the "MU-MLC" technique). The effectiveness of this strategy was compared with two other simpler on-line image guidance strategies that are more commonly used in the clinic. To this end, five prostate cancer patients, with at least 15 treatment CT studies each, were analyzed. The prescription dose was set to the maximum dose that did not violate the rectum and bladder dose-volume constraints, and hence, was unique to each patient. Deformable image registration and dose-tracking was performed on each CT image to obtain the cumulative treatment dose distributions. From this, maximum, minimum, and mean dose, as well as generalized equivalent uniform dose (gEUD) were calculated for each image guidance strategy. As expected, some dosimetric differences in the clinical target volume (CTV) were observed between the three image guidance strategies investigated. For example, up to +/-2% discrepancy in prostate minimum dose were observed among the techniques. Of them, only the "MU -MLC" technique did not reduce the prostate minimum dose for all patients (i.e., > or = 100%). However, the differences were clinically not significant to indicate the preference of one strategy over another, when using a uniform 5 mm margin size. For the organ-at-risks (OARs), the large rectum sparing effect (< or =5.7 Gy, gEUD) and bladder overdosing effect (< or = 16 Gy, gEUD) were observed. This was likely due to the use of bladder contrast during CT simulation studies which was not done during the treatment CT studies. Therefore, ultimately, strategies to maintain relatively constant rectum and bladder volumes, throughout the treatment course, are required to minimize this effect. In conclusion, the results here suggest that simple translational corrections based on three-dimensional (3D) images is adequate to maintain target coverage, for margin sizes at least as large as 5 mm. In addition, due to large fluctuations in OAR volumes, innovative image guidance strategies are needed to minimize dose and maintain consistent sparing during the whole course of radiation therapy.

Delineation of moving targets with slow MVCT scans: implications for adaptive non-gated lung tomotherapy.

Accurate imaging is a prerequisite for adaptive radiation therapy of mobile tumours. We present an evaluation of the performance of slow computed tomography (CT) for mapping and delineating the excursion boundary of a moving object using a tumour phantom scanned with the helical MVCT scanner of a tomotherapy unit. A spherical test object driven by sinusoidal motion in both the lateral and cranial-caudal directions was used to determine how well MVCT images depict the true envelope of the motion. Such information is useful in interpreting the CT images relative to the static object case when radiotherapy gating is to be used or in determining the internal target volume (ITV) when beam gating is not possible. A computer simulation of the CT imaging process was developed which incorporates the third generation fan beam geometry and helical acquisition technique of the tomotherapy MVCT system. Motion artefacts are mainly characterized by the parameter alpha=Tgantry/Trespiration which is interpreted as the period of the gantry rotation (Tgantry) in units of the respiratory period (Trespiration). Experimental tests were performed using a fixed gantry period of 10 s per full rotation and respiratory period ranging from 4.0 (alpha=2.5) to 1.0 (alpha=10) s. These cases represent typical clinical imaging conditions on the tomotherapy unit, as well as an extreme test case where the gantry period is intentionally set to be much greater than the respiratory period (termed an 'ultra-slow' scan). The accuracy of target (ITV) delineation is evaluated by comparing volumes generated using iso-density contours on the MVCT images to the true motion envelope, known a priori in this phantom study. As expected, motion artefacts are present in clinical MVCT images and they are not averaged over the slow gantry period of rotation. Furthermore, artefacts are not significantly affected by scanning with different helical pitch values. Greater distortions from the true density distribution are observed for lateral motion compared to cranial-caudal motion. Volumes generated by iso-density contours yield better agreement with the motion envelope for scans performed under ultra-slow conditions (alpha=10) compared to typical clinical imaging conditions (alpha=2.5). If the MVCT gantry cannot be rotated very quickly due to engineering constraints in order to achieve ultra-fast CT, we suggest an opposite approach as an interim measure for mapping the ITV. Adjusting MVCT scan conditions to a very slow acquisition (alpha=10) may be a good compromise for determining the ITV for non-gated adaptive tomotherapy of moving lung tumours.