Hereby, we proposed a protocol to illustrate the effect of aggregate surface morphology on the ITZ microstructure. The SEM-BSE image were quantitatively analyzed to obtain ITZ's porosity gradient via digital image processing and a K-means clustering algorithm was further employed to establish a relationship between porosity gradient and surface roughness.