Quantification of Local Metabolic Tumor Volume Changes by Registering Blended PET-CT Images for Prediction of Pathologic Tumor Response

Deformable registration with intensity correction for CESM monitoring response to Neoadjuvant Chemotherapy

This paper proposes a robust longitudinal registration method for Contrast Enhanced Spectral Mammography in monitoring neoadjuvant chemotherapy. Because breast texture intensity changes with the treatment, a non-rigid registration procedure with local intensity compensations is developed. The approach allows registering the low energy images of the exams acquired before and after the chemotherapy. The measured motion is then applied to the corresponding recombined images. The difference of registered images, called residual, makes vanishing the breast texture that did not changed between the two exams. Consequently, this registered residual allows identifying local density and iodine changes, especially in the lesion area. The method is validated with a synthetic NAC case where ground truths are available. Then the procedure is applied to 51 patients with 208 CESM image pairs acquired before and after the chemotherapy treatment. The proposed registration converged in all 208 cases. The intensity-compensated registration approach is evaluated with different mathematical metrics and through the repositioning of clinical landmarks (RMSE: 5.9 mm) and outperforms state-of-the-art registration techniques.

Inter- and intrafraction motion assessment and accumulated dose quantification of upper gastrointestinal organs during magnetic resonance-guided ablative radiation therapy of pancreas patients

Background and purpose

Stereotactic body radiation therapy (SBRT) of locally advanced pancreatic cancer (LAPC) is challenging due to significant motion of gastrointestinal (GI) organs. The goal of our study was to quantify inter and intrafraction deformations and dose accumulation of upper GI organs in LAPC patients.

Materials and methods

Five LAPC patients undergoing five-fraction magnetic resonance-guided radiation therapy (MRgRT) using abdominal compression and daily online plan adaptation to 50 Gy were analyzed. A pre-treatment, verification, and post-treatment MR imaging (MRI) for each of the five fractions (75 total) were used to calculate intra and interfraction motion. The MRIs were registered using Large Deformation Diffeomorphic Metric Mapping (LDDMM) deformable image registration (DIR) method and total dose delivered to stomach_duodenum, small bowel (SB) and large bowel (LB) were accumulated. Deformations were quantified using gradient magnitude and Jacobian integral of the Deformation Vector Fields (DVF). Registration DVFs were geometrically assessed using Dice and 95th percentile Hausdorff distance (HD95) between the deformed and physician’s contours. Accumulated doses were then calculated from the DVFs.

Results

Median Dice and HD95 were: Stomach_duodenum (0.9, 1.0 mm), SB (0.9, 3.6 mm), and LB (0.9, 2.0 mm). Median (max) interfraction deformation for stomach_duodenum, SB and LB was 6.4 (25.8) mm, 7.9 (40.5) mm and 7.6 (35.9) mm. Median intrafraction deformation was 5.5 (22.6) mm, 8.2 (37.8) mm and 7.2 (26.5) mm. Accumulated doses for two patients exceeded institutional constraints for stomach_duodenum, one of whom experienced Grade1 acute and late abdominal toxicity.

Conclusion

LDDMM method indicates feasibility to measure large GI motion and accumulate dose. Further validation on larger cohort will allow quantitative dose accumulation to more reliably optimize online MRgRT.

Multi-level multi-modality (PET and CT) fusion radiomics: prognostic modeling for non-small cell lung carcinoma

We developed multi-modality radiomic models by integrating information extracted from¹⁸F-FDG PET and CT images using feature- and image-level fusions, toward improved prognosis for non-small cell lung carcinoma (NSCLC) patients. Two independent cohorts of NSCLC patients from two institutions (87 and 95 patients) were cycled as training and testing datasets. Fusion approaches were applied at two levels, namely feature- and image-levels. For feature-level fusion, radiomic features were extracted individually from CT and PET images and concatenated. Alternatively, radiomic features extracted separately from CT and PET images were averaged. For image-level fusion, wavelet fusion was utilized and tuned with two parameters, namely CT weight and Wavelet Band Pass Filtering Ratio. Clinical and combined clinical + radiomic models were developed. Gray level discretization was performed at 3 different levels (16, 32 and 64) and 225 radiomics features were extracted. Overall survival (OS) was considered as the endpoint. For feature reduction, correlated (redundant) features were excluded using Spearman's correlation, and best combination of top ten features with highest concordance-indices (via univariate Cox model) were selected in each model for further multivariate Cox model. Moreover, prognostic score's median, obtained from the training cohort, was used intact in the testing cohort as a threshold to classify patients into low- versus high-risk groups, and log-rank test was applied to assess differences between the Kaplan-Meier curves. Overall, while models based on feature-level fusion strategy showed limited superiority over single-modalities, image-level fusion strategy significantly outperformed both single-modality and feature-level fusion strategies. As such, the clinical model (C-index = 0.656) outperformed all models from single-modality and feature-level strategies, but was outperformed by certain models from image-level fusion strategy. Our findings indicated that image-level fusion multi-modality radiomics models outperformed single-modality, feature-level fusion, and clinical models for OS prediction of NSCLC patients.

Generalizable cone beam CT esophagus segmentation using physics-based data augmentation

Automated segmentation of the esophagus is critical in image-guided/adaptive radiotherapy of lung cancer to minimize radiation-induced toxicities such as acute esophagitis. We have developed a semantic physics-based data augmentation method for segmenting the esophagus in both planning CT (pCT) and cone beam CT (CBCT) using 3D convolutional neural networks. One hundred and ninety-one cases with their pCTs and CBCTs from four independent datasets were used to train a modified 3D U-Net architecture and a multi-objective loss function specifically designed for soft-tissue organs such as the esophagus. Scatter artifacts and noises were extracted from week-1 CBCTs using a power-law adaptive histogram equalization method and induced to the corresponding pCT were reconstructed using CBCT reconstruction parameters. Moreover, we leveraged physics-based artifact induction in pCTs to drive the esophagus segmentation in real weekly CBCTs. Segmentations were evaluated using the geometric Dice coefficient and Hausdorff distance as well as dosimetrically using mean esophagus dose and D _5cc. Due to the physics-based data augmentation, our model trained just on the synthetic CBCTs was robust and generalizable enough to also produce state-of-the-art results on the pCTs and CBCTs, achieving Dice overlaps of 0.81 and 0.74, respectively. It is concluded that our physics-based data augmentation spans the realistic noise/artifact spectrum across patient CBCT/pCT data and can generalize well across modalities, eventually improving the accuracy of treatment setup and response analysis.

Quantification of accumulated dose and associated anatomical changes of esophagus using weekly Magnetic Resonance Imaging acquired during radiotherapy of locally advanced lung cancer

•

MRI is suited for tracking volumetric changes/accumulating doses in the esophagus.

•

Introduced medial axis of esophagus to calculate inter-fraction positional uncertainty.

•

Planned and accumulated esophagus dose-volume parameter differences are significant.

•

Longitudinal expansion of esophagus may link to acute esophagitis.

Background and purpose

Minimizing acute esophagitis (AE) in locally advanced non-small cell lung cancer (LA-NSCLC) is critical given the proximity between the esophagus and the tumor. In this pilot study, we developed a clinical platform for quantification of accumulated doses and volumetric changes of esophagus via weekly Magnetic Resonance Imaging (MRI) for adaptive radiotherapy (RT).

Material and methods

Eleven patients treated via intensity-modulated RT to 60–70 Gy in 2–3 Gy-fractions with concurrent chemotherapy underwent weekly MRIs. Eight patients developed AE grade 2 (AE2), 3–6 weeks after RT started. First, weekly MRI esophagus contours were rigidly propagated to planning CT and the distances between the medial esophageal axes were calculated as positional uncertainties. Then, the weekly MRI were deformably registered to the planning CT and the total dose delivered to esophagus was accumulated. Weekly Maximum Esophagus Expansion (MEex) was calculated using the Jacobian map. Eventually, esophageal dose parameters (Mean Esophagus Dose (MED), V_90% and D_5cc) between the planned and accumulated dose were compared.

Results

Positional esophagus uncertainties were 6.8 ± 1.8 mm across patients. For the entire cohort at the end of RT: the median accumulated MED was significantly higher than the planned dose (24 Gy vs. 21 Gy p = 0.006). The median V_90% and D_5cc were 12.5 cm³ vs. 11.5 cm³ (p = 0.05) and 61 Gy vs. 60 Gy (p = 0.01), for accumulated and planned dose, respectively. The median MEex was 24% and was significantly associated with AE2 (p = 0.008).

Conclusions

MRI is well suited for tracking esophagus volumetric changes and accumulating doses. Longitudinal esophagus expansion could reflect radiation-induced inflammation that may link to AE.

LDeform: Longitudinal deformation analysis for adaptive radiotherapy of lung cancer

PURPOSE

Conventional radiotherapy for large lung tumors is given over several weeks, during which the tumor typically regresses in a highly nonuniform and variable manner. Adaptive radiotherapy would ideally follow these shape changes, but we need an accurate method to extrapolate tumor shape changes. We propose a computationally efficient algorithm to quantitate tumor surface shape changes that makes minimal assumptions, identifies fixed points, and can be used to predict future tumor geometrical response.

METHODS

A novel combination of nonrigid iterative closest point (ICP) and local shape-preserving map algorithms, LDeform, is developed to enable visualization, prediction, and categorization of both diffeomorphic and nondiffeomorphic tumor deformations during an extended course of radiotherapy.

RESULTS

We tested and validated our technique on 31 longitudinal CT/MRI subjects, with five to nine time points each. Based on this tumor deformation analysis, regions of local growth, shrinkage, and anchoring are identified and tracked across multiple time points. This categorization in turn represents a rational biomarker of local response. Results demonstrate useful predictive power, with an averaged Dice coefficient and surface mean-squared error of 0.85 and 2.8 mm, respectively, over all images.

CONCLUSIONS

We conclude that the LDeform algorithm can facilitate the adaptive decision-making process during lung cancer radiotherapy.

Multi-Level Multi-Modality Fusion Radiomics: Application to PET and CT Imaging for Prognostication of Head and Neck Cancer

To characterize intra-tumor heterogeneity comprehensively, we propose a multi-level fusion strategy to combine PET and CT information at the image-, matrix-and feature-levels towards improved prognosis. Specifically, we developed fusion radiomics in the context of 3 prognostic outcomes in a multi-center setting (4 centers) involving 296 head & neck cancer patients. Eight clinical parameters were first utilized to build a (1) clinical model. We also built models by extracting 127 radiomics features from (2) PET images alone; (3-8) PET and CT images fused via wavelet-based fusion (WF) using CT-weights of 0.2, 0.4, 0.6 and 0.8, gradient transfer fusion (GTF), and guided filtering-based fusion (GFF); (9) fused matrices (sumMat); (10-11) fused features constructed via feature averaging (avgFea) and feature concatenation (conFea); and finally, (12) CT images alone; above models were also expanded to include both clinical and radiomics features. Seven variations of training and testing partitions were investigated. Highest performance in 5, 6 and 5 partitions was achieved by image-level fusion strategies for RFS, MFS and OS prediction, respectively. Among all partitions, WF0.6 and WF0.8 showed significantly higher performance than CT model for RFS (C-index: 0.60 ± 0.04 vs. 0.56 ± 0.03, p-value: 0.015) and MFS (C-index: 0.71 ± 0.13 vs. 0.62 ± 0.08, p-value: 0.020) predictions, respectively. In partition CER 23 vs. 14, WF0.6 significantly outperformed Clinical model for RFS prediction (C-index: 0.67 vs. 0.53, p-value: 0.003); both avgFea and WF0.6 showed C-index of 0.64 and significantly higher than that of PET only (C-index: 0.51, p-value: 0.018 and 0.031, respectively) for OS prediction. Fusion radiomics modeling showed varying improvements compared to single modality models for different outcome predictions in different partitions, highlighting the importance of generalizing radiomics models. Image-level fusion holds potential to capture more useful characteristics.