Modeling flow and transport in porous media using pore‐scale modeling is reliant on rock properties derived from digital rock images using segmentation techniques. These digital rock images obtained using computed tomography incorporate the variation in the intensity of phases depending on the attenuation of X‐rays. A standard technique is the segmentation of tomographic images based on user‐selected grayscale thresholding, allowing the identification of different phases. This threshold is subjective based on the operator and results in loss of essential information about the grayscale variation after segmentation. This paper implements the gray‐level co‐occurrence matrix (GLCM) incorporating the full range of grayscale information. The GLCM captures the relative occurrence of grayscale values in a spatial map. These maps show visually connected/disconnected populations of different phases such as pore space, quartz grains, minerals, and other features. We show that each rock has its own GLCM signature depending on the variations in gray‐level intensities. Several statistical measures are calculated: (1) GLCM contrast describing local variation in the gray‐level intensities, (2) GLCM angular second moment, describing the rock homogeneity; (3) GLCM mean, describing weighted average of the probability of occurrence of features based on their location on the GLCM map; and (4) GLCM correlation, measuring the linear dependencies of grayscale values and the degree of (an) isotropy in the micro–computed tomographic images of each of the rock types. The GLCM method provides a pathway to alleviate user biases and allow automation of micro–computed tomography analyses.