Integrating local and distant radiation-induced lung injury: Development and validation of a predictive model for ventilation loss

BACKGROUND

Investigations on radiation-induced lung injury (RILI) have predominantly focused on local effects, primarily those associated with radiation damage to lung parenchyma. However, recent studies from our group and others have revealed that radiation-induced damage to branching serial structures such as airways and vessels may also have a substantial impact on post-radiotherapy (RT) lung function. Furthermore, recent results from multiple functional lung avoidance RT trials, although promising, have demonstrated only modest toxicity reduction, likely because they were primarily focused on dose avoidance to lung parenchyma. These observations emphasize the critical need for predictive dose-response models that effectively incorporate both local and distant RILI effects.

PURPOSE

We develop and validate a predictive model for ventilation loss after lung RT. This model, referred to as P+A, integrates local (parenchyma [P]) and distant (central and peripheral airways [A]) radiation-induced damage, modeling partial (narrowing) and complete (collapse) obstruction of airways.

METHODS

In an IRB-approved prospective study, pre-RT breath-hold CTs (BHCTs) and pre- and one-year post-RT 4DCTs were acquired from lung cancer patients treated with definitive RT. Up to 13 generations of airways were automatically segmented on the BHCTs using a research virtual bronchoscopy software. Ventilation maps derived from the 4DCT scans were utilized to quantify pre- and post-RT ventilation, serving, respectively, as input data and reference standard (RS) in model validation. To predict ventilation loss solely due to parenchymal damage (referred to as P model), we used a normal tissue complication probability (NTCP) model. Our model used this NTCP-based estimate and predicted additional loss due radiation-induced partial or complete occlusion of individual airways, applying fluid dynamics principles and a refined version of our previously developed airway radiosensitivity model. Predictions of post-RT ventilation were estimated in the sublobar volumes (SLVs) connected to the terminal airways. To validate the model, we conducted a k-fold cross-validation. Model parameters were optimized as the values that provided the lowest root mean square error (RMSE) between predicted post-RT ventilation and the RS for all SLVs in the training data. The performance of the P+A and the P models was evaluated by comparing their respective post-RT ventilation values with the RS predictions. Additional evaluation using various receiver operating characteristic (ROC) metrics was also performed.

RESULTS

We extracted a dataset of 560 SLVs from four enrolled patients. Our results demonstrated that the P+A model consistently outperformed the P model, exhibiting RMSEs that were nearly half as low across all patients (13 ± 3 percentile for the P+A model vs. 24 ± 3 percentile for the P model on average). Notably, the P+A model aligned closely with the RS in ventilation loss distributions per lobe, particularly in regions exposed to doses ≥13.5 Gy. The ROC analysis further supported the superior performance of the P+A model compared to the P model in sensitivity (0.98 vs. 0.07), accuracy (0.87 vs. 0.25), and balanced predictions.

CONCLUSIONS

These early findings indicate that airway damage is a crucial factor in RILI that should be included in dose-response modeling to enhance predictions of post-RT lung function.

Systematic Evaluation of the Impact of Lung Segmentation Methods on 4-Dimensional Computed Tomography Ventilation Imaging Using a Large Patient Database

PURPOSE

A novel form of lung functional imaging applied for functional avoidance radiation therapy has been developed that uses 4-dimensional computed tomography (4DCT) data and image processing techniques to calculate lung ventilation (4DCT-ventilation). Lung segmentation is a common step to define a region of interest for 4DCT-ventilation generation. The purpose of this study was to quantitatively evaluate the sensitivity of 4DCT-ventilation imaging using different lung segmentation methods.

METHODS AND MATERIALS

The 4DCT data of 350 patients from 2 institutions were used. Lung contours were generated using 3 methods: (1) reference segmentations that removed airways and pulmonary vasculature manually (Lung-Manual), (2) standard lung contours used for planning (Lung-RadOnc), and (3) artificial intelligence (AI)-based contours that removed the airways and pulmonary vasculature (Lung-AI). The AI model was based on a residual 3-dimensional U-Net and was trained using the Lung-Manual contours of 279 patients. We compared the Lung-RadOnc or Lung-AI with Lung-Manual contours for the entire 4DCT-ventilation functional avoidance process including lung segmentation (surface Dice similarity coefficient [Surface DSC]), 4DCT-ventilation generation (correlation), and subanalysis of 10 patients on a dosimetric endpoint (percentage of high functional volume of lung receiving ≥20 Gy [fV20{%}]).

RESULTS

Surface DSC comparing Lung-Manual/Lung-RadOnc and Lung-Manual/Lung-AI contours was 0.40 ± 0.06 and 0.86 ± 0.04, respectively. The correlation between 4DCT-ventilation images generated with Lung-Manual/Lung-RadOnc and Lung-Manual/Lung-AI were 0.48 ± 0.14 and 0.85 ± 0.14, respectively. The difference in fV20[%] between 4DCT-ventilation generated with Lung-Manual/Lung-RadOnc and Lung-Manual/Lung-AI was 2.5% ± 4.1% and 0.3% ± 0.5%, respectively.

CONCLUSIONS

Our work showed that using standard planning lung contours can result in significantly variable 4DCT-ventilation images. The study demonstrated that AI-based segmentations generate lung contours and 4DCT-ventilation images that are similar to those generated using manual methods. The significance of the study is that it characterizes the lung segmentation sensitivity of the 4DCT-ventilation process and develops methods that can facilitate the integration of this novel imaging in busy clinics.

Deep learning-based combination of [18F]-FDG PET and CT images for producing pulmonary perfusion image

BACKGROUND

The main application of [18F] FDG-PET (18 FDG-PET) and CT images in oncology is tumor identification and quantification. Combining PET and CT images to mine pulmonary perfusion information for functional lung avoidance radiation therapy (FLART) is desirable but remains challenging.

PURPOSE

To develop a deep-learning-based (DL) method to combine 18 FDG-PET and CT images for producing pulmonary perfusion images (PPI).

METHODS

Pulmonary technetium-99 m-labeled macroaggregated albumin SPECT (PPISPECT ), 18 FDG-PET, and CT images obtained from 53 patients were enrolled. CT and PPISPECT images were rigidly registered, and registration displacement was subsequently used to align 18 FDG-PET and PPISPECT images. The left/right lung was separated and rigidly registered again to improve the registration accuracy. A DL model based on 3D Unet architecture was constructed to directly combine multi-modality 18 FDG-PET and CT images for producing PPI (PPIDLM ). 3D Unet architecture was used as the basic architecture, and the input was expanded from a single-channel to a dual-channel to combine multi-modality images. For comparative evaluation, 18 FDG-PET images were also used alone to generate PPIDLPET . Sixty-seven samples were randomly selected for training and cross-validation, and 36 were used for testing. The Spearman correlation coefficient (rs ) and multi-scale structural similarity index measure (MS-SSIM) between PPIDLM /PPIDLPET and PPISPECT were computed to assess the statistical and perceptual image similarities. The Dice similarity coefficient (DSC) was calculated to determine the similarity between high-/low- functional lung (HFL/LFL) volumes.

RESULTS

The voxel-wise rs and MS-SSIM of PPIDLM /PPIDLPET were 0.78 ± 0.04/0.57 ± 0.03, 0.93 ± 0.01/0.89 ± 0.01 for cross-validation and 0.78 ± 0.11/0.55 ± 0.18, 0.93 ± 0.03/0.90 ± 0.04 for testing. PPIDLM /PPIDLPET achieved averaged DSC values of 0.78 ± 0.03/0.64 ± 0.02 for HFL and 0.83 ± 0.01/0.72 ± 0.03 for LFL in the training dataset and 0.77 ± 0.11/0.64 ± 0.12, 0.82 ± 0.05/0.72 ± 0.06 in the testing dataset. PPIDLM yielded a stronger correlation and higher MS-SSIM with PPISPECT than PPIDLPET (p < 0.001).

CONCLUSIONS

The DL-based method integrates lung metabolic and anatomy information for producing PPI and significantly improved the accuracy over methods based on metabolic information alone. The generated PPIDLM can be applied for pulmonary perfusion volume segmentation, which is potentially beneficial for FLART treatment plan optimization.

Deep learning-based lung image registration: A review

Lung image registration can effectively describe the relative motion of lung tissues, thereby helping to solve series problems in clinical applications. Since the lungs are soft and fairly passive organs, they are influenced by respiration and heartbeat, resulting in discontinuity of lung motion and large deformation of anatomic features. This poses great challenges for accurate registration of lung image and its applications. The recent application of deep learning (DL) methods in the field of medical image registration has brought promising results. However, a versatile registration framework has not yet emerged due to diverse challenges of registration for different regions of interest (ROI). DL-based image registration methods used for other ROI cannot achieve satisfactory results in lungs. In addition, there are few review articles available on DL-based lung image registration. In this review, the development of conventional methods for lung image registration is briefly described and a more comprehensive survey of DL-based methods for lung image registration is illustrated. The DL-based methods are classified according to different supervision types, including fully-supervised, weakly-supervised and unsupervised. The contributions of researchers in addressing various challenges are described, as well as the limitations of these approaches. This review also presents a comprehensive statistical analysis of the cited papers in terms of evaluation metrics and loss functions. In addition, publicly available datasets for lung image registration are also summarized. Finally, the remaining challenges and potential trends in DL-based lung image registration are discussed.

Metrics of dose to highly ventilated lung are predictive of radiation-induced pneumonitis in lung cancer patients

PURPOSE

To identify metrics of radiation dose delivered to highly ventilated lung that are predictive of radiation-induced pneumonitis.

METHODS AND MATERIALS

A cohort of 90 patients with locally advanced non-small cell lung cancer treated with standard fractionated radiation therapy (RT) (60-66 Gy in 30-33 fractions) were evaluated. Regional lung ventilation was determined from pre-RT 4-dimensional computed tomography (4DCT) using the Jacobian determinant of a B-spline deformable image registration to estimate lung tissue expansion during respiration. Multiple voxel-wise population- and individual-based thresholds for defining high functioning lung were considered. Mean dose and volumes receiving dose ≥ 5-60 Gy were analyzed for both total lung-ITV (MLD,V5-V60) and highly ventilated functional lung-ITV (fMLD,fV5-fV60). The primary endpoint was symptomatic grade 2+ (G2+) pneumonitis. Receiver operator curve (ROC) analyses were used to identify predictors of pneumonitis.

RESULTS

G2+ pneumonitis occurred in 22.2% of patients, with no differences between stage, smoking status, COPD, or chemo/immunotherapy use between G<2 and G2+ patients (P≥ 0.18). Highly ventilated lung was defined as voxels exceeding the population-wide median of 18% voxel-level expansion. All total and functional metrics were significantly different between patients with and without pneumonitis (P≤ 0.039). Optimal ROC points predicting pneumonitis from functional lung dose were fMLD ≤ 12.3 Gy, fV5 ≤ 54% and fV20 ≤ 19 %. Patients with fMLD ≤ 12.3 Gy had a 14% risk of developing G2+ pneumonitis whereas risk significantly increased to 35% for those with fMLD > 12.3 Gy (P = 0.035).

CONCLUSIONS

Dose to highly ventilated lung is associated with symptomatic pneumonitis and treatment planning strategies should focus on limiting dose to functional regions. These findings provide important metrics to be used in functional lung avoidance RT planning and designing clinical trials.

Radiation pneumonitis prediction model with integrating multiple dose-function features on 4DCT ventilation images

PURPOSE

Radiation pneumonitis (RP) is dose-limiting toxicity for non-small-cell cancer (NSCLC). This study developed an RP prediction model by integrating dose-function features from computed four-dimensional computed tomography (4DCT) ventilation using the least absolute shrinkage and selection operator (LASSO).

METHODS

Between 2013 and 2020, 126 NSCLC patients were included in this study who underwent a 4DCT scan to calculate ventilation images. We computed two sets of candidate dose-function features from (1) the percentage volume receiving > 20 Gy or the mean dose on the functioning zones determined with the lower cutoff percentile ventilation value, (2) the functioning zones determined with lower and upper cutoff percentile ventilation value using 4DCT ventilation images. An RP prediction model was developed by LASSO while simultaneously determining the regression coefficient and feature selection through fivefold cross-validation.

RESULTS

We found 39.3 % of our patients had a ≥ grade 2 RP. The mean area under the curve (AUC) values for the developed models using clinical, dose-volume, and dose-function features with a lower cutoff were 0.791, and the mean AUC values with lower and upper cutoffs were 0.814. The relative regression coefficient (RRC) on dose-function features with upper and lower cutoffs revealed a relative impact of dose to each functioning zone to RP. RRCs were 0.52 for the mean dose on the functioning zone, with top 20 % of all functioning zone was two times greater than that of 0.19 for these with 60 %-80 % and 0.17 with 40 %-60 % (P < 0.01).

CONCLUSIONS

The introduction of dose-function features computed from functioning zones with lower and upper cutoffs in a machine learning framework can improve RP prediction. The RRC given by LASSO using dose-function features allows for the quantification of the RP impact of dose on each functioning zones and having the potential to support treatment planning on functional image-guided radiotherapy.

A Transfer Learning Framework for Deep Learning-Based CT-to-Perfusion Mapping on Lung Cancer Patients

Purpose

Deep learning model has shown the feasibility of providing spatial lung perfusion information based on CT images. However, the performance of this method on lung cancer patients is yet to be investigated. This study aims to develop a transfer learning framework to evaluate the deep learning based CT-to-perfusion mapping method specifically on lung cancer patients.

Methods

SPECT/CT perfusion scans of 33 lung cancer patients and 137 non-cancer patients were retrospectively collected from two hospitals. To adapt the deep learning model on lung cancer patients, a transfer learning framework was developed to utilize the features learned from the non-cancer patients. These images were processed to extract features from three-dimensional CT images and synthesize the corresponding CT-based perfusion images. A pre-trained model was first developed using a dataset of patients with lung diseases other than lung cancer, and subsequently fine-tuned specifically on lung cancer patients under three-fold cross-validation. A multi-level evaluation was performed between the CT-based perfusion images and ground-truth SPECT perfusion images in aspects of voxel-wise correlation using Spearman’s correlation coefficient (R), function-wise similarity using Dice Similarity Coefficient (DSC), and lobe-wise agreement using mean perfusion value for each lobe of the lungs.

Results

The fine-tuned model yielded a high voxel-wise correlation (0.8142 ± 0.0669) and outperformed the pre-trained model by approximately 8%. Evaluation of function-wise similarity indicated an average DSC value of 0.8112 ± 0.0484 (range: 0.6460-0.8984) for high-functional lungs and 0.8137 ± 0.0414 (range: 0.6743-0.8902) for low-functional lungs. Among the 33 lung cancer patients, high DSC values of greater than 0.7 were achieved for high functional volumes in 32 patients and low functional volumes in all patients. The correlations of the mean perfusion value on the left upper lobe, left lower lobe, right upper lobe, right middle lobe, and right lower lobe were 0.7314, 0.7134, 0.5108, 0.4765, and 0.7618, respectively.

Conclusion

For lung cancer patients, the CT-based perfusion images synthesized by the transfer learning framework indicated a strong voxel-wise correlation and function-wise similarity with the SPECT perfusion images. This suggests the great potential of the deep learning method in providing regional-based functional information for functional lung avoidance radiation therapy.

Prediction of radiation pneumonitis with machine learning using 4D-CT based dose-function features

In this article, we highlight the fundamental importance of the simultaneous use of dose-volume histogram (DVH) and dose-function histogram (DFH) features based on functional images calculated from 4-dimensional computed tomography (4D-CT) and deformable image registration (DIR) in developing a multivariate radiation pneumonitis (RP) prediction model. The patient characteristics, DVH features and DFH features were calculated from functional images by Hounsfield unit (HU) and Jacobian metrics, for an RP grade ≥ 2 multivariate prediction models were computed from 85 non-small cell lung cancer patients. The prediction model is developed using machine learning via a kernel-based support vector machine (SVM) machine. In the patient cohort, 21 of the 85 patients (24.7%) presented with RP grade ≥ 2. The median area under curve (AUC) was 0.58 for the generated 50 prediction models with patient clinical features and DVH features. When HU metric and Jacobian metric DFH features were added, the AUC improved to 0.73 and 0.68, respectively. We conclude that predictive RP models that incorporate DFH features were successfully developed via kernel-based SVM. These results demonstrate that effectiveness of the simultaneous use of DVH features and DFH features calculated from 4D-CT and DIR on functional image-guided radiotherapy.

Quantifying pulmonary perfusion from noncontrast computed tomography

Purpose

Computed tomography (CT)‐derived ventilation methods compute respiratory induced volume changes as a surrogate for pulmonary ventilation. Currently, there are no known methods to derive perfusion information from noncontrast CT. We introduce a novel CT‐Perfusion (CT‐P) method for computing the magnitude mass changes apparent on dynamic noncontrast CT as a surrogate for pulmonary perfusion.

Methods

CT‐Perfusion is based on a mass conservation model which describes the unknown mass change as a linear combination of spatially corresponding inhale and exhale HU estimated voxel densities. CT‐P requires a deformable image registration (DIR) between the inhale/exhale lung CT pair, a preprocessing lung volume segmentation, and an estimate for the Jacobian of the DIR transformation. Given this information, the CT‐P image, which provides the magnitude mass change for each voxel within the lung volume, is formulated as the solution to a constrained linear least squares problem defined by a series of subregional mean magnitude mass change measurements. Similar to previous robust CT‐ventilation methods, the amount of uncertainty in a subregional sample mean measurement is related to measurement resolution and can be characterized with respect to a tolerance parameter τ. Spatial Spearman correlation between single photon emission CT perfusion (SPECT‐P) and the proposed CT‐P method was assessed in two patient cohorts via a parameter sweep of τ. The first cohort was comprised of 15 patients diagnosed with pulmonary embolism (PE) who had SPECT‐P and 4DCT imaging acquired within 24 h of PE diagnosis. The second cohort was comprised of 15 nonsmall cell lung cancer patients who had SPECT‐P and 4DCT images acquired prior to radiotherapy. For each test case, CT‐P images were computed for 30 different uncertainty parameter values, uniformly sampled from the range [0.01, 0.125], and the Spearman correlation between the SPECT‐P and the resulting CT‐P images were computed.

Results

The median correlations between CT‐P and SPECT‐P taken over all 30 test cases ranged between 0.49 and 0.57 across the parameter sweep. For the optimal tolerance τ = 0.0385, the CT‐P and SPECT‐P correlations across all 30 test cases ranged between 0.02 and 0.82. A one‐sample sign test was applied separately to the PE and lung cancer cohorts. A low Spearmen correlation of 15% was set as the null median value and two‐sided alternative was tested. The PE patients showed a median correlation of 0.57 (IQR = 0.305). One‐sample sign test was statistically significant with 96.5 % confidence interval: 0.20–0.63, P < 0.00001. Lung cancer patients had a median correlation of 0.57(IQR = 0.230). Again, a one‐sample sign test for median was statistically significant with 96.5 percent confidence interval: 0.45–0.71, P < 0.00001.

Conclusion

CT‐Perfusion is the first mechanistic model designed to quantify magnitude blood mass changes on noncontrast dynamic CT as a surrogate for pulmonary perfusion. While the reported correlations with SPECT‐P are promising, further investigation is required to determine the optimal CT acquisition protocol and numerical method implementation for CT‐P imaging.

An assessment of the correlation between robust CT-derived ventilation and pulmonary function test in a cohort with no respiratory symptoms

OBJECTIVE

To evaluate CT-ventilation imaging (CTVI) within a well-characterized, healthy cohort with no respiratory symptoms and examine the correlation between CTVI and concurrent pulmonary function test (PFT).

METHODS

CT scans and PFTs from 77 Caucasian participants in the NORM dataset (clinicaltrials.gov NCT00848406) were analyzed. CTVI was generated using the robust Integrated Jacobian Formulation (IJF) method. IJF estimated total lung capacity (TLC) was computed from CTVI. Bias-adjusted Pearson's correlation between PFT and IJF-based TLC was computed.

RESULTS

IJF- and PFT-measured TLC showed a good correlation for both males and females [males: 0.657, 95% CI (0.438-0.797); females: 0.667, 95% CI (0.416-0.817)]. When adjusting for age, height, smoking, and abnormal CT scan, correlation moderated [males: 0.432, 95% CI (0.129-0.655); females: 0.540, 95% CI (0.207-0.753)]. Visual inspection of CTVI revealed participants who had functional defects, despite the fact that all participant had normal high-resolution CT scan.

CONCLUSION

In this study, we demonstrate that IJF computed CTVI has good correlation with concurrent PFT in a well-validated patient cohort with no respiratory symptoms.

ADVANCES IN KNOWLEDGE

IJF-computed CTVI's overall numerical robustness and consistency with PFT support its potential as a method for providing spatiotemporal assessment of high and low function areas on volumetric non-contrast CT scan.

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.

Variations Between Dose-Ventilation and Dose-Perfusion Metrics in Radiation Therapy Planning for Lung Cancer

Purpose

Currently, several active clinical trials of functional lung avoidance radiation therapy using different imaging modalities for ventilation or perfusion are underway. Patients with lung cancer often show ventilation-perfusion mismatch, whereas the significance of dose-function metric remains unclear. The aim of the present study was to compare dose-ventilation metrics with dose-perfusion metrics for radiation therapy plan evaluation.

Methods and Materials

Pretreatment 4-dimensional computed tomography and ^99mTc-macroaggregated albumin single-photon emission computed tomography perfusion images of 60 patients with lung cancer treated with radiation therapy were analyzed. Ventilation images were created using the deformable image registration of 4-dimensional computed tomography image sets and image analysis for regional volume changes as a surrogate for ventilation. Ventilation and perfusion images were converted into percentile distribution images. Analyses included Pearson’s correlation coefficient and comparison of agreements between the following dose-ventilation and dose-perfusion metrics: functional mean lung dose and functional percent lung function receiving 5, 10, 20, 30, and 40 Gy (fV₅, fV₁₀, fV₂₀, fV₃₀, and fV₄₀, respectively).

Results

Overall, the dose-ventilation metrics were greater than the dose-perfusion metrics (ie, fV₂₀, 26.3% ± 9.9% vs 23.9% ± 9.8%). Correlations between the dose-ventilation and dose-perfusion metrics were strong (range, r = 0.94-0.97), whereas the agreements widely varied among patients, with differences as large as 6.6 Gy for functional mean lung dose and 11.1% for fV₂₀. Paired t test indicated that the dose-ventilation and dose-perfusion metrics were significantly different.

Conclusions

Strong correlations were present between the dose-ventilation and dose-perfusion metrics. However, the agreement between the dose-ventilation and dose-perfusion metrics widely varied among patients, suggesting that ventilation-based radiation therapy plan evaluation may not be comparable to that based on perfusion. Future studies should elucidate the correlation of dose-function metrics with clinical pulmonary toxicity metrics.

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.

A deep learning method for producing ventilation images from 4DCT: First comparison with technegas SPECT ventilation

PURPOSE

The purpose of this study is to develop a deep learning (DL) method for producing four-dimensional computed tomography (4DCT) ventilation imaging and to evaluate the accuracy of the DL-based ventilation imaging against single-photon emission-computed tomography (SPECT) ventilation imaging (SPECT-VI). The performance of the DL-based method is assessed by comparing with density change- and Jacobian-based (HU and JAC) methods.

MATERIALS AND METHODS

Fifty patients with esophagus or lung cancer who underwent thoracic radiotherapy were enrolled in this study. For each patient, 4DCT scans paired with 99mTc-Technegas SPECT/CT were acquired before the first radiotherapy treatment. 4DCT and SPECT/CT were first rigidly registered using MIMvista and converted to data matrix using MATLAB, and then transferred to a DL model based on U-net for correlating 4DCT features and SPECT-VI. Two forms of 4DCT dataset [(a) ten phases and (b) two phases of peak-exhalation and peak-inhalation] as input are studied. Tenfold cross-validation procedure was used to evaluate the performance of the DL model. For comparative evaluation, HU and JAC methodologies are used to calculate specific ventilation imaging based on 4DCT (CTVI) for each patient. The voxel-wise Spearman's correlation was evaluated over the whole lung between each of CTVI and corresponding SPECT-VI. The SPECT-VI and produced CTVIs were segmented into high, median, and low functional lung (HFL, MFL, and LFL) regions. The spatial overlap of corresponding HFL, MFL, and LFL for each CTVI against SPECT-VI was also evaluated using the dice similarity coefficient (DSC). The averaged DSC of functional lung regions was calculated and statistically analyzed with a one-factor ANONA model among different methods.

RESULTS

The voxel-wise Spearman r_s values were (0.22 ± 0.31), (-0.09 ± 0.18), and (0.73 ± 0.16)/(0.71 ± 0.17) for the CTVI_HU , CTVI_JAC , and CTVI_DL(1) /CTVI_DL(2) . These results showed the DL method yielded the strongest correlation with SPECT-VI. Using the DSC as the spatial overlap metric, we found that the CTVI_HU , CTVI_JAC , and CTVI_DL(1) /CTVI_DL(2) methods achieved averaged DSC values for all patients to be (0.45 ± 0.08), (0.33 ± 0.04), and (0.73 ± 0.09)/(0.71 ± 0.09), respectively. The results demonstrated that the DL method yielded the highest similarity with SPECT-VI with the prominently significant difference (P < 10^-7 ).

CONCLUSIONS

This study developed a DL method for producing CTVI and performed a validation against SPECT-VI. The results demonstrated that DL method can derive CTVI with greatly improved accuracy in comparison to HU and JAC methods. The produced ventilation images can be more accurate and useful for lung functional avoidance radiotherapy and treatment response modeling.

Robust HU-based CT ventilation from an integrated mass conservation formulation

Computed tomography (CT) ventilation algorithms estimate volume changes induced by respiratory motion. Existing Hounsfield Unit (HU) methods approximate volume change from the measured HU variations between spatially corresponding voxel locations within a temporally resolved CT image pair, assuming that volume changes are caused solely by changes in air content. Numerical implementations require a deformable image registration to determine the inhale/exhale spatial correspondence, a preprocessing lung volume segmentation, a preprocessing high-intensity vessel segmentation, and a post-processing smoothing applied to the raw volume change estimates obtained for each lung tissue voxel.

PURPOSE

We introduce the novel mass-conserving volume change (MCVC) method for estimating voxel volume changes from the HU values within an inhale/exhale CT image pair. MCVC is based on subregional volume change estimates that possess quantitatively characterized and controllable levels of uncertainty. MCVC is therefore robust to small variations in DIR solutions and the resulting ventilation images are overall more reproducible. In contrast to existing HU methods, MCVC does not require a preprocessing lung vessel segmentation or pre/post-processing Gaussian smoothing.

METHODS

Subregional volume change estimates are defined in terms of mean density ratios. As such, the corresponding uncertainty is characterized using Gaussian statistics and standard error analysis of the sample density means. A numerical solution is obtained from the MCVC formulation by solving a constrained linear least squares problem defined by a series of subregional volume change estimates. Reproducibility of the MCVC method with respect to DIR solution was assessed using expert-determined landmark point pairs and inhale/exhale phases from 10 four-dimensional CTs (4DCTs) available on www.dir-lab.com. MCVC was also compared to the robust Integrated Jacobian Formulation (IJF), a transformation-based ventilation method.

RESULTS

The ten Dir-Lab 4DCT cases were registered twice with the same DIR algorithm, but using different degrees of freedom (DIR 1 and DIR 2). Standard HU ventilation (HUV) and MCVC ventilation images were computed for both solutions. The average spatial errors (300 landmarks per case) for DIR 1 ranged between 0.74 and 1.50 mm, whereas for DIR 2 they ranged between 0.68 and 1.18 mm. Despite the differences in spatial errors between the two DIR solutions, the average Pearson correlation between the corresponding HUV images was 0.94 (0.03), while for the MCVC images it was 1.00 (0.00). The average correlation between MCVC and IJF ventilation over the ten test cases was 0.81 (0.14), whereas for HUV and IJF it was 0.56 (1.11).

CONCLUSION

While HUV is robust to DIR solution, its implementation depends on heuristic Gaussian smoothing and vessel segmentation. MCVC is based on subregional volume change measurements with quantifiable and controllable levels of uncertainty. The subregional approach eliminates the need for Gaussian smoothing and lung vasculature segmentation. Numerical experiments are consistent with the underlying mathematical theory and indicate that MCVC ventilation is more reproducible with respect to DIR algorithm than standard HU methods. MCVC results are also more consistent with the robust IJF method, which suggests that incorporating robustness leads to more consistent results across both DIRs and ventilation algorithms.

Deformable image registration for radiation therapy: principle, methods, applications and evaluation

Background: Deformable image registration (DIR) is increasingly used in the field of radiation therapy (RT) to account for anatomical deformations. The aims of this paper are to describe the main applications of DIR in RT and discuss current DIR evaluation methods. Methods: Articles on DIR published from January 2000 to October 2018 were extracted from PubMed and Science Direct. Our search was restricted to articles that report data obtained from humans, were written in English, and address DIR methods for RT. A total of 207 articles were selected from among 2506 identified in the search process. Results: At planning, DIR is used for organ delineation using atlas-based segmentation, deformation-based planning target volume definition, functional planning and magnetic resonance imaging-based dose calculation. In image-guided RT, DIR is used for contour propagation and dose calculation on per-treatment imaging. DIR is also used to determine the accumulated dose from fraction to fraction in external beam RT and brachytherapy, both for dose reporting and adaptive RT. In the case of re-irradiation, DIR can be used to estimate the cumulated dose of the two irradiations. Finally, DIR can be used to predict toxicity in voxel-wise population analysis. However, the evaluation of DIR remains an open issue, especially when dealing with complex cases such as the disappearance of matter. To quantify DIR uncertainties, most evaluation methods are limited to geometry-based metrics. Software companies have now integrated DIR tools into treatment planning systems for clinical use, such as contour propagation and fraction dose accumulation. Conclusions: DIR is increasingly important in RT applications, from planning to toxicity prediction. DIR is routinely used to reduce the workload of contour propagation. However, its use for complex dosimetric applications must be carefully evaluated by combining quantitative and qualitative analyses.

Patch-based lung ventilation estimation using multi-layer supervoxels

Patch-based approaches have received substantial attention over the recent years in medical imaging. One of their potential applications may be to provide more anatomically consistent ventilation maps estimated on dynamic lung CT. An assessment of regional lung function may act as a guide for radiotherapy, ensuring a more accurate treatment plan. This in turn, could spare well-functioning parts of the lungs. We present a novel method for lung ventilation estimation from dynamic lung CT imaging, combining a supervoxel-based image representation with deformations estimated during deformable image registration, performed between peak breathing phases. For this we propose a method that tracks changes of the intensity of previously extracted supervoxels. For the evaluation of the method we calculate correlation of the estimated ventilation maps with static ventilation images acquired from hyperpolarized Xenon129 MRI. We also investigate the influence of different image registration methods used to estimate deformations between the peak breathing phases in the dynamic CT imaging. We show that our method performs favorably to other ventilation estimation methods commonly used in the field, independently of the image registration method applied to dynamic CT. Due to the patch-based approach of our method, it may be physiologically more consistent with lung anatomy than previous methods relying on voxel-wise relationships. In our method the ventilation is estimated for supervoxels, which tend to group spatially close voxels with similar intensity values. The proposed method was evaluated on a dataset consisting of three lung cancer patients undergoing radiotherapy treatment, and this resulted in a correlation of 0.485 with XeMRI ventilation images, compared with 0.393 for the intensity-based approach, 0.231 for the Jacobian-based method and 0.386 for the Hounsfield units averaging method, on average. Within the limitation of the small number of cases analyzed, results suggest that the presented technique may be advantageous for CT-based ventilation estimation. The results showing higher values of correlation of the proposed method demonstrate the potential of our method to more accurately mimic the lung physiology.

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.

Robust CT ventilation from the integral formulation of the Jacobian

Computed tomography (CT) derived ventilation algorithms estimate the apparent voxel volume changes within an inhale/exhale CT image pair. Transformation‐based methods compute these estimates solely from the spatial transformation acquired by applying a deformable image registration (DIR) algorithm to the image pair. However, approaches based on finite difference approximations of the transformation's Jacobian have been shown to be numerically unstable. As a result, transformation‐based CT ventilation is poorly reproducible with respect to both DIR algorithm and CT acquisition method.

Purpose

We introduce a novel Integrated Jacobian Formulation (IJF) method for estimating voxel volume changes under a DIR‐recovered spatial transformation. The method is based on computing volume estimates of DIR mapped subregions using the hit‐or‐miss sampling algorithm for integral approximation. The novel approach allows for regional volume change estimates that (a) respect the resolution of the digital grid and (b) are based on approximations with quantitatively characterized and controllable levels of uncertainty. As such, the IJF method is designed to be robust to variations in DIR solutions and thus overall more reproducible.

Methods

Numerically, Jacobian estimates are recovered by solving a simple constrained linear least squares problem that guarantees the recovered global volume change is equal to the global volume change obtained from the inhale and exhale lung segmentation masks. Reproducibility of the IJF method with respect to DIR solution was assessed using the expert‐determined landmark point pairs and inhale/exhale phases from 10 four‐dimensional computed tomographies (4DCTs) available on http://www.dir-lab.com. Reproducibility with respect to CT acquisition was assessed on the 4DCT and 4D cone beam CT (4DCBCT) images acquired for five lung cancer patients prior to radiotherapy.

Results

The ten Dir‐Lab 4DCT cases were registered twice with the same DIR algorithm, but with different smoothing parameter. Finite difference Jacobian (FDJ) and IFJ images were computed for both solutions. The average spatial errors (300 landmarks per case) for the two DIR solution methods were 0.98 (1.10) and 1.02 (1.11). The average Pearson correlation between the FDJ images computed from the two DIR solutions was 0.83 (0.03), while for the IJF images it was 1.00 (0.00). For intermodality assessment, the IJF and FDJ images were computed from the 4DCT and 4DCBCT of five patients. The average Pearson correlation of the spatially aligned FDJ images was 0.27 (0.11), while it was 0.77 (0.13) for the IFJ method.

Conclusion

The mathematical theory underpinning the IJF method allows for the generation of ventilation images that are (a) computed with respect to DIR spatial accuracy on the digital voxel grid and (b) based on DIR‐measured subregional volume change estimates acquired with quantifiable and controllable levels of uncertainty. Analyses of the experiments are consistent with the mathematical theory and indicate that IJF ventilation imaging has a higher reproducibility with respect to both DIR algorithm and CT acquisition method, in comparison to the standard finite difference approach.

Comparison of CT ventilation imaging and hyperpolarised gas MRI: effects of breathing manoeuvre

Image registration of lung CT images acquired at different inflation levels has been proposed as a surrogate method to map lung 'ventilation'. Prior to clinical use, it is important to understand how this technique compares with direct ventilation imaging modalities such as hyperpolarised gas MRI. However, variations in lung inflation level have been shown to affect regional ventilation distributions. Therefore, the aim of this study was to evaluate the impact of lung inflation levels when comparing CT ventilation imaging to ventilation from ³He-MRI. Seven asthma patients underwent breath-hold CT at total lung capacity (TLC) and functional residual capacity (FRC). ³He-MRI and a same-breath ¹H-MRI were acquired at FRC+1L and TLC. Percentage ventilated volumes (%VVs) were calculated for FRC+1L and TLC ³He-MRI. TLC-CT and registered FRC-CT were used to compute a surrogate ventilation map from voxel-wise intensity differences in Hounsfield unit values, which was thresholded at the 10th and 20th percentiles. For direct comparison of CT and ³He-MRI ventilation, FRC+1L and TLC ³He-MRI were registered to TLC-CT indirectly via the corresponding same-breath ¹H-MRI data. For ³He-MRI and CT ventilation comparison, Dice similarity coefficients (DSCs) between the binary segmentations were computed. The median (range) of %VVs for FRC+1L and TLC ³He-MRI were 90.5 (54.9-93.6) and 91.8 (67.8-96.2), respectively (p   =  0.018). For MRI versus CT ventilation comparison, statistically significant improvements in DSCs were observed for TLC ³He MRI when compared with FRC+1L, with median (range) values of 0.93 (0.86-0.93) and 0.86 (0.68-0.92), respectively (p   =  0.017), for the 10-100th percentile and 0.87 (0.83-0.88) and 0.81 (0.66-0.87), respectively (p   =  0.027), for the 20-100th percentile. Correlation of CT ventilation imaging and hyperpolarised gas MRI is sensitive to lung inflation level. For ventilation maps derived from CT acquired at FRC and TLC, a higher correlation with gas ventilation MRI can be achieved if the MRI is acquired at TLC.

The VAMPIRE challenge: A multi-institutional validation study of CT ventilation imaging

PURPOSE

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

METHODS

The 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 r S , and Dice similarity coefficients evaluated for low function lung ( DSC low ) and high function lung ( DSC high ).

RESULTS

A 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 r S 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 r S , 0.52 (0.36-0.67) for DSC low , and 0.45 (0.28-0.62) for DSC high . All other algorithms exhibited at least one negative r S value, and/or one DSC value less than 0.5.

CONCLUSIONS

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

An initial investigation of hyperpolarized gas tagging magnetic resonance imaging in evaluating deformable image registration-based lung ventilation

BACKGROUND

Deformable image registration (DIR)-based lung ventilation mapping is attractive due to its simplicity, and also challenging due to its susceptibility to errors and uncertainties. In this study, we explored the use of 3D Hyperpolarized (HP) gas tagging MRI to evaluate DIR-based lung ventilation.

METHOD AND MATERIAL

Three healthy volunteers included in this study underwent both 3D HP gas tagging MRI (t-MRI) and 3D proton MRI (p-MRI) using balanced steady-state free precession pulse sequence at end of inhalation and end of exhalation. We first obtained the reference displacement vector fields (DVFs) from the t-MRIs by tracking the motion of each tagging grid between the exhalation and the inhalation phases. Then, we determined DIR-based DVFs from the p-MRIs by registering the images at the two phases with two commercial DIR algorithms. Lung ventilations were calculated from both the reference DVFs and the DIR-based DVFs using the Jacobian method and then compared using cross correlation and mutual information.

RESULTS

The DIR-based lung ventilations calculated using p-MRI varied considerably from the reference lung ventilations based on t-MRI among all three subjects. The lung ventilations generated using Velocity AI were preferable for the better spatial homogeneity and accuracy compared to the ones using MIM, with higher average cross correlation (0.328 vs 0.262) and larger average mutual information (0.528 vs 0.323).

CONCLUSION

We demonstrated that different DIR algorithms resulted in different lung ventilation maps due to underlining differences in the DVFs. HP gas tagging MRI provides a unique platform for evaluating DIR-based lung ventilation.

Interim Analysis of a Two-Institution, Prospective Clinical Trial of 4DCT-Ventilation-based Functional Avoidance Radiation Therapy

PURPOSE

Functional imaging has been proposed that uses 4DCT images to calculate 4DCT-based lung ventilation (4DCT-ventilation). We have started a 2-institution, phase 2 prospective trial evaluating the feasibility, safety, and preliminary efficacy of 4DCT-ventilation functional avoidance. The trial hypothesis is that the rate of grade ≥2 radiation pneumonitis could be reduced to 12% with functional avoidance, compared with a 25% rate of pneumonitis with a historical control. The trial employed a Simon 2-stage design with a planned futility analysis after 17 evaluable patients. The purpose of this work is to present the trial design and implementation, dosimetric data, and clinical results for the planned futility analysis.

METHODS AND MATERIALS

Eligible patients were patients with lung cancer who were prescribed doses of 45 to 75 Gy. For each patient, the 4DCT data were used to generate a 4DCT-ventilation image using the Hounsfield unit technique along with a compressible flow-based image registration algorithm. Two intensity modulated radiation therapy treatment plans were generated: (1) a standard lung plan and (2) a functional avoidance treatment plan that aimed to reduce dose to functional lung while meeting target and normal tissue constraints. Patients were treated with the functional avoidance plan and evaluated for thoracic toxicity (presented as rate and 95% confidence intervals [CI]) with a 1-year follow-up.

RESULTS

The V20 to functional lung was 21.6% ± 9.5% (mean ± standard deviation) with functional avoidance, representing a decrease of 3.2% (P < .01) relative to standard, nonfunctional treatment plans. The rates of grade ≥2 and grade ≥3 radiation pneumonitis were 17.6% (95% CI, 3.8%-43.4%) and 5.9% (95% CI, 0.1%-28.7%), respectively.

CONCLUSIONS

Dosimetrically, functional avoidance achieved reduction in doses to functional lung while meeting target and organ at risk constraints. On the basis of Simon's 2-stage design and the 17.6% grade ≥2 pneumonitis rate, the trial met its futility criteria and has continued accrual.

Evaluation of lung toxicity risk with computed tomography ventilation image for thoracic cancer patients

Background

Four-dimensional computed tomography (4D-CT) ventilation is an emerging imaging modality. Functional avoidance of regions according to 4D-CT ventilation may reduce lung toxicity after radiation therapy. This study evaluated associations between 4D-CT ventilation-based dosimetric parameters and clinical outcomes.

Methods

Pre-treatment 4D-CT data were used to retrospectively construct ventilation images for 40 thoracic cancer patients retrospectively. Fifteen patients were treated with conventional radiation therapy, 6 patients with hyperfractionated radiation therapy and 19 patients with stereotactic body radiation therapy (SBRT). Ventilation images were calculated from 4D-CT data using a deformable image registration and Jacobian-based algorithm. Each ventilation map was normalized by converting it to percentile images. Ventilation-based dosimetric parameters (Mean Dose, V5 [percent lung volume receiving ≥5 Gy], and V20 [percent lung volume receiving ≥20 Gy]) were calculated for highly and poorly ventilated regions. To test whether the ventilation-based dosimetric parameters could be used predict radiation pneumonitis of ≥Grade 2, the area under the curve (AUC) was determined from the receiver operating characteristic analysis.

Results

For Mean Dose, poorly ventilated lung regions in the 0–30% range showed the highest AUC value (0.809; 95% confidence interval [CI], 0.663–0.955). For V20, poorly ventilated lung regions in the 0–20% range had the highest AUC value (0.774; 95% [CI], 0.598–0.915), and for V5, poorly ventilated lung regions in the 0–30% range had the highest AUC value (0.843; 95% [CI], 0.732–0.954). The highest AUC values for Mean Dose, V20, and V5 were obtained in poorly ventilated regions. There were significant differences in all dosimetric parameters between radiation pneumonitis of Grade 1 and Grade ≥2.

Conclusions

Poorly ventilated lung regions identified on 4D-CT had higher AUC values than highly ventilated regions, suggesting that functional planning based on poorly ventilated regions may reduce the risk of lung toxicity in radiation therapy.

Correlation of Regional Lung Ventilation and Gas Transfer to Red Blood Cells: Implications for Functional-Avoidance Radiation Therapy Planning

PURPOSE

To investigate the degree to which lung ventilation and gas exchange are regionally correlated, using the emerging technology of hyperpolarized (HP)-129Xe magnetic resonance imaging (MRI).

METHODS AND MATERIALS

Hyperpolarized-129Xe MRI studies were performed on 17 institutional review board-approved human subjects, including 13 healthy volunteers, 1 emphysema patient, and 3 non-small cell lung cancer patients imaged before and approximately 11 weeks after radiation therapy (RT). Subjects inhaled 1 L of HP-129Xe mixture, followed by the acquisition of interleaved ventilation and gas exchange images, from which maps were obtained of the relative HP-129Xe distribution in three states: (1) gaseous, in lung airspaces; (2) dissolved interstitially, in alveolar barrier tissue; and (3) transferred to red blood cells (RBCs), in the capillary vasculature. The relative spatial distributions of HP-129Xe in airspaces (regional ventilation) and RBCs (regional gas transfer) were compared. Further, we investigated the degree to which ventilation and RBC transfer images identified similar functional regions of interest (ROIs) suitable for functionally guided RT. For the RT patients, both ventilation and RBC functional images were used to calculate differences in the lung dose-function histogram and functional effective uniform dose.

RESULTS

The correlation of ventilation and RBC transfer was ρ = 0.39 ± 0.15 in healthy volunteers. For the RT patients, this correlation was ρ = 0.53 ± 0.02 before treatment and ρ = 0.39 ± 0.07 after treatment; for the emphysema patient it was ρ = 0.24. Comparing functional ROIs, ventilation and RBC transfer demonstrated poor spatial agreement: Dice similarity coefficient = 0.50 ± 0.07 and 0.26 ± 0.12 for the highest-33%- and highest-10%-function ROIs in healthy volunteers, and in RT patients (before treatment) these were 0.58 ± 0.04 and 0.40 ± 0.04. The average magnitude of the differences between RBC- and ventilation-derived functional effective uniform dose, fV20Gy, fV10Gy, and fV5Gy were 1.5 ± 1.4 Gy, 4.1% ± 3.8%, 5.0% ± 3.8%, and 5.3% ± 3.9%, respectively.

CONCLUSION

Ventilation may not be an effective surrogate for true regional lung function for all patients.

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.

Dosimetric feasibility of 4DCT-ventilation imaging guided proton therapy for locally advanced non-small-cell lung cancer

Background

The principle aim of this study is to incorporate 4DCT ventilation imaging into functional treatment planning that preserves high-functioning lung with both double scattering and scanning beam techniques in proton therapy.

Methods

Eight patients with locally advanced non-small-cell lung cancer were included in this study. Deformable image registration was performed for each patient on their planning 4DCTs and the resultant displacement vector field with Jacobian analysis was used to identify the high-, medium- and low-functional lung regions. Five plans were designed for each patient: a regular photon IMRT vs. anatomic proton plans without consideration of functional ventilation information using double scattering proton therapy (DSPT) and intensity modulated proton therapy (IMPT) vs. functional proton plans with avoidance of high-functional lung using both DSPT and IMPT. Dosimetric parameters were compared in terms of tumor coverage, plan heterogeneity, and avoidance of normal tissues.

Results

Our results showed that both DSPT and IMPT plans gave superior dose advantage to photon IMRTs in sparing low dose regions of the total lung in terms of V5 (volume receiving 5Gy). The functional DSPT only showed marginal benefit in sparing high-functioning lung in terms of V5 or V20 (volume receiving 20Gy) compared to anatomical plans. Yet, the functional planning in IMPT delivery, can further reduce the low dose in high-functioning lung without degrading the PTV dosimetric coverages, compared to anatomical proton planning. Although the doses to some critical organs might increase during functional planning, the necessary constraints were all met.

Conclusions

Incorporating 4DCT ventilation imaging into functional proton therapy is feasible. The functional proton plans, in intensity modulated proton delivery, are effective to further preserve high-functioning lung regions without degrading the PTV coverage.

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.

Regional Lung Function Profiles of Stage I and III Lung Cancer Patients: An Evaluation for Functional Avoidance Radiation Therapy

PURPOSE

The development of clinical trials is underway to use 4-dimensional computed tomography (4DCT) ventilation imaging to preferentially spare functional lung in patients undergoing radiation therapy. The purpose of this work was to generate data to aide with clinical trial design by retrospectively characterizing dosimetric and functional profiles for patients with different stages of lung cancer.

METHODS AND MATERIALS

A total of 118 lung cancer patients (36% stage I and 64% stage III) from 2 institutions were used for the study. A 4DCT-ventilation map was calculated using the patient's 4DCT imaging, deformable image registration, and a density-change-based algorithm. To assess each patient's spatial ventilation profile both quantitative and qualitative metrics were developed, including an observer-based defect observation and metrics based on the ventilation in each lung third. For each patient we used the clinical doses to calculate functionally weighted mean lung doses and metrics that assessed the interplay between the spatial location of the dose and high-functioning lung.

RESULTS

Both qualitative and quantitative metrics revealed a significant difference in functional profiles between the 2 stage groups (P<.01). We determined that 65% of stage III and 28% of stage I patients had ventilation defects. Average functionally weighted mean lung dose was 19.6 Gy and 5.4 Gy for stage III and I patients, respectively, with both groups containing patients with large spatial overlap between dose and high-function regions.

CONCLUSION

Our 118-patient retrospective study found that 65% of stage III patients have regionally variant ventilation profiles that are suitable for functional avoidance. Our results suggest that regardless of disease stage, it is possible to have unique spatial interplay between dose and high-functional lung, highlighting the importance of evaluating the function of each patient and developing a personalized functional avoidance treatment approach.

A complete 4DCT-ventilation functional avoidance virtual trial: Developing strategies for prospective clinical trials

Introduction

4DCT‐ventilation is an exciting new imaging modality that uses 4DCT data to calculate lung‐function maps. Because 4DCTs are acquired as standard of care for lung cancer patients undergoing radiotherapy, 4DCT‐ventiltation provides functional information at no extra dosimetric or monetary cost to the patient. The development of clinical trials is underway to use 4DCT‐ventilation imaging to spare functional lung in patients undergoing radiotherapy. The purpose of this work was to perform a virtual trial using retrospective data to develop the practical aspects of a 4DCT‐ventilation functional avoidance clinical trial.

Methods

The study included 96 stage III lung cancer patients. A 4DCT‐ventilation map was calculated using the patient's 4DCT‐imaging, deformable registration, and a density‐change‐based algorithm. Clinical trial inclusion assessment used quantitative and qualitative metrics based on the patient's spatial ventilation profile. Clinical and functional plans were generated for 25 patients. The functional plan aimed to reduce dose to functional lung while meeting standard target and critical structure constraints. Standard and dose‐function metrics were compared between the clinical and functional plans.

Results

Our data showed that 69% and 59% of stage III patients have regional variability in function based on qualitative and quantitative metrics, respectively. Functional planning demonstrated an average reduction of 2.8 Gy (maximum 8.2 Gy) in the mean dose to functional lung.

Conclusions

Our work demonstrated that 60–70% of stage III patients would be eligible for functional planning and that a typical functional lung mean dose reduction of 2.8 Gy can be expected relative to standard clinical plans. These findings provide salient data for the development of functional clinical trials.

Assessing the use of 4DCT-ventilation in pre-operative surgical lung cancer evaluation

PURPOSE

A primary treatment option for lung cancer patients is surgical resection. Patients who have poor lung function prior to surgery are at increased risk of developing serious and life-threatening complications after surgical resection. Surgeons use nuclear medicine ventilation-perfusion (VQ) scans along with pulmonary function test (PFT) information to assess a patient's pre-surgical lung function. The nuclear medicine images and pre-surgery PFTs are used to calculate percent predicted postoperative (%PPO) PFT values by estimating the amount of functioning lung tissue that would be lost with surgical resection. Nuclear medicine imaging is currently considered the standard of care when evaluating the amount of ventilation that would be lost due to surgery. A novel lung function imaging modality has been developed in radiation oncology that uses 4-Dimensional computed tomography data to calculate ventilation maps (4DCT-ventilation). Compared to nuclear medicine, 4DCT-ventilation is cheaper, does not require a radioactive contrast agent, provides a faster imaging procedure, and has improved spatial resolution. In this work we perform a retrospective study to assess the use of 4DCT-ventilation as a pre-operative surgical lung function evaluation tool. Specifically, the purpose of our study was to compare %PPO PFT values calculated with 4DCT-ventilation and %PPO PFT values calculated with nuclear medicine ventilation-perfusion imaging.

METHODS

The study included 16 lung cancer patients that had undergone 4DCT imaging, nuclear medicine imaging, and had Forced Expiratory Volume in 1 second (FEV₁ ) acquired as part of a standard PFT. The 4DCT datasets, spatial registration, and a density-change-based model were used to compute 4DCT-ventilation maps. Both 4DCT-ventilation and nuclear medicine images were used to calculate %PPO FEV₁ . The %PPO FEV₁ was calculated by scaling the pre-surgical FEV₁ by (1-fraction of total resected ventilation); where the resected ventilation was determined using either the 4DCT-ventilation or nuclear medicine imaging. Calculations were done assuming both lobectomy and pneumonectomy resections. The %PPO FEV₁ values were compared between the 4DCT-ventilation-based calculations and the nuclear medicine-based calculations using correlation coefficients, average differences, and Receiver Operating Characteristic (ROC) analysis.

RESULTS

Overall the 4DCT-ventilation derived %PPO FEV₁ values agreed well with nuclear medicine-derived %PPO FEV₁ data with correlations of 0.99 and 0.81 for lobectomy and pneumonectomy, respectively. The average differences between the 4DCT-ventilation and nuclear medicine-based calculation for %PPO FEV₁ were less than 5%. ROC analysis revealed predictive accuracy that ranged from 87.5% to 100% when assessing the ability of 4DCT-ventilation to predict for nuclear medicine-based %PPO FEV₁ values.

CONCLUSIONS

4DCT-ventilation is an innovative technology developed in radiation oncology that has great potential to translate to the surgical domain. The high correlation results when comparing 4DCT-ventilation to the current standard of care provide a strong rationale for a prospective clinical trial assessing 4DCT-ventilation in the clinical setting. 4DCT-ventilation can reduce the cost and imaging time for patients while providing improved spatial accuracy and quantitative results for surgeons.

The numerical stability of transformation-based CT ventilation

Abstract

Computed tomography (CT)-derived ventilation imaging utilizes deformable image registration (DIR) to recover respiratory-induced tissue volume changes from inhale/exhale 4DCT phases. While current strategies for validating CT ventilation rely on analyzing its correlation with existing functional imaging modalities, the numerical stability of the CT ventilation calculation has not been characterized.

Purpose

The purpose of this study is to examine how small changes in the DIR displacement field can affect the calculation of transformation-based CT ventilation.

Methods

First, we derive a mathematical theorem, which states that the change in ventilation metric induced by a perturbation to single displacement vector is bounded by the perturbation magnitude. Second, we introduce a novel Jacobian constrained optimization method for computing user-defined CT ventilation images.

Results

Using the Jacobian constrained method, we demonstrate that for the same inhale/exhale CT pair, it is possible to compute two DIR transformations that have similar spatial accuracies, but generate ventilation images with significantly different physical characteristics. In particular, we compute a CT ventilation image that perfectly correlates with a single-photon emission CT perfusion scan.

Conclusion

The analysis and experiments indicate that while transformation-based CT ventilation is a promising modality, small changes in the DIR displacement field can result in large relative changes in the ventilation image. As such, approaches for improving the reproducibility of CT ventilation are still needed.

Evaluation of the ΔV 4D CT ventilation calculation method using in vivo xenon CT ventilation data and comparison to other methods

Ventilation distribution calculation using 4D CT has shown promising potential in several clinical applications. This study evaluated the direct geometric ventilation calculation method, namely the ΔV method, with xenon-enhanced CT (XeCT) ventilation data from four sheep, and compared it with two other published meth-ods, the Jacobian and the Hounsfield unit (HU) methods. Spearman correlation coefficient (SCC) and Dice similarity coefficient (DSC) were used for the evaluation and comparison. The average SCC with one standard deviation was 0.44 ± 0.13 with a range between 0.29 and 0.61 between the XeCT and ΔV ventilation distributions. The average DSC value for lower 30% ventilation volumes between the XeCT and ΔV ventilation distributions was 0.55 ± 0.07 with a range between 0.48 and 0.63. Ventilation difference introduced by deformable image registration errors improved with smoothing. In conclusion, ventilation distributions generated using ΔV-4D CT and deformable image registration are in reasonably agreement with the in vivo XeCT measured ventilation distribution.

Comparison of CT-based Lobar Ventilation with 3He MR Imaging Ventilation Measurements

PURPOSE

To compare lobar lung ventilation computed from expiratory and inspiratory computed tomographic (CT) data with direct measurements of ventilation at hyperpolarized helium 3 ((3)He) magnetic resonance (MR) imaging by using same-breath hydrogen 1 ((1)H) MR imaging examinations to coregister the multimodality images.

MATERIALS AND METHODS

The study was approved by the national research ethics committee, and written patient consent was obtained. Thirty patients with asthma underwent breath-hold CT at total lung capacity and functional residual capacity. (3)He and (1)H MR images were acquired during the same breath hold at a lung volume of functional residual capacity plus 1 L. Lobar segmentations delineated by major fissures on both CT scans were used to calculate the percentage of ventilation per lobe from the change in inspiratory and expiratory lobar volumes. CT-based ventilation was compared with (3)He MR imaging ventilation by using diffeomorphic image registration of (1)H MR imaging to CT, which enabled indirect registration of (3)He MR imaging to CT. Statistical analysis was performed by using the Wilcoxon signed-rank test, Pearson correlation coefficient, and Bland-Altman analysis.

RESULTS

The mean ± standard deviation absolute difference between the CT and (3)He MR imaging percentage of ventilation volume in all lobes was 4.0% (right upper and right middle lobes, 5.4% ± 3.3; right lower lobe, 3.7% ± 3.9; left upper lobe, 2.8% ± 2.7; left lower lobe, 3.9% ± 2.6; Wilcoxon signed-rank test, P < .05). The Pearson correlation coefficient between the two techniques in all lobes was 0.65 (P < .001). Greater percentage of ventilation was seen in the upper lobes with (3)He MR imaging and in the lower lobes with CT. This was confirmed with Bland-Altman analysis, with 95% limits of agreement for right upper and middle lobes, -2.4, 12.7; right lower lobe, -11.7, 4.6; left upper lobe, -4.9, 8.7; and left lower lobe, -9.8, 2.8.

CONCLUSION

The percentage of regional ventilation per lobe calculated at CT was comparable to a direct measurement of lung ventilation at hyperpolarized (3)He MR imaging. This work provides evidence for the validity of the CT model, and same-breath (1)H MR imaging enables regional interpretation of (3)He ventilation MR imaging on the underlying lung anatomy at thin-section CT.

Clinical validation of 4-dimensional computed tomography ventilation with pulmonary function test data

PURPOSE

A new form of functional imaging has been proposed in the form of 4-dimensional computed tomography (4DCT) ventilation. Because 4DCTs are acquired as part of routine care for lung cancer patients, calculating ventilation maps from 4DCTs provides spatial lung function information without added dosimetric or monetary cost to the patient. Before 4DCT-ventilation is implemented it needs to be clinically validated. Pulmonary function tests (PFTs) provide a clinically established way of evaluating lung function. The purpose of our work was to perform a clinical validation by comparing 4DCT-ventilation metrics with PFT data.

METHODS AND MATERIALS

Ninety-eight lung cancer patients with pretreatment 4DCT and PFT data were included in the study. Pulmonary function test metrics used to diagnose obstructive lung disease were recorded: forced expiratory volume in 1 second (FEV1) and FEV1/forced vital capacity. Four-dimensional CT data sets and spatial registration were used to compute 4DCT-ventilation images using a density change-based and a Jacobian-based model. The ventilation maps were reduced to single metrics intended to reflect the degree of ventilation obstruction. Specifically, we computed the coefficient of variation (SD/mean), ventilation V20 (volume of lung ≤20% ventilation), and correlated the ventilation metrics with PFT data. Regression analysis was used to determine whether 4DCT ventilation data could predict for normal versus abnormal lung function using PFT thresholds.

RESULTS

Correlation coefficients comparing 4DCT-ventilation with PFT data ranged from 0.63 to 0.72, with the best agreement between FEV1 and coefficient of variation. Four-dimensional CT ventilation metrics were able to significantly delineate between clinically normal versus abnormal PFT results.

CONCLUSIONS

Validation of 4DCT ventilation with clinically relevant metrics is essential. We demonstrate good global agreement between PFTs and 4DCT-ventilation, indicating that 4DCT-ventilation provides a reliable assessment of lung function. Four-dimensional CT ventilation enables exciting opportunities to assess lung function and create functional avoidance radiation therapy plans. The present work provides supporting evidence for the integration of 4DCT-ventilation into clinical trials.

A Voxel-by-Voxel Comparison of Deformable Vector Fields Obtained by Three Deformable Image Registration Algorithms Applied to 4DCT Lung Studies

BACKGROUND

Commonly used methods of assessing the accuracy of deformable image registration (DIR) rely on image segmentation or landmark selection. These methods are very labor intensive and thus limited to relatively small number of image pairs. The direct voxel-by-voxel comparison can be automated to examine fluctuations in DIR quality on a long series of image pairs.

METHODS

A voxel-by-voxel comparison of three DIR algorithms applied to lung patients is presented. Registrations are compared by comparing volume histograms formed both with individual DIR maps and with a voxel-by-voxel subtraction of the two maps. When two DIR maps agree one concludes that both maps are interchangeable in treatment planning applications, though one cannot conclude that either one agrees with the ground truth. If two DIR maps significantly disagree one concludes that at least one of the maps deviates from the ground truth. We use the method to compare 3 DIR algorithms applied to peak inhale-peak exhale registrations of 4DFBCT data obtained from 13 patients.

RESULTS

All three algorithms appear to be nearly equivalent when compared using DICE similarity coefficients. A comparison based on Jacobian volume histograms shows that all three algorithms measure changes in total volume of the lungs with reasonable accuracy, but show large differences in the variance of Jacobian distribution on contoured structures. Analysis of voxel-by-voxel subtraction of DIR maps shows differences between algorithms that exceed a centimeter for some registrations.

CONCLUSION

Deformation maps produced by DIR algorithms must be treated as mathematical approximations of physical tissue deformation that are not self-consistent and may thus be useful only in applications for which they have been specifically validated. The three algorithms tested in this work perform fairly robustly for the task of contour propagation, but produce potentially unreliable results for the task of DVH accumulation or measurement of local volume change. Performance of DIR algorithms varies significantly from one image pair to the next hence validation efforts, which are exhaustive but performed on a small number of image pairs may not reflect the performance of the same algorithm in practical clinical situations. Such efforts should be supplemented by validation based on a longer series of images of clinical quality.

Evaluation of accuracy of B-spline transformation-based deformable image registration with different parameter settings for thoracic images

Deformable image registration (DIR) is fundamental technique for adaptive radiotherapy and image-guided radiotherapy. However, further improvement of DIR is still needed. We evaluated the accuracy of B-spline transformation-based DIR implemented in elastix. This registration package is largely based on the Insight Segmentation and Registration Toolkit (ITK), and several new functions were implemented to achieve high DIR accuracy. The purpose of this study was to clarify whether new functions implemented in elastix are useful for improving DIR accuracy. Thoracic 4D computed tomography images of ten patients with esophageal or lung cancer were studied. Datasets for these patients were provided by DIR-lab (dir-lab.com) and included a coordinate list of anatomical landmarks that had been manually identified. DIR between peak-inhale and peak-exhale images was performed with four types of parameter settings. The first one represents original ITK (Parameter 1). The second employs the new function of elastix (Parameter 2), and the third was created to verify whether new functions improve DIR accuracy while keeping computational time (Parameter 3). The last one partially employs a new function (Parameter 4). Registration errors for these parameter settings were calculated using the manually determined landmark pairs. 3D registration errors with standard deviation over all cases were 1.78 (1.57), 1.28 (1.10), 1.44 (1.09) and 1.36 (1.35) mm for Parameter 1, 2, 3 and 4, respectively, indicating that the new functions are useful for improving DIR accuracy, even while maintaining the computational time, and this B-spline-based DIR could be used clinically to achieve high-accuracy adaptive radiotherapy.

Estimation of lung motion fields in 4D CT data by variational non-linear intensity-based registration: A comparison and evaluation study

Accurate and robust estimation of motion fields in respiration-correlated CT (4D CT) images, usually performed by non-linear registration of the temporal CT frames, is a precondition for the analysis of patient-specific breathing dynamics and subsequent image-supported diagnostics and treatment planning. In this work, we present a comprehensive comparison and evaluation study of non-linear registration variants applied to the task of lung motion estimation in thoracic 4D CT data. In contrast to existing multi-institutional comparison studies (e.g. MIDRAS and EMPIRE10), we focus on the specific but common class of variational intensity-based non-parametric registration and analyze the impact of the different main building blocks of the underlying optimization problem: the distance measure to be minimized, the regularization approach and the transformation space considered during optimization. In total, 90 different combinations of building block instances are compared. Evaluated on proprietary and publicly accessible 4D CT images, landmark-based registration errors (TRE) between 1.14 and 1.20 mm for the most accurate registration variants demonstrate competitive performance of the applied general registration framework compared to other state-of-the-art approaches for lung CT registration. Although some specific trends can be observed, effects of interchanging individual instances of the building blocks on the TRE are in general rather small (no single outstanding registration variant existing); the same level of accuracy is, however, associated with significantly different degrees of motion field smoothness and computational demands. Consequently, the building block combination of choice will depend on application-specific requirements on motion field characteristics.

Comparison of 4-dimensional computed tomography ventilation with nuclear medicine ventilation-perfusion imaging: a clinical validation study

PURPOSE

Four-dimensional computed tomography (4DCT) ventilation imaging provides lung function information for lung cancer patients undergoing radiation therapy. Before 4DCT-ventilation can be implemented clinically it needs to be validated against an established imaging modality. The purpose of this work was to compare 4DCT-ventilation to nuclear medicine ventilation, using clinically relevant global metrics and radiologist observations.

METHODS AND MATERIALS

Fifteen lung cancer patients with 16 sets of 4DCT and nuclear medicine ventilation-perfusion (VQ) images were used for the study. The VQ-ventilation images were acquired in planar mode using Tc-99m-labeled diethylenetriamine-pentaacetic acid aerosol inhalation. 4DCT data, spatial registration, and a density-change-based model were used to compute a 4DCT-based ventilation map for each patient. The percent ventilation was calculated in each lung and each lung third for both the 4DCT and VQ-ventilation scans. A nuclear medicine radiologist assessed the VQ and 4DCT scans for the presence of ventilation defects. The VQ and 4DCT-based images were compared using regional percent ventilation and radiologist clinical observations.

RESULTS

Individual patient examples demonstrate good qualitative agreement between the 4DCT and VQ-ventilation scans. The correlation coefficients were 0.68 and 0.45, using the percent ventilation in each individual lung and lung third, respectively. Using radiologist-noted presence of ventilation defects and receiver operating characteristic analysis, the sensitivity, specificity, and accuracy of the 4DCT-ventilation were 90%, 64%, and 81%, respectively.

CONCLUSIONS

The current work compared 4DCT with VQ-based ventilation using clinically relevant global metrics and radiologist observations. We found good agreement between the radiologist's assessment of the 4DCT and VQ-ventilation images as well as the percent ventilation in each lung. The agreement lessened when the data were analyzed on a regional level. Our study presents an important step for the integration of 4DCT-ventilation into thoracic clinical practice.