Skip to main content

18F-FET PET radiomics-based survival prediction in glioblastoma patients receiving radio(chemo)therapy



Quantitative image analysis based on radiomic feature extraction is an emerging field for survival prediction in oncological patients. 18F-Fluorethyltyrosine positron emission tomography (18F-FET PET) provides important diagnostic and grading information for brain tumors, but data on its use in survival prediction is scarce. In this study, we aim at investigating survival prediction based on multiple radiomic features in glioblastoma patients undergoing radio(chemo)therapy.


A dataset of 37 patients with glioblastoma (WHO grade 4) receiving radio(chemo)therapy was analyzed. Radiomic features were extracted from pre-treatment 18F-FET PET images, following intensity rebinning with a fixed bin width. Principal component analysis (PCA) was applied for variable selection, aiming at the identification of the most relevant features in survival prediction. Random forest classification and prediction algorithms were optimized on an initial set of 25 patients. Testing of the implemented algorithms was carried out in different scenarios, which included additional 12 patients whose images were acquired with a different scanner to check the reproducibility in prediction results.


First order intensity variations and shape features were predominant in the selection of most important radiomic signatures for survival prediction in the available dataset. The major axis length of the 18F-FET-PET volume at tumor to background ratio (TBR) 1.4 and 1.6 correlated significantly with reduced probability of survival. Additional radiomic features were identified as potential survival predictors in the PTV region, showing 76% accuracy in independent testing for both classification and regression.


18F-FET PET prior to radiation provides relevant information for survival prediction in glioblastoma patients. Based on our preliminary analysis, radiomic features in the PTV can be considered a robust dataset for survival prediction.


Glioblastoma is the most aggressive primary brain tumor with a median survival of less than two years [1,2,3,4]. Combining surgery, chemotherapy and radiotherapy is the current standard of care, but prognosis remains poor because of the heterogeneous and complex pathogenesis [2]. Therefore, survival prediction based on images acquired before treatment may represent an important step towards a more personalized approach. In recent years, a lot of research has been done in this area with a strong focus on quantitative medical image analysis [5].

There is an increasing interest towards the use of 18F-Fluorethyltyrosine (18F-FET) Positron Emission Tomography (PET) imaging for the characterization of aggressive brain tumors such as glioblastoma. 18F-FET PET images were reported to provide valuable information, able to effectively support treatment response assessment in glioblastoma [6,7,8,9,10]. In the framework of the recent randomized multicenter phase II trial ARTE, higher 18F-FET PET TBR (tumor-to-background ratio) on pre-treatment scans and persistent 18F-FET PET signal of longer contrast-enhancing tumor was found to associate with inferior overall survival [10]. However, dedicated studies are needed to assess the predictive power related to 18F-FET PET radiomics for survival prediction in glioblastoma. The effective use of radiomics to predict overall survival as an indicator of treatment response comes with a large variability of different strategies regarding image acquisition, tumor segmentation, image processing and feature extraction, data handling and implementation of machine learning models [5, 11, 12]. For example, reliable tumor segmentation methods should reduce the inherent intra- and inter-observer variability [13]. At the same time, segmented contours should be readily available to facilitate their integration in the treatment planning process. In addition, site-specific image acquisition protocols introduce variations in imaging scanners, image reconstruction and even SUV intensities, which can have an effect on the extracted radiomic features [12, 14, 15]. This calls for specific studies to evaluate the reproducibility of radiomic features, in order to maximize their use in survival prediction models.

The objective of this work was to identify key radiomic features in 18F-FET PET images acquired before treatment for glioblastoma patients undergoing radiochemotherapy, and to assess their reproducibility, aiming at consistent prediction accuracy of overall survival.


Radiomic features were extracted from four different volumes of interest (VOIs) to identify the VOIs with the most reliable information in terms of survival prediction. After features were extracted, principal component analysis (PCA) was employed to reduce dimensionality of the feature space, working either as data reduction or feature selection method. Four different datasets were generated with the PCA and used as input data for the machine learning algorithms. Overall, three different prediction models were implemented: (1) a univariate Cox regression model to predict survival time, providing preliminary investigation of potentially relevant features, (2) a random forest classification model to predict overall survival at 1 year, (3) a random forest regression model to predict survival time.

Clinical dataset

The retrospective study included 37 patients with glioblastoma (WHO grade 4) diagnosed between 2009 and 2017. Survival time was defined, in days, as the time between the day of treatment planning CT acquisition and the day of decease. Due to the high malignancy of glioblastoma and the related poor survival prognosis, the event of decease occurred in all of the cases. The following inclusion criteria were applied for patient selection:

  • Evidence of macroscopic tumor before the start of treatment planning

  • PET image acquired not more than 3 weeks before the treatment planning CT to avoid significant anatomical changes and tumor growth/shrinkage between PET imaging and start of treatment

  • IDH-wildtype glioblastoma, in cases of available IDH mutational status

  • Follow-up of at least one year.

18F-FET PET images were acquired at the Department of Nuclear Medicine, University Hospital of Munich prior to primary radiation treatment. 18F-FET PET images were used to delineate the tumor based on standardized uptake values (SUVs) and subsequently extract radiomic features within the selected volumes of interest (VOIs). Planning target volumes (PTVs) were delineated within the treatment planning software Oncentra External Beam (version 4.5, Nucletron, Elekta AB, Sweden). When available (32 out of 37 patients), the reconstructed treatment planning CTs and structures were restored from the archive using the treatment planning software.

18F-FET PET images acquired with a Siemens ECAT Exact HR+ (Siemens Healthineers, Erlangen, Germany) scanner according to the standard protocol of the department, described in detail previously [16], were used as initial training set (N = 25 cases). A transmission scan was performed during 15 min with a 68Ge rod source. Thereafter, approximately 180 MBq of 18F-FET were injected as an intravenous bolus. Once the patient was in position, a 40 min dynamic emission recording in 3D mode consisting of 16 frames was started. The ordered subset expectation maximization (OSEM) algorithm was used for image reconstruction slice per slice (OSEM-2D). The OSEM-2D reconstruction consisted of 6 iterations and used 16 projection subsets. The voxel size in images reconstructed with the described OSEM-2D algorithm was 2.035 × 2.035 × 2.425 mm3, with a resulting image matrix of size 128 × 128 × 63.

The remaining 12 patients were acquired in list mode with a Biograph 64 PET/CT (Siemens Healthineers, Erlangen, Germany), and OSEM-3D was used to reconstruct images. OSEM-3D images were reconstructed with 4 iterations, 21 subsets, a 5 mm Gaussian post-reconstruction filter and standard corrections (attenuation, random and scattered coincidences, dead time decay). The reconstruction matrix was set to 336 × 336 × 109, with voxel size 1.018 × 1.018 × 2.027 mm3.

The key features of the available clinical dataset are summarized in Table 1.

Table 1 Features of the clinical dataset

Feature extraction

Four different volumes of interest (VOIs) were used to compare the prognostic power of the radiomics analysis when limited to a specific region. The PMOD software (PMOD Technologies, Zurich, Switzerland) was used to generate PET-based VOIs and to register PTV contours to the PET images. PET-based VOIs were defined using a semi-automatic segmentation method, as described by Unterrainer et al. [16]. A tumor-to-background ratio (TBR) threshold of 1.6 has been considered to be optimal to differentiate tumor tissue and surrounding healthy tissue [16, 17]. Besides 1.6, two additional TBR thresholds were chosen, one smaller value of 1.4 and one larger value of 1.8. This procedure resulted in a set of three PET-based VOIs for every patient, further on referred to as voi14, voi16 and voi18 corresponding to the TBR threshold 1.4, 1.6 and 1.8, respectively.

The fourth VOI was generated using the PTV, which was defined by experienced radiation oncologists during the treatment planning process. The gross tumor volume (GTV) was contoured based on the contrast enhanced T1 weighted MRI, complementing the treatment planning CT. T2 and Fluid Attenuated Inversion Recovery (FLAIR) MRI images were used in addition if available. The clinical target volume (CTV) was then created by adding a 2 cm margin to the GTV [18]. The final planning target volume (PTV) was obtained by adding 0.3–0.5 cm to the CTV. As the goal was to extract radiomics features from the PET image, the PTV structures had to be registered to the PET images. The PET image was loaded into PMOD, serving as reference image, whereas the treatment planning CT image was iteratively registered via rigid matching. The resulting transformation matrix was applied to the PTV contour to be overlaid onto the PET image. This provided PET intensity values within the PTV with no need to interpolate raw image intensities. The workflow of the image registration procedure is depicted in Fig. 1.

Fig. 1
figure 1

Workflow of the image registration procedure that was implemented to overlay the PTV and PET-based VOIs onto the FET-PET image

The PTV, when defined on the treatment planning CT, may include bone and other structures with high CT numbers that do not provide any information about tracer uptake in tumor tissue and surrounding brain tissue. As those regions would cause noise in the radiomic features and impair identification of intensity related patterns in tumor regions, these structures were excluded from the PTV. As shown schematically in Fig. 1, the PTV was much larger compared to the metabolically active tumor, consistently with previously reported studies [19]. This provided a way to assess the role of areas featuring a low intensity 18F-FET PET signal, included in the PTV but outside the metabolically active area (TBR-based VOIs).

Image intensities were rebinned using a fixed bin width approach [11, 20]. To determine the appropriate bin width, the mean TBR intensity range of the initial 25 patients was calculated from voi14 and subsequently divided by the number of bins. Such a number was selected aiming at optimal feature robustness, as investigated by Tsujikawa et al. for PET imaging [21], and therefore set to 64. According to this strategy, a bin width of 0.04 was obtained for TBR-based resampling, by dividing the mean TBR intensity range (voi14) in 64 bins.

Following rebinning, the python package pyradiomics was used to extract quantitative image features for the selected VOIs (i.e. voi14, voi16, voi18 and PTV) in the 18F-FET PET image [22, 23]. In total 107 radiomic features were extracted from the 3D PET images, including the following feature classes: First Order (FO), Shape (SH), Gray Level Co-occurrence Matrix (GLCM), Gray Level Size Zone Matrix (GLSZM), Gray Level Run Length Matrix (GLRLM), Neighbouring Gray Tone Difference Matrix (NGTDM), Gray Level Dependence Matrix (GLDM).

Dimensionality reduction and feature selection

Unsupervised methods were applied to reduce the data dimensionality and perform feature selection. For both tasks we applied PCA based on the singular value decomposition (SVD) of the feature set, as provided within the scikit-learn package (version 0.22.2) [24,25,26]. PCA was utilized following the standardization of radiomic features, which resulted in feature values with zero mean and unit variance.

Four different types of feature datasets were considered:

  • Top: the dataset consisted of 10 features per principal component and was generated using PCA-based feature selection: the features with the highest PCA loadings were retained, focusing on the components that explained 90–95% of the total variance in the data

  • PC: the dataset consisted of the principal component values of each patient, calculated with the entire feature set of 107 features, aiming at the same total variance applied for the Top features

  • SH-FO: the dataset consisted of all the shape and first order features (32 features overall)

  • PC(SH-FO): the dataset consisted of the principal component values of each patient, calculated with shape and first order features (32 overall), with the same goal in terms of total variance explained.

Identification of relevant radiomic features

The extracted radiomic features were analyzed to highlight the relevant ones for the purpose of survival prediction.

The Cox proportional hazard model was applied to identify relevant radiomic features in survival prediction based on the calculated hazard ratios. The Cox regression was implemented using the lifelines package in python (version 0.24.4) [27]. Univariate Cox models were fitted to the four feature datasets and potentially relevant features were identified based on the reported p-value of the Wald statistics. A threshold corresponding to 95% significance (i.e. 5% p-value) was selected to classify whether the corresponding radiomic features had a statistically relevant impact on survival prediction.

Survival prediction based on multiple features

We used random forests for both classification and regression analysis based on multiple radiomic features. Random forest is an ensemble method, which means that it combines a number of base estimators built with decision trees in order to reduce variance and improve generalizability [28, 29]. Random forest algorithms were implemented with the python package Scikit-learn [29]. Survival prediction was applied for classification (categorize survival at one-year follow-up) and regression (prediction of survival time). In Scikit-learn, the random forest algorithm for classification provided the use of different class weights to adjust for the imbalance in class frequencies.

To select the best hyper-parameter settings in the models, a coarse to fine strategy was applied using a randomized search cross-validation on training data. To reduce computation time, the randomized search cross-validation was applied twice: (1) in a very broad range of hyper-parameters to identify the order of magnitude of the best performing hyper-parameters and (2) once the magnitude was assessed, multiple values within this magnitude were included in the final parameter grid.

The performance in survival prediction was quantified on the optimal model (i.e. with the determined hyperparameters): we rated survival classification at 12 months based on the Area Under the Receiver Operating Characteristic Curve (AUC), and survival time prediction based on the Concordance Index (C-index). We considered different testing scenarios, where data from the 25 OSEM-2D cases were mixed with the 12 OSEM-3D cases:

  1. 1.

    Only OSEM-3D data in the test dataset, only OSEM-2D data in the training dataset

  2. 2.

    A fivefold train and test cross-validation using the entire dataset (OSEM-2D and OSEM-3D data combined).


Patients and treatment

Patients were treated with conventionally fractionated radiochemotherapy with temozolomide (TMZ) based on the EORTC26981/22981-NCIC CE3 protocol [30] in 34 cases, in one case with conventionally fractionated radiochemotherapy with TMZ and lomustine according to the NOA-09 protocol [31], in one case with hypofractionated radiotherapy only and in one case with conventionally fractionated radiotherapy only due to a contraindication for chemotherapy. Radiotherapy dosage was 30 × 2 Gy in 32 cases, 29 × 2 Gy due to discontinuation in one case, 33 × 1.8 Gy in two cases, 12 × 1.8 Gy/2.1 Gy on a simultaneous integrated boost (SIB) due to discontinuation in one case and with 10 × 3 Gy plus 1 × 4 Gy in one case.

Identification of relevant radiomic features

When considering the initial 25 patients, a total number of 5–7 principal components turned out to be sufficient to explain 90–95% of the total variance in the data. The Top dataset included the 10 features with the highest loadings for the first 5–7 principal components: Fig. 2 shows the features that appeared most frequently in the Top dataset across the 4 VOIs (voi14, voi16, voi18 and PTV). Shape and first order features were predominantly selected (87.5% and 55.6% average selection frequency, respectively) compared to all other feature types (37.3% average selection frequency). This indicates that shape and first order features have larger variations across the analyzed glioblastoma patients, which might be relevant for outcome prediction.

Fig. 2
figure 2

Circle plot showing the frequency of feature selection in the Top dataset. The height of each color bar in the circular plot is proportional to the selection frequency across the 4 VOIs. Color codes denote different radiomic features types, as indicated by text labels in the inner part of the circle

When analyzing the statistical relevance of single radiomic features via the Cox proportional hazard model, very few radiomic features turned out to be significant in the Top dataset. Figure 3 depicts graphically which radiomic features were classified as significant in the different regions at 95% confidence in univariate Cox proportional hazard models.

Fig. 3
figure 3

Statistically relevant features from the univariate Cox model fit for the Top radiomics dataset. All of the listed features were associated to a hazard ratio larger than 1, meaning that increasing values reduce the survival probability

Figure 4 shows the effect on survival of the shape feature Major Axis Length in voi14 and voi16 which resulted in significantly different prediction performance for the univariate Cox proportional hazard model. The log-rank test comparison between the Kaplan–Meier curves, where patients were stratified according to the median value, exhibited a p-value smaller than 0.05. This indicates the ability of the radiomic feature Major Axis Length to stratify patients in different risk classes, connected to a significantly different probability of survival at 95% statistical confidence.

Fig. 4
figure 4

Kaplan–Meier survival curves showing the effect of the shape feature Major Axis Length in voi14 (left panel) and voi16 (right panel). Patient stratification based on the median value is depicted, with the corresponding p-value of the log-rank test

Survival prediction based on multiple features

Table 2 reports the average results for the two testing scenarios. The best performance was obtained when considering the PTV as VOI for independent testing on OSEM-3D images. Conversely, results for the CV case were in general lower than independent testing. This indicates that features selected on OSEM-2D images can adequately predict survival in OSEM-3D, whereas data mixing in training does not allow the implemented models to achieve satisfactory prediction performance.

Table 2 Results on the testing experiments, expressed as average over all scenarios (1–4)


We analyzed the reproducibility of radiomic features extracted from 18F-FET PET images, and their use for survival analysis in glioblastoma patients. Two different patient cohorts were available which included cases with PET images reconstructed via an OSEM-2D algorithm, as well as more recent cases with OSEM-3D reconstruction. The OSEM-2D cohort served as a homogeneous training dataset, providing a consistent source of data to optimize the performance of survival prediction models based on multiple radiomic features. The random forests survival models, featuring the parameters optimized on the first patient cohort, were then tested in multiple scenarios including newer data. This provided a framework to check the reproducibility of radiomic features extracted from 18F-FET PET images to achieve accurate survival prediction in glioblastoma patients.

When analyzing the role of single radiomic features, shape and first order features extracted from 18F-FET PET images were highlighted as the most relevant in the analyzed dataset. This was confirmed both in unsupervised PCA feature selection (Fig. 2) and in the analysis of hazard ratios in univariate Cox proportional models for PET-based VOIs (Fig. 3). As visible in Fig. 3, different features were deemed significant when analyzing the univariate Cox proportional models in multiple regions (VOIs). Univariate Cox models were used to assess significant variables, instead of multivariate analysis. This choice was justified by the intrinsic limitations in Cox proportional models for reduced number of events per variable, as discussed in [32].

A closer look at these results highlights that increasing values of the shape feature Major Axis Length (i.e. the maximal extension of the biologically active area in the 18F-FET PET image) in voi14 and voi16 are significantly correlated to a reduced overall survival time. This means that larger regions of biologically active tumor in the 18F-FET PET image, when the image is thresholded at 1.4 or 1.6 TBR, are clear indicators of lower survival probability. When extending the analysis to the PTV, the Elongation was also found to be significantly correlated to survival. The Elongation is defined as the square root of the ratio between the minor and major axis of the PTV, with values ranging between 0 (maximal elongation) and 1 (sphere-like shape) [22]. As increasing values are associated to poorer survival, this means that spherical PTVs seem less favorable in terms of survival compared to elongated PTVs. It should be stressed here that Elongation values in our dataset were relatively high (above 0.57), meaning that the minor axis of the PTV is at least 33% of the major axis. Specific further studies are therefore required to check performance in case of elongated PTVs, i.e. corresponding to small Elongation values.

Higher order features were not linked to significant survival prediction in TBR-based VOIs. Conversely, increasing values of GLCM Joint Entropy in the PTV, which encodes variability in intensity, were shown to decrease survival. The existence of multiple significant features calls for survival prediction models that can better handle multiple radiomics signatures simultaneously, such as the implemented random forest algorithms.

The generalization capabilities were analyzed by investigating the reproducibility of radiomics-based predictions in two testing scenarios, where images reconstructed with the OSEM-3D algorithm were also included. In independent testing, significantly better results were obtained when considering the PTV as VOI (Table 2). In this setting, feature sets that included shape and first order 18F-FET PET features obtained more consistent results, thus confirming the findings of univariate Cox models. The Cox regression models resulted in general lower performance than the implemented Random Forest predictors. When testing the univariate Cox models on OSEM-3D data, the C-index for the Major Axis Length in VOI16 and VOI14 reached 0.58 and 0.56, respectively.

Conversely, the performance in CV when mixing different images (OSEM-2D vs. OSEM-3D) were inferior, with the exception of the PC(SH-FO) dataset for survival classification in the PTV (0.7 AUC). This can be partly explained by the lower degree of reproducibility between VOIs segmented in the OSEM-2D reconstructed images versus the OSEM-3D ones. Results from independent testing on OSEM-3D images indicated that features selected on OSEM-2D images can adequately predict survival on different images. Overall best testing performance was found for the SH-FO dataset in classification for the PTV, thus confirming the importance of shape and first order features.

The optimal performance in the PTV shows that also areas of the 18F-FET PET image at lower image intensities, i.e. outside of the high intensity region at 1.4–1.6 TBR, may play a key role in radiomics-based survival prediction. This also strengthens the relevance of expert-based contouring for PTV delineation, which can potentially highlight relevant areas of low/medium biological activity, due to the intrinsic relevance in survival prediction.

A comparative analysis of our results with other studies is challenging, due to the lack of extensive literature on the topic. 18F-FET PET images have been mostly used for diagnostic purposes, and/or to differentiate progression and treatment-related changes [5, 9]. Carles et al. [8] have reported an analysis of relevant 18F-FET PET features in recurrent glioblastoma patients, including also overall survival as clinical endpoint. More complex features were highlighted as significant for recurrent glioblastoma patients, whereas our results suggest that shape and first order features are more predictive for overall survival in primary glioblastoma. Conversely, best results in our study were achieved for the PTV, where the GLCM Joint Entropy turned out to be significant in univariate analysis. This confirms that higher order features are indeed relevant for primary glioblastoma patients, especially when considered within the PTV region.

MR images have been used more frequently to predict survival in a radiomics-based approach. The performance reported in our study based on 18F-FET PET images are in line with previously reported outcomes on MRI-based radiomics in glioblastoma [33,34,35,36,37,38]. These studies include also more recent extensions from traditional radiomics to deep learning models, designed to include traditional risk factors as potential predictors [35]. Therefore, 18F-FET PET images are suggested as valuable candidates for the prediction of overall survival in primary glioblastoma patients. The complimentary use of MRI and 18F-FET PET images has been recently shown in recurrent glioblastoma patients to derive independent biomarkers of response to treatment [38].

The main limitation in our study is the reduced number of patients, also in reason of the unique characteristics of the analyzed patient cohort, which consists solely of glioblastoma patients receiving 18F-FET PET prior to primary radio(chemo)therapy. This motivated the use of stringent criteria for feature selection, but hindered a more extensive analysis of the relevance of single features in survival prediction. For this reason, the univariate analysis was restricted to the Top features dataset, which might have underestimated the relevance of complex textural features, as these latter were under-represented in the PCA feature selection procedure (Fig. 2). Reported results indeed confirm the appropriateness of such a selection, as PCA including only SH-FO features outperformed PCA including more complex features (Table 2). As further limitation, no comparison with MRI-based radiomic features has been carried out in this work. Different studies have investigated the use of quantitative features from MRI and 18F-FET PET for treatment response and survival [6, 10, 38], highlighting the complimentary role of these imaging modalities. A comparative radiomics-based analysis is therefore suggested to highlight the need of dedicated imaging protocols for glioblastoma survival prediction.


The findings of feature analysis and prediction performance based on 18F-FET PET images show that shape features in voi14, voi16, and PTV are significant predictors of overall survival in glioblastoma patients. When radiomics-based prediction is enlarged to the PTV a more robust prediction of survival is achieved, with 76% accuracy in independent testing both for classification and regression analyses. These preliminary results should be extended to a larger patient cohort for further validation of our findings.

Availability of data and materials

Due to privacy and ethical concerns, the source data cannot be made available.



Area under the curve


Computerized tomography




Fluid attenuated inversion recovery


First order


Gray level co-occurrence matrix


Gray level dependence matrix


Gray level run length matrix


Gray level size zone matrix


Isocitrate DeHydrogenase


Mean absolute error


Magnetic resonance imaging


Mean square error


Neighbouring gray tone difference matrix


Ordered subset expectation maximization


Principal component


Principal component analysis


Positron emission tomography


Planning target volume


Standard deviation




Standardized uptake value


Tumor to BackgRound




Volume of interest


  1. Wirsching HG, Galanis E, Weller M. Glioblastoma. In: Berger MS, Weller M, editors. Handbook of clinical neurology, vol. 134. Amsterdam: Elsevier; 2016. p. 381–97.

    Google Scholar 

  2. Hanif F, Muzaffar K, Perveen K, Malhi SM, Simjee SU. Glioblastoma multiforme: a review of its epidemiology and pathogenesis through clinical presentation and treatment. Asian Pac J Cancer Prev. 2017;18(1):3–9.

    PubMed  PubMed Central  Google Scholar 

  3. Tamimi AF, Juweid M. Epidemiology and outcome of glioblastoma. In: De Vleeschouwer S, editor. Glioblastoma [internet]. Brisbane: Codon Publications; 2017.

    Chapter  Google Scholar 

  4. Fernandes C, Costa A, Osório L, et al. Current standards of care in glioblastoma therapy. In: De Vleeschouwer S, editor., et al., Glioblastoma [internet]. Brisbane: Codon Publications; 2017.

    Chapter  Google Scholar 

  5. Kocher M, Ruge MI, Galldiks N, Lohmann P. Applications of radiomics and machine learning for radiotherapy of malignant brain tumors. Strahlenther Onkol. 2020;196(10):856–67.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Galldiks N, Dunkl V, Ceccon G, Tscherpel C, Stoffels G, Law I, et al. Early treatment response evaluation using FET PET compared to MRI in glioblastoma patients at first progression treated with bevacizumab plus lomustine. Eur J Nucl Med Mol Imaging. 2018;45(13):2377–86.

    Article  CAS  PubMed  Google Scholar 

  7. Pyka T, Gempt J, Hiob D, Ringel F, Schlegel J, Bette S, et al. Textural analysis of pre-therapeutic [18F]-FET-PET and its correlation with tumor grade and patient survival in high-grade gliomas. Eur J Nucl Med Mol Imaging. 2016;43(1):133–41.

    Article  PubMed  Google Scholar 

  8. Carles M, Popp I, Starke MM, Mix M, Urbach H, Schimek-Jasch T, Eckert F, Niyazi M, Baltas D, Grosu AL. FET-PET radiomics in recurrent glioblastoma: prognostic value for outcome after re-irradiation? Radiat Oncol. 2021;16(1):46.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Suchorska B, Jansen NL, Linn J, Kretzschmar H, Janssen H, Eigenbrod S, Simon M, Pöpperl G, Kreth FW, la Fougere C, Weller M, Tonn JC, Network GG. Biological tumor volume in 18FET-PET before radiochemotherapy correlates with survival in GBM. Neurology. 2015;84(7):710–9.

    Article  CAS  PubMed  Google Scholar 

  10. Wirsching HG, Roelcke U, Weller J, Hundsberger T, Hottinger AF, von Moos R, Caparrotti F, Conen K, Remonda L, Roth P, Ochsenbein A, Tabatabai G, Weller M. MRI and 18FET-PET predict survival benefit from bevacizumab plus radiotherapy in patients with isocitrate dehydrogenase wild-type glioblastoma: results from the randomized ARTE trial. Clin Cancer Res. 2021;27(1):179–88.

    Article  CAS  PubMed  Google Scholar 

  11. Leijenaar RTH, Nalbantov G, Carvalho S, Van Elmpt WJC, Troost EGC, Boellaard R, et al. The effect of SUV discretization in quantitative FDG-PET radiomics: the need for standardized methodology in tumor texture analysis. Sci Rep. 2015;5(5):11075.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Altazi BA, Zhang GG, Fernandez DC, Montejo ME, Hunt D, Werner J, et al. Reproducibility of F18-FDG PET radiomic features for different cervical tumor segmentation methods, gray-level discretization, and reconstruction algorithms. J Appl Clin Med Phys. 2017;18(6):32–48.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Lee M, Woo B, Kuo MD, Jamshidi N, Kim JH. Quality of radiomic features in glioblastoma multiforme: impact of semi-automated tumor segmentation software. Korean J Radiol. 2017;18(3):498–509.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Liu Z, Wang S, Dong D, Wei J, Fang C, Zhou X, et al. The applications of radiomics in precision diagnosis and treatment of oncology: opportunities and challenges. Theranostics. 2019;9(5):1303–22.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Rizzo S, Botta F, Raimondi S, Origgi D, Fanciullo C, Morganti AG, et al. Radiomics: the facts and the challenges of image analysis. Eur Radiol Exp. 2018;2(1):36.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Unterrainer M, Vettermann F, Brendel M, Holzgreve A, Lifschitz M, Zähringer M, et al. Towards standardization of 18F-FET PET imaging: do we need a consistent method of background activity assessment? EJNMMI Res. 2017;7:1–8.

    Article  CAS  Google Scholar 

  17. Vomacka L, Unterrainer M, Holzgreve A, Mille E, Gosewisch A, Brosch J, et al. Voxel-wise analysis of dynamic 18F-FET PET: a novel approach for non-invasive glioma characterisation. EJNMMI Res. 2018;8(1):91.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Niyazi M, Brada M, Chalmers AJ, Combs SE, Erridge SC, Fiorentino A, Grosu AL, Lagerwaard FJ, Minniti G, Mirimanoff RO, Ricardi U, Short SC, Weber DC, Belka C. ESTRO-ACROP guideline “target delineation of glioblastomas.” Radiother Oncol. 2016;118(1):35–42.

    Article  PubMed  Google Scholar 

  19. Fleischmann DF, Unterrainer M, Schön R, Corradini S, Maihöfer C, Bartenstein P, Belka C, Albert NL, Niyazi M. Margin reduction in radiotherapy for glioblastoma through 18F-fluoroethyltyrosine PET? A recurrence pattern analysis. Radiother Oncol. 2020;145:49–55.

    Article  CAS  PubMed  Google Scholar 

  20. Branchini M, Zorz A, Zucchetta P, Bettinelli A, De Monte F, Cecchin D, et al. Impact of acquisition count statistics reduction and SUV discretization on PET radiomic features in pediatric 18F-FDG-PET/MRI examinations. Phys Med. 2019;59:117–26.

    Article  PubMed  Google Scholar 

  21. Tsujikawa T, Tsuyoshi H, Kanno M, Yamada S, Kobayashi M, Narita N, et al. Selected PET radiomic features remain the same. Oncotarget. 2018;9(29):20734–46.

    Article  PubMed  PubMed Central  Google Scholar 

  22. PyRadiomics Radiomic Features.

  23. van Griethuysen J, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, Beets-Tan R, Fillion-Robin JC, Pieper S, Aerts H. Computational radiomics system to decode the radiographic phenotype. Can Res. 2017;77(21):e104–7.

    Article  Google Scholar 

  24. Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, et al. Scikit-learn: machine learning in Python. J Mach Learn Res. 2011;12(Oct):2825–30.

    Google Scholar 

  25. Geron A. Dimensionality reduction. In: Géron A, editor. Hands-on machine learning with scikit-learn and TensorFlow. 1st ed. Beijing: O’Reilly Media; 2017. p. 207–27.

    Google Scholar 

  26. sklearn.decomposition.PCA | scikit-learn.

  27. Davidson-Pilon C, Kalderstam J, Jacobson N, Zivich P, Kuhn B, Williamson M, et al. CamDavidsonPilon/lifelines: v0.24.4. 2020.

  28. Fratello M, Tagliaferri R. Decision trees and random forests. In: Ranganathan S, Gribskov M, Nakai K, Schoenbach C, editors. Encyclopedia of bioinformatics and computational biology, vol. 1–3. Elsevier: Amsterdam; 2019. p. 374–83.

    Chapter  Google Scholar 

  29. Ensemble methods | scikit-learn 0.22.2 documentation.

  30. EORTC-NCIC 26981/22981.

  31. Herrlinger U, Tzaridis T, Mack F, Steinbach JP, Schlegel U, Sabel M, Hau P, Kortmann RD, Krex D, Grauer O, Goldbrunner R, Schnell O, Bähr O, Uhl M, Seidel C, Tabatabai G, Kowalski T, Ringel F, Schmidt-Graf F, Suchorska B, Neurooncology Working Group of the German Cancer Society. Lomustine-temozolomide combination therapy versus standard temozolomide therapy in patients with newly diagnosed glioblastoma with methylated MGMT promoter (CeTeG/NOA-09): a randomised, open-label, phase 3 trial. Lancet (London, England). 2019;393(10172):678–88.

    Article  CAS  PubMed  Google Scholar 

  32. Ramia-de-Cap M, Kaul S. Limitations of multivariate survival analysis. Mod Pathol. 2021.

    Article  PubMed  Google Scholar 

  33. Chaddad A, Sabri S, Niazi T, et al. Prediction of survival with multi-scale radiomic analysis in glioblastoma patients. Med Biol Eng Comput. 2018;56:2287–300.

    Article  PubMed  Google Scholar 

  34. Vils A, Bogowicz M, Tanadini-Lang S, Vuong D, Saltybaeva N, Kraft J, Wirsching HG, Gramatzki D, Wick W, Rushing E, Reifenberger G, Guckenberger M, Weller M, Andratschke N. Radiomic analysis to predict outcome in recurrent glioblastoma based on multi-center MR imaging from the prospective DIRECTOR trial. Front Oncol. 2021;14(11):636672.

    Article  Google Scholar 

  35. Lao J, Chen Y, Li ZC, Li Q, Zhang J, Liu J, Zhai G. A deep learning-based radiomics model for prediction of survival in glioblastoma multiforme. Sci Rep. 2017;7(1):10353.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Bae S, Choi YS, Ahn SS, Chang JH, Kang SG, Kim EH, Kim SH, Lee SK. Radiomic MRI phenotyping of glioblastoma: improving survival prediction. Radiology. 2018;289(3):797–806.

    Article  PubMed  Google Scholar 

  37. Wang B, Zhang S, Wu X, Li Y, Yan Y, Liu L, Xiang J, Li D, Yan T. Multiple survival outcome prediction of glioblastoma patients based on multiparametric MRI. Front Oncol. 2021;25(11):778627.

    Article  Google Scholar 

  38. Lombardi G, Spimpolo A, Berti S, Campi C, Anglani MG, Simeone R, Evangelista L, Causin F, Zorzi G, Gorgoni G, Caccese M, Padovan M, Zagonel V, Cecchin D. PET/MR in recurrent glioblastoma patients treated with regorafenib: [18F]FET and DWI-ADC for response assessment and survival prediction. Br J Radiol. 2022;95(1129):20211018.

    Article  PubMed  Google Scholar 

Download references




Open Access funding enabled and organized by Projekt DEAL.

Author information

Authors and Affiliations



TW, DF, MN and MR established the study design. MN, SC and CB were involved the preparation of the radiation treatment plans. TW and DF reviewed the radiation treatment plans included into the study. DF collected and reviewed the clinical data. AH, MU, NA and PB reviewed the 18F-FET PET images included into the study. TW and LK extracted the radiomic features. TW and MR prepared the analysis for the identification of the most relevant features for survival prediction. GL, MI and KP provided support regarding the physical-technical implementation of the analysis. TW, DF, MN and MR prepared the manuscript. Manuscript review and discussion was provided by SC, IP and AG. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Marco Riboldi.

Ethics declarations

Ethics approval and consent to participate

The retrospective analysis was approved by the local ethics committee of the LMU Munich (Study Number 20-255).

Consent for publication

Not applicable.

Competing interests

The Department of Radiation Oncology of the University Hospital of LMU Munich has research agreements with Brainlab, Elekta and Viewray.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wiltgen, T., Fleischmann, D.F., Kaiser, L. et al. 18F-FET PET radiomics-based survival prediction in glioblastoma patients receiving radio(chemo)therapy. Radiat Oncol 17, 198 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Quantitative image analysis
  • Radiomics
  • Survival analysis
  • Glioblastoma
  • Radiotherapy