Quantifying interfraction and intrafraction tumor motion in lung stereotactic body radiotherapy using respiration-correlated cone beam computed tomography

Quantitative assessment of intertarget position variations based on 4D-CT and 4D-CBCT simulations in single-isocenter multitarget lung stereotactic body radiation therapy

BACKGROUND

In single-isocenter multitarget stereotactic body radiotherapy (SBRT), geometric miss risks arise from uncertainties in intertarget position. However, its assessment is inadequate, and may be interfered by the reconstructed tumor position errors (RPEs) during simulated CT and cone beam CT (CBCT) acquisition. This study aimed to quantify intertarget position variations and assess factors influencing it.

METHODS

We analyzed data from 14 patients with 100 tumor pairs treated with single-isocenter SBRT. Intertarget position variation was measured using 4D-CT simulation to assess the intertarget position variations (ΔD) during routine treatment process. Additionally, a homologous 4D-CBCT simulation provided RPE-free comparison to determine the impact of RPEs, and isolating purely tumor motion induced ΔD to evaluate potential contributing factors.

RESULTS

The median ΔD was 4.3 mm (4D-CT) and 3.4 mm (4D-CBCT). Variations exceeding 5 mm and 10 mm were observed in 31.1% and 5.5% (4D-CT) and 20.4% and 3.4% (4D-CBCT) of fractions, respectively. RPEs necessitated an additional 1-2 mm safety margin. Intertarget distance and breathing amplitude variability showed weak correlations with variation (Rs = 0.33 and 0.31). The ΔD differed significantly by locations (upper vs. lower lobe and right vs. Left lung). Notably, left lung tumor pairs exhibited the highest risk.

CONCLUSIONS

This study provide a reliable way to assess intertarget position variation by using both 4D-CT and 4D-CBCT simulation. Consequently, single-isocenter SBRT for multiple lung tumors carries high risk of geometric miss. Tumor motion and RPE constitute a substantial portion of intertarget position variation, requiring correspondent strategies to minimize the intertarget uncertainties.

Markerless liver online adaptive stereotactic radiotherapy: feasibility analysisCervantes

Objective. Radio-opaque markers are recommended for image-guided radiotherapy in liver stereotactic ablative radiotherapy (SABR), but their implantation is invasive. We evaluate in thisin-silicostudy the feasibility of cone-beam computed tomography-guided stereotactic online-adaptive radiotherapy (CBCT-STAR) to propagate the target volumes without implanting radio-opaque markers and assess its consequence on the margin that should be used in that context.Approach. An emulator of a CBCT-STAR-dedicated treatment planning system was used to generate plans for 32 liver SABR patients. Three target volume propagation strategies were compared, analysing the volume difference between the GTVPropagatedand the GTVConventional, the vector lengths between their centres of mass (lCoM), and the 95th percentile of the Hausdorff distance between these two volumes (HD95). These propagation strategies were: (1) structure-guided deformable registration with deformable GTV propagation; (2) rigid registration with rigid GTV propagation; and (3) image-guided deformable registration with rigid GTV propagation. Adaptive margin calculation integrated propagation errors, while interfraction position errors were removed. Scheduled plans (PlanNon-adaptive) and daily-adapted plans (PlanAdaptive) were compared for each treatment fraction.Main results.The image-guided deformable registration with rigid GTV propagation was the best propagation strategy regarding tolCoM(mean: 4.3 +/- 2.1 mm), HD95 (mean 4.8 +/- 3.2 mm) and volume preservation between GTVPropagatedand GTVConventional. This resulted in a planning target volume (PTV) margin increase (+69.1% in volume on average). Online adaptation (PlanAdaptive) reduced the violation rate of the most important dose constraints ('priority 1 constraints', 4.2 versus 0.9%, respectively;p< 0.001) and even improved target volume coverage compared to non-adaptive plans (PlanNon-adaptive).Significance. Markerless CBCT-STAR for liver tumours is feasible using Image-guided deformable registration with rigid GTV propagation. Despite the cost in terms of PTV volumes, daily adaptation reduces constraints violation and restores target volumes coverage.

How much rotational error is clinically acceptable for single-isocenter/two-lesion lung SBRT treatment on halcyon ring delivery system (RDS)?

PURPOSE

SBRT treatment of two separate lung lesions via single-isocenter/multi-target (SIMT) plan on Halcyon RDS could improve patient comfort, compliance, patient throughput, and clinic efficiency. However, aligning two separate lung lesions synchronously via a single pre-treatment CBCT scan on Halcyon can be difficult due to rotational patient setup errors. Thus, to quantify the dosimetric impact, we simulated loss of target(s) coverage due to small, yet clinically observable rotational patient setup errors on Halcyon for SIMT treatments.

METHODS

Seventeen previously treated 4D-CT based SIMT lung SBRT patients with two separate lesions (total 34 lesions, 50 Gy in five fractions to each lesion) on TrueBeam (6MV-FFF) were re-planned on Halcyon (6MV-FFF) using a similar arc geometry (except couch rotation), dose engine (AcurosXB algorithm), and treatment planning objectives. Rotational patient setup errors of [± 0.5⁰ to ± 3.0⁰] on Halcyon were simulated via Velocity registration software in all three rotation axes and recalculated dose distributions in Eclipse treatment planning system. Dosimetric impact of rotational errors was evaluated for target coverage and organs at risk (OAR).

RESULTS

Average PTV volume and distance to isocenter were 23.7 cc and 6.1 cm. Average change in Paddick's conformity indexes were less than -5%, -10%, and -15% for 1°, 2°, and 3°, respectively for yaw, roll, and pitch rotation directions. Maximum drop off of PTV(D100%) coverage for 2° rotation was -2.0% (yaw), -2.2% (roll), and -2.5% (pitch). With ±1° rotational error, no PTV(D100%) loss was found. Due to anatomical complexity: irregular and highly variable tumor sizes and locations, highly heterogenous dose distribution, and steep dose gradient, no trend for loss of target(s) coverage as a function of distance to isocenter and PTV size was found. Change in maximum dose to OAR were acceptable per NRG-BR001 within ±1.0° rotation, but were up to 5 Gy higher to heart with 2° in the pitch rotation axis.

CONCLUSION

Our clinically realistic simulation results show that rotational patient setup errors up to 1.0° in any rotation axis could be acceptable for selected two separate lung lesions SBRT patients on Halcyon. Multivariable data analysis in large cohort is ongoing to fully characterize Halcyon RDS for synchronous SIMT lung SBRT.

Deep learning-based motion compensation for four-dimensional cone-beam computed tomography (4D-CBCT) reconstruction

BACKGROUND

Motion-compensated (MoCo) reconstruction shows great promise in improving four-dimensional cone-beam computed tomography (4D-CBCT) image quality. MoCo reconstruction for a 4D-CBCT could be more accurate using motion information at the CBCT imaging time than that obtained from previous 4D-CT scans. However, such data-driven approaches are hampered by the quality of initial 4D-CBCT images used for motion modeling.

PURPOSE

This study aims to develop a deep-learning method to generate high-quality motion models for MoCo reconstruction to improve the quality of final 4D-CBCT images.

METHODS

A 3D artifact-reduction convolutional neural network (CNN) was proposed to improve conventional phase-correlated Feldkamp-Davis-Kress (PCF) reconstructions by reducing undersampling-induced streaking artifacts while maintaining motion information. The CNN-generated artifact-mitigated 4D-CBCT images (CNN enhanced) were then used to build a motion model which was used by MoCo reconstruction (CNN+MoCo). The proposed procedure was evaluated using in-vivo patient datasets, an extended cardiac-torso (XCAT) phantom, and the public SPARE challenge datasets. The quality of reconstructed images for XCAT phantom and SPARE datasets was quantitatively assessed using root-mean-square-error (RMSE) and normalized cross-correlation (NCC).

RESULTS

The trained CNN effectively reduced the streaking artifacts of PCF CBCT images for all datasets. More detailed structures can be recovered using the proposed CNN+MoCo reconstruction procedure. XCAT phantom experiments showed that the accuracy of estimated motion model using CNN enhanced images was greatly improved over PCF. CNN+MoCo showed lower RMSE and higher NCC compared to PCF, CNN enhanced and conventional MoCo. For the SPARE datasets, the average (± standard deviation) RMSE in mm^-1 for body region of PCF, CNN enhanced, conventional MoCo and CNN+MoCo were 0.0040 ± 0.0009, 0.0029 ± 0.0002, 0.0024 ± 0.0003 and 0.0021 ± 0.0003. Corresponding NCC were 0.84 ± 0.05, 0.91 ± 0.05, 0.91 ± 0.05 and 0.93 ± 0.04.

CONCLUSIONS

CNN-based artifact reduction can substantially reduce the artifacts in the initial 4D-CBCT images. The improved images could be used to enhance the motion modeling and ultimately improve the quality of the final 4D-CBCT images reconstructed using MoCo.

An integrated ultrasound imaging and abdominal compression device for respiratory motion management in radiation therapy

BACKGROUND

Radiotherapy to tumors in the abdomen is challenging because of the significant organ movement and tissue deformation caused by respiration.

PURPOSE

A motion management strategy that integrated ultrasound (US) imaging with abdominal compression was developed and evaluated, where US was used to real-time monitor organ motion after abdominal compression.

METHODS

A device that combined a US imaging system and an abdominal compression plate (ACP) was developed. Twenty-one healthy volunteers were involved to evaluate the motion management efficacy. Each volunteer was immobilized on a flat bench by the device. Abdominal US data were successively collected with and without ACP compression and experiments were repeated three times to verify the imaging reproducibility. A template matching algorithm based on normalized cross correlation (NCC) was implemented to track the targets (vessels in the liver, pancreas and stomach) automatically. The matching algorithm was validated by comparing with the manual references. Automatic tracking was judged as failed if the center of mass difference from manual tracking was beyond a failure threshold. Based on the locations obtained through the template matching algorithm, the motion correlation between liver and pancreas/stomach was investigated using Pearson correlation test. Paired Student's t-test was used to analyze the difference between the results without and with ACP compression.

RESULTS

The liver motion amplitude over all 21 volunteers was significantly (p<0.001) reduced from 14.9 ± 5.5/3.4 ± 1.8 mm in superior-inferior (SI)/anterior-posterior (AP) direction before ACP compression to 7.3 ± 1.5/1.6 ± 0.7 mm after ACP compression. The mean liver motion standard deviation before compression was on average 2.8/1.4 mm in SI/AP direction, and was significantly (p<0.001) reduced to 0.9/0.4 mm after compression. The failure rates of automatic tracking for liver, pancreas and stomach were reduced for failure thresholds of 1-5 mm after applying ACP. The Pearson correlation coefficients between liver and pancreas/stomach were 0.98/0.97 without ACP and 0.96/0.94 with ACP in SI direction, and were 0.68/0.68 and 0.43/0.42 in AP direction. The motion prediction errors for pancreas/stomach with ACP have significantly (p<0.001) reduced to 0.45 ± 0.36/0.52 ± 0.43 mm from 0.69 ± 0.56/0.71 ± 0.66 mm without ACP in SI direction, and to 0.38 ± 0.33/0.39 ± 0.27 mm from 0.44 ± 0.35/0.61 ± 0.59 mm in AP direction.

CONCLUSIONS

The proposed strategy that combines real-time US imaging and abdominal compression has the potential to reduce the abdominal organ motion while improving both target tracking reliability and motion reproducibility. Furthermore, the observed correlation between liver and pancreas/stomach motion indicates the possibility of indirect pancreas/stomach tracking using liver markers as tracking surrogates. The strategy is expected to provide an alternative for respiratory motion management in the radiation treatment of abdominal tumors. This article is protected by copyright. All rights reserved.

Treatment Time Optimization in Single Fraction Stereotactic Ablative Radiation Therapy: A 10-Year Institutional Experience

Purpose

Stereotactic ablative radiation therapy (SABR) delivered in a single fraction (SF) can be considered to have higher uncertainty given that the error probability is concentrated in a single session. This study aims to report the variation in technology and technique used and its effect on intrafraction motion based on a 10 years of experience in SF SABR.

Methods and Materials

Records of patients receiving SF SABR delivered at our instruction between 2010 and 2019 were included. Treatment parameters were extracted from the patient management database by using an in-house script. Treatment time was defined as the time difference between the first image acquisition to the last beam off of a single session. The intrafraction variation was measured from the 3-dimensional couch displacement measured after the first cone beam computed tomography (CBCT) acquired during a treatment.

Results

The number of SF SABR increased continuously from 2010 to 2019 and were mainly lung treatments. Treatment time was minimized by using volumetric modulated arc therapy, flattening filter-free dose rate, and coplanar field (24 ± 9 min). Treatment time increased as the number of CBCTs per session increased. The most common scenario involved both 2 and 3 CBCTs per session. On the average, a CBCT acquisition added 6 minutes to the treatment time. All treatments considered, the average intrafraction variation was 1.7 ± 1.6 mm.

Conclusions

SF SABR usage increased with time in our institution. The intrafraction motion was acceptable and therefore a single fraction is an efficacious treatment option when considering SABR.

A standardized workflow for respiratory-gated motion management decision-making

Purpose

Motion management of tumors within the lung and abdomen is challenging because it requires balancing tissue sparing with accuracy of hitting the target, while considering treatment delivery efficiency. Physicists can play an important role in analyzing four‐dimensional computed tomography (4DCT) data to recommend the optimal respiratory gating parameters for a patient. The goal of this work was to develop a standardized procedure for making recommendations regarding gating parameters and planning margins for lung and gastrointestinal stereotactic body radiotherapy (SBRT) treatments. In doing so, we hoped to simplify decision‐making and analysis, and provide a tool for troubleshooting complex cases.

Methods

Factors that impact gating decisions and planning target volume (PTV) margins were identified. The gating options included gating on exhale with approximately a 50% duty cycle (Gate3070), exhale gating with a reduced duty cycle (Gate4060), and treating for most of respiration, excluding only extreme inhales and exhales (Gate100). A standard operating procedure was developed, as well as a physics consult document to communicate motion management recommendations to other members of the treatment team. This procedure was implemented clinically for 1 year and results are reported below.

Results

Identified factors that impact motion management included the magnitude of motion observed on 4DCT, the regularity of breathing and quality of 4DCT data, and ability to observe the target on fluoroscopy. These were collated into two decision tables—one specific to lung tumors and another for gastrointestinal tumors—such that a physicist could answer a series of questions to determine the optimal gating and PTV margin. The procedure was used clinically for 252 sites from 213 patients treated with respiratory‐gated SBRT and standardized practice across our 12‐member physics team.

Conclusion

Implementation of a standardized procedure for respiratory gating had a positive impact in our clinic, improving efficiency and ease of 4DCT analysis and standardizing gating decision‐making amongst physicists.

Simulation of dosimetric consequences of intrafraction variation of tumor drift in lung cancer stereotactic body radiotherapy

Objective

The purpose of this study was to investigate the target dose discrepancy caused by intrafraction variation during stereotactic body radiotherapy (SBRT) for lung cancer.

Methods

Intensity-modulated radiation therapy (IMRT) plans were designed based on average computed tomography (AVG CT) utilizing the planning target volume (PTV) surrounding the 65% and 85% prescription isodoses in both phantom and patient cases. Variation was simulated by shifting the nominal plan isocenter along six directions from 0.5 mm to 4.5 mm with a 1-mm step size to produce a series of perturbed plans. The dose discrepancy between the initial plan and the perturbed plans was calculated as the percentage of the initial plan. Dose indices, including ΔD₉₉ for internal target volume (ITV) and gross tumor volume (GTV), were adopted as endpoint samples. The mean dose discrepancy was calculated under the 3-dimensional space distribution.

Results

We found that motion can lead to serious dose degradation of the target and ITV in lung SBRT, especially during SBRT with PTV surrounding the lower isodose line. Lower isodose line may lead to larger dose discrepancy, while make steeper dose fall-off gradient. This phenomenon was compromised when 3-dimensional space distribution was considered.

Discussion

This result may provide a prospective reference for target dose degradation due to motion during lung SBRT treatment.

Tumour motion analysis from planning to end of treatment course for a large cohort of peripheral lung SBRT targets

BACKGROUND

The aim is to quantify and analyse tumour motion during a course of treatment for lung SBRT patients.

MATERIAL AND METHODS

Peak-to-peak motion of 483 tumours in 441 patients treated with peripheral lung SBRT at a single institution over a two year period was measured on planning CT and at all treatment fractions. Planning 4D-CT scans were analysed using our clinical workflow involving deformable propagation of the delineated target to all phases. Similarly, acquisition of the 4D-CBCT data followed the clinical workflow based on XVI 5.0 available on Elekta linacs. Differences and correlations of the peak-to-peak motion on the planning CT and at treatment were analysed.

RESULTS

On the planning CT, a total of 81.4% of the tumours had a peak-to-peak motion <10 mm, and 96.1% had <20 mm. The largest motion was observed in the CC direction, with largest amplitude for tumours located in the caudal posterior part of the lung. The difference in amplitude in CC between planning CT and first fraction had a mean and standard deviation of 0.3 mm and 3.5 mm, respectively, and the largest differences were observed in the caudal posterior part of the lung. Patients with a difference in tumour motion amplitude exceeding two standard deviations (>7 mm) at the first fraction were evaluated individually, and they all had poor 4DCT image quality. The difference between the first and second/third fractions had a mean and standard deviation of 0.4 mm/0.5 mm and 2.0 mm/1.9 mm.

CONCLUSION

Tumour motion at first treatment was similar to motion at planning, and motion at subsequent treatments was very similar to motion at first treatment. Large tumour motions are located towards the caudal posterior tumour locations. Patients with poor 4D-CT image quality should be closely followed at the first treatment to verify the motion.

Incorporation of tumor motion directionality in margin recipe: The directional MidP strategy

PURPOSE

Planning target volume (PTV) definition based on Mid-Position (Mid-P) strategy typically integrates breathing motion from tumor positions variances along the conventional axes of the DICOM coordinate system. Tumor motion directionality is thus neglected even though it is one of its stable characteristics in time. We therefore propose the directional MidP approach (MidP dir), which allows motion directionality to be incorporated into PTV margins. A second objective consists in assessing the ability of the proposed method to better take care of respiratory motion uncertainty.

METHODS

11 lung tumors from 10 patients with supra-centimetric motion were included. PTV were generated according to the MidP and MidP dir strategies starting from planning 4D CT.

RESULTS

PTV_{MidP dir} volume didn't differ from the PTV_MidP volume: 31351 mm³ IC95% [17242-45459] vs. 31003 mm³ IC95% [ 17347-44659], p = 0.477 respectively. PTV_{MidP dir} morphology was different and appeared more oblong along the main motion axis. The relative difference between 3D and 4D doses was on average 1.09%, p = 0.011 and 0.74%, p = 0.032 improved with directional MidP for D_99% and D_95%. D_2% was not significantly different between both approaches. The improvement in dosimetric coverage fluctuated substantially from one lesion to another and was all the more important as motion showed a large amplitude, some obliquity with respect to conventional axes and small hysteresis.

CONCLUSIONS

Directional MidP method allows tumor motion to be taken into account more tightly as a geometrical uncertainty without increasing the irradiation volume.

Novel Clinically Weight-Optimized Dynamic Conformal Arcs (WO-DCA) for Liver SBRT: A Comparison with Volumetric Modulated Arc Therapy (VMAT)

PURPOSE

To evaluate the feasibility of shortening the duration of liver stereotactic radiotherapy (SBRT) without jeopardizing dosimetry or conformity by utilizing weight-optimized dynamic conformal arcs (WO-DCA) as opposed to volumetric modulated arc therapy (VMAT) for tumors away from critical structures.

METHODS

Nineteen patients with liver metastasis were included, previously treated with 50 Gy in 4 fractions with VMAT technique using two partial coplanar arcs of 6 MV beams delivered in high-definition multi-leaf collimator (HD-MLC). Two coplanar partial WO-DCA were generated on Pinnacle treatment planning system (TPS) for each patient; and MLC aperture around the planning target volume (PTV) was automatically generated at different margins for both arcs and maintained dynamically around the target during arc rotation. Weight of the two arcs using optimization method was adjusted between the arcs to maximize tumor coverage and protect organs at risk (OAR) based on the RTOG-0438 protocol.

RESULTS

The WO-DCA plans successfully "agreed" with the standard VMAT for OAR (liver, spinal cord, stomach, duodenum, small bowel, and heart) and PTV (Dmean, D98%, D2%, CI, and GI), with superior mean quality assurance (QA) pass rate (97.06 vs 93.00 for VMAT; P < 0.001 and t = 8.87). Similarly, the WO-DCA technique additionally reduced the beam-on time (3.26 vs 4.43; P < 0.001) and monitor unit (1860 vs 2705 for VMAT; P < 0.001) values significantly.

CONCLUSION

The WO-DCA plans might minimize small-field dosimetry errors and defeat patient-specific VMAT QA requirements due to the omission of MLC beam modulation through the target volume. The WO-DCA plans may additionally enable faster treatment delivery times and lower OAR without sacrificing target doses in SBRT of liver tumors away from critical structures.

Lung Stereotactic Body Radiotherapy (SBRT) Using Spot-Scanning Proton Arc (SPArc) Therapy: A Feasibility Study

Purpose

We developed a 4D interplay effect model to quantitatively evaluate breathing-induced interplay effects and assess the feasibility of utilizing spot-scanning proton arc (SPArc) therapy for hypo-fractionated lung stereotactic body radiotherapy (SBRT). The model was then validated by retrospective application to clinical cases.

Materials and Methods

A digital lung 4DCT phantoms was used to mimic targets in diameter of 3cm with breathing motion amplitudes: 5, 10, 15, and 20 mm, respectively. Two planning groups based on robust optimization were generated: (1) Two-field Intensity Modulated Proton Therapy (IMPT) plans and (2) SPArc plans via a partial arc. 5,000 cGy relative biological effectiveness (RBE) was prescribed to the internal target volume (ITV) in five fractions. To quantitatively assess the breathing induced interplay effect, the 4D dynamic dose was calculated by synchronizing the breathing pattern with the simulated proton machine delivery sequence, including IMPT, Volumetric repainting (IMPT_volumetric), iso-layered repainting (IMPT_layer) and SPArc. Ten lung patients’ 4DCT previously treated with VMAT SBRT, were used to validate the digital lung tumor model. Normal tissue complicated probability (NTCP) of chestwall toxicity was calculated.

Result

Target dose were degraded as the tumor motion amplitude increased. The 4D interplay effect phantom model indicated that motion mitigation effectiveness using SPArc was about five times of IMPT_volumetric or IMPT_layer using maximum MU/spot as 0.5 MU at 20 mm motion amplitude. The retrospective study showed that SPArc has an advantage in normal tissue sparing. The probability of chestwall’s toxicity were significantly improved from 40.2 ± 29.0% (VMAT) (p = 0.01) and 16.3 ± 12.0% (IMPT) (p = 0.01) to 10.1 ± 5.4% (SPArc). SPArc could play a significant role in the interplay effect mitigation with breathing-induced motion more than 20 mm, where the target D99 of 4D dynamic dose for patient #10 was improved from 4,514 ± 138 cGy [RBE] (IMPT) vs. 4,755 ± 129 cGy [RBE] (SPArc) (p = 0.01).

Conclusion

SPArc effectively mitigated the interplay effect for proton lung SBRT compared to IMPT with repainting and was associated with normal tissue sparing. This technology may make delivery of proton SBRT more technically feasible and less complex with fewer concerns over underdosing the target compared to other proton therapy techniques.

Application of discrete cosine transform to assess the effect of tumor motion variations on the definition of ITV in lung and liver SBRT

PURPOSE

To use Discrete Cosine Transform to include tumor motion variations on ITV definition of SBRT patients.

METHODS

Data from 66 patients was collected. 2D planar fluoroscopy images (FI) were available for 54 patients. Daily CBCT projections (CBCTp) from 29 patients were employed to measure interfraction amplitude variability. Systematic amplitude variations were obtained from 17 patients with data from both FI and CBCTp. Tumor motion curves obtained from FI were characterized with a Cosine model (CM), based on cosine functions to the power of 2, 4 or 6, and DCT. Performance of both models was evaluated by means of R² coefficient and by comparing their results on Internal Target Volume (ITV) margins against those calculated from original tumor motion curves. Amplitude variations from CBCTp, as well as estimations of baseline shift variations were added to the DCT model to account for their effect on ITV margins.

RESULTS

DCT replicated tumor motion curves with a mean R² values for all patients of 0.86, 0.91 and 0.96 for the lateral (LAT), anterior-posterior (AP) and cranio-caudal (CC) directions respectively. CM yielded worst results, with R² values of 0.64, 0.61 and 0.74 in the three directions. Interfraction amplitude variation increased ITV margins by a 9%, while baseline shift variability implied a 40% and 80-100% increase for normalized values of baseline shift of 0.2 and 0.4 respectively.

CONCLUSIONS

Probability distribution functions of tumor positions can be successfully characterized with DCT. This permits to include tumor motion variablilities obtained from patient population into patient specific ITVs.

Feasibility of a TPS-integrated method to incorporate tumor motion in the margin recipe

BACKGROUND AND PURPOSE

There are several alternatives to the widespread ITV strategy in order to account for breathing-induced motion in PTV margins. The most sophisticated one includes the generation of a motion-compensated CT scan with the CTV placed in its average position - the mid-position approach (MidP). In such configuration, PTV margins integrate breathing as another random error. Despite overall irradiated volume reduction, such approach is barely used in clinical practice because of its dependence to deformable registration and its unavailability in commercial treatment planning systems. As an alternative, the mid-ventilation approach (MidV) selects the phase in the 4D-CT scan that is the closest to the MidP, with a residual error accounted for in the PTV margin. We propose a treatment planning system-integrated strategy, aiming at better approximating the MidP approach without its drawbacks: Hybrid MidV-MidP approach, i.e., the delineation on the MidV-CT and translation at the mid-position coordinates using treatment planning system built-in capabilities.

MATERIAL AND METHODS

Forty-five lung lesions treated with stereotactic radiotherapy were selected. PTV was defined using MidP, MidV, Hybrid MidV-MidP and ITV strategies. Margin definitions were adapted and resulting PTVs were compared.

RESULTS

Hybrid MidV-MidP showed similar target volume and location than the MidP and confirmed that margin-incorporated tumor motion strategies lead to significantly smaller PTVs than the ITV with mean reduction of 26 ± 7%.

CONCLUSION

We report on the successful implementation of a pseudo-MidP solution without its inherent drawbacks. It answers the need for TPS-embedded tumor motion range identification and related margin's component calculation.

A method using 4D dose accumulation to quantify the interplay effect in lung stereotactic body radiation therapy

The purpose of this study was to devise and evaluate a method to quantify the dosimetric uncertainty produced by the interplay between the movement of multileaf collimator and respiratory motion in lung stereotactic body radiation therapy. The method calculates the dose distribution for all control points from a dynamic treatment in all respiratory phases. The methodology includes some characteristics of a patient's irregular breathing patterns. It selects, for each control point, the phases with maximum and minimum mean dose over the tumor and their corresponding adjacent phases, whenever necessary. According to this selection, the dose matrices from each control point are summed up to obtain two dose distributions in each phase, which are accumulated in the reference phase subsequently by deformable image registration (DIR). D ₉₅ and [Formula: see text] were calculated over those accumulated dose distributions for Gross Tumor Volume (GTV), Planning Target Volume-based on Internal Target Volume approach-and Evaluation Target Volume (ETV), a novel concept that applies to 4D dose accumulation. With the ETV, DIR and interplay uncertainties are separated. The methodology also evaluated how variations in the breathing rate and field size affects the mean dose received by the GTV. The method was applied retrospectively in five patients treated with intensity modulated radiotherapy-minimum area defined by the leaves configuration at any control point was at least 4 cm². Uncertainties in tumor coverage were small (in most patients, changes on D ₉₅ and [Formula: see text] were below 2% for GTV and ETV) but significant over- and under-dosages near ETV, which can be accentuated by highly irregular breathing. Uncertainties in mean dose for GTV tended to decrease exponentially with increasing field size and were reduced by an increase of breathing rate. The implementation of this method would be helpful to assess treatment quality in patients with irregular breathing. Furthermore, it could be used to study interplay uncertainties when small field sizes are used.

Risk of target coverage loss for stereotactic body radiotherapy treatment of synchronous lung lesions via single-isocenter volumetric modulated arc therapy

Treating multiple lung lesions synchronously via single‐isocenter volumetric modulated arc therapy (VMAT) stereotactic body radiation therapy (SBRT) improves treatment efficiency and patient compliance. However, aligning multiple lung tumors accurately on single pretreatment cone beam CTs (CBCTs) can be problematic. Tumors misaligned could lead to target coverage loss. To quantify this potential target coverage loss due to small, clinically realistic setup errors, a novel simulation method was developed. This method was used on 26 previously treated patients with two metastatic lung lesions. Patients were treated with 4D CT‐based, highly conformal noncoplanar VMAT plans (clinical VMAT) with 6MV‐flattening filter free (FFF) beam using AcurosXB dose calculation algorithm with heterogeneity corrections. A single isocenter was placed approximately between the lesions to improve patient convenience and clinic workflow. Average isocenter to tumor distance was 5.9 cm. Prescription dose was 54 Gy/50 Gy in 3/5 fractions. For comparison, a plan summation (simulated VMAT) was executed utilizing randomly simulated, clinically relevant setup errors, obtained from pretreatment setup, per treatment fraction, in Eclipse treatment planning system for each of the six degrees of freedom within ± 5.0 mm and ± 2°. Simulations yielded average deviations of 27.4% (up to 72% loss) (P < 0.001) from planned target coverage when treating multiple lung lesions using a single‐isocenter plan. The largest deviations from planned coverage and desired biological effective dose (BED10, with α/β = 10 Gy) were seen for the smallest targets (<10 cc), some of which received < 100 Gy BED10. Patient misalignment resulted in substantial decrease in conformity and increase in the gradient index, violating major characteristics of SBRT. Statistically insignificant differences were seen for normal tissue dose. Although, clinical follow‐up of these patients is ongoing, the authors recommend an alternative treatment planning strategy to minimize the probability of a geometric miss when treating small lung lesions synchronously with single‐isocenter VMAT SBRT plans.

Cleaning the dose falloff in lung SBRT plan

PURPOSE

To investigate a planning technique that can possibly reduce low-to-intermediate dose spillage (measured by R50%, D2cm values) in lung SBRT plans.

MATERIALS AND METHODS

Dose falloff outside the target was studied retrospectively in 102 SBRT VMAT plans of lung tumor. Plans having R50% and/or D2cm higher than recommended tolerances in RTOG protocols 0813 and 0915 were replanned with new optimization constraints using novel shell structures and novel constraints. Violations in the RTOG R50% value can be rectified with a dose constraint to a novel shell structure ("OptiForR50"). The construction of structure OptiForR50% and the novel optimization criteria translate the RTOG goals for R50% into direct inputs for the optimizer. Violations in the D2cm can be rectified using constraints on a 0.5 cm thick shell structure with inner surface 2cm from the PTV surface. Wilcoxon signed-rank test was used to compare differences in dose conformity, volume of hot spots, R50%, D2cm of the target in addition to the OAR doses. A two-sided P-value of 0.05 was used to assess statistical significance.

RESULTS

Among 102 lung SBRT plans with PTV sizes ranging from 5 to 179 cc, 32 plans with violations in R50% or D2cm were reoptimized. The mean reduction in R50% (4.68 vs 3.89) and D2cm (56.49 vs 52.51) was statistically significant both having P < 0.01. Target conformity index, volume of 105% isodose contour outside PTV, normal lung V20, and mean dose to heart and aorta were significantly lowered with P < 0.05.

CONCLUSION

The novel planning methodology using multiple shells including the novel OptiForR50 shell with precisely calculated dimensions and optimizer constraints lead to significantly lower values of R50% and D2cm and lower dose spillage in lung SBRT plans. All plans were successfully brought into the zone of no RTOG violations.

High versus low dose irradiation for tumor immune reprogramming

Local administration of ionizing radiation to tumors can promote anticancer immune responses that lead to the abscopal regression of distant metastases, especially in patients receiving systemic immune-checkpoint inhibitors. Growing preclinical evidence indicates that high-dose irradiation administered locally to destroy malignant lesions, can promote the release of danger-associated molecular patterns that lead to the recruitment of immune cells, thus inducing a systemic response against tumor antigens that protects against local disease relapse and also mediates distant antineoplastic effects. An accumulating body of preclinical evidence supports also the implementation of low-dose irradiation to induce tumor immune reprogramming. Here, we provide the rationale for a clinical research agenda to refine future clinical practice based on innovative combinations of radiation-immunotherapy.

On the use of single‐isocenter VMAT plans for SBRT treatment of synchronous multiple lung lesions: Plan quality, treatment efficiency, and early clinical outcomes

Cone‐beam computed tomography (CT)‐guided volumetric‐modulated arc therapy (VMAT) plans for stereotactic body radiotherapy (SBRT) treatment of synchronous multiple lung lesions with a flattening filter‐free (FFF) beam is a safe and highly effective treatment option for oligometastases lung cancer patients. Fourteen patients with metastatic non–small‐cell lung cancer (NSCLC) lesions (two to five) received a single‐isocenter VMAT SBRT treatment in our clinic. Four‐dimensional (4D) CT‐based treatment plans were generated using advanced AcurosXB‐based dose calculation algorithm using heterogeneity corrections with a single isocenter placed between/among the lesions. Compared to 10X‐FFF and traditional flattened 6X (6X‐FF) beams, 6X‐FFF beam produced highly conformal radiosurgical dose distribution to each target volume, reduced dose to adjacent organs at risk (OAR), and significantly reduced the lung SBRT fraction duration to < 3.5 min/fraction for 54/50 Gy treatments in 3/5 fractions — significantly improving patient convenience and clinic workflow. Early follow‐up CT imaging (mean, 9 months) results show high local control rates (100%) with no acute lung or rib toxicity. Longer clinical follow up in a larger patient cohort is ongoing to further validate the outcomes of this treatment approach.

The effect of megavoltage field size on intrafraction cone-beam CT image quality

To investigate the effects of scatter from a megavoltage treatment beam on intrafraction cone beam CT (CBCT) image quality. The effects of treatment beam field size and phantom geometry were investigated as well as the clinical success of IFI. Intrafraction imaging (IFI) was performed on four phantoms with four different MV field sizes using a 6 MV FFF source. The image quality of the intrafraction CBCT images was compared to that of a baseline CBCT (i.e. with no treatment beam on) and quantified using noise and low contrast visibility. Increasing the kV tube current was explored as a possible method to reduce noise induced by the MV photon scatter in the intrafraction-CBCTs. The clinical success of all IFI patients over a 2 month period was reviewed. Intrafraction-CBCT image quality and low-contrast visibility deteriorated as MV field size increased. The extent of image degradation was found to depend on the mass of the phantom resulting in a more pronounced effect for a pelvic phantom than a thoracic phantom. While increasing the tube current could reduce the noise in the intrafraction-CBCT images, increasing the current by a factor of 4 failed to reach baseline image quality. Anatomy was found to be the primary indication of clinical IFI failure with all observed failures occurring during abdominal treatments. Image quality was found to decrease with increasing MV field size and decrease with increasing treatment anatomy mass. When considering intrafraction imaging clinically, the primary indicator of IFI failure is treatment anatomy. IFI can be used during chest treatments with high success rates but care must be taken for abdominal treatments and failures should be expected.

A novel and clinically useful dynamic conformal arc (DCA)-based VMAT planning technique for lung SBRT

Purpose

Volumetric modulated arc therapy (VMAT) is gaining popularity for stereotactic treatment of lung lesions for medically inoperable patients. Due to multiple beamlets in delivery of highly modulated VMAT plans, there are dose delivery uncertainties associated with small‐field dosimetry error and interplay effects with small lesions. We describe and compare a clinically useful dynamic conformal arc (DCA)‐based VMAT (d‐VMAT) technique for lung SBRT using flattening filter free (FFF) beams to minimize these effects.

Materials and Methods

Ten solitary early‐stage I‐II non‐small‐cell lung cancer (NSCLC) patients were treated with a single dose of 30 Gy using 3–6 non‐coplanar VMAT arcs (clinical VMAT) with 6X‐FFF beams in our clinic. These clinically treated plans were re‐optimized using a novel d‐VMAT planning technique. For comparison, d‐VMAT plans were recalculated using DCA with user‐controlled field aperture shape before VMAT optimization. Identical beam geometry, dose calculation algorithm, grid size, and planning objectives were used. The clinical VMAT and d‐VMAT plans were compared via RTOG‐0915 protocol compliances for conformity, gradient indices, and dose to organs at risk (OAR). Additionally, treatment delivery efficiency and accuracy were recorded.

Results

All plans met RTOG‐0915 requirements. Comparing with clinical VMAT, d‐VMAT plans gave similar target coverage with better target conformity, tighter radiosurgical dose distribution with lower gradient indices, and dose to OAR. Lower total number of monitor units and small beam modulation factor reduced beam‐on time by 1.75 min (P < 0.001), on average (maximum up to 2.52 min). Beam delivery accuracy was improved by 2%, on average (P < 0.05) and maximum up to 6% in some cases for d‐VMAT plans.

Conclusion

This simple d‐VMAT technique provided excellent plan quality, reduced intermediate dose‐spillage, and dose to OAR while providing faster treatment delivery by significantly reducing beam‐on time. This novel treatment planning approach will improve patient compliance along with potentially reducing intrafraction motion error. Moreover, with less MLC modulation through the target, d‐VMAT could potentially minimize small‐field dosimetry errors and MLC interplay effects. If available, d‐VMAT planning approach is recommended for future clinical lung SBRT plan optimization.

A simple, yet novel hybrid-dynamic conformal arc therapy planning via flattening filter-free beam for lung stereotactic body radiotherapy

Purpose

Due to multiple beamlets in the delivery of highly modulated volumetric arc therapy (VMAT) plans, dose delivery uncertainties associated with small‐field dosimetry and interplay effects can be concerns in the treatment of mobile lung lesions using a single‐dose of stereotactic body radiotherapy (SBRT). Herein, we describe and compare a simple, yet clinically useful, hybrid 3D‐dynamic conformal arc (h‐DCA) planning technique using flattening filter‐free (FFF) beams to minimize these effects.

Materials and Methods

Fifteen consecutive solitary early‐stage I‐II non‐small‐cell lung cancer (NSCLC) patients who underwent a single‐dose of 30 Gy using 3–6 non‐coplanar VMAT arcs with 6X‐FFF beams in our clinic. These patients’ plans were re‐planned using a non‐coplanar hybrid technique with 2–3 differentially‐weighted partial dynamic conformal arcs (DCA) plus 4–6 static beams. About 60–70% of the total beam weight was given to the DCA and the rest was distributed among the static beams to maximize the tumor coverage and spare the organs‐at‐risk (OAR). The clinical VMAT and h‐DCA plans were compared via RTOG‐0915 protocol for conformity and dose to OAR. Additionally, delivery efficiency, accuracy, and overall h‐DCA planning time were recorded.

Results

All plans met RTOG‐0915 requirements. Comparison with clinical VMAT plans h‐DAC gave better target coverage with a higher dose to the tumor and exhibited statistically insignificance differences in gradient index, D_2cm, gradient distance and OAR doses with the exception of maximal dose to skin (P = 0.015). For h‐DCA plans, higher values of tumor heterogeneity and tumor maximum, minimum and mean doses were observed and were 10%, 2.8, 1.0, and 2.0 Gy, on average, respectively, compared to the clinical VMAT plans. Average beam on time was reduced by a factor of 1.51. Overall treatment planning time for h‐DCA was about an hour.

Conclusion

Due to no beam modulation through the target, h‐DCA plans avoid small‐field dosimetry and MLC interplay effects and resulting in enhanced target coverage by improving tumor dose (characteristic of FFF‐beam). The h‐DCA simplifies treatment planning and beam on time significantly compared to clinical VMAT plans. Additionally, h‐DCA allows for the real time target verification and eliminates patient‐specific VMAT quality assurance; potentially offering cost‐effective, same or next day SBRT treatments. Moreover, this technique can be easily adopted to other disease sites and small clinics with less extensive physics or machine support.

Flattening filter free VMAT for a stereotactic, single-dose of 30 Gy to lung lesion in a 15-min treatment slot

Cone‐beam CT‐guided single dose of lung stereotactic body radiotherapy (SBRT) treatment with a flattening filter free (FFF) beam and volumetric modulated arc therapy (VMAT) is a safe and highly effective treatment modality for selective small lung lesions. Four‐dimensional (4D) CT‐based treatment plans were generated using advanced AcurosXB algorithm for heterogeneity corrections. 6X‐FFF beam produced highly conformal radiosurgical dose distribution to the target and reduced lung SBRT fraction duration to less than 10 min for a single dose of 30 Gy, significantly improving patient comfort and clinic workflow. Early follow‐up CT imaging results (mean, 8 months) show high local control rates (100%) with no acute lung or rib toxicity. Longer clinical follow‐up in a larger patient cohort managed in this fashion is underway to further validate this treatment approach.

A new respiratory monitor system for four-dimensional computed tomography by measuring the pressure change on the back of body

OBJECTIVE

A novel respiratory monitoring method based on the periodical pressure change on the patient's back was proposed and assessed by applying to four-dimensional CT (4DCT) scanning.

METHODS

A pressure-based respiratory monitoring system is developed and validated by comparing to real-time position management (RPM) system. The pressure change and the RPM signal are compared with phase differences and correlations calculated. The 4DCT images are reconstructed by these two signals. Internal and skin artifacts due to mismatch between CT slices and respiratory phases are evaluated.

RESULTS

The pressure and RPM signals shows strong consistency (R = 0.68±0.19 (1SD)). The time shift is 0.26 ± 0.51 (1SD) s and the difference of breath cycle is 0.02 ± 0.17 (1SD) s. The quality of 4DCT images reconstructed by two signals is similar. For both methods, the number of patients with artifacts is eight and the maximum magnitudes of artifacts are 20 mm (internal) and 10 mm (skin). The average magnitudes are 8.8 mm (pressure) and 8.2 mm (RPM) for internal artifacts, and 5.2 mm (pressure) and 4.6 mm (RPM) for skin artifacts. The mean square gray value difference shows no significant difference (p = 0.52).

CONCLUSION

The pressure signal provides qualified results for respiratory monitoring in 4DCT scanning, demonstrating its potential application for respiration monitoring in radiotherapy.

ADVANCES IN KNOWLEDGE

Pressure change on the back of body is a novel and promising method to monitor respiration in radiotherapy, which may improve treatment comfort and provide more information about respiration and body movement.

FFF-VMAT for SBRT of lung lesions: Improves dose coverage at tumor-lung interface compared to flattened beams

Purpose

To quantify the differences in dosimetry as a function of ipsilateral lung density and treatment delivery parameters for stereotactic, single dose of volumetric modulated arc therapy (VMAT) lung stereotactic body radiation therapy (SBRT) delivered with 6X flattening filter free (6X‐FFF) beams compared to traditional flattened 6X (6X‐FF) beams.

Materials/methods

Thirteen consecutive early stage I–II non‐small‐cell‐lung cancer (NSCLC) patients were treated with highly conformal noncoplanar VMAT SBRT plans (3–6 partial arcs) using 6X‐FFF beam and advanced Acuros‐based dose calculations to a prescription dose of 30 Gy in one fraction to the tumor margin. These clinical cases included relatively smaller tumor (island tumors) sizes (2.0–4.2 cm diameters) and varying average ipsilateral lung densities between 0.14 g/cc and 0.34 g/cc. Treatment plans were reoptimized with 6X‐FF beams for identical beam/arc geometries and planning objectives. For same target coverage, the organs‐at‐risk (OAR) dose metrics as a function of ipsilateral lung density were compared between 6X‐FFF and 6X‐FF plans. Moreover, monitor units (MU), beam modulation factor (MF) and beam‐on time (BOT) were evaluated.

Results

Both plans met the RTOG‐0915 protocol compliance. The ipsilateral lung density and the tumor location heavily influenced the treatment plans with 6X‐FFF and 6X‐FF beams, showing differences up to 12% for the gradient indices. For similar target coverage, 6X‐FFF beams showed better target conformity, lower intermediate dose‐spillage, and lower dose to the OAR. Additionally, BOT was reduced by a factor of 2.3 with 6X‐FFF beams compared to 6X‐FF beams.

Conclusion

While prescribing dose to the tumor periphery, 6X‐FFF VMAT plans for stereotactic single‐dose lung SBRT provided similar target coverage with better dose conformity, superior intermediate dose‐spillage (improved dose coverage at tumor interface), and improved OAR sparing compared to traditional 6X‐FF beams and significantly reduced treatment time. The ipsilateral lung density and tumor location considerably affected dose distributions requiring special attention for clinical SBRT plan optimization on a per‐patient basis. Clinical follow up of these patients for tumor local‐control rate and treatment‐related toxicities is in progress.

Imaging as Part of a Quality Assurance Program: Predictors of Interobserver Variability for Pretreatment Image Registration for Lung SBRT

Purpose:

To evaluate the magnitude of interobserver variability in pretreatment image registration for lung stereotactic body radiation therapy patients in aggregate and within 3 clinical subgroups and to determine methods to identify patients expected to demonstrate larger variability.

Methods and Materials:

Retrospective image registration was performed for the first and last treatment fraction for 10 lung stereotactic body radiation therapy patients by 16 individual observers (5 physicians, 6 physicists, and 5 therapists). Registration translation values were compared within and between subgroups overall and between the first and the last fractions. Four metrics were evaluated as possible predictors for large interobserver variability.

Results:

The mean 3-dimensional displacement vector for all patients over all comparisons was 2.4 ± 1.8 mm. Three patients had mean 3-dimensional vector differences >3 mm. This cohort of patients showed a significant interfraction difference in variance (P value = .01), increasing from first fraction to last. A significant difference in interobserver variability was observed between physicians and physicists (P value < .01) and therapists and physicists (P value < .01) but not between physicians and therapists (P value = .07). Three of the 4 quantities evaluated as potential predictive metrics showed statistical correlation with increased interobserver variation, including target excursion and local target/lung contrast.

Conclusion:

Variability in pretreatment image guidance represents an important treatment consideration, particularly for stereotactic body radiation therapy, which employs small margins and a small number of treatment fractions. As a result of the data presented here, we have initiated weekly “registration rounds” to familiarize all staff physicians with the target and normal anatomy for each stereotactic body radiation therapy patient and minimize interobserver variations in image registration prior to treatment. The metrics shown here are capable of identifying patients for which large interobserver variations would be anticipated. These metrics may be used in the future to develop thresholds for additional interventions to mitigate registration variations.

Reconstruction of a high-quality volumetric image and a respiratory motion model from patient CBCT projections

PURPOSE

To develop and evaluate a method of reconstructing a patient- and treatment day- specific volumetric image and motion model from free-breathing cone-beam projections and respiratory surrogate measurements. This Motion-Compensated Simultaneous Algebraic Reconstruction Technique (MC-SART) generates and uses a motion model derived directly from the cone-beam projections, without requiring prior motion measurements from 4DCT, and can compensate for both inter- and intrabin deformations. The motion model can be used to generate images at arbitrary breathing points, which can be used for estimating volumetric images during treatment delivery.

METHODS

The MC-SART was formulated using simultaneous image reconstruction and motion model estimation. For image reconstruction, projections were first binned according to external surrogate measurements. Projections in each bin were used to reconstruct a set of volumetric images using MC-SART. The motion model was estimated based on deformable image registration between the reconstructed bins, and least squares fitting to model parameters. The model was used to compensate for motion in both projection and backprojection operations in the subsequent image reconstruction iterations. These updated images were then used to update the motion model, and the two steps were alternated between. The final output is a volumetric reference image and a motion model that can be used to generate images at any other time point from surrogate measurements.

RESULTS

A retrospective patient dataset consisting of eight lung cancer patients was used to evaluate the method. The absolute intensity differences in the lung regions compared to ground truth were 50.8 ± 43.9 HU in peak exhale phases (reference) and 80.8 ± 74.0 in peak inhale phases (generated). The 50th percentile of voxel registration error of all voxels in the lung regions with >5 mm amplitude was 1.3 mm. The MC-SART was also applied to measured patient cone-beam projections acquired with a linac-mounted CBCT system. Results from this patient data demonstrate the feasibility of MC-SART and showed qualitative image quality improvements compared to other state-of-the-art algorithms.

CONCLUSION

We have developed a simultaneous image reconstruction and motion model estimation method that uses Cone-beam computed tomography (CBCT) projections and respiratory surrogate measurements to reconstruct a high-quality reference image and motion model of a patient in treatment position. The method provided superior performance in both HU accuracy and positional accuracy compared to other existing methods. The resultant reference image and motion model can be combined with respiratory surrogate measurements to generate volumetric images representing patient anatomy at arbitrary time points.

Evaluation of the use of abdominal compression of the lung in stereotactic radiation therapy

The goal of this retrospective study was to determine the benefit in using abdominal compression to reduce tumor motion for patients treated with lung stereotactic body radiotherapy. Forty-four lung lesions (n = 44) from 37 patients (N = 37) treated at the University of Toledo's Dana Cancer Center were assessed by determining the overall tumor displacement along with possible surrogates such as change in tidal volume and diaphragm displacement, with and without abdominal compression. Measurements of lung capacity were acquired from the 4DCT at maximum and minimum respiration in order to determine the tidal volume, with and without abdominal compression. Tumor centroid and diaphragm apex motion was then assessed in 3 dimensions from phase 0 to phase 50. This was measured in centimeters using the ruler method on MIM software, both with and without the compression belt. Change in overall tumor movement was 0.61 cm ± 0.09 cm with compression, and 0.60 cm ± 0.09 cm without the compression belt. Delta tumor motion was reduced in 5 cases, increased (made worse) in 6 cases, and did not clinically impact the remaining 33 cases. Average tidal volume with abdominal compression was 379.7 mL or 12.0% ± 0.724% of total lung volume while average tidal volume without abdominal compression was 337.7 mL or 10.5% ± 0.649% of total lung volume. Change in diaphragm position throughout the breathing cycle was 1.21 cm ± 0.10 cm with compression, and 1.28 ± 0.13 cm without the compression belt. These findings indicate that abdominal compression may not be an effective method in the reduction of respiratory motion, and can even negatively impact tumor motion by increasing its displacement. Compression decreased tumor motion in 5 out of the 44 cases studied. The 5 cases that benefitted tended to be lesions close to the diaphragm but these 5 corresponded to less than half of the inferior lesions, suggesting that even inferior lung lesions may not be prime candidates for abdominal compression.

Simulation of dosimetry impact of 4DCT uncertainty in 4D dose calculation for lung SBRT

Background

Due to the heterogeneity of patient’s individual respiratory motion pattern in lung stereotactic body radiotherapy (SBRT), treatment planning dose assessment using a traditional four-dimensional computed tomography (4DCT_traditional) images based on a uniform breathing curve may not represent the true treatment dose delivered to the patient. The purpose of this study was to evaluate the accumulated dose discrepancy between based on the 4DCT_traditional and true 4DCT (4DCT_true) that incorporated with the patient’s real entire breathing motion. The study also explored a novel 4D robust planning strategy to compensate for such heterogeneity respiratory motion uncertainties.

Methods

Simulated and measured patient specific breathing curves were used to generate 4D targets motion CT images. Volumetric-modulated arc therapy (VMAT) was planned using two arcs. Accumulated dose was obtained by recalculating the plan dose on each individual phase image and then deformed the dose from each phase image to the reference image. The “4 D dose” (D^4D) and “true dose” (D^true) were the accumulated dose based on the 4DCT_traditional and 4DCT_true respectively. The average worse case dose discrepancy (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{\Delta D} $$\end{document}ΔD¯) between D^4D and D^true in all treatment fraction was calculated to evaluate dosimetric /planning parameters and correlate them with the heterogeneity of respiratory-induced motion patterns. A novel 4D robust optimization strategy for VMAT (4D Ro-VMAT) based on the probability density function(pdf) of breathing curve was proposed to improve the target coverage in the presence of heterogeneity respiratory motion. The data were assessed with a paired t-tests.

Results

With increasing breathing amplitude from 5 to 20 mm, target \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{\Delta {D}_{99}} $$\end{document}ΔD99¯, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{\Delta {D}_{95}} $$\end{document}ΔD95¯ increased from 1.59,1.39 to 10.15%,8.66% respectively. When the standard deviation of breathing amplitude increased from 15 to 35% of the mean amplitude, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{\Delta {D}_{99}} $$\end{document}ΔD99¯, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{\Delta {D}_{95}} $$\end{document}ΔD95¯ increased from 4.06,3.48 to 10.25%,6.63% respectively. The 4D Ro-VMAT plan significantly improve the target dose compared to VMAT plan.

Conclusion

When the breathing curve amplitude is more than 10 mm and standard deviation of amplitude is higher than 25% of mean amplitude, special care is needed to choose an appropriated dose accumulation approach to evaluate lung SBRT plan target coverage robustness. The proposed 4D Ro_VMAT strategy based on the pdf of patient specific breathing curve could effectively compensate such uncertainties.

Evaluation of plan quality and treatment efficiency for single-isocenter/two-lesion lung stereotactic body radiation therapy

Purpose/objectives

To evaluate the plan quality and treatment delivery efficiency of single‐isocenter/two‐lesions volumetric modulated arc therapy (VMAT) lung stereotactic body radiation therapy (SBRT).

Materials/methods

Eight consecutive patients with two peripherally located early stage nonsmall‐cell‐lung cancer (NSCLC) lung lesions underwent single‐isocenter highly conformal noncoplanar VMAT SBRT treatment in our institution. A single‐isocenter was placed between the two lesions. Doses were 54 or 50 Gy in 3 and 5 fractions respectively. Patients were treated every other day. Plans were calculated in Eclipse with AcurosXB algorithm and normalized to at least 95% of the planning target volume (PTV) receiving 100% of the prescribed dose. For comparison, two‐isocenter plans (isocenter placed centrally in each target) were retrospectively created. Conformity indices (CIs), heterogeneity index (HI), gradient index (GI), gradient distance (GD), and D_2cm were calculated. The normal lung V5, V10, V20, mean lung dose (MLD) and other organs at risk (OARs) doses were evaluated. Total number of monitor units (MUs), beam‐on time, and patient‐specific quality assurance (QA) results were recorded.

Results

The mean isocenter to tumor distance was 6.7 ± 2.3 cm. The mean combined PTV was 44.0 ± 23.4 cc. There was no clinically significant difference in CI, HI, GD, GI, D_2cm, and V20 including most of the OARs between single‐isocenter and two‐isocenter lung SBRT plans, evaluated per RTOG guidelines. However, for single‐isocenter plans as the distance between the lesions increased, the V5, V10, and MLD increased, marginally. The total number of MUs and beam‐on time was reduced by a factor of 1.5 for a single‐isocenter plan compared to a two‐isocenter plan. The single‐isocenter/two‐lesions VMAT lung SBRT QA plans demonstrated an accurate dose delivery of 98.1 ± 3.2% for clinical gamma passing rate of 3%/3 mm.

Conclusion

The SBRT treatment of two peripherally located lung lesions with a centrally placed single‐isocenter was dosimetrically equivalent to two‐isocenter plans. Faster treatment delivery for single‐isocenter treatment can improve patient compliance and reduce the amount of intrafraction motion errors for well‐suited patients.

A method for optimizing planning target volume margins for patients receiving lung stereotactic body radiotherapy

Lung stereotactic-body radiotherapy (SBRT) places additional requirements on targeting accuracy over standard approaches. In treatment planning, a tumour volume is geometrically expanded and the resulting planning target volume (PTV) is covered with the prescribed dose. This ensures full dose delivery despite various uncertainties encountered during treatment. We developed a retrospective technique for optimizing the PTV expansion for a patient population. The method relies on deformable image registration (DIR) of the planning CT to a treatment cone-beam CT (CBCT). The resulting transformation is used to map the planned target onto the treatment geometry, allowing the computation of the achieved target/PTV overlap. Basic validation of the method was performed using an anthropomorphic respiratory motion phantom. A self-validation technique was also implemented to allow estimation of the DIR error for the data being analyzed. Our workflow was used to retrospectively optimize PTV margin for 25 patients treated over 93 fractions. Targets for these patients were contoured on 4D CT images. SBRT delivery followed CBCT acquisition and a couch correction. A post-treatment CBCT was also acquired in some cases. Our basic validation demonstrated that the DIR-based technique is capable of transforming target volumes from planning CTs to treatment CBCTs with sub-mm accuracy. Our clinical analysis showed that the minimum percentages of target volumes covered for 3, 4, and 5 mm PTV margins were 92.1, 97.6, and 99.2, respectively. Analyzing data acquired before and just after treatment demonstrated that margins exceeding 5 mm did not significantly improve coverage. Finally, a 5 mm PTV margin achieved  ⩾95% target volume coverage with  ⩾95% probability. Our technique is accurate, automated, self-validating, and incorporates complex ITV shapes/deformations to allow PTV margin optimization. The analysis of clinical data indicates a 5 mm PTV margin is optimal for our process. This approach is generalizable to other disease sites and treatment strategies.

Data-driven respiratory motion compensation for four-dimensional cone-beam computed tomography (4D-CBCT) using groupwise deformable registration

PURPOSE

To demonstrate the feasibility of using a purely data-driven, a posteriori respiratory motion modeling and reconstruction compensation method to improve 4D-CBCT image quality under clinically relevant image acquisition conditions.

METHODS

Evaluated workflows that utilized a combination of groupwise deformable image registration and motion-compensated image reconstruction algorithms. Groupwise registration is an approach that simultaneously registers all temporal frames of a 4D image to a common reference instead of one at a time so as to minimize the influence of any individual time point on the global smoothness or accuracy of the resulting deformation model. Four-dimensional cone-beam CT (4D-CBCT) Feldkamp-Davis-Kress (FDK) reconstructions were registered to either iteratively computed mean respiratory phase (mean-frame) or preselected respiratory phase (fixed-frame) reference images to model respiratory motion. The resulting 4D transformations were used to deform projection data during the FDK backprojection operation to create motion-compensated reconstructions. Tissue interface sharpness (TIS) was defined as the slope of a sigmoid curve fit to a mobile tissue boundary and was used to evaluate image quality in regions susceptible to motion artifacts. Image quality improvement was assessed for 19 clinical cases by evaluating mitigation of view aliasing artifacts, TIS, image noise reduction, and contrast for implanted fiducial markers.

RESULTS

Average (standard deviation) diaphragm TIS recovery relative to initial 4D-CBCT reconstructions was observed to be 87% (46%) using fixed-frame registration alone; 87% (47%) using fixed frame with motion-compensated reconstruction; 101% (68%) using mean-frame registration alone; and 99% (65%) using mean frame with motion-compensated reconstruction. Noise was reduced in sampled soft tissue ROIs by 58% for both fixed-frame registration and registration with motion compensation and by 57% and 58% on average for the corresponding mean-frame methods, respectively. Average improvement in local CNR was observed to be respectively 93% and 98% for fixed-frame registration and registration with motion compensation methods and 116% and 111% for the corresponding mean-frame methods.

CONCLUSION

Data-driven groupwise registration and motion-compensated reconstruction offer a feasible means of improving the quality of 4D-CBCT images acquired under clinical conditions. The addition of motion compensation reconstruction after groupwise registration visibly reduced the impact of view aliasing artifacts for the clinical image datasets studied.

The dosimetric differences in calculating lung SBRT plans on different image data sets: Comparison of the free breathing scan to both the average intensity projection scan and to the sum of calculations on each respiratory phase of the 4DCT scan

The purpose of this study is to evaluate the dosimetric differences to the lung internal target volume (ITV) in stereotactic body radiation therapy (SBRT) when calculated on the free breathing (FB) scan in comparison to calculations on the average intensity projection (AIP) scan as well as the sum of dose calculated on each of the treated respiratory phases. The clinical treatment plan data for 16 SBRT lung patients were retrospectively chosen for this study, 5 of which were 30% to 70% respiratory phase gated. All patients had ITV contours and fixed monitor units from the treatment approved plan copied over to each scan on which calculations were to be made. The results of this study yielded 6 out of 16 patients with greater than or equal to 2% difference in ITV maximum dose (D₀) and 2 of 16 patients with greater than or equal to 2% difference in ITV minimum dose (D₁₀₀). The range of ITV dose differences for these 8 patients was 2% to 4.7% with no patients exceeding a 5% difference in D₀ or D₁₀₀. None of the patients had greater than or equal to 2% difference in ITV mean dose (D_mean). Sixty-three percent of patients with greater than 2% ITV dose difference in any of the categories reviewed were those patients with greater than 1 cm gross tumor motion. This study concluded that in order to reduce uncertainty in dose to the ITV, tumor motion should be assessed at simulation and limited to below 1 cm in any direction if possible.

Advances in the use of motion management and image guidance in radiation therapy treatment for lung cancer

The development of advanced radiation technologies, including intensity-modulated radiation therapy (IMRT), stereotactic body radiation therapy (SBRT) and proton therapy, has resulted in increasingly conformal radiation treatments. Recent evidence for the importance of minimizing dose to normal critical structures including the heart and lungs has led to incorporation of these advanced treatment modalities into radiation therapy (RT) for non-small cell lung cancer (NSCLC). While such technologies have allowed for improved dose delivery, implementation requires improved target accuracy with treatments, placing increasing importance on evaluating tumor motion at the time of planning and verifying tumor position at the time of treatment. In this review article, we describe issues and updates related both to motion management and image guidance in the treatment of NSCLC.

Evaluation of the 4D RADPOS dosimetry system for dose and position quality assurance of CyberKnife

PURPOSE

The Synchrony respiratory motion tracking of the CyberKnife system purports to provide real-time tumor motion compensation during robotic radiosurgery. Such a complex delivery system requires thorough quality assurance. In this work, RADPOS applicability as a dose and position quality assurance tool for CyberKnife treatments is assessed quantitatively for different phantom types and breathing motions, which increase in complexity to more closely resemble clinical situations.

METHODS

Two radiotherapy treatment experiments were performed where dose and position were measured with the RADPOS probe housed within a Solid Water phantom. For the first experiment, a Solid Water breast phantom was irradiated using isocentric beam delivery while stationary or moving sinusoidally in the anterior/posterior direction. For the second experiment, a phantom consisting of a Solid Water tumor in lung equivalent material was irradiated using isocentric and non-isocentric beam delivery while either stationary or moving. The phantom movement was either sinusoidal or based on a real patient's breathing waveform. For each experiment, RADPOS dose measurements were compared to EBT3 GafChromic film dose measurements and the CyberKnife treatment planning system's (TPS) Monte Carlo and ray-tracing dose calculation algorithms. RADPOS position measurements were compared to measurements made by the CyberKnife system and to the predicted breathing motion models used by the Synchrony respiratory motion compensation.

RESULTS

For the static and dynamic (i.e., sinusoidal motion) cases of the breast experiment, RADPOS, film and the TPS agreed at the 2.0% level within 1.1 σ of estimated combined uncertainties. RADPOS position measurements were in good agreement with LED and fiducial position measurements, where the average standard deviation (SD) of the differences between any two of the three position datasets was ≤0.5 mm for all directions. For the 10 mm peak to peak amplitude sinusoidal motion of the breast experiment, the average Synchrony correlation errors were ≤0.2 mm, indicative of an accurate predictive model. For all the cases of the lung experiment, RADPOS and film measurements agreed with each other at the 2.0% level within 1.5 σ of estimated experimental uncertainties provided that the measurements were corrected for imaging dose. The measured dose for RADPOS and film were 4.0% and 3.4% higher, respectively, than the TPS for the most complex dynamic cases (i.e., irregular motion) considered for the lung experiment. Assessment of the Synchrony correlation models by RADPOS showed that model accuracy declined as motion complexity increased; the SD of the differences between RADPOS and model position data measurements was ≤0.8 mm for sinusoidal motion but increased to ≤2.6 mm for irregular patient waveform motion. These results agreed with the Synchrony correlation errors reported by the CyberKnife system.

CONCLUSIONS

RADPOS is an accurate and precise QA tool for dose and position measurements for CyberKnife deliveries with respiratory motion compensation.

The impact of technology on the changing practice of lung SBRT

Stereotactic body radiotherapy (SBRT) for lung tumours has been gaining wide acceptance in lung cancer. Here, we review the technological evolution of SBRT delivery in lung cancer, from the first treatments using the stereotactic body frame in the 1990's to modern developments in image guidance and motion management. Finally, we discuss the impact of current technological approaches on the requirements for quality assurance as well as future technological developments.

Intra- and inter-fractional liver and lung tumor motions treated with SBRT under active breathing control

Purpose

To assess intra‐ and inter‐fractional motions of liver and lung tumors using active breathing control (ABC).

Methods and Materials

Nineteen patients with liver cancer and 15 patients with lung cancer treated with stereotactic body radiotherapy (SBRT) were included in this retrospective study. All patients received a series of three CTs at simulation to test breath‐hold reproducibility. The centroids of the whole livers and of the lung tumors from the three CTs were compared to assess intra‐fraction variability. For 15 patients (8 liver, 7 lung), ABC‐gated kilovoltage cone‐beam CTs (kV‐CBCTs) were acquired prior to each treatment, and the centroids of the whole livers and of the lung tumors were also compared to those in the planning CTs to assess inter‐fraction variability.

Results

Liver intra‐fractional systematic/random errors were 0.75/0.39 mm, 1.36/0.97 mm, and 1.55/1.41 mm at medial‐lateral (ML), anterior‐posterior (AP), and superior‐inferior (SI) directions, respectively. Lung intra‐fractional systematic/random errors were 0.71/0.54 mm (ML), 1.45/1.10 mm (AP), and 3.95/1.93 mm (SI), respectively. Substantial intra‐fraction motions (>3 mm) were observed in 26.3% of liver cancer patients and in 46.7% of lung cancer patients. For both liver and lung tumors, most inter‐fractional systematic and random errors were larger than the corresponding intra‐fractional errors. However, these inter‐fractional errors were mostly corrected by the treatment team prior to each treatment based on kV CBCT‐guided soft tissue alignment, thereby eliminating their effects on the treatment planning margins.

Conclusions

Intra‐fractional motion is the key to determine the planning margins since inter‐fractional motion can be compensated based on daily gated soft tissue imaging guidance of CBCT. Patient‐specific treatment planning margins instead of recipe‐based margins were suggested, which can benefit mostly for the patients with small intra‐fractional motions.

Radiobiological analysis of stereotactic body radiation therapy for an evidence-based planning target volume of the lung using multiphase CT images obtained with a pneumatic abdominal compression apparatus: a case study

This study evaluated the efficiency of stereotactic body radiation therapy of lung (SBRT-Lung) in generating a treatment volume using conventional multiple-phase three-dimensional computed tomography (3D-CT) of a patient immobilized with pneumatic abdominal compression. The institutional protocol for SBRT-Lung using the RapidArc technique relied on a planning target volume (PTV) delineated using 3D-CT and accounted for linear and angular displacement of the tumor during respiratory movements. The efficiency of the institutional protocol was compared with that of a conventional method for PTV delineation based on radiobiological estimates, such as tumor control probability (TCP) and normal tissue complication probability (NTCP), evaluated using dose-volume parameters. Pneumatic abdominal compression improved the TCP by 15%. This novel protocol improved the TCP by 0.5% but reduced the NTCP for lung pneumonitis (0.2%) and rib fracture (1.0%). Beyond the observed variations in the patient's treatment setup, the institutional protocol yielded a significantly consistent TCP (p < 0.005). The successful clinical outcome of this case study corroborates predictions based on radiobiological evaluation and deserves validation through an increased number of patients.