Engineering topics

Li, Ruijiang

Publications and source records attributed to Li, Ruijiang.

14 records

DOE OSTI2017

Robust Estimation of Electron Density From Anatomic Magnetic Resonance Imaging of the Brain Using a Unifying Multi-Atlas Approach

Purpose: To develop a reliable method to estimate electron density based on anatomic magnetic resonance imaging (MRI) of the brain. Methods and Materials: We proposed a unifying multi-atlas approach for electron density estimation based on standard T1- and T2-weighted MRI. First, a composite atlas was constructed through a voxelwise matching process using multiple atlases, with the goal of mitigating effects of inherent anatomic variations between patients. Next we computed for each voxel 2 kinds of conditional probabilities: (1) electron density given its image intensity on T1- and T2-weighted MR images; and (2) electron density given its spatial location in a reference anatomy, obtained by deformable image registration. These were combined into a unifying posterior probability density function using the Bayesian formalism, which provided the optimal estimates for electron density. We evaluated the method on 10 patients using leave-one-patient-out cross-validation. Receiver operating characteristic analyses for detecting different tissue types were performed. Results: The proposed method significantly reduced the errors in electron density estimation, with a mean absolute Hounsfield unit error of 119, compared with 140 and 144 (P<.0001) using conventional T1-weighted intensity and geometry-based approaches, respectively. For detection of bony anatomy, the proposed method achieved an 89% area under the curve, 86% sensitivity, 88% specificity, and 90% accuracy, which improved upon intensity and geometry-based approaches (area under the curve: 79% and 80%, respectively). Conclusion: The proposed multi-atlas approach provides robust electron density estimation and bone detection based on anatomic MRI. If validated on a larger population, our work could enable the use of MRI as a primary modality for radiation treatment planning.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2016

Quantitative Analysis of {sup 18}F-Fluorodeoxyglucose Positron Emission Tomography Identifies Novel Prognostic Imaging Biomarkers in Locally Advanced Pancreatic Cancer Patients Treated With Stereotactic Body Radiation Therapy

Purpose: To identify prognostic biomarkers in pancreatic cancer using high-throughput quantitative image analysis. Methods and Materials: In this institutional review board–approved study, we retrospectively analyzed images and outcomes for 139 locally advanced pancreatic cancer patients treated with stereotactic body radiation therapy (SBRT). The overall population was split into a training cohort (n=90) and a validation cohort (n=49) according to the time of treatment. We extracted quantitative imaging characteristics from pre-SBRT {sup 18}F-fluorodeoxyglucose positron emission tomography, including statistical, morphologic, and texture features. A Cox proportional hazard regression model was built to predict overall survival (OS) in the training cohort using 162 robust image features. To avoid over-fitting, we applied the elastic net to obtain a sparse set of image features, whose linear combination constitutes a prognostic imaging signature. Univariate and multivariate Cox regression analyses were used to evaluate the association with OS, and concordance index (CI) was used to evaluate the survival prediction accuracy. Results: The prognostic imaging signature included 7 features characterizing different tumor phenotypes, including shape, intensity, and texture. On the validation cohort, univariate analysis showed that this prognostic signature was significantly associated with OS (P=.002, hazard ratio 2.74), which improved upon conventional imaging predictors including tumor volume, maximum standardized uptake value, and total legion glycolysis (P=.018-.028, hazard ratio 1.51-1.57). On multivariate analysis, the proposed signature was the only significant prognostic index (P=.037, hazard ratio 3.72) when adjusted for conventional imaging and clinical factors (P=.123-.870, hazard ratio 0.53-1.30). In terms of CI, the proposed signature scored 0.66 and was significantly better than competing prognostic indices (CI 0.48-0.64, Wilcoxon rank sum test P<1e-6). Conclusion: Quantitative analysis identified novel {sup 18}F-fluorodeoxyglucose positron emission tomography image features that showed improved prognostic value over conventional imaging metrics. If validated in large, prospective cohorts, the new prognostic signature might be used to identify patients for individualized risk-adaptive therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2016

Robust Intratumor Partitioning to Identify High-Risk Subregions in Lung Cancer: A Pilot Study

Purpose: To develop an intratumor partitioning framework for identifying high-risk subregions from {sup 18}F-fluorodeoxyglucose positron emission tomography (FDG-PET) and computed tomography (CT) imaging and to test whether tumor burden associated with the high-risk subregions is prognostic of outcomes in lung cancer. Methods and Materials: In this institutional review board–approved retrospective study, we analyzed the pretreatment FDG-PET and CT scans of 44 lung cancer patients treated with radiation therapy. A novel, intratumor partitioning method was developed, based on a 2-stage clustering process: first at the patient level, each tumor was over-segmented into many superpixels by k-means clustering of integrated PET and CT images; next, tumor subregions were identified by merging previously defined superpixels via population-level hierarchical clustering. The volume associated with each of the subregions was evaluated using Kaplan-Meier analysis regarding its prognostic capability in predicting overall survival (OS) and out-of-field progression (OFP). Results: Three spatially distinct subregions were identified within each tumor that were highly robust to uncertainty in PET/CT co-registration. Among these, the volume of the most metabolically active and metabolically heterogeneous solid component of the tumor was predictive of OS and OFP on the entire cohort, with a concordance index or CI of 0.66-0.67. When restricting the analysis to patients with stage III disease (n=32), the same subregion achieved an even higher CI of 0.75 (hazard ratio 3.93, log-rank P=.002) for predicting OS, and a CI of 0.76 (hazard ratio 4.84, log-rank P=.002) for predicting OFP. In comparison, conventional imaging markers, including tumor volume, maximum standardized uptake value, and metabolic tumor volume using threshold of 50% standardized uptake value maximum, were not predictive of OS or OFP, with CI mostly below 0.60 (log-rank P>.05). Conclusion: We propose a robust intratumor partitioning method to identify clinically relevant, high-risk subregions in lung cancer. We envision that this approach will be applicable to identifying useful imaging biomarkers in many cancer types.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2015

Optimization approaches to volumetric modulated arc therapy planning

Volumetric modulated arc therapy (VMAT) has found widespread clinical application in recent years. A large number of treatment planning studies have evaluated the potential for VMAT for different disease sites based on the currently available commercial implementations of VMAT planning. In contrast, literature on the underlying mathematical optimization methods used in treatment planning is scarce. VMAT planning represents a challenging large scale optimization problem. In contrast to fluence map optimization in intensity-modulated radiotherapy planning for static beams, VMAT planning represents a nonconvex optimization problem. In this paper, the authors review the state-of-the-art in VMAT planning from an algorithmic perspective. Different approaches to VMAT optimization, including arc sequencing methods, extensions of direct aperture optimization, and direct optimization of leaf trajectories are reviewed. Their advantages and limitations are outlined and recommendations for improvements are discussed.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2015

Simultaneous beam sampling and aperture shape optimization for SPORT

Purpose: Station parameter optimized radiation therapy (SPORT) was recently proposed to fully utilize the technical capability of emerging digital linear accelerators, in which the station parameters of a delivery system, such as aperture shape and weight, couch position/angle, gantry/collimator angle, can be optimized simultaneously. SPORT promises to deliver remarkable radiation dose distributions in an efficient manner, yet there exists no optimization algorithm for its implementation. The purpose of this work is to develop an algorithm to simultaneously optimize the beam sampling and aperture shapes. Methods: The authors build a mathematical model with the fundamental station point parameters as the decision variables. To solve the resulting large-scale optimization problem, the authors devise an effective algorithm by integrating three advanced optimization techniques: column generation, subgradient method, and pattern search. Column generation adds the most beneficial stations sequentially until the plan quality improvement saturates and provides a good starting point for the subsequent optimization. It also adds the new stations during the algorithm if beneficial. For each update resulted from column generation, the subgradient method improves the selected stations locally by reshaping the apertures and updating the beam angles toward a descent subgradient direction. The algorithm continues to improve the selected stations locally and globally by a pattern search algorithm to explore the part of search space not reachable by the subgradient method. By combining these three techniques together, all plausible combinations of station parameters are searched efficiently to yield the optimal solution. Results: A SPORT optimization framework with seamlessly integration of three complementary algorithms, column generation, subgradient method, and pattern search, was established. The proposed technique was applied to two previously treated clinical cases: a head and neck and a prostate case. It significantly improved the target conformality and at the same time critical structure sparing compared with conventional intensity modulated radiation therapy (IMRT). In the head and neck case, for example, the average PTV coverage D99% for two PTVs, cord and brainstem max doses, and right parotid gland mean dose were improved, respectively, by about 7%, 37%, 12%, and 16%. Conclusions: The proposed method automatically determines the number of the stations required to generate a satisfactory plan and optimizes simultaneously the involved station parameters, leading to improved quality of the resultant treatment plans as compared with the conventional IMRT plans.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

Accuracy of surface registration compared to conventional volumetric registration in patient positioning for head-and-neck radiotherapy: A simulation study using patient data

Purpose: 3D optical surface imaging has been applied to patient positioning in radiation therapy (RT). The optical patient positioning system is advantageous over conventional method using cone-beam computed tomography (CBCT) in that it is radiation free, frameless, and is capable of real-time monitoring. While the conventional radiographic method uses volumetric registration, the optical system uses surface matching for patient alignment. The relative accuracy of these two methods has not yet been sufficiently investigated. This study aims to investigate the theoretical accuracy of the surface registration based on a simulation study using patient data. Methods: This study compares the relative accuracy of surface and volumetric registration in head-and-neck RT. The authors examined 26 patient data sets, each consisting of planning CT data acquired before treatment and patient setup CBCT data acquired at the time of treatment. As input data of surface registration, patient’s skin surfaces were created by contouring patient skin from planning CT and treatment CBCT. Surface registration was performed using the iterative closest points algorithm by point–plane closest, which minimizes the normal distance between source points and target surfaces. Six degrees of freedom (three translations and three rotations) were used in both surface and volumetric registrations and the results were compared. The accuracy of each method was estimated by digital phantom tests. Results: Based on the results of 26 patients, the authors found that the average and maximum root-mean-square translation deviation between the surface and volumetric registrations were 2.7 and 5.2 mm, respectively. The residual error of the surface registration was calculated to have an average of 0.9 mm and a maximum of 1.7 mm. Conclusions: Surface registration may lead to results different from those of the conventional volumetric registration. Only limited accuracy can be achieved for patient positioning with an approach based solely on surface information.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

An initial study on the estimation of time-varying volumetric treatment images and 3D tumor localization from single MV cine EPID images

Purpose: In this work the authors develop and investigate the feasibility of a method to estimate time-varying volumetric images from individual MV cine electronic portal image device (EPID) images. Methods: The authors adopt a two-step approach to time-varying volumetric image estimation from a single cine EPID image. In the first step, a patient-specific motion model is constructed from 4DCT. In the second step, parameters in the motion model are tuned according to the information in the EPID image. The patient-specific motion model is based on a compact representation of lung motion represented in displacement vector fields (DVFs). DVFs are calculated through deformable image registration (DIR) of a reference 4DCT phase image (typically peak-exhale) to a set of 4DCT images corresponding to different phases of a breathing cycle. The salient characteristics in the DVFs are captured in a compact representation through principal component analysis (PCA). PCA decouples the spatial and temporal components of the DVFs. Spatial information is represented in eigenvectors and the temporal information is represented by eigen-coefficients. To generate a new volumetric image, the eigen-coefficients are updated via cost function optimization based on digitally reconstructed radiographs and projection images. The updated eigen-coefficients are then multiplied with the eigenvectors to obtain updated DVFs that, in turn, give the volumetric image corresponding to the cine EPID image. Results: The algorithm was tested on (1) Eight digital eXtended CArdiac-Torso phantom datasets based on different irregular patient breathing patterns and (2) patient cine EPID images acquired during SBRT treatments. The root-mean-squared tumor localization error is (0.73 ± 0.63 mm) for the XCAT data and (0.90 ± 0.65 mm) for the patient data. Conclusions: The authors introduced a novel method of estimating volumetric time-varying images from single cine EPID images and a PCA-based lung motion model. This is the first method to estimate volumetric time-varying images from single MV cine EPID images, and has the potential to provide volumetric information with no additional imaging dose to the patient.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

Assessing the Dosimetric Impact of Real-Time Prostate Motion During Volumetric Modulated Arc Therapy

Purpose: To develop a method for dose reconstruction by incorporating the interplay effect between aperture modulation and target motion, and to assess the dosimetric impact of real-time prostate motion during volumetric modulated arc therapy (VMAT). Methods and Materials: Clinical VMAT plans were delivered with the TrueBeam linac for 8 patients with prostate cancer. The real-time target motion during dose delivery was determined based on the 2-dimensional fiducial localization using an onboard electronic portal imaging device. The target shift in each image was correlated with the control point with the same gantry angle in the VMAT plan. An in-house-developed Monte Carlo simulation tool was used to calculate the 3-dimensional dose distribution for each control point individually, taking into account the corresponding real-time target motion (assuming a nondeformable target with no rotation). The delivered target dose was then estimated by accumulating the dose from all control points in the plan. On the basis of this information, dose–volume histograms and 3-dimensional dose distributions were calculated to assess their degradation from the planned dose caused by target motion. Thirty-two prostate motion trajectories were analyzed. Results: The minimum dose to 0.03 cm{sup 3} of the gross tumor volume (D{sub 0.03cc}) was only slightly degraded after taking motion into account, with a minimum value of 94.1% of the planned dose among all patients and fractions. However, the gross tumor volume receiving prescription dose (V{sub 100%}) could be largely affected by motion, dropping below 60% in 1 trajectory. We did not observe a correlation between motion magnitude and dose degradation. Conclusions: Prostate motion degrades the delivered dose to the target in an unpredictable way, although its effect is reduced over multiple fractions, and for most patients the degradation is small. Patients with greater prostate motion or those treated with stereotactic body radiation therapy would benefit from real-time prostate tracking to reduce the margin.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

Nonisocentric Treatment Strategy for Breast Radiation Therapy: A Proof of Concept Study

Purpose: To propose a nonisocentric treatment strategy as a special form of station parameter optimized radiation therapy, to improve sparing of critical structures while preserving target coverage in breast radiation therapy. Methods and Materials: To minimize the volume of exposed lung and heart in breast irradiation, we propose a novel nonisocentric treatment scheme by strategically placing nonconverging beams with multiple isocenters. As its name suggests, the central axes of these beams do not intersect at a single isocenter as in conventional breast treatment planning. Rather, the isocenter locations and beam directions are carefully selected, in that each beam is only responsible for a certain subvolume of the target, so as to minimize the volume of irradiated normal tissue. When put together, the beams will provide an adequate coverage of the target and expose only a minimal amount of normal tissue to radiation. We apply the nonisocentric planning technique to 2 previously treated clinical cases (breast and chest wall). Results: The proposed nonisocentric technique substantially improved sparing of the ipsilateral lung. Compared with conventional isocentric plans using 2 tangential beams, the mean lung dose was reduced by 38% and 50% using the proposed technique, and the volume of the ipsilateral lung receiving ≥20 Gy was reduced by a factor of approximately 2 and 3 for the breast and chest wall cases, respectively. The improvement in lung sparing is even greater compared with volumetric modulated arc therapy. Conclusions: A nonisocentric implementation of station parameter optimized radiation therapy has been proposed for breast radiation therapy. The new treatment scheme overcomes the limitations of existing approaches and affords a useful tool for conformal breast radiation therapy, especially in cases with extreme chest wall curvature.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2013

Clinical Implementation of Intrafraction Cone Beam Computed Tomography Imaging During Lung Tumor Stereotactic Ablative Radiation Therapy

Purpose: To develop and clinically evaluate a volumetric imaging technique for assessing intrafraction geometric and dosimetric accuracy of stereotactic ablative radiation therapy (SABR). Methods and Materials: Twenty patients received SABR for lung tumors using volumetric modulated arc therapy (VMAT). At the beginning of each fraction, pretreatment cone beam computed tomography (CBCT) was used to align the soft-tissue tumor position with that in the planning CT. Concurrent with dose delivery, we acquired fluoroscopic radiograph projections during VMAT using the Varian on-board imaging system. Those kilovolt projections acquired during millivolt beam-on were automatically extracted, and intrafraction CBCT images were reconstructed using the filtered backprojection technique. We determined the time-averaged target shift during VMAT by calculating the center of mass of the tumor target in the intrafraction CBCT relative to the planning CT. To estimate the dosimetric impact of the target shift during treatment, we recalculated the dose to the GTV after shifting the entire patient anatomy according to the time-averaged target shift determined earlier. Results: The mean target shift from intrafraction CBCT to planning CT was 1.6, 1.0, and 1.5 mm; the 95th percentile shift was 5.2, 3.1, 3.6 mm; and the maximum shift was 5.7, 3.6, and 4.9 mm along the anterior-posterior, left-right, and superior-inferior directions. Thus, the time-averaged intrafraction gross tumor volume (GTV) position was always within the planning target volume. We observed some degree of target blurring in the intrafraction CBCT, indicating imperfect breath-hold reproducibility or residual motion of the GTV during treatment. By our estimated dose recalculation, the GTV was consistently covered by the prescription dose (PD), that is, V100% above 0.97 for all patients, and minimum dose to GTV >100% PD for 18 patients and >95% PD for all patients. Conclusions: Intrafraction CBCT during VMAT can provide geometric and dosimetric verification of SABR valuable for quality assurance and potentially for treatment adaptation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2013

Automatic Prostate Tracking and Motion Assessment in Volumetric Modulated Arc Therapy With an Electronic Portal Imaging Device

Purpose: To assess the prostate intrafraction motion in volumetric modulated arc therapy treatments using cine megavoltage (MV) images acquired with an electronic portal imaging device (EPID). Methods and Materials: Ten prostate cancer patients were treated with volumetric modulated arc therapy using a Varian TrueBeam linear accelerator equipped with an EPID for acquiring cine MV images during treatment. Cine MV images acquisition was scheduled for single or multiple treatment fractions (between 1 and 8). A novel automatic fiducial detection algorithm that can handle irregular multileaf collimator apertures, field edges, fast leaf and gantry movement, and MV image noise and artifacts in patient anatomy was used. All sets of images (approximately 25,000 images in total) were analyzed to measure the positioning accuracy of implanted fiducial markers and assess the prostate movement. Results: Prostate motion can vary greatly in magnitude among different patients. Different motion patterns were identified, showing its unpredictability. The mean displacement and standard deviation of the intrafraction motion was generally less than 2.0 ± 2.0 mm in each of the spatial directions. In certain patients, however, the percentage of the treatment time in which the prostate is displaced more than 5 mm from its planned position in at least 1 spatial direction was 10% or more. The maximum prostate displacement observed was 13.3 mm. Conclusion: Prostate tracking and motion assessment was performed with MV imaging and an EPID. The amount of prostate motion observed suggests that patients will benefit from its real-time monitoring. Megavoltage imaging can provide the basis for real-time prostate tracking using conventional linear accelerators.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2012

4D cone beam CT via spatiotemporal tensor framelet

Purpose: On-board 4D cone beam CT (4DCBCT) offers respiratory phase-resolved volumetric imaging, and improves the accuracy of target localization in image guided radiation therapy. However, the clinical utility of this technique has been greatly impeded by its degraded image quality, prolonged imaging time, and increased imaging dose. The purpose of this letter is to develop a novel iterative 4DCBCT reconstruction method for improved image quality, increased imaging speed, and reduced imaging dose. Methods: The essence of this work is to introduce the spatiotemporal tensor framelet (STF), a high-dimensional tensor generalization of the 1D framelet for 4DCBCT, to effectively take into account of highly correlated and redundant features of the patient anatomy during respiration, in a multilevel fashion with multibasis sparsifying transform. The STF-based algorithm is implemented on a GPU platform for improved computational efficiency. To evaluate the method, 4DCBCT full-fan scans were acquired within 30 s, with a gantry rotation of 200°; STF is also compared with a state-of-art reconstruction method via spatiotemporal total variation regularization. Results: Both the simulation and experimental results demonstrate that STF-based reconstruction achieved superior image quality. The reconstruction of 20 respiratory phases took less than 10 min on an NVIDIA Tesla C2070 GPU card. The STF codes are available at https://sites.google.com/site/spatiotemporaltensorframelet . Conclusions: By effectively utilizing the spatiotemporal coherence of the patient anatomy among different respiratory phases in a multilevel fashion with multibasis sparsifying transform, the proposed STF method potentially enables fast and low-dose 4DCBCT with improved image quality.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2011

Bridging the gap between IMRT and VMAT: Dense angularly sampled and sparse intensity modulated radiation therapy

Purpose: To propose an alternative radiation therapy (RT) planning and delivery scheme with optimal angular beam sampling and intrabeam modulation for improved dose distribution while maintaining high delivery efficiency. Methods: In the proposed approach, coined as dense angularly sampled and sparse intensity modulated RT (DASSIM-RT), a large number of beam angles are used to increase the angular sampling, leading to potentially more conformal dose distributions as compared to conventional IMRT. At the same time, intensity modulation of the incident beams is simplified to eliminate the dispensable segments, compensating the increase in delivery time caused by the increased number of beams and facilitating the plan delivery. In a sense, the proposed approach shifts and transforms, in an optimal fashion, some of the beam segments in conventional IMRT to the added beams. For newly available digital accelerators, the DASSIM-RT delivery can be made very efficient by concatenating the beams so that they can be delivered sequentially without operator's intervention. Different from VMAT, the level of intensity modulation in DASSIS-RT is field specific and optimized to meet the need of each beam direction. Three clinical cases (a head and neck (HN) case, a pancreas case, and a lung case) are used to evaluate the proposed RT scheme. DASSIM-RT, VMAT, and conventional IMRT plans are compared quantitatively in terms of the conformality index (CI) and delivery efficiency. Results: Plan quality improves generally with the number and intensity modulation of the incident beams. For a fixed number of beams or fixed level of intensity modulation, the improvement saturates after the intensity modulation or number of beams reaches to a certain level. An interplay between the two variables is observed and the saturation point depends on the values of both variables. For all the cases studied here, the CI of DASSIM-RT with 15 beams and 5 intensity levels (0.90, 0.79, and 0.84 for the HN, pancreas, and lung cases, respectively) is similar with that of conventional IMRT with seven beams and ten intensity levels (0.88, 0.79, and 0.83) and is higher than that of single-arc VMAT (0.75, 0.75, and 0.82). It is also found that the DASSIM-RT plans generally have better sparing of organs-at-risk than IMRT plans. It is estimated that the dose delivery time of DASSIM-RT with 15 beams and 5 intensity levels is about 4.5, 4.4, and 4.2 min for the HN, pancreas, and lung case, respectively, similar to that of IMRT plans with 7 beams and 10 intensity levels.Conclusion: DASSIS-RT bridges the gap between IMRT and VMAT and allows optimal sampling of angular space and intrabeam modulation, thus it provides improved conformity in dose distributions while maintaining high delivery efficiency.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2011

A Bayesian approach to real-time 3D tumor localization via monoscopic x-ray imaging during treatment delivery

Purpose: Monoscopic x-ray imaging with on-board kV devices is an attractive approach for real-time image guidance in modern radiation therapy such as VMAT or IMRT, but it falls short in providing reliable information along the direction of imaging x-ray. By effectively taking consideration of projection data at prior times and/or angles through a Bayesian formalism, the authors develop an algorithm for real-time and full 3D tumor localization with a single x-ray imager during treatment delivery. Methods: First, a prior probability density function is constructed using the 2D tumor locations on the projection images acquired during patient setup. Whenever an x-ray image is acquired during the treatment delivery, the corresponding 2D tumor location on the imager is used to update the likelihood function. The unresolved third dimension is obtained by maximizing the posterior probability distribution. The algorithm can also be used in a retrospective fashion when all the projection images during the treatment delivery are used for 3D localization purposes. The algorithm does not involve complex optimization of any model parameter and therefore can be used in a ''plug-and-play'' fashion. The authors validated the algorithm using (1) simulated 3D linear and elliptic motion and (2) 3D tumor motion trajectories of a lung and a pancreas patient reproduced by a physical phantom. Continuous kV images were acquired over a full gantry rotation with the Varian TrueBeam on-board imaging system. Three scenarios were considered: fluoroscopic setup, cone beam CT setup, and retrospective analysis. Results: For the simulation study, the RMS 3D localization error is 1.2 and 2.4 mm for the linear and elliptic motions, respectively. For the phantom experiments, the 3D localization error is < 1 mm on average and < 1.5 mm at 95th percentile in the lung and pancreas cases for all three scenarios. The difference in 3D localization error for different scenarios is small and is not statistically significant. Conclusions: The proposed algorithm eliminates the need for any population based model parameters in monoscopic image guided radiotherapy and allows accurate and real-time 3D tumor localization on current standard LINACs with a single x-ray imager.

62 RADIOLOGY AND NUCLEAR MEDICINE↗