Impact of Preprocessing Parameters in Medical Imaging-Based Radiomic Studies: A Systematic Review

BACKGROUND

Lately, radiomic studies featuring the development of a signature to use in prediction models in diagnosis or prognosis outcomes have been increasingly published. While the results are shown to be promising, these studies still have many pitfalls and limitations. One of the main issues of these studies is that radiomic features depend on how the images are preprocessed before their computation. Since, in widely known and used software for radiomic features calculation, it is possible to set these preprocessing parameters before the calculation of the radiomic feature, there are ongoing studies assessing the stability and repeatability of radiomic features to find the most suitable preprocessing parameters for every used imaging modality.

MATERIALS AND METHODS

We performed a comprehensive literature search using four electronic databases: PubMed, Cochrane Library, Embase, and Scopus. Mesh terms and free text were modeled in search strategies for databases. The inclusion criteria were studies where preprocessing parameters' influence on feature values and model predictions was addressed. Records lacking information on image acquisition parameters were excluded, and any eligible studies with full-text versions were included in the review process, while conference proceedings and monographs were disregarded. We used the QUADAS-2 (Quality Assessment of Diagnostic Accuracy Studies 2) tool to investigate the risk of bias. We synthesized our data in a table divided by the imaging modalities subgroups.

RESULTS

After applying the inclusion and exclusion criteria, we selected 43 works. This review examines the impact of preprocessing parameters on the reproducibility and reliability of radiomic features extracted from multimodality imaging (CT, MRI, CBCT, and PET/CT). Standardized preprocessing is crucial for consistent radiomic feature extraction. Key preprocessing steps include voxel resampling, normalization, and discretization, which influence feature robustness and reproducibility. In total, 44% of the included works studied the effects of an isotropic voxel resampling, and most studies opted to employ a discretization strategy. From 2021, several studies started selecting the best set of preprocessing parameters based on models' best performance. As for comparison metrics, ICC was the most used in MRI studies in 58% of the screened works.

CONCLUSIONS

From our work, we highlighted the need to harmonize the use of preprocessing parameters and their values, especially in light of future studies of prospective studies, which are still lacking in the current literature.

Radiomic feature robustness in CT scans affected by fiducial marker induced streak artifacts for patients with hepatocellular carcinoma

BACKGROUND

Stereotactic body radiation therapy for hepatocellular carcinoma necessitates the implantation of gold fiducial markers in the liver, resulting in artifacts on computed tomography (CT) scans, which affect radiomic feature values.

PURPOSE

This report aims to assess the effect of these artifacts on radiomic features and how removing CT slices affects radiomic features extracted from 3D regions of interest (ROI).

METHODS

First, the range variation in 38 tumor contours unaffected by artifacts was assessed after sequentially and randomly removing 25%, 50%, 75% of slices. Subsequently, the agreement of feature values before and after removing ROI slices containing artifacts from 186 patients' CT scans was assessed with Lin's concordance correlation coefficient (CCC) and a Wilcoxon signed-rank test.

RESULTS

In artifact-free tumor volumes, at least 71% of features remain robust with up to 50% of slices removed, while 56% remains robust with up to 75% of slices removed. When comparing contours before and after removing slices containing artifacts, around a third of features in the tumor contour and surrounding area remain robust (CCC > 0.9), compared to 44% in the healthy liver. Concerning the tumor, 13% (Gray Level Size Zone Matrix) to 61% (first order) of the features remain robust (CCC > 0.9). Over 90% of features differ significantly as assessed by Wilcoxon signed-rank test, however.

CONCLUSIONS

This study demonstrates that removing slices containing artifacts is a feasible solution for the CT fiducial problem in this patient population and provides insight into which features are affected.

Automated Lung Cancer Diagnosis Applying Butterworth Filtering, Bi-Level Feature Extraction, and Sparce Convolutional Neural Network to Luna 16 CT Images

Accurate prognosis and diagnosis are crucial for selecting and planning lung cancer treatments. As a result of the rapid development of medical imaging technology, the use of computed tomography (CT) scans in pathology is becoming standard practice. An intricate interplay of requirements and obstacles characterizes computer-assisted diagnosis, which relies on the precise and effective analysis of pathology images. In recent years, pathology image analysis tasks such as tumor region identification, prognosis prediction, tumor microenvironment characterization, and metastasis detection have witnessed the considerable potential of artificial intelligence, especially deep learning techniques. In this context, an artificial intelligence (AI)-based methodology for lung cancer diagnosis is proposed in this research work. As a first processing step, filtering using the Butterworth smooth filter algorithm was applied to the input images from the LUNA 16 lung cancer dataset to remove noise without significantly degrading the image quality. Next, we performed the bi-level feature selection step using the Chaotic Crow Search Algorithm and Random Forest (CCSA-RF) approach to select features such as diameter, margin, spiculation, lobulation, subtlety, and malignancy. Next, the Feature Extraction step was performed using the Multi-space Image Reconstruction (MIR) method with Grey Level Co-occurrence Matrix (GLCM). Next, the Lung Tumor Severity Classification (LTSC) was implemented by using the Sparse Convolutional Neural Network (SCNN) approach with a Probabilistic Neural Network (PNN). The developed method can detect benign, normal, and malignant lung cancer images using the PNN algorithm, which reduces complexity and efficiently provides classification results. Performance parameters, namely accuracy, precision, F-score, sensitivity, and specificity, were determined to evaluate the effectiveness of the implemented hybrid method and compare it with other solutions already present in the literature.

A comparative analysis of preclinical computed tomography radiomics using cone-beam and micro-computed tomography scanners

Background and purpose

Radiomics analysis extracts quantitative data (features) from medical images. These features could potentially reflect biological characteristics and act as imaging biomarkers within precision medicine. However, there is a lack of cross-comparison and validation of radiomics outputs which is paramount for clinical implementation. In this study, we compared radiomics outputs across two computed tomography (CT)-based preclinical scanners.

Materials and methods

Cone beam CT (CBCT) and µCT scans were acquired using different preclinical CT imaging platforms. The reproducibility of radiomics features on each scanner was assessed using a phantom across imaging energies (40 & 60 kVp) and segmentation volumes (44-238 mm3). Retrospective mouse scans were used to compare feature reliability across varying tissue densities (lung, heart, bone), scanners and after voxel size harmonisation. Reliable features had an intraclass correlation coefficient (ICC) > 0.8.

Results

First order and GLCM features were the most reliable on both scanners across different volumes. There was an inverse relationship between tissue density and feature reliability, with the highest number of features in lung (CBCT=580, µCT=734) and lowest in bone (CBCT=110, µCT=560). Comparable features for lung and heart tissues increased when voxel sizes were harmonised. We have identified tissue-specific preclinical radiomics signatures in mice for the lung (133), heart (35), and bone (15).

Conclusions

Preclinical CBCT and µCT scans can be used for radiomics analysis to support the development of meaningful radiomics signatures. This study demonstrates the importance of standardisation and emphasises the need for multi-centre studies.

Increasing the efficiency of cone-beam CT based delta-radiomics using automated contours to predict radiotherapy-related toxicities in prostate cancer

Extracting longitudinal image quantitative data, known as delta-radiomics, has the potential to capture changes in a patient's anatomy throughout the course of radiation treatment for prostate cancer. Some of the major challenges of delta-radiomics studies are contouring the structures for individual fractions and accruing patients' data in an efficient manner. The manual contouring process is often time consuming and would limit the efficiency of accruing larger sample sizes for future studies. The problem is amplified because the contours are often made by highly trained radiation oncologists with limited time to dedicate to research studies of this nature. This work compares the use of automated prostate contours generated using a deformable image-based algorithm to make predictive models of genitourinary and changes in total international prostate symptom score in comparison to manually contours for a cohort of fifty patients. Area under the curve of manual and automated models were compared using the Delong test. This study demonstrated that the delta-radiomics models were similar for both automated and manual delta-radiomics models.

Investigating the role of imaging factors in the variability of CT-based texture analysis metrics

OBJECTIVE

This study assesses the robustness of first-order radiomic texture features namely interquartile range (IQR), coefficient of variation (CV) and standard deviation (SD) derived from computed tomography (CT) images by varying dose, reconstruction algorithms and slice thickness using scans of a uniform water phantom, a commercial anthropomorphic liver phantom, and a human liver in-vivo.

MATERIALS AND METHODS

Scans were acquired on a 16 cm detector GE Revolution Apex Edition CT scanner with variations across three different nominal slice thicknesses: 0.625, 1.25, and 2.5 mm, three different dose levels: CTDIvol of 13.86 mGy for the standard dose, 40% reduced dose and 60% reduced dose and two different reconstruction algorithms: a deep learning image reconstruction (DLIR-high) algorithm and a hybrid iterative reconstruction (IR) algorithm ASiR-V50% (AV50) were explored, varying one at a time. To assess the effect of non-linear modifications of images by AV50 and DLIR-high, images of the water phantom were also reconstructed using filtered back projection (FBP). Quantitative measures of IQR, CV and SD were extracted from twelve pre-selected, circular (1 cm diameter) regions of interest (ROIs) capturing different texture patterns across all scans.

RESULTS

Across all scans, imaging, and reconstruction settings, CV, IQR and SD were observed to increase with reduction in dose and slice thickness. An exception to this observation was found when using FBP reconstruction. Lower values of CV, IQR and SD were observed in DLIR-high reconstructions compared to AV50 and FBP. The Poisson statistics were more stringently noted in FBP than DLIR-high and AV50, due to the non-linear nature of the latter two algorithms.

CONCLUSION

Variation in image noise due to dose reduction algorithms, tube current, and slice thickness show a consistent trend across phantom and patient scans. Prospective evaluation across multiple centers, scanners and imaging protocols is needed for establishing quality assurance standards of radiomics.

Technical and clinical considerations of a physical liver phantom for CT radiomics analysis

OBJECTIVE

This study identifies key characteristics to help build a physical liver computed tomography (CT) phantom for radiomics harmonization; particularly, the higher-order texture metrics.

MATERIALS AND METHODS

CT scans of a radiomics phantom comprising of 18 novel 3D printed inserts with varying size, shape, and material combinations were acquired on a 64-slice CT scanner (Brilliance 64, Philips Healthcare). The images were acquired at 120 kV, 250 mAs, CTDIvol of 16.36 mGy, 2 mm slice thickness, and iterative noise-reduction reconstruction (iDose, Philips Healthcare, Andover, MA). Radiomics analysis was performed using the Cancer Imaging Phenomics Toolkit (CaPTk), following automated segmentation of 3D regions of interest (ROI) of the 18 inserts. The findings were compared to three additional ROI obtained of an anthropomorphic liver phantom, a patient liver CT scan, and a water phantom, at comparable imaging settings. Percentage difference in radiomic metrics values between phantom and tissue was used to assess the biological equivalency and <10% was used to claim equivalent.

RESULTS

The HU for all 18 ROI from the phantom ranged from -30 to 120 which is within clinically observed HU range of the liver, showing that our phantom material (T3-6B) is representative of biological CT tissue densities (liver) with >50% radiomic features having <10% difference from liver tissue. Based on the assessment of the Neighborhood Gray Tone Difference Matrix (NGTDM) metrics it is evident that the water phantom ROI show extreme values compared to the ROIs from the phantom. This result may further reinforce the difference between a structureless quantity such as water HU values and tissue HU values found in liver.

CONCLUSION

The 3-D printed patterns of the constructed radiomics phantom cover a wide span of liver tissue textures seen in CT images. Using our results, texture metrics can be selectively harmonized to establish clinically relevant and reliable radiomics panels.

Novel Functional Radiomics for Prediction of Cardiac Positron Emission Tomography Avidity in Lung Cancer Radiotherapy

PURPOSE

Traditional methods of evaluating cardiotoxicity focus on radiation doses to the heart. Functional imaging has the potential to provide improved prediction for cardiotoxicity for patients with lung cancer. Fluorine-18 (18F) fluorodeoxyglucose (FDG)-positron emission tomography (PET)/computed tomography (CT) imaging is routinely obtained in a standard cancer staging workup. This work aimed to develop a radiomics model predicting clinical cardiac assessment using 18F-FDG PET/CT scans before thoracic radiation therapy.

METHODS

Pretreatment 18F-FDG PET/CT scans from three study populations (N = 100, N = 39, N = 70) were used, comprising two single-institutional protocols and one publicly available data set. A clinician (V.J.) classified the PET/CT scans per clinical cardiac guidelines as no uptake, diffuse uptake, or focal uptake. The heart was delineated, and 210 novel functional radiomics features were selected to classify cardiac FDG uptake patterns. Training data were divided into training (80%)/validation (20%) sets. Feature reduction was performed using the Wilcoxon test, hierarchical clustering, and recursive feature elimination. Ten-fold cross-validation was carried out for training, and the accuracy of the models to predict clinical cardiac assessment was reported.

RESULTS

From 202 of 209 scans, cardiac FDG uptake was scored as no uptake (39.6%), diffuse uptake (25.3%), and focal uptake (35.1%), respectively. Sixty-two independent radiomics features were reduced to nine clinically pertinent features. The best model showed 93% predictive accuracy in the training data set and 80% and 92% predictive accuracy in two external validation data sets.

CONCLUSION

This work used an extensive patient data set to develop a functional cardiac radiomic model from standard-of-care 18F-FDG PET/CT scans, showing good predictive accuracy. The radiomics model has the potential to provide an automated method to predict existing cardiac conditions and provide an early functional biomarker to identify patients at risk of developing cardiac complications after radiotherapy.

Quantitation of Oncologic Image Features for Radiomic Analyses in PET

Radiomics is an emerging and exciting field of study involving the extraction of many quantitative features from radiographic images. Positron emission tomography (PET) images are used in cancer diagnosis and staging. Utilizing radiomics on PET images can better quantify the spatial relationships between image voxels and generate more consistent and accurate results for diagnosis, prognosis, treatment, etc. This chapter gives the general steps a researcher would take to extract PET radiomic features from medical images and properly develop models to implement.

Radiotranscriptomics of non-small cell lung carcinoma for assessing high-level clinical outcomes using a machine learning-derived multi-modal signature

BACKGROUND

Multi-omics research has the potential to holistically capture intra-tumor variability, thereby improving therapeutic decisions by incorporating the key principles of precision medicine. The purpose of this study is to identify a robust method of integrating features from different sources, such as imaging, transcriptomics, and clinical data, to predict the survival and therapy response of non-small cell lung cancer patients.

METHODS

2996 radiomics, 5268 transcriptomics, and 8 clinical features were extracted from the NSCLC Radiogenomics dataset. Radiomics and deep features were calculated based on the volume of interest in pre-treatment, routine CT examinations, and then combined with RNA-seq and clinical data. Several machine learning classifiers were used to perform survival analysis and assess the patient's response to adjuvant chemotherapy. The proposed analysis was evaluated on an unseen testing set in a k-fold cross-validation scheme. Score- and concatenation-based multi-omics were used as feature integration techniques.

RESULTS

Six radiomics (elongation, cluster shade, entropy, variance, gray-level non-uniformity, and maximal correlation coefficient), six deep features (NasNet-based activations), and three transcriptomics (OTUD3, SUCGL2, and RQCD1) were found to be significant for therapy response. The examined score-based multi-omic improved the AUC up to 0.10 on the unseen testing set (0.74 ± 0.06) and the balance between sensitivity and specificity for predicting therapy response for 106 patients, resulting in less biased models and improving upon the either highly sensitive or highly specific single-source models. Six radiomics (kurtosis, GLRLM- and GLSZM-based non-uniformity from images with no filtering, biorthogonal, and daubechies wavelets), seven deep features (ResNet-based activations), and seven transcriptomics (ELP3, ZZZ3, PGRMC2, TRAK1, ATIC, USP7, and PNPLA2) were found to be significant for the survival analysis. Accordingly, the survival analysis for 115 patients was also enhanced up to 0.20 by the proposed score-based multi-omics in terms of the C-index (0.79 ± 0.03).

CONCLUSIONS

Compared to single-source models, multi-omics integration has the potential to improve prediction performance, increase model stability, and reduce bias for both treatment response and survival analysis.

A radiomic- and dosiomic-based machine learning regression model for pretreatment planning in 177 Lu-DOTATATE therapy

BACKGROUND

Standardized patient-specific pretreatment dosimetry planning is mandatory in the modern era of nuclear molecular radiotherapy, which may eventually lead to improvements in the final therapeutic outcome. Only a comprehensive definition of a dosage therapeutic window encompassing the range of absorbed doses, that is, helpful without being detrimental can lead to therapy individualization and improved outcomes. As a result, setting absorbed dose safety limits for organs at risk (OARs) requires knowledge of the absorbed dose-effect relationship. Data sets of consistent and reliable inter-center dosimetry findings are required to characterize this relationship.

PURPOSE

We developed and standardized a new pretreatment planning model consisting of a predictive dosimetry procedure for OARs in patients with neuroendocrine tumors (NETs) treated with 177 Lu-DOTATATE (Lutathera). In the retrospective study described herein, we used machine learning (ML) regression algorithms to predict absorbed doses in OARs by exploiting a combination of radiomic and dosiomic features extracted from patients' imaging data.

METHODS

Pretreatment and posttreatment data for 20 patients with NETs treated with 177 Lu-DOTATATE were collected from two clinical centers. A total of 3412 radiomic and dosiomic features were extracted from the patients' computed tomography (CT) scans and dose maps, respectively. All dose maps were generated using Monte Carlo simulations. An ML regression model was designed based on ML algorithms for predicting the absorbed dose in every OAR (liver, left kidney, right kidney, and spleen) before and after the therapy and between each therapy session, thus predicting any possible radiotoxic effects.

RESULTS

We evaluated nine ML regression algorithms. Our predictive model achieved a mean absolute dose error (MAE, in Gy) of 0.61 for the liver, 1.58 for the spleen, 1.30 for the left kidney, and 1.35 for the right kidney between pretherapy 68 Ga-DOTATOC positron emission tomography (PET)/CT and posttherapy 177 Lu-DOTATATE single photon emission (SPECT)/CT scans. Τhe best predictive performance observed was based on the gradient boost for the liver, the left kidney and the right kidney, and on the extra tree regressor for the spleen. Evaluation of the model's performance according to its ability to predict the absorbed dose in each OAR in every possible combination of pretherapy 68 Ga-DOTATOC PET/CT and any posttherapy 177 Lu-DOTATATE treatment cycle SPECT/CT scans as well as any 177 Lu-DOTATATE SPECT/CT treatment cycle and the consequent 177 Lu-DOTATATE SPECT/CT treatment cycle revealed mean absorbed dose differences ranges from -0.55 to 0.68 Gy. Incorporating radiodosiomics features from the 68 Ga-DOTATOC PET/CT and first 177 Lu-DOTATATE SPECT/CT treatment cycle scans further improved the precision and minimized the standard deviation of the predictions in nine out of 12 instances. An average improvement of 57.34% was observed (range: 17.53%-96.12%). However, it's important to note that in three instances (i.e., Ga,C.1 → C3 in spleen and left kidney, and Ga,C.1 → C2 in right kidney) we did not observe an improvement (absolute differences of 0.17, 0.08, and 0.05 Gy, respectively). Wavelet-based features proved to have high correlated predictive value, whereas non-linear-based ML regression algorithms proved to be more capable than the linear-based of producing precise prediction in our case.

CONCLUSIONS

The combination of radiomics and dosiomics has potential utility for personalized molecular radiotherapy (PMR) response evaluation and OAR dose prediction. These radiodosiomic features can potentially provide information on any possible disease recurrence and may be highly useful in clinical decision-making, especially regarding dose escalation issues.

The role of parametric feature maps to correct different volume of interest sizes: an in vivo liver MRI study

BACKGROUND

Different volume of interest (VOI) sizes influence radiomic features. This study examined if translating images into feature maps before feature sampling could compensate for these effects in liver magnetic resonance imaging (MRI).

METHODS

T1- and T2-weighted sequences from three different scanners (two 3-T scanners, one 1.5-T scanner) of 66 patients with normal abdominal MRI were included retrospectively. Three differently sized VOIs (10, 20, and 30 mm in diameter) were drawn in the liver parenchyma (right lobe), excluding adjacent structures. Ninety-three features were extracted conventionally using PyRadiomics. All images were also converted to 93 parametric feature maps using a pretested software. Agreement between the three VOI sizes was assessed with overall concordance correlation coefficients (OCCCs), while OCCCs > 0.85 were rated reproducible. OCCCs were calculated twice: for the VOI sizes of 10, 20, and 30 mm and for those of 20 and 30 mm.

RESULTS

When extracted from original images, only 4 out of the 93 features were reproducible across all VOI sizes in T1- and T2-weighted images. When the smallest VOI was excluded, 5 features (T1-weighted) and 7 features (T2-weighted) were reproducible. Extraction from parametric maps increased the number of reproducible features to 9 (T1- and T2-weighted) across all VOIs. Excluding the 10-mm VOI, reproducibility improved to 16 (T1-weighted) and 55 features (T2-weighted). The stability of all other features also increased in feature maps.

CONCLUSIONS

Translating images into parametric maps before feature extraction improves reproducibility across different VOI sizes in normal liver MRI.

RELEVANCE STATEMENT

The size of the segmented VOI influences the feature quantity of radiomics, while software-based conversion of images into parametric feature maps before feature sampling improves reproducibility across different VOI sizes in MRI of normal liver tissue.

KEY POINTS

• Parametric feature maps can compensate for different VOI sizes. • The effect seems dependent on the VOI sizes and the MRI sequence. • Feature maps can visualize features throughout the entire image stack.

Spatially-aware clustering improves AJCC-8 risk stratification performance in oropharyngeal carcinomas

OBJECTIVE

Evaluate the effectiveness of machine learning tools that incorporate spatial information such as disease location and lymph node metastatic patterns-of-spread, for prediction of survival and toxicity in HPV+ oropharyngeal cancer (OPC).

MATERIALS & METHODS

675 HPV+ OPC patients that were treated at MD Anderson Cancer Center between 2005 and 2013 with curative intent IMRT were retrospectively collected under IRB approval. Risk stratifications incorporating patient radiometric data and lymph node metastasis patterns via an anatomically-adjacent representation with hierarchical clustering were identified. These clusterings were combined into a 3-level patient stratification and included along with other known clinical features in a Cox model for predicting survival outcomes, and logistic regression for toxicity, using independent subsets for training and validation.

RESULTS

Four groups were identified and combined into a 3-level stratification. The inclusion of patient stratifications in predictive models for 5-yr Overall survival (OS), 5-year recurrence free survival, (RFS) and Radiation-associated dysphagia (RAD) consistently improved model performance measured using the area under the curve (AUC). Test set AUC improvements over models with clinical covariates, was 9 % for predicting OS, and 18 % for predicting RFS, and 7 % for predicting RAD. For models with both clinical and AJCC covariates, AUC improvement was 7 %, 9 %, and 2 % for OS, RFS, and RAD, respectively.

CONCLUSION

Including data-driven patient stratifications considerably improve prognosis for survival and toxicity outcomes over the performance achieved by clinical staging and clinical covariates alone. These stratifications generalize well to across cohorts, and sufficient information for reproducing these clusters is included.

Bayesian feature selection for radiomics using reliability metrics

Introduction: Imaging of tumors is a standard step in diagnosing cancer and making subsequent treatment decisions. The field of radiomics aims to develop imaging based biomarkers using methods rooted in artificial intelligence applied to medical imaging. However, a challenging aspect of developing predictive models for clinical use is that many quantitative features derived from image data exhibit instability or lack of reproducibility across different imaging systems or image-processing pipelines.

Methods: To address this challenge, we propose a Bayesian sparse modeling approach for image classification based on radiomic features, where the inclusion of more reliable features is favored via a probit prior formulation.

Results: We verify through simulation studies that this approach can improve feature selection and prediction given correct prior information. Finally, we illustrate the method with an application to the classification of head and neck cancer patients by human papillomavirus status, using as our prior information a reliability metric quantifying feature stability across different imaging systems.

Radiomic feature repeatability and its impact on prognostic model generalizability: A multi-institutional study on nasopharyngeal carcinoma patients

BACKGROUND AND PURPOSE

To investigate the radiomic feature (RF) repeatability via perturbation and its impact on cross-institutional prognostic model generalizability in Nasopharyngeal Carcinoma (NPC) patients.

MATERIALS AND METHODS

286 and 183 NPC patients from two institutions were included for model training and validation. Perturbations with random translations and rotations were applied to contrast-enhanced T1-weighted (CET1-w) MR images. RFs were extracted from primary tumor volume under a wide range of image filtering and discretization settings. RF repeatability was assessed by intraclass correlation coefficient (ICC), which was used to equally separate the RFs into low- and high-repeatable groups by the median value. After feature selection, multivariate Cox regression and Kaplan-Meier analysis were independently employed to develop and analyze prognostic models. Concordance index (C-index) and P-value from log-rank test were used to assess model performance.

RESULTS

Most textural RFs from high-pass wavelet-filtered images were susceptible to image perturbations. It was more prominent when a smaller discretization bin number was used (e.g., 8, mean ICC = 0.69). Using high-repeatable RFs for model development yielded a significantly higher C-index (0.63) in the validation cohort than when only low-repeatable RFs were used (0.57, P = 0.024), suggesting higher model generalizability. Besides, significant risk stratification in the validation cohort was observed only when high-repeatable RFs were used (P < 0.001).

CONCLUSION

Repeatability of RFs from high-pass wavelet-filtered CET1-w MR images of primary NPC tumor was poor, particularly when a smaller bin number was used. Exclusive use of high-repeatable RFs is suggested to safeguard model generalizability for wide-spreading clinical utilization.

Artificial intelligence for differentiating COVID-19 from other viral pneumonias on CT: comparative analysis of different models based on quantitative and radiomic approaches

BACKGROUND

To develop a pipeline for automatic extraction of quantitative metrics and radiomic features from lung computed tomography (CT) and develop artificial intelligence (AI) models supporting differential diagnosis between coronavirus disease 2019 (COVID-19) and other viral pneumonia (non-COVID-19).

METHODS

Chest CT of 1,031 patients (811 for model building; 220 as independent validation set (IVS) with positive swab for severe acute respiratory syndrome coronavirus-2 (647 COVID-19) or other respiratory viruses (384 non-COVID-19) were segmented automatically. A Gaussian model, based on the HU histogram distribution describing well-aerated and ill portions, was optimised to calculate quantitative metrics (QM, n = 20) in both lungs (2L) and four geometrical subdivisions (GS) (upper front, lower front, upper dorsal, lower dorsal; n = 80). Radiomic features (RF) of first (RF1, n = 18) and second (RF2, n = 120) order were extracted from 2L using PyRadiomics tool. Extracted metrics were used to develop four multilayer-perceptron classifiers, built with different combinations of QM and RF: Model1 (RF1-2L); Model2 (QM-2L, QM-GS); Model3 (RF1-2L, RF2-2L); Model4 (RF1-2L, QM-2L, GS-2L, RF2-2L).

RESULTS

The classifiers showed accuracy from 0.71 to 0.80 and area under the receiving operating characteristic curve (AUC) from 0.77 to 0.87 in differentiating COVID-19 versus non-COVID-19 pneumonia. Best results were associated with Model3 (AUC 0.867 ± 0.008) and Model4 (AUC 0.870 ± 0.011. For the IVS, the AUC values were 0.834 ± 0.008 for Model3 and 0.828 ± 0.011 for Model4.

CONCLUSIONS

Four AI-based models for classifying patients as COVID-19 or non-COVID-19 viral pneumonia showed good diagnostic performances that could support clinical decisions.

External validation of a CT-based radiomics signature in oropharyngeal cancer: Assessing sources of variation

BACKGROUND AND PURPOSE

Radiomics is a high-throughput approach that allows for quantitative analysis of imaging data for prognostic applications. Medical images are used in oropharyngeal cancer (OPC) diagnosis and treatment planning and these images may contain prognostic information allowing for treatment personalization. However, the lack of validated models has been a barrier to the translation of radiomic research to the clinic. We hypothesize that a previously developed radiomics model for risk stratification in OPC can be validated in a local dataset.

MATERIALS AND METHODS

The radiomics signature predicting overall survival incorporates features derived from the primary gross tumor volume of OPC patients treated with radiation +/- chemotherapy at a single institution (n = 343). Model fit, calibration, discrimination, and utility were evaluated. The signature was compared with a clinical model using overall stage and a model incorporating both radiomics and clinical data. A model detecting dental artifacts on computed tomography images was also validated.

RESULTS

The radiomics signature had a Concordance index (C-index) of 0.66 comparable to the clinical model's C-index of 0.65. The combined model significantly outperformed (C-index of 0.69, p = 0.024) the clinical model, suggesting that radiomics provides added value. The dental artifact model demonstrated strong ability in detecting dental artifacts with an area under the curve of 0.87.

CONCLUSION

This work demonstrates model performance comparable to previous validation work and provides a framework for future independent and multi-center validation efforts. With sufficient validation, radiomic models have the potential to improve traditional systems of risk stratification, treatment personalization and patient outcomes.

Cone-beam CT delta-radiomics to predict genitourinary toxicities and international prostate symptom of prostate cancer patients: a pilot study

For prostate cancer (PCa) patients treated with definitive radiotherapy (RT), acute and late RT-related genitourinary (GU) toxicities adversely impact disease-specific quality of life. Early warning of potential RT toxicities can prompt interventions that may prevent or mitigate future adverse events. During intensity modulated RT (IMRT) of PCa, daily cone-beam computed tomography (CBCT) images are used to improve treatment accuracy through image guidance. This work investigated the performance of CBCT-based delta-radiomic features (DRF) models to predict acute and sub-acute International Prostate Symptom Scores (IPSS) and Common Terminology Criteria for Adverse Events (CTCAE) version 5 GU toxicity grades for 50 PCa patients treated with definitive RT. Delta-radiomics models were built using logistic regression, random forest for feature selection, and a 1000 iteration bootstrapping leave one analysis for cross validation. To our knowledge, no prior studies of PCa have used DRF models based on daily CBCT images. AUC of 0.83 for IPSS and greater than 0.7 for CTCAE grades were achieved as early as week 1 of treatment. DRF extracted from CBCT images showed promise for the development of models predictive of RT outcomes. Future studies will include using artificial intelligence and machine learning to expand CBCT sample sizes available for radiomics analysis.

A Framework for Harmonization of Radiomics Data for Multicenter Studies and Clinical Trials

PURPOSE

Variability in computed tomography images intrinsic to individual scanners limits the application of radiomics in clinical and research settings. The development of reproducible and generalizable radiomics-based models to assess lesions requires harmonization of data. The purpose of this study was to develop, test, and analyze the efficacy of a radiomics data harmonization model.

MATERIALS AND METHODS

Radiomic features from biopsy-proven untreated hepatic metastasis (N = 380) acquired from 167 unique patients with pancreatic, colon, and breast cancers were analyzed. Radiomic features from volume-match 551 samples of normal liver tissue and 188 hepatic cysts were included as references. A novel linear mixed effect model was used to identify effects associated with lesion size, tissue type, and scanner model. Six separate machine learning models were then used to test the effectiveness of radiomic feature harmonization using multivariate analysis.

RESULTS

Proposed model identifies and removes scanner-associated effects while preserving cancer-specific functional dependence of radiomic features on the tumor size. Data harmonization improves the performance of classification models by reducing the scanner-associated variability. For example, the multiclass logistic regression model, LogitBoost, demonstrated the improvement in sensitivity in the range from 15% to 40% for each type of liver metastasis, whereas the overall model accuracy and the kappa coefficient increased by 5% and 8% accordingly.

CONCLUSION

The model removed scanner-associated effects while preserving cancer-specific functional dependence of radiomic features.

Uniqueness of radiomic features in non-small cell lung cancer

PURPOSE

The uniqueness of radiomic features, combined with their reproducibility, determines the reliability of radiomic studies. This study is to test the hypothesis that radiomic features extracted from a defined region of interest (ROI) are unique to the underlying structure (e.g., tumor).

APPROACH

Two cohorts of non-small cell lung cancer (NSCLC) patients were retrospectively retrieved from a GE and a Siemens CT scanner. The lung nodules (ROI) were delineated manually and radiomic features were extracted using IBEX. The same ROI was then translocated randomly to four other tissue regions of the same set of images: adipose, heart, lung beyond nodule, and muscle for radiomic feature extraction. Coefficient of variation (CV) within different ROIs and concordance correlation coefficient (CCC) between lung nodule and a given tissue region were calculated to test to determine feature uniqueness. The radiomic features were considered nonunique when (1) the CV < 10% and CCC > 0.85 for over 50% of patients; and (2) the CCC > 0.85 appeared in ≥2 tissue regions beyond the defined region.

RESULTS

In total, 14 patients from GE and 18 patients from Siemens are analyzed. The results show that 12 features fall below the 10% CV threshold for over 50% of patients in the GE cohort and 29 features in the Siemens cohort. According to CCC, 18 radiomic features in GE and 16 features in Siemens are identified as nonunique, with 11 overlapping features. Combining CV and CCC, 9 of 123 calculated features (7.3%) are identified as nonunique to a defined ROI.

CONCLUSIONS

Radiomic feature uniqueness should be considered to improve the reliability of radiomics study.

Impact of image quality on radiomics applications

Radiomics features extracted from medical images have been widely reported to be useful in the patient specific outcome modeling for variety of assessment and prediction purposes. Successful application of radiomics features as imaging biomarkers, however, is dependent on the robustness of the approach to the variation in each step of the modeling workflow. Variation in the input image quality is one of the main sources that impacts the reproducibility of radiomics analysis when a model is applied to broader range of medical imaging data. The quality of medical image is generally affected by both the scanner related factors such as image acquisition/reconstruction settings and the patient related factors such as patient motion. This article aimed to review the published literatures in this field that reported the impact of various imaging factors on the radiomics features through the change in image quality. The literatures were categorized by different imaging modalities and also tabulated based on the imaging parameters and the class of radiomics features included in the study. Strategies for image quality standardization were discussed based on the relevant literatures and recommendations for reducing the impact of image quality variation on the radiomics in multi-institutional clinical trial were summarized at the end of this article.

Image resampling and discretization effect on the estimate of myocardial radiomic features from T1 and T2 mapping in hypertrophic cardiomyopathy

Radiomics is emerging as a promising and useful tool in cardiac magnetic resonance (CMR) imaging applications. Accordingly, the purpose of this study was to investigate, for the first time, the effect of image resampling/discretization and filtering on radiomic features estimation from quantitative CMR T1 and T2 mapping. Specifically, T1 and T2 maps of 26 patients with hypertrophic cardiomyopathy (HCM) were used to estimate 98 radiomic features for 7 different resampling voxel sizes (at fixed bin width), 9 different bin widths (at fixed resampling voxel size), and 7 different spatial filters (at fixed resampling voxel size/bin width). While we found a remarkable dependence of myocardial radiomic features from T1 and T2 mapping on image filters, many radiomic features showed a limited sensitivity to resampling voxel size/bin width, in terms of intraclass correlation coefficient (> 0.75) and coefficient of variation (< 30%). The estimate of most textural radiomic features showed a linear significant (p < 0.05) correlation with resampling voxel size/bin width. Overall, radiomic features from T2 maps have proven to be less sensitive to image preprocessing than those from T1 maps, especially when varying bin width. Our results might corroborate the potential of radiomics from T1/T2 mapping in HCM and hopefully in other myocardial diseases.

Building reliable radiomic models using image perturbation

Radiomic model reliability is a central premise for its clinical translation. Presently, it is assessed using test–retest or external data, which, unfortunately, is often scarce in reality. Therefore, we aimed to develop a novel image perturbation-based method (IPBM) for the first of its kind toward building a reliable radiomic model. We first developed a radiomic prognostic model for head-and-neck cancer patients on a training (70%) and evaluated on a testing (30%) cohort using C-index. Subsequently, we applied the IPBM to CT images of both cohorts (Perturbed-Train and Perturbed-Test cohort) to generate 60 additional samples for both cohorts. Model reliability was assessed using intra-class correlation coefficient (ICC) to quantify consistency of the C-index among the 60 samples in the Perturbed-Train and Perturbed-Test cohorts. Besides, we re-trained the radiomic model using reliable RFs exclusively (ICC > 0.75) to validate the IPBM. Results showed moderate model reliability in Perturbed-Train (ICC: 0.565, 95%CI 0.518–0.615) and Perturbed-Test (ICC: 0.596, 95%CI 0.527–0.670) cohorts. An enhanced reliability of the re-trained model was observed in Perturbed-Train (ICC: 0.782, 95%CI 0.759–0.815) and Perturbed-Test (ICC: 0.825, 95%CI 0.782–0.867) cohorts, indicating validity of the IPBM. To conclude, we demonstrated capability of the IPBM toward building reliable radiomic models, providing community with a novel model reliability assessment strategy prior to prospective evaluation.

Discrimination of Tumor Texture Based on MRI Radiomic Features: Is There a Volume Threshold? A Phantom Study

Radiomics is emerging as a promising tool to extract quantitative biomarkers—called radiomic features—from medical images, potentially contributing to the improvement in diagnosis and treatment of oncological patients. However, technical limitations might impair the reliability of radiomic features and their ability to quantify clinically relevant tissue properties. Among these, sampling the image signal in a too-small region can reduce the ability to discriminate tissues with different properties. However, a volume threshold guaranteeing a reliable analysis, which might vary according to the imaging modality and clinical scenario, has not been assessed yet. In this study, an MRI phantom specifically developed for radiomic investigation of gynecological malignancies was used to explore how the ability of radiomic features to discriminate different image textures varies with the volume of the analyzed region. The phantom, embedding inserts with different textures, was scanned on two 1.5T and one 3T scanners, each using the T2-weighted sequence of the clinical protocol implemented for gynecological studies. Within each of the three inserts, six cylindrical regions were drawn with volumes ranging from 0.8 cm3 to 29.8 cm3, and 944 radiomic features were extracted from both original images and from images processed with different filters. For each scanner, the ability of each feature to discriminate the different textures was quantified. Despite differences observed among the scanner models, the overall percentage of discriminative features across scanners was >70%, with the smallest volume having the lowest percentage of discriminative features for all scanners. Stratification by feature class, still aggregating data for original and filtered images, showed statistical significance for the association between the percentage of discriminative features with VOI sizes for features classes GLCM, GLDM, and GLSZM on the first 1.5T scanner and for first-order and GLSZM classes on the second 1.5T scanner. Poorer results in terms of features’ discriminative ability were found for the 3T scanner. Focusing on original images only, the analysis of discriminative features stratified by feature class showed that the first-order and GLCM were robust to VOI size variations (>85% discriminative features for all sizes), while for the 1.5T scanners, the GLSZM and NGTDM feature classes showed a percentage of discriminative features >80% only for volumes no smaller than 3.3 cm3, and equal or larger than 7.4 cm3 for the GLRLM. As for the 3T scanner, only the GLSZM showed a percentage of discriminative features >80% for all volume sizes above 3.3 cm3. Analogous considerations were obtained for each filter, providing useful indications for feature selection in this clinical case. Similar studies should be replicated with suitably adapted phantoms to derive useful data for other clinical scenarios and imaging modalities.

Predictive Value of Delta-Radiomics Texture Features in 0.35 Tesla Magnetic Resonance Setup Images Acquired During Stereotactic Ablative Radiotherapy of Pancreatic Cancer

Purpose

The purpose of this work is to explore delta-radiomics texture features for predicting response using setup images of pancreatic cancer patients treated with magnetic resonance image guided (MRI-guided) stereotactic ablative radiotherapy (SBRT).

Methods

The total biological effective dose (BED) was calculated for 30 patients treated with MRI-guided SBRT that delivered physical doses of 30–60 Gy in three to five fractions. Texture features were then binned into groups based upon BED per fraction by dividing BED by the number of fractions. Delta-radiomics texture features were calculated after delivery of 20 Gy BED (BED20 features) and 40 Gy BED (BED40 features). A random forest (RF) model was constructed using BED20 and then BED40 features to predict binary outcome. During model training, the Gini Index, a measure of a variable’s importance for accurate prediction, was calculated for all features, and the two features that ranked the highest were selected for internal validation. The two features selected from each bin were used in a bootstrapped logistic regression model to predict response and performance quantified using the area under the receiver operating characteristic curve (AUC). This process was an internal validation analysis.

Results

After RF model training, the Gini Index was highest for gray-level co-occurrence matrix-based (GLCM) sum average, and neighborhood gray tone difference matrix-based (NGTDM) busyness for BED20 features and gray-level size zone matrix-based (GLSZM) large zones low gray-level emphasis and gray-level run length matrix-based (GLRLM) run percentage was selected from the BED40-based features. The mean AUC obtained using the two BED20 features was AUC = 0.845 with the 2.5 percentile and 97.5 percentile values ranging from 0.794 to 0.856. Internal validation of the BED40 delta-radiomics features resulted in a mean AUC = 0.567 with a 2.5 and 97.5 percentile range of 0.502–0.675.

Conclusion

Early changes in treatment quantified with the BED20 delta-radiomics texture features in low field images acquired during MRI-guided SBRT demonstrated better performance in internal validation than features calculated later in treatment. Further analysis of delta-radiomics texture analysis in low field MRI is warranted.

Radiomics of Patients with Locally Advanced Rectal Cancer: Effect of Preprocessing on Features Estimation from Computed Tomography Imaging

The purpose of this study was to investigate the effect of image preprocessing on radiomic features estimation from computed tomography (CT) imaging of locally advanced rectal cancer (LARC). CT images of 20 patients with LARC were used to estimate 105 radiomic features of 7 classes (shape, first-order, GLCM, GLDM, GLRLM, GLSZM, and NGTDM). Radiomic features were estimated for 6 different isotropic resampling voxel sizes, using 10 interpolation algorithms (at fixed bin width) and 6 different bin widths (at fixed interpolation algorithm). The intraclass correlation coefficient (ICC) and the coefficient of variation (CV) were calculated to assess the variability in radiomic features estimation due to preprocessing. A repeated measures correlation analysis was performed to assess any linear correlation between radiomic feature estimate and resampling voxel size or bin width. Reproducibility of radiomic feature estimate, when assessed through ICC analysis, was nominally excellent (ICC > 0.9) for shape features, good (0.75 < ICC ≤ 0.9) or moderate (0.5 < ICC ≤ 0.75) for first-order features, and moderate or poor (0 ≤ ICC ≤ 0.5) for textural features. A number of radiomic features characterized by good or excellent reproducibility in terms of ICC showed however median CV values greater than 15%. For most textural features, a significant (p < 0.05) correlation between their estimate and resampling voxel size or bin width was found. In CT imaging of patients with LARC, the estimate of textural features, as well as of first-order features to a lesser extent, is appreciably biased by preprocessing. Accordingly, this should be taken into account when planning clinical or research studies, as well as when comparing results from different studies and performing multicenter studies.

Overview of radiomics in prostate imaging and future directions

Recent advancements in imaging technology and analysis methods have led to an analytic framework known as radiomics. This framework extracts comprehensive high-dimensional features from imaging data and performs data mining to build analytical models for improved decision-support. Its features include many categories spanning texture and shape; thus, it can provide abundant information for precision medicine. Many studies of prostate radiomics have shown promising results in the assessment of pathological features, prediction of treatment response, and stratification of risk groups. Herein, we aimed to provide a general overview of radiomics procedures, discuss technical issues, explain various clinical applications, and suggest future research directions, especially for prostate imaging.

Inter-Reader Variability of Volumetric Subsolid Pulmonary Nodule Radiomic Features

OBJECTIVE

To evaluate the inter-observer consistency for subsolid pulmonary nodule radiomic features.

MATERIALS AND METHODS

Subsolid nodules were selected by reviewing radiology reports of CT examinations performed December 1, 2015 to April 1, 2016. Patients with CTs at two time points were included in this study. There were 55 patients with subsolid nodules, of whom 14 had two nodules. Of 69 subsolid nodules, 66 were persistent at the second time point, yielding 135 lesions for segmentation. Two thoracic radiologists and an imaging fellow segmented the lesions using a semi-automated volumetry algorithm (Syngo.via Vb20, Siemens). Coefficient of variation (CV) was used to assess consistency of 91 quantitative measures extracted from the subsolid nodule segmentations, including first and higher order texture features. The accuracy of segmentation was visually graded by an experienced thoracic radiologist. Influencing factors on radiomic feature consistency and segmentation accuracy were assessed using generalized estimating equation analyses and the Exact Mann-Whitney test.

RESULTS

Mean patient age was 71 (38-93 years), with 39 women and 16 men. Mean nodule volume was 1.39mL, range .03-48.2mL, for 135 nodules. Several radiomic features showed high inter-reader consistency (CV<5%), including entropy, uniformity, sphericity, and spherical disproportion. Descriptors such as surface area and energy had low consistency across inter-reader segmentations (CV>10%). Nodule percent solid component and attenuation influenced inter-reader variability of some radiomic features. The presence of contrast did not significantly affect the consistency of subsolid nodule radiomic features. Near perfect segmentation, within 5% of actual nodule size, was achieved in 68% of segmentations, and very good segmentation, within 25% of actual nodule size, in 94%. Morphologic features including nodule margin and shape (each p <0.01), and presence of air bronchograms (p = 0.004), bubble lucencies (p = 0.02) and broad pleural contact (p < 0.01) significantly affected the probability of near perfect segmentation. Stroke angle (p = 0.001) and length (p < 0.001) also significantly influenced probability of near perfect segmentation.

CONCLUSIONS

The inter-observer consistency of radiomic features for subsolid pulmonary nodules varies, with high consistency for several features, including sphericity, spherical disproportion, and first and higher order entropy, and normalized non-uniformity. Nodule morphology influences the consistency of subsolid nodule radiomic features, and the accuracy of subsolid nodule segmentation.

Role of hepatic metastatic lesion size on inter-reader reproducibility of CT-based radiomics features

OBJECTIVES

To evaluate the effect of hepatic metastatic lesion size on inter-reader reproducibility of CT-based 2D radiomics imaging features.

METHODS

Computerized tomography (CT) scans of 59 liver metastases from 34 patients with colorectal cancer were evaluated. Image segmentation was performed manually by three readers blinded to each other's results. For each radiomics feature, we created two datasets by sorting measurements according to size, i.e., (i) from the smallest to the largest lesion and (ii) from the largest to the smallest lesion. The Lin concordance correlation coefficient (CCC) was employed to analyze the reproducibility of radiomics features. In particular, the CCC was computed as a function of a number of elements in the dataset, by gradually adding lesions from each sorted dataset. To evaluate the effect of lesion size, we analyzed the difference between these two functions thus assessing the contribution of small and large lesions into the reproducibility of radiomics features.

RESULTS

Inter-reader reproducibility of CT-based 2D radiomics features assessed using Lin's CCC demonstrates tumor-size dependence. For example, the Lin CCC for GLCM contrast equals 0.88 (95% C.I. 0.84 to 0.92, p < 0.003) and could change by an additional + / - 0.06 depending on the presence of large or small lesions.

CONCLUSIONS

Groups of "large" and "small" lesions show different inter-reader reproducibility. The inter-reader reproducibility from the mixed group consisting of "large" and "small" lesions depends on the lesion-size distribution and can vary widely. This finding could partially explain variability in reproducibility of radiomics features in the literature.

KEY POINTS

• Groups of "large" and "small" lesions show different inter-reader reproducibility. • The inter-reader reproducibility from the mixed group consisting of "large" and "small" lesions depends on the lesion-size distribution.

Radiomics for Predicting Lung Cancer Outcomes Following Radiotherapy: A Systematic Review

Lung cancer's radiomic phenotype may potentially inform clinical decision-making with respect to radical radiotherapy. At present there are no validated biomarkers available for the individualisation of radical radiotherapy in lung cancer and the mortality rate of this disease remains the highest of all other solid tumours. MEDLINE was searched using the terms 'radiomics' and 'lung cancer' according to the Preferred Reporting Items for Systematic Reviews and Met-Analyses (PRISMA) guidance. Radiomics studies were defined as those manuscripts describing the extraction and analysis of at least 10 quantifiable imaging features. Only those studies assessing disease control, survival or toxicity outcomes for patients with lung cancer following radical radiotherapy ± chemotherapy were included. Study titles and abstracts were reviewed by two independent reviewers. The Radiomics Quality Score was applied to the full text of included papers. Of 244 returned results, 44 studies met the eligibility criteria for inclusion. End points frequently reported were local (17%), regional (17%) and distant control (31%), overall survival (79%) and pulmonary toxicity (4%). Imaging features strongly associated with clinical outcomes include texture features belonging to the subclasses Gray level run length matrix, Gray level co-occurrence matrix and kurtosis. The median cohort size for model development was 100 (15-645); in the 11 studies with external validation in a separate independent population, the median cohort size was 84 (21-295). The median number of imaging features extracted was 184 (10-6538). The median Radiomics Quality Score was 11% (0-47). Patient-reported outcomes were not incorporated within any studies identified. No studies externally validated a radiomics signature in a registered prospective study. Imaging-derived indices attained through radiomic analyses could equip thoracic oncologists with biomarkers for treatment response, patterns of failure, normal tissue toxicity and survival in lung cancer. Based on routine scans, their non-invasive nature and cost-effectiveness are major advantages over conventional pathological assessment. Improved tools are required for the appraisal of radiomics studies, as significant barriers to clinical implementation remain, such as standardisation of input scan data, quality of reporting and external validation of signatures in randomised, interventional clinical trials.

Radiomics in Oncology: A 10-Year Bibliometric Analysis

Objectives

To date, radiomics has been applied in oncology for over a decade and has shown great progress. We used a bibliometric analysis to analyze the publications of radiomics in oncology to clearly illustrate the current situation and future trends and encourage more researchers to participate in radiomics research in oncology.

Methods

Publications for radiomics in oncology were downloaded from the Web of Science Core Collection (WoSCC). WoSCC data were collected, and CiteSpace was used for a bibliometric analysis of countries, institutions, journals, authors, keywords, and references pertaining to this field. The state of research and areas of focus were analyzed through burst detection.

Results

A total of 7,199 pieces of literature concerning radiomics in oncology were analyzed on CiteSpace. The number of publications has undergone rapid growth and continues to increase. The USA and Chinese Academy of Sciences are found to be the most prolific country and institution, respectively. In terms of journals and co-cited journals, Scientific Reports is ranked highest with respect to the number of publications, and Radiology is ranked highest among co-cited journals. Moreover, Jie Tian has published the most publications, and Phillipe Lambin is the most cited author. A paper published by Gillies et al. presents the highest citation counts. Artificial intelligence (AI), segmentation methods, and the use of radiomics for classification and diagnosis in oncology are major areas of focus in this field. Test-retest statistics, including reproducibility and statistical methods of radiomics research, the relation between genomics and radiomics, and applications of radiomics to sarcoma and intensity-modulated radiotherapy, are frontier areas of this field.

Conclusion

To our knowledge, this is the first study to provide an overview of the literature related to radiomics in oncology and may inspire researchers from multiple disciplines to engage in radiomics-related research.

A deep look into radiomics

Radiomics is a process that allows the extraction and analysis of quantitative data from medical images. It is an evolving field of research with many potential applications in medical imaging. The purpose of this review is to offer a deep look into radiomics, from the basis, deeply discussed from a technical point of view, through the main applications, to the challenges that have to be addressed to translate this process in clinical practice. A detailed description of the main techniques used in the various steps of radiomics workflow, which includes image acquisition, reconstruction, pre-processing, segmentation, features extraction and analysis, is here proposed, as well as an overview of the main promising results achieved in various applications, focusing on the limitations and possible solutions for clinical implementation. Only an in-depth and comprehensive description of current methods and applications can suggest the potential power of radiomics in fostering precision medicine and thus the care of patients, especially in cancer detection, diagnosis, prognosis and treatment evaluation.

Radiomics of hepatocellular carcinoma: promising roles in patient selection, prediction, and assessment of treatment response

Radiomics refers to the process of conversion of conventional medical images into quantifiable data ("features") which can be further mined to reveal complex patterns and relationships between the voxels in the image. These high throughput features can potentially reflect the histology of biologic tissues at macroscopic and microscopic levels. Several studies have investigated radiomics of hepatocellular carcinoma (HCC) before and after treatment. HCC is a heterogeneous disease with diverse phenotypical and genotypical landscape. Due to this inherent heterogeneity, HCC lesions can manifest variable aggressiveness with different response to treatment options, including the newer targeted therapies. Hence, radiomics can be used as a potential tool to enable patient selection for therapies and to predict response to treatments and outcome. Additionally, radiomics may serve as a tool for earlier and more efficient assessment of response to treatment. Radiomics, radiogenomics, and radio-immunoprofiling and their potential roles in management of patients with HCC will be discussed and critically reviewed in this article.

Risk Score Generated from CT-Based Radiomics Signatures for Overall Survival Prediction in Non-Small Cell Lung Cancer

Simple Summary

Despite recent advancements in lung cancer treatment, individuals with lung cancer have a dismal 5-year survival rate of only 15%. In patients with non-small cell lung cancer (NSCLC), medical images have lately been employed as a valuable marker for predicting overall survival. The primary goal of this study was to develop a risk score based on computed tomography (CT) based radiomics feature signatures that may be used to predict survival in NSCLC patients. After analyzing 577 NSCLC patients from two data sets, we discovered that the risk score model’s prediction ability as a prognostic indicator was superior to other clinical indicators (age, stage, and gender), and the possibility of patient risk stratification with survival was evaluated using a risk score representation of 10 radiomics signatures. According to this study, the risk score generated using CT-based radiomics signatures promises to predict overall survival in NSCLC patients.

Abstract

This study aimed to create a risk score generated from CT-based radiomics signatures that could be used to predict overall survival in patients with non-small cell lung cancer (NSCLC). We retrospectively enrolled three sets of NSCLC patients (including 336, 84, and 157 patients for training, testing, and validation set, respectively). A total of 851 radiomics features for each patient from CT images were extracted for further analyses. The most important features (strongly linked with overall survival) were chosen by pairwise correlation analysis, Least Absolute Shrinkage and Selection Operator (LASSO) regression model, and univariate Cox proportional hazard regression. Multivariate Cox proportional hazard model survival analysis was used to create risk scores for each patient, and Kaplan–Meier was used to separate patients into two groups: high-risk and low-risk, respectively. ROC curve assessed the prediction ability of the risk score model for overall survival compared to clinical parameters. The risk score, which developed from ten radiomics signatures model, was found to be independent of age, gender, and stage for predicting overall survival in NSCLC patients (HR, 2.99; 95% CI, 2.27–3.93; p < 0.001) and overall survival prediction ability was 0.696 (95% CI, 0.635–0.758), 0.705 (95% CI, 0.649–0.762), 0.657 (95% CI, 0.589–0.726) (AUC) for 1, 3, and 5 years, respectively, in the training set. The risk score is more likely to have a better accuracy in predicting survival at 1, 3, and 5 years than clinical parameters, such as age 0.57 (95% CI, 0.499–0.64), 0.552 (95% CI, 0.489–0.616), 0.621 (95% CI, 0.544–0.689) (AUC); gender 0.554, 0.546, 0.566 (AUC); stage 0.527, 0.501, 0.459 (AUC), respectively, in 1, 3 and 5 years in the training set. In the training set, the Kaplan–Meier curve revealed that NSCLC patients in the high-risk group had a lower overall survival time than the low-risk group (p < 0.001). We also had similar results that were statistically significant in the testing and validation set. In conclusion, risk scores developed from ten radiomics signatures models have great potential to predict overall survival in NSCLC patients compared to the clinical parameters. This model was able to stratify NSCLC patients into high-risk and low-risk groups regarding the overall survival prediction.

Oropharyngeal cancer patient stratification using random forest based-learning over high-dimensional radiomic features

To improve risk prediction for oropharyngeal cancer (OPC) patients using cluster analysis on the radiomic features extracted from pre-treatment Computed Tomography (CT) scans. 553 OPC Patients randomly split into training (80%) and validation (20%), were classified into 2 or 3 risk groups by applying hierarchical clustering over the co-occurrence matrix obtained from a random survival forest (RSF) trained over 301 radiomic features. The cluster label was included together with other clinical data to train an ensemble model using five predictive models (Cox, random forest, RSF, logistic regression, and logistic-elastic net). Ensemble performance was evaluated over the independent test set for both recurrence free survival (RFS) and overall survival (OS). The Kaplan-Meier curves for OS stratified by cluster label show significant differences for both training and testing (p val < 0.0001). When compared to the models trained using clinical data only, the inclusion of the cluster label improves AUC test performance from .62 to .79 and from .66 to .80 for OS and RFS, respectively. The extraction of a single feature, namely a cluster label, to represent the high-dimensional radiomic feature space reduces the dimensionality and sparsity of the data. Moreover, inclusion of the cluster label improves model performance compared to clinical data only and offers comparable performance to the models including raw radiomic features.

Impact of Lesion Delineation and Intensity Quantisation on the Stability of Texture Features from Lung Nodules on CT: A Reproducible Study

Computer-assisted analysis of three-dimensional imaging data (radiomics) has received a lot of research attention as a possible means to improve the management of patients with lung cancer. Building robust predictive models for clinical decision making requires the imaging features to be stable enough to changes in the acquisition and extraction settings. Experimenting on 517 lung lesions from a cohort of 207 patients, we assessed the stability of 88 texture features from the following classes: first-order (13 features), Grey-level Co-Occurrence Matrix (24), Grey-level Difference Matrix (14), Grey-level Run-length Matrix (16), Grey-level Size Zone Matrix (16) and Neighbouring Grey-tone Difference Matrix (five). The analysis was based on a public dataset of lung nodules and open-access routines for feature extraction, which makes the study fully reproducible. Our results identified 30 features that had good or excellent stability relative to lesion delineation, 28 to intensity quantisation and 18 to both. We conclude that selecting the right set of imaging features is critical for building clinical predictive models, particularly when changes in lesion delineation and/or intensity quantisation are involved.

Intrinsic radiomic expression patterns after 20 Gy demonstrate early metabolic response of oropharyngeal cancers

PURPOSE

This study investigated the prognostic potential of intra-treatment PET radiomics data in patients undergoing definitive (chemo) radiation therapy for oropharyngeal cancer (OPC) on a prospective clinical trial. We hypothesized that the radiomic expression of OPC tumors after 20 Gy is associated with recurrence-free survival (RFS).

MATERIALS AND METHODS

Sixty-four patients undergoing definitive (chemo)radiation for OPC were prospectively enrolled on an IRB-approved study. Investigational ¹⁸ F-FDG-PET/CT images were acquired prior to treatment and 2 weeks (20 Gy) into a seven-week course of therapy. Fifty-five quantitative radiomic features were extracted from the primary tumor as potential biomarkers of early metabolic response. An unsupervised data clustering algorithm was used to partition patients into clusters based only on their radiomic expression. Clustering results were naïvely compared to residual disease and/or subsequent recurrence and used to derive Kaplan-Meier estimators of RFS. To test whether radiomic expression provides prognostic value beyond conventional clinical features associated with head and neck cancer, multivariable Cox proportional hazards modeling was used to adjust radiomic clusters for T and N stage, HPV status, and change in tumor volume.

RESULTS

While pre-treatment radiomics were not prognostic, intra-treatment radiomic expression was intrinsically associated with both residual/recurrent disease (P = 0.0256, χ 2 test) and RFS (HR = 7.53, 95% CI = 2.54-22.3; P = 0.0201). On univariate Cox analysis, radiomic cluster was associated with RFS (unadjusted HR = 2.70; 95% CI = 1.26-5.76; P = 0.0104) and maintained significance after adjustment for T, N staging, HPV status, and change in tumor volume after 20 Gy (adjusted HR = 2.69; 95% CI = 1.03-7.04; P = 0.0442). The particular radiomic characteristics associated with outcomes suggest that metabolic spatial heterogeneity after 20 Gy portends complete and durable therapeutic response. This finding is independent of baseline metabolic imaging characteristics and clinical features of head and neck cancer, thus providing prognostic advantages over existing approaches.

CONCLUSIONS

Our data illustrate the prognostic value of intra-treatment metabolic image interrogation, which may potentially guide adaptive therapy strategies for OPC patients and serve as a blueprint for other disease sites. The quality of our study was strengthened by its prospective image acquisition protocol, homogenous patient cohort, relatively long patient follow-up times, and unsupervised clustering formalism that is less prone to hyper-parameter tuning and over-fitting compared to supervised learning.

Observing deep radiomics for the classification of glioma grades

Deep learning is a promising method for medical image analysis because it can automatically acquire meaningful representations from raw data. However, a technical challenge lies in the difficulty of determining which types of internal representation are associated with a specific task, because feature vectors can vary dynamically according to individual inputs. Here, based on the magnetic resonance imaging (MRI) of gliomas, we propose a novel method to extract a shareable set of feature vectors that encode various parts in tumor imaging phenotypes. By applying vector quantization to latent representations, features extracted by an encoder are replaced with a fixed set of feature vectors. Hence, the set of feature vectors can be used in downstream tasks as imaging markers, which we call deep radiomics. Using deep radiomics, a classifier is established using logistic regression to predict the glioma grade with 90% accuracy. We also devise an algorithm to visualize the image region encoded by each feature vector, and demonstrate that the classification model preferentially relies on feature vectors associated with the presence or absence of contrast enhancement in tumor regions. Our proposal provides a data-driven approach to enhance the understanding of the imaging appearance of gliomas.

Uncertainty measurement of radiomics features against inherent quantum noise in computed tomography imaging

OBJECTIVES

Quantum noise is a random process in X-ray-based imaging systems. We addressed and measured the uncertainty of radiomics features against this quantum noise in computed tomography (CT) images.

METHODS

A clinical multi-detector CT scanner, two homogeneous phantom sets, and four heterogeneous samples were used. A solid tumor tissue removed from a male BALB/c mouse was included. We the placed phantom sets on the CT scanning table and repeated 20 acquisitions with identical imaging settings. Regions of interest were delineated for feature extraction. Statistical quantities-average, standard deviation, and percentage uncertainty-were calculated from these 20 repeated scans. Percentage uncertainty was used to measure and quantify feature stability against quantum noise. Twelve radiomics features were measured. Random noise was added to study the robustness of machine learning classifiers against feature uncertainty.

RESULTS

We found the ranges of percentage uncertainties from homogeneous soft tissue phantoms, homogeneous bone phantoms, and solid tumor tissue to be 0.01-2138%, 0.02-15%, and 0.18-16%, respectively. Overall, it was found that the CT features ShortRunHighGrayLevelEmpha (SRHGE) (0.01-0.18%), ShortRunLowGrayLevelEmpha (SRLGE) (0.01-0.41%), LowGrayLevelRunEmpha (LGRE) (0.01-0.39%), and LongRunLowGrayLevelEmpha (LRLGE) (0.02-0.66%) were the most stable features against the inherent quantum noise. The most unstable features were cluster shade (1-2138%) and max probability (1-16%). The impact of random noise to the prediction accuracy by different machine learning classifiers was found to be between 0 and 12%.

CONCLUSIONS

Twelve features were used for uncertainty measurements. The upper and lower bounds of percentage uncertainties were determined. The quantum noise effect on machine learning classifiers is model dependent.

KEY POINTS

• Quantum noise is a random process and is intrinsic to X-ray-based imaging systems. This inherent quantum noise creates unpredictable fluctuations in the gray-level intensities of image pixels. Extra cautions and further validations are strongly recommended when unstable radiomics features are selected by a predictive model for disease classification or treatment outcome prognosis. • We addressed and used the statistical quantity of percentage uncertainty to measure the uncertainty of radiomics features against the inherent quantum noise in computed tomography (CT) images. • A clinical multi-detector CT scanner, two homogeneous phantom sets, and four heterogeneous samples were used in the stability measurement. A solid tumor tissue removed from a male BALB/c mouse was included in the heterogeneous sample.

Repeatability of CBCT radiomic features and their correlation with CT radiomic features for prostate cancer

PURPOSE

Radiomic features of cone-beam CT (CBCT) images have potential as biomarkers to predict treatment response and prognosis for patients of prostate cancer. Previous studies of radiomic feature analysis for prostate cancer were assessed in a variety of imaging modalities, including MRI, PET, and CT, but usually limited to a pretreatment setting. However, CBCT images may provide an opportunity to capture early morphological changes to the tumor during treatment that could lead to timely treatment adaptation. This work investigated the quality of CBCT-based radiomic features and their relationship with reconstruction methods applied to the CBCT projections and the preprocessing methods used in feature extraction. Moreover, CBCT features were correlated with planning CT (pCT) features to further assess the viability of CBCT radiomic features.

METHODS

The quality of 42 CBCT-based radiomic features was assessed according to their repeatability and reproducibility. Repeatability was quantified by correlating radiomic features between 20 CBCT scans that also had repeated scans within 15 minutes. Reproducibility was quantified by correlating radiomic features between the planning CT (pCT) and the first fraction CBCT for 20 patients. Concordance correlation coefficients (CCC) of radiomic features were used to estimate the repeatability and reproducibility of radiomic features. The same patient dataset was assessed using different reconstruction methods applied to the CBCT projections. CBCT images were generated using 18 reconstruction methods using iterative (iCBCT) and standard (sCBCT) reconstructions, three convolution filters, and five noise suppression filters. Eighteen preprocessing settings were also considered.

RESULTS

Overall, CBCT radiomic features were more repeatable than reproducible. Five radiomic features are repeatable in > 97% of the reconstruction and preprocessing methods, and come from the gray-level size zone matrix (GLSZM), neighborhood gray-tone difference matrix (NGTDM), and gray-level-run length matrix (GLRLM) radiomic feature classes. These radiomic features were reproducible in > 9.8% of the reconstruction and preprocessing methods. Noise suppression and convolution filter smoothing increased radiomic features repeatability, but decreased reproducibility. The top-repeatable iCBCT method (iCBCT-Sharp-VeryHigh) is more repeatable than the top-repeatable sCBCT method (sCBCT-Smooth) in 64% of the radiomic features.

CONCLUSION

Methods for reconstruction and preprocessing that improve CBCT radiomic feature repeatability often decrease reproducibility. The best approach may be to use methods that strike a balance repeatability and reproducibility such as iCBCT-Sharp-VeryLow-1-Lloyd-256 that has 17 repeatable and eight reproducible radiomic features. Previous radiomic studies that only used pCT radiomic features have generated prognostic models of prostate cancer outcome. Since our study indicates that CBCT radiomic features correlated well with a subset of pCT radiomic features, one may expect CBCT radiomics to also generate prognostic models for prostate cancer.

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.

Identification of robust and reproducible CT-texture metrics using a customized 3D-printed texture phantom

OBJECTIVE

The objective of this study was to evaluate the robustness and reproducibility of computed tomography-based texture analysis (CTTA) metrics extracted from CT images of a customized texture phantom built for assessing the association of texture metrics to three-dimensional (3D) printed progressively increasing textural heterogeneity.

MATERIALS AND METHODS

A custom-built 3D-printed texture phantom comprising of six texture patterns was used to evaluate the robustness and reproducibility of a radiomics panel under a variety of routine abdominal imaging protocols. The phantom was scanned on four CT scanners (Philips, Canon, GE, and Siemens) to assess reproducibility. The robustness assessment was conducted by imaging the texture phantom across different CT imaging parameters such as slice thickness, field of view (FOV), tube voltage, and tube current for each scanner. The texture panel comprised of 387 features belonging to 15 subgroups of texture extraction methods (e.g., Gray-level Co-occurrence Matrix: GLCM). Twelve unique image settings were tested on all the four scanners (e.g., FOV125). Interclass correlation two-way mixed with absolute agreement (ICC3) was used to assess the robustness and reproducibility of radiomic features. Linear regression was used to test the association between change in radiomic features and increased texture heterogeneity. Results were summarized in heat maps.

RESULTS

A total of 5612 (23.2%) of 24 090 features showed excellent robustness and reproducibility (ICC ≥ 0.9). Intensity, GLCM 3D, and gray-level run length matrix (GLRLM) 3D features showed best performance. Among imaging variables, changes in slice thickness affected all metrics more intensely compared to other imaging variables in reducing the ICC3. From the analysis of linear trend effect of the CTTA metrics, the top three metrics with high linear correlations across all scanners and scanning settings were from the GLRLM 2D/3D and discrete cosine transform (DCT) texture family.

CONCLUSION

The choice of scanner and imaging protocols affect texture metrics. Furthermore, not all CTTA metrics have a linear association with linearly varying texture patterns.

MRI Texture Analysis for the Prediction of Stereotactic Radiosurgery Outcomes in Brain Metastases from Lung Cancer

This study aims to evaluate the utility of texture analysis in predicting the outcome of stereotactic radiosurgery (SRS) for brain metastases from lung cancer. From 83 patients with lung cancer who underwent SRS for brain metastasis, a total of 118 metastatic lesions were included. Two neuroradiologists independently performed magnetic resonance imaging (MRI)-based texture analysis using the Imaging Biomarker Explorer software. Inter-reader reliability as well as univariable and multivariable analyses were performed for texture features and clinical parameters to determine independent predictors for local progression-free survival (PFS) and overall survival (OS). Furthermore, Harrell’s concordance index (C-index) was used to assess the performance of the independent texture features. The primary tumor histology of small cell lung cancer (SCLC) was the only clinical parameter significantly associated with local PFS in multivariable analysis. Run-length non-uniformity (RLN) and short-run emphasis were the independent texture features associated with local PFS. In the non-SCLC (NSCLC) subgroup analysis, RLN and local range mean were associated with local PFS. The C-index of independent texture features was 0.79 for the all-patients group and 0.73 for the NSCLC subgroup. In conclusion, texture analysis on pre-treatment MRI of lung cancer patients with brain metastases may have a role in predicting SRS response.

Harmonization strategies for multicenter radiomics investigations

Carrying out large multicenter studies is one of the key goals to be achieved towards a faster transfer of the radiomics approach in the clinical setting. This requires large-scale radiomics data analysis, hence the need for integrating radiomic features extracted from images acquired in different centers. This is challenging as radiomic features exhibit variable sensitivity to differences in scanner model, acquisition protocols and reconstruction settings, which is similar to the so-called 'batch-effects' in genomics studies. In this review we discuss existing methods to perform data integration with the aid of reducing the unwanted variation associated with batch effects. We also discuss the future potential role of deep learning methods in providing solutions for addressing radiomic multicentre studies.

Impact of quantization algorithm and number of gray level intensities on variability and repeatability of low field strength magnetic resonance image-based radiomics texture features

PURPOSE

The purpose of this work was to investigate the impact of quantization preprocessing parameter selection on variability and repeatability of texture features derived from low field strength magnetic resonance (MR) images.

METHODS

Texture features were extracted from low field strength images of a daily image QA phantom with four texture inserts. Feature variability over time was quantified using all combinations of three quantization algorithms and four different numbers of gray level intensities. In addition, texture features were extracted using the same combinations from the low field strength MR images of the gross tumor volume (GTV) and left kidney of patients with repeated set up scans. The impact of region of interest (ROI) preprocessing on repeatability was investigated with a test-retest study design.

RESULTS

The phantom ROIs quantized to 64 Gy level intensities using the histogram equalization method resulted in the greatest number of features with the least variability. There was no clear method that resulted in the highest repeatability in the GTV or left kidney. However, eight texture features extracted from the GTV were repeatable regardless of ROI processing combination.

CONCLUSION

Low field strength MR images can provide a stable basis for texture analysis with ROIs quantized to 64 Gy levels using histogram equalization, but there is no clear optimal combination for repeatability.

Impact of inter-reader contouring variability on textural radiomics of colorectal liver metastases

BACKGROUND

Radiomics is expected to improve the management of metastatic colorectal cancer (CRC). We aimed at evaluating the impact of liver lesion contouring as a source of variability on radiomic features (RFs).

METHODS

After Ethics Committee approval, 70 liver metastases in 17 CRC patients were segmented on contrast-enhanced computed tomography scans by two residents and checked by experienced radiologists. RFs from grey level co-occurrence and run length matrices were extracted from three-dimensional (3D) regions of interest (ROIs) and the largest two-dimensional (2D) ROIs. Inter-reader variability was evaluated with Dice coefficient and Hausdorff distance, whilst its impact on RFs was assessed using mean relative change (MRC) and intraclass correlation coefficient (ICC). For the main lesion of each patient, one reader also segmented a circular ROI on the same image used for the 2D ROI.

RESULTS

The best inter-reader contouring agreement was observed for 2D ROIs according to both Dice coefficient (median 0.85, interquartile range 0.78-0.89) and Hausdorff distance (0.21 mm, 0.14-0.31 mm). Comparing RF values, MRC ranged 0-752% for 2D and 0-1567% for 3D. For 24/32 RFs (75%), MRC was lower for 2D than for 3D. An ICC > 0.90 was observed for more RFs for 2D (53%) than for 3D (34%). Only 2/32 RFs (6%) showed a variability between 2D and circular ROIs higher than inter-reader variability.

CONCLUSIONS

A 2D contouring approach may help mitigate overall inter-reader variability, albeit stable RFs can be extracted from both 3D and 2D segmentations of CRC liver metastases.

A systematic review and meta-analysis of the prognostic value of radiomics based models in non-small cell lung cancer treated with curative radiotherapy

BACKGROUND AND PURPOSE

Radiomics allows extraction of quantifiable features from imaging. This study performs a systematic review and meta-analysis of the performance of radiomics based prognostic models in non-small cell lung cancer (NSCLC).

MATERIALS AND METHODS

A literature review was performed following PRISMA guidelines. Medline, EMBASE and Cochrane databases were searched for articles investigating radiomics features predictive of overall survival (OS) in NSCLC treated with curative intent radiotherapy. A random-effects meta-analysis of Harrell's Concordance Index (C-index) was performed on the performance of radiomics models.

RESULTS

Of the 2746 articles retrieved, 40 studies of 55 datasets and 6223 patients were eligible for inclusion in the systematic review. There was significant heterogeneity in the methodology for feature selection and model development. Twelve datasets reported the C-index of radiomics based models in predicting OS and were included in the meta-analysis. The C-index random effects estimate was 0.57 (95% CI 0.53-0.62). There was significant heterogeneity (I² = 70.3%).

CONCLUSIONS

Based on this review, radiomics based models for lung cancer have to date demonstrated modest prognostic capabilities. Future research should consider using standardised radiomics features, robust feature selection and model development, and deep learning techniques, absolving the need for pre-defined features, to improve imaging-based models.

The role of radiomics in prostate cancer radiotherapy

"Radiomics," as it refers to the extraction and analysis of a large number of advanced quantitative radiological features from medical images using high-throughput methods, is perfectly suited as an engine for effectively sifting through the multiple series of prostate images from before, during, and after radiotherapy (RT). Multiparametric (mp)MRI, planning CT, and cone beam CT (CBCT) routinely acquired throughout RT and the radiomics pipeline are developed for extraction of thousands of variables. Radiomics data are in a format that is appropriate for building descriptive and predictive models relating image features to diagnostic, prognostic, or predictive information. Prediction of Gleason score, the histopathologic cancer grade, has been the mainstay of the radiomic efforts in prostate cancer. While Gleason score (GS) is still the best predictor of treatment outcome, there are other novel applications of quantitative imaging that are tailored to RT. In this review, we summarize the radiomics efforts and discuss several promising concepts such as delta-radiomics and radiogenomics for utilizing image features for assessment of the aggressiveness of prostate cancer and its outcome. We also discuss opportunities for quantitative imaging with the advance of instrumentation in MRI-guided therapies.