Reproducibility of Tumor Motion Probability Distribution Function in Stereotactic Body Radiation Therapy of Lung Cancer

Data-driven rapid 4D cone-beam CT reconstruction for new generation linacs

Objective.Newer generation linear accelerators (Linacs) allow 20 s cone-beam CT (CBCT) acquisition which reduces radiation therapy treatment time. However, the current clinical application of these rapid scans is only 3DCBCT. In this paper we propose a novel data-driven rapid 4DCBCT reconstruction method for new generation linacs.Approach.This method relies on estimating the magnitude of the diaphragm motion from an initial 3D reconstruction. This estimated motion is used to linearly approximate a deformation vector field (DVF) for each respiration phase. These DVFs are then used for motion compensated Feldkamp-Davis-Kress (MCFDK) reconstructions. This method, named MCFDK Data Driven (MCFDK-DD), was compared to a MCFDK reconstruction using a prior motion model (MCFDK-Prior), a 3D-FDK reconstruction, and a conventional acquisition (4 mins) conventional reconstruction 4DCBCT (4D-FDK). The data used in this paper were derived from 4DCT volumes from 12 patients from The Cancer Imaging Archives. Image quality was quantified using RMSE of line plots centred on the tumour, tissue interface width (TIW), the mean square error (MSE) and structural similarity index measurement (SSIM).Main Results.The tumour line plots in the Superior-Inferior direction showed reduced RMSE for the MCFDK-DD compared to the 3D-FDK method, indicating the MCFDK-DD method provided a more accurate tumour location. Similarly, the TIW values from the MCFDK-DD reconstructions (median 8.6 mm) were significantly reduced for the MCFDK-DD method compared to the 3D-FDK reconstructions (median 14.8 mm, (p< 0.001). The MCFDK-DD, MCFDK-Prior and 3D-FDK had median MSE values of1.08×10-6mm-1,1.11×10-6mm-1and1.17×10-6mm-1respectively. The corresponding median SSIM values were 0.93, 0.92 and 0.92 respectively indicating the MCFDK-DD had good agreement with the conventional 4D-FDK reconstructions.Significance.These results demonstrate the feasibility of creating accurate data-driven 4DCBCT images for rapid scans on new generation linacs. These findings could lead to increased clinical usage of 4D information on newer generation linacs.

A Simulation Study of Tolerance of Breathing Amplitude Variations in Radiotherapy of Lung Cancer Using 4DCT and Time-Resolved 4DMRI

As patient breathing irregularities can introduce a large uncertainty in targeting the internal tumor volume (ITV) of lung cancer patients, and thereby affect treatment quality, this study evaluates dose tolerance of tumor motion amplitude variations in ITV-based volumetric modulated arc therapy (VMAT). A motion-incorporated planning technique was employed to simulate treatment delivery of 10 lung cancer patients' clinical VMAT plans using original and three scaling-up (by 0.5, 1.0, and 2.0 cm) motion waveforms from single-breath four-dimensional computed tomography (4DCT) and multi-breath time-resolved 4D magnetic resonance imaging (TR-4DMRI). The planning tumor volume (PTV = ITV + 5 mm margin) dose coverage (PTV D95%) was evaluated. The repeated waveforms were used to move the isocenter in sync with the clinical leaf motion and gantry rotation. The continuous VMAT arcs were broken down into many static beam fields at the control points (2°-interval) and the composite plan represented the motion-incorporated VMAT plan. Eight motion-incorporated plans per patient were simulated and the plan with the native 4DCT waveform was used as a control. The first (D95% ≤ 95%) and second (D95% ≤ 90%) plan breaching points due to motion amplitude increase were identified and analyzed. The PTV D95% in the motion-incorporated plans was 99.4 ± 1.0% using 4DCT, closely agreeing with the corresponding ITV-based VMAT plan (PTV D95% = 100%). Tumor motion irregularities were observed in TR-4DMRI and triggered D95% ≤ 95% in one case. For small tumors, 4 mm extra motion triggered D95% ≤ 95%, and 6-8 mm triggered D95% ≤ 90%. For large tumors, 14 mm and 21 mm extra motions triggered the first and second breaching points, respectively. This study has demonstrated that PTV D95% breaching points may occur for small tumors during treatment delivery. Clinically, it is important to monitor and avoid systematic motion increase, including baseline drift, and large random motion spikes through threshold-based beam gating.

Image-guidance technique comparison on respiratory reproducibility and dose indexes for stereotactic body radiotherapy in lung tumor

We investigated respiratory reproducibility from position errors of gold internal fiducial markers for breath-hold (BH) and real-time tumor tracking (RTT) techniques for stereotactic body radiotherapy in lung tumors. The relationship between position errors and dose indexes was checked for both techniques. The stereotactic body radiotherapy plan in lung tumors was planned for 29 patients. The tumor positioning was arranged using 1.5 mm diameter gold internal fiducial markers. First, CT images were acquired to analyze position errors of gold markers for BH and RTT techniques. The offset plans for both techniques were calculated by displacing the mean position errors. The dose indexes (D98, D95, D2, mean dose) in a planning target volume were evaluated from dose volume histograms for the original plan, BH, and RTT offset plans. The relationship between position errors and dose indexes was analyzed using the root mean square (RMS) for both techniques. For the BH, the RMS was 3.29 mm at the lower lobe. Similarly, it was 1.34 mm for the RTT. The difference for D98 by position error for BH was -7.0 ± 10.8% at the lower lobe and the difference of all dose indexes for the RTT was less than 1%. The D2 and mean dose for both techniques were nearly the same as those of the original plan. In conclusion, the adaptation of the BH technique should be ≤2 mm RMS. If the position error is >2 mm RMS, the RTT technique should be used instead of the BH technique.

Introduction of a pseudo demons force to enhance deformation range for robust reconstruction of super-resolution time-resolved 4DMRI

PURPOSE

The purpose of this study was to enhance the deformation range of demons-based deformable image registration (DIR) for large respiration-induced organ motion in the reconstruction of time-resolved four-dimensional magnetic resonance imaging (TR-4DMRI) for multi-breath motion simulation.

METHODS

A demons-based DIR algorithm was modified to enhance the deformation range for TR-4DMRI reconstruction using the super-resolution approach. A pseudo demons force was introduced to accelerate the coarse deformation in a multi-resolution (n = 3) DIR approach. The intensity gradient of a voxel was applied to its neighboring (5 × 5 × 5) voxels with a weight of Gaussian probability profile (σ = 1 voxel) to extend the demons force, especially on those voxels that have little intensity gradience but high-intensity difference. A digital 4DMRI phantom with 3-8 cm diaphragmatic motions was used for DIR comparison. Six volunteers were scanned with two high-resolution (highR: 2 × 2 × 2 mm³ ) breath-hold (BH) 3DMR images at full inhalation (BHI) and full exhalation (BHE) and low-resolution (lowR: 5 × 5 × 5 mm³ ) free-breathing (FB) 3DMR cine images (2 Hz) under an IRB-approved protocol. A cross-consistency check (CCC) (BHI→FB←BHE), with voxel intensity correlation (VIC) and inverse consistency error (ICE), was introduced for cross-verification of TR-4DMRI reconstruction.

RESULTS

Using the digital phantom, the maximum deformable magnitude is doubled using the modified DIR from 3 to 6 cm at the diaphragm. In six human subjects, the first 15-iteration DIR using the pseudo force deforms 200 ± 150% more than the original force, and succeeds in all 12 cases, whereas the original demons-based DIR failed in 67% of tested cases. Using the pseudo force, high VIC (>0.9) and small ICE (1.6 ± 0.6 mm) values are observed for DIR of BHI&BHE, BHI→FB, and BHE→FB. The CCC identifies four questionable cases, in which two cases need further DIR refinement, without missing true negative.

CONCLUSIONS

The introduction of a pseudo demons force enhances the largest deformation magnitude up to 6 cm. The cross-consistency check ensures the quality of TR-4DMRI reconstruction. Further investigation is ongoing to fully characterize TR-4DMRI for potential multi-breathing-cycle radiotherapy simulation.

Accelerating volumetric cine MRI (VC-MRI) using undersampling for real-time 3D target localization/tracking in radiation therapy: a feasibility study

PURPOSE

To accelerate volumetric cine MRI (VC-MRI) using undersampled 2D-cine MRI to provide real-time 3D guidance for gating/target tracking in radiotherapy.

METHODS

4D-MRI is acquired during patient simulation. One phase of the prior 4D-MRI is selected as the prior images, designated as MRIprior. The on-board VC-MRI at each time-step is considered a deformation of the MRIprior. The deformation field map is represented as a linear combination of the motion components extracted by principal component analysis from the prior 4D-MRI. The weighting coefficients of the motion components are solved by matching the corresponding 2D-slice of the VC-MRI with the on-board undersampled 2D-cine MRI acquired. Undersampled Cartesian and radial k-space acquisition strategies were investigated. The effects of k-space sampling percentage (SP) and distribution, tumor sizes and noise on the VC-MRI estimation were studied. The VC-MRI estimation was evaluated using XCAT simulation of lung cancer patients and data from liver cancer patients. Volume percent difference (VPD) and Center of Mass Shift (COMS) of the tumor volumes and tumor tracking errors were calculated.

RESULTS

For XCAT, VPD/COMS were 11.93  ±  2.37%/0.90  ±  0.27 mm and 11.53  ±  1.47%/0.85  ±  0.20 mm among all scenarios with Cartesian sampling (SP  =  10%) and radial sampling (21 spokes, SP  =  5.2%), respectively. When tumor size decreased, higher sampling rate achieved more accurate VC-MRI than lower sampling rate. VC-MRI was robust against noise levels up to SNR  =  20. For patient data, the tumor tracking errors in superior-inferior, anterior-posterior and lateral (LAT) directions were 0.46  ±  0.20 mm, 0.56  ±  0.17 mm and 0.23  ±  0.16 mm, respectively, for Cartesian-based sampling with SP  =  20% and 0.60  ±  0.19 mm, 0.56  ±  0.22 mm and 0.42  ±  0.15 mm, respectively, for radial-based sampling with SP  =  8% (32 spokes).

CONCLUSIONS

It is feasible to estimate VC-MRI from a single undersampled on-board 2D cine MRI. Phantom and patient studies showed that the temporal resolution of VC-MRI can potentially be improved by 5-10 times using a 2D cine image acquired with 10-20% k-space sampling.

Assessment of Intrafraction Breathing Motion on Left Anterior Descending Artery Dose During Left-Sided Breast Radiation Therapy

PURPOSE

To use 4-dimensional computed tomography (4D-CT) imaging to predict the level of uncertainty in cardiac dose estimates of the left anterior descending artery that arises due to breathing motion during radiation therapy for left-sided breast cancer.

METHODS AND MATERIALS

The fast helical CT (FH-CT) and 4D-CT of 30 left-sided breast cancer patients were retrospectively analyzed. Treatment plans were created on the FH-CT. The original treatment plan was then superimposed onto all 10 phases of the 4D-CT to quantify the dosimetric impact of respiratory motion through 4D dose accumulation (4D-dose). Dose-volume histograms for the heart, left ventricle (LV), and left anterior descending (LAD) artery obtained from the FH-CT were compared with those obtained from the 4D-dose.

RESULTS

The 95% confidence interval of 4D-dose and FH-CT differences in mean dose estimates for the heart, LV, and LAD were ±0.5 Gy, ±1.0 Gy, and ±8.7 Gy, respectively.

CONCLUSION

Fast helical CT is a good approximation for doses to the heart and LV; however, dose estimates for the LAD are susceptible to uncertainties that arise due to intrafraction breathing motion that cannot be ascertained without the additional information obtained from 4D-CT and dose accumulation. For future clinical studies, we suggest the use of 4D-CT-derived dose-volume histograms for estimating the dose to the LAD.

Megavoltage conebeam CT cine as final verification of treatment plan in lung stereotactic body radiotherapy

INTRODUCTION

To analyse the clinical impact of megavoltage conebeam computed tomography (MV-CBCT) cine on internal target volume (ITV) coverage in lung stereotactic body radiotherapy (SBRT).

METHODS

One hundred and six patients received lung SBRT. All underwent 4D computed tomography simulation followed by treatment via image guided 3D conformal or intensity modulated radiation. Prior to SBRT, all patients underwent MV-CBCT cine, in which raw projections are displayed as beam's-eye-view fluoroscopic series with the planning target volume (PTV) projected onto each image, enabling verification of tumour motion relative to the PTV and assessment of adequacy of treatment margin.

RESULTS

Megavoltage conebeam computed tomography cine was completed 1-2 days prior to SBRT. Four patients (3.8%) had insufficient ITV coverage inferiorly at cine review. All four plans were changed by adding 5 mm on the PTV margin inferiorly. The mean change in PTV volumes was 3.9 cubic centimetres (cc) (range 1.85-6.32 cc). Repeat cine was performed after plan modification to ensure adequate PTV coverage in the modified plans.

CONCLUSIONS

PTV margin was adequate in the majority of patients with this technique. MV-CBCT cine did show insufficient coverage in a small subset of patients. Insufficient PTV margins may be a function of 4D CT simulation inadequacies or deficiencies in visualizing the ITV inferior border in the full-inhale phase. MV-CBCT cine is a valuable tool for final verification of PTV margins.

A technique for estimating 4D-CBCT using prior knowledge and limited-angle projections

PURPOSE

To develop a technique to estimate onboard 4D-CBCT using prior information and limited-angle projections for potential 4D target verification of lung radiotherapy.

METHODS

Each phase of onboard 4D-CBCT is considered as a deformation from one selected phase (prior volume) of the planning 4D-CT. The deformation field maps (DFMs) are solved using a motion modeling and free-form deformation (MM-FD) technique. In the MM-FD technique, the DFMs are estimated using a motion model which is extracted from planning 4D-CT based on principal component analysis (PCA). The motion model parameters are optimized by matching the digitally reconstructed radiographs of the deformed volumes to the limited-angle onboard projections (data fidelity constraint). Afterward, the estimated DFMs are fine-tuned using a FD model based on data fidelity constraint and deformation energy minimization. The 4D digital extended-cardiac-torso phantom was used to evaluate the MM-FD technique. A lung patient with a 30 mm diameter lesion was simulated with various anatomical and respirational changes from planning 4D-CT to onboard volume, including changes of respiration amplitude, lesion size and lesion average-position, and phase shift between lesion and body respiratory cycle. The lesions were contoured in both the estimated and "ground-truth" onboard 4D-CBCT for comparison. 3D volume percentage-difference (VPD) and center-of-mass shift (COMS) were calculated to evaluate the estimation accuracy of three techniques: MM-FD, MM-only, and FD-only. Different onboard projection acquisition scenarios and projection noise levels were simulated to investigate their effects on the estimation accuracy.

RESULTS

For all simulated patient and projection acquisition scenarios, the mean VPD (±S.D.)∕COMS (±S.D.) between lesions in prior images and "ground-truth" onboard images were 136.11% (±42.76%)∕15.5 mm (±3.9 mm). Using orthogonal-view 15°-each scan angle, the mean VPD∕COMS between the lesion in estimated and "ground-truth" onboard images for MM-only, FD-only, and MM-FD techniques were 60.10% (±27.17%)∕4.9 mm (±3.0 mm), 96.07% (±31.48%)∕12.1 mm (±3.9 mm) and 11.45% (±9.37%)∕1.3 mm (±1.3 mm), respectively. For orthogonal-view 30°-each scan angle, the corresponding results were 59.16% (±26.66%)∕4.9 mm (±3.0 mm), 75.98% (±27.21%)∕9.9 mm (±4.0 mm), and 5.22% (±2.12%)∕0.5 mm (±0.4 mm). For single-view scan angles of 3°, 30°, and 60°, the results for MM-FD technique were 32.77% (±17.87%)∕3.2 mm (±2.2 mm), 24.57% (±18.18%)∕2.9 mm (±2.0 mm), and 10.48% (±9.50%)∕1.1 mm (±1.3 mm), respectively. For projection angular-sampling-intervals of 0.6°, 1.2°, and 2.5° with the orthogonal-view 30°-each scan angle, the MM-FD technique generated similar VPD (maximum deviation 2.91%) and COMS (maximum deviation 0.6 mm), while sparser sampling yielded larger VPD∕COMS. With equal number of projections, the estimation results using scattered 360° scan angle were slightly better than those using orthogonal-view 30°-each scan angle. The estimation accuracy of MM-FD technique declined as noise level increased.

CONCLUSIONS

The MM-FD technique substantially improves the estimation accuracy for onboard 4D-CBCT using prior planning 4D-CT and limited-angle projections, compared to the MM-only and FD-only techniques. It can potentially be used for the inter/intrafractional 4D-localization verification.