Drawing the Line: Artificial Intelligence-Based Glioblastoma Segmentation and Machine Learning Assessment—An Approach to Tumor Volume Quantification and Survival Prediction

1. Introduction

Brain tumors encompass a wide range of abnormal intracranial cellular growth patterns and include some of the most aggressive types of cancer. Glioblastoma (GB) remains the most common primary brain tumor in adults, characterized by its malignant and fatal nature, with an incidence of 2–3 cases per 100,000 people [1,2,3].

The latest Classification of Tumors of the Central Nervous System (CNS) published by the World Health Organization (WHO) in 2021 establishes a histopathologically distinct diagnosis for GB. Currently, GB is identified as a diffuse astrocytic tumor in adults that is IDH-wildtype [4,5,6].

Evidence suggests that GB is predominantly prevalent in the elderly population, with a median diagnostic age of 64 years. Moreover, it is most common in male patients and it is characterized by increased cell division, invasiveness, neoangiogenesis, and central necrosis [7,8,9]. Consequently, GB has an extremely poor prognosis, being an clinical entity with a median survival time of 15-16 months and a 5-year survival rate of approximately 5-10% [10]. One of the primary factors contributing to the incurability of GB is its high rate of recurrence, which is likely due to the unfeasibility of a complete surgical removal of the tumor mass [11].

Nowadays, recent reports have analyzed comprehensive genomic and epigenomic characteristics to understand long-term survival in patients with GB [12,13]. Long-term survivors (LTS) of GB are defined as those who live longer than 2 years post-diagnosis [14]. It is generally accepted that there is a subset of GB patients who live longer than 2 years and are classified as LTS, with 5% to 13% surviving an exceptional 5 years [15]. The rarity of LTS is indicated by the sharp decline in survival rates over time, with only about 6% of patients surviving up to 5 years [16]. Thus, despite the aggressive nature of GB, there is a subset of patients who achieve long-term survival.

Despite significant advancements in defining molecular subgroups, the etiology of GB remains unknown. Less than 20% of patients exhibit a familial predisposition to cancer. Specific risk factors, such as being exposed to radiation, having a genetic predisposition, and having specific inherited disorders, may contribute to the development of the disease [17].

A hallmark of its biological aggression is that GB develops by infiltrating the apparently normal tissue surrounding the tumor, namely the peritumoral brain zone [11,18]. The ability of GB to invade surrounding brain tissue and expand quickly are its most distinctive characteristics [19]. However, reflecting the underlying mass effect and increased intracranial pressure, symptoms frequently include headaches, convulsions, altered cognition, and neurological impairments [20]. The majority of the patients with GB are diagnosed by contrast-enhanced magnetic resonance imaging (MRI) and subsequently confirmed by histological diagnosis, but less than 60% of patients over the age of 70 receive the pathological confirmation [21].

According to current neuro-oncological guidelines, the standard treatment for this tumor type is maximal safe resection, followed by radiotherapy and concomitant chemotherapy with temozolomide, as well as adjuvant treatment consisting of 6 cycles of temozolomide. Among the patients expected not to tolerate standard radiotherapy, hypofractionated radiotherapy may be used [22]. However, the current conventional therapeutic modalities offer only a limited amelioration in the clinical outcomes and unfortunately recurrence is predicted in 90% of patients within 12 months following the diagnosis [23,24,25].

Therefore, in clinical practice, early evaluation of GB requires neuroimaging analysis obtained through contrast-enhanced MRI. The classic MRI finding of GB is an irregularly shaped, rim-enhancing lesion surrounding a central hypointense necrotic core. Involvement of the deep white matter and the corpus callosum is common. Depending on the tumor volume and perilesional edema, the intra-axial expansive process can displace midline structures and present a significant mass effect on normal brain matter (Chen et al., 2022), (Alexander et al., 2017) [17,26].

In the modern era of neuroimaging, studies in the field place emphasis on the importance of an early and accurate GB imaging diagnosis, as well as on the detailed analysis of imaging pattern, in order to improve the prognosis paradigm and permit an efficient patient management [27]. Therefore, GB imaging holds the ability to contribute with precise and significant data for a highly needed progress in the treatment alternatives for these patients as well as for the monitoring of the tumor progression [28].

However, craniocerebral contrast-enhanced MRI remains the pivotal paraclinical investigation in GB, providing essential parameters which can be investigated in relation to overall survival (OS). These include the tumor volume itself, the volume of the contrast-enhancing and non-enhancing tumor, the volume of the perilesional edema, the fractal parameters of necrosis, and geometric measures of the tumoral rim. Consequently, the volumetric parameters play a crucial role in predicting outcomes for GB patients, aiding in the identification of high-risk patients requiring treatment adjustments and further seeking alternative therapies for better results [29,30].

For this reason, pre-treatment data obtained via MRI has the potential to offer significant insights into understanding the tumor kinetics [31]. Moreover, artificial intelligence (AI) is a dynamically advancing field, with exciting results based on MRI images. Numerous AI models have been the subject of current research and show the possibility of being integrated into clinical practice [32]. Thus, the application of deep learning (DL) and machine learning (ML) algorithms for MRI-based tumor segmentation and grading has yielded state-of-the-art results. These models exhibit exceptional performance, serving as critical tools for enhancing surgical strategies and augmenting therapeutic response [33].

Within this framework, our study was designed to quantify the potential correlations between the specific volumetric components of GB — including the solid tumor volume, the volume of the tumor necrosis, and the volume of the perilesional edema — as well as the OS rates in patients diagnosed with this highly aggressive cerebral pathology.

This approach seeks to evaluate the extent to which the tumor’s internal architecture and its micro- or macro-environmental impact —assessed via AI- imaging techniques—can serve as reliable predictive or prognostic markers.

Also, the goal of this research is to employ a multi-model approach, training multiple ML models to predict survival outcomes in patients with GB.

2. Materials and Methods

A retrospective analysis was conducted on patients with GB who received treatment at the “Prof. Dr. N. Oblu” Emergency Clinical Hospital (Iași, Romania) from 2019 to 2023. The final study cohort was selected based on well-defined inclusion and exclusion criteria, which are detailed below.

Inclusion criteria:

- patients with a confirmed diagnosis of GB;

- availability of high-quality, pre-treatment MRI scans;

- patients who underwent standard-of-care treatment (e.g., maximal safe resection followed by radiotherapy and/or chemotherapy);

- complete clinical records and documented survival data (OS), with clear dates of diagnosis and death or last known follow-up;

- adult patients (aged > 18 years) at the time of diagnosis.

Exclusion criteria:

- patients who received any form of anti-tumor therapy prior to the baseline MRI scan;

- poor MRI quality that precludes reliable volumetric analysis;

- patients lost to follow-up or those for whom survival data could not be accurately quantified;

- presence of other active primary cancers;

- patients who died within 30 days of surgery from non-oncological complications.

2.1. Data Collection, Data Cleaning and Preprocessing

Clinical and radiological data were systematically compiled into a dedicated database. Collected variables included demographics (age at diagnosis and sex) and temporal data (dates of diagnosis and death). Radiographic parameters encompassed the topographic localization, the manual volumetric measurements, and the automated outputs from the AI model, including the total tumor volume and the sub-compartmental volumes (contrast-enhancing component, peritumoral edema, and necrosis).

All the collected parameters were systematically organized and stored in a digital database for subsequent statistical analysis. Columns were standardized and renamed for consistency prior to analysis.

The dates of the histopathological diagnosis and the times of death were used to determine the expected survival time and the date April 1^st 2025 was used as a placeholder for missing death dates. Duplicate patient identifiers were removed to ensure patient-level analyses. All analyses were performed on the cleaned dataset, by specific exclusion reasons (missing or inconsistent dates, incomplete imaging measures).

The primary outcome for survival analyses was OS (months) measured from histopathologic diagnosis to death. Survival durations are expressed in months and an event indicator was set to 1 for observed deaths and 0 for censored observations.

We derived multiple imaging biomarkers from the raw variables. Manual tumor volumes were normalized using the scalar π/6 to convert raw measurements into a standardized volume estimate. AI-derived measurements included the total estimated volume and the compartmental volumes (contrast, necrosis, edema). Ratios and differences were also computed. To reduce extreme influences, ratio variables were truncated at predefined caps.

2.2. Manual Volumetric Assessment of the Tumor

The tumor volume was quantified using the ellipsoidal method, based on the principle that the morphology of brain pathology is most accurately approximated by an elliptical geometry rather than a cubic one [34].

The protocol involved a two-stage process. First, the 3D T1 weighted contrast-enhanced (CE-T1W) acquisitions underwent rigorous spatial normalization. This entailed aligning the medio-sagittal line in the coronal and axial planes and the bi-commissural axis in the sagittal plane to ensure anatomical consistency. Secondly, the maximum orthogonal diameters were measured across all three planes. The product of these measurements was subsequently adjusted by the ellipsoidal constant of 0.52 to derive the final volumetric data.

In order to guarantee the robustness of the findings and minimize subjective bias, all measurements were performed independently by two senior neuroradiologists.

2.3. Automated Volumetric Segmentation of the Tumor

Automated volumetric analysis was performed using the mediaire software solutions mdbrain (distributed in Romania via Supermedical), a class II CE-certified medical device. The software employs a DL algorithm based on a 3D convolutional neural network with a U-Net architecture.

To facilitate a comprehensive segmentation, CE-T1W and T2-weighted Fluid-Attenuated Inversion Recovery (T2W FLAIR) sequences were co-registered and integrated.

This multi-sequence approach allows for the simultaneous characterization of the solid enhancing component and the central necrotic core (identifiable via CE-T1W), alongside the peripheral vasogenic edema (characterized by hyperintense signal on FLAIR without contrast enhancement).

The total tumoral volume was thus partitioned into these three distinct sub-compartments, with the resulting Digital Imaging and Communication in Medicine (DICOM)-processed outputs illustrated in Figure 1.

2.4. Statistical Analysis

The distributions of continuous variables were inspected using histograms and kernel density estimates. Pairwise scatterplots and correlation assessments were used to identify strong associations between imaging-derived features (for example, AI edema vs AI necrosis vs volume measures) and with survival. Summary statistics (mean, median, standard deviation, minimum, maximum and counts) were computed overall and stratified by key categorical variables (tumor side, precise localization, gender).

We compared continuous outcomes between groups using parametric or nonparametric tests depending on the number of groups and distributional properties. For two-group comparisons we applied the Student’s t-test or Mann–Whitney U test; for series of more than two groups we used one-way ANOVA or Kruskal–Wallis tests. The value of p ≤ 0.05 was considered statistically significant.

2.5. Survival Analysis Methods

The Kaplan–Meier estimator was used to describe survival functions for the overall cohort and for the stratified groups [35]. Survival curves and median survivals with 95% confidence intervals were computed and plotted. Log-rank tests were used for unadjusted comparisons where appropriate [36]. To assess the proportional hazards assumption graphically, complementary log–log (cloglog) plots of log(-log(S(t))) versus log(time) were generated for major strata and visually inspected.

Multivariable Cox proportional hazards models were fit to estimate hazard ratios (HR) associated with demographic and imaging-derived covariates [37]. Categorical predictors were dummy-encoded (drop-first), and models were fit only when the sample size for the covariate grouping exceeded a minimum threshold to avoid unstable estimates. Duration and event variables were constructed as described above; hazard ratios with 95% confidence intervals and p-values are reported.

2.6. Predictive Modelling

The integration of AI into medical research has revolutionized the field of oncology. In the context of GB, by leveraging high-dimensional imaging data and clinical information, ML algorithms can extract valuable insights into patient prognosis. In this study, we harnessed the power of AI to explore the potential correlations between tumor volumetry and survival outcomes in GB patients, laying the groundwork for the development of more accurate predictive models that can inform clinical decision-making.

2.6.1. Regression – Survival in Months

We modeled the problem as a regression estimator. For handling the large number of features compared to the total patients in the study, we chose random forest (RF) models.

For creating the RF model, the dataset is split into sub-datasets. Each sub-dataset is used to construct a decision tree (DT), in the end forming a forest of many trees. The strength of RF lies in its assembly of multiple DTs, enabling it to achieve higher and more robust performance by better generalizing across the dataset.

To predict continuous survival (months) we compared four model families: RF Regressor, XGBoost Regressor, k-Nearest Neighbors (KNN) Regressor, and a fully connected Neural Network implemented in Keras [38,39,40,41].

The feature set used for model development is: gender, age at diagnosis, manual computed volume normalised, AI model estimated volume, AI model contrast volume, AI model necrosis volume, AI model edema volume, all volume ratios, tumor location, and tumor side.

A cross-validation holdout was used with a 20% testing and 80% training split. Model performance was evaluated on the held-out test set using mean absolute error (MAE), mean squared error (MSE) and${\mathrm{R}}^2$.

2.6.2. Classification – 6-Month Threshold

For binary prediction of survival ≥6 months we compared Random Forest Classifier, XGBoost Classifier, k-Nearest Neighbors Classifier, and a Neural Network trained with class weights. The same train/test split strategy described above was used. The class ratio split was used in training to adapt to the imbalance. Classification performance metrics included accuracy, precision, recall, F1 score and area under the ROC curve (AUC). Confusion matrices were reported to illustrate the trade-off between sensitivity and specificity.

2.7. Model Explainability

Model-agnostic local explanations were produced using LIME (local interpretable model-agnostic explanations) for tabular data. An explainer was initialized on the training set and explanations were generated for a sample of test instances; feature contributions were aggregated across instances to produce mean absolute importances per model.

2.8. Software, packages and reproducibility

Analyses were conducted in Python 3.13 using pandas, NumPy, SciPy, scikit-learn, xgboost, Keras/TensorFlow, lifelines, seaborn, matplotlib, lime and TableOne. Random seeds were set where models support deterministic seeds.

3. Results

The preliminary screening identified 695 patients with a primary diagnosis of GB treated at the “Prof. Dr. N. Oblu” Emergency Clinical Hospital between January 2019 and December 2023. A significant number of cases were ineligible for inclusion in the study due to missing pre-operative imaging data. More specifically, the initial MRI investigations were conducted externally and could not be accessed for volumetric processing. Furthermore, numerous patients were excluded from the study because their MRI scans were unsuitable for AI processing owing to the presence of artifacts.

A total of 79 patients were ultimately included in the study after meeting all the eligibility criteria. The baseline characteristics of the study cohort are summarized in Table 1.

SD: standard deviation; N: number; N/A: not applicable.

3.1. Group Comparisons

Histograms of manual normalized tumor volume and AI-estimated volume are shown in Figure 2, demonstrating right-skewness and overlapping.

Pairwise scatterplots (Figure 3) showed a strong positive correlation between the AI-estimated volume and the manual normalized volume (Spearman ρ = 0.93, p value = 2.57e-38) and low correlation between the AI edema and the AI necrosis (p value = 0.26). These relationships motivated inclusion of both volume and compartment ratios in downstream models.

The survival outcomes were stratified by tumor laterality (left, center, right), as detailed in Table 2. The statistical tests indicated a significant difference in median survival between tumor sides (ANOVA p value = 0.0022), central tumors being the most aggressive.

Max: maximum value; Min: minimum value; OS: overall survival; SD: standard deviation.

The comparison of edema values derived from the AI based model according to tumor location (left/central/right) demonstrated a statistically significant difference among groups (one-way ANOVA, p = 0.049), with centrally located tumors showing lower AI-derived edema values compared with laterally located tumors (Figure 4).

3.2. Survival Analysis

The Kaplan–Meier survival curves for the overall study population stratified by sex are presented in Figure 5, while the ones stratified according to the precise tumor localization are shown in Figure 6. The median overall survival for the cohort was 8.16 months (95%, confidence interval 0.44–0.62).

A multivariable Cox proportional hazards analysis was performed, adjusting for age at diagnosis, manual normalized volume, and AI-estimated volumes. Significant independent predictors associated with decreased survival rates included higher manual normalized volume (HR = 0.01051, p = 0.01906), higher AI model contrast volume (HR = 0.01794, p = 0.01106), and higher AI model estimated volume (HR = 0.01052, p = 0.03193).

3.3. Regression Models Predicting Continuous Survival Performance

The four trained regression models, as described in the Materials and Methods section, were used to predict continuous survival outcomes expressed in months. Performance on the held-out test set is summarized in Table 3. The XGBoost regressor produced MAE (mean absolute error) = 5.342 months, MSE (mean squared error) = 36.268; The Random Forest and KNN models demonstrated slightly inferior predictive performance, whereas the Neural Network exhibited evidence of overfitting, as reflected by suboptimal MSE performance.

KNN: k-Nearest Neighbors; MAE: mean absolute error; MSE: mean squared error.

The classification model, a simplified approach to the regression model, is less sensitive to large errors in exact month estimates, which degrades regression performance. In practice, this simplification improved predictive performance and achieved higher discrimination (AUC) than the regression models' continuous-error metrics. We therefore present the classification results alongside regression; however, this approach misses information regarding finer-grained survival differences, which require a prespecified threshold (6 months here) chosen for its clinical relevance.

3.4. Classification Models (6-Months Survival) Performance

For binary prediction of survival ≥6 months we compared Random Forest Classifier, XGBoost Classifier, k-Nearest Neighbors Classifier, and a Neural Network trained with class weights. Classification performance metrics included accuracy, precision, recall, F1 score and area under the ROC curve (AUC). Confusion matrices with absolute counts and row percentages were reported to illustrate the trade-off between sensitivity and specificity (see Figure 7 and Figure 8).

3.5. Feature Importance and Explainability

Feature importances indicate that the top predictors included AI-model necrosis, manually normalized volume, and AI edema volume. Aggregated LIME explanations corroborated the importance of these features across multiple model families. The most relevant features for each model are also consistent with our intuition, as shown in Table 4.

AI: artificial intelligence; KNN: k-Nearest Neighbors; LIME: local interpretable model-agnostic explanations.

4. Discussion

GB constitutes the most prevalent primary brain malignancy in adults, distinguished by hallmark features such as extensive necrosis and diffuse infiltration into perilesional parenchyma. These pathological attributes validate its aggressive clinical profile, marked by pronounced resistance to standard-of-care therapies and a high recurrence rate. Consequently, the prognosis remains dismal, underscoring the imperative for expedited diagnosis and tailored therapeutic interventions to optimize patient management [29,42,43].

Nowadays, progress is being made in identifying prognostic factors derived from MRI features of GB, promising a targeted, personalized approach to achieve positive long-term outcomes for GB patients. In this regard, Sanghani and colleagues [44] evaluated 13 tumor shape features for predicting OS in GB patients. Significant 3D shape features, including the bounding ellipsoid volume ratio, sphericity, and spherical disproportion, derived from both FLAIR and contrast-enhanced T1-weighted masks, were prognostic indicators of OS in GB patients. These features were more effective when derived from the contrast-enhanced T1-weighted sequence compared to the FLAIR sequence. Similar to our study, their research used univariate and multivariate Cox regression analyses, along with Kaplan-Meier survival curves, in order to assess the prognostic value of these shape features for OS in GB patients. Their findings suggest that irregular tumor shapes are associated with poor survival, whereas regular tumor shapes are associated with prolonged survival.

However, the quantification of tumor volume as a prognostic marker for GB was reported by Reeves and Marks in 1979 [45]. Their research indicated that larger GB are associated with shorter post-treatment survival. Even today, the volume of GB measured by MRI is considered one of the most significant prognostic factors in consequence of its strong correlation with disease progression and patient survival. Our findings reveal that larger tumors are usually associated with a poorer prognosis, as they indicate a higher tumor burden and a greater risk of impairing essential brain functions [46,47].

Regarding tumor volume, a recent study [48] focuses on the preoperative MRI evaluation and provides insights into potential prognostic indicators and the importance of accurate GB volume measurements for predicting patient outcomes. The cohort consisted of 44 patients, predominantly male, with a median age of 58.15 years. Palpan Flores et al. used various volumetric methods in order to measure tumor size and compartments, showing that GB volumes were optimal neuroimaging biomarkers for predicting survival in GB patients. It is remarkable that OS was associated with age ≥ 65 years and a tumor volume measured in FLAIR-T2 of ≥ 2,000 mm² or ≥ 60 cm³.

Moreover, understanding the dynamics of contrast-enhancing tumor volume over time can provide essential awareness of the aggressiveness of the disease and guide personalized treatment strategies for patients with GB.

In this regard, Auer and colleagues [49] proposed turning the spotlight on quantitative volumetric analysis, recognizing the limitations of the existing methods and the need for more robust prognostic tools. Their research included a cohort of 49 patients diagnosed with GB, with the objective of examining the possibility of 3D volumetric measurements of tumor enhancement as a predictive radiomic imaging biomarker. Results indicate a substantial correlation between high contrast-enhanced tumor volume and shortened survival, providing a baseline for considering 3D volumetric approaches as predictive biomarkers in GB. Thus, this study provides a glimpse into the prospective resolution of the ongoing discussion on the effectiveness of different measurement approaches, illustrating the potential of 3D volumetric quantification as a clinically relevant tool for GB prognosis.

The GB volume at diagnosis provides prognostically relevant information and supports clinicians in making appropriate, individualized management decisions. In addition to tumor volume, the presence and extent of necrosis—an imaging hallmark of GB — also influence prognosis and therapeutic planning and therefore warrant particular attention.

Necrosis is a diagnostic hallmark of GB, as its presence is identified in over 85% of the cases [50]. Researchers have long been intrigued by the relationship between necrosis and GB progression because understanding this connection could lead to advances in treatment options for patients. Studies have shown that necrosis within GB tumors can impact patient survival rates as well as their response to therapy, highlighting the importance of further investigating this phenomenon [51,52].

Unquestionably, necrotic patterns in GB that are linked to patient survival have been the focus of extensive research in recent years. Understanding the relationship between necrotic areas within GB tumors and patient outcomes is vitally important for developing more appropriate treatment options. In this regard, Liu and colleagues [53] investigated the prognostic role of necrotic patterns in GB using fractal dimension (FD) and lacunarity analyses derived from MRI data. The authors analyzed clinical and MRI data of 95 GB patients and calculated the FD and lacunarity of the necrosis by fractal analysis. The study revealed significant correlations between lower FD values and larger lacunarity values and a decline in PFS (p=0.006 and p=0.012, respectively) and OS (p=0.008 and p=0.005, respectively) in patients with GB. This demonstrates that the fractal parameters of necrosis in GB can predict patient survival and are accordingly linked to the biological processes concerning tumor necrosis. In addition, lacunarity was found to be associated with the suppression of apoptosis – and necrosis-related biological processes, shedding light on the importance of necrosis in GB prognosis and the molecular pathways that lead to it.

As a further matter, Curtin and collaborators [54] explored the use of morphological metrics, specifically lacunarity and FD, for the analysis of GB abnormalities on MRI and their connection to patient survival. The study found significant correlations between these morphological metrics and patient outcomes, with the shape of T2/FLAIR (fluid-attenuated inversion recovery) abnormalities showing the strongest link to OS.

Most recently, a study published in 2024 by Ma et al. [55], involving a cohort of 150 patients, demonstrated that imaging-defined necrosis was strongly associated with glioblastoma and with specific genetic alterations, including 1p/19q non-codeletion and homozygous deletion of CDKN2A/B, and functioned as an independent prognostic marker for tumor outcome. This research emphasizes the critical role of necrosis, particularly imaging necrosis, in evaluating GB, providing valuable insights for an accurate diagnosis and prognosis. Additionally, in this study, dynamic contrast-enhanced-MRI metrics demonstrated a high diagnostic efficiency in identifying tumor necrosis.

However, some studies discuss the anatomical features of GB on imaging, including cysts, necrotic changes, and the extent of edema, and emphasize the importance of complete resection for better outcomes. Other studies focus on the prognostic value of morphological MRI-based imaging biomarkers, such as tumor volume and surface regularity, for predicting survival in GB patients. The common direction is to highlight the significance of specific imaging characteristics in understanding and predicting outcomes for patients with GB. For example, in recent years, numerous studies have focused on the link between perilesional edema and the prognosis in GB patients (Table 5).

CE-T1WI: contrast-enhanced T1-weightedimage, FLAIR: fluid attenuated inversion recovery, GB: glioblastoma, MRI: magnetic resonance imaging, OS: overall survival, PFS: progression-free survival, PTBE: peritumoral brain edema, T1WI: T1 weighted image, T2WI: T2 weighted image.

Overall, by incorporating tumor volumes into ML models for prognostic evaluation, our study revealed that the Neural Network achieved superior discriminative performance, while the Random Forest displayed robust clinical relevance in confusion matrix assessments. In particular, the Random Forest achieved a favorable balance between true-positive and true-negative classifications in the held-out dataset, suggesting significant utility in clinical contexts where specific error profiles are critical, such as minimizing the risk of misclassifying high-risk patients.

From a practical perspective, the Neural Network may be considered the primary model for survival ranking and prediction, while the Random Forest represents a complementary alternative when greater interpretability of decision thresholds or increased robustness to noise in tabular features is desired.

The findings presented in our research are likely to make a significant impact on the diagnosis and management of GB patients. However, our study also has some limitations. This is a retrospective analysis and has a restricted number of subjects in the cohort provided. Furthermore, the unicentre nature of this investigation has a detrimental impact on its validity. Therefore, this study requires further confirmation through multicentred and prospective research in order to be applicable, and so that its potential findings may be effectively used in the current clinical practice.

5. Conclusions

As a summary, our study contributes to the development of advanced computational tools for neuro-oncology by integrating AI-based segmentation and multi-model ML approaches, and the results highlight the importance of imaging biomarkers in understanding the prognosis of GB.

Future work may reasonably explore model combination strategies, such as ensemble or stacked approaches, or selection of the optimal model based on predefined clinical priorities (e.g., optimization of sensitivity versus specificity), prior to external validation.