A polar-coordinate-based pencil beam algorithm for VMAT dose computation with high-resolution gantry angle sampling.

PURPOSE: Most commercially available treatment planning systems (TPSs) approximate the continuous delivery of volumetric modulated arc therapy (VMAT) plans with a series of discretized static beams for treatment planning, which can make VMAT dose computation extremely inefficient. In this study, we developed a polar-coordinate-based pencil beam (PB) algorithm for efficient VMAT dose computation with high-resolution gantry angle sampling that can improve the computational efficiency and reduce the dose discrepancy due to the angular under-sampling effect. METHODS AND MATERIALS: 6 MV 1 × 1 m m 2 $1 \times 1{\rm{\;m}}{{\rm{m}}^2}$ pencil beams were simulated on a uniform cylindrical phantom under an EGSnrc Monte Carlo (MC) environment. The MC-generated PB kernels were collected in the polar coordinate system for each bixel on a 40 × 40 c m 2 $40 \times 40{\rm{\;c}}{{\rm{m}}^2}$ fluence map and subsequently fitted via a series of Gaussians. The fluence was calculated using a detectors' eye view with off-axis and MLC transmission factors corrected. Doses of VMAT arc on the phantom were computed by summing the convolution results between the corresponding PB kernels and fluence for each bixel in the polar coordinate system. The convolution was performed using fast Fourier transform to expedite the computing speed. The calculated doses were converted to the Cartesian coordinate system and compared with the reference dose computed by a collapsed cone convolution (CCC) algorithm of the TPS. A heterogeneous phantom was created to study the heterogeneity corrections using the proposed algorithm. Ten VMAT arcs were included to evaluate the algorithm performance. Gamma analysis and computation complexity theory were used to measure the dosimetric accuracy and computational efficiency, respectively. RESULTS: The dosimetric comparisons on the homogeneous phantom between the proposed PB algorithm and the CCC algorithm for 10 VMAT arcs demonstrate that the proposed algorithm can achieve a dosimetric accuracy comparable to that of the CCC algorithm with average gamma passing rates of 96% (2%/2mm) and 98% (3%/3mm). In addition, the proposed algorithm can provide better computational efficiency for VMAT dose computation using a PC equipped with a 4-core processor, compared to the CCC algorithm utilizing a dual 10-core server. Moreover, the computation complexity theory reveals that the proposed algorithm has a great advantage with regard to computational efficiency for VMAT dose computation on homogeneous medium, especially when a fine angular sampling rate is applied. This can support a reduction in dose errors from the angular under-sampling effect by using a finer angular sampling rate, while still preserving a practical computing speed. For dose calculation on the heterogeneous phantom, the proposed algorithm with heterogeneity corrections can still offer a reasonable dosimetric accuracy with comparable computational efficiency to that of the CCC algorithm. CONCLUSIONS: We proposed a novel polar-coordinate-based pencil beam algorithm for VMAT dose computation that enables a better computational efficiency while maintaining clinically acceptable dosimetric accuracy and reducing dose error caused by the angular under-sampling effect. It also provides a flexible VMAT dose computation structure that allows adjustable sampling rates and direct dose computation in regions of interest, which makes the algorithm potentially useful for clinical applications such as independent dose verification for VMAT patient-specific QA.

Tumor phase recognition using cone-beam computed tomography projections and external surrogate information.

PURPOSE: Directly extracting the respiratory phase pattern of the tumor using cone-beam computed tomography (CBCT) projections is challenging due to the poor tumor visibility caused by the obstruction of multiple anatomic structures on the beam's eye view. Predicting tumor phase information using external surrogate also has intrinsic difficulties as the phase patterns between surrogates and tumors are not necessary to be congruent. In this work, we developed an algorithm to accurately recover the primary oscillation components of tumor motion using the combined information from both CBCT projections and external surrogates. METHODS: The algorithm involved two steps. First, a preliminary tumor phase pattern was acquired by applying local principal component analysis (LPCA) on the cropped Amsterdam Shroud (AS) images. In this step, only the cropped image of the tumor region was used to extract the tumor phase pattern in order to minimize the impact of pattern recognition from other anatomic structures. Second, by performing multivariate singular spectrum analysis (MSSA) on the combined information containing both external surrogate signal and the original waveform acquired in the first step, the primary component of the tumor phase oscillation was recovered. For the phantom study, a QUASAR respiratory motion phantom with a removable tumor-simulator insert was employed to acquire CBCT projection images. A comparison between LPCA only and our method was assessed by power spectrum analysis. Also, the motion pattern was simulated under the phase shift or various amplitude conditions to examine the robustness of our method. Finally, anatomic obstruction scenarios were simulated by attaching a heart model, PVC tubes, and RANDO® phantom slabs to the phantom, respectively. Each scenario was tested with five real-patient breathing patterns to mimic real clinical situations. For the patient study, eight patients with various tumor locations were selected. The performance of our method was then evaluated by comparing the reference waveform with the extracted signal for overall phase discrepancy, expiration phase discrepancy, peak, and valley accuracy. RESULTS: In tests of phase shifts and amplitude variations, the overall peak and valley accuracy was -0.009 ± 0.18 sec, and no time delay was found compared to the reference. In anatomical obstruction tests, the extracted signal had 1.6 ± 1.2 % expiration phase discrepancy, -0.12 ± 0.28 sec peak accuracy, and 0.01 ± 0.15 sec valley accuracy. For patient studies, the extracted signal using our method had -1.05 ± 3.0 % overall phase discrepancy, -1.55 ± 1.45% expiration phase discrepancy, 0.04 ± 0.13 sec peak accuracy, and -0.01 ± 0.15 sec valley accuracy, compared to the reference waveforms. CONCLUSIONS: An innovative method capable of accurately recognizing tumor phase information was developed. With the aid of extra information from the external surrogate, an improvement in prediction accuracy, as compared with traditional statistical methods, was obtained. It enables us to employ it as the ground truth for 4D-CBCT reconstruction, gating treatment, and other clinic implementations that require accurate tumor phase information.

Beam flatness modulation for a flattening filter free photon beam utilizing a novel direct leaf trajectory optimization model.

Flattening filter free (FFF) linear accelerators produce a fluence distribution that is forward peaked. Various dosimetric benefits, such as increased dose rate, reduced leakage and out of field dose has led to the growth of FFF technology in the clinic. The literature has suggested the idea of vendors offering dedicated FFF units where the flattening filter (FF) is removed completely and manipulating the beam to deliver conventional flat radiotherapy treatments. This work aims to develop an effective way to deliver modulated flat beam treatments, rather than utilizing a physical FF. This novel optimization model is an extension of the direct leaf trajectory optimization (DLTO) previously developed for volumetric modulated radiation therapy (VMAT) and is capable of accounting for all machine and multileaf collimator (MLC) dynamic delivery constraints, using a combination of linear constraints and a convex objective function. Furthermore, the tongue and groove (T&G) effect was also incorporated directly into our model without introducing nonlinearity to the constraints, nor nonconvexity to the objective function. The overall beam flatness, machine deliverability, and treatment time efficiency were assessed. Regular square fields, including field sizes of 10 × 10 cm2 to 40 × 40 cm2 were analyzed, as well as three clinical fields, and three arbitrary contours with "concave" features. Quantitative flatness was measured for all modulated FFF fields, and the results were comparable or better than their open FF counterparts, with the majority having a quantitative flatness of less than 3.0%. The modulated FFF beams, due to the included efficiency constraint, were able to achieve acceptable delivery time compared to their open FF counterpart. The results indicated that the dose uniformity and flatness for the modulated FFF beams optimized with the DLTO model can successfully match the uniformity and flatness of their conventional FF counterparts, and may even provide further benefit by taking advantage of the unique FFF beam characteristics.

A self-checking treatment couch coordinate calculation system in radiotherapy.

PURPOSE: Traditionally, the treatment couch coordinates (TCCs) for patients undergoing radiotherapy can only be determined at the time of treatment, placing pressure on the treating therapists and leaving several pathways for errors such as wrong-site treatment or wrong treatment table shift from a reference point. The purpose of this work is to propose an accurate, robust, and streamlined system that calculates TCC in advance. METHODS: The proposed system combines the advantages of two different calculation methods that use an indexed immobilization device. The first method uses an array of reference ball bearings (BBs) embedded in the CT scanner's couch-top. To obtain the patient-specific TCC, the spatial offset of the treatment planning isocenter from the reference BB is used. The second method performs a calculation using the one-to-one mapping relationship between the CT scanner's DICOM (Digital Imaging and Communications in Medicine) coordinate system and the TCC system. Both methods use a reference point in the CT coordinate system to correlate a point in the TCC system to perform the coordinate transfer between the two systems. Both methods were used to calculate the TCC and the results were checked against each other, creating an integrated workflow via automated self-checking. The accuracy of the calculation system was retrospectively evaluated with 275 patients, where the actual treatment position determined with cone-beam CT was used as a reference. RESULTS: An efficient workflow transparent to the therapists at both CT simulation and treatment was created. It works with any indexed immobilization device and can be universally applied to all treatment sites. The two methods had comparable accuracy, with 95% of the calculations within 3 mm. The inter-fraction variation was within ± 1.0 cm for 95% of the coordinates across all the treatment sites. CONCLUSIONS: A robust, accurate, and streamlined system was implemented to calculate TCCs in advance. It eases the pressure on the treating therapists, reduces patient setup time, and enhances the patient safety by preventing setup errors.

Angular under-sampling effect on VMAT dose calculation: An analysis and a solution strategy.

PURPOSE: Most VMAT algorithms compute the dose on discretized apertures with small angular separations for practical reasons. However, machines deliver the VMAT dose with a continuously moving MLC and gantry and a continuously changing dose rate. The computed dose can deviate from the delivered dose, especially if no, or loose, MLC movement constraints are applied for the VMAT optimization. The goal of this paper is to establish a simplified mathematical model to analyze the discrepancy between the VMAT plan calculation dose and the delivered dose and to provide a reasonable solution for clinical implementation. METHODS: A simplified metric is first introduced to describe the discrepancy between doses computed with discretized apertures and a continuous delivery model. The delivery fluences were formed separately for six different leaf movement scenarios. The formula was then rewritten in a more general form. The correlation between discretized and continuous fluence is summarized using this general form. The Fourier analysis for the impacts from three separate factors - dose kernel width, aperture width, aperture distance - to the dose discrepancy is also presented in order to provide insight into the dose discrepancy caused by under-sampling in the frequency domain. Finally, a weighting-based interpolation (WBI) algorithm, which can improve the aperture interpolation efficiency, is proposed. The associated evaluation methods and criteria for the proposed algorithm are also given. RESULTS: The comparisons between the WBI algorithm and the equal angular interpolation (EAI) method suggested that the proposed algorithm has a great advantage with regard to aperture number efficiency. To achieve a 90% gamma passing rate using the dose computed with apertures generated with 0.5° EAI, with the initial optimization apertures as the standard for the comparison, the WBI needs only 66% and 54% of the aperture numbers that the EAI method needs for a 2° and a 4° angular separation of the VMAT optimization, respectively. The results also suggested that the weighted dose error index value, Θ, can be used as a stopping criterion for an interpolation algorithm, e.g., WBI or EAI, or as an indicator for sampling level evaluations. The phantom results indicate that the gamma passing rate decreases with increasing depth, from the phantom surface to the iso center, for the plans computed with under-sampled apertures. No obvious variation trends were observed for the plans computed with well-sampled apertures. CONCLUSIONS: The mathematical analysis suggests that the dose discrepancies due to under-sampling are strongly correlated with the aperture width, the distance between apertures, and the width of the dose kernel. The WBI algorithm proves to be an efficient aperture interpolation strategy and is useful for dose computation of VMAT plans.

Image-Guided High-Dose Rate Brachytherapy in Cervix Carcinoma Using Balloon Catheter and Belt Immobilization System.

PURPOSE: The efficacy of image-guided high-dose rate brachytherapy for cervical cancer is limited by the ineffective rectal sparing devices available commercially and the potential applicator movement. We developed a novel device using a balloon catheter and a belt immobilization system, serving for rectal dose reduction and applicator immobilization purposes, respectively. METHODS: The balloon catheter is constructed by gluing a short inflatable tube to a long regular open-end catheter. Contrast agent (10) cm3 is injected into the inflatable end, which is affixed to the tandem and ring applicator, to displace the posterior vaginal wall. The belt immobilization system consists of a specially designed bracket that can hold and fix itself to the applicator, a diaper-like Velcro fastener package used for connecting the patient's pelvis to the bracket, and a buckle that holds the fasteners to stabilize the whole system. The treatment data for 21 patients with cervical cancer using both balloon catheter and belt immobilization system were retrospectively analyzed. Computed tomography and magnetic resonance images, acquired about 30 minutes apart, were registered to evaluate the effectiveness of the immobilization system. RESULTS: In comparison with a virtual rectal blade, the balloon decreased the rectal point dose by 34% ± 4.2% (from 276 ± 57 to 182 ± 38 cGy), corresponding to an extra sparing distance of 7.9 ± 1.1 mm. The maximum sparing distance variation per patient is 1.4 ± 0.6 mm, indicating the high interfractional reproducibility for rectum sparing. With the immobilization system, the mean translational and rotational displacements of the applicator set are <3 mm and <1.5°, respectively, in all directions. CONCLUSIONS: The rectal balloon provides significant dose reduction to the rectum and it may potentially minimize patient discomfort. The immobilization system permits almost no movement of the applicator during treatment. This work has the potential to be promoted as a standardized solution for high-dose rate treatment of cervical cancer.

Parameterization of photon beam dosimetry for a linear accelerator.

PURPOSE: In radiation therapy, accurate data acquisition of photon beam dosimetric quantities is important for (1) beam modeling data input into a treatment planning system (TPS), (2) comparing measured and TPS modeled data, (3) the quality assurance process of a linear accelerator's (Linac) beam characteristics, (4) the establishment of a standard data set for comparison with other data, etcetera. Parameterization of the photon beam dosimetry creates a data set that is portable and easy to implement for different applications such as those previously mentioned. The aim of this study is to develop methods to parameterize photon beam dosimetric quantities, including percentage depth doses (PDDs), profiles, and total scatter output factors (S(cp)). METHODS: S(cp), PDDs, and profiles for different field sizes, depths, and energies were measured for a Linac using a cylindrical 3D water scanning system. All data were smoothed for the analysis and profile data were also centered, symmetrized, and geometrically scaled. The S(cp) data were analyzed using an exponential function. The inverse square factor was removed from the PDD data before modeling and the data were subsequently analyzed using exponential functions. For profile modeling, one halfside of the profile was divided into three regions described by exponential, sigmoid, and Gaussian equations. All of the analytical functions are field size, energy, depth, and, in the case of profiles, scan direction specific. The model's parameters were determined using the minimal amount of measured data necessary. The model's accuracy was evaluated via the calculation of absolute differences between the measured (processed) and calculated data in low gradient regions and distance-to-agreement analysis in high gradient regions. Finally, the results of dosimetric quantities obtained by the fitted models for a different machine were also assessed. RESULTS: All of the differences in the PDDs' buildup and the profiles' penumbra regions were less than 2 and 0.5 mm, respectively. The differences in the low gradient regions were 0.20% ± 0.20% (<1% for all) and 0.50% ± 0.35% (<1% for all) for PDDs and profiles, respectively. For S(cp) data, all of the absolute differences were less than 0.5%. CONCLUSIONS: This novel analytical model with minimum measurement requirements was proved to accurately calculate PDDs, profiles, and S(cp) for different field sizes, depths, and energies.

Impact of scanning parameters and breathing patterns on image quality and accuracy of tumor motion reconstruction in 4D CBCT: a phantom study.

Four-dimensional, cone-beam CT (4D CBCT) substantially reduces respiration-induced motion blurring artifacts in three-dimension (3D) CBCT. However, the image quality of 4D CBCT is significantly degraded which may affect its accuracy in localizing a mobile tumor for high-precision, image-guided radiation therapy (IGRT). The purpose of this study was to investigate the impact of scanning parameters hereinafter collectively referred to as scanning sequence) and breathing patterns on the image quality and the accuracy of computed tumor trajectory for a commercial 4D CBCT system, in preparation for its clinical implementation. We simulated a series of periodic and aperiodic sinusoidal breathing patterns with a respiratory motion phantom. The aperiodic pattern was created by varying the period or amplitude of individual sinusoidal breathing cycles. 4D CBCT scans of the phantom were acquired with a manufacturer-supplied scanning sequence (4D-S-slow) and two in-house modified scanning sequences (4D-M-slow and 4D-M-fast). While 4D-S-slow used small field of view (FOV), partial rotation (200°), and no imaging filter, 4D-M-slow and 4D-M-fast used medium FOV, full rotation, and the F1 filter. The scanning speed was doubled in 4D-M-fast (100°/min gantry rotation). The image quality of the 4D CBCT scans was evaluated using contrast-to-noise ratio (CNR), signal-to-noise ratio (SNR), and motion blurring ratio (MBR). The trajectory of the moving target was reconstructed by registering each phase of the 4D CBCT with a reference CT. The root-mean-squared-error (RMSE) analysis was used to quantify its accuracy. Significant decrease in CNR and SNR from 3D CBCT to 4D CBCT was observed. The 4D-S-slow and 4D-M-fast scans had comparable image quality, while the 4D-M-slow scans had better performance due to doubled projections. Both CNR and SNR decreased slightly as the breathing period increased, while no dependence on the amplitude was observed. The difference of both CNR and SNR between periodic and aperiodic breathing patterns was insignificant (p > 0.48). At end-exhale phases, the motion blurring was negligible for both periodic and aperiodic breathing patterns; at mid-inhale phase, the motion blurring increased as the period, the amplitude or the amount of cycle-to-cycle variation on amplitude increased. Overall, the accuracy of localizing the moving target in 4D CBCT was within 2 mm under all studied cases. No difference in the RMSEs was noticed among the three scanning sequences. The 4D-M-fast scans, free of volume truncation artifacts, exhibited comparable image quality and accuracy in tumor motion reconstruction as the 4D-S-slow scans with reduced imaging dose (0.60 cGy vs. 0.99 cGy) due to the use of faster gantry rotation and the F1 filter, suggesting its suitability for clinical use.

Priori mask guided image reconstruction (p-MGIR) for ultra-low dose cone-beam computed tomography.

Recently, the compressed sensing (CS) based iterative reconstruction method has received attention because of its ability to reconstruct cone beam computed tomography (CBCT) images with good quality using sparsely sampled or noisy projections, thus enabling dose reduction. However, some challenges remain. In particular, there is always a tradeoff between image resolution and noise/streak artifact reduction based on the amount of regularization weighting that is applied uniformly across the CBCT volume. The purpose of this study is to develop a novel low-dose CBCT reconstruction algorithm framework called priori mask guided image reconstruction (p-MGIR) that allows reconstruction of high-quality low-dose CBCT images while preserving the image resolution. In p-MGIR, the unknown CBCT volume was mathematically modeled as a combination of two regions: (1) where anatomical structures are complex, and (2) where intensities are relatively uniform. The priori mask, which is the key concept of the p-MGIR algorithm, was defined as the matrix that distinguishes between the two separate CBCT regions where the resolution needs to be preserved and where streak or noise needs to be suppressed. We then alternately updated each part of image by solving two sub-minimization problems iteratively, where one minimization was focused on preserving the edge information of the first part while the other concentrated on the removal of noise/artifacts from the latter part. To evaluate the performance of the p-MGIR algorithm, a numerical head-and-neck phantom, a Catphan 600 physical phantom, and a clinical head-and-neck cancer case were used for analysis. The results were compared with the standard Feldkamp-Davis-Kress as well as conventional CS-based algorithms. Examination of the p-MGIR algorithm showed that high-quality low-dose CBCT images can be reconstructed without compromising the image resolution. For both phantom and the patient cases, the p-MGIR is able to achieve a clinically-reasonable image with 60 projections. Therefore, a clinically-viable, high-resolution head-and-neck CBCT image can be obtained while cutting the dose by 83%. Moreover, the image quality obtained using p-MGIR is better than the quality obtained using other algorithms. In this work, we propose a novel low-dose CBCT reconstruction algorithm called p-MGIR. It can be potentially used as a CBCT reconstruction algorithm with low dose scan requests.

A method of surface marker location optimization for tumor motion estimation in lung stereotactic body radiation therapy.

PURPOSE: Accurately localizing lung tumor localization is essential for high-precision radiation therapy techniques such as stereotactic body radiation therapy (SBRT). Since direct monitoring of tumor motion is not always achievable due to the limitation of imaging modalities for treatment guidance, placement of fiducial markers on the patient's body surface to act as a surrogate for tumor position prediction is a practical alternative for tracking lung tumor motion during SBRT treatments. In this work, the authors propose an innovative and robust model to solve the multimarker position optimization problem. The model is able to overcome the major drawbacks of the sparse optimization approach (SOA) model. METHODS: The principle-component-analysis (PCA) method was employed as the framework to build the authors' statistical prediction model. The method can be divided into two stages. The first stage is to build the surrogate tumor matrix and calculate its eigenvalues and associated eigenvectors. The second stage is to determine the "best represented" columns of the eigenvector matrix obtained from stage one and subsequently acquire the optimal marker positions as well as numbers. Using 4-dimensional CT (4 DCT) and breath hold CT imaging data, the PCA method was compared to the SOA method with respect to calculation time, average prediction accuracy, prediction stability, noise resistance, marker position consistency, and marker distribution. RESULTS: The PCA and SOA methods which were both tested were on all 11 patients for a total of 130 cases including 4 DCT and breath-hold CT scenarios. The maximum calculation time for the PCA method was less than 1 s with 64 752 surface points, whereas the average calculation time for the SOA method was over 12 min with 400 surface points. Overall, the tumor center position prediction errors were comparable between the two methods, and all were less than 1.5 mm. However, for the extreme scenarios (breath hold), the prediction errors for the PCA method were not only smaller, but were also more stable than for the SOA method. Results obtained by imposing a series of random noises to the surrogates indicated that the PCA method was much more noise resistant than the SOA method. The marker position consistency tests using various combinations of 4 DCT phases to construct the surrogates suggested that the marker position predictions of the PCA method were more consistent than those of the SOA method, in spite of surrogate construction. Marker distribution tests indicated that greater than 80% of the calculated marker positions fell into the high cross correlation and high motion magnitude regions for both of the algorithms. CONCLUSIONS: The PCA model is an accurate, efficient, robust, and practical model for solving the multimarker position optimization problem to predict lung tumor motion during SBRT treatments. Due to its generality, PCA model can also be applied to other imaging guidance system whichever using surface motion as the surrogates.

A further investigation of the centroid-to-centroid method for stereotactic lung radiotherapy: a phantom study.

PURPOSE: Our previous study [B. Lu et al., "A patient alignment solution for lung SBRT setups based on a deformable registration technique," Med. Phys. 39(12), 7379-7389 (2012)] proposed a deformable-registration-based patient setup strategy called the centroid-to-centroid (CTC) method, which can perform an accurate alignment of internal-target-volume (ITV) centroids between averaged four-dimensional computed tomography and cone-beam computed tomography (CBCT) images. Scenarios with variations between CBCT and simulation CT caused by irregular breathing and/or tumor change were not specifically considered in the patient study [B. Lu et al., "A patient alignment solution for lung SBRT setups based on a deformable registration technique," Med. Phys. 39(12), 7379-7389 (2012)] due to the lack of both a sufficiently large patient data sample and a method of tumor tracking. The aim of this study is to thoroughly investigate and compare the impacts of breathing pattern and tumor change on both the CTC and the translation-only (T-only) gray-value mode strategies by employing a four-dimensional (4D) lung phantom. METHODS: A sophisticated anthropomorphic 4D phantom (CIRS Dynamic Thorax Phantom model 008) was employed to simulate all desired respiratory variations. The variation scenarios were classified into four groups: inspiration to expiration ratio (IE ratio) change, tumor trajectory change, tumor position change, tumor size change, and the combination of these changes. For each category the authors designed several scenarios to demonstrate the effects of different levels of breathing variation on both of the T-only and the CTC methods. Each scenario utilized 4DCT and CBCT scans. The ITV centroid alignment discrepancies for CTC and T-only were evaluated. The dose-volume-histograms (DVHs) of ITVs for two extreme cases were analyzed. RESULTS: Except for some extreme cases in the combined group, the accuracy of the CTC registration was about 2 mm for all cases for both the single and the combined scenarios. The performance of the CTC method was insensitive to region-of-registration (ROR) size selections, as suggested by the comparable accuracy between 1 and 2 cm expansions of the ROR selections for the method. The T-only method was suitable for some single scenarios, such as trajectory variation, position variation, and size variation. However, for combined scenarios and/or a large variation in the IE ratio, the T-only method failed to produce reasonable registration results (within 3 mm). The discrepancy was close to, or even greater than, 1 cm. In addition, unlike the CTC method, the T-only method was sensitive to the ROR size selection. The DVH analysis suggested that a large ITV to PTV margin should be considered if a breathing pattern variation is observed. CONCLUSIONS: The phantom study demonstrated that the CTC method was reliable for scenarios in which breathing pattern variation was involved. The T-only gray value method worked for some scenarios, but not for scenarios that involved an IE ratio variation. For scenarios involving position variation, the T-only method worked only with a careful selection of the ROR, whereas the CTC method was independent of ROR size as long as the ITVs were included in the ROR. One indication of the dose consequence analysis was that a large ITV to PTV margin should be considered if a breathing pattern variation is observed.

A study of IMRT planning parameters on planning efficiency, delivery efficiency, and plan quality.

PURPOSE: To improve planning and delivery efficiency of head and neck IMRT without compromising planning quality through the evaluation of inverse planning parameters. METHODS: Eleven head and neck patients with pre-existing IMRT treatment plans were selected for this retrospective study. The Pinnacle treatment planning system (TPS) was used to compute new treatment plans for each patient by varying the individual or the combined parameters of dose∕fluence grid resolution, minimum MU per segment, and minimum segment area. Forty-five plans per patient were generated with the following variations: 4 dose∕fluence grid resolution plans, 12 minimum segment area plans, 9 minimum MU plans, and 20 combined minimum segment area∕minimum MU plans. Each plan was evaluated and compared to others based on dose volume histograms (DVHs) (i.e., plan quality), planning time, and delivery time. To evaluate delivery efficiency, a model was developed that estimated the delivery time of a treatment plan, and validated through measurements on an Elekta Synergy linear accelerator. RESULTS: The uncertainty (i.e., variation) of the dose-volume index due to dose calculation grid variation was as high as 8.2% (5.5 Gy in absolute dose) for planning target volumes (PTVs) and 13.3% (2.1 Gy in absolute dose) for planning at risk volumes (PRVs). Comparison results of dose distributions indicated that smaller volumes were more susceptible to uncertainties. The grid resolution of a 4 mm dose grid with a 2 mm fluence grid was recommended, since it can reduce the final dose calculation time by 63% compared to the accepted standard (2 mm dose grid with a 2 mm fluence grid resolution) while maintaining a similar level of dose-volume index variation. Threshold values that maintained adequate plan quality (DVH results of the PTVs and PRVs remained satisfied for their dose objectives) were 5 cm(2) for minimum segment area and 5 MU for minimum MU. As the minimum MU parameter was increased, the number of segments and delivery time were decreased. Increasing the minimum segment area parameter decreased the plan MU, but had less of an effect on the number of segments and delivery time. Our delivery time model predicted delivery time to within 1.8%. CONCLUSIONS: Increasing the dose grid while maintaining a small fluence grid allows for improved planning efficiency without compromising plan quality. Delivery efficiency can be improved by increasing the minimum MU, but not the minimum segment area. However, increasing the respective minimum MU and∕or the minimum segment area to any value greater than 5 MU and 5 cm(2) is not recommended because it degrades plan quality.

A patient alignment solution for lung SBRT setups based on a deformable registration technique.

PURPOSE: In this work, the authors propose a novel registration strategy for translation-only correction scenarios of lung stereotactic body radiation therapy setups, which can achieve optimal dose coverage for tumors as well as preserve the consistency of registrations with minimal human interference. METHODS: The proposed solution (centroid-to-centroidor CTC solution) uses the average four-dimensional CT (A4DCT) as the reference CT. The cone-beam CT (CBCT) is deformed to acquire a new centroid for the internal target volume (ITV) on the CBCT. The registration is then accomplished by simply aligning the centroids of the ITVs between the A4DCT and the CBCT. Sixty-seven cases using 64 patients (each case is associated with separate isocenters) have been investigated with the CTC method and compared with the conventional gray-value (G) mode and bone (B) mode registration methods. Dosimetric effects among the tree methods were demonstrated by 18 selected cases. The uncertainty of the CTC method has also been studied. RESULTS: The registration results demonstrate the superiority of the CTC method over the other two methods. The differences in the D99 and D95 ITV dose coverage between the CTC method and the original plan is small (within 5%) for all of the selected cases except for one for which the tumor presented significant growth during the period between the CT scan and the treatment. Meanwhile, the dose coverage differences between the original plan and the registration results using either the B or G method are significant, as tumor positions varied dramatically, relative to the rib cage, from their positions on the original CT. The largest differences between the D99 and D95 dose coverage of the ITV using the B or G method versus the original plan are as high as 50%. The D20 differences between any of the methods versus the original plan are all less than 2%. CONCLUSIONS: The CTC method can generate optimal dose coverage to tumors with much better consistency compared with either the G or B method, and it is especially useful when the tumor position varies greatly from its position on the original CT, relative to the rib cage.

An approach for online evaluations of dose consequences caused by small rotational setup errors in intracranial stereotactic radiation therapy.

PURPOSE: The purpose of this work is to investigate the impact of small rotational errors on the magnitudes and distributions of spatial dose variations for intracranial stereotactic radiotherapy (SRT) treatment setups, and to assess the feasibility of using the original dose map overlaid with rotated contours (ODMORC) method as a fast, online evaluation tool to estimate dose changes (using DVHs) to clinical target volumes (CTVs) and organs-at-risks (OARs) caused by small rotational setup errors. METHODS: Fifteen intracranial SRT cases treated with either three-dimensional conformal radiation therapy (3DCRT) or intensity-modulated radiation therapy (IMRT) techniques were chosen for the study. Selected cases have a variety of anatomical dimensions and pathologies. Angles of ±3° and ±5° in all directions were selected to simulate the rotational errors. Dose variations in different regions of the brain, CTVs, and OARs were evaluated to illustrate the various spatial effects of dose differences before and after rotations. DVHs accounting for rotations that were recomputed by the treatment planning system (TPS) and those generated by the ODMORC method were compared. A framework of a fast algorithm for multicontour rotation implemented by ODMORC is introduced as well. RESULTS: The average values of relative dose variations between original dose and recomputed dose accounting for rotations were greater than 4.0% and 10.0% in absolute mean and in standard deviation, respectively, at the skull and adjacent regions for all cases. They were less than 1.0% and 2.5% in absolute mean and in standard deviation, respectively, for dose points 3 mm away from the skull. The results indicated that spatial dose to any part of the brain organs or tumors separated from the skull or head surface would be relatively stable before and after rotations. Statistical data of CTVs and OARs indicate the lens and cochleas have the large dose variations before and after rotations, whereas the remaining ROIs have insignificant dose differences. DVH comparisons suggest that the ODMORC method is able to estimate the DVH of CTVs fairly accurately (within 1.5% of relative dose differences for evaluation volumes). The results also show that most of the OARs including the brain stem, spinal cord, chiasm, hippocampuses, optic nerves, and retinas, which were relatively distal from the skull and surface, had good agreement (within 2.0% of relative dose differences for 0.1 cc of the volumes ) between the ODMORC method and the recomputation, whereas OARs more proximate to the bone-tissue interface or surface, such as the lenses and cochlea, had larger dose variations (greater than 5.0%) for some cases due to the incapability of the ODMORC to account for scatter contribution variations proximate to interfaces and intrinsic dose calculation uncertainties for ROIs with small volumes. CONCLUSIONS: The ODMORC method can be implemented as an online evaluation system for rotation-induced dose changes of CTVs and most OARs and for other related dose consequence analyses.

Feasibility study of performing IGRT system daily QA using a commercial QA device.

The purpose of this study was to investigate the feasibility of using a single QA device for comprehensive, efficient daily QA of a linear accelerator (Linac) and three image-guided stereotactic positioning systems (IGSPSs). The Sun Nuclear Daily QA 3 (DQA3) device was used to perform daily dosimetry and mechanical accuracy tests for an Elekta Linac, as well as daily image geometric and isocenter coincidence accuracy tests for three IGSPSs: the AlignRT surface imaging system; the frameless SonArray optical tracking System (FSA) and the Elekta kV CBCT. The DQA3 can also be used for couch positioning, repositioning, and rotational tests during the monthly QA. Based on phantom imaging, the Linac coordinate system determined using AlignRT was within 0.7 mm/0.6° of that of the CBCT system. The difference is attributable to the different calibration methods that are utilized for these two systems. The laser alignment was within 0.5 mm of the isocenter location determined with the three IGSPSs. The ODI constancy was ± 0.5 mm. For gantry and table angles of 0°, the mean isocenter displacement vectors determined using the three systems were within 0.7 mm and 0.6° of one another. Isocenter rotational offsets measured with the systems were all ≤ 0.5°. For photon and electron beams tested over a period of eight months, the output was verified to remain within 2%, energy variations were within 2%, and the symmetry and flatness were within 1%. The field size and light-radiation coincidence were within 1mm ± 1 mm. For dosimetry reproducibility, the standard deviation was within 0.2% for all tests and all energies, except for photon energy variation which was 0.6%. The total measurement time for all tasks took less than 15 minutes per QA session compared to 40 minutes with our previous procedure, which utilized an individual QA device for each IGSPS. The DQA3 can be used for accurate and efficient Linac and IGSPS daily QA. It shortens QA device setup time, eliminates errors introduced by changing phantoms to perform different tests, and streamlines the task of performing dosimetric checks.

Characterization of a multi-axis ion chamber array.

PURPOSE: The aim of this work was to characterize a multi-axis ion chamber array (IC PROFILER; Sun Nuclear Corporation, Melbourne, FL, USA) that has the potential to simplify the acquisition of LINAC beam data. METHODS: The IC PROFILER (or panel) measurement response was characterized with respect to radiation beam properties, including dose, dose per pulse, pulse rate frequency (PRF), and energy. Panel properties were also studied, including detector-calibration stability, power-on time, backscatter dependence, and the panel's agreement with water tank measurements [profiles, fractional depth dose (FDD), and output factors]. RESULTS: The panel's relative deviation was typically within (+/-) 1% of an independent (or nominal) response for all properties that were tested. Notable results were (a) a detectable relative field shape change of approximately 1% with linear accelerator PRF changes; (b) a large range in backscatter thickness had a minimal effect on the measured dose distribution (typically less than 1%); (c) the error spread in profile comparison between the panel and scanning water tank (Blue Phantom, CC13; IBA Schwarzenbruck, DE) was approximately (+/-) 0.75%. CONCLUSIONS: The ability of the panel to accurately reproduce water tank profiles, FDDs, and output factors is an indication of its abilities as a dosimetry system. The benefits of using the panel versus a scanning water tank are less setup time and less error susceptibility. The same measurements (including device setup and breakdown) for both systems took 180 min with the water tank versus 30 min with the panel. The time-savings increase as the measurement load is increased.

Characterization of a real-time surface image-guided stereotactic positioning system.

PURPOSE: The AlignRT3C system is an image-guided stereotactic positioning system (IGSPS) that provides real-time target localization. This study involves the first use of this system with three camera pods. The authors have evaluated its localization accuracy and tracking ability using a cone-beam computed tomography (CBCT) system and an optical tracking system in a clinical setting. METHODS: A modified Rando head-and-neck phantom and five patients receiving intracranial stereotactic radiotherapy (SRT) were used to evaluate the calibration, registration, and position-tracking accuracies of the AlignRT3C system and to study surface reconstruction uncertainties, including the effects due to interfractional and intrafractional motion, skin tone, room light level, camera temperature, and image registration region of interest selection. System accuracy was validated through comparison with the Elekta kV CBCT system (XVI) and the Varian frameless SonArray (FSA) optical tracking system. Surface-image data sets were acquired with the AlignRT3C daily for the evaluation of pretreatment and interfractional and intrafractional motion for each patient. Results for two different reference image sets, planning CT surface contours (CTS) and previously recorded AlignRT3C optical surface images (ARTS), are reported. RESULTS: The system origin displacements for the AlignRT3C and XVI systems agreed to within 1.3 mm and 0.7 degrees. Similar results were seen for AlignRT3C vs FSA. For the phantom displacements having couch angles of 0 degrees, those that utilized ART_S references resulted in a mean difference of 0.9 mm/0.4 degrees with respect to XVI and 0.3 mm/0.2 degrees with respect to FSA. For phantom displacements of more than +/- 10 mm and +/- 3 degrees, the maximum discrepancies between AlignRT and the XVI and FSA systems were 3.0 and 0.4 mm, respectively. For couch angles up to +/- 90 degrees, the mean (max.) difference between the AlignRT3C and FSA was 1.2 (2.3) mm/0.7 degrees (1.2 degrees). For all tests, the mean registration errors obtained using the CT_S references were approximately 1.3 mm/1.0 degrees larger than those obtained using the ART_S references. For the patient study, the mean differences in the pretreatment displacements were 0.3 mm/0.2 degrees between the AlignRT3C and XVI systems and 1.3 mm/1 degrees between the FSA and XVI systems. For noncoplanar treatments, interfractional motion displacements obtained using the ART_S and CT_S references resulted in 90th percentile differences within 2.1 mm/0.8 degrees and 3.3 mm/0.3 degrees, respectively, compared to the FSA system. Intrafractional displacements that were tracked for a maximum of 14 min were within 1 mm/1 degrees of those obtained with the FSA system. Uncertainties introduced by the bite-tray were as high as 3 mm/2 degrees for one patient. The combination of gantry, aSi detector panel, and x-ray tube blockage effects during the CBCT acquisition resulted in a registration error of approximately 3 mm. No skin-tone or surface deformation effects were seen with the limited patient sample. CONCLUSIONS: AlignRT3C can be used as a nonionizing IGSPS with accuracy comparable to current image/marker-based systems. IGSPS and CBCT can be combined for high-precision positioning without the need for patient-attached localization devices.

Wide field array calibration dependence on the stability of measured dose distributions.

PURPOSE: The aim of this work was to simulate the effect of dose distribution changes on detector array calibrations and to explore compensatory methods that are used during calibration measurements. METHODS: The array calibration technique that was investigated is known as wide field (WF) calibration. Using this method, a linear array [y-axis (65 detectors) of the IC PROFILER (Sun Nuclear Corporation, Melbourne, FL)] is calibrated with three measurements (alpha, theta, and lamda); each measurement uses the same radiation field, which is larger than the array. For measurement configuration theta, the array is rotated by 180 degrees from its position in a; for lamda, the array is shifted by one detector from its position in theta. The relative detector sensitivities are then determined through ratios of detector readings at the same field locations (using theta and lamda). This method results in error propagation that is proportional to the number of detectors in the array. During the procedure, the calibration protocol operates under three postulates, which state that (a) the beam shape does not change between measurements; (b) the relative sensitivities of the detectors do not change; and (c) the scatter to the array does not change as the array is moved. The WF calibration's sensitivity to a postulate (a) violation was quantified by applying a sine shaped perturbation (of up to 0.1%) to a, theta, or lamda, and then determining the change relative to a baseline calibration. Postulate (a) violations were minimized by using a continuous beam and mechanized array movement during theta and lamda. A continuously on beam demonstrated more stable beam symmetry as compared to cycling the beam on and off between measurements. Additional side-scatter was also used to satisfy postulate (c). RESULTS: Simulated symmetry perturbations of 0.1% to theta or lamda resulted in calibration errors of up to 2%; alpha was relatively immune to perturbation (<0.1% error). Wide field calibration error on a linear accelerator with similar symmetry variations was +/- 1.6%. Using a continuous beam during theta and lamda with additional side-scatter reduced the calibration error from +/- 1.6% to +/- 0.48%. CONCLUSIONS: This work increased the reproducibility of WF calibrations by limiting the effect of measurement perturbations primarily due to linear accelerator symmetry variations. The same technique would work for any array using WF calibration.

Quality assessment of frameless fractionated stereotactic radiotherapy using cone beam computed tomography.

PURPOSE: A quality assessment of intracranial stereotactic radiotherapy was performed using cone beam computed tomography (CBCT). Setup errors were analyzed for two groups of patients: (1) those who were positioned using a frameless SonArray (FSA) system and immobilized with a bite plate and thermoplastic (TP) mask (the bFSA group); and (2) those who were positioned by room laser and immobilized using a TP mask (the mLAS group). METHODS AND MATERIALS: A quality assurance phantom was used to study the system differences between FSA and CBCT. The quality assessment was performed using an Elekta Synergy imager (XVI) (Elekta Oncology Systems, Norcross, GA) and an On-Board Imager (OBI) (Varian Medical Systems, Palo Alto, CA) for 25 patients. For the first three fractions, and weekly thereafter, the FSA system was used for patient positioning, after which CBCT was performed to obtain setup errors. RESULTS: (1) Phantom tests: The mean differences in the isocenter displacements for the two systems was 1.2 ± 0.7 mm. No significant variances were seen between the XVI and OBI units (p~0.208). (2)Patient tests: The mean of the displacements between FSA and CBCT were independent of the CBCT system used; mean setup errors for the bFSA group were smaller (1.2 mm) than those of the mLAS group (3.2 mm) (p < 0.005). For the mLAS patients, the 90th percentile and the maximum rotational displacements were 3° and 5°, respectively. A 4-mm drift in setup accuracy occurred over the treatment course for 1 bFSA patient. CONCLUSIONS: System differences of less than 1 mm between CBCT and FSA were seen. Error regression was observed for the bFSA patients, using CBCT (up to 4 mm) during the treatment course. For the mLAS group, daily CBCT imaging was needed to obtain acceptable setup accuracies.

Assessment of the setup dependence of detector response functions for mega-voltage linear accelerators.

PURPOSE: Accurate modeling of beam profiles is important for precise treatment planning dosimetry. Calculated beam profiles need to precisely replicate profiles measured during machine commissioning. Finite detector size introduces perturbations into the measured profiles, which, in turn, impact the resulting modeled profiles. The authors investigate a method for extracting the unperturbed beam profiles from those measured during linear accelerator commissioning. METHODS: In-plane and cross-plane data were collected for an Elekta Synergy linac at 6 MV using ionization chambers of volume 0.01, 0.04, 0.13, and 0.65 cm3 and a diode of surface area 0.64 mm2. The detectors were orientated with the stem perpendicular to the beam and pointing away from the gantry. Profiles were measured for a 10 x 10 cm2 field at depths ranging from 0.8 to 25.0 cm and SSDs from 90 to 110 cm. Shaping parameters of a Gaussian response function were obtained relative to the Edge detector. The Gaussian function was deconvolved from the measured ionization chamber data. The Edge detector profile was taken as an approximation to the true profile, to which deconvolved data were compared. Data were also collected with CC13 and Edge detectors for additional fields and energies on an Elekta Synergy, Varian Trilogy, and Siemens Oncor linear accelerator and response functions obtained. Response functions were compared as a function of depth, SSD, and detector scan direction. Variations in the shaping parameter were introduced and the effect on the resulting deconvolution profiles assessed. RESULTS: Up to 10% setup dependence in the Gaussian shaping parameter occurred, for each detector for a particular plane. This translated to less than a +/- 0.7 mm variation in the 80%-20% penumbral width. For large volume ionization chambers such as the FC65 Farmer type, where the cavity length to diameter ratio is far from 1, the scan direction produced up to a 40% difference in the shaping parameter between in-plane and cross-plane measurements. This is primarily due to the directional difference in penumbral width measured by the FC65 chamber, which can more than double in profiles obtained with the detector stem parallel compared to perpendicular to the scan direction. For the more symmetric CC13 chamber the variation was only 3% between in-plane and cross-plane measurements. CONCLUSIONS: The authors have shown that the detector response varies with detector type, depth, SSD, and detector scan direction. In-plane vs. cross-plane scanning can require calculation of a direction dependent response function. The effect of a 10% overall variation in the response function, for an ionization chamber, translates to a small deviation in the penumbra from that of the Edge detector measured profile when deconvolved. Due to the uncertainties introduced by deconvolution the Edge detector would be preferable in obtaining an approximation of the true profile, particularly for field sizes where the energy dependence of the diode can be neglected. However, an averaged response function could be utilized to provide a good approximation of the true profile for large ionization chambers and for larger fields for which diode detectors are not recommended.