N Wen

Dose delivered from Varian's CBCT to patients receiving IMRT for prostate cancer

With the increased use of cone beam CT (CBCT) for daily patient setup, the accumulated dose from CBCT may be significantly higher than that from simulation CT or portal imaging. The objective of this work is to measure the dose from daily pelvic scans with fixed technical settings and collimations. CBCT scans were acquired in half-fan mode using a half bowtie and x-rays were delivered in pulsed-fluoro mode. The skin doses for seven prostate patients were measured on an IRB-approved protocol. TLD capsules were placed on the patients skin at the central axis of three beams: AP, left lateral (Lt Lat) and right lateral (Rt Lat). To avoid the ring artefacts centred in the prostate, the treatment couch was dropped 3 cm from the patients tattoo (central axis). The measured AP skin doses ranged 3-6 cGy for 20-33 cm separation. The larger the patient size the less the AP skin dose. Lateral doses did not change much with patient size. The Lt Lat dose was approximately 4.0 cGy, which was approximately 40% higher than the Rt Lat dose of approximately 2.6 cGy. To verify this dose asymmetry, surface doses on an IMRT QA phantom (oval shaped, 30 cm x 20 cm) were measured at the same three sites using TLD capsules with 3 cm table-drop. The dose asymmetry was due to: (1) kV source rotation which always starts from the patients Lt Lat and ends at Lt Lat. Gantry rotation gets much slower near the end of rotation but dose rate stays constant and (2) 370 degrees scan rotation (10 degrees scan overlap on the Lt Lat side). In vivo doses were measured inside a Rando pelvic heterogeneous phantom using TLDs. The left hip (femoral head and neck) received the highest doses of approximately 10-11 cGy while the right hip received approximately 6-7 cGy. The surface and in vivo doses were also measured for phantoms at the central-axis setup. The difference was less than approximately 12% to the table-drop setup.

Dosimetric verification and clinical evaluation of a new commercially available Monte Carlo-based dose algorithm for application in stereotactic body radiation therapy (SBRT) treatment planning.

Modern cancer treatment techniques, such as intensity-modulated radiation therapy (IMRT) and stereotactic body radiation therapy (SBRT), have greatly increased the demand for more accurate treatment planning (structure definition, dose calculation, etc) and dose delivery. The ability to use fast and accurate Monte Carlo (MC)-based dose calculations within a commercial treatment planning system (TPS) in the clinical setting is now becoming more of a reality. This study describes the dosimetric verification and initial clinical evaluation of a new commercial MC-based photon beam dose calculation algorithm, within the iPlan v.4.1 TPS (BrainLAB AG, Feldkirchen, Germany). Experimental verification of the MC photon beam model was performed with film and ionization chambers in water phantoms and in heterogeneous solid-water slabs containing bone and lung-equivalent materials for a 6 MV photon beam from a Novalis (BrainLAB) linear accelerator (linac) with a micro-multileaf collimator (m(3) MLC). The agreement between calculated and measured dose distributions in the water phantom verification tests was, on average, within 2%/1 mm (high dose/high gradient) and was within +/-4%/2 mm in the heterogeneous slab geometries. Example treatment plans in the lung show significant differences between the MC and one-dimensional pencil beam (PB) algorithms within iPlan, especially for small lesions in the lung, where electronic disequilibrium effects are emphasized. Other user-specific features in the iPlan system, such as options to select dose to water or dose to medium, and the mean variance level, have been investigated. Timing results for typical lung treatment plans show the total computation time (including that for processing and I/O) to be less than 10 min for 1-2% mean variance (running on a single PC with 8 Intel Xeon X5355 CPUs, 2.66 GHz). Overall, the iPlan MC algorithm is demonstrated to be an accurate and efficient dose algorithm, incorporating robust tools for MC-based SBRT treatment planning in the routine clinical setting.

Combining scatter reduction and correction to improve image quality in cone‐beam computed tomography (CBCT)

PURPOSE The authors propose a combined scatter reduction and correction method to improve image quality in cone-beam computed tomography (CBCT). Although using a beam-block approach similar to previous techniques to measure the scatter, this method differs in that the authors utilize partially blocked projection data obtained during scatter measurement for CBCT image reconstruction. This study aims to evaluate the feasibility of the proposed approach. METHODS A 1D grid, composed of lead septa, was placed between the radiation source and the imaging object for scatter measurement. Image data were collected from the grid interspace regions while the scatter distribution was measured in the blocked regions under the grid. Scatter correction was performed by subtracting the measured scatter from the imaging data. Image information in the penumbral regions of the grid was derived. Three imaging modes were developed to reconstruct full CBCT images from partial projection data. The single-rotation half-fan mode uses interpolation to fill the missing data. The dual-rotation half-fan mode uses two rotations, with the grid offset by half a grid cycle, to acquire two complementary sets of projections, which are then merged to form complete projections for reconstruction. The single-rotation full-fan mode was designed for imaging a small object or a region of interest. Full-fan projection images were acquired over a 360 degrees scan angle with the grid shifting a distance during the scan. An enlarged Catphan phantom was used to evaluate potential improvement in image quality with the proposed technique. An anthropomorphic pelvis phantom was used to validate the feasibility of reconstructing a complete set of CBCT images from the partially blocked projections using three imaging modes. Rigid-body image registration was performed between the CBCT images from the single-rotation half-fan mode and the simulation CT and the results were compared to that for the CBCT images from dual-rotation mode and conventional CBCT images. RESULTS The proposed technique reduced the streak artifact index from 58% to 1% in comparison with the conventional CBCT. It also improved CT number linearity from 0.880 to 0.998 and the contrast-to-noise ratio (CNR) from 4.29 to 6.42. Complete sets of CBCT images with overall improved image quality were achieved for all three image modes. The longitudinal resolution was slightly compromised for the single-rotation half-fan mode. High resolution was retained for the dual-rotation half-fan and single-rotation full-fan modes in the longitudinal direction. The registration error for the CBCT images from the single-rotation half-fan mode was 0.8 +/- 0.3 mm in the longitudinal direction and negligible in the other directions. CONCLUSIONS The proposed method provides combined scatter correction and direct scatter reduction. Scatter correction may eliminate scatter artifacts, while direct scatter reduction may improve the CNR to compensate the CNR degradation due to scatter correction. Complete sets of CBCT images are reconstructed in all three imaging modes. The single-rotation mode can be used for rigid-body patient alignment despite degradation in longitudinal resolution. The dual-rotation mode may be used to improve CBCT image quality for soft tissue delineation in adaptive radiation therapy.

Evaluation of multiple image-based modalities for image-guided radiation therapy (IGRT) of prostate carcinoma: A prospective study

PURPOSE Setup errors and prostate intrafraction motion are main sources of localization uncertainty in prostate cancer radiation therapy. This study evaluates four different imaging modalities 3D ultrasound (US), kV planar images, cone-beam computed tomography (CBCT), and implanted electromagnetic transponders (Calypso/Varian) to assess inter- and intrafraction localization errors during intensity-modulated radiation therapy based treatment of prostate cancer. METHODS Twenty-seven prostate cancer patients were enrolled in a prospective IRB-approved study and treated to a total dose of 75.6 Gy (1.8 Gy/fraction). Overall, 1100 fractions were evaluated. For each fraction, treatment targets were localized using US, kV planar images, and CBCT in a sequence defined to determine setup offsets relative to the patient skin tattoos, intermodality differences, and residual errors for each patient and patient cohort. Planning margins, following van Herks formalism, were estimated based on error distributions. Calypso-based localization was not available for the first eight patients, therefore centroid positions of implanted gold-seed markers imaged prior to and immediately following treatment were used as a motion surrogate during treatment. For the remaining 19 patients, Calypso transponders were used to assess prostate intrafraction motion. RESULTS The means (μ), and standard deviations (SD) of the systematic (Σ) and random errors (σ) of interfraction prostate shifts (relative to initial skin tattoo positioning), as evaluated using CBCT, kV, and US, averaged over all patients and fractions, were: [μ CBCT = (-1.2, 0.2, 1.1) mm, Σ CBCT = (3.0, 1.4, 2.4) mm, σ CBCT = (3.2, 2.2, 2.5) mm], [μkV = (-2.9, -0.4, 0.5) mm, Σ kV = (3.4, 3.1, 2.6) mm, σ kV = (2.9, 2.0, 2.4) mm], and [μ US = (-3.6, -1.4, 0.0) mm, Σ US = (3.3, 3.5, 2.8) mm, σ US = (4.1, 3.8, 3.6) mm], in the anterior-posterior (A/P), superior-inferior (S/I), and the left-right (L/R) directions, respectively. In the treatment protocol, adjustment of couch was guided by US images. Residual setup errors as assessed by kV images were found to be: μ residual = (-0.4, 0.2, 0.2) mm, Σ residual = (1.0, 1.0,0.7) mm, and σ residual = (2.5, 2.3, 1.8) mm. Intrafraction prostate motion, evaluated using electromagnetic transponders, was: μ intrafxn = (0.0, 0.0, 0.0) mm, Σ intrafxn = (1.3, 1.5, 0.6) mm, and σ intrafxn = (2.6, 2.4, 1.4) mm. Shifts between pre- and post-treatment kV images were: μ kV(post-pre) = (-0.3, 0.8, -0.2), Σ kV(post-pre) = (2.4, 2.7, 2.1) mm, and σ kV(post-pre) = (2.7, 3.2, 3.1) mm. Relative to skin tattoos, planning margins for setup error were within 10-11 mm for all image-based modalities. The use of image guidance was shown to reduce these margins to less than 5 mm. Margins to compensate for both residual setup (interfraction) errors as well as intrafraction motion were 6.6, 6.8, and 3.9 mm in the A/P, S/I, and L/R directions, respectively. CONCLUSIONS Analysis of interfraction setup errors, performed with US, CBCT, planar kV images, and electromagnetic transponders, from a large dataset revealed intermodality shifts were comparable (within 3-4 mm). Interfraction planning margins, relative to setup based on skin marks, were generally within the 10 mm prostate-to-planning target volume margin used in our clinic. With image guidance, interfraction residual planning margins were reduced to approximately less than 4 mm. These findings are potentially important for dose escalation studies using smaller margins to better protect normal tissues.

Implementation of a Novel Algorithm For Generating Synthetic CT Images From Magnetic Resonance Imaging Data Sets for Prostate Cancer Radiation Therapy

PURPOSE To describe and evaluate a method for generating synthetic computed tomography (synCT) images from magnetic resonance simulation (MR-SIM) data for accurate digitally reconstructed radiograph (DRR) generation and dose calculations in prostate cancer radiation therapy. METHODS AND MATERIALS A retrospective evaluation was performed in 9 prostate cancer patients who had undergone MR-SIM in addition to CT simulation (CT-SIM). MR-SIM data were used to generate synCT images by using a novel, voxel-based weighted summation approach. A subset of patients was used for weight optimization, and the number of patients to use during optimization was determined. Hounsfield unit (HU) differences between CT-SIM and synCT images were analyzed via mean absolute error (MAE). Original, CT-based treatment plans were mapped onto synCTs. DRRs were generated, and agreement between CT and synCT-generated DRRs was evaluated via Dice similarity coefficient (DSC). Dose was recalculated, and dose-volume metrics and gamma analysis were used to evaluate resulting treatment plans. RESULTS Full field-of-view synCT MAE across all patients was 74.3 ± 10.9 HU with differences from CTs of 2.0 ± 8.1 HU and 11.9 ± 46.7 HU for soft tissue structures (prostate, bladder, and rectum) and femoral bones, respectively. Calculated DSCs for anterior-posterior and lateral DRRs were 0.90 ± 0.04 and 0.92 ± 0.05, respectively. Differences in D99%, mean dose, and maximum dose to the clinical target volume from CT-SIM dose calculations were 0.75% ± 0.35%, 0.63% ± 0.34%, and 0.54% ± 0.33%, respectively, for synCT-generated plans. Gamma analysis (2%/2 mm dose difference/distance to agreement) revealed pass rates of 99.9% ± 0.1% (range, 99.7%-100%). CONCLUSION Generated synCTs enabled accurate DRR generation and dose computation for prostate MR-only simulation. Dose recalculated on synCTs agreed well with original planning distributions. Further validation using a larger patient cohort is warranted.

Using patient-specific phantoms to evaluate deformable image registration algorithms for adaptive radiation therapy

The quality of adaptive treatment planning depends on the accuracy of its underlying deformable image registration (DIR). The purpose of this study is to evaluate the performance of two DIR algorithms, B‐spline‐based deformable multipass (DMP) and deformable demons (Demons), implemented in a commercial software package. Evaluations were conducted using both computational and physical deformable phantoms. Based on a finite element method (FEM), a total of 11 computational models were developed from a set of CT images acquired from four lung and one prostate cancer patients. FEM generated displacement vector fields (DVF) were used to construct the lung and prostate image phantoms. Based on a fast‐Fourier transform technique, image noise power spectrum was incorporated into the prostate image phantoms to create simulated CBCT images. The FEM‐DVF served as a gold standard for verification of the two registration algorithms performed on these phantoms. The registration algorithms were also evaluated at the homologous points quantified in the CT images of a physical lung phantom. The results indicated that the mean errors of the DMP algorithm were in the range of 1.0~3.1mm for the computational phantoms and 1.9 mm for the physical lung phantom. For the computational prostate phantoms, the corresponding mean error was 1.0–1.9 mm in the prostate, 1.9–2.4 mm in the rectum, and 1.8–2.1 mm over the entire patient body. Sinusoidal errors induced by B‐spline interpolations were observed in all the displacement profiles of the DMP registrations. Regions of large displacements were observed to have more registration errors. Patient‐specific FEM models have been developed to evaluate the DIR algorithms implemented in the commercial software package. It has been found that the accuracy of these algorithms is patient‐dependent and related to various factors including tissue deformation magnitudes and image intensity gradients across the regions of interest. This may suggest that DIR algorithms need to be verified for each registration instance when implementing adaptive radiation therapy. PACS numbers: 87.10.Kn, 87.55.km, 87.55.Qr, 87.57.nj

Clinical commissioning and use of the Novalis Tx linear accelerator for SRS and SBRT

The purpose of this study was to perform comprehensive measurements and testing of a Novalis Tx linear accelerator, and to develop technical guidelines for commissioning from the time of acceptance testing to the first clinical treatment. The Novalis Tx (NTX) linear accelerator is equipped with, among other features, a high‐definition MLC (HD120 MLC) with 2.5 mm central leaves, a 6D robotic couch, an optical guidance positioning system, as well as X‐ray‐based image guidance tools to provide high accuracy radiation delivery for stereotactic radiosurgery and stereotactic body radiation therapy procedures. We have performed extensive tests for each of the components, and analyzed the clinical data collected in our clinic. We present technical guidelines in this report focusing on methods for: (1) efficient and accurate beam data collection for commissioning treatment planning systems, including small field output measurements conducted using a wide range of detectors; (2) commissioning tests for the HD120 MLC; (3) data collection for the baseline characteristics of the on‐board imager (OBI) and ExacTrac X‐ray (ETX) image guidance systems in conjunction with the 6D robotic couch; and (4) end‐to‐end testing of the entire clinical process. Established from our clinical experience thus far, recommendations are provided for accurate and efficient use of the OBI and ETX localization systems for intra‐ and extracranial treatment sites. Four results are presented. (1) Basic beam data measurements: Our measurements confirmed the necessity of using small detectors for small fields. Total scatter factors varied significantly (30% to approximately 62%) for small field measurements among detectors. Unshielded stereotactic field diode (SFD) overestimated dose by ~ 2% for large field sizes. Ion chambers with active diameters of 6 mm suffered from significant volume averaging. The sharpest profile penumbra was observed for the SFD because of its small active diameter (0.6 mm). (2) MLC commissioning: Winston Lutz test, light/radiation field congruence, and Picket Fence tests were performed and were within criteria established by the relevant task group reports. The measured mean MLC transmission and dynamic leaf gap of 6 MV SRS beam were 1.17% and 0.36 mm, respectively. (3) Baseline characteristics of OBI and ETX: The isocenter localization errors in the left/right, posterior/anterior, and superior/inferior directions were, respectively, −0.2±0.2 mm, −0.8±0.2 mm, and −0.8±0.4 mm for ETX, and 0.5±0.7 mm, 0.6±0.5 mm, and 0.0±0.5 mm for OBI cone‐beam computed tomography. The registration angular discrepancy was 0.1±0.2°, and the maximum robotic couch error was 0.2°. (4) End‐to‐end tests: The measured isocenter dose differences from the planned values were 0.8% and 0.4%, measured respectively by an ion chamber and film. The gamma pass rate, measured by EBT2 film, was 95% (3% DD and 1 mm DTA). Through a systematic series of quantitative commissioning experiments and end‐to‐end tests and our initial clinical experience, described in this report, we demonstrate that the NTX is a robust system, with the image guidance and MLC requirements to treat a wide variety of sites — in particular for highly accurate delivery of SRS and SBRT‐based treatments. PACS numbers: 87.55.Qr, 87.53.Ly, 87.59.‐e

Correlation of dose computed using different algorithms with local control following stereotactic ablative radiotherapy (SABR)-based treatment of non-small-cell lung cancer

PURPOSE To retrospectively compute dose distributions for lung cancer patients treated with SABR, and to correlate dose distributions with outcome using a tumor control probability (TCP) model. METHODS Treatment plans for 133 NSCLC patients treated using 12 Gy/fxn × 4 (BED=106 Gy), and planned using a pencil-beam (1D-equivalent-path-length, EPL-1D) algorithm were retrospectively re-calculated using model-based algorithms (including convolution/superposition, Monte Carlo). 4D imaging was performed to manage motion. TCP was computed using the Marsden model and associations between dose and outcome were inferred. RESULTS Mean D95 reductions of 20% (max.=33%) were noted with model-based algorithms (relative to EPL-1D) for the smallest tumors (PTV<20 cm(3)), corresponding to actual delivered D95 BEDs of ≈ 60-85 Gy. For larger tumors (PTV>100 cm(3)), D95 reductions were ≈ 10% (BED>100 Gy). Mean lung doses (MLDs) were 15% lower for model-based algorithms for PTVs<20 cm(3). No correlation between tumor size and 2-year local control rate was observed clinically, consistent with TCP calculations, both of which were ≈ 90% across all PTV bins. CONCLUSION Results suggest that similar control rates might be achieved for smaller tumors using lower BEDs relative to larger tumors. However, more studies with larger patient cohorts are necessary to confirm this possible finding.

Characteristics of a novel treatment system for linear accelerator–based stereotactic radiosurgery

The purpose of this study is to characterize the dosimetric properties and accuracy of a novel treatment platform (Edge radiosurgery system) for localizing and treating patients with frameless, image‐guided stereotactic radiosurgery (SRS) and stereotactic body radiotherapy (SBRT). Initial measurements of various components of the system, such as a comprehensive assessment of the dosimetric properties of the flattening filter‐free (FFF) beams for both high definition (HD120) MLC and conical cone‐based treatment, positioning accuracy and beam attenuation of a six degree of freedom (6DoF) couch, treatment head leakage test, and integrated end‐to‐end accuracy tests, have been performed. The end‐to‐end test of the system was performed by CT imaging a phantom and registering hidden targets on the treatment couch to determine the localization accuracy of the optical surface monitoring system (OSMS), cone‐beam CT (CBCT), and MV imaging systems, as well as the radiation isocenter targeting accuracy. The deviations between the percent depth‐dose curves acquired on the new linac‐based system (Edge), and the previously published machine with FFF beams (TrueBeam) beyond Dmax were within 1.0% for both energies. The maximum deviation of output factors between the Edge and TrueBeam was 1.6%. The optimized dosimetric leaf gap values, which were fitted using Eclipse dose calculations and measurements based on representative spine radiosurgery plans, were 0.700 mm and 1.000 mm, respectively. For the conical cones, 6X FFF has sharper penumbra ranging from 1.2−1.8 mm (80%‐20%) and 1.9−3.8 mm (90%‐10%) relative to 10X FFF, which has 1.2−2.2 mm and 2.3−5.1 mm, respectively. The relative attenuation measurements of the couch for PA, PA (rails‐in), oblique, oblique (rails‐out), oblique (rails‐in) were: −2.0%, −2.5%, −15.6%, −2.5%, −5.0% for 6X FFF and −1.4%, −1.5%, −12.2%, −2.5%, −5.0% for 10X FFF, respectively, with a slight decrease in attenuation versus field size. The systematic deviation between the OSMS and CBCT was −0.4±0.2 mm, 0.1±0.3 mm, and 0.0±0.1 mm in the vertical, longitudinal, and lateral directions. The mean values and standard deviations of the average deviation and maximum deviation of the daily Winston‐Lutz tests over three months are 0.20±0.03 mm and 0.66±0.18 mm, respectively. Initial testing of this novel system demonstrates the technology to be highly accurate and suitable for frameless, linac‐based SRS and SBRT treatment. PACS number: 87.56.J‐

Contouring variability of human- and deformable-generated contours in radiotherapy for prostate cancer

This study was designed to evaluate contouring variability of human-and deformable-generated contours on planning CT (PCT) and CBCT for ten patients with low-or intermediate-risk prostate cancer. For each patient in this study, five radiation oncologists contoured the prostate, bladder, and rectum, on one PCT dataset and five CBCT datasets. Consensus contours were generated using the STAPLE method in the CERR software package. Observer contours were compared to consensus contour, and contour metrics (Dice coefficient, Hausdorff distance, Contour Distance, Center-of-Mass [COM] Deviation) were calculated. In addition, the first day CBCT was registered to subsequent CBCT fractions (CBCTn: CBCT2-CBCT5) via B-spline Deformable Image Registration (DIR). Contours were transferred from CBCT1 to CBCTn via the deformation field, and contour metrics were calculated through comparison with consensus contours generated from human contour set. The average contour metrics for prostate contours on PCT and CBCT were as follows: Dice coefficient-0.892 (PCT), 0.872 (CBCT-Human), 0.824 (CBCT-Deformed); Hausdorff distance-4.75 mm (PCT), 5.22 mm (CBCT-Human), 5.94 mm (CBCT-Deformed); Contour Distance (overall contour)-1.41 mm (PCT), 1.66 mm (CBCT-Human), 2.30 mm (CBCT-Deformed); COM Deviation-2.01 mm (PCT), 2.78 mm (CBCT-Human), 3.45 mm (CBCT-Deformed). For human contours on PCT and CBCT, the difference in average Dice coefficient between PCT and CBCT (approx. 2%) and Hausdorff distance (approx. 0.5 mm) was small compared to the variation between observers for each patient (standard deviation in Dice coefficient of 5% and Hausdorff distance of 2.0 mm). However, additional contouring variation was found for the deformable-generated contours (approximately 5.0% decrease in Dice coefficient and 0.7 mm increase in Hausdorff distance relative to human-generated contours on CBCT). Though deformable contours provide a reasonable starting point for contouring on CBCT, we conclude that contours generated with B-Spline DIR require physician review and editing if they are to be used in the clinic.