Automated testing platform for radiotherapy treatment planning scripts.

Realizing the potential of user-developed automation software interacting with a treatment planning system (TPS) requires rigorous testing to ensure patient safety and data integrity. We developed an automated test platform to allow comparison of the treatment planning database before and after the execution of a write-enabled script interacting with a commercial TPS (Eclipse, Varian Medical Systems, Palo Alto, CA) using the vendor-provided Eclipse Scripting Application Programming Interface (ESAPI). The C#-application known as Write-Enable Script Testing Engine (WESTE) serializes the treatment planning objects (Patient, Structure Set, PlanSetup) accessible through ESAPI, and then compares the serialization acquired before and after the execution of the script being tested, documenting identified differences to highlight the changes made to the treatment planning data. The first two uses of WESTE demonstrated that the testing platform could acquire and analyze the data quickly (<4 s per test case) and facilitate the clinical implementation of write-enabled scripts.

Implementation of free breathing respiratory amplitude-gated treatments.

PURPOSE: The purpose of this study was to provide guidance in developing and implementing a process for the accurate delivery of free breathing respiratory amplitude-gated treatments. METHODS: A phase-based 4DCT scan is acquired at time of simulation and motion is evaluated to determine the exhale phases that minimize respiratory motion to an acceptable level. A phase subset average CT is then generated for treatment planning and a tracking structure is contoured to indicate the location of the target or a suitable surrogate over the planning phases. Prior to treatment delivery, a 4DCBCT is acquired and a phase subset average is created to coincide with the planning phases for an initial match to the planning CT. Fluoroscopic imaging is then used to set amplitude gate thresholds corresponding to when the target or surrogate is in the tracking structure. The final imaging prior to treatment is an amplitude-gated CBCT to verify both the amplitude gate thresholds and patient positioning. An amplitude-gated treatment is then delivered. This technique was commissioned using an in-house lung motion phantom and film measurements of a simple two-field 3D plan. RESULTS: The accuracy of 4DCBCT motion and target position measurements were validated relative to 4DCT imaging. End to end testing showed strong agreement between planned and film measured dose distributions. Robustness to interuser variability and changes in respiratory motion were demonstrated through film measurements. CONCLUSIONS: The developed workflow utilizes 4DCBCT, respiratory-correlated fluoroscopy, and gated CBCT imaging in an efficient and sequential process to ensure the accurate delivery of free breathing respiratory-gated treatments.

Validation of clinical acceptability of deep-learning-based automated segmentation of organs-at-risk for head-and-neck radiotherapy treatment planning.

Introduction: Organ-at-risk segmentation for head and neck cancer radiation therapy is a complex and time-consuming process (requiring up to 42 individual structure, and may delay start of treatment or even limit access to function-preserving care. Feasibility of using a deep learning (DL) based autosegmentation model to reduce contouring time without compromising contour accuracy is assessed through a blinded randomized trial of radiation oncologists (ROs) using retrospective, de-identified patient data. Methods: Two head and neck expert ROs used dedicated time to create gold standard (GS) contours on computed tomography (CT) images. 445 CTs were used to train a custom 3D U-Net DL model covering 42 organs-at-risk, with an additional 20 CTs were held out for the randomized trial. For each held-out patient dataset, one of the eight participant ROs was randomly allocated to review and revise the contours produced by the DL model, while another reviewed contours produced by a medical dosimetry assistant (MDA), both blinded to their origin. Time required for MDAs and ROs to contour was recorded, and the unrevised DL contours, as well as the RO-revised contours by the MDAs and DL model were compared to the GS for that patient. Results: Mean time for initial MDA contouring was 2.3 hours (range 1.6-3.8 hours) and RO-revision took 1.1 hours (range, 0.4-4.4 hours), compared to 0.7 hours (range 0.1-2.0 hours) for the RO-revisions to DL contours. Total time reduced by 76% (95%-Confidence Interval: 65%-88%) and RO-revision time reduced by 35% (95%-CI,-39%-91%). All geometric and dosimetric metrics computed, agreement with GS was equivalent or significantly greater (p<0.05) for RO-revised DL contours compared to the RO-revised MDA contours, including volumetric Dice similarity coefficient (VDSC), surface DSC, added path length, and the 95%-Hausdorff distance. 32 OARs (76%) had mean VDSC greater than 0.8 for the RO-revised DL contours, compared to 20 (48%) for RO-revised MDA contours, and 34 (81%) for the unrevised DL OARs. Conclusion: DL autosegmentation demonstrated significant time-savings for organ-at-risk contouring while improving agreement with the institutional GS, indicating comparable accuracy of DL model. Integration into the clinical practice with a prospective evaluation is currently underway.

Quality assurance for image-guided radiation therapy utilizing CT-based technologies: a report of the AAPM TG-179.

PURPOSE: Commercial CT-based image-guided radiotherapy (IGRT) systems allow widespread management of geometric variations in patient setup and internal organ motion. This document provides consensus recommendations for quality assurance protocols that ensure patient safety and patient treatment fidelity for such systems. METHODS: The AAPM TG-179 reviews clinical implementation and quality assurance aspects for commercially available CT-based IGRT, each with their unique capabilities and underlying physics. The systems described are kilovolt and megavolt cone-beam CT, fan-beam MVCT, and CT-on-rails. A summary of the literature describing current clinical usage is also provided. RESULTS: This report proposes a generic quality assurance program for CT-based IGRT systems in an effort to provide a vendor-independent program for clinical users. Published data from long-term, repeated quality control tests form the basis of the proposed test frequencies and tolerances. CONCLUSION: A program for quality control of CT-based image-guidance systems has been produced, with focus on geometry, image quality, image dose, system operation, and safety. Agreement and clarification with respect to reports from the AAPM TG-101, TG-104, TG-142, and TG-148 has been addressed.

Novel dosimetric phantom for quality assurance of volumetric modulated arc therapy.

The objective of this work is to assess the suitability and performance of a new dosimeter system with a novel geometry for the quality assurance (QA) of volumetric modulated arc therapy (VMAT). The new dosimeter system consists of a hollow cylinder (15 and 25 cm inner and outer diameters) with 124 diodes embedded in the phantom's cylindrical wall forming four rings of detectors. For coplanar beams, the cylindrical geometry and the ring diode pattern offer the advantage of invariant perpendicular incidence on the beam central axis for any gantry angle and also have the benefit of increasing the detector density as both walls of the cylinder sample the beam. Other advantages include real-time readout and reduced weight with the hollow phantom shape. A calibration method taking into account the variation in radiation sensitivity of the diodes as a function of gantry angle was developed and implemented. In this work, the new dosimeter system was used in integrating mode to perform composite dose measurements along the cylindrical surface supporting the diodes. The reproducibility of the dosimeter response and the angular dependence of the diodes were assessed using simple 6 MV photon static beams. The performance of the new dosimeter system for VMAT QA was then evaluated using VMAT plans designed for a head and neck, an abdominal sarcoma, and a prostate patient. These plans were optimized with 90 control points (CPs) and additional versions of each plan were generated by increasing the number of CPs to 180 and 360 using linear interpolation. The relative dose measured with the dosimeter system for the VMAT plans was compared to the corresponding TPS dose map in terms of relative dose difference (% deltaD) and distance to agreement (DTA). The dosimeter system's sensitivity to gantry rotation offset and scaling errors as well as setup errors was also evaluated. For static beams, the dosimeter system offered good reproducibility and demonstrated small residual diode angular dependence after calibration. For VMAT deliveries, the agreement between measured and calculated doses was good with > or = 86.4% of the diodes satisfying 3% of % deltaD or 2 mm DTA for the 180 CP plans. The phantom offered sufficient sensitivity for the detection of small gantry rotation offset (3 degrees) and scaling errors (1 degree) as well as phantom setup errors of 1 mm, although the results were plan dependent. With its novel geometry, the dosimeter system was also able to experimentally demonstrate the discretization effect of the number of CPs used in the TPS to simulate a continuous arc. These results demonstrate the suitability of the new dosimeter system for the patient-specific QA of VMAT plans and suggest that the dosimeter system can be an effective tool in the routine QA and commissioning of treatment machines capable of VMAT delivery and cone-beam CT image guidance.

2D-3D registration for prostate radiation therapy based on a statistical model of transmission images.

PURPOSE: In external beam radiation therapy of pelvic sites, patient setup errors can be quantified by registering 2D projection radiographs acquired during treatment to a 3D planning computed tomograph (CT). We present a 2D-3D registration framework based on a statistical model of the intensity values in the two imaging modalities. METHODS: The model assumes that intensity values in projection radiographs are independently but not identically distributed due to the nonstationary nature of photon counting noise. Two probability distributions are considered for the intensity values: Poisson and Gaussian. Using maximum likelihood estimation, two similarity measures, maximum likelihood with a Poisson (MLP) and maximum likelihood with Gaussian (MLG), distribution are derived. Further, we investigate the merit of the model-based registration approach for data obtained with current imaging equipment and doses by comparing the performance of the similarity measures derived to that of the Pearson correlation coefficient (ICC) on accurately collected data of an anthropomorphic phantom of the pelvis and on patient data. RESULTS: Registration accuracy was similar for all three similarity measures and surpassed current clinical requirements of 3 mm for pelvic sites. For pose determination experiments with a kilovoltage (kV) cone-beam CT (CBCT) and kV projection radiographs of the phantom in the anterior-posterior (AP) view, registration accuracies were 0.42 mm (MLP), 0.29 mm (MLG), and 0.29 mm (ICC). For kV CBCT and megavoltage (MV) AP portal images of the same phantom, registration accuracies were 1.15 mm (MLP), 0.90 mm (MLG), and 0.69 mm (ICC). Registration of a kV CT and MV AP portal images of a patient was successful in all instances. CONCLUSIONS: The results indicate that high registration accuracy is achievable with multiple methods including methods that are based on a statistical model of a 3D CT and 2D projection images.

Accuracy of automatic couch corrections with on-line volumetric imaging.

The purpose of this study was to characterize automatic remote couch adjustment and to assess the accuracy of automatic couch corrections following localization with cone-beam CT (CBCT). Automatic couch movement was evaluated through passive reflector markers placed on a phantom, tracked with an optical tracking system (OTS). Repeated couch movements in the lateral, cranial/caudal, and vertical directions were monitored through the OTS to assess velocity and response time. In conjunction with CBCT, remote table movement for patient displacements following initial setup was available on four linear accelerators (Elekta Synergy). After the initial CBCT scan assessment, patients with isocenter displacements that exceeded clinical protocol tolerances were corrected using remote automatic couch movement. A verification CBCT scan was acquired after any remote movements. These verification CBCT datasets were assessed for the following time periods: one month post clinical installation, and six months later to monitor remote couch correction stability. Residual error analysis was evaluated using the verification scans. The mean +/- standard deviations (mu +/- sigma) of couch movement based on phantom measurements with the OTS were 0.16 +/- 0.48 mm, 0.32 +/- 0.30 mm, 0.11 +/- 0.12 mm in the L/R, A/P, and S/I couch directions, respectively. The fastest maximum velocity was observed in the inferior direction at 10.5 mm/s, and the slowest maximum velocity in the left direction at 3.6 mm/s. From 1134 verification CBCT registrations for 207 patients, the residual error for each translational direction from each month evaluated are reported. The mu was less than 0.3 mm in all directions, and sigma was in the order of 1 mm. At a 3 mm threshold, 21 of the 1134 fractions (2%) exceeded tolerance, attributed to patient intrafraction movement. Remote automatic couch movement is reliable and effective for adjusting patient position with a precision of approximately 1mm. Patient residual error observed in this study demonstrates that displacement is minimal after remote couch adjustment.

A quality assurance program for image quality of cone-beam CT guidance in radiation therapy.

The clinical introduction of volumetric x-ray image-guided radiotherapy systems necessitates formal commissioning of the hardware and image-guided processes to be used and drafts quality assurance (QA) for both hardware and processes. Satisfying both requirements provides confidence on the system's ability to manage geometric variations in patient setup and internal organ motion. As these systems become a routine clinical modality, the authors present data from their QA program tracking the image quality performance of ten volumetric systems over a period of 3 years. These data are subsequently used to establish evidence-based tolerances for a QA program. The volumetric imaging systems used in this work combines a linear accelerator with conventional x-ray tube and an amorphous silicon flat-panel detector mounted orthogonally from the accelerator central beam axis, in a cone-beam computed tomography (CBCT) configuration. In the spirit of the AAPM Report No. 74, the present work presents the image quality portion of their QA program; the aspects of the QA protocol addressing imaging geometry have been presented elsewhere. Specifically, the authors are presenting data demonstrating the high linearity of CT numbers, the uniformity of axial reconstructions, and the high contrast spatial resolution of ten CBCT systems (1-2 mm) from two commercial vendors. They are also presenting data accumulated over the period of several months demonstrating the long-term stability of the flat-panel detector and of the distances measured on reconstructed volumetric images. Their tests demonstrate that each specific CBCT system has unique performance. In addition, scattered x rays are shown to influence the imaging performance in terms of spatial resolution, axial reconstruction uniformity, and the linearity of CT numbers.

Statistical process control for IMRT dosimetric verification.

Patient-specific measurements are typically used to validate the dosimetry of intensity-modulated radiotherapy (IMRT). To evaluate the dosimetric performance over time of our IMRT process, we have used statistical process control (SPC) concepts to analyze the measurements from 330 head and neck (H&N) treatment plans. The objectives of the present work are to: (i) Review the dosimetric measurements of a large series of consecutive head and neck treatment plans to better understand appropriate dosimetric tolerances; (ii) analyze the results with SPC to develop action levels for measured discrepancies; (iii) develop estimates for the number of measurements that are required to describe IMRT dosimetry in the clinical setting; and (iv) evaluate with SPC a new beam model in our planning system. H&N IMRT cases were planned with the PINNACLE treatment planning system versions 6.2b or 7.6c (Philips Medical Systems, Madison, WI) and treated on Varian (Palo Alto, CA) or Elekta (Crawley, UK) linacs. As part of regular quality assurance, plans were recalculated on a 20-cm-diam cylindrical phantom, and ion chamber measurements were made in high-dose volumes (the PTV with highest dose) and in low-dose volumes (spinal cord organ-at-risk, OR). Differences between the planned and measured doses were recorded as a percentage of the planned dose. Differences were stable over time. Measurements with PINNACLE3 6.2b and Varian linacs showed a mean difference of 0.6% for PTVs (n=149, range, -4.3% to 6.6%), while OR measurements showed a larger systematic discrepancy (mean 4.5%, range -4.5% to 16.3%) that was due to well-known limitations of the MLC model in the earlier version of the planning system. Measurements with PINNACLE3 7.6c and Varian linacs demonstrated a mean difference of 0.2% for PTVs (n=160, range, -3.0%, to 5.0%) and -1.0% for ORs (range -5.8% to 4.4%). The capability index (ratio of specification range to range of the data) was 1.3 for the PTV data, indicating that almost all measurements were within +/-5%. We have used SPC tools to evaluate a new beam model in our planning system to produce a systematic difference of -0.6% for PTVs and 0.4% for ORs, although the number of measurements is smaller (n=25). Analysis of this large series of H&N IMRT measurements demonstrated that our IMRT dosimetry was stable over time and within accepted tolerances. These data provide useful information for assessing alterations to beam models in the planning system. IMRT is enhanced by the addition of statistical process control to traditional quality control procedures.

Automated 2D-3D registration of portal images and CT data using line-segment enhancement.

In prostate radiotherapy, setup errors with respect to the patient's bony anatomy can be reduced by aligning 2D megavoltage (MV) portal images acquired during treatment to a reference 3D kilovoltage (kV) CT acquired for treatment planning purposes. The purpose of this study was to evaluate a fully automated 2D-3D registration algorithm to quantify setup errors in 3D through the alignment of line-enhanced portal images and digitally reconstructed radiographs computed from the CT. The line-enhanced images were obtained by correlating the images with a filter bank of short line segments, or "sticks" at different orientations. The proposed methods were validated on (1) accurately collected gold-standard data consisting of a 3D kV cone-beam CT scan of an anthropomorphic phantom of the pelvis and 2D MV portal images in the anterior-posterior (AP) view acquired at 15 different poses and (2) a conventional 3D kV CT scan and weekly 2D MV AP portal images of a patient over 8 weeks. The mean (and standard deviation) of the absolute registration error for rotations around the right-lateral (RL), inferior-superior (IS), and posterior-anterior (PA) axes were 0.212 degree (0.214 degree), 0.055 degree (0.033 degree) and 0.041 degree (0.039 degree), respectively. The corresponding registration errors for translations along the RL, IS, and PA axes were 0.161 (0.131) mm, 0.096 (0.033) mm, and 0.612 (0.485) mm. The mean (and standard deviation) of the total registration error was 0.778 (0.543) mm. Registration on the patient images was successful in all eight cases as determined visually. The results indicate that it is feasible to automatically enhance features in MV portal images of the pelvis for use within a completely automated 2D-3D registration framework for the accurate determination of patient setup errors. They also indicate that it is feasible to estimate all six transformation parameters from a 3D CT of the pelvis and a single portal image in the AP view.

Improving image-guided target localization through deformable registration.

PURPOSE: To quantify the improvements in online target localization using kV cone beam CT (CBCT) with deformable registration. METHODS AND MATERIAL: Twelve patients treated under a 6 fraction liver cancer radiation therapy protocol were imaged in breath hold using kV CBCT at each treatment fraction. The images were imported into the treatment planning software and rigidly registered by fitting the liver, identified on the daily kV CBCT image, into the liver contours, previously drawn on the planning CT. The liver was then manually contoured on each CBCT image. Deformable registration was automatically performed, aligning the CT liver to the liver on each CBCT image using MORFEUS, a biomechanical model based deformable registration algorithm. The tumor, defined on planning CT, was mapped onto the CBCT, through MORFEUS. The center of mass (COM) displacement of the tumor was computed. RESULTS: The mean (SD) displacement magnitude (absolute value) of the COM following deformable registration was 0.08 (0.07), 0.10 (0.11), and 0.10 (0.17) cm in the left-right (LR), anterior-posterior (AP), and superior-inferior (SI) directions, respectively. The maximum displacement of the COM was 0.34, 0.65, and 0.97 cm in the LR, AP, and SI directions, respectively. Fifteen percent of the treatment fractions had a COM displacement of greater than 0.3 cm and 33% of patients had at least 1 fraction with a displacement of greater than 0.3 cm. The deformable registration, excluding the manual contouring of the liver, was performed in less than 1 minute, on average. DISCUSSION: Rigid registration of the liver volume between planning CT and verification kV CBCT localizes the tumor to within 0.3 cm for the majority (66%) of patients; however, larger offsets in tumor position can be observed due to liver deformation.

2D-3D registration for cranial radiation therapy using a 3D kV CBCT and a single limited field-of-view 2D kV radiograph.

PURPOSE: We present and evaluate a fully automated 2D-3D intensity-based registration framework using a single limited field-of-view (FOV) 2D kV radiograph and a 3D kV CBCT for 3D estimation of patient setup errors during brain radiotherapy. METHODS: We evaluated two similarity measures, the Pearson correlation coefficient on image intensity values (ICC) and maximum likelihood measure with Gaussian noise (MLG), derived from the statistics of transmission images. Pose determination experiments were conducted on 2D kV radiographs in the anterior-posterior (AP) and left lateral (LL) views and 3D kV CBCTs of an anthropomorphic head phantom. In order to minimize radiation exposure and exclude nonrigid structures from the registration, limited FOV 2D kV radiographs were employed. A spatial frequency band useful for the 2D-3D registration was identified from the bone-to-no-bone spectral ratio (BNBSR) of digitally reconstructed radiographs (DRRs) computed from the 3D kV planning CT of the phantom. The images being registered were filtered accordingly prior to computation of the similarity measures. We evaluated the registration accuracy achievable with a single 2D kV radiograph and with the registration results from the AP and LL views combined. We also compared the performance of the 2D-3D registration solutions proposed to that of a commercial 3D-3D registration algorithm, which used the entire skull for the registration. The ground truth was determined from markers affixed to the phantom and visible in the CBCT images. RESULTS: The accuracy of the 2D-3D registration solutions, as quantified by the root mean squared value of the target registration error (TRE) calculated over a radius of 3 cm for all poses tested, was ICCAP : 0.56 mm, MLGAP : 0.74 mm, ICCLL : 0.57 mm, MLGLL : 0.54 mm, ICC (AP and LL combined): 0.19 mm, and MLG (AP and LL combined): 0.21 mm. The accuracy of the 3D-3D registration algorithm was 0.27 mm. There was no significant difference in mean TRE for the 2D-3D registration algorithms using a single 2D kV radiograph with similarity measure and image view point. There was no significant difference in mean TRE between ICCLL , MLGLL , ICC (AP and LL combined), MLG (AP and LL combined), and the 3D-3D registration algorithm despite the smaller FOV used for the 2D-3D registration. While submillimeter registration accuracy was obtained with both ICC and MLG using a single 2D kV radiograph, combining the results from the two projection views resulted in a significantly smaller (P≤0.05) mean TRE. CONCLUSIONS: Our results indicate that it is possible to achieve submillimeter registration accuracy with both ICC and MLG using either single or dual limited FOV 2D kV radiographs of the head in the AP and LL views. The registration accuracy suggests that the 2D-3D registration solutions presented are suitable for the estimation of patient setup errors not only during conventional brain radiation therapy, but also during stereotactic procedures and proton radiation therapy where tighter setup margins are required.

A magnetic resonance imaging study of prostate deformation relative to implanted gold fiducial markers.

PURPOSE: To describe prostate deformation during radiotherapy and determine the margins required to account for prostate deformation after setup to intraprostatic fiducial markers (FM). METHODS AND MATERIALS: Twenty-five patients with T1c-T2c prostate cancer had three gold FMs implanted. The patients presented with a full bladder and empty rectum for two axial magnetic resonance imaging (MRI) scans using a gradient recalled echo (GRE) sequence capable of imaging the FMs. The MRIs were done at the time of radiotherapy (RT) planning and a randomly assigned fraction. A single observer contoured the prostate surfaces. They were entered into a finite element model and aligned using the centroid of the three FMs. RESULTS: During RT, the prostate volume decreased by 0.5%/fraction (p = 0.03) and the FMs in-migrated by 0.05 mm/fraction (p < 0.05). Prostate deformation was unrelated to differential bladder and bowel filling, but was related to a transurethral resection of the prostate (TURP) (p = 0.003). The standard deviation for systematic uncertainty of prostate surface contouring was 0.8 mm and for FM centroid localization was 0.4 mm. The standard deviation of random interfraction prostate deformation was 1.5 mm and for FM centroid variability was 1.1 mm. These uncertainties from prostate deformation can be incorporated into a margin recipe to determine the total margins required for RT. CONCLUSIONS: During RT, the prostate exhibited: volume decrease, deformation, and in-migration of FMs. Patients with TURPs were prone to prostate deformation.

Comparison of localization performance with implanted fiducial markers and cone-beam computed tomography for on-line image-guided radiotherapy of the prostate.

PURPOSE: The aim of this work was to assess the accuracy of kilovoltage (kV) cone-beam computed tomography (CBCT)-based setup corrections as compared with orthogonal megavoltage (MV) portal image-based corrections for patients undergoing external-beam radiotherapy of the prostate. METHODS AND MATERIALS: Daily cone-beam CT volumetric images were acquired after setup for patients with three intraprostatic fiducial markers. The estimated couch shifts were compared retrospectively to patient adjustments based on two orthogonal MV portal images (the current clinical standard of care in our institution). The CBCT soft-tissue based shifts were also estimated by digitally removing the gold markers in each projection to suppress the artifacts in the reconstructed volumes. A total of 256 volumetric images for 15 patients were analyzed. RESULTS: The Pearson coefficient of correlation for the patient position shifts using fiducial markers in MV vs. kV was (R2 = 0.95, 0.84, 0.81) in the left-right (LR), anterior-posterior (AP), and superior-inferior (SI) directions, respectively. The correlation using soft-tissue matching was as follows: R2 = 0.90, 0.49, 0.51 in the LR, AP and SI directions. A Bland-Altman analysis showed no significant trends in the data. The percentage of shifts within a +/-3-mm tolerance (the clinical action level) was 99.7%, 95.5%, 91.3% for fiducial marker matching and 99.5%, 70.3%, 78.4% for soft-tissue matching. CONCLUSIONS: Cone-beam CT is an accurate and precise tool for image guidance. It provides an equivalent means of patient setup correction for prostate patients with implanted gold fiducial markers. Use of the additional information provided by the visualization of soft-tissue structures is an active area of research.

A frequency-based approach to locate common structure for 2D-3D intensity-based registration of setup images in prostate radiotherapy.

In many radiotherapy clinics, geometric uncertainties in the delivery of 3D conformal radiation therapy and intensity modulated radiation therapy of the prostate are reduced by aligning the patient's bony anatomy in the planning 3D CT to corresponding bony anatomy in 2D portal images acquired before every treatment fraction. In this paper, we seek to determine if there is a frequency band within the portal images and the digitally reconstructed radiographs (DRRs) of the planning CT in which bony anatomy predominates over non-bony anatomy such that portal images and DRRs can be suitably filtered to achieve high registration accuracy in an automated 2D-3D single portal intensity-based registration framework. Two similarity measures, mutual information and the Pearson correlation coefficient were tested on carefully collected gold-standard data consisting of a kilovoltage cone-beam CT (CBCT) and megavoltage portal images in the anterior-posterior (AP) view of an anthropomorphic phantom acquired under clinical conditions at known poses, and on patient data. It was found that filtering the portal images and DRRs during the registration considerably improved registration performance. Without filtering, the registration did not always converge while with filtering it always converged to an accurate solution. For the pose-determination experiments conducted on the anthropomorphic phantom with the correlation coefficient, the mean (and standard deviation) of the absolute errors in recovering each of the six transformation parameters were Theta(x):0.18(0.19) degrees, Theta(y):0.04(0.04) degrees, Theta(z):0.04(0.02) degrees, t(x):0.14(0.15) mm, t(y):0.09(0.05) mm, and t(z):0.49(0.40) mm. The mutual information-based registration with filtered images also resulted in similarly small errors. For the patient data, visual inspection of the superimposed registered images showed that they were correctly aligned in all instances. The results presented in this paper suggest that robust and accurate registration can be achieved with intensity-based methods by focusing on rigid bony structures in the images while diminishing the influence of artifacts with similar frequencies as soft tissue.

Image guidance: treatment target localization systems.

Highly conformal radiation therapy tailors treatment to match the target shape and position, minimizing normal tissue damage to a greater extent than previously possible. Technological advances such as intensity-modulated radiation therapy, introduced a decade ago, have yielded significant gains in tumor control and reduced toxicity. Continuing advances have focused on the characterization and control of patient movement, organ motion, and anatomical deformation, which all introduce geometric uncertainty. These sources of uncertainty limit the effectiveness of high-precision treatment. Target localization, performed using appropriate technologies and frequency, is a critical component of treatment quality assurance. Until recently, the target position with respect to the beams has been inferred from surface marks on the patient's skin or through an immobilization device, and verified using megavoltage radiographs of the treatment portal. Advances in imaging technologies have made it possible to image soft tissue volumes in the treatment setting. Real-time tracking is also possible using a variety of technologies, including fluoroscopic imaging and radiopaque markers implanted in or near the tumor. The capacity to acquire volumetric soft tissue images in the treatment setting can also be used to assess anatomical changes over a course of treatment. Enhancing localization practices reduces treatment errors, and gives the capacity to monitor anatomical changes and reduce uncertainties that could influence clinical outcomes. This review presents the technologies available for target localization, and discusses some of the considerations that should be addressed in the implementation of many new clinical processes in radiation oncology.

Feasibility of a novel deformable image registration technique to facilitate classification, targeting, and monitoring of tumor and normal tissue.

PURPOSE: To investigate the feasibility of a biomechanical-based deformable image registration technique for the integration of multimodality imaging, image guided treatment, and response monitoring. METHODS AND MATERIALS: A multiorgan deformable image registration technique based on finite element modeling (FEM) and surface projection alignment of selected regions of interest with biomechanical material and interface models has been developed. FEM also provides an inherent method for direct tracking specified regions through treatment and follow-up. RESULTS: The technique was demonstrated on 5 liver cancer patients. Differences of up to 1 cm of motion were seen between the diaphragm and the tumor center of mass after deformable image registration of exhale and inhale CT scans. Spatial differences of 5 mm or more were observed for up to 86% of the surface of the defined tumor after deformable image registration of the computed tomography (CT) and magnetic resonance images. Up to 6.8 mm of motion was observed for the tumor after deformable image registration of the CT and cone-beam CT scan after rigid registration of the liver. Deformable registration of the CT to the follow-up CT allowed a more accurate assessment of tumor response. CONCLUSIONS: This biomechanical-based deformable image registration technique incorporates classification, targeting, and monitoring of tumor and normal tissue using one methodology.

Characterization of scattered radiation in kV CBCT images using Monte Carlo simulations.

Kilovoltage (kV) cone beam computed tomography (CBCT) images suffer from a substantial scatter contribution. In this study, Monte Carlo (MC) simulations are used to evaluate the scattered radiation present in projection images. These predicted scatter distributions are also used as a scatter correction technique. Images were acquired using a kV CBCT bench top system. The EGSnrc MC code was used to model the flat panel imager, the phantoms, and the x-ray source. The x-ray source model was validated using first and second half-value layers (HVL) and profile measurements. The HVLs and the profile were found to agree within 3% and 6%, respectively. MC simulated and measured projection images for a cylindrical water phantom and for an anthropomorphic head phantom agreed within 8% and 10%. A modified version of the DOSXYZnrc MC code was used to score phase space files with identified scattered and primary particles behind the phantoms. The cone angle, the source-to-detector distance, the phantom geometry, and the energy were varied to determine their effect on the scattered radiation distribution. A scatter correction technique was developed in which the MC predicted scatter distribution is subtracted from the projections prior to reconstruction. Preliminary testing of the procedure was done with an anthropomorphic head phantom and a contrast phantom. Contrast and profile measurements were obtained for the scatter corrected and noncorrected images. An improvement of 3% for contrast between solid water and a liver insert and 11% between solid water and a Teflon insert were obtained and a significant reduction in cupping and streaking artifacts was observed.

The stability of mechanical calibration for a kV cone beam computed tomography system integrated with linear accelerator.

The geometric accuracy and precision of an image-guided treatment system were assessed. Image guidance is performed using an x-ray volume imaging (XVI) system integrated with a linear accelerator and treatment planning system. Using an amorphous silicon detector and x-ray tube, volumetric computed tomography images are reconstructed from kilovoltage radiographs by filtered backprojection. Image fusion and assessment of geometric targeting are supported by the treatment planning system. To assess the limiting accuracy and precision of image-guided treatment delivery, a rigid spherical target embedded in an opaque phantom was subjected to 21 treatment sessions over a three-month period. For each session, a volumetric data set was acquired and loaded directly into an active treatment planning session. Image fusion was used to ascertain the couch correction required to position the target at the prescribed iso-center. Corrections were validated independently using megavoltage electronic portal imaging to record the target position with respect to symmetric treatment beam apertures. An initial calibration cycle followed by repeated image-guidance sessions demonstrated the XVI system could be used to relocate an unambiguous object to within less than 1 mm of the prescribed location. Treatment could then proceed within the mechanical accuracy and precision of the delivery system. The calibration procedure maintained excellent spatial resolution and delivery precision over the duration of this study, while the linear accelerator was in routine clinical use. Based on these results, the mechanical accuracy and precision of the system are ideal for supporting high-precision localization and treatment of soft-tissue targets.

Patient dose from kilovoltage cone beam computed tomography imaging in radiation therapy.

Kilovoltage cone-beam computerized tomography (kV-CBCT) systems integrated into the gantry of linear accelerators can be used to acquire high-resolution volumetric images of the patient in the treatment position. Using on-line software and hardware, patient position can be determined accurately with a high degree of precision and, subsequently, set-up parameters can be adjusted to deliver the intended treatment. While the patient dose due to a single volumetric imaging acquisition is small compared to the therapy dose, repeated and daily image guidance procedures can lead to substantial dose to normal tissue. The dosimetric properties of a clinical CBCT system have been studied on an Elekta linear accelerator (Synergy RP, XVI system) and additional measurements performed on a laboratory system with identical geometry. Dose measurements were performed with an ion chamber and MOSFET detectors at the center, periphery, and surface of 30 and 16-cm-diam cylindrical shaped water phantoms, as a function of x-ray energy and longitudinal field-of-view (FOV) settings of 5,10,15, and 26 cm. The measurements were performed for full 360 degrees CBCT acquisition as well as for half-rotation scans for 120 kVp beams using the 30-cm-diam phantom. The dose at the center and surface of the body phantom were determined to be 1.6 and 2.3 cGy for a typical imaging protocol, using full rotation scan, with a technique setting of 120 kVp and 660 mAs. The results of our measurements have been presented in terms of a dose conversion factor fCBCT, expressed in cGy/R. These factors depend on beam quality and phantom size as well as on scan geometry and can be utilized to estimate dose for any arbitrary mAs setting and reference exposure rate of the x-ray tube at standard distance. The results demonstrate the opportunity to manipulate the scanning parameters to reduce the dose to the patient by employing lower energy (kVp) beams, smaller FOV, or by using half-rotation scan.