First-order statistics are histogram-based metrics that quantify the frequency distribution of voxel intensity values within a delineated region of interest (ROI). These descriptors—including mean, median, standard deviation, skewness, kurtosis, and entropy—treat the image as a collection of independent observations, ignoring the spatial arrangement of pixels. They provide a foundational characterization of overall tissue density or signal intensity on modalities such as CT, MRI, and PET.
Glossary
First-Order Statistics

What is First-Order Statistics?
First-order statistics describe the distribution of individual voxel intensity values within a region of interest without considering spatial relationships or pixel interactions.
Because they disregard spatial context, first-order features are computationally efficient but cannot capture textural heterogeneity or morphological patterns. They are sensitive to acquisition parameters and require intensity discretization and Z-score normalization to ensure reproducibility across scanners. Despite their simplicity, these metrics serve as critical inputs to radiomic signatures and are standardized under the Image Biomarker Standardisation Initiative (IBSI) guidelines.
Core First-Order Metrics
First-order statistics describe the distribution of individual voxel intensity values within a region of interest (ROI) without considering spatial relationships or pixel interactions. These metrics form the foundational layer of radiomic analysis, quantifying tumor heterogeneity through simple histogram properties.
Mean Intensity
The arithmetic average of all voxel intensity values within the segmented ROI. This metric provides a baseline measure of overall tissue density or signal activity.
- Clinical relevance: Distinguishes hypo-intense from hyper-intense lesions
- Calculation: Sum of all voxel intensities divided by total voxel count
- Limitation: Highly sensitive to outliers and does not capture heterogeneity
- Example: A mean Hounsfield Unit of 35 in a lung nodule suggests ground-glass opacity rather than solid tissue
Variance and Standard Deviation
Variance measures the spread of intensity values around the mean, while standard deviation is its square root. These metrics quantify the degree of dispersion within the ROI.
- High variance indicates heterogeneous tissue with mixed densities
- Low variance suggests homogeneous tissue composition
- Clinical application: Tumor heterogeneity often correlates with malignancy and treatment resistance
- Standard deviation is expressed in the original intensity units, making it more interpretable than variance
Skewness
A measure of the asymmetry of the intensity histogram distribution. Skewness indicates whether the tail of the distribution extends toward higher or lower intensity values.
- Positive skew: Tail extends toward higher intensities; mean > median
- Negative skew: Tail extends toward lower intensities; mean < median
- Zero skew: Symmetric distribution around the mean
- Clinical insight: Skewness shifts can indicate the presence of enhancing or necrotic sub-regions within tumors
Kurtosis
A measure of the 'peakedness' or tailedness of the intensity histogram relative to a normal distribution. Kurtosis quantifies the concentration of values around the mean versus the extremes.
- High kurtosis (leptokurtic): Sharp peak with heavy tails; indicates a dense core with outlier voxels
- Low kurtosis (platykurtic): Flat distribution; suggests uniform heterogeneity
- Excess kurtosis is often reported, where a normal distribution equals zero
- Example: A necrotic tumor core surrounded by an enhancing rim produces high kurtosis values
Entropy
A measure of the inherent randomness or unpredictability in the intensity distribution. Higher entropy indicates greater heterogeneity and disorganization within the ROI.
- Calculation: Based on Shannon's information theory applied to the intensity histogram
- Maximum entropy occurs when all intensity bins are equally probable
- Minimum entropy occurs when all voxels share the same intensity value
- Clinical significance: High entropy in glioblastomas correlates with poorer overall survival and genetic heterogeneity
Percentiles and Range
Order statistics that describe the intensity values at specific cumulative distribution thresholds, along with the overall intensity range.
- 10th and 90th percentiles: Robust alternatives to minimum and maximum, less sensitive to outlier voxels
- Interquartile range (IQR): Difference between 25th and 75th percentiles; robust measure of dispersion
- Intensity range: Maximum minus minimum; highly sensitive to noise and artifacts
- Median: The 50th percentile; robust central tendency measure unaffected by extreme values
How First-Order Statistics Are Computed
First-order statistics are computed directly from the frequency distribution of voxel intensity values within a segmented region of interest, without any consideration for spatial relationships or pixel interactions.
The computation begins with intensity discretization, where continuous Hounsfield Unit or signal intensity values are binned into a finite number of discrete gray levels—typically between 8 and 256 bins—to reduce noise sensitivity. A histogram is then constructed by tallying the frequency of each discretized intensity value across all voxels in the region of interest (ROI). From this single-variable distribution, scalar metrics are derived: the mean calculates the average intensity, variance measures the spread around the mean, skewness quantifies distribution asymmetry, and kurtosis evaluates the peakedness or tailedness relative to a normal distribution.
Additional descriptors include energy, which sums the squared values of histogram bin probabilities to measure uniformity, and entropy, which quantifies randomness using Shannon's information theory formula. Percentile-based metrics such as the 10th, 25th, 50th (median), 75th, and 90th percentiles capture distribution extremes without parametric assumptions. The uniformity metric directly sums the squared frequencies of each intensity bin, while the interquartile range provides a robust measure of dispersion. All calculations are standardized by the Image Biomarker Standardisation Initiative (IBSI) to ensure reproducibility across different software implementations and imaging platforms.
Frequently Asked Questions
Clear, technical answers to common questions about histogram-based radiomic features and their role in quantifying tumor heterogeneity.
First-order statistics are quantitative metrics that describe the distribution of individual voxel intensity values within a segmented region of interest (ROI) without considering spatial relationships or pixel interactions. These histogram-based features include the mean, median, standard deviation, skewness, kurtosis, energy, and entropy of the gray-level frequency distribution. Unlike second-order texture matrices such as GLCM or GLRLM, first-order statistics treat the ROI as a bag of voxels, analyzing only the statistical properties of the intensity histogram. They serve as the foundational layer of radiomic feature extraction, providing baseline descriptors of tumor density, heterogeneity, and overall signal characteristics that often correlate with cellularity, necrosis, or treatment response.
Enabling Efficiency, Speed & Accuracy
Intelligent Analysis, Decision & Execution
We build AI systems for teams that need search across company data, workflow automation across tools, or AI features inside products and internal software.
Talk to Us
Search across company data
Give teams answers from docs, tickets, runbooks, and product data with sources and permissions.
Useful when people spend too long searching or get different answers from different systems.

Automate internal workflows
Use AI to route work, draft outputs, trigger actions, and keep approvals and logs in place.
Useful when repetitive work moves across multiple tools and teams.

Add AI to products and internal tools
Build assistants, guided actions, or decision support into the software your team or customers already use.
Useful when AI needs to be part of the product, not a separate tool.
First-Order vs. Higher-Order Statistics
Comparative analysis of statistical feature families based on spatial complexity and information captured from medical imaging regions of interest.
| Feature | First-Order Statistics | Second-Order Statistics | Higher-Order Statistics |
|---|---|---|---|
Spatial Relationship Consideration | |||
Computational Complexity | Low | Moderate | High |
Primary Input Data | Histogram of voxel intensities | Pairs of voxels at defined offsets | Runs, zones, or neighborhoods of voxels |
Key Texture Information Captured | Global intensity distribution | Local spatial arrangement and contrast | Regional homogeneity and structural patterns |
Typical Feature Count | 14-18 features | 22-26 features (per matrix) | 16-20 features (per matrix) |
Sensitivity to ROI Rotation | |||
Example Metrics | Mean, Variance, Skewness, Kurtosis, Entropy | Contrast, Correlation, Energy, Homogeneity (GLCM) | Short Run Emphasis (GLRLM), Zone Percentage (GLSZM) |
Primary Clinical Utility | Baseline tumor intensity characterization | Tissue texture and heterogeneity quantification | Advanced structural pattern and roughness analysis |
Related Terms
Explore the core statistical and methodological concepts that complement first-order histogram analysis in the radiomics pipeline.
Shape Features
Morphological descriptors that quantify the three-dimensional geometric properties of a segmented region of interest. Unlike first-order statistics, these metrics are independent of intensity values and focus purely on form.
- Volume: Total number of voxels multiplied by voxel dimensions
- Sphericity: How closely the ROI approximates a perfect sphere
- Compactness: Ratio of volume to surface area
- Maximum 3D Diameter: The largest pairwise Euclidean distance between surface mesh vertices
Gray-Level Co-occurrence Matrix (GLCM)
A second-order statistical method that quantifies texture by calculating the frequency of specific pairs of pixel intensities occurring at a defined spatial offset. While first-order statistics ignore spatial relationships, GLCM captures directional patterns.
- Generates matrices for multiple angles (0°, 45°, 90°, 135°)
- Derives features like contrast, homogeneity, and correlation
- Sensitive to the chosen distance offset and discretization bin width
Entropy
A first-order statistical measure of the randomness or inherent unpredictability in the distribution of voxel intensity values. High entropy indicates a heterogeneous, noisy region with many intensity levels; low entropy suggests a homogeneous, uniform texture.
- Calculated as:
-Σ p(i) * log₂(p(i))where p(i) is the probability of intensity i - Often correlates with tumor heterogeneity and treatment resistance
- Highly dependent on the number of discretization bins chosen
Intensity Discretization
The process of binning continuous voxel intensity values into a finite number of discrete gray levels. This is a critical preprocessing step that directly impacts the reproducibility of first-order and texture features.
- Fixed Bin Width: Bins of constant Hounsfield Unit width (e.g., 25 HU)
- Fixed Bin Count: A constant number of bins across all ROIs
- IBSI guidelines recommend reporting discretization parameters for reproducibility
Image Biomarker Standardisation Initiative (IBSI)
An independent international collaboration providing consensus-based reference values and standardized nomenclature for radiomic feature computation. IBSI ensures that first-order statistics calculated in different software packages yield identical results.
- Publishes benchmark datasets with ground-truth feature values
- Defines mathematical formulas for 174 radiomic features
- Critical for multi-center clinical trial harmonization
Z-Score Normalization
A feature scaling technique that standardizes radiomic feature values by centering them to a mean of zero and scaling to a standard deviation of one. Essential before feeding first-order statistics into machine learning classifiers.
- Formula:
z = (x - μ) / σ - Preserves the shape of the distribution while removing scale effects
- Prevents features with larger numeric ranges from dominating model training

About the author
Prasad Kumkar
CEO & MD, Inference Systems
Prasad Kumkar is the CEO & MD of Inference Systems and writes about AI systems architecture, LLM infrastructure, model serving, evaluation, and production deployment. Over 5+ years, he has worked across computer vision models, L5 autonomous vehicle systems, and LLM research, with a focus on taking complex AI ideas into real-world engineering systems.
His work and writing cover AI systems, large language models, AI agents, multimodal systems, autonomous systems, inference optimization, RAG, evaluation, and production AI engineering.
Partnered with leading AI, data, and software stack.
How We Work
Custom AI workflows for your Business
One-fit-all AI don't work for modern businesses. At Inferensys, we aim to understand your business & custom requirements; which we use to define most efficient agentic workflows, the data, and the tools for your business.
01
Review the use case
We understand the task, the users, and where AI can actually help.
Read more02
Pick the right approach
We define what needs search, automation, or product integration.
Read more03
Build the first useful version
We implement the part that proves the value first.
Read more04
Improve from there
We add the checks and visibility needed to keep it useful.
Read moreThe first call is a practical review of your use case and the right next step.
Talk to Us