Search

PurposeCT ventilation imaging (CTVI) is being used to achieve functional avoidance lung cancer radiation therapy in three clinical trials (NCT02528942, NCT02308709, NCT02843568). To address the need for common CTVI validation tools, we have built the Ventilation And Medical Pulmonary Image Registration Evaluation (VAMPIRE) Dataset, and present the results of the first VAMPIRE Challenge to compare relative ventilation distributions between different CTVI algorithms and other established ventilation imaging modalities.MethodsThe VAMPIRE Dataset includes 50 pairs of 4DCT scans and corresponding clinical or experimental ventilation scans, referred to as reference ventilation images (RefVIs). The dataset includes 25 humans imaged with Galligas 4DPET/CT, 21 humans imaged with DTPA-SPECT, and 4 sheep imaged with Xenon-CT. For the VAMPIRE Challenge, 16 subjects were allocated to a training group (with RefVI provided) and 34 subjects were allocated to a validation group (with RefVI blinded). Seven research groups downloaded the Challenge dataset and uploaded CTVIs based on deformable image registration (DIR) between the 4DCT inhale/exhale phases. Participants used DIR methods broadly classified into B-splines, Free-form, Diffeomorphisms, or Biomechanical modeling, with CT ventilation metrics based on the DIR evaluation of volume change, Hounsfield Unit change, or various hybrid approaches. All CTVIs were evaluated against the corresponding RefVI using the voxel-wise Spearman coefficient rS , and Dice similarity coefficients evaluated for low function lung ( DSClow ) and high function lung ( DSChigh ).ResultsA total of 37 unique combinations of DIR method and CT ventilation metric were either submitted by participants directly or derived from participant-submitted DIR motion fields using the in-house software, VESPIR. The rS and DSC results reveal a high degree of inter-algorithm and intersubject variability among the validation subjects, with algorithm rankings changing by up to ten positions depending on the choice of evaluation metric. The algorithm with the highest overall cross-modality correlations used a biomechanical model-based DIR with a hybrid ventilation metric, achieving a median (range) of 0.49 (0.27-0.73) for rS , 0.52 (0.36-0.67) for DSClow , and 0.45 (0.28-0.62) for DSChigh . All other algorithms exhibited at least one negative rS value, and/or one DSC value less than 0.5.ConclusionsThe VAMPIRE Challenge results demonstrate that the cross-modality correlation between CTVIs and the RefVIs varies not only with the choice of CTVI algorithm but also with the choice of RefVI modality, imaging subject, and the evaluation metric used to compare relative ventilation distributions. This variability may arise from the fact that each of the different CTVI algorithms and RefVI modalities provides a distinct physiologic measurement. Ultimately this variability, coupled with the lack of a "gold standard," highlights the ongoing importance of further validation studies before CTVI can be widely translated from academic centers to the clinic. It is hoped that the information gleaned from the VAMPIRE Challenge can help inform future validation efforts.

PurposeLung functional image guided radiation therapy (RT) that avoids irradiating highly functional regions has potential to reduce pulmonary toxicity following RT. Tumor regression during RT is common, leading to recovery of lung function. We hypothesized that computed tomography (CT) ventilation image-guided treatment planning reduces the functional lung dose compared to standard anatomic image-guided planning in 2 different scenarios with or without plan adaptation.Methods and materialsCT scans were acquired before RT and during RT at 2 time points (16-20 Gy and 30-34 Gy) for 14 patients with locally advanced lung cancer. Ventilation images were calculated by deformable image registration of four-dimensional CT image data sets and image analysis. We created 4 treatment plans at each time point for each patient: functional adapted, anatomic adapted, functional unadapted, and anatomic unadapted plans. Adaptation was performed at 2 time points. Deformable image registration was used for accumulating dose and calculating a composite of dose-weighted ventilation used to quantify the lung accumulated dose-function metrics. The functional plans were compared with the anatomic plans for each scenario separately to investigate the hypothesis at a significance level of 0.05.ResultsTumor volume was significantly reduced by 20% after 16 to 20 Gy (P = .02) and by 32% after 30 to 34 Gy (P

Background and purposeThe impact of audiovisual (AV) biofeedback on four dimensional (4D) positron emission tomography (PET) and 4D computed tomography (CT) image quality was investigated in a prospective clinical trial (NCT01172041).Material and methods4D-PET and 4D-CT images of ten lung cancer patients were acquired with AV biofeedback (AV) and free breathing (FB). The 4D-PET images were analyzed for motion artifacts by comparing 4D to 3D PET for gross tumor volumes (GTVPET) and maximum standardized uptake values (SUVmax). The 4D-CT images were analyzed for artifacts by comparing normalized cross correlation-based scores (NCCS) and quantifying a visual assessment score (VAS). A Wilcoxon signed-ranks test was used for statistical testing.ResultsThe impact of AV biofeedback varied widely. Overall, the 3D to 4D decrease of GTVPET was 1.2±1.3cm(3) with AV and 0.6±1.8cm(3) for FB. The 4D-PET increase of SUVmax was 1.3±0.9 with AV and 1.3±0.8 for FB. The 4D-CT NCCS were 0.65±0.27 with AV and 0.60±0.32 for FB (p=0.08). The 4D-CT VAS was 0.0±2.7.ConclusionThis study demonstrated a high patient dependence on the use of AV biofeedback to reduce motion artifacts in 4D imaging. None of the hypotheses tested were statistically significant. Future development of AV biofeedback will focus on optimizing the human-computer interface and including patient training sessions for improved comprehension and compliance.

PurposeTo investigate the hypothesis that CT ventilation functional image-based IMRT plans designed to avoid irradiating highly-functional lung regions are comparable to single-photon emission CT (SPECT) ventilation functional image-based plans.Methods and materialsThree IMRT plans were created for eight thoracic cancer patients using: (1) CT ventilation functional images, (2) SPECT ventilation functional images, and (3) anatomic images (no functional images). CT ventilation images were created by deformable image registration of 4D-CT image data sets and quantitative analysis. The resulting plans were analyzed for the relationship between the deviations of CT-functional plan metrics from anatomic plan metrics (ΔCT-anatomic) and those of SPECT-functional plans (ΔSPECT-anatomic), and moreover for agreements of various metrics between the CT-functional and SPECT-functional plans.ResultsThe relationship between ΔCT-anatomic and ΔSPECT-anatomic was strong (e.g., R=0.94; linear regression slope 0.71). The average differences and 95% limits of agreement between the CT-functional and SPECT-functional plan metrics (except for monitor units) for various structures were mostly less than 1% and 2%, respectively.ConclusionsThis study demonstrated a reasonable agreement between the CT ventilation functional image-based IMRT plans and SPECT-functional plans, suggesting the potential for CT ventilation imaging to serve as a surrogate for SPECT ventilation in functional image-guided radiotherapy.

Purpose: 4-dimensional computed tomography (4D-CT)-based pulmonary ventilation imaging is an emerging functional imaging modality. The purpose of this study was to investigate the physiological significance of 4D-CT ventilation imaging by comparison with pulmonary function test (PFT) measurements and single-photon emission CT (SPECT) ventilation images, which are the clinical references for global and regional lung function, respectively. Methods and Materials: In an institutional review board-approved prospective clinical trial, 4D-CT imaging and PFT and/or SPECT ventilation imaging were performed in thoracic cancer patients. Regional ventilation (V 4DCT) was calculated by deformable image registration of 4D-CT images and quantitative analysis for regional volume change. V 4DCT defect parameters were compared with the PFT measurements (forced expiratory volume in 1second (FEV1; % predicted) and FEV1/forced vital capacity (FVC; %). V 4DCT was also compared with SPECT ventilation (V SPECT) to (1) test whether V 4DCT in V SPECT defect regions is significantly lower than in nondefect regions by using the 2-tailed t test; (2) to quantify the spatial overlap between V 4DCT and V SPECT defect regions with Dice similarity coefficient (DSC); and (3) to test ventral-to-dorsal gradients by using the 2-tailed t test. Results: Of 21 patients enrolled in the study, 18 patients for whom 4D-CT and either PFT or SPECT were acquired were included in the analysis. V 4DCT defect parameters were found to have significant, moderate correlations with PFT measurements. For example, V4DCTHU defect volume increased significantly with decreasing FEV1/FVC (R=-0.65, P

PurposeThis study aims to evaluate the potential and feasibility of positron emission tomography for dynamic lung tumor tracking during radiation treatment. The authors propose a center of mass (CoM) tumor tracking algorithm using gated-PET images combined with a respiratory monitor and investigate the geometric accuracy of the proposed algorithm.MethodsThe proposed PET dynamic lung tumor tracking algorithm estimated the target position information through the CoM of the segmented target volume on gated PET images reconstructed from accumulated coincidence events. The information was continuously updated throughout a scan based on the assumption that real-time processing was supported (actual processing time at each frame ≈ 10 s). External respiratory motion and list-mode PET data were acquired from a phantom programmed to move with measured respiratory traces (external respiratory motion and internal target motion) from human subjects, for which the ground truth target position was known as a function of time. The phantom was cylindrical with six hollow sphere targets (10, 13, 17, 22, 28, and 37 mm in diameter). The measured respiratory traces consisted of two sets: (1) 1D-measured motion from ten healthy volunteers and (2) 3D-measured motion from four lung cancer patients. The authors evaluated the geometric accuracy of the proposed algorithm by quantifying estimation errors (Euclidean distance) between the actual motion of targets (1D-motion and 3D-motion traces) and CoM trajectories estimated by the proposed algorithm as a function of time.ResultsThe time-averaged error of 1D-motion traces over all trajectories of all targets was 1.6 mm. The error trajectories decreased with time as coincidence events were accumulated. The overall error trajectory of 1D-motion traces converged to within 2 mm in approximately 90 s. As expected, more accurate results were obtained for larger targets. For example, for the 37 mm target, the average error over all 1D-motion traces was 1.1 mm; and for the 10 mm target, the average error over all 1D-motion traces was 2.8 mm. The overall time-averaged error of 3D-motion traces was 1.6 mm, which was comparable to that of the 1D-motion traces. There were small variations in the errors between the 3D-motion traces, although the motion trajectories were very different. The accuracy of the estimates was consistent for all targets except for the smallest.ConclusionsThe authors developed an algorithm for dynamic lung tumor tracking using list-mode PET data and a respiratory motion signal, and demonstrated proof-of-principle for PET-guided lung tumor tracking. The overall tracking error in phantom studies is less than 2 mm.

PurposeFour-dimensional (4D) computed tomography (CT) ventilation imaging is a novel promising technique for lung functional imaging. The current standard 4D CT technique using phase-based sorting frequently results in artifacts, which may deteriorate the accuracy of ventilation imaging. The purpose of this study was to quantify the variability of 4D CT ventilation imaging due to 4D CT sorting.Methods4D CT image sets from nine lung cancer patients were each sorted by the phase-based method and anatomic similarity-based method, designed to reduce artifacts, with corresponding ventilation images created for each method. Artifacts in the resulting 4D CT images were quantified with the artifact score which was defined based on the difference between the normalized cross correlation for CT slices within a CT data segment and that for CT slices bordering the interface between adjacent CT data segments. The ventilation variation was quantified using voxel-based Spearman rank correlation coefficients for all lung voxels, and Dice similarity coefficients (DSC) for the spatial overlap of low-functional lung volumes. Furthermore, the correlations with matching single-photon emission CT (SPECT) ventilation images (assumed ground truth) were evaluated for three patients to investigate which sorting method provides higher physiologic accuracy.ResultsAnatomic similarity-based sorting reduced 4D CT artifacts compared to phase-based sorting (artifact score, 0.45 ± 0.14 vs 0.58 ± 0.24, p = 0.10 at peak-exhale; 0.63 ± 0.19 vs 0.71 ± 0.31, p = 0.25 at peak-inhale). The voxel-based correlation between the two ventilation images was 0.69 ± 0.26 on average, ranging from 0.03 to 0.85. The DSC was 0.71 ± 0.13 on average. Anatomic similarity-based sorting yielded significantly fewer lung voxels with paradoxical negative ventilation values than phase-based sorting (5.0 ± 2.6% vs 9.7 ± 8.4%, p = 0.05), and improved the correlation with SPECT ventilation regionally.ConclusionsThe variability of 4D CT ventilation imaging due to 4D CT sorting was moderate overall and substantial in some cases, suggesting that 4D CT artifacts are an important source of variations in 4D CT ventilation imaging. Reduction of 4D CT artifacts provided more physiologically convincing and accurate ventilation estimates. Further studies are needed to confirm this result.

Rationale and objectivesA novel ventilation imaging method based on four-dimensional (4D) computed tomography (CT) has been applied to the field of radiation oncology. Understanding its reproducibility is a prerequisite for clinical applications. The purpose of this study was to quantify the reproducibility of 4D CT ventilation imaging over different days and the same session.Materials and methodsTwo ventilation images were created from repeat 4D CT scans acquired over the average time frames of 15 days for 6 lung cancer patients and 5 minutes for another 6 patients. The reproducibility was quantified using the voxel-based Spearman rank correlation coefficients for all lung voxels and Dice similarity coefficients (DSC) for the spatial overlap of segmented high-, moderate-, and low-functional lung volumes. Furthermore, the relationship between the variation in abdominal motion range as a measure of the depth of breathing and variation in ventilation was evaluated using linear regression.ResultsThe voxel-based correlation between the two ventilation images was moderate on average (0.50 ± 0.15). The DSCs were also moderate for the high- (0.60 ± 0.08), moderate- (0.46 ± 0.06), and low-functional lung (0.58 ± 0.09). No patients demonstrated strong correlations. The relationship between the motion range variation and ventilation variation was found to be moderate and significant.ConclusionsWe investigated the reproducibility of 4D CT ventilation imaging over the time frames of 15 days and 5 minutes and found that it was only moderately reproducible. Respiratory variation during 4D CT scans was found to deteriorate the reproducibility. Improvement of 4D CT imaging is necessary to increase the reproducibility of 4D CT ventilation imaging.

Background: During prostate stereotactic body radiation therapy (SBRT), prostate tumor translational motion may deteriorate the planned dose distribution. Most of the major advances in motion management to date have focused on correcting this one aspect of the tumor motion, translation. However, large prostate rotation up to 30° has been measured. As the technological innovation evolves toward delivering increasingly precise radiotherapy, it is important to quantify the clinical benefit of translational and rotational motion correction over translational motion correction alone. Purpose: The purpose of this work was to quantify the dosimetric impact of intrafractional dynamic rotation of the prostate measured with a six degrees‐of‐freedom tumor motion monitoring technology. Methods: The delivered dose was reconstructed including (a) translational and rotational motion and (b) only translational motion of the tumor for 32 prostate cancer patients recruited on a 5‐fraction prostate SBRT clinical trial. Patients on the trial received 7.25 Gy in a treatment fraction. A 5 mm clinical target volume (CTV) to planning target volume (PTV) margin was applied in all directions except the posterior direction where a 3 mm expansion was used. Prostate intrafractional translational motion was managed using a gating strategy, and any translation above the gating threshold was corrected by applying an equivalent couch shift. The residual translational motion is denoted as Tres$T_{res}$. Prostate intrafractional rotational motion Runcorr$R_{uncorr}$ was recorded but not corrected. The dose differences from the planned dose due to Tres$T_{res}$ + Runcorr$R_{uncorr}$, ΔD(Tres$T_{res}$ + Runcorr$R_{uncorr}$) and due to Tres$T_{res}$ alone, ΔD(Tres$T_{res}$), were then determined for CTV D98, PTV D95, bladder V6Gy, and rectum V6Gy. The residual dose error due to uncorrected rotation, Runcorr$R_{uncorr}$ was then quantified: ΔDResidual$\Delta D_{Residual}$ = ΔD(Tres$T_{res}$ + Runcorr$R_{uncorr}$) ‐ ΔD(Tres${T}_{\textit{res}}$). Results: Fractional data analysis shows that the dose differences from the plan (both ΔD(Tres$T_{res}$ + Runcorr$R_{uncorr}$) and ΔD(Tres$T_{res}$)) for CTV D98 was less than 5% in all treatment fractions. ΔD(Tres$T_{res}$ + Runcorr$R_{uncorr}$) was larger than 5% in one fraction for PTV D95, in one fraction for bladder V6Gy, and in five fractions for rectum V6Gy. Uncorrected rotation, Runcorr$R_{uncorr}$ induced residual dose error, ΔDResidual$\Delta D_{Residual}$, resulted in less dose to CTV and PTV in 43% and 59% treatment fractions, respectively, and more dose to bladder and rectum in 51% and 53% treatment fractions, respectively. The cumulative dose over five fractions, ∑D(Tres$T_{res}$ + Runcorr$R_{uncorr}$) and ∑D(Tres$T_{res}$), was always within 5% of the planned dose for all four structures for every patient. Conclusions: The dosimetric impact of tumor rotation on a large prostate cancer patient cohort was quantified in this study. These results suggest that the standard 3–5 mm CTV‐PTV margin was sufficient to account for the intrafraction prostate rotation observed for this cohort of patients, provided an appropriate gating threshold was applied to correct for translational motion. Residual dose errors due to uncorrected prostate rotation were small in magnitude, which may be corrected using different treatment adaptation strategies to further improve the dosimetric accuracy. [ABSTRACT FROM AUTHOR]