Evaluation of four-dimensional computed tomography (4D-CT)-based pulmonary ventilation: The high correlation between 4D-CT ventilation and 81mKr-planar images was found

Development of a dynamic deformable thorax phantom for the quality management of deformable image registration

The purpose of this study was to develop a novel dynamic deformable thorax phantom for deformable image registration (DIR) quality assurance (QA) and to verify as a tool for commissioning and DIR QA. The phantom consists of a base phantom, an inner phantom, and a motor-derived piston. The base phantom is an acrylic cylinder phantom with a diameter of 180 mm. The inner phantom consists of deformable, 20 mm thick disk-shaped sponges. To evaluate the physical characteristics of the phantom, we evaluated its image quality and deformation. DIR accuracies were evaluated using the three types of commercially DIR software (MIM, RayStation, and Velocity AI) to test the feasibility of this phantom. We used different DIR parameters to test the impact of parameters on DIR accuracy in various phantom settings. To evaluate DIR accuracy, a target registration error (TRE) was calculated using the anatomical landmark points. The three locations (i.e., distal, middle, and proximal positions) had different displacement amounts. This result indicated that the inner phantom was not moved but deformed. In cases with different phantom settings and marker settings, the ranges of the average TRE were 0.63-15.60 mm (MIM). In cases with different DIR parameters settings, the ranges of the average TRE were as follows: 0.73-7.10 mm (MIM), 8.25-8.66 mm (RayStation), and 8.26-8.43 mm (Velocity). These results suggest that our phantom could evaluate the detailed DIR behaviors with TRE. Therefore, this is indicative of the potential usefulness of our phantom in DIR commissioning and QA.

Evaluation of four-dimensional cone beam computed tomography ventilation images acquired with two different linear accelerators at various gantry speeds using a deformable lung phantom

We evaluated four-dimensional cone beam computed tomography (4D-CBCT) ventilation images (VI_CBCT) acquired with two different linear accelerator systems at various gantry speeds using a deformable lung phantom. The 4D-CT and 4D-CBCT scans were performed using a computed tomography (CT) scanner, an X-ray volume imaging system (Elekta XVI) mounted in Versa HD, and an On-Board Imager (OBI) system mounted in TrueBeam. Intensity-based deformable image registration (DIR) was performed between peak-exhale and peak-inhale images. VI_CBCT- and 4D-CT-based ventilation images (VI_CT) were derived by DIR using two metrics: one based on the Jacobian determinant and one on changes in the Hounsfield unit (HU). Three different DIR regularization values (λ) were used for VI_CBCT. Correlations between the VI_CBCT and VI_CT values were evaluated using voxel-wise Spearman's rank correlation coefficient (r). In case of both metrics, the Jacobian-based VI_CBCT with a gantry speed of 0.6 deg/sec in Versa HD showed the highest correlation for all the gantry speeds (e.g., λ = 0.05 and r = 0.68). Thus, the r value of the Jacobian-based VI_CBCT was greater or equal to that of the HU-based VI_CBCT. In addition, the ventilation accuracy of VI_CBCT increased at low gantry speeds. Thus, the image quality of VI_CBCT was affected by the change in gantry speed in both the imaging systems. Additionally, DIR regularization considerably influenced VI_CBCT in both the imaging systems. Our results have the potential to assist in designing CBCT protocols, incorporating VI_CBCT imaging into the functional avoidance planning process.

The impact of 4DCT-ventilation imaging-guided proton therapy on stereotactic body radiotherapy for lung cancer

Functional lung avoidance during radiotherapy can help reduce pulmonary toxicity. This study assessed the potential impact of four-dimensional computed tomography (4DCT)-ventilation imaging-guided proton radiotherapy (PT) on stereotactic body radiotherapy (SBRT) by comparing it with three-dimensional conformal radiotherapy (3D-CRT) and volumetric modulated arc therapy (VMAT), which employ photon beams. Thirteen lung cancer patients who received SBRT with 3D-CRT were included in the study. 4DCT ventilation was calculated using the patients' 4DCT data, deformable image registration, and a density-change-based algorithm. Three functional treatment plans sparing the functional lung regions were developed for each patient using 3D-CRT, VMAT, and PT. The prescribed doses and dose constraints were based on the Radiation Therapy Oncology Group 0618 protocol. We evaluated the region of interest (ROI) and functional map-based dose-function metrics for 4DCT ventilation and the irradiated dose. Using 3D-CRT, VMAT, and PT, the percentages of the functional lung regions that received ≥ 5 Gy (fV5) were 26.0%, 21.9%, and 10.7%, respectively; the fV10 were 14.4%, 11.4%, and 9.0%, respectively; and fV20 were 6.5%, 6.4%, and 6.6%, respectively, and the functional mean lung doses (fMLD) were 5.6 Gy, 5.2 Gy, and 3.8 Gy, respectively. These results indicated that PT resulted in a significant reduction in fMLD, fV5, and fV10, but not fV20. The use of PT reduced the radiation to highly functional lung regions compared with those for 3D-CRT and VMAT while meeting all dose constraints.

Evaluation of the effect of user-guided deformable image registration of thoracic images on registration accuracy among users

User-guided deformable image registration (DIR) has allowed users to actively participate in the DIR process and is expected to improve DIR accuracy. The purpose of this study was to evaluate the time required for and effect of user-guided DIR on registration accuracy for thoracic images among users. In this study, 4-dimensional computed tomographic images of 10 thoracic cancer patients were used. The dataset for these patients was provided by DIR-Lab (www.dir-lab.com) and included a coordinate list of anatomical landmarks (300 bronchial bifurcations). Four medical physicists from different institutions performed DIR between peak-inhale and peak-exhale images with/without the user-guided DIR tool, Reg Refine, implemented in MIM Maestro (MIM software, Cleveland, OH). DIR accuracy was quantified by using target registration errors (TREs) for 300 anatomical landmarks in each patient. The average TREs with user-guided DIR in the 10 images by the 4 medical physicists were 1.48, 1.80, 3.46, and 3.55 mm, respectively, whereas the TREs without user-guided DIR were 3.28, 3.45, 3.56, and 3.28 mm, respectively. The average times taken by the 4 physicists to use the user-guided DIR were 10.0, 6.7, 7.1, and 8.0 min, respectively. This study demonstrated that user-guided DIR can improve DIR accuracy and requires only a moderate amount of time (<10 min). However, 2 of the 4 users did not show much improvement in DIR accuracy, which indicated the necessity of training prior to use of user-guided DIR.

CT-based ventilation imaging in radiation oncology

A form of lung function imaging is emerging that uses phase-resolved four-dimensional CT (4DCT or breath-hold CT) images along with image processing techniques to generate lung function maps that provide a surrogate of lung ventilation. CT-based ventilation (referred to as CT-ventilation) research has gained momentum in Radiation Oncology because many lung cancer patients undergo four-dimensional CT simulation as part of the standard treatment planning process. Therefore, generating CT-ventilation images provides functional information without burdening the patient with an extra imaging procedure. CT-ventilation has progressed from an image processing calculation methodology, to validation efforts, to retrospective demonstration of clinical utility in Radiation Oncology. In particular, CT-ventilation has been proposed for two main clinical applications: functional avoidance radiation therapy and thoracic dose-response assessment. The idea of functional avoidance radiation therapy is to preferentially spare functional portions of the lung (as measured by CT-ventilation) during radiation therapy with the hypothesis that reducing dose to functional portions of the lung will lead to reduced rates of radiation-related thoracic toxicity. The idea of imaging-based dose-response assessment is to evaluate pre- to post-treatment CT-ventilation-based imaging changes. The hypothesis is that early, imaging-change-based response can be an early predictor of subsequent thoracic toxicity. Based on the retrospective evidence, the clinical applications of CT-ventilation have progressed from the retrospective setting to on-going prospective clinical trials. This review will cover basic CT-ventilation calculation methodologies, validation efforts, presentation of clinical applications, summarize on-going clinical trials, review potential uncertainties and shortcomings of CT-ventilation, and discuss future directions of CT-ventilation research.

Changes in Regional Ventilation During Treatment and Dosimetric Advantages of CT Ventilation Image Guided Radiation Therapy for Locally Advanced Lung Cancer

PURPOSE

Lung 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 MATERIALS

CT 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.

RESULTS

Tumor volume was significantly reduced by 20% after 16 to 20 Gy (P = .02) and by 32% after 30 to 34 Gy (P < .01) on average. In both scenarios, the lung accumulated dose-function metrics were significantly lower in the functional plans than in the anatomic plans without compromising target volume coverage and adherence to constraints to critical structures. For example, functional planning significantly reduced the functional mean lung dose by 5.0% (P < .01) compared to anatomic planning in the adapted scenario and by 3.6% (P = .03) in the unadapted scenario.

CONCLUSIONS

This study demonstrated significant reductions in the accumulated dose to the functional lung with CT ventilation image-guided planning compared to anatomic image-guided planning for patients showing tumor regression and changes in regional ventilation during RT.

Evaluation of functionally weighted dose-volume parameters for thoracic stereotactic ablative radiotherapy (SABR) using CT ventilation

For the purpose of reducing radiation pneumontisis (RP), four-dimensional CT (4DCT)-based ventilation can be used to reduce functionally weighted lung dose. This study aimed to evaluate the functionally weighted dose-volume parameters and to investigate an optimal weighting method to realize effective planning optimization in thoracic stereotactic ablative radiotherapy (SABR). Forty patients treated with SABR were analyzed. Ventilation images were obtained from 4DCT using deformable registration and Hounsfield unit-based calculation. Functionally-weighted mean lung dose (fMLD) and functional lung fraction receiving at least x Gy (fVx) were calculated by two weighting methods: thresholding and linear weighting. Various ventilation thresholds (5th-95th, every 5th percentile) were tested. The predictive accuracy for CTCAE grade ≧ 2 pneumonitis was evaluated by area under the curve (AUC) of receiver operating characteristic analysis. AUC values varied from 0.459 to 0.570 in accordance with threshold and dose-volume metrics. A combination of 25th percentile threshold and fV₃₀ showed the best result (AUC: 0.570). AUC values with fMLD, fV₁₀, fV₂₀, and fV₄₀ were 0.541, 0.487, 0.548 and 0.563 using a 25th percentile threshold. Although conventional MLD, V₁₀, V₂₀, V₃₀ and V₄₀ showed lower AUC values (0.516, 0.477, 0.534, 0.552 and 0.527), the differences were not statistically significant. fV₃₀ with 25th percentile threshold was the best predictor of RP. Our results suggested that the appropriate weighting should be used for better treatment outcomes in thoracic SABR.

Technical Note: Correction for the effect of breathing variations in CT pulmonary ventilation imaging

PURPOSE

The accuracy and precision of computed tomography (CT) pulmonary ventilation imaging with conventional CT scanners are limited by breathing variations. We propose a method to correct for the effect of breathing variations in CT ventilation imaging based on external respiratory signals acquired throughout a scan.

METHODS

The proposed method is based on: (a) calculating voxel-by-voxel abdominal surface motion ranges using four-dimensional (4D) CT image datasets spatiotemporally correlated with external respiratory monitor data, and (b) applying the correction factor, which is defined as the ratio of the overall mean of the abdominal surface motion range in the lungs to that of each voxel, to the CT ventilation value. The performance of the proposed method was investigated by comparing voxel-wise correlations of the uncorrected and corrected CT ventilation images with single-photon emission CT (SPECT) ventilation images as a ground truth for nine patients. CT ventilation images were calculated by deformable image registration of the 4D-CT image datasets, followed by calculation of regional volume changes. A Steiger's Z-test was used to determine the statistical significance of the difference between the correlations for the uncorrected and corrected CT ventilation images.

RESULTS

The proposed correction method resulted in significant increases (P < 0.05) in the correlation between CT and SPECT ventilation in three patients, trends toward significant increase (P: 0.13-0.18) in two patients, no significant differences in three patients, and a significantly decreased correlation in one patient. The average standard deviation of the abdominal surface motion range in three patients showing significant increases was 0.27 (range 0.10-0.37), which was greater than 0.17 (range 0.07-0.38) in the other six patients.

CONCLUSIONS

The proposed method to correct for the effect of breathing variations could be readily implemented and has the potential to improve the accuracy of CT ventilation imaging as demonstrated in a nine-patient study.