Dosimetric evaluation of irradiation geometry and potential air gaps in an acrylic miniphantom used for external audit of absolute dose calibration for a hybrid 1.5 T MR-linac system.

INTRODUCTION: To investigate the impact of partial lateral scatter (LS), backscatter (BS) and presence of air gaps on optically stimulated luminescence dosimeter (OSLD) measurements in an acrylic miniphantom used for dosimetry audit on the 1.5 T magnetic resonance-linear accelerator (MR-linac) system. METHODS: The following irradiation geometries were investigated using OSLDs, A26 MR/A12 MR ion chamber (IC), and Monaco Monte Carlo system: (a) IC/OSLD in an acrylic miniphantom (partial LS, partial BS), (b) IC/OSLD in a miniphantom placed on a solid water (SW) stack at a depth of 1.5 cm (partial LS, full BS), (c) IC/OSLD placed at a depth of 1.5 cm inside a 3 cm slab of SW/buildup (full LS, partial BS), and (d) IC/OSLD centered inside a 3 cm slab of SW/buildup at a depth of 1.5 cm placed on top of a SW stack (full LS, full BS). Average of two irradiated OSLDs with and without water was used at each setup. An air gap of 1 and 2 mm, mimicking presence of potential air gap around the OSLDs in the miniphantom geometry was also simulated. The calibration condition of the machine was 1 cGy/MU at SAD = 143.5 cm, d = 5 cm, G90, and 10 × 10 cm2 . RESULTS: The Monaco calculation (0.5% uncertainty and 1.0 mm voxel size) for the four setups at the measurement point were 108.2, 108.1, 109.4, and 110.0 cGy. The corresponding IC measurements were 109.0 ± 0.03, 109.5 ± 0.06, 110.2 ± 0.02, and 109.8 ± 0.03 cGy. Without water, OSLDs measurements were ∼10% higher than the expected. With added water to minimize air gaps, the measurements were significantly improved to within 2.2%. The dosimetric impacts of 1 and 2 mm air gaps were also verified with Monaco to be 13.3% and 27.9% higher, respectively, due to the electron return effect. CONCLUSIONS: A minimal amount of air around or within the OSLDs can cause measurement discrepancies of 10% or higher when placed in a high b-field MR-linac system. Care must be taken to eliminate the air from within and around the OSLD.

Efficiency and safety increases after the implementation of a multi-institutional automated plan check tool at our institution.

PURPOSE: The plan check tool (PCT) is the result of a multi-institutional collaboration to jointly develop a flexible automated plan checking framework designed with the versatility to be shared across collaborating facilities while supporting the individual differences between practices. We analyze the effect that PCT has had on the efficiency and effectiveness of initial chart checks at our institution. METHODS AND MATERIALS: Data on errors identified during initial chart checks were acquired during two time periods: before the introduction of PCT in the clinic (6/24/2015 to 7/31/2015, 187 checks) and post-clinical release (4/14/2016 to 5/2/2016, 186 checks). During each time period, human plan checkers were asked to record all issues that they either manually detected or that were detected by PCT as well as the amount of time, less breaks, or interruptions, it took to check each plan. RESULTS: After the clinical release of PCT, there was a statistically significant decrease in the number of issues recorded by the human plan checkers both related to checks explicitly performed by PCT (13 vs 50, P < 0.001) and in issues identified overall (127 vs 200, P < 0.001). The mean and medium time for a plan check decreased by 20%. CONCLUSIONS: The use of a multi-institutional, configurable, automated plan checking tool has resulted in both substantial gains in efficiency and moving error detection to earlier points in the planning process, decreasing their likelihood that they reach the patient. The sizeable startup effort needed to create this tool from scratch was mitigated by the sharing, and subsequent co-development, of software code from a peer institution.

A phantom study to evaluate three different registration platform of 3D/3D, 2D/3D, and 3D surface match with 6D alignment for precise image-guided radiotherapy.

PURPOSE: To evaluate two three-dimensional (3D)/3D registration platforms, one two-dimensional (2D)/3D registration method, and one 3D surface registration method (3DS). These three technologies are available to perform six-dimensional (6D) registrations for image-guided radiotherapy treatment. METHODS: Fiducial markers were asymmetrically placed on the surfaces of an anthropomorphic head phantom (n = 13) and a body phantom (n = 8), respectively. The point match (PM) solution to the six-dimensional (6D) transformation between the two image sets [planning computed tomography (CT) and cone beam CT (CBCT)] was determined through least-square fitting of the fiducial positions using singular value decomposition (SVD). The transformation result from SVD was verified and was used as the gold standard to evaluate the 6D accuracy of 3D/3D registration in Varian's platform (3D3DV), 3D/3D and 2D/3D registration in the BrainLab ExacTrac system (3D3DE and 2D3D), as well as 3DS in the AlignRT system. Image registration accuracy from each method was quantitatively evaluated by root mean square of target registration error (rmsTRE) on fiducial markers and by isocenter registration error (IRE). The Wilcoxon signed-rank test was utilized to compare the difference of each registration method with PM. A P < 0.05 was considered significant. RESULTS: rmsTRE was in the range of 0.4 mm/0.7 mm (cranial/body), 0.5 mm/1 mm, 1.0 mm/1.5 mm, and 1.0 mm/1.2 mm for PM, 3D3D, 2D3D, and 3DS, respectively. Comparing to PM, the mean errors of IRE were 0.3 mm/1 mm for 3D3D, 0.5 mm/1.4 mm for 2D3D, and 1.6 mm/1.35 mm for 3DS for the cranial and body phantoms respectively. Both of 3D3D and 2D3D methods differed significantly in the roll direction as compared to the PM method for the cranial phantom. The 3DS method was significantly different from the PM method in all three translation dimensions for both the cranial (P = 0.003-P = 0.03) and body (P < 0.001-P = 0.008) phantoms. CONCLUSION: 3D3D using CBCT had the best image registration accuracy among all the tested methods. 2D3D method was slightly inferior to the 3D3D method but was still acceptable as a treatment position verification device. 3DS is comparable to 2D3D technique and could be a substitute for X-ray or CBCT for pretreatment verification for treatment of anatomical sites that are rigid.

Impact of varying air cavity on planning dosimetry for rectum patients treated on a 1.5 T hybrid MR-linac system.

PURPOSE: To investigate the dosimetric impact of magnetic (B) field on varying air cavities in rectum patients treated on the hybrid 1.5 T MR-linac. METHODS: Artificial air cavities of varying diameters (0.0, 1.0, 1.5, 2.0, 2.5, 3.0, and 5.0 cm) were created for four rectum patients (two prone and two supine). A total of 56 plans using a 7 MV flattening filter-free beam were generated with and without B-field. Reference intensity-modulated radiation therapy treatment plans without air cavity in the presence and absence of B-field were generated to a total dose of 45/50 Gy. The reference plans were copied and recalculated for the varying air cavities. D95 (PTV45 -PTV50 ), D95 (PTV50 -aircavity), V50 (PTV50 -aircavity), Dmax (PTV50 -aircavity), and V110% (PTV50 -aircavity) were extracted for each patient. Annulus rings of 1-mm-diameter step size were generated for one of the air cavity plans (3.0 cm) for all four patients to determine Dmax (%) and V110% (cc) within each annulus. RESULTS: In the presence of B-field, hot spots at the cavity interface start to become visible at ~1 cm air cavity in both supine and prone positioning due to electron return effect (ERE). In the presence of B-field Dmax and V110% varied from 5523 ± 49 cGy and 0.09 ± 0.16 cc for 0 cm air cavity size to 6050 ± 109 cGy and 11.6 ± 6.7 cc for 5 cm air cavity size. The hot spots were located within 3 mm inside the rectal-air interface, where Dmax increased from 110.4 ± 0.5% without B-field to 119.2 ± 0.8 % with B-field. CONCLUSIONS: Air cavities inside rectum affects rectum plan dosimetry due ERE. Location and magnitude of hot spots are dependent on the size of the air cavity.

Dosimetric evaluation of an atlas-based synthetic CT generation approach for MR-only radiotherapy of pelvis anatomy.

PURPOSE: To investigate the potential of an atlas-based approach in generation of synthetic CT for pelvis anatomy. METHODS: Twenty-three matched pairs of computed tomography (CT) and magnetic resonance imaging (MRI) scans were selected from a pool of prostate cancer patients. All MR scans were preprocessed to reduce scanner- and patient-induced intensity inhomogeneities and to standardize their intensity histograms. Ten (training dataset) of 23 pairs were then utilized to construct the coregistered CT-MR atlas. The synthetic CT for a new patient is generated by appropriately weighting the deformed atlas of CT-MR onto the new patient MRI. The training dataset was used as an atlas to generate the synthetic CT for the rest of the patients (test dataset). The mean absolute error (MAE) between the deformed planning CT and synthetic CT was computed over the entire CT image, bone, fat, and muscle tissues. The original treatment plans were also recomputed on the new synthetic CTs and dose-volume histogram metrics were compared. The results were compared with a commercially available synthetic CT Software (MRCAT) that is routinely used in our clinic. RESULTS: MAE errors (±SD) between the deformed planning CT and our proposed synthetic CTs in the test dataset were 47 ± 5, 116 ± 12, 36 ± 6, and 47 ± 5 HU for the entire image, bone, fat, and muscle tissues respectively. The MAEs were 65 ± 5, 172 ± 9, 43 ± 7, and 42 ± 4 HU for the corresponding tissues in MRCAT CT. The dosimetric comparison showed consistent results for all plans using our synthetic CT, deformed planning CT and MRCAT CT. CONCLUSION: We investigated the potential of a multiatlas approach to generate synthetic CT images for the pelvis. Our results demonstrate excellent results in terms of HU value assignment compared to the original CT and dosimetric consistency.

Dual-input tracer kinetic modeling of dynamic contrast-enhanced MRI in thoracic malignancies.

Pulmonary perfusion with dynamic contrast-enhanced (DCE-) MRI is typically assessed using a single-input tracer kinetic model. Preliminary studies based on perfusion CT are indicating that dual-input perfusion modeling of lung tumors may be clinically valuable as lung tumors have a dual blood supply from the pulmonary and aortic system. This study aimed to investigate the feasibility of fitting dual-input tracer kinetic models to DCE-MRI datasets of thoracic malignancies, including malignant pleural mesothelioma (MPM) and nonsmall cell lung cancer (NSCLC), by comparing them to single-input (pulmonary or systemic arterial input) tracer kinetic models for the voxel-level analysis within the tumor with respect to goodness-of-fit statistics. Fifteen patients (five MPM, ten NSCLC) underwent DCE-MRI prior to radiotherapy. DCE-MRI data were analyzed using five different single- or dual-input tracer kinetic models: Tofts-Kety (TK), extended TK (ETK), two compartment exchange (2CX), adiabatic approximation to the tissue homogeneity (AATH) and distributed parameter (DP) models. The pulmonary blood flow (BF), blood volume (BV), mean transit time (MTT), permeability-surface area product (PS), fractional interstitial volume (vI ), and volume transfer constant (KTrans ) were calculated for both single- and dual-input models. The pulmonary arterial flow fraction (γ), pulmonary arterial blood flow (BFPA ) and systemic arterial blood flow (BFA ) were additionally calculated for only dual-input models. The competing models were ranked and their Akaike weights were calculated for each voxel according to corrected Akaike information criterion (cAIC). The optimal model was chosen based on the lowest cAIC value. In both types of tumors, all five dual-input models yielded lower cAIC values than their corresponding single-input models. The 2CX model was the best-fitted model and most optimal in describing tracer kinetic behavior to assess microvascular properties in both MPM and NSCLC. The dual-input 2CX-model-derived BFA was the most significant parameter in differentiating adenocarcinoma from squamous cell carcinoma histology for NSCLC patients.

Diffusion-weighted MRI of the lung at 3T evaluated using echo-planar-based and single-shot turbo spin-echo-based acquisition techniques for radiotherapy applications.

PURPOSE: To compare single-shot echo-planar (SS-EPI)-based and turbo spin-echo (SS-TSE)-based diffusion-weighted imaging (DWI) in Non-Small Cell Lung Cancer (NSCLC) patients and to characterize the distributions of apparent diffusion coefficient (ADC) values generated by the two techniques. METHODS: Ten NSCLC patients were enrolled in a prospective IRB-approved study to compare and optimize DWI using EPI and TSE-based techniques for radiotherapy planning. The imaging protocol included axial T2w, EPI-based DWI and TSE-based DWI on a 3 T Philips scanner. Both EPI-based and TSE-based DWI sequences used three b values (0, 400, and 800 s/mm2 ). The acquisition times for EPI-based and TSE-based DWI were 5 and 8 min, respectively. DW-MR images were manually coregistered with axial T2w images, and tumor volume contoured on T2w images were mapped onto the DWI scans. A pixel-by-pixel fit of tumor ADC was calculated based on monoexponential signal behavior. Tumor ADC mean, standard deviation, kurtosis, and skewness were calculated and compared between EPI and TSE-based DWI. Image distortion and ADC values between the two techniques were also quantified using fieldmap analysis and a NIST traceable ice-water diffusion phantom, respectively. RESULTS: The mean ADC for EPI and TSE-based DWI were 1.282 ± 0.42 × 10-3 and 1.211 ± 0.31 × 10-3  mm2 /s. The average skewness and kurtosis were 0.14 ± 0.4 and 2.43 ± 0.40 for DWI-EPI and -0.06 ± 0.69 and 2.89 ± 0.62 for DWI-TSE. Fieldmap analysis showed a mean distortion of 13.72 ± 8.12 mm for GTV for DWI-EPI and 0.61 ± 0.4 mm for DWI-TSE. ADC values obtained using the diffusion phantom for the two techniques were within 0.03 × 10-3  mm2 /s with respect to each other as well as the established values. CONCLUSIONS: Diffusion-weighted turbo spin-echo shows better geometrical accuracy compared to DWI-EPI. Mean ADC values were similar with both acquisitions but the shape of the histograms was different based on the skewness and kurtosis values. The impact of differences in respiratory technique on ADC values requires further investigation.

Developing a MLC modifier program to improve fiducial detection for MV/kV imaging during hypofractionated prostate volumetric modulated arc therapy.

PURPOSE: To develop an Eclipse plug-in (MLC_MODIFIER) that automatically modifies control points to expose fiducials obscured by MLC during VMAT, thereby facilitating tracking using periodic MV/kV imaging. METHOD: Three-dimensional fiducial tracking was performed during VMAT by pairing short-arc (3°) MV digital tomosynthesis (DTS) images to triggered kV images. To evaluate MLC_MODIFIER efficacy, two cohorts of patients were considered. For first 12 patients, plans were manually edited to expose one fiducial marker. Next for 15 patients, plans were modified using MLC_MODIFIER script. MLC_MODIFIER evaluated MLC apertures at appropriate angles for marker visibility. Angles subtended by control points were compressed and low-dose "imaging" control points were inserted and exposed one marker with 1 cm margin. Patient's images were retrospectively reviewed to determine rate of MV registration failures. Failure categories were poor DTS image quality, MLC blockage of fiducials, or unknown reasons. Dosimetric differences in rectum, bladder, and urethra D1 cc, PTV maximum dose, and PTV dose homogeneity (PTV HI) were evaluated. Statistical significance was evaluated using Fisher's exact and Student's t test. RESULT: Overall MV registration failures, failures due to poor image quality, MLC blockage, and unknown reasons were 33% versus 8.9% (P < 0.0001), 8% versus 6.4% (P < 0.05), 13.6% versus 0.1% (P < 0.0001), and 7.6% versus 2.4% (P < 0.0001) for manually edited and MLC_MODIFIER plans, respectively. PTV maximum and HI increased on average from unmodified plans by 2.1% and 0.3% (P < 0.004) and 22.0% and 3.3% (P < 0.004) for manually edited and MLC_MODIFIED plans, respectively. Changes in bladder, rectum, and urethra D1CC were similar for each method and less than 0.7%. CONCLUSION: Increasing fiducial visibility via an automated process comprised of angular compression of control points and insertion of additional "imaging" control points is feasible. Degradation of plan quality is minimal. Fiducial detection and registration success rates are significantly improved compared to manually edited apertures.

Impact of daily soft-tissue image guidance to prostate on pelvic lymph node (PLN) irradiation for prostate patients receiving SBRT.

PURPOSE: To determine the impact of using fiducial match for daily image-guidance on pelvic lymph node (PLN) coverage for prostate cancer patients receiving stereotactic body radiation therapy (SBRT). METHODS: Thirty patients underwent SBRT treatment to the prostate and PLN from 2014 to 2016. Each patient received either 800cGy × 5 or 500cGy × 5 to the prostate and 500cGy × 5 to the PLN. A 5 mm clinical target volume (CTV)-to-planning target volume (PTV) margin around the PLN was used for planning. Two registrations with planning computed tomography (PCT) for each of the daily cone beam CTs (CBCTs) were performed: a rigid registration to fiducials and to the bony anatomy. The average translational difference between fiducial and bony match as well as percentage of fractions with differences > 5mm were calculated. Changes in bladder and rectal volume as well as center-of-mass (COM) position from simulation parameters, and their correlation with translational difference were also evaluated. The dosimetric impact of the translational differences was calculated by shifting the plan isocenter. RESULTS: The average translational difference between fiducial and bony match was 0.06 ± 0.82, 2.1 ± 4.1, -2.8 ± 4.3, and 5.5 ± 4.2 mm for lateral, vertical, longitudinal, and vector directions. The average change in bladder and rectal volume from simulation was -67.2 ± 163.04 cc (-12 ± 52%) and -1.6 ± 18.75 (-2 ± 30%) cc. The average change in COM of bladder from the simulation position was 0.34 ± 2.49, 4.4 ± 8.1, and -3.9 ± 7.5 mm along the LR, AP, and SI directions. The corresponding COM change for the rectum was 0.17 ± 1.9, 1.34 ± 3.5, and -0.6 ± 5.2 mm. CONCLUSIONS: The 5 mm margin covered ~75% of fractions receiving PLN irradiation with SBRT, daily CBCT and fiducial-guided setup. The dosimetric impact on PLN coverage was significant in 19% of fractions or 25% of patients. A larger translational shift was due to variation in rectal volume and changes in COM position of the bladder and rectum. A consistent bladder positioning and/or rectum filling compared with presimulation volume were essential for adequate coverage of PLN in a hypofractionated treatment regime.

Clinical evaluation of 4D MRI in the delineation of gross and internal tumor volumes in comparison with 4DCT.

PURPOSE: To evaluate clinical utility of respiratory-correlated (RC) four-dimensional magnetic resonance imaging (4DMRI) for lung tumor delineation and motion assessment, in comparison with the current clinical standard of 4D computed tomography (4DCT). METHODS AND MATERIALS: A prospective T2-weighted (T2w) RC-4DMRI technique was applied to acquire coronal 4DMRI images for 14 lung cancer patients (16 lesions) during free breathing (FB) under an IRB-approved protocol, together with a breath-hold (BH) T1w 3DMRI and axial 4DMRI. Clinical simulation CT and 4DCT were acquired within 2 h. An internal navigator was applied to trigger amplitude-binned 4DMRI acquisition whereas a bellows or real-time position management (RPM) was used in the 4DCT reconstruction. Six radiation oncologists manually delineated the gross and internal tumor volumes (GTV and ITV) in 399 3D images using programmed clinical workflows under a tumor delineation guideline. The ITV was the union of GTVs within the breathing cycle without margin. Average GTV and motion range were assessed and ITV variation between 4DMRI and 4DCT was evaluated using the Dice similarity index, mean distance agreement (MDA), and volume difference. RESULTS: The mean tumor volume is similar between 4DCT (GTV4DCT  = 1.0, as the reference) and T2w-4DMRI (GTVT2wMR  = 0.97), but smaller in T1w MRI (GTVT1wMR  = 0.76), suggesting possible peripheral edema around the tumor. Average GTV variation within the breathing cycle (22%) in 4DMRI is slightly greater than 4DCT (17%). GTV motion variation (-4 to 12 mm) and ITV variation (∆VITV =-25 to 95%) between 4DCT and 4DMRI are large, confirmed by relatively low ITV similarity (Dice = 0.72 ± 0.11) and large MDA = 2.9 ± 1.5 mm. CONCLUSION: Average GTVs are similar between T2w-4DMRI and 4DCT, but smaller by 25% in T1w BH MRI. Physician training and breathing coaching may be necessary to reduce ITV variability between 4DMRI and 4DCT. Four-dimensional magnetic resonance imaging is a promising and viable technique for clinical lung tumor delineation and motion assessment.

Institutional experience with SRS VMAT planning for multiple cranial metastases.

BACKGROUND AND PURPOSE: This study summarizes the cranial stereotactic radiosurgery (SRS) volumetric modulated arc therapy (VMAT) procedure at our institution. MATERIALS AND METHODS: Volumetric modulated arc therapy plans were generated for 40 patients with 188 lesions (range 2-8, median 5) in Eclipse and treated on a TrueBeam STx. Limitations of the custom beam model outside the central 2.5 mm leaves necessitated more than one isocenter pending the spatial distribution of lesions. Two to nine arcs were used per isocenter. Conformity index (CI), gradient index (GI) and target dose heterogeneity index (HI) were determined for each lesion. Dose to critical structures and treatment times are reported. RESULTS: Lesion size ranged 0.05-17.74 cm3 (median 0.77 cm3 ), and total tumor volume per case ranged 1.09-26.95 cm3 (median 7.11 cm3 ). For each lesion, HI ranged 1.2-1.5 (median 1.3), CI ranged 1.0-2.9 (median 1.2), and GI ranged 2.5-8.4 (median 4.4). By correlating GI to PTV volume a predicted GI = 4/PTV0.2 was determined and implemented in a script in Eclipse and used for plan evaluation. Brain volume receiving 7 Gy (V7 Gy ) ranged 10-136 cm3 (median 42 cm3 ). Total treatment time ranged 24-138 min (median 61 min). CONCLUSIONS: Volumetric modulated arc therapy provide plans with steep dose gradients around the targets and low dose to critical structures, and VMAT treatment is delivered in a shorter time than conventional methods using one isocenter per lesion. To further improve VMAT planning for multiple cranial metastases, better tools to shorten planning time are needed. The most significant improvement would come from better dose modeling in Eclipse, possibly by allowing for customizing the dynamic leaf gap (DLG) for a special SRS model and not limit to one DLG per energy per treatment machine and thereby remove the limitation on the Y-jaw and allow planning with a single isocenter.

Dose calculation for hypofractionated volumetric-modulated arc therapy: approximating continuous arc delivery and tongue-and-groove modeling.

Hypofractionated treatments generally increase the complexity of a treatment plan due to the more stringent constraints of normal tissues and target coverage. As a result, treatment plans contain more modulated MLC motions that may require extra efforts for accurate dose calculation. This study explores methods to minimize the differences between in-house dose calculation and actual delivery of hypofractionated volumetric-modulated arc therapy (VMAT), by focusing on arc approximation and tongue-and-groove (TG) modeling. For dose calculation, the continuous delivery arc is typically approximated by a series of static beams with an angular spacing of 2°. This causes significant error when there is large MLC movement from one beam to the next. While increasing the number of beams will minimize the dose error, calculation time will increase significantly. We propose a solution by inserting two additional apertures at each of the beam angle for dose calculation. These additional apertures were interpolated at two-thirds' degree before and after each beam. Effectively, there were a total of three MLC apertures at each beam angle, and the weighted average fluence from the three apertures was used for calculation. Because the number of beams was kept the same, calculation time was only increased by about 6%-8%. For a lung plan, areas of high local dose differences (> 4%) between film measurement and calculation with one aperture were significantly reduced in calculation with three apertures. Ion chamber measurement also showed similar results, where improvements were seen with calculations using additional apertures. Dose calculation accuracy was further improved for TG modeling by developing a sampling method for beam fluence matrix. Single element point sampling for fluence transmitted through MLC was used for our fluence matrix with 1 mm resolution. For Varian HDMLC, grid alignment can cause fluence sampling error. To correct this, transmission volume averaging was applied. For three paraspinal HDMLC cases, the average dose difference was greatly reduced in film and calculation comparisons with our new approach. The gamma (3%, 3 mm) pass rates have improved significantly from 74.1%, 90.0%, and 90.4% to 99.2%, 97.9%, and 97.3% for three cases, for calculation without volume averaging and calculation with volume averaging, respectively. Our results indicate that more accurate MLC leaf position and transmission sampling can improve accuracy and agreement between calculation and measurement, and are particularly important for hypofractionated VMAT that consists of large MLC movement.

Simultaneous MV-kV imaging for intrafractional motion management during volumetric-modulated arc therapy delivery.

The purpose of this study was to evaluate the accuracy and clinical feasibility of a motion monitoring method employing simultaneously acquired MV and kV images during volumetric-modulated arc therapy (VMAT). Short-arc digital tomosynthesis (SA-DTS) is used to improve the quality of the MV images that are then combined with orthogonally acquired kV images to assess 3D motion. An anthropomorphic phantom with implanted gold seeds was used to assess accuracy of the method under static, typical prostatic, and respiratory motion scenarios. Automatic registra-tion of kV images and single MV frames or MV SA-DTS reconstructed with arc lengths from 2° to 7° with the appropriate reference fiducial template images was performed using special purpose-built software. Clinical feasibility was evaluated by retrospectively analyzing images acquired over four or five sessions for each of three patients undergoing hypofractionated prostate radiotherapy. The standard deviation of the registration error in phantom using MV SA-DTS was similar to single MV images for the static and prostate motion scenarios (σ = 0.25 mm). Under respiratory motion conditions, the standard deviation of the registration error increased to 0.7mm and 1.7 mm for single MV and MV SA-DTS, respectively. Registration failures were observed with the respiratory scenario only and were due to motion-induced fiducial blurring. For the three patients studied, the mean and standard deviation of the difference between automatic registration using 4° MV SA-DTS and manual registration using single MV images results was 0.07±0.52mm. The MV SA-DTS results in patients were, on average, superior to single-frame MV by nearly 1 mm - significantly more than what was observed in phantom. The best MV SA-DTS results were observed with arc lengths of 3° to 4°. Registration failures in patients using MV SA-DTS were primarily due to blockage of the gold seeds by the MLC. The failure rate varied from 2% to 16%. Combined MV SA-DTS and kV imaging is feasible for intratreatment motion monitoring during VMAT of anatomic sites where limited motion is expected, and improves registration accuracy compared to single MV/kV frames. To create a clinically robust technique, further improvements to ensure visualization of fiducials at the desired control points without degradation of the treatment plan are needed.

Clinical experience with two frameless stereotactic radiosurgery (fSRS) systems using optical surface imaging for motion monitoring.

The purpose of this study was to compare two clinical immobilization systems for intracranial frameless stereotactic radiosurgery (fSRS) under the same clinical procedure using cone-beam computed tomography (CBCT) for setup and video-based optical surface imaging (OSI) for initial head alignment and intrafractional motion monitoring. A previously established fSRS procedure was applied using two intracranial immobilization systems: PinPoint system (head mold and mouthpiece) and Freedom system (head mold and open face mask). The CBCT was used for patient setup with four degrees of freedom (4DOF), while OSI was used for 6DOF alignment prior to CBCT, post-CBCT setup verification at all treatment couch angles (zero and nonzero), and intrafractional motion monitoring. Quantitative comparison of the two systems includes residual head rotation, head restriction capacity, and patient setup time in 25 patients (29 lesions) using PinPoint and 8 patients (29fractions) using Freedom. The maximum possible motion was assessed in nine volunteers with deliberate, forced movement in Freedom system. A consensus-based comparison of patient comfort level and clinical ease of use is reported. Using OSI-guided corrections, the maximum residual rotations in all directions were 1.1° ± 0.5° for PinPoint and 0.6° ± 0.3° for Freedom. The time spent performing rotation corrections was 5.0 ± 4.1 min by moving the patient with PinPoint and 2.7 ± 1.0min by adjusting Freedom couch extension. After CBCT, the OSI-CBCT discrepancy due to different anatomic landmarks for alignment was 2.4 ± 1.3 mm using PinPoint and 1.5 ± 0.7 mm using Freedom. Similar results were obtained for setup verification at couch angles (< 1.5 mm) and for motion restriction: 0.4± 0.3 mm/0.2° ± 0.2° in PinPoint and 0.6 ± 0.3 mm/0.3° ± 0.1° in Freedom. The maximum range of forced head motion was 2.2 ± 1.0 mm using Freedom. Both intracranial fSRS immobilization systems can restrict head motion within 1.5 mm during treatment as monitored by OSI. Setting a motion threshold for beam-hold ensures that head motion is constrained within the treatment margin during beam-on periods. The capability of 6D setup is useful to improve treatment accuracy. Patient comfort and clinical workflow should play a substantial role in system selection, and Freedom system outperforms PinPoint system in these two aspects.

Incorporation of treatment plan spatial and temporal dose patterns into a prostate intrafractional motion management strategy.

PURPOSE: Periodic MV∕KV radiographs taken during volumetric modulated arc therapy (VMAT) for hypofractionated treatment provide guidance in intrafractional motion management. The choice of imaging frequency and timing are key components in delivering the desired dose while reducing associated overhead such as imaging dose, preparation, and processing time. In this project the authors propose a paradigm with imaging timing and frequency based on the spatial and temporal dose patterns of the treatment plan. METHODS: A number of control points are used in treatment planning to model VMAT delivery. For each control point, the sensitivity of individual target or organ-at-risk dose to motion can be calculated as the summation of dose degradations given the organ displacements along a number of possible motion directions. Instead of acquiring radiographs at uniform time intervals, MV∕KV image pairs are acquired indexed to motion sensitivity. Five prostate patients treated via hypofractionated VMAT are included in this study. Intrafractional prostate motion traces from the database of an electromagnetic tracking system are used to retrospectively simulate the VMAT delivery and motion management. During VMAT delivery simulation patient position is corrected based on the radiographic findings via couch movement if target deviation violates a patient-specific 3D threshold. The violation rate calculated as the percentage of traces failing the clinical dose objectives after motion correction is used to evaluate the efficacy of this approach. RESULTS: Imaging indexed to a 10 s equitime interval and correcting patient position accordingly reduces the violation rate to 19.5% with intervention from 44.5% without intervention. Imaging indexed to the motion sensitivity further reduces the violation rate to 12.1% with the same number of images. To achieve the same 5% violation rate, the imaging incidence can be reduced by 40% by imaging indexed to motion sensitivity instead of time. CONCLUSIONS: The simulation results suggest that image scheduling according to the characteristics of the treatment plan can improve the efficiency of intrafractional motion management. Using such a technique, the accuracy of delivered dose during image-guided hypofractionated VMAT treatment can be improved.

Impact of a Centralized Database System on Radiation Therapy Quality Assurance Management at a Large Health Care Network: 5 Years' Experience.

PURPOSE: This study reports the impact of using a centralized database system for major equipment quality assurance (QA) at a large institution. METHODS AND MATERIALS: A centralized database system has been implemented for radiation therapy machine QA in our institution at 6 campuses with 11 computed tomographies and 22 linear accelerators (LINACs). The database system was customized to manage monthly and annual computed tomography and LINAC QA. This includes providing the same set of QA procedures across the enterprise, digitally storing all measurement records, and generating trend analyses. Compared with conventional methods (ie, paper forms), the effectiveness of the database system was quantified by changes in the compliance of QA tests and perceptions of staff to the efficiency of data retrieval and analyses. An anonymized questionnaire was provided to physicists enterprise-wide to assess workflow changes. RESULTS: With the implementation of the database system, the compliance of QA test completion improved from 80% to >99% for the entire institution. This resonates with the 56% of physicists who found the database system helpful in guiding them through QA, and 25% of physicists found the contrary, and 19% reported no difference (n = 16). Meanwhile, 40% of physicists reported longer times needed to record data using the database system compared with conventional methods, and another 40% suggested otherwise. In addition, 87% and 80% found the database more efficient to analyze and retrieve previous data, respectively. This was also reflected by the shorter time taken to generate year-end QA statistics using the software (5 vs 30 min per LINAC). Overall, 94% of physicists preferred the centralized database system over conventional methods and endorsed continued use of the system. CONCLUSIONS: A centralized database system is useful and can improve the effectiveness and efficiency of QA management in a large institution. With consistent data collection and proper data storage using a database, high-quality data can be obtained for failure modes and effects analyses as per TG 100.

Automated and Clinically Optimal Treatment Planning for Cancer Radiotherapy.

Each year, approximately 18 million new cancer cases are diagnosed worldwide, and about half must be treated with radiotherapy. A successful treatment requires treatment planning with the customization of penetrating radiation beams to sterilize cancerous cells without harming nearby normal organs and tissues. This process currently involves extensive manual tuning of parameters by an expert planner, making it a time-consuming and labor-intensive process, with quality and immediacy of critical care dependent on the planner's expertise. To improve the speed, quality, and availability of this highly specialized care, Memorial Sloan Kettering Cancer Center developed and applied advanced optimization tools to this problem (e.g., using hierarchical constrained optimization, convex approximations, and Lagrangian methods). This resulted in both a greatly improved radiotherapy treatment planning process and the generation of reliable and consistent high-quality plans that reflect clinical priorities. These improved techniques have been the foundation of high-quality treatments and have positively impacted over 4,000 patients to date, including numerous patients in severe pain and in urgent need of treatment who might have otherwise required longer hospital stays or undergone unnecessary surgery to control the progression of their disease. We expect that the wide distribution of the system we developed will ultimately impact patient care more broadly, including in resource-constrained countries.

Determination of action thresholds for electromagnetic tracking system-guided hypofractionated prostate radiotherapy using volumetric modulated arc therapy.

PURPOSE: Hypofractionated prostate radiotherapy may benefit from both volumetric modulated are therapy (VMAT) due to shortened treatment time and intrafraction real-time monitoring provided by implanted radiofrequency(RF) transponders. The authors investigate dosimetrically driven action thresholds (whether treatment needs to be interrupted and patient repositioned) in VMAT treatment with electromagnetic (EM) tracking. METHODS: VMAT plans for five patients are generated for prescription doses of 32.5 and 42.5 Gy in five fractions. Planning target volume (PTV) encloses the clinical target volume (CTV) with a 3 mm margin at the prostate-rectal interface and 5 mm elsewhere. The VMAT delivery is modeled using 180 equi-spaced static beams. Intrafraction prostate motion is simulated in the plan by displacing the beam isocenter at each beam assuming rigid organ motion according to a previously recorded trajectory of the transponder centroid. The cumulative dose delivered in each fraction is summed over all beams. Two sets of 57 prostate motion trajectories were randomly selected to form a learning and a testing dataset. Dosimetric end points including CTV D95%, rectum wall D1cc, bladder wall D1cc, and urethra Dmax, are analyzed against motion characteristics including the maximum amplitude of the anterior-posterior (AP), superior-inferior (SI), and left-right components. Action thresholds are triggered when intrafraction motion causes any violations of dose constraints to target and organs at risk (OAR), so that treatment is interrupted and patient is repositioned. RESULTS: Intrafraction motion has a little effect on CTV D95%, indicating PTV margins are adequate. Tight posterior and inferior action thresholds around 1 mm need to be set in a patient specific manner to spare organs at risk, especially when the prescription dose is 42.5 Gy. Advantages of setting patient specific action thresholds are to reduce false positive alarms by 25% when prescription dose is low, and increase the sensitivity of detecting dose limits violations by 30% when prescription dose is high, compared to a generic 2 mm action box. The sensitivity and specificity calculated from the testing dataset are consistent to the learning set, which indicates that the patient specific approach is reliable and reproducible within the scope of the prostate database. CONCLUSIONS: This work introduces a formalism for ensuring a VMAT delivery meets the most clinically important dose requirements by using patient specific and dosimetric-driven action thresholds to hold the beam and reposition the patient when necessary. Such methods can provide improved sensitivity and specificity compared to conventional methods, which assume directionally symmetric action thresholds.

Deep learning auto-segmentation and automated treatment planning for trismus risk reduction in head and neck cancer radiotherapy.

BACKGROUND AND PURPOSE: Reducing trismus in radiotherapy for head and neck cancer (HNC) is important. Automated deep learning (DL) segmentation and automated planning was used to introduce new and rarely segmented masticatory structures to study if trismus risk could be decreased. MATERIALS AND METHODS: Auto-segmentation was based on purpose-built DL, and automated planning used our in-house system, ECHO. Treatment plans for ten HNC patients, treated with 2 Gy × 35 fractions, were optimized (ECHO0). Six manually segmented OARs were replaced with DL auto-segmentations and the plans re-optimized (ECHO1). In a third set of plans, mean doses for auto-segmented ipsilateral masseter and medial pterygoid (MIMean, MPIMean), derived from a trismus risk model, were implemented as dose-volume objectives (ECHO2). Clinical dose-volume criteria were compared between the two scenarios (ECHO0 vs. ECHO1; ECHO1 vs. ECHO2; Wilcoxon signed-rank test; significance: p < 0.01). RESULTS: Small systematic differences were observed between the doses to the six auto-segmented OARs and their manual counterparts (median: ECHO1 = 6.2 (range: 0.4, 21) Gy vs. ECHO0 = 6.6 (range: 0.3, 22) Gy; p = 0.007), and the ECHO1 plans provided improved normal tissue sparing across a larger dose-volume range. Only in the ECHO2 plans, all patients fulfilled both MIMean and MPIMean criteria. The population median MIMean and MPIMean were considerably lower than those suggested by the trismus model (ECHO0: MIMean = 13 Gy vs. ≤42 Gy; MPIMean = 29 Gy vs. ≤68 Gy). CONCLUSIONS: Automated treatment planning can efficiently incorporate new structures from DL auto-segmentation, which results in trismus risk sparing without deteriorating treatment plan quality. Auto-planning and deep learning auto-segmentation together provide a powerful platform to further improve treatment planning.

Dose correction strategy for the optimization of volumetric modulated arc therapy.

PURPOSE: Dose calculation during optimization of volumetric modulated arc therapy (VMAT) is necessarily simplified to keep computation time manageably low; however the approximations used in the scatter dose calculation lead to discrepancy with more accurate dose calculation following optimization. The purpose of this study is to develop a dose correction strategy in optimization that can minimize the disagreement. METHODS: VMAT delivery is modeled using a number of static equispaced beams. Dose correction factors (C(ij)) are associated with each beam i and point j inside the region of interest. C(ij) is calculated as the ratio of dose obtained from the full scatter dose calculation over that from the partial scatter dose calculation in optimization. VMAT optimization algorithm is a multiple resolution approach. The dose correction factors are calculated at the beginning of each resolution and applied as multiplicative corrections to the partial scatter dose during optimization. Clinical cases for brain, prostate, paraspinal, and esophagus are utilized to evaluate the method. Treatment plans created with and without the correction scheme are normalized such that the complication rates of organs at risk (OARs) are comparable. The resulting planning target volume (PTV) mean doses are used to compare plan quality. RESULTS: The difference between the dose calculated at the end of optimization and at the end of the final forward dose calculation is reduced from 7% and 5% for the PTV and OAR mean doses without correction to approximately 1% with correction. Applying dose correction during optimization saves planners 2-4 h in average in treatment planning, and has a positive impact on plan quality, evidenced by a noticeably higher PTV mean dose: 2.1%, 2.4%, 0.5%, and 9.3% of the corresponding prescription dose in the brain, esophagus, prostate, and paraspinal cases, respectively. CONCLUSIONS: When dose correction is applied during optimization, dose discrepancies between optimization and full dose calculation are reduced. Integrating dose correction in VMAT optimization allows planners to adjust the optimization constraints more easily and confidently during optimization and has the potential to improve plan quality.