Technical Note: Proof of concept for radiomics-based quality assurance for computed tomography.

PURPOSE: Routine quality assurance (QA) testing to identify malfunctions in medical imaging devices is a standard practice and plays an important role in meeting quality standards. However, current daily computed tomography (CT) QA techniques have proven to be inadequate for the detection of subtle artifacts on scans. Therefore, we investigated the ability of a radiomics phantom to detect subtle artifacts not detected in conventional daily QA. METHODS: An updated credence cartridge radiomics phantom was used in this study, with a focus on two of the cartridges (rubber and cork) in the phantom. The phantom was scanned using a Siemens Definition Flash CT scanner, which was reported to produce a subtle line pattern artifact. Images were then imported into the IBEX software program, and 49 features were extracted from the two cartridges using four different preprocessing techniques. Each feature was then compared with features for the same scanner several months previously and with features from controlled CT scans obtained using 100 scanners. RESULTS: Of 196 total features for the test scanner, 79 (40%) from the rubber cartridge and 70 (36%) from the cork cartridge were three or more standard deviations away from the mean of the controlled scan population data. Feature values for the artifact-producing scanner were closer to the population mean when features were preprocessed with Butterworth smoothing. The feature most sensitive to the artifact was co-occurrence matrix maximum probability. The deviation from the mean for this feature was more than seven times greater when the scanner was malfunctioning (7.56 versus 1.01). CONCLUSIONS: Radiomics features extracted from a texture phantom were able to identify an artifact-producing scanner as an outlier among 100 CT scanners. This preliminary analysis demonstrated the potential of radiomics in CT QA to identify subtle artifacts not detected using the currently employed daily QA techniques.

Correction to: A predictive model for distinguishing radiation necrosis from tumour progression after gamma knife radiosurgery based on radiomic features from MR images.

The original version of this article, published on 24 November 2017, unfortunately contained a mistake.

A predictive model for distinguishing radiation necrosis from tumour progression after gamma knife radiosurgery based on radiomic features from MR images.

OBJECTIVES: To develop a model using radiomic features extracted from MR images to distinguish radiation necrosis from tumour progression in brain metastases after Gamma Knife radiosurgery. METHODS: We retrospectively identified 87 patients with pathologically confirmed necrosis (24 lesions) or progression (73 lesions) and calculated 285 radiomic features from four MR sequences (T1, T1 post-contrast, T2, and fluid-attenuated inversion recovery) obtained at two follow-up time points per lesion per patient. Reproducibility of each feature between the two time points was calculated within each group to identify a subset of features with distinct reproducible values between two groups. Changes in radiomic features from one time point to the next (delta radiomics) were used to build a model to classify necrosis and progression lesions. RESULTS: A combination of five radiomic features from both T1 post-contrast and T2 MR images were found to be useful in distinguishing necrosis from progression lesions. Delta radiomic features with a RUSBoost ensemble classifier had an overall predictive accuracy of 73.2% and an area under the curve value of 0.73 in leave-one-out cross-validation. CONCLUSIONS: Delta radiomic features extracted from MR images have potential for distinguishing radiation necrosis from tumour progression after radiosurgery for brain metastases. KEY POINTS: • Some radiomic features showed better reproducibility for progressive lesions than necrotic ones • Delta radiomic features can help to distinguish radiation necrosis from tumour progression • Delta radiomic features had better predictive value than did traditional radiomic features.

Cost-effective immobilization for whole brain radiation therapy.

To investigate the inter- and intra-fraction motion associated with the use of a low-cost tape immobilization technique as an alternative to thermoplastic immobilization masks for whole-brain treatments. The results of this study may be of interest to clinical staff with severely limited resources (e.g., in low-income countries) and also when treating patients who cannot tolerate standard immobilization masks. Setup reproducibility of eight healthy volunteers was assessed for two different immobilization techniques. (a) One strip of tape was placed across the volunteer's forehead and attached to the sides of the treatment table. (b) A second strip was added to the first, under the chin, and secured to the table above the volunteer's head. After initial positioning, anterior and lateral photographs were acquired. Volunteers were positioned five times with each technique to allow calculation of inter-fraction reproducibility measurements. To estimate intra-fraction reproducibility, 5-minute anterior and lateral videos were taken for each technique per volunteer. An in-house software was used to analyze the photos and videos to assess setup reproducibility. The maximum intra-fraction displacement for all volunteers was 2.8 mm. Intra-fraction motion increased with time on table. The maximum inter-fraction range of positions for all volunteers was 5.4 mm. The magnitude of inter-fraction and intra-fraction motion found using the "1-strip" and "2-strip" tape immobilization techniques was comparable to motion restrictions provided by a thermoplastic mask for whole-brain radiotherapy. The results suggest that tape-based immobilization techniques represent an economical and useful alternative to the thermoplastic mask.

Comparison of reconstructed prompt gamma emissions using maximum likelihood estimation and origin ensemble algorithms for a Compton camera system tailored to proton range monitoring.

Compton-based prompt gamma (PG) imaging is being investigated by several groups as a potential solution for in vivo range monitoring in proton therapy. The performance of this technique depends on the detector system as well as the ability of the reconstruction method to obtain good spatial resolution to establish a quantitative correlation between the PG emission and the proton beam range in the patient. To evaluate the feasibility of PG imaging for range monitoring, we quantitatively evaluated the emission distributions reconstructed by a Maximum Likelihood Expectation Maximization (MLEM) and a Stochastic Origin Ensemble (SOE) algorithm. To this end, we exploit experimental and Monte Carlo (MC) simulation data acquired with the Polaris-J Compton Camera (CC) prototype. The differences between the proton beam range (RD) defined as the 80% distal dose fall-off and the PG range (RPG), obtained by fitting the distal end of the reconstructed profile with a sigmoid function, were quantified. A comparable performance of both reconstruction algorithms was found. For both experimental and simulated irradiation scenarios, the correlation between RD and RPG was within 5 mm. These values were consistent with the ground truth distance (RD-RPGg≈ 3 mm) calculated by using the expected PG emission available from MC simulation. Furthermore, shifts of 3 mm in the proton beam range were resolved with the MLEM algorithm by calculating the relative difference between the RPG for each reconstructed profile. In non-homogeneous targets, the spatial changes in the PG emission due to the different materials could not be fully resolved from the reconstructed profiles; however, the fall-off region still resembled the ground truth emission. For this scenario, the PG correlation (RD-RPG) varied from 0.1 mm to 4 mm, which is close to the ground truth correlation (3 mm). This work provides a framework for the evaluation of the range monitoring capabilities of a CC device for PG imaging. The two investigated image reconstruction algorithms showed a comparable and consistent performance for homogeneous and heterogeneous targets.

The effects of Compton camera data acquisition and readout timing on PG imaging for proton range verification.

The purpose of this study was to determine how the characteristics of the data acquisition (DAQ) electronics of a Compton camera (CC) affect the quality of the recorded prompt gamma (PG) interaction data and the reconstructed images, during clinical proton beam delivery. We used the Monte-Carlo-plus-Detector-Effect (MCDE) model to simulate the delivery of a 150 MeV clinical proton pencil beam to a tissue-equivalent plastic phantom. With the MCDE model we analyzed how the recorded PG interaction data changed as two characteristics of the DAQ electronics of a CC were changed: (1) the number of data readout channels; and (2) the active charge collection, readout, and reset time. As the proton beam dose rate increased, the number of recorded PG single-, double-, and triple-scatter events decreased by a factor of 60× for the current DAQ configuration of the CC. However, as the DAQ readout channels were increased and the readout/reset timing decreased, the number of recorded events decreased by <5× at the highest clinical dose rate. The increased number of readout channels and reduced readout/reset timing also resulted in higher quality recorded data. That is, a higher percentage of the recorded double- and triple-scatters were "true" events (caused by a single incident gamma) and not "false" events (caused by multiple incident gammas). The increase in the number and the quality of recorded data allowed higher quality PG images to be reconstructed even at the highest clinical dose rates.

Applications of Machine Learning to Improve the Clinical Viability of Compton Camera Based in vivo Range Verification in Proton Radiotherapy.

We studied the application of a deep, fully connected Neural Network (NN) to process prompt gamma (PG) data measured by a Compton camera (CC) during the delivery of clinical proton radiotherapy beams. The network identifies 1) recorded "bad" PG events arising from background noise during the measurement, and 2) the correct ordering of PG interactions in the CC to help improve the fidelity of "good" data used for image reconstruction. PG emission from a tissue-equivalent target during irradiation with a 150 MeV proton beam delivered at clinical dose rates was measured with a prototype CC. Images were reconstructed from both the raw measured data and the measured data that was further processed with a neural network (NN) trained to identify "good" and "bad" PG events and predict the ordering of individual interactions within the good PG events. We determine if NN processing of the CC data could improve the reconstructed PG images to a level in which they could provide clinically useful information about the in vivo range and range shifts of the proton beams delivered at full clinical dose rates. Results showed that a deep, fully connected NN improved the achievable contrast to noise ratio (CNR) in our images by more than a factor of 8x. This allowed the path, range, and lateral width of the clinical proton beam within a tissue equivalent target to easily be identified from the PG images, even at the highest dose rates of a 150 MeV proton beam used for clinical treatments. On average, shifts in the beam range as small as 3 mm could be identified. However, when limited by the amount of PG data measured with our prototype CC during the delivery of a single proton pencil beam (~1 × 109 protons), the uncertainty in the reconstructed PG images limited the identification of range shift to ~5 mm. Substantial improvements in CC images were obtained during clinical beam delivery through NN pre-processing of the measured PG data. We believe this shows the potential of NNs to help improve and push CC-based PG imaging toward eventual clinical application for proton RT treatment delivery verification.

Secondary Particle Interactions in a Compton Camera Designed for in vivo Range Verification of Proton Therapy.

The purpose of this study was to determine the types, proportions, and energies of secondary particle interactions in a Compton camera (CC) during the delivery of clinical proton beams. The delivery of clinical proton pencil beams ranging from 70 to 200 MeV incident on a water phantom was simulated using Geant4 software (version 10.4). The simulation included a CC similar to the configuration of a Polaris J3 CC designed to image prompt gammas (PGs) emitted during proton beam irradiation for the purpose of in vivo range verification. The interaction positions and energies of secondary particles in each CC detector module were scored. For a 150-MeV proton beam, a total of 156,688(575) secondary particles per 108 protons, primarily composed of gamma rays (46.31%), neutrons (41.37%), and electrons (8.88%), were found to reach the camera modules, and 79.37% of these particles interacted with the modules. Strategies for using CCs for proton range verification should include methods of reducing the large neutron backgrounds and low-energy non-PG radiation. The proportions of interaction types by module from this study may provide information useful for background suppression.

Radiomics feature robustness as measured using an MRI phantom.

Radiomics involves high-throughput extraction of large numbers of quantitative features from medical images and analysis of these features to predict patients' outcome and support clinical decision-making. However, radiomics features are sensitive to several factors, including scanning protocols. The purpose of this study was to investigate the robustness of magnetic resonance imaging (MRI) radiomics features with various MRI scanning protocol parameters and scanners using an MRI radiomics phantom. The variability of the radiomics features with different scanning parameters and repeatability measured using a test-retest scheme were evaluated using the coefficient of variation and intraclass correlation coefficient (ICC) for both T1- and T2-weighted images. For variability measures, the features were categorized into three groups: large, intermediate, and small variation. For repeatability measures, the average T1- and T2-weighted image ICCs for the phantom (0.963 and 0.959, respectively) were higher than those for a healthy volunteer (0.856 and 0.849, respectively). Our results demonstrated that various radiomics features are dependent on different scanning parameters and scanners. The radiomics features with a low coefficient of variation and high ICC for both the phantom and volunteer can be considered good candidates for MRI radiomics studies. The results of this study will assist current and future MRI radiomics studies.

Computed Tomography Radiomics Kinetics as Early Imaging Correlates of Osteoradionecrosis in Oropharyngeal Cancer Patients.

Osteoradionecrosis (ORN) is a major side-effect of radiation therapy in oropharyngeal cancer (OPC) patients. In this study, we demonstrate that early prediction of ORN is possible by analyzing the temporal evolution of mandibular subvolumes receiving radiation. For our analysis, we use computed tomography (CT) scans from 21 OPC patients treated with Intensity Modulated Radiation Therapy (IMRT) with subsequent radiographically-proven ≥ grade II ORN, at three different time points: pre-IMRT, 2-months, and 6-months post-IMRT. For each patient, radiomic features were extracted from a mandibular subvolume that developed ORN and a control subvolume that received the same dose but did not develop ORN. We used a Multivariate Functional Principal Component Analysis (MFPCA) approach to characterize the temporal trajectories of these features. The proposed MFPCA model performs the best at classifying ORN vs. Control subvolumes with an area under curve (AUC) = 0.74 [95% confidence interval (C.I.): 0.61-0.90], significantly outperforming existing approaches such as a pre-IMRT features model or a delta model based on changes at intermediate time points, i.e., at 2- and 6-month follow-up. This suggests that temporal trajectories of radiomics features derived from sequential pre- and post-RT CT scans can provide markers that are correlates of RT-induced mandibular injury, and consequently aid in earlier management of ORN.

Computational model for detector timing effects in Compton-camera based prompt-gamma imaging for proton radiotherapy.

This paper describes a realistic simulation of a Compton-camera (CC) based prompt-gamma (PG) imaging system for proton range verification for a range of clinical dose rates, and its comparison to PG measured data with a pre-clinical CC. We used a Monte Carlo plus Detector Effects (MCDE) model to simulate the production of prompt gamma-rays (PG) and their energy depositions in the CC. With Monte Carlo, we simulated PG emission resulting from irradiation of a high density polyethylene phantom with a 150 MeV proton pencil beam at dose rates of 5.0 × 108, 2.6 × 109, and 4.6 × 109 p+ s-1. Realistic detector timing effects (e.g. delayed triggering time, event-coincidence, dead time, etc,) were added in post-processing to allow for flexible count rate variations. We acquired PG emission measurements with our pre-clinical CC during irradiation with a clinical 150 MeV proton pencil beam at the same dose rates. For simulations and measurements, three primary changes could be seen in the PG emission data as the dose rate increased: (1) reduction in the total number of detected events due to increased dead-time percentage; (2) increase in false-coincidence events (i.e. multiple PGs interacting, rather than a single PG scatter); and (3) loss of distinct PG emission peaks in the energy spectrum. We used the MCDE model to estimate the quality of our measured PG data, primarily with regards to true and false double-scatters and triple-scatters recorded by the CC. The simulation results showed that of the recorded double-scatter PG interactions 22%, 57%, and 70% were false double-scatters and for triple-scatter interactions 3%, 21%, and 35% were false events at 5.0 × 108, 2.6 × 109, and 4.6 × 109 p+ s-1, respectively. These false scatter events represent noise in the data, and the high percentage of these events in the data represents a major limitation in our ability to produce usable PG images with our prototype CC.

Use of uniform shots for robust planning of mask-based treatment in Gamma Knife Icon.

PURPOSE: To verify whether Icon automatic correction is robust in preserving plan quality. MATERIALS/METHODS: An end-to-end phantom was used to verify Icon's correction accuracy qualitatively. For quantitative assessment, two plans, a composite- and a uniform-shot-only, were created for an elliptical- (E) and a sausage-shaped (S) lesion inside a PseudoPatient head phantom with a film insert. The phantom was irradiated in the planned and three other positions under each plan: 14° pitch (B); 14° rotation + 8° pitch (C); 95° rotation + 4-cm shift (D). RESULTS: Icon accurately corrects the locations of the shots. For the uniform-shot plans: all gamma index passing rates were >97%, and the differences between the planned and the delivery doses (minimum, maximum, and mean) were all ≤0.1 Gy. For the composite-shot plans, however, the dose differences increased as the phantom was shifted through positions B-D, with a gamma index passing rate of 61% for lesion-E in position D, and 92%, 79%, and 45% for lesion-S in positions B, C, and D, respectively. CONCLUSIONS: Plans using only uniform shots are more robust to deviations in treatment position. The tolerance for such deviations may be lower for plans using composite shots.

Matching and Homogenizing Convolution Kernels for Quantitative Studies in Computed Tomography.

The sharpness of the kernels used for image reconstruction in computed tomography affects the values of the quantitative image features. We sought to identify the kernels that produce similar feature values to enable a more effective comparison of images produced using scanners from different manufactures. We also investigated a new image filter designed to change the kernel-related component of the frequency spectrum of a postreconstruction image from that of the initial kernel to that of a preferred kernel. A radiomics texture phantom was imaged using scanners from GE, Philips, Siemens, and Toshiba. Images were reconstructed multiple times, varying the kernel from smooth to sharp. The phantom comprised 10 cartridges of various textures. A semiautomated method was used to produce 8 × 2 × 2 cm regions of interest for each cartridge and for all scans. For each region of interest, 38 radiomics features from the categories intensity direct (n = 12), gray-level co-occurrence matrix (n = 21), and neighborhood gray-tone difference matrix (n = 5) were extracted. We then calculated the fractional differences of the features from those of the baseline kernel (GE Standard). To gauge the importance of the differences, we scaled them by the coefficient of variation of the same feature from a cohort of patients with non-small cell lung cancer. The noise power spectra for each kernel were estimated from the phantom's solid acrylic cartridge, and kernel-homogenization filters were developed from these estimates. The Philips C, Siemens B30f, and Toshiba FC24 kernels produced feature values most similar to GE Standard. The kernel homogenization filters reduced the median differences from baseline to less than 1 coefficient of variation in the patient population for all of the GE, Philips, and Siemens kernels except for GE Edge and Toshiba kernels. For prospective computed tomographic radiomics studies, the scanning protocol should specify kernels that have been shown to produce similar feature values. For retrospective studies, kernel homogenization filters can be designed and applied to reduce the kernel-related differences in the feature values.

Imaging-Genomic Study of Head and Neck Squamous Cell Carcinoma: Associations Between Radiomic Phenotypes and Genomic Mechanisms via Integration of The Cancer Genome Atlas and The Cancer Imaging Archive.

PURPOSE: Recent data suggest that imaging radiomic features of a tumor could be indicative of important genomic biomarkers. Understanding the relationship between radiomic and genomic features is important for basic cancer research and future patient care. We performed a comprehensive study to discover the imaginggenomic associations in head and neck squamous cell carcinoma (HNSCC) and explore the potential of predicting tumor genomic alternations using radiomic features. METHODS: Our retrospective study integrated whole-genome multiomics data from The Cancer Genome Atlas with matched computed tomography imaging data from The Cancer Imaging Archive for the same set of 126 patients with HNSCC. Linear regression and gene set enrichment analysis were used to identify statistically significant associations between radiomic imaging and genomic features. Random forest classifier was used to predict the status of two key HNSCC molecular biomarkers, human papillomavirus and disruptive TP53 mutation, on the basis of radiomic features. RESULTS: Widespread and statistically significant associations were discovered between genomic features (including microRNA expression, somatic mutations, and transcriptional activity, copy number variations, and promoter region DNA methylation changes of pathways) and radiomic features characterizing the size, shape, and texture of tumor. Prediction of human papillomavirus and TP53 mutation status using radiomic features achieved areas under the receiver operating characteristic curve of 0.71 and 0.641, respectively. CONCLUSION: Our exploratory study suggests that radiomic features are associated with genomic characteristics at multiple molecular layers in HNSCC and provides justification for continued development of radiomics as biomarkers for relevant genomic alterations in HNSCC.

Radiomics features of the primary tumor fail to improve prediction of overall survival in large cohorts of CT- and PET-imaged head and neck cancer patients.

Radiomics studies require many patients in order to power them, thus patients are often combined from different institutions and using different imaging protocols. Various studies have shown that imaging protocols affect radiomics feature values. We examined whether using data from cohorts with controlled imaging protocols improved patient outcome models. We retrospectively reviewed 726 CT and 686 PET images from head and neck cancer patients, who were divided into training or independent testing cohorts. For each patient, radiomics features with different preprocessing were calculated and two clinical variables-HPV status and tumor volume-were also included. A Cox proportional hazards model was built on the training data by using bootstrapped Lasso regression to predict overall survival. The effect of controlled imaging protocols on model performance was evaluated by subsetting the original training and independent testing cohorts to include only patients whose images were obtained using the same imaging protocol and vendor. Tumor volume, HPV status, and two radiomics covariates were selected for the CT model, resulting in an AUC of 0.72. However, volume alone produced a higher AUC, whereas adding radiomics features reduced the AUC. HPV status and one radiomics feature were selected as covariates for the PET model, resulting in an AUC of 0.59, but neither covariate was significantly associated with survival. Limiting the training and independent testing to patients with the same imaging protocol reduced the AUC for CT patients to 0.55, and no covariates were selected for PET patients. Radiomics features were not consistently associated with survival in CT or PET images of head and neck patients, even within patients with the same imaging protocol.

Correction: Harmonizing the pixel size in retrospective computed tomography radiomics studies.

[This corrects the article DOI: 10.1371/journal.pone.0178524.].

Guidelines and Experience Using Imaging Biomarker Explorer (IBEX) for Radiomics.

Imaging Biomarker Explorer (IBEX) is an open-source tool for medical imaging radiomics work. The purpose of this paper is to describe how to use IBEX's graphical user interface (GUI) and to demonstrate how IBEX calculated features have been used in clinical studies. IBEX allows for the import of DICOM images with DICOM radiation therapy structure files or Pinnacle files. Once the images are imported, IBEX has tools within the Data Selection GUI to manipulate the viewing of the images, measure voxel values and distances, and create and edit contours. IBEX comes with 27 preprocessing and 132 feature choices to design feature sets. Each preprocessing and feature category has parameters that can be altered. The output from IBEX is a spreadsheet that contains: 1) each feature from the feature set calculated for each contour in a data set, 2) image information about each contour in a data set, and 3) a summary of the preprocessing and features used with their selected parameters. Features calculated from IBEX have been used in studies to test the variability of features under different imaging conditions and in survival models to improve current clinical models.

Prognostic value of combining a quantitative image feature from positron emission tomography with clinical factors in oligometastatic non-small cell lung cancer.

BACKGROUND AND PURPOSE: Oligometastatic non-small cell lung cancer (NSCLC) is a heterogeneous condition with few known risk stratification factors. A quantitative imaging feature (QIF) on positron emission tomography (PET), gray-level co-occurrence matrix energy, has been linked with outcome of nonmetastatic NSCLC. We hypothesized that GLCM energy would enhance the ability of models comprising standard clinical prognostic factors (CPFs) to stratify oligometastatic patients based on overall survival (OS). MATERIALS AND METHODS: We assessed 79 patients with oligometastatic NSCLC (≤3 metastases) diagnosed in 2007-2015. The primary and largest metastases at diagnosis were delineated on pretreatment scans with GLCM energy extracted using imaging biomarker explorer (IBEX) software. Iterative stepwise elimination feature selection based on the Akaike information criterion identified the optimal model comprising CPFs for predicting OS in a multivariate Cox proportional hazards model. GLCM energy was tested for improving prediction accuracy. RESULTS: Energy was a significant predictor of OS (P = 0.028) in addition to the selected CPFs. The c-indexes for the CPF-only and CPF + Energy models were 0.720 and 0.739. CONCLUSIONS: Incorporating Energy strengthened a CPF model for predicting OS. These findings support further exploration of QIFs, including markers of the primary tumor vs. those of the metastatic sites.

Lung tumor segmentation methods: Impact on the uncertainty of radiomics features for non-small cell lung cancer.

PURPOSE: To evaluate the uncertainty of radiomics features from contrast-enhanced breath-hold helical CT scans of non-small cell lung cancer for both manual and semi-automatic segmentation due to intra-observer, inter-observer, and inter-software reliability. METHODS: Three radiation oncologists manually delineated lung tumors twice from 10 CT scans using two software tools (3D-Slicer and MIM Maestro). Additionally, three observers without formal clinical training were instructed to use two semi-automatic segmentation tools, Lesion Sizing Toolkit (LSTK) and GrowCut, to delineate the same tumor volumes. The accuracy of the semi-automatic contours was assessed by comparison with physician manual contours using Dice similarity coefficients and Hausdorff distances. Eighty-three radiomics features were calculated for each delineated tumor contour. Informative features were identified based on their dynamic range and correlation to other features. Feature reliability was then evaluated using intra-class correlation coefficients (ICC). Feature range was used to evaluate the uncertainty of the segmentation methods. RESULTS: From the initial set of 83 features, 40 radiomics features were found to be informative, and these 40 features were used in the subsequent analyses. For both intra-observer and inter-observer reliability, LSTK had higher reliability than GrowCut and the two manual segmentation tools. All observers achieved consistently high ICC values when using LSTK, but the ICC value varied greatly for each observer when using GrowCut and the manual segmentation tools. For inter-software reliability, features were not reproducible across the software tools for either manual or semi-automatic segmentation methods. Additionally, no feature category was found to be more reproducible than another feature category. Feature ranges of LSTK contours were smaller than those of manual contours for all features. CONCLUSION: Radiomics features extracted from LSTK contours were highly reliable across and among observers. With semi-automatic segmentation tools, observers without formal clinical training were comparable to physicians in evaluating tumor segmentation.

Effect of tube current on computed tomography radiomic features.

Variability in the x-ray tube current used in computed tomography may affect quantitative features extracted from the images. To investigate these effects, we scanned the Credence Cartridge Radiomics phantom 12 times, varying the tube current from 25 to 300 mA∙s while keeping the other acquisition parameters constant. For each of the scans, we extracted 48 radiomic features from the categories of intensity histogram (n = 10), gray-level run length matrix (n = 11), gray-level co-occurrence matrix (n = 22), and neighborhood gray tone difference matrix (n = 5). To gauge the size of the tube current effects, we scaled the features by the coefficient of variation of the corresponding features extracted from images of non-small cell lung cancer tumors. Variations in the tube current had more effect on features extracted from homogeneous materials (acrylic, sycamore wood) than from materials with more tissue-like textures (cork, rubber particles). Thirty-eight of the 48 features extracted from acrylic were affected by current reductions compared with only 2 of the 48 features extracted from rubber particles. These results indicate that variable x-ray tube current is unlikely to have a large effect on radiomic features extracted from computed tomography images of textured objects such as tumors.