Pretreatment CT and PET radiomics predicting rectal cancer patients in response to neoadjuvant chemoradiotherapy.

BACKGROUND: The purpose of this study was to characterize pre-treatment non-contrast computed tomography (CT) and 18F-fluorodeoxyglucose positron emission tomography (PET) based radiomics signatures predictive of pathological response and clinical outcomes in rectal cancer patients treated with neoadjuvant chemoradiotherapy (NACR T). MATERIALS AND METHODS: An exploratory analysis was performed using pre-treatment non-contrast CT and PET imaging dataset. The association of tumor regression grade (TRG) and neoadjuvant rectal (NAR) score with pre-treatment CT and PET features was assessed using machine learning algorithms. Three separate predictive models were built for composite features from CT + PET. RESULTS: The patterns of pathological response were TRG 0 (n = 13; 19.7%), 1 (n = 34; 51.5%), 2 (n = 16; 24.2%), and 3 (n = 3; 4.5%). There were 20 (30.3%) patients with low, 22 (33.3%) with intermediate and 24 (36.4%) with high NAR scores. Three separate predictive models were built for composite features from CT + PET and analyzed separately for clinical endpoints. Composite features with α = 0.2 resulted in the best predictive power using logistic regression. For pathological response prediction, the signature resulted in 88.1% accuracy in predicting TRG 0 vs. TRG 1-3; 91% accuracy in predicting TRG 0-1 vs. TRG 2-3. For the surrogate of DFS and OS, it resulted in 67.7% accuracy in predicting low vs. intermediate vs. high NAR scores. CONCLUSION: The pre-treatment composite radiomics signatures were highly predictive of pathological response in rectal cancer treated with NACR T. A larger cohort is warranted for further validation.

Superficial and peripheral dose in compensator-based FFF beam IMRT.

Flattening filter-free (FFF) beams produce higher dose rates. Combined with compensator-based intensity modulated radiotherapy (IMRT) techniques, the dose delivery for each beam can be much shorter compared to the flattened beam MLC-based or flattened beam compensator-based IMRT. This 'snap shot' IMRT delivery is beneficial to patients for tumor motion management. Due to softer energy, superficial doses in FFF beam treatment are usually higher than those from flattened beams. Due to no flattening filter, thus less photon scattering, peripheral doses are usually lower in FFF beam treatment. However, in compensator-based IMRT using FFF beams, the compensator is in the beam pathway. Does it introduce beam hardening effects and scattering such that the superficial dose is lower and peripheral dose is higher compared to FFF beam MLC-based IMRT? This study applied Monte Carlo techniques to investigate the superficial and peripheral doses in compensator-based IMRT using FFF beams and compared it to the MLC-based IMRT using FFF beams and flattened beams. Besides varying thicknesses of brass slabs to simulate varying thicknesses of compensators, a simple cone-shaped compensator was simulated to mimic a clinical application. The dose distribution in water phantom by the cone-shaped compensator was then simulated by multiple MLC-defined FFF and flattened beams with varying apertures. After normalization to the maximum dose, Dmax, the superficial and peripheral doses were compared between the FFF beam compensator-based IMRT and FFF/flattened beam MLC-based IMRT. The superficial dose at the central 0.5 mm depth was about 1% (of Dmax) lower in the compensator-based 6 MV FFF (6FFF) IMRT compared to the MLC-based 6FFF IMRT, and about 8% higher than the flattened 6 MV MLC-based IMRT dose. At 8 cm off-axis at depth of central maximum dose, dmax, the peripheral dose between the 6FFF and flattened 6 MV MLC demonstrated similar doses, while the compensator dose was about 1% (of Dmax) higher. Compensators reduce the superficial doses slightly compared to open FFF beams, but increases the peripheral doses due to scatter in the compensator.

Reproducibility of F18-FDG PET radiomic features for different cervical tumor segmentation methods, gray-level discretization, and reconstruction algorithms.

Site-specific investigations of the role of radiomics in cancer diagnosis and therapy are emerging. We evaluated the reproducibility of radiomic features extracted from 18 Flourine-fluorodeoxyglucose (18 F-FDG) PET images for three parameters: manual versus computer-aided segmentation methods, gray-level discretization, and PET image reconstruction algorithms. Our cohort consisted of pretreatment PET/CT scans from 88 cervical cancer patients. Two board-certified radiation oncologists manually segmented the metabolic tumor volume (MTV1 and MTV2 ) for each patient. For comparison, we used a graphical-based method to generate semiautomated segmented volumes (GBSV). To address any perturbations in radiomic feature values, we down-sampled the tumor volumes into three gray-levels: 32, 64, and 128 from the original gray-level of 256. Finally, we analyzed the effect on radiomic features on PET images of eight patients due to four PET 3D-reconstruction algorithms: maximum likelihood-ordered subset expectation maximization (OSEM) iterative reconstruction (IR) method, fourier rebinning-ML-OSEM (FOREIR), FORE-filtered back projection (FOREFBP), and 3D-Reprojection (3DRP) analytical method. We extracted 79 features from all segmentation method, gray-levels of down-sampled volumes, and PET reconstruction algorithms. The features were extracted using gray-level co-occurrence matrices (GLCM), gray-level size zone matrices (GLSZM), gray-level run-length matrices (GLRLM), neighborhood gray-tone difference matrices (NGTDM), shape-based features (SF), and intensity histogram features (IHF). We computed the Dice coefficient between each MTV and GBSV to measure segmentation accuracy. Coefficient values close to one indicate high agreement, and values close to zero indicate low agreement. We evaluated the effect on radiomic features by calculating the mean percentage differences (d¯) between feature values measured from each pair of parameter elements (i.e. segmentation methods: MTV1 -MTV2 , MTV1 -GBSV, MTV2 -GBSV; gray-levels: 64-32, 64-128, and 64-256; reconstruction algorithms: OSEM-FORE-OSEM, OSEM-FOREFBP, and OSEM-3DRP). We used |d¯| as a measure of radiomic feature reproducibility level, where any feature scored |d¯| ±SD ≤ |25|% ± 35% was considered reproducible. We used Bland-Altman analysis to evaluate the mean, standard deviation (SD), and upper/lower reproducibility limits (U/LRL) for radiomic features in response to variation in each testing parameter. Furthermore, we proposed U/LRL as a method to classify the level of reproducibility: High- ±1% ≤ U/LRL ≤ ±30%; Intermediate- ±30% < U/LRL ≤ ±45%; Low- ±45 < U/LRL ≤ ±50%. We considered any feature below the low level as nonreproducible (NR). Finally, we calculated the interclass correlation coefficient (ICC) to evaluate the reliability of radiomic feature measurements for each parameter. The segmented volumes of 65 patients (81.3%) scored Dice coefficient >0.75 for all three volumes. The result outcomes revealed a tendency of higher radiomic feature reproducibility among segmentation pair MTV1 -GBSV than MTV2 -GBSV, gray-level pairs of 64-32 and 64-128 than 64-256, and reconstruction algorithm pairs of OSEM-FOREIR and OSEM-FOREFBP than OSEM-3DRP. Although the choice of cervical tumor segmentation method, gray-level value, and reconstruction algorithm may affect radiomic features, some features were characterized by high reproducibility through all testing parameters. The number of radiomic features that showed insensitivity to variations in segmentation methods, gray-level discretization, and reconstruction algorithms was 10 (13%), 4 (5%), and 1 (1%), respectively. These results suggest that a careful analysis of the effects of these parameters is essential prior to any radiomics clinical application.

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

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

Is wax equivalent to tissue in electron conformal therapy planning? A Monte Carlo study of material approximation introduced dose difference.

With CT-based Monte Carlo (MC) dose calculations, material composition is often assigned based on the standard Hounsfield unit ranges. This is known as the density threshold method. In bolus electron conformal therapy (BolusECT), the bolus material, machineable wax, would be assigned as soft tissue and the electron density is assumed equivalent to soft tissue based on its Hounsfield unit. This study investigates the dose errors introduced by this material assignment. BEAMnrc was used to simulate electron beams from a Trilogy accelerator. SPRRZnrc was used to calculate stopping power ratios (SPR) of tissue to wax, SPR (tissue) (wax), and tissue to water, SPR(tissue) (water), for 6, 9, 12, 15, and 18 MeV electron beams, of which 12 and 15MeV beams are the most commonly used energies in BolusECT. DOSXYZnrc was applied in dose distribution calculations in a tissue phantom with either flat wax slabs of various thicknesses or a wedge-shaped bolus on top. Dose distribution for two clinical cases, a chest wall and a head and neck, were compared with the bolus material treated as wax or tissue. The SPR(tissue) (wax) values for 12 and 15MeV beams are between 0.935 and 0.945, while the SPR(tissue) (water) values are between 0.990 and 0.991. For a 12 MeV beam, the dose in tissue immediately under the bolus is overestimated by 2.5% for a 3 cm bolus thickness if the wax bolus is treated as tissue. For 15 MeV beams, the error is 1.4%. However, in both clinical cases the differences in the PTV DVH is negligible. Due to stopping power differences, dose differences of up to 2.5% are observed in MC simulations if the bolus material is misassigned as tissue in BolusECT dose calculations. However, for boluses thinner than 2 cm that are more likely encountered in practice, the error is within clinical tolerance.

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

Ventilation imaging using 4D CT is a convenient and low-cost functional imaging methodology which might be of value in radiotherapy treatment planning to spare functional lung volumes. Deformable image registration (DIR) is needed to calculate ventilation imaging from 4D CT. This study investigates the dependence of calculated ventilation on DIR methods and ventilation algorithms. DIR of the normal end expiration and normal end inspiration phases of the 4D CT images was used to correlate the voxels between the two respiratory phases. Three different DIR algorithms, optical flow (OF), diffeomorphic demons (DD), and diffeomorphic morphons (DM) were retrospectively applied to ten esophagus and ten lung cancer cases with 4D CT image sets that encompassed the entire lung volume. The three ventilation extraction methods were used based on either the Jacobian, the change in volume of the voxel, or directly calculated from Hounsfield units. The ventilation calculation algorithms used are the Jacobian, ΔV, and HU method. They were compared using the Dice similarity coefficient (DSC) index and Bland-Altman plots. Dependence of ventilation images on the DIR was greater for the ΔV and the Jacobian methods than for the HU method. The DSC index for 20% of low-ventilation volume for ΔV was 0.33 ± 0.03 (1 SD) between OF and DM, 0.44 ± 0.05 between OF and DD, and 0.51 ± 0.04 between DM and DD. The similarity comparisons for Jacobian were 0.32 ± 0.03, 0.44 ± 0.05, and 0.51 ± 0.04, respectively, and for HU they were 0.53 ± 0.03, 0.56 ± 0.03, and 0.76 ± 0.04, respectively. Dependence of extracted ventilation on the ventilation algorithm used showed good agreement between the ΔV and Jacobian methods, but differed significantly for the HU method. DSC index for using OF as DIR was 0.86 ± 0.01 between ΔV and Jacobian, 0.28 ± 0.04 between ΔV and HU, and 0.28 ± 0.04 between Jacobian and HU, respectively. When using DM or DD as DIR, similar values were obtained when comparing the different ventilation calculation methods. The similarity values for the 20% high-ventilation volume were close to those found for the 20% low-ventilation volume. The results obtained with DSC index were confirmed when using the Bland-Altman plots for comparing the ventilation images. Our data suggest that ventilation calculated from 4D CT depends on the DIR algorithm employed. Similarities between ΔV and Jacobian are higher than between ΔV and HU, and Jacobian and HU.

Comparing dose in the build-up region between compensator- and MLC-based IMRT.

The build-up dose in the megavoltage photon beams can be a limiting factor in intensity-modulated radiation therapy (IMRT) treatments. Excessive surface dose can cause patient discomfort and treatment interruptions, while underdosing may lead to tumor repopulation and local failure. Dose in the build-up region was investigated for IMRT delivery with solid brass compensator technique(compensator-based IMRT) and compared with that of multileaf collimator (MLC)-based IMRT. A Varian Trilogy linear accelerator equipped with an MLC was used for beam delivery. A special solid brass step-wise compensator was designed and built for testing purposes. Two step-and-shoot MLC fields were programmed to produce a similar modulated step-wise dose profile. The MLC and compensator dose profiles were measured and adjusted to match at the isocenter depth of 10 cm. Build-up dose in the 1-5 mm depth range was measured with an ultrathin window, fixed volume parallel plate ionization chamber. Monte Carlo simulations were used to model the brass compensator and step-and-shoot MLC fields. The measured and simulated profiles for the two IMRT techniques were matched at the isocenter depth of 10 cm. Different component contributions to the shallow dose, including the MLC scatter, were quantified. Mean spectral energies for the open and filtered beams were calculated. The compensator and MLC profiles at 10 cm depth were matched better than ± 1.5%. The build-up dose was up to 7% lower for compensator IMRT compared to MLC IMRT due to beam hardening in the brass. Low-energy electrons contribute 22% and 15% dose at 1 mm depth for compensator and MLC modalities, respectively. Compensator-based IMRT delivers less dose in the build-up region than MLC-based IMRT does, even though a compensator is closer to the skin than the MLC.

Providing a fast conversion of total dose to biological effective dose (BED) for hybrid seed brachytherapy.

Optimization of permanent seed implant brachytherapy plans for treatment of prostate cancer should be based on biological effective dose (BED) distributions, since dose does not accurately represent biological effects between different types of sources. Currently, biological optimization for these plans is not feasible due to the amount of time necessary to calculate the BED distribution. This study provides a fast calculation method, based on the total dose, to calculate the BED distribution. Distributions of various numbers of hybrid seeds were used to calculate total dose distributions, as well as BED distributions. Hybrid seeds are a mixture of different isotopes (in this study (125)I and (103)Pd). Three ratios of hybrid seeds were investigated: 25/75, 50/50, and 75/25. The total dose and BED value from each voxel were coupled together to produce graphs of total dose vs. BED. Equations were then derived from these graphs. The study investigated four types of tissue: bladder, rectum, prostate, and other normal tissue. Equations were derived from the total dose - BED correspondence. Accuracy of conversion from total dose to BED was within 2 Gy; however, accuracy of conversion was found to be better for high total dose regions as compared to lower dose regions. The method introduced in this paper allows one to perform fast conversion of total dose to BED for brachytherapy using hybrid seeds, which makes the BED-based plan optimization practical. The method defined here can be extended to other ratios, as well as other tissues that are affected by permanent seed implant brachytherapy (i.e., breast).

Pretreatment CT and 18 F-FDG PET-based radiomic model predicting pathological complete response and loco-regional control following neoadjuvant chemoradiation in oesophageal cancer.

INTRODUCTION: To develop a radiomic-based model to predict pathological complete response (pCR) and outcome following neoadjuvant chemoradiotherapy (NACRT) in oesophageal cancer. METHODS: We analysed 68 patients with oesophageal cancer treated with NACRT followed by esophagectomy, who had staging 18F-fluorodeoxyglucose (18 F-FDG) positron emission tomography (PET) and computed tomography (CT) scans performed at our institution. An in-house data-chjmirocterization algorithm was used to extract 3D-radiomic features from the segmented primary disease. Prediction models were constructed and internally validated. Composite feature, Fc  = α * FPET  + (1 - α) * FCT , 0 ≤ α ≤ 1, was constructed for each corresponding CT and PET feature. Loco-regional control (LRC), recurrence-free survival (RFS), metastasis-free survival (MFS) and overall survival (OS) were estimated by Kaplan-Meier analysis, and compared using log-rank test. RESULTS: Median follow-up was 59 months. pCR was achieved in 34 (50%) patients. Five-year RFS, LRC, MFS and OS were 67.1%, 88.5%, 75.6% and 57.6%, respectively. Tumour Regression Grade (TRG) 0-1 indicative of complete response or minimal residual disease was significantly associated with improved 5-year LRC [93.7% vs 71.8%; P = 0.020; HR 0.19, 95% CI 0.04-0.85]. Four sepjmirote pCR predictive models were built for CT alone, PET alone, CT+PET and composite. CT, PET and CT+PET models had AUC 0.73 ± 0.08, 0.66 ± 0.08 and 0.77 ± 0.07, respectively. The composite model resulted in an improvement of pCR predicting power with AUC 0.87 ± 0.06. Stratifying patients with a low versus high radiomic score showed clinically relevant improvement in 5-year LRC favouring low-score group (91.1% vs. 80%, 95% CI 0.09-1.77, P = 0.2). CONCLUSION: The composite CT/PET radiomics model was highly predictive of pCR following NACRT. Validation in larger data sets is warranted to determine whether the model can predict clinical outcomes.

CT-based radiomic features to predict pathological response in rectal cancer: A retrospective cohort study.

INTRODUCTION: Innovative biomarkers to predict treatment response in rectal cancer would be helpful in optimizing personalized treatment approaches. In this study, we aimed to develop and validate a CT-based radiomic imaging biomarker to predict pathological response. METHODS: We used two independent cohorts of rectal cancer patients to develop and validate a CT-based radiomic imaging biomarker predictive of treatment response. A total of 91 rectal cancer cases treated from 2009 to 2018 were assessed for the tumour regression grade (TRG) (0 = pathological complete response, pCR; 1 = moderate response; 2 = partial response; 3 = poor response). Exploratory analysis was performed by combining pre-treatment non-contrast CT images and patterns of TRG. The models built from the training cohort were further assessed using the independent validation cohort. RESULTS: The patterns of pathological response in training and validation groups were TRG 0 (n = 14, 23.3%; n = 6, 19.4%), 1 (n = 31, 51.7%; n = 15, 48.4%), 2 (n = 12, 20.0%; n = 7, 22.6%) and 3 (n = 3, 5.0%; n = 3, 9.7%), respectively. Separate predictive models were built and analysed from CT features for pathological response. For pathological response prediction, the model including 8 radiomic features by random forest method resulted in 83.9% accuracy in predicting TRG 0 vs TRG 1-3 in validation. CONCLUSION: The pre-treatment CT-based radiomic signatures were developed and validated in two independent cohorts. This imaging biomarker provided a promising way to predict pCR and select patients for non-operative management.

Dose mapping: validation in 4D dosimetry with measurements and application in radiotherapy follow-up evaluation.

The purpose of this paper is to validate a dose mapping program using optical flow method (OFM), and to demonstrate application of the program in radiotherapy follow-up evaluation. For the purpose of validation, the deformation matrices between four-dimensional (4D) CT data of different simulated respiration phases of a phantom were calculated using OFM. The matrices were then used to map doses of all phases to a single-phase image, and summed in equal time weighting. The calculated dose should closely represent the dose delivered to the moving phantom if the deformation matrices are accurately calculated. The measured point doses agreed with the OFM calculations better than 2% at isocenters, and dose distributions better than 1mm for the 50% isodose line. To demonstrate proof-of-concept for the use of deformable image registration in dose mapping for treatment evaluation, the treatment-planning CT was registered with the post-treatment CT image from the positron emission tomography (PET)/CT resulting in a deformation matrix. The dose distribution from the treatment plan was then mapped onto the restaging PET/CT using the deformation matrix. Two cases in which patients had thoracic malignancies are presented. Each patient had CT-based treatment planning for radiotherapy and restaging fluorodeoxy glucose (FDG)-PET/CT imaging 4-6 weeks after completion of treatments. Areas of pneumonitis and recurrence were identified radiographically on both PET and CT restaging images. Local dose and standard uptake values for pneumonitis and recurrence were studied as a demonstration of this method. By comparing the deformable mapped dose to measurement, the treatment evaluation method which is introduced in this manuscript proved to be accurate. It thus provides a more accurate analysis than other rigid or linear dose-image registration when used in studying treatment outcome versus dose.

Accounting for reconstruction kernel-induced variability in CT radiomic features using noise power spectra.

Large variability in computed tomography (CT) radiomics feature values due to CT imaging parameters can have subsequent implications on the prognostic or predictive significance of these features. Here, we investigated the impact of pitch, dose, and reconstruction kernel on CT radiomic features. Moreover, we introduced correction factors to reduce feature variability introduced by reconstruction kernels. The credence cartridge radiomics and American College of Radiology (ACR) phantoms were scanned on five different scanners. ACR phantom was used for 3-D noise power spectrum (NPS) measurements to quantify correlated noise. The coefficient of variation (COV) was used as the variability assessment metric. The variability in texture features due to different kernels was reduced by applying the NPS peak frequency and region of interest (ROI) maximum intensity as correction factors. Most texture features were dose independent but were strongly kernel dependent, which is demonstrated by a significant shift in NPS peak frequency among kernels. Percentage improvement in robustness was calculated for each feature from original and corrected %COV values. Percentage improvements in robustness of 19 features were in the range of 30% to 78% after corrections. We show that NPS peak frequency and ROI maximum intensity can be used as correction factors to reduce variability in CT texture feature values due to reconstruction kernels.

Investigating multi-radiomic models for enhancing prediction power of cervical cancer treatment outcomes.

Quantitative image features, also known as radiomic features, have shown potential for predicting treatment outcomes in several body sites. We quantitatively analyzed 18Fluorine-fluorodeoxyglucose (18F-FDG) Positron Emission Tomography (PET) uptake heterogeneity in the Metabolic Tumor Volume (MTV) of eighty cervical cancer patients to investigate the predictive performance of radiomic features for two treatment outcomes: the development of distant metastases (DM) and loco-regional recurrent disease (LRR). We aimed to fit the highest predictive features in multiple logistic regression models (MLRs). To generate such models, we applied backward feature selection method as part of Leave-One-Out Cross Validation (LOOCV) within a training set consisting of 70% of the original patient cohort. The trained MLRs were tested on an independent set consisted of 30% of the original cohort. We evaluated the performance of the final models using the Area under the Receiver Operator Characteristic Curve (AUC). Accordingly, six models demonstrated superior predictive performance for both outcomes (four for DM and two for LRR) when compared to both univariate-radiomic feature models and Standard Uptake Value (SUV) measurements. This demonstrated approach suggests that the ability of the pre-radiochemotherapy PET radiomics to stratify patient risk for DM and LRR could potentially guide management decisions such as adjuvant systemic therapy or radiation dose escalation.

The clinical outcome correlations between radiation dose and pretreatment metabolic tumor volume for radiotherapy in head and neck cancer: A retrospective analysis.

This study was to investigate the clinical outcomes between radiation dose and pretreatment metabolic tumor volume (MTV) in patients with head and neck cancer treated with definitive chemoradiotherapy.Thirty-four patients received pretreatment F-fluorodeoxyglucose (F-FDG) positron emission tomography-computed tomography (PET/CT) were recruited for this study. The CT-based volume (gross tumor volume of the primary [GTVp]) and 4 types of MTVs were measured on the basis of either a maximal standardized uptake value (SUVmax) of 2.5 (MTV2.5), 3.0 (MTV3.0), or a fixed threshold of 40% (MTV40%), 50% (MTV50%). F-FDG PET-CT images before treatment, and data including response to treatment, local recurrence, death due to the cancer, disease-free survival (DFS) and primary relapse-free survival (PRFS), were collected for analysis.The Wilcoxon rank test showed that all values determined by the different delineation techniques were significantly different from the GTVp (P < .05). Tumor volume and the homogeneity of target dose of MTV2.5, MTV3.0, MTV40%, and MTV50% were significantly different between the 2 groups of patients through treatment outcomes (P < .05).The survival curves for DFS and PRFS demonstrated that the homogeneity of the target dose in MTVs was a good indicator. The homogeneity of target dose in the tumor is a potential indicator of DSF and PRFS in patients with head and neck cancer who underwent radiotherapy.

Ventilation Series Similarity: A Study for Ventilation Calculation Using Deformable Image Registration and 4DCT to Avoid Motion Artifacts.

The major problem with ventilation distribution calculations using DIR and 4DCT is the motion artifacts in 4DCT. Quite often not all phases would exhibit mushroom motion artifacts. If the ventilation series similarity is sufficiently robust, the ventilation distribution can be calculated using only the artifact-free phases. This study investigated the ventilation similarity among the data derived from different respiration phases. Fifteen lung cancer cases were analyzed. In each case, DIR was performed between the end-expiration phase and all other phases. Ventilation distributions were then calculated using the deformation matrices. The similarity was compared between the series ventilation distributions. The correlation between the majority phases was reasonably good, with average SCC values between 0.28 and 0.70 for the original data and 0.30 and 0.75 after smoothing. The better correlation between the neighboring phases, with average SCC values between 0.55 and 0.70 for the original data, revealed the nonlinear property of the dynamic ventilation. DSC analysis showed the same trend. To reduce the errors if motion artifacts are present, the phases without serious mushroom artifacts may be used. To minimize the effect of the nonlinearity in dynamic ventilation, the calculation phase should be chosen as close to the end-inspiration as possible.

Fiducial markers coupled with 3D PET/CT offer more accurate radiation treatment delivery for locally advanced esophageal cancer.

BACKGROUND AND AIMS: The role of three-dimensional positron emission tomography/computed tomography (3 D PET/CT) in esophageal tumors that move with respiration and have potential for significant mucosal inflammation is unclear. The aim of this study was to determine the correlation between gross tumor volumes derived from 3 D PET/CT and endoscopically placed fiducial markers. METHODS: This was a retrospective, IRB approved analysis of 40 patients with esophageal cancer with fiducials implanted and PET/CT. The centroid of each fiducial was identified on PET/CT images. Distance between tumor volume and fiducials was measured using axial slices. Image features were extracted and tested for pathologic response predictability. RESULTS: The median adaptively calculated threshold value of the standardized uptake value (SUV) to define the metabolic tumor volume (MTV) border was 2.50, which corresponded to a median 23 % of the maximum SUV. The median distance between the inferior fiducial centroid and MTV was - 0.60 cm (- 3.9 to 2.7 cm). The median distance between the superior fiducial centroid and MTV was 1.25 cm (- 4.2 to 6.9 cm). There was no correlation between MTV-to-fiducial distances greater than 2 cm and the gastroenterologist who performed the fiducial implantation. Eccentricity demonstrated statistically significant correlations with pathologic response. CONCLUSIONS: There was a stronger correlation between inferior fiducial location and MTV border compared to the superior extent. The etiology of the discordance superiorly is unclear, potentially representing benign secondary esophagitis, presence of malignant nodes, inflammation caused by technical aspects of the fiducial placement itself, or potential submucosal disease. Given the concordance inferiorly and the ability to more precisely set up the patient with daily image guidance matching to fiducials, it may be possible to minimize the planning tumor volume (PTV) margin in select patients, thereby, limiting dose to normal structures.

Intrinsic dependencies of CT radiomic features on voxel size and number of gray levels.

PURPOSE: Many radiomics features were originally developed for non-medical imaging applications and therefore original assumptions may need to be reexamined. In this study, we investigated the impact of slice thickness and pixel spacing (or pixel size) on radiomics features extracted from Computed Tomography (CT) phantom images acquired with different scanners as well as different acquisition and reconstruction parameters. The dependence of CT texture features on gray-level discretization was also evaluated. METHODS AND MATERIALS: A texture phantom composed of 10 different cartridges of different materials was scanned on eight different CT scanners from three different manufacturers. The images were reconstructed for various slice thicknesses. For each slice thickness, the reconstruction Field Of View (FOV) was varied to render pixel sizes ranging from 0.39 to 0.98 mm. A fixed spherical region of interest (ROI) was contoured on the images of the shredded rubber cartridge and the 3D printed, 20% fill, acrylonitrile butadiene styrene plastic cartridge (ABS20) for all phantom imaging sets. Radiomic features were extracted from the ROIs using an in-house program. Features categories were: shape (10), intensity (16), GLCM (24), GLZSM (11), GLRLM (11), and NGTDM (5), fractal dimensions (8) and first-order wavelets (128), for a total of 213 features. Voxel-size resampling was performed to investigate the usefulness of extracting features using a suitably chosen voxel size. Acquired phantom image sets were resampled to a voxel size of 1 × 1 × 2 mm3 using linear interpolation. Image features were therefore extracted from resampled and original datasets and the absolute value of the percent coefficient of variation (%COV) for each feature was calculated. Based on the %COV values, features were classified in 3 groups: (1) features with large variations before and after resampling (%COV >50); (2) features with diminished variation (%COV <30) after resampling; and (3) features that had originally moderate variation (%COV <50%) and were negligibly affected by resampling. Group 2 features were further studied by modifying feature definitions to include voxel size. Original and voxel-size normalized features were used for interscanner comparisons. A subsequent analysis investigated feature dependency on gray-level discretization by extracting 51 texture features from ROIs from each of the 10 different phantom cartridges using 16, 32, 64, 128, and 256 gray levels. RESULTS: Out of the 213 features extracted, 150 were reproducible across voxel sizes, 42 improved significantly (%COV <30, Group 2) after resampling, and 21 had large variations before and after resampling (Group 1). Ten features improved significantly after definition modification effectively removed their voxel-size dependency. Interscanner comparison indicated that feature variability among scanners nearly vanished for 8 of these 10 features. Furthermore, 17 out of 51 texture features were found to be dependent on the number of gray levels. These features were redefined to include the number of gray levels which greatly reduced this dependency. CONCLUSION: Voxel-size resampling is an appropriate pre-processing step for image datasets acquired with variable voxel sizes to obtain more reproducible CT features. We found that some of the radiomics features were voxel size and gray-level discretization-dependent. The introduction of normalizing factors in their definitions greatly reduced or removed these dependencies.

Variability of Image Features Computed from Conventional and Respiratory-Gated PET/CT Images of Lung Cancer.

Radiomics is being explored for potential applications in radiation therapy. How various imaging protocols affect quantitative image features is currently a highly active area of research. To assess the variability of image features derived from conventional [three-dimensional (3D)] and respiratory-gated (RG) positron emission tomography (PET)/computed tomography (CT) images of lung cancer patients, image features were computed from 23 lung cancer patients. Both protocols for each patient were acquired during the same imaging session. PET tumor volumes were segmented using an adaptive technique which accounted for background. CT tumor volumes were delineated with a commercial segmentation tool. Using RG PET images, the tumor center of mass motion, length, and rotation were calculated. Fifty-six image features were extracted from all images consisting of shape descriptors, first-order features, and second-order texture features. Overall, 26.6% and 26.2% of total features demonstrated less than 5% difference between 3D and RG protocols for CT and PET, respectively. Between 10 RG phases in PET, 53.4% of features demonstrated percent differences less than 5%. The features with least variability for PET were sphericity, spherical disproportion, entropy (first and second order), sum entropy, information measure of correlation 2, Short Run Emphasis (SRE), Long Run Emphasis (LRE), and Run Percentage (RPC); and those for CT were minimum intensity, mean intensity, Root Mean Square (RMS), Short Run Emphasis (SRE), and RPC. Quantitative analysis using a 3D acquisition versus RG acquisition (to reduce the effects of motion) provided notably different image feature values. This study suggests that the variability between 3D and RG features is mainly due to the impact of respiratory motion.

Monte Carlo study of radiation dose enhancement by gadolinium in megavoltage and high dose rate radiotherapy.

MRI is often used in tumor localization for radiotherapy treatment planning, with gadolinium (Gd)-containing materials often introduced as a contrast agent. Motexafin gadolinium is a novel radiosensitizer currently being studied in clinical trials. The nanoparticle technologies can target tumors with high concentration of high-Z materials. This Monte Carlo study is the first detailed quantitative investigation of high-Z material Gd-induced dose enhancement in megavoltage external beam photon therapy. BEAMnrc, a radiotherapy Monte Carlo simulation package, was used to calculate dose enhancement as a function of Gd concentration. Published phase space files for the TrueBeam flattening filter free (FFF) and conventional flattened 6MV photon beams were used. High dose rate (HDR) brachytherapy with Ir-192 source was also investigated as a reference. The energy spectra difference caused a dose enhancement difference between the two beams. Since the Ir-192 photons have lower energy yet, the photoelectric effect in the presence of Gd leads to even higher dose enhancement in HDR. At depth of 1.8 cm, the percent mean dose enhancement for the FFF beam was 0.38±0.12, 1.39±0.21, 2.51±0.34, 3.59±0.26, and 4.59±0.34 for Gd concentrations of 1, 5, 10, 15, and 20 mg/mL, respectively. The corresponding values for the flattened beam were 0.09±0.14, 0.50±0.28, 1.19±0.29, 1.68±0.39, and 2.34±0.24. For Ir-192 with direct contact, the enhanced were 0.50±0.14, 2.79±0.17, 5.49±0.12, 8.19±0.14, and 10.80±0.13. Gd-containing materials used in MRI as contrast agents can also potentially serve as radiosensitizers in radiotherapy. This study demonstrates that Gd can be used to enhance radiation dose in target volumes not only in HDR brachytherapy, but also in 6 MV FFF external beam radiotherapy, but higher than the currently used clinical concentration (>5 mg/mL) would be needed.

Monte Carlo comparison of superficial dose between flattening filter free and flattened beams.

This study investigates the superficial dose from FFF beams in comparison with the conventional flattened ones using a Monte Carlo (MC) method. Published phase-space files which incorporated real geometry of a TrueBeam accelerator were used for the dose calculation in phantom and clinical cases. The photon fluence on the central axis is 3 times that of a flattened beam for a 6 MV FFF beam and 5 times for a 10 MV beam. The mean energy across the field in air at the phantom surface is 0.92-0.95 MeV for the 6 MV FFF beam and 1.18-1.30 MeV for the corresponding flattened beam. At 10 MV, the values are 1.52-1.72 and 2.15-2.87 MeV for the FFF and flattened beams, respectively. The phantom dose at the depth of 1 mm in the 6 MV FFF beam is 6% ± 2.5% (of the maximum dose) higher compared to the flattened beam for a 25 × 25 cm(2) field and 14.6% ± 1.9% for the 2 × 2 cm(2) field. For the 10 MV beam, the corresponding differences are 3.4% ± 1.5% and 10.7% ± 0.6%. The skin dose difference at selected points on the patient's surface between the plans using FFF and flattened beams in the head-and-neck case was 6.5% ± 2.3% (1SD), and for the breast case it was 6.4% ± 2.3%. The Monte Carlo simulations showed that due to the lower mean energy in the FFF beam, the clinical superficial dose is higher without the flattening filter compared to the flattened beam.