Integration of on-line imaging, plan adaptation and radiation delivery: proof of concept using digital tomosynthesis

Prospects for daily online adaptive radiotherapy via ethos for prostate cancer patients without nodal involvement using unedited CBCT auto-segmentation

Purpose

 Implementing new online adaptive radiation therapy technologies is challenging because extra clinical resources are required particularly expert contour review. Here, we provide the first assessment of Varian's Ethos™ adaptive platform for prostate cancer using no manual edits after auto‐segmentation to minimize this impact on clinical efficiency.

Methods

Twenty‐five prostate patients previously treated at our clinic were re‐planned using an Ethos™ emulator. Clinical target volumes (CTV) included intact prostate and proximal seminal vesicles. The following clinical margins were used: 3 mm posterior, 5 mm left/right/anterior, and 7 mm superior/inferior. Adapted plans were calculated for 10 fractions per patient using Ethos's auto‐segmentation and auto‐planning workflow without manual contouring edits. Doses and auto‐segmented structures were exported to our clinical treatment planning system where contours were modified as needed for all 250 CTVs and organs‐at‐risk. Dose metrics from adapted plans were compared to unadapted plans to evaluate CTV and OAR dose changes.

Results

Overall 96% of fractions required auto‐segmentation edits, although corrections were generally minor (<10% of the volume for 70% of CTVs, 88% of bladders, and 90% of rectums). However, for one patient the auto‐segmented CTV failed to include the superior portion of prostate that extended into the bladder at all 10 fractions resulting in under‐contouring of the CTV by 31.3% ± 6.7%. For the 24 patients with minor auto‐segmentation corrections, adaptation improved CTV D98% by 2.9% ± 5.3%. For non‐adapted fractions where bladder or rectum V90% exceeded clinical thresholds, adaptation reduced them by 13.1% ± 1.0% and 6.5% ± 7.3%, respectively.

Conclusion

 For most patients, Ethos's online adaptive radiation therapy workflow improved CTV D98% and reduced normal tissue dose when structures would otherwise exceed clinical thresholds, even without time‐consuming manual edits. However, for one in 25 patients, large contour edits were required and thus scrutiny of the daily auto‐segmentation is necessary and not all patients will be good candidates for adaptation.

Offline multiple adaptive planning strategy for concurrent irradiation of the prostate and pelvic lymph nodes

PURPOSE

Concurrent irradiation of the prostate and pelvic lymph nodes (PLNs) can be challenging due to the independent motion of the two target volumes. To address this challenge, the authors have proposed a strategy referred to as Multiple Adaptive Planning (MAP). To minimize the number of MAP plans, the authors' previous work only considered the prostate motion in one major direction. After analyzing the pattern of the prostate motion, the authors investigated a practical number of intensity-modulated radiotherapy (IMRT) plans needed to accommodate the prostate motion in two major directions simultaneously.

METHODS

Six patients, who received concurrent irradiation of the prostate and PLNs, were selected for this study. Nine MAP-IMRT plans were created for each patient with nine prostate contours that represented the prostate at nine locations with respect to the PLNs, including the original prostate contour and eight contours shifted either 5 mm in a single anterior-posterior (A-P), or superior-inferior (S-I) direction, or 5 mm in both A-P and S-I directions simultaneously. From archived megavoltage cone beam CT (MV-CBCT) and a dual imaging registration, 17 MV-CBCTs from 33 available MV-CBCT from these patients showed large prostate displacements (>3 mm in any direction) with respect to the pelvic bones. For each of these 17 fractions, one of nine MAP-IMRT plans was retrospectively selected and applied to the MV-CBCT for dose calculation. For comparison, a simulated isocenter-shifting plan and a reoptimized plan were also created for each of these 17 fractions. The doses to 95% (D95) of the prostate and PLNs, and the doses to 5% (D5) of the rectum and bladder were calculated and analyzed.

RESULTS

For the prostate, D95 > 97% of the prescription dose was observed in 16, 16, and 17 of 17 fractions for the MAP, isocenter-shifted, and reoptimized plans, respectively. For PLNs, D95 > 97% of the prescription doses was observed in 10, 3, and 17 of 17 fractions for the three types of verification plans, respectively. The D5 (mean ± SD) of the rectum was 45.78 ± 5.75, 45.44 ± 4.64, and 44.64 ± 2.71 Gy, and the D5 (mean ± SD) of the bladder was 45.18 ± 2.70, 46.91 ± 3.04, and 45.67 ± 3.61 Gy for three types of verification plans, respectively.

CONCLUSIONS

The MAP strategy with nine IMRT plans to accommodate the prostate motions in two major directions achieved good dose coverage to the prostate and PLNs. The MAP approach can be immediately used in clinical practice without requiring extra hardware and software.

Digital tomosynthesis (DTS) for verification of target position in early stage lung cancer patients

PURPOSE

The ability to verify intrafraction tumor position is clinically useful for hypofractionated treatments. Short arc kV digital tomosynthesis (DTS) could facilitate more frequent target verification. The authors used DTS combined with triangulation to determine the mean temporal position of small-volume lung tumor targets treated with stereotactic radiotherapy. DTS registration results were benchmarked against online clinical localization using registration between free-breathing cone-beam computed tomography (CBCT) and the average intensity projection (AvIP) of the planning 4DCT.

METHODS

In this retrospective study, 76 sets of kV-projection images from online CBCT scans of 13 patients were used to generate DTS image slices (CB-DTS) with nonclinical research software (DTS Toolkit, Varian Medical Systems). Three-dimensional tumor motion was 1.3-4 mm in six patients and 6.1-25.4 mm in seven patients on 4DCT (significant difference in the mean of the groups, P < 0.01). The 4DCT AvIP was used to digitally reconstruct the Reference-DTS. DTS registration and DTS registration combined with triangulation were investigated. Progressive shortening of total DTS arc lengths from 95° to 35° around 0° gantry position was evaluated for different scenarios: DTS registration using the entire arc; DTS registration plus triangulation using two nonoverlapping arcs; and for 55° and 45° total gantry rotation, DTS registration plus triangulation using two overlapping arcs. Finally, DTS registration plus triangulation performed at eight gantry angles, each separated by 45° was evaluated using full fan kV projection data for one patient with an immobile tumor and five patients with mobile tumors.

RESULTS

For DTS registration alone, shortening arc length did not influence accuracy in X- and Y-directions, but in Z-direction, mean deviations from online CBCT localization systematically increased for shorter arc length (P < 0.05). For example, using a 95° arc mean DTS-CBCT difference was 0.8 mm (1 SD = 0.6 mm) and for a 35° arc the mean was 2.4 mm (1 SD = 1.7 mm). DTS plus triangulation using nonoverlapping-arcs increased accuracy in Z-direction for tested arc lengths ≤55° (P < 0.01). Overlapping arcs increased accuracy in Y-direction for tumors with motion >4 mm (P < 0.02) but increased Z-direction accuracy was only observed with 55° total gantry rotation. The 95th percentile deviations with this overlapping technique in X-, Y-, and Z-directions were 1.3, 2.0, and 2.5 mm, respectively. For the five patients with mobile tumors where DTS + triangulation was performed with 45° intervals, the pooled deviation from online CBCT correction showed, for X-, Y-, and Z-directions, mean of 1.1 mm, standard deviations (SD) of 0.9, 1.0, and 0.9 mm, respectively. The mean + 2 SD was <3 mm for each direction.

CONCLUSIONS

Short-arc DTS verification of time averaged lung tumor position is feasible using free-breathing kV projection data and the AvIP of the 4DCT as a reference. Observed differences between DTS and online CBCT registration with AvIP were ≤3 mm (mean + 2 SD), however, the increased temporal resolution of DTS + triangulation also identified short period deviations from the average target position on the CBCT. Short-arc DTS appears promising for intrafraction tumor position monitoring during stereotactic lung radiotherapy delivered with a rotational technique.

Respiration-phase-matched digital tomosynthesis imaging for moving target verification: a feasibility study

PURPOSE

To develop a respiration-phase-matched digital tomosynthesis (DTS) technique to monitor moving targets, and to evaluate its accuracy for various imaging parameters and anatomical characteristics.

METHODS

Previously developed 3D-DTS techniques, registering onboard DTS (OB-DTS, reconstructed from onboard projections) to reference DTS (R-DTS, reconstructed from DRRs of 3D reference CT), are inadequate to monitor moving targets. The authors' proposed respiration-phase-matched DTS technique registers OB-DTS to R-DTS reconstructed from DRRs generated by the same phase images of 4D reference CT as the corresponding onboard projections. To evaluate the improved accuracy of the author's technique, the authors performed thoracic phantom studies using (1) simulation with the 4D digital extended-cardiac-torso (XCAT) phantom, and (2) experiments with an anthropomorphic motion phantom. The studies were performed for various: respiratory cycle (RC), scan angle, and fraction of RC contained therein. Also, the authors assessed the accuracy of their technique relative to target size/location, and respiration inconsistencies from the R-DTS to OB-DTS.

RESULTS

In both simulation and experimental studies, the respiration-phase-matched DTS technique is significantly more accurate in determining moving target positions. For 324 different scenarios simulated by XCAT, the respiration-phase-matched DTS technique localizes the 3D target position to errors of 1.07 ± 0.57 mm (mean ± S.D.), as compared to (a) 2.58 ± 1.37 and (b) 7.37 ± 4.18 mm, for 3D-DTS using 3D reference CT of (a) average intensity projection and (b) free-breathing CT. For 60 scenarios evaluated through experimental study, the uncertainties corresponding to those above are 1.24 ± 0.87, 2.42 ± 1.80, and 5.77 ± 6.45 mm, respectively. For a given scan angle, the accuracy of respiration-phase-matched DTS technique is less dependent on RC and the fraction of RC included in the scan. Increasing scan angle improves its accuracy. For different target locations, the targets near the chest wall or in the middle of lung provide higher registration accuracy compared to those near the mediastinum and diaphragm. Larger targets provide higher registration accuracy than small targets. Different respiratory cycle inconsistencies from R-DTS to OB-DTS minimally affect the registration accuracy. Increasing the respiratory amplitude inconsistencies will decrease the accuracy.

CONCLUSIONS

The respiration-phase-matched DTS is more accurate and robust in determining moving target positions than 3D-DTS. It has potential application in pretreatment setup, post-treatment analysis, and intrafractional target verification.

Incorporation of treatment plan spatial and temporal dose patterns into a prostate intrafractional motion management strategy

PURPOSE

Periodic MV∕KV radiographs taken during volumetric modulated arc therapy (VMAT) for hypofractionated treatment provide guidance in intrafractional motion management. The choice of imaging frequency and timing are key components in delivering the desired dose while reducing associated overhead such as imaging dose, preparation, and processing time. In this project the authors propose a paradigm with imaging timing and frequency based on the spatial and temporal dose patterns of the treatment plan.

METHODS

A number of control points are used in treatment planning to model VMAT delivery. For each control point, the sensitivity of individual target or organ-at-risk dose to motion can be calculated as the summation of dose degradations given the organ displacements along a number of possible motion directions. Instead of acquiring radiographs at uniform time intervals, MV∕KV image pairs are acquired indexed to motion sensitivity. Five prostate patients treated via hypofractionated VMAT are included in this study. Intrafractional prostate motion traces from the database of an electromagnetic tracking system are used to retrospectively simulate the VMAT delivery and motion management. During VMAT delivery simulation patient position is corrected based on the radiographic findings via couch movement if target deviation violates a patient-specific 3D threshold. The violation rate calculated as the percentage of traces failing the clinical dose objectives after motion correction is used to evaluate the efficacy of this approach.

RESULTS

Imaging indexed to a 10 s equitime interval and correcting patient position accordingly reduces the violation rate to 19.5% with intervention from 44.5% without intervention. Imaging indexed to the motion sensitivity further reduces the violation rate to 12.1% with the same number of images. To achieve the same 5% violation rate, the imaging incidence can be reduced by 40% by imaging indexed to motion sensitivity instead of time.

CONCLUSIONS

The simulation results suggest that image scheduling according to the characteristics of the treatment plan can improve the efficiency of intrafractional motion management. Using such a technique, the accuracy of delivered dose during image-guided hypofractionated VMAT treatment can be improved.

Optimal parameters for clinical implementation of breast cancer patient setup using Varian DTS software

Digital tomosynthesis (DTS) was evaluated as an alternative to cone‐beam computed tomography (CBCT) for patient setup. DTS is preferable when there are constraints with setup time, gantry‐couch clearance, and imaging dose using CBCT. This study characterizes DTS data acquisition and registration parameters for the setup of breast cancer patients using nonclinical Varian DTS software. DTS images were reconstructed from CBCT projections acquired on phantoms and patients with surgical clips in the target volume. A shift‐and‐add algorithm was used for DTS volume reconstructions, while automated cross‐correlation matches were performed within Varian DTS software. Triangulation on two short DTS arcs separated by various angular spread was done to improve 3D registration accuracy. Software performance was evaluated on two phantoms and ten breast cancer patients using the registration result as an accuracy measure; investigated parameters included arc lengths, arc orientations, angular separation between two arcs, reconstruction slice spacing, and number of arcs. The shifts determined from DTS‐to‐CT registration were compared to the shifts based on CBCT‐to‐CT registration. The difference between these shifts was used to evaluate the software accuracy. After findings were quantified, optimal parameters for the clinical use of DTS technique were determined. It was determined that at least two arcs were necessary for accurate 3D registration for patient setup. Registration accuracy of 2 mm was achieved when the reconstruction arc length was > 5° for clips with HU ≥ 1000°; larger arc length (≥ 8°) was required for very low HU clips. An optimal arc separation was found to be ≥ 20° and optimal arc length was 10°. Registration accuracy did not depend on DTS slice spacing. DTS image reconstruction took 10–30 seconds and registration took less than 20 seconds. The performance of Varian DTS software was found suitable for the accurate setup of breast cancer patients. Optimal data acquisition and registration parameters were determined.

PACS numbers: 87.57.‐s, 87.57.nf, 87.57.nj

Ultra-Fast Digital Tomosynthesis Reconstruction Using General-Purpose GPU Programming for Image-Guided Radiation Therapy

The purpose of this work is to demonstrate an ultra-fast reconstruction technique for digital tomosynthesis (DTS) imaging based on the algorithm proposed by Feldkamp, Davis, and Kress (FDK) using standard general-purpose graphics processing unit (GPGPU) programming interface. To this end, the FDK-based DTS algorithm was programmed “in-house” with C language with utilization of 1) GPU and 2) central processing unit (CPU) cards. The GPU card consisted of 480 processing cores (2 × 240 dual chip) with 1,242 MHz processing clock speed and 1,792 MB memory space. In terms of CPU hardware, we used 2.68 GHz clock speed, 12.0 GB DDR3 RAM, on a 64-bit OS. The performance of proposed algorithm was tested on twenty-five patient cases (5 lung, 5 liver, 10 prostate, and 5 head-and-neck) scanned either with a full-fan or half-fan mode on our cone-beam computed tomography (CBCT) system. For the full-fan scans, the projections from 157.5°–202.5° (45°-scan) were used to reconstruct coronal DTS slices, whereas for the half-fan scans, the projections from both 157.5°–202.5° and 337.5°–22.5° (2 × 45°-scan) were used to reconstruct larger FOV coronal DTS slices. For this study, we chose 45°-scan angle that contained ~80 projections for the full-fan and ~160 projections with 2 × 45°-scan angle for the half-fan mode, each with 1024 × 768 pixels with 32-bit precision. Absolute pixel value differences, profiles, and contrast-to-noise ratio (CNR) calculations were performed to compare and evaluate the images reconstructed using GPU- and CPU-based implementations. The time dependence on the reconstruction volume was also tested with (512 × 512) × 16, 32, 64, 128, and 256 slices. In the end, the GPU-based implementation achieved, at most, 1.3 and 2.5 seconds to complete full reconstruction of 512 × 512 × 256 volume, for the full-fan and half-fan modes, respectively. In turn, this meant that our implementation can process > 13 projections-per-second (pps) and > 18 pps for the full-fan and half-fan modes, respectively. Since commercial CBCT system nominally acquires 11 pps (with 1 gantry-revolution-per-minute), our GPU-based implementation is sufficient to handle the incoming projections data as they are acquired and reconstruct the entire volume immediately after completing the scan. In addition, on increasing the number of slices (hence volume) to be reconstructed from 16 to 256, only minimal increases in reconstruction time were observed for the GPU-based implementation where from 0.73 to 1.27 seconds and 1.42 to 2.47 seconds increase were observed for the full-fan and half-fan modes, respectively. This resulted in speed improvement of up to 87 times compared with the CPU-based implementation (for 256 slices case), with visually identical images and small pixel-value discrepancies (< 6.3%), and CNR differences (< 2.3%). With this achievement, we have shown that time allocation for DTS image reconstruction is virtually eliminated and that clinical implementation of this approach has become quite appealing. In addition, with the speed achievement, further image processing and real-time applications that was prohibited prior due to time restrictions can now be tempered with.

On-line adaptive radiation therapy: feasibility and clinical study

The purpose of this paper is to evaluate the feasibility and clinical dosimetric benefit of an on-line, that is, with the patient in the treatment position, Adaptive Radiation Therapy (ART) system for prostate cancer treatment based on daily cone-beam CT imaging and fast volumetric reoptimization of treatment plans. A fast intensity-modulated radiotherapy (IMRT) plan reoptimization algorithm is implemented and evaluated with clinical cases. The quality of these adapted plans is compared to the corresponding new plans generated by an experienced planner using a commercial treatment planning system and also evaluated by an in-house developed tool estimating achievable dose-volume histograms (DVHs) based on a database of existing treatment plans. In addition, a clinical implementation scheme for ART is designed and evaluated using clinical cases for its dosimetric qualities and efficiency.

A novel solid-angle tomosynthesis (SAT) scanning scheme

PURPOSE

Digital tomosynthesis (DTS) recently gained extensive research interests in both diagnostic and radiation therapy fields. Conventional DTS images are generated by scanning an x-ray source and flat-panel detector pair on opposite sides of an object, with the scanning trajectory on a one-dimensional curve. A novel tomosynthesis method named solid-angle tomosynthesis (SAT) is proposed, where the x-ray source scans on an arbitrary shaped two-dimensional surface.

METHODS

An iterative algorithm in the form of total variation regulated expectation maximization is developed for SAT image reconstruction. The feasibility and effectiveness of SAT is corroborated by computer simulation studies using three-dimensional (3D) numerical phantoms including a 3D Shepp-Logan phantom and a volumetric CT image set of a human breast.

RESULTS

SAT is able to cover more space in Fourier domain more uniformly than conventional DTS. Greater coverage and more isotropy in the frequency domain translate to fewer artifacts and more accurately restored features in the in-plane reconstruction.

CONCLUSIONS

Comparing with conventional DTS, SAT allows cone-shaped x-ray beams to project from more solid angles, thus provides more coverage in the spatial-frequency domain, resulting in better quality of reconstructed image.