Investigation of four-dimensional computed tomography-based pulmonary ventilation imaging in patients with emphysematous lung regions

Establishing a 4D-CT lung function related volumetric dose model to reduce radiation pneumonia

In order to study how to use pulmonary functional imaging obtained through 4D-CT fusion for radiotherapy planning, and transform traditional dose volume parameters into functional dose volume parameters, a functional dose volume parameter model that may reduce level 2 and above radiation pneumonia was obtained. 41 pulmonary tumor patients who underwent 4D-CT in our department from 2020 to 2023 were included. MIM Software (MIM 7.0.7; MIM Software Inc., Cleveland, OH, USA) was used to register adjacent phase CT images in the 4D-CT series. The three-dimensional displacement vector of CT pixels was obtained when changing from one respiratory state to another respiratory state, and this three-dimensional vector was quantitatively analyzed. Thus, a color schematic diagram reflecting the degree of changes in lung CT pixels during the breathing process, namely the distribution of ventilation function strength, is obtained. Finally, this diagram is fused with the localization CT image. Select areas with Jacobi > 1.2 as high lung function areas and outline them as fLung. Import the patient's DVH image again, fuse the lung ventilation image with the localization CT image, and obtain the volume of fLung different doses (V60, V55, V50, V45, V40, V35, V30, V25, V20, V15, V10, V5). Analyze the functional dose volume parameters related to the risk of level 2 and above radiation pneumonia using R language and create a predictive model. By using stepwise regression and optimal subset method to screen for independent variables V35, V30, V25, V20, V15, and V10, the prediction formula was obtained as follows: Risk = 0.23656-0.13784 * V35 + 0.37445 * V30-0.38317 * V25 + 0.21341 * V20-0.10209 * V15 + 0.03815 * V10. These six independent variables were analyzed using a column chart, and a calibration curve was drawn using the calibrate function. It was found that the Bias corrected line and the Apparent line were very close to the Ideal line, The consistency between the predicted value and the actual value is very good. By using the ROC function to plot the ROC curve and calculating the area under the curve: 0.8475, 95% CI 0.7237-0.9713, it can also be determined that the accuracy of the model is very high. In addition, we also used Lasso method and random forest method to filter out independent variables with different results, but the calibration curve drawn by the calibration function confirmed poor prediction performance. The function dose volume parameters V35, V30, V25, V20, V15, and V10 obtained through 4D-CT are key factors affecting radiation pneumonia. Establishing a predictive model can provide more accurate lung restriction basis for clinical radiotherapy planning.

Prospective Trial of Functional Lung Avoidance Radiation Therapy for Lung Cancer: Quality of Life Report

PURPOSE

A novel form of lung function imaging has been developed that uses 4-dimensional computed tomography (4DCT) data to generate lung ventilation images (4DCT-ventilation). Functional avoidance uses 4DCT-ventilation to reduce doses to functional lung with the aim of reducing pulmonary side effects. A phase 2, multicenter 4DCT-ventilation functional avoidance clinical trial was completed. The purpose of this work was to quantify changes in patient-reported outcomes (PROs) for patients treated with functional avoidance and determine which metrics are predictive of PRO changes.

MATERIALS AND METHODS

Patients with locally advanced lung cancer receiving curative-intent radiation therapy were accrued. Each patient had a 4DCT-ventilation image generated using 4DCT data and image processing. PRO instruments included the Functional Assessment of Cancer Therapy-Lung (FACT-L) questionnaire administered pretreatment; at the end of treatment; and at 3, 6, and 12 months posttreatment. Using the FACT-Trial Outcome Index and the FACT-Lung Cancer Subscale results, the percentage of clinically meaningful declines (CMDs) were determined. A linear mixed-effects model was used to determine which patient, clinical, dose, and dose-function metrics were predictive of PRO decline.

RESULTS

Of the 59 patients who completed baseline PRO surveys. 83% had non-small cell lung cancer, with 75% having stage 3 disease. The median dose was 60 Gy in 30 fractions. CMD FACT-Trial Outcome Index decline was 46.3%, 38.5%, and 26.8%, at 3, 6, and 12 months, respectively. CMD FACT-Lung Cancer Subscale decline was 33.3%, 33.3%, and 29.3%, at 3, 6, and 12 months, respectively. Although an increase in most dose and dose-function parameters was associated with a modest decline in PROs, none of the results were significant (all P > .053).

CONCLUSIONS

The current work presents an innovative combination of use of functional avoidance and PRO assessment and is the first report of PROs for patients treated with prospective 4DCT-ventilation functional avoidance. Approximately 30% of patients had clinically significant decline in PROs at 12 months posttreatment. The study provides additional data on outcomes with 4DCT-ventilation functional avoidance.

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.

Assessing four-dimensional CT stress maps derived from patient-specific biomechanical models of the lung with pulmonary function test data in lung cancer patients

OBJECTIVE

This study aims to retrospectively compare the stress map of the lung with pulmonary function test (PFT) results in lung cancer patients and to evaluate the potential of the stress map as an imaging biomarker for chronic obstructive pulmonary disease (COPD).

METHODS

25 lung cancer patients with pre-treatment four-dimensional CT (4DCT) and PFT data were retrospectively analysed. PFT metrics were used to diagnose obstructive lung disease. For each patient, forced expiratory volume in 1 s (FEV1 % predicted) and the ratio of FEV1 and forced vital capacity (FEV1/FVC) were recorded. 4DCT and biomechanical model-deformable image registration (BM-DIR) were used to obtain the lung stress map. The relationship between the mean of the total lung stress and PFT data was evaluated, and the COPD classification grade was also evaluated.

RESULTS

The mean values of the total lung stress and FEV1 % predicted showed a significant strong correlation [R = 0.833, (p < 0.001)]. The mean values and FEV1/FVC showed a significant strong correlation [R = 0.805, (p < 0.001)]. For the total lung stress, the area under the curve and the optimal cut-off value were 0.94 and 510.8 Pa for the classification of normal or abnormal lung function, respectively.

CONCLUSION

This study has demonstrated the potential of lung stress maps based on BM-DIR to accurately assess lung function by comparing them with PFT data.

ADVANCES IN KNOWLEDGE

The derivation of stress map directly from 4DCT is novel method. The BM-DIR-based lung stress map can provide an accurate assessment of lung function.

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.

CT-derived vessel segmentation for analysis of post-radiation therapy changes in vasculature and perfusion

Vessel segmentation in the lung is an ongoing challenge. While many methods have been able to successfully identify vessels in normal, healthy, lungs, these methods struggle in the presence of abnormalities. Following radiotherapy, these methods tend to identify regions of radiographic change due to post-radiation therapytoxicities as vasculature falsely. By combining texture analysis and existing vasculature and masking techniques, we have developed a novel vasculature segmentation workflow that improves specificity in irradiated lung while preserving the sensitivity of detection in the rest of the lung. Furthermore, radiation dose has been shown to cause vascular injury as well as reduce pulmonary function post-RT. This work shows the improvements our novel vascular segmentation method provides relative to existing methods. Additionally, we use this workflow to show a dose dependent radiation-induced change in vasculature which is correlated with previously measured perfusion changes (R² = 0.72) in both directly irradiated and indirectly damaged regions of perfusion. These results present an opportunity to extend non-contrast CT-derived models of functional change following radiation therapy.

Measuring Indirect Radiation-Induced Perfusion Change in Fed Vasculature Using Dynamic Contrast CT

Recent functional lung imaging studies have presented evidence of an “indirect effect” on perfusion damage, where regions that are unirradiated or lowly irradiated but that are supplied by highly irradiated regions observe perfusion damage post-radiation therapy (RT). The purpose of this work was to investigate this effect using a contrast-enhanced dynamic CT protocol to measure perfusion change in five novel swine subjects. A cohort of five Wisconsin Miniature Swine (WMS) were given a research course of 60 Gy in five fractions delivered locally to a vessel in the lung using an Accuray Radixact tomotherapy system with Synchrony motion tracking to increase delivery accuracy. Imaging was performed prior to delivering RT and 3 months post-RT to yield a 28−36 frame image series showing contrast flowing in and out of the vasculature. Using MIM software, contours were placed in six vessels on each animal to yield a contrast flow curve for each vessel. The contours were placed as follows: one at the point of max dose, one low-irradiated (5−20 Gy) branching from the max dose vessel, one low-irradiated (5−20 Gy) not branching from the max dose vessel, one unirradiated (<5 Gy) branching from the max dose vessel, one unirradiated (<5 Gy) not branching from the max dose vessel, and one in the contralateral lung. Seven measurements (baseline-to-baseline time and difference, slope up and down, max rise and value, and area under the curve) were acquired for each vessel’s contrast flow curve in each subject. Paired Student t-tests showed statistically significant (p < 0.05) reductions in the area under the curve in the max dose, and both fed contours indicating an overall reduction in contrast in these regions. Additionally, there were statistically significant reductions observed when comparing pre- and post-RT in slope up and down in the max dose, low-dose fed, and no-dose fed contours but not the low-dose not-fed, no-dose not-fed, or contralateral contours. These findings suggest an indirect damage effect where irradiation of the vasculature causes a reduction in perfusion in irradiated regions as well as regions fed by the irradiated vasculature.

Assessment of a computed tomography-based radiomics approach for assessing lung function in lung cancer patients

PURPOSE

We aimed to assess radiomics approaches for estimating three pulmonary function test (PFT) results (forced expiratory volume in one second [FEV1], forced vital capacity [FVC], and the ratio of FEV1 to FVC [FEV₁/FVC]) using data extracted from chest computed tomography (CT) images.

METHODS

This retrospective study included 85 lung cancer patients (mean age, 75 years ±8; 69 men) who underwent stereotactic body radiotherapy between 2012 and 2020. Their pretreatment chest breath-hold CT and PFT data before radiotherapy were obtained. A total of 107 radiomics features (Shape: 14, Intensity: 18, Texture: 75) were extracted using two methods: extraction of the lung tissue (<-250 HU) (APPROACH 1), and extraction of small blood vessels and lung tissue (APPROACH 2). The PFT results were estimated using the least absolute shrinkage and selection operator regression. Pearson's correlation coefficients (r) were determined for all PFT results, and the area under the curve (AUC) was calculated for FEV₁/FVC (<70 %). Finally, we compared our approaches with the conventional formula (Conventional).

RESULTS

For the estimated FEV₁/FVC, the Pearson's r were 0.21 (P =.06), 0.69 (P <.01), and 0.73 (P <.01) for Conventional, APPROACH 1, and APPROACH 2, respectively; the AUCs for FEV₁/FVC (<70 %) were 0.67 (95 % confidence interval [CI]: 0.55, 0.79), 0.82 (CI: 0.72, 0.91; P =.047) and 0.86 (CI: 0.78, 0.94; P =.01), respectively.

CONCLUSIONS

The radiomics approach performed better than the conventional equation and may be useful for assessing lung function based on CT images.

Combining Serial and Parallel Functionality in Functional Lung Avoidance Radiation Therapy

PURPOSE

Functional lung avoidance (FLA) radiation therapy (RT) aims to minimize post-RT pulmonary toxicity by preferentially avoiding dose to high-functioning lung (HFL) regions. A common limitation is that FLA approaches do not consider the conducting architecture for gas exchange. We previously proposed the functionally weighted airway sparing (FWAS) method to spare airways connected to HFL regions, showing that it is possible to substantially reduce risk of radiation-induced airway injury. Here, we compare the performance of FLA and FWAS and propose a novel method combining both approaches.

METHODS

We used breath-hold computed tomography (BHCT) and simulation 4-dimensional computed tomography (4DCT) from 12 lung stereotactic ablative radiation therapy patients. Four planning strategies were examined: (1) Conventional: no sparing other than clinical dose-volume constraints; (2) FLA: using a 4DCT-based ventilation map to delineate the HFL, plans were optimized to reduce mean dose and V_13.50 in HFL; (3) FWAS: we autosegemented 11 to 13 generations of individual airways from each patient's BHCT and assigned priorities based on the relative contribution of each airway to total ventilation. We used these priorities in the optimization along with airway dose constraints, estimated as a function of airway diameter and 5% probability of collapse; and (4) FLA + FWAS: we combined information from the 2 strategies. We prioritized clinical dose constraints for organs at risk and planning target volume in all plans. We performed the evaluation in terms of ventilation preservation accounting for radiation-induced damage to both lung parenchyma and airways.

RESULTS

We observed average ventilation preservation for FLA, FWAS, and FLA + FWAS as 3%, 8.5%, and 14.5% higher, respectively, than for Conventional plans for patients with ventilation preservation in Conventional plans <90%. Generalized estimated equations showed that all improvements were statistically significant (P ≤ .036). We observed no clinically relevant improvements in outcomes of the sparing techniques in patients with ventilation preservation in Conventional plans ≥90%.

CONCLUSIONS

These initial results suggest that it is crucial to consider the parallel and the serial nature of the lung to improve post-radiation therapy lung function and, consequently, quality of life for patients.

Deriving Pulmonary Ventilation Images From Clinical 4D-CBCT Using a Deep Learning-Based Model

Purpose

The current algorithms for measuring ventilation images from 4D cone-beam computed tomography (CBCT) are affected by the accuracy of deformable image registration (DIR). This study proposes a new deep learning (DL) method that does not rely on DIR to derive ventilation images from 4D-CBCT (CBCT-VI), which was validated with the gold-standard single-photon emission-computed tomography ventilation image (SPECT-VI).

Materials and Methods

This study consists of 4D-CBCT and 99mTc-Technegas SPECT/CT scans of 28 esophagus or lung cancer patients. The scans were rigidly registered for each patient. Using these data, CBCT-VI was derived using a deep learning-based model. Two types of model input data are studied, namely, (a) 10 phases of 4D-CBCT and (b) two phases of peak-exhalation and peak-inhalation of 4D-CBCT. A sevenfold cross-validation was applied to train and evaluate the model. The DIR-dependent methods (density-change-based and Jacobian-based methods) were used to measure the CBCT-VIs for comparison. The correlation was calculated between each CBCT-VI and SPECT-VI using voxel-wise Spearman's correlation. The ventilation images were divided into high, medium, and low functional lung regions. The similarity of different functional lung regions between SPECT-VI and each CBCT-VI was evaluated using the dice similarity coefficient (DSC). One-factor ANONA model was used for statistical analysis of the averaged DSC for the different methods of generating ventilation images.

Results

The correlation values were 0.02 ± 0.10, 0.02 ± 0.09, and 0.65 ± 0.13/0.65 ± 0.15, and the averaged DSC values were 0.34 ± 0.04, 0.34 ± 0.03, and 0.59 ± 0.08/0.58 ± 0.09 for the density change, Jacobian, and deep learning methods, respectively. The strongest correlation and the highest similarity with SPECT-VI were observed for the deep learning method compared to the density change and Jacobian methods.

Conclusion

The results showed that the deep learning method improved the accuracy of correlation and similarity significantly, and the derived CBCT-VIs have the potential to monitor the lung function dynamic changes during radiotherapy.

Predicting radiation pneumonitis with fuzzy clustering neural network using 4DCT ventilation image based dosimetric parameters

Background

To develop a fuzzy clustering neural network to predict radiation-induced pneumonitis (RP) using four-dimensional computed tomography (4DCT) ventilation image (VI) based dosimetric parameters for thoracic cancer patients.

Methods

The VI were retrospectively calculated from pre-treatment 4DCT data using a deformable image registration (DIR) and an improved VI algorithm. Similar to dose-volume histogram (DVH) of intensity modulated radiotherapy (IMRT), dose-function histogram (DFH) was derived from dose distribution and VI. Then, the dose-function metrics were calculated from DFH. For comparison, the dose-volume metrics were calculated from DVH. Correspondingly, two sets of feature vectors were formed from the dose-volume metrics and the dose-function metrics, respectively. For the feature vectors of each set, they were first pre-processed by principal component analysis (PCA) to reduce feature dimensions. Then, they were grouped to few clusters determined by fuzzy c-means (FCM) algorithm. Next, the neural network was trained to correlate the dosimetric parameters with RP based on the feature vectors of each cluster. Finally, the occurrence of RP was predicted by the neural network on the test data.

Results

Through PCA analysis, the top 5 principal components were selected. Their contribution is more than 98%, which is adequate to represent the original feature space of input data. Based on the clustering validity indexes, the optimal number of clusters is 4 and used for subsequent fuzzy clustering of the input data. After network training, the areas under the curve (AUC) of the prediction model is 0.77 using VI-based dosimetric parameters and 0.67 using structure-based dosimetric parameters.

Conclusions

Compared to the structure-based dosimetric features, the VI-based dosimetric features are more relevant to lung function and presented higher prediction accuracy of RP. The fuzzy clustering neural network improved the prediction accuracy of RP compared to the conventional neural network. The combination of VI-based dose-function metrics and fuzzy clustering neural network provides an effective predictive model for assessing lung toxicity risk after radiotherapy.

Radiation-induced airway changes and downstream ventilation decline in a swine model

Purpose.To investigate indirect radiation-induced changes in airways as precursors to atelectasis post radiation therapy (RT).Methods.Three Wisconsin Miniature Swine (WMSTM) underwent a research course of 60 Gy in 5 fractions delivered to a targeted airway/vessel in the inferior left lung. The right lung received a max point dose <5 Gy. Airway segmentation was performed on the pre- and three months post-RT maximum inhale phase of the four-dimensional (4D) computed tomography (CT) scans. Changes in luminal area (Ai) and square root of wall area (WA) for each airway were investigated. Changes in ventilation were assessed using the Jacobian ratio and were measured in three different regions: the inferior left lung <5 Gy (ILL), the superior left lung <5 Gy (SLL), and the contralateral right lung <5 Gy (RL).Results.Airways (n = 25) in the right lung for all swine showed no significant changes (p = 0.48) in Ai post-RT compared to pre-RT. Airways (n = 28) in the left lung of all swine were found to have a significant decrease (p < 0.001) in Ai post-RT compared to pre-RT, correlated (Pearson R = -0.97) with airway dose. Additionally,WAdecreased significantly (p < 0.001) with airway dose. Lastly, the Jacobian ratio of the ILL (0.883) was lower than that of the SLL (0.932) and the RL (0.955).Conclusions.This work shows that for the swine analyzed, there were significant correlations between Ai andWAchange with radiation dose. Additionally, there was a decrease in lung function in the regions of the lung supplied by the irradiated airways compared to the regions supplied by unirradiated airways. These results support the hypothesis that airway dose should be considered during treatment planning in order to potentially preserve functional lung and reduce lung toxicities.

Dynamic chest radiography: clinical validation of ventilation and perfusion metrics derived from changes in radiographic lung density compared to nuclear medicine imaging

Background

Dynamic chest radiography (DCR) is a type of non-contrast-enhanced functional lung imaging with a dynamic flat-panel detector (FPD). This study aimed to assess the clinical significance of ventilation and perfusion metrics derived from changes in radiographic lung density on DCR in comparison to nuclear medicine imaging-derived metrics.

Methods

DCR images of 42 lung cancer patients were sequentially obtained during respiration using a dynamic FPD imaging system. For each subdivided lung region, the maximum change in the averaged pixel value (Δ_max), i.e., lung density, due to respiration and cardiac function was calculated, and the percentage of Δ_max relative to the total of all lung regions (Δ_max%) was computed for ventilation and perfusion, respectively. The Δ_max% was compared to the accumulation of radioactive agents such as Tc-99m gas and Tc-99m macro-aggregated albumin (radioactive agents%) on ventilation and perfusion scans in the subdivided lung regions, by Spearman's correlation coefficient (r) and the Dice similarity coefficients (DSC). To facilitate visual evaluation, Δ_max% was visualized as a color scaling, where larger Δ_max values were indicated by higher color intensities.

Results

We found a moderate correlation between Δ_max% and radioactive agents% on ventilation and perfusion scans, with perfusion metrics (r=0.57, P<0.001) showing a higher correlation than ventilation metrics (r=0.53, P<0.001). We also found a good or strong correlation (r≥0.5) in 80.9% (34/42) of patients for perfusion metrics (r=0.60±0.16) and in 52.4% (22/42) of patients for ventilation metrics (r=0.53±0.16). DSC indicated a moderate correlation for both metrics. Decreased pulmonary function was observed in the form of reduced color intensities on color-mapping images.

Conclusions

DCR-derived ventilation and perfusion metrics correlated reasonably well with nuclear medicine imaging findings in lung subdivisions, suggesting that DCR could provide useful information on pulmonary function without the use of radioactive contrast agents.

In Silico Ventilation Within the Dose-Volume is Predictive of Lung Function Post-radiation Therapy in Patients with Lung Cancer

Lung cancer is a leading cause of death worldwide. Radiation therapy (RT) is one method to treat this disease. A common side effect of RT for lung cancer is radiation-induced lung damage (RILD) which leads to loss of lung function. RILD often compounds pre-existing smoking-related regional lung function impairment. It is difficult to predict patient outcomes due to large variability in individual response to RT. In this study, the capability of image-based modelling of regional ventilation in lung cancer patients to predict lung function post-RT was investigated. Twenty-five patient-based models were created using CT images to define the airway geometry, size and location of tumour, and distribution of emphysema. Simulated ventilation within the 20 Gy isodose volume showed a statistically significant negative correlation with the change in forced expiratory volume in 1 s 12-months post-RT (p = 0.001, R = - 0.61). Patients with higher simulated ventilation within the 20 Gy isodose volume had a greater loss in lung function post-RT and vice versa. This relationship was only evident with the combined impact of tumour and emphysema, with the location of the emphysema relative to the dose-volume being important. Our results suggest that model-based ventilation measures can be used in the prediction of patient lung function post-RT.

Functionally weighted airway sparing (FWAS): a functional avoidance method for preserving post-treatment ventilation in lung radiotherapy

Recent changes to the guidelines for screening and early diagnosis of lung cancer have increased the interest in preserving post-radiotherapy lung function. Current investigational approaches are based on spatially mapping functional regions and generating regional avoidance plans that preferentially spare highly ventilated/perfused lung. A potentially critical, yet overlooked, aspect of functional avoidance is radiation injury to peripheral airways, which serve as gas conduits to and from functional lung regions. Dose redistribution based solely on regional function may cause irreparable damage to the 'supply chain'. To address this deficiency, we propose the functionally weighted airway sparing (FWAS) method. FWAS (i) maps the bronchial pathways to each functional sub-lobar lung volume; (ii) assigns a weighting factor to each airway based on the relative contribution of the sub-volume to overall lung function; and (iii) creates a treatment plan that aims to preserve these functional pathways. To evaluate it, we used four cases from a retrospective cohort of SAbR patients treated for lung cancer. Each patient's airways were auto-segmented from a diagnostic-quality breath-hold CT using a research virtual bronchoscopy software. A ventilation map was generated from the planning 4DCT to map regional lung function. For each terminal airway, as resolved by the segmentation software, the total ventilation within the sub-lobar volume supported by that airway was estimated and used as a function-based weighting factor. Upstream airways were weighted based on the cumulative volumetric ventilation supported by corresponding downstream airways. Using a previously developed model for airway radiosensitivity, dose constraints were determined for each airway corresponding to a <5% probability of airway collapse. Airway dose constraints, ventilation scores, and clinical dose constraints were input to a swarm optimization-based inverse planning engine to create a 3D conformal SAbR plan (CRT). The FWAS plans were compared to the patients' prescribed CRT clinical plans and the inverse-optimized clinical plans. Depending on the size and location of the tumour, the FWAS plan showed superior preservation of ventilation due to airflow preservation through open pathways (i.e. cumulative ventilation score from the sub-lobar volumes of open pathways). Improvements ranged between 3% and 23%, when comparing to the prescribed clinical plans, and between 3% and 35%, when comparing to the inverse-optimized clinical plans. The three plans satisfied clinical requirements for PTV coverage and OAR dose constraints. These initial results suggest that by sparing pathways to high-functioning lung subregions it is possible to reduce post-SAbR loss of respiratory function.

Comparison of dynamic flat-panel detector-based chest radiography with nuclear medicine ventilation-perfusion imaging for the evaluation of pulmonary function: A clinical validation study

PURPOSE

Dynamic chest radiography (DCR) is a flat-panel detector (FPD)-based functional lung imaging technique capable of measuring temporal changes in radiographic lung density due to ventilation and perfusion. The aim of this study was to determine the diagnostic performance of DCR in the evaluation of pulmonary function based on changes in radiographic lung density compared to nuclear medicine lung scans.

METHODS

This study included 53 patients with pulmonary disease who underwent DCR and nuclear medicine imaging at our institution. Dynamic chest radiography was conducted using a dynamic FPD system to obtain sequential chest radiographs during one breathing cycle. The maximum change in the average pixel value (Δ_max ) was measured, and the percentage ofΔ_max in each lung region, calculated relative to the sum of all lung regions (Δ_max %), was calculated for each factor (ventilation and perfusion). The Δ_max % was compared with the accumulation of radioactive agents (radioactive agents%) on ventilation and perfusion scans in each lung and lung region using correlation coefficients and scatter plots. The ratio of ventilation to perfusion Δ_max % was calculated and compared with nuclear medicine ventilation-perfusion (V/Q) findings in terms of sensitivity and specificity for V/Q mismatch in each lung region.

RESULTS

There was a high correlation between Δ_max % and radioactive agents% for each lung (Ventilation: r = 0.81, perfusion: r = 0.87). However, correlation coefficients were lower (0.37 to 0.80) when comparing individual lung regions, with the upper lung regions showing the lowest correlation coefficients. The sensitivity and specificity of DCR for V/Q mismatch were 63.3% (19/30) and 60.1% (173/288), respectively. Motion artifacts occasionally increased Δ_max %, resulting in false negatives.

CONCLUSIONS

Ventilation and perfusion Δ_max % correlated reasonably with radioactive agents% on ventilation and perfusion scans. Although the regional correlations were lower and the detection performance for V/Q mismatch was not enough for clinical use at the moment, these results suggest the potential for DCR to be used as a functional imaging modality that can be performed without the use of radioactive contrast agents. Further technical improvement is required for the implementation of DCR-based V/Q studies.

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.

Accounting for respiratory motion in small serial structures during radiotherapy planning: proof of concept in virtual bronchoscopy-guided lung functional avoidance radiotherapy

Respiratory motion management techniques in radiotherapy (RT) planning are primarily focused on maintaining tumor target coverage. An inadequately addressed need is accounting for motion in dosimetric estimations in smaller serial structures. Accurate dose estimations in such structures are more sensitive to motion because respiration can cause them to move completely in or out of a high dose-gradient field. In this work, we study three motion management strategies (m1-m3) to find an accurate method to estimate the dosimetry in airways. To validate these methods, we generated a 'ground truth' digital breathing model based on a 4DCT scan from a lung stereotactic ablative radiotherapy (SAbR) patient. We simulated 225 breathing cycles with  ±10% perturbations in amplitude, respiratory period, and time per respiratory phase. A high-resolution breath-hold CT (BHCT) was also acquired and used with a research virtual bronchoscopy software to autosegment 239 airways. Contours for planning target volume (PTV) and organs at risk (OARs) were defined on the maximum intensity projection of the 4DCT (CT_MIP) and transferred to the average of the 10 4DCT phases (CT_AVG). To design the motion management methods, the RT plan was recreated using different images and structure definitions. Methods m1 and m2 recreated the plan using the CT_AVG image. In method m1, airways were deformed to the CT_AVG. In m2, airways were deformed to each of the 4DCT phases, and union structures were transferred onto the CT_AVG. In m3, the RT plan was recreated on each of the 10 phases, and the dose distribution from each phase was deformed to the BHCT and summed. Dose errors (mean [min, max]) in airways were: m1: 21% (0.001%, 93%); m2: 45% (0.1%, 179%); and m3: 4% (0.006%, 14%). Our work suggests that accurate dose estimation in moving small serial structures requires customized motion management techniques (like m3 in this work) rather than current clinical and investigational approaches.

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.

Availability of a simplified lung ventilation imaging algorithm based on four-dimensional computed tomography

PURPOSE

It is still not conclusive which four-dimensional computed tomography (4DCT)-based ventilation imaging algorithm is most accurate and efficient. In this study, we proposed a simplified algorithm (VIAAVG) which only requires the average computed tomography (AVG CT) as input, and quantitatively compared its accuracy and efficiency with three other popular algorithms.

MATERIAL AND METHODS

Fifty patients with lung or esophageal cancer who underwent radiotherapy were enrolled. Single photon emission computed tomography (SPECT) ventilation images (VI-SPECT) and 4DCT were acquired 1-3 days before the first treatment session. The end of exhalation and the end of inhalation CT were registered to derive deformable vector field (DVF) using MIMvista. 4DCT-based ventilation images (CTVI) were first calculated respectively by means of four algorithms (VIAJAC, VIAHU, VIAPRO and VIAAVG). The computation times were compared using paired t-test. The corresponding CTVIs (CTVIJAC, CTVIHU, CTVIPRO and CTVIAVG) and VI-SPECT were segmented into three equal sub-volumes (high, medium and low function lung, respectively) after smoothing and normalization. The Dice Similarity Coefficients (DSCs) were calculated for each sub-volume between each CTVI and VI-SPECT. The average DSCs for high, medium and low function lung in different CTVIs for each patient were compared using paired t-test.

RESULTS

The mean DSCs for CTVIJAC, CTVIHU, CTVIPRO and CTVIAVG were 0.3255, 0.4465, 0.5865 and 0.5958, respectively. The average computation times for CTVIJAC, CTVIHU, CTVIPRO and CTVIAVG were 18.3 s, 24.2 s, 144.8 s and 15.0 s.

CONCLUSION

VIAAVG is available for clinical use because of its high accuracy, improved efficiency and less input requirement compared to the other 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.

SPECT V/Q in Lung Cancer Radiotherapy Planning

Curative-intent lung cancer radiation therapy either alone (RT) or combined with immuno-chemotherapy is associated with potential risk of serious radiation-induced lung injury. This review provides a summary of the role of SPECT ventilation perfusion (V/Q) imaging as an emerging adjunct to lung cancer RT planning and treatment dosimetry. Denoted "functional lung avoidance RT" it is hypothesized that preferential dosimetric avoidance of physiologically functional lung may reduce the frequency of radiation-induced lung injury. SPECT V/Q imaging datasets available during the planning process allows the prioritization (or "personalization') of RT dose to minimize the volume of functional lung probabilistically exposed to injurious radiation dose. Selective escalation of target dose and adaptive planning and replanning is also enabled. The emergent importance of the tumor-lung microenvironment and its biologic relationship to local immune effectors in lung cancer provides further incentive to individualize RT planning and delivery. This review examines important normal tissue dosimetric constraints that are part of current standards-of-care and the new dosimetric parameters associated with functional lung avoidance RT. SPECT V/Q has been a valuable tool in investigating the feasibility and efficacy of functional lung avoidance RT but is yet to become main stream due to the lack of large clinical trials. It is encouraging however that functional lung avoidance is feasible in RT dose-target delineation and some of the more promising studies are discussed.

Virtual Bronchoscopy-Guided Treatment Planning to Map and Mitigate Radiation-Induced Airway Injury in Lung SAbR

PURPOSE

Radiation injury to the bronchial tree is an important yet poorly understood potential side effect in lung stereotactic ablative radiation therapy (SAbR). We investigate the integration of virtual bronchoscopy in radiation therapy planning to quantify dosage to individual airways. We develop a risk model of airway collapse and develop treatment plans that reduce the risk of radiation-induced airway injury.

METHODS AND MATERIALS

Pre- and post-SAbR diagnostic-quality computerized tomography (CT) scans were retrospectively collected from 26 lung cancer patients. From each scan, the bronchial tree was segmented using a virtual bronchoscopy system and registered deformably to the planning CT. Univariate and stepwise multivariate Cox regressions were performed, examining factors such as age, comorbidities, smoking pack years, airway diameter, and maximum point dosage (D_max). Logistic regression was utilized to formulate a risk function of segmental collapse based on D_max and diameter. The risk function was incorporated into the objective function along with clinical dosage volume constraints for planning target volume (PTV) and organs at risk (OARs).

RESULTS

Univariate analysis showed that segmental diameter (P = .014) and D_max (P = .007) were significantly correlated with airway segment collapse. Multivariate stepwise Cox regression showed that diameter (P = .015), D_max (P < .0001), and pack/years of smoking (P = .02) were significant independent factors associated with collapse. Risk management-based plans enabled significant dosage reduction to individual airway segments while fulfilling clinical dosimetric objectives.

CONCLUSION

To our knowledge, this is the first systematic investigation of functional avoidance in lung SAbR based on mapping and minimizing doses to individual bronchial segments. Our early results show that it is possible to substantially lower airway dosage. Such dosage reduction may potentially reduce the risk of radiation-induced airway injury, while satisfying clinically prescribed dosimetric objectives.

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.

Functional-guided radiotherapy using knowledge-based planning

BACKGROUND AND PURPOSE

There are two significant challenges when implementing functional-guided radiotherapy using 4DCT-ventilation imaging: (1) lack of knowledge of realistic patient specific dosimetric goals for functional lung and (2) ensuring consistent plan quality across multiple planners. Knowledge-based planning (KBP) is positioned to address both concerns.

MATERIAL AND METHODS

A KBP model was created from 30 previously planned functional-guided lung patients. Standard organs at risk (OAR) in lung radiotherapy and a ventilation contour delineating areas of high ventilation were included. Model validation compared dose-metrics to standard OARs and functional dose-metrics from 20 independent cases that were planned with and without KBP.

RESULTS

A significant improvement was observed for KBP optimized plans in V20Gy and mean dose to functional lung (p = 0.005 and 0.001, respectively), V20Gy and mean dose to total lung minus GTV (p = 0.002 and 0.01, respectively), and mean doses to esophagus (p = 0.005).

CONCLUSION

The current work developed a KBP model for functional-guided radiotherapy. Modest, but statistically significant, improvements were observed in functional lung and total lung doses.

A prospective gating method to acquire a diverse set of free-breathing CT images for model-based 4DCT

Breathing motion modeling requires observation of tissues at sufficiently distinct respiratory states for proper 4D characterization. This work proposes a method to improve sampling of the breathing cycle with limited imaging dose. We designed and tested a prospective free-breathing acquisition protocol with a simulation using datasets from five patients imaged with a model-based 4DCT technique. Each dataset contained 25 free-breathing fast helical CT scans with simultaneous breathing surrogate measurements. Tissue displacements were measured using deformable image registration. A correspondence model related tissue displacement to the surrogate. Model residual was computed by comparing predicted displacements to image registration results. To determine a stopping criteria for the prospective protocol, i.e. when the breathing cycle had been sufficiently sampled, subsets of N scans where 5  ⩽  N  ⩽  9 were used to fit reduced models for each patient. A previously published metric was employed to describe the phase coverage, or 'spread', of the respiratory trajectories of each subset. Minimum phase coverage necessary to achieve mean model residual within 0.5 mm of the full 25-scan model was determined and used as the stopping criteria. Using the patient breathing traces, a prospective acquisition protocol was simulated. In all patients, phase coverage greater than the threshold necessary for model accuracy within 0.5 mm of the 25 scan model was achieved in six or fewer scans. The prospectively selected respiratory trajectories ranked in the (97.5  ±  4.2)th percentile among subsets of the originally sampled scans on average. Simulation results suggest that the proposed prospective method provides an effective means to sample the breathing cycle with limited free-breathing scans. One application of the method is to reduce the imaging dose of a previously published model-based 4DCT protocol to 25% of its original value while achieving mean model residual within 0.5 mm.

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.

Evaluating the Toxicity Reduction With Computed Tomographic Ventilation Functional Avoidance Radiation Therapy

PURPOSE

Computed tomographic (CT) ventilation imaging is a new modality that uses 4-dimensional (4D) CT information to calculate lung ventilation. Although retrospective studies have reported on the reduction in dose to functional lung, no work to our knowledge has been published in which the dosimetric improvements have been translated to a reduction in the probability of pulmonary toxicity. Our work estimates the reduction in toxicity for CT ventilation-based functional avoidance planning.

METHODS AND MATERIALS

Seventy previously treated lung cancer patients who underwent 4DCT imaging were used for the study. CT ventilation maps were calculated with 4DCT deformable image registration and a density change-based algorithm. Pneumonitis was graded on the basis of imaging and clinical presentation. Maximum likelihood methods were used to generate normal tissue complication probability (NTCP) models predicting grade 2 or higher (2+) and grade 3+ pneumonitis as a function of dose (V5 Gy, V10 Gy, V20 Gy, V30 Gy, and mean dose) to functional lung. For 30 patients a functional plan was generated with the goal of reducing dose to the functional lung while meeting Radiation Therapy Oncology Group 0617 constraints. The NTCP models were applied to the functional plans and the clinically used plans to calculate toxicity reduction.

RESULTS

By the use of functional avoidance planning, absolute reductions in grade 2+ NTCP of 6.3%, 7.8%, and 4.8% were achieved based on the mean fV20 Gy, fV30 Gy, and mean dose to functional lung metrics, respectively. Absolute grade 3+ NTCP reductions of 3.6%, 4.8%, and 2.4% were achieved with fV20 Gy, fV30 Gy, and mean dose to functional lung. Maximum absolute reductions of 52.3% and 16.4% were seen for grade 2+ and grade 3+ pneumonitis for individual patients.

CONCLUSION

Our study quantifies the possible toxicity reduction from CT ventilation-based functional avoidance planning. Reductions in grades 2+ and 3+ pneumonitis were 7.1% and 4.7% based on mean dose-function metrics, with reductions as high as 52.3% for individual patients. Our work provides seminal data for determining the potential toxicity benefit from incorporating CT ventilation into thoracic treatment planning.

Evaluating Which Dose-Function Metrics Are Most Critical for Functional-Guided Radiation Therapy

PURPOSE

Four-dimensional (4D) computed tomography (CT) ventilation imaging is increasingly being used to calculate lung ventilation and implement functional-guided radiation therapy in clinical trials. There has been little exhaustive work evaluating which dose-function metrics should be used for treatment planning and plan evaluation. The purpose of our study was to evaluate which dose-function metrics best predict for radiation pneumonitis (RP).

METHODS AND MATERIALS

Seventy lung cancer patients who underwent 4D CT imaging and pneumonitis grading were assessed. Pretreatment 4D CT scans of each patient were used to calculate ventilation images. We evaluated 3 types of dose-function metrics that combined the patient's 4D CT ventilation image and treatment planning dose distribution: (1) structure-based approaches; (2) image-based approaches using the dose-function histogram; and (3) nonlinear weighting schemes. Log-likelihood methods were used to generate normal tissue complication probability models predicting grade 3 or higher (ie, grade 3+) pneumonitis for all dose-function schemes. The area under the curve (AUC) was used to assess the predictive power of the models. All techniques were compared with normal tissue complication probability models based on traditional, total lung dose metrics.

RESULTS

The most predictive models were structure-based approaches that focused on the volume of functional lung receiving ≥20 Gy (AUC, 0.70). Probabilities of grade 3+ RP of 20% and 10% correspond to V20 (percentage of volume receiving ≥20 Gy) to the functional subvolumes of 26.8% and 9.3%, respectively. Imaging-based analysis with the dose-function histogram and nonlinear weighted ventilation values yielded AUCs of 0.66 and 0.67, respectively, when we evaluated the percentage of functionality receiving ≥20 Gy. All dose-function metrics outperformed the traditional dose metrics (mean lung dose, AUC of 0.55).

CONCLUSIONS

A full range of dose-function metrics and functional thresholds was examined. The calculated AUC values for the most predictive functional models occupied a narrow range (0.66-0.70), and all showed notable improvements over AUC from traditional lung dose metrics (0.55). Identifying the combinations most predictive of grade 3+ RP provides valuable data to inform the functional-guided radiation therapy process.

A four-dimensional computed tomography comparison of healthy and asthmatic human lungs

The purpose of this study was to explore new insights in non-linearity, hysteresis and ventilation heterogeneity of asthmatic human lungs using four-dimensional computed tomography (4D-CT) image data acquired during tidal breathing. Volumetric image data were acquired for 5 non-severe and one severe asthmatic volunteers. Besides 4D-CT image data, function residual capacity and total lung capacity image data during breath-hold were acquired for comparison with dynamic scans. Quantitative results were compared with the previously reported analysis of five healthy human lungs. Using an image registration technique, local variables such as regional ventilation and anisotropic deformation index (ADI) were estimated. Regional ventilation characteristics of non-severe asthmatic subjects were similar to those of healthy subjects, but different from the severe asthmatic subject. Lobar airflow fractions were also well correlated between static and dynamic scans (R²>0.84). However, local ventilation heterogeneity significantly increased during tidal breathing in both healthy and asthmatic subjects relative to that of breath-hold perhaps because of airway resistance present only in dynamic breathing. ADI was used to quantify non-linearity and hysteresis of lung motion during tidal breathing. Non-linearity was greater on inhalation than exhalation among all subjects. However, exhalation non-linearity among asthmatic subjects was greater than healthy subjects and the difference diminished during inhalation. An increase of non-linearity during exhalation in asthmatic subjects accounted for lower hysteresis relative to that of healthy ones. Thus, assessment of non-linearity differences between healthy and asthmatic lungs during exhalation may provide quantitative metrics for subject identification and outcome assessment of new interventions.

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.

Improving the Accuracy of Registration-Based Biomechanical Analysis: A Finite Element Approach to Lung Regional Strain Quantification

Tissue deformation plays an important role in lung physiology, as lung parenchyma largely deforms during spontaneous ventilation. However, excessive regional deformation may lead to lung injury, as observed in patients undergoing mechanical ventilation. Thus, the accurate estimation of regional strain has recently received great attention in the intensive care community. In this work, we assess the accuracy of regional strain maps computed from direct differentiation of B-Spline (BS) interpolations, a popular technique employed in non-rigid registration of lung computed tomography (CT) images. We show that, while BS-based registration methods give excellent results for the deformation transformation, the strain field directly computed from BS derivatives results in predictions that largely oscillate, thus introducing important errors that can even revert the sign of strain. To alleviate such spurious behavior, we present a novel finite-element (FE) method for the regional strain analysis of lung CT images. The method follows from a variational strain recovery formulation, and delivers a continuous approximation to the strain field in arbitrary domains. From analytical benchmarks, we show that the FE method results in errors that are a fraction of those found for the BS method, both in an average and pointwise sense. The application of the proposed FE method to human lung CT images results in 3D strain maps are heterogeneous and smooth, showing high consistency with specific ventilation maps reported in the literature. We envision that the proposed FE method will considerably improve the accuracy of image-based biomechanical analysis, making it reliable enough for routine medical applications.

Imaging regional lung function: a critical tool for developing inhaled antimicrobial therapies

Alterations in regional lung function due to respiratory infection have a significant effect on the deposition of inhaled treatments. This has consequences for treatment effectiveness and hence recovery of lung function. In order to advance our understanding of respiratory infection and inhaled treatment delivery, we must develop imaging techniques that can provide regional functional measurements of the lung. In this review, we explore the role of functional imaging for the assessment of respiratory infection and development of inhaled treatments. We describe established and emerging functional lung imaging methods. The effect of infection on lung function is described, and the link between regional disease, function, and inhaled treatments is discussed. The potential for lung function imaging to provide unique insights into the functional consequences of infection, and its treatment, is also discussed.

Quantitative assessment of pulmonary perfusion using dynamic contrast-enhanced CT in patients with chronic obstructive pulmonary disease: correlations with pulmonary function test and CT volumetric parameters

BACKGROUND

Pulmonary function test (PFT) is commonly used to help diagnose chronic obstructive pulmonary disease (COPD) and other lung diseases. However, it cannot be used to evaluate regional function and morphological abnormalities.

PURPOSE

To quantitatively evaluate pulmonary perfusion imaging using dynamic contrast-enhanced (DCE) computed tomography (CT) and observe its correlations with PFT and CT volumetric parameters in COPD patients.

MATERIAL AND METHODS

PFT and CT pulmonary perfusion examination were performed in 63 COPD patients. Perfusion defects were quantitated by calculating the CT value ratio (RHU) between perfusion defects (HUdefect) and normal lung (HUnormal). Volumetric CT data were used to calculate emphysema index (EI), total lung volume (TLV), and total emphysema volume (TEV). Emphysematous parenchyma was defined as the threshold of lung area lower than -950 HU. Correlations between RHU and TLV, TEV, EI, and PFT were assessed using Spearman correlation analysis.

RESULTS

The positive rate of perfusion defects on CT perfusion images was higher than that of emphysema on CT mask images (χ(2) = 17.027, P < 0.001). The Spearman correlation test showed that RHU was positively correlated with FEV1 (R = 0.59, P < 0.001), FEV1% Predicted (R = 0.61, P < 0.001), FVC (R = 0.47, P = 0.002), and FEV1/FVC (R = 0.65, P < 0.001), and negatively correlated with EI (R = -0.67, P < 0.001).

CONCLUSION

CT perfusion imaging is more sensitive in detecting emphysema that is inconspicuous on CT images. RHU is correlated with PFT and CT volumetric parameters, suggesting that it is more sensitive in detecting early COPD changes and may prove to be a potential predictor of focal lung function.

Dosimetric impact of 4-dimensional computed tomography ventilation imaging-based functional treatment planning for stereotactic body radiation therapy with 3-dimensional conformal radiation therapy

PURPOSE

The purpose of this study was to clarify the dosimetric impact of 4-dimensional computed tomography (4D-CT)-derived, ventilation-guided functional avoidance for stereotactic body radiation therapy (SBRT) with 3-dimensional conformal radiation therapy.

METHODS AND MATERIALS

Eleven lung cancer patients with peripheral tumors no greater than 5 cm in size were studied (average planning target volume, 42.4 ± 32.5 cm(3)). Four-dimensional-CT ventilation imaging was performed using deformable image registration for spatial mapping of the peak-exhale 4D-CT image to the peak-inhale 4D-CT image and computation of the Jacobian-based ventilation metric. For each patient, anatomical and functional plans were created using 7 to 9 noncoplanar beams for SBRT (40-56 Gy/4-8 fractions). The anatomical plan was generated without incorporating ventilation information. In the functional plan, functional dose-volume constraints were applied in planning to spare the high-functional lung that was defined as the 90th percentile functional volume. The beam directions of the 2 plans were automatically determined by beam angle optimization.

RESULTS

The percentage of volume receiving a dose of ≥5 Gy (V5), V10, V20, and mean dose to the high-functional lung were 20.5%, 15.6%, 7.8%, and 4.6 Gy, respectively, for the anatomical plan, whereas they were 12.3%, 8.2%, 4.6%, and 3.2 Gy, respectively, for the functional plan. No significant differences in minimum dose, maximum dose, and conformity index of the planning target volume and in all dosimetric parameters for normal tissues between the anatomical and functional plans were seen.

CONCLUSIONS

We compared anatomical and functional plans for SBRT with 3-dimensional conformal radiation therapy for the first time. Our results demonstrated that a functional plan for SBRT reduced the dose in the high-functional regions without a significant change in the total lung or planning target volume even if the radiation technique cannot modulate beam intensity. Thus, this technique can be safely used in functional planning for SBRT.

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.

Review of automatic pulmonary lobe segmentation methods from CT

The computational detection of pulmonary lobes from CT images is a challenging segmentation problem with important respiratory health care applications, including surgical planning and regional image analysis. Several authors have proposed automated algorithms and we present a methodological review. These algorithms share a number of common stages and we consider each stage in turn, comparing the methods applied by each author and discussing their relative strengths. No standard method has yet emerged and none of the published methods have been demonstrated across a full range of clinical pathologies and imaging protocols. We discuss how improved methods could be developed by combining different approaches, and we use this to propose a workflow for the development of new algorithms.

Paired inspiratory/expiratory volumetric CT and deformable image registration for quantitative and qualitative evaluation of airflow limitation in smokers with or without copd

RATIONALE AND OBJECTIVES

To evaluate paired inspiratory/expiratory computed tomography (CT; iCT/eCT) and deformable image registration for quantitative and qualitative assessment of airflow limitation in smokers.

MATERIALS AND METHODS

Paired iCT/eCT images acquired from 35 smokers (30 men and 5 women) were coregistered and subtraction images (air trapping CT images [aCT]) generated. To evaluate emphysema quantitatively, the percentage of low-attenuation volume (LAV%) on iCT was calculated at -950 HU, as were mean and kurtosis on aCT for quantitative assessment of air trapping. Parametric response maps of emphysema (PRMe) and of functional small airways disease (PRMs) were also obtained. For qualitative evaluation of emphysema, low-attenuation areas on iCT were scored by consensus of two radiologists using Goddard classification. To assess air trapping qualitatively, the degree of air trapping on aCT was scored. For each quantitative and qualitative index, the Spearman rank correlation coefficient for forced expiratory flow in 1 second was calculated, and differences in correlation coefficients were statistically tested.

RESULTS

The correlation coefficients for the indices were as follows: mean on aCT, 0.800; kurtosis on aCT, -0.726; LAV%, -0.472; PRMe, -0.570; PRMs, -0.565; addition of PRMe and PRMs, -0.653; emphysema score, -0.502; air trapping score, -0.793. The indices showing significant differences were as follows: mean on aCT and addition of PRMe and PRMs (P = 1.43 × 10(-8)); air trapping score and emphysema score (P = .0169).

CONCLUSIONS

Air trapping images yielded more accurate quantitative and qualitative evaluation of airflow limitation than did LAV%, PRMe, PRMs, and Goddard classification.

Comparison between multivolume CT-based surrogates of regional ventilation in healthy subjects

RATIONALE AND OBJECTIVES

The assessment of regional ventilation is of critical importance when investigating lung function during disease progression and planning of pulmonary interventions. Recently, different computed tomography (CT)-based parameters have been proposed as surrogates of lung ventilation. The aim of the present study was to compare these parameters, namely variations of density (ΔHU), specific volume (sVol), and specific gas volume (ΔSVg) between different lung volumes, in relation to their topographic distribution within the lung.

MATERIALS AND METHODS

Ten healthy volunteers were scanned via high-resolution CT at residual volume (RV) and total lung capacity (TLC); ΔHU, sVol, and ΔSVg were mapped voxel by voxel after registering TLC onto RV. Variations of the three parameters along the vertical and horizontal directions were analyzed.

RESULTS

Along the vertical direction (from ventral to dorsal regions), a strong dependence on gravity was found in ΔHU and sVol, with greater values in the dorsal regions of the lung (P < .001), whereas ΔSVg was more homogeneously distributed within the lung. Conversely, along the caudocranial direction (from lung bases to apexes) where no gravitational gradient is present, the three parameters behaved similarly, with lower values at the apices.

CONCLUSIONS

ΔHU, sVol, and ΔSVg behave differently along the gravity direction. As the greater amount of air delivered to the dependent portion of the lung supplies a larger number of alveoli, the amount of gas delivered to alveoli compared to the mass of tissue is not gravity dependent. The minimization of gravity dependence in the distribution of ventilation when using ΔSVg suggests that this parameter is more reliable to discriminate healthy from pathologic regions.

Optimal beam arrangement for pulmonary ventilation image-guided intensity-modulated radiotherapy for lung cancer

Background

The principal aim of this study was to evaluate the feasibility of incorporating four-dimensional (4D)-computed tomography (CT)-based functional information into treatment planning and to evaluate the potential benefits of individualized beam setups to better protect lung functionality in patients with non-small cell lung cancer (NSCLC).

Methods

Peak-exhale and peak-inhale CT scans were carried out in 16 patients with NSCLC treated with intensity-modulated radiotherapy (IMRT). 4D-CT-based ventilation information was generated from the two sets of CT images using deformable image registration. Four kinds of IMRT plans were generated for each patient: two anatomic plans without incorporation of ventilation information, and two functional plans with ventilation information, using either five equally spaced beams (FESB) or five manually optimized beams (FMOB). The dosimetric parameters of the plans were compared in terms of target and normal tissue structures, with special focus on dose delivered to total lung and functional lung.

Results

In both the anatomic and functional plans, the percentages of both the functional and total lung regions irradiated at V₅, V₁₀, and V₂₀ (percentage volume irradiated to >5, >10 and >20 Gy, respectively) were significantly lower for FMOB compared with FESB (P < 0.05), but there was no significant difference for V₃₀ (P > 0.05). Compared with FESB, a greater degree of sparing of the functional lung was achieved in functional IMRT plans with optimal beam arrangement, without compromising target volume coverage or the irradiated volume of organs at risk, such as the spinal cord, esophagus, and heart.

Conclusions

Pulmonary ventilation image-guided IMRT planning with further optimization of beam arrangements improves the preservation of functional lung in patients with NSCLC.

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.

Evaluation of various deformable image registration algorithms for thoracic images

We evaluated the accuracy of one commercially available and three publicly available deformable image registration (DIR) algorithms for thoracic four-dimensional (4D) computed tomography (CT) images. Five patients with esophagus cancer were studied. Datasets of the five patients were provided by DIR-lab (dir-lab.com) and consisted of thoracic 4D CT images and a coordinate list of anatomical landmarks that had been manually identified. Expert landmark correspondence was used for evaluating DIR spatial accuracy. First, the manually measured displacement vector field (mDVF) was obtained from the coordinate list of anatomical landmarks. Then the automatically calculated displacement vector field (aDVF) was calculated by using the following four DIR algorithms: B-spine implemented in Velocity AI (Velocity Medical, Atlanta, GA, USA), free-form deformation (FFD), Horn-Schunk optical flow (OF) and Demons in DIRART of MATLAB software. Registration error is defined as the difference between mDVF and aDVF. The mean 3D registration errors were 2.7 ± 0.8 mm for B-spline, 3.6 ± 1.0 mm for FFD, 2.4 ± 0.9 mm for OF and 2.4 ± 1.2 mm for Demons. The results showed that reasonable accuracy was achieved in B-spline, OF and Demons, and that these algorithms have the potential to be used for 4D dose calculation, automatic image segmentation and 4D CT ventilation imaging in patients with thoracic cancer. However, for all algorithms, the accuracy might be improved by using the optimized parameter setting. Furthermore, for B-spline in Velocity AI, the 3D registration error was small with displacements of less than ∼10 mm, indicating that this software may be useful in this range of displacements.