---
title: Handcrafted Radiomics Features
url: https://www.emergentmind.com/topics/handcrafted-radiomics-features-hrf
type: topic
---

# Handcrafted Radiomics Features

Handcrafted radiomics features (HRFs), also called engineered features, are quantitative descriptors computed from predefined image-derived measurements extracted from a segmented region or volume of interest. They convert imaging data into structured, high-dimensional feature representations while retaining closed-form mathematical definitions for intensity histograms, shape indices, texture descriptors, filter responses, and related biomarkers. In contrast to deep radiomics, HRFs are user-driven, mathematically explicit, and typically more interpretable, but their reliability depends on harmonized acquisition, segmentation, preprocessing, and validation protocols [2602.00102] [1808.07954] [2407.13813].

## 1. Definition and conceptual position

In the standardized radiomics pipeline, HRFs fall into three major families—first-order statistics, shape descriptors, and texture features—each encoding distinct aspects of image data. First-order statistics summarize the intensity distribution inside the ROI without regard to spatial arrangement; shape descriptors quantify geometry independent of intensity; texture features encode spatial relationships between gray levels through statistical matrices such as the gray-level co-occurrence matrix (GLCM). Broader implementations additionally include higher-order features derived from wavelets, Laplacian-of-Gaussian or Gabor filtering, as well as moment-invariant and diagnostic features in software frameworks such as PySERA [2602.00102] [1808.07954] [2511.15963].

HRFs occupy a distinct methodological position relative to deep radiomics. Handcrafted radiomics is a user-driven workflow in which an ROI is defined, a fixed set of mathematically defined features is computed, and feature selection and modeling follow; deep radiomics embeds feature extraction, selection, and modeling in one or more neural-network layers. The trade-off reported across reviews is consistent: HRFs are transparent, interpretable, and workable with smaller datasets, whereas deep radiomics can capture more complex patterns but demands larger training sets and yields less explicit feature semantics [2407.13813] [1808.07954].

A common misconception is that mathematical explicitness alone guarantees robustness. The survey literature instead emphasizes that reliable HRF computation hinges on preprocessing choices and parameter settings, and that scanner- and segmentation-induced variability can overwhelm true biological signal when protocols are not standardized [2602.00102]. This suggests that the defining property of HRFs is not simply interpretability, but interpretability under controlled extraction conditions.

## 2. Feature families and mathematical structure

First-order HRFs probe the intensity histogram of the ROI. For voxel intensities $\{x_i\}_{i=1}^N$, representative descriptors include the mean $\mu=\frac{1}{N}\sum_{i=1}^N x_i$, the variance $\sigma^2=\frac{1}{N}\sum_{i=1}^N (x_i-\mu)^2$, skewness, kurtosis, entropy, energy, uniformity, percentiles, and related summary statistics. Across reviews and application papers, these measures are used to characterize attenuation, signal intensity, enhancement behavior, and heterogeneity at the histogram level [2602.00102] [1808.07954].

Shape descriptors quantify lesion or organ geometry independently of gray-level intensity. Frequently used examples are volume, surface area, compactness, sphericity, surface-to-volume ratio, maximum diameter, and axis lengths. A representative sphericity definition is
$$
\Psi=\frac{\pi^{1/3}(6V)^{2/3}}{A},
$$
where $V$ is volume and $A$ is surface area. Such features are explicitly linked in the literature to irregularity, deviation from spherical morphology, and margin complexity [2507.08546] [2602.00102].

Texture HRFs model spatial relationships between gray levels. For a normalized GLCM $P(i,j)$, common statistics include contrast,
$$
\mathrm{Contrast}=\sum_{i,j}(i-j)^2P(i,j),
$$
correlation, energy $\sum_{i,j}P(i,j)^2$, homogeneity, dissimilarity, and entropy. Other matrix families extend the same logic to different spatial organizations: GLRLM measures runs of identical gray levels, GLSZM measures connected zones, NGTDM measures local gray-tone differences, and GLDM or NGLDM measure dependence structure [2602.00102] [1808.07954] [2511.15963].

Recent task-specific papers also extend the classical family structure. RadiomicNet uses Gray-Level Co-occurrence Matrix and Local Binary Pattern features, with LBP defined on a circular neighborhood by
$$
\mathrm{LBP}_{P,R}(x_c)=\sum_{p=0}^{P-1}s(I(x_p)-I(x_c))\,2^p,\qquad
s(u)=
\begin{cases}
1,&u\ge 0\\
0,&u<0
\end{cases}
$$
using $P=24$ and $R=3$; its 13-dimensional radiomics vector comprises 5 GLCM, 4 LBP, and 4 first-order intensity moments [2607.02185]. A different extension is Image Frequency Analysis, which computes Pearson-style correlations between adjacent rows or columns and then applies a one-dimensional Fourier transform to the resulting sequences, yielding a directional frequency descriptor of the image [2111.05855]. This suggests that HRFs are not limited to conventional ROI-level histogram and texture matrices, but also include explicitly designed descriptors of cross-line structure and local coding patterns.

## 3. Extraction pipeline, preprocessing, and parameterization

The canonical handcrafted-radiomics workflow comprises image acquisition and reconstruction, segmentation, preprocessing, feature extraction, and feature reduction or model building. Segmentation may be manual, semi-automatic, or fully automatic, and segmentation quality directly impacts feature stability. The preprocessing stage typically includes intensity normalization, spatial resampling to isotropic voxels, grayscale discretization, and, when required, denoising or artifact correction [1808.07954] [2602.00102].

Parameterization is not ancillary; it changes the feature being computed. Standard choices reported across surveys include spatial resampling to isotropic voxels such as $1\times1\times1\ \mathrm{mm}^3$, normalization by z-score, histogram matching, or Nyúl’s mapping, and gray-level discretization by either fixed bin width or fixed bin count. For texture matrices, the distance $d$, directions $\theta$, and aggregation strategy are all consequential: longer distances capture coarser patterns, while multi-directional averaging improves rotational invariance at the cost of interpretability [2602.00102] [1808.07954].

Concrete implementations make these dependencies explicit. RadiomicNet computes GLCM features over distances $d\in\{1,3\}$ and angles $\theta\in\{0^\circ,45^\circ,90^\circ,135^\circ\}$ with 64 gray-level quantization, and averages the five statistics over all $(d,\theta)$ pairs; its LBP uses radius $R=3$, $P=24$ neighbors, and uniform encoding [2607.02185]. In PySERA, standardized preprocessing includes resampling, discretization, normalization, optional range-based re-segmentation, optional outlier removal, and partial-volume binarization; the framework computes 557 handcrafted features, including 487 IBSI-compliant features, 10 moment-invariant descriptors, and 60 diagnostic features [2511.15963].

Software ecosystems formalize these workflows. The survey literature explicitly names PyRadiomics, IBEX, and CERR as radiomics software environments aligned with IBSI definitions, while PySERA extends this landscape with a Python-native, object-oriented implementation benchmarked against PyRadiomics, MITK, and LIFEx [2602.00102] [2511.15963]. In practice, the software choice is inseparable from the definition of the extracted feature set, because implementation details, defaults, and compliance level affect reproducibility.

## 4. Reproducibility, stability, and dimensionality control

Feature instability is one of the central methodological problems in handcrafted radiomics. Major sources of variability include scanner manufacturer, reconstruction kernel, slice thickness, and ROI delineation; segmentation uncertainty propagates directly into feature values. The survey literature identifies two standard reproducibility measures: the intraclass correlation coefficient (ICC), with features above $0.80$ deemed reproducible, and the coefficient of variation (CoV) [2602.00102].

Harmonization is therefore a structural component of HRF analysis rather than a post hoc convenience. Feature-domain methods such as ComBat and Bayesian ComBat are reported as ways to remove batch effects while retaining biological variance, and IBSI provides reference definitions and test-case images to standardize implementation across software. At the same time, modality-specific results show that no single harmonization method is uniformly optimal. In an echocardiography study on left-ventricle hypertension disease phenotyping, ComBat yielded mean Jensen–Shannon divergence $0.0735\pm0.0617$ versus $0.0585\pm0.0746$ on original images, whereas a self-supervised convolutional-kernel filter reduced it to $0.0470\pm0.0367$ and increased AUC from $0.637$ to $0.868$ [2310.08897].

Because HRF studies often operate in the $p\gg n$ regime, feature selection and dimensionality reduction are equally central. Reported strategies include filter methods such as variance thresholding, correlation analysis, mutual information, and mRMR; wrapper methods such as recursive feature elimination; embedded methods such as LASSO and Elastic Net; and dimensionality reduction by PCA, kernel PCA, or autoencoders. A recurring methodological requirement is that all selection and reduction steps must reside inside nested cross-validation loops to avoid data leakage and optimistic bias [2602.00102].

The instability induced by multicollinearity has motivated dedicated statistical frameworks. One example is the four-step approach of redundancy filtering, penalized maximum-likelihood estimation of a filtered correlation matrix, maximum-likelihood factor analysis, and factor scoring, producing a compact latent-feature space for regression-based prediction. This approach was explicitly proposed to address the high dimensionality, extreme multicollinearity, and instability of standard penalized classifiers in radiomics data [1903.11696]. A plausible implication is that stable HRF modeling often depends as much on covariance structure control as on the semantic content of the original descriptors.

## 5. Hybridization with deep learning and task-specific deployment

A major recent development is the direct integration of HRFs into deep architectures rather than their use as separate tabular covariates. RadiomicNet is a two-stream hybrid segmentation network in which a radiomics embedding modulates skip-connection attention through a Radiomics Attention Gate, and a Radiomics Consistency Loss aligns texture complexity with prediction entropy. With 13 handcrafted radiomics features as input to the radiomics branch, the model achieved a Dice Similarity Coefficient of $0.763\pm0.231$ on BUSI and $0.854\pm0.112$ on Kvasir-SEG, reduced Expected Calibration Error on BUSI from $0.142$ to $0.118$, and used only $3.27$M parameters. Gradient-based feature importance identified GLCM dissimilarity ($15.24\%$), GLCM energy ($14.56\%$), and LBP entropy ($11.49\%$) as the dominant radiomics cues [2607.02185].

Hybridization can also occur at both global and voxel levels. In pancreatic ductal adenocarcinoma detection, a two-stage nnUNet-based workflow first extracted 1,486 global radiomic descriptors over the pancreas mask, then selected ten features by a Pearson correlation test with false discovery rate control and recursive feature elimination wrapped around a linear SVM, and finally injected eight non-shape features as voxel-level parametric maps while using the ten-dimensional global vector in radiomics-aware cross-attention at the bottleneck. The resulting method achieved AUC $=0.96$ and AP $=0.84$ in cross-validation on PANORAMA, and AUC $=0.95$ and AP $=0.78$ on an external in-house cohort [2602.17986].

Patch-level fusion has been developed for transformer architectures. RE-ViT extracts 91 patch-level HRFs from six families and fuses them by early fusion with linearly projected image patches before ViT encoding. Reported results were AUC $=0.950\pm0.011$ on BUSI, $0.989\pm0.004$ on ChestXray2017, and $0.986\pm0.001$ on Retinal OCT, with statistically significant degradation when the radiomics stream was removed [2504.10916].

Not all effective HRF deployments require large generic feature banks. In hepatocellular carcinoma classification on 3D CT, three handcrafted radiological features—arterial-phase hyper-enhancement, non-peripheral wash-out, and enhancing capsule—were derived from lesion and parenchyma measurements motivated by LI-RADS criteria and fused with deep outputs in an $L_2$-regularized logistic regression. On the transfer set $D_2$, HF + size alone achieved $83.0\pm0.2\%$ AUC and the combined DLF + HF + size model achieved $83.2\pm0.3\%$, outperforming the baseline “Tiny” net at $77.8\pm1.8\%$ [2501.08097].

Prognostic and radiogenomic settings show a more mixed interaction between handcrafted and learned representations. In lung-cancer survival prediction, PET_HRF with semi-supervised learning reduced MAE from $2.11\pm0.46$ years to $1.55\pm0.25$ years, and CT_HRF with PCA plus Component-Wise Gradient Boosting Survival Analysis achieved external c-index $0.656\pm0.02$ [2502.01661]. In multicenter NSCLC radiogenomics for EGFR/KRAS prediction, IBSI-compliant HRF-only pipelines generalized best, reaching external AUC $0.77\pm0.07$ and accuracy $0.77\pm0.00$, whereas fusion improved robustness over DFR-only models but did not consistently outperform HRFs [2603.24922]. In multicenter NSCLC survival modeling, ComBat-harmonized tumor radiomics combined with clinical covariates achieved C-index $0.7552$ and 5-year t-AUC $0.8820$ [2505.17893].

HRFs have also been extended to retrieval and generation. RadiomicsRetrieval uses 72 tumor-level HRFs aligned with promptable image embeddings and anatomical positional embedding; radiomics-based retrieval correlations reached $0.9663$ on BraTS and $0.9448$ on NSCLC for Top-5 HRF queries [2507.08546]. Radiomics-conditioned tumor synthesis uses 16 shape features and 74 texture features as conditioning vectors in a GAN-based shape generator and a diffusion-based texture generator, while RadiomicsFill conditions glioma synthesis on 67 HRFs to control tumor geometry and heterogeneity [2509.24182] [2311.02586].

## 6. Interpretability, misconceptions, and open problems

The principal interpretive appeal of HRFs is that downstream decisions can be traced to named image phenotypes. RadiomicNet describes its attention modulation as a direct, differentiable function of named radiomics features with clear clinical meaning, and its learned scalar gate and channel weights are tied to specific texture descriptors [2607.02185]. Patient-specific radiomic feature selection for knee MRI retains fully interpretable radiomic features while learning a sparse per-patient subset, with case studies linking top-weighted features such as entropy, interquartile range, compactness, energy, uniformity, and the 10th percentile to clinically meaningful anatomy and pathology [2503.13131]. In multicenter CT radiogenomics, SHAP highlighted morphology- and heterogeneity-related radiomic phenotypes such as Morph-VolDens-AABB, coefficient of variation, moment invariants, GLDZM-SDHGE, and GLRLM-LRE as key predictors [2603.24922].

Interpretability, however, should not be conflated with automatic validity. Reviews emphasize persistent challenges of feature instability, limited reproducibility, validation bias, domain shift, and restricted clinical translation. The survey literature calls for harmonized protocols, unbiased validation frameworks, standardized multi-center benchmarks, federated learning, multimodal fusion, and regulatory-compliant, end-to-end reproducible workflows with uncertainty quantification [2602.00102]. The PET/SPECT review similarly reports that only $48\%$ of studies explicitly adhered to IBSI definitions, that class imbalance was addressed in only $59\%$, that missing data were discussed in $29\%$, and that low population diversity remained common; across 226 studies, deep radiomics achieved the highest mean accuracy ($0.862$) while fusion yielded the highest mean AUC ($0.861$), indicating that handcrafted radiomics is neither obsolete nor universally dominant [2507.16065].

A second misconception is that fusion with deep learning invariably improves performance. The empirical record is more conditional. Hybrid methods can be highly effective when HRFs inject heterogeneity priors, spatially resolved parametric maps, or clinically grounded descriptors into the network, but multicenter radiogenomics results show that HRF-only models may still generalize better than DFR-only or fused systems under strong site variation [2602.17986] [2603.24922]. This suggests that the enduring value of HRFs lies less in competing with learned features on representational capacity than in providing stable, inspectable, and domain-grounded structure for modeling under limited data, heterogeneous acquisition, and translational constraints.

Source: https://www.emergentmind.com/topics/handcrafted-radiomics-features-hrf