Image quality comparisons of coil setups in 3T MRI for brain and head and neck radiotherapy simulations.

PURPOSE: MRI is increasingly used for brain and head and neck radiotherapy treatment planning due to its superior soft tissue contrast. Flexible array coils can be arranged to encompass treatment immobilization devices, which do not fit in diagnostic head/neck coils. Selecting a flexible coil arrangement to replace a diagnostic coil should rely on image quality characteristics and patient comfort. We compared image quality obtained with a custom UltraFlexLarge18 (UFL18) coil setup against a commercial FlexLarge4 (FL4) coil arrangement, relative to a diagnostic Head/Neck20 (HN20) coil at 3T. METHODS: The large American College of Radiology (ACR) MRI phantom was scanned monthly in the UFL18, FL4, and HN20 coil setup over 2 years, using the ACR series and three clinical sequences. High-contrast spatial resolution (HCSR), image intensity uniformity (IIU), percent-signal ghosting (PSG), low-contrast object detectability (LCOD), signal-to-noise ratio (SNR), and geometric accuracy were calculated according to ACR recommendations for each series and coil arrangement. Five healthy volunteers were scanned with the clinical sequences in all three coil setups. SNR, contrast-to-noise ratio (CNR) and artifact size were extracted from regions-of-interest along the head for each sequence and coil setup. For both experiments, ratios of image quality parameters obtained with UFL18 or FL4 over those from HN20 were formed for each coil setup, grouping the ACR and clinical sequences. RESULTS: Wilcoxon rank-sum tests revealed significantly higher (p < 0.001) LCOD, IIU and SNR, and lower PSG ratios with UFL18 than FL4 on the phantom for the clinical sequences, with opposite PSG and SNR trends for the ACR series. Similar statistical tests on volunteer data corroborated that SNR ratios with UFL18 (0.58 ± 0.19) were significantly higher (p < 0.001) than with FL4 (0.51 ± 0.18) relative to HN20. CONCLUSIONS: The custom UFL18 coil setup was selected for clinical application in MR simulations due to the superior image quality demonstrated on a phantom and volunteers for clinical sequences and increased volunteer comfort.

Feasibility and clinical implementation of MRI-guided surface brachytherapy.

Purpose: Best practices for high-dose-rate surface applicator brachytherapy treatment (SABT) have long relied on computed tomography (CT)-based imaging to visualize diseased sites for treatment planning. Compared with magnetic resonance (MR)-based imaging, CT provides insufficient soft tissue contrast. This work described the feasibility of clinical implementation of MR-based imaging in SABT planning to provide individualized treatment optimization. Material and methods: A 3D-printed phantom was used to fit Freiberg flap-style (Elekta, The Netherlands) applicator. Images were taken using an optimized pointwise encoding time reduction with radial acquisition (PETRA) MR sequence for catheter visualization, and a helical CT scan to generate parallel treatment plans. This clinical study included three patients undergoing SABT for Dupuytren's contracture/palmar fascial fibromatosis imaged with the same modalities.SABT planning was performed in Oncentra Brachy (Elekta Brachytherapy, The Netherlands) treatment planning software. A geometric analysis was conducted by comparing CT-based digitization with MR-based digitization. CT and MR dwell positions underwent a rigid registration, and average Euclidean distances between dwell positions were calculated. A dosimetric comparison was performed, including point-based dose difference calculations and volumetric segmentations with Dice similarity coefficient (DSC) calculations. Results: Euclidean distances between dwell positions from CT-based and MR-based plans were on average 0.68 ±0.05 mm and 1.35 ±0.17 mm for the phantom and patients, respectively. The point dose difference calculations were on average 0.92% for the phantom and 1.98% for the patients. The D95 and D90 DSC calculations were both 97.9% for the phantom, and on average 93.6% and 94.2%, respectively, for the patients. Conclusions: The sub-millimeter accuracy of dwell positions and high DSC's (> 0.95) of the phantom demonstrated that digitization was clinically acceptable, and accurate treatment plans were produced using MR-only imaging. This novel approach, MRI-guided SABT, will lead to individualized prescriptions for potentially improved patient outcomes.

Edge roughness quantifies impact of physician variation on training and performance of deep learning auto-segmentation models for the esophagus.

Manual segmentation of tumors and organs-at-risk (OAR) in 3D imaging for radiation-therapy planning is time-consuming and subject to variation between different observers. Artificial intelligence (AI) can assist with segmentation, but challenges exist in ensuring high-quality segmentation, especially for small, variable structures, such as the esophagus. We investigated the effect of variation in segmentation quality and style of physicians for training deep-learning models for esophagus segmentation and proposed a new metric, edge roughness, for evaluating/quantifying slice-to-slice inconsistency. This study includes a real-world cohort of 394 patients who each received radiation therapy (mainly for lung cancer). Segmentation of the esophagus was performed by 8 physicians as part of routine clinical care. We evaluated manual segmentation by comparing the length and edge roughness of segmentations among physicians to analyze inconsistencies. We trained eight multiple- and individual-physician segmentation models in total, based on U-Net architectures and residual backbones. We used the volumetric Dice coefficient to measure the performance for each model. We proposed a metric, edge roughness, to quantify the shift of segmentation among adjacent slices by calculating the curvature of edges of the 2D sagittal- and coronal-view projections. The auto-segmentation model trained on multiple physicians (MD1-7) achieved the highest mean Dice of 73.7 ± 14.8%. The individual-physician model (MD7) with the highest edge roughness (mean ± SD: 0.106 ± 0.016) demonstrated significantly lower volumetric Dice for test cases compared with other individual models (MD7: 58.5 ± 15.8%, MD6: 67.1 ± 16.8%, p < 0.001). A multiple-physician model trained after removing the MD7 data resulted in fewer outliers (e.g., Dice ≤ 40%: 4 cases for MD1-6, 7 cases for MD1-7, Ntotal = 394). While we initially detected this pattern in a single clinician, we validated the edge roughness metric across the entire dataset. The model trained with the lowest-quantile edge roughness (MDER-Q1, Ntrain = 62) achieved significantly higher Dice (Ntest = 270) than the model trained with the highest-quantile ones (MDER-Q4, Ntrain = 62) (MDER-Q1: 67.8 ± 14.8%, MDER-Q4: 62.8 ± 15.7%, p < 0.001). This study demonstrates that there is significant variation in style and quality in manual segmentations in clinical care, and that training AI auto-segmentation algorithms from real-world, clinical datasets may result in unexpectedly under-performing algorithms with the inclusion of outliers. Importantly, this study provides a novel evaluation metric, edge roughness, to quantify physician variation in segmentation which will allow developers to filter clinical training data to optimize model performance.

Left Anterior Descending Coronary Artery Radiation Dose Association With All-Cause Mortality in NRG Oncology Trial RTOG 0617.

PURPOSE: A left anterior descending (LAD) coronary artery volume (V) receiving 15 Gy (V15 Gy) ≥10% has been recently observed to be an independent risk factor of major adverse cardiac events and all-cause mortality in patients with locally advanced non-small cell lung cancer treated with radiation therapy. However, this dose constraint has not been validated in independent or prospective data sets. METHODS AND MATERIALS: The NRG Oncology/Radiation Therapy Oncology Group (RTOG) 0617 data set from the National Clinical Trials Network was used. The LAD coronary artery was manually contoured. Multivariable Cox regression was performed, adjusting for known prognostic factors. Kaplan-Meier estimates of overall survival (OS) were calculated. For assessment of baseline cardiovascular risk, only age, sex, and smoking history were available. RESULTS: There were 449 patients with LAD dose-volume data and clinical outcomes available after 10 patients were excluded owing to unreliable LAD dose statistics. The median age was 64 years. The median LAD V15 Gy was 38% (interquartile range, 15%-62%), including 94 patients (21%) with LAD V15 Gy <10% and 355 (79%) with LAD V15 Gy ≥10%. Adjusting for prognostic factors, LAD V15 Gy ≥10% versus <10% was associated with an increased risk of all-cause mortality (hazard ratio [HR], 1.43; 95% confidence interval, 1.02-1.99; P = .037), whereas a mean heart dose ≥10 Gy versus <10 Gy was not (adjusted HR, 1.12; 95% confidence interval, 0.88-1.43; P = .36). The median OS for patients with LAD V15 Gy ≥10% versus <10% was 20.2 versus 25.1 months, respectively, with 2-year OS estimates of 47% versus 67% (P = .004), respectively. CONCLUSIONS: In a reanalysis of RTOG 0617, LAD V15 Gy ≥10% was associated with an increased risk of all-cause mortality. These findings underscore the need for improved cardiac risk stratification and aggressive risk mitigation strategies, including implementation of cardiac substructure dose constraints in national guidelines and clinical trials.

Observer preference of artificial intelligence-generated versus clinical prostate contours for ultrasound-based high dose rate brachytherapy.

BACKGROUND: For trans-rectal ultrasound (TRUS)-based high dose rate (HDR) prostate brachytherapy, prostate contouring can be challenging due to artifacts from implanted needles, bleeding, and calcifications. PURPOSE: To evaluate the geometric accuracy and observer preference of an artificial intelligence (AI) algorithm for generating prostate contours on TRUS images with implanted needles. METHODS: We conducted a retrospective study of 150 patients, who underwent HDR brachytherapy. These patients were randomly divided into training (104), validation (26) and testing (20) sets. An AI algorithm was trained/validated utilizing the TRUS image and reference (clinical) contours. The algorithm then provided contours for the test set. For evaluation, we calculated the Dice coefficient between AI and reference prostate contours. We then presented AI and reference contours to eight clinician observers, and asked observers to select their preference. Observers were blinded to the source of contours. We calculated the percentage of cases in which observers preferred AI contours. Lastly, we evaluate whether the presence of AI contours improved the geometric accuracy of prostate contours provided by five resident observers for a 10-patient subset. RESULTS: The median Dice coefficient between AI and reference contours was 0.92 (IQR: 0.90-0.94). Observers preferred AI contours for a median of 57.5% (IQR: 47.5, 65.0) of the test cases. For resident observers, the presence of AI contours was associated with a 0.107 (95% CI: 0.086, 0.128; p < 0.001) improvement in Dice coefficient for the 10-patient subset. CONCLUSION: The AI algorithm provided high-quality prostate contours on TRUS with implanted needles. Further prospective study is needed to better understand how to incorporate AI prostate contours into the TRUS-based HDR brachytherapy workflow.

Prospective deployment of an automated implementation solution for artificial intelligence translation to clinical radiation oncology.

Introduction: Artificial intelligence (AI)-based technologies embody countless solutions in radiation oncology, yet translation of AI-assisted software tools to actual clinical environments remains unrealized. We present the Deep Learning On-Demand Assistant (DL-ODA), a fully automated, end-to-end clinical platform that enables AI interventions for any disease site featuring an automated model-training pipeline, auto-segmentations, and QA reporting. Materials and methods: We developed, tested, and prospectively deployed the DL-ODA system at a large university affiliated hospital center. Medical professionals activate the DL-ODA via two pathways (1): On-Demand, used for immediate AI decision support for a patient-specific treatment plan, and (2) Ambient, in which QA is provided for all daily radiotherapy (RT) plans by comparing DL segmentations with manual delineations and calculating the dosimetric impact. To demonstrate the implementation of a new anatomy segmentation, we used the model-training pipeline to generate a breast segmentation model based on a large clinical dataset. Additionally, the contour QA functionality of existing models was assessed using a retrospective cohort of 3,399 lung and 885 spine RT cases. Ambient QA was performed for various disease sites including spine RT and heart for dosimetric sparing. Results: Successful training of the breast model was completed in less than a day and resulted in clinically viable whole breast contours. For the retrospective analysis, we evaluated manual-versus-AI similarity for the ten most common structures. The DL-ODA detected high similarities in heart, lung, liver, and kidney delineations but lower for esophagus, trachea, stomach, and small bowel due largely to incomplete manual contouring. The deployed Ambient QAs for heart and spine sites have prospectively processed over 2,500 cases and 230 cases over 9 months and 5 months, respectively, automatically alerting the RT personnel. Discussion: The DL-ODA capabilities in providing universal AI interventions were demonstrated for On-Demand contour QA, DL segmentations, and automated model training, and confirmed successful integration of the system into a large academic radiotherapy department. The novelty of deploying the DL-ODA as a multi-modal, fully automated end-to-end AI clinical implementation solution marks a significant step towards a generalizable framework that leverages AI to improve the efficiency and reliability of RT systems.

Clinical validation of deep learning algorithms for radiotherapy targeting of non-small-cell lung cancer: an observational study.

BACKGROUND: Artificial intelligence (AI) and deep learning have shown great potential in streamlining clinical tasks. However, most studies remain confined to in silico validation in small internal cohorts, without external validation or data on real-world clinical utility. We developed a strategy for the clinical validation of deep learning models for segmenting primary non-small-cell lung cancer (NSCLC) tumours and involved lymph nodes in CT images, which is a time-intensive step in radiation treatment planning, with large variability among experts. METHODS: In this observational study, CT images and segmentations were collected from eight internal and external sources from the USA, the Netherlands, Canada, and China, with patients from the Maastro and Harvard-RT1 datasets used for model discovery (segmented by a single expert). Validation consisted of interobserver and intraobserver benchmarking, primary validation, functional validation, and end-user testing on the following datasets: multi-delineation, Harvard-RT1, Harvard-RT2, RTOG-0617, NSCLC-radiogenomics, Lung-PET-CT-Dx, RIDER, and thorax phantom. Primary validation consisted of stepwise testing on increasingly external datasets using measures of overlap including volumetric dice (VD) and surface dice (SD). Functional validation explored dosimetric effect, model failure modes, test-retest stability, and accuracy. End-user testing with eight experts assessed automated segmentations in a simulated clinical setting. FINDINGS: We included 2208 patients imaged between 2001 and 2015, with 787 patients used for model discovery and 1421 for model validation, including 28 patients for end-user testing. Models showed an improvement over the interobserver benchmark (multi-delineation dataset; VD 0·91 [IQR 0·83-0·92], p=0·0062; SD 0·86 [0·71-0·91], p=0·0005), and were within the intraobserver benchmark. For primary validation, AI performance on internal Harvard-RT1 data (segmented by the same expert who segmented the discovery data) was VD 0·83 (IQR 0·76-0·88) and SD 0·79 (0·68-0·88), within the interobserver benchmark. Performance on internal Harvard-RT2 data segmented by other experts was VD 0·70 (0·56-0·80) and SD 0·50 (0·34-0·71). Performance on RTOG-0617 clinical trial data was VD 0·71 (0·60-0·81) and SD 0·47 (0·35-0·59), with similar results on diagnostic radiology datasets NSCLC-radiogenomics and Lung-PET-CT-Dx. Despite these geometric overlap results, models yielded target volumes with equivalent radiation dose coverage to those of experts. We also found non-significant differences between de novo expert and AI-assisted segmentations. AI assistance led to a 65% reduction in segmentation time (5·4 min; p<0·0001) and a 32% reduction in interobserver variability (SD; p=0·013). INTERPRETATION: We present a clinical validation strategy for AI models. We found that in silico geometric segmentation metrics might not correlate with clinical utility of the models. Experts' segmentation style and preference might affect model performance. FUNDING: US National Institutes of Health and EU European Research Council.

Knowledge-based inverse treatment planning for low-dose-rate prostate brachytherapy.

PURPOSE: Permanent low-dose-rate brachytherapy is a widely used treatment modality for managing prostate cancer. In such interventions, treatment planning can be a challenging task and requires experience and skills of the planner. We developed a novel knowledge-based (KB) optimization method based on previous treatment plans. The purpose of this method was to generate clinically acceptable plans that do not require extensive manual adjustments in clinical scenarios. METHODS: Objective functions used in current inverse planning methods are preferably based on spatial invariant dose objectives rather than spatial dose distributions. Therefore, they are prone to return suboptimal plans resulting in time consuming plan adjustments. To overcome this limitation, a KB approach is introduced. The KB model uses the dose distributions of previous clinical plans projected onto a standardized geometry. From those standardized distributions a template plan is generated. The treatment plans were optimized with an in-house developed planning system by solving a constraint inverse optimization problem that mimics the projected template dose plan constraint to DVH metrics. The method is benchmarked under an IRB-approved retrospective study by comparing optimization time, dosimetric performance, and clinical acceptability against current clinical practice. The quality of the KB model is evaluated with a Turing test. RESULTS: The KB model consists of five high-quality treatment plans. Those plans were selected by one of our experts and showed all desired dosimetric features. After generating the model treatment plans were created with one run of the optimizer for the remaining 20 patients. The optimization time including needle optimization ranged from 6 to 29 s. Based on a Wilcoxon signed rank test the new plans are dosimetrically equivalent to current clinical practice. The Turing test showed that the proposed method generates plans that are equivalent to current clinical practice and that the dose prediction drives the optimization to achieve high-quality treatment plans. CONCLUSIONS: This study demonstrated that the proposed KB model was able to capture user-specific features in isodose lines which can be used to generate acceptable treatment plans with a single run of the optimization engine in under a minute. This could potentially reduce the time in the operating room and the time a patient is under anesthesia.

Development and clinical implementation of semi-automated treatment planning including 3D printable applicator holders in complex skin brachytherapy.

PURPOSE: High-dose-rate brachytherapy (HDR-BT) is a treatment option for malignant skin diseases compared to external beam radiation therapy, HDR-BT provides improved target coverage, better organ sparing, and has comparable treatment times. This is especially true for large clinical targets with complex topologies. To standardize and improve the quality and efficacy of the treatments, a novel streamlined treatment approach in complex skin HDR-BT was developed and implemented. This approach consists of auto generated treatment plans and a 3D printable applicator holder (3D-AH). MATERIALS AND METHODS: The in-house developed planning system automatically segments computed tomography simulation images (a), optimizes a treatment plan (b), and generates a model of the 3D-AH (c). The 3D-AH is used as an immobilization device for the flexible Freiburg flap applicator used to deliver treatment. The developed, automated planning is compared against the standard clinical plan generation process for a flat 10 × 10 cm2 field, curved fields with radii of 4, 6, and 8 cm, and a representative clinical case. The quality of the 3D print is verified via an additional CT of the flap applicator latched into the holder, followed by an automated rigid registration with the original planning CT. Finally, the methodology is implemented and tested clinically under an IRB approval. RESULTS: All automatically generated plans were reviewed and accepted for clinical use. For the clinical workflow, the coverage achieved at a prescription depth for the flat 4, 6, and 8 cm applicator was (100.0 ± 4.9)%, (100.0 ± 4.9)%, (96.0 ± 0.3)%, and (98.4 ± 0.3)%, respectively. For auto planning, the coverage was (99.9 ± 0.3)%, (100.0 ± 0.2)%, (100.0 ± 0.3)%, and (100.1 ± 0.2)%. For the clinical test case, the D90 for the clinical workflow and auto planning was found to be 93.5% and 100.29% of the prescribed dose, respectively. Processing of the patient's CT to generate trajectories and positions as well as the 3D model of the applicator took <5 min. CONCLUSION: This workflow automates time intensive catheter digitizing and treatment planning. Compared to printing full applicators, the use of 3D-AH reduces the complexity of the 3D prints, the amount of the material to be used, the time of 3D printing, and amount of quality assurance required. The proposed methodology improves the overall treatment plan quality in complex HDR-BT and impact patient treatment outcomes potentially.

Combined stereotactic biopsy and stepping-source interstitial irradiation of glioblastoma multiforme.

Patients diagnosed with glioblastoma multiforme receiving stereotactic biopsy only either due to tumor localization or impaired clinical status face a devastating prognosis with very short survival times. One strategy to provide an initial cytoreductive and palliative therapy at the time of the stereotactic biopsy is interstitial irradiation through the pre-defined trajectory of the biopsy channel. We designed a novel treatment planning system and evaluated the treatment potential of a fixed-source and a stepping-source algorithm for interstitial radiosurgery on non-spherical glioblastoma in direct adjacency to risk structures. Using both setups, we show that radiation doses delivered to 100% of the gross tumor volume shifts from sub-therapeutic (10-12 Gy) to sterilizing single doses (25-30 Gy) when using the stepping source algorithm due to improved sparing of organs-at-risk. Specifically, the maximum doses at the brain stem were 100% of the PTV dose when a fixed central source and 38% when a stepping-source algorithm was used. We also demonstrated precision of intracranial target points and stability of superficial and deep trajectories using both a phantom and a body donor study. Our setup now for the first time provides a basis for a clinical proof-of-concept trial and may widen palliation options for patients with limited life expectancy that should not undergo time-consuming therapies.

A fast inverse treatment planning strategy facilitating optimized catheter selection in image-guided high-dose-rate interstitial gynecologic brachytherapy.

PURPOSE: Interstitial high-dose rate (HDR) brachytherapy is an important therapeutic strategy for the treatment of locally advanced gynecologic (GYN) cancers. The outcome of this therapy is determined by the quality of dose distribution achieved. This paper focuses on a novel yet simple heuristic for catheter selection for GYN HDR brachytherapy and their comparison against state of the art optimization strategies. The proposed technique is intended to act as a decision-supporting tool to select a favorable needle configuration. MATERIALS: The presented heuristic for catheter optimization is based on a shrinkage-type algorithm (SACO). It is compared against state of the art planning in a retrospective study of 20 patients who previously received image-guided interstitial HDR brachytherapy using a Syed Neblett template. From those plans, template orientation and position are estimated via a rigid registration of the template with the actual catheter trajectories. All potential straight trajectories intersecting the contoured clinical target volume (CTV) are considered for catheter optimization. Retrospectively generated plans and clinical plans are compared with respect to dosimetric performance and optimization time. RESULTS: All plans were generated with one single run of the optimizer lasting 0.6-97.4 s. Compared to manual optimization, SACO yields a statistically significant (P ≤ 0.05) improved target coverage while at the same time fulfilling all dosimetric constraints for organs at risk (OARs). Comparing inverse planning strategies, dosimetric evaluation for SACO and "hybrid inverse planning and optimization" (HIPO), as gold standard, shows no statistically significant difference (P > 0.05). However, SACO provides the potential to reduce the number of used catheters without compromising plan quality. CONCLUSION: The proposed heuristic for needle selection provides fast catheter selection with optimization times suited for intraoperative treatment planning. Compared to manual optimization, the proposed methodology results in fewer catheters without a clinically significant loss in plan quality. The proposed approach can be used as a decision support tool that guides the user to find the ideal number and configuration of catheters.

A fast multitarget inverse treatment planning strategy optimizing dosimetric measures for high-dose-rate (HDR) brachytherapy.

PURPOSE: In this study, we introduce a novel, fast, inverse treatment planning strategy for interstitial high-dose-rate (HDR) brachytherapy with multiple regions of interest solely based on dose-volume-histogram-related dosimetric measures (DMs). METHODS: We present a new problem formulation of the objective function that approximates the indicator variables of the standard DM optimization problem with a smooth logistic function. This problem is optimized by standard gradient-based methods. The proposed approach is then compared against state-of-the-art optimization strategies. RESULTS: All generated plans fulfilled prescribed DMs for all organs at risk. Compared to clinical practice, a statistically significant improvement (p=0.01) in coverage of target structures was achieved. Simultaneously, DMs representing high-dose regions were significantly reduced (p=0.01). The novel optimization strategies run-time was (0.8 ± 0.3) s and thus outperformed the best competing strategies of the state of the art. In addition, the novel DM-based approach was associated with a statistically significant (p=0.01) increase in the number of active dwell positions and a decrease in the maximum dwell time. CONCLUSIONS: The generated plans showed a clinically significant increase in target coverage with fewer hot spots, with an optimization time approximately three orders of magnitude shorter than manual optimization currently used in clinical practice. As optimization is solely based on DMs, intuitive, interactive, real-time treatment planning, which motivated the adoption of manual optimization in our clinic, is possible.

Feasibility of using single photon counting X-ray for lung tumor position estimation based on 4D-CT.

PURPOSE: In stereotactic body radiation therapy of lung tumors, reliable position estimation of the tumor is necessary in order to minimize normal tissue complication rate. While kV X-ray imaging is frequently used, continuous application during radiotherapy sessions is often not possible due to concerns about the additional dose. Thus, ultra low-dose (ULD) kV X-ray imaging based on a single photon counting detector is suggested. This paper addresses the lower limit of photons to locate the tumor reliably with an accuracy in the range of state-of-the-art methods, i.e. a few millimeters. METHOD: 18 patient cases with four dimensional CT (4D-CT), which serves as a-priori information, are included in the study. ULD cone beam projections are simulated from the 4D-CTs including Poisson noise. The projections from the breathing phases which correspond to different tumor positions are compared to the ULD projection by means of Poisson log-likelihood (PML) and correlation coefficient (CC), and template matching under these metrics. RESULTS: The results indicate that in full thorax imaging five photons per pixel suffice for a standard deviation in tumor positions of less than half a breathing phase. Around 50 photons per pixel are needed to achieve this accuracy with the field of view restricted to the tumor region. Compared to CC, PML tends to perform better for low photon counts and shifts in patient setup. Template matching only improves the position estimation in high photon counts. The quality of the reconstruction is independent of the projection angle. CONCLUSIONS: The accuracy of the proposed ULD single photon counting system is in the range of a few millimeters and therefore comparable to state-of-the-art tumor tracking methods. At the same time, a reduction in photons per pixel by three to four orders of magnitude relative to commercial systems with flatpanel detectors can be achieved. This enables continuous kV image-based position estimation during all fractions since the additional dose to the patient is negligible.