# Radiomics-Based Classification of Tumor and Healthy Liver on Computed Tomography Images

Vincent Béni Sèna Zossou, Freddy Houéhanou Rodrigue Gnangnon, Olivier

Biaou, Florent de Vathaire, Rodrigue S. Allodji, Eugène C. Ezin

► To cite this version:

Vincent Béni Sèna Zossou, Freddy Houéhanou Rodrigue Gnangnon, Olivier Biaou, Florent de Vathaire, Rodrigue S. Allodji, et al. Radiomics-Based Classification of Tumor and Healthy Liver on Computed Tomography Images. Cancers, 2024, 16 (6), 10.3390/cancers16061158. [https://doi.org/10.3390/cancers16061158](https://doi.org/10.3390/cancers16061158)

HAL Id: hal-04831583

https://hal.science/hal-04831583v1

Submitted on 16 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Radiomics-Based Classification of Tumor and Healthy Liver on Computed Tomography Images

Vincent-Béni Sèna Zossou <fnref n="1" />,<fnref n="2" />,<fnref n="3" />,<fnref n="4" />,*, Freddy Houéhanou Rodrigue Gnangnon <fnref n="5" />, Olivier Biaou <fnref n="6" />, Florent de Vathaire <fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Rodrigue S. Allodji <fnref n="1" />,<fnref n="2" />,<fnref n="3" />, and Eugène C. Ezin <fnref n="7" />,<fnref n="8" />

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Zossou, V.B.S.; Gnangnon, F.H.R.; Biaou, O.; de Vathaire, F.; Allodji, R.S.; Ezin, E.C. Radiomics-Based Classification of Tumor and Healthy Liver on Computed Tomography Images. *Cancers* **2024**, *16*, 1158. https://doi.org/10.3390/cancers16061158

Academic Editor: Dania Cioni

Received: 15 February 2024

Revised: 10 March 2024

Accepted: 11 March 2024

Published: 14 March 2024

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Simple Summary:** Liver malignancies, particularly hepatocellular carcinoma, and metastases stand as prominent contributors to cancer mortality. Within abdominal computed tomography imaging, much of data remains underused by radiologists. Radiomics uses advanced image analysis to extract quantitative features from medical scans for deeper diagnosis, treatment, and prognosis insights. Machine learning algorithms enable analyzing these features, facilitating an automatic, rapid, and efficient medical management process. We used these algorithms to train models that can distinguish between healthy livers and those with tumors, as well as between malignant and benign tumors, using CT images from the electronic medical record of the Centre National Hospitalier Universitaire Hubert Koutoukou Maga (CNHU-HKM) in Benin. The high correlation scores suggest that the radiomics signature is a prognostic biomarker for hepatic tumor screening.

**Abstract:** Liver malignancies, particularly hepatocellular carcinoma and metastasis, stand as prominent contributors to cancer mortality. Much of the data from abdominal computed tomography images remain underused by radiologists. This study explores the application of machine learning in differentiating tumor tissue from healthy liver tissue using radiomics features. Preoperative contrast-enhanced images of 94 patients were used. A total of 1686 features classified as first-order, second-order, higher-order, and shape statistics were extracted from the regions of interest of each patient's imaging data. Then, the variance threshold, the selection of statistically significant variables using the Student's t-test, and lasso regression were used for feature selection. Six classifiers were used to identify tumor and non-tumor liver tissue, including random forest, support vector machines, naive Bayes, adaptive boosting, extreme gradient boosting, and logistic regression. Grid search was used as a hyperparameter tuning technique, and a 10-fold cross-validation procedure was applied. The area under the receiver operating curve (AUROC) assessed the performance. The AUROC scores varied from 0.5929 to 0.9268, with naive Bayes achieving the best score. The radiomics features extracted were classified with a good score, and the radiomics signature enabled a prognostic biomarker for hepatic tumor screening.

**Keywords:** radiomic features; classification; machine learning; liver lesions; hepatocellular carcinoma; metastasis

<footnote n="1">¹ Université Paris-Saclay, UVSQ, Univ. Paris-Sud, CESP, Équipe Radiation Epidemiology, 94805 Villejuif, France; rodrigue.allodji@gustaveroussy.fr (R.S.A.)</footnote>

<footnote n="2">² Department of Clinical Research, Radiation Epidemiology Team, Gustave Roussy, 94805 Villejuif, France</footnote>

<footnote n="3">³ Centre de Recherche en Épidémiologie et Santé des Populations (CESP), U1018, Institut National de la Santé et de la Recherche Médicale, 94805 Villejuif, France</footnote>

<footnote n="4">⁴ Ecole Doctorale Sciences de l'Ingénieur, Université d'Abomey-Calavi, Abomey-Calavi 384, Benin</footnote>

<footnote n="5">⁵ Department of Visceral Surgery, CNHU-HKM, Cotonou 229, Benin</footnote>

<footnote n="6">⁶ Department of Radiology, CNHU-HKM, Cotonou 229, Benin</footnote>

<footnote n="7">⁷ Institut de Formation et de Recherche en Informatique, Université d'Abomey-Calavi, Abomey-Calavi 384, Benin; eugene.ezin@uac.bj</footnote>

<footnote n="8">⁸ Institut de Mathématiques et de Sciences Physiques, Université d'Abomey-Calavi, Dangbo 384, Benin</footnote>

<footnote>* Correspondence: vincent.zossou@gustaveroussy.fr</footnote>
# 1. Introduction

Liver cancer is the sixth-most common and third-deadliest cancer worldwide, with projections suggesting an increase to over a million cases by 2030. It primarily affects males, especially in regions with high hepatitis B prevalence, such as Asia and Sub-Saharan Africa. Hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma make up 99% of liver cancer cases [1]. Early diagnosis, crucial for managing liver diseases, relies significantly on diagnostic imaging to detect primary tumors or metastases and monitor tumor progression.

High-resolution imaging modalities such as magnetic resonance imaging (MRI), ultrasound, and computed tomography (CT) help detect liver tumors. Radiologists and oncologists most frequently use computed tomography (CT) for liver lesion evaluation and staging [2]. Radiomics is a non-invasive method for determining the relationship between quantitative medical picture characteristics and underlying biological events [3]. The use of radiomics-based computer-aided diagnosis (RCAD) [4] from CT images significantly enhances diagnostic processes by providing an objective and rapid assessment, which is particularly beneficial for less experienced or less confident radiologists. Additionally, it is especially relevant in settings where CT scans are a primary tool for routine checkups, making it a valuable asset for medical facilities with varying resources.

The approach to liver cancer treatment varies based on the tumor's phenotype [5,6]. Initially, radiologists employ MRI or CT scans to assess the phenotype, and various manual grading criteria for liver tumors have been established [7]. However, the visual characteristics of these lesions on scans can greatly differ due to histological variations [8], leading to subjective interpretations and a lack of consensus on their exact definitions [7]. While tissue biopsy offers definitive histological confirmation, it is not routinely required for HCC diagnosis thanks to advancements in non-invasive methods utilizing characteristic imaging features and biomarkers.

Many RCAD studies have been conducted and published. In 2007, S. G. Mougiakakou and others [9] used the gray-level histogram (GLH), the gray-level co-occurrence matrix (GLCM), and other methods to extract texture features from CT images. They then built models using support vector machine (SVM), K-nearest neighbors, the probabilistic neural network (PNN), and other methods to help with the diagnosis of normal liver, hepatic cyst, hepatic hemangioma (HEM), and hepatocellular carcinoma (HCC). In 2013, S. S. Kumar et al. [10] used the GLH, GLCM, wavelet transform, contourlet transform, and their fusions to extract texture data from CT images and then constructed PNN models for classifying HCC and HEM. J. J. Qiu et al. [11] investigated the same classification task, applying SVM models to data extracted from plain CT images with the methods of the GLH, GLCM, etc. In 2017, Y. Zhou et al. [12] used the Laplacian of Gaussian (LoG) spatial band-pass filter, GLH, and GLCM to extract textural data from arterial- and portal venous-phase CT images and, then, constructed logistic regression (LR) models to identify early recurrence of HCC within a year. In 2020 and 2021, Nie et al. [13,14] designed a nomogram based on CT radiomics to differentiate HCC from focal nodular hyperplasia (FNH) and hepatocellular adenoma (HA) in the normal non-cirrhotic liver with high accuracy. In 2022, S. Lysdahlgaard used morphological and statistical methods along with LR, random forest (RF), and SVM models to obtain radiomic features like the size, shape, and location of the tumor and healthy liver tissue [15].

Recent research has demonstrated significant potential in using radiomics to identify variables associated with clinical outcomes [16,17]. However, only a few studies have used the entire liver as the region of interest (ROI) in their research. Additionally, no African investigations have been conducted on applying radiomics or machine learning to classify liver lesions. This study evaluates the effectiveness of radiomic features extracted from CT scans of the liver in distinguishing between healthy and cancerous tissues and differentiating benign tumors from malignant ones, thereby leveraging radiomics for nuanced diagnostic insights into liver health. To our knowledge, this is the first study in Africa.
## 2. Materials and Methods

### 2.1. Proposed System

A machine learning (ML) model has been developed to accurately classify healthy and tumor livers using radiomic features extracted from CT images. There are five stages in the workflow process of the proposed system: segmentation of ROIs, features extraction, statistical analysis, feature selection, and classification. The entire liver was interactively evaluated as the ROI. Figure 1 shows the block diagram of the proposed system. The work consists of two main phases: classifying healthy and tumor livers, then, among those with tumors, classifying them as benign or malignant.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** The block diagram of the proposed system.</figcaption></figure>

### 2.2. Dataset Description

In this study, 94 CT images obtained from the electronic medical record of the Centre National Hospitalier Universitaire Hubert Koutoukou Maga (CNHU-HKM) for January 2018 to December 2021 were retrospectively included. The dataset reveals a diverse array of tumor sizes, with dimensions ranging from a minimum of 3.082 mm to a maximum of 144.4 mm and a median value of 11.18 mm. The cohort comprises 38 healthy patients, 23 confirmed to have HCC and 33 having metastatic tumors. Each scanner folder is an abdominal CT scan including four injection phases of contrast enhancement: 15–20 s for the early arterial phase, 45 s for the late arterial phase, 80 s for the portal venous phase, and 4 min for the equilibrium delay phase. Only the portal phase was used because it can show complex structural details, improve vascular clarity, boost tissue contrast, and make it easier to find lesions. The imaging dataset's details are reported in Table 1.

<table><caption>**Table 1.** Imaging dataset's details.</caption><thead><tr><td>Image Size</td><td>Slice Thickness</td><td>Pixel Spacing</td><td>Slice Number</td></tr></thead><tbody><tr><td>256 × 256</td><td>1.5 mm</td><td>0.6172 mm</td><td>1200–1700</td></tr></tbody></table>

### 2.3. ROI Segmentation

Image segmentation is the procedure by which an image is partitioned into significant regions, with distinct attributes assigned to individual pixels. The performance of image-segmentation methods is contingent upon various factors, including the type of image, manner of application, size, and color intensity. Two radiologists, each with three years of experience in abdominal CT diagnosis, used the Radiology Informatics Laboratory Contour
(RIL-Contour) software [18] to segment the ROI in the images manually. Figure 2 shows a sample CT scan of a liver tumor and the segmentation.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** The segmentation of liver. (left) Original image with HCC masses inside. (middle) Liver segmented. (right) Tumor segmented.</figcaption></figure>

## 2.4. *Intensity Normalization*

Variations in CT image intensity can occur, even when the same CT scanner and identical scanning parameters are used. These variations may affect the process of extracting features from the images. Furthermore, the primary image may be blurred, leading to potential confusion between the ROI and other organs or regions present in the image. The dataset was preprocessed to improve the contrast of the images. The Hounsfield unit (HU) is a relative quantitative measurement of radiodensity (a property of a tissue to attenuate/absorb an X-ray beam) that radiologists use to interpret CT images. The HU values typically vary between 1000 and 3000 [19]. The HU values of the image were set to the range [150, 300] in this study.

## 2.5. *Feature Extraction*

The main goal of radiomics is to create a machine learning algorithm that can categorize outcomes based on objective criteria, utilizing quantitative data from medical images. This work comprehensively analyzed 1686 radiomics features per patient, classified into four categories: first-order, second-order, higher-order, and shape features (Table 3). First-order features can accurately depict the spatial distribution of various voxel intensities, irrespective of the three-dimensional configuration [20]. The attributes of the second-order statistics can be identified by analyzing a density histogram as follows:

(i) The gray-level co-occurrence matrix (GLCM) that quantifies an image's texture by computing the frequency of pairs of pixels with given values and a defined spatial connection [21].

(ii) The gray-level run length matrix (GLRLM), a matrix that can extract texture features for texture analysis. The GLRLM technique is utilized to extract advanced statistical texture properties [22].

(iii) The gray-level size zone matrix (GLSZM) measures the number of gray-level zones in an image. A gray-level zone is a group of connected voxels with the same gray-level intensity [23].

(iv) The gray-level dependence matrix (GLDM), which measures the relationships between gray levels in an image. There is a gray-level dependency when a count of connected voxels within a distance $\delta$ of the center voxel depends on it [24].

Based on the first-order and second-order features, filter grids were applied to the images, and transformed features, namely higher-order statistical features, were obtained. Six filters were applied: exponential, square, square root, logarithm, Laplacian of Gaussian (LoG), and wavelet. Each filter was selected for its capacity to enhance or reveal specific aspects of imaging data, contributing to a more comprehensive and precise analysis. The exponential filter intensifies low-amplitude variations in the image, enabling the better distinction of subtle regions that might be significant for radiomic analysis. The
square filter amplifies high intensities, while the square root filter mitigates high intensities, making low-contrast areas more discernible. The logarithmic filter compresses the dynamic range of high-dynamic images, facilitating the visualization and analysis of fine details in bright and dark regions. The LoG is particularly effective for edge detection and contour enhancement. It accentuates regions with rapid changes in intensity, which is crucial for identifying lesion boundaries or other relevant structures. The sigma values for LoG were 0.5, 1.0, 2.0, 3.0, 4.0, and 5.0. Wavelet analysis decomposes the image into components at different scales, allowing for extracting pertinent information that other methods might miss. This is particularly useful for capturing localized and multi-scale features in imaging data. The wavelet transforms and decomposes an image into a series of sub-images, each representing different frequency components (low (L) and high (H)) along different orientations (horizontal, vertical, and diagonal). L denotes low frequency (smooth parts of the image), and H denotes high frequency (edges and fine details). Based on these frequencies, eight levels of wavelet decomposition have been obtained: wavelet-LHL, wavelet-LHH, wavelet-HLL, wavelet-LLH, wavelet-HLH, wavelet-HHH, wavelet-HHL, and wavelet-LLL). The extracted features are shown in Table 2.

<table><caption>**Table 2.** Radiomics features in the radiomics analysis.</caption><thead><tr><th>Types</th><th>Features</th></tr></thead><tbody><tr><td>Shape (n = 14)</td><td>Elongation, Flatness, Least Axis Length, Major Axis Length, Maximum 2D Diameter Column, Maximum 2D Diameter Row, Maximum 2D Diameter Slice, Maximum 3D Diameter, Mesh Volume, Minor Axis Length, Sphericity, Surface Area, Surface Volume Ratio, Voxel Volume</td></tr><tr><td>First-order statistics (n = 18)</td><td>10 Percentile, 90th Percentile, Energy, Entropy, Interquartile Range, Kurtosis, Maximum, Mean Absolute Deviation, Mean, Median, Minimum, Range, Robust Mean Absolute Deviation, Root Mean Square, Skewness, Total Energy, Uniformity, Variance</td></tr><tr><td>Second-order statistics</td><td></td></tr><tr><td>GLCM (n = 24)</td><td>Auto-Correlation, Cluster Prominence, Cluster Shade, Cluster Tendency, Contrast, Correlation, Difference Average, Difference Entropy, Difference Variance, Id, Idm, Idmn, Idn, Imc1, Imc2, Inverse Variance, Joint Average, Joint Energy, Joint Entropy, MCC, Maximum Probability, Sum Average, Sum Entropy, Sum of Squares</td></tr><tr><td>GLRLM (n = 16)</td><td>Gray-Level Non-Uniformity, Gray-Level Non-Uniformity Normalized, Gray-Level Variance, High Gray-Level Run Emphasis, Long-Run Emphasis, Long-Run High Gray-Level Emphasis, Long-Run Low Gray-Level Emphasis, Low Gray-Level Run Emphasis, Run Entropy, Run Length Non-Uniformity, Run Length Non-Uniformity Normalized, Run Percentage, Run Variance, Short-Run Emphasis, Short-Run High Gray-Level Emphasis, Short-Run Low Gray-Level Emphasis</td></tr><tr><td>GLSZM (n = 16)</td><td>Gray-Level Non-Uniformity, Gray-Level Non-Uniformity Normalized, Gray-Level Variance, High Gray-Level Zone Emphasis, Large Area Emphasis, Large Area High Gray-Level Emphasis, Large Area Low Gray-Level Emphasis, Low Gray-Level Zone Emphasis, Size Zone Non-Uniformity, Size Zone Non-Uniformity Normalized, Small Area Emphasis, Small Area High Gray-Level Emphasis, Small Area Low Gray-Level Emphasis, Zone Entropy, Zone Percentage, Zone Variance</td></tr></tbody></table>
<table><caption>**Table 2.** *Cont.*</caption><thead><tr><th>Types</th><th>Features</th></tr></thead><tbody><tr><td>GLDM (n = 14)</td><td>Dependence Entropy, Dependence Non-Uniformity, Dependence Non-Uniformity Normalized, Dependence Variance, Gray-Level Non-Uniformity, Gray-Level Variance, High Gray-Level Emphasis, Large Dependence Emphasis, Large Dependence High Gray-Level Emphasis, Large Dependence Low Gray-Level Emphasis, Low Gray-Level Emphasis, Small Dependence Emphasis, Small Dependence High Gray-Level Emphasis, Small Dependence Low Gray-Level Emphasis</td></tr><tr><td>High-order statistics (n = 1584)</td><td>First-Order and Second-Order Features Are Transformed by LoG, Exponential, Square, Square Root, Logarithm, Wavelet (Wavelet-LHL, Wavelet-LHH, Wavelet-HLL, Wavelet-LLH, Wavelet-HLH, Wavelet-HHH, Wavelet-HHL, Wavelet-LLL)</td></tr></tbody></table>

GLCM: gray-level co-occurrence matrix; GLRLM: gray-level run length matrix; GLSZM: gray-level size zone matrix; L; low, H: high.

## 2.6. Statistical Analysis and Feature Selection

Three steps were followed to select the most relevant features. In the first step, a variance criterion was set to decrease the number of features that did not meet the consistency threshold, which was set to 0.8 [20]. Variance thresholds remove features whose values do not change much from observation to observation (i.e., their variance falls below a threshold). These features provide little value. In the second step, based on the Student's t-test, the features that did not show statistical differences ($p > 0.05$) were removed. The t-test is any statistical hypothesis test in which the test statistic follows a Student's t-distribution under the null hypothesis [25]. This process can be viewed as preconditioning the predictive model [20]. The last step was penalization with the lasso algorithm [26]. The best value of lambda was identified, and the variance inflation factor (VIF) was calculated to obtain the most relevant features. The VIF measures the severity of multicollinearity in multiple regression models. It represents the ratio of the estimator variance of the regression coefficient to the variance when no linear correlation between the independent variables is assumed. The VIF can be calculated as follows:

$$VIF = \frac{1}{1 - R_i^2} \qquad (1)$$

where $R_i$ is the negative correlation coefficient of the regression analysis for other independent variables. The larger the VIF, the greater the possibility of collinearity between the independent variables. Generally, multicollinearity is assumed when the VIF value is greater than five. Thus, removing the radiomic features with a VIF value greater than five is necessary. Statistical outcomes for the variable of age were derived from a t-test analysis, while the results for the variable of sex were obtained using a *chi-squared* test. To build the classification model, the segmentation was labeled as two categories: tumor tissue and normal tissue. All steps were implemented using the R software version 4.3.2.

## 2.7. Models and Evaluation

Six supervised machine learning techniques were used for the classification, logistic regression (LR), random forests (RF), support vector machine (SVM), extreme gradient boosting (XGBoost), adaptive boosting (AdaBoost), and naive Bayes. The 10-fold cross-validation approach was used with random data splitting (70% for training and 30% for test). For the hyper-parameters' optimization, a grid search was used for each model [27]. The optimal parameters were found as follows:

* SVM: "C" = 0.05050505;

* Naive Bayes: "laplace" = 0, "adjust" = 1.5, "usekernel" = TRUE;
* XGBoost: "rounds" = 100, "max_depth" = 3, "eta" = 0.01, "gamma" = 0.2, "colsample_bytree" = 0.5, "min_child_weight" = 3, "subsample" = 1;

* RF: "mtry" = 21;

* LR: "alpha" = 1, "lambda" = 0.1;

* AdaBoost: "iter" = 150, "maxdepth" = 5, "nu" = 1.

The models' performances were evaluated using the area under the receiver operating curve (AUROC), sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and Matthews correlation coefficient (MCC). The best radiomics model was then screened. The AUROC predictive model's performance is presented as 95% corrected confidence intervals (CIs). The normal liver class was used as the positive class by the models.

## 3. Results

### 3.1. Clinical Factors of Patients

According to the patient's clinical records, sex and age were used as clinical characteristics. Table 3 shows the statistical results between the tumor and non-tumor groups.

<table><caption>**Table 3.** Statistical resultsof clinical characteristics between tumor and non-tumor groups.</caption><thead><tr><th>Clinical Characteristics</th><th>Tumor</th><th>Non-Tumor</th><th>p-Value</th></tr></thead><tbody><tr><td>Sex</td><td></td><td></td><td>0.42</td></tr><tr><td>Male</td><td>34 (36%)</td><td>21 (22%)</td><td></td></tr><tr><td>Female</td><td>20 (21%)</td><td>19 (20%)</td><td></td></tr><tr><td>Age</td><td>55.6 ± 16</td><td>51.3 ± 15</td><td>0.182</td></tr></tbody></table>

### 3.2. Radiomic Feature Selection

The dataset comprised diverse tumor contrast levels, variations in tissue size abnormalities, and various quantities of lesions. Validation was performed for the segmentations of each CT scan. There were 1686 features, including original, exponential, square, square root, logarithm, LoG, and wavelet transformations, extracted, as presented in Table 2. After calculations based on the variance threshold, 804 features underwent the first dimensionality reduction. Then, features with a *p*-value > 0.05 were excluded by computing the Student's *t*-test, and 487 features were selected. Among them, 28 features with non-zero coefficients were retained after lasso logistic regression analysis. Their VIF was calculated, and the variables with a *VIF* > 5 were removed to avoid model over-fitting. A total of seven radiomics features were retained as a result. The VIFs were all smaller than 5, indicating no multicollinearity among the seven radiomic features. Table 4 presents these variables with their VIF values. Figure 3 compares the selected variables between liver and tumor tissue. A correlation heat map of the selected radiomic features is presented in Figure 4.

<table><caption>**Table 4.** Model collinearity analysis</caption><thead><tr><th>Features</th><th>VIF</th></tr></thead><tbody><tr><td>log.sigma.0.5.mm.3D_firstorder_Maximum</td><td>1.226227</td></tr><tr><td>log.sigma.2.0.mm.3D_glszm_SizeZoneNonUniformity</td><td>4.81969</td></tr><tr><td>log.sigma.5.0.mm.3D_glszm_SizeZoneNonUniformity</td><td>3.135544</td></tr><tr><td>square_glszm_LargeAreaEmphasis</td><td>2.184934</td></tr><tr><td>logarithm_gldm_DependenceVariance</td><td>1.688747</td></tr><tr><td>exponential_glcm_ClusterProminence</td><td>3.00102</td></tr><tr><td>exponential_glcm_ClusterShade</td><td>1.598127</td></tr></tbody></table>
exponential_glcm_ClusterShade

<figure><img src="image_6.png" /></figure>

log sigma 0.5 mm 3D firstorder Maximum

<figure><img src="image_7.png" /></figure>

logsigma 5.0.mm.3D_glszm
_sizeZoneNonUniformity

<figure><img src="image_8.png" /></figure>

square_glszm_LargeAreaEmphasis

<figure><img src="image_9.png" /></figure>

logarithm_gldm_DependenceVariance

<figure><img src="image_10.png" /></figure>

exponential_glcm_ClusterProminence

<figure><img src="image_11.png" /></figure>

logsigma 2.0.mm.3D_glszm_SizeZoneNonUniformity

<figure><img src="image_12.png" /></figure>

Healthy liver

Tumor liver

<caption>**Figure 3.** Selected radiomic features compared between tumor and non-tumor tissue.</caption>
<figure><img src="image_13.png" /><figcaption>**Figure 4.** Heat map correlation of the selected features.</figcaption></figure>

### 3.3. Construction and Evaluation of the Classification Model of Tumor Presence or Not

The naive Bayes model achieved the best performance among the other algorithms, with an AUROC of 0.9268 (95% CI: 0.8224–1), sensitivity of 0.9091, specificity of 0.9444, NPV of 0.9444, PPV of 0.9091, and MCC of 0.8536. AdaBoost and XGBoost performed similarly, with an AUROC of 0.8813 (95% CI: 0.75–1). SVM was the worst algorithm, with an AUROC of 0.8258 (95% CI: 0.677–0.9745). All scores obtained during the test phase are presented in Table 5. The ROC curves of all models are plotted in Figure 5. Twenty-two radiomic features showed a higher importance correlation and are presented in Table 6.

<table><caption>**Table 5.** Metrics' measures for each ML algorithm.</caption><thead><tr><th>Algorithm</th><th>AUROC (95% CI)</th><th>Sensitivity</th><th>Specificity</th><th>NPV</th><th>PPV</th><th>MCC</th></tr></thead><tbody><tr><td>SVM</td><td>0.8258 (0.677–0.9745)</td><td>0.8182</td><td>0.8333</td><td>0.75</td><td>0.8824</td><td>0.6419</td></tr><tr><td><strong>Naive Bayes</strong></td><td><strong>0.9268 (0.8224–1)</strong></td><td><strong>0.9091</strong></td><td><strong>0.9444</strong></td><td><strong>0.9444</strong></td><td><strong>0.9091</strong></td><td><strong>0.8536</strong></td></tr><tr><td>XGBoost</td><td>0.8813 (0.75–1)</td><td>0.8182</td><td>0.9444</td><td>0.8947</td><td>0.9</td><td>0.7785</td></tr><tr><td>RF</td><td>0.8535 (0.7126–0.9945)</td><td>0.8182</td><td>0.8889</td><td>0.8889</td><td>0.8182</td><td>0.7071</td></tr><tr><td>Logistic</td><td>0.8359 (0.6875–0.9842)</td><td>0.7273</td><td>0.9444</td><td>0.8889</td><td>0.85</td><td>0.7045</td></tr><tr><td>AdaBoost</td><td>0.8813 (0.75–1)</td><td>0.8182</td><td>0.9444</td><td>0.8947</td><td>0.9</td><td>0.7785</td></tr></tbody></table>

CI: confidence interval.
<figure><img src="image_14.png" /><figcaption>**Figure 5.** The ROC curves of the classification model of tumor presence or not.</figcaption></figure>

<table><caption>**Table 6.** The most important radiomic features in tumor or no tumor classification.</caption><thead><tr><td>Feature</td><td>Importance</td></tr></thead><tbody><tr><td>log.sigma.2.0.mm.3D_firstorder_Median</td><td>0.198357587</td></tr><tr><td>wavelet.LLL_glcm_Autocorrelation</td><td>0.139106233</td></tr><tr><td>log.sigma.2.0.mm.3D_firstorder_Kurtosis</td><td>0.066871147</td></tr><tr><td>log.sigma.2.0.mm.3D_firstorder_InterquartileRange</td><td>0.061976577</td></tr><tr><td>logarithm_glcm_ClusterShade</td><td>0.055246304</td></tr></tbody></table>
<table><caption>**Table 6.** *Cont.*</caption><thead><tr><td>Feature</td><td>Importance</td></tr></thead><tbody><tr><td>wavelet.LHL_firstorder_Maximum</td><td>0.049298670</td></tr><tr><td>wavelet.LLH_gldm_SmallDependenceHighGrayLevelEmphasis</td><td>0.043135931</td></tr><tr><td>exponential_firstorder_MeanAbsoluteDeviation</td><td>0.033085765</td></tr><tr><td>log.sigma.1.0.mm.3D_glcm_SumAverage</td><td>0.029689345</td></tr><tr><td>original_glszm_HighGrayLevelZoneEmphasis</td><td>0.028640204</td></tr><tr><td>log.sigma.1.0.mm.3D_glszm_HighGrayLevelZoneEmphasis</td><td>0.024491582</td></tr><tr><td>wavelet.LLH_glrlm_LongRunHighGrayLevelEmphasis</td><td>0.023100819</td></tr><tr><td>wavelet.LLL_firstorder_Kurtosis</td><td>0.022962333</td></tr><tr><td>wavelet.HHH_glszm_LargeAreaLowGrayLevelEmphasis</td><td>0.022767771</td></tr><tr><td>log.sigma.2.0.mm.3D_glszm_SizeZoneNonUniformity</td><td>0.022381763</td></tr><tr><td>log.sigma.2.0.mm.3D_firstorder_MeanAbsoluteDeviation</td><td>0.019902419</td></tr><tr><td>wavelet.LHL_glcm_ClusterShade</td><td>0.018137160</td></tr><tr><td>original_firstorder_Mean</td><td>0.016972684</td></tr><tr><td>log.sigma.0.5.mm.3D_glrlm_LongRunHighGrayLevelEmphasis</td><td>0.016380692</td></tr><tr><td>log.sigma.5.0.mm.3D_glszm_SmallAreaHighGrayLevelEmphasis</td><td>0.013439397</td></tr><tr><td>log.sigma.1.0.mm.3D_glszm_GrayLevelNonUniformity</td><td>0.013436358</td></tr><tr><td>original_glrlm_RUNLengthNonUniformity</td><td>0.011691586</td></tr></tbody></table>

### 3.4. Construction and Evaluation of the Classification Model of Benign or Malignancy Tumor

The naive Bayes model achieved the best performance again, among the other algorithms, with an AUROC of 0.8571 (95% CI: 0.6764–1), sensitivity of 0.7143, specificity of 1, NPV of 0.8333, PPV of 1, and MCC of 0.7715. It was followed by the lasso regression with an AUROC of 0.7569 (95% CI: 0.5472–0.9667). All other models performed similarly, with an AUROC of 0.5929 (95% CI: 0.3873–0.7985). All scores obtained during the test phase are presented in Table 7. The ROC curves of all models are plotted in Figure 6. Eighteen radiomic features showed a higher importance correlation and are presented in Table 8.

<table><caption>**Table 7.** Metrics' measures for each ML algorithm in benign or malignancy tumor classification.</caption><thead><tr><td>Algorithm</td><td>AUROC (95% CI)</td><td>Sensitivity</td><td>Specificity</td><td>NPV</td><td>PPV</td><td>MCC</td></tr></thead><tbody><tr><td>SVM</td><td>0.5929 (0.3873–0.7985)</td><td>0.2857</td><td>0.9</td><td>0.6429</td><td>0.6667</td><td>0.2398</td></tr><tr><td><strong>Naive Bayes</strong></td><td><strong>0.8571 (0.6764–1)</strong></td><td><strong>0.7143</strong></td><td><strong>1</strong></td><td><strong>0.8333</strong></td><td><strong>1</strong></td><td><strong>0.7715</strong></td></tr><tr><td>XGBoost</td><td>0.5929 (0.3873–0.7985)</td><td>0.2857</td><td>0.9</td><td>0.6429</td><td>0.6667</td><td>0.2398</td></tr><tr><td>RF</td><td>0.5929 (0.3873–0.7985)</td><td>0.2857</td><td>0.9</td><td>0.6429</td><td>0.6667</td><td>0.2398</td></tr><tr><td>Logistic</td><td>0.7569 (0.5472–0.9667)</td><td>0.6250</td><td>0.8889</td><td>0.7273</td><td>0.8333</td><td>0.5367</td></tr><tr><td>AdaBoost</td><td>0.5929 (0.3873–0.7985)</td><td>0.2857</td><td>0.9</td><td>0.6429</td><td>0.6667</td><td>0.2398</td></tr></tbody></table>

CI: confidence interval.

<table><caption>**Table 8.** The most important radiomic features in benign or malignancy tumor classification.</caption><thead><tr><td>Feature</td><td>Importance</td></tr></thead><tbody><tr><td>exponential_glrlm_HighGrayLevelRunEmphasis</td><td>0.3296813462</td></tr><tr><td>originalShape_Maximum2DDiameterColumn</td><td>0.1837066150</td></tr><tr><td>originalShape_Maximum2DDiameterRow</td><td>0.0666429693</td></tr><tr><td>originalShape_Maximum3DDiameter</td><td>0.0603642515</td></tr><tr><td>exponential_glrlm_GrayLevelVariance</td><td>0.0464352022</td></tr><tr><td>square_firstorder_RobustMeanAbsoluteDeviation</td><td>0.0389244507</td></tr><tr><td>exponential_glcm_ClusterTendency</td><td>0.0284419206</td></tr><tr><td>exponential_glszm_HighGrayLevelZoneEmphasis</td><td>0.0253803621</td></tr><tr><td>exponential_firstorder_RootMeanSquared</td><td>0.0226305816</td></tr><tr><td>squareroot_glszm_HighGrayLevelZoneEmphasis</td><td>0.0206921639</td></tr><tr><td>wavelet.LHL_firstorder_Range</td><td>0.0195810157</td></tr><tr><td>wavelet.LLH_gldm_GrayLevelNonUniformity</td><td>0.0187361912</td></tr><tr><td>exponential_glszm_GrayLevelNonUniformity</td><td>0.0173371526</td></tr></tbody></table>
<table><caption>**Table 8.** *Cont.*</caption><thead><tr><td>Feature</td><td>Importance</td></tr></thead><tbody><tr><td>original_gldm_GrayLevelNonUniformity</td><td>0.0164583230</td></tr><tr><td>log.sigma.3.0.mm.3D_glszm_LargeAreaLowGrayLevelEmphasis</td><td>0.0151196504</td></tr><tr><td>square_glrlm_ShortRunHighGrayLevelEmphasis</td><td>0.0127392186</td></tr><tr><td>squareroot_firstorder_90Percentile</td><td>0.0121420080</td></tr><tr><td>logarithm_firstorder_90Percentile</td><td>0.0119382485</td></tr></tbody></table>

XGBoost ROC curve, AUC = 0.5929

<figure><img src="image_15.png" /></figure>

AdaBoost ROC curve, AUC = 0.5929

<figure><img src="image_16.png" /></figure>

Logistic regression ROC curve, AUC = 0.7569

<figure><img src="image_17.png" /></figure>

Naive Bayes ROC curve, AUC = 0.8571

<figure><img src="image_18.png" /></figure>

Random forest ROC curve, AUC = 0.5929

<figure><img src="image_19.png" /></figure>

SVM ROC curve, AUC = 0.5929

<figure><img src="image_20.png" /><figcaption>**Figure 6.** The ROC curves of the classification model of benign or malignancy tumor.</figcaption></figure>
## 4. Discussion

The study shows that machine learning algorithms can capture the difference between normal and tumor liver tissue and can detect benign or malignancy tumors in CT images at the portal phase using radiomic features. In radiomics, a significant array of features, often numbering in the hundreds, is generated to delineate a specified ROI through various methodologies [28]. These features undergo evaluation for their potential as prognostic indicators. Additionally, the critical process of feature selection demands a focus on their consistency and responsiveness concerning the delineation methodology for them to be deemed clinically viable. To date, no studies have systematically assessed the application of various machine learning algorithms in differentiating between tumor and non-tumor tissues using the entire liver as the ROI. Furthermore, no African recorded studies employ machine learning to classify liver lesions.

The approach used six machine learning algorithms, enhanced with hyper-parameter tuning through grid search and validated via 10-fold cross-validation to identify the most effective model based on the extracted features. This method ensures the selection of the best model by systematically exploring a range of parameter settings and rigorously evaluating model performance across multiple subsets of the data. The study revealed that the naive Bayes algorithm outperformed the others, achieving the highest marks across the AUROC, sensitivity, specificity, NPV, PPV, and MCC metrics for the two classification tasks. This demonstrates the reliability of the data and the stability of the naive Bayes model in this study. In detecting the presence of a tumor or not, AdaBoost and XGBoost were closely matched as the second-best performing algorithms, demonstrating equivalent outcomes in their performance evaluations. However, all models performed well, with no drastic difference in the scores.

In contrast, only naive Bayes and lasso regression yielded satisfactory scores for detecting malignant or benign tumors. The other models produced low scores. In total, 1686 features were extracted from portal-phase CT images. Twenty-two features showed a higher correlation in the presence of tumor or not classification, while there were eighteen in the classification of benign or malignancy tumor.

Applying lasso regression for variable selection identified seven variables with a VIF below 5. The correlation map in Figure 4 demonstrates no correlation between features. The boxplot in Figure 3 demonstrates that most radiomic features extracted from tumor tissue cover a broader spectrum than those derived from normal liver tissue (5 features out of 7). The study of Lysdahlgaard observed the same phenomenon [15]. Zhou et al. [29] conducted a similar study, classifying gross tumor volume and normal liver tissue in HCC. They also identified seven variables with a VIF below 5. However, their findings differed in model performance; they identified naive Bayes as the least effective model, whereas XGBoost emerged as the superior algorithm in terms of AUC performance.

Age and sex were added to the selected features. It should be noted that there was no statistically significant difference in age and sex between healthy livers and those with tumors (*p* > 0.05). However, age had a *p*-value of 0.182, suggesting that it may merit more examination in future studies. This study demonstrates the feasibility of using the entire liver as the region of interest for segmentation rather than extracting the local lesion for observations. It can be useful for automatically and swiftly identifying ill patients. It is particularly interesting because liver and lesion segmentation research yields more accurate results when segmenting the liver than the lesions themselves.

Nevertheless, it remains important to distinguish between lesions such as HCC, hepatic cyst, HEM, FNH, HA, etc., to administer appropriate treatments and avoid unnecessary procedures. Stratifying tumors for specific lesion types would have been ideal; however, the study's small sample size limited this approach. Additionally, complementing hepatic CT with ultrasound or MRI could further enhance diagnostics, potentially unveiling new radiomic features from these modalities and improving lesion characterization and treatment accuracy. A study conducted by Yao et al. explored 2560 radiomic features from 177 patients' multimodal ultrasound images and found that radiomic models could be
useful for evaluating hepatic tumors, especially for diagnosis, differential diagnosis, and clinical prognosis [30]. By comparing the radiomic properties of ultrasound, MRI, and computed tomography, it may be possible to predict and differentiate between various types of hepatic lesions.

## 5. Conclusions

This research performed a preliminary radiomics investigation using CT scans. A correlation has been established between radiomic features and the distinction between healthy and tumor liver tissue, on the one hand, and between malignant and benign tumors, on the other hand. The naive Bayes model showed great results for the two tasks. It achieved an AUROC of 0.9268 (95% CI: 0.8224–1) in the distinction between healthy and tumor liver tissue and an AUROC of 0.8571 (95% CI: 0.6764–1) in the distinction between malignant and benign tumors. It demonstrated that the developed radiomics signature correlated statistically strongly with healthy and tumor liver tissue using the entire liver as the ROI. However, there is a need for larger retrospective and prospective research on liver CT scans that examines potential prognostic markers and has strict reference criteria.

**Author Contributions:** Conceptualization, R.S.A.; methodology, R.S.A. and F.H.R.G.; software, V.-B.S.Z.; validation, R.S.A., F.d.V. and E.C.E.; formal analysis, V.-B.S.Z. and F.H.R.G.; investigation, V.-B.S.Z. and F.H.R.G.; resources, O.B.; data curation, V.-B.S.Z.; writing—original draft preparation, V.-B.S.Z.; writing—review and editing, E.C.E. and R.S.A.; visualization, R.S.A., O.B. and F.H.R.G.; supervision, R.S.A., O.B. and E.C.E.; project administration, R.S.A., O.B. and E.C.E.; funding acquisition, R.S.A. and V.-B.S.Z. All authors have read and agreed to the published version of the manuscript.

**Funding:** The French Government funded this research through its Embassy in the Benin Republic.

**Institutional Review Board Statement:** The study was conducted following the Declaration of Helsinki and approved by the Ethics Committee of Université de Parakou in Benin Republic (protocol code 0510/CLERB-UP/P/SP/R/SA of 28 February 2022).

**Informed Consent Statement:** Patient consent was waived due to the nature of the retrospective study.

**Data Availability Statement:** The data supporting the conclusions of this article are available upon request from the corresponding author. These data are not publicly accessible due to restrictions on data dissemination imposed by the Ethics Committee.

**Conflicts of Interest:** The authors declare no conflicts of interest.

## Abbreviations

The following nomenclature are used in this manuscript:

<table><tr><td>AdaBoost</td><td>adaptive boosting</td></tr><tr><td>AUROC</td><td>area under the receiver operating curve</td></tr><tr><td>CI</td><td>confidence interval</td></tr><tr><td>CNHU-HKM</td><td>Centre National Hospitalier Universitaire Hubert Koutoukou Maga</td></tr><tr><td>CT</td><td>computed tomography</td></tr><tr><td>FNH</td><td>focal nodular hyperplasia</td></tr><tr><td>GLCM</td><td>gray-level co-occurrence matrix</td></tr><tr><td>GLDM</td><td>gray-level dependence matrix</td></tr><tr><td>GLH</td><td>gray-level histogram</td></tr><tr><td>GLRLM</td><td>gray-level run length matrix</td></tr><tr><td>GLSZM</td><td>gray-level size zone matrix</td></tr><tr><td>H</td><td>high</td></tr><tr><td>HA</td><td>hepatocellular adenoma</td></tr><tr><td>HCC</td><td>hepatocellular carcinoma</td></tr><tr><td>HEM</td><td>hepatic hemangioma</td></tr><tr><td>HU</td><td>Hounsfield unit</td></tr><tr><td>L</td><td>low</td></tr><tr><td>LoG</td><td>Laplacian of Gaussian</td></tr></table>
<table><tr><td>LR</td><td>logistic regression</td></tr><tr><td>MCC</td><td>Matthews correlation coefficient</td></tr><tr><td>ML</td><td>machine learning</td></tr><tr><td>MRI</td><td>magnetic resonance imaging</td></tr><tr><td>NPV</td><td>negative prediction value</td></tr><tr><td>PNN</td><td>probabilistic neural network</td></tr><tr><td>PPV</td><td>positive prediction value</td></tr><tr><td>RCAD</td><td>radiomics-based computer-aided</td></tr><tr><td>RF</td><td>random forests</td></tr><tr><td>RIL-Contour</td><td>Radiology Informatics Laboratory Contour</td></tr><tr><td>ROI</td><td>region of interest</td></tr><tr><td>SVM</td><td>support vector machine</td></tr><tr><td>VIF</td><td>variance inflation factor</td></tr><tr><td>XGBoost</td><td>extreme gradient boosting</td></tr></table>

References

1. Sung, H.; Ferlay, J.; Siegel, R.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. *CA Cancer J. Clin.* **2021**, *71*, 209–249. [CrossRef]

2. Hann, L.; Winston, C.; Brown, K.; Akhurst, T. Diagnostic imaging approaches and relationship to hepatobiliary cancer staging and therapy. *Semin. Surg. Oncol.* **2000**, *19*, 94–115. [CrossRef]

3. Lambin, P.; Rios-Velazquez, E.; Leijenaar, R.; Carvalho, S.; Van Stiphout, R.; Granton, P.; Zegers, C.; Gillies, R.; Boellard, R.; Dekker, A.; et al. Radiomics: Extracting more information from medical images using advanced feature analysis. *Eur. J. Cancer* **2012**, *48*, 441–446. [CrossRef]

4. Gillies, R.; Kinahan, P.; Hricak, H. Radiomics: Images Are More than Pictures, They Are Data. *Radiology* **2016**, *278*, 563–577. [CrossRef] [PubMed]

5. Balogh, J.; Victor, D.; Asham, E.; Burroughs, S.; Boktour, M.; Saharia, A.; Li, X.; Ghobrial, M.; Monsour, H. Hepatocellular carcinoma: A review. *J. Hepatocell. Carcinoma* **2016**, *3*, 41–53. [CrossRef] [PubMed]

6. Grazioli, L.; Bondioni, M.P.; Haradome, H.; Motosugi, U.; Tinti, R.; Frittoli, B.; Gambarini, S.; Donato, F.; Colagrande, S. Hepatocellular adenoma and focal nodular hyperplasia: Value of gadoxetic acid-enhanced MR imaging in differential diagnosis. *Radiology* **2012**, *262*, 520–529. [CrossRef] [PubMed]

7. Mitchell, D.; Bruix, J.; Sherman, M.; Sirlin, C. LI-RADS (Liver Imaging Reporting and Data System): Summary, discussion, and consensus of the LI-RADS Management Working Group and future directions. *Hepatology* **2015**, *61*, 1056–1065. [CrossRef] [PubMed]

8. Starmans, M.; Klein, S.; Van Der Voort, S.; Thomeer, M.; Miclea, R.; Niessen, W. Classification of malignant and benign liver tumors using a radiomics approach. In *Medical Imaging 2018: Image Processing*; SPIE: Bellingham, WA, USA, 2018; Volume 10574, p. 48.

9. Mougiakakou, S.; Valavanis, I.; Nikita, A.; Nikita, K. Differential diagnosis of CT focal liver lesions using texture features, feature selection and ensemble driven classifiers. *Artif. Intell. Med.* **2007**, *41*, 25–37. [CrossRef] [PubMed]

10. Kumar, S.; Moni, R.; Rajeesh, J. An automatic computer-aided diagnosis system for liver tumours on computed tomography images. *Comput. Electr. Eng.* **2013**, *39*, 1516–1526. [CrossRef]

11. Qiu, J.; Wu, Y.; Hui, B.; Chen, J.; Ji, L.; Wang, M. A Novel Texture Analysis Method Based on Reverse Biorthogonal Wavelet and Co-Occurrence Matrix Applied for Classification of Hepatocellular Carcinoma and Hepatic Hemangioma. *J. Med. Imaging Health Inform.* **2018**, *8*, 1835–1843. [CrossRef]

12. Zhou, Y.; He, L.; Huang, Y.; Chen, S.; Wu, P.; Ye, W.; Liu, Z.; Liang, C. CT-based radiomics signature: A potential biomarker for preoperative prediction of early recurrence in hepatocellular carcinoma. *Abdom. Radiol.* **2017**, *42*, 1695–1704. [CrossRef]

13. Nie, P.; Yang, G.; Guo, J.; Chen, J.; Li, X.; Ji, Q.; Wu, J.; Cui, J.; Xu, W. A CT-based radiomics nomogram for differentiation of focal nodular hyperplasia from hepatocellular carcinoma in the non-cirrhotic liver. *Cancer Imaging* **2020**, *20*, 20. [CrossRef]

14. Nie, P.; Wang, N.; Pang, J.; Yang, G.; Duan, S.; Chen, J.; Xu, W. CT-Based Radiomics Nomogram: A Potential Tool for Differentiating Hepatocellular Adenoma From Hepatocellular Carcinoma in the Noncirrhotic Liver. *Acad. Radiol.* **2021**, *28*, 799–807. [CrossRef]

15. Lysdahlgaard, S. Comparing Radiomics features of tumour and healthy liver tissue in a limited CT dataset: A machine learning study. *Radiography* **2022**, *28*, 718–724. [CrossRef]

16. Thawani, R.; McLane, M.; Beig, N.; Ghose, S.; Prasanna, P.; Velcheti, V.; Madabhushi, A. Radiomics and radiogenomics in lung cancer: A review for the clinician. *Lung Cancer* **2018**, *115*, 34–41. [CrossRef]

17. Wakabayashi, T.; Ouhmich, F.; Gonzalez-Cabrera, C.; Felli, E.; Saviano, A.; Agnus, V.; Savadjiev, P.; Baumert, T.; Pessaux, P.; Marescaux, J.; et al. Radiomics in hepatocellular carcinoma: A quantitative review. *Hepatol. Int.* **2019**, *13*, 546–559. [CrossRef] [PubMed]

18. Philbrick, K.; Weston, A.; Akkus, Z.; Kline, T.; Korfiatis, P.; Sakinis, T.; Kostandy, P.; Boonrod, A.; Zeinoddini, A.; Takahashi, N.; et al. RIL-Contour: A Medical Imaging Dataset Annotation Tool for and with Deep Learning. *J. Digit. Imaging* **2019**, *32*, 571–581. [CrossRef] [PubMed]
19. DenOtter, T.; Schubert, J. *Hounsfield Unit*; StatPearls: Treasure Island, FL, USA, 2024.

20. Wu, J.; Liu, A.; Cui, J.; Chen, A.; Song, Q.; Xie, L. Radiomics-based classification of hepatocellular carcinoma and hepatic haemangioma on precontrast magnetic resonance images. *BMC Med. Imaging* **2019**, *19*, 23. [CrossRef] [PubMed]

21. Thibault, G.; Angulo, J.; Meyer, F. Advanced statistical matrices for texture characterization: Application to cell classification. *IEEE Trans. Bio-Med. Eng.* **2014**, *61*, 630–637. [CrossRef]

22. Kairuddin, W.; Mahmud, W. Texture Feature Analysis for Different Resolution Level of Kidney Ultrasound Images. *IOP Conf. Ser. Mater. Sci. Eng.* **2017**, *226*, 012136. [CrossRef]

23. Thibault, G.; Fertil, B.; Navarro, C.; Pereira, S.; Ley, N.; Sequeira, J.; Mari, J. Texture Indexes and Gray Level Size Zone Matrix Application to Cell Nuclei Classification. In Proceedings of the Pattern Recognition and Information Processing (PRIP), Minsk, Belarus, 19–21 May 2009; pp. 140–145.

24. Khalil, M.M. Radiomics Feature Selection in 18F-FDG PET Imaging: Investigation of the Most robust and Reproducible Candidates. *Egypt. J. Nucl. Med.* **2020**, *21*, 80–102.

25. Student. The probable error of a mean. *Biometrika* **1908**, *6*, 1–25. [CrossRef]

26. Lindstrom, J.; Dahl, F. Model Selection with Lasso in Multi-group Structural Equation Models. *Struct. Equ. Model. Multidiscip. J.* **2019**, *27*, 33–42. [CrossRef]

27. Lecun, Y.; Bottou, L.; Bengio, Y.; Haffner, P. Gradient-based learning applied to document recognition. *Proc. IEEE* **1998**, *86*, 2278–2324. [CrossRef]

28. Owens, C.; Peterson, C.; Tang, C.; Koay, E.; Yu, W.; Mackin, D.; Li, J.; Salehpour, M.; Fuentes, D.; Court, L.; et al. Lung tumor segmentation methods: Impact on the uncertainty of radiomics features for non-small cell lung cancer. *PLoS ONE* **2018**, *13*, e0205003. [CrossRef] [PubMed]

29. Zhang, H.; Huang, D.; Wang, Y.; Zhong, H.; Pang, H. CT radiomics based on different machine learning models for classifying gross tumor volume and normal liver tissue in hepatocellular carcinoma. *Cancer Imaging* **2024**, *24*, 20. [CrossRef] [PubMed]

30. Yao, Z.; Dong, Y.; Wu, G.; Zhang, Q.; Yang, D.; Yu, J.; Wang, W. Preoperative diagnosis and prediction of hepatocellular carcinoma: Radiomics analysis based on multi-modal ultrasound images. *BMC Cancer* **2018**, *18*, 1089. [CrossRef] [PubMed]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.