Here, we present a comprehensive method to illustrate the uneven distribution of the interfacial transition zone (ITZ) around the aggregate and the effect of aggregate surface morphology on the formation of ITZ. First, a model concrete sample is prepared with a spherical ceramic particle in roughly the central part of the cement matrix, acting as a coarse aggregate used in common concrete/mortar. After curing until the designed age, the sample is scanned by X-ray computed tomography to determine the relative location of the ceramic particle inside the cement matrix. Three locations of the ITZ are chosen: above the aggregate, on the side of the aggregate, and below the aggregate. After a series of treatments, the samples are scanned with a SEM-BSE detector. The resultant images were further processed using a digital image processing method (DIP) to obtain quantitative characteristics of the ITZ. The surface morphology is characterized at the pixel level based on the digital image. Thereafter, K-means clustering method is used to illustrate the effect of surface roughness on ITZ formation.