Dependence of ventilation image derived from 4D CT on deformable image registration and ventilation algorithms

Incorporation of Functional Lung Imaging Into Radiation Therapy Planning in Patients With Lung Cancer: A Systematic Review and Meta-Analysis

Our purpose was to provide an understanding of current functional lung imaging (FLI) techniques and their potential to improve dosimetry and outcomes for patients with lung cancer receiving radiation therapy (RT). Excerpta Medica dataBASE (EMBASE), PubMed, and Cochrane Library were searched from 1990 until April 2023. Articles were included if they reported on FLI in one of: techniques, incorporation into RT planning for lung cancer, or quantification of RT-related outcomes for patients with lung cancer. Studies involving all RT modalities, including stereotactic body RT and particle therapy, were included. Meta-analyses were conducted to investigate differences in dose-function parameters between anatomic and functional RT planning techniques, as well as to investigate correlations of dose-function parameters with grade 2+ radiation pneumonitis (RP). One hundred seventy-eight studies were included in the narrative synthesis. We report on FLI modalities, dose-response quantification, functional lung (FL) definitions, FL avoidance techniques, and correlations between FL irradiation and toxicity. Meta-analysis results show that FL avoidance planning gives statistically significant absolute reductions of 3.22% to the fraction of well-ventilated lung receiving 20 Gy or more, 3.52% to the fraction of well-perfused lung receiving 20 Gy or more, 1.3 Gy to the mean dose to the well-ventilated lung, and 2.41 Gy to the mean dose to the well-perfused lung. Increases in the threshold value for defining FL are associated with decreases in functional parameters. For intensity modulated RT and volumetric modulated arc therapy, avoidance planning results in a 13% rate of grade 2+ RP, which is reduced compared with results from conventional planning cohorts. A trend of increased predictive ability for grade 2+ RP was seen in models using FL information but was not statistically significant. FLI shows promise as a method to spare FL during thoracic RT, but interventional trials related to FL avoidance planning are sparse. Such trials are critical to understanding the effect of FL avoidance planning on toxicity reduction and patient outcomes.

Advances in CT-based lung function imaging for thoracic radiotherapy

The objective of this review is to examine the potential benefits and challenges of CT-based lung function imaging in radiotherapy over recent decades. This includes reviewing background information, defining related concepts, classifying and reviewing existing studies, and proposing directions for further investigation. The lung function imaging techniques reviewed herein encompass CT-based methods, specifically utilizing phase-resolved four-dimensional CT (4D-CT) or end-inspiratory and end-expiratory CT scans, to delineate distinct functional regions within the lungs. These methods extract crucial functional parameters, including lung volume and ventilation distribution, pivotal for assessing and characterizing the functional capacity of the lungs. CT-based lung ventilation imaging offers numerous advantages, notably in the realm of thoracic radiotherapy. By utilizing routine CT scans, additional radiation exposure and financial burdens on patients can be avoided. This imaging technique also enables the identification of different functional areas of the lung, which is crucial for minimizing radiation exposure to healthy lung tissue and predicting and detecting lung injury during treatment. In conclusion, CT-based lung function imaging holds significant promise for improving the effectiveness and safety of thoracic radiotherapy. Nevertheless, challenges persist, necessitating further research to address limitations and optimize clinical utilization. Overall, this review highlights the importance of CT-based lung function imaging as a valuable tool in radiotherapy planning and lung injury monitoring.

4DCT ventilation function image-based functional lung protection for esophageal cancer radiotherapy

BACKGROUND

4DCT (four-dimensional computed tomography) can effectively obtain functional lung ventilation images for patients and integrate them into radiotherapy treatment planning. Studies have not been performed on esophageal cancer, and there is no clear consensus on the optimal functional lung threshold for functional lung.

METHODS

Functional lung images were generated for 11 patients with esophageal cancer. The correlation between the dose-volume parameters of functional lung (FL) as defined by different thresholds and the change of PFT/PDFT (pulmonary [diffusion] function test) metrics before and after radiotherapy were evaluated. FL-sparing planning was generated for each patient to preserve the functional lung and compared to conventional anatomical CT (non-sparing) planning.

RESULTS

There was a significant positive correlation between the FL0.8 (defined Jacobian value ≤ 0.8), FL0.84, and FL0.9 dose-volume parameters and ΔFEV1/FVC (reduction before and after radiotherapy), and the FL0.8‑V₃₀ correlation was the strongest (r = 0.819, P < 0.01). The FL-sparing planning had a target area conformity index and homogeneity index comparable to the non-sparing planning (P > 0.05). For FL, the FL-sparing planning achieved lower FL-MLD (6.30 ± 2.14 Gy vs. 7.83 ± 2.70 Gy), V₁₀ (17.13 ± 7.70% vs. 27.40 ± 9.48%), and V₂₀ (6.96 ± 3.85% vs. 11.63 ± 7.19%) compared to the non-sparing planning (P < 0.05), while heart and spinal cord doses were not significantly different between the two planning groups.

CONCLUSION

The 4DCT-based FL irradiation dose for esophageal cancer was significantly associated with a decrease in FEV1/FVC. The optimal FL defined as a Jacobian value ≤ 0.8 or about 21% of the whole lung volume may be a good choice. FL-sparing planning significantly reduced the FL dose without compromising target area coverage.

Functional lung imaging in thoracic tumor radiotherapy: Application and progress

Radiotherapy plays an irreplaceable and unique role in treating thoracic tumors, but the occurrence of radiation-induced lung injury has limited the increase in tumor target doses and has influenced patients' quality of life. However, the introduction of functional lung imaging has been incorporating functional lungs into radiotherapy planning. The design of the functional lung protection plan, while meeting the target dose requirements and dose limitations of the organs at risk (OARs), minimizes the radiation dose to the functional lung, thus reducing the occurrence of radiation-induced lung injury. In this manuscript, we mainly reviewed the lung ventilation or/and perfusion functional imaging modalities, application, and progress, as well as the results based on the functional lung protection planning in thoracic tumors. In addition, we also discussed the problems that should be explored and further studied in the practical application based on functional lung radiotherapy planning.

Different functional lung-sparing strategies and radiotherapy techniques for patients with esophageal cancer

Background

Integration of 4D-CT ventilation function images into esophageal cancer radiation treatment planning aimed to assess dosimetric differences between different functional lung (FL) protection strategies and radiotherapy techniques.

Methods

A total of 15 patients with esophageal cancer who had 4D-CT scans were included. Lung ventilation function images based on Jacobian values were obtained by deformation image registration and ventilation imaging algorithm. Several different plans were designed for each patient: clinical treatment planning (non-sparing planning), the same beam distribution to FL-sparing planning, three fixed-beams FL-sparing intensity-modulated radiation therapy (IMRT) planning (5F-IMRT, 7F-IMRT, 9F-IMRT), and two FL-sparing volumetric modulated arc therapy (VMAT) planning [1F-VMAT (1-Arc), 2F-VMAT (2-Arc)]. The dosimetric parameters of the planning target volume (PTV) and organs at risk (OARs) were compared and focused on dosimetric differences in FL.

Results

The FL-sparing planning compared with the non-sparing planning significantly decreased the FL-Dmean, V5-30 and Lungs-Dmean, V10-30 (Vx: volume of receiving ≥X Gy), although it slightly compromised PTV conformability and increased Heart-V40 (P&lt; 0.05). The 5F-IMRT had the lowest PTV-conformability index (CI) but had a lower Lungs and Heart irradiation dose compared with those of the 7F-IMRT and 9F-IMRT (P&lt; 0.05). The 2F-VMAT had higher PTV-homogeneity index (HI) and reduced irradiation dose to FL, Lungs, and Heart compared to those of the 1F-VMAT planning (P&lt; 0.05). The 2F-VMAT had higher PTV conformability and homogeneity and decreased FL-Dmean, V5-20 and Lungs-Dmean, V5-10 but correspondingly increased spinal cord-Dmean compared with those of the 5F-IMRT planning (P&lt; 0.05).

Conclusion

In this study, 4D-CT ventilation function image-based FL-sparing planning for esophageal cancer can effectively reduce the dose of the FL. The 2F-VMAT planning is better than the 5F-IMRT planning in reducing the dose of FL.

Artificial intelligence in functional imaging of the lung

Artificial intelligence (AI) is transforming the way we perform advanced imaging. From high-resolution image reconstruction to predicting functional response from clinically acquired data, AI is promising to revolutionize clinical evaluation of lung performance, pushing the boundary in pulmonary functional imaging for patients suffering from respiratory conditions. In this review, we overview the current developments and expound on some of the encouraging new frontiers. We focus on the recent advances in machine learning and deep learning that enable reconstructing images, quantitating, and predicting functional responses of the lung. Finally, we shed light on the potential opportunities and challenges ahead in adopting AI for functional lung imaging in clinical settings.

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.

Technical Note: Deriving ventilation imaging from 4DCT by deep convolutional neural network

PURPOSE

Ventilation images can be derived from four-dimensional computed tomography (4DCT) by analyzing the change in HU values and deformable vector fields between different respiration phases of computed tomography (CT). As deformable image registration (DIR) is involved, accuracy of 4DCT-derived ventilation image is sensitive to the choice of DIR algorithms. To overcome the uncertainty associated with DIR, we develop a method based on deep convolutional neural network (CNN) to derive ventilation images directly from the 4DCT without explicit image registration.

METHODS

A total of 82 sets of 4DCT and ventilation images from patients with lung cancer were used in this study. In the proposed CNN architecture, the CT two-channel input data consist of CT at the end of exhale and the end of inhale phases. The first convolutional layer has 32 different kernels of size 5 × 5 × 5, followed by another eight convolutional layers each of which is equipped with an activation layer (ReLU). The loss function is the mean-squared-error (MSE) to measure the intensity difference between the predicted and reference ventilation images.

RESULTS

The predicted images were comparable to the label images of the test data. The similarity index, correlation coefficient, and Gamma index passing rate averaged over the tenfold cross validation were 0.880 ± 0.035, 0.874 ± 0.024, and 0.806 ± 0.014, respectively.

CONCLUSIONS

The results demonstrate that deep CNN can generate ventilation imaging from 4DCT without explicit deformable image registration, reducing the associated uncertainty.

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